t0 = Z1^2
t1 = Y1^2
t2 = r*t0
t3 = t2-t1
t4 = t3^2
V = s*t4
t5 = Z1^2
t6 = Y1^2
t7 = r*t5
t8 = t7+t6
W = t8^2
Y3 = W-V
Z3 = W+V