This paper describes a critical evaluation of Non Divergent Load flow methods in well, ill and unsolvable conditioned systems. The comparison studies deals Multiple Load Flow Solution (MLFS) based Second-Order Load-Flow (SOLF)in polar coordinate and Continuation Load Flow (CLF). The analytical bases, ability consideration of theses methods to return operation of power system from unsolvable to solvable region solution. Attention is given to the problems and techniques of to provide optimal recommendations of the parameters, that using in these Non-Divergent Load flow methods to lead to solvable region, based on inequality constraints of power system. A part of the survey, this paper also presents the comparison of numerical result using different type of aforesaid load flow methods for well and iII-conditioned systems. Accordingly, load flow simulation has been solved using the C++ programming.A t first of 70's, ability of load flow methods were considered to over come, outage security assessment, optimization and stability [1 ], [31]. Thus, the load-flow studies required the less computer memory and far less execution time. These issues were addressed in the early 1970's when Stott and Alsac presented Fast Decoupled Load Flow (FDLF) model [I]. The FDLF was the usual choice in transmission applications. But, during 70's power system operation was confronted to ill conditioning. The ill conditioning for FDLF was high ratios of lines rlx, connections of very low and very high impedance lines at a bus that could cause FDLF to converge slowly or was divergent [2]. However, when reliability and accuracy, rather than speed of response, was a concern, or when the decoupling principle did not hold, the Newton-Raphson (NR) was the preferred .But conventional NR was poor convergence for radial distribution systems. There were heavy loading at some buses and results in low voltages at these [12], [36]. On other meaning that as the system loading approaches critical loading the Jacobian matrix tends to become singular[3], [39]. Under this condition, even the NR was encountered difficulties to reach a solution. Because, efficient sparsity oriented implementation of NR decreases and no solution from initial estimate increases [7],[17]. This issue has 978-1-4244-5940-711 0/$26.00©20 10 IEEE motivated the development of alternative methodologies, based on the NR iterative scheme, specifically [34]. At end of 70's, second-order load-flow (SOLF) methods began to appear [5], [20]. First second order load fl ow technique based on the Taylor series expansion of a load flow equation was in polar coordinate form. In many cases, this second order required lesser iterations, had better convergence characteristics in almost all the load flows attempted by the authors than NR technique. Moreover, it had also been shown that the elements of the second order coefficient matrix need not be stored separately. Rectangular forms of second order method as a fast load flow method retaining nonlinearly (Iwamoto's Method) was introduced in 1978 [5],...