We give a mathematically rigorous derivation of Ehrenfests equations for the evolution of position and momentum expectation values, under general and natural assumptions which include atomic and molecular Hamiltonians with Coulomb interactions.