def Jacobian_re_im(self, pars):
partials = []<EOL>partials.append(self.dre_drho0(pars)[:, np.newaxis])<EOL>partials.append(self.dre_dm(pars))<EOL>partials.append(self.dre_dtau(pars))<EOL>partials.append(self.dre_dc(pars))<EOL>partials.append(self.dim_drho0(pars)[:, np.newaxis])<EOL>partials.append(self.dim_dm(pars))<EOL>partials.append(self.dim_dtau(pars))<EOL>partials.append(self.dim_dc(pars))<EOL>print('<STR_LIT>')<EOL>for x in partials:<EOL><INDENT>print(x.shape)<EOL><DEDENT>J = np.concatenate(partials, axis=<NUM_LIT:1>)<EOL>return J<EOL>
r"""
        :math:`J`

        >>> import sip_models.res.cc as cc
        >>> import numpy as np
        >>> f = np.logspace(-3, 3, 20)
        >>> pars = [100, 0.1, 0.04, 0.8]
        >>> obj = cc.cc(f)
        >>> J = obj.Jacobian_re_im(pars)