thenon1.m

来自「计算动力学系统的分岔图」· M 代码 · 共 24 行

epsb=1e-7;Arc=7.5;dk=0.01;amax=0.69;amin=0.2;dmin=0.001;damax=0.07;damin=0.001;dmax=0.2;  %1

p0=[0.46661702380495;0.46661702380495];p=[0.3;-1.057;-0.5;0];%alpha,beta,R,S
lamb=2.76261564458262; %FP2
del=1e-1;
v=[0.51553133957551;-0.85687072415592];%Fp2
n=2;
MM=Umanifold(p0,Arc,amax,amin,dmax,dmin,damax,damin, dk,p,epsb,lamb,del,v,n);
global M
for i=1:size(M,2) hold on, plot(M(1,i),M(2,i),'-b.'), end
 hold on
 plot(0.46661702380495,0.46661702380495,'--sg'),hold on,

del=-1e-3;
MM=Umanifold(p0,Arc,amax,amin,dmax,dmin,damax,damin, dk,p,epsb,lamb,del,v,n);
global M
for i=1:size(M,2) hold on, plot(M(1,i),M(2,i),'-b.'), end
 hold on
 plot(0.46661702380495,0.46661702380495,'--sg'),hold on,