def dim_dc(self, pars):
self._set_parameters(pars)<EOL>nom1a = - self.m * np.log(self.w * self.tau) * self.otc *np.sin(self.ang)<EOL>nom1b = - self.m * self.otc * (np.pi / <NUM_LIT>) * np.cos(self.ang)<EOL>term1 = (nom1a + nom1b) / self.denom<EOL>nom2 = (self.m * self.otc * np.sin(self.ang)) *(<NUM_LIT:2> * np.log(self.w * self.tau) * self.otc * np.cos(self.ang) -<EOL><NUM_LIT:2> * self.otc * (np.pi / <NUM_LIT>) * np.sin(self.ang) +<EOL><NUM_LIT:2> * np.log(self.w * self.tau) * self.otc2)<EOL>term2 = nom2 / self.denom ** <NUM_LIT:2><EOL>result = term1 + term2<EOL>result *= self.rho0<EOL>return result<EOL>
r"""
        :math:`\frac{\partial \hat{\rho''}(\omega)}{\partial c} = \rho_0
        \frac{-m sin(\frac{c \pi}{2}) ln(\omega \tau)(\omega \tau)^c - m
        (\omega \tau)^c \frac{\pi}{2} cos(\frac{\pi}{2}}{1 + 2 (\omega \tau)^c
        cos(\frac{c \pi}{2}) + (\omega \tau)^{2 c}} + \rho_0 \frac{\left[-m
        (\omega \tau)^c cos(\frac{c \pi}{2}) \right] \cdot \left[ -2 ln(\omega
        \tau) (\omega \tau)^c cos(\frac{c \pi}{2}) + 2 (\omega \tau)^c
        \frac{\pi}{2} cos(\frac{c \pi}{2}) \right] + \left[2 ln(\omega \tau)
        (\omega \tau)^{2 c}\right]}{\left[1 + 2 (\omega \tau)^c cos(\frac{c
        \pi}{2}) + (\omega \tau)^{2 c}\right]^2}`