from random import*
proba_initiale_A = 0.3
proba_initiale_B = 1 - proba_initiale_A
pAA = 0.9
pAB = 0.1
pBA = 0.4
pBB = 0.6
def etat_suivant(etat):
if etat == "A":
if random() < pAB:
return "B"
else:
return "A"
if etat == "B":
if random() < pBA:
return "A"
else:
return "B"
def etat_asymptotique():
if random() < proba_initiale_A:
etat = "A"
else:
etat = "B"
for i in range(1000):
etat = etat_suivant(etat)
return etat
def simule_distribution(n):
total = 0
for i in range(n):
if etat_asymptotique() == "A":
total = total + 1
resultat = total/n
return resultat
Nos manuels sont collaboratifs, n'hésitez pas à nous en faire part.
j'ai une idée !
Oups, une coquille