use assert_approx_eq::assert_approx_eq;
use sequential_integration::calculate_single_integral_simpson;

#[test]
fn calculate_single_integral_simpson_not_const_equation() {
    let equation = |x: f64| (1. - x.powf(2.)).sqrt();

    let first_integral_begin = -1.;
    let first_integral_end = 1.;
    let first_integral_step = 0.05;

    let result = calculate_single_integral_simpson(
        equation,
        first_integral_begin,
        first_integral_end,
        first_integral_step,
    )
    .unwrap();

    println!(
        "result: {}, expected: {}, diff: {}",
        result,
        std::f64::consts::FRAC_PI_2,
        result - std::f64::consts::FRAC_PI_2
    );

    assert_approx_eq!(result, std::f64::consts::FRAC_PI_2, 1e-2);
}