请教一个ODE
帮别人问个问题,请各位指点:下面的程序哪里除了问题?
function dy=myfun_3(t,y)
globalm 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 = ;
y0 = ;
=ode45(@myfun_3, Tspan, y0);% 是命令行
plot(tout,yout(:,1),'-')
hold on
plot(tout,yout(:,2),'--')
legend('x','y')
xlabel('time')
ylabel('density') 没有看见有啥错误，下面的是你想要得吗? 只是改了一下顺序。注意下面的m文件的文件名是"myabc.m".
function myabc
Tspan = ;
y0 = ;
=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));
页:
[1]