130 LETL=0
150 PRINT
160 PRINT"   S P H E R I C A L  T R I A N G L E  S O L U T I O N"
170 PRINT
180 PRINT
185 GOTO2035
200 LETP=3.14159265#
210 C=180/P
250 LETL=L+1
260 PRINT"CASE NUMBER";L
270 PRINT
280 PRINT
290 PRINT"POSITION 1:"
300 PRINT
310 PRINTABS(D0);"DEG  ";M0;"MIN  ";
320 IFD0<0THEN350
330 PRINT"NORTH ";
340 GOTO360
350 PRINT"SOUTH ";
360 PRINT"LATITUDE"
370 PRINTABS(T0);"DEG  ";N0;"MIN  ";
380 IFT0<0THEN410
390 PRINT"WEST ";
400 GOTO420
410 PRINT"EAST ";
420 PRINT"LONGITUDE"
430 PRINT
440 PRINT
450 PRINT"POSITION 2:"
460 PRINT
470 PRINTABS(D1);"DEG  ";M1;"MIN  ";
480 IFD1<0THEN510
490 PRINT"NORTH ";
500 GOTO520
510 PRINT"SOUTH ";
520 PRINT"LATITUDE"
530 PRINTABS(T1);"DEG  ";N1;"MIN  ";
540 IFT1<0THEN570
550 PRINT"WEST ";
560 GOTO580
570 PRINT"EAST ";
580 PRINT"LONGITUDE"
590 PRINT
600 PRINT
610 LETM0=M0/60
620 LETN0=N0/60
630 LETM1=M1/60
640 LETN1=N1/60
650 IFD0>=0THEN680
660 LETD0=(ABS(D0)+M0+90)/C
670 GOTO690
680 LETD0=(90-(D0+M0))/C
690 IFD1>=0THEN720
700 LETD1=(ABS(D1)+M1+90)/C
710 GOTO730
720 LETD1=(90-(D1+M1))/C
730 IFT0>=0THEN760
740 LETT0=-(ABS(T0)+N0)/C
750 GOTO770
760 LETT0=(T0+N0)/C
770 IFT1>=0THEN800
780 LETT1=-(ABS(T1)+N1)/C
790 GOTO810
800 LETT1=(T1+N1)/C
810 LETT=ABS(T1-T0)
820 LETF=FNC(D0)*FNC(D1)+FNS(D0)*FNS(D1)*FNC(T)
830 LETF1=SQR(ABS(1-F*F))
831 LETK4=ABS(F1/F)
832 IFK4<-1THEN838
833 IFK4>1THEN836
834 LETZ=FNA(K4)
835 GOTO850
836 LETZ=FNT(K4)
837 GOTO850
838 LETZ=FNN(K4)
850 IFF>=0THEN870
860 LETZ=P-Z
870 LETG=(FNC(D1)-F*FNC(D0))/(F1*FNS(D0))
871 LETK4=ABS(SQR(ABS(1-G*G))/G)
872 IFK4<-1THEN878
873 IFK4>1THEN876
874 LETB0=FNA(K4)
875 GOTO890
876 LETB0=FNT(K4)
877 GOTO890
878 LETB0=FNN(K4)
890 IFG>=0THEN910
900 LETB0=P-B0
910 IF(T0-T1)>=0THEN930
920 LETB0=2*P-B0
930 LETH=(FNC(D0)-F*FNC(D1))/(F1*FNS(D1))
931 LETK4=ABS(SQR(ABS(1-H*H))/H)
932 IFK4<-1THEN938
933 IFK4>1THEN936
934 LETB1=FNA(K4)
935 GOTO950
936 LETB1=FNT(K4)
937 GOTO950
938 LETB1=FNN(K4)
950 IFH>=0THEN970
960 LETB1=P-B1
970 IF(T1-T0)>=0THEN990
980 LETB1=2*P-B1
990 PRINT"LOCAL HOUR ANGLE (AT NORTH POLE):"
1000 PRINT
1010 PRINTINT(10*T*C+.5)/10;"DEG"
1020 PRINTINT(T*C);"DEG  ";
1030 PRINTINT(600*(T*C-INT(T*C))+.5)/10;"MIN"
1040 LETH7=T*C/15
1050 LETM7=(H7-INT(H7))*60
1060 LETS7=(M7-INT(M7))*60
1070 PRINTINT(H7);"HRS  ";INT(M7);"MIN  ";INT(10*S7+.5)/10;"SEC"
1080 PRINT
1090 PRINT
1100 PRINT"ZENITH (GREAT CIRCLE) DISTANCES:"
1110 PRINT
1120 PRINTINT(10*Z*C+.5)/10;"DEG"
1130 PRINTINT(Z*C);"DEG  ";
1140 PRINTINT(600*(Z*C-INT(Z*C))+.5)/10;"MIN"
1150 PRINTINT(600*Z*C+.5)/10;"NAUTICAL MILES"
1160 PRINTINT(600*Z*C*6080.2/5280+.5)/10;"STATUTE MILES"
1170 PRINT
1180 PRINT
1190 PRINT"TRUE BEARINGS (GREAT CIRCLE COURSES):"
1200 PRINT
1210 PRINT" POSITION 2 FROM POSITION 1:"
1220 PRINTINT(10*B0*C+.5)/10;"DEG"
1230 PRINTINT(B0*C);"DEG  ";INT(600*(B0*C-INT(B0*C))+.5)/10;"MIN"
1240 PRINT
1250 PRINT" POSITION 1 FROM POSITION 2:"
1260 PRINTINT(10*B1*C+.5)/10;"DEG"
1270 PRINTINT(B1*C);"DEG  ";INT(600*(B1*C-INT(B1*C))+.5)/10;"MIN"
1280 PRINT
1290 PRINT
1300 PRINT"ALTITUDE (REMOTE CELESTIAL POSITION ABOVE LOCAL HORIZON):"
1310 PRINT
1330 LETA7=90-Z*C
1340 LETA8=ABS(A7)
1350 IFA7<0THEN1380
1360 PRINTINT(10*A7+.5)/10;
1370 GOTO1390
1380 PRINT-1*INT(10*A8+.5)/10;
1390 PRINT"DEG"
1400 IFA7<0THEN1440
1410 PRINTINT(A7);"DEG  ";
1420 PRINTINT(600*(A7-INT(A7))+.5)/10;"MIN"
1430 GOTO1460
1440 PRINT-1*INT(A8);"DEG  ";
1450 PRINTINT(600*(A8-INT(A8))+.5)/10;"MIN"
1460 PRINT
1470 PRINT
1480 IF(1+H0)*(1+J0)=1THEN1730
1490 PRINT"OBSERVED ALTITUDE:"
1500 PRINT
1510 PRINTH0;"DEG  ";J0;"MIN"
1520 LETH1=H0+J0/60
1530 PRINTINT(100*H1+.5)/100;"DEG"
1540 PRINT
1550 PRINT
1560 PRINT"LINE OF POSITION:"
1570 PRINT
1580 LETI4=60*(H1-A7)
1590 PRINTINT(10*ABS(I4)+.5)/10;"MILES ";
1600 IFI4>0THEN1630
1610 PRINT"AWAY ";
1620 GOTO1640
1630 PRINT"TOWARDS ";
1640 PRINT"ON LINE BEARING  ";
1650 IFI4>0THEN1720
1660 LETB9=B0*C+180
1670 IFB9>360THEN1700
1680 PRINTINT(10*B9+.5)/10;"DEGREES TRUE ";
1690 GOTO1730
1700 LETB9=B9-360
1710 GOTO1680
1720 PRINTINT(10*B0*C+.5)/10;"DEGREES TRUE"
1730 PRINT:PRINT
1790 PRINT"ANOTHER PROBLEM  ";
1800 INPUT A$
1806 PRINT
1810 IF LEFT$(A$,1)="Y" THEN 2130
1820 GOTO 3000
2030 REM THESE SERIES REPLACE TRIG FUNCTIONS. GET ORIG. AND
2031 REM PUT THE FUNCTIONS BACK --- GOOD GRIEF !!!
2035 DEF FNV(X)=X^5/120
2040 DEF FNW(X)=X^6/720
2050 DEF FNX(X)=X^7/5040
2060 DEF FNY(X)=X^8/40320!
2070 DEF FNZ(X)=X^9/362880!
2080 DEF FNS(X)=X-X^3/6+FNV(X)-FNX(X)+FNZ(X)
2090 DEF FNC(X)=1-X^2/2+X^4/24-FNW(X)+FNY(X)
2100 DEF FNA(X)=X-X^3/3+X^5/5-X^7/7
2110 DEF FNT(X)=3.14159/2-1/X+1/(3*X^2)-1/(5*X^5)
2120 DEF FNN(X)=-3.14159/2-1/X+1/(3*X^2)-1/(5*X^5)
2130 PRINT
2150 PRINT"LATITUDE OF POSITION 1 (DEG,MIN) ";
2160 INPUTD0,M0
2170 PRINT"LONGITUDE OF POSITION 1 (DEG,MIN) ";
2180 INPUTT0,N0
2190 PRINT"LATITUDE OF POSITION 2 (DEG,MIN) ";
2200 INPUTD1,M1
2210 PRINT"LONGITUDE OF POSITION 2 (DEG,MIN) ";
2220 INPUTT1,N1
2230 LETH0=0
2240 LETJ0=0
2250 PRINT
2260 PRINT
2270 GOTO200
3000 END
ITUDE OF POSITION 2 (DEG,MIN) ";
2220 INPUTT1,N1
2230 LETH0=0
2240 LETJ0=0
2250 PRINT
2260 PRINT
2270 GOTO200
3000 EN