####################
# Exercise 4.1.1
####################
rm(list = ls()) # Clear work space

# Number of samples
N <- 100

# Mean
mu <- 17

# Standard deviation
s <- 2

# Number of bins in histogram
NBins <- 20

# Generate samples from the Normal distribution
X <- rnorm(N, mean = mu, sd = s)

# Plot a histogram
{
  par(mfrow = c(1, 2))
  plot(X, main = "Data")
  hist(X, breaks = NBins, main = "Histogram of Data")
}