微分方程为：在附件中。
我编的程序如下：
x(1)=x;x(2)=y;z=[x(1),x(2)]';
function xdot=shier(t,z)
r=1;d=0.5;a=0.1;b=0.02;
xdot=[r-a*x(2),0;0,-d+b*x(1)]*z;
ts=0:0.1:15;
x0=[25,2]';
[t,z]=ode45(@shier,ts,x0);
[t,z]
plot(t,z),grid,
gtext('\fontsize{12}x(t)'),gtext('\fontsize{12}y(t)'),
pause,plot(z(:,1),z(:,2)),grid,
xlabel('x'),ylabel('y')

运行后出现：
??? Error: File: shier.m Line: 2 Column: 1
Function definitions are not permitted at the prompt or in scripts.
这是为什么？

[此贴子已经被作者于2007-7-20 23:46:41编辑过]