# 二項分布

# パラメータ
n <- 10
p <- 0.5

# 確率密度関数
x <- 0 : n
f <- dbinom(x, n, p)

# グラフ
title <- sprintf("Bi(%d,%g)", n, p)
barplot(f, names.arg = x, xlab = "x", ylab = "f", space = 0, main = title)

# 値の出力
print(sprintf("f(%d)=%f", x, f), quote = FALSE)