home *** CD-ROM | disk | FTP | other *** search
/ Simtel MSDOS - Coast to Coast / simteldosarchivecoasttocoast2.iso / astrnomy / ephem421.zip / OBLIQ.C < prev    next >
C/C++ Source or Header  |  1990-09-13  |  442b  |  22 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. obliquity (mjd, eps)
  8. double mjd;
  9. double *eps;
  10. {
  11.     static double lastmjd = -10000, lasteps;
  12.  
  13.     if (mjd != lastmjd) {
  14.         double t;
  15.         t = mjd/36525.;
  16.         lasteps = degrad(2.345229444E1
  17.             - ((((-1.81E-3*t)+5.9E-3)*t+4.6845E1)*t)/3600.0);
  18.         lastmjd = mjd;
  19.     }
  20.     *eps = lasteps;
  21. }
  22.