%!demo
%! ## Compute the interquartile range for a Generalized Pareto distribution
%! pd = makedist ("GeneralizedPareto", "k", 0.5, "sigma", 1, "theta", 0)
%! iqr_value = iqr (pd)
%!
%! ## Use this to calculate the interquartile range, which measures the spread
%! ## of the middle 50% of the distribution, useful for understanding variability
%! ## in exceedance values.
