Solution for exercise from Langtangen

import numpy as np

m= 0.0
s= 2.0
x= 1.0
print(1/np.sqrt(2*np.pi)/s*np.exp(-0.5*((x-m)/s)**2))