compute lambda = (Y1 + Y2)/(X1 + X2)
compute X3 = lambda^2 + lambda + X1 + X2 + a2
compute Y3 = lambda(X1 + X3) + X3 + Y1