分岔图程序出图对称

Arno

function ydot=xiezhenqi(t,y,flag,vb)
b=5e-8;
k1=5;
k3=15e12;
m=5e-12;
d=2e-6;
c0=1.875e-18;
w=0.5;
w0=sqrt(k1/m);
vac=0.21;
F=0.5*c0*(vb+vac*sin(w*w0*t)).^2/(d-y(1)).^2-0.5*c0*vb^2/(d+y(1)).^2;
ydot=[y(2)
    (-b*y(2)-k1*y(1)-k3*y(1).^3+F)/m];
end


clear
vb=0:0.1:10;
c=0.5;
T=2*pi/c;
hold on
for j=1:length(vb)
   [t,y]=ode45('xiezhenqi',[0:T/10:T],[0,0],[],vb(j));
    plot(vb(j),y(:,2),'linewidth',5)
end
不知道我的这个程序问题出在哪里？  分岔图运行出的图关于x轴对称。。。  

  对于频闪法画分岔图还望各位大大请教！！