N <- c(0:20)
P <- data.frame( Prob = choose(20, N)* 0.15^(20-N) *0.85^N )
P$N <- round(2000*P$Prob)

#画图：
barplot(P$N,col = "red")

plot(P$Prob, type = "b")

什么意思啊，里面的参数不是很懂啊

加了注释你再看看吧，不知道你到底要哪样的图。

谢谢 爱迪生_陈！借花献佛，稍微改造了一下程序。
m = 20  # 钉板层数，也是最底层钉子数
p = 0.15  # 落在钉子左侧的概率
n = 5000  # 玻璃球数
space <- c(0:m)  # 最下面21个格子
# 计算掉到每个格子的概率，排列问题
result <- data.frame(Prob = choose(m, space)* p^(m-space) *(1-p)^ space)
# 四舍五入，5000个球掉入每个格子球数的数学期望
result$Num <- round(n*result$Prob)
            
#画图：
barplot(result$Num, col = "red")  # 条形图
plot(result$Prob, type = "b")     # 点线图