import matplotlib.pyplot as plt
import math
import numpy as np

def f(n):
    fact=1
    if n != 0:
        for i in range(1,n+1):
            fact=fact*i
    return fact
        
def k_parmi_n(k,n):
    return f(n)/(f(k)*f(n-k))
    
def proba_binom(n,p,k):
    return k_parmi_n(k,n)*p**k*(1-p)**(n-k)
    

def graph(k,p):
    plt.close()
    X=[n for n in range(10,k+1)]
    plt.plot(X,p2(k,p),'r.')
    plt.plot([10,k],[p*(1-p),p*(1-p)],'b-')
    plt.show()