• Remix
  • Share
  • New Entry

  • MvLevi

  • /
  • Personlized Lyapunov Fractal

on 10 Oct 2022
  • 57
  • 350
  • 0
  • 2
  • 269
C=[-9:9e-3:9];
D=[-9:9e-3:9];
for q=1:2001
for j=1:2001
X=.5;
for i=1:5
if mod(i,2)==0
X(i+1)=X(i)-C(q)*(.5+.3*cos(X(i)))^-1;
else
X(i+1)=X(i)-D(j)*(.5+.3*cos(X(i)))^-1;
end
end
P=diff(X);
L(q,j)=mean(log(abs(P)));
end
end
pcolor(C,D,-L)
shading flat
colormap(twilight)
axis off
caxis([-3.5 3.5])
Remix Tree