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.
Oups, une coquille
j'ai une idée !