{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "e2dbdbb0",
   "metadata": {},
   "source": [
    "# Incertitude type composée pour la mesure d'une distance le long du banc d'optique."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7111e6ac",
   "metadata": {},
   "source": [
    "Pour exploiter la méthode d'autocollimation, on doit estimer la distance séparant l'objet de la lentille, le long d'un banc d'optique gradué.\n",
    "\n",
    "Si on note $X_{objet}$ la position de l'objet et $X_{L}$ la position de la lentille le long du banc d'optique, la distance focale $f'$ de la lentille est évaluée avec la relation : $f'=X_{L}-X_{objet}$\n",
    "\n",
    "Pour obtenir une évaluation de $f'$, on doit donc déterminer la valeur mesurée $f'_{mes}$ de $f'$ ainsi que l'incertitude type $u(f')$ à partir des valeurs mesurées directement pour $X_{objet}$ la position de l'objet obtenue avec l'incertitude type $u\\left(X_{objet}\\right)$ et $X_{L}$ la position de la lentille connue avec l'incertitude type $u\\left(X_{L}\\right)$.\n",
    "\n",
    "Dans cette situation où on évalue une grandeur à partir de la combinaison de plusieurs mesures, on dit qu'il faut réaliser la compisition des incertitudes pour obtenir l'incertitude type sur la grandeur. Pour la formule de type somme-différence qui fait le lien entre $f'$, $X_{objet}$ et $X_{L}$, la formule théorique de composition des incertitudes est la suivante :\n",
    "\n",
    "estimation de la différence :$f'=X_{objet}-X_{L}$ et estimation de l'incertitude : $u^2 (f')=u^2\\left(X_{objet}\\right)+u^2\\left(X_{L}\\right)$\n",
    "\n",
    "On se propose de la vérifier par simulation de Monte Carlo.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bad803d7",
   "metadata": {},
   "source": [
    "## Evaluation de la distance focale de la lentille convergente par la méthode d'autocollimation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "281c0e61",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAPKElEQVR4nO3db8yddX3H8fdnFPEPLGC4YYWW3WiqWVm2Yu50biwLG4syNFYfYEo2VzOW+gAz3UwmaDL1QRPMFLcH06WKs9tQbBQncW6KTGN8INgiYktlq9LBbTtap5u4JbiW7x6cq/FYzv33nHP/+d3vV3Jyn/M71++c75fTfM51fue6DqkqJElt+ZnlLkCSNHqGuyQ1yHCXpAYZ7pLUIMNdkhq0brkLALjwwgtrcnJyucuQpFVl//7936uqiUH3rYhwn5ycZN++fctdhiStKkn+fab7XJaRpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGrYgzVCVp1CZv/scZ7zty6yuWsJLl4Z67JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNmjPck2xM8sUkh5IcTPKmbvydSb6b5MHucl3fnFuSHE7ySJKXj7MBSdIzzeckppPAW6rqgSTnAfuT3NPd976qek//xkk2A9uBK4BLgC8keVFVnRpl4f1mOllhLZyoIEmDzLnnXlXHquqB7vqTwCHg0lmmbAPurKqnqupR4DCwdRTFSpLmZ0E/P5BkErgSuA+4Cnhjkt8H9tHbu/8BveD/at+0aQa8GSTZCewEuOyyyxZTuyStGCvt5w7m/YVqknOBTwJvrqofAh8AXghsAY4B7z296YDp9YyBqt1VNVVVUxMTEwutW5I0i3ntuSc5m16w31FVdwFU1RN9938Q+Ex3cxrY2Dd9A3B0JNUu0GLfSVfaO7AkLdSc4Z4kwO3Aoaq6rW98fVUd626+BjjQXb8b+GiS2+h9oboJuH+kVWvV8I1SWh7z2XO/Cngd8M0kD3ZjbwNuSLKF3pLLEeANAFV1MMle4GF6R9rcNM4jZbT8ZgtwSctjznCvqq8weB39s7PM2QXsGqIuSSPip6fltxyvgWeoSlKD/D8xSdI8raYlSMN9hPz4u7L5+oyO/y1XPpdlJKlBhrskNchlGUmr1jjWwFfTuvps3HOXpAYZ7pLUIMNdkhrkmvsqtti1QQ9VW14eRqilYLhLGinfvFYGw10agkGmlcpwl+awUg6Na+GNxKXEpeMXqpLUIPfcpQaslE8XWjkMdwnDUe0x3DUvS32adwtrrKv9DWO11z+blns7zTV3SWqQe+76KWthj0ZaCwz3Fc6wlbQYhrtWpNbX46VxM9ylFcRPahoVv1CVpAat2T33tbyHtJZ7l9YK99wlqUGGuyQ1aM0uyyw1j/6QtJQM9wVyvVrSauCyjCQ1yHCXpAbNuSyTZCPwt8DPAU8Du6vqL5M8H/g4MAkcAV5bVT/o5twC3AicAv6oqj43luq1Jvn9hTS3+ay5nwTeUlUPJDkP2J/kHuD1wL1VdWuSm4Gbgbcm2QxsB64ALgG+kORFVXVqPC1IP+F3IlLPnMsyVXWsqh7orj8JHAIuBbYBe7rN9gCv7q5vA+6sqqeq6lHgMLB1xHVLkmaxoDX3JJPAlcB9wMVVdQx6bwDARd1mlwKP902b7sYkSUtk3odCJjkX+CTw5qr6YZIZNx0wVgMebyewE+Cyyy6bbxmS1iCX2xZuXnvuSc6mF+x3VNVd3fATSdZ3968Hjnfj08DGvukbgKNnPmZV7a6qqaqampiYWGz9kqQB5gz39HbRbwcOVdVtfXfdDezoru8APt03vj3JOUkuBzYB94+uZEnSXOazLHMV8Drgm0ke7MbeBtwK7E1yI/AYcD1AVR1Mshd4mN6RNjd5pIwkLa05w72qvsLgdXSAa2aYswvYNURdkqQheIaqJDXIHw6T1jCPQmmXe+6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGuTPD6wAngIuadTcc5ekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGzRnuST6c5HiSA31j70zy3SQPdpfr+u67JcnhJI8kefm4CpckzWw+e+4fAa4dMP6+qtrSXT4LkGQzsB24opvz/iRnjapYSdL8zBnuVfVl4PvzfLxtwJ1V9VRVPQocBrYOUZ8kaRGGWXN/Y5KHumWbC7qxS4HH+7aZ7sYkSUtoseH+AeCFwBbgGPDebjwDtq1BD5BkZ5J9SfadOHFikWVIkgZZVLhX1RNVdaqqngY+yE+WXqaBjX2bbgCOzvAYu6tqqqqmJiYmFlOGJGkGiwr3JOv7br4GOH0kzd3A9iTnJLkc2ATcP1yJkqSFWjfXBkk+BlwNXJhkGngHcHWSLfSWXI4AbwCoqoNJ9gIPAyeBm6rq1FgqlyTNaM5wr6obBgzfPsv2u4BdwxQlSRqOZ6hKUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBs0Z7kk+nOR4kgN9Y89Pck+Sf+v+XtB33y1JDid5JMnLx1W4JGlm89lz/whw7RljNwP3VtUm4N7uNkk2A9uBK7o5709y1siqlSTNy5zhXlVfBr5/xvA2YE93fQ/w6r7xO6vqqap6FDgMbB1NqZKk+VrsmvvFVXUMoPt7UTd+KfB433bT3ZgkaQmN+gvVDBirgRsmO5PsS7LvxIkTIy5Dkta2xYb7E0nWA3R/j3fj08DGvu02AEcHPUBV7a6qqaqampiYWGQZkqRBFhvudwM7uus7gE/3jW9Pck6Sy4FNwP3DlShJWqh1c22Q5GPA1cCFSaaBdwC3AnuT3Ag8BlwPUFUHk+wFHgZOAjdV1akx1S5JmsGc4V5VN8xw1zUzbL8L2DVMUZKk4XiGqiQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDVo3zOQkR4AngVPAyaqaSvJ84OPAJHAEeG1V/WC4MiVJCzGKPfffrKotVTXV3b4ZuLeqNgH3drclSUtoHMsy24A93fU9wKvH8BySpFkMG+4FfD7J/iQ7u7GLq+oYQPf3okETk+xMsi/JvhMnTgxZhiSp31Br7sBVVXU0yUXAPUm+Nd+JVbUb2A0wNTVVQ9YhSeoz1J57VR3t/h4HPgVsBZ5Ish6g+3t82CIlSQuz6HBP8rwk552+DrwMOADcDezoNtsBfHrYIiVJCzPMsszFwKeSnH6cj1bVPyf5GrA3yY3AY8D1w5cpSVqIRYd7VX0H+OUB4/8JXDNMUZKk4XiGqiQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWrQ2MI9ybVJHklyOMnN43oeSdIzjSXck5wF/BXwO8Bm4IYkm8fxXJKkZxrXnvtW4HBVfaeqfgzcCWwb03NJks6wbkyPeynweN/taeBX+jdIshPY2d38UZJHxlTLOF0IfG+5i1hi9rw2rLWel63fvHuo6T8/0x3jCvcMGKufulG1G9g9pudfEkn2VdXUctexlOx5bVhrPbfY77iWZaaBjX23NwBHx/RckqQzjCvcvwZsSnJ5kmcB24G7x/RckqQzjGVZpqpOJnkj8DngLODDVXVwHM+1zFb1stIi2fPasNZ6bq7fVNXcW0mSVhXPUJWkBhnuktQgw32AJM9Ocn+SbyQ5mORd3fj13e2nk8x42FSS85N8Ism3khxK8qtLV/3CjaDfP+62O5DkY0mevXTVL84sPf9597o9lORTSc6fYf6q+3mNYXpOsjHJF7t/zweTvGnJG1iEYV/nbtuzknw9yWeWrPBRqCovZ1zoHad/bnf9bOA+4KXALwAvBr4ETM0yfw/wh931ZwHnL3dP4+qX3glrjwLP6W7vBV6/3D0N0fPLgHXd+LuBdw+YexbwbeAF3ev7DWDzcvc05p7XAy/prp8H/GvrPfc9xp8AHwU+s9z9LOTinvsA1fOj7ubZ3aWq6lBVzXombZKfBX4DuL17rB9X1X+Ns95hDdNvZx3wnCTrgOeyCs5pmKXnz1fVyW78q/TO0TjTqvx5jWF6rqpjVfVAd/1J4BC9N/YVbcjXmSQbgFcAHxp7sSNmuM+g+yj2IHAcuKeq7pvn1BcAJ4C/6T7KfSjJ88ZV56gstt+q+i7wHuAx4Bjw31X1+bEVOkLz6PkPgH8aMHXQz2us+KCDoXruf4xJ4Ep6e8Er3pA9/wXwp8DTYytwTAz3GVTVqaraQu8dfWuSX5zn1HXAS4APVNWVwP8AK35NdrH9JrmA3l7r5cAlwPOS/N7YCh2h2XpO8nbgJHDHgKlz/rzGSjVEz6e3ORf4JPDmqvrhmMsdicX2nOSVwPGq2r9UtY6S4T6HbknlS8C185wyDUz37R18gl7YrwqL6Pe3gUer6kRV/R9wF/Br46luPM7sOckO4JXA71a36HqGVf/zGovomSRn0wv2O6rqrqWpdHQW0fNVwKuSHKG39PZbSf5+SYodAcN9gCQTp789T/IcegH2rfnMrar/AB5P8uJu6Brg4XHUOSrD9EtvOealSZ6bJPT6PTSWQkdopp6TXAu8FXhVVf3vDNNX5c9rDNNz99reDhyqqtuWqOShDdNzVd1SVRuqapLea/wvVbUqPpUCHi0z6AL8EvB14CHgAPBn3fhr6O21PQU8AXyuG78E+Gzf/C3Avm7+PwAXLHdPY+73XfTeDA4Afwecs9w9DdHzYXrr6Q92l7+eoefr6B0x8m3g7cvdz7h7Bn6d3tLTQ33bXbfcPY37de57nKtZZUfL+PMDktQgl2UkqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWrQ/wPcep1c6n/aPQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAPrElEQVR4nO3cb4xcV33G8e/TBGgpSCTEiYzj1BEyEgkCI61cqlRV2khNGpBCKqVyXtBURDVSEylpqYoNUqGq3IZS/rwpVKZEpC2QWgqIiKaAiaABCQh2MMGOceOSlBi7sSlIBCGltfPri7kWgz27szuzszN7/P1Io7n33Hvn/nw8+8yZM3cmVYUkqS2/MO0CJEnLz3CXpAYZ7pLUIMNdkhpkuEtSg86fdgEAF110UW3YsGHaZUjSqrJ3794fVNWaQdtmItw3bNjAnj17pl2GJK0qSf5rvm1Oy0hSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoNm4huq0mq1Ydu/zrvtybtev4KVSD/PcNeiGGLS6uK0jCQ1yHCXpAY5LSPhtJPa48hdkhpkuEtSgwx3SWrQ0HBPsj7JF5McTHIgyR1d+7uSfD/Jvu52fd8x25McTnIoybWT/AdIks62mA9UTwJvrapHkrwY2Jtkd7ft/VX1t/07J7kC2AJcCbwM+EKSV1TVqeUsXJI0v6Ej96o6VlWPdMvPAAeBdQsccgNwb1U9W1VPAIeBzctRrCRpcZZ0KWSSDcBrga8DVwG3J/l9YA+90f2P6AX/1/oOO8KAF4MkW4GtAJdddtkotUvSilpNl8wuOtyTvAi4D7izqn6c5EPAXwLV3b8XeDOQAYfXWQ1VO4GdAHNzc2dtl7R4qyl0tDIWdbVMkufRC/aPVdUnAarq6ao6VVXPAR/mZ1MvR4D1fYdfChxdvpIlScMs5mqZAB8BDlbV+/ra1/btdiOwv1u+H9iS5AVJLgc2Ag8vX8mSpGEWMy1zFfAm4NtJ9nVtbwduTrKJ3pTLk8BbAKrqQJJdwGP0rrS5zStlJGllDQ33qvoKg+fRH1jgmB3AjjHqkiSNoYkfDpvvwyQ/SJJ0rvLnBySpQYa7JDWoiWkZScvPa+dXN0fuktQgR+4rxFHQyrCfNS2z9txz5C5JDXLkvkSz9uoszRr/RmaDI3dJapAjd/2chUZdoxwziZHaKDVOwyhfrlst/zbNPsNdWmEGuFaC4d4o5z2lc9s5G+4thJ8jQGl+LfyNj+OcDXdJq2eA4I8DLp3hLg2xWgJwPqu9fo3GSyElqUGO3CUtme8GZp8jd0lqkOEuSQ1yWuYc5FtqnetW+m9gGpdlGu4z4Fy/HncQX4Ck8RjumihfuKTpMNwlzQTfrS0vw30Z+eSUFubfyMox3CWtWr5YzM9LISWpQY7cZ5wjE0mjcOQuSQ0y3CWpQU7L6JzhFJfOJUNH7knWJ/likoNJDiS5o2u/MMnuJI939xf0HbM9yeEkh5JcO8l/gCTpbIsZuZ8E3lpVjyR5MbA3yW7gD4AHq+quJNuAbcDbklwBbAGuBF4GfCHJK6rq1GT+CcvPEZ6k1W5ouFfVMeBYt/xMkoPAOuAG4Oput3uALwFv69rvrapngSeSHAY2A19d7uKHMaQlnauW9IFqkg3Aa4GvA5d0wX/6BeDibrd1wFN9hx3p2iRJK2TR4Z7kRcB9wJ1V9eOFdh3QVgMeb2uSPUn2nDhxYrFlSJIWYVHhnuR59IL9Y1X1ya756SRru+1rgeNd+xFgfd/hlwJHz3zMqtpZVXNVNbdmzZpR65ckDTB0zj1JgI8AB6vqfX2b7gduAe7q7j/d1/7xJO+j94HqRuDh5SxabfAzEWlyFnO1zFXAm4BvJ9nXtb2dXqjvSnIr8D3gJoCqOpBkF/AYvSttbltNV8pIUgsWc7XMVxg8jw5wzTzH7AB2jFGXJGkM/vyAJDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJatDQcE9yd5LjSfb3tb0ryfeT7Otu1/dt257kcJJDSa6dVOGSpPktZuT+UeC6Ae3vr6pN3e0BgCRXAFuAK7tjPpjkvOUqVpK0OEPDvaoeAn64yMe7Abi3qp6tqieAw8DmMeqTJI1gnDn325M82k3bXNC1rQOe6tvnSNcmSVpBo4b7h4CXA5uAY8B7u/YM2LcGPUCSrUn2JNlz4sSJEcuQJA0yUrhX1dNVdaqqngM+zM+mXo4A6/t2vRQ4Os9j7KyquaqaW7NmzShlSJLmMVK4J1nbt3ojcPpKmvuBLUlekORyYCPw8HglSpKW6vxhOyT5BHA1cFGSI8A7gauTbKI35fIk8BaAqjqQZBfwGHASuK2qTk2kcknSvIaGe1XdPKD5IwvsvwPYMU5RkqTx+A1VSWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1KCh4Z7k7iTHk+zva7swye4kj3f3F/Rt257kcJJDSa6dVOGSpPktZuT+UeC6M9q2AQ9W1UbgwW6dJFcAW4Aru2M+mOS8ZatWkrQoQ8O9qh4CfnhG8w3APd3yPcAb+9rvrapnq+oJ4DCweXlKlSQt1qhz7pdU1TGA7v7irn0d8FTffke6NknSClruD1QzoK0G7phsTbInyZ4TJ04scxmSdG4bNdyfTrIWoLs/3rUfAdb37XcpcHTQA1TVzqqaq6q5NWvWjFiGJGmQUcP9fuCWbvkW4NN97VuSvCDJ5cBG4OHxSpQkLdX5w3ZI8gngauCiJEeAdwJ3AbuS3Ap8D7gJoKoOJNkFPAacBG6rqlMTql2SNI+h4V5VN8+z6Zp59t8B7BinKEnSePyGqiQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTp/nIOTPAk8A5wCTlbVXJILgX8BNgBPAr9XVT8ar0xJ0lIsx8j9N6tqU1XNdevbgAeraiPwYLcuSVpBk5iWuQG4p1u+B3jjBM4hSVrAuOFewOeT7E2ytWu7pKqOAXT3Fw86MMnWJHuS7Dlx4sSYZUiS+o015w5cVVVHk1wM7E7yncUeWFU7gZ0Ac3NzNWYdkqQ+Y43cq+pod38c+BSwGXg6yVqA7v74uEVKkpZm5HBP8stJXnx6GfhtYD9wP3BLt9stwKfHLVKStDTjTMtcAnwqyenH+XhVfTbJN4BdSW4FvgfcNH6ZkqSlGDncq+q7wGsGtP8PcM04RUmSxuM3VCWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDJhbuSa5LcijJ4STbJnUeSdLZJhLuSc4D/g74HeAK4OYkV0ziXJKks01q5L4ZOFxV362q/wXuBW6Y0LkkSWc4f0KPuw54qm/9CPCr/Tsk2Qps7VZ/kuTQiOe6CPjBiMdO0qzWBbNbm3UtjXUtzUzWlXePVdevzLdhUuGeAW31cytVO4GdY58o2VNVc+M+znKb1bpgdmuzrqWxrqU51+qa1LTMEWB93/qlwNEJnUuSdIZJhfs3gI1JLk/yfGALcP+EziVJOsNEpmWq6mSS24HPAecBd1fVgUmci2WY2pmQWa0LZrc261oa61qac6quVNXwvSRJq4rfUJWkBhnuktSgmQz3JOcl+WaSz3TrNyU5kOS5JPNeMpTkjiT7u33v7Gu/MMnuJI939xfMSF3vSvL9JPu62/XLVNd7knwnyaNJPpXkJfMcN/AnIibYX+PWNe3+ujvJ8ST7z2ifdn/NV9fU+ivJ+iRfTHKwe97f0bdtav01pK5p9tcvJnk4ybe6uv6ib9tI/TWT4Q7cARzsW98P/C7w0HwHJHkV8If0vh37GuANSTZ2m7cBD1bVRuDBbn0W6gJ4f1Vt6m4PLFNdu4FXVdWrgf8Atg+oa6GfiJhUf41bF0ypvzofBa4b0D61/hpSF0yvv04Cb62qVwKvA26bhefXkLpgev31LPBbVfUaYBNwXZLXddtG6q+ZC/cklwKvB/7hdFtVHayqYd9gfSXwtar6aVWdBP4duLHbdgNwT7d8D/DGGalrbPPU9fnuXABfo/c9gzMt9BMRk+qvcesa2xh1UVUPAT8csGma/bVQXWMbta6qOlZVj3TLz9ALu3Xd5qn115C6xjZGXVVVP+lWn9fdTl/tMlJ/zVy4Ax8A/gx4bonH7Qd+I8lLk7wQuJ6ffZHqkqo6Br3/XODiGakL4Pbu7drdI749HVbXm4F/G9A+6CciTj/JV6K/RqkLptdfC5lmfw0z9f5KsgF4LfD1rmkm+mtAXTDF/uqmc/YBx4HdVTVWf81UuCd5A3C8qvYu9diqOgi8m95boM8C36L3FmyW6/oQ8HJ6b8OOAe9dzrqSvKM718cGbR5U7lLOP4W6ptlfEzPBuqbeX0leBNwH3FlVP17K+adQ11T7q6pOVdUmeiP7zd2U7uiqamZuwF/TG6k9Cfw38FPgn/u2fwmYW+Rj/RXwR93yIWBtt7wWODQLdZ3RvgHYv1x1AbcAXwVeOM+xvwZ8rm99O7B90v01Tl3T7K+FzjvN/lpsf0yjv+hNLXwO+JNZ6q/56pp2f53xWO8E/nSc/lp04St9A64GPnNG25dYIESBi7v7y4DvABd06+8BtnXL24C/mZG61vbt88fAvctRF70P1x4D1iyw//nAd4HLgefTe0dx5ST7axnqmlp/9R131h/9NPtrSF3TfH4F+EfgAwO2TfP5tVBd0+yvNcBLuuVfAr4MvGGc/hqp8JW4ndE5N9J7RXwWeJpuZAe8DHig75gvd534LeCavvaX0vuU+fHu/sIZqeufgG8Dj9L77Z21y1TXYXrz1vu629/PU9f19D69/0/gHSvQX+PWNe3++gS9t+v/1/2/3zoj/TVfXVPrL+DX6U2nPdq33/XT7q8hdU2zv14NfLM7937gz8f9e/TnBySpQTP1gaokaXkY7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalB/w99+XKucIISCgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAARyklEQVR4nO3df5Dcd13H8efLFspPpbXXTkhSUzCgrTOkGAuKMpUytraOKQ5lUkeNTp3oWLQoo7Q4DvhHZorDrz+gOIVWMgotGSg0IgOUCFZnlJpigaShEmlsQ2ITBOSHM4WEt3/sN7K97N3t3d5uNp88HzM3u/vZ73f3lcvd6z732e9+L1WFJKktP3CiA0iSlp/lLkkNstwlqUGWuyQ1yHKXpAadfqIDAJx99tm1Zs2aEx1Dkk4q991331eqambQfVNR7mvWrGHnzp0nOoYknVSS/Odc97ksI0kNstwlqUGWuyQ1yHKXpAZZ7pLUIMtdkhpkuUtSgyx3SWqQ5S5JDZqKd6hKWtiaG/5uzvv23XTlBJPoZGC5S1NkvgKXFsNy11RyliqNxnLX48xVqhaqdHLxBVVJapAzd50wri9L4+PMXZIa5MxdY+XsXDoxLHdpAR65o5OR5S5NmL/NaBIsdwkLV+2x3DUyi1GaPh4tI0kNcuYujcAXWzWtnLlLUoOcuUuN87eLU5Mzd0lq0IIz9yRPAu4Bzui2f39VvS7JWcD7gDXAPuAVVfW1bp8bgWuBo8AfVNXHxpJepyRnopPh5/nkNszM/THgJVX1PGAdcHmSFwI3ADuqai2wo7tNkguAjcCFwOXAzUlOG0N2SdIcFiz36vlWd/MJ3UcBG4Ct3fhW4Kru+gbgjqp6rKoeAvYCFy9naEnS/IZ6QbWbed8H/Cjw9qr6dJJzq+ogQFUdTHJOt/lK4F/6dt/fjc1+zM3AZoDzzjtv6f8CqY9LCVLPUC+oVtXRqloHrAIuTvIT82yeQQ8x4DFvqar1VbV+ZmZmqLCSpOEs6miZqvo68Cl6a+mPJlkB0F0e6jbbD6zu220VcGDUoJKk4S1Y7klmkjyju/5k4KXAF4DtwKZus03AXd317cDGJGckOR9YC9y7zLklSfMYZs19BbC1W3f/AWBbVX04yT8D25JcCzwMXA1QVbuTbAMeAI4A11XV0fHElyQNsmC5V9XngIsGjP83cOkc+2wBtoycTtJYeUbPdvkOVUlqkOeWkcZkkrNiZ+CazXLXUFoojxb+DdKwXJaRpAY5cz8FOYOV2ufMXZIaZLlLUoMsd0lqkGvukpaVZ+acDs7cJalBlrskNchlGUmL5uG008+ZuyQ1yHKXpAa5LNMof22WTm3O3CWpQZa7JDXIcpekBlnuktQgy12SGuTRMicxj4iRNBdn7pLUoAXLPcnqJJ9MsifJ7iTXd+OvT/LlJPd3H1f07XNjkr1JHkxy2Tj/AZKk4w2zLHMEeHVVfSbJ04H7ktzd3feWqnpj/8ZJLgA2AhcCzwQ+keQ5VXV0OYNLkua24My9qg5W1We6698E9gAr59llA3BHVT1WVQ8Be4GLlyOsJGk4i1pzT7IGuAj4dDf0yiSfS3JbkjO7sZXAI3277WfAD4Mkm5PsTLLz8OHDi08uSZrT0OWe5GnAB4BXVdU3gHcAzwbWAQeBNx3bdMDuddxA1S1Vtb6q1s/MzCw2tyRpHkOVe5In0Cv291TVnQBV9WhVHa2q7wHv5PtLL/uB1X27rwIOLF9kSdJChjlaJsCtwJ6qenPf+Iq+zV4G7Oqubwc2JjkjyfnAWuDe5YssSVrIMEfLvAj4deDzSe7vxl4LXJNkHb0ll33A7wBU1e4k24AH6B1pc51HykgC/3j2JC1Y7lX1TwxeR//IPPtsAbaMkEuSNALfoSpJDbLcJalBlrskNchyl6QGWe6S1CDLXZIa5B/rmHL+QQ5JS+HMXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBnmcu6Sp4Lnel5czd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGrRguSdZneSTSfYk2Z3k+m78rCR3J/lid3lm3z43Jtmb5MEkl43zHyBJOt4wM/cjwKur6seBFwLXJbkAuAHYUVVrgR3dbbr7NgIXApcDNyc5bRzhJUmDLVjuVXWwqj7TXf8msAdYCWwAtnabbQWu6q5vAO6oqseq6iFgL3DxMueWJM1jUWvuSdYAFwGfBs6tqoPQ+wEAnNNtthJ4pG+3/d3Y7MfanGRnkp2HDx9eQnRJ0lyGLvckTwM+ALyqqr4x36YDxuq4gapbqmp9Va2fmZkZNoYkaQhDlXuSJ9Ar9vdU1Z3d8KNJVnT3rwAOdeP7gdV9u68CDixPXEnSMBY8K2SSALcCe6rqzX13bQc2ATd1l3f1jb83yZuBZwJrgXuXM7SkU8tcZ4z0bJFzG+aUvy8Cfh34fJL7u7HX0iv1bUmuBR4Grgaoqt1JtgEP0DvS5rqqOrrcwVsy36lOJWkpFiz3qvonBq+jA1w6xz5bgC0j5JIkjcB3qEpSgyx3SWqQ5S5JDbLcJalB/oHsCfGIGEmT5MxdkhpkuUtSgyx3SWqQ5S5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQg/1jHMvIPckiTNd/33L6brpxgkumz4Mw9yW1JDiXZ1Tf2+iRfTnJ/93FF3303Jtmb5MEkl40ruCRpbsMsy7wbuHzA+Fuqal338RGAJBcAG4ELu31uTnLacoWVJA1nwXKvqnuArw75eBuAO6rqsap6CNgLXDxCPknSEozyguork3yuW7Y5sxtbCTzSt83+buw4STYn2Zlk5+HDh0eIIUmabanl/g7g2cA64CDwpm48A7atQQ9QVbdU1fqqWj8zM7PEGJKkQZZU7lX1aFUdrarvAe/k+0sv+4HVfZuuAg6MFlGStFhLKvckK/puvgw4diTNdmBjkjOSnA+sBe4dLaIkabEWPM49ye3AJcDZSfYDrwMuSbKO3pLLPuB3AKpqd5JtwAPAEeC6qjo6luSSpDktWO5Vdc2A4Vvn2X4LsGWUUJKk0Xj6AUlqkKcfkNSkU/3UBJb7Inn+GEknA5dlJKlBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQb5DVdIp51Q4NYEzd0lqkOUuSQ2y3CWpQZa7JDXIF1QlqU8rL7Y6c5ekBlnuktQgy12SGmS5S1KDFiz3JLclOZRkV9/YWUnuTvLF7vLMvvtuTLI3yYNJLhtXcEnS3IaZub8buHzW2A3AjqpaC+zobpPkAmAjcGG3z81JTlu2tJKkoSx4KGRV3ZNkzazhDcAl3fWtwKeA13Tjd1TVY8BDSfYCFwP/vEx5J2K+Q6Ek6WSw1DX3c6vqIEB3eU43vhJ4pG+7/d3YcZJsTrIzyc7Dhw8vMYYkaZDlfkE1A8Zq0IZVdUtVra+q9TMzM8scQ5JObUst90eTrADoLg914/uB1X3brQIOLD2eJGkpllru24FN3fVNwF194xuTnJHkfGAtcO9oESVJi7XgC6pJbqf34unZSfYDrwNuArYluRZ4GLgaoKp2J9kGPAAcAa6rqqNjyi5JmsMwR8tcM8ddl86x/RZgyyihJEmj8R2qktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUGWuyQ1yHKXpAZZ7pLUoAXfoSpJ6pnvbz3su+nKCSZZmDN3SWqQ5S5JDTpll2X8U3qSWubMXZIaZLlLUoMsd0lqkOUuSQ2y3CWpQZa7JDXIcpekBlnuktQgy12SGjTSO1ST7AO+CRwFjlTV+iRnAe8D1gD7gFdU1ddGiylJWozlmLn/fFWtq6r13e0bgB1VtRbY0d2WJE3QOJZlNgBbu+tbgavG8BySpHmMWu4FfDzJfUk2d2PnVtVBgO7ynBGfQ5K0SKOeFfJFVXUgyTnA3Um+MOyO3Q+DzQDnnXfeiDEkSf1GmrlX1YHu8hDwQeBi4NEkKwC6y0Nz7HtLVa2vqvUzMzOjxJAkzbLkck/y1CRPP3Yd+AVgF7Ad2NRttgm4a9SQkqTFGWVZ5lzgg0mOPc57q+qjSf4V2JbkWuBh4OrRY0qSFmPJ5V5VXwKeN2D8v4FLRwm1XPxrS5JOVb5DVZIaZLlLUoMsd0lqkOUuSQ2y3CWpQaO+Q1WSxPxH5+276coJJulx5i5JDbLcJalBlrskNchyl6QGWe6S1CDLXZIaZLlLUoMsd0lqUBNvYvLUvpL0eM7cJalBlrskNchyl6QGWe6S1KAmXlCVpGl2Is4Y6cxdkhpkuUtSgyx3SWrQ2Mo9yeVJHkyyN8kN43oeSdLxxlLuSU4D3g78InABcE2SC8bxXJKk441r5n4xsLeqvlRV3wHuADaM6bkkSbOM61DIlcAjfbf3Ay/o3yDJZmBzd/NbSR5cwvOcDXxlSQkn72TKCuYdN/OO10mTN28Alp73R+a6Y1zlngFj9bgbVbcAt4z0JMnOqlo/ymNMysmUFcw7buYdL/OOb1lmP7C67/Yq4MCYnkuSNMu4yv1fgbVJzk/yRGAjsH1MzyVJmmUsyzJVdSTJK4GPAacBt1XV7jE81UjLOhN2MmUF846becfrlM+bqlp4K0nSScV3qEpSgyx3SWrQ1Jd7kuuT7EqyO8mr5tnup5IcTfLyCcYblGPBvEkuSXJ/t80/TDji7Czz5k3yQ0n+Nslnu21+a8L5bktyKMmuvrGzktyd5Ivd5Zlz7DvxU2AsNW+S1Uk+mWRP93m+fprz9m17WpJ/S/Lhac+b5BlJ3p/kC93n+aenPO8fdl8Lu5LcnuRJi3ryqpraD+AngF3AU+i9+PsJYO2A7U4D/h74CPDyac4LPAN4ADivu33OlOd9LfCG7voM8FXgiRPM+GLg+cCuvrG/AG7ort9wLN+Ar4n/AJ4FPBH4LHDBFOddATy/u/504N+nOW/ftn8EvBf48DR/PXT3bQV+u7v+ROAZ05qX3htBHwKe3N3eBvzmYp572mfuPw78S1X9b1UdAf4BeNmA7X4f+ABwaJLhBhgm768Cd1bVwwBVdSIzD5O3gKcnCfA0euV+ZFIBq+qe7jn7baD3jUp3edWAXU/IKTCWmreqDlbVZ7rr3wT20PsGH6sRPr8kWQVcCbxrXPlmW2reJD9Ir2hv7R7nO1X19bEF7Yzy+aU34XpyktPpTcAW9V6haS/3XcCLk/xwkqcAV/D4N0eRZCW9QvrLE5BvtgXzAs8BzkzyqST3JfmNiaf8vmHyvo3eD4EDwOeB66vqe5ONeZxzq+og9EoROGfANoNOgTH2spzDMHn/X5I1wEXAp8cfbaBh874V+BPgZPh6eBZwGPirbhnpXUmeOsmQfRbMW1VfBt4IPAwcBP6nqj6+mCeZ6nKvqj3AG4C7gY/S+9V69qzxrcBrquroZNMdb8i8pwM/SW/GcxnwZ0meM8mcxwyZ9zLgfuCZwDrgbd0saNoteAqMaZTkafR+C31VVX3jROeZS5JfAg5V1X0nOsuQTqe3PPKOqroI+Da9JZGp1K3DbwDOp/e999Qkv7aYx5jqcgeoqlur6vlV9WJ6v958cdYm64E7kuwDXg7cnOSqyab8viHy7gc+WlXfrqqvAPcAz5t0zmOGyPtb9JaRqqr20lsH/LFJ55zl0SQrALrLQUtb03QKjGHykuQJ9Ir9PVV15wTzzTZM3hcBv9x9390BvCTJ30wu4uMM+/Wwv6qO/Tb0fnplfyIMk/elwENVdbiqvgvcCfzMYp5k6ss9yTnd5XnArwC3999fVedX1ZqqWkPvP+z3qupDk855zEJ5gbuAn0tyercU8gJ666snxBB5HwYu7bY5F3gu8KVJZhxgO7Cpu76J3ud0tmk6BcaCebvXNG4F9lTVmyeYbZAF81bVjVW1qvu+2wj8fVUtama5jIbJ+1/AI0me2w1dSu/AhhNhmK/fh4EXJnlK97VxKYvtiXG/WjzqB/CP9P4TPgtc2o39LvC7A7Z9NyfwaJlh8wJ/3G2zi96v31Obl96vhB+nt96+C/i1Cee7nd6a43fpzb6uBX4Y2EHvt4wdwFl9WT/St+8V9I46+Q/gT6c5L/Cz9JaNPkdvGex+4IppzTvrMS5hckfLjPL1sA7Y2X2OPwScOeV5/xz4Qvd999fAGYt5bk8/IEkNmvplGUnS4lnuktQgy12SGmS5S1KDLHdJapDlLkkNstwlqUH/B4P63h7jXslsAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "la valeur moyenne sur les tirages réalisés est égale à  10.09704365923575\n",
      "l'incertitude associée à la mesure est égale à  0.31466617556376014\n",
      "l'incertitude théorique sur la distance focale est  0.3109126351029605\n"
     ]
    }
   ],
   "source": [
    "#on commence classiquement par importer la library numpy (sous l'alias np) pour le calcul numérique\n",
    "#son sous module numpy.random pour effectuer les tirages aléatoires selon des lois bien controlées\n",
    "#et la library matplotlib.pyplot (sous l'alias pl) pour la réalisation de graphique.\n",
    "import numpy as np\n",
    "import numpy.random as rd\n",
    "import matplotlib.pyplot as pl\n",
    "\n",
    "#on créée une liste de 10000 tirages aléatoires uniformément répartis sur l'intervalle étudié.\n",
    "#ce qu'on appelle généralement simulation Monte-Carlo.\n",
    "N=10000     # nombre de tirage aléatoire utilisé\n",
    "\n",
    "#Liste des valeurs pour X(objet)\n",
    "X_O=32.0    \n",
    "l_XO=0.5   \n",
    "X_O_MC=X_O+rd.uniform(-l_XO,l_XO,N)\n",
    "\n",
    "#Liste des valeurs pour X(Lentille)\n",
    "X_L=42.1    \n",
    "l_XL=0.2   \n",
    "X_L_MC=X_L+rd.uniform(-l_XL,l_XL,N)\n",
    "\n",
    "#Liste des valeurs pour f'\n",
    "f_MC=X_L_MC-X_O_MC\n",
    "\n",
    "#on peut alors réaliser les histogrammes pour visualiser les distributions.\n",
    "pl.figure(1)\n",
    "pl.hist(X_O_MC,bins='rice')\n",
    "pl.figure(2)\n",
    "pl.hist(X_L_MC,bins='rice')\n",
    "pl.figure(3)\n",
    "pl.hist(f_MC,bins='rice')\n",
    "#et en demander l'affichage.\n",
    "pl.show()\n",
    "\n",
    "#on évalue alors la distance focale pour cette expérience\n",
    "moyenne_f=np.average(f_MC)\n",
    "print(\"la valeur moyenne sur les tirages réalisés est égale à \",moyenne_f)\n",
    "#on évalue alors l'inceritude type \n",
    "u_f=np.std(f_MC,ddof=1)\n",
    "#on affiche l'incertitude associée\n",
    "print(\"l'incertitude associée à la mesure est égale à \",u_f)\n",
    "\n",
    "#on compare alors l'incertitude à la valeur numérique théorique. \n",
    "#Evaluation de type B de l'incertitude\n",
    "u_XO=l_XO/np.sqrt(3)\n",
    "u_XL=l_XL/np.sqrt(3)\n",
    "u_f_theo=np.sqrt(u_XO**2+u_XL**2)\n",
    "print(\"l'incertitude théorique sur la distance focale est \",u_f_theo)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7637365f",
   "metadata": {},
   "source": [
    "**Conclusion : On retiendra qu'on peut utiliser la formule de propagation des incertitude pour évaluer l'incertitude associée à une grandeur s'exprimant à l'aide d'une relation de type somme-différence en fonction des paramètres expérimentaux directement mesurés**\n",
    "\n",
    "**Pour Y s'exprimant par la formule générale de type somme-différence $Y=a*X_{1}+b*X_{2}$, l'incertitude type composée est exprimée par la relation $u^2\\left(Y\\right)=a^2*u^2\\left(X_{1}\\right)+b^2*u^2\\left(X_{2}\\right)$**\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "488ced46",
   "metadata": {},
   "source": [
    "## Evaluation de la distance focale de la lentille divergente par la méthode d'autocollimation\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3f9a28e3",
   "metadata": {},
   "source": [
    "Pour la lentille divergente, on a déjà mis en place la méthode d'autocolimation pour le doublet ce qui donne."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "f1f00e5f",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "la distance focale obtenue pour le doublet est alors 22.0\n",
      "l'incertitude associée est 0.3109126351029605\n"
     ]
    }
   ],
   "source": [
    "#pour l'objet\n",
    "X_Od=32.0             #position de l'objet\n",
    "l_XOd=0.5             #demi largeur d'intervalle\n",
    "u_XOd=l_XOd/np.sqrt(3)#incertitude type \n",
    "\n",
    "#pour le doublet de lentille\n",
    "X_Ld=54.0             #poisition du doublet\n",
    "l_Ld=0.2             #demi largeur d'intervalle\n",
    "u_Ld=l_Ld/np.sqrt(3) #incertitude type \n",
    "\n",
    "#pour la distance focale\n",
    "f_d=X_Ld-X_Od\n",
    "u_f_d=np.sqrt(u_XOd**2+u_Ld**2)\n",
    "\n",
    "print (\"la distance focale obtenue pour le doublet est alors\", f_d)\n",
    "print(\"l'incertitude associée est\", u_f_d)\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4500fca2",
   "metadata": {},
   "source": [
    "On a donc obtenue une estimation de la distance focale de la lentille convergente et de la distance focale du doublet, on peut donc obtenir une estimation de la distance focale de la lentille divergente en composant les incertitudes à l'aide de la formule $V_{doublet}=V_{conv}+V_{div}$ ce qui se traduit par $f'_{div}=\\frac{f'_{doublet}*f'_{conv}}{f'_{conv}-f'_{doublet}}$ qui n'est pas une formule simple. \n",
    "On va devoir composer les incertitudes avec la méthode Monte Carlo.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "a3a0472e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAASjUlEQVR4nO3df4xd513n8fcHN3W7aREJnWRd/9hExag4BVwxMpX6T2gK9jYIt0iBCVBZIsJFSqV2xWqJqRDlD0tB9AfSLqnkblOs3VKvBVSxws/UtCoFUuMUN43jmowaN5naik3Zqg0Ig93v/jHH2ht7Zu6duXPnjp95v6TRPec559z7nUf2Zx499/xIVSFJast3jbsASdLyM9wlqUGGuyQ1yHCXpAYZ7pLUIMNdkho0cLgnWZfk75M82q3fnOSxJM90rzf17LsvyXSS00l2jqJwSdL8FjNyfw9wqmf9AeBoVW0FjnbrJNkGTAF3ALuAh5KsW55yJUmDGCjck2wC7gb+Z0/zbuBgt3wQeHtP+6GqulhVzwLTwI5lqVaSNJCXDbjf7wD/DXh1T9utVXUOoKrOJbmla98IPN6z30zX9hJJ9gJ7AW688cYfef3rX7+4yiVpjXviiSf+saom5trWN9yT/CRwvqqeSHLnAJ+XOdquucdBVR0ADgBMTk7W8ePHB3hrSdIVSb4237ZBRu5vBn4qyduAVwDfneR/Ay8k2dCN2jcA57v9Z4DNPcdvAs4urXRJ0lL0nXOvqn1VtamqbmP2i9K/rKpfAI4Ae7rd9gCPdMtHgKkk65PcDmwFji175ZKkeQ065z6XB4HDSe4DngPuAaiqk0kOA08Dl4D7q+ry0JVKkgaW1XDLX+fcJWnxkjxRVZNzbfMKVUlqkOEuSQ0y3CWpQYa7JDXIcJekBg1zKqQ0Frc98Mfzbjvz4N0rWIm0ejlyl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIi5ik64QXb2kxDHdpFVkowKXFMNylEZkvqB1layU45y5JDeob7klekeRYki8lOZnkN7v29yf5epIT3c/beo7Zl2Q6yekkO0f5C0iSrjXItMxF4C1V9WKSG4DPJ/nTbtuHq+oDvTsn2QZMAXcArwU+neT7fUi2WuQcuVarviP3mvVit3pD97PQU7V3A4eq6mJVPQtMAzuGrlSSNLCB5tyTrEtyAjgPPFZVX+g2vTvJk0keTnJT17YReL7n8JmuTZK0QgYK96q6XFXbgU3AjiRvAD4CvA7YDpwDPtjtnrne4uqGJHuTHE9y/MKFC0soXZI0n0WdCllV30zyWWBX71x7ko8Cj3arM8DmnsM2AWfneK8DwAGAycnJhaZ5pJHzAiG1pm+4J5kA/r0L9lcCbwV+K8mGqjrX7fYO4Klu+Qjw+0k+xOwXqluBY8tfuq53Bqo0OoOM3DcAB5OsY3Ya53BVPZrkfyXZzuyUyxngXQBVdTLJYeBp4BJwv2fKSOPjH9G1qW+4V9WTwBvnaH/nAsfsB/YPV5okaam8/YC0wjw3XivB2w9IUoMMd0lqkNMyWpVGMXXR8nRIy7+blsaRuyQ1yJG7muIIVprlyF2SGuTIXSPlSPr65cVP1zfDXVrD/OPbLqdlJKlBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOe5a2ieKy2tPo7cJalBhrskNahvuCd5RZJjSb6U5GSS3+zab07yWJJnutebeo7Zl2Q6yekkO0f5C0iSrjXIyP0i8Jaq+mFgO7AryZuAB4CjVbUVONqtk2QbMAXcAewCHkqybgS1S5Lm0Tfca9aL3eoN3U8Bu4GDXftB4O3d8m7gUFVdrKpngWlgx3IWLUla2EBny3Qj7yeA7wN+t6q+kOTWqjoHUFXnktzS7b4ReLzn8Jmu7er33AvsBdiyZcvSfwNpxDwbSNejgb5QrarLVbUd2ATsSPKGBXbPXG8xx3seqKrJqpqcmJgYqFhJ0mAWdbZMVX0T+Cyzc+kvJNkA0L2e73abATb3HLYJODtsoZKkwQ1ytsxEku/pll8JvBX4CnAE2NPttgd4pFs+AkwlWZ/kdmArcGyZ65YkLWCQOfcNwMFu3v27gMNV9WiSvwUOJ7kPeA64B6CqTiY5DDwNXALur6rLoylfkjSXvuFeVU8Cb5yj/RvAXfMcsx/YP3R1kqQl8d4yGohnjEjXF28/IEkNMtwlqUFOy0hatIWm6c48ePcKVqL5OHKXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBnueul/A2A1IbDHdJy8oLnFYHp2UkqUGGuyQ1yGkZSSvGKZuV48hdkhpkuEtSgwZ5QPbmJJ9JcirJySTv6drfn+TrSU50P2/rOWZfkukkp5PsHOUvIEm61iBz7peAX6mqLyZ5NfBEkse6bR+uqg/07pxkGzAF3AG8Fvh0ku/3IdmStHL6jtyr6lxVfbFb/jZwCti4wCG7gUNVdbGqngWmgR3LUawkaTCLmnNPchvwRuALXdO7kzyZ5OEkN3VtG4Hnew6bYY4/Bkn2Jjme5PiFCxcWX7kkaV4Dh3uSVwF/CLy3qr4FfAR4HbAdOAd88Mqucxxe1zRUHaiqyaqanJiYWGzdkqQFDBTuSW5gNtg/UVV/BFBVL1TV5ar6DvBR/v/UywywuefwTcDZ5StZktTPIGfLBPgYcKqqPtTTvqFnt3cAT3XLR4CpJOuT3A5sBY4tX8mSpH4GOVvmzcA7gS8nOdG1/Rpwb5LtzE65nAHeBVBVJ5McBp5m9kyb+z1TRpJWVt9wr6rPM/c8+p8scMx+YP8QdUmShuAVqpLUIMNdkhpkuEtSgwx3SWqQ4S5JDfJhHZJWBR/ksbwcuUtSgxy5r0ELjZAktcGRuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDfI890Z5Lru0thnu1zEDXNJ8nJaRpAYN8oDszUk+k+RUkpNJ3tO135zksSTPdK839RyzL8l0ktNJdo7yF5AkXWuQkfsl4Feq6geANwH3J9kGPAAcraqtwNFunW7bFHAHsAt4KMm6URQvSZpb33CvqnNV9cVu+dvAKWAjsBs42O12EHh7t7wbOFRVF6vqWWAa2LHMdUuSFrCoOfcktwFvBL4A3FpV52D2DwBwS7fbRuD5nsNmurar32tvkuNJjl+4cGEJpUuS5jPw2TJJXgX8IfDeqvpWknl3naOtrmmoOgAcAJicnLxmu2Z5RoykpRho5J7kBmaD/RNV9Udd8wtJNnTbNwDnu/YZYHPP4ZuAs8tTriRpEIOcLRPgY8CpqvpQz6YjwJ5ueQ/wSE/7VJL1SW4HtgLHlq9kSVI/g0zLvBl4J/DlJCe6tl8DHgQOJ7kPeA64B6CqTiY5DDzN7Jk291fV5eUuXJI0v77hXlWfZ+55dIC75jlmP7B/iLokSUPwClVJapDhLkkNMtwlqUHeFVLSqjff9R5nHrx7hSu5fjhyl6QGGe6S1CCnZSRdtxa6Pcdan7Jx5C5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDVokAdkP5zkfJKnetren+TrSU50P2/r2bYvyXSS00l2jqpwSdL8Bhm5/x6wa472D1fV9u7nTwCSbAOmgDu6Yx5Ksm65ipUkDaZvuFfV54B/GvD9dgOHqupiVT0LTAM7hqhPkrQEw8y5vzvJk920zU1d20bg+Z59Zrq2ayTZm+R4kuMXLlwYogxJ0tWWGu4fAV4HbAfOAR/s2jPHvjXXG1TVgaqarKrJiYmJJZYhSZrLkh7WUVUvXFlO8lHg0W51Btjcs+sm4OySq1sjFnrggCQtxZJG7kk29Ky+A7hyJs0RYCrJ+iS3A1uBY8OVKElarL4j9ySfBO4EXpNkBvgN4M4k25mdcjkDvAugqk4mOQw8DVwC7q+qyyOpXJI0r77hXlX3ztH8sQX23w/sH6YoSdJwfEC2pCat9Ydne/sBSWqQ4S5JDTLcJalBhrskNcgvVFeIFypJWkmO3CWpQYa7JDXIcJekBjnnLmnNWQsXODlyl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDWob7gneTjJ+SRP9bTdnOSxJM90rzf1bNuXZDrJ6SQ7R1W4JGl+g4zcfw/YdVXbA8DRqtoKHO3WSbINmALu6I55KMm6ZatWkjSQQR6Q/bkkt13VvBu4s1s+CHwW+NWu/VBVXQSeTTIN7AD+dpnqXdW8ra+k1WKp95a5tarOAVTVuSS3dO0bgcd79pvp2q6RZC+wF2DLli1LLGPlGeCSrgfL/YVq5miruXasqgNVNVlVkxMTE8tchiStbUsN9xeSbADoXs937TPA5p79NgFnl16eJGkplhruR4A93fIe4JGe9qkk65PcDmwFjg1XoiRpsfrOuSf5JLNfnr4myQzwG8CDwOEk9wHPAfcAVNXJJIeBp4FLwP1VdXlEtUuS5jHI2TL3zrPprnn23w/sH6YoSdJwvEJVkhpkuEtSgwx3SWqQ4S5JDVrqFapN8ypUae1a6P//mQfvXsFKhuPIXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUFD3RUyyRng28Bl4FJVTSa5Gfg/wG3AGeBnqur/DlemJI3f9XTHyOUYuf9YVW2vqslu/QHgaFVtBY5265KkFTSKaZndwMFu+SDw9hF8hiRpAcOGewF/keSJJHu7tlur6hxA93rLkJ8hSVqkYZ/E9OaqOpvkFuCxJF8Z9MDuj8FegC1btgxZxuL5tCVJLRtq5F5VZ7vX88CngB3AC0k2AHSv5+c59kBVTVbV5MTExDBlSJKusuRwT3JjkldfWQZ+AngKOALs6XbbAzwybJGSpMUZZlrmVuBTSa68z+9X1Z8l+TvgcJL7gOeAe4YvU5K0GEsO96r6KvDDc7R/A7hrmKIkScPxClVJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSg4a9cZgkidX3IA9H7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBTZwK6fNQJemlmgh3SVrNxnEOvNMyktQgw12SGjSycE+yK8npJNNJHhjV50iSrjWScE+yDvhd4D8D24B7k2wbxWdJkq41qpH7DmC6qr5aVf8GHAJ2j+izJElXGdXZMhuB53vWZ4Af7d0hyV5gb7f6YpLTwGuAfxxRTa2wj/qzjwZjP/U38j7Kbw11+H+ab8Oowj1ztNVLVqoOAAdeclByvKomR1RTE+yj/uyjwdhP/V3PfTSqaZkZYHPP+ibg7Ig+S5J0lVGF+98BW5PcnuTlwBRwZESfJUm6ykimZarqUpJ3A38OrAMerqqTAxx6oP8ua5591J99NBj7qb/rto9SVf33kiRdV7xCVZIaZLhLUoPGHu5JfjvJV5I8meRTSb6na//eJJ9J8mKS/zHmMsduvn7qtu3rbvNwOsnOMZY5VknuSXIyyXeSTPa0vzzJx5N8OcmXktw5virHa4E+uiHJwa6PTiXZN846x2mBPvr5JCd6fr6TZPsYS13Q2MMdeAx4Q1X9EPAPwJV/VP8K/DrwX8dV2CozZz91t3WYAu4AdgEPdbd/WIueAn4a+NxV7b8EUFU/CPw48MEkq+Hf/jjM10f3AOu7PvoR4F1Jblvh2laLOfuoqj5RVdurajvwTuBMVZ1Y+fIGM/Z/4FX1F1V1qVt9nNlz4qmqf66qzzMb8mvefP3E7G0dDlXVxap6Fphm9vYPa05Vnaqq03Ns2gYc7fY5D3wTuC4vTBnWAn1UwI1JXga8Evg34FsrWtwqsUAf9boX+ORK1LNUYw/3q/wi8KfjLuI60NtPc93qYeOKV7S6fQnYneRlSW5ndmS6uc8xa80fAP8MnAOeAz5QVf803pJWtZ9llYf7ijyJKcmngf84x6b3VdUj3T7vAy4Bn1iJmlajJfZT31s9tGSQPprDw8APAMeBrwF/w2wfNmmJfbQDuAy8FrgJ+Kskn66qr46ozLFaYh9dOfZHgX+pqqdGUtwyWZFwr6q3LrQ9yR7gJ4G7ag2feL/EflpTt3ro10fzHHMJ+C9X1pP8DfDMcta1miylj4CfA/6sqv4dOJ/kr5mdumoy3JfYR1dMscpH7bAKpmWS7AJ+FfipqvqXcdezWi3QT0eAqSTruymHrcCxcdS4WiX5D0lu7JZ/HLhUVU+PuazV5jngLZl1I/Am4CtjrmnV6b6Iv4fZ25ivamO/QjXJNLAe+EbX9HhV/XK37Qzw3cDLmf0S7CfW6n/KPv30Pmbn4S8B762qNfm9RZJ3AP8dmGD238uJqtrZnfXx58B3gK8D91XV18ZV5zgt0EevAj7O7JfPAT5eVb89tkLHaL4+6rbdCTxYVW8aV32DGnu4S5KW39inZSRJy89wl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ36f4OWVmialm+1AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "la valeur moyenne sur les tirages réalisés est égale à  -18.67597685318012\n",
      "l'incertitude associée à la mesure est égale à  0.8864607127420258\n"
     ]
    }
   ],
   "source": [
    "#on a déjà créer les listes de tirages aléatoires pour la lentille convergente.\n",
    "#on crée maintenant les listes des positionspour le doublet\n",
    "X_Ld_MC=X_Ld+rd.uniform(-l_Ld,l_Ld,N)\n",
    "\n",
    "#Liste des valeurs pour la distance focale de la lentille divergente\n",
    "f_div_MC=(X_Ld_MC-X_O_MC)*(X_L_MC-X_O_MC)/((X_L_MC-X_O_MC)-(X_Ld_MC-X_O_MC))\n",
    "\n",
    "pl.hist(f_div_MC,bins='rice')\n",
    "pl.show()\n",
    "\n",
    "#on évalue alors la distance focale pour cette expérience\n",
    "moyenne_f_div=np.average(f_div_MC)\n",
    "print(\"la valeur moyenne sur les tirages réalisés est égale à \",moyenne_f_div)\n",
    "#on évalue alors l'inceritude type \n",
    "u_f_div=np.std(f_div_MC,ddof=1)\n",
    "#on affiche l'incertitude associée\n",
    "print(\"l'incertitude associée à la mesure est égale à \",u_f_div)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1e027311",
   "metadata": {},
   "source": [
    "On pet observer qu'on obtient bien une distance focale négative ce qui est logique puisque la lentille est divergente. On peut noter aussi l'allure de la distribution des valeurs aléatoire qui n'est plus du tout uniforme mais se centre de plus en plus sur la valeur moyenne. On peut noter enfin que l'incertitude sur cette valeur est grande."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4e4173bc",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
