%!demo
%! ## Compute the interquartile range for a Poisson distribution
%! data = poissrnd (5, 10000, 1);
%! pd = fitdist (data, "Poisson");
%! iqr_value = iqr (pd)
%!
%! ## Use this to calculate the interquartile range, measuring the spread of the
%! ## middle 50% of the distribution, helpful for understanding variability in count data.
