中边桩坐标计算程序：ZBJS主程序X:Y（这是立仪点坐标）:D“X1”：E“Y1”（前一交点坐标）：F“X2”：G“Y2”（曲线所在交点坐标）：H“X3”：V“Y3”（后一交点坐标）：L“LS”（缓和曲线长）：W“J”（所在曲线交点桩号）I=D-F：J=E-G：Prog“C”（计算前一直线的极坐标）：O（字母O）=J（前一直线方位角）：I=H-F：J=V-G：Prog“C”：U=J（后一直线极坐标及方位角M=O（字母O）-U（计算交点夹角）：M<0（数字O）M=M+360◣M<180M=180-M:C=1:M=M-180:C=-1◣P=L^2/(24R):Q=L/2-L^3/(240R^2):T=Q+(R+P)TAN(M/2)◢S=RMπ/180+L◢M=MC:M“J0”◢M=W-T：N=M：N“ZH”◢N=M+L：N“HY”◢N=M+S/2：N“QZ”◢N=M+S-L：N“YH”◢N=M+S：N“HZ”◢（计算和显示曲线主点桩号，计算坐标时，可以不用显示，将黑三角形改为：就行了）Prog“B”B（子程序、主要是进行桩号判断，进行分段处理）Lbl0：Fixm:{KZB}:I=K-W+T:J=Z:B（输入K：计算的桩号，Z：与中线的夹角、B：边桩的距离，左负右正，这是为计算边桩坐标用的）I<S-LI≦0M=I+BcosZ:N=BsinZ:Goto2◣C=1J=180-Z◣I≦LProg“D”:Prog“F”◣Lbl2:J=O（字母O）:Prog“G”:Goto0◣I=S-I:I≦0M=I-BcosZ:N=-BsinZ:C1J=180-Z◣Prog“D”:N=-N◣J=U:Prog“G”：Goto0C（子程序、主要是进行坐标转换，转换为极坐标）Pol（I，J）：J<0J=J+360◣D（子程序、主要计算缓和曲线上的X和Ｙ）A=RL：M=I-I^5/(40A^2)+I^9/(3456A^4):N=I^3/(6A)-I^7/(336A^3)+I^11/(42240A^5):A=90I^2/(Aπ):Prog“E”E（子程序，主要是计算边桩独立坐标）A=A-J:M=M-CBcosA:N=CN-BsinAF子程序（计算圆曲线上独立坐标）A=90(2I-L)/(Rπ):M=RsinA+Q:N=P+R(1-cosA):Prog“E”G（子程序，计算和显示坐标）A=cosJ:J=sinJ:I=F+AT-MA+NJ:J=G+TJ-MJ-NA:M=I:N=J:I=I-X:J=J-YProg“C”:M”X=”◢N”Y=”◢（计算和显示坐标）I◢（放样时显示的距离）J◢（放样时显示的角度，在计算坐标时，可以不用显示，将黑三角形改为：就行了）说明：1、由于显示的原因，所有“◣”均为空心的，不知怎么回事，全变成实心的了。2、程序源于网络，只进行了部分修改和完善。3、所有（）中的均为说明，程序中不要出现。为便于看懂，增加一个操作实例：1、先运行程序，会出现XY，X=1000，Y=1000（这是立仪点坐标，可以任意输，没有影响）。X1=3046919.844，Y1=450373.316，X2=3046843.992，Y2=450328.693，X3=3046776.665，Y3=450088.808，这是三个点，确定了两条直线和一夹角。LS=40（缓和曲线长），J=48668.004（这是计算的平曲线的交点桩号，即为48K+668.004），R=142（平曲线半径）。这是基本的资料，都要输的。2、计算结果：T=77.338（切线长）。S=148.688，这是平曲线长，J0=43.85461096（这是转角值，正为右，负为左，换成度分秒就是43度51分16.6秒，先按黄色SHIF，再按。，，，就出来了）。ZH=48590.666（直缓点桩号，下面4个就不说了，你应该知道），HY=48630.666，QZ=48665.010，YH=48699.354，HZ=48739.354。这是计算平曲线的基本数据，在程序中可以不要显示，只要将相应的黑色三角换成回车就行了。3、K？输入要计算的点的桩号，如48600（K48+600）。Z？是与中线的夹角，正交就输90。一般都输90就行了。B？输入左右边桩的距离，左负右正，中桩输0就行。得出：X=3046902.618，Y=450363.154，这就是K48+600中桩的坐标。完成计算。4、如果还要计算，直接按EXE键，再次要求输入桩号，如果要算桩，直接按就行了，刚已经输过48600了，Z？也直接按，B？输入8，表示右边桩8米的坐标，X=3046906.727，Y=450356.2903。5、因为是计算器，所以只能一个个的输，一个个的算，边