Moon 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

scatter([0,.5],[0,0],5e5,[1;0]*[1,1,1],'.');

% Add the panes

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

FaceA=.4);

% Pretty...

shading faceted

axis equal off

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