关于2-D微分方程的MATLAB求解，请大家赐教

mingpangzi

我编了2个m文件，求解二维的偏微分方程dx/dt=x-y;dy/dt=1-x*y2; m文件如下：
m1文件f.m：
function dxdt=f(t,x)
dxdt1=x(1)-x(2);
dxdt2=1-x(1)*x(2)^2;
dxdt=[dxdt1;dxdt2];
end end
m2文件simu.m：
function [t,x]= simu(t,x);
t0=0; te=60;
x0=[1;3];
[t,x] = ode45(@f, [t0 te],x0,[]); plot(t,x(1));
plot(t,x(2));
plot(x(1),x(2));
不管我怎么调初始值x0，得到的x1和x2对时间变化的图形总是一条水平直线，即不随时间变化。
为什么？请达人赐教！

[ 本帖最后由 mingpangzi 于 2009-9-30 04:43 编辑 ]

mingpangzi

诚意请教，请会的人指点一二！