coins <- dbinom(0:8, 8, 1/2)
sum(coins)
## [1] 1
plot(0:8, coins, type = "h")
par(mfrow = c(5, 2))
for (i in 1:10) plot(0:i, dbinom(0:i, i, 1/2), type = "h", xlab = paste("0:",
i, sep = ""), ylab = paste("dbinom(0:", i, ",", i, ",", "1/2", ")", sep = ""))