On the Numerical Integration of Discontinuous Functions

The numerical integration of discontinuous functions is best carried out by integrating each continuous portion of the integrand separately. Through the use of a simple example, it is demonstrated that integrating through a discontinuity (by blindly applying a standard technique such as Simpson´s rule) increases the computation time and the sensitivity to round-off error.