一、程式：
Defm 1：N=12：
Lbl 1：S“L1=”：A“X1=”：B“Y1=”：C“A1=”：O=C：I“R1=”：D=I：K“L2=”：T“R2=”：F=T：I=0＝＞T=0＝＞J=1       J“L-R=”
K－S＞0＝＞Z=1：≠＞Z=－1   D≠0＝＞D=D-1    F≠0＝＞F=F-1       
Lbl 2                                                   ‘将半径转换为曲率
Z[1]=0：Q=0：U=0:V=0：P=N+1：{L}：L“LC=”
R=1＝＞Goto3
L＞S＝＞Goto1      ‘当反算时，若计算点在线元里程范围外，则转至Lbl 1重新输入线元要素
L＜K＝＞Goto1
Goto4
Lbl 3               ‘当正算时，若计算点在线元里程范围外，则转至Lbl 1重新输入线元要素
L＜S＝＞Goto1
L＞K＝＞Goto1
Lbl 4
H=（L－S）÷N：D≠0＝＞Goto5   F= 0＝＞Goto 8
Lbl 5                                判断线元是否为直线，是则转至直线计算程式Lab 8
Dsz  P：Goto 6：Goto 7
Lbl 6                       ‘计算各点曲率及方位角
E=P÷2：G=D+（F－D）÷（K－S）×H×P：M=C+（G+D）×H×P×90÷π×J×Z
P=N＝＞O=M   E≠Int E＝＞Z[1]=Z[1]+CosM：Q=Q+SinM：≠＞U=U+CosM：V=V+SinM
Goto 5     
X=A+Abs H÷3×（Cos C+4×Z[1]+2×U－CosO））▲
Y= B+Abs H÷3×（Sin C+4×Q+2×V－Sin O）▲
Z=－1＝＞Z=Z+180
O＞360＝＞O=O－360：≠＞O＜0＝＞O=O+360
O”A=”▲
Prog “LRZB”
Goto 2
Lbl 8                     ‘直线上坐标计算程式
H=（L—S）×Z
X=A+H×Cos O▲          ‘线路中心点坐标
Y=B+H×Sin O▲          ‘线路中心点坐标
O”A=”▲                  ‘线路中心点切线方位角
Prog “LRZB”
Goto 2
二、变量说明：
N—-曲线元N值   S—曲线元起点里程   A—-起点X坐标  B—-起点Y坐标  C—-起点切线方位角
I、D—-起点半径   K—曲线元终点里程    T、F—-终点半径   J—-线元左右偏判别（1右-1左）
R—-曲线正反算判别（1为正算-1为反算）   注：正算即从起点算至终点，反算即从终点算至起点
L—-计算点里程  Z—-坐标反算判别
子程式“LRZB”
‘计算线路中线左右两侧点坐标
Lbl1：W=-1：{WR}：W“ANG=”：W=1＝＞Goto 2   R“D=”：O=O+W↓
X“[X]”=X+R×Cos O▲       ‘线路侧点坐标
Y“[Y]”=Y+R×Sin O▲        ‘线路侧点坐标
O＞360＝＞O=O－360：≠＞O＜0＝＞O=O+360
O“A=”▲     ‘侧点计算基点至该侧点连线方方位角
Goto 1
Lbl 2
说明：
本子程式计算曲线两侧点坐标,可以无限计算连续点坐标,前提是当提示”ANG”时,输入转向角度就是了.如果输入”-1”则回到计算中线坐标上来.重新计算下一点坐标
W—-夹角（相当于曲线里面的转向角，为前一直线（或切线）的延长线至计算点的夹角
R—-前一点至计算点的直线长度
前方交会是分别在已知坐标点A、B安置经纬仪向待定点P观测水平角 和检查角 ，以确定待定点P的坐标。为保证交会定点的精度，在选定P点时，应使交会角 位于30°~120°之间，最好近于90
变量对照表
数学模型变量 fx-4800P变量 
 P °′″ A点的水平角观测值
  Q °′″ B点的水平角观测值
  A  
  B  
  C °′″ 检查角
  E,F m A点坐标
  K,L m B点坐标
  M,N m C点坐标
  X,Y m P点坐标
  W ″ 检查角差值
程序名：QFJH
E: F: K: L
{PQ}
A=1÷tan P: B=1÷tan Q
X=(EB+KA+L-F)÷(A+B) 
Y=(FB+LA+E-K)÷(A+B) 
G=X-E: H=Y-F:Pol(G,H)
J 0 Z=J+360:  Z=J 
M: N
G=M-E: H=N-F:Pol(G,H)
J 0 J=J+360 
{C}
W=Z-J:W 0 W=W+360 
W=3600(W-C) 
案例：
前方交会数据
点名  坐标(m)  坐标(m) 、 水平角
A 781.227 1323.385   53°33′42″
B 769.862 1428.989   50°34′55″
C 866.294 1306.117   54°03′08″
程序步骤：
步骤 显示 按键 注释
1 E? ××× 781.227  输入A点坐标
2 F? ××× 1323.385  
3 K? ××× 769.862  输入B点坐标
4 L？××× 1428.989  
5 P? ××× 53 33 42   输入 
6 Q？××× 50 34 55   输入 
7 X=843.5338463   P点坐标
8 Y=1380.640577   
9 J=42.58080857   中间计算结果
10 M? ××× 866.294  输入C点坐标
11 N？××× 1306.117  
12 C? ××× 54 3 08   输入检查角 
13 W=11.91808111   检查角差值