x = (−1)s × 1.f × 2e;   ulp(x) = 2e−p   (normal, 2e ≤ |x| < 2e+1);   ε32 = 2−23, ε64 = 2−52