pangjialai.m

来自「混沌映射庞加莱界面Matlab程序」· M 代码 · 共 13 行

clear
ff=0.222;
options=odeset('RelTol',1e-7);
tt=2*pi/1.2
[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff);
plot(x(2000:end,1),x(2000:end,2),'-')
%pojialai
hold on
i=2000:100:3000
plot(x(i,1),x(i,2),'*')