求助：用龙格库塔法求微分方程

[liuwei7854]

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

[anbcjys]

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

[sunjq]

请试试
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]);