We propose several extensions of the Starobinsky model of inflation, which obey all observational constraints on the inflationary parameters, by demanding that both the inflaton scalar potential in the Einstein frame and the F (R) gravity function in the Jordan frame have the explicit dependence upon fields and parameters in terms of elementary functions. All our models are continuously connected to the original Starobinsky model via changing the parameters. We modify the Starobinsky (R + R 2 ) model by adding an R 3 -term, an R 4term, and an R 3/2 -term, respectively, and calculate the scalar potentials, the inflationary observables and the allowed limits on the deformation parameters in all these cases. We also deform the scalar potential of the Starobinsky model in the Einstein frame, in powers of y = exp − 2 3 φ/M P l , where φ is the canonical inflaton (scalaron) field, calculate the corresponding F (R) gravity functions in the two new cases, and find the restrictions on the deformation parameters in the lowest orders with respect to the variable y that is physically small during slow-roll inflation.