thenon6.m

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

epsb=1e-12;Arc=8;dk=1e-3;amax=0.5;amin=0.2; 
dmin=0.01;damax=0.7;damin=1e-3;dmax=1e-2;
p0=[ -1.62114638486435; -1.62114638486435];p=[-0.4;1.03;-0.1;0];
p1=p0-1e-2*[ -0.96357221418883  ;-0.26744829041001  ];
n=1;
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