Description: Bias-precession-nutation matrix given CIP.
The source code for the ANSI C routine Bpn2xy is available: