Beta-Binomialfunktion in Python

Ich möchte die Wahrscheinlichkeit berechnen, die durch eine Binomialverteilung für vorgegebene x (Erfolge), n (Versuche) und p (Wahrscheinlichkeit) gegeben ist - deren letzterer durch eine Wahrscheinlichkeitsmassenfunktion Beta (a, b) gegeben ...