intsim function

Numerical integration using Simpson's rule