clear all

p=.4; %fraction of cells that survive each reproduction

%Define vector n(t):  number of cells per generation
n(1)=1000; %initial number of cells (generation 0)
generation_end=20;  %predict up to this generation in the future

for t=2:generation_end
    n(t)=4*n(t-1)*p ;
end

figure
plot(n,'.','MarkerSize',24)
ylabel('n(t)','FontSize',20)
xlabel('t','FontSize',20)
set(gca,'FontSize',20)