Shell

# initialize
n <- 25
r <- (1:n)^1.5
x <- r*cos(seq(-pi, pi, length=n) )
y <- r*sin(seq(-pi, pi, length=n) )
color <- rgb(0.133, 1, 0.133, alpha=0.25)

# plot
par(bg="black", mar=c(1, 1, 1, 1))
plot(x, y, type="n", asp=1)
for (i in 1:n) for (j in 1:n) lines(x[c(i, j)], y[c(i, j)], col=color)