pub fn gaussian(x: f64, mu: f64, sigma: f64) -> f64
Gaussian function
N(x|μ,σ) = 1/√(2πσ^2) exp(-(x-μ)^2/(2σ^2))