A := 0;
B := 2;
C := 1 + (3^(1/2))*I;
E := 1 + tan(Pi/6)*I;x1:= solve (2*(3+d^2)=(3+d)^2,d)[1];
x2:= solve (2*(3+d^2)=(3+d)^2,d)[2];with(geometry):
draw([circ(cA,pA,A,1),circ(cB,pB,B,1),circ(cC,pC,C,1),circ(cE,pE,E,1/x1),circ(cE2,pE2,E,1/abs(x2))]);circ := (cname,pname,z,r) -> circle(cname,[c2p(pname,z),r]):
c2p := (name,z) -> point(name,Re(z),Im(z)):simplify(2*(A^2+ B^2+ C^2+ (E/x)^2) - (A+B+C+E/x)^2);