A smooth and hypnotic dot-tunnel animation created using Nelson

[Nelson-numerical-software]

% dot_tunnel_image2.m  (NELSON)

%
% PARAMETERS
%
nRings        = 60;     % number of rings
nDotsPerRing  = 55;     % dots on each ring
maxR          = 120;    % outer radius
minR          = 1.5;    % center threshold
speedShrink   = 0.94;   % shrink factor per frame
rotSpeed      = 0.015;  % rotation per frame (set 0 for static rings)

%
% RADIAL DISTRIBUTION
% (non-linear spacing → deeper/tunnel feel)
%
radial = linspace(0, 1, nRings);
radial = radial .^ 1.4;     % curve exponent (tweak 1.2–1.8)
RR0 = minR + (maxR - minR) * radial;   % starting radii

%
% DOT ANGLES (evenly spaced)
%
theta = linspace(0, 2*pi, nDotsPerRing+1);
theta(end) = [];  % remove duplicate

[TH, RR] = meshgrid(theta, RR0);

%
% INITIAL COORDINATES
%
X = RR .* cos(TH);
Y = RR .* sin(TH);

%
% FIGURE SETUP
%
figure(Color='k');
axis('equal')
axis('off')
axis([-maxR maxR -maxR maxR]);

% grayscale intensity by radius → bright in center
brightness = 1 - (RR(:)./maxR);
brightness = max(0.2, brightness);
col = [brightness brightness brightness];

h = scatter(X(:), Y(:), 18, col, 'filled');

%
% MAIN LOOP
%
while ishandle(h)

    % Slight rotation for “flow” effect (optional)
    TH = TH + rotSpeed;

    % Shrink radii toward center
    RR = RR * speedShrink;

    % Respawn rings that reached center
    idx = RR < minR;
    if any(idx(:))
        RR(idx) = maxR;
    end

    % Update coordinates
    X = RR .* cos(TH);
    Y = RR .* sin(TH);

    % Update brightness (Image #2 glow effect)
    brightness = 1 - (RR(:) ./ maxR);
    brightness = max(0.15, brightness);
    col = [brightness brightness brightness];

    %
    % UPDATE SCATTER
    %
    set(h, ...
        'XData', X(:), ...
        'YData', Y(:), ...
        'CData', col);

    drawnow;
end

https://youtu.be/tzA4O23hRgs