% 
% Written by:
% -- 
% John L. Weatherwax, Ph.D.         2002-11-16
% MIT Lincoln Laboratory
% 244 Wood Street
% Lexington, MA 02420-9176 USA
% 
% email: weatherwax@ll.mit.edu
% Voice: (781) 981-5370       Fax: (781) 981-0721
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

close all;

tK = linspace( 50, 100, 100 ) + 273;
L  = 41.790; R = 8.314;
P0 = 0.1234; 

C0 = P0*exp( L/( R*373 ) );

pBar = C0*exp( -L./( R*tK ) );

figure;

plot( tK, pBar, '-o' );
xlabel( 'temp (Kelvin)' ); ylabel( 'vapor pressure (Bar)' ); title( 'vapor pressure v.s. t' );

grid on;

heightsM = [ 1430 3090 4420 8850 ];
mGOverRT = 9.822*10^(-5);
pHeightAtm = exp( -heightsM*mGOverRT )