由本站论坛　bg6iam　提供，谢谢



Casio fx-4800 直线、曲线、螺旋线、匝道中边桩坐标计算程序

Lb1  0：{DEFG}：A“XA”：B“YA”：C“CA”：D“1÷RA”：E“1÷RB”：F“DKA”：G“DKB”
Lb1  1：｛HOR｝：H“DKI”：O“DL”：R“DR”
H>G=>Goto 0⊿
H<F=>Goto 0⊿
P=（E－D）÷Abs（G－F）：Q= Abs(H－F)：I=P×Q
J=C＋（I＋2×D）×Q×90÷π
M=C＋（I÷4＋2×D）×Q×45÷（2×π）
N=C＋（3÷4×I＋2D）×Q×135÷（2×π）
K=C＋（I÷2＋2×D）×Q×45÷π
X=A＋Q÷12×（cosC＋4×（cosM＋cosN）＋2×cosK＋cosJ）
Y=B＋Q÷12×（sinC＋4×（sinM＋sinN）＋2×sinK＋sinJ）

U“XL”=X＋O×cos（J－90）◢
V“YL”=Y＋O×sin（J－90） ◢
W“XR”=X＋R×cos（J＋90）◢
Z“YR”=Y＋R×sin（J＋90） ◢
Goto 1
Lb1 3：A＝X：B＝Y：D＝E：F＝G：Goto 0

本程序是根据复化辛普生公式（n=4）编制的
计算器显示及输入方法：
XA？、YA？：路线起算点的X、 Y坐标。
CA？路线起算点切线方位角。
1÷RA、1÷RB：曲线元起点及终点半径，以倒数输入，如起点、终点半径分别为60和40米则输入1÷60，1÷40，
注意：遇ZH，HZ，GQ点时半径都输入0。左传为负、右转为正，直线输0。
DKA？、 DKB：路线的起点里程和终点里程。
DKI？：计算点的里程。
DL？、DR？：输入左、右边桩的距离。可以0输入。
DKI= ：显示计算点里程。
X= 、 Y=  显示要求里程点的X、Y坐标。
XL= 、 YL=  ：显示计算点的左边桩X、Y坐标。
XR= 、 YR=  ：显示计算点的右边桩X、Y坐标。

说明：
1、DKI？再次出现输入下一个要求的点的桩号，如果输入数据大于DKB？所输入的数据，则会再次出现RB？这时输入下一段线型的半径，直线输0，DKB？也会随后出现，这次出现输入下一段线型的终点里程桩号。则可以正常运算其它桩号点的坐标。
2、本程序可以反复计算下一线型。
3、上述说明中线型是指曲线上的某一段，例如直缓点到缓圆点则为一段的线型。