import random
import matplotlib.pyplot as plt


def prop(i,j):
    return i/(i+j)
    
eff_total=int(input("Effectif population ? "))  
eff_jv=eff_total*0.7
eff_pas_jv=eff_total*0.3
p=prop(eff_jv,eff_pas_jv)

liste_proba=[]

for i in range(20):
    liste_proba.append(p)
    if random.random() <= p:
        eff_jv=eff_jv-1
    else:
        eff_pas_jv=eff_pas_jv-1
    p=prop(eff_jv,eff_pas_jv)
    
    
plt.plot([i for i in range(20)],liste_proba,'b.')
plt.show()