Numerical integration using Simpson's method
