"Mais calculé avec la **méthode** des (aiguilles de Buffon)[https://fr.wikipedia.org/wiki/Aiguille_de_Buffon], on obitendrait comme **approximation** :"
"Mais calculé avec la **méthode** des [aiguilles de Buffon](https://fr.wikipedia.org/wiki/Aiguille_de_Buffon), on obitendrait comme **approximation** :"
]
]
},
},
{
{
...
@@ -78,6 +78,82 @@
...
@@ -78,6 +78,82 @@
"2/(sum((x+np.sin(theta))>1)/N)"
"2/(sum((x+np.sin(theta))>1)/N)"
]
]
},
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Avec un arguement \"fréquentiel\" de surface"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Sinon, une méthode plus simple à comprendre et ne faisant pas intervenir d'appel à la fonction sinus se base sur le fait que si $X\\sim U(0,1)$ et $Y\\sim U(0,1)$ alors $P[X^2 + Y^2 \\leq 1] = \\pi/4$ (voir [méthode de Monte Carlo sur Wikipédia](https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Monte-Carlo#D%C3%A9termination_de_la_valeur_de_%CF%80). Le code suivant illustre ce fait :"
