Function var_coefficients 

fn var_coefficients(
    trains_binned: &[Vec<f64>],
    order: usize,
) -> (Vec<f64>, Vec<f64>)

Fit VAR(order) model. Returns (beta [order*d × d, row-major], sigma [d×d, row-major]).