Numerical analysis of some integral equations with singularities
In this thesis we consider new approaches to the numerical solution of a class of Volterra integral equations, which contain a kernel with singularity of non-standard type. The kernel is singular in both arguments at the origin, resulting in multiple solutions, one of which is differentiable at the origin. We consider numerical methods to approximate any of the (infinitely many) solutions of the equation. We go on to show that the use of product integration over a short primary interval, combined with the careful use of extrapolation to improve the order, may be linked to any suitable standard method away from the origin. The resulting split-interval algorithm is shown to be reliable and flexible, capable of achieving good accuracy, with convergence to the one particular smooth solution.