bernoulli.py

#!/usr/bin/env python2
# -*- coding: utf-8 -*-

from scipy.special import binom

def bernoulli(p, n, k):
    return binom(n, k) * (p ** k) * ((1.0 - p) ** (n - k))



def main():
    p = 1.0/50.0      #prawdopodobienstwo wylosowania tego elementu co trzeba w jednym losowaniu
    n = 1000          #liczba losowan
    k = 10            #ilosc potrzebnych takich samych elementow
    
    print bernoulli(p, n, k)



if __name__ == "__main__":
    main()