Systems of quasilinear elliptic equations with dependence on the gradient via subsolution-supersolution method