from math import sqrt from random import random from matplotlib.pyplot import * def succes(p): if random()<p: res=True else: res=False return res def frequence_succes(n,p): nb_succes=0 for i in range(n): if succes(p): nb_succes+=1 f_succes=nb_succes/n return f_succes def graph_sim(N,n,p): axis((-0.1*N,1.1*N,-0.1,1.1)) text(0,1.1,"N={} et n={}".format(N,n)) plot([-0.1*N,1.1*N],[p,p], color="red") plot([-0.1*N,1.1*N], [p-1/sqrt(n),p-1/sqrt(n)], color="green") plot([-0.1*N,1.1*N], [p+1/sqrt(n),p+1/sqrt(n)], color="green") Lx=list(range(1,N+1)) Ly=[] for i in range(N): Ly.append(frequence_succes(n,p)) scatter(Lx,Ly) show()