|
发表于 2016-4-2 11:09:31
|
显示全部楼层
简单的做法就是用Matlab解好方程,然后把解用C语言写出来就行了
>> syms x y
>> solve(y==-0.8*x^3+20000*x,x)
ans =
25000/(3*(((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)) + (((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)
- 12500/(3*(((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)) - (((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)/2 + (3^(1/2)*(25000/(3*(((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)) - (((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3))*i)/2
- 12500/(3*(((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)) - (((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)/2 - (3^(1/2)*(25000/(3*(((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3)) - (((25*y^2)/64 - 15625000000000/27)^(1/2) - (5*y)/8)^(1/3))*i)/2
|
|