using this python script could you solve this question in its entirety

import sys

from pathlib import Path

# Add parent directory to Python path

sys$.$path.append$($str$($Path$(\_\_$file$\_\_).$resolve$().$parent.parent$))$

import numpy as np

import equations as eq $$# Import equations module

# Parameters

# Generator G$1$

G$1=[49\mathrm{.}0$e$6,13\mathrm{.}8$e$3,0\mathrm{.}01,0\mathrm{.}11]$# $[$MVA$,$ kV$,$ Rpu, Xpu$]$

# Transformer T$1$

T$1=[47\mathrm{.}0$e$6,13\mathrm{.}8$e$3,120$e$3,0\mathrm{.}01,0\mathrm{.}053]$# $[$MVA$,$ kVpri, kVsec, Rpu, Xpu$]$

# Transformer T$2$

T$2=[47\mathrm{.}0$e$6,120$e$3,14$e$3,0\mathrm{.}011,0\mathrm{.}063]$# $[$MVA$,$ kVpri, kVsec, Rpu, Xpu$]$

# Line L$1$

Zline $=16+75$j $$# Line impedance in ohms

# Load SLoad

Sload $=(43+9$j$)*1$e$6/3$ # Load power in VA

Sbase $=$ G$1[0]/3$

t$1=$ T$1[1]/$T$1[2]$

t$2=$ T$2[1]/$T$2[2]$

#pu table

Va$,$ Ia$,$ Za $=$ eq$.$table$($G$1[1]/$np$.$sqrt$(3),$ Sbase, $1)$

Vb$,$ Ib$,$ Zb $=$ eq$.$table$($Va$,$ Sbase, t$1)$

Vc$,$ Ic$,$ Zc $=$ eq$.$table$($Vb $,$Sbase, t$2)$

#impedance calcs

#Zone A

Zg$1\_$pu $=$ G$1[2]+1$j$*$G$1[3]$

#Zone B

Zt$1\_$pu $=$ eq$.$pu$\_$new$($T$1[3]+1$j$*$T$1[4],$ T$1[2]/$ np$.$sqrt$(3),$ Vb$,$ T$1[0]/3,$ Sbase$)$

Zl$1\_$pu $=($Zline$.$real $+1$j$*$Zline$.$imag$)/$ Zb

Zt$2\_$pu $=$ eq$.$pu$\_$new$($T$2[3]+1$j$*$T$2[4],$ T$2[1]/$ np$.$sqrt$(3),$ Vb$,$ T$2[0]/3,$ Sbase$)$

#Zone C

ZLoad$1\_$pu $=($Vc$/$np$.$conj$($Sload$/$Vc$))/$Zc

#matrix

Zline$1\_$pu $=$ Zt$1\_$pu $+$ Zl$1\_$pu $+$ Zt$2\_$pu

#Question $1$: Bus admittance matrix

Y$11=1/$Zg$1\_$pu $+1/$Zline$1\_$pu

Y$12=-1/$Zline$1\_$pu

Y$21=-1/$Zline$1\_$pu

Y$22=1/$ZLoad$1\_$pu $+1/$Zline$1\_$pu

Ybus $=$ np$.$array$([[$Y$11,$ Y$12],[$Y$21,$ Y$22]])$

#Question $2$: Bus impedance matrix

Zbus $=$ np$.$linalg.inv$($Ybus$)$

Ig$1\_$pu $=1/$ Zg$1\_$pu

Imatrix $=$ np$.$array$([$Ig$1\_$pu$,0])$

Busmatrix $=$ np$.$linalg.solve$($Ybus$,$ Imatrix$)$

#Question $3$: Vb$1$: Voltage at bus $1[$mag$,$ phs$][$kV$]$

Vbus$1=$ Busmatrix$[0]*$ Va

Vb$1=$ eq$.$format$\_$phasor$($eq$.$phasor$($Vbus$1*$np$.$sqrt$(3)/1$e$3))$

#Question $4$: IG$1$: Current from G$1[$mag$,$ phs$][$A$]$

Ig$1=($G$1[1]/$np$.$sqrt$(3)-$ Vbus$1)/($Zg$1\_$pu$*$Za$)$

IG$1=$ eq$.$format$\_$phasor$($eq$.$phasor$($Ig$1))$

#Question $5$: Vb$2$: Voltage at bus $2[$mag$,$ phs$][$kV$]$

Vbus$2=$ Busmatrix$[1]*$ Vc

Vb$2=$ eq$.$format$\_$phasor$($eq$.$phasor$($Vbus$2*$np$.$sqrt$(3)/1$e$3))$

#Question $6$: IL$1$: Current through L$1[$mag$,$ phs$][$A$]$

Iline $=($Vbus$1/$t$1-($Vbus$2*$t$2))/($Zline $+($Zt$1\_$pu$+$Zt$2\_$pu$)*$Zb$)$

IL$1=$ eq$.$format$\_$phasor$($eq$.$phasor$($Iline$))$

#Question $7$: ILoad$1$: Current to Load$1[$mag$,$ phs$][$A$]$

Iload $=$ Vbus$2/($ZLoad$1\_$pu$*$Zc$)$

ILoad$1=$ eq$.$format$\_$phasor$($eq$.$phasor$($Iload$))$

#Question $8$: Pline: Power lost on L$1[$kW$]$

Sline $=3*$ Iline $*$ np$.$conj$($Iline$)*$ Zline

Pline $=$ Sline.real $/1$e$3$

#Question $9$: SLoad$\_$actual: Actual power delivered to load $[$MVA$]$

Sload$\_$actual $=3*$ Vbus$2*$ np$.$conj$($Iload$)/1$e$6$

print$("1$ Bus admittance matrix: $","["+",".$join$([$f$"[\{\text{'},\text{'}.$join$($f$\text{'}\{$x:$\mathrm{.}5\}\text{'}$ for x in row$)\}]"$ for row in Ybus$])+"]")$

print$("2$ Bus impedance matrix: $","["+",".$join$([$f$"[\{\text{'},\text{'}.$join$($f$\text{'}\{$x:$\mathrm{.}5\}\text{'}$ for x in row$)\}]"$ for row in Zbus$])+"]")$

print$($f$"3$ Vb$1$: Voltage at Bus $1[$kV$]$: $\{$Vb$1\}")$

print$($f$"4$ IG$1$: Current from G$1[$A$]$: $\{$IG$1\}")$

print$($f$"5$ Vb$2$: Voltage at Bus $2[$kV$]$: $\{$Vb$2\}")$

print$($f$"6$ IL$1$: Current through L$1[$A$]$: $\{$IL$1\}")$

print$($f$"7$ ILoad$1$: Current to Load$1[$A$]$: $\{$ILoad$1\}")$

print$($f$"8$ Pline: Power lost on L$1[$kW$]$: $\{$Pline:$\mathrm{.}4\}")$

print$($f$"9$ SLoad$\_$actual: Actual power delivered to load $[$MVA$]$: $\{$Sload$\_$actual:$\mathrm{.}4\}")$