Simpson's Rule

Code used for approximating certain integrals with Simpson's Rule. Normally around 20 intervals should suffice to get 6 places of accuracy.