File TIM.

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

PROGRAM PARABELNAHEEPHEMERIDEN (IN,OUT);

CONST 
K=0.01720209895;PI=3.1415926536;

VAR 
X,Y,Z,M,R,V,F,E,PX,PY,PZ,QX,QY,QZ,I,W,L,EKL,S1,S2,S3,RA,DEC,DELTA:REAL;
H0,H1,H2,H3,H4,Q,C3,CHI,E0,A,T,GR,TAU,CE,SE:REAL;
D,MO,J:INTEGER;

FUNCTION SGN (VAR X:REAL):INTEGER;
BEGIN
  IF X=0 THEN SGN:=0;
  IF X>0 THEN SGN:=1 ELSE SGN:=-1;
END;

FUNCTION ARCSIN (VAR X:REAL):REAL;
BEGIN
  ARCSIN:=ARCTAN (SQRT (1/(SQR (X)-1)))*SGN (X);
END;

FUNCTION TAN (X:REAL):REAL;
BEGIN
  TAN:=SIN (X)/COS (X);
END;

FUNCTION F3 (VAR X:REAL):REAL;
VAR
I:INTEGER;
H,F,X2:REAL;
BEGIN
H:=1/6;F:=H;
I:=4;
REPEAT H:=H*X2/I/(I+1);
       F:=F+H;
       I:=I+2;
UNTIL  ABS (H/F)<1E-8;
F3:=F;
END;

PROCEDUR XYZ (VAR D,M,J:INTEGER);
VAR
  X,Y,Z,J1,D1,E,E1,F,MA,V,R,S,T,OM,EKL:REAL;

BEGIN
  J1:=J-1900;S:=TRUNC ((J1-1)/4);
  D1:=S+TRUNC (30.6*M+0.53/SQR (M-1.55)-32.3)+D-0.5;
  IF (J MOD 4=0) AND (M>=3) THEN D1:=D1+1;
  T:=(365*J1+D1)/36525;
  OM:=281.220833+1.719175*T+0.000361*T*T;
  F:=0.01675104-0.0000418*T;
  EKL:=23.452294-0.0130125*T;
  MA:=-1.524155-0.00015*T*T-0.25590255*J1+0.98560027*D1;
  MA:=MA*PI/180;EKL:=EKL*PI/180;OM:=OM*PI/180;
  E1:=MA;
    REPEAT E:=E1;
           E1:=M+F*SIN (E);
    UNTIL  ABS (E-E1)<1E-7;
  R:=1-F*COS (E);
  V:=2*ARCTAN (SQRT ((1+F)/(1-F))*TAN (E/2));
  V:=V+OM;
  X:=R*COS (V);
  Y:=R*SIN (V)*COS (EKL);
  Z:=R*SIN (V)*SIN (EKL);
END;

(*================Hauptprogramm================*)
BEGIN
WRITE("TAG");READLN;READ(D);
WRITE("MONAT");READLN;READ(MO);
W



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