Thanks a lot for explanations. I have tried my own code to plot the bifurcation diagram but still have one problem
% define the vector of colors, to plot the dat aof reach value of r in % different color color_vec = ['b']; r=0.05; for rho = 0: 0.001:10
% do a loop of values of r for j = 1 : size(r_vec,2) % define delay r = r_vec(j); % define coefficient rho %rho = 0.5; % define the number of deiscrete times in interval [k*r, (k+1)*r] n = 50; % define the nimber of iterations k = 1,.....,Nit Nit = 300; % define the initial vector which is a vector of size n x0 = 0.3*ones(1,n); % define vectors x_next and x_previous x_next = zeros(1,n); x_previous = zeros(1,n); %x_n x_pp=zeros(1,n); %x_n-1 x_ppp=zeros(1,n); %x_n-2 x_pppp=zeros(1,n); %x_n-3 time = zeros(1,n);
hold on for i = 1:Nit time = linspace(i*r,(i+1)*r,Nit); x_next = x_previous+(time(i)-i.*r_vec(j))*rho*x_pppp.*(1-x_ppp); x_previous = x_next; end end plot(rho,x_next(1),'.', 'Color', char(color_vec(j)))
I always get a straight line. What went wrong please? Thanks a lot for your help. I really appreciate it.