g = Table[
  ListPlot[
    Table[{k, Binomial[k + r - 1, r - 1] (r/(10 + r))^r (1 - r/(10 + r))^k}, {k, 0, 25}],
    PlotRange -> 0.118, 
    Filling -> Axis,
    PlotStyle -> PointSize[0.02], 
    FillingStyle -> Directive[Hue[0.6], Thickness[0.01], Opacity[0.3]],
    Epilog -> Inset[Style["\!\(\*StyleBox[\"r\",\nFontSlant->\"Italic\"]\) = " <> ToString[r], 16], {23, 0.08}], 
    Prolog -> {Opacity[0.5], Green, Thin, Line[{{10 - Sqrt[(10 (10 + r))/r], .1}, {10 + Sqrt[(10 (10 + r))/r], .1}}], 
               Dashed, Thick, Orange, Line[{{10, 0}, {10, 1}}]},
    ImageSize -> 200,
    LabelStyle -> Directive[Larger],
    AxesLabel -> {"\!\(\*StyleBox[\"k\",\nFontSlant->\"Italic\"]\)"},
    AxesOrigin -> {-1, 0}
    ], 
  {r, {1, 2, 3, 4, 5, 10, 20, 40}}
  ]
Export["c:/users/stpasha/documents/negbinomial.gif", g, "DisplayDurations" -> {.5, .5, .5, .5, .5, .5, .5, .15}, "TransparentColor" -> White, ImageSize -> 250]

我，此作品的版權所有人，釋出此作品至公共領域。此授權條款在全世界均適用。 這可能在某些國家不合法，如果是的話： 我授予任何人有權利使用此作品於任何用途，除受法律約束外，不受任何限制。