  In the present paper we study the geodesic motion of test particles and light
rays in the spacetime of a static charged black hole in $f(R)$ gravity. The
complete set of analytic solutions of the geodesic equations in the spacetime
of this black hole are presented. The geodesic equations are solved in terms of
Weierstrass elliptic $\wp$ function and derivatives of Kleinian $\sigma$
function. With the help of parametric diagrams and effective potentials we
analyze the geodesic motion and give a list of all possible orbit types. The
different types of the resulting orbits are characterized in terms of the
conserved energy, angular momentum, charge and cosmological constant.
