缓和曲线与原曲线任意点坐标计算程序
Z=90S÷(R∏):U=S2÷(24R)-S^3÷(2688R^3):Q=S÷2-S^3÷(240R2):T=(R＋U)tan(N÷2)＋Q▲
L=R(N-2Z)∏÷180▲
Lbl0:{G}:K=Abs(G-F):S=0=＞Goto2: ≠＞Goto1⊿
Lbl1:K≤S=＞C=K-K^5÷(40R2S2) +K^9÷(3456R^4S^4):V=K^3÷(6RS)-K^7÷(336R^3S^3)+K^11÷(42240R^5S^5):O=tan-1(V÷C): Goto3: ≠＞Goto2⊿
Lbl2: O=180(K-S)÷(R∏)+Z:C=RsinO+Q:V=R(1-cosO)+U
Lbl3:Rec(C,B):X=D+I:Y=P+J: Rec(V,B):X=X+I:Y=Y+J▲
{H,W}:A=B+MO:ReC(H,A+90W):X=X+I▲Y=Y+J▲
Goto0
S：缓和段长度
R：圆曲线
∏：圆周率
T：缓和曲线切线长
N：缓和曲线偏角（输入绝对值）
L：圆曲线长
K：距离缓直点或圆起点距离
G：计算里程桩号F：缓直点（直缓点）或直圆点（圆直点）里程桩号
D：缓直点或圆起点X坐标P：缓直点（直缓点）或圆起点（圆终点）Y坐标
B：前切线方位角（计算缓直至圆缓或直圆至圆直用）或后切线方位角（计算直缓至缓圆或圆直至直圆用）
M：缓和曲线偏角方向控制（左输“ -1” ，右输“ 1” ）（反向计算时相应反号）
X：计算点X值
Y：计算点Y值
H：边桩至中桩距离W：计算左边桩坐标时，输“ -1” 。
适用于计算：1、从缓直点至圆缓点或直缓点至缓圆点坐标
2、从直圆点至圆直点或圆直点至直圆点坐标，此时S=0。
本人编的4800计算器程序，有什么不到之处，恳请各位朋友给出宝贵的意见，本人一定予以改进。在此多谢各位！
本文来自：中国路桥网（www.9to.com），原文地址：http://www.9to.com/Article/Technique/daolu/200902/29227.html