#calculate the Volume
#Initialization of variables
V1=0.284 #cm^3 /g
V2=1.43 #cm^3 /g
P1=142.4 #mm
P2=760. #mm
#calculations
z=(1/V1 - 1/V2)/(1/P1 - 1/P2)
invVm=1/V2 - z/P2
Vm=1/invVm
#results
print '%s %.1f %s' %("Volume =",Vm,"cm^3/g")
print 'The answer in the textbook is a bit different due to rounding off error.'
#calculate the total area
#Initialization of variables
Vm=2.86 #cc/g
P=1 #atm
R=82.06 #cm^3 atm/deg mol
T=273.2 #deg
N=6.023*10**23
sigma=16.2*10**-16 #cm^2 /molecule
#calculations
n=P*Vm/(R*T)
A=N*n*sigma
#results
print '%s %.2e %s' %("total area =",A," cm^2 (g catalyst)^-1")