求助：用龙格库塔法求微分方程 - MATLAB爱好者论坛-LabFans.com

liuwei7854 · 2009-09-06, 22:56

用龙格库塔法（四阶）求问分方程组
dy/dx=1/(z-x);
dz/dx=-1/y+1;
y(0)=1,z(0)=1,x在[0,1]。
希望能得到帮助，向高手致敬！！

anbcjys · 2009-09-07, 08:37

按照龙格库塔法给出的微分方程表示形式来描述你的方程就行了

sunjq · 2009-09-23, 17:26

请试试
function dx=myf(t,x)
dx=zeros(3,1)
dx(1)=1;
dx(2)=1/(x(3)-x(1));
dx(3)=-1/x(2)+1;

[t,y]=ode45('myf',[0 1],[0 x0 y0]);

2009-09-06, 22:56	#1
liuwei7854 初级会员注册日期: 2009-05-02 年龄: 39 帖子: 2 声望力: 0	求助：用龙格库塔法求微分方程用龙格库塔法（四阶）求问分方程组 dy/dx=1/(z-x); dz/dx=-1/y+1; y(0)=1,z(0)=1,x在[0,1]。希望能得到帮助，向高手致敬！！

2009-09-07, 08:37	#2
anbcjys 高级会员注册日期: 2008-09-14 年龄: 43 帖子: 351 声望力: 24	回复: 急，求助：用龙格库塔法求微分方程按照龙格库塔法给出的微分方程表示形式来描述你的方程就行了 __________________ qq604443022

2009-09-23, 17:26	#3
sunjq 初级会员注册日期: 2009-09-23 帖子: 2 声望力: 0	回复: 急，求助：用龙格库塔法求微分方程请试试 function dx=myf(t,x) dx=zeros(3,1) dx(1)=1; dx(2)=1/(x(3)-x(1)); dx(3)=-1/x(2)+1; [t,y]=ode45('myf',[0 1],[0 x0 y0]);