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.