  In this work, we analyze the wormhole solutions in $f(R)$ gravity.
Specifically we sought for wormhole geometry solutions for the following three
shape functions: (i)
$b(r)=r_{0}+\rho_{0}r_{0}^{3}\ln\left(\frac{r_{0}}{r}\right)$, (ii)
$b(r)=r_{0}+\gamma r_{0}\left(1-\dfrac{r_{0}}{r}\right)$, and (iii)
$b(r)=\alpha +\beta r$, under some legitimate physical conditions on the
parameters as well as constants involved here with in the shape functions. It
is observed from the graphical plots that the behaviour of the physical
parameters are interesting and viable.
