Die Programme für ATARI ST Computer laufen mit dem ATARI STEmulator  (mit ARARI  TOS ROM 2.06) unter Windows 9x/ME (s. Abschnitt Kompatibilität). Kopieren Sie die ATARI ST Progr. aus dem Ordner GFA/GFA_PROG. und den GFA-BASIC-Interpreter auf das Festplatten-Laufwerk >C<.

REM ATARI ST GFA-BASIC46  ASTROGRAPH II FÜR GROSSMASSTÄBIGE FELDSTERNAUFNAHMEN
DEFFN gzb(x) = ATN(1 / (1 + 0.006739501819 * (6378140 / (6378140 + eh))) * TAN(x))
DEFFN ro(x) = 6356755.288157 / SQR(1 - 0.006694384999591 * COS(ATN(0.993305615 * TAN(x))) ^ 2) + eh
DEFFN r(x) = x - INT(x / (2 * PI)) * (2 * PI)
DIM ar(20),dek(20),ar1(20),dek1(20),x(20),y(20),x1(20),y1(20),ko(15),ko1(15),koo(15),p1(15,15),p(15,15),xo(20),yo(20),n(20),aa(20),dd(20),w1(20),w2(20),w3(20)
GOSUB astrograph
PROCEDURE astrograph
  REM -- ASTROGRAPH II - PHOTOGRAPHISCHE ORTSBESTIMMUNG AUF GROSSFELDAUFNAHMEN --
  k1 = RAD(20.49552 / 3600)
  neu:
  CLS
  PRINT AT(1,1);"- ASTRONOMISCHE ORTSBESTIMMUNG II -"
  PRINT AT(1,3);"DATUM (TT.MM.JJJJ)...........................: ";
  REM  INPUT "",a$
  a$="1.1.2000" !EINTRAG TESTDATEN-------------------
  a1 = ABS(LEN(a$) - 2)
  a2 = VAL(LEFT$(a$,2))
  a3 = VAL(MID$(a$,4 + (a2 < 10),2))
  a4 = VAL(RIGHT$(a$,a1 + (a2 < 10) + (a3 < 10) + (a2 > 9) * 2 + (a3 > 9) * 2))
  IF a2 = 0 OR a1 > 10 OR ABS(a2) > 31 OR ABS(a3) > 12 THEN
    GOTO neu
  ENDIF
  GOSUB jd
  PRINT AT(1,4);"UHRZEIT (UT AUFNAHMEMITTE h,m,s).............: ";
  REM  INPUT "",h,m,s
  h = 12
  m=0 !-------------------------
  s = 0
  ut = h + m / 60 + s / 3600
  PRINT AT(1,5);"KORREKTIONSWERT FÜR EPHEMERIDENZEIT (SEK)....: ";
  REM  INPUT "",ddt  !ddt=TDT-UT
  ddt=0 !---------------
  ddt = ddt / 3600
  PRINT AT(1,6);"HÖHE ÜBER N.N. (METER).......................: ";
  REM   INPUT "",eh
  eh = 0
  PRINT AT(1,7);"GEOGRAPHISCHE BREITE (o,','')................: ";
  REM   INPUT "",h$,m,s
  h$ = "50"
  m=0 !----------------------
  s = 0
  GOSUB grad        !KOORDINATEN BEOBACHTRUNGSSTATION
  br=FN gzb(g)      !FÜR TÄGL. ABERRATION
  p1 = FN ro(g) / 6378140
  PRINT AT(1,8);"GEOGRAPHISCHE L¥NGE (o,','').................: ";
  REM   INPUT "",h$,m,s
  h$ = "6"
  m=0 !---------------------
  s = 0
  GOSUB grad
  lgeo = DEG(g)
 jdbeo=jd+(ut+ddt)/24 !JULIAN. DATUM DER AUFNAHMEMITTE
  PRINT AT(1,9);"KATALOGEPOCHE REFERENZSTERNE (1950.0, 2000.0): ";
  REM  INPUT "",ep
  ep=2000 !-----------------
  PRINT AT(1,10);"BESSEL/FK4 ODER JULIAN./FK5  B,J.............: ";
  REM   INPUT "",ep$
  ep$="J" !------------------------
  ep$=UPPER$(ep$)
  IF ep$ = "B" THEN
    jdep = 2415020.31352 + 365.242198781 * (ep - 1900)
  ENDIF
  IF ep$ = "J" THEN
    jdep = 2451545 + 365.25 * (ep - 2000)
  ENDIF
  t1 = (jd - 2451545) / 36525
  t = (jdbeo - 2451545) / 36525
  GOSUB nut
  GOSUB osz
  jdd=(jdbeo-jdep)/365.25 !DIFFERENZ AUFNAHMENITTE-EPOCHE IN JULIAN. JAHREN
  PRINT AT(1,11);"PRIMÄRBRENNWEITE DER AUFNAHMEOBTIK (mm)......: ";  !IN;MILIIMETER
  REM  INPUT "",f
  fo=5000 !------------------
  PRINT AT(1,12);"ANZAHL REFERENZSTERNE........................: ";
  REM  INPUT "",anz
  anz=10 !---------------
  CLS
  i = 1
  PRINT AT(1,2);"KATALOG-REKTASZENSION STERN AUFNAHMEMITTE (h,m,s): ";  !AR;DER;OPTISCHEN;ACHSE
  REM  INPUT "",h$,m,s
  h$ = "0"
  m=20 !----------------------
  s = 0
  GOSUB grad
  ar(i) = g * 15
  REM EIGENBEWEGUNG BEZOGEN AUF DIESELBE EPOCHE WIE AR UND DEKLIN.
  PRINT AT(1,3);"EIGENBEWEGUNG STERN AR AUFNAHMEMITTE (ZEITSEK.)..: "; !FALLS STERN INMIITEN DER AUFNAHME SONST;pm = 0
  REM   INPUT "",pm
  pm=0 !--------------------
  ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
  PRINT AT(1,4);"KATALOG-DEKLINATION STERN AUFNAHMEMITTE (o,','').: ";  !DEKLIN. OPTISCHE ACHSE
  REM   INPUT "",h$,m,s
  h$ = "-15"
  m=40 !-----------------
  s = 0
  GOSUB grad
  dek(i) = g
  PRINT AT(1,5);"EIGENBEWEGUNG STERN DEKLIN. AUFNAHMEMITTE ('')...: ";
  REM   INPUT "",pm
  pm=0 !------------------
  dek(i) = dek(i) + RAD(pm / 3600) * jdd
  PRINT AT(1,6);"J¥HRL. PARALLAXE ('')............................: ";  !FALLS;STERN;INMITTEN;DER;AUFNAHME, SONST;p = 0
  REM   INPUT "",p
  p=0 !---------------------
  GOSUB appapos
  do = dek(i)
  aro = ar(i)
  CLS
  GOTO jump !----------------
  FOR i=1 TO anz  !EINGABE REFERENZSTERNE
  PRINT AT(1,2);i;") KATALOG-REKTASZENSION REFERENZSTERN (h,m,s): ";
  INPUT "",h$,m,s
  GOSUB grad
  ar(i) = g * 15
  REM EIGENBEWEGUNG BEZOGEN AUF DIESELBE EPOCHE WIE AR UND DEKLIN.
  PRINT AT(1,3);i;") EIGENBEWEGUNG AR REFERENZSTERN (ZEITSEK.)..: ";
  INPUT "",pm
  ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
  PRINT AT(1,4);i;") KATALOG-DEKLINATION REFERENZSTERN (o,','').: ";
  INPUT "",h$,m,s
  GOSUB grad
  dek(i) = g
  PRINT AT(1,5);i;") EIGENBEWEGUNG DEKLIN. REFERENZSTERN ('')...: ";
  INPUT "",pm
  dek(i) = dek(i) + RAD(pm / 3600) * jdd
  PRINT AT(1,6);i;") J¥HRL. PARALLAXE ('')......................: ";
  INPUT "",p
  GOSUB appapos
  PRINT AT(1,7);i;") x-,y-MESSUNG REFERENZSTERN (mm)............: "; !IN MILLIMETER
  INPUT "",x(i),y(i)
  REM STRENGE STANDARDKOORDINATEN X,Y
  j = SIN(dek(i)) * SIN(do) + COS(dek(i)) * SIN(do) * COS(ar(i) - aro)
  x1(i) = (COS(dek(i)) * SIN(ar(i) - AR)) / j
  y1(i) = (SIN(dek(i)) * COS(do) - COS(dek(i)) * SIN(do) * COS(ar(i) - aro)) / j
  CLS
