code: def calc_bias_correction_abc(acceleration, total_curvature):
      
      inner_arg = 2 * norm.cdf(acceleration) * norm.cdf(-1 * total_curvature)
      bias_correction = norm.ppf(inner_arg)
      return bias_correction
docs: Calculate the bias correction constant for the approximate bootstrap
      confidence (ABC) intervals.
  
      Parameters
      ----------
      acceleration : 1D ndarray of scalars.
          Should contain the ABC intervals' estimated acceleration constants.
      total_curvature : 1D ndarray of scalars.
          Should denote the ABC intervals' computred total curvature values.
  
      Returns
      -------
      bias_correction : 1D ndarray of scalars.
          Contains the computed bias correction for the MLE estimates that the
          ABC interval is being computed for.
  
      References
      ----------
      Efron, Bradley, and Robert J. Tibshirani. An Introduction to the Bootstrap.
          CRC press, 1994. Section 22.6, Equation 22.40, line 1.