Abstract
The product integration method is used for the numerical solution of weakly singular integral equations of the second kind. These equations often have solutions which have derivative singularities at the endpoints of the range of integration. Therefore, the order of convergence results of de Hoog and Weiss for smooth solutions do not hold in general. In this paper it is shown that their results may be regained for the general case by using an appropriate nonuniform mesh. The spacing of the knot points is defined by the behavior of the solution at the endpoints. If the solution is smooth enough the mesh becomes uniform. Numerical examples are given.