solve (x-d1)^2 + (sqrt(r2^2-x^2)-d2)^2 = r1^2