pub fn probability_density(x: f64, avg: f64, stddev: f64) -> f64
probability_density normalize x using the mean and the standard deviation and return the PDF https://en.wikipedia.org/wiki/Probability_density_function
