Global stability analysis of the WASP-148 planetary system. We fixed all orbital parameters of the stable MAP solution (Table 4), and we varied the inclination, i_c, and the longitude of the ascending node, Ω_c, of planet-c. The step size is 0.25° in both axes. For each initial condition, the system was integrated over 50 kyr and a stability criterion was derived with the frequency analysis of the mean longitude. White dashed curves give the isolines of constant mutual inclination I = 10^°, 20°, 30°, and 40°. The white dot marks the position of the stable MAP solution from Table 4. The colour bar corresponds to the one in Fig. 11.