A method to incorporate the steady-state model of the generalized power flow controller into a Newton-Raphson power flow algorithm is disclosed. The disclosed method 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 disclosed method only incorporates the control variables of the shunt voltage sourced converter into the state vector of Newton-Raphson power flow algorithm. The increment of 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.