{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"provenance":[],"authorship_tag":"ABX9TyNr1ejFfkSqlcUC1MoKJS9k"},"kernelspec":{"name":"python3","display_name":"Python 3"},"language_info":{"name":"python"}},"cells":[{"cell_type":"markdown","metadata":{"id":"BgHcopIjinRi"},"source":["activité 7.2 Cinétiques d’hydrolyse du 2-clhoro-2-méthylpropane\n","\n","3°)"]},{"cell_type":"code","metadata":{"colab":{"base_uri":"https://localhost:8080/","height":508},"id":"1sBdQPZacMId","executionInfo":{"status":"ok","timestamp":1695277461591,"user_tz":-120,"elapsed":990,"user":{"displayName":"Lionnel Malara","userId":"16342666118884215780"}},"outputId":"7905732b-a013-4b0c-e98d-5279d2810067"},"source":["import matplotlib.pyplot as plt #Pour tracer des graphiques\n","import numpy as np #Pour faire divers calculs\n","\n","## Données\n","\n","t=np.array([0, 30,60,80,100,120]) # valeur du temps en s\n","y=np.array([0,0.34,0.66,0.89,1.13,1.33]) # valeur de la fonction y\n","\n","# Courbe d'évolution de la fonction y=f(t) et régression linéaire\n","\n","\n","plt.figure(1)\n","p=np.polyfit(t,y,1)                     #Régression linéaire de y en fonction de t : équation de la forme y = p[0]t + p[1]\n","plt.plot(t,np.polyval(p,t),'r')         #Tracé de la droite de régression en rouge ; polyval permet de placer les ordonnées y prévues par la droite modèle.\n","plt.xlabel('t en s') #Légende de l’axe des abscisses\n","plt.ylabel('y') #Légende de l’axe des ordonnées\n","plt.plot(t,y,'b+',label='y=f(t)') #Représente y en fonction de x avec des croix bleues avec une étiquette pour la courbe\n","plt.title('vérification ordre 1') #Titre du graphe\n","plt.legend() #Affiche l’étiquette de la courbe\n","plt.grid() #Affiche le quadrillage\n","plt.show() #Affiche le graphique\n","\n","print('pente = ',format(p[0],\"#.2e\"),'/s ',' y_0 = ',format(p[1],\"#.2f\"))\t#Affiche le résultat\n","\n","t2=t[1:]\n","y2=y[1:]\n","k2=y2/t2\n","k_moy=np.mean(k2) #moyenne de k\n","u_kmoy= np.std(k2, ddof=1)/5**(1/2) #incertitude-type de k_moy\n","print (\"k_moy =\",format(k_moy,\"#.4e\"),\" avec u(k_moy) = \",format(u_kmoy,\"#.1e\"),'par calcul direct') #Affiche résultat et incertitude-type\n","\n"],"execution_count":2,"outputs":[{"output_type":"display_data","data":{"text/plain":["<Figure size 640x480 with 1 Axes>"],"image/png":"iVBORw0KGgoAAAANSUhEUgAAAjcAAAHHCAYAAABDUnkqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABZHElEQVR4nO3deVhUdf//8eewS4bkxuZaVqYVIqaZWlq43S4/u60MN7T0NpdvFmlFd4q2mVam3bepqbmU5dJiZqnhlpYmbtwtpmZpeiug3kaoJAzM+f1xcopAHRQ4M8PrcV1edc6cOfOedyqvPp/POcdmGIaBiIiIiJfwsboAERERkdKkcCMiIiJeReFGREREvIrCjYiIiHgVhRsRERHxKgo3IiIi4lUUbkRERMSrKNyIiIiIV1G4EfEydrudF198kY8//tjqUkRELKFwI+JlnnzySWbPns2tt97q8nvmzZuHzWbj4MGDhfa/9NJLXH311fj6+tKkSRMA6tWrx4ABA0qvYBeNGzcOm81W7p9bls7XdxG5PAo3Il7ko48+4u2332bVqlXUqFHjss712Wef8fjjj9OqVSvmzp3LCy+8UEpVnl9OTg7jxo1jw4YNZf5ZFU1qairDhg0jNjYWf39/rwuKIn/mZ3UBIlJ6Dh48yMqVK2nQoEGJ3tevXz/uv/9+AgMDnfvWrVuHj48Pc+bMISAgwLl/7969+PiUzf8X5eTkMH78eADatm1b6LWnn36aJ598skw+tyL49NNPmT17NjfffDNXX301+/bts7okkTKjkRsRLzJy5EiaNm3q8vFnzpwBwNfXl6CgoEL/N3/s2DEqVapUKNgABAYG4u/vXzoFl4Cfnx9BQUHl/rmXw+FwcPbs2cs6h2EY/Pbbb5ddy9ChQ/n111/Zvn077du3v+zzibgzhRsRN/Tee+9hs9n4/PPPi7w2c+ZMbDYb3377rXPfnj17uOeee6hatSpBQUE0a9aM5cuXF3rfufUdn3/+OcOGDaNmzZrUqlWr0Gvn1n7YbDbmzp3LmTNnsNls2Gw25s2bBxS/5iYrK4tHH32UevXqERgYSK1atejfvz8nTpwAIC8vj7FjxxIbG0uVKlW44ooraNOmDevXr3ee4+DBg86ptPHjxzs/d9y4cUDxa27y8/N59tlnueaaawgMDKRevXo89dRT5ObmFjquXr16dO3alS+++ILmzZsTFBTE1VdfzYIFC1z4r2GGwMcee4zatWsTGBjI9ddfz8svv4xhGIWOs9lsjBgxgoULF9K4cWMCAwNZtWoVAN999x133nknlSpVolatWjz33HM4HI4in3Wu1tWrV9OsWTMqVarEzJkznX1+5JFHnHU0aNCAiRMnFnuevwoLC6NSpUoufV8RT6dpKRE31KVLFypXrsySJUu44447Cr22ePFiGjduzI033giYPzRbtWpFVFQUTz75JFdccQVLliyhR48evP/++9x9992F3j9s2DBq1KjB2LFjnSM3f/XWW2/xxhtvkJqayuzZswG47bbbij329OnTtGnThu+//54HHniApk2bcuLECZYvX85///tfqlevTnZ2NrNnzyY+Pp7Bgwdz6tQp5syZQ8eOHUlNTaVJkybUqFGD6dOnM3ToUO6++27+/ve/A3DzzTeft0+DBg1i/vz53HPPPTz22GNs3bqVCRMm8P333/Phhx8WOnb//v3cc889PPjggyQkJPDmm28yYMAAYmNjady48Xk/wzAMunfvzvr163nwwQdp0qQJq1evZvTo0Rw5coRXX3210PHr1q1jyZIljBgxgurVq1OvXj0yMjJo164d+fn5zv9Gb7zxxnnDxt69e4mPj2fIkCEMHjyY66+/npycHO644w6OHDnCkCFDqFOnDps3byYpKYn09HSmTJly3u8gUuEYIuKW4uPjjZo1axr5+fnOfenp6YaPj4/xzDPPOPfdddddxk033WScPXvWuc/hcBi33Xabce211zr3zZ071wCM1q1bFzrnn187cOCAc19CQoJxxRVXFKmrbt26RkJCgnN77NixBmB88MEHRY51OByGYRhGfn6+kZubW+i1X375xQgLCzMeeOAB577jx48bgJGcnFzkXMnJycaf/8pKS0szAGPQoEGFjhs1apQBGOvWrStUM2Bs3LjRue/YsWNGYGCg8dhjjxX5rD9btmyZARjPPfdcof333HOPYbPZjP379zv3AYaPj4/x3XffFTr2kUceMQBj69athT6/SpUqRfp+rtZVq1YVOsezzz5rXHHFFca+ffsK7X/yyScNX19f49ChQxf8Hn82fPhwQ3/9izfTtJSIm+rVqxfHjh0rdOXQe++9h8PhoFevXgCcPHmSdevWcd9993Hq1ClOnDjBiRMn+N///kfHjh354YcfOHLkSKHzDh48GF9f31Kr8/333yc6OrrICBHgnEby9fV1rt1xOBycPHmS/Px8mjVrxs6dOy/pcz/99FMAEhMTC+1/7LHHAPjkk08K7W/UqBFt2rRxbteoUYPrr7+en3766aKf4+vry8MPP1zkcwzDYOXKlYX233HHHTRq1KjIOW699VaaN29e6PP79OlT7GfWr1+fjh07Ftq3dOlS2rRpw1VXXeX873zixAni4uIoKChg48aNF/weIhWJpqVE3FSnTp2oUqUKixcv5q677gLMKakmTZpw3XXXAeZUi2EYjBkzhjFjxhR7nmPHjhEVFeXcrl+/fqnW+eOPP9KzZ8+LHjd//nxeeeUV9uzZg91uv+x6fv75Z3x8fIpcGRYeHk5oaCg///xzof116tQpco6rrrqKX3755aKfExkZyZVXXllo/w033OB8/c+K+z4///wzLVq0KLL/+uuvL/YzizvHDz/8wNdff33eS/yPHTtW/BcQqYAUbkTcVGBgID169ODDDz/k9ddfJzMzky+//LLQ/WbOLSQdNWpUkf/TP+evP/ytWFT69ttvM2DAAHr06MHo0aOpWbMmvr6+TJgwgR9//PGyzu3q/VrON1pl/GVR8OUqjf4Wdw6Hw0H79u15/PHHi33PucArIgo3Im6tV69ezJ8/n7Vr1/L9999jGIZzSgrg6quvBsDf35+4uDhLarzmmmsKXblVnPfee4+rr76aDz74oFAYSU5OLnRcSW4sV7duXRwOBz/88INzFAUgMzOTrKws6tat6/K5LvY5a9as4dSpU4VGb/bs2eN83ZVz/PDDD0X279271+U6rrnmGk6fPm3Zf2cRT6I1NyJuLC4ujqpVq7J48WIWL15M8+bNC01Z1KxZk7Zt2zJz5kzS09OLvP/48eNlXmPPnj35z3/+U+TqJPhjVOTcqMmfR0m2bt3Kli1bCh0fHBwMmJc8X8zf/vY3gCJXCU2ePBkwrzgrDX/7298oKCjg3//+d6H9r776Kjabjc6dO7t0jq+++orU1FTnvuPHj7Nw4UKX67jvvvvYsmULq1evLvJaVlYW+fn5Lp9LxNtp5EbEjfn7+/P3v/+dRYsWcebMGV5++eUix0ybNo3WrVtz0003MXjwYK6++moyMzPZsmUL//3vf/nPf/5TpjWOHj2a9957j3vvvZcHHniA2NhYTp48yfLly5kxYwbR0dF07dqVDz74gLvvvpsuXbpw4MABZsyYQaNGjTh9+rTzXJUqVaJRo0YsXryY6667jqpVq3LjjTc6L3v/s+joaBISEnjjjTfIysrijjvuIDU1lfnz59OjRw/atWtXKt+vW7dutGvXjn/+858cPHiQ6OhoPvvsMz766CMeeeQRrrnmmoue4/HHH+ett96iU6dOjBw50nkpeN26dfn6669dqmP06NEsX76crl27Oi9hP3PmDN988w3vvfceBw8epHr16ud9/88//8xbb70FwPbt2wF47rnnAHNkqV+/fi7VIeIRrLxUS0QuLiUlxQAMm81mHD58uNhjfvzxR6N///5GeHi44e/vb0RFRRldu3Y13nvvPecx5y733rZtW5H3X86l4IZhGP/73/+MESNGGFFRUUZAQIBRq1YtIyEhwThx4oRhGOYl4S+88IJRt25dIzAw0IiJiTFWrFhhJCQkGHXr1i10rs2bNxuxsbFGQEBAocvC/3opuGEYht1uN8aPH2/Ur1/f8Pf3N2rXrm0kJSUVuiz+XM1dunQp8l3uuOMO44477iiy/69OnTplPProo0ZkZKTh7+9vXHvttcZLL73kvNT9HMAYPnx4sef4+uuvjTvuuMMICgoyoqKijGeffdaYM2dOsZeCF1fruTqSkpKMBg0aGAEBAUb16tWN2267zXj55ZeNvLy8C36H9evXG0Cxv1zpgYgnsRlGKa+mExEREbGQ1tyIiIiIV1G4EREREa+icCMiIiJeReFGREREvIql4Wbjxo1069aNyMhIbDYby5Ytc/m9X375JX5+fjRp0qTM6hMRERHPY2m4OXPmDNHR0UybNq1E78vKyqJ///7O5+2IiIiInOM2l4LbbDY+/PBDevTocdFj77//fq699lp8fX1ZtmwZaWlpLn+Ow+Hg6NGjXHnllSW61buIiIhYxzAMTp06RWRkJD4+Fx6b8bg7FM+dO5effvqJt99+23l3zQvJzc0lNzfXuX3kyBEaNWpUliWKiIhIGTl8+DC1atW64DEeFW5++OEHnnzySTZt2oSfn2ulT5gwgfHjxxfZP3v2bOdzbERERMS95eTkMGjQoEIPsD0fjwk3BQUF9O7dm/Hjx3Pddde5/L6kpCQSExOd29nZ2dSuXZsePXoQEhJSqjXa7XZSUlJo3749/v7+pXpub6NeuU69cp165Tr1qmTUL9eVVa+ys7MZNGiQS0tKPCbcnDp1iu3bt7Nr1y5GjBgBmOtnDMPAz8+Pzz77jDvvvLPI+wIDAwkMDCyy39/fv8x+g5blub2NeuU69cp16pXr1KuSUb9cV9q9Ksm5PCbchISE8M033xTa9/rrr7Nu3Tree+896tevb1FlIiIi4k4sDTenT59m//79zu0DBw6QlpZG1apVqVOnDklJSRw5coQFCxbg4+PDjTfeWOj9NWvWJCgoqMh+ERERqbgsDTfbt2+nXbt2zu1za2MSEhKYN28e6enpHDp0yJLaCgoKsNvtJXqP3W7Hz8+Ps2fPUlBQUEaVeQddhi8iImXF0nDTtm1bLnSbnXnz5l3w/ePGjWPcuHGlWpNhGGRkZJCVlXVJ7w0PD+fw4cP64X0RNpvtovcpEBERuRQes+amvJwLNjVr1iQ4OLhEIcXhcHD69GkqV66sH9wX4HA4OHLkCKGhoRcMtyIiIpdC4eZPCgoKnMGmWrVqJX6/w+EgLy+PoKAghZuLqFGjBr/++qum70REpNTpJ/CfnFtjo5v7lT1/f39sNpvCjYiIlDqFm2JovUzZO9djTUuJiEhpU7gRERERr6JwI0Xk5OTQs2dPQkJCsNlszivH9u7dS3h4OKdOnbrg+2+99Vbef//9cqhURESkKIUbKWL+/Pls2rSJzZs3k56eTpUqVQDzOV3/93//53xo2bx58wgNDS3y/qeffponn3wSh8NRnmWLiIgACjdlKj0dxo0z/+lJfvzxR2644QZuvPFGwsPDsdlsHDp0iBUrVjBgwICLvr9z586cOnWKlStXln2xIiLiVjJW7GDJvHqW/uxTuClD6ekwfnzZh5sFCxZQrVo1cnNzC+3v0aMH/fr1K9G52rZtyyuvvMLGjRux2Wy0bdsWgCVLlhAdHU1UVBQAGzZsYODAgfz666/YbDZsNpvzhoq+vr787W9/Y9GiRZf93URExEM4HDBpEsfu+T/eWRZNxlHrLhhRuPEC9957LwUFBSxfvty579ixY3zyySc88MADbNq0icqVK1/w18KFCwH44IMPGDx4MC1btiQ9PZ0PPvgAgE2bNtGsWTPn+W+77TamTJlCSEgI6enppKenM2rUKOfrzZs3Z9OmTeXUARERsdTx49C1KzzxBDbH77f4yMuzrBzdxK+UZWTY2L8ffHxg505z37l/AkREmL9KU6VKlejduzdz587l3nvvBeDtt9+mTp06tG3blrNnz5KWlnbBc4SFhQFQtWpVgoODCQgIIDw83Pn6zz//XCjcBAQEUKVKFWw2W6HjzomMjOTw4cM4HA7d0FBExJtt3Ej6fSNJz7RBwK1su3cCLIRd3wfhV9k8pCx+9l2Iwk0pmzcvkIkTC/8wHzz4j39PTjbX4ZS2wYMHc8stt3DkyBGioqKYN28eAwYMwGazUalSJRo0aHBZ5//tt98ICgpy+fhKlSrhcDjIzc2lUqVKl/XZIiLihgoKYMIESE5mpmMM4xkHeYA5EcBDD/0RMcrqZ9/5KNyUsgEDcrnnngB8fHzYudMMNrNmQdOm5utllVxjYmKIjo5mwYIFdOjQge+++45PPvkEMKeUOnfufMH3z5w5kz59+pz39erVq/PLL7+4XM/Jkye54oorFGxERLxRZib06QNr1wIw5J6TdB+ZA8HBbNuWz0MP+TFjRj633GLGjPIctQGFm1IXHm4QEmJOS53TtOkf4aYsDRo0iClTpnDkyBHi4uKoXbs2AM2aNXN5Wup8YmJi2L17d6F9AQEB5318wrfffktMTIzrxYuIiGdYu9YMNpmZEBwM06YRMWAA5/JLfr65kDgmxiiXn33FUbjxIr1792bUqFHMmjWLBQsWOPeXxrRUx44dGTRoEAUFBfj6+gJQr149Tp8+zdq1a4mOjiY4ONj5XK5NmzbRoUOHy/pMERFxIwUF8Mwz8OyzYBjQuDEsWQKNGlldWRFa6VmGIiLMecbyGo6rUqUKPXv2pHLlyvTo0aNUz925c2f8/PxYs2aNc99tt93GQw89RK9evahRowaTJk0C4MiRI2zevJmBAweWag0iImKRo0fhrrvMcGMYMGgQpKYWG2zCw6FXrz0Uc61JudHITRmKiCjfBVRgBos+ffoQGBh4yeeYMmVKkX1+fn489dRTTJ48mY4dOzr3T58+nenTpxc69rXXXmPAgAHUqlXrkmsQERE3sXo19O0LJ05A5cowcyb07n3ewyMiID5+LxER15RjkYUp3HiJX375hQ0bNrBhwwZef/31MvmMIUOGkJWVxalTp5yPYChOzZo1SUxMLJMaRESknOTnw5gx8OKL5nZ0tDkNdd111tblAoUbLxETE8Mvv/zCxIkTuf7668vkM/z8/PjnP/950eMee+yxMvl8EREpJ4cPQ3w8fPmluT10KEyeDCW4JYiVFG68xMGDB60uQUREvMGKFZCQACdPQkiIeT+T++6zuqoS0YJiERERMR+X8Nhj0K2bGWxiY81b7HtYsAGN3BTLMKx72FdFca7HNpvN4kpERISDB6FXL/MKKICRI2HiRLiMi1OspJGbP/H39wcgJyfH4kq8n91uxzAM5z1zRETEIh9+CDExZrAJDTW3p0zx2GADGrkpxNfXl9DQUI4dOwZAcHBwiUYWHA4HeXl5nD17Vg+LvACHw8Hx48fJyclRuBERsUpuLoweDf/6l7ndogUsWgT16llaVmlQuPmLc0+4PhdwSsIwDH777TcqVaqk6ZaLsNls/Prrr+qTiIgV9u83p6F27jS3R42CF16A32cwPJ3CzV/YbDYiIiKoWbMmdru9RO+12+1s3LiR22+/3TnFJcWz2Wzs3bvX6jJERCqeJUvMOwyfOgXVqsH8+dCli9VVlSqFm/Pw9fUt8ZSJr68v+fn5BAUFKdxcREmDo4iIXKbffoPERJgxw9xu3RrefRe88G7yWhgiIiLi7fbuhVtv/SPYJCXB+vVeGWxAIzciIiLe7e234aGH4MwZqFED3noL/vSMQG+kkRsRERFvlJMDDz4I/fqZwaZtW0hL8/pgAwo3IiIi3mf3brjlFnjzTbDZIDkZ1qyByEirKysXmpYSERHxFoYB8+bB8OHmAuLwcFi4EO680+rKypXCjYiIiDc4fRqGDTPX1ADExZnrbcLCrK3LApqWEhER8XRffw3NmpnBxscHnnsOVq+ukMEGNHIjIiLiuQwDZs2Chx82H6cQFWXeu6ZNG6srs5TCjYiIiCfKzoYhQ8znQQF07gwLFkD16tbW5QY0LSUiIuJpdu6E2Fgz2Pj6wqRJsGKFgs3vNHIjIiLiKQwDpk2Dxx6DvDyoU8cMOC1bWl2ZW1G4ERER8QRZWeZN+T74wNzu3h3mzoWqVS0tyx1ZOi21ceNGunXrRmRkJDabjWXLll3w+A8++ID27dtTo0YNQkJCaNmyJatXry6fYkVERKySmgoxMWaw8feHKVNg2TIFm/OwNNycOXOG6Ohopk2b5tLxGzdupH379nz66afs2LGDdu3a0a1bN3bt2lXGlYqIiFjAMODVV80neB88CPXrw5dfwsiR5p2HpViWTkt17tyZzp07u3z8lClTCm2/8MILfPTRR3z88cfExMSUcnUiIiIWOnkSBgyAjz82t3v2hNmzITTUyqo8gkdfLeVwODh16hRVNSwnIiLeZPNmaNLEDDYBAeYi4qVLFWxc5NELil9++WVOnz7Nfffdd95jcnNzyc3NdW5nZ2cDYLfbsdvtpVrPufOV9nm9kXrlOvXKdeqV69Srkim3fjkc+LzyCj5jx2IrKMBo0ID8hQvN9Tb5+WX72aWkrHpVkvPZDMMwSvXTL5HNZuPDDz+kR48eLh3/zjvvMHjwYD766CPi4uLOe9y4ceMYP358se8PDg6+1HJFRERKVcCvv9J06lTCdu4E4L9t2vCfYcPIr1TJ4srcQ05ODr179+bXX38lJCTkgsd6ZLhZtGgRDzzwAEuXLqVLly4XPLa4kZvatWtz4sSJizanpOx2OykpKbRv3x5/f/9SPbe3Ua9cp165Tr1ynXpVMmXdL9umTfj264ft6FGMoCAKXn0V44EHPHLRcFn1Kjs7m+rVq7sUbjxuWurdd9/lgQceYNGiRRcNNgCBgYEEBgYW2e/v719mf6DL8tzeRr1ynXrlOvXKdepVyZR6vwoKYMIESE4GhwMaNsS2ZAl+N91Uep9hkdLuVUnOZWm4OX36NPv373duHzhwgLS0NKpWrUqdOnVISkriyJEjLFiwADCnkhISEpg6dSotWrQgIyMDgEqVKlGlShVLvoOIiMglyciAvn1h7Vpzu39/c+Fw5crW1uUFLL1aavv27cTExDgv405MTCQmJoaxY8cCkJ6ezqFDh5zHv/HGG+Tn5zN8+HAiIiKcv0aOHGlJ/SIiIpdk7Vrzaqi1ayE4GObNg/nzFWxKiaUjN23btuVCS37mzZtXaHvDhg1lW5CIiEhZKiiAZ56BZ581b9DXuDEsWQKNGlldmVfxuDU3IiIiHunoUejdGz7/3NweNAimTjVHbqRUKdyIiIiUtVWroF8/OHHCnHqaOdMMOlImPPoOxSIiIm4tPx+SkqBzZzPYREfDjh0KNmVMIzciIiJl4fBhiI83H3QJMGwYvPIKBAVZW1cFoHAjIiJS2lasgIQE8+GXISHmAy/vvdfqqioMTUuJiIiUlrw8eOwx6NbNDDaxsbBzp4JNOdPIjYiISGk4cADuvx9SU83tkSNh4kQo5i75UrY0ciMiInIR6enw7rvXk55+ngM+/NB8cndqKoSGmttTpijYWEThRkRE5CIyMmDx4ob8/tSfP+TmwsMPw9//Dr/+Ci1aQFoauPAQaCk7CjciIiKXYv9+uO02+Ne/zO1Ro2DTJqhb19q6RGtuREREipOejnMaatcum/Offn7AZ58R8dxwIs7sh2rVzOdCdeliXbFSiMKNiIhIMWbOhPHjz22ZPy4feujcj80OJNOHca3XwrvvQq1aVpQo56FpKRERkWIMGWLeTHjHDpgxIx+AN6LGsYOm7CCWIQ8Hwfr1CjZuSCM3IiIixYiIMH8B5C//FOhO7JHlNK3xX3j7bejQwdL65PwUbkRERM4nJwf+7//we3MX0B1H01j4eAVERlpdmVyAwo2IiEhxvvsO7rsPdu8mnAiGN36Pmu/9CyL1bCh3p3AjIiLyZ4YB8+bB8OHw228QHk7N+fNo/9tvRNTytbo6cYEWFIuIiJxz+rT5wMsHHjCDTfv2kJaG0a6d1ZVJCSjciIiIAHz9NTRrBm+9BT4+8NxzsGoVhIVZXZmUkKalRESkYjMMeOMN80GXubkQFWXeu6ZNG6srk0ukcCMiIhVXdjb84x+weLG53bkzLFgA1atbW5dcFk1LiYhIxbRzJ8TGmsHG1xcmTYIVKxRsvIBGbkREpGIxDJg2DR57DPLyoE4dWLQIWra0ujIpJQo3IiJScWRlwYMPwgcfmNvdu8PcuVC1qqVlSenStJSIiFQMqakQE2MGG39/mDIFli1TsPFCCjciIuLdDANefRVat4aDB6F+ffjyS/PqKJvN6uqkDGhaSkREvNfJkzBgAHz8sbndsyfMng2hoVZWJWVMIzciIuKdNm+GJk3MYBMQYC4iXrpUwaYCULgRERHv4nDAxIlw++1w+DA0aABffQXDhmkaqoLQtJSIiHiP48fNZ0OtXGlux8fDzJlw5ZXW1iXlSuFGRES8w8aNZpg5ehSCguC112DQII3WVECalhIREc9WUGA+5LJdOzPYNGxoXvY9eLCCTQWlkRsREfFcmZnQpw+sXWtu9+9vLhyuXNnausRSCjciIuKZ1q41g01mJgQHw+uvm+ttpMLTtJSIiHiWggJITob27c1g07gxbNumYCNOGrkRERHPcfQo9O4Nn39ubg8aBFOnmiM3Ir9TuBEREc+wejX07QsnTphrambONIOOyF9oWkpERNxbfj4kJUGnTmawiY6GHTsUbOS8NHIjIiLu6/Bh8941X35pbg8bBq+8Yt7HRuQ8LB252bhxI926dSMyMhKbzcayZcsu+p4NGzbQtGlTAgMDadCgAfPmzSvzOkVExAIrVpjPhvrySwgJgSVLzMu8FWzkIiwNN2fOnCE6Oppp06a5dPyBAwfo0qUL7dq1Iy0tjUceeYRBgwaxevXqMq5URETKTV4ePPYYdOtmPtU7NhZ27oR777W6MvEQlk5Lde7cmc6dO7t8/IwZM6hfvz6vvPIKADfccANffPEFr776Kh07diyrMkVEpLwcPAj33w9bt5rbI0eaD8EMDLS0LPEsHrWgeMuWLcTFxRXa17FjR7Zs2WJRRSIiUmo+/BBiYsxgExpqbk+ZomAjJeZRC4ozMjIICwsrtC8sLIzs7Gx+++03KlWqVOQ9ubm55ObmOrezs7MBsNvt2O32Uq3v3PlK+7zeSL1ynXrlOvXKdW7Vq9xcfJ58Et/flyg4mjenYOFCqFsX3KE+3Kxfbq6selWS83lUuLkUEyZMYPz48UX2f/bZZwSX0U2fUlJSyuS83ki9cp165Tr1ynVW9yo4PZ1bXn6Z0B9/BOCHHj34vm9fjO++g+++s7S24ljdL09S2r3Kyclx+ViPCjfh4eFkZmYW2peZmUlISEixozYASUlJJCYmOrezs7OpXbs2HTp0ICQkpFTrs9vtpKSk0L59e/z9/Uv13N5GvXKdeuU69cp17tAr29Kl+D7+OLZTpzCqVaNgzhzq/e1v1LOkmgtzh355irLq1bmZF1d4VLhp2bIln376aaF9KSkptGzZ8rzvCQwMJLCY+Vp/f/8y+w1aluf2NuqV69Qr16lXrrOkV7/9BomJMGOGud26NbZ338WvVq3yreMS6PeW60q7VyU5l6ULik+fPk1aWhppaWmAeal3Wloahw4dAsxRl/79+zuPf+ihh/jpp594/PHH2bNnD6+//jpLlizh0UcftaJ8EREpqb174dZb/wg2SUmwfj14QLARz2HpyM327dtp166dc/vc9FFCQgLz5s0jPT3dGXQA6tevzyeffMKjjz7K1KlTqVWrFrNnz9Zl4CIinmDhQhgyBM6cgRo14K23QH9/SxmwNNy0bdsWwzDO+3pxdx9u27Ytu3btKsOqRESkVOXkwMMPw5w55nbbtmbQiYy0tCzxXh51nxsREfEwu3dD8+ZmsLHZIDkZ1qxRsJEy5VELikVExEMYBsybB8OHmwuIw8PN0Zo777S6MqkAFG5ERKR0nT5tPr37rbfM7fbtzX//y01YRcqKpqVERKT0fP013HKLGWZ8fOC552DVKgUbKVcauRERkctnGDBrlvmgy7NnISoK3n0X2rSxujKpgBRuRETk8mRnm5d4L1pkbnfuDAsWQPXq1tYlFZampURE5NLt2gWxsWaw8fWFSZNgxQoFG7GURm5ERKTkDANef918jEJeHtSpYwacCzwOR6S8KNyIiEjJZGXBoEHw/vvmdvfuMHcuVK1qaVki52haSkREXJeaCjExZrDx94cpU2DZMgUbcSsKNyIicnGGAa++Cq1bw8GDUL8+fPmleXWUzWZ1dSKFaFpKREQu7ORJGDgQli83t3v2hNmzITTU0rJEzkcjNyIicn6bN0OTJmawCQiAadNg6VIFG3FrCjciIlKUw2Fe1n377XD4MDRoAF99ZT5WQdNQ4uY0LSUiIoUdPw4JCbBypbkdHw8zZ8KVV1pbl4iLFG5EROQPGzeaYeboUQgKgtdeMy/71miNeBBNS4mICBQUmA+5bNfODDYNG5qXfQ8erGAjHkfhRkSkAkpPh3ffvZ70dCAzEzp1gjFjzLU2/fvDtm1w001WlylySRRuREQqoIwMWLy4IRkfb4foaFizBoKDYd48mD8fKle2ukSRS6Y1NyIiFVGBAwDfESOATGjcGJYsgUaNrK1LpBQo3IiIVBDp6eYvjh8n7eF5wGh2EYOtRw8YNYqIqyoRYW2JIqVC4UZEpIKYORPGjweoAYwGYDCzYRmwDJKTYdw4q6oTKT0KNyIiFUF+PkP+N4nuvAfAjqju/OPIOGbMyOeWW8wfBREathEvoXAjIuLtDh+G+HgivvzSnHYaOpT8+NFwO8TEGDRtanWBIqVLV0uJiHizTz4xnw315ZcQEgKLF8Prr0NggNWViZQZhRsREW+UlwejRkHXruZTvWNjYedOuO8+AMLDoVevPYSHW1ynSBlQuBER8TYHD5oPvHzlFXN75Ehz5Oaaa5yHRERAfPxerbMRr6Q1NyIi3mTZMhg4ELKyIDQU5s6FHj2srUmknGnkRkTEG+TmmiM0d99tBpsWLSAtTcFGKiSFGxERT/fjj9CqlfkEbzDX2mzaBHXrWluXiEU0LSUi4smWLIFBg+DUKahWzXwuVJcuVlclYimN3IiIeKKzZ2HoUOjVyww2rVub01AKNiIKNyIiHmffPrj1Vpgxw9xOSoL166FWLWvrEnETmpYSEfEkCxfCkCFw5gzUqAFvvQUdO1pdlYhb0ciNiIgnyMkx19b07WsGm7ZtzWkoBRuRIhRuRETc3e7d0Lw5zJkDNpv5+O41ayAy0urKRNySpqVERNzZvHkwfLg5chMebk5L3Xmn1VWJuDWN3IiIuKPTpyEhwbzbcE4OxMWZ01AKNiIXpXAjIuJuvv4abrkFFiwAHx947jlYvRrCwqyuTMQjaFpKRMRdGAbMmmU+RuHsWYiKgnffhTZtrK5MxKNYPnIzbdo06tWrR1BQEC1atCA1NfWCx0+ZMoXrr7+eSpUqUbt2bR599FHOnj1bTtWKiJSR7Gzo3du8zPvsWejc2ZyGUrARKTFLw83ixYtJTEwkOTmZnTt3Eh0dTceOHTl27Fixx7/zzjs8+eSTJCcn8/333zNnzhwWL17MU089Vc6Vi4iUol27IDYWFi0CX1+YNAlWrIDq1a2uTMQjWRpuJk+ezODBgxk4cCCNGjVixowZBAcH8+abbxZ7/ObNm2nVqhW9e/emXr16dOjQgfj4+IuO9oiIuCXDgGnTzLsN798PdeqYD7wcPdpcayMil8SyNTd5eXns2LGDpKQk5z4fHx/i4uLYsmVLse+57bbbePvtt0lNTaV58+b89NNPfPrpp/Tr1++8n5Obm0tubq5zOzs7GwC73Y7dbi+lb4PznH/+p5yfeuU69cp1HtWrrCx8hwzB58MPAXB07UrB7NlQtSqUQ/0e1Ss3oH65rqx6VZLz2QzDMEr101109OhRoqKi2Lx5My1btnTuf/zxx/n888/ZunVrse977bXXGDVqFIZhkJ+fz0MPPcT06dPP+znjxo1j/PjxRfa/8847BAcHX/4XEREpodAffqDZyy9zRWYmDj8/vktI4KeuXc0b9IlIsXJycujduze//vorISEhFzzWo66W2rBhAy+88AKvv/46LVq0YP/+/YwcOZJnn32WMWPGFPuepKQkEhMTndvZ2dnUrl2bDh06XLQ5JWW320lJSaF9+/b4+/uX6rm9jXrlOvXKdW7fK8PA51//wuepp7DZ7Rj16+NYuJCGzZrRsJxLcfteuRn1y3Vl1atzMy+usCzcVK9eHV9fXzIzMwvtz8zMJDw8vNj3jBkzhn79+jFo0CAAbrrpJs6cOcM//vEP/vnPf+JTzBx1YGAggYGBRfb7+/uX2W/Qsjy3t1GvXKdeuc4te3XypHlDvuXLze2ePbHNno1faKilZbllr9yY+uW60u5VSc5l2Yq1gIAAYmNjWbt2rXOfw+Fg7dq1haap/iwnJ6dIgPH19QXAotk1EZGL27wZmjQxg01AgLmIeOlSsDjYiHgrS6elEhMTSUhIoFmzZjRv3pwpU6Zw5swZBg4cCED//v2JiopiwoQJAHTr1o3JkycTExPjnJYaM2YM3bp1c4YcERG34XDAyy/DU09BQQE0aABLlkBMjNWViXg1S8NNr169OH78OGPHjiUjI4MmTZqwatUqwn6/xfihQ4cKjdQ8/fTT2Gw2nn76aY4cOUKNGjXo1q0bzz//vFVfQUSkeMePm8+GWrnS3I6Ph5kz4corra1LpAKwfEHxiBEjGDFiRLGvbdiwodC2n58fycnJJCcnl0NlIiKXaONGM8wcPQpBQfDaazBokK6GEiknukuUiEhpKSgwH3LZrp0ZbBo2hNRUGDxYwUakHFk+ciMi4hUyM6FvX1izxtzu399cOFy5srV1iVRACjciIpdr3TrzoZeZmRAcDK+/bq63ERFLaFpKRORSFRRAcjLExZnBpnFj2LZNwUbEYhq5ERG5FEePQp8+cO7Ch0GDYOpUc+RGRCylcCMiUlKrV0O/fubl3pUrm5d49+5tdVUi8jtNS4mIuCo/H5KSoFMnM9hER8OOHQo2Im5GIzciIq44fNi8d82XX5rbQ4fC5MnmfWxExK0o3IiIXMwnn5iXdp88CSEhMGsW3Hef1VWJyHloWkpE5Hzy8mDUKOja1Qw2sbGwc6eCjYib08iNiEhxDh6E+++HrVvN7ZEjYeJECAy0tCwRuTiFGxGRv1q2DAYOhKwsCA2FuXOhRw9raxIRl2laSkTknNxcc4Tm7rvNYNOiBezapWAj4mEUbkREAH78EVq1Mp/gDeZam02boF49S8sSkZLTtJSIyJIl5h2GT52CatVg/nzo0sXqqkTkEmnkRkQqrrNnzfvV9OplBpvWrSEtTcFGxMMp3IhIxbRvH9x6K8yYYW4nJcH69VCrlrV1ichl07SUiFQ8CxfCkCFw5gzUqAFvvQUdO1pdlYiUEo3ciEjFkZNjrq3p29cMNm3bmtNQCjYiXkXhRkQqht27oXlzmDMHbDZIToY1ayAy0urKRKSUaVpKRLzfvHkwfLg5chMebk5L3Xmn1VWJSBnRyI2IeK/TpyEhwbzbcE4OxMWZ01AKNiJeTeFGRLzT11/DLbfAggXg4wPPPQerV0NYmNWViUgZK3G4SUhIYOPGjWVRi4jI5TMMbLNnm49O2LMHoqJgwwb45z/NkCMiXq/Ef9J//fVX4uLiuPbaa3nhhRc4cuRIWdQlIlJi6ftO8fnIXzg+7BnzBn2dO5vTUG3aWF2aiJSjEoebZcuWceTIEYYOHcrixYupV68enTt35r333sNut5dFjSIiF7drF8c7DeTVQwM56lMLJk2CFSugenWrKxORcnZJY7Q1atQgMTGR//znP2zdupUGDRrQr18/IiMjefTRR/nhhx9Ku04RkeIZBkybBrfeiu2/hwEomDkTRo/WNJRIBXVZf/LT09NJSUkhJSUFX19f/va3v/HNN9/QqFEjXn311dKqUUSkeFlZpHcdzM4Rc9iZ15jtNw8EYFfBzezcCTt3Qnq6xTWKSLkrcbix2+28//77dO3albp167J06VIeeeQRjh49yvz581mzZg1LlizhmWeeKYt6RURM27ZB06bM/LQWsewklp0M+fr/AHjoIT9iYyE2FmbOtLhOESl3Jb6JX0REBA6Hg/j4eFJTU2nSpEmRY9q1a0doaGgplCci8heGAVOnwuOPg93OkNor6f78vdC4Mdu25fPQQ37MmJHPLbeYf71FRFhcr4iUuxKHm1dffZV7772XoKCg8x4TGhrKgQMHLqswEZEiTp40b8i3fLm53bMnEbNnE/H7/0zl5xsAxMQYNG1qUY0iYrkSh5t+/fqVRR0iIhe2eTPcfz8cPgwBAfDqqzB0qPmcKBGRP9GlBCLi3hwO87Lu2283g02DBvDVVzBsWJFgEx4OvXrtITzcolpFxC0o3IiI+zp+HLp2hSeegIICiI83L4GKiSn28IgIiI/fq3U2IhWcngouIu5p40YzzBw9CkFB8NprMGiQpqFE5KI0ciMi7qWgwHzIZbt2ZrBp2BBSU2HwYAUbEXGJRm5ExH1kZkLfvrBmjbndv7959+HKla2tS0Q8isKNiLiHdeugTx/IyIDgYDPUDBhgdVUi4oE0LSUi1ioogORkiIszg03jxubdhxVsROQSWR5upk2bRr169QgKCqJFixakpqZe8PisrCyGDx9OREQEgYGBXHfddXz66aflVK2IlKqjR81Q88wz5p2HBw0y19c0amR1ZSLiwSydllq8eDGJiYnMmDGDFi1aMGXKFDp27MjevXupWbNmkePz8vJo3749NWvW5L333iMqKoqff/5Zj3oQ8USrV0O/fubl3pUrmw+B6t3b6qpExAtYGm4mT57M4MGDGTjQfJLvjBkz+OSTT3jzzTd58sknixz/5ptvcvLkSTZv3oy/vz8A9erVK8+SReRy5efDmDHw4ovmdnQ0LFkC111nbV0i4jUsCzd5eXns2LGDpKQk5z4fHx/i4uLYsmVLse9Zvnw5LVu2ZPjw4Xz00UfUqFGD3r1788QTT+Dr61vse3Jzc8nNzXVuZ2dnA+bTze12eyl+I5znK+3zeiP1ynVe1avDh/Ht1w+fzZsBKBgyBMdLL5n3sSmF7+dVvSpj6lXJqF+uK6teleR8loWbEydOUFBQQFhYWKH9YWFh7Nmzp9j3/PTTT6xbt44+ffrw6aefsn//foYNG4bdbic5ObnY90yYMIHx48cX2f/ZZ58RHBx8+V+kGCkpKWVyXm+kXrnO03sVtn07TadOxefUKezBwaQNG8bR1q3Nq6RKmaf3qjypVyWjfrmutHuVk5Pj8rEedSm4w+GgZs2avPHGG/j6+hIbG8uRI0d46aWXzhtukpKSSExMdG5nZ2dTu3ZtOnToQEhISKnWZ7fbSUlJoX379s5pMymeeuU6j+9VXh4+Y8bg++qrADiaNoWFC2lyzTU0KeWP8vhelSP1qmTUL9eVVa/Ozby4wrJwU716dXx9fcnMzCy0PzMzk/DzPPUuIiICf3//QlNQN9xwAxkZGeTl5REQEFDkPYGBgQQGBhbZ7+/vX2a/Qcvy3N5GvXKdR/bq4EHzSd5bt5rbI0fiM3EiPsX8mSxNHtkri6hXJaN+ua60e1WSc1l2KXhAQACxsbGsXbvWuc/hcLB27VpatmxZ7HtatWrF/v37cTgczn379u0jIiKi2GAjIhZatsx8wOXWrRAaCh9+CFOmQBkHGxERS+9zk5iYyKxZs5g/fz7ff/89Q4cO5cyZM86rp/r3719owfHQoUM5efIkI0eOZN++fXzyySe88MILDB8+3KqvICJ/lZsLI0fC3XdDVha0aAG7dkGPHlZXJiIVhKVrbnr16sXx48cZO3YsGRkZNGnShFWrVjkXGR86dAgfnz/yV+3atVm9ejWPPvooN998M1FRUYwcOZInnnjCqq8gIn/244/Qqxfs2GFujxoFL7wAGsYXkXJk+YLiESNGMGLEiGJf27BhQ5F9LVu25KuvvirjqkSkxJYuNe8wnJ0N1arB/PnQpYvVVYlIBWT54xdExMOdPQvDhsF995nBpnVrSEtTsBERyyjciMil27cPbr0Vpk83t5OSYP16qFXL2rpEpEKzfFpKRDzUwoUwZAicOQM1asBbb0HHjlZXJSKikRsRKaGcHHNtTd++ZrBp29achlKwERE3oXAjIq7bvRuaN4c5c8Bmg+RkWLMGIiOtrkxExEnTUiLimnnzYPhwc+QmPNyclrrzTqurEhEpQiM3InJhp09DQgIMHGgGm7g4cxpKwUZE3JTCjYic39dfwy23wIIF4OMDzz0Hq1fD7zfaFBFxR5qWEpGiDANmzTIfo3D2LERFwbvvQps2VlcmInJRCjciUlh2tnmJ96JF5nbnzubITfXq1tYlIuIiTUuJyB927YLYWDPY+PrCpEmwYoWCjYh4FI3ciIg5DfX665CYCHl5UKeOGXBatrS6MhGRElO4EanosrLMm/K9/7653b07zJ0LVataWpaIyKXStJRIRbZtGzRtagYbf3+YMgWWLVOwERGPppEbkYrIMGDqVHj8cbDboX59WLzYvOxbRMTDKdyIVDQnT5o35Fu+3Nzu2RNmz4bQUEvLEhEpLZqWEqlINm+GJk3MYBMQANOmwdKlCjYi4lUUbkQqAofDvKz79tvh8GFo0AC++gqGDTMfgCki4kU0LSXi7Y4fN58NtXKluR0fDzNnwpVXWluXiEgZUbgR8WYbN5ph5uhRCAqC114zL/vWaI2IeDFNS4l4o4IC8yGX7dqZwaZhQ0hNhcGDFWxExOtp5EbE22RmQt++sGaNud2/v7lwuHJla+sSESknCjci3mTdOujTBzIyIDjYDDUDBlhdlYhIudK0lIg3KCiA5GSIizODTePG5t2HFWxEpALSyI2Ipzt61Byt2bDB3B40yLz7cHCwpWWJiFhF4UbEk61eDf36mZd7V65sXuLdu7fVVYmIWErTUiKeKD8fkpKgUycz2ERHw44dCjYiImjkRsTzHD5s3rvmyy/N7aFDYfJk8z42IiKicCPiUT75xLy0++RJCAmBWbPgvvusrkpExK1oWkrEA9jsdnyeeAK6djWDTWws7NypYCMiUgyN3Ii4u4MHaf3Pf+K7b5+5PXIkTJwIgYHW1iUi4qYUbkTc2bJl+A0cSNWsLIzQUGxz50KPHlZXJSLi1jQtJeKOcnPNEZq778aWlcXJ664jPzVVwUZExAUauRFxNz/+CL16mZd2AwWJiXxx2210rlfP2rpERDyERm5E3MnSpdC0qRlsqlWDFStwvPgihp/+P0RExFUKNyLu4OxZGDbMvPopOxtat4a0NOjSxerKREQ8jsKNiNX27YNbb4Xp083tpCRYvx5q1bK2LhERD6WxbhErLVwIQ4bAmTNQowa89RZ07Gh1VSIiHk0jNyJWyMkxn97dt68ZbNq2NaehFGxERC6bW4SbadOmUa9ePYKCgmjRogWpqakuvW/RokXYbDZ66PJY8SS7d0Pz5jBnDthskJwMa9ZAZKTVlYmIeAXLw83ixYtJTEwkOTmZnTt3Eh0dTceOHTl27NgF33fw4EFGjRpFmzZtyqlSkVIwbx7ccgt89x2Eh5uhZtw48PW1ujIREa9hebiZPHkygwcPZuDAgTRq1IgZM2YQHBzMm2++ed73FBQU0KdPH8aPH8/VV19djtWKXKLTpyEhAQYONKek4uLMaag777S6MhERr2PpguK8vDx27NhBUlKSc5+Pjw9xcXFs2bLlvO975plnqFmzJg8++CCbNm264Gfk5uaSm5vr3M7OzgbAbrdjt9sv8xsUdu58pX1eb1ShevX11/j16YNt714MHx8cyck4nngCfHzAhe9foXp1mdQr16lXJaN+ua6selWS81kabk6cOEFBQQFhYWGF9oeFhbFnz55i3/PFF18wZ84c0tLSXPqMCRMmMH78+CL7P/vsM4KDg0tcsytSUlLK5LzeyKt7ZRjU/ewzbpozB1teHr9Vq8b2xERONm4Mq1aV+HRe3atSpl65Tr0qGfXLdaXdq5ycHJeP9ahLwU+dOkW/fv2YNWsW1atXd+k9SUlJJCYmOrezs7OpXbs2HTp0ICQkpFTrs9vtpKSk0L59e/z9/Uv13N7G63uVnY3vsGH4LFkCgKNTJ/zefJNbXfx9+2de36tSpF65Tr0qGfXLdWXVq3MzL66wNNxUr14dX19fMjMzC+3PzMwkPDy8yPE//vgjBw8epFu3bs59DocDAD8/P/bu3cs111xT6D2BgYEEBgYWOZe/v3+Z/QYty3N7G6/s1a5d5p2G9+83FwpPmIDPY4/h43N5S9y8sldlRL1ynXpVMuqX60q7VyU5l6ULigMCAoiNjWXt2rXOfQ6Hg7Vr19KyZcsixzds2JBvvvmGtLQ056/u3bvTrl070tLSqF27dnmWL1KYYcC0aebdhvfvhzp1YNMmGD3aXF8jIiLlwvJpqcTERBISEmjWrBnNmzdnypQpnDlzhoEDBwLQv39/oqKimDBhAkFBQdx4442F3h8aGgpQZL9IucrKMm/K9/775nb37jB3LlStamlZIiIVkeXhplevXhw/fpyxY8eSkZFBkyZNWLVqlXOR8aFDhy57OF+kTG3bBr16wYED4O8PL70EDz9s3qBPRETKneXhBmDEiBGMGDGi2Nc2bNhwwffOmzev9AsScYVhwNSp8Pjj5iXd9evD4sXmTfpERMQybhFuRDzOyZPmDfmWLze3e/aE2bPh92lSERGxjuZ7REpq82Zo0sQMNgEB5iLipUsVbERE3ITCjYirHA6YNAluvx0OH4YGDeCrr2DYMK2vERFxI5qWEnHF8ePms6FWrjS34+Nh5ky48kpr6xIRkSIUbkQuZuNGM8wcPQpBQfDaa+Zl3xqtERFxS5qWEvmT9HQYN878JwUF8Nxz0K6dGWwaNoTUVBg8WMFGRMSNKdyI/El6OowfD+nf/g86dYIxY8y1Nv37m/ezuekmq0sUEZGL0LSUSHHi4+F/ayA42LwaasAAqysSEREXKdxIhZee/sc01M5nVgD/j53/qwNX3wMTJxLR6moirC5SRERcpnAjFd7MmeZUFPgC/w+AwcyGn4B7ITnZXIcjIiKeQeFGKrwh162ne+hzkPULOwNvY3Duv5k1C5o2NV+P0LCNiIhHUbiRiis/H8aMIeLFF81pp+hoSB4NfzeDzblwIyIinkVXS0nFdPgwtG0LL75obg8dat5tuG5dS8sSEZHLp5EbqXg++cS8tPvkSQgJgVmz4L77AHMKKjlZU1EiIp5M4UYqjrw8eOopeOUVczs2FhYvhmuucR4SEaHFwyIink7hRiqGgwfh/vth61Zze+RImDgRAgMtLUtEREqfwo14v2XLYOBAyMqC0FCYOxd69LC2JhERKTNaUCzeKzfXHKG5+24z2LRoAbt2KdiIiHg5hRvxTj/+CK1amU/wBhg1CjZtgnr1LC1LRETKnqalxPssXQqDBkF2NlSrBvPnQ5cuVlclIiLlRCM34j3OnoVhw8zLurOzoXVrSEtTsBERqWAUbsQ77NsHt94K06eb20lJsH491KplbV0iIlLuNC0lnm/hQhgyBM6cgRo14K23oGNHq6sSERGLaORGPFdOjrm2pm9fM9i0bWtOQynYiIhUaAo34pl274bmzWHOHLDZzGcmrFkDkZFWVyYiIhbTtJR4nnnzYPhwc+QmPNyclrrzTqurEhERN6GRG/Ecp09DQoJ5t+GcHIiLM6ehFGxERORPFG7EM3zzDdxyCyxYAD4+8NxzsHo1hIVZXZmIiLgZTUuJezMMmD0bHn7YvI9NVBS8+y60aWN1ZSIi4qYUbsR9ZWebl3gvWmRud+5sjtxUr25tXSIi4tY0LSXuadcuiI01g42vL0yaBCtWKNiIiMhFaeRG3IthwOuvQ2Ii5OVBnTpmwGnZ0urKRETEQyjciPvIyjJvyvf+++Z29+4wdy5UrWppWSIi4lk0LSXuYds2aNrUDDb+/jBlCixbpmAjIiIlpnAj1jIMM8i0agUHDkD9+vDllzBypHnnYRERkRLStJRYxv/UKXx79jQXCgP07Gle9h0aamldIiLi2RRuxBK2r76i7aOP4nPiBAQEwKuvwtChGq0REZHLpmkpKV8OB0yahG+7dgSfOIHRoAF89RUMG6ZgIyIipUIjN1J+jh83nw21ciU24L9t2hC2bBn+WjQsIiKlyC1GbqZNm0a9evUICgqiRYsWpKamnvfYWbNm0aZNG6666iquuuoq4uLiLni8uImNG6FJE1i5EoKCyJ8+nR2JiXDllVZXJiIiXsbycLN48WISExNJTk5m586dREdH07FjR44dO1bs8Rs2bCA+Pp7169ezZcsWateuTYcOHThy5Eg5Vy4ucTjg+eehXTs4ehQaNoTUVIwHH9Q0lIiIlAnLw83kyZMZPHgwAwcOpFGjRsyYMYPg4GDefPPNYo9fuHAhw4YNo0mTJjRs2JDZs2fjcDhYu3ZtOVcuF5WZCZ06wdNPmyGnf3/zfjY33WR1ZSIi4sUsDTd5eXns2LGDuLg45z4fHx/i4uLYsmWLS+fIycnBbrdTVes23Mu6deY0VEoKBAfDvHkwfz5Urmx1ZSIi4uUsXVB84sQJCgoKCAsLK7Q/LCyMPXv2uHSOJ554gsjIyEIB6c9yc3PJzc11bmdnZwNgt9ux2+2XWHnxzp2vtM/rUQoK8HnuOXxeeAGbYWA0akT+O+9Ao0bwp76oV65Tr1ynXrlOvSoZ9ct1ZdWrkpzPo6+WevHFF1m0aBEbNmwgKCio2GMmTJjA+PHji+z/7LPPCA4OLpO6UlJSyuS87i7o5EmaTp5MjW+/BeBg+/Z8O2gQBQcPwsGDxb6novbqUqhXrlOvXKdelYz65brS7lVOTo7Lx1oabqpXr46vry+ZmZmF9mdmZhIeHn7B97788su8+OKLrFmzhptvvvm8xyUlJZGYmOjczs7Odi5CDgkJubwv8Bd2u52UlBTat2+Pv79/qZ7b3dlSUvB94glsx49jVK5MwbRpRMXHE3We4ytyr0pKvXKdeuU69apk1C/XlVWvzs28uMLScBMQEEBsbCxr166lR48eAM7FwSNGjDjv+yZNmsTzzz/P6tWradas2QU/IzAwkMDAwCL7/f39y+w3aFme2+3k58PYsTBhgrkdHY1tyRL8rrvOpbdXqF5dJvXKdeqV69SrklG/XFfavSrJuSyflkpMTCQhIYFmzZrRvHlzpkyZwpkzZxg4cCAA/fv3Jyoqigm///CcOHEiY8eO5Z133qFevXpkZGQAULlyZSprsWr5OnwY4uPNB12CeZfhV16B80wRioiIlAfLw02vXr04fvw4Y8eOJSMjgyZNmrBq1SrnIuNDhw7h4/PHRV3Tp08nLy+Pe+65p9B5kpOTGTduXHmWXrF98ol5affJkxASYj7w8t57ra5KRETE+nADMGLEiPNOQ23YsKHQ9sHzLEyVcmK3w1NPwcsvm9uxsbB4MVxzjbV1iYiI/M4two14iIMH4f77YetWc3vkSJg4EYpZ0yQiImIVhRtxzbJlMHAgZGVBaCjMnQu/LwIXERFxJ5Y/fkHcXG6uOUJz991msGnRAtLSFGxERMRtKdzI+f34I7RqBa+9Zm6PGgWbNkHdutbWJSIicgGalpLiLV0KgwZBdjZUq2Y+F6pLF6urEhERuSiN3EhhZ8+a96u57z4z2LRubU5DKdiIiIiHULiRP+zbB7feCtOng81mXvK9fj3UqmV1ZSIiIi7TtJSYFi6EIUPgzBmoUQPefhs6dLC6KhERkRLTyE1Fl5Njrq3p29cMNm3bmtNQCjYiIuKhFG4qst27oXlzmDPHnIZKToY1ayAy0urKRERELpmmpSqqefNg+HBz5CY83JyWuvNOq6sSERG5bBq5qWhOn4aEBPNuwzk50L69OQ2lYCMiIl5C4aYi+eYbuOUWWLAAfHzguedg1Sr4/QnsIiIi3kDTUhWBYcDs2fDww+Z9bKKi4N13oU0bqysTEREpdQo33i4727zEe9Eic7tzZ3Pkpnp1a+sSEREpI5qW8ma7dkFsrBls/Pxg0iRYsULBRkREvJpGbryRYZh3GX70UcjLgzp1zIDTsqXVlYmIiJQ5hRtvk5UFgwfDe++Z2927w9y5ULWqpWWJiIiUF01LeZNt26BpUzPY+PvDlCmwbJmCjYiIVCgaufEGhgFTp8Ljj4PdDvXrw+LF5mXfIiIiFYzCjac7edK8Id/y5eZ2z57mZd+hoZaWJSIiYhVNS3myLVsgJsYMNgEBMG0aLF2qYCMiIhWawo0ncjjMy7rbtIFDh6BBA/jqKxg2zHwApoiISAWmaSlPc+IE9O8PK1ea2/HxMHMmXHmltXWJiIi4CYUbT7JpkxlmjhyBoCD417/gwQc1WiMiIvInmpbyBA4HPP88tG1rBpuGDSE1FQYNUrARERH5C43cuLvMTOjXD1JSzO3+/c2Fw5UrW1uXiIiIm1K4cWfr1kGfPpCRAcHB8PrrkJBgdVUiIiJuTdNS7qigAJKTIS7ODDY33mjefVjBRkRE5KI0cuNujh41R2s2bDC3Bw0y7z4cHGxpWSIiIp5C4cadfPYZ9O0Lx4+ba2pmzoTeva2uSkRExKNoWsod5OfDU09Bx45msImOhh07FGxEREQugUZurPbf/5r3rvniC3N72DB45RXzPjYiIiJSYgo3VvrkE3OR8P/+ByEh5gMv773X6qpEREQ8mqalrGC3w+jR0LWrGWxiY2HnTgUbERGRUqCRm/J28CDcfz9s3WpujxwJEydCYKClZYmIiHgLhZvytGwZDBwIWVkQGgpz50KPHtbWJCIi4mU0LVUecnPhkUfg7rvNYHPrrZCWpmAjIiJSBhRuytqPP0KrVuaN+ABGjYKNG6FuXWvrEhER8VKalipLS5eadxjOzoZq1WD+fOjSxeqqREREvJpbjNxMmzaNevXqERQURIsWLUhNTb3g8UuXLqVhw4YEBQVx00038emnn5ZTpReWng7vvns96QdzzfvV3HefGWxatzanoRRsREREypzl4Wbx4sUkJiaSnJzMzp07iY6OpmPHjhw7dqzY4zdv3kx8fDwPPvggu3btokePHvTo0YNvv/22nCsvKiMDFi9uyPG/PQjTp4PNZt55eP16qFXL6vJEREQqBMvDzeTJkxk8eDADBw6kUaNGzJgxg+DgYN58881ij586dSqdOnVi9OjR3HDDDTz77LM0bdqUf//73+VceTFWrQbAtv8HqFEDVq2C558HP83+iYiIlBdLf+rm5eWxY8cOkpKSnPt8fHyIi4tjy5Ytxb5ny5YtJCYmFtrXsWNHli1bVuzxubm55ObmOrezs7MBsNvt2O32y/wG5lRURgbw8Qq+fu4ToCvbr72f/Ff7QWh1wg/ZiYi47I/xOud6Xxr/DbydeuU69cp16lXJqF+uK6teleR8loabEydOUFBQQFhYWKH9YWFh7Nmzp9j3ZGRkFHt8RkZGscdPmDCB8ePHF9n/2WefERwcfImV/+Hdd69n8eKGwN2//4IhPzwOXc3Xe/XaQ3z83sv+HG+VkpJidQkeQ71ynXrlOvWqZNQv15V2r3Jyclw+1uvnS5KSkgqN9GRnZ1O7dm06dOhASEjIZZ8/JgYSE800uX1LHsMfuYJp03Jp1syc8QsPv4aIiGsu+3O8jd1uJyUlhfbt2+Pv7291OW5NvXKdeuU69apk1C/XlVWvzs28uMLScFO9enV8fX3JzMwstD8zM5Pw8PBi3xMeHl6i4wMDAwks5tEG/v7+pdL0OnXMX3/WrJkPzZvrN78rSuu/Q0WgXrlOvXKdelUy6pfrSrtXJTmXpQuKAwICiI2NZe3atc59DoeDtWvX0rJly2Lf07Jly0LHgzn0db7jRUREpGKxfFoqMTGRhIQEmjVrRvPmzZkyZQpnzpxh4MCBAPTv35+oqCgmTJgAwMiRI7njjjt45ZVX6NKlC4sWLWL79u288cYbVn4NAMLDzTU24eGahhIREbGK5eGmV69eHD9+nLFjx5KRkUGTJk1YtWqVc9HwoUOH8PH5Y4Dptttu45133uHpp5/mqaee4tprr2XZsmXceOONVn0Fp4gIiI/fqzU2IiIiFrI83ACMGDGCESNGFPvahg0biuy79957uffee8u4KhEREfFElt/ET0RERKQ0KdyIiIiIV1G4EREREa+icCMiIiJeReFGREREvIrCjYiIiHgVhRsRERHxKgo3IiIi4lUUbkRERMSruMUdisuTYRhAyR6d7iq73U5OTg7Z2dl6auxFqFeuU69cp165Tr0qGfXLdWXVq3M/t8/9HL+QChduTp06BUDt2rUtrkRERERK6tSpU1SpUuWCx9gMVyKQF3E4HBw9epQrr7wSm81WqufOzs6mdu3aHD58mJCQkFI9t7dRr1ynXrlOvXKdelUy6pfryqpXhmFw6tQpIiMjCz1QuzgVbuTGx8eHWrVqlelnhISE6De/i9Qr16lXrlOvXKdelYz65bqy6NXFRmzO0YJiERER8SoKNyIiIuJVFG5KUWBgIMnJyQQGBlpdittTr1ynXrlOvXKdelUy6pfr3KFXFW5BsYiIiHg3jdyIiIiIV1G4EREREa+icCMiIiJeReFGREREvIrCTSmZNm0a9erVIygoiBYtWpCammp1SZabMGECt9xyC1deeSU1a9akR48e7N27t9AxZ8+eZfjw4VSrVo3KlSvTs2dPMjMzLarYfbz44ovYbDYeeeQR5z716g9Hjhyhb9++VKtWjUqVKnHTTTexfft25+uGYTB27FgiIiKoVKkScXFx/PDDDxZWbJ2CggLGjBlD/fr1qVSpEtdccw3PPvtsoefzVNR+bdy4kW7duhEZGYnNZmPZsmWFXnelLydPnqRPnz6EhIQQGhrKgw8+yOnTp8vxW5SPC/XKbrfzxBNPcNNNN3HFFVcQGRlJ//79OXr0aKFzlGevFG5KweLFi0lMTCQ5OZmdO3cSHR1Nx44dOXbsmNWlWerzzz9n+PDhfPXVV6SkpGC32+nQoQNnzpxxHvPoo4/y8ccfs3TpUj7//HOOHj3K3//+dwurtt62bduYOXMmN998c6H96pXpl19+oVWrVvj7+7Ny5Up2797NK6+8wlVXXeU8ZtKkSbz22mvMmDGDrVu3csUVV9CxY0fOnj1rYeXWmDhxItOnT+ff//4333//PRMnTmTSpEn861//ch5TUft15swZoqOjmTZtWrGvu9KXPn368N1335GSksKKFSvYuHEj//jHP8rrK5SbC/UqJyeHnTt3MmbMGHbu3MkHH3zA3r176d69e6HjyrVXhly25s2bG8OHD3duFxQUGJGRkcaECRMsrMr9HDt2zACMzz//3DAMw8jKyjL8/f2NpUuXOo/5/vvvDcDYsmWLVWVa6tSpU8a1115rpKSkGHfccYcxcuRIwzDUqz974oknjNatW5/3dYfDYYSHhxsvvfSSc19WVpYRGBhovPvuu+VRolvp0qWL8cADDxTa9/e//93o06ePYRjq1zmA8eGHHzq3XenL7t27DcDYtm2b85iVK1caNpvNOHLkSLnVXt7+2qvipKamGoDx888/G4ZR/r3SyM1lysvLY8eOHcTFxTn3+fj4EBcXx5YtWyyszP38+uuvAFStWhWAHTt2YLfbC/WuYcOG1KlTp8L2bvjw4XTp0qVQT0C9+rPly5fTrFkz7r33XmrWrElMTAyzZs1yvn7gwAEyMjIK9apKlSq0aNGiwvUK4LbbbmPt2rXs27cPgP/85z988cUXdO7cGVC/zseVvmzZsoXQ0FCaNWvmPCYuLg4fHx+2bt1a7jW7k19//RWbzUZoaChQ/r2qcA/OLG0nTpygoKCAsLCwQvvDwsLYs2ePRVW5H4fDwSOPPEKrVq248cYbAcjIyCAgIMD5m/+csLAwMjIyLKjSWosWLWLnzp1s27atyGvq1R9++uknpk+fTmJiIk899RTbtm3j4YcfJiAggISEBGc/ivszWdF6BfDkk0+SnZ1Nw4YN8fX1paCggOeff54+ffoAqF/n4UpfMjIyqFmzZqHX/fz8qFq1aoXu3dmzZ3niiSeIj493PjizvHulcCPlYvjw4Xz77bd88cUXVpfilg4fPszIkSNJSUkhKCjI6nLcmsPhoFmzZrzwwgsAxMTE8O233zJjxgwSEhIsrs79LFmyhIULF/LOO+/QuHFj0tLSeOSRR4iMjFS/pNTZ7Xbuu+8+DMNg+vTpltWhaanLVL16dXx9fYtctZKZmUl4eLhFVbmXESNGsGLFCtavX0+tWrWc+8PDw8nLyyMrK6vQ8RWxdzt27ODYsWM0bdoUPz8//Pz8+Pzzz3nttdfw8/MjLCxMvfpdREQEjRo1KrTvhhtu4NChQwDOfujPpGn06NE8+eST3H///dx0003069ePRx99lAkTJgDq1/m40pfw8PAiF47k5+dz8uTJCtm7c8Hm559/JiUlxTlqA+XfK4WbyxQQEEBsbCxr16517nM4HKxdu5aWLVtaWJn1DMNgxIgRfPjhh6xbt4769esXej02NhZ/f/9Cvdu7dy+HDh2qcL276667+Oabb0hLS3P+atasGX369HH+u3platWqVZFbCuzbt4+6desCUL9+fcLDwwv1Kjs7m61bt1a4XoF5JYuPT+G/6n19fXE4HID6dT6u9KVly5ZkZWWxY8cO5zHr1q3D4XDQokWLcq/ZSueCzQ8//MCaNWuoVq1aodfLvVelvkS5Alq0aJERGBhozJs3z9i9e7fxj3/8wwgNDTUyMjKsLs1SQ4cONapUqWJs2LDBSE9Pd/7KyclxHvPQQw8ZderUMdatW2ds377daNmypdGyZUsLq3Yff75ayjDUq3NSU1MNPz8/4/nnnzd++OEHY+HChUZwcLDx9ttvO4958cUXjdDQUOOjjz4yvv76a+P//b//Z9SvX9/47bffLKzcGgkJCUZUVJSxYsUK48CBA8YHH3xgVK9e3Xj88cedx1TUfp06dcrYtWuXsWvXLgMwJk+ebOzatct5hY8rfenUqZMRExNjbN261fjiiy+Ma6+91oiPj7fqK5WZC/UqLy/P6N69u1GrVi0jLS2t0N/3ubm5znOUZ68UbkrJv/71L6NOnTpGQECA0bx5c+Orr76yuiTLAcX+mjt3rvOY3377zRg2bJhx1VVXGcHBwcbdd99tpKenW1e0G/lruFGv/vDxxx8bN954oxEYGGg0bNjQeOONNwq97nA4jDFjxhhhYWFGYGCgcddddxl79+61qFprZWdnGyNHjjTq1KljBAUFGVdffbXxz3/+s9APnYrar/Xr1xf7d1RCQoJhGK715X//+58RHx9vVK5c2QgJCTEGDhxonDp1yoJvU7Yu1KsDBw6c9+/79evXO89Rnr2yGcafblMpIiIi4uG05kZERES8isKNiIiIeBWFGxEREfEqCjciIiLiVRRuRERExKso3IiIiIhXUbgRERERr6JwIyIiIl5F4UZE3ELbtm155JFHrC5DRLyAwo2IiIh4FYUbEbHcgAED+Pzzz5k6dSo2mw2bzcbBgweLPTY3N5dRo0YRFRXFFVdcQYsWLdiwYYPz9Xnz5hEaGsrq1au54YYbqFy5Mp06dSI9Pf28n//LL7/Qp08fatSoQaVKlbj22muZO3duKX9LESkvflYXICIydepU9u3bx4033sgzzzwDQI0aNYo9dsSIEezevZtFixYRGRnJhx9+SKdOnfjmm2+49tprAcjJyeHll1/mrbfewsfHh759+zJq1CgWLlxY7DnHjBnD7t27WblyJdWrV2f//v389ttvZfNlRaTMKdyIiOWqVKlCQEAAwcHBhIeHn/e4Q4cOMXfuXA4dOkRkZCQAo0aNYtWqVcydO5cXXngBALvdzowZM7jmmmsAMxCdC03nO29MTAzNmjUDoF69eqX0zUTECgo3IuIxvvnmGwoKCrjuuusK7c/NzaVatWrO7eDgYGewAYiIiODYsWPnPe/QoUPp2bMnO3fupEOHDvTo0YPbbrut9L+AiJQLhRsR8RinT5/G19eXHTt24OvrW+i1ypUrO//d39+/0Gs2mw3DMM573s6dO/Pzzz/z6aefkpKSwl133cXw4cN5+eWXS/cLiEi5ULgREbcQEBBAQUHBBY+JiYmhoKCAY8eO0aZNm1L9/Bo1apCQkEBCQgJt2rRh9OjRCjciHkrhRkTcQr169di6dSsHDx6kcuXKVK1aFR+fwhd0XnfddfTp04f+/fvzyiuvEBMTw/Hjx1m7di0333wzXbp0uaTPHjt2LLGxsTRu3Jjc3FxWrFjBDTfcUBpfS0QsoEvBRcQtjBo1Cl9fXxo1akSNGjU4dOhQscfNnTuX/v3789hjj3H99dfTo0cPtm3bRp06dS75swMCAkhKSuLmm2/m9ttvx9fXl0WLFl3y+UTEWjbjQhPRIiIiIh5GIzciIiLiVRRuRERExKso3IiIiIhXUbgRERERr6JwIyIiIl5F4UZERES8isKNiIiIeBWFGxEREfEqCjciIiLiVRRuRERExKso3IiIiIhXUbgRERERr/L/AYFPW79NGxvbAAAAAElFTkSuQmCC\n"},"metadata":{}},{"output_type":"stream","name":"stdout","text":["pente =  1.11e-02 /s   y_0 =  0.00\n","k_moy = 1.1168e-02  avec u(k_moy) =  6.4e-05 par calcul direct\n"]}]},{"cell_type":"markdown","metadata":{"id":"lPoC5sRmSa9h"},"source":["activité 7.3 Décomposition de l’éthylamine en phase gazeuse\n","\n","1°)"]},{"cell_type":"code","metadata":{"colab":{"base_uri":"https://localhost:8080/","height":508},"id":"xtPRfY8tSjEn","executionInfo":{"status":"ok","timestamp":1695277822632,"user_tz":-120,"elapsed":683,"user":{"displayName":"Lionnel Malara","userId":"16342666118884215780"}},"outputId":"6b74f9a6-9d56-48b6-f86a-2e263ffe0dba"},"source":["import matplotlib.pyplot as plt #Pour tracer des graphiques\n","import numpy as np #Pour faire divers calculs\n","\n","## Données\n","\n","t=np.array([0, 1,4,10]) # valeur du temps en min\n","P_T=np.array([7240,7900,9480,11720]) # valeur de la pression totale en Pa\n","P_0=7240  # valeur de la pression initiale en Pa\n","PE=2*P_0-P_T # valeur de la pression partielle de l'éthyleamine en Pa\n","y =np.log(PE) # valeur de la fonction\n","\n","# Courbe d'évolution de la fonction y=f(t) et régression linéaire\n","\n","plt.figure(1)\n","p=np.polyfit(t,y,1)                     #Régression linéaire de y en fonction de t : équation de la forme y = p[0]t + p[1]\n","plt.plot(t,np.polyval(p,t),'r')         #Tracé de la droite de régression en rouge ; polyval permet de placer les ordonnées y prévues par la droite modèle.\n","plt.xlabel('t en min') #Légende de l’axe des abscisses\n","plt.ylabel('y') #Légende de l’axe des ordonnées\n","plt.plot(t,y,'b+',label='y=f(t)') #Représente y en fonction de x avec des croix bleues avec une étiquette pour la courbe\n","plt.title('vérification ordre 1') #Titre du graphe\n","plt.legend() #Affiche l’étiquette de la courbe\n","plt.grid() #Affiche le quadrillage\n","plt.show() #Affiche le graphique\n","\n","print('pente = ',format(p[0],\"#.2e\"),'/min ',' y_0 = ',format(p[1],\"#.2f\"))\t#Affiche le résultat\n","\n","t2=t[1:]\n","y2=-y[1:]+np.log(P_0)\n","k2=y2/t2\n","k_moy=np.mean(k2) #moyenne de k\n","u_kmoy= np.std(k2, ddof=1)/3**(1/2) #incertitude-type de k_moy\n","print (\"k_moy =\",format(k_moy,\"#.2e\"),\" avec u(k_moy) = \",format(u_kmoy,\"#.1e\"),'par calcul direct') #Affiche résultat et incertitude-type\n","\n"],"execution_count":5,"outputs":[{"output_type":"display_data","data":{"text/plain":["<Figure size 640x480 with 1 Axes>"],"image/png":"iVBORw0KGgoAAAANSUhEUgAAAjcAAAHHCAYAAABDUnkqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABTpklEQVR4nO3deZyNdf/H8deZM4sZW8kyY52xL8mMNVvIGolKCTeDuBUqRKFbyJbdjSIqVNZEKZXmln0b20hlabHFMGQZTDTmXL8/vj9TkxkGM+eaOfN+Ph7zuOdcc13X+ZxvbvN2fTeHZVkWIiIiIh7Cy+4CRERERNKSwo2IiIh4FIUbERER8SgKNyIiIuJRFG5ERETEoyjciIiIiEdRuBERERGPonAjIiIiHkXhRsTDxMfH8+abb/L555/bXYqIiC0UbkQ8zMCBA3n33Xd58MEHU33N3LlzcTgcHD58OMnx8ePHU7x4cZxOJ6GhoQAEBwfTuXPntCs4lYYNG4bD4XD7+6anlNpdRO6Owo2IB/nss8/46KOP+Prrr8mXL99d3eubb77hlVdeoXbt2syZM4fRo0enUZUpi4uLY9iwYaxduzbd3yuriYyMpGfPnlSpUgUfHx+PC4oif+dtdwEiknYOHz7MV199RcmSJW/ruo4dO/LMM8/g5+eXeOzbb7/Fy8uL9957D19f38TjBw4cwMsrff5dFBcXx/DhwwGoX79+kp/95z//YeDAgenyvlnBl19+ybvvvssDDzxA8eLFOXjwoN0liaQbPbkR8SAvvfQSlStXTvX5ly9fBsDpdJItW7Yk/5qPiYnB398/SbAB8PPzw8fHJ20Kvg3e3t5ky5bN7e97N1wuF1euXLmre1iWxR9//HHXtTz//PNcuHCBHTt20Lhx47u+n0hGpnAjkgEtXboUh8PBunXrbvjZO++8g8Ph4Pvvv088tn//ftq0aUOePHnIli0bVatWZcWKFUmuuz6+Y926dfTs2ZP8+fNTuHDhJD+7PvbD4XAwZ84cLl++jMPhwOFwMHfuXCD5MTfnz5+nb9++BAcH4+fnR+HChenUqRNnzpwB4M8//+T111+nSpUq5M6dm+zZs1O3bl3WrFmTeI/Dhw8ndqUNHz488X2HDRsGJD/m5tq1a4wYMYISJUrg5+dHcHAwgwcP5urVq0nOCw4O5tFHH2Xjxo1Ur16dbNmyUbx4cT744INU/NcwIfDll1+mSJEi+Pn5UaZMGSZMmIBlWUnOczgc9O7dm/nz51OhQgX8/Pz4+uuvAfjhhx94+OGH8ff3p3DhwowcORKXy3XDe12vddWqVVStWhV/f3/eeeedxHbu06dPYh0lS5Zk7Nixyd7nnwoUKIC/v3+qPq9IZqduKZEMqEWLFuTIkYMlS5ZQr169JD9bvHgxFSpU4P777wfML83atWtTqFAhBg4cSPbs2VmyZAmtW7fmk08+4fHHH09yfc+ePcmXLx+vv/564pObf/rwww+ZNWsWkZGRvPvuuwDUqlUr2XMvXbpE3bp12bdvH127dqVy5cqcOXOGFStW8Ntvv5E3b15iY2N59913adeuHd27d+fixYu89957NG3alMjISEJDQ8mXLx8zZszg+eef5/HHH+eJJ54A4IEHHkixnbp168a8efNo06YNL7/8Mtu2bWPMmDHs27eP5cuXJzn3559/pk2bNjz77LOEh4fz/vvv07lzZ6pUqUKFChVSfA/LsnjsscdYs2YNzz77LKGhoaxatYoBAwZw/PhxJk+enOT8b7/9liVLltC7d2/y5s1LcHAwJ0+epEGDBly7di3xv9GsWbNSDBsHDhygXbt29OjRg+7du1OmTBni4uKoV68ex48fp0ePHhQtWpTNmzczaNAgoqOjmTJlSoqfQSTLsUQkQ2rXrp2VP39+69q1a4nHoqOjLS8vL+uNN95IPNawYUOrYsWK1pUrVxKPuVwuq1atWlapUqUSj82ZM8cCrDp16iS5599/dujQocRj4eHhVvbs2W+oq1ixYlZ4eHji69dff90CrGXLlt1wrsvlsizLsq5du2ZdvXo1yc/OnTtnFShQwOratWvisdOnT1uANXTo0BvuNXToUOvvf2VFRUVZgNWtW7ck5/Xv398CrG+//TZJzYC1fv36xGMxMTGWn5+f9fLLL9/wXn/36aefWoA1cuTIJMfbtGljORwO6+eff048BlheXl7WDz/8kOTcPn36WIC1bdu2JO+fO3fuG9r9eq1ff/11knuMGDHCyp49u3Xw4MEkxwcOHGg5nU7r6NGjN/0cf9erVy9Lf/2LJ1O3lEgG1bZtW2JiYpLMHFq6dCkul4u2bdsCcPbsWb799luefvppLl68yJkzZzhz5gy///47TZs25aeffuL48eNJ7tu9e3ecTmea1fnJJ59QqVKlG54QAYndSE6nM3Hsjsvl4uzZs1y7do2qVauya9euO3rfL7/8EoB+/folOf7yyy8DsHLlyiTHy5cvT926dRNf58uXjzJlyvDrr7/e8n2cTicvvvjiDe9jWRZfffVVkuP16tWjfPnyN9zjwQcfpHr16knev0OHDsm+Z0hICE2bNk1y7OOPP6Zu3brce++9if+dz5w5Q6NGjUhISGD9+vU3/RwiWYm6pUQyqGbNmpE7d24WL15Mw4YNAdMlFRoaSunSpQHT1WJZFkOGDGHIkCHJ3icmJoZChQolvg4JCUnTOn/55ReefPLJW543b948Jk6cyP79+4mPj7/reo4cOYKXl9cNM8MCAwO55557OHLkSJLjRYsWveEe9957L+fOnbvl+xQsWJCcOXMmOV6uXLnEn/9dcp/nyJEj1KhR44bjZcqUSfY9k7vHTz/9xHfffZfiFP+YmJjkP4BIFqRwI5JB+fn50bp1a5YvX87bb7/NqVOn2LRpU5L1Zq4PJO3fv/8N/9K/7p+//O0YVPrRRx/RuXNnWrduzYABA8ifPz9Op5MxY8bwyy+/3NW9U7teS0pPq6x/DAq+W2nRvsndw+Vy0bhxY1555ZVkr7keeEVE4UYkQ2vbti3z5s1j9erV7Nu3D8uyErukAIoXLw6Aj48PjRo1sqXGEiVKJJm5lZylS5dSvHhxli1bliSMDB06NMl5t7OwXLFixXC5XPz000+JT1EATp06xfnz5ylWrFiq73Wr9/nf//7HxYsXkzy92b9/f+LPU3OPn3766YbjBw4cSHUdJUqU4NKlS7b9dxbJTDTmRiQDa9SoEXny5GHx4sUsXryY6tWrJ+myyJ8/P/Xr1+edd94hOjr6hutPnz6d7jU++eST7Nmz54bZSfDXU5HrT03+/pRk27ZtbNmyJcn5AQEBgJnyfCvNmzcHuGGW0KRJkwAz4ywtNG/enISEBKZPn57k+OTJk3E4HDzyyCOpusfWrVuJjIxMPHb69Gnmz5+f6jqefvpptmzZwqpVq2742fnz57l27Vqq7yXi6fTkRiQD8/Hx4YknnmDRokVcvnyZCRMm3HDOW2+9RZ06dahYsSLdu3enePHinDp1ii1btvDbb7+xZ8+edK1xwIABLF26lKeeeoquXbtSpUoVzp49y4oVK5g5cyaVKlXi0UcfZdmyZTz++OO0aNGCQ4cOMXPmTMqXL8+lS5cS7+Xv70/58uVZvHgxpUuXJk+ePNx///2J097/rlKlSoSHhzNr1izOnz9PvXr1iIyMZN68ebRu3ZoGDRqkyedr2bIlDRo04LXXXuPw4cNUqlSJb775hs8++4w+ffpQokSJW97jlVde4cMPP6RZs2a89NJLiVPBixUrxnfffZeqOgYMGMCKFSt49NFHE6ewX758mb1797J06VIOHz5M3rx5U7z+yJEjfPjhhwDs2LEDgJEjRwLmyVLHjh1TVYdIpmDnVC0RubWIiAgLsBwOh3Xs2LFkz/nll1+sTp06WYGBgZaPj49VqFAh69FHH7WWLl2aeM716d7bt2+/4fq7mQpuWZb1+++/W71797YKFSpk+fr6WoULF7bCw8OtM2fOWJZlpoSPHj3aKlasmOXn52eFhYVZX3zxhRUeHm4VK1Ysyb02b95sValSxfL19U0yLfyfU8Ety7Li4+Ot4cOHWyEhIZaPj49VpEgRa9CgQUmmxV+vuUWLFjd8lnr16ln16tW74fg/Xbx40erbt69VsGBBy8fHxypVqpQ1fvz4xKnu1wFWr169kr3Hd999Z9WrV8/Kli2bVahQIWvEiBHWe++9l+xU8ORqvV7HoEGDrJIlS1q+vr5W3rx5rVq1alkTJkyw/vzzz5t+hjVr1lhAsl+paQORzMRhWWk8mk5ERETERhpzIyIiIh5F4UZEREQ8isKNiIiIeBSFGxEREfEoCjciIiLiURRuRERExKNkuUX8XC4XJ06cIGfOnLe11LuIiIjYx7IsLl68SMGCBfHyuvmzmSwXbk6cOEGRIkXsLkNERETuwLFjxyhcuPBNz8ly4eb6xnfHjh0jV65caXrv+Ph4vvnmG5o0aYKPj0+a3lv+onZ2D7Wze6id3Udt7R7p1c6xsbEUKVIkyQa2Kcly4eZ6V1SuXLnSJdwEBASQK1cu/R8nHamd3UPt7B5qZ/dRW7tHerdzaoaUaECxiIiIeBSFGxEREfEoCjciIiLiUbLcmBsREZE7lZCQQHx8vN1lZGjx8fF4e3tz5coVEhISbutaX1/fW07zTg2FGxERkVuwLIuTJ09y/vx5u0vJ8CzLIjAwkGPHjt32enJeXl6EhITg6+t7VzUo3IiIiNzC9WCTP39+AgICtAjsTbhcLi5dukSOHDlu6ynM9UV2o6OjKVq06F21scKNiIjITSQkJCQGm/vuu8/ucjI8l8vFn3/+SbZs2W67iylfvnycOHGCa9eu3dU0cg0oFhERuYnrY2wCAgJsrsTzXe+Out2xOv+kcCMiIpIK6opKf2nVxgo3IiIi4lEUbkREROSW4uLiePLJJ8mVKxcOhyNx5tiBAwcIDAzk4sWLN73+wQcf5JNPPnFDpQo3aSo6GhYuLEN0tN2ViIiIpK158+axYcMGNm/eTHR0NLlz5wZg0KBBvPDCC4kbWs6dO5dixYrdcP1//vMfBg4ciMvlSvdaFW7S0Mlv9rJ4cVlOnrS7EhERyaiio2HYMDLdP4R/+eUXypUrx/33309gYCAOh4OjR4/yxRdf0Llz51te/8gjj3Dx4kW++uqrdK9V4SatHD2Ks18/8/3uKFtLERGRjCs6GoYPT/9w88EHH3Dfffdx9erVJMdbt25Nx44db+te9evXZ+LEiaxfvx6Hw0H9+vUBWLJkCZUqVaJQoUIArF27lmeffZbY2FicTicOh4Nhw4YB4HQ6ad68OYsWLbrrz3YrCjd3KToadu2CXdsT2BVQG4A9z89mV6/32LU9IdMlcxER8QxPPfUUCQkJrFixIvFYTEwMK1eupGvXrmzYsIEcOXLc9Gv+/PkALFu2jO7du1OzZk2io6NZtmwZABs2bKBq1aqJ969VqxaTJ08mZ86cHD9+nOjoaPr375/48+rVq7Nhw4Z0/+xaxO8uvfOOSeAQAowC4N/WLHgbeBuGvnyJYRNy2FihiIjYLTr6ryc1u3Yl/V+AoCDzlZb8/f1p3749c+bM4amnngLgo48+omjRotSvX58rV64QFRV103sUKFAAgDx58hAQEICvry+BgYGJPz9y5EiScOPr65s44DgwMPCGRfwKFizIsWPHcLlcabKHVEoUbu5Sjx7w2GPm++3br/Hcc97M7LSBaksGwpU/CJp3DZpNgkaN7C1URERs89c/hP/Svftf3w8dasbhpLXu3btTrVo1jh8/TqFChZg7dy6dO3fG4XDg7+9PyZIl7+r+f/zxB9myZUv1+f7+/rhcLq5evYq/v/9dvffNqFvqLgUFQeXK5isszAIgrNeDVN71LpUrXiPozF5o0gT+8x+4ds3makVExA49esDOneZr9mxzbPbsv4716JE+7xsWFkalSpX44IMP2LlzJz/88EPi4N/b6ZZKSd68eTl37lyq6zl79izZs2dP12ADenKTfsqVg23boG9fE9lHjYJ162DhQihc2O7qRETEjZLrdrr+D+P01q1bN6ZMmcLx48dp1KgRRYoUAaBq1aqp7pZKSVhYGD/++GOSY76+vilO9/7+++8JCwtLffF3SE9u0lBgILRtu5/E7kh/f5g5ExYtgpw5YeNGCA2FlSvtLFNERLKQ9u3b89tvvzF79my6du2aePx6t9TNvq6vXZOSpk2bsmXLliR7QQUHB3Pp0iVWr17NmTNniIuLS/zZhg0baNKkSdp/yH9QuElDQUHQrt2BGweFtW0Lu3dDlSrw++/w6KPQvz/8+actdYqIiH2CgswYm7QeQJyS3Llz8+STT5IjRw5at26dpvd+5JFH8Pb25n//+1/isVq1atGlSxfatWtHvnz5GDduHADHjx9n8+bNdOnSJU1rSI7CjbuUKAGbNkGfPub1xIlQty4cOmRrWSIi4l5BQWbwsLvCDZhg0aFDB/z8/O74HlOmTGHt2rVJjnl7ezN48GAmTZqU5PikSZOIiYnBsqzEdW6mTp1K586dKeyGoRkKN+7k5weTJ8Onn8K990JkJISFgZv22hARkazl3LlzLF++nLVr19KrV690eY8ePXrw0EMP3XJvqfz58zNixIh0qeGfFG7s0KoVREVBzZpw4QK0aQO9e8OVK3ZXJiIiHiQsLIzOnTszduxYypQpky7v4e3tzWuvvXbL8Tkvv/zyLQcopxWFG7sULWpmTw0caF6/9ZYJOwcP2luXiIh4jMOHD3PhwoUkqwRnBQo3dvLxgTFj4KuvIF8+8zSnShVYsMDuykRERDIthZuMoFkzE2zq14dLl6BDB+jWDf42fU5EROxlWZbdJXi8tGpjhZuMomBB+N//zBB6hwPeew+qVYMffrC7MhGRLM3HxwcgyXotkj7+/P8lUpxO513dRysUZyROp1n84KGHzNObH380AWfaNOja1YQeERFxK6fTyT333ENMTAwAAQEBOPT3cYpcLhd//vknV65cua3NMV0uF6dPnyYgIABv77uLJwo3GVGDBqabqlMnWLXKdFF9+61Z7fgWo9FFRCTtXd8J+3rAkZRZlsUff/yBv7//bYdALy8vihYtetfhUeEmo8qfH778EsaPh9deM4OMIyNhyRKzNo6IiLiNw+EgKCiI/PnzEx8fb3c5GVp8fDzr16/noYceSuzSSy1fX9/betqTEoWbjMzLC1591axk/Mwz8PPP8OCDMGkS9OypbioRETdzOp13PR7E0zmdTq5du0a2bNluO9ykFQ0ozgxq1TLdVI89Zvaj6t3bLPx3/rzdlYmIiGQ4CjeZRZ48ZtuGKVPM+jjLlpnuqW3b7K5MREQkQ1G4yUwcDnjpJdi8GYoXh8OHoU4dmDABXC67qxMREckQFG4yo6pVYdcuePppuHYNBgyAli3hzBm7KxMREbGdwk1mlTs3LFpkpof7+ZmZVaGhsH693ZWJiIjYSuEmM3M4oEcPM0W8TBk4ftyskTNyJCQk2F2diIiILRRuPMEDD8COHWbRP5cLhgyBpk3h5Em7KxMREXE7hRtPkSMHzJsHc+dCQACsXg2VKkFEhN2ViYiIuJXCjacJDzdPcSpWhJgY8wTntdfMwGMREZEswNZwk5CQwJAhQwgJCcHf358SJUowYsSIW255Pn/+fCpVqkRAQABBQUF07dqV33//3U1VZwLlypn1b3r0AMuC0aPNWJxjx+yuTEREJN3ZGm7Gjh3LjBkzmD59Ovv27WPs2LGMGzeOadOmpXjNpk2b6NSpE88++yw//PADH3/8MZGRkXTv3t2NlWcC/v5mJtWiRWazzY0bzWyqL76wuzIREZF0ZWu42bx5M61ataJFixYEBwfTpk0bmjRpQmRkZIrXbNmyheDgYF588UVCQkKoU6cOPXr0uOk1WVrbtmZNnCpV4OxZsx7Oyy+bbRxEREQ8kK3hplatWqxevZqDBw8CsGfPHjZu3MgjjzyS4jU1a9bk2LFjfPnll1iWxalTp1i6dCnNmzd3V9mZT8mSsGmTWd0YzMabderAoUNER8OwYRAdbWuFIiIiacbWXcEHDhxIbGwsZcuWxel0kpCQwKhRo+jQoUOK19SuXZv58+fTtm1brly5wrVr12jZsiVvvfVWsudfvXqVq1evJr6OjY0FzJbsab1t/fX7pfV904SXF4wfj6NOHZz//jeO7duxwsI49upihg9vSvPm8eTNa3eRqZOh29mDqJ3dQ+3sPmpr90ivdr6d+zmsW43eTUeLFi1iwIABjB8/ngoVKhAVFUWfPn2YNGkS4eHhyV7z448/0qhRI/r27UvTpk2Jjo5mwIABVKtWjffee++G84cNG8bw4cNvOL5gwQICAgLS/DNlBv6nT1N1wgTyHDjALsKowi4mvxlBSNk4u0sTERFJVlxcHO3bt+fChQvkypXrpufaGm6KFCnCwIED6dWrV+KxkSNH8tFHH7F///5kr+nYsSNXrlzh448/Tjy2ceNG6taty4kTJwgKCkpyfnJPbooUKcKZM2du2Ti3Kz4+noiICBo3boyPj0+a3jstRUfDyWPXcMyaxZ4P9tKdd5lVaBiVJneEYkUJDIR/NGOGklnaObNTO7uH2tl91NbukV7tHBsbS968eVMVbmztloqLi8PLK+mwH6fTiesmO1zHxcXh7Z20bKfTCZDsFHI/Pz/8/PxuOO7j45Nuf7jT895p4f33YfhwH+ClxGP/Pj4MnjbfDx1qxuFkdBm9nT2F2tk91M7uo7Z2j7Ru59u5l60Dilu2bMmoUaNYuXIlhw8fZvny5UyaNInHH3888ZxBgwbRqVOnJNcsW7aMGTNm8Ouvv7Jp0yZefPFFqlevTsGCBe34GJlOjx6wc6f5mj3bHJtdejw7qcxOKtNjf1+4fNneIkVERO6QrU9upk2bxpAhQ+jZsycxMTEULFiQHj168PrrryeeEx0dzdGjRxNfd+7cmYsXLzJ9+nRefvll7rnnHh5++GHGjh1rx0fIlIKCbux2qvxRPyqvvAxvvAGLd8N3q2DJErj/fnuKFBERuUO2hpucOXMyZcoUpkyZkuI5c+fOveHYCy+8wAsvvJB+hWVFTqfpi6pXDzp0gH37oFo1mDYNnn3W7EAuIiKSCWhvqSwuKMiMsUl8ktOgAURFmT2prlyB7t1N2Pn/KfQiIiIZncJNFhcUZB7YJOmmyp8fvvwS3nzTPNFZuNCscLxrl11lioiIpJrCjSTPywtefRXWr4ciReDnn6FmTZg+3WzGKSIikkEp3MjN1apluqkee8zsR/XCC/Dkk3DunN2ViYiIJEvhRm4tTx749FOYMgV8fGD5cggLg61b7a5MRETkBgo3kjoOh9l4c/NmKF4cjhyBunVhwgS4yaKLIiIi7qZwI7enalUzsPjpp+HaNRgwAFq2hDNn7K5MREQEULiRO5E7NyxaBDNngp+fmVlVqZIZfCwiImIzhRu5Mw6H2cchMhLKlIETJ8waOSNHQkKC3dWJiEgWpnAjd+eBB2DHDggPN2NvhgwxCwCePGl3ZSIikkUp3Mjdy5ED5s41XwEBsHq16aaKiLC7MhERyYIUbiTthIebrcYrVoSYGPME57XXzMBjERERN1G4kbRVtixs22bG41gWjB5txuIcO2Z3ZSIikkUo3Eja8/c3M6kWLYKcOWHjRggNhc8/t7syERHJAhRuJP20bQu7d5tNN8+eNVs49OtntnEQERFJJwo3kr5KlIBNm6BPH/N68mSoUwd+/dXWskRExHMp3Ej68/MzoebTT+Hee2H7drM31dKldlcmIiIeSOFG3KdVK7PDeM2aEBsLTz0FPXvClSt2VyYiIh5E4Ubcq2hRWLcOBg40r2fMgAcfhAMH7K1LREQ8hsKNuJ+PD4wZA199BfnywZ49ZtDxRx/ZXZmIiHgAhRuxT7Nmppuqfn24fBk6doSuXc33IiIid0jhRuxVsCD8738wbJjZjHPOHKhWDb7/3u7KREQkk1K4Efs5nTB0qNmTKigI9u0zAefdd80qxyIiIrdB4UYyjgYNTDdVs2ZmBlX37tChg5lZJSIikkoKN5Kx5M8PK1fC2LHmic7ChWaw8a5ddlcmIiKZhMKNZDxeXvDKK7B+vZk6/vPPZm2c6dPVTSUiIrekcCMZV61aZm+qVq3MflQvvABPPgnnztldmYiIZGAKN5Kx5ckDy5fDf/9r1sdZvhzv6tW5V4v+iYhIChRuJONzOODFF2HzZiheHMeRI9QZPBiviRPB5bK7OhERyWAUbiTzqFoVdu3C1aYNXgkJOAcNgkcfhdOn7a5MREQyEIUbyVxy5yZh/nyinn8eK1s2s4VDaKjZr0pERASFG8mMHA6ONG3KtY0boWxZOHECHn4YRoyAhAS7qxMREZsp3Ejm9cADsH07hIebsTevvw5NmkB0tN2ViYiIjRRuJHPLkQPmzjVfAQHw7bemmyoiwubCRETELgo34hnCw2HnTqhYEWJioGlTeO01uHbN7spERMTNFG7Ec5QtC9u2wXPPmZWMR4+G+vXh2DG7KxMRETdSuBHP4u8PM2bA4sWQKxds2mS6qT7/3O7KRETETRRuxDM9/bTZbLNKFTh7Fh57DPr1M9s4iIiIR1O4Ec9VooR5ctOnj3k9eTLUqQO//mprWSIikr4UbsSz+fmZUPPZZ3DvvWbqeFgYLF1qd2UiIpJOFG4ka3jsMYiKMjuNx8bCU09Bz55w5YrdlYmISBpTuJGso2hRWLsWBg40r2fMgAcfBO0wLiLiURRuJGvx8YExY+DrryFfPtizxww6/ugjuysTEZE0onAjWVPTpqabqkEDuHwZOnaErl3N9yIikqkp3EjWVbCg2aZh+HDw8oI5c6BaNdi71+7KRETkLijcSNbmdJoNN1evhqAg2LcPqleH2bPNKsciIpLpKNyIgNmmISoKmjUzM6j+/W9o397MrBIRkUxF4Ubkuvz5YeVKGDvWPNFZtAgqVzYrHYuISKahcCPyd15e8MorsH69mTr+yy9QsyZMm6ZuKhGRTELhRiQ5tWrB7t3QurXZj+rFF+GJJ+DcObsrExGRW1C4EUlJnjywbBlMnQq+vvDpp2brhq1b7a5MRERuQuFG5GYcDnjhBdi82WzEeeQI1K0L48eDy2V3dSIikgyFG5HUqFLFDCxu2xauXTPjch59FE6ftrsyERH5B4UbkdTKlQsWLoRZsyBbNvjqKwgNhXXr7K5MRET+RuFG5HY4HNC9O0RGQtmycOIEPPwwjBgBCQl2VyciIijciNyZihVhxw7o3NmMvXn9dWjSBKKj7a5MRCTLU7gRuVPZs5v9qObNM99/+63ppoqIsLsyEZEsTeFG5G516mSe4jzwAMTEmB3HX3vNDDwWERG3U7gRSQtly5r1b557zqxkPHq02a/q2DG7KxMRyXJsDTcJCQkMGTKEkJAQ/P39KVGiBCNGjMC6xTL3V69e5bXXXqNYsWL4+fkRHBzM+++/76aqRVLg7w8zZsDixWZm1aZNppvq88/trkxEJEvxtvPNx44dy4wZM5g3bx4VKlRgx44ddOnShdy5c/Piiy+meN3TTz/NqVOneO+99yhZsiTR0dG4tKCaZBRPP23WxWnbFnbuhMceg7594c03zUrHIiKSrmwNN5s3b6ZVq1a0aNECgODgYBYuXEhkZGSK13z99desW7eOX3/9lTx58iReJ5KhlChhntwMHAhTpsDkybBhg3mqU7y43dWJiHg0W8NNrVq1mDVrFgcPHqR06dLs2bOHjRs3MmnSpBSvWbFiBVWrVmXcuHF8+OGHZM+enccee4wRI0bg7+9/w/lXr17l6tWria9jY2MBiI+PJz4+Pk0/z/X7pfV9JalM085eXjBuHI66dXF264Zjxw6ssDASZs7EatPG7upuKdO0cyandnYftbV7pFc73879HNatBrikI5fLxeDBgxk3bhxOp5OEhARGjRrFoEGDUrymWbNmrF27lkaNGvH6669z5swZevbsSYMGDZgzZ84N5w8bNozhw4ffcHzBggUEBASk6ecRSYn/6dNUmTiR+/bvB+BQs2Z837UrLnVTiYikSlxcHO3bt+fChQvkypXrpufaGm4WLVrEgAEDGD9+PBUqVCAqKoo+ffowadIkwsPDk72mSZMmbNiwgZMnT5I7d24Ali1bRps2bbh8+fINT2+Se3JTpEgRzpw5c8vGuV3x8fFERETQuHFjfHx80vTe8pdM287x8XgNH45z3DgArIoVubZgAZQpY3Nhycu07ZzJqJ3dR23tHunVzrGxseTNmzdV4cbWbqkBAwYwcOBAnnnmGQAqVqzIkSNHGDNmTIrhJigoiEKFCiUGG4By5cphWRa//fYbpUqVSnK+n58ffn5+N9zHx8cn3f5wp+e95S+Zrp19fGDsWLNdQ8eOOPbuxefBB80Mq44d7a4uRZmunTMptbP7qK3dI63b+XbuZetU8Li4OLy8kpbgdDpvOvOpdu3anDhxgkuXLiUeO3jwIF5eXhQuXDjdahVJM02bwp490KABXL5sFgHs0sV8LyIid83WcNOyZUtGjRrFypUrOXz4MMuXL2fSpEk8/vjjiecMGjSITp06Jb5u37499913H126dOHHH39k/fr1DBgwgK5duyY7oFgkQwoKMts0DB9uBh7PnQvVqsHevXZXJiKS6dkabqZNm0abNm3o2bMn5cqVo3///vTo0YMRI0YknhMdHc3Ro0cTX+fIkYOIiAjOnz9P1apV6dChAy1btmTq1Kl2fASRO+d0mg03V682YWffPqheHWbPNqsci4jIHbF1zE3OnDmZMmUKU6ZMSfGcuXPn3nCsbNmyRGhzQvEU9eubbqpOneDrr+Hf/zabcL7zjlnpWEREbov2lhLJCPLlg5UrzYBjpxMWLYLKlWHXLrsrExHJdBRuRDIKLy945RWzknHRovDLL1CzJkybpm4qEZHboHAjktHUrAm7d0Pr1vDnn/Dii/DEE3DunN2ViYhkCgo3IhlRnjywbBlMnWo22/z0UwgLgy1b7K5MRCTDU7gRyagcDnjhBdi82WzEeeQI1K0L48bBTdaCEhHJ6hRuRDK6KlXMwOJnnoGEBHj1VWjRAk6ftrsyEZEMSeFGJDPIlQsWLDBr4GTLZqaMh4bCunV2VyYikuEo3IhkFg4HdOsGkZFQrhycOGH2qXrjDfNER0REAIUbkcynYkXYvt3sR+VywdCh0LgxREfbXZmISIagcCOSGWXPDu+/Dx98YL5fswYqVYJvvrG7MhER2ynciGRmHTvCzp3wwANmgHHTpjB4MFy7ZndlIiK2UbgRyezKlIGtW+H5583rMWPMflXHjtlaloiIXRRuRDyBvz+8/TYsWWJmVm3aZGZTff653ZWJiLidwo2IJ3nqKbN1Q9WqcPYsPPYY9OtntnEQEckiFG5EPE3x4ubJTd++5vXkyVC7Nvz6q711iYi4icKNiCfy9YVJk+Czz+Dee2HHDrM31ccf212ZiEi6U7gR8WSPPQZRUVCrFsTGwtNPQ8+ecOWK3ZWJiKQbhRsRT1e0KKxdC4MGmdczZkCNGnDggK1liYikF4UbkazAxwdGj4ZVqyB/fvjuO7Mh54cfpnhJdDQsXFhGCx+LSKajcCOSlTRpYrqpHn4YLl+GTp3MNg6XL99w6smTsHhxWU6edH+ZIiJ3Q+FGJKsJCjLbNLzxBnh5wdy5Zur43r12VyYikia87S5ARGzgdMKQIfDQQ9C+PezfD9WrE/3GbKIf7gAOB7t3OwDYvduB9///TREUZL5ERDIyhRuRrKxePdNNFR4OX33FO6/8zHAc//9D89fDc8/99dfE0KEwbJjbqxQRuS3qlhLJ6vLlgy++gHHj6OF8j51UZmfhVswcZBb9mznzGjt3mv05e/SwuVYRkVTQkxsRMWNvBgwgqE4dgp55Bo7uxhofDUQSFmpRubLdBYqIpJ6e3IjIX2rWNN1UrVvjuHYNAK9XXjH7VImIZBIKNyKS1L33wrJl5B/2HEMcb1Bo/SKzdcOWLXZXJiKSKgo3InIjh4MCg7vw8PjsBJbIDkePQt26MG4cuFx2VyciclMKNyKSogslS3Jt2zZ45hlISIBXX4UWLeD0abtLExFJkcKNiNxcrlywYAHMng3ZssHXX0NoKKxbZ3dlIiLJUrgRkVtzOKBbN4iMhHLl4MQJs4XDG2+YJzoiIhmIwo2IpF7FirB9u9mPyuUyq/o1box21xSRjEThRkRuT/bs8P778MEH5vs1a6BSJbNflYhIBqBwIyJ3pmNHs2xxpUpmgHHTpjBoEMTH212ZiGRxCjcicufKlIGtW6FnT/P6zTehfn0zdVxExCYKNyJyd7Jlg7fego8/NjOrNm82s6lWrLC7MhHJohRuRCRttGkDu3dDtWpw7hy0agV9+8Kff9pdmYhkMQo3IpJ2iheHjRuhXz/zesoUqF0bfvnF1rJEJGtRuBGRtOXrCxMnmm6pPHlgxw6oXBmWLLG7MhHJIhRuRCR9tGxpdhivXRtiY6FtW3j+efjjD7srExEPp3AjIumnSBFYuxYGDzarHM+cCQ8+CPv3212ZiHgwhRsRSV/e3jBqlNmTKn9++O47qFoVPvzQ7spExEMp3IiIezRpYrqpHn4YLl+GTp3MNg6XL9tdmYh4GIUbEXGfoCCzTcMbb4CXF8yda57i7N1rd2Ui4kEUbkTEvZxOGDIEvv0WChY042+qV4fZs8Gy7K5ORDyAwo2I2KNePdNN9cgjcOUK/Pvf0L69mVklInIXFG5ExD758sEXX8C4cWbg8aJFZk2cnTvtrkxEMjGFGxGxl5cXDBgAGzZAsWJmNeOaNWHqVHVTicgdUbgRkYzhwQfN3lSPPw7x8fDSS+b7s2ftrkxEMhmFGxHJOO69Fz75BKZNM9s4fPYZhIXBli12VyYimYjCjYhkLA4H9O5tAk3JknD0KNStC2PHgstld3Uikgko3IhIxnR9YPEzz0BCAgwcCM2bQ0yM3ZWJSAancCMiGVeuXLBggVkDJ1s2WLUKQkPNflUiIilQuBGRjM3hgG7dYPt2KFcOoqOhYUMYPtw80RER+QeFGxHJHO6/3wScLl3M2Jthw6BxYzhxwu7KRCSDUbgRkcwje3Z4/32zo3j27LBmjemmWrXK7spEJANRuBGRzOdf/zKDjStVgtOnoVkzGDTIrI8jIlnebYeb8PBw1q9fnx61iIikXpkysHUr9OxpXr/5JtSvb6aOi0iWdtvh5sKFCzRq1IhSpUoxevRojh8/fsdvnpCQwJAhQwgJCcHf358SJUowYsQIrFQuub5p0ya8vb0JDQ294xpEJBPLlg3eegs+/tjMrNq82XRTrVhhd2UiYqPbDjeffvopx48f5/nnn2fx4sUEBwfzyCOPsHTpUuJv85Hw2LFjmTFjBtOnT2ffvn2MHTuWcePGMW3atFtee/78eTp16kTDhg1v9yOIiKdp08Zs3VCtGpw7B61aQd++8OefdlcmIja4ozE3+fLlo1+/fuzZs4dt27ZRsmRJOnbsSMGCBenbty8//fRTqu6zefNmWrVqRYsWLQgODqZNmzY0adKEyMjIW1773HPP0b59e2rWrHknH0FEPE3x4rBxI/TrZ15PmQK1a5uNOEUkS/G+m4ujo6OJiIggIiICp9NJ8+bN2bt3L+XLl2fcuHH07dv3ptfXqlWLWbNmcfDgQUqXLs2ePXvYuHEjkyZNuul1c+bM4ddff+Wjjz5i5MiRNz336tWrXL16NfF1bGwsAPHx8bf9pOlWrt8vre8rSamd3SNTtrPDAW++iaNuXZzPPotjxw6sypVJmDED66mn7K4uWZmynTMptbV7pFc73879HFZqB7j87eYrVqxgzpw5fPPNNzzwwAN069aN9u3bkytXLgCWL19O165dOXfu3E3v5XK5GDx4MOPGjcPpdJKQkMCoUaMYNGhQitf89NNP1KlThw0bNlC6dGmGDRvGp59+SlRUVLLnDxs2jOHDh99wfMGCBQQEBKT+g4tIppLt9GmqTprEffv2AXCoaVO+79oVl5+fzZWJyJ2Ii4ujffv2XLhwITFvpOS2n9wEBQXhcrlo164dkZGRyQ7mbdCgAffcc88t77VkyRLmz5/PggULqFChAlFRUfTp04eCBQsSHh5+w/kJCQm0b9+e4cOHU7p06VTVO2jQIPpdf0yNeXJTpEgRmjRpcsvGuV3x8fFERETQuHFjfHx80vTe8he1s3t4RDt36EDCG2/gNXYsIatWEXz8ONcWLICyZe2uLJFHtHMmobZ2j/Rq5+s9L6lx2+Fm8uTJPPXUU2TLli3Fc+655x4OHTp0y3sNGDCAgQMH8swzzwBQsWJFjhw5wpgxY5INNxcvXmTHjh3s3r2b3r17A+bpj2VZeHt788033/Dwww8nucbPzw+/ZP6l5uPjk25/uNPz3vIXtbN7ZOp29vGBMWPg4YfhX//C8f33+Dz4IMyYAZ062V1dEpm6nTMZtbV7pHU73869bjvcdOzY8XYvSVFcXBxeXknHNDudTlwuV7Ln58qVi7179yY59vbbb/Ptt9+ydOlSQkJC0qw2EfEgjRvDnj3QoQN8+y2Eh5v/nT4dcuSwuzoRSWN3NaD4brVs2ZJRo0ZRtGhRKlSowO7du5k0aRJdu3ZNPGfQoEEcP36cDz74AC8vL+6///4k98ifPz/ZsmW74biISBKBgfDNN+ZJztChMG8ebNsGixfDAw/YXZ2IpCFbt1+YNm0abdq0oWfPnpQrV47+/fvTo0cPRowYkXhOdHQ0R7XiqIikBacT/vMfsydVwYKwfz/UqAGzZsHtza0QkQzM1nCTM2dOpkyZwpEjR/jjjz/45ZdfGDlyJL6+vonnzJ07l7Vr16Z4j2HDhqU4U0pEJFkPPQRRUfDII3DlCvToAe3awW0MWBSRjEsbZ4pI1pQvH3zxBYwfD97epnuqcmWzIaeIZGoKNyKSdXl5Qf/+sGEDFCtmVjOuWROmTlU3lUgmpnAjIvLgg2Zvqscfh/h4eOkl8/3Zs3ZXJiJ3QOFGRATg3nvhk09g2jTw9YXPPoOwMNiyxe7KROQ2KdyIiFzncEDv3ibQlCwJR49C3bowdiyksP6WiGQ8CjciIv90fWBxu3aQkAADB0Lz5hATY3dlIpIKCjciIsnJlQvmz4d33wV/f1i1CkJD4SZLU4hIxqBwIyKSEocDnn0WIiOhfHmIjoaGDWH4cPNER0QyJIUbEZFbuf9+E3C6djVjb4YNg0aN4MQJuysTkWQo3IiIpEb27PDee/Dhh+b7tWtNN9WqVXZXJiL/oHAjInI7/vUv2LULKlWC06ehWTMz4Dg+3u7KROT/KdyIiNyu0qVh61bo2dO8HjsW6tUzU8dFxHYKNyIidyJbNnjrLfj4YzOzassW00312Wd2VyaS5SnciIjcjTZtzNYN1arBuXPQujX06QNXr9pdmUiWpXAjInK3iheHjRuhXz/z+r//hdq1zUacIuJ2CjciImnB1xcmToTPP4c8ecwKx2FhsGSJ3ZWJZDkKNyIiaenRRyEqCurUgYsXoW1beO45+OMPuysTyTIUbkRE0lqRIrBmDbz2mlnl+J13oEYN2L/f7spEsgSFGxGR9ODtDSNHmkX+8ueHvXuhShX44AO7KxPxeAo3IiLpqXFj2LPH7EkVFwfh4dC5M1y6ZHdlIh5L4UZEJL0FBponOCNHgpcXzJtnpo5/953dlYl4JIUbERF3cDrNGJw1a6BQITP+pnp1vGbPBsuyuzoRj6JwIyLiTg89ZGZTNW8OV6/i7NWLqhMmwIULdlcm4jEUbkRE3C1vXrMezoQJWN7eFNq0Ce8aNWDHDrsrE/EICjciInbw8oKXXyZh7Vri8uXD8euvUKuWWd1Y3VQid0XhRkTERlb16qydPBlX69YQH2/2pWrdGs6etbkykcxL4UZExGbxOXKQsHgxTJ9utnFYscLsML55s92liWRKCjciIhmBwwG9esHWrVCyJBw7ZgYfjx0LLpfd1YlkKgo3IiIZSVgY7NoF7dpBQgIMHGhmVsXE2F2ZSKahcCMiktHkzAnz58O774K/v1kAMDQU1q61uzKRTEHhRkQkI3I44NlnITISypeH6GizhcPw4eaJjoikSOFGRCQju/9+E3C6djVjb4YNg0aN4MQJuysTybAUbkREMrrs2eG99+CjjyBHDtM9FRpquqtE5AYKNyIimUWHDrBzpwk2p09Ds2ZmwHF8vN2ViWQoCjciIplJ6dKwZYuZNg5mqni9enDkiL11iWQgCjciIplNtmxmwb+lSyF3bhN2wsLgs8/srkwkQ1C4ERHJrJ58EnbvhurV4dw5s21Dnz5w9ardlYnYSuFGRCQzCwmBDRvg5ZfN6//+12zA+fPP9tYlYiOFGxGRzM7XFyZMgM8/hzx5zArHlSvD4sV2VyZiC4UbERFP8eijsGcP1KkDFy/CM89Ajx7wxx92VybiVgo3IiKepHBhWLMGXnvNrHI8axbUqAH79tldmYjbKNyIiHgab28YOdIs8pc/P+zdC1Wrwrx5dlcm4hYKNyIinqpxY9NN1bAhxMVB584QHg6XLtldmUi6UrgREfFkgYHmCc7IkeDlBR98YJ7ifPed3ZWJpBuFGxERT+d0mjE4a9ZAoUJw4IBZG+edd8Cy7K5OJM0p3IiIZBUPPQRRUdC8uVno77nnzIyqCxfsrkwkTSnciIhkJXnzmvVwJkwwA4+XLDFr4uzYYXdlImlG4UZEJKvx8jIrGm/cCMHB8OuvZlXj//5X3VTiERRuRESyqho1zN5UTzwB8fFmX6rWreHsWbsrE7krCjciIlnZPfeY3cXfests47BiBYSGwqZNdlcmcscUbkREsjqHA3r2hK1boVQpOHYM6tWDN98El8vu6kRum8KNiIgYYWGwcye0bw8JCTBoEDzyCMTE2F2ZyG1RuBERkb/kzAkffQTvvQf+/vDNN1CpklkjRySTULgREZGkHA7o2hW2b4fy5eHkSbOFw7Bh5omOSAancCMiIsmrUMEEnGefNVPEhw+HRo3gxAm7KxO5KYUbERFJWUAAvPuu6arKkQPWrjXdVF9/bXdlIilSuBERkVvr0MEMNg4NhTNnzEDjgQPN+jgiGYzCjYiIpE7p0rBlC/TqZV6PHWumjB85Ym9dIv+gcCMiIqmXLRtMn24W/sud24SdsDD47DO7KxNJZGu4SUhIYMiQIYSEhODv70+JEiUYMWIE1k32Nlm2bBmNGzcmX7585MqVi5o1a7Jq1So3Vi0iIjz5pNm6oXp1OHfObNvQp4/ZbVzEZraGm7FjxzJjxgymT5/Ovn37GDt2LOPGjWPatGkpXrN+/XoaN27Ml19+yc6dO2nQoAEtW7Zk9+7dbqxcREQICYENG8wmnGA23qxVC37+2d66JMvztvPNN2/eTKtWrWjRogUAwcHBLFy4kMjIyBSvmTJlSpLXo0eP5rPPPuPzzz8nLCwsPcsVEZF/8vWFCROgQQMID4ddu6ByZZg1C555xu7qJIuyNdzUqlWLWbNmcfDgQUqXLs2ePXvYuHEjkyZNSvU9XC4XFy9eJE+ePMn+/OrVq1z922PS2NhYAOLj44lP41H+1++X1veVpNTO7qF2dg+PaecmTWD7dpydOuG1cSO0a0fC6tW4Jk40Kx1nAB7T1hlcerXz7dzPYd1sgEs6c7lcDB48mHHjxuF0OklISGDUqFEMGjQo1fcYN24cb775Jvv37yd//vw3/HzYsGEMHz78huMLFiwgICDgruoXEZGkHAkJlFm0iNJLl+KwLGKLFmX7gAFcKlLE7tIkk4uLi6N9+/ZcuHCBXLly3fRcW8PNokWLGDBgAOPHj6dChQpERUXRp08fJk2aRHh4+C2vX7BgAd27d+ezzz6jUaNGyZ6T3JObIkWKcObMmVs2zu2Kj48nIiKCxo0b4+Pjk6b3lr+ond1D7ewentrOjtWrcXbujOPUKayAABKmTsXq1MnWmjy1rTOa9Grn2NhY8ubNm6pwY2u31IABAxg4cCDP/H+/bMWKFTly5Ahjxoy5ZbhZtGgR3bp14+OPP04x2AD4+fnh5+d3w3EfH590+8OdnveWv6id3UPt7B4e187NmkFUFPzrXzhWr8a7WzdYtw7eftusdGwjj2vrDCqt2/l27mXrbKm4uDi8vJKW4HQ6cblcN71u4cKFdOnShYULFyYORhYRkQwmMBBWrYKRI8HLCz78EKpWhT177K5MPJyt4aZly5aMGjWKlStXcvjwYZYvX86kSZN4/PHHE88ZNGgQnf72KHPBggV06tSJiRMnUqNGDU6ePMnJkye5cOGCHR9BRERuxumE114ze1IVKgQHDkCNGjBzptmMUyQd2Bpupk2bRps2bejZsyflypWjf//+9OjRgxEjRiSeEx0dzdGjRxNfz5o1i2vXrtGrVy+CgoISv1566SU7PoKIiKRG3bqmm6pFC7PQ3/PPQ9u2oH+YSjqwdcxNzpw5mTJlyg1r1/zd3Llzk7xeu3ZtutYkIiLpJG9eWLECJk82m25+/LHZjHPRIqhWze7qxINobykREXEfLy+zovHGjRAcDL/+CrVrw5Qp6qaSNKNwIyIi7lejhtmb6oknID4e+vY1+1OdPWt3ZeIBFG5ERMQe99xjdhd/6y2zjcOKFRAaCps22V2ZZHIKNyIiYh+HA3r2hK1boVQpOHYM6tWDN9+EWywLIpIShRsREbFfWJgZXNyhAyQkwKBB8MgjEBNjd2WSCSnciIhIxpAzp1no7/33zWab33wDlSrBt9/aXZlkMgo3IiKScTgc0KULbN8OFSrAyZPQqBEMHWqe6IikgsKNiIhkPBUqQGQkdOtmpoi/8QY0bAjHj9tdmWQCCjciIpIxBQTA7Nkwf77ZbHPdOjOb6quv7K5MMjiFGxERydjat4ddu0ywOXMGmjeHV1816+OIJEPhRkREMr5SpWDLFujd27weNw4eegiOHLG3LsmQFG5ERCRzyJYNpk2DTz6B3LnN2jihofDpp3ZXJhmMwo2IiGQuTzxhtm6oXh3On4fHH4eXXjK7jYugcCMiIplRSAhs2AD9+5vXU6dCrVrw88/21iUZgsKNiIhkTr6+MH48fPEF3HefGXRcuTIsWmR3ZWIzhRsREcncWrSAqCioWxcuXoR27aBHD/jjD7srE5so3IiISOZXuLDZpuE//zGrHM+aZcbk7Ntnd2ViA4UbERHxDN7eMGKE2ZOqQAH4/nuoWhXmzbO7MnEzhRsREfEsjRqZbqpGjSAuDjp3hk6d4NIluysTN1G4ERERzxMYCKtWwahR4OUFH36I94MPkuvQIbsrEzdQuBEREc/k5QWDB8PatVCoEI6DB3nolVfwmjXLbMYpHkvhRkREPFvduhAVhat5c5zx8Th794a2beHCBbsrk3SicCMiIp4vb14Sli/n+y5dsLy94eOPISwMtm+3uzJJBwo3IiKSNTgc/NKqFQnr1kFwMBw6BLVrw+TJ6qbyMAo3IiKSpVjVqpm9qZ58EuLjoV8/aNUKfv/d7tIkjSjciIhI1nPPPaZr6q23zDYOn39udhjftMnuyiQNKNyIiEjW5HBAz56wbRuUKgW//Qb16sGYMeBy2V2d3AWFGxERydpCQ2HnTujQARISzPTxRx6BU6fsrkzukMKNiIhIzpzw4Yfw/vvg72+2cAgNNftVSaajcCMiIgKmm6pLFzM9vEIFOHnSbOEwdKh5oiOZhsKNiIjI31WoAJGR0K2bmSL+xhvQsCEcP253ZZJKCjciIiL/FBAAs2fD/PmQIwesW2e6qb76yu7KJBUUbkRERFLSvj3s2mVWMz5zBpo3h1dfNevjSIalcCMiInIzpUrBli3wwgvm9bhx8NBDcPiwrWVJyhRuREREbsXPD6ZOhU8+MQsAbt1qnuYsX253ZZIMhRsREZHUeuIJs3VDjRpw/rx5/eKLcPWq3ZXJ3yjciIiI3I7gYNiwAQYMMK+nTYNateDnn20tS/6icCMiInK7fHzM2JuVK+G++8yg48qVYdEiuysTFG5ERETuXPPmEBUFdevCxYvQrh38+98QF2d3ZVmawo2IiMjdKFzYbNMwZIhZ5Xj2bDMm58cf7a4sy1K4ERERuVve3mYl44gIKFAAvv8eqlWDuXPNKsfiVgo3IiIiaaVhQ9izx+xJFRdn9qoKD4dLl+yuLEtRuBEREUlLBQrAqlUwahR4eZndxqtUMaFH3ELhRkREJK15ecHgwbB2LRQqBAcPmnE4M2eqm8oNFG5ERETSS926ZjZVixZmob/nn4e2beHCBbsr82gKNyIiIukpb174/HOYONEMPP74Y7N1w/btdlfmsRRuRERE0pvDAf36waZNEBIChw5B7dowebK6qdKBwo2IiIi7VK9uVjNu0wbi403gadUKfv/d7so8isKNiIiIO91zDyxZAm+/bXYb//xzCA2FjRvtrsxjKNyIiIi4m8NhBhdv3QqlS8Nvv0H9+jBmDLhcdleX6SnciIiI2CU0FHbsgA4dICHBTB9v1gxOnbK7skxN4UZERMROOXOahf7efx/8/c0WDqGhsHq13ZVlWgo3IiIidnM4zFYNO3ZAhQpw8iQ0bgyvvw7XrtldXaajcCMiIpJRlC8PkZHQrZuZIj5ihNmv6vhxuyvLVBRuREREMpKAAJg9G+bPhxw5YP1600315Zd2V5ZpKNyIiIhkRO3bmzVxwsLgzBmzhcMrr5j1ceSmFG5EREQyqlKlYMsWeOEF83r8eHjoITh82NayMjqFGxERkYzMzw+mToVPPjELAG7dap7mLF9ud2UZlq3hJiEhgSFDhhASEoK/vz8lSpRgxIgRWLfYZ2Pt2rVUrlwZPz8/SpYsydy5c91TsIiIiF2eeAJ274YaNeD8efP6xRfNbuOShK3hZuzYscyYMYPp06ezb98+xo4dy7hx45g2bVqK1xw6dIgWLVrQoEEDoqKi6NOnD926dWPVqlVurFxERMQGwcGwYQMMGGBeT5sGtWrBzz/bWlZG423nm2/evJlWrVrRokULAIKDg1m4cCGRkZEpXjNz5kxCQkKYOHEiAOXKlWPjxo1MnjyZpk2buqVuERER2/j4wLhxZruGTp3MoOPKlWHWLHjmGburyxBsDTe1atVi1qxZHDx4kNKlS7Nnzx42btzIpEmTUrxmy5YtNGrUKMmxpk2b0qdPn2TPv3r1Klf/9sguNjYWgPj4eOLTeMT59ful9X0lKbWze6id3UPt7D4e19aNG8OOHTg7dcJrwwZo1w5XRAQJkyaZ6eQ2Sa92vp372RpuBg4cSGxsLGXLlsXpdJKQkMCoUaPo0KFDitecPHmSAgUKJDlWoEABYmNj+eOPP/D390/yszFjxjB8+PAb7vPNN98QkE7/8SMiItLlvpKU2tk91M7uoXZ2H09ra0efPpQuWJAyS5bg9f77XPrf/9gxYAAXixSxta60bue4uLhUn2truFmyZAnz589nwYIFVKhQIXEMTcGCBQkPD0+T9xg0aBD9+vVLfB0bG0uRIkVo0qQJuXLlSpP3uC4+Pp6IiAgaN26Mj49Pmt5b/qJ2dg+1s3uond3Ho9u6ZUsSunTB2bkzuY4epcErr5AwdSpWp05mawc3Sq92vt7zkhq2hpsBAwYwcOBAnvn/PsKKFSty5MgRxowZk2K4CQwM5NQ/dks9deoUuXLluuGpDYCfnx9+fn43HPfx8Um3P9zpeW/5i9rZPdTO7qF2dh+PbeumTSEqCjp2xBERgXf37rBuHbz9ttmc083Sup1v5162zpaKi4vDyytpCU6nE5fLleI1NWvWZPU/dkqNiIigZs2a6VKjiIhIplGgAHz9NYweDU4nfPQRVK1qQk8WYmu4admyJaNGjWLlypUcPnyY5cuXM2nSJB5//PHEcwYNGkSnTp0SXz/33HP8+uuvvPLKK+zfv5+3336bJUuW0LdvXzs+goiISMbi5QWDBsHatVC4MBw8CA8+CDNmmM04swBbw820adNo06YNPXv2pFy5cvTv358ePXowYsSIxHOio6M5evRo4uuQkBBWrlxJREQElSpVYuLEibz77ruaBi4iIvJ3deqYJzaPPmoW+uvZE55+2iwA6OFsHXOTM2dOpkyZwpQpU1I8J7nVh+vXr8/u3bvTrzARERFPcN99sGIFTJkCr74KS5fCzp2weDFUq2Z3delGe0uJiIh4MocD+vaFjRshJAQOHYLatWHyZI/tplK4ERERyQqqVzerGbdpA/Hx0K8fPPYY/P673ZWlOYUbERGRrOKee2DJEjM93M8PvvgCQkPNUx0PonAjIiKSlTgc8PzzsHUrlC4Nv/1m9qkaMwZushRLZqJwIyIikhWFhsKOHdChAyQkwODB0KwZ/GOh3MxI4UZERCSrypkTPvwQ3n/fbLYZEWFCzz8Wy81sFG5ERESyMocDunSB7dvh/vvh5Emz4/jrr8O1a3ZXd0cUbkRERATKl4dt26B7dzNFfMQIaNgQjh+3u7LbpnAjIiIiRkAAzJoFCxZAjhywfj1UqgRffml3ZbdF4UZERESSatfOrIkTFmbWwWnRAgYMMOvjZAIKNyIiInKjUqVgyxZ44QXzesIEqFsXDh+2tazUULgRERGR5Pn5wdSpsGyZWQBw2zbzNGfZMrsruymFGxEREbm5xx+H3buhRg2zq/iTT5onOleu2F1ZshRuRERE5NaCg2HDBjP2BmD6dKhVC376ydaykqNwIyIiIqnj4wPjxpnZU3nzmqc5lSvDwoV2V5aEwo2IiIjcnkcegagoeOghuHQJ2rc36+PExdldGaBwIyIiIneiUCGzTcPrr5tVjt99F6pXJ3rdQRYuLEN0tH2lKdyIiIjInfH2huHDzZ5UgYHwww+cfrQ7ixeX5aTCjYiIiGRaDRuabqrGjXFc/f8ZVKdP21aOt23vLCIiIh4hOhqiowvA6K/ZEfAlfAa7j+fHe5f5eVCQ+XIXhRsRERG5K++8Y3qnTIfQowA899xfEWPoUBg2zH31KNyIiIjIXenRAx57zHy/ffs1nnvOm5kzr1GtmokZ7nxqAwo3IiIicpf+3u107ZoFQFiYReXK9tSjAcUiIiLiURRuREREJM0EBkLbtvsJDLSvBoUbERERSTNBQdCu3QG3j7P5O4UbERER8SgKNyIiIuJRFG5ERETEoyjciIiIiEdRuBERERGPonAjIiIiHkXhRkRERDyKwo2IiIh4FIUbERER8SgKNyIiIuJRstyu4JZldiuNjY1N83vHx8cTFxdHbGwsPj4+aX5/MdTO7qF2dg+1s/uord0jvdr5+u/t67/HbybLhZuLFy8CUKRIEZsrERERkdt18eJFcufOfdNzHFZqIpAHcblcnDhxgpw5c+JwONL03rGxsRQpUoRjx46RK1euNL23/EXt7B5qZ/dQO7uP2to90qudLcvi4sWLFCxYEC+vm4+qyXJPbry8vChcuHC6vkeuXLn0fxw3UDu7h9rZPdTO7qO2do/0aOdbPbG5TgOKRURExKMo3IiIiIhHUbhJQ35+fgwdOhQ/Pz+7S/Foamf3UDu7h9rZfdTW7pER2jnLDSgWERERz6YnNyIiIuJRFG5ERETEoyjciIiIiEdRuBERERGPonCTRt566y2Cg4PJli0bNWrUIDIy0u6SPM6YMWOoVq0aOXPmJH/+/LRu3ZoDBw7YXZbHe/PNN3E4HPTp08fuUjzO8ePH+de//sV9992Hv78/FStWZMeOHXaX5VESEhIYMmQIISEh+Pv7U6JECUaMGJGq/Ynk5tavX0/Lli0pWLAgDoeDTz/9NMnPLcvi9ddfJygoCH9/fxo1asRPP/3kltoUbtLA4sWL6devH0OHDmXXrl1UqlSJpk2bEhMTY3dpHmXdunX06tWLrVu3EhERQXx8PE2aNOHy5ct2l+axtm/fzjvvvMMDDzxgdyke59y5c9SuXRsfHx+++uorfvzxRyZOnMi9995rd2keZezYscyYMYPp06ezb98+xo4dy7hx45g2bZrdpWV6ly9fplKlSrz11lvJ/nzcuHFMnTqVmTNnsm3bNrJnz07Tpk25cuVK+hdnyV2rXr261atXr8TXCQkJVsGCBa0xY8bYWJXni4mJsQBr3bp1dpfikS5evGiVKlXKioiIsOrVq2e99NJLdpfkUV599VWrTp06dpfh8Vq0aGF17do1ybEnnnjC6tChg00VeSbAWr58eeJrl8tlBQYGWuPHj088dv78ecvPz89auHBhutejJzd36c8//2Tnzp00atQo8ZiXlxeNGjViy5YtNlbm+S5cuABAnjx5bK7EM/Xq1YsWLVok+bMtaWfFihVUrVqVp556ivz58xMWFsbs2bPtLsvj1KpVi9WrV3Pw4EEA9uzZw8aNG3nkkUdsrsyzHTp0iJMnTyb5+yN37tzUqFHDLb8bs9zGmWntzJkzJCQkUKBAgSTHCxQowP79+22qyvO5XC769OlD7dq1uf/+++0ux+MsWrSIXbt2sX37drtL8Vi//vorM2bMoF+/fgwePJjt27fz4osv4uvrS3h4uN3leYyBAwcSGxtL2bJlcTqdJCQkMGrUKDp06GB3aR7t5MmTAMn+brz+s/SkcCOZUq9evfj+++/ZuHGj3aV4nGPHjvHSSy8RERFBtmzZ7C7HY7lcLqpWrcro0aMBCAsL4/vvv2fmzJkKN2loyZIlzJ8/nwULFlChQgWioqLo06cPBQsWVDt7MHVL3aW8efPidDo5depUkuOnTp0iMDDQpqo8W+/evfniiy9Ys2YNhQsXtrscj7Nz505iYmKoXLky3t7eeHt7s27dOqZOnYq3tzcJCQl2l+gRgoKCKF++fJJj5cqV4+jRozZV5JkGDBjAwIEDeeaZZ6hYsSIdO3akb9++jBkzxu7SPNr13392/W5UuLlLvr6+VKlShdWrVycec7lcrF69mpo1a9pYmeexLIvevXuzfPlyvv32W0JCQuwuySM1bNiQvXv3EhUVlfhVtWpVOnToQFRUFE6n0+4SPULt2rVvWMrg4MGDFCtWzKaKPFNcXBxeXkl/1TmdTlwul00VZQ0hISEEBgYm+d0YGxvLtm3b3PK7Ud1SaaBfv36Eh4dTtWpVqlevzpQpU7h8+TJdunSxuzSP0qtXLxYsWMBnn31Gzpw5E/ttc+fOjb+/v83VeY6cOXPeMI4pe/bs3HfffRrflIb69u1LrVq1GD16NE8//TSRkZHMmjWLWbNm2V2aR2nZsiWjRo2iaNGiVKhQgd27dzNp0iS6du1qd2mZ3qVLl/j5558TXx86dIioqCjy5MlD0aJF6dOnDyNHjqRUqVKEhIQwZMgQChYsSOvWrdO/uHSfj5VFTJs2zSpatKjl6+trVa9e3dq6davdJXkcINmvOXPm2F2ax9NU8PTx+eefW/fff7/l5+dnlS1b1po1a5bdJXmc2NhY66WXXrKKFi1qZcuWzSpevLj12muvWVevXrW7tExvzZo1yf6dHB4eblmWmQ4+ZMgQq0CBApafn5/VsGFD68CBA26pzWFZWqZRREREPIfG3IiIiIhHUbgRERERj6JwIyIiIh5F4UZEREQ8isKNiIiIeBSFGxEREfEoCjciIiLiURRuRESSUb9+ffr06WN3GSJyB7SIn4iku/r16xMaGsqUKVPsLiXVzp49i4+PDzlz5rS7FBG5TdpbSkQkGXny5LG7BBG5Q+qWEpF01blzZ9atW8d///tfHA4HDoeDw4cPJ3vu1atX6d+/P4UKFSJ79uzUqFGDtWvXJv587ty53HPPPaxatYpy5cqRI0cOmjVrRnR0dIrvv3btWhwOB6tWrSIsLAx/f38efvhhYmJi+OqrryhXrhy5cuWiffv2xMXFJV73z26p4OBgRo8eTdeuXcmZMydFixbVJpciGZTCjYikq//+97/UrFmT7t27Ex0dTXR0NEWKFEn23N69e7NlyxYWLVrEd999x1NPPUWzZs346aefEs+Ji4tjwoQJfPjhh6xfv56jR4/Sv3//W9YxbNgwpk+fzubNmzl27BhPP/00U6ZMYcGCBaxcuZJvvvmGadOm3fQeEydOpGrVquzevZuePXvy/PPPc+DAgdtrEBFJdwo3IpKucufOja+vLwEBAQQGBhIYGIjT6bzhvKNHjzJnzhw+/vhj6tatS4kSJejfvz916tRhzpw5iefFx8czc+ZMqlatSuXKlenduzerV6++ZR0jR46kdu3ahIWF8eyzz7Ju3TpmzJhBWFgYdevWpU2bNqxZs+am92jevDk9e/akZMmSvPrqq+TNm/eW14iI+2nMjYhkCHv37iUhIYHSpUsnOX716lXuu+++xNcBAQGUKFEi8XVQUBAxMTG3vP8DDzyQ+H2BAgUICAigePHiSY5FRkam+h4Oh4PAwMBUvbeIuJfCjYhkCJcuXcLpdLJz584bnuzkyJEj8XsfH58kP3M4HKRm0uffr3M4HMnex+Vypfoeqb1GRNxP4UZE0p2vry8JCQk3PScsLIyEhARiYmKoW7eumyoTEU+kMTciku6Cg4PZtm0bhw8f5syZM8k+7ShdujQdOnSgU6dOLFu2jEOHDhEZGcmYMWNYuXKlDVWLSGalcCMi6a5///44nU7Kly9Pvnz5OHr0aLLnzZkzh06dOvHyyy9TpkwZWrduzfbt2ylatKibKxaRzEwrFIuIiIhH0ZMbERER8SgKNyIiIuJRFG5ERETEoyjciIiIiEdRuBERERGPonAjIiIiHkXhRkRERDyKwo2IiIh4FIUbERER8SgKNyIiIuJRFG5ERETEoyjciIiIiEf5P+vBsTduMrh1AAAAAElFTkSuQmCC\n"},"metadata":{}},{"output_type":"stream","name":"stdout","text":["pente =  -9.64e-02 /min   y_0 =  8.89\n","k_moy = 9.49e-02  avec u(k_moy) =  1.2e-03 par calcul direct\n"]}]}]}