docs: def _get_taulny(self, C, mag):
          """
          Returns the inter-event random effects coefficient (tau)
          Equation 28.
          """
          if mag <= 4.5:
              return C["tau1"]
          elif mag >= 5.5:
              return C["tau2"]
          else:
              return C["tau2"] + (C["tau1"] - C["tau2"]) * (5.5 - mag)