clear; popsize=2000; % population size tmax=50000; % number of seasons penalty=0.4; % fecundity cost of resampling, 0.4 i standard cost reprcost=0.01; % survival cost of reproductive effort mu0=0.05; % baseline mortality sampl_init=logspace(-2,1,7); %starting values for sampling effort repr_init=linspace(1,7,3); % starting values for reproductive effort mutations=[1 1]; %where do mutations occur pvec=linspace(0.4,0.9,6); % environmental autocorrelation values determined by p evec=linspace(0.01,0.21,3); % sampling error values % Initialise matrices to hold data sampleffort=nan(length(sampl_init),tmax,length(repr_init),length(pvec),length(evec)); repreffort=nan(length(sampl_init),tmax,length(repr_init),length(pvec),length(evec)); clutchsize=nan(length(sampl_init),tmax,length(repr_init),length(pvec),length(evec)); maladaptation=nan(length(sampl_init),tmax,length(repr_init),length(pvec),length(evec)); PopDist=nan(length(sampl_init),popsize,6,length(repr_init),length(pvec),length(evec)); MisMatch=nan(length(sampl_init),length(repr_init),length(pvec),length(evec)); A=nan(length(sampl_init),length(repr_init),length(pvec),length(evec)); countdown=length(sampl_init)*length(repr_init)*length(pvec)*length(evec); % Run the simulation for all different parameter values for p=1:length(pvec) for epsilon=1:length(evec) for i=1:length(sampl_init) for j=1:length(repr_init) [sampleffort(i,:,j,p,epsilon),repreffort(i,:,j,p,epsilon),clutchsize(i,:,j,p,epsilon),maladaptation(i,:,j,p,epsilon),MisMatch(i,j,p,epsilon),A(i,j,p,epsilon)]=evomemory(popsize,tmax,pvec(p),evec(epsilon),penalty,reprcost,mu0,[1 sampl_init(i) repr_init(j)],mutations); countdown=countdown-1 %this is just to print how many simulation runs are left end end end end %end save lifeplast