y = (1/sigma(sq rt(2 pi))) e -(((x-mu)^2)/(2sigma^2))