Tspan = linspace(0,20,1000);

k=1;
%for c0 = 0:0.2:2
   %c0
   [t u] = ode15s('diff', Tspan, [0 5 0 0 0]);
%   max_m(k) = max(u(:,3));
%   figure(3)
%   hold on 
%   plot(t,u(:,3))
%   k=k+1;
%end

%b)
figure(1)
plot(t,u)
xlabel('Tid (s)')
ylabel('Koncentration (M)')
grid
legend('C_6H_6','O_2','C_4H_2O_3','CO_2','H_2O')

%b) ii)
%figure(2)
%plot(0:0.2:2,max_m)
%grid
