Dependence of the heavy quark wave function renormalization on the c-quark mass in QCD

With use of numerical methods for calculation of higher-order corrections to the perturbation theory, results for the renormalization constant of the wave function of a heavy quark are obtained with allowance for contribution from c-quark loops. In calculations the Laporta’s automated algorithm, the sector decomposition techniques, and the Mellin-Barnes representation have been used.