We consider a system of N bosons interacting through a two-body potential with, possibly, Coulomb-type singularities. We show that the difference between the many-body Schrodinger evolution in the mean-field regime and the effective nonlinear Hartree dynamics is at most of the order 1/N, for any fixed time. The N-dependence of the bound is optimal.