Rossler系统

function dx=Rossler(t,x);
% Rossler方程[不显含时间t的自治系统]
% 方程如下:
% dx/dt=-y-z
% dy/dt=x+a*y
% dz/dt=b+z*(x-c)
%
% Example:
% [T,Y]=ode45('Rossler',[0,400],[1;1;1]);
% plot3(Y(:,1),Y(:,2),Y(:,3));
%
% Author's email: ustb03-07@yahoo.com.cn
%
a=.15;
b=.2;
c=10;
dx=zeros(3,1);

dx(1)=-x(2)-x(3);
dx(2)=x(1)+a*x(2);
dx(3)=b+x(3)*(x(1)-c);