clear all
close all
clc

% dati di riferimento
y0 = 50;
m_vec = 0.01:0.01:0.9;
b = 0;
x = zeros(1,500); % tempi fino a 500!!
n = 0:length(x);

% inizializzo vettore anno estinzione
anno_estinzione = zeros(1,length(m_vec));

% calcolo l'anno di estinzione al variare del parametro m
for k = 1:length(m_vec)
    % faccio girare il sistema
    y = [y0, popolazione(y0, m_vec(k), b, x)];
    % trovo l'anno di estinzione
    anno = n(find(y<1,1));
    if isempty(anno)
        anno = max(n); % di default, se non si estingue 
                       % metto l'ultimo anno considerato
    end
    anno_estinzione(k) = anno;
end

% rappresento il risultato
figure
plot(m_vec, anno_estinzione,'b-')
grid on
xlabel('tasso di mortalità m')
ylabel('anno di estinzione')