home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / x / volume19 / xephem / part21 / obliq.c < prev    next >
Encoding:
C/C++ Source or Header  |  1993-05-15  |  426 b   |  23 lines
  1. #include <stdio.h>
  2. #include "astro.h"
  3.  
  4. /* given the modified Julian date, mjd, find the obliquity of the
  5.  * ecliptic, *eps, in radians.
  6.  */
  7. void
  8. obliquity (mjd, eps)
  9. double mjd;
  10. double *eps;
  11. {
  12.     static double lastmjd = -10000, lasteps;
  13.  
  14.     if (mjd != lastmjd) {
  15.         double t;
  16.         t = mjd/36525.;
  17.         lasteps = degrad(2.345229444E1
  18.             - ((((-1.81E-3*t)+5.9E-3)*t+4.6845E1)*t)/3600.0);
  19.         lastmjd = mjd;
  20.     }
  21.     *eps = lasteps;
  22. }
  23.