Mesurer l'epsilon machine

eps = 1.0
q = 0
while 1.0 + eps / 2 != 1.0:
    eps /= 2
    q += 1
print(f"Puissance : 2**-{q}, Valeur : {eps}")

Le programme s'arrête pour $q=52$, renvoyant eps $=2^{-52}\approx 2.22\cdot 10^{-16}$. C'est la valeur limite de la précision de la mantisse de 52 bits du format double précision autour de la valeur $1$.