Standard dynamical system analysis of Einstein-Maxwell equation in $f(R)$ theories is considered in this work. We investigate cosmological dynamics of a uniform magnetic field in the Orthogonal Spatially Homogeneous (OSH) Bianchi type I universe with viable $f(R)$ models of gravity. In this work, the $f(R) = R -\alpha R^n$ and $f(R) = \left( R^b - \Lambda\right)^c$ models are examined by using our dynamical system analysis. Our results show that both of two $f(R)$ models have a viable cosmological consequence identical to the analysis present in Ref.\cite{Amendola:2007nt} for the FLRW background. Contrary to Ref.\cite{Amendola:2007nt}, we discover in our models that there is an additional anisotropic and non-zero cosmological magnetic fields fixed point emerging before the present of the standard matter epoch. This means that the universe has initially isotropic stage with the intermediated epoch as the anisotropic background and it ends up with the isotropic late-time acceleration. The primordial magnetic fields play a crucial role of the shear evolutions obtained from these two models which have the same scaling of the cosmic time as $\sigma\sim t^{-\frac13}$, instead of $\sigma\sim t^{-1}$ for the absence of the primordial magnetic cases.

We consider a dynamical system of phantom scalar field under exponential potential in background of loop quantum cosmology. In our analysis, there is neither stable node nor repeller unstable node but only two saddle points, hence no Big Rip singularity. Physical solutions always possess potential energy greater than magnitude of the negative kinetic energy. We found that the universe bounces after accelerating even in the domination of the phantom field. After bouncing, the universe finally enters oscillatory regime.

In this paper, we investigate and analyze the cosmological dynamics of the universe, with an effect of modified $f(R)$ gravity emerging at cosmological scale. We choose the Einstein frame as a physical frame. We consider phase portraits of the universe at the late time from modified $f(R)$ gravity model. This result gives our universe an acceleration phase expansion without introducing existence of dark energy dominating our universe.