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')
Leave a Comment