ex4_1_4.py

# exercise 4.1.4

import numpy as np

# Number of samples
N = 1000

# Mean
mu = np.array([13, 17])

# Covariance matrix
S = np.array([[4,3],[3,9]])

# Generate samples from the Normal distribution
X = np.random.multivariate_normal(mu, S, N)

print('Ran Exercise 4.1.4')