%º¯Êý
function dxdt=differential(t,x)
dxdt=[x(1)*(1-x(1))*(5-6*x(2));x(2)*(1-x(2))*(1-6*x(1))]; end
%Ö÷³ÌÐò
clear
%y-x
for i=0:0.1:1
for j=0:0.1:1
[T,Y]=ode45('differential',[0 5],[i j]); figure(1)
grid on
plot(Y(:,1),Y(:,2));
hold on
end
end
%dx/dt-t
for i=0:0.1:1
for j=0:0.1:1
[T,Y]=ode45('differential',[0 5],[i j]);
figure(2)
grid on
plot(T,Y(:,1));
hold on
end
end
%dy/dt-t
for i=0:0.1:1
for j=0:0.1:1
[T,Y]=ode45('differential',[0 5],[i j]); figure(3)
grid on
plot(T,Y(:,2));
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余1页未读,立即下载