XX = X1^2
YY = Y1^2
ZZ = Z1^2
D = XX+YY
DD = D^2
t0 = ccd*ZZ
t1 = D-t0
E = 4*t1
t2 = YY-XX
t3 = D*t2
H = 2*t3
t4 = XX*E
P = DD-t4
t5 = YY*E
Q = DD-t5
QQ = Q^2
t6 = Q+X1
t7 = t6^2
t8 = H+Q
t9 = t7-QQ
t10 = t9-XX
X3 = t8*t10
t11 = H-P
t12 = P*Y1
t13 = t11*t12
Y3 = 2*t13
t14 = Q+Z1
t15 = t14^2
t16 = t15-QQ
t17 = t16-ZZ
Z3 = P*t17