I prefer to calculate in polar coordinates because there is a problem with calculating \theta only from (x,y). For example if I use \theta=\mbox{atan2}(y,x) then if there is any oscillation outside of [-\pi,\pi] the angle will wrap around and strange things happen.
Ok, I just tried that. If...