We present molecular dynamics computations of the time-dependent auto-correlation function of the single-particle density in a classical one-component plasma for three thermodynamic states in the range of intermediate and strong coupling. The deviations from the Gaussian approximation are calculated and the data are analyzed by the standard memory function formalism.