NEXT i
jump: !------------------
REM EINGABE REFERENZSTERNE
REM (BEI VERWENDUNG OBIGER EINGABESCHLEIFE DIESEN TESTTEIL LÖSCHEN)--
i=1 !NR. 1
h$="0"  !REKTASZENSION STUNDEN
m=7     !MINUTEN
s=46.923 !SEKUNDEN
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$="-22" !DEKLINATION
m=30     !BOGENMINUTEN
s=30.96  !BOGENSEKUNDEN
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=2 !NR. 2
h$ = "0"
m = 6
s = 24.982
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-18"
m = 2
s = 17.81
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=3 !NR. 3
h$ = "0"
m = 17
s = 8.883
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-15"
m = 39
s = 27.94
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=4 ! NR. 4
h$ = "0"
m = 34
s = 45.931
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-24"
m = 31
s = 40.79
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=5 ! NR. 5
h$ = "0"
m = 30
s = 45.105
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-17"
m = 42
s = 54.03
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=6 ! NR. 6
h$ = "0"
m = 43
s = 35.372
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-17"
m = 59
s = 11.82
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ---------------------------
i = 7
h$ = "0"
m = 30
s = 22.661
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-23"
m = 47
s = 15.72
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i = 8
h$ = "0"
m = 24
s = 49.643
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-18"
m = 28
s = 18.66
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=9 !NR. 9
h$ = "0"
m = 18
s = 16.954
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-21"
m = 8
s = 18.22
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=10 ! NR. 10
h$ = "0"
m = 20
s = 2.973
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-17"
m = 42
s = 1.67
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=11 ! NR. 11
h$ = "0"
m = 15
s = 46
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-28"
m = 40
s = 0
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
i=12 ! NR. 12
h$ = "2"
m = 43
s = 35
GOSUB grad
ar(i) = g * 15
pm = 0
ar(i) = FN r(ar(i) + RAD((pm * 15) / 3600) * jdd)
h$ = "-5"
m = 59
s = 11.82
GOSUB grad
dek(i) = g
pm = 0
dek(i) = dek(i) + RAD(pm / 3600) * jdd
p = 0
GOSUB appapos
GOSUB messung
x(i) = x
y(i) = y
REM ----------------------------------------------
REM STRENGE STANDARDKOORDINATEN X,Y AUSGLEICHUNG BEOBACHTUNG MINUS RECHNUNG
FOR i = 1 TO anz
  j = SIN(dek(i)) * SIN(do) + COS(dek(i)) * COS(do) * COS(ar(i) - aro)
  x1(i) = (COS(dek(i)) * SIN(ar(i) - aro)) / j
  y1(i) = (SIN(dek(i)) * COS(do) - COS(dek(i)) * SIN(do) * COS(ar(i) - aro)) / j
NEXT i
REM ---------------------
GOSUB ausgleich
CLS
PRINT AT(1,1);"BESTIMMUNG DER PROJIZIERTEN/REPRODUZIERTEN NEGATIVVERGRÖSSERUNG
PRINT AT(1,2);"NR. ZWEIER REFERENZSTERNE GEMESS. WINKELDISTANZ A,B: ";
REM  INPUT "",nr1,nr2
nr1 = 2
nr2=4 !---------------
PRINT AT(1,3);"GEMESSENE WINKELDISTANZ STERN NR. ";nr1;" ZU ";nr2;" (mm)......: "; !MILLIMETER
REM  INPUT "",dis
dis=101.65 !--------------
d = ACOS(SIN(dek(nr1)) * SIN(dek(nr2)) + COS(dek(nr1)) * COS(dek(nr2)) * COS(ar(nr2) - ar(nr1)))
d=DEG(d)*3600 !DISTANZ DER ZWEI STERNE IN BOGENSEKUNDEN
f1=fo/206264.8062471 !1'' (BOGENSEK.) AM HIMMEL = f1 MILLIMETER AUF DEM NEGATIV
v=dis/(d*f1) !VERGRÖSSERUNG DER DIAPROJEKTION/KOPIE GEGENÜBER DEM NEGATIV
r=fo*v !EFFEKTIVE BRENNWEITE IN MILLIMETER
f2=r/206264.8062471 !MASSTAB 1'' AUF DER DIAPROJEKTION/KOPIE = f2 MILLIMETER
f3=1/f2 !1 mm DER DIAPROJEKTION/KOPIE = f3 BOGENSEKUNDEN
CLS
PRINT AT(1,1);"ANZAHL MESSUNGEN OBJEKT: ";  !GERADE MEHRFACH HIN UND ZURÜCK MESSEN
REM  INPUT "",ns
ns=1 !------------------
PRINT AT(1,2);SPC(50)
j = 0
GOTO ik !----------------
FOR i = 1 TO ns
  j = j + 1
  PRINT AT(1,j);j;") MESSUNG OBJEKT x,y (mm): ";  !EINGABE;MILLIMETER
  INPUT "",x(i),y(i)
