Figure 5.25 (page 114):

Adsorbed oxygen concentration versus gas-phase oxygen concentration.

Code for Figure 5.25

Text of the GNU GPL.

main.m

%% Copyright (C) 2001, James B. Rawlings and John G. Ekerdt
%%
%% This program is free software; you can redistribute it and/or
%% modify it under the terms of the GNU General Public License as
%% published by the Free Software Foundation; either version 2, or (at
%% your option) any later version.
%%
%% This program is distributed in the hope that it will be useful, but
%% WITHOUT ANY WARRANTY; without even the implied warranty of
%% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
%% General Public License for more details.
%%
%% You should have received a copy of the GNU General Public License
%% along with this program; see the file COPYING.  If not, write to
%% the Free Software Foundation, 59 Temple Place - Suite 330, Boston,
%% MA 02111-1307, USA.

%
% This m-file is for Spring 1995 homework problem number 18.
% This m-file is for prob18 in ~/che372/problems/kinetics
%
% Make some data with a little bit of noise for oxygen adsorption
%
% gas constant and experimental conditions
R = 0.082;
T = 400;
Ndata = 12;
P = [linspace(.1,.5,5) linspace(2,5,Ndata-5)];

% Langmuir model parameters
cm = 9.6e-6;
K = 12;
c_gas = [P/(R*T)]';


% Surface concentration
tmp  =  sqrt(K*c_gas);
c_surf = cm*tmp  ./ (1 + tmp);

%
% Put some noise in the data
%
randn('seed',0)
c_surf_noise = c_surf + randn(Ndata,1)*0.25e-6;
table = [c_gas c_surf_noise c_surf 1./sqrt(c_gas) 1./c_surf_noise 1./c_surf];

plot(table(:,1),table(:,2),'o');
axis ([-0.01,0.16,0,6e-6]);
title ('Figure 5.25')