complex_numbers.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<figure>\n",
    "  <IMG SRC=\"gfx/Logo_norsk_pos.png\" WIDTH=100 ALIGN=\"right\">\n",
    "</figure>\n",
    "\n",
    "# Complex numbers\n",
    "*Roberto Di Remigio*, *Luca Frediani*\n",
    "\n",
    "Complex numbers are extremely important in mathematics and physics. The concept of a\n",
    "complex number was introduced as early as the 16th century in order to solve algebraic equations\n",
    "without a real solution, such as:\n",
    "\n",
    "\\begin{equation}\n",
    "    x^2 + 1 = 0\n",
    "\\end{equation}\n",
    "\n",
    "If we _define_ the _imaginary unit_ as:\n",
    "\n",
    "\\begin{equation}\n",
    "    \\mathrm{i}^2 = -1\n",
    "\\end{equation}\n",
    "\n",
    "then the quadratic equation above admits the two, purely imaginary, solutions $+\\mathrm{i}$ and $-\\mathrm{i}$.\n",
    "\n",
    "In general, a complex number is any entity that can be written in the form:\n",
    "\n",
    "\\begin{equation}\n",
    " z = a +\\mathrm{i}b = \\Re(z) +\\mathrm{i}\\Im(z)\n",
    "\\end{equation}\n",
    "\n",
    "where $a, b$ are real numbers and are called the _real_ and _imaginary_ part of the complex number $z$.\n",
    "Apart from all the usual operations, i.e. addition, subtraction, multiplication and division, we can furthermore define _complex conjugation_:\n",
    "\n",
    "\\begin{equation}\n",
    " z^* = a - \\mathrm{i} b = \\Re(z) - \\mathrm{i}\\Im(z)\n",
    "\\end{equation}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Exercise 1\n",
    "\n",
    "Write down explicitly the form of the sum, difference, multiplication and division of two complex numbers $z = a + \\mathrm{i}b$ and $w=c + \\mathrm{i}d$.\n",
    "Prove also that:\n",
    "\n",
    " \\begin{alignat}{2}\n",
    "  \\Re(z) = a = \\frac{z+z^*}{2};\\quad&\\quad \\Im(z) = b = \\frac{z-z^*}{2\\mathrm{i}}\n",
    " \\end{alignat}\n",
    "\n",
    "We can draw a powerful analogy between points in the plane and complex numbers. Let's rename the $x$-axis the _real_ axis and the $y$-axis the _imaginary_ axis. Any complex\n",
    "number can be represented in this plane, simply drawing the point with $x$ coordinate equal to its real part and $y$ coordinate equal to its imaginary part. This amounts to\n",
    "specifying the components of a 2-dimensional vector:\n",
    "\n",
    "\\begin{equation}\\label{eq:vector}\n",
    " \\overrightarrow{v} = \\begin{pmatrix}\n",
    "      \\Re(z) \\\\\n",
    "      \\Im(z)\n",
    "     \\end{pmatrix}\n",
    "     = \\begin{pmatrix}\n",
    "      a \\\\\n",
    "      b\n",
    "     \\end{pmatrix}\n",
    "\\end{equation}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Exercise 2\n",
    "\n",
    "Draw the following numbers in the complex plane: $z=\\frac{\\sqrt{2}}{2} + \\mathrm{i}\\frac{\\sqrt{2}}{2}$, $w=\\frac{\\sqrt{3}}{2} + \\mathrm{i}\\frac{1}{2}$ and their complex conjugates $z^*$ and $w^*$.\n",
    "The following example shows how to draw arrows using `matplotlib` and the `arrow` function. This function takes two four arguments, the $x$ and $y$ coordinates of the arrow start and the $x$ and $y$ coordinates of the arrow end."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEACAYAAABWLgY0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEFBJREFUeJzt3W2MpXV5x/HvtWy3whZpgIToIpt2kWwkRUtkJJW0BzV1\nIaQrxtaFBCvWdlOL8R2DL8xO0poFXhhijTVLVlJNzNpCE7dlCRjkhNDyMBp5EHdlkQb3AZEHbQRC\nMuxefXEOu2eHmTn3nLnP0/98P8lJ5p7zn/tc/DPz2x/3OWcmMhNJUplWDXsASVL/GPKSVDBDXpIK\nZshLUsEMeUkqmCEvSQXrGvIRsTMino+Ix5dY89WI2B8Rj0bE++odUZLUqypN/jbgo4vdGRGXARsy\n893AVuAbNc0mSVqhriGfmQ8Av15iyWbgW+21DwOnRcRZ9YwnSVqJOq7JrwMOdBwfan9OkjRkPvEq\nSQVbXcM5DgHv6jg+u/25t4gIf1GOJPUgM6OXr6va5KN9W8hu4FMAEXEx8JvMfH6xE2Wmt0y2bds2\n9BlG5eZeuBfuxdK3leja5CPiO0ADOCMifgFsA9a08jp3ZOaeiLg8Ip4GXgWuXdFEkqTadA35zLy6\nwprr6hlHklQnn3gdkkajMewRRoZ7cZx7cZx7UY9Y6fWeZT1YRA7y8SSpBBFB9vmJV0nSGDLkJalg\nhrwkFcyQl6SCGfKSVDBDXpIKZshLUsEMeUkqmCEvSQUz5CWpYIa8JBXMkJekghnyklQwQ16SCmbI\nS1LBDHlJKpghL0kFM+QlqWCGvCQVzJCXpIIZ8pJUMENekgpmyEtSwQx5SSqYIS9JBTPkJalghrwk\nFcyQl6SCGfKSVDBDXpIKZshLUsEMeUkqmCEvSQUz5CWpYJVCPiI2RcS+iHgqIqYXuP+MiLgrIh6N\niCci4tO1TypJWrbIzKUXRKwCngI+DBwGZoEtmbmvY8024G2Z+cWIOBP4GXBWZr4x71zZ7fEkSSeK\nCDIzevnaKk1+Ctifmc9m5hywC9g8b80vgVPbH58KvDQ/4CVJg7e6wpp1wIGO44O0gr/TrcC9EXEY\n+D3gk/WMJ0laiSohX8UXgccy89KI2AB8PyIuyMxX5i+cmZk59nGj0aDRaNQ0giSVodls0mw2azlX\nlWvyFwMzmbmpfXwDkJl5U8eaPcCXM/O/28f3AtOZ+cN55/KavCQtU7+vyc8C50bE+ohYA2wBds9b\nsxf4SHuYs4DzgGd6GUiSVJ+ul2sy80hEXAfcQ+sfhZ2ZuTcitrbuzh3AduC2iHgMCOD6zHy5n4NL\nkrrrermm1gfzco0kLVu/L9dIksaUIS9JBTPkJalghrwkFcyQl6SCGfKSVDBDXpIKZshLUsEMeUkq\nmCEvSQUz5CWpYIa8JBXMkJekghnyklQwQ16SCmbIS1LBDHlJKpghL0kFM+QlqWCGvCQVzJCXpIIZ\n8pJUMENekgpmyEtSwQx5SSqYIS9JBTPkJalghrwkFcyQl6SCGfKSVDBDXpIKZshLUsEMeUkqmCEv\nSQUz5CWpYJVCPiI2RcS+iHgqIqYXWdOIiB9HxE8i4r56x5Qk9SIyc+kFEauAp4APA4eBWWBLZu7r\nWHMa8D/An2fmoYg4MzNfXOBc2e3xJEkniggyM3r52ipNfgrYn5nPZuYcsAvYPG/N1cAdmXkIYKGA\nlyQNXpWQXwcc6Dg+2P5cp/OA0yPivoiYjYhr6hpQktS71TWe50LgQ8Ba4MGIeDAzn67p/JKkHlQJ\n+UPAOR3HZ7c/1+kg8GJmvg68HhH3A+8F3hLyMzMzxz5uNBo0Go3lTSxJhWs2mzSbzVrOVeWJ15OA\nn9F64vU54BHgqszc27FmI/DPwCbgd4GHgU9m5k/nncsnXiVpmVbyxGvXJp+ZRyLiOuAeWtfwd2bm\n3ojY2ro7d2Tmvoi4G3gcOALsmB/wkqTB69rka30wm7wkLVu/X0IpSRpThrwkFcyQl6SCGfKSVDBD\nXpIKZshLUsEMeUkqmCEvSQUz5CWpYIa8JBXMkJekghnyklQwQ17qwl+qp3FmyEsLeO2117j99tu5\n/PK/YuPGP+aNN94Y9khST+r683/S2HvttdfYs2cP3/zmv/GDH9zNmjVTvPLKw3zuc3/L6tX+qGg8\n+fvkNdHeGuwf4Le//UvgSlat+lfe855/50c/up81a9YMe1RNsJX8PnlDXhPr6NGjvP3tp3PkyNm8\n/voXgCuBM9v3Psipp36MJ554hPXr1w9xSsk/GiL1JCK4+eabOOmk51m16gDw++17Xubkk7fw7W/f\nasBr7BnymlgRwUUXXcirr77I0aP/yNq1fwI8yckn/zWf+cwn2Lz5L4Y9orRihrwm0tzcHFdcsZmp\nqSkADh06xM03/w1ve9sH2bDhBb7yle1DnlCqh9fkNXFmZ2ePhfv27TcyPX09Ea3LnS+88AKnnHIK\na9euHeaI0gl84lWqYG5ujiuv/AR33rkbgMOHD/OOd7xjyFNJ3fnEq9TF7Owsa9as4c47d7N9+40c\nPXrUgNdE8B0eKprtXZPOJq9i2d4lm7wKZHuXjrPJqyi2d+lENnkVwfYuLcwmr7Fne5cWZ5PX2LK9\nS93Z5DWWbO9SNTZ5jRXbu7Q8NnmNDdu7tHw2eY0827vUO5u8RprtXVoZm7xGku1dqkelJh8RmyJi\nX0Q8FRHTS6y7KCLmIuLj9Y2oSWN7l+rTtclHxCrga8CHgcPAbER8LzP3LbDuRuDufgyq8tnepfpV\nafJTwP7MfDYz54BdwOYF1n0euB34VY3zaULY3qX+qHJNfh1woOP4IK3gPyYi3gl8LDMvjYgT7pOW\nYnuX+quuV9fcAnReq+/pz1Rpstjepf6r0uQPAed0HJ/d/lyn9wO7ovXXkM8ELouIuczcPf9kMzMz\nxz5uNBo0Go1ljqxxZ3uXltZsNmk2m7Wcq+sf8o6Ik4Cf0Xri9TngEeCqzNy7yPrbgP/MzP9Y4D7/\nkPeEm52dZWqqdUVv+/YbmZ6+nlY3kLSYlfwh765NPjOPRMR1wD20Lu/szMy9EbG1dXfumP8lvQyi\nstnepeHo2uRrfTCb/ESyvUsr09cmL/XK9i4Nn7+7Rn3hK2ek0WCTV61s79JoscmrNrZ3afTY5LVi\ntndpdNnktSK2d2m02eTVE9u7NB5s8lo227s0Pmzyqsz2Lo0fm7wqsb1L48kmryXZ3qXxZpPXomzv\n0vizyestbO9SOWzyOoHtXSqLTV6A7V0qlU1etnepYDb5CWZ7l8pnk59QtndpMtjkJ4ztXZosNvkJ\nYnuXJo9NfgLY3qXJZZMvnO1dmmw2+ULZ3iWBTb5ItndJb7LJF8T2Lmk+m3whbO+SFmKTH3O2d0lL\nscmPMdu7pG5s8mPI9i6pKpv8mLG9S1oOm/yYsL1L6oVNfgzY3iX1yiY/wmzvklbKJj+ibO+S6mCT\nHzG2d0l1ssmPENu7pLpVavIRsQm4hdY/Cjsz86Z5918NTLcPfwv8fWY+UeegJbO9S+qXrk0+IlYB\nXwM+CpwPXBURG+ctewb408x8L/BPwK11D1oq27ukfqrS5KeA/Zn5LEBE7AI2A/veXJCZD3WsfwhY\nV+eQJbK9SxqEKtfk1wEHOo4PsnSIfxa4ayVDlc72LmlQan11TURcClwLXLLYmpmZmWMfNxoNGo1G\nnSOMNNu7pCqazSbNZrOWc0VmLr0g4mJgJjM3tY9vAHKBJ18vAO4ANmXmzxc5V3Z7vFLNzs4yNTUF\nwPbtNzI9fT0RMeSpJI2DiCAzewqMKk1+Fjg3ItYDzwFbgKvmDXAOrYC/ZrGAn1S2d0nD1PWafGYe\nAa4D7gGeBHZl5t6I2BoRf9de9iXgdODrEfHjiHikbxOPEa+9Sxq2rpdran2wCblcY3uXVKeVXK7x\nHa81s71LGiX+7pqa2N4ljSKbfA1s75JGlU1+BWzvkkadTb5HtndJ48Amv0y2d0njxCa/DLZ3SePG\nJl+B7V3SuLLJd2F7lzTObPKLsL1LKoFNfgG2d0mlsMl3sL1LKo1Nvs32LqlEE9/kbe+SSjbRTd72\nLql0E9nkbe+SJsXENXnbu6RJMjFN3vYuaRJNRJO3vUuaVEU3edu7pElXbJO3vUtSgU3e9i5JxxXV\n5G3vknSiIpq87V2SFjb2Td72LkmLG9smb3uXpO7Gssnb3iWpmrFq8rZ3SVqesWnytndJWr6Rb/K2\nd0nq3Ug3edu7JK3MSDZ527sk1WPkmrztXZLqMzJN3vYuSfUbiSZve5ek/qgU8hGxKSL2RcRTETG9\nyJqvRsT+iHg0It5X5bxzc3NcccVmpqamgFZ7v+GGaSKi8n+AJGlxXUM+IlYBXwM+CpwPXBURG+et\nuQzYkJnvBrYC3+h23klv781mc9gjjAz34jj34jj3oh5VmvwUsD8zn83MOWAXsHnems3AtwAy82Hg\ntIg4a6GT2d5b/AY+zr04zr04zr2oR5WQXwcc6Dg+2P7cUmsOLbAGYKLbuyQN2lBeXeMrZyRpMCIz\nl14QcTEwk5mb2sc3AJmZN3Ws+QZwX2Z+t328D/izzHx+3rmWfjBJ0oIys6dr2lWa/CxwbkSsB54D\ntgBXzVuzG/gH4LvtfxR+Mz/gVzKkJKk3XUM+M49ExHXAPbSu4e/MzL0RsbV1d+7IzD0RcXlEPA28\nClzb37ElSVV0vVwjSRpffXnHa7/ePDWOuu1FRFwdEY+1bw9ExB8NY85BqPJ90V53UUTMRcTHBznf\nIFX8GWlExI8j4icRcd+gZxyUCj8jZ0TEXe2seCIiPj2EMfsuInZGxPMR8fgSa5afm5lZ643WPxxP\nA+uB3wEeBTbOW3MZcGf74w8AD9U9xyjcKu7FxcBp7Y83TfJedKy7F/gv4OPDnnuI3xenAU8C69rH\nZw577iHuxTZg+5v7ALwErB727H3Yi0uA9wGPL3J/T7nZjyZf65unxlzXvcjMhzLz/9qHD7HI+wsK\nUOX7AuDzwO3ArwY53IBV2YurgTsy8xBAZr444BkHpcpe/BI4tf3xqcBLmfnGAGcciMx8APj1Ekt6\nys1+hHytb54ac1X2otNngbv6OtHwdN2LiHgn8LHM/Beg5FdiVfm+OA84PSLui4jZiLhmYNMNVpW9\nuBU4PyIOA48BXxjQbKOmp9wcmV81POki4lJar0q6ZNizDNEtQOc12ZKDvpvVwIXAh4C1wIMR8WBm\nPj3csYbii8BjmXlpRGwAvh8RF2TmK8MebBz0I+QPAed0HJ/d/tz8Ne/qsqYEVfaCiLgA2AFsysyl\n/ndtnFXZi/cDu6L1i4zOBC6LiLnM3D2gGQelyl4cBF7MzNeB1yPifuC9tK5fl6TKXnwQ+DJAZv48\nIv4X2Aj8cCATjo6ecrMfl2uOvXkqItbQevPU/B/S3cCn4Ng7ahd881QBuu5FRJwD3AFck5k/H8KM\ng9J1LzLzD9u3P6B1Xf5zBQY8VPsZ+R5wSUScFBGn0Hqibe+A5xyEKnuxF/gIQPsa9HnAMwOdcnCC\nxf8PtqfcrL3Jp2+eOqbKXgBfAk4Hvt5usHOZOTW8qfuj4l6c8CUDH3JAKv6M7IuIu4HHgSPAjsz8\n6RDH7ouK3xfbgdsi4jFaAXh9Zr48vKn7IyK+AzSAMyLiF7ReVbSGFeamb4aSpIKNxJ//kyT1hyEv\nSQUz5CWpYIa8JBXMkJekghnyklQwQ16SCmbIS1LB/h/ZCfrRTDND5QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9004881850>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "\n",
    "ax = plt.axes()\n",
    "# Use the arrow function\n",
    "# The arrow start at x=0.0, y=0.0 and ends at x=0.5, y=0.5\n",
    "ax.arrow(x=0, y=0, dx=0.5, dy=0.5)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The analogy between complex numbers and 2-dimensional vectors can be brought a bit further. Using Pythagoras' theorem we can calculate the length of a 2-dimensional vector\n",
    "(its _norm_ or its _magnitude_) in terms of its components:\n",
    "\n",
    "\\begin{equation}\n",
    " v = \\sqrt{a^2 + b^2} = \\sqrt{\\Re(z)^2 + \\Im(z)^2}\n",
    "\\end{equation}\n",
    "\n",
    "This lets us introduce a third way of representing a complex number, the _polar_ representation. Instead of specifying the cartesian components of our point in the complex plane\n",
    "(cartesian or _rectangular_ representation), we can specify the _length_ of the vector, $\\rho$, and its _angle_ with respect to the real axis  $\\theta$. We assume that\n",
    "angles increase in an anti-clockwise sense, with the zero corresponding to the $x$-axis.\n",
    "\n",
    "## Exercise 3\n",
    "\n",
    "Using the Figure below\n",
    "<figure>\n",
    "  <IMG SRC=\"gfx/trigonometry.png\" WIDTH=400>\n",
    "</figure>\n",
    "prove the following sets of relations between the rectangular and polar representations:\n",
    "\n",
    " \\begin{alignat}{2}\n",
    "  \\left\\{\n",
    "  \\begin{array}{lr}\n",
    "    \\rho = \\sqrt{a^2 + b^2} \\\\\n",
    "    \\tan\\theta = \\frac{b}{a}\n",
    "  \\end{array}\n",
    "\\right.\n",
    "\\quad&\\quad\n",
    "\\left\\{\n",
    "  \\begin{array}{lr}\n",
    "    a = \\rho\\cos\\theta\\\\\n",
    "    b = \\rho\\sin\\theta\n",
    "  \\end{array}\n",
    "\\right.\n",
    " \\end{alignat}\n",
    "The tangent function is defined as $\\tan\\alpha=\\frac{\\sin\\alpha}{\\cos\\alpha}$."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using the second set among the equations above we have:\n",
    "\n",
    "\\begin{equation}\n",
    " z = \\rho\\cos\\theta + \\mathrm{i}\\rho\\sin\\theta = \\rho(\\cos\\theta + \\mathrm{i}\\sin\\theta)\n",
    "\\end{equation}\n",
    "\n",
    "which we will rewrite as:\n",
    "\n",
    "\\begin{equation}\n",
    " z = \\rho \\mathrm{e}^{\\mathrm{i}\\theta}\n",
    "\\end{equation}\n",
    "\n",
    "The exponential above can be expanded as:\n",
    "\n",
    "\\begin{equation}\n",
    "  \\mathrm{e}^{\\mathrm{i}\\theta} = \\cos\\theta + \\mathrm{i}\\sin\\theta\n",
    "\\end{equation}\n",
    "\n",
    "known as _Euler's formula_ and $\\mathrm{e}=2.71828\\ldots$ is the _Napier's constant_. This formula is extremely important both from a theoretical and practical point of view.\n",
    "\n",
    "## Exercise 4\n",
    "\n",
    "Given $z=\\rho \\mathrm{e}^{\\mathrm{i}\\theta}$ how would you write the Euler representation of its complex conjugate $z^*$?\n",
    "\n",
    "    Hint: what happens to $a+\\mathrm{i} b$ when taking the complex conjugate?\n",
    "   \n",
    "Draw some more vectors in the complex plane using `matplotlib`: $\\mathrm{e}^{\\mathrm{i}\\theta}$ where $\\theta$ is $30^\\circ$, $150^\\circ$, $210^\\circ$ and $330^\\circ$.\n",
    "What do all these vectors have in common?\n",
    "You can write polar vectors in `matplotlib` in two ways. Either by using Euler's formula to find the arrow start and end coordinates or by using the proper syntax for polar plots.\n",
    "We invite you to use both and define the appropriate functions for the calculations needed when using the Cartesian representation.\n",
    "What follows is an example of polar plot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAENCAYAAAAha/EUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd8lFX2/z93MkmmJhBQehVBEQERBFFXVLoIFlDEVVd3\nXXd1XXT9quv+3LUgIsq6qBRBUKRXQUCpUkMJLSGhl5DeQ8r0en5/THESMsnM5LnzzITn/XrdV+aZ\neZ5zTpLnOXPvueeey4gIEhISEkIhE9sACQmJpoXkVCQkJARFcioSEhKCIjkVCQkJQZGcioSEhKBI\nTkVCQkJQJKci4RfG2GTGWIa7/d39XnPG2HbG2HnG2DbGWKLP+QsZY6mMsdHiWS0hNpJTkagTxtht\nAP4IoD+AvgDGMMZuAvBPADuJqAeAXQDe9Tk/x33+86IYLRERSE5Fwh+3AkghIgsROQDsA/A4gLEA\nfnCf8wOAR92vHQDUAOIASBmV1zGSU5HwxykA97mHOyoAowF0ANCKiIoBgIiKALRyvz4HIBbAXgBz\nxDFZIhKQi22ARGRCROcYY9MB7ACgB5AKV2+kNk6fa94Ik3kSEYzUU5HwCxF9T0T9iWgIgEoA5wEU\nM8ZaAQBjrDWAEhFNlIhAJKci4RfG2A3unx0BPAZgOYCNAP7gPuV5AD+JYpxExMKkVcoS/mCM7QOQ\nBMAG4A0i2sMYSwKwGq74SjaAJ4moUkQzJSIMyalISEgIijT8kZCQEBTJqUhISAiK5FQkJCQERXIq\nEhISgiI5FQkJCUGRMmolAoIxxuC6XzyNAbB7GhE567lc4jpCmlK+TnE7iSQAbQG08fxUq9Vd4+Pj\nuwBobbPZWplMpkS73S4HwGQymZMx5oyJiSHGGJxOJ7Pb7TIikjHGSCaTOeLi4gzx8fElMTExeTab\nLae6uvoygEJ3K3D/LHUvUpRogkhO5TqAMdYCwJ0ymWxAQkLCgzab7Xaz2ZwUHx9v12q1tlatWlGP\nHj1Yp06dFElJSfJOnTqhY8eOaNOmDVq1agWFQoGYmBi4/NC1EBGcTicsFgtKS0tRWFiIgoICFBYW\n4vTp087i4mJTaWmpPT8/nxUVFcWZTKZYlUpVpVAozlRVVe212+1HARwHkE/SDRn1SE6lieHrQJo1\na3a/2WzuZ7fbE3r37m289957VX379o0dOHAgOnXqBKVSKYqNVqsVBQUFSEtLw+bNmx2nTp0ynjlz\nJtZut9sUCsVpnU73q+RoohfJqUQ5jLFmAEY2a9Zsot1uv8dmsyV07tzZOmrUqPi77ror9s4770S3\nbt0gkwkbk9+zZw+GDBkimDwiQl5eHlJSUnD06FHHoUOHDGlpaXFWq9WhUqmOVFRUrACwmYgKBVMq\nwQXJqUQhjLEujLGxzZs3f0av1/e55ZZb7K+88orqgQce4OJA6kJop1IXHkezb98+rF271rB161Z5\nbGxsvslk+sFut28AkCH1YiIPyalEAYwxGYABCoXi8bi4uKfsdnvrxx9/3Dl+/Hjl0KFDoVarxTYx\nLNhsNuzbtw/r16+3rlmzxmYwGEwANhgMhlUA9hGRVWwbJSSnErG4Z2cGJiQkvGqxWB5t3bo1Pfnk\nk8rHHntMftdddyEmJkZsE0WFiHD69GksX77cuWnTJv3FixdjFQrFrqqqqjkAtkmzS+IhOZUIgzGm\nBvB0QkLCP9VqdevJkycrx48fL7vpppvENq0G4Rj+BENRUREWLVqEb775xlheXq43m83/s9vtC4io\nTGzbrjuISGoR0AD0UKvVcxQKheHuu+82bNiwgRwOB0Uqu3fvFtuEOnE6nZSSkkJPP/20MS4uzqLV\natcCGAT3F6jU+DeppyIijDE5gLGJiYnvOhyOXq+++mrMX//619hOnTqJbVqToLy8HN99953ziy++\nMJlMpsKqqqrpAFYQkUFs25oyklMRAcZY89jY2NdjY2Nf69Gjh/ytt97SPv7444iPjxfbtCaJ0+nE\ntm3b8O6775rOnTtHcrn8B4PB8CkR5YhtW5NE7K7S9dQAqOLi4v6lUCj09957ryUtLY2ilUgd/jTE\nlStX6M0337QqlUqjSqWaDaAlRcC90ZSatEo5DDDG5DKZ7GWlUpk/bNiwf6Wmpqr3798f16dPH7FN\nu+7o3LkzZsyYEXvp0iXl448//qJSqcxSKpUfMMY0YtvWVJCGPxxxTws/odFoZnbv3r3ZrFmz1Hff\nfbfYZkn4cOnSJbzzzjvGX375xWGz2f7lcDjmk5Tv0igkp8IJxthDGo1mbrt27dp89dVXmmHDhvld\nkCchPsePH8cbb7xhSE1N1ev1+n8AWElSOYfQEHv81dQagDsSExMPtGnTRr9gwYKInhZuDNEaU2mI\nXbt2Uffu3fVqtToTwAiKgHsq2poUUxEIxphCpVL9V61WH/zkk0/uzs7OVv/xj38MyzocCeF44IEH\ncO7cOfXixYu7tGnTZl1CQsI6915HEgEiDX8EgDF2l0ajWT1w4MAbli5dqmrdurXYJkkIgMFgwNtv\nv2357rvvLGaz+TkiknZjDADJqTQCxli8Uqn8RC6X/3XevHmKiRMnMilu0vTYtWsXnn32WaNOp9ui\n0+n+TERXxbYpkpH65iHCGBugVqvP9+vX75WLFy8qn3766evKoezZs0dsE8LGgw8+iAsXLqieffbZ\nMQqF4jJjbKzYNkUyklMJEsZYvEqlmpGQkLD322+/7bh//35Fq1atxDZLgjNqtRqzZ8+O37RpU7NW\nrVqtSEhIWCvFWvwgdqQ42AZgIYBiAOk+730E4CSANAA7AbR3v98JgBHACXeb43PNGPc184PQPUCt\nVmf379/fXFRURBLXJ3q9nl555RWzUqm8CmAshXYfywCkAtjoPn4fQJ7PvTqSat7zqQBGh6Ir3E10\nA0L4Z9wLoG8tp6Lxef0agAX0m1NJ9yNnpfsf+xGAng3plcvlf9ZoNMbly5c7nU4nNWV0Oh2ZTCbv\n8Z49eygnJ8d7vGHDBrp8+bL3+JdffqnxeVFREZnN5vAYKyI///wzNW/e3KRSqb4AEEPB3cdvAFha\ny6n8o47zbnN/FgNgVTA6xGpRN/whomQAFbXe0/scqgH41tDwF+hgAOIAqADY/OljjMVqNJpvW7Ro\n8eXx48ebROyEiGCz/fYrr1mzBqdOnfIe7969G4WFv5WCvemmm9CiRQvv8UMPPYTMzEzv8cCBA2t8\nnp6ejuLiYu/xokWLcOXKFe9xYWEhnM7ozysbPXo0Lly4oOjdu/eftVrtDsZYQiDXMcbaAxgNYEHt\nj+o43QHXPR0HIDpmVcT2aqE01NEDAfAxgBwA5wAk+pyng6s7uRvAvT7nDwVwDMD0evS00Gq1h4cM\nGWKorKykaKWiooIKCgq8xxs2bKD09PRGyWxM8tvGjRspPz+/hizfnlG0YbVa6cUXXzSrVKpcADdT\nw/fvGrh62/ejZk/lClxD+AUAmvmc/z8ARwD8riHZkdBENyAko+sf1rwD4Hv36zgAzd2v+7mdjiZA\nHb1UKlXhP/7xD4vdbq99H0U0RqORsrKyvMdHjx6lCxcuiGhR/Rw/fpyMRiMRuYosrV+/PiozkWfP\nnu1QqVTVAIaR//vqYQCz3K+H+DiVG/BbisfHABb6kxHpTXQDQjK6fqfSAa4q63V9thtAvwDkj1Mo\nFIbZs2dHTfDEtyd1+vRpSk1NFdGa0HE6nZSWlkaeuJXRaGx0ryqc7N27lxISEoxxcXH/QB3V5gB8\n4v5yy4Rrt0Y9gMUU4P0dDU10A0IyGujs6zgAdPN5/RqAJe7XLQHI3K+7Asj17VbWIZcplcr3k5KS\njCkpKXXdMxFJbm4urVq1Kqw6w7X2x2Qy0c6dO73H0RAkz8zMpLZt25o0Gs0yAPHk/37zHf609nn/\nDQDL/V0X6U10A4I2GFgO1568FrfHfwHAWgAZcE27rQNwo/vcxwGccsdUjqGeKTkAKo1Gs6l79+5G\n3/F+JOJ0OunDDz8UdYgg1oLC1NRUWrdunSi6g0Gn09HDDz9s0Gq1aQBaUcNOZTGAdHdMZYO/a6Kh\nSWn6ABhjWq1Wu/uBBx7ouXLlSqVY24HWR1FREeRyOVq2bCm2KRHF5s2b0a9fP7Rt21ZsU67B6XTi\nvffes82cOfOqyWS6i66T8pVRN6UsNIyxZhqNJvmJJ57otX79+oh0KACQkZFx3e/1UxcPPPCAaHtC\nN4RMJsMnn3wS+9Zbb7VUq9VHGWNdxbYpLIjdVRKzAWihVqvPP/nkk5ZIG6ufP3+eVq9eLbYZfonE\neioOh0P0YaE/3DNDZQB6UATc+zyb6AaI9osDN2o0mktvvPFGxDkUIiKz2RzRQclIdCq1MZlMEfU3\n/Pzzzx1yeZwJwG0UAc8Ar3ZdDn8YYy3UavXBv/3tbx3/+9//xkVKhuz06dNhMpkAAPHx8RFdfjKS\ndif0x9mzZ7Ft2zaxzfDSrl07md1uVcjlimOMsR5i28OL6y5Q646hHBozZsxNy5cvjwVcxXg0GvGL\nqRNRRDuSaKeiogLNmzcXRffixYvx/PPP4+GHH8HAgQOc06ZNq3QHby+LYhBPxO4qhbMB0Gq12pN/\n/etfzZ5usclkoi+//JLEID09ndasWSOK7sYSDcMfX5xOJ3399ddks9nCrnv58uUEgB56aJj3vY8/\n/tihUCjKAHSkCHg2hGyiGxC2XxRQajSaYy+88IIpUsbZkWJHKESbU6lNuJzLsmXLCACNGfPINZ/9\n97//tavV6gIAbSgCnhGhmugGhOWXBJhGo1k/aNAgc30zA+Xl5VRSUuL3cyFYuHBhjcV9EuKwfft2\nOnToEFcdP/zwg1+H4uE///mPVaPRnAKgpAh4VoRoohsQjqZQKN677bbb9J5Fa/6orq6mDRs21HtO\nY7FYLFzlS0QGvw15Hqr3PKfTSXfffbdZrVavRR1rhaKxiW4A918QGNO8eXNDXl6e338sb7Zt2xbV\nQ526iPbhjy+5ubn0/fffCybP41Dq66H4YjAY6JZbbtErFIp3KAKemcY20Q3g+ssBPRUKhWH//v3+\n/6N+SE1NFWQo5HA46ODBg42WE2k0JadC5KqJIgSBDHnqIicnh5o1a2ZElJSMrK+JbgC3XwxIUqvV\nBYsWLQqpi1BdXd2o8gFC3aQS4cVut1NycnJI1wY65PFHcnIyxcfHGwDcQhHwDIXaRDeAyy8FyDUa\nzYHJkyeLUii1vLycZs2aJYZqCQHYsWNH0NcEO+Txx4IFCxwajSYX7uJi0dhEN4BHU6lUc3v16mUW\nqmLb+vXruc8KRRtNbfjTGEId8vjj1VdfNWs0mn0A5BQBz1OwTXQDhG5yufxPHTp0MFRUVDT4zwuU\n6upqunr1ar3nWCyWkL7hogW9Xl9j5mrt2rXU0GxatON0OmnKlCn1Btk9PZShQ4f5PSdYbDYb9erV\ny6RSqWZRBDxTwTbRDRD0lwH6aTQaw7lz5xr+zwlMcXFxjW0roo3Lly9TeXm593jRokWUmZnpPV65\nciX5zqBt2bKFCgsLvccLFy6sURd3/vz5lJub6z2OxJXDgVBfb1eoIU9dlJWVUdu2bQ2MsScpAp6t\nYJroBgj2iwBxKpUqe9asWVznbmfOnEl6vZ6nirBw4MCBGrVfDx48SKWlpdz0zZo1q8YQMhoTAH17\nv55M2YcfHsNN35EjR0ipVFYjyqrAiW6AUE2lUn06dOhQA+98EN+SBFVVVTR37lyu+oTi/PnzNTJI\ni4qKGjVD1diYynfffefdlsPpdEb8Fh0Oh4NmzpxJRL/1UEaNGs1d7+TJky1arXZ7NCXGiW6AIL8E\ncKdWqzWG+9vPZrORwWAIq85AcTqdNYYrBQUFgvawhAzUmkwm+uKLLwSTx5OlS5dyG/LUhdlspo4d\nO+oZYxMpAp61QJroBjT6F3APexYvXhy2lFWn00lOp5M++OCDiM2UzczMpO3bt4ttRkicO3eOli1b\nJrYZ1+DpoQwfPjys24ZE2zBIdAMa21Qq1af3338/92GPB0/JwkjD6XTSRx99FLUB0dr4/j/z8/NF\nd96+QVmn00k///xzWPW/9dZbFq1Wuy0ahkGiG9Ao40Ua9tS+wZ1OJ4mxLarT6ayxEXo4H7xw5qnU\nDiqHG6HzUELBbDZTmzZtjAAifhgUteUkGWNxarV67ezZsxVt2rThrs/pdHoc2TXV2cxmM5YvX87d\nhtps2LChxkbpTbVq3ODBg3H77bcDAOx2OzZs2OD9X/BmxYoVeP755/HQQ0OxadPGaz4/deoUKisr\nudsRHx+Pn376SalUKuczxlpxV9gYxPZqoTalUvlp//79w1Zwac6cOVRcXBwWXfURCTaIzcmTJ8Oi\nxzPkeeSRsX7PKS8vpyNHjoTFHiLvMGgrIngYJLoBIRkN9BZj2BMIRUVF3FL69Xq9oEv0mwJ79uzh\nkvMSCUOeujCbzdSyZUszgPEUAc9iXU10A0JpCQkJe7/88suwRCSD7QlVV1cLHsTzjZtECpGy9ken\n01FOTo6gMkNdbXzs2DFB7fDHjh07SKPRFACIpQh4Hmu3qIupMMbuZ4wNfPnll7nbXlJSgrlz5wZ1\njVarxejRowWzYceOHUhPTxdMXlNDo9GgQ4cOAIDq6mpcvty44vQrVqzApEmTMGbMI9i5c2dQ1xYW\nFsLhcDRKfyAMHToUffr0SZDJZH/iriwUxPZqwTQALCEhIWPRokXBu/cQacwUbUpKirS6OYxYLBba\nvHlzyNd7hjz1xVAihWPHjpFCoagEoKYIeDZ9m+gGBGUsMK5bt266aMnFqK6uplOnTgV93TfffMN1\nHc71QjD3SWMLLPlSWloaliHrAw88YIyLi3uPIuDZ9G2iGxCwoUCMRqPJacw3UaAsXryYysrKuOvx\nRzQUx46UmIo/nE4nffjhhwHFxIRebZyZmUlbt24VRFZ9nD9/npRKpR5AEkXAM+ppohsQsKHAC7fc\ncktYMmd9SwAIxerVq+vtfWRnZ9e7zD7SiHSnEiiROssTKC+++KJJpVJ9Sdc+LwsBFANIr/2Zzzlf\nAbgIIA1AX3/nBdtEdxYBGQko1Gp12YEDBxr1DxCT6upq0ul0fj9fvHixKLvnXS/s3bv3mvdWrFjh\nXsszgpte3l+C+fn5FBcXZwbQnmo+M/cC6OvPqQAYBeBn9+uBAA7XdV4oLSpmf+Li4v52zz33KAYP\nHsxVz969ez1/cMHRarX17tf87LPPQi6Xc9Et4frytNls3uOVK1fi6aefxpgxj2Dbtq1cdJrNZnz+\n+edcZHto27Yt/v73vzOtVvup7/tElAygop5LxwFY7D43BUCiYJm6QnknXg1AYnx8vCGUgGew/Prr\nr9x1EBF9/vnnZDQao7okYzQPfzwxFJ4FljyEY0hbUVFBGo3GBOBWqvnsdIL/nsomAIN9jncC6FfX\nucG2iO+pxMbGvvLwww+z2267jbuuBx98kLsOAJg8eTKUSiV69+4NpVIZFp0SLpYvX45JkyZh1KjR\n2Lx5E3d9MTEx3HU0a9YMb775ZqxWq/2Yu7JAEMIz8WpwzfiUHD16tLHOvF7ErDoWTcHZaMdTAnL4\n8OFh1817O93S0lKKi4uzAGhJgfVUvgHwlM/xOQhUryXSeyoj27dvr+jfvz83BU6nEzNnzuQmvzab\nNm3CpUuXALgc+ieffOL5p0pwZOXKlXjmmWcwatRobNu2Lez6PVm/vGjZsiWeeOIJh1wu982yZe5W\nFxsBPAcAjLFBACqJqFgQY4TwTLxaQkJC8ldffdU4Fx5hCLl1CG+2bNlCRUVF3uO5c+d6F+/t3r27\nxjER0bx580TN7/GHZ5andqZsXl4eLV++XCSrhCclJYU0Gk0xgBgAywEUALAAyAHwAoCXAfyZfuud\nzAJwCcBJCBRPIYrgKWUAXTUajSlaA5nBYrfbG9xbiDdLliyhCxcueI9zcnL8FseuK1BrtVprDOdm\nzJhBVVVVgtsZDA2l3ofbPp47MTidTurYsaMBIu/HLLrz8NeUSuX/Xn/9da6ppatXr+Yp3ktFRQX9\n8MMP9Z6j1+tpwYIFYbHHw86dO2tU2OeJ3W4POMNVKDw9lIceGho2nQ3x2Wef1fn+li1bqEePHnTz\nzTfTp59+es3nZWVlNHLkSOrTpw/16tXLbwmM2bNnU2Ji4h6SnEotowBlfHy8/uLFi3X+4YTA6XSG\ntC4nFCwWi+jf2ESu3zk/P997LOZygCtXrtCWLVu4yfc35PHHkiVLRNvPyeFw0E033URZWVlktVqp\nT58+dPbs2RrnfPDBB/TPf/6TiFxB2aSkpDqTJQ0GA6nVahOALiTS8xupgdoJPXv2ZN26deOmgDGG\ncExTA0BcXBwSEhICPj8vLw+lpaWC23HlyhWcOnWqhl2hsmfPnkbZ0qlTJ/Tu3btRMvyxZMkSPP30\n03jkkbHYuPGngK4ZNWqUaMmHR44cwc0334xOnTohNjYWEydOxE8/1bS7devW0Ol0AACdTocWLVrU\naa9KpcKLL74oUygUfwuL8XUQkU4lMTHxnQ8++EDFS77dbuclugZ79+5FWVlZ0NclJibi5MmTgtjw\nyy+/wGAwAAC6du2K4cOHCyK3sTDG0LZtW+/x6tWrcebMmUbLXblyJZ577jk89NDQgB0KALRo0QLx\n8fGN1h8IJpMJu3fv9h7n5+fXmB1q37498vPza1zz0ksv4fTp02jbti369OmDL7/80q/8V199NY6I\nXmaMKYS3vmEizqkwxu6Mi4vr/PDDD3PT8dlnn8HpdHKT76FZs2Zo2bJl0NdptVoMHTpUEBs6dOgA\nlUp4/zxkyBBB5T355JO49dZbGyXDN/V+584dIcnIzs7GxYsXG2VHQyiVyqD/J9OmTUOfPn1QUFCA\n1NRUvPrqq9Dr9XWe26NHD/Tv358AjBfA3KCJOKeiUqn+OmHChHiemYjvvvsuZDL+v3qfPn0aLePA\ngQNBDYUsFgv279/vPb799tujpsq+x87y8nJ8/fXXQV3rO+Spq+p9oLRt2xbl5eUhXx8oAwcO9L5u\n164dcnJyvMd5eXlo165djfMPHDiACRMmAABuuukmdOnSBefOnfMr/4033tA0b958ssBmB4ZYwZy6\nGgCmUqnKxNzjRQgyMzMFm+Worq6uMc3bEHl5eXT58mVBdNcH77U/wfz9PEHZoUOHcbSID06nk+x2\nuzdQa7FYqE+fPnTmzJka5/3jH/+gDz74gIhcxdXbt29fb4kOg8HgybBtTuF+jsOtsF5jgL6tW7fW\n8Zp2LCkpoczMTC6yfRGj4n24Z3LCuaBw+fLlfqu4BTvLEygOh4P7BnE7duzwTulv2bKFunfvTt26\ndaNp06YRkasC4Lx584jINeMzZswY6t27N91+++0BJe317dvXAGASXc9ORS6Xf/Daa69xezqOHj3K\nZTuHcLFs2bI6M1Z//fVXSk5OFsGi8HDx4sU6ncqSJUu4FViyWCw0a9Ysv583lFdC5HK8ffv2pdtu\nu42GDBlyzed2u53r2q9vv/2WmjVrtomuZ6ei1Wozd+3aJcCfs2mi0+lEXfwYCXh6sZ4eyogRI8Nu\nQyB5JZWVldSzZ0/Ky8sjIhKl5nBhYSHFx8cbAcRRGJ/jiAnUMsba2e32tvfee6/YpoTMjBkzuMrX\naDRQKFyzhDqdTtStOxqbpxIqU6dOxdKlS71B2a1bt4TdhkDySpYvX44nnnjCG3CtbxYwNzeXi52t\nW7dGx44d7QDu46LADxHjVACMGT16tD02NpaL8CVLlnCR68vLL7/MXYeHf/3rXyFNV0c7N910E559\n9lmMGfNIUHkooWKxWLBq1aoa7wWSV3LhwgVcvXoVDzzwAAYMGFDv/bd1K5/KcwAwbtw4tUqlmsBN\nQR1EjFNp3rz57ydMmKDmJf/uu+/mJdqLVqvlrsPDzJkzaySPhRuh81QCYeXKld6NvjZt2gij0chd\nZ3x8PO68886gr7Pb7Thx4gS2bNmCrVu3YsqUKd6SF7V56aWXGmumX5577jmZTCZ7nIUxryAinApj\nTKPT6e4aOXIkNx08U/5zcnJq1D/lhcPhwC+//AKgZkWxcOgWG8+QZ+xYVx6KzWbDN9980yiZTqcz\noIzn2vdOIHkl7du3x4gRI6BQKNCiRQv87ne/EyxLOhh69eoFlSvTrlfYlIYzgOOvAXjs7rvvFn/F\nXYgsWbKkUTsZBkplZSWdP3++xnsOh4M+/vhj7rprE84pZU9QduRI4YKyNpuNJkx4liZNeiGg861W\nq7cMRCB5JWfPnqWhQ4eS3W4ng8FAvXr1otOnT/uVn5yczO0eevrppy1yufzfdD3N/iQkJKyYOXMm\nl+SU8vJymjt3Lg/R1zXhcioehzJ2rP88lGCnZS0WC3Xo0Jnk8jgyGAwBXZORkVFjS9WG8kqIXAXO\ne/bsSbfffjs1VGwsJSWFWwGvjRs3UrNmzc7SdeZU8tPS0gT6E9bE6XSGZQtKnly8eDGgAk4Wi4XL\nRmhisXjx4gYdCpGr/mtqampAMq1WK/3rX+8RAAIgemGscGA0Gik2NtYKQEHXg1MBoJXL5TZ/FcYi\nmdLS0rB8Y2/ZsiWgjcZ0Ol2DxaCiBU8PZdiwxqfeexLVWrdu7XUmffv29+aQeDhy5AjJ5XJat25d\no3VGGl27dq0CMICuE6fyu9tuu41LPrTVauWaLFZYWEjFxcXc5EcyPJ1pIEOeQHE4HNSuXTuvM5HJ\nZPTzzz/Xed6DDz5IDz/8cL1OhefODidPnrwmiU4oRowYYQLwFwrDMx0Jsz93Dh48mEvdhxMnTiAl\nJYWHaACu5KIbb7yRm/zGzOpkZmZyKfTEG89q47Fjx16TUNYQeXl5NXI+zp07h5iYGG8OyYEDBzB1\n6lRkZGRcc+3XX3+N8ePHN/j/LCws9HwZCk7Hjh3r3cWyMQwfPlyh1WrDk1kaDs9VX0tMTPwp3LVZ\no4VPPvkk5LUh1dXVtGfPHoEtIpo/f36NgOKqVasEK8P422rj0GvK5ubmUmFhId1yy23e3snQoUO9\n6f1Lliyh1157rcY1+fn53rU5f/jDH5rk8OfQoUPUrFmzS3Q99FScTucAnvv68OKLL76Aw+HgquPd\nd98NeYeW2ZEJAAAgAElEQVQ7rVaL+++/v9E2zJkzB8XFv20HM2nSJCQmJnqPBw8e7C1LSUSYOnVq\nSD0sT4GlsWPHYseO0Aos6fV6/P3vb6BNmzY4d+40vvjif1i9ejW6dOlSb02Z119/HdOnT/ceE6ee\niJj06dMHBoOhY1iqwYXDc/lrALSxsbFWHkFas9lMJSUlgsv1oNPpuMkWml27dgX8t9i6dWtAcQN/\nMRXfXAu9Xh9QLybQWR5/2Gw2+uc//+Xtmbz00sveBXyHDh2iESNGeM+dNm3aNauKu3TpQl26dKHO\nnTuTRqOhVq1a0U8//eRX386dO7lN/yYnJ9ebz9IYWrdubUAYgrViO5Xf9ezZk0uQ9tSpU1FbDuDg\nwYNUXV0tmLzq6mq6cuVKQOcGGtgOJFBbUlJCK1eurPecxszyOJ1O+vbbBV5nMnz4KKqsrKTs7Gz6\n5ZdfiCiwRDVfAhn+5OXlcVt1rNPpuE1zjx8/3oAwBGvFdipvTJw4MeqSSHiXH9i/f3/Y9sex2Wx+\n96IRmtqJZo2Z5dm5c6fXmXTs2LneOjmBJKp5eOGFF5pkTIWIaNasWZSQkLCUmrJTSUxM3DBjxgwB\n/2zhYerUqWKbEDKLFi26JkEuXDVa/ve//3kTEZcuXVqnQ2mo+NFnn33mdSYA6h2mSNQkXMFasZ1K\n9okTJwT8s/1G7TUyEi50Oh1ZrdYaeySHQmPyVDw9lFGjRtV4v77iR0VFRXTjja28zmTPnj20ZcsW\nGjhwoF89vmn1QpOcnEy8ainPmzePeGz3azQaSSaT2QHIieNzLersj9lsvqFjx46Cy3U4HDhx4oTg\ncnlTWVmJ77//nqsOjUYDp9OJzZs3c9XjD88sz+DBg70rrj3UVfxozZo1eOCBh9C6dWuUlBRj6dKl\ncDqduP/++zFo0KBr6pj40rZtW24ruPv27VujpoqQPPXUU+BRV0ipVEKtVlsB8EuuAsTrqQBQxsTE\n2MO5t64QNPYbvj6sVqugAdqG+PjjjwNK/xcKTw9lzJgxtGfPnmtycNauXUsvvfQSEbliPQ89NNTb\nM/nooynXrOL9/PPPvedLBEa3bt0qAfSnpjj8AdC1efPmgS0RjSBmz54ttgkhU1FRUaPqfjjKNXj4\nreq9/yLVa9eupT/96U80c+aXXmfSo0ePOncK2LVrF/Xs2fO6WBAoJPfff38lgEeoiTqVe3v16sVl\nOjka9w2y2+3cH/IlS5b43XYi2JXcwcRUPFXv65rlycrK8vaWZsyY4XUm/fr1p/fff7/OYO3Jkyep\nW7dudOnSpQZ179+/nzIyMgK2NRh++OEHblPL06dP5yJ3+PDhJgB/pibqVCaMHDmSS2GmNWvW8BDL\nld27d9O+fftE0W2z2fxuM+GPQJ3Kb1XvR9T5+ZkzZ7wzQQCIMUYpKSl+c0qys7OpW7du3v1yGsJg\nMHDbv0ev13MbPgZa5yVY3n33XSeAD6mJOpXJkyZNiqr9JoqLi7n9s4mC25UvGmgoDyU/P5/kcrnX\noaSnpzeYU/KnP/2JkpKS6I477qC+ffvSgAEDwvb7NAXmzJlDiYmJi6kpOpXY2NjP3333XUH/YLxZ\nv359VBZBWr9+fcBp+kajUZAufX2p93q9nnr37uN1Jlu3bm20PonA2LBhA7Vo0WIfcXy2RZtS1mg0\nXW+++WbB5TqdTm774Tz66KNISkoSXK7T6URVVZXgcj30798fN9xwQ0Dn2u32gBb01bfvz8qVK/Hc\nc89h2LBhNcoXOBwOTJr0DDQaDdLTT2Lu3LkgIowYMQJLly6FxWIJyMZgWbt2LTIzM7nInjZtGhe5\n69atw/nz5wWX26ZNG9hsNuHzOHzh6bHqay1atDjC4xtKr9fTtm3bBJfLk9zcXFq9erXYZgSFv5iK\nZ8gzbtw473tOp5OmTJni7Zm88cY/rglKFxUVccvsNZlM3GIfPJLUiFzpBTwC97m5uaRSqSqoKQ5/\ntFpt3smTJwX9g/EmkNmGSKMxcZpz584FtdK7riHPqlWrvM5kxIiRYd9IXqImVquVZDKZA0AMNbXh\nj81ma9aqVSux1IfEgQMHxDYhaKZMmRLytW3btsXFixcDOtcz5Bk6dCh++uknHD58GIwxPPXUU2jX\nrj0qKiqwdesWb+0VCXGIjY2FUqm0AmjBTQkvb9VQi4uLM/GY6svOzhZlM+zGwLPObaiV4xrCd/jj\nO8tz5cqVGgv+Ai25QORaq8Or98prEejy5cu59GBtNpt39kto1Gq1GUB7Emv4A2AhgGIA6T7vDQBw\nBECq+2d/n8/eBXARwFkAw33eHwPgJID5RITY2FirUGUIfTl+/DhlZ2cLLpcns2bNEtuEBtm2bVuN\noZDHqXiGPKNHj65RYPrw4cNB67BardycIK9tWnjFPoiI21CxRYsWegCdyf8zPxLAOQAXALzjfq8L\ngBQAOwEk+ruWAnQq9wLoW8up7PY4DACjAOx2v+7pdjRyAJ0BXALA3J+thGub1Y8A9GSM2aNpPx6b\nzUaZmZlimxEUQt7w1dXVlJ+fX+M9Tw8lKSnJ60yiLeB8PdK6dWsdgG5U9/Mucz+3nQDEup/nWwF8\n7n7vAQCv1nWtpzUYUyGiZAAVtd4uBOApVNoMgGep6FgAK4nITkRZ7h7LXe7PGIA4ACoANiKS8ViJ\nyQuTyYS0tDSxzQiKxYsXo6SkRBBZWq22xobwK1aswNNPPw0AuHr1KqZOnQoiwoQJEwTRJ8GPmJgY\nguuLvy7uAnCRiLKJyAZXZ2AcADsAjbtZ61VQn8fx8V6dULOn0hFALoAc988O7ve/BjDJ57wFAB53\nvx4K4BiA6QBkjDEuXjgjIyPqdiTkufKZB++++663Z/L73/9e0O4/r32hp0yZwkXukSNHuCXv8bK5\nZcuWBgC9qO5n/Qm4QxTu498D+ApAOwB7AGwAoKrrWk/z560aYiGA14hoA2NsPIDvAAxrwHntBNAf\nABhjNUrEexKphgwZ0ujjwsJC5OXlQaFQCCIvHMfTpk3Do48+GjH2NHRcUVEBtVqNdevWYcSIEYLq\neP7557nYnJiYiD179gj+N2nbti1iYmK4/M0LCgrgQUj57p0Fgpr5JaJ8AEMCOdcT76j/JMY6AdhE\nRL3dx9VElODzeSURNWOM/dOln6a7398K4H0iumZHL8aY0+FwMJlM9F1CAsJmsyE/Px+dO3cW25SA\nsdlsiImJAY+/se8DKhFdtGnTRldUVDSQiM7W/owxNgjAB0Q00n1c45kOhEDvNuZuHi4yxu53K30I\nrtgJAGwEMJExFscY6wKgG1yzQ9cKZMxpt9sDtVN07HY7jh07JrYZQfHTTz/h0qVLgsiyWCy4evWq\n99jjUCwWC8rLywXRIRE2GAB/JfGOAujGGOvEGIsDMBGu5zpgGnQqjLHlAA4C6M4Yy2GMvQDgzwA+\nY4ylAvjYfQwiOgNgNYAzAH4B8Ar56QrFxMQ4eZT6y87O5nKTK5VKjB8/XnC5AGps1iUk48ePR/fu\n3QWRtWPHDhgMhmvet9ls15SFDJXLly+HvJFYQ2zfvp2L3NzcXFRXVwsul4i4yAUAu93O4Aq81qXX\nAeBvALYDOA3XxMs1PZp6qS/gwrPxSn5LS0sLKuEqEpgzZ47YJgQNjw3a9Xo9t9onvPaA2rdv3zVT\n7UJgMBi43RcajcYEjslvAcVUeKBQKAw5OTkqnhucC83FixfBY2U1TyoqKtC8efOQrrVYLDhz5gzu\nuOOOaz6rK6Zy+fJlJCQkBLwiWkIcNBqN2WAwdCYiLl1k0aKkcXFxVULlUISLI0fqDA9FNEuXLg15\nz+cLFy6gRYu6l4jUFaS98cYbce7cuZB0SYQHu90Ok8kUB4BbIExMp1JQWFgouFyLxYKMjAzB5QLA\nM888w0Wu3W4Hj78FALz22mshb/J+++23I5gtVLRaLe67776QdJ0+fRrr168P6dqGyMzMxNmzwYUF\nAuXUqVNc5BqNRi71ZUpKSqBQKHRExG2WRDSn4nQ683zn4YVCLpcjNzdXcLm82bp1q9gmAHA55Q0b\nNjR4Xn1FmgBg165dKC0tDVjvrbfeikceeSTg84NBJpNxKa4FgFvPbM+ePcjKyhJcbkFBAeRyOb+K\nYPCfqssdg8GQyaOyVUxMDEaPHi24XMA1s9SyZUuo1WpB5crlcrzwwguCyvSlvLwcFRUV6NatW4Pn\nms1m3HnnnY3WOWDAAFRU1F7d4R+ZTMYlnwYA19wiXjOCvO7hwsJCyOXyPC7C3YjWU7FarbkFBQV8\n6gdyori4mNv0L08SEhICLqeYmJgY0M57DSW+abXagIdOQuXSSDRMQUEBHA5HFk8dYqazFpSWlta/\nMClEeAVU77rrLnTt2pWL7Kqqqnq38GwMsbGxGD58uN/P7XY7vvzySy66AeD777/322txOp1ISbkm\n4VowXMWh+Awts7OzUVZWxkU2r4TCgoIC0ul0fAr2uhHTqRTm5eU5eQgOZiwfSfAq2O2L03ntn1wu\nl+MPf/hDUHIaiqn48uSTT0Kr1db5mUwm4xYA98jv3bs3F9n5+fnglZKxfPlyLnKPHTtmcTqdfL69\n3IjaU8nNzeVS++Dhhx/mIRYAkJqaykVuYmIiRo0axUW2h6qqKr89ksTExDrfFwK1Wg25/NrwXahT\n3cGQmJhYo2SDkAwePJhbTs5rr73GRa5Op7PAVbqEG6L2VKqrq+PESr4LlezsbLFNCJnExES8/vrr\n3uOPP/4Yoa6/CnUxoUdnRkZGje07JMKDOybI1amIllELAPHx8aaioiJFqBmf/iAiJCcnh5wzIRZE\nhJ07d2LYsHqrSAiCw+EAY4zbjIs/nE5n2HRu2rQJ3bt3R48ePQSXXVlZiezsbPTp00dw2WVlZdBq\ntYiPjxdcdvPmzY2VlZU9iIjbDJCodQeUSmUZj5wSxhisVi4xYK4wxqBUKrnqsFqtcDqd2LlzZ6MC\n2sHEVHyprKz0vjabzSHrD4QhQ4YItqCyNtXV1VCpVFxkb9u2DTwW21osFuh0uni4ak5zQ1SnwhjL\n2LVrFxfZDz30EBe5ALB7925usu+9915usgHghx9+QHl5OUaMGIFBgwZx1VWb4uJibNu2DYBrxumr\nr77iqk+r1XoKEglOx44dua0De+YZ1y6OQpORkQGNRpNHrjKR3BDVqVRWVu4+efJk1HUpeN2o4eCl\nl166Jrh46tSpoHstocRUWrVq5a1rK5fL8fbbbwctI1CibV1ZODh+/DjgqojPFbHLrh0/ceKEiYfg\nnJwcnD59modo7hXP9u3bh8OHDwsmz2Kx1Ju0d9ttt4HXxm7bt29HcnJyvecYjUZB8z1sNhu3dUQA\ncOXKFW6zgDqdjts6sEOHDpmqqqr2cxHug9hO5cS5c+dUPKYWk5KS6szJiAbuueceQQOAe/furXdx\nGmMMnTp18h7Pnj27RpW3uvAXUyGiGmuvhg8f3uCQzuFw4Ndff633nGCIjY3Fyy+/LJi82sTExHBL\n/b906RK3xLft27cTgONchPvCq1BLoE2tVpecOnWq0YVnws2aNWua7L7AdrvdWyHfbrfT999/f82e\nzJ4iTXa7nXz/fyUlJbRhw4aw2SoRGGazmeRyuQ0NVMIXoondU0FsbOxR91gvqhg4cGDIOR6Bcvbs\n2ZATxCwWS8izO77FsmUyGR588EFvHMlqtWLq1KneIaDD4aixdueGG27AuHHjQtILuFb9NiYjmlcm\narTjDtLmE5GRty7RnUpVVdWelJQULsHa8+fPcxv7dujQgduUoger1RryYrusrCy0bt260TYwxmos\nDIyLi8P/+3//r8ZxY5xIbdq1a4fLly+HdC0RCbLCuj7mzp3Lpc4J4Jpuv3DhAhfZx44dA9WxqwUP\nRE1+AwDG2IN9+vT5MS0tTfA8cYPBgLKyshrxAqEhoqieDQqV63WLjsaU52yIixcvwuFw4JZbbhFc\n9rBhwyw7d+78PyKaJbjwWojeUwFw4vTp02oewVq1Ws3VoVRWVmLWLO7/IxBRQEFni8WCNWvWcLcn\nXGzbti3godDFixcbPkkAeDkUALj55pu5OBQAyM7OtiAcQVpEgFMhokqFQlEWjbVNmzVrhr/85S/c\n9Vy8eBGrVq1q8DyHw4HBgwdztwfgP60OuBbsBRK3qqqq4lYusraeaMRisSArK0sF4GQ49InuVAAg\nJibm4P79fKbPDQYD5s+fz0U24Jq+5E337t0xceLEBs9TqVRo164dd3vChVarRZs2bRo8LzExEWPH\njuVqy9WrV/Hjjz9yk5+SksKtANiBAwegUqmywhGkBSLEqVRVVa1etWoVl52T1Go1JkyYwEO0l0uX\nLnGrq+HBX9zG6XTiv//9L1fddRHq2p9QmT9//jU9BYPBwH0GzkNSUhLXkp9xcXFo2bIlF9nz58+3\nG43Ghru6AhERTgXA1gMHDij1ej0X4TzHwYBrlzpeVdtq89lnn9WYfZDJZGEZgonN73//+2sKPa1e\nvZrbLn7h5o477gh514OGSE5ONttsNn4pxrUQffbHQ1JSUsrChQvveuyxx7jIt1gsXJaShxu73Q65\nXH7dzjoB4Z9x27ZtG0aMGBE2fUJy/vx59OvXr9JoNCZRmB72SOmpoLKycunatWu5jfnEGCLwwFNB\nbdq0adxLB0QiDocDb731VliXYPgrhSkUPGcQV65cSTKZbEO4HAoA8dP0PQ1AZ41GY7Lb7SEmIovP\nwoULyWazhUWX1WqlKVOmhEVXXfDYSzkQ9u/fT6dPnxZFNy9KS0u5ye7Ro4cewCgK47McMcMfANBq\ntbnbtm1rH65pUaEpKSlBUlJSnfVYhcBisSAmJoab/GCIlOQ3o9HILbPZZrOFZXaPF2VlZWjfvr3Z\nYrE0J6KwdWsjZvgDADabben69eu5FZDJzc0NaoOrYLnxxhu5PvArV66scwWrw+HAyZNhSUHwEk6H\nkpmZiX379l3zvtVqxbx587jpnT59OjfZgGvmzmTiUvkDAPDLL79ApVLtD6dDASIoUAsAjLG7Onbs\nuDM7O5vLIDY/Px95eXkYOHAgD/FeLl68yK0qWF0QEX766Sc8+uijYdMZTnJzc9G6deuo7jXUxdGj\nR6HT6fDggw9ykX/fffcZk5OTXyWiRVwU+CHSnIpMqVReTUtLS+RVWzQcbNy4Effff78g215YLBaU\nlpaiffv2AV8TjuLSvIc/RqPRUxojoPOrqqpgtVq5bZkRbVgsFiQkJNisVms7IgrrRlgRNfwhImdM\nTMzm77//PjqrK7kZO3asYPvoHDx4MOiZjlmzZnEd5oWDTZs2wWAwBHy+TCaDEFnZDocDS5YsabQc\nsdm5cydUKtX5cDsUAJEz++NpAAa3a9dO5ykSxINFixZROGaZahc2EoNIsCFQqqurxTaBrFYr5eTk\ncNezZ88ervKHDRumY4z9iUR4hiOqp+LmkE6nKxayvGBthg4dGpY8h927d9cZYGwIi8WCAwcOCGLD\nqVOnsG7dOkFk8USn0wlWYCk9PT3kQk+xsbEBbVDfGBwOh+cLlAvZ2dnYu3dvDBGt4KakPsTwZA01\nAH8eOHCgoXG+OjIIpadw+fJlys7O5mANUW5uriC9FyHyVM6cOUNFRUWNllOb6upqOnbsWFDXXL16\nlaKxrGldvP3221aVSvUNifT8RmJPBQCWp6WlMR4bjfmi0+m4ygdC286ja9euNaqtCUl6enrY1ik1\nREVFBRISEgSXq9Vqg64Ad+nSJW57LocTi8WCefPm2Y1G4//EsiEinQoR6eVy+eK5c+dy3fSIZ46D\nL0TUYM6DxWLBsmXLuNsyevRo70ySw+HARx99FFJXPJSZn+TkZOzcudN7PHjwYO47Mm7atCmgodCA\nAQO4LzwFgE8++YTr0GfdunVgjGUQ0XluShogoqaUfWGM3ZqYmHi8pKREGRcXJ7Y5jaahzE+LxYLy\n8nJRvy3Ly8uxdOlSTJ48GYBrapoxFlRvy2w2Q6FQAHAV7k5PT8dTTz3Fxd5A0Ol0sFgsfssKnDp1\nCr169QqbPbwXtnbt2tV05cqV3xMRv+IvDRCRPRUAIKKzRHQuHN/e4aChVPL4+HjRu98tWrTwOhQA\nKCoqwoIFC7zHOTk5WLhwIQBXnorvMeDaZGvDhg3e41tvvVVUhwK4hkL+HIrT6US4Kw7ydCjp6eko\nKioyA9jITUkgiBXMCaQBGH/nnXdyn2dcsmQJbxVe0tLS6OeffyYiVxB32rRpUTXt60GsBYWN4euv\nvxZt2vrgwYPc/89//OMfTXFxcR+SyM9txA5/AIAxFqtSqUpSUlKa8eyiXrp0Cd26deMmvzYmk8kb\nS/B9LcEXs9mM+Ph4XLhwAa1btxYsQTEQtm7dipEjR3KTX11djVatWlnMZnNXIirgpigAInb4AwBE\nZHM4HLM/+eQTrguiwulQAHhjDkQkOZQwolAowBjDlStXAk7/FwqeDgUAZs2a5YyLi9sttkMBItyp\nAIDFYpn5448/hrzBVDCEQwfgWv1qMpmwceNGpKWlhUWn0IS7Rq2QjBgxAlOnTuU6C+PBZuM6gQkA\n0Ov1+OyzzyzV1dX/4q4sACLeqRBRGWPss7fffpt7JfBDhw6FvM1oMLzzzjtQKpUYN26coBuxS/gn\nOTnZu/sfYwz/+c9/uJekdDgc+Pzzz7nqAIAvvvjCTkQ7iIjPdpxBEtExFQ+MMY1KpcpLTk5OvOOO\nO8Q2JyQ8xar9Rf8zMzPRuXNn7quLr1eysrLQqVOnaxwJEUGv13MvGcmLsrIydOjQwWI2m3sRUWh7\n5ApMVNzBRKS3WCz/fuGFF/hVtKmpT3CZa9eurbPAkofKykqcOXNGcL3XM06n07vGq3PnznX2TCwW\nCxYvXhxu0wRjypQplpiYmOWR4lAARPaUsm8DEKdWq4vDMZX5ww8/0OXLl7nriWaiYUp5wYIFlJeX\nF3a9hYWFtHr1au56srOzSalUGgG0pgh4Rj0tKoY/Hhhjk3r16jUvPT1dw3M87PmbNFaHxWJBfn4+\nunbtGtR1e/fuhUwmw3333dco/TyJlBq1QlJWVgYianShJ6PRFf7jVTvXw/jx400///zzVyaT6Z9c\nFQVJVAx/fFiZlZVVtHr1aq5Kgk1N98exY8dCqll7//33cy952Vgi1aHMnz8/5BmX+Ph4HD58uNE2\nqFQq7g7l9OnT2Lx5M5nN5mlcFYWC2F2lYBuAUTfccIMpHFthGAwGmjVrFnc99VFcXEyHDx8W1YZo\ngkcphUD5+eefyWAIT8WOESNG6GNjY9+iCHgma7do66kAwFar1Xp60aJF3MdtKpUKzz//fFDXWCwW\nQXM4WrZsGRFbctQmUvJUrFYr9u7d6z1u1aqVIHKPHz8edKGnjh07cu+hAK4So+6tTL/mriwEos6p\nEBFVVVW9+tZbb5mvXr3KXZ9Gownq/NLS0qBjKPUhk8lq1AZJTU31xnwkgKtXr3JZiNm9e3cUFhYG\ndU04Vjvb7XZMnDjRbDAY/o/CvPVGoESdUwEAIkqxWq0/PPbYY5aGzxaGdevWIS8vr8Hz2rdvz63A\nEuCqGl9VVcVNfqCIGVM5deoUCgpc2eitW7fmsh2KVqtF7969GzzPYDCEtdc2ffp0e0VFRSqAH8Km\nNFjEHn+F2gCoVSpV4YYNGwIfiDYCk8lEOp2uzs/MZjMtWrQoLHb4cuXKFUpNTQ27XrE5fPgwWa3W\nsOlbt24dlZSU1PlZVlZW2OI4GRkZpFKp9AA6UQQ8g/5aVE0p14Yxdl/z5s23Xrp0SZWUlCSaHXa7\nHeXl5YKN54PRe+bMmYC+UYUmnFPK+fn5SE5OFq02i06ng9PpDOuq5trY7Xb07NnTeOXKlX/YbLbw\nlCwMkagc/nggov0Wi2Xxs88+G9ax5cyZM2vENeRyedgdikevr0NZtWpVk8nK9f37Nm/eHI899pho\ntmi12mscim8xqnDw6aef2gsLCzPsdvv8sCoOBbG7So1tANRKpbIoXMMgIlfldafTSVOnTo24Aku+\n9qSmpkacfYHgdDrpo48+CsveTMHyxRdfkE6no7S0tLDpjJZhj6dF9fDHg1jDIKvVitjYWO6rXUNl\n8+bNGDFiBGJjY+F0OkFEiImJEdusOlm2bBkGDx6MLl26iG1KvVitVoSzZrLdbkefPn0MFy5ceDPS\nhz0eonr444GI9lut1sUvv/wy9wWHvpuQMcYwdepU3ipDZsyYMd5NzfV6PWbNmuX9zONkQqWxMx6H\nDh3CsWPHvMcTJkyIeIcyc+bMGk45HBvSvffee/acnJzoGPZ4ELurJFSDaxhUvGrVqob7k41gxowZ\npNfrvcfROLwgcs0c+c5YlZeXU35+fsDXB7ugcP/+/bRz507vcVlZGfHc2pYHVVVV3tdOp5M+/PBD\nrv9/n2FPR4qAZyzQ1iSGPx4YY/clJiZuS09PV/LMFakLh8OBiooKv5XbI538/HxkZWXhnnvuAeBK\nssvOzsajjz4KwFXHt6qqypuIl5OTg6qqKtx+++0AgJMnTyI3NxdjxowBAKSkpKCsrAwPP/wwANe3\nejTWirHb7aJkNOv1evTr18/gnu2Jnl4K0HR6Kp4WHx//do8ePfS+vYnGYjabG1zTYTQa6bvvvhNM\nZ6RRWVlJubm53uPi4mLKzMz0Hlut1qjreQTC9OnTyWQy1XuOw+GgiooKwXQ6HA4aMGCAWa1Wr4C7\nkFo0NdENEPwXAphGo1k9ePBgs1Bd09WrV4tSlyOSiYZ6KuHCYDDQvHnzBJP373//26rRaE4CiKcI\neKaCbaIbwOWXAhQajSbjww8/DF/apQ8Wi4UWLFgghuqw0ZSdyo4dO0Sbzl67di2pVKpSRFjhpWBa\nk4qp+MIYa6tSqdKXL1/eYty4cUFfb7FYkJWVhR49eoSkv7CwEG3atAnpWglxOXjwIAYPHhzStYWF\nhZDL5SEVejpx4gTuueces9lsvpeIjodkQAQQfZGzACGiAqPROGrixImWkydPBn39yZMnG7Unj69D\nCb1d+m4AAA51SURBVHa1q0T4sVqt3tehOhTAtar9xIkTQV9XWlqKUaNGGc1m8wvR7FCAJuxUAICI\njlqt1j+NHDnSWF/R6bq46667BFttvGXLlho3bVMgUuqpCIHBYMCcOXMEkaXVajFixIigrrFarRgz\nZoxBr9fPIaKVghgiIk3aqQCAw+FYWl1dPf+RRx4xNFRm0GKxYPv27YLb8OKLL4Y1C1MiONRqNV5/\n/XXB5R4+fLjBQk9EhHHjxlnPnDlz2Gg0viO4ESLQ5J0KABiNxv/LyMg4NmbMGFt9MaTKykrccsst\nXG1Zs2YNTp06xVVHOIjUGrWBcuLECfz8889cddx22231bssCAF9//bVz//79hXq9/nEi4p+iGwaa\nbKC2NoyxRLVafeSVV17pPH369LhIXa8jER5sNpt3CYNYLFq0iF599dUKo9E4gIgyRTVGQK6LngoA\nEFGVwWC4Z+7cuXkffPCBdxxktVrx7bffimJTQUGBaLobSzTGVL7//ntUVlYCQNgdyqpVq1BWVuY9\nnjVrFv31r3+tNhqN9zQlhwKgaeap1NcAtFKpVHmTJ0+2EbmyF/1V9QoHkbi8PxCiMU/l6tWroumu\nrq72ZmX/+OOPpFQqqwDcThHwTAjdRDdAlF8aaKdSqQpmzJgRUU/0vHnzpMxdAcnKyqJly5aJbUYN\nNm3aRCqVSgfgToqAZ4FHu25iKrVhjHVSq9WH33vvvZb//Oc/I2YPDCKK2Pos0YbFYkFsbGzELGRc\nsmQJ/vznPxvMZvNDRJQitj28iIy/tggQUbbBYOg/ZcqUktdff73eWaFw4nEohYWF+Oabb0S2xj+R\nGlOZMWMGdDodANeOg5HiUFasWEEvv/xytdls/l1TdigArs/hj28D0FqtVme+/fbblkiujXL+/Pmw\n7X4XCJESU7HZbFRaWiq2GfWyYMECp0qlqkATjaHUbpHhxkWEiIoMBsNdc+bMyXrxxRctdrtdbJPq\nRCaT4dKlS2Kb4SVS8lSSk5NRUlIithl1QkSYMWOG429/+1uV0Wi8m4gyxLYpLIjt1YJtANoD2AXg\nNIAMAK+53/8MwFkAaQDWAUhwv98JgBHACXeb4yNrDICTAOYDaKbRaA4NGjTIIOYsQaB88cUXos5m\niEVubi7Nnz9fbDMaxGKx0B/+8AeTRqO5DKAHgBQAqe779hNy3X/jAZwC4ADQj367LwO6Z4njc9aY\nJroBQRsMtAbQ1/1aA+A8gFsADAUgc7//KYBpPv+gdD+yVsIVV/oIQE8AcrVaPat9+/aGM2fOBHDr\nRAZGo5HWr18fVp3hGv44nU7avXu3t2yj0+mM+BKexcXFdOuttxq1Wu02AFpy3Wsq988YAIcB3ON2\nNje7vyRrO5WA7tm6zhG7Rd3wh4iKiCjN/VoPV++kHRHtpN/SnA/D1aPx4G86hQGIA6ACYCMiu16v\n/1thYeFrd955p3ndunWcfgthUSgUNZYXlJeX4/z58yJa1DhycnJgNv+2lZNvsJUxFtGzY2lpaejV\nq5cxKyvrK51ON4qIdABAREb3KfFwOYUKIjpPRBdR9/0Z0D0rsPnCILZXa0wD0BlAFgBNrfc3Apjk\n4/V1cHUjdwO41+e8oQCOAZheh+y7lUrl1WnTptki/ZuxNtXV1ZScnOw9zs3NpcuXL4toUf1kZ2fX\nGMr9+OOPJGQ50HCxatUqp1qtNjDGJtC195MMruFPNYDPan22G9f2VIK+ZyOliW5AyIa7hj7HAIyr\n9f7/A7DO5zgWQHP3634Acmo7oXp0tNdqtWefeOIJo9ForOs+igpKSkpqbH6VkpJCx44dE82eo0eP\n0vnz573Hu3btivgZnPpwOBz0+uuvWxUKRTmAO6j+eyoBrp70/T7v1XYqId+zkdBENyAkowE5gK0A\nJtd6/w8ADqCe2p61/4EB6FJpNJqN3bp1MzSVbFeLxUJlZWXe419//ZX27NnjPT5x4kSNh95gMFyz\nIbpvTMVms5HZbPYep6en17h+06ZNdPToUe9xVlZWRE2PNwadTkejRo0yaDSaVACtKLB76t8A3qQA\n78lg71mxm+gGhGQ0sBjAF7XeGwlXZL1Frfdb4rcAblcAuQCaBamPKRSK95KSkgzbt2+npk5JSQkV\nFxd7j/ft20cpKSne461bt9LXX3/tPd6xY0eNnk9OTk4Np9VUSUlJoQ4dOhg0Gs1SAHHk//5pCSDR\n/VoJYB+Ah3w+3w2ftH0h7lkxm+gGBG2wK2rugGvqONU97hwF4CKAbNSahgPwOFzTdifcw6XRjdA9\nTKVSlY0ePdriu7GUxPWF1Wql999/36pQKPQymexFNLCNBoDb3fdfKlzTwf/nfv9Rt8MwASgEsIUE\nvmfFaNft2p9QcddlmaVSqR5ftmyZatiwYWKbJBFG0tLSMHHiRENhYeHx6urq3xNRrtg2RRpRN6Us\nNkRUpdfrny0tLX3s0UcfLRs7dqylurpabLPCTqSu/eGFzWbDBx98YBs0aJDx4sWLk6urq4dIDqVu\nJKcSIkS03Wg0dtu1a9e6bt26GXfs2CG2SRKcSE9PR+/evQ3/+9//DlksllscDsdCkrr4fpGGPwLA\nGBuuVquXjRkzRjN//nxFQkKC2CZJCIDNZsOUKVNs06dPt9vt9tecTud3kjNpGMmpCIQ71jJHoVA8\nOnfuXNX48eMjOvNTon7279+Pl156SYqdhILYkeKm1gA8qNVqz99yyy26zZs3U1MlUkofCM3Jkydp\nwIABRqVSWQrgGUThBuliNymmIjBEtEun09167ty5Pz355JPF9913n/748ajecO664MqVK3jyySeN\ngwYNqk5NTf2XyWRqT0TLiEjqygeJNPzhCGMsViaTvaRQKKbefffd8XPmzFF2795dbLMkfCguLsaH\nH35oWbhwITHGZlgsls/IvQhQIjQkpxIGGGPquLi4N2Uy2duTJk2KmTJliqJt27Zim3VdU11djU8/\n/dT+5Zdf2hljiwwGw/tEFJnVnqINscdf11MD0EKtVn+pUCiMjz/+uDUnJ4eilWiNqZSXl9P06dMd\narXapNVq1wHoTBFwbzSlJsVUwggRlev1+slms7nHtm3bFnXv3t04atQo/Y4dO+B0NokdLyOWI0eO\nYPTo0eZ27dqZp06d+pPBYBhYXV39BBFliW1bU0Ma/ogIY0zDGJukVqvfValUN77zzjuKF154Qda8\neXOxTWsSmEwmrFq1Cp9//rkuKyvLbLVav7Tb7fOJqP5d0yUaheRUIgDmSmi5OyEh4f+sVuuoBx98\nEO+//77irrvuEtu0qOTSpUuYPXu29ZtvvkF8fHxKVVXVZ3At1nOIbdv1gORUIgzG2I1yufyluLi4\nyV26dFG8+eab2nHjxiEpKUls02qwZ8+eiKmoDwAGgwErVqzAt99+a8jIyCAAC0wm0ywiuiy2bdcb\nklOJUBhjMQBGJSYmTjYajff16dPH8tRTT2nHjRvHbr75ZrHNiwinUlBQgI0bN2L16tW6gwcPxiuV\nyozKysqvAKwhIpOoxl3HSE4lCmCMKeHK1H3SbreP1Wg0cc8991zsY489Fjto0CDExMSIbWJYICJk\nZGRgw4YNjpUrVxoyMzPlCoVid1VV1VIAW4moUmwbJSSnEnUwxmQA7oyPj39MoVA8Zbfb2w0ePJj+\n8pe/KIYMGRJxw6TGotfrcfDgQaxfv96yevVqu9lsNjPGVhsMhrUA9hNRZFaUv46RnEqUwxjrBOCR\n5s2bP2MwGO5Qq9XOXr16OUeOHKkaMGAA69evH1q0aCG4Xh7DH71ej9TUVBw/fhwHDhwwHD582FlY\nWKhMSEg4r9PpVtjt9g0AzpB000Y0klNpQrh7MTcD6K9UKgcplcr7dDrdrQkJCfaBAwc67rnnHs0d\nd9zBBgwYgBYtWjRqFXVjnYpOp0NaWhq2bNmCc+fOGY4ePeosKipSqFSqXIfDscdgMBwEcBzAaak3\nEl1ITqWJ4+No7lQqlXfHx8c/YDAYbmKMyVu0aGFSqVSym266iXr06BHXvn37uBtvvBEdOnRA27Zt\n0aZNGyQmJgblfPR6PQoLC1FQUIDCwkKcPXsWlZWV9vz8fHNOTo4jMzMz1mAwyO12u0yr1V62Wq2H\nDAZDMiQH0mSQnMp1CmNMC6AtgDaenwqFoqNCobhZJpO1sdvtrcxmc5LD4ZDL5XKHXC53xsTEkNPp\nZPHx8faYmBjY7XZmsVhi3OnZMrvdLnM6nVAqlVfj4uJKGGN5RqMxx2QyZcJV2LnA52elNIxpmkhO\nRaJeGGMquLbZlMO1yZXc3WRwbbtpA2B3NxsAg+Qsrm8kpyIhISEo0oJCCQkJQZGcioSEhKBITkVC\nQkJQJKciISEhKJJTkZCQEBTJqUh4YYy1Z4ztYoydZoxlMMZec7+/kjF2wt2uMMZO+FzzLmPsImPs\nLGNsuM/7Yxhj/7+9O1aJKwijOP4/KBbBN1gDCcRCK5V0IWVqIcVq0CewESytfAAbW4U0gSSENFoF\nsmWaCCEGYS0EO8UyhaVyUtxZc7NJkeK6u+D5VbMfwzK3Ocy3zJ39IWl3GM8SwzM+7AXESLkGNmwf\nSZoEvkn6bHu5N0HSNvCzjGeANjADTAEdSdPlnMoqMA9sSZq13R30w8RwZKcSt2xf2j4q4yvgBGj1\nTWsDb8t4EXhv+7rc9XoK9K6rE9WhuQdUh+LinkioxD9JegTMAV9rtefApe2zUmoB9b8DPed3CO0B\nX4Ab26d3vd4YHWl/4i+l9fkIrJcdS88r4N3/fIftDvD0DpYXIy6hEn+QNE4VKG9s79fqY8BLYKE2\n/Rx4WPs8VWpxj6X9iX6vqS5C2umrvwBObF/UagfAsqQJSY+BJ8DhgNYZIyo7lbgl6RmwAhxL+g4Y\n2LT9CViir/Wx3ZX0AehS/Ri7ljeUI28pR0Sj0v5ERKMSKhHRqIRKRDQqoRIRjUqoRESjEioR0aiE\nSkQ06hdaPUpRSxQqGwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f90048e5710>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Enable polar coordinates\n",
    "ax = plt.axes(polar=True)\n",
    "# Use the arrow function\n",
    "# The arrow forms a 45° angle with the x axis and has length 1\n",
    "ax.arrow(x=45/180.*np.pi, y=0.0, dx=0, dy=1)\n",
    "# We can still use the Cartesian format\n",
    "ax.arrow(x=0.0, y=0.0, dx=0.5, dy=0.5)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Exercise 5\n",
    "\n",
    "The Taylor expansion of $\\mathrm{e}^{\\mathrm{i}\\theta}$ is given below. In this exercise we will try to ``prove'', by a convincing yet non-formal argument, that Euler's formula does indeed hold true.\n",
    "\n",
    "\\begin{equation}\n",
    "\\begin{aligned}\n",
    " e^{\\mathrm{i}\\theta} &= \\sum_{n = 0}^\\infty \\frac{(\\mathrm{i}\\theta)^n}{n!} \\\\\n",
    " &= 1 + \\mathrm{i}\\theta + \\frac{(\\mathrm{i}\\theta)^2}{2!} + \\frac{(\\mathrm{i}\\theta)^3}{3!} + \\frac{(\\mathrm{i}\\theta)^3}{3!} + \\ldots \\\\\n",
    " &= \\left( 1 - \\frac{\\theta^2}{2!} + \\frac{\\theta^4}{4!} - \\frac{\\theta^6}{6!} +\\ldots\\right)\n",
    " + \\mathrm{i}\\left( \\theta - \\frac{\\theta^3}{3!} + \\frac{\\theta^5}{5!} - \\frac{\\theta^7}{7!} + \\ldots \\right) \\\\\n",
    " &= \\cos\\theta + \\mathrm{i}\\sin\\theta\n",
    " \\end{aligned}\n",
    "\\end{equation}\n",
    "\n",
    "Using `matplotlib`, you will produce two plots. \n",
    "The former will contain the plot of $\\cos\\theta$, $1$, $\\left( 1 - \\frac{\\theta^2}{2!}\\right)$,\n",
    "$\\left( 1 - \\frac{\\theta^2}{2!} + \\frac{\\theta^4}{4!}\\right)$ and $\\left( 1 - \\frac{\\theta^2}{2!} + \\frac{\\theta^4}{4!} - \\frac{\\theta^6}{6!} \\right)$.\n",
    "The latter will contain the plot of $\\sin\\theta$, $\\theta$, $\\left(\\theta - \\frac{\\theta^3}{3!}\\right)$,\n",
    "$\\left(\\theta - \\frac{\\theta^3}{3!} + \\frac{\\theta^5}{5!}\\right)$ and $\\left(\\theta - \\frac{\\theta^3}{3!} + \\frac{\\theta^5}{5!} - \\frac{\\theta^7}{7!} \\right)$.\n",
    "What is it you notice?\n",
    "\n",
    "As a quick reminder, you should use **both** `numpy` and `matplotlib` to generate these plots. Defining functions will also help ;)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Bonus\n",
    "\n",
    "The general expression for a Taylor expansion at $x_0$ is:\n",
    "\n",
    "\\begin{equation}\n",
    " f(x) = \\sum_{n=0}^\\infty\\frac{(x-x_0)^n}{n!}f^{(n)}(x_0)\n",
    "\\end{equation}\n",
    "\n",
    "where $f^{(n)}(x_0) = \\frac{\\mathrm{d}^n f}{\\mathrm{d}x^n}(x_0)$ is the $n$-th order derivative function of $f(x)$ evaluated at $x_0$.\n",
    "Write down explicitly the first few terms for the expansion of $\\cos(x)$ and $\\sin(x)$ for $x_0=0$ and verify that the identification made in the Taylor expansion above is correct."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}