A method to incorporate the steady-state model of the generalized power flow controller into a Newton-Raphson power flow algorithm adopts a flexible steady-state model of the generalized power flow controller, which can be applied to calculate the power flow solution of a power grid embedded with STATCOM, UPFC, GUPFC and the generalized power flow controller in a single framework. The method only incorporates the control variables of the shunt voltage sourced converter into the state vector of the Newton-Raphson power flow algorithm. The increment of the state variables due to incorporating the generalized power flow controller is less than the prior art. Further, the method can preserve the quadratic convergence characteristic of the Newton-Raphson power flow algorithm after embedding the generalized power flow controller into a power grid.