We present a solution of the infinite dimensional Hubbard model obtained by the numerical renormalization group method. This method gives directly the dynamical response in real frequency at zero temperature, and also gives reliable estimation of static quantities such as the electronic specific heat near absolute zero. We note technical points which need due care in applying this method.