code: def dre_dsigmai(self, pars):
          r
          self._set_parameters(pars)
          terms = self.m * self.num / self.denom
          specs = np.sum(terms, axis=1)
          result = 1 - specs
  
          return result
docs: r"""
          :math:Add formula