A = X1^2
B = A^2
C = Y1^2
D = C^2
E = Z1^2
t0 = E^2
F = d1*t0
t1 = B+D
G = d2d1*t1
H = A*E
I = C*E
J = H+I
t2 = d2*J
K = G+t2
t3 = F+J
Z3 = t3+G
t4 = K+H
X3 = t4+D
t5 = K+I
Y3 = t5+B