Flat, radiation universes with quadratic corrections and asymptotic   analysis
  It was shown long ago by T. V. Ruzmaikina and A. A. Ruzmaikin that within the framework of a homogeneous and isotropic cosmological model quadratic corrections of the gravitational field cannot provide solutions that are both regular initially and go over to Friedmann type at later times. We find here, by applying a dynamical systems approach, the general form of the solution to this class of models in the neighborhood of the initial singularity under the above conditions. 
Spiros Cotsakis and Antonios Tsokaros