function z=stk(t,y)
z(1,:)=y(1)-0.015*y(1).*y(2);
z(2,:)=-y(2)+0.01*y(1).*y(2);
function text
Y0=[100;20];
[t,Y]=ode23(@stk,[0,20],Y0);
x=Y(:,1);y=Y(:,2);
plot(x,y)
错误是
>> stk
错误使用 stk (line 2)
输入参数的数目不足。
由于自定义的微分方程函数 z=stk(t,y) 定义不完善而出现错误。
完整的函数应为
function z=stk(t,y)
z(1,:)=y(1)-0.015*y(1).*y(2);
z(2,:)=-y(2)+0.01*y(1).*y(2);
z=z(:);
end
更改后运行,可得到其图形。