updated gpdRangeFit to let user tell the function how to compute covariance. This is because in some cases the observed info is singular and whilst the numerical approximation is probably poor it at least lets us produce the required plot