This paper considers the numerical evaluation of the Principal Value of the improper integral over a region [a, b] which may include q. Methods of Newton–Cotes and modified Gaussian type are considered. Some of these require the evaluation of f(q), while others do not. The methods are illustrated by some numerical examples.