function [A, b] = chemreaction()
%CHEMREACTION Initialize the linear system describibg the particular chemical reaction
A = zeros(4,4);
b = zeros(4,1);
A(1,1) =  1.;
A(1,3) = -1.;
A(2,1) =  2.;
A(2,2) =  3.;
A(2,3) = -6.;
A(2,4) = -1.;
A(3,1) =  2.;
A(3,2) =  1.;
A(3,4) = -2.;
A(4,2) =  1.;
A(4,3) = -2.;
end

