import numpy.random as rd
import numpy as np
import matplotlib.pyplot as plt

def phi(x):
    return 1/np.sqrt(2*np.pi)*np.exp(-x**2/2)
x = np.linspace(-5,5,50)
plt.close()
plt.plot(x,phi(x))
N = 10000
X = ....
plt.hist(X,range=(-5,5),bins=50,density=True,color="yellow",edgecolor="red")
plt.show()