{ "metadata": { "kernelspec": { "codemirror_mode": { "name": "ipython", "version": 2 }, "display_name": "IPython (Python 2)", "language": "python", "name": "python2" }, "name": "", "signature": "sha256:edca853a48a4921f446560c028cc5b5d2722e33c45a91bbb29b9c59e5bed1529" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "The purpose of this assignment is to show that one can obtain the Legendre Polynomials numerically from the QR decompostion" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The Legendre polynomials are orthogonal on [-1,1] and have many applications. Please see \n", "http://en.wikipedia.org/wiki/Legendre_polynomials" ] }, { "cell_type": "code", "collapsed": false, "input": [ "# Let x be a vector of equally spaced points from -1 to 1\n", "# Let the columns of A be the powers of x\n", "# One can think of A as having first column 1, second column x, third column x^2, etc.\n", "N = 1000\n", "h = 1/N\n", "k = 5\n", "A = zeros(2*N+1,k+1)\n", "x = -1: 1/N : 1\n", "for i=1:2N+1,j=1:k\n", " A[i,j] = x[i]^(j-1)\n", "end" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "# Compute the QR decomposition of A (Gram-Schmidt)\n", "(q,r) = qr(A)\n", "q=q/sqrt(h); # Normalize q to integrate to 1;" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "You are asked to compare graphically with the first five columns of Q with the Legendre polynomials to see if they line up\n", "Note: You must compare with the nth Legendre polynomial times sqrt(n+.5). Do you see why? Note that it is possible\n", "you may be off by a factor of -1. Can you explain why?" ] }, { "cell_type": "code", "collapsed": false, "input": [ "using PyPlot;" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 23 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(x, -q[:,5], linewidth=2,\"r\")\n", "plot(x, (35*x.^4 - 30*x.^2 + 3)/8 * sqrt(4.5))\n", "legend({\"QR with h=$h\",\"Quartic Legendre Polynomial\"})" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAArYAAAIQCAYAAAB0Ri0fAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XlclOX+//HXPeygQqigsiOQhnuLllqYqZVt5smyr1lpq9W3+p1zbNMSy/J0bPmWeeqYiZ3KFo9aWWppmOaSqZWmogjDJigKKIoL2/37gyVHQBEZhuX9fDzmcU733Pd1f2YYxjfXXPO5DdM0TUREREREmjiLowsQEREREakPCrYiIiIi0iwo2IqIiIhIs6BgKyIiIiLNgoKtiIiIiDQLCrYiIiIi0iwo2IqIiIhIs6BgKyIiIiLNgoKtiIiIiDQLdg+2K1eu5O677yYqKgovLy8CAwO55ZZb2LJly1mPjYuLw2KxVHvLzs62d+kiIiIi0oQ42/sE7733HgcOHODJJ58kOjqaAwcO8Nprr9GvXz+WL1/OoEGDzjpGXFwcXbp0sdnm6+trr5JFREREpAkyTNM07XmC7Oxs/Pz8bLYVFBQQERFBt27d+P7772s8Ni4ujnHjxrFp0yb69OljzzJFREREpImz+1KE00MtgJeXF127diUjI6NWY9g5e4uIiIhIM+CQL48dPnyYLVu2EB0dXav9b7jhBpydnWnbti0jR45k+/btdq5QRERERJoauy9FqM6YMWP44osv2LBhA717965xv+XLl7N27Vr69etHmzZt2Lp1K9OnTycvL49169bRvXv3BqxaRERERBqzBg+2kydPZtq0acycOZMJEyac8/Gpqal0796dwYMHs2jRomr3ycrKIisr63xLFRERERE76dixIx07dqzfQc0GNGXKFNMwDPOVV145r3GuvfZas0OHDtXel5mZaXbq1MkEdNNNN91000033XRrpLdOnTqZmZmZ55UJT2f3dl8VYmNjK29PP/30eY9nGEa127OyssjMzOSjjz6ia9eu530ekTN54oknePPNNx1dhrQAeq1JQ9FrTRrCzp07GTNmDFlZWfU6a9sgwfbFF18kNjaWyZMnM3ny5PMaKzk5mTVr1jBs2LAz7te1a1e1CBO78/Hx0etMGoRea9JQ9FqTpszuwfa1117jhRde4Nprr+X6669nw4YNNvf369cPgPHjx/Phhx+SnJxMUFAQAEOGDOHqq68mOjqaVq1asW3bNl599VWcnZ158cUX7V26iIiIiDQhdg+2S5YswTAMli1bxrJly2zuMwyDkpISAEpLSyktLbXpWdu9e3c+/vhj0tPTOX78OH5+flxzzTVMnjyZiIgIe5cuIiIiIk2I3YNtfHx8rfabO3cuc+fOtdn2+uuv26MkEREREWmGHHKBBpHmYvTo0Y4uQVoIvdakoei1Jk2Zgq3IedA/ANJQ9FqThqLXmjRlCrYiIiIi0iw0WB9bEZH6lpiYyJEjRxxdhoiInKJ169ZERkY65NwKtiLSJCUmJhIVFeXoMkREpBq7d+92SLhVsBWRJqliplZXGRQRaTwqrijmqE/TFGxFpEnTVQZFRKSCvjwmIiIiIs2Cgq2IiIiINAsKtiIiIiLSLCjYioiIiEizoGArIiIiIs2Cgq2IiIiINAsKtiIiIiLSLCjYiohIoxITE4PFcm7/PMXFxWGxWJg3b955n3/VqlVYLBZiY2PPeywRaVgKtiIizcSmTZu49957CQ8Px9PTE29vb3r06MHEiRPJzMys9piKEHfqzdXVlYCAAEaOHMmaNWsa+FGAYRgYhlFtnWcLm6cfd751NLSSkhLeeOMNevTogaenJ23btmX48OGsX7++QcbKzc3liSeeIDQ0FHd3dwICAhg/fjx79+6tdv85c+bw4IMP0rdvXzw9PbFYLEyePPmcaxWpL7rymIhIM/DUU0/xz3/+ExcXF4YMGcLtt99OYWEha9euZcaMGcyaNYt58+YxcuTIao8PDQ3lnnvuAeDYsWNs2rSJRYsWsXjxYj799FNuu+22BnssH374IcePH6/2PkeEzYZimiZ33HEH//3vf+nSpQuPPfYYOTk5fPbZZ1x55ZX897//5aabbrLbWDk5OVxxxRUkJiYyePBg7rzzTnbu3MncuXP55ptvWL9+PWFhYTbH/PWvfyU/Px9fX18CAgJISkpq1j8jaQLMZmbz5s0mYG7evNnRpYiIHel3/U9TpkwxDcMww8PDzR07dlS5/7///a/p4eFhOjs7m/Hx8Tb3xcfHm4ZhmIMGDapy3PTp003DMMzQ0FB7lV5rFXVOmTKl2vvnzp1rGoZhzps3r97OFRsbe95jnYtPPvnENAzDHDBggHny5MnK7b/88ovp5uZm+vn5mfn5+XYb64EHHjANwzD/9re/2Wx/6623TMMwzGuvvbbKeZYvX26mpaWZpmmacXFxpmEY5uTJk2v9mKX5qe17s73ew7UUQUSkCbNarbz00ku4urry1Vdf0bVr1yr73HrrrbzxxhuUlJTw8MMPY5pmrca+7777AEhLSyMnJ+es+x89ehRXV1cGDBhgs72goABXV1csFgsfffSRzX3/+te/sFgsxMXFVW47fY3tPffcw9VXXw1AbGyszbKJ1atX24xnmibx8fHExMTQpk0bvL29GT58ODt27KjVYz7db7/9xvDhw/Hx8cHLy4urrrqKdevW1Wmss5k1axYA06ZNw9XVtXL7JZdcwu23386BAwdYsGCBXcY6evQo//nPf2jdujVTpkyxGevRRx8lJCSE5cuXY7Vabe4bOnQoQUFBALV+XYnYk4KtiLQsl1wCgYENf7vkErs8nLlz51JSUsItt9xCdHR0jfvdd999dOjQgV27dvHjjz/WauySkpLK/+/m5nbW/Vu1akXfvn3ZtGkTBQUFldvXrFlDcXExACtXrrQ5ZuXKlRiGweDBg222n/px9ogRI7j77ruBstA7ZcqUyltoaKjNcUuWLGHYsGH4+Pjw8MMPM3DgQJYuXUpMTAwHDx6s1eOu8Msvv9C/f38KCwt54IEHuOGGG1i7di2DBw8mISHhnMY6mxMnTrB+/Xq8vLwYOHBglfuvu+46oGytsT3G2rBhAydOnKB///54eXnZ7G8YBsOGDQMgPj6+tg9JxCG0xlZEWpZ9+6CGL8I0RT/99BMAQ4YMOeN+Tk5ODBo0iPnz57N27VpiYmLOOvacOXMA6N27N61atapVPVdffTVr165l9erVlQFq5cqVODs7M3DgQJtgW1paSnx8POHh4ZWzftW5+eab8fb2Zt68ecTExPD888/XuO+XX37J8uXLGTRoUOW2Z599lunTp/PBBx8wceLEWj0OgG+++YZ58+Zx1113VW7797//zUMPPcRbb71VOSsKkJKSYjPrXBv33nsvISEhACQlJVFaWkp4eHi1a1QjIiIA2L1791nHrctYu3btAiAyMrLaMSuOSUxMPOv5RRxJwVZEWpYOHZrVefft2wdwxmBYITAw0OaYU1mt1sqPoI8dO8bmzZuJj4/Hz8+P999/v9b1DB48mBdffJGVK1faBNu+fftyyy238Pjjj5OYmEhkZCS//fYbeXl59frFtNGjR9uEWoAHHniA6dOns2nTpnMaa+DAgTahFmDcuHE88sgjVcZKSUlh6tSptR7bMAyuvvrqymB7+PBhALy9vavdv2L7oUOHzjp2Xcaqz/OLOJKCrYi0LOcYbpqK2nwTvWKfU5cYVEhNTa0SzNq1a0d8fHy163Zrcvnll+Ph4cEPP/wAlLWP2rp1K5MmTapcbvDDDz8QGRlZuU/F+tn6cEk1Sz4qAn1eXt55j+Xs7Iy/v3+VsWJiYigtLT2n8UWk/mmNrYhIE9ahfCY4LS3trPtmZGQAEBAQUOW+imBWWlpKbm4u7777Lvn5+QwbNozs7Oxa1+Pi4kL//v35/fffycnJIT4+ntLSUgYPHsxFF11Ehw4dKpcjrFy5EovFUq/BtroZR2fnsjmc6gL9uY5VMd65jlXbc1XMnJ6uYruPj49dxqrP84s4kmZsRUSasIEDB7Jq1SpWrFjB+PHja9yvpKSk8os//fv3P+OYPj4+PPDAAxQXF/Poo48yYcKEWn8bH8qWI6xYsYL4+Hh++OEHvLy8uPzyy4Gy2dlly5ZRWFjImjVriI6Opl27drUeu7E63zW2nTt3xmKxkJycTElJCU5OTjb7VqxtjYqKOuu4dRmrS5cuQM1reM/l/CKOpGArItKE3Xvvvbz88sssXLiQnTt31rhs4IMPPiArK4uoqCiuuuqqWo390EMP8a9//YuFCxeybt06rrjiilodd+qSg/j4eAYMGFA5azp48GA++eQTZs2axbFjx6p0Q6hJRTir75nS+nK+a2zd3d3p378/a9asYc2aNVW+3Ld06VKgdss26jJWv379cHd3Z+3atRw9etTmy4KlpaV89913GIZRZf2ySGOjpQgiIk1YaGgozz77LEVFRdx4443s3Lmzyj6LFy/m8ccfx8PDg3nz5tX6ylCnXsL2ueeeq3VNvXv3xtvbm4ULF7Jr1y6b8Frx/6dPnw7Ufn1t27ZtgdotuXCEU5dy1OZWUlLClVdeaTPGww8/DMCkSZM4efJk5fZffvmFzz77DD8/vypXjktLSyMhIaHKldrOdSwvLy/Gjh3L0aNHq/SxnTlzJqmpqQwbNqxKezWRxkYztiIiTdyUKVMoKCjg9ddfp2fPngwbNoyLLrqIoqIi1q1bx8aNG/Hz82POnDn07dv3nMa+9dZb6dWrFz/++CPfffcdQ4cOPesxTk5OxMTE8OWXXwLYBNvg4GA6d+5MUlISzs7ONc4en97sv0uXLgQGBvLpp5/i4uJCUFAQhmEwduxYgoODz+kxNVZ33HEHCxcuZMGCBfTu3Zsbbrih8jK4pmkye/bsKm3Xxo4dy+rVq4mPj7d5Lusy1ssvv8yqVat4/fXX+e2337j00kvZuXMnX331Ff7+/rzzzjtVan7//fcrW87t2bMHgK+++qryD5CuXbvy1FNP1evzJHImCrYiIk2cYRjMmDGDUaNGMWvWLH788UdWrFhROVMXGhrKr7/+Wu2XoWozezt16lRuuukmJk+eXKtgC2Vh9ssvv8TX15fevXtXuS8pKYmLL76Y1q1bV1vT6XVZLBYWLlzI008/zeeff86RI0cAuPLKKyuDbW1nos+XPc8zf/58rrjiCj744ANmzpyJh4cHMTExTJo0iX79+lVbS3XPV13G8vX1Zf369cTGxrJ48WLWrFlDu3btGDduHFOnTqVTp05Vjlm7di0ffvhh5fkNw2Dbtm1s3boVwzC46qqrFGylQRlmM7sG3pYtW7j44ovZvHkzffr0cXQ5ImIn+l0/u6NHjzJgwAC2bt3K1KlTmTRpkqNLEpFmrrbvzfZ6D9caWxGRZqpVq1YsWbKEjh078vzzz/PJJ584uiQREbvSUgQRkWYsMDCQpUuXsmjRItLT0yksLMTV1dXRZYmI2IWCrYhIM9ejRw969Ojh6DJEROxOSxFEREREpFlQsBURERGRZkHBVkRERESaBQVbEREREWkWFGxFREREpFlQsBURERGRZkHBVkRERESaBQVbEREREWkWFGxFREREpFlQsBURkUbPYrEwaNAgR5fRIoSGhhIWFuboMuxi1apVWCwWYmNjHV1Kg6qvxx0XF4fFYmHevHn1VFn9U7AVEWkmNm3axL333kt4eDienp54e3vTo0cPJk6cSGZmpqPLO6PahCnDMOxybovFgsWifw5PZa/nujbuueeeyp9Jxa1Vq1Z0796d5557jsOHD5/3ORz5+Bypvh53Y37+nB1dgIiInL+nnnqKf/7zn7i4uDBkyBBuv/12CgsLWbt2LTNmzGDWrFnMmzePkSNHOrrUGp3pH8uEhAQ8PT0dcm5xjFtuuYVevXoBkJWVxddff80rr7zCggUL2LhxI97e3g6usOno27cvCQkJtGvXztGl2J2CrYhIExcbG8s///lPwsLCWLJkCV27drW5f+HChYwZM4Y77riD77//npiYGMcUeh6ioqIcXYI0sFtuuYWxY8dW/veMGTPo27cvO3bs4O2332bSpEkOrK5p8fDwaHS/Q0lJ9hm3+X72sn69oysQEbE7q9XKSy+9hKurK1999VWVUAtw66238sYbb1BSUsLDDz+MaZqV902ZMgWLxcLq1aurHJeSkoLFYuHee++12b57926efvppLrnkEtq3b4+7uzuhoaE88MADpKenVxnn1PV9GzZs4Nprr+WCCy6oXKtnsVhIS0urPF/F7dTz1rTGtqSkhHfffZf+/fvj7e2Np6cnkZGR3H///ezZs+ecnsvamj9/PoMGDcLHxwcPDw8uuugipk2bRmFhYbX7f/zxx/Tp0wdPT0/8/f0ZO3YsmZmZxMTE1LgEYvny5Vx//fW0a9cOd3d3IiIimDhxYrUfw1cs4zh27Bh///vfCQ4Oxt3dncjISP7xj3/U+DhmzpxJdHQ0Hh4eBAYG8uijj9b4Mf+payu/+eYbrrzyStq0aWNTf3FxMbNmzaJfv360adMGLy8v+vTpwzvvvGPzmqsrLy8v7r77bgA2btxYuf3QoUM8/fTTREVF4eHhga+vL8OGDWPFihW1Gvfyyy/HycmJ1NTUau+fMWMGFouF119/vXJbXZ/zzz77jIEDB1a+Vrt3784rr7zCyZMnq+xbcY6CggKefPJJgoKC8PT0pFevXixcuBAoe85feuklIiMj8fDwICIignfeeafKWDWtsd28eTOPP/44PXv2pG3btpUB+K9//St5eXm1ev7qJD+fb57fYJehm+2MbfHGzfCIo6sQEbGvuXPnUlJSwsiRI4mOjq5xv/vuu4/Y2Fh27drFjz/+eE6ztqd/TL9w4ULee+89rr76agYMGICrqytbt25lzpw5fPXVV2zevJmAgIAq46xbt46XX36ZK6+8kvvvv5/s7GyioqJ44YUXePPNNwF48sknK/ev+Bi6pjoKCwu54YYbWLFiBSEhIdx11120bt0aq9XK4sWLGThwIBEREbV+nLUxbtw44uLiCA4OZtSoUXh7e7N+/XomT57MypUr+e6773B2/vOf1ldffZWnn34aX19f7rnnHry9vfnuu+8YMGAA3t7e1S6BiI2NJTY2lnbt2nHjjTfi5+fH77//zowZM/j2229Zt24dbdq0sXleioqKGDp0KFlZWQwfPhxnZ2cWLVrEM888w/Hjx5kyZYrNOR5//HHefvttOnXqxIMPPoizszNffvklGzdupKioCDc3t2of/xdffMGyZcsYPnw4EyZMwGq1AlBUVMSNN97Id999R9euXRkzZgzu7u788MMPPPbYY2zYsIH//Oc/5/38VwTkiv/Ny8vjiiuuYNeuXfTt25e//OUvHDhwgM8//5xhw4Yxc+ZMHn744TOOOWHCBH7++Wdmz57NSy+9VOV8//73v3F3d+eee+6p3F6X57xiuZCfnx933XUXXl5efPPNNzz33HMsW7aMFStW4OLiUuUc11xzDYcOHWLEiBGcPHmS+fPnc9ttt/Htt9/yr3/9i82bN3P99dfj5ubGF198wWOPPUb79u0ZNWpUlcd6+utt9uzZLF68mJiYGIYOHUpJSQm//PILb7zxBt9++y2//PILrVq1OuvP5ZylppKdkFv/4wKYzczmzZtNwPzmqlGOLkVE7Kjid33z5s3ndFxBgWlu3tzwt4IC+zwPgwYNMg3DMN9///2z7nvnnXeahmGYL730UuW2F154wTQMw/zxxx+r7G+1Wk3DMMx7773XZvvevXvNwsLCKvsvXbrUdHJyMh966CGb7fHx8aZhGKZhGOa///3vamsLCQkxw8LCaqzdMAxz0KBBNtueeeYZ0zAM8+abb65ST2FhoXngwIEaxzt9bIvFctb95s6daxqGYY4aNco8efKkzX1TpkwxDcMw33zzzcptSUlJprOzs+nn52dmZGTY7D969Ohqz/vDDz+YhmGYAwcONPPz823ui4uLMw3DMJ944gmb7SEhIaZhGObw4cPNEydOVG7Pzs42fXx8TG9vb7OoqKhy+9q1a03DMMzIyEgzLy+vcvuJEyfMyy+/3DQMo8rPouKxOzk5mcuXL6/y3FS8jp588kmztLS0cntJSYk5fvx40zAMc/HixVWOq87dd99tGoZhzps3z2Z7QUGBGR0dbRqGYb744oumaZrm/fffbxqGYT7yyCM2++7atcts06aN6erqalqt1srtFa/F2NhYm8fdrl07s2PHjmZxcbHNOCtXrjQNwzDHjBljs/1cn/Offvqp8nk99XVZXFxsDh8+3DQMw5w2bVq157jppptsXt9r1qwxDcMw27RpY1522WXm4cOHK+9LSUkxXV1dzd69e9uMVd3jNk3TTE1Ntfl5VXjvvfdMwzDM6dOn22yveB2c/rM5Va3em7/91uzDe3V6Dz+bZhts50b9xdGliIgd1TXYbt5smtDwt3p+767UtWtX0zCMasPG6SZOnGgahmE++uijldvqEmzPpFu3bmZ4eLjNtop/VPv06VPjcecabIuLi01vb2/Ty8vLzMrKqnV9NY1dm2Dbq1cv083NzSZInFpPu3btzEsvvbRy24svvmgTwk6VmppqOjk5VTnvLbfcYhqGYe7cubPGGvz8/Gy2hYSEmBaLxUxOTq6y/9ixY03DMMzt27dXbrvvvvtMwzDMuLi4KvuvWrXqjMF2xIgRVY4pKSkxfX19zYCAgGpDUl5enmmxWMzbbrut2sd0uopge8stt5gvvPCC+cILL5gPPfSQGRAQUBnIDx06ZJ48edL08PAw27RpYx46dKjKOJMmTTINwzCnTp1aua2mgFfxu7Fw4UKb7aNGjTINwzB/+uknm+3n+pxXhPs5c+ZU2X/Xrl2mk5NTld+binOcGswrhIWFmRaLxYyPj69y36BBg0xXV1ebn0VNj7smpaWlZps2bczBgwfbbK+3YPvvf5uBLLZLsG22SxH2H2i+y4dFpO66dIHNmx1zXnuqzbf6K/YpKSk57/N99NFHxMXF8fvvv3Po0CGbMWv6GPuyyy477/NWSEhIID8/n379+tGhQ4d6G7cmx44d4/fff6d9+/Y2ay1P5erqSkJCQuV///rrrwAMGDCgyr7BwcEEBQWRlpZms339+vW4uLjw2WefVbsutbCwkAMHDpCXl8cFF1xQud3Hx6fadmlBQUEANuslt2zZgmEYXHXVVVX279+//xlbn/Xt27fKtt27d5OXl0e7du1q7JPq7u5u89zUxpdffsmXX34JgKenJ6GhoYwZM4ann34ab29vtm3bxokTJ7j00kur7ZAwePBgpk2bVvlzOJOHHnqIGTNm8N577zFixAgAsrOzWbRoEdHR0fTv37/KMefynP/6668YhlHtOvGoqCgCAgJISUkhPz/fZpmJj48PoaGhVY7p1KkTqampXHzxxdXeV1RUxL59++jYseMZH3dRURHvvfcen376KTt27CA/P5/S0tLK+/fu3XvG4+vKTM9gP93sMnazDbbZh93KJkrUwkVETuHpCX36OLqK+tOhQwcSEhKqBKTqZGRkAFS7/vVcPPnkk/zf//0fnTp14rrrriMgIAAPDw+gbM1vTbXUZwA9dOgQcP6PpbYqQsqBAweYOnVqjfud+gdGxRex/P39q93X39+/ynOVk5NDSUnJGRvpG4bB0aNHbYJtTa2vKtb7nvqHx5nqcnZ2PmNLqOp+hjk5OQAkJibW+NwYhkFBQUGN41YnLi7OpivC6SoeR02vq4rttel7GxYWxrBhw1i+fDkpKSmEhoYyd+5ciouLefDBB6s9pi7PeU21duzYkYyMDA4fPmwTbM92jtatW9d4X1FRUbXHnur2229n8eLFdO7cmREjRtChQwfc3NwwTZM333yz2i+11Yec5MMUUf0fwOer2Qbb/aVt4eBBaN/e0aWIiNjNwIEDWbVqFStWrGD8+PE17ldSUkJ8fDyAzexTxexccXFxlWMqwuOpsrOzeeutt+jevTvr1q3Dy8vL5v6PP/64xhrqs1esj48PYL8ZpdNVBIw+ffqwadOmWh1TEVD2799fbbeK/fv313iegwcP1rXUs6o4x759+6rMOBYXF3Pw4EGCg4OrPba6n2HFeLfeeisLFiyo52prdurjqE5WVpbNfmczYcIEli1bxuzZs5k2bRqzZ8/G09PzjOH6XGvNysoiPDz8vGutD5s2bWLx4sVcc801LF26FCcnp8r7TNM8Y3eH85WRXH0HkfrQbD+vz8YfymcnRESaq3vvvRcnJycWLlzIzp07a9zvgw8+ICsri6ioKJuPoCtm/aqbZa0uwCUnJ2OaJkOHDq0SajMyMkhOTq7T43BycjqnJRJdu3bF29ub33//vcZgU59atWpFdHQ0f/zxR63bIPUp/2hgzZo1Ve5LTU2ttjXa5ZdfTm5uLjt27Di/gs/g4osvxjRNfvzxxyr3/fTTTzYfRddG165d8fHxYf369dX+gWQvXbp0wcPDg99//73aWdmKP+T61PIjmuHDhxMSEsLcuXP59ttvSU5OZtSoUTYzqHXVp08fTNNk1apVVe7bs2cPGRkZhIWF1cu5aquiHd7NN99sE2oBfv75Z06cOGG3c2fstd+n6c022O7HD6p50xARaU5CQ0N59tlnK9stVRduFy9ezOOPP46Hhwfz5s2zmXXr168f8GfbsArp6enVfqxcMcO3Zs0amwB09OhR7r///jqv323bti3Z2dm1/ujTYrHwyCOPcPz4cSZMmFDlY9fCwsJ6n/X8f//v/1FYWMi4ceOqDVJ5eXk26znvvPNOnJ2defvttyuXgUDZbNgzzzxTbYCsaHd2//33V87inaqgoICff/75vB5HRduqadOm2YT0EydO8Mwzz5zzeE5OTjz22GNkZWXxv//7v9UGoqysrDP+4VUXLi4ujBkzhvz8fCZPnmxzX1JSEm+99Raurq7cddddtRrPMAwefPBB9u3bx/333w+Urb2tD+PGjQPgpZdesnldlpSU8Le//Q3TNM/4iYs9VPwuV/wBUCE7O5tHHrFvv9SMbFcs2OePoGa7FKFsxrb6ZssiIs3JlClTKCgo4PXXX6dzJSHRAAAgAElEQVRnz54MGzaMiy66iKKiItatW8fGjRvx8/Njzpw5Vb78c+mllxITE8OqVau47LLLGDRoEPv372fJkiUMGzaMzz//3GZ/f39/7rjjDj799FN69erFkCFDOHz4MN9//31l8/jffvvtnB/DNddcw6ZNm7juuusqe+P26tWLG264ocZjXnjhBX7++WcWL15MVFQU119/Pa1btyY9PZ3vv/+eGTNm1PpjZNM0bfqUnsowDGbNmsW9997L5s2bmTVrFp07d2bYsGEEBQWRm5uL1WplzZo1jBs3jlmzZgEQHh7O1KlTefbZZ+nZsye33347bdq04fvvvycvL4+ePXuybds2m3NdffXVTJ8+nWeeeYbIyEiuv/56QkNDOXr0KKmpqaxevZqBAwfy7bff1u6JrcYVV1zBY489xttvv023bt0YOXIkLi4ufPnll7Rt25aOHTue8wUVJk+ezO+//867777L119/zaBBgwgICCA7O5vExMTKHsbVLck4H9OnT2fNmjXMnDmTX375hZiYGA4ePMjnn39OQUEBM2fOJCQkpNbjjR8/nilTppCVlUXPnj3r7QuPl19+ORMnTuTVV1+lW7du/OUvf8HT05OlS5eyfft2Bg4cyN///vdzGvNcf0anu/TSS+nfvz8LFy6kf//+9O/fn/3797Ns2TK6dOlCp06d6uXCGlUcOcLeE7605wBVF+PUg3rtsXCaFStWmGPHjjUjIyNNT09PMyAgwLz55ptr3dph//795t133222a9fO9PT0NC+//HJz5cqVZzymos2EEz+bxU89Wx8PQ0Qaobq2+2rOfv75Z/Puu+82Q0NDTXd398resWFhYdW2Q6pw+PBh88EHHzT9/PxMNzc3s3v37ubs2bPNlJSUatt9HTt2zHzuuefMiIgI093d3QwODjYfffRRMycnx4yJianSwqo2rYYKCgrMhx9+2AwMDDSdnZ1Ni8Vic97q+tiaZlmbrZkzZ5qXXXaZ2apVK9PLy8uMiooyH3zwQTMpKalWz1tFu6+K56viVrHNYrHYtPhasmSJecMNN5h+fn6mq6ur2bFjR7Nv377m5MmTzV27dlUZ/z//+Y/Zu3dv093d3fTz8zPvuusuMzMz04yOjjZ9fX2rremnn34yR40aZXbq1Ml0dXU1/fz8zN69e5t//etfq7zmQ0NDa2yVNmXKFNNisVTbzm3mzJlm165dTTc3NzMgIMB89NFHzcOHD1c7XlxcnGmxWM7Y5qnisQ4ePNj09fU1XV1dzcDAQHPgwIHmK6+8UqWXb03uueeeWp2rwqFDh8ynnnrKjIyMNN3c3MwLLrjAHDp0qPn9999X2bc2r8URI0aYhmGY7777bo371PU5//TTT80BAwaYrVu3Nt3d3c1u3bqZL7/8cpW+yGc7R3W/ZxUqnr/U1NTKbTU97tzcXHPChAmV7xkRERHmc889Zx47dqzOr4Ozvjfv3Gnewwdmdz6wy3u4YZr2iONlRo0axYEDBxg1ahTR0dEcOHCA1157jU2bNrF8+fJq215UOHnyJJdccgn5+flMnz4dPz8/Zs6cyTfffMOKFSu48sorqz1uy5Yt5e0vNpM5ci4dF7xtp0cnIo5U8bu+efPmWq+ha2mOHj3KgAED2Lp1K1OnTmXSpEmOLknK5efn4+/vT58+fVi7dq2jy5FypaWldO7cmZycHDIzM+1z1a1m7qzvzStWMGSICWxmBc/U+3u4XZcizJw5Ez8/P5tt1157LREREbz88stnDLZz5sxh+/btrF+/vvKjs5iYGHr27MnEiRPZsOHs1xhOTynhzB3cRESar1atWrFkyRL69u3L888/T3h4OHfeeaejy2pRDh48iLe3t82lUouLi/nrX//KyZMnGTlypAOrk9N99tlnpKamMmHCBIVae9m7lwwuow9L7TK8XYPt6aEWwMvLi65du9ospK/OokWL6NKli816MCcnJ8aMGcOzzz5LVlbWWRsPZ2RaqL924CIiTU9gYCBLly5l0aJFpKenU1hYiKurq6PLajEWLFjA888/z5AhQwgMDCQ3N5fVq1eTmJjIxRdfzKOPPuroEoWyL3Xl5uYyZ84cWrduXacv0UntmOkZpDOS6+yzwrbhvzx2+PBhtmzZwjXXXHPG/f74449qr4rSvXt3ALZv337GYOvCSTIO6CINIiI9evSgR48eji6jRerXrx8xMTGsW7eO7OxsoOzb6JMnT+app57SHxmNxPPPP4+rqyvR0dG89tprDXbhj5YoP/kgBbTCjwN2Gb/Bg21Fe5bnnnvujPvl5ubi6+tbZXvFtoorndTEn2wyiv0hNxfatq17wSIiInXUq1evKp0lpPE51969UncVF2fws9OMbYP2sZ08eTKffPIJb7zxBr1797brufzYTwaBukiDiIiISCORkV7Ws8DfcuYJyrpqsBnb2NhYpk2bxssvv8yECRPOun/btm3Jzc2tsr1iW9uzzMKm8zoptOWmB4qh/HrYo0ePZvTo0XWoXkRERETqYv78+cyfPx+AtNQ/gBuZ7pwDdriyboME29jY2Mrb008/XatjunfvztatW6tsr2hm3a1btzMefz1DWM2TfDXue3jwwXMvWkRERETOW+XE4okTxHq8wn4e4v+6XMvF1eS882X3pQgvvvgisbGxTJ48ucol785kxIgRJCQksHHjxsptxcXFfPTRR/Tr148OHTqc8Xh/stlLAKXpe+tcu4iIiIjUk8xMMggkkIzKT9Prm12D7WuvvcYLL7zAtddey/XXX8+GDRtsbhXGjx+Pi4sL6enpldvGjRtHdHQ0t912G/Pnz2fFihWMGjWKxMRE/vGPf5z13H7spxA3DibmnXVfEREREbGzjIw/g201LWHrg12XIixZsgTDMFi2bBnLli2zuc8wDEpKSoCybyOWlpbaXJPY1dWVlStXMnHiRB577DGOHTtG7969Wbp0KQMHDjzruf0pa6uSYS3CPk+diDQGO3fudHQJIiJS7ozvyRkZ7KUbV/Fj0wy28fHxtdpv7ty5zJ07t8p2Pz8/4uLi6nRuP/d8OAHpey3oYpsizU/r1q0BGDNmjIMrERGR01W8R9vIyCCDa5vujK0j+fq74JxaRMZBd12kQaQZioyMZPfu3Rw5csTRpYiIyClat25NZGRkle0F1mzy8C1fY9vFLudutsHW4u9HQOpeMgrbw+HD4OPj6JJEpJ5V98YpIiKN0949xwHsOmPboBdoaFD+/gSSoYs0iIiIiDQCGWllV3gLZC+0b2+XczTfYOvnRxDpZcH2lG4LIiIiItLwMvaVLRQI6FgKzvZZNNB8g61mbEVEREQah5MnychvTVsO4h5inx620EKCrZmuYCsiIiLiMOU9bAPYC0FBdjtNsw+2J/AgN0kXaRARERFxmPR00gkimDQF2zrx8yv71h1lF2kQEREREQcpD7ZBpCvY1kmbNgS65wBaYisiIiLiUOnppBGsYFtnhkGHQGecKCYj29XR1YiIiIi0WAXJ+8nDt2wpQnCw3c7TfIMt4BQcQEeySD/ZHvLzHV2OiIiISIuUvuckgGZsz0tgoFp+iYiIiDhYWqoJQLBzlt2uOgYKtiIiIiJiZ+n7XTEoJSDQAIv94qeCrYiIiIjYz9GjpB9vSwf24RLc0a6navbBNoiy9hK6SIOIiIiIA5R3RLB3D1to7sE2qKwR8DG8yE3McXQ1IiIiIi3PqT1s7dgRAVpIsAVITyp0cDEiIiIiLVAD9bCF5h5sfX0J8iibqU1LNxxcjIiIiEjLY6amNcjldKG5B1vDwD/YDRcKSct2B9N0dEUiIiIiLUrunlyO46kZ2/pgCSlb05FW1AFycx1djoiIiEiLkpZUBNj/4gzQAoItwWXfwksjGNLTHV2NiIiISIuSvrdsOWiwx0G44AK7nqtlBdu0NEdXIyIiItJymCbp2e64UIhfsDsY9v3OU4sJtukEKdiKiIiINKS8PNKKOhBIBpbgQLufrsUE20w6UWTVRRpEREREGkxaw3VEgBYSbINIpxQnMncfdXQ1IiIiIi3HqRdnULCtB4GBlRdpSLOWOLgYERERkRakAS/OAC0h2Lq5EeRXdtWxtExnBxcjIiIi0nKUpGawl4CySUY7X04XWkKwBVqH+HIBuaTntYKiIkeXIyIiItIi7Es8QgnOmrGtV5Utv4Jg715HVyMiIiLSIlQsA1WwrU/qZSsiIiLS4NL3lkXNYO988PKy+/laTLANIl3BVkRERKShlJaSnuNJK47gHezdIKdsMcFWM7YiIiIiDWj/ftJKAwgiHSPE/l8cgxYWbA/jQ/6ebEdXIyIiItL8lfewbaiLM0ALC7YA6YknHFyMiIiISAvQwBdngJYSbNu3J9h1P6CVCCIiIiINIi2tQS/OAC0l2BoGHYNdcKKYtP2uYJqOrkhERESkWTuZlEE2/mWfmoeENMg5W0awBZxDAuhEJmkn/eHwYUeXIyIiItKsZew+BpT3sFWwrWendkZIT3d0NSIiIiLNWuXFGSyZ0KlTg5yzxQXbdIK00FZERETEztKyXAAI6lQCzs4Ncs4WF2zVy1ZERETEzgoKSC1oix/78Qjr0GCnbTnBNqisj1oGgZSkaCmCiIiIiN2kppJKCCGkNtj6WmhJwbZ8xrYIV/bv1pfHREREROxGwdbOgsobBANpycUOLkZERESkGVOwtTNPT4J9CwBI2+vk4GJEREREmq/SlLLvNSnY2pFPiDetOEJ6rhcUa9ZWRERExB727zpEIW5lwTY4uMHO26KCrRFS3hnBDISsLEeXIyIiItIspSaX9bBVsLUntfwSERERsbvUjLJlnyFtC8DTs8HOq2ArIiIiIvWnqIjUvNZ4cwifsAsa9NQtMtimEqJgKyIiImIPGRmkmg3/xTFogcE2hFRyaEfBHq2xFREREal3Dmr1BS002AKk7j7p4GJEREREmiEF2wbi70+oa9lMbUqKY0sRERERaY7MlFRSCG3wjgjQ0oKtxUKnEBecKSJ1nxuYpqMrEhEREWlW8vbkcJTWmrFtCE5hwQSRTkphR8jNdXQ5IiIiIs1KxXJPBduGEBJCKCmkEKr1CCIiIiL1LDW9LF6GeB6EC5phu6+jR48yceJEhg4dSvv27bFYLMTGxtbq2Li4OCwWS7W37Ozscy8mNJRQUspafinYioiIiNSf0lJSsz1w5zh+oZ5gGA16eueGOMnBgweZPXs2vXr1YsSIEbz//vsY5/hA4+Li6NKli802X1/fcy8mNJQQdvEt10PK2nM/XkRERESql51NakkAwaRhhDbsMgRooGAbGhpKXl4eADk5Obz//vvnPEa3bt3o06dPfRRDKMvZTweOJ2Xicf4jioiIiAjYtvpq4I4I4IA1tmYdOxHU9bgqytfYAqQlHKufMUVERETEoT1soQl9eeyGG27A2dmZtm3bMnLkSLZv3163gTp2JMQ5E4AUq9p9iYiIiNSbtDRSKZ9EdECwbZClCOejY8eOTJo0iX79+tGmTRu2bt3K9OnT6devH+vWraN79+7nNqDFQmCIE05JxaRklfeybeCFzSIiIiLNUcGeLA7SvnzGdniDn7/RB9thw4YxbNiwyv8eMGAAw4cPp3v37jz//PMsWrTonMd0Dg0kMCmD1JP+kJcHdfkSmoiIiIjYSNt1HHBMD1toQksRThUSEkL//v3ZsGFD3QYob/mVQiikptZrbSIiIiItVWpq2TLPEOdM6Nixwc/f6Gdsz+RMLcOeeOIJfHx8bLaNHj2a0aNHl7f8SiWRyLJetr1727lSERERkeYvNcsNJ4rpFGgBS9n86fz585k/f77NfocOHbLL+ZtksE1OTmbNmjU2SxRO9+abb9bcHiw0lFCSWME1kPKznaoUERERaUEOHyb1hB+BZOAcGli5uXJi8RRbtmzh4osvrvcSGizYLl26lIKCAo4cOQLA9u3bWbBgAQDDhw/Hw8OD8ePH8+GHH5KcnExQUBAAQ4YM4eqrryY6OppWrVqxbds2Xn31VZydnXnxxRfrVkxICKGsJJMATu5Jx61eHqGIiIhIC+bgVl/QgMF2woQJpJavZzUMgy+++IIvvvgCwzCwWq0EBwdTWlpKaWmpTc/a7t278/HHH5Oens7x48fx8/PjmmuuYfLkyURERNStmPKlCFC2yDnyvB+diIiISAtntZJKCOEkQ1iYQ0posGBrtVrPus/cuXOZO3euzbbXX3+9/ovp1IlQpwwogVRrqYKtiIiIyPmyWkmlD4OId1iwbZJdEc6bkxOBwRYslJCS5eroakRERESavKI9qewloOxTcQXbhuUa2olOZJJy3B/s9M08ERERkZYiLeEYJpayq46FhjqkhhYbbCt62aYSol62IiIiIufJmlwKQJhzBnTq5JAaWnywTSG0rJetiIiIiNSNaWLNdMdCCcEhBjg5OaSMFh1sQyhrS6FgKyIiInIecnOxnuxIIBm4hAc5rIyWG2xDQgglhb0EUJiU7uhqRERERJouqxUrYYRhddj6WmjJwbZ8KUIpTmTsKnB0NSIiIiJN16nB1kEdEaAlB9uAAEItZTO1qcklDi5GREREpAlTsHUwZ2fKr9pLSqZ62YqIiIjUVUFiJtn4aymCI7mFdqQTe0k51h4OH3Z0OSIiIiJNUsrO4wCasXWoU1t+qZetiIiISJ1YUwwAQt33g5+fw+pQsCUFK2Fq+SUiIiJSF6ZJyn4PXDlJpzA3MAyHldLig204yQq2IiIiInW1bx/W4kBCSMUSFuLQUlp2sA0PJ5xk9hLAid1pjq5GREREpOlpJB0RQMGWMKyYWEjdoV62IiIiIucsJUXBtlHo1IlwlwwArEmlDi5GREREpAlqJFcdg5YebC0WAkJdcKGQ5CwPME1HVyQiIiLSpOQl7OcwPpqxbQycOocSQirWogDIznZ0OSIiIiJNinVXIeD4HragYFu5zjaZcEhOdnQ1IiIiIk2KNbUsToa1zoELLnBoLQq2YWF/tvyyWh1djYiIiEjTUVKC9WArWnGEtmFtHF2Ngm1Fy68kOmMmacZWREREpNYyMrCWhhCGFSPcscsQQMG2cilCPt7kJex3dDUiIiIiTUd5q69QUhy+vhYUbCuXIsCfi59FREREpBYa0cUZQMEWvL0J8zkEQHKqng4RERGR2jKTraQQ2ih62IKCLQAXdPbFm0NYD7aBQs3aioiIiNTGvh25nMBDM7aNidG5ouVXGKSmOrocERERkSbBuqcEQDO2jUr5OttkwtXyS0RERKSWrGlOAIT55kOrVg6uRsG2THnLLythukiDiIiISG2cOIE1z5u2HKR1ZAdHVwMo2JYpb/mVSgglezRjKyIiInJWp3ZE6NzZ0dUACrZlypciFOHK3p35jq5GREREpPFLSlKwbZSCgwkzyr40lpxY4uBiRERERJqApCSSKVvOSXi4o6sBFGzLuLgQElSKQSnWDBdHVyMiIiLS6BUmppJOEJ1J0oxtY+MeEUgAe0k+3gHy8hxdjoiIiEijlvLHUUpxUrBtlMLCynvZquWXiIiIyNkkJZX9b2e3vdCxo2OLKadgW0Etv0RERERqp7SUpH2euFBIYLgrGIajKwIUbP8UHv7njK2CrYiIiEjNMjNJKg4hDCtOEY6/lG4FBdsK5S2/9tOBY7szHF2NiIiISOOVlEQSnRvV+lpQsP1T+YwtgHXnCQcXIyIiItKIKdg2cu3aEe65H4Bka+NYJyIiIiLSGJXuSSaZ8LJg20h62IKC7Z8Mgw6dvXDnOEn7W0GJLtQgIiIiUp2s7bmcwEMzto2ZJTyUziSRVBoK6emOLkdERESkUUraXTYB2JlkCA11bDGnULA9VUQEEexhDxF/NmcTERERERtJ6a4AhAUUgpubg6v5k4LtqSIiymZs6Qx79ji6GhEREZHG59Ahkgr8CSADj8hAR1djQ8H2VOUztlbCKE7U1cdEREREqkhObpQdEUDB1lZ5sC3GhbRthx1djYiIiEjj00hbfYGCra2gIDo7pwF/LooWERERkVOcGmwbUasvULC15eREcJgTzhSxZ68HmKajKxIRERFpVA7tzCKXtpqxbQqcI8MIw8qeomDYt8/R5YiIiIg0Kkk7TgIo2DYJp7b8UmcEERERERtJKU4AdPbOgQsucHA1thRsT6eWXyIiIiLVKywkKccbH/LwjfB1dDVVKNiernzGNonOlCbqIg0iIiIilVJTSTLDG+UyBFCwrapzZyLYwwk8yPwj19HViIiIiDQeSUkko2DbdISG0tkouzjDnl1q+SUiIiJSqRH3sAUF26pcXQkLKcWgtOw6yGr5JSIiIgLAyV0ppBPUKHvYgoJttdwigwkmjT3HO0FOjqPLEREREWkUUv44iomlLNhGRDi6nCoUbKujll8iIiIiVSQllgLQ2W0vBAQ4uJqq7B5sjx49ysSJExk6dCjt27fHYrEQGxtb6+Ozs7O55557aN++PV5eXlxxxRX88MMPdqyYypZfe4iAJHVGEBEREaGkhKQsT9w4QUCEB1ga3/yo3Ss6ePAgs2fPpqioiBEjRgBgGEatjj158iSDBw8mPj6et956i6+++gp/f3+uvfZaVq9ebb+iyzsjJNEZM1EztiIiIiKkpbGnJJQwrFiiGt8yBABne58gNDSUvLw8AHJycnj//fdrfeycOXPYvn0769evp2/fvgDExMTQs2dPJk6cyIYNG+xSc9lShDiO0IYD27Pxs89ZRERERJqOxER2E0UUuyEy0tHVVKtB55DNc+wwsGjRIrp06VIZagGcnJwYM2YMGzduJCsrq75LLBNe3p8N2LOzyD7nEBEREWlKEhNJJJJIEhVs6+KPP/6gR48eVbZ3794dgO3bt9vnxB4edO50AoCkNBf7nENERESkCSlKSCKFUAXbusrNzcXXt+p1iCu25dixFZdXVAAdyWTPET84fNhu5xERERFpCqxbj1CCs4Jtk1T+BTJ1RhARERGBxD1lX/6P9NgLHTs6uJrq2f3LY+ejbdu25ObmVtlesa1t27Y1HvvEE0/g4+Njs2306NGMHj26dicvb/m1g4vKetn26VP7wkVERESak+Jidu9rgwfHCIj0hFp2uAKYP38+8+fPt9l26NCh+q4QaOTBtnv37mzdurXK9m3btgHQrVu3Go9988036XM+YTQiggh+5WtuhD0r6z6OiIiISFOXkkJiaTgR7DnnVl/VTSxu2bKFiy++uD4rBBr5UoQRI0aQkJDAxo0bK7cVFxfz0Ucf0a9fPzp06GC/k5dffSyHduTtsFP3BREREZGmoAl0RIAGCrZLly5lwYIFfP3110BZN4MFCxawYMECjh8/DsD48eNxcXEhPT298rhx48YRHR3Nbbfdxvz581mxYgWjRo0iMTGRf/zjH/YtunPnsh8ekLi90L7nEhEREWnMmkiwbZClCBMmTCA1NRUou+rYF198wRdffIFhGFitVoKDgyktLaW0tNSm162rqysrV65k4sSJPPbYYxw7dozevXuzdOlSBg4caN+iW7cm0i8fsmG31YXL7Hs2ERERkUbrxE4raQSXB9sBji6nRg0SbK1W61n3mTt3LnPnzq2y3c/Pj7i4ODtUdXatuwTQMTuT3Yf9ID8f2rRxSB0iIiIijpS07RgmlkY/Y9uo19g6XFTZZeN2EwWJiY6uRkRERMQhKlp9RXllgp+fg6upmYLtmZwabHfvdnQ1IiIiIg2vsJDE7Da04gj+Ud7n1OqroSnYnskpwdbcpWArIiIiLVBKColmBJEkYkSeW6uvhqZgeyblwbaAVmT9tt/R1YiIiIg0vCbSEQEUbM8sPJwoYw8Au3eWOLgYEREREQdITGQ3UQq2TZ6bG+GhpVgoYXeqG5zSikxERESkJSjYkUomAUSxW8G2qXO9MIwwrOw+GQz7tRxBREREWpY928oupqUZ2+ZAnRFERESkBUtMKouLka33Q7t2Dq7mzBRsz+bCCxVsRUREpGU6eZLdBy/AhzzaRrVt1K2+QMH27MpnbJPoTHHCHkdXIyIiItJwkpP/bPUV1biXIYCC7dmVB9tiXEj5Nc/R1YiIiIg0nPJWX03hi2OgYHt2gYFEuaUBWokgIiIiLUxCwp89bC+80NHVnJWC7dlYLARGeuDOcXZntoLiYkdXJCIiItIgDm9LIxt/BdvmxHJh2V8qu0s7Q2qqo8sRERERaRCJW09p9RUV5eBqzk7BtjbU8ktERERaoF3JLgBc2PEItG7t4GrOTsG2NhRsRUREpKXJyWHX0U50JJM2XQMcXU2tKNjWRnmwTSeYY9utjq5GRERExP527SKBLlzIriaxvhYUbGunPNgC7Nl6zMHFiIiIiDSA8mDbhQQF22alXTuifA4AsHuPnjIRERFp/kp27iaRyLIZ2y5dHF1OrSil1VLbC9txAbnszvGFY5q1FRERkeYt7dccTuDRpGZsnR1dQFNhXBhF1M/lXyDbswd69HB0SSLSnJWWwtGjkJ8Px4+DkxM4O4ObG/j6gouLoysUkWZuV4IJQBdXKwQHO7ia2lGwra3ydba7uLCsM4KCrYicL9OEvXthyxbYsoWSnbvZt+coqdZS9uZ5kk9rjtCaY3hioRRXCnGlkLbk4N/qGP7tSwnr4oZndFjZx4QXXwzdupUFYBGR81FcTEJmG9w5TnCUO1iaxof8everragoLuQ3lnADZsIPGI6uR0SapsxM+P57Sr5byW/f7eeng134jV78yi3s4CKKcLXZ3ZMCvCigFAtFuHACdwpxg6PAUTCspYQuTeEidtCHr+jv9gL9+pp4D+0Lw4dDz55g6B1LRM6R1cqu0gii2I2la9NYhgAKtrV34YV04VPy8OXA75n4OboeEWk6MjNhwQJyPlrK4l868Q3DWcWb5OGLGyfozjYuYyP38T5hWAlud5zAQGjT1gUn71bg4VG2NKGkBI4f5+j+Avbtg337ILEwmB1cxA4uYhYTePFkO1CM2h0AACAASURBVIzVpVyyehM3TvovN/lNpMdtF2KMvQsuvVQhV0RqJyGhybX6AgXb2ouKoisJAOzcWqRgKyJnVlIC33xDwVtz+HxlWz7ldlYyARODK1jH4/wfg93XcVl/F1wv6wV9+kC3RyE0FNzdzzh0KyACiCgtZcDevbBrF2zbhrn+Q3av3sdP+yNYzjD+yd95PrsN4e8kMfadDxkb9gxh918D998P7do1xLMgIk3Vrl0k8D8MZI2CbbPk7k5EWAlO1mISUty5yjQ18yEiVeXnw7/+xfY3v+fdfTfzIR9yhNbEsIq3eYxbL9yB/21XwtCh0Pc5cHU9+5g1sVggKKjsds01GE/ChcCFVivjv/2Wwq/H8OMPJcwvGskM/sYUayyDnv2Bx194mBvuugCnJ/+3bE2uiMhpDm9LYx8dy2dsr3N0ObXWNFYCNxKuF0XQmSR2FoZDRoajyxGRxiQvD6ZM4aeA27nu6R5027eCL7iNx3gba+CV/PD8jzz8x6P4J/wIL74IAweeX6g9k7AweOQRXJd9xZADn/DBHNg/8DY+4n84iRu3FH3BhR9M5O3u73H8ptth61b71CEiTdaurScBmlSrL1CwPTddu9KVneykK+zc6ehqRKQxOH4cpk0jPmAMMbExDDy6lAwC+Zj/Ie26h3jpmz6EpPwIsbEQHd3w9Xl7w7hxeK5exv/smcraJxawwfNqLmETT/IGnb9+g3d6vkfhX0ZDQkLD1ycijdKupLIP9aP886FNGwdXU3sKtueiS9ll5RLoon8ARFq60lL4+GN2hg/nhkk9ufr4NxylFYstI/j9nje5c/cUXL9dDNdfX9aDtjHo3BneeIO++77k0zf3/3/27ju65vt/4PjzZgqCik3tkciwt6ANQq2i+MZqQ+0VI7YkktgrRqitSmPVrr1TW4wMe5bGjhUSkdz7+yPaX7XaGrl53/F6nOP0nMtNnpxTefnkPbiQqxZ12UUfZlLyp3H84Dwerc8AePxYdakQQqVHjzj/LB/5uYW9UwHVNe9FBtv38fqJ7a8UIj7yquoaIYQqUVHEVWlIz/ZPcL2zg7OUZqXmfxzvPJdml6disXghlCihuvKf2dtDv34Uu76bpSGPiHaoQ0VO0FG7BPfpLTlZpCUsWJA6vAshzM+FC1yglNEtQwAZbN/P6ye2ABdPPVccI4RId4mJ6IaPYEXZ8TidWMpy2jGeoZzz7E/raD80C+anrm81FnZ20K8fpW9sZU1ADHtsGvCErFR8vJMeXV7xpFaT1JsWhRDm5cKF/z/qy9FRdc17kcH2fWTPjmOOhwCcuyzXWQphVg4f5kbphjQaVwMv7XJqcYBzxZowaHt9bLdtgNKlVRd+uEyZwN+fzy7N49RXY5lGf5bTDpeD37HVeRBMm5Z6fJkQwiyknL3AJUrIE1tzkNW5AHmJ5fzTvKm7oIUQpi05Gd3oQJbVmIPbtfVE48JGy+as9osmX8zO1GO7TEXBgliv/pF+O5sQXaAhTpzji6T1dBqQlcfuTeDmTdWFQoh0cP3UI5KwNbrLGUAG2/cnJyMIYT6uX+dxzca0DShBB91SmrKRqAreNDkTnHrKga2t6kL9qFuXgue2s73Pz8ynCz/RknKHQznm0gk2bFBdJ4TQs/PndAA42l6HQoWUtrwvGWzfl5OTnIwghDnYtInDLl1wOzqPrTQkTNOWHwKukvXIdjXHdqW3zJnRzJjOtwc6cjpfI3JzlxpPtzDlywNoe/eFly9VFwoh9CEpiQu/ZcaOFxRwzGw4p7q8Ixls35ejI06c4yIlSY65oLpGCJHWtFp0fv6ENt1G7ec/U5Bfify0Mf/7pTf4+4OVmV3Y6O5OkehNhH85lf5MYxBTaBLqSVyNJhAbq7pOCJHWLl3ivK4kJbmIRWnj2jgGMti+v9dLEV5hw7WTssZWCJPy6BEvvviKr4OK0ZtQejKbvS1mUTB6C1SvrrpOnU8+wXrtSibOycIWq6YcoSqVI2YTU6YtHDmiuk4IkZbOpS63LM1Zo9wUK4Pt+ypQAEe7XwE4d16jOEYIkWYuX+Zm+WbU2D6KNXzFck17Qia+wnpNmFHduqM3Gg10707DY6M5ka8ZdiRQ9cEmNrpPgkWLVNcJIdKILuYsMTjLYGs2NBryOWXFnqecu/sJJCaqLhJCfKyDBzlZsStVrq/gEZ9wJGsD2u70Bl/f1IFO/L9y5Shyeh2Ha/hSj518mbyasZ0voxs+AnQ61XVCiI9079RvPCJ76mDr5KQ6573JYPsBNKVfbyDTlYKLF1XnCCE+xqpVbK4zmVpPNlKAWxwp+TVuZ34ADw/VZYYrZ04y793Emp578SOQEYyl67jCJHfwhqQk1XVCiI9wNjIZgNKWF6F4ccU1708G2w8hR34JYRqmTiW0zX6aJa+hPjvY91kgeY5tNLrjbZSwtsYidCYBM3OwGG8W403z5S153qAlPH2quk4I8SGSkzl7MzPWJFG8pAVYG99lVDLYfojXV+uexxHdOTnySwijo9OhGzmKoIGP6E0o/ZjOmk5bybh9HWTNqrrOuPTuzTdrm7LZugV7+QyPvSO4X+NLuH9fdZkQ4n1du8bZ5JKU4gJWzsZ1McPvZLD9EK+f2D4hG3dO3VZdI4R4H1otur79GDQmG34EMZZhTAmIx2LBPKN8OmEQmjenwb6h7M/SlGsUoVZ0KLHVv4Lb8vejEEbl7FnOUtpo19eCDLYfpnhxHC0vA3A+OllxjBDinSUnk+L9Ld/OKsNUBhJKT4bNyIfG3082iX2s6tWpcHQ2B3O35DmZqHV5Ib9WawO//qq6TAjxrv482BrhiQggg+2HsbamWDGw4hXnfs0EKSmqi4QQ/yU5meR2X9N2qSff8zU/aDrSc3Fl6NNHdZnpcHSk+OEfOJC/LVosqHVjKVertYUrV1SXCSHewYNTN7lHbhlszZF16RKU4BLnkovDjRuqc4QQ/yYlheSOnWi/qglracFqSy/ar24G33yjusz0FClC4SMrOFDkG2xIolbsCi5U95bhVggjcO5M6qkmpTXnoWRJxTUfRgbbD+XkhDMxxOAMZ8+qrhFC/BOtlpROXfg6zJM1fMVKy3Y03/ANtGypusx0FShAgUOr2F+qG1l5wmf3VnCllrcsSxDCkGm1nL1mhyXJlCiaAhkyqC76IDLYfihnZ5yJIRoXiIlRXSOEeButlpQu3fFeWoeVtCHMoj0t1raHRo1Ul5m+PHnIG76KPSV7YM8zPGKX8mut9hAbq7pMCPE2v/7K2VfFKcElbJxLqK75YDLYfigXF5yJ4T65uH9CliIIYXB0OnR9+9FtUWWW045lmo60WtUKmjZVXWY+cuYk9/5V7C78LRp01L2xgDt1/gf37qkuE0L81blzRr9xDGSw/XClSuGiSV2CEHP6leIYIcTfBAczNLQAC/mWJZpO/C+smSw/UCFPHgoc+JHd+b/mBRmpe2k2Dz5vDU+eqC4TQvyZCZyIADLYfrgMGSheQoM1ScRcl5MRhDAoc+cy2e8JExlCCP3osMQD2rRRXWW+Pv2UouHfsyt3e+6Tk4Yxk4hv4gUvX6ouE0K89vjUNWLJb9Rn2IIMth/F2qUUpbhATHJJuHpVdY4QAuCnn/i++2F8mcwIguk3uSB07Ki6ShQpguOBeWzP9j8uUIpW4X141d4btFrVZUII/nQiAmfB0VFxzYeTwfZjyAYyIQzL3r1s/t8yOrOALswjaNBTGDhQdZX4XcmSlN02nnW2XuzGg65r6qHrPwB0OtVlQpg3nY6zV2yxIIWSnyZC5syqiz6YDLYf4/UGshic0UVFq64RwrxFR3O0cRCtkn+kGRuY0/EImokTVFeJv6pSBY+1vVis6cwSvPGb4QCTJ6uuEsK83b7N2YTCFOUqds5FVdd8FBlsP4azMy5EE4cDdyNuqa4Rwnzdvcv1Bt1p+iKMCkSwvOFyLBfOk2tyDdUXX9Bu4edMxJdgRvHd4CuwYoXqKiHMV3Q0ZymNMzHg4qK65qOky2AbHx+Pj48P+fPnx87OjnLlyrFy5cr/fN+SJUuwsLB46497hnBcTIkSOFtdBCAmUjaPCaFEQgJPGrWl0W9zyUw868v4k2HNMrCyUl0m/o23N4OCstGPEHoRytYOP8LRo6qrhDBPUVH/fyKCkQ+26fI3f4sWLThx4gQTJkygZMmSLF++HC8vL7RaLV5eXv/5/iVLluD4l4XM2bNn11fuu7OxoVhJS2zPJhJzIzMer16BtbXqKiHMh1bLqw6daBUxhFjycTh3c3Js+REyZlRdJt6BZsRwplztytXFP9MmeRmHGzXB+eQPULCg6jQhzMqzU5f5lUI4cQ5cPFXnfBS9D7Zbtmxh165dhIWF0eb1cTu1a9fmxo0b+Pr60qZNGyws/v3BsYuLC+XLl9d36gexdHHC8ex5orVOcPmyUR+RIYSx0Y3yo/dPddjLZ+zI0AzH7dMhXz7VWeJdaTRYfhfK8ovNqHGwME0eLubYFx3JcWSzUW9eEcLYRJ9IBMCVaKM+wxbSYSnCunXrsLe3p1WrVm+87u3tTWxsLEff4VtPOkPeMeviggvRxOAsJyMIkZ6WLydk7HPm0Y15mu58tqYXlCmjukq8Lxsb7DcsY2OhvsSTmZYxo0ny+lrOBhcivaSkEHU1E5Yk41jsFdjZqS76KHofbKOjo3FycvrbU1lXV1cAYt5hGGzcuDFWVlY4ODjQsmXLd3pPunl95JecjCBEOjp5kt3eyxjEZAYxCe+QMtCokeoq8aEcHCi87TvWZerAEarSY/MX6IaPUF0lhHm4do3oVyUpwSUyuJVUXfPR9D7YPnz48K3rYX9/7eHDh//43rx58zJy5EgWLlzIvn37CAoK4vjx41StWpWoqCi9Nb+X10d+PSEbsRG3VdcIYfoePOB6kz60efUDHuxmfKdL0KeP6irxsRwdqbFuEPM03VlEZ2ZMTIDVq1VXCWH6oqKIwhVXouD1Q0djZtDHfXl6ehIYGMgXX3xBzZo16dmzJ+Hh4Wg0Gvz8/FTnpSpWDGeby4CcjCCE3iUn8+KrjjSPnUUWnrKi4hQsZ8+UY71MRb16fB1SjgFMYRCTCe84H86eVV0lhEnTRUX//2Br5CciQDoMtg4ODm99KhsXF/fHz7+PQoUKUaNGDY4cOZImfR/N0pIijrbY8YLoW9nk7nMh9Eg3eAjf7m/PRUqy3uFbsm9YDLa2qrNEWurThwltI6nBQVolLiW2STd48kR1lRAm686JWzwkBy5Em8Rgq/dTEdzc3AgLC0Or1b6xzvb3pQQuH/iHqPmPJzQ+Pj5ky5btjde8vLze6Xix92Xh6kzpyLPE6Jzg4kWTeJQvhMH58UemToMw2rLSsi1uG4LkBARTpNFgNX8OK083psLZpbS6Op69HTpjs34V/McJOkKI9xd1OvW7za7WF6BECb18jrCwMMLCwt547fHjx3r5XHofbJs3b878+fNZs2YNrVu3/uP1JUuWkD9/fqpUqfJeH+/q1auEh4fj6fnv56yFhISk3xFhf9pARnS0DLZCpLXoaPZ5f89gtjCE8bSeVQtq1FBdJfQlY0Zyb1rA6rKdqf1sE4M21WLGuHEwQjaUCZGmXr4k+lY2MvKcok62ervY5m0PFk+ePEmFChXS/HPpfbBt0KAB9erVo0ePHjx9+pRixYoRFhbGjh07WL58+R9PXjt37szSpUu5evUqn376KQD16tXj888/x9nZmcyZMxMVFcXEiROxsrIiKChI3+nvzsUFZ8JZSwu00T8b9sJlIYxNfDx3m3fHK2k1tdnPGO8r0G2e6iqhb0WLUm1Vf0Ia9qcXoVQe2YH2lXZA/fqqy4QwHefPE6VLfThn4Wb8yxAgnW4eW7t2LSNGjMDPz4+4uDicnJxYsWLFG09wtVotWq32jTNrXV1dWb58OTdv3iQhIYFcuXJRt25dRo0aRfHixdMj/d04O+NGKPHYc+PoHYqo7hHCVOh0pHTrSdvLo9Gh4cfSY7Cc/bNsFjMXDRrQIyiCo6O+pytzKfc/T5xjXCFvXtVlQpiG6NSNY25EmsT6WkinwTZTpkyEhIQQEhLyj79m8eLFLF68+I3Xpk6dqu+0tFG4MG4Zr8ALiIzSyGArRFpZuJCgH4uyjzrstGtGnvXfQYYMqqtEOtIMH8acX77i5PbytH70Hcf/14mMezaDpaXqNCGMXkpkDDG0oB3LwaWu6pw0Id81TwsWFuR1y4kDD4i8l1t28AqRFiIj2dVzLYH44c9oPl/cQW8bG4QBs7Ag47J5rMzZh2sUoe+BlhAcrLpKCJNw9dgDErEzmTNsQQbbNKMpWwY3IjlDmdQNZEKID/fsGbHNe9Hu1WLqsosR3R5Cmzaqq4QqOXJQek0goZo+LORblgdcgr17VVcJYfSiYlKXdblmugav9zcZOxls04qbG25EEokbREaqrhHCeOl0JHftSdurQViRzDKXCViGTFFdJVSrVYtvAovSnh/ozhwutRkJ9+6prhLCeD19StT9POTgPrlcc5vM3gUZbNPK68H2MsV5fuKc6hohjNcPPzBuRWHCcefHjF3ItX6erKsVAGiGDWV2ndXk5Tat788isW0n0GpVZwlhnGJi/rhxTONmGssQQAbbtOPqShnOoMOCmGPPVdcIYZyuXOFI9yWMxp/hjKX24m+gWDHVVcJQWFpiv2I+q7L34BxODNrdAKbI03whPkh0NNG4mMxVur+TwTatZMlC6UIvsCCFyEt28hRBiPeVnMwzr660S5hPRU7g1/EG/OlIQCEAyJ2bsquGM4VBhNKbzcMOwpkzqquEMDoJJ89xiRIy2Ip/Zle2FCW5SOTLknD9uuocIYxLcDB9j7fnHrlY/ukwrGdNU10kDJWHBz0HZaQRm+mcMpd7rXtDQoLqKiGMyrmjT9FiiQvRUKaM6pw0I4NtWvrzBjJ5giDEuzt0iFWB51mCN7M0fSm2ahzY26uuEgZMExzEQudp6NDw7UVfdEOGqk4SwnjodESdS73KwDn/E8ieXXFQ2pHBNi39abDVnZGTEYR4J0+f8msbX7rp5tCalXQMKApVq6quEobO1pbcq2aywLonm2jKgpkvYPt21VVCGIfr14lMLEERrmJfzoBuck0DMtimpTKpZ9k+Iju/Hb2lukYIo5DSsw8db43Bnmd8V2UJmuHDVCcJY1G6NE2n1KYL8/AhhEvtR8ODB6qrhDB8Z85wmrKU45RJLUMAGWzTVtGiuGW4BEDkadk8JsR/WrWKKctzc4Ba/JCxO5+smANW6XLTtzAVvXszte5W8hFL+wfTePVtD9DpVFcJYdB0p05zmrKU5TSULas6J03JYJuWLC0p6JqVLDwh8k5OiI9XXSSE4bp7l5iu0xlFEAOZQu25baFwYdVVwthoNGReOptlWXsTQQXGbHCGxYtVVwlh0H47eos4HCjDGXliK/6dpoybXK0rxH/R6XjVpSdfP5lOMa4Q9OVJaNdOdZUwVnnzUmVJD0YRRDAjOdpnGdy8qbpKCIN1+nTqf8vaXTS5s8JlsE1rr9fZytW6QvyLZcsYt8mZ05RlabZ+ZJg3w2SucxSKfPklIzreojwn+eZFKInesiRBiLd68oTTd/PwCXF8WiY7WJjWKGhavxtD8PpkhAuUIjEiRnWNEIbnt9841XM+QYxiGOOouLAH5MypukqYAKvpU1icayhXKYr/bndYsEB1khCGJzL1u8plOIOmrGktQwAZbNOeqytuRJKCFeeOPVNdI4Rh0el42akHX8fPwpkYRrW5BC1aqK4SpiJbNpyX+BJAAJMZxBGfFXDjhuoqIQzL6xMRynLa5NbXggy2ae+TT3DJ/xiAM+dt5VthQvzZokUE7qjCeRz5PscgbGaHqC4SpqZhQ3y9H1KBCLxfzJIlCUL8xbPj57lMCZM8EQFksNUL+3LFKcFFTiU6ytW6Qvzuxg2O9l3OeIbiz2jKLPYxqdtuhOGwmjaJJblfL0nYWxvmzVOdJITBiDyaev10Wc6Aq6vimrQng60+lCtHOU5xinJw8qTqGiHU0+lI9O7BNy9CqUAEQzregcaNVVcJU5U1K6W/H8Jo/FOXJPRfKQ8ZhABITubMlcxYk4RTiWTIlEl1UZqTwVYfXg+2pymL9uRp1TVCqLd4McF7q3OVoizJPRSrGVNVFwlT5+nJoM6PqcgJvBNCSfy6G2jl4hxh5i5d4nSyM6U5i03Z0qpr9EIGW314Pdg+IwtXDt5RXSOEWrdvc8ZnMRMYwgjGUHqxL2TNqrpKmAGrqRNZnHcEVymK3wEPWLhQdZIQap0+bdIbx0AGW/0oVIhyWa8BcOq0nM0pzFtyr350fjYNR84z9H83oGFD1UnCXGTJ8seShCkMJGLAMoiNVV0lhDLJp6KIwjX1xjET3DgGMtjqh0ZDrgqfkp9bnHpSBO7IU1thptatI2RdQU5SngVZB2Ezc4rqImFu6tVjYIf7uBLFt/EhvOrlo7pICGUuHX5AInbyxFZ8gD9vIDt1SnWNEOnv8WOudJuIH4H0YzpVZn8NOXKorhJmyHraRBZk8yUSN6atLwzr1qlOEkKJ01GWAJTJ9ivkz6+4Rj9ksNWX8uX/GGx1EXIygjA/Ot/BdL0fTG7uElw/HLy8VCcJc+XgQMU5nfEhBH9Gc7nbJHj8WHWVEOnr9m1OPylMQW6QvUIRk73GXAZbfSlXjvKc5B65uX34uuoaIdLXvn0sWpDCHjyYm8GHTPNDTPYvUWEk2rQh0PMQebhD9/uB6IYMVV0kRPqKiPjjKl0qVFBdozcy2OpLyZKUy3AegJMRcuuNMCMJCdz2Hs5ApvA1S6g/uT4ULKi6Spg7jYZM86bxXYb+7KYu389LhAMHVFcJkW50ESc5RbnU9bXly6vO0RsZbPXF0pKCZbPzCXGcuptXvu0lzMfo0fS+PhBbXjK18kro0UN1kRCpChbEc6IH7VjGAKZy13soJCaqrhIiXdw6eIN75KYCEfLEVnwYTfk/bSA7LRc1CDMQFcWmSedZS0umWw4k+/fTwEL+mhEGpGdPplVYjgVafK72geBg1UVCpIuIiNT/VrC/BMWKqY3RI/mKo0+v19nK1brCLGi1PO/iQx9tCPXZTpuRJcDRUXWVEG+ytCTn95OZZunLCrzYMu4MxMSorhJCv+7eJSKuMLm4S/7yuU16z4MMtvr0+mSE6xTh0ZELqmuE0K9Fiwg6Wo875GF24UlohsnmHGGgnJ1pP+xT6rOdntqZPO/aH3SyF0KYsJMniaACFYhAU9F0lyGADLb65exMOcsoAE4fS1IcI4Qe3b9P9MDFTGEgIwmm2IJhYGurukqIf6QZMZzQghO5Qx6CD30G33+vOkkIvdFFpA62FTlh0hvHQAZb/bK1paSzNRl5zslfHeDFC9VFQuiF1ncI3Z9OoBhX8P3fLfDwUJ0kxL/LkIHi8wYzjHFMZhBn+8+Hhw9VVwmhF78dvG4WG8dABlu9s6xQljKc4ZSuLERFqc4RIu3t38/i7zUcpCbfZRqEbcgE1UVCvBtPT4a0vEJhrtPz8Rh0Q4epLhJCLyJOpC61qZDxPJQoobhGv2Sw1bfXG8giqCAbyITpSUrifpfhDGYiHfmeOpMbQ+7cqquEeGcZZkxktt0g9lOHZQsS4NAh1UlCpK0HD4h4UPD/N46Z+Ek1pv27MwQVUte0XKAUTw9Fq64RIm1NmYLvpS7o0DCp/Aro2lV1kRDvJ18+6o33oA0rGMgUHnUZDMnJqquESDtmtHEMZLDVv7JlqWRxEh0WnDyYoLpGiLRz9Sr7/PfyPd8wQTOMXAvHmfyTAGGievZkqssiEsnA8LPtYMYM1UVCpJnfN45VIMLkN46BDLb6lzEjjs6WZCKe49dzwPPnqouE+Hg6HUk9fejxajrVOUjnfpmhbFnVVUJ8GCsr8i0MJphRzKUbx0ZsgJs3VVcJkSZiD17jLnnMYuMYyGCbLiwrV6A8JzmhqwCnTqnOEeLjrV3L1O2luUQJvsvlj0VggOoiIT5O5cr07JZCGc7QPXEaKf0GqC4SIk2cOP5641iGs1CqlOIa/ZPBNj1UqkQljnOcSnD8uOoaIT7O8+fc7DORIEbRlxm4zukJ9vaqq4T4aFbjg/nuk+Gcpiyz1+WBLVtUJwnxcR49IuJeAXJyjwLlcoKlpeoivZPBNj1UqkRFTnCNojz45bzqGiE+zrhxDLo9AHueEVD3IDRvrrpIiLSRLRtVZrSjK/MYSTC3ewZBYqLqKiE+nJltHAMZbNOHqyuVbM4AEHHkleIYIT7C5cvsmXCcVbRhkuUwssweb9J3jgsz1K4d46pvxoYkhtzoAVOmqC4S4oPpjh3//41jZrC+FmSwTR/W1hQrm4VsPOJ4bD54/Fh1kRAf5FXfgfRJnkoNfqG9b16TP+hbmCGNhk++G8c4zQh+oCMHA3fDr7+qrhLig8QeuMxd8qRepVuliuqcdCGDbTrRVKpIRU6krrM9cUJ1jhDvb/NmZm4txnkcmZUrCM3IEaqLhNAPV1c69c5IRY7TO2kKKQN9VRcJ8UFOHNMCUCHzRShZUnFN+pDBNr3IBjJhzBITud0rmAAC6MEcys7sDJkyqa4SQm8sAgOYlW0UpynH/DWfwJ49qpOEeD+//caxuGLk4TYFKuczm3PGzeN3aQheD7a3yUfsgcuqa4R4P5MnM+TXntjyksCaO6FVK9VFQuhXtmxUmdIabxYxgjE87DESXskeCWFEjh3jKFWowlE0VSqrrkk3Mtiml1KlqJjxHCAPbIWRuXGDX4L28gMdGacZQfbvxsqGMWEevvmGceVWk4wVIy92gFmzVBcJ8c60R1O/S1yFo2azvhZksE0/lpYUqJiH3Nzh+MMicOeO6iIh3klyf196J02hEsfo1M8enJ1VJwmRPiwsyD03kNEEMJdunBq1Vv7uFkbj/L47PCUrlTkGleWJrdADTeXU5QgnqCgbyIRxpiUlNgAAIABJREFU2LmTuetyEokbodn9sAjwU10kRPqqVIlenRIozVl6Px+PbshQ1UVC/LeUFI6dsUGDlkr5YiFvXtVF6UYG2/T0+qKG41RCd0zWIwgDl5TE/Z7+jCSYziyk0lQvyJpVdZUQ6c56fBAzMw3jEDVYtjQFDh1SnSTEv7twgaOJZXDiHFmqllZdk65ksE1PFStSiePE4cC1AzdV1wjx72bMYNjlTmjQMa7iWujQQXWREGrkzMlnExrQilUMZiJPew6FlBTVVUL8sz9tHDOn9bUgg236KlKEyp+knohwJMIadDrFQUL8g9hYIvw2sJBvCcKPHHPHmM1RMUK8VbduTHZcyFOyEHSmCSxYoLpIiH+UcOgUkbiZ3fpakME2fWk05KhWghJc5Ei8M1yWY7+EYdING45PwlhciKJbVx2UL686SQi1rKwoOG8kwxlLCD6cH7IY4uJUVwnxVicPxJOCFVU4ZjZX6f5OBtv0Vq0aVTnCYarB4cOqa4T4u4gI1ix9zi+4MzWTH1ZjA1UXCWEY3N0Z2OY3CnGDfk9GowsYrbpIiL9LSODopezY8QKX0lqwt1ddlK5ksE1v1apRjcOcpiwJv0SorhHiTTodiX0HM5iJNGEj9cbUAQcH1VVCGIwMk4OZYjOcHXjyc+h1OHdOdZIQbzp9mqPaipTnJNZVzO+7bTLYprfKlanKUZKxJmLvU9U1QrxpzRpCDlXiFgWYXGQ29OypukgIw1KgAE1HuOLBLgZoJ5HUz1f2SwjDcuwYx6hslhvHIB0G2/j4eHx8fMifPz92dnaUK1eOlStXvtN77927xzfffEPOnDnJlCkT1atXZ4+x39dtb4+ri46MPOfIlRwQH6+6SIhUiYncGTCRMYygN7MoObMPWFurrhLC4Gh8BzEtz0SuUIxZO0vCli2qk4T4w73957hOkdTB1sw2jkE6DLYtWrRg6dKlBAQEsG3bNipVqoSXlxdhYWH/+r6XL1/i4eHB3r17mTFjBhs3biR37tw0aNCAAwcO6Dtbr6yqV6YSxzmsqyr36wrDMW0aI291w4Yk/OqEwxdfqC4SwjDZ2eEa0pluzCUQP+73DYKkJNVVQgBw9FAyAFVsz4CLi+Ka9KfXwXbLli3s2rWLOXPm0KVLF2rXrs28efOoV68evr6+aLXaf3zvwoULiYmJYdWqVXh5eeHh4cGaNWsoWbIkgwcP1me2/r1eZ3uYaugOH1FdIwTcucOpoM0sohOjNQF8MjMQNBrVVUIYrtatCayyBQ06Rl39BkJDVRcJAbdvc+xuIXJxl4JV8prld930OtiuW7cOe3t7WrVq9cbr3t7exMbGcvTo0X99r6OjI1X+tD7E0tKS9u3bc+zYMW7fvq23br17PdjeJh8391xSXSMEuhEj6Z8wBkfO062Lziz/lS/Ee9FoyBE6Gn8CmU8XIv3WwP37qquEuTt8mCNUpQpH0dSorrpGCb0OttHR0Tg5OWHxl4PdXV1dAYiJifnX97q5uf3t9Xd5r8ErUYKqWc8DcPiYpWw8EGqdOsX6RXHspw5TM47COthfdZEQxqFCBXp985wSXKJ/fCC6UX6qi4SZSw5PHWxrcBCqy2Cb5h4+fEj27Nn/9vrvrz18+PAf3xsXF/fB7zV4Fhbkql6colzhyLPScPWq6iJhrnQ6XvYbzCAm0YCtNAisDjlzqq4SwmhYjwtkqt1I9uDBhnl3ITJSdZIwY1G77xGPPdU5BFWrqs5RQo77UkUuahCGYP16ZoaX4QaFmFJwBvTpo7pICOOSJw9f+FeiAVsZpJvIy75y/JdQJDGRQzFZsSaJiiWeQo4cqouUsNLnB3dwcHjrk9W419cQOvzLwe8ODg5//Lr3fS+Aj48P2bJle+M1Ly8vvLy8/rM7XVSrRjXWs4avePlLGLbt26suEubm5Uvu9R9HEDvpwRxKz+gONjaqq4QwPj4+TJ3VBNdbW5i+vwyD16+H5s1VVwlzExHBQW1VKhCBXU3DukY3LCzsb6dhPX78WC+fS6+DrZubG2FhYWi12jfW2UZFRQHg8i8bVFxdXYl8y7d03uW9ACEhIZQ35PvtK1emGkNJwpaTe59QTXWPMD8zZuB3oxMWaAlw3wNN16ouEsI42driNKMHvVqEEsxIvvapS+4vvgBbW9VlwpwcOsRBWvMVawxufe3bHiyePHmSChXSfgDX61KE5s2bEx8fz5o1a954fcmSJeTPn/+NEw/e9t7z589z7NixP15LTk5m2bJlVK1alTx58uitO11kyYKbcwp2vODIJQd4/lx1kTAn9+4RFfAT8+mCP4E4zAyQ472E+Bhffol/zT1Y84oRv3aFkBDVRcLM3Np9gV8pZNYbx0DPg22DBg2oV68ePXr0YMGCBezdu5euXbuyY8cOJk6ciOb1F9LOnTtjbW3NzZs3/3hvp06dcHZ2plWrVoSFhbFr1y5at27NpUuXmDBhgj6z04119UpU4jgHddXkogaRrnSj/Oj/IpjiXKZnp0QoU0Z1khDGTaMhe2gQgRp/FtGJU6M3wp07qquEudDpOHQ4daaqniUGHB0VB6mj981ja9eupUOHDvj5+dGwYUOOHz/OihUr3ngkrdVq0Wq16P604N7Gxobdu3fz2Wef0adPH5o2bcrdu3fZunUr7u7u+s5OH9WrU5NfCMcdXfgvqmuEuYiMZPP82+ymLlPsRmEzNkB1kRCmwc2Nbl10lOYs/RLGoRs2XHWRMBdXrnDwqQtFuUKeGsXAwnzPBtD77zxTpkyEhIQQGxtLYmIip06donXr1m/8msWLF5OSkkLBggXfeD1XrlwsWbKEBw8e8OLFCw4ePMjnn3+u7+T0U6sW7oRzj9xc2i5Hfol0oNOR5DOYgbpJ1GUnjfwrQu7cqquEMBlWwQFMyziScGqxZkk8RESoThLm4NAhDlHd7JchgBz3pVaRIlTPcw0LUgg/YQfJyaqLhKnbtInZex25QjGm5p+Kxqef6iIhTEvOnNQLrk0TNuLLRBL7Dpbjv4TePd9/glOUk8EWGWzV0mjIUrscZThD+MtKcPq06iJhypKSeNB/DKPxpyvzcJ3+rezaFkIfevViSuFZxJKPaYcqw182UAuR1o7tiScFK6prjkDlyqpzlJLBVjV39z/W2RIerrpGmLLQUAKudkCLBYFVt0KLFqqLhDBNNjaUmNmXPsxkLMO53X8iJCaqrhKmKi6OQ9fzkoUnOJe1hsyZVRcpJYOtau7uuBPOVYoRuyNadY0wVQ8ecNZvBd/RnVEEkzM0QI73EkKfGjViVJ1fsOUlI3/rDtOmqS4SpuqXXzhIDapxGIs6tVTXKCeDrWouLrhnSb2I4peDGlmLJfTD35+B8QEU5jp9Oj4BQ768RAhToNGQbWbq8V+L8eZk0M9w+7bqKmGCUvYe4CA1cCccaslgK4OtahYW5KlVkuJcIvxZGbhwQXWRMDUxMWz97gbbaMhk25HYjh+tukgI8+DiQtduFjhxjv4JY9CNGKm6SJig09vv8pSs1GY/1KypOkc5GWwNwevlCOG4w4EDqmuEKdHpeOXjywDtJOqwl2aj3CBvXtVVQpgNqyB/pmUcyQFqs3bxEzh5UnWSMCXPnrH/fC4ykEAlp+eQI4fqIuVksDUErwfbSNx4vFvOPBRpaOtW5u4qygVKMS3vJDQD+qsuEsK85MhB/eBaNGKzHP8l0t6hQxzQuVONw9h+Zt7HfP1OBltDUKECNW1PoMOCQ/teqq4RpuLVK+L6jcaf0XRmIWWnfQ12dqqrhDA/vXoxudAsbvIp0w9WgJ9+Ul0kTIR2f+p3e2uzX9bXviaDrSGwsaF41Rzk5g7h90rBzZuqi4QpmDOHwMteJGFDUMWN8Jcb/4QQ6cTGBseZvehFKGMYwd0BE+T4L5EmorfdIg4HanEA3N1V5xgEGWwNhKbWn9bZynm24mPFxXFh1DJC6cUIxpBntp8c7yWESo0b41f7ANa8YtTNLhASorpIGLuEBPZHfoINL6la9D7ky6e6yCDIYGsoatXCnXCOU4mEPYdV1whjN3o0g56OogC38Gn3ACpVUl0khHnTaMg+czSjNQEs4FtOB26EO3dUVwljdvQo+1NqUJlj2NWporrGYMhgayiqVuUzy3CSsOXQ9meqa4QxO3eOHbMuspkmTLIZSYYJcryXEAbB1ZVuXXQ4cl6O/xIfTbf/AAeolbq+tnZt1TkGQwZbQ5E5M86VM5GD++y9VRxiY1UXCSOVPGAwA7STcOcALYeXgvz5VScJIV6zDvZnasZR7OMzNix6CKdOqU4SRur8tuvcJ5dsHPsLGWwNiIXHZ3zGXvbwOezdqzpHGKPt21mwLT9nKc20XOPR+A5SXSSE+LOcOWkQVIOGbGEQk3jZ11eO/xLv7+VL9kdkxopXVM93AwoVUl1kMGSwNSSff87n7OE4lXi2/ZDqGmFskpN53M+fUQTxNd9TYVp7yJhRdZUQ4q9692ZKwRlcpzAzfykLa9eqLhLG5sgR9r+qRkVOkKluNdkc/Ccy2BqSatX4zPogyVjzy/bnqmuEsZk3j+ALX5GAHWPK/QReXqqLhBBvY2OD08ye9GQ2QYzi3oDxcvyXeC+63XvYT+3UZQiff646x6DIYGtIMmSgZI2c5OM39t4rDdeuqS4SxuLRIy6NWMIM+jKU8eSbPVL+BS+EIWvSBP9a+7AkBb9fO8P06aqLhBE5t/kKt8nH5+wBDw/VOQZFBlsDo/H4XNbZivcXHMzgx8PIwx0Gtr4FVauqLhJC/BuNBodZownQjGY+XYgMXC/Hf4l3Ex/PrjM5seElNYvfhQIFVBcZFBlsDc3rdbanKMejrUdU1whjcPEie0IiWU9zJlqPxG5ykOoiIcS7cHWlx7fJlOASA14EoRs5SnWRMAYHDrBbW4fqHCJjvRqqawyODLaGplIlPrc7ghZLDux+JbtlxX9KGTSE/trJVOMQbYYUhk8/VZ0khHhH1mMCmGo3kt3UZdPCe3D6tOokYeCSd+1jH3Woyy5ZhvAWMtgaGmtrCtcpTGGusfdRGbhwQXWRMGS7d7N4kwORlGFajrFohgxWXSSEeB85c9IwsBqebGMgk0nqO0geaIh/dWLzHZ6SFQ/2QJ06qnMMjgy2huizz/icPanrbPfsUV0jDFVKCk/7jmQEY2jPD1SZ0hoyZ1ZdJYR4T5q+fZjy6XSuUYRZ4W6wbp3qJGGoHjxg16WCZOEJFcsmg4OD6iKDI4OtIXq9zjYKN+5tOaG6RhiqhQsZe7YZz7BnnNsKaN9edZEQ4kPY2OA8szvdmEsgftzvPxZevlRdJQzR3r3sxoM67MOqbh3VNQZJBltDVLYsHlkjANi1zxJSUhQHCYPz5AlXh81nGv0ZzEQKhA4DC/nfWQij1bQpo2vuAsD/105y/Jd4qxfbwzlEdTzYLetr/4F8JTRElpbkqedKGU6z43kNOH5cdZEwNGPGMCRuMDm5j2+Lq1CzpuoiIcTH0GjIEToaf00gc+lG9Oif4O5d1VXCwPyy9RlJ2FLXch+4u6vOMUgy2BoqT0/qs4Md1Ee3bbvqGmFILl/mwNQTrKEV46z8yDQlUHWRECItuLnRq/NLinFFjv8Sf3fjBrtinchLLE5Vs0KmTKqLDJIMtobK0xNPtnObfEStu6y6RhiQlEFD8EmZTGWO0s43HxQurDpJCJFGbMb4M8VuFDupz5YFsXL8l/h/27ezGw882I2mgafqGoMlg62h+vRTapR6iB0v2BGVBx4/Vl0kDMGePXy/ISunKE9I9iAshg9VXSSESEu5ctF4dCXqspMBTOFVPzn+S6S6v+EQpyiXen5tgwaqcwyWDLYGLEODOtRhH9t19WH3btU5QrWUFJ72GcFwxuLFj1Sb8pUc7yWECdL07cPUAtO4THFmH3CG9etVJwnVXr1i+x5rdFjgmf0ElC+vushgyWBryF6vsw3HnRc/71VdI1RbuJBxZ5vylCxMcF0OHTuqLhJC6IOtLa4zu9KVeQQQwMP+wXL8l7k7fJhtibUpx0nyNCwnp+D8C/mTMWS1a+NpvZeXZODA5qfy7Shz9vp4r6kMwJdJfDpbjvcSwqQ1a0ZgjR1osSDgxjcwY4bqIqGQdut2tuNJQ7aCp6yv/TfyldGQZcyIY61cFOAmO+6Xlet1zVlwMIPjhpCT+wxuKcd7CWHyNBpyhgbgpwliDj04G7BKjv8yYyfW/soDctKAbVC/vuocgyaDrYHTeNbHk9R/qbFjh+ococKlS+yfdpKf+IrxVqPkeC8hzEWZMvTp9IIiXGPgi0Dw81NdJFS4c4dtF4uQlcdUK/cScudWXWTQZLA1dK/X2Z7FmZvrI1TXCAVSBg2hX8oUqnCEtoMLQKFCqpOEEOnEZmwAkzOMYhsN2Tr/Fpw5ozpJpLcdO9hKQ+qxE6uG9VTXGDwZbA2dqyv1ckdhSTJbf7GHhATVRSI97d7N4o3ZOUNZQhyCsRg2RHWRECI95cpF09EV+JzdDNBNluO/zNDDDb9wlCqp62vlmK//JIOtodNo+KRxDWpwkM2v6sOePaqLRHpJTuZp35GMYAztWEbVKa3keC8hzJCmX1+m5Z/CRUry3X5H2LBBdZJILykp7NyhSz3mK/MhqFpVdZHBk8HWGDRpQmM2s4u6JKzbprpGpJeFCxlz9kueYc84txXQoYPqIiGECra2uM3swrcswJ/RxPUPkuO/zMXx42yNr4kbZ8hfrzRYW6suMngy2BqDunVpbL2DBDKyZ/0T+TaUOXj8mCvDFhCCD0OYIMd7CWHuvvySoOrbSMaK0dc7wsyZqotEOtBu2MQ2GqQuQ2jcWHWOUZCvlMYgUyYcPfJTjMtsflhNNg+Yg+BgfB8NIxf38G15DWrUUF0khFBJoyFXqD8jGUMovTgXsBLu3VNdJfTsxKqr3CM3DdkGjRqpzjEKMtgaCU2TxjRmM5tpjG7jJtU5Qp8uXWJvyBnW0YIJ1qPIOCVIdZEQwhCULUu/Ts8oxA0GPQ+Q479M3bVrbLzqTHYeUqNKshzz9Y5ksDUWjRrRmM3c4lMiV8tFDaYsZeBgfFImU5XDeA3+VI73EkL8wXasP5My+LGFRmyfdwMiI1UnCX3ZtIkNNKMRP2PVTJ7WvisZbI1FoULUcnlEZp6xOboQ3Lmjukjow7ZtLNyUk0jKEOIQjGaoHO8lhPiT3LlpHlCG2uxjgG4yyT5y/JepurryONG40owN0KSJ6hyjIYOtEbFp2gBPtrOJJrBli+ockdaSkojr7cdwxtKBpVQJ8ZLjvYQQf6Px6ce0/FM4hxPz9haHjRtVJ4m09uQJG4/kwoaX1C90EZydVRcZDRlsjcnrY7+OUZm7a8JV14i0Nn06flc6koQNEyr9BO3aqS4SQhgiW1vKzfCmE4vwI5BHPqPl+C9Ts20bG7WN8GA39l96gEajushoyGBrTCpV4guHYwBs3m0nt5CZkthYzvivZQ498COIvHMD5C8yIcQ/a96c4GpbeIktQdfbw6xZqotEGopbvZsD1KIpG2UZwnuSwdaYWFqSq1k13AlnbVIj2LFDdZFII7rBQ+iTMIESXKJvlwQoV051khDCkGk05AkdxQjGMpM+XPRfDvfvq64SaeHVK7ZuhRSsaGK/H2rVUl1kVGSwNTYtW9KCteyiLk/CZJ2tSTh4kBXLkwmnFjPsR2AzbrTqIiGEMShXDp+vH1GAW3L8lykJD2fDi7pU5Dj5G5WV28bekwy2xsbDgxaZd5KELT9vTIGkJNVF4mOkpBDfw5dBTOZL1lF/Ql1wcFBdJYQwEhkmjGZSBj820ZSdc69CVJTqJPGRXq7awFYapp6G0KyZ6hyjI4OtsbG15dOm5ajEMdYmNIA9e1QXiY8xfz5jopoQR3amOs6Hrl1VFwkhjEnu3LT0d6EW++mnm0ZSXzn+y6ilpLB71QPisaepzXa5bewDyGBrjFq2pCU/sZWGvFi1WXWN+FAPH3Jp6EKmMJDBTKTIvGFgaam6SghhZDQ+/ZiZfwIXKcn0fW6werXqJPGhDh1i9aO6lOQCrg0LgL296iKjI4OtMWrQgBYZtvKCTGz/KR5SUlQXiQ/h54fPkwDycpshra+Du7vqIiGEMcqQAbfZ3enDTEbjz62+E+HZM9VV4gMkrVzHer6kFavRtPpKdY5RksHWGGXMSIlGJXElkp+e1oVwOdPW6Jw5w+Y5N9lCI6baDifj1GDVRUIIY9a0KQGeR8hMPAPv+kJgoOoi8b60WnavuM9jPqGV1Xpo3Fh1kVGSwdZYtWhBC9ayiSYkrVqvuka8D52OxF4D8dFNpS47aRHgBvnzq64SQhi5rLPHMdl6GKtow66pkRATozpJvI+jR1n98DNKcBG3+nkga1bVRUZJBltj1bgxLa028pSs7F75ALRa1UXiXS1dytSDlblBIWYUmoqmv4/qIiGEKShalHYjiuDOAfpoQ0jq0U82khmRV6tkGUJa0PtgGx8fj4+PD/nz58fOzo5y5cqxcuXKd3rvkiVLsLCweOuPe/fu6bncwGXJgkv9fJTkAivi6sHBg6qLxLuIi+N6/+kEM5J+TMdprg/Y2qquEkKYCM2QwYQWGM8lShASXh5+/FF1kngXOh17frzDI7LTynIdNG2qusho6X2wbdGiBUuXLiUgIIBt27ZRqVIlvLy8CAsLe+ePsWTJEo4cOfLGj+zZs+ux2jho2nrRlh9ZR3MSlsouWKMwfDh9HwWQnTj8v4wET0/VRUIIU5IhA67z+tCHmQTix83+U+HJE9VV4r8cOcLqe7UoziXKeOQAmXE+mF4H2y1btrBr1y7mzJlDly5dqF27NvPmzaNevXr4+vqifcdvn7u4uFC5cuU3flhZWekz3Tg0a4aX7TqekYWfVzyVyxoM3dGjbJh7h000ZXqGIdjPGqe6SAhhiho2JKDRCex5xsD7Q8DfX3WR+A9JS1ewjuapyxDaeqnOMWp6HWzXrVuHvb09rVq1euN1b29vYmNjOXr06Dt9HJ2sEXq7zJkp2dyZihznx/imsGOH6iLxT5KTed7Fh75MpyFbaDGmomwYE0LoTdbQsUy2Hs5qWrNzxjk4c0Z1kvgnr16x/ceHxOGAl81aaN5cdZFR0+tgGx0djZOTExYWb34aV1dXAGLeccdm48aNsbKywsHBgZYtW77z+8xC27a05Ud+phGPl8jpCAZr9myCoppxj1zMKjULTd8+qouEEKasUCHa+pegFvvprZvByx4+ssnYUO3axbKnTXDjDK7NikKWLKqLjJpeB9uHDx++dS3s7689fPjwX9+fN29eRo4cycKFC9m3bx9BQUEcP36cqlWrEiX3Yafy9KRN1u28wpp1m6wgPl51kfir2Fhihi1jCgMZzliKLhoJspRGCKFnmkEDCS00iSsUY/Lh6rB0qeok8RZPl6xlI01pzzJo1051jtF758F23759/3hCwV9/REZGpkmcp6cngYGBfPHFF9SsWZOePXsSHh6ORqPBz88vTT6H0bOxIV8bd+qwjx+TWsLGjaqLxF/o+g+g54tJFOUqg70fQPXqqpOEEObA1haX+f0YwFSCGMXl/qHw4IHqKvFnz5+zdr0FL7HFK8sWaNBAdZHRe+fHRo6OjixYsOCdfm3BggUBcHBweOtT2bi4uD9+/n0VKlSIGjVqcOTIkfd+r8lq146285bSjbncXvQNedu2VV0kfrdzJz+ssuEAtdmVpTm2k97t/yEhhEgT9erh33IZq366Q4/HY9kxYCCapd+rrhK/27SJZUmtqMM+CrSpIcc/poF3Hmzz5MlDp06d3uuDu7m5ERYWhlarfWOd7e/LCFxcXN7r4/2ZRqP515/38fEhW7Zsb7zm5eWFl5cJ7jasWZOW+XrSMzaFlXtz4XP/PuTMqbpKJCYS1304g/gZL37EY1oT+IB/zAkhxMfINGsCc7b14Yvnq1n+w/e0/3o3eHiozhLAbwu2sofFLOBbaNtRdY7ehIWF/e2Y18ePH+vnk+n0aOvWrTqNRqNbuXLlG697enrqChQooNNqte/9Ma9cuaLLlCmTrkWLFm/9+YiICB2gi4iI+KBmo+Xrq2vBGp0bp3XaqdNU1widTqcbOVLXhbm6LDzWxVZqqtOlpKguEkKYq3nzdG0I0+Xgnu5B4Qo63YsXqotEbKxuksZXZ0uC7nE+J7P7GqGveU2vm8caNGhAvXr16NGjBwsWLGDv3r107dqVHTt2MHHixDeeunbu3Blra2tu3rz5x2v16tVj3LhxbNy4kT179jB9+nTc3d2xsrIiKChIn+nGx9ubTiwikjKcmn1YrlFULSqK/eMOMp+ujLcYQd5FY8BCbrAWQijSuTMhlcNIxgrf6z1BvoYqp1v6A0t17WnCJrJ6t5CvEWlE73+Ka9eupUOHDvj5+dGwYUOOHz/OihUr/rYkQKvVotVq3ziz1tXVleXLl9OhQwcaNGjApEmTqFu3LidOnKB06dL6TjcuTk54Vn1KXmJZfLkmRESoLjJfKSkkdu5Fl5TvqMEvdBvuAB+x7EYIIT6ahQV5loxnguVwFtOJfROOgpwupI5Ox4nZx4jCjU4sAm9v1UUmQ6PTmdajvZMnT1KhQgUiIiIoX7686pz0tWABQ7s8YB5die0SQIZ5M1QXmacZMxjR7xmTGcTpws1xOrcWMmRQXSWEEGj9R1Mr0IMH5OBMpS7YHt4Hlpaqs8zPoUN0qxHFFr7geu1vsNy3W3VRutPXvCbPvU1J69Z4Z1jBI7KzcdlTSEhQXWR+btwgcuhyJjKYEYzBaekwGWqFEAbDYvhQ5hYez1WKMu64B3z3neoks/R87jLC8MKbxVh2/kZxjWmRwdaUZMlCqTZlqc5BFiX8D9atU11kXnQ6Unr05tuEmZTiAkO7xIG7u+oqIYT4f7a2OC/xZTATGccwzg1eDL/9prrKvMTHs3qllngy0ynTKmjZUnWRSZHB1tR06kQnFrGD+twMlcst1jWTAAAcJ0lEQVQa0tWKFczcWowTVGR+juHYTBqjukgIIf6udm1Gfn2Lwlyn04uZpHTrKRuO09Pq1Sx82Q4PdlO4XQ3ImFF1kUmRwdbUuLvTumgEmXjOgkNOcPGi6iLz8PAh13tPZiTB9CKUavM7QdasqquEEOKtMkwdy6JPBnGUKkz/uRgsX646yWycn7aVX3DnWxbIpjE9kMHW1Gg02HdvRwd+YB5deRU6T3WRWdD16Uv3uDF8wiPGNj0KX36pOkkIIf5Z9uzUWNSZfkxnBGO41CsEbt9WXWX6TpxgQVRlsvOQL12vQpUqqotMjgy2pqhTJ3raLOQOeVm34CG8eKG6yLStXcvCMDu204C5mQZi/90k1UVCCPHfvvySMa3OkJ/f6PR0GtpuPWRJgp69mLGARaQuGbTt0xX+r717j++xbvw4/vpucxjNaRvblOggbpuImKmUcw4xDGMZQxkawlLJnZLNccMMbQ5JY0zIqYMzSW4kh3TnTo5z3CxmGbPr94c79+2Xuzb7btf2/b6fj4fHI5++313vqz4+e7t2fT/XXzxFVXJPxdYWubri3aM2z7GV2IxgSEgwO5HtuniRE/3H8TpTCWEubWa/BJ6eZqcSEcmRUrGTmVtuBDt4lpjVVeD/PfZUrCglhYTFFtIox0CXRdCjh9mJbJKKra0aNIiBxLKV5zk8eb3+Fp5PjIGD6Js6kXKkMbXNRujZ0+xIIiI55+ZGk7m9GEQMbxLBzwOnwLlzZqeySca8+cRkvUo71lAt5AUoXdrsSDZJxdZW1a+Pf/3TVOIcs/75AnzzjdmJbM/SpcxJqsBGmjPXZRhl50Xpx0oiUvR06kRk5z1U5AL9fp1M9gDtkmB12dl8Hb2b76nDYGJg4ECzE9ksFVsbVvy1V+lPHB8RzJWouWbHsS3nz3Ps1QmMYDKvMpsWcV2hUiWzU4mI3JcHZk8mvuwItvACs1Z56pYEa/v8c2Ykd+YJfqR5cwtUr252IpulYmvLunYltHwimZQgbnkFOHnS7ES2wTDIfjWUkLQpuHGJSR2+hq5dzU4lInL/3NxoNrcHocQykkn8MzQaTp82O5XNSP5gPp/SiUHMxGGwrtbmJxVbW1ayJF5hXehBAtFGGDenTDc7kW34+GNmrqrMVp5nftlhuMRN1S0IIlL0de7MpC67eZDTBF2Zyc1efSE72+xURd/evczY+RTO/EbwIzugXTuzE9k0FVtbN2gQw4vHcJqHWDrnMqSlmZ2oaDt2jB9CZxDORAYzgxfmBoG7u9mpRESsovSHUSyqOJzvqMu4zX4QHW12pCLvyvgYZhHKAGZTJnwAODqaHcmmqdjaOnd3fPo2oDXrmZw5GGP2HLMTFV1ZWWT26EOPjDiq8QsTex7QM75FxLaUL0+DxcMYw3t8wNvsemMFHDhgdqqi6/hx5qxwJ4NSDKmwCHr1MjuRzVOxtQfDhjGCKeynLhsn7YPMTLMTFU3jxvHWty9xhJokVA7HedZUsxOJiFhf06a8New6T/MPgrLmk969H1y/bnaqIilz0nSijKG8zMdUHtIFnJ3NjmTzVGztweOP09S/LHX4jkmpIXom+P3YuZOv3vuGqQwnwvI2dZJGg4uL2alERPKFU8T7fPzEB5zDg9eP9IM33zQ7UtGTksKi+OucxYuRJWZoi68ComJrJyzhIxlFJF/Sim/HrIWsLLMjFR1XrnCp+2CCjfm04EuG/r0s+PqanUpEJP+UKMFjSZFEOYUTxyusiD4OX3xhdqoiJXtqNJNuhNGBldTo9wy4uZkdyS6o2NoLX1+6NEujJj8w9kxf+PhjsxMVDYaBMSCUfqfGcIPiLKg/E4e3deVCROyAtzf9Jj1BZ5IIYR7HA9+EM2fMTlU0pKSQNPUk/6QGoxwnQ3i42YnshoqtHXEcO4YxvMd62vDt6NVw86bZkQq/+HhiFldgFR2Jdw7Da9k0cHIyO5WISIGwhL1GfIullCONbpdncaN7L/3ELwduTYnm3etv0Jr1+Pb3gSpVzI5kN1Rs7UnjxgQ0//dV2+R+umr7V/bvZ/egjxjOFIYQTceP/KFqVbNTiYgUHAcHyiXEklgxjO+oy1s7XoR33zU7VeGWkkJiVDJH+BtjHd/X/ckFTMXWzuiqbQ5duUJqp350vbmIunzHxNDjEBBgdioRkYLn5kaDT0cxwfImUxjBmg++1/22fyJrcjRjr4fTjtU06P+krtYWMBVbe+PnR0CLX6nJD7xzNhTi481OVPgYBkbffvT+ZQxXKMNSn3EUj5pgdioREfM0bszQSA/asZpgFnAqMFz3295LcjKfTDnHTzzBWCddrTWDiq0dcnz/XcbzFl/Rki/e3AJXr5odqXCJiWFKUhVW8xILS4Xy8KrpUKKE2alERExlGTGcBS0SKM01Ai7PITMgCG7cMDtWofLb2+MYc3M0nVjOU4Mb62qtCVRs7VHDhnToUpxn2M7IX9/mVuQksxMVHtu2sXXYCkYRSTgTaPdJIFSrZnYqERHzOTjgujiG5RUHsp86vPZNIAwdanaqwuPwYaZ9VJZkvIh4YDyMHm12IrukYmunLJERTHYcxUFqs3DSeTh92uxI5jt5khP+Q+lyK5Hn2Ma4YanQsaPZqURECg9XV55e83dmOYURxyvMmXUL4uLMTlUoXBw2nvHGmwwkluqju4Krq9mR7JKKrb169FEahjWkK4mMvjmGa+FjzU5krowMMl7qjn9qPA+QztLnZ1Fs0nizU4mIFD5PP02fOD8GMpPXmMHO0I9h506zU5nriy949ys/HMjmHc+5EBZmdiK7pWJrz0aPJsIlghRc+WBxNdixw+xE5vj3h8X6fv8a/+QJVj74Gm7L54Cjo9nJREQKp969iRr0Mw35ls63EknuONB+P0yWmcmhV6Yzh1d5mw9wmzASnJ3NTmW3VGztWYUKPBLRn1FEMpkR/Bgy0T63/4qIYNKSB1lCIPNLhPLk+kioUMHsVCIihVrxqAks84vGkVt0vPghGe27wbVrZscqcNmTphB6chSP8S/C/PZCUJDZkeyaiq29GzCAN+p8yUOcYvDRMIxp081OVLASElj+9l5GEclbfEDXxf7g7W12KhGRwq9YMTxWzeEzj1c5TC2CvnudW4FBcOuW2ckKzvHjfPTeCXbwLLEOr1FiVjRYLGansmsqtvbO0RHnOdHMIIyNNCdx9EE4edLsVAVj61a+CZ5NEIvoRiLvj3MAf3+zU4mIFB1ubjz1ZSRLnENYSUfeWN349k4JhmF2svxnGKT0e4ORNz+gJ4toGuYNtWubncruqdgKNGhAmwFV6EwSYZkTudhruO0vSkeOcLT967TP+pT67GF+n+04vDXK7FQiIkWPjw/tV79CtMNwpjCCWTFZEBVldqr8t2ABwze+SBZOTKk4UY8aLiRUbOW2iAhmeowjGwcGbQ2AOXPMTpR/kpO51Konba4uwY1LrGo6nZIfTtePj0RE7lezZoTNq0MY0xhMDKuHb4GlS81OlX9On2bVoC/5iN5EMYxK8yKgbFmzUwkqtvK7cuWotGACsQxkGV1ZOnQnHDtmdirru3iRKy90oM2p2fxKWdbVHEGFlfPAycnsZCIiRVtwMFPfSaMjKwlgKVsC58DatWansj7D4GLwCF75LZr2fEbvl7OhbVuzU8m/qdjKf7RqRddXyhPAUgZmTuVMwFDb2iUhLY2MFh1o99MUfqI6X3iF8MiGD8HFxexkIiI2wXHsGBL6bOA5ttE+eyV7/D+AzZvNjmVVRlQ0oZu6cAtHPqw0Bsu0aLMjyX9RsZW7TZ5M7EORlOQ6gftGkDXKRh4JeO0amS92pNP3Y9jHU6xzC6bujhng5WV2MhER22GxUCIuhk87J+DNIVrf/Iwf2o6EXbvMTmYdu3fz4cijLKcLswjFY0EklC9vdir5Lyq2cjcXF9ySZrPEMYid+PH3qWVg1SqzU+XN1avcaNWewF1hbOF5VpXphd+2SKhWzexkIiK2x9GRBxbHsa7VdLxIpulva/iheVjRfzpZWhp7O75PWHYUocQSEP4ItG5tdir5f1Rs5Y8aNOCZqM58wNuM523W9fwEfvrJ7FT3Jy2N683a0uXroaylLUmlgmm25R2oWdPsZCIitqtYMcqvnM/GZ8dSifM8f20NB5oPg02bzE52f7KyuNylPwFnp1GbA0T5LoVx48xOJfegYiv3NngwIzv/Qns+I/BaHIebD4GUFLNT5c6lS2Q834aX/jGar2jBqgeCaLctHOrWNTuZiIjtK1kS988/ZlOT93iIU7zw2zr2vfg2rF9vdrJcuzlkBN039uMy5VlWtj8lEhdCsWJmx5J7ULGVe7NYcJgXzyd/G09VjtP21CzOt+kDmZlmJ8uZn38mrWErXvw+gp34sb5cD1p//Q7Uq2d2MhER+1GqFK6ff8LG1pN5jH/R9MZ6trafDIsWmZ0sx4yYmQyMrcUmmrLcsRtVV02DKlXMjiX/g4qt/G9lyuCyfimr3fuSSQk67H6L9O79ICvL7GR/btcuTjzdhcbHFnIQH75068nz30ToiTAiImYoWZJyny3kK/9Z1GMvLW+tY8nLa+C99wr/w4CWLGHCa6eJpz/x9KNpXCA0aWJ2KvkTKrby56pUocq62awuEcBhavHSyj5kBPYtvOV22TL2NBlOw8vr+Q1nvnn0Zfz2TIcaNcxOJiJiv4oVo8yyuazv9yndSCSQJUz6+1WMXsGF9yeBq1Yxq8d23iSCMYwleJQX9Oljdir5Cyq28tfq16f+ytGsc+rAtzTEP6kH11/uX7j2uL15E0aMYHHXT2ly40uqcpxdjUfwxJ5P4OGHzU4nIiKOjhT/MIaPJl5gNOMIZxLBi5qT0agZHD9udrq7JSUR1/lzBhozGUI07/ZPhvHjzU4lOaBiKznTujXPrnid1Y7+bOM52i4J4tc2gZCebnYySE7mRtPWhE2pQg8W05GVbOo5j4qblkC5cmanExGR31ksWEaO4P3lf+OTYr1ZTmcafTeTn2v7w+rVZqe7bc4cogJ28sqtWQwihqiee7HMnqXHrhcRKraSc+3a0XT5ID53as9e6vHchndI9usCp06Zk8cwICGBf9VsT5Md45jNAGIcwlg0I41SH8+B4sXNySUiIn+uUyd67B7KrgcDyKAU9a5uJvGlRTB4sHkXTLKyyH59BMMHpPM6U3mDSGb03odlwXxwUF0qKvR/SnKnQweabHiHHS5tSKUC9Q/OY7t3aMFv35KcTHZAN2J77uDJK9u4QEW2uXdh0I5ALIMH6W/WIiKFXZ06+BxazJ7279GKL+hOIt1mPsulvz0HGzcWbJZz50h9oTPto14gimHMYDCRI1OxzJsLTk4Fm0XyRMVWcq9JE7y/iWP3g515nKO8cGUlE9ps4VboYLhyJX+PnZkJEyZw8DF/mi0PZRCxBPMR33d+H98j86FRo/w9voiIWE/ZspRdtZDEGRdZUuxlNtAc71Pr+Kj5QrIDusGxY/l7fMOAJUv4pnowT+2Yxi58WevYgcGxtWDiRF0kKYJUbOX+1KqF54Ev2Ng2iuFM4U0i8Jv9Moce94eEBMjOtu7xbtyA+HhSqjdi8KjS1PltJ8l48aVLZ2IT3XggaQG4ulr3mCIikv8sFhg8mG5H3uWQb3+asJXefIRf0jB2PREMr78OZ85Y/7iHDpHewp+hgedofHU9Hpxjn3trXtw6CkJDrX88KRAqtnL/ypfHafUKJkSX5OsSzbiKC3UufMGAnldIrtUCliy5XUjz4uxZiIjg3CN+hPdP5eGT21hILyZa3uDgqzNp8cuH0LWrdc5HRETM8+ijeH6dRGLcVbaU68h1StIoazsvRrXg64d7QEgIfPtt3ve+3buXm4G9mF07luobY/mQV5jMCL4OmMbDP6yHxo2tcz5iChVbyRuLBYYModGReXzX4g0m8AbLCODRH9fQN/Aa+zzawNCht++XyknJNYzb277ExpLdug0bHuxN4FtVefjM18xmAEOYxs9N+jJ8XxDFZ0/XVVoREVvi4AD9+tHkxEL2jkpicbFenOIhnrm1lfrzBzLbdz6/Vn8axoyB7dtztu2kYcBPP8HkyVyq24KJ9ROpvmQsA40YmrKJw14teX2ZH45LF4ObW/6fo+Qri2EU9sd+5M6+ffuoV68ee/fu5amnnjI7jn0xDNi4kV/fGM/Mfb7MIpTTPERNfuAlPqOp03bqeWfi+uSD4OEB5cvffk9mJiQnc+v4KY7u+ZU9qdX4kpZ8TmsuUpGa/EBf5hLS8jTlxw4DX1+zz1RERArC+fNkT5vB2uk/E3ctkLW0xZFbPMt22rCORsX38WStLEr7PHL7+0qFCrff9+/vK5nHznDoHxnsTPsbq2nPFp7HgWy6kcjwsnOp/WZbGDIESpY09zztUH71NRVbsT7DgG3byIqZzfpPf2NF9kusoR0XqQhARc7jRTLluYwFgxsU5yyenKEy13EG4En204Z1vOS+i4b9a2Pp1xeqVTPzrERExCxXr0JSEmc+XMvKXZVYS1s28wLXccZCNp6cxYtkXLiKBYMMSnEWT5Lx4ibFKcYNnmcLHVhFN+8fcBvcHV5+GUqVMvvM7JaKbQ6p2BYyaWnw+edkr1nHzzvPs+cXV/7Fo5zFk8uUx4KBE1l4cA4vkvEpc4I6dR1wbV4X2rWDJ5/Up1JFROQ/zpyBDRvI/HwzP+xI5bvT7pygCmfxJJ0HAChBJp6c5SFO8ZTLv6jd2AXn5xtCp07w+OMmn4BA/vU1bc4m+atcOejeHYfu3XkcePzqVTh5Es6du701mMUCxYqBpzd4eUGlSiqyIiLyv1WuDMHBlAgOpi5QNyPj9rZgly5Baurt+3SLlQX3WlC1Kri76/uKHVGxlYLl4gK1at3+JSIiklelSoG3t9kppJDQrggiIiIiYhNUbEVERETEJuRrsU1PTyc8PJyWLVvi7u6Og4MDY8eOzdXXuHDhAr1798bd3Z3SpUvj5+fHpk2b8imxiIiIiBRV+VpsL126RFxcHDdv3sTf3x8ASy5u4M7MzKRZs2Zs3ryZ6dOn89lnn1GpUiVat27Ntm3b8iu2SI4tXrzY7AhiJzTXpKBorklRlq/FtmrVqly+fJnNmzcTERGR6/fPnTuXw4cPs3TpUgIDA2nWrBlJSUlUr16d8PDwfEgskjv6BiAFRXNNCormmhRlBXaP7f1sl7tixQpq1KhBw4YN74w5OjoSFBTE7t27OXv2rDUjioiIiEgRVqg/PHbo0CFq1679h3EfHx8ADh8+XNCRRERERKSQKtTFNjU1lQq/P/f5v/w+lpKSUtCRRERERKSQyvEDGrZs2ULTpk1z9Nr9+/ff80prQTpy5Iipxxf7kJaWxr59+8yOIXZAc00KiuaaFIT86mk5LrY1atQgPj4+R6996KGH7jvQf3N1dSU1NfUP47+Pubq6/uHfeXp64uXlRVBQkFUyiPyVevXqmR1B7ITmmhQUzTUpCF5eXnh6elr1a+a42Hp4eBASEmLVg/8VHx8fDhw48IfxgwcPAuB9j0foeXp6smfPHn2wTERERKQQ8/T0NK/YmsHf35+BAweye/duGjRoAEBWVhaLFi3C19cXDw+Pe74vP/5DiYiIiEjhlu/Fdv369Vy7do2rV68Ct3cySEpKAqBt27Y4OzsD0LdvXxYuXMixY8fu3MoQEhLCzJkzCQgIIDIyEnd3d2JjYzl69CgbNmzI7+giIiIiUoRYjPvZYDYXqlWrxokTJ24fzGK5s5+txWLhl19+oUqVKgD06dOHhQsX3jUGtx+pGx4ezpo1a8jIyKBu3bq8//77Of4gm4iIiIjYh3wvtiIiIiIiBaFQ72P7V9LT0wkPD6dly5a4u7vj4ODA2LFjc/U1Lly4QO/evXF3d6d06dL4+fmxadOmfEosRV16ejpDhw6lcuXKODs7U7duXRITE3P03gULFuDg4HDPXxcuXMjn5FIY5WU+ae2S3LjfuaZ1S3Ijr73MGutaof7w2F+5dOkScXFx1KlTB39/f+Lj47FYLDl+f2ZmJs2aNePKlStMnz6dihUrEhMTQ+vWrdmwYQPPPfdcPqaXoqhTp07s2bOHCRMmUL16dT755BMCAwPJzs4mMDAwR19jwYIF1KhR466xez2IRGzf/c4nrV2SW3ldu7RuSU7kpZdZbV0zbMSlS5cMi8VijB07NsfvmTlzpmGxWIxdu3bdGcvKyjJq1aplNGzYMD9iShG2du1aw2KxGEuWLLlrvGXLlkblypWNW7du/en758+fb1gsFmPv3r35GVOKiLzMJ61dkht5mWtat+R+5baXWWtdK9K3Ivw34z5uFV6xYgU1atSgYcOGd8YcHR0JCgpi9+7d2gtX7rJixQpcXFwICAi4a7xPnz4kJyfz7bff5ujr3M9cFduTl/mktUtywxprl9Ytya3czhlrrWs2U2zvx6FDh+756F8fHx/g9tZkIr87dOgQNWvWxMHh7j82uZ0v7dq1w8nJCVdXVzp37qx5ZqfyMp+0dkluWGPt0rol+c1a61qRvsc2r1JTU+95j9DvYykpKQUdSQqxlJQUHnvssT+M53S+eHp6Mnr0aHx9fSlTpgwHDhwgMjISX19fdu7ceecPr9iHvMwnrV2SG3mZa1q3pKBYa10rNMV2y5YtOd6bdv/+/fds9SI5ZcZ8a9WqFa1atbrz+2eeeYa2bdvi4+PDmDFjWLFiRZ6PISJiTVq3pKgpNMW2Ro0axMfH5+i1vz+ZLK9cXV1JTU39w/jvY66urlY5jhQ+uZlvvz8wxNXV9Z5/Y8zLfHn44Ydp3Lgxu3btyvV7pWjLy3zS2iW5Ye21S+uW5AdrrWuFpth6eHgQEhJSoMf08fHhwIEDfxg/ePAgAN7e3gWaRwrO/cy32rVrs3jxYrKzs++6V80a8yU329SJbcjLfNLaJbmRX2uX1i2xJmuta3b94TF/f39+/PFHdu/efWcsKyuLRYsW4evri4eHh4nppLDx9/cnPT2dpKSku8YXLFhA5cqV7/okZ04dO3aM7du306hRI2vFlCIiL/NJa5fkhrXXLq1bkh+stq7lalOyQmjdunXGsmXLjHnz5hkWi8Xo2rWrsWzZMmPZsmVGRkbGndeFhIQYTk5OxsmTJ++MZWZmGt7e3kaVKlWMhIQE46uvvjL8/f2N4sWLG9u2bTPjdKSQa9mypVGhQgUjLi7O2LRpk9G/f3/DYrEYCQkJd73uXvOtefPmxvjx441Vq1YZGzduNKKjow0vLy+jbNmyxuHDhwv6VKQQyMl80tol1nC/c03rluRWTnpZfq5rRb7YVq1a1bBYLIbFYjEcHBzu+ucTJ07ceV3v3r3/MGYYhnH+/HkjODjYcHV1NZydnQ0/Pz9j48aNBX0aUkSkp6cbQ4YMMTw9PY0SJUoYderUMRITE//wunvNt2HDhhm1atUyypQpYxQrVsyoXLmy0atXL+Po0aMFeQpSiORkPmntEmu437mmdUtyKye9LD/XNYthaNdlERERESn67PoeWxERERGxHSq2IiIiImITVGxFRERExCao2IqIiIiITVCxFRERERGboGIrIiIiIjZBxVZEREREbIKKrYiIiIjYBBVbEREREbEJKrYiIiIiYhNUbEVERETEJqjYioiIiIhN+D/2sNuLzT1/mQAAAABJRU5ErkJggg==", "text": [ "Figure(PyObject )" ] }, { "metadata": {}, "output_type": "pyout", "prompt_number": 22, "text": [ "PyObject " ] } ], "prompt_number": 22 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Please submit pictures up to the quartic like the one above and answer the two questions." ] }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }