公路逐桩坐标计算程序
(可以计算对称、不对称缓和曲线)
Lb1 0
Z=？  
V=？
W=V+2：Fixm
｛K｝
Lb1 1
K>Z[W+5Z+4]=>W=W+1:Goto 1⊿    (判断桩号在哪个交点范围，就是该交点曲线起点至下一交点曲线起点)
S=K-Z[W+5Z+3]                              （计算该桩号与曲线起点的距离） 
R=Z[W+2Z+2]：L=Z[W+3Z+2]:E=Z[W+4Z+2]    （读取该交点曲线要素R、Ls1 、Ls2）
Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]）       （计算该交点与下一交点直线方位角）
J<0=>J=J+360⊿
A=J
Pol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]）      （计算该交点与上一交点直线方位角）
J<0=>J=J+360⊿
C=A-J:A=J                            （计算偏角）
W=V+2=>Goto2⊿                     （如果桩号在起点与第一交点曲线起点之间，则转 Lb1 2 ）
I=Abs（tan（c÷2））
M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2
P=L^2÷6R-L^4÷336R^3-R（1-cos（90L÷πR））
Q=E^2÷6R-E^4÷336R^3-R（1-cos（90E÷πR））
D=（P-Q）I÷2 : F=（P+Q+2R）I÷2
M=F+M-D:Q=F+N+D
N=πRAbsC÷180+（L+E）÷2
X=Z[W-1]-McosA
Y=Z[W+Z+1]-MsinA
M=Z[W-1]+Qcos（A+C）
V=Z[W+Z+1]+Qsin（A+C）       
Q=AbsC÷C
S≤L=>P=0:Goto3⊿    （如果桩号在第一缓和曲线内，则转 Lb1 3）
S≤N-E=>S=S-L:Goto4⊿（如果桩号在圆曲线内，则转 Lb1 4）
S≤N=>S=N-S
Q=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3 ⊿   （如果桩号在第二缓和曲线内，则转 Lb1 3）
P=A+C:S=S-N:D=M+ScosP:F=V+SsinP
Goto6                                           （如果桩号在直线内，则转 Lb1 6）
Lb1 2
P=A+C
D=Z[W-1]+ScosP
F=Z[W+Z+1]+SsinP:Goto6
Lb1 3
I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4
J=Q（S^3÷6RL-S^7÷336R^3÷L^3）
P=P+A+90QS^2÷πRL:Goto5
Lb1 4
M=90（2S+L）÷πR
I=RsinM+L÷2-L^3÷240R^2
J=Q（L^2÷24R+R（1-cosM））
P=A+QM
Lb1 5
D=X+IcosA-jsinA:F=Y+JcosA+IsinA 
Lb1 6
D″X=″◢                （结果显示X坐标）
F″Y=″◢                （结果显示Y坐标）
P″AT=″◢                 （结果显示该桩号方位角）
{BO}：B″S″O″⊿″      （输入边桩距离，交角）
P=P+O
L″XB″=D+BcosP◢             （结果显示边桩X坐标）
M″YB″=F+BsinP◢             （结果显示边桩Y坐标）
以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。
     该程序是最基本的，如再加几条语句就可以处理断链。
    该程序还需变量输入，把整条线路曲线要素一次输入到4800中，就可以计算线路任意中、边桩坐标。
如有需要，需RMB ￥200，附有变量输入方法及详细使用说明。联系：13858423737