function[CV reactions rev_reactions init species start stop increment] = chemMonad()
%----Modify the following code to change the model-----


reactions = {'B --> 2B' [.013]};

rev_reactions = {};

species = {'B';
           };


%InitialValues
%B =
init(:,:,1) = [8.1612e+003;
               ];
           


CV = {};
           
           
start = 0;
stop = 200;
increment = 1;
                 
         
