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 3x Interpreter auf das Festplatten-Laufwerk >C<. Der Atari ST Interpreter akzeptiert lediglich die Einfügung von ASCII-Dateien mit LST-Extension unter MERGE (Dateiname.LST - z. B. mit WordPad, Editor o.a. erstellt).

REM ATARI ST GFA-BASIC48 PROGR. FIREBALL
DIM hr(11),hl(11),hb(11),el(11),eb(11),rr(11),x(11),nam$(10),s(50,2),c(50,2),i(3),ia(3),p(20,20),ko(20),k(50),k1(50),na$(4)
DIM v(20),w(20),y(20),xx(600),yy(600),ll2(600),ll1(600),up1(600),e$(600),ent1(600),h2$(180),h3$(90),h4$(90),c1(600),d(600),xx1(600),yy1(600),zz1(600),xax(60 0),yay(600),zaz(600),ent(600),ho(600),az(600),s$(600),rx(600),d1(200),a1(200),d2(200
),a2(200)
DEFFN z(x)=x-INT(x/(PI*2))*(PI*2)
DEFFN gr(g,m,s)=gg+gm/60+gs/3600
DEFFN gr1(g,m,s)=gg*i+gm*i/60+gs*i/3600
DEFFN sinh(x)=(EXP(x)-EXP(-x))/2
DEFFN cosh(x)=(EXP(x)+EXP(-x))/2
DEFFN tanh(x)=-EXP(-x)/(EXP(x)+EXP(-x))*2+1
DEFFN arctanh(x)=LOG((1+x)/(1-x))/2
DEFFN aco(x)=ATN(SQR(1-x*x)/x)
DEFFN ggb(x)=ATN((1+0.006739501819*(6378140/(6378140+eh)))*TAN(x))
DEFFN gzb(x)=ATN(1/(1+0.006739501819*(6378140/(6378140+eh)))*TAN(x))
DEFFN ro(x)=6356.755288158/SQR(1-0.006694384999592*COS(FN gzb(x))^2)+eh/1000
SETCOLOR 0,0
xo=320  !GRAPHIKAUFLÖSUNG 640x400
yo=200
ae=149597870 !=1 ASTRONOMISCHE EINHEIT 1 AE=149597870 KM
xx$="6"
CLS
BOX 1,1,619,399
DEFTEXT 1,0,0,32
TEXT 220,40,"FIREBALL"
PRINT AT(3,7);
REM INPUT "DATUM (Z. B. 1.2.1999)....................";a$
a$="28.7.1978"
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))
GOSUB jd
jd1=jd
tz=(jd1-2378496.5)/36525
ddt=46/3600 !EINGABE SEK. EPHEMERIDENZEITKORREKTUR (Z. B. ddt=42 SEK. FÜR 1972)
PRINT AT(3,9);
REM INPUT "UHRZEIT UT (Z. B. 12h3m1s)................";gg,gm,gs  !WELTZEIT UT
gg=22
gm=2
gs=29
ut=FN gr(g,m,s)
jd1=jd+(ut+ddt)/24
t=(jd1-2451545)/36525  !JULIAN. JAHRHUNDERTE AB EPOCHE J2000
te=(2451545-jd1)/36525
fo6=t/10 !JULIAN. JAHRTAUSENDE AB EPOCHE J2000
t3=(jd-2451545)/36525
f62=(jd1-2451545)/365.25
PRINT AT(3,11);
REM INPUT "MASSTAB (1 = 1:1, 0.1 = 1:10).............";ska
ska=0.2 !----------
PRINT AT(23,10);"Taste >1< - METEORBAHNBERECHNUNG"
PRINT AT(23,11);"Taste >2< - STERNKARTE"
PRINT AT(23,12);"Taste >3< - PLANETARIUM"
PRINT AT(23,13);"Taste >4< - NEUSTART"
PRINT AT(23,14);"Taste >S< - ENDE"
DO
  w$=UPPER$(INKEY$)
  IF w$="1" THEN
    GOSUB metbahn
  ENDIF
  IF w$="2" THEN
    GOSUB karte  !GNOMONISCHE STERNKARTE
  ENDIF
  IF w$="3" THEN
    GOSUB rium    !PLANETARIUM
    w$="4"
  ENDIF
  IF w$="4" THEN
    RUN
  ENDIF
  EXIT IF w$="S"
