subm_calibrate_colicin.m

function z=calibrate_colicin

opt = optimset('MaxIter', 5000, 'MaxFunEvals', 5000);

%addpath('../Data');
subm_fetch_data;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%SSE=1225.146 %eps=0.001, cov=0.01, 
params=[0.114471500852052
   1.177541496882635
   0.211658150775181
   9.666505932271139
   0.855708199615597
   0.317377470424927
   0.421198591219959
   0.223911224561197
   0.001100072834240
   0.006155849122120
   0.001713291744643
   0.007084442015885
   0.002343229576080
   0.009261500126116
   0.002221571067967
   0.006992461503957
   0.003102900711591
   0.003007528602092
   0.003979439128267
   0.004125429804641
   0.005001246130118
   0.002519631469717
   0.004162155718298
   0.002635664694760
   0.004996158765916
   0.001611042006734
   0.002589779939195
   0.000000201235590
   0.001032964357047
   0.000000249393925
   0.002290538203115
   0.000000125346725];





sizep=length(params);



lb=0*ones(1,sizep); ub=100*ones(1,sizep); % Set lower and upper bounds (lb and ub)

% Comment out one of fmincon or simulannealbnd
z=fmincon(@subm_run_colicin,params,[],[],[],[],lb, ub,[], opt);
%z=simulannealbnd(@subm_run_colicin,params,lb,ub, opt);
format long
z %report output


end