function [x,y,r]=cercle3pts(P1,P2,P3)
[a,b,c]=eqdroite(P1,P2);
[d,e,f]=eqdroite(P2,P3);
delta=a*e-d*b;
x=(c*e-b*f)/delta;
y=(a*f-d*c)/delta;
r=sqrt((P1(1)-x)^2+(P1(2)-y)^2);