Figure 13 (page 26):

Measurements of species concentrations versus time for parallel reactions.

Code for Figure 13

Text of the GNU GPL.

main.m

clear('all'); close('all');
%
% Create data for the A->B->C model
% jbr,  11/2008

clear all
close all
global theta

k1 = 2.5;
k2 = 0.5;
ca0 = 1;
cb0 = 0;
cc0 = 0;
c0 = [ca0; cb0; cc0];
theta = [k1; k2];
tfinal = 5;
nplot = 100;



measure.states = [1,2,3];
measure.time = linspace(0, tfinal, 20)';

%create the measurements by solving the model and adding noise
[tsolver, cabc] = ode15s(@paraode, measure.time, c0);
randn('seed',0);
y  = cabc(:,measure.states);
measure.data = y + 0.02*randn(size(y));
myfile = fopen('parallel_data.dat', 'w');
for i = 1: size(y,1)
  fprintf(myfile, '%8.2f', measure.time(i), measure.data(i,:));
  fprintf(myfile, '\n');
end
fclose(myfile);

plot(measure.time, measure.data, 'o');

paraode.m

function xdot = paraode(t, x)
  global theta
  ca  = x(1);
  cb  = x(2);
  cc  = x(3);
  k1  = theta(1);
  k2  = theta(2);
  r1 = k1*ca;
  r2 = k2*ca;
  xdot = [-r1-r2; r1; r2];