  We study the primordial magnetic field generated by the simple model $f^2 FF$
in Starobinsky, $R^2$-inflationary, model. The scale invariant PMF is achieved
at relatively high power index of the coupling function, $\left| \alpha \right|
\approx 7.44$. This model does not suffer from the backreaction problem as long
as, the rate of inflationary expansion, $H$, is in the order of or less than
the upper bound reported by Planck ($\le 3.6 \times 10^{-5} M_\rm{Pl}$) in both
de Sitter and power law expansion, which show similar results. We calculate the
lower limit of the reheating parameter, $R_\rm{rad} > 6.888$ in
$R^2$-inflation. Based on the upper limit obtained from CMB, we find that the
upper limits of magnetic field and reheating energy density as,
$\left(\rho_{B_\rm{end}} \right)_\rm{CMB} < 1.184 \times 10^{-20} M_\rm{Pl}^4$
and $\left(\rho_\rm{reh} \right)_\rm{CMB} < 8.480 \times 10^{-22} M_\rm{Pl}^4$.
All of foregoing results are well more than the lower limit derived from WMAP7
for both large and small field inflation. By using the Planck inflationary
constraints, 2015 in the context of ${R^2}$-inflation, the upper limit of
reheating temperature and energy density for all possible values of $w
_\rm{reh}$ are respectively constrained as, $T_\rm{reh} < 4.32 \times 10^{13}
\rm{GeV}$ and $\rho_\rm{reh} < 3.259 \times 10^{-18} M_\rm{Pl}^4$ at $n_\rm{s}
\approx 0.9674$. This value of spectral index is well consistent with Planck,
2015 results. Adopting $T_\rm{reh}$, enables us to constrain the reheating
e-folds number, $N_\rm{reh}$ on the range $1 < N_\rm{reh} < 8.3$, for $- 1/3 <
w_\rm{reh} < 1$. By using the scale invariant PMF generated by $f^2 FF$, we
find that the upper limit of present magnetic field, $B_0 < 8.058 \times
10^{-9} \rm{G}$.
