- /
-
Personlized Lyapunov Fractal
on 10 Oct 2022
- 57
- 424
- 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])