Untitled

function [xc,yc,R] = circfit(x,y)
    A = [x y ones(size(x))];
    b = -(x.^2 + y.^2);
    c = A\b;
    xc = -c(1)/2;
    yc = -c(2)/2;
    R = sqrt((c(1)^2+c(2)^2)/4-c(3));
end




(x-xc)^2 + (y-yc)^2 = R^2

x^2 + y^2 -2xcx - 2ycy + xc^2 + yc^2 - R^2=0

x^2 + y^2 + Ax + By + C = 0

A = -2xc
B = -2yc
C = xc^2 + yc^2 - R^2


A = x1    y1    1
    x2    y2    1
    .     .
    .     .     1

b = -(x^2+y^2)

Ac = b, dove c = [A B C]^T

xc = -A/2, yc = -B/2