clear;
clc;
r=6e6;
v1=11e3;
v2=66e3;
xg=.1;
xt=.09;
z=4+(1*%i);
zb=v2^2/r;
zpu=z/zb;
E=1;
Ifault=E/(zpu+((xg+xt)*%i));
Ir=real(Ifault);
Ii=imag(Ifault);
Imod=sqrt((Ir^2)+(Ii^2));
Ib=r/(sqrt(3)*v2);
i=Imod*Ib;
igb=r/(sqrt(3)*v1);
ig=igb*Imod;
i=fix(i);
ig=fix(ig);
mprintf("the base current on HT side = %dA\n the current from generator = %dA",i,ig);