Code for Figure 1

main.m

clear('all'); close('all');
% 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.

npts  = 251;
x = linspace(-5, 7, npts)';
sigma = [1/2, 1, 2];
y = [ exp(-0.5*(x-1).^2./sigma(1)^2)/sqrt(2*pi)/sigma(1), ...
exp(-0.5*(x-1).^2./sigma(2)^2)/sqrt(2*pi)/sigma(2), ...
exp(-0.5*(x-1).^2./sigma(3)^2)/sqrt(2*pi)/sigma(3) ];

table = [x, y];

save -ascii normal1d.dat table;

plot (x, y)
title ('Figure 1')