第二次修改：有关路基和隧道的CASIO 5800P 万能曲线 计算程序，可以正、反坐标，放边坡和隧道圆心，计算开挖半径
曲线任意里程中边桩坐标正反算
(CASIO fx-5800P计算器)程序
一、程序功能
本程序由一个主程序(ZHUCHENXU)和几个子程序——正算子程序(SUB1)、反算子程序(SUB2)、数据库(SUB3)、算方位角程序(FA)、使用方位角算坐标(NE)、边坡放样程序(BIANPO)、隧道圆心放样程序（SDYX）构成，可以根据直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。
二、源程序
增加变量22→DimZ
1.主程序(ZHUCHENXU)
“1.SZ → XY，2.XY → SZ，3.SDYX，4. BIANPO “？N：
Lbl 1：”ZHUANGHAO=”?S： Prog “SUB3″： N≥2=>Goto 2： Abs(S-O)→W：”BIANZHU=”?→Z:Prog “SUB1″： “XS=”:X→X◢”YS=”:Y→Y◢”FS=”:(F-90) →F:F►DMS◢Prog “FA”：Goto 1:
Lbl 2：”CX=”?X: X→I：”CY=”?Y: Y→J：”DMG=”?→Z[7]：Prog “SUB2″： “S=”:(O+W) →S◢”Z=”：Z→Z◢If N=3:Then Prog”SDYX”: IfEnd ：If N=4:Then Prog”BIANPO”:IfEnd:Goto 2
2. 正算子程序(SUB1)
0.1739274226 →Z[1]： 0.3260725774→B： 0.0694318442→K： 0.3300094782→L： (1-L) →F： (1-K)→M： U+W(Z[1]cos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Z[1]cos(G+QEMW(C+MWD))) →X： V+W(Z[1]sin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))
+Bsin(G+QEFW(C+FWD))+Z[1]sin(G+QEMW(C+MWD))) →Y： G+QEW(C+WD)+90→F： X+ZcosF→X： Y+ZsinF→Y
3. 反算子程序(SUB2)
G-90 →Z[2]： Abs((Y-V)cosZ[2]-(X-U)sinZ[2]) →W： 0→Z：Lbl 0：Prog “SUB1″： Z[2]+QEW(C+WD) →L： (J-Y)cosL-(I-X)sinL→Z：
If Abs(Z)<1E-6:Then Goto1：Else W+Z→W：Goto 0: IfEnd
Lbl 1： 0→Z：Prog “SUB1″： (J-Y)÷sinF→Z
4. 数据库(SUB3)
If S<226100:Then Cls:Stop: IfEnd
If S<226255.833:Then 226100→O： 99037.9736→U： 4282.3590→V： 201.3266304→G： 1045→P： 1045→R： 253.543→H： 0→Q: Goto 1:IfEnd
……
……
If S>236005.651:Then Cls:Stop: IfEnd
Lbl 1: 1÷P→C：(P-R)÷(2HPR) →D：180÷π→E：
5. 算方位角程序(FA)
LbI0：“XC=”? A : “YC=”? T :X：Y
0→J： POI（X-A，Y-T）→I：“BC=”：I◢If J＜0： Then 360+J→J：ElseJ→J :IfEnd
“FW=”：J►DMS◢ Goto0
6. 使用方位角算坐标(NE)
LbI0：“X=”? X: “Y=”? Y: “L=”? L:“F0”?F： L×Cos（F）+X→N: “N=”N◢L×Sin（F）+Y→E: “E=”E◢
Goto0
7. 边坡放样程序(BIANPO)
Prog ” TCD “：Prog ” BHD “：IF S<226255.731 : then 232.257→ Z[17]:225844.269→ Z[18]:-0.69→Z[19]:11789.745→ Z[22]: Goto 3: IfEnd
IF S<226351.25 : then 230.837→ Z[17]:226050→ Z[18]:2.8→Z[19]:10^45→ Z[22]: Goto 3: IfEnd
IF S<226688.75 : then 239.272→ Z[17]:226351.25→ Z[18]:2.8→Z[19]:-45000→ Z[22]: Goto 3: IfEnd
IF S<227318 : then 243.997→ Z[17]:226520→ Z[18]:2.05→Z[19]: 10^45→ Z[22]: Goto 3: IfEnd
IF S<227982 : then 260.356→ Z[17]:227318→ Z[18]:2.05→Z[19]:-16000→ Z[22]: Goto 3: IfEnd
Lbl 3:Z[17]+(S- Z[18]) ×Z[19]/100+ (S- Z[18])^2/2/ Z[22] – (13+ Z[20] )×0.02→ Z[4]: “QPDGD=”：Z[4] ◢
Z[7]-Z[4] → Z[11]：If Z[11]〈0 Then Abs（Z[11] ）→Z[11]: “T=”：Z[11] ◢13+ Z[20]+ Z[21]→ Z[3]: 8→ Z[5]: 0→ Z[6]: 1.5→ Z[8]: 1.75→ Z[9]:Else“W=”：Z[11] ◢15.6+ Z[20]+ Z[21]→ Z[3]: 10→ Z[5]: 2→ Z[6]: 1.0→ Z[8]: 1.25→ Z[9]:IfEnd
Z[11]÷Z[5] →Z[12]：If Z[12]〈1 Then Z[8]×Z[11] →Z[13]：Else Z[8]×Z[5]+Z[9]×Z[5]×（Z[12]-1）→Z[13]：IfEnd: AbsZ-（Z[13]+Z[6] ×Int Z[12]+Z[3]）→Z[14]:“BP+C-J=”: Z[14] ◢
说明： 只适用二级边坡。
Z[14]=“BP+C-J=” ………………边坡超欠挖宽度（+表超，-表欠）
8. 隧道圆心放样程序（SDYX）
IF S<229000:THEN 243.997+0.0205×(S-226520)-(S-227318)^2/2/16000+1.503 →Z[10]:“NGC=”：Z[10] ◢ELSE 286.385-0.023×(S-230500)- (S-230593.753)^2/2/165000+1.503→Z[10]:“YLJ=”：Z[10] ◢IfEnd
Z[10]- Z[7] →Z[15]:“+SG,-JD”：Z[15] ◢√(（Z[15]^2+(ABS(Z)-6.8)^2）)→Z[16]: “R=”:Z[16] ◢
9. TCD（停车带）
IF S<227020 : then O→ Z[20]:reture: IfEnd
IF S<227035 : then (S-227020) ×3.5/15→ Z[20]: reture: IfEnd
IF S<227065 : then 3.5→ Z[20]: reture: IfEnd
IF S<227080 : then (S-227065) ×3.5/15→ Z[20]: Else 0→Z[20] :IfEnd
10. BHD（变化段）
IF S<227300 : then O→ Z[21]: reture: IfEnd
IF S<227450 : then (S-227300) ×1.85/300→ Z[21]: Else 0.925→Z[21] :IfEnd
三、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
2、输入与显示说明
SZ → XY，2.XY → SZ，3.SDYX，4. BIANPO
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算 里程和边距。输入3表示隧道圆心放样和计算实测半径；输入4表示进行路基边坡放样。
3、盖板涵斜交操作说明：
（以左前点为例）
先把盖板涵斜交角度减９０后＝Ａ
再盖板涵左长度*ＣｏｓＡ＝Ｂ，盖板涵左长度*ＳｉｎＡ＝Ｃ，盖板涵身宽度/ＣｏｓＡ＝Ｄ
利用YQXJS程序计算此桩号正交时Ｂ宽度的坐标值，及法线方位角。
再转到NE程序，Ｘ，Ｙ，Ｆ三项不要修改，只改Ｌ＝Ｃ＋Ｄ/２（计算后点Ｄ前为—）
4、SUB3数据库说明：
S待求桩号，O=曲线起点桩号，U曲线起点X坐标，V曲线起点Y坐标，G曲线起点切线方位角，
P曲线起点半径， R曲线终点半径， H=曲线长度，Q转向（直线0，左转-1，右转1）