By using a model for the non-perturbative gluon propagator, suggested by analyticity, we determine the heavy quark potential and calculate the mass spectrum of the ccBAR and bbBAR systems. The same model reproduces the results of lattice calculations for the gluon condensate.