NEXT i
ik: !-------
x(ns) = ROUND(x(1),1)
y(ns) = ROUND(y(1),1)
FOR i = 1 TO ns
  xo(i)=koo(1)+koo(2)*x(i)+koo(3)*y(i)+koo(4)*x(i)*y(i)+koo(5)*x(i)*x(i)+koo(6)*y(i)*y(i)+koo(7)*x(i)*y(i)*y(i)+koo(8)*x(i)*x(i)*y(i)+koo(9)*x(i)*x(i)*x(i)+k oo(10)*y(i)*y(i)*y(i)  !AUSGEGLICHENE STANDARDKOORDINATEN
  yo(i)=ko1(1)+ko1(2)*x(i)+ko1(3)*y(i)+ko1(4)*x(i)*y(i)+ko1(5)*x(i)*x(i)+ko1(6)*y(i)*y(i)+ko1(7)*x(i)*y(i)*y(i)+ko1(8)*x(i)*x(i)*y(i)+ko1(9)*x(i)*x(i)*x(i)+ko 1(10)*y(i)*y(i)*y(i)  !DES OBJEKTS
  j=COS(do)-yo(i)*SIN(do)
  ar1(i)=FN r(aro+ATN(xo(i)/j))  !REKTASZENSION OBJEKT
  dek1(i)=ATN((SIN(do)+yo(i)*COS(do))/SQR(xo(i)^2+j^2)) !DEKLINATION OBJEKT
NEXT i
FOR i=1 TO ns
  n(i)=xo(i)
NEXT i
GOSUB sig
PRINT AT(1,1);"AUSGEGLICHENE STANDARDKOORDINATEN DES OBJEKT:"
PRINT AT(1,2);USING "X: ##.######### ",mw;
PRINT "/ m.F. ";mf;"/ Streuung ";sig
FOR i=1 TO ns
  n(i)=yo(i)
NEXT i
GOSUB sig
PRINT AT(1,3);USING "Y: ##.######### ",mw;
PRINT "/ m.F. ";mf;"/ Streuung ";sig
FOR i=1 TO ns
  n(i)=DEG(ar1(i))/15
NEXT i
GOSUB sig
arap=RAD(mw)*15
GOSUB sexa
IF ep$="J" THEN
  PRINT AT(1,4);"SCHEINBARER ORT OBJEKT, WAHRES ÄQUINOKTIUM DES DATUM/FK5:"
ELSE
  PRINT AT(1,4);"SCHEINBARER ORT OBJEKT, WAHRES ÄQUINOKTIUM DES DATUM/FK4:"
ENDIF
PRINT AT(1,5);USING "AR.....: ## h ## m ##.### s",h,m,s;
PRINT "/ m.F. ";mf;"/ Streuung ";sig
FOR i=1 TO ns
  n(i)=DEG(dek1(i))
NEXT i
GOSUB sig
dekap=RAD(mw)
GOSUB sexa
PRINT AT(1,6);"DEKLIN.: ";z$;USING "## o ## ' ##.## ''",h,m,s;
PRINT "/ m.F. ";mf;"/ Streuung ";sig
GOSUB app  !SUBROUTINE MITTL. ORT
arj=ar !MITTL. AR
dj=d   !MITTL. DEKLIN.
mw=DEG(AR)/15
GOSUB sexa
IF ep$="J" THEN
  PRINT AT(1,7);"ASTROMETR. ORT OBJEKT, MITTL. ÄQUINOKTIUM J2000/FK5"
ELSE
  PRINT AT(1,7);"ASTROMETR. ORT OBJEKT, MITTL. ÄQUINOKTIUM B1950/FK4"
ENDIF
PRINT AT(1,8);USING "AR.....:  ## h ## m ##.### s",h,m,s
mw=DEG(d)
GOSUB sexa
PRINT AT(1,9);"DEKLIN.: ";z$;USING "## o ## ' ##.## ''",h,m,s
IF ep$="J" THEN
  GOSUB conv1
ELSE
  GOSUB conv
ENDIF
mw=DEG(AR)/15
GOSUB sexa
IF ep$="J" THEN
  PRINT AT(1,10);"ASTROMETR. ORT OBJEKT, MITTL. ÄQUINOKTIUM B1950/FK4"
ELSE
  PRINT AT(1,10);"ASTROMETR. ORT OBJEKT, MITTL. ÄQUINOKTIUM J2000/FK5"
ENDIF
PRINT AT(1,11);USING "AR.....: ## h ## m ##.### s",h,m,s
mw=DEG(d)
GOSUB sexa
PRINT AT(1,12);"DEKLIN.: ";z$;USING "## o ## ' ##.## ''",h,m,s
PRINT AT(1,16);"PRIMÄRBRENNWEITE DER AUFNAHMEOPTIK: ";fo;" mm"
PRINT AT(1,17);"EFFEKTIVE BRENNWEITE: ";r;" mm"
PRINT AT(1,18);"VERGRÖSSERUNGSFAKTOR: ";v;"x"
PRINT AT(1,19);"MASSTAB: 1''  = ";f1;" mm AUF DEM NEGATIV"
PRINT AT(1,20);"MASSTAB: 1''  = ";f2;" mm AUF DER KOPIE/DIAPROJEKTION"
PRINT AT(1,21);"MASSTAB: 1 mm = ";f3;" '' AUF DER KOPIE/DIAPROJEKTION"
PRINT AT(1,23);"TASTE >1<"
PRINT AT(1,24);"TASTE >2<"
PRINT AT(1,25);"TASTE >S<"
DO
  w$=UPPER$(INKEY$)
  IF w$="1" THEN
    GOSUB tafel
    IF w$="1" THEN
      GOSUB tafel
    ENDIF
  ENDIF
  IF w$="2" THEN
    GOSUB astrographZ
  ENDIF
  EXIT IF w$="S"
LOOP
RETURN
PROCEDURE grad
  j=1
  IF LEFT$(h$,1)="-" THEN
    j=-1
  ENDIF
  g=RAD(ABS(VAL(h$))+m/60+s/3600)*j
