% bifurkációs diagram

kezdet=0.2; %0.5; %mindegy honnan indítok kb ugyanaz jön ki;
iteraciok=500;
mu=0.1:0.001:4;% egész diagram
%mu=3.8:0.00001:3.9; %periodikus ablak 

%{
%Ricker
a=0.1;
mu=0.1:0.01:30; %ez most a "b" paraméter
%}
Population=zeros(length(mu), iteraciok); %adatok ide: sorok: mu-k, oszlopok: iterációk
Population(:,1)=kezdet*ones(length(mu),1); %kezdet minden mu-höz

for i=1:(iteraciok-1) %megcsinálom az iterációt
    Population(:, i+1)=mu'.*Population(:, i).*(1-Population(:, i));
    %Population(:,i+1)=mu'.*Population(:,i).*exp(-a*Population(:,i)); %Ricker
end

figure;
plot(mu, Population(:, 0.8*iteraciok:iteraciok), 'b.', 'MarkerSize', 1);
title('Bikurkációs diagram');
xlabel('mu');
ylabel('populáció');
