Deg：Fix4 定义四位有效小数
“XJ=”？W 输入右夹角，垂直为90°、三都河为75° Lbl 1：“DKI=”？H：“L=”？O：Prog“ROAD-DATA1” 第一循环开始，输入变量桩号与边距，进入子程序调取参数
（E-D）÷Abs（G-F）→P：Abs（H-F）→Q：P×Q→I 计算参数P、Q、I
“J=”：C+(I+2D) ×Q×90÷π→J◢ 计算切线方位角并显示方位角 C+（I÷4+2D）Q×45÷(2π) →M：C+（3I÷4+2D）Q×135÷(2π) →N 计算参数M、N C+（I÷2+2D）Q×45÷π→K 计算参数K
A+（Q÷12）（cosC+4(cosM+cosN)+2cosK+cosJ）→X：B+（Q÷12）（sinC+4(sinM+sinN)+2sinK+sinJ）→Y： “XL=”：X+Ocos（J+W）→U◢
“YL=”：Y+Osin（J+W）→V◢
计算中桩坐标，显示相应边桩坐标 Goto1 循环结束
“ROAD-DATA1” 数据库子程序
84994.506→F：86081.185→G：If H≤G：Then 3094440.1858→A：498527.6831→B：264°45′56.83″→C：0→D：0→E：Return：IfEnd
上述的意思是，定义84994.506为起点桩号，86081.185为终点桩号，假如输入的桩号H小于等于86081.185时，则3094440.1858为起点的X坐标，498527.6831为起点的Y坐标，264°45′56.83″为起点的切线方位角，0为起点的曲率，0为终点的曲率，此段结束。（说明这段为直线上面）
G→F：86326.185→G：If H≤G：Then3094341.0511→A：497445.5354→B：264°45′56.83″→C：0→D：-1÷2000→E：Return：IfEnd
上述的意思是，此段的起点为上一段在终点，86326.185为该段的终点桩号，假如输入的桩号H小于等于86326.185时，则3094341.0511为起点的X坐标，497445.5354为起点的Y坐标，264°45′56.83″为起点的切线方位角，0为起点的曲率，-1÷2000为终点的曲率，此段结束。（说明该段位于缓和曲线上面，该段由直线接到一个半径为2000的圆曲线上面。该圆曲线是往路线的左边转的，假如向右转则为正的曲率） G→F：86823.097→G：If H≤G：Then3094313.729→A：497202.105→B：261°15′23.11″→C：-1÷20000→D：-1÷2000→E：Return：IfEnd
上述的意思是，此段的起点为上一段在终点，86823.097为该段的终点桩号，假如输入的桩号H小于等于86823.097时，则3094313.729为起点的X坐标，497202.105为起点的Y坐标，261°15′23.11″为起点的切线方位角，-1÷2000为起点的曲率，-1÷2000为终点的曲率，此段结束。（说明该段位于半径为2000、并向左转圆曲线上面，上述的起点坐标和方位角可以按照前一段计算出来，及为前一段的终点坐标和方位交） „„„„
直到把本标段的桩号编辑完毕，在此就不一一写出来了
在程序使用时，只需要输入待算点桩号，边距，就能算出准确的算出坐标。要算左边的则把边距输入负的边距，右边反之，当边距输入0时，则为中桩坐标。
任意点标高计算程序的扩充程序是为了对路基精加工编辑的 主程序“BGJS”
“GC=”？A 输入距设计路面的理论高差（见附表1）
Lbl 1：“DKI=”？X：Prog“SQXSJK” 循环开始，输入桩号，进入中桩标高数据库调取参数 AI（T-Abs(Q-X)）²÷（2R）→Y：B（O-T×M）→H：H+（T-Abs（Q-X））×M+K+Y→Z 计算设计中桩标高
“D=”？D：“Z=”？P 输入边距和测量标高 If D＞0：Then Goto 0：Else Goto 2：IfEnd 判断左边还是右边
Lbl 0：D-1→S：Prog “YBPB-HP” 执行右边计算 扣除中间分离带1m，并进入右边坡比数据库（超高） “CZ=”：Z-V×S-A-P→W◢ 计算测量标高距94区顶标高。若为正，则为填土不够。若为负则反之 Lbl 2：AbsD-1→S：Prog “ZBPB-HPB”执行左边计算 扣除中间分离带1m，并进入左边边坡比数据库（超高）
“CZ=”：Z-V×S-A-P→W◢ 计算测量标高距94区顶标高。若为正，则为填土不够。若为负则反之 Goto 1 循环结束 （若算其他填筑区和挖方区，则把1.55改为其他相应的高度） 中桩标高数据库程序“SQXSJK”
85380→F：85500→G：If X≤G：Then0.005→M：2000→R：1→I：1→B：120→T：213.538→O：85500→Q：0→K：Return：IfEnd
上述的意思是：当桩号在85380到85500之间时，该段是在凹曲线的变坡点之前，在任意中桩点标高的计算为（T-Abs(Q-X)）²÷（2R）=Y、（O-T×M）=H、H+（T-Abs（Q-X））×M +Y=Z 该段的M为在该段竖曲线前的坡比，即为0.005
G→F：85620→G：If X≤G：Then-0.02→M：2000→R：1→I：1→B：120→T：213.538→O：85500→Q：0→K：Return：IfEnd
上述的意思是：当桩号在85500到85620之间时，该段是在凹曲线的变坡点之后，在任意中桩点标高的计算为（T-Abs(Q-X)）²÷（2R）=Y、（O-T×M）=H、H+（T-Abs（Q-X））×M +Y=Z 该段的M为在该段竖曲线后的负坡比，即为-0.02
G→F：85980→G：If X≤G：Then0.02→M：9000→R：-1→I：1→B：180→T：219.538→O：85800→Q：0→K：Return：IfEnd
上述的意思是：当桩号在85620到85980之间时，该段是在凸曲线上，在任意中桩点标高的计算为（T-Abs(Q-X)）²÷（2R）=Y、（O-T×M）=H、H+（T-Abs（Q-X））×M +Y=Z 该段的特点是曲线线前后直线的坡比的绝对值相等，符号相反。因此不曾在坡比变化，M都是0.02
G→F：86348→G：If X≤G：Then0.02→M：0→B：86348→T： 0→Q：208.578→K：Return：IfEnd 上述的意思是：当桩号在85980到86348之间时，该段位于直线上面，在任意点中桩标高计算则直接为坡比乘以距离加上起算点的标高 „„„.
综上所述：当位于竖曲线上时， B=1、K=0
在凹曲线上时，I=1、 在凸曲线上时，I=-1
在变坡点前时，M=变坡前的坡比 在变坡点之后时，M=负的边坡后的坡比
在没有竖曲线时， B=0、Q=0、T=终点桩号、K=终点标高
竖曲线长度计算：L=R(M-N) 式中R为竖曲线半径，M为前坡比，N为后坡比，对于编辑竖
曲线桩号有一定的作用
“YBPB-HP”右边超高要素程序
85560→E：86081.185→C：If X≤C ：Then 0.02→V：Return：IfEnd
C→E：86196.185→C：If X≤C ：Then 0.02-（4÷11500）×（X-86081.185）→V：Return：IfEnd „„„„„
其中0.02-（4÷11500）×（X-86081.185）的有得为：按照超高的渐变型式，建立数学模型得到
左边超高要素子程序编程方式相同，就不在此介绍了
注意：在操作此程序时要注意的是，边距输入为左负右正，当边距小于1的时候，输入1或者-1（原因为中央分隔带没有坡比），
任意点桩号计算程序“ZHJS”
主程序“ZHJS”
Deg：Fix3：10→Dimz 定义为三位小数，扩充10为变量 Lbl 1 第一循环开始
“N=”？S：“E=”？T：“KP=”？H 输入变量全站仪显示的N E坐标，输入大致桩号KP （减少计算器的运行量，减少运算时间）
Lbl 2 第二循环开始 Prog “XXGS” 进入坐标运算程序
J-90→R 计算垂直切线方位交的直线方位角 （T-Y）cosR-（S-X）sinR→W 主线上该点距起点的距离
If AbsW ＞0.0001：Then H+W→H：Goto 2：IfEnd 满足上述条件后W再加上大致桩号，进入上述循环计算，把这次循环的H做为下次计算的大致桩号，当桩号重合后进入下面运算 “DKI=”：H◢ 显示桩号
“D=”：（Y-T）÷sin R →Z[2] ◢ 计算并显示距中桩之间的距离（左负右正） Goto 1 循环完成
“XXGS”坐标计算程序
Prog“ROAD-DATA1” 进入子程序调取参数 （E-D）÷Abs（G-F）→P：Abs（H-F）→Q：P×Q→I 计算参数P、Q、I
“J=”：C+(I+2D) ×Q×90÷π→J◢ 计算切线方位角并显示方位角 C+（I÷4+2D）Q×45÷(2π) →M：C+（3I÷4+2D）Q×135÷(2π) →N 计算参数M、N C+（I÷2+2D）Q×45÷π→K 计算参数K
A+（Q÷12）（cosC+4(cosM+cosN)+2cosK+cosJ）→X：B+（Q÷12）（sinC+4(sinM+sinN)+2sinK+sinJ）→Y 计算中桩坐标
“ROAD-DATA1”数据库
跟坐标第一节介绍的一样，不作介绍（为了简单输入，可以只把坐标编写到千米位）
还有两个4800计算器计算程序（防止在半径较小的曲线中构筑物出现喇叭形），1、 C “X1”：D “Y1”：E “X2”：F“Y2”：Fixm：Pol（E-C，F-D：I“S1-2=” ◢ J≤0 ＝＞ J=J+360⊿ J “A1-2=” ◢
2、 C “X1”：D “Y1”：S“L”：R“A”：Fixm ：X“X2”=C+Rec（S，R）◢ Y “Y2”=D+J◢
因篇幅问题不能全部显示，请点此查看更多更全内容