function [y,density] = randnorm(mu,sigma)
%
% y = randnorm(mu,sigma)
%
% Genera y distribuito come una v.a. normale di media mu e deviazione
% standard sigma.
%
% Calcola anche la densita` per mu(1),lambda(1). Puo` essere testato con
%
% testva(@randnorm,mu*ones(1,10000),sigma*ones(1,10000))
y = mu+sqrt(2)*sigma.*erfinv(2*rand(size(mu))-1);
%y = mu+sigma.*randn(size(mu));
density = @(x) 1/sqrt(2*pi*sigma(1)^2)*exp(-(x-mu(1)).^2./(2*sigma(1)^2));