Ab initio coupled Hartree-Fock calculation of the 31P NMR shielding constants in P_4, P_2, and PN