From 1a0fd2e8cb900604c6c347a8bff8a813755e9c03 Mon Sep 17 00:00:00 2001 From: 934efbaf7a6731c72b5e03e55b33b979 <934efbaf7a6731c72b5e03e55b33b979@app-learninglab.inria.fr> Date: Wed, 2 Sep 2020 09:23:53 +0000 Subject: [PATCH] First part of exercice --- module2/exo1/toy_notebook_fr.ipynb | 98 +++++++++++++++++++++++++++++- 1 file changed, 95 insertions(+), 3 deletions(-) diff --git a/module2/exo1/toy_notebook_fr.ipynb b/module2/exo1/toy_notebook_fr.ipynb index 0bbbe37..adb3213 100644 --- a/module2/exo1/toy_notebook_fr.ipynb +++ b/module2/exo1/toy_notebook_fr.ipynb @@ -1,5 +1,98 @@ { - "cells": [], + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# toy_notebook_fr" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## March 28, 2019" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "1. A propos du calcul de $\\pi$" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + " 1. En demendant à la lib Maths" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Mon ordinateur indique que $\\pi$ vaut *approximativement*:" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "3.141592653589793\n" + ] + } + ], + "source": [ + "from math import *\n", + "print(pi)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + " 2. En utilisant la méthode des [aiguilles de Buffon](https://fr.wikipedia.org/wiki/Aiguille_de_Buffon), on obtiendrait comme **approximation**:" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "3.128911138923655" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "import numpy as np\n", + "np.random.seed(seed=42)\n", + "N=10000\n", + "x=np.random.uniform(size=N, low=0, high=1)\n", + "theta=np.random.uniform(size=N, low=0, high=pi/2)\n", + "2/(sum((x+np.sin(theta))>1)/N)" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], "metadata": { "kernelspec": { "display_name": "Python 3", @@ -16,10 +109,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.3" + "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 } - -- 2.18.1