%!demo
%! ## Compute the interquartile range for a Birnbaum-Saunders distribution
%! pd = makedist ("BirnbaumSaunders", "beta", 1, "gamma", 0.5)
%! 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 failure times.
