Abstract In this study, we analyze f (R) gravity employing the f (R) = R+ ₁ Rᵖ - ₂ R^-q f (R) = R + α 1 R p - α 2 R - q model within the Friedmann–Lemaître–Robertson–Walker (FLRW) background. We derive the Friedmann equations via modified gravity action and subsequently reexpress them in terms of the standard Friedmann equations. Our investigation explores the behavior of a bouncing cosmological model within this modified gravity context, offering a potential solution to the singularity problem encountered in the standard Big Bang cosmology. We represent cosmological parameters as functions of cosmic time and scrutinize the conditions for a cosmic bounce. Additionally, we reconstruct the f (R) gravity model using the redshift parameter and graphically present cosmological parameters as functions of redshift. These plots reveal an accelerated nature of Universe expansion. Furthermore, we reconstruct f (R) gravity models for scale factors, a (t) = e^e^{ t} a (t) = e e λ t which exhibits the bouncing behavior. Finally, we investigated the stability scenario by examining the sound speed. The resulting graph of sound speed plotted against redshift confirms late-time stability.
Ilyas et al. (Fri,) studied this question.