Numerical integration using Simpson's method