5800全线任意坐标计算程序
1. 正算主程序(ZHCX) （不运行）
8→DimZ
1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E
“Z=”?Z:”YJJ=”?A:Abs(S-O)→W
0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ]
U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X:
V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y：
G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd
”X=”:X→X◢
”Y=”:Y→Y◢
If F﹤0:Then F+360→F:IfEnd
”QX FWJ=”:F▼DMS◢
“C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标
Lbi 1
X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3
Lbi 2
X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3
Lbi 3
“QH-X=”: X →X◢
“QH-Y=”: Y →Y◢
Prog “FY”
2 . 参数子程序（直接运行）
M(主线) 一条线路一个名称
“S=”?S
If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd
…
…
If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd
Lbi 1
Prog “ZBJS”
3. 放样程序（FY）（不运行）
“X0=”?M:“Y0=”？N
Pol((X-M, Y-N)
If J﹤0:Then J+360→J:IfEnd
“FWJ=”:J▲DMS◢
“JU-LI=”:I◢
说明：
一、程序功能及原理
1.功能说明：
本程序由一个主程序和两个子程——正算主程序(ZH)、参数子程序（M）、放样程序（FY），可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径、左偏或右偏）及里程边距或坐标，对该线元段范围内任意里程中边桩坐标进行正算。

本程序可以在CASIO fx-4800P计算器及CASIO fx-4850P计算器上运行。

2．计算原理：
利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。

利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。

二、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1。

(2) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(3) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。

(4) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

(5) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、运行参数子程序（要做哪条线就直接运行哪条线的参数子程序）
输入
S （S=?）为任意里程桩号
Z (Z=?) 为斜距（当所求点位于中线时，Z=0；当位于中线左侧时，Z 取负值；当位于中线右侧时，Z取正值）
A (YJJ=?) 为右夹角（0≦W≦180）
输出
X=********* 计算的X坐标值
Y=********* 计算的Y坐标值
QXFWJ= *°*′*″计算的里程桩号K的切线方位角
输入
C=? 当C=1时，为计算斜交斜做的桥涵坐标；当C=2时，为计算斜交正做的桥涵坐标
QHJU=? 相对于刚计算完坐标的那个点的大桩号（+）或小桩号（-）的距离输出
QH-X=*******.*** 计算的斜交斜做或斜交正做的X坐标值
QH-Y=********* 计算的斜交斜做或斜交正做的Y坐标值输入
X0=? 测站的X值
Y0=? 测站的Y值
输出
FWJ= *°*′*″计算的X 、Y对于测站的方位角
JU-LI=******* 计算的X 、Y对于测站的距离。