def response(self, parameters):
<EOL>self._set_parameters(parameters)<EOL>terms = self.m * (<NUM_LIT:1> - (<NUM_LIT:1> / (<NUM_LIT:1> + (<NUM_LIT> * self.w * self.tau) ** self.c)))<EOL>specs = np.sum(terms, axis=<NUM_LIT:1>)<EOL>rcomplex = self.rho0 * (<NUM_LIT:1> - specs)<EOL>response = sip_response.sip_response(self.f, rcomplex=rcomplex)<EOL>return response<EOL>
r"""Complex response of the Cole-Cole model::
        :math:`\hat{\rho} = \rho_0 \left(1 - \sum_i m_i (1 - \frac{1}{1 + (j
        \omega \tau_i)^c_i})\right)`

        Parameters
        ----------
        parameters: list or tuple or numpy.ndarray
            Cole-Cole model parameters: rho0, m, tau, c (all linear)

        Returns
        -------
        response: :class:`sip_models.sip_response.sip_response`
            model response object