SS1 = S1^2
CC1 = 1-SS1
t0 = S1+C1
t1 = t0^2
t2 = t1-SS1
SC1 = t2-CC1
E = D1^2
F = SC1^2
G = a*F
t3 = 4*E
Z3 = t3+G
t4 = 4*E
D3 = t4-G
t5 = CC1^2
t6 = 4*t5
t7 = F+t6
t8 = 2*t7
C3 = t8-Z3
t9 = SC1+D1
t10 = t9^2
t11 = t10-E
t12 = t11-F
S3 = 2*t12