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