Generate the matrix of precession between two epochs, using the old, pre-IAU1976, Bessel-Newcomb model, using Kinoshita's formulation (double precision)
BEP0 double beginning Besselian epoch BEP1 double ending Besselian epoch
RMATP double[3][3] precession matrix The matrix is in the sense v(bep1) = rmatp * v(bep0)
Kinoshita, H. (1975) 'Formulas for precession', SAO Special Report No. 364, Smithsonian Institution Astrophysical Observatory, Cambridge, Massachusetts.
slaDeuler Defined in slamac.h: DAS2R P.T.Wallace Starlink 30 October 1993