Working in Einstein frame, we introduce, in order to avoid singularities, holonomy corrections to the f(R)=R+aR2 model. We perform a detailed analytical and numerical study when holonomy corrections are taken into account ...