Code covered by the BSD License  

Highlights from
Generation of Random Variates

image thumbnail

Generation of Random Variates



generates random variates from over 870 univariate distributions

stacy_pdf(x, alpha, bet, theta, xtrem)
% stacy_pdf.m - evaluates a Stacy Probability Density.
%   See "The Amoroso Distribution", G. E. Crooks, Tech. Note 003v5 (2010-01-14)
%   Vector form of the PDF!!!
%  Created by Jim Huntley,  02/17/10

function [pdf] = stacy_pdf(x, alpha, bet, theta, xtrem)

%persistent coef

    coef = exp(log(abs(bet/theta)) - gammaln(alpha));

argx = x ./ theta;
pdf = coef .* argx.^(alpha*bet-1) .* exp(-(argx.^bet));


Contact us