code: def cdf_abramowitz_stegun(x):
      
      if x >= 0:
          # if x > 7.0:
          #    return 1 - ONE_OVER_SQRT_OF_TWO_PI * exp(-0.5*x*x)/sqrt(1.0+x*x)
          result = 1.0 / (1.0 + 0.2316419 * x)
          ret = 1.0 - ONE_OVER_SQRT_OF_TWO_PI * exp(-0.5 * x * x) * (
              result * (0.31938153 +
                        result * (-0.356563782 +
                                  result * (1.781477937 +
                                            result * (-1.821255978 +
                                                      result * 1.330274429)))))
          return ret
  
      if x <= 0:
          # if x < -7.0:
          #    return 1 - one_over_sqrt_of_two_pi() * exp(-0.5*x*x)/sqrt(1.0+x*x)
          result = 1.0 / (1.0 - 0.2316419 * x)
          ret = ONE_OVER_SQRT_OF_TWO_PI * exp(-0.5 * x * x) * (
              result * (0.31938153 +
                        result * (-0.356563782 +
                                  result * (1.781477937 +
                                            result * (-1.821255978 +
                                                      result * 1.330274429)))))
  
          return ret
docs: The cumulative distribution function of the standard normal distribution.
      The standard implementation, following Abramowitz/Stegun, (26.2.17).