Binomijakauma

# Binomijakauma 1000 kolikkoa heitetään
# Mikä on todennäköisyys, että klaavoja vähintään 495, enintään 510
# Juhani Kaukoranta 22.5.2018

from decimal import *
import math


def binomi(n,k):
    a = math.factorial(n)//(math.factorial(n-k)*math.factorial(k))
    return(a)

tulos = 0
for k in range(495,511):
    tulos = tulos+0.5**k*0.5**(1000-k)*binomi(1000, k)

print("Tulos = ",tulos)