thenon5.m

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

epsb=1e-15;Arc=9;dk=1e-3;amax=0.5;amin=0.1; dmin=0.1;damax=0.3;damin=0.1;dmax=0.2; %-3

p0=[-1.62114638486435;-1.62114638486435];p=[-0.4;1.03;-0.1;0];%alpha,beta,R,S

lamb=3.12684825232441; 
del=-1e-3;
v=[-0.30461229725090;-0.95247642929551];
n=1;
 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),'-r.'), end
 hold on
 plot(-1.62114638486435,-1.62114638486435,'--sg'),hold on,