epsilon = 1;
it = 0;
maxit = 10000;
while it < maxit
    epsilon = epsilon/2;
    b = 1 + epsilon;
    if b==1
        break;
    end
    it=it+1;
end