§ Как галактика выглядит

Я искал этот алгоритм годами и найти не мог. Сегодня я смог его найти и немного видоизменил, чтобы красиво выглядело.

§ Исходящий код

1SCREEN 13
2
3' Palette
4FOR i = 1 TO 63: OUT 968, i: OUT 969, i: OUT 969, i: OUT 969, i: NEXT
5
6FOR i = 0 TO 80000
7
8  ' Logaritmics distribution
9   s = -LOG(RND(1))
10
11  ' Spiro Ellipsis
12  t = s * 200
13  u = s * 20 * SIN(t)
14  v = s * 25 * COS(t)
15  p = s * 3
16
17  ' Precalculate
18  c = COS(p): s = SIN(p)
19
20  ' Rotate (x,y) to P angle, r-dist
21  x = u * c + v * s
22  y = v * c - u * s
23  r = x * x + y * y
24
25  ' Bald
26  z = 40 * EXP(-r / 1000) * (RND(1) - .5)
27  x = x + RND(1) * 5 - z * .7
28  y = y + RND(1) * 5 + z
29
30  ' Rotation
31  x = x * 1.4 + y * .6
32  y = x * .2 + y * .8
33
34  ' Set
35  x = x + 160
36  y = y + 100
37
38  IF ABS(x) < 320 AND ABS(y) < 200 THEN
39    k = POINT(x, y) + 1
40    IF k > 63 THEN k = 63
41    PSET (x, y), k
42  END IF
43
44NEXT i