docs: def ad_unif_inf(statistic):
      """
      Approximates the limiting distribution to about 5 decimal digits.
      """
      z = statistic
      if z < 2:
          return (exp(-1.2337141 / z) / sqrt(z) *
                                  (2.00012 + (.247105 - (.0649821 - (.0347962 -
                                  (.011672 - .00168691 * z) * z) * z) * z) * z))
      else:
          return exp(-exp(1.0776 - (2.30695 - (.43424 - (.082433 -
                                  (.008056 - .0003146 * z) * z) * z) * z) * z))