Clifford:=proc()local a,b,C,i,j,P,k:
for i from 1 to 4 do
C[i]:=x^2+y^2+a[i]*x+b[i]*y:
od:
for i from 1 to 4 do
for j from i+1 to 4 do
P[i,j]:=
OtherCeCe
(C[i],C[j],[0,0]):
od:
od:
for i from 1 to 4 do
for j from i+1 to 4 do
for k from j+1 to 4 do
C[i,j,k]:=
Ce
(P[i,j],P[i,k],P[j,k]):
od:od:od:
ConcurrentCes
(C[1,2,3],C[1,2,4],C[1,3,4],C[2,3,4]):
end: