
u = 1;
for i=0:999
  if (1+u == 1)
    break;
  end
  fprintf(1,' %4d u = %20.15e\n',i,u); 
  u=u/2;
end
