Develop a Matlab code that solves the power flow problem in Lecture $6-7$ for both PQ and PV buses and the slack bus using the Gauss$-$Seidel method. You are expected to develop a code with a loop and obtain the converged results. Please create aRun ModeRun Mode

Solution Animation Running

AC

Viewing Present

Solution Animation Running

AC

Viewing Present matlab code and a document recording your results $($P$,$ Q$,$ V$,$ delta for each bus$).$ Your results should match with the results shown in PowerWorld.

Results after the $1$st iteration:

P$1=2\mathrm{.}3792$

Q$1=-0\mathrm{.}4256$

v$2$abs $=0\mathrm{.}9613$

v$2$angle $=-16\mathrm{.}5430$

v$3$abs $=1\mathrm{.}0500$

v$3$angle $=2\mathrm{.}0766$

v$4$abs $=1\mathrm{.}0354$

v$4$angle $=0\mathrm{.}3529$

v$5$abs $=1\mathrm{.}0061$

v$5$angle $=-2\mathrm{.}7473$

You should get the results for these variables after they converge. The convergence criteria is the changes of the values for v$2,$ v$4,$ and v$5$ between two iterations is less than $1$e$-5.$

Your correct answer should be very close to the numbers shown in the following figure: