Harmonic and biharmonic boundary value problems (BVP) arising in physical situations in fluid mechanics are, in general, intractable by analytic techniques. In the last decade there has been a rapid increase in the application of integral equation techniques for the numerical solution of such problems 1,2,3]. One such method is the boundary integral equation method (BIE) which is based on Green's Formula 4] and enables one to reformulate certain BVP as integral equations. The reformulation has the effect of reducing the dimension of the problem by one. Because discretisation occurs only on...