A = X1^2
t0 = Y1*Z1
B = A+t0
C = X1*Z1
D = C^2
t1 = B^2
t2 = a2*D
t3 = B*C
t4 = t1+t3
E = t4+t2
X3 = C*E
t5 = B+C
t6 = A^2
t7 = t6*C
t8 = t5*E
Y3 = t8+t7
Z3 = C*D