This paper is concerned with the problem of finding all the real solutions (all components of the solution vector must be real values) of load flow equations. Solutions in which some of the components are complex values are of no interest as they have no physical significance as a load flow solution. This problem issignificant not only because of its theoretical challenge but also, its relationship with several system behavior related issues. Approaches suggested so far for solving this problem are rather ad hoc, computationally demanding and have been demonstrated only on very small systems. Further, it has been subsequently shown by others that many of these methods are not capable of finding all solutions. In this work a new approach is proposed which is more systematic and seems to have the potential to handle even large problems. We show that for any system it is possible to find the multiple load flow solutions (MLFS) corresponding to a given operating point extremely easily, starting from a set of points that are referred to as zero load solutions (ZLS) in this paper.It is shown that the complete set of ZLS is unique for a system and MLFS for any other operating point can be obtained starting from these ZLS using only the Newton’s load flow method. The set of procedures for implementing the proposed scheme are illustrated and their features are highlighted by considering several sample systems.