Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: why I am not getting a bifurcation diagram?
Replies: 3   Last Post: Feb 20, 2013 1:20 PM

 Messages: [ Previous | Next ]
 Sanaa Posts: 171 Registered: 3/20/12
why I am not getting a bifurcation diagram?
Posted: Feb 20, 2013 7:20 AM

Hi,
I am trying to plot a bifurcation diagram for Rossler discrete-time system.
What's wrong in my code?
I am only getting a line!!!!

itermax=300;min=itermax-9;% That is we plot from 291:300 i.e. 10 values of x only.
r=0.1;b=0.2;c=4;
for a=0:0.001:0.4
x=0.1;
x0=x;
y=0.2;
y0=y;
z=0.3;
z0=z;
for n=1:itermax
xn=x0-r*(y0+z0);
yn=y0+r*(x0+a*y0);
zn=z0+r*(b+x0*z0-c*z0);
x=[x xn];
x0=xn;
y=[y yn];
y0=yn;
z=[z zn];
z0=zn;
end
plot(a*ones(10),y(min:itermax),'b.','linewidth',2)
hold on
end
xminmax = xlim;
line(xminmax, [0 0], 'Linewidth', 1,'Color','k')
fsize=15;
xlabel('a','FontSize',fsize)
ylabel('\ity','FontSize',fsize)
hold off
%print(gcf, '-djpeg', '-zbuffer', 'bif01_standard.png');

Any help will be highly appreciated.
Thanks

Date Subject Author
2/20/13 Sanaa
2/20/13 dpb
2/20/13 Sanaa
2/20/13 dpb