T 'T 就是 N 的取值 范围 2<=T<=12 且必须为偶数
Lb10:{EG}:A”XA”:B”YA”:C”CA”:D”1/RA”:E”1/RB”:F”DKA”:G”DKB”
Lb1 1:{HOR}:H”DKI”:O”DL”:R”DR”:H ＞ G= ＞ Goto 3 Δ
P=(E-D)/ABS(G-F)
Q=ABS(H-F)
I=P×Q
J=C+(I+2D)Q×90/π
U=0:V=0:W=0:Z=0
S=T
Lbl 4:
U = U + Cos(C + (I * (2 * S - 1) / (2 * T) + 2 * D) * (Q / 2) * (2 * S - 1) / (2 * T)*180/π)
V = V + Sin(C + (I * (2 * S - 1) / (2 * T) + 2 * D) * (Q / 2) * (2 * S - 1) / (2 * T) *180/π)
DSZ S
Goto 4
S=T-1
Lbl 5
W = W + Cos(C + (I * S / T + 2 * D) * (Q / 2) * (S / T) *180/π)
Z = Z + Sin(C + (I * S / T + 2 * D) * (Q / 2) * (S / T) *180/π)
DSZ S
Goto 5
X=A+Q/(6T)*(cosC+4U+2W+cosJ) ▲
Y=B+Q/(6T)*(sinC+4V+2Z+sinJ) ▲
U”XL”=X+Ocos(J-90) ▲ V”YL”= Y+Osin(J-90) ▲
W”XR”=X+Rcos(J+90) ▲ Z”YR”=Y+Rsin(J+90) ▲
Goto 1
Lbl 3:A=X:B=Y:D=E:F=G:C=J:Goto 0
其他参数的意义，参见李少元编写的程序，在此不提及了。
我建议大家不要一味追求精度，具体工程采用具体的实际计算方法。当 N ＝ 12 时，在 Fx4800 上的计算比较慢的。如果不是那种缓和曲线特别长，圆曲线半径特别小的线形， N ＝ 2 就可以了。