y'+y-y^(2/3)=0; x=0:0.1:1 怎么用matlab求解?

问题描述:

y'+y-y^(2/3)=0; x=0:0.1:1 怎么用matlab求解?

这是个微分方程,需要有初始条件才能求解.

假设初始条件y(0)=0.5,则代码如下:

dy=@(x,y)-y+y^(2/3);
ode45(dy,0:0.1:1,0.5)

如果我想输出y值的矩阵呢?

把第二句改成

[t,y]=ode45(dy,0:0.1:1,0.5);