这种椭圆该怎么用宏程序编？

不完整椭圆的结构尺寸如图所示，椭圆的方程由图可知为X2/122+Z2/182=1。从图中可知自变量2的起始坐标分别为A（-10）和B（-39）。Z轴方向上椭圆的原点O1与编程原点O不重合，两者相差10+14.5=24.5。椭圆图形的起点为Z=14.5（自变量初始值），程序循环的起点为A，坐标ZA=14.5-24.5=-10。椭圆图形的终点为Z=-14.5（自变量终点值），程序循环的终点为B坐标ZB=-14.5-24.5=-39。

编程时把Z作为自变量，由椭圆的方程可知因变X2=12×√（1-Z2/182）。从图中可知X轴方向上的原点O1与编程原点O不重合，两者相差12mm，所以在指令机床移动时，X的坐标为X=2×（12+X2）。编制的椭圆加工程序如下：

N1M03S800T0101

N2G00X54Z3

N3G71U1R0.5P5Q16X0.5Z0F100

N4M03S1000T0101

N5G00G42X0

N6G01Z0

N7G01X20

N8Z-10

N9#1=14.5

N10WHILE#1GE[-14.5]

N11#2=12*SQRT[1-#1*#1/[18*18]]

N12G01X[24+2*#2]Z[#1-24.5]

N13#1=#1-0.5

N14ENDW

N15W-8

N16G40X54

N17G00Z200

N18M30