A new method for finding the optimum gas injection rates for a group of continuous gas lift wells to maximize the total oil production rate is established. The new method uses a quasi-Newton nonlinear optimization technique which is incorporated with the gradient projection method. The method is capable of accommodating restrictions to the gas injection rates. The only requirement for fast convergence is that a reasonable estimate of the gas injection rates must be supplied as an initial point to the optimization method. A method of estimating the gas injection rates is developed for that purpose. A computer program is developed capable of implementing the new optimization method as well as generating the initial estimate of the gas injection rates. This program is then successfully tested on field data under both unlimited and limited gas supply. The new optimization technique demonstrates superior performance, faster convergence, and greater application.