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