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))