This paper completes and improves the previous one by the author and J. A. Belward
[Comput. Math. Appl. 30, No. 7, 5-14 (1995; Zbl 0834.65004
)] and the one by the author [ibid. 30, No. 7, 15-19 (1995; Zbl 0834.65005
)], dealing with the application of a particular version of the Lanczos
-methods to approximate the Bessel functions
in the complex plane. That version uses Faber polynomials up to degree 15 as the perturbation terms. The author introduces symbolic representation of the scaled Faber polynomials and the appropriately modified (automated)
-method. Numerical comparisons, showing the accuracy improvements achieved by this new version of the
-method, are given and discussed.