Energy Density Inhomogeneities with Polynomial $f(R)$ Cosmology

In this paper, we study the effects of polynomial $f(R)$ model on the stability of homogeneous energy density in self-gravitating spherical stellar object. For this purpose, we construct couple of evolution equations which relate the Weyl tensor with matter parameters. We explore different factors responsible for density inhomogeneities with non-dissipative dust, isotropic as well as anisotropic fluids and dissipative dust cloud. We find that shear, pressure, dissipative parameters and $f(R)$ terms affect the existence of inhomogeneous energy density.