%!demo
%! ## Compute the interquartile range for a Half-normal distribution
%! data = abs (randn (10000, 1));
%! pd = fitdist (data, "HalfNormal");
%! iqr_value = iqr (pd)
%!
%! ## Use this to calculate the interquartile range, which measures the spread
%! ## of the middle 50% of the distribution, helpful for understanding central
%! ## variability in positive data.
