求大神帮忙看看程序哪里有问题？

%计算均匀直线阵方向系数D随阵元N程序
clc;clear;
global n posai sita d;
sita=(0:pi/600:pi);k=2*pi;
nn=(2:19/100:20);d=0.25;
for jj=1:length(nn);
n=nn(jj);
beta=0;
%beta=k*d;
%beta=beta=-k*d-pi/n;
posai=beta+k*d*cos(sita);
jifen=0;
f=zxz(sita);
for i=1:length(sita);
jifen=jifen+f(i)*pi/600;
end;
fxxs(jj)=2/jifen;
end;
plot(nn,fxxs);hold on