//Model Lotki-Volterry
function Ndot = LV(t,y)
    Ndot = [b*y(1)-p1*y(1)*y(2);
            -d*y(2)+p2*y(1)*y(2)];
endfunction
b = 1; //Dzietność
d = 3; //Umieralność
p1 = 0.8;
p2 = 0.5;
t = 0:0.1:20; //Oś czasu
y0 = [100;200];

k = 0.1; //Wsp. konkurencji

y = ode(y0,t(1),t,LV);
plot(t,y(1,:),'blue',t,y(2,:),'red');

xtitle('Model Lotki-Volterry','Czas [min]','Liczebność');

//Model Verhulsta