Nonlinear spectroscopy with density matrix propagation