laplace.m

function y  = randlap(mu, sigma)
%   mean mu and standard deviation sigma. 
%   mu      : mean
%   sigma   : standard deviation

% Generate Laplacian noise
u = rand(1, 1)-0.5;
b = sigma / sqrt(2);