x=rand(1,100)*5
y=rand(1,100)*5
z=rand(1,100)*5

[xc, yc, zc, R] = sphere_fit(x,y,z);

% Draw
[xs, ys, zs] = sphere;

xs = xs * R + xc;
ys = ys * R + yc;
zs = zs * R + zc;

surf(xs, ys, zs);
alpha 0.2
hold
plot3(x,y,z, 'o')