“…The integration is divided in two parts with the first one, during the f (R) era, performed in terms of the cosmic time and the second one, during the mGCG and ΛCDM eras, performed in terms of the scale factor. Following previous works, [31,48,49], we stop the integration when k 2 ≫ z ′′ /z, i.e. when the mode is well inside the horizon, since in this regime X and Y have a highly oscillatory sinusoidal behaviour while |β| 2 is virtually constant.…”