In this paper we give details of the numerical analysis of an integral equation of the form (1) y ( t ) = f ( t ) + ∫ 0 t s μ - 1 t μ y ( s ) d s μ > 0 . The distinctive feature of the equation is the presence of a singularity at t = 0 for all values of μ > 0 and at s = 0 for all values of t > 0 for 0 < μ < 1. This means that conventional analytical and numerical theory does not apply. In fact, for 0 < μ < 1, the equation has an infinite family of solutions. As background, we give details of the analytical results on existence and uniqueness of solutions to (1) and we give new results on the use of numerical schemes that will yield approximations to any specific solution. We conclude with some numerical results that show that our methods enable us to find approximations of any chosen order to any of the infinite class of exact solutions to (1).