急需一条公路4800圆曲线计算程序

急需一条公路4800圆曲线计算程序
我是个初学者
刚在实习中
对这些还不是很明白
希望那位大哥帮帮我
给条程序让我做参考?

有缓和曲线或无缓和曲线任意桩号中、边桩坐标4800计算程序
有缓和曲线或无缓和曲线任意桩号中、边桩坐标4800计算程序
四、QXZB(曲线坐标)
A＂ZH(ZH)=″:B＂X(ZH)=″:C＂Y(ZH)=″ ＂X(HZ)=″:E＂Y(HZ)=″:F＂FWJ-1=″:G＂FWJ-2=″:H＂R=″:I＂L1=″:J＂L=″:K＂L2=″:L＂ZJ(Z=1,Y=-1)=″
O=180/pi (pi为圆周率,电脑中找不到此符号以此代替)
LbI A
{MN}
M＂ZHUANG HAO=″
P=M-A
MM>(A+I+J+K)?GOTO A⊿
N＂ZHUANG JU(Z+,Y-)=″
M>(A+I+J)?GOTO D⊿
M>(A+I)?GOTO C⊿
LbI B(第一缓和曲线)
X=B+(P-(P5÷40÷H2÷I2))COSF+(P3÷6÷H÷I-P7÷336÷H3÷I3)LSINF
Y=C+(P-(P5÷40÷H2÷I2))SINF-(P3÷6÷H÷I-P7÷336÷H3÷I3)LCOSF
X=X+NCOS(F-(P2÷2÷H÷I)OL-90)◢
Y=Y+NSIN(F-(P2÷2÷H÷I)OL-90)◢
GOTO A
LbI C(圆曲线)
X=B+(HSIN(((P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2))COSF-(H(1-COS(((P-I)÷H+I÷2÷H)O))+ I2÷24÷H)(-1)LSINF
Y=C+(HSIN(((P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2))SINF+(H(1-COS(((P-I)÷H+I÷2÷H)O))+I2÷24÷H)(-1)LCOSF
X=X+NCOS(F-((P-I)÷H+I÷2÷H)OL-90)◢
Y=Y+NSIN(F-((P-I)÷H+I÷2÷H)OL-90)◢
GOTO A
LbI D(第二缓和曲线)
Q=A+I+J+K-M
X=D-(Q-Q5÷40÷H2÷K2)COSG+(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LSING
Y=E-(Q-Q5÷40÷H2÷K2)SING-(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LCOSG
X=X+NCOS(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90)◢
Y=Y+NSIN(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90)◢
GOTO A
变量说明：
以上变量依次为：
ZH(ZH)=直缓点桩号,无缓和曲线则为直圆点；
X(ZH)=直缓点X坐标；无缓和曲线则为直圆点；
Y(ZH)=直缓点Y坐标,无缓和曲线则为直圆点；
X(HZ)=缓直点X坐标,无缓和曲线则为圆直点；
Y(HZ)=缓直点Y坐标,无缓和曲线则为圆直点；
FWJ-1=第一切线方位角；
FWJ-2=第二切线方位角；
R=半径；
L1=第一缓和曲线长；
L=圆曲线长；
L2=第二缓和曲线长；
ZJ(Z=1,Y=-1)=转角形式(左=1,右=-1)
ZHUANG HAO=所求点桩号
ZHUANG JU(Z+,Y-)=边桩桩距,左边桩为正,右为负,中桩0
我已经用过了,没发现什么问题,请大家提意见.