{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np, matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[0. 1.67203]\n", " [0.1 1.79792]\n", " [0.2 2.37791]\n", " [0.3 2.66408]\n", " [0.4 2.11245]\n", " [0.5 2.43969]\n", " [0.6 1.88843]\n", " [0.7 1.59447]\n", " [0.8 1.79634]\n", " [0.9 1.0781 ]\n", " [1. 0.21066]]\n" ] } ], "source": [ "# Os dados são lidos dos arquivos \"t_data.txt\" (tempo) e \"h_data.txt\" (altura)\n", "t_data = np.loadtxt(\"t_data.txt\"); h_data = np.loadtxt(\"h_data.txt\"); # Leitura\n", "n = t_data.size; t_data = t_data.reshape(n,1); h_data = h_data.reshape(n,1); # Reshape\n", "print(np.hstack([t_data, h_data])) # Dados" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAXWklEQVR4nO3de7SddX3n8ffHGCUWNGpiJQGJjDSriGj0iLi81NZLlFKgrTNeoV5aVlE7tnXSAbXV0pnRrrSujqXK0MKyeEGrzcR4m9QLFXXEciBIuJjKKEgSKhEMoB5tgO/8sXfgEE7O2SfZz97Z+3m/1jorez/Ps/f+/nLOPp/z/H6//XtSVUiS2utBwy5AkjRcBoEktZxBIEktZxBIUssZBJLUcgaBJLWcQSA1IMkNSV4w7DqkXhgEEvf+4p5KcmeSnUn+b5LfTeJ7RGPPH3LpPr9WVYcARwDvBv4rcP5wS5KaZxBIe6iq26tqA/Ay4LeSHJPkV5NsSnJHkpuSvHP6Y5KcmuTGJLcmedse+x6a5K+SbO9+/VWSh3b3LUny6e5ZyG1JvuJZiAbNHzhpL6rqX4CtwHOAHwOnAYuBXwXOSHIKQJKjgfcDpwLLgEcDh017qrcBxwNPAZ4MHAe8vbvvLd3XWAr8PPBWwHVfNFAGgTS77cCjquqfq2pzVd1TVVcBFwG/1D3mpcCnq+qSqvoZ8MfAPdOe41XA2VV1S1XtAP6UTmgA7AIOBY6oql1V9ZVyATANmEEgzW45cFuSZyS5OMmOJLcDvwss6R6zDLhp9wOq6sfArdOeYxlw47T7N3a3AawFrgf+Kcl3kpzZUDukvTIIpL1I8nQ6QfBV4CPABuDwqnoEcC6Q7qE3A4dPe9zD6HQP7badzgD0bo/rbqOq7qyqt1TVkcBJwB8meX4zLZJmZhBIe0jy8CQnAh8FPlRVm4FDgNuq6qdJjgNeOe0hnwBOTPLsJA8Bzub+762LgLcnWZpkCfAnwIe6r3VikickCXA7cDf371aSGvfgYRcgHUA+leQuOr+IrwXeQ+cvf4A3AH+Z5Bzgy8A/0Bk4pqquSfJGOmcNP9d93NZpz/vfgIcDV3Xvf7y7DeAo4Bw6g8U/BN5XVRc30Thpb+K4lCS1m11DktRyBoEktZxBIEktZxBIUsuN3KyhJUuW1IoVK4ZdhiSNlMsvv/wHVbV0pn0jFwQrVqxgcnJy2GVI0khJcuPe9tk1JEktZxBIUssZBJLUcgaBJLWcQSBJLTdys4Z04Fu/aRtrN25h+84pli1exJrVKzll1fJhlyVpLwwC9dX6Tds4a91mpnbdDcC2nVOctW4zgGEgHaDsGlJfrd245d4Q2G1q192s3bhlSBVJmotBoL7avnNqXtslDZ9BoL5atnjRvLZLGj6DQH21ZvVKFi1ccL9tixYuYM3qlUOqSNJcHCxWX+0eEHbWkDQ6DAL13SmrlvuLXxohdg1JUssZBJLUcgaBJLWcQSBJLWcQSFLLNRYESQ5PcnGSa5Nck+TNMxzzvCS3J7my+/UnTdUjSZpZk9NH7wLeUlVXJDkEuDzJ56vq2j2O+0pVndhgHZKkWTR2RlBVN1fVFd3bdwLXAU4ul6QDzEDGCJKsAFYB35hh9zOTfDPJ55I8cS+PPz3JZJLJHTt2NFmqJLVO40GQ5GDgH4Hfr6o79th9BXBEVT0Z+Gtg/UzPUVXnVdVEVU0sXbq00XolqW0aXWIiyUI6IfDhqlq35/7pwVBVn03yviRLquoHTdYl9ZtXZdMoaywIkgQ4H7iuqt6zl2MeC3y/qirJcXTOUG5tqiapCV6VTaOuyTOCZwGnApuTXNnd9lbgcQBVdS7wUuCMJHcBU8DLq6oarEnqu9muymYQaBQ0FgRV9VUgcxxzDnBOUzVIg+BV2TTq/GSxtJ+8KptGnUEwptZv2saz3v0lHn/mZ3jWu7/E+k3bhl3S2PKqbBp1XphmDDl4OVhelU2jziAYQw5eDp5XZdMos2toDDl4KWk+DIIx5OClpPkwCMaQg5eS5sMxgjHk4KWk+TAIxpSDl5J6ZdeQJLWcQSBJLWcQSFLLGQSS1HIGgSS1nEEgSS1nEEhSyxkEktRyBoEktZxBIEktZxBIUssZBJLUcgaBJLWcQSBJLecy1Bob6zdt8xoM0j4wCDQW1m/axlnrNjO1624Atu2c4qx1mwEMA2kOdg1pLKzduOXeENhtatfdrN24ZUgVSaPDINBY2L5zal7bJd3HINBYWLZ40by2S7qPQaCxsGb1ShYtXHC/bYsWLmDN6pVDqkgaHQ4WayzsHhB21pA0fwaBxsYpq5b7i1/aB3YNSVLLNRYESQ5PcnGSa5Nck+TNMxyTJO9Ncn2Sq5I8tal6JEkza7Jr6C7gLVV1RZJDgMuTfL6qrp12zEuAo7pfzwDe3/1XkjQgjZ0RVNXNVXVF9/adwHXAnh24JwMXVselwOIkhzZVkyTpgQYyRpBkBbAK+MYeu5YDN027v5UHhgVJTk8ymWRyx44djdUpSW3UeBAkORj4R+D3q+qOfXmOqjqvqiaqamLp0qX9LVCSWq7RIEiykE4IfLiq1s1wyDbg8Gn3D+tukyQNSJOzhgKcD1xXVe/Zy2EbgNO6s4eOB26vqpubqkmS9EBNzhp6FnAqsDnJld1tbwUeB1BV5wKfBU4Argd+Ary2wXokSTNoLAiq6qtA5jimgDc2VYMkaW4uMSGNMK/Kpn4wCKQR5VXZ1C+uNSSNKK/Kpn4xCKQR5VXZ1C8GgTSivCqb+sUgkEaUV2VTvzhYLI0or8qmfjEIGub0PjWpjVdl8z3VfwZBg5zeJ/WX76lmOEbQIKf3Sf3le6oZBkGDnN4n9ZfvqWYYBA1yep/UX76nmmEQNMjpfVJ/+Z5qhoPFDXJ6n9Rfvqeakc5K0KNjYmKiJicnh12GJI2UJJdX1cRM++wakqSW67lrKMkjgWXAFHBDVd3TWFWSpIGZNQiSPILOFcReATwE2AEcBPx8kkuB91XVxY1XuZ/8JKIk7d1cZwSfAC4EnlNVO6fvSPI04NQkR1bV+Q3Vt9/8JKIkzW7WIKiqF86y73Lg8r5X1GezfRLRIJCk+Y0RHAusmP6YqlrXQE195ScRJWl2PQVBkguAY4FrgN2DxAUc8EGwbPEits3wS99PIkpSR69nBMdX1dGNVtKQNatX3m+MAPwkoiRN12sQfD3J0VV1baPVNMBPIkr950y88dJrEFxIJwz+DfgZEKCq6tjGKuujNl68Q2qKM/HGT69BcD5wKrCZ+8YIJLWQM/HGT69BsKOqNjRaiaSR4Ey88dNrEGxK8hHgU3S6hoDRmD4qqb+ciTd+el10bhGdAHgR8GvdrxObKkrSgctrAoyfns4Iquq1TRciaTQ4E2/8zLXo3NvpLCx32172/wrwsKr6dBPFSTowORNvvMx1RrAZ+FSSnwJXcN/qo0cBTwG+APyPmR7Y/TTyicAtVXXMDPufB3wS+G5307qqOnveLZAk7Ze5Fp37JPDJJEcBzwIOBe4APgScXlWzTRP4AHAOnc8g7M1XqsqxBkkaol7HCL4NfHs+T1xVlyRZsS9FSZIGZ9iXqnxmkm8m+VySJ+7toCSnJ5lMMrljx45B1idJY2+YQXAFcERVPRn4a2D93g6sqvOqaqKqJpYuXTqo+iSpFYYWBFV1R1X9qHv7s8DCJEuGVY8ktVWv1yM4CHg98EQ6s4YAqKrX7esLJ3ks8P2qqiTH0QmlW/f1+SRJ+6bXJSY+CHwLWA2cDbwKuG62ByS5CHgesCTJVuAdwEKAqjoXeClwRpK7gCng5VVV+9AGSdJ+SC+/e5NsqqpVSa6qqmOTLKQz9fP45ku8v4mJiZqcnBz0y0rSSEtyeVVNzLSv1zGCXd1/dyY5BngE8Jh+FCdJGq5eu4bOS/JI4O3ABuBg4I8bq0qSNDBzBkGSBwF3VNUPgUuAIxuvSpI0MHN2DVXVPcAfDaAWSdIQ9DpG8IUk/yXJ4Uketfur0cokSQPR6xjBy7r/vnHatsJuIkkaeb0uOvf4pguRJA1Hr58sPm2m7VU12xLTkqQR0GvX0NOn3T4IeD6dReMMAkkacb12Df3e9PtJFgMfbaIgSdJg7evqoz8GHDeQpDHQ6xjBp+jMEoJOeBwNfLypoiRJg9PrGMFfTLt9F3BjVW1toB5J0oD12jV0QlV9ufv1taramuTPG61MkjQQvQbBC2fY9pJ+FiJJGo5Zu4aSnAG8ATgyyVXTdh0CfK3JwiRJgzHXGMFHgM8B7wLOnLb9zqq6rbGqJEkDM1cQLADu4P5rDAGQ5FGGgSSNvrmC4HLumzaaPfa56JwkjYFZg8DF5iRp/PX6OQK6l6o8is5aQwBU1SVNFCVJGpxeP1n828CbgcOAK4Hjga8Dv9JYZZKkgej1cwRvprMC6Y1V9cvAKmBnU0VJkgan1yD4aVX9FCDJQ6vqW8DK5sqSJA1Kr2MEW7tLT68HPp/kh8CNTRUlSRqcXq9H8Ovdm+9McjHwCOD/NFaVJB1g1m/axtqNW9i+c4plixexZvVKTlm1fNhl9UXPs4Z2q6ovN1GIJB2o1m/axlnrNjO1624Atu2c4qx1mwHGIgz29cI0ktQaazduuTcEdpvadTdrN24ZUkX9ZRBI0hy275ya1/ZRYxBI0hyWLV40r+2jxiCQpDmsWb2SRQsX3G/booULWLN6PGbRNxYESS5IckuSq/eyP0nem+T6JFcleWpTtUjS/jhl1XLe9RtPYvniRQRYvngR7/qNJ43FQDHsw6yhefgAcA5w4V72v4TO2kVHAc8A3t/9V5IOOKesWj42v/j31NgZQXdButmuV3AycGF1XAosTnJoU/VIkmY2zDGC5cBN0+5v7W57gCSnJ5lMMrljx46BFCdJbTESg8VVdV5VTVTVxNKlS4ddjiSNlWEGwTbg8Gn3D+tukyQN0DCDYANwWnf20PHA7VV18xDrkaRWamzWUJKLgOcBS5JsBd4BLASoqnOBzwInANcDPwFe21QtkqS9aywIquoVc+wv4I1Nvb4kqTcjMVgsSWqOQSBJLWcQSFLLGQSS1HIGgSS1nEEgSS1nEEhSyxkEktRyBoEktZxBIEktZxBIUssZBJLUcgaBJLWcQSBJLWcQSFLLGQSS1HIGgSS1nEEgSS1nEEhSyxkEktRyBoEktZxBIEktZxBIUssZBJLUcgaBJLWcQSBJLWcQSFLLGQSS1HIGgSS1nEEgSS1nEEhSyzUaBElenGRLkuuTnDnD/tck2ZHkyu7XbzdZjyTpgR7c1BMnWQD8DfBCYCtwWZINVXXtHod+rKre1FQdkqTZNXlGcBxwfVV9p6r+HfgocHKDrydJ2gdNBsFy4KZp97d2t+3pN5NcleQTSQ5vsB5J0gyGPVj8KWBFVR0LfB74+5kOSnJ6kskkkzt27BhogZI07poMgm3A9L/wD+tuu1dV3VpVP+ve/TvgaTM9UVWdV1UTVTWxdOnSRoqVpLZqMgguA45K8vgkDwFeDmyYfkCSQ6fdPQm4rsF6JEkzaGzWUFXdleRNwEZgAXBBVV2T5Gxgsqo2AP85yUnAXcBtwGuaqkeSNLNU1bBrmJeJiYmanJwcdhmSNFKSXF5VEzPtG/ZgsSRpyAwCSWo5g0CSWs4gkKSWMwgkqeUMAklqOYNAklrOIJCkljMIJKnlDAJJajmDQJJaziCQpJYzCCSp5QwCSWo5g0CSWq6xC9NIkvpj/aZtrN24he07p1i2eBFrVq/klFXL+/b8BoEkHcDWb9rGWes2M7XrbgC27ZzirHWbAfoWBnYNSdIBbO3GLfeGwG5Tu+5m7cYtfXsNg0CSDmDbd07Na/u+MAgk6QC2bPGieW3fFwaBJB3A1qxeyaKFC+63bdHCBaxZvbJvr+FgsSQdwHYPCDtrSJJa7JRVy/v6i39Pdg1JUssZBJLUcgaBJLWcQSBJLWcQSFLLpaqGXcO8JNkB3LiPD18C/KCP5YwC29wOtrkd9qfNR1TV0pl2jFwQ7I8kk1U1Mew6Bsk2t4Ntboem2mzXkCS1nEEgSS3XtiA4b9gFDIFtbgfb3A6NtLlVYwSSpAdq2xmBJGkPBoEktdxYBkGSFyfZkuT6JGfOsP+hST7W3f+NJCuGUGZf9dDmP0xybZKrknwxyRHDqLOf5mrztON+M0klGfmphr20Ocl/6n6vr0nykUHX2G89/Gw/LsnFSTZ1f75PGEad/ZLkgiS3JLl6L/uT5L3d/4+rkjx1v1+0qsbqC1gA/D/gSOAhwDeBo/c45g3Aud3bLwc+Nuy6B9DmXwYe1r19Rhva3D3uEOAS4FJgYth1D+D7fBSwCXhk9/5jhl33ANp8HnBG9/bRwA3Drns/2/xc4KnA1XvZfwLwOSDA8cA39vc1x/GM4Djg+qr6TlX9O/BR4OQ9jjkZ+Pvu7U8Az0+SAdbYb3O2uaourqqfdO9eChw24Br7rZfvM8CfAX8O/HSQxTWklzb/DvA3VfVDgKq6ZcA19lsvbS7g4d3bjwC2D7C+vquqS4DbZjnkZODC6rgUWJzk0P15zXEMguXATdPub+1um/GYqroLuB149ECqa0YvbZ7u9XT+ohhlc7a5e8p8eFV9ZpCFNaiX7/MvAL+Q5GtJLk3y4oFV14xe2vxO4NVJtgKfBX5vMKUNzXzf73PyCmUtk+TVwATwS8OupUlJHgS8B3jNkEsZtAfT6R56Hp2zvkuSPKmqdg6zqIa9AvhAVf1lkmcCH0xyTFXdM+zCRsU4nhFsAw6fdv+w7rYZj0nyYDqnk7cOpLpm9NJmkrwAeBtwUlX9bEC1NWWuNh8CHAP8c5Ib6PSlbhjxAeNevs9bgQ1Vtauqvgv8K51gGFW9tPn1wD8AVNXXgYPoLM42rnp6v8/HOAbBZcBRSR6f5CF0BoM37HHMBuC3urdfCnypuqMwI2rONidZBfwvOiEw6v3GMEebq+r2qlpSVSuqagWdcZGTqmpyOOX2RS8/2+vpnA2QZAmdrqLvDLDGfuulzd8Dng+Q5BfpBMGOgVY5WBuA07qzh44Hbq+qm/fnCceua6iq7kryJmAjnRkHF1TVNUnOBiaragNwPp3Tx+vpDMq8fHgV778e27wWOBj4eHdc/HtVddLQit5PPbZ5rPTY5o3Ai5JcC9wNrKmqkT3b7bHNbwH+Nskf0Bk4fs0o/2GX5CI6Yb6kO+7xDmAhQFWdS2cc5ATgeuAnwGv3+zVH+P9LktQH49g1JEmaB4NAklrOIJCkljMIJKnlDAJJajmDQK2SZHGSNwy7jt2SLEry5SQLZjnmC0keOci61C4GgdpmMZ3VZw8UrwPWVdXdsxzzQQ6smjVmDAK1zbuB/5DkyiRrAZKsSXJZd233P+1uW5HkW0k+kORfk3w4yQu6i7l9O8lx3ePemeSDSb7e3f473e1JsjbJ1Uk2J3nZXup5FfDJ7mMOTXJJt7arkzyne8wGOuvpSI0Yu08WS3M4Ezimqp4CkORFdNbiOY7O+u4bkjyXzrIFTwD+I52/2i8DXgk8GzgJeCtwSvc5j6WzltHPAZuSfAZ4JvAU4Ml01r25LMkl05cC6C6ZcGRV3dDd9EpgY1X9925X0cMAquqH6VxM6dGj/ClhHbgMArXdi7pfm7r3D6YTDN8DvltVmwGSXAN8saoqyWZgxbTn+GRVTQFTSS6mEyrPBi7qdvl8P8mXgadz/3VylgA7p92/DLggyUJgfVVdOW3fLcAyRntxRB2g7BpS2wV4V1U9pfv1hKo6v7tv+gqt90y7fw/3/yNqz3Vael23ZYrOAmmdB3UuSPJcOitJfiDJadOOPah7vNR3BoHa5k46S1TvthF4XZKDAZIsT/KYeT7nyUkOSvJoOouFXQZ8BXhZkgVJltL5Bf8v0x/UvYrYgiQHdV/7COD7VfW3wN/RuVwh3avnPRa4YZ51ST2xa0itUlW3dgd8rwY+V1VruksXf727KuuPgFfTWbmzV1cBF9Pp6vmzqtqe5H/TGSf4Jp0zhD+qqn+b4bH/RKcb6Qt0QmRNkl3dOnafETwNuLR7NT2p71x9VNoPSd4J/Kiq/mIfH/9U4A+q6tRZjvmfdC4288V9q1KanV1D0hBV1RXAxbN9oAy42hBQkzwjkKSW84xAklrOIJCkljMIJKnlDAJJajmDQJJa7v8DhvTo8uHDbewAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.scatter(t_data, h_data)\n", "plt.title(\"Dados\")\n", "plt.xlabel(\"tempo (s)\"); plt.ylabel(\"altura (m)\");" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[1. 0. 0. ]\n", " [1. 0.1 0.01]\n", " [1. 0.2 0.04]\n", " [1. 0.3 0.09]\n", " [1. 0.4 0.16]\n", " [1. 0.5 0.25]\n", " [1. 0.6 0.36]\n", " [1. 0.7 0.49]\n", " [1. 0.8 0.64]\n", " [1. 0.9 0.81]\n", " [1. 1. 1. ]]\n" ] } ], "source": [ "A = np.hstack([np.ones((n,1)), t_data, t_data**2]) # Matriz do sistema sobredeterminado\n", "print(A)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[11. 5.5 3.85 ]\n", " [ 5.5 3.85 3.025 ]\n", " [ 3.85 3.025 2.5333]]\n" ] } ], "source": [ "G = np.dot(A.transpose(), A) # Matriz do sistema normal\n", "print(G)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[19.63208 ]\n", " [ 8.386632]\n", " [ 4.995481]]\n" ] } ], "source": [ "d = np.dot(A.transpose(), h_data) # Lado direito do sistema normal\n", "print(d)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 1.65415252 3.90262443 -5.20208625]\n" ] } ], "source": [ "a = np.linalg.solve(G, d) # Solução do sistema normal\n", "a = np.reshape(a, (3,)) # Formato mais fácil para se trabalhar\n", "print(a)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "10.404172494172505\n" ] } ], "source": [ "g = -2*a[2]; print(g) # Aproximação para a aceleração da gravidade" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "def h(t): return a[0] + a[1]*t + a[2]*t**2 # Altura em função do tempo usando os coeficientes calculados" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAyzElEQVR4nO3dd3gU5drH8e8dCCUkFOnSAopIB+lSVBApShCwRzyIiBQRRbHhqx456EEUy+FYIqAHCUUUBCkiKBqRGqSDCtKVEpESCCWQ+/1jFgwYSEJ2Mpvs/bmuvbI7MzvzG0j23pln5nlEVTHGGBO8QrwOYIwxxltWCIwxJshZITDGmCBnhcAYY4KcFQJjjAlyVgiMMSbIWSEwxpggZ4XABDQRuUdE4kXkiIjsFpE5ItLC61yBQkQ+EpF/eZ3D5GxWCEzAEpFBwJvAy0BpoCLwDtD5EtaV16/hjMlFrBCYgCQiRYCXgP6qOlVVj6pqsqp+oaqDfcuc821YRK4XkV2pXm8TkadEZA1w1Pf80/O285aIvO17fr+IbBSRRBHZIiIPpZOxp2/5AyIyV0QqXWC5SBFREektIr/7jmyeSDU/v4i86Zv3u+95/tT7JCKPi8g+33vv983rDUQDT/qOmL7wTa8uIt+KyEERWS8iUam21VFENvj28bfUOUzwskJgAlUzoAAwLYvruRu4GSgKTAI6ikgEgIjkAe4AJviW3QfcAhQG7gfeEJFr0lqpiHQGngW6AiWB74GJ6WS5AagK3AQ8JSI3+qYPAZoC9YC6QGPguVTvKwMUAcoBDwD/FZFiqhoDxAKvqmq4qnYSkVDgC+AroBQwAIgVkWq+dY0BHlLVCKAW8E06mU0QsEJgAlVx4A9VPZXF9bytqjtV9Ziqbgd+BLr45rUGklR1CYCqzlLVX9XxHc6HacsLrLcP8IqqbvRlfBmod6GjAp9/+o5s1gIf4hQpcL7Vv6Sq+1Q1Afgn0D3V+5J985NVdTZwBKhG2poC4cC/VfWkqn4DzEy1rWSghogUVtUDqvrjRfKaIGGFwASq/UAJP5zb33ne6wn89aF4D38dDSAiHURkiYj8KSIHgY5AiQustxLwlu/0y0HgT0BwvrVnJMt24HLf88t9r9OaB7D/vIKYhPNhn5bLgZ2qmnLe+s7k6oazX9tF5DsRaXaRvCZIWCEwgWoxcAK49SLLHAXCUr0uk8Yy53evOwW4XkTK4xwZTADnPD3wGfAaUFpViwKzcT7c07IT5xRL0VSPgqq66CJ5K6R6XhH43ff8d5zCkta89Jy/f78DFUQk9d92ReA3AFVdrqqdcU4bfQ58ksHtmFzMCoEJSKp6CHge53z4rSISJiKhvm/tr/oWW4Vzzv8yESkDPJqB9SYA3+Kcmtmqqht9s/IB+YEE4JSIdMA5l38h7wHPiEhNcBq3ReT2dDb/f779qInTBjHZN30i8JyIlBSREr79Hp/evvjsBaqker0U54jhSd+/1/VAJ2CSiOQTkWgRKaKqycBhIOX8FZrgY4XABCxVfR0YhNNwmoDzLfxhnG+yAB8Dq4FtOOfzJ/9tJWmbANxIqtNCqpoIPILzDfkAzmmjGRfJNg0YjvMBexhYB3RIZ7vfAZuBr4HXVPUr3/R/AfHAGmAtTjtGRu8NGINzzv+giHyuqidxPvg7AH/gXG57n6r+5Fu+O7DNl7kPTvuECXJiA9MY4y4RiQS2AqF+aPw2xu/siMAYY4KcFQJjjAlydmrIGGOCnB0RGGNMkMtxHXGVKFFCIyMjvY5hjDE5yooVK/5Q1ZJpzctxhSAyMpL4+HivYxhjTI4iItsvNM9ODRljTJCzQmCMMUHOCoExxgS5HNdGYIwxF5KcnMyuXbs4fvy411E8U6BAAcqXL09oaGiG32OFwBiTa+zatYuIiAgiIyMRuVDHsbmXqrJ//3527dpF5cqVM/w+OzVk/C42FiIjISTE+Rkb63UiEyyOHz9O8eLFg7IIAIgIxYsXz/QRkR0RGL+KjYXevSEpyXm9fbvzGiDa+rk02SBYi8AZl7L/dkRg/GrIkL+KwBlJSc50Y0xgskJg/GrHjsxNNyY3e/HFF3nttdcyvHx4+IVGIHWXFQLjVxUrZm66MV6y9iyHFQLjV8OGQVjYudPCwpzpxgSSM+1Z27eD6l/tWVktBsOGDeOqq66iRYsW/PzzzwB88MEHNGrUiLp169KtWzeSfOdPt27dSrNmzahduzbPPffc2XWoKoMHD6ZWrVrUrl2byZOdwfd2795Nq1atqFevHrVq1eL777/PWlgfKwTGr6KjISYGKlUCEednTIw1FJvA40Z71ooVK5g0aRKrVq1i9uzZLF++HICuXbuyfPlyVq9eTfXq1RkzZgwAAwcOpG/fvqxdu5ayZcueXc/UqVNZtWoVq1evZv78+QwePJjdu3czYcIE2rVrd3ZevXr1Lj1sKnbVkPG76Gj74DeBz432rO+//54uXboQ5jssjoqKAmDdunU899xzHDx4kCNHjtCuXTsAfvjhBz777DMAunfvzlNPPQXAwoULufvuu8mTJw+lS5fmuuuuY/ny5TRq1IiePXuSnJzMrbfe6rdCYEcExpiglJ3tWT169GDUqFGsXbuWF1544Zzr/DNzuWerVq2Ii4ujXLly9OjRg3HjxvklnxUCY0xQcqM9q1WrVnz++eccO3aMxMREvvjiCwASExMpW7YsycnJxKZqhGjevDmTJk0COGd6y5YtmTx5MqdPnyYhIYG4uDgaN27M9u3bKV26NA8++CC9evXixx9/vPSwqdipIWNMUDpz+nLIEOd0UMWKThHIymnNa665hjvvvJO6detSqlQpGjVqBMDQoUNp0qQJJUuWpEmTJiQmJgLw1ltvcc899zB8+HA6d+58dj1dunRh8eLF1K1bFxHh1VdfpUyZMvzvf/9jxIgRhIaGEh4e7rcjghw3ZnHDhg3VBqYxxqRl48aNVK9e3esYnkvr30FEVqhqw7SWt1NDxhgT5FwrBCJSQUQWiMgGEVkvIgPTWOZ6ETkkIqt8j+fdymOMMSZtbrYRnAIeV9UfRSQCWCEi81R1w3nLfa+qt7iYwxhjzEW4dkSgqrtV9Uff80RgI1DOre0ZY4y5NNnSRiAikUB9YGkas5uJyGoRmSMiNS/w/t4iEi8i8QkJCW5GNcaYoON6IRCRcOAz4FFVPXze7B+BSqpaF/gP8Hla61DVGFVtqKoNS5Ys6WpeY4wJNq4WAhEJxSkCsao69fz5qnpYVY/4ns8GQkWkhJuZjHGD9WJpcjI3rxoSYAywUVVHXmCZMr7lEJHGvjz73cpkjBvc6sXS5Ezbtm2jVq1aac579NFHiYuLA+DNN9882wspwI033siBAweyJeP53DwiaA50B1qnujy0o4j0EZE+vmVuA9aJyGrgbeAuzWl3uJmgZ6OymYzYv38/S5YsoVWrVsDfC0H37t155513PMnm2uWjqroQuGhvSqo6ChjlVgZjsoONyhaYHv3yUVbtWeXXddYrU48327+Z7nKnT5/mwQcfZNGiRZQrV47p06fz2Wef0b59ewDefvttfv/9d2644QZKlCjBggULiIqKomXLlgzx4BuE3VlsTBbZqGzmfJs2baJ///6sX7+eokWL8tlnn/HDDz/QoEEDAB555BEuv/xyFixYwIIFCwAoVqwYJ06cYP/+7D87bp3O5VKxsf7tTMtc2LBhTptA6tNDNiqb9zLyzd0tlStXPjtWQIMGDdi2bRu7d+8mvaseS5Uqxe+//07x4sWzIeVfrBDkQmcaL898MJ1pvIScWQxOp5wmKTmJY6eOoaqICOI765gnJA+FQguRP29+z/K50Yulydny5//r9zFPnjwcO3aMggULnjMOQVqOHz9OwYIF3Y73N1YIcqGLNV4GyofTgWMH2HpwK1sPbGXbwW3sPrKbhKQE9h3dx76j+/gj6Q8STySSlJzEidMn0l1faEgo4fnCCc8XTpECRShVqBRlwstQulBpyoSXoVxEOaoUq8IVl11BybCSmRoMJCNsVDaTnurVq7N582auv/56ACIiIkhMTKRECeeKeVVlz549REZGZns2KwS5UKA0Xqoqu4/sZu3etazbt451CetYt28dm/Zv4tCJQ+csWyBvAUoXKk3JQiUpE16G2qVqUzh/YcJCw84+CuYtiIhw5sIyRTmdcpqjyUc5cvIIR04eIfFkIgePH2Tvkb0s2bWEPUf2kJR8blUMzxdOlWJVuLrE1dQpVYc6pZ1HxSIV/V4gjDnj5ptv5v3336dXr14A9O7dm/bt259tK1ixYgVNmzYlb97s/1i28QhyochI53TQ+SpVgm3b3Ntu4olE4n+PZ8muJSz5bQlLdi1h39F9Z+eXCS9DrVK1qFa8GpWLVqZyscpnfxbJX8S1D+EjJ4+w89BOthzYwq8HfuXXP3/l1wO/svGPjWw5sOXscoXzF6ZB2QY0K9+MaytcS9PyTSkelr3nak3WBPp4BC1atGDmzJkULVr0b/MGDhxIVFQUbdq0yfJ2MjsegR0R5ELZ1Xh59ORRFu5YyPwt85m/dT5r9q4hRVMAqFa8Gh2u7EDDyxtSu1RtapaqSYkwb24aD88XTvWS1ale8u8fEIknElm3bx1r9q5h9d7VLPttGcN/GM5pPQ04+3Fdpetoe0VbWlduzWUFL8vu+CYXef3119mxY0eahaBWrVp+KQKXwo4Icik3rhpSVdYnrGf6T9OZt2Uei3ct5uTpk+TLk4/mFZrTqlIrmpVvRuNyjSlWsJh/dsQDR08eJf73eBbvWsyinYv4bvt3HD5xGEFoeHlD2lZpS4eqHWhWvhl5QvJ4HdeksnHjRq6++uqgPsWnqvz000+ZOiKwQmAu6nTKaRbtXMTnP33O9J+n8+uBXwGoX6Y+N1a5kRur3EiLii0ICw1LZ00516mUUyz7bRnzfp3HvC3zWLJrCaf1NKUKlaJztc50uboLrSu39vTKJePYunUrERERFC9ePCiLgaqyf/9+EhMTqVy58jnzrBCYTFFVlv62lPFrxvPJ+k9ISEogX558tKnchluvvpVOV3WibERZr2N65tDxQ3y5+Uum/jSV2Ztmc+TkESLyRdCpWifurX0vba9oS94QO+vqheTkZHbt2pXuZZq5WYECBShfvjyhoaHnTLdCYDJk0/5NxK6NZfya8fx64FcK5C1AVLUoulXvRvsr21M4f2GvIwac46eO8/WWr5n20zSmbpzKgeMHKF2oNHfXupvudbtTv0z9oPxmagKPFQJzQSdOnWDqxqm8t+I94rbHIQitK7fm3jr30rV6V/vwz4QTp04we9NsPl7zMTN/mUlySjI1S9akd4Pe3Ff3PooWKOp1RBPErBCYv9n852ZiVsTw4aoP+SPpD64odgUPXvMg99a5l3KFbUTRrPrz2J98sv4Txq4cy/LflxMWGsY9te6hb6O+XFP2Gq/jmSBkhcAAzrn/BdsW8Nqi15izeQ55JA+dr+5MnwZ9aFOlDSFifRC6YcXvK3g3/l0mrJ3AsVPHaFyuMY81fYzbatxmbQkm21ghCHLJp5OZsmEKry16jZV7VlKqUCn6NezHgw0e5PKIy72OFzQOHj/IuNXjGLVsFJv+3ERk0Ugea/oYD9R/gEL5Cnkdz+RyVgiCVFJyEjErYhi5eCQ7D+/k6hJX83izx7m3zr0UyFvA63hBK0VTmPHzDEYsGsGinYsoVqAY/Rr1Y2CTgZQsZGNyG3dYIQgySclJvB//PsN/GM7eo3tpVakVg68dTMeqHe30T4BZtHMRIxaNYPpP0wkLDePhxg/zxLVPeHYXtsm9rBAEifMLQOvKrXnhuhdoVamV19FMOjYmbGRo3FAmrZtEoXyFGNB4AI83e9z6OjJ+Y4UglzuVcoqxK8fywrcvsOfIHtpUbsML171Ay0otvY5mMmlDwgZe+u4lPln/CeH5whnUbBBPXPsE4fnCvY5mcriLFQI7T5CDqSozfp5B7Xdr89DMh7ii2BXE9Yhj/n3zrQjkUDVK1mDSbZNY23ct7a5sxz+/+ydV/1OVmBUxnEo55XU8k0tZIcihlu5aynUfXUfnSZ1RVT6/83O+v/97KwC5RM1SNZly+xSWPLCEKy+7kodmPkSdd+sw85eZ5LSjeBP4rBDkMHuO7OEfn/+DpmOa8sv+X3jv5vdY128dna/ubF0Z5EJNyjchrkcc0+6cxqmUU3Sa2IkbP76R9fvWex3N5CJWCHKI5NPJvLH4Da76z1VMWjeJZ1o8w6YBm3io4UN2U1IuJyLcevWtrO+3nv90+A8rd6+k3vv1eOKrJ0g8keh1PJMLWCHIAb7Z+g313q/HoK8G0aJiC9b1XcfLbV4mIn+E19ECSmysMzpbSIjzMzbW60T+FZonlIcbP8wvA37h/nr3M3LxSKqNqsaEtRPsdJHJEisEASzhaALRU6NpM64Nx5KPMf2u6cy6ZxZVi1f1OlrAiY11RmXbvh1UnZ+9e+e+YgBQIqwEMZ1iWNJrCeUKlyN6ajStx7Xml/2/eB3N5FB2+WgAUlVi18by6JePcvjEYZ5p8QxPt3iagqEFvY4WsLwap9lrp1NOM2blGJ6a/xTHko/x4vUv8nizxwnNE5r+m01QsfsIcpDtB7fTd1Zf5myeQ9PyTRndaTQ1S9X0OlbACwlxjgTOJwIpKdmfJ7vtTtzNgDkD+GzjZ9QvU5/RUaOtl1NzDruPIAdI0RRGLRtFzXdqErc9jrfav8XC+xdaEcigihUzNz23KRtRlk/v+JTP7viM3Ud20/iDxjw9/2mOJR/zOprJAawQBIAdh3bQ9uO2DJgzgBYVW7C+33oeafKIDYyeCcOGQdh5wyaHhTnTg0nX6l3Z0G8DPer1YPgPw2kQ04Afd//odSwT4KwQeEhVGb9mPHXercPSXUv5oNMHzImeQ6WilbyOluNER0NMjNMmIOL8jIlxpgebYgWLMTpqNHPvncuhE4doMroJr3z/CqdTTnsdzQQoayPwyP6k/fSd1ZcpG6ZwbYVrGXfrOK647AqvY5lc5s9jf9J3Vl8+Wf8JzSs0Z1yXcVQpVsXrWMYD1kYQYOb9Oo/a79Zm2k/TeLn1y8T1iLMiYFxxWcHLmNRtErFdY1m3bx1136vL2JVj7b4Dcw7XCoGIVBCRBSKyQUTWi8jANJYREXlbRDaLyBoRydWXOZxKOcWQr4fQbnw7ihYoyrJey3im5TPWFmBcJSLcU/se1vZdS+NyjXlgxgNET422u5LNWW4eEZwCHlfVGkBToL+I1DhvmQ5AVd+jN/Cui3k8tfPQTq7/6HpeXvgyPev3JL53PPXL1vc6lgkiFYpUYF73eQxrPYzJ6yfTIKYBq/as8jqWCQCuFQJV3a2qP/qeJwIbgXLnLdYZGKeOJUBRESnrViavzPplFvXer8fqvasZ32U8o6NGExYalv4bjfGzEAnh2ZbPsuAfCziafJSmo5vyzvJ37FRRkMuWNgIRiQTqA0vPm1UO2Jnq9S7+XiwQkd4iEi8i8QkJCa7l9Lfk08k88dUT3DLxFioUrsCK3iuIrhOEl7GYgNOqUitWPbSK1pVb0392f+749A4OHT/kdSzjEdcLgYiEA58Bj6rq4UtZh6rGqGpDVW1YsmTOGNx739F9tP24La8vfp2+DfuypNcSrip+ldexjDmrZKGSzLxnJsNvHM60jdNoPLoxGxI2eB3LeMDVQiAioThFIFZVp6axyG9AhVSvy/um5WjLf1tOg5gGLP1tKR93+Zh3bn6HAnkLeB3LmL8JkRCebP4k3/zjGw4eP0iT0U2YujGtP1WTm7l51ZAAY4CNqjryAovNAO7zXT3UFDikqrvdypQdxq4cS8sPW5JH8rCo5yLurXOv15GMSVerSq1Y0XsFNUrWoNsn3Xj262ftBrQg4uYRQXOgO9BaRFb5Hh1FpI+I9PEtMxvYAmwGPgD6uZjHVSdPn6TfrH48MOMBWlZqaVcFmRynfOHyxPWIo1f9Xryy8BVunnAzfx770+tYJhvYncV+kHA0gS6Tu/DDzh948tonGdZmmI0aZnK0mBUxPDz7YSoWqcjMe2ZydYmrvY5kssjuLHbRhoQNNBndhBW7VzD5tskMbzvcioDJNm6Nyta7QW++6/EdiScTaTq6KfO3zPfPik1AskKQBV/9+hXNxjTj2KljxPWI446ad3gdyQQRt0dla1ahGUt7LaVCkQq0H9+e9+Pf98+KTcCxQnCJ3l3+Lh1jOxJZNJKlvZbSqFwjryOZIDNkCCQlnTstKcmZ7i+RRSP5oecP3HTFTfSZ1YdBcwdZI3IuZIUgk06nnGbgnIH0m92PDlU7sPD+hVQsEiSjn5iAsmNH5qZfqsL5CzPj7hk80vgR3ljyBrdOvtX6KcplrBBkQlJyEl0md+HtZW/zWNPH+PzOz4nIH+F1LBOksnNUtrwheXmrw1v8t+N/mbNpDq0+asWeI3v8vyHjCSsEGbQ/aT9txrVh5i8z+W/H/zKy3UjrNdR4yotR2fo16scXd3/BL/t/4dox1/LL/l/c25jJNlYIMmD7we00H9uclbtX8ukdn9KvUY693cHkIl6NytahagcW/GMBiScTaT62Oct+W+buBo3rrBCkY83eNTQb04y9R/cyr/s8ulbvmqn3u3V5nzHgfOhv2wYpKc7P7Bqas3G5xizquYiIfBHc8L8bmL1pdvZsGPubcoMVgov4dtu3tPywJSESwvf3f0/LSi0z9X63L+8zxktVi1dl0QOLqFa8GlETo/ho1Ueub9P+ptxhdxZfwLSN07jrs7u4otgVfHnvl5d0ZVBkpPOLer5KlZxvb8bkBodPHKbr5K58vfVrXr3xVQY3H+zatuxv6tLZncWZNH7NeG6fcjvXlL2GhT0v/fLQ7Lq8zxgvFc5fmNnRs7mj5h08Of9Jnl/wvGsD3djflDusL4TzvB//Pn1n9eWGyjcw/a7phOcLv+R1VayY9rcXNy7vM8ZL+fLkY0LXCYSHhjM0bihHTh7h9Ztex+mE2H/sb8oddkSQysjFI+kzqw8dq3Zk5t0zs1QEwJvL+4zxSp6QPHwQ9cHZG88emvmQ3+9Ctr8pd9gRAaCqDI0bygvfvsDtNW5nfNfx5MuTL8vrPXMFx5AhzqFrxYrOL2x2XdlhTHYLkRDebP8mEfkjGPb9MI4mH+Wjzh8RmifUL+u3vyl3BH1jsary1PynGLFoBP+o+w9GR4223kON8YN/L/w3z3z9DJ2rdWbybZPJnze/15GCmjUWX4Cq8tjcxxixaAT9G/VnbOexVgSM8ZOnWzzNqA6jmP7zdLp90o0Tp054HclcQIYLgYgUE5GaIlJFRHJ8AVFVBs0dxFtL32Jgk4H8p8N/CMn5u2VMQOnfuD/v3fweszbNsmIQwC76ySciRUTkWRFZCywB3gc+AbaLyBQRuSE7QmbV+Xcijh+vPPHVE7y59E0GNhnIG+3e8PvVDcYYx0MNHzpbDG6bcpsVgwCU3nmQT4FxQEtVPZh6hog0ALqLSBVVHeNSviw7cyfimX7bt29X7p84mFONRzKg8QArAsZkg4caPoSi9J3Vl9un3M6U26dYm0EAyfWNxefeiajQ9iloPoKIDQ9zaNLbVgSMyUbvLn+XfrP70emqTlYMstnFGosz3DIqInWAyNTvUdWpWU7nsr/uOFS48WloPgKW9SNxztvIZCsCxmSnvo36AtBvdj/u+PQOptw+xS+XapusyVAhEJGxQB1gPZDim6xAwBeCs3ci1vsIWrwKy/vC7FFUqmRFwBgv9G3UF0XpP7s/3ad1Z0LXCTa2h8cyekTQVFVruJrEJcOG+doI1t0FoUkQ35ewMLE7EY3xUL9G/UhKTmLwvMGEh4bzQdQHdtWehzL6L79YRHJkITg7eEe5gkh8fypVDMmWwTuMyc38MSbAE9c+wfOtnmfsqrEMmjvItY7qTPoyekQwDqcY7AFOAAKoqtZxLZkfRUfbB78x/vL3K/Gc15D5v7MXr3+RwycO8+bSNymcvzAv3fCSf8OaDMloIRgDdAfW8lcbgTEmCA0Z8lcROCMpyZme2UIgIoxsN5IjJ48wNG4oEfkiXB3PwKQto4UgQVVnuJrEGJMj+HtMABHhvVve40jyEZ6c/yQR+SPo07DPpQc0mZbRQrBSRCYAX+CcGgJyxuWjxhj/cmNMgDwheRh36ziOnjxKv1n9KBFWgttq3HbpKzSZktHG4oI4BeAmoJPvcYtboYwxgcutMQFC84Qy+bbJXFvhWqKnRvPttm+ztkKTYbn+zmJjjP/Fxro3JsCfx/6k5Yct2XV4F3E94qhbpq5/VhzkLrkbahF5TkQuu8j81iJiRwbGBJnoaGew+JQU56c/r8q7rOBlfBn9JYXzF6Z9bHu2Hdzmv5WbNKV3amgt8IWIfC0iI0TkSRF5XkQ+9vVI2glYmtYbRWSsiOwTkXUXmH+9iBwSkVW+x/NZ2xVjTG5RoUgFvoz+kuOnjtNufDsSjiZ4HSlXu2ghUNXpqtoc6IPTvUQe4DAwHmisqo+p6oX+hz4C2qez/e9VtZ7vYRcQG2POqlmqJjPvnsmOQzu4ZeItHD151OtIuVaGrhpS1U3ApsysWFXjRCTyUkIZYwxA84rNmXzbZLpM7sLtU25nxt0zbBRBF3jduUczEVktInNEpOaFFhKR3iISLyLxCQl2iGhMMImqFsW7N7/LnM1zGDB7gHVF4QIvS+uPQCVVPSIiHYHPgappLaiqMUAMOFcNZVtCY0xA6N2gN1sObGH4D8OpWrwqg5oN8jpSruLZEYGqHlbVI77ns4FQESnhVR5jTGB7uc3L3FbjNp746gmmbZzmdZxcJaPjERQAHgBqAgXOTFfVnpe6YREpA+xVVRWRxjhFaf+lrs8Yk7uFSAjjbh3HzkM7iZ4azXc9vqNRuUZex8oVMnpE8DFQBmgHfAeUBxIv9gYRmQgsBqqJyC4ReUBE+ojImU5EbgPWichq4G3gLrWTf8aYiygYWpAZd8+gdHhpOk3sxPaDafR1YTItQ3cWi8hKVa0vImtUtY6IhOJc+tnU/YjnsjuLjTEbEzbSbEwzyhcuzw89f6BIgSJeRwp4l3xncSrJvp8HRaQWUAQo5Y9wxhiTWdVLVmfqnVP5ef/P3D7ldk6lnPI6Uo6W0UIQIyLFgOeAGcAGYLhrqYwxJh2tK7cm5pYY5m2Zx+NzH/c6To6WbmOxiIQAh1X1ABAHVHE9lTHGZMD99e9n3b51jFwyktqla9Prml5eR8qR0j0iUNUU4MlsyGKMMZk2vO1w2l3Rjn6z+rFwx0Kv4+RIGT01NF9EnhCRCiJy2ZmHq8mMMSYD8obkZdJtk6hcrDJdJ3dlx6FLHCotiGW0ENwJ9Mc5NbTC97BLd4wxAaFogaLMuGsGJ0+fpPOkztZBXSZlqBCoauU0HtZWYIwJGNVKVGNit4ms3rOaHtN7WJ9EmZChQiAi96X1cDucMcZkRoeqHXi17at8uuFT/hX3L6/j5BgZ7XQu9X3cBYA2OJ3GjfN7ImOMyYLHmz3Omr1reP7b56lXph6dqnXyOlLAy+h4BANSvxaRosAkNwIZY0xWiAgxnWLYkLCB7tO6E987nisvu9LrWAHtUnsfPQpU9mcQY4zxlwJ5C/DpHZ+SJyQPXSZ3scbjdGS0jeALEZnhe8wEfsYZP8AYYwJSZNFIJnabyPp96+k9s7c1Hl9ERtsIXkv1/BSwXVV3uZDHGGP85qYrbmLoDUN5bsFzNCnXhEeaPOJ1pICU0VNDHVX1O9/jB1XdJSLW15AxJuA90/IZoqpF8fhXj9udxxeQ0ULQNo1pHfwZxBhj3HBmQJvIopHcPuV2difu9jpSwLloIRCRviKyFmdwmTWpHluBNdkT0RhjsqZIgSJMu3Mah08c5o5P7yD5dHL6bwoi6R0RTAA64XQ93SnVo4Gq3utyNmOM8ZtapWoxJmoMC3cs5Kn5T3kdJ6CkVwjyAIdx+hlKTPXAOp0zxuQ0d9W6iwGNB/DGkjeY/tN0r+MEjPSuGloBnLnmSs6bp9jYBMaYHGZE2xEs2rmIHtN7sLLMSiKLRnodyXMXPSI407mc72Gdzhljcrz8efPzye2fkKIp3PnpnZw8fdLrSJ7L8J3FIlJMRBqLSKszDzeDGWOMW6oUq8LYqLEs+20ZT89/2us4nsvoncW9cMYimAv80/fzRfdiGWOMu7rV6GbtBT4ZPSIYiNMD6XZVvQGoDxx0K5QxxmSHEW1H0KBsA3pM78G2g9u8juOZjBaC46p6HEBE8qvqT0A192IZY4z7rL3AkdFCsMvX9fTnwDwRmQ5sdyuUMcZkF2svyPhQlV1U9aCqvgj8HzAGuNXFXMYYk2261ejGw40e5o0lbzB70+w0l4mNhchICAlxfsbGZmtEV0lO65q1YcOGGh8f73UMY0wuc/zUcRp/0Ji9R/eyps8aSoeXPjsvNhZ694akpL+WDwuDmBiIjvYg7CUQkRWq2jCteZc6MI0xxuQqBfIWYGK3iRw+cZge03uQoiln5w0Zcm4RAOf1kCHZHNIlVgiMMcanZqmajLxpJF9u/pK3l759dvqOHWkvf6HpOY0VAmOMSaVPwz5EVYviqflPsWrPKgAqVkx72QtNz2msEBhjTCoiwpioMRQvWJy7P7ubpOQkhg1z2gRSCwuDYcO8yehvrhUCERkrIvtEZN0F5ouIvC0im31jHFzjVhZjjMmMEmEl+LjLx/z8x88MmjuI6GinYbhSJRBxfuakhuL0uHlE8BHQ/iLzOwBVfY/ewLsuZjHGmExpU6UNg68dzPsr3mfaxmlER8O2bZCS4vzMLUUAXCwEqhoH/HmRRToD49SxBCgqImXdymOMMZk1tPVQGpRtQK8verHr8C6v47jGyzaCcsDOVK93+ab9jYj0FpF4EYlPSEjIlnDGGJMvTz4mdpvI8VPH6Tm9JzntvquMyhGNxaoao6oNVbVhyZIlvY5jjAkiVYtX5fWbXmfelnm8s/wdr+O4wstC8BtQIdXr8r5pxhgTUB5q8BDtr2zP4HmD+fmPn72O43deFoIZwH2+q4eaAodUdbeHeYwxJk1nLiktGFqQ+z6/j1Mpp7yO5FduXj46EVgMVBORXSLygIj0EZE+vkVmA1uAzcAHQD+3shhjTFZdHnE57978Lst+W8Yr37/idRy/Sm/w+kumqnenM1+B/m5t3xhj/O2Omncw/efpvBT3Eh2rdqTB5Q28juQXOaKx2BhjAsWoDqMoXag03ad151jyMa/j+IUVAmOMyYRiBYvxYecP2fjHRp79+lmv4/iFFQJjjMmktle05eFGD/Pm0jf5Zus3XsfJMisExhhzCYa3Hc5Vxa/i/un3k3gi0es4WWKFwBhjLkFYaBgfdv6QnYd28tT8p7yOkyVWCIwx5hJdW+FaHmv6GO/Gv5ujTxFZITDGmCwY2nooV152Jb1m9OLIySNex7kkVgiMMSYLwkLDGBs1lm0Ht/HM/Ge8jnNJrBAYY0wWtazUkkeaPMKo5aP4btt3XsfJNCsExhjjB8NaD6NKsSo8MOMBjp486nWcTLFCYIwxflAoXyHGRo3l1wO/MuSbIV7HyRQrBMYY4yfXRV5H/0b9eXvp2yzcsdDrOBlmhcAYY/zo3zf+m0pFK9Fzek+SkpO8jpMhVgiMMcaPwvOFMyZqDJv+3MSL377odZwMsUJgjDF+1rpya3rV78XIxSNZuXul13HSZYXAGGNc8GrbVykRVoJeX/QK+BHNrBAYY4wLihUsxqiOo/hx94+8ueRNr+NclBUCY4xxSbfq3YiqFsXzC55ny4EtXse5ICsExhjjEhHhvx3/S96QvPSZ2QdnhN7AY4XAGGNcVL5weV5p8wrztsxj/JrxXsdJkxUCY4xxWd9GfWlWvhmPzX2MhKMJXsf5GysExhjjshAJYXTUaA6fOMygrwZ5HedvrBAYY0w2qFGyBs+0eIbxa8Yzd/Ncr+OcwwqBMcZkk2dbPsvVJa6mz6w+AdVDqRUCY4zJJvnz5ifmlhi2HdzGv+L+5XWcs6wQGGNMNmpZqSX317uf1xa/xvp9672OA1ghMMaYbPdq21cpnL8w/Wb3C4h7C6wQGGNMNisRVoLhNw4nbnsc41aP8zqOFQJjjPFCz/o9aVa+GU/Me4I/j/3paRYrBMYY44EQCeG9W97jwLEDPDP/GW+zeLp1Y4wJYnVK12Fgk4HE/BjDkl1LPMvhaiEQkfYi8rOIbBaRp9OY30NEEkRkle/Ry808xhgTaF68/kXKRZSjz8w+no1b4FohEJE8wH+BDkAN4G4RqZHGopNVtZ7vMdqtPMYYE4gi8kfwVvu3WL13NaOWjfIkg5tHBI2Bzaq6RVVPApOAzi5uzxhjcqSu1bvS4coO/N+C/+O3w79l+/bdLATlgJ2pXu/yTTtfNxFZIyKfikgFF/MYY0xAEhFGdRzFqZRTPDb3sWzfvteNxV8AkapaB5gH/C+thUSkt4jEi0h8QkLgdeFqjDFZVaVYFZ5t8SxTNkzh6y1fZ+u23SwEvwGpv+GX9007S1X3q+oJ38vRQIO0VqSqMaraUFUblixZ0pWwxhjjtcHNB1OlWBUe+fIRkk8nZ9t23SwEy4GqIlJZRPIBdwEzUi8gImVTvYwCNrqYxxhjAlqBvAV4o90bbEjYkK0Nx64VAlU9BTwMzMX5gP9EVdeLyEsiEuVb7BERWS8iq4FHgB5u5THGmJyg01Wd6HBlB1787kX2HNmTLduUQOjwKDMaNmyo8fHxXscwxhjX/LL/F2q9U4t7at/DR7d+5Jd1isgKVW2Y1jyvG4uNMcac56riVzGo2SD+t/p/LN652PXtWSEwxpgA9Fyr57g84nIGzBnA6ZTTrm7LCoExxgSg8HzhvNb2NVbsXsGYlWNc3ZYVAmOMCVB31bqLVpVa8ezXz7raVbUVAmOMCVAiwn86/IcDxw/wf9/8n2vbsUJgjDEBrE7pOvRv1J/3VrzHqj2rXNmGFQJjjAlw/7z+n1xW8DL+tyrNXniyLK8razXGGOM3xQoWY1mvZUQWjXRl/VYIjDEmB6hcrLJr67ZTQ8YYE+SsEBhjTJCzQmCMMUHOCoExxgQ5KwTGGBPkrBAYY0yQs0JgjDEBLjYWIiMhJMT5GRvr3/XbfQTGGBPAYmOhd29ISnJeb9/uvAaIjvbPNuyIwBhjAtiQIX8VgTOSkpzp/mKFwBhjAtiOHZmbfimsEBhjTACrWDFz0y+FFQJjjAlgw4ZBWNi508LCnOn+YoXAGGMCWHQ0xMRApUog4vyMifFfQzHYVUPGGBPwoqP9+8F/PjsiMMaYIGeFwBhjgpwVAmOMCXJWCIwxJshZITDGmCAnqup1hkwRkQRg+yW+vQTwhx/j5AS2z8HB9jk4ZGWfK6lqybRm5LhCkBUiEq+qDb3OkZ1sn4OD7XNwcGuf7dSQMcYEOSsExhgT5IKtEMR4HcADts/BwfY5OLiyz0HVRmCMMebvgu2IwBhjzHmsEBhjTJDLlYVARNqLyM8isllEnk5jfn4Rmeybv1REIj2I6VcZ2OdBIrJBRNaIyNciUsmLnP6U3j6nWq6biKiI5PhLDTOyzyJyh+//er2ITMjujP6Wgd/tiiKyQERW+n6/O3qR019EZKyI7BORdReYLyLytu/fY42IXJPljapqrnoAeYBfgSpAPmA1UOO8ZfoB7/me3wVM9jp3NuzzDUCY73nfYNhn33IRQBywBGjode5s+H+uCqwEivlel/I6dzbscwzQ1/e8BrDN69xZ3OdWwDXAugvM7wjMAQRoCizN6jZz4xFBY2Czqm5R1ZPAJKDzect0Bv7ne/4p0EZEJBsz+lu6+6yqC1T1zBDYS4Dy2ZzR3zLy/wwwFBgOHM/OcC7JyD4/CPxXVQ8AqOq+bM7obxnZZwUK+54XAX7Pxnx+p6pxwJ8XWaQzME4dS4CiIlI2K9vMjYWgHLAz1etdvmlpLqOqp4BDQPFsSeeOjOxzag/gfKPIydLdZ98hcwVVnZWdwVyUkf/nq4CrROQHEVkiIu2zLZ07MrLPLwL3isguYDYwIHuieSazf+/pshHKgoyI3As0BK7zOoubRCQEGAn08DhKdsuLc3roepyjvjgRqa2qB70M5bK7gY9U9XURaQZ8LCK1VDXF62A5RW48IvgNqJDqdXnftDSXEZG8OIeT+7MlnTsyss+IyI3AECBKVU9kUza3pLfPEUAt4FsR2YZzLnVGDm8wzsj/8y5ghqomq+pW4BecwpBTZWSfHwA+AVDVxUABnM7ZcqsM/b1nRm4sBMuBqiJSWUTy4TQGzzhvmRnAP3zPbwO+UV8rTA6V7j6LSH3gfZwikNPPG0M6+6yqh1S1hKpGqmokTrtIlKrGexPXLzLyu/05ztEAIlIC51TRlmzM6G8Z2ecdQBsAEamOUwgSsjVl9poB3Oe7eqgpcEhVd2dlhbnu1JCqnhKRh4G5OFccjFXV9SLyEhCvqjOAMTiHj5txGmXu8i5x1mVwn0cA4cAUX7v4DlWN8ix0FmVwn3OVDO7zXOAmEdkAnAYGq2qOPdrN4D4/DnwgIo/hNBz3yMlf7ERkIk4xL+Fr93gBCAVQ1fdw2kE6ApuBJOD+LG8zB/97GWOM8YPceGrIGGNMJlghMMaYIGeFwBhjgpwVAmOMCXJWCIwxJshZITBBRUSKikg/r3OcISIFReQ7EclzkWXmi0ix7MxlgosVAhNsiuL0PhsoegJTVfX0RZb5mMDKbHIZKwQm2PwbuEJEVonICAARGSwiy319u//TNy1SRH4SkY9E5BcRiRWRG32duW0Skca+5V4UkY9FZLFv+oO+6SIiI0RknYisFZE7L5AnGpjue09ZEYnzZVsnIi19y8zA6U/HGFfkujuLjUnH00AtVa0HICI34fTF0xinf/cZItIKp9uCK4Hbcb61LwfuAVoAUcCzwK2+ddbB6cuoELBSRGYBzYB6QF2cfm+Wi0hc6q4AfF0mVFHVbb5J9wBzVXWY71RRGICqHhBnMKXiOfkuYRO4rBCYYHeT77HS9zocpzDsALaq6loAEVkPfK2qKiJrgchU65iuqseAYyKyAKeotAAm+k757BWR74BGnNtPTgngYKrXy4GxIhIKfK6qq1LN2wdcTs7uHNEEKDs1ZIKdAK+oaj3f40pVHeObl7qH1pRUr1M490vU+f20ZLTflmM4HaQ5b3IGJGmF05PkRyJyX6plC/iWN8bvrBCYYJOI00X1GXOBniISDiAi5USkVCbX2VlECohIcZzOwpYD3wN3ikgeESmJ8wG/LPWbfKOI5RGRAr5tVwL2quoHwGic4QrxjZ5XBtiWyVzGZIidGjJBRVX3+xp81wFzVHWwr+vixb5eWY8A9+L03JlRa4AFOKd6hqrq7yIyDaedYDXOEcKTqronjfd+hXMaaT5OERksIsm+HGeOCBoAS3yj6Rnjd9b7qDFZICIvAkdU9bVLfP81wGOq2v0iy7yFM9jM15eW0piLs1NDxnhIVX8EFlzshjJgnRUB4yY7IjDGmCBnRwTGGBPkrBAYY0yQs0JgjDFBzgqBMcYEOSsExhgT5P4fhJA7qFLKE1EAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Gráfico dos dados e da curva obtida\n", "t = np.linspace(0.0, 1.0, 1000)\n", "plt.plot(t_data, h_data, 'bo', label = \"dados\")\n", "plt.plot(t, h(t), 'g-', label = \"h(t)\")\n", "plt.title(\"Curva e pontos\")\n", "plt.xlabel(\"tempo (s)\")\n", "plt.ylabel(\"altura (m)\")\n", "plt.legend() # Show labels" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.9" } }, "nbformat": 4, "nbformat_minor": 2 }