Description: Bias-precession-nutation matrix given CIP.
The source code for the Fortran routine BPN2XY is available: