Normal Probability Distributions

> ND := (mu,sigma,x) -> 1/sqrt(2*Pi*sigma^2)*exp(-1/2*((x-mu)/sigma)^2);