请教一个ODE

lzh1975

帮别人问个问题,请各位指点:
下面的程序哪里除了问题?
function dy=myfun_3(t,y)
global  m k a b
m=0.5;k=1;a=1;b=1;
dy=zeros(2,1);
dy(1,1)=y(1).*(1-y(1))-y(1).*y(2).^m./(1+k*k*y(1).*y(1));
dy(2,1)=a*y(2).*(1-b*y(2)./y(1));
Tspan = [0 100];
y0 = [0.1; 0.2];
[tout, yout]=ode45(@myfun_3, Tspan, y0);% 是命令行
plot(tout,yout(:,1),'-')
hold on
plot(tout,yout(:,2),'--')
legend('x','y')
xlabel('time')
ylabel('density')

xyz999

没有看见有啥错误，下面的是你想要得吗? 只是改了一下顺序。注意下面的m文件的文件名是"myabc.m".

function myabc
Tspan = [0 100];
y0 = [0.1; 0.2];
[tout, yout]=ode45(@myfun_3, Tspan, y0);
plot(tout,yout(:,1),'-')
hold on
plot(tout,yout(:,2),'--')
legend('x','y')
xlabel('time')
ylabel('density')

function dy=myfun_3(t,y)

m=0.5;k=1;a=1;b=1;
dy=zeros(2,1);
dy(1,1)=y(1).*(1-y(1))-y(1).*y(2).^m./(1+k*k*y(1).*y(1));
dy(2,1)=a*y(2).*(1-b*y(2)./y(1));