Code
#reset plot to defaults
.pardefault <- par()
# generate data
theta = 1:100
x = sin(theta)
y = cos(theta)
#create color scheme
color_list=c("#08519c", "#3182bd", "#6baed6", "#bdd7e7")
# set up graphic parameters & open plot
op = par(bg = 'black', mar = rep(0.5, 4))
plot.new()
plot.window(xlim = c(-1, 1), ylim = c(-1, 1), asp = 1)
# add circles one at a time
lines(x, y, col=color_list[1])
lines(0.8 * x, 0.8 * y, col=color_list[2])
lines(0.6 * x, 0.6 * y, col=color_list[3])
lines(0.4 * x, 0.4 * y, col=color_list[4])
