The scattering solution from N scatterers, each of whose scattered field is approximated by M harmonics is an NM unknown problem. A straight-forward solution to this problem can be obtained by casting it into an NM linear algebraic equation. The solution of the linear algebraic equation will involve order N3M3 flouting point operations. However, via the use of a recursive algorithm and the translation formula, an order N2M3 algorithm to solve such a problem is possible.