H:=proc(A,B):(A[1]-B[1])^2+(A[2]-B[2])^2: end:
P:=proc(t):[(t+1/t)/2,(t-1/t)/2/I]:end:
T:=proc(A,B,C,D):evalb(expand(((A+B-C-D)^2-4*(A*B+C*D))^2-64*A*B*C*D)=0):end:
U:=proc(A,B,C):H(A,B)*H(B,C)*H(C,A):end:
S:=proc() local K,t,i:K:={seq(P(t[i]),i=1..4)}:T(seq(U(op(K minus {P(t[i])})),i=1..4)):end: print(S()):