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

小弟初学Matlab，按照公式敲一个泰勒线源综合的程序，但是仿真出来的结果不对啊，求教各位大神，程序哪里出错了呢？ @"|i"Hk^  
Mkv|TyC  
SLL=input('请输入SLL=');n=input('请输入n='); X-JV'KE}^z  
R0=10^(SLL/-20);A=(1/pi)*acosh(R0);n0=2*A^2+4; Z1,gtl ?  
x=ones(1,n);lamda=8.57;L=(n-1)*lamda/2; d=lamda/2;n1=6; lDc;__}Ws  
for i=1:n P $>`  
    if i<=n1-1 BC/_:n8O  
        x(i)=(n1)^2*(A^2+(i-1/2)^2)/(A^2+(n1-1/2)^2); Mm=Mz  
    else gPEqjj  
         x(i)=i^2; z?ucIsbR  
    end Jb(Y,LO^  
     }zMf7<C  
end |4b)>8TL/  
for j=1:n 8_we: 9A  
    if j<n1 '2^ Yw  
       for k=1:n1-1 \pY^^ l*  
           B(j,k)=1-j^2/x(k); #y; yN7W  
       end; #a9_
~\s  
       B1=B(j,:); ;IZ
?19Q  
       s(j)=(factorial(n1-1))^2/(factorial(n1-1+j)*factorial(n1-1-j))*prod(B1); 'Ud|Ex@A9  
    else 0 bSA_  
        s(j)=0; Hs?zq  
    end k#)
Ad*t  
end fCX*R"  
if mod(n,2)>0;     svhrf;3:  
    for o=1:(n-1)/2 8cVzFFQP  
        z(o)=o*d;z(n+1-o)=-o*d; &M"ouy Zo9  
    z((n-1)/2+1)=((n-1)/2+1)*d; fA^Em)cs2  
    end O\beKBT;  
else H[W
eu  
    for o=1:n/2 H\G{
3.T.9  
        z(o)=(2*o+1)/2*d;z(n+1-o)=-(2*o+1)/2*d; F $yO  
    end uV]ULm#,i
 
end P-Su5F  
for l=1:n [Ot,q/hBJ  
    for m=1:n1-1 (#Mp 5C'X  
        z0(m)=s(m)*cos(2*m*pi*z(l)/L); ,Og4 ?fS  
        y(l)=0+z0(m); RKkGITDk  
    end M##';x0  
    I(l)=1+2*y(l); ]~c+'E`  
end xtfRrX^  
BWq/TG=>  
rT7^-B*  
附正确结果：SLL=-20，n=19时，相应激励为 YD[H