Numerical integration using Simpson's rule