Stability analysis of the nominal fit (Table 2) of the GJ 163 planetary system. For fixed initial conditions, the phase space of the system is explored by varying the semi-major axis a_p and eccentricity e_p of each planet, b, c, and d, respectively. The step size is 10^-5 AU in semi-major axis and 10^-2 in eccentricity. For each initial condition, the system is integrated over 200 yr and a stability criterion is derived with the frequency analysis of the mean longitude (Laskar 1990, 1993). As in Correia et al. (2005, 2009, 2010), the chaotic diffusion is measured by the variation in the frequencies. The red zone corresponds to highly unstable orbits, while the dark blue region can be assumed to be stable on a billion-year timescale. The contour curves indicate the value of χ² obtained for each choice of parameters.