% Script machine_eps.m
% find the machine epsilon for base 2

e = 1;
while 1+e > 1
   e = e / 2;
end
e = e * 2;
disp(e);
disp(eps);