{
  "metadata": {
    "language_info": {
      "codemirror_mode": {
        "name": "python",
        "version": 3
      },
      "file_extension": ".py",
      "mimetype": "text/x-python",
      "name": "python",
      "nbconvert_exporter": "python",
      "pygments_lexer": "ipython3",
      "version": "3.8"
    },
    "kernelspec": {
      "name": "python",
      "display_name": "Python (Pyodide)",
      "language": "python"
    }
  },
  "nbformat_minor": 4,
  "nbformat": 4,
  "cells": [
    {
      "cell_type": "markdown",
      "source": "# Calculando la gravedad con el movimiento de un pendulo",
      "metadata": {}
    },
    {
      "cell_type": "code",
      "source": "import numpy as np\nimport matplotlib.pyplot as plt\nfrom scipy.interpolate import interp1d\nfrom scipy.optimize import curve_fit\nimport math",
      "metadata": {
        "trusted": true
      },
      "execution_count": 38,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "L = np.loadtxt('Longitudes.txt')\nM1= np.loadtxt('Tiempo 1.txt')",
      "metadata": {
        "trusted": true
      },
      "execution_count": 9,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "plt.figure()\nplt.plot(M1,L, '.') \nplt.xlabel(r'$Tiempo$')\nplt.ylabel(r'$Longitud$')\nplt.show()",
      "metadata": {
        "trusted": true
      },
      "execution_count": 44,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": "<Figure size 640x480 with 1 Axes>",
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAG0CAYAAADO5AZFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAvC0lEQVR4nO3de3RU5aH+8WcngWSIZBACkxlIICgXIUC4GSOIVeLhKLJEaoVKherReMEeo1wKVvBYkShVqniBo0vBLkFoe8C6PBqLUVEUIkISIypySSBKEhs0MzGEYJP9+8PD/By5lAyT7M3O97PWXqt5Z8/Ok7c68/juPXsM0zRNAQAAOESU1QEAAAAiiXIDAAAchXIDAAAchXIDAAAchXIDAAAchXIDAAAchXIDAAAcJcbqAK2tqalJBw4cUMeOHWUYhtVxAADAKTBNU7W1tfL5fIqKOvnaTJsrNwcOHFBycrLVMQAAQBjKy8vVo0ePk+7T5spNx44dJf0wOQkJCRanAQAApyIQCCg5OTn4Pn4yba7cHD0VlZCQQLkBAOAMcyqXlHBBMQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAcBTKDQAAiIgKf70+2FOtCn+9pTna3BdnAgCAyFu7db/mrStRkylFGVLupEGaPDLFkiys3AAAgNNS4a8PFhtJajKle9Z9YtkKDuUGAACcltLqumCxOarRNFVWfciSPJQbAABwWlIT4xVlhI5FG4Z6JXawJA/lBgAAnBav26XcSYMUbfzQcKINQ4smpcnrdlmShwuKAQDAaZs8MkVj+nZVWfUh9UrsYFmxkSg3AAAgQrxul6Wl5ihOSwEAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEexVblpbGzU/PnzlZqaKpfLpXPOOUcPPPCATPP/f2GFaZpasGCBvF6vXC6XsrKytGvXLgtTAwAAO7FVuXn44Ye1bNkyPfnkk/rss8/08MMPa/HixXriiSeC+yxevFhLly7V8uXLVVBQoPj4eI0bN06HDx+2MDkAALALw/zxsojFrrzySnk8Hj333HPBsZ///OdyuVx68cUXZZqmfD6fZs6cqVmzZkmS/H6/PB6PVq5cqSlTphxzzIaGBjU0NAR/DgQCSk5Olt/vV0JCQsv/UQAA4LQFAgG53e5Tev+21crNhRdeqPz8fH3xxReSpOLiYm3atEmXX365JKm0tFSVlZXKysoKPsftdisjI0ObN28+7jFzc3PldruDW3Jycsv/IQAAwDK2+m6puXPnKhAIqH///oqOjlZjY6MefPBBTZ06VZJUWVkpSfJ4PCHP83g8wcd+at68ebr77ruDPx9duQEAAM5kq3Lz5z//WatWrdLq1as1cOBAFRUVKScnRz6fT9OnTw/rmLGxsYqNjY1wUgAAYFe2KjezZ8/W3Llzg9fODBo0SPv27VNubq6mT5+upKQkSVJVVZW8Xm/weVVVVUpPT7ciMgAAsBlbXXNz6NAhRUWFRoqOjlZTU5MkKTU1VUlJScrPzw8+HggEVFBQoMzMzFbNCgAA7MlWKzcTJkzQgw8+qJSUFA0cOFCFhYVasmSJbrzxRkmSYRjKycnRwoUL1adPH6Wmpmr+/Pny+XyaOHGiteEBAIAt2KrcPPHEE5o/f75uv/12ff311/L5fLrlllu0YMGC4D5z5sxRXV2dsrOzVVNTo9GjRysvL09xcXEWJgcAAHZhq/vctIbmfE4eAADYwxl7nxsAAIDTRbkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAMAGKvz1+mBPtSr89VZHOePFWB0AAIC2bu3W/Zq3rkRNphRlSLmTBmnyyBSrY52xWLkBAMBCFf76YLGRpCZTumfdJ6zgnAbKDQAAFiqtrgsWm6MaTVNl1YesCeQAlBsAACyUmhivKCN0LNow1CuxgzWBHIByAwCAhbxul3InDVK08UPDiTYMLZqUJq/bZXGyMxcXFAMAYLHJI1M0pm9XlVUfUq/EDhSb00S5AQDABrxuF6UmQjgtBQAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHIVyAwAAHMVW5aZXr14yDOOYbcaMGZIk0zS1YMECeb1euVwuZWVladeuXRanBgAAdmKrcrN161ZVVFQEtw0bNkiSfvGLX0iSFi9erKVLl2r58uUqKChQfHy8xo0bp8OHD1sZGwAA2IhhmqZpdYgTycnJ0auvvhpcnfH5fJo5c6ZmzZolSfL7/fJ4PFq5cqWmTJlySscMBAJyu93y+/1KSEhosewAACBymvP+bauVmx87cuSIXnzxRd14440yDEOlpaWqrKxUVlZWcB+3262MjAxt3rz5hMdpaGhQIBAI2QAAgHPZtty8/PLLqqmp0a9//WtJUmVlpSTJ4/GE7OfxeIKPHU9ubq7cbndwS05ObrHMAADAerYtN88995wuv/xy+Xy+0zrOvHnz5Pf7g1t5eXmEEgIAADuKsTrA8ezbt09vvvmm1q1bFxxLSkqSJFVVVcnr9QbHq6qqlJ6efsJjxcbGKjY2tsWyAgAAe7Hlys2KFSvUrVs3jR8/PjiWmpqqpKQk5efnB8cCgYAKCgqUmZlpRUwAAGBDtlu5aWpq0ooVKzR9+nTFxPz/eIZhKCcnRwsXLlSfPn2Umpqq+fPny+fzaeLEidYFBgAAtmK7cvPmm29q//79uvHGG495bM6cOaqrq1N2drZqamo0evRo5eXlKS4uzoKkAADAjmx9n5uWwH1uAAA48zjiPjcAAADhoNwAAABHodwAAABHodwAAABHodwAAABHodwAABBhFf56fbCnWhX+equjtEm2u88NAABnsrVb92veuhI1mVKUIeVOGqTJI1OsjtWmsHIDAECEVPjrg8VGkppM6Z51n7CC08ooNwAAREhpdV2w2BzVaJoqqz5kTaA2inIDAECEpCbGK8oIHYs2DPVK7GBNoDaKcgMAQIR43S7lThqkaOOHhhNtGFo0KU1et8viZG0LFxQDABBBk0emaEzfriqrPqReiR0oNhag3AAAEGFet4tSYyFOSwEAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEeh3AAAAEexXbn56quv9Ktf/UpdunSRy+XSoEGD9NFHHwUfN01TCxYskNfrlcvlUlZWlnbt2mVhYgAAYCe2KjfffvutRo0apXbt2un111/Xp59+qkcffVRnn312cJ/Fixdr6dKlWr58uQoKChQfH69x48bp8OHDFiYHAAB2YZimaVod4qi5c+fq/fff13vvvXfcx03TlM/n08yZMzVr1ixJkt/vl8fj0cqVKzVlypR/+TsCgYDcbrf8fr8SEhIimh8AALSM5rx/22rl5pVXXtGIESP0i1/8Qt26ddPQoUP17LPPBh8vLS1VZWWlsrKygmNut1sZGRnavHnzcY/Z0NCgQCAQsgEAAOeyVbnZu3evli1bpj59+uiNN97Qbbfdpv/8z//UCy+8IEmqrKyUJHk8npDneTye4GM/lZubK7fbHdySk5Nb9o8AAACWslW5aWpq0rBhw7Ro0SINHTpU2dnZuvnmm7V8+fKwjzlv3jz5/f7gVl5eHsHEAADAbmxVbrxerwYMGBAydt5552n//v2SpKSkJElSVVVVyD5VVVXBx34qNjZWCQkJIRsAAM1V4a/XB3uqVeGvtzoK/gVblZtRo0Zp586dIWNffPGFevbsKUlKTU1VUlKS8vPzg48HAgEVFBQoMzOzVbMCANqOtVv3a9RDb+m6Zws06qG3tHbrfqsj4SRsVW7uuusubdmyRYsWLdLu3bu1evVqPfPMM5oxY4YkyTAM5eTkaOHChXrllVdUUlKiadOmyefzaeLEidaGBwA4UoW/XvPWlajp/z5b3GRK96z7hBUcG4uxOsCPjRw5UuvXr9e8efP0+9//XqmpqXrsscc0derU4D5z5sxRXV2dsrOzVVNTo9GjRysvL09xcXEWJgcAOFVpdV2w2BzVaJoqqz4kr9tlTSiclK3uc9MauM8NAKA5Kvz1GvXQWyEFJ9owtGnuJZSbVnTG3ucGAAC78bpdyp00SNGGIemHYrNoUhrFxsZsdVoKAAA7mjwyRWP6dlVZ9SH1SuxAsbE5yg0AAKfA63ZRas4QnJYCAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACO0qxvBc/MzNTQoUOVnp6u9PR0DR48WHFxcS2VDQAAoNmaVW7Gjx+vjz/+WI8++qj27NkjwzDUp0+fYNk5unXr1q2l8gIAAJyUYZqmGc4TP/zwQ02cOFGjR49Wu3btVFhYqM8//1yGYcjj8ejAgQORzhoRgUBAbrdbfr9fCQkJVscBAACnoDnv381aufmx2267TU899ZSuvvrq4Nhrr72m7OxsTZ8+PdzDAgAAnJawLyj+7LPPlJ6eHjJ2xRVX6Omnn9YHH3xwurkAAADCEna5GTlypF544YVjxgcNGqQPP/zwtEIBAACEK+zTUkuWLNGll16qffv26a677lJaWpqOHDmiRx99VImJiZHMCAAAcMrCLjfDhw9XQUGBZsyYofT0dLVr105NTU2KiYnRc889F8mMAAAApyzscvPNN9+of//+ys/P1/79+1VUVKSoqCgNHz5cXq83khkBAABOWdjlJjExUd27d9eQIUNCtqSkpEjmAwAAaJawy01JSYmKiopUXFysrVu36plnntE333yjuLg4paWlqaCgIJI5AQAATknY5WbgwIEaOHCgpk6dKkkyTVN5eXn6zW9+o7Fjx0YsIAAAQHNE7IszDcPQ5ZdfrhdffFGVlZWROiwAAECzRPxbwS+44AK9/fbbkT4sAADAKQn7tNRZZ52lQYMGaciQIRo8eLCGDBmi/v37a+vWraqtrY1kRgAAgFMWdrn561//qqKiIhUVFenxxx/Xnj17ZJqmDMPQAw88EMmMAAAApyzsbwX/qUOHDqm0tFRdunSx9cfB+VZwAADOPM15/w77mpsLL7xQgUAg+HOHDh00cOBAWxcbAADgfGGXmy1btujw4cPHjAcCAf32t789rVAAAADhana5ueaaa/TQQw/JMAx9/fXXxzxeV1enRx55JCLhAAAAmqvZFxSnpKTo1VdflWmaGjJkiLp06RL86oX09HTt3LmT75YCAEiSKvz1Kq2uU2pivLxul9Vx0EY0u9wsWbJEktS+fXu9//77OnDggAoLC1VUVKT169erqalJixcvjnhQAMCZZe3W/Zq3rkRNphRlSLmTBmnyyBSrY6ENCPvTUt9//71iYmJkGEakM7UoPi0FAC2vwl+vUQ+9paYfvcNEG4Y2zb2EFRyEpcU+LXX99dervr5eklRRUXHGFRsAQOsora4LKTaS1GiaKqs+ZE0gtCnNKjfx8fFqaGiQJPXq1UtdunTRJZdcorvuuksrV65UUVGRvv/++7DD/Nd//ZcMwwjZ+vfvH3zcNE0tWLBAXq9XLpdLWVlZ2rVrV9i/DwDQMlIT4xX1k//+jTYM9UrsYE0gtCnNKjfLly9Xp06dJEmlpaV6/vnndckll2j//v164IEHNHz4cJ111lkaMmRI2IEGDhyoioqK4LZp06bgY4sXL9bSpUu1fPlyFRQUKD4+XuPGjTvuR9IBANbxul3KnTRI0f+3wh9tGFo0KY1TUmgVYX/9Ql1dnSZMmKCrrroqOFZbW6uioiJ9/PHH4QeKiTnujQBN09Rjjz2me++9N/g7//SnP8nj8ejll1/WlClTwv6dAIDImzwyRWP6dlVZ9SH1SuxAsUGrCbvcpKWlKS4uTgMGDAj5KPjgwYN10UUXhR1o165d8vl8iouLU2ZmpnJzc5WSkqLS0lJVVlYqKysruK/b7VZGRoY2b958wnLT0NAQPJUmKeSuygCAluV1uyg1aHVh36F448aNSkhIUPfu3VVbW6tnn31WP/vZz9SlSxf169dP8+fPV01NTbOOmZGRoZUrVyovL0/Lli1TaWmpLrroItXW1qqyslKS5PF4Qp7j8XiCjx1Pbm6u3G53cEtOTm723woAAM4cYX8UfNiwYZo/f76uvvrq4Fh+fr5uueUWTZ8+XW+++ab279+vDz/8UF27dg0rXE1NjXr27KklS5bovPPO06hRo3TgwIGQmwRee+21MgxDa9euPe4xjrdyk5yczEfBAQA4g7TKF2d+/vnnGjhwYMjY2LFj9cc//lHFxcV65513NGLECN1zzz3h/gp16tRJffv21e7du4PX4VRVVYXsU1VVddIv64yNjVVCQkLIBgAAnCvscjN8+HCtWrXqmPG0tDT9/e9/l2EYmj17tt58882ww3333Xfas2ePvF6vUlNTlZSUpPz8/ODjgUBABQUFyszMDPt3AAAAZwm73DzyyCNasmSJrr/+en3++eeSpCNHjuiPf/yjOnfuLEnq2rXrMSstJzNr1ixt3LhRZWVl+uCDD3T11VcrOjpav/zlL2UYhnJycrRw4UK98sorKikp0bRp0+Tz+TRx4sRw/wwAAOAwYX9a6uinlO68804NGDBAsbGx+uc//6mYmBitWLFCklRYWCifz3fKx/zyyy/1y1/+UgcPHlTXrl01evRobdmyJXjNzpw5c1RXV6fs7GzV1NRo9OjRysvLU1xcXLh/BgAAcJiwLyj+sX379qm4uFjR0dEaPnx48BqY9957T1VVVbrmmmtOO2ik8N1SAACceZrz/h2RcnMmodwAAHDmac77d9inpbZu3aq5c+fqH//4h84991ylp6cHt5QUvtIeAABYI+wLiq+//npFR0crOztbqamp2rhxo2644YbgF2oCAABYIeyVm/Lycv3v//6vzjnnnJDxffv2qaio6HRzAQAAhCXscpOZmamvvvrqmHLTs2dP9ezZ87SDAQAAhCPs01J33XWXfv/73+ubb76JZB4AAIDTEvbKzYQJE2QYhvr27aurrrpKmZmZGjp0qAYNGqT27dtHMiMAAMApC7vc7N69W8XFxcFt0aJFKisrU7t27dSvXz99/PHHkcwJAABwSsIuN71791bv3r1DvhU8EAiouLiYYgMAACzDTfwAAIDtNef9O+wLiv/5z3/qwQcfVGZmpoYNG6bp06drw4YN4R4OAAAgIsIuN3PnztXTTz+tsWPHauLEiWpoaNCVV16pG264QW1sMQgAANhI2NfcrF69WmvWrNGYMWOCY6Wlpbryyiv1yCOPaPbs2REJCAAA0Bxhr9zU1dWpR48eIWOpqal64okn9Mwzz5x2MAAAgHCEXW5Gjx6tF1544Zjx1NRUHThw4LRCAQAAhCvs01IPP/ywRo0apW+//Va/+c1v1KdPH33//fd64oknNGDAgEhmBAAAOGVhl5u0tDS98847ys7O1pNPPqn27dursbFRnTp10vr16yOZEQAA4JSFXW4kaejQodq6dat27typHTt2qGPHjsrIyOD+MQAAwDKnVW6O6tevn/r16ydJ2rVrlyZMmKCNGzdG4tAAAADNEvYFxSdy5MgRbdq0KdKHBQAAOCURLzcAAABWavZpqVtvvVXDhw/X0KFDNXjwYLVv374lcgEAAISl2eWmpKREq1atUl1dndq1a6cBAwZo2LBhGj58uIYNG6aoKBaDAACAdZpdbt5//32ZpqmdO3dq+/btwW39+vWqqamRJBmGEemcAAAApySsT0sZhqH+/furf//+uu6664Lje/fu1bZt21RYWBixgAAAAM1hmG3sK7wDgYDcbrf8fj/34wEA4AzRnPdvLpABAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOQrkBAACOYtty89BDD8kwDOXk5ATHTNPUggUL5PV65XK5lJWVpV27dlkXEgAA2I4ty83WrVv13//93xo8eHDI+OLFi7V06VItX75cBQUFio+P17hx43T48GGLkgIAALuxXbn57rvvNHXqVD377LM6++yzg+Omaeqxxx7Tvffeq6uuukqDBw/Wn/70Jx04cEAvv/yydYEBAICt2K7czJgxQ+PHj1dWVlbIeGlpqSorK0PG3W63MjIytHnz5hMer6GhQYFAIGQDAADOFWN1gB9bs2aNtm/frq1btx7zWGVlpSTJ4/GEjHs8nuBjx5Obm6v7778/skEBAIBt2Wblpry8XHfeeadWrVqluLi4iB133rx58vv9wa28vDxixwYAAPZjm3Kzbds2ff311xo2bJhiYmIUExOjjRs3aunSpYqJiQmu2FRVVYU8r6qqSklJSSc8bmxsrBISEkI2AADgXLYpN2PHjlVJSYmKioqC24gRIzR16lQVFRWpd+/eSkpKUn5+fvA5gUBABQUFyszMtDA5AACwE9tcc9OxY0elpaWFjMXHx6tLly7B8ZycHC1cuFB9+vRRamqq5s+fL5/Pp4kTJ1qQGAAA2JFtys2pmDNnjurq6pSdna2amhqNHj1aeXl5Eb1GBwAAnNkM0zRNq0O0pkAgILfbLb/fz/U3AACcIZrz/m2ba24AAAAigXIDAAAchXIDAAAchXIDAAAchXIDAAAchXIDAK2kwl+vD/ZUq8Jfb3UUwNHOqPvcAMCZau3W/Zq3rkRNphRlSLmTBmnyyBSrYwGOxMoNALSwCn99sNhIUpMp3bPuE1ZwgBZCuQGAFlZaXRcsNkc1mqbKqg9ZEwhwOMoNALSw1MR4RRmhY9GGoV6JHawJBDgc5QYAWpjX7VLupEGKNn5oONGGoUWT0uR1uyxOBjgTFxQDQCuYPDJFY/p2VVn1IfVK7ECxAVoQ5QYAWonX7aLUAK2A01IAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRKDcAAMBRbFVuli1bpsGDByshIUEJCQnKzMzU66+/HnzcNE0tWLBAXq9XLpdLWVlZ2rVrl4WJAQCA3diq3PTo0UMPPfSQtm3bpo8++kiXXnqprrrqKu3YsUOStHjxYi1dulTLly9XQUGB4uPjNW7cOB0+fNji5AAAwC4M0zRNq0OcTOfOnfWHP/xBN954o3w+n2bOnKlZs2ZJkvx+vzwej1auXKkpU6ac0vECgYDcbrf8fr8SEhJaMjoAAIiQ5rx/22rl5scaGxu1Zs0a1dXVKTMzU6WlpaqsrFRWVlZwH7fbrYyMDG3evPmEx2loaFAgEAjZAACAc9mu3JSUlOiss85SbGysbr31Vq1fv14DBgxQZWWlJMnj8YTs7/F4go8dT25urtxud3BLTk5u0fwAAMBatis3/fr1U1FRkQoKCnTbbbdp+vTp+vTTT8M+3rx58+T3+4NbeXl5BNMCAAC7ibE6wE+1b99e5557riRp+PDh2rp1qx5//HH99re/lSRVVVXJ6/UG96+qqlJ6evoJjxcbG6vY2NgWzQwAAOzDdis3P9XU1KSGhgalpqYqKSlJ+fn5wccCgYAKCgqUmZlpYUIAAGAntlq5mTdvni6//HKlpKSotrZWq1ev1jvvvKM33nhDhmEoJydHCxcuVJ8+fZSamqr58+fL5/Np4sSJVkcHAAA2Yaty8/XXX2vatGmqqKiQ2+3W4MGD9cYbb+iyyy6TJM2ZM0d1dXXKzs5WTU2NRo8erby8PMXFxVmcHAAA2IXt73MTadznBgCAM48j7nMDoGVU+Ov1wZ5qVfjrrY4CAC3CVqelALSstVv3a966EjWZUpQh5U4apMkjU6yOBQARxcoN0EZU+OuDxUaSmkzpnnWfsIIDwHEoN0AbUVpdFyw2RzWapsqqD1kTCABaCOUGaCNSE+MVZYSORRuGeiV2sCYQALQQyg3QRnjdLuVOGqRo44eGE20YWjQpTV63y+JkABBZXFAMtCGTR6ZoTN+uKqs+pF6JHSg2AByJcgO0MV63i1IDwNE4LQUAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAAByFcgMAABzFVuUmNzdXI0eOVMeOHdWtWzdNnDhRO3fuDNnHNE0tWLBAXq9XLpdLWVlZ2rVrl0WJAQCA3diq3GzcuFEzZszQli1btGHDBn3//ff6t3/7N9XV1QX3Wbx4sZYuXarly5eroKBA8fHxGjdunA4fPmxhcgAAYBeGaZqm1SFO5B//+Ie6deumjRs3asyYMTJNUz6fTzNnztSsWbMkSX6/Xx6PRytXrtSUKVP+5TEDgYDcbrf8fr8SEhJa+k8AAAAR0Jz3b1ut3PyU3++XJHXu3FmSVFpaqsrKSmVlZQX3cbvdysjI0ObNm497jIaGBgUCgZANAAA4l23LTVNTk3JycjRq1CilpaVJkiorKyVJHo8nZF+PxxN87Kdyc3PldruDW3JycssGBwAAlrJtuZkxY4Y++eQTrVmz5rSOM2/ePPn9/uBWXl4eoYQAAMCOYqwOcDx33HGHXn31Vb377rvq0aNHcDwpKUmSVFVVJa/XGxyvqqpSenr6cY8VGxur2NjYFs0LAADsw1YrN6Zp6o477tD69ev11ltvKTU1NeTx1NRUJSUlKT8/PzgWCARUUFCgzMzM1o4LAABsyFYrNzNmzNDq1av1t7/9TR07dgxeR+N2u+VyuWQYhnJycrRw4UL16dNHqampmj9/vnw+nyZOnGhteAAAYAu2KjfLli2TJP3sZz8LGV+xYoV+/etfS5LmzJmjuro6ZWdnq6amRqNHj1ZeXp7i4uJaOS0AALAjW9/npiVwnxsAAM48jrnPDQAAQHNRbgAAgKNQbgAAgKNQbgAAgKNQbgAAgKNQbiKowl+vD/ZUq8Jfb3UUAADaLFvd5+ZMtnbrfs1bV6ImU4oypNxJgzR5ZIrVsQAAaHNYuYmACn99sNhIUpMp3bPuE1ZwAACwAOUmAkqr64LF5qhG01RZ9SFrAgEA0IZRbiIgNTFeUUboWLRhqFdiB2sCAQDQhlFuIsDrdil30iBFGz80nGjD0KJJafK6XRYnAwCg7eGC4giZPDJFY/p2VVn1IfVK7ECxAQDAIpSbCPK6XZQaAAAsxmkpAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKJQbAADgKG3uu6VM05QkBQIBi5MAAIBTdfR9++j7+Mm0uXJTW1srSUpOTrY4CQAAaK7a2lq53e6T7mOYp1KBHKSpqUkHDhxQx44dZRiG1XEUCASUnJys8vJyJSQkWB3HVpibE2NuTo75OTHm5sSYm5Ozen5M01Rtba18Pp+iok5+VU2bW7mJiopSjx49rI5xjISEBP5lOgHm5sSYm5Njfk6MuTkx5ubkrJyff7VicxQXFAMAAEeh3AAAAEeh3FgsNjZW9913n2JjY62OYjvMzYkxNyfH/JwYc3NizM3JnUnz0+YuKAYAAM7Gyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyk0reOqpp9SrVy/FxcUpIyNDH3744Sk97/3331dMTIzS09NbNqCFmjs3DQ0N+t3vfqeePXsqNjZWvXr10vPPP99KaVtXc+dm1apVGjJkiDp06CCv16sbb7xRBw8ebKW0refdd9/VhAkT5PP5ZBiGXn755X/5nHfeeUfDhg1TbGyszj33XK1cubLFc1qhuXOzbt06XXbZZeratasSEhKUmZmpN954o3XCWiCcf3aOcvrrcThzY+fXY8pNC1u7dq3uvvtu3Xfffdq+fbuGDBmicePG6euvvz7p82pqajRt2jSNHTu2lZK2vnDm5tprr1V+fr6ee+457dy5Uy+99JL69evXiqlbR3Pn5v3339e0adP0H//xH9qxY4f+8pe/6MMPP9TNN9/cyslbXl1dnYYMGaKnnnrqlPYvLS3V+PHjdckll6ioqEg5OTm66aabHPkm3ty5effdd3XZZZfptdde07Zt23TJJZdowoQJKiwsbOGk1mju/BzVFl6Pw5kbW78em2hR559/vjljxozgz42NjabP5zNzc3NP+rzJkyeb9957r3nfffeZQ4YMaeGU1mju3Lz++uum2+02Dx482FoRLdPcufnDH/5g9u7dO2Rs6dKlZvfu3Vs0p9UkmevXrz/pPnPmzDEHDhwYMjZ58mRz3LhxLZjMeqcyN8czYMAA8/777498IJtpzvy0hdfjHzuVubH76zErNy3oyJEj2rZtm7KysoJjUVFRysrK0ubNm0/4vBUrVmjv3r267777WiOmJcKZm1deeUUjRozQ4sWL1b17d/Xt21ezZs1SfX19a8VuFeHMTWZmpsrLy/Xaa6/JNE1VVVXpr3/9q6644orWim1bmzdvDplLSRo3btxJ/x1sq5qamlRbW6vOnTtbHcU22sLrcTjs/nrc5r44szVVV1ersbFRHo8nZNzj8ejzzz8/7nN27dqluXPn6r333lNMjHP/7wlnbvbu3atNmzYpLi5O69evV3V1tW6//XYdPHhQK1asaI3YrSKcuRk1apRWrVqlyZMn6/Dhw/rnP/+pCRMmNHv53YkqKyuPO5eBQED19fVyuVwWJbOfRx55RN99952uvfZaq6PYQlt5PQ6H3V+PWbmxkcbGRl133XW6//771bdvX6vj2E5TU5MMw9CqVat0/vnn64orrtCSJUv0wgsv2Oa/Fqzy6aef6s4779SCBQu0bds25eXlqaysTLfeeqvV0XCGWL16te6//379+c9/Vrdu3ayOYzlej0/O7q/HVNEWlJiYqOjoaFVVVYWMV1VVKSkp6Zj9a2tr9dFHH6mwsFB33HGHpB/+ATJNUzExMfr73/+uSy+9tFWyt7Tmzo0keb1ede/ePeQr78877zyZpqkvv/xSffr0adHMrSWcucnNzdWoUaM0e/ZsSdLgwYMVHx+viy66SAsXLpTX623x3HaVlJR03LlMSEhg1eb/rFmzRjfddJP+8pe/HHMKr61qS6/H4bD76zErNy2offv2Gj58uPLz84NjTU1Nys/PV2Zm5jH7JyQkqKSkREVFRcHt1ltvVb9+/VRUVKSMjIzWjN+imjs30g+nXg4cOKDvvvsuOPbFF18oKipKPXr0aPHMrSWcuTl06JCiokL/dY6OjpYkmW386+MyMzND5lKSNmzYcMK5bGteeukl3XDDDXrppZc0fvx4q+PYRlt6PQ6H7V+PrbyauS1Ys2aNGRsba65cudL89NNPzezsbLNTp05mZWWlaZqmOXfuXPP6668/4fOdfHV+c+emtrbW7NGjh3nNNdeYO3bsMDdu3Gj26dPHvOmmm6z6E1pMc+dmxYoVZkxMjPn000+be/bsMTdt2mSOGDHCPP/88636E1pMbW2tWVhYaBYWFpqSzCVLlpiFhYXmvn37TNM8dm727t1rdujQwZw9e7b52WefmU899ZQZHR1t5uXlWfUntJjmzs2qVavMmJgY86mnnjIrKiqCW01NjVV/Qotq7vz8lJNfj5s7N3Z/PabctIInnnjCTElJMdu3b2+ef/755pYtW4KPTZ8+3bz44otP+Fwn/8tkms2fm88++8zMysoyXS6X2aNHD/Puu+82Dx061MqpW0dz52bp0qXmgAEDTJfLZXq9XnPq1Knml19+2cqpW97bb79tSjpmmz59ummax5+bt99+20xPTzfbt29v9u7d21yxYkWr524NzZ2biy+++KT7O004/+z8mJNfj8OZGzu/Hhum2cbXrAEAgKNwzQ0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0AAHAUyg0A27nllls0depUq2MAOEPx3VIAWs0bb7yhf//3f/+X+4wYMUKxsbGKj49vpWQAnIRyA6DV1NfXy+/3B39OS0vT7bffrttvvz041rVrV0VHR1sRD4BDcFoKQKtxuVxKSkpSUlKSGhsbdfDgQV100UXBsaSkJJWXl8swDJWVlQWft3//fl133XU6++yz1blzZ02dOlXffvutJKmsrEyGYeh//ud/NGbMGLlcLo0cOVL79+/Xe++9pwsuuEAdOnTQ2LFjVVNTEzxmZWWlDMPQ448/rqFDhyouLk4DBw7Upk2bQjJ/8sknuuKKK5SQkKCkpCTNnDlTR44caY3pAhAmyg0ASxQWFkqShg0bFjJeXFysTp06qVevXpKk3bt3a/jw4Tr33HO1ZcsWbdiwQbt379bs2bOD+0vSsmXLtGjRIn3wwQeqqqrSr371Kz300EN68skn9fbbb6u4uFgrVqwI/p6ioiJJ0vPPP6/HHntMRUVFSklJ0dSpU9XU1BTMeOGFF2rYsGHavn271qxZo5deekkPP/xwS04NgNMUY3UAAG3T9u3blZycrC5duoSMFxUVafDgwcGfj562uv/++4Njc+bMCZaboqIide7cWWvXrg0e6+KLL9amTZu0Y8cOdejQQZI0cuRIVVZWBo9RXFysdu3a6W9/+1uwSC1cuFAjRozQV199peTkZN188826/vrrtXDhQknSueeeqxtuuEGvvvqq5s+fH/lJARARlBsAlti+ffsxqzbSD6UjPT1dkrRv3z5t2LBBmzZt0qOPPhrcp7GxUcnJycH9r7766pCStH//fk2ePDlYbI6OXXXVVcGfi4qKNGnSpGCxkaSEhITg//7888+1bds2vfjiiyH52rdvr4aGhvD+aACtgtNSACxxonJTVFSkIUOGSPqhuHTu3Fkff/yxioqKgltJSYnefvvt4P4ZGRkhxyguLtYFF1wQ/Pnw4cPauXNn8LhHn3e0RB21efNmJSYmqnv37tqxY4fatWunvn37huzz6aefatCgQaf1twNoWazcAGh11dXVKi8vP6bcBAIBlZWVBUtHu3btVFtbK5/PF7IK89P9hw4dGhwrLS2V3+8PGSspKZFpmsFSUl9fr127dqmxsTG4T1NTkx577DFNnz5dUVFR6tixoxobG/X9998rNjY2eOz169frlVdeidhcAIg8Vm4AtLrt27dLOv7FxNHR0Ro4cKAkKSMjQwkJCZo2bZqKi4u1e/du5eXlKScnJ2T/tLS04DGOXoPTs2fPkLFzzjlHZ511lqQfyo5hGHrxxRe1efNmffbZZ5o8ebJqamp07733Bn93p06dNHfuXO3du1dvvfWWxo8frylTpvzLe/UAsBblBkCrKywslMfjkc/nCxkvLi5W//79gyslnTt31muvvaaDBw9qzJgxGjZsmH73u9+pd+/ewf379eunuLi4kGP8eNXm6NhPT0n1799f99xzj37+859rxIgRamxs1MaNG9WpUydJktvt1ssvv6x3331XAwcO1M0336xp06aFfOIKgD1xEz8Abc6MGTP07bffavXq1VZHAdACWLkB0Ob89OPmAJyFcgOgTTFNUyUlJZQbwME4LQUAAByFlRsAAOAolBsAAOAolBsAAOAolBsAAOAolBsAAOAolBsAAOAolBsAAOAolBsAAOAolBsAAOAolBsAAOAo/w9J0Dw9uWFp3QAAAABJRU5ErkJggg=="
          },
          "metadata": {}
        }
      ]
    },