C1Z1 = C1*Z1 C1D1 = C1*D1 D1Z1 = D1*Z1 D1S1 = D1*S1 U = C1Z1^2 V = C1D1^2 W = D1Z1^2 UV = U-V t0 = C1Z1*D1S1 S3 = 2*t0 t1 = U-W C3 = t1+V D3 = W-UV Z3 = W+UV