RETURN
PROCEDURE sexa
  IF mw<0 THEN
    z$="-"
  ELSE
    z$="+"
  ENDIF
  g=ABS(mw)
  h=FIX(g)
  m=FIX((g-h)*60)
  s=((g-h)*60-m)*60
RETURN
PROCEDURE tafel
  CLS
  PRINT AT(15,1);"       EINMESSUNG EINES KOORDINATENKREUZES    "
  PRINT AT(15,2);"STRENGE STANDARDKOORDINATEN X,Y DER REFERENZSTERNE"
  PRINT AT(15,3);"           WAHRS ÄQUINOKTIUM DES DATUMS       "
  u=4
  FOR i=1 TO anz
    u=u+1
    PRINT AT(1,u);i;") X: ";-r*x1(i);" mm, Y: ";r*y1(i);" mm"
  NEXT i
  PRINT AT(1,24);"TASTE >W<"
  REPEAT
  UNTIL UPPER$(INKEY$)="W"
  CLS
  PRINT AT(15,1);"SCHEINBARE REKTASZENSION UND DEKLINATION"
  PRINT AT(15,2);"     WAHRES ÄQUINOKTIUM DES DATUMS"
  u=3
  FOR i=1 TO anz
    u=u+1
    mw=DEG(ar(i))/15
    GOSUB sexa
    PRINT AT(1,u);i;USING ") AR:  ## h ## m ##.### s",h,m,s;
    mw=DEG(dek(i))
    GOSUB sexa
    PRINT ",  DEKLIN.: ";z$;USING "## o ## ' ##.## ''",h,m,s
  NEXT i
  PRINT AT(1,23);"TASTE >1<"
  PRINT AT(1,24);"TASTE >2<"
  REPEAT
  UNTIL UPPER$(INKEY$)="1" OR UPPER$(INKEY$)="2"
  DO
    w$=UPPER$(INKEY$)
    IF w$="1" OR w$="2" THEN
      GOTO j
    ENDIF
  LOOP
  j:
RETURN
PROCEDURE jd
  jd=1720994.5+a2+FIX(30.6001*((a3-12*(a3<3))+1))+FIX(f*365.245*(a4+(a3<3)))
RETURN
PROCEDURE ausgleich
  REM AUSGLEICHSRECHNUNG NACH DER METHODE DER KLEINSTEN QUADRATE
  REM ----- NORMALGLEICHUNG MATRIX 10 X 10 ----------
  rr=0
  br=0
  cr=0
  dr=0
  er=0
  fr=0
  gr=0
  hr=0
  ir=0
  jr=0
  FOR i=1 TO anz
    rr=rr+x1(i)
    br=br+x(i)*x1(i)
    cr=cr+y(i)*x1(i)
    dr=dr+x(i)*y(i)*x1(i)
    er=er+x(i)*x(i)*x1(i)
    fr=fr+y(i)*y(i)*x1(i)
    gr=gr+x(i)*y(i)*y(i)*x1(i)
    hr=hr+x(i)*x(i)*y(i)*x1(i)
    ir=ir+x(i)*x(i)*x(i)*x1(i)
    jr=jr+y(i)*y(i)*y(i)*x1(i)
  NEXT i
  GOSUB mess
  m=10
  n=10
  FOR i=1 TO m
    FOR k=1 TO n+1
      p(i,k)=p1(i,k)
    NEXT k
  NEXT i
  GOSUB elim
  FOR i=1 TO m
    koo(i)=ko(i)
  NEXT i
  i=1
  GOSUB gl
  PRINT
  PRINT "TASTE >W<"
  REPEAT
  UNTIL UPPER$(INKEY$)="W"
  CLS
  rr=0
  br=0
  cr=0
  dr=0
  er=0
  fr=0
  gr=0
  hr=0
  ir=0
  jr=0
  FOR i=1 TO anz
    rr=rr+y1(i)
    br=br+x(i)*y1(i)
    cr=cr+y(i)*y1(i)
    dr=dr+x(i)*y(i)*y1(i)
    er=er+x(i)*x(i)*y1(i)
    fr=fr+y(i)*y(i)*y1(i)
    gr=gr+x(i)*y(i)*y(i)*y1(i)
    hr=hr+x(i)*x(i)*y(i)*y1(i)
    ir=ir+x(i)*x(i)*x(i)*y1(i)
    jr=jr+y(i)*y(i)*y(i)*y1(i)
  NEXT i
  GOSUB mess
  m=10
  n=10
  FOR i=1 TO m
    FOR k=1 TO n+1
      p(i,k)=p1(i,k)
    NEXT k
  NEXT i
  GOSUB elim
  FOR i=1 TO m
    ko1(i)=ko(i)
  NEXT i
  i=0
  GOSUB gl
  PRINT
  PRINT "TASTE >W<"
  REPEAT
  UNTIL UPPER$(INKEY$)="W"
RETURN
PROCEDURE elim
  FOR j=1 TO n-1     !GAUSS ELIMINATION
  nr=j
  no=ABS(p(j,j))
  FOR i=j+1 TO n    STEP !ZEILENPIVOT
    noo=ABS(p(i,j))
    EXIT IF (noo-no)<0
    no=noo
    nr=i
  NEXT i
  IF nr=j THEN
    GOTO jum1
  ENDIF
  FOR i=j TO m+1
    no=p(nr,i)
    p(nr,i)=p(j,i)
    p(j,i)=no
  NEXT i
  jum1:
  FOR i=j+1 TO m+1   STEP !ELIMINATION
    p(j,i)=p(j,i)/p(j,j)
  NEXT i
  FOR i=j+1 TO n
    FOR k=j+1 TO m+1
      p(i,k)=p(i,k)-p(j,k)*p(i,j)
    NEXT k
  NEXT i
NEXT j
ko(n)=p(n,n+1)/p(n,n)  !RÜCKSUBSTITUTION
FOR j=n-1 TO 1 STEP -1
  ko(j)=p(j,n+1)
  FOR i=j+1 TO n
    ko(j)=ko(j)-p(j,i)*ko(i)
  NEXT i
