平曲線坐標、反算樁號計算程序
prog "PQX"
{mn}:M"X0":N"Y0":Defm 1:Z[1]=0
Lbl 1:{ABGHVWK}:K"JL":G"JX":H"JY":V"I0":W"J0":A"R":B"LS":Fixm
P=B^2/24/A:Q=B/2-B^3/240/A^2:T=(A+P)tan(Abs W/2)+Q:C=K-T:D=C+B:F=D+π*A*AbsW/180:E=F-B
w<0=>s=-1:≠=>s=1△
Lbl 2:{L}:L:L=0=>GOTO 1△
Lbl 6:L≤C=>O=K-L:R=G-OcosV:U=H-OsinV:Z=V:GOTO 3△
L≥F=>O=L-F+T:Z=V+W:R=G+OcosZ:U=H+OsinZ:GOTO 3△
L<D=>O=L-C:prog "XY":Z=V+SR:R=G+(I-T)cosV-SJsinV:U=H+(I-T)sinV+JScosV:GOTO 3△
L>E=>O=F-L:prog "XY":O=V+W:Z=O-SR:R=G+(T-I)cosO-SJsinO:U=H+(T-I)sinO+SJcosO:GOTO 3△
O=180(L-D+B/2)/π/A:I=AsinO+Q:J=P+A(1-cosO):Z=V+SO:R=G-(T-I)cosV-SJsinV:U=H-(T-I)sinV+SJcosV
Lbl 3:Z[1]=0.002=>GOTO 5:≠=>O=0△R"XL"◢U"YL"◢Z"ZL"◢
Lbl 4:{O}:O"YC":O=0=>GOTO 2△
O=-1=>GOTO 5△
X=R-OsinZ:Y=U+OcosZ:pol(X-M,Y-N:J<0=>J=J+360△
X"XZ"◢Y"YZ"◢J"A0"◢I"DD"◢GOTO 4
Lbl 5:{XY}:X"XF":X=-1=>Z[1]=0:GOTO 4△Y"YF":pol(X-R,Y-U:O=Icos(J-Z):AbsO≤0.001=>O"DL"◢L+O◢O=Isin(J-Z):O"YC"◢GOTO 5:≠=>Z[1]=0.002:L=L+O:GOTO 6△
prog "XY"
U=AB:I=O-O^5/40/U^2:J=O^3/6/U-O^7/336/U^3:R=90O^2/π/U
變量名稱說明：
M"X0":N"Y0"為測站點坐標
K"JL":G"JX":H"JY":V"I0":W"J0":A"R":B"LS"分別為交點樁號、坐標、起始方位角、轉角(左-右+)、半徑、緩長
L、R"XL"、U"YL"、Z"ZL" 分別為輸入樁號、求出中樁坐標、即時方位角
O"YC"：X"XZ"◢Y"YZ"◢J"A0"◢I"DD" 分別為輸入邊樁距離(左-右+)、求出邊樁坐標、邊樁點到測站點的方位角和平距
O"YC"=0返回計算中樁坐標
O"YC"=-1進入計算樁號過程(輸入坐標X"XF" 、Y"YF"分別為求出樁號L+O,及至中樁邊距O"YC"，O"DL"為樁號計算的誤差距)
O"YC"=約0.00001可求中樁點到測站點方位角和平距
X"XF"=-1返回計算坐標過程
計算樁號過程時，需多次確認輸入的坐標X"XF"、Y"YF"直到出現O"DL"樁號誤差距較小時即可得出準確的樁號
=>前不用冒號
≠=>前有冒號
△前后不用冒號