Though an accurate discretization approach for gas flow dynamics, the method of characteristics (MOC) is liable to instability for inappropriate step sizes. This letter addresses the numerical stability limitation of MOC, in the context of IEGS’s optimal scheduling. Specifically, the proposed method enables flexible temporal step sizes without sacrificing accuracy, significantly reducing non-convergence due to numerical oscillations. The effectiveness of the proposed method is validated through case studies in different simulation settings.