NEXT j
RETURN
PROCEDURE nut
  REM DELAUNAY-ELEMENTE
  d = FN r(5.198468 + 7771.37714 * t)      !Mittl. Elongation Mond-Sonne
  ll = FN r(6.24006 + 628.3019555 * t)      !Mittl. Anomalie Sonne (Erde)
  l = FN r(2.3555482 + 8328.691422 * t)   !Mittl. Anomalie Mond
  f = FN r(1.627903 + 8433.46615813 * t) !Mittl. Länge Mond ab aufsteigenden Mondknoten
  af = FN r(2.182438 - 33.757045 * t)        !Mittl. Länge aufst. Mondknoten
  REM NUTATION IN LÄNGE
  nu = -172000 * SIN(af) + 2060 * SIN(2 * af)  - 1320 * SIN(2 * f - 2 * d + 2 * af) + 1430 * SIN(ll)
  nu = nu / 10000
  REM NUTATION IN SCHIEFE
  nu1 = 92000 * COS(af) - 895 * COS(2 * af) + 5740 * COS(2 * f - 2 * d + 2 * af) + 220 * COS(ll + 2 * f - 2 * d + 2 * af)
  nu1 = nu1 / 10000
   f7=0.409093-0.000227*t !MITTL. EKLIPTISCHIEFE
  ec=f7+RAD(nu1/3600) !WAHRE EKLIPTIKSCHIEFE
  REM lms,ms,ls = geometr. mittl. L„nge der Sonne. mittl. Anomalie der Sonne, geometr. wahre L„nge der Sonne
  lms=FN r(4.89506+628.33197*t)
  ms=FN r(6.24006+628.30196*t)
  c=0.03342*SIN(ms)+0.00035*SIN(2*ms)
  l=FN r(lms+c)
  v=FN r(ms+c)
  REM ds, ars = Deklin. und Rektaszension der Sonne Äquinoktium des Datums JD
  ds=ASIN(SIN(ec)*SIN(l))  !Deklin. Sonne
  ys=SIN(l)*COS(ec)
  xs=COS(l)
  r1=SQR(xs^2+ys^2)
  x1=xs/r1
  y1=ys/r1
  ars=FN r(ATN(y1/(1+x1))*2) !Rektaszension Sonne
  e=0.0167086
  pe=RAD(102.94+0.71953*t)
 
rs = 10 ^ (((3040 - 15 * t) + (-727600 + 1810 * t) * COS(ms) + (-9090 + 50 * t) * COS(2 * ms)) / 100000000) //Radius Vektor Sonne
RETURN
PROCEDURE appapos
  d1=dek(i)
  a1=ar(i)
  t2=(jdbeo-jdep)/36525
  IF ep$="J" THEN
    GOSUB praez
  ELSE
    t2=(jdbeo-jdep)/36524.2198782
    GOSUB praez1
  ENDIF
  GOSUB reduk
  dek(i)=d+ded+dpa+dnu+dab+di
  ar(i)=FN r(AR+ard+apa+anu+aab+ai)
RETURN
PROCEDURE app
  d=dekap
  AR=arap
  p=0
  GOSUB reduk
  dekm=dekap-(dpa+ded+dnu+dab+di)
  arm=FN r(arap-(apa+ard+anu+aab+ai))
  REM ----- PRÄZESSION ----------
  t2=(jdep-jdbeo)/36525
  d1=dekm
  a1=arm
  IF ep$="J" THEN
    GOSUB praez
  ELSE
    t2=(jdep-jdbeo)/36524.2198782
    GOSUB praez1
  ENDIF
RETURN
PROCEDURE praez
  REM ----- PRÄZESSION FK5 ----------
  to=(jdep-2451545)/36525
  w1=RAD(((2306.218+1.4*to)*t2+0.3*t2*t2)/3600)
  w2=RAD(((2306.218+1.4*to)*t2+1.1*t2*t2)/3600)
  w3=RAD(((2004.311-0.85*to)*t2-0.43*t2*t2)/3600)
  d=ASIN(SIN(w3)*COS(d1)*COS(a1+w1)+COS(w3)*SIN(d1))
  x=(COS(w3)*COS(d1)*COS(a1+w1)-SIN(w3)*SIN(d1))/COS(d)
  y=(COS(d1)*SIN(a1+w1))/COS(d)
  AR=FN r(ATN(y/(1+x))*2+w2)
RETURN
PROCEDURE praez1
  REM ----- PRÄZESSION FK4----------
  to=(jdep-2433282.423459)/36524.2198781
  w1=RAD(((2304.952+1.4*to)*t2+0.3*t2*t2)/3600)
  w2=RAD(((2304.952+1.4*to)*t2+1.1*t2*t2)/3600)
  w3=RAD(((2004.259-0.85*to)*t2-0.43*t2*t2)/3600)
  d=ASIN(SIN(w3)*COS(d1)*COS(a1+w1)+COS(w3)*SIN(d1))
  x=(COS(w3)*COS(d1)*COS(a1+w1)-SIN(w3)*SIN(d1))/COS(d)
  y=(COS(d1)*SIN(a1+w1))/COS(d)
  AR=FN r(ATN(y/(1+x))*2+w2)
RETURN
PROCEDURE reduk
  REM -------- LICHTABLENKUNG ----------------
  REM elo = Elongationswinkel bzw. Winkeldistanz Stern-Sonne.
  abl=RAD(0.00407/3600)
  elo=ACOS(SIN(ds)*SIN(d)+COS(ds)*COS(d)*COS(AR-ars))
  ard=abl*COS(ds)*SIN(AR-ars)/(1-COS(elo)*COS(d))
  ded=abl*(SIN(d)*COS(ds)*COS(AR-ars)-COS(d)*SIN(ds))/(1-COS(elo))
  REM ---------- JÄHRLICHE PARALLAXE  ----------
  apa=RAD((rs*p*(COS(AR)*COS(ec)*SIN(l)-SIN(AR)*COS(l))*(1/COS(d)))/3600)
  dpa=RAD((rs*p*(COS(d)*SIN(ec)*SIN(l)-COS(AR)*SIN(d)*COS(l)-SIN(AR)*SIN(d)*COS(ec)*SIN(l)))/3600)
  REM ---------- NUTATION --------------------
  anu=(COS(ec)+SIN(ec)*SIN(AR)*TAN(d))*RAD(nu/3600)-COS(AR)*TAN(d)*RAD(nu1/3600)
  dnu=(SIN(ec)*COS(AR)*RAD(nu/3600)+SIN(AR)*RAD(nu1/3600))
  REM --------- JÄHRL. ABERRATION (FIXSTERNABERRATION) ------
  aab=-k1*((COS(AR)*COS(l)*COS(ec)+SIN(AR)*SIN(l))/COS(d))+e*k1*((COS(AR)*COS(pe)*COS(ec)+SIN(AR)*SIN(pe))/ COS(d))
 dab=-k1*(COS(l)*COS(ec)*(TAN(ec)*COS(d)-SIN(AR)*SIN(d))+COS(AR)*SIN(d)*SIN(l))+e*k1*(COS(pe)*COS(ec)*(TA N(ec)*COS(d)-SIN(AR)*SIN(d))+COS(AR)*SIN(d)*SIN(pe))
  REM ------- TÄGL. ABERRATION -----
  di=RAD(-0.319/3600)*p1*COS(br)*SIN(osz-AR)*SIN(d)
  ai=RAD(-0.319/3600)*((p1*COS(br)*COS(osz-AR))/COS(d))
