Symbolic computation on complex polynomial solution of differential equations

A symbolic computation scheme, based on the Lanczos tau-method, is proposed for obtaining exact polynomial solutions to some perturbed differential equations with suitable boundary conditions. The automated tau-method uses symbolic Faber polynomials as the perturbation terms for arbitrary circular sections of the complex plane and has advantages of avoiding rounding error and easy manipulation over the numerical counterpart. The method is illustrated by applying it to the modified Bessel function of the first kind I-0(z) and the quality of the approximation is discussed. (C) 1996 Academic Press Limited