x0=0.01;
a=3.8;
N=100;
x=zeros(N,1);
x(1)=x0;
for i=2:N
    x(i)=f_logistic(x(i-1),a);
end

% take the first three digits
y0=round(x(10)*10^3)/10^3;
y=zeros(N-10,1);
y(1)=y0;
for i=2:N-10
    y(i)=f_logistic(y(i-1),a);
end

disp('Lorentz experiment:');
disp(['x(10)=', num2str(x(10)),', y(1)=',num2str(y(1))]);
disp(['x(20)=', num2str(x(20)),', y(10)=',num2str(y(10))]);
disp(['x(50)=', num2str(x(50)),', y(40)=',num2str(y(40))]);
disp(['x(100)=', num2str(x(100)),', y(90)=',num2str(y(90))]);