谁帮我看看这个程序哪里错了？？

I=2300;
    H=5290;
    G=2020;
R=1950;
P=3350;
C=2400;
A=3000;
K=sqrt(I^2+(H-G)^2)
a=asin(I/K)
n=980/31.886
w=2*pi*n/60
b=0:2*pi/100:2*pi;
b2=2*pi-b+a
L=sqrt(R*R+K*K-2*R*K.*cos(b2))
c=asin(R./L.*sin(b2))
b3=acos((P*P+L.*L-C*C)/2/P./L)-c
b4=acos((P*P-L.*L-C*C)/2/C./L)-c
x=acos((C*C+L.*L-P*P)/2./L/C)
m=R.*b2./P.*sin(b4-b2)./sin(b3-b4)
n=R/C.*b2.*sin(b3-b2)./sin(b3-b4)
f=n.*((-w-m).*cot(b2-b3)-(m-n).*cot(b3-b4))
g=n*A
h=f*A
Emax=acos((C*C+K*K-(R+P)*(R+P))/2/C/K)
Emin=acos((C*C+K*K-(P-R)*(P-R))/2/C/K)
S=(Emax-Emin)*A