fn solve_ols(x: &[Vec<f64>], y: &[f64], p: usize) -> Vec<f64>
Solve OLS via normal equations with regularization.