Hi all,

happy new year!

I noted that the Planck distribution is a geometric distribution with a different parametrization, see Issue #9359:

import numpy as np
from scipy.stats import planck, geom

a = 0.5
k = np.arange(20)
sum(abs(geom.pmf(k, 1-np.exp(-a), loc=-1) - planck.pmf(k, a))) # 1.30e-18

I don't know if there is a specific reason to have the Planck distribution in addition to the geometric. If not, I would propose to deprecate it.

Any views? Thanks

Christoph