The well established power flow methods -Gauss-Seidel, NewtonRaphson, and the Fast Decoupled Load Wow - are all based on major, classical methodologies of applied mathematics. The Krylov Subspace Power How (KSPF) presented in this paper uses a newer, very successful approach the Krylov subspace methodology - developed in applied linear algebra for the iterative solution of large, sparse systems of linear equations. The method has been adapted to nonlinear equations and used for the solution of the power Bow problem with either an approximation of the Jacobian, as in the Fast Decoupled Load Flow, or in a direct Newton-like manner but without explicitly forming the Jacobian. Convergence rates are from linear to almost quadratic. First the general methodology is described in the paper and then its application to the power flow problem. '1 ne main advantage of KSPF is that no matrix factorizations, only sparse matrix-vector multiplications or evaluations of residuals are used. Preliminary tests suggest that KSPF may become a competitive alternative to existing methods, especially in the case of large systems. © 1995 IEEE.