High order methods for weakly singular integral equations with nonsmooth input functions

To solve one-dimensional linear weakly singular integral equations on bounded intervals, with input functions which may be smooth or not, we propose to introduce first a simple smoothing change of variable, and then to apply classical numerical methods such as product-integration and collocation based on global polynomial approximations. The advantage of this approach is that the order of the methods can be arbitrarily high and that the associated linear systems one has to solve are very well-conditioned.