Chen系统

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

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