function rysunek = pasozyty(a,c)%%,k,r)
%% funkcja rysuje zaleznosc ilosci nosicieli i pasozytow od predkosci
%% reprodukcji nosicieli
lambda=0.1:0.1:1000;
P=zeros(10000,1);
N=zeros(10000,1);
id=1:10000;

P(round(abs(lambda*10)))=-1/a*log(1./lambda);
for indx=1:10000,
    if P(indx)<=0,
        P(indx)=0;
    end 
end 
N(round(abs(lambda*10)))=lambda./(c*a*(1.-lambda)).*log(1./lambda);

if N(id)<=0,
    N(id)=0;
end 
NN=zeros(100,100);
%%NN2=zeros(100,100);
index=1:100;
index2=1:100;
for index3=1:100,
    for index4=1:100,
NN(index3,index4)=lambda(10*index4)/(c*(lambda(10*index4)-1))*P(index3);
if P(index3)<=0
    NN(index3,index4)=lambda(10*index4)/(c*(lambda(10*index4)-1))*P(index3);
end
    %%NN2(index3,index4)=lambda(index4)/(c*(lambda(index4)-1))*P(index3);
if NN(index3,index4)<=0,
    NN(index3,index4)=0;
end
%%if NN2(index3,index4)<=0,
%%    NN2(index3,index4)=0;
%%end
    end
end
subplot(2,2,1);

plot(lambda,P,'-r');
title('P*(lambda)');

subplot(2,2,3);

plot(lambda,N,'-b');
title('N*(lambda)');
subplot(2,2,[2;4]);
%%hold on;
surfl(P(index),lambda(10*index2),NN);
zlabel('N*');
xlabel('P');
ylabel('lambda');
title('N*(P*,lambda)');
%%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
%% function [element,en]=el(c,n)
%% element=0.3e-4*c*n;
%% element=1000/6*c*n;
%% en=n+12/1000-element;
%% element=2/5*c*n;
%% en=n+5-element;
%% koniec tej funkcji