function lsc5-3-14()
x = [-1 0 1.1 2.4]';
y = [-2 2.4 -4 -1.6]';
A = [2*x 2*y ones(4,1)];
b = x.^2+y.^2;
B = A'*A;
C = inv(B);
z = C*A'*b
R = sqrt(z(3)+z(1)^2+z(2)^2)
theta = linspace(0,360,361)*pi/180;
xx = z(1) + R*cos(theta);
yy = z(2) + R*sin(theta);
plot(x,y)

