mathematica求解微分方程的数值解,做出图像,R = 0.008;Y = 47*2*3.14;z = 0.18;m = 0.14;k = 20;g = 10;f = 20;i = 0.143;e = 28;NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) + 0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0.5*k*z^2 - k*z*Sqrt[R*Y^2 + z^2 - R*y[t]^2] + 0.5*f*(y[t] - Y)^2 - m*g*(Sqrt[R*Y^2 + z^2 - R*y[t]^2] - z) == e,y[0] == 0,y'[0] == 1},y[t],{t,10}];Plot[Evaluate[y[t] /.%],{t,0,10}]这是代码,我用的mathematica 8 看不出错误来.

问题描述:

mathematica求解微分方程的数值解,做出图像,
R = 0.008;
Y = 47*2*3.14;
z = 0.18;
m = 0.14;
k = 20;
g = 10;
f = 20;
i = 0.143;
e = 28;
NDSolve[{0.5*m*(R^4 *y[t]^2 *y'[t]^2)/(R*Y^2 + z^2 - R*y[t]^2) +
0.5*i*y'[t]^2 + 0.5*k*(R*Y^2 + z^2 - R*y[t]^2) + 0.5*k*z^2 -
k*z*Sqrt[R*Y^2 + z^2 - R*y[t]^2] + 0.5*f*(y[t] - Y)^2 -
m*g*(Sqrt[R*Y^2 + z^2 - R*y[t]^2] - z) == e,y[0] == 0,
y'[0] == 1},y[t],{t,10}];
Plot[Evaluate[y[t] /.%],{t,0,10}]
这是代码,我用的mathematica 8 看不出错误来.