
function [circleCenter,radius] = ThreePointCircle(obj,x,y,z)
A=[x(1)-y(1),x(2)-y(2);z(1)-y(1),z(2)-y(2)];
B=[x(1)^2-y(1)^2+x(2)^2-y(2)^2;z(1)^2-y(1)^2+z(2)^2-y(2)^2];
ab=A\B;
a=ab(1)/2;
b=ab(2)/2;
circleCenter = [a,b];
c2 = (x(1)-a)^2+(x(2)-b)^2;
radius = sqrt(c2);
end