# L-13 MCS 260 Wed 10 Feb 2010 : npdf

# The name of the function and module are the same.

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