Untitled

global alpha beta gamma delta0

alpha=0.05; %[1/giorno]
beta=0.00114; %[1/cat/giorno]
gamma=0.000001; %[1/rat/giorno]

% 1-G/70 = limite della popolazione di gatti che il territorio può
% sostenere

% delta = tasso di mortalità dei gatti

delta0=0.001; %[1/giorno]
deltafin=0.015; %[1/giorno] %dopo l'irrorazione con l'insetticida


function dydt=Lotka_Volterra(t,y)

global alpha beta gamma delta0

T=y(1);
G=y(2);

dTdt=alpha*T-beta*T*G;
dGdt=gamma*(1-G/70)*T*G-delta0*G;

dydt=[dTdt; dGdt];

end

T0=20;
G0=20;

y0=[T0; G0];

tspan=[0 7300];

[t, y]=ode45(@Lotka_Volterra, tspan, y0);

figure(1)
plot(t,y(:,1), '-b', 'LineWidth', 1.5);
grid on
xlabel('Tempo(giorni)')
ylabel('Topi')
title ('Numero di topi nel tempo')


figure(2)
plot(t,y(:,2), '-r', 'LineWidth', 1.5);
grid on
xlabel('Tempo(giorni)')
ylabel('Gatti')
title ('Numero di gatti nel tempo')