Simplified the numerical scheme in phugoid.py

Instead of using a rotation over angle dtheta around the center of
the osculating circle we just translate a distance ds along the tangent.