Function mean_se 

pub fn mean_se(x: &[f64], mult: f64) -> (f64, f64, f64)

(y, ymin, ymax) = mean ± mult standard errors (mean_se(x, mult)).