t0 = X1+Y1
U1 = t0^2
U2 = Y1^2
S1 = X1^2
V1 = S1^2
aS1 = a2*S1
t1 = U2-V1
T = t1-aS1
t2 = U1-U2
X3 = t2-S1
t3 = T+V1
t4 = U2+aS1
t5 = 4*V1
t6 = t3*t4
Y3 = t6+t5
Z3 = T-V1