def npdf ( x, mu = 0, sigma = 1 ):
   "normal probability density function"
   import math
   f = math.exp(-(x - mu)**2/(2*sigma**2))
   f = f/(sigma*math.sqrt(2*math.pi))
   return f