In this paper, an approach is presented for the numerical evaluation of weakly, strongly, hyper- and super-singular boundary integrals which exist in the Cauchy principal value sense in two-dimensional problems. In this approach, the singularities involved in integration kernels are analytically removed by expressing the nonsingular parts of the integration kernels as polynomials of the distance r. A self-contained Fortran code is listed and described for implementation of the proposed approach. The attached code is also able to evaluate general regular integrals using Gaussian quadrature, which enables the code to evaluate any two-dimensional boundary integral. Some examples are provided to verify the correctness of the presented formulations and the included code.