pub fn polyfit<T: RealField + Copy>( x_values: &[T], y_values: &[T], polynomial_degree: usize, ) -> Result<Vec<T>, &'static str>