Previously, the AC-case negative bias temperature instability has been modelled with iterative method to account for its recovery effect. A significant disadvantage of the iterative method is its low time-efficiency. A new algorithm is proposed, which requires a one-time simulation; besides, it gives identical result as the iterative method. The time-efficiency could thus be improved by orders.