%%HP: T(3)A(D)F(.);
@ ANALEMMA by Craig A. Finseth
DIR
  GO
    \<< RCLF DEG
ERASE { # 0h # 0h }
PVIEW FRAME BASE
TICK LBL 7 FREEZE
STOF
    \>>
  BASE
    \<< 0 PLACE 0 365
      FOR d d PLACE
DUP ROT LINE 5
      STEP DROP
    \>>
  E -3.19866112964
  FIGE
    \<< L0 .0057183 -
\Ga - \GDp \Ge COS * +
'E' STO
    \>>
  FIGJD
    \<< 78 +
2451545.5 + 'JD'
STO
    \>>
  FIGL0
    \<< T 10 / \-> \Gt
      \<< 280.4664567
360007.698278 \Gt * +
.03032028 \Gt SQ * +
360 MOD 'L0' STO
      \>>
    \>>
  FIGT
    \<< JD 2451545 -
36525 / 'T' STO
    \>>
  FIG\Ga
    \<< 357.5291
35999.0503 T * +
360 MOD \-> M
      \<< 1.9146
.004817 T * -
.000014 T SQ * - M
SIN * .019993
.000101 T * - M 2 *
SIN * + .00029 M 3
* SIN * + L0 + \-> OD
        \<< 125.04
1934.136 T * - 360
MOD \-> \GW
          \<< OD
.00569 - .00478 \GW
SIN * - \-> \Gl
            \<< \Ge
.00256 \GW COS * +
COS \Gl SIN * \Gl COS
SWAP R\->C ARG 360 +
360 MOD
            \>>
          \>>
        \>>
      \>> '\Ga' STO
    \>>
  FIG\GD
    \<< 125.04452
1934.136261 T * -
360 MOD 280.4665
36000.7698 T * +
360 MOD 218.3165
481267.8813 T * +
360 MOD \-> \GW L LP
      \<< -.00172
HMS\-> \GW SIN *
.000132 HMS\-> 2 L *
SIN * - .000023
HMS\-> 2 LP * SIN * -
.000021 HMS\-> 2 \GW *
SIN * + '\GDp' STO
.00092 HMS\-> \GW COS *
.000057 HMS\-> L 2 *
COS * + .00001 HMS\->
LP 2 * COS * +
.000009 HMS\-> \GW 2 *
COS * - '\GD\Ge' STO
      \>>
    \>>
  FIG\Gd
    \<< T 36525 * 78
- 365 / 360 * 360
MOD SIN \Ge * '\Gd' STO
    \>>
  FIG\Ge
    \<< 23.2621448
HMS\-> .0046815 HMS\->
T * - .000000059
HMS\-> T SQ * - \GD\Ge +
'\Ge' STO
    \>>
  FRAME
    \<< (-51.17,-20)
(51.17,-20) LINE
(-51.17,-10)
(51.17,-10) LINE
(-51.17,0)
(51.17,0) LINE
(-51.17,10)
(51.17,10) LINE
(-51.17,20)
(51.17,20) LINE
(0,-25) (0,25) LINE
    \>>
  JD 2451967.5
  L0 336.902466384
  LBL
    \<< PICT (-32,1)
"Mar" 2 \->GROB GXOR
PICT (-32,16) "Aug"
2 \->GROB GXOR PICT
(-32,25) "Jul" 1
\->GROB GXOR PICT
(27,25) "Jun" 1
\->GROB GXOR PICT
(27,19) "May" 1
\->GROB GXOR PICT
(27,13) "Apr" 1
\->GROB GXOR PICT
(27,4) "Sep" 2
\->GROB GXOR PICT
(27,-6) "Oct" 2
\->GROB GXOR PICT
(27,-17) "Nov" 1
\->GROB GXOR PICT
(27,-22) "Dec" 1
\->GROB GXOR PICT
(-32,-19) "Jan" 1
\->GROB GXOR PICT
(-32,-9) "Feb" 2
\->GROB GXOR
    \>>
  PLACE
    \<< FIGJD FIGT
FIG\GD FIG\Ge FIGL0
FIG\Ga FIGE FIG\Gd E 4
* \Gd R\->C
    \>>
  PPAR {
(-51.17,-25)
(51.17,25) X 0
(0,0) PARAMETRIC Y
}
  T
1.15674195756E-2
  TICK
    \<< { 0 10 40 71
91 101 132 163 184
193 224 254 274 285
316 344 } DUP SIZE
\-> d n
      \<< 1 n
        FOR j d j
GET PLACE PIXOFF
        NEXT
      \>>
    \>>
  \Ga 340.091306964
  \GD\Ge
-4.10799632761E-4
  \GDp
-4.47118519728E-3
  \Gd -8.10071364255
  \Ge 23.4387298868
END