Function integral 

pub fn integral<F>(f: F, a: f64, b: f64, n: usize) -> f64
where
    F: Fn(f64) -> f64,

Numerical integral using trapezoidal rule

Arguments

• f - Function to integrate
• a - Lower bound
• b - Upper bound
• n - Number of intervals (default: 1000)

Example

use dodecet_encoder::calculus;

let f = |x: f64| x * x;
let integral = calculus::integral(f, 0.0, 2.0, 1000);
assert!((integral - 8.0/3.0).abs() < 0.01); // ∫x²dx = x³/3, from 0 to 2 = 8/3