
    {
      "cell_type": "code",
      "source": "pi = math.pi\n",
      "metadata": {
        "trusted": true
      },
      "execution_count": 41,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "def modelo_periodo(L, g):\n    return 2*pi*math.sqrt(L/g)",
      "metadata": {
        "trusted": true
      },
      "execution_count": 45,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "Y = np.log(np.abs(L))\nX = np.log(M1)",
      "metadata": {
        "trusted": true
      },
      "execution_count": 46,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "plt.plot (X, Y)\nplt.xlabel(r'$\\ln(l)$')\nplt.ylabel(r'$\\ln(y)$')",
      "metadata": {
        "trusted": true
      },
      "execution_count": 48,
      "outputs": [
        {
          "execution_count": 48,
          "output_type": "execute_result",
          "data": {
            "text/plain": "Text(0, 0.5, '$\\\\ln(y)$')"
          },
          "metadata": {}
        },
        {
          "output_type": "display_data",
          "data": {
            "text/plain": "<Figure size 640x480 with 1 Axes>",
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGyCAYAAAACgQXWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAABJlUlEQVR4nO3deVzUdf4H8NdwDQjMIAgCgkDiBQjeOlqZiUe5prWtmgaU5mY/W3VrK9laXdoKq7XN7TC7VPIgNY9a77wzRERQFK9CDpHDiwEEBpj5/P4gJlFRrpnvHK/n48Ef853vF16fMOfld77zfcuEEAJEREREFspG6gBEREREhsSyQ0RERBaNZYeIiIgsGssOERERWTSWHSIiIrJoLDtERERk0Vh2iIiIyKLZSR1AajqdDpcuXYKrqytkMpnUcYiIiKgJhBAoKyuDr68vbGzufu7G6svOpUuX4O/vL3UMIiIiaoG8vDz4+fnddR+rLzuurq4A6v5jKRQKidMQERFRU5SWlsLf31//On43Vl926t+6UigULDtERERmpimXoPACZSIiIrJoLDtERERk0Vh2iIiIyKKx7BAREZFFY9khIiIii8ayQ0RERBaNZYeIiIgsGssOERERWTSWHSIiIrJoLDtERERk0Vh2iIiIyKKx7BAREZFFY9khIiIigzlTWIrcqxWSZmDZISIiIoP434lLePyTn/H8ylRUVmsly2GyZWfhwoWQyWSYO3duk/Y/dOgQ7Ozs0Lt3b4PmIiIiorur1eoQv/U0XlydhsoaLTycHVBdq5Msj51kP/kuUlJSsHTpUoSHhzdp/5KSEkRHR2PEiBEoKioycDoiIiJqzLUb1fjLmmM49MtVAMDMYV3wyujusLWRSZbJ5M7slJeXY+rUqfjiiy/Qvn37Jh0zc+ZMTJkyBSqVysDpiIiIqDEn89UY99FPOPTLVbRzsMXHU/pg3iM9JC06gAmWnVmzZmHs2LGIjIxs0v7Lli1DVlYWFixY0KT9NRoNSktLG3wRERFR62w4dhF/XPIz8ksqEeDRDhv/byj+EO4rdSwAJvY2VmJiIo4dO4aUlJQm7X/+/HnMmzcPBw8ehJ1d05YSHx+PuLi41sQkIiKi39RodXh7y2ks/zkbADC8uyc+nNQHynb20ga7icmc2cnLy8OcOXOwatUqODo63nN/rVaLKVOmIC4uDt26dWvyz4mNjYVardZ/5eXltSY2ERGR1bpcpsHUL5P1RWf2w8H4KmaASRUdAJAJIYTUIQBg06ZNePzxx2Fra6vfptVqIZPJYGNjA41G0+C5kpIStG/fvsE2nU4HIQRsbW2xc+dOPPzww/f8uaWlpVAqlVCr1VAoFG27KCIiIguVnleCmd+korC0Ci5yOyyaGIHRod5G+/nNef02mbexRowYgYyMjAbbnn32WfTo0QOvvfZag1IDAAqF4rb9P/30U+zZswfr169HUFCQwTMTERFZo29TcvGPTadQrdXhPk9nfB7VH8FeLlLHapTJlB1XV1eEhYU12Obs7AwPDw/99tjYWOTn5yMhIQE2Nja37e/l5QVHR8fbthMREVHraWq1iPshE6uTcwEAo0I6YtHECLg6mtbbVrcymbLTFAUFBcjNzZU6BhERkdUpKq3CCytTcSy3BDIZ8FJkN8waHgwbiT9W3hQmc82OVHjNDhER0d0dzb6GF1Ydw+UyDVwd7fDfyX0wvIeXpJnM8podIiIiMi1CCKw8nIO4HzJRqxPo3tEVS6P6IbCDs9TRmoVlh4iIiG5TVaPFPzadxLrUiwCAsb188N6T4XCWm191ML/EREREZFCXSioxc2UqTlxUw0YGvDamB/784H2QyUz/+pw7YdkhIiIivaRfr+LF1cdw9UY13NrZ46On+uCBrp5Sx2oVlh0iIiKCEAJfH8rGO1tPQ6sTCPFRYGlUP/i7t5M6Wqux7BAREVm5ymot5m04gc3plwAAE3r7Iv6JcDg52N7jSPPAskNERGTF8q5V4PlvUpFZUApbGxn+/mhPTBsaaLbX59wJyw4REZGVOnj+Mv6yJg0lFTXwcHbAx1P6QtXFQ+pYbY5lh4iIyMoIIbD0QBbe234GOgGE+ynx2dP94OvmJHU0g2DZISIisiI6ncCcb9Pxw/G663P+1M8P/5oQBkd7y7g+505YdoiIiKzIgfOX8cPxS7CzkWHBY6F4elBni7o+505YdoiIiKzIzswiAMDEAf6IGhwgcRrjsJE6ABERERmHTifw429lZ2RIR4nTGA/LDhERkZXIyFejuEwDZwdbDLHAT101hmWHiIjISuz67azOsO6ekNtZ7gXJt2LZISIishK7rPAtLIBlh4iIyCrkXq3A2aIy2NrIMLy7l9RxjIplh4iIyArszCwEAAwIbA+3dg4SpzEulh0iIiIr8OPp+rewvCVOYnwsO0RERBaupKIaKdnXAQCjrOx6HYBlh4iIyOLtOVMMrU6gh7cr/N3bSR3H6Fh2iIiILJy1fgqrHssOERGRBauq0WL/ucsAgMieLDtERERkYZKyrqKiWouOCjl6dVJKHUcSLDtEREQWrP4trMieHWFjY9nTzRvDskNERGShrHXw561YdoiIiCzUiZsGf6qsaPDnrVh2iIiILNSu3+6abG2DP2/FskNERGShfswsBmDdb2EBLDtEREQWyZoHf96KZYeIiMgC1Q/+HBjobnWDP2/FskNERGSBrP2uyTcz2bKzcOFCyGQyzJ07t9F9NmzYgJEjR8LT0xMKhQIqlQo7duwwXkgiIiITdP1GNVKyrwFg2QFMtOykpKRg6dKlCA8Pv+t+Bw4cwMiRI7F161akpqZi+PDhGDduHNLS0oyUlIiIyPTsPVsMnYDVDv68lZ3UAW5VXl6OqVOn4osvvsBbb711130//PDDBo/feecdbN68GT/88AP69Olzx2M0Gg00Go3+cWlpaaszExERmZKD568AAEb0tO4Lk+uZ3JmdWbNmYezYsYiMjGz2sTqdDmVlZXB3d290n/j4eCiVSv2Xv79/a+ISERGZnPq3sAYFWe+NBG9mUmUnMTERx44dQ3x8fIuO//e//43y8nJMnDix0X1iY2OhVqv1X3l5eS2NS0REZHIK1VW4eL0SNjKgT2c3qeOYBJN5GysvLw9z5szBrl274Ojo2OzjV69ejbi4OGzevBleXo2ftpPL5ZDL5a2JSkREZLKO5tSd1enpo4Cro73EaUyDyZSd1NRUFBcXo2/fvvptWq0WBw4cwMcffwyNRgNb2zvf6joxMRHPPfcc1q1b16K3v4iIiCzF0ezrAID+Ae0lTmI6TKbsjBgxAhkZGQ22Pfvss+jRowdee+21RovOmjVrMG3aNCQmJmLs2LHGiEpERGSy6s/s9A9s/PpVa2MyZcfV1RVhYWENtjk7O8PDw0O/PTY2Fvn5+UhISABQ99ZVTEwMFi9ejEGDBqGwsO5ukU5OTlAqlcZdABERkcTKNbXIvFT3KeP+gTyzU8+kLlC+l4KCAuTm5uoff/7556itrcWsWbPg4+Oj/5ozZ46EKYmIiKSRnlsCnQA6uTnBR+kkdRyTYTJndu5k3759DR4vX778rs8TERFZs/qPnA/gWZ0GzOrMDhERETWu/nqdfrxepwGWHSIiIgtQq9UhLbcEAM/s3Iplh4iIyAKcLihDRbUWro526OblKnUck8KyQ0REZAH0b2EFtIeNjUziNKaFZYeIiMgC1N9McACv17kNyw4REZGZE0LoP4nVj3dOvg3LDhERkZm7eL0SxWUa2NvKEOHnJnUck8OyQ0REZObqz+qEdVLCyeHO45WsGcsOERGRmUvh8M+7YtkhIiIyc6kc/nlXLDtERERmrKSiGueKygHwzE5jWHaIiIjMWGpO3VtY93VwhoeLXOI0pollh4iIyIwd/a3s9OeIiEax7BAREZmxo799Eqt/AK/XaQzLDhERkZnS1Gpx/KIaAM/s3A3LDhERkZk6ma9Gda0OHs4OCOrgLHUck8WyQ0REZKbq76/TL6A9ZDIO/2wMyw4REZEZEkLg0C9XAHD4572w7BAREZkZIQTe/F8mDp6vKzv3d+0gcSLTZid1ACIiImo6nU7g9U0nseZILgDgXxPC0NNHIXEq08ayQ0REZCZqtTq8uv4ENqTlQyYD3v1jOCb295c6lslj2SEiIjIDNVod5iamY0tGAWxtZPhgYgTG9+4kdSyzwLJDRERk4qpqtHhx9TH8eLoY9rYyfPRUX4wJ85Y6ltlg2SEiIjJhldVa/Pmbozh4/grkdjb4LKofhnf3kjqWWWHZISIiMlHlmlpMX56C5AvX4GRvi69i+mNIMD951VwsO0RERCZIXVmDZ5YdQVpuCVzkdlj+7AD05/10WoRlh4iIyMRcu1GNqK+ScepSKZRO9kiYNhAR/m5SxzJbLDtEREQmpLisClFfHsHZojJ4ODvgm+mDEOLL++i0BssOERGRiShQV2LqF8nIunIDXq5yrJ4xCMFerlLHMnssO0RERCYg71oFpnx5GHnXKtHJzQmrnhuEQE4ybxMsO0RERBLLulyOqV8mo0BdhQCPdlj13CD4tW8ndSyLYbKDQBcuXAiZTIa5c+fedb99+/ahb9++kMvlCA4OxvLly42Sj4iIqC2cLSzDxKWHUaCuQhdPZ6x9XsWi08ZMsuykpKRg6dKlCA8Pv+t+Fy5cwNixYzF8+HCkp6dj7ty5eO6557Bjxw4jJSUiImq5k/lqTP48CVfKNejh7Ypvn1eho8JR6lgWx+TKTnl5OaZOnYovvvgC7du3v+u+n332GYKCgrBo0SL07NkTL774Ip588kn85z//MVJaIiKiljmWex1PfXEY1ytqEO6nROKfB6ODi1zqWBbJ5MrOrFmzMHbsWERGRt5z36SkpNv2Gz16NJKSkho9RqPRoLS0tMEXERGRMSVnXUXUl8koq6pF/4D2WPncILi1c5A6lsUyqQuUExMTcezYMaSkpDRp/8LCQnTs2LHBto4dO6K0tBSVlZVwcnK67Zj4+HjExcW1SV4iIqLmOnj+MmYkHEVVjQ5Dunjgi+j+cJab1MuxxTGZMzt5eXmYM2cOVq1aBUdHw71fGRsbC7Varf/Ky8sz2M8iIiK62e7TRZi+oq7oPNTdE18/M4BFxwhM5r9wamoqiouL0bdvX/02rVaLAwcO4OOPP4ZGo4GtrW2DY7y9vVFUVNRgW1FRERQKxR3P6gCAXC6HXM73RImIyLi2ZhRg9po01OoERod2xH+f6gO5ne29D6RWM5myM2LECGRkZDTY9uyzz6JHjx547bXXbis6AKBSqbB169YG23bt2gWVSmXQrERERM2xMe0iXl57HDoBPBbhi0UTI2BvazJvrlg8kyk7rq6uCAsLa7DN2dkZHh4e+u2xsbHIz89HQkICAGDmzJn4+OOP8eqrr2LatGnYs2cP1q5diy1bthg9PxER0Z0kHslF7MYMCAH8qZ8fFv4xHLY2MqljWRWzqpUFBQXIzc3VPw4KCsKWLVuwa9cuREREYNGiRfjyyy8xevRoCVMSERHVWX7oAuZtqCs6UYMD8C6LjiRkQgghdQgplZaWQqlUQq1WQ6HgVFkiImobn+3/FQu3nQEAzHggCH9/tCdkMhadttKc12+TeRuLiIjIEggh8OGP57F493kAwOyHg/HXkd1YdCTEskNERNRGhBBYuP0Mlu7PAgC8Mro7Zg0PljgVsewQERG1AZ1OIO6HU1iRlAMAmP+HEEy7P0jiVASw7BAREbWaVifw+sYMJKbkQSYD3p7QC1MGdZY6Fv2GZYeIiKgVarU6/G3dcWxKvwQbGfD+kxH4Yz8/qWPRTVh2iIiIWqi6Voc5iWnYdrIQdjYyLJ7cB2PDfaSORbdg2SEiImoBIQTmfltXdBxsbfDJ1L4YGdLx3geS0bHsEBERtcD3xy9ha0Zd0fkipj+GdfOUOhI1wqzuoExERGQKSiqq8a//ZQIAXnw4mEXHxLHsEBERNdPCbWdwpbwawV4umDmsi9Rx6B5YdoiIiJohOesqElPyAADxT/SCgx1fSk0df0NERERNpKnVInZjBgDgqYGdMSDQXeJE1BQsO0RERE306d5fkXX5Bjq4yDFvTA+p41ATsewQERE1wS/FZViy71cAwD8fC4Gynb3EiaipWHaIiIjuQacT+PuGk6jW6jC8uyfG9uKNA80Jyw4REdE9rD2ahyPZ1+Bkb4t/TQiDTCaTOhI1A8sOERHRXRSXVeGdracBAC+P6ga/9u0kTkTNxbJDRER0F//632mUVtUirJMCzwwJlDoOtQDLDhERUSP2ni3GD8frppkvfCIcdrZ82TRH/K0RERHdQUV1Ld7YeBIAMG1oEMI6KSVORC3FskNERHQHH/54Hvkllejk5oS/juwmdRxqBZYdIiKiW5zMV+Orny4AAN6aEAZnuZ3Eiag1WHaIiIhuotUJxG7IgFYnMDbcB8N7eEkdiVqJZYeIiOgmy3/ORka+Gq6OdlgwLkTqONQGWHaIiIh+k19SiUU7zwIAYh/pCS9XR4kTUVtg2SEiIgIghMD8TSdRUa1F/4D2mDzAX+pI1EZYdoiIiABsO1mI3WeKYW8rQ/wTvWBjw5EQloJlh4iIrJ66sgYLvj8FAHhhWBd07egqcSJqSyw7RERk9d7bfgaXyzS4r4Mz/m94sNRxqI2x7BARkVU7mn0Nq5JzAQBvP94Ljva2EieitsayQ0REVqu6VofYDRkAgIn9/aDq4iFxIjIEkyo7S5YsQXh4OBQKBRQKBVQqFbZt23bXY1atWoWIiAi0a9cOPj4+mDZtGq5evWqkxEREZM4+P/ArzheXw8PZAX9/tKfUcchATKrs+Pn5YeHChUhNTcXRo0fx8MMPY/z48Th16tQd9z906BCio6Mxffp0nDp1CuvWrcORI0cwY8YMIycnIiJzk3W5HP/d8wsAYP64ELi1c5A4ERmKSQ37GDduXIPHb7/9NpYsWYLDhw8jNDT0tv2TkpIQGBiI2bNnAwCCgoLw/PPP49133zVKXiIiMk9CCLy+8SSqa3V4oGsHPBbhK3UkMiCTOrNzM61Wi8TERNy4cQMqleqO+6hUKuTl5WHr1q0QQqCoqAjr16/Ho48+2uj31Wg0KC0tbfBFRETWZX3qRSRlXYWjvQ3entALMhnvqWPJTK7sZGRkwMXFBXK5HDNnzsTGjRsREnLn2SRDhw7FqlWrMGnSJDg4OMDb2xtKpRKffPJJo98/Pj4eSqVS/+XvzztkEhFZk6vlGry99TQAYG5kN3T2aCdxIjI0kys73bt3R3p6OpKTk/HCCy8gJiYGmZmZd9w3MzMTc+bMwfz585Gamort27cjOzsbM2fObPT7x8bGQq1W67/y8vIMtRQiIjJBb205jZKKGvTwdsX0+4OkjkNGIBNCCKlD3E1kZCS6dOmCpUuX3vZcVFQUqqqqsG7dOv22n376CQ888AAuXboEHx+fe37/0tJSKJVKqNVqKBSKNs1ORESm5eD5y4j66ghkMmDj/w1Fb383qSNRCzXn9dvkzuzcSqfTQaPR3PG5iooK2Ng0XIKtbd3NoEy8wxERkZFVVmvx+saTAIAYVSCLjhUxqU9jxcbG4pFHHkHnzp1RVlaG1atXY9++fdixY4f++fz8fCQkJACo+/TWjBkzsGTJEowePRoFBQWYO3cuBg4cCF9fXllPRES/+++e88i9VgEfpSP+Nrq71HHIiEyq7BQXFyM6OhoFBQVQKpUIDw/Hjh07MHLkSABAQUEBcnNz9fs/88wzKCsrw8cff4yXX34Zbm5uePjhh/nRcyIiauB0QSm+OJAFAIh7LBQucpN6+SMDM/lrdgyN1+wQEVk2rU7gj0t+RnpeCcaEeuOzqH5SR6I2YFHX7BAREbXGquQcpOeVwEVuh38+dvsNasnysewQEZHFKlRX4b3tZwEAr43pDm+lo8SJSAosO0REZLEWfH8S5Zpa9OnshqmDAqSOQxJh2SEiIou041Qhdpwqgp2NDPFP9IKNDUdCWCuWHSIisjhlVTVYsPkUAODPD96HHt78AIo1Y9khIiKL8+8dZ1FYWoUAj3aYPaKr1HFIYiw7RERkUdJyryPhcA4A4O0JveBobytxIpIayw4REVmMGq0OsRsyIATwRJ9OuL9rB6kjkQlg2SEiIovx5cELOFNYhvbt7PH62J5SxyETwbJDREQWIefqDSzefQ4A8PrYEHi4yCVORKaCZYeIiMyeEAJvbDqJqhodhnTxwB/7dpI6EpkQlh0iIjJ7m9Mv4eD5K3Cws8Hbj/eCTMZ76tDvWHaIiMisXb9RjTf/lwkAmDOiK4I6OEuciEwNyw4REZm1d7aexrUb1ejW0QUzHrhP6jhkglh2iIjIbP386xWsS70ImQyIfyIcDnZ8WaPb8U8FERGZpaoaLV7feBIAMHVQZ/QLaC9xIjJVdq05uKamBoWFhaioqICnpyfc3d3bKhcREdFdfbL3F1y4cgNernK8OqaH1HHIhDX7zE5ZWRmWLFmCYcOGQaFQIDAwED179oSnpycCAgIwY8YMpKSkGCIrERERAOB8URk+2/8rACDusVAoHO0lTkSmrFll54MPPkBgYCCWLVuGyMhIbNq0Cenp6Th37hySkpKwYMEC1NbWYtSoURgzZgzOnz9vqNxERGSldDqB2A0ZqNEKRPb0wpgwb6kjkYlr1ttYKSkpOHDgAEJDQ+/4/MCBAzFt2jR89tlnWLZsGQ4ePIiuXTltloiI2s6alFwczbkOZwdbvDk+jPfUoXuSCSFESw7My8uDv79/W+cxutLSUiiVSqjVaigUCqnjEBHRXRSXVmHEB/tRVlWL+X8IwbT7g6SORBJpzut3iz+N1aNHD8yfPx8VFRUt/RZERETNEvdDJsqqahHup0TMkECp45CZaHHZ2bVrF3bs2IGuXbti+fLlbRiJiIjodrtPF2FLRgFsbWSIf6IXbG349hU1TYvLzpAhQ5CcnIz4+Hj84x//QL9+/XDw4MG2zEZERAQAuKGpxfzNpwAAz90fhFBfpcSJyJy0+qaC0dHROHv2LMaOHYtHHnkETz75JC5cuNAW2YiIiAAAH+w6h/ySSvi1d8KcSH7whZqnze6gPGrUKDz33HPYuHEjQkJC8Oqrr6K8vLytvj0REVmpjItqLDtU94/otyaEoZ1Dq+6HS1aoxX9iPvvsM6SkpCAlJQWnT5+GjY0NwsLCMHPmTERERCAxMREhISHYsGED+vfv35aZiYjIStRqdZi34QR0AngswhcPdfeSOhKZoRZ/9Nzf3x+DBg3C4MGDMXjwYPTr1w9OTk4N9nnnnXewevVqnDx5sk3CGgI/ek5EZLq+OJCFt7eehsLRDrtffgiernKpI5GJaM7rd4vLTlMUFRXB19cXWq3WUD+i1Vh2iIhMU961Coz6zwFU1mjx7h97YdKAzlJHIhNilPvsNIWXlxf27NljyB9BREQWSAiBf2w+icoaLQYGuWNif/O/iS1Jp1llJzc3t1nf/NKlSxg2bFizjiEiIvrfiQLsO3sZDrY2eOfxXhwJQa3SrLIzYMAAPP/883edaq5Wq/HFF18gLCwM3333XasDEhGRdVFX1CDuh0wAwP8N74JgLxeJE5G5a1bZyczMhLOzM0aOHAlvb2+MHTsWM2bMwF/+8hc8/fTT6Nu3L7y8vPD111/jvffew+zZs5sVZsmSJQgPD4dCoYBCoYBKpcK2bdvueoxGo8Hrr7+OgIAAyOVyBAYG4uuvv27WzyUiItNQVlWDl9cdx5VyDbp4OuOFh7pIHYksQIsuUK6srMSWLVvw008/IScnB5WVlejQoQP69OmD0aNHIywsrEVhfvjhB9ja2qJr164QQmDFihV4//33kZaW1uik9fHjx6OoqAhvvfUWgoODUVBQAJ1Oh6FDhzbpZ/ICZSIi05Cacw1zv01H3rVK2NrIsGbGYAwMcpc6Fpkok/k0Vltwd3fH+++/j+nTp9/23Pbt2zF58mRkZWXB3b1p/0NoNBpoNBr949LSUvj7+7PsEBFJpFarw3/3/IKP95yHTgB+7Z3w4aTe6B/IokONM5lPY7WGVqtFYmIibty4AZVKdcd9vv/+e/Tv3x/vvfceOnXqhG7duuFvf/sbKisrG/2+8fHxUCqV+i9/f17hT0QklZyrN/CnpUn47+66ovNEn07YOucBFh1qU6265/bu3buxe/duFBcXQ6fTNXiupdfNZGRkQKVSoaqqCi4uLvrxE3eSlZWFn376CY6Ojti4cSOuXLmC//u//8PVq1exbNmyOx4TGxuLl156Sf+4/swOEREZjxAC3x3Lx4LNJ3GjWgtXRzu8NSEM43t3kjoaWaAWl524uDi8+eab6N+/P3x8fNrsY4Hdu3dHeno61Go11q9fj5iYGOzfv/+OhUen00Emk2HVqlVQKusm4H7wwQd48skn8emnn952R2cAkMvlkMt5B04iIqmoK2rw940Z2JJRAAAYGOSODyZGwK99O4mTkaVq1Wys5cuXIyoqqi3zwMHBAcHBwQCAfv36ISUlBYsXL8bSpUtv29fHxwedOnXSFx0A6NmzJ4QQuHjxIrp25WRcIiJT8vOvV/Dy2uMoUFfBzkaGv47shpnDusDWhvfRIcNpcdmprq7GkCFD2jLLHel0ugYXFN9s6NChWLduHcrLy+HiUncfhnPnzsHGxgZ+fn4Gz0ZERE1TXavDB7vOYemBXyEEENTBGYsn90a4n5vU0cgKtPgC5eeeew6rV69uyyyIjY3FgQMHkJ2djYyMDMTGxmLfvn2YOnWq/vno6Gj9/lOmTIGHhweeffZZZGZm4sCBA3jllVcwbdq0O76FRURExvdLcTmeWHIIn+2vKzpPDfTH//5yP4sOGU2Lz+xUVVXh888/x48//ojw8HDY29s3eP6DDz5o9vcsLi5GdHQ0CgoKoFQqER4ejh07dmDkyJEAgIKCggYjK1xcXLBr1y785S9/Qf/+/eHh4YGJEyfirbfeaumyiIiojQghsCo5F29tyURVjQ5u7eyx8IlwjAnzljoaWZkW32dn+PDhd31+7969LQpkbLypIBFR27tarsFr353Aj6eLAQD3B3fAookR6KhwlDgZWYrmvH63+MyOuZQZIiIyrn1ni/HK+hO4XKaBg60NXh3THdOGBsGGFyGTRJpVdm6+P83dyGQyLFq0qEWBiIjIPFXVaLFw2xks/zkbANDVywWLJ/dBiC/PmpO0mlV20tLSmrRfW91zh4iIzMOZwlLMWZOOs0VlAIAYVQBiH+0JR3tbiZMRNbPs8K0rIiK6mU4nsPznbCzcfgbVtTp0cHHA+09GYHgPL6mjEem1alwEERFZr+LSKry87jgOnr8CABjRwwvvPhmODi68Sz2ZFpYdIiJqtp2nCvHadydwvaIGcjsbvPGHEDw9qDMvYyCTxLJDRERNVlFdi7e2nMbq5Lp7noX4KPDfp3oj2MtV4mREjWPZISKiJsm4qMacxDRkXbkBAHj+wfvw0qhukNvxImQybSw7RER0V1qdwOcHsrBo51nU6gS8FY5YNDECQ4M7SB2NqElYdoiIqFGXSirx12/TkXzhGgDgkTBvvPN4L7R3dpA4GVHTsewQEdEd/e/EJfx9QwZKq2rRzsEW/xwXij/19+NFyGR2WHaIiKiBsqoa/PP7THx37CIAIMLfDYsn9UZgB2eJkxG1DMsOERHppeZcx9xv05B3rRI2MuDF4cH4y4iusLe1kToaUYux7BAREWq1Ony89xd8tOcXaHUCndyc8OHk3hgQ6C51NKJWY9khIrJyuVcrMPfbNBzLLQEATOjtizcnhEHhaC9tMKI2wrJDRGSlhBDYcCwfC74/hXJNLVzldnjr8TCM791J6mhEbYplh4jICqkravD6pgz870QBAGBAYHt8MLE3/N3bSZyMqO2x7BARWZmkX6/i5bXpuKSugp2NDH8d2Q0zh3WBrQ0/Uk6WiWWHiMhKVNfq8J8fz+Gz/b9CCCDQox0+nNwHvf3dpI5GZFAsO0REVuDXy+WYk5iGk/mlAIBJ/f0xf1wInOV8GSDLxz/lREQWTAiBNUfy8K//ZaKyRgu3dvZY+EQvjAnzkToakdGw7BARWair5RrM25CBXZlFAIChwR5Y9Kfe8FY6SpyMyLhYdoiILNCBc5fx8rrjuFymgYOtDV4Z3R3T7w+CDS9CJivEskNEZEGqarR4b/tZfH3oAgAg2MsFiyf3RqivUuJkRNJh2SEishBnC8swJzENZwrLAADRqgDEPtITTg62EicjkhbLDhGRmRNCYPnP2YjfdgbVtTp0cHHAe0+G4+EeHaWORmQSWHaIiMxYcVkVXll3AvvPXQYADO/uifeejICnq1ziZESmg2WHiMhM7coswmvfncC1G9WQ29ng9bE9ETU4ADIZL0ImuhnLDhGRmams1uKtLZlYlZwLAOjpo8B/J/dG146uEicjMk0sO0REZuRkvhqzE9OQdfkGAGDGA0H42+jukNvxImSixrDsEBGZAZ1O4PODWVi08yxqtAIdFXIs+lNv3N+1g9TRiEyejdQBbrZkyRKEh4dDoVBAoVBApVJh27ZtTTr20KFDsLOzQ+/evQ0bkojIyC6VVGLql8lYuO0MarQCo0M7YvucB1l0iJrIpM7s+Pn5YeHChejatSuEEFixYgXGjx+PtLQ0hIaGNnpcSUkJoqOjMWLECBQVFRkxMRGRYW05UYDYDSdQWlULJ3tb/POxEEzs78+LkImaQSaEEFKHuBt3d3e8//77mD59eqP7TJ48GV27doWtrS02bdqE9PT0Jn//0tJSKJVKqNVqKBSKNkhMRNR65Zpa/PP7U1ifehEAEOGnxIeT+yCog7PEyYhMQ3Nev03qzM7NtFot1q1bhxs3bkClUjW637Jly5CVlYWVK1firbfeuuf31Wg00Gg0+selpaVtkpeIqK0cy72OuYnpyL1WAZkMmPVQMOZEdoW9rUldeUBkNkyu7GRkZEClUqGqqgouLi7YuHEjQkJC7rjv+fPnMW/ePBw8eBB2dk1bSnx8POLi4toyMhFRm6jV6vDJ3l/x3z3nodUJdHJzwn8m9cbAIHepoxGZNZP7Z0L37t2Rnp6O5ORkvPDCC4iJiUFmZuZt+2m1WkyZMgVxcXHo1q1bk79/bGws1Gq1/isvL68t4xMRtUjetQpM+vww/vPjOWh1AuN7+2LrnAdYdIjagMlfsxMZGYkuXbpg6dKlDbaXlJSgffv2sLX9/d4SOp0OQgjY2tpi586dePjhh+/5/XnNDhFJSQiBjWn5mL/5FMo1tXCV2+FfE8IwoU8nqaMRmTSLuGannk6na3CNTT2FQoGMjIwG2z799FPs2bMH69evR1BQkLEiEhG1iLqyBm9sOokfjl8CAPQPaI//TOoNf/d2EicjsiwmVXZiY2PxyCOPoHPnzigrK8Pq1auxb98+7NixQ/98fn4+EhISYGNjg7CwsAbHe3l5wdHR8bbtRESmJjnrKl5aexz5JZWwtZFh7oiueOGhLrDjRchEbc6kyk5xcTGio6NRUFAApVKJ8PBw7NixAyNHjgQAFBQUIDc3V+KUREQtV12rw4c/nsOS/b9CCCDAox0+nNQbfTq3lzoakcUy+Wt2DI3X7BCRsWRdLsfcb9Nx4qIaADCxvx/mjwuFi9yk/t1JZBYs6podIiJzJ4RAYkoe3vwhE5U1Wiid7BH/RC882stH6mhEVoFlh4jIgK7dqMa8705gZ2bdKJshXTywaGIEfJROEicjsh4sO0REBnLw/GW8vPY4iss0sLeV4ZXR3fHc/ffBxoZzrYiMiWWHiKiNVdVo8f6Os/jqpwsAgC6ezlg8uQ/COiklTkZknVh2iIja0LmiMsxek4YzhWUAgKjBAfj7oz3h5GB7jyOJyFBYdoiI2oAQAit+zsY7286gulYHD2cHvPdkOEb07Ch1NCKrx7JDRNQKWp3AnjPF+PJgFpIvXAMAPNTdE+8/GQFPV7nE6YgIYNkhImqRazeq8W1KHlYezkF+SSUAwMHOBq8/2hPRqgDIZLwImchUsOwQETXD8bwSJCTl4IcTl1BdqwMAuLWzx6T+/nh6cADnWhGZIJYdIqJ7qKrRYsuJAiQczsHxvBL99l6dlIhWBWBchC8c7XkBMpGpYtkhImrExesVWJWci29T8nDtRjUAwMHWBmPDfRCtCkBvfze+XUVkBlh2iIhuotMJHPr1Clb8nIM9Z4qg+216oK/SEVMHB2DSAH90cOGFx0TmhGWHiAhAaVUNvku9iG+ScpB15YZ++/3BHRClCsCIHl6ws7WRMCERtRTLDhFZtTOFpUhIysGmtHxUVGsBAC5yOzzZzw9PDw5AsJeLxAmJqLVYdojI6tRoddhxqhAJSTk48tu9cQCgW0cXRKkC8XifTnCR869HIkvB/5uJyGoUl1Zh9ZFcrDmSi6JSDQDA1kaG0aEdETU4EIPvc+cFx0QWiGWHiCyaEAIp2deRkJSN7ScLUfvbFccdXOSYMtAfTw3qDB+lk8QpiciQWHaIyCJVVNdiU9olJCRl64dyAkD/gPaIUgXgkTAfONjxgmMia8CyQ0QW5cKVG/gmKQfrUvNQVlULAHC0t8GE3p0QpQpAqK9S4oREZGwsO0Rk9rQ6gb1nipFwOAcHzl3Wbw/waIeowQH4Uz9/KNvZS5iQiKTEskNEZuv6jWp8e7RuGOfF63XDOGUyYHh3L0SrAvBgV0/Y2PCCYyJrx7JDRGbnxMW6YZzfH799GOfUQQHo7MFhnET0O5YdIjILVTVabM0owIqkhsM4wzopEK0KxGMcxklEjWDZISKTdrdhnFGqAPThME4iugeWHSIyOUIIHPrlKlYkZWP3aQ7jJKLWYdkhIpOhH8Z5OAdZl38fxjk02ANRgwMR2ZPDOImo+Vh2iEhyZwvLkJCUjY23DOP8Y9+6e+MEe7lKnJCIzBnLDhFJokarw85TRUhIykbyTcM4u3q5IFoVgMf7+nEYJxG1Cf5NQkRGVVxahTVH8rD6SE6DYZyjQjoiWsVhnETU9lh2iMjghBA4mnMdCUk52JZR0GAY51MD/TGFwziJyIBYdojIYCqqa7E5/RJW/NxwGGe/gPaI5jBOIjISlh0ianMXrtzAysM5WHuUwziJSHom9U+qJUuWIDw8HAqFAgqFAiqVCtu2bWt0/w0bNmDkyJHw9PTU779jxw4jJiaielqdwO7TRYj++giG/3sfvvrpAsqqahHg0Q5vjO2J5NhILPxjOIsOERmdSZ3Z8fPzw8KFC9G1a1cIIbBixQqMHz8eaWlpCA0NvW3/AwcOYOTIkXjnnXfg5uaGZcuWYdy4cUhOTkafPn0kWAGR9bl+oxprj+bhmzsM44xSBWAYh3ESkcRkQgghdYi7cXd3x/vvv4/p06c3af/Q0FBMmjQJ8+fPv+PzGo0GGo1G/7i0tBT+/v5Qq9VQKBRtkpnIGtQP4/zh+CVofhvGqXSyx6QB/niawziJyMBKS0uhVCqb9PptUmd2bqbVarFu3TrcuHEDKpWqScfodDqUlZXB3d290X3i4+MRFxfXVjGJrIqmVostJwqQkJSD9FuHcQ4OxLgIXzg5cBgnEZkWkzuzk5GRAZVKhaqqKri4uGD16tV49NFHm3Tse++9h4ULF+LMmTPw8vK64z48s0PUfPkllVh1OAffpuTh6m/DOO1tZRjbywfRQwI5jJOIjM6sz+x0794d6enpUKvVWL9+PWJiYrB//36EhITc9bjVq1cjLi4OmzdvbrToAIBcLodczgGCRPdSP4wzISkbP940jNNH6Yipgzpj0oDO8HTl/0tEZPpM7szOrSIjI9GlSxcsXbq00X0SExMxbdo0rFu3DmPHjm3W929OMySyBmU3DeP89aZhnEO6eCBaFYDInh05jJOIJGfWZ3ZupdPpGrztdKs1a9Zg2rRpSExMbHbRIaLfnSuqG8a54djvwzidHWzxx35+iBocgK4dOYyTiMyTSZWd2NhYPPLII+jcuTPKysqwevVq7Nu3T3/vnNjYWOTn5yMhIQFA3VtXMTExWLx4MQYNGoTCwkIAgJOTE5RK3suD6F5qtDrsyizCip8bDuMM9nJBDIdxEpGFMKm/xYqLixEdHY2CggIolUqEh4djx44dGDlyJACgoKAAubm5+v0///xz1NbWYtasWZg1a5Z+e0xMDJYvX27s+ERmo7isColH8rAq+fZhnFGqAKju8+AFx0RkMUz+mh1D4zU7ZC2EEEjNuY4VSTnYfrIANdr6YZwOeGpgZw7jJCKzYlHX7BBR69QP40xIysHpglL99vphnGPCvCG3471xiMhysewQWajsKzfwzeEcrDuah9KbhnGOj6gbxhnWide1EZF1YNkhsiBancD+c8VY8XMO9p+7rN/e2b0dogYH4E/9/eDWzkHChERExseyQ2QB6odxrkzOQd6134dxPtTNE9GqQAzrxmGcRGS9WHaIzFjGRTUSkrLx/S3DOCf298PTgwMQ4OEscUIiIumx7BCZGU2tFlsz6oZxpuWW6LeH+ioQo+IwTiKiW7HsEJmJ/JJKrE7OQeKRhsM4H+3lg2hVIPp25jBOIqI7YdkhMmFCCPz861Ws+JnDOImIWoplh8gEcRgnEVHbYdkhMiH1wzg3HsvHDQ7jJCJqEyw7RBKrrR/GmZSNw1kNh3FGqwLweJ9OcHW0lzAhEZF5Y9khkkj9MM7VybkoLK0CANjIgFEh3ohWBUDVhcM4iYjaAssOkRHVD+NMSMrBtluGcU4eUDeM09eNwziJiNoSyw6REVRWa7E5PR8JSTnIvGkYZ9/ObohWBeKRXhzGSURkKCw7RAaUfeUGVh7OwdqbhnHK7WwwvrcvolWBHMZJRGQELDtEbax+GGdCUg72nW04jPPpwZ0xsb8/h3ESERkRyw5RGympqBvG+c3hhsM4h3XzRAyHcRIRSYZlh6iVTuarseLnhsM4FY52mDTAn8M4iYhMAMsOUQtoarXYllGIFUnZDYZxhvgoEDMkAI9FdOIwTiIiE8GyQ9QMl0oqsarRYZwB6Nu5Pe+NQ0RkYlh2iO6hfhhnQlI2dmX+PozTW1E3jHPyQA7jJCIyZSw7RI0oq6rBhmP5+OZwDn4pLtdvV91XN4xzZAiHcRIRmQOWHaJbnC8qQ0JSDjYcu9hgGOcTff0QpQpANw7jJCIyKyw7RPh9GGdCUg6Ssq7qt3fxdEa0KhBP9OUwTiIic8WyQ1btcpkGiUdyseqWYZwjQzoiWhWIIRzGSURk9lh2yOoIIXAst24Y59aM34dxejg7YPJAf0wZFIBOHMZJRGQxWHbIalRWa/H98Xys+LnhMM4+nd0Qw2GcREQWi2WHLF7O1Rv4JikH61IvQl1ZA6BuGOdjEXXDOHv5cRgnEZElY9khi6TTCew/dxkrkrKx/9xliN/ujePv7oSowQH4Uz9/tHfmME4iImvAskMWpX4Y58rDuci9VqHf/lB3T0SrAjCsmxdsOYyTiMiqsOyQRTiZr0ZCUjY2pzccxjmxf90wzsAOHMZJRGStTOr2r0uWLEF4eDgUCgUUCgVUKhW2bdt212P27duHvn37Qi6XIzg4GMuXLzdOWJKcplaLTWn5eOLTQ/jDRz9h7dGL0NTqEOKjwMIneiH575F44w8hLDpERFbOpM7s+Pn5YeHChejatSuEEFixYgXGjx+PtLQ0hIaG3rb/hQsXMHbsWMycOROrVq3C7t278dxzz8HHxwejR4+WYAVkDJdKKrE6OReJKbm4Uv77MM5HwuqGcfYL4DBOIiL6nUyI+ks3TZO7uzvef/99TJ8+/bbnXnvtNWzZsgUnT57Ub5s8eTJKSkqwffv2Jn3/0tJSKJVKqNVqKBSKNstNbUsIgaRfryIhKQe7ThdB+9s0Tm+FI6YM6ozJA/3h5eoocUoiIjKW5rx+m9SZnZtptVqsW7cON27cgEqluuM+SUlJiIyMbLBt9OjRmDt3bqPfV6PRQKPR6B+XlpY2ui9Jr1xTiw3HLiIhqeEwzsH3uSNaFYiRIR1hz2GcRER0FyZXdjIyMqBSqVBVVQUXFxds3LgRISEhd9y3sLAQHTt2bLCtY8eOKC0tRWVlJZycbr8Lbnx8POLi4gySndrO+aIyfHM4B9+l/j6Ms52DLZ7o2wnRqkAO4yQioiYzubLTvXt3pKenQ61WY/369YiJicH+/fsbLTzNFRsbi5deekn/uLS0FP7+/m3yval1arU6/Hi6CCt+bjiM8z5PZ8RwGCcREbWQyZUdBwcHBAcHAwD69euHlJQULF68GEuXLr1tX29vbxQVFTXYVlRUBIVCccezOgAgl8shl8vbPji12OUyDb5NqRvGWaD+fRhnZM+OiBnCYZxERNQ6Jld2bqXT6RpcY3MzlUqFrVu3Nti2a9euRq/xIdNRN4yzBAlJ2RzGSUREBmVSZSc2NhaPPPIIOnfujLKyMqxevRr79u3Djh079M/n5+cjISEBADBz5kx8/PHHePXVVzFt2jTs2bMHa9euxZYtW6RcBt1F/TDOhKQcnLr0+8Xhvf3dEDMkAI/28uEwTiIialMmVXaKi4sRHR2NgoICKJVKhIeHY8eOHRg5ciQAoKCgALm5ufr9g4KCsGXLFvz1r3/F4sWL4efnhy+//JL32DFBOVdvYOXhHKw9ymGcRERkXCZ/nx1D4312DKd+GGdCUjb23TKM8+lBAZjYn8M4iYioZSziPjtkvkoqqrHu6EWsTM5BztXfh3EO61Y3jPOh7hzGSURExsOyQ23mZL4a3yTlYPPxfFTV/D6M80+/DeMM4owqIiKSAMsOtUp1rQ7bThYgISkHqTnX9dt7+igQrQrA+N6+aOfAP2ZERCQdvgpRixSo64Zxrjny+zBOOxsZHunlgxgO4yQiIhPCskNNJoRAUtZVfJOUg52Zvw/j7KiQY8rAADw10B9eCg7jJCIi08KyQ/dUrqnFxt+GcZ6/aRjnoKC6YZyjQjmMk4iITBfLDjXql+IyJCTlYMOxfJRragH8PowzanAguntzGCcREZk+lh1qoG4YZzESkrLx868Nh3FGDw7AE/38oOAwTiIiMiMsOwQAuFKuQeKRXKxOzsWlW4ZxRqsCMTSYwziJiMg8sexYsfphnN8kZWNrRiGqtXX3xnF3dsDkAf6YOpjDOImIyPyx7Fihqhotvk+/hITD2TiZ33AYZ7Sqbhinoz2HcRIRkWVg2bEiuVcrsDI5B2uP5qGkom4Yp4N+GGcAwv3cpA1IRERkACw7Fk6nE9h//jK+ScrB3rPF+mGcfu2d8PTgumGc7hzGSUREFoxlx0KpK2qwLjUP3xxuOIzzwW6eiB4cgOE9OIyTiIisA8uOhTl1qW4Y56b034dxujra4U/9/PH04M64z9NF4oRERETGxbJjARobxtnD2xUxQwI5jJOIiKwaXwHNWKG6CquTc7D6SB6ulGsA1A3jHBPmjZghgejPYZxEREQsO+amsWGcXq5yTB3EYZxERES3YtkxExzGSURE1DIsOybul+JyfJOUje9uGcb5eJ9OiFIFoIe3QuKEREREpo1lxwTVD+P85nA2Dv1y0zDODs6IUgXgjxzGSURE1GQsOybkSrkG36bkYdXhnAbDOEf07IhoVQCGdukAG94bh4iIqFlYdiQmhEBaXgm+ScrBlhMFDYZxThrgj6mDOsOvfTuJUxIREZkvlh2JVNVo8f3xS/gmKQcZ+Wr99gh/N0QPDsDYcA7jJCIiagssO0aWd60CKw/n4NtbhnGOC68bxhnh7yZtQCIiIgvDsmMEOp3AgfOXkXDLMM5ObnXDOCcN4DBOIiIiQ2HZMaD6YZwrD+cg+6ZhnA907YAYVSCHcRIRERkBy46BHMu9jilfHOYwTiIiIomx7BhIqK8Czg52CPSQI1oViAl9OIyTiIhICnz1NRC5nS1++Mv98FE6chgnERGRhFh2DMjXzUnqCERERFbPpCZHxsfHY8CAAXB1dYWXlxcmTJiAs2fP3vO4VatWISIiAu3atYOPjw+mTZuGq1ev3vM4IiIisnwmVXb279+PWbNm4fDhw9i1axdqamowatQo3Lhxo9FjDh06hOjoaEyfPh2nTp3CunXrcOTIEcyYMcOIyYmIiMhUmdTbWNu3b2/wePny5fDy8kJqaioefPDBOx6TlJSEwMBAzJ49GwAQFBSE559/Hu+++67B8xIREZHpM6kzO7dSq+vGKLi7uze6j0qlQl5eHrZu3QohBIqKirB+/Xo8+uijd9xfo9GgtLS0wRcRERFZLpMtOzqdDnPnzsXQoUMRFhbW6H5Dhw7FqlWrMGnSJDg4OMDb2xtKpRKffPLJHfePj4+HUqnUf/n7+xtqCURERGQCTLbszJo1CydPnkRiYuJd98vMzMScOXMwf/58pKamYvv27cjOzsbMmTPvuH9sbCzUarX+Ky8vzxDxiYiIyETIhKif1GQ6XnzxRWzevBkHDhxAUFDQXfeNiopCVVUV1q1bp9/2008/4YEHHsClS5fg4+Nz1+NLS0uhVCqhVquhUCjaJD8REREZVnNev03qzI4QAi+++CI2btyIPXv23LPoAEBFRQVsbBouw9bWVv/9iIiIyLqZVNmZNWsWVq5cidWrV8PV1RWFhYUoLCxEZWWlfp/Y2FhER0frH48bNw4bNmzAkiVLkJWVhUOHDmH27NkYOHAgfH19pVgGERERmRCT+uj5kiVLAAAPPfRQg+3Lli3DM888AwAoKChAbm6u/rlnnnkGZWVl+Pjjj/Hyyy/Dzc0NDz/8MD96TkRERABM9JodY+I1O0RERObHbK/ZISIiImprLDtERERk0Uzqmh0p1L+LxzspExERmY/61+2mXI1j9WWnrKwMAHgnZSIiIjNUVlYGpVJ5132s/gJlnU6HS5cuwdXVFTKZzKA/q7S0FP7+/sjLy7OKi6G5XsvG9Vo2a1svYH1rNvf1CiFQVlYGX1/f2+63dyurP7NjY2MDPz8/o/5MhUJhln+wWorrtWxcr2WztvUC1rdmc17vvc7o1OMFykRERGTRWHaIiIjIorHsGJFcLseCBQsgl8uljmIUXK9l43otm7WtF7C+NVvTeq3+AmUiIiKybDyzQ0RERBaNZYeIiIgsGssOERERWTSWHSIiIrJoLDsG9vbbb2PIkCFo164d3NzcmnSMEALz58+Hj48PnJycEBkZifPnzxs2aBu5du0apk6dCoVCATc3N0yfPh3l5eV3Paa8vBwvvvgi/Pz84OTkhJCQEHz22WdGStw6LVkvAJw+fRqPPfYYlEolnJ2dMWDAAOTm5hohceu0dL31Zs6cCZlMhg8//NBwIdtQc9dbU1OD1157Db169YKzszN8fX0RHR2NS5cuGTF1033yyScIDAyEo6MjBg0ahCNHjtx1/3379qFv376Qy+UIDg7G8uXLjRO0jTRnvRs2bMDIkSPh6ekJhUIBlUqFHTt2GDFt6zX391vv0KFDsLOzQ+/evQ0b0JgEGdT8+fPFBx98IF566SWhVCqbdMzChQuFUqkUmzZtEsePHxePPfaYCAoKEpWVlYYN2wbGjBkjIiIixOHDh8XBgwdFcHCweOqpp+56zIwZM0SXLl3E3r17xYULF8TSpUuFra2t2Lx5s5FSt1xL1vvLL78Id3d38corr4hjx46JX375RWzevFkUFRUZKXXLtWS99TZs2CAiIiKEr6+v+M9//mPYoG2kuestKSkRkZGR4ttvvxVnzpwRSUlJYuDAgaJfv35GTN00iYmJwsHBQXz99dfi1KlTYsaMGcLNza3RP4dZWVmiXbt24qWXXhKZmZnio48+Era2tmL79u1GTt4yzV3vnDlzxLvvviuOHDkizp07J2JjY4W9vb04duyYkZO3THPXW+/69evivvvuE6NGjRIRERHGCWsELDtGsmzZsiaVHZ1OJ7y9vcX777+v31ZSUiLkcrlYs2aNARO2XmZmpgAgUlJS9Nu2bdsmZDKZyM/Pb/S40NBQ8eabbzbY1rdvX/H6668bLGtbaOl6J02aJJ5++mljRGxTLV2vEEJcvHhRdOrUSZw8eVIEBASYRdlpzXpvduTIEQFA5OTkGCJmiw0cOFDMmjVL/1ir1QpfX18RHx9/x/1fffVVERoa2mDbpEmTxOjRow2as600d713EhISIuLi4gwRr821dL2TJk0Sb7zxhliwYIFFlR2+jWViLly4gMLCQkRGRuq3KZVKDBo0CElJSRImu7ekpCS4ubmhf//++m2RkZGwsbFBcnJyo8cNGTIE33//PfLz8yGEwN69e3Hu3DmMGjXKGLFbrCXr1el02LJlC7p164bRo0fDy8sLgwYNwqZNm4yUuuVa+vvV6XSIiorCK6+8gtDQUGNEbRMtXe+t1Go1ZDJZk9/GNobq6mqkpqY2+HvGxsYGkZGRjf49k5SU1GB/ABg9erTJ/70EtGy9t9LpdCgrK4O7u7uhYraZlq532bJlyMrKwoIFC4wR06hYdkxMYWEhAKBjx44Ntnfs2FH/nKkqLCyEl5dXg212dnZwd3e/a/aPPvoIISEh8PPzg4ODA8aMGYNPPvkEDz74oKEjt0pL1ltcXIzy8nIsXLgQY8aMwc6dO/H444/jiSeewP79+40Ru8Va+vt99913YWdnh9mzZxs6Yptq6XpvVlVVhddeew1PPfWUSQ1avHLlCrRabbP+niksLLzj/qWlpaisrDRY1rbQkvXe6t///jfKy8sxceJEQ0RsUy1Z7/nz5zFv3jysXLkSdnaWNyOcZacF5s2bB5lMdtevM2fOSB2zzRh6vR999BEOHz6M77//HqmpqVi0aBFmzZqFH3/8sQ1X0XSGXK9OpwMAjB8/Hn/961/Ru3dvzJs3D3/4wx8kuyjbkOtNTU3F4sWLsXz5cshksjZO3jLG+v+3pqYGEydOhBACS5YsaYPkJJXVq1cjLi4Oa9euva0AWwKtVospU6YgLi4O3bp1kzqOQVhefTOCl19+Gc8888xd97nvvvta9L29vb0BAEVFRfDx8dFvLyoqkuzK+Kau19vbG8XFxQ2219bW4tq1a/p13aqyshJ///vfsXHjRowdOxYAEB4ejvT0dPz73/++7bS5MRhyvR06dICdnR1CQkIabO/Zsyd++umnVuVuKUOu9+DBgyguLkbnzp3127RaLV5++WV8+OGHyM7Obm38ZjPkeuvVF52cnBzs2bPHpM7qAHV/Dm1tbVFUVNRge1FRUaNr8/b2vuP+CoUCTk5OBsvaFlqy3nqJiYl47rnnsG7dOkn+PmqJ5q63rKwMR48eRVpaGl588UUAdf8wE0LAzs4OO3fuxMMPP2yU7IbCstMCnp6e8PT0NMj3DgoKgre3N3bv3q0vN6WlpUhOTsYLL7xgkJ95L01dr0qlQklJCVJTU9GvXz8AwJ49e6DT6TBo0KA7HlNTU4OamhrY2DQ8yWhra6s/C2Jshlyvg4MDBgwYgLNnzzbYfu7cOQQEBLQ+fAsYcr1RUVF3vM4jKioKzz77bOvDt4Ah1wv8XnTOnz+PvXv3wsPDo82ytxUHBwf069cPu3fvxoQJEwDUvbjt3r1b/2J3K5VKha1btzbYtmvXLqhUKkPHbbWWrBcA1qxZg2nTpiExMVH/jzFz0Nz1KhQKZGRkNNj26aefYs+ePVi/fj2CgoKMEduwJL5A2uLl5OSItLQ0ERcXJ1xcXERaWppIS0sTZWVl+n26d+8uNmzYoH+8cOFC4ebmJjZv3ixOnDghxo8fb1YfPe/Tp49ITk4WP/30k+jatettH9W9db3Dhg0ToaGhYu/evSIrK0ssW7ZMODo6ik8//dTY8ZutJevdsGGDsLe3F59//rk4f/68/iO8Bw8eNHb8ZmvJem9lLp/GEqL5662urhaPPfaY8PPzE+np6aKgoED/pdFopFhCoxITE4VcLhfLly8XmZmZ4s9//rNwc3MThYWFQggh5s2bJ6KiovT713/0/JVXXhGnT58Wn3zyidl99Lw56121apWws7MTn3zySYPfY0lJiVRLaJbmrvdWlvZpLJYdA4uJiREAbvvau3evfh8AYtmyZfrHOp1O/OMf/xAdO3YUcrlcjBgxQpw9e9b44Vvg6tWr4qmnnhIuLi5CoVCIZ599tkGxE+L29RYUFIhnnnlG+Pr6CkdHR9G9e3exaNEiodPpjJy++VqyXiGE+Oqrr0RwcLBwdHQUERERYtOmTUZM3XItXe/NzKnsNHe9Fy5cuOP/77f+P28qPvroI9G5c2fh4OAgBg4cKA4fPqx/LiYmRgwbNqzB/nv37hW9e/cWDg4O4r777rvr79kUNWe9w4YNu+PvMSYmxvjBW6i5v9+bWVrZkQkhhFFPJREREREZET+NRURERBaNZYeIiIgsGssOERERWTSWHSIiIrJoLDtERERk0Vh2iIiIyKKx7BAREZFFY9khIiIii8ayQ0RERBaNZYeIiIgsGssOEZmVhx56CHPnzm3WMVevXoWXlxeys7MbbP/b3/6GCRMmYPLkyVi0aFHbhSQik8KyQ0QW7+2338b48eMRGBjYYHt6ejp69+6NN954A2+//TbUarU0AYnIoFh2iMiiVVRU4KuvvsL06dNve+748ePo3bs3wsLC0KVLF6xcuVKChERkaCw7RGS2HnroIcyePRuvvvoq3N3d4e3tjX/+858N9tm6dSvkcjkGDx7cYPvFixdx5coVREREAADGjRuHxMREY0UnIiNi2SEis7ZixQo4OzsjOTkZ7733Ht58803s2rVL//zBgwfRr1+/245LT0+HUqlEUFAQAGDgwIE4cuQINBqN0bITkXGw7BCRWQsPD8eCBQvQtWtXREdHo3///ti9e7f++ZycHPj6+t52XHp6uv6sDgD4+vqiuroahYWFRslNRMbDskNEZi08PLzBYx8fHxQXF+sfV1ZWwtHR8bbjbi07Tk5OAOqu8SEiy8KyQ0Rmzd7evsFjmUwGnU6nf9yhQwdcv379tuPqP4lV79q1awAAT09PwwQlIsmw7BCRRevTpw8yMzMbbCsrK0NWVlaDsnPy5En4+fmhQ4cORk5IRIbGskNEFm306NE4depUg7M7x48fh62tLUJDQ/XbDh48iFGjRkkRkYgMjGWHiCxar1690LdvX6xdu1a/LT09HT169IBcLgcAVFVVYdOmTZgxY4ZUMYnIgGRCCCF1CCIiQ9qyZQteeeUVnDx5EjY2t/8bb8mSJdi4cSN27twpQToiMjQ7qQMQERna2LFjcf78eeTn58Pf3/+25+3t7fHRRx9JkIyIjIFndoiIiMii8ZodIiIismgsO0RERGTRWHaIiIjIorHsEBERkUVj2SEiIiKLxrJDREREFo1lh4iIiCwayw4RERFZNJYdIiIismj/D2GQq17noidYAAAAAElFTkSuQmCC"
          },
          "metadata": {}
        }
      ]
    },
    {
      "cell_type": "code",
      "source": "m, b =np.polyfit (X,Y,1)\n",
      "metadata": {
        "trusted": true
      },
      "execution_count": 49,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "",
      "metadata": {},
      "execution_count": null,
      "outputs": []
    }
  ]
}