Rossler Attractor Bifurcation Diagram

Hello,

I am trying to create a Bifurcation Diagram for the Rossler Attractor system, but for some reason my code only keeps giving me a point. I am not sure where I am going wrong. Any help is much appreciated.

Here is my code so far.

function [ T,X ] = test( a, b0, b1, c, n, Tmid, t, initV, T, eps )

for b = b0:((b1 - b0)/n):b1; options = odeset('RelTol',eps,'AbsTol',[eps eps eps/10]); [T,X] = ode45(@(tt,xx) F(tt, xx, a, b, c), T, initV, options); end

Xmid = X(Tmid:length(T),1); Tmid = T(Tmid:end);

X1 = Xmid(:,1); l = length(Tmid); A = [];

    for j=1:l-t;
        [A] = [A: b max(X1(j:j+t,1))];
    end

figure(1) plot(A(:,1),A(:,2),'ro','MarkerEdgeColor','b','MarkerFaceColor','b','MarkerSize',1.5)

xlabel('b','fontsize',18); ylabel('max(X)','fontsize',18); end

function dx = F( T,X, a, b, c)

dx = zeros(3,1); dx(1) = -X(2) - X(3); dx(2) = X(1) + a*X(2); dx(3) = b + X(3)*(X(1) - c); return

end

Answer 1

[Sanaa]

See this plz http://homepages.ulb.ac.be/~dgonze/INFO/matlab/rossler.m

[Umesh Prajapati]

I tried your code. Result is wrong. We don't get what we need.