Solve for y₁ d = square root of (x_2-x_1)^2+(y_2-y_1)^2