

%%% Baseline Run
        
alpha   =  5; 
beta    =  3;
gamma   =  0.3;
dead    =  0.3;
mu1     =  13;
mu2     =  20;
maimed  =  0.05;
phi     = 0.005;

plist = [alpha beta gamma dead mu1 mu2 maimed phi]; 

RunBam18('Baseline-d',plist)

%%% Mu1-Min Run

alpha   =  1;        % Change only  Mu1 Min
beta    =  3;
gamma   =  0.3;
dead    =  0.3;
mu1     =  7;
mu2     =  20;
maimed  =  0.05;
phi     = 0.005;

plist = [alpha beta gamma dead mu1 mu2 maimed phi]; 

RunBam18('Mu1-Min',plist)

%%%  Mu1-Max Run

alpha   =  10;      % Change only to Mu1 Max 
beta    =  3;
gamma   =  0.3;
dead    =  0.3;
mu1     =  17;
mu2     =  20;
maimed  =  0.05;
phi     = 0.005;

plist = [alpha beta gamma dead mu1 mu2 maimed phi]; 

RunBam18('Mu1-Max',plist)


%%% Mu1-Min Run

alpha   =  1;        % Change only  Mu1 Min
beta    =  3;
gamma   =  0.3;
dead    =  0.3;
mu1     =  13;
mu2     =  16;
maimed  =  0.05;
phi     = 0.005;

plist = [alpha beta gamma dead mu1 mu2 maimed phi]; 

RunBam18('Mu2-Min',plist)

%%%  Mu1-Max Run

alpha   =  10;      % Change only to Mu1 Max 
beta    =  3;
gamma   =  0.3;
dead    =  0.3;
mu1     =  13;
mu2     =  24;
maimed  =  0.05;
phi     = 0.005;

plist = [alpha beta gamma dead mu1 mu2 maimed phi]; 

RunBam18('Mu2-Max',plist)