#calculate the Increase in kinetic energy
#Initialization of variables
m1=1.008142
m2=1.008982
#calculations
dm=m1-m2
dt=abs(dm) *931
#results
print '%s %.3f %s' %("Increase in kinetic energy =",dt,"Mev")
#calculate the Thickness
#Initialization of variables
import math
d=8.642 #g/cc
M=112.41 #g/mol
ratio=0.01/100
nb=2400
#calculations
n=d*6.02*10**23 /M
sigma=nb*10**-24
x=-2.303*math.log10(ratio) /(sigma*n)
#results
print '%s %.3f %s' %("Thickness =",x,"cm")
#calculate the Threshold
#Initialization of variables
M1=4
M2=14
E=-1.2 #Mev
#calculations
R1=1.5*10**-13 *(M1)**(1/3.)
R2=1.5*10**-13 *(M2)**(1/3.)
V1=2*7*(4.8*10**-10)**2 /(R1+R2)
V2=V1/(1.6*10**-6)
x=(M1+M2)*V2/M2
#results
print '%s %.1f %s' %("Threshold =",x," Mev")
#calculate the Time taken
#Initialization of variables
import math
t=1622. #years
per=1. #percent
#calculations
Nratio=1-per/100.
x=t*math.log10(Nratio) / math.log10(0.5)
#results
print '%s %.1f %s' %("Time taken =",x," years")