mathpack::fitting
pub fn weighted_least_squares( f: Function, params: &[f64], vars: &[Vec<f64>], y: &[f64], weights: &[f64]) -> f64