楼主你好!
这个方程是解不出解析解的。
我用ODE45给您算了一下。
先赋个初值y(0)=10000
编写m文件
function f=exam2fun (t,y)
f=3.9464/(y*(0.8-y*0.2679))*(3.8873e3*y+7.6004e-006+0.0299*sqrt(0.8333*y));
f=f(:);
再在主窗口输入
[t,y]=ode45('exam2fun',[0:1:1000],10000)
plot(t,y)
可以求得0到1000内的y(t)
最后有一部分为虚数,若增大初值,虚数分界线变大
得到的图像如图,最后一部分虚数没有显示。
您可以选择自己的初值,然后选择范围求解并画图。
温馨提示:答案为网友推荐,仅供参考