  We determine, by numerical simulations on a lattice, the gauge-invariant
two-point correlation functions of the gauge field strengths in the QCD vacuum
with four flavours of dynamical staggered fermions.
