The Numerical Renormalization Group applied to the Bose Hubbard Model