autocatalator.m

计算3维系统分叉点

% Author: Housam Binous

% Three variable autocatalator

% National Institute of Applied Sciences and Technology, Tunis, TUNISIA

% Email: [e-mail]

function xdot=autocatalator(t,x)

global nu

kappa=10;
delta=2e-2;
sigma=5e-3;

xdot(1)=kappa+nu*x(3)-x(1)*x(2)^2-x(1);
xdot(2)=(x(1)*x(2)^2+x(1)-x(2))/sigma;
xdot(3)=(x(2)-x(3))/delta;

xdot=xdot';

end