[讨论]求大神看看泰勒线阵综合的程序哪里出错了

小弟初学Matlab，按照公式敲一个泰勒线源综合的程序，但是仿真出来的结果不对啊，求教各位大神，程序哪里出错了呢？ e,EK,,iY5  
ovo/!YJ2  
SLL=input('请输入SLL=');n=input('请输入n='); 5QAdcEcN@O  
R0=10^(SLL/-20);A=(1/pi)*acosh(R0);n0=2*A^2+4; y
B.G
=90  
x=ones(1,n);lamda=8.57;L=(n-1)*lamda/2; d=lamda/2;n1=6; >L^xlm%7o  
for i=1:n QS^~77q  
    if i<=n1-1 +fM&su=wl  
        x(i)=(n1)^2*(A^2+(i-1/2)^2)/(A^2+(n1-1/2)^2); BYS lKTh  
    else 0|ZVA+  
         x(i)=i^2; {{32jU7<  
    end RoX &+~  
     RL6Vkd?  
end q%>7L<r  
for j=1:n Hk@r5<{  
    if j<n1 7skljw(  
       for k=1:n1-1 ,I ZqLA  
           B(j,k)=1-j^2/x(k); 4:D:|r  
       end; J7:9_/e0T  
       B1=B(j,:); 7!pLK&_  
       s(j)=(factorial(n1-1))^2/(factorial(n1-1+j)*factorial(n1-1-j))*prod(B1); y=w`w>%  
    else rOW;yJ[  
        s(j)=0; 3 -tO;GKb  
    end R<|ejw  
end D z@1rc<B  
if mod(n,2)>0;     W@^J6sH  
    for o=1:(n-1)/2 fmFzW*,E  
        z(o)=o*d;z(n+1-o)=-o*d; sm1;MF]/u  
    z((n-1)/2+1)=((n-1)/2+1)*d; @M=$qO_$9  
    end {=3B)+N  
else P'sfi>A  
    for o=1:n/2 !n eo\  
        z(o)=(2*o+1)/2*d;z(n+1-o)=-(2*o+1)/2*d; 8"x9#kyU<3  
    end #EJP(wXa  
end R"Kz!NTB  
for l=1:n ur6e&bTp  
    for m=1:n1-1 RwW$O@0  
        z0(m)=s(m)*cos(2*m*pi*z(l)/L); FqwIJ|ct  
        y(l)=0+z0(m); :s"2Da3B  
    end w
FvT0  
    I(l)=1+2*y(l); %%&e"&7HE  
end #1[z;Mk0  
8.yCA  
HMV)U
{  
附正确结果：SLL=-20，n=19时，相应激励为 @TysXx