def dre_drho0(self, pars):
self._set_parameters(pars)<EOL>numerator = self.m * self.otc * (np.cos(self.ang) + self.otc)<EOL>term = numerator / self.denom<EOL>specs = np.sum(term, axis=<NUM_LIT:1>)<EOL>result = <NUM_LIT:1> - specs<EOL>return result<EOL>
r""" Compute partial derivative of real parts with respect to
        :math:`\rho_0`

        :math:`\frac{\partial \hat{\rho'}(\omega)}{\partial \rho_0} = 1 -
        \frac{m (\omega \tau)^c cos(\frac{c \pi}{2}) + (\omega \tau)^c}{1 + 2
        (\omega \tau)^c cos(\frac{c \pi}{2}) + (\omega \tau)^{2 c}}`

        Note that partial derivatives towards :math:`\rho_0` are 1D, in
        contrast to the other parameter derivatives, which usually return 2D
        arrays!

        Returns
        -------
        dre_drho0: :class:`numpy.ndarray`
            Size N (nr of frequencies) array with the derivatives