code: def atanh(x):
      
      if isinstance(x, UncertainFunction):
          mcpts = np.arctanh(x._mcpts)
          return UncertainFunction(mcpts)
      else:
          return np.arctanh(x)
docs: Inverse hyperbolic tangent