1 = (x2+y2)/(r2+a2) + z2/r2
r2(r2+a2) = r2(x2+y2) + (r2+a2)z2
r4 + r2a2 = r2x2 + r2y2 + r2z2 + a2z2
r4 + (a2-x2-y2-z2)r2 - a2z2 = 0
r2 = ((x2+y2+z2-a2) ± √((x2+y2+z2-a2)2+4a2z2))/2