This is just the Gamma distribution with shape parameter α=n an integer.
Implementation: scipy.stats.erlang