%主程序
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) %图(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));
hold on
end
end
评论1