Form the matrix of nutation for a given date (IAU 1980 theory). (double precision)
Final report of the IAU working group on nutation, chairman P.K.Seidelmann, 1980. Kaplan, G.H., 1981, USNO circular no. 163, pA3-6.
date double TDB (loosely ET) as Modified Julian Date (=JD-2400000.5)
rmatn double[3][3] nutation matrix The matrix is in the sense v(true) = rmatn * v(mean) .
slaNutc, slaDeuler P.T.Wallace Starlink 14 July 1993