# =============================================================
# Blue feathers
# =============================================================
# generate pairs of x-y values
x1 <- c(seq(0, pi, length = 50), seq(pi, 2*pi, length = 50))
y1 <- cos(x1) / sin(x1)
x2 <- seq(1.02 * 2 * pi + pi/2, 4*pi + pi/2, length = 50)
y2 <- tan(x2)
png("blue_feathers.png", width = 700, height = 400)
# set graphical parameters
op <- par(bg="black", mar=rep(.5,4))
# plot
plot(c(x1, x2), c(y1, y2), type = "n", ylim = c(-11, 11))
for (i in seq(-10, 10, length = 100))
{
lines(x1, y1 + i, col = hsv(runif(1, 0.65, 0.7), 1, 1, runif(1, 0.7)),
lwd = 4 * runif(1, 0.3))
lines(x2, y2 + i, col = hsv(runif(1, 0.65, 0.7), 1, 1, runif(1, 0.7)),
lwd = 4 * runif(1, 0.3))
}
# signature
legend("bottomright", legend = "© Gaston Sanchez", bty = "n",
text.col = "gray70")
# reset par
par(op)
dev.off()