% FILE d011hlp.m begins.
x=1; a=3;
disp('    n     x(n)            Error');
for j=1:6
    disp([j-1, x, x-sqrt(3)]);
%    fprintf('   %2d%16.10f   %12.4e\n',j-1,x,x-sqrt(3));
    x=0.5*(x+a/x);
end
% FILE d011hlp.m ends.