Commit fe7b6ca5 authored by Olivia Guillin's avatar Olivia Guillin

Avec_la_figure

parent 9ab8e45f
#+TITLE: À propos du calcul de $\pi$ ho#+TITLE: À propos du calcul de $\pi$
#+LANGUAGE: fr #+LANGUAGE: fr
...@@ -41,8 +41,7 @@ Sinon une méthode plus simple à comprendre et ne faisant pas ...@@ -41,8 +41,7 @@ Sinon une méthode plus simple à comprendre et ne faisant pas
intervenir d'appel à la fonction sinus se base sur le fait que si 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 \le1] = \pi/4$ (voir [[https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Monte-Carlo#D%C3%A9termination_de_la_valeur_de_%CF%80][méthode $X\sim U(0,1)$ et $Y\sim U(0,1)$ alors $P[X^2+Y^2 \le1] = \pi/4$ (voir [[https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Monte-Carlo#D%C3%A9termination_de_la_valeur_de_%CF%80][méthode
de Monte Carlo sur wikipedia]]. Le code suivant illustre ce fait : de Monte Carlo sur wikipedia]]. Le code suivant illustre ce fait :
#+begin_src python :results file :session :var matplot_lib_filename=(org-babel-temp-file "figure" ".png") :exports both
#+begin_src python :results output file :var matplot_lib_filename="figure_pi_mc2.png" :exports both :session *python*
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
np.random.seed(seed=42) np.random.seed(seed=42)
...@@ -59,15 +58,11 @@ ax.scatter(x[reject], y[reject], c='r', alpha=0.2, edgecolor=None) ...@@ -59,15 +58,11 @@ ax.scatter(x[reject], y[reject], c='r', alpha=0.2, edgecolor=None)
ax.set_aspect('equal') ax.set_aspect('equal')
plt.savefig(matplot_lib_filename) plt.savefig(matplot_lib_filename)
print(matplot_lib_filename) matplot_lib_filename
#+end_src #+end_src
#+RESULTS: #+RESULTS:
[[file:Traceback (most recent call last): [[file:/var/folders/21/wtvjlyxj21b3gj703l4k0j_40000gp/T/babel-0OJhDQ/figureX0cFtl.png]]
File "<stdin>", line 1, in <module>
File "/var/folders/21/wtvjlyxj21b3gj703l4k0j_40000gp/T/babel-4uvhkJ/python-bYwGgQ", line 4, in <module>
np.random.seed(seed=42)
NameError: name 'np' is not defined]]
Il est alors aisé d'obtenir une approximation (pas terrible) de $\pi$ en Il est alors aisé d'obtenir une approximation (pas terrible) de $\pi$ en
comptant combien de fois, en moyenne, $X^2+Y^2$ est inférieur à 1 : comptant combien de fois, en moyenne, $X^2+Y^2$ est inférieur à 1 :
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment