“…It turns then out that these models can generate a viable sequence 5 of radiation dominated, matter dominated and accelerating era matching with the constraints. This has been shown for various parameterizations of the function f (R), most often with some power-law forms (with one, two or three powers of R), but also with square-root, logarithmic and exponential forms for the curvature correction terms [53,50,54,55,56]. As an example we show, in the left panel of Fig.…”