Concept:
Bus Admittance Matrix (YBus):
Consider a small power system network consisting of two generating stations, three transmission lines, one load, and a static capacitor connected to load bus 3.
Assumed that the network is symmetrical and operating under the balance conditions.

The node voltage equation of the system can be written as,
I1=(y12+y31)V1−y12V2−y31V3
I2=−y12V1+(y12+y23)V2−y23V3
−I3=−y31V1−y23V2+(y31+y23+y30)V3
Where,
y12=z121,;y23=z231,;y31=z311
The above equation can be written in form of a matrix,
\(\left[ {\begin{array}{{20}{c}} {{I_1}}\ {{I_2}}\ { - {I_3}} \end{array}} \right] = \left[ {\begin{array}{{20}{c}} {{y_{12}} + {y_{31}}}&{ - {y_{12}}}&{ - {y_{31}}}\ { - {y_{12}}}&{{y_{12}} + {y_{23}}}&{ - {y_{23}}}\ { - {y_{31}}}&{ - {y_{23}}}&{{y_{31}} + {y_{23}} + {y_{30}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{V_1}}\ {{V_2}}\ {{V_3}} \end{array}} \right]\)
\(\left[ {\begin{array}{{20}{c}} {{I_1}}\ {{I_2}}\ { - {I_3}} \end{array}} \right] = \left[ {\begin{array}{{20}{c}} {{Y_{11}}}&{{Y_{12}}}&{{Y_{13}}}\ {{Y_{21}}}&{{Y_{22}}}&{{Y_{23}}}\ {{Y_{31}}}&{{Y_{32}}}&{{Y_{33}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{V_1}}\ {{V_2}}\ {{V_3}} \end{array}} \right]\)
Where,
Y11 = y12 + y31, Y22 = y21 + y23, Y33 = y30 + y31 + y32
Y12 = Y21 = - y12, Y23 = Y32 = -y23, Y31 = Y13 = -y31
The element Y11, Y22, Y33 is termed as self admittances.
The elements Y12, Y13, Y21, Y23, Y31, Y32 termed as mutual admittance.
Calculation:
Given circuit can be drawn as,

Applying KCL at node V,
I1 + I2 + I3 + I4 = 0
j1E1−V+j1E2−V+j1E3−V+j1−V=0
4V = E1 + E2 + E3
I1=j1E1−V=j1E1−(4E1+E2+E3)
I1=j4−3E1+j41E2+j41E3
I2=j1E2−V=j1E2−(4E1+E2+E3)
I2=j41E1+j4−3E2+j41E3
I3=j1E3−V=j1E3−(4E1+E2+E3)
I1=j41E1+j41E2+j4−3E3
The equation of current can be written in form of the matrix,
\(\left[ {\begin{array}{{20}{c}} {{I_1}}\ {{I_2}}\ {{I_3}} \end{array}} \right] = j\left[ {\begin{array}{{20}{c}} {\frac{{ - 3}}{4}}&{\frac{1}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{{ - 3}}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{1}{4}}&{\frac{{ - 3}}{4}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{E_1}}\ {{E_2}}\ {{E_3}} \end{array}} \right]\)
\({Y_{bus}} = j\left[ {\begin{array}{*{20}{c}} {\frac{{ - 3}}{4}}&{\frac{1}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{{ - 3}}{4}}&{\frac{1}{4}}\ {\frac{1}{4}}&{\frac{1}{4}}&{\frac{{ - 3}}{4}} \end{array}} \right]\)
\(Y_{bus}=\left[ \begin{array}{20 {c*}} {-\frac 3 4 j}&{\frac 1 4 j}&{\frac 1 4 j}\\ {\frac 1 4 j}&{-\frac 3 4 j}&{\frac 1 4j}\\ {\frac 1 4 j}&{\frac 1 4 j}&{-\frac 3 4 j} \end{array} \right]\)