matlab 运算出错

matlab 运算出错
function [w1,w2,VoA,VoB,VoC,VoD,VoE,VA1,VB1,VC1,VD1,VE1,VA2,VB2,VC2,VD2,VE2,AVA,AVB,AVC,AVD,AVE,RA1,RB1,RC1,RD1,RE1,RA2,RB2,RC2,RD2,RE2,CpA,CpB,CpC,CpD,CpE]=mesh(z1,z2);
m=2;
h=1;
T=6400;
n=1000/60;
i=2;
a=20;
b=8;
e=210000
r1=m*z1/2;
r2=m*z2/2;
ra1=r1+h*m;
ra2=r2+h*m;
a1=acos(r1*cos(a*pi/180)/ra1);
a2=acos(r2*cos(a*pi/180)/ra2);
w1=2*pi*n;
w2=w1*r1/r2;
p=T/(r1*cos(a*pi/180)*b);
CHD=(z1*(tan(a1*pi/180)-tan(a*pi/180))+z2*(tan(a2*pi/180)-tan(a*pi/180)))/(2*pi);
pb=pi*m*cos(a*pi/180);
AE=pb*CHD;
K=AE*(pb-1)/pb;
C1C2=z1+z2;
R2=ra1;
R1=ra2;
AC1=sqrt(R1^2+C1C2^2-2*R1*C1C2*cos((a2-a)*pi/180));
BC1=sqrt(R2^2+pb^2-2*R2*pb*sin(a1*pi/180));
C1C=r1;
C1D=sqrt(R2^2+K^2-2*R2*K*sin(a1*pi/180));
EC2=sqrt(C1C2^2+R2^2-2*C1C2*R2*cos((a1-a)*pi/180));
DC2=sqrt(pb^2+R1^2-2*pb*R1*sin(a2*pi/180));
CC2=r2;
BC2=sqrt(R1^2+K^2-2*R1*K*sin(a2*pi/180));
RA1=sqrt(AC1^2-(r1*cos(a*pi/180))^2);
N1N2=AE+RA1+sqrt(EC2^2-(2*r2*cos(pi/9)^2));
RA2=N1N2-RA1;
RB1=RA1+K;
RB2=N1N2-RB1;
RC1=sqrt(C1C^2-(r1*cos(a*pi/180))^2);
RC2=N1N2-RC1;
RD1=RA1+pb;
RD2=N1N2-RD1;
RE1=RA1+AE;
RE2=N1N2-RE1;
VA1=w1*AC1;
VA2=w2*R1;
VoA=VA2*cos(a*pi/180);
AVA=(w1+w2)*(K+RC1-RB1);
VB1=w1*BC1;
VB2=w2*BC2;
VoB=VB2*cos(a*pi/180);
AVB=(w1+w2)*(RC1-RB1);
VC1=w1*C1C;
VC2=w2*CC2;
VoC=VC2*cos(a*pi/180);
AVC=(w1+w2)*0
VD1=w1*C1D;
VD2=w2*DC2;
VoD=VD2*cos(a*pi/180);
AVD=(w1+w2)*(RD1-RC1);
VE1=w1*R2;
VE2=w2*EC2;
VoE=VE2*cos(a*pi/180);
AVE=(w1+w2)*(RD1-RC1+K);
CpA=0.418*sqrt(p*e*(RA1+RA2)/(RA1*RA2));
CpB=0.418*sqrt(p*e*(RB1+RB2)/(2*RB1*RB2));
CpC=0.418*sqrt(p*e*(RC1+RC2)/(RC1*RC2));
CpD=0.418*sqrt(p*e*(RD1+RD2)/(2*RD1*RD2));
CpE=0.418*sqrt(p*e*(RE1+RE2)/(RE1*RE2));
最上边方括号里面都是我想要计算的值 为什么赋值之后 结果会是下面的情况
[w1,w2,VoA,VoB,VoC,VoD,VoE,VA1,VB1,VC1,VD1,VE1,VA2,VB2,VC2,VD2,VE2,AVA,AVB,AVC,AVD,AVE,RA1,RB1,RC1,RD1,RE1,RA2,RB2,RC2,RD2,RE2,CpA,CpB,CpC,CpD,CpE]=mesh(21,42);
e =
210000
AVC =
0
只出现了两个值 并且e是已知的