Load flow is an important tool for studying, designing, and analyzing power systems. It allows power system engineers to determine whether the operation and configuration of the power system is safe under varying loading conditions. It is necessary to model and simulate such a system in order to determine the power flow and losses. This research paper focuses on using numerical methods such as Newton Raphson and Gauss Seidel power flow equations for load flow analysis to calculate bus voltage magnitudes, phase angles, real and reactive power of each bus of an IEEE 9-bus test system. Newton Raphson's computation offers fast, accurate convergence but demands complex implementation, whereas Gauss Siedel is simpler but converges slower with lower accuracy. The analysis was carried out using a MATLAB program. By manipulating variables such as power injections, voltage magnitudes, and phase angles, it solves nonlinear equations iteratively to establish stable operating points which aids in enhancing power system analysis. The line losses for the two methods are compared and the system's total load and generation power are also displayed. The consideration of line losses and assessment of total load generation is crucial for maintaining system efficiency, reliability and preventing voltage instability and equipment damage. The results are also used to generate a directed graph which shows the interconnected nature of the power system, aiding engineers in understanding power flow paths, identifying potential issues, and making informed decisions about system operations. The Newton Raphson method yields the lowest loss, with 4.585MW and 10.789Mvar. In comparison, the Gauss Seidel method achieved 4.809MW and 10.798Mvar.