function y=NO2(t,u,k11,k12,k21,k22)
% Oxidation av NO till NO2.
% Atkins och Jones sid 581-583.  

% reaktionshastigheter
r11=k11*u(1)^2;  
r12=k12*u(3);
r21=k21*u(2)*u(3);
r22=k22*u(4)^2;

y=zeros(size(u));

% hastighetsekvation
y(1)=-2*r11+2*r12;     % NO
y(2)=-r21+r22;         % O2
y(3)=r11-r12-r21+r22;  % N2O2
y(4)=2*r21-2*r22;      % NO2