RETURN
PROCEDURE sig
  k=0
  vv=0
  FOR i=1 TO ns
    k=k+n(i)
  NEXT i
  mw=k/ns  !ARITHMERTISCHER MITTELWERT
  FOR i=1 TO ns
    vv=vv+(n(i)-mw)^2 !FEHLERQUADRATSUMME
  NEXT i
  IF ns>1 THEN
    sig=SQR(vv/(ns-1)) !STANDARDABWEICHUNG
    mf=SQR(vv/(ns*(ns-1))) !MITTL. FEHLER DES MITTELSWERTES
  ENDIF
RETURN
PROCEDURE osz
  REM ORTSSTERNZEIT
  sz=6.697374558333+0.05133690722222*t1
  sz1=2400*t1
  sz1=sz1-INT(sz1/24)*24
  sz=sz+sz1
  sz=sz-INT(sz/24)*24
  sz=sz+ut*1.002737909+lgeo/15+((nu*COS(ec))/15)/3600
  sz=sz-INT(sz/24)*24
  osz=RAD(sz*15)
RETURN
PROCEDURE conv
  REM astrometrisch B1950/FK4 zu astrometrisch J2000/FK5
  x=COS(dj)*COS(arj)
  y=COS(dj)*SIN(arj)
  z=SIN(dj)
  REM Subktraktion der E-Terme
  x1=x-(-1.6E-06)+(x*(-1.6E-06))*x
  y1=y-(-3.2E-07)+(y*(-3.3E-07))*y
  z1=z-(-1.4E-07)+(z*(-1.4E-07))*z
  xo=0.99993*x1-0.01118*y1-0.00486*z1
  yo=0.01118*x1+0.99994*y1-0.00003*z1
  zo=0.00486*x1-0.00003*y1+0.99999*z1
  t=(jdbeo-2433282.423459)/36525
  x=xo+((-0.00265*x1-1.154*y1+2.11114*z1)/1000000)*t
  y=yo+((1.1541*x1-0.0129*y1+0.0236*z1)/1000000)*t
  z=zo+((-2.1113*x1-0.0056*y1+0.01026*z1)/1000000)*t
  ro=SQR(x^2+y^2+z^2)
  d=ASIN(z/ro)
  x=x/(ro*COS(d))
  y=y/(ro*COS(d))
  AR=FN r(ATN(y/(1+x))*2)
RETURN
PROCEDURE conv1
  REM astrometrisch J2000/FK5 zu astrometrisch B1950/FK4
  x=COS(dj)*COS(arj)
  y=COS(dj)*SIN(arj)
  z=SIN(dj)
  xo=0.99993*x+0.01118*y+0.00486*z
  yo=-0.01118*x+0.99994*y-0.00003*z
  zo=-0.00486*x-0.00003*y+0.99999*z
  t=(jdbeo-2451545)/36525
  x=xo+((-0.00265*x-1.154*y+2.11114*z)/1000000)*t
  y=yo+((1.1541*x-0.0129*y+0.0236*z)/1000000)*t
  z=zo+((-2.1113*x-0.0056*y+0.01026*z)/1000000)*t
  REM Addition der E-Terme
  x1=x+(-1.6E-06)-(x*(-1.6E-06))*x
  y1=y+(-3.2E-07)-(y*(-3.3E-07))*y
  z1=z+(-1.4E-07)-(z*(-1.4E-07))*z
  ro=SQR(x1^2+y1^2+z1^2)
  d=ASIN(z1/ro)
  x1=x1/(ro*COS(d))
  y1=y1/(ro*COS(d))
  AR=FN r(ATN(y1/(1+x1))*2)
RETURN
PROCEDURE gl
  IF i=1 THEN
    PRINT "PLATTENKONSTANTEN x:"
  ELSE
    PRINT "PLATTENKONSTANTEN y:"
  ENDIF
  PRINT "A: ";ko(1);"; B: ";ko(2);"; C: ";ko(3)
  PRINT "D: ";ko(4);"; E: ";ko(5);"; F: ";ko(6)
  PRINT "G: ";ko(7);"; H: ";ko(8);"; I: ";ko(9)
  PRINT "J: ";ko(10)
  PRINT "LINKE UND RECHTE SEITE NORMALGLEICHUNG:"
  PRINT anz*ko(1)+b1*ko(2)+c1*ko(3)+d1*ko(4)+e1*ko(5)+f1*ko(6)+g1*ko(7)+h1*ko(8)+i1*ko(9)+j1*ko(10),rr
  PRINT b1*ko(1)+bb*ko(2)+cb*ko(3)+db*ko(4)+eb*ko(5)+fb*ko(6)+gb1*ko(7)+hb*ko(8)+ib*ko(9)+jb*ko(10),br
  PRINT c1*ko(1)+bc*ko(2)+cc*ko(3)+dc*ko(4)+ecc*ko(5)+fc*ko(6)+gc*ko(7)+hc*ko(8)+ic*ko(9)+jc*ko(10),cr
  PRINT d1*ko(1)+bd*ko(2)+cd*ko(3)+dd*ko(4)+ed*ko(5)+fd*ko(6)+gd*ko(7)+hd*ko(8)+id*ko(9)+jd*ko(10),dr
  PRINT e1*ko(1)+be*ko(2)+ce*ko(3)+de*ko(4)+ee*ko(5)+fe*ko(6)+ge*ko(7)+he*ko(8)+ie*ko(9)+je*ko(10),er
  PRINT f1*ko(1)+bf*ko(2)+cf*ko(3)+df*ko(4)+ef*ko(5)+ff*ko(6)+gf*ko(7)+hf*ko(8)+if*ko(9)+jf*ko(10),fr
  PRINT g1*ko(1)+bg*ko(2)+cg*ko(3)+dg*ko(4)+eg*ko(5)+fg*ko(6)+gg*ko(7)+hg*ko(8)+ig*ko(9)+jg*ko(10),gr
  PRINT h1*ko(1)+bh*ko(2)+ch*ko(3)+dh*ko(4)+eh*ko(5)+fh*ko(6)+gh*ko(7)+hh*ko(8)+ih*ko(9)+jh*ko(10),hr
  PRINT i1*ko(1)+bi*ko(2)+ci*ko(3)+di*ko(4)+ei*ko(5)+fi*ko(6)+gi*ko(7)+hi*ko(8)+ii*ko(9)+ji*ko(10),ir
  PRINT j1*ko(1)+bj*ko(2)+cj*ko(3)+dj*ko(4)+ej*ko(5)+fj*ko(6)+gj*ko(7)+hj*ko(8)+ij*ko(9)+jj*ko(10),jr
