def dre_dc(self, pars):
self._set_parameters(pars)<EOL>num1a = np.log(self.w * self.tau) * self.otc * np.sin(self.ang)<EOL>num1b = self.otc * np.cos(self.ang) * np.pi / <NUM_LIT><EOL>term1 = (num1a + num1b) / self.denom<EOL>num2 = self.otc * np.sin(self.c / np.pi) * <NUM_LIT:2><EOL>denom2 = self.denom ** <NUM_LIT:2><EOL>term2 = num2 / denom2<EOL>num3a = <NUM_LIT:2> * np.log(self.w * self.tau) * self.otc * np.cos(self.ang)<EOL>num3b = <NUM_LIT:2> * ((self.w * self.tau) ** <NUM_LIT:2>) * np.pi / <NUM_LIT> * np.sin(self.ang)<EOL>num3c = <NUM_LIT:2> * np.log(self.w * self.tau) * self.otc2<EOL>term3 = num3a - num3b + num3c<EOL>result = self.sigmai * self.m * (term1 + term2 * term3)<EOL>return result<EOL>
r"""
        :math:Add formula