matlab数值解法求解二阶微分方程 ODE45函数y‘’+(2/x)y'+y^(3/2)=0初值x=0时 y=1 y'=0 解的区间随意 就取【0,20】吧说下具体怎么用这个函数求解...或者帮忙看下错误吧：m文件function dydt = eq1(t,y)dydt(1)=

Question

matlab数值解法求解二阶微分方程 ODE45函数
y‘’+(2/x)y'+y^(3/2)=0
初值x=0时 y=1 y'=0
解的区间随意 就取【0,20】吧
说下具体怎么用这个函数求解...
或者帮忙看下错误吧：
m文件
function dydt = eq1(t,y)
dydt(1)=y(2);
dydt(2)=-y(1)^(3/2)-(2/t)*y(2);
dydt=dydt';
end
命令窗口：
[t,y]=ode45('eq1',[0 15],[1 0]);
为什么y的值全部是NaN呢?

Answer 1

因为你x=0时2/x是无穷大呀,然后y'又是0,然后(2/x)y'就是nan了,所以后面算的全是nan了.