Commit 5f910f09 authored by Lucie Leonarski's avatar Lucie Leonarski

Mon deuxième commit magit pour l'exercice 2

parent 0a39ce05
...@@ -13,13 +13,14 @@ ...@@ -13,13 +13,14 @@
* En demandant à la lib maths * En demandant à la lib maths
Mon ordinateur m'indique que \pi vaut approximativement: Mon ordinateur m'indique que $\pi$ vaut approximativement:
#+begin_src python :results output :session :exports both #+begin_src python :results output :session :exports both
from math import * from math import *
pi pi
#+end_src #+end_src
#+RESULTS:
* En utilisant la méthode des aiguilles de Buffon * En utilisant la méthode des aiguilles de Buffon
...@@ -34,10 +35,11 @@ theta = np.random.uniform(size=N, low=0, high=pi/2) ...@@ -34,10 +35,11 @@ theta = np.random.uniform(size=N, low=0, high=pi/2)
2/(sum((x+np.sin(theta))>1)/N) 2/(sum((x+np.sin(theta))>1)/N)
#+end_src #+end_src
#+RESULTS:
* Avec un argument "fréquentiel" de surface * Avec un argument "fréquentiel" de surface
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 ∼ U(0,1) et Y ∼ U(0,1) alors P[X^2+Y^2 ≤1]=\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 : 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 ≤1]=\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 :
#+begin_src python :results output file :session :var matplot_lib_filename="./random_uniform_pi.png" :exports results #+begin_src python :results output file :session :var matplot_lib_filename="./random_uniform_pi.png" :exports results
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
......
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