%!demo
%! ## Compute confidence intervals for parameters of a fitted Birnbaum-Saunders
%! ## distribution
%! pd = makedist ("BirnbaumSaunders", "beta", 1, "gamma", 0.5)
%! randg ("seed", 21);
%! data = random (pd, 1000, 1);
%! pd_fitted = fitdist (data, "BirnbaumSaunders")
%! ci = paramci (pd_fitted, "Alpha", 0.05)
%!
%! ## Use this to obtain confidence intervals for the estimated parameters (beta
%! ## and gamma), providing a range of plausible values given the data.
