from random import random from matplotlib.pyplot import * def OM2(x,y): return x**2+y**2 def KM2(x,y): return (x-1)**2+(y-1)**2 def simul(n): NbSuccès=0 for i in range(n): x=random() y=random() if OM2(x,y)<=1 and KM2(x,y)<=1: NbSuccès=NbSuccès+1 return NbSuccès/n def graphique(N,n): plot([0,N],[0.5,0.5],color='black') Lx=list(range(1,N+1)) Ly=[] for i in range(N): f=simul(n) Ly.append(f) scatter(Lx,Ly,color='red') axis((-0.05*N,N+0.05*N,-0.05,1.05)) show() graphique(200,100)