- /
-
Galaxy
on 12 Oct 2021
- 3
- 27
- 0
- 0
- 278
figure('color','k')
R=(randn(1,1e6))+2;
i=find(abs(R)<4);
R=R(i);
T=3*sqrt(pi)*randn(1,1e6);
X=R.*cos(T(i));
Y=R.*sin(0.5*T(i)).*cos(0.5*T(i));
Z=exp(-R.^2/64).*cos(-9*exp(-R/30+3)-2*angle(Y-X*1i))/pi.^0.5;
n=.08*randn(1,numel(i));
scatter3(X+n,Y,Z+n,1,abs(Z+n))
axis equal off
colormap jet
view(2)