%!demo
%! ## Generate random samples from a Burr distribution
%! pd = makedist ("Burr", "alpha", 1, "c", 2, "k", 3)
%! rand ("seed", 21);
%! samples = random (pd, 50, 1);
%! hist (samples, 90)
%! title ("Histogram of 50 random samples from Burr(alpha=1, c=2, k=3)")
%! xlabel ("values in x")
%! ylabel ("Frequency")
%!
%! ## This generates random samples from a Burr distribution, useful for
%! ## simulating non-negative data like household income under specific parameters.
