function F = modelGradacji(t,N) 
% N – populacja szkodnikow
% t - czas
rB = 0.6; % wspolczynnik rozrodczosci
KB = 0.5; % pojemno?? srodowiska zwiazana z ilo?ci? igiel na drzewach
% funkcja opisujaca drapieznictwo ze strony ptakow
A = sqrt(0.3);
B = 4;% parametry A i B okre?laj? drapie?nictwo.


F = rB*N.*(1-N/KB)-B*N.^2./(A^2+N.^2);