source 2007 Duquesne, page 103, Table 1, plus V=R/2
parameter b
assume b = -2 a
parameter half
assume half 2=1
compute A1 = XX1
compute A2 = XX2
compute A3 = half R1
compute A4 = half R2
compute A5 = ZZ1
compute A6 = ZZ2
compute A7 = Y1
compute A8 = Y2
compute A9 = A7 A8
compute A7 = A7+A3
compute A8 = A8+A4
compute A3 = A3 A4
compute A7 = A7 A8
compute A7 = A7-A9
compute A7 = A7-A3
compute A4 = A1 A2
compute A8 = A5 A6
compute A1 = A1+A5
compute A2 = A2+A6
compute A5 = A1 A2
compute A5 = A5-A4
compute A5 = A5-A8
compute A1 = A8-A4
compute A2 = A8+A4
compute A6 = b A3
compute A6 = A9-A6
compute A6 = A6 A2
compute A3 = 2 A3
compute A3 = A5 A3
compute A8 = A6+A3
compute A2 = A7^2
compute A4 = A1 A7
compute A6 = A1^2
compute X3 = A7
compute Y3 = A8
compute Z3 = A1
compute XX3 = A2
compute R3 = 2 A4
compute ZZ3 = A6