The derivation of the Hartree equation from many-body systems of Bosons in the mean field limit has been very intensively studied in the last couple of years. However, very few results exist showing convergence of the k-th marginal of the N-body density matrix to the projection to the k-fold tensor product of the solution of the Hartree equation in stronger trace norms like the energy trace norm, see [MS], [Lu]. This issue is from a physical view point very important. The reason is that one can then approximate expectation values of certain observables of the N-body system by means of the Hartree equation, with relaxation of the very restrictive assumption that the observables are bounded operators. Here we consider the non-relativistic case. We prove, assuming only H1-regularity of the initial data, convergence in the energy trace norm without rates, and convergence in any other weaker Sobolev trace norm with rates. Our proof is simple and uses the functional aN introduced by Pickl in [Pi].