  A fully selfconsistent treatment for gap generation and Coulomb screening in
excitonic insulators is presented. The method is based on the equations of
motion for the relevant dynamical variables combined with a variational
approach. Applying the theory for a model system of bilayer graphene, an
excitonic groundstate with a gap exceeding 10 meV is predicted.
