The Cauchy Principal Value Integral of a user supplied function (also known as the Finte Hilbert Transform) is evaluated by a Double Exponential quadrature. This is just a trapezoidal integration after a change of variable. The quadrature points are chosen symmetrically around the point of evaluation in the new variable. Coding is due to Mohankumar and Natarajan.
