$$a_{1} = \frac{x^{2} + y^{2} - \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} + 8}{4 x}$$
=
$$\frac{x^{2} + y^{2} - \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} + 8}{4 x}$$
=
0.25*(8 + x^2 + y^2 - 8*(1 + 0.015625*x^4 + 0.015625*y^4 + 0.25*y^2 - 0.25*x^2 + 0.03125*x^2*y^2)^0.5)/x
$$b_{1} = \frac{x^{2} + y^{2} - \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} - 8}{4 y}$$
=
$$\frac{x^{2} + y^{2} - \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} - 8}{4 y}$$
=
0.25*(-8 + x^2 + y^2 - 8*(1 + 0.015625*x^4 + 0.015625*y^4 + 0.25*y^2 - 0.25*x^2 + 0.03125*x^2*y^2)^0.5)/y
$$a_{2} = \frac{x^{2} + y^{2} + \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} + 8}{4 x}$$
=
$$\frac{x^{2} + y^{2} + \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} + 8}{4 x}$$
=
0.25*(8 + x^2 + y^2 + 8*(1 + 0.015625*x^4 + 0.015625*y^4 + 0.25*y^2 - 0.25*x^2 + 0.03125*x^2*y^2)^0.5)/x
$$b_{2} = \frac{x^{2} + y^{2} + \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} - 8}{4 y}$$
=
$$\frac{x^{2} + y^{2} + \sqrt{x^{4} + 2 x^{2} y^{2} - 16 x^{2} + y^{4} + 16 y^{2} + 64} - 8}{4 y}$$
=
0.25*(-8 + x^2 + y^2 + 8*(1 + 0.015625*x^4 + 0.015625*y^4 + 0.25*y^2 - 0.25*x^2 + 0.03125*x^2*y^2)^0.5)/y