Multiloop Calculations of Beta-Function of \({N=1}\) Supersymmetric Theories, Regularized by Higher Derivatives

Using a special C++ program we investigate contributions to the \(\beta\) -function in \({N}=1\) SQED with \(N_{f}\) flavors from vacuum graphs in the form of integrals of double total derivatives in four loop approximation. We compare this result with the three-loop anomalous dimension calculated earlier. In this way NSVZ equation was checked in this order of perturbation theory.