- /
-
City lights at night
on 23 Oct 2021
- 0
- 12
- 0
- 0
- 246
col=20;
m=300;
cx=0;
cy=0;
l=1;
x=linspace(cx-l,cx+l,m);
y=linspace(cy-l,cy+l,m);
[X,Y]=meshgrid(x,y);
Z=X +i*Y;
R=Z;
c=-0.5-i*0.8660254;
for k=1:col;
Z=2/3*Z + 1/3*1./(eps+Z.^2);
end
W=abs(Z-c);
A=angle(Z);
colormap copper(256);
mesh(W-A);
shading flat;
axis('square','equal','off');