%This m-file is used to plot the evolution of the mean telomere lenght
%through a user defined number of population doublings

function[pd,M]= length_vs_pd(N_max)

M= zeros(1,N_max);

for N= 1:1:N_max
    [L,p]= load_telomere_data(N,0);
    M(N)= mean((mean(L,2)));
end

pd= 1:N;

figure; hold on;
xlabel('population doubling');
ylabel('mean telomere length (kbp)');
plot(M);

end