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


reactions = {'B + P --> I1'      [ 0.1] ;    
             'I1 + P --> I2'     [0.1]  ; 
             'I2 + P --> I3'     [0.1]  ; 
             'I1 --> L'          [0.06667]  ;              
             'I1 --> K'          [ 0.6]  ;              
             'I2 --> L'          [0.06667]  ;              
             'I2 --> K'          [0.15556]  ;              
             'I3 --> L'          [0.06667 ]  ;              
             'I3 --> K'          [0.06667]  ;              
             'K --> 100P'        [0.025]; 
             'B + FA --> 2B'      [0.05];
             'L + FA --> 2L'      [0.05];
             'FT --> FA'         [0.05]};

rev_reactions = {};

species = {'B';
             'P';
             'I1';
             'I2';
             'I3';
             'L';
             'K';
             'FA';
             'FT'};


%InitialValues
%B =
init(:,:,1) = [1000 1000 ;
               1000 1000 ];
           
%P =
init(:,:,2) = [100 0;
               0 0 ];

%I1 =
init(:,:,3) = [0 0 ;
               0 0 ];           
%I2 =
init(:,:,4) = [0 0 ;
               0 0 ];
%I3 =
init(:,:,5) = [0 0 ;
               0 0 ];
%L =
init(:,:,6) = [0 0 ;
               0 0 ];
%K =
init(:,:,7) = [0 0 ;
               0 0 ];
%FA =
init(:,:,8) = [0 0 ;
               0 0 ];
           
%FT =
init(:,:,9) = [1e+6 1e+6 ;
               1e+6 1e+6 ];
               

CV = {'FA' 'FT'};
           
           
start = 0;
stop = 120;
increment = 1;
                 
         
