#calculate Reactance of coil
#Given data :
import math
f=50.##Supply frequency in Hz
C=4.5*10.**-6##in Farad
Omega_L=1./3./2./math.pi/f/C##in ohm
print '%s %.2f' %("Reactance of coil (ohm) :",Omega_L)#
#calculate Necessary Inductance of peterson coil,Rating of supressor coil in MVA
#Given data :
import math
V=132.*1000.##V
f=50.##Hz
r=10./1000.##m
d1=4.##m
d2=4.##m
d3=d1+d2##m
epsilon_o=8.854*10.**-12##constant
l_tl=192.*1000.##length of transmission line in m
C=2.*math.pi*epsilon_o/math.log((d1*d2*d3)**(1./3.)/r)*l_tl##in Farad
L=1./3./(2.*math.pi*f)**2./C##H
print '%s %.3f' %("Necessary Inductance of peterson coil in H : ",L)#
VP=V/math.sqrt(3.)##V
IL=VP/(2.*math.pi*f)/L##A
Rating=VP*IL/1000.##kVA
print '%s %.2f' %("Rating of supressor coil in MVA :",Rating/1000)#