Given z_1 we can calculate z_2 according to
4z_2^2+z_1^2-2iz_1z_2=0 or
z_2 =i/4(1pm sqrt(5)) z_1
which means that given z_1 to obtain z_2 we need two transformations
1) Scaling by 1/4(1pm sqrt5)
2) Rotating pi/2 counterclockwise as associated to the product by i
so with z_0 = 0 we can build two triangles
[z_0,z_1,z_2^a] and [z_0,z_1,z_2^b]
with
z_2^a=i/4(1+ sqrt(5)) z_1
z_2^b=i/4(1- sqrt(5)) z_1
so z_2^a,z_0,z_2^b are aligned and
[z_2^a,z_2^b] is perpendicular to [z_0,z_1]
and given (z_1)_k all triangles [z_2^a,z_1,z_2^b]_k are similar.
The least angles at [z_2^a,z_0,z_1] and [z_2^b,z_0,z_1] are respectively
alpha = arctan((1-sqrt(5))/4), beta = arctan((1+sqrt(5))/4)
and
cot(alpha)+cot(beta) = 4/(1 - sqrt[5]) + 4/(1 + sqrt[5])=-2