##Chapter 10, Example 10.1, Page 280
import math
## Thermal utilization factor
Summation = ((0.0055*103.4)+(0.720*687)+(99.2745*2.73))/100
sigma = 0.0034
f = 7.662/(7.662+(sigma*450.))
print'%s %.2f %s'%("Total thermal macroscopic = ",Summation," N^U cm^1\n")
print'%s %.2f %s'%(" f = \n",f,"")
## Answer may vary due to round off error
##Chapter 10, Example 10.2, Page 280
import math
## Thermal fission factor
neeta = (2.42*587.)/(687. +(2.73*0.98/0.02))
print'%s %.2f %s'%(" Thermal fission factor = \n",neeta,"")
## Answer may vary due to round off error
##Chapter 10, Example 10.3, Page 282
import math
## Find the probability
P = math.exp(-6.85*10**-4*368.)
Pnl = 1./(1.+(578.*6.85*10**-4))
print'%s %.2f %s'%("Fast-neutron nonleakage probability = \n",P,"")
print'%s %.2f %s'%(" Thermal-neutron nonleakage probability = \n",Pnl,"")
## Answer may vary due to round off error
##Chapter 10, Example 10.4, Page 283
import math
## k of a homogeneous
f = 687./(687. +(0.0034*40000.))
k = 2.07*f
print'%s %.2f %s'%("f =\n",f,"")
print'%s %.2f %s'%(" k = \n",k,"")
##Answer may vary due to round off error
##Chapter 10, Example 10.5, Page 284
import math
##Calculate radius R
L = 578.
T = 368.
Bc = 6.358*10**-4
R = math.sqrt(math.pi**2/Bc)
print'%s %.2f %s'%(" R = ",R," cm \n")
##Answer may vary due to round off error
##Chapter 10, Example 10.6, Page 285
import math
## mass of U235
m = (((4./3.)*math.pi*125**3*1.60)*235.)/(40000.*12.)
print'%s %.2f %s'%(" m = ",m*10**-3," kg \n")
##Answer may vary due to round off error
##Chapter 10, Example 10.7, Page 285
import math
## Keff
Keff = 1./(1.-0.0065*0.1)
print'%s %.2f %s'%(" Keff = \n",Keff,"")
##Answer may vary due to round off error
##Chapter 10, Example 10.8, Page 293
import math
##Resulting reactor period
bt = 0.0065
dt = 0.00065
T = (bt*12.8)/dt
Pt = 10000.
P0 = 10.
t = T*math.log(Pt/P0)
print'%s %.2f %s'%(" Resulting reactor period = ",T," sec \n")
print'%s %.2f %s'%(" t = ",t," sec\n")
##Answer may vary due to round off error