超混沌分岔图

clc
clear all
XA=100;XC=67;XD=100;XE=70;XF=90;XR=69;XG=80;
Z=[];
for XB=linspace(42,70,10);
    options = odeset('RelTol',1e-6,'AbsTol',[0.1 0.1 0.1 0.1]);
    [T,X]=ode45('chen',[0,50],[0.1 0.1 0.1 0.1],options,XA,XB,XC,XD,XE,XF,XR,XG);
  
    n=length(X);
    for k=round(n/2):n
      if abs(X(k,1))<1
        Z=[Z,XB+abs(X(k,2))*i];
      end
    end
end
figure
plot(Z,'.','markersize',1);
xlabel('b'),ylabel('|x|');