%
% figure 2-2b
%
load influent_cod.dat
t = [0:length(influent_cod)-1]' * 0.1 / 24;
%
plot(t,influent_cod);
set( gca, 'fontsize', 20 );
xlabel( 'Time in days' );
ylabel( 'Influent COD in mg/l' );
%
print -f1 -deps2 fig2-2b.eps
%
