const = 100;
deltaE = 8.314*77*Log[const]
pVap = 1;
monolayerCap = 1;
amtAdsorbed = 
        monolayerCap*const*
          pressure/((pVap - 
                    pressure)*(1 + (const - 1)*
                      pressure/pVap));
Plot[amtAdsorbed, {pressure, 0, 0.6}];

n = 2;
amtAdsorbed2 = 
        monolayerCap*const*
          relP/(1 - relP)*(1 - (n + 1)*relP^n + 
                n*relP^(n + 1))/(1 + (const - 1)*relP - 
                const*relP^(n + 1));
Plot[amtAdsorbed2, {relP, 0, 1}];

molWt = 28;
density = 0.751;
MolArea = 
      4*0.866*(molWt/(4*2^0.5*6.022*10^23*density))^(2/
              3)

soln1 = 
      Solve[{(c - 1)/vM/c == 
            2*10^-3/(0.205 - 0.065), 1/vM/c == 
            0.2*10^-3}, {vM, c}]

vM_norm = soln1[[1, 1, 2]]/0.606
cNew = soln1[[1, 2, 2]]

deltaE = 1.987*90.1*Log[cNew]

surfaceArea = 116.2/22400*6.022*10^23*16.2*10^-20

Calculation of cm^3/m^2 (see Langmuir 2004 paper citing Brunauer et al., 1938).

(1/0.808)*28/6.022*10^23/16.2*10^-20

molVol = 82.06*273.15 // N

molVol/(6.022*10^23)/(16.2*10^-20)

lrr3@psu.edu (08/29/07)