The cosmological dynamics for the model f(R) = R − μRc [1 − exp (− 2R∕Rc)] was studied in [396].