# 11.
# a)
# m=0, sigma^2=1
x <- rnorm(n)
m.hat <- mean(x)
m.hat
sd2.hat <- sum((x - mean(x)) ^ 2) / n
sd2.hat
sd.hat<-sqrt(sd2.hat)
sd.hat
# ili
(n - 1) / n * var(x)

fitdistr(x, "normal")