RETURN
PROCEDURE mess
  b1=0  !PARAMETER AUF NULL SETZEN
  c1=0
  d1=0
  e1=0
  f1=0
  g1=0
  h1=0
  i1=0
  j1=0
  REM ----------
  bb=0
  cb=0
  db=0
  eb=0
  fb=0
  gb=0
  hb=0
  ib=0
  jb=0
  REM ------------
  bc=0
  cc=0
  dc=0
  ecc=0
  fc=0
  gc=0
  hc=0
  ic=0
  jc=0
  REM ------------
  bd=0
  cd=0
  dd=0
  ed=0
  fd=0
  gd=0
  hd=0
  id=0
  jd=0
  REM ------------
  be=0
  ce=0
  de=0
  ee=0
  fe=0
  ge=0
  he=0
  ie=0
  je=0
  REM ------------
  bf=0
  cf=0
  df=0
  ef=0
  ff=0
  gf=0
  hf=0
  if=0
  jf=0
  REM ------------
  bg=0
  cg=0
  dg=0
  eg=0
  fg=0
  gg=0
  hg=0
  ig=0
  jg=0
  REM ------------
  bh=0
  ch=0
  dh=0
  eh=0
  fh=0
  gh=0
  hh=0
  ih=0
  jh=0
  REM ------------
  bi=0
  ci=0
  di=0
  ei=0
  fi=0
  gi=0
  hi=0
  ii=0
  ji=0
  REM ------------
  bj=0
  cj=0
  dj=0
  ej=0
  fj=0
