Modified homotopy perturbation method (HPM) was used
to solve the hypersingular integral equations (HSIEs) of the first kind on the interval [−1,1] with the assumption that the kernel of the hypersingular integral is constant on the diagonal of the domain. Existence of inverse of hypersingular integral operator leads to the convergence of HPM in certain cases. Modified HPM and its norm convergence are obtained in Hilbert space. Comparisons between modified HPM, standard HPM, Bernstein polynomials approach Mandal and Bhattacharya (Appl Math Comput 190:1707−1716, 2007), Chebyshev expansion method Mahiub et al. (Int J Pure Appl Math 69(3):265–274, 2011) and reproducing kernel Chen and Zhou (Appl Math Lett 24:636–641, 2011) are made by solving five examples. Theoretical and practical examples revealed that the modified HPM dominates the standard HPM and others. Finally, it is found that the modified HPM is exact, if the solution of the problem is a product of weights and polynomial functions. For rational solution the absolute error decreases very fast by increasing the number of collocation points.