#-------------------------------------------------------
# semplice generatore di numeri pseudocasuali gaussiani
# che fa uso di
# - proprietà di media e varianza di combinazioni lineari
#   -> da cui media a varianza di media aritmetica 
# - teorema del limite centrale
#
# GdA, giugno 2023 
#--------------------------------------------------------

myrnorm <- function (n, mu=0, sigma=1) {
    x <- rep(-6,n)
    for (i in 1:12) {
        x <- x + runif(n)
    }
    x <- sigma*x + mu
    return(x)
}

x <- myrnorm(100000, 2, 3)
hist(x, nc=100, col='cyan', prob=TRUE) # anche se in nome 'prob' fa schifo...
print(mean(x))
print(sd(x))