A church window

on 19 Oct 2021

% Draw circle

m=666;

t=linspace(0,2*pi,m) ;

R=.5;

% Generate random numbers

x=rand(m,1);y=rand(m,1);

% Get random points inside the circle

i=inpolygon(x,y,R+R*cos(t),R+R*sin(t));

% Generate delaunay

d=delaunayTriangulation(x(i),y(i));

% z=rand(nnz(i),3);

% plot

figure(color='k')

patch(faces=d.ConnectivityList,vertices=d.Points,facevertexc=rand(nnz(i),3))

shading faceted

axis equal off

material([1 1 R 15 R])

camlight