function X = gauss_real(mu,sigma)  


U = rand(1);
  
V = 2*U-1;
S = V*V;
X = mu + (V*sqrt(-2*log(S)/S)) * sqrt(sigma);
end

 