Black hole solutions in F(R) gravity with conformal anomaly

In this paper, we consider $F(R)=R+f(R)$ theory instead of Einstein gravity with conformal anomaly and look for its analytical solutions. Depending on the free parameters, one may obtain both uncharged and charged solutions for some classes of $F(R)$ models. Calculation of Kretschmann scalar shows that there is a singularity located at $r=0$, which the geometry of uncharged (charged) solution is corresponding to the Schwarzschild (Reissner-Nordstr\"om) singularity. Further, we discuss the viability of our models in details. We show that these models can be stable depending on their parameters and in different epoches of the universe.