def response(self, parameters):
<EOL>self._set_parameters(parameters)<EOL>terms = self.m / (<NUM_LIT:1> + (<NUM_LIT> * self.w * self.tau) ** self.c)<EOL>specs = np.sum(terms, axis=<NUM_LIT:1>)<EOL>ccomplex = self.sigmai * (<NUM_LIT:1> - specs)<EOL>response = sip_response.sip_response(self.f, ccomplex=ccomplex)<EOL>return response<EOL>
r"""Return the forward response in base dimensions
        :math:`\hat{\sigma }(\omega ) = \sigma _\infty \left(1 - \sum_i \frac
        {m_i}{1 + (j \omega \tau_i)^c_i}\right)`

        Parameters
        ----------
        pars:

        Returns
        -------
        response: Nx2 array, first axis denotes frequencies, seconds real and
                  imaginary parts