A numerical method for hypersingular integrals of the first kind