# Chapter 17, Electromagnetic waves¶

## Example 1, page 550¶

In :
from __future__ import division
from numpy import pi
# magnitude
#given data :
R=7*10**8 # in m
P=3.8*10**26 # in Watt
S=P/(4*pi*R**2)
print "Magnitude of poynting vector, S = %0.3e W/m^2 " %S

Magnitude of poynting vector, S = 6.171e+07 W/m^2


## Example 2, page 551¶

In :
from numpy import pi
# Poynting vector
#given data :
R=1.5*10**11 # in m
P=3.8*10**26 # in Watt
S=P/(4*pi*R**2) # in W/m**2
Se=round(S*60/(4.2*10**4))
print "Poynting vector, Se = %0.2f cal/cm^2-m " %Se

Poynting vector, Se = 2.00 cal/cm^2-m


## Example 3, page 560¶

In :
from numpy import sqrt
# Amplitude and magnetic field
#given data :
S=2 # in cal/cm**2- min
EH=S*4.2*10**4/60 # joule/m**2 sec
mu0=4*pi*10**-7
epsilon0=8.85*10**-12
EbyH=sqrt(mu0/epsilon0)
E=sqrt(EH*EbyH)
H=EH/E
E0=E*sqrt(2)
H0=H*sqrt(2)
print "E = %0.2f V/m "%E
print "H  = %0.3f Amp-turn/m "%H
print "Amplitude of electric fields of radiation, E0 = %0.f V/m  " %E0
print "Magnetice field of radition, H0 = %0.2f Amp-turn/m " %H0

E = 726.32 V/m
H  = 1.928 Amp-turn/m
Amplitude of electric fields of radiation, E0 = 1027 V/m
Magnetice field of radition, H0 = 2.73 Amp-turn/m


## Example 4, page 560¶

In :
from numpy import pi
# electric and magnetic field
#given data :
r=2 # in m
mu0=4*pi*10**-7
epsilon0=8.85*10**-12
EbyH=sqrt(mu0/epsilon0)
EH=1000/(4*r**2*pi**2) # in W/m**2
E=sqrt(EH*EbyH)
H=(EH/E)
print "Intensities of electric, E = %0.2f V/m" %E
print "Magnetic field of radiation, H = %0.4f Amp-turn/m  " %H

Intensities of electric, E = 48.85 V/m
Magnetic field of radiation, H = 0.1296 Amp-turn/m


## Example 5, page 593¶

In :
from math import degrees, pi, asin, sin, tan
# Degree of polarization
#given data :
thetai=45 # in degree
n=1.5 #/ index
thetar= degrees(thetar)
Rl=sin((thetai-thetar)*pi/180)**2/sin((thetai+thetar)*pi/180)**2
Rp=tan(thetai-thetar*pi/180)**2/tan((thetai+thetar)*pi/180)**2
D=((Rl-Rp)/(Rl+Rp))*100
print "Degree of polarization, D = %0.2f %%" %D
# answer is wrong in the textbook

Degree of polarization, D = 49.44 %


## Example 6, page 594¶

In :
# Frequency
#given data :
Del=1 # in m
mu=4*pi*10**-7 # in H/m
sigma=4 # in siemen/m
v=1*10**-3/(pi*Del**2*mu*sigma)
print "Frequency, v = %0.1f kHz " %v

Frequency, v = 63.3 kHz