docs: def L(self, x, tau):
          """
          Logarithm that appears frequently
          :param x: r/Rs
          :param tau: t/Rs
          :return:
          """
  
          return np.log(x * (tau + np.sqrt(tau ** 2 + x ** 2)) ** -1)