pub fn polynomial_regression(x: &[f64], y: &[f64], degree: usize) -> Vec<f64>
Polynomial regression of given degree. Returns coefficients [a0, a1, …, a_deg].