Skip to content
Projects
Groups
Snippets
Help
Loading...
Help
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
M
mooc-rr
Project
Project
Details
Activity
Releases
Cycle Analytics
Repository
Repository
Files
Commits
Branches
Tags
Contributors
Graph
Compare
Charts
Issues
0
Issues
0
List
Board
Labels
Milestones
Merge Requests
0
Merge Requests
0
CI / CD
CI / CD
Pipelines
Jobs
Schedules
Charts
Wiki
Wiki
Snippets
Snippets
Members
Members
Collapse sidebar
Close sidebar
Activity
Graph
Charts
Create a new issue
Jobs
Commits
Issue Boards
Open sidebar
4370a735ad7d84eb72e83b8c07159a69
mooc-rr
Commits
2c85c94a
Commit
2c85c94a
authored
Jul 29, 2020
by
Pierre Ayoub
Browse files
Options
Browse Files
Download
Email Patches
Plain Diff
[MODULE2/EXO1] Achieved
parent
95d00c65
Changes
3
Hide whitespace changes
Inline
Side-by-side
Showing
3 changed files
with
449 additions
and
83 deletions
+449
-83
pi_monte_carlo.png
module2/exo1/pi_monte_carlo.png
+0
-0
toy_document_orgmode_python_fr.html
module2/exo1/toy_document_orgmode_python_fr.html
+377
-0
toy_document_orgmode_python_fr.org
module2/exo1/toy_document_orgmode_python_fr.org
+72
-83
No files found.
module2/exo1/pi_monte_carlo.png
0 → 100644
View file @
2c85c94a
79.1 KB
module2/exo1/toy_document_orgmode_python_fr.html
0 → 100644
View file @
2c85c94a
<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN"
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html
xmlns=
"http://www.w3.org/1999/xhtml"
lang=
"fr"
xml:lang=
"fr"
>
<head>
<!-- 2020-07-29 mer. 14:28 -->
<meta
http-equiv=
"Content-Type"
content=
"text/html;charset=utf-8"
/>
<meta
name=
"viewport"
content=
"width=device-width, initial-scale=1"
/>
<title>
À propos du calcul de \(\pi\)
</title>
<meta
name=
"generator"
content=
"Org mode"
/>
<meta
name=
"author"
content=
"Pierre AYOUB"
/>
<style
type=
"text/css"
>
<!
--
/*
--
><![
CDATA
[/*><!
--
*/
.title
{
text-align
:
center
;
margin-bottom
:
.2em
;
}
.subtitle
{
text-align
:
center
;
font-size
:
medium
;
font-weight
:
bold
;
margin-top
:
0
;
}
.todo
{
font-family
:
monospace
;
color
:
red
;
}
.done
{
font-family
:
monospace
;
color
:
green
;
}
.priority
{
font-family
:
monospace
;
color
:
orange
;
}
.tag
{
background-color
:
#eee
;
font-family
:
monospace
;
padding
:
2px
;
font-size
:
80%
;
font-weight
:
normal
;
}
.timestamp
{
color
:
#bebebe
;
}
.timestamp-kwd
{
color
:
#5f9ea0
;
}
.org-right
{
margin-left
:
auto
;
margin-right
:
0px
;
text-align
:
right
;
}
.org-left
{
margin-left
:
0px
;
margin-right
:
auto
;
text-align
:
left
;
}
.org-center
{
margin-left
:
auto
;
margin-right
:
auto
;
text-align
:
center
;
}
.underline
{
text-decoration
:
underline
;
}
#postamble
p
,
#preamble
p
{
font-size
:
90%
;
margin
:
.2em
;
}
p
.verse
{
margin-left
:
3%
;
}
pre
{
border
:
1px
solid
#ccc
;
box-shadow
:
3px
3px
3px
#eee
;
padding
:
8pt
;
font-family
:
monospace
;
overflow
:
auto
;
margin
:
1.2em
;
}
pre
.src
{
position
:
relative
;
overflow
:
visible
;
padding-top
:
1.2em
;
}
pre
.src
:before
{
display
:
none
;
position
:
absolute
;
background-color
:
white
;
top
:
-10px
;
right
:
10px
;
padding
:
3px
;
border
:
1px
solid
black
;
}
pre
.src
:hover:before
{
display
:
inline
;}
/* Languages per Org manual */
pre
.src-asymptote
:before
{
content
:
'Asymptote'
;
}
pre
.src-awk
:before
{
content
:
'Awk'
;
}
pre
.src-C
:before
{
content
:
'C'
;
}
/* pre.src-C++ doesn't work in CSS */
pre
.src-clojure
:before
{
content
:
'Clojure'
;
}
pre
.src-css
:before
{
content
:
'CSS'
;
}
pre
.src-D
:before
{
content
:
'D'
;
}
pre
.src-ditaa
:before
{
content
:
'ditaa'
;
}
pre
.src-dot
:before
{
content
:
'Graphviz'
;
}
pre
.src-calc
:before
{
content
:
'Emacs Calc'
;
}
pre
.src-emacs-lisp
:before
{
content
:
'Emacs Lisp'
;
}
pre
.src-fortran
:before
{
content
:
'Fortran'
;
}
pre
.src-gnuplot
:before
{
content
:
'gnuplot'
;
}
pre
.src-haskell
:before
{
content
:
'Haskell'
;
}
pre
.src-hledger
:before
{
content
:
'hledger'
;
}
pre
.src-java
:before
{
content
:
'Java'
;
}
pre
.src-js
:before
{
content
:
'Javascript'
;
}
pre
.src-latex
:before
{
content
:
'LaTeX'
;
}
pre
.src-ledger
:before
{
content
:
'Ledger'
;
}
pre
.src-lisp
:before
{
content
:
'Lisp'
;
}
pre
.src-lilypond
:before
{
content
:
'Lilypond'
;
}
pre
.src-lua
:before
{
content
:
'Lua'
;
}
pre
.src-matlab
:before
{
content
:
'MATLAB'
;
}
pre
.src-mscgen
:before
{
content
:
'Mscgen'
;
}
pre
.src-ocaml
:before
{
content
:
'Objective Caml'
;
}
pre
.src-octave
:before
{
content
:
'Octave'
;
}
pre
.src-org
:before
{
content
:
'Org mode'
;
}
pre
.src-oz
:before
{
content
:
'OZ'
;
}
pre
.src-plantuml
:before
{
content
:
'Plantuml'
;
}
pre
.src-processing
:before
{
content
:
'Processing.js'
;
}
pre
.src-python
:before
{
content
:
'Python'
;
}
pre
.src-R
:before
{
content
:
'R'
;
}
pre
.src-ruby
:before
{
content
:
'Ruby'
;
}
pre
.src-sass
:before
{
content
:
'Sass'
;
}
pre
.src-scheme
:before
{
content
:
'Scheme'
;
}
pre
.src-screen
:before
{
content
:
'Gnu Screen'
;
}
pre
.src-sed
:before
{
content
:
'Sed'
;
}
pre
.src-sh
:before
{
content
:
'shell'
;
}
pre
.src-sql
:before
{
content
:
'SQL'
;
}
pre
.src-sqlite
:before
{
content
:
'SQLite'
;
}
/* additional languages in org.el's org-babel-load-languages alist */
pre
.src-forth
:before
{
content
:
'Forth'
;
}
pre
.src-io
:before
{
content
:
'IO'
;
}
pre
.src-J
:before
{
content
:
'J'
;
}
pre
.src-makefile
:before
{
content
:
'Makefile'
;
}
pre
.src-maxima
:before
{
content
:
'Maxima'
;
}
pre
.src-perl
:before
{
content
:
'Perl'
;
}
pre
.src-picolisp
:before
{
content
:
'Pico Lisp'
;
}
pre
.src-scala
:before
{
content
:
'Scala'
;
}
pre
.src-shell
:before
{
content
:
'Shell Script'
;
}
pre
.src-ebnf2ps
:before
{
content
:
'ebfn2ps'
;
}
/* additional language identifiers per "defun org-babel-execute"
in ob-*.el */
pre
.src-cpp
:before
{
content
:
'C++'
;
}
pre
.src-abc
:before
{
content
:
'ABC'
;
}
pre
.src-coq
:before
{
content
:
'Coq'
;
}
pre
.src-groovy
:before
{
content
:
'Groovy'
;
}
/* additional language identifiers from org-babel-shell-names in
ob-shell.el: ob-shell is the only babel language using a lambda to put
the execution function name together. */
pre
.src-bash
:before
{
content
:
'bash'
;
}
pre
.src-csh
:before
{
content
:
'csh'
;
}
pre
.src-ash
:before
{
content
:
'ash'
;
}
pre
.src-dash
:before
{
content
:
'dash'
;
}
pre
.src-ksh
:before
{
content
:
'ksh'
;
}
pre
.src-mksh
:before
{
content
:
'mksh'
;
}
pre
.src-posh
:before
{
content
:
'posh'
;
}
/* Additional Emacs modes also supported by the LaTeX listings package */
pre
.src-ada
:before
{
content
:
'Ada'
;
}
pre
.src-asm
:before
{
content
:
'Assembler'
;
}
pre
.src-caml
:before
{
content
:
'Caml'
;
}
pre
.src-delphi
:before
{
content
:
'Delphi'
;
}
pre
.src-html
:before
{
content
:
'HTML'
;
}
pre
.src-idl
:before
{
content
:
'IDL'
;
}
pre
.src-mercury
:before
{
content
:
'Mercury'
;
}
pre
.src-metapost
:before
{
content
:
'MetaPost'
;
}
pre
.src-modula-2
:before
{
content
:
'Modula-2'
;
}
pre
.src-pascal
:before
{
content
:
'Pascal'
;
}
pre
.src-ps
:before
{
content
:
'PostScript'
;
}
pre
.src-prolog
:before
{
content
:
'Prolog'
;
}
pre
.src-simula
:before
{
content
:
'Simula'
;
}
pre
.src-tcl
:before
{
content
:
'tcl'
;
}
pre
.src-tex
:before
{
content
:
'TeX'
;
}
pre
.src-plain-tex
:before
{
content
:
'Plain TeX'
;
}
pre
.src-verilog
:before
{
content
:
'Verilog'
;
}
pre
.src-vhdl
:before
{
content
:
'VHDL'
;
}
pre
.src-xml
:before
{
content
:
'XML'
;
}
pre
.src-nxml
:before
{
content
:
'XML'
;
}
/* add a generic configuration mode; LaTeX export needs an additional
(add-to-list 'org-latex-listings-langs '(conf " ")) in .emacs */
pre
.src-conf
:before
{
content
:
'Configuration File'
;
}
table
{
border-collapse
:
collapse
;
}
caption
.t-above
{
caption-side
:
top
;
}
caption
.t-bottom
{
caption-side
:
bottom
;
}
td
,
th
{
vertical-align
:
top
;
}
th
.org-right
{
text-align
:
center
;
}
th
.org-left
{
text-align
:
center
;
}
th
.org-center
{
text-align
:
center
;
}
td
.org-right
{
text-align
:
right
;
}
td
.org-left
{
text-align
:
left
;
}
td
.org-center
{
text-align
:
center
;
}
dt
{
font-weight
:
bold
;
}
.footpara
{
display
:
inline
;
}
.footdef
{
margin-bottom
:
1em
;
}
.figure
{
padding
:
1em
;
}
.figure
p
{
text-align
:
center
;
}
.inlinetask
{
padding
:
10px
;
border
:
2px
solid
gray
;
margin
:
10px
;
background
:
#ffffcc
;
}
#org-div-home-and-up
{
text-align
:
right
;
font-size
:
70%
;
white-space
:
nowrap
;
}
textarea
{
overflow-x
:
auto
;
}
.linenr
{
font-size
:
smaller
}
.code-highlighted
{
background-color
:
#ffff00
;
}
.org-info-js_info-navigation
{
border-style
:
none
;
}
#org-info-js_console-label
{
font-size
:
10px
;
font-weight
:
bold
;
white-space
:
nowrap
;
}
.org-info-js_search-highlight
{
background-color
:
#ffff00
;
color
:
#000000
;
font-weight
:
bold
;
}
.org-svg
{
width
:
90%
;
}
/*]]>*/
--
>
</style>
<link
rel=
"stylesheet"
type=
"text/css"
href=
"http://www.pirilampo.org/styles/readtheorg/css/htmlize.css"
/>
<link
rel=
"stylesheet"
type=
"text/css"
href=
"http://www.pirilampo.org/styles/readtheorg/css/readtheorg.css"
/>
<script
src=
"https://ajax.googleapis.com/ajax/libs/jquery/2.1.3/jquery.min.js"
></script>
<script
src=
"https://maxcdn.bootstrapcdn.com/bootstrap/3.3.4/js/bootstrap.min.js"
></script>
<script
type=
"text/javascript"
src=
"http://www.pirilampo.org/styles/lib/js/jquery.stickytableheaders.js"
></script>
<script
type=
"text/javascript"
src=
"http://www.pirilampo.org/styles/readtheorg/js/readtheorg.js"
></script>
<script
type=
"text/javascript"
>
/*
@licstart The following is the entire license notice for the
JavaScript code in this tag.
Copyright (C) 2012-2019 Free Software Foundation, Inc.
The JavaScript code in this tag is free software: you can
redistribute it and/or modify it under the terms of the GNU
General Public License (GNU GPL) as published by the Free Software
Foundation, either version 3 of the License, or (at your option)
any later version. The code is distributed WITHOUT ANY WARRANTY;
without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the GNU GPL for more details.
As additional permission under GNU GPL version 3 section 7, you
may distribute non-source (e.g., minimized or compacted) forms of
that code without the copy of the GNU GPL normally required by
section 4, provided you include this license notice and a URL
through which recipients can access the Corresponding Source.
@licend The above is the entire license notice
for the JavaScript code in this tag.
*/
<!--
/*--><![CDATA[/*><!--*/
function
CodeHighlightOn
(
elem
,
id
)
{
var
target
=
document
.
getElementById
(
id
);
if
(
null
!=
target
)
{
elem
.
cacheClassElem
=
elem
.
className
;
elem
.
cacheClassTarget
=
target
.
className
;
target
.
className
=
"code-highlighted"
;
elem
.
className
=
"code-highlighted"
;
}
}
function
CodeHighlightOff
(
elem
,
id
)
{
var
target
=
document
.
getElementById
(
id
);
if
(
elem
.
cacheClassElem
)
elem
.
className
=
elem
.
cacheClassElem
;
if
(
elem
.
cacheClassTarget
)
target
.
className
=
elem
.
cacheClassTarget
;
}
/*]]>*/
//-->
</script>
<script
type=
"text/x-mathjax-config"
>
MathJax
.
Hub
.
Config
({
displayAlign
:
"center"
,
displayIndent
:
"0em"
,
"HTML-CSS"
:
{
scale
:
100
,
linebreaks
:
{
automatic
:
"false"
},
webFont
:
"TeX"
},
SVG
:
{
scale
:
100
,
linebreaks
:
{
automatic
:
"false"
},
font
:
"TeX"
},
NativeMML
:
{
scale
:
100
},
TeX
:
{
equationNumbers
:
{
autoNumber
:
"AMS"
},
MultLineWidth
:
"85%"
,
TagSide
:
"right"
,
TagIndent
:
".8em"
}
});
</script>
<script
type=
"text/javascript"
src=
"https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.0/MathJax.js?config=TeX-AMS_HTML"
></script>
</head>
<body>
<div
id=
"content"
>
<h1
class=
"title"
>
À propos du calcul de \(\pi\)
</h1>
<div
id=
"table-of-contents"
>
<h2>
Table des matières
</h2>
<div
id=
"text-table-of-contents"
>
<ul>
<li><a
href=
"#org9e00bf0"
>
1. En demandant à la lib maths
</a></li>
<li><a
href=
"#org7830e5e"
>
2. En utilisant la méthodes des aiguilles de Buffon
</a></li>
<li><a
href=
"#orgec609c5"
>
3. Avec un argument "fréquentiel" de surface
</a></li>
</ul>
</div>
</div>
<div
id=
"outline-container-org9e00bf0"
class=
"outline-2"
>
<h2
id=
"org9e00bf0"
><span
class=
"section-number-2"
>
1
</span>
En demandant à la lib maths
</h2>
<div
class=
"outline-text-2"
id=
"text-1"
>
<p>
Mon ordinateur m'indique que \(\pi\) vaut approximativement :
</p>
<div
class=
"org-src-container"
>
<pre
class=
"src src-python"
><span
style=
"color: #4d9391;"
>
from
</span>
math
<span
style=
"color: #4d9391;"
>
import
</span>
*
pi
</pre>
</div>
<pre
class=
"example"
>
3.141592653589793
</pre>
</div>
</div>
<div
id=
"outline-container-org7830e5e"
class=
"outline-2"
>
<h2
id=
"org7830e5e"
><span
class=
"section-number-2"
>
2
</span>
En utilisant la méthodes des aiguilles de Buffon
</h2>
<div
class=
"outline-text-2"
id=
"text-2"
>
<p>
Mais calculé avec la
<b>
méthode
</b>
des
<a
href=
"https://fr.wikipedia.org/wiki/Aiguille_de_Buffon"
>
aiguilles de Buffon
</a>
, on obtiendrait comme
<b>
approximation
</b>
:
</p>
<div
class=
"org-src-container"
>
<pre
class=
"src src-python"
><span
style=
"color: #4d9391;"
>
import
</span>
numpy
<span
style=
"color: #4d9391;"
>
as
</span>
np
np.random.seed
<span
style=
"color: #6bd9db;"
>
(
</span>
seed=
<span
style=
"color: #cd5c60;"
>
42
</span><span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #807f96;"
>
N
</span>
=
<span
style=
"color: #cd5c60;"
>
10000
</span>
<span
style=
"color: #807f96;"
>
x
</span>
= np.random.uniform
<span
style=
"color: #6bd9db;"
>
(
</span>
size=N, low=
<span
style=
"color: #cd5c60;"
>
0
</span>
, high=
<span
style=
"color: #cd5c60;"
>
1
</span><span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #807f96;"
>
theta
</span>
= np.random.uniform
<span
style=
"color: #6bd9db;"
>
(
</span>
size=N, low=
<span
style=
"color: #cd5c60;"
>
0
</span>
, high=pi/
<span
style=
"color: #cd5c60;"
>
2
</span><span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #cd5c60;"
>
2
</span>
/
<span
style=
"color: #6bd9db;"
>
(
</span><span
style=
"color: #80bcb6;"
>
sum
</span><span
style=
"color: #ab98b5;"
>
(
</span><span
style=
"color: #5D8272;"
>
(
</span>
x + np.sin
<span
style=
"color: #807f96;"
>
(
</span>
theta
<span
style=
"color: #807f96;"
>
)
</span><span
style=
"color: #5D8272;"
>
)
</span>
>
<span
style=
"color: #cd5c60;"
>
1
</span><span
style=
"color: #ab98b5;"
>
)
</span>
/ N
<span
style=
"color: #6bd9db;"
>
)
</span>
</pre>
</div>
<pre
class=
"example"
>
3.128911138923655
</pre>
</div>
</div>
<div
id=
"outline-container-orgec609c5"
class=
"outline-2"
>
<h2
id=
"orgec609c5"
><span
class=
"section-number-2"
>
3
</span>
Avec un argument "fréquentiel" de surface
</h2>
<div
class=
"outline-text-2"
id=
"text-3"
>
<p>
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 Y(0, 1)\) alors \(P[X^2 + Y^2 \le 1] = \pi / 4\) (voir méthode de
<a
href=
"https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Monte-Carlo#D%C3%A9termination_de_la_valeur_de_%CF%80"
>
Monte
Carlo sur Wikipedia
</a>
). Le code suivant illustre ce fait :
</p>
<div
class=
"org-src-container"
>
<pre
class=
"src src-python"
><span
style=
"color: #4d9391;"
>
import
</span>
matplotlib.pyplot
<span
style=
"color: #4d9391;"
>
as
</span>
plt
np.random.seed
<span
style=
"color: #6bd9db;"
>
(
</span>
seed=
<span
style=
"color: #cd5c60;"
>
42
</span><span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #807f96;"
>
N
</span>
=
<span
style=
"color: #cd5c60;"
>
1000
</span>
<span
style=
"color: #807f96;"
>
x
</span>
= np.random.uniform
<span
style=
"color: #6bd9db;"
>
(
</span>
size=N, low=
<span
style=
"color: #cd5c60;"
>
0
</span>
, high=
<span
style=
"color: #cd5c60;"
>
1
</span><span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #807f96;"
>
y
</span>
= np.random.uniform
<span
style=
"color: #6bd9db;"
>
(
</span>
size=N, low=
<span
style=
"color: #cd5c60;"
>
0
</span>
, high=
<span
style=
"color: #cd5c60;"
>
1
</span><span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #807f96;"
>
accept
</span>
=
<span
style=
"color: #6bd9db;"
>
(
</span>
x*x+y*y
<span
style=
"color: #6bd9db;"
>
)
</span>
<
=
<span
style=
"color: #cd5c60;"
>
1
</span>
<span
style=
"color: #807f96;"
>
reject
</span>
= np.logical_not
<span
style=
"color: #6bd9db;"
>
(
</span>
accept
<span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #807f96;"
>
fig
</span>
,
<span
style=
"color: #807f96;"
>
ax
</span>
= plt.subplots
<span
style=
"color: #6bd9db;"
>
(
</span><span
style=
"color: #cd5c60;"
>
1
</span><span
style=
"color: #6bd9db;"
>
)
</span>
ax.scatter
<span
style=
"color: #6bd9db;"
>
(
</span>
x
<span
style=
"color: #ab98b5;"
>
[
</span>
accept
<span
style=
"color: #ab98b5;"
>
]
</span>
, y
<span
style=
"color: #ab98b5;"
>
[
</span>
accept
<span
style=
"color: #ab98b5;"
>
]
</span>
, c=
<span
style=
"color: #6fb593;"
>
'b'
</span>
, alpha=
<span
style=
"color: #cd5c60;"
>
0
</span>
.
<span
style=
"color: #cd5c60;"
>
2
</span>
, edgecolor=
<span
style=
"color: #ab98b5;"
>
None
</span><span
style=
"color: #6bd9db;"
>
)
</span>
ax.scatter
<span
style=
"color: #6bd9db;"
>
(
</span>
x
<span
style=
"color: #ab98b5;"
>
[
</span>
reject
<span
style=
"color: #ab98b5;"
>
]
</span>
, y
<span
style=
"color: #ab98b5;"
>
[
</span>
reject
<span
style=
"color: #ab98b5;"
>
]
</span>
, c=
<span
style=
"color: #6fb593;"
>
'r'
</span>
, alpha=
<span
style=
"color: #cd5c60;"
>
0
</span>
.
<span
style=
"color: #cd5c60;"
>
2
</span>
, edgecolor=
<span
style=
"color: #ab98b5;"
>
None
</span><span
style=
"color: #6bd9db;"
>
)
</span>
ax.set_aspect
<span
style=
"color: #6bd9db;"
>
(
</span><span
style=
"color: #6fb593;"
>
'equal'
</span><span
style=
"color: #6bd9db;"
>
)
</span>
plt.savefig
<span
style=
"color: #6bd9db;"
>
(
</span>
matplot_lib_filename
<span
style=
"color: #6bd9db;"
>
)
</span>
<span
style=
"color: #4d9391;"
>
print
</span><span
style=
"color: #6bd9db;"
>
(
</span>
matplot_lib_filename
<span
style=
"color: #6bd9db;"
>
)
</span>
</pre>
</div>
<div
class=
"figure"
>
<p><img
src=
"pi_monte_carlo.png"
alt=
"pi_monte_carlo.png"
/>
</p>
</div>
<p>
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 :
</p>
<div
class=
"org-src-container"
>
<pre
class=
"src src-python"
><span
style=
"color: #cd5c60;"
>
4
</span>
* np.mean
<span
style=
"color: #6bd9db;"
>
(
</span>
accept
<span
style=
"color: #6bd9db;"
>
)
</span>
</pre>
</div>
<pre
class=
"example"
>
3.112
</pre>
</div>
</div>
</div>
<div
id=
"postamble"
class=
"status"
>
<p
class=
"date"
>
Date: 2020-07-29 mer. 00:00
</p>
<p
class=
"author"
>
Auteur: Pierre AYOUB
</p>
<p
class=
"date"
>
Created: 2020-07-29 mer. 14:28
</p>
<p
class=
"validation"
><a
href=
"http://validator.w3.org/check?uri=referer"
>
Validate
</a></p>
</div>
</body>
</html>
module2/exo1/toy_document_orgmode_python_fr.org
View file @
2c85c94a
#+TITLE:
Votre titre
#+TITLE:
À propos du calcul de $\pi$
#+AUTHOR:
Votre nom
#+AUTHOR:
Pierre AYOUB
#+DATE:
La date du jour
#+DATE:
[2020-07-29 mer.]
#+LANGUAGE: fr
#+LANGUAGE: fr
# #+PROPERTY: header-args :eval never-export
# #+PROPERTY: header-args :eval never-export
...
@@ -11,83 +11,72 @@
...
@@ -11,83 +11,72 @@
#+HTML_HEAD: <script type="text/javascript" src="http://www.pirilampo.org/styles/lib/js/jquery.stickytableheaders.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="http://www.pirilampo.org/styles/lib/js/jquery.stickytableheaders.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="http://www.pirilampo.org/styles/readtheorg/js/readtheorg.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="http://www.pirilampo.org/styles/readtheorg/js/readtheorg.js"></script>
* Quelques explications
* En demandant à la lib maths
Ceci est un document org-mode avec quelques exemples de code
Mon ordinateur m'indique que $\pi$ vaut approximativement :
python. Une fois ouvert dans emacs, ce document peut aisément être
exporté au format HTML, PDF, et Office. Pour plus de détails sur
#+BEGIN_SRC python :results value :exports both :session *sess_python*
org-mode vous pouvez consulter https://orgmode.org/guide/.
from math import *
pi
Lorsque vous utiliserez le raccourci =C-c C-e h o=, ce document sera
#+END_SRC
compilé en html. Tout le code contenu sera ré-exécuté, les résultats
récupérés et inclus dans un document final. Si vous ne souhaitez pas
#+RESULTS:
ré-exécuter tout le code à chaque fois, il vous suffit de supprimer
: 3.141592653589793
le # et l'espace qui sont devant le ~#+PROPERTY:~ au début de ce
document.
* En utilisant la méthodes des aiguilles de Buffon
Comme nous vous l'avons montré dans la vidéo, on inclue du code
Mais calculé avec la *méthode* des [[https://fr.wikipedia.org/wiki/Aiguille_de_Buffon][aiguilles de Buffon]], on obtiendrait comme
python de la façon suivante (et on l'exécute en faisant ~C-c C-c~):
*approximation* :
#+begin_src python :results output :exports both
#+BEGIN_SRC python :results value :exports both :session *sess_python*
print("Hello world!")
import numpy as np
#+end_src
np.random.seed(seed=42)
N = 10000
#+RESULTS:
x = np.random.uniform(size=N, low=0, high=1)
: Hello world!
theta = np.random.uniform(size=N, low=0, high=pi/2)
2 / (sum((x + np.sin(theta)) > 1) / N)
Voici la même chose, mais avec une session python, donc une
#+END_SRC
persistance d'un bloc à l'autre (et on l'exécute toujours en faisant
~C-c C-c~).
#+RESULTS:
#+begin_src python :results output :session :exports both
: 3.128911138923655
import numpy
x=numpy.linspace(-15,15)
* Avec un argument "fréquentiel" de surface
print(x)
#+end_src
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
#+RESULTS:
\sim Y(0, 1)$ alors $P[X^2 + Y^2 \le 1] = \pi / 4$ (voir méthode de [[https://fr.wikipedia.org/wiki/M%25C3%25A9thode_de_Monte-Carlo#D%25C3%25A9termination_de_la_valeur_de_%25CF%2580][Monte
#+begin_example
Carlo sur Wikipedia]]). Le code suivant illustre ce fait :
[-15. -14.3877551 -13.7755102 -13.16326531 -12.55102041
-11.93877551 -11.32653061 -10.71428571 -10.10204082 -9.48979592
#+HEADER: :var matplot_lib_filename="pi_monte_carlo.png"
-8.87755102 -8.26530612 -7.65306122 -7.04081633 -6.42857143
#+BEGIN_SRC python :results output file :exports both :session *sess_python*
-5.81632653 -5.20408163 -4.59183673 -3.97959184 -3.36734694
import matplotlib.pyplot as plt
-2.75510204 -2.14285714 -1.53061224 -0.91836735 -0.30612245
0.30612245 0.91836735 1.53061224 2.14285714 2.75510204
np.random.seed(seed=42)
3.36734694 3.97959184 4.59183673 5.20408163 5.81632653
N = 1000
6.42857143 7.04081633 7.65306122 8.26530612 8.87755102
x = np.random.uniform(size=N, low=0, high=1)
9.48979592 10.10204082 10.71428571 11.32653061 11.93877551
y = np.random.uniform(size=N, low=0, high=1)
12.55102041 13.16326531 13.7755102 14.3877551 15. ]
#+end_example
accept = (x*x+y*y) <= 1
reject = np.logical_not(accept)
Et enfin, voici un exemple de sortie graphique:
#+begin_src python :results output file :session :var matplot_lib_filename="./cosxsx.png" :exports results
fig, ax = plt.subplots(1)
import matplotlib.pyplot as plt
ax.scatter(x[accept], y[accept], c='b', alpha=0.2, edgecolor=None)
ax.scatter(x[reject], y[reject], c='r', alpha=0.2, edgecolor=None)
plt.figure(figsize=(10,5))
ax.set_aspect('equal')
plt.plot(x,numpy.cos(x)/x)
plt.tight_layout()
plt.savefig(matplot_lib_filename)
print(matplot_lib_filename)
plt.savefig(matplot_lib_filename)
#+END_SRC
print(matplot_lib_filename)
#+end_src
#+RESULTS:
[[file:pi_monte_carlo.png]]
#+RESULTS:
[[file:./cosxsx.png]]
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 :
Vous remarquerez le paramètre ~:exports results~ qui indique que le code
#+BEGIN_SRC python :results value :exports both :session *sess_python*
ne doit pas apparaître dans la version finale du document. Nous vous
4 * np.mean(accept)
recommandons dans le cadre de ce MOOC de ne pas changer ce paramètre
#+END_SRC
(indiquer ~both~) car l'objectif est que vos analyses de données soient
parfaitement transparentes pour être reproductibles.
#+RESULTS:
: 3.112
Attention, la figure ainsi générée n'est pas stockée dans le document
org. C'est un fichier ordinaire, ici nommé ~cosxsx.png~. N'oubliez pas
de le committer si vous voulez que votre analyse soit lisible et
compréhensible sur GitLab.
Enfin, n'oubliez pas que nous vous fournissons dans les ressources de
ce MOOC une configuration avec un certain nombre de raccourcis
claviers permettant de créer rapidement les blocs de code python (en
faisant ~<p~, ~<P~ ou ~<PP~ suivi de ~Tab~).
Maintenant, à vous de jouer! Vous pouvez effacer toutes ces
informations et les remplacer par votre document computationnel.
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment