Path: news.mathworks.com!not-for-mail
From: <HIDDEN>
Newsgroups: comp.soft-sys.matlab
Subject: Incomplete diagram!
Date: Fri, 12 Jul 2013 12:32:20 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 67
Message-ID: <krot0k$q4g$1@newscl01ah.mathworks.com>
Reply-To: <HIDDEN>
NNTP-Posting-Host: www-02-blr.mathworks.com
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: newscl01ah.mathworks.com 1373632340 26768 172.30.248.47 (12 Jul 2013 12:32:20 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Fri, 12 Jul 2013 12:32:20 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 3391687
Xref: news.mathworks.com comp.soft-sys.matlab:799237

Hi, 
I have a problem that I wish to plot what is called a bifurcation diagram. 
The map under study is given by
x_((n+1)r)= 1 - ru *x^2_(nr), 
where r is a positive constant may take 0.1,0.2,...
ru is a positive parameter, and the time t runs from nr to (n+1)r, n=0,1,2,...

I have the matlab code
%%%%%%%%%

% define the vector of values that r takes
r_vec = [0.2];
% define the vector of colors, to plot the dat aof reach value of r in
% different color
color_vec = ['b'];

for rho = 0: 0.005: 2


% 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 = 10;
%  define the nimber of iterations k = 1,.....,Nit
Nit = 1000;
% 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);
time = zeros(1,n);

% initialize x_previous
x_previous = x0;
x_next = x_previous;
time = linspace(0,r,n);
    
  % set the figure details
    figure(1)
    xlabel('t')
    ylabel('x')
    %plot(time,x_next,'*')
    hold on
for i = 1:Nit
    
    x_next = 1 - rho*x_previous.^2;
    x_previous = x_next;
    time = linspace(i*r,(i+1)*r,n);
    
 end
 end


      plot(rho,x_next(1),'.', 'Color', char(color_vec(j)),'linewidth',10)
      
end

hold off
% save the figure
% print -depsc map.eps
%%%%
The code works but I am wondering why I don't get a complete picture of the diagram, it seems that it is divided into 2 parts!!!!
Any help please.