tp-fisica-fibra-optica.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Ejercicio de Optica\n",
    "## • Fisica en la transmision a travez de una fibra optica"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Ejercicio:\n",
    "####     Una fibra tiene un índice de refracción de 1.6 para el núcleo y 1.4 para el revestimiento. Si el rayo incidente toca la frontera con un ángulo de incidencia de 30º, calcule el ángulo derefracción con el que se desvía el rayo en el revestimiento. ###"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Snell's Law"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "\\begin{eqnarray}\n",
    "    µ_1 \\sin\\delta_1 &=&  \\sin\\delta_2 \\\\\n",
    "    \\sin\\delta_1 &=& µ_2/µ_1 \\sin\\delta_2 \\\\\n",
    "    \\delta_1 &=& \\arcsin (µ_2/µ_1 \\sin\\delta_2) \\\\\n",
    "\\end{eqnarray}\n",
    "## El angulo critico se calcula para \n",
    "\\begin{eqnarray}\n",
    "\\delta_2 &=& 90 \n",
    "\\end{eqnarray}\n",
    "## Por lo tanto \n",
    "\\begin{eqnarray}\n",
    "\\delta_1 &=& \\arcsin (µ_2/µ_1 \\sin90)\\\\\n",
    "\\delta_1 &=& \\arcsin (µ_2/µ_1)\\\\\n",
    "\\end{eqnarray}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Datos:\n",
    "\\begin{eqnarray}\n",
    "µ_1 &=& 1.6 \\\\\n",
    "µ_2 &=& 1.4 \\\\\n",
    "α &=& 30° \\\\\n",
    "\\\\\n",
    "\\delta_1 &=& \\arcsin (1.4/4.6) = 17.71° \\\\\n",
    "\\end{eqnarray}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Importo las librerias necesarias"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.lines import Line2D    \n",
    "import math\n",
    "from matplotlib.patches import Arc  \n",
    "from matplotlib.widgets import Slider, Button, RadioButtons"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Defino las constantes que vamos a usar"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "ejex = 0\n",
    "ejey = 0\n",
    "core_radio = 2\n",
    "cladding_radio = 3\n",
    "cable_long = 20\n",
    "core_u = 1.6\n",
    "cladding_u = 1.4\n",
    "angle_ray = 30"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Defino los objetos y las funciones que voy a necesitar"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "class Linea:\n",
    "    def __init__(self, cable, angle, point, color):\n",
    "        self.cable = cable\n",
    "        self.angle = angle\n",
    "        self.x0 = point[0]\n",
    "        self.y0 = point[1]\n",
    "\n",
    "        if point[1] < core_radio and point[1] > -core_radio:\n",
    "            if angle == 0:\n",
    "                self.xf = 50\n",
    "                self.yf = 0\n",
    "            else:\n",
    "                # Si el inicio es en el core\n",
    "                if angle > 0:\n",
    "                    self.yf = core_radio\n",
    "                if angle < 0:\n",
    "                    self.yf = -core_radio\n",
    "                self.xf = core_radio/math.tan(math.radians(angle))\n",
    "\n",
    "\n",
    "        elif point[1] == core_radio:\n",
    "            if angle <= self.cable.get_angulo_critico():\n",
    "                # Rebota\n",
    "                self.yf = -core_radio\n",
    "                self.xf = (point[0] + \n",
    "                           math.fabs((core_radio + core_radio)/\n",
    "                                     math.tan(math.radians(angle))))\n",
    "                self.angle = -angle\n",
    "            else:\n",
    "                # Refracta\n",
    "                self.yf = cladding_radio\n",
    "                self.xf = (point[0] + \n",
    "                           math.fabs(\n",
    "                               core_radio/\n",
    "                               math.tan(math.radians(\n",
    "                                   self.cable.angulo_refraccion(angle)))))\n",
    "                self.angle = self.cable.angulo_refraccion(angle)\n",
    "\n",
    "        elif point[1] == -core_radio:\n",
    "            if angle <= self.cable.get_angulo_critico():\n",
    "                # Rebota\n",
    "                self.yf = core_radio\n",
    "                self.xf = (point[0] + math.fabs((core_radio + core_radio)/\n",
    "                                                math.tan(math.radians(angle))))\n",
    "                self.angle = -angle\n",
    "            else:\n",
    "                # Refracta\n",
    "                self.yf = -cladding_radio\n",
    "                self.xf = (point[0] +\n",
    "                           math.fabs(core_radio /\n",
    "                                     math.tan(math.radians(\n",
    "                                         self.cable.angulo_refraccion(angle)))))\n",
    "                self.angle = self.cable.angulo_refraccion(angle)\n",
    "            \n",
    "        self.color = color\n",
    "        self.draw_line()\n",
    "\n",
    "    def draw_line(self):\n",
    "        ax.add_line(Line2D([self.x0, self.xf],\n",
    "                           [self.y0, self.yf], color=self.color))\n",
    "\n",
    "    def rebotar(self, cable):\n",
    "        if (self.angle > 0 and self.angle < 180 and\n",
    "            self.angle > cable.get_angulo_critico()):\n",
    "            new_line = Linea(cable.angulo_refraccion(self.angle),\n",
    "                             [self.xf, self.yf], 'blue')\n",
    "        else:\n",
    "            new_line = Linea(self.angle, [self.xf, self.yf], 'green')\n",
    "        return new_line\n",
    "\n",
    "class Cable:\n",
    "    def __init__(self, core_radio, cladding_radio, core_u, cladding_u,\n",
    "                 core_color='0.5', cladding_color='#2ca02c'):\n",
    "        self.core_radio = core_radio\n",
    "        self.cladding_radio = cladding_radio\n",
    "        self.core_color = core_color\n",
    "        self.cladding_color = cladding_color\n",
    "        self.core_u = core_u\n",
    "        self.cladding_u = cladding_u\n",
    "\n",
    "    def draw_cable(self):\n",
    "        plt.axhspan(self.core_radio*(-1), self.core_radio, xmin=0, xmax=1,\n",
    "                    facecolor=self.core_color, alpha=0.4)\n",
    "        plt.axhspan(self.core_radio, self.cladding_radio, xmin=0, xmax=1,\n",
    "                    facecolor=self.cladding_color, alpha=0.4)\n",
    "        plt.axhspan(self.core_radio*(-1), self.cladding_radio*(-1), xmin=0,\n",
    "                    xmax=1, facecolor='#2ca02c', alpha=0.4)\n",
    "    \n",
    "        \n",
    "    def angulo_refraccion(self, angle):\n",
    "        return (90-math.degrees(math.asin(math.sin(math.radians(90 - angle)) *\n",
    "                                          (self.core_u/self.cladding_u))))\n",
    "\n",
    "        \n",
    "    def get_angulo_critico(self):\n",
    "        return 90 - math.degrees(math.asin((self.cladding_u/self.core_u)))\n",
    "    \n",
    "    def draw_critic_line(self):\n",
    "        linea_critica = Linea(self, self.get_angulo_critico(), [0,0], 'green')\n",
    "\n",
    "\n",
    "\n",
    "class Grafico:\n",
    "    def __init__(self, x0, y0, xf, yf):\n",
    "        self.x0 = x0\n",
    "        self.y0 = y0\n",
    "        self.xf = xf\n",
    "        self.yf = yf\n",
    "    \n",
    "    def draw_grafico(self):\n",
    "        ax.set_aspect('equal')\n",
    "        ax.set_aspect('auto')\n",
    "        ax.set_xlim(self.x0, self.xf)\n",
    "        ax.set_ylim(self.y0, self.yf)\n",
    "        plt.grid(True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Manos a la obra"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xl8VOXd/vHPNxuBAAmQsIZNSRAFZFVAQSoCsoj7VqsoWuR5fCoqRWr7q9X2qVgVURBRZHNBrAURa0ERJY/UigpKZSvighg22UMChExy//6YZEwgAZIMnMzher9eeZk5Z84914yTizP3nDljzjlERMSforwOICIiJ49KXkTEx1TyIiI+ppIXEfExlbyIiI+p5EVEfEwlLyedmTUzs2wziy68nGFmd3id60SYWU8zW+/B7c40s/89xnpnZq1OZSaJTCp5CRsz22hmBwsLveinsXNuk3OupnMu3+uMx3NkeTrnljrnWldwrHQz+5uZ7TSzfWb2pZndV/SPncipoJKXcLussNCLfracrBuyoCr5HDazM4FPgB+Ads65ROBaoAtQy8tscnqpkn8g4i9m1qJwDzmm2OIzzexTM8sys/lmVrfY9buZ2b/MbK+Z/dvMehdbl2Fmfzazj4ADwBml3F6bwuvtNbM1Zjak2LqZZvacmb1nZvvN7P/MrHnhug8Lr/bvwlch15tZbzPLLLZ9UzN7w8x2mNkuM3umjLv9MPAv59x9zrmtAM659c65nzvn9haO9Tcz21a4l/+hmZ1zxBjJpeUs5f5WM7MnzGyTmW0vvH/Vy8glpxmVvHjlFmAY0AgIABMAzKwJ8A/gf4G6wK+BuWaWUmzbm4HhBPeIvy8+qJnFAn8HFgH1gV8Bs8ys+JTLTcCfgGRgJTALwDnXq3D9uYWvQv56xNjRwNuFt9kCaAK8Vsb9uwSYc5zHYCGQVpjz86Icx8tZikeBdKAD0Kow14PHuW05XTjn9KOfsPwAG4FsYG/hz5uFy1sADogpvJwBPFpsu7OBw0A0MAZ4+Yhx3wWGFtv2j8fI0BPYBkQVWzYbeKjw95nAa8XW1QTygaaFlx3Qqtj63kBm4e/dgR1F9+M4j0UecGk5HrukwttOLE9OwIAc4Mxi1+0OfOf180E/VeOn+MtnkXC4wjm3+ASu90Ox378HYgnusTYHrjWzy4qtjwWWlLHtkRoDPzjnCo4Yv0lp2zvnss1sd9F2x8ncFPjeORc4zvUAdhF8lVKqwlcFfyY4T58CFOVNBvaVI2cKUANYYWah4Qn+gymi6RrxTNNivzcjuOe7k2CBveycSyr2k+Cce7TY9Y916tQtQNMj3pBtBmwu7bbNrCbBaaETeYP4B6DZEe8tlGUxcPUx1v8cuJzgtE4iwVc7ECzo8uTcCRwEzin2eCU652qeQEY5DajkxSu/MLOzzawG8EdgjgseYvkKcJmZ9TezaDOLL3zzM/UEx/2E4Buy95tZbOGbtpdRcu58oJldaGZxBOe8lznnivaOt1PKm7mFPgW2Ao+aWUJhtgvKuO4fgB5m9riZNQQws1Zm9oqZJRF8PyGX4B5/DeCRUsY4Vk4ACl+xvACMN7P6hbfTxMz6l5FLTjMqefHKywTnnbcB8cDdAIUldjnwW4Lz3z8AoznB56pz7jDBUh9AcC/3WeAW59x/il3tVYIlvBvoDPyi2LqHgBcLj8y57oix8wvHbgVsAjKB68vI8Q3BufEWwBoz2wfMBZYD+4GXCE4jbQbWAstKGeZYOYsbA3wNLDOzLIKvIip0bL/4jzmnLw2R04eZzST4Rur/8zqLyKmgPXkRER9TyYuI+Jima0REfEx78iIiPubJh6GSkpJcq1aRe5bUnJwcEhISvI5RYZGcP5Kzg/J7LdLzr1ixYqdzLuX41/yJJyXfoEEDli9f7sVNh0VGRga9e/f2OkaFRXL+SM4Oyu+1SM9vZt8f/1olabpGRMTHVPIiIj6mkhcR8bEqcxbKvLw8MjMzOXTokNdRjisxMZF169Z5HaPCIjl/JGcHqFmzJnl5ecTGxnodRU4TVabkMzMzqVWrFi1atKDYKVOrpP3791OrVuR+g1sk54/k7M45MjMzyczMpGXLll7HkdNElZmuOXToEPXq1avyBS9SUWZGYmJiRLxaFf+oMiUPqODF9/Qcl1OtSpW8iIiEl0q+mJo1j/4yneeee46XXnrJgzQiIpVXZd54rapGjBjhdQQRkQrTnvxxPPTQQzzxxBMA9O7dmzFjxtC7d2/S09NZunQpAPn5+YwePZquXbvSvn17nn/+eQCys7Pp06cPnTp1ol27dsyfP9+z+yEip6cquSf/8N/XsHZLVljHPLtxbf5w2TmVHicQCJCRkcHSpUt5+OGHWbx4MdOmTSMxMZHPPvuM3NxcLrjgAvr160fTpk2ZN28etWvXZufOnXTr1o0hQ4bozTcROWWqZMlXZVdddRUAnTt3ZuPGjQAsWrSIL7/8kjlz5gCwb98+NmzYQGpqKr/97W/58MMPiYqKYvPmzWzfvp2GDRt6FV9ETjNVsuTDscd9slSrVg2A6OhoAoEAEPyQy8SJE+nfv3+J686cOZMdO3awYsUKYmNjadGihY6RFpFTSnPyYdC/f38mT55MXl4eAF999RU5OTns27eP+vXrExsby5IlS/j++3KfJVREpFKq5J68Vw4cOEBqamro8n333XdC291xxx1s3LiRTp064ZwjJSWFN998k5tuuonLLruMdu3a0aVLF84666yTFV1EpFQq+WIKCgqOuT4jIwMInj8lOTk5NCcfFRXFI488wiOPPHLUNh9//HG4Y4qInLCwTdeYWbSZfWFmb4drTBERqZxwzsmPBCL3HLAiIj4UlpI3s1RgEDA1HOOJiEh4hGtO/ingfqDME32b2XBgOEBKSkpofrtIYmIi+/fvD1Ockys/Pz9ispYmkvNHcnYI5j906NBRz/9IkZ2dHbHZIfLzV0SlS97MBgM/OudWmFnvsq7nnJsCTAFo3bq1O/Ib09etWxcxXwYRyV9cAZGdP5KzQzB/fHw8HTt29DpKhWRkZHDk324kifT8FRGO6ZoLgCFmthF4DbjYzF4Jw7giIlJJlS5559wDzrlU51wL4AbgA+fcLyqdzAPR0dF06NCBtm3bctlll7F3797jbjNhwgTatGnDTTfdxMyZM/mf//mfCt22TmksIieDPvFaTPXq1Vm5ciWrV6+mbt26TJo06bjbPPvss7z33nvMmjWrUrc9YsQIbrnllkqNISJypLCWvHMuwzk3OJxjeqV79+5s3rw5dPnxxx8PnUr4z3/+MxAs5m+//ZYBAwYwfvz4Etvv2LGDq6++mq5du9K1a1c++ugjAEaOHMkf//hHAN5991169epFQUFBiVMar1y5km7dutG+fXuuvPJK9uzZA/x0quPzzjuvxKmORUTKUjU/8brwN7BtVXjHbNgOBjx6QlfNz8/n/fff5/bbbweCZ5ncsGEDn376Kc45Bg4cyIcffshzzz3HO++8w5IlS0hOTmbmzJmhMUaOHMm9997LhRdeyKZNm+jfvz/r1q1j7NixdO3alZ49e3L33XezYMECoqJK/lt7yy23MHHiRC666CIefPBBHn74YZ566ikgeKrjTz/9lAULFoROdSwiUpaqWfIeOXjwIB06dGDz5s20adOGvn37AsGSX7RoUeiIiKysLDZs2ECvXr3KHGvx4sWsXbs2dDkrK4vs7Gxq1qzJCy+8QK9evRg/fjxnnnlmie327dvH3r17ueiiiwAYOnQo1157bWh9aac6FhEpS9Us+RPc4w63ojn5AwcO0L9/fyZNmsTdd9+Nc44HHniAO++8Ezixw/gKCgpYtmwZ8fHxR61btWoV9erVY8uWLeXOWNqpjkVEyqI3XktRo0YNJkyYwLhx4wgEAvTv35/p06eTnZ0NwJYtW/jxxx+POUa/fv2YOHFi6PLKlSsB+P777xk3bhxffPEFCxcu5JNPPimxXWJiInXq1AnNt7/88suhvXoRkfKqmnvyVUDHjh1p3749s2fP5uabb2bdunV0794dCO7xz549m/r165e5/YQJE7jrrrto3749gUCAXr16MXnyZG6//XaeeOIJGjduzLRp07j11lv57LPPSmz74osvMmLECA4cOMAZZ5zBjBkzTup9FRH/MufcKb/R1q1bu/Xr15dYtm7dOtq0aXPKs1SEHz51Gan5Izk7BPNnZmZGzHP9SJH+idFIz29mK5xzXcqzjaZrRER8TCUvIuJjKnkRER9TyYuI+JhKXkTEx1TyIiI+ppIvZtu2bdxwww2ceeaZdO7cmYEDB/LVV18ddb1LLrkEgI0bN/Lqq6+Gli9fvpy777670jlGjx7NWWedFTpBWdEpj/Py8hg6dCjt2rWjTZs2jB07ttTtP/jgAzp16kTbtm0ZOnRo6JOxGRkZJCYmcsEFF9ChQ4fQidIAWrRoQbt27ejQoQNduvx0hNbu3bvp27cvaWlp9O3bN3SytFmzZtG+fXvatWtHjx49+Pe//10iQ35+Ph07dmTw4J/OV3f77bdz7rnn0r59e6655prQh8vuvfdeOnToQIcOHUhPTycpKSm0zaWXXkpSUlKJcQDef/99OnXqRIcOHbjwwgv5+uuvAZg5cyYpKSmh8aZO/ekbKceMGUPbtm1p27Ytf/3rXys81pIlS0LLOnToQHx8PG+++SYAa9asoXv37gwdOpSCgoIy/g+LnELOuVP+k56e7o60du3ao5adSgUFBa5bt25u8uTJoWUrV650H374YehyXl6ec865rKws55xzS5YscYMGDQp7lnfffTd0W/fff7+7//77nXPOzZo1y11//fXOOedycnJc8+bN3XfffVdi2/z8fJeamurWr1/vnHPu97//vZs6dWqJvEX5i2vevLnbsWPHUctHjx7txo4d65xzbuzYsaEsH330kdu9e7dzzrkFCxa48847r8R248aNczfeeGOJx2ffvn2h3++9997QuMVNmDDB3XbbbaHLixcvdm+99VZonKLsaWlpoefMpEmT3NChQ51zzs2YMcPdddddR4379ttvu0suucTl5eW57Oxs16VLl1Ce8o5V3K5du1ydOnVcTk6Oc865YcOGuR9//NFNmDDBLVy48KjrZ2Vlef5cr4wlS5Z4HaFSIj0/sNyVs2+1J19oyZIlxMbGMmLEiNCyc889l/z8fHr27MmQIUM4++yzAWjUqBEAv/nNb1i6dCkdOnRg/PjxZGRkhPY4s7Ozue2222jXrh3t27dn7ty5AMyePZt27drRtm1bxowZU2qWfv36ERMT/DByt27dyMzMBMDMyMnJIRAIcPDgQeLi4qhdu3aJbXft2kVcXBzp6ekA9O3bN3TbFTF//nyGDh0KBE+WVrTH2qNHD+rUqXNURoDMzEz+8Y9/cMcdd5QYqyirc46DBw9iZkfd3uzZs7nxxhtDl/v06VPqh5/MjKysLCB4UrfGjRsf836sXbuWXr16ERMTQ0JCAu3bt+edd96p0FjFzZkzhwEDBlCjRg0g+ArGzIiKisJ58EFDkSNVydMa/OXTv/Cf3f8J65hn1T2LMeeVXqoAq1evpnPnzqWu+/zzz1m9ejUtW7YssfzRRx/liSee4O233wYo8QXBf/rTn0hMTGTVquApk/fs2cOWLVsYM2YMK1asoE6dOvTr148333yTK664gjvuuIMRI0aUmCoBmD59Otdffz0A11xzDfPnz6dRo0YcOHCA8ePHU7du3RLXT05OJhAIsHz5crp06cKcOXP44YcfQus//vhjevToQWpqKk888QTnnHMOECy6fv36YWbceeedDB8+HIDt27eH/lFr2LAh27dvP+rxmTZtGgMGDAhdvueee3jsscdK/cLt2267jQULFnD22Wczbty4Euu+//57vvvuOy6++OKjtjvS1KlTGThwINWrV6d27dosW7YstG7u3Ll8+OGHpKenM378eJo2bcq5557Lww8/zKhRozhw4ABLliwJ/aNd3rGKe+2117jvvvtCl0eOHMmgQYNIS0srscMg4hVPSj7f5fO3r/5WYlnLvJbsPrQbgEOBQwQKwnuGxUOBQ6HxS5OTl1PqdbIOZ9GpSycSGyWWWLf70G6yDmeRV5AXWl788rvvvcsLL73w0zbV4YPFH9C9Z3eia0WTFcjiiuuuYNEHi+h1aS8ee+ax0LhFxv1lHAVWwICrB7D70G4++dcn5JPP6m9Xs3fPXgZdMoguPbvQomWLEpmnvDiFX438Fbm5ufzskp9BVHDcFme34Iv1X1ArthaLlixiyOVD+Gx18Lw5f1/8dxo3acyOH3dw9eCraXxGY3pc2AOHK/mYWMmMS/9vKVOmTmHB+wuC93vBu9SuV5uW57Tknx/+s8TjAzBu8jgey3+MMfeNYdor07jplptC66a/PJ3BVwxmX94+yCv5/6BoHCswdh/azV+e+Auz582my3ldmPjkRO4aeRdPT36aC/teyOdXfk61atWYOXUmP7/558x/Zz5denWhd7/enN/9fOol16PzeZ3JLcit0FhFtm3dxpervqTrRV1D97F5m+Ys/L+FAEfdDwArMHLyco56/keKmNyYiM0OkZ+/IqrknvyvOv3qlN/mWWefxd/n/b3UdTUSapziNPDqy6+yaMEi5i2cF5rWmPP6HC7udzGxsbGk1E/h/O7ns3LFyqNKvmu3rvzj/X8AsGTxEr7Z8A3w03QJudD30r6MHjmaXTt3US+5Ho2bBKcoUuqnMGjIID7/7HN6XNiDlPopbNu6jYaNGrJt6zaSU5JDt7Nm1Rru+a97+Ov8v1K3XvAVxScff8I7b7/D4ncWk5uby/6s/dx52508P+P50HbR0dFcde1VTHxyYomSnzdnHo+Nf+y4j83OHTtZs2oNXc4Lvuq58porufby4Dn3i3IA3HzbzTz0u4dCl0eNGcWoMaMAGD50OGemnVnhsQDmz53PoCGDiI2NPW5mEa9oTr5Qr969yM3N5cVpL4aWrVm1hmUfLStzm5o1a5K9P7vUdRddfBHTnpsWurx3z146denEv5b+i107d5Gfn8/c1+fSo2ePo7Z9f9H7THxyIrPmzArN9QKkpqayNCN4CuKcnByWf7qctNZpR22/48cdAOTm5vL0uKe59Ze3ArB92/bQPPGKz1ZQUFBA3Xp1ycnJCU2t5OTksGTxEtqcEzyB1oBBA3jtldcAeO2V1xg4eCAAmZsyGXrDUCZPm0yrtFah237wTw+y+pvVrFy/khdeeoGevXvy/Izncc7x7TffAsE5+Xfefoe09J+yf7X+K/bu2UvXbl3LfLyLJNVJIisri683BI+Cyfggg/TWwfcgtm3dFrrewrcXhpbn5+eze1dwb3vNqjWsWb2Gn13yswqNVWTu63O56rqrjptXxEtVck/eC2bGS399id+N/h0Txk2gWnw1mjVvxsDLBpa5zTntziEqOope5/Xixl/cSLsO7ULrRv1mFPffcz8XdL6A6OhoRv92NJddcRkP/ulBLr/0cpxz9L20b2j8kf81klvvuJWOnTsy5t4x5ObmcvXgqwHocl4Xxk0cx+0jbudXw39Fj049cM7x85t/zjntgnPq119xPU89+xSNGjfimfHP8O7CdykoKGDYL4fRq3fwG6zemvcWM16YQUxUDPE14pn60lTMjB0/7uCW64NfIh4IBLj6+qvp069PMNevRzLsF8OY9eIsUpulMv2V6QA8PvZxdu/ezeh7RgMQHRPNBx99UOZj5ZzjrjvuYv/+/TjnaNuuLY9PeDy0ft7f5nHltVce9WbsoD6D2PDVBnKyc2h7ZlsmTJzAxQMvZvyk8dx6461ERUWRlJTEhOcnADDl2Sm88493iImJoU6dOjzzwjNA8PDTQZcMAqBWrVo8N/250Jvb5R0LYNP3m9icuZkLel5Q5n0WqQo8OdVwq/RWbuzbJY/xbnmoJWekn3HKs1SE5RquWuQeORHJ+SM5OwTzf/P9N3wX/53XUSok5rsYAi0j9xvJIj3/da2v06mGRUTkJyp5EREfq1Jz8s65Uj8gI+IX+oCUlMUKAsTmHyAuP5u4QA6x+TnE5ecQG8ghLj+b2MCBCo1bZUo+13LZv2c/terUUtGLLznn2L9vP7mW63UUOUmiCg7/VNCB7KNLu9Rlwf/GFhw65tj5VrFDdatMyW+L3Qa7oNrOal5HOS4LGC4mcvfIIjl/JGcHyDuUx7ba245/RfGGc0QX5JZawnGB7OCedYllhXvZhb9Hu7xjDp8XFU9eTAKHoxPIi04gp1p99tSoWWLZ4ZgEDkcfvawgKg64rtx3qcqUfEFUAVuqbfE6xgmJ9HfoIzl/JGcHiNkaQ0GSzk55UrkCYvMPllrM8QcPEp2ZdYyyPkCUyy97aIy86BrB8i0s4ay4VA5HFxVzjRIFXWJZdA1c1Kmv3CpT8iIiRczlF85FF5uTzj/wU2mXWFY0Z114/fwDGGW/2is4EFVYwsHiPRxTk5y4+sfemw6tqwEWWcerqORF5KSIKsgruTcdKuFSlhX+t6isT2R++nBMYSFHJ3AoNoms+CYcjqkZWhbamy62zG1O5FDLGDiN3vdTyYtI6ULz0z/tOVc/fIionVkllh19JEiwxI83Px2Iqsbh6Jqhss6plsLe6BbFirlGiYIuviw4P11+MRYDFrnTfRWhkhfxM1dATP6hYBkfcTheyYIuWlZ0BEiwrEudny52BunD0TUKpzeCBbw/vvFPxVxsT7u0ZV7MT5+O9CiLVHWu4Dh7zqUcCVLscL1jzk8TddRc9IG45J/mpIsXc0wCBdsTOdgsPmLnp09HKnmRU8AKAsHpjvxDRGVnHecY6mJ704EcYgsOHnPsfIsJ7U0fjk4gNzaR/fGNj3k4XtGyQFR8ueanY3bFEKh2ek13RDqVvMiJcI7ogsOhqY3yHkMdU3D4p7H2Hj18cH46IXTUR061FPKim5c8HK9oTvqIZQUWe1q9kSjlU+mSN7OmwEtAA8ABU5xzT1d2XJGwc46YgoPHKeaiuesDJT7kEpeffczjpwHyoqsH30gsLOv98Y2OOF46gfzdtTnUKP7ooz40Py0nSTieWQFglHPuczOrBawws/ecc2vDMLZISS6f2EB2GcdLH/8Y6mPNTzusxN50aH66jMPxjtzDPpH56ZicGAKJmu6QU6fSJe+c2wpsLfx9v5mtA5oAKnkpkxUE2LPmIM1Tt1E9KpvovCyy92fSqFr1I44EOeIY6vyDUPZX9VJg0SUOy8uNqR3coy5+ON4RH4QpWhaIqq5pD/GdsL5GNLMWQEfgk3COK/6S+U0tPpqewJwre8E3R68PRFULzT0fjkngYFw99sU0Iy86gcD+WhxKrnHUUR9FR4LkR8WpqEWKCVvJm1lNYC5wj3Muq5T1w4HhACkpKSR8nRCumz7lAoEA1b6u+idSK4tX+bdvr86sWWksWdKEpFr7ebDR/WxIfZmNbiuJcc3p0+Aa2ta+EFfaB13ygz+B6AAx+49+2kYX/lR1eu54K9LzV0RYSt7MYgkW/Czn3BulXcc5NwWYApCWluYaNGgQjpv2xPbt21H+E7dzZwwzZjRi7txkoqLg/MHL2NjxTv5kqzgr7izuTL+fC+tfeEKnmNZj7y3ljzzhOLrGgGnAOufck5WPJH6RlRXNyy83YPbs+uTlRdGhz+ds6fJffBz3Ga1rt+bX6ePo1aCXvj9A5CQKx578BcDNwCozW1m47LfOuQVhGFsi0MGDUbz2Wn1eeqkB+/fH0K7nanZ2v4sVNT4kvXY6o9Kf4KIGF6ncRU6BcBxd809Af63C4cPGvHnJTJ/eiF27Ymnd5Sv29bybVYnvklYrjcfTH+eihhcRpY/Ci5wy+gSGVFp+PixYUI8pUxqxdWs1WrT9ntgb72N98hu0qtWKx9Ifo3fD3ip3EQ+o5KXCnIMPPkjiueca89131WnSaisNfjmGjY1f5szaZ/Jo2qNc3OhilbuIh1TyUm7OwSef1GLSpCasW5dASupO6t8yis0tJ3NGrTMYmz6WPo36qNxFqgCVvJTLl18m8MwzTfj881ok1d9HyvX3syP9SVomNuOR9Efo06gP0RYJR6yLnB5U8nJCNmyozrPPNmbp0iRqJuWQfMUf2dn2z7RIaszItIfp27ivyl2kClLJyzFt2lSN559vzKJFdYivkUu9gePY1eFB6tVJ4Z70B1XuIlWcSl5KtX17LFOnNuKtt5KJiglQ5+Ln2d3lAVKSazMy7Tf0b9Jf5S4SAVTyUsLevdHMnNmQ11+vT36Bo1b3Wew9/9ckNIhjZPq99G/cnxid+1wkYuivVQDIzo5i1qwGvPpqAw4cNGp3ns++HvdSs0k+96T9F5c2uVTlLhKB9Fd7mjt0yJgzJ4UZMxqyb18stdovwl14DzVbZHFP2h0MaDJA5S4SwfTXe5oKBOCtt5KZOrURP/4YR0Lrj+C6kdRO28Y9rYYxKHWQyl3EB/RXfJopKICMjEa89lobfvghnuotVsLQe0g8ewP3pg1jcOpglbuIj+iv+TThHPzzn4k8+2xjNmyoQXyjr+DGUdQ+dzn3pd3O4KaPEhsV63VMEQkzlfxpYPnymkya1IRVq2oSl7wJrh5D7c4Z3Nv6NoY0/Z3KXcTHVPI+tnZtDZ59tjHLliUSm7gdBt9HUvcFXNv4Cm46+w3iokv5mj0R8RWVvA999108kyc35oMP6hCdsAf6jSKp5xyGtbmRy5v+jT0796jgRU4TKnkf2bIljilTGvGPBXWx2INw0UMk9X6JYW2v4Ypms6kWfXp9gbGIqOR9YefOGKZPb8TcN+pRQAB3/pPU6fMCw84dwpXNXlG5i5zGVPIRrOiLsme9mszhPIOOU0nq+wy3d+rPlc2nEx8d73VEEfGYSj4CFX1R9owXkzmQXQ3avkpS/6cYdn5Prmr+nMpdREJU8hHk8GHjjTeSmTItmaw9NSD979Tu/zjDenbmmubjVe4ichSVfAQo+qLsZ55LZtf2mtA8g1rXPcrtF7fhmhZjVe4iUiaVfBVW9EXZT0+qx5ZNSdBoOQm3PcKw/k25rsX/o3pMda8jikgVp5Kvgoq+KPuJCXXZ+FUyJK+jxo13c9uQRK5veS81Ymp4HVFEIoTuY0zMAAAK6ElEQVRKvor5978TePzpJP7zZUNI3Ej1qx/ktqtiueHMX6rcRaTcVPJVxIYN1fnL07VZuSwVErYRP3g0Q687zM/TbyAhJsHreCISoVTyHtu0qRqPP1OTj5c0g2r7qNbvYW6+cQ+/OOtKasbW9DqeiEQ4lbxHtm+PZdyz1flg4RkQfYi4Xk9y001bGdpuiMpdRMJGJX+K7d0bzZPPx/POvDMoKIDY817ghlu+5faOg1TuIhJ2KvlTJDs7ionTY3nz9Zbk58YT02E21w1dw53n96dWbFev44mIT6nkT7JDh4znXoniry+3JC8nkehz3uSaW77gv3v2pnZcG6/jiYjPqeRPkkAApr/ueHF6M3L3phDdajFDhn7MyD4XkBg3xOt4InKaUMmHWUEBzHorwAtTUjnwYxOimi7j0v+ezq8HdyEpbqDX8UTkNKOSDxPnYO57h3jm2UZkZ56BNfiSS0bNYcxVbalT7RKv44nIaUolHwZvL83mqWdS2PtNZ6zuN/T878n8/ufp1I2/0OtoInKaC0vJm9mlwNNANDDVOfdoOMat6t5fvo/HJySyc+1FUHsz3W6bwR9ubUZKwnleRxMRAcJQ8mYWDUwC+gKZwGdm9pZzbm1lx66qln+dw+TfH2Lb5xdDjZ10vvFVHvplAxrVbu91NBGREsKxJ38e8LVz7lsAM3sNuBzwXclv23OQEX/IJfPjqyD2AO2ueIOH/zuRZnVbex1NRKRU4Sj5JsAPxS5nAucfeSUzGw4MB0hJSWH79u1huOlTKy9QwI8bW3LmxW8xamiA5vUSIC834u5LIBCIuMxFIjk7KL/XIj1/RZyyN16dc1OAKQBpaWmuQYMGp+qmwyrjjRz27KpGgwbNvI5SYdu3bydSH/9Izg7K77VIz18RUWEYYzPQtNjl1MJlvhQXE+11BBGRExaOkv8MSDOzlmYWB9wAvBWGcUVEpJIqPV3jnAuY2f8A7xI8hHK6c25NpZOJiEilhWVO3jm3AFgQjrFERCR8wjFdIyIiVZRKXkTEx1TyIiI+ppIXEfExlbyIiI+p5EVEfEwlLyLiYyp5EREfU8mLiPiYSl5ExMdU8iIiPqaSFxHxMZW8iIiPqeRFRHxMJS8i4mMqeRERH1PJi4j4mEpeRMTHVPIiIj6mkhcR8TGVvIiIj6nkRUR8TCUvIuJjKnkRER9TyYuI+JhKXkTEx1TyIiI+ppIXEfExlbyIiI+p5EVEfEwlLyLiYyp5EREfU8mLiPiYSl5ExMcqVfJm9riZ/cfMvjSzeWaWFK5gIiJSeZXdk38PaOucaw98BTxQ+UgiIhIulSp559wi51yg8OIyILXykUREJFzCOSc/DFgYxvFERKSSYo53BTNbDDQsZdXvnHPzC6/zOyAAzDrGOMOB4QApKSls2rmpQoGrAguY8nskkrOD8nst0vNXxHFL3jl3ybHWm9mtwGCgj3POHWOcKcAUgFbprVygZaCsq1Z5Md/FoPzeiOTsoPxei/T8FXHckj8WM7sUuB+4yDl3IDyRREQkXCo7J/8MUAt4z8xWmtlzYcgkIiJhUqk9eedcq3AFERGR8NMnXkVEfEwlLyLiYyp5EREfU8mLiPiYSl5ExMdU8iIiPqaSFxHxMZW8iIiPqeRFRHxMJS8i4mMqeRERH1PJi4j4mEpeRMTHVPIiIj6mkhcR8TGVvIiIj1XqS0MqKtqiuTb9Wi9uOiwytmTQO7231zEqLJLzR3J2UH6vRXr+itCevIiIj6nkRUR8TCUvIuJjKnkRER9TyYuI+JhKXkTEx1TyIiI+ppIXEfExlbyIiI+p5EVEfEwlLyLiYyp5EREfU8mLiPiYSl5ExMdU8iIiPqaSFxHxMZW8iIiPqeRFRHwsLCVvZqPMzJlZcjjGExGR8Kh0yZtZU6AfsKnycUREJJzCsSc/HrgfcGEYS0REwsicq3g3m9nlwMXOuZFmthHo4pzbWcZ1hwPDAVJSUjq//vrrFb5dr2VnZ1OzZk2vY1RYJOeP5Oyg/F6L9Pw/+9nPVjjnupRnm+OWvJktBhqWsup3wG+Bfs65fccr+eJat27t1q9fX56cVUpGRga9e/f2OkaFRXL+SM4Oyu+1SM9vZuUu+ZjjXcE5d0kZN9YOaAn828wAUoHPzew859y28oQQEZGT47glXxbn3CqgftHl8uzJi4jIqaHj5EVEfKzCe/JHcs61CNdYIiISHtqTFxHxMZW8iIiPqeRFRHxMJS8i4mMqeRERH1PJi4j4mEpeRMTHVPIiIj6mkhcR8TGVvIiIj6nkRUR8TCUvIuJjKnkRER9TyYuI+JhKXkTEx1TyIiI+dtwv8j4pN2q2H4jcb/KGZCCSv+YwkvNHcnZQfq9Fev7Wzrla5dkgbN8MVU7ry/uN41WJmS1Xfm9EcnZQfq/5IX95t9F0jYiIj6nkRUR8zKuSn+LR7YaL8nsnkrOD8nvttMvvyRuvIiJyami6RkTEx1TyIiI+5lnJm9njZvYfM/vSzOaZWZJXWU6UmV1qZuvN7Gsz+43XecrDzJqa2RIzW2tma8xspNeZKsLMos3sCzN72+ss5WVmSWY2p/B5v87MunudqTzM7N7C585qM5ttZvFeZzoWM5tuZj+a2epiy+qa2XtmtqHwv3W8zFiWMrJXqDO93JN/D2jrnGsPfAU84GGW4zKzaGASMAA4G7jRzM72NlW5BIBRzrmzgW7AXRGWv8hIYJ3XISroaeAd59xZwLlE0P0wsybA3UAX51xbIBq4wdtUxzUTuPSIZb8B3nfOpQHvF16uimZydPYKdaZnJe+cW+ScCxReXAakepXlBJ0HfO2c+9Y5dxh4Dbjc40wnzDm31Tn3eeHv+wkWTBNvU5WPmaUCg4CpXmcpLzNLBHoB0wCcc4edc3u9TVVuMUB1M4sBagBbPM5zTM65D4HdRyy+HHix8PcXgStOaagTVFr2inZmVZmTHwYs9DrEcTQBfih2OZMIK8kiZtYC6Ah84m2ScnsKuB8o8DpIBbQEdgAzCqebpppZgtehTpRzbjPwBLAJ2Arsc84t8jZVhTRwzm0t/H0b0MDLMJVwwp15UkvezBYXzt8d+XN5sev8juBUwqyTmUWCzKwmMBe4xzmX5XWeE2Vmg4EfnXMrvM5SQTFAJ2Cyc64jkEPVnSo4SuHc9eUE/7FqDCSY2S+8TVU5Lnj8eMQdQ17ezjyp565xzl1yrPVmdiswGOjjqv4B+5uBpsUupxYuixhmFkuw4Gc5597wOk85XQAMMbOBQDxQ28xecc5FStFkApnOuaJXT3OIoJIHLgG+c87tADCzN4AewCuepiq/7WbWyDm31cwaAT96Hag8KtKZXh5dcynBl95DnHMHvMpRDp8BaWbW0sziCL7p9JbHmU6YmRnB+eB1zrknvc5TXs65B5xzqc65FgQf+w8iqOBxzm0DfjCz1oWL+gBrPYxUXpuAbmZWo/C51IcIeuO4mLeAoYW/DwXme5ilXCramZ594tXMvgaqAbsKFy1zzo3wJMwJKtyLfIrgkQXTnXN/9jjSCTOzC4GlwCp+mtP+rXNugXepKsbMegO/ds4N9jpLeZhZB4JvGscB3wK3Oef2eJvqxJnZw8D1BKcKvgDucM7lepuqbGY2G+hN8PTC24E/AG8CrwPNgO+B65xzR74567kysj9ABTpTpzUQEfGxqnJ0jYiInAQqeRERH1PJi4j4mEpeRMTHVPIiIj6mkhcR8TGVvIiIj/1/jG8CUQTRQdMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Inicializo el grafico\n",
    "fig = plt.figure()\n",
    "ax = fig.add_subplot(111)\n",
    "\n",
    "# Dibujo el grafico\n",
    "grafico = Grafico(-2, -5, 12, 5)\n",
    "grafico.draw_grafico()\n",
    "\n",
    "# Dibujo el cable\n",
    "cable = Cable(core_radio, cladding_radio, core_u=core_u, cladding_u=cladding_u)\n",
    "cable.draw_cable()\n",
    "cable.draw_critic_line()\n",
    "\n",
    "linea_nueva = Linea(cable, angle_ray, [0,0], 'blue')\n",
    "\n",
    "linea_nueva_refrac = Linea(cable, linea_nueva.angle ,\n",
    "                           [linea_nueva.xf, linea_nueva.yf], 'orange')\n",
    "if not linea_nueva_refrac.angle > 0:\n",
    "    linea_nueva_refrac2 = Linea(cable, linea_nueva_refrac.angle ,\n",
    "                                [linea_nueva_refrac.xf, linea_nueva_refrac.yf],\n",
    "                                'orange')\n",
    "\n",
    "# Genero el grafico\n",
    "\n",
    "ax.plot(0 , 0, label='Linea')\n",
    "ax.plot(0 , 0, label='Reflexion')\n",
    "ax.plot(0 , 0, label='Critico:' + str(cable.get_angulo_critico()) + '°')\n",
    "\n",
    "plt.title('Fiber optic Cable')\n",
    "ax.legend()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Cuando no supera el angulo critico, se ve como refracta en un angulo menor al que tenia de incidencia.\n",
    "#### Si utilizamos un angulo de incidencia menor (20°), al superarlo, se produce una completa reflexion interna, por lo que el rayo de luz rebota"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Angulo de entrada critico: 28.955024371859857\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEICAYAAABMGMOEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xl4VdXV+PHvykwGkhACIQlzFCdABKnYanEoDq9j+1q1zq1Sf7VVQSuIvs4zKEhVkBartdahWsc6K6itVRFHlAJhTgIhCSEDIfP6/XFPkpuQ4ZLcc4dkfZ4nD/eeYZ+V/Vyy7jl7n3VEVTHGGGO6EhHsAIwxxoQHSxjGGGN8YgnDGGOMTyxhGGOM8YklDGOMMT6xhGGMMcYnljBMSBORYSJSKSKRzvvlInJpsOPyhYgcJSJrgnDcx0Xkjk7Wq4jkBDIm0ztYwjAhQUQ2icgeJzk0/WSq6hZVTVTVhmDH2JW2f4hV9SNVHdPNtvYXkb+LSLGIlInINyIysylxGhMMljBMKDnVSQ5NPwVuHUg8QvLzLyKjgU+BrcBYVU0GzgImAUnBjM30bSH5H8aYJiIywvnmHuW1eLSIfCYi5SLysogM8Nr+CBH5WER2icjXIjLVa91yEblTRP4NVAGj2jnegc52u0TkOxE5zWvd4yKyWETeEZEKEflARIY76z50NvvaOTs6W0Smikie1/5DReQfIlIkIiUi8lAHv/atwMeqOlNVtwGo6hpV/YWq7nLa+ruIbHfOPj4UkYPbtDGwvTjb+X1jRWSeiGwRkULn9+vXQVymj7OEYcLRhcAvgSFAPbAQQESygH8CdwADgGuBF0Qk3WvfC4DpeL6pb/ZuVESigVeBt4FBwO+Ap0TE+7LSecDtwEDgK+ApAFU92lk/3jk7erZN25HAa84xRwBZwDMd/H7HA8930QdvAPs5cX7RFEdXcbbjHmB/4FAgx4nrpi6ObfoqVbUf+wn6D7AJqAR2OT8vOctHAApEOe+XA/d47XcQUAtEArOAJ9u0+xZwkde+t3USw1HAdiDCa9nTwC3O68eBZ7zWJQINwFDnvQI5XuunAnnO6ylAUdPv0UVf1AEn7kPfpTjHTt6XOAEBdgOjvbadAmwM9ufBfkLzx/s035hgO0NV3/Vhu61erzcD0Xi+SQ8HzhKRU73WRwPLOti3rUxgq6o2tmk/q739VbVSRHY27ddFzEOBzapa38V2ACV4zp7a5Zyt3IlnXCMdaIp3IFC2D3GmA/HAShFpbh5P8jVmL3ZJyoSjoV6vh+H5Rl6M54/hk6qa4vWToKr3eG3fWXnmAmBom8HwYUB+e8cWkUQ8l758GZzfCgxrMxbTkXeBn3Wy/hfA6XguXSXjOQsDzx/7fYmzGNgDHOzVX8mqmuhDjKYPsoRhwtH5InKQiMQDtwHPq2fa7V+BU0XkBBGJFJE4Z+A528d2P8UzGH6diEQ7A+an0nqs4WQR+ZGIxOAZI/hEVZu+tRfSzkC64zNgG3CPiCQ4sf2wg21vBo4UkbkikgEgIjki8lcRScEz/lKD50wkHrirnTY6ixMA50zqj8B8ERnkHCdLRE7oIC7Tx1nCMOHoSTzX6bcDccCVAM4fxNOBOXjGC7YCv8fHz7mq1uJJECfh+fb9CHChqv7Xa7O/4fmDvhOYCJzvte4W4AlnhtXP27Td4LSdA2wB8oCzO4hjPZ6xhBHAdyJSBrwAfA5UAH/Bc6ksH/ge+KSdZjqL09ssIBf4RETK8ZzddOveEdP7iao9QMkYX4jI43gGsW8MdizGBIOdYRhjjPGJXxKGiDwmIjtEZJXXsgHOjUPrnH9TO9j3ImebdSJykT/iMcYY439+uSQlIkfjmUP/F1U9xFl2H7BTVe8RkdlAqqrOarPfADzXZSfhmb2yEpioqqU9DsoYY4xf+eUMQ1U/xDO45u104Ann9RPAGe3segLwjqrudJLEO8CJ/ojJGGOMf7l5495gderg4JnNMridbbJofSNRHq1vkmomItPxlHQgLi5u4rBhw/wYavhqbGwkIsKGosD6wpv1RQvrixZr164tVtX0rrdsX0Du9FZVFZEeXftS1SXAEoAxY8bomjUBf8xASFq+fDlTp04NdhghwfqihfVFC+uLFiKyueutOuZm2i0UkSEAzr872tkmn9Z37WbT+q5aY4wxIcLNhPEK0DTr6SLg5Xa2eQuYJiKpziyqac4yY4wxIcZf02qfBv4DjBGRPBH5FZ6yyT8RkXV4at7c42w7SUT+BKCqO/GULVjh/NzmLDPGGBNi/DKGoarndrDquHa2/Ry41Ov9Y8Bj/ojDGBPe6urqyMvLo7q62m9tJicns3r1ar+1Fw7i4uLIzs4mOjrar+1aeXNjTMjIy8sjKSmJESNG4FVyvUcqKipISuo7T7ZVVUpKSsjLy2PkyJF+bdvmmhljQkZ1dTVpaWl+SxZ9kYiQlpbm17O0JpYwjDEhxZJFz7nVh5YwjDHG+MQShjHGeElM3PuBg4sXL+Yvf/lLEKIJLTbobYwxXbj88suDHUJIsDMMY4zpwi233MK8efMAmDp1KrNmzWLy5Mnsv//+fPTRRwA0NDTw+9//nsMPP5xx48bx6KOPAlBZWclxxx3HYYcdxtixY3n55fbuYQ4PdoZhjAlJt776Hd8XlPe4nYaGBiIjIwE4KLM/N596cI/brK+v57PPPuP111/n1ltv5d1332Xp0qUkJyezYsUKampq+OEPf8i0adMYOnQoL774Iv3796e4uJgjjjiC0047LSwH9y1hGGPMPvrpT38KwMSJE9m0aRMAb7/9Nt988w3PP/88AGVlZaxbt47s7GzmzJnDhx9+SEREBPn5+RQWFpKRkRGs8LvNEoYxJiT540wA3LlxLzY2FoDIyEjq6+sBzw1zf/jDHzjhhBNabfv4449TVFTEypUriY6OZsSIEa7cIxEINoZhjDF+cMIJJ7Bo0SLq6uoAWLt2Lbt376asrIxBgwYRHR3NsmXL2Ly5RxXGg8rOMIwxxktVVRXZ2dnN72fOnOnTfpdeeimbNm3isMMOQ1VJT0/npZde4rzzzuPUU09l7NixTJo0iQMOOMCt0F1nCcMYY7w0NjZ2un758uXNrwcOHNg8hhEREcFdd93FXXfdtdc+//nPf/wZYtDYJSljjDE+sYRhjDHGJ5YwjDHG+MQShjHGGJ9YwjDGGOMTVxOGiIwRka+8fspF5Oo220wVkTKvbW5yMyZjjDHd4+q0WlVdAxwKICKRQD7wYjubfqSqp7gZizHG+CIyMpKxY8dSX1/PyJEjefLJJ0lJSel0n4ULF7Jo0SIOO+wwfvKTn/D555/z0EMP7fOxFy9eTHx8PBdeeGF3w3dVIC9JHQesV9Xwvc3RGNPr9evXj6+++opVq1YxYMAAHn744S73eeSRR3jnnXd46qmnenTsyy+/PGSTBQQ2YZwDPN3Buiki8rWIvCEi/ikgY4wxPTRlyhTy8/Ob38+dO7e5fPnNN98MeP7Ib9iwgZNOOon58+e32r+oqIif/exnHH744Rx++OH8+9//BuCqq67itttuA+Ctt97i6KOPprGxsVUZ9a+++oojjjiCcePGceaZZ1JaWgp0XF49EAJyp7eIxACnAde3s/oLYLiqVorIycBLwH7ttDEdmA6Qnp7e6m7LvqyystL6wmF90SJc+yI5OZmKigoAYpfdTMSO73rcZj+FeqeSeOOgg6k55tYu96moqKChoYG33nqLCy64gIqKCt577z2+//573nvvPVSVs88+mzfffJO5c+fyxhtv8Oqrr5KWlsZTTz1FbW0tFRUV/OY3v+HXv/41U6ZMYevWrZx55pl8/vnnzJkzh6lTpzJx4kSuvvpqnn/+eXbv3k1NTQ3R0dFUVFRw/vnnM3fuXH70ox9xxx13cMMNN3DvvffS0NBAVVUV7733Hm+99RY33XQTr7zyyl6/Q3V1td8/A4EqDXIS8IWqFrZdoarlXq9fF5FHRGSgqha32W4JsARgzJgxOnXqVJdDDg/Lly/H+sLD+qJFuPbF6tWrWyrLRsdAZM//RNU31BPV1E50DDFdVK7ds2cPRx11FPn5+Rx44IGcfvrpREZG8q9//Ytly5Zx9NFHA56knJ+fT1JSEiJCYmIiSUlJxMXFERMTQ1JSEh988AHr1q1rbruyshIRYfDgwSxdupSjjz6a+fPnM378eMBTBTc2NpbGxkbKy8s56aSTAJg+fTpnnXUWSUlJREZGcs4555CUlMRRRx3F7Nmz263GGxcXx4QJE3rcf94ClTDOpYPLUSKSARSqqorIZDyXyUoCFJcxJlSddI9fmtmzj+XNm8YwqqqqOOGEE3j44Ye58sorUVWuv/56fv3rX/vcVmNjI5988glxcXF7rfv2229JS0ujoKDA5/aatFdePRBcH8MQkQTgJ8A/vJZdLiJND8n9X2CViHwNLATOUVV1Oy5jjOlMfHw8Cxcu5P7776e+vp4TTjiBxx57jMrKSgDy8/PZsWNHp21MmzaNP/zhD83vv/rqKwA2b97M/fffz5dffskbb7zBp59+2mq/5ORkUlNTm8cnnnzySX784x/789frFtfPMFR1N5DWZtlir9cPAfs+/8wYY1w2YcIExo0bx9NPP80FF1zA6tWrmTJlCgCJiYn89a9/ZdCgQR3uv3DhQq644grGjRtHfX09Rx99NIsWLeJXv/oV8+bNIzMzk6VLl3LxxRezYsWKVvs+8cQTXH755VRVVTFq1Cj+/Oc/u/q7+kLC8cv8mDFjdM2aNcEOIySE67VqN1hftAjXvli9ejUHHnigX9t044l74aC9vhSRlao6qbttWmkQY4wxPrGEYYwxxieWMIwxxvjEEoYxxhifWMIwxhjjE0sYxhhjfGIJwxhjvGzfvp1zzjmH0aNHM3HiRE4++WTWrl2713ZHHnkkAJs2beJvf/tb8/LPP/+cK6+8MmDxBpIlDGOMcagqZ555JlOnTmX9+vWsXLmSu+++m8LCljJ4TaU4Pv74Y2DvhDFp0iQWLlwY2MADxBKGMcY4li1bRnR0NJdffnnzsvHjx9PQ0MBRRx3FaaedxkEHHQR47vQGmD17Nh999BGHHnoo8+fPZ/ny5Zxyiud5cJWVlVxyySWMHTuWcePG8cILLwDw9NNPM3bsWA455BBmzZoV4N+y+wJVfNAYY/bJvZ/dy393/rfH7TQ0NBAZGQnAAQMOYNbkjv9Ar1q1iokTJ7a77osvvmDVqlWMHDmy1fJ77rmHefPm8dprrwG0Kil+++23k5yczLfffgtAaWkpBQUFzJo1i5UrV5Kamsq0adN46aWXOOOMM3ryawaEnWEYY4wPJk+evFey6Mq7777LFVdc0fw+NTWVFStWMHXqVNLT04mKiuK8887jww8/9He4rrAzDGNMSOrsTGBf7EstqYMPPpjnn3++3XUJCQl+iSec2RmGMcY4jj32WGpqaliyZEnzsm+++abTx6AmJSU1PyWwrZ/85CetngleWlrK5MmT+eCDDyguLqahoYGnn346JEqX+8IShjHGOESEF198kXfffZfRo0dz8MEHc/3115ORkdHhPuPGjSMyMpLx48fv9UzvG2+8kdLSUg455BDGjx/PsmXLGDJkCPfccw/HHHMM48ePZ+LEiZx++ulu/2p+YZekjDHGS2ZmJs8999xeyy+77LJW75sepBQdHc3777/fal1TafnExESeeOKJvdo699xzOffcc/0UceDYGYYxxhifWMIwxhjjE0sYxpiQEo5PAQ01bvWh6wlDRDaJyLci8pWIfN7OehGRhSKSKyLfiMhhbsdkjAlNcXFxlJSUWNLoAVWlpKSEuLg4v7cdqEHvY1S1uIN1JwH7OT8/ABY5/xpj+pjs7Gzy8vIoKiryW5vV1dWu/PEMZXFxcWRnZ/u93VCYJXU68Bf1fKX4RERSRGSIqm4LdmDGmMCKjo7e57upu7J8+XImTJjg1zb7qkAkDAXeFhEFHlXVJW3WZwFbvd7nOctaJQwRmQ5MB0hPT29Vr6Uvq6ystL5wWF+0sL5oYX3hP4FIGD9S1XwRGQS8IyL/VdV9LpziJJolAGPGjNGmec593fLly7G+8LC+aGF90cL6wn9cH/RW1Xzn3x3Ai8DkNpvkA0O93mc7y4wxxoQQVxOGiCSISFLTa2AasKrNZq8AFzqzpY4Aymz8whhjQo/bl6QGAy+KSNOx/qaqb4rI5QCquhh4HTgZyAWqgEtcjskYY0w3uJowVHUDML6d5Yu9XitwRdttjDHGhBa709sYY4xPLGEYY4zxiSUMY4wxPrGEYYwxxieWMIwxxvjEEoYxxhifWMIwxhjjE0sYxhhjfGIJwxhjjE8sYRhjjPGJJQxjjDE+sYRhjDHGJ5YwjDHG+MQShjHGGJ9YwjDGGOMT8TyOIrzk7J+jd792d7DDCAlRG6OoH1kf7DBCgvVFC+uLFtYXLX4+5ucrVXVSd/e3MwxjjDE+sYRhjDHGJ64lDBEZKiLLROR7EflORK5qZ5upIlImIl85Pze5FY8xxpiecfOZ3vXANar6hYgkAStF5B1V/b7Ndh+p6ikuxmGMMcYPXDvDUNVtqvqF87oCWA1kuXU8Y4wx7nLzDKOZiIwAJgCftrN6ioh8DRQA16rqdx20MR2YDpCenk7UxoCEHvKkRqwvHNYXLawvWlhf+I/rvSgiicALwNWqWt5m9RfAcFWtFJGTgZeA/dprR1WXAEvAM63Wpsl52JTBFtYXLawvWlhf+I+rs6REJBpPsnhKVf/Rdr2qlqtqpfP6dSBaRAa6GZMxxpjucXOWlABLgdWq+kAH22Q42yEik514StyKyRhjTPe5eYbxQ+AC4FivabMni8jlInK5s83/AqucMYyFwDnq51vP38l9h/zy/L2Wl39ZQs3y1f48VO/T2MiR6+6jf9WWYEcS0lJ3r2dK7v0QhlUTAmnCpj+StfOTYIcR0mJqyzhq7e1ENNQEO5R2uTaGoar/AqSLbR4CHtrXtkWFURWjutxuQ8UGlq5YylkjzuL3h/y+eXlBQQwJX97A8Qe9ybr8h6nuf9i+hhAyChsKGVYxzJW2M4seZkjZ5wzcvZXvRr6ARsS5chx/cbMvOhLRWMUhGy4numEnUzZ/RGHaxQE9fkeC0RedSa5YzujidxhZ8hGrZBr1UYG78hxqfdEhbWTMll+SWL2KY9c9xoasucGOaC+9+k7vB1c/SHxUPJftf1nzssrKCGbMyGH6Y49RHZ3NmMLfE1u7NYhRhqbkyg/J2Pk4lXGHEFtfyIjtt9k36La0gZHbbiSqoZTdsQeSVbyIpN3tTQTs2+Jqchm5/Waqo4cBQk7+TKSxOthhhZzM4kUkVq+ivN9hpFYuI7302WCHtJdemzA+LfqUf+/4N7/c75ekxKQAUF8Pc+aMYtOmOObcspPNI+aDQE7eVUQ2tJ3A1Xf1q17LyIIbqIo9gLVDF5M38LcMqHibISVLgh1aSMkqeoiUyg/ZOuga1g5dRHXMCEYVzCK2ZlOwQwsZUfU7ycmfSYP0Y93QR9iYeSfx1asZue0m0MZghxcyBpS9xpCdf6Yo+UzWZT/CrsSjGbrjfvrv/jjYobXSKxNGgzbw4OoHGdJvCGePOLt5+YIF2Xz8cTKzZm3hBz+ooDYmm/WZ84ipK2BUwSxQm3oXVV9MTv4MGiKTWJ/9ABoRR+GAiyjufyqZJX8ktfzNYIcYEtJ2vURG6ZPsSDmLotRzaIxMJDd7Piox5ORfTWT9rmCHGHTSWMPo/GuJri9hffYD1EVnUJZ4NHnpV5Fa+T6ZxYuCHWJISKz6kuHb76A8/nC2DJ4FEVFsHHIHe2JzGFVwPXE1ucEOsVmvTBiv573O2vK1/PaA3xIbGQvA3/+ezjPPDOYXvyjkpz8tbt62Mn4CmzNupH/VCoYV3tunL7tIYzU5+TOJbCgjN+sB6qLSnRXClow5VPQ7jBHbbyNhz9fBDTTIknavYHjh3ZTFT2HroGual9dGZ7I+ax4x9TsYXfB7pLE2iFEGmSojtt9GYvU3bBxyK1VxBzWv2pF6HkXJZzJk558ZUPZaEIMMvtjarYwuuJbamCw2ZN4L4hlWboyIJzdrPg3Sj5z8mUTV7wxypB69LmHsqd/DI/99hINTDmZa5jQA/vOf/sybN5SjjtrFVVfl7bXPzuRT2DbgEtLLXmRQ6VOBDjk0aCMjt91MfPVqNmbeyZ64A1qvlmjWZ91HbdRgRudfS0xdQZACDa7Y2k2MKriO6pjhbMi8u/k/eJPd/caxKeNmkvZ8ybDCu/rsF5AhJX9iQMVb5A+8gl1Jx7deKcKWwbMoj5/M8O13kFj1RXCCDLLIhnJy8q8GhdysBTRE9m+1vi56MOuzHyC6voTR+dcijcGfOdXrEsZTG56iqKaIGQfNQERYvz6O2bNHMXr0Hu68cyORke3vVzDw/1GaeBzZRQ+SXPlBYIMOAZnFi0mtfI+89KsoS/xxu9s0RKaQmz0f0Xpy8q4moqEywFEGV2TDLnLyZqASRW7WfBojE9vdrrT/CRSk/ZqB5a+RsfPxwAYZAlLL3yKz5FGK+5/C9gEXt7+RRLEh8x5qY7IYnd8HJ55oPaMKZhNTm8/6rLnUxAxtd7OquIPYOMRzphYKE096VcIori7mifVPcEzGMRw64FB27oxixowc4uIaeeCBXOLjOxlkkwjn1PlARhbcSL/qNYELPMg8A26PUZR8BjtSz+t025qYEWzIvI+42s2MKri+z4z7iNYxOv86YuoLWZ81j9qYzutobku7lJKkE8kqfpiUincDFGXwJez5hhHbb6Wi3wS2DJ4D0vHM+obI/uRmLQDwjPv0lYknqgwrvJf+VZ+xOeNGKuM7n9a/K+k48gdewYCKtxhS8scABdm+XpUwHl37KHWNdVx54JXU1AjXXjuakpJoHnggl4yMui7314g4crMeoCGyPzn5M4iuLwpA1MHVesBtdqf/wZtUJBzO5sHXk1z1H4buaPcm/t5FlWHb7yJpzxdsyriJ3f3Gd72PCJsz/o/KuHGeS3172q2p2avE1BUwOv9aaqMGsT5zLhoR0+U+NTFDWZ81l5jafEYVzO4TX0AGlf6N9LIX2TbgEnYm+/Zkh+0DLqa4/ylkliwJ6sSTXpMwcstzeXnLy5w14iyy44dy++3D+eabRG69dSMHH1zlczv1UQPJzZpPZEMFo/Ov6dXzxWNq8zwDbtGZrQbcfFGScgbbU89n0K7nQnK+uD8N3vkEA8tfpSDtMkr7n+jzfhoRy/qsedRFpZGTP5Pouu0uRhlcEQ2V5OTNQLTOcz0+KsXnfSvjD3MmnnzW6yeeJFd+SHbRAkoTj6Vg4P/zfUcRtgyeQ0W/Cc7Ek2/cC7ITvSZhLFy9kIToBC7d71KWLs3gzTfT+M1v8jn++H2f3rgnbn+v+eL/1yvni0c2VLQMuGXvPeDmi/z034XsfHF/Sal4n+zih9iZNI1tadP3ef/6qAHkZs0nQj0z0CIaff/yEja0nlHb5hBXu4kNmfdSEztin5voCxNPvO9v2jjkNpB9+/OrETGsz5xLbdSgoE086RUJ45OiT/i46GN+td+v+HT5CBYvzuJ//qeESy7p/jc6z3zxq0mtXNb75otrvecGs9o81mfdR01MN8smSGTIzhf3h6YvDJVxY9mUcbNPl+vaUx07mg2Zd9OvJpeRBTeANvg50uAaumM+ybs/ZsvgWVQkTO52O56JJ8c6E08+9GOEwddyf1P/5vubuqMhKoXcrAWI1pGTNyPgE0/CPmE0aAMLvl9AVnwWB5X/kltuGcGhh1Zwww2bu/v/u9mO1F80zxdPK3vVPwEHmyrDCu+jf9VnbMm4gcr4iT1qrtV88bwZITNfvKei6woZnTeDusgBrM+ah0bE9qi98oQj2TroWlJ2f0RW0R/8FGXwpZc+x6Bdz1KYeh7FKT/tWWMSwcYhtzkTT27oNRNPWu5vKm99f1M31cSOYEPmvcTVbmLUtjkBHfcJ+4Tx6tZXya3I5fwBNzD7ugNIT69j3rz1xMT44Tqo13zxYdvv7BXzxQeVPk162T/YNuBiSpJP9UubzfPFG3aGzHzxnohorCInfwaRuofc7PnUR6X5pd2i1J+zI+VsMkr/ysBdL/qlzWDqv/s/DN1xP7sSjiIv/Uq/tNl64slMouqLu94plLW6v+kO9sSN8UuzFQmT2TJ4Nsm7P2bojvl+adMXYZ0wquqrWLxmMQf3m8Lzd19Eba3w4IPrSEnx4yl/L5ovnlz5EdlF8ylNPIaCgb/xa9ue+eK3h8x88W7TBkZu+z/61eSyYchdVMfm+LX5rYNmUJZwJMMK7yFp92d+bTuQ4mrWM6pgNntiR7Mx806QDm5w6gbPxJMHnBvbwrtQYWbxoi7vb+qu4pQzKUw9j0G7niW99Dm/tt2RsE4YT65/kuKqUiKef45Nm/px770bGDHC/99uGyL7sy7rwbAuVOgZcJtDVewBbOrGgJsvdiUd6xQqDP588e7yFBT8gK2DZlKe+EP/H0CiPIkoZjijCq4Ly0KFnoKCM2gUz9lAY0S834+xJ25M2Bcq9C4o2NX9Td2Vl34luxKOciae/MeVY3gL24RRVF3EkxueZOi//863K7K57jpPQUG3hHOhwrYFBRsj+rl2LE+hwuDPF++OtgUF3eIpVLggLAsVehcUzM3yFBR0SzgXKtyroGBPB1Q7IpGeUj6xoxlVMJu4mvXuHMcRtglj0ZpF1H5yKVvfP5Nf/KKQn/3M/Wud4ViosMOCgq4dUNiScYNXocLgzBffVx0VFHRLWBYqVGX49ttbCgr2O9j1Q4ZjocKe3N/UHZ6JJw84hQrdnXgSlgmjprGGV5bVoG8s6LCgoFvCar64NjJi+y2eAbchd+xVUNC1w0o06zPvc+aLXxPyhQpjazd3WlDQLZ5ChTeFTaHCISV/Iq3iTfIH/mbvgoJuaZ54crhTqPDLwBy3m1rf3/Rgt+5v6o666IyAFCp0PWGIyIkiskZEckVkdjvrY0XkWWf9pyIyoqs2i3ZXwt+fZeSoPdxxR8cFBd0SLoUKM0seZUDFu+SnX0lZ0tSAHrshKoXc7AUhX6jQU1Dw6i4LCrqltP+JFKRNZ2D5awze+URAj72Cc0pwAAAUyklEQVQvWhcUvCSwB5coNmTe65l4UnBt6E48aXV/U8cFBd3imXhyK4nV3zB8++2ufAFxNWGISCTwMHAScBBwrogc1GazXwGlqpoDzAfu7ard6uIs4uOEhQs2kJAQhMGwMChUOKDsnwwpWUpR8hkUpp4flBg8hQrv9RQqDPB8cV/sa0FBt2xLu4ydSSeQXfwQKRXvBSWGziTs+dbngoJuaS5UqCFaqHCv+5s6Lyjoll1Jx5M/8ArSKt5kSMmf/N6+22cYk4FcVd2gqrXAM8DpbbY5HWj6avU8cJxIF5/IhmgWLtjsU0FBt4TyfPGEqq8YXngH5fGT2OrmgJsPWs8XD6FChd0pKOgWETZl3OQUKryJ+OrvgxdLG56CgtfsU0FBt4RyoUJPQUH/3t/UXS2FCh8ltfwtv7bt9sXaLMD7/DEP+EFH26hqvYiUAWlAq7/AIjIdmA7Qr38mQwZtpbDQrbB9V97vRqZU/J7hm37HJ4n30Cjdu+W/u+rr6yn06oj4hm2MrZxJlQzik6hrqd8R/DuvCzmShtgzGbXrOYpqBrA51p3/UG37ojOjqv/OwOpXWRd7Luv2TIA9wf8w7YyZxZG1Mxm55So+TppPdUT3JyjsS190JEqrmFJxLTTW8Gn8XewuqQGC20+FZFEd/zvGV80nfdMtrOr32y6/EPmjLzozqO5TsncvYHv0kXzZ8DNC4Q9TkVzK5MiNDN92C9sr4tgV5Z/xy8CM7vmBqi4BlgDst99+Onjw4CBH1GQwmyrvYnT+Nfyg4WE2ZN7jyj0OHSksLKSpLyIbKhiz5Q4iRNg0/A+kdbdGlAtKdTa7Cko4qHIJMWkHUZ5wpN+P4d0XnUmpeJ/Rux5nZ9I0yofMZHAQz8BaG8zGmgc5YMuvOKLmLtYM+1O373HwtS86pPXk5N9JYmMe67IXkpgwicCO7nSsnvPYVrSLYTv/jCQfyI4Bnd/j0OO+6ES/6rWM2TKXqtgDKBh2H4O7WSPKDVvqH+SALRdz+J47+e/wx6mNzuxxm27/ZcsHvEd+sp1l7W4jIlFAMlDiclx+FRLzxZ0nePW4oKBbQqRQYeuCgjcF9XJde6pjc1oKFW67MWiFCrN3LPAqKNj2okDwtRQqXBC0QoVt72/qbkFBt7hRqNDthLEC2E9ERopIDHAO8EqbbV4BLnJe/y/wvmqIzy9sR1Dni6syrHAu/as+ZUvGnB4XFHRL6/nigX+wfXTdDkbnz/QqKBha/8GbeAoVXkNK5YdBKVSYXvp3Bu96hsLUX/S8oKBb9ipUuDawhw/0/U3dtFehwh5yNWGoaj3wW+AtYDXwnKp+JyK3ichpzmZLgTQRyQVmAntNvQ0LQXyw/aBdz5Be9gLbB1xESfJpXe8QRIGaL95Wc0HBxiq/FhR0S1Hq2exI+XnACxV6CgrOcwoKXhWw43ZH2ydkBmziSauCgncG7P6m7vJMPJlFsh+eWeP6xXZVfV1V91fV0ap6p7PsJlV9xXldrapnqWqOqk5W1Q1ux+SaIBQqTK/7jOwdnoKC+QOvcP14/hDwB9s3FxRc50pBQbdsHTQzoIUK3Swo6JZgFCrMLF7sWkFBtxSn/JRCP9SzCss7vUNZq0KFLs8X71e9jgm776Mqdn/XCgq6JZAPtne9oKBbWhUqnOVqocJAFBR0i6dQ4R3O+NTNrhYq9BQUfIyi5DNcKyjoFn+UoA+fvzBhxFOosGm+uDuFCj0DbldTJ/Gsz57vakFBtwTiwfatCgqmnO3KMdzUUqgwyvkC4v9ChYEsKOiWssQfOxNP3nNt4knrgoKzQ27CRJf8cMZoCcMlLQ+293+hQs+A2zVENpSxMuGmkB1w65LLD7ZPrPq8dUHBcPsP7vAUKrzfU6gw/zr/FioMQkFBt7g58STQBQVDlSUMF7lSqLC5oOD3bBxyB+VR4XE9viNuPdg+tnYzo/MDX1DQLS2FCr/wa6HCjJKlgS8o6BaXChW2Lii4IGAFBUORJQyX+btQYUtBwd8FvKCgW/w9X7yloGBkUAoKusVTqPAyvxUqTC1/m6ySxZT0/5/AFxR0i78LFbYqKBiC9zcFmCUMt/mxUOGAstcZUrKU4uTTKUy9wI9BBp+/HmzfUlBwe1ALCrplW9p0r0KF73e7HU9BwVuo6DeBzYNvCNvLde1pW6gwqrGbD1bbq6BgaN7fFEiWMALAH4UKPQUFb6c8flJ4Drj5oMcPtvcqKLg52AUF3dKqUOH/datQYUzdNkbnX0NdVHrQCwq6xbtQ4WFVd3frC8ig0qdDpqBgqLCEESCe+eLzuzVfvHnALWoIGzLvRSXaxUiDqycPth9V8wIDy1+lIO0ydvY/yaUIg08jYlmfNZe6yAGMzptJdJ3vxe4iGirJybsa0VrWZT1IQ1SKi5EGV9PEk4H1X+/zxJPkyo/ILvLc31Qw8DcuRhleLGEE0J64/ff5wfaeAbcZiDY6A27JAYg0uFoebD+P/j7enZpSsYwx1Z6CgtvSprscYfDVR6WRmz2fSN1DTv4MIhqrut5J6xm1bQ5xtZvYkHkvNbEjXI8z2HYmn0Ju7M+diSd/82mfftVrGVkwh6rYA8Lu/ia3WU8E2D4VKmwuKLi1bw24NT/YvqlQYecPtvck4BvZFTkmJAsKumVfCxW2FBS8LiQLCrplbdwFPhcqbFtQMBzvb3KTJYwg8Gm++F4FBScFNsggaypU2ChxnT7YvqWgYCorE24M2YKCbmldqPChDrdrXVDwZwGMMAQ0FSqMPaDTQoXhUlAwmCxhBIMPhQpbCgpeGPIFBd1SF51BblbHhQpbFxRcQG1EapAiDS5PocKzyCh9krRdL+21PpwKCrpFI+JYn91JocLm+5tWe8rwh3hBwWCxhBEsnRQq7F/5L6+Cgr8NYpDBV9Xv4PYfbK+NYVlQ0C1bB11DWcKRDC+8m6TdK5qXx9VscAoKjmJj5h1hUVDQLXVR6R0WKmy5v+nKXnN/kxssYQRRe4UK+1WvY1TBnLAsKOiW9h5sn1X0h/AsKOiWVoUKryOhIY+o+lKvgoLzaYxICHaUQddeocIBZf9kSMlSipLPoDD1/GCHGNLCu15CL9BUqHC/rb9hdP5MYuu20RCZwPosG3Dztn3AxcTWbiaz5FGi6wtJL3spbAsKuqUxMpHcrPkcsOViJlXejNalEV1fzJqhj4ZlQUG3NBUqHFq0gOHbbmFA5TuUx09i6+BZfWbCRHfZ19cQUBl/GJsHzyZpz1dE1xeTmzWfuuhBwQ4rtDiFCqtichhY9hKVcYeEdUFBt9TGZLE+ax79tJDE6u/YmHETVf0OCXZYIWdH6nmUJE1jYMXr1Ef0Z0Pmfb36/iZ/CcszDBVlQ1L4PmepPRsSxxJfdQglCfvxXXoM4NvvF1Uc1ev6ojM7ci7nR7n38UHOdKritrRa19f6okNJSUQVnkS/uO18kZkD0rf7pKPPxZb4czl+9bd8OfRiClOKgQA9sS+MhWXC6JVE+HDMTcGOIuRVxA/jjXEdTx81Ht8mXkr9SP8/h6U3qY9K4M2xDwY7jLDiSsIQkbnAqUAtsB64RFX3evKLiGwCKoAGoF5V+9bNBsYYE0bcGsN4BzhEVccBa4HrO9n2GFU91JKFMcaENlcShqq+rdpcHvITINuN4xhjjAmcQIxh/BJ4toN1CrwtIgo8qqpLOmpERKYD0wHS09OJ2mjDLwBSI9YXDuuLFtYXLawv/KfbvSgi7wLtTe6+QVVfdra5AagHOno+6Y9UNV9EBgHviMh/VbXd6mBOMlkCkLN/jtqAnkfUxigb3HRYX7SwvmhhfeE/3U4YqtrpA4BF5GLgFOA41fYL0atqvvPvDhF5EZgMdF5O0hhjTFC4MoYhIicC1wGnqWq7hfpFJEFEkppeA9OAVW7EY4wxpufcmiX1EJCE5zLTVyKyGEBEMkXkdWebwcC/RORr4DPgn6r6pkvxGGOM6SFXRoJUtd3SoapaAJzsvN4A9MKHLhtjTO8UllMHIiWSs/Y/K9hhhITlBcuZuv/UYIcREqwvWlhftLC+8B8rPmiMMcYnljCMMcb4xBKGMcYYn1jCMMYY4xNLGMYYY3xiCcMYY4xPLGEYY4zxiSUMY4wxPrGEYYwxxieWMIwxxvjEEoYxxhifWMIwxhjjE0sYxhhjfGIJwxhjjE8sYRhjjPGJJQxjjDE+sYRhjDHGJ64lDBG5RUTynWd6fyUiJ3ew3YkiskZEckVktlvxGGOM6Rm3H9E6X1XndbRSRCKBh4GfAHnAChF5RVW/dzkuY4wx+yjYl6QmA7mqukFVa4FngNODHJMxxph2uH2G8VsRuRD4HLhGVUvbrM8Ctnq9zwN+0F5DIjIdmA6Qnp7O8uXL/R9tGKqsrLS+cFhftLC+aGF94T89Shgi8i6Q0c6qG4BFwO2AOv/eD/yyu8dS1SXAEoAxY8bo1KlTu9tUr7J8+XKsLzysL1pYX7SwvvCfHiUMVT3el+1E5I/Aa+2sygeGer3PdpYZY4wJMW7Okhri9fZMYFU7m60A9hORkSISA5wDvOJWTMYYY7rPzTGM+0TkUDyXpDYBvwYQkUzgT6p6sqrWi8hvgbeASOAxVf3OxZiMMcZ0k2sJQ1Uv6GB5AXCy1/vXgdfdisMYY4x/BHtarTHGmDBhCcMYY4xPLGEYY4zxiSUMY4wxPrGEYYwxxieWMIwxxvjEEoYxxhifWMIwxhjjE0sYxhhjfGIJwxhjjE8sYRhjjPGJJQxjjDE+sYRhjDHGJ5YwjDHG+MQShjHGGJ9YwjDGGOMTSxjGGGN8YgnDGGOMT1x5RKuIPAuMcd6mALtU9dB2ttsEVAANQL2qTnIjHmOMMT3nSsJQ1bObXovI/UBZJ5sfo6rFbsRhjDHGf1xJGE1ERICfA8e6eRxjjDHuc3sM4yigUFXXdbBegbdFZKWITHc5FmOMMT0gqtq9HUXeBTLaWXWDqr7sbLMIyFXV+ztoI0tV80VkEPAO8DtV/bCDbacD0wHS09MnPvfcc92Ku7eprKwkMTEx2GGEBOuLFtYXLawvWhxzzDErezJW3O2E0WXDIlFAPjBRVfN82P4WoFJV53W17ZgxY3TNmjU9D7IXWL58OVOnTg12GCHB+qKF9UUL64sWItKjhOHmJanjgf92lCxEJEFEkppeA9OAVS7GY4wxpgfcTBjnAE97LxCRTBF53Xk7GPiXiHwNfAb8U1XfdDEeY4wxPeDaLClVvbidZQXAyc7rDcB4t45vjDHGv+xOb2OMMT6xhGGMMcYnljCMMcb4xBKGMcYYn1jCMMYY4xNLGMYYY3xiCcMYY4xPLGEYY4zxiSUMY4wxPrGEYYwxxieWMIwxxvjEEoYxxhifWMIwxhjjE0sYxhhjfGIJwxhjjE8sYRhjjPGJJQxjjDE+sYRhjDHGJ5YwjDHG+KRHCUNEzhKR70SkUUQmtVl3vYjkisgaETmhg/1HisinznbPikhMT+Ixxhjjnp6eYawCfgp86L1QRA4CzgEOBk4EHhGRyHb2vxeYr6o5QCnwqx7GY4wxxiU9ShiqulpV17Sz6nTgGVWtUdWNQC4w2XsDERHgWOB5Z9ETwBk9iccYY4x7olxqNwv4xOt9nrPMWxqwS1XrO9mmmYhMB6Y7b2tEZJWfYg13A4HiYAcRIqwvWlhftLC+aDGmJzt3mTBE5F0go51VN6jqyz05+L5Q1SXAEiemz1V1Uhe79AnWFy2sL1pYX7SwvmghIp/3ZP8uE4aqHt+NdvOBoV7vs51l3kqAFBGJcs4y2tvGGGNMiHBrWu0rwDkiEisiI4H9gM+8N1BVBZYB/+ssuggI2BmLMcaYfdPTabVnikgeMAX4p4i8BaCq3wHPAd8DbwJXqGqDs8/rIpLpNDELmCkiuXjGNJb6eOglPYm7l7G+aGF90cL6ooX1RYse9YV4vugbY4wxnbM7vY0xxvjEEoYxxhifhFXCEJETnVIjuSIyO9jxBJKIDBWRZSLyvVOO5Spn+QAReUdE1jn/pgY71kARkUgR+VJEXnPe98lSMyKSIiLPi8h/RWS1iEzpq58LEZnh/P9YJSJPi0hcX/lciMhjIrLD+x61jj4H4rHQ6ZNvROQwX44RNgnDKS3yMHAScBBwrlOCpK+oB65R1YOAI4ArnN9/NvCequ4HvOe87yuuAlZ7ve+rpWYeBN5U1QOA8Xj6pM99LkQkC7gSmKSqhwCReEoU9ZXPxeN4SjF56+hzcBKe2av74bkhepEvBwibhIGntEiuqm5Q1VrgGTwlSPoEVd2mql84ryvw/FHIwtMHTzib9ZnyKiKSDfwP8CfnfZ8sNSMiycDRODMMVbVWVXfRRz8XeO4t6yciUUA8sI0+8rlQ1Q+BnW0Wd/Q5OB34i3p8gueeuCFdHSOcEkYWsNXrfaelRHozERkBTAA+BQar6jZn1XZgcJDCCrQFwHVAo/N+n0rN9CIjgSLgz87luT+JSAJ98HOhqvnAPGALnkRRBqykb34umnT0OejW39NwShgGEJFE4AXgalUt917n3AzZ6+dJi8gpwA5VXRnsWEJAFHAYsEhVJwC7aXP5qQ99LlLxfHMeCWQCCex9iabP8sfnIJwShi/lRno1EYnGkyyeUtV/OIsLm04lnX93BCu+APohcJqIbMJzafJYPNfxU5xLEdB3Ph95QJ6qfuq8fx5PAumLn4vjgY2qWqSqdcA/8HxW+uLnoklHn4Nu/T0Np4SxAtjPmfEQg2cw65UgxxQwzjX6pcBqVX3Aa9UreMqqQB8pr6Kq16tqtqqOwPM5eF9Vz6MPlppR1e3AVhFpqkJ6HJ4KC33uc4HnUtQRIhLv/H9p6os+97nw0tHn4BXgQme21BFAmdelqw6F1Z3eInIynmvXkcBjqnpnkEMKGBH5EfAR8C0t1+3n4BnHeA4YBmwGfq6qbQe+ei0RmQpcq6qniMgoPGccA4AvgfNVtSaY8QWCiByKZ/A/BtgAXILny2Cf+1yIyK3A2XhmFX4JXIrn2nyv/1yIyNPAVDzl3AuBm4GXaOdz4CTUh/BcsqsCLlHVLivZhlXCMMYYEzzhdEnKGGNMEFnCMMYY4xNLGMYYY3xiCcMYY4xPLGEYY4zxiSUMY4wxPrGEYYwxxif/H9XNem9ryHnSAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Inicializo el grafico\n",
    "fig = plt.figure()\n",
    "ax = fig.add_subplot(111)\n",
    "\n",
    "# Dibujo el grafico\n",
    "grafico = Grafico(0, -10, 100, 10)\n",
    "grafico.draw_grafico()\n",
    "\n",
    "# Dibujo el cable\n",
    "cable = Cable(core_radio, cladding_radio, core_u=core_u, cladding_u=cladding_u)\n",
    "print('Angulo de entrada critico:', cable.get_angulo_critico())\n",
    "cable.draw_cable()\n",
    "cable.draw_critic_line()\n",
    "\n",
    "\n",
    "linea_nueva = Linea(cable, 20, [0,0], 'blue')\n",
    "\n",
    "linea_nueva_refrac = Linea(cable, linea_nueva.angle ,\n",
    "                           [linea_nueva.xf, linea_nueva.yf],\n",
    "                           'orange')\n",
    "if not linea_nueva_refrac.angle > 0:\n",
    "    \n",
    "    x_inicial = linea_nueva_refrac.xf\n",
    "    y_inicial = linea_nueva_refrac.yf\n",
    "    angulo_inicial = linea_nueva_refrac.angle\n",
    "    while x_inicial < 100:\n",
    "        linea_nueva_refrac_aux = Linea(cable, angulo_inicial ,\n",
    "                                       [x_inicial, y_inicial],\n",
    "                                       'orange')\n",
    "        x_inicial =linea_nueva_refrac_aux.xf\n",
    "        y_inicial = linea_nueva_refrac_aux.yf\n",
    "        angulo_inicial = linea_nueva_refrac_aux.angle\n",
    "\n",
    "\n",
    "            \n",
    "            \n",
    "ax.plot(0 , 0, label='Linea')\n",
    "ax.plot(0 , 0, label='Reflexion')\n",
    "ax.plot(0 , 0, label='Critico')\n",
    "\n",
    "plt.title('Fiber optic Cable')\n",
    "ax.legend()\n",
    "\n",
    "# Genero el grafico\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " Fuentes\n",
    "\n",
    " - http://www.ciscopress.com/articles/article.asp?p=170740\n",
    "\n",
    " - https://www.youtube.com/watch?v=pTakwkDtBDw\n",
    "\n",
    " - https://www.youtube.com/watch?v=MrJswUU143M\n",
    "\n",
    " - https://www.youtube.com/watch?v=02wPSDOXMhc\n",
    "\n",
    " - https://www.samm.com/calculating-fiber-loss-and-maximum-distance-estimates\n",
    "\n",
    " - https://www.youtube.com/watch?v=W8uoqsEWX1M\n",
    "\n",
    " - https://www.seas.es/blog/automatizacion/reflexion-y-refraccion-de-la-luz-en-transmisiones-de-fibra-optica/"
   ]
  }
 ],
 "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}