use numlib::integrate::*;

fn main() {
    println!("Demonstration of the trapezoid quadrature and simpson quadrature of x^2 over [-1, 1]. The exact solution is 2/3.");
    println!("Trapezoid solution {}", composite_trapezoid(-1.0, 1.0, 100, square));
    println!("Simpson's solution {}", simpson_quadrature(-1.0, 1.0, 100, square)); 
    println!("3/8's Simpson's solution {}", three_eighths_simpson(-1.0, 1.0, 100, square));
}

fn square(x : f64) -> f64{
    return x*x;
}