用Runge-Kutta法解非线性微分方程

songchangsha

我的问题是用Runge-Kutta法解非线性微分方程，但没法看数值结果。
如微分方程为y'=y－x/y
M文件为function y=ff1(x,y)
       y=y－x/y
命令为：
options=odeset('RelTol',1e-4,'AbsTol',[1e-4]);
ode45=('ff1',[0,1],1,options
现在可以看解对应的图，就是没法看x 具体结果,也没法看x对应的y值。
谢谢！！！

catny

ode45可以输出结果的
[t,y]=ode45(...)

yydeep

应该做两个函数
如：
function xdot=sys.m
xdot=(y,y-xy)
function simu.m
[t,y]=ode45('sys',t0,tfinal,[x,y],tol)