function f = Laplace_pdf(X,muX,sigmaX)
% f = Laplace_pdf(X,muX,sigmaX)
% return the PDF of the Laplace distribution with mean muX
% and std.dev sigmaX, at points X
 
%    Reference: http://en.wikipedia.org/wiki/Laplace_distribution

  sr2 = sqrt(2);

  f =  1/(sr2*sigmaX) * exp(-sr2*abs(X-muX) / sigmaX);

% ------------------------------------------------------------- Laplace_pdf
