home *** CD-ROM | disk | FTP | other *** search
/ Collection of Education / collectionofeducationcarat1997.iso / SCIENCE / AA_51.ZIP / DEFLEC.C < prev    next >
C/C++ Source or Header  |  1993-02-09  |  610b  |  33 lines
  1. /* Correct for light deflection due to solar gravitation.
  2.  * AA page B37
  3.  *
  4.  * SE, pq, etc. were computed earlier by angles.c.
  5.  */
  6.  
  7. #include "kep.h"
  8.  
  9. int relativity( p, q, e )
  10. /* unit vector from earth to object:
  11.  */
  12. double p[];
  13. /* heliocentric ecliptic rectangular coordinates
  14.  * of earth and object:
  15.  */
  16. double q[], e[];
  17. {
  18. double C;
  19. int i;
  20.  
  21. C = 1.974e-8/(SE*(1.0+qe));
  22. for( i=0; i<3; i++ )
  23.     {
  24.     dp[i] = C*(pq*e[i]/SE - ep*q[i]/SO);
  25.     p[i] += dp[i];
  26.     }
  27. if( prtflg )
  28.     printf( "elongation from sun %.2f degrees, ", acos( -ep )/DTR );
  29. showcor( "light defl.", p, dp );
  30. return(0);
  31. }
  32.  
  33.