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!