Usually Newton method is applied to solve these equations, which requires the solution of a linear system at each iteration.
The formulation of the power flow problem is modified by increasing the number of variables to be computed, such that an under-determined linear system is solved at each iteration.