<gh_stars>0
from math import *
def f(mu,sigma2,x):
    return 1/sqrt(2.*pi*sigma2)*exp(-0.5*(x-mu)**2 /sigma2)

print(f(10.,4.,8.))

# maximum gaussian will be when x =mu
from math import *
def f(mu,sigma2,x):
    return 1/sqrt(2.*pi*sigma2)*exp(-0.5*(x-mu)**2 /sigma2)

print(f(10.,4.,8.))

# maximum gaussian will be when x =mu
# maximum gaussian will be when x =mu