Iterated solutions of linear operator equations with the Tau method

The Tau Method produces polynomial approximations of solutions of differential equations. The purpose of this paper is (i) to extend the recursive formulation of this method to general linear operator equations defined in a separable Hilbert space, and (ii) to develop an iterative refinement procedure which improves on the accuracy of Tau approximations. Applications to Fredholm integral equations demonstrate the effectiveness of this technique.