library(ggplot2)
bws <- exp(seq(log(0.1), log(2), length.out = 50))

for (bw in bws)
    plot(density(MASS::geyser$duration, bw = bw), ylim = c(0, 1.0))