Skip to content

Instantly share code, notes, and snippets.

@cdeil
Created February 26, 2014 15:23
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save cdeil/9231500 to your computer and use it in GitHub Desktop.
Save cdeil/9231500 to your computer and use it in GitHub Desktop.
{
"metadata": {
"name": "",
"signature": "sha256:d8744dae97c9a36724599c528b9d822cfe5123464e3a0337e686961ddc67b733"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"%pylab inline\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from astropy.modeling.models import LogParabola1D\n",
"from gammapy.spectrum.crab import crab_flux"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"crab_ctools = LogParabola1D(amplitude=1e6 * 4.73232e-17, x_0=1, alpha=2.46315, beta=0.138834)\n",
"energy = np.logspace(-1, 1, 100)\n",
"plt.plot(energy, energy ** 2 * crab_ctools(energy), label='crab_ctools')\n",
"plt.plot(energy, energy ** 2 * crab_flux(energy, reference='meyer'), label='crab_meyer')\n",
"plt.plot(energy, energy ** 2 * crab_flux(energy, reference='hess_ecpl'), label='crab_hess_ecpl')\n",
"plt.loglog()\n",
"plt.ylim(5e-12, 1e-10)\n",
"plt.xlabel('Energy (TeV)')\n",
"plt.ylabel('E^2 dN/dE (TeV cm^-1 s^-1)')\n",
"plt.legend(loc='best');"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEWCAYAAACjYXoKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYjekbB/DvaUEqadO0KFTKUoNs2cq+bxkkpYQwJPzs\njMJgMmMZ+y6yjDUakiVCIjvZs0Qqa6G0n/P8/nhGlLZTnaW6P9d1Luq87/s8cc65e5/lvgWMMQZC\nCCGkBBRk3QFCCCFlHwUTQgghJUbBhBBCSIlRMCGEEFJiFEwIIYSUGAUTQgghJUbBhBBCSIlRMCGE\nEFJichtMnj9/jpEjR2LgwIEAgC9fvsDV1RUeHh7YvXu3jHtHCCHke3IbTGrXro3Nmzdnf33o0CEM\nGjQIGzduRGBgoAx7RgghJDeJBxN3d3fo6enBysoqx/eDg4NhaWkJc3Nz+Pr6Fnqd2NhY1KxZEwCg\nqKgokb4SQggpHokHk+HDhyM4ODjH94RCIcaPH4/g4GDcv38fe/bswYMHD+Dv749JkyYhLi7uh+sY\nGRkhJiYGACASiSTdbUIIIWKQeDBp27YtNDU1c3zvypUrMDMzQ61ataCsrAxHR0ccOXIELi4uWL58\nOQwMDJCQkIAxY8bg5s2b8PX1hYODAw4ePIhff/0Vffr0kXS3CSGEiEFJFo1+P2QF8LuOiIiIHMdo\naWlh/fr1Ob63devWQq9tZmaGp0+flk5HCSGkgjA1NcWTJ0+Kfb5MJuAFAoHErv306VMwxsr8w9vb\nu9y0W9JrFud8cc4p6rFFOa6gY2T1fyqJhyx+lvLy2hT3vJK+7op6TEl/CZdJMDE0NMye/wCAmJgY\nGBkZyaIrcsve3r7ctFvSaxbnfHHOKeqxRTlOVv9v0iaLn7O8vDbFPa+0XneS/j8TMMYkXhwrOjoa\nvXv3RmRkJAAgKysLFhYWCAkJgYGBAZo3b449e/agXr16JW5LIBBACj8SIWLz8fGBj4+PrLtBSJ5K\n+tkp8TuTIUOGoFWrVnj8+DFq1qyJbdu2QUlJCatXr0bXrl1Rv359DB48uFQCCSHyrKLctZCKSSp3\nJtJEdyaEECK+kn52ymQ1l8Q5OQFjxgBt2wISnOwnRN5paWkhMTFR1t0gckRTUxMJCQmlft3yGUxa\ntuTBRCDgf7q4ANWry7pXhEhdYmIi3amTHCS1mlZuc3OVhE9CAkLXrAHWrgXCw4HatYERI4Br12Td\nNUIIkSuhoaGlsjCkYsyZvH0LbN0KbNgAaGvzu5UhQwBVVdl0khApoTlEklt+r4mSvlYqRjD5SiQC\nTpwA1q0DLl78NrfSoIF0O0mIlFAwIblJKpiUy2GufCkoAN27A4GBwM2bgKYm0Lkz0K4dsHs3kJ4u\n6x4SQkiZVLGCyfeMjYH584EXLwAvLz4MZmwMzJgBPHsm694RQoqoVq1aCAkJkXU38uTj4wMXFxdZ\nd0MqKm4w+UpZGRgwADh9GrhwAcjMBFq0+HYHIxTKuoeEkAIIBIJSXaFkb2+PLVu2lMq1JJmHUN5Q\nMPle3brA0qXAy5d8gn7xYr4SbMECID5e1r0jpMIRyuCXudIMABVpvoqCSV5UVIBhw4BLl4AjR4DY\nWKB+fWDgQCAkBKhALxBCJCUmJgYODg6oUaMGdHR04Onpie3bt6N169aYPHkydHR0MG/ePDx79gwd\nOnSAjo4OdHV14ezsjE+fPuW41pUrV9CgQQNoaWnB3d0d6UWY/zxy5AgaNWoEDQ0NmJmZ4cSJE5g9\nezYuXLiA8ePHQ11dHRMmTAAAhIeHo1mzZqhevTqaN2+OS5cuZV8nLi4Offr0gba2NszNzXOUG/9e\nWloanJ2doaOjA01NTTRv3hxv374twb+gnGHlDADm7e3Nzp49W7oX/vSJsTVrGGvYkLG6dRlbtoyx\nhITSbYOQUlbYW5z/ZlTyh7iysrKYtbU1mzx5MktJSWFpaWksLCyM+fn5MSUlJbZ69WomFApZamoq\ne/LkCTt9+jTLyMhg7969Y+3atWMTJ07MvpaJiQmzsrJir169YgkJCax169Zszpw5BbYfERHBNDQ0\n2OnTpxljjMXGxrKHDx8yxhizt7dnW7ZsyT72w4cPrHr16mznzp1MKBSyPXv2ME1NTZbw3/u/bdu2\nbNy4cSw9PZ3dunWL6erqsjNnzjDGGPP29mbOzs6MMcbWr1/PevfuzVJTU5lIJGI3btxgnz9/Fv8f\nr4RyvybOnj3LvL29C32tFHrdEp0thyQeH0Uixi5cYMzJiTENDcbc3BiLiODfJ0TOyOvvi+Hh4UxX\nV5cJhcIc39+2bRszNjYu8NyAgADWuHHj7K9r1arFNmzYkP11UFAQMzU1LfAaHh4ebPLkyXk+Z29v\nzzZv3pz99Y4dO1iLFi1yHGNra8v8/PzYy5cvmaKiIktOTs5+bubMmczNzY0xljOYbN26lbVq1Yrd\nuXOnwL5JWn6viZK+VmiYS1wCAdCmDbBrFxAVBVhaAo6OQNOmwObNwJcvsu4hIXIvJiYGJiYmUFD4\n8SPo+yqsAPDmzRs4OjrCyMgIGhoacHFxwYcPH/I9x9jYGHFxcQW2/+rVK5iamub7/PfzJnFxcTA2\nNs7xvImJCeLi4hAfHw8tLS2ofrcB2tjYGLGxsT9c08XFBV27doWjoyMMDQ0xffp0ZGVlFdjPsoSC\nSUno6gLTpwNPngC//85Xfxkb86XGDx7IuneEyK2aNWvi5cuXeU6w554AnzVrFhQVFXH37l18+vQJ\n/v7+EIlEOY55+fJljr8bGBgU2n5+JWpzt29oaIgXL17k+N6LFy9gaGgIAwMDJCQkIDk5OUf7eRX7\nU1JSwty5c3Hv3j2Eh4fj6NGj2LFjR4H9LEsomJSG7zdD3rgBqKkB7dvzx/79fLkxISRbixYtoK+v\njxkzZiAlJQVpaWm4ePFinscmJydDVVUV1apVQ2xsLP78888czzPGsGbNGsTGxiIhIQELFy6Eo6Nj\nge2PGDEC27Ztw5kzZyASiRAbG4tHjx4BAPT09HKUsO3RowceP36MPXv2ICsrC3v37sXDhw/Rq1cv\nGBkZoVWrVpg5cybS09Nx584dbN26Fc7Ozj+0GRoaisjISAiFQqirq0NZWRmKiori/tPJrxINkskh\nufmR0tMZ++cfxtq1Y0xfn7HffmMsJkbWvSIVjNy8H/Lw8uVL1q9fP6atrc10dHSYl5cX8/PzY23b\nts1x3L1795iNjQ1TU1NjjRs3ZkuXLmU1a9bMfr5WrVrsjz/+YPXr12fVq1dnbm5uLDU1tdD2AwIC\nmLW1NVNXV2dmZmbs5MmTjDHGLl26xOrWrcs0NTWZl5cXY4yxsLAwZmNjwzQ0NFjTpk3ZxYsXs6/z\n6tUr1qtXL6alpcVMTU1zzN/4+PgwFxcXxhhje/bsYRYWFkxVVZXp6ekxLy+vH+aMpCG/10RJXysV\nKzeXrNy7B6xfz+dZ7OyAsWOBTp34HQ0hEiSX7wciU5TosYjk+s2TnMwDyrp1fKJ+zBhg+HBAS0vW\nPSPllFy/H4hMUKLH8kBNDRg9mieZ3LGD/1mnDg8oV6/KuneElCuLFi2Curr6D4+ePXvKumvlEt2Z\nyNq7d8C2bXwYTFsb+PVXYPBgoGpVWfeMlANl7v1AJI6GuYpIIBDA29sb9vb2sLe3l1g7KSnAhw9A\nQgKQmMgfnz9/e3z58u2Rmsqz26el8T+zsnj+yKwsXmKFMUDARLD9FIxf3q5Fg+TLOKk3DEEmY/Fe\n0xyVKwNVqvAsL1Wr8oe6+rdH9er8oanJHzo6/JgKlGOO5IOCCckt92siNDQUoaGhmDdvHgWT7wkE\nAugs0YFmFU1oqmhCW0UbOlV1oK2ijRqqNVBDtQb01PSgp6oHA3UD6KnpQUlBCQD/UH//Hnj16tsj\nLg54/ZrneXz9mt9IvHvHg4CODp/u0NTkH+YaGkC1avwDXk2NF3JUVeVB4GtAqFSJJypWUgIUFfmD\n95u3LxQCyq+eQ+/weugf34aPdRojqvOviK7fEykZSkhJ4YEsKenb49MnHsw+fuQB7v17fk0dHUBP\n79tDXx8wNAQMDAAjI74lRleXgk55RsGE5EZ3JkUkEAjwJvkNPqZ9REJqAj6kfMD7lPd4n/Ie71Le\nIfbTG7x49xavPsXjfXo8vojeQylTG4IkI2R+MIJySk1oKRnjpyomMNEwgZlubZj9pAt9fQH09IAa\nNfgHsKqqFD6E09KAAweANWt4ssnRo4GRI3lkKERKCg96b98Cb97wR1wcf8TG8kAZE8PXBNSsCdSq\nxadvatcGTE0BMzP+UFOT8M9IJIqCCcmNgkkRff0HefcOuHuXb0T/+nj8mH/A1qnDPyhNTQGT2lmo\nbvQGlXViIVR9hfcZMXj56SWiP0XjxccXeP7xOdKz0lFHsw7qaNaBmZYZzLXMYa5tDgttCxioG0in\nZsHNm8DatTy4dOsGjBsHtG5d4oj25QvPuB8dDTx/zh9PnvBMMU+f8rsuS0v+qFePJ09u0IDHM7qj\nkX8UTEhuFEyKSCAQoEYNhowM/qFXvz7/EKxXD7Cw4EM74m46/ZT2Cc8/PsfThKeISohC1IcoRCVE\n4fGHx0jJTIGFjgUsdSxRX6c+GtRogAa6DVBbszYUBBJYLPfxI+Dnx5cXV6nCJ+yHDpXILYRIxO9e\nHj0CHj7kj/v3eZBmDLCyAqytgZ9/Bho14v/eVaqUejdICVSEYFKrVi1s2bIFHTt2lHVXygQKJkUk\nEAgQE8NgaCid35wTUxPx6MMjPHj3APff3ce9d/dw7909fEj5gAY1GsC6hjWs9azR6KdG+Pmnn1Gt\ncrXSaVgkAs6c4UNg584Bzs58M2S9eqVz/QIwxofNIiOB27f549YtfidjYQE0acLzXjZvzgNOpUoS\n7xLJR0UIJrVr18aWLVvQoUMHWXelTKBgUkTy8ub5mPYRd9/eReSbSNx6fQu33tzC3bd3oa+mDxsD\nGzTVbwobAxvY6NtAo4pGyRqLiQE2bOBZixs04ENgffrwWX4pSk3lAeb6deDaNeDKFeDZM3730rIl\nYGvLH7mSwhIJkpf3Q3EJhcJC81eVp2AiEonyzKRcmiQVTOQ3cU8xyfOPlCnMZHff3GXbb21nE4Im\nsFZbWjHVharMcrUlGxYwjK25soZdj7vOMoWZxWsgPZ2x3bsZa92aMUNDxubPZyw+vnR/CDElJTEW\nGsrY4sWM9enDmI4OY8bGvBzM2rWM3b1LpWAkSZ7fDy9fvmT9+/dnurq6TFtbm40fP575+fmxVq1a\nsUmTJjFtbW3222+/sadPn7L27dtn5/AaOnQo+/jxY/Z1atWqxRYvXszq16/PNDU12fDhw1laWlqB\nbZ89e5YZGhqyJUuWMF1dXaavr88CAgLYsWPHmLm5OdPS0mKLFy/OPl4kErHFixczU1NTpq2tzQYN\nGpRdHKtHjx5s1apVOa5vZWXFDh8+zBhj7MGDB6xTp05MS0uLWVhYsH379mUf5+rqysaMGcO6d+/O\nVFVVWUhISIn/XQuT32uipK8V+X2lFZM8v3nykinMZDfjb7IN1zYw98PurN7qekxtkRqz22bHZofM\nZkGPg1hiaqL4F755k7FRoxirXp2xIUMYCwuTi09tkYixR48Y27yZ1xWrXZsxXV3GfvmFsdWrGXvw\nQC66WW7I6/tB1pUWz549y5SUlNiCBQtYVlYW27RpE9PW1mZOTk4sOTmZ3bt3j6moqLDo6GjGGGMr\nVqxgtra2LDY2lmVkZLDRo0ezIUOGMMYY27dvX47iWbdu3WLa2tosMzOTJScnMyMjI+bn58eEQiG7\nefMm09HRYffv32eM8WCioaHBwsPDGWOs0CBYGiQVTGiYSw4lpiYiIjYCF19eRFhMGK7FXYOppina\nmbSDnYkd2pq0RQ3VGkW72NcJ+7Vr+U7G8eMBJye52mH/8iWf9jl7Fjh9mu+16dgR6NwZ6NKlSCuh\nST4Kez8I5pXOxCLzFu89d+nSJfTt2xevX7/OMazj5+cHb2/vH+qHfO/w4cOYP38+bty4AYAPc82c\nORMeHh4AgOPHj8PT0zPfeiUA36jXo0cPfPnyBQKBAElJSdDQ0EBERASaNWsGAGjatCnmzp2LPn36\noF69elizZk32UFp8fDxMTEyQlpaGjIwMGBgY4OrVqzA1NcWUKVOQlpaG1atXY+/evVizZg3Onz+f\n3fbo0aNhaGiIuXPnws3NLfvnlhZJDXNJd1CdFImmiia6mXVDN7NuAIBMYSZuxN/AuRfnsO3WNowI\nHAGjakboULsD2tdqD/ta9tBU0cz7YtWrAxMnAhMmAKdO8Qn7GTMAV1c+YW9mJsWfLG/GxoCLC38w\nxpcmh4QAAQGApyff+9K1K9CjB59zUVaWdY/LD3GDQGkRt9Kil5cXwsLCkJSUBJFIBK1cyVHFrbQI\nANra2tnL+lVUVADwWiZfqaioZBe9evHiBfr375+jv0pKSnjz5g309fUxaNAg+Pv7w9vbG//88w8O\nHjyYfV5ERAQ0Nb+9P7OysjBs2DAA/AM8r0JaZREFkzJAWVEZLYxaoIVRC0xrPQ1ZoizcjL+JM8/P\nYMP1DXA97ApLHUt0qtMJnep0QuuarVFZqXLOiygo8E/krl35ZpL164FWrQAbG3630q2b+GumJUAg\nAMzN+WPMGF5X7MoV4PhxYNIkPqHfuTPQqxcPLjo6su4xKY7vKy3mnmAvqNJi9erVcfjwYXh6euY4\nRtxKi+IyNjbGtm3bYGtrm+fzrq6uGDZsGFq3bo2qVauiRYsW2efZ2dnh5MmTpdofeVQuswb7+Pgg\nNDRU1t2QGCUFJTQzbIbpbaYj2DkY76a+w5+d/4SiQBEzQ2ZC909d9NrdC6siVuFJQh63+rVrA76+\nwIsXvH69jw9Qty7w11882ZgcUVbmezN//52vErt/n8e9gAC+6bRdO97t7wrjkTJA1pUWxTVmzBjM\nmjUrO2i9e/cOgYGB2c/b2tpCIBBgypQp2XcdANCrVy88fvwYO3fuRGZmJjIzM3H16lU8fPgwu++y\nFhoaCh8fn5JfqEQzLnKoHP5IYvuQ8oHtvbuXuR92Z/p/6TOzlWbMM8iTBUcFs9TMfCrQRUQw5uLC\nJ+xHjOAT+HIuNZWxY8cY8/BgTE+PMSsrxubOZezOHZrE/0qe3w+yrLR49uzZHNfIzMxkCgoK7MWL\nF9nfa9OmDdu1axdjjK/mWrZsGbOwsGDq6urM1NSUzZ49O8c1FyxYwAQCAXv+/HmO7z969Ij17Nkz\ne9Vax44d2e3btxljjLm5ubHffvut6P9opSC/10RJXys0AV/OMcZw+81tBEUFISgqCJFvI9Ghdgf0\nMu+F3ha9f5zIf/uW71dZtw4wMeFDYA4Ocr/zUCgELl8GDh0CDh7kiTV/+QUYNIjvc6moqV/o/SA9\n/v7+2LRpU47Jdnkkk02LmZmZOHnyJM6fP4/o6GgIBAKYmJigXbt26Nq1K5SkvCmuKOjNU7D3Ke8R\n/CQYgY8CcfLpSTSo0QB9Lfqin2U/1NWu++3ArCwgMBBYtYrnUxk9GvDw4KmH5RxjfNPk/v3Avn08\na7OjI39YWMi6d9JF7wfpSElJQYcOHTB+/Hg4OzvLujsFknowWbBgAQ4ePAhbW1s0b94cBgYGEIlE\niI+Px5UrV3D58mX88ssvmDNnTrEblwR68xRdelY6QqNDcfjhYRx5dASaKprob9kfA+oNQKOfGn2b\nCL17l68C27uXT1iMH8+XVZWBX/cZAyIieNf37uXp952ceGAp5TlauVSR3w+LFi3C4sWLf/h+u3bt\ncOzYsVJr58SJExgwYAA6d+6MgwcPSnwHe0lJPZgEBgaid+/e+WbEFYlEOHr0KPr06VPsxiWhIr95\nSkLERIh4FYFDDw7h0MNDAIBf6v2CX+r/gqYGTfnr4OuelTVreOGW8eP5p/J/yyrlnVDI97Ls2gUc\nPsxzh7m6Av36ydW2m1JF7weSG+XmKiJ685QcYwy3Xt/CwQcHsf/+fmQKMzGowSAMbjCY37EwBpw4\nwYfArl0DRozge1aMjWXd9SJLSQGOHAF27OBzLQMGAO7uZeaGq8jo/UByk6tg4uHhgY0bNxa7UUmi\nN0/p+jqBv/fuXuy9txfKisoY0nAIhjQcAgsdC174ZM0awN8fsLfnuwzt7MrUJ3J8PO/+1q386+HD\nATe38rHznt4PJDepB5OEfPYbMMZgbW2N2NjYYjcqSfTmkRzGGK7GXcWeyD3Ye28vflL7Cc7WznBs\n6AgDqPNP5FWr+OYQT09eZ6UMjR8xBly6BGzZwleFdejA1xx07sz3fJZF9H4guUk9mCgoKMDExCTP\nk2JjY5GRkVHsRiWJ3jzSIRQJERodil2RuxDwMABNDZpimPUw9LfsB7Xzl3lQuXSJ/5r/66+8LnAZ\n8vkzsGcPsHEjkJjIF7MNH87LNpcl9H4guUk9mJiZmSEkJCTPgFKzZk3ExMQUu1FJojeP9KVmpiLw\nUSD87/gj7GUY+lr2hevPrrAXGUNh7To+aW9nx+9W7O3L1BAYAFy9yrPPHDrEU7iMH8/rs5SFH4Pe\nDyQ3SQWTfG/eJ06ciMTExDyfmzp1arEbJOWPirIKBjccjKNOR/Fo/CM00muESScmoc6/neDdSw3P\nb4XysaJx43iN302b+Ax4GdGsGR/6evaM/33YMF5Jcts2IC1N1r0jtWrVQkhIiFjnhIaG/pBQsjxR\nUFDAs2fPpNtmfk+MHz8ejRo1yvO5CRMmSKxDpGzTU9PDJNtJuD3mNgIGB+Bj2kc0/6cDOqoewO5/\nZiNjyWLg33/57voZM3j++TJCU5MnYH70CFi4kG+KNDEBfvsNKEKSWiIhAoEg3y0MRHrEmlb8Wi9A\n3pX3RI9lRWP9xvi7+994NekVRtuMxvY7O2BwxxUTxtbCo3/9+K/1jRsDAwcCFy7wGfAyQEGB790M\nCgLOn+dzKg0a8BT6N2/Kunfli1AolHUXyr3SSvQoVjC5evVqiRuUBh8fH9jb28u6G+Q/lZUqY1CD\nQTjhfALXPa5Ds4omOoaPhq1VBPyPLEBGG1u+V8XGBti+HUhPl3WXi8zCAli9mmf1t7ICevfmq8CO\nHQNEIln3Tr7FxMTAwcEBNWrUgI6ODjw9PbF9+3a0bt0akydPho6ODubNm4dnz56hQ4cO0NHRga6u\nLpydnfHp06cc17py5QoaNGgALS0tuLu7I72Ir6Fly5ZBT08PBgYGOQpUpaenY8qUKTAxMcFPP/2E\nsWPHIu2/Mc3379+jV69e0NTUhLa2Ntq1a5d9nq+vL4yMjFCtWjVYWlrizJkzBbbPGMMff/wBMzMz\n6OjoYPDgwTmmF8LCwtCqVStoamrC2NgYO3bsAAC4ublhzJgx6NKlC6pVqwZ7e/scafjFYW9vL/2s\nwV27dhXncJkQ80ciMpIpzGRHHh5hPXf1ZNq+2szrqCd7sXsdY1268BTAc+fKvH59cWRkMLZzJ2ON\nGjHWoAFjfn6MpafLrj/y+n6Ql7K93t7eLCsriwUFBbGqVatm15afOHEi69u3L0tMTGRJSUmsd+/e\nbObMmYwxxmbMmMHGjBnDsrKyWFZWFgsLC2OMMfbw4UNWs2ZNFv/f6/bFixfs6dOnBfajoHLA0dHR\nTF1dnf3zzz8sKyuLffjwgd26dYsxxsv9qqurswsXLrD09HTm5eXF2rRpk31dgUCQb9v5vSZK+lqR\nz1daCcjrm4fk73niczY7ZDb76a+fmL2fPQs+soxljfbg6fBdXBi7dk3WXRSbSMTYyZOMderEmJER\nY8uXM5acLP1+FPp+4IOLJX+IKTw8nOnq6jKhUJjj+9u2bWPGxsYFnhsQEMAaN26c/XWtWrXYhg0b\nsr8OCgpipqamBV7j7NmzTEVFJUf7NWrUYBEREUwkEjFVVdUcH8bh4eGsdu3ajDHG5s6dy/r27cue\nPHmS45pRUVGsRo0a2YGvKOrVq8dCQkKyv46Li2PKysosKyuLLVq0iDk4OOR5nqura3bQYYyx5ORk\npqioyF69esUYk00wKXSYKzIysuS3P4QUoFb1Wvi9w+94MfEFxtiMwaKEwzCuexS+2z2QZG4C9O8P\ntG3Lc8tnZcm6u0UiEPAFbKdO8TxgFy/ymmQLFvA5FrlRWuFETOKW7XV0dISRkRE0NDTg4uKCDx8+\n5HuOOGV7v2+/atWqSE5Oxrt375CSkgIbGxtoampCU1MT3bt3x/v37wHw1axmZmbo0qULTE1N4evr\nC4Bvp1ixYgV8fHygp6eHIUOGID4+vsA+REdHo3///tnt1K9fP7sc8KtXr1CnTp08z8td7ldVVRVa\nWlpF+rklpcBgEhISgrFjx0qrL6SCq6RYCYMbDsY5t3MIHhqM5wqfYFxpNVyWtkaUUzdg2TJes/6v\nv3jSyTLCxoav/Dp/ni8vNjMDZs0C/vtsqpC+L9ubW0Flez99+gR/f3+Ick1IlWbZXh0dHaioqOD+\n/ftITExEYmIiPn78iM+fPwMA1NTU8Ndff+Hp06cIDAzEsmXLsudGhgwZggsXLuDFixcQCASYPn16\ngW0ZGxsjODg4u53ExESkpKTAwMAANWvWxNN8SogyxnLs9UtOTkZCQkKplysWR77BZOfOnZgyZQoC\nAgKk2R9CAABWelZY32s9nk14hsZGzdA1YwtaDM9C8OIREN64DtSpwzdBRkXJuqtFZmnJ96Zcv87v\nTiwsgClTgDdvZN0z6ZPnsr0KCgoYNWoUJk6ciHfv3gHgWT++1nE/duwYnjx5AsYYqlWrBkVFRSgq\nKuLx48c4c+YM0tPTUblyZVSpUuWH+va5FVQOeOjQoTh9+jT279+PrKwsfPjwAbdv384+NygoCBcv\nXkRGRgZ+++032NrawtDQsNg/d0nlG0zc3d1x8OBB6OrqSrM/hOSgqaKJybaTEeUZhVltZmFJxlkY\nNz2PlZtHIUVFGWjVii+hOnu2zCwtrlWLF7K8fRvIyADq1wemTuVFLisKBQUF/Pvvv3jy5AmMjY1R\ns2ZN7N+/P889I97e3rhx4wY0NDTQu3dvDBgwIMcxAoEAQ4cOzR52Mjc3L1KdpYL2pvj6+sLMzAwt\nW7aEhoaW3YG3AAAgAElEQVQGOnfujMePHwMAoqKi0LlzZ6irq6NVq1YYN24c7OzskJ6ejpkzZ0JX\nVxf6+vp4//59nvVUvufl5YU+ffpkr8qytbXFlStXAPC7t6CgICxduhTa2tpo3Lgx7ty5k913Jycn\nzJs3D9ra2rh58yZ27txZpJ9NUvJNp7J48WKcPHkSQUFBUCkj9SoASh9REdx5cwcrLq9AwMMAOJs6\nYNZLE+hv/oeXFp44ERgyhNftLSNiYwFfX15nZcQIYNo0QEendK5N74fyafjw4TAyMsKCBQvEPlfq\n6VRmzpwJd3d39O3bt9gXJ0QSrPWssbXvVjwc9xBa2kZoJFqDvrNNETlpKNiePfxX/wULgP+GKOSd\noSGwciW/U0lO5sNfc+eWqWkhImXy+AtCgRPwLi4umDZtmrT6QohY9NT0MK/9PER7RaOHRS8M+LwJ\nLX5JxIn1UyF68QKoW5fnkL9/X9ZdLRIjI2DtWl5v7NUrwNyc37GUoTRmcmXRokVQV1f/4dGzZ0+p\n9qN79+559uOPP/4o9jXlMYUMVVok5YaIiRD4KBBLLi7B2y9vMcfCA0PDk6C8YRNP2zJ5MtCpU9lI\n9wvg4UOe9ys8HJgzBxg5kpeKEQe9H0huUh/m+mrfvn3ZS+IWLFiA/v3748aNG8VukBBJURAooJ9l\nP4SPCMf2fttxKCEMNdU3YcnOX5HStwcwaRLPWuznVyZStlha8iXFgYFAQADP/7V/f5lZZ0AqmELv\nTKysrBAZGYmwsDDMmTMHU6ZMwYIFCxARESGtPoqFfhMj37v39h6WhC/B0cdHMbLRCExPaQytdX5A\nZCRPiT9mDKCtLetuFsmpU8D06YCSEvDnn7xETGHo/UByk9mdydd10kePHsWoUaPQq1cvua2ySEhu\nDWo0wPZ+23HD4wZSslJh9mgcxk0wQ9z+bcCTJ3wH4bhx/O9yrnNnPp8ycSKvUd+3Lx8KI0QeFBpM\nDA0N4eHhgb1796Jnz55IS0v7YfcpIfLOpLoJVvVYhQfjHqBa5WqwuuiE4f2AZ2FHgerVAVtbwMGB\nT1DIMQUFwMkJePCAZ5hp25bHwvx202tqamZP1tKDHgKBAJqamhJ5bRY6zPXlyxcEBwfD2toa5ubm\niI+PR2RkJLp06SKRDpWUQEC39aRwiamJWBmxEquvrkbnOp3xm81k1Dt6GVi+nBd6nzqV/+pfyA5m\nWfvwAZg3j9ernzGDlxQuQ1tsiBwp6WdnuVzN5e3tDXt7e6ppQgr1Of0z1l5di+WXl8POxA6/tZ4F\nq0tPef6vd+/4CjA3N6BqVVl3tUAPHvDULI8eAUuXAn36lJlFa0TGQkNDERoainnz5lEw+R7dmZDi\nSM5Ixrqr67D00lK0MW6Due1+g/XTZD7THR4OjB3Lf+2X8/RCJ07wRWtGRsCKFTxVCyFFUdLPTrEq\nLRJSXqlVUsPU1lPxdMJT2BrZosvOrhgYtwKRGxbwksKvX/NNkGPHyvVkfdeufCd9z558tdekSUCu\nooSESAQFE0K+o1pJFf9r9T88nfAULQxboLN/ZwyOnIsHv0/kS6d0dPhk/S+/AP8l5JM3ysqAlxff\n+P/lC1CvHq+GTOtmiCQVa5jr694TeUTDXKQ0JWckY/WV1Vh2aRm6mnXF3HZzYV5ZH9i6lU9O1K7N\nMzN27y63kxRXr/IROgUFnq6lcWNZ94jII4lNwB88eDDfxkaPHp1ddUzeUDAhkvA5/TNWXF6BlREr\n0c+yH+bazYVxVX1g3z5gyRL+a/+0aYCjo/g5T6RAJOLxb/Zs3sX58wENDVn3isgTiQUTZWVlODk5\n/VBWkzGGAwcOIDk5udiNShIFEyJJCakJ+Cv8L2y4vgFDrYZidtvZ0FOtwWe+lywBnj7lK8BGjgRU\nVWXd3R+8fw/MnAkcO8ZvrBwd5faGikiZxIJJkyZNsH37dlhZWf3wXM2aNXOUjJQnFEyINLxJfoPF\nYYvhf8cfY5uOxZRWU1C9SnU+j7JkCa/R++uvfHyptIqTlKJLl3gmGT09PvRlZibrHhFZk9hqrhUr\nVqBatWp5Pnfo0KFiN0hIeaCnpocV3VbghscNxCfFo+6qulhycQlSG1sBBw4AYWG86lXdujz/yXc1\nyuWBrS0vH9y1K9CyJS//QlmSSEnQPhNCSsHD9w8x+8xsXIm9Am87b7g1coOSghIQF8d31W/dyncS\nTp/O0wHLkZcveUqWZ8+AjRuB1q1l3SMiC1LdZ9KkSZNiN0RIeWapY4mDgw7iwMAD2BW5Cw3XNkTA\ngwAwfX2+8fHJE8DUlG/+GDCAZ2yUE8bGPM29jw8wcCDfSkN7U4i4xAom9Bs/IQVrYdQCZ4adwYpu\nK+BzzgdttrVB2MswQFOTV7h69owHlP79gS5dgLNn5aJAiUDAA8m9e4BQCDRsCPz7r6x7RcoSsYKJ\ntMtdElIWCQQCdDPrhhseNzDaZjSGHhqKfv/0w6P3j/gKrwkT+KovR0dg9Gg+rnT0qFwEFU1NPtS1\nYwffPT9kCE9RRkhhaM6EEAlLy0rDqohVWBK+BAPrD4S3nTf01PT4k0IhcPAgsGgRDyazZvHd9XKQ\nrTglBfD2Bvz9gb//BgYNomXE5ZnE50y2bdtW7IsTQoAqSlUwtfVUPBz3EJUVK6PB2gZYdGERUjNT\nedAYNAi4eZMHlL//5tkZ/fyAzEyZ9rtqVT7dc+QI3+Q4YABPUUZIXgoMJvPnz8fp06el1RdCyjXt\nqtpY3m05Lo+8jJuvb8JitQX8b/tDxET8V/6ePYGLF4H16/ntgLk5/7uM69W3aAHcuMFzfP38M7B7\nt1yMyBE5k+8w16hRo5CUlITdu3f/sAtentEwFykrLr68iMknJ0PERFjWZRnamrTNeUB4OLBwIU8D\nPHUqMGqUzOuqXLsGuLry1c3r1vE6YqR8kNgOeHV1dURERKB+GSuIQMGElCUiJsLeu3sxI2QGmhk0\nw5LOS1BHs07Og65fB37/Hbh8Gfjf//jWdTU12XQYQFoaX0bs5wesWcOHv0jZJ7E5k3///RcDBw7E\nEzmu3UBIWacgUMAQqyF4OO4hGv/UGM02NcP0U9PxOf3zt4NsbICAAJ7/6+pVvl9l8WIgKUkmfa5S\nBfjjD+DwYZ7ny8UF+PhRJl0hciTfYGJvb4+9e/fC2dlZmv0hpEJSUVbB7HazcXfsXbxNeQuL1RbY\nfGMzhCLht4OsrYG9e/nelLt3eVD5/XeZ7TBs2ZKvG9DQAKysgFOnZNINIicKXRocGxsLQ0NDafWn\nxGiYi5QH1+OuwyvYCymZKfi7298/zqcAvOD7778Dx4/zvSteXjLLK3/qFODuDjg48LsWFRWZdIOU\ngMSXBpelQPKVj48PQkNDZd0NQorNxsAGF4ZfwNRWU+F0yAmOBxwR8ylXpm4LC77qKzycb4I0NeVr\neGVwp9K5M18n8Po10LQpcOuW1LtAiik0NBQ+Pj4lvk6hdybPnj3DqlWrEB0djaysLH6SQIDAwMAS\nNy4JdGdCypsvGV/ge9EXa66uwcQWEzGl1RSoKOfxq39UFF/9dezYtzuVfDJ/SwpjwK5dvKTLtGn8\nzzK0GLRCk9hqrq+sra0xcuRINGzYMHuJsEAggJ2dXbEblSQKJqS8ep74HFNOTcGN+BtY1mUZ+ln2\ngyCvLelRUTynfHAwT3/v6Qmoq0u1r9HRgLMzn6zfvh0ogwMcFY7Eg0nz5s1x5cqVYjcgbRRMSHkX\n8iwEE4InwKiaEVZ2WwkLHYu8D3z0iA97nT7NlxSPGyfV6o9ZWXzR2erVwIYNQL9+UmuaFIPEg4m/\nvz+ePn2Krl27onLlytnfl9d09BRMSEWQKczE6iursfDCQrg3dsdv7X6DeuV87j7u3wfmzQPOneP1\nVMaMkeoM+aVLwNChQPfuwF9/0eS8vJJ4MJkxYwb8/f1hZmaWYyf82bNni92oJFEwIRXJ6+TXmH56\nOkKehWBpl6UY1GBQ3kNfAHDnDjB3Lt/GPmsWr1NfqZJU+vnpE0+QfO8esGcPT3FP5IvEg4mpqSke\nPHiASlJ60ZUUBRNSEYW9DMO4oHHQqaqD1d1Xo55uvfwPvnoV+O034OFDnhbYxQVQUpJ4Hxnju+an\nTeM5LUeOpCzE8kTiS4OtrKyQmJhY7AYIIZLXxrgNrntcRz+Lfmjn1w7TT01HckZy3gc3a8Yn5/39\n+ad7w4Z8M6RIJNE+CgTA8OHA+fPAqlWAkxPw+XPh55GyodA7Ezs7O9y5cwfNmjXLnjOhpcGEyK/X\nya8x9dRUnIs+h+Vdl8OhnkP+Q1+M8R2Hc+bwGfOFC4Fu3SR+y5CaypcNnzoF7NsHyOkUbIUi8WGu\nvDb/0dJgQuTfuehz+DXoV5homGBV91Uw1TLN/2DGeLKtOXMALS2+DKtNG4n3cd8+vshs3jxee56G\nvWRH4sHk2bNn0NfXh8p/SzBSU1Px+vVr1K5du9iNShIFE0K+yRBmYPml5fgz/E94tfDCtNbTUFmp\ncv4nCIXAzp18LqVhQz65YW0t0T5GRfH6YHXrAps2SX2fJfmPxOdMBg4cCMXvSogqKChg0KBBxW6Q\nECI9lRQrYXqb6bjucR3X4q/Ber01zjw/k/8Jioq8YMmjR0CXLvzh7Aw8fy6xPpqb8+XDmpo8Fcud\nOxJrikhQocFEKBTmWMlVuXJlZGRkSLRThJDSZVLdBEccj2BJpyUYfmQ4XAJc8PbL2/xPqFyZp2SJ\nigLMzPinvJcX8O6dRPpXpQovKuntDXTsyNcFkLKl0GCio6ODI0eOZH995MgR6OjoSLRThBDJ6GvZ\nF/d+vQd9NX00XNsQm29s5mWD86OuzithPXjA51Xq1eOZir98kUj/hg7leyt9fXlhybQ0iTRDJKDQ\nOZMnT55g6NChiIuLAwAYGRllb2KURzRnQkjR3H59G6OPjoaSghI29NqABjUaFH7S06fA7NnAhQv8\nNsLdXSJ7VJKT+T6UqCjg4EGgVq1Sb4LkIvEJ+K+S/qvqpi7lhHHiomBCSNEJRUJsuL4B3qHeGG0z\nGnPazUEVpSqFn3j1Kt99+OYNv43o1avUl2IxBqxcyReW+fnxFctEcqQWTMoKCiaEiC8uKQ5ewV64\n9foWNvTagA61OxR+EmM83f20aUCNGjzxVtOmpd63sDDA0RHw8OArlymlvWRQMMmFggkhxffvo38x\nLmgcOtbpiL86/wXtqtqFn5SVBWzbxoe9OnTgy4mNjUu1X/HxwMCBgLY2sGOHzApKlmsSWxqcmZlZ\n7IsSQsqm3ha9ce/Xe1CvpI6G6xpiT+Sewj9glJT4bPnjx0CdOkDjxsDMmaWaK0VfHzhzBqhZk2eD\nuXev1C5NSkm+wcTQ0BAjR45ESEgI/aZPSAWiXlkdK7uvRMDgACwKW4Ree3rh5aeXhZ+opsbrp9y+\nDcTF8bLCmzbxjZCloFIlXhtl9mzA3h4ICCiVy5JSkm8wuX//Ppo2bYoFCxbAyMgIXl5euHz5sjT7\nRgiRoZZGLXHd4zpaGbVCkw1NsCpiFYSiIgQGIyNeXvHoUZ5MsnFjICSk1Prl6goEBfFtL97eEs9P\nSYqoSHMmcXFx2LdvH/bu3Yu3b99i8ODBWLRokTT6JzaaMyGk9D18/xAe/3ogU5SJLX22oL5u/aKd\nyBhw6BAwdSpPy/LXX3wTZCl48wb45ReeSszfn9KwlJTE06kAgIGBAUaMGIExY8ZATU0NmzdvLnaD\nhJCyx1LHEqFuoRhmPQx2fnaYf24+MoRFyIQhEAADBvBqjy1b8se0aaUyn6Knx2949PUBW1u+BYbI\nToHBJDU1Ffv27YODgwPMzMxw5swZ+Pr6Zm9gJIRUHAoCBYxtNhY3PG4gIjYCTTc2xdXYq0U7uUoV\nYMYM4O5d4P17wNKSrwAr4RhVpUo8Dcv48UCrVqU6mkbElO8wl5OTE06dOgU7OzsMGTIEPXr0yM4c\nLM9omIsQyWOMYXfkbkw+ORmuP7tinv08qCiL8flw9SrP/ZWVxStltWxZ4j6FhvL9KHPm8LT2lM5e\nPBLbZ7J9+3Y4ODjI/Y733CiYECI9b7+8hedxT9yMv4ktfbagrUnbop8sEgG7dvE7li5dgD/+4GNX\nJfD8OdC7N9C2Ld89r6xcostVKBKbM3F1dYW6ujpev36NESNGoNt/uQzu37+PLVu2FLtBQkj5UUO1\nBvb+she+nXzheNARE45PyL9ccG4KCrz+/IMHgI4Or5+yYgW/Wymm2rWB8HAgJoanX0lIKPaliJgK\nnYB3c3NDly5dsudJzM3NsXz5col3jBBSdvSv1x+RYyPxMe0jrNcVUjMlt2rVgD//5Mkjjx3jS4nP\nny92X6pVA44c4aWAW7TgpVmI5BUaTN6/f4/BgwdnF8hSVlaGkgSyhBJCyjYtFS3s6L8Dq7qvguth\nV4w9OhZJ6UlFv4ClJXDyJN884uzMH69fF6svioo8Ps2YAbRrx3fPE8nKN5hk/Xerqaamhg8fPmR/\n//Lly9CgxDiEkHz0rNsTkWMjkSHMgPV6a4Q8E2OJlUDAN488eMA3P1pZ8Qn6Yg59jRgB7N0LDBnC\nN+MTycl3Ar5Jkya4ceMGrl+/Dk9PT9y7dw8NGjTAu3fvcODAAfz888/S7muR0AQ8IfIj+EkwPP71\nQA/zHviz859Qryzmgp779/nSrI8f+RrgFi2K1Y+oKKBnT6BvX54xnzIP/0hiE/BfL2pjY4Pz588j\nPDwcGzduxP379+U2kHzl4+OD0NBQWXeDkAqvm1k3RI6NRKYws/D683mpX5+PUU2ZAvTrB/z6Kw8s\nYvpaZ/7KFZ59OCVF7EuUW6GhofDx8SnxdfK9MzEyMsLkyZPBGPshYgkEAkyePLnEjUsC3ZkQIp+C\nooIw+uho9KnbB76dfaFWSU28CyQmArNm8dn1pUv5phIxN5Okp/MKjo8eAYGBwE8/ideF8kxidyZC\noRBJSUlITk7O/vPr42vVRUIIKaoe5j0QOTYSXzK/4Of1P+P8CzFXbGlqAuvW8VxfixcD3bvzjSVi\nqFyZ10Pp2ZOnYHnwQLwukPzle2fSuHFj3Lx5U9r9KTG6MyFE/gU+CsTYY2MxsP5ALOq4CFWVq4p3\ngcxMYNkyvmRr2jRg0iSxdyhu385P3bcPsLMTr/nySCqJHgkhpDT1seiDO2Pu4O2Xt2iyoQkiXkWI\ndwFlZWD6dD4Jcvo00Lw5cP26WJdwdQV27+ZzKLt2idc8+VG+dyYvX76EcSGlN5OSkuQu3QrdmRBS\ntuy/tx+exz3h3tgd3nbeqKxUWbwLMAbs3Mkn6YcNA+bNA6oW/U7n7l0+7DV2LI9PFTWnl8Ryc3Xq\n1AkWFhbo27cvmjZtCi0tLQDAhw8fcO3aNRw+fBhRUVE4ffp0sRuXBAomhJQ9b5LfwOOoB54nPod/\nf3/8/FMxVoy+fQtMnAhERACbNwPt2xf51NhYoEcPoE0bntPrvz3aFYrEggkAnDlzBrt378bFixez\n06kYGBigTZs2GDp0KOzt7YvdsKRQMCGkbGKMYfvt7Zh6aiomt5yMqa2nQkmhGNk2jh7ltxk9egBL\nlgBF3GT9+TPg4MCrD+/eLdbNTbkg0WBSFlEwIaRse/HxBYYfGY60rDTs6L8DZlrFqMz46RMfszp2\nDNiwgQeWIsjI4Lvmnz4F/v0X0NYWv+myiibgCSHlikl1E5wedhqDGwyG7RZbrL+2XvwPOQ0NvmPe\nz4/voB8+nO9TKUSlSnyVV9u2/PHyZfF+hoqIggkhRO4oCBTg1dIL593OY/ONzeixuwfikopR4bVj\nR+DOHT5mZWUFBAUV3rYCT7ni4QG0bg1ERhbjB6iAKJgQQuRWPd16uDTiEpobNEfjDY2x/95+8S+i\nrg6sWQP4+/O7lBEjilSDfuJEvo2lUycgLKwYna9g8g0mZ77L2fw81y7TQ4cOSa5HhBDyHWVFZcxr\nPw+BjoGYfWY2XAJc8DFN/PxcaN+e36UoKfG7lCIUjHd05HtQHBx4+hWSv3yDyf/+97/svzs4OOR4\nbsGCBZLrESGE5KGFUQvcHH0T6pXU8fP6nxEaHSr+RdTV+YT8xo2Amxvg5VVo1sdOnfjo2OjRwNat\nxep6hUDDXISQMkO1kirW9lyLdT3XwemgE6aenIr0rHTxL9S1K79Lef+el2S8erXAw5s2Bc6dAxYs\n4PMp5EcUTAghZU4P8x64PeY2ohKi0Hxzc9x9e1f8i2hq8jGs+fOBXr34nwUU4apbl8+d7NjBVx3T\nDoSc8t1noqGhATs7OzDGcOHCBbRt2zb7uQsXLuBjMWoKSAPtMyGk4mCMYevNrZgRMgNz2s6BZwtP\nKAiK8TtybCwf9kpO5qlZTE3zPfTDB75txdqarz4uL7vlJbZpsbDiUvK4+x2gYEJIRfQk4QmcDzlD\no4oGtvXdBgN1A/EvIhLxEsG//853zru55ZuoKzkZ6N8fqF6dx57KYqYTk0e0Az4XCiaEVExZoiz8\nfv53rL+2Hut6rkP/ev2Ld6HISMDJiVd5XL+eD4flIT2dH5aczEusqKqWoPNyQGLBxMrKqsBG79y5\nU+xGJYmCCSEV26WYS3AOcEaHWh2wvNty8Ss6AkBqKi92EhjIbz2+G+b/XlYWMGoU8Pgxz9xSvXoJ\nOy9DEgsm0dHRAIC1a9cCAFxcXMAYw67/Ev/7yumSBgomhJCk9CRMCJ6Aiy8vYpfDLjQzbFa8Cx09\nyuv8/vorMHt2nhMkIhGvzXX+PHDiBFCjRgk7LyMSH+Zq1KgRbt26leN78lyFkYIJIeSrfff2YXzQ\neEy2nYypraZCUaEYs+VxcYCzM48au3YBhoY/HMIYMHcucOAAr9WVxyFyT+KJHhljCPsul8DFixfp\nw5oQUiYMajAI1zyu4fiT4+jk3wkxn2LEv4iBAXDqFNC5M2Bjk2d+L4GA70EZPhxo107s0vTlQqF3\nJtevX8fw4cPx6dMnAED16tWxbds2NGnSRCodFBfdmRBCchOKhPC96Iu/I/7G2h5rMaD+gOJdKCyM\nz7o7OgILF+ZZd37tWuCPP4CTJwFLyxJ2XIokNswVHh4OW1tbCP5bGvd1X0l1OZ9homBCCMlPxKsI\nDD00FB1qd8DyrsuhWqkYS7Dev+cF5BMTgX/+AfIob759OzBzJp9DKWAtk1yR2DDXjh070KRJEwwe\nPBh+fn5IS0uT+0BCCCEF+ZrfK0OYAZuNNrgZX4y5Xx0dXjmrXz+geXMgOPiHQ1xdgeXL+cjYjRul\n0PEyoNBhrgcPHuD48eM4efIkPn78iA4dOqBbt25o3bo1FOVw6yfdmRBCimJP5B54BXthZpuZ8Grp\nVbyd8xcuAEOGAO7ugLf3D6u9Dh/mCSKPHAFatiyljkuIVDctpqSk4OzZszh+/DguXbqE69evF7th\nSaFgQggpqueJz+F0yAnVq1SHX18/6KnpiX+RN294QFFQAPbsAXR1czx9/Di/Uzl0CGjTppQ6LgG0\nAz4XCiaEEHFkCjMx79w8bL25FX79/NDFtIv4FxEKgTlz+NLhAwf48Nd3Tp/m8/b79gFymolKcsFE\nTU0te/L9e0KhEOnp6RAKhcVuVJIomBBCiuPs87NwCXDBkIZDsLDjQlRSrCT+RQIC+LjW77/zrfHf\nfYaGhgKDBvF407lz6fW7tEhsAj45ORlJSUnZj/j4eMyePRt6enqYOHFisRskhBB51L52e9wacwuP\nPjxC662t8SThifgX6d+fLx/++28eTNK/1Vqxt+dDXUOH5jlnX+YVOuP08eNH+Pj4wMrKCklJSbh2\n7RqWLl0qjb4RQohU6VTVwRHHIxhmPQy2W2yxO3K3+BepWxe4fBn4+BGws+Pp7f/Tpg2flB82LM+9\nj2VavsHk3bt3mDFjBho3bgxFRUXcunULCxcuhLa2tjT7RwghUiUQCODZwhOnXE5h/rn5GH5kOJIz\nksW7iLo6sH8/0KcPnz+5eDH7qVateP5INzeeHLK8yHfORFVVFTo6OnB3d8+eP/l6qEAgwOTJk6Xa\n0aKiORNCSGlJzkiG53FPXIq5hL2/7MXPP/0s/kW+LudatIgnjfxPRATQuzewZQv/U9YkNgHv4+OT\n3cD3GGMQCATw9vYudqOSRMGEEFLadt7ZiUknJsHHzge/Nvs1z8VJBXr0iN+ldO0KLFsGKCkBAK5c\n4RWDt27lf8oSLQ3OhYIJIUQSoj5EwfGgI0w0TLClzxZoquRdNCtfiYl8P0pWFl8jrKUF4FtA2bYN\n6NlTAh0vIokFE09Pzzwb+RqRV65cWexGJYmCCSFEUtKz0jHt1DQceXQEewbsgW1NW/EukJXFi24d\nPcofdesC+DbkJcuAIrGlwTY2NrCxsUF6ejpu3LiBunXrwtzcHDdv3kRGRkaxGySEkLKqslJl/N39\nb6zsvhL99vaDb5gvRExU9AsoKfFhrqlTefXGM2cAAC1a8HRfw4eX3WXDhQ5ztWjRAmFhYVD+L9Vy\nZmYm2rRpg4iICKl0UFx0Z0IIkYaYTzEYcnAI1CqpYUf/HaihKmaJxbNneSr7rxscAVy6BPTtC+ze\nDXTqJIFOF0DixbE+fvyIz58/Z3+dlJSUnY6eEEIqqpoaNXHW9Sya6DdBkw1NEBodKt4F2rfnGxyX\nLAFmzABEItja8o2NTk481pQlhd6ZbNu2DT4+Pmjfvj0YYzh37hx8fHzg5uYmpS6Kh+5MCCHSduLJ\nCbgdccPYpmMxu+1s8coDv3/P09nr6wM7dgAqKjh3Dhg4EDh4kI+GSYNUVnPFx8cjIiICAoEAzZs3\nh76+frEblDQKJoQQWYhLioPTQScoKihil8Mu/KT2U9FPTkvjaeyfP+eTJzo6CAnhi78CA6WTvp6W\nBudCwYQQIitZoizMPzcfm29shn9/f3Ss07HoJzMGzJ7Nsw4fPw6Ymmbvdzx+nJeflyQKJrlQMCGE\nyFrIsxC4BLhgtM1ozGk3R7xhr3XrgAULeBKv5s1x5AhPRHzyJGBtLbk+UzDJhYIJIUQexCfFY+ih\noYSRsYsAABAQSURBVBAIBOIPewUGAiNGAH5+QM+e2LcPmDiRryS2tJRMfyW+mosQQoj49NX1ccrl\nFNrUbAObjTY48/xM0U/u04fPnYwYAWzfjkGDgMWLeR2UZ88k1+eSyDeY3LlzBy1btoSRkRE8PDyQ\nmJiY/VzzXFXECCGE/EhRQRHz2s/D9n7b4XzIGQvOLYBQVMTCgi1b8vXB3t7AkiVwHcYwaxbQsSMQ\nEyPZfhdHvsFk7Nix8PHxQWRkJOrWrYvWrVvjyRNeLCYzM1PiHXv+/DlGjhyJgQMH5vk1IYSUFZ3q\ndMI1j2sIeR6C7ru64+2Xt0U7sV49vhdlxw7gf//D2NEieHryDY1v3ki2z+LKN5gkJSWhW7du0NTU\nxJQpU7BmzRp069YNly9flkrHateujc2bN+f7NSGElCUG6gY4Pew0mhk0g81GG1x4caFoJxoZARcu\n8IJbI0Zg8oQsDB0KdOkCJCRIts/iyDeYCAQCfPr0Kfvr9u3b49ChQ3B2dsbLly+L3IC7uzv09PRg\nZWWV4/vBwcGwtLSEubk5fH19i9F1QggpW5QUlLCw40Js7LURA/cPLHpuL01N4NQpIC4OGDQIv01N\nQ+fOQPfuQFKS5PtdFPkGk2nTpuH+/fs5vmdtbY0zZ87AwcGhyA0MHz4cwbkylwmFQowfPx7BwcG4\nf/8+9uzZgwcPHsDf3x+TJk1CXFycmD8GIYSUHd3Nu+PqqKs4/Ogw+uzpg4TUItxiqKryVV5KShD0\n7oU/vZPx8898rj41VfJ9LoxUlgZHR0ejd+/eiIyMBABcunQJ8+bNyw4yf/zxBwBgxowZ2eckJCRg\n1qxZCAkJwciRIzFq1KgcX0+fPj3PtnIX7rK3t4e9vb2EfjJCCCm+DGEGZpyegYCHAdj3yz40M2xW\n+ElCIeDhATx8CGHgMbh4Vsfnz0BAAPBfPt4iCQ0NRWhoaPbX8+bNk+w+k8jIyB+GqMSVO5gcOHAA\nJ06cwKZNmwAAO3fuREREBFatWlWidgDaZ0IIKXsOPTiEMUfHwNvOu2iVHEUiYNIk4MIFZB47if6j\ndKChAfj7AwrF3PAh0X0mISEhGDt2bLEvnh+xS14SQkg55lDPAeEjwrH55mY4HXJCUnohEyEKCsCK\nFUC3blDuZIf9K+Px6hUwfjzPyiIL+QaTnTt3YsqUKQgICCj1Rg0NDRHz3ULpmJgYGBkZlXo7hBBS\nVphpmSHcPRxqympotqkZ7r29V/AJAgGwaBHg5ASV7vY4uiEWERHAnDnS6W9u+QYTd3d3HDx4ELq6\nuqXeaNOmTREVFYXo6GhkZGRg79696NOnT6m3QwghZYmKsgo29dmEmW1mwn67PXbe2Vn4SbNnAyNG\nQL2XHU5ufomDB4HlyyXf1x+wfCxatIjZ29uzlJSU/A4pEkdHR6avr88qVarEjIyM2NatWxljjAUF\nBbG6desyU1NTtmjRohK18b0CfiRCCCkz7ry+w8xXmrMx/45haZlphZ+wbBljtWuzV2HPmbExY35+\n4rVX0s/OAifg/f394e/vj5MnT0ovupUQTcATQsqLT2mf4B7ojpefXmL/wP2oVb1WwSesWgUsW4an\nm8+ijXMtbNjAlw4XhcSzBp8+fRqdpF2MuAQomBBCyhPGGJZfXg7fi77w6+uH7ubdCz5h1Spg+XLc\nWRmKTu7GOHAAaNeu8HYknjW4LAWSr3x8fHKsnyaEkLJKIBBgsu1kHBh4AKP+HYW5Z+cWnCzS0xOY\nMAHWXu1x6O8YDBwI3L6d/+GhoaHw8fEpeUfzG//y9fXN/vu+fftyPDdz5swSja1JUgE/EiGElGnx\nSfHM3s+edfHvwt59eVfwwUuXMmZqygLXv2IGBow9fVrw4SX97Mz3zmTPnj3Zf1+0aFGO544fP17y\nKEYIIUQsP6n9hFMup9BIrxFsNtrgSuyV/A+ePBkYNQq9l3fEwglv0KWLZDMNU3EsQggpQ5QUlODb\n2Rd/d/sbvXb3wvpr6/Of65g+HXB0hNuuzhjl8AE9ekguMSQFE0IIKYP6WfbDRfeLWHN1DdyOuCEl\nMyXvA729ge7dMS2kC9pafYSDA5CRUfr9yXc1l6KiIqpWrQoASE1NhYqKSvZzqampyMrKKv3elAJa\nzUUIqUi+ZHzB6KOjcfftXRwcdBCmWqY/HsQY4OUFdu06HLVOQklD9Yc8XhJfGlzWUDAhhFQ0jDGs\nvboW887Nw9a+W9Grbq8fDxKJAHd3CGPj0elLIJrYVsbSpd+epmCSCwUTQkhFdSnmEgYdGIThjYbD\n284bigqKOQ/IygIGDkSGUBFNn/wDt5FKmDyZP0XBJBcKJoSQiuxN8hsMPjAYKsoq2OWwC1oqWjkP\nSE8HevVCspYx6l3cjD//EsDRUQqbFssi2rRICKmo9NT0cHrYaTTQbQCbjTa4EX8j5wGVKwMBAVCL\nvovrXWdhzJhQuLr6lLhdujMhhJByat+9fRgXNA5/df4Lro1ccz75/j3Qpg2iOo5B6/0T8e4dDXPl\nQMGEEEK+uff2Hhz2OaBj7Y5Y0W0FKilW+vbkixdAmzZ47P4HLOY7UzD5HgUTQgjJ6VPaJ7gedsWb\nL29wYOABGFYz/Pbk3btAx44QvH1LcyaEEELyp1FFA4cGH0Lvur3RfHNzXHhx4duTDRsCVwpIy1JE\ndGdCCCEVyIknJzDs8DDMbjsbns09IRAIgP+3d68hUa17HMd/s1VCQyqjIjWyUtSmPGpSUkblhSJO\nxhlqzCBCiygpyCCCQirCF0UURWK7pAuBpW2FGgoxMN90MRwMQqPMnVhpdDuBWOR4OS88e2jS9k6X\nOup8P9CbZ61nzV97WD+ftdY8Szwa3AthAgB/78///ilLkUXzp83X7//+XX4+fjwaDADon9mTZuv+\nlvvq7u7WkgtL1Pi50fAxCRMA8EB+Pn668p8r2vyvzYoviDd8PC5zAYCHszfbFRcUx2WuH/ENeAD4\nNZWVlbKdsxk+DjMTAAA34AEA7keYAAAMI0wAAIYRJgAAwwgTAIBhhAkAwDDCBABgGGECADCMMAEA\nGDYmw4TlVADg11RWVurQoUOGj8NyKgAAllMBALgfYQIAMIwwAQAYRpgAAAwjTAAAhhEmAADDCBMA\ngGGECQDAMMIEGCasyoCxjDABhglhgrGMMBmh3HXiGYrPNXrMgfTvT59f3fdX9vOUwHDHzzlWxmZ/\n+w3WuBvq/zPCZIQiTIz1J0yGFmFirP9YDJMxt9BjaGioGhoa3F0GAIwqc+bM0YsXLwbcf8yFCQBg\n+HGZCwBgGGECADCMMAEAGEaYAAAMI0wAAIZ5RJi8fPlSW7du1fr1691dCuDU1tamzZs3a9u2bSos\nLHR3OYDTQM6ZHhEms2bNUkFBgbvLAFyUlpbKarXq3LlzunnzprvLAZwGcs4cVWGSmZmpadOmaf78\n+S7tZWVlioiIUFhYmI4ePeqm6oD+jdE3b95oxowZkiQvL69hrxWeZajPn6MqTDIyMlRWVubS1tnZ\nqZ07d6qsrEx1dXW6evWqnj59qitXrig7O1vNzc1uqhaeqD9jNDg4WK9evZIkdXV1uaNceJD+jM2B\nGFVhsnTpUk2aNMml7dGjRwoNDVVISIh8fHy0YcMG3bhxQ5s2bdLJkycVGBioT58+afv27Xr8+DEz\nFwyp/oxRi8WikpISZWVlKTU11U0Vw1P0Z2wO5JzpPRRFD6fvLxVIUnBwsKqqqlz2CQgI0NmzZ4e7\nNEDSz8eon5+fLly44MbK4Ol+NjYHcs4cVTOTvphMJneXAPwtxihGqsEcm6M+TIKCgpzXnSXp1atX\nCg4OdmNFgCvGKEaqwRyboz5M4uLiVF9fr8bGRrW3t6uoqIjrzxhRGKMYqQZzbI6qMElPT9fixYv1\n/PlzzZgxQxcvXpS3t7fOnDmjlStXau7cuUpLS1NkZKS7S4WHYoxipBrqscn7TAAAho2qmQkAYGQi\nTAAAhhEmAADDCBMAgGGECQDAMMIEAGAYYQIAMIwwgcfx8vJSTEyM89+xY8fcXZJTcnKympqanLVN\nnz5dwcHBiomJUWxsrDo6Olz2v3z5sjZu3OjS9uHDB02dOlXt7e2yWq16+fLlcP4I8FB8aREex9/f\nX62trYN6zI6ODnl7G1uEu6KiQiUlJcrLy3O2HT58WP7+/tqzZ0+ffVpbWzV79mw1NTXJ19dXknT2\n7FlVV1eroKBAd+7ckc1m0+nTpw3VBvwTZibA/4WEhOjQoUNasGCBoqKi9OzZM0k972rPzMzUokWL\nFBsb63zF7qVLl5SamqqkpCSlpKTo69evslqtMpvNslgsio+Pl91u18WLF5Wdne38nPPnz/cZDoWF\nhVq7dm2v9r/+3rPb7Vq+fLni4uK0atUqvX37Vv7+/lq2bJlsNptz/2vXrik9PV2StHz5ct2+fXvw\nfknATxAm8Dhfv351ucx1/fp1ST3LcU+ZMkV2u107duzQ8ePHJUm5ublKSkpSVVWVKioqtHfvXn35\n8kWSVFNTo5KSEt29e1d5eXmaPHmyamtrdeTIEdntdplMJlmtVtlsNnV2dkrqCaEtW7b0quvevXuK\ni4vr1W4ymdTR0aFdu3appKRE1dXVysjI0IEDByT1rLl07do1SVJzc7Pq6+uVmJgoSfLx8VFQUNCA\n354H/KpR/3IsoL98fX1VU1PT5zaLxSJJio2NVWlpqSSpvLxcNpvNGS7fvn1TU1OTTCaTUlJSNHHi\nREk9YbB7925JktlsVlRUlCRp/PjxSkxMlM1mU0REhBwOh8xmc6/Pbm5uVkBAQK/27u5uPXv2TLW1\ntUpOTpbU87rVwMBASdLq1auVlZWl1tZWFRcXa926dS7vqQgMDFRjYyOLS2JIESbAd8aNGyep5yb9\n9ze7S0tLFRYW5rJvVVWVxo8f79L2s1uQW7duVW5uriIjI5WZmdnvurq7u2U2m3X//v1e23x9fbVq\n1SqVlpaqqKhIJ0+e7NX3t9+4CIGhxQgD/sHKlStdbmD/Nav5MTiWLFmi4uJiSVJdXZ2ePHni3LZw\n4UK9fv1ahYWFzvsZPwoMDNTHjx/73BYeHq7379/r4cOHkiSHw6G6ujrn9vT0dJ04cULv3r1TfHy8\nS9+WlhbNnDnzV39cYEAIE3icH++Z7N+/v9c+JpPJeakoJydHDodDUVFRmjdvng4ePNhrH0nKysrS\n+/fvZTablZOTI7PZrAkTJji3W61WJSQkuLR9LyEhQdXV1X3W4uPjoz/++EP79u1TdHS0YmJi9ODB\nA+c+ycnJamlpUVpamktfh8Oh169fKyIioh+/IaD/eDQYGCRdXV1yOBwaN26cGhoalJKSoufPnzsf\nGV6zZo327NmjFStW9Nm/srJSRUVFys/PH7SaysvLdevWLZ06dWrQjgn0hZkJMEja2tqUkJCg6Oho\nWSwW5efny9vbW58/f1Z4eLj8/Px+GiRSz2O89fX1g/odmIKCApfHkoGhwswEAGAYMxMAgGGECQDA\nMMIEAGAYYQIAMIwwAQAY9j/leS/XaJlfCQAAAABJRU5ErkJggg==\n",
"text": [
"<matplotlib.figure.Figure at 0x10eee0390>"
]
}
],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"flux_ratio = crab_ctools(energy) / crab_flux(energy, reference='meyer')\n",
"plt.plot(energy, flux_ratio)\n",
"plt.xlabel('Energy (TeV)')\n",
"plt.ylabel('Flux ratio (ctools / meyer)');"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEPCAYAAABCyrPIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4zWf+//HnIakGQWpJbbULIoi1Siq2arUxlao2HV+a\n0BqKKtPl234V02HoOrTGWLrQDUVLWk3Vj1iqocSaqK1SS1BBJARN4vz+uEmpJCeJc84ny+txXeca\nydleydU579yf+77ft81ut9sREZESr5TVAUREpHBQQRAREUAFQURErlJBEBERQAVBRESuUkEQERHA\nxQUhIiICX19fAgICsr0/OjqaihUrEhgYSGBgIP/85z9dGUdERHLh4coXDw8PZ+TIkQwcODDHx3Tp\n0oXly5e7MoaIiOSBS0cIQUFB+Pj45PoY7YsTESkcLJ1DsNlsbNy4kZYtW9K7d2/i4+OtjCMiUqK5\n9JKRI61bt+bIkSOULVuWb7/9locffph9+/ZZGUlEpMSytCB4e3tn/fuBBx5g+PDhnDlzhjvuuOOG\nxzVs2JCDBw+6O56ISJHWoEEDDhw4kOfHW3rJ6OTJk1lzCJs3b8Zut99UDAAOHjyI3W7XzW5n/Pjx\nlmcoLDf9LvS70O8i91t+/5B26QghLCyMtWvXkpSURO3atZk4cSLp6ekADB06lMWLFzNz5kw8PDwo\nW7YsCxYscGUcERHJhUsLwueff57r/c888wzPPPOMKyOIiEgeaadyERMcHGx1hEJDv4s/6HfxB/0u\nCs5mt9sL/UYAm81GEYgpIlKo5PezUyMEEREBVBBEROQqFQQREQFUEERE5CoVBBERAVQQRETkKhUE\nEREBVBBEROQqFQQREQFUEERE5CoVBBERAVQQRETkKhUEEREBVBBEROQqS89ULknOnYMdO8ztl1/g\n2DE4cQLS0swN4Pbbza1yZfD1herVoX59aNQIGjeGatWs/RlEpHjTeQguYrfD5s2wdCmsWgX79kHz\n5tCypfmAr1kT7rwTypUDLy+w2eDyZVMcTp+Gkyfh+HE4eNA8d+9e87g2bcytUyfo2BHKl7f6JxWR\nwiq/n50qCE525gzMng1z5oCHB/TvD716Qfv2cNttBX9dux1+/RW2boUtW2DDBti2zRSZ7t2hd2/o\n0MG8p4gIqCBY5uxZmDwZ3n8f+vSBESPMX/I2m+ve8+JFiImB77+HFSvgyBFTfEJD4YEHzOhDREou\nFQQ3u3IF/vtfGD/efBC/+qq5HGSFY8fgm29gyRJTKLp3h8ceMwXKy8uaTCJiHRUENzp0CAYNgowM\nc4nI39/qRH84fRqWL4fPP4effoK+fWHAAAgOhlJaWyZSIqgguElkJAwZAi+8AKNHQ+nSVifKWWKi\nKQzz50NKCgweDOHh1o1kRMQ9VBBczG6HN9+E6dNh0SKz0qeosNshNtaMZhYtMiuVhg838w4aNYgU\nPyoILmS3m9HAmjVmErdWLasTFdyFC7BwIcyYYfZIDB9uRg0+PlYnExFnye9np/4uzCO7HUaNgk2b\nYN26ol0MwKxAiogwS1g/+cSMHOrXh2HDzJ4HESl5VBDy6OWXTTH47juoVMnqNM5js8Hdd5ui8PPP\nZof0vfdCSAisXm0KoYiUDLpklAezZsFbb8HGjVClimUx3ObiRVMg3n7bLFd98UV45BFtehMpagrV\nJaOIiAh8fX0JCAjI9XE//fQTHh4eLF261JVxCmTtWrPHYMWKklEMwBSBp56CuDj4xz/g3XfBz8/s\nt7h0yep0IuIqLi0I4eHhREVF5fqYzMxMXnzxRe6///5CMXF8vePH4YknYN48aNjQ6jTuV6oUPPSQ\naZMxb57Z9FavHrz+OqSmWp1ORJzNpQUhKCgIHwfLVt5991369etH1apVXRkl365cMZvOhgwxyzJL\nus6dzd6LlStND6UGDeC11yA52epkIuIslk4qHzt2jGXLljFs2DDAXO8qLGbMMJu4xo2zOknhEhBg\nNrmtX286sTZsCP/3f6apn4gUbZYWhNGjRzNlypSsiY/Ccsnol19g4kT4+GNNpObEzw8++si0xfjt\nN9PSW4VBpGiz9ONu69atPP744wAkJSXx7bff4unpSZ8+fW567IQJE7L+HRwcTHBwsEsy2e1mk9bz\nz5sPOcldvXqm3ffLL8OkSeZ3Nnw4jBmjTW4i7hYdHU10dHSBn+/yZacJCQmEhISwa9euXB8XHh5O\nSEgIoaGhN93nzmWnS5bAhAlmo5anp1veslg5dMgUhq++gpEj4bnnoEIFq1OJlEyFatlpWFgY99xz\nD3v37qV27dp88MEHzJo1i1mzZrnybQvs8mXTrO7f/1YxKKh69WDuXLOJ75dfzBzD1KmmVYaIFG7a\nmHadt94y+w6WL3f5W5UYe/aYEde6dfDSSzB0qDk3WkRcT83tCiglxfw1Gx0NzZq59K1KpB07zIqt\n7dvNRr9BgzRhL+JqheqSUVEybZrZb6Bi4BotW5qR18KFpi1G8+bwxRdmv4eIFA4aIWB23davb3oV\naWWR69nt5hzol182BWHyZFOMC9E2FJFiQZeMCuDtt80k6MKFLnsLyYbdblZ1vfIK1KgBU6ZAhw5W\npxIpPlQQ8ik93ayMiYyEwECXvIU4kJEBH35oNgN26GCWrTZpYnUqkaJPcwj5tHgxNG6sYmAlDw/T\nXXXfPmjfHoKCzNfHjlmdTKRkKfEFYfp0cxKaWK9sWXP2wr59cMcd0KKFmWdQAz0R98hzQbh06RKX\nL192ZRa3i42FxERzOpgUHj4+ZjPb9u1w8qQZwb3zjtk4KCKuk2NBuHLlCkuXLuXRRx+lZs2a1KtX\njzp16lCzZk369evHl19+WWia0RXUnDmmvXXp0lYnkezUrg3vv2+O8ly92jTU++QTLVUVcZUcJ5Xv\nvfdegoKC6NOnD61ataJMmTIAXL58mW3btrF8+XI2bNjAunXrXB/SBZPKFy6YD5ydO6FWLae+tLjI\nunWmtcjly+aQnp49rU4kUrg5bZXR5cuXs4pATvLyGGdwRUH4+GPT13/FCqe+rLjYtaWq//u/Zu/I\n1KnQqpXVqUQKJ6etMipTpgwZGRk0yWX9nzuKgat8/LFpnyBFi80G/fpBfDz06QP33w8DB8Kvv1qd\nTKToy3VS2cPDAz8/P34tZv9vS0yELVvMB4oUTZ6e8MwzZkVSnTrQurW5nHT2rNXJRIouh6uMzpw5\ng7+/P926dSMkJISQkJBsD7ApSr74whQDLy+rk8itqlDBnO28a5cpBn5+WpEkUlAOdypnd/qOzWaj\nS5cursqU7fs5cw6hUyfTLqF3b6e9pBQScXGmzXZcnOmR1L8/lCrxu22kpHJJ64qEhAQOHDhAjx49\nSEtLIyMjgwpuPAbLmQXhyBEzCXn8ONx2m1NeUgqhNWvMMailSsEbb4Ab/34RKTSc3rpi9uzZPPro\nowwdOhSAo0eP0rdv34IntNiyZfDQQyoGxV3XrrB5sznC88knzSXC+HirU4kUbg4LwowZM9iwYUPW\niKBx48b89ttvLg/mKsuWwcMPW51C3KFUKQgLg59/huBgcxs6FE6csDqZSOHksCCUKVPmhuWlGRkZ\n2Ipo4/rkZNPm+r77rE4i7lSmDIwZA3v3grc3+Pubzqrnz1udTKRwcVgQunTpwqRJk0hLS+P777/n\n0UcfJaSINv9ZudJ00ixXzuokYgUfH3jzTbPkeN8+0yNp9mzTfltE8jCpfOXKFebOncvKlSsB6NWr\nF0OGDHHrKMFZk8pPPmnaKw8ffuuZpOjbssVMPJ88aXY8P/SQTm2T4sXpq4wiIyN58MEHKWXh2j1n\nFIQrV6B6dYiJMQfiiIBphbFihdnUVrWqWZHUrp3VqUScw+mrjBYsWEDDhg154YUX+Pnnn28pnJV2\n7IBKlVQM5EY2Gzz4oPnvY8AAs+AgLAx++cXqZCLu57AgfPrpp2zbto369evz5JNP0rFjR2bPnk1q\naqo78jnN999rMlly5uFhWqHv2wfNmplLi2PGwJkzVicTcZ88XQeqWLEi/fr147HHHiMxMZEvv/yS\nwMBApk+f7up8TrNypdoli2PlysG4cWan88WLphXGG2/ApUtWJxNxPYdzCMuWLeOjjz5i//79DBw4\nkCeffJJq1aqRlpZGs2bNSEhIcH3IW5xDuHwZKlc2Z/RWrOjEYFLs7d1rWm1v3Wp6Jg0YoFYYUnQ4\nfVJ50KBBDB48mHvvvfem+1atWkWPHj3ynzKfbrUgrFsHY8fCTz85MZSUKD/8YFYkpaWZEYNGm1IU\nqJdRNv7xD7MJ6fXXnRhKShy7HZYuNc3z6tUz/z3pcB4pzNTLKBtr10I2AxyRfLHZ4JFHTE+kv/zF\nHM7zP/+jw3mk+HBpL6OIiAh8fX0JCAjI9v5ly5bRsmVLAgMDadOmDatXr85H9LxJTzdNzjp1cvpL\nSwl1/eE89eqZw3nGjtWKJCn6XNrLKDw8nKioqBzv79GjBzt27GDbtm189NFHPP3003l63fzYts38\nn9bHx+kvLSVchQrmcuTu3XDhglmRNHWqWZ0kUhS5tJdRUFAQPrl8Epe7rqnQ+fPnqVKlSp5eNz9+\n+EGjA3Gt6tXhv/+FDRvMaNTPDz78EDIzrU4mkj8OC8KUKVOoWrUqAQEBzJo1i969e/PPf/7TaQG+\n+uormjZtygMPPOCSfQ0//ggdOzr9ZUVu4ucHS5bAggXw/vtmwvmbb8xktEhRkKdVRrciISGBkJAQ\ndu3alevj1q9fz5AhQ9i7d+9N99lsNsaPH5/1dXBwMMHBwXl6/7vugv/3/6BRo3zFFrkldjtERpo9\nDFWqmEtJd99tdSop7qKjo2849njixInOb2736quvkpCQQMbVPsE2m42UlJQ8vUFeCwJAgwYN2Lx5\nM5UrV74xZAGXnSYmQkAAJCWpi6VYIyMD5s+H8eNN07zJk6FJE6tTSUnh9GWno0ePZt68eZw+fZrU\n1FRSU1PzXAwcOXjwYFbY2NhYgJuKwa3YvNn0pFExEKt4eEBEhFmR1LGjOY/j6afNrnmRwsbD0QNq\n1aqFv79/gdpfh4WFsXbtWpKSkqhduzYTJ04kPT0dgKFDh7JkyRLmz5+Pp6cn5cuXZ8GCBfn/CXLx\n00+mIIhYzcvL7HQeMgSmTIEWLUxhePFF04VXpDBweMkoJiaGV199la5du3Lb1ZPpbTYbY8aMcUvA\na+9XkEtG990Ho0aZg09ECpMjR2DCBDPP8PzzMGKEKRoizuT0S0bjxo2jfPnyXLp0ifPnz3P+/Pki\n0frabjcNydq0sTqJyM1q1zYrkdauNYc2NWoEc+boOE+xlsMRQvPmzdm9e7e78mSrICOEX38112wT\nE10USsSJNm0yK5KOHYNJk0yLDM19ya1y+gihd+/efPfdd7cUygqxsRAYaHUKkbzp0MEsj373XbMS\nqX1787WIOzkcIZQvX560tDRuu+02PD09zZPysezUGQoyQnj1VXOOshP30Im4xZUrsHgxvPIK1KkD\n//qXznmWgnH6COH8+fNcuXKFS5cuOX3ZqSvt2KHWxFI0lSoF/fubrqqPPgp9+0JoqPlaxJWK7dlP\nO3ZAy5ZWpxApOE9PGDoU9u83/bi6doVBg+DQIauTSXFVLAtCcjKcPg0NGlidROTWeXmZ9tr790Pd\nutC2rVmmevy41cmkuCmWBWH3bmjWTGffSvFSoQJMnAg//wxlykDz5mZl0tmzVieT4iLHj8w2bdrw\n7LPPEhUVxaVLl9yZ6Zbt3m16GIkUR1WrwltvwfbtZiTcuLFZPFEEtgdJIZdjQYiJieHhhx9mzZo1\ndOnShQceeIBp06axb98+d+YrkLg489eTSHFWuzbMng0bN8KePWZz29tv64AeKbg8t78+duwYUVFR\nfPfddxw4cIC7776b//znP67OB+R/6VTXrvDyy9CzpwtDiRQyu3aZrqqbNpklq0OGwNVuM1JC5fez\ns0DnIWRmZhITE0MnNx1Flt8fytfXHJ1Zo4YLQ4kUUlu2wLhxZq7h1Vfhf/7HdF2VksctBcHd8vND\nnT4N9eublUba+i8l2fr18H//BydOmMno/v210KKkcfrGtKJmzx6zwkjFQEq6oCCIjob33oN33jEb\nNb/8Ukd6Ss6KZUHQiVQihs1m5tJiYkyPpNdeMx2AIyNVGORmDgvC888/T0pKCunp6XTv3p0qVarw\n8ccfuyNbgfz8swqCyJ/ZbOZckK1bzfzCK6+YhnpRUSoM8geHBWHlypVUqFCBr7/+mrp163Lw4EHe\neOMNd2QrkL17wc/P6hQihZPNZnojbd9uDuYZM8a0xVi1SoVB8lAQMq6e2PH111/Tr18/KlasiK0Q\nX6BXQRBxrFQp0zhv1y4YORKeeQa6dDFzDlJyOSwIISEhNGnShK1bt9K9e3d+++03br/9dndky7f0\ndDh82KwyEhHHSpeGsDCzmXPIEHPr1g3WrbM6mVghT8tOT58+TaVKlShdujQXLlwgNTWVO++80x35\ngLwvndq/H3r1gl9+cUMokWIoPR0++cS0wqhTx2x069LF6lRSUPlddprjdpUlS5bcdGno2gvbbDZC\nQ0MLGNF19u832/dFpGA8PSE8HAYMgE8/hYgIuOsumDBBhaEkyLEgREZG5jpXUFgLQsOGVqcQKfo8\nPeHJJ+GvfzWFYfBg0ztp/HgIDrY6nbhKsdqpPGoU1KsHzz3nhlAiJUhGhikMr70GtWqZEYMKQ+Hn\n9J3KycnJPPfcc7Rp04Y2bdowduxYzp07d0shXeXgQR2KI+IKHh7mtLaffzaXlIYMMZeQVq/WctXi\nxGFBiIiIoEKFCnzxxRcsWrQIb29vwsPD3ZEt31QQRFzr+sIQEQHDhkHnztrgVlw4vGTUsmVLduzY\n4fB7rpSXYc+VK1C2LJw5Y/5XRFwvMxMWLTKrksqWNbugQ0LUS6ywcPolIy8vL9avX5/19YYNGyhb\nCD9xExPBx0fFQMSdru1j2LULXnrJtNsODITFi80faVK0OBwhbN++nYEDB2bNG/j4+DBv3jxatmzp\nloCQtyq3fj28+KI5PUpErGG3w9dfm8nntDTTM6l/f1M4xP1cdh5CSkoKABUqVMjzi0dERPDNN99Q\nrVo1du3addP9n376Ka+//jp2ux1vb29mzpxJixYtbg6Zhx/q44/h22/hs8/yHE9EXMRuh5Ur4R//\ngKQkc4LhE0+Y5aziPi5bZdS1a1e6du2ar1VG4eHhREVF5Xh//fr1WbduHTt37mTcuHE8/fTTeQ7+\nZwkJZsmpiFjPZjNdAzZsgJkzYd48aNzY/PvSJavTSU5cusooKCgIHx+fHO/v2LEjFStWBKBDhw4c\nPXo0j7FvlpBgttqLSOFhs5neSKtXm9H7ihWm19ibb0JqqtXp5M8cFoSDBw8yceJE6tevT4MGDZgw\nYQIHDx50epD333+f3r17F/j5v/4Kdes6L4+IOFfHjuZgnm+/Nec+169vNridOWN1MrnG4dHb11YZ\nBQUFAa5ZZbRmzRo++OADfvjhhxwfM2HChKx/BwcHE/ynbZK//qoRgkhR0LIlLFgA+/bB1Kmm3czg\nweZshurVrU5XtEVHRxN9Cz3MXb7KKCEhgZCQkGwnlQF27txJaGgoUVFRNMyhEZGjiRHtQRApug4f\nNpeQPvkEHn8cXnhBo31ncfqkcoUKFdi5c2fWbfv27Xh7e99SyGsOHz5MaGgon3zySY7FIC9OnYIK\nFVQMRIqiu+6C6dPN7ueKFc2ZzwMHwu7dVicreRyOEAIDA9m2bdsN32vTpg1bt251+OJhYWGsXbuW\npKQkfH19mThxIunp6QAMHTqUIUOG8OWXX3LXXXcB4OnpyebNm28O6aDK/fQT/O1v5rxYESnakpPN\naqRp06BdO7PhrVMnq1MVTU7bh7Bnzx7i4+N5/vnnefPNN7Hb7dhsNlJSUnjjjTeIi4tzWmiHIR38\nUEuWmH0IX33ltkgi4mIXL8JHH5nLSdWrm8LQu7c5/lPyxmkH5Ozbt4/IyEjOnTtHZGRk1ve9vb2Z\nM2fOraV0sqNHzbBTRIoPLy/TPO+pp8wffePGmaLw4otmrkGb3JzP4SWjH3/8kY4dO7orT7YcVbm/\n/x2qVTOTUSJSPNnt8P33MGWK6Ww8dqxZnVSunNXJCi+nTyrPnDmT5OTkrK/Pnj1LREREwdK5yNGj\n5tAOESm+bDa47z6zye2LL2DtWtOdYOJEOH3a6nTFg8OCsHPnTipVqpT1tY+PD7GxsS4NlV9Hj0LN\nmlanEBF3ad/eXEZavx6OHDFnqY8YAb/8YnWyos1hQbDb7Zy5bivhmTNnyMzMdGmo/Dp2TCMEkZLI\nzw/mzoW4OLP0vH17ePRRyGaxouSBwzmE+fPnM2nSJPr374/dbueLL77glVdeYeDAge7KmOt1MLsd\nbr/dLFXz8nJbJBEphM6fh/ffh3feMQtNnn8eHnyw5K5Mckn767i4OFavXo3NZqNbt240a9bslkLm\nV24/VFKS6aKofigick1Ghjmk5403zLkMY8fCgAHmj8eSxGkFITU11eGO5Lw8xhly+6F27jR91rWr\nUUT+zG6H6GizlyE21swzDBsGd9xhdTL3cNoqo759+/LMM8+wcuXKm+YQvvvuO4YNG0bfvn1vLa0T\nJCZCjRpWpxCRwshmg65d4ZtvYNUqs1y1YUMYNUoT0NnJsSCsWrWKRx55hEWLFtGpUycqVqxIxYoV\n6dSpE4sXL+axxx5j1apV7syarcREdUgUEcf8/eGDD8zVhHLlzAR0aKhZqZSPP6KLtTwfoWml3IY9\nkyaZiaR//cvNoUSkSLtwwZzk9u9/mxVKY8aYFUrFaQe00zemFXbHj2uEICL5V64cDB9uuqyOH2+W\nr9arZ3ZCl9RFKkW+IJw4oYIgIgVXqhSEhJgd0JGRsGcPNGgAzzxjDvEpSYpFQbjzTqtTiEhxEBho\nLiPFxYGPj2m73acPrFlTMuYZ8jSHsH37dtavX4/NZiMoKCjPp6U5S27XwRo2NGe0Nmrk1kgiUgKk\npZnW+u+8Yza+jh5tOq2WKWN1srxx+hzCtGnTGDBgAKdOneLkyZMMGDCA6dOn31JIZzpxAnx9rU4h\nIsVR2bIwdCjEx5sFLJ99Zs5uf/VVs8KxuHE4QggICCAmJoZyV3vMXrhwgbvvvjvHM5JdIacqd+EC\nVKliqrjN5rY4IlKC7dkD774Ln39uDuwZNQo6dLA6VfZcssqo1HWNQEoVoqYgJ0+a0YGKgYi4S9Om\n8J//wKFD5vznsDBTED79FH7/3ep0tybHE9OuCQ8Pp0OHDoSGhmK32/nqq68KzXkI1wqCiIi7Vapk\n9i48+yx8/TVMn26a6f3tb+YyU1H8bMrTpPLWrVvZsGFD1qRyYGCgO7JlyWnYs2yZWTt83QmfIiKW\n2b3bXE5atMisTho1yowirOK05nYpKSlUqFAhq4/RtYfZrl6fucON3aFy+qHmzIFNm0xREBEpLM6c\nMZ9LM2aYs1pGjIBHHoHbbnNvDqcVhAcffJBvvvmGunXrZhWB6x06dKjgKfMppx9q0iQzsTx5stui\niIjkWUaGuZIxY4aZjB4yxFxOcteBXi45D8FqOf1Qzz5rtpqPHm1BKBGRfIiPN5PRn30G3bqZndDB\nwa5dFOP0VUbdu3fP0/es8NtvULWq1SlERBxr1gzeew9+/dUUhBEjoHlzM3pISbE6nZFjQbh48SKn\nT5/m1KlTnDlzJuuWkJDAsWPH3JkxR6dOqSCISNHi7W2a6u3ebYpBdDTUrWtGDHFx1mbLcdnprFmz\nmDZtGomJibS5bprc29ubESNGuCWcI0lJKggiUjTZbOaSUXAwHDsGs2dDz57g52eKw1/+4v5W3A7n\nEKZPn86oUaPclSdbOV0Hq1nTrDJy1wSNiIgrpafDl1+akcPBg/D00zB4sPmsKwiXTCrv3r2b+Ph4\nLl26lPW9gQMHFixhAWT3Q9nt5sDs5GTTdEpEpDjZtctMQi9caEYRQ4eaEUR+mkU4fVJ5woQJjBw5\nkhEjRrBmzRpeeOEFli9fnqcXj4iIwNfXl4CAgGzv//nnn+nYsSO33347b731Vp5DgzklzdNTxUBE\niqeAAJg5Ew4fhgcegP/9X9PVeepUs6DGFRwWhMWLF7Nq1SqqV6/Ohx9+yI4dO0hOTs7Ti4eHhxMV\nFZXj/ZUrV+bdd9/l73//e94TX5WUZBrbiYgUZ+XLw1NPwdatsGCBObTHz8+04Y6Odu45DQ4LgpeX\nF6VLl8bDw4Nz585RrVo1jhw5kqcXDwoKwsfHJ8f7q1atStu2bfEswMyJCoKIlCQ2G7RrB++/bxrr\nde5slq42bWrOa3DGsZ8OC0K7du04e/YsTz31FG3btiUwMJB77rnn1t/5FiUlQeXKVqcQEXG/SpVM\nMdi1y7TIiI2F+vVh0CDYuLHgo4Zcu53a7XZeeuklfHx8+Nvf/kavXr1ISUlx+4lpYOYyrgkODub0\n6WAVBBEp0Ww2M1Lo3Nn8kTxuXDQhIdF4eEDbtvl/PYftr3v37s3u3bsBqFevXv7fwUmuLwhgWs2q\nIIiIGFWqwMyZwcyYEUx0tGn+CRPz9Rq5XjKy2Wy0adOGzZs3FzxlHhSkndLp0yoIIiJ/VqqUaY3x\n+ef5f67DfQh+fn4cOHCAOnXqZB2jabPZ2Llzp8MXDwsLY+3atSQlJeHr68vEiRNJT08HYOjQoZw4\ncYJ27dqRkpJCqVKl8Pb2Jj4+nvLly98YMpu1tCNHmiVYFu+ZExEptJy+MS0hISHb79etWzc/uW5J\ndj/UE0/Agw/CX//qthgiIkVKfguCwzkEd37w58eZM+DGM3pERIq9fGyCLlxUEEREnKtIF4Rc9ryJ\niEg+OSwI8fHxN30vOjraFVny5exZFQQREWdyWBD69+/P1KlTsdvtpKWlMXLkSF566SV3ZMvRlStw\n7pzZrSciIs7hsCBs2rSJI0eO0LFjR9q3b0/16tXZuHGjO7LlKCUFypZ1/+ERIiLFmcOC4OHhgZeX\nFxcvXuTSpUvUr1+fUvlpyO0CyckaHYiIOJvDT/b27dtz++23s2XLFtavX89nn33Go48+6o5sOdL8\ngYiI8zmQGfR9AAAOLklEQVTchzB37lzatWsHQPXq1Vm+fDnz5893ebDcaIQgIuJ8DguCr68vhw8f\nvuF7Xbp0cVmgvEhO1ghBRMTZ8tTt1GazAXDp0iUOHTqEn58fcXFxLg+Xk+RkqFjRsrcXESmWHBaE\na62vr4mNjWXGjBkuC5QXGiGIiDhfvpcLtW7dmk2bNrkiS55phCAi4nwORwhvvfVW1r+vXLlCbGws\nNWvWdGkoR86dg7vusjSCiEix47AgpKamZs0heHh48NBDD/HII4+4PFhukpOhRQtLI4iIFDsOC8Kf\nj64sDM6d0yUjERFny7EghISE5Pgkm83G8uXLXRIoL9THSETE+XIsCGPHjs3xSdcuIVlFIwQREefL\nsSDUq1ePOnXquDNLnqWkqCCIiDhbjstOH3744ax/Wz2J/GcpKVChgtUpRESKlzztQ/jll19cnSNf\nzp1TQRARcbYid4Tm779DZibcfrvVSUREiheb3W63Z3dH6dKlKVu2LAAXL17Ey8vrjyfZbKSkpLgn\n4dX3uxYzKQn8/OD0abe9vYhIkXT9Z2de5DipnJmZ6ZRAzqb5AxER1yhyl4xUEEREXKPIFYTUVPD2\ntjqFiEjxU+QKQkqKCoKIiCsUuYKgEYKIiGu4tCBERETg6+tLQEBAjo8ZNWoUjRo1omXLlmzbts3h\na6ogiIi4hksLQnh4OFFRUTnev2LFCg4cOMD+/fuZPXs2w4YNc/iaKggiIq7h0oIQFBSETy5nXS5f\nvpxBgwYB0KFDB5KTkzl58mSur3n+vAqCiIgrWDqHcOzYMWrXrp31da1atTh69Giuz9EIQUTENSyf\nVP7zLjpHrbU1QhARcQ2HJ6a5Us2aNTly5EjW10ePHs3xvOZrJ7f98AP4+AQDwS7PJyJSlERHRxMd\nHV3g5+fYy8hZEhISCAkJYdeuXTfdt2LFCt577z1WrFhBTEwMo0ePJiYm5uaQ1/Xj6NsXBg2C67pz\ni4hINpzWy8gZwsLCWLt2LUlJSdSuXZuJEyeSnp4OwNChQ+nduzcrVqygYcOGlCtXjg8//NDha54/\nD+XKuTK1iEjJ5PIRgjNcX+U6doS33zb/KyIiOcvvCMHySeX80ghBRMQ1imRBKF/e6hQiIsWPCoKI\niABFsCBcuKBLRiIirlCkCkJmJly6BNed5ikiIk5SpApCWhqULQulilRqEZGioUh9tOpykYiI66gg\niIgIoIIgIiJXqSCIiAhQxArCtUllERFxviJVEDRCEBFxnSJVEKpUgXvvtTqFiEjxVOS6nYqISN4U\n+26nIiLiGioIIiICqCCIiMhVKggiIgKoIIiIyFUqCCIiAqggiIjIVSoIIiICqCCIiMhVKggiIgKo\nIIiIyFUqCCIiAqggiIjIVSoIIiICuLggREVF0aRJExo1asTUqVNvuv/s2bP07duXli1b0qFDB+Li\n4lwZR0REcuGygpCZmcmIESOIiooiPj6ezz//nD179tzwmMmTJ9O6dWt27NjB/PnzefbZZ10Vp9iI\njo62OkKhod/FH/S7+IN+FwXnsoKwefNmGjZsSN26dfH09OTxxx9n2bJlNzxmz549dO3aFQA/Pz8S\nEhI4deqUqyIVC/qP/Q/6XfxBv4s/6HdRcC4rCMeOHaN27dpZX9eqVYtjx47d8JiWLVuydOlSwBSQ\nX3/9laNHj7oqkoiI5MJlBcFmszl8zEsvvURycjKBgYG89957BAYGUrp0aVdFEhGR3Nhd5Mcff7T3\n6tUr6+vJkyfbp0yZkutz6tata09NTb3p+w0aNLADuummm2665ePWoEGDfH1ue+Aibdu2Zf/+/SQk\nJFCjRg0WLlzI559/fsNjzp07h5eXF7fddhtz5syhS5culC9f/qbXOnDggKtiiojIVS4rCB4eHrz3\n3nv06tWLzMxMBg8eTNOmTZk1axYAQ4cOJT4+nieffBKbzUbz5s15//33XRVHREQcsNntdrvVIURE\nxHqFeqeyo41tJcWRI0fo2rUr/v7+NG/enOnTp1sdyXKZmZkEBgYSEhJidRRLJScn069fP5o2bUqz\nZs2IiYmxOpJl/vWvf+Hv709AQABPPPEEly9ftjqS20RERODr60tAQEDW986cOUPPnj1p3Lgx9913\nH8nJyQ5fp9AWhLxsbCspPD09eeedd4iLiyMmJoYZM2aU2N/FNdOmTaNZs2Z5Ws1WnD377LP07t2b\nPXv2sHPnTpo2bWp1JEskJCQwZ84cYmNj2bVrF5mZmSxYsMDqWG4THh5OVFTUDd+bMmUKPXv2ZN++\nfXTv3p0pU6Y4fJ1CWxDysrGtpLjzzjtp1aoVAOXLl6dp06YkJiZanMo6R48eZcWKFQwZMoSSfMXz\n3LlzrF+/noiICMDM21WsWNHiVNaoUKECnp6epKWlkZGRQVpaGjVr1rQ6ltsEBQXh4+Nzw/eWL1/O\noEGDABg0aBBfffWVw9cptAUhLxvbSqKEhAS2bdtGhw4drI5imeeee4433niDUqUK7X++bnHo0CGq\nVq1KeHg4rVu35qmnniItLc3qWJa44447GDt2LHfddRc1atSgUqVK9OjRw+pYljp58iS+vr4A+Pr6\ncvLkSYfPKbT/jyrplwKyc/78efr168e0adOyXZ5bEnz99ddUq1aNwMDAEj06AMjIyCA2Npbhw4cT\nGxtLuXLl8nRZoDg6ePAg//73v0lISCAxMZHz58/z6aefWh2r0LDZbHn6TC20BaFmzZocOXIk6+sj\nR45Qq1YtCxNZKz09nUceeYQBAwbw8MMPWx3HMhs3bmT58uXUq1ePsLAwVq9ezcCBA62OZYlatWpR\nq1Yt2rVrB0C/fv2IjY21OJU1tmzZwj333EPlypXx8PAgNDSUjRs3Wh3LUr6+vpw4cQKA48ePU61a\nNYfPKbQF4fqNbb///jsLFy6kT58+VseyhN1uZ/DgwTRr1ozRo0dbHcdSkydP5siRIxw6dIgFCxbQ\nrVs35s+fb3UsS9x5553Url2bffv2AbBq1Sr8/f0tTmWNJk2aEBMTw8WLF7Hb7axatYpmzZpZHctS\nffr0Yd68eQDMmzcvb39I5mtfs5utWLHC3rhxY3uDBg3skydPtjqOZdavX2+32Wz2li1b2lu1amVv\n1aqV/dtvv7U6luWio6PtISEhVsew1Pbt2+1t27a1t2jRwt63b197cnKy1ZEsM3XqVHuzZs3szZs3\ntw8cOND++++/Wx3JbR5//HF79erV7Z6envZatWrZP/jgA/vp06ft3bt3tzdq1Mjes2dP+9mzZx2+\njjamiYgIUIgvGYmIiHupIIiICKCCICIiV6kgiIgIoIIgIiJXqSCIiAiggiDFROnSpQkMDMy6vf76\n61ZHytKjRw8OHz6cla169erUqlWLwMBAWrduTUZGxg2PnzdvHk888cQN30tKSqJatWr8/vvv9O/f\nn0OHDrnzR5ASQvsQpFjw9vYmNTXVqa+ZkZGBh8etHSq4evVqlixZwowZM7K+N3HiRLy9vRkzZky2\nz0lNTaV+/focPnwYLy8vAP773/+yZcsW5s6dy/fff09kZKTOxRCn0whBirW6desyYcIE2rRpQ4sW\nLdi7dy8AFy5cICIigg4dOtC6dWuWL18OwEcffUSfPn3o3r07PXv25OLFi/Tv3x9/f39CQ0O5++67\n2bp1Kx9++CHPPfdc1vvMmTMn2w/4zz77jL/85S83ff/a32Fbt24lODiYtm3bcv/993PixAm8vb3p\n0qULkZGRWY9fsGABYWFhAAQHB7NixQrn/ZJErlJBkGLh4sWLN1wy+uKLLwDT5bFq1aps3bqVYcOG\n8eabbwIwadIkunfvzqZNm1i9ejXPP/98Vuvobdu2sWTJEtasWcOMGTOoXLkycXFxvPbaa2zduhWb\nzUb//v2JjIwkMzMTMIVk8ODBN+X64YcfaNu27U3ft9lsZGRkMHLkSJYsWcKWLVsIDw/nlVdeASAs\nLCzrgJfExET2799Pt27dAHNgUs2aNUv8IUnifLc2HhYpJLy8vNi2bVu294WGhgLQunVrli5dCsDK\nlSuJjIzMKhCXL1/m8OHD2Gw2evbsSaVKlQDzgX6toaC/vz8tWrQAoFy5cnTr1o3IyEiaNGlCenp6\nto3lEhMTueOOO276vt1uZ+/evcTFxWX17c/MzKRGjRoA9O7dm+HDh5OamsqiRYvo16/fDe2La9So\nQUJCQok9IU1cQwVBir0yZcoAZuL5+gncpUuX0qhRoxseu2nTJsqVK3fD93KaZhsyZAiTJk2iadOm\nWaeW5Yfdbsff3z/bNs1eXl7cf//9LF26lIULF/LOO+/c9NySfkCQOJ/+i5ISqVevXjdMyl4bXfz5\nw79Tp04sWrQIgPj4eHbt2pV1X/v27Tl69CifffZZ1vX9P6tRowanT5/O9j4/Pz9OnTpFTEwMYM68\niI+Pz7o/LCyMt99+m99++4277777huceP36cOnXq5PXHFckTFQQpFv48h/Dyyy/f9JjrT40aN24c\n6enptGjRgubNmzN+/PibHgMwfPhwTp06hb+/P+PGjcPf3/+Gc4v79+9P586dczzLuHPnzmzZsiXb\nLJ6enixevJgXX3yRVq1aERgYyI8//pj1mB49enD8+HEee+yxG56bnp7O0aNHadKkST5+QyKOadmp\nSC6uXLlCeno6ZcqU4eDBg/Ts2ZN9+/ZlLUcNCQlhzJgxdO3aNdvnR0dHs3DhQmbOnOm0TCtXruSb\nb75h2rRpTntNEdAIQSRXFy5coHPnzrRq1YrQ0FBmzpyJh4cHycnJ+Pn5UbZs2RyLAZglovv373fq\nHom5c+fesORVxFk0QhAREUAjBBERuUoFQUREABUEERG5SgVBREQAFQQREblKBUFERAD4/zvoHFwo\nVjziAAAAAElFTkSuQmCC\n",
"text": [
"<matplotlib.figure.Figure at 0x10fa21210>"
]
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment