S=SU ←┘ Goto 0 ←┘ Lbi 2 ←┘ Z=H-A ←┘
O=90Z2÷(∏RV) ←┘
X=Z-Z^5÷(40R2V2)+Z^9÷(3456R^4V^4) ←┘
Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ←┘ Lbi 6 ←┘ Rec(X,F) ←┘
X=Z[1]+I:Y=Z[2]+J ←┘ Rec(Z,F+90G) ←┘ X=X+I:Y=Y+J ←┘ Rec(S,F+OG+EU) ←┘ \ \ S=SU ←┘ Goto 0 ←┘ Lbi 3 ←┘ Z=H-A-V÷2 ←┘ O=180Z÷(R∏) ←┘ X=Rsin O+M ←┘ Z=R(1-cos O)+P ←┘ Goto 6 ←┘
Lbi 4 ←┘ Z=D-H ←┘
O=90Z2÷(VR∏) ←┘
X=Z-Z^5÷(40R2V2)+Z^9÷(3456R^4V^4) ←┘
Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ←┘ Rec(X,F+GN+180) ←┘ X=Z[3]+I:Y=Z[4]+J ←┘ Rec(Z,F+GN+180-90G) ←┘ X=X+I:Y=Y+J ←┘
Rec(S,F+GN+180-OG-(180-E)U) ←┘ \ \ S=SU ←┘ Goto 0 ←┘ Lbi 5 ←┘
Rec(H-D+T,F+GN) ←┘ X=W+I:Y=K+J ←┘ Rec(S,F+GN+EU) ←┘ \ \ S=SU ←┘ Goto 0 ←┘
请注意 :∏ 这是圆周率 3.14 的附号 ,10E 这个 E 请按 \键输入 .
不明白者联系 QQ326275120 注： XY 程序单代号
Q ——交点里程桩号（去 K 后以纯数字表示输入） W ——交点 X 坐标值 K ——交点 Y 坐标值 R ——圆曲线半径 F ——第一直线段方位角 N ——交点偏角 V ——缓和曲线长度
G ——如果为左偏角输入“ -1 ” ，如果为右偏角输入“ 1 ” S ——边桩至中桩的距离，求中桩时输入 0 H ——所求点的里程桩号
E ——左边桩输入“ -90 ” ，右边桩输入“ 90 ” XY ——所求点的坐标值 二、水平测量 程序名 :SPCL Lbi 0 ←┘ M:A ←┘ {B} ←┘ \B ◢ Goto 0 ←┘ M ——水准点高程 A ——后视水准点计数
B ——前视所求点计数 H ——所求点高程
卡西欧5800P计算器曲线任意里程中边桩坐标正反算
一、程序功能
本程序由一个主程序(ZHUCHENXU)和几个子程序——正算子程序(SUB1)、反算子程序(SUB2)、数据库(SUB3)、算方位角程序(FA)、使用方位角算坐标(NE)、边坡放样程序(BIANPO)、隧道圆心放样程序（SDYX）构成，可以根据直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。
二、源程序
增加变量16→DimZ 1.主程序(ZHUCHENXU)
\3.SDYX，4. BIANPO \？N：
Lbl 1：”ZHUANGHAO=”?S： Prog \： N≥2=>Goto 2： Abs(S-O)→W：”BIANZHU=”?→Z:Prog \
\-90) →F:F►DMS◢Prog \Goto 1: Lbl 2：”CX=”?X: X→I：”CY=”?Y: Y→J：”DMG=”?→Z[7]：Prog \： \→S◢\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 \： 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 \： (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)
“QPDJL=”? → Z[3]: “QPDGD=”? → Z[4]: “PG=” ? → Z[5]：2→Z[6]： “P1=” ? → Z[8]： “P2=” ? → Z[9]： “P3=” ? → Z[10] (没有的话就删)： Z[7]-Z[4] → Z[11]：If Z[11]〈0 Then AbsZ[11] →Z[11]: “T=”：Z[11] ◢Else“W=”：Z[11] ◢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[3]= QPDJL ?????? 起坡点离中桩距离 Z[4]=“QPDGD” ?????起坡点设计标高 Z[5]=“PG” ?????? 坡高
Z[6]=“PTK” ?????? 平台宽度 Z[7]=“DMG” ??????实测地面高度 Z[8]=“P1” ?????? 一级边坡坡度 Z[9]=“P2” ?????? 二级边坡坡度
Z[10]=“P3” ?????? 三级边坡坡度(没有的话就删)
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] ◢ 三、使用说明 1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