Function polymers::math::integrate_2d_grid

pub fn integrate_2d_grid(
    f: &dyn Fn(&f64, &f64) -> f64,
    grid_x: &[f64],
    grid_y: &[f64],
    da: &f64,
) -> f64