symulacje

Generowanie wartości z rozkładu normalnego uciętego – metoda eliminacji

# liczba generowanych wartości
n=100
# parametry rozkładu uciętego
z=c(); a=-2; b=2; d=dnorm(0)
licz=0
repeat
{
u1=runif(1,a,b)
u2=runif(1,0,d)
if (u2<dnorm(u1)) z=c(z,u1); licz=licz+1
if (licz>=n) break
}

# wyświetlenie wyniku
z