% Matlab's epsilon estimation: eps is the smallest number for which 
%    1<1+eps (in floating point representation)
e=1;
while e+1>1
    e=e/2;
end
fprintf('eps=%6.2e\n',2*e)