ash/freed   oct 1968    subroutine monlib
koziel formulation of moon rotation-libration
modified sept 1972 by r.w.king to calculate angles for lunar
rotation integration

mnrt1= 0 lunorb or rotmat is calling routine
= 1 mrtfn is calling routine
ilib= 0 input angles are tau,rho,isig and derivatives
= 1 input angles are psi,theta,phi (mod) and derivatives
inert = 0 ecliptic rates are mod, noninertial
inert = 1 ecliptic rates are inertial w.r.t. epoch
moon rotation-libration matrix (output of this routine)