Hai... Being a professor and researcher in Power Systems Engineering, I can deliver your project in 3 days using MATLAB code. Basically the convergence characteristics and computational time of load flow methods (GS, NR, DL and FD) are dependent on system size, topology, loading level and contingency. I can give the solution even for a larger system size, which can suit for simulating different operating conditions/contingency and even for research.