A polynomial approximation to Bessel functions that arises from an electromagnetic scattering problem is examined. The approximation is extended to Bessel functions of any integer order, and the relationship to the Taylor series is derived. Numerical calculations show that the polynomial approximation and the Taylor series truncated to the same order have similar accuracies.