File GRKRSS.FC (FOCAL source file)

Directory of image this file is from
This file as a plain text file

C-PS/8 FOCAL, 1971
 1.01 C DIESES PROGRAMM === GRKRSS ===  BERECHNET DIE ENTFERNUNG
 1.02 C ZWEIER PUNKTE AUF DEM ERDELLIPSOID. DAMIT MAN GEKNICKTE
 1.04 C SOWIE DREIECKSKURSE NICHT EINZELN ZUSAMMENSETZEN MUSS, 
 1.05 C SIEHT DAS PROGRAMM AUCH DIESES VOR.
 1.09 C BEI DEM CDC-PROGRAMM WIRD MIT MITTLEREM ERDRADIUS
 1.10 C UND DER ANNAHME, DIE ERDE SEI EINE KUGEL, GERECHNET. DAHER ER-
 1.11 C GEBEN SICH UNTERSCHIEDE CDC-RECHNUNG UND DIESER (CA.1%).
 1.12 C CDC-VERSION IN VORBEREITUNG.          STIEMER 18.1.1974
 1.31 T " WENN SIE EINE ERKLAERUNG MOECHTEN :<CTRL/P,W 1,RETURN> "!!
 1.32 T " HABEN SIE AUCH UND NUR DIE FEATURES 3,4
 SONST NEUSTART PFOCAL!"!!
 1.35 T !!!" ENTFERNUNG ZWISCHEN ZWEI PUNKTEN AUF DEM ERDELLIPSOID",!!!
 2.01 E
 2.02 C ABFRAGEN
 2.10 A !!"STRECKE(1)
 GEKNICKTE BAHN(2)
 DREIECK(3)
"I1,!
 2.12 IF (I1-2)2.19,2.18,2.18
 2.18 A "LAENGE "LG(1),LM(1),LS(1)," BREITE "BG(1),BM(1),BS(1),!
 2.19 A "LAENGE "LG(2),LM(2),LS(2)," BREITE "BG(2),BM(2),BS(2),!
 2.20 A "LAENGE "LG(3),LM(3),LS(3)," BREITE "BG(3),BM(3),BS(3),!!
 3.01 C ENTSCHEIDUNGEN WELCHER KURS
 3.10 I (I1-2) 4.01,6.50,7.01
 4.01 C STRECKE
 4.04 S I2=1.5
 4.10 S J=2;D 8;S L1=L(J);S B1=B(J);S J=3;D 8;S L2=L(J);S B2=B(J)
 4.12 D 9;D 10;D 12;T !;S S8=S9;S P8=S8;S I=2;D 13.25;S I=1;D 13.15;D 13.16
 4.13 S I=3;D 13.25;S I=2;D 13.15
 4.90 GOTO 2.01
 5.01 C ZIEL IN GEKNICKTER BAHN
 5.10 S J=1;D 8;S L1=L(J);S B1=B(J)
 5.11 S J=2;D 8;S L2=L(J);S B2=B(J)
 5.12 S J=3;D 8;S L3=L(J);S B3=B(J)
 5.15 D 9;D 10;S S8=S9
 5.17 S L1=L2;S L2=L3;S B1=B2;S B2=B3
 5.19 D 9;D 10;S S7=S9;S P8=S8;S P7=S7;S PS=P8+P7
 5.20 D 13
 5.25 I (I1-2)20.01,6.05,7.01
 6.02 D 13
 6.05 GOTO 20.01
 6.50 S I2=1.5
 6.51 GOTO 5.01
 7.01 C DREIECK
 7.04 S I2=2
 7.10 D 5;S P8=2*P8;S P7=2*P7
 7.11 S L1=L(3);S B1=B(3);S L2=L(1);S B2=B(1)
 7.18 D 9;D 10;S S6=S9;S P6=2*S6
 7.20 S I=4;D 13.25
 7.25 T %7.1,"                                       ",S6,"   ",P6,!
 7.30 S I=1;D 13.25;S I=4;D 13.15;S I2=1
 7.31 GOTO 20.01
 8.01 C UMRECHNUNG IN BOGENMASS
 8.10 S LU(J)=LM(J)/60;S LV(J)=LS(J)/3600
 8.12 S L(J)=LG(J)+LU(J)+LV(J)
 8.15 S L(J)=L(J)*(3.1415926536/180)
 8.20 S BU(J)=BM(J)/60;S BV(J)=B(J)/3600
 8.22 S B(J)=BG(J)+BU(J)+BV(J)
 8.25 S B(J)=B(J)*(3.1415926536/180)
 9.01 C KONSTANTEN
 9.07 S X=B;D 14;S CB=X
 9.08 S X=B;D 15;S SB=X
 9.09 S C5=0.006768469
 9.10 S B=(B1+B2)/2;S V2=1+C5*CB
 9.12 S LD=(L2-L1);S BD=(B2-B1);S N=6400-22*CB
 9.15 S C6=SB;S D5=LD^2;S D6=BD^2
10.01 C ENTFERNUNG ZWEIER PUNKTE
10.10 S S1=(N/V2)*BD*<1-((2+C6)/24)*D5>
10.15 S S2=N*CB*LD*<1+D6/24-(C6/24)*D5>
10.20 S S9=FSQT(S1^2+S2^2)
11.01 C PUNKTEBERECHNUNG
12.01 C AUSDRUCK
12.10 T !!"  NR.         LAENGE      BREITE        STRECKE      PUNKTE"!
13.01 C AUSDRUCK
13.10 D 12.10
13.11 T !
13.14 S I=1;D 13.25
13.15 T %2,I,"         "A1,A2,A3,"  ",A4,A5,A6,!
13.16 T %7.1,"                                       ",S8,"   ",P8*I2,!
13.17 S I=2;D 13.25;D 13.15
13.19 T %7.1,"                                       ",S7,"   ",P7*I2,!
13.20 S I=3;D 13.25;D 13.15
13.25 S A1=LG(I);S A2=LM(I);S A3=LS(I);S A4=BG(I);S A5=BM(I);S A6=BS(I)
14.01 C COS-FUNKTION
14.05 S X=1-X^2/2+X^4/24-X^6/720+X^8/(720*56)
15.01 C SIN-FUNKTION
15.05 S X=X-X^3/6+X^5/120-X^7/6040+X^9/44880
20.01 C SCHLUSS
20.05 S SP=I2*(P6+P7+P8);S SS=S6+S7+S8
20.06 T %7.1,!,"                                 SUMME ",SS,"   ",SP,!!
20.30 GOTO 2.01
25.10 O O LPT:,
25.11 W A



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search