Left: comparison between analytical and numerical solutions of the symmetron model for two different redshifts. The density distribution is given by a sphere of constant density. Different colours correspond to different levels of refinement. The thin line is the analytic solution for redshift z = 0 (upper line) and z = 1 (lower line). Right: same for the f(R) solver. See text for details.