Question

微分方程求解,跪求！！

y^3dx+2(x^2-xy^2)dy=0，当x=1时y=1

Answer 1

>> f=inline('1/6-y/30','t','y');
>> [t,y]=ode45(f,[0,5],[0]);
>> plot(t,y)

另外，由高数的知识，可以得到解析解为：y=5-5*e^(-t/30)
把上述的解析值也放到上图中。

>> hold on
>> plot(t,5-5*exp(-t/30),'r*')
可以看到二者相当吻合。

追问

话说，没看懂...可以用求微分方程的方法来说不啦/(ㄒoㄒ)/~~谢谢