[MATLAB数学相关] 请教二阶非线性微分方程求解 - MATLAB爱好者论坛-LabFans.com

lxcwym · 2009-08-10, 15:29

y-1=d^2(m^2/2/y^2)/dx^2+d^2(lny)/dx^2
y(inf)=1
y'(inf)=0

lxcwym · 2009-08-10, 16:11

function dy=plas(x,y)
dy=zeros(5,1);
y(1)=1/(2*y(3)^2);
y(2)=log(y(3));
dy(1)=y(4);
dy(2)=y(5);
dy(4)=y(3)-1-dy(5);

>>[x,y]=ode45(@plas,[-inf,inf],[1,0])

我这样编写，问题出现在哪里呢？谢谢！

2009-08-10, 15:29	#1
lxcwym 初级会员注册日期: 2009-08-08 年龄: 45 帖子: 8 声望力: 0	请教二阶非线性微分方程求解 y-1=d^2(m^2/2/y^2)/dx^2+d^2(lny)/dx^2 y(inf)=1 y'(inf)=0

2009-08-10, 16:11	#2
lxcwym 初级会员注册日期: 2009-08-08 年龄: 45 帖子: 8 声望力: 0	回复: 请教二阶非线性微分方程求解 function dy=plas(x,y) dy=zeros(5,1); y(1)=1/(2*y(3)^2); y(2)=log(y(3)); dy(1)=y(4); dy(2)=y(5); dy(4)=y(3)-1-dy(5); >>[x,y]=ode45(@plas,[-inf,inf],[1,0]) 我这样编写，问题出现在哪里呢？谢谢！