
use gkquad::prelude::{integral, integral2, DynamicY};

fn main() {
    println!("1D: {}", 2.0 * 1.0);

    println!(
        "2D: {}",
        2.0 * integral(|x: f64| (1.0 - x * x).sqrt(), -1.0..1.0)
            .unwrap()
            .estimate
    );
