rosser.m

rossler分岔代码

function u=rosser(t,x)
global c
a=0.2;b=0.2;
x0=[0.5,0.4,0.3]';
u=[-x(2)-x(3),x(1)+a*x(2),b+x(3)*(x(1)-c)]';
% function dy=rosser(t,y)
% global c
% dy=zeros(3,1);
% dy(1)=-(y(2)+y(3));
% dy(2)=y(1)+0.2*y(2);
% dy(3)=0.2+y(3)*(y(1)-c);