def dim_drho0(self, pars):
self._set_parameters(pars)<EOL>result = np.sum(- self.m * self.otc * np.sin(self.ang) / self.denom,<EOL>axis=<NUM_LIT:1>)<EOL>return result<EOL>
r"""
        :math:`\frac{\partial \hat{\rho}''(\omega)}{\partial \rho_0} = -
        \frac{m (\omega \tau)^c sin(\frac{c \pi}{2})}{1 + 2
        (\omega \tau)^c cos(\frac{c \pi}{2}) + (\omega \tau)^{2 c}}`