Regular black holes in $f(R)$ gravity coupled to nonlinear electrodynamics

We obtain a class of regular black hole solutions in four-dimensional $f(R)$ gravity, $R$ being the curvature scalar, coupled to a nonlinear electromagnetic source. The metric formalism is used and static spherically symmetric spacetimes are assumed. The resulting $f(R)$ and nonlinear electrodynamics functions are characterized by a one-parameter family of solutions which are generalizations of known regular black holes in general relativity coupled to nonlinear electrodynamics. The related regular black holes of general relativity are recovered when the free parameter vanishes, in which case one has $f(R)\propto R$. We analyze the regularity of the solutions and also show that there are particular solutions that violate only the strong energy condition