求大牛帮忙怎么分岔图就是直线

wotaisa

tic;
uu=100:1:1000;
t0=0;tf=10;
hold on
for j=1:length(uu)
[t,y]=ode45('zf1',[t0,tf],[0,0,0,0,0,0,0,0,0,0,0,0],[],uu(j));
plot(uu(j),y(1:1:10,6),'linewidth',5)
title('bifurcation')
xlabel('参数 uf')
ylabel('w位移')
end
toc
zf1.m
function rdot=zf(t,y,flag,uu)
a1=-3.2341*10^11;a2=1/78.5*(10.511e11);a3=-1/78.5*(1.066e11);b1=1/78.5*(1.0511e12);b2=1/78.5*(5.8546e12);b3=-1/78.5*(1.1319e12);c1=1/497*(1.0611e11);c2=1/497*(1.1319e12);
c3=(2.277e11+141850*(uu)^2)/497;
rdot=[a1*y(4)+a2*y(5)+a3*y(6);b1*y(4)-b2*y(5)-b3*y(6);c1*y(4)-c2*y(5)+c3*y(6);y(1);y(2);y(3)];
end
运算后只是一条直线，没得到分岔啊。求大牛们帮忙看看。不胜感激。

wotaisa

zf1.m
function rdot1=zf1(t,y,flag,uu)
a1=-3.2341*10^11;a2=1/78.5*(10.511e11);a3=-1/78.5*(1.066e11);b1=1/78.5*(1.0511e12);b2=1/78.5*(5.8546e12);b3=-1/78.5*(1.1319e12);c1=1/497*(1.0611e11);c2=1/497*(1.1319e12);
c3=(2.277e11+141850*(uu)^2)/497;c4=15347*uu/497;d1=-4.2127e12/78.5;d2=2.1022e12/78.5;d3=-2.1325e11/78.5;e1=2.1022e12/78.5;e2=-6.4405e12/78.5;e3=1.132e12/78.5;f1=2.1325e11/162.269;f2=-1.132e12/162.269;
f3=(2.2857e11+237390*uu^2)/162.269;f4=-6420.5*uu/162.269;
rdot1=[a1*y(4)+a2*y(5)+a3*y(6);b1*y(4)-b2*y(5)-b3*y(6);c1*y(4)-c2*y(5)+c3*y(6)+c4*y(9);y(1);y(2);y(3);d1*y(10)+d2*y(11)+d3*y(12);e1*y(10)+e2*y(11)+e3*y(12);f1*y(10)+f2*y(11)+f3*y(12)+f4*y(3);y(7);y(8);y(9)];

wotaisa

plot(uu(j),y(1:1:10,6),'linewidth',5)
想让y坐标值更小，怎么设置？y(1:110)