Variational band theory of vibronic polarons in crystals. III. Numerical   calculations
  In the preceding Part II, we derived variational equations for the phonon Fourier amplitudes and for the Fourier amplitudes of the fractional contribution of the electronic bands to the trial variational state. These equations are now solved by means of iterations for each value of the total momentum in order to obtain the energy vs. momentum relation for the ground state. Another result is mapping out the phonon and band Fourier amplitudes in the parameter space of the mixing constant and the electron hopping energy. 
M. Georgiev, A.G. Andreev and M.D. Ivanovich