# Credit: https://stackoverflow.com/a/57579379 import sys eps = 1.0 while eps + 1 > 1: eps /= 2 eps *= 2 sys_eps = sys.float_info.epsilon print(f"The calculated epsilon: {eps}") print(f" The system epsilon: {sys_eps}")