format long e
j=0
while c(1+2^(-j))>1
    [j c(1+2^(-j))]
    j=j+1
end
j
epsilon_de_maquina= eps
2^(-j)