# Variables
import math 
m = 5.313*10**-23 			#g
k = 1.38*10**-16 
T = 298 			#K
R = 82.06 			#ml-atm /mol K
h = 6.624*10**-27 			#J /mol
# Calculations
Qt = (2*math.pi*m*k*T)**1.5*R*T/h**3
# Results
print  'Qt  = %.2e '%(Qt)
# Variables
import math 
Qe = 4.029
Qe1 = -37.02
Qe2 = 4.695*10**4
T = 300 			#K
R = 1.98 			#cal /mol K
Qe3 = 4.158
Qe4 = -200.8
Qe5 = 2.546*10**5
T1 = 500 			#K
# Calculations
Ce = R*((Qe2/Qe)-(Qe1/Qe)**2)/T**2
Ce1 = R*((Qe5/Qe3)-(Qe4/Qe3)**2)/T1**2
# Results
print  'electronic contribution  = %.3f cal deg**-1.g.atom**-1'%(Ce)
print  ' electronic contribution  = %.3f cal deg**-1.g.atom**-1'%(Ce1)
# Variables
import math 
I = 0.459*10**-40 			#g cm**2
k = 1.38*10**-16
T = 300 			#K
h = 6.624*10**-27 			#J/mol
I1 = 245*10**-40 			# g cm**2
# Calculations
Qr = I*k*T*8*math.pi**2*0.5/h**2
Qr1 = I1*k*T*8*math.pi**2/h**2
# Results
print  'Rotational Partition  = %.2f '%(Qr)
print  ' Rotational Partition  = %.f '%(Qr1)
# note : rounding off error. please check.
# Variables
import math 
h = 1.439
T = 300. 			#K
w = 4405. 			#cm**-1
w1 = 565. 			#cm**-1
# Calculations
Qv1 = (1-math.e**(-h*w/T))**-1
Qv2 = (1-math.e**(-h*w1/T))**-1
# Results
print  'Vibrational Partition  = %.3f '%(Qv1)
print  ' Vibrational Partition  = %.3f '%(Qv2)
# Variables
h = 1.439
T = 300 			#K
w = 565 			#cm**-1
R = 1.98 			#cal /mol K
n = 0.56
# Calculations
Qr = h*w/T
Cv = n*R
# Results
print  'Vibrational Partition  = %.2f '%(Qr)
print  ' Cv  = %.2f cal deg**-1 mole**-1'%(Cv)
# Variables
R = 1.986 			#cal deg**-1 mole**-1
cv1 = 0.392 			#cal deg**-1 mole**-1
cv2 = 0.004 			#cal deg**-1 mole**-1
cv3 = 0.003 			#cal deg**-1 mole**-1
cv4 = 1.265 			#cal deg**-1 mole**-1
cv5 = 0.247 			#cal deg**-1 mole**-1
cv6 = 0.225 			#cal deg**-1 mole**-1
# Calculations
Cv  = 3*R+cv1+cv2+cv3
Cv1 = 3*R+cv4+cv5+cv6
# Results
print  'Total heat capacity  = %.2f cal deg**-1 mole**-1'%(Cv)
print  ' Total heat capacity  = %.2f cal deg**-1 mole**-1'%(Cv1)
# Variables
R = 1.98 			#cal/mol K
# Calculations
Cv = 2.856*R
# Results
print  'Cv  = %.2f cal deg**-1 g.atom**-1'%(Cv)
# Variables
R = 1.98 			#cal/mol K
n = 3
# Calculations
Cv = n*R*0.8673
# Results
print  'Cv  = %.2f cal deg**-1 g.atom**-1'%(Cv)
# Variables
R = 1.98 			#cal/mol K
n = 3
# Calculations
Cv = n*R*0.904
# Results
print  'Cv  = %.2f cal deg**-1.g.atom**-1'%(Cv)