- /
-
Vortex
on 12 Oct 2021
- 10
- 21
- 0
- 0
- 225
c = 1.68;
d = 2*pi/c;
theta = 0:d:1000;
r = sqrt(theta);
x = r.*cos(theta);
y = r.*sin(theta);
[v,c] = voronoin([x',y']);
for i = 1:length(c)
if all(c{i}~=1)
xp = v(c{i},1);
yp = v(c{i},2);
patch(xp,yp,i);
end
end
axis square off
axis(20*[-1 1 -1 1])