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


reactions = {'X --> 2B' [.05];
             'B + P --> I1'      [ 0.00000001] ;
             'X + P --> I1'      [0.00000001];
             'I1 --> K'          [ 0.9]  ;
             'K --> 100P + S'        [0.025];
             'I1 --> L'          [ .0005];
             'Y --> 2L'         [.05]};
             

rev_reactions = {'B + S <--> X' [5.12829e-7] [5e-7];
                 'L + S <--> Y' [5.12829e-7] [5e-7]};
                  
                 

species = {'B';
           'S';
           'X';
           'P';
           'K';
           'I1';
           'L';
           'Y';
            };


%InitialValues
%B =
init(:,:,1) = [1.7952e+004;
               ];
           
%S =
init(:,:,2) = [6.5e4;
               ];

%X =
init(:,:,3) = [0;
               ];     
 
%P = 
init(:,:,4) = [100000];

%K = 
init(:,:,5) = [0];

%I1 = 
init(:,:,6) = [0];

%L = 
init(:,:,7) = [0];

%Y = 
init(:,:,8) = [0];

% %ZS = 
% init(:,:,9) = [3e4];
% 



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