  A Bose-Hubbard model, describing bosons in a harmonic trap with a
superimposed optical lattice, is studied using a fast and accurate variational
technique (MF+NRG): the Gutzwiller mean-field (MF) ansatz is combined with a
Numerical Renormalization Group (NRG) procedure in order to improve on both.
Results are presented for one, two and three dimensions, with particular
attention to the experimentally accessible momentum distribution and possible
satellite peaks in this distribution. In one dimension, a comparison is made
with exact results obtained using Stochastich Series Expansion.
