n = 10000000
x = rnorm(n,0,1)
a = x[x >= -1 & x <= 1]
sigma1 = length(a)/length(x)
print(sigma1)
a = x[x >= -2 & x <= 2]
sigma2 = length(a)/length(x)
print(sigma2)
a = x[x >= -3 & x <= 3]
sigma3 = length(a)/length(x)
print(sigma3)
n = 10000000
x = rnorm(n,0,1)
a = x[x >= -1 & x <= 1]
sigma1 = length(a)/length(x)
print(sigma1)
a = x[x >= -2 & x <= 2]
sigma2 = length(a)/length(x)
print(sigma2)
a = x[x >= -3 & x <= 3]
sigma3 = length(a)/length(x)
print(sigma3)