Argh - Doppler cooling is too complicated

Since it looks like I will be spending the rest of my life of this, I might present the problem: Basically, we have a nice 1D model of Doppler cooling in a linear Paul trap. When we wanted to check that it was OK to ignore the remaining dimensions, we entered a world of pain (or rather: I entered a world of pain). – Even the arm waiving seems to take for ever!