integSimpson(f, a, b, rel_eps, max_steps)

Calculates the integral of a function using the Simpson's method.

fThe function to be integrated.
aL.h.s. of the integration interval.
bR.h.s. of the integration interval.
rel_eps(Optional, defaults to 1e-10) precision parameter.
max_steps(Optional, defaults to 20) Max number of simpson interation, it implies a maximal 2^max_steps function evaluation.
return valueThe numerical integral.