n0 = [50; 100];   
D = [ 1  -1   0    
      0   1  -1 ];
c = [1 .005 0.6];

L = [1 1 0
     0 1 1];               % Stoich coeffs of the reactants
K = sum(L);                    % molecularity
tf = 30;                      % simulation stop time in seconds 
runs = 1;                     % number of realizations(simulation runs)