clear;
clc;
V=3000e3;
r1=30;
r=5000e3;
vb2=11e3;
vb3=33e3;
x=.2;
Xt=.05*r/V;
Xl=r1*r/(vb3^2);
xtotal=(x+Xt+Xl)*%i;
MVA=r*%i*1e-6/xtotal;
Ifault=MVA*1e6/(sqrt(3)*vb3*%i);
Ir=real(Ifault);
Ii=imag(Ifault);
Imod=sqrt((Ir^2)+(Ii^2));
Iangle=atand(Ir/Ii)-90;
Imod=round(Imod);
MVA=round(MVA*10)/10;
mprintf("the value of falut current = %d/_%d Amp \n fault MVA =%f MVA",Imod,Iangle,MVA);