# exercise 4.1.1

import numpy as np
from matplotlib.pyplot import figure, hist, plot, show, subplot, title

# Number of samples
N = 200

# Mean
mu = 17

# Standard deviation
s = 2

# Number of bins in histogram
nbins = 20

# Generate samples from the Normal distribution
X = np.random.normal(mu, s, N).T
# or equally:
X = np.random.randn(N).T * s + mu

# Plot the samples and histogram
figure(figsize=(12, 4))
title("Normal distribution")
subplot(1, 2, 1)
plot(X, ".")
subplot(1, 3, 3)
hist(X, bins=nbins)
show()