- /
-
Sieve of Eratosthenes
on 4 Oct 2021
- 2
- 38
- 1
- 0
- 259
nOffset = [0 1e6 1e12];
for i=1:3
nX = 90;
nY = 135;
A = nOffset(i) + (1:nX*nY);
A = reshape(A, [nX nY])';
sieveA = isprime(A);
set(gcf,'Color', 'w')
subplot(1,3,i);
imagesc(sieveA); colormap('gray')
axis image; set(gca, 'XTick', [], 'YTick', [])
end