\begin{Schunk}
\begin{Sinput}
 N <- 437
 p <- 0.63
 x1 <- max(0, N * p - 4 * sqrt(p * (1 - p) * N))
 x2 <- min(N * p + 4 * sqrt(p * (1 - p) * N), N)
 x <- as.integer(x1):as.integer(x2 + 1)
 pseq <- dbinom(x, N, p)
 plot(x, pseq, type = "h", xlab = "k", ylab = paste("Prob(k, N=", N, ", p=", 
      p, ")"))
 points(x, pseq, pch = 18, cex = 0.5)
\end{Sinput}
\end{Schunk}
