thenon2.asv

计算动力学系统的分岔图

% 
epsb=1e-9;Arc=5;dk=1e-3;amax=0.5;amin=0.2; dmin=0.01;damax=0.7;damin=1e-3;dmax=1e-2;
p0=[-0.42861702380495;-0.42861702380495];
p=[0.3;-1.057;-0.5;0];
p1=p0-1e-2*[-0.81439328674458;0.58031334165721];
n=2;
MM=Smanifold(p0,p1,Arc,amax,amin,dmax,dmin,damax,damin, dk,p,epsb,n);
global SM
for i=1:size(SM,2) hold on, plot(SM(1,i),SM(2,i),'-r.'), end
hold on
plot(-0.42861702380495,-0.42861702380495,'--sg'),
pause
epsb=1e-15;Arc=5;dk=0.5e-3;amax=0.4;amin=0.1; dmin=0.01;damax=0.07;damin=1e-5;dmax=0.0001;
p1=p0+1e-2*[-0.81439328674458;0.58031334165721];%p0+1e-3
MM=Smanifold(p0,p1,Arc,amax,amin,dmax,dmin,damax,damin, dk,p,epsb,n);
global SM
for i=1:size(SM,2) hold on, plot(SM(1,i),SM(2,i),'-b.'), end
hold on