par(mfrow=c(2,2))

n<-25 # ½ÃÇàȽ¼ö

p_list<-c(0.1, 0.2, 0.5, 0.7) # ¹ß»ýÈ®·ü

for (i in 1:length(p_list)) {
    p_x<-dbinom(x=0:n, n, p_list[i])
    plot(x=0:n, p_x, xlab="x", ylab="P(X=x)", 
         ylim=c(0, 0.3), xlim=c(0,n), main=paste("p=", p_list[i]))
    x_seq<-seq(0,n,1)
    lines(x_seq, p_x, type="h", col="blue")
}