Orientation of the earth by numerical integration

A fundamental problem is the determination of the orientation of the earth in the celestial coordinate system. Classical reductions for precession and nutation can be expected to be consistent with the present-day observations, however, corrections to the classical theory are difficult to model because of the large number of coefficients involved. Consequently, a portion of the research has been devoted to numerically integrating the Eulerian equations of motion for a rigid earth and considering the six initial conditions of the integration as unknowns. Comparison of the three adjusted Eulerian angles from the numerical integration over 1000 days indicates agreement with classical theory to within 0.003 seconds of arc.