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))