CASIO fx-4800P計算器程序
DA DI——N E (公路彎道坐標計算)
K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”
M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（240R2）▲
N“HY”=M＋S▲
O“QZ”=M＋（πRAbs F÷180＋S）÷2▲
Z[1]“YH”=M＋πRAbs F÷180▲
Q“HZ”= Z[1]＋S▲
T=K－M▲
L=Q－M▲
E=（R＋S2÷（24R））÷cos（F÷2）－R▲
Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△
{U}：U“I——B OUT＋IN－”
Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋ S2÷（24R）－Ucos H：△
F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsin Abs F－Tsin Abs F－Wcos F：△△
C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△
I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲
J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲
G=G＋20
Goto 0
程序運行
輸入：JD？ ——交點樁號
R？ ——圓曲線半徑
LS？ ——緩和曲線長度
PJ——L－R＋？ ——偏角 右偏為正
ZH——N？ ——ZH點N坐標
ZH——E？ ——ZH點E坐標
JD——N？ ——JD點N坐標
JD——E？ ——JD點E坐標
輸出：ZH= ——ZH點樁號
HY= ——HY點樁號
QZ= ——QZ點樁號
YH= ——YH點樁號
HZ= ——HZ點樁號
T= ——切線長
L= ——曲線長
E= ——外矢距
再輸入：ZHUANG HAO？ ——輸入任一點樁號
I——B OUT＋IN－？ ——距中線距離（外側為正）
再輸出：COM——N= ——計算出N坐標
COM——E= ——計算出E坐標
以此循環——再輸入——再輸出
程序中各字母所代表的含義：
A——ZH點N坐標
B——ZH點E坐標
C——JD點N坐標
D——JD點E坐標
E——外矢距
F——偏角（右偏為正）
G——彎道內任一點樁號
H——彎道內任一點切線角
I——計算出的N坐標
J——計算出的E坐標
K——交點樁號
L——曲線長
M——ZH點樁號
N——HY點樁號
O——QZ點樁號
P——ZH--JD方位角
Q——HZ點樁號
R——圓曲線半徑
S——緩和曲線長度
T——切線長
U——距中線距離（外側為正）
V——支距
W——支距
X——彎道內坐標
Y——彎道內坐標
Z——任一點曲線長度
Z[1]——YH點樁號
作者：路更群
