docs: def dre_dm(self, pars):
          r"""
          :math:Add formula
          """
          self._set_parameters(pars)
          terms = self.num / self.denom
          result = - self.sigmai * terms
  
          return result