求大神看看我的程序哪里不对
 

小女菜鸟一枚，只懂皮毛。这是根据立铣刀铣削力模型编出来的程序，要出铣削力的仿真图形。一直调不出来，望各位大神赐教！！！:o:o:o:o
clear
ap=3;ae=8;
c=0.1;
n=10000;
D=16;N=2;beta=pi*30/180;
fi_st=0;fi_ex=pi/3;
Ktc=675.5;Krc=277.9;Kac=195.8;
Kte=6.6;Kre=0.43;Kae=7.5;
deta_fi=pi/100;
deta_ap=0.01;
fi_p=2*pi/N;
K=2*pi/deta_fi;
L=ap/deta_ap;
for i=1:K
Fx(i)=0;Fy(i)=0;Fz(i)=0;Ft(i)=0;fi(i)=0;
fi(i)=fi_st+i*deta_fi;
for k=1:N;
fi_1=fi(i)+(k-1)*fi_p;
fi_2=fi_1;
for j=1:L;
a(j)=j*deta_ap;
fi_2=fi_1-2*tan(beta)*a(j)/D;
if fi_st<=fi_2<=fi_st;
h(j)=c*sin(fi_2);
deta_Ft=deta_ap*(Ktc*h(j)+Kte);
deta_Fr=deta_ap*(Kre*h(j)+Kre);
deta_Fa=deta_ap*(Kac*h(j)+Kae);
deta_Fx=-deta_Ft*cos(fi_2)-deta_Fr*sin(fi_2);
deta_Fy=deta_Ft*sin(fi_2)-deta_Fr*cos(fi_2);
deta_Fz=deta_Fa;
Fx(i)=Fx(i)+deta_Fx;
Fy(i)=Fy(i)+deta_Fy;
Fz(i)=Fz(i)+deta_Fz;
else ap>ap_st;
ap=ap+deta_ap;
end
F(i)=(Fx(i)^2+Fy(i)^2)^(1/2);
Tc(i)=Ft(i)*D/2;
end
end
end
plot(i,Fx);
plot(i.Fy);
plot(i.Fz);
plot(i,Tc)