%!demo
%! ## Compute the median for different Half-normal distributions
%! data1 = 0.5 * abs (randn (10000, 1));
%! data2 = 1.0 * abs (randn (10000, 1));
%! pd1 = fitdist (data1, "HalfNormal");
%! pd2 = fitdist (data2, "HalfNormal");
%! median1 = median (pd1)
%! median2 = median (pd2)
%!
%! ## Use this to find the median value, which splits the distribution
%! ## into two equal probability halves, robust to skewness in positive data.
