作业帮 > 数学 > 作业

matlab数值解法求解二阶微分方程 ODE45函数

来源:学生作业帮 编辑:拍题作业网作业帮 分类:数学作业 时间:2024/05/01 23:51:36
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呢?
因为你x=0时2/x是无穷大呀,然后y'又是0,然后(2/x)y'就是nan了,所以后面算的全是nan了.