{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Rješenja zadataka iz predmeta Akvizicija i obrada eksperimentalnih podataka\n", "#### Novembar 2018" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### ZADATAK 1\n", "Koristeći Newtonov metod naći korjen jednačine $x^4-2x^3+5x^2-6=0$ na intervalu [1,2]. Newtonov metod koristi rekurzivnu formulu $x_{n+1}=x_n-\\frac{f\\left(x_n\\right)}{f'\\left(x_n\\right)}$, gdje se u prvoj iteraciji n=1 uzima neki proizvoljan broj iz zadatog intervala. Uporediti dobijene vrijednosti za 3, 5 i 50 iteracija do na 15. decimalu (prikazati u tabeli). Nacrtati funkciju u datom intervalu." ] }, { "cell_type": "code", "execution_count": 89, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Unijeti broj iteracija: 5\n", "Unijeti pocetnu vrijednost x: 1.1\n", "1.21756215475076\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEPCAYAAABBUX+lAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XdUFPf+P/4nu9KrIB0UK0VBFIzEq0awKyiaWJJojDGaxHKjKf6MKZZ4kktivmpiQdOjxnhNoqjYo7EXNCiIFAsg0pvSy+7O7w+vfGJsLAw7C/N8nMM5LMy89/XenX0yvOc9MwaCIAggIqIWTyF1AUREpBsMfCIimWDgExHJBAOfiEgmGPhERDLBwCcikgkGPhGRTDDwW6jdu3cjKChI0hoKCgowceJETJo0CS+99BLy8vIkrSc2NhYTJkzApEmT8NZbb6G2tlbSeoh0zYAnXrU8Go0Gb775Jm7duoXt27dLVodarYaBgQEUCgV+//135OTkYObMmZLVk5ubC2tra5iYmGDFihXw9vbGsGHDJKuHSNe4h99AISEhOHXqlFbrjBw5EmfPnn3ge7Ht2rULQ4cOhYGBQZO0X19KpRIKxd1NrLy8HJ07d5a0HkdHR5iYmAC4vza6X1Num2JryOcQADw9PeHv748VK1Y0QVXASy+9BF9fXzz//PNN0n5DtfgtPiQkBH5+fujRo0fdV25uriS1REdHo3fv3g98Lya1Wo29e/dixIgRorRXU1ODhQsXIjg4GD169EB4eDiOHj1a7/UTExMxbtw4bNq0CT4+PqLUNHnyZPj6+ta9n0OHDtVq/YyMDBw7dgwDBgwQpR7g7vs5fPhw+Pv7Y9CgQTh//nyj27x9+zZmzZoFf39/BAcHY9euXSJU+mTabJsNDVx9EBUVhXnz5tU93rRpE8aOHYtu3bphwYIFj1yvPp+Jn376CUuWLGmy2huqldQF6EJkZCT69OkjdRk6sXPnTgwfPly0vVeVSgVnZ2ds3LgRLi4uOHr0KObOnYtdu3bBzc3tiet7e3tj27Zt2LNnD9avX4+lS5c+ctmvvvoKADBnzpwntvvRRx9h3Lhx9e/I/5SVlWHBggX47LPPYGRkpPX6D3Py5EksX74cK1asgJ+fH/Lz80Vpd+nSpTA0NMTJkyeRmJiI1157DV5eXpL/pyQWlUqFVq30J4IcHBwwc+ZMHD9+HNXV1Y9crrGfCSm1+D38R/H09ER6enrd4wULFtz3711ISAi+/fZbhIWFISAgAHPnzn3kRnD9+nWEhIQgOjoaAJCdnY3Zs2cjKCgIvXv3rgu5v+8N/XPPaMOGDRg0aBB69OiBESNG4ODBg4+s/bPPPsOsWbPqHkdERGDKlCmora3FtWvXEBUVhWnTpiE9PR3Lli1rwKvzf8zMzDBnzhy4ublBoVAgODgYbm5uSEhIeGI9NTU1dT+ztLSEqalpo2qpj8e9NiqVCm+99RZmz56NDh06iPacX331FWbOnAl/f38oFAo4OjrC0dGxUW1WVFTgwIEDePPNN2Fubo7AwECEhIQgKirqocs/aXvVZvv653b6qHbfffddZGVl4fXXX0ePHj3w9ddfA7h7rGTOnDkICgpCSEgIfvrpp/va3rBhA8LCwuDv74+1a9fi3//+933Pv2zZsrrtVpu6n/SZfpIhQ4Zg0KBBsLGxeexy9flM6Cv9+fOqh/bu3YtvvvkGxsbGeP755/H7778/MCaXkJCAWbNmYdGiRQgODoZarcZrr72GoKAgHD58GEqlEvHx8U98Lnd3d2zevBn29vbYt28f3n33XRw4cAAODg4PLDt9+nQMHjwYiYmJuHjxIk6cOIGff/4ZhoaGePfdd+uWGzt2LD744IPGvxB/U1BQgLS0NHTq1OmJ9SQkJGD58uVQKBQwNjbGJ598IlodX3zxBZYvX4727dtj3rx5dUMQj3ttduzYgUuXLmHt2rVYu3Ytnn/++UYPfanValy+fBkhISEYPHgwqqurMWjQIMyfP7/ueME9r732Gi5cuPDQdgICArB+/fq6x2lpaVAoFGjfvn3dz7y8vBATE/PIWh63vWqzfdW33c8//xwXLlzAsmXL6v6D1mg0eOONNxASEoIvvvgCubm5ePnll9G+fXv069cPwN0how0bNqB169YoLCzE+vXrUVZWBgsLC6jVauzbtw+rV69udN268rDPhL6SReDPmjULSqUSAPDUU09h7dq19Vpv8uTJdXtqwcHBSExMvO/358+fx6+//orPPvusbgpkXFwc8vLyMH/+/Lp/VwMDA5/4XMOHD6/7fsSIEVi/fj3i4uIwaNCgB5Zt3bo1pkyZgvnz56OsrAw///wzLC0tH1ju999/r1c/66u2thbvvPMOxowZg44dOz6xnh49emDz5s2i1gAA77zzDjp27AgjIyNER0fj9ddfR1RUFNq2bfvY1yY8PBzh4eGi1lJQUIDa2lrs27cPmzdvRqtWrTBz5kysW7fuvvFhAPcF+pNUVFQ88J5aWlqivLz8kes8bnvVZvvSpt1/io+PR1FREWbPng3gbmCPHz8ee/bsqQv8yZMnw9nZGQDg6uoKHx8fHDp0COHh4Thz5gxMTEzg7+/f6Lp14VGfCX0li8Bfs2ZNg8bw7e3t6743NTV9YB75L7/8gl69et033z07OxsuLi5aj03u2LED33//PTIzMwHc/cAXFxc/cnlvb2+sXr0ay5cvr/vw1NfkyZNx7ty5h/6uZ8+e2LJlywM/12g0mD9/PgwNDfHhhx+KVs/f93rvDRX8+OOPAB7c672ne/fudd+PGTMGu3fvxtGjRzF58uRG1QJo/9rc24ufPHly3V7n1KlTHxr42jAzM0NZWdl9PysrK4O5ufkj13nc9qrt9lXfdv8pMzMTeXl59+3kqNXq+x7/8z0JDQ3F7t27ER4ejt27dyM0NFSUupvakz4T+kgWgf8wpqamqKysrHucn5+v9bjrkiVL8PXXX+OTTz7BwoULAdzdmLOzs7U6IJWZmYkPPvgAP/zwA3r06AGlUonRo0c/cvnk5GQsXrwYY8aMwW+//YawsDCt6t64caNWywuCgPfffx8FBQX4+uuvYWhoKFo9fw90bQ7a/p2BgQHunU6i69fG2toaTk5O9ZoC++qrrz52SOebb76pe+zh4QG1Wo20tDR4eHgAAJKSkho0bKDt9tUYzs7OcHNzw4EDBx65zD9fq+HDhyMiIgI5OTk4ePAgtm7d2qC6xfhM19eTPhP6SrYHbb28vLB7926o1WocO3bssWOjj2Jubo5vvvkG58+fx/LlywEAfn5+sLe3xxdffIGKigpUV1c/8kN+T2VlJQwMDGBrawsA+O2333D16tWHLpubm4s33ngDS5YswaJFi5CSktLkc6YXLVqE69evIzIy8oFxaV3XU1JSUjeLQqVSYefOnTh//jz69u0ryWsD3D1WsnHjRhQWFuLOnTv48ccfHzrl85tvvkFsbOxDv/4e9sDdPfzBgwfjyy+/REVFBS5cuIA//vijQUGtzfalrTZt2iAjI6PusZ+fHywsLLBhwwZUVVVBrVYjJSUFcXFxj2zD1tYWTz31FN577z24ubnVDY1oW3djP9MqlQrV1dXQaDRQq9V129g9CxYsqJuu+bjPhD6TbeC///77OHLkCAIDA7Fr164GjwlaWVnhu+++w7Fjx7By5UoolUpERkYiPT0dwcHB6N+/P/bu3fvYNjp16oRXXnkFEydORJ8+fZCSkoKePXs+sFxZWRmmT5+Ol19+GQMHDoSpqSmmTZvWZCePAHf3srZu3YrExET07du3bu77zp07JalHpVJh5cqVCAoKQlBQEDZt2oQ1a9bAwcFB57XcM3PmTPj6+mLo0KEYMWIEfHx88MYbbzS63UWLFqGqqgp9+vTB22+/jcWLFzdoSmZ9t6+GmDFjBtatW4fAwEB8++23UCqVWLduHZKSkjBw4EAEBQXhgw8+eGB46p9CQ0Nx6tSp+4ZztK27sZ/pdevWwc/PDxs2bMDOnTvh5+eHdevW1f0+OzsbPXv2fOxnQt/x0goSGTBgAD7//HP06tVL6lKIZMfX1xdGRkaYPHky5s6d+8Tla2pqMHr0aOzcubNewzdTp07FxYsX4efnV3dMSh8w8CVQVFSEAQMGYN++fXBxcZG6HCKSCdkO6UglLi4OQ4YMwaRJkxj2RKRT3MMnIpIJ7uETEckEA5+ISCYY+EREMsHAJyKSCb25tEJxcTk0Gu2PH9vZWaCw8PEndbQ07LM8sM8tX2P6q1AYoHXrR19b6WH0JvA1GqFBgX9vXblhn+WBfW75dNlfDukQEckEA5+ISCYY+EREMqFV4EdERCAkJASenp5ISUl54PerV69+5O+IiEhaWgX+wIEDsXnzZri6uj7wu4SEBFy8eJHXhyEi0lNaBX5gYOBDbxlXU1ODpUuXYtGiRfW68w8REemeKGP4q1atwqhRo+Du7i5Gc0RELd6JuGzM+PQQqmvVOnvORs/Dj42NRXx8PN55551GtWNnZ9Hgde3tLRv13M0R+ywP7HPLdD4xFz/sS0L3Tm3g6myts5GRRgd+TEwMbty4gYEDBwIAcnJyMG3aNHz66afo27dvvdspLCxr0AkI9vaWyM8v1Xq95ox9lgf2uWVKzS7BZz/Hws3eHAum9EJBQcPPtNV2R7nRgT9jxgzMmDGj7nFISAgiIyPRpUuXxjZNRNSi5N+uxKptl2Bhaoi547rDzMQQ5aVVOnt+rcbwly1bhv79+yMnJwdTp07FyJEjm6ouIqIWpaSiBv9v60WoNQLemtAdNhbGOq9Bb+54xSGd+mOf5YF9bjmqalT4fEssbuWX492JPdDJzRpA4/rbkCEdnmlLRNSEVGoN1u1IQFpOKV4f3bUu7KXAwCciaiKCIODHvUmIv1GIyUM90aOzvaT1MPCJiJrIr0ev4+TlHIz6lwcG+D94hQJdY+ATETWBA+duYu+ZmxjQwxWj+7aXuhwADHwiItGdTsjBL4evIaCLPSYN7qI3l5xh4BMRiSjueiG+i06Ep7sNZozygUKhH2EPMPCJiERz7dYdrN0eD1d7c/z7OT8YtlJKXdJ9GPhERCK4lV+GldsuobWlMd4a7w9TY725ZXgdBj4RUSPl367EF1svwshQgbcn+MPK3Ejqkh6KgU9E1Ai3y6qx/JdYqFQavD3BH21sTKUu6ZEY+EREDVRWWYsvtl5ESXkt5o33h6t9wy/zrgsMfCKiBqiqUWHVtkvILarAnGd90cHFSuqSnoiBT0SkpVqVGl/9Fo/U7FK8NqobfDxspS6pXhj4RERaUKk1iIxKQGJ6MV4Z6YUAT2mvj6MNBj4RUT1pBAHf7UlE7NUCTBrSBX26OUtdklYY+ERE9SAIAjbuT8aZhFw8+0wHhPR0k7okrWkV+BEREQgJCYGnpydSUlIAAMXFxZg+fTqGDh2KsLAwzJ49G0VFRU1SLBGRFARBwNbD13D0YhZGPt0OI5/2kLqkBtEq8AcOHIjNmzfD1fX/LvNpYGCAV199Ffv378euXbvg7u6O5cuXi14oEZFUdhxPxYGYDAwKcMPY/h2kLqfBtAr8wMBAODvfP2ZlY2OD3r171z329/dHVlaWONUREUks+nQadp1KQ18/Z0wc1FlvrnzZEKJe7EGj0WDLli0ICQnRel1t7834d/b2lg1et7lin+WBfZZW1LHr+O3oDTzTww3zXugJZRNc+VKX/RU18D/++GOYmZlh0qRJWq/Lm5jXH/ssD+yztP6MzcRP+5MR4GmPSYM7oaiwTPTn0PVNzEUL/IiICKSnpyMyMhIKBSf/EFHzdfxSFn7anwy/jnZ4bVRXKFtIpokS+CtWrMDly5exYcMGGBnp51XiiIjq4/TlHPywNwld29ti1phuaKVsGWEPaBn4y5Ytw4EDB1BQUICpU6fCxsYGK1euRGRkJDw8PDBx4kQAgJubG9asWdMkBRMRNZVzibn4JvoKvNq1xpyxvnp3A5PGMhAEQfuB8ybAMfz6Y5/lgX3WrfNJeYiMSkAnN2vMG9cdxkZNH/a6HsNvOf+rEBE10IXku2HfwdUKbz7np5OwlwIDn4hk7a+UfERGJaC9iyXmjeuul7cmFAsDn4hkKzYlH+t2XEY7J0u9vQ+tmBj4RCRLsSn5WLvjMto6yiPsAZFPvCIiag7uhf29PXszE3lEoTx6SUT0PxeS8xEZJb+wBxj4RCQj96Zetne2xDyZhT3AwCcimTiXmIsNO6+gg6tVi5+N8yjy6zERyc7phBx8s/sKOrvZYO44P5gYyTP65NlrIpKN43FZ+GFPEjzb2uDN53RzBq2+YuATUYt17xLHXdvbYvZYXxgbyjfsAQY+EbVQB89nYMuhq/DraIdZY7q1uAuhNQQDn4hanOjTafjt6A307GKP10d3bVGXOG4MBj4RtRiCIGDH8VTsOpWGIB9HTAv1bjE3LxEDA5+IWgRBEPDfI9ew/1wG+vo54+VhXlA0wT1omzOt/vRFREQgJCQEnp6eSElJqft5amoqJkyYgKFDh2LChAlIS0sTu04iokfSaAT8tD8Z+89lYGBPN7w8nGH/MFoF/sCBA7F582a4urre9/NFixbhhRdewP79+/HCCy/go48+ErVIIqJHUWs0+Cb6Co5ezMLIp9vhhcGdoTBg2D+MVoEfGBgIZ2fn+35WWFiIK1euIDQ0FAAQGhqKK1euoKioSLwqiYgeolalxtrtl3EmIRfPPtMBzz7TEQYM+0dq9NGM7OxsODo6Qqm8O+VJqVTCwcEB2dnZjS6OiOhRqmpUWLktDrFXC/Di4C4Y+bSH1CXpPb05aKvtvRn/zt7eUsRKmgf2WR7Y54crrahBxM+xuHrrNuY93xMhge46qKxp6PI9bnTgOzs7Izc3F2q1GkqlEmq1Gnl5eQ8M/TwJb2Jef+yzPLDPD1dcWo3/99+LyC2qwMzwbvBtZ9NsX6dmdxNzOzs7eHt7Y/fu3QCA3bt3w9vbG7a2to1tmojoPrnFFfh00wUU3KnCvHHd0bOLvdQlNSta7eEvW7YMBw4cQEFBAaZOnQobGxtER0dj8eLFWLBgAdauXQsrKytEREQ0Vb1EJFMZeWX4YutFaDQC5j/fA+2draQuqdkxEARB+3GUJsAhnfpjn+WBff4/KRm3serXOJgYKfH2BH+4tDGXoDrx6XpIR28O2hIRPUzs1XxERiXAzsoEb0/wh521idQlNVsMfCLSWyfisvHD3iS0c7LA3HHdYWlmJHVJzRoDn4j0jiAI2HMmHb8dvYGuHq0xc4yvLG9JKDa+gkSkVzSCgF8OXcWhC7fQ28cR00Z68/LGImHgE5HeqFWpsWFnAs4l5mFwoDsmDOzE6+KIiIFPRHqhoqoWKzacQfz1Aowb0BHDerfldXFExsAnIskVl1ZjxX8vIqeoAtPDfPB0VyepS2qRGPhEJKnM/DKs2HYJFVUqLHo1CK6tTaUuqcXikRAikkxiejE+2fQX1BoBC17sCf8uDlKX1KJxD5+IJHEmIQffRifC0dYM88Z15wlVOsDAJyKdEgQB0afT8fuxG/B0t8HsZ31hbmIodVmywMAnIp1RqTXYdCAZxy5lI8jHEVNHeMOwFUeWdYWBT0Q6UVmtwrodl3E5tQihfTwwpl97TrvUMQY+ETW5opIqrNx2CdmFFXh5uBf6d3eRuiRZYuATUZNKyynBql/jUFOrxtzx3dHVgzdHkgoDn4iazF8p+diwKwGWpkZ4Z5I/XO0bfu9qajzRAv/IkSNYtWoVBEGARqPBnDlzMGTIELGaJ6JmRBAE7Dt3E78euY72LlaYM9YX1hbGUpcle6IEviAImD9/PjZv3owuXbogKSkJzz//PAYNGgSFgkfgieTk7zNxenk5YNpIbxgZKqUuiyDiHr5CoUBp6d1bdZWWlsLBwYFhTyQzZZW1WLs9Hkk3byO0jwfC+7Xn1S71iCiBb2BggJUrV2LmzJkwMzNDeXk51q9fL0bTRNRMZBeWY9W2OBSVVvMCaHpKlJuYq1QqvPrqq5gzZw4CAgJw4cIFvP3224iOjoa5ecu42TARPdpfyXn47KcYGLZS4v2pT8GLM3H0kih7+ImJicjLy0NAQAAAICAgAKamprh+/Tr8/Pzq1UZhYRk0Gu3/9jTmru/NFfssD82hz4Ig4ND5W/jl8FW4trHAv5/zhZ25YYPrbg59FlNj+qtQGMDOTrtZT6IEvpOTE3JycnDjxg106NAB169fR0FBAdq2bStG80Skh/5+cLZH5zaYHuYDEyPO9NZnorw79vb2WLx4Md588826U6U//fRT2NjYiNE8EemZO+U1WLM9Htdu3eHB2WZEtD/Ho0aNwqhRo8Rqjoj0VHpOKb76PQ5lFbV4fXRXPOXtKHVJVE/8/4uI6u3slVx8vycRFmaGeG9SANo5WUpdEmmBgU9ET6TRCPjt6HXsPXsTnd2sMXOML6zNjaQui7TEwCeixyqvqsX6qARcTi1CcA9XPD+oM1opeVJlc8TAJ6JHupVXhq9+j0NRSTWmDPPEM/6uUpdEjcDAJ6KHOpeYi+/2JMLUuBX+vxd7opOrtdQlUSMx8InoPmqNBr/9eQP7zt1EJzdrzAzvBhte6bJFYOATUZ2S8hpERl1G0s3bCO7piucHcry+JWHgExEA4EZWCdZsj0dZZS2mjfTGv3ydpS6JRMbAJ5I5QRDwZ2wmfj50Fa0tjbGQ8+tbLAY+kYxV16rx075knE7IgV9HO7wa6gMLU0Opy6ImwsAnkqmcogqs3R6PzPxyhPdrj9A+HrweTgvHwCeSofNJefhuTyJaKRWYN747unWwk7ok0gEGPpGMqNQa/PfINRw6fwsdXazwRng32FqZSF0W6QgDn0gmCu5UYt2OBKRml2BQgBvGh3TilEuZYeATycDFawX4dvcVaAQBM8O7IdDLQeqSSAIMfKIWTKXW4Lej17H/XAbaOljgjTHd4NjaTOqySCKiBX51dTU++eQTnD59GsbGxvD398fHH38sVvNEpKWCO5WIjErAjawSBPd0xcSQTjBspZS6LJKQaIH/+eefw9jYGPv374eBgQEKCgrEapqItHQhOR/f70mEAAFvhHdDLw7hEEQK/PLycuzYsQNHjx6tu6dtmzZtxGiaiLRQq1Jj6+FrOPxXJjycLPH66K5w4BAO/Y8ogZ+RkQEbGxusXr0aZ8+ehbm5Od58800EBgaK0TwR1UN2YTnWRyXgZl4ZhvRyx3MDOnIWDt3HQBAEobGNXL58Gc8++yyWL1+OsLAwXLp0Ca+//joOHjwICwsLMeokokcQBAF/xNxE5PZ4GBsqMXdiD/TycZK6LNJDouzhu7i4oFWrVggNDQUAdO/eHa1bt0Zqaip8fX3r1UZhYRk0Gu3/9tjbWyI/v1Tr9Zoz9lke6tPniioVftqfhHOJefBqa4PpYV3R2tK42b5WcnufG9NfhcIAdnba7VCLEvi2trbo3bs3Tp48ib59+yI1NRWFhYVo166dGM0T0UNcu3UH63cmoLi0GmP7d8CIoHZQKHgtHHo00WbpLFmyBAsXLkRERARatWqFzz77DFZWVmI1T0T/o9ZoEH0qHVEnU2FnZYL3JvdERxfefpCeTLTAd3d3x8aNG8VqjogeIv92Jb7edQXXMu8gqKsjJg/xhKkxz5+k+uGWQtQMCIKA0wk52HQgBQYGwIwwHwR15YFZ0g4Dn0jPlVXWYuP+ZMQk5aGzmzWmh/qgjY2p1GVRM8TAJ9JjCWlF+C46ESXlNTwwS43GwCfSQzW1anwdFY+dx27A2c4Mc54NgIcTJ0FQ4zDwifRMek4pNuxKQHZhBQb2dMNzwR1hbMiLnlHjMfCJ9IRao0H06XTsOpkGK3MjLJ3xNNxsOVZP4mHgE+mBrIJyfBt9BanZpQjyccSLQ7rAw91WVmedUtNj4BNJSKMRcOh8Bn47dgPGhkrejYqaFAOfSCK5xRX4LjoRV2/dgX+nNpgyzBPWFsZSl0UtGAOfSMc0goAjf2Vi25/XoFQoMG2kN/p0c6q7lwRRU2HgE+lQXnEFvt+ThOSM2+jWwRYvD/OCrZWJ1GWRTDDwiXRAIwg4fOEWfj16HUqFAaaO8EJfX2fu1ZNOMfCJmlhOUQW+25OIa7fuwLeDHaYM8+RePUmCgU/URNQaDQ6cy8COE6kwVHKsnqTHwCdqAjdzS/H93iSk55SiR+c2mDzUEzacgUMSY+ATiahWpcauU+nYeyYd5iatMDO8GwI87blXT3pB9Fvar169Gp6enkhJSRG7aSK9lpJxG4u+i8HuU2no7eOIZdODEOjlwLAnvSHqHn5CQgIuXrwIFxcXMZsl0msVVbX49c/r+PNiFtpYm+CtCd3Rrb2d1GURPUC0wK+pqcHSpUuxfPlyTJkyRaxmifSWIAi4kJyPzYdSUFJegyG93DGmXwcYG/HKlqSfRAv8VatWYdSoUXB3d2/Q+nZ2Fg1+bnt7ywav21yxz9LKK67Ahu3xOJuQgw6u1lj86tPo5G4j+vPoU591RW591mV/RQn82NhYxMfH45133mlwG4WFZdBoBK3Xs7e3lN0VBdln6ag1Gvxx/ha2H0+FAAHjgzthcC83KBUK0evTlz7rktz63Jj+KhQGWu8oixL4MTExuHHjBgYOHAgAyMnJwbRp0/Dpp5+ib9++YjwFkeRuZJXgp/1JuJlbBr+Odpg0pAvaWPN69dR8iBL4M2bMwIwZM+oeh4SEIDIyEl26dBGjeSJJVVTV4rdjN/DnX5mwtjDiVEtqtjgPn+gRBEHAmSu52Hr4GkorajAo0B3h/drD1JgfG2qemmTLPXz4cFM0S6QzWQXl2HQgGUk3b6O9sxXmjeuOdk7yOphILQ93VYj+pqpGhZ0n03AwJgMmRkq8NMwT/bu7QMHhG2oBGPhEuDt8cz45H7/8cRXFpdXo6+eM5wZ0hJWZkdSlEYmGgU+yl1lQjp8PpiAxvRhtHSzwRng3dHK1lrosItEx8Em2KqpU2HkyFX9cuAUTIyUmDemCAf6uUCg4fEMtEwOfZEcjCDgZn43f/ryO0opa9Ovugmef6QBLDt9QC8fAJ1m5nnUHPx+8itTsEnR0tcKb47qjvbOV1GUR6QQDn2ShuLQav/7FI6AXAAAQ1UlEQVR5HacTcmBtboRXQ70R1NWJs29IVhj41KLV1KpxICYD0afTodZoMCKoHUY+3Y4nT5EscaunFkkQBMQk5WHbkesoLKlCj85tMCGkExxam0ldGpFkGPjU4tzIKsEvf1zFtcw7cHewwCsje8C7XWupyyKSHAOfWoyikir8evQ6ziTkwsrcCFOGeaKfnwunWRL9DwOfmr3KahX2nEnHgZgMCAIw8ul2GBHEcXqif+IngpotlVqD45eysONEKkorahHU1RHP9u8IO2sTqUsj0ksMfGp2BEFA7NUC/PrndeQUVaCLuw3mjuvE+fRET8DAp2blWuYdbDtyDVdv3YGznRnmPOsL/05teDMSonoQJfCLi4sxf/583Lx5E0ZGRmjXrh2WLl0KW1tbMZonwq28Uny9PR5/peTD2twILw31RL/uzlAqFFKXRtRsiBL4BgYGePXVV9G7d28AQEREBJYvX45PPvlEjOZJxopLqxF1IhUn4rNh2EqB8H7tMbRXWxgbKaUujajZESXwbWxs6sIeAPz9/bFlyxYxmiaZKq+qxZ4z6Th0/hY0GgEj/9UeA/1dYGXOC5wRNZToY/gajQZbtmxBSEiI2E2TDFTXqHHwfAb2nr2JqmoVend1xJh+HeDT2QH5+aVSl0fUrBkIgiCI2eCSJUuQm5uL1atXQ8HxVaqnWpUa+8+kY+uhFNwurcZTPk6YNNwL7V14IxIisYga+BEREUhOTkZkZCSMjLT717uwsAwajfal2Ntbym7PryX1Wa3R4NTlHOw8kYbCkip4uttg7DMd0NnN5r7lWlKf64t9bvka01+FwgB2dhZarSPakM6KFStw+fJlbNiwQeuwJ/nRCALOJeYi6kQacosq4OFkiZeHe8HHozWnWBI1EVEC/+rVq4iMjISHhwcmTpwIAHBzc8OaNWvEaJ5aEI0gIDYlHztOpCIzvxyu9uaYPdYXPTpzLj1RUxMl8Dt37ozk5GQxmqIWShAEXLpWiB3Hb+BmXhmcbM0wY5QPnvJ25E1IiHSEZ9pSkxIEAXHXCxF1IhVpOaVwsDG9e7cpHydexZJIxxj41CTuBf3Ok6lIzS5FG2sTTB3hhT7dnHh2LJFEGPgkqntDNztP3t2jb2NtgqnDvfB0Nye0UjLoiaTEwCdR3DsYu+tUGm7mlsHBxhRTR3jh6a4MeiJ9wcCnRtFoBJxLykX06XRk5pfDsbUppo30RlBXRw7dEOkZBj41iEqtwemEHOw5nY7c4ko42/1v1o2XIw/GEukpBj5ppaZWjeNx2dh3Nh2FJdVo62CBmeHd0NPTntMrifQcA5/qpbJahSOxmTgQk4GS8hp0dLXCi0M80b2jHU+YImomGPj0WCXlNTh4PgOH/8pEZbUKXT1aI3R0V3Rxt2HQEzUzDHx6qLzbldh/9iZOxGdDpdIgwNMeI55uBw8n3jeWqLli4NN90nJKsO/sTcQk5UGpMECfbs4Y1rstnGzNpC6NiBqJgU8QBAGXU4uw7+xNJKYXw9RYiWFPtcXgXu6wsTCWujwiEgkDX8ZUag3OXsnFvnM3kZlfDhsLI4wP7oRn/F1gasxNg6il4adahsoqa3H0YiYOXbiFO2U1cLM3x7SR3ujt48izYolaMAa+jOQWV+BQzC0cj89CTa0GXT1aY9oIb3Rtb8sZN0QywMBv4QRBQErGbRyIycDFqwVQKg3Q28cRQ3u1hZuDdrdHI6LmTbTAT01NxYIFC3D79m3Y2NggIiICHh4eYjVPWqpVaXAuMRcHYzJwM68MFqaGCO3jgZCerrDmgVgiWRIt8BctWoQXXngBo0ePRlRUFD766CP89NNPYjVP9XS7rBp/xmbiz9hMlFTUwrWNOV4e7oUgH0cYGSqlLo+IJCRK4BcWFuLKlSv4/vvvAQChoaH4+OOPUVRUBFtbWzGegh5DEATcyCrBHxduISYpDxqNAL+OdhjUyx0+7XhTcCK6S5TAz87OhqOjI5TKu3uQSqUSDg4OyM7OZuA3oVqVBjFJuTh0/hbSckphYqREcE9XDAxwg2NrnihFRPfTm4O2dnYNP4Bob28pYiX6L7+4EntjMnDgbPrdaZUOFnh9jC+CA91hZmIodXlNRm7vM8A+y4Eu+ytK4Ds7OyM3NxdqtRpKpRJqtRp5eXlwdnaudxuFhWXQaAStn9ve3hL5+aVar9fcaAQBienFOHzhFi5dK4AAwL9TG4T0dIOPx91hm/LSKpSXVkldapOQy/v8d+xzy9eY/ioUBlrvKIsS+HZ2dvD29sbu3bsxevRo7N69G97e3hzOEUF5VS1OxmXjyMUs5BZVwMLUEGODO+MpzzZoY20qdXlE1IyINqSzePFiLFiwAGvXroWVlRUiIiLEalp2BEHAjewS/BmbiXOJeahVadDR1QqjQn0Q6OUAF2drWe0FEZE4RAv8jh07Ytu2bWI1J0uV1SqcvZKLP2MzcTOvDMaGSvyrmxMG9HBFW0d5jWsSkfj05qCtnKXllODP2CycvZKL6lo13B0sMHmoJ4J8HHkRMyISDdNEIhVVKpy9koOjl7JwM7cMRq0UeMrbEc/4u6CDixXnzhOR6Bj4OiQIAq5l3sGxS1mIScpDTa0G7g4WeHFwFzzd1QlmJnw7iKjpMGF04E55DU5fzsHxuCxkF1bA2EiJIB8nPOPvAg8nS+7NE5FOMPCbiFqjQfyNIhy/lIW464VQawR0dLXC1BFe6OXlABMjvvREpFtMHZFlFZTjZHw2Tl3OwZ3yGliZGWJwL3f09XWGSxtzqcsjIhlj4IugvKoW5xLzcDI+GzeySqBUGMCvox36+jrDt6Md7yJFRHqBgd9Aao0Gl28U4eTlHFy8WgCVWgNXe3NMCOmEoK5OsDY3krpEIqL7MPC1IAgCbuaW4XRCDs5cyUVJeQ0sTA3xjL8L/uXrhHaOPABLRPqLgV8PRSVVOHslF6cScpCZXw6lwgDdO7VBn25O8OOQDRE1Ewz8R6ioUuFCSh7OJOQiKb0YAoCOLlaYPKQLenk7wsK05V6GmIhaJgb+39SqNIi/UYgzCTm4eK0QKrUGDjamGNW3PYK6OvKmIkTUrMk+8DUaAUk3i3H2Si7OJ+ejsloFS7O74/JBXR3RwZmXOSCilkGWgX/vHrBnr+QiJikPd8prYGykRM/O9ni6qyO8PVpDqeC4PBG1LLIJ/HszbM4l5SImMQ8Fd6rQSmkAv45t0NvHEX4d7WBsqJS6TCKiJtOiA18QBGTml+NcUh5iEnORW1wJhYEBfNq3xui+7dGjsz0vWEZEstHotFuyZAlOnz4NIyMjmJmZ4f3334evr68YtTWIIAjIKihHTFIeYpLykF1YAQMDwKttawzt3RYBXexhacaToohIfhod+P3798fChQthaGiII0eOYN68eTh06JAYtdWbIAi4lV+O80l5OJ/8v5AH4NnWBoMC3NDT04FnvhKR7DU68IODg+u+9/f3R05ODjQaDRQ6OOhZXFqNPecycDz2FnKLK2FgAHi622BggBsCutjD2sK4yWsgImouDARBEMRqbPXq1UhKSsLq1avFavKx1m+Pw55TafDr1AZ9/FwQ1M0JrS1NdPLcRETNzRMDf8yYMcjKynro706dOgWl8u7MlujoaHz55ZfYvHkz2rRpo3UhhYVl0Gi0+9uj0QiwtDZFeWmV1s/XnNnbWyI/v1TqMnSKfZYHufW5Mf1VKAxgZ2eh1TpPHNLZvn37Exs5ePAgVqxYgR9++KFBYd9QCoUBzEwMZRf4REQN0egx/CNHjuDTTz/F999/Dzc3NzFqIiKiJtDowH/vvfdgaGiIf//733U/++GHH9C6devGNk1ERCJqdOCfOXNGjDqIiKiJ8YIxREQywcAnIpIJBj4RkUzozZXDFIqGX3O+Mes2V+yzPLDPLV9D+9uQ9UQ905aIiPQXh3SIiGSCgU9EJBMMfCIimWDgExHJBAOfiEgmGPhERDLBwCcikgkGPhGRTDDwiYhkolkEfkREBEJCQuDp6YmUlJSHLqNWq7FkyRIMGjQIgwcPxrZt23RcpXjq0981a9Zg5MiRGDVqFMaOHYvjx4/ruEpx1afP99y4cQPdu3dHRESEjqprGvXt8549exAWFobQ0FCEhYWhoKBAh1WKqz59LiwsxIwZMxAWFoZhw4Zh8eLFUKlUOq5UPMXFxZg+fTqGDh2KsLAwzJ49G0VFRQ8sV1lZiblz52Lw4MEYNmwYjhw5In4xQjMQExMjZGVlCcHBwUJycvJDl9m+fbvwyiuvCGq1WigsLBT69esnZGRk6LhScdSnv8eOHRMqKioEQRCExMREISAgQKisrNRlmaKqT58FQRBUKpUwadIk4a233hL+85//6LBC8dWnz3FxccLw4cOFvLw8QRAEoaSkRKiqqtJlmaKqT5+XLVtW997W1NQIzz33nBAdHa3LMkVVXFwsnDlzpu7xf/7zH+G99957YLmvvvpKWLhwoSAIgpCamir06dNHKCsrE7WWZrGHHxgYCGdn58cus2fPHowbNw4KhQK2trYYNGgQ9u3bp6MKxVWf/vbr1w+mpqYAAE9PTwiCgNu3b+uivCZRnz4DwIYNGzBgwAB4eHg0fVFNrD59/uGHH/DKK6/A3t4eAGBpaQljY2NdlNck6tNnAwMDlJeXQ6PRoKamBrW1tXB0dNRRheKzsbFB79696x77+/sjKyvrgeX27t2LiRMnAgA8PDzQrVs3HDt2TNRamkXg10d2djZcXFzqHjs7OyMnJ0fCinRnx44daNu2LZycnKQupUklJSXhxIkTePnll6UuRWeuX7+OjIwMvPjiixgzZgzWrl0LoYVf73DmzJlITU1F3759674CAgKkLksUGo0GW7ZsQUhIyAO/y8rKgqura93jpsiwFhP4cnXu3DmsWrUKX3zxhdSlNKna2lp8+OGHWLJkCZRKpdTl6IxarUZycjK+//57bNy4EceOHUNUVJTUZTWpffv2wdPTEydOnMCxY8dw/vz5Zvvf+j99/PHHMDMzw6RJkyR5/hYT+M7Ozvf9m5Sdnd3i93hjY2Px7rvvYs2aNejQoYPU5TSp/Px83Lx5EzNmzEBISAh+/PFH/Pe//8WHH34odWlNysXFBcOGDYORkREsLCwwcOBAxMXFSV1Wk9q0aRNGjRoFhUIBS0tLhISE4OzZs1KX1WgRERFIT0/HypUroVA8GL0uLi7IzMyse9wUGdZiAn/YsGHYtm0bNBoNioqKcOjQIQwdOlTqsppMXFwc5s2bhy+//BJdu3aVupwm5+LigrNnz+Lw4cM4fPgwpkyZgvHjx+Pjjz+WurQmFRoaihMnTkAQBNTW1uLMmTPw8vKSuqwm5ebmVjd2XVNTg9OnT6Nz584SV9U4K1aswOXLl7FmzRoYGRk9dJlhw4Zh69atAIC0tDTEx8ejX79+otbRLAJ/2bJl6N+/P3JycjB16lSMHDkSADB9+nTEx8cDAEaPHg03NzcMGTIE48ePx6xZs+Du7i5l2Q1Wn/4uWbIEVVVV+OijjzB69GiMHj0aycnJUpbdKPXpc0tTnz6PHDkSdnZ2GDFiBMLDw9GpUyc899xzUpbdKPXp88KFC3HhwgWEhYUhPDwcHh4eGD9+vJRlN8rVq1cRGRmJvLw8TJw4EaNHj8asWbMA3M2t3NxcAMC0adNQUlKCwYMH47XXXsPSpUthYWEhai284xURkUw0iz18IiJqPAY+EZFMMPCJiGSCgU9EJBMMfCIimWDgExHJBAOfiEgmGPhERDLx/wNR2ibGEZM3AwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "\n", "def funkcija(x):\n", " return x**4-2*x**3+5*x**2-6\n", "\n", "def fprime(x):\n", " return 4*x**3-6*x**2+10*x\n", " \n", "n=input('Unijeti broj iteracija: ')\n", "x0=input('Unijeti pocetnu vrijednost x: ')\n", "\n", "x0=float(x0)\n", "n=int(n)\n", "\n", "for k in range (0,n):\n", " x=x0-funkcija(x0)/fprime(x0)\n", " x0=x\n", "\n", "print(x)\n", "\n", "x=np.linspace(1,2,100)\n", "\n", "import seaborn as sns\n", "sns.set()\n", "plt.plot(x,funkcija(x))\n", "plt.title('Funkcija $x^4-2x^3+5x^2-6=0$ na intervalu [1,2]')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Zadatak 3 \n", "Napisati program koji učitava matricu kolonu X, a zatim izračunava standardnu devijaciju učitane matrice. Ispisati vrijednost standardne devijacije pozivom na funkciju. Standardna devijacija se računa prema formuli $\\sigma=\\sqrt{\\frac{1}{N-1}\\Sigma_i\\left(x_i-\\overline{x}\\right)^2}$, gdje je N  broj elemenata matrice, a $\\overline{x}$ srednja vrijednost svih elemenata." ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Broj clanova niza 5\n", "Unesi element niza 1\n", "Unesi element niza 2\n", "Unesi element niza 3\n", "Unesi element niza 4\n", "Unesi element niza 5\n", "Standardna devijacija niza elemenata [1. 2. 3. 4. 5.] je jednaka 1.5811388300841898\n" ] } ], "source": [ "import numpy as np\n", "\n", "N=int(input('Broj clanova niza '))\n", "A=np.zeros(N)\n", "Sum=0.\n", "\n", "for k in range(0,N):\n", " A[k]=input('Unesi element niza ')\n", " Sum=Sum+A[k] #Suma elemenata niza\n", "\n", "xbar=Sum/N #Srednja vrijednost\n", "\n", "Sum_i=0. #Suma ispod korijena\n", "for j in range(0,N):\n", " Sum_i=Sum_i+(A[j]-xbar)**2\n", "\n", "std_dev=np.sqrt(1./(N-1)*Sum_i) #Standardna devijacija\n", "\n", "print('Standardna devijacija niza elemenata ',A, ' je jednaka ',std_dev)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Zadatak 4\n", "Korištenjem operatora dvotačka definisati vremenski interval od 0 do 15 sekundi sa korakom 0.05 sekundi, a zatim u tom intervalu na jednom grafu nacrtati funkcije  $x\\left(t\\right)=\\cos\\left(2\\pi ft\\right)$  i $v\\left(t\\right)=-2\\pi f\\sin\\ \\left(2\\pi ft\\right)$. Uzeti da je f=0.2 Hz. Obavezno napisati naziv grafa (\"Zavisnost elongacije i brzine od vremena\") i označiti na grafu pripadnost linija odgovarajućim funkcijama, te nacrtati gridlines." ] }, { "cell_type": "code", "execution_count": 88, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEcCAYAAAD6GqKbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsnXl4VdXV8H/n3CnzSOYBEhJmkECCIDKGQRFRURAtWIcXxXmoX2urVWnVvtiW+qmtn7VWa7XOigo4FEVQUBREGUMCZCDzPCd3OGd/f9zcC4Ek3CT33CTk/p4nz5PknnP2Onefs9fea629liSEEHjx4sWLFy9uRu5rAbx48eLFy7mJV8F48eLFixdN8CoYL168ePGiCV4F48WLFy9eNMGrYLx48eLFiyZ4FYwXL168eNEEr4IZpBQXF5OWloaiKH0tSrcpLCxk5MiR2Gy2vhalQ3bv3s3ChQudf19yySXs2rWr19c5nZEjR5Kfn98jGU+lvzwL7733Htdcc02fyuDFvej7WgAvnfPhhx/yyCOPnPH/5uZm7rzzTu64444eXzs2Npa9e/f2RrweU1hYSGZmJgcPHkSvP/cewfT0dD799FPn35s2bXLLdbSiL58FL+c2597bfQ6xZMkSlixZ0u5/b7/9Nk899RTLly/vI6m8DERsNts5qcw7YjDda3/HayIbQBw6dIgnnniC9evXExkZCcC7777LxRdfTFpaGpmZmbzxxhvO4y+++GK2bt3q/Ntms3H++edz8ODBM8xM7733HpmZmaSlpTF37lw+/PBD5/+vueYa1q1bR0ZGBnPnzmXbtm3Oa5aVlbFmzRqmTJnC/Pnzeeutt5yf7du3j6VLlzJp0iQuuOAC/vCHPwCwcuVKADIyMkhLS+tw9qyqKn//+9+ZN28e559/PnfffTe1tbUdfi9dyfDMM89w991388tf/pK0tDQuueQS9u/f7/z84MGDXH755aSlpXHXXXdxzz338Je//AWAuro6brnlFqZOnUpGRga33HILpaWlznNra2v59a9/zYUXXkhGRga33XYbALt27WLmzJnO4+bOncvOnTu7fV+nX6cjtm3bRmZmJueffz7r1q1DVVXA3m8rVqzgiSeeYMqUKTzzzDMsWbKEtLQ058/IkSPZtWvXGc/CqlWreOqpp1ixYgVpaWnceOONVFdXO9v88ccfWbFiBenp6SxZsqRL89+xY8dYtWoV6enpXHLJJXz++efOz2pqalizZg2TJk3iqquuoqCgoNPr3HTTTbz66qvt/rdkyRI+++wzwG4ufO2111iwYAELFixwtn3DDTcwZcoUFi5cyObNm53nPvDAAzz66KP8z//8D2lpaaxYsYKKigoef/xxMjIyuOiiizh06JDz+LKyMu68806mTp3K3LlzeeWVV5yfne0Zc/R3WloaixYt4r///W+n93nOIbwMCOrq6kRmZqZ4/vnn2/1/69atIj8/X6iqKnbt2iUmTJggDhw4IIQQ4plnnhH33Xdfu2MXLlwohBDixIkTYsSIEcJqtYqmpiaRlpYmjh07JoQQoqysTGRnZwshhHj33XfFmDFjxJtvvilsNpt47bXXxPTp04WqqkIIIX72s5+JRx55RLS2topDhw6J888/X+zcuVMIIcTy5cvF+++/L4QQorGxUezdu/eMtjvjpZdeEsuWLRMlJSXCbDaL3/72t+Lee+/t8PyuZHj66afFuHHjxJdffilsNpv405/+JJYtWyaEEMJsNovZs2eLl19+WVgsFvHpp5+KsWPHivXr1wshhKiurhaffPKJaG5uFg0NDeLOO+8Ut956q1PG1atXi7vvvlvU1tYKi8Uidu3aJYQQ4ttvvxUzZsxwHjdnzhyxY8eOs97X6Zx+ndMZMWKEWLlypaipqRFFRUViwYIF4q233nL22+jRo8Urr7wirFaraGlpaXfuG2+8IRYuXCgaGhrO+D5XrlwpMjMzxfHjx0VLS4tYuXKl+OMf/yiEEKK0tFRMmTJFfPnll0JRFPH111+LKVOmiKqqqjPks1gsYt68eeK5554TZrNZ7Ny5U0ycONH5nN1zzz3irrvuEk1NTeLIkSPiwgsvFCtWrOjwXt9//31x9dVXO//OyckRkydPFmaz2fldXH/99aKmpka0tLSIpqYmMXPmTPHOO+8Iq9UqDhw4IKZMmeJ8rn/1q1+JKVOmiP3794vW1laxatUqMWfOHPH+++8Lm80m1q9fL1auXCmEEEJRFHHFFVeIZ555RpjNZlFQUCDmzp0rtm/fLoTo+hkTQojNmzeL0tJSoSiK2LRpkzjvvPNEWVlZp/16LuFdwQwAhBD86le/IjU1ldWrV7f7bPbs2SQmJiJJElOmTGH69Ons3r0bgEsvvZQvvviClpYWAD766CMWL17cYRuyLJOTk0NrayuRkZGkpqY6P4uNjWX58uXodDquuOIKKioqqKyspKSkhD179nD//fdjMpkYPXo0y5Yt44MPPgBAr9dTUFBAdXU1/v7+TJw40eV7fvPNN7n33nuJjo7GaDRyxx138Omnn57h2D+bDACTJ09m1qxZ6HQ6LrvsMrKysgD46aefsNlsXHfddRgMBhYsWMD48eOd54WGhrJw4UJ8fX0JCAjg1ltv5fvvvwegvLyc7du3s3btWoKDgzEYDEyZMsVt9+Uqq1evJiQkhNjYWK677jo2btzo/CwyMpJVq1ah1+vx8fFx/n/37t089dRTPPfccwQEBHR43aVLl5KUlISPjw8XXXQRhw8fBuCDDz5g5syZzJo1C1mWmT59OuPGjWu3qnXw008/0dzczM0334zRaGTatGnMmTOHTZs2oSgKn332GXfddRd+fn6MGDGCK664otP7nDdvHllZWRQVFQH2Z3n+/PkYjUbnMTfffDMhISH4+Pjw5ZdfEhcXx5VXXoler2fs2LEsXLiwnU9r/vz5jBs3DpPJxPz58zGZTFx++eXodDoWLVrkvOf9+/dTXV3NHXfcgdFoJCEhgeXLl7dbEXX2jIHdkhAVFYUsyyxatIihQ4eyb9++Tu/1XMJrqBwAvPDCC+Tk5PDee+8hSVK7z7Zt28Zf//pX8vLyUFWV1tZWRowYAcDQoUMZPnw4W7duZc6cOXzxxRds2LDhjOv7+fnxl7/8hX/+8588+OCDTJo0iV/96lcMHz4cgCFDhjiP9fX1BeyBBrW1tQQHB7cbpGJjYzlw4AAAjz/+OE8//TQXX3wx8fHx3HHHHcyZM8eley4uLub2229Hlk/OgWRZpqqqqt1x5eXlXcpwuvw+Pj6YzWZsNhvl5eVERUW1+05jYmKcv7e0tPCHP/yBr776irq6OgCamppQFIXS0lKCg4MJDg526X5cua+oqKhuXet0eePi4igvL3f+HR0dfcbxJSUl3HPPPfzv//4vSUlJnV43IiLC+buvry/Nzc1O+T/55JMOTa+nU15eTnR0dLt7jY2NpaysjOrqamw2Wzv5Y2NjO5UnICCAWbNmsWnTJm6++WY2bdrE73//+3bHnHqtoqIi9u3bR3p6uvN/iqK082mGh4c7f/fx8TnjOXHcc1FREeXl5Wdc69S/O3vG9Ho9GzZs4KWXXnIqx+bmZmpqajq913MJr4Lp5+zatYvnnnuO1157jaCgoHafWSwW7rrrLtatW0dmZiYGg4HbbrsNcUqC7MWLF7Nx40ZUVSUlJYWhQ4d22M6MGTOYMWMGra2tPPXUU/z2t7/lP//5T5eyRUZGUldXR2Njo3OALykpcQ6Uw4YNY/369aiq6pyt7tq16wwl2RHR0dE88cQTTJ48+YzPCgsLXZahKyIiIigrK0MI4ZSppKSEhIQEAP75z3+Sm5vLW2+9RUREBIcPH+byyy9HCEF0dDR1dXXU19ef0S89va+eUFJS4lxtFhcXO31zwBnfc2trK7fffjs///nPmTVrVo/ai4mJ4bLLLuOxxx4767GRkZGUlpaiqqpTyZSUlDBs2DDCwsLQ6/WUlJQ4JzIlJSVdXm/x4sU8++yzZGRk0NraeoZSO32ikJGRwUsvvdTdWzyDmJgY4uPjnf6e7lBUVMRDDz3Eyy+/TFpamnOFM1jwmsj6MeXl5dx333385je/YcyYMWd8brFYsFgszpd127Zt7Nixo90xixYtYseOHbz++uudmscqKyv5/PPPaW5uxmg04ufnh06nO6t8MTExpKWlsX79esxmM1lZWbzzzjtceumlgN2cUl1djSzLzkFYp9MRFhaGLMucOHGi02tfc801PPXUU85ZX3V1NVu2bOm2DF0xceJEdDodr776KjabjS1btrRzzjY1NWEymQgKCqK2tpZnn33W+VlkZCQzZ85k7dq11NXVYbVaneazrnD1vlzlxRdfpK6ujpKSEl555RUWLVrU6bG/+c1vSEpKOsPM2h2WLFnC1q1b+eqrr1AUBbPZzK5du9oFPziYMGECvr6+/OMf/8BqtbJr1y6++OILFi1ahE6nY/78+Tz77LO0tLRw9OhR3n///S7bnjVrFsXFxTz99NMsWrSo3crodGbPnk1eXh4bNmzAarVitVrZt28fx44d6/Y9T5gwgYCAAP7+97/T2tqKoihkZ2e7ZOZqaWlBkiTCwsIAe1BOTk5Ot2UYqHgVTD/mrbfeorKykieeeKJd9E9aWhoPP/wwAQEBPPTQQ9xzzz1kZGSwceNG5s6d2+4akZGRTJw4kb1793Y6+KiqyksvvcSMGTOYMmUK33//fYf7bzpi/fr1FBUVMWPGDO644w7uvPNOpk+fDsBXX33FJZdcQlpaGo8//jh/+ctfMJlM+Pr6smbNGq655hrS09P58ccfz7juddddx9y5c7nxxhtJS0tj+fLlnb7QXcnQFUajkWeeeYZ33nmHjIwMPvzwQ2bPnu206//85z/HbDYzdepUrr76ambMmNHu/CeffBK9Xs/FF1/MBRdcwL/+9a+zttmd+3KFzMxMli5dyuWXX87s2bO56qqrOj1206ZNbNmypd1z5PDXuUpMTAx/+9vfeP7555k2bRqzZs3ixRdfdEavnYrRaOS5555j+/btTJ06lbVr1/Lkk086VywPP/wwzc3NTJ8+nQceeIClS5d22bbRaGT+/Pns3Lmz08mSg4CAAF588UU2b97MjBkzuPDCC/nTn/6ExWLp1v2CfVL03HPPkZWVRWZmJlOnTuWhhx6isbHxrOempKRw4403smLFCi644AKys7OZNGlSt2UYqEhCeAuOefHiYNmyZaxYsYIrr7zSbdecPXs2f/zjH8nIyHDbNb14GQh4VzBeBjXfffcdFRUV2Gw23n//fY4cOXLGSqU3VFdXU11dTVxcnNuu6cXLQMHr5PcyqMnNzeWee+6hubmZhIQEnn766XaO8t6wb98+brzxRlauXNllhJQXL+cqXhOZFy9evHjRBK+JzIsXL168aIJXwXjx4sWLF03wKhgvXrx48aIJg9LJX1PThKr2zPUUHh5AVdXZ49/7iv4uH3hldAf9XT7wyugu+oOMsiwRGurf7fMGpYJRVdFjBeM4vz/T3+UDr4zuoL/LB14Z3cVAkLEjvCYyL168ePGiCV4F48WLFy9eNMGrYLx48eLFiyZ4FYwXL168eNEEr4Lx4sWLFy+a4FUwXrx48eJFEwZlmLIXbVFqi7F89w62osNIegP6pAyM6Zcj+wT2tWiDGlvRISx7NqBU5CL5hWAYPQvj+IVIOkNfizZoEULFlr0Dy75PUetLkUNiMJ63CP3wqS5Vfu3veFcwXtyKregQze89iq34MIaUqehiR2M9vJXm93+H2lDR1+INWiwH/kvLpidRm6oxjJ6NHBSJ5bt3aPl4PcLa2tfiDUqEqtK67SVat70IegOG0fZiga1fPI9552sIcWYRt4GGdwXjxW0oNcW0fPp/kQMj8V30C2T/UPv/x82n+eP1tGz+M35LHwW8KxlPYs3djXnna+iHTcJn7hokvb1ipzV7B63b/kHr1hfwmX/HOTFjHkhYdr+HLfsrjJOWYJx8BZIkIVQV8643se7/FMkvBFNa15U7+zveFYwXtyBUG62fP4ekN7ZTLgC6qBR859+BWl+Geed/+lDKwYfaVEPrtheRI5LwybzNqVwADCOmY5qyHFveHqyHt/ahlIMPW/FhLD9uxDBqJqb0pU7lLskypqkr0CdPwbL7XZTyY30sae/wKhgvbsF68HPU6hOYZl7fTrk40MeNwTjhYqxHttNalN0HEg5OzLveBMWKb+atSLozDRaGCRehix2N+ft3UVsb+kDCwYdQbZh3/BspMALTBT8743NJkvCZeQOSbzCtO14d0KYyr4Lx0muEuQnzng3oEsajHzqp0+OMk5Yg+YVQteVlvHXutEcpP47t6LcYJ1yMHNRxlU5JkjBdsBIsLVh++NDDEg5OrFnbUWuK8Zl2LZLe1OExktEX0/nLUStyaTz4tYcldB9eBeOl11gOfg6WFkwZV3Vpx5cMPhgnLcFceASlJMuDEg5OLHs/ApM/xvMWdXmcLiwOfep0rIe3oTbXeUi6wYlQbVh+3IQclYJu6MQuj9WnTEUOi6d2x7sDdhXjVTBeeoWwWbAe+C+6xPPQDRl61uMNIy5EFxCKZe9GD0g3eFFqirDl78U4bj6S0fesx5smXgKqFeuB/3pAusGL7eguRGMVprTFZw2qkCQZ48TFWCsLseXt9ZCE7sWrYLz0Ctvx7xGtDRgnXOTS8ZLeSNDki1CKDqLWlmos3eDFeugL0OkxjM106Xg5JBp94kSsWdsQilVj6QYvlkOfIwdHo0s4z6Xj9ckZ6IKG2PtzAOJVMF56heXQ58ghMehiRrl8TuDETJB1WAboS9PfEdZWrNk70CdldGtzq2FsJqK1AVvubg2lG7wolXmo5ccxjM10OSRcknUEpc1vm5CVaCyh+/EqGC89Rqkusr8wo2d3aw+FPiAU/bBJ2HJ2IlSbhhIOTmy5u8HaimHMnG6dp4sbgxQYgfXIwHUq92esR74CnQFD6gXdOi9wYiZIMtbsHRpJph1eBeOlx9iOfgOSjH741G6fa0idjjA3opw4oIFkgxtrzjdIgRHoolK7dZ4kyRhSL0ApOoTaVKORdIMTodqwHfsO/dA0JFP3Sg/rA0LRxY/DevSbAefs9yoYLz1CCBXr0W/QxY9D9gvu9vm6hHFIpgCsOTs1kG7wojbVoBQdwpB6QY925htSpgEC29Fv3S/cIEY5cQDR2tDt1YsDQ+o0RGMVSsnA2kPmVTBeeoRamY9orMIw/PwenS/JevTJ6dgKfkLYLG6WbvBiy98LCPQ97Bc5JBp5yFCseXvcK9ggx5q7B4x+6BLG9eh8/dBJoDNgG2D94lUwXnqELX8vSBL6RNeiYTpCP3QS2MzePTFuxJa/Fyk4CjkkpsfX0A+dhFp2zLsnxk0IVUUp+BF94gQkuWfpHyWDCV3cGGz5ewfUJmWvgvHSI2x5e9FFj0DyCejxNXSxo0BvGrAx/v0NYWlBKTpst/P3InGlfuhEQKAU/OQ+4QYxSvkxRGsD+qFpvbqOfmgaoqEStbrQTZJpj1fBeOk2akMFavWJXr8wkt6IPmH8gJuV9VdshQdAtfW6X+TwRKSA8DZzm5feouTvBUmHPmF8r66jb9v5P5D6pV8pmHXr1jF37lxGjhxJdnbHzixFUVi7di3z5s1j/vz5vP322x6W0ost/0fg5APfG/RD0xDNtaiV+b2+1mDHlr8XTP7oolJ6dR1JktAPnYit8CDCZnaTdIMXW/5edLGjkIx+vbqO7BeCHJnsVTA9JTMzk9dee424uLhOj/noo48oKCjgs88+48033+SZZ56hsHDgLBnPBWz5e5FDYpGDo3t9LV3iBJAkbPk/uEGywYtQFWwFP6FPnIgk63p9Pf3QNFAsKEWH3CDd4EWtLUWtLen1qtKBfmgaakXugAkj71cKJj09nZiYrp2TmzdvZtmyZciyTFhYGPPmzeOTTz7xkIRehKUFpfiIW1YvALJPILroEc5VkZeeoZQdBXOT2/pFFzMKDL7efukltgL3rfbt15nUdt2B4R8bcBUtS0pKiI2Ndf4dExNDaan2Oa1aLTb++v4BggJMTE4dQlrqkEFZAVApPQJCQddLe/Kp6OLGYtn9HqK10eWgAZuisutQGd9nldPQbCE23J/Zk+IYHtv9PTnnAkrRIXtUX9wYt1xP0unRx47CVny4W+e1mG18ubeIA7nVCCAhIoDM9HgiQ86ecPNcxFZ0CDk4GjlwiFuuJ4fGIvmF2Pt79GyXz6uub+XLH4sZMzSUUUPPrNekFQNOwbiD8PDuRz61WmxEhfuzL6eCb/aXcP7YaO67dhJ+PgYNJOwdERHalSSu+ukYrToDUWMnIp9SHbG7nCpj69h0ine/h19THgEJ084uQ10Lf3llN4fzqokZ4k9UqB8/Hqtix4FSls5O4eeXjEGWe6/8tfwe3cGp8hVXZGOKTiYyvvdmSwd1I9Oo+mwvIcYWDMEd15M5laMnavnDv76nvKaF4fHB+Jr0bN1bxLYfi7jtqvPIzEh0m2zuRKt+FoqNvLJsAsfNYkgv2zhVxvLkCTQf/5EhQwJcmuR+/VMRT7+5F4tVZXhCqEef6wGnYGJiYiguLmbChAnAmSsaV6iqakRVux+1tHJeKmHLJ/LGp1m8vfUoDz23g/uunojJ0Hubt7uIiAikokK7yoRNR39CjkqhqsYM9MwBfLqMwhAJBh9qDv9AS3jXG9Hqmy2se+0HquvNrL50DFPHRCFJEq0WG299cZT3vjxKdV0LqxaM6NUKU+vvsbecKp+wmmktysY4fqFbZVaCkgGo2L8bw8gZXR5bUNbAuv/sxc+k4zcrJ5MSH0xERCDZxyv5x8ZDPPXGXmrrWph5XvfeVa3Rsp+VsqMISyuWsJRetXG6jLbwVNQD2ynLzkIXFt/lubsOlfH3Dw+SHBfELZeOZUiIb49kkWWpRxPzfuWDcYWLLrqIt99+G1VVqa6uZsuWLSxcuNBj7et0MgsyElh96RiOFtbx6mdHPNZ2XyNaG1GrCtDFjnbrdSVZjy56BMpZzDGqKnj+g4NU1rVyz7IJTBsb7VQiPkY9qxaO5OLzE/lybxFb9xa5Vcb+jFKWA6ri9n6RQ+OQfALPaiZrbLHy9Lv78DXpeOBnduXiIDTQxL3Lz2NcchivfHKEo4WDZ/OmrS1AQhfreqZxV3Bc72zvS0FZAy9uOkRqfDD3r0hjSB+YKfuVgnnssceYOXMmpaWl3HDDDVxyySUArF69mv379wNw2WWXER8fz4IFC1i+fDm33347CQkJHpd1yugoFl8wjB37S/nucJnH2+8LbG077vVuHsgc11RrS1Cbazs95tPvCzicX8PK+SMYmXimHVmSJK6cPZzxyeG88XkOxZVNbpezP6IUHwZJhy56hFuvK0kSutjRKMVZXe5T+venR6hrtHD7FeMJD/Y543O9TmbNknGEB5t4/sMDtFoGRwZtpSQLOTyhWyUTXEEOjEAKjOhSwVhtKv/vg4ME+Bq4fen4PrOy9CsF89BDD7F9+3YOHTrEjh072LRpEwAvvPAC48fbnco6nY61a9eyZcsWtmzZwtVXX91n8i65cBhDowJ5/fMcWszn/kujFB0GvQk5Msnt19bF2ZVWZy9NZV0LH3ydS1rqEC6c0HmkoSxJ3HTJaIx6Ha9+dmRQbOC0FR1GFzUcydBxfffeoIsdjWiqRtR3PInaf7yK77PKWXJhEkkxQZ1ex89Hz+rFY6mqN/Phjjy3y9nfEDYLSmkOulj3BF2cjj52NLbiLITacXblT3blU1rdzI2LRhPo13NfaW/pVwpmoKGTZVYtHEldo4XN3577GwWV4sPoYkb2OJ9SV8hhiWDytyuxDnhv23EQcM281LP6VoL8jVw5K5msglp+yK50u6z9CWFpRq3Mdbt5zIG+TfHbOugXVRW88XkOUWF+XDTl7A78lPhgLpwQw3+/P0FFbYvbZe1PKOXHQLFqstqHtgmZpRm1quCMz2obzWz6Jp/0kRGMSw7XpH1X8SqYXpIcG8SU0ZFs2V1IfdO5mxVYba5DrS1G72Z7sgNJltHHjOzQ3l9Y0ciuQ2VkpsczJNg1O/LMibFEhfmx4evjPQroGCgoJdkghNvt/A6koCgk/9AOV5bfHCylpKqZK2cmY9C7NpRcMSMZWZb48Otcd4var1CKs0CS0MW412zpwDGh6KhfNu3MR1EFV84erknb3cGrYNzAZRcmYbEpfPr9mbOJcwWl7CiA2+38p6KLTkU0VJyRxXfjzjxMRh0Xnz/U9WvJMlfMSKKoookfsivcLWq/QSk7ave/RCZrcn1JktBFpTr734GqCjbuzCMxKoDJIyNcvl5ooIm5k+LYebCU8nN4FaOUHUUOS+h1epjOkP1C7H6Y0/qlttHMtp+KmD4+hqhQbdruDl4F4wZiwv2ZPDKSL/cWn7O+GKXsKMh65CGuD/LdxVGBUSk/+dJU1rbwfVY5syfGEeDbvT1H6SMjiQzx5ZPvCs5ZX4xSdhR5SCKS3v3+Fwe66FREUzVqY7Xzfz8eraSspoVFU4d2Oxx8QUYisiTx3+9PuFvUfoFQVZTyY73OCXc2dFEp9lDoU57tz/cUoiiCi6f2jz1HXgXjJi6akkiL2caO/SV9LYomqGVHkSOGIem021gqDxkKsh6l9KSC2bKnEFmSmJfedbx/h9eTJRZMSeB4cT3Hi+vdKWq/QKgKSsVx7QeySLup5VTF/9n3JxgS7NOt1YuD0EATU8dG8dW+Yppbz70JmVpTBNZWjygY0VKHaLD7GS1WhS/3FjFpRES/WL1ANxRMU1MTpaWlNDUNjtDP7pIcG8Sw6EC2/1Ryzs2WhWJDqczV/IWRdAbkiGGobct+m6Ky80ApE1OHEBZ0ZvirK0wbG43JoOOrfcXuFLVfoFafAJtF836RhySCzuhU/KXVzWSfqGV2Whw6uWdz1LmT4rFY1XMyxN9pTvaAgoGTiv+HnAqaWm3MndR5smBP0+XTkZ2dze9//3syMzNJT09n9uzZpKenM2/ePH73u99x5Mjg2WToCjMmxFBY0Uh+Wf/dAd4T1Kp8UGzOmayW6CKHo1TmIhQbP+ZU0thiZUYXYclnw9ekJ2NUJLsOl5+THV4iAAAgAElEQVRz+y8cA77mil/Wo4sY5hzIvt5XgixJXDCu52lphkUHEhfhz1f7zr0Vv1J2FMk3CCmw+6u77iCHxYPe5HwOvvqphCHBPoz0YK6xs9Gpgrnvvvv4xS9+QUREBH/84x/59ttvOXDgAN9++y1PPvkkUVFR3H///dx7772elLdfc/6YKAx6+Zx7aTw1kDnbUGyoVfl8vb+E0EAT45J6F2o547wYzBaF3VnnlrNfKT+K5B+K5B+meVu6qBTUynxsllZ2HChhwvBwQgJ67veRJIkZE2LJLamnsKLRjZL2PUr5UXRRKZonw5Vke3CHUn6UytoWDufXcOGEGOR+lIS3UwWzePFiPvroI9asWcOkSZMIDg5Gr9cTHBzMpEmTuOWWW/joo49YsmSJJ+Xt1/j5GJg8MoJvD5ZhsSp9LY7bUMpykAKHIPtrPzNyKLGG/Cz2H69i+vjoXieuTIkLJirM75wzkyllnhnIoC0AQ1U4um8/dY2WLje7usrUsVHoZImvfjp3JmRqSz2irswjkzFoU/xVJ/jmp3wkYPq43veLO+lUwcydO9elC8yZM8dtwpwLzBgfQ4vZds6Exgoh7ANZpGdeGNk/FCkgnOrjBxECLhzf+xfGPluOIaewjtLqZjdI2ffYGmoQDZWe65cou3m0NGs/QX4GJgzv/Qa+ID8jE1OH8M3BUmxKxzvSBxpq2TEAZA8qGIRKwaH9jEkK6zBVT1/ikodu48aNHDtm/+KOHz/Oz372M6677jrn/7ycZOTQUEIDTXyfVd7XorgF0VSNaK5FF+W5TVu6qBR86vJJiQsm0k3RMNPG2v0F50q/mIvsJcU91S+ybxAERmCsy2PK6Cj0OvcEoF4wLprGFitZ+QOjQuPZUMqP2fclDRnmkfYcftEwczEXjHVfqQZ34dJT8tRTTxEcbM+Q+uSTTzJhwgQyMjJYu3atpsINRGRJYvKICA7kVp8TTmWlwr7jWquNfB3RHBBPIE1MS3bfbCw00MTwuCD2HDlHFEyJfYOlHO65/Q51PnHEy5Wkjzp7bRhXGZcUhsmoY/eRc2PFr1TkIofHI/WiVlJ3kHwCaDKEkqiv4rwU9xQ1cycuKZjq6mqGDBmC2Wxmz5493Hvvvdx+++1kZWVpLd+AZPLICKw2lX3HqvpalF6jVuTZB7Iwz2Wszmq0T2bOC3VvSPzkEZEUlDWeEzvIzaXHkMM8N5AB5LSEEKprZrgbXXEGvY7zhoezN6diwKf0EUKgVOahG+L+ZLBdtZlrDiPZpwY/n/5X3sslBRMWFkZ+fj7bt29n/PjxGI1GzGbzObffw12kxocQ5GdgzzkwK1Mq85DD4jw6kH11Qo8A/Fvc65R3bAr8YYD3ixACc8kxdBHDPNamxarwfbk9D5yodm9i10kjImhotpJT2HmphoGAaKgAcxOyB/uloKyRo60hBKgNqC39bzOxSwrmtttuY+nSpTz44IPcdNNNAHzzzTeMGqVNgr2BjixLpI2IYN+xKqy2gRtNJoRAqcj16EBWXd9KdmkrzaYIp3nOXUSE+JIYFcCe7IFtJhMNlagtjcgRnpspH8ytJrc1BIGEUpHn1mtPGB6OQS8PeDOZ43vRebBf9mSXU6jYTWOqm/vFHXSpYFpa7KaEpUuX8vXXX7Nt2zamT58OwHnnncf69eu1l3CAMnlEBGarwoHc6rMf3E8RDZX2GZkHl/x72qLvjFHJqBV5bl8lTx4ZybGiemoaelbuuT+gVLb5xTyo+HcfqUDv44scEu12xe9j1DMuKYwfsitQB7BVRKnIBZ0eOdRzO+n3HKnANyYJkJzPRX+iSwUzZ84cbrzxRl5++WXKy8vx9T2ZKj08PJyICG13qg5kRg0NxdekZ2/OwK1H0hcD2Y85lcSE+xEQ35ZnqYsKlz1h8gj7M/vj0YHbL2pFXttA1v38bD1qTxXsO1bJeSlD0EUkoVbmub2NSSMiqGkwUzCAs2ColXn2DMo6z/hCyqqbKalqZvzIeOSQ6IG3gvnqq69YvXo1paWlrFmzhgULFvDYY4/x9ddfY7Gcu7VP3IFeJzNmWCgHc6sHrK9KrcgDWWdPSeEBWi02cgprmTA83Bnm6e6XJibcj/AgHw4cH7gBGEpFLqbIoR4byHJL6mlqtdn7JSIJ0VyL2uTesGJHYawDxwfmil8IFaUiz6Pmsf1tz/D44eHIQ4a5fWXpDrpUMAaDgWnTpvHAAw/w8ccf8+KLLzJs2DD+9a9/MXPmTNasWcPrr79OZeXAnQ1qyfjkcGoazBRVDMwEoUpFbtuMTLsMyqeSVVCLTRGMTw63J1iU3L/slySJ8cPDOZRfMyA39wmholTmYYzx3L6k/cerkCQYMyxMM8Uf7G9kaHSgc9AcaIi6crC2eGz/C8D+49VEhfkRGeKrmeLvLd3aLZWQkMDKlSt54YUX2LZtGytWrCA7O5utW7dqJd+AZlySPUfU/tyB99I4BjJPmsf2H6/CZNCRGh+CpDchh8a53aEMMD4pDLNF4Whh3dkP7meI+nKwtGDyqIKpJjk2iABfg2aKH2B8chjHiuppbrW6/dpa4/g+PBV4YbEqHCmoYXzbGONoVwvzZW/oloJpbGykrKyMsrIyamtrGT16NI888gjLli3TSr4BTViQD/ER/uwfgPthRH0FWFqQPTQjE0Kw/1gVo4eGOsvvykOSUCty3W5iHDU0FJ0sDcjZskPhmqI9o2Dqmy3kldQzvi3hqJaKf1xSOKoQHMrrX7NwV1Aq8kBnQA6N9Uh72SdqsdhUp2lRF96m+PuZH8YlBbNz504yMzPJyMhg1qxZzp/Zs2drLN7AZ1xyODmFdQOu0qVSZS//rNOwguWplNe0UFnXyrjkk5mBdUMSEa0NiBb3rjR8TXpS44PZPwDt/WpVAcg6jBGe8Ysdyq1GYLfzO5DDE+1yuJnhcUH4mvQDUvGrVQX2ja+yziPt7T9ejUEvMyoxBADJYEIOjtakX3qDSwrmwQcfZM2aNezevZuDBw86fw4cOKC1fAOe8cnhKKogq2BgzcrUqgKQJI+FXDoGFceMDHBmD9DipRmfHE5hReOAC1dWqk8gh8Z6zC+2/3g1Ab4GhkYHOv+nC0+w2/tb3RvxpZNlxg4L5cAAC4wRQqBWnUAX7rlsFwdyqxiZEILRcFKhyWEJKNX9qwy1SwrGbDazdOlS/P390el07X68dE1qfDAmg27A7YdRqk4gB8d4bAf/gdxqokJ9iQw5GQrveGEVDRSMQ5EdHGD9olYWIId5Jv+YEIKDedWMTQprV2PE0b5a5f7BbFxbYExx5cAJjBHNtQhzo8f6paqulZKqZqeP14EcnohoqESY+89355KCuf766/nHP/4xoGYV/QW9TiY1IXjAZYtVqwo8lkhRUVWyT9Qyelj7F0Yy+SMFhGsykMVF+BPgaxhQK0u1pR7RUuexmXJxVTP1TRZGn1YhUQ53rCzdmzIGcLaVVTBw0sY4vgfZQ/3ieGZPf1+cE7LqQo/I4QouKZgFCxbw1ltvMXnyZDIzM9v9eDk7oxNDKalqprZxYJhjRGsjoqnaYy9MXmkDrRbFaU8+FV14oiYKRpYkRiWGcDi/ZsBMnBymQk8pfsekaNTpCsY3CMkvBEWDfokI8SU8yIfDA2hC5vgePKX4D+fXEOBrIC7Cv93/Hc9Ff/LDuLRT66677iI9PZ2LLroIH5/+VdBmIDDKOSurYeqY/lez4XQcdlxPvTDOgSzxzDS9cngCtoIfETaL2811o4eGsvtIBeW1LUS5qe6MlqjOgcxzCiY8yIeIDopYyRopfrD3y94ce9qY/lT+tzPUqhNIgRFIRu2fISHs/txRiSFnfDeSXwiST6Bm/dITXFIwhYWFbNiwAVl2T5GhwcbQqEB8TXqy8msHhIJxPKAemykX1BIX4U+Q/5kKRA5LACFQa4rcvkvaqfjzawaEglGqCpD8Q5F8AjRvS20byCamDumwJLMuPAFL0UGEYnN7RoFRQ0P4en8JheWNJEYFnv2EPkatKvDYZKyitoXqejOLpp45GZMkCTm8fzn6XdIYmZmZfPvtt1rLcs4iyxIjE0IGjL1fqTqB5BOI5BuseVs2RSWnsLbD1QucnK1r4eiPDvMj2N84YOz9atUJj9XlKSxvpKnV1mm/yGEJoCqote4tqQAnV7IDwUwmbGbU+jKP9YvjWe2qX9TqQoTaP7K4uzT1sFgs3HrrraSnpxMe3r4W95NPPqmJYOcaoxJD+PFoJdX1rYQF9W8zo1ptd/B3NHN1N8eL67FY1U5fGCkoAgw+mtiVJUli1NBQstr8MJ64354iFCtqbQnGoRM90p7DbHm6g9/BSXv/Cbeb7MKCfIgM9SUrv4aFUzxXsbMnqNVFIIRH/WLB/kZiwjtecevCE7EqVtS6MnQe2vTZFS4pmNTUVFJTU7WW5ZzGYY45nF/D9PExfSxN5wjVhlpThGHsPI+0l1VQgwSM7MDBDyBJMnJYvKb2/l2HyiipaiZ2iP/ZT+gj1JpiEIoHI5VqiQzx7XQyJAdHgc6AUlWAgelub3/00FC+O1yGoqro+rFp3rkh2QP9IoTgcH4NIxNDOp0MnYzwKxg4CuaOO+7QWg4AcnNzeeCBB6itrSUkJIR169YxbNiwdsc888wz/Oc//yEy0l4XfNKkSTzyyCMeka83xEcG4O+j50hBbb9WMGptKSg2jzmSjxTUkhAZQIBv5xsHdeGJWI9+o8kqw6HYjhTU9G8F44wg034gU4XgyIlaMkZ1Xo5DasuyrVXE0sjEELb9WExBWSNJMUGatOEO1KoCMPggBQ7RvK2ymhbqmixnRPWdihwSC7LOLlfKVM1lOhudTg2ysrJcuoCrx7nCI488wrXXXsunn37Ktddey8MPP9zhcZdffjkffPABH3zwwYBQLmAPi02NDyGnqH8nWPTkQGZTVI4X15Oa0PHqxYEclgCWFkSj+7N2R4b4Euxv7Pf9olSdAJ0ROUj7IJHiiiZazDZS47vuF11YAmrVCU3CvEe0td3fE5I6TISSpP0qK+eE3f8yoot+kXR65JDYfuPo7/RbWbt2LatXr2bjxo2UlZW1+6y8vJxNmzaxevVqfv/737tFkKqqKg4dOsTixYsBWLx4MYcOHaK6emDttO6KlPhgyqqbqW/uv7V0lKoTIOuRQ7RfZZ0ob8RsVUiN7zqYQMsd/ZIkkRIfPAAGMkeuKw8MZIX2gexs/SKHJyLMjW4vCgd2P0x4kMkpS39ECNWeusdDDv6cojr8ffREd+J/caBlCHl36dRE9vrrr7N161beeOMNHnzwQWRZxt/fn6YmexqCadOmsXLlSmbNmuUWQUpKSoiKinKmn9HpdERGRlJSUkJYWPsdq5s2beLrr78mIiKCO++8k7S0tG61FR7euzDPiIiehU6mj43hnS+PUdFgYfjQ8LOf0EN6Kh9ASWMxUkQCkVGdL8PdQUREIN8cLgfg/AlxDDklRczpqMGjyUPCt6Wc0F7cW2dMHBnFniMVyEY94cEn5ejN9+hOhBDk1xTiP2pqO5m0ku9EZTMhgSbGpEZ2aZJsGT6Skp0QaKvAL6Jjk2pvZBw7fAgHjlUyZEiApgEYPZXRWlNKo7WV4GGpBGn8rEREBJJb0sCY5HCiIrs2GdYOTaE6Zwdhfio6f+0jQbuiSx/MnDlzmDNnDlarlby8PBoaGggODmbo0KHo9Z6ppnc6K1asYM2aNRgMBnbs2MFtt93G5s2bCQ11fUCsqmpEVXu2rI+ICKSiomdJ/kJ9deh1EnsOlTI8Spu9DL2RD6C1JBddwoReXeNsOGTce6Sc8CATwmo7a3tScBQNJ45i00Cu2FC7I3vXvmIyRkW2k7E/oDZWobY2YvGPdsqkpXwHjlUyPDaIysrGLo8TOvskqfp4Nk3BI874vLcyJkb4s31vEYePVhDRxQSkN/RGRmuu3T3QbIzErPH7ciy/iqKKRqaNjTqrvDZTFADl2YfRx491iwyyLPVoYu7Sevvo0aOkpqYyadIkhg8frolyiYmJoaysDEWxx28rikJ5eTkxMe1NNRERERgMdofw9OnTiYmJIScnx+3yaIFBr2NYdFC/XfarzXWIlnqPRcQcLawl5Sx2fge68ARNTGQACZEBGA2y08bd3/DkxteaBjOVda2kxp195isZ/ZACh2jm6E9pk6G/mi+dGcfDtM84fqztO0hxoV+ckWTVfZ8yxiUFc8MNN7Bo0SL+9re/ceKENra98PBwRo8ezcaNGwHYuHEjo0ePPsM8dqo/6PDhwxQVFZGU5Lk62L0lJT6Y/NIGrLb+sRHqVNS2JHmesClX1bVS22g5q53fgRyWgGioQFhb3S6LXieTHBPUbx39juSFujDta8AcbfsOXFb8bRv7tCA+IgBfk67/TsiqC5GCopD0Js3byimqQ6+TSIo5uylO9glsyxXX90kvXVIwO3bs4Je//CXHjx/nsssu4+qrr+bf//43VVXuLQz06KOP8uqrr7Jw4UJeffVV1q5dC8Dq1avZv38/AOvXr2fx4sUsWbKEhx56iCeffJKIiM7DKfsbqXHB2BRBbkn/ML+cilpTBOCRGjA53ZiRAc5Zolrj/p3jYFf8J8oaabX0v8Jwak0Rkn+YR3Jd5ZyoxaiXSXTRhCuHxqHWlSEU95c5lmWJ5Njgfqv41ZoidB6ql3S0sI6h0YEY9K6VSJFD45zvc1/ikq1Lp9Mxe/ZsZs+eTWtrK59//jmvv/4669atc2vRseHDh/P222+f8f8XXnjB+fu6devc1l5fMLxtxp5TWMuIs4Tnehq1pgjJFIDkq/2+g5yiOnxNOuIjXBvIdKH22btaU4QuMtnt8qTGh6CKfI4X1zPmtDTofY1aU+SxUrw5RXUkxwah17kWrSaHxYNQUOtK0Wmw8k2ND+aDr3JpbrXi5+OZImuuIBQran05+uQMzduyWBXySuuZl+769yuHxWM9tBUhVI+EUHcqR3cONpvNbN26lc2bN3PgwAHS09O1kuucJcjPSHSYX7+0Kys1RchhcR5JmXK0sJbk2GBk2bW2pMAI+85xjcwxw2ODkeh/9n6hqqi1JR5ZVbZabJwoayTFRbMlnFztqtXazJZT44IRwNGiek2u31PU2lIQqmdW+ydqsSnCJb+YAzk0FhQLosH9e8e6g0srmG3btvHRRx/xxRdfkJKSwqJFi3j00UcHlGmqP5ESH8ze7P6VjlwIgVpTjMEDu38bW6wUVTSRPjLS5XMkWUYOidVs2e/noycuwr/fmWNEQwUoVo+YYnKL61GFICXO9ZW1HBINkqxZvyTHBiNLEkeLapkwXLvQ/u7iNCd7wMF/OM++F3B4NxS/7hTFLwe5/p65G5dWMOvWrSM5OZkNGzbw1ltvcf3113uVSy9IjQumqdVGaVVzX4viRDTXgqXZnmpCY7LyqhGcfSPf6chh2tqVU+JDOFZU1+MQdi1QnH4x7fslp7AOCUiJc91EKukMyMFRmvWLyagjISqAnBP9S/GrNUUgycjB2mdWOJxbTVSYH0F+rtdDcqyslJq+dfS7pGA2b97MbbfdRmJi/85sOlBwpEbpT9Exnp6RyZLdgdsd5NA4RFONZjXHU+ODabUoFFZ0vf/Dk3g08KKojrgI/277OuTQOKci1ILU+GByS+qxKapmbXQXtaYYOSgSSaetX0gVgsN5Vd2ejElGX3u58T529PffNKXnMFGhvgT6GfqVvd+TA1lWXjUJUQGYjK5FxDjQhTlmZdpEkjls3Dn9rF/sEWTabDR0tiMEx4vrXA5PPhU5NA5RV46waZMCKTU+BItNpaCs/yh+pabII+9KWXUzDc3WbvlfHMihcZr5xlyWoU9bH6RIkkRyTBDHS/qP41KtKbZHkPlom/JCVQU5J2oZHtv9SDWnQ1mjWVl4sA9B/kZy+1m/eGJVWVrVTItZ6Vm/hMUBArW2xP2CgVOm/tIvQrEi6ss80i/Hi+33nNxTBVNb0qfFx7wKpo9Ijg2itKqZ5tb+se/CUxFkJdXNtJhtPUrBLgWE24uPaaRgnIq/uJ8MZI4IMg/4xRyD97Ae9IvsCCHXKMIvNNBEcICR48X9Y2VpjyATHusXX5OOmLDu74HShcWDakOtLzv7wRrhVTB9RFJsEALILe37wUy01bz3xJI/t23w7pGCkdoiyTQayMDeL6XVzTT2g4zXoqHcYxFkx0vq8TH2bCCTgyPtNUgGieL3pL8yt6SelPhQl8P5T0XrEHKXZOjpib/+9a955513nLnDvHQPxwCb2w9eGnsEWYtHIpVyS+vxNZ095Xhn6DSOJEtuM8f0h7xkigcHsrySeoZFB/ZoIJPayjto6ehPjg2irKaFxhb3ZwzoLp6KILPaVE6UNzKik2qvZ0MOjQGkPnX092oFs3HjRpYsWeIuWQYV/j4GosL8+sWszJMO/ty2AmM93f8jh8YhWupRW7VJtZMUbfdBZRfUaHL97uBIi6O1Kcba5kDvTeVIrVOTOGTL6wd+GLWmGDk4SvMIssKKRvsGy4Selc6Q9CakoIg+VTA9Tov8hz/8AQCrte9nFAOV5JggDuZVa1IKuDt4SsFYbQonyhu5fHRUj69x6rJfjh3lLtGc+PkYiAn3I7uglrkT+7amuaciyE6UN6KootcKxnZsF8LaimTwcaN0dpJigpCwm/LGJffthkvFQznIHJPP1MQQ6GFyXF0f5yTrtQ/GkTrfS/dJjg2ivslCdb25T+VQa4qQfAKRNc5BVtA2kI1I7HkxMznsZE4yrUiOCSK7oEaTUsDdQW0LvNAah4M/uQcRZA60Tkbqa9ITM8S/z1f8wmbxWARZbkk9Qf7GXtXCsUeSlSGUvgkm6nIFU11dzQcffMCXX35JVlYWjY2NBAQEMGrUKGbOnMkVV1xxRjp9L67jeKGPl9QTHuz+WZ+rKDXFHvG/5LVlkB6RGIqw9uyBl/xCwOirrTkmNogdB0qpqmvtstKmljgiyAxx7ikY1RW5JfUE+xsJDex52nndKSHkWiQjBUiKCeSno1V9uuJX69oiyDxhTi6pJyk6sFf3KofFnZKMVPtyD2e039kHf/7zn7niiivIzc3lqquu4qWXXmLz5s289NJLXHXVVeTn53PFFVfwpz/9yZPynlMkRAag10l96uj3ZATZ8WL7QNYbZSpJErrQeI84+vtyn5I9gszmmRxkJfV2E1QvBjIpMFLTZKRgz0vW2GKlss79NYFcxekX07hfmttSSSX1YlUJ2oeQn41OVzCRkZH897//xWg8M//NmDFjuPTSSzGbzR2m1/fiGnqdTGJUYJ/G93s0gswNAxnYX25r7veazWTjIwIw6GWOF9czpRf+ot7gqRxkza02SqqamTqmd/epdTJSsJsuwT5R0aqE8tk4GUGm7XORX1qP4OQ99xStk5Getf3OPli1alWHyuVUTCYTK1eudLtQg4nkmCDyyhpQ1L7Js+QpB39zq5XS6t7PyKBt2W9uQrRoo5j1OpnhccF9unPcsXdB6wiy/LZ9WO7qFy0HsrgIfwx6uW/7pabIIxFkuaV2c3JPNr6eitbJSM+Gy1FkhYWFHDlyhObm9hmAL730UrcLNZhIig1iy55CiiqaSIzSNk1LRzgHMo0VTF7bC9PbGRmcFknmp03RthGJoXzyTR42RXW5+JY7UWuKkQLCNY8gc5gBh0W7p19sOTsR5iYkk79L5wghqKmpwGJpBc4eVHHf5fa+Ly3N742oTsrLZdRuTO7U4TNB1tPgpvY7IyXCyv+5Mp7GumIOHiztloyno2b8DFQbdS7IrNPpCQgIwdfXtf47Gy4pmOeff56//vWvpKSk4ONz0n4uSZJXwfSSU+39faJgaj0TQXYyFUnv77FdTrJ4bZzgIxJD+fCr4xRVNDE0um/6xTOO5AYiQ30J8O39jNzhL1JqitFHp7p0TmNjHZIkERUV71LlRWN9Kw3NVqKiAtxiHtXrZWw21wZvIVTUShXJLwTZv+eRkK5gLW/Ex6gjIsS3WzJ2hNpUg2iuRR6S0OV3LITAarVQW1sB4BYl45KC+ec//8l7771HSkpKrxv00p7IEF/8ffTkFtcze6Jn6nufiqciyI4X1xMV6ou/G8reSr5BSD6Bmi77Rw61DyC5JfUeVzCejiAb6abS3SdDlYvARQXT0tJIWFiUy2V9TQYd9cKCxaZiMnQvG3evUdr2/GlsHrMpKoqiYjS4Xv+lS/Rt17FZwdB5pKAkSRiNJkJCIqirq3SLgnGpV0NCQoiL8/zgNxiQJImk2L7JrCyEsJuZPBWp5AY7P9i/Mzk0VtPUJFFhfgT4Gvpk34Wo90wEWU2DmZoGc682WJ6KFBAOelO3IpZUVUGnc32/t7FNqZgtnk9RJWx2BSPp3TTwd4LFar83k8E9plmHv0gorm2KNxiMKG7aN+PSHfzmN7/ht7/9Lfv376e4uLjdj5fekxwTRHFFEy1mz26GEk01YNU+gqymwUxto8VtAxnYwy/V6iLNNkNKkkRybFCfOJQVDwVeONKuuE/xyz1KGdMdU5deJyHLEmZrH+RAVCyApPkKxnFvRr2bVmg6AyCBi/V63BmZ6ZKCsVqt7Nixg2XLljF37lznT2ZmptsEGcwkt2VWzi/VJr9WZ6i1nonpdwzSblUwYXFgbUE0VbvtmqeTFBNEcaXnFb8zsi8kRtN2jpfUo5MlEiMD3HZNrXOSSZKEyaBzzvLdzVVXXcq1117Jz39+DatWLWfLlk9Pfmizgk5/xgD8+OOP8u67b7pNBrNVxWjQkZ19mLVrH+r2+aqqcuutN1Febk/T//HHGzlRWu5cwXzzzdc8+eTjbpO3K1xSMGvXruW+++5jz549HDx40Plz4MABreUbFDgGXk+byTwVQZar0UAG2qUmgZOK39MJFj0VQZZXUk9chL/T7OQOdGGxiJZ6RKt21SdNRh1Wm4qiarN6feyxdfzrX6/z29/+jiee+B21tfbM2kKxaoI7s0AAACAASURBVG4eE0JgsSoYDTKjRo3hkUce6/Y1tm7dQlJSMpGR9r06mzd/RGFpmXMFM23ahWRlHaaoSPvNly4pGEVRWLp0Kf7+/uh0unY/XnpPoJ+RyBBfj9v7PZWD7HhxPfERAW4dyBxmPU9k8PW44vdAZgVVCHJLGtwSNn4qcogjkky7fnE497VaxTgYMWIUfn5+lJQUYbNZ+duLL/LzW29j1arlPPvsUx2WKmlubuaJJ9ayatVyVq1azquvvgxATU01119/rfPn4ovn8vTTfwZg7dqHuOmmVVx33dX8+tf301Bfj8mg44cfdnPTTavOaENVVe6993befvsNAHJzj3PllYudK5YPP3yf+fMvAmDTpg85cuQw//f557np7nv4/vtvAZg7dx4bN37g9u/sdFzyrt144438/e9/Z82aNX2a9fdcJik2iGwP1yCx1xXX1v+iCkFeaQPnj45063Vln0Ak3yBNFUyAr4HIUF9ySzxnuhSqYo8g0yj82kF5TQvNZluvN/KdTrtIMtK7ff6O/SV8va/r0ssCgdmioNfJ3dqjdOGEGKaPd93s+MMPu7FYLMTHJ/Lhhnc5mpvHP/7fP5B9/Ln//rv48MP3ueKKq9qd8/LL/0BVVV555U2am5u45ZYbGT48lWnTpvPyy/8BIDs7iwcf/BVXXnk1AHfffT8hIfZIvmf/9iwfbXid++65t1O5ZFnm4Yd/z803X8+IESP585//l1/84gEiI6Ow2Wzs37+PMWPsz88llyzh4483suKqZUwbm+LcuDtu3AT+9renXf4ueopLCubf//43lZWVPP/8884vwsGXX36phVyDjqToQHYdKqO20UxIQM+TDrqKPQdZMYbUaZq2U9aLEslnQw6NQ9HQRAb2VYwnFb+orwBV+wgyR/47d69gJP8wTctaA0hISJKEqlGAx0MP/Qqj0YS/vz+PP76OwMBAdu/+jovmzMbo64ekN7Bo0aVs3771DAWze/d33H33/UiShL9/APPmLWD37u+YNm06AGVlpTz00K945JHfExdnzxP2yScb+eyzT7DZrDQ1NxMVHY9B37XiDA0N49e/fpi77lrDVVet4IILLgSgtrYWg0GPyXRavj/ZPtQLxYpkMBEWFk5FRbk7vq4ucUnB/PGPf9RajkGPI5Int6SetNQIzds7GUHmIQe/myKVTkUOjcWavVPT7LpJMUHsOlRGTYO5V9mGXcVTEWS5JXYzTOwQ9+zYduAIIe+pb2z6eNdWGRW1LbRaFBLc6Ndz8Nhj60hObr/nTwjF/oydEkHW8TMnOP3fjuOamhr55S/v5dZb72TcuAkA/PTTXjZseJfnnvsnoaGhvPXeBrZ8ttGl5zk7O4uQkJB2isJkMmGxdBAtJusAqS0SDiwWCyaT9s+zS+vLKVOmdPrjxT0kRgUiS5LHwmI9lYMst7jBPpCFu3cggzbZNY4kS/ZwJUXVQ0kuHRtIe1Ii+WzIIdoXuTIZdCiKik3xTA6/jPPO4+Ot21AUBZvNxscfbyQ9/czxLz39fDZu/AAhBM3NTXz++Wekp0/BZrPx0EO/YuHCRcyZM895fENDA/7+AQQHB2M2m/n8v5txpUsOHTrAu+++zcsvv05tbQ0bNrwDQGBgIKGhYZSUnFTw/v7+NDc32TdZti368vNzSUlxbTNsb+hUwbzyyisda8JTsFgsvPLKK24XajBiMuiIj/D3WOp+jw1kpRoOZKemjNGIxKgAZEnymKNfrSmyR5BpUBXSgU1RyS9rJMkNaXs6whFJpjRr950Z2zYhemo/zOL5mQxPTuaGG67lhhuuZfjwVC699Iozjrv++v9BCMF1113NLbfcwMKFi5g69QL27/+JPXu+57PPPnY6+v/zn1eYOvUC4uLiufbaq/jF/XcxLCm1XTnxjlYyDQ0NrF37EA8++AihoWE8/PBj/PvfL5OTcwSAmTPnsGvXN87jlyxZyssv/4P/ufcX7DmcDcCuXd8wa9Zcd39NZ9CpiayyspL58+cza9YsMjIySEpKwt/fn6amJvLy8vjuu+/Yvn07l112meZCDhaGxQSxO6vcIwWV1JpizSPIbIpKQVkD8yYnaHJ93amhygkTNGnD2Kb4PbeCKdZ8VVlU0YRNUTXxi8FJxW+pPAG+iZq0YdTrQLJHkrkj/ZCDd9756Iz/CaGiQ+WONbd1mIPswQcfdf7u5+fX7m8HaWmT2b79uw7b/N3v7OXnG5otVNW1EhdhN/vV1NQQFBR8xvGBgYG8+eYG599Dhgzh3Xc3Ov9evvwaHn30QS67bCmSJDF9+gymT5/h/LyurpYjR7L4P//nNx3K4046VTD33Xcf119/Pe+//z7vvPMO2dnZNDQ0EBQUxMiRI5k1axb33nsvoaHuS/qWm5vLAw88QG1tLSEhIaxbt45hw4a1O0ZRFB577DG++uorJEni5ptvZtmyZW6ToS9Jjg1i+0/FlNe0EBXmp2lbigdCYU+UN2JThCb+FwDJJ0DzSDKw+4++P1yOKkS72aW78VQEmWM15m4Hv4P/3957h8lRXfn7b1V1mJxzlEY5SyiDECAJBaxok5YVNrYMNtl4+RmMMSLYa4PXJnox9rL+EtYY2BUCRBISQUggIUBCo6zR5Jxz6O6q+v3R6taMpmemZ6aqpmdU7/PwPEx39a2jrq576p5z7ud4flfOqiLI0MfBiKKAzSLR4TQgRObyaJDpuwemwykjigIWSWDr1v/jn/98mZ///J5+j5OSksq1126kpqaauLju+dzS0hLuvvteQ9rd95rkj4mJYdOmTWzatEl3QwA2b97Mddddx7p163jzzTd54IEHuoXg3n77bQoLC9m+fTv19fWsX7+ehQsXkpZmfDtQrem870JPB3O2guxC3c4BnXfw6ycW6a4k09nBJEfw6UG340/S87o0VhpWQRYeYtWtTbenksxRVQT6+BfAHSZraXfpvuL37IAXLDpLxDhk7FYJQRBYv/57rF//vQGPtWTJsh7fmzRJfxFVD34l+fPy8ny+/vXXX2tmSE1NDUePHmX16tUArF69mqNHj1Jb2zWB++6773LVVVchiiIxMTEsW7aM999/XzM7ekN1tOk6fkpcCDYDGioZpUGWV3ZmIovQL5/gqVjSS5MMzjp+va+LbFA7Xq06i/aEp5LMUa3vTnG7VUJVVP0T/QZokCmKitOlaLoZORDwy8FcffXV/OMf//D+7XQ6+cMf/sDtt9+umSFlZWUkJiZ61QEkSSIhIYGysrJux6WknJ0Yk5OTKS8v18yOnlBa6mj+f7fQXnxct3NIokhmUrjuE5lhFWRlTbpOZOCpJGvXtZIsJS4Em1XUvQDDiMKLtg4XpdUtuuVfPEjRqTiri3Q9h2dHv+5hsh40yLREawXlQMHvjZa/+MUv+OSTT/jhD3/I7373OxISEti6dWvfHw5AYmP7Xzsvhwm0oNJefIL4BRN1sMrNlDFxvLMnj+iY0AF3UoyP7z0kVZ9bTRuQMG4CUog+4avWdidlNS1cOjvdpz192egvbaPHUbYbwpU6QuJHaTKmh842jkuPpqi6RTO7fVHRWoESGU9Cin/7oAZiS/bpalRg5sREXf8t9elZ1J74jJhQFSmkZ2dWWSli6WNTYU9IkoggCDhdyoDHAPr8rEN2IlrtgzpHXzjb3Cvw0GAbktTdkel5bl+IoqjJ78MvBzNx4kRef/11rrzySn70ox/x3e9+l9/+Vls1zuTkZCoqKpBlGUmSkGWZyspKkpOTux1XWlrK9OnuqqFzVzT+UFPTjDIAoTwhJApHVRFVVfpJhyRFBeF0KRw8Wj6gRlfx8eF92tdelIsQFE5tiwAt+vxbjhfUoaqQEGHvZo8/NvqLKriLTOryT9ESoV1DvHNtTIsLYefXJZSVN+jWQrmtvAAhItmv72ag3+GBY+7VfkyIRdffscsaC0DlqRNYUnp+IFMUZVDdGm1WkbYO14DH6KtbpKoq7hyMPXRQdvZFW7sLiySiqiouV9e5abAdLQeCoihdfh+iKAzowdyvO6WiooKf/OQnWK1WfvWrX7Fz507+9Kc/4XJpJ2MeGxvLpEmT2LbNXW63bds2Jk2aRExMTJfjVq5cyeuvv46iKNTW1rJjxw5WrFihmR29IUanuCtjdGSUAfF+2YBSWCMS/HC2kkyu1V8yxiUrlFS16DK+u4Ks3JCwZVxkEOEh+lZEefco1et7XexWCYdL0S8H56kgM6DJ2EjLv4CfDmbdunXMnDmT1157jY0bN/Lmm29y+PBhvve9gVc5+OLBBx/k5ZdfZsWKFbz88ss89NBDANx4441kZ2d7bUlLS2P58uVcffXV3HrrraSn67PP4lzE6FQc1cW6JpTjI4PcnRR1cjCeCjK9J7Lcskbio/SfyOBMD5J6/SvJQD9lZW8FWYz+FWRZOpWNd0YIjUGwh3hbQuiF3SqBquLQ6Anf0w/mhhuu47rrvsejj/07LpfL2xXSH6qrq7j99p/4fbxHkWCg+Zf29nY2bbqetjZ3EdJrr/2DurqzOcktW17npZf+PqCxB4tfIbJnn32WWbNmef9OTEzkv//7vzXfxT9mzBhef/31bq//7W9/8/6/JElex2M0YnQqqrMdtbkGITxOl3MIgsDoZP06Kaotte4KMr0nsrJGxqZ23ySmB2J0Ks6Tu3UtV4074/jzShu5bJb2351XgyxKvwR/Q4uDmsZ2ls3Rv6RfEARscWm4dC4h79xC2a7RCsCjRSbLMrf89Ifs2vsly9Z3Db8qioIgCD5/b3Fx8Tz99HN+n6/Dm+AfmP3/+7//5NJLlxAc7O4f9NprrzBnzjyio93Rn7VrN3Dddd/ju9+9itBQ7bXbesMvB9PZuXTm+9//vqbGBDqde5CIOjkYcIeVDufV0O5wEWTzv1+5PyhGTGTNHdQ2djB6jv5PynDmupypJBPCYnU5h9fxl+vj+I2oINOjs2hv2OLScJzY36/POE/uwXlil9/Hq0CEQ0YVBVr7SIRbJyzGOv4iv8d2OBw4HB2ER7grIZ9//jlKSoppa2ulpKSY++9/iH//94e9x5eUFHHjjbdw8cWX8OMfX8877+wEYNGiOdx00y3s2vUJDQ0N3HrrHVx6qbsb8EMP3U9efh7t7Q5Gj8rkl798gIiIrtfn4MFv+OMfH+Pvfz9byfujH23k9tvvYtas2bz11hs89dRfAHjhheeprq7yKkJv3vwbRo/OYt68Bezc+SFr13aXt9ETbWevEU4XaZKMGbqdZ3RyBKrqbqE8IUM7pQQ42wFSzxWMp3+KUROZN95fW4Kok4OBM44/t4a2DhfBdq0dfylCeJyuGmT5ZY0IAmQm6psX82CNT0f99iOUtkbdJIkEQNRYut8zOZeUFDN35nTmzZ3vfe/gwW/47//+H2/bEk+Pl717P+fpp//E8uWraGtr7TZmaGgo//VfL3Lo0EEeeOCXXgdz551306HYcMkK77zxIv/zPy9w881dt3/MnHkBbW1t5OScYuzYceTm5tDc3MTMmRdQUVFOe3s7SUnuYqgf/GATb7+9tZsi9JQp09i7d4/pYAIZwR6KFBZtyM5xcE/U2juYEoTgCMQg/SaZ3LJGREEwbCLzOv76EsjQR5MMzrZQLqzQ57oYkRdLjQvDbjMmmWyLc+dGlbpSvx2MdfxF/VplADiaOmho7nArkmsgquqZnNvb27j/njt47c23ufZ6t5rJwoUXdeuJderUCf74x9/zpz89Q1RUlE8Hs3SpuxBpypRpVFdX0dHRgd1u5733tvHOe++gyDJORwfp6b6lD1au/A7vvfc2t9/+c955522uuGINgiBQVVXpDYX1RmxsHJWV+vd/OZeRtavHAGxxabprX0WE2oiLDNIlD2OEBlnemV7vRk1k7kqySN0ryUZ1cvxa4q0g0zFsqarqmQS/MU4fwBbvniz1vl9sOrVQtksiC+fM5quDB7yvBQd3lQqqrKzg/vvv4de/frhH5wBgs7mLXTwbyWVZ9vaC+cV9j/Lc317mxhtvxuHo8Pn5lStXs2PHB3R0dLBjxwesXPkdt412e4+f6YzD0WFI/5dz8cvB3HzzzT5fv+222zQ1ZjhgjU/XXZoE3JOZ1g7GW0Gm80SWX9aoe3nyubglY3R2/CFux691JZnSWKF7BVlVfRst7dq3SO4NKTwGrMEDbj7mL3adpPtlZwffHjlKenqmz/dbW1v4xS/u4qabbmX69Jn9Hr+pqYmQkFDCwiNAlXnnnbd6PDYpKYlRo7J44on/YNSoLG9ILCMjk5qami6tVUJDQ2lubu7y+fz8PMaOHd9vGweLXw5m3759Pl//8kvf8tMjGVtcOrg6UJtrdD1PVnIE1Q3tNLb03pOnP7gryNp1zb9U1rknsqwUYyrIPLhLlfV3/KOTIzSXjFEM0CDTW0HZF2e7W+rr+CVJxCKJmjmY+++/hxtuuI4f/PgGFEXlhz+6yedxn376MYWFBbz00t+9PV62b3/P7/MsWHAhiUkp/OJnN/CrX97FhAkTej3+iivW8Pbbb3DFFau9r9ntQVxwwWwOHDirC3nlldfy7//+MDfccB15ebkAfPnlXi69VP/+L+fSaw7mySefBNzaY57/91BUVNTvHfQjAVu8J66sfyUZuMNNM8Zqcx4jNMhyS42tVPLQWZNMr0oycP+79h+vpLHFQUSoNnt8jKjsyy1txGYVSY3XvrNob0jRqbgKDvR94CCx2SQ6HIN3MJ37wSgNFaiyC+lMT5ZNm7rubVm1ajWrVq3GF54KMoDdu7/q8l7nv++8ezMAyX50fF2x4gpWrLii2+v/+q8/4B//eJH58xcCsGbNetasWe99v6AgH0VRvG2ajaTXFUx5eTnl5eWoqur9f89/ycnJ3ZzO+YC1U+JSTzKTwhEEbXf0G1EKm3um13uqxr3e+8KzKtN7Y19nx68VSm3JmQoy/WLkeWWNjEoMRxIN1rSKTkVtb0Jp01co1G4RNW+hrMpOXSX6VVXF4VQGvX9n2rQZXHjhxd6NludSWVnOv/3bvYM6x0DpdQXzu9+5O63NmjWLq6++2hCDAh0pOAwhJEr3SrIgm4WUuFBNE8pyrQEVZKX6tUjuDSnq7B4lPSvJOjt+LVeWeq4qXbJCQXkzS2frW9zhi7N7x/yvJBsInoISh1PWRCtOVRWQnWDXb2OiR+JGC4mYziuWc5k7d8Ggxx8ofl2JCy64gOrqagBaWlp46qmneOaZZ3r0mCMdI+LKgHdHv1Z5BaWuBDFGv13cTpdCUWWTIVIk5+KtJDPA8afGhWqW6FdlF0p9OZKO18XdWVQxPC8GnfYo6V1JZtFYut+rQabfCsYT0tNKgSAQ8cvB/Nu//RuNje4b6tFHH2X//v0cPHiQBx54QFfjAhUxOtWQSrKs5Aia25xUN7QPeixVVXR/Uva0SDYykdwZMSZV99AluB1/flmTJtdfaagAVdb1uhglPOoLITT6TCVZzw5Gi+9RFAWsFkmzRL8qu4tr+qNB1l86t0gOJFRVwb2FdfD4tdGypKSErKwsVFVlx44dbNu2jaCgIJYuXaqJEcMNMTrVW0mmlyYZdO2kGB8VPKix1KYacDl0TvA3AAzJCgbOaJKd+Ez3FrqjkyP47FAZVQ3tJAzyuhhVeBERatO1s2hPCIJwxvH7djAWi42WlkZCQwffmM5uE2nVqoWyy4neXSwdzrMtkgMBVVWRZRdNTXXYbNr8VvxyMDabjebmZk6fPk1SUhIxMTG4XC46Ovre4DMSMUqTLDXe3XQst7SReZMSBzWW5wbXs997XlkjkWE2osON39AFZ6qwdBYjhU6Ov7RRGwcjCIhRyX0fPEDyyhrJ0rmzaG9IUSk9VpJFR8dTV1dFc3P9oM/T7pBpbnNS6rL7bNrlC1EUUZTuYTWlrREUBVEuHLRdvlBUqGtsJ8RuodzR+zTck416IIoSwcFhhIVpE071y8GsXr2aH/zgB7S0tLBx40YAjh49Slqa/qqsgYhRmmQWSSQzKUyTiiW5zt0fXc89MLmlQzuReSvJDHD8VotIXlkj8ycP0vHXFiNEJCLo1G/E3Vm0lQVTknQZ3x/E6FTUE7t8apJJkoW4OG2ca2FFE3/6x35uWjuZBZP9+/f21Lit+R//hpQ0juAlP9XEtnM5ml/Ln944zc+vmUFWUu9l9Vo26DMavxzMfffdx+7du7FYLCxY4K5IEASBX/7yl7oaF6gI9lBDKskARidFsOtQKbKiDKrEVKktcffosIX0ffAAaG5zUlHXxqLp+j2J94WRjj8jURvHr9SV6LyqdE9MQxW2hHMcv46VZClxbsefX9bkt4PxhepoQ22uQYy+TEPrujJU+8WMxu8Za9GiRWRkZHDgwAFKS0uZNm0aCxcu1NO2gMawSrKUCBxOhdLq7gJ6/cFdQabfRJY/BDvFz8VQx58cQUF5E/IgQheqy4HSWKHvqtKT4B9A+22tMKqSzCKJZCaGD7rCz6hwcmJMCKFB+uV4AgG/HExlZSUbN25k+fLl3H777SxfvpyNGzdSUVGht30Bi5GVZDC4jX2qoqDU69vFMre0EQEM1bryhVGOPys5AodrcC2UlYZyUFV9K8hKG0mODSFkCCcyISTKEE0ycDv+wkE6fm/zN50cv6qqZ8LJQ+f0jcIvB/Pggw8yceJEvvzyS3bv3s2XX37JxIkT2bx5s972BSydK8n0JCE6mBC7ZXAOprESZJeuT2S5ZY0kx4Vq3ielv5x1/PomRUefCTnllw88Nq7UnsmL6XRdVFUlt6xxyMMwfVWSacnolPDBO/7aErDYdCsUqWvqoKHFMST7kozGLwfz9ddfc8899xAS4o7fh4SE8Itf/IIDB/TXGApUOleS6Ym7k2L4oAQWZZ1LYc8+kQ19PPms46/t++BBkBAVTGiQxRtLHwhKXSkIEmKkPgn4mka3WOpQ5l88SNEphq1gYHArfqWuBDEqBUHQR1bH85sJhOuiN359g5GRkZw+fbrLa7m5ud1ae55PdEko68zolAiKq1oGvIlM8VSQ6aRBVt3QTnOb0/tUP5QYFe8XBGHQLRXk2mLEqEQESZ9Vn9GdRXvDKE0yj+MftIPROS9mkQTS4vWToQkU/Ppl//jHP+aGG27gyiuvJCUlhdLSUrZs2cKdd96pt30Bi9EJZUVVKapoZmxa/5fVbjHFeN3a8XqfyAJgIpM6ryx1rCQD93V594sCOs5smOsvSl0JUvxoHSxzk1vagEUSSU8Y+omss+PXs5LMveKPGLCGn9rRgtpar284ubSRjMRwrJaR3+/Rr3/h1VdfzeOPP05dXR0ff/wxdXV1/PGPf+Saa67R276AxkhNMmDA1TFKXam+CsqljVgtxkvB+8JYxx+OoqoUVvR/MlOdHahN1boXXmQmhmki/jhYjFpZgrvQpKSqZUDy/WfDyfrs8ZMVhfzyoc+LGYXfa/OFCxee12XJvhCjU3Ee36W7NElUmJ3ocPuAlv2q4kJpKMOWqd/TfF6ZW0E5ECYyOJvo15vOO/rHpUX1cXRXlPpSQNUtFCMrCgXlTSyeGRg9m4SQKLAZU0mWdWbFX1DRxPj0fl6XWn03JJdWt+JwKudF/gX8dDAOh4M33niDY8eO0dradT/GY489pothwwGjNMnAfdN4tL76g9JQAYp+YoouWaGgoonLZhkvBd8TYnTKGcev6JaoBbfjj4mwD2hlqfdei+LKFhyuwJnI3N0tjaokO7PiL23sv4OpKwFrEEJojB6mcXqI9fqMxi8Hc++993L8+HEuu+wy4uL0nUiHE0ZpkgGMSY3k65NVNDR3EBnmv9aXpwGXXg6moLwJp0thbGrglFx2dfzxup5rbGokp4r77/jl2mKQLAgRCTpYBTklbpsC6bpI0Sm48vWvPI0MtREfFeT9DvqDUutWHNcrIpFT3EB4iHXQGnbDBb8czGeffcbOnTvP66oxX3iePuXaEiw6J5THnUnunypuYM5E/yclvcUUPZPruAEUH+hFl4SyAQ7my2OV1DS0ExvpfxGFUleKGJWMIOrTC+RUcT3R4fYhUVDuCTE6FfW4b00yrRmbGsWRvJp+h6+VuhIsmbN0syunuIGxqZEBo6CsN37FD5KTk3E4HHrbMuwQ7KEIoTEotUW6nyszyV110t+nMr3FFHNKGoiPCurXqkpvPJVkcq3+8X5P7uVUSf/UgJXaYt0SyeC+LoE2kRmZ6B+XFkljq5PKev+bIiptjajtTbrlXxqaO6isb+t3vm4449cKZv369dxyyy18//vfJza2q/Ln+Z74F2PSvIlBPbFIIqOTIzhV3M+JTEcxRVVVySmuZ8ro3tVgjcZTSabU6z+RpSWEYrdJnCpu8FtgUXW0obbU6ha2rGlop7axg5XzAmdVCec4mJRJup7Lu+IvaiAx2j+BV7178wTial9v/HIwL7/8MgB/+tOfurwuCAI7d+7U3qphhBSbjqPkCKrs0m3DnIdxaZG8v6/Q730XbjHFSixj5uliT2VdG42tTsalB94NY1QlmSSKjEmJIKcfeRhvgl+nJ2XPairQnpSNrCRLjgslxG4hp6Teb4Vvb75Sp/bVOSUNWC0imUMoPGo0fs2IH330kd52DFvEmHRQZJSGMqSYdF3PNTY1EllRySttZGJmdJ/Hu8UUFf2fyAIokezBXUn2qe6VZOC+Lm9/nk9bh8svLTa9pXtyihuwWyXSEoZ+X1JnjKwkEwWBsWn9K8BQ6orBHooQrM/v+VRxA6MDqJzfCALiX9rW1sbPfvYzLr/8clauXMnHH3/s87h9+/YxY8YM1q1bx7p167jqqqsMtrQ74hmnotTon4fx7OI/5WceRu8lf05JPSF2C8lxgTWRgaeSzOFuFa0z49KjUNWzJah9obeYYk5xA2NSIwbVP0gvpOhUlNoS3VXIwb3iL6tppbnN6dfxSl0pkk4VZB1OmcKKJsb1s2x6uNPj49Yll1zi1xf9ySefDNqI559/VuQQMwAAIABJREFUntDQUD788EPy8/P513/9V7Zv305oaPeJa8yYMWzZsmXQ59QKMSoRRIsheZjQICupcaF+52GUmiIQ9RNTPFXcwNi0SMQASiR7kM6EOZTaYsQIfSvJ3F083fH+qX7ko5TaojOlsNo7gLYOF0VVzay5cJTmY2uBGJOGevxT1LYGd8hMRzwl2jnFDcwc17szV1UVuaYI6zh9csp5pY3IihpQZeNG0KOD+cMf/uD9/+zsbLZu3cr111/v1SJ7+eWXWb9+vSZGvPfee/z+978HYNSoUUydOpVdu3axatUqTcbXE0G0IEYnu/c1GMDYNHdZrKKqfU7sck0hYnSKLrmh5jZ3K96FQ9iKtzfccXQBubYQyyj9yk4Bgu0W0hPC/Krwc09khVhHz9bFltOlDahq4OVfPJxd8Rci6uxgRidHIIkCp0rq+3YwTdXgbPPapzWeqMMY08G4mTfvbGL44Ycf5vnnnycx8Wz/8cWLF/PjH/+YH/3oR4M2orS0lNTUs2Gc5ORkysvLfR6bn5/Phg0bsFgsXHfddWzYsGHQ5x8sYkw6cukxQ841NjWSTw+WUlrVQlofIoZKTRFS2lRd7MgJ8IoYwRqEEJGAUl1oyPnGpUaxO7sMl6z0GmNXW+qgowUxJkMXO3KKGxCEwN0pLsW6J3C5phBL+nRdz2WzSmQmhftVgCHXun8nUpx+1yUlLpSw4JHdwfJc/Hq0rays9PaC8RASEuJ3R8sNGzZQWuq7cuTzzz/3awyAKVOm8OmnnxIeHk5RURE//OEPSUxM5MILL/R7DIDY2MGpy8bHd60Cqc8YS+2pz4kJBSlE3wqRBTNEnn/nGGUN7cya4rs6Jj4+HLmlgaa2BiIyxxEVr71NJfsKsUgCc6enDkhJ+NzvUA+UlCw6KvIGfK7+fG725CR2flNMs1NhXFLPTre1/iQtQOzYiQQN8jvwZV9+RTOjUyLJSOu7CMQIutsYTntkPNbmMkN+A9PHxfPOnjyiokOwWnz/TuPjw6k7Vkk7AgnjJiLatN2cqigquaUNLJqZashvMZDwy8EsWbKEm2++mZtvvpmkpCTKysp47rnnWLJkiV8neeONN3p9PyUlhZKSEmJi3Po/ZWVlzJ8/v9txYWFnHUN6ejrLli3jm2++6beDqalpRlEGlmSMjw+nqqqreq7L7t5ZX3nqGBad6/tFVSUy1MaB4xXM9bHs99jnKj4KQJs9AWfVwLsu9sShU1VkJIbTWN/a98Hn4Os71ANnWDKu419QWVKJYOufNEd/bUyIcG9k3Z9dRlRQz7dVR94JABrFWJoG8R34sk9WFE4U1LFoWrIh329f9PgdRqXRVppniI2pMSE4XQpfZZf5bHXhsbGtKAchMpGaBifgX1GAvxRXNtPS7iItNmRA/2aj7pfeEEVhQA/mfmUZH3roIWbOnMnmzZvZsGEDDz74IDNmzOChhx7q9wl9sXLlSl599VXAHQLLzs7m4osv7nZcZWWlt/qkvr6ePXv2MHHiRE1sGAxip4Sy3gie8sui3pf9nqo2KVb7Jb/TpZBX1hTwCUtPOMaI6xITEURshL3PAgylpsjdm6efDs8fCiua6XDKA+oZZCRibDpKQxmqS391kLMSS71fF7mm0Pt70RrPuQM1nKwnfq1g7HY7d999N3fffbcuRmzatIl7772Xyy+/HFEUefjhh72rlSeffJKEhAT+5V/+he3bt/PKK69gsViQZZl169axbNkyXWzqD0JwJEJQuCGlygAT0qP4+kQV1fVtxPUgmifXFCKERiMEad9sKre0AZesMCEjMBPJHsQzzlWuLUJKGqf7+canR3O4D/0rRceJ7ESheyIL+OsSkw6qqnvDNYCIUBvJsSEcL6xn1YJMn8eojjbUpirECd0farXgeKFbFy7+PBG47Izf5UV79uzhnXfeoba2lr/85S9kZ2fT3NysiVRMSEgITz31lM/3OnfN3LhxIxs3bhz0+bRGEATEmDTDKsk8myyPFdZxcQ8/WqW2SLeKmGMFdQiC29EFMkJoDNhCUGqMSfRPzIziiyPllFS3+GyHqzo7UBoqsIzpHv7VguOFdSTFhBAVQLpwvvCsquWaQt0dDMDEjGg+P1LeYwGGfEZLUI8Ev6qqnCisY8romIDShTMKv0JkL730Eg8++CCjRo1i//79AAQFBfHkk0/qatxwQoxJR6krRlUU3c+VGhdKeIiV4wV1Pt9XZSdKXZku4TFwP5FlJoYTEhTYFTGCICDFpiMbtLKclHHG8fdwXZS6YkD1rqy0xCUrnCiq90vhYagRIuLBGmTYin9iZjQdDpn8ct95DM8DiB6VfaXVLTS2OpmYEfjXRQ/8cjAvvPACf//737npppsQz+wOzsrKIi8vT1fjhhNSbPqZneOVup9LEAQmZkRzvLDe545opa4UVBlRh1BMh1PmdEnDsJjIwB0mU2qLDHH8cVHBxEUG9ej4ZW9eTPvrUlDeRIdDZtIwuC6CILpFYg1aWXpChj1dF6WmyC0RE6r9d+d52BgO10UP/HIwLS0tJCe7S2I9yzyXy4XVGthPsEbiSfQb9bQ8MTOauqYOKuu6y5F7blw9VjA5JQ3Iijpsnsik2Ay342/U3/GD+7qcLKpH8eX4awrBGqyLRMzxQvdEFuj5Fw9SbAZybZEhkjERITbS4kO939G5yDVFSLEZuoSwjhfWExcZ1GOudKTjl4OZO3cuf/3rX7u89uKLL/osJT5fEaNTQRCNi/efmUiO+bhp5JoikGwIEYnd3hssxwvqkERh2FTEeFZxno10ejMpI5qWdhdFFc3d3lNqipBi03WRiDleUEdqfCgRIfr0/dEaMSYdHG2ozdWGnG9iRjQ5xQ04XV1Xsqoiu+WEdMhXKmfyL8PlYUwP/Pql33///Xz44YcsWbKElpYWVqxYwfvvv8+9996rt33DBsFiQ4xKRjbIwSTFhBAZZvO57Hcn+NMQdBA7PF5Yx6jkcL9UgwMBMSrljOM3bmUJdHtaVlUFWafCC5escKq4YVhNZGd39Bt3XRwuhbyyxi6vO2vLQHbokuD37H+ZmDk8VpV64NcskZCQwP/93/9x6NAhSktLSU5OZvr06d58jIkbMS4TueSoIecSBIFJGdEcza/tokt2VutqjubnbOtwkVfaxKoF+hQP6IHRjj863E5idDDHCupYMe/s9+TWumrXJS+WW9qIw6UMKwfjdrSC2/GPukD3803IiEIAjubXMr5T9aOjsqCTPdriyb8Mp+uiNX57CEEQmDFjBqtWrWLmzJmmc/GBFJeJ2lqP0tq/rpMDZcroGBpbnV3CMXJTrVvrSoeJ7FhBHYqqMnV0jOZj64kYm27YCgbc1+V4YR1Ol+x9TdZx42t2bg2iIDBpGD0pC1Y7QmSCYSHl0CAro5IjOJJX2+V1R0U+CBLimTbbWnI4t4bk2BBiIrSVnhlO9LiCWbVqFe+99x7Qu3S/1WolLS2N2267jdmz9VGIHS6IcaMAUKoLEQ1Itnom+sN5Nd4ueR0V7so+PUphD+fWEGSThp0irBiTgStnL2p7sy4bT89lalYsH31TwsniBqaMcl8jpaYQBEGXfu+Hc2sZkxoR8GXj5yLFpCNXFxh2vmlZMbz9eT7NbU6v6GRHRT5idDKCpO131+GQOVFUz5IL9OmOOVzo0cE88sgj3v/vLN1/LoqicPLkSe655x527NihrXXDDO8Gsup8LBn6KsUCRIbZyUgMIzu3lu8sHAWceSIDzbtrqqpKdm4Nk0fFDLuOfJ74ulxTiCV1su7nm5QRjUUSyD5d08XBiBGJCBZtN0E2NHdQUNHEdxdnaTquEYixGbjyvkJ1tCLYQvr+wCCZlhXLW3vyOZJXy/zJ7gIYR0UeYpL2clPHC+twySrTsvruDzSS6dHBzJlzNobfWbrfFwsWLKCoyLgQRKAi2IIRIhJRDH0qi+W9vYW0trsICbLQUXYaITJJc62rsppWaho7+M6Fwys8Bu7cGIBSnQ8GOBi7TWJ8ehSHO4Vj5Op8pGTtJzLPOYbjRCbFjwJAri7QXSQW3P1hQoMsHM6tYf7kRJSWOuTmOuw6qAkczq3FZhUZnz68Vvta06ODefbZZ7n55psBet2x75Fy+dWvfqWxacMTKS4Tucq4DajTsmJ554sCjhXUMntCAh1lp5EStdfdys51tx4ebvkXADEoHCE8Drkq37BzTsuK5dWPcqhpaCfa2o7aUuedULUkO7eGiFAb6Yn6h/60xhtSrsoHAxyMKApMGR1Ddp67MEapzne/roODyc6tYWJGdI8tAs4Xeox1dG74VV5e3uN/Jl0R4zJRm6pQO1oMOV9WSgTBdons3FqU1nrkphqkOD2eyNwJy7jI4blhTIobZajjn3pmRZGdV6PbRKYoKkfyapk6OiYg21b3hRgcgRAWa/gDWWOLg6KKZvcDhyBqXnhRUddKZX3bsFxVak2PKxiPFL+iKKxdu5bZs2djsw2PTVxDiXQmHCNXFxgS77dIIpMzY8jOrUGe3A6AqPGTsjth2cCSC7RPUBuFGD/aHe83KNGfEhtCTISd7NM1LJTzQRA0n8jyyhppaXcxNWv4rSo9SPGjkc84YCPwrMAP5dYQ35iPNS4NwaptXiz79JnV/jC+LlrRZ7ZWFEVuueUW07n4ydl4v3F5mBlj46hr6qAm7wQgeJ2cVmTn1uCSFWaMGb5PZB7VXqMmM0EQmDE2jiN5tTgr8xCjUhCs2parfnOqCkkUmDp6+F4XMX4UamOlYSv+yDA7o5PDOXiyEqUqD3vyGM3PceBUNUkxISRG61+4EOj4LRVz8OBBvW0ZEZyN9+cads6Z4+IQBYHWkhyscamaT2Rfn6wiLNjK+GGic+UL78rSwDzM7PHxOFwyzopcb75BK1RV5esTVUzMiBrWfd6lM9+LkdflgvHx1FZUoLY1Yk/StvquqdXBicJ6Zk+I13Tc4YpfO/lTUlK48cYbWbp0KUlJSV32xHTu12LiRorPQq48bdj5woKtTMyIJLipBLvGO/idLoVvc6qZOzEBaRhvrhXsoQgRCSgGxvsnZESREuzE4mjSPMGfX9ZIZV0bK+cNH1UFX3gdTHUelrQphpxzzoQETn7uDmPZU8aiZV/NA6eqUVSVORMSNBx1+OKXg+no6PB2jqyoqNDVoJGAlJCFK/dLlNYGxBBjyhQXZAUR9m0rLaFpaJlhOJpfS7tDZvYIuGGkuFGGOn5JFFmU5oAaUGO0DVt+kV2GAMwaP7yflIWgMITweHclmUEkxoQwJaIRRRWwJWRCvXYu5puTVcRFBpExDKv69KBXB6MoCqIo8rvf/c4oe0YEYoI7rqtU5SJmzjLknNMi3c2UDjeEsUDDcb8+UUWw3cLkUcNfT0mKH+12/G2NiMERhpxzcngjcrXAyeYwpmk47ueHShmXHkVk6PDPjUrxow11/ADjQxooa4gitk27dgGt7S6O5NWybE7aedm90he9xjwWL17MY489xsmTJ42yZ0QgxWWAICJXGpeHCWoqRkHgo1ztftiyonDgVBUzx8YOu937vvBU1xkZJot1llOpRLE/p0GzMctrWykobxoxcX4pfhRqcw1KW2PfB2uAqqpEO8spdMWy93CZZuN+e7oaWVFHxGpfK3qdNR588EGKi4u58sor2bBhAy+88AK1tbW9fcQEECx2xJh0Qx2MXJFDW3ASOeVtFFd270UyELJza2lpd42YeLIUPxoEAbkix5DzqaqCUnWalvAMvj5RhcMp9/0hP9h7pBxBcBcRjAS8K/4KY1YxakMFoqOFuqBUPvmmWLNx9x6pIDrcTlaKMavj4UCvDmbZsmU89dRT7N69m2uuuYb333+fSy65hJ/+9Kd88MEHOJ1Oo+wcdkgJo5GrclFV/Vv1qoqMXJVHaMZELJLIZ4e0eSr77NtSIkKsTBvG5cmdEaxBiDEZhoVjlPoycLQRkzWZtg4X35ysGvyYisru7DJmjU8YMSq9bscvIVca4/g950kcP5WjebWU1Qy+RLq2sZ3DeTVcNC15WG561Qu/4h4RERFce+21vPLKK7z77rtMnTqV3/3udyxatEhv+4YtUsIYd8e+Bv2LIpTaYnC2E5QyjvlTk/jiSHm3zn39paHFwaHTNVw4NXlEhMc8SIljkCtzURVtVhO94VkppU6aTlxkkCaO/2hBLbWNHSwb5tVjnREsNncvJYNWlnJ5DthCmDV7KqIosDt78Ndlz+FyVBUWTUvSwMKRQ79mDofDQXZ2NocOHaK6uprx48frZdewR4x319cbESbz3JhS4jiWz8ukuc3JtzmDa0X7xeFyZEVl0fRkLUwMGKTEseBsdztlnZHLcxDsYUhRSVw8PZljBXVU1bcNaszPvi0jNMjCgqkjayJzO/48VMWl+7nkihykxDFEhgczd1Iin2eXIysDfyBTVJXdh0qZmBFFgrm5sgt+OZivvvqKX//611x44YU8+eSTzJgxgw8++ICXXnpJb/uGLWJUMliDDXkqkytyEIIjEMLjmDE+nuhw+6CellVV5bNDpYxJjSAlLlRDS4ceKXEsgCFhMqUyBzFxDIIgcNG0ZARgzyCelpvbnBw4VcXCKUkjTkRRShwHskP3xnCqoxWlrgQpwf07WDYvg4YWB9mnB55bPlVUT1V9OxdP175p2XCnVwfz9NNPs2zZMq+q8nPPPccHH3zArbfeSmrq8NWlMgJBFJGSxiKX61+B534iG4cgCEiiwKJpyRzOq6FygE/LJ4vqKatpHZE3jBAejxAcgVx+StfzqO3NKPVlXmXrmIggpmTFsOvbUlzywJ6Wdx8qwyWPvFUldHL8Oj+QuSMKqvd8cyYlEhFq4+MDJQMe8+MDJQTbLVwwQqr6tKRXB3Pw4EHuuusudu/ezSOPPHLed6zsL1LSBJS6EtR2baq6fKG0NqA2VXlvGIBLZ6UiCgIf7h/Y0+AHXxYRFmz1NmUaSQiCgJQ4TvcVjGf8ztfl8jnp1Dc72He0/3k5l6zw4VdFTMyIIiMxXDM7AwUxLAYhNEZ/B1OR4xYeTXCHsC2SyJILUsnOraGkqv/3aXV9G18dr+KSmSnYrSNrVakFvTqY559/nu985zvY7dqqjZ4vSMnuHJVLx1WMpyJG7DSRRYfbWTAlkc8OldLY2r9dyiVVzRzMqWbJBakj9oaREsegNlbquu/CPZGJXpFNcCv5psaH8v6+QhS1fxv89h2toK6pgxUjKLl/LlLiWEMcjBid1qUh32WzUrFZRN7bV9jv8T7YX4QgwLLZ53dr5J4YOeVBAYgUPxoki65hMrn8FIhSNwXlKxZk4nQpvPtF/1Sdt+7OI8gmsXQE3zDimbCVnmEyufwUYmx6Fyl4QRD4zsJMSqpb2H+s0u+xXLLCW3vySE8IGzEl476QEse6N1w267PXTlUU5IrTSIldFZTDQ2xcOiuVL46U96tkuaahnU8PlnDh1KQRUzKuNaaD0RFBsrqFL8tO6HYOufQ4UsIYBEtXyZDk2FAunJLER9+UUO1nLuZ0aQNfn6hi+dx0wkOGvwRJT0jxo0CyIZce02V81eVArszx2SJ53qREUuNDeWNXrt+l5J8eLKWqvp0Ni7NG9B4LKXkCgG7XRanOB2ebz+tyxYJMbBaJLZ/6X/W59TP3sWsv0r7B30jBdDA6IyVPQKkuQHW2az622tGCUlOA1EO72Q2Ls5BEgX/s6PtJXVFUXvrgBFFhthEdhoEzjj9pHHLpcV3GlytyQHZhSe1+XURB4JolY6msb+O9fX2vLhtbHLyxK5dJmdHDuh+PP4ix6WAPxaXTdXGdcVy+7peIUBtXLMzk65NV3vbgvXGquJ49h8u5fE46sZHm6qUnTAejM1LSeFAVXWLLrrIToKpIPXTOjIkIYu2iURzMqe6zPHbb5/kUVjRz7dJxBNv9Etke1kipk1DqilFatdMI8yCXHnPnX5Im+Hx/6uhY5kxMYNvn+RSUN/U4jqqq/P3dYzhcMhuXjx/xAoqCIGJJmYRcehS1nzkqf5BLjyFGp/aocL5yXgZJMSH8v/eO09zWs0pJW4eL5985RmyE3Vy99EFAOJg333yTNWvWMHnyZF5++eVej33ttde4/PLLWbZsGQ8//DDKIDZIGYGUONatf6VDmEwuPQaSzVsR44sVczOYkB7FS9tPkFvqO6l98FQ1b+7JY+GUROZOHBm6Y31hSXE7ZblM+6dlV+kxxPjRXRLJ53L98vGEh9j4z63ZNDR3+Dzm7T35fHu6hqsuG0ty7Mjaj9QTUsok1OYa1KbBy+p0RpVdyOUne1ztA1gtIjetnUxji4Pn3jzsM4QpKwr/te0oVfVt3LhmCnbbyCyE0YqAcDCTJk3i8ccfZ/Xq1b0eV1RUxDPPPMOrr77K9u3bKSgo4K233jLIyoEh2IIRE8bgKj6i+dhyyTGkpHEIUs8dDUVR4CfrphARYuPx1w5yJK9rAnXvkXL+c+thMhPD2bh8woh/SvYgxmW6N8KWaBvvVx1tKJV5WHqZyMCdWL5lw1QaW5w89sqBLsllRVHZsus0W3fnsXBK0nlVoSSdCSu6NM7DyFW54HJ4x++JUUkRfH/lBI7k1/H0lkNdVjJtHS6e3XqEA6equXbpOManD98Or0YRELEQj+SM2EfHxA8++IBly5YRExMDwFVXXcWWLVtYv3697jYOBkvaVBxfv4na3owQpE0jIqW1AaWuGNu4vru/RIXZ+f/+ZRZP/e8h/vjqQSZlRpMSF0p+WSOnSxsZlxbJ7d+bfl6ExjwIooSUPEH7iaz8FKhyj2HLzoxJieRnV03nz28c5oHnv2TmuDgiQ20czqulsq6NRdOT+cHK88fpA4iRyQghUW7HP/ESzcZ1P0gIWHwk+M/l4ukpKIrKy9tPcu9fvuCC8fGIInxzsprWdhfXLh3H5XPSNbNtJDOsZpSysjJSUs7uLk9JSaGsrP/SG7Gxg5vk4+P7t9Gtfdp8Sr/eSkjjacLStREIbT56iBYgbvIcgs6xx5d98fHhPPHzS3nzs9N88nUxnx8uIyk2lJ9smMaqhaOQDBa07O93qAcNE2ZR8+FBou0dWCLiur0/EBtrDuXQJllInDIT0dr3/rH4+HAmjY3nf3eeYt/RclranIxOieTH66axcFrvO/YD4Tvsi4HYWJk1nba8Q8TFhWnmXEurTmJLyiIhrbuGmy8br7x8IrOnJPPajpMcPl2DoqpMGxvHlUvGMT7D+OZ7w+Fa+8IQB7NhwwZKS0t9vvf5558jScbGMWtqmlGUgSUR4+PDqarqOTHrC9WSCPZQao/upy1hxoDOey7tx74CazCNlniaOtnTl31LZqSwZEZXCZja2sHLlfeHgXyHeiBHuHNXld/uwzpxcZf3Bmpjy6kDSAljqKl3QD+6vX/34tF89+KuCePezh8o32FvDNRGV8xY5MO7qDh5DClm8CsF1dlOe/EJrFOXdbOnNxvDrCI/WtV9xWP09x4I11oUhQE9mBviYN544w1NxklOTu7iqEpLS0lODnxdJkEUsaROQS4+jKqqg34qU1UFV8FBLOnTEEQzyThQxJg0hNAYXIUHuzmYgaA0VaPUFGGff40G1p2/SBnTAXAVHNTEwbiKskFxYcmYOeixTPpHQCT5/WXFihXs2LGD2tpaFEXh9ddfZ9WqVUNtll9Y0qehttZrIhOvVOWhtjViyTRvmMEgCAKWzJm4ig+juvonqeMLV8FBACyZswY91vmMGBKFGJ+Fq+CAJuO5Cg6CPRQpaZwm45n4T0A4mG3btrF48WLef/99nnzySRYvXkxOjnvfyJNPPskrr7wCQHp6OrfccgtXX301y5cvJy0tjbVr1w6l6X4jpU0FwFV0aNBjufIPgCBiSZ8+6LHOdyyZs8DlQC45OuixXAUHECOTEKNGVq+WocCSOROlMheltX5Q46iKjFz4LZb06eZqfwgIiCT/6tWreyxRvvPOO7v8fe2113LttdcaYZamiKHRiLGZuPK/xj7zO4May1VwACl5gmYVaeczUspEsAbhKjgwqBWh6mhFLjuObdoKDa07f7GMmoXjqy24Cg5im3TpgMeRK3JQO5qxjLpAO+NM/CYgVjDnC5Yx89xPZY0D30SmNFai1JWY4TGNECQrlvRpuAoOoqoD37TrjvPLSGZ4TBPE6DSE8LhBh8lcBQdAtGA5E0EwMRbTwRiINWseAM7cLwc8hivffcOZcX7tsGTOQm1rQBlEe2tXwQGEoHCkhDF9H2zSJ+782CzkkiMD1vFTVdW92k+Z2Kuqgol+mA7GQMSIeHfy8vTAHYwz90vEmHTEiPND0sUILBkzQJRw5u4f0OdVV4e7qi9zFkIfm4VN/Mcy6gKQXd7iif6i1BSiNlRgGWU2ShwqzLvBYKxj5qHUFKA0lPf7s0pDBUrlaSxjF+pg2fmLYA/Fkj4dV85e1AFo27nyD4CzHcs487poiZQ8ASE0BmfOFwP6vPPU5yBKWLPmamyZib+YDsZgLJ4w2el9/f6s89TngIB1bN/yMCb9wzLuQtS2BuSS/mvGOU99jhAa4+1nYqINgiBiHbsAuSi7391HVUXGlbPXXT1mFsMMGaaDMRgxLAYpZRLOE5/162lZVWScJz5DSp2MGBajo4XnJ5aMGQj2MJzHP+3X55SmauTibKzjLkQQzNtJayzjF4Gq4DzxWb8+5yo8iNrWgGWCNtJMJgPDvCOGAOvky1CbqpGL/d8T4yo8iNpSi3XKEh0tO38RLDYsEy7Glf8NSkud359zHvsEAOsgSmlNekaKTkFKnojz2Mf9eiBzHv0YITTG3L0/xJgOZgiwjLoAISQKR/aHfn/Gmb3dvGF0xjb5MlBVnEd2+HW86uzAefxTLBkzEcO7i2WaaIN1yhLUpmpcBd/4dbxcW4JcfBjrpEvNzZVDjOlghgBBtGCbthy55Ahy5ek+j3eVn0QuO4Ft+krzhtERMSIBS9ZcHEd2Irc193m88/gq0o0sAAANOUlEQVQnqO1NWGdcYYB15y+WUbMRIhJxHHjbr06XjoNvg8WObbK52h9qTAczRFgnXQb2UDr2b+n1plFVFcf+LQhB4Vgnadcfw8Q3tllrwNlOw943ez1OdbThOPguUvIELKbGla4IooR95ndQqgtw5fe+ipFri3Gd3od18mVmcj8AMB3MECHYgrHPXo9ccqTXm8aVu98tQTJnA4Kl7/4iJoNDik3HMnYh9fveQmmo6PG4jm/eQm1rMJWTDcIy/kLE6FQ69r7SozCpqqp0fP4/YAvBPrP37rgmxmA6mCHEOvkyxJg0Ona/4FPUT2mpo2PPS4ixmVgnXmq8gecp9vlXI0hW2j7+K6ri6va+q+wEzuz3sU64GCkhawgsPP8QRAv2izaiNlXT8cUrPo9xHv4QufQY9rnfM1cvAYLpYIYQQbQQtOSnqI522j54ErX9bNxfaW+i7YMnUF0Ogpb+xNwhbiBiaDTx37kZpfI07Z88j6rI3vfk2mLad/wZISIB+8LrhtDK8w9LyiSs01fiPPYxjkPvd3nPmf81HfteRcqY6Q4/mwQEAaGmfD4jxaQRvOxm2j78My3/9wC2qZejqirOIx+itjcRfPntSFEpfQ9koilhky+ivigfx1dbaG0oxzp+EWpLLY4jOxGsQYSs+JmpbzUE2Odd6V7F7P0nctkJpIwZKJWncZ7cjRg/muDLbtSszbLJ4DEdTABgyZxFyJp7ad/zMh37XgVw3yyX32aKJw4h9gvWIkYk0LHvNTr2vAQIWEbNwr7wOrMseYgQRAtBS2/G8e27OA6+61VLtk5e6g6NmU4/oDAdTIAgJY4l9LsPeiUxxOCIIbbIBMA6dgGWMfNRW2oRbCHmBBYACKKEfdYabDOuQG2uRQiNQpCsQ22WiQ9MBxNgmI4l8BAEASEsdqjNMDkHQZQQIuKH2gyTXjAzxyYmJiYmumA6GBMTExMTXTAdjImJiYmJLpgOxsTExMREF0wHY2JiYmKiC6aDMTExMTHRhfOyTFkUB7fTd7Cf15tAtw9MG7Ug0O0D00atGGobB3p+QfWnwYKJiYmJiUk/MUNkJiYmJia6YDoYExMTExNdMB2MiYmJiYkumA7GxMTExEQXTAdjYmJiYqILpoMxMTExMdEF08GYmJiYmOiC6WBMTExMTHTBdDAmJiYmJrpgOhg/ycvL45prrmHFihVcc8015OfnD7VJXairq+PGG29kxYoVrFmzhttuu43a2tqhNssnzzzzDBMmTODkyZNDbUo3Ojo62Lx5M8uXL2fNmjX8+te/HmqTuvHxxx+zfv161q1bx5o1a9i+ffuQ2vPoo4+yZMmSbtc0kO4ZXzYG2j3T0/foIZDvmx5RTfzi+uuvV7du3aqqqqpu3bpVvf7664fYoq7U1dWpe/fu9f79+9//Xv3lL385hBb55vDhw+qmTZvUSy+9VD1x4sRQm9ONRx55RP3tb3+rKoqiqqqqVlVVDbFFXVEURZ0zZ473uzt27Jg6c+ZMVZblIbNp//79amlpqXrZZZd1uaaBdM/4sjHQ7pmevkdVDfz7pifMFYwf1NTUcPToUVavXg3A6tWrOXr0aECtEKKiopg/f77375kzZ1JaWjqEFnXH4XDw8MMPs3nzZgQh8AQGW1pa2Lp1K3feeafXvri4uCG2qjuiKNLU1ARAU1MTCQkJiOLQ3cpz5swhOTm5y2uBds/4sjHQ7hlfNkLg3ze9cV6qKfeXsrIyEhMTkSQJAEmSSEhIoKysjJiYmCG2rjuKovDKK6+wZMmSoTalC08++SRr164lPT19qE3xSVFREVFRUTzzzDPs27eP0NBQ7rzzTubMmTPUpnkRBIEnnniCW265hZCQEFpaWnjuueeG2qxumPeMdgT6fdMb5gpmBPLII48QEhLCxo0bh9oULwcOHCA7O5vrrrtuqE3pEZfLRVFREZMnT2bLli3cfffd3H777TQ3Nw+1aV5cLhfPPfcc//mf/8nHH3/Ms88+y1133UVLS8tQmzasCcR7BobHfdMbpoPxg+TkZCoqKpBlGQBZlqmsrPS5nB1qHn30UQoKCnjiiSeGNGxyLvv37yc3N5elS5eyZMkSysvL2bRpE7t37x5q07ykpKRgsVi8YZ0ZM2YQHR1NXl7eEFt2lmPHjlFZWcns2bMBmD17NsHBwZw+fXqILeuKec9ow3C4b3ojsL7NACU2NpZJkyaxbds2ALZt28akSZMCbqn/+OOPc/jwYf785z9js9mG2pwu3HTTTezevZuPPvqIjz76iKSkJJ5//nkWLVo01KZ5iYmJYf78+ezZswdwV0HV1NSQmZk5xJadJSkpifLycnJzcwE4ffo01dXVZGRkDLFlXTHvGW0YDvdNb5gNx/zk9OnT3HvvvTQ2NhIREcGjjz5KVlbWUJvl5dSpU6xevZpRo0YRFBQEQFpaGn/+85+H2DLfLFmyhL/85S+MHz9+qE3pQlFREffddx/19fVYLBZ+9rOfcckllwy1WV146623+Nvf/uZN+N5xxx0sW7ZsyOz5zW9+w/bt26muriY6OpqoqCjeeeedgLpnfNn4xBNPBNQ909P32JlAvW96wnQwJiYmJia6YIbITExMTEx0wXQwJiYmJia6YDoYExMTExNdMB2MiYmJiYkumA7GxMTExEQXTAdjYjIAHnjgAW8561dffcWKFSuG2KKe+ec//8lvf/vbPo+77bbb2LVrlwEWmZwvmGXKJuc1mzZtYvr06dx5551dXt+xYwebN2/m008/xWIZHpJ9Tz/9NAUFBfzHf/yH9zWHw8Hll1/Oa6+9RmJiYq+fP3ToEA8++CBbtmzR21ST8wRzBWNyXrNhwwbefPNNzn3Oeuutt1izZo1P5+KRPxkO7Ny5k6ysrD6dC8D06dNpbm4mOzvbAMtMzgdMB2NyXrNs2TIaGhr46quvvK81NDR4m3oB3HvvvWzevJkbb7yRmTNnsm/fPu69914ef/xxAPbt28fixYu9n6+oqOD2229nwYIFLFmyhBdffNH73tNPP80dd9zB3XffzaxZs1izZg15eXk899xzLFy4kEsuuaSLzlRTUxP33XcfixYt4uKLL+bxxx/36eB27drFc889x3vvvcesWbNYu3at9/W5c+d6j+vo6ODuu+9m/vz5zJkzh+9973tUV1d73583bx6ffvrpYL9WExPAdDAm5zlBQUGsWrWKrVu3el977733yMrKYuLEid7Xtm3bxk9/+lO++eYbr9CkLxRF4eabb2bChAns2rWLF154gRdeeIHPPvvMe8zHH3/MunXr2L9/P5MmTWLTpk0oisKuXbu49dZbeeCBB7zH3nPPPVgsFrZv387WrVvZs2cPr7/+erfzLl68mJ/85CesWrWKAwcO8NZbbwFw8uRJRo8e7T3ujTfeoLm5mU8++YR9+/bx0EMPeWVSAMaMGcPx48f7+S2amPjGdDAm5z3r16/n/fffp729HYCtW7eyYcOGLscsXbqU2bNnI4oidru9x7Gys7Opra3ltttuw2azkZ6eztVXX827777rPWbOnDlcfPHFWCwWVq5cSV1dHTfddBNWq5UrrriCkpISGhsbqa6uZteuXdx3332EhIQQGxvLDTfc0E2fqjeampoIDQ31/m2xWKivr6egoABJkpg6dSphYWHe90NDQ2lsbPR7fBOT3hge2UsTEx2ZM2cOMTEx7Ny5k+nTp3P48GGeeeaZLsf4KzNfUlJCZWVllyZlsix3+Ts2Ntb7/0FBQURHR3sbc3lWE62trVRWVuJyuboo5yqK0i/J+4iIiC69YtatW0d5eTk///nPaWxsZO3atdx1111YrVbA3dUzIiLC7/FNTHrDdDAmJrgn3q1bt5KXl8dFF1004FbJycnJpKWlsX379kHblJSUhM1mY+/evX5VsvlqpzthwgTy8/O9f1utVm677TZuu+02iouLuemmmxg9ejRXXXUV4FYN7xwaNDEZDGaIzMQEd5jsiy++4LXXXvMm9wfC9OnTCQsL469//Svt7e3IsszJkyc5dOhQv8dKSEjgoosu4ve//z3Nzc0oikJhYSFffvmlz+NjY2MpKSlBURTva5dccgn79+/3/r13715OnDiBLMuEhYVhsVi8qydwN7jqXLBgYjIYTAdjYoK7D8isWbNoa2tj6dKlAx5HkiSeffZZjh8/ztKlS1mwYAH333//gNsuP/bYYzidTq644grmzp3LHXfcQVVVlc9jV65cCcD8+fO9OaTLLruM3NxcKioqAKiuruaOO+5g9uzZXHHFFcybN89bcXbo0CGCg4OZPn36gGw1MTkXc6Olickg+eKLL7j//vvZuXPnUJvik1dffZWcnBx+9atf9Xrc7bffzpVXXhlwDdZMhi+mgzExGSQvvvgiO3fu5IUXXhhqU0xMAgozyW9iMgh+85vf8NFHH/Hoo48OtSkmJgGHuYIxMTExMdEFM8lvYmJiYqILpoMxMTExMdEF08GYmJiYmOiC6WBMTExMTHTBdDAmJiYmJrpgOhgTExMTE134/wGhe6JlwjxFeAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "t=np.arange(0,15,0.05)\n", "f=0.2\n", "omega=2*np.pi*f\n", "x=np.cos(omega*t)\n", "v=-omega*np.sin(omega*t)\n", "\n", "plt.plot(t,x,label='Polozaj x(t)')\n", "plt.plot(t,v,label='Brzina v(t)')\n", "plt.legend()\n", "plt.title('Zavisnost elongacije i brzine od vremena')\n", "plt.xlabel('Vrijeme t(s)')\n", "plt.ylabel('Vrijednost x(m), v(m/s)')\n", "\n", "plt.show()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Zadatak 5\n", "Napisati skriptu koja učitava vrijednost temperature u stepenima celzijusa, a zatim nudi korisniku da pritiskom na F (Farenheit) ili K (Kelvin) prikaže vrijednost temperature u odabranoj skali $\\left(F=\\frac{9}{5}C+32\\right)$,$\\left(K=C+273.15\\right)$." ] }, { "cell_type": "code", "execution_count": 50, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Program za konverziju °C u Kelvin ili Fahrenheit\n", "Unesi temperaturu u Celzijusima76.3\n", "Odaberi K za Kelvin, F za Fahrenheit k\n", "Temperatura od 76.3 °C odgovara 349.45 Kelvina\n" ] } ], "source": [ "import numpy as np\n", "\n", "print('Program za konverziju °C u Kelvin ili Fahrenheit')\n", "C=float(input('Unesi temperaturu u Celzijusima '))\n", "\n", "izbor=input('Odaberi K za Kelvin, F za Fahrenheit ')\n", "izbor=izbor.upper() #ako korisnik unese malo slovo\n", "\n", "K=C+273.15; F=9/5*C+32\n", "\n", "if izbor=='K':\n", " print('Temperatura od ',C,' °C odgovara ',K,' Kelvina')\n", "elif izbor == 'F':\n", " print('Temperatura od ',C,' °C odgovara ',F,' Fahrenheita')\n", "else:\n", " print('Nepravilan izbor znaka za konverziju')\n", " \n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Zadatak 2\n", "U nekom fizičkom eksperimentu je mjerena zavisnost fizičke veličine Y od druge fizičke veličine X. Rezultati mjerenja su snimljeni u datoteci rezultati_mjerenja.csv. Učitati podatke i izvršiti fit metodom najmanjih kvadrata (linearna metoda). Rješenje prikazati na grafikonu (podaci, dobijena prava i navesti vrijednosti koeficijenata). \n" ] }, { "cell_type": "code", "execution_count": 56, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "\n", "data=pd.read_csv('rezultati_mjerenja.csv', header=None)" ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
01
05.110117.5920
15.52779.1302
28.518613.6620
37.003211.8540
45.85986.8233
\n", "
" ], "text/plain": [ " 0 1\n", "0 5.1101 17.5920\n", "1 5.5277 9.1302\n", "2 8.5186 13.6620\n", "3 7.0032 11.8540\n", "4 5.8598 6.8233" ] }, "execution_count": 57, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data.head()" ] }, { "cell_type": "code", "execution_count": 91, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEcCAYAAAA2g5hwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3XtUVPX6P/D3zOCgBjaCgKOSGieNMMvAzMww1BLDNE8mmVpej5VaGRidPGCG1WhLy9Q49k3TltnJNAqUvOQpzbCTl0rB289l6pGbAioooMzs3x+cmRjYw56B2Xtu79darRUze/Z+GLafZ3/uKkEQBBARETVB7eoAiIjI/TFZEBGRJCYLIiKSxGRBRESSmCyIiEgSkwUREUlisiBZbd68GU899ZSrw0CfPn1w7tw5V4chKSUlBUuXLrXr2NTUVKxYsULmiICff/4ZDz74YIvP8+ijj+Lnn392QkTkCn6uDoA8X1xcHC5evAiNRmN57fHHH0dqaqoLo7J26NAhl1w3Li4O6enpuP/++51+7gULFjj9nHLasmWLq0OgFmCyIKfIyMiQpUBsqdraWvj58TYnaik2Q5GiDAYDnnrqKVRUVODMmTMYP348oqOj0a9fP7z00kuW43r27Il169Zh8ODB6NevHwwGA0wmk+X9L7/8EvHx8ejbty+mTJmC8+fPW312/fr1ePjhh/Hwww9bXjtz5gyAuqae+fPnY+rUqejTpw8SExNx4cIFLFy4EH379sWwYcOQn59vOV9xcTFmzZqF++67D3FxcVi3bp3lvQ8++AAvvvgi5s6diz59+uDRRx/F4cOHAQDJyckoKCjAjBkz0KdPH3z00UcAgNmzZ2PAgAGIjo7G008/jZMnT4p+V+bmn9WrV6N///544IEHsGnTJsv7TTVZbd68GYmJiXjzzTcRHR2NYcOGITc31+p3mjFjBu69914MHToUX3zxheW96upqpKSkoG/fvhg+fLjl9zFbtWoVhgwZgj59+mD48OHYsWOH1ftffPEF4uPjLe/n5eUBqKtl/fTTT6LxkvtjsiBFmEwmzJs3DydOnMDq1asRGBiI999/HwMGDMAvv/yC3bt3Y/z48Vaf2bFjBzZt2oSvvvoKu3btshSUO3fuxD//+U8sX74cubm5iI6OxiuvvGL12Z07d+KLL77A1q1bRePJycnBSy+9hH379kGr1WLs2LGIiorCvn378Mgjj+Dtt9+2xP3cc8+hZ8+e2L17N9auXYu1a9diz549lnPt2rULjz76KPbv34+4uDi8+eabAIDFixejU6dOyMjIwKFDhzBt2jQAwIMPPoht27YhNzcXd9xxB5KSkmx+bxcvXkRFRQV2796NhQsXYsGCBbh8+bJd3/nvv/+O8PBw7Nu3D7Nnz8bMmTNx6dIlAMArr7yCjh07Ys+ePVi2bBmWLFliSSbLly/H2bNnsWPHDnz88cfIzMy0Om94eDjWr1+PAwcOYObMmUhOTkZJSYnle/3ggw9gMBhw8OBBfPjhh9DpdHbFS+6NyYKc4oUXXkBMTIzlv/pPqrW1tZgzZw4uX76MDz/8EG3atAEA+Pn5oaCgACUlJfD390dMTIzVOadNmwadTodOnTph4sSJyM7OBgB8/vnnmD59OiIiIuDn54cZM2bg6NGjVrWL6dOnQ6fToXXr1qLxDh06FL169YK/vz+GDh0Kf39/jBo1ChqNBsOHD8fRo0cBAIcPH0ZZWRlmzpwJrVaL8PBwPPnkk1ZJKDo6GrGxsdBoNBg5ciSOHTvW5Hf1xBNPICAgAFqtFrNmzcKxY8dQUVEheqyfnx9eeOEFtGrVCrGxsWjbti1Onz7d5PnNgoKC8Mwzz6BVq1YYPnw4unfvju+//x6FhYU4cOAAkpKS4O/vj8jISIwZMwZff/01gLoCf8aMGdDpdNDr9ZgwYYLVeePj4xEWFga1Wo3hw4eja9eu+P333wHU1fimTp2K3r17Q6VSoWvXrujcubNd8ZJ7Y2MuOcWKFSts9lmcPXsWx44dw8aNG6HVai2vJycn4/3338cTTzyBm2++GZMmTcITTzxheV+v11v+v3Pnzpan14KCArz11lswGAyW9wVBQHFxsaVgqv9ZMcHBwZb/b926NTp06GD187Vr1wAA58+fR0lJiVUiMxqNVj83/GxNTY3NvhKj0YilS5fi22+/RVlZGdTquue18vJyBAYGNjpep9NZnadNmzaW2KSEhYVBpVJZfu7UqRNKSkpQUlKCm2++GQEBAVbvHTlyBABQUlJi9f116tTJ6ryZmZlYs2aNJTlfu3YN5eXlAIDCwkLccsstdsVHnoXJgmR366234umnn8a0adOwdu1a3HrrrQCAkJAQpKenAwD279+PSZMmoW/fvujatSuAuoLntttuA1CXIEJDQwHUJYIZM2bgscces3nN+oVkS+j1enTp0gXbt293yvmysrLw3XffYc2aNejSpQsqKirQt29fyLH4c3FxMQRBsHwXhYWFiIuLQ2hoKC5fvozKykpLwigsLERYWBiAur9L/e++sLDQcs7z589j3rx5+OSTT9CnTx9LbcpMr9fj7NmzTv9dyPXYDEWKSEhIwJw5czBp0iRLYZKTk4OioiIAwM033wyVSmV50gaAjz/+GJcvX0ZhYSHWrVuH4cOHAwASExOxatUqS8dwRUUFcnJyZIm7d+/eCAgIwKpVq1BdXQ2j0YgTJ05Yml2kdOjQwWp+x9WrV6HVatG+fXtUVVVhyZIlssQNAGVlZVi3bh1u3LiBnJwcnDp1CrGxsdDr9ejTpw+WLFmCmpoaHDt2DF9++SVGjBgBoK6ZadWqVbh8+TKKiorw6aefWs5ZVVUFlUqFoKAgAMCmTZusOuifeOIJrF69GkeOHIEgCDhz5oxV8yB5LtYsyClmzJhhNc/i/vvvbzRh7PHHH8eNGzfwzDPP4NNPP8Xhw4fx1ltvobKyEsHBwXj99dcRHh5uOX7w4MEYPXo0Kisr8fjjj1uaqIYOHYqrV69izpw5OH/+PAIDA3H//fcjPj7e6b+XRqPBhx9+CIPBgMGDB+P69evo3r271citpkyfPh3p6elYvHgxnnvuOSQmJuLHH3/EwIEDodPp8OKLL2LDhg1OjxuoS3RnzpzBfffdhw4dOmDZsmVo3749AGDJkiVIS0vDwIED0a5dO8yaNQsDBgwAAMycORNpaWkYPHgwQkNDMXr0aMsIsL/85S+YPHkyEhMToVKpMGrUKNxzzz2Wa8bHx+PSpUt45ZVXUFJSgs6dO2PRokXst/ACKm5+RO6oZ8+e2L59u6VJihyzefNmbNy4UbZERL6HzVBERCSJyYKIiCSxGYqIiCSxZkFERJKYLIiISJIiQ2fLy8sxd+5cnD17FlqtFl27dsWCBQsQFBSEnj17okePHpbx9YsWLULPnj2VCIuIiOykSJ/FpUuXcPz4cfTr1w9A3cqjly9fxltvvYWePXvi4MGDuOmmm5p9/vLyqzCZHP81goMDUFpa2ezrKo3xysvT4gU8L2bGKy9741WrVWjf3rEyV5GahU6nsyQKALj77rudOv7bZBKalSzMn/UkjFdenhYv4HkxM155yRWv4jO4TSYTNmzYgLi4OMtrEyZMgNFoxIMPPohZs2ZZLTZHRESup/jQ2TfeeAPFxcVYvnw51Go1CgsLodfrUVlZieTkZPTo0QMvv/yykiEREZEERWsWBoMBZ86cQUZGhqVD27wUckBAAMaMGYM1a9Y4fN7S0spmVb1CQgJx4YL4PgLuiPHKy9PiBTwvZsYrL3vjVatVCA4OkDzO6jPNDcpRS5cuxZEjR7BixQpLM9Ply5dRXV0NoG6DnG3btiEyMlKpkIiIyE6K1CxOnjyJjIwMdOvWDYmJiQCALl26YOrUqUhNTYVKpUJtbS369OmDF198UYmQiIg8Um5eETb/cAqlV2oQ3M4fo2Mj0D+qo+zXVSRZ3HbbbTh+/Ljoe1lZWUqEQETk8XLzirA25xiu15oAAKVXarA2p24bX7kTBmdwExF5iM0/nLIkCrPrtSZs/uGU7NdmsiAi8hClV2ocet2ZmCyIiDxEcDt/h153JiYLIiIPMTo2Alo/62Jb66fG6NgI2a/NPbiJiDyEuRPba0dDERGRc/SP6qhIcmiIzVBERCSJyYKIiCQxWRARkST2WRCRV3DVMhi+gsmCiDyeK5fB8BVshiIij+fKZTB8BZMFEXk8Vy6D4SuYLIjI47lyGQxfwWRBRB7Plctg+Ap2cBORx1NiGQxfH23FZEFEXkHOZTA42orNUEREkjjaismCiEgSR1sxWRARSeJoKyYLIiJJHG3FDm4iIknmTuzPdhzH1WojAEDbyreetX3rtyUiaoEbtYLl/yurarE25xhy84pcGJFymCyIiOzg6yOimCyIiOzg6yOimCyIiOzg6yOimCyIiOzg6yOiOBqKiMgOSqw/5c4USRbl5eWYO3cuzp49C61Wi65du2LBggUICgrCr7/+itTUVNTU1KBz585YvHgxgoODlQiLiMghcq4/5e4UaYZSqVSYOnUqtm3bhqysLISHh+Pdd9+FIAhITk5Gamoqtm3bhpiYGLz77rtKhERERA5QJFnodDr069fP8vPdd9+NgoICHD58GP7+/oiJiQEAJCYm4ttvv1UiJCIicoDifRYmkwkbNmxAXFwcCgsL0alTJ8t7QUFBMJlMuHTpEnQ6nd3nDA4OaHY8ISGBzf6sKzBeeXlavIDnxcx45SVXvIonizfffBNt27bF+PHjsWPHDqecs7S0EiaTIH1gAyEhgbhwocIpMSiB8crL0+IFPC9mxisve+NVq1UOP2QrmiwMBgPOnDmDjIwMqNVq6PV6FBQUWN4vKyuDSqVyqFZBRETyU2yexdKlS3HkyBGsWLECWq0WANCrVy9UV1dj//79AIDPP/8c8fHxSoVERER2UqRmcfLkSWRkZKBbt25ITEwEAHTp0gUrVqzAokWLkJaWZjV0loiI3IsiyeK2227D8ePHRd+75557kJWVpUQYRETUTFzug4iIJDFZEBGRJCYLIiKSxGRBRESSmCyIiEgSkwUREUlisiAiIknc/IiISEG5eUUeuYESkwURkUJy84qwNucYrteaAAClV2qwNucYALh9wvDpZPH9gXP4JDvP4zI8EXmmzT+csiQKs+u1Jmz+4ZTblz0+myxy84qw7tvjqLlhBOBZGZ6IPFPplRqHXncnPtvBvfmHU5ZEYWbO8EREcghu5+/Q6+7EZ5OFJ2d4IvJMo2MjoPWzLna1fmqMjo1wUUT289lmqOB2/qKJwRMyPBGJc/eRRuZY3DlGW3w2WYyOjbDqswA8J8MTUWOeMtKof1RHt4rHXj7bDNU/qiNmjrnLUpMIbuePZ+Jv98g/IhE1PdKIWs5naxYAMCg6HFG3cL9vIm/Afkh5+WzNgoi8iyePNPIETBZE5BU8eaSRJ/DpZigi8h6ePNLIEzBZEJHX8NSRRp6AzVBERCSJyYKIiCQxWRARkST2WRApxN2XoiBqCpMFkQI8ZSkKIlvYDEWkAC5FQZ6OyYJIAVyKgjydYs1QBoMB27Ztw/nz55GVlYUePXoAAOLi4qDVauHvXzclPykpCQMHDlQqLCJFcEl88nSKJYvBgwdj4sSJePrppxu9t2zZMkvyIPJGo2MjrPosAC5FQZ5FsWQRExOj1KWIZNWcUU1cioI8nUoQBEHJC8bFxSEjI8OqGSogIACCICA6Ohpz5sxBu3btlAyJyG7fHziH5Rt/s9o0y7+VBjPH3IVB0eEujIxIXi4fOrt+/Xro9Xpcv34dCxcuxIIFC/Duu+86dI7S0kqYTI7nvJCQQFy4UOHw51yF8crLnng/yc6zShQAUHPDiE+y81yyN4o3fsfuxFvjVatVCA4OcOjcLh8NpdfrAQBarRbjxo3DwYMHXRwRkW0c1US+yqXJ4tq1a6ioqMuCgiBg69atiIyMdGVIRE3iBjvkqxRrhkpPT8f27dtx8eJFTJo0CTqdDhkZGZg1axaMRiNMJhMiIiKQlpamVEhEDuOoJvJViiWLefPmYd68eY1ez8zMVCoEohbjqCbyVS7v4CbyNNxgh3wRk4UP4aqnRNRcTBY+gqueOo7JlehPLh86S8rgqqeOMSdX85BYc3LNzStycWRErsFk4SM4P8AxTK5E1pgsfATnBziGyZXIGpOFjxgdGwGtn/Wfm/MDbGNyJbLGZOEj+kd1xDPxt1sKu+B2/ngm/nZ22NrA5EpkjaOhfAjnB9iPk++IrDFZENnA5Er0JzZDERGRJCYLIiKSxGYoahJnMRMRwGThcxwp/LlEiPtp+Pd7NiHKJTv0ke9hsnAxuZ/c65//ptYa1NwwodZYtwWtVOHf1CxmJgvliSXv5Rt/w8RhPfn3INmxz8KF5F5/qOH5r1YbLYnCrKklLDiL2b2IJe+aG0YuQUKKYLJwIbnXHxI7vxhbhT9nMbsXJm9yJSYLF5L7H7+957FV+HMWs3th8iZXYrJwIbn/8dtznqYKfy4R4l7Ekrd/Kw2TNymCHdwuNDo2wqrDEnDuk7vY+TUqoE1rP1RW1drVoc5ZzO5DbAkSjoYipTBZuJDc6w9xfSPv0zB5h4QE4sKFChdGRL6CycLF5H5yZ82AiJyBfRZERCTJ7mTx9ttv4+jRo3LGQkREbsruZqja2lpMmTIFQUFBeOyxx/DYY4+hY0c2bxAR+QK7axb/+Mc/sGfPHrzyyis4duwY4uPj8eyzzyIzMxNXr16VM0YiInIxh/osNBoNHnroISxZsgRffPEFysrKkJKSggceeACvv/46iouL5YqTiIhcyKFkUVlZiY0bN2LChAkYP3487rrrLqxfvx5bt25F27ZtMXXqVNHPGQwGxMXFoWfPnjhx4oTl9dOnT2Ps2LF45JFHMHbsWPzxxx8t+mWInCk3rwjJK/di8ju7kLxyr9PW7CLyRHb3WcyePRt79uxB37598dRTT2HIkCHQarWW91977TVER0eLfnbw4MGYOHEinn76aavX09LSMG7cOIwcORJff/01UlNTsW7dumb+KkTOw+XZiazZXbO46667sH37dqxatQrDhw+3ShQAoFar8dNPP4l+NiYmBnq93uq10tJS5OfnIyEhAQCQkJCA/Px8lJWVOfo7EDmd3Is8Enkau2sWU6ZMkTymTZs2dl+4sLAQYWFh0Gg0AOr6Q0JDQ1FYWIigoCC7z0MkB67wSmTNK2ZwBwcHNPuzISGBToxEfoxXXuZ4Q9q3wYXyqsbvt2/jdr+Tu8UjhfHKS654XZYs9Ho9iouLYTQaodFoYDQaUVJS0qi5yh6lpZUwmQTpAxvwtHV1GK+86sc76oHuoos8jnqgu1v9Tp78HXsCb41XrVY5/JDtsuU+goODERkZiezsbABAdnY2IiMj2QRFboHLsxNZU6RmkZ6eju3bt+PixYuYNGkSdDodtmzZgvnz5yMlJQUrV65Eu3btYDAYlAjHLnLvjU3uj4swEv1JJQiC4+03bsbZzVANh00CdU0Qrn6y9NYqsbtobryufLDwle/YVbw1Xo9qhnJnHDZJ9jI/WJhHSZnnY3ACH3kbJgsRHDZJ9uKDBfkKJgsRcu+NTd6DDxbkK5gsRIyOjYDWz/qrcebe2OQ9+GBBvoLJQgSHTZK9+GBBvsIrZnDLgcMmyR7me4TDrMnbMVkQtRAfLMgXMFmQ1+BESiL5MFmQ7JQoxLn/BJG82MFNslJq0hrnOxDJizULarGmag5NFeLOfOLnfAciebFmQS0iVXNQqhDnfAcieTFZ+IDcvCIkr9yLye/sQvLKvU5tApJq/lGqEOd8ByJ5sRnKQZ424kbujl+pmsPo2AjRFXydXYhzvgORvJgsHOCJI27k7jMIbucvmjDMNYfmFuLNScqc70AkHyYLByjVWetMcvcZ2FNzcLQQ98SkTOTtmCwc4MoRN7l5Rcj8MRcXyqssT9qA9BO71JO/1DWlzi9Vc6h/DrUKMAmQrCl4YlIm8nZMFg4wF3Zir4sVrIBz2tDFnrTXbD0KwSTA+L94bD19N7fPwJGne1s1h+8PnLM6h0kiVjMOgyVyPxwN5QBbO7eaBDQaPrpm61Gszs53ymQ0sSftWuOficJMbBJac1fQdcYkt3U5Rxudw55zcRgskfthzcIBtpp01CqIFuYNNbcpxZEnarFjm9Px64yn+4vlVc26hlIjqIjIfqxZOMDWWH5bNQ4xzWlKceSJ2llP3854uu/Qvk2zrsH9RIjcD5OFA2wVYo4UoOb+DUeIJSk/jQoalfVxznz6dsYkt4nxkY3OYe+5+kd1xOLnB2B1ShwWPz+AiYLIxdgM5SBbTToNm038NCqrDmgzc/+G+Vz2XhMAMn887fBoqOZyxiS3QdHhuFJR7fBoKCJyP0wWTmCrYAWAj7PzGzVTNafvon9URzw26DZcuFAhem0xLZ1t7oxJbpwoR+QdfDJZmAvRsis1CHLSU66tQvGjrHzR4+UeBsqJbUTkTD6XLJQoRMUmojUk9zBQTmwjImfyuWQhdyHaMBmJJYqGnbtyLE4ox8Q2d19E0d3jI/JkPpcs5J4dvGHnCZsT0QDgptYajBva02o5DDlqOi1Z5kOMuzdruXt8RJ7OLZJFXFwctFot/P3rCrKkpCQMHDhQlmvZKkQD2rT8q8jNK0JlVW2Tx7TW+lkVXnLVdOyd2Gbv07i7N2u5e3xEns4tkgUALFu2DD169JD9OqNjI7Bm69FGM6yrqmuRm1fUooLFnqUwGiYqR2o63x84h0+y8+xqZrFn6KsjT+Puvl6Tu8dH5OncJlkopX9UR3y24zhqjUar140CWvwUak/B1LAZyFZNBwCSV+61FPC5eUVY9+1x1NwwWq4l1cwiNWzVkadxZzdrOZu7x0fk6dxmBndSUhJGjBiB+fPn48qVK7Je62q1UfT1lj6FShVMYs1AYjOl68djXnxw8w+nLInCzNGF/cTOb+/r7r5tqbvHR+TpVIIgOLCykTwKCwuh1+tx/fp1LFy4EFevXsW7774r2/Ump2/HBZFF7kLat8HqeQ9bvfb9gXNYl3MUF8ur0KF9G0yMj8Sg6HDR835/4Bze/9ch0UUEA9u2wvRRd4p+1nwNsZjMcV0sr4LYH0oFYM64e+yOsT5Hvof6cTp6HaW4e3xEnswtkkV9x48fx3PPPYddu3bZ/ZnS0kqYHFjNr2FbPVD3FNpwsTp7j2t47g07T1g6uhuOfmrK5Hds/85Ndcxfv2FyKMb6sTr6+zkiJCSw0Yxzd+Zp8QKeFzPjlZe98arVKgQHBzh0bpf3WVy7dg1GoxGBgYEQBAFbt25FZGSkrNe0tdaS2O5uDUmNsGnJ8hZN9V+Iva5RAYIg2N3vIDby6Zn42zk3gYgkuTxZlJaWYtasWTAajTCZTIiIiEBaWprs17W11pLY03ZDco2wERvu2hSVWtVk/0v90V22Rj49E387Fj8/wDm/ABF5LZcni/DwcGRmZro6DAuxEUINyTXCpuFwVym1RsHmciKA9eq2LZmHwJnRROTyZOEquXlFyPwxt1EzlFQh7awRNrYKYHMhbG8NwyTUxSR2bP1k0Nx5CJwZTUSAGw2dVZK5ADSPBKo/RLWpWoOzdmwzX9/W/tz21G4axmSL+RrN3fnOGXtxE5Hn88lk0VQBaGu8/rQRdzi0Y1tuXhGSV+7F5Hd2IXnlXqvd8aQKYEf6RMw1Eqlk0Nx5CJwZTUSAjyaLpgpAc0Gu/t+Wpc2pTUjVHKQKYEf6RMxxSSWD5u5r7Yy9uInI8/lkn4U9Q1RNQt3Q1JobRnyUlW+pddiTNKQ6k5tamiI3r6jRTO2mfg8ze9aCas6wXnsXJCQi7+aTycLeIapGAZbJdaVXavBRVj4+yspHQBs/PDWkR6MJfFKjmMzv2SqAe0cEi8al9VPBJMBqZrhYgS3HFqbO2IubiDyfTyYLsUl5jrTBV1bVYs3Wo5Zz2TM3A/izJmCrALbVsR3YVovRsRE2JxHKjftoE5FPJgug8aS85JV7HUoYtUbB0qxkz+ilhjUBsQK4qf26bU0iJCJSgs8mi4bzLHpHBGPv4SK7h6wCsOrAboq9NQEus01E7sonk4XYRLO9h4sw4M6O+P1Uqd01jPoji2wV8o4spcHOZCJyVz45dNbWaKXfT5Vi8fMDsDolTvJp3k+jshTitvak6B0R7FBczR3eSkQkN5+sWdgz0aypEVMqFTCwt95SiPeP6oj/999L+PehAqvj9h4uwl+66Bwq7NmZTETuyCdrFvZMNGv4lF+fINQlgvqzsn8/VdroOC6LQUTewieThb3NRv2jOmLx8wNEE0bDRMBlMYjIm/lkM5StZqN/HyrAvw8VNBq9ZE8i4EgmIvJmPlmzAMSbjcxKr9Rgzdajlmamm1prbB5rPqa5C/UREXkCn6xZANLNQ7VGARt2ngAA1NywPfdidXbdRDoui0FE3sxnk4U9S3xUVtVi8w+nrNZkasgowDKTmyOZiMhb+WyyGB0bgXXfHpdc4dWeDur6x3ALUiLyRj6bLPpHdUS7wNb4JDvPZkK4qbUGrbV+kglDrfqz74JbkBKRN/LZZAEAg6LDEXWLDrl5RVidnY/6rU0aFTBuaE8A0vthm4S6Y1r5qZrcx4KIyFP5dLIwa9g5fVNrDWqNf64Cq/VTIaCNHyqraqFS1U3Ka+h6rQnXa8XPb09TFpuviMid+ezQ2YbME/CmjbgDNTdMVn0Z12sFVFXXYtqIO/Dxq3EOn1tqroXUNqxERK7GZNGArdFPRgH4ODsfuXlFNudd+LdSNWuuRVPbsBIRuQOfbob6/sA5Swe3uemnqSYjc9+ESqUSfb+VnwZPDemBz3Ycx9XqupqJtpV0PuZSIUTk7nw2WeTmFVkNnS29UoP/s7FTXX1NdXSb9+u+UStYvSY1IopLhRCRu/PZZqjNP5xqNMfC9tQ7+32cne9wkxKXCiEid+ezNQu5mnhMNjJOU9fjUiFE5O7cIlmcPn0aKSkpuHTpEnQ6HQwGA7p16ybrNe1Z7sPZ12uKHEuFcDguETmLWzRDpaWlYdy4cdi2bRvGjRuH1NRU2a/pjCYee/sUXNGkxOG4RORMLk8WpaWlyM/PR0JCAgAgISEB+fmvCEIwAAALdklEQVT5KCsrk/W6LX3CNicAWwlD/b8BU67aR5vDcYnImVzeDFVYWIiwsDBoNHVzFzQaDUJDQ1FYWIigoCBZrx3Svg0ulFc1el2tEu97ML/esEmn4XIgWj+1SxJEfRyOS0TO5PJk4QzBwQHN+tzE+Egs3/ib1ago/1YaDI7pgu/2/7fR6zPH3IVB0eFW53hsUCDaBbbGupyjuFhehQ7t22BifGSj45wlJCTQvuNsJMKQ9m3sPoczKHktZ/C0eAHPi5nxykuueF2eLPR6PYqLi2E0GqHRaGA0GlFSUgK9Xm/3OUpLK2GyNQypCYOiw3Glolq0E7hzh5savR51iw4XLlQ0Ok/ULToY/tbf6jWx41oqJCTQ7vOOeqC7aI1n1APdZYlNjCPxugNPixfwvJgZr7zsjVetVjn8kO3yZBEcHIzIyEhkZ2dj5MiRyM7ORmRkpOxNUFI8fSMjDsclImdyebIAgPnz5yMlJQUrV65Eu3btYDAYFLnu9wfOefX+E56e8IjIfbhFsoiIiMDGjRsVv+66nKPcf4KIyA4uHzrrShdFOoABjhgiImrIp5NFh/ZtRF8PaOMWFS4iIrfh06XixPhIvP+vQ432r6iqrkVuXlGzm6K4zAYReRufrlkMig6Hv8h+E0YBzZ7pzGU2iMgb+XSyAGDZpKih5vZbcJkNIvJGPp8sbK3t1NyNh7jMBhF5I59PFs7eeMjZyYeIyB34fLLoH9URz8TfbinMW7pKLHe9IyJv5NOjocycOdOZy2wQkTdispABl9kgIm/j881QREQkzSdrFuZJc2VXahDEZiIiIkk+lyzMk+a8daVZIiI5+FwzFCfNERE5zueSBSfNERE5zueSBSfNERE5zueSBSfNERE5zuc6uOtPmuNoKCIi+/hcsgD+nDQXEhKICxcqXB0OEZHb87lmKCIichyTBRERSWKyICIiSUwWREQkySs6uNVqlUs+6wqMV16eFi/geTEzXnnZE29zfieVIAhCcwIiIiLfwWYoIiKSxGRBRESSmCyIiEgSkwUREUlisiAiIklMFkREJInJgoiIJDFZEBGRJCYLIiKS5BXLfTQlLi4OWq0W/v5126YmJSVh4MCBVsdUVVXhtddeQ15eHjQaDV599VU89NBDisf63//+Fy+88ILl54qKClRWVuI///mP1XEffPABPvvsM4SGhgIA7rnnHqSlpSkSo8FgwLZt23D+/HlkZWWhR48eAIDTp08jJSUFly5dgk6ng8FgQLdu3Rp93mg0Ij09HXv27IFKpcL06dMxZswYxWMuLy/H3LlzcfbsWWi1WnTt2hULFixAUFBQo8+npKTgp59+Qvv27QEAw4YNw3PPPadovIB99zKg/P0sFq+99zKg/P3c1N/+119/RWpqKmpqatC5c2csXrwYwcHBjc6h5HdsK97Lly8jNTUVFy5cgJ+fH+68806kpaWhdevWjc4xYcIEFBQUICAgAAAwceJE/PWvf3UsEMHLPfTQQ8Lx48ebPOaDDz4Q/v73vwuCIAinT58W7r//fqGyslKJ8JqUnp4uvPHGG41eX7ZsmfDOO++4ICJB+OWXX4SCgoJG3+uECROEzMxMQRAEITMzU5gwYYLo57/66ith8uTJgtFoFEpLS4WBAwcK586dUzzm8vJyYd++fZZj3nnnHeG1114T/fyrr74qfPrpp7LGWJ+t79iee1kQlL+fbcVbn617WRCUv59t/e1NJpMwZMgQ4ZdffhEEQRBWrFghpKSkiJ5Dye/YVrznzp0T8vLyBEEQBKPRKLz44ovC8uXLRc8xfvx4YdeuXS2Kg81QAHJycpCYmAgA6NatG3r16oXdu3e7NKbr168jKyvL8ewvs5iYGOj1eqvXSktLkZ+fj4SEBABAQkIC8vPzUVZW1ujzW7duxZgxY6BWqxEUFIQhQ4bg22+/VTxmnU6Hfv36WX6+++67UVBQIGsc9hKL1xFK389S8brbvWzrb3/48GH4+/sjJiYGAJCYmGjz3lTyO7YVb5cuXXDHHXcAANRqNXr37i3rPewTySIpKQkjRozA/PnzceXKlUbvFxQUoHPnzpaf9Xo9ioqKlAyxkV27diEsLAxRUVGi72/ZsgUjRozA5MmTcejQIYWjs1ZYWIiwsDBoNBoAgEajQWhoKAoLC0WP7dSpk+Vnd/iuTSYTNmzYgLi4OJvHrFmzBiNGjMDzzz+PU6dOKRidNal7GXC/+1nqXgZcdz/X/9s3vDeDgoJgMplw6dKlRp9z1Xds616trq7Gpk2bmryHFy1ahBEjRiApKQnFxcUOX9vrk8X69evxzTffYNOmTRAEAQsWLHB1SHbZtGmTzSexxMREfPfdd8jKysKUKVPw/PPPo7y8XOEIvcebb76Jtm3bYvz48aLvv/zyy9ixYweysrLw8MMPY+rUqTAajQpH6Z33MuDa+1nqb+9uxOKtra3Fyy+/jPvuuw+DBw8W/dyiRYuQk5ODzMxM3HrrrXjppZccvrbXJwtz9Vir1WLcuHE4ePBgo2M6deqE8+fPW34uLCxEx44dFYuxoeLiYvzyyy8YMWKE6PshISFo1aoVAGDAgAHQ6/U4efKkkiFa0ev1KC4uthSgRqMRJSUlok0Ter3eqqrs6u/aYDDgzJkzeO+996BWi/9zCAsLs7w3atQoXLt2zSVP6vbcy4B73c9S9zLguvu54d++4b1ZVlYGlUoFnU7X6LOu+I7F7lWj0YikpCTcfPPNmDdvns3Pmu8djUaDiRMn4rfffoPJZHLo+l6dLK5du4aKigoAgCAI2Lp1KyIjIxsdN2zYMPzrX/8CAPzxxx84fPiw6CgTpXz11VeIjY21jL5pqH4V8ujRozh//jy6d++uVHiNBAcHIzIyEtnZ2QCA7OxsREZGio4sGjZsGDZu3AiTyYSysjLs3LkTjzzyiNIhAwCWLl2KI0eOYMWKFdBqtTaPq/9979mzB2q1GmFhYUqEaGHvvQy41/0sdS8Drrmfxf72vXr1QnV1Nfbv3w8A+PzzzxEfHy/6eaW/Y7F4TSYTUlJSoNFosHDhQqhU4hsa1dbW4uLFi5aft2zZgh49eth8OLLFqzc/OnfuHGbNmgWj0QiTyYSIiAjMmzcPoaGhGDlyJFatWoWwsDBcu3YNKSkpOHr0KNRqNZKTkzFkyBCXxf3II4/g9ddfx4MPPmh5bdq0aZg9ezbuvPNOvPrqq8jLy4NarUarVq0we/ZsxMbGKhJbeno6tm/fjosXL6J9+/bQ6XTYsmULTp06hZSUFFy5cgXt2rWDwWDArbfe2ih2o9GIBQsWYO/evZb3xo4dq3jM7733HhISEtCtWzfLUMMuXbpgxYoVAGB1fzz77LMoLS2FSqVCQEAA5s6di7vvvlvReDMyMmzeyw3jVfp+tnVPAOL3MuDa+/nkyZM2//YHDx5EWlqa1dDZDh06AHDdd2wr3jFjxuBvf/ubVcFvHnZcXFyM6dOn4+uvv8a1a9cwfvx43LhxAwAQGhqK119/3fLv015enSyIiMg5vLoZioiInIPJgoiIJDFZEBGRJCYLIiKSxGRBRESSmCyIiEgSkwUREUlisiAiIklMFkROdPbsWdx7773Iy8sDULeURb9+/fDzzz+7ODKilmGyIHKiW265BUlJSUhKSkJVVRX+/ve/Y/To0Vb7ERB5Ii73QSSDGTNmWFYl3bRpU5MLFRJ5AtYsiGTw5JNP4sSJE5gwYQITBXkF1iyInOzq1asYOXIk+vXrh927dyMrK0t0TwQiT8KaBZGTLVy4EFFRUVi4cCEGDRqEtLQ0V4dE1GJMFkROtHPnTuzZswdvvPEGACAlJQX5+fn45ptvXBwZUcuwGYqIiCSxZkFERJKYLIiISBKTBRERSWKyICIiSUwWREQkicmCiIgkMVkQEZEkJgsiIpLEZEFERJL+P54YrVbfOeVCAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x=data[0]; y=data[1]\n", "plt.scatter(x,y)\n", "plt.xlabel('x'); plt.ylabel('y')\n", "plt.title('Eksperimentalni podaci')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Linearna regresija:\n", "$y=ax+b$ gdje je \n", "$a=\\dfrac{\\Sigma_i^N (x_i-\\bar{x})(y_i-\\bar{y})}{\\Sigma_i^N (x_i-\\bar{x})^2}$ i $b=\\bar{y}-a\\bar{x}.$" ] }, { "cell_type": "code", "execution_count": 99, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEcCAYAAAA2g5hwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3XlYVGX7wPHvzLBvwiACiruCqKm4pCjuOyKab6a5pmZlppVpWb/UNK3Uyt5Sc0ktW+zNMmMRcSf33HLf9x1lhn2AWc7vD2IS2ZFhZuD5XFfXJXPOPOcehp77nGeVSZIkIQiCIAiFkJs7AEEQBMHyiWQhCIIgFEkkC0EQBKFIIlkIgiAIRRLJQhAEQSiSSBaCIAhCkUSyEKxGQEAA169fL/NyZ86cyZIlS8q83KJEREQwduxY489BQUHcvHkTgOnTp7No0aJilXPr1i0CAgLQ6XQmiTPHwYMH6dSpk0mvIVgukSyEInXr1o2mTZuiUqlyvT5gwAACAgK4detWkWVYckUzZ84cJk6cWO7XDQ8PZ/Xq1cafjx07Rs2aNcs9jvKyYcMGnn/+eXOHIZSSSBZCsdSoUYPo6Gjjz+fPnycjI8OMEQmWxNRPNYL5iWQhFMuAAQPYuHGj8eeNGzcycODAXOdkZWUxf/58unTpQvv27Zk5cyYZGRmkp6czfvx44uPjCQoKIigoiPv375OVlcW8efMICQkhJCSEefPmkZWVZSzvm2++MR779ddfC41v5MiRfPHFFwwdOpSgoCDGjh2b60lo8uTJdOjQgVatWjF8+HAuXrxoPPZok0/OE9DKlSsJDg4mJCSEbdu2ERcXR+/evXn66adZtmyZ8b0nTpxgyJAhtG7dmpCQEObMmZPrMwQEBLBu3Tp69epFmzZtmD17NjmLJjx+p/14M1tycjIvvfQSQUFBDB48mBs3bhT+Jf0jNjaWbt26ceHCBcaNG8cPP/yQ63h4eDhbtmwBYO7cuXTu3JmWLVsyaNAgDh8+bDwvIyOD6dOn06ZNG0JDQzl58mSucrp168aKFSvo378/LVq0QKfTsWLFCnr06EFQUBChoaFs3boVgMuXLzNr1iz+/vtvgoKCaN26NQC7du1i4MCBtGzZks6dO/PVV18V6zMKZiAJQhG6du0q7d27V+rVq5d06dIlSafTSZ06dZJu3bol+fv7Szdv3pQkSZLmzp0rvfzyy5JarZZSUlKkl19+Wfr0008lSZKkAwcOSB07dsxV7hdffCENHjxYevjwoZSQkCANGTJEWrRokSRJkhQXFycFBwdL58+fl9LS0qQpU6ZI/v7+0rVr1/KNccSIEVL37t2lK1euSBqNRhoxYoS0cOFC4/H169dLKSkpUmZmpjR37lwpPDzceOydd96RPv/8c2OcgYGB0ldffSVlZWVJ//vf/6S2bdtKU6ZMkVJSUqQLFy5ITZs2lW7cuCFJkiSdPHlSOnbsmKTVaqWbN29Kffr0kdasWWMs29/fX3rppZekpKQk6fbt21Lbtm2luLg4SZIk6bfffpOGDh2a69ycz/fOO+9Ibdq0kY4fPy5ptVppypQp0htvvJHvZ79586bk7+8vabVa6ddff5V69OhhLOf333+XhgwZYjz34sWLUqtWraTMzExJkiRp48aNkkqlkrRarbRq1Sqpffv2UkZGhiRJkrRw4ULp+eefl9RqtXTnzh2pX79+ub7Drl27SuHh4dKdO3ckjUYjSZIkbdq0Sbp3756k1+ul6OhoqXnz5tL9+/fz/bw5v+9z585Jer1eOnv2rBQcHCxt3bo1388pmJd4shCKLefpYu/evdSrVw9vb2/jMUmSWL9+Pe+99x7u7u64uLjw8ssv52q6elxkZCQTJ07E09MTpVLJxIkTiYiIACAmJoZBgwbh7++Pk5MTr732WpHxDRo0iLp16+Lg4ECfPn04e/as8dizzz6Li4sLdnZ2TJo0iXPnzpGSkpJvOTY2NkyYMAFbW1tCQ0NRq9WMGjUKFxcXGjZsSMOGDTl//jwATZs2pUWLFtjY2ODn58eQIUM4dOhQrvLGjx+Pm5sb1atXp23btpw7d67IzwLQs2dPmjVrho2NDeHh4bk+T36+++47Vq1axffff0/t2rUB6NGjB+fOneP27dtA9u+8Z8+e2NnZAdnfqYeHBzY2NowdO5asrCyuXr0KZH8Hr7zyCu7u7vj6+jJy5Mg81xw5ciS+vr44ODgA0LdvX7y9vZHL5YSGhlK7dm1OnDhRYMxt27YlICAAuVxOo0aN6NevH3/99Vexfj9C+bIxdwCC9RgwYAAjRozg1q1bDBgwINcxlUqFRqNh0KBBxtckScJgMBRYXnx8PNWrVzf+XL16deLj443HmjZtajxWo0aNIuPz8vIy/tvR0ZH09HQA9Ho9ixYtYvPmzahUKuTy7HsktVqNq6trnnLc3d1RKBQAxkrQ09PTeNze3p60tDQArl69yieffMKpU6fQaDTo9XqaNGlSaFw57y1K1apVjf92cHAwfp6CrFq1iokTJ+Lj42N8zcXFhc6dOxMdHc1LL71EdHQ0H374ofH46tWrWb9+PfHx8chkMlJTU1Gr1UD2d+Dr62s899HvKsejxyG7eXLNmjXG5JSenm4sLz/Hjx/n008/5eLFi2i1WrKysujTp0+hn1MwD/FkIRRbjRo18PPzIy4ujl69euU65uHhgYODA9HR0Rw+fJjDhw9z5MgRjh07BoBMJstTXrVq1bhz547x57t371KtWjXjsbt37xqPPXpeSUVGRrJ9+3bWrFnDkSNH2LFjB4Cx7+BJfPDBB9SrV4/Y2FiOHj3Km2++WSbllsbq1av5+uuviY2NzfV6WFgY0dHRHDt2jIyMDNq2bQvA4cOHWblyJV988QWHDh3i8OHDuLq6GuP38vLK9R08+u8cj36vt2/f5v3332fGjBkcPHiQw4cP07Bhw3zPzfHWW2/RvXt34uLiOHLkCEOHDjXb708onEgWQonMmzeP7777Dicnp1yvy+VyBg8ezEcffURCQgIA9+/fZ/fu3UD2nXliYmKupp9+/frx9ddfo1KpUKlULFmyhP79+wPQp08ffv/9dy5duoRGo2Hx4sWljjktLQ07Ozs8PDzQaDR8/vnnpS4rv7KdnZ1xdnbm8uXLrFu3rszKLqkGDRrwzTffMGfOHLZv3258vXPnzty5c4cvv/yS0NBQ45NVWloaCoUCpVKJTqdj8eLFpKamGt/Xt29fVqxYQVJSEvfu3eP7778v9PoajQaZTIZSqQTgt99+yzWQwNPT0ziwIUdaWhpVqlTB3t6eEydOEBUVVSa/C6HsiWQhlEitWrV46qmn8j02bdo0ateuzXPPPUfLli154YUXjO3f9evXp1+/fvTo0YPWrVtz//59Xn31VZo2bUp4eDjh4eE0adKEV199Fciu4EaPHs3o0aPp2bMn7dq1K3XMAwcOpHr16nTs2JF+/frRokWLUpf1uHfeeYeoqChatmzJjBkzCA0NLbOyS6NRo0YsW7aMGTNmEBcXB4CdnR09e/Zk3759hIWFGc8NCQmhU6dO9O7dm27dumFvb5+rWem1116jevXqdO/enbFjx+ZpenxcgwYNGDt2LEOHDqV9+/ZcuHCBli1bGo+3a9eOBg0aEBISYny6mTVrFl9++SVBQUEsWbKEvn37luWvQyhDMkk88wmV3Ntvv02tWrWK1Yleln799VciIiJYu3ZtuV5XEEpDPFkIlZpOp+Pq1av4+fmV+7UvXbpklusKQmmIZCFUah06dMDNzY3evXuX63VfffVVdu/enWttKEGwZKIZShAEQSiSeLIQBEEQiiSShSAIglCkcpnBrVarefvtt7lx4wZ2dnbUrl2bOXPmoFQqCQgIwN/f3zj2e8GCBQQEBJRHWIIgCEIxlUufRWJiIufPnzeOrZ4/fz5JSUl89NFHBAQEcPToUZydnUtdvlqdhsFQ8o/h6elCQkJq0SdaCBGvaVlbvGB9MYt4Tau48crlMjw8SlbnlsuThbu7uzFRALRo0aJMZ7oaDFKpkkXOe62JiNe0rC1esL6YRbymZap4y30hQYPBwLp16+jWrZvxtZEjR6LX6+nUqROTJk0yrogpCIIgWIZyHzo7e/Zs7t+/z+LFi5HL5dy9exdfX19SU1OZNm0a/v7+vPnmm+UZkiAIglCEcn2ymD9/PtevX2fZsmXGDu2ctWhcXFwYPHgwa9asKXG5CQmpuR69JElCrX5AVlYGUHAulMvlhS6hbWksO14ZdnYOeHh4GVcX9fJy5cGD/PeMsETWFi9YX8wiXtMqbrxyuQxPT5cSlV1uyWLRokWcOnWKFStWGJuZkpKSsLe3x8HBAZ1OR2xsLIGBgU98rdTUJGQyGd7efshkBY8OtrGRo9NZauWblyXHK0kGEhMfkpqahKuru7nDEQShjJVLsrh48SLLli2jTp06DB06FAA/Pz9efPFFZs6ciUwmQ6fTERQUxOuvv/7E19NoUlEqvQtNFELZksnkuLp6oFLdF8lCEExo/+l7bIi7TEJyJp5u9gzqXJ/gJj5Fv/EJlUuyeHQbysdFRkaW+fUMBj0KhdgEsLwpFDYYDHpzhyEIFdb+0/f4LuYcWf+0MCQkZ/JdTPY2vaZOGBX21ju/XbkE0xK/c0EwrQ1xl42JIkeWzsCGuMsmv3aFTRaCIAgVTUJyZoleL0siWQiCIFgJTzf7Er1elkSyEIqUnJzEgAHlu9+DIAh5DepcHzub3NW2nY2cQZ3rm/zaIllUcs8+27/Ic86ePUOjRo3LIRpBEAoT3MSH0X0bGZ8kPN3sGd23UcUZDVXZqdUqhg59hpiYncbJiG+9NZmwsHC6du1RqjJPnTrJkiVfkJmZQXp6Os89N4xBgwYD8NFHs/Hx8WXs2Je4efMG06a9wezZHxEQ0KhU1zp37gxublV4882J3LlzG3//RsyaNRcbG/HnIwjlLbiJT7kkh8dViv/bk/ftJWnPn3lel8lkPOlqJ1VCOuHWvkOh53h4KFEqPbly5TINGjRk+/atyGTkShTvv/82t27dyvf9y5evxt7eIddrtWrVZvHiFSgUClJTU3n22f7GZDF+/ATGjRtJSEgn5syZwbvvzix1ogA4d+4sWm0Wn3zyOQqFgokTx7N7965SJzpBEKxPpUgWlqB58yBOnTqOr68vK1YsYdGiJbmOz527oETl7d+/h8jIjaSkpGAw6I1PLABeXtXo06cfEyeOZ+7cBTRv3iLXe995503u378PwMOHD3jhhWEAKBQKVq36Ps+1zp8/y+LFK7C3z370bdjQH7VaXaJ4BUGwbpUiWbi175Dv3X95Lp/RvHkQR44c4urVK/TrF0716jVyHS/Jk8XevbvZuPFX5s5dgKdnVXbs2MYff2wwHlerVRw4sBcHB0eqVfPOU978+YuM/3722f58++1PBcatVqvIysrMFe+ZM6cJDS26r0MQhIqjUiQLS9CsWQu++moRVatWZdWqH/IcL8mTxYUL5wgICMTTsyoPHsSzcuVSOnXqCkBKSgpvvTWZsWNfIjFRzddff8WCBYuKKLFgZ8+eITU1lXv37uHj48OGDetxd3enceOmpS5TEEzBXMtgVBYiWZQTHx9ftFotb775Nra2tk9UVu/eoUyfPoUXXhhGQEAjfH1rEBAQSEZGBm+//QaDBj1Lly7d0el0/PLLOo4ePUzLlq1Lda3z58/yzDODmTt3JomJaho08Gf27I+fKH5BKGvmXAajsij3/SxM4fElyu/du46PT+0i31eezVA///wD165dZfr0GaUuw5JXnc3x6O++oi7vbEmsLWZTxTtt6d58ZzF7utmz8NXCB6AUpqL+fkuzRLmYZ2Fi169fY/jwZ/nrr4NMnvyWucMRhArJnMtgmIMhMxNDhqZcrymaoUysdu06/Pjjr+YOQxAqNE83+wKfLCoSfXoaidu3od66BYfadfB7a1q5XVskC0EQrN6gzvVz9VlA+S2DUR70qamot8aSuGMbBo0G5xZBVH3mP+Uag0gWgiBYvZxObFOOhjLHaCtdcjLqLZtJ3LkdKTMTl1atUfbrj0Otovtky5pIFoIgVAimXAajvEdb6RITUcXGkBS3E0mrxbVNW5T9+mNfo0bRbzYRkSwEQRCKUNimQ2WZLLSqBNSbN5H0ZxySwYBb22CU/cKw8/Ets2uUlkgWgiAIRTD1aCvtwweoYqJJ2rMbALfgDihDw7CrVq1Myi8LIlkIgiAUwVSjrbLu30e1KYrkA/uQyWRU6dgZZd9QbD2rPlG5piCShSAIQhHKerRV5p07qKIjSfnrADIbG9y7dMOjTyi2Hh5lFXKZE8nCShw9epilS//LN9/kXRW2OPbsieP48b+ZOPH1Mo5MECq+nH6Jn7aeJy1DD4CdbcnnNGfeuklCVCSpRw4hs7XFo1dvPHr1waaKe5nGawoiWTyiIi9EFhLSmZCQzuYOQxCsmlb377JCqRpdsUdEZdy4jioygtRjR5DZO+DRJxSPXr2xcXUzabxlSSSLf5hyaFxISGvGjBnPoUMHSUpK5OWXJ9KlS3cADhzYx/LlizEYDLi7ezBt2nv4+dUEYMWKpWzfvgUvr2oEBjYxlpeQ8JAPPvg/0tLSyMrKon37Drz6avYTg1arZfnyJRw8uA+5XEH16jX4+ONP2bQpkn37dpd43wxBELKVZkSU5soVVFF/kHbiOHJHR5Rh4Xj06IXCpWTrMlkCkSz+YeqhcXK5nGXLVnPjxjVeeWUczZsHATB37ky++moFdevWIypqI7Nnv8/Kld+xZ8+f7N37J2vW/IS9vT3vvjvVWJaLiyvz5y/CyckJnU7HlCmvceDAPtq1a8/336/hzp3brF79I7a2tiQmJj5x7IIglGxElObiRRKi/iD99Cnkzs54DhyEe7ceKJycTB2myYhk8Q9TD40LCxsAQK1adfD3D+D06ZOAjPr1/albtx4AoaHhfPbZfNLT0zh27DDduvXE6Z8/rrCwAaxduwoAg8HA0qX/5eTJE4BEQkICFy9eoF279uzbt4fXXnvDuAy6u7vlt4UKgjUoakSUJElozp8jISoCzbmzKFxdqfqf53Dv2hW5g2N5h1vmRLL4R3kuRJa9KLwMkJDJCjqn4JXj//e/H0lJSWbFim+xt7dn/vx5ZGVlFvk+QRBKr8ARUZ3qkXb6FKqoCDQXL6CoUgWv54ZSpXNX5PYVZyFDsUT5PwZ1ro+dTe5fR1kuRBYdHQHAzZs3uHTpPE2aNKVJk2ZcunSB69evARATE0XDhgE4OTnTqtXT7NixDY1Gg16vZ9OmCGNZKSkpeHpWxd7engcP4tmzJ854rEOHjvzyyzq0Wi2AaIYShDIS3MSH0X0bGW8gPV3teDEAakSs5PaiT9E+fIDX88Op+/FCPHr1qVCJAsrpyUKtVvP2229z48YN7OzsqF27NnPmzEGpVPL3338zc+ZMMjMzqVGjBgsXLsTT07M8wsrF1AuR2dnZMWHCWBITE5k27T08PJQAvP/+HGbP/j/0ej3u7h7MnPkhkF3pnzp1gjFjhlG1qhdBQa14+PABAIMHD2XGjHcYM2YY1ap506pVG+N1Rox4geXLFzNmzDBsbGzx8/MTndqCUEaCm/jQLrAaqX8fQxUVQeax6+g8Pak2cjRu7UOQP+EumJasXHbKS0xM5Pz587Rt2xaA+fPnk5SUxLx58+jVqxcff/wxrVu3ZunSpdy8eZOPPy7Ztp2WvlNeSEhrtmz509j/UFpipzzTsrZ4wfpituZ4JYOB1COHSYiKIOv2LWy9qqHsF4Zbu/bIbCyjRd+UO+WVyyd0d3c3JgqAFi1asG7dOk6ePIm9vT2tW2fvDz106FC6d+9e4mQhCIJgKpLBQMqhg6iiIsm6ewdbHx98xr2E69NtkSkU5g6v3JR7OjQYDKxbt45u3bpx9+5dqlevbjymVCoxGAwkJiaWaBTP4xkyPl6OjU3xumOKe96TOHDgaJmVVR7xPgm5XI6Xl6vx50f/bQ2sLV6wvpitJV6DTsf97Tu49esGMu7cxal2LepMnULV9u0sOkmY6vdb7sniww8/xMnJiREjRrB169YyKfPxZiiDwYBWq0dW0FCjf1hDs86jLD1eSZIwGAzGx2BrbnKwFtYWszXEK+l0JO3bg3pTNNqHD7CvWQvfCa/hEtQS5HIeqtLNHWKBrL4ZKsf8+fO5fv06y5YtQy6X4+vry507d4zHVSoVMpnsiecGyOUK9HodNjYVt7PJEun1OuRyy73jEoTCGLRZJO/ZjSomGp1KhX2dugS+PA5dnYAibzwrg3JLFosWLeLUqVOsWLECOzs7AJo2bUpGRgaHDx+mdevW/Pzzz/Tt2/eJr+Xo6EJKSiLu7p7IZJbdbFNRSJKBlBQ1jo7Wt4yBULkZMjNJ2h2HavMm9ImJONRvgPeoMTg1aYqympvFPwmVl3JJFhcvXmTZsmXUqVOHoUOHAuDn58eSJUtYsGABs2bNyjV09km5uFRBrX7A/fu3gIIHe8nlcgwGy23WeZxlxyvDzs4BF5cq5g5EEIrFkJFBYtxO1Jtj0Kck4+gfgOe4l3BsFCieJPJRLsmiYcOGnD9/Pt9jLVu2JDIyskyvJ5PJUCqL3mHKGtpPH2Vt8QqCJdJrNCTu2IZ6ayyG1FScApug7B+Ok3+AuUOzaJYxOFgQBMHE9GlpJG7finrbFgzp6Tg/1QxlWDiO9RuYOzSrIJKFIAgVmj4lBfXWWBJ3bseg0eDcIgjPsHAc6tQ1d2hWRSQLQRAqJF1SEuotm0nctQMpKwuXlq3wDBuAfc2a5g7NKolkIQhChaJLVKPaHEPSn7uQtFpcn26Lsl9/7KvXMHdoVk0kC0EQKgStKgFVzCaSd8chGQy4tQtGGdofO5+KsTWyuYlkIQiCVdM+eIAqJpqkvbsBcGvfAWVoGHZeRY+IFIpPJAtBEKxS1v17qKKjSD6wD5lcTpVOnVH26YetGbY4qAxEshAEwapk3rmDKjqSlL8OILOxwb1bdzx6h2Lr4WHu0Ipl/+l7Jts3x5REshAEwSpk3rxJQnQEqUcOI7O1xaNXbzx69cGmivXsM7//9L1cW7MmJGfyXcw5AItPGJU6Wew6cpNvo05bXYYXhMok4/o1EqIiSDt2FLmDA8q+/fDo2RuFq3Usdf6oDXGXc+3hDZClM7Ah7rLF1z2VNlnsP32PtZvPk6nVA9aV4QWhMtBcuYwqKoK0E8eROzmh7D8Aj+49UbhY72KVCcmZJXrdklTaZLEh7rIxUeSwlgwvCBWZ5uIFEqIiSD99CrmzM54DB+HerQeKJ9yW2BJ4utnnmxg83ezNEE3JVNpkYc0ZXhAqGkmS0Jw/R0LkH2jOn0Ph6krV/zyHe9euyB0czR1emRnUuX6uPgsAOxs5gzrXN2NUxVNpk4U1Z3hBqCgkSSL9zGlUURFoLl5AUcUdr+eep0rnLsjtS/7/oqWPNMqJxZJjLEilTRaDOtfP1WcB1pPhBcHaSZJE2onjqKIiyLh6BRsPJdWGjcCtYyfktnalKtNaRhoFN/GxqHiKq9Imi+AmPri5OojRUIJQjiSDgZSjR1BFRZB54zo2VatSbeQLuLXvgNz2ybZBtuaRRtag0iYLgC6tatKklvWM0RYEayUZDKQeOcytzVGkX7+BbTVvvF8Yh1u7YGQ2ZVMNiX5I06rUyUIQBNOS9HpS/jqIKjqSrHt3cfSrgc+LL+Hapi0yhaJMryX6IU1LJAtBEMqcpNORfGA/qk1RaOPvY1fDD9+XX6Vu7y48VKWb5JrWPNLIGohkIQhCmTFotSTv24sqJgrdw4fY16qN76uTcGkRhEwuL/OniUdZ80gjayCShSAIT8ygzSJp95+oYzahU6twqFuPasNG4PxUc2QyWbnFYa0jjayBSBaCIJSaITOTpLhdqGI3oU9KwrGhP94vjMWpcZNyTRKC6YlkIQhCiRkyNCTu3Il6y2b0Kck4NgrEc/wrOAY0EkmighLJQhCEYtOnp5O4YxvqrbEY0tJwatIUz7BwHBv6mzs0wcREshCEcmLpS1EURp+ainr7VhK3b8WQno5zs+Yow8JxrCdGGlUWIlkIQjmwlqUoHqdPSUG9NZbEHdswZGTgHNQSz7BwHGrXMXdoQjkTyUIQyoG1LUWhS0pCvSWGxJ07kLRaXFq1xrNfOPY1a5o7NMFMRLIQhHJgLUtRaNVq1Js3kfTnLiSdDten26Hs1x/76tXNHZpgZuWWLObPn09sbCy3b98mMjISf//sDrFu3bphZ2eH/T/LEU+dOpWOHTuWV1iCUC4sfSkKbcJDVDGbSN7zJ5LBgFtwB5Sh/bDztrynHsE8yi1ZdO/enVGjRjF8+PA8x7788ktj8hCEishSl6LIehCPalMUyfv2AlClQ0eUffth6+Vl1rgEy1NuyaJ169bldSlBMKnSjGqytKUosu7dQ7UpkuQD+5HJ5VTp1AVl31BslZ5miUewfDJJkqTyvGC3bt1YtmxZrmYoFxcXJEmiVatWTJkyBTc3t/IMSRCKbdeRmyxefzzXpln2tgpeG9ycLq0sv/M3/cZNbq7/jYd79iK3scG7dy9qPDMAe0+luUMTLJzZO7h//PFHfH19ycrKYt68ecyZM4dPP/20RGUkJKRiMJQ853l5ufLgQUqJ32cuIl7TKk6830adzpUoADK1er6NOm2WvVGK+zvOvHmDhKgIUo8eQWZnh0evPnj07I1NlSokG4By+p4q4t+EJSluvHK5DE9PlxKVbfZk4evrC4CdnR3Dhg1jwoQJZo5IEApmLaOacmRcu0ZC1B+k/X0MuaMjytAwPHr0QuHqau7QBCtj1mSRnp6OXq/H1dUVSZLYtGkTgYGB5gxJEApl6aOacmguXyIhMoL0UyeQOznhGT4Q9+49UTg7mzs0wUqVW7KYO3cuW7Zs4eHDh4wZMwZ3d3eWLVvGpEmT0Ov1GAwG6tevz6xZs8orJEEoMUsd1ZQj/cJ5VJERpJ89jdzFhaqDnqVK1+4oHB3NHZpg5cq9g9sURJ+FZaqo8VrSGk9eXq7ExyejOXeWhMg/0Fw4j8LVDY8+fXHv3BW5g4NZ4ipIRf2bsBQVus9CEKyNpWywI0kS6qPHuPnDz2RcvoRSfAY/AAAgAElEQVTC3R2vocOo0rEzcnvLahYTrJ9IFpWIJd0RC6UnSRJpx/8mISqCzGtXsVEqqTZ8JG4hHZHb2pk7PKGCEsmikrDWVU/NydKSq2QwkHrsCKqoSDJv3sC2qhf1J05A/lQrZDbif2XBtMRfWCVhbauempslJVfJYCDl8F+ooiLJunMbW29vvMeMw61tMNV8PayqTV2wXiJZVBLWNj/A3CwhuUp6PSl/HSAhOhLtvXvY+VbHZ/zLuLZpi0wuL5cYBCGHSBaVhLXMD7AU5kyukk5H8v69qDZFo30Qj51fTXxfmYhLy1YiSQhmI5JFJWHp8wMsjTmSq0GrJXnvblQx0egSErCvVZvqEyfh3DxIJAnB7ESyqCQsbdVTS1eeydWQlUXS7jjUmzehU6txqFePasNH4fxUM2QyWZlfTxBKQySLSsRS5gdYg/JIrobMTJLidqKKjUGflIRjQ3+8x7yIU2BjkSQEiyOShSAUwFTJ1ZChIXHnDtSxm9GnpuDYKBDP8a/g1EisiyZYLpEsBKGc6NPTSNy+DfW2LRjS0nBq0hTPsAE4Nmxo7tAEoUgiWQiCielTU1Fv20Li9q0YNBqcm7dA2S8cx3r1zB2aIBSbSBZCoSxtFrM10aUko94SS+KO7UiZGbi0bIUyLByHWrXNHZoglJhIFpVMSSp/S5rFbE10iYmoY2NIjNuJpNXi2uZplP36Y1/D74nLfvz7eyGsiVl26BMqH5EszMzUd+6Plu/soCBTa0Cnz17OvajK3xJmMVsTrUqFevMmknbHIel0uLYLxjM0DDvf6mVSfn7Je/H644zqEyC+D8HkRLIwI1PfuT9eflqGPs85hVX+YomQ4tEmPES1KZrkvbuRJAm34PYo+4Zh5+1dptfJL3lnavUieQvlQiQLMzL1nXt+5eenoMpfLBFSuKz4eFSbIknevw+AKiEdUfbth21VL5NcTyRvwZxEsjAjU//PX9xyCqr8xRIh+cu6d5eE6EhSDh5AJpfj3rkLHn1CsVV6mvS6InkL5iSShRmZ+n/+gsp/VGGVv1giJLfM27dQRUeScugvZLa2eHTviUfvvti4l08Hc37J295WUemTt1A+RLIwI1PfuedXvkIGjg42pGp0xar8xRIhkHHjOqqoCFKPHkFm74BH77549OqDjZtbucaRX/IWo6GE8iKShRmZ+s5dPBk8mYyrV0iIiiDt+N/IHR1RhvXHo0dvFC4l2+i+LD2evL28XMXmR0K5EMnCzEx95y6eDEpOc/kSCZF/kH7qJHInZzwHPIN79x4onJzNHZogmI1IFoLwj6RTp7n1w8+knz2DwsWVqoOepUrX7igcHc0dmiCYXbGTxccff8zAgQMJDBQrYwoVhyRJpJ89gyoqAs2F8yjc3Kg6eAjuXbohtxejjAQhR7GThU6nY9y4cSiVSsLDwwkPD8fHRzRvCNZJkiTST50kISqCjMuXsPHwoO6LY1G0bIfczs7c4QmCxZFJkiQV92S9Xs+ff/5JZGQkO3fupHnz5gwcOJCePXvi7Gy+9tyEhFQMhmJ/DCNr6xwU8T45SZJI+/sYCdGRZF67io3SE2VoP9w6dMS7utLi4i2KJf6OCyPiNa3ixiuXy/D0LNlAjRIli0ddvHiRt956iwsXLuDo6EhoaCiTJ0/Gu4yXOCgOkSwskyXFKxkMpB49gio6gsybN7H18kIZGoZbcAdkNtkP2JYUb3FZW8wiXtMyZbIoUQd3amoqMTExREREcOHCBXr16sWsWbOoXr06q1ev5sUXXyQyMjLP++bPn09sbCy3b98mMjISf39/AK5evcr06dNJTEzE3d2d+fPnU6dOnRJ9AEEojGQwkHLoL1TREWTduYOttw8+Y8fj2rYdMoWi0PeK5dkF4V/FThaTJ09m9+7dtGnThueff54ePXpg90jb7rvvvkurVq3yfW/37t0ZNWoUw4cPz/X6rFmzGDZsGAMGDOCPP/5g5syZrF27tpQfRRD+Jen1pBw8QEJ0JNr797CrXh2fl17BtfXTyOTyIt8vlmcXhNyKnSyaN2/OjBkz8PLKf5E0uVzOvn378j3WunXrPK8lJCRw5swZ1qxZA0BYWBgffvghKpUKpVJZ3LAEIRdJpyN5315UMVFoHzzAvmZNfCdMxCWoVbGSRA6xPLsg5FbsZDFu3Lgiz3EswXj0u3fv4u3tjeKfpgCFQkG1atW4e/euSBZCiRm0WpL37EYVE41OlYB9nbpUHzIM5+YtkMlkJS5PrPAqCLlViEl5Je2oeZSXl2sZRmJ6It7c9JmZ3N+yjdsbNpKlUuEaEEDNSRNwDypdksiJ18vDkQdqTd7jHo4W9x1YWjxFEfGalqniNVuy8PX15f79++j1ehQKBXq9nvj4eHx9fUtclhgNZZlMGa8hI4PEuJ2oY2PQJyfj6B+A35gXcWwUiE4m4+HD1BKX+Wi8A0Pq5rvI48CQuhb1HYi/CdOqqPGafDRUWfL09CQwMJCoqCgGDBhAVFQUgYGBoglKKJReoyFp53bUW2LRp6bgFNgE5SvhOPkHlOl1xCKMgpBbuSSLuXPnsmXLFh4+fMiYMWNwd3cnOjqaDz74gOnTp7N06VLc3NyYP39+eYRTLGLYpGXRp6eRuH0b6q1bMKSn4dS0GZ79w3Gs38Bk1xSLMArCv0o9Kc+SlHUz1OPDJiG7CWJ030ZmrTwq6iNxYfSpqai3xZK4fRsGjQbnFkF4hoXjUKduGUX5r9LGa84bi8r4N1GeKmq8VtUMZcnEsEnz0yUno96ymcSdO5AyM3Bp1Rplv/441Kpt7tByEfMxhMpCJIt8iGGT5qNLTEQVG0NS3E4krRbXNm1R9uuPfY0a5g4tX+LGQqgsRLLIh6n3xhby0qoSUG/eRNKfcUgGA25tg1H2C8POp+Sj48qTuLEQKguRLPJh6r2xhX9pHz5AtSmapL27AXAL7oAyNAy7atXMHFnxiBsLobIQySIfYtik6WXdv49qUxTJB/Yhk8mo0rEzyr6h2HpWNXdoJSJuLITKQiSLAohhk6aRdfcOCdGRpBw8gMzGBvcu3fDoE4qth4e5QysVcWMhVBYiWQjlIvP2LVRREaQcPoTM1haPnr3x6N0Hmyru5g7tiYkbC6EyEMlCMKmMG9dRRUaQeuwIMnsHPPqE4tGrNzaubmV+LTGRUhBMRyQLwSQ0V65wZnk06kNHMNg7cMynJbsdG+J8341BN9IJblK2yULMdxAE0xLJQihTmosXSYj6g/TTp7BxdSE9uCdrHlYlRbIFIMNElbiY7yAIpiWShfDE9p+6y/7oPTx16wi1NfcwODpT7T+Dqf/sAF5cGEeKlHtoqSkqcTHfQRBMSyQLodQkSeJIzB40MZGEaeJJVTiyvWprTisDGF7rKQKcHMutEhfzHQTBtESyqATKuuNXkiTSTh5HFRWB25UroHBia9U2HHdriE5uA4bsZqHwLg3LrRIX8x0EwbREsighaxtxU5Ydv5LBQNrxYyRERpB54zo2np5s9mrHSbf66GWKXOfmJIjyqsTFfAdBMC2RLErAGkfclEXHr2QwkHr0MAlRkWTduomtVzW8XxiLW7v23FxxEH0hTw6lrcRLk5TFfAdBMB2RLErAGkfcPEmfgWQwkHLoIKqoSLLu3sHWxwefceNxfbodMkX2k0RxnhxKWolbY1IWhIpOJIsSMOeIm/2n77Fxz34eqDXGO20o+o69NH0Gkk5H8sH93Pn9DxSJD4m3c+dknW60GNCDuk9Vz3VuUU8Ojz4hyGVgkCjyScEak7IgVHQiWZRATmWX3+v5NZtA2bSh53envWbTWSSDhP6feAq6+y5Jn4Gk05G0bw/qTdFoHz4gwV7JHp/OXHCuBTIZx2MvgFye5zMU9OSw68jNXNc2FBFrDjEMVhAsj0gWJVDQzq0GiVJX5sWR3522Tp83mPzuvovTZ2DQZpG8ZzeqmGh0KhX2deoS6xbEMZk3yGSFll+YtTFn88RdnLLEMFhBsDwiWZRAQZWYXEapK/PiKMkddX7nFnTnb8jMJOnPXag2x6BPSsShfgO8R43BqUlTjs3f+cSxPFRrShwriGGwgmCJRLIogYIqsYLunvNTmqaUgpJUQecWxZCRQeKuHahjN6NPScYxoBGeL76EY6NAZP88SZTF3X1VD0ceFJIwCipLDIMVBMsjkkUJFFSJ5fxcHDn9GyWp+PJLUjYKWa5mLij67luv0ZC4YxvqrbEYUlNxCmyCsn84Tv4BxbpmSe/uR/UN5Ktf/s43mRZVlhgGKwiWRSSLEiqoEitOZQ7/9m/klFXcawJs3HO1xKOhAPRpaai3bSFx+1YM6ek4P9UMZVg4jvUbFHnNJ7m779KqJskpGSUeDSUIguURyaIMFFSxAqyKOpOnY7w0fRfBTXwI79KQBw9S8r12fg4cvsz13yNpHH8ae0lLVv3GNHh+MA516hb7mk9aoYsnBEGoGCplssgZ5qpKzkRZRne5BVWKKyPP5Hu+KYeB6pKSOP3Tb7gc3UcLScc5l9rs82hGkq0no9McCTbZlQVBqKgqXbIoj9nB+U1Ee5wphoHqEtWoNseQ9Ocu7LO0nHWpwz7lUyTY/bN1qZjYJghCKVW6ZGHq2cGPJ6P8EsXjnbtPujihNiEB1eZoknf/iWQw4NYumIW3vFDb5d2N7kmeaCx9EUVLj08QrFmlSxamnh28btuFQofSOjsoGNYzINdyGKV90tE+eIAqJoqkvXsAcGvfAWVoGHZe1ZAv3QtlOLHN0tdrsvT4BMHaWUSy6NatG3Z2dtjbZ1dkU6dOpWPHjia5VkHzB1wcn/xXsf/0PVI1ukLPcbCzyVV5leZJJ+v+PVTRUSQf2IdMLqdKp84o+/TD1tPTeE5xh74W927c0tdrsvT4BMHaWUSyAPjyyy/x9/c3+XUGda7Pmk1n88yw1mToSjz/4XEb4i4Xec7jiaokTzq7Nh/i2s+/Ul99GYNMjrZFMI2HPYuth0eec4sz9LUkd+OWvl6TpccnCNbOYpJFeQlu4sNPW8+j0+tzva6XeOK70OJUTI83AxU2O3va0r0M6lyflm5aLvz0C3YXT1JbZsMh90AOujdBp3Fm9J1MgvPmCqDoYasluRu39PWaLD0+QbB2cnMHkGPq1Kn079+fDz74gOTkZJNeKy1Dn+/rT3oXWlTFlF8z0KDO9bGzyf9rsIm/g3rlUq7PnoHs8jn2ezzF13UGsbNqa9JtHI0Ve2mV5G48vzgtab0mS49PEKydTJKkAtZSLT93797F19eXrKws5s2bR1paGp9++qnJrjd27pZ81yzy8nBk9fu9cr2268hN1sac5aFaQ1UPR0b1DaRLq5r5lrvryE3++79j+S4i6Opky0sDn8r3vTnXyInJN+MBHVQnaJB+mwy5HWd8nuJP+wZkKPImIxkwZVjLYsdY2t/Do3GW9DrlxdLjEwRrZhHJ4lHnz59nwoQJ7Nixo9jvSUhIxVDQ+uH5eLytHrLvQkf3bVRom35B5z1e9rptF4wd3Y+PfirMzFnr6KA6QV3NXTRyO/5yb8zRKo3IVNgV2jGfpTWUKMYn+Xwl4eXlmmfGuSWztnjB+mIW8ZpWceOVy2V4erqUqGyz91mkp6ej1+txdXVFkiQ2bdpEYGCgSa9Z0FpL+e3u9riiRtiUdHkLSZLQnD9HQuQfjLh9jjSFAzs8W3KsSgBaua3xvPxiUciy31/cfof8Rj6N7ttIzE0QBKFIZk8WCQkJTJo0Cb1ej8FgoH79+syaNcvk1y1oraX87rYfVxYjbCRJIv30KRKiIsi4dBFFFXfSOoex6q476YbidSXJ5LJC+18eHd1V0Min0X0bsfDVDk/8eQRBqNjMnixq1qzJxo0bzR2GUX4jhB73JCNsJEki7cRxVFERZFy9go2HkmrDRuDWsRNyWzuGF/JU8zidXipwORHIvbrtk8xDEDOjBUEwe7Iwl/2n77Fxz/48zVBFVdKlHWEjGQykHjuKKjqSzBvX0bt5sK9WCAdsa+N+zpFB3qpcTVhFPd3kMEgFb8D0aDIo7TwEMTNaEASopMmisAqwsHkPpbmrlgwGUg8fIiE6kqzbt7Ct5k1qr2f55roTOS1Ij1fAxXm6eTymola3Le08BDEzWhAEqKTJorAKsKBlMko6Qmj/yduc2LiVp+4co6o2CZ2yGn4vvoRrm7YsXX6ADH3uirs4TwH5yUleBTVd5SSD0u58J2ZGC4IAFjQprzwVVgHmJBJ59lbUeLrZlyhRSDodR3+JQr70E7pe34VBJmOjdye+8urDadd6yBSKIivgkvSJ5MRV1KS04CY+jO7byFh2cT9XQbGImdGCULlUyieLwpqacl43SNlDUzO1elZGnjE+dRRUuRq0WpL37UEVE43Lw4fcs1eywacLF5xrgkwG+n+XEymsSWj/6XtkavMf4ZTf+TmKsxZUaXatK4u9uAVBsH6VMlnkVwHmRy9hnFyXkJzJysgzrIw8g4ujDc/38Ce4iQ8GbRbJu//kbkQkitQk7thXZa9vNy471chOEo/ISRAFVcDN6nvmG5edjQyDRK6Z4flV2KbYwrQs9uKubB4+fMjChfM4d+4sCQkPWb8+Al/f6gWev3Ll1+zevYvr168xatRYxo172Xhs7drVfP/9GuPPer0BrTaLyMituLu7069fP27fvm08npWVRdu27VmwYNE/5+tZtWo50dERpKen4+fnx5dfLsfV1TVXDJMnv8LRo4fZtesANjbZ1cLFi+dZtGghly9fxMnJmfDwZxgzZrzxPdu3b2X16uXEx8fj7e3NSy9NpFOnLsbj//vfj/z441oyMzPo3LkbU6e+azw2adLLXL16mawsLb6+1XnxxZfp2DH7vfv27eH779dw9epl7Ozs6dChI5MmvYmTk7Px/YcOHeTrr7/kxo3ruLlVYeLEN+jevSc3blxn6dL/curUCfR6A4GBjXnjjanUqlWn0Ljs7OyMx3/5ZR3r169DrVZRvXp15s5dSK1atYv8LpYu/S/btm0hNTUVV1dXwsOfYfTocQV+79amUiaL/CbllaQNPlWj4/uokzgc3Y3b33vRJyVy17Eae6r34Jqjb54kkSPnSaCgCrigjm1XJzsGda5f4CRCUxP7aJeMXC6jbdv2jBw5hldeGVvk+X5+NZkwYTJ//PFbnmOjRo1l1Kh/y1i1ajnHjx/D3T1798Po6GjjXCFJknjuuYF07do91/mnTp1g+fI1eHv7/FMB2+W6xpYtMej1eZ9mZ89+n44du/DVV8u5d+8ur746joYN/QkJ6cyDB/F8+OEMPv74M9q1a8/+/XuZMeMdfv01Eg8PJQcP7ueHH77jv//9mqpVvXjvvamsWrWcmTPfA+D116dSp05dbGxsOH36FG+88Srr1m2gatWqpKamMnr0OFq0aElWVhazZ7/PkiX/Zdq07PdevXqF2bPf5//+7wPatGlLWloqqampAKSmphAS0on33puFk5Mza9asZPr0t/jpp+zfbUFxTZgwCYDIyI1ER//BggVfUKdOXTIyEtFq5cX6LsLCBjBmzEs4Ojry4EE8b775GnXq1KVz525F/g1Yg0qZLCDvpLxpS/cWK2HYGbQEJZ3n6cQzOF/MwK5RIOurdeCUQVlgkoC8TwL5VcCFjWgqaBKhUPa+//5bIiN/R61W4+3tzfjxr9K5c9div1+p9GTQoMHodIXvbZKjb98wALZu3VzoeZIkERu7Kdfd/aP+/vsoiYkqunTJThbJycmsX7+Ob79dh4+PLwD16jXI9Z7U1FRWr17J++/P5pVXxuQ6dvfuHXr16otCoaBGDT+eeqoFV69eISSkM/Hx8bi4uBIcnD2hs337EBwdHbl9+xYeHkpiYqIICxtAvXrZf/MvvPAic+bMMCaLBg0aGq8jk4FeryM+/h5Vq1alV68+xmMODg6Ehw9k1arlxte++24VAwYMMl67ShV3qlTJrrAbN25K48ZNjecOGTKM775bRVJSIlWquBcY14QJkzAYDKxZs5L33ptF3br1AKhVq1a+/8/l9108+vQC2TcNt27dzOebsk6VNlk8Ps+iWX1P9p68V2DTlL0+i1ZJ52iTeAZHQxZXHH3Zp2zGnKnPc+qTHdkr+hWguE8CYplty1Cjhh9Ll36DUunJzp3b+PDDGTRpspHbt2/xzjtvFPi++fO/oHnzFiaL6/jxY6jVqgLvVGNioujSpTuOjo4AXLlyCYVCwc6d2/jll3U4OTkzePBQ/vOf54zvWb58Cc888x88H9k4K8dzzw1j8+Zoxo+fwJ07tzh9+iTDh48CoFGjQOrUqcuePXEEB4ewd+9ubG3tqF8/OwlcvXqFjh07G8tq0MAflSoBtVpNTrXz9ttvcPjwX2RlZfH008E0atQ438/199/HjJU3wOnTp6hRw49Ro4aQlJRIq1ZP88YbU3Fzq5Lvez09PY3JpKC4kpIS0WgyiI+/z5Url/noo9koFAoGDXqGIUNGI5fnHjxS0Hfx/fffsnbtKjQaDb6+NejZsw8VRaVMFvnNs9h78h4dnvLhxOWEXBW2gz6T1olnaZ10FgeDlktOfuxVNuOuQ9VcI4sKquRLspSG6Ey2DN269TD+u3v3Xnz//becPXuKjh27sHnzLrPFlZMMnJyc8hzLyMhg164dfPLJZ8bX4uPjSU1N5ebNG6xf/wc3b97kjTcmUKtWLdq0ace5c2c4efI4r7/+Fg8exOcps337EObOncXPP/+AXq9nzJjxBAY2AUChUNCnTyizZ79PVlYWNjY2fPjhfGOi0mjScXb+d6E6F5fsf6elpWFvn12pL1jwBTqdjkOHDnLjxrU8FTLAoUMHiImJYsWKb42vPXhwn9jYTXz++WKqVvVi3rxZLFq0kFmz5uZ6b3z8fT7/fD6vvfam8bWC4kpPT+fhwwfGa3733c+kpqYwbdpknJ3dCQ9/pljfxciRLzBixGguXjzP7t1xxvIrgko5dLageRYnLiew8NUOrJ7ejepOEp0SjjLh2gZC1Ce47ujLGr9+/Fq9G3cdqmKjkBkr8YL2pGhWP+/dWmFKO7xVKFsxMVG88MIw+vTpQp8+Xbh69TJJSYlmjSkzM4OdO7cbm6weFxe3Azc3N4KCWhlfy9mmeMyY8djbO9CgQUO6d+/F/v17MRgMfPbZJ7z++lvGDu1HJScn8dZbk3nhhRfZvn0vGzZEc/DgfjZsWA9kdzAvXfoVX321nJ0797N48Qrmz/+QixfPA+Do6ERaWpqxvLS07D4FZ2fnXNexsbEhOLgDBw8eYM+euFzHTp06yezZ7zN37nxq1aqd63OFhvanVq3aODk5MXLkWPbv35vrvWq1milTXuOZZ57NdXdfUFxOTk7G39ewYaNwdXXF17c6Q4YMyVN2Ud+FTCbD378Rdnb2uZrPrF2lfLIobNisLikJ9ZYYRp7ZDjot51xqs8+jGQ/s/92OTiaDjs18jZV4cBMfLt1KZOexO7nK23vyHg383EtU2YvOZPO6d+8uCxbM44svvqZp06dQKBS88MIwJCm76WHq1MkFvvfTT7+kefMgk8QVF7czTzJ4VExMFH369EP2SL/Zv/0CedtI09LSOHfuLLNmZfchGAzZHdyDBvXjww8/wd7eHrlcbqwQq1XzpkePXhw4sJdBgwZz8eIFmjcPMjYdBQY2oXHjphw69BcNGwZQt249Ll26QPfuPQG4dOkiSqUnHh4e+fYB6PV6bt++Zfz5woVzTJ8+henTZ9K69dO5zq1fv2Guz/m45ORkpkyZSIcOnfKMRiooripV3LG3d8DW1rbQsqHo7+Lfz6TL9ZmsXaV8ssivD8BFl05Y0hGuTp+Kekssbq1akzRqCnv8e+VKFACSlJ0I9p++Z3ztxOWEPGU+6U52QvnTaDTIZDI8PHJGG0Vw9Wr2d9i8eRBbt+4u8L9HE0VmZiZarRYArTaLzMyCB0/odDoyMzMxGAzo9XoyMzPzjE7KLxnkiI+/z7FjR/Lc6dao4Ufz5kGsXbuKrKwsrl27yvbtW2nfviMuLi5s3BjDmjU/smbNjyxc+F8AVq36nsaNm1KzZi1AYsuWzRgMBhISHrJ9+1Zjn0RgYGNOnDhmfJK4cOEcx4//bUxQffr0++d3d4Xk5GS++26VMb7r16+xf/9eMjMz0Ol0xMZu4vjxo7RokV35Xrlyibfemsybb04jJKRTns8bGtqfTZsiuX37FhkZGfz443d06BACZD8pvPXWazz1VHPjCKdHFRaXg4MD3br15Kef1pKenkZ8/H1++eUXY9mFfRcGg4GNG38jOTkZSZI4c+YUv/++nlatcic6a2Zxmx+VxpNsfuSmTaWd+hTNki8hl4F7+w4oQ/th5/3v3X1BI6Ue7ZMY+0nBmzWtnl42Q+cq6kYsliIn3uXLl7Bx42/I5TL69OnH+fPn6N07lP79Bxa7rJCQ1nle27PnMAALF34EYBwKOm/eB8TEROU69733ZhEa2h+ABw/iefbZ/vz446/4+eXe+c/Ly5XPP/+S/fv3snTpN3mu+eBBPB9//CEnT/6Nu7uS4cNHMXDgf/Kcd/fuHQYPDs81z+LIkUN8/fVX3Lx5HXt7Bzp06Mjrr0/FwcEBgN9++x+//LIOlUqFu7sHgwYN5vnnRxjL/PnnH/6Zz5BJly7Z8xlq1PDk0KETzJv3AdeuXUWhkOPnV4uRI8cYR5x99NFsYmKijNcB8Pb25YcffjH+vGrVcn7/PbtJrG3bYF5/fRpubm7ExEQxb94HODg45KrMv/9+PT4+PgXGlTOcOC0tlQUL5rFv315cXV0ZMuQ5nntulLGsgr4Lg8HA1Kmvc/bsaXQ6LZ6eXoSGhjFy5Jgin1TKkik3P6qUyQLg542Hkf+5maeSLwMyTrjV54BHU5JsXfOMXipOIihOQnlS1lr5WgtrixesL2YRr2mZMllUymYoAPmh3TRNucLfVfxZVvsZYqsFk2SbPas1ITmTNf4iZigAABjASURBVJvOGpuZnB0UBZaTc05RazMJgiBYs0rZwQ2w1aUpO5wDyZLb5Xtcp5dYt+0CAJnagpcFWR2VPZFOLIshCEJFVmmThXsVpyJnbKdqdGyIu5xrTabH6aV/FwgUI5kEQaioKm2yGNS5Pms3ny9yhdfiLAHy6DliC1JBECqiSpssgpv44ObqwLdRpwtMCM4OChzsbIpMGHLZv30XYgtSQRAqokqbLAC6tKpJk1ru7D99j9VRZ3i0tUkhg2E9A4Ci98M2SNnn2NrIxBakgiBUSJU6WeR4vHPa2UGBTv/vKrB2NjJcHG1I1eiQybIn5T0uS2cgq4BFRovTlCWarwRBsGSVdujs44Kb+LDw1Q6M79+YTK0hV19Glk5Ck6FjfP/GrHqn5BPsilo1NmeSYE5SyWm+enSGuFA+nn22P4cOHTR3GIJgccSTxWMKGv2kl2DVP8NknR0UpGXk7Ri3t5UhSbISrxpb0MKGovmqYpgzZwZHjvyFRpOBUunJ8OGjCpwNvnDhR2zZEmP8WafTYWNjy9atf5KVlcVnn33C4cN/kZycjJ+fH9OmTaVx45bG8zMyMli8+At27tyKTqejQQN/lixZCWTPel67dnWuzY++/XYdNWr4kZiYyLvvTuH69esYDHpq167LxImv06xZiyLjguyFBz/++EMOHTpAlSruvPzya8Z9KSRJYu3a1fzxxwbS0lJp1649b7/9f8bVXx88iOezzz7h+PG/cXBwYPTosQwc+CyQ/3pcGo2GuXPn06VLdyRJYuXKr9m0KRKNJp2GDQOYMuUd434VyclJfPrpJxw58hcATz8dzNSp043XLmzHPoDIyEgWLvyUpKRE2rRpy7vvzjQuhV7c73X16hWsXr2CRYuW0KZN23y/d2tQqZPFriM3jR3cOU0/hTUZ5fRNFDR939ZGwfM9/Plp63ljMrGzLfrhraj9wAXrNmLEC0yfPgM7OzuuX7/GpEkv07BhAI0aBeY5d9q094zLgED2UiA5f296vZ5q1bxZvHgF3t4+7N+/lzfeeIPvvvvZuG3rggXz0Ot1/PDDr7i5uXHx4oVc5Xfv3ouZMz/Mc11HR0fefXcmfn61kMlk7N4dxzvvTCEycgs2NjaFxgXw2WfzsbW1JSJiCxcvXuDtt1+nQYOG1KtXn82bo4mN3cTXX6+ibt3qTJ78BosWLeT992cD2ZVugwb+zJ27gKtXrzB58ivUqlWHli1bG9fjynH06GHeeWcKbdu2B2DHjm1ER0ewdOk3+Pj4snLl18ydO5PVq38EYMWKr0lJSeaXX/5AkuD//u9tVq9ewaRJU4DCd+y7cuUyM2fOZMGCL/D3b8SCBfP47LNPmD3742J/r7dv32LXru14elYt8u/E0lXaZqj9p++xeP3xXE0/3xSwU92jsnSGAofb5uzXrdVJuV4rqkmpoGYqsemReZw7d4YRIwbTp09XPvpodqGLABZHvXr1jXfzMln2f8VZjVSj0bBr1w7jQneOjo6MG/cyvr7VkcvldOjQET8/P86fPwvAjRvX2LPnT95++//w8PBAoVDkm5DyY29vT61adZDL5UiShFwuJyUlmZSU5CLj0mg0xMXt4MUXX8HJyYnmzVsQEtKJ2NhNAOzd+ydhYQPw9vbB2dmZ4cNHs2PHVjIyMkhPT+fYsSOMGjUWGxsbGjb0p0uX7kRH/5FvnJs3R9O1678bPN29e5tmzVpQo4YfCoWCXr36cu3aVeP5d+/eoWPHLjg7u+Dy/+3de1xUdfrA8Q/DLQiN+0XUvGwmoWXqZm0ZSSCosHhZLz8VNS1d3c0yUUZTES8vxUtZZtv6Wq1Md+vlmhqgaHhJt620vCRgmxooAaIC5oWLMHN+f0xMDjPDgMIMwvP+b871mXO+zMP5nnO+j5sbzz77HNnZP+nn/+53D+nHwrq9Yh/oKheGhobSo0dPXF1deemlqXzxxQFKS3VDnNflvL7xxgqmTn0ZR0fHOp2HpqzFXll8+sU5ox/9hhgka0NKFjWHqbLUpSRFj5qWvXt3s3r1WlxcXIiPn8GHH27gj38cyoQJo8yu89praoNyoDWtWrWc3buTqaiooEuXh/UlQWtz8OA+3N096NGjp8n5xcVF5OTk0LGjrp1kZmbg7+/Phg1/Z8+eXXh5eTNx4mR9mVXQ/XAPGBCKl5c3w4aNYMiQPxlsc/z4UZw/n0NVVRXR0YPx8PC0GFdu7nlUKnuDmhOdO3fhxIljgK4b6vYh6BRF4datW+TmXiAwMLB66m17UPjpJ+PRmsvLdXUkkpLe0E8LC4tg//7PuXDhPG3aBJKWlkKfPk/p5w8dOpzt27cSFhYB6Op+1BzJ1lzFvuzsc/Tp83v9coGBbXF0dOTChQv6JFzbed2/Px1HRweeeuoZIMno+9xrWmyyaKwuHnPjGda2PxkqpGkZNmwEfr+OOjxu3ETefHMlkydPu6sqeXFxambMmEVGximOH//W4L6BOWlpqURGDjTZ7VlVVUVi4nyGDBnCgw92AHR9/z/9dI6QkFB27EgjI+N7Zs9+lQ4dOtGhQ0dCQ8OJiRmKh4cnWVkZzJs3Gzc3N4PiQB9++DEVFRUcOnTAbA3xmnGVlZUZVYRzc3PT/wf+5JNPs2XLJkJDw7nvvjZs2fIhoCsi5Op6P927P8YHH/yDadNeIScnm4MH9+Pu7m60X12ScjeoI+Hl5c2jjz7O6NHDsLe3x9fXj7fe+pt+/sMPd6WysopBg3QJs1ev3zNkyHCD7Zqr2FdWVkarVq0Mlr3//t++F5g/r6Wlpaxfv4433njH5DG8FzWJbqjs7GxGjhxJREQEI0eOJCcnp9H3ae0uHkv7q34aa6M6lJXTnm6QRPFV5kVmvfslE5fvZ9a7X8rTVXXk6/vbsffzC+DKlSt1XnfmzOmEh/clPLyvwQ1h0JUifeyxHly+fInt2/9d63YKCy9y4sQxIiMHGc3TarUsXjwfR0cH5s+fr5/u7OyMg4MD48dPwtHRkccf78Xjj/fmyJGvAV3hH29vH+zt7ene/TH+9Kf/4+DBfUbbd3Z2Jjw8ks2bPzC652EqLhcXF33FuWo3b97E1VVXFW/QoD8SHh7Byy9PYdCgQfTsqRu+3cfHF4CEhCUUFOQzbNggVq1aRv/+kfj6+hnFtXt3qlEdiY0b1/PDD1l8+mkq+/Z9yQsvvMQrr0ylvLwcgPnz42nXrj179x5iz54vCAxsy6JF8422bapin4uLCzdumP9e1Uyd1w0b/k5ExEDatAmkuWgSVxYJCQmMHj2amJgYdu7cyYIFC9i0aVOj7nNoSGf9exR3ylzt7Zps0aVkqs64vE1eN9V91qD7cfT29ubixYvExg43u86sWXPp338Aq1e/bXH7NavCmZKWlkq3bo8SGNjWYLqiKCxfvpji4mJWrXrr175w3Q9jdWGiujL3zlC1qqoq8vPzeOihLrXG1a7dg2g0GnJzL/xaNAnOnv2Rjh07AaBSqZg0aQqTJk3Bx6cVqamf4+Pjq08W/v4BrFixRr+9hQtf19f6rqZLUt8xe/Zcg+nnzp0hNDRcn1wGDozm7bdXk5PzE127PsLZs2d47TW1/h7H4MHDmDbtRbPf+fZz07FjZ3744Qf9vLy8n6msvEX79u0trvvdd0e5fLlQnzyuXi1hwYI5jBkzjrFjJ5jdf1Nm8yuLoqIisrKyiIrS3SyLiooiKyuL4uLiRt3v3f5gVicAc1cMql//+bFVHe3aHscVtfv0061culTItWu/8NFH7/P88/3x9/evtUpe//4DTG6rpKSY9PQ9lJaWotFo+Oabr0hP30OvXsbFkW6XlpZqssbzqlXLyMnJJinpTZyd7zOY16NHT/z8/Nm8+QOqqqr4/vsTHD/+nb4P//DhgwaV3P7970/o2zcE0NW7PnnyBJWVlVRUlLN58wcUFxcTHNzNYlwuLi6EhPTjH/94j7KyMr7//gT/+c8XREQMBHSPr+bl/YyiKJw9e5a1a99gwoQX9d09OTnZlJbepLKykj17dnH06NeMHDnGYB979uwymTy7dn2EAwfSKS4uQqvVkpaWSlVVFYGB7fTzU1J2UFFRTkVFOTt3btdX87NUsS88PJIDBw5w8uRxysrK2LDh74SE9MPV9X6L5/Wtt95l06ZP9JUIvb19mDVrLkOHjqj1vDdlNr+yKCgowM/PD3t7Xc0IXb+jLwUFBXh6Gt9ca0g+Hi5cLikzmq6yM33voXp6zXsKpm5O2yJB3E4ex71z4eGRvPbaX7ly5TLPPBNiVMe5fuzYsWMbq1YtQ6tV8Pf3Z/r0mfpn+auvWG6v5JaR8T2XL18iNDTMYEsXLxawc+enODk5EROju2FrZ2dHXNwc+vcfgIODA8uWrSYpaQmbN3+Av38A8+Yl6u9ppKfvZdmyxVRW3sLHx5cxY8brf/grK2+xZs0q8vPzcHBwoHPn37Fy5Rq8vX30+zcXF8DMmWqWLVtEdHQ4rVs/wMyZc/TvOly9epX4+BlculSIl5cXQ4eOICZmqH7dI0e+YtOmjZSXl/PQQw+zevVaPDwMSxmnpaUyenSs0X7HjBlPSUkJL7wwmrKyctq2bcuSJSv09xrmzFnAmjWrGDJkEIqiEBQUzNy5CwHdVdrGjetZsOC3in2Jict4+OGugO5pp8TERBIT53Ht2i/07v0Ec+cm1Om8PvCA4T0XlUpFq1atcHV1NfoO9wqbV8rLyMggPj6e1NRU/bSBAweycuVKgoODa1nz7h38Lpd3tp40eCrK2dGe53u3Zd+3PxtN/+vwx3iuVzuT29m0+zRXSsrw9nBh3IAgk8tZ08Qle00mQh8PFzbO62+DiIQQ9zKbX1kEBARQWFiIRqPB3t4ejUbDpUuXCAgIqPM27qSsKugGErx2vdzkU0iB3vcbTQ9u726yZGFwe3eSpjxlMK0xSjHWp8Tj4Gc6mrziGfxMR6uViWyuJSmbknstZom3cTVmWVWbJwsvLy+CgoJISUkhJiaGlJQUgoKCGr0LypJ7vZCRPI4rhGhINk8WAAsXLkStVvPuu+/SunVrkpKs8wLLwe9ym/UTQ/d6whNCNB1NIll07tyZrVu3Wn2/m3aflgH8hBCiDmz+6KwtXTFxAxjkiSEhhKipRScLbw8Xk9PdXJrEBZcQQjQZLfpXcdyAIN765LhR/Yqy8iq+yrx4x11RUvVOCNHctOgri+d6tcPZRL0JjcIdv+ksVe+EEM1Ri04WgMmKd3Dn9y1kmA0hRHPU4pNFQxcekmE2hBDNUYtPFkNDOuPkYHgY7maUWKl6J4Rojlp8sngq2J/xA7rqf8zvdpTYhk4+QgjRFLTop6GqNeSbzjLMhhCiOZJk0QhkmA0hRHPT4ruhhBBCWNYiryyqX5orvlaBp3QTCSGERS0uWUhtaiGEqL8W1w0lL80JIUT9tbhkIS/NCSFE/bW4ZCEvzQkhRP21uGQhL80JIUT9tbgb3Le/NCdPQwkhRN20uGQBv7005+PTisuXr9s6HCGEaPJaXDeUEEKI+pNkIYQQwiJJFkIIISySZCGEEMKiZnGDW6Wys8m6tiDxNq57LV6492KWeBtXXeK9k+9kpyiKcicBCSGEaDmkG0oIIYRFkiyEEEJYJMlCCCGERZIshBBCWCTJQgghhEWSLIQQQlgkyUIIIYRFkiyEEEJYJMlCCCGERc1iuI/ahIaG4uTkhLOzrmxqXFwcffv2NVimrKyMOXPmkJmZib29PfHx8fTr18/qsf7888/85S9/0X++fv06N27c4MiRIwbLrV27ln/+85/4+voC0LNnTxISEqwSY1JSEnv27CEvL4/k5GS6dOkCQHZ2Nmq1mqtXr+Lu7k5SUhIdOnQwWl+j0bBkyRIOHz6MnZ0dkydPZvjw4VaPuaSkhNmzZ3PhwgWcnJx48MEHWbRoEZ6enkbrq9Vq/vvf/+Lh4QFAZGQkU6dOtWq8ULe2DNZvz6birWtbBuu359rO/YkTJ1iwYAEVFRUEBgaycuVKvLy8jLZhzWNsLt5ffvmFBQsWcPnyZRwcHOjevTsJCQncd999RtuIjY0lPz8fNzc3AMaNG8ewYcPqF4jSzPXr10/53//+V+sya9euVebOnasoiqJkZ2crf/jDH5QbN25YI7xaLVmyRElMTDSa/vbbbyvLly+3QUSKcvToUSU/P9/ouMbGxio7duxQFEVRduzYocTGxppcf/v27crEiRMVjUajFBUVKX379lVyc3OtHnNJSYny9ddf65dZvny5MmfOHJPrx8fHKx999FGjxng7c8e4Lm1ZUazfns3FeztzbVlRrN+ezZ17rVarhIWFKUePHlUURVHWrVunqNVqk9uw5jE2F29ubq6SmZmpKIqiaDQa5ZVXXlHeeecdk9sYO3assn///ruKQ7qhgN27dzNq1CgAOnToQLdu3Th06JBNY7p16xbJycn1z/6NrHfv3gQEBBhMKyoqIisri6ioKACioqLIysqiuLjYaP1du3YxfPhwVCoVnp6ehIWFkZaWZvWY3d3d6dOnj/5zjx49yM/Pb9Q46spUvPVh7fZsKd6m1pbNnftTp07h7OxM7969ARg1apTZtmnNY2wu3rZt2/LII48AoFKpePTRRxu1DbeIZBEXF0d0dDQLFy7k2rVrRvPz8/MJDAzUfw4ICODixYvWDNHI/v378fPzIzg42OT81NRUoqOjmThxIsePH7dydIYKCgrw8/PD3t4eAHt7e3x9fSkoKDC5bJs2bfSfm8Kx1mq1/Otf/yI0NNTsMu+//z7R0dFMmzaNc+fOWTE6Q5baMjS99mypLYPt2vPt575m2/T09ESr1XL16lWj9Wx1jM211fLycrZt21ZrG16xYgXR0dHExcVRWFhY7303+2SxZcsWPvvsM7Zt24aiKCxatMjWIdXJtm3bzP4nNmrUKPbt20dycjKTJk1i2rRplJSUWDnC5mPx4sW4uroyduxYk/NnzJjB559/TnJyMv379+fFF19Eo9FYOcrm2ZbBtu3Z0rlvakzFW1VVxYwZM3jyySd5/vnnTa63YsUKdu/ezY4dO+jUqROvvvpqvffd7JNF9eWxk5MTo0eP5tixY0bLtGnThry8PP3ngoIC/P39rRZjTYWFhRw9epTo6GiT8318fHB0dATg6aefJiAggDNnzlgzRAMBAQEUFhbqf0A1Gg2XLl0y2TUREBBgcKls62OdlJTE+fPnWbNmDSqV6T8HPz8//bzBgwdTWlpqk//U69KWoWm1Z0ttGWzXnmue+5pts7i4GDs7O9zd3Y3WtcUxNtVWNRoNcXFxPPDAA8ybN8/sutVtx97ennHjxnHy5Em0Wm299t+sk0VpaSnXr18HQFEUdu3aRVBQkNFykZGRfPLJJwDk5ORw6tQpk0+ZWMv27dsJCQnRP31T0+2XkKdPnyYvL4+OHTtaKzwjXl5eBAUFkZKSAkBKSgpBQUEmnyyKjIxk69ataLVaiouLSU9PJyIiwtohA/Dmm2+SkZHBunXrcHJyMrvc7cf78OHDqFQq/Pz8rBGiXl3bMjSt9mypLYNt2rOpc9+tWzfKy8v59ttvAfj4448ZMGCAyfWtfYxNxavValGr1djb27N06VLs7EwXNKqqquLKlSv6z6mpqXTp0sXsP0fmNOviR7m5ubz88stoNBq0Wi2dO3dm3rx5+Pr6EhMTw/r16/Hz86O0tBS1Ws3p06dRqVTMmjWLsLAwm8UdERHB66+/zrPPPquf9tJLLzF9+nS6d+9OfHw8mZmZqFQqHB0dmT59OiEhIVaJbcmSJezdu5crV67g4eGBu7s7qampnDt3DrVazbVr12jdujVJSUl06tTJKHaNRsOiRYv48ssv9fNGjhxp9ZjXrFlDVFQUHTp00D9q2LZtW9atWwdg0D4mTJhAUVERdnZ2uLm5MXv2bHr06GHVeN977z2zbblmvNZuz+baBJhuy2Db9nzmzBmz5/7YsWMkJCQYPDrr7e0N2O4Ym4t3+PDhTJkyxeCHv/qx48LCQiZPnszOnTspLS1l7NixVFZWAuDr68vrr7+u//usq2adLIQQQjSMZt0NJYQQomFIshBCCGGRJAshhBAWSbIQQghhkSQLIYQQFkmyEEIIYZEkCyGEEBZJshBCCGGRJAshGtCFCxd44oknyMzMBHRDWfTp04dvvvnGxpEJcXckWQjRgNq3b09cXBxxcXGUlZUxd+5chg4dalCPQIh7kQz3IUQj+POf/6wflXTbtm21DlQoxL1AriyEaAQjRozgxx9/JDY2VhKFaBbkykKIBnbz5k1iYmLo06cPhw4dIjk52WRNBCHuJXJlIUQDW7p0KcHBwSxdupTnnnuOhIQEW4ckxF2TZCFEA0pPT+fw4cMkJiYCoFarycrK4rPPPrNxZELcHemGEkIIYZFcWQghhLBIkoUQQgiLJFkIIYSwSJKFEEIIiyRZCCGEsEiShRBCCIskWQghhLBIkoUQQgiLJFkIIYSw6P8Bob+m3Ojrqr0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "N=len(x)\n", "\n", "xbar= np.mean(x); ybar=np.mean(y)\n", "USum=0.; DSum=0.\n", "for k in range (0,N):\n", " USum+=(x[k]-xbar)*(y[k]-ybar)\n", " DSum+=(x[k]-xbar)**2\n", " \n", "a=USum/DSum; b=ybar-a*xbar\n", "\n", "xfit=np.linspace(min(x),max(x),100)\n", "yfit=a*xfit+b\n", "\n", "import seaborn as sns\n", "sns.set()\n", "plt.scatter(x,y,label='podaci')\n", "plt.plot(xfit,yfit,'r',label='$y=ax+b$')\n", "plt.xlabel('x'); plt.ylabel('y')\n", "plt.title('Metod najmanjih kvadrata')\n", "plt.legend()\n", "plt.annotate('a='+str(a),(12.5,2.5))\n", "plt.annotate('b='+str(b),(12.5,0))\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.0" } }, "nbformat": 4, "nbformat_minor": 2 }