为什么我的分岔图出不来

ymsncx

function dx=yang(t,x)
global A ;
dx=zeros(3,1);
dx(1)=6.3*x(2)-1.8*x(1)-A*(abs(x(1)-1)-abs(x(1)+1));
dx(2)=0.7*x(1)-0.7*x(2)+x(3);
dx(3)=-7.0*x(2);

function chua
clear;
global A
range=[0.1: 0.0001: 0.35];
period=2*pi;
k=0;
YY1=[];
Step=2*pi/100;
for A=range
  x0=[0 0.001 0];
  A
  k=k+1;     %出去前面60个周期的数据，并将最后的结果作为下一次积分的初值。
  tspan=[0:step:60*period];
  [t,Y]=ode45('yang',tspan ,x0);
  y0=Y(end,:);
  j=1;
  for i=60:200
    tspan=[i*period:step(i+1)*period];
    [t,Y]ode45('yang',tspan ,x0);
    YY1(k, j)=Y(end,1);
    j=j+1;    %取出每一个周期内的第一个解的最后一个值。
    y0=Y(end,:);
   end
end
bifdata=YY1(:,end-51:end);
plot(range,bifdata,'k','markersize',1)；
各位请都给点意见吧

咕噜噜

[t,Y]ode45('yang',tspan ,x0);
这里有问题
不过你的程序有问题吧

function chua
clear;
global A
range=[0.1: 0.0001: 0.35];
period=2*pi;
k=0;
YY1=[];
Step=2*pi/100;
for A=range
  x0=[0 0.001 0];
  A
  k=k+1;     %出去前面60个周期的数据，并将最后的结果作为下一次积分的初值。
  tspan=[0:step:60*period];
  [t,Y]=ode45('yang',tspan ,x0);
  y0=Y(end,:);
  j=1;
  for i=60:200
    tspan=[i*period:step(i+1)*period];
    [t,Y]ode45('yang',tspan ,x0);
    YY1(k, j)=Y(end,1);
    j=j+1;    %取出每一个周期内的第一个解的最后一个值。
    y0=Y(end,:);
   end
end
bifdata=YY1(:,end-51:end);
plot(range,bifdata,'k','markersize',1)；
为什么出现两次 [t,Y]=ode45('yang',tspan ,x0);

ymsncx

我是按照liliangbiao提供的Duffing系统的分岔图的Matlab实现中的程序编的，实现不了