We construct an approximate solution to the cosmological perturbation theory around Einstein–de Sitter background up to the fourth-order perturbations. This could be done with the help of the specific symmetry condition imposed on the metric, from which follows that the model density forms an infinite, cubic lattice. To verify the convergence of the perturbative construction, we express the resulting metric as a polynomial in the perturbative parameter and calculate the exact Einstein tensor. In our model, it seems that physical quantities averaged over large scales overlap with the respective Einstein–de Sitter prediction, while local observables could differ significantly from their background counterparts. As an example, we analyze the behavior of the local measurements of the Hubble constant and compare them with the Hubble constant of the homogeneous background model. A difference between these quantities is important in the context of a current Hubble tension problem.