Lorsque vous utiliserez le raccourci =C-c C-e h o=, ce document sera
Lorsque vous utiliserez le raccourci =C-c C-e h o=, ce document sera
compilé en html. Tout le code contenu sera ré-exécuté, les résultats
compilé en html. Tout le code contenu sera ré-exécuté, les résultats
* En utilisant la méthode des aiguilles de Buffon
récupérés et inclus dans un document final. Si vous ne souhaitez pas
récupérés et inclus dans un document final. Si vous ne souhaitez pas
Mais calculé avec la *méthode* des [[https://fr.wikipedia.org/wiki/Aiguille_de_Buffon][aiguilles de Buffon]], on obtiendrait
ré-exécuter tout le code à chaque fois, il vous suffit de supprimer
ré-exécuter tout le code à chaque fois, il vous suffit de supprimer
comme *approximation* :
le # et l'espace qui sont devant le ~#+PROPERTY:~ au début de ce
le # et l'espace qui sont devant le ~#+PROPERTY:~ au début de ce
document.
document.
#+begin_src R :results output :session *R* :exports both
Comme nous vous l'avons montré dans la vidéo, on inclue du code
set.seed(42)
python de la façon suivante (et on l'exécute en faisant ~C-c C-c~):
Comme nous vous l'avons montré dans la vidéo, on inclut du code
N = 100000
#+begin_src python :results output :exports both
R de la façon suivante (et on l'exécute en faisant ~C-c C-c~):
x = runif(N)
theta = pi/2*runif(N)
#+begin_src R :results output :exports both
2/(mean(x+sin(theta)>1))
print("Hello world!")
print("Hello world!")
#+end_src
#+end_src
#+end_src
#+RESULTS:
#+RESULTS:
: Hello world!
#+RESULTS:
: [1] 3.14327
Voici la même chose, mais avec une session python, donc une
: [1] "Hello world!"
* Avec un argument "fréquentiel" de surface
Voici la même chose, mais avec une session R (c'est le cas le
Sinon, une méthode plus simple à comprendre et ne faisant pas
plus courant, R étant vraiment un langage interactif), donc une
intervenir d'appel à la fonction sinus se base sur le fait que si $X\sim
persistance d'un bloc à l'autre (et on l'exécute toujours en faisant
persistance d'un bloc à l'autre (et on l'exécute toujours en faisant
U(0,1)$ et $Y\sim U(0,1)$ alors $P[X^2+Y^2\leq 1] = \pi/4$ (voir [[https://fr.wikipedia.org/wiki/M%25C3%25A9thode_de_Monte-Carlo#D%25C3%25A9termination_de_la_valeur_de_%25CF%2580][méthode de
~C-c C-c~).
~C-c C-c~).
#+begin_src python :results output :session :exports both
Monte Carlo sur Wikipedia]]). Le code suivant illustre ce fait :
import numpy
#+begin_src R :results output :session *R* :exports both
x=numpy.linspace(-15,15)
#+begin_src R :results output graphics :file figure_pi_mc1.png :exports both :width 600 :height 400 :session *R*
