5800计算器程序，求高人修改下

线路坐标正反算
1. 主程序(TYQXJS)
Lbl 4："1.SZ => XY"："2.XY => SZ"：？N：？S：Prog“SUB0”↙
1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1： Goto2：↙
Lbl 1:? Z： Abs(S-O) →W：Prog "SUB1"： "XS="：X◢
"YS="：Y◢F-90→F： “FS=”：F▲DMS◢：Goto4↙
Lbl 2：？X： ？Y：X→I：Y→J： Prog“SUB2”：O+W→S：“S=”：S◢“Z=”：Z◢Goto4↙

2. 正算子程序(SUB1)
0.1739274226→A： 0.3260725774→B： 0.0694318442→K： 0.3300094782→L： 1-L→F： 1-K→M：
U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW
(C+FWD))+Acos(G+QEMW(C+MWD))) →X： V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) →Y：
G+QEW(C+WD)+90→F： X+Zcos（F）→X： Y+Zsin（F）→Y
2. 反算子程序(SUB2)
G-90→T： Abs((Y-V)cosT-(X-U)sinT) →W： 0→Z：Lbl 0：Prog "SUB1"：
T+QEW(C+WD) →L： (J-Y)cos（L）-(I-X)sin（L）→Z：Abs（Z）<1E-6=>Goto1： W=W+Z：Goto 0Δ←┘
Lbl 1：Z=0：Prog "SUB1"： (J-Y)÷sin（F）→Z
SUB0 数据库子程序
Goto1↙ 同时保存多个曲线时的指针
Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd↙
Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd
………………………..
为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。
三、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右
侧时，Z取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆
弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45
次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的
值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、输入与显示说明
输入部分：
1. SZ => XY
2. XY = > SZ
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。
S ？正算时所求点的里程，反算时为所求点的近似里程，反算时输入的近似里程不能夸线元，若夸线元的话，需进行第二次反算，并以第一次的结果做为近似里程，否则反算的里程及边距不对， Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
显示部分：
XS=××× 正算时，计算得出的所求点的X坐标
YS=××× 正算时，计算得出的所求点的Y坐标
FS=××× 正算时，所求点对应的中线点的切线方位角
S=××× 反算时，计算得出的所求点的里程
Z=××× 反算时，计算得出的所求点的边距