We show how Mathematica can be used to obtain numerical solutions of integral equations by exploiting a combination of iteration and interpolation. The efficacy of the method is demonstrated by considering three classical integral equations of applied mathematics: Love's equation for the condenser problem, Theodorsen's equation associated with conformal mapping, and Nekrasov's equation arising in the theory of water waves. The success of the approach depends on the use of an appropriate method for the interpolation.