Z2sqr = Z2^2
Z1sqr = Z1^2
S1 = Z2sqr * Z2
S2 = Z1sqr * Z1
S1 = Y1 * S1
S2 = Y2 * S2
R = S2 - S1
U1 = X1 * Z2sqr
U2 = X2 * Z1sqr
H = U2 - U1
Rsqr = R^2
Z3 = H * Z1
Hsqr = H^2
Z3 = Z3 * Z2
Hcub = Hsqr * H
U2 = U1 * Hsqr
Hsqr = 2 * U2
X3 = Rsqr - Hsqr
X3 = X3 - Hcub
Y3 = U2 - X3
S2 = S1 * Hcub
Y3 = R * Y3
Y3 = Y3 - S2