LOOP
END
PROCEDURE metbahn
  REM -------------------- METEORBAHN ----------------
  REM SUB-METEOR ORT (GROUNDTRACK) --------------------------
  REM  METHODE I
  eh=0 !EINTRAG SEEHÖHE STATION A (41 N.N.)
  b=50+38/60+0/3600  !EINTRAG GEOGRAPH. BREITE (+50o31'30.12'') STATION A
  b=RAD(b)
  b1=FN gzb(b) !GEOZENTRISCHE BREITE STATION A
  r1=FN ro(b)  !GEOZENTR. DISTANZ STATION A + SEEH™HE
  l1=+(6+47/60+0/3600) !EINTRAG GEOGRAPH. LÄNGE STATION A
  da=RAD(8) !DEKLINATION ANFANGSPUNKT DER METEORSPUR STATION A
  ara=RAD(292.5) !REKTASZENSION ANFANGSPUNKT STATION A
  de=RAD(12) !DEKLINATION ENDPUNKT DER METEORSPUR STATION A
  are=RAD(294) !REKTASZENSION ENDPUNKT STATION A
  REM ----------------------------------
  CLS
  eh=0 !EINTRAG SEEHÖHE STATION B
  bb1=50+48/60+0/3600  !EINTRAG GEOGRAPH. BREITE STATION B
  bb1=RAD(bb1)
  b2=FN gzb(bb1)
  r2=FN ro(bb1)
  l2=+(7+14/60+0/3600) !EINTRAG GEOGRAPH. LÄNGE STATION B
  dba=RAD(2.48171)   !DEKLINATION ANFANGSPUNKT DER METEORSPUR STATION B
  arba=RAD(278.64886) !REKTASZENSION ANFANGSPUNKT STATION B
  dbe=RAD(3.58468)   !DEKLINATION ENDPUNKT DER METEORSPUR STATION B
  arbe=RAD(273.99534) !REKTASZENSION ENDPUNKT STATION B
  REM ----------------------------------
  x1=r1*COS(b1)*COS(RAD(l1)) !KARTESISCHE STATIONSKOORDINATEN
  y1=r1*COS(b1)*SIN(RAD(l1))
  z1=r1*SIN(b1)
  x2=r2*COS(b2)*COS(RAD(l2))
  y2=r2*COS(b2)*SIN(RAD(l2))
  z2=r2*SIN(b2)
  dx=x2-x1
  dy=y2-y1
  dz=z2-z1
  REM ----- GREENWICH ORTSSTERNZEIT
  lgeo=0
  GOSUB sternz
  gsz=RAD(ru*15)
  REM SUB-METEOR ENDPUNKT ----------------------------
  dbve=50 !EVTL. VARIIEREN
  dao=de
  arao=are
  REPEAT
    x=dbve*COS(dbe)*COS(arbe-gsz)
    y=dbve*COS(dbe)*SIN(arbe-gsz)
    z=dbve*SIN(dbe)
    x3=dx+x
    y3=dy+y
    z3=dz+z
    dave=SQR(x3^2+y3^2+z3^2) !ENTFERNUNG ENDPUNKT METEORSPUR VON STATION A
    de1=ASIN(z3/dave)
    x=x3/(dave*COS(de1))
    y=y3/(dave*COS(de1))
    are1=FN z(ATN(y/(1+x))*2+gsz)
    ww1=DEG(ABS(de1-dao))
    ww=DEG(ABS(are1-arao))
    dbve=dbve+ww  !ENTFERNUNG ENDPUNKT DER METEORSPUR VON STATION B
    PRINT AT(1,1);USING "ITERATION: ####.############",ww1
    PRINT AT(30,1);USING "####.############",ww
    UNTIL ww<1.0E-10 !Schleife
    k(1)=ROUND(dave,3)
    k(2)=ROUND(dbve,3)
    REM ---------------
    gx=r1*COS(b1)*COS(RAD(l1))+dave*COS(de)*COS(are-gsz)
    gy=r1*COS(b1)*SIN(RAD(l1))+dave*COS(de)*SIN(are-gsz)
    gz=r1*SIN(b1)+dave*SIN(de)
    gd=SQR(gx^2+gy^2+gz^2) !geozentr. Vektor des Meteorpunkts md(da,ara)
    gdek=ASIN(gz/gd) !GEOZENTR. BREITE SUB-METEOR PUNKT (ENDPUNKT)
    x=gx/(gd*COS(gdek))
    y=gy/(gd*COS(gdek))
    gar=FN z(ATN(y/(1+x))*2) !GEOGRAPH. LÄNGE SUB-METEOR ENDPUNKT
    eh=0
    gb1=FN ggb(gdek) !GEOGRAPHISCHE BREITE SUB-METEOR ENDPUNKT
    gr1=FN ro(gb1)   !GEOZENTR. ENTFERNUNG SUB-METEOR ENDPUNKT
    hm=gd-gr1          !HÖHE SUB-METEOR ENDPUNKT šBER N.N.
    REM ---------------
    br=gdek !(=geozentrische Breite=Deklination)
    r=gd    !(=geozentrische Entfernung Endpunkt Meteor)
    FOR i=1 TO 10  STEP !ERDABLATTUNG
      gr=6356.755288158/SQR(1-0.006694384999592*COS(d1)^2)
      lat=ATN(1.00673950182*TAN(d1))
      h=SQR(r^2-gr^2*SIN(lat-d1)^2)-gr*COS(lat-d1)
      d1=br-ASIN((h/r)*SIN(lat-d1))
    NEXT i
    k(3)=ROUND(DEG(lat),5)
    k(4)=ROUND(DEG(d1),5)
    k(5)=ROUND(DEG(gar),5)
    k(6)=ROUND(gr,3)
    k(7)=ROUND(h,3)
    REM SUB-METEOR ANFANGSPUNKT ------------------------
    dbva=50 !EVTL. VARIIEREN
    dao=da
    arao=ara
    REPEAT
      x=dbva*COS(dba)*COS(arba-gsz)
      y=dbva*COS(dba)*SIN(arba-gsz)
      z=dbva*SIN(dba)
      x3=dx+x
      y3=dy+y
      z3=dz+z
      dava=SQR(x3^2+y3^2+z3^2) !ENTFERNUNG ANFANGSPUNKT VON STATION A
      da1=ASIN(z3/dava)
      x=x3/(dava*COS(da1))
      y=y3/(dava*COS(da1))
      ara1=FN z(ATN(y/(1+x))*2+gsz)
      ww1=DEG(ABS(da1-dao))
      ww=DEG(ABS(ara1-arao))
      dbva=dbva+ww  !ENTFERNUNG ANFANGSPUNKT DER METEORSPUR VON STATION B
      PRINT AT(1,1);USING "ITERATION: ####.############",ww1
      PRINT AT(30,1);USING "####.############",ww
      UNTIL ww<1.0E-10 !Schleife
      CLS
      k(8)=dava
      k(9)=dbva
      REM ---------------
      gx1=r1*COS(b1)*COS(RAD(l1))+dava*COS(da)*COS(ara-gsz)
      gy1=r1*COS(b1)*SIN(RAD(l1))+dava*COS(da)*SIN(ara-gsz)
      gz1=r1*SIN(b1)+dava*SIN(da)
      gd=SQR(gx1^2+gy1^2+gz1^2) !geozentr. Vektor des Meteorpunkts md(da,ara)
      gdek=ASIN(gz1/gd) !GEOZENTR. BREITE SUB-METEOR PUNKT (ANFANGSPUNKT)
      x=gx1/(gd*COS(gdek))
      y=gy1/(gd*COS(gdek))
      gar=FN z(ATN(y/(1+x))*2)
      eh=0
      gb1=FN ggb(gdek) !GEOGRAPHISCHE BREITE DES SUB-METEOR PUNKTES
      gr=FN ro(gb1)   !GEOZENTR. ENTFERNUNG DES SUB-METEOR PUNKTES
      hm=gd-gr         !HÖHE SUB-METEOR PUNKT ÜBER N.N.
      REM ---------------
      br=gdek !(=geozentrische Breite=Deklination)
      r=gd    !(=geozentrische Entfernung des Meteors)
      FOR i=1 TO 10 STEP !ERDABPLATTUNG
        gr=6356.755288158/SQR(1-0.006694384999592*COS(d1)^2)
        lat=ATN(1.00673950182*TAN(d1))
        h=SQR(r^2-gr^2*SIN(lat-d1)^2)-gr*COS(lat-d1)
        d1=br-ASIN((h/r)*SIN(lat-d1))
      NEXT i
      k(10)=ROUND(DEG(lat),5)
      k(11)=ROUND(DEG(d1),5)
      k(12)=ROUND(DEG(gar),5)
      k(13)=ROUND(gr,3)
      k(14)=ROUND(h,3)
      REM ------------------
      lx=gx1-gx
      ly=gy1-gy
      lz=gz1-gz
      l=SQR(lx^2+ly^2+lz^2) !WAHRE LÄNGE DER METEORSPUR IN DER ATMOSPHÄRE (km)
      k(15)=ROUND(l,4)
      REM SPHÄRISCHE KOORDINATEN DES SCHEINBAREN RADIANTEN (ANFANG- u. ENDPUNKT METEORSPUR)
      y=SIN(are)/TAN(de)-SIN(ara)/TAN(da)
      x=COS(are)/TAN(de)-COS(ara)/TAN(da)
      s=SQR(x^2+y^2)
      x=x/s
      y=y/s
      arp1=FN z(ATN(-(x/(1+y)))*2-PI)
      dp1=ATN(-COS(arp1-are)/TAN(de))
      REM ----------------------------
      y=SIN(arbe)/TAN(dbe)-SIN(arba)/TAN(dba)
      x=COS(arbe)/TAN(dbe)-COS(arba)/TAN(dba)
      s=SQR(x^2+y^2)
      x=x/s
      y=y/s
      arp2=FN z(ATN(-(x/(1+y)))*2-PI)
      dp2=ATN(-COS(arp2-arbe)/TAN(dbe))
      REM -------------------------------
      y=SIN(arp1)/TAN(dp1)-SIN(arp2)/TAN(dp2)
      x=COS(arp1)/TAN(dp1)-COS(arp2)/TAN(dp2)
      arr=ATN(1/-(y/x))  !REKTASZENSION RADIANT, ÄQUINOKTIUM DES DATUMS
      s=SQR(x^2+y^2)
      x=x/s
      y=y/s
      arr=FN z(ATN(-(x/(1+y)))*2-PI)
      dr=ATN(-COS(arr-arp1)/TAN(dp1)) !DEKLINATION RADIANT
      dro=dr
      arro=arr
      GOSUB pr !SCHEINBARER RADIANT ÄQUINOKTIUM J2000
      k(25)=ROUND(DEG(do1),5) !DEKLIN. RADIANT J2000
      k(26)=ROUND(DEG(aro1)/15,5) !AR RADIANT J2000
      lgeo=l1
      GOSUB sternz
      osz1=RAD(ru*15) !ORTSSTERNZEIT STATION A
      hr1=ASIN(SIN(b)*SIN(dr)+COS(b)*COS(dr)*COS(osz1-arr)) !HÖHE RADIANT
      IF hr1<0 THEN
        dr=-dr
        arr=FN z(arr+PI)
      ENDIF
      y=(COS(dr)*SIN(osz1-arr))/COS(hr1)
      x=(-COS(b)*SIN(dr)+SIN(b)*COS(dr)*COS(osz1-arr))/COS(hr1)
      azr=FN z(ATN(y/(1+x))*2) !AZIMUT RADIANT STATION A
      k(16)=ROUND(DEG(dr),5)
      k(17)=ROUND(DEG(arr/15),5)
      sw=ACOS(SIN(de)*SIN(da)+COS(de)*COS(da)*COS(are-ara))
      PRINT AT(1,1);
      PRINT "SUB-METEOR ORT (AUFLEUCHTPUNKT):"
      PRINT "GEOGR. BREITE: ";k(10);" GRAD"
      PRINT "GEOZ.  BREITE: ";k(11);" GRAD"
      PRINT "GEOGR.  LÄNGE: ";k(12);" GRAD"
      PRINT "GEOZ.   ENTF.: ";k(13);" km"
      PRINT "Höhe über N.N: ";k(14);" km"
      PRINT "SUB-METEOR ORT (VERLÖSCHPUNKT):"
      PRINT "GEOGR. BREITE: ";k(3);" GRAD"
      PRINT "GEOZ.  BREITE: ";k(4);" GRAD"
      PRINT "GEOGR.  LÄNGE: ";k(5);" GRAD"
      PRINT "GEOZ.   ENTF.: ";k(6);" km"
      PRINT "Höhe über N.N: ";k(7);" km"
      PRINT "DEKLIN. SCHEINB. RADIANT: ";k(16);" GRAD"
      PRINT "AR      SCHEINB. RADIANT: ";k(17);" Std."
      PRINT "WAHRE LÄNGE METEORSPUR..: ";k(15);" km"
      PRINT "LÄNGE METEORSPUR IN GRAD: ";ROUND(DEG(sw),5);" GRAD"
      PRINT "SCHEINB. RADIANT HÖHE...: ";ROUND(DEG(hr1),4);" GRAD"
      PRINT "SCHEINB. RADIANT AZIMUT.: ";ROUND(DEG(azr),4);" GRAD"
      PRINT AT(60,13);k(25);" GRAD J2000"
      PRINT AT(60,14);k(26);" STD. J2000"
      REM METHODE II ----------------------------
      REM SUB-METEOR ORT (VERLÖSCHPUNKT DER METEORBAHN) ------------------
      h1=ASIN(SIN(b)*SIN(de)+COS(b)*COS(de)*COS(osz1-are)) !HÖHE šBR HORIZONT ENDPUNKT METEOR
      y=(COS(de)*SIN(osz1-are))/COS(h1)
      x=(-COS(b)*SIN(de)+SIN(b)*COS(de)*COS(osz1-are))/COS(h1)
      az1=FN z(ATN(y/(1+x))*2) !AZIMUT ENDPUNKT METEOR STATION A
      REM -----------------------------------
      lgeo=l2
      GOSUB sternz
      osz2=RAD(ru*15) !ORTSSTERNZEIT STATION B
      h2=ASIN(SIN(bb1)*SIN(dbe)+COS(bb1)*COS(dbe)*COS(osz2-arbe)) !HÖHE ÜBER HORIZONT ENDPUNKT METEOR
      y=(COS(dbe)*SIN(osz2-arbe))/COS(h2)
      x=(-COS(bb1)*SIN(dbe)+SIN(bb1)*COS(dbe)*COS(osz2-arbe))/COS(h2)
      az2=FN z(ATN(y/(1+x))*2) !AZIMUT ENDPUNKT METEOR STATION B
      REM STATION A-----------------------------------
      x=COS(az1)
      y=SIN(az1)*SIN(b)
      j1=SQR(x^2+y^2)
      x=x/j1
      y=y/j1
      n1=FN z(RAD(l1)-ATN(y/(1+x))*2)
      j1=ASIN(j1)
      a11=SIN(n1)*SIN(j1)
      b11=-COS(n1)*SIN(j1)
      c11=-COS(j1)
      REM STATION B -------------
      x=COS(az2)
      y=SIN(az2)*SIN(bb1)
      j2=SQR(x^2+y^2)
      x=x/j2
      y=y/j2
      n2=FN z(RAD(l2)-ATN(y/(1+x))*2)
      j2=ASIN(j2)
      a12=SIN(n2)*SIN(j2)
      b12=-COS(n2)*SIN(j2)
      c12=-COS(j2)
      REM AKKUMULATION
      a=a11*a11+a12*a12
      ab=a11*b11+a12*b12
      bb=b11*b11+b12*b12
      ac=a11*c11+a12*c12
      bc=b11*c11+b12*c12
      REM NORMALGLEICHUNG ---------
      m=2 !ANZAHL GLEICHUNGEN
      n=2 !ANZAHL UNBEKANNTE
      p(1,1)=a
      p(1,2)=ab
      p(1,3)=ac
      p(2,1)=ab
      p(2,2)=bb
      p(2,3)=bc
      GOSUB elim
      xo1=ko(1)
      yo1=ko(2)
      s=SQR(xo1^2+yo1^2)
      y=yo1/s
      x=xo1/s
      eh=0
      lm=FN z(ATN(y/(1+x))*2) !GEOGRAPISCHE LÄNGE ENDPUNKT METEORSPUR
      bm=ATN(1/s)   !GEOGRAPHISCHE BREITE ENDPUNKT METEORSPUR
      rm=FN ro(bm) !GEOZENTR. ENTFERNUNG DES SUB-METEOR PUNKTES (ENDPUNKT)
      lm1=lm
      bm1=FN gzb(bm)
      rm1=rm
      lgeo=DEG(lm)
      GOSUB sternz
      oszm=RAD(ru*15)
      REM ---------------------
      k1(1)=ROUND(DEG(bm),5)
      k1(2)=ROUND(DEG(FN gzb(bm)),5) !GEOZENTR. BREITE
      k1(3)=ROUND(DEG(lm),5)
      k1(4)=ROUND(rm,3)
      REM u=ACOS(SIN(b)*SIN(bm)+COS(b)*COS(bm)*COS(RAD(l1)-lm))
      u=ASIN(SQR(SIN((b-bm)/2)^2+COS(b)*COS(bm)*SIN((RAD(l1)-lm)/2)^2))*2
      hme=2*6378.14*SIN(u/2)*(SIN(h1+u/2)/COS(h1+u))
      hme1=hme
      k1(5)=ROUND(hme,3) !HÖHE VERLÖSCHPUNKT ÜBER N.N.
      x=(6378.14+hme)*COS(bm)*COS(oszm)-6378.14*COS(b)*COS(osz1)
      y=(6378.14+hme)*COS(bm)*SIN(oszm)-6378.14*COS(b)*SIN(osz1)
      z=(6378.14+hme)*SIN(bm)-6378.14*SIN(b)
      dve=SQR(x*x+y*y+z*z)  !ENTFERNUNG DES VERLÖSCHPUNKTES VON STATION A
      dee=ASIN(z/dve) !DEKLIN. VERLÖSCHPUNKT STATION A
      x=x/(dve*COS(dee))
      y=y/(dve*COS(dee))
      aree=FN z(ATN(y/(1+x))*2) !AR VERLÖSCHPUNKT STATION A
      REM SUB-METEOR ORT (AUFLEUCHTPUNKT DER METEORBAHN) ------------------
      lgeo=l1
      GOSUB sternz
      osz=RAD(ru*15) !ORTSSTERNZEIT STATION A
      h1=ASIN(SIN(b)*SIN(da)+COS(b)*COS(da)*COS(osz-ara)) !HÖHE ÜBR HORIZONT ENDPUNKT METEOR
      y=(COS(da)*SIN(osz-ara))/COS(h1)
      x=(-COS(b)*SIN(da)+SIN(b)*COS(da)*COS(osz-ara))/COS(h1)
      az1=FN z(ATN(y/(1+x))*2) !AZIMUT ENDPUNKT METEOR STATION A
      REM -----------------------------------
      lgeo=l2
      GOSUB sternz
      osz=RAD(ru*15) !ORTSSTERNZEIT STATION B
      h2=ASIN(SIN(bb1)*SIN(dba)+COS(bb1)*COS(dba)*COS(osz-arba)) !HÖHE ÜBER HORIZONT ENDPUNKT METEOR
      y=(COS(dba)*SIN(osz-arba))/COS(h2)
      x=(-COS(bb1)*SIN(dba)+SIN(bb1)*COS(dba)*COS(osz-arba))/COS(h2)
      az2=FN z(ATN(y/(1+x))*2) !AZIMUT ENDPUNKT METEOR STATION B
      REM STATION A-----------------------------------
      x=COS(az1)
      y=SIN(az1)*SIN(b)
      j1=SQR(x^2+y^2)
      x=x/j1
      y=y/j1
      n1=FN z(ATN(y/(1+x))*2-RAD(l1))
      j1=ASIN(j1)
      a11=SIN(n1)*SIN(j1)
      b11=-COS(n1)*SIN(j1)
      c11=-COS(j1)
      REM STATION B -------------
      x=COS(az2)
      y=SIN(az2)*SIN(bb1)
      j2=SQR(x^2+y^2)
      x=x/j2
      y=y/j2
      n2=FN z(ATN(y/(1+x))*2-RAD(l2))
      j2=ASIN(j2)
      a12=SIN(n2)*SIN(j2)
      b12=-COS(n2)*SIN(j2)
      c12=-COS(j2)
      REM AKKUMULATION
      a=a11*a11+a12*a12
      ab=a11*b11+a12*b12
      bb=b11*b11+b12*b12
      ac=a11*c11+a12*c12
      bc=b11*c11+b12*c12
      REM NORMALGLEICHUNG ---------
      m=2 !ANZAHL GLEICHUNGEN  !LEAST SQUARE FIT
      n=2 !ANZAHL UNBEKANNTE
      p(1,1)=a
      p(1,2)=ab
      p(1,3)=ac
      p(2,1)=ab
      p(2,2)=bb
      p(2,3)=bc
      GOSUB elim
      xo1=ko(1)
      yo1=ko(2)
      s=SQR(xo1^2+yo1^2)
      x=xo1/s
      y=yo1/s
      he=0
      lm=FN z(PI*2-ATN(y/(1+x))*2) !GEOGRAPISCHE LÄNGE ENDPUNKT METEORSPUR
      bm=ATN(1/s)   !GEOGRAPHISCHE BREITE ENDPUNKT METEORSPUR
      rm=FN ro(bm) !GEOZENTR. ENTFERNUNG DES SUB-METEOR PUNKTES (ENDPUNKT)
      k1(6)=ROUND(DEG(bm),5)
      k1(7)=ROUND(DEG(FN gzb(bm)),5) !GEOZENTR. BREITE
      k1(8)=ROUND(DEG(lm),5)
      k1(9)=ROUND(rm,3)
      REM ----------------------
      REM u=ACOS(SIN(b)*SIN(bm)+COS(b)*COS(bm)*COS(RAD(l1)-lm))
      u=ASIN(SQR(SIN((b-bm)/2)^2+COS(b)*COS(bm)*SIN((RAD(l1)-lm)/2)^2))*2
      hme=2*6378.14*SIN(u/2)*(SIN(h1+u/2)/COS(h1+u))
      k1(10)=ROUND(hme,3)
      REM -----------------------
      lgeo=DEG(lm)
      GOSUB sternz
      oszm=RAD(ru*15)
      z=(6378.14+hme)*SIN(bm)-6378.14*SIN(b)
      x=(6378.14+hme)*COS(bm)*COS(oszm)-6378.14*COS(b)*COS(osz1)
      y=(6378.14+hme)*COS(bm)*SIN(oszm)-6378.14*COS(b)*SIN(osz1)
      dva=SQR(x*x+y*y+z*z)  !ENTFERNUNG DES VERLÖSCHPUNKTES VON STATION A
      daa=ASIN(z/dva) !DEKLIN. VERLÖSCHPUNKT STATION A
      x=x/(dva*COS(daa))
      y=y/(dva*COS(daa))
      araa=FN z(ATN(y/(1+x))*2) !AR VERLÖSCHPUNKT STATION A
      REM -----------------------
      REM  SPHÄRISCHE KOORDINATEN DES SCHEINBAREN RADIANTEN
      REM STATION A
      y=TAN(da)
      x=(TAN(de)-TAN(da)*COS(are-ara))/SIN(are-ara)
      g=SQR(x^2+y^2)
      x=x/g
      y=y/g
      q1=FN z(ara-ATN(y/(1+x))*2)
      g1=ATN(g)
      REM STATION B
      y=TAN(dba)
      x=(TAN(dbe)-TAN(dba)*COS(arbe-arba))/SIN(arbe-arba)
      g=SQR(x^2+y^2)
      x=x/g
      y=y/g
      q2=FN z(arba-ATN(y/(1+x))*2)
      g2=ATN(g)
      a11=SIN(q1)*SIN(g1)
      b11=-COS(q1)*SIN(g1)
      c11=-COS(g1)
      a12=SIN(q2)*SIN(g2)
      b12=-COS(q2)*SIN(g2)
      c12=-COS(g2)
      REM AKKUMULATION
      a=a11*a11+a12*a12
      ab=a11*b11+a12*b12
      bb=b11*b11+b12*b12
      ac=a11*c11+a12*c12
      bc=b11*c11+b12*c12
      REM NORMALGLEICHUNG ---------
      m=2 !ANZAHL GLEICHUNGEN
      n=2 !ANZAHL UNBEKANNTE
      p(1,1)=a
      p(1,2)=ab
      p(1,3)=ac
      p(2,1)=ab
      p(2,2)=bb
      p(2,3)=bc
      GOSUB elim
      xo1=ko(1)
      yo1=ko(2)
      s=SQR(xo1^2+yo1^2)
      x=xo1/s
      y=yo1/s
      arr=FN z(ATN(y/(1+x))*2)
      dr=ATN(1/s)
      hr=ASIN(SIN(b)*SIN(dr)+COS(b)*COS(dr)*COS(osz1-arr)) !HÖHE RADIANT STATION A
      IF hr<0 THEN
        dr=-dr
        arr=FN z(arr+PI)
      ENDIF
      k1(11)=ROUND(DEG(dr),5) !DEKLIN. SCHEINBARER RADIANT
      k1(12)=ROUND(DEG(arr/15),5) !AR SCHEINB. RADIANT
      REM --------------------------------------
      sw=ACOS(SIN(de)*SIN(da)+COS(de)*COS(da)*COS(are-ara))
      k1(13)=ROUND(DEG(sw),5) !SCHEINB. LÄNGE DER METEORSPUR IN GRAD
      dmw=ACOS(-SIN(dr)*SIN(de)-COS(dr)*COS(de)*COS(arr-are))
      l1=(dave*SIN(sw))/SIN(sw+dmw)
      k1(14)=ROUND(l1,4)
      ao=6378.14^2
      a1=2*6378.14*SIN(h1)
      y=(6378.14+hme)^2
      dmev=-a1/2+SQR(a1^2/4-(ao-y)) !ENTFERNUNG STATION A vom ENDPUNKT BAHNSPUR
      dmav=(dave*SIN(dmw))/SIN(sw+dmw) !ENTFERNUNG STATION A vom ANFANGSPUNKT BAHNSPUR
      PRINT AT(44,1);"METHODE II:"
      PRINT AT(44,2);k1(6);" GRAD"
      PRINT AT(44,3);k1(7);" GRAD"
      PRINT AT(44,4);k1(8);" GRAD"
      PRINT AT(44,5);k1(9);" km"
      PRINT AT(44,6);k1(10);" km"
      PRINT AT(44,8);k1(1);" GRAD"
      PRINT AT(44,9);k1(2);" GRAD"
      PRINT AT(44,10);k1(3);" GRAD"
      PRINT AT(44,11);k1(4);" km"
      PRINT AT(44,12);k1(5);" km"
      PRINT AT(44,13);k1(11);" GRAD"
      PRINT AT(44,14);k1(12);" STD."
      PRINT AT(44,15);k1(14);" km"
      PRINT AT(44,16);k1(13);" GRAD"
      REM bma,lma,hma,=geograph. Breite,Länge und Höhe über N.N. des Aufleuchtpunktes
      hre=ASIN(SIN(bm1)*SIN(dr)+COS(bm1)*COS(dr)*COS(oszm-arr))
      y=(COS(dr)*SIN(oszm-arr))/COS(hre)
      x=(-COS(bm1)*SIN(dr)+SIN(bm1)*COS(dr)*COS(oszm-arr))/COS(hre)
      azre=FN z(ATN(y/(1+x))*2)
      hee=rm1+hme1
      i=ATN((l1*COS(hre))/(hee+l1*SIN(hre)))
      o=(l1*SIN(hre+i/2))/COS(i/2)
      hma=hme1+o
      bma=ASIN(COS(i)*SIN(bm1)-SIN(i)*COS(bm1)*COS(azre))
      x=(COS(i)*COS(bm1)+SIN(i)*SIN(bm1)*COS(azre))/COS(bma)
      y=(SIN(i)*SIN(azre))/COS(bma)
      lma=FN z(lm1-ATN(y/(1+x))*2)
      REM KORREKTUR DES SCHEINBAREN RADIANTEN ---------------
      zd=ACOS(SIN(b)*SIN(dr)+COS(b)*COS(dr)*COS(osz1-arr)) !ZENITDISTANZ RADIANT STATION A
      y=(COS(b)*SIN(osz1-arr))/SIN(zd) !osz1,b=Ortsstenzeit u. geograph. Br. Station A
      x=(SIN(b)*COS(dr)-COS(b)*SIN(dr)*COS(osz1-arr))/SIN(zd)
      c=FN z(ATN(y/(1+x))*2) !PARALLAKTISCHER WINKEL RADIANT STATION A
      PRINT AT(1,21);
      INPUT "Geschwindigkeit (km/s) für Mitte der Meteorbahn ";vo
      INPUT "Verlangsamung   (km/s) für Mitte der Meteorbahn ";dv
      PRINT AT(1,21);SPC(60);!60;LEERZEICHEN
      PRINT AT(1,22);SPC(60)
      v=vo*(1-(8.77193*dv)/(vo^2*COS(zd))) !METORGEXWINDIGKEIT AUSSERHALB DER ATMOSPHÄRE
      o=FN z(osz1+PI/2)
      w=0.4561*COS(b)
      p=ACOS(COS(dr)*COS(arr-o))
      vgo=SQR(v^2+w^2-2*v*w*COS(p)) !BAHNGESCHWIDIGKEIT FšR NICHTROTIERENDE ERDE
      vg=SQR(vgo^2-11.16944^2) !GEOZENTRISCHE BAHNGECHWINDIGKEIT
      dz=ATN(((vgo-vg)/(vgo+vg))*TAN(zd/2))*2 !ZENITATTRAKTION RADIANT
      dar=-dz*SIN(c)*(1/COS(dr))
      dd=-dz*COS(c)
      ar1=FN z(arr+dar) !AR RADIANT FšR ZENITTATTRAKTION KORRIGIERT
      dr1=dr+dd !DEKLIN. RADIANT FšR ZENITATTRATKTION KORRIGIERT
      dar=-(0.46383/vgo)*COS(b)*COS(osz1-arr)*(1/COS(dr))
      dd=-(0.46383/vgo)*COS(b)*SIN(osz1-arr)*SIN(dr)
      ar2=FN z(ar1+dar) !AR RADIANT VON TÄGL. ABERRATION BERFREIT
      dr2=dr1+dd        !DEKLIN. RADIANT VON TÄGL. ABERRATION BEFREIT
      k(15)=ROUND(DEG(dr2),4)
      k(16)=ROUND(DEG(ar2)/15,4)
      lms=FN z(4.8951+628.33197*t) !mittl. ekl.Länge der Sonne
      ms=FN z(6.24006+628.301955*t) !mittl. Anomalie
      cc=0.03342*SIN(ms)+0.00035*SIN(2*ms) !Mittelpunktsgleichung
      ls=FN z(lms+cc) !genäherte ekl. ä„nge der Sonne (0.01 Grad)
      e=0.016708617 !Mittl. Erdbahnexentrizität
      p=FN z(1.7966+0.0125582*t) !Mittl. Länge Perihelion Erde
      tex=RAD((e/0.000290888)/60)*SIN(ls-p)
      al=FN z(ls-PI/2+tex) !APEX DER ERDBEWEGUNG
      wa=FN z(ms+cc) !wahre Anomalie
      rs = 10 ^ (((3040 - 15 * t) + (-727600 + 1810 * t) * COS(ms) + (-9090 + 50 * t) * COS(2 * ms)) / 100000000) !Entf. Erde-Sonne in AE
      bq=ASIN(SIN(dr2)*COS(ec)-COS(dr)*SIN(ec)*SIN(ar2)) !EKL. BR. RADIANT ar2
      x=(COS(dr2)*COS(ar2))/COS(bq)
      y=(SIN(dr2)*SIN(ec)+COS(dr2)*COS(ec)*SIN(ar2))/COS(bq)
      lq=FN z(ATN(y/(1+x))*2) !EKL. L¥NGE RADIANT dr2
      k(17)=ROUND(DEG(bq),4)
      k(18)=ROUND(DEG(lq),4)
      dro=dr2
      arro=ar2
      loo=ls
      GOSUB pr !KORRIGIERT RADIANT ÄQUINOKTIUM J2000
      k(27)=ROUND(DEG(do1),5) !DEKLIN. RADIANT J2000 KORRIGIERT
      k(28)=ROUND(DEG(aro1)/15,5) !AR RADIANT J2000 KORRIGIERT
      k(29)=ROUND(DEG(loo),5) !EKL. LÄNGE DER SONNE J2000
      vt=SQR(29.7847^2*(2/rs-1)) !BAHNGESCHWINDIGKEIT ERDE
      vh=SQR(vg^2+vt^2-2*vg*vt*COS(bq)*COS(lq-al)) !HELOZENTR. BAHNGESCHWINDIGKEIT METEOR
      z=ACOS(COS(bq)*COS(lq-al))
      y=SIN(bq)/SIN(z)
      x=(COS(bq)*SIN(lq-al))/SIN(z)
      s=FN z(ATN(y/(1+x))*2)
      mm=ASIN((vt/vh)*SIN(z))+z
      br=ASIN(SIN(mm)*SIN(s)) !EKLIPTIKALE BREITE DES WAHREN RADIANTEN, ÄQUINOKTIUM DES DATUMS
      x=COS(mm)/COS(br)
      y=(SIN(mm)*COS(s))/COS(br)
      lr=FN z(ATN(y/(1+x))*2+al) !EKL. LÄNGE DES WAHREN RADIANTEN, ÄQUINOKTIUM DES DATUMS
      k(19)=ROUND(DEG(br),4)
      k(20)=ROUND(DEG(lr),4)
      k(21)=ROUND(vt,3)
      k(22)=ROUND(vg,3)
      k(23)=ROUND(vh,3)
      k(24)=ROUND(DEG(ls),4)
     PRINT AT(1,19);"DEKLIN. RADIANT KORRIGIERT: ";k(15);" GRAD"
PRINT AT(1,20);"AR      RADIANT KORRIGIERT: ";k(16);" Std."
PRINT AT(1,21);"AR      RADIANT KORRIGIERT: ";ROUND(DEG(ar2),4);" GRAD"
PRINT AT(1,22);"EKL.BR. RADIANT KORRIGIERT: ";k(17);" GRAD"
PRINT AT(1,23);"EKL. L. RADIANT KORRIGIERT: ";k(18);" GRAD"
PRINT AT(1,24);"EKL. LÄNGE DER SONNE......: ";k(24);" GRAD"
PRINT AT(45,19);k(27);" GRAD J2000"
PRINT AT(45,20);k(28);" GRAD J2000"
PRINT AT(45,24);k(29);" GRAD J2000"
PRINT AT(70,24);"LEERTASTE"
REPEAT
UNTIL INKEY$=" "
REM BAHNELEMENTE ------
IF br>0 THEN
  kn=ls !EKL. LÄNGE AUFSTEIGENDER BAHNKNOTEN EKLIPTIK-METEORBAHN
  y=SIN(br)/SIN(u)
ELSE
  kn=FN z(ls+PI)
ENDIF
u=ACOS(COS(br)*COS(lr-ls))
y=SGN(br)*SIN(br)/SIN(u)
x=(COS(br)*SIN(lr-ls))/SIN(u)
i=FN z(ATN(y/(1+x))*2) !INCLINATION EKLITIK-METEORBAHN
k2=29.784691826
k1=0.01720209895
ko=2*k2^2/rs
IF vh^2<ko THEN
  na$(1)="ELLIPTISCHE METEORBAHN"
  nr=1
  a=1/(2/rs-vh^2/k2^2) !GROSSE HALBACHSE DER BAHN IN AE
  eo=ACOS(SGN(br)*(((rs*SIN(br))/(SQR(a)*SIN(i)))*(vh/k2)))
  p=a*COS(eo)^2 !BAHNPARAMETER
  e=SIN(eo) !EXENTRIZITÄT
  q=a*(1-e) !PERIHELDISTANZ
  x=p/rs-1
  y=(p/rs)*(1/TAN(u))
  e1=SQR(x^2+y^2)
  x=x/e1
  y=y/e1
  wa=FN z(ATN(y/(1+x))*2)
  wp=FN z(PI+ls-wa) !LÄNGE DES PERIHELS
  w=FN z(wp-kn) !PERIHELARGUMENT
  um=a^1.5 !UMLAUFZEIT SIDERISCHE JAHRE
  n=((PI*2)/365.242189)*(1/um)
  ex=ATN(TAN(wa/2)/SQR((1+e)/(1-e)))*2
  m1=ex-e*SIN(ex)
  m1=FN z(m1)
  tp=jd1-m1/n !JULIAN. DATUM PERIHELZEIT
  jd=tp
ENDIF
IF ABS(vh^2-ko)<0.01 THEN
  na$(2)="PARABOLISCHE METEORBAHN"
  nr=2
  q=rs*SIN(u)^2 !HELIOZENTR. PERIHELDISTANZ IN AE
  p=q/0.5
  wp=FN z(ls+2*u) !LÄNGE PERIHELION
  w=FN z(wp-kn) !PERIHELARGUMENT
  wa=ATN((1/TAN(lr-ls))*COS(i))*2 !WAHRE ANOMALIE
  tp=jd1-((1/3)*TAN(wa/2)^3+TAN(wa/2))*(SQR(2)/k1)*q^(3/2)
  n=0
  um=0
  jd=tp !PERIHELZEIT
ENDIF
IF vh^2>ko THEN
  na$(3)="HYPERBOLISCHE METEORBAHN"
  nr=3
  a=1/(2/rs-vh^2/k2^2)
  p=((rs*vh*COS(br)*SIN(lr-ls))/(k*COS(i)))^2 !BAHNPARAMETER
  e1=SQR(1+p/ABS(a))
  x=p/rs-1
  y=(p/rs)*(1/TAN(u))
  e=SQR(x^2+y^2) !EXENTRIZTITÄT
  x=x/e
  y=y/e
  wa=FN z(ATN(y/(1+x))*2) !WARE ANOMALIE
  q=a*(e-1)
  wp=FN z(PI+ls-wa) !LÄNGE PERIHELION
  w=FN z(wp-kn) !PERIHELARGUMENT
  um=0
  n=SQR(k1^2/ABS(a)^3) !MITTL. TÄGL. BEWEGUNG (RAD)
  h=FN arctanh(SQR((e-1)/(e+1))*TAN(wa/2))*2
  tp=jd1-SQR(ABS(a)^2/k1^2)*(e*FN sinh(h)-h)
  jd=tp
ENDIF
oa=a
oe=e
oi=i
ok=kn
ow=w
CLS
PRINT AT(1,1);"GESCHWINDIGKEIT ERDE......: ";k(21);" km/s"
PRINT AT(1,2);"GEOZENTR.   GESCHW. METEOR: ";k(22);" km/s"
PRINT AT(1,3);"HELIOZENTR. GESCHW. METEOR: ";k(23);" km/s"
IF nr=1 THEN
  PRINT AT(16,4);na$(nr)
ENDIF
IF nr=2 THEN
  PRINT AT(16,4);na$(nr)
ENDIF
IF nr=3 THEN
  PRINT AT(16,4);na$(nr)
ENDIF
PRINT AT(1,4);"BAHNELEMENTE:"
PRINT AT(1,5);"PERIHELZEIT (JD): ";tp;" / ";"TAG: ";FIX(tag);".";mon;".";jah;"  ";ROUND(FRAC(tag)*24,3);" STD."
PRINT AT(1,6);"UMLAUFZEIT......: ";ROUND(um,5);" SID. JAHRE"
PRINT AT(1,7);"Mittl. Bewegung  n: ";ROUND(DEG(n),6);" GRAD/TAG"
PRINT AT(1,8);"GROSSE HALBACHSE a: ";ROUND(a,5);" AE"
PRINT AT(1,9);"PERIHELDISTANZ   q: ";ROUND(q,4);" AE"
PRINT AT(1,10);"BAHNPARAMETER   p: ";ROUND(p,4);" AE"
PRINT AT(1,11);"EXENTRIZITÄT    e: ";ROUND(e,4)
PRINT AT(1,12);"PERIHELION.....wp: ";ROUND(DEG(wp),4);" GRAD"
PRINT AT(1,13);"Inclination     i: ";ROUND(DEG(i),4);" GRAD"
PRINT AT(1,14);"Perihelargument w: ";ROUND(DEG(w),4);" GRAD"
PRINT AT(1,15);"Aufst. Konten   k: ";ROUND(DEG(kn),4);" GRAD"
REM ---------------------
PRINT AT(1,17);"PLANETARIUM BAHNELEMENTE INTERN šBERNEHMEN - TASTE 5"
PRINT AT(1,18);"PLANETARIUM BAHNELEMENTE EXTERN EINGEBEN   - TASTE 6"
DO
  xx$=UPPER$(INKEY$)
  IF xx$="5" OR xx$="6" THEN
    GOTO jump
  ENDIF
LOOP
jump:
REM --------------
PRINT AT(23,21);"Taste >1< - METEORBAHNBERECHNUNG"
PRINT AT(23,22);"Taste >2< - STERNKARTE"
PRINT AT(23,23);"Taste >3< - PLANETARIUM"
PRINT AT(23,24);"Taste >4< - NEUSTART"
RETURN
PROCEDURE rium
  CLS
  BOX 1,1,619,399
  DEFTEXT 1,0,0,32
  TEXT 200,40,"PLANETARIUM"
  REM EKLIPTIKALE KOORDINATEN DES BEOBACHTERS --------------------
  PRINT AT(3,7);
  INPUT "ENTFERNUNG ZUR SONNE (ASTRONOM. EINHEITEN)";d
  PRINT AT(3,8);
  INPUT "EKLIPTIKALE BREITE DES BEOBACHTERS (GRAD) ";br
  PRINT AT(3,9);
  INPUT "EKLIPTIKALE L¥NGE  DES BEOBACHTERS (GRAD) ";le
  CLS
  br=RAD(br)+1.0E-07  !EKLIPTIKALE BREITE DES BEOBACHTERS
  le=RAD(le)+1.0E-07   !EKL. LÄNGE DES BETRACHTERS
  rb=d*ae
  fp=DEG(1)*32
  rd=DEG(1/rb)*32
  r1=rd*696000*ska  !SONNENRADIUS KM IN PIX
  IF r1<4 THEN
    DEFFILL 1,2,8
    PCIRCLE xo,yo,3
  ENDIF
  REM ---- RICHTUNG ARIES-PUNKT (RAM HORN) NULL-LÄNGENGRAD ---------
  ra=3*ae
  e=0
  i=0
  k=0
  w=0
  nr=6.29
  GOSUB pb
  LINE 320,200,x,y
  REM -------------------
  GOSUB kr
  REM MITTL. BAHNELEMENTE PLANETEN, MITTL. ÄQUINOKTIUM DES DATUMS
  RESTORE el
  nr=0.1 !EINTRAG SCHRITTWEITE
  FOR pl=1 TO 8
    READ e$(pl),i,k,w,e,ra
    PRINT AT(1,1);SPC(20)
    PRINT AT(1,1);e$(pl)
    i=RAD(i)
    k=RAD(k)
    w=RAD(w)
    ra=ra*ae
    GOSUB pb
    EXIT IF rb<3 AND pl=4
    EXIT IF rb<9 AND pl=8 AND ABS(br)>0.35
  NEXT pl
  n:
  REM BAHNELEMENTE SPORADISCHE METEORE/PERIODISCHE METEORSCHWÄRME ---------
  PRINT AT(1,1);SPC(20)
  PRINT AT(1,1);"OBJEKTBAHN"
  REM EXTERNE EINGABE BAHNELEMENTE (KOMET,ASTEROID USW.)
  IF xx$="6" THEN
    epo=2000      !EINTRAG JAHR ÄUINOKTIALEPOCHE BAHNELEMENTE
    jd3=2415020.31352+365.242198781*(epo-1900)
    ra=2.7883*ae      !EINTRAG GROSSE BAHNHALBACHSE (a) in AE
    e=0.7166        !EINTRAG EXENTRIZITÄT (PARABEL e=1)
    q=0.76*ae     !EINTRAG PERIHELDISTANZ BEI PARABOLISCHER BAHN
    ii=RAD(8.7391)   !EINTRAG INCLINATION EPOCHE
    kk=RAD(125.4597)   !EINTRAG AUFST. KNOTENLÄNGE EPOCHE
    ww=RAD(242.04)   !EINTRAG PERIHELARGUMENT EPOCHE
    i=ii
    k=kk
    w=ww
    GOSUB praez  !BAHNELEMENTE, ÄQUINOKTIUM DES DATUMS
    a2=10.565    !EINTRAG TAG+BRUCHTEIL PERIHELDURCHGANG
    a3=1     !EINTRAG MONAT PERIHELDURCHGANG
    a4=1974    !EINTRAG JAHR  PERIHELDURCHGANG
    GOSUB jd
    jd2=jd !JULIAN. DATUM DER PERIHELDURCHGANGSZEIT
  ENDIF
  REM INTERN BERECHNETE BAHNELEMENTE
  IF xx$="5" THEN
    i=oi  !BAHNELEMENTE, ÄQUINOKTIUM DES DATUMS
    k=ok
    w=ow
    e=oe
    ra=oa*ae
    jd2=tp
  ENDIF
  nr=0.05    !EINTRAG SCHRITTWEITE ---------------------------
  GOSUB pb   !BAHN, ¥QUINOKTIUM DES DATUMS
  gau=0.01720209895
  IF e<1 THEN !ELLIPTISCHE BAHN
  n=(gau/(ra/ae)^(3/2)) !Mittl. tägl. Bewegung in RAD
  ma=n*(jd1-jd2) !Mittl. Anomalie
ENDIF
IF e=1 THEN !PARABOLISCHE BAHN
aa1=(3*gau*(jd1-jd2))/(2*SQR(2)*(q/ae)^(3/2))
ENDIF
IF e>1 THEN !HYPERBOLISCHE BAHN
n=SQR(gau^2/ABS(ra/ae)^3) !MITTL. TÄGL. BEWEGUNG (RAD)
ma=n*(jd1-jd2) !MITTL. ANOMALIE
ENDIF
GOSUB pb1
nam$(1)="Objekt"
el(1)=el
eb(1)=eb
rr(1)=rr
GOSUB pb2
x(1)=INT(x)
y(1)=INT(y)
REM ------------------------------------
'  GOSUB solu
GOSUB exit
DO
SHOWM !MAUSZEIGER EIN - BAHNPUNKT MIT MAUSZEIGER ANKLICKEN
w$=UPPER$(INKEY$)
IF w$="S" THEN
RUN  "!TASTE >S< NEUSTART"
ENDIF
IF w$="P" THEN
BSAVE "f1",XBIOS(2),32000 !TASTE >P< GRAPHIK AUF DISKETTE SPEICHERN (ATARI ST)
ENDIF
x=MOUSEX
y=MOUSEY
z=MOUSEK
SHOWM
FOR is=1 TO 10
IF INT(x)=x(is) AND INT(y)=y(is) THEN
  GOSUB info
  SETMOUSE x(is)+1,y(is)+1 !MAUSZEIGER versetzen
ENDIF
NEXT is
IF w$="R" THEN
GOTO ju
ENDIF
LOOP
ju:
RETURN
PROCEDURE info
  PRINT AT(1,1);SPC(30)
  PRINT AT(1,2);SPC(30)
  PRINT AT(1,3);SPC(30)
  PRINT AT(1,4);SPC(30)
  PRINT AT(1,1);nam$(is)
  PRINT AT(1,2);"Heliozentr. Breite: ";ROUND(DEG(eb(is)),5)
  PRINT AT(1,3);"Heliozentr. L„nge: ";ROUND(DEG(el(is)),5)
  PRINT AT(1,4);"Heliozentr. Entf.: ";ROUND(rr(is)/ae,7);" AE"
RETURN
PROCEDURE kr  !RADIUS SONNE - ZENTRALKÖRPER DES SONNENSYSTEMS
ns=1
FOR i1=0 TO PI*2 STEP 0.01
  x=xo+r1*COS(i1)
  y=yo+r1*SIN(i1)
  IF ns=1 THEN
    ns=0
    PLOT x,y
  ENDIF
  DRAW  TO x,y
NEXT i1
RETURN
PROCEDURE pb  !ORBIT
ns=1
IF e<1 THEN
nn=-nr
ENDIF
IF e=1 THEN
nn=-50  !EINTRAG BAHNLÄNGE PARABEL
jj=ABS(nn)
ENDIF
IF e>1 THEN
nn=-8  !EINTRAG BAHNLÄNGE HYPERBEL
jj=ABS(nn)
ENDIF
DO
ADD nn,nr
REM TRANSFORMATION BAHNEBENE ZU EKLIPTIKEBENE
IF e<1 THEN
  ra1=ra*(1-e*COS(nn))
  v=ATN(SQR((1+e)/(1-e))*TAN(nn/2))*2
ENDIF
IF e=1 THEN
  v=ATN((SQR(nn^2+1)+nn)^(1/3)-(SQR(nn^2+1)-nn)^(1/3))*2
  ra1=(2*q)/(1+COS(v))
ENDIF
IF e>1 THEN !HYPERBOLIDE BAHN
v=ATN(SQR((e+1)/(e-1))*FN tanh(nn/2))*2
ra1=ABS(ra)*(e*FN cosh(nn)-1)
ENDIF
x=ra1*(COS(k)*COS(w+v)-SIN(k)*SIN(w+v)*COS(i))
y=ra1*(SIN(k)*COS(w+v)+COS(k)*SIN(w+v)*COS(i))
zv=ra1*(SIN(w+v)*SIN(i))
rr=SQR(x^2+y^2+zv^2)  !HELIOZENTR. ENTFERNUNG
eb=ASIN(zv/rr)  !HELIOZENTR. BREITE
o=COS(eb)
IF o==0 THEN
o=1.0E-20
ENDIF
x2=x/(rr*o)
y2=y/(rr*o)
IF x2<=-1 THEN
el=PI
ELSE
el=FN z(ATN(y2/(1+x2))*2) !HELIOZENTR. LÄNGE
ENDIF
REM TRANSFORMATION HELIOZENTR. EKLIPTIKEBENE ZUM BEOBACHTUNGSZENTRUM
x=rr*COS(eb)*COS(el)+rb*COS(-br)*COS(le+PI)
y=rr*COS(eb)*SIN(el)+rb*COS(-br)*SIN(le+PI)
z=rr*SIN(eb)+rb*SIN(-br)
efb=SQR(x^2+y^2+z^2) !ENTFERNUNG GESTIRN-BEOBACHTER IN KM
bb=z/efb          !EKL. BREITE DES GESTIRNS FÜR DAS BEOBACHTUNGSZENTRUM
IF ABS(bb)>=1 THEN
bb=0.9999999999*SGN(bb)
ENDIF
bb=ASIN(bb)
x1=x/(efb*COS(bb))
y1=y/(efb*COS(bb))
IF x1=<-0.9999999999 THEN
ll=PI
ELSE
ll=FN z(ATN(y1/(1+x1))*2) !EKL. LÄNGE DES GESTIRNS FÜR DAS VEHIKEL- BZW. BEOBACHTUNGSZENTRUM
ENDIF
dd=SIN(bb)*SIN(-br)+COS(bb)*COS(-br)*COS(ll-(le+PI))
IF ABS(dd)>1 THEN
dd=0.99999999*SGN(dd)
ENDIF
REM GNOMONISCHE PROJEKTION
s=ATN(TAN(bb)/COS(ll-(le+PI)))
x1=-fp*((COS(s)*TAN(ll-(le+PI)))/COS(s-(-br)))*ska
y1=-fp*TAN(s-(-br))*ska
IF dd<0 THEN
x1=5000
y1=5000
PLOT xo+x1,yo+y1
ENDIF
rx=SQR(x1*x1+y1*y1)
x=xo+x1
y=yo+y1
IF ns=1 THEN
ns=0
PLOT x,y
ENDIF
no=FN z(el-le)
IF rx<r1 AND no>PI/2 AND no<PI*1.5 THEN
c=POINT(x,y)
COLOR c
PLOT x,y
ENDIF
COLOR 1
IF x>-10 AND x<650 AND y>-10 AND y<410 THEN !AUF DER 640x400 GRAFIKFLÄCHE LINIEN, SONST PUNKTE
DRAW  TO x,y
ELSE
PLOT x,y
ENDIF
IF e<1 THEN
EXIT IF nn>6.29
ENDIF
IF e=1 OR e>1 THEN
EXIT IF nn>jj
ENDIF
LOOP
RETURN
PROCEDURE pb1
  REM TRANSFORMATION BAHNEBENE ZU EKLIPTIKEBENE
  IF e<1 THEN
    ex=ma
    FOR is=1 TO 200
      ex=ex-(ex-e*SIN(ex)-ma)/(1-e*COS(ex))
    NEXT is
    r=ra*(1-e*COS(ex))
    v=ATN(SQR((1+e)/(1-e))*TAN(ex/2))*2
  ENDIF
  IF e=1 THEN
    v=ATN((SQR(aa1^2+1)+aa1)^(1/3)-(SQR(aa1^2+1)-aa1)^(1/3))*2
    r=(2*q)/(1+COS(v))
  ENDIF
  IF e>1 THEN !HYPERBOLIDE BAHN
  IF ma<0 THEN
    ex=-(-6*ma)^(1/3)
  ELSE
    ex=(6*ma)^(1/3)
  ENDIF
  w1=1.0E-10
  REPEAT
    ww=ex
    ex=ex-(e*FN sinh(ex)-ex-ma)/(e*FN cosh(ex)-1)
  UNTIL ABS(ex-ww)<w1
  v=ATN(SQR((e+1)/(e-1))*FN tanh(ex/2))*2
  r=ABS(ra)*(e*FN cosh(ex)-1)
ENDIF
x=r*(COS(k)*COS(w+v)-SIN(k)*SIN(w+v)*COS(i)) !HELIOZENTR. KART. KOORDINATEN
y=r*(SIN(k)*COS(w+v)+COS(k)*SIN(w+v)*COS(i))
zv=r*(SIN(w+v)*SIN(i))
rr=SQR(x^2+y^2+zv^2)  !HELIOZENTR. ENTFERNUNG
eb=ASIN(zv/rr)  !HELIOZENTR. BREITE
o=COS(eb)
IF o==0 THEN
  o=1.0E-20
ENDIF
x2=x/(rr*o)
y2=y/(rr*o)
IF x2<=-1 THEN
  el=PI
ELSE
  el=FN z(ATN(y2/(1+x2))*2) !HELIOZENTR. LÄNGE
ENDIF
RETURN
PROCEDURE pb2
  REM TRANSFORMATION HELIOZENTR. EKLIPTIKEBENE ZUM BEOBACHTUNGSZENTRUM
  x=rr*COS(eb)*COS(el)+rb*COS(-br)*COS(le+PI)
  y=rr*COS(eb)*SIN(el)+rb*COS(-br)*SIN(le+PI)
  z=rr*SIN(eb)+rb*SIN(-br)
  efb=SQR(x^2+y^2+z^2) !ENTFERNUNG GESTIRN-BEOBACHTER IN KM
  bb=z/efb          !EKL. BREITE DES GESTIRNS FšR DAS BEOBACHTUNGSZENTRUM
  IF ABS(bb)>=1 THEN
    bb=0.9999999999*SGN(bb)
  ENDIF
  bb=ASIN(bb)
  x1=x/(efb*COS(bb))
  y1=y/(efb*COS(bb))
  IF x1=<-0.9999999999 THEN
    ll=PI
  ELSE
    ll=FN z(ATN(y1/(1+x1))*2) !EKL. LÄNGE DES GESTIRNS FÜR DAS VEHIKEL- BZW. BEOBACHTUNGSZENTRUM
  ENDIF
  dd=SIN(bb)*SIN(-br)+COS(bb)*COS(-br)*COS(ll-(le+PI))
  IF ABS(dd)>1 THEN
    dd=0.99999999*SGN(dd)
  ENDIF
  REM GNOMONISCHE PROJEKTION
  s=ATN(TAN(bb)/COS(ll-(le+PI)))
  x1=-fp*((COS(s)*TAN(ll-(le+PI)))/COS(s-(-br)))*ska
  y1=-fp*TAN(s-(-br))*ska
  IF dd<0 THEN
    x1=5000
    y1=5000
    PLOT xo+x1,yo+y1
  ENDIF
  rx=SQR(x1*x1+y1*y1)
  x=xo+x1
  y=yo+y1
  no=FN z(el-le)
  IF rx<r1 AND no>PI/2 AND no<PI*1.5 THEN
    c=POINT(x,y)
    COLOR c
    PLOT x,y
  ENDIF
  DEFFILL 1
  PCIRCLE x,y,3
RETURN
PROCEDURE jd
  jd=1720994.5+a2+FIX(30.6001*((a3-12*(a3<3))+1))+FIX(365.2425*(a4+(a3<3)))
RETURN
PROCEDURE praez !TRANSFORMATION PRÄZESSION BAHNELEMENTE
IF epo<1950 THEN !REDUKTION AUF ÄQUINOKTIUM B1950/FK4
to=(jd3-2451545)/36525
t2=(2433282.4235-jd3)/36525
ip=RAD(0.0131*t2)
lp=RAD(174.876+0.9137*to-0.2416*t2)
pn=RAD(1.397*t2)
i1=ACOS(COS(i)*COS(ip)+SIN(i)*SIN(ip)*COS(k-lp))
x=(-SIN(ip)*COS(i)+COS(ip)*SIN(i)*COS(k-lp))/SIN(i1)
y=(SIN(i)*SIN(k-lp))/SIN(i1)
k1=FN z(ATN(y/(1+x))*2+lp+pn)
x=(SIN(i)*COS(ip)-COS(i)*SIN(ip)*COS(k-lp))/SIN(i1)
y=(-SIN(ip)*SIN(k-lp))/SIN(i1)
w1=FN z(ATN(y/(1+x))*2+w)
i=i1
k=k1
w=w1
epo=1950
ENDIF
IF epo=1950 THEN !REDUKTION ÄQUINOKTIUM B1950/FK4 AUF J2000/FK5
j=RAD(0.00652)
l=RAD(5.19856)
l1=RAD(4.50002)
ij=ACOS(COS(i)*COS(j)-SIN(i)*SIN(j)*COS(l+k))
x=(SIN(i)*COS(j)+COS(i)*SIN(j)*COS(l+k))/SIN(ij)
y=(SIN(j)*SIN(l+k))/SIN(ij)
wj=FN z(ATN(y/(1+x))*2+w)
x=(COS(i)*SIN(j)+SIN(i)*COS(j)*COS(l+k))/SIN(ij)
y=(SIN(i)*SIN(l+k))/SIN(ij)
kj=FN z(ATN(y/(1+x))*2-l1)
i=ij
w=wj
k=kj
epo=2000
ENDIF
IF epo=2000 THEN !REDUKTION J2000/FK5 AUF ÄQUINOKTIUM DES DATUMS/FK5
ip=RAD(0.0131*t/3600)
lp=RAD(174.876-0.2416*t)
pa=RAD(1.397*t)
id=ACOS(COS(i)*COS(ip)+SIN(i)*SIN(ip)*COS(k-lp))
x=(-SIN(ip)*COS(i)+COS(ip)*SIN(i)*COS(k-lp))/SIN(id)
y=(SIN(i)*SIN(k-lp))/SIN(id)
kd=FN z(ATN(y/(1+x))*2+lp+pa)
x=(SIN(i)*COS(ip)-COS(i)*SIN(ip)*COS(k-lp))/SIN(id)
y=(-SIN(ip)*SIN(k-lp))/SIN(id)
wd=FN z(ATN(y/(1+x))*2+w)
i=id
w=wd
k=kd
ENDIF
RETURN
PROCEDURE pr  !PRÄZESSION J2000/FK5
w1=RAD(((2306.22+1.4*t)*te)/3600)
w2=RAD(((2306.22+1.4*t)*te)/3600)
w3=RAD(((2004.311-1*t)*te)/3600)
do1=ASIN(SIN(w3)*COS(dro)*COS(arro+w1)+COS(w3)*SIN(dro))
x=(COS(w3)*COS(dro)*COS(arro+w1)-SIN(w3)*SIN(dro))/COS(do1)
y=(COS(dro)*SIN(arro+w1))/COS(do1)
aro1=FN z(ATN(y/(1+x))*2+w2)
loo=FN z(loo+RAD((1.397+0.0006173*t)*te))
RETURN
el: !MITTL. BAHNELEMENTE J2000 = genäherte Bahnelenemte Äquinoktium des Datums für einige Jahrzehnte
DATA Erdbahn,0,0,102.9,0.0167,1.000001
DATA Merkurbahn,7,48.3,29.125,0.2056,0.387
DATA Venusbahn,3.395,76.68,54.9,0.0068,0.723
DATA Marsbahn,1.85,49.56,286.5,0.0934,1.52368
DATA Jupiterbahn,1.3,100.464441,273.86687,0.0485,5.2026,Saturnbahn,2.489,113.6655
DATA 339.39126,0.05508,9.554909,Uranusbahn,0.7732,74.005,98.9
DATA 0.046295,19.218446,Neptunbahn,1.76995,131.7841,276.3396,0.324062
DATA 0.00899,30.110386
PROCEDURE karte
  r=32*DEG(1)
  lya=3.26156378
  raa=(180/PI)*3600
  rab=raa/15
  rac=977792.2217
  neu:
  CLS
  BOX 1,1,619,399
  DEFTEXT 1,0,0,32
  TEXT 220,40,"STERNKARTE"
  DO
    PRINT AT(21,12);"Taste >1< - Sphärische  Sternkarte"
    PRINT AT(21,13);"Taste >2< - Gnomonische Sternkarte"
    w$=UPPER$(INKEY$)
    IF w$="1" OR w$="2" THEN
      GOTO ee
    ENDIF
  LOOP
  ee:
  PRINT AT(21,12);SPC(40)
  PRINT AT(21,13);SPC(40)
  nr1=VAL(w$)
  w$=""
  DO
    PRINT AT(20,12);"Taste >1< - Koordinaten Horizontebene"
    PRINT AT(20,13);"Taste >2< - Koordinaten Äquatorebene"
    w$=UPPER$(INKEY$)
    IF w$="1" OR w$="2" THEN
      GOTO ee1
    ENDIF
  LOOP
  w$=""
  ee1:
  PRINT AT(20,12);SPC(40)
  PRINT AT(20,13);SPC(40)
  numr=VAL(w$)
  GOSUB sternb
  GOSUB jd
  PRINT AT(2,14);"Geograph. Breite (o,',''): ";
  INPUT "",gg$,gm,gs
  GOSUB dec
  bo=RAD(FN gr1(gg,gm,gs))
  PRINT AT(2,16);"Geograph. L„nge  (o,',''): ";
  INPUT "",gg$,gm,gs
  GOSUB dec
  lgeo=FN gr1(gg,gm,gs)
  GOSUB sternz
  zm=RAD(ru*15)
  GOSUB do
  CLS
  GOSUB j
  GOSUB alig
  CLS
  GOSUB mil
  GOSUB stern
  IF alig=1 THEN
    alig=0
    GOSUB stern
  ENDIF
  GOSUB metspur  !METEORBAHN EINZEICHNEN
  GOSUB metspur1
  GOSUB um !RADIANT MARKIEREN
  GOSUB ta
  DO
    w$=UPPER$(INKEY$)
    IF w$="H" AND numr=1 THEN
      GOSUB ho
      GOSUB ta
    ENDIF
    IF w$="N" THEN
      GOSUB netz
      GOSUB ta
    ENDIF
    IF w$="E" THEN
      GOSUB ekl
      GOSUB ta
    ENDIF
    IF w$="G" THEN
      GOSUB gal
      GOSUB ta
    ENDIF
    IF w$="S" THEN
      RUN
    ENDIF
    IF w$="I" THEN
      GOSUB info1
      GOSUB ta
    ENDIF
    IF w$="K" THEN
      GOSUB koo
      GOSUB ta
    ENDIF
  LOOP
RETURN
PROCEDURE metspur
  RESTORE met
  anz=22 !EINTRAG ANZAHL METEORSPUREN
  FOR i=1 TO anz
    READ d1,a1,d2,a2 !DEKLIN.,AR ANFANG, DEKLIN., AR ENDE DER METEORSPUR
    d1(i)=RAD(d1) !GRAD
    a1(i)=RAD(a1) !GRAD
    d2(i)=RAD(d2)
    a2(i)=RAD(a2)
  NEXT i
  FOR i=1 TO anz
    ll2=d1(i)
    ll1=a1(i)
    IF numr=1 THEN
      GOSUB hoh
      d=ho
      ar=az
    ELSE
      d=ll2
      ar=ll1
    ENDIF
    GOSUB mat
    x2=xx
    y2=yy
    h2=h
    ll2=d2(i)
    ll1=a2(i)
    IF numr=1 THEN
      GOSUB hoh
      d=ho
      ar=az
    ELSE
      d=ll2
      ar=ll1
    ENDIF
    GOSUB mat
    IF h2>0 OR h>0 THEN
      LINE x2,y2,xx,yy
    ENDIF
  NEXT i
  met:
  REM d1,ar1,d2,ar2,d1,ar1,d2,ar2 EINTRAG DEKLINATION U. REKTASZENSION ANFANG/ENDE METEORSPUREN
  DATA 59,69,59,88,58,75,56,92,62,65,64,81,63,38,68,32,57,28,57,14,52,72,48,83,42,71,36,79,38,54,30,57,36,40,29,38,38,26,29,17
  DATA 46,21,35,7,49,14,45,6,33,67,29,70,28,47,19,48,53,16,50,2,44,7,41,1,31,35,25,33,31,45,23,45
  DATA 30,21,20,14,27,13,22,9,26,9,24,7,45,84,42,91,52,87,48,95
RETURN
PROCEDURE metspur1
  d1(1)=da
  a1(1)=ara
  d2(1)=de
  a2(1)=are
  d1(2)=dba
  a1(2)=arba
  d2(2)=dbe
  a2(2)=arbe
  FOR i=1 TO 2
    ll2=d1(i)
    ll1=a1(i)
    IF numr=1 THEN
      GOSUB hoh
      d=ho
      ar=az
    ELSE
      d=ll2
      ar=ll1
    ENDIF
    GOSUB mat
    x2=xx
    y2=yy
    h2=h
    ll2=d2(i)
    ll1=a2(i)
    IF numr=1 THEN
      GOSUB hoh
      d=ho
      ar=az
    ELSE
      d=ll2
      ar=ll1
    ENDIF
    GOSUB mat
    x3=xx
    y3=yy
    IF h2>0 THEN
      LINE x2,y2,x3,y3
    ENDIF
  NEXT i
RETURN
PROCEDURE do
  IF numr=2 THEN
    PRINT AT(2,18);"Eintrag Deklination Plattenmitte (o,','') z. B. Radiant: ";
  ELSE
    PRINT AT(2,18);"Eintrag H”he  Plattenmitte  (o,','') z. B. Radiant: ";
  ENDIF
  INPUT "",gg$,gm,gs
  GOSUB dec
  do=RAD(FN gr1(gg,gm,gs))+1.0E-07
  IF numr=2 THEN
    PRINT AT(2,19);"Eintrag          AR Plattenmitte (o,','') z. B. Radiant: ";
  ELSE
    PRINT AT(2,19);"Eintrag Azimut Plattenmitte (o,','') z. B. Radiant: ";
  ENDIF
  INPUT "",gg$,gm,gs
  GOSUB dec
  IF numr=1 THEN
    aro=2*PI-(RAD(FN gr1(gg,gm,gs))+1.0E-07)
  ELSE
    aro=(RAD(FN gr1(gg,gm,gs))+1.0E-07)
  ENDIF
RETURN
PROCEDURE alig
  BOX 160,130,430,190
  DO
    PRINT AT(24,10);"Taste >1< - mit  Alignemente"
    PRINT AT(24,11);"Taste >2< - ohne Alignemente"
    w$=UPPER$(INKEY$)
    IF w$="1" OR w$="2" THEN
      GOTO jump1
    ENDIF
  LOOP
  jump1:
  alig=0
  IF w$="1" THEN
    alig=1
  ENDIF
RETURN
PROCEDURE dec
  i=1
  IF LEFT$(gg$,1)="-" THEN
    i=-1
  ENDIF
  gg=ABS(VAL(gg$))
RETURN
PROCEDURE ho
  COLOR 1
  FOR az=0 TO PI*2 STEP RAD(10)
    n=1
    FOR h1=PI/2 TO -PI/2-0.05 STEP -0.05
      d=h1
      ar=az
      GOSUB mat
      IF n=1 THEN
        n=0
        c=POINT(xx,yy)
        COLOR 0
        PLOT xx,yy
      ENDIF
      IF xx<-10 OR xx>650 OR yy<-10 OR yy>420 THEN  !Eckpunkte Graphikfl„che
      PLOT xx,yy
    ENDIF
    COLOR 1
    IF h>=0 THEN
      DRAW  TO xx,yy
    ENDIF
  NEXT h1
NEXT az
FOR h1=PI/2 TO -PI/2 STEP -RAD(10)
  n=1
  FOR az=0 TO PI*2+0.05 STEP 0.05
    d=h1
    ar=az
    GOSUB mat
    IF n=1 THEN
      n=0
      c=POINT(xx,yy)
      COLOR 0
      PLOT xx,yy
    ENDIF
    IF xx<-10 OR xx>650 OR yy<-10 OR yy>420 THEN
      PLOT xx,yy
    ENDIF
    COLOR 1
    IF h>=0 THEN
      DRAW  TO xx,yy
    ENDIF
  NEXT az
NEXT h1
RETURN
.
.
.
Fortsetzung
 
 
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