clear, clc
x=[0 4 6.5]';
y=[1.42 6.18 4.75]';

A=[ones(size(x)) x y];
b=[x.^2+y.^2];
c=A\b;
x0=c(2)/2
y0=c(3)/2
R=sqrt(c(1)+x0^2+y0^2)
fi=0:pi/50:2*pi;
plotx=x0+R*cos(fi);
ploty=y0+R*sin(fi);
plot(plotx,ploty,x,y,'*')
axis equal