Methods and systems are provided for controlling gas lift in at least one production well, which determine a quadratic function that relates flow rate or value of oil produced from a production well to gas flow rate for gas injected into the production well based on flow rate measurements for gas injected into the production well and corresponding flow rate measurements of oil produced from the production well. The quadratic function is used to determine an optimal flow rate for gas injected into the production well. The optimal flow rate for gas injected into the production well is used to control the flow rate of gas injected into the production well, and an oil flow rate produced from the production well is measured at the optimal gas flow rate for iterative processing to refine the quadratic function and determination of optimal gas flow rate if need be.