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