  We extend the gauge invariant variational approach of Phys. Rev. D52 (1995)
3719, hep-th/9408081, to theories with fermions. As the simplest example we
consider the massless Schwinger model in 1+1 dimensions. We show that in this
solvable model the simple variational calculation gives exact results.
