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.