Ulam's night sky

Eugen Prel
on 30 Oct 2021
- 5
- 45
- 0
- 0
- 279
n = 137;
u = isprime(spiral(n));
figure('Color',[0, 0, 0.08],'Renderer','painters'); hold on;
m = 'ph';
for k = 1:n
for l = 1:n
if u(l,k) ~= 0
scatter(l,k,'SizeData',randi(13),'Marker',m(randi(2)),'MarkerEdgeColor','none','MarkerFaceColor','w','MarkerFaceAlpha',rand);
end
end
end
axis off equal