import numpy as np
from scipy.stats import norm
import matplotlib.pyplot as plt
def norm_dist_prob(theta):
y = norm.pdf(theta,loc=10,scale=5)
return y
x = np.arange(1,20,0.1)
y = [norm_dist_prob(x[i]) for i in range(len(x))]
plt.scatter(x,y)
plt.show()