\displaystyle r_{1, 2} = \frac{-(x_0 + y_0) \pm \sqrt{(x_0 + y_0)^2 + 4 (1 - x_0 y_0)}}{2}