% Oxidation av NO till NO2.
% Atkins och Jones sid 581-583.  

T=1;   % stopptid

% två-stegs reaktionskonstanterna 
k11=10;    % snabb
k12=10;    % snabb
k21=0.01;  % långsam
k22=0;     % långsam

% begynnelsevärden 
u0=[0.5;1.0;0;0];

[t,U]=ode45(@(t,u)NO2(t,u,k11,k12,k21,k22),[0,T],u0);   % MATLABs standard ODE-lösare

figure(1)
clf
plot(t,U)

k=2*k11*k21/k12;  % reaktionskonstanten
x=NO2a(U,k11,k12,k21,k22)/k;  % skall vara = 1 

hold on 
plot(t,x,'k:')
hold off
legend('[NO]','[O_2]','[N_2O_2]','[NO_2]','u_4''/(u_2 u_1^2)')
axis([0 t(end) 0 2])
