Solve ((A - (r^2 - (y - C)^2)^(1/2)) - B)^2/(a^2) + ((y - D)^2)/(b^2) = 1 for y