fn gauss_legendre_grids(
    score: [f64; 2],
    num_grids: usize
) -> (Vec<f64>, Vec<f64>)