{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np # Paquete para manejar arreglos de datos\n", "import matplotlib.pyplot as plt # Paquete para graficar\n", "from scipy.optimize import curve_fit\t# Se importa metodo de fiteo" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAO1UlEQVR4nO3db2hd933H8c9Xvuq9xZXXsRuVKDe2Css8gqHyuNgww6hCV+w02H5UGrt+1E1PZkjYht2AnxT0ZE9KEOkDi7aModqmkIaNLG5rew5xIJX/xE4Wx3Vlgqt6Dsi3XmTpgW52re8e6EZxXUU6wvfc8/U97xcI3yOdnPv9SdGbw7l/ZO4uAEBcXVkPAABYHqEGgOAINQAER6gBIDhCDQDBFdI4aLlc9v7+/mX3aTQaKhRSufvQWHe+sO58eZB1X7hwoebujyz1tVS+k/39/Tp//vyy+9RqNZXL5TTuPjTWnS+sO18eZN1m9tvP+hqXPgAgOEINAMERagAIjlADQHCEGgCCI9QA0AJTd+b098fe19TMXMuPTagBoAVGTk3o0o0ZjZycaPmx8/eMdABooY2HjqvemF/cHhuf1Nj4pIqFLl0d3tGS++CMGgAewJkDg9o50KdS90JOS91d2jXQpzMHB1t2H4QaAB5A77qSeooF1Rvz+twaU70xr55iQb09pZbdB5c+AOAB1Wbr2rt1g57+ix699psZ3WrxA4qEGgAe0OF9VUkL7/Xx109uaPnxufQBAMElOqM2s+uSZiTdldRw92qaQwEAPrWaSx+D7l5LbRIAwJK49AEAwSU9o3ZJvzQzl3TY3Ufv38HMhiQNSVKlUlGttvzJ9/T09CpH7QysO19Yd76kte6kod7m7jfNrFfSCTP7tbu/ce8OzXiPSlK1WvUkf+Ugj38BQmLdecO68yWNdSe69OHuN5v/Tkl6RdKWlk8CAFjSiqE2s7Vm1vPJbUlfl/Re2oMBABYkufTxJUmvmNkn+x9x95+nOhUAYNGKoXb3DyR9pQ2zAACWwNPzACA4Qg0AwRFqAAiOUANAcIQaAIIj1AAQHKEGgOAINQAER6gBIDhCDQDBEWoACI5QA0BwhBoAgiPUABAcoQaA4Ag1AARHqAEgOEINAMERagAIjlADQHCEGgCCI9QAEByhBoDgCDUABEeoASA4Qg0AwRFqAAiOUANAcIlDbWZrzOyimb2a5kAAgD+0mjPq5yRdSWsQAMDSEoXazCqSviHph+mOAwC4XyHhfi9KOiCp57N2MLMhSUOSVKlUVKvVlj3g9PR0wrvuLKw7X1h3vqS17hVDbWbPSJpy9wtm9tXP2s/dRyWNSlK1WvVyubzinSfZpxOx7nxh3cubujOn/Ucv6qU9m9XbU0p5qvSl8fNOculjm6SdZnZd0jFJT5nZWMsnAZBLI6cmdO76bY2cnMh6lLBWPKN29xckvSBJzTPqf3b3b6c8F4AOt/HQcdUb84vbY+OTGhufVLHQpavDOzKcLB6eRw0gE2cODGrnQJ9K3QsZKnV3addAn84cHMx4sniSPpgoSXL31yW9nsokAHKld11JPcWC6o15FQtdqjfm1VMsdMR16lZbVagBoJVqs3Xt3bpBe7as15Gzk7o1M5f1SCERagCZObyvunh7ePemDCeJjWvUABAcoQaA4Ag1AARHqAEgOEINAMERagAIjlADQHCEGgCCI9QAEByhBoDgCDUABEeoASA4Qg0AwRFqAAiOUANAcIQaAIIj1AAQHKEGgOAINQAER6gBIDhCDQDBEWoACI5QA0BwhBoAgiPUABAcoQaA4FYMtZmVzOysmb1jZpfN7HvtGAwAsKCQYJ+6pKfcfdbMuiW9aWbH3f1XKc8GAFCCULu7S5ptbnY3PzzNoQAAn0pyRi0zWyPpgqQ/l/QDdx9fYp8hSUOSVKlUVKvVlj3m9PT0qoftBKw7X1h3vqS17kShdve7kgbM7IuSXjGzTe7+3n37jEoalaRqterlcnnF4ybZpxOx7nxJuu6pO3Paf/SiXtqzWb09pZSnSh8/79ZZ1bM+3P0jSa9L2t7ySYCcGzk1oXPXb2vk5ETWoyCYFc+ozewRSf/n7h+Z2eclfU3Sv6Q+GZATGw8dV70xv7g9Nj6psfFJFQtdujq8I8PJEEWSM+pHJZ02s3clnZN0wt1fTXcsID/OHBjUzoE+lboXfh1L3V3aNdCnMwcHM54MUSR51se7kja3YRYgl3rXldRTLKjemFex0KV6Y149xUJHXKdGayR6MBFAumqzde3dukF7tqzXkbOTujUzl/VICIRQAwEc3lddvD28e1OGkyAi3usDAIIj1AAQHKEGgOAINQAER6gBIDhCDQDBEWoACI5QA0BwhBoAgiPUABAcoQaA4Ag1AARHqAEgOEINAMERagAIjlADQHCEGgCCI9QAEByhBoDgCDUABEeoASA4Qg0AwRFqAAiOUGNVpu7M6ZuH39LUzFzWowC5QaixKiOnJnTu+m2NnJzIehQgNwpZD4CHw8ZDx1VvzC9uj41Pamx8UsVCl64O78hwMqDzrXhGbWaPm9lpM7tiZpfN7Ll2DIZYzhwY1M6BPpW6F/6XKXV3addAn84cHMx4MqDzJTmjbkj6J3d/28x6JF0wsxPu/n7KsyGQ3nUl9RQLqjfmVSx0qd6YV0+xoN6eUtajAR1vxVC7+4eSPmzenjGzK5Iek0Soc6Y2W9ferRu0Z8t6HTk7qVs8oAi0xaquUZtZv6TNksbTGAaxHd5XXbw9vHtThpMA+ZI41Gb2BUkvS3re3e8s8fUhSUOSVKlUVKvVlj3e9PT06ibtEKw7X1h3vqS17kShNrNuLUT6J+7+s6X2cfdRSaOSVK1WvVwur3jcJPt0ItadL6w7X9JYd5JnfZikH0m64u7fb/kEAIBlJXnByzZJ+yQ9ZWaXmh9PpzwXAKApybM+3pRkbZgFALAEXkIOAMERagAIjlADQHCEGgCCI9QAEByhBoDgCDUABEeoASA4Qg0AwRFqAAiOUANAcIQaAIIj1AAQHKEGgOAINQAER6gBIDhCDQDBEWoACI5QA0BwhBoAgiPUABAcoQaA4Ag1AARHqAEgOEINAMERagAIjlADQHCEGgCCI9QAENyKoTazH5vZlJm9146BHnZTd+b0zcNvaWpmLutRAHSIJGfU/yppe8pzdIyRUxM6d/22Rk5OZD0KgA5RWGkHd3/DzPrTH+XhtvHQcdUb84vbY+OTGhufVLHQpavDOzKcDMDDbsVQJ2VmQ5KGJKlSqahWqy27//T0dKvuOoR//7uv6MXXJ3X62v+q3phXsdClwSf+VM9/df0ffC86bd1Jse58Yd2t1bJQu/uopFFJqlarXi6XV/xvkuzzsCiXpT/7k9/r47u/V7HQpY/vzqu8bq3+ckPfEvt2zrpXg3XnC+tunZaFGlJttq69Wzdoz5b1OnJ2Urd4QBFACxDqFjq8r7p4e3j3pgwnAdBJkjw976iktyRtNLMbZvad9McCAHwiybM+nm3HIACApfHKRAAIjlADQHCEGgCCI9QAEByhBoDgCDUABEeoASA4Qg0AwRFqAAiOUANAcIQaAIIj1AAQHKEGgOAINQAER6gBIDhCDQDBEWoACI5QA0BwhBoAgiPUABAcoQaA4Ag1AARHqAEgOEINAMERagAIjlADQHCEGgCCI9QAEFyiUJvZdjO7ambXzOy7aQ8FAPjUiqE2szWSfiBph6QnJT1rZk+mPRgAYEGSM+otkq65+wfu/rGkY5J2pTsWAOAThQT7PCbpd/ds35C09f6dzGxI0pAkVSoV1Wq1ZQ86PT2dfMoOwrrzhXXnS1rrThJqW+Jz/kefcB+VNCpJ1WrVy+Xyige+f5+pO3Paf/SiXtqzWb09pQSjPZySfG86EevOF9bdOkkufdyQ9Pg92xVJN1s+iaSRUxM6d/22Rk5OpHF4AHgoJTmjPifpCTP7sqT/kfQtSXtaOcTGQ8dVb8wvbo+NT2psfFLFQpeuDu9o5V0BwENnxTNqd29I2i/pF5KuSPqpu19u5RBnDgxq50CfSt0L45S6u7RroE9nDg628m4A4KGU5Ixa7v6apNfSGqJ3XUk9xYLqjXkVC12qN+bVUyx09HVqAEgqUajboTZb196tG7Rny3odOTupWzNzWY8EACGECfXhfdXF28O7N2U4CQDEwnt9AEBwhBoAgiPUABAcoQaA4Ag1AARHqAEgOHP/o/dXevCDmt2S9NsVditLWv4t9joT684X1p0vD7LuDe7+yFJfSCXUSZjZeXevrrxnZ2Hd+cK68yWtdXPpAwCCI9QAEFyWoR7N8L6zxLrzhXXnSyrrzuwaNQAgGS59AEBwhBoAgmt7qM1su5ldNbNrZvbddt9/Vszsx2Y2ZWbvZT1Lu5jZ42Z22syumNllM3su65naxcxKZnbWzN5prv17Wc/ULma2xswumtmrWc/STmZ23cz+28wumdn5lh67ndeozWyNpN9I+lst/NHcc5Kedff32zZERszsbyTNSvo3d8/FG26b2aOSHnX3t82sR9IFSbtz8vM2SWvdfdbMuiW9Kek5d/9VxqOlzsz+UVJV0jp3fybredrFzK5Lqrp7y1/o0+4z6i2Srrn7B+7+saRjkna1eYZMuPsbkm5nPUc7ufuH7v528/aMFv7m5mPZTtUevmC2udnd/Oj4R+7NrCLpG5J+mPUsnaTdoX5M0u/u2b6hnPzi5p2Z9UvaLGk820nap3kJ4JKkKUkn3D0Pa39R0gFJ81kPkgGX9Eszu2BmQ608cLtDbUt8ruPPMvLOzL4g6WVJz7v7naznaRd3v+vuA5IqkraYWUdf8jKzZyRNufuFrGfJyDZ3/ytJOyT9Q/NyZ0u0O9Q3JD1+z3ZF0s02z4A2al6ffVnST9z9Z1nPkwV3/0jS65K2ZzxK2rZJ2tm8VntM0lNmNpbtSO3j7jeb/05JekULl3pbot2hPifpCTP7spl9TtK3JP1Hm2dAmzQfUPuRpCvu/v2s52knM3vEzL7YvP15SV+T9Otsp0qXu7/g7hV379fC7/Z/ufu3Mx6rLcxsbfMBc5nZWklfl9SyZ3i1NdTu3pC0X9IvtPDA0k/d/XI7Z8iKmR2V9JakjWZ2w8y+k/VMbbBN0j4tnFldan48nfVQbfKopNNm9q4WTlBOuHuunq6WM1+S9KaZvSPprKT/dPeft+rgvIQcAILjlYkAEByhBoDgCDUABEeoASA4Qg0AwRFqAAiOUANAcP8PV3Dnjt5pbdoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "x = [0.1, 0.9, 1.8, 3.2, 4.2, 5]\n", "y = [-0.1, 1, 2.1, 2.9, 3.9, 5.2]\n", "plt.plot(x, y, '*')\n", "plt.grid(alpha=0.3)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([ 0.99439421, -0.01913201])" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "z = np.polyfit(x, y, 1)\n", "z" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deVjVZR738fdtGDY+OmZqG5VN+YBboqJpmWWWS6ZpuaTWTOVoc8Als3EZbccoxy0iTNN03C8rl5Q0wiUU11BREdc01CxEc8tAgfv5Q8enKUuEc/id5fO6Li/PkXN+53N79HN9r/v8+GGstYiIiPcq5XQAERH5YypqEREvp6IWEfFyKmoRES+nohYR8XJBnjhopUqVbNWqVYv03Ly8PIKCPBLLa2nN/i/Q1gta85VKTU3NttZWvtTXPPK3WLVqVb7++usiPTc7O5tKlSq5OZF305r9X6CtF7TmK2WM+fb3vqatDxERL6eiFhHxcipqEREvp6IWEfFyKmoRES+nohYRcYOskzn0nL2drFM5bj+2ilpExA1il+5m88FTxCbtdvuxA+tsdBERNwsdtpjcvIKL96evy2T6ukyCg0qxM7q1W15DE7WISDGsHNiMduE3Uab0+TotU7oUj4XfxMpBzdz2GipqEZFiqFK+DOWCg8jNK+Dqqwy5eQWUCw6iSrkybnsNbX2IiBRT9ulcut99G4/833J8vusUR9z8gaKKWkSkmMY/HQGcv9bHPTVuc/vxtfUhIuLlCjVRG2P2A6eAfCDPWhvhyVAiIvL/XcnWRzNrbbbHkoiIyCVp60NExMsVdqK2QKIxxgLjrbUTfv0AY0wvoBdASEgI2dlFG75PnDhRpOf5Mq3Z/wXaeiHw1pyVlcXatWtp166d249d2KK+11r7nTGmCvClMWaHtTb5lw+4UN4TACIiImxxfrJDoP1UCNCaA0GgrRcCY835+flMmDCBf/3rXxhj6NKlC2XLlnXraxRq68Na+92F37OAeUBDt6YQEfFBX3/9NY0aNSIyMpJ69eqRkJDg9pKGQhS1MaasMabcf28DLYBtbk8iIuIjfvzxRyIjI2nYsCEHDx5k5syZJCUlUa1aNY+8XmG2Pq4H5hlj/vv4mdbaJR5JIyLixay1TJs2jZdeeomjR4/Sp08f3njjDf785z979HUvW9TW2m+AOh5NISLi5dLT04mMjCQ5OZlGjRrxxRdfULdu3RJ5bZ2eJyLyB06fPs3AgQMJDw9n27ZtfPjhh6SkpJRYSYOu9SEicknWWubNm0e/fv04ePAgPXr04O2333bkTBZN1CIiv7J3717atGnDE088QcWKFUlJSWHixImOnW6oohYRuSAnJ4c33niDmjVrsnLlSsaMGUNqair33HOPo7m09SEiAiQmJhIVFcWePXvo0qULo0aN4uabb3Y6FqCJWkQC3KFDh+jcuTMtW7bEGENiYiKzZ8/2mpIGFbWIBKhz584xevRowsLCWLhwIW+++SZbt27l4Ycfdjrab2jrQ0QCTkpKCi6Xi61bt/LII4/w3nvv8Ze//MXpWL9LE7WIBIwjR47w3HPP0aRJE44fP868efNYtGiRV5c0qKhFJAAUFBQwYcIEQkNDmTZtGoMGDSIjI4P27dtz4fIYXk1bHyLi1zZu3IjL5WL9+vXcf//9xMfHU6NGDadjXRFN1CLil06cOEHfvn1p0KAB+/fvZ9q0aSxfvtznSho0UYuIn7HWMmvWLF588UWysrKIjIwkOjqaChUqOB2tyFTUIuI3MjIyiIqKYvny5TRo0ICEhATq16/vdKxi09aHiPi8n376iSFDhlCnTh02bdrEuHHjWLNmjV+UNGiiFhEft2DBAvr27UtmZiZ/+9vfGDFiBFWqVHE6lltpohYRn7Rv3z7atWtH+/btKV++PMnJyUyZMsXvShpU1CLiY3Jzc3nrrbeoWbMmy5YtY+TIkWzcuJH77rvP6Wgeo60PEfEZS5cuJSoqip07d9KxY0fGjBlDSEiI07E8ThO1iHi97777jq5du/LQQw+Rl5fH4sWL+fjjjwOipEFFLSJeLC8vj3fffZewsDDmzZvHq6++yrZt22jVqpXT0UqUtj5ExCutWbMGl8tFWloaLVu2JC4ujjvvvNPpWI7QRC0iXuXo0aP07NmTe+65h+zsbD755BMWL14csCUNKmoR8RIFBQVMmjSJ0NBQJk+ezEsvvURGRgZPPPGET1zhzpO09SEijktLS8PlcrFmzRqaNGlCfHw8tWvXdjqW19BELSKOOXnyJP3796d+/frs3r2bKVOmkJycrJL+FU3UIlLirLXMmTOH/v378/333/P8888zfPhwKlas6HQ0r6SiFpEStWvXLqKiokhKSqJevXrMnz+fhg0bOh3LqxV668MYc5UxZpMxZpEnA4mIf/r55595+eWXqV27Nhs2bCAuLo7169erpAvhSibqfkAGUN5DWUTETyUkJNCnTx/27dvHU089xb///W9uuOEGp2P5jEJN1MaYEKANMNGzcUTEn3z77bd06NCBRx99lDJlyrB8+XKmTZumkr5ChZ2oxwIDgXK/9wBjTC+gF0BISAjZ2dlFCnTixIkiPc+Xac3+L9DWe/bsWcaOHcv7778PwMsvv8w//vEPrr766iJ3gy/w1Pt82aI2xjwKZFlrU40xD/ze46y1E4AJABEREbZSpUpFDlWc5/oqrdn/Bcp6V6xYQWRkJBkZGbRv356xY8dy2223/eFzsk7m0HvWJuK61aVKuTIllNQzPPE+F2br416gnTFmPzAbeNAYM93tSUTEp33//fc8/fTTNGvWjJycHGbMmMG8efMuW9IAsUt3s2H/MWKTdpdAUt9z2YnaWjsEGAJwYaJ+yVr7lIdziYiPyM/P54MPPmDo0KH8/PPPDBs2jCFDhnDmzJnLPjd02GJy8wou3p++LpPp6zIJDirFzujWnoztU/SdiSJSZP89va537940aNCALVu28Oabb/KnP/2pUM9fObAZ7cJvokzp81VUpnQpHgu/iZWDmnkyts+5oqK21q6w1j7qqTAi4huOHTuGy+WiUaNGHD58mNmzZ5OYmEhoaOgVHadK+TKUCw4iN6+A4KBS5OYVUC44yOf3qd1N35koIoVmreU///kP//znPzl27Bj9+vXj9ddfp3z5on97RfbpXLrffRvdGt7KzPWZHDmV48bE/kFFLSKFsnXrViIjI1m1ahWNGzcmPj6e8PDwYh93/NMRF29Ht69V7OP5I+1Ri8gfOn36NP/85z+pW7cuGRkZTJo0iVWrVrmlpKVwNFGLyCVZa5k7dy79+vXj0KFD9OzZk5iYGK677jqnowUcTdQi8ht79uzhkUceoWPHjlSqVInVq1czYcIElbRDVNQiclFOTg6vvfYatWrVIiUlhXfffZevv/6axo0bOx0toGnrQ0QAWLJkCb1792bv3r08+eSTjBo1iptuusnpWIImapGAd+DAATp27Ejr1q0JCgoiKSmJWbNmqaS9iIpaJECdO3eOkSNHUr16dRISEhg+fDhpaWk0b97c6WjyK9r6EAlAK1euxOVykZ6ezqOPPkpsbCy3336707Hkd2iiFgkgWVlZPPPMMzRt2pRTp06xYMECFi5cqJL2cipqkQDw3yvchYaGMnPmTIYMGcL27dtp166d09GkELT1IeLnUlNTcblcbNiwgWbNmvH+++9TvXp1p2PJFdBELeKnjh8/fvHyo5mZmcyYMYOlS5eqpH2QilrEz1hrmT59OqGhoYwbN47evXuzc+dOunXrhjHG6XhSBNr6EPEj27dvJzIykq+++oqGDRuyePFi6tWr53QsKSZN1CJ+4KeffmLw4MHUqVOHLVu2MH78eNasWaOS9hOaqEV8mLWWBQsW0LdvXw4cOMCzzz7LO++8Q+XKlZ2OJm6kiVrER33zzTe0bduWDh06UKFCBVauXMlHH32kkvZDKmoRH5Obm0t0dDQ1a9bkq6++YtSoUaSmptKkSROno4mHaOtDxId8+eWXREVFsXv3bjp16sTo0aMJCQlxOpZ4mCZqER9w6NAhnnzySVq0aIG1liVLljBnzhyVdIBQUYt4sby8PMaMGUNYWBjz58/n9ddfZ+vWrbRs2dLpaFKCtPUh4qVSUlKIjIxky5YttG7dmvfee4877rjD6VjiAE3UIl4mOzubHj160KRJE44dO8bcuXNJSEhQSQcwFbWIlygoKGDixImEhoYydepUBg4cSEZGBh06dNC3fgc4bX2IeIHNmzfjcrlYu3YtTZs2JT4+npo1azodS7yEJmoRB504cYJ+/fpRv3599u7dy9SpU1mxYoVKWv7HZSdqY0wZIBkIvvD4T6y1r3o6mIg/s9Yye/ZsXnzxRX744QdcLhfR0dFce+21TkcTL1SYrY9c4EFr7WljTGlglTFmsbV2rYezifilHTt2EBUVxbJly4iIiOCzzz6jQYMGTscSL3bZrQ973ukLd0tf+GU9mkrED505c4ahQ4dy1113kZqaSnx8PGvXrlVJy2UV6sNEY8xVQCpwJ/C+tXbdJR7TC+gFEBISQnZ2dpECnThxokjP82Vas/+bO3cu0dHRHDhwgM6dO/Pqq69SpUoVfvzxR6ejeUygvcfguTUXqqittflAuDGmAjDPGFPLWrvtV4+ZAEwAiIiIsJUqVSpyqOI811dpzf5p//799OvXj88++4waNWqwYsUK7r///ss+L+tkDr1nbSKuW12qlCtTAkk9IxDe41/zxJqv6KwPa+1xYAXQyu1JRPzI2bNniYmJoUaNGiQlJfHqq6+yefPmQpU0QOzS3WzYf4zYpN0eTiq+oDBnfVQGzllrjxtjrgEeAt7xeDIRH7Vs2TKioqLYsWMHjz/+OGPGjOFPf/oTpUuXvuxzQ4ctJjev4OL96esymb4uk+CgUuyMbu3J2OLFCjNR3wgsN8ZsATYAX1prF3k2lojvOXz4MN27d6d58+acPXuWhIQEPv30U2699dZCH2PlwGa0C7+JMqXP/9csU7oUj4XfxMpBzTwVW3zAZSdqa+0WoG4JZBHxSXl5eYwbN45hw4aRk5PDK6+8wuDBg7nmmmuu+FhVypehXHAQuXkFBAeVIjevgHLBQT69Ty3Fp28hFymGdevW4XK52LRpEy1atCAuLo5q1aoV65jZp3PpfvdtdGt4KzPXZ3LkVI6b0oqvUlGLFMHRo0f517/+xYcffsiNN97InDlz6Nixo1sunjT+6YiLt6Pb1yr28cT36VofIlegoKCAyZMnExYWxqRJk+jfvz87duygU6dOusKdeIwmapFC2rJlC5GRkaSkpHDvvfcSHx/PXXfd5XQsCQCaqEUu49SpUwwYMIB69eqxc+dOJk+eTHJyskpaSowmapHfYa3lk08+4YUXXuDw4cP07NmTmJgYKlas6HQ0CTCaqEUuYffu3bRq1YrOnTtz/fXXs3r1asaPH6+SFkeoqEV+4eeff+aVV16hVq1arF27ltjYWNavX0+jRo2cjiYBTFsfIhd8/vnn9O7dm3379tGtWzdGjhzJjTfe6HQsEU3UIpmZmTz++OO0adOG4OBgli5dyowZM1TS4jVU1BKwzp07x4gRI6hevTpLliwhJiaGtLQ0HnzwQaejifwPbX1IQEpOTsblcrF9+3Yee+wxxo4dS9WqVZ2OJXJJmqgloPzwww/87W9/4/777+fMmTN89tlnzJ8/XyUtXk1FLQEhPz+f+Ph4QkNDmTVrFkOHDiU9PZ22bds6HU3ksrT1IX5vw4YNuFwuUlNTad68OXFxcYSFhTkdS6TQNFGL3/rxxx+JjIzk7rvv5tChQ8yaNYsvv/xSJS0+R0Utfsday9SpUwkNDWX8+PH07duXHTt28OSTT+oKd+KTtPUhfiU9PZ3IyEiSk5Np1KgRiYmJhIeHOx1LpFg0UYtfOH36NAMHDiQ8PJxt27bx4YcfkpKSopIWv6CJWnyatZZ58+bRr18/Dh48SI8ePXj77bepVKmS09FE3EYTtfisvXv30qZNG5544gkqVqxISkoKEydOVEmL31FRi8/JycnhjTfeoGbNmqxcuZIxY8aQmprKPffc43Q0EY/Q1of4lMTERKKiotizZw9dunRh1KhR3HzzzU7HEvEoTdTiEw4dOkTnzp1p2bIlxhgSExOZPXu2SloCgopavNq5c+cYPXo0YWFhLFy4kDfffJOtW7fy8MMPOx1NpMSoqKVYsk7m0Hn8GrJO5bj92KtWraJ+/foMGDCApk2bkp6ezrBhwwgODnb7a4l4MxW1FEvs0t1s2H+M2KTdbjvmkSNHeO6557jvvvs4fvw48+bNY9GiRfzlL39x22uI+BJ9mChFEjpsMbl5BRfvT1+XyfR1mQQHlWJndOsiHbOgoICJEycyePBgTp06xeDBgxk2bBhly5Z1V2wRn3TZidoYc4sxZrkxJsMYk26M6VcSwcS7rRzYjHbhN1Gm9Pl/QmVKl+Kx8JtYOahZkY63ceNGGjduzPPPP0+dOnVIS0sjJiZGJS1C4bY+8oAB1trqQCMgyhhTw7OxxNtVKV+GcsFB5OYVEBxUity8AsoFB1GlXJkrOs7x48fp06cPDRo0YP/+/UybNo1ly5ZRo4b+iYn812W3Pqy1h4HDF26fMsZkADcD2z2cTbxc9ulcut99G90a3srM9ZkcuYIPFK21zJgxgwEDBpCVlUVkZCTR0dFUqFDBg4lFfNMV7VEbY6oCdYF1nggjvmX80xEXb0e3r1Xo52VkZNCzZ09SUlJo0KABCQkJ1K9f3xMRRfxCoYvaGPN/gE+BF6y1Jy/x9V5AL4CQkBCys7OLFOjEiRNFep4vC5Q1//TTT4wePZpx48ZxzTXXMHLkSJ566imuuuqqIv978RWB8h7/ktbsPoUqamNMac6X9Axr7dxLPcZaOwGYABAREWGLc2GcQLyojr+vecGCBfTt25fMzEyeeeYZBg4cSPXq1Z2OVaL8/T2+FK3ZPQpz1ocBJgEZ1trRbk8gfm3fvn20a9eO9u3bU758eZKTk5k8eTKVK1d2OpqIzyjMWR/3Ak8DDxpjNl/49YiHc4mPy83NZfjw4dSoUYNly5YxcuRINm7cyH333ed0NBGfU5izPlYB+kFzUmhJSUlERUWxa9cuOnbsyJgxYwgJCXE6lojP0reQi9t89913dO3alYcffpj8/HwWL17Mxx9/rJIWKSYVtRRbXl4e7777LmFhYcybN4/XXnuNbdu20apVK6ejifgFXetDimXNmjW4XC7S0tJo2bIlcXFx3HnnnU7HEvErmqilSI4ePUrPnj255557yM7O5pNPPmHx4sUqaREPUFHLFSkoKGDSpEmEhoYyefJkXnrpJTIyMnjiiSc4fyaniLibtj6k0NLS0nC5XKxZs4YmTZoQHx9P7dq1nY4l4vc0UctlnTx5kv79+1OvXj327NnDlClTSE5OVkmLlBBN1PK7rLXMmTOH/v378/333/P8888zfPhwKlas6HQ0kYCiopZL2rVrF1FRUSQlJVGvXj3mz59Pw4YNnY4lEpC09SH/48yZM7z88svUrl2bDRs2EBcXx/r161XSIg7SRC0XLVq0iD59+rB//36eeuop/v3vf3PDDTc4HUsk4GmiFr799ls6dOhA27Ztueaaa1i+fDnTpk1TSYt4CRV1ADt79izvvPMONWrUIDExkXfeeYfNmzfzwAMPOB1NRH5BWx8BasWKFURGRpKRkUGHDh0YO3Yst956q9OxROQSNFEHmO+//56nn36aZs2akZOTw6JFi5g7d65KWsSLqagDRH5+PnFxcYSFhTFnzhyGDRvGtm3baNOmjdPRROQytPURANavX4/L5WLjxo089NBDxMXFERoa6nQsESkkTdR+7NixY/zjH/+gUaNGHD58mNmzZ5OYmKiSFvExKmo/ZK1lypQphIaGMnHiRF544QV27NhBly5ddIU7ER+kovYzW7dupWnTpjz77LNUq1aN1NRURo8eTfny5Z2OJiJFpKL2E6dOneKll16ibt26ZGRkMGnSJFatWkWdOnWcjiYixaQPE32ctZZPP/2UF154gUOHDtGzZ09iYmK47rrrnI4mIm6iidqH7dmzh9atW9OpUycqV67M6tWrmTBhgkpaxM+oqH1QTk4Or732GrVq1WL16tW8++67bNiwgcaNGzsdTUQ8QFsfPmbJkiX07t2bvXv30rVrV0aNGsWNN97odCwR8SBN1D7iwIEDdOzYkdatWxMUFERSUhIzZ85USYsEABW1lzt37hwjR46kevXqJCQkMHz4cNLS0mjevLnT0USkhGjrw4utXLkSl8tFeno6jz76KLGxsdx+++1OxxKREqaJ2gtlZWXxzDPP0LRpU06fPs2CBQtYuHChSlokQF22qI0xHxljsowx20oikL/KOplD5/FryDqV87uPyc/P54MPPiA0NJSZM2cyZMgQ0tPTadeuXQkmFRFvU5iJegrQysM5/F7s0t1s2H+M2KTdl/x6amoqjRs3xuVyUbduXdLS0njrrbcoW7ZsCScVEW9z2T1qa22yMaaq56P4p9Bhi8nNK7h4f/q6TKavyyQ4qBQ7o1tz/PhxBg0axOTJk6lSpQozZsyga9euuniSiFzktg8TjTG9gF4AISEhZGdnF+k4J06ccFckr7Dg73UYuyKT5Xt+JDevgOCgUjSrdi397r+FcePG8eqrr3L06FF69OjBkCFDKF++PEePHnU6tsf52/t8OYG2XtCa3cltRW2tnQBMAIiIiLCVKlUq8rGK81xvU6kSXPfno5zNP0pwUCnO5hdQKv8srme689VXX9GwYUNmzZoVkKfb+dP7XBiBtl7Qmt1FZ32UgOzTuXS/+zZmPluPO+x3fJqQyJYtWxg/fjxr1qzRFe5E5A/pPOoS8MFT9VmwYAEdmj3CgQMHeO6553h7504qV67sdDQR8QGFOT1vFrAGCDXGHDTG9PB8LP/xzTff0LZtWzp06ECFChVYtWoVkyZNUkmLSKEV5qyPriURxN/k5uYyYsQI3nrrLYKCghg1ahR9+vShdOnSTkcTER+jrQ8P+PLLL4mKimL37t106tSJ0aNHExIS4nQsEfFR+jDRjQ4dOkSXLl1o0aIF1lq++OIL5syZo5IWkWJRUbtBXl4eY8aMISwsjAULFvDGG2+wdetWWrRo4XQ0EfED2vooppSUFCIjI9myZQutW7fmvffe44477nA6loj4EU3URZSdnU2PHj1o0qQJx44dY+7cuSQkJKikRcTtVNRXqKCggA8//JDQ0FCmTp3KwIEDycjIoEOHDro+h4h4hLY+rsCmTZtwuVysW7eOpk2bEh8fT82aNZ2OJSJ+ThN1IZw4cYJ+/foRERHBvn37mDp1KitWrFBJi0iJ0ET9B6y1zJ49mxdffJEffvgBl8tFdHQ01157rdPRRCSAqKh/x44dO4iKimLZsmVERESwcOFCIiIinI4lIgFIWx+/cubMGYYOHcpdd91Famoq8fHxrF27ViUtIo7RRP0LCxcupE+fPnz77bf89a9/ZcSIEVx//fVOxxKRAKeJGti/fz+PPfYY7dq1o2zZsqxYsYL//Oc/KmkR8QoBXdRnz54lJiaGGjVqkJSUxIgRI9i8eTP333+/09FERC4K2K2PZcuWERUVxY4dO3j88ccZO3Yst9xyi9OxRER+I+Am6sOHD9O9e3eaN2/O2bNnSUhI4NNPP1VJi4jXCpiizsvLIzY2lrCwMD755BNeeeUVtm3bxiOPPOJ0NBGRPxQQWx9r167F5XKxefNmWrRoQVxcHNWqVXM6lohIofj1RH306FF69epF48aNOXLkCHPmzGHJkiUqaRHxKX5Z1AUFBUyePJmwsDA++ugjBgwYQEZGBp06ddIV7kTE5/jd1seWLVuIjIwkJSWFe++9l3HjxlG7dm2nY4mIFJnfTNSnTp1iwIAB1KtXj507dzJ58mSSk5NV0iLi83x+orbW8vHHH9O/f38OHz5Mz549iYmJoWLFik5HExFxC5+eqHft2kXLli3p0qUL119/PatXr2b8+PEqaRHxKz5Z1D///DOvvPIKtWvXZt26dcTGxrJ+/XoaNWrkdDQREbfzua2Pzz//nN69e7Nv3z66d+/OyJEjueGGG5yOJSLiMT4zUWdmZvL444/Tpk0bgoODWbZsGdOnT1dJi4jf8/qiPnv2LCNGjKB69eosWbKEmJgY0tLSaNasmdPRRERKRKGK2hjTyhiz0xizxxgz2NOh/uurr76ibt26DBo0iIcffpjt27czePBgrr766pKKICLiuMsWtTHmKuB9oDVQA+hqjKnhyVA//PADf/3rX3nggQc4c+YMn332GfPnz6dq1aqefFkREa9UmIm6IbDHWvuNtfYsMBt4zBNh8vPz+eijjwgNDWX27NkMHTqU9PR02rZt64mXExHxCYU56+Nm4MAv7h8E7v71g4wxvYBeACEhIWRnZ19RkOPHj9OxY0fS0tJo2rQp77zzDnfeeSdnzpzhzJkzV3QsX3PixAmnI5S4QFtzoK0XtGZ3KkxRX+oqRvY3f2DtBGACQEREhK1UqdIVBbnuuusIDw8nKiqKv//97/9z8aSskzn0nrWJuG51qVKuzBUd11dc6d+XPwi0NQfaekFrdpfCbH0cBH75409CgO/cHcQYw5QpU+jQocNvrnAXu3Q3G/YfIzZpt7tfVkTE6xVmot4AVDPG3A4cAp4Eunk01QWhwxaTm1dw8f70dZlMX5dJcFApdka3LokIIiKOu+xEba3NA3oDXwAZwBxrbbqngwGsHNiMduE3Uab0+ZhlSpfisfCbWDlI51CLSOAo1LeQW2s/Bz73cJbfqFK+DOWCg8jNKyA4qBS5eQWUCw7y231qEZFL8fprfWSfzqX73bfRreGtzFyfyZFTOU5HEhEpUV5f1OOfjrh4O7p9LQeTiIg4w+uv9SEiEuhU1CIiXk5FLSLi5VTUIiJeTkUtIuLlVNQiIl7OWPub6ysV/6DGHAG+LeLTKwFXduk936c1+79AWy9ozVfqNmtt5Ut9wSNFXRzGmK+ttRGXf6T/0Jr9X6CtF7Rmd9LWh4iIl1NRi4h4OW8s6glOB3CA1uz/Am29oDW7jdftUYuIyP/yxolaRER+QfR8wiYAAAK2SURBVEUtIuLlvKaojTGtjDE7jTF7jDGDnc5TEowxHxljsowx25zOUhKMMbcYY5YbYzKMMenGmH5OZ/I0Y0wZY8x6Y0zahTW/7nSmkmKMucoYs8kYs8jpLCXBGLPfGLPVGLPZGPO1W4/tDXvUxpirgF3Aw5z/YbobgK7W2u2OBvMwY0xT4DQw1Vrr9xfbNsbcCNxord1ojCkHpALt/fl9Nud/UnNZa+1pY0xpYBXQz1q71uFoHmeMeRGIAMpbax91Oo+nGWP2AxHWWrd/k4+3TNQNgT3W2m+stWeB2cBjDmfyOGttMnDM6RwlxVp72Fq78cLtU5z/GZw3O5vKs+x5py/cLX3hl/PTkYcZY0KANsBEp7P4A28p6puBA7+4fxA//w8c6IwxVYG6wDpnk3jehS2AzUAW8KW11u/XDIwFBgIFTgcpQRZINMakGmN6ufPA3lLU5hJ/5vdTR6Ayxvwf4FPgBWvtSafzeJq1Nt9aGw6EAA2NMX69zWWMeRTIstamOp2lhN1rra0HtAaiLmxtuoW3FPVB4JZf3A8BvnMoi3jQhX3aT4EZ1tq5TucpSdba48AKoJXDUTztXqDdhT3b2cCDxpjpzkbyPGvtdxd+zwLmcX5L1y28pag3ANWMMbcbY64GngQ+cziTuNmFD9YmARnW2tFO5ykJxpjKxpgKF25fAzwE7HA2lWdZa4dYa0OstVU5/395mbX2KYdjeZQxpuyFD8gxxpQFWgBuO5vLK4raWpsH9Aa+4PwHTHOstenOpvI8Y8wsYA0Qaow5aIzp4XQmD7sXeJrzE9bmC78ecTqUh90ILDfGbOH8QPKltTYgTlcLMNcDq4wxacB6IMFau8RdB/eK0/NEROT3ecVELSIiv09FLSLi5VTUIiJeTkUtIuLlVNQiIl5ORS0i4uVU1CIiXu7/AenaoX7Zu4r3AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "def lineal(x, m, a):\n", " x=np.array(x)\n", " y= m*x +a\n", " return y\n", "\n", "x_lineal=[0,5]\n", "y_lineal=lineal(x_lineal, z[0], z[1])\n", "plt.plot(x_lineal, y_lineal, 'k')\n", "plt.plot(x, y, '*')\n", "plt.grid(alpha=0.3)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Se obtienen datos del espectro (velocidad y temperatura)\n", "v,T = np.genfromtxt('sdf_111_111', unpack = True, skip_header=108)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 0, 'Vel [km/s]')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3xb5fX48c+RZcl7xSPDzl6EEEIWFMKmjLLaMkuh0NJS2tLSRSnw7aCU0n5p+X1LC6W0FCgjzFBGCiSMACFk7z2cxInjeA95ypae3x/3SrET25ETS7Kt8369/LJ0daV7bElHR+c+97lijEEppVTscEQ7AKWUUpGliV8ppWKMJn6llIoxmviVUirGaOJXSqkYo4lfKaVijCZ+pZSKMZr41YAlIrtFpElE6tv9/DUKcfxaRJ6N9HaV6ooz2gEoFWaXGmPei3YQ3RERAcQY4492LCo2aMWvYo6IjBWRj0SkVkQqROTFdrcZEfmBiBTatz0oIo52t39DRDaLSLWIvCsiI9rddryILBCRKhEpFZG7ReRC4G7gGvsbx1p73YUicr+IfAo0AqNFZKiIvGHff4eIfCty/xUVS7TiV7HoPmA+cDbgAmYccvuX7GUpwHvAVuCfIvJFrCR+KbAd+DkwBzhVRFLtdf9o3x4PTDLGLBWR3wFjjTHXH7KdG4CL7McX+/4bgaHARGCBiBQaY97vxb9dKa341YD3HxGpaffzLaAVGAEMNcY0G2MWHXKfPxhjqowxRcD/AV+xl38beMAYs9kY0wb8DphqV/2XAAeMMX+yH9NjjFl6hNieMsZstB9rMDAbuNO+/xrgn1gfDkr1Kk38aqD7ojEmo93PP4CfYVXYy0Rko4h845D77G13eQ9WBQ7Wh8WfAx8iQJX9OMOAAmBnD2Nrv52hQJUxxnPItof18DGVOiJN/CrmGGMOGGO+ZYwZilXFPyoiY9utUtDu8nBgv315L/DtQz5IEo0xi+3bxnS1yRCW7wey7JZR+20Xh/hnKRUyTfwq5ojIVSKSb1+txkrAvnar3CEimSJSANwOBHb+PgbcJSLH24+TLiJX2be9BQwWkR+KiFtEUkXkZPu2UmBk+53EhzLG7AUWAw+ISIKITAFuBp479r9YqY408auB7s1DxvG/BswElopIPfAGcLsxZle7+7wOrATWAPOAJwCMMa8BfwBeEJE6YAPWzlnsFs3nsXbsHsDa+Xu2/Xgv278rRWRVN7F+BRiJVf2/BvzKGLPgWP54pTojeiIWpQ4SEQOMM8bsiHYsSoWLVvxKKRVjNPErpVSM0VaPUkrFGK34lVIqxvSLKRuys7PNyJEjox2GUkr1KytXrqwwxuQcurxfJP6RI0eyYsWKaIehlFL9iojs6Wy5tnqUUirGaOJXSqkYo4lfKaVijCZ+pZSKMZr4lVIqxmjiV0qpGKOJXymlYowmfqUGkHc2HKCktinaYag+ThO/UgNEmaeZW59dye0vrIl2KKqP08Sv1ACxeEclAK0+f5QjUX2dJn6lBohPtlcAUJCZFOVIVF8XtsQvIgUi8qGIbBaRjSJyu708S0QWiMh2+3dmuGJQKlYYY/h0h5X4G1raohyN6uvCWfG3AT8xxhwHnAJ8T0QmAT8H3jfGjAPet68rpY7BrooGDtQ1A+Bp1sSvuhe2xG+MKTHGrLIve4DNwDDgcuBpe7WngS+GKwalYsXmEg8AwzISqWtujXI0qq+LSI9fREYCJwFLgTxjTAlYHw5AbiRiUGog21bqwSEwtSBDK351RGFP/CKSArwK/NAYU9eD+90iIitEZEV5eXn4AlRqANhe5mF4VhI5qW48WvGrIwhr4heReKyk/5wxZq69uFREhti3DwHKOruvMeZxY8wMY8yMnJzDTiCjlGpn6wEP4/NSSU1w4mlpw+/Xc2mrroVzVI8ATwCbjTEPtbvpDeBG+/KNwOvhikGpWNDS5mN3ZWMw8RsDDV5t96iuhfPUi6cBNwDrRSRwKOHdwO+Bl0TkZqAIuCqMMSg14O2qaMDnN4zLS6HR6wOskT2pCfFRjkz1VWFL/MaYRYB0cfO54dquUrFme2k9AOPzUtlZbl3WHbyqO3rkrlL9XJmnBYAh6Qmk2VW+7uBV3dHEr1Q/V9tkJfnUhHhSE6wv8Vrxq+5o4leqn6tt9JKW4CTOIcG+vh7EpbqjiV+pfq62qZX0JCvhp2nFr0KgiV+pfq62qZWMRBeAVvwqJJr4lernappaSU+0En5CvAOnQ7TiV93SxK9UP9e+1SMipCXG66ge1S1N/Er1c7WNByt+wJq2QSt+1Q1N/Er1Y8YYu8d/MPFnJMZT1eCNYlSqr9PEr1Q/1uD10eY3HSr+/Kwk9lY1RjEq1ddp4leqHwscvJWRdDDxj8hKYl91E2160nXVBU38SvVjtY1W4m9f8Y8clEyb37C/pjlaYak+ThO/Uv1YTZPVy0+3x/EDDB+UBMCeqoaoxKT6Pk38SvVjdU2HV/wjAom/Uvv8qnOa+JXqx2oaD+/x56Um4HI6KNIdvKoLmviV6sdqO6n4HQ5heFYSeyq11aM6p4lfqX6spqmV+DghyRXXYfmIrCRt9aguaeJXqh+rafSSnujCOsX1QcMHJVFU1YgxetJ1dThN/Er1Y+WeFnJT3YctH5GVRKPXR3l9SxSiUn2dJn6l+rEyTws5nSX+QckAFGm7R3VCE79S/VhZXRcVvw7pVN3QxK9UP+X3GyrqO6/48zOTcAjs0SGdqhOa+JXqp6obvbT5TacVv8vpYEh6IkU6pFN1QhO/Uv1UYMdtTmpCp7ePGJSkFb/qlCZ+pfqpsjor8eemHV7xg5X4deeu6owmfqX6qTKPXfGndJ74h2clU9ng1dMwqsNo4leqnyr3dF/xD82wWkCldTo9s+pIE79S/VSZp5kUt5Mkl7PT23Pt3n+gJaRUgCZ+pfqp8i4O3grIs78JlHq04lcdaeJXqp8q97R02d8HyE0LtHq04lcdaeJXqp+qbWrtMA//oVLcTlLcTu3xq8No4leqn/I0t5Ga0HXiB8hNdQdH/ygVoIlfqX6qrqmV1ITOd+wG5Ka5KdOKXx1CE79S/ZDfb6j3tpGW2H3Fn5eWoD1+dRhN/Er1Q/XeNoyBtCNU/Fbib9YTsqgONPEr1Q/V2efaPWKrJ9VNS5ufuua2SISl+glN/Er1Qx47kacdaeduWuAgLu3zq4PClvhF5F8iUiYiG9ot+7WIFIvIGvvnC+HavlIDWSDxhzKqB3Qsv+oonBX/U8CFnSz/f8aYqfbPf8O4faUGrFBbPYFvBDpRm2ovbInfGPMxUBWux1cqlnlarER+pFE9KW7rg6G+RXv86qBo9PhvE5F1disoMwrbV6rfO9jq6b7iT3bHAdCgiV+1E+nE/zdgDDAVKAH+1NWKInKLiKwQkRXl5eWRik+pfiHUVk+yXfE3eH1hj0n1HxFN/MaYUmOMzxjjB/4BzOpm3ceNMTOMMTNycnIiF6RS/YCnuQ2304HbGdftem6nA6dDtNWjOoho4heRIe2ufgnY0NW6Sqmu1YUwTw+AiJDsdmqrR3XQ/ffEYyAic4CzgGwR2Qf8CjhLRKYCBtgNfDtc21dqIKtrbj3iUbsBKW6nVvyqg7AlfmPMVzpZ/ES4tqdULLFm5gzt7ZvsjtOKX3WgR+4q1Q95mluPOJQzwGr16M5ddZAmfqX6oZ5U/NrqUYfq8pUjIg+FcP86Y8yvey8cpVQo6ppaSXWHWPG7nByo1bl61EHdlQxXAL85wv1/Cvy616JRSoWkZz1+HdWjOurulfMXY0y3O2P1yFulIq+51UdTq4/MZFdI66e447TVozrorsc/p6sbROQiAGPMH3s9IqVUtyobvAAMCjHxJ7udNHh9ejIWFdRd4n9PRIYfulBEvgY8Er6QlFLdqbYTf1aoFX+CE5/f0NLmD2dYqh/pLvH/DCv5jw4sEJE7gDuxDsxSSkVBZU8Tv87QqQ7RZY/fGPOmiLQA74rI5cA3gNnAGcaYykgFqJTqqKrBOqlKqIk/2WVP1NbSRnaKO2xxqf6j23H8xpj5wLeAj4HjgLM16SsVXVUN1sycg5JDS+LJWvGrQ3Q3jr8aa04dAZKwqv1iERHAGGOyIhOiUqq9qoYW4hzSowO4AD16VwV198rJjlgUSqmQVTV4yUxy4XBISOvryVjUobrr8Wt5oFQfVNXgDXkoJ+jOXXW4Lnv8IrLsSHcOZR2lVO+qavCSmRzadA3Q7ixcmviVrbtWzwkisqqb2wUY1MvxKKWOoLLBy3GD00JeX3fuqkN1l/gnh3B/fSUpFWHVDd6Qh3KC1eoRsc7apRR03+PfGclAlFJH5vMbappaQ56nByDOIaQlxFPT6A1jZKo/0fn4lepHqhu9GBP6PD0BWckuqhtbwxSV6m808SvVjwTm6elJxQ+QkaQVvzpIE79S/Uhds1W1p4d42sWAzCQXVQ2a+JXliIlfRGaKyBIRqRWRZhFpEZG6SASnlOoosIM21KN2A6yKX1s9yhLKq+dR4HrgBWAWcBNQEMaYlFJdqGuykndaQs8q/qwkF9Xa6lG2UFo9DmPMVsBpjGk1xvwDOC/McSmlOhGo+NMSe1bxZya7aPT6aG7VA/JVaBV/g4i4gLUi8jugBEgJb1hKqc4cbcWfkWStX9PYyuD0uF6PS/UvoVT8N9nr3Qb4gHHAlWGMSSnVBU9zG644BwnxPUvemUnWKCBt9yg4QsUvInHAr4wxNwLNwC8iEpVSqlN1za09bvOAJn7V0ZFOxOIDhohIz75XKqXCoq6ptcdtHiA4qVt1g47sUaH1+AuBT0TkdaAhsNAY83DYolJKdaquuY3UHo7hB634VUehJP5yYAHWWbiSwhuOUqo7VsXf81bPwZ27mvhVCInfGKN9faX6CE9zK8MyEnt8P7czjmRXnM7Xo4AQEr+ILMA6924HxpjzwxKRUqpLdc1tR7VzFyAjyRWc60fFtlBeQf/T7nICcAXQEp5wlFLdOdqduwBpifE6J78CQmv1LD1k0Uci8lGY4lFKdaG51UdLm7/H8/QEpLjjqG/RVo8KrdXT/hxvDmA6MCRsESmlOuUJTtdwdBV/ittJRb22elRorZ6NWD1+wTrV4i7gW+EMSil1OE/z0U3XEJDsdrKnsrE3Q1L9VCiJf7QxpsP3QxE5uu+aSqmjdrQTtAWkuJ16wnUFhDZXz6E9foBlvR2IUqp7gQnaUo+h4m/QxK/opuIXkVysXn6iiJyA1eoBSCOEA7lE5F/AJUCZMWayvSwLeBEYCewGrjbGVB9D/ErFjNpg4j+6ij/Z7aTB68PvNzgccuQ7qAGru4r/YuCvQD7WyVgesX/uJrTJ2p4CLjxk2c+B940x44D37etKqRAU1zQBMCS95wdwgTWqB6BR5+SPeV2WDsaYJ4EnReRqY8xLPX1gY8zHIjLykMWXA2fZl58GFgJ39vSxlYpFuysaGJTs6vH5dgOS3dbbvb65jRS37qaLZaGM439JRC4Ajsc6gCuw/HdHsb08Y0yJff8Su52klArBrooGRmYnH/X9A8led/CqUE62/ihwI/BjIBHr/LtjwxwXInKLiKwQkRXl5eXh3pxSfd6uigZG9ULi1x28KpRRPbONMdcBlfaEbSdj9f2PRqmIDAGwf5d1taIx5nFjzAxjzIycnJyj3JxSA0NDSxtlnpZjSvzJmviVLZTE3xz4LSKD7esjj3J7b2B9e8D+/fpRPo5SMWV3pXUqjN6o+LXVo0JJ/P8VkQzgj8AarGGYrxzpTiIyB/gMmCAi+0TkZuD3wOdFZDvwefu6UuoIdlVYiX/koF6o+L2a+GPdkc656wDeNsbUAC+LyFtAojGm6kgPbIz5Shc3ndvzMJWKbbsDiT/76M+FlGwP56xv0eGcse5I59z1A39ud70plKSvlOpdxTXNZKe4SHId/TDMlHbDOVVsC6XVs0BELg97JEqpLpXVNZOTmnDkFbuRGB+HQ3TnrgptkrbbgHQRaQGasKZuMMaYrLBGppQKKvO0kJfmPqbHEBGSdaI2RWiJPzvsUSilulXmaea4IanH/DgpOlGbIoRWjzHGB1wF3GlfHgJMDXdgSimLz28o97SQe4ytHghM1KaJP9aFcuTuX4GzgRvsRY3AY+EMSil1UGVDC34DucfY6gHsVo+O6ol1obR6TjXGTBOR1QDGmCoRcYU5LqWUrayuBaBXKv4Ud5y2elRIo3pa7fH8BkBEBgH+sEallAoq81gHz/dGxZ/idupwThVS4n8EeBXIEZF7gUXAH8IalVIq6GDF3xuJPz547l4Vu0KZlvnfIrISOM9edJUxZkN4w1JKBZR5rMSf0wuJPyMpPngmLxW7Qqn4AeKAVsDbg/sopXpBaV0zmUnxuJ1xx/xYGYnxNHh9eNu0WxvLQhnVcw8wBxiKNR3z8yJyV7gDU0pZynppKCdYFT+gVX+MC2VUz/XAdGNMI4CI3A+sBB4IZ2BKKUtVg5dBKb0zkC49yXqc2iZvr7SOVP8USttmDx0/IJxAYXjCUUodqqbRG6zUj1WGfb7emkat+GNZKBV/I7BRRN7FGtJ5PrBIRB4CMMb8OIzxKRXzaptaSU/snYo/8AGiiT+2hZL459k/AUvCFItS6hDGGGqbWnux4rc+QGq0xx/TQhnO+UQkAlFKHa7R66PVZ0hP7J3Enx6s+L298niqfwplVM+FIrJcRMpEpEpEqkVET8aiVAQEKvOMXkr8qW4nDtFRPbEulFbPX4GrgfXoVA1KRVSt3YvvrVaPwyGkJ8Zrjz/GhZL49wFr7NMwKqUiqKbJasmk9VLFD5CR5NIef4wLJfH/DHhTRBYCLYGFxpiHwxWUUspSF2z19N6EuFbFrz3+WBZK4r8Xa7qGDLTVo1REBVoy6b3U6gGrbVTVoIk/loWS+HONMdPDHolS6jC9vXM38FiF5Q299niq/wnlyN33ReScsEeilDpMbVMr8XFCkuvYJ2gLyEhy6aieGBdK4v8W8J6I1OtwTqUiq6axlfTEeESk1x4zLTGeuuZWfH7Ta4+p+pdQWj3ZYY9CKdWpuqbWXjt4KyAzKR5jrG8TWcl6FtVYdMSK3xjjA64C7rQvDwGmhjswpZQ1nDMjqXeTcyDZV+vInpgVypG7fwXOBm6wFzUCj4UzKKWUJdDq6U2Z9geJjuyJXaH0+E81xnwbaAYwxlQB+v1QqQiobWrt1RE9cLDi18Qfu0JJ/K0i4sCakhkRGYSO51cqIqobvL06hh/atXo08cesLhO/iAR2/D4CvArkiMi9wCLgDxGITamYVt/SRoPXR15a75x2MSDY6tEef8zqblTPMmCaMebfIrISOA8Q4CpjzIaIRKdUDCutawYgL613T5GY6IojMT5OK/4Y1l3iDw4cNsZsBDaGPxylVEBpbSDx927FD1a7p6pBD+KKVd0l/hwR6fK0isaYh8IQj1LKVuoJX+LPTI6nqqHlyCuqAam7xB8HpNCu8ldKRU5pnZWYw5L4k1xU6Zz8Mau7xF9ijPlNxCJRSnVwoLaZFLeTFHcoB9j3zKBkF3sqG3v9cVX/0N1wTq30lYqiMk9zr+/YDchMdunO3RjWXSlxbrg2KiK7AQ/gA9qMMTPCtS2l+qsDtc1hafMAZCW58LS04W3z43KGcjiPGki6fMbtI3TD6WxjzFRN+kp1rrSuhcFhSvyZ9kFceiau2KQf9Ur1QX6/oczTTG64Kn478VdquycmRSvxG2C+iKwUkVuiFINSfVZ1o5dWn2FwuHr8STptQyzr/eECoTnNGLNfRHKBBSKyxRjzcfsV7A+EWwCGDx8ejRiVipoyjzWUMyc1PBV/WqL11q9rbgvL46u+LSoVvzFmv/27DHgNmNXJOo8bY2YYY2bk5OREOkSloqqy3qrEs1PCMxFuWoI18ZunWcfyx6KIJ34RSRaR1MBl4HxA5/5Rqp2Keqviz04NT6snkPi14o9N0Wj15AGv2ecQdQLPG2PeiUIcSvVZwcSfEp7En5Jgt3r0pOsxKeKJ3xhTCJwY6e0q1Z9U1HtxxTlISwjPWzTOIaS6ndRpqycm6XBOpfqgivoWBqW4sL8Zh0VqghOPtnpikiZ+pfqgivqWsLV5AtIS47XVE6M08SvVBwUq/nBKS4jXVk+M0sSvVB9UWe8Ne8WvrZ7YpYlfqT7GGBORxJ+WqBV/rNLEr1QfU9fUhtfnD9vBWwFpCU7qmrTij0Wa+JXqYyoawjuGPyAtMR5PcyvGmLBuR/U9mviV6mMqPJFJ/KkJTvwGGry+sG5H9T2a+JXqYwITtGWnhn9UD+jRu7FIE79SfUxRlXUu3ILMpLBuJy0xMF+PJv5Yo4lfqT5mT2UD2SluksNwkvX2Uu3pIHRIZ+zRxK9UH7OnspERg8Jb7YO2emKZJn6l+piiqkZGZEUg8WurJ2Zp4leqD2lu9XGgrpnhEan4rVZPbaMm/lgTrVMvKqUO8cd3t7J0VyXGEJFWT3piPCJQpYk/5mjiV6qPmL/pANtK6wEYMSg57NtzxjnISnJRbg8fVbFDWz1K9QFtPj+7KxqD1yPR4wfISXUHz/alYodW/Er1AXurm/D6/Jw1IQeArOTwHrwVkJ2iiT8WaeJXqg/YWWa1eH5w7jimDc+M2HZzUt3s3t0Qse2pvkFbPUr1ATvLrcQ/JjslotvNTnFRUd+iE7XFGE38SvUBO8vryU5xk54UH9Ht5qS6aW71U9+iR+/GEk38SvUBO8sbGJMT/pE8hwrMAFpR7434tlX0aOJXKsqMMewoq2dMbmTbPNA+8esO3liiiV+pKKtq8FLb1MrYnMgn/pxUK/HrWP7YoolfqSjbYY/o0YpfRYomfqWibGe5NZwyGj3+rGQXDtGKP9Zo4lcqynaW15MQ72BoemLEtx3nEAaluJmzbC+PfLgj4ttX0aGJX6ko21lez+jsFBwOicr27zh/AplJ8Tz64Q7afP6oxKAiSxO/UlG2s7yesVHo7wdcPbOA284ZS4PXx5YDnqjFoSJHE79SUdTc6mNfdRNjojCip72ZI7MAWL67KqpxqMjQxK9UFO2qaMAYGJMb+R277Q3NSGRYRiIrdldHNQ4VGZr4lYqi4FDOKFf8ADNHZrJsd5XO2xMDNPErFUU7y+sRgVHZ0a34AaaNyKTc00JJbXNw2dxV+3h68e7oBaXCQhO/UlG0s7yBgswkEuLjoh0Kk4elA7C+uDa47JEPd/DH+Vt1tM8Ao4k/wqoavGw5UBftMFQP1Ta14m3r/eS3s6w+KgdudWbSkDTiHMKK3VV88+kVLNxaxs7yBjzNbR0+DPqzhpY27n5tPWV1zUdeeQAb8Im/0du3ppv9w9tbuOqxz/pFBbW0sJLnlu6Jdhi97unFu1laWNnpbX6/6dDjNsbw5Uc/5cR75/PNf68ILttQXHtMvfAlhZV87/lVFFbU94n+PkBCfBzjclN4Zske3ttcyk9fXhe8bfHOzv9f/c0HW8p4fmkRb60r6fF9P9leHtwn098N6MT/4LtbuPJvn9EQwlzjq4qqKapsPOJ6x2r5nio8zW1sOeChtYvk3+rzc/0/l/LRtvIeP35RZSM7yg6Oxe5qG5WdzM1ijGHh1rLgff7xSSH3vbUJv//oE1xheT3VDX1nyt9Gbxu/eWsTT366GwjMjHnw/3XbnFV8f87q4PWKei+rimrIz0zk423lbD3g4dkle7jkL4v492cdPxTb/592lHkormnqMo53Nhxg3roSmlv9UZmjpyvHD02nudV6/ivqW3A6hNE5ySzaXnHE+364tYz9h/zNxhje2VByTHMBNXl9/GvRLlrafCGtv3ZvDTc9uazTcwwE/o7Ve2t6FMPeqka+8dRyHlqwtUf366uikvhF5EIR2SoiO0Tk5+HazoyRWWw5UMePXlyDMYYnFu3id//dzA1PLOVHL64Jrvff9SVc+bfF/PSVtZ0+Tl1zK7WNrcccT02jl0J7XpbnlhYx+VfvsnjH4W+oTfvrWLSjgg+3lPHYRzu57flVXT7mrooGbvn3Cqrs5Pqd51Zy47+W4/cbyj0tnHjvfJ75bDd3v7aeB97eDFhjtWfe/x6vrd7X4bHeXFfCTU8u57GFOwHYWuqhudXP3upGNhTXctGfP+H9zaXB9Y0xrOhmFEhNo5fL/vopv3xjY+j/pF7Q6vPz4xfXdFrVr9lbg89v2G4n+/c2l3HeQx+ztLCSktom3t5wgKW7Do5lD6z3swsn4nI6+NUbG/jdf7cA1gdj4JvbOxtKmPiLd3howTZ8fsPXnljGpX9Z1GWFuLfqYJExPi+1d/7wXnDCsDQAzp2YC8Dxw9I5f9Jgluyq5K656/j1GxspqT38A62qwcvNTy3nO8+uxNfuA/ChBdu49dlVPPjOwYRZ7mnhmr9/1uEDtzuvrynmN29t4tWVxZ3evnhnBTc/tZwbnljKYx/t5McvrWHh1nKW7apke6knWPgZY1hkv99WF/Vs2Oof52+l1WfYXdFIYXk9//ykMKyjn6oavLy0Ym/YthHxxC8iccAjwEXAJOArIjIpHNs6e0IuP71gAvM3lTJvfQn3vbWJf35SyOqiGl5bXUxpXTMV9S388MU1xMc5WLmnmprGjtXpPz4u5OT73+f6J5bS5vPz/ubSbp+MB/67mTteXsuLy4vYUVbfYd1AlSECc5YV0dLm57llRYc9xoo91ouysKKBdzceYN76kk4/eIwx/PzVdcG/b+sBDxv311Fc08SSXZWs3FNFo9fHL17fyPNLi3huSREtbT7ufXMjfgN/W7gTYwytPj8ltU385f3tAPz940L2VjWyt8p6g3+4pYxr/v4Zm0vqmLv64JvvxeV7ufKxz7r8ZvLkp7upb2nj/c2lNLeGVq2FYl91I498uKNDhe3zm2DCWbarirmri/ne86so83Ts5QbGqe+pbKTV52fdPus5eW5pEXNXFWOMlZgCr4NA4j55VBZXTMtnSWEVQzISuO/y49lXbX1QACwprMLr8/Pw+9v528Id7K9tpqrB2+WHdlFVI+dOzOX5b53MtOEZvfa/OVbnHpfHKaOz+N2XT2D22GwunTKE758zlmtnFjBn2V6eWrybuaus18Cm/XXc8MRSLn74E+YsK8JvYO2+Wp5dYn0TenXlPk+iCLAAABqmSURBVP7ywQ6SXHF8sLUs+Hwt21XF0l1V/P5t6wN0b1Ujm/Yf3O/l9xveXLufbzy1nF+/sZFP7Cr9qcW7DnvvLd5ZwXX/WMr64loq6738/u0t7CxvQAQWba/ksr9+yveeX8XKPVXcP28zxTVNjM5JZl91Ew/N38qrKzsWPwG1jVax52lu5ZtPL+f1NftJdsWxp7KBZ5cU8dt5m9lVEb5zFf/jk0J+9so6VhX17JtJqKJxsvVZwA5jTCGAiLwAXA5sCsfGrp05nAff3coDdpX23o/PpKXNz0V//oSPtpbjaWnD2+bnD1ecwJ2vruejbeVcNHkILqeDljYfD87fitMhrC+u5Zkle7j3zU08ct00dlc24HY6+OrJI0h0WSMy6ppb+fvHhTgdwsv2Cyovzc0vLpnERZOHsHpPNQ6Bsybk8sGWMlxxDhZsKuWFZUVsL6tnwuBUrp5RwMo9VsW5q6KemoZWjLGq9PMm5QGweEcFia44VhfVsHRXFfFxwvyNB9hX3UicQ0hwOnhtVTGDUtzExwm5qQkku+PYVlrPb9/azIbiOs47Lpf3Npfx0bZyPtlewROLdgFw+7nj+MsH27n7tfXB/+FjHxXS4PUxfUQmSwur8PsNtU2tPGx/UCzcWs5ZE3I7/N8bvW08+eku8jMT2VfdxHefW0Vzq49/f2MWzriD9UZzq49yTwsFWUkhP6ePLtzJ80uLODE/g9njsgH4zZsbWV9cy9zvnsaCTaW4nA48zW385KW1PP31WcF5cAIfqm1+w57KBrbaUxS8s+EASworSYyPo6nVx46yemaMzGJ7aT2pCU5yU938+rJJ3H7uOPLS3BgDTy7ezd8/3sklU4awvczDCcPSKalt4lH7G9P1pwzn2SVFFFU2MnzQwb/PGENRVSNnjs/h1DHZIf/dkVCQlcQLt3wOgGe/eXJw+QNfnsK9l03moj9/zCr7f/jH+VtZu7eGRq+PzSV1DElPYFxeKr/772ZcTgf3vbWJk0dlceX0fO54ZR1Pf7abIekJHb5t/e87W3hmyR7i4xy8fOvnuOnJZfh8hv21zaS4nXywpYzE+DgGJbvYVlrPksIqPjdmEC8uL+LlFfso87RQkJXIuz88gySXk6WFlRyoa+bxjwuZs6yIplYfC7dar3Gf3xDnEL5/zlh+9OJaHv7AmpQuziF88aRhAHiaW7n9hTV8uLWM3FQ3Z0/I5f0tZdx10UTiHMJv521m8U7rg+jjbeWM7mL/TG1jK2mJTvwG5m88QF1zK9fMHA7AKyv3UVrXzMKtZTS0+Hj9ttOoqG8hJ8UdfG98uKUMgDfWFDN9RGZvPsVAdBL/MGBvu+v7gJO7WPeYZSW7OKkgg1VFNRRkJQbHS+elufloWzm7Kho4YVg6V00v4MF3t3Lnq+v45esb+fiOs9lW5sHb5ufG00fxj092BWcvfODtzeyrtqrhv39cyF0XTeTL0/LZYI98+OeNM8jPTGLlniqeW1rEbc+vxhW3FoNhwuA0zploJf5fXTaJe17bwM/nrsftdNDS5rfbJ9YbK1BxAyzdVcl5k/Jo9fm55ZmVNHjbMAbOOy6XUdnJPLV4N+uLazlzfA7ZKS7mrSthdE4Kk4ak8cp3TqWhpY2T7lvAM0v2MDY3hUe+Oo3Tfv8hc5YVsX5fLVPy07nsxKF8/bRRrCqqDlZZLqeDA3XNjM5O5srp+dw1dz3XPP4Zy+0YB6cl8NG2ch6av5X9tc1kJbtIT4xnWEYidc1tPHb9dL77/Co+sF/In+6s5MzxOQB42/zc9OQylhRWcemJQ/nTVSficlov/G2lHvbXNB32geJt8zPP3jH3ysq9zB6XTZvPzxtr91PT1EpDSxsLNpVyxrhszj0uj7vmruexj3fy3bPG4vMbVu+p5sSCDNburWF7aT1bSz1MHJzKjrJ6UtxO/ueSSfxgzmq224l/W6mHcbkpiAhuZxyD060PeRH41umjuWvuej7bWcn20nrOGJ/D5GFpzFm2l+wUN18/bRTPLinio+3l3DBoRPBvKPe00NLm7/Bh0B+4nA5OGp7JB1vKaGhpY9GOCq4/eQQtbT6eW1rE5yflcds5Y7n44UXcNXc9BVmJ/Pnak3A5HTgE7n1zE26ngzPH5zAkPYFBKS4eXbiTZFccVc1efjfPel9dMGkwPzk/j7Mm5HDq7z+gqdXHry+bxC9e38gHW0oxGO5+bQOuOAdNrT6evGkmSS4rlZ08ehBgFUob99eRnhjP+LwUvG1+HvnqNHx+Q15aAr/4z0ZOHpWFp7mNX76+gYunDKG+uY0bn1zGpv11fOO0UTy9eDcvLN/LldPz+faZY/hgi9XmDMxn9Mn2Cm46bZR9uZyiqka+MnM4a/bVcM3fP2PSkDSqGr3B9/GYnBSGZiTy05etlnJ6Yjy1Ta388MU1zFtXQqrbyf1fPoHpIzLZcsCDy+ngrXUl/OKSSR2Kpd4QjcTf2RSEh/VOROQW4BaA4cOHH9MGz56Qy6qiGs4cn4OItfkzx+cwd1UxbX7DvZcdj8MhXDEtn7mriyn3tPDOxpLgeUhvnj2apxfvoaLei9vpYF91E3lpbv501VT+990t/PTltVxw/GDW77MS/5T8DLKSXYzNTeGKafnMW1/C5hIPza0+zjsuj1mjsjgxP4PJw9KoqvcyIjuZiyYP5htPLefOV61Ke9rwjODXvGRXHEsKrW8B6/bVUN/SxsmjsijISuL+L01mQ3Et//hkF26ng/+5+DjqW9p4acU+1hfX8rXPjSA+zkFGkospw9JZu6+WW84YjdsZx6UnDuGpxbsxBm4/b1ywIrl0ylA+2V5BQryDmSOz+GR7BWeMz2HWqMB8LtVceuJQZo8dRH2Lj/ve2sTDH+wgNcFJk9dHm9+Qk+pmWEYip4wexO3njmNPZSOvrtrHqyv3UdPoZXVRDSv3VLO+uJbzJ+Xx5tr9XHzCEC6cPJh560r48UtraPX5eeU7p/LuxgPccMoIEuPjeG11MbVNrYzLTeHtDQc4f30JyW4n1XYrbO7qYoprmvj+OWO5ZmYBH2wp49EPd/L1U0exuqgaT0sb1588nLV7a1hfXEtRVSM/PHc8Xzm5gEHJbgS485V1vLpyH6+s3Me2Ax4unjKk09fVl04axp/mb+P/vbeNMk8L43JTGD84lTnL9nLy6CxGZydTkJXIvHX7cTsdXD51KG5nHEV2f78n33L6imnDM3ll5T6eWbIHb5ufz0/KY0xOMptL6rh6RgG5qQk89fWZrCqq4arp+cHjE75wwhA27a+jsKKBhVvLOX1cNk/cNJPSumaMgVMeeJ/3t5QxfUQmj90wPbi9L08bxovL93LOxDxeWbmP5burWbqrivzMRN64bTZ1Ta2d/h9PKsjk2SVFnD0hhz9dPRWHEHzvA7z/kzMZlOxiwaZSvvPcKj7ZXs4f3t7KrsoG/n7DdM49Lo/cVDf/XLSLn5w/HoDhWQeH3SbGx/FZYSXeNj8up4M/zt/G2r01vLPhAMU1TWQmuahrbmNwWgJ3XjiRe9/cxB/e2cLn7W/t834wm/F5qXzp0U+Zt66E44daQ2nvnrueq2cUANa37wff3dqhWOot0Uj8+4CCdtfzgf2HrmSMeRx4HGDGjBnHtIfjgsmDefiD7Xxh8sE38HUnj2BvVRMzRmZyzUwrnLu+cBw/v2giZ/9xIW+s3Y/T4WBcbgqD0xOYkp/Oij3V/PT8Cfzlg+387IKJzB6XTaN3LLc8s5JtpR7WFdeSn5lIVrIruB1nnIPLpw7j8qkdYzoh3zpY5vvnjgsu+9v103l68W5W7anm6pkFfPuZlQB8eVo+zy3dw8o9VXy6oxIReOz66WTa25k2PJP/vWIKnxszKPgmmDUyi2W7qzgx/2D/+LKpw2hq9XH51KHW9ROHBke3nN2usr7g+MHc85/1jMtNZeLgVDvxZzM6O5mcVDeJ8XE8eOUUEuLj2FFWz3329ubccgrGGC5+eBFbSz3ceuYYHA7h63ZV1OT18eKKvbyxdj+J8XGMH5zKfV+czLUzC5h1/3u8vaGEM8Znc/dr64NV+LV/X4LX5+ft9QeobvTiaW4jJ9XNn64+kS89upjvPrcKV5xVUfoNPPrhDkSsXrWI8M3Zo1iwqZT/ri/h050VpLqdXHriUP7vve28veEAxsCEwSnkpiYE//4xucnBlhDQ5cyZCfFxXDk9n8c+slo74/JSOHXMIGaMyOSLU4chIpw5PodnlxSxpLCK/TVN/PC88cHEP7w/Jv4R1uvp0Q93kJ4Yz8yRmTjjHMz97mnBdY4fms7xQ9M73O+v102jor6FGb99D6/Pz9g863+al2b93ycNSWNTSV0wMQbcc/EkrpyeT06qmxkjs3j840J8fsOdF04kPTGe9MT4TuM8eXQWLqf13ovrZLrrwHZPG5dNnEP42SvrqKj38szNszh9nJVkv33mGL55+ujg/QuyEhEBY+C6k4cHB4v84pJJbDvgYUJeKhuKa6lpauVfN83s8J6qbWrlntc2sKG4jin5B/8/d190HA+8vYVHrpsGwAX/9zH/+nQXJ+an883TR5EQH8ekIWmhPDU9Eo3EvxwYJyKjgGLgWuC6cG5wfF4qa355Psnug3/u1IIM5txyymHrigiXnTiUv364A5fTwRXT8gH43JhBrNlbw5enDePm2aOCPePj7CdlywEPG4qtlsnRSnE7+d7ZYwGCOxezkl386PPj+WR7Od98egWZyS6OH5oWTPqBmK+eWdDhsb579hg2PFfLyaOzgstunj2Km2eP6vA/GDEoiYzEeHLTDia+9KR4fnDOOHJS3YzKTmbt3lo+NzobEeHxG6aTmhAfrOTG5qbw52uncuqYbPsNIvzq0kl89/lVXDk9v0NMXzt1BOuKa7n1zNFcOmVoh/nnz580mHnrS3h+aRG1Ta384pJJfLy9goff387XTxvJc0uLmJCXyh0XTGB8XiqD0xNY/cvP8+GWMu54eR2nj8th4/5aSmqbOWl4RvBcsrNGZTFyUBL/+nQXheUNfPGkYSTExzElPz24Y3bC4I5vrHG5qWworuPOCyeyfHfVYe2m9r48bdjBxJ+bitsZxyvfOTV4+7dOH016YjybSzz8beFOrpiWT1FVIyIwLCPyJ145VuNyUxmU7KLR6+OuL0zoUQsiO8XN6JxkCssbGJfbcSTT6eOz2VRSx3nHdUz8KW4n00dYr+FZI7P4m73/5MLJg7vdVn5mEmt/eX5w/1tX0hLimT7cmqPovONyg0k/oP2HhtsZx5C0BPbXNvPtM0cD8MSiXaS4nTS1+rj59FFcMS2fck8Lg9MTOjzOdbOGs6G4ljnL9nLZiUODy08dm82b358dvP7vm2dR3eDljPE5uJ1xHd6vvSniid8Y0yYitwHvAnHAv4wxYR/v1z7pH8lVMwp4f0sZfmN9nQe49cwxXDh5MIPsc5QGDMtIJMXt5LOdleypbAx+ezhWGUkuspJdjMtNISvZxVNfn8VNTy6jsLyBb58x+oj3P2tCLht+fUG3J/cQEZ68aSZOx+Fv3vbfRF669XPByycNP3xH0+VTh3W4furYbNb88vzD1jt+aDpv3356p7FcPGUIL67YywNvb+HE/HSmj8jkpOGZXDJlCOPzUvnuWWPJSIonvl2iSUuI5/Kpw5g8LJ3UBCd3vLyOj7aVd6gaRYTrTxnBb+dtRgSutZ+f31w+mfqWNnZVNBxWeV8zs4C8tARuPXM03zlrTKfxBozPS2XysDR2ljV0mshHDErmjgsmWvsrHlzIk5/upqqhhcFpCX1imoaeinMIb/1gNknxTtKTOq+2uzNzRJad+Dt+i7r1jDGcVJDZ7XkJpo3IRAQm5KWGNLfRkZJ+wNkTc1m2u4rvnDX2iOsOH5REc5ufnBQ3/3Pxcbyz4QBPfmoNjJg4OJU4hxyW9MF6Hf72iyfw+Ul5zB7bddsmMD12uEl/mIlvxowZZsWKFdEOo0tX/G0xq4qqMQbe+v7s4Jwnx+rpxdYoiPOPt6qbRm8bLy3fy8VThgYr2oHCGMMLy/eyYnc1184qOKo3wIPvbuGRD3ey4EdnMK7d2HhjDBX1XgymQ0sHCI70OBbLd1exq7zhsG9dh/r2MytYsbua5lYfFxw/mIeumdrt+gPRou0V3P/fzcz9zqkhJ+b2fvvWJk4syODSdlXzsWpu9bFxf23wm0V35m88QKmnhRtOsXbW3/Paep5bWoRDYOO9Fx7V3xROIrLSGDPjsOWa+I9d4Mk/YVh6h69tKrLKPS18uqMiODSvr3l344Hgfpt5P5h9WB9c9T/vbSrlm/9eweicZD74yVnRDucwXSX+aPT4B5yJdp//K7OObfSROjY5qe4+m/TB2oGeleziuCGpmvQHiFPHDsIV52Di4L5z9HUoNPH3gosmD2ZPRUNwf4BSnXE5Hbxy6+e6HImi+p8kl5O/XncSI/vA+RR6Qls9Sik1QHXV6hnQs3MqpZQ6nCZ+pZSKMZr4lVIqxmjiV0qpGKOJXymlYowmfqWUijGa+JVSKsZo4ldKqRjTLw7gEpFyYM8xPkw2cPiZzaOvL8bVF2OCvhlXX4wJNK6e6IsxQe/ENcIYc9h0oP0i8fcGEVnR2RFs0dYX4+qLMUHfjKsvxgQaV0/0xZggvHFpq0cppWKMJn6llIoxsZT4H492AF3oi3H1xZigb8bVF2MCjasn+mJMEMa4YqbHr5RSyhJLFb9SSik08SulVMyJmcQvIj8VESMi2fZ1EZGHRWSHiKwTkWkRjOU+e5trRGS+iAyNdkz29h8UkS32tl8TkYx2t91lx7VVRC6IYExXichGEfGLyIxDbotKTO22f6G97R0i8vNIb79dHP8SkTIR2dBuWZaILBCR7fbvzAjHVCAiH4rIZvv5u72PxJUgIstEZK0d17328lEistSO60URcUUyLjuGOBFZLSJvhT0mY8yA/wEKgHexDgLLtpd9AXgbEOAUYGkE40lrd/kHwGPRjsne/vmA0778B+AP9uVJwFrADYwCdgJxEYrpOGACsBCY0W551GKytx9nb3M04LJjmRTJ56tdLGcA04AN7Zb9L/Bz+/LPA89lBGMaAkyzL6cC2+znLNpxCZBiX44HltrvtZeAa+3ljwHficLz+GPgeeAt+3rYYoqViv//AT8D2u/Jvhz4t7EsATJEZEgkgjHG1LW7mtwurqjFZMc13xjTZl9dAuS3i+sFY0yLMWYXsAOYFaGYNhtjtnZyU9Riss0CdhhjCo0xXuAFO6aIM8Z8DFQdsvhy4Gn78tPAFyMcU4kxZpV92QNsBob1gbiMMabevhpv/xjgHOCVaMUlIvnAxcA/7esSzpgGfOIXkcuAYmPM2kNuGgbsbXd9n70sUnHdLyJ7ga8Cv+wLMR3iG1jfPqBvxRUQ7Ziivf0jyTPGlICVhIHcaAUiIiOBk7Cq66jHZbdU1gBlwAKsb2417YqeaDyX/4dVnPrt64PCGZOztx4omkTkPWBwJzfdA9yN1cI47G6dLOu1sa3dxWSMed0Ycw9wj4jcBdwG/CrcMYUSl73OPUAb8FzgbuGMK5SYOrtbOGMKQbS33y+ISArwKvBDY0ydVchGlzHGB0y192G9htVOPGy1SMUjIpcAZcaYlSJyVmBxOGMaEInfGHNeZ8tF5ASs/u9a+wWXD6wSkVlYn6AF7VbPB/aHO6ZOPA/Mw0r8YY0plLhE5EbgEuBcYzcXwx1XD/5X7YX9f9XHt38kpSIyxBhTYrcLyyIdgIjEYyX954wxc/tKXAHGmBoRWYjV488QEaddYUf6uTwNuExEvgAkAGlY3wDCFtOAbvUYY9YbY3KNMSONMSOx3qzTjDEHgDeAr9kjaU4BagNfQcNNRMa1u3oZsMW+HLWY7LguBO4ELjPGNLa76Q3gWhFxi8goYBywLFJxdSHaMS0HxtkjL1zAtXZMfcUbwI325RuBrr45hYXdo34C2GyMeagPxZUTGK0mIonAeVj7Hz4EroxGXMaYu4wx+XaOuhb4wBjz1bDGFOk919H8AXZzcFSPAI9g9ffW027ESATieBXYAKwD3gSGRTsme/s7sPrWa+yfx9rddo8d11bgogjG9CWsD+wWoBR4N9oxtdv+F7BGq+zEaktFdPvt4pgDlACt9v/qZqwe8fvAdvt3VoRjmo3VmljX7vX0hT4Q1xRgtR3XBuCX9vLRWIXDDuBlwB2l5/IsDo7qCVtMOmWDUkrFmAHd6lFKKXU4TfxKKRVjNPErpVSM0cSvlFIxRhO/UkrFGE38SikVYzTxqwFLRBYeOl2ziPxQRB49wv3qO1k2UkSa7DleAtc3HH7vHsX3ORH5Rxe3JYo1bbdX7KnEleotmvjVQDYH60jI9q61lx+NncaYqccWUgcXAu90doMxpsneVl+aBkINEJr41UD2CnCJiLghOEvkUGCRff0OEVku1oln7j3ajYjIaPsEGjNF5CYR+Y+IvCkiu0TkNhH5sX37EhHJanfXc4H3ROR4++Qga+xYxnW1LaV6gyZ+NWAZYyqxDnm/0F50LfCiMcaIyPlY8/vMAqYC00XkjJ5uQ0QmYE3B8XVjzHJ78WTgOvux7wcajTEnAZ8BX7Pvlw20GmNqgVuBP9sV/gysaReUChtN/Gqga9/uad/mOd/+WQ2sAiZifRD0RA7WxFnXG2PWtFv+oTHGY4wpB2qx5mMCa/6lke22P9++/Blwt4jcCYwwxjT1MA6lekQTvxro/gOcK9b5ixONfVYorAnxHjDGTLV/xhpjnujhY9diTWp32iHLW9pd9re77ufgVOgXYff3jTHPY83S2gS8KyLn9DAOpXpEE78a0Ix1mr2FwL/ouFP3XeAb9olCEJFhItLTs0F5sU6H9zURuS7UO9lTFk/BmrESERkNFBpjHsaatnhKD+NQqkcGxIlYlDqCOcBc2o3wMcbMF5HjgM/sk/TUA9fTwxODGGMa7DMoLRCRhhDvNh1YbQ5OjXsNcL2ItAIHgN/0JAalekqnZVYqBPaIoLeMMZN74bH+B+tE7S+EsO5urPMyVBzrdpUK0MSvVAhEpABYDFT28lj+rraXiLXTNwc4wRhTFe5tqtihiV8ppWKM7txVSqkYo4lfKaVijCZ+pZSKMZr4lVIqxvx/SEz7AV8pBbwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Se plotea espectro\n", "plt.plot(v,T)\n", "plt.title('Espectro')\n", "plt.ylabel('Temperatura [K]')\n", "plt.xlabel('Vel [km/s]')" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Valores fiteados: (t0,M,S) = (19.710920444950002, 9.480355666763051, 1.8797736472691602)\n" ] } ], "source": [ "# Se define funcion gaussiana a fitear\n", "def f_gauss(x,T0,mean,stdv):\n", " return T0*np.exp(-((x-mean)**2)/(2*(stdv**2)))\n", "\n", "# First Guess: T0=20, m=10, s=1 obtenidos al inspeccionar espectro ploteado\n", "fg = [20, 10, 1] \n", "coefs,cov = curve_fit(f_gauss,v,T, p0=fg) # Se fitea\n", "\n", "t0,M,S = coefs[0],coefs[1],coefs[2] # Se extraen los coeficientes fiteados\n", "\n", "print ('Valores fiteados: (t0,M,S) =',(t0,M,S))" ] }, { "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.7.4" } }, "nbformat": 4, "nbformat_minor": 2 }