Have attempted to add a bit to the MAP() routine that evaluates
the Hessian matrix at the end of the likelihood optimisation.
Not really sure if it's working properly though as the uncertainties
on the parameters that it implies don't always seem reliable...