следует также посмотреть на , что эта ссылка
делает два списка для k, например xk и yk
xk
yk
from scipy import stats custm = stats.rv_discrete(name='test', values=(xk, yk)) custm.rvs(size=1)