Higher derivative corrections of f(R) gravity with varying equation of state in the case of variable G and Λ

In this paper, we study three models of f(R) modified gravity including higher order terms based on different equation of state parameters. We also assume variable G and Λ. Using numerical analysis, we obtain the behavior of some important cosmological parameters, like the Hubble expansion H parameter and the deceleration parameter q. Statefinder diagnostics are also performed for all models. (author)