[laosam280]

代码我写好了：
%这就是一个例子，初值我都自己定义了，你要计算可以修改。把上面的X1'=AX0+(B+1)X1+X1^2*X2；X2'=X0+BX1-X1^2*X2这个方程组写成一个函数文件，存为myfunfun.m

t_final=100; x0=[1;1];
[t,x]=ode45('myfunfun',[0,t_final],x0);
plot(t,x)

可惜了，因为有非线性项，我电脑没有在短时间内算出来。明天我去实验室用大型机给你算一下，再找个时间把结果发给你。