- /
-
Inverted Sunspots
on 26 Oct 2022
- 19
- 71
- 0
- 2
- 280
N=1e8;
x=single(zeros(N,1)); % required for speed up
y=x;
F=@(x,y)[,];
a=-16;b=16;c=1;d=-1;
% Clifford Attractor - discrete system
for k=1:N-1
x(k+1)=sin(a*y(k))+cos(a*x(k));
y(k+1)=sin(b*x(k))-cos(b*y(k));
end
I=zeros(900);
E=1;
D=size(I,1)-2*E;
for k=1:N
W=ceil(D*(x(k)+2)/4)+E;
U=ceil(D*(y(k)+2)/4)+E;
I(W,U)=I(W,U)+.002;
end
imshow(I,'Colormap',hot);
set(gcf,'Color','k');
%axis tight;