m(meas_point)=macro/x,y,z
f(po01)=feat/point,cart,x,y,z,0,0,1
meas/point,f(po01),1
ptmeas/cart,x,y,z,0,0,1
endmes
output/fa(po01)
endmac
DECL/INTGR,i
decl/double,x,y,z
decl/double,Rz,e,rz,k1
Rz=25
e=1.25
k1=0.5
do/i,0,360,0.1
x=Rz*sin(i/180*3.1415926)-e*sin(Zb*i/180*3.1415926)
y=Rz*cos(i/180*3.1415926)-e*cos(Zb*i/180*3.1415926)
call/m(meas_point),x,y,10
enddo
上面程序可运行,并输出相应结果
但在其x=Rz*sin(i/180*3.1415926)-e*sin(Zb*i/180*3.1415926)式子后加上2.5*(k1*sin
(Zb*i/180*3.1415926)-sin(i/180*3.1415926))/sqrt(1+k1*k1-2*k1*cos(9*i/180*3.1415926))
后运行,系统提示这句出了问题,请问各位这是怎么回事
谢谢!