i=5;
r=.1;
va=i**2*r;
j=10+2*va;
print"the Ct of VA and A may be used",j,i
iS=5.;
pr=2.;
ir=2.5;
pe=pr*(iS/ir)**2
print"the burden on transformer Pe=VA",pe
from math import sqrt
ct=2000./5.;
i=40e3;
r1=.31;
a=28.45e-4;
r2=2.;
iS=i/ct;
e=iS*(r1+r2);
f=50.;
B=e/(4.4*f*ct*a);
C=B/sqrt(2);
C=round(C*10)/10;
print"saturation magnetic field max=Wb\t rms value=Wb",B,C
r1=.1;
r2=.4;
r=r1+r2;
i=1e3/10;
ip=100*5/50;
ie=10.;
e=45.;
y=i-ie;
per=-(ie*y-(10*i))/(i*10);
print"the percentage R.E at 1000A =percent",per*100