%FUNKCJA PIERWSZA
clear, clc
r1 = 0.2;
r2 = 0.9;
r3 = 1.1;
r4 = 1.5;
x1 = 0.23;
r = r2;
suma = 0;
for x=2:1:500
   if x1(x-1)<=0.5
      x1 = [x1 r*x1(x-1)];
      suma = suma + log(abs(r));
   else
      x1 = [x1 -r*x1(x-1)+r];
      suma = suma + log(abs(r));
 end
end
N = 500;
lambda = 1/N*suma %"-" - NIE CHAOS, "+" - CHAOS, im wieksa lambda tym WIEKSZY CHAOS!Q!N!@#!!!
x=x1;
plot (x1);



%FUNKCJA DRUGA
clear, clc
r1 = 1.5;
r2 = 3.4;
r3 = 3.5;
r4 = 3.6;
x1 = 0.501;
r = r4;
suma = 0;
for x=2:1:500
   x1 = [x1 r*(x1(x-1))*(1-x1(x-1))];
   suma = suma + log(abs(r-2*r*x1(x-1)));
end
N = 500;
lambda = 1/N*suma %"-" - NIE CHAOS, "+" - CHAOS, im wieksa lambda tym WIEKSZY CHAOS!Q!N!@#!!!
x=x1;
plot (x1);



%FUNKCJA HENNONA
clear, clc
x = 0.5; %od 0 do 1
y = 0.5 %od 0 do 1
a = 0.85; %od -1 do 1
b = -0.56; %od -1 do 1
for i=2:1:500
    y = [y b*x(i-1)];
    x = [x 1-a*(x(i-1))^2+y(i-1)];
end
N = 500;
plot (y,x);


FUNKCJA HENNONA