求达人帮忙：为啥while嵌套while出不来哦

[liusanshao]

是哪里的问题呢
全部程序如下：
bc=20;g=9.81;y1=6;y3=[0.8 1.0 1.5 2.0 3.0 4.0 5.0];a=5.2;r=7.8;
w=[0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0];Q=zeros(8,7);y2=zeros(8,7);
Re11=zeros(8,7);Kc11=zeros(8,7);Ecorr1=zeros(8,7);M11=zeros(8,7);
Re13=zeros(8,7);Kc13=zeros(8,7);M13=zeros(8,7);
yw1=zeros(8,7);Fw1=zeros(8,7);F31=zeros(8,7);E0=zeros(8,7);Ecorr3=zeros(8,7);Re31=zeros(8,7);Kc31=zeros(8,7);M31=zeros(8,7);
Re33=zeros(8,7);Kc33=zeros(8,7);M33=zeros(8,7);yw3=zeros(8,7);Fw3=zeros(8,7);F33=zeros(8,7);EM=zeros(8,7);Cd=zeros(8,7);
M01=zeros(8,7);M03=zeros(8,7);Ecorr0=zeros(8,7);Re01=zeros(8,7);Re03=zeros(8,7);Kc01=zeros(8,7);Kc03=zeros(8,7);
o=acos((a-w)/r);p=0.643;
Ks=1.001-0.2349*o-0.1843*o.^2+0.1133*o.^3;
yj=Ks.*w;
for i=1:8
for j=1:7
%-----------------------------------------------------从判断一个闸门状态的自由或淹没出流开始求Q------------------------------------------
if y1>=0.81*y3(j)*(y3(j)/w(i))^0.72
Cd(i,j)=Ks(i)/sqrt(1+Ks(i)*w(i)/y1);
Q(i,j)=Cd(i,j)*w(i)*bc*sqrt(2*g*y1);
else
%-----------------------------------------------------求y2开始------------------------------------------------------------------------------
y21=0.001; y22=6.000;
Ecorr1(i,j)=(y21-yj(i))*(0.52-0.34*atan(7.98*((y21-yj(i))/yj(i))-0.83));
R1=bc*y1/(bc+2*y1);
%--------------------------------------------------------------开始求对应y21的流量----------------------
Q11=0.1; Q12=500.0;
Re11(i,j)=Q11*R1/(bc*w(i)*0.000001306);
Kc11(i,j)=0.15*exp((-0.000005)*Re11(i,j));
M11(i,j)=(Ks(i)*w(i)*bc*sqrt(2*g*(y1-y21+Ecorr1(i,j))/(1+Kc11(i,j))))/Q11-1;

while abs(Q11-Q12)>0.0001
Q13=0.5*(Q11+Q12);
Re13(i,j)=Q13*R1/(bc*w(i)*0.000001306);
Kc13(i,j)=0.15*exp(-0.000005*Re13(i,j));
M13(i,j)=(Ks(i)*w(i)*bc*sqrt(2*g*(y1-y21+Ecorr1(i,j))/(1+Kc13(i,j))))/Q13-1;
if M11(i,j)*M13(i,j)<0
Q12=Q11;
end
Q11=Q13;
M11(i,j)=M13(i,j);
end
Q1=Q13;
%--------------------------------------------------------结束求对应y21的流量---------------
yw1(i,j)=p*y21+(1-p)*y3(j);
Fw1(i,j)=0.5*g*bc*yw1(i,j)^2;
F31(i,j)=0.5*g*bc*y3(j)^2;
E0(i,j)=(Q1*sqrt((Q1/(bc*yj(i)))^2-2*g*Ecorr1(i,j))+0.5*bc*g*y21^2+Fw1(i,j))-(Q1*Q1/(bc*y3(j))+F31(i,j));
while abs(y21-y22)>0.001
y23=0.5*(y22+y21);
Ecorr3(i,j)=(y23-yj(i))*(0.52-0.34*atan(7.98*((y23-yj(i))/yj(i))-0.83));
R1=bc*y1/(bc+2*y1);
%-------------------------------------开始求对应y23的流量----------------------------
Q31=0.1;Q32=500.0;
Re31(i,j)=Q31*R1/(bc*w(i)*0.000001306);
Kc31(i,j)=0.15*exp((-0.000005)*Re31(i,j));
M31(i,j)=Q31-(Ks(i)*w(i)*bc*sqrt(2*g*(y1-y23+Ecorr3(i,j))/(1+Kc31(i,j))));
while abs(Q31-Q32)>0.001
Q33=0.5*(Q11+Q12);
Re31(i,j)=Q33*R1/(bc*w(i)*0.000001306);
Kc31(i,j)=0.15*exp(-0.000005*Re31(i,j));
M33(i,j)=Q33-(Ks(i)*w(i)*bc*sqrt(2*g*(y1-y23+Ecorr3(i,j))/(1+Kc31(i,j))));
if M31(i,j)*M33(i,j)<0
Q32=Q31;
end
Q31=Q33;
M31(i,j)=M33(i,j);
end
Q3=Q33;
%--------------------------------------结束求对应y23的流量---------------------------
yw3(i,j)=p*y23+(1-p)*y3(j);
Fw3(i,j)=0.5*g*bc*yw3(j)^2;
F33(i,j)=0.5*g*bc*y3(j)^2;
EM(i,j)=(Q3*sqrt((Q3/(bc*yj(i)))^2-2*g*Ecorr3(i,j))+0.5*bc*g*y23^2+Fw1(i,j))-(Q3*Q3/(bc*y3(j))+F33(i,j));
if E0(i,j)*EM(i,j)<0
y22=y21;
end
y21=y23;
E0(i,j)=EM(i,j);
end
y2(i,j)=y23;
%-----------------------------------------------------求y2结束------------------------------------------------------------------------------
%-----------------------------------------------------求Q开始--------------

Ecorr0(i,j)=(y2(i,j)-yj(i))*(0.52-0.34*atan(7.98*((y2(i,j)-yj(i))/yj(i))-0.83));
R1=bc*y1/(bc+2*y1);
Q01=0.1;Q02=500.0;
Re01(i,j)=Q01*R1/(bc*w(i)*0.000001306);
Kc01(i,j)=0.15*exp((-0.000005)*Re01(i,j));
M01(i,j)=(Ks(i)*w(i)*bc*sqrt(2*g*(y1-y2(i,j)+Ecorr0(i,j))/(1+Kc01(i,j))))/Q01-1;
while abs(Q01-Q02)>0.001
Q03=0.5*(Q01+Q02);
Re03(i,j)=Q03*R1/(bc*w(i)*0.000001306);
Kc03(i,j)=0.15*exp(-0.000005*Re03(i,j));
M03(i,j)=(Ks(i)*w(i)*bc*sqrt(2*g*(y1-y2(i,j)+Ecorr0(i,j))/(1+Kc03(i,j))))/Q03-1;
if M01(i,j)*M03(i,j)<0
Q02=Q01;
end
Q01=Q03;
M01(i,j)=M03(i,j);
end
Q(i,j)=Q03;
%-----------------------------------------------------求Q结束------------------------------------------------------------------------------
end
%-----------------------------------------------------结束这一一个闸门状态的自由或淹没出流求Q，进入下一个状态------------------------------
end
end

先谢谢各位回答我这个傻傻的问题了

[laosam280]

缺少了end

[liusanshao]

作者: laosam280

缺少了end

没有少吧
敢问在哪个地方加呢