{ "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$\n", "## 1.1 En demandant à la lib maths\n", "Mon ordinateur m'indique que $\\pi$ vaut _approximativement_" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "3.141592653589793\n" ] } ], "source": [ "from math import *\n", "print(pi)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1.2 En utilisant la méthode des aiguilles de Buffon\n", "Mais calculé avec la __méthode__ des aiguilles 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": [ "## 1.3 Avec un argument \"fréquentiel\" de surface\n", "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 " ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }