- /
-
Infinity
on 4 Oct 2021
- 7
- 32
- 1
- 0
- 242
C = [flipud(gray); hot];
C = C(16:end-64,:);
N = 1000;
t = 2*pi*(1:N)/N;
x = cos(t)./(1+sin(t).^2);
z = sin(t);
y = x.*z;
r = 1e3*(.5 + abs(x)).^2;
C = C(1+round((length(C)-1)*abs(x)),:);
scatter3(x,y,z,r,C,'filled')
view(0,90)
axis equal off
set(gcf,'color','k')