Bradford Distribution

There is one shape parameter, $c>0$, and the support is $x\in [0,1]$. Define $k=\log \left(1+c\right)$.

where $\mathrm{Ei}\left(z\right)$ is the exponential integral function.

Implementation: scipy.stats.bradford