support <- 0:40
plot(support+3, dnbinom(support,3,0.5), type="l")
points(support+3, dnbinom(support,3,0.3), type="l")
points(support+3, dnbinom(support,3,0.1), type="l")