由本人自网络收集,在此谢谢原作者!

主:
PQX Lbi 0：prog“TAD”Goto 1↙
Lbi 1:Z[1]“ZH-X”=M+Tcos（F+180）:Z[2]“ZH-Y”=N+Tsin（F+180）：V[1]“HZ-X”=M+Tcos（F+Z）：V[2]“HZ-Y”=N+Tsin（F+Z）：H=0=>Goto 2:≠>D=H:prog“QXY”:A“HY-X”= Z[1]+I cos(F+SZ÷AbsZ) ：C“HY-Y”= Z[2]+I sin(F+SZ÷AbsZ):Q“YH-X”= V[1]+I cos(F+Z+180-SZ÷AbsZ ) ：E“YH-Y”= V[2]+ Isin(F+Z+180-SZ÷AbsZ )↙
Lbi 2:K≤G=>D=Abs(G-K):O=M:P=N:I=D+T:J=F+180:W=F:Goto 4⊿↙ H=0→Goto 3⊿↙ K≤G+H=>D=Abs(K-G):O=Z[1]:P=Z[2]:prog“QXY” :J=F+SZ÷AbsZ:W=F+3SZ÷AbsZ: Goto 4⊿↙ Lbi 3:K≤G+L-H=>D=Abs(K-G-H):H=0=>H=1:O=Z[1]:P=Z[2]:≠>H=H:O=A:P=C:prog“QXY”H=1=>Q=0⊿I=D:J=F+(Q+E)Z÷AbsZ:W=F+(Q+2E)Z÷AbsZ: Goto 4⊿↙
K≤V=>D=Abs(V-K):O=V[1]:P=V[2]:prog“QXY”:F=F+Z+180: :J=F-SZ÷AbsZ:W=F-3SZ÷AbsZ:U=180+U: Goto 4⊿↙
Lbi 4:X=O+IcosJ+Bcos(W+U)▲ Y=P+IsinJ+Bsin(W+U)▲ I=0：J=0：I=POL（X-x0，Y-y0：J<0→J=J+360△J“FWJ”=J▲I“r”=I▲ （x0、y0为置仪点坐标）
子:QXY
A=D-D^5/40R²H²+ D^9/3456R^4H^4- D^13/599040R^6H^6+ D^17/17542600R^8H^8+1E-99:C=D³/6RH-D^7/336R³H³+D^11/42240R^5H^5-D^15/976800R^7H^7+ D^19/3530097000R^9H^9: S=tan¯1(C÷A):I=√(A²+C²):Q=90H÷лR:E=90D÷лR:D=2RsinE
子:TAD
Lbi 0:{KBU}:KBU↙
K≤xxx=>M=xxx:N=xxx: F=xxx: Z=xxx: R=xxx: H=xxx: T=xxx: L=xxx: G=xxx: V=xxx:≠>K≤xxx=>M=xxx:N=xxx: F=xxx: Z=xxx: R=xxx: H=xxx: T=xxx: L=xxx: G=xxx: V=xxx:≠>K≤xxx=>M=xxx … … … … … ⊿
注:
K：输入计算点桩号;
B：输入距中桩距离,左－,右＋;
U：输入与前进方向偏角计算中桩坐标时B、U输入0,M－交点X坐标;N－交点Y坐标;
F－路线前进方向方位角;
Z－转角,左转时连同“-”号输入;
R－半径;H-缓和曲线长,无缓和曲线时为0;T－切线长；
L－曲线总长;
G－直缓点桩号;
V－缓直的桩号在Lbi1行处,把各主点后的“:”改为“▲”，如Z[1]“ZH-X”=M+Tcos（F+180）:改成Z[1]“ZH-X”=M+Tcos（F+180）▲可直接显示主点桩号坐标,在子程“TAD”内可输入全线角桩资料,K≤本转角缓直点桩号时输入本转角曲线要素。

程序有些繁锁,重复语句多,对于它的准确性和可靠性望大家给予验证,修正.fengdisheng123@126.com