main_autocatalator.m

计算3维系统分叉点

% Author: Housam Binous

% Three variable autocatalator

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

% Email: [e-mail]

clc

clear

global t x y z nu

figure(2)

for nu=0.10:0.0001:0.20

nu

[t x]=ode45(@autocatalator,0:0.001:14,[0.01 0 0]);

y=diff(log10(x(4000:end,3)))/0.001;

z=diff(y)/0.001;

k=1;

clear a

for i=1:100
    t0=4+0.08*i;
    option=optimset('display','off');
    zer=fsolve(@diff_autocatalator,t0,option);
    if interp1(t(4000:end-2),z,zer)>0 
        a(k)=interp1(t(4000:end),log10(x(4000:end,3)),zer); k=k+1;
    end
end
kmax=k-1;
h=plot(nu.*ones(1,kmax),a,'r.');
hold on
set(h,'MarkerSize',0.1);
end