clear;
clc;
rb=20e6;
r=10e6;
v1=11e3;
v2=66e3;
x1=5;
X1=x1*rb/r;
xa=20;
xb=20;
xc=20;
xd=20;
xbus=25;
xtr=X1;
xcd=inv(inv(xc)+inv(xd));
xab=inv(inv(xa)+inv(xb));
xcdbus=xcd+xbus;
xn=inv(inv(xab)+inv(xcdbus));
xth=xtr+xn;
mva=rb/xth*100;
i=mva/(1.745*v2);
i=round(i);
printf("the SC MVA=%fMVA \n the SC current=%dA",mva/1e6,i);