solution.m

来自「对于Rossler方程组这样的非线形微分方程组」· M 代码 · 共 17 行

global a; 
global b; 
global c; 
b=2; 
c=4; 
t0=[0,200]; 
for a=0:0.02:0.65 
    [t,x]=ode45('rossler',t0,[0,0,0]); 
    a 
    subplot(1,2,1);
    plot(t,x(:,1),'r',t,x(:,2),'g',t,x(:,3),'b');         
    title('x(红色),y( 绿色),z(篮色)随t变化情况');xlabel('t')
    subplot(1,2,2); 
        plot3(x(:,1),x(:,2),x(:,3))     
        title('相图');xlabel('x');ylabel('y');zlabel('z'); 
        pause 
end