This solution does not rely on trigonometry. Instead, the stereographic projection is used: $$ (x,y) = \left(\frac{1-q^2}{1+q^2},\frac{2q}{1+q^2}\right) $$

The differential equations were obtained with Mathematica. They are considerably large.