function rysunek = giu(r)
%% funkcja rysuje zaleznosc ilosci nosicieli i pasozytow od predkosci
%% reprodukcji nosicieli
g=zeros(100,1);
f=zeros(100,1);
for R=0:10,
for q=0.1:10,    
for u=0.1:0.1:10,
g(roun(10*u))=R*u*(1-u/q)
f(roun(10*u)=u*u/1+u

end

plot (0.1:0.1:10,g,0.1:0.1:10,f)
end
end
N=zeros(100,50,50);
for k=0.1:0.1:5,
for No = 0.1:0.1:5,
%%N=zeros(100,50,50);
for t=1:100,

N(t,round(k*10),round(10*No))=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