palEcmat palEcmat - Form the equatorial to ecliptic rotation matrix - IAU 2006 precession model

Description:
The equatorial to ecliptic rotation matrix is found and returned. The matrix is in the sense V(ecl) = RMAT $*$ V(equ); the equator, equinox and ecliptic are mean of date.

Invocation:
palEcmat( double date, double rmat[3][3] )

Arguments:

date = double (Given)
TT as Modified Julian Date (JD-2400000.5). The difference between TT and TDB is of the order of a millisecond or two (i.e. about 0.02 arc-seconds).
rmat = double[3][3] (Returned)
Rotation matrix