- /
-
Sacks prime spiral
on 31 Oct 2021
- 6
- 63
- 0
- 0
- 200
figure('Color','k');
m=200;
r=1:m*m-1;
t=[];
for n=2:m
l=linspace(0,2*pi,2*n);
t=[t,l(1:end-1)];
end
y=isprime(r);
t=t(y);
r=r(y);
polarplot(t,r,'Marker','.','LineStyle','none','MarkerSize',2.2,'Color','#FFB976');
axis off