{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Sherpa Quick Start\n",
    "==========\n",
    "\n",
    "This tutorial shows some basic Sherpa features.\n",
    "\n",
    "Workflow:\n",
    " - create synthetic data: a parabola with noise and error bars\n",
    " - load data in Sherpa\n",
    " - plot data using matplotlib\n",
    " - set, inspect, edit a model to fit the data\n",
    " - fit the data\n",
    " - compute the confidence intervals for the parameters\n",
    " - explore the parameter space"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First of all, let's activate the inline matplotlib mode. Sherpa seamlessly uses matplotlib to provide immediate visual feedback to the user. Support for matplotlib in Sherpa requires the matplotlib package to be installed."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The following commands just avoid some unnecessary logging duplication when using Sherpa:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import logging\n",
    "logging.getLogger('sherpa').propagate = 0"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, let's create a simple synthetic dataset, using numpy: a parabola between x=-5 and x=5, with some randomly generated noise"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "\n",
    "x = np.arange(-5, 5.1)\n",
    "y = x*x + 23.2 + np.random.normal(size=x.size)\n",
    "e = np.ones(x.size)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's import Sherpa:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "WARNING: failed to import sherpa.astro.xspec; XSPEC models will not be available\n"
     ]
    }
   ],
   "source": [
    "from sherpa.astro import ui as sherpa"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Depending on how you installed Sherpa, certain special features may be enabled or disabled. Sherpa prints some warning messages when it cannot find some of its modules, as shown above. These warnings are benign. You can refer to the Sherpa documentation to find out what additional features you can enable in Sherpa and how to enable them."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's load and plot the data we just created. Notice we are assigning the ID `mydata` to the dataset we are loading. We will use this ID to refer to the same dataset in the rest of the tutorial. Sherpa can deal with multiple datasets, fit them simultaneously with the same model, and even link parameters between models. Sherpa can read ASCII table and FITS files (provided the `pyfits` package is installed)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEPCAYAAABFpK+YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAEd5JREFUeJzt3X+sZGV9x/H3h10R0epCbdatUKGJxN+iFULV1usPyEIU\n",
       "qX/4g1TRWkNsRUIqFTDWNSa1WkWbNjFQsFmtRgkqBQOVFbka22oXu4CChtaAkQYWaqHVYKuWb/+4\n",
       "Z3cv9wd7L8yZMzPP+5WQOefMmT3f4d77mWee8zznpKqQJLXjgKELkCSNl8EvSY0x+CWpMQa/JDXG\n",
       "4Jekxhj8ktSYjX0fIMltwH8D/wf8vKqOTXIo8FngScBtwKur6t6+a5EkjafFX8BcVT2nqo7ttp0D\n",
       "7Kiqo4BrunVJ0hiMq6snS9ZPBrZ3y9uBU8ZUhyQ1b1wt/i8nuS7JW7ptm6tqd7e8G9g8hjokSYyh\n",
       "jx94QVXdkeRXgB1Jvrf4yaqqJF43QpLGpPfgr6o7use7k3wBOBbYneQJVXVnki3AXUtf54eBJD00\n",
       "VbW0e/0Beg3+JAcDG6rqx0keDZwAvBe4HDgN+ED3eNlKr99f8dMsybaq2jZ0HX2Y5fcGvr9p18D7\n",
       "22+jue8W/2bgC0n2HOtTVXV1kuuAS5K8mW44Z891SNLYJMwBc93qHDDfLc9X7V0eTK/BX1W3Akev\n",
       "sP0/gZf1eWxJGkoX7vMACVW190NgIjhzdzjzQxfQo/mhC+jZ/NAF9Gx+6AJ6Nj90AUPLpN6IJUnN\n",
       "ch+/pDZ0Lf6xZdlastMWvyQ1xuCXpMYY/JLUmHHM3J1Ykz7kStJ0S7iwe7wSOLWKibgKsSd39x5v\n",
       "vCdgJM2+hHngRd3qJVW8pv9jenJXkoZ0X/e4Ezh9yEIWs8W/93i2+CWNVsIm4B7gkHF186wlOw3+\n",
       "vccz+CWNnuP4JUmDM/h54Jn37quZJM0sg3/BUd3jicAFQxYiSX1rehz/IhN55l3SdFoyR+irCdu6\n",
       "5YmYI+TJXYY58y5JfXBUz7qO56geSdPPUT2SpGUMfklqjMEvSY0x+CWpMU2f3PWyzJJmjaN6JKkx\n",
       "juqRJC1j8EtSYwx+SWqMwS9JjTH4JakxBr8kNab34E+yIcmuJFd069uS3N5t25Vka981SJL2Gcf1\n",
       "+M8EbgZ+qVsv4PyqOn8Mx5YkLdFriz/JYcBJwEWw95LHWbQsSRqzvrt6PgKcDdy/aFsBZyS5IcnF\n",
       "SbzHrSSNUW/Bn+TlwF1VtYsHtvA/BhwJHA3cAXy4rxokScv12cf/fODkJCcBBwGPTfKJqnrDnh2S\n",
       "XARcsdo/kGTbotX5qprvqVZJmkpJ5th3scm1vWYcF2lL8iLgHVX1iiRbquqObvtZwDFVdeoKr/Ei\n",
       "bZK0TmvJznGM6oGFrp49nzAfTPLsbv1W4PQx1SBJwssyS9JM8bLMkqRlDH5JaozBL0mNMfglqTEG\n",
       "vyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BL\n",
       "UmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1\n",
       "pvfgT7Ihya4kV3TrhybZkeSWJFcn2dR3DZKkfcbR4j8TuBmobv0cYEdVHQVc061Lksak1+BPchhw\n",
       "EnARkG7zycD2bnk7cEqfNUiSHqjvFv9HgLOB+xdt21xVu7vl3cDmnmuQJC2ysa9/OMnLgbuqaleS\n",
       "uZX2qapKUis91/0b2xatzlfV/EiLlKQp1+Xr3LpeU7Vq7j7cYv4UeD3wC+Ag4LHA54FjgLmqujPJ\n",
       "FuDaqnrKCq+vqsrS7ZKk1a0lO3vr6qmq86rq8Ko6Engt8JWqej1wOXBat9tpwGV91SBJWm6c4/j3\n",
       "fLX4M+D4JLcAL+nWJUlj0ltXz8NlV4+kUUmYY18/+Bww3y3PV+1dnglryU6DX1JTEqqKmc2WQfv4\n",
       "JUmTyeCXpMYY/JLUmN4mcGm5lk4wSZpcntwdyKyfYJImUcKFwFuAq4BTq7h34JJGzpO7kvRAR3WP\n",
       "JwIXDFnIkAx+SS25r3vcCZw+ZCFDsqtnIHb1SOOXsAm4BzhkFrt5wAlcE83gl4Yx63979vFLkpYx\n",
       "+AfQjSwg4cruq6ckjY3BPwxHFkgajME/DEcWSBqMJ3cH0MLIAmmStDRr3lE9E2zWRxZIGoajeiRJ\n",
       "yxj8ktQYg1+SGmMf/xi1dIJJ0jA8uStJjfHkriRpGYNfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8k\n",
       "NabX4E9yUJJvJrk+yc1J3t9t35bk9iS7uv+29lmHJGmf3idwJTm4qu5LshH4OvAO4KXAj6vq/Ad5\n",
       "nRO4JGmdRjKBK8nbkxzyUIuoqj03HTkQ2MDCdegBL0ksSUNYS1fPZmBnkkuSbE2yrsBOckCS64Hd\n",
       "wLVVdVP31BlJbkhycRLvOytJY7Kmrp4kBwAnAG8EngdcAlxcVd9f84GSxwFfAs4Bbgbu7p56H7Cl\n",
       "qt68ZP8C3rto03xVza/1eJLUgiRz7Lv4I8B7RnaRtiRHA28CtgJfAY4DvlxVZ6+jwHcDP62qDy3a\n",
       "dgRwRVU9c8m+9vFL0jqNqo//zCTfAj4I/APwjKp6K/AbwKv289rH7+nGSfIo4HhgV5InLNrtd4Bv\n",
       "768OSdJobFzDPocCr6qqHyzeWFX3J3nFfl67BdjedRUdAHyyqq5J8onuG0QBtwKnP4TaJUkPgdfj\n",
       "l6QZ4vX4JUnLGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGrOW\n",
       "i7RpCiXMse8a3XPAfLc8X7V3WVKDvEhbAxKqyltdSi1YS3ba4pc0dn4jHZYt/gbY4tck8/dztLws\n",
       "syRpGYNfkhpj8M+4hAu7xysTNg1dj6ThGfyz76ju8UTggiELkTQZDP7Zd1/3uBNvaq8J4zfSYRj8\n",
       "s+/U7vGEKu4dtBJpOb+RDsDgn3F7wt7Q14TyG+kAHMffAMdJa1J13Tv3AIfYOBmNtWSnwd8Ag1+T\n",
       "zN/P0XIClyRpGYNfkhpjV8+M8iJYmhZ29YyWffySJp7BP1qDBn+Sg4CvAo8EDgT+rqrOTXIo8Fng\n",
       "ScBtwKuratnZfINfml1+I+3P4C3+JAdX1X1JNgJfB94BnAz8R1V9MMk7gUOq6pwVXmvwS9I6DT6q\n",
       "p6r2TM44ENjAwnjdk4Ht3fbtwCl91iBJeqBegz/JAUmuB3YD11bVTcDmqtrd7bIb2NxnDZKkB+r1\n",
       "1otVdT9wdJLHAV9K8uIlz1eSVfuakmxbtDpfVfO9FCpJUyrJHPvOl6ztNeMa1ZPk3cBPgd8H5qrq\n",
       "ziRbWPgm8JQV9rePX5LWadA+/iSPT7KpW34UcDywC7gcOK3b7TTgsr5qkCQt12dXzxZge5IDWPiA\n",
       "+WRVXZNkF3BJkjfTDefssQZJ0hJO4JKkGTL4cE5J0uQx+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPw\n",
       "S1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8k\n",
       "Ncbgl6TGbBy6AM2GhDlgrludA+a75fmqvcuSJkCqaugaVrSWO8VrMiVUFf7spAGsJTvt6pGkxtjV\n",
       "I8muusbY1aORs6tnuvnzm25ryU5b/JpKtlClh84Wv0Zu3C1GW6ij5f/P6ebJXY1dwoXd45UJm4au\n",
       "R9JyvQZ/ksOTXJvkpiTfSfL2bvu2JLcn2dX9t7XPOjRWR3WPJwIXDFmIpJX13cf/c+Csqro+yWOA\n",
       "byXZARRwflWd3/PxNX73dY87gdOHLETSynpt8VfVnVV1fbf8E+C7wBO7p+1DnE2ndo8nVHHvoJVo\n",
       "3eyqa8PY+viTHAE8B/hGt+mMJDckuTiJv2AzYk/YG/pTy666Bowl+LtunkuBM7uW/8eAI4GjgTuA\n",
       "D4+jDs0eW6gjZ1ddA3ofzpnkEcAXgauq6qMrPH8EcEVVPXPJ9gLeu2jTfFXN91epRmWcwwET5oEX\n",
       "dauXVPGacRx3VnUfnvcAh/itbTokmWPfnBaA9+xvOGevwZ8kwHbgR1V11qLtW6rqjm75LOCYqjp1\n",
       "yWsdxz+lxhz8V7LQLbGTGTqvMOQENcfxT7e1ZGffwf9C4GvAjSyM5AE4D3gdC908BdwKnF5Vu5e8\n",
       "1uCfIkMF1bhaqC0FscE/3QYP/ofD4NdazXowzvrxNFrO3JUkLWOLX1Nv1lvE4zieF72bHXb1qAmz\n",
       "GMRDHk/TzeDXzPJkq7Qyg1/qwZiHq14IvAW4Cjh1Voarqj+e3JVGbICZwl5CQSNn8EvrM+4g9hIK\n",
       "GjmDX1qfcQexVzvVyNnHL63DENey8eSu1sOTu1IPHNWjSebJXUnSMga/JDXGrh5pnbyEgiaZffzS\n",
       "iBjEmhYGvyQ1xpO7kqRlDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqM\n",
       "wS9JjTH4JakxBr8kNabX4E9yeJJrk9yU5DtJ3t5tPzTJjiS3JLk6yaY+65Ak7dN3i//nwFlV9XTg\n",
       "OOAPkzwVOAfYUVVHAdd0601JMjd0DX2Z5fcGvr9pN+vvby16Df6qurOqru+WfwJ8F3gicDKwvdtt\n",
       "O3BKn3VMqLmhC+jR3NAF9Gxu6AJ6Njd0AT2bG7qAoY2tjz/JEcBzgG8Cm6tqd/fUbmDzuOqQpNaN\n",
       "JfiTPAb4HHBmVf148XO1cAuwybwNmCTNoN5vvZjkEcAXgauq6qPdtu8Bc1V1Z5ItwLVV9ZQlr/PD\n",
       "QJIegv3denFjnwdPEuBi4OY9od+5HDgN+ED3eNnS13q/XUnqR68t/iQvBL4G3Mi+7pxzgX8GLgF+\n",
       "DbgNeHVV3dtbIZKkvXrv6pEkTZaJnrmb5Iwk3+0mf31g6Hr6kOSPktyf5NChaxmlJH/e/exuSPL5\n",
       "JI8buqZRSLI1yfeS/GuSdw5dzyitNuFyliTZkGRXkiuGrmXUkmxKcmn3d3dzkuNW23digz/Ji1kY\n",
       "7/+sqnoG8KGBSxq5JIcDxwM/GLqWHlwNPL2qng3cwkIX31RLsgH4K2Ar8DTgdd2ExFmx2oTLWXIm\n",
       "cDOzOZLwL4Arq+qpwLNYmDe1ookNfuCtwPur6ucAVXX3wPX04Xzgj4cuog9VtaOq7u9WvwkcNmQ9\n",
       "I3Is8G9VdVv3e/kZ4JUD1zQyq0y4/NVhqxqdJIcBJwEXATM1eKT7Rv1bVfVxgKr6RVX912r7T3Lw\n",
       "Pxn47STfSDKf5HlDFzRKSV4J3F5VNw5dyxj8HnDl0EWMwBOBHy5av73bNnOWTLicFR8Bzgbu39+O\n",
       "U+hI4O4kf5PkX5L8dZKDV9u51+Gc+5NkB/CEFZ56Fwu1HVJVxyU5hoVRQL8+zvoerv28v3OBExbv\n",
       "PpaiRuhB3t95VXVFt8+7gJ9V1afHWlw/ZrF7YJluwuWlLEy4/MnQ9YxCkpcDd1XVrhm9Vs9G4LnA\n",
       "26pqZ5KPsnANtD9ZbefBVNXxqz2X5K3A57v9dnYnQH+5qn40tgIfptXeX5JnsPAJfcPCVAcOA76V\n",
       "5NiqumuMJT4sD/bzA0jyRha+Wr90LAX179+BwxetH85Cq39mdBMuPwf8bVUtm18zxZ4PnJzkJOAg\n",
       "4LFJPlFVbxi4rlG5nYUehJ3d+qU8yMUvJ7mr5zLgJQBJjgIOnKbQfzBV9Z2q2lxVR1bVkSz80J47\n",
       "TaG/P0m2svC1+pVV9T9D1zMi1wFPTnJEkgOB17AwGXEmPMiEy6lXVedV1eHd39trga/MUOhTVXcC\n",
       "P+yyEuBlwE2r7T9oi38/Pg58PMm3gZ8BM/NDWsEsdiH8JXAgsKP7VvNPVfUHw5b08FTVL5K8DfgS\n",
       "sAG4uKpWHTkxhV4A/C5wY5Jd3bZzq+rvB6ypL7P4N3cG8KmuUfJ94E2r7egELklqzCR39UiSemDw\n",
       "S1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfWqMkx3Q3lnlkkkd3Nyt52tB1SevlzF1p\n",
       "HZK8j4WLfD0K+GFVzeSd4TTbDH5pHbqrV14H/BT4zfIPSFPIrh5pfR4PPBp4DAutfmnq2OKX1iHJ\n",
       "5cCnWbgp0JaqOmPgkqR1m+TLMksTJckbgP+tqs8kOQD4xyRzVTU/cGnSutjil6TG2McvSY0x+CWp\n",
       "MQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5Jasz/A0VOMDlt3mZsAAAAAElFTkSuQmCC\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sherpa.load_arrays(\"mydata\", x, y, e)\n",
    "sherpa.plot_data(\"mydata\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can set the model we want to fit to the data using the `set_model` call. There are different ways to instantiate the models: in this case, we just use the string `polynom1d` to refer to a 1D polynomial. The name of the model will be `poly`, and will be accessible as a Python variable. One can use more object oriented patterns to access and instantiate built-in models. Also, new models can be added by the user as Python functions or from tabular data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "sherpa.set_model(\"mydata\", \"polynom1d.poly\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Several Sherpa commands can be used to inspect the model. In this case we just use a simple `print` to get a summary of the model and its components."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "polynom1d.poly\n",
      "   Param        Type          Value          Min          Max      Units\n",
      "   -----        ----          -----          ---          ---      -----\n",
      "   poly.c0      thawed            1 -3.40282e+38  3.40282e+38           \n",
      "   poly.c1      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c2      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c3      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c4      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c5      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c6      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c7      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.c8      frozen            0 -3.40282e+38  3.40282e+38           \n",
      "   poly.offset  frozen            0 -3.40282e+38  3.40282e+38           \n"
     ]
    }
   ],
   "source": [
    "print(poly)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "By default, only the first component (the intercept) is __thawed__, i.e. is free to vary in the fit. This corresponds to a constant function. In order to fit a parabola, we need to __thaw__ the coefficients of the first two orders in the polynomial, as shown below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "sherpa.thaw(poly.c1, poly.c2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We are going to fit the dataset using the default settings. However, Sherpa has a number of optimization algorithms, each configurable by the user, and a number of statistics that can be used to take into account the error and other characteristics of data being fitted."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dataset               = mydata\n",
      "Method                = levmar\n",
      "Statistic             = chi2\n",
      "Initial fit statistic = 11908.3\n",
      "Final fit statistic   = 7.5272 at function evaluation 8\n",
      "Data points           = 11\n",
      "Degrees of freedom    = 8\n",
      "Probability [Q-value] = 0.48096\n",
      "Reduced statistic     = 0.9409\n",
      "Change in statistic   = 11900.8\n",
      "   poly.c0        22.6647     \n",
      "   poly.c1        0.0391891   \n",
      "   poly.c2        1.00155     \n"
     ]
    }
   ],
   "source": [
    "sherpa.fit(\"mydata\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Notice that Sherpa used a Levemberg Marquadt minimization strategy, and the $\\chi^2$ error function. \n",
    "The best fit values of the parameters $c_0=22.9, c_1=-0.1, c_2=1.1$ are close to the ones defined when we generated the dataset, i.e. $23.2, 0, 1$ "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In order to get immediate feedback of the fit results, we can plot the fit and the residuals. Again, Sherpa facilitates the creation of the plots, but users can harvest the power of matplotlib directly if they want to."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEPCAYAAACqZsSmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYXHXZ//H3JyAg8CCJCMkDoVhQKQqogKKyCsTQiwKC\n",
       "SBVQaQ8qTR6dPfLYUIr9EglKE4006RCRFbBggFADEhEU+NGUrqBA7t8f5yRsll22ZOZ8z5zzeV3X\n",
       "Xrszc2bPfTLZe77zbbciAjMza5ZxqQMwM7PyOfmbmTWQk7+ZWQM5+ZuZNZCTv5lZAzn5m5k10KJl\n",
       "nETSvcBTwIvA8xGxvqQJwM+BVYB7gZ0i4oky4jEza7qyWv4B9ETEuhGxfnHfkcCMiFgduLK4bWZm\n",
       "JSiz20cDbm8DnFr8fCqwXYmxmJk1Wpkt/19Jul7SvsV9K0TEw8XPDwMrlBSLmVnjldLnD2wUEQ9K\n",
       "eh0wQ9Kd/R+MiJDkfSbMzEpSSvKPiAeL749KOg9YH3hY0sSIeEjSJOCRgc/zG4KZ2dhExMCu9gV0\n",
       "PPlLWhJYJCKelrQUMAXIgAuAPYCvF9/PH+z5w11AN5PUGxG9qePolDpfX52vDXx97TkHPUBPcbMH\n",
       "6Ct+7ouY/3OHzj18w7mMlv8KwHmS5p3vzIi4QtL1wHRJ+1BM9SwhFjOzUhQJvg9AIiLmvxFUQseT\n",
       "f0TcA6wzyP2PAZt2+vxmZvZyZQ34VlbKj2b9zlVXfakD6KC+1AF0WF/qADqsr6wTSZxUfL8E2DWC\n",
       "SixmVZWLuUiKMvv8i49mtR1jMLPySfQBGxc3p0ewc+fPOXzu9N4+ZmadNPm3i7DEEwAzgf0TRzOf\n",
       "k39Bmjv/o5nEsqnjMbPup0yrsPf73siKfwSYUpUuH3Dyf8n2e27BO34IsDnww8TRmFmX06SbtuDp\n",
       "iX/kxr1v5u4pvwH+R6K3GGdMzn3+88418aY+Pj51Yy7+7p+44yMbVukd2sy6izItClxIPo3909Eq\n",
       "N9GOJHc2frbPfA+vsx3Tz36cvd7/WhQrgbeXNrMxO4G8Z+XgshP/SDW+20eiR6IX+B/+9t7fcO2R\n",
       "V/Ps+Gu1+sU7pI7NzLqPMh0EfBDYKVrxfOp4huJun8HOm+mLwFZAT7TiX2Wf38y6kzJtAUwD3hOt\n",
       "uCdZHFWa6ilpEUmzJF1Y3O6VdH9x3yxJU8uKZQSOAe4ETlemxn86MrPhKdPbgJ8AH06Z+EeqzMR2\n",
       "CDCbfG9/iu/HF9W91o2Iy0qM5RUVfXT7AssBX0scjplVnDJNJB/gPSRa8bvU8YxEKclf0krAFsDJ\n",
       "vFTRS/1+rpxoxb+BHYDtlGm/1PGYWTUp06uBXwLTohVnpY5npMpq+Z8AHAbM7XdfAAdJulnSNEmV\n",
       "W1gVrfgHsCXwJWWakjoeM6uWolv4VGAOeXdx1+h48pe0FfBIRMxiwZb+D4DVyHf8fBA4rtOxjEW0\n",
       "Yg6wI3CGMq2VOh4zq5RjgBWBT1R1SudQOj7bR9JXgI8DLwBLAMsA50TE7v2OWRW4MCLWHvDcIC/8\n",
       "Mk9fRPR1NOAhKNNu5C/0htGaX3vYzBpKmfYAvkieEx5NGovUAwusHG4NN9un1KmekjYGPhcRW0ua\n",
       "NK+8o6RDgXdFxK4Djk8y1XMoytQiH7v4gKeAmjWXMr0f+AX5dPA7UsczUKWmehbES7N9jpV0i6Sb\n",
       "ybc7PbTkWMbiS8BdwGmeAmrWTMr0JmA68LEqJv6R8iKvUVKmxYEZwO+jFUekjsfMyqNME4DfA8dF\n",
       "K05KHc9Qqtjy73rFFNDtgR2Uad/U8ZhZOZRpMeAc4KIqJ/6RcvIfg2IK6BbAMcq0Wep4zKyzlEnk\n",
       "MxSfAg5PHE5bOPmPUTEFdCfgTGVaM3U8ZtZRhwHrkffzv5g6mHZw8l8I0Yqrgc8AFynTCqnjMbP2\n",
       "U6YdgIOBraMVz6SOp12c/BdStOIM8hV+FyjTkqnjMbP2UaZ3klf22zZacX/qeNrJyb89MvLl3Z4C\n",
       "alYTyjQZOB/YN1pxQ+p42s2Jqg2KZd37ACsAX0kcjpktJGVamnyXzhOjFeenjqcTnPzbpN8U0A8r\n",
       "0ydSx2NmY6NMiwBnATOp6J5j7eDk30bRir+T7wL6ZWXaNHU8ZjYm3wSWJEHh9TKVtZ//wCpeEyTN\n",
       "kHSXpCuquJ3zWEUr7iLfBfSnyrRG6njMbOSU6ZPA5sBHqlx/tx3KavkPrOJ1JDAjIlYHrixu10Yx\n",
       "BfRzwMWeAmrWHYqaHb3AltGKxxOH03Fl7Oc/WBWvbcinR1J8367TcZQtWnEacDrwy6LSj5lVVLFQ\n",
       "8wxgx2jF3anjKUMZLf/BqnitEDF/T/yHyWfJ1FELuBtPATWrLGVannxmz2eiFdekjqcsHU1Ir1DF\n",
       "a77ItxWt5aBKvymgk4AvJw7HzAZQpiXI5/KfWSzYbIxFO/z73wNsI2kLiipekk4HHpY0MSIekjQJ\n",
       "eGSoXyCpt9/NZJW8xipa8ZwybQf8QZnmRCtOSR2Tmc3frO0U4G/kn9K71iCVvIZ/Tln7+Q+o4nUs\n",
       "8I+I+LqkI4FlI+Jlg75V3M9/rJTpzcDVwK7RiitTx2PWdMrUC3wI+GC04tnE4bRVFffzn/dO8zVg\n",
       "M0l3AR8sbtdatOJP5LuAnuUpoGblk+iR6JXo1Qbfmc2z4w/he7ddS29skDq2FFzJq2TKtDv5dLIN\n",
       "oxVDdneZWWco0/v45+uuZqlH145W3JY6nk6oYsvfeuNvzN7hPh5f5QYt/tTV81siGl1/nZmNnjJN\n",
       "5d//dQXnnAm9caxEbRaYjpZb/gkUA00/4a/v3Z1Vrp3QhAUlZqkp0y7AiZxxyQP8efN1i7unR7Bz\n",
       "yrg6wS3/iiqmgO7Fg+8AuFqZVkwcklmtKdOBwLHAJvx584eKu2cC+6eLKi23/BOROAliX6Z89k9s\n",
       "+K0lGTd3SrTiztRxmdVJ8Sm7F9gFmBKtuLfo6nkcGB/BEynj6xS3/KttdRBccfybuepLDwJ9yrR+\n",
       "6qDM6qLYmvl7wFbAe6MV9wLMS/h1Tfwj5ZZ/IhKXkO8eOBOYQq/eR77g5OPRisuSBmfW5ZRpcfK9\n",
       "tZYDtotWPFVMqugpDukB+oqf+yLm/1wLI8mdTv6JDPbRU5k2As4FPtu0peZm7aJM/0X+d/QU8LFo\n",
       "xXOJQyqdu30qbLCPntGK35IvevuKMn0mVWxm3UqZliPfJv4eYKcmJv6Rcsu/ZCP56KlMKwOXk+80\n",
       "eESdqwmZtUvxd3MFcA7wv03+u3G3TxdTptcCFwF3AZ+oe1Uhs4WhTG8lbzAdH604MXU8qVWi20fS\n",
       "EpKuk3STpNmSvlrc3yvp/qK84yxJUzsdSzeJVvwD2JR8wOp8ZVoqcUhmlaRMGwBXAUc78Y9cKS1/\n",
       "SUtGxL8kLQpcS17icBPg6Yg4/hWe19iW/zzK9CryKmirA1sVbwpmBijTh8grcO0Zrbg4dTxV0baW\n",
       "v6SDJY0fayAR8a/ix8WARchnucAQBV7sJUV3z57ANcC1Rb+mWeMp00eB08incjrxj9JIu31WAGZK\n",
       "mi5pqqRRJW1J4yTdRF6y8aqIuL146CBJN0uaJqmxGywNJ1oR0YrDgR+RvwGsmToms5SU6QDgm8Cm\n",
       "xSw5G6URd/tIGgdMIW+FvhOYDkyLGHmxY0mvIR+UORKYDTxaPHQMMCki9hlwfABZv7u6rpJXuynT\n",
       "bsBxwA7+T29NU2zX0AJ2pdiuIW1E1TBIJa9WW2f7SFoH2AuYCvwa2BD4VUQcNorf8QXg2Yj4Zr/7\n",
       "VgUujIi1Bxzb+D7/wSjTVPLVi3tHKy5MHY9ZGYrtGr4NvBuY6noYQ2tnn/8hkm4g3xXvt8BaEfEp\n",
       "4B3ADsM8d7l5XTqSXg1sBsySNLHfYdsDt44kFoNi+4etgB8p016p4zHrtGK7hp8CawA9TvwLb6QF\n",
       "3CcAO0TEX/vfGRFzJW09zHMnAacW3UbjgNMj4kpJpxWfJIJ8NV5jt1Ydi2jFdcq0MXC5Mk0Evtbk\n",
       "RS1WX8q0NHAe8DSwuVfttocXeXW5ohbAZeRL2j8TrZibOCSztim2a7gEuAX4ZLTihcQhdYVKLPKy\n",
       "zopWPAC8n7wL7gxlWixxSGZtUUxrvgb4FbCvE397ueVfE8r0auAsYCnymUBPJw7JbMz6bddwQrTi\n",
       "hNTxdBu3/BskWvEs8BHgXuDXyrR82ojMxqbfdg3/68TfOW7510wxD/pLwM7Ah+iNVWhIAQvrfso0\n",
       "BTgT2CtacVHqeLqVd/VssKJg9ZHAltGKmwEkIsJbalg1vGx7854vPsq7j9+cG/Y/Ki4/7jvpIut+\n",
       "Tv4Np0w7Ad8FdqQ3PgbsC1wK7Nr0+qVWLVr/+8GWBzwAbBGtuCV1PN3Off4NF62YDuwCnM26095T\n",
       "3L058MN0UZm9RJmW0i7b3caGJ8D3b51Db/wtdUxN4eRfc9GKK4EPsdnhb2Szw2Gxp2/AC+qsApTp\n",
       "g8At6MXlOPkP8MhaPbhhUhp3+zSEtt7vjSzx1BzW+MU9jJu7T7TiqtQxWTMp0zLkW8VsCXyS3jiA\n",
       "/BPpTGCKuyQXXiX6/CUtAfwGWJx8P/9fRsRRkiYAPwdWIZ+euFNEPDHguU7+bSQR9Gpb4PvkqyYP\n",
       "j1b4D81Ko0zzuh0vAw6LVjwpsSx5jY/xTvztUYnkXwQyWCWvbYC/R8Sxko4AxkfEkQOe5+S/kAYt\n",
       "GL/MfYuzxyZr8to56wEHRCt+mSg8awhlmgAcD2xMvlr3Vws87plobVWZ5D//ZNKS5J8C9gTOATaO\n",
       "iIeLHT77IuItA4538u+gYmO4k4EbgYOjFQ8nDslqSJm2A74HnAscFa14BoZomOS8BmUhVSb5Fzt6\n",
       "3gi8AfhBRBwu6fGIGF88LuCxebf7Pc/Jv8OKbSF6yd+QDwNO9+6g1g7K9DrgO8B6wD7RimsSh9QY\n",
       "lUn+80/2UiWvo4Bz+yd7SY9FxIQBx7uSV0mU6R3ANPJSm/u7QpKNVbHKfGfgRPKiQ61oza/jbR3Q\n",
       "8Upe7TCvkhfwCaAnIh6SNIm8tq+7fRJSpleRj8d8lnyLiO9FK15MG5V1E2WaRD6hYHXySnPXJQ6p\n",
       "kSqxyGuoSl7ABcAexWF7AOd3OhZ7ZdGK56MVXwU2AnYkLxa/RuKwrAsok5RpD+Bm4DZgPSf+aitj\n",
       "qufawKnkbzTzKnl9o5jqOR1YGU/1rBxlGgd8kvwTwLfJK4X9J21UVkXFvvs/JK/at1e0YlbikBqv\n",
       "cn3+o+Xkn54yTQZ+QL4eY59oxR8Th2QVUTQQ9gX+D/gW8PVoxfNpozJw8rc2KQbwdiGfp30m8MVo\n",
       "xT/TRmUpKdPryacJL03e2r89cUjWj5O/tVUxde9EYENgv2LfIGsQZVoEOBD4AvA14ESXV6weJ3/r\n",
       "CGXakrwraAbwuWjF44lDshIo01vIpwPPJe8CvCtxSDaESsz2sfqJVlwMrAU8B9ymTDskDsk6SJkW\n",
       "VaYjyLdmOQvY2Im/+7nlbwtFmd5H3vd7K3BgtOKhxCFZGynT2sApwJPke/LckzgkGwF3+1gplGkJ\n",
       "4IvkC/eOAH7iLSK6mzItRl4G9CDyFfnT/Jp2Dyd/K5UyrUPeJ/wY+RYRf0kcko1SMaA7FfgKcD/5\n",
       "63h/2qhstJz8rXRa4olN2PyQI1nzFxtx95Sn+NPWlzBrr7/BOO/UWGHFQq29gX2A/wecAPzcrf3u\n",
       "5ORvySjT8lz+zYf50OfuABYl7zc+NVrxYOLQrFDs5bQlsB+wAflg7o+iFTcnDcwWWiWSv6TJwGnA\n",
       "8kAAJ0XEtyX1kvcRP1ocelREXDbguU7+XayoHDaOfF3A3sBHgGvIu4Yu8WrQNIoFWp8g38b7L8BJ\n",
       "wNn0xvp4f/1aqErynwhMjIibJC0N3ABsB+wEPB0Rx7/Cc538u5TESeRL/y8Fdo3gCWVamnzDuH2A\n",
       "N5I3Ck6JVtw5xnP04GQ1IsUA7nbkr8k65FstnxytmD3o8a6s1dUqkfxfdkLpfOC75DtHPhMRx73C\n",
       "sU7+XUqij7xkH8D0CHZe4PFMbyb/NLA7eetzGjB9XpWnMZzPyWoQyrQ6ecLfHZhN3so/L1rx3JDP\n",
       "GeSNu4xYrX0ql/wlrUpexnFN8j3j9yKfP3w98Fnv6lkfEpcAmwMzgSlDJZCi33kL8jeC95OX+psG\n",
       "/H6kg41OVgsqpt5+mPzf5C3ku+qeHK2YM6LnD/PGbdVXqeRfdPn0Af8XEedLWp6X+vuPASZFxD4D\n",
       "nuNKXl1KYlngcWD8SJNxUQhkd/I3grnkg8SnDVdb2Mkqp0xrkif83ci7V38EXDDarbhH+sZt1VHZ\n",
       "Sl6SXgVcBFwaEScO8viqwIURsfaA+93y72Jj7YopdhHdiHxsYHvgKvJPA5cNtolYXZPVSMY0lGlJ\n",
       "8vGz/ci33f4x+YKsMa/EHcsbt1VLJVr+RXH2U4F/RMSh/e6fFJFP+5N0KPCuiNh1wHOd/LtMuwdh\n",
       "lWkZ8uS2D3lyO5V8kHh+F0ZZySrlAPPAN9JiQd1+wEeB35G38i9u1w6bHkPpblVJ/u8FrgZuIZ/q\n",
       "CfB58v3h1ynuuwfYP2LBj/dO/tZf0a2xN/Bx4A7ybqGzoxX/LDtZpTgfvVqGPNnvB6xA/mnolGjF\n",
       "fW06Rw+ePVULlUj+C8PJ3wZTTFvciqdXOIwlnlqH+ze4nd8dtgiPrHU5T678HB1OVmUMMBdVsiYC\n",
       "b+Cq1jEs88DGrPWzF1jk+ctZ9N/fA66IVrzY7vNaPTj5W+0p00rkg8Q7kq8deJG8JvSgX9GKhU7U\n",
       "7RpgVqbFgdWA1wNvKL7m/bwa8DRwN3dt+Qbu3Xh5btkNnpnU2AFtGzknf2uUYqB4ArDqEF+rAc/z\n",
       "ym8OTw57npFPYxUwnpcn9nk/Lw/cR77O4e7ia97Pf5m35qGuA9rWOU7+Zv0sxJvDPbz05vBU/wFm\n",
       "evU0MJnBW+9vAMRLib1/cr8buH8kA7SefWOj5eRvNgoD3hxWY/A3iP8A9/KXTdbl9VfOAVYmX68y\n",
       "WHL/C/DYWHfG9ACsjZWTv1kbadwLPYy/Zwsm3bgssAHPTriGB971BP9e9ldOxlYlTv5mZg3kAu5m\n",
       "ZjYoJ38zswZy8jcza6COJ39JkyVdJel2SbdJOri4f4KkGZLuknSFpGU7HUvVFDvx1Vadr6/O1wa+\n",
       "viYoo+X/PHBoRKxJXs7vAElvBY4EZkTE6sCVxe2m6UkdQIf1pA6gg3pSB9BhPakD6LCe1AGk1vHk\n",
       "HxEPRcRNxc/PkG/ItSKwDfkOjRTft+t0LGZmliu1z7/Yt39d4DpghX67eD5MvkuhmZmVoOxKXr8B\n",
       "jikqeT0eEeP7Pf5YREwY8JzqLkIwM6uw4eb5L1pGEEUlr3OA0yPi/OLuhyVNjIiHJE0CHhn4PC/w\n",
       "MjPrjI4n/6KS1zRg9oASjhcDf5T0GDCJfB8UMzMrQRl9/huRF5T+gKRZxddU4MvAn4Elyat8LVpU\n",
       "/TIzsw6rxN4+kpYkHw/YIyJmp47HzKzukq7wlTRO0k3ks32ucuI3MytHKQO+Q4mIucA6kl4DXC6p\n",
       "JyL65j3u2T5mZmNTidk+w4mIJyVdDLwTFtwXvc4zfiT1RkRv6jg6pc7XV+drA19ftxtJwzlZt4+k\n",
       "5ebt5yPp1cBmwKxU8ZiZNUnKlv8k4FRJ48jfhE6PiCsTxmNm1hjJkn9E3Aqsl+r8FdGXOoAO60sd\n",
       "QAf1pQ6gw/pSB9BhfakDSK0SUz2H4jKOZmaj5zKOZmY2KCd/M7MGcvI3M2sgJ38zswZy8jczayAn\n",
       "f6sFid7UMZh1E0/1tFqQiAj8f8UMT/U0M7MhpNzbZ7KkqyTdLuk2SQenisXMrGmSdftImghMjIib\n",
       "iuLuNwDbRcQd/Y5xt4+NiLt9zF5S6W6fiHgoIm4qfn4GuAP471TxmJk1SSX6/CWtCqwLXJc2EjOz\n",
       "ZkhezKXo8jkbOKT4BDDw8d5+N/v6V/oyS0WiN8LTS60aJPUAPaN6TsqpnpJeBVwEXBoRJw7yuPv8\n",
       "bUTK7vP3GINVWaX7/CUJmAbMHizxm5lZ56Ts898I2A34gKRZxdfUhPGYmTWGV/haLbjbx+wlle72\n",
       "MTNriiruPeWWv9WCW/5WZeX//3TL38zMBuHkb2bWQE7+ZmYN5ORvZtZATv7W9SROKr5fIrFs6njM\n",
       "uoGTf4NUcbpZm6xefN8c+GHKQMy6hZN/PzVOjvO0UgfQIf8qvs8E9k8ZiFm3SL2x2ynAlsAjEbH2\n",
       "II+XOs+/7nO363p9RVfP48D4CJ4o6Zy1/Le0zvA8/5f7MeD9fGyhzEv4ZSV+szpImvwj4hryFpuZ\n",
       "mZUodcvfrOukmF3UgPEoK1nySl7DcSUvq6CBs4t2LuGcLfAbgA2u6yp5wfz6vRd6wLfz6nx9ZV6b\n",
       "xCXkiX8mMKWMsYY6v3ZN4AFfs3rYtfheSuI364SkyV/SWcDvgNUl3Sdpr5TxmI2EZxe1n8c0ype8\n",
       "2+eVuNunvep8fXXfz7/Orx3U+/qKCQL7ApcCu5bTTehuHzOz1Cq5/YiTv5lZZ1Vy+xEnfzNrnJLH\n",
       "GCo5QcB9/gucr779jlDv66t7H3ydXzuo/7+np3qamVklOPmbmTWQk39DuNqVmfXn5F9oQHKs5HSz\n",
       "NspSB2DWTZz8X1L35FjJ6WbtEuEVomajkXp7h6mS7pQ0R9IRKWMhQXL0dDMzSyXZVE9JiwB/AjYF\n",
       "HiBPurtExB39jiltqmcTSgHWfbpgmfzatVfd/z091XNB6wN/joh7I+J54GfAtqmC8WZdVlUNGI+y\n",
       "BIZN/pIOljS+A+deEbiv3+37i/usJrxTY9vUfTzKEhhJJa8VgJmSbgROAS6P9vQVjeh3SCqxXyp8\n",
       "vvafr1Xe+cpU5r/lxcAWwB+BKTtJT+5UznnL1Ii/hUptpzBs8o+IoyV9AZgC7Al8V9J0YFpE3L0Q\n",
       "534AmNzv9mTy1v/A85fcL+fzdev5yiTRGxG9JZ2rGI9af3zEE7Xslqz7/83yzzf8G82I+vwjYi7w\n",
       "EPAw8CIwHjhb0jcWIr7rgTdJWlXSYuR1UC9YiN9nVpoyp5bWfTzKYxppjKTP/xBJNwDHAr8F1oqI\n",
       "TwHvAHYY64kj4gXgQOByYDbw8/4zfcysMTymkcBI+vwnADtExF/73xkRcyVtvTAnj4hLyavbmFlz\n",
       "1XoBYlV5S+cFzlf3ub/1Pl+d1fnfsiFrbHrL7Cqs+jx/M7PSxzRSjDFUcfsRJ38zaxqPMeDkb2bN\n",
       "4zEG3Oc/4Hz17hOv+/nqrO7/lmVeX4oxhrK5z9/MxqTOW3PUfd3ESLnlv8D56t0yrvv56qzur13d\n",
       "z1c2t/xtIFe7MjPAyT8ZTzczs5SSJH9JO0q6XdKLktZLEcMQymwZe7qZmSWTquV/K7A9cHWi8w+q\n",
       "5Jaxp5uZWTIj2dun7SLiTgCptuMtI7Er+XQz19Q1s9K5zz+Ruk838za9befBemurjrX8Jc0AJg7y\n",
       "0Ocj4sJOndcqY+CYxs4JY+l6Hqy3dutY8o+IzdrxeyT19rvZFxF97fi91nEe0zAriaQeoGdUz0m5\n",
       "yEvSVcDnIuKGIR6vbRlAqPdCkyYsoa+zui+6qvPfHlR4kZek7SXdB2wIXCzJBV1qpu5jGmbdLtVs\n",
       "n/OA81Kc28wqyQPaJfPePgnV/6Nnva+vzur+2tX/+ira7WNmZmk5+ZtZEzW+m8ndPgnV/6Nnva+v\n",
       "zvzadTd3+5jZqHl1djM4+ZvZQN5xtgGc/M1sIK/ObgD3+SdU937Vul9fXXl1dvdzn7+ZjZpXZzdD\n",
       "qu0dviHpDkk3SzpX0mtSxGFm1lSpWv5XAGtGxNuBu4CjEsWRWuPnGptZGkmSf0TMiIi5xc3rgJVS\n",
       "xJGa92g3s1Sq0Oe/N3BJ6iDMzJokaSUvSUcD/4mIn3YqDkvK3VpmFZVsqqekPYF9gU0i4rkhjgkW\n",
       "TCCu5GVWAk/T7S6DVPJqDTfVM0nylzQVOA7YOCL+/grH1Xqev1lVOfl3t5HkzlTJfw6wGPBYcdfv\n",
       "I+LTgxzn5G+WgJN/d6ts8h8pJ3+zNJz8u5tX+JqZ2aCc/M3MGsjJ38ysgZz8zcwayMnfzKyBnPzN\n",
       "bDBenV1znuppZlYznuppZmaDcvI3M2ugVJW8jimqeN0k6UpJk1PEYWbWVKla/sdGxNsjYh3gfKCV\n",
       "KI6kip34aqvO11fnawNfXxOkquT1dL+bSwND7uxZcz2pA+iwntQBdFBP6gA6rCd1AB3WkzqA1DpW\n",
       "zGU4kr4MfBz4F7BhqjjMzJqoYy1/STMk3TrI19YAEXF0RKwM/AQ4oVNxmJnZyyWf5y9pZeCSiFhr\n",
       "kMequwjBzKzChpvnn6TbR9KbImJOcXNbYNZgx3mBl5lZZ6Sq5HU28GbgReBu4FMR8UjpgZiZNVTy\n",
       "bh8zMytf5Vf4SjpI0h2SbpP09dTxdIKkz0qaK2lC6ljaSdI3itfuZknnSnpN6pjaQdJUSXdKmiPp\n",
       "iNTxtJOkyZKuknR78Td3cOqY2k3SIpJmSbowdSztJmlZSWcXf3ezJQ05k7LSyV/SB4BtgLcVA8Lf\n",
       "TBxS2xWrmzcD/po6lg64AlgzIt4O3AUclTiehSZpEeC7wFRgDWAXSW9NG1VbPQ8cGhFrkk/BPqBm\n",
       "1wdwCDAbqGO3x7fIJ9C8FXgbcMdQB1Y6+QOfAr4aEc8DRMSjiePphOOBw1MH0QkRMSMi5hY3rwNW\n",
       "ShlPm6wP/Dki7i3+X/6MfNJCLUTEQxFxU/HzM+TJ47/TRtU+klYCtgBOhnoVqC8+Wb8vIk4BiIgX\n",
       "IuLJoY6vevJ/E/B+SX+Q1CfpnakDaidJ2wL3R8QtqWMpwd7AJamDaIMVgfv63b6/uK92JK0KrEv+\n",
       "xl0XJwCHAXOHO7ALrQY8KunHkm6U9CNJSw51cLIVvvNImgFMHOSho8njGx8RG0p6FzAdeH2Z8S2s\n",
       "Ya7vKGBK/8NLCaqNXuH6Ph8RFxbHHA38JyJ+WmpwnVHHroKXkbQ0cDZwSPEJoOtJ2gp4JCJm1XRv\n",
       "n0WB9YADI2KmpBOBI4EvDnVwUhGx2VCPSfoUcG5x3MxiUPS1EfGP0gJcSENdn6S1yN+pb5YEeZfI\n",
       "DZLW76Zpr6/0+gFI2pP8Y/YmpQTUeQ8A/XehnUze+q8NSa8CzgHOiIjzU8fTRu8BtpG0BbAEsIyk\n",
       "0yJi98RwiN4PAAABu0lEQVRxtcv95D0JM4vbZ5Mn/0FVvdvnfOCDAJJWBxbrpsT/SiLitohYISJW\n",
       "i4jVyF+49bop8Q9H0lTyj9jbRsRzqeNpk+uBN0laVdJiwM7ABYljahvlLZFpwOyIODF1PO0UEZ+P\n",
       "iMnF39tHgV/XKPETEQ8B9xW5EmBT4Pahjk/e8h/GKcApkm4F/gPU5oUaRB27E74DLAbMKD7d/D4i\n",
       "Pp02pIUTES9IOhC4HFgEmBYRQ86o6EIbAbsBt0iat/L+qIi4LGFMnVLHv7mDgDOLhsndwF5DHehF\n",
       "XmZmDVT1bh8zM+sAJ38zswZy8jczayAnfzOzBnLyNzNrICd/M7MGcvI3M2sgJ38zswZy8jcbBUnv\n",
       "KorTLC5pqaLgyRqp4zIbLa/wNRslSceQbwz2auC+iKhlhTmrNyd/s1Eqdr28HngWeHf4j8i6kLt9\n",
       "zEZvOWApYGny1r9Z13HL32yUJF0A/JS8sNCkiDgocUhmo1b1LZ3NKkXS7sC/I+JnksYBv5PUExF9\n",
       "iUMzGxW3/M3MGsh9/mZmDeTkb2bWQE7+ZmYN5ORvZtZATv5mZg3k5G9m1kBO/mZmDeTkb2bWQP8f\n",
       "8ihYr2EIClMAAAAASUVORK5CYII=\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sherpa.plot_fit_resid(\"mydata\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can now compute the confidence intervals for the free parameters in the fit:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dataset               = mydata\n",
      "Confidence Method     = confidence\n",
      "Iterative Fit Method  = None\n",
      "Fitting Method        = levmar\n",
      "Statistic             = chi2gehrels\n",
      "confidence 1-sigma (68.2689%) bounds:\n",
      "   Param            Best-Fit  Lower Bound  Upper Bound\n",
      "   -----            --------  -----------  -----------\n",
      "   poly.c0           22.6647    -0.455477     0.455477\n",
      "   poly.c1         0.0391891   -0.0953463    0.0953463\n",
      "   poly.c2           1.00155   -0.0341394    0.0341394\n",
      "poly.c0 lower bound:\t-0.455477\n",
      "poly.c1 lower bound:\t-0.0953463\n",
      "poly.c2 lower bound:\t-0.0341394\n",
      "poly.c0 upper bound:\t0.455477\n",
      "poly.c1 upper bound:\t0.0953463\n",
      "poly.c2 upper bound:\t0.0341394\n"
     ]
    }
   ],
   "source": [
    "sherpa.conf(\"mydata\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Sherpa allows to inspect the parameter space. In the cell below we ask sherpa to show us the projection of the confidence regions for the `c0` and `c1` parameters. The contours are configurable by the user: by default they show the confidence curves at 1, 2, and 3 $\\sigma$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEZCAYAAACTsIJzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYHGW5/vHvnbCIBGTxSDAkRBABUWQA48TjMm4YRMCj\n",
       "LIdFCSKLPxcEVEDUSutRgXME3A8CmiCygwgHZREYFaXZHEAIYZVdQDYVFIXw/P54q01Np7unZ6a7\n",
       "36rq53Ndc2W6u7r7mZpM3f3UW2+VzAznnHOuZkrsApxzzuWLB4NzzrlRPBicc86N4sHgnHNuFA8G\n",
       "55xzo3gwOOecG8WDwfUFSbMk/VWSYtfSiqTvSfp8UV7XlZN8HoPLG0n3AC8DlgLPAJcCHzOzv8Ss\n",
       "qxlJC4AjgGeB54HFwCFmVo1Uz3xgHzN7c4z3d8XnHYPLIwPea2arAa8DXgvk+dOuAael9f4bcCVw\n",
       "bqMFJfnfnMs9/0/qcs3MHgEuATYDkDQo6beSnpR0g6S31paV9ApJv5L0F0mXSvqOpB+lj82W9EJt\n",
       "wyzp5ZLOl/S4pDskfSTzOgsknSlpUfpaN0vaqkWZSr8ws+eBk4HpktaWtDDdjfMzSU8DQ5I2lTSc\n",
       "/gw3S9o+894LJX05c/u96c/5pKTfSHpt5rGZks6V9KikxyR9S9ImwP8Cc9NdZ080ed1905/7cUk/\n",
       "lbRu5rEXJO0v6fb0fb89rl+aKzwPBpdXApC0HjAPuFrSDOD/gC+Z2ZrAp4FzJK2dPudUoAqsBSwA\n",
       "9iR8mm/kdOA+YF1gJ+Crkt6WeXx74DTgJcD5QFsbR0krA/OB+8zs8fTu3YAvm9k04FrgAuAiQnfx\n",
       "CeDHkl6VLmu1miUNACcB+6Y/0/HA+ZJWlDQ1XRd/ANYHZhC6liXA/sBVZraama3V4HXfDnwV2Dn9\n",
       "+e9N10fWdsDWwObALpLe3c7P78rBg8HlkYDzJP2FsPG+C/gKYUP/MzO7CMDMfgFcB2wnaRZhQ/ZF\n",
       "M3vezH5D2KAvN9gsaSbwRuBQM/unmd0InAh8KLPYr83sIguDcKcQdmm1soukJ9N6B4D/yDx2npld\n",
       "lX6/BbCqmR2Z1nkFYQO/W4PX3A843syuteBk4B/AXGAOYaP+GTP7u5n9w8x+m1l/rewBnGRmN5jZ\n",
       "P4HDCR3GrMwyR5rZX8zsfuCKtG7XJzwYXB4ZsKOZrQ4MAW8HtiJ8Mt453b3xZLoh/ndgOvBy4Akz\n",
       "ezbzOvc3ef3ass9k7ruP8Km75pHM938DXiRpiqQ90l00f5V0YWaZM8xsTTNbx8zeaWYjmZ/lgbr3\n",
       "rq/r3vT+eusDh9T9vOsRAmEmcK+ZvdDkZ2yl1iWEAsN6eJzRP//Dme//BkybwPu4glohdgHOtWJm\n",
       "v5L0LeAo4GLgR2a2X/1yktYH1pK0ipn9Pb17Fo13JT2ULjvNzJ7OLPtAg2Xr6/kx8OP6u2n9KT1b\n",
       "w0PATEmyZYcErg8safC8+4CvmNlX6x+QNBeYJWmqmS1t8X6NPATMzrzWqsDawINjPM/1Ce8YXBEc\n",
       "R9h1ciWwvaRtJE2V9CJJQ5JmmNm9hN1KC9J98HOB99JgI5nuHvkt8DVJK0vaHPgwYZfRRLQKhfrH\n",
       "qoRP4J9N6xxK6zy9wXNOAA6QNEfBqpK2kzQNuBr4I3CkpBen6+KN6fMeAdaTtGLda9Ze9zRgb0mv\n",
       "S8dEvgpUzey+Cfx8roQ8GFzumdljwCLgYGAH4HPAo4RP1Iew7P/xHoT9748DXwbOAP6ZfanM97sR\n",
       "PjU/RDi09ItmdnlmufpAafUpvNHyDR8zs+cIA9vbAn8iDGp/0MxuX+6JZtcTBp6/DTwB3EE6DpLu\n",
       "QtoeeCVhPdwP7JI+9TLgFuBhSY/W12FmlwFfAM5Jf/5XAP/Z4mdt9fO5EvIJbq60JJ0BLDazSuxa\n",
       "2iVpEXCHmf1X7Fpc//KOwZWGpK0lbZgOEm9L6C7Oi11XuyStAGxMOATVuWh88NmVyXTCbqG1CbtW\n",
       "DkgPRS2KhwnzHM6JXYjrb74ryTnn3Ci+K8k559wopdiVJMnbHuecGycza3wosplF+yKcA2cJ4TC8\n",
       "Q1ss93rC6Yzf3+Rxm+D7L4j5809ivXndXrfXnbOvotXdarsZbVdSehKwbxPC4dXAbpI2bbLcUYST\n",
       "jvlEG+ec67KYYwxzgDvN7B4Lk35OB3ZssNwngLMJk4Gcc851WcxgmMHok4k9wOiTeJGeZnlH4Hvp\n",
       "XZ0eSxju8Ov1ynDsAiZoOHYBEzQcu4AJGo5dwAQNxy5ggoZjF9ApMYOhnY38ccBhFnaIZc/10pkC\n",
       "zIY7+Xq94nX3ltfdW153fDGPSnqQcOrgmpksf3bLrYDTFa7f/lJgW0nPmdn59S+mcN3dmuEy/ZKc\n",
       "c26y0hM2DrW1bDo63XPp9P/bgHcQTuR1DbCbmd3aZPkfAheY2XLX0k3PYOwD084516ZW281oHYOZ\n",
       "PS/p44Rz7E8lXFHqVkn7p48fH6s255zrZ6U4JYZ3DM45Nz6ttpt+SgznnHOjeDA455wbxYPBOefc\n",
       "KB4MzjnnRvFgcM45N4oHg3POuVFKcT0G5yaqIk0hzLrfLP16NbAB+fvQtBS4G1gM3JJ+3Z+U4Xhz\n",
       "lzs+j8H1hSYBsBmwKfAXlm1sFwN3Eq7/kScrAhsyuv7VgFtZVrcHhmtbq+2mB4MrlXEGwC3A4sTs\n",
       "qTjVTk5FWpNlP1/t380IgfGvnw8PDNeAB4MrnX4KgPHywHDt8GBwheUB0DkeGC7Lg8HlngdAPB4Y\n",
       "/cmDweVORXoJsD3htOseADnUZmBcBlyQmP0lVp1uYjwYXC5UpNWBHYBdCBcMGQYuBH6PB0BhZALj\n",
       "tcB7gTcDVwBnAecnZn+NWJ5rkweDi6oibQgcCOwJXMmyDcifoxbmOqIirUG4NvvOwBuBRcA3E7M/\n",
       "RC3MteTB4HquEq7H+ibgIOAtwAnAtxOzB6MW5rqqIs0EPg7sQ+gIjwGu8jGJ/PFgcD1TkVYEdgIO\n",
       "BtYAjgUWJWbPRC3M9VRFmgbsDXwK+BMhIM5NzPI2cbBveTC4rkt3J+wLfBK4i7Ah+L/E7IWohbmo\n",
       "KtJUwrjSQcD6wDeBE303YnweDK5r6sYPLgSOTcx+F7cql0cV6fWEgJiHj0NE58HgOsrHD9xk+DhE\n",
       "PuQ2GCTNA44DpgInmtlRdY/vCHwJeCH9+oyZXd7gdTwYesDHD1wn+ThEXLkMBoV9j7cB7wQeBK4F\n",
       "djOzWzPLrGrpRkfSa4GfmNkrG7yWB0MX+fiB6yYfh4gjr8EwF0jMbF56+zAAMzuyxfLHmtlgg8c8\n",
       "GLrAxw9cr9WNQ5wMfMPHIbqj1XYz5oV6ZgD3Z24/ALyhfiFJ7wO+BqwLbNOb0vpXk/GD1/r4geuF\n",
       "xOxaYPfMOMR1FekKfByip2J2DB8A5pnZvuntPYE3mNknmiz/ZsI4xMYNHvOOYZJ8/MDlUd04xGOE\n",
       "gDjHxyEmL6+7kgaBBZldSYcDL9QPQNc95y5gjpk9Xne/AZXMXcNmNtz5qsunyfjBhYnZ0qiFOZfh\n",
       "4xCTJ2mIcI6ymiSPwbACYfD5HcBDwDUsP/i8IXC3mZmkLYGzzGzDBq/lHcM41Y0f/IwwfnB93Kqc\n",
       "G5uPQ3RGLjsGAEnbsuxw1ZPM7GuS9gcws+MlfRb4EPAc8DRwsIV9kPWv48HQBp9/4MokMw7xEcLZ\n",
       "XX0cYhxyGwyd4sEwtkoI4S8Rxg+OI4wfPB23Kucmr8E4xBcSs0viVpV/Hgx9LD13/rGEDuFgwkVV\n",
       "fPzAlU46DrEjoXP4BXCIj0E012q7OaXXxbjeqUjbES6C8wyweWJ2noeCK6vEbGlidi7hAkLPA7+v\n",
       "SO+OXFYhecdQQnVdwj6J2RWRS3Ku5yrSO4ET8e6hIe8Y+kiDLsFDwfWlxOwXePcwId4xlIR3CRM3\n",
       "QFXABsDc9GszwpFyebKUEPhXAVXgDyMMFv+Pt0e8e1ieDz6XXNolHA/8FDjUjzZqbYDqqsDWLAuC\n",
       "QcKnyqvSrxuAf0YrsLGVgS1YVvMUQkDUar5uhMG/xSsv/yrSasB/A+8B9k3MLo5cUlQeDCXlXcLY\n",
       "GnQDc4GNWfbp+yrgqhEG72/6IjmT/kyzWBZqc4HXALeyrKO4Cu8qGvLuIfBgKCHvEhpLu4HXs2yD\n",
       "OZcwQfKqzNfvRhh8NlqRXTBAdRVgS0YH4FRG/9zeVaQq0uqE7mFb+rR78GAoEe8SliljN9ApY3QV\n",
       "2V1Qfd1V9HP34MFQEuns5RPo0y5hgOoUwqnZ30ofdQOd0kZXMQxc029BUdc97JOYXRq5pJ7wYCiB\n",
       "irQH4T/vHv3UJaRhMAjsQjgt+J+BS+jjbqBTMl1FraN4N7AqcDZwJnB1P4VERXoX8CPgwMTsjNj1\n",
       "dJsHQ8FlQuFdidktsevptjQM5gI7sywMzgTOGmFwcczayiwNis0I631nYBohJM4ihETpL+VakTYn\n",
       "fPAofTh4MBRYv4SCh0G+9HNI9Es4eDAUVNlDId34vBEPg1xrEBKrEQLiLKBaxt1N/RAOHgwF1Aeh\n",
       "MBs4CVgPOBUPg8IYoLoZYcxnN+APwL4jDN4Xt6rOK3s4eDAUTJlDIf30uT/wZeB/gK+PMOjX7y2g\n",
       "AaorAp8lXAfhcOCksnUPZQ4HD4YCKXkozCZ0CasB871DKIcBqq8BFgKPU8Luoazh4GdXLYiyhsIA\n",
       "VQ1QPQC4lvAH9kYPhfIYYfBmwoEDvwKuH6D6kbQzLIXE7CZgG+AbFWnX2PX0gncMOVHiUJiNdwl9\n",
       "o8zdQ9k6B+8Ycq6MoeBdQn8qc/fQT52DdwyRlTQUZuNdQt8ra/dQls4htx2DpHmSlki6Q9KhDR7f\n",
       "Q9KNkm6S9BuFX0hplC0UvEtwWWXtHvqhc4jWMUiaCtwGvBN4kLAx2c3Mbs0sMxdYbGZ/ljQPWGBm\n",
       "gw1eq3AdQwlDYTbeJbgmytg9FL1zyGvHMAe408zuMbPngNOBHbMLmNlVtuw0uFcTJkMVXplCwbsE\n",
       "144ydg9l7hxiBsMMIHtmzAfS+5rZB/hZVyvqgZKFwmzCeew/DLx1hMGjfLKaa2aEwedGGPwK8Dbg\n",
       "AOCiAaqzIpc1KWUNhxUivnfb+7AkvY2w8fn3FsssyNwcNrPhCVfWJWUJBZ+97CZjhMGbB6jOJcya\n",
       "vn6A6ueAE4s6azoxu6kibQNcUpHI624lSUPAUFvLRhxjGCSMGcxLbx8OvGBmR9UttzlwLjDPzO5s\n",
       "8lq5H2MoUSjMxscSXIeUaeyhaGMOeR1juA7YSNJsSSsBuwLnZxeQNIsQCns2C4UiKEMo+FiC64YG\n",
       "Yw/7FnXsoUy7laLOY1C4VOVxhMsLnmRmX5O0P4CZHS/pROA/gNqniOfMbE6D18ltx1CSUJiNdwmu\n",
       "y8rSPRSlc/CT6EVS9FDwsQTXa3VnbC3s2EMRwsGDIYIShMJsvEtwkZShe8h7OHgw9FhF2g34OgUM\n",
       "Be8SQKpMBXLz/yllZsnS2EX0Uhm6h0w4fDwxOzt2PVkeDD1UkbYALgXelpjdHLue8ejnLkGqTAfe\n",
       "T7h05ZvIYTAAvyZc+vRcs+TRyPX0TNG7h4r0OsJ8n1xtEzwYeqQirU442ipJzE6LXU+7+rVLqAuD\n",
       "LYALCdcxvtgseTZmbfWkyirAPEKt7wGuJ9TaFyFR9O6hIu0FHAa8PjF7OnY94MHQMxXpFOCZxGz/\n",
       "2LW0a4DqWoRPoavTB12CVBHwn8B+5DwMmmkSEt8HzjRLiv8H3UJd97DrCINPxa2ofRXphwCJ2d6x\n",
       "awEPhp6oSK8ArgFmJWZ/j1lLuwaoTiHMHbkP+GTZuwSpMoOwAZ0BLAAuKkoYNJMJiQrwB+AAs+SP\n",
       "cavqrrR7+C6wNvCBonQOFWlVwt/aFonZ/WMt3215neBWNh8CTitKKKQOJvxxHVjmUJAqkirzgRsI\n",
       "E/TmmCXnFT0UAMySv5slPwFeD9wE3ChV9kw7o1IaYfA54OPALOATkctpW2L2DKE7/WDsWsbiHUMH\n",
       "VKQpwJ3ATonZ72LVMR4DVN8I/ASYM8LgvbHr6Za6LmG+WXJD5JK6SqpsRdjVcjcl7x4GqG4AVIHt\n",
       "Rhi8NnY97aiEUwGdDGycRN74esfQfW8GngZGYhfSjgGqawOnAR8payg06RJKHQoAZsn1wNb0Qfcw\n",
       "wuDdwEeBMwaorhG7njZdDbxAOA1IbnnH0AHpoNLvE7NjYtXQrsy4wpIRBj8du55u6LcuoZl+6R4G\n",
       "qH6TcK2WQow3VKTDgA0Ss/1i1uEdQxdVpGnA+4Afx66lTbVxhcNjF9Jp/dolNNNH3cNnKNZ4w4+A\n",
       "nSrSi2MX0ox3DJOUHp/8gcRshxjvPx5lHlfwLqG1sncPRRtvqEgXAScnZqfGqsE7hu6aT/iDy7Wy\n",
       "jit4l9CesncPBRxvWEjYduSSdwyTkJm7MCMx+2ev3388BqieBPx1hMFPxa6lU7xLmJgydw8DVL8L\n",
       "TBlh8IDYtbRSkVYhXM442pwG7xi6pzZ3Ie+hsAGwI2ESVOF5lzA5Je8evgjsMkB1/diFtJLOdzqT\n",
       "nM5p8I5hgoo0dyHtFh4YYTCJXctkpV3C8YSjULxLmKQydg8DVL8GrFmAruENhIHoKHMavGPojkLM\n",
       "Xch0C8fFrmUy0i5hL8L6vg7vEjqirnu4QarsUYLu4esUoGsg7IZeSg7nNHgwTNx8YGHs2YttOAL4\n",
       "zgiDT8YuZKLSLuEC4CBgG7NkgVmS6913RWKW/MMs+QLhhHyHAedJlXUjlzVhIww+Rugqc31Idrrt\n",
       "WEgOB6E9GCagKHMXit4teJfQWyXrHorSNZxCDuc0eDBMzAeAXydmj8QuZAyF7Ra8S4ijLN1DgbqG\n",
       "BwmnyXhf7FqyogaDpHmSlki6Q9KhDR7fRNJVkp6VdEiMGpuYT87nLhS1W/AuIR9K0j0UpWtYSM52\n",
       "J0U7KknSVOA24J3Ag4TDDnczs1szy/wbsD4hTZ80s683ea2eHZVUlLkLRTwSyY84yqciH7lUhCOU\n",
       "KtKLCNvAns5pyOtRSXOAO83sHjN7Djid8An3X8zsT2Z2HfBcjAKbyP3chaJ1C94l5FvBu4fcdw2J\n",
       "2bPkbE5DzGCYAWTT8YH0vtxK5y7sRc53I1GgsQUfSyiGoo49FGWsgXR3UkXKReDGDIaO7sOStCDz\n",
       "NdTJ187I/dyFInULUuXNeJdQKA26h9wdg99A7rsGejCnQdJQdjvZatkVulVEGx4EZmZuzyR0DRNi\n",
       "ZgsmW1Ab5pP/uQuF6BakyjqEk/rtZZb8PHY9rn1myT+AL0iVa4AzpMqWZsljsetqZoTBxwao1rqG\n",
       "XI41JGZWkRYStjG/7cZ7mNkwMFy7Lanp+GPMjuE6YCNJsyWtBOxKuIBMI9HbqyLMXShKtyBVphKO\n",
       "317ooVBcZskFhLHBRVIl74e+F6FryM2chgn9MiVN+o/ZzJ4nXND7YmAxcIaZ3Sppf0n7p+8zXdL9\n",
       "hP3Pn5d0n8IGOoYizF0oRLdA+OS2ErAgch1u8o4A1gRyfTXAIow15GlOQ9PDVSVt2ew5wIVmNr1r\n",
       "VY1TLw5XrUhXAN9KzM7t5vtMVNotXANslOdgkCpDhF1IW5klD0Uux3WAVJlFONz8/WbJb2LX08wA\n",
       "1ZcCtwMDeb0mSUXaFdgnMdum2+/VarvZKhiWAr9q8pqDZrZKh+qbtG4HQ0WaSRgknZGY/aNb7zMZ\n",
       "6XnoHxth8Iuxa2lGqqwM3AV82Cy5JHY9rnOkynbA94ANzZI8HV4+ygDVI4FVRhg8MHYtjWTmNLwm\n",
       "MevqfJGJzmNYAuxvZm+r/wJyO9DUJS8D7s1rKKRmETqGPFsdWNlDoXzMkgsJv99VY9cyhqsJfyu5\n",
       "lM5peICwzYmmVTAsaPH4JztfinPOuTxoeriqmZ3V4rGfdKcc55xzsY15VJKkr0laM3N7TUn/1d2y\n",
       "nHPOxdLO4arbmtm/jnJJv9+ueyU555yLqZ1gmKIwUg6ApFUIx6A755wroXZOifFj4DJJPyDMYdgb\n",
       "OLmrVTnnnItmzGAws6Mk3QS8I73rS2Z2cXfLcs45F0tbJ9Ezs58Dfk4b55zrAxM9V9IJnS7EuTKS\n",
       "GIpdg3Pj1c7hqsdI2qzu7uO7VI9zZTMUuwDnxqudjuFW4PuSrpF0gKSXpJfbdM45V0LtDD6fAJwg\n",
       "aRPCRSR+L+lK4AQzu6LL9TlXOOnuo6H0ZpK5WOOw2bILpTiXV20NPkuaCmwCbAr8CbgROFjSAWa2\n",
       "axfrc65w0o3/MIAEZn7dCVcsYwaDpGOB7YHLga+YWe0MnkdJuq2bxTnnnOu9djqGm4DPm9kzDR57\n",
       "Q4frca5shmMX4Nx4NQ0GSVsBRgiGjbVsR6kAM7PfmdlT3S/RueLyMQVXRK06hq8TgqGZt3W4Fuec\n",
       "cznQ6noMQz2swznnXE60M8FtJUkHSjon/fqEpBU78eaS5klaIukOSYc2Weab6eM3ShroxPs655xr\n",
       "rp0Jbt8DtgS+A3wX2Cq9b1LSQ2C/DcwDXg3sJmnTumXeA7zSzDYC9uvE+zrnnGutnaOSXm9mm2du\n",
       "X5aebXWy5gB3mtk9AJJOB3YkzLSu2QFYBGBmV0taQ9I6ZvZIB97fOedcA+10DM9LemXthqQNgec7\n",
       "8N4zgPsztx9I7xtrmfU68N7OOeeaaKdj+AxwuaS7CYeqziZcrGeyWh3xlKW62+0+zznn3AS0Ewy/\n",
       "Bb4PvB14Crg4vW+yHgRmZm7PJHQErZZZL71vOZIWZG4Om9nw5Ev8l+eBaR18vW54DlgrdhFjeA54\n",
       "sVR5qVnyWOxiXOdIlXUIl/ztxN6EblqL8P8wlyphwtg0ulCjpCHaPNtvO7uSTgZeAXyZMFi8AfCj\n",
       "CdaWdR2wkaTZklYCdgXOr1vmfOBDAJIGgaeajS+Y2YLM13AH6su6GVi1Ir22w6/bSScDnxqgWt9h\n",
       "5YZZ8hTwLeBkqTKha4G4/JEqU4FTgP8xS56OXU8zA1SnAJ8iHbfMqS0I2+UlnX5hMxvObidbLdvO\n",
       "H+dmZraPmV1hZpeb2UeA+uszTKTI54GPEzqQxcAZZnarpP0l7Z8u8zPgbkl3Eq4B8f8m+74TkZgt\n",
       "JWx494rx/m36KeH3uX3sQsbwBWB1wi5KVw6fI+x9qMQuZAzvB54Ffha7kBbmA4sSsxdiFiGz1rvs\n",
       "JZ0CfMfMrkpvDwIfM7MP9qC+tkgyM+vqJ+WKtDHwS2BmYpbLVnSA6n8QNrxbjTCY27EYqTITuBbY\n",
       "ySy5MnY9buKkytuAU4GtzJKHYtfTTNot3AgcNsLghbHraaQS9pw8AAwmZnd3+/1abTfb6Ri2Bn4j\n",
       "6V5J9xDGF7aW9PsOHbZaCInZbcDdwLtj19JCIboGs+R+YB/gVKny0tj1uIlJxxVOAfbKcyikitAt\n",
       "bAcs7kUojKWdYJhHGFd4K2HgYgNgW8LGZ4euVZZPCwmtXi6NMPgCoZ1fkOexBgCz5ELgNOBCqfLK\n",
       "sZZ3+SJVNgYuBH5ollwSu55W0m4hARbkuZMmbFsWRq4BaCMYzOyeVl89qDFPzgTeVZHWjl1IC4Xo\n",
       "GlKHE8KhKlUO9AHp/JMqU6XKIcBvCBuxL8atqC257xYq0jqED99nx64F2usYXCoJpxm/ENgtdi3N\n",
       "FKxreMEsOQ6YC+wMDHv3kF9pl/ArwoeON5gl3zZLog6SjqVA3cLuwHmJWS6O6vJgGL+F5Hh3UqpI\n",
       "XQNmyR2ET0vn4t1D7tR1CacBbzdL7opcVruK0C2IMGl4YeRS/sX/+MbvMmB6nuc0FKlrqDFLlnr3\n",
       "kD9F7BJqCtQtbAGsRljPueDBME4FmdMABesaarx7yIeCdwk1ue8WUvPJwdyFrDHnMRRBL+YxZBVh\n",
       "TgMUZ15DM1JlI+CHwAvAh82SOyOX1BfSLuEHhNMy7FPAQCjEvAXo/dyFrMnOY3B1CjKnAQraNdR4\n",
       "99BbJekSaorSLeRm7kKW/5FN3EJyPghdxLGGeg3GHn7pYw+dV+SxhHoFGluAHM1dyPJgmLgizGmA\n",
       "gncNNZnu4Ry8e+iYknUJNYXoFvI2dyHLxxgmoSKdCvw2Mft2r997PIo+1lAvM/ZgwN4+9jAxZRhL\n",
       "qFeUsQWAinQQ8LrEbH6M9/cxhu5ZSM53J6VqXcP7YhfSCd49TE5Ju4SanSlGt5C7uQtZ/sc0Obmf\n",
       "0wD/Gms4GPjOANXpsevpBB97mJgyjSXUG6C6HvAN4OACdMa5m7uQ5cEwCQWa08AIg5cDJwCnDlCd\n",
       "GrueTvHuoT0l7xIYoLoC4ef61giDv45dTxvmk7O5C1n+BzR5i4A9K9KKsQtpw5cI19D+QuxCOsm7\n",
       "h9bK3CVkfAn4O/C12IWMJZ27sDvhQ2UueTBMUoHmNDDC4FLCf8j9Bqi+I3Y9nebdw2hl7xJqBqjO\n",
       "I1wCeM90t2ne5XLuQlbf/tF02EKKMQjNCIN/JPwR/ags4w1Z3j0EfdIl1MYVFgJ7jDD4aORy2jWf\n",
       "nA4613gwdEZR5jQAMMLgLyjheENWv3YP/dIlwHLjCr+MXU878jx3Iav0fyi9UITrNDRQyvGGrH7r\n",
       "HvqlS8gozLhCxu7ATxOzv8YupBUPhs5ZSEF2J0H5xxuyGnQPn5IqpemU+qlLqCnguELu5y5kRZn5\n",
       "LGkt4AxgfeAeYBcLn7rrl/sBYaDmUTNrOlcg1sznrIo0FbgdOCgxOz9mLeORhsLpwOeB7xfg+O9J\n",
       "SWdNfx/YhBAUZwFXmiVLoxY2TmmwvYXQCb0fuAXYrw8CYQrwUcI5wD5QlF1IABVpJ+ArwKZ5OEy1\n",
       "1XYzVjAcDTxmZkdLOhRY08wOa7Dcm4GngZPzHgwAFWkucB4wJzG7N3Y97RqguhnhFBN/Bj4ywmBh\n",
       "ap8oqfIqwkZ1Z2AdChASDcLgIULNZ/XDaUEGqL6CcAqPFwF7jzC4JHJJbatIGwBVYLvE7NrY9UA+\n",
       "g2EJ8FYze0TSdGDYzDZpsuxs4IIiBANARfo0sBPwlsTsn7HraVc6kPdp4BD6pHuoaRASvwYei1rU\n",
       "8lYlHBLdV2EAy3UJRwHHpLtCC6EirUzYzXdyYvbN2PXU5DEYnjSzNdPvBTxRu91g2dkUKximED59\n",
       "Pgfsnpg9H7mkcenH7iErDYk5wOqxa6nzPHB5v4RBTZG7BIB04usZhN/frkmOzloaJRgkXQo0Ok7+\n",
       "CGBRNggKtQzvAAAOOElEQVQkPWFmazV5ndm0EQyETxM1w2Y2PIGyO6IivYhwcZmnKWY49G334PKh\n",
       "6F0CjAqFFYCdE7N/xKxH0hAwlLkryVvHsAQYMrOHJa0LXFGWXUk1RQ8H8O7BxVH0LgHyFwqN5PG0\n",
       "2+ez7MRzexEGbEslMXuWMEA4DTi1Iq0QuaRxG2HwFuCNhLPIXjdAdf+iXgnO5d8A1SkDVD8GXEs4\n",
       "bfabPBTiiHm46pnALDKHq0p6OXCCmW2XLnca4fjztYFHgS+a2Q8bvF7uOoaaMnQO4N2D664ydAlQ\n",
       "rFDI3eBzp+U5GKBU4eBjD66jyjCWUFOkUAAPhlwoSziAdw+uM8rSJUDxQgHyOcbQd8ow5lDjYw9u\n",
       "MsoyllBTxFAYi3cMPVamzgG8e3DjU6YuAYodCt4x5EiZOgfw7sG1p65LuJCCdwlQ7FAYi3cMkZSt\n",
       "cwDvHlxjZesSoByh4B1DDpWtcwDvHtxoaZfw/4BrKEmXAOUIhbF4xxBZGTsH8O6h36VdwknAKpSk\n",
       "S4ByhYJ3DDlWxs4BvHvoV3VdQuGPOMoqUyiMxTuGnChr5wDePfSLsnYJUM5Q8I6hAMraOYB3D2VX\n",
       "5i4ByhkKY/GOIWfK3DmAdw9lU+YuAcodCt4xFEiZOwdo2D18doDqjMhluXEaoLreANXDKGmXAOUO\n",
       "hbF4x5BTZe8cAAaovhr4LLADsJhwxt1zRhh8MGphrqEBqusRLlu7M7AJ4fT5R5UtEKA/QsFPoldQ\n",
       "/RAOAANUVwbeSdjgeEjkSIMw+CnhmtOXjTBYmGuaj0c/hAJ4MBRav4RDjYdEfP0YBjX9EgrgwVB4\n",
       "deHwwTL/Z81qEhJnAZcAt40w+ELE8kojvSbCJsA29GEY1KR/Z6fSB6EAHgylkP6nPQXYGJifmF0f\n",
       "uaSeyoTEToSr+q0BXA1cBVSBq0cY/HO8CotjgOoawBuAucBg+v3jwC+Bs+mjMKipSHMIR8v9Htir\n",
       "7KEAHgylUZEE7A4cA5wAfLkf/gM3MkB1HZZt2OYCWxEuE3tV5qvvu4pMNzA38zULuJ5l66k6wuCj\n",
       "0YqMKP3AlQAfBg4EzkjKsFFsgwdDyVSkdYH/BTagD7uHRgaorghszugNYN91FU26gScYHZg3jTBY\n",
       "6rGqdmS6hNuAjyZmj0Quqac8GErIu4exDVCdzrKOYi6wJSXqKrwbmJh+7hKychcMktYijPyvT/hD\n",
       "3cXMnqpbZiZwMvAywIDvm9k3m7xe3wVDjXcP7St6V+HdwOT1e5eQlcdgOBp4zMyOlnQosKaZHVa3\n",
       "zHRgupndIGka4VPQ+8zs1gav17fBAN49TEaTruIx4G8x62pgZcKHJO8GJsC7hOXlMRiWAG81s0fS\n",
       "ABg2s03GeM55wLfM7LIGj/V1MNR49zB5aVcxi7AhzpOlwN0jDD4Xu5Ci8S6hsTwGw5Nmtmb6vYAn\n",
       "arebLD+bcCjdZmb2dIPHPRhS3j04F3iX0FqUYJB0KTC9wUNHAIuyQSDpCTNbq8nrTAOGgf8ys/Oa\n",
       "LGNAJXPXsJkNT7D0UvDuwfUz7xKWJ2kIGMrcleStY1gCDJnZwwobsCsa7UpSmJ7+f8DPzey4Fq/n\n",
       "HUMD3j24fuNdQvvyuCvpaOBxMztK0mHAGg0GnwUsSpc7aIzX82BoIe0evgtsDXwLOCExezJuVc51\n",
       "TiUc6bg/8DHC0WUf8y6htTwGw1qEE6PNInO4qqSXAyeY2XaS3gT8CriJcLgqwOFmdlGD1/NgaENF\n",
       "2gI4CNge+DHwjcTszrhVOTdxFelVhM5gN8L5nY5LzG6MW1Ux5C4YOs2DYXwqIYA/BuwHXEnY1XSl\n",
       "t9yuCNJdpG8lfMiZCxwPfDcx+2PUwgrGg8E1VJFWBT4EfAr4KyEgzkrM/JBIlzsVaSVgF+BgwqVE\n",
       "jwVOSczyNuekEDwYXEsVaQrwHsIf3Eb4OITLkbrxgyWEDzAXJWaFPJVJXngwuLb5OITLCx8/6C4P\n",
       "BjduPg7hYvDxg97xYHAT5uMQrhd8/KD3PBjcpPk4hOsGHz+Ix4PBdVRFGiC0+u/FxyHcBPj4QXwe\n",
       "DK4rfBzCjYePH+SLB4Prqsw4xEHAnwkBcbaPQzjw8YO88mBwPZGOQ2xHCAgfh+hzPn6Qbx4Mrufq\n",
       "xiFOIYxD3BW3KtcLPn5QDB4MLpoG4xAnAZf66b/LJT3d9TbAPvj4QSF4MLjo0nGIPYE9gNcAFwBn\n",
       "Ab9IzJ6NWZubmDQM3g3sTNiFeCPhKLUf+/hB/nkwuFxJu4j3EwYktwbuBhYDt6Rfi4E7fPA6H9LB\n",
       "442AzdKvV6f/vgK4mhDw5yRmD0cr0o2bB4PLrfRT58Ysv9GZCdyFB0bPjBEA97Hsd1D7fdzu3V5x\n",
       "eTC4whkjMO5mdFjcggdG2zwAHHgwuBLxwGifB4BrxYPBlV4/B0ZdANR+bg8A15IHg+tbZQoMDwDX\n",
       "SR4MztUZIzAeBPIWDisCM/AAcB2Su2BQmCp/BrA+cA+wi5k9VbfMi4BfAisDKwE/NbPDm7yeB4Pr\n",
       "iDQw1gemxK6lzlLgPg8A1yl5DIajgcfM7GhJhwJrmtlhDZZ7sZn9TdIKhFmznzazKxss58HgnHPj\n",
       "0Gq7GetT0Q7AovT7RcD7Gi1ky2ZPrgRMBZ7ofmnOOdffYgXDOmb2SPr9I8A6jRaSNEXSDekyV5jZ\n",
       "4l4V6Jxz/WqFbr2wpEuB6Q0eOiJ7w8xMUsP9WRZOz7uFpJcAF0saMrPhJu+3IHNzuNlyzjnXjyQN\n",
       "AUNtLRtpjGEJMGRmD0tal9ANbDLGc74A/N3M/qfBYz7G4Jxz45DHMYbzgb3S7/cCzqtfQNJLJa2R\n",
       "fr8K8C5gpGcVOudcn4p5uOqZwCwyh6sqnHXzBDPbTtLmwEJCeE0BfmRm/93k9bxjcM65ccjd4aqd\n",
       "5sHgnHPjk8ddSc4553LKg8E559woHgzOOedG8WBwzjk3Sl8HQzrho3C87t7yunvL646vr4OBNmcB\n",
       "5tBQ7AImaCh2ARM0FLuACRqKXcAEDcUuYIKGYhfQKf0eDM455+p4MDjnnBulNBPcYtfgnHNFU+qZ\n",
       "z8455zrHdyU555wbxYPBOefcKKUMBkkzJV0h6RZJN0v6ZHr/zul9SyVt2eL58yQtkXRHek3qotR9\n",
       "j6SbJI1IuiZyzf8t6VZJN0o6N73YUqPn521dt1t3z9f1GHV/Oa35BkmXSZrZ5Pl5W9/t1p2r9Z15\n",
       "/BBJL6RnjG70/Cjre9LMrHRfhCvHbZF+Pw24DdgU2AR4FXAFsGWT504F7gRmAysCNwCb5r3u9Dl/\n",
       "ANbKybp+FzAlvf9I4MiCrOsx6461rseoe7XMMp8ATizI+h6z7jyu7/T2TOCiZrXFXN+T/Splx2Bm\n",
       "D5vZDen3TwO3Ai83syVmdvsYT58D3Glm95jZc8DpwI7drTiYZN01PT39eIuaL7VwaVaAq4H1Gjw9\n",
       "j+u6nbpren6q9xZ1/zWz2DTgsQZPz+P6bqfumtys7/ThY4DPtnh6tPU9WaUMhixJs4EBwh95O2YA\n",
       "92duP5De11MTqBvAgF9Iuk7Svt2oq5UWNX8Y+FmDp+R9XTerGyKva1i+bklfkXQf4aqIRzZ4Si7X\n",
       "dxt1Q87Wt6QdgQfM7KYWT8nF+p6IUgeDpGnA2cCBadq3I/rxuxOsG+DfzWwA2Bb4mKQ3d6XABprV\n",
       "LOkI4J9mdmqDp+V2XY9RN0Rc19C4bjM7wsxmEa58eGyDp+VyfbdRN+RofQMvAJ8DkuwiDZ4WfX1P\n",
       "VGmDQdKKwDnAKWa23DWlW3iQsO+wZiYh6XtiEnVjZn9M//0T8BNCK9t1zWqWNB94D7BHk6fmcl23\n",
       "UXe0dQ1t/R85FXh9g/tzub4zmtWdt/W9IWHc4EZJfyDsbrxe0svqnhp1fU9K7EGObnwR0vtk4Ngm\n",
       "j18BbNXksRWAuwi/+JXo7QDdZOp+MelAHrAq8Btgm1g1A/OAW4CXtnhu7tZ1m3VHWddj1L1R5vtP\n",
       "EK6RXoT13U7duVvfdcs0G3yOtr4n/XPHLqBLv8w3Edq9G4CR9Gtb4H2EfX5/Bx4Gfp4u/3Lgwszz\n",
       "tyUcfXAncHgR6gY2SJ93A3Bzr+puUfMdwL2Z+75bkHU9Zt2x1vUYdZ8N/D69/xzgZQVZ32PWncf1\n",
       "XbfM3aTBkJf1PdkvPyWGc865UUo7xuCcc25iPBicc86N4sHgnHNuFA8G55xzo3gwOOecG8WDwTnn\n",
       "3CgeDM51iKRhSVtN4vlrSbpU0u2SLpG0Rifrc65dHgzOdY4xufPjHAZcamavAi5LbzvXcx4MzjUh\n",
       "aXZ6kZVTJC2WdJakVSS9Q9Lv0gvHnCRppdFP096Sjs3csa+kY+pe+9D0+TdI+mp69w7AovT7RYQZ\n",
       "7871nAeDc629CviOmb0a+AtwCPBDYBcz25xwPpyPZpY34Exge0lT0/vmAyfVFpC0LSEE5pjZFsDR\n",
       "6UPrmNkj6fePAOt05SdybgweDM61dr+ZXZV+fwrwduBuM7szvW8R8JbsE8zsGeByQjhsAqxoZrdk\n",
       "FnkH8AMzezZd/qn6N7Vwrho/X42LwoPBudayG2cBTzH63PvNrip2IrA3oVv4QYPHGz3vEUnTASSt\n",
       "Czw63mKd6wQPBudamyVpMP1+d+A6YLakDdP7PggMZ5YXgJldQzhP/+7AaXWveSmwt6RVACStmd5/\n",
       "PuEqZqT/jut6HM51igeDc63dRrhi2GLgJYTr/O4NnCXpJuB54H8zy2c7jDOBK83sz5K2lnQCgJld\n",
       "TAiB6ySNEMYtIFzW8l2Sbifssmp2mUvnuspPu+1cE+k1fi8ws9dO8PkXAMeY2RWdrMu5bvOOwbnW\n",
       "xv3JSdIakm4D/uah4IrIOwbnnHOjeMfgnHNuFA8G55xzo3gwOOecG8WDwTnn3CgeDM4550bxYHDO\n",
       "OTfK/weWtLLMB1vVbgAAAABJRU5ErkJggg==\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sherpa.reg_proj(poly.c0, poly.c1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can also directlty inspect the parameter space. For instance, in the plot below Sherpa displays the Interval Projection of the `c0` parameter, i.e. a plot of the error for each value of the parameter, around the minimum found by the optimization method during the fit."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecXGXd/vHPRULvEaUJRjpq6P1ByUMRUEEFFUhUuvwA\n",
       "AR/lERCQoCKCgFQRQhGEoDQVHlAEYUGKQDCEZpDei4AkdAL5/v64zyaTZWd3ZqfcU67367UvdmbO\n",
       "nLnmZJnv3OfcRRGBmZlZf+bIHcDMzFqXi4SZmZXlImFmZmW5SJiZWVkuEmZmVpaLhJmZleUiYdaH\n",
       "pBmSlmvC64yVdE277Ne6k4uENY2kxyVtVsF2PZJ2b0amakkaXRSR1yRNkzRF0i5D2VdEXBgRW9aY\n",
       "Z2SRZ+b/y/XYr1kvFwlrpih+KtluyEo/MBvkmYhYMCIWAg4CxktatZ8cwxqcY7aXa+JrWRdxkbBm\n",
       "k6RdJN0s6eeSXpH0qKStigePAj4NnFp8Wz+5uH8VSddKern49v7Vkh3+WtLpkq6S9DpwoKTnSouF\n",
       "pC9Lmlz8vp6k2yT9R9Kzkk6RNOdQ3kxE/BH4D7Bq8b5ukXSCpJeAIyQtJOl8SS8WLalDJanIsYuk\n",
       "v5VkHOg9zivp+GIfr0q6SdI8wE3FJq8WLZsN+tnvRpLuLJ53h6QNSx7rkfSj4t9jmqRrJH1oKMfC\n",
       "OpOLhDVbbythPWAK8CHgWOBsgIg4FPgbsG/xbX1/SfMD1wIXAB8GdgR+2efb+07ATyJiAeAk4A1g\n",
       "05LHxwAXFr+/BxxQvPaGwGbAPtW+EUlzSPoysAhwb8n7egT4CPBT4FRgQeDjwCbAN4Fd+9nXYO/x\n",
       "OGDNIu8I4PvADFJBBVg4IhaKiL/32e8I4CrgxOJ5JwBXSVq0ZLOdgF2KzHMBB1Z7LKxzuUhYLk9E\n",
       "xNmRJg87H1hS0kdKHi89ffIF4LGIOC8iZkTE3cDlwFdLtvlDRNwGEBHvABeRPvyQtCCwdXEfEfGP\n",
       "iLij2NcTwJmkD/BKLSXpP8C/gcOBr0fEQ8Vjz0bEaRExA5gO7AAcEhFvFK91PPCNfvZZ9j0WLaJd\n",
       "gQMi4rni8b9HxLsMfprp88CDxXWKGRHxW1Jx3rZ4PIBzI+LhiHgbuBhYo4pjYR1ueO4A1rWe7/0l\n",
       "It4szsAsALzYe3fJth8D1i8+mHsNJxWX3m2f7rP/CcCtkvYGtgPuioinACStRPpGvTYwX7GviX0D\n",
       "SloWuH9WzFio+P3ZiFimzPt6quT3xYA5gSdK7nsSWLqf5w30Hj8EzENqoVRrqeI1Sz1R3N/r+ZLf\n",
       "3yL9O5gBbklYa+p74fpJ4MaIWLTkZ8GI2LfsDiL+Sfow3Jp0qmlCycOnAw8AK0TEwsCh9PP/QkQ8\n",
       "WbzOgiUFoprsL5FaEyNL7luWDxY0GPg9vgy8DawwyOv15xlSASr1seJ+s0G5SFizVdIL5wVg+ZLb\n",
       "/wesJOnrkuYsftaVtMog+5wAfId03v6SkvsXAF4D3iz2sXdV76BCEfE+6fTNUZIWkPQx4H9I1x36\n",
       "uooy77E4dXUOcIKkJSUNk7ShpLlIp7xmMPvxKvWnYr87SRouaQdgFdIx7eWeUVaWi4Q1W9B/V9jS\n",
       "2ycBXyl6Pp0YEa8DnyVdzH0GeA44mnSRtXSffV0EfAb4a0S8UnL/gaTWxTTS9Yjf9nn+YN/Oyz3e\n",
       "X479SBfRHyVdkL8QOPcDT4x4jYHf44Gki+N3kloWRwOKiDeBo4BbiuO1fmmOiHiZdL3je6SWzYHA\n",
       "F/ocj77v3YvM2Exq5KJDkg4Bvk76pnMvsGtxUbF0m9HAL0jnbl+KiNENC2TWQiTtBoyNiEEHGJrl\n",
       "0rCWhKSRwJ7AWhExChhG+pZUus0iwGnANhHxKeArjcpj1oI+SWphmLWsRvZumka6aDefpPdJvUj6\n",
       "XiwbA1wWEU8DRMRLDcxj1jIk/YF0HeGrg21rllPDWhLFOc/jSb02ngVejYjr+my2IjBC0g2SJkrq\n",
       "r/+4WceJiC9FxKiImJI7i9lAGnm6aXlSz5KRpD7ZC0ga22ezOYG1gM8BWwKHS1qxUZnMzKw6jTzd\n",
       "tA5wa9G7AkmXAxsxa2oESAOPXoqIt4C3JN0ErA48VLojSe5tYWY2BBFRUxfnRhaJKaSWwbykgUCb\n",
       "A3f02eaPpInchgFzA+uTRsJ+QK1vtFNIGhcR43LnaAUNORYT1APAmPbqZee/i1l8LGapxxfshhWJ\n",
       "iJgs6XzSdAczgH+QplTeq3j8jIiYIunPwD3FNuMj4oFGZTIzs+o0dO6miDiWNMNnqTP6bHMcaYZL\n",
       "MzNrMR5x3X56cgdoIT25A7SQntwBWkhP7gCdxEWizURET+4MrcLHYhYfi1l8LOrLRcLMzMpq6NxN\n",
       "9SIp3LvJzKw69fjsdEvCzMzKcpEwM7OyXCTMzKwsFwkzMyurrYqExKoSa+bOYWY2FBKjJEblzlGN\n",
       "tioSpEkDj8kdwjrYBPXMnL/JrP5+TJrotG20W5G4GBgl8YncQczMqiGxHLAxcEHuLNVoqyIRwTuk\n",
       "uZ/2z53FzKxK3wbOieCN3EGq0VZFovArYAeJEbmDmJlVQmJBYGfgtNxZqtV2RSKC54ErgT1yZzEz\n",
       "q9DOwA0RPJE7SLXarkgUTgL2lRo71bmZWa0k5gD2I31utZ22nbtJ4hbghAguyxTLzGxQElsDPwXW\n",
       "iqCpH7jdPnfTScABuUOYmQ3iAOCkZheIemnnlsScwKPAthFMypPMzKw8iVVIiyCNjODt5r9+F7ck\n",
       "IphO6ing7rBm1qr2A87MUSDqpW1bEul+PgQ8DKwcwYvNT2Zm1j+JRYDHgE9E8FyeDF3ckgCI4GXg\n",
       "UmCv3FnMzPrYHbg6V4Gol7ZuSaTHGAVcQzrn925zk1nH6Z23aUyMzprD2prEMNJZjh0iuCNfji5v\n",
       "SQBEcC/wT+CrubOYmRW2BZ7PWSDqpe2LROEkfAHbzFrH/rTp4Lm+OqVIXAV8WGKD3EHMrLtJrAas\n",
       "BJ0x0LcjikQE7wOn4MF1ZpbfAcAvi276ba8jikThHGBLiaVzBzGz7iTxYWA74MzcWeql7Xs3zb4d\n",
       "pwBTIzisCbHMzGYj8QNg+Qh2z50F2qB3k6RDJN0v6V5JEyTNXWa7dSW9J2m7Gl/yVGBPiXlq3I+Z\n",
       "WVWKqYL2AU7OnaWeGlYkJI0E9gTWiohRwDBgx362G0Zat/rPQE0VL4IHgbuAMbXsx8xsCLYHHo5g\n",
       "cu4g9dTIlsQ0YDown6ThwHzAM/1stx9p1PS/6/S6JwEHSLUVHDOzKh1Ah3R7LdWwIhERrwDHA08C\n",
       "zwKvRsR1pdtIWhr4InB679Pq8NJ/AeYCNqnDvszMBiWxHrAkcEXuLPXWyNNNywPfAUYCSwELSBrb\n",
       "Z7MTgYMjXT0XNZ5uAijmbD8FD64zs+bZHzi16I7fURrWu0nSDsAWEbFHcfsbwAYRsW/JNo8yqzAs\n",
       "BrwJ7BkRV/TZVwBHltzVExE95V+bBYDHgXUieLzmN2Pdw3M3WZUklgQeAJaL4D95s2g0MLrkriNq\n",
       "7d3UyDWipwCHS5oXeBvYHGafxyQiluv9XdK5wJV9C0TJtuMqfeEIXpf4NfBt4MCqk5uZVW5v4KLc\n",
       "BQKg+PLc03tb0hG17rOR1yQmA+cDE4F7irvHS9pLUjOm9j4V2KVoVZiZ1V3R3X4vOqzba6mOGkz3\n",
       "wedxOXBtxMwL42YD8+kmq4LEzsBOEWyVO0t/Wn4wXQs4Gdhf6vj3aWZNVnSzP4AObkVA5xeJG4F3\n",
       "gS1yBzGzjrMxMD9pIHDH6ujTTem57AZ8JYLP1TmWmXUxiUuBnghOzZ2lnHqcbuqGIjEPaUDfp4tp\n",
       "O8zMaiLxMeAfpGWTX8udpxxfk6hABG8D40ndYc3M6mEf4PxWLhD10vEtifR8lgbuBT4ewdT6JTOz\n",
       "biMxP2mw7voRPJo5zoDckqhQBM8A1wC75c5iZm3v68CtrV4g6qUrWhJpH2wATABW7MT5Vcys8Ypu\n",
       "r/cB+0Vwfe48g3FLojq3Ay8Bn88dxFrYBPXMHFBn9kGbATOAG3IHaZauKRLF7LAnAN/3WhNmNkTf\n",
       "B04sPk+6QtcUicIlwEfwWhNmViWJ9YFVgN/kztJMXVUkimsRRwOH585iZm3nMOCYCN7NHaSZuqpI\n",
       "FC4AlpfYKHcQM2sPEmsCawNn587SbF1XJCKYDvyM9K3AzKwShwLHFYNzu0rXdIGdfX/MDTwCfCmC\n",
       "ifXar5l1HolPAteTVp57I3eeargL7BBF8A5wLOnbgZnZQH4A/KLdCkS9dGVLIu2T+UitiS0jZq6c\n",
       "Z2Y2k8SKwK3A8hFMy52nWm5J1CCCN0njJtyaMLNyfgCc2o4Fol66tiWR9suCpNbEZyKYUu/9m1n7\n",
       "khgJ3AWsEMF/MscZErckalRM83sy6duCmVmpg4Ez2rVA1EtXtyTSvlmY1JpYP4JHGvEa1kZ6520a\n",
       "E6Oz5rCsJD4K3AOsHMG/c+cZKrck6qBYX+J00rcGMzOA/wXObecCUS9dXyQKJwLbSyybO4iZ5SWx\n",
       "BPAN4LjcWVqBiwQQwcvAWaQZHs2su30XmBDBc7mDtAIXiVmOB8ZILJk7iJnlIbEYsAdwTO4srcJF\n",
       "ohDBC8D5wIG5s5hZNgcAl0bwVO4graLrezfN/josDdxLm/doMLPqSSwCPAys1ynrV7t3U51F8Azw\n",
       "O9I5STPrLvsBV3VKgagXtyQ+8FozR1muGMErzXhNM8urmH3hUWDjCB7Mnade2qIlIekQSfdLulfS\n",
       "BElz93l8rKTJku6RdIuk1RqdaSARPA78Edg/Zw4za6q9gb92UoGol4a2JCSNJM3DvmpEvCPpd8DV\n",
       "EXFeyTYbAg9ExFRJWwHjImKDPvtpWksivR4rAreR5o/v2om9zLpBMSP0o8AWEdybO089tUNLYhow\n",
       "HZhP0nBgPuCZ0g0i4raImFrcvB34aIMzDSqCh4BrgH1zZzGzhtsTuK3TCkS9NLRIRMQrpPEHTwLP\n",
       "Aq9GxHUDPGV34OpGZqrCUcB3JObPHcSaaIJ6Zs7fZB1PYh7SFBw/yZ2lVQ1v5M4lLQ98BxgJTAUu\n",
       "kTQ2Ii7sZ9v/BnYD/qvMvsaV3OyJiJ565y0VwQMSNwF7kdadMLPOswswOYK7cgepB0mjgdF13WeD\n",
       "r0nsAGwREXsUt78BbBAR+/bZbjXgcmCriHi4n/009ZrErNdlDVLLZvkI3mr261sGngW2a0jMCTwE\n",
       "7BTBbbnzNEI7XJOYAmwgaV5JAjYHHijdQNKypALx9f4KRE4R3A1MJLVwzKyzfB14uFMLRL0MWiQk\n",
       "zS/pcEnji9srSvpCJTuPiMmkqS4mwsx1pMdL2kvSXsXtHwKLAqdLmiTpjqrfRWP9BDhIYq7cQcys\n",
       "PiSGkxYb87WIQQx6uknSxaTBZd+MiE9Kmh+4NSJWb0bAIkOW002zXp9rSPO5jM+VwZrEp5u6gsRY\n",
       "4P+Rli5u/RHFQ1SPz85KisRdEbG2pEkRsWZx3+QuKxIbk1pEK0XwXq4cZlY7iTmA+4DvRPCX3Hka\n",
       "qVnXJN6RNG/Jiy4PvFPLi7abCG4mdeMdkzuLmdVsO+B14NrcQdpBJUViHPBn4KOSJpBGUB/UyFAt\n",
       "6sfADySG5Q5iZkMjIeAw4MedfJqpnirqAitpMaB3qoy/R8RLDU31wdfPeropZUDArcCJEfwuZxYz\n",
       "GxqJbUhf+NbshiLRrGsSmwAB9L5QAETETbW8cDVaoUikHGwJnAx8KoLpufOYWeWKswD/AI6M4PLc\n",
       "eZqhHp+dlYy4/l+YWXHnAdYj9XbatJYXblN/IV2b+BZwWuYsZladbwKvAb/PHaSdVD3iWtIywEkR\n",
       "sV1jIvX7mi3RkgCQWJ1ULFaKYOpg21ubcRfYjlTMwfYgsH0Et+fO0yy5Rlw/Daxay4u2swgmA1cB\n",
       "h+TOYmYV+x5wczcViHoZ9HSTpFNKbs4BrAGdMRlWDQ4H7pE4PYIncocxs/IklgQOANbJnaUdVXJN\n",
       "orQgvAdcFBE3NyhPW4jgGYlTgZ8CY3PnMbMB/Qg4J4LHcgdpR4MWiYj4dRNytKOfA/+SWDeCO3OH\n",
       "MbMPkhgFbAusnDtLuypbJCQNtEpTRETWtahzi+B1iR8Cx0ts0g19rs3a0M+BoyJ4NXeQdlW2d1Ox\n",
       "PnVZEfF4/eP0r5V6N5Uq+l3fDRwewR9y5zGzWYpxTacCn4zg3dx5cmjKYLpW0KpFAmb+IZ5CGmDX\n",
       "lX+IZq2m+AI3CTgionvHRTSlC6ykDSXdKekNSdMlzZA0rZYX7SQRXAM8Rlrm1Mxawy7Aq+AWfq0q\n",
       "mioc2BG4mNSF7JvAyhFxcOPjzczQsi0JmHlx7DpgZZ/7NMtLYgHSwLkvR9Bqi5g1VdMG00XEQ8Cw\n",
       "iHg/Is4FtqrlRTtNBPcCV5JWujKzvA4Ebuz2AlEvlYyTeEPS3MBkSccCzzNrsj+b5XDgPolfRvB4\n",
       "7jBm3UhiKWA/YO3cWTpF2ZaEpHWLX79RbPdt4E3go8D2jY/WXiJ4jjRD7NG5s1gNJqhn5vxN1o5+\n",
       "DJzlL2r1M1BLYrykBYCLSKOsHyAtQGTlHUcaYLe+54gxay6J1YAv4IFzdVW2JRERa5AO+PvApZLu\n",
       "kXTwYOMnulkEb5BOOx1fLFJkZk1Q/P92HGnFOXceqaMBL1xHxJSIGBcRnyD1aloYuF7SrU1J157O\n",
       "AxYEvpw7iFkX2RL4GHBG7iCdpqLeTZLmAD4CLA7MD7zQyFDtLIL3Sb0rjpGYK3ces04nMZzUivi+\n",
       "V4ysvwGLhKTPSPolaQ2JA4G/AStFhL8lDyCCa4GHgb1zZzHrArsCLwNX5A7SiQaau+kp0lKdFwGX\n",
       "RES21kOrD6brj8SngOtJA+z+kzuPWScqBs79C9g2gom587Sahs7dJGlkMyfxG0g7FgkAiTOBaREc\n",
       "mDuLWSeS+BGwXARfz52lFXmCvxYnsQRwP7BuBI/mzmPWSSSWBu4B1vIKkf3Ltca1VSiC54ET8QA7\n",
       "s0b4CXCmC0RjNbRISDpE0v2S7pU0oZjeo+82J0t6SNJkSWs2Mk8mJwD/JbFh7iBmnUJiDWBr4Ge5\n",
       "s3S6SqYKP1/SoiW3R0g6p4LnjQT2BNaKiFHAMNJssqXbfA5YISJWBL4FnF5V+jZQDLA7DDjOA+zM\n",
       "alcycO5HEUzNnafTVdKSWC0iZvbOiYhXgLUqeN40YDown6ThwHzAM3222ZY0+IyIuB1YRNLilQRv\n",
       "M78hjS/xnFetznM3tYOtSXPIjc8dpBtUUiQkaUTJjRGkVsGAimJyPKkb7bPAqxFxXZ/NlgaeKrn9\n",
       "NOkfv6N4gJ1ZfRQD536OB841TSVThR8P3CbpYtIU4V8FjhrsSZKWB74DjASmApdIGhsRF/bdtM/t\n",
       "frtbSRpXcrMnInoqyN4yIrhOYgppNt0Tcucxa1O7Ay+S1m+xPiSNBkbXc5+DFomIOL9YnW5T0gf4\n",
       "l4sZYQezDnBrRLwMIOlyYCOgtEg8AyxTcvujfPCUVG+OcRW8Zqv7HvA3iUsjeDJ3GLN2UnQp/zHw\n",
       "2Yj+v0x2u+LLc0/vbUlH1LrPgdaTWKj47wjgOWACafT186WnnwYwBdhA0rySBGwO9C0uV5AmDkTS\n",
       "BqRTUh07L1QEU4CTgF/5IrZZ1U4lrRVxd+4g3WSglsRFwOeBf9D/KaCPD7TjiJgs6XxgIjCj2M94\n",
       "SXsVj58REVdL+pykh4E3SHOwdLpjgTuBscAFmbOYtQWJ7YFR4JHVzeYR1xlIrANcBYyK4MXcecxa\n",
       "mcSipJkLvhbBzbnztJOmjLiW9NdK7rPKFRORnUda7tTMBnY8cLkLRB5lTzdJmpc0tuHDfa5BLETq\n",
       "umq1OQK4R+KLEfwxdxizViSxBanTzKjcWbrVQNck9gIOAJYC7iq5/zXSBSSrQQRvSewBXChxo5dc\n",
       "NJtdMQ34mcD/i+C13Hm61aDXJCTtFxGnNClPuQwddU2ilMTpwPAI9sydxayVSJwELBLBzrmztKtm\n",
       "zQL7gqQFixc8XNLlkiqZlsMqcxCwpcRmuYOYtQqJjYCvAd/NnaXbVVIkDo+I1yRtDGwGnAP8qrGx\n",
       "ukcE00jLnJ4pMX/uPF3PczdlJzEPcDawfwQv587T7SopEu8X//0CMD4i/g+Ys3GRuk8EVwF/B36U\n",
       "O4tZCziUNBj30txBrLIi8YykM4EdgKskzVPh86w6BwBjJdbPHcQsF4nVSZ1m9vHUG62hkg/7rwHX\n",
       "AJ+NiFeBRYH/bWiqLhTBS6QJEc/2TLHWjYoZXs8GDo7gudx5LBl07iZgbuAG4OVivMQ7pKk2rP5+\n",
       "BzwK/CB3ELMMvgu8CpybO4jN0rC5m6x6EYTEPsAkicsiuDd3JrNmkFgJ+D6wnk8ztRbP3dSCJL4F\n",
       "7AFsWCxYZNaxJOYgna34fQQn5s7TSTx3U+caD7xOupht1um+BcwFZB20a/3z3E0tqDjt9C3g7xJ/\n",
       "jOCR3JnMGkFiGdJCQpu41dyaPHdTi4rgYYmjgfESm/k8rXWaYuGt04GTIz6wIJm1iErmbto/IrJO\n",
       "ad1t1yR6FV0CbwXOjOCs3HnM6kliDHAwsE4E7+bO04nq8dlZ0YVrSZ8CPgHM03tfRJxfywtXo1uL\n",
       "BIDEKOB6YI2I/tf/Nms3Eh8B7gG2ieDO3Hk6VVOKhKRxwCbAJ0mrqW0N3BwRX6nlhavRzUUCQOJI\n",
       "YE3giz7t1GC98zaNidFZc3Q4iYuApyM8MLeRmjUL7FeAzYHnImJXYHVgkVpe1Kr2U2A50uh3s7Ym\n",
       "sS2wDmnhLWtxlRSJtyLifeA9SQsDLwLLNDaWlYrgHWB34ESJxXLnMRsqiYWB04A9I3gzdx4bXCVF\n",
       "YqKkRUl99ycCk0gXU62JIridNAr+1KJXiFk7+gVwVQQ9uYNYZQbqAgtAROxd/PorSdcAC0XE5MbG\n",
       "sjIOA24D9iF9GzNrGxK7AhsC6+XOYpWrasR1RDwWEZM94jqPonm+HXCExAa585hVSmJN4FhgO69X\n",
       "3V484rrNRPCIxO7AJRJrR/Bi7kwdxb2a6k5iBHAZsG8E/8ydx6rjEddtKIIri8WJfivx2Qjey53J\n",
       "rD/F5H2/Af4QwcW581j1POK6TUkMA/4E/COCg3PnMeuPxA+BLYBNI5ieO0+3aehgOknrAk9HxHPF\n",
       "7Z2B7YHHgXER8UotL1xVSBeJfhXdYScC/xPB73PnMSslsRVwFrCuV5rLo9GD6c4krUKHpM8APwPO\n",
       "A6YVj1lmxZKnXwXOKBZtMWsJEiNJnxc7uUC0t4GKxBwlrYUdgDMi4rKIOAxYsZKdS1pZ0qSSn6mS\n",
       "9u+zzWKS/izpbkn3SdplSO+kSxXz3hwOXC4xf+48ZhLzkC5U/yyCv+XOY7UZ6HTTfcCaETFd0oPA\n",
       "tyLixuKx+yPik1W9kDQH8AywXkQ8VXL/OGDuiDhE0mLAg8DiEfFeyTY+3TSAYnDdOaT1yMd6fqca\n",
       "eO6mmkmcBSwI7Oi/xbwafbrpIuBGSVcAb0L6RiBpRdJi5dXaHHiktEAUniN1q6X478ulBcIGV/yP\n",
       "uA+wKvDtzHGsixXdszcCdneB6Axlu8BGxFGSrgeWAP4SETOKhwTsN4TX2hGY0M/944HrJT1L+vbh\n",
       "SeyGIIK3JLYHbpO4K8JTp1hzSawNHA18JoLXc+ex+hhwWo6IuK2f+/5V7YtImgvYBjion4d/ANwd\n",
       "EaMlLQ9cK2n1iHitzz7GldzsiYieanN0uggeldgNuLgYaPdC7kzWHSQ+BFwK7BPBlNx5upWk0cDo\n",
       "eu5z0Lmb6mRr4K6I+Hc/j20EHAUQEY9IegxYmdS1c6aIGNfokJ0ggqskziENtNvCA+2s0YoxOxcA\n",
       "l0Vwae483az48tzTe1tSzdOxVzILbD3sRLrG0Z8ppOsVSFqcVCAebVKuTnUkqfvyT3MHsa5wODAv\n",
       "eFBnJ6po+dKaXkCaH3gC+HjvKSRJewFExBlFj6ZzgWVJRevoiJjQZx/u3VSlovk/ETgwgsty57HO\n",
       "JPE50ripdSJ4Pncem13T1rjOzUViaCTWAa4GPh3Bg7nzWGeR+Djwd9LMrrfkzmMf1KzlS61NRTCR\n",
       "1DHgcokFcuexziExL2nA3FEuEJ3NLYkuIHE2MD9pioTW/we3llYM3jwbmAcP3mxpbklYpb5Nmkpl\n",
       "/8E2NKvAHqTV5b7lAtH53JLoEiXnj7eP4Obceaw9SawLXAVsHEHVY6asudySsIpF8BiwM2lFu0/k\n",
       "ztOyJqhn5vxNNhuJFYDLgb1cILqHi0QXieDPwIHAXzy1uFWjaIn+FTjSa5d0l2aNuLYWEcGFEnMD\n",
       "10n8dwSP5M5krU1iWVKBODaCs3LnseZykehCEZwjMRfwV4lNIngidyZrTRJLkQrEKRGcljuPNZ+L\n",
       "RJeK4FdFobi+KBRP585krUVicVKBOCuCX+TOY3m4SHSxCE4uaVGM9jKT1kviw6QCcVEEx+TOY/m4\n",
       "C6whcSgwFhgdwYu581heEiOA60ldXQ/zWIj2VY/PTrckjAiO6nMx++XcmSwPiYWBa4DrcIEw3JKw\n",
       "QjHVws9I07ZvFjGkJWqtjUksCPwFuBM4wAWi/XkWWKurolCcQFoIaosIpmWOZE0iMT/wJ+ABYG8X\n",
       "iM7gImF1VxSKU4E1gC29VnHnK2Z0/T/Sui97RDBjkKdYm3CRsIaQmAM4A1gB+HwEb2aOZA0iMQ/w\n",
       "B+AlYOcI3s8cyerIRcIapigUvwaWALaN4O28iZqkd96mMTE6a44mKLo/Xwa8BYzxeuidxxP8WcMU\n",
       "pxx2BV4GLit6P1mHkJgT+C3wPmlNCBcI65eLhJVVnHr4JvA28Lvig8XanMRw4DfA3MAOEUzPHMla\n",
       "mIuEDaj4ANmJ9LcyofiAsTYlMQw4BxhBWlvkncyRrMW5SNigIngX+CqwIHBh0RvG2kxxyvBsYBng\n",
       "S11znclq4iJhFSm+cX6ZdA77VonlM0eyKkh8DPgbqdBv4x5rVin3brKqFOMo9gGOAPaM4I+ZI9kg\n",
       "JLYi9VT7OXCCB8p1D3eBtWwk1gcuJvWQOdS9Y1pPcf3hh8DuwE4R/C1zJGsyFwnLSmIx4EJgLtKH\n",
       "0POZI1mh5N9mbmBH/9t0J4+TsKwieAn4HHAjMFHi05kjGTNbeXcBk4DNXSCsFm5JWF2UnPc+Djje\n",
       "572br7gJYKKZAAALaElEQVRetC/pFNMeEVyROZJl5tNN1lKKHjSXAE8Bu0UwNXOkriGxADAeWJU0\n",
       "/uGRzJGsBbT06SZJK0uaVPIzVdL+/Ww3unj8PqmYN8faUgRPAJ8GngfulFgtc6TqTVDPzPmb2oTE\n",
       "qsAdwJvAhi4QVk8NGz0bEQ8CawJImgN4Bvh96TaSFgFOA7aMiKclLdaoPNYcxXiKfSXGktbOPjCC\n",
       "83Ln6lQSOwKnAAdFcE7uPNZ5mjXFwubAIxHxVJ/7xwCXRcTTABHxUpPyWINFcKHE3aTJATcirXTm\n",
       "Eb51UszgehzweeCzEUzKHMk6VLN6N+0ITOjn/hWBEZJukDRR0jealMeaIIL7gXWBRYFbJD6eOVJH\n",
       "kFiG1KPsY8DaLhDWSA0vEpLmArYhXdDsa05gLVI3yi2BwyWt2OhM1jwRvAbsQJp19HaJbTNHamsS\n",
       "W5LWoP4D8GWvRW6N1ozTTVsDd0XEv/t57CngpYh4C3hL0k3A6sBDfTeUNK7kZk9E9DQgqzVA0R32\n",
       "RIk7gfMk9gcOj+C2zNHahsQ6wI9JvZd2iuCGzJGsBUkaDYyu6z4b3QVW0m+BP0XEBy5eSlqFtJ7y\n",
       "lqSRobcDO0TEA322cxfYDlGsSbEzcDhwH/DDCO7Km6p1FT3EfgSsAxwFnF3Myms2qJbuAgsgaX7S\n",
       "RevLS+7bS9JeABExBfgzcA+pQIzvWyCss0QwPYKzgJWAq4ErJC6XGJU5WkuRWFXid8A1QA+wYgSn\n",
       "u0BYs3kwnWVVrE2xN/B94AZgXAQP5k2VTzEF+xHAVsDxwKkRvJE3lbWrlm9JmA0mgrciOAFYgdSi\n",
       "vFni1xLLZY7WVBLLSowntagfBlaI4BgXCMvNRcJaQgSvR3A0qVg8DtwhcUbR3bNjSSwlcSppMr4X\n",
       "gZUi+FEE0zJHMwNcJKzFRDA1gnHAysArwN0SJ0ssmTdZfUl8ROJ44F7gbWCVCA6N4JXM0cxm4yJh\n",
       "LSmClyM4BPgE8B5wn8TPJT7c0Bdu8NxNEiMkfgr8k7QOx6gIDoygvy7iZtm5SFhLi+CFCL4LrAbM\n",
       "Bzwscb3EYRIbFV1qW5bEcIkNJA6RuJZ0Km0xYM0I9ovg2bwJzQbm3k3WViQWAjYGNi1+VgBuAa4n\n",
       "9Y6aFMH7Q36B3lbEmBg9xHxzkAaEbgr8N2lW3MeKbNcDN3kKdWuWenx2NmuCP7O6KC7oXl38IDEC\n",
       "2IT0oXwesJTEjcz6UL4/ghmNylMs9LMqs4rWJqQL0DeQFmHapVjBz6wtuUhYWysu9P6++EFiCdK0\n",
       "BJsC+wELSdzArKLxUC2r5hVFYXlSK6G3tfBmsf/LgG/7FJJ1krY53cQ4bux7fxzR/ykBHdn/hUdv\n",
       "333bSyzL/yxzDW8vsijvLLQIaA40I7Uszr71A2s/33LERktM//B9c42e9KkPjmxO289LupDee3rr\n",
       "Bsbp3Ebl9/bevpbtfbrJbBARPKkjn36BhZ9+gQDem2cuYo7eDhs79t3+v468FX4w/29ZfHJ/u9sR\n",
       "mA48W9oa0ZGNSG7WGtqmJeEL12Zm1fG0HGZm1lAuEmZmVpaLhJmZleUiYWZmZblImJVq8NxNZu3G\n",
       "RcLMzMpykTAzs7JcJMzMrCwXCTMzK8tFwszMyvK0HGZmHcrTcpiZWUO5SJiZWVkuEmZmVpaLhJmZ\n",
       "leUiYWZmZblImJXy3E1ms2lYkZC0sqRJJT9TJe1fZtt1Jb0nabtG5TEzs+o1rEhExIMRsWZErAms\n",
       "DbwJ/L7vdpKGAccAfwY8FmIQkkbnztAqfCxm8bGYxceivpp1umlz4JGIeKqfx/YDLgX+3aQs7W50\n",
       "7gAtZHTuAC1kdO4ALWR07gCdpFlFYkdgQt87JS0NfBE4vbir9Yd/m5l1kYYXCUlzAdsAl/Tz8InA\n",
       "wZHmBhE+3WRm1lIaPneTpC8Ce0fEVv089iizCsNipOsWe0bEFX22cwvDzGwIap27aXi9ggxgJ+Ci\n",
       "/h6IiOV6f5d0LnBl3wJRbOcWhplZBg093SRpftJF68tL7ttL0l6NfF0zM6uPtpgq3MzM8sg64lrS\n",
       "MpJukHS/pPt6B9tJ+rmkf0qaLOlySQtX+tx2VMtxKNnHsGLQ4pXNS15/tR4LSYtIurTY9gFJGzT3\n",
       "HdRPHY7FIcVz75U0QdLczX0H9TPAsfhxcRzulvRXScuUef5WkqZIekjSQc1NX1+1HIshfW5GRLYf\n",
       "YAlgjeL3BYAHgVWBLYA5ivt/Bvys0ufmfD85jkPJPr4LXAhckfv95DwWwHnAbsXvw4GFc7+nHMcC\n",
       "GAk8Csxd3P4dsHPu99SAY7FgyTb7AWf189xhwMPFMZkTuLtdPyvqcCyq/tzM2pKIiOcj4u7i99eB\n",
       "fwJLRcS1ETGj2Ox24KOVPrc5yeurluMAIOmjwOeAs2jzbsS1HIviG/WnI+Kc4vnvRcTUJkWvuxr/\n",
       "LqYB04H5JA0H5gOeaULshhjgWLxWstkCwEv9PH094OGIeDwipgO/JY3Paku1HIuhfG62zAR/kkYC\n",
       "a5L+6EvtBlw9xOe2nSEeh18A/wvMKPN4WxrCsfg48G9J50r6h6TxkuZrbMrmqPZYRMQrwPHAk8Cz\n",
       "wKsRcV1jUzZH32Mh6ShJTwI7k1pWfS0NlM728HRxX9sbwrEo+9xyWqJISFqANDXHAUV1673/UODd\n",
       "iPjAaO3BntuOhnIcJH0BeDEiJtHmrYhSQ/ybGA6sBfwyItYC3gAObkbeRhri38XywHdIp1iWAhaQ\n",
       "NLY5iRunv2MREYdGxLLAr0lfmPrqyN45QzwWZZ9bTvYiIWlO4DLggoj4Q8n9u5BOoZT9wy733HZU\n",
       "w3HYCNhW0mOk8SibSjq/wXEbqoZj8TTwdETcWdy+lFQ02lYNx2Id4NaIeDki3iN1Q9+owXEbqoL/\n",
       "3ycA6/Zz/zNA6UXcZUh/K22rhmNR/edm5gswAs4HftHn/q2A+4HFqn1uO/7Uchz6bL8JaUBi9veU\n",
       "61gANwErFb+PA47J/Z5yHAtgdeA+YN5iP+cB++Z+Tw04FiuW/L4f8Jt+njsceITUqpqL9r9wXcux\n",
       "qPpzM/eb3Zh0Hv1uYFLxszXwEPBEyX2/LLZfCrhqgOdulfsfsNnHoc9+NqH9ezfVdCyKD8c7gcmk\n",
       "b8/t3Lup1mPx/aKY3FsUiTlzv6cGHItLi/d3N+nb8UfKHIutST15HgYOyf1+ch2LoXxuejCdmZmV\n",
       "lf2ahJmZtS4XCTMzK8tFwszMynKRMDOzslwkzMysLBcJMzMry0XCrAqSeiStXcPzR0i6VtK/JP1F\n",
       "0iL1zGdWby4SZtUJapsL6GDg2ohYCfgrHTC3lHU2FwnrapJGFovRXFAsUnSJpHklbVbMJHuPpLMl\n",
       "zTX707SrpF+U3LGnpBP67Pug4vl3S/ppcfe2pNHPFP/9UkPfoFmNXCTMYCXgtIj4BGkdhu8B5wJf\n",
       "i4jVSHP/7F2yfQAXA9tIGlbctwtwdu8GkrYmFYT1ImIN4NjiocUj4oXi9xeAxRvyjszqxEXCDJ6K\n",
       "iNuK3y8ANgUejYiHi/vOAz5T+oSIeAO4nlQoViHNi3R/ySabAedExNvF9q/2fdFIc+J4XhxraS4S\n",
       "ZrN/UAt4ldnX5ii3TsdZwK6kVsQ5/Tze3/NekLQEgKQlgRerDWvWTC4SZrCspA2K38cAE4GRxcI9\n",
       "AN8Aekq2F0BE3EFaOnQMaS2PUtcCu0qaF0DSosX9V5BWDaP4b1uvg2Kdz0XCLE0hva+kB4CFgRNI\n",
       "LYRLJN0DvAf8qmT70pbHxcDNETFV0jqSxgNExDWkgjBR0iTSdQ5IS0puIelfpNNaAy4xaZabpwq3\n",
       "rlas83tlRIwa4vOvBE6IiBvqmcusVbglYTaEi8eSFpH0IPCmC4R1MrckzMysLLckzMysLBcJMzMr\n",
       "y0XCzMzKcpEwM7OyXCTMzKwsFwkzMyvr/wO1CDMUOGq07QAAAABJRU5ErkJggg==\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sherpa.int_proj(poly.c0)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}
