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