We study the shape of the optical luminosity function of Quasi Stellar Objects (QSOs) from the Sloan Digital Sky Survey Data Release Seven (SDSS DR7) over the redshift range 0.3 ≤ ≤ 2.4. By using the Levenberg-Marquardt method of nonlinear least square fit, the observed QSO luminosity function is fitted by a double power-law model with luminosity evolution characterized by a second order polynomial in redshift. For a flat universe with Ω m =0.3 and Ω Λ =0.7, we determine the best-fitting optical luminosity function model parameters.