\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}
