Code covered by the BSD License  

Highlights from
Generation of Random Variates

image thumbnail

Generation of Random Variates

by

 

generates random variates from over 870 univariate distributions

relbw_pdf(x, M, Gamm)
% relbw_pdf.m - evaluates a Relativistic Breit Wigner Probability Density.
%   See "Relativistic Breit-Wigner Distribution", Wikipedia 
%
%  Created by Jim Huntley,  08/64/11
%

function [pdf] = relbw_pdf(x, M, Gamm)

%persistent Msq Gammsq gam k

%if(isempty(Msq))
    Msq = M^2;
    Gammsq = Gamm^2;
    gam = sqrt(Msq*(Msq+Gammsq));
    k = 2*sqrt(2)*M*Gamm*gam / (pi*(sqrt(Msq+gam)));
%end

pdf = k ./ ((x.^2-Msq).^2+Msq*Gammsq);

return

Contact us