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