AbstractIt is suggested that the effects of change in the celestial coordinates due to general precession be computed by solving numerically the differential equations for the instataneous rates of change in right ascension and declination. An approximate solution in closed from is given. Similar differential equations may be solved for the complete rotation matrix, requiring only the knowledge of the variations of the two functionmandnwith time. This suggests the existence of a relation between the classical precessional elementsz, ζ0and Θ, and this relation is derived. Some numerical examples are also presente