Sun through the flower window

on 20 Oct 2022

p = rand(70); % Radii

g = rand(70)*7; % Angles

p(p<(cos(g*4)+.1).^2/1.5+.1) = 0; % Flower function

% Triangulate

x=p.*cos(g);

y=p.*sin(g);

p = delaunay(x,y);

Warning: Duplicate data points have been detected and removed.
Some point indices will not be referenced by the triangulation.

% Add the sun

plot(0,0,'y.','markersize',650)

% Add the panes

h = patch('faces',p,'vertices',[x(:) y(:)], 'FaceVertexCData', hsv(size(p, 1)),...

'FaceAlpha', .55);

% Pretty...

shading faceted

axis equal off

set(gcf, 'color', 'k');