def dim_dc(self, pars):
self._set_parameters(pars)<EOL>num1a = self.m * np.sin(self.ang) * np.log(self.w * self.tau)* self.otc<EOL>num1b = self.m * self.otc * np.pi / <NUM_LIT:2> * np.cos(np.pi / <NUM_LIT:2>)<EOL>term1 = self.sigma0 * (-num1a - num1b) / self.denom<EOL>num2a = -self.m * self.otc * np.cos(self.ang)<EOL>num2b = -<NUM_LIT:2> * np.log(self.w * self.tau) * self.otc * np.cos(self.ang)<EOL>num2c = <NUM_LIT:2> * self.otc * np.pi / <NUM_LIT:2> * np.cos(self.ang)<EOL>num2d = <NUM_LIT:2> * np.log(self.w * self.tau) * self.otc2<EOL>numerator = num2a * (num2b + num2c) + num2d<EOL>term2 = self.sigma0 * numerator / (self.denom ** <NUM_LIT:2>)<EOL>result = term1 + term2<EOL>return result<EOL>
r"""
        :math:Add formula