1原文程序
function dx=paper6 (t,x)
dx=[x(2);
x(3);
5.5*x(1)-3.5*x(2)-1*x(3)-1*x(1)*x(1)*x(1)];
%%%%%%%%%%%%%%%%%%%%%%%%%%
x0=[0.1; 0.3; 0.5];
[t,y]=ode45('paper6',[0,300],x0);
for i=1:3
figure(i)
plot(t,y(:,i))
end
figure(4)
plot(y(:,1),y(:,2))
figure(5)
plot(y(:,1),y(:,3))
figure(6)
plot(y(:,2),y(:,3))
figure(7)
plot3(t,y(:,1),y(:,2))%,grid,
figure(8)
plot3(t,y(:,1),y(:,3))%,grid,