A semi-analytical method is presented for the design of gradient index (GRIN) flat lenses. Closed-form expressions are derived to define the refractive index distribution of the lens, for several cases: collimating lenses with on-axis feed, collimating lenses with off-axis feed, lenses converting spherical wavefronts with different wavenumbers, lenses changing the focal number of a quasi-optical system, and Fresnel zone lenses. The design equations are validated by ray-tracing simulations in inhomogeneous media, implemented by numerical solution of the Eikonal equation.