  We introduce a numerical method for generating the approximating polynomials
used in fermionic calculations with smeared link actions. We investigate the
stability of the algorithm and determine the optimal weight function and the
optimal type of discretization. The achievable order of polynomial
approximation reaches several thousands allowing fermionic calculations using
the Hypercubic Smeared Link action even with physical quark masses.
