Path: utzoo!utgpu!news-server.csri.toronto.edu!cs.utexas.edu!sdd.hp.com!spool.mu.edu!uwm.edu!psuvax1!rutgers!modus!otello!gear!wolf!Tom.Metcalf
From: Tom.Metcalf@wolf.fidonet.sublink.org (Tom Metcalf)
Newsgroups: comp.sys.handhelds
Subject: Celestial Navigation with HP-48SX
Message-ID: <3897.27A7599E@wolf.fidonet.sublink.org>
Date: 31 Jan 91 00:33:47 GMT
Organization: Rekursive Labs sdf- Pisa, Italy
Lines: 206

MESS
            END
          UNTIL 'HP
<1.2'
\v/ / 57.3 H1 H2 + 2
/ COS * * 225 D1 D2
+ 2 / COS SQ * DH
DT / SQ - \v/ / "ERR"
\->TAG
      \>>
    \>>
  DRLAT
35.1633333337
    \>>
  LOP
    \<< \-> g d a l
      \<<
        IF 'ABS(l)\=/
90'
        THEN 'g+
ACOS((SIN(a)-SIN(l)
*SIN(d))/(COS(l)*
COS(d)))' \->NUM
        ELSE g
        END DUP IM
        IF 0 \=/
        THEN DROP g
        END
        IF 'ABS(l)>
90-ABS(d)+a'
        THEN 180 +
        END LON
RANGE l R\->C
      \>>
    \>>
  CST { SOLVE ADDOB
SETUP INIT ADV
ADDDR DR PLOTP SAIL
WVIEW ERROR TIME }
  REFRACT
    \<< 0 \-> h rp
      \<< '1/TAN(h+
7.31/(h+4.4))' \->NUM
'rp' STO 'rp*((
PRESS-80)/930)/(1+
.00008*(rp+39)*(
TMPTR-10))' \->NUM 60
/
      \>>
    \>>
  MESS
    \<< 3 DISP 7
FREEZE 0 WAIT DROP
    \>>
  PPAR {
(66.9101324376,39.9691000001)
(66.5090675652,40.1191000001)
X 0 (0,0) FUNCTION
Y }
  T\Gg 6
  G\Gg 231.103333334
  PRESS 1010
  TMPTR 10
  a0 '-(G12*G23-G13
*G22)*G13+(G11*G23-
G12*G13)*G23-(G11*
G22-G12^2)*G33'
  a1 'G11*G22-G12^2
+G11*G33-G13^2+G22*
G33-G23^2'
  TF 51.8827999999
  CRS 320
  SPD 4.06448990652
  EV3 '-2*\v/Q*COS((\Gh
+360)/3)+N/3'
  EV2 'N-\Ga1-\Ga3'
  EV1 '-2*\v/Q*COS(\Gh/
3)+N/3'
  \Gm
6.00499443433E-3
  \Gb3 5.86314838322
  \Gb2 -.662967085927
  \Gb1
-1.64470052966E-2
  E3
[ .493568602538 .624304833881 .605502691141 ]
  E2
[ .153840003259 .62256501051 -.767297895922 ]
  E1
[ -.855992574723 .471864686196 .211235484564 ]
  INTERP
    \<< \-> T V1 V2
      \<< V1 V2 V1 -
T2 T1 - / T T1 - *
+
      \>>
    \>>
  GSUM
    \<< \-> DS DC GS GC
HS
      \<< 0 'G11' STO
0 'G12' STO 0 'G13'
STO 0 'G22' STO 0
'G23' STO { 3 } 0
CON 'R' STO OBS
OBJ\-> OBJ\-> DROP DROP
'N' STO 1 N
        START SIN
'HS' STO DUP SIN
'DS' STO COS 'DC'
STO DUP SIN 'GS'
STO COS 'GC' STO DS
SQ 'G11' STO+ DS DC
GC * * 'G12' STO+
DS DC GS * * 'G13'
STO+ DC SQ GC SQ *
'G22' STO+ DC SQ GS
GC * * 'G23' STO+ R
OBJ\-> DROP DC GS HS
* * + ROT DS HS * +
ROT DC GC HS * * +
ROT { 3 } \->ARRY 'R'
STO
        NEXT N G11
G22 + - 'G33' STO
      \>>
    \>>
  OUT
    \<< OBJ\-> DROP \-> U
V W
      \<<
        IF 'ABS(U)>
1'
        THEN U SIGN
'U' STO
        END U ASIN
V W R\->C ARG \->FMT
"LON" \->TAG SWAP
\->FMT "LAT" \->TAG
      \>>
    \>>
  UVW
    \<< \Gb1 \Ga1 \Gm - /
E1 * \Gb2 \Ga2 \Gm - / E2
* \Gb3 \Ga3 \Gm - / E3 *
+ +
    \>>
  EIGEN
    \<< \-> EV
      \<< 'G12*G23-
G13*G22+G13*EV'
\->NUM 'G13*G12-G11*
G23+G23*EV' \->NUM '
G11*G22-SQ(G12)-(
G11+G22)*EV+SQ(EV)'
\->NUM { 3 } \->ARRY
DUP ABS
        IF DUP 0 \=/
        THEN /
        ELSE DROP
        END
      \>>
    \>>
  \Ga2 2.6356149168
  \Ga3 6.29496706443
  \Ga1 .06941801877
  \Gh 'ACOS(R1/Q^1.5)
'
  R1 'A0/2+N/3*(A1/
6-Q)'
  Q '(N/3)^2-A1/3'
  N 9
  A0 -1.1517219228
  A1 17.2110524034
  G33 3.86275122638
  R
[ 2.80595360979 3.23989100557 4.05537118357 ]
  G23 1.1275173595
  G22 3.49049143543
  G13 1.55763290776
  G12 2.16410219514
  G11 1.64675733819
  GHA2
60.5550000011
  DEC2
22.0816666668
  T2 12
  GHA1
45.5566666678
  DEC1
22.0750000002
  T1 11
  LU 1
  SEMI .26333333386
  HP .9333333352
  HGT 3.048
  INDX 0
  BODY "S"
END
''
--------------------------------- CUT HERE ------------------
--  
 WolfNet BBS Pisa (Italy) Tel. +39-50-589050 300-14.4K Baud Matrix 2:332/602.0
 Tom Metcalf - via FidoNet node 2:332/602
 UUCP: ...!gear!wolf!Tom.Metcalf
 ARPA: Tom.Metcalf@wolf.fidonet.sublink.org