求助

xzy919

<p>qb=1<br/>jj=5<br/>s=50<br/>fx=0<br/>x1=1<br/>x2=1<br/>i=1<br/>pi=3.1415<br/>gamn=45*pi/180<br/>gamm=135*pi/180<br/>fa0=acos(((qb+x1)*(qb+x1)-x2*x2+jj*jj)/(2*(qb+x1)*jj))<br/>pu0=acos(((qb+x1)*(qb+x1)-x2*x2-jj*jj)/(2*x2*jj))<br/>fai=fa0+0.5*pi*i/s<br/>pui=pu0+2*(fai-fa0)*(fai-fa0)/(3*Pi)<br/>ri=sqrt(qb*qb+jj*jj-2*qb*jj*cos(fai))<br/>alfi=acos((ri*ri+x2*x2-x1*x1)/(2*ri*x2))<br/>bati=acos((ri*ri+jj*jj-qb*qb)/(2*ri*jj))<br/>g1=x1*x1+x2*x2-(jj-qb)*(jj-qb)-2*x1*x2*cos(gamn)<br/>g2=x1*x1-x2*x2+(jj+qb)*(jj+qb)-2*x1*x2*cos(gamn)<br/>*if,fai,gt,0,then<br/>&nbsp;*elseif,fai,le,pi<br/>&nbsp;psi=pi-alfi-bati<br/>&nbsp;*elseif,fi,gt,pi,lt,2*pi<br/>&nbsp;psi=pi-alfi+bati<br/>*endif<br/>fx=fx+(pui-psi)*(pui-psi)</p><p>帮我看一下程序有没有错误</p>