请高手看看这分岔程序

earth

clear all
global c
zhang=[];
M=[0:0.01:50];
counter=1;
for counter=1:length(M)
    c=M(counter);
    i=2;
xmax=0;
xmaxold=0;
frmdata=[];
error=0;
tspan=[0 300];
var=1;
y0=[0.1;0.1;0.1;0.1];
[t,x]=ode45(@chaoLu2,tspan,y0);
while i < size(x,1)
       if x(i-1,var) < x(i,var) & x(i+1,var) <= x(i,var)
            xmax=x(i,var);
            if xmaxold ~= 0
                frmdata=[frmdata ; xmax xmaxold];
            end
                xmaxold=xmax;
       end
        i=i+1;
end
r= length(frmdata)-20:length(frmdata);
fradata1(1,r)=frmdata(r,1);
zhang=[zhang;fradata1];%这个地方老是显示错误，这是错误信息{All rows in the bracketed expression must have the same
number of columns.}%
end

plot(M,zhang,'k.','markersize',1);
xlabel(sprintf('c'));
ylabel(sprintf('x'))

function dx=chaoLu2(t,X)

x=X(1);

y=X(2);

z=X(3);
w=X(4);

global c;
dx=zeros(4,1);

dx(1)=-35*(x-y-yz);

dx(2)=-x*z+w+14*y;

dx(3)=-3*z+y*x;
dx(4)=-c*x
请高手能给我指导下，谢谢！

咕噜噜

搜索论坛，用最大值法求分叉图
你这个就是最大值法

无水1324

是程序出了问题还是什么？