%
% Calcola l'epsilon di macchina
%
meps = 1.0;
t=0;
while ((1.0+meps)>1.0)
  meps=meps/2;
  t = t+1;
end
meps = 2*meps
t = t-1

sprintf('eps calcolato:   %.15E\n eps precalcolato: %.15E\n',...
meps,eps)
sprintf('meps-eps = %.15E\n',meps-eps)