diff --git a/doc/pub/week40/ipynb/week40.ipynb b/doc/pub/week40/ipynb/week40.ipynb index 0da89492e..c457584b6 100644 --- a/doc/pub/week40/ipynb/week40.ipynb +++ b/doc/pub/week40/ipynb/week40.ipynb @@ -611,7 +611,15 @@ "execution_count": 5, "id": "17635e19", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "gamma_j after 500 epochs: 9.97108e-05\n" + ] + } + ], "source": [ "import numpy as np \n", "\n", @@ -654,7 +662,34 @@ "execution_count": 6, "id": "a21b3bf0", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Own inversion\n", + "[[4.00444539]\n", + " [3.05946757]]\n", + "Eigenvalues of Hessian Matrix:[0.33390381 4.08547851]\n", + "theta from own gd\n", + "[[4.00444539]\n", + " [3.05946757]]\n", + "theta from own sdg\n", + "[[3.93407215]\n", + " [3.12257566]]\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAHFCAYAAADxOP3DAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAABV40lEQVR4nO3deXwTZeIG8Ce9uVpooaUXtCBSoNy3UAHlkKMBKx6gCOzqT+USwVVYXAVd7MqKoqKyeIBXwYOCKSqK0gIKuJxynxYohVLK0dJCr/T9/ZFNIDRpjiaZmczz/Xzywc5MJu8kHfP0PTVCCAEiIiIilfKRugBEREREUmIYIiIiIlVjGCIiIiJVYxgiIiIiVWMYIiIiIlVjGCIiIiJVYxgiIiIiVWMYIiIiIlVjGCIiIiJVYxgikqHly5dDo9GYHn5+foiMjMRDDz2EY8eOSVauuXPnQqPRSPb6SnDy5EloNBq8/vrrUheFiOzkJ3UBiMi6ZcuWISEhAaWlpfjtt98wf/58ZGZm4vDhw2jUqJHUxSMi8goMQ0QylpiYiG7dugEA+vfvD71ej5deeglr1qzBxIkTJS4dSamiosJUa0hEtcNmMiIFMQaj8+fPm7aVlpZi5syZ6NSpE0JCQhAaGorevXvj22+/rfZ8jUaDKVOm4LPPPkObNm1Qt25ddOzYEWvXrq127HfffYdOnTohMDAQ8fHxVpt9SktLMXv2bMTHxyMgIADR0dGYPHkyrly5YnZcXFwcRowYgbVr16Jz586oU6cO2rRpY3rt5cuXo02bNqhXrx569OiBHTt22Hw/jM2JmZmZeOqpp9C4cWOEhYUhJSUFZ8+erXbtc+fOrXaOuLg4TJgwodo5N2zYgMcffxxhYWEIDg7Go48+ipKSEuTl5eGBBx5Aw4YNERkZiWeffRYVFRXVzltVVYX58+ejWbNmCAoKQrdu3fDLL79UO+7YsWMYO3YswsPDERgYiDZt2uDdd981OyYrKwsajQafffYZZs6ciejoaAQGBuL48eO4du0ann32WcTHxyMoKAihoaHo1q0bVqxYYfP9IyID/klBpCDZ2dkAgNtvv920raysDJcuXcKzzz6L6OholJeX4+eff0ZKSgqWLVuGRx991Owc3333HbZv346XX34Z9evXx4IFC3DvvffiyJEjaNGiBQDgl19+wciRI9G7d2+sXLkSer0eCxYsMAthACCEwKhRo/DLL79g9uzZSEpKwt69e/HSSy9h69at2Lp1KwIDA03H//HHH5g9ezbmzJmDkJAQzJs3DykpKZg9ezZ++eUXvPrqq9BoNHj++ecxYsQIZGdno06dOjbfl8ceewzDhw9HWloacnJy8Le//Q2PPPIINmzY4PR7/dhjjyElJQUrV67E7t278fe//x2VlZU4cuQIUlJS8H//93/4+eef8dprryEqKgozZswwe/7ixYvRvHlzLFq0CFVVVViwYAGGDh2KjRs3onfv3gCAgwcP4o477kCzZs2wcOFCNG3aFD/++COmTZuGgoICvPTSS2bnnD17Nnr37o0lS5bAx8cH4eHhmDFjBj777DP885//ROfOnVFSUoL9+/fj4sWLTl87keoIIpKdZcuWCQBi27ZtoqKiQly9elWsW7dONG3aVNx5552ioqLC6nMrKytFRUWF+Otf/yo6d+5stg+AiIiIEEVFRaZteXl5wsfHR6Smppq29ezZU0RFRYnr16+bthUVFYnQ0FBx8/821q1bJwCIBQsWmL3Ol19+KQCIpUuXmrY1b95c1KlTR5w5c8a0bc+ePQKAiIyMFCUlJabta9asEQCETqez632aNGmS2fYFCxYIAOLcuXNm1/7SSy9VO0fz5s3F+PHjq51z6tSpZseNGjVKABBvvPGG2fZOnTqJLl26mH7Ozs4WAKy+fwMHDjRtGzJkiIiJiRGFhYVm55wyZYoICgoSly5dEkIIkZmZKQCIO++8s1r5ExMTxahRo6ptJyL7sZmMSMZ69eoFf39/NGjQAPfccw8aNWqEb7/9tlo/ka+//hp9+vRB/fr14efnB39/f3z00Uc4dOhQtXMOGDAADRo0MP0cERGB8PBwnDp1CgBQUlKC7du3IyUlBUFBQabjGjRogOTkZLNzGWtebm5mAoD7778f9erVq9Ys1KlTJ0RHR5t+btOmDQBDf6i6detW224sky1ardbs5w4dOjj0fEtGjBhh9rOxTMOHD6+23dLrWHv/Nm3aBL1ej9LSUvzyyy+49957UbduXVRWVpoew4YNQ2lpKbZt22Z2zvvuu6/a6/To0QM//PADZs2ahaysLFy/ft3payZSK4YhIhn79NNPsX37dmzYsAFPPPEEDh06hDFjxpgdk56ejgceeADR0dH4/PPPsXXrVmzfvh1/+ctfUFpaWu2cYWFh1bYFBgaavkQvX76MqqoqNG3atNpxt267ePEi/Pz80KRJE7PtGo0GTZs2rdZUExoaavZzQEBAjdstld+SW6/J2DRXm2DgSFktldPa+1deXo7i4mJcvHgRlZWVeOedd+Dv72/2GDZsGACgoKDA7PmRkZHVzvn222/j+eefx5o1azBgwACEhoZi1KhRkk7BQKQ07DNEJGNt2rQxdZoeMGAA9Ho9PvzwQ3zzzTcYPXo0AODzzz9HfHw8vvzyS7M5gMrKypx6zUaNGkGj0SAvL6/avlu3hYWFobKyEhcuXDALREII5OXloXv37k6VwR0CAwMtvifu6ltj7f0LCAhA/fr14e/vD19fX4wbNw6TJ0+2eI74+Hizny3N8VSvXj3MmzcP8+bNw/nz5021RMnJyTh8+LBrLobIy7FmiEhBFixYgEaNGuHFF19EVVUVAMMXZEBAgNkXZV5ensXRZPYwjuZKT083q/G4evUqMjIyzI69++67ARgC2c1WrVqFkpIS0345iIuLw969e822bdiwAcXFxW55PWvvX1JSEnx9fVG3bl0MGDAAu3fvRocOHdCtW7dqD0u1eDWJiIjAhAkTMGbMGBw5cgTXrl1z9WUReSXWDBEpSKNGjTB79mw899xzSEtLwyOPPIIRI0YgPT0dkyZNwujRo5GTk4NXXnkFkZGRTjeVvPLKK7jnnnswaNAgzJw5E3q9Hq+99hrq1auHS5cumY4bNGgQhgwZgueffx5FRUXo06ePaTRZ586dMW7cOFddeq2NGzcO//jHP/Diiy+iX79+OHjwIBYvXoyQkBC3vJ6vry8GDRqEGTNmoKqqCq+99hqKioowb9480zFvvfUW+vbti6SkJDz11FOIi4vD1atXcfz4cWRkZNg1Gq5nz54YMWIEOnTogEaNGuHQoUP47LPP0Lt3b7N+WERkHcMQkcJMnToVixcvxssvv4wxY8Zg4sSJyM/Px5IlS/Dxxx+jRYsWmDVrFs6cOWP2xeuIQYMGYc2aNXjhhRfw4IMPomnTppg0aRKuX79udk6NRoM1a9Zg7ty5WLZsGebPn4/GjRtj3LhxePXVV82G1Uvtb3/7G4qKirB8+XK8/vrr6NGjB7766iuMHDnSLa83ZcoUlJaWYtq0acjPz0e7du3w3XffoU+fPqZj2rZti127duGVV17BCy+8gPz8fDRs2BCtWrUy9Ruy5a677oJOp8Obb76Ja9euITo6Go8++ijmzJnjlusi8kYaIYSQuhBEREREUmGfISIiIlI1hiEiIiJSNYYhIiIiUjWGISIiIlI1hiEiIiJSNYYhIiIiUjVVzDNUVVWFs2fPokGDBhansyciIiL5EULg6tWriIqKgo+P++pvVBGGzp49i9jYWKmLQURERE7IyclBTEyM286vijDUoEEDAIY3Mzg4WOLSEBERkT2KiooQGxtr+h53F1WEIWPTWHBwMMMQERGRwri7iws7UBMREZGqMQwRERGRqjEMERERkaoxDBEREZGqMQwRERGRqjEMERERkaoxDBEREZGqMQwRERGRqjEMERERkaoxDBEREZGqMQwRERGRqjEMERERkaoxDBEREZGqMQwRERGRqjEMERERkaoxDBEREZGqMQwRERGRqkkehjZt2oTk5GRERUVBo9FgzZo1Vo994oknoNFosGjRIo+Vj4iIiLyb5GGopKQEHTt2xOLFi2s8bs2aNfj9998RFRXloZIRERGRGvhJXYChQ4di6NChNR6Tm5uLKVOm4Mcff8Tw4cM9VDIiIiJSA8lrhmypqqrCuHHj8Le//Q3t2rWTujhERETkZSSvGbLltddeg5+fH6ZNm2b3c8rKylBWVmb6uaioyB1FIyIiIi8g65qhnTt34q233sLy5cuh0Wjsfl5qaipCQkJMj9jYWDeWkoiIiJRM1mFo8+bNyM/PR7NmzeDn5wc/Pz+cOnUKM2fORFxcnNXnzZ49G4WFhaZHTk6O5wpNREREiiLrZrJx48Zh4MCBZtuGDBmCcePGYeLEiVafFxgYiMDAQHcXj4iIiLyA5GGouLgYx48fN/2cnZ2NPXv2IDQ0FM2aNUNYWJjZ8f7+/mjatClat27t6aISERGRF5I8DO3YsQMDBgww/TxjxgwAwPjx47F8+XKJSkVERERqIXkY6t+/P4QQdh9/8uRJ9xWGiIiIVEfWHaiJiIiI3I1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFTNT+oCEBERkYP0emDzZuDcOSAyEkhKAnx9pS6VYkleM7Rp0yYkJycjKioKGo0Ga9asMe2rqKjA888/j/bt26NevXqIiorCo48+irNnz0pXYCIiIimlpwNxccCAAcDYsYZ/4+IM28kpkoehkpISdOzYEYsXL66279q1a9i1axf+8Y9/YNeuXUhPT8fRo0eh1WolKCkREZHE0tOB0aOBM2fMt+fmGrYzEDlFI4QQUhfCSKPRYPXq1Rg1apTVY7Zv344ePXrg1KlTaNasmV3nLSoqQkhICAoLCxEcHOyi0hIREXmQXm+oAbo1CBlpNEBMDJCd7TVNZp76/pa8ZshRhYWF0Gg0aNiwodVjysrKUFRUZPYgIiJStM2brQchABACyMkxHEcOUVQYKi0txaxZszB27NgaE2JqaipCQkJMj9jYWA+WkoiIyA3OnXPtcWSimDBUUVGBhx56CFVVVXjvvfdqPHb27NkoLCw0PXJycjxUSiIiIjeJjHTtcWSiiKH1FRUVeOCBB5CdnY0NGzbYbDcMDAxEYGCgh0pHRETkAUlJhj5BubmGJrFbGfsMJSV5vmwKJ/uaIWMQOnbsGH7++WeEhYVJXSQiIiLP8/UF3nrL8N8ajfk+48+LFnlN52lPkjwMFRcXY8+ePdizZw8AIDs7G3v27MHp06dRWVmJ0aNHY8eOHfjiiy+g1+uRl5eHvLw8lJeXS1twIiIiT0tJAb75BoiONt8eE2PYnpIiTbkUTvKh9VlZWRgwYEC17ePHj8fcuXMRHx9v8XmZmZno37+/Xa/BofVERORVVDIDtae+vyXvM9S/f3/UlMdkNA0SERGRPPj6AnZWCJBtkjeTEREREUmJYYiIiIhUjWGIiIiIVE3yPkNERESkchJ3CGcYIiIiIumkpwNPP22+7lpMjGFOpYEDPVIEhiEiIiKSRno6MHp09Rm1c3MN2z/91CPFYJ8hIiIi8jy93lAjZGkKHeO2WbM8UhSGISIiIvK8zZvNm8ZuJYShhsgDGIaIiIjI886dk7oEJgxDRERE5HmRkVKXwIRhiIiIiDwvKckwakyjsbxfo6m+IK2bMAwRERGR5/n6GobPA9UDkfHnf/3LI0VhGCIiIiJppKQA33xTvQYoJsawXav1SDE4zxARERFJJyUFGDnS8gzURUUeKQLDEBEREUnL1xfo31+yl2cYIiIiohskXidMCgxDREREZFDTOmEpKdKVy83YgZqIiIhurBN266zQxnXC0tOlKZcHMAwRERGpnT3rhE2fbjjOCzEMERERqZ0964Tl5BiO80IMQ0RERGpn7zphMlpPzJUYhoiIiNTO3nXCZLSemCsxDBEREamdPeuExcYajvNCDENERERqZ886YYsWee18QwxDREREZHudMC+eZ4iTLhIREZFBTeuEeTGGISIi8l4qXFqi1iReJ0wKDENEROSdVLq0BDmOfYaIiMj7qHhpCXIcwxAREXkXlS8tQY5jGCIiIu+i8qUlyHHsM0RERN5F5UtLyILCOq4zDBERkXdR+dISklNgx3U2kxERkXdR+dISklJox3WGISIi8i4qX1pCMgruuM4wRERE3se4tERUlPn26GivX1pCMvPnK7bjOvsMERGRsjjSOddaU5mSybFzcno68NJL9h0rw47rrBkiIiLlSE8H4uKAAQOAsWMN/8bFVe+LotC+KzbZe/2uptcDWVnAihWGf29u6jI2j9lLhh3XJQ9DmzZtQnJyMqKioqDRaLBmzRqz/UIIzJ07F1FRUahTpw769++PAwcOSFNYIiKSjr0BR8F9V2okVcCzFcBszet0Mzs7rl85VYiV07bgr4nbnC62IyQPQyUlJejYsSMWL15scf+CBQvwxhtvYPHixdi+fTuaNm2KQYMG4erVqx4uKRERScaRgOONky5KFfDsCWCONHs99JDVJr3sTTl4K2UjBobuQpO4uhjzzh34JqdXLQpvP8n7DA0dOhRDhw61uE8IgUWLFmHOnDlI+V9nt08++QQRERFIS0vDE0884cmiEhGRVBwJON446aIj1++qFedtBTCNxhDAli2z/5yvvw706gWkpKCqsgrbPzkI3ccF0O2Mxv6yVgBiTYe2DTyOwYlHsGhnra/EJsnDUE2ys7ORl5eHwYMHm7YFBgaiX79+2LJli9UwVFZWhrKyMtPPRUVFbi8rERG5kSMBxxsnXczNte84VwY8ewMYYJjXKTfXcnC6+SkArk+chKdnhSLjRFucr0o07fNFJZIa7oO2XyGSp8bjtrtvQ1FROBaFuOBabJC8mawmeXl5AICIiAiz7REREaZ9lqSmpiIkJMT0iI2NtXosEREpgCMBx5smXdTrgZdfBp56yr7jXRnw7A1W+fk35nWyQSME6hadx7FjwPmqcASjEA/GbsEXk37DhT+LkXm5M55Z0x+33d28FgV3nKzDkJHmll9oIUS1bTebPXs2CgsLTY8cY3IlIiJlciTgeMuki+npQESEYci6rX6y7gh4jgRQ47xOoaF2PWVC7Aasf20XLlytg5Wn78DYd/ugUXxD58taS7IOQ02bNgWAarVA+fn51WqLbhYYGIjg4GCzBxERyVxNw7cdDTjGL+foaPNjY2KUMemisePyxYv2P8fVAc+BAFpeXI6fj8fhnZA5dp16wqd3YeBzXRBQP8B15a0FWYeh+Ph4NG3aFOvXrzdtKy8vx8aNG3HHHXdIWDIiInIpe+bPcTTgpKQAJ08CmZlAWprh3+xs+QehmjouW9KkiXsCno0AKgBsavskHor/HU0aXMeg57tgevbTyEEMqqCsJkrJO1AXFxfj+PHjpp+zs7OxZ88ehIaGolmzZpg+fTpeffVVtGrVCq1atcKrr76KunXrYuzYsRKWmoiIXMZYC3Lrl79x+PbNX/QpKcDIkfbPwOzr69zoKilneXZk3h4AePNN9wU8YwC9ZRX68wjHZPEOVv14v2lbhE8+Rtx2BHmJkxGz+u8ANOafqYybKCUPQzt27MCAAQNMP8+YMQMAMH78eCxfvhzPPfccrl+/jkmTJuHy5cvo2bMnfvrpJzRo0ECqIhMRkavYO3x75MgbX6DOBhx7padX+/JHTIyhlsQTtUqOjgi7tabMhfTlevy3oBUyopfj7Pk8lFX44BwisVkkoQq+SAw8Bm3XXGj/0hjdx7eFj184gCQg/XbL7+GiRbKsmdMIYW89nHIVFRUhJCQEhYWF7D9ERCQnWVmGJjFbMjPdG4CMrNVSGWs1PNHfyN73BDA0OWVnu7SmpSS/BOvf3A/dqnKsPZ6AC6KJaZ8fKtCv0T4k97+K5GnxaNG/mfUTuaB2zVPf35LXDBERkYrJaYJEZ2qp3MHYcdmOeXtc1eR0dlceMhYeRcbPQfg5vwPK0NO0LwSFGNZ8P7RaDe6Z2Q4Nm3ex76TursFzIYYhIiKSjpwmSJRilmdLjB2XR482BDBLgSgsDFi61OlaKlElsPebo9D95xx028Kx41pbAE1N++P9TmNk+2wkPxyMpKcS4V+3j5MXowwMQ0REJB1btSAajWG/J0YfyamWykrHZYSFAdOmAXPmOFwjVFZUho3v7oduRTF0B29Djr41gNYAAA2q0LP+AWjvuAjtpBi0TW4JjU8NTWBehmGIiIikU1MtiKdHH8mplgpwfOScBRePXcIPbxyCbq0P1p1ph6voatpXB9cwuOk+aO8px/CZCYhIbO+Oq1AEdqAmIiLpWRrBFRvr2dFHer1hbiNrTWXGWioXd1h2tWPrT0L39knoNjfEr4XtUYUbZY30yUPy7UeRfH8Q7p7eHnVC60hYUtvYgZqIiNTDBbUgtebrC4wZA/z739aPkeEcOfpyPbZ+eAC65ZeQ8UcsDpe3BBBn2t8h6Ai03c5B+1g4uj6cAB+/plbPpVYMQ0REJA9Sjz5KTwdef936/meflc0cOcV5xfjpjf3QpVfiuz8TUCA6mPb5oxz9Q/dBe1cxkqe3RPM+N/oGkWVsJiMiIlJAE9mZ7eeQ8cYx6H6ugw0FHVCOQNO+RprLGNb8ILSjNBjyTDuENAuRpIyuxmYyIiIiT5HLsPqbX7JKYPeKw8j48Dx0v0dg1/U2AG503m7pdwojO2ZDO64h+jyRCL8g7x7+7k4MQ0RERDIZVl9WVIbMt/dBt/IaMg7fhjP6NgDaADAMf7+jwX5o+15C8qRYJAxrAY1Pc7eWxyYp13BzIYYhIiIiCYfVFxy5iO9ePwTd93746Ww7FKObaV9dlGBI5D5oh1Vi2DOtEd6uQw1n8jCp13BzIfYZIiIiMvYZsjX5o4v6DB354U/oFp+G7tdG2FKUaDb8PcrnHLQJx6B9sA4GTGuPoIZBtX49l/PQGm6e+v5mGCIiIgJufMEDlid/rMUXfGVpJbZ8cAC6Ty5Dt7c5jlXEm+3vVOcwtN3zoH08Al3GJkDjo3HqdTzCg53N2YGaiIjIk6wtgRET49Tkj0VnivDTmwegW6PHd9ltcUl0NO3zRznuCtsL7cASjHj6NjTrnQAgwTXX4W4y7GxeWwxDRERERrWc/PH01lxkvHkcug31kHmxAyrQ27QvVHMJI+IPQpvii8HPJKJBVLcaziRjMuls7koMQ0RERDdzYPLHqsoq7Eo7jIyP8qHb3hR7ricAiDbtv90/G9qOp6Ad3wi9H2sHv6C+7imzJ8ltDTcXYBgiIiJyQOmVUmx4ax90X15HxpFWOFvVFkBbAIAP9OgTvB/apMtIntwMrYe2ABBf4/kUJynJ0HRoq7N5UpLny+YkhiEiIiIb8g9cwHcLD0O3zh8/nWuPa+hu2lcfVzEkej+0w/QYNrMNGrfuaPkkXjInD3x9DcPnR482BB9Lnc1luIZbTRiGiIiIbiGqBA6tPYGM989A91sotl5NhMCNmo4Y37PQtjkG7UP10H9qewQG967hbPCqOXkAuLyzudQ4tJ6ISEnsrV3wlloID6osrcSvS/ZD99kV6PbG40Sl+ezOXesehLZnPpIfb4pOD7a2f/i7h+bkkYSbf884z5ALMQwRkVewt3bB22oh3KjwdCHWLTyADJ3A96fa4rJoZNoXgDLc3WQvtAOvY8QzrRDT3YkOwQpYAFbOGIZciGGIiBTP3toFb66FsMbB2omTv55BxqIT0GXWR9alDqiEv2lfY00BRrQ8BG2KPwY9k4j6TevXrmxZWcCAAbaPy8xUzJw8nsRJF4mIyECvN9T0WPrbVQhD0Jk+HRgxwr7jRo6sXS2EnJrg7KgFq6qswo7PDkH30QXodkZhX+ntAGJMhycEnIC2Uw60j4ag1+2X4FuQb7iuJnVqXz4vnJPHGzEMERHJnb0z/r73nvtnBpZTE5y1WrDcXIjRo7HjvlR8sL83Mo62Rl5VO9NuX1Sib8h+aPtdQfKUOLQa1BJI/wN4epzrr8sL5+TxRgxDRERyZ2+twYkTrj3frWoIHxg92rNNcDZqywQ0iPjmXXyEZ1EFXzRAEYbG7od2hMDQmW0R2rLTjePdeV1eOCePN/KRugBERGSDvbUGLVu69nw3s9VUBxia4PR6x8/tDBu1ZT4QaIYcvBH/Dn5K3YmCq0H48vQdePi9PghteaOTtNuvyzgnD3Cj35aRQufk8UYMQ0REcmesXbj1y9RIowFiY4FJk+w7zplaCEcW53SjimsV+OXfu/DphA12Hf/0/AgMmtUVAfUDLB/giesyzskTHW2+PSbGOzu0KxDDEBGRXOn1htFIX30FPP74jU7QN7u5diEgwH21EBJ2BL5yqhArpm7BmOZb0KTeNQx8rgs+PnWXfU+2VQvmqetKSQFOnjSMGktLM/ybnc0gJBPsM0REJEeWOiqHhRn+vXjxxrZbZ/x118zAHu4I/GfWaejeykbGxgbYdLk9KnGHaV8TzQXc1hK4fj4cQcUXoKlNXxxPXpcDC8CSZ3GeISIiubE1V9DcuUCrVp6dgdo4eaCtjsBOTh5YVVmF/y4/CN3HBdDtisaBslZm+9sFHoO2Sy6SJ4Shx4S28A3wvfE+AZbXx7KnCcrN10W1w0kXXYhhiIgk4UwgkfOMxa4IHzcpyS/Bz4v2Q/dNOdYeT0C+aGLa54tK3NlwH7T9i5A8NQ4t72pu+SSWatBiYx2rBXPxdZHrMAy5EMMQEXmcs/PxyH3G4lqGj3N7zmPtG0eg+zEIP+e3RyluTGzYEJfwTHgaBna6iPYTuqLBA0PtC3yuqAVzRagil2MYciGGISLyqNosibFiBTB2rO3XSEsDxoypXTmd5UD4EFUC+9KPQbfkLHTbmmB7STuz/XF+OdC2+xN/7bAd7TcsgiY398ZOT0/mKKeZtQkAw5BLMQwRkcfUtplL7jVDdigvLsfGxfugSytGxsGWOKWPMdvfs95+aO8ogPapaLQbeRs0a1arbz01sgvDkAsxDBGRx9Q2zCi0Q++lE5fxw8KD0K3VYF1OOxQhxLSvDq5hUNN9SB5cjhEzW6Nph/AbT5RzHymSHBdqJSJyJU81gdR23hrjjMWjRxuCgKUOvTKZsfj4L6egezsbuo0N8WthIvToY9rX1Oc8RrQ6Au39Qbj76UTUbdzT8kkcmfRQpjVhpHwMQ0Tk/Ty5uKij89ZYCmnumiuolvTlevy+7CB0yy5CtzsWh8pbArgxyqt90FFou56F9q9N0G1cG/j4Rdg+KVd1JxlgGCIi7+bpxUUdWZjTVkgbOVLyDr3FecVY/+Z+6NIrsPZEGxSI9qZ9fqhA/9C9SO5fjOSnWyD+ztsB3O7YC3BVd5IB2fcZqqysxNy5c/HFF18gLy8PkZGRmDBhAl544QX4+Ni3mgj7DBGplFT9UeyZtwaQbafh3B3nkPHGMWT8XAe/XGiPMgSZ9jXUXMGwZgegHanBPTPbIaRZSA1nsoNC+0iRZ7DP0P+89tprWLJkCT755BO0a9cOO3bswMSJExESEoKnn35a6uIRkZxJ1R/FVjPXyJGGAGBtpXSNxrBS+siRHgkAokrgj6+PQvefc9D9Ho6d19oCuFET08LvFEZ2yEbywyHo+2Qi/Ov2sX4yRymhjxSH3Hs92YehrVu3YuTIkRg+fDgAIC4uDitWrMCOHTskLhkRyZ6U/VFqaubKypK803BZURmy3tkH3coSZBy6DTn61gBaAwA0qEKv+geg7XMR2kkxaDOiJTQ+VmaAdgWZ9pEC4Nn+ZiQZ2Yehvn37YsmSJTh69Chuv/12/PHHH/j111+xaNEiqYtGRHIndX8UawtzShTSLh67hO8XHoRurS/W5SaiGN1M++qiBIMj90F7TwWGz0xAeLv2NZzJDWTSR8qMp/ubkWRkH4aef/55FBYWIiEhAb6+vtDr9Zg/fz7G1DDzallZGcrKykw/FxUVeaKoRCQ3jnRm9iQPhrSjP2ZD984p6DY3wm9FiahC3xun98lD8u1HoX2wDu6alog6ob2qn8CTTURyWtVdrzfUCMmkKZPcS/Zh6Msvv8Tnn3+OtLQ0tGvXDnv27MH06dMRFRWF8ePHW3xOamoq5s2b5+GSErkY+ynUnlz7o7gxpFWWVmLrhweQ8ell6P5ohiPlLQDEm/Z3DDoCbfdz0D4Wji5jE+Dj19T6ydTcRMT5j9RFyFxMTIxYvHix2bZXXnlFtG7d2upzSktLRWFhoemRk5MjAIjCwkJ3F5fINVatEiImRgjD/3INj5gYw3ZynKX3MzZW2vdz1SohNBrD4+ZyGbc5ULai3CLxzbNbxKMtNoswTYHZ6fxRJgaHbRfvjM4SJ3/Ncbx8N5/MyfIpUlpa9Wu39EhLk7qkXq2wsNAj39+yrxm6du1atSH0vr6+qKqqsvqcwMBABAYGurtoRO7BfgquJ8f+KLXsNJzz+1lkvHkMup/rIfNie5Sjt2lfqOYShsUdhHaUL4bMaIfgmG41nMkCe5qInnwSuH4diI6W/r10hD01rno9cP68feerqSmTtbvK4dao5QLjx48X0dHRYu3atSI7O1ukp6eLxo0bi+eee87uc3gqWRLVWmVl9RqMW/8qj401HEfeobJSiMxMQw1DZqbVz7ZKXyV2fn5QvNQvU3Suc7Dar8Zt/tliZtdMsfHtPaLiekXtypSZaV+tiNJqLe2pcbV0jDP3Imt3XcJT39+yn3Tx6tWr+Mc//oHVq1cjPz8fUVFRGDNmDF588UUEBATYdQ5OukiK4ckVy/lXq+yVXilF5tv7oPvyOjIOt0Ju1Y1aCB/ocUfwfiT3uQztlGZofU88ND4a17zwihXA2LH2Hy+DiSJtslbjas9EmLeydb32vJZc3yeZ4ar1LsQwRIph75dQWhpQw4hKm9TcMVbmLhwqwHcLDyPjBz/8eDYRJahv2lcPxRgStR/aYZUYNiMBTdo0dk8h7A3lN5PzTNH2zEQeHW0IL7m5ts8XG2u9KVOqWc+9FGegJlIjTwy5lkufJCXWTLmhzKJK4Mi6bOgWn4bu11BsuZoIcdPw92ifc9C2OYbkB+pgwLT2CGpoYfi7q9ka7WaJnEdX2TMyrKb9N3vzTWDqVOufO0ehOc/S/eUhDENEcuLueXHkMneKEmumXFjmytJK/Paf/dB9dgW6vXE4XtECQAvT/s51DkHb4zy0j0eg85gEaHw8vEhpTVMS2CLH1eVdWaaIiJrvDSlnPVcya/dXaqpnXt+RDkanT592R78lt2MHajezswMo2cmFQ66rsbdjbGamq66mOiUO2XZBmQtzCsVXz2wRj8RvFo00l8xOE4BScU/j/4r3HtooTm/L9cAF2cnezsSe+t1xlqMdwmtzfXK4x5SmhvurEPDI97dDYahu3brihRdeEMXFxe4qj1swDLkRR0y4h7vmxZF67hQljparRZlP/poj3hmdJQaF7hD+KDN7WpimQIxvuVl88+wWUZRbJMGF2cn4x87nnwvRuLGyPjsj42do6QvXWPaYGCGio2s+xp7rs+e15Po+ScHG/SXLMPTbb7+JHj16iMjISPHxxx+7q0wuxzDkJkr8C19J3FHjJvVfrVK/vjMcKLO+Qi/+u/yAeKFvpugYdLjaIa0DToi/dc8Um9/9Q1SWKfDL0J21lu5mT9lddX1Kfp88zcb9JcswZPTJJ5+ImJgY0alTJ5Epp/9pWcEw5AZK/AufpP+rVeqaKWfYWeYPov4hIn3OmW32QaW4M2S3eH1Epjiy7k+pr8Q15Dibt73sKburrk/J75Mn2bi/PBWGnB5af/36daSmpmLhwoUYPHgw/v3vf+O2225zVVcml+LQejfw5Hw45FrG0WSA4X83Rp6YA0WJvzd2lrk/MrER/VEfV3FP9H5oR+gxbGZbhLUKdX8ZPU2JIwGN7J2B2hXXp+T3yVNs3F9FAEIA+c4zdO3aNezatQurVq3C22+/DX9/f0yePBlz585FgwYNXF3OWmEYcgNPzYdD7mFp5EZNc6e4inEOFluj5eQ0B4teDxEXB5zJhQbVyywAXEYjvJT4DZLHhqDf5EQEBnM5ICK72Ph/gqfCkEND65csWYLt27dj+/btOHToEHx9fdGhQwdMnjwZnTp1whdffIG2bdti9erV6NbNwbVwSFk8MR8OuY9Ua3XJdRV5CyquVeDXJfuh+7wIpWdfwHt4EgKApTmeG+Ey3pl3BUi5y8OlJFI4W/9PcK6+xmEO1QzFxsaiV69epke3bt2qLYj66quvIi0tDfv373d5YZ3FmiE3kPIvfFY9K5+lmqkmTYCHHzaENIk+0yunCrFu4QFk6AS+P90OV0RDAIalL84jHGG4ZDEMybJGi0hJrNRWF736KkLGjZNvM5k158+fR1RUFPR6vStPWysMQ24iRd8TJU7WR5YZQ+233wKffw4UFNzY58HPNHtTDjLe+hO6zAbYeLk9KuFv2tdEcwEjbjuMv3Tbh74rJts+mZz6OhEpjYU/dItKSpS5HEd4eDg2bNjg6tOSHKWkGAKPpXDijr4ncllGglzD1xe4dMkQejz4mVZVVmH7Jweh+7gAup3R2F/WCkCsaX+bgBPQds6BdmIYek5sC9+AJGDFGWCFHSfnrMJEzvP1leyPCS7USrXniWYrLn7ofTz4mV4ruIaf39yHjFVlyDiWgPNV4aZ9vqhEUsN90PYrRPLUeNx2d/PqJ1DiKDgiL8CFWkk5PJHm5bL4IfsruY6bP9O8vflYu/AIdD8GYP35DihFT9O+YBRiaOwBaJMF7pnRFqEtO9d8MnevGXcr/p4ReRTDECmDHBY/ZH8l13LxZyqqBPavPgbdkrPI2NoYv5ckArhRA9Tc9wy07U5A+3AD3DkpEQH177C/rJ4cBcffMyKPYxgiZZB6KD/7K7meCz7T8uJybH5/P3RfXIXuQAucrLwdwO2m/T3q7Ye2dwG0T0Yh8d5W0PjEOF9eT/SRc+b3jLVIRLXGPkOkDFIP5Wd/Jddz8jO9nH0FP7x+ALoMDX7IaYcihJj2BeE6Bobvg3ZIKUbMaI3IThHuKbc7woczv2esRSIvxz5DRDeTcrI+ufRXkhNXBAIHPtMTG04h452T0GUFY9OV9tCjj+nQcM0FJLc6DO3oAAx8pj3qNu5R++uzVW53fM6O/p6xtpLIZRiGSDmcbaao7Re3HPoryYkrayOsfKYiOhpHh8/A8n+HQjf2OA6W3QbgxiivxMBj0HbNRfLExugxoS18/FzUcVlKjvye6fWG98xSjZoQhjA5fbphAkvWVhLZxDBEyuLoMhKu+OKWur+SnLiyNsIYUsvKgOXLcf3ydez/bDcyd4bgzTMPIO8/TU2H+qECdzbaB23/q0ieFo8W/VsBaOW665IDR37PWFtJ5FLsM0Tey9oXt6MzZCtxcVF3cKRPC1BzYLUQUs8gGtPwNlbD8JmEoBDDmu+HVqvBPTPboWHzEHg1R37PvvqKCyWTKnjq+9vHbWcmcoZeb5jgbsUKw7/OLutiqxkBMDQj2HN+Y98W4EaQMpLZ4qJuZW9txPz5hi/1AQMMX9gDBhh+Tk+HqBI4OfMdiPvuQ9Ut54rCWXyD0fgo7hX88vpuXCipi7STffDQ23d4fxACHPs9Y20lkWsJFSgsLBQARGFhodRFoZqsWiVETIwQhq9VwyMmxrDdUZmZ5uex9sjMrF35YmOdK58SpaXZ955aeFT97/GEZok4jRiht3acRmN4Tysrpb5a6djze1ZZaThGo7H8nvN9JC/hqe9v9hkieXD1yBh3dHp2tL+St6lFLYMGQBU0mCf+gQhcsH4c+7rY93sm5ehKIi/EMETSc8fIGHc1I0i4kKDkbCxJIWAIPdb4QNQYhMyoZWSeNfb8nnl6oWQiL8Y+Q2rgqn447uLIyBh7Gb+4axIb67q1pNSghj4tVdDApSMx2NfFPikpwMmThgVi09IM/2ZnMwgROYhhyNulp1vtzCoJS8HMHU1avr62R9E89BCbERxUfMdgbNO+iks+jc22n0EM5uFF+07SpEn1DsJGGg1DqqOMtUhjxhj+5e80AfL/I1hmGIa8mbEfzq21LsZ+OJ4KRMab8plnDH/x3xrMjh2z7zyO1Bbo9Yb/CdRk5Ur+D8IOZ7afw/tjNmFY+HaERfqj97ez0ER/Dv2RicfwAeZHvI3fp32BGX8+baiNsxV03nvvxs+37gfY14WotuT2R7ASuLV7tkyocjSZcbSJtRE+nhptYmlkzK3l0GiECAtz7cgYd4wmU4kqfZXYlXZIzO2fKbrUOVjtLWvpd1LM6JopMt/cLcpLys2fvGrVjc/U0udsHBGl9pF5RO5ivAet/b9WYfcYR5NR7chhhlprI8RuLYexRsD4364YGcMlNBxSVlSGzLf3QbfyGjIO34Yz+gQACQAADarQu8F+aPtcgnZyLBKGtYDGp7nlE9nbqVftI/OI3IHLtDiNYchbSR0GaropbyUEcPEicN99wMaNQEHBjX3OjozxxknpXLxaesGRi/h+4SHovvfFj7mJKEY30766KMGQyH3QDqvEsGdaI7xdB/tPbG/QUfPIPCJ3kMMfwQrFMOStpA4Dtm5KS1atMvzbpAnw8MOGL1Rnv/BtDAM3LW2glI66Lloc9cgPf0K3+DR0vzbClqJEVKGvaV+Uzzkktz4G7YN1cNfT7RHUsJfz5WXQIfI8qf8IVjCGIW8ldRiozc1WUGD4kq9NzYc3TUpXiwkpK0srseWDA9B9chkZe5vjaEULAC1M+zvVOQxt9zxoH49Al7EJ0PgoqKaMiMxJ/UewgnGhVm9m/BIFLIcBR2d1dkRWlmEEg7NctfippRqV2Fj3T0rnqiYtRxZH/d/5r569ih8X7odujR7fZbfFJRFqOtwf5bgrbC+S7ypB8jO3oVnvaCcujohkyQsXlfbU9zfDkLeTMgzUdFPaKzOz9s0tLu5rY5OLmrQA2B0qzy/+Ct9sDIduQz1kXWyPcgSa9oVqLmF4/EFoR/li8DPtEByjsntASp7+3SOS8o9gN2AYciFVhyFAuv8hW7spHZGWZnvyRDmx1qTl7P+IVqwwzBNiwxikYSVuvE+t/LMxsuMpaMc3Qu/H2sEviC3iHufKUEzkCKn+CHYDhiEXUn0YkpKlm7JJE6BfP0MwsMUVNUOe4kSTlk121gzdhZ9RGdwYyX0vQzulGVoPbWHzOR6nploSV4diIkd5yf3GMORCDEMSs3RTAl7Xtm13PykHAl7+3jzU690Bda5dsDhdfBWA6/XCcf2/+9C4bbgDhXVCbf7nqqZaEneEYiKV8tT3N5fjUAOp16ixtHZSDYt+Km60l5G9I+hyc63uElUCh9aewGtDs9AneC+adgzHuGtLYJj60Px9EhoNfDQa1Pv0ffcHodpM7y+XZWE8xR0LDxORWykiDOXm5uKRRx5BWFgY6tati06dOmHnzp3ufVGpA4SryHmNGuNsxdG3jGiKiVFmM4K9w1Wfecbs/a8srUTWoj2Y0TULtwedQtvklpi1rj+2XO0AAR+cqtMGX7ebi8rQCLPTaDz1PtUmzNiaERcwzIir1PvLEs71QqQ4sm8mu3z5Mjp37owBAwbgqaeeQnh4OE6cOIG4uDi0bNnSrnM4XM3mLVX6Sum34CVt23aPoNNoIABsGf5PvL+/H74/1RaXRSPT7gCU4e4me6EdeB0jnmmFmO6RN87v6feptk0+L78MvPSS7ddRUt8wW9zQXEqkVuwz9D+zZs3Cb7/9hs21qFJ26M1USoCwhf0WpGHPemwAqqDBGcQgHtmogi8aawowvMVhaFP8MHhGIuo3re+hAttQmy/29HTDEiv2UNqowZp44VwvRFJhn6H/0el06NatG+6//36Eh4ejc+fO+OCDD2p8TllZGYqKiswedvGmKn32W5CGsemvceMaD/OBQDPk4L2Ed/Dre3uRV9oIy4/3RcqCXvIJQoDzTT7Ge8le3jYj7uOPWw9CgPL6wxF5OdmHoT///BPvv/8+WrVqhR9//BFPPvkkpk2bhk8//dTqc1JTUxESEmJ6xMbG2vdi3hQg2G9BEtcvXcfaPTH40P8Ju45/4sUI9HmqA3wDZPrF6Oz0/o6sTefra744ryPk1rfP2EfPWtOgUvvDEXk7IXP+/v6id+/eZtumTp0qevXqZfU5paWlorCw0PTIyckRAERhYWHNL5aWJoQh8tT8SEtzxaW5V2amfdeSmSl1SRUvb1+++HD8JqFtuk3UQYkAhOiHTO94/ysrhYiJEUKjsVx+jUaI2FjDcTez9166+TyrVjlWtlWrDGW7+TwxMY6fx1VWrbL+PgFCzJtX/X0iohoVFhba9/1dS7KvGYqMjETbtm3NtrVp0wanT5+2+pzAwEAEBwebPex8MdceJyXjQq23Dls30mgMM5JKuWq73P6qt5OoEjjw7XGkDslC7wb7ENk+DI99kgRdXk9cR1008z2DDolVKA0Jh5Dz+28PZ6dAcOYecaQJWm7D9WtqYgcM79WHH3q2TERkN9mHoT59+uDIkSNm244ePYrmzZu7/sWUECDsJfd5fOQ85N+CimsV2LBwN6Z33oiWgTlIHHUb/v5Tf2wrbg8BH3SvdwAv35WFPV8ewcnyaLy97y4Effy+YWYgOb7/jnBmCgRb99KtHGmClmPfPm9qYidSI7fWO7nAf//7X+Hn5yfmz58vjh07Jr744gtRt25d8fnnn9t9Doeq2YxV3bdWdxu3SVUF7yxLTQmxsdJeh7XmBJm9x5dPXhFpU34TY5r/KkJwxayogbguhof/Lv7z8EaRu/Oc9ZPI8f13VmWloVkvLc3wr60mH2v3Um2boOXYBOxNTexEMuKpZjLZhyEhhMjIyBCJiYkiMDBQJCQkiKVLlzr0fIffTG/6AhPC8S8xd5fl1vfWnj4oHnIi85RYdG+WuKvRTuGHcrOiNdHki7+02iRWz9omis8X239SOb3/nmbpXqptgJFj8JBjQCPyAp4KQ7KfZ8gVnJqnwFsmAnSEJ65ZZhPSVVVW4b/LD0L3cQF0u6JxoKyV2f62gceh7XwG2vEN0aPlRfgW5Kvn98FVjH3DHngAuHTJ8jGOzL0js98hAJxbiMhNPDXPkJ/bzqx0xvW01MJTs27LYMj/tYJr+PnNfcj4+jouHr+MQFGKc4jEIbSALypxZ8N90PYvQvLUOLS86zYgfS/w9Djlz0guFV9f4O67gQ8+MHRuBswDg6N9qIz9kWwFD0/27TP20Rs92vD6tbk+IvI8t9Y7yYSnqtkUy5N9eCRqTji7O08sHbdRjAj/XQThmrgXq8RpmDffFNdrIq6+96n5ExXSv0kxXNUELde+fd7WxE4kMTaTuZCnqtkUydPLdnioOUFUCexLP4aM/5yFbmtj/Lck0bTvXqTjG4yGBsJ8Hfhbl1zx5HujpmZZV12rpdrM2FhDDYyUNXZq+iyJ3Ixrk7kQw1ANpOh/YZwjBrDcnODkDL3lxeXY9N5+6L64Ct2BljiljzHb37PefozsfR4zdz8C/4t5sDjo++aAs3mzY++Ns1+C3rIwsBQYPIi8GvsMkWdI0YfHOG+NpQDg4F/1l05cxg8LD0K3VoN1Oe1QhC6mfUG4jkERe6EdUo7hz9yOyE6JQFYBMCDP+glvng/GkffG2UBjbWFX4+SBXLqhZmrr20dEbsEwpHZSzbqdkgKMHOnUX/XHfzmFjHeyodsYgs1X2kOPPqZ9ET75SG51GNr7g3D304mo27in+ZMdCTj2XvOxY8DcuY4HGluTB2o0hskDR45kbQcRkRuxmcyb2dOEoIAhwfpyPX5fdhC6ZReh2x2LQ+Utzfa3DzoKbdezSJ7YGN3Ht4WPXw0TqzvSLJiUZPu9Mc7K7Ey/IjkOEScikhE2k1Ht2NtsI9MhwcV5xVj/5n5krK7A2uMJuCDam/b5oQL9Gu2DdsBVJD/dAvF33g7gdvtO7MiwbHvem8cft75COWDe7HZroJHBNANERKSAtcnICY4uYunM2lNukLvjHP7z8CYMD9+OxpF+SFnQC8uOJeGCaIKGmisY2/w3rJy2BRdOXsPPl7pg2qp+iL8z1rEXcXTNNlvvTSvzSRqtshRovGlhYCIiBWMzmbepzXBwD4/MEVUCf3x9FLr/nIPu93DsvNbWbH8Lv1PQtj8J7SPB6PtkIvzr+rvuxR0dlm3tvalNU5cCmiiJiKTEofUupKowJPN+KGVFZch6Zx8yviyB7uBtyNHfqHHRoAq96h+Ats9FaCfFoM2IltD42LnquTNcEf5qG2jcNM0AEZE3YJ8hco4M+6FcPHYJ3y88CN1aX6zLTUQxupn21UUJBkfug/aeCgyfmYDwdu1rOJOLuWJYdm37XLlwmgEiInIOw5C3cUU/FBfUmBz9MRu6d04h49eG+LWwParQ98ZL++Qh+faj0D5YB3dNS0Sd0F4OnVt2ahtoajHNAABOPEhEVEtsJvM2rmi2cWLyQH25Hls/PADd8kvQ/dEMR8pbmO3vGHQE2u7noH0sHF3GJtQ8/F2ppAglnL2aiLwY+wy5kKrCEGD4grzvvurbbfVDsTYbspXnXT17FT+9uR+6dD2+y26DiyLMtM8f5egfug/au4qRPL0lmvcxXxqDXMDBz4uISGnYZ4hqJywMuHjRfFtoKLB0aa1mQz4T1QMZb/0J3c91sKGgA8rR23RYI81lDI87AO0oXwyZ0Q7BMV1dfFFkwtmriYhchmHI21irLQCAS5esP2/+fOvD8QHT5IGP9D6Ojehv2nyb/0loO5yEdlxD9HkiEX5BfW80F21mHxa32bzZrs/L4mSPRERkhmHIm9RUW2BkqbYgPb3mWZRvEoVc9GmwF9q+l6Cd0gyt74mHxifO/Fzsw+J+rh41yE7YRKRiDEPexJnagv8FKAHAnhl93lteHw3Hd7C8kyuwe44rZ69mgCUilfPCIT0q5kBtgagSOPz9n1hxxzvAmTN2BSHExqLhIyMs77PVhwUw1Erp9faVkWpmXGPt1iVFjDQaw4zaSUk1n8fRpVuIiLwQw5A3CQ+367D3555H66CTaDO8BXT/jbD//DVNHuhIrZQn6fWGWblXrDD86y1hzNE11ixhgCUiAsAw5D3S04Hx42s8pAoanEYsphydimMV8QhAGWKDC+07/7x5NTeZyHDma6SnG+ZcGjAAGDvW8G9cnPfUdtR2gV25BlgiIg9jGFI6vR54+WXDvEK5uVYPq/pfQ9g/8AoeabEV3zy7FQW55Vhw6fGam1sAw/45c2ouh9xWYFdL809KCnDypGGtubQ0w7/Z2fb19ZFjgCUikgA7UCuZpY6vVhT5N8a5CbPw8eJH4BtwS9OJrbW13nrL9sgiYx8WWzNf2+rD4ihLo6AAdc3B4+waa3ILsEREEmHNkEKVfbIC4r77IOwIQgDQ8IcVaLN0RvUgBNS+uQVwTR8WR1lrBrNzziTVN/+4qhM2EZHCsWZIQc7vv4DvFh7Gd+t88Fbe3xAFB9Jsfn7N+2u7WKjxHLYWLHXVfDY1DeO3c86kWjf/KH1uHmOAralW0NUBlohIhhiGZExUCRxaewK6985A91sYthW3g0AS+iELMbDeP8gie5o6nG1uuVlNocpV89nYMwrKHrVp/vGWuXnsCbBERF6OC7XKTMW1Cvy6ZD90nxchY18cTlQ2N9vfte5BzI5bgfsO/tO+E9papd5TXLmoaFaWoUnMWbV9T7xxgVSl13IRkVfiqvUuJPcwVHi6EOsWHoDuW4HvT7fDFdHQtC8Qpbi7yT5oB13HiGdaIbpbpONhYN48oFUr6b7k9HpDXx5r/XgcDScrVhj6CNnDWvOPs4HF1ddCRERWcdV6L5e9Kcew+ntmA2y83B6VuMO0r4nmAoa3PAxtij8GPZOI+k27mz/Z1sgto7Aww78396GRoinH1YuK2tu8NW8e8MEHrm3+4QKpREReh2HIQ6oqq7D9k4PIWFYA3c4o7Cu9HUCsaX+bgBPQds6BdmIYek5sC9+AGkbw1NTx1ejBB4GvvpLHOmGuns/G3mH8c+YYHq5s/uHcPEREXodhyI2uFVzDL2/th+7rUmQcS8D5qkTTPl9UIqnhPiTfWYjkKXFoNaglgJb2n9xax9fYWGDhQmDGDPnMs+Pq+WwcHQXlyhoazs1DROR12GfIxfL25mPtwiPQ/RiAn8+3x3XUNe1rgCIMjd0P7QiBoTPbIrRlo9q/oKWOr5s329enKDPTM005xn42tmpyHO1nY2lEV2yse0dBuetaiIioGvYZUghRJXDg2+PQvZ8L3ZbG+L0kEcCNBVOb+56Btt0JaB9ugDsnJSKg/h3WT+YMS8Ph5daU4675bFwxN5KjODcPEZHXYRhyQsW1Cmx6dx90X1yF7kALnKxsBaCVaX+PevuR3KsA2iej0D6lFTQ+MZ4toBybctw1n40r5kZyFOfmISLyKmwms9Pl7Cv44fUDyFirwQ+n26EQIaZ9QbiOgeH7oB1SihEzWiOyU4Sriu4cOTfleNN8Nt50LUREMsRmMhk4seEUMt45CV1WMDZdaQ89+pj2hWsuILnVYSTfF4CB0xNRL7yHhCW9hZybcqSoyXEXb7oWIiIVYxi6ib5cj/8uPwjdsovI2B2NA2WtANyYAbpd4DFou+RC+5fG6DGhLXz8ZLyAJZtyiIiI7KK4ZrLU1FT8/e9/x9NPP41FixbZ9ZyaqtlK8kvw86L90H1TjrXHE5Avmpj2+aIS/Rrthbb/VSRPi0eL/s1ceSmewaYcIiJSKDaTWbB9+3YsXboUHTp0qNV5zu7Kw9o3jkK3Pgi/5LdHKXqa9oWgEEObHYBWC9wzoy0axXepbbGlpdSmHE+HOIZGIiLVUkwYKi4uxsMPP4wPPvgA//ynnYuU3mKBdjN+2hGH7SXtADQ1bW/pm43pcd9iQK/raP1Id/gNGsAvQks8FRg8vSK8t6xAT0RETlFMM9n48eMRGhqKN998E/3790enTp2sNpOVlZWhrKzM9HNRURFiY2MBFAIIhgZV6Fn/AJJ7X8QjnfcjNu01aPhFWDNPBQZPrwjvjSvQExF5CU81k/m47cwutHLlSuzatQupqal2HZ+amoqQkBDTwxCEgKHhO/Dh+M04+0cBtl5tj78/eQnN/j3NPAgBN9bvSk939aUokzEwuPt90usNgcvaMiKAYRkRvV6Zr0dERLIk+5qhnJwcdOvWDT/99BM6duwIAE7XDJklS+NcPNZWIOeyCgaefJ+ysjy7jIinX4+IiBzCmqH/2blzJ/Lz89G1a1f4+fnBz88PGzduxNtvvw0/Pz/oLfzVHhgYiODgYLNHNZs3W/+CBww1Azk5huPcRa83fCGvWGH4V441EJ58nzy9jIjcli0hIiJJyL4D9d133419+/aZbZs4cSISEhLw/PPPw9fZ2gipvwiV0mnXk++Tp5cRkeOyJURE5HGyD0MNGjRAYmKi2bZ69eohLCys2nabvv4aaNnSMApKyi9Ca512jX1w5NRp15PvU1KSIRDaWkYkyUWTXXr69YiISJZk30zmUo89ZugjEhcHXLhg+KIzjhq6lUYDxMa6/otQaZ12jYHBE++TcRkR43lvfR3AtcuIePr1iIhIlhQZhrKysuyefdqi3FzgwQeBMWMMP3vyi1AOfZUc4enAYFxGJDrafHtMjHtqzDz9ekREJDuKDEO1ZqyBWbkS+OorICrKfH90tPu+CKXuq+QMKQLKyZOGUVxpaYZ/s7PdF0w8/XpERCQrsu8z5DbGGpiDB603AbmDUjvtpqQAI0d6bskKTy8jotRlS4iIqNZkP8+QK5jmKQBg9ywFGo17aj2M8/bY6rSr9vmNiIhI9TjPkNSEAP7v/1zbkdm4tpdxJBk77RIREUlOvc1k9rh4EZg/H3jxxdqfy9K8Qj4+5mErJsYQhNxRG2WteYurtRMRkcqps5lMo7HcRGVJaCiQn1+7gGBtXiGj6dMN/XHcEURqmtwRUMbEj0REpEqeaiZTZxiKjTXMOfTSS/adoDZrU0m5BlpNK7Jb+9i5WjsREckE+wy5w4cf3hg2PWeOodbHHrUZ5i7VvEL2TO5orTyAvCZ+JCIiciN1haH77zfU8Pj6Gh5PP23f82ozzF2qeYVshbCayG3iRyIiIjdSVxi61Zw5QFiY9f2uWGpCqnmFXBGu5DTxIxERkZuoOwz5+gJLl1re56ph7p5c2+tmrghXcpv4kYiIyA3UHYYAQyfhVasMgeVmrlpqQqrFQG2FsJq4K6ARERHJEMMQ4P61qaRYDNSeEFbTPk78SEREKqGuofVuHppnkxQTHFqaZyg21hB2AOv7OKyeiIgkxnmGXEg2YUgqnIGaiIgUyFPf31yOQw1qWpGdq7UTEZHKsc8QERERqRrDEBEREakawxARERGpGsMQERERqRrDEBEREakawxARERGpGsMQERERqRrDEBEREakawxARERGpGsMQERERqRrDEBEREakawxARERGpGsMQERERqRrDEBEREakawxARERGpGsMQERERqRrDEBEREakawxARERGpGsMQERERqRrDEBEREakawxARERGpGsMQERERqZrsw1Bqaiq6d++OBg0aIDw8HKNGjcKRI0ekLhYRERF5CdmHoY0bN2Ly5MnYtm0b1q9fj8rKSgwePBglJSVSF42IiIi8gEYIIaQuhCMuXLiA8PBwbNy4EXfeeaddzykqKkJISAgKCwsRHBzs5hISERGRK3jq+9vPbWd2k8LCQgBAaGio1WPKyspQVlZm+rmoqMjt5SIiIiJlkn0z2c2EEJgxYwb69u2LxMREq8elpqYiJCTE9IiNjfVgKYmIiEhJFNVMNnnyZHz33Xf49ddfERMTY/U4SzVDsbGxbCYjIiJSEDaT3WLq1KnQ6XTYtGlTjUEIAAIDAxEYGOihkhEREZGSyT4MCSEwdepUrF69GllZWYiPj5e6SERERORFZB+GJk+ejLS0NHz77bdo0KAB8vLyAAAhISGoU6eOxKUjIiIipZN9nyGNRmNx+7JlyzBhwgS7zsGh9URERMrDPkP/I/OsRkRERAqnqKH1RERERK7GMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqxjBEREREqsYwRERERKrGMERERESqppgw9N577yE+Ph5BQUHo2rUrNm/eLHWRiIiIyAsoIgx9+eWXmD59OubMmYPdu3cjKSkJQ4cOxenTp6UuGhERESmcRgghpC6ELT179kSXLl3w/vvvm7a1adMGo0aNQmpqqs3nFxUVISQkBIWFhQgODnZnUYmIiMhFPPX9LfuaofLycuzcuRODBw822z548GBs2bJFolIRERGRt/CTugC2FBQUQK/XIyIiwmx7REQE8vLyLD6nrKwMZWVlpp8LCwsBGBImERERKYPxe9vdjViyD0NGGo3G7GchRLVtRqmpqZg3b1617bGxsW4pGxEREbnPxYsXERIS4rbzyz4MNW7cGL6+vtVqgfLz86vVFhnNnj0bM2bMMP185coVNG/eHKdPn3brmyk3RUVFiI2NRU5Ojqr6SvG6ed1qwOvmdatBYWEhmjVrhtDQULe+juzDUEBAALp27Yr169fj3nvvNW1fv349Ro4cafE5gYGBCAwMrLY9JCREVb9ERsHBwbxuFeF1qwuvW13Uet0+Pu7t4iz7MAQAM2bMwLhx49CtWzf07t0bS5cuxenTp/Hkk09KXTQiIiJSOEWEoQcffBAXL17Eyy+/jHPnziExMRHff/89mjdvLnXRiIiISOEUEYYAYNKkSZg0aZJTzw0MDMRLL71ksenMm/G6ed1qwOvmdasBr9u9162ISReJiIiI3EX2ky4SERERuRPDEBEREakawxARERGpGsMQERERqZoiw9B7772H+Ph4BAUFoWvXrti8eXONx2/cuBFdu3ZFUFAQWrRogSVLllQ7ZtWqVWjbti0CAwPRtm1brF692l3Fd5oj152eno5BgwahSZMmCA4ORu/evfHjjz+aHbN8+XJoNJpqj9LSUndfikMcue6srCyL13T48GGz47zt854wYYLF627Xrp3pGCV83ps2bUJycjKioqKg0WiwZs0am8/xhvvb0ev2lvvb0ev2lvvb0ev2lvs7NTUV3bt3R4MGDRAeHo5Ro0bhyJEjNp/niXtccWHoyy+/xPTp0zFnzhzs3r0bSUlJGDp0KE6fPm3x+OzsbAwbNgxJSUnYvXs3/v73v2PatGlYtWqV6ZitW7fiwQcfxLhx4/DHH39g3LhxeOCBB/D777976rJscvS6N23ahEGDBuH777/Hzp07MWDAACQnJ2P37t1mxwUHB+PcuXNmj6CgIE9ckl0cvW6jI0eOmF1Tq1atTPu88fN+6623zK43JycHoaGhuP/++82Ok/vnXVJSgo4dO2Lx4sV2He8t97ej1+0t97ej122k9Pvb0ev2lvt748aNmDx5MrZt24b169ejsrISgwcPRklJidXneOweFwrTo0cP8eSTT5ptS0hIELNmzbJ4/HPPPScSEhLMtj3xxBOiV69epp8feOABcc8995gdM2TIEPHQQw+5qNS15+h1W9K2bVsxb94808/Lli0TISEhriqiWzh63ZmZmQKAuHz5stVzquHzXr16tdBoNOLkyZOmbUr4vG8GQKxevbrGY7zl/r6ZPddtiRLv75vZc93ecn/fzJnP2xvubyGEyM/PFwDExo0brR7jqXtcUTVD5eXl2LlzJwYPHmy2ffDgwdiyZYvF52zdurXa8UOGDMGOHTtQUVFR4zHWzulpzlz3raqqqnD16tVqi90VFxejefPmiImJwYgRI6r9ZSml2lx3586dERkZibvvvhuZmZlm+9TweX/00UcYOHBgtVna5fx5O8Mb7m9XUOL9XRtKvr9dwVvu78LCQgCocRFWT93jigpDBQUF0Ov11Varj4iIqLaqvVFeXp7F4ysrK1FQUFDjMdbO6WnOXPetFi5ciJKSEjzwwAOmbQkJCVi+fDl0Oh1WrFiBoKAg9OnTB8eOHXNp+Z3lzHVHRkZi6dKlWLVqFdLT09G6dWvcfffd2LRpk+kYb/+8z507hx9++AGPPfaY2Xa5f97O8Ib72xWUeH87wxvu79rylvtbCIEZM2agb9++SExMtHqcp+5xxSzHcTONRmP2sxCi2jZbx9+63dFzSsHZMq5YsQJz587Ft99+i/DwcNP2Xr16oVevXqaf+/Tpgy5duuCdd97B22+/7bqC15Ij1926dWu0bt3a9HPv3r2Rk5OD119/HXfeeadT55SKs2Vcvnw5GjZsiFGjRpltV8rn7Shvub+dpfT72xHedH87y1vu7ylTpmDv3r349ddfbR7riXtcUTVDjRs3hq+vb7W0l5+fXy0VGjVt2tTi8X5+fggLC6vxGGvn9DRnrtvoyy+/xF//+ld89dVXGDhwYI3H+vj4oHv37rL5S6I2132zXr16mV2TN3/eQgh8/PHHGDduHAICAmo8Vm6ftzO84f6uDSXf366itPu7Nrzl/p46dSp0Oh0yMzMRExNT47GeuscVFYYCAgLQtWtXrF+/3mz7+vXrcccdd1h8Tu/evasd/9NPP6Fbt27w9/ev8Rhr5/Q0Z64bMPzFOGHCBKSlpWH48OE2X0cIgT179iAyMrLWZXYFZ6/7Vrt37za7Jm/9vAHDaI3jx4/jr3/9q83Xkdvn7QxvuL+dpfT721WUdn/XhtLvbyEEpkyZgvT0dGzYsAHx8fE2n+Oxe9zurtYysXLlSuHv7y8++ugjcfDgQTF9+nRRr149U6/6WbNmiXHjxpmO//PPP0XdunXFM888Iw4ePCg++ugj4e/vL7755hvTMb/99pvw9fUV//rXv8ShQ4fEv/71L+Hn5ye2bdvm8euzxtHrTktLE35+fuLdd98V586dMz2uXLliOmbu3Lli3bp14sSJE2L37t1i4sSJws/PT/z+++8evz5rHL3uN998U6xevVocPXpU7N+/X8yaNUsAEKtWrTId442ft9EjjzwievbsafGcSvi8r169Knbv3i12794tAIg33nhD7N69W5w6dUoI4b33t6PX7S33t6PX7S33t6PXbaT0+/upp54SISEhIisry+z39tq1a6ZjpLrHFReGhBDi3XffFc2bNxcBAQGiS5cuZsPyxo8fL/r162d2fFZWlujcubMICAgQcXFx4v333692zq+//lq0bt1a+Pv7i4SEBLObSy4cue5+/foJANUe48ePNx0zffp00axZMxEQECCaNGkiBg8eLLZs2eLBK7KPI9f92muviZYtW4qgoCDRqFEj0bdvX/Hdd99VO6e3fd5CCHHlyhVRp04dsXTpUovnU8LnbRw6be331lvvb0ev21vub0ev21vub2d+z73h/rZ0zQDEsmXLTMdIdY9r/ldAIiIiIlVSVJ8hIiIiIldjGCIiIiJVYxgiIiIiVWMYIiIiIlVjGCIiIiJVYxgiIiIiVWMYIiIiIlVjGCIiIiJVYxgiIiIiVWMYIiIiIlVjGCIiRVqxYgWCgoKQm5tr2vbYY4+hQ4cOKCwslLBkRKQ0XJuMiBRJCIFOnTohKSkJixcvxrx58/Dhhx9i27ZtiI6Olrp4RKQgflIXgIjIGRqNBvPnz8fo0aMRFRWFt956C5s3b2YQIiKHsWaIiBStS5cuOHDgAH766Sf069dP6uIQkQKxzxARKdaPP/6Iw4cPQ6/XIyIiQuriEJFCsWaIiBRp165d6N+/P959912sXLkSdevWxddffy11sYhIgdhniIgU5+TJkxg+fDhmzZqFcePGoW3btujevTt27tyJrl27Sl08IlIY1gwRkaJcunQJffr0wZ133on//Oc/pu0jR45EWVkZ1q1bJ2HpiEiJGIaIiIhI1diBmoiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVI1hiIiIiFSNYYiIiIhUjWGIiIiIVO3/ATyZ5axZC3zBAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "%matplotlib inline\n", "\n",