Comparison of model function fitting algorithms for the quantification of strongly overlapping peaks in 31P NMR spectra