n=200;%stevilo tock
x0=2;%parametri kroznice, ki jih uporabimo, da generiramo podatke
y0=3;
r=4;
faktor=0.1;%faktor, ki doloca sum v podatkih
t=linspace(0,pi,n)';
x=r*cos(t)+x0+randn(n,1)*faktor;%generiramo podatke
y=r*sin(t)+y0+randn(n,1)*faktor;
clf
hold on
plot(x,y,'rx');

A=[2*x 2*y ones(n,1)];%leva stran predolocenega sistema
p=x.^2+y.^2;%desna stran

x=A\p;%resimo sistem z vgrajeno metodo

x0=x(1)%nazaj izracunamo parametre
y0=x(2)
r=sqrt(x(3)+x0^2+y0^2)

t=linspace(0,2*pi,n);
plot(r*cos(t)+x0,r*sin(t)+y0);%narisemo dobljeno kroznico

axis equal