Cold Dark Stars made of self-gravitating fermions in the degenerate limit are constructed in General Relativity and in R-square gravity, $f(R)=R+\alpha R^2$. The properties of the resulting Cold Dark Stars in both theories of gravity are studied. It is found that the same gravitational potential is generated for different election of the parameters of the model, such as the mass of the fermion, the self-interacting strength or the value of $\alpha$, thus, a possible confusion in the determination of the dark matter properties and the favored theory of gravity might arise.

In this paper we investigate, from the dynamical systems perspective, the evolution of a Kantowski-Sachs metric in a generic class of $f(R)$ models. We present conditions (i. e., differentiability conditions, existence of minima, monotony intervals, etc.) for a free input function related to the $f(R)$, that guarantee the asymptotic stability of well-motivated physical solutions, specially, self-accelerated solutions, allowing to describe both inflationary- and late-time acceleration stages of the cosmic evolution. We discuss which $f(R)$ theories allows for a cosmic evolution with an acceptable matter era, in correspondence to the modern cosmological paradigm. We find a very rich behavior, and amongst others the universe can result in isotropized solutions with observables in agreement with observations, such as de Sitter, quintessence-like, or phantom solutions. Additionally, we find that a cosmological bounce and turnaround are realized in a part of the parameter-space as a consequence of the metric choice.