R negative binomial distribution

> par(mfrow=c(5,2))
> for (i in (1:10)/10) {
+   barplot(dbinom(0:10, 10, i), main=sprintf('prob=%s',i))
+ }