Numerical integration: The Cavalieri-Simpson method