a0 = S1 a1 = C1 a2 = D1 m = a1^2 l2 = a0*a2 n = l2^2 t0 = a1*a2 t1 = t0^2 l3 = 2*t1 r3 = m+n t2 = a1+l2 t3 = t2^2 r0 = t3-r3 r1 = m-n r2 = l3-r1 S3 = r0 C3 = r1 D3 = r2 Z3 = r3