t0 = 1/X1
t1 = Y1*t0
lambda = X1+t1
t2 = lambda^2
t3 = t2+lambda
X3 = t3+a2
t4 = X1+X3
t5 = lambda*t4
t6 = t5+X3
Y3 = t6+Y1