##Chapter 9, Example 9.1, Page 241
import math
## Iron kerma and absorbed dose rates
Sp = 10**14
r = 100.
mew = 0.03031
mtr = 0.02112 ## mew/pro
men = 0.01983 ## mew/pro
p0 = 10**-6*Sp*math.exp(-mew*r)/(4.*math.pi*r**2)
K0 = 1.602*10**-10*mtr*p0
D0 = 1.602*10**-10*men*p0
print'%s %.2e %s'%("p0 = ",p0," cm^-2s^-1\n")
print'%s %.2e %s'%(" K0 = ",K0," Gy/s\n")
print'%s %.2e %s'%(" D0 = ",D0," Gy/s\n")
## Answers may vary due to round off error
##Chapter 9, Example 9.2, Page 242
import math
## kerma rate
fsMs = (0.6022/18.)*((2.*12.8*0.5)+(3.5*0.1107))
K = 1.602*10**-10*fsMs*10**10*0.1
print'%s %.2f %s'%("fsUs/p = ",fsMs," cm^2/g\n")
print'%s %.2f %s'%(" K = ",K," Gy/s\n")
## Answers may vary due to round off error
##Chapter 9, Example 9.3, Page 245
import math
##Find fluence and H
Sp = 10**9
dt = 600.
r = 1500.
E = 0.03103
phi = Sp*dt/(4.*math.pi*r**2)
H = 1.602*10**-10*E*phi
print'%s %.2f %s'%("fluence = ",phi," cm^2\n")
print'%s %.2f %s'%(" H = ",H*10**8," microSv\n")
## Answer may vary due to round off error