clear
close all

a=1.0;
while (1+a ~= 1)
  e=a;
  a=a/2;
end

fprintf('double eps = %e\n',e)


a=single(1.0);
while (1+a ~= 1)
  e=a;
  a=a/2;
end

fprintf('single eps = %e\n',e)

%>> a4_geteps
%double eps = 2.220446e-16
%single eps = 1.192093e-07