gj = 0
hj = 0
ij = 0
jj = 0
FOR i=1 TO anz  !AUFSUMMIERUNG PARAMETER
b1 = b1 + x(i)
c1 = c1 + y(i)
d1 = d1 + x(i) * y(i)
e1 = e1 + x(i) * x(i)
f1 = f1 + y(i) * y(i)
g1 = g1 + x(i) * y(i) * y(i)
h1 = h1 + x(i) * x(i) * y(i)
i1 = i1 + x(i) * x(i) * x(i)
j1 = j1 + y(i) * y(i) * y(i)
REM -------------
bb = bb + x(i) * x(i)
cb = cb + y(i) * x(i)
db = db + x(i) * y(i) * x(i)
eb = eb + x(i) * x(i) * x(i)
fb = fb + y(i) * y(i) * x(i)
gb1 = gb1 + x(i) * y(i) * y(i) * x(i)
hb = hb + x(i) * x(i) * y(i) * x(i)
ib = ib + x(i) * x(i) * x(i) * x(i)
jb = jb + y(i) * y(i) * y(i) * x(i)
REM -------------
bc = bc + x(i) * y(i)
cc = cc + y(i) * y(i)
dc = dc + x(i) * y(i) * y(i)
ecc = ecc + x(i) * x(i) * y(i)
fc = fc + y(i) * y(i) * y(i)
gc = gc + x(i) * y(i) * y(i) * y(i)
hc = hc + x(i) * x(i) * y(i) * y(i)
ic = ic + x(i) * x(i) * x(i) * y(i)
jc = jc + y(i) * y(i) * y(i) * y(i)
REM -------------
bd = bd + x(i) * x(i) * y(i)
cd = cd + y(i) * x(i) * y(i)
dd = dd + x(i) * y(i) * x(i) * y(i)
ed = ed + x(i) * x(i) * x(i) * y(i)
fd = fd + y(i) * y(i) * x(i) * y(i)
gd = gd + x(i) * y(i) * y(i) * x(i) * y(i)
hd = hd + x(i) * x(i) * y(i) * x(i) * y(i)
id = id + x(i) * x(i) * x(i) * x(i) * y(i)
jd = jd + y(i) * y(i) * y(i) * x(i) * y(i)
REM -------------
be = be + x(i) * x(i) * x(i)
ce = ce + y(i) * x(i) * x(i)
de = de + x(i) * y(i) * x(i) * x(i)
ee = ee + x(i) * x(i) * x(i) * x(i)
fe = fe + y(i) * y(i) * x(i) * x(i)
ge = ge + x(i) * y(i) * y(i) * x(i) * x(i)
he = he + x(i) * x(i) * y(i) * x(i) * x(i)
ie = ie + x(i) * x(i) * x(i) * x(i) * x(i)
je = je + y(i) * y(i) * y(i) * x(i) * x(i)
REM -------------
bf = bf + x(i) * y(i) * y(i)
cf = cf + y(i) * y(i) * y(i)
df = df + x(i) * y(i) * y(i) * y(i)
ef = ef + x(i) * x(i) * y(i) * y(i)
ff = ff + y(i) * y(i) * y(i) * y(i)
gf = gf + x(i) * y(i) * y(i) * y(i) * y(i)
hf = hf + x(i) * x(i) * y(i) * y(i) * y(i)
if = if + x(i) * x(i) * x(i) * y(i) * y(i)
jf = jf + y(i) * y(i) * y(i) * y(i) * y(i)
REM -------------
bg = bg + x(i) * x(i) * y(i) * y(i)
cg = cg + y(i) * x(i) * y(i) * y(i)
dg = dg + x(i) * y(i) * x(i) * y(i) * y(i)
eg = eg + x(i) * x(i) * x(i) * y(i) * y(i)
fg = fg + y(i) * y(i) * x(i) * y(i) * y(i)
gg = gg + x(i) * y(i) * y(i) * x(i) * y(i) * y(i)
hg = hg + x(i) * x(i) * y(i) * x(i) * y(i) * y(i)
ig = ig + x(i) * x(i) * x(i) * x(i) * y(i) * y(i)
jg = jg + y(i) * y(i) * y(i) * x(i) * y(i) * y(i)
REM -------------
bh = bh + x(i) * x(i) * x(i) * y(i)
ch = ch + y(i) * x(i) * x(i) * y(i)
dh = dh + x(i) * y(i) * x(i) * x(i) * y(i)
eh = eh + x(i) * x(i) * x(i) * x(i) * y(i)
fh = fh + y(i) * y(i) * x(i) * x(i) * y(i)
gh = gh + x(i) * y(i) * y(i) * x(i) * x(i) * y(i)
hh = hh + x(i) * x(i) * y(i) * x(i) * x(i) * y(i)
ih = ih + x(i) * x(i) * x(i) * x(i) * x(i) * y(i)
jh = jh + y(i) * y(i) * y(i) * x(i) * x(i) * y(i)
REM -------------
bi = bi + x(i) * x(i) * x(i) * x(i)
ci = ci + y(i) * x(i) * x(i) * x(i)
di = di + x(i) * y(i) * x(i) * x(i) * x(i)
ei = ei + x(i) * x(i) * x(i) * x(i) * x(i)
fi = fi + y(i) * y(i) * x(i) * x(i) * x(i)
gi = gi + x(i) * y(i) * y(i) * x(i) * x(i) * x(i)
hi = hi + x(i) * x(i) * y(i) * x(i) * x(i) * x(i)
ii = ii + x(i) * x(i) * x(i) * x(i) * x(i) * x(i)
ji = ji + y(i) * y(i) * y(i) * x(i) * x(i) * x(i)
REM -------------
bj = bj + x(i) * y(i) * y(i) * y(i)
cj = cj + y(i) * y(i) * y(i) * y(i)
dj = dj + x(i) * y(i) * y(i) * y(i) * y(i)
ej = ej + x(i) * x(i) * y(i) * y(i) * y(i)
fj = fj + y(i) * y(i) * y(i) * y(i) * y(i)
gj = gj + x(i) * y(i) * y(i) * y(i) * y(i) * y(i)
hj = hj + x(i) * x(i) * y(i) * y(i) * y(i) * y(i)
ij = ij + x(i) * x(i) * x(i) * y(i) * y(i) * y(i)
jj = jj + y(i) * y(i) * y(i) * y(i) * y(i) * y(i)
REM ----------------
NEXT i
m=10  !EINTRAG ANZAHL GLEICHUNGEN
n=10  !EINTRAG ANZAHL UNBEKANNTE
p1(1,1)=anz  !anz=ANZAHL DER MESSUNGEN BZW. REFERENZPUNKTE
p1(1,2)=b1   !p(m,n)=m=ZEILEN,n+1=SPALTEN; n+1=UNBEKANNTE+RESIDUUM
p1(1,3) = c1
p1(1,4) = d1
p1(1,5) = e1
p1(1,6) = f1
p1(1,7) = g1
p1(1,8) = h1
p1(1,9) = i1
p1(1,10) = j1
p1(1,11) = rr
REM -------------------------
p1(2,1) = b1
p1(2,2) = bb
p1(2,3) = cb
p1(2,4) = db
p1(2,5) = eb
p1(2,6) = fb
p1(2,7) = gb1
p1(2,8) = hb
p1(2,9) = ib
p1(2,10) = jb
p1(2,11) = br
REM -------------------------
p1(3,1) = c1
p1(3,2) = bc
p1(3,3) = cc
p1(3,4) = dc
p1(3,5) = ecc
p1(3,6) = fc
p1(3,7) = gc
p1(3,8) = hc
p1(3,9) = ic
p1(3,10) = jc
p1(3,11) = cr
REM -------------------------
p1(4,1) = d1
p1(4,2) = bd
p1(4,3) = cd
p1(4,4) = dd
p1(4,5) = ed
p1(4,6) = fd
p1(4,7) = gd
p1(4,8) = hd
p1(4,9) = id
p1(4,10) = jd
p1(4,11) = dr
REM -------------------------
p1(5,1) = e1
p1(5,2) = be
p1(5,3) = ce
p1(5,4) = de
p1(5,5) = ee
p1(5,6) = fe
p1(5,7) = ge
p1(5,8) = he
p1(5,9) = ie
p1(5,10) = je
p1(5,11) = er
REM -------------------------
p1(6,1) = f1
p1(6,2) = bf
p1(6,3) = cf
p1(6,4) = df
p1(6,5) = ef
p1(6,6) = ff
p1(6,7) = gf
p1(6,8) = hf
p1(6,9) = if
p1(6,10) = jf
p1(6,11) = fr
REM -------------------------
p1(7,1) = g1
p1(7,2) = bg
p1(7,3) = cg
p1(7,4) = dg
p1(7,5) = eg
p1(7,6) = fg
p1(7,7) = gg
p1(7,8) = hg
p1(7,9) = ig
p1(7,10) = jg
p1(7,11) = gr
REM -------------------------
p1(8,1) = h1
p1(8,2) = bh
p1(8,3) = ch
p1(8,4) = dh
p1(8,5) = eh
p1(8,6) = fh
p1(8,7) = gh
p1(8,8) = hh
p1(8,9) = ih
p1(8,10) = jh
p1(8,11) = hr
REM -------------------------
p1(9,1) = i1
p1(9,2) = bi
p1(9,3) = ci
p1(9,4) = di
p1(9,5) = ei
p1(9,6) = fi
p1(9,7) = gi
p1(9,8) = hi
p1(9,9) = ii
p1(9,10) = ji
p1(9,11) = ir
REM -------------------------
p1(10,1) = j1
p1(10,2) = bj
p1(10,3) = cj
p1(10,4) = dj
p1(10,5) = ej
p1(10,6) = fj
p1(10,7) = gj
p1(10,8) = hj
p1(10,9) = ij
p1(10,10) = jj
p1(10,11) = jr
RETURN
PROCEDURE messung
  r = 5000 !EFFEKTIVE KAMERABRENNWEITE
  z = ATN(TAN(dek(i)) * COS(ar(i) - aro)) 
  x = -r * COS(z) * TAN(ar(i) - aro) / COS(z - do)  //x-Wert Paßpunktstern der Aufnahme
  y = r * TAN(z - do)      !y-Wert Paßpunktsterrn
RETURN
 

Alle Rechte vorbehalten (all rights reserved), auch die der fotomechanischen Wiedergabe und der Speicherung in elektronischen Medien, Translation usw. Dasselbe gilt für das Recht der öffentlichen Wiedergabe. Copyright © by H. Schumacher, Spaceglobe

Sternbeobachter - Sterntagebuch - Produktinformation - www.spaceglobe.de