solve d2^2=(sqrt(d1^2 - (y3 - y1)^2) + x1 -x2)^2 + (y3 - y2)^2 for y3