The inputs must be numeric and a single value.
For this function to be well defined, xi<1.
All parameters must be provided.
This function can be optimised to find the fitted value of the location parameter mu for the GPD if the scale and shape parameters are known.