\begin{Schunk}
\begin{Sinput}
 par(mfcol = c(1, 2))
 x <- 0:10
 y <- dbinom(x, p = 0.66, size = 10)
 plot(x, y, type = "h", lty = 4, xlab = "10 Flips with a Biased Coin", ylab = "Chance of Observing x Heads")
 points(x, y, pch = 16)
 y <- c(y[1], y)
 x <- c(-1, x)
 plot(x + 0.5, y, type = "s", lty = 4, xlab = "10 Flips with a Biased Coin", 
      ylab = "Chance of Observing x Heads")
 par(mfcol = c(1, 1))
\end{Sinput}
\end{Schunk}
