function rysunek = verhulst(r)
%% funkcja rysuje zaleznosc ilosci nosicieli i pasozytow od predkosci
%% reprodukcji nosicieli
N=zeros(100,50,50);
for k=1:0.1:5,
for No = 1:0.1:5,
%%N=zeros(100,50,50);
for t=1:100,

N(round(k*10),round(10*No),t)=k/(1+(k/No-1)*exp(-r*t));
end
end
end
k=1:50;
No=1:50;
surfl(k,No,N(1));
zlabel('N');
xlabel('k');
ylabel('N0');
title('N(k,No,1)');
%%plot3(-1/a.*log(1./exp(r.*(1-NN./k))),exp(r*(1.-NN./k)),NN);
%%hold off;
%%subplot(2,2,4);
%%surfl(P(index),lambda(index2),NN2);
%%zlabel('N*');
%%xlabel('P');
%%ylabel('lambda');
%%title('N*(P*,lambda) zoom w poblizu lambda ~= 0');
%% koniec tej funkcji