Ab initio calculation of the 31P NMR shielding tensor for the series PO_a F_b^{-(2a+b-5)}, a+b=4 and for HPO_4^{2-}