4800程序请高手修改，适合水利工程的放样！

圆曲线和直线正反算,非常适合于水利工程中渠道的测量放样,但不知哪个语句有问题,圆曲线右偏时,计算结果错误,左偏和直线段都是对的,请大师百忙之中能抽出点空修改一下,不胜感谢!!!
CASIO fx-4800P程序编制（里程中边桩计算坐标及反算）
程序编制
一、源程序(SZ-XY-SZ) "1.SZ => XY"："2. XY => SZ"：NU"X0"V"Y0"O"S0"G"F0"Q：Q=0=>R=1E-45：≠> R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙ Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：F=G+B：E=F+B+90：C=LK+2MRsinQB：X"XS" = U+Rec(C,F)+ZcosE◢ Y"YS"=V+J+ZsinE◢ Goto 1←┙ Lbl 2：W=G+90Q：E=U+Rec(R,W)：F=V+J←┙ Lbl 3：{XY}：XY：D=Pol(X-E,Y-F):T=J：J<0=>T=J+360△B=T-W-180M：P=Abs B：S"S"=O+ DKcosP+ MπRP÷180◢ Z"Z"=DKsinB+MQ(R-D)◢ Goto 3←┙ 二、使用说明 1.由里程、中边桩计算坐标 启动程序后按以下提示输入数据： N ？1 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩) X0 ？142.353 (直线或圆曲线起点的X坐标） Y0 ？368.166 (直线或圆曲线起点的Y坐标） S0 ？155.252 (直线或圆曲线起点的里程） F0 ？60°54′03″(直线或圆曲线起点的切线方位角） Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且 右偏时Q=1) R ？ 100 (圆曲线的半径) S ？ 197.016 (中线里程) Z ？ -5.251 (中线左侧5.251m的边桩；在中线上Z=0，在左侧<0，在右侧Z>0) 显示计算结果 XS=172.7467767 YS=395.2316807 2.由坐标反算里程、中边桩 启动程序后按以下提示输入数据： N ？2 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩) X0 ？142.353 (直线或圆曲线起点的X坐标） Y0 ？368.166 (直线或圆曲线起点的Y坐标） S0 ？155.252 (直线或圆曲线起点的里程） F0 ？60°54′03″(直线或圆曲线起点的切线方位角） Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且 右偏时Q=1) R ？ 100 (圆曲线的半径) X ？ 172.7467767 (待求里程及中边桩的点的X坐标) Y ？ 395.2316807 (待求里程及中边桩的点的Y坐标) 显示计算结果 S=197.016 Z=-5.251 三、注意事项 1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等 曲线的计算； 2.程序中“R=1E-45”的意思为10的-45次方，输入程序时切莫输错。
在圆曲线右偏时反算计算错误,请大师能修改一下程序?不胜感谢!!!
我计算了一个数据:
上切点坐标:x0=408.6357,y0=689.9914,s0=21406.73,F0=194度14分21.96秒,Q=1,R=250.时下切点(x=220.2193,y=441.6333)的桩号和轴距,正确答案应该是:s=21743.32,z=0.01700961962.
右偏时此程序桩号s计算错误,轴距z正确,请大师详查后修改