%
% Figure 2.1a
%
load figure_2_1a_data
%
figure;
hax = plotyy( dat(:,1), dat(:,2), ...
        dat(:,1), dat(:,3) );
axesyy( hax, [ 0 3 -200 200 0 300 ], ...
          'Time in days', ...
          'Flowrate in Ml/d', ...
          'COD in mg/l', 2, ...
          [0 100 200], [0 50 100 150] );
%
printsto -f1 -deps2 fig2-1a.eps
%
       
