%1. RUN MC SIMULATION BY CALLING:  mc_static(rep,N,T,R0,Rmax,beta0,outputfile)
%   ===> OUTPUT: outputfile.mat 

%2. CALCULATE SUMMARY STATISTICS FOR EACH MC SIMULATION BY CALLING: create_statistics
%   ===> OUTPUT: outputfile-stats.cls 


mc_static(10000,300,10,2,5,1,'results_T10')
create_statistics('results_T10')

mc_static(10000,300,30,2,5,1,'results_T30')
create_statistics('results_T30')

mc_static(10000,300,100,2,5,1,'results_T100')
create_statistics('results_T100')

mc_static(10000,300,300,2,5,1,'results_T300')
create_statistics('results_T300')






