logisticcircle.m

混沌与分岔

clear;clc;
for x1=0.737:0.00001:0.999;
    i=i+1;x2=x1;
    for k=1:100
        x(k)=3.8*x2*(1-x2);
        x2=x(k);
        if k>3
            for j=1:k-1
                if (abs(x1-x(j))<0.000001)
                    x1;
                    j;
                    break;
                end
            end
        end
    end
end