clc
    clear all
    close all    
%     N=1024;
%     fs=1000;     %采樣頻率
%     f=10;      %信號頻率
%     dt=1/fs;   %采樣周期
%     t=0:dt:(N-1)*dt; %采樣時間間隔
%     ff=0:fs/N:(N-1)*fs/N;
%     data=sin(2*pi*f*t);
    lambda = 3.6;
    x = 0.4*ones(1,length(lambda));
    
    k1 = 4000;                   % 前面的迭代點數
    k2 = 1024;                   % 后面的迭代點數
    
    f = zeros(k1+k2,length(lambda));
    for i = 1:k1+k2
        x = lambda .* x .* (1 - x);
        f(i,:) = x;
    end
    f = f(k1+1:end,:);
    X=f';

    le=largest_lyapunov_exponent(X,1000);