Untitled

%ESERCIZIO 2
%LOTKA-VOLTERRA
alfa=0.05; %[1/giorno]
beta=0.00114; %[1/gatto/giorno]
gamma=0.000001; %[1/ratto/giorno]
%(1-G/70) limite della popolazione di gatti che il territorio può sostenere 
%delta è il tasso di mortalità dei gatti
delta0=0.001; %[1/giorno]
deltafin=0.015; %[1/giorno] (dopo l'irrorazione con insetticida)
%calcola numero di gatti G e di topi T
tspan=[0 7300];
function dydt=lotka_volterra(t,y)
alfa=0.05; 
beta=0.00114; 
gamma=0.000001;
delta0=0.001;
T=y(1);
G=y(2);
T0=10;
G0=10;
dTdt=alfa*T-beta*T*G;
dGdt=gamma*(1-G/70)*T*G-delta0*G;
dydt=[dTdt dGdt];
end