Function sum_squared_deviations 

pub fn sum_squared_deviations(values: &[f64], mean: f64) -> f64

Compute the sum of (x - mean)^2 for all elements.