use numlib::integrate::*;
use std::f64::consts::PI;
fn main() {
println!("Demonstration of composite trapezoid rule and simpson rule over periodic intervals. We integrate sin(x) over 0 to 4pi. The solution should be 0.");
println!("Trapezoid solution: {}", composite_trapezoid(0.0, PI*4.0, 100, sin));
println!("Simpson's solution: {}", simpson_quadrature(0.0, PI*4.0, 100, sin));
}
fn sin(x : f64) -> f64 {
return x.sin();
}