For m = 1…2000 and n = 1…1200 the pixel in row n, column m gets the colour below. Every symbol is a function of the point (x, y) unless it only depends on s. This page evaluates exactly these lines, then uses the lobe functions Ls to give the picture depth.
PixelScreen coordinates map to the plane with 420 pixels per unit.
\[\operatorname{rgb}\!\Big(F\big(H_0(\tfrac{m-1000}{420},\tfrac{601-n}{420})\big),\;F\big(H_1(\cdot)\big),\;F\big(H_2(\cdot)\big)\Big)\]
ClampA smooth way of writing ⌊255·min(max(x,0),1)⌋.
\[F(x)=\Big\lfloor 255\,e^{-e^{-1000x}}\,|x|^{e^{-e^{1000(x-1)}}}\Big\rfloor\]
ColourNebula, core glow and stars, one channel per v ∈ {0, 1, 2}.
\[H_v(x,y)=\tfrac{11}{10}\big(1-W\big)K_vA+\tfrac{4+v^2-v}{2}\,W+T_v\]
Nebula50 filament layers, each with its own tint.
\[K_v=\sum_{s=1}^{50}I_s\Big(\tfrac{19}{20}\Big)^{s}\frac{12-4v+v^2+(v-1)\cos(2s^2)+8\cos\big((7+v)s^2\big)}{50}\]
Filaments
\[I_s=45\,C_{1,s}+6\,C_{0,s}\]\[C_{v,s}=e^{-e^{-\frac{1+15v}{4}\left(\cos\left(\frac15\left(\frac{23}{20}\right)^{s}\left(\cos(15s^2)S+\sin(15s^2)Q_s\right)+2\cos(27s^2)\right)\cos\left(\frac15\left(\frac{23}{20}\right)^{s}\left(\cos(15s^2)Q_s-\sin(15s^2)S\right)+2\cos(28s^2)\right)-\frac54+2A+\frac{E}{7}\right)}}\]
Stars30 rotated lattices; warm and blue tints alternate with s.
\[T_v=\sum_{s=1}^{30}\frac{v^2-2v+4+(v-1)(-1)^s}{4}\Big(4e^{-e^{200\left(M_s^2+N_s^2-\frac{1}{800}-\frac{B_s}{200}\right)}}+e^{-e^{20M_s^2+20N_s^2-\frac{7}{50}}}\Big)\]\[B_s=e^{-e^{5\cos\left(20\arctan\frac{M_s}{N_s}+2\cos\left(9\arctan\frac{M_s}{N_s}+s^2\right)\right)+\frac{15}{4}}}\]\[M_s=O\big(6^s5^{-s}2(\cos(19s^2)P_s+\sin(19s^2)Q_s)+2\cos(27s^2)\big)\]\[N_s=O\big(6^s5^{-s}2(\cos(19s^2)Q_s-\sin(19s^2)P_s)+2\cos(28s^2)\big)\]
Shell lightBrightest just inside each lobe wall.
\[A=\frac14\sum_{s=1}^{27}\Big(\prod_{u=0}^{s-1}1-J_u\Big)J_s\,e^{-e^{\frac{3}{20}(s-23)}-e^{-3L_s}}\]
Core & noise
\[O(x)=\arccos(\cos x),\qquad W=e^{-e^{10\sqrt{x^2+y^2}-1+\frac{E}{4}}},\qquad E=\sum_{s=1}^{50}\Big(\tfrac{19}{20}\Big)^{s}D_s\]\[\begin{aligned}D_s={}&\cos\!\Big(5^s4^{-s}\big(\cos(7s)S+\sin(7s)Q_s+2\cos(17s)\big)+4\cos\!\big(5^s4^{-s}(\cos(4s)S+\sin(4s)Q_s)\big)+2\cos(5s)\Big)\\&\times\cos\!\Big(5^s4^{-s}\big(\cos(7s)Q_s-\sin(7s)S+2\cos(15s)\big)+4\cos\!\big(5^s4^{-s}(\cos(8s)S+\sin(8s)Q_s)\big)+2\cos(7s)\Big)\end{aligned}\]
Lobes27 nested bipolar shapes; S is a signed distance to the first one containing the point.
\[S=2\sum_{s=1}^{27}\Big(\prod_{u=0}^{s-1}1-J_u\Big)J_sL_s,\qquad J_s=e^{-e^{25-50s}-e^{10L_s}}\]\[L_s=\sqrt{U_s^2+\Big(2R(s)^{\frac{3}{10}}|U_s|^{-\frac{3}{10}}\Big(y-\big(\tfrac{3}{20}+\tfrac{\cos(4s^2)}{5}\big)x\Big)\Big)^2}-R(s)\]
Frames
\[P_s=y\cos(15s^2)-x\sin(15s^2),\qquad Q_s=x\cos(15s^2)+y\sin(15s^2)\]\[R(s)=\tfrac{s}{10}+\tfrac{3}{50}\cos(5s^2),\qquad U_s=x+\big(\tfrac{3}{20}+\tfrac{\cos(3s^2)}{5}\big)y+10^{-4}\]