10 CLS:PRINT "Resection Program:":CLS:SET N
20 PRINT "Input the co-ords of three knownpoints in a CLOCKWISE direction":CLS
30 PRINT "x1:";x1;:INPUT x1:PRINT "y1:";y1;:INPUT y1:CLS
40 PRINT "x2:";x2;:INPUT x2:PRINT "y2:";y2;:INPUT y2:CLS
50 PRINT "x3:";x3;:INPUT x3:PRINT "y3:";y3;:INPUT y3:CLS
60 PRINT "Azimuth 1:";DMS$(AZ01);:INPUT AZ01:CLS
70 PRINT "Azimuth 2:";DMS$(AZ02);:INPUT AZ02:CLS
80 PRINT "Azimuth 3:";DMS$(AZ03);:INPUT AZ03:CLS:PRINT "Wait...";
90 Z=POL((x2-x1),(y2-y1)):AZ12=90-Y:Z=POL((x3-x2),(y3-y2)):AZ23=90-Y:Z=POL((x1-x3),(y1-y3)):AZ31=90-Y
100 A=AZ31+180-AZ12:B=AZ12+180-AZ23:C=AZ23+180-AZ31
110 c=AZ02-AZ01+360:b=AZ01-AZ03+360:a=AZ03-AZ02+360
120 K1=1/((1/TAN(A))-(1/TAN(a)))
130 K2=1/((1/TAN(B))-(1/TAN(b)))
140 K3=1/((1/TAN(C))-(1/TAN(c)))
150 EP=((K1*x1)+(K2*x2)+(K3*x3))/(K1+K2+K3)
160 NP=((K1*y1)+(K2*y2)+(K3*y3))/(K1+K2+K3)
170 Z=POL((x1-EP),(y1-NP):brg1=90-Y:IF brg1<0 THEN brg1=brg1+360
180 Z=POL((x2-EP),(y2-NP):brg2=90-Y:IF brg2<0 THEN brg2=brg2+360
190 Z=POL((x3-EP),(y3-NP):brg3=90-Y:IF brg3<0 THEN brg3=brg3+360
200 CLS:BEEP 0:PRINT "E: ";ROUND(EP,-5);TAB(32);"N: ";ROUND(NP,-5):CLS
210 BEEP 0:PRINT "Set Bearings";TAB(16);"1:";DMS$(brg1);TAB(32);"2:";DMS$(brg2);TAB(48);"3:";DMS$(brg3)
220 CLS:GOTO 20