Autoregressive Moving Average (ARMA): Artificial data
=====================================================


.. _tsa_arma_1_notebook:

`Link to Notebook GitHub <https://github.com/statsmodels/statsmodels/blob/master/examples/notebooks/tsa_arma_1.ipynb>`_

.. raw:: html

   
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[1]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="kn">from</span> <span class="nn">__future__</span> <span class="kn">import</span> <span class="n">print_function</span>
   <span class="kn">import</span> <span class="nn">numpy</span> <span class="kn">as</span> <span class="nn">np</span>
   <span class="kn">import</span> <span class="nn">statsmodels.api</span> <span class="kn">as</span> <span class="nn">sm</span>
   <span class="kn">import</span> <span class="nn">pandas</span> <span class="kn">as</span> <span class="nn">pd</span>
   <span class="kn">from</span> <span class="nn">statsmodels.tsa.arima_process</span> <span class="kn">import</span> <span class="n">arma_generate_sample</span>
   <span class="n">np</span><span class="o">.</span><span class="n">random</span><span class="o">.</span><span class="n">seed</span><span class="p">(</span><span class="mi">12345</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>Generate some data from an ARMA process:</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[2]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">arparams</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="o">.</span><span class="mi">75</span><span class="p">,</span> <span class="o">-.</span><span class="mi">25</span><span class="p">])</span>
   <span class="n">maparams</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">array</span><span class="p">([</span><span class="o">.</span><span class="mi">65</span><span class="p">,</span> <span class="o">.</span><span class="mi">35</span><span class="p">])</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p>The conventions of the arma_generate function require that we specify a 1 for the zero-lag of the AR and MA parameters and that the AR parameters be negated.</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[3]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">arparams</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">r_</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="o">-</span><span class="n">arparams</span><span class="p">]</span>
   <span class="n">maparam</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">r_</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="n">maparams</span><span class="p">]</span>
   <span class="n">nobs</span> <span class="o">=</span> <span class="mi">250</span>
   <span class="n">y</span> <span class="o">=</span> <span class="n">arma_generate_sample</span><span class="p">(</span><span class="n">arparams</span><span class="p">,</span> <span class="n">maparams</span><span class="p">,</span> <span class="n">nobs</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing text_cell rendered">
   <div class="prompt input_prompt">
   </div>
   <div class="inner_cell">
   <div class="text_cell_render border-box-sizing rendered_html">
   <p> Now, optionally, we can add some dates information. For this example, we&#39;ll use a pandas time series.</p>
   </div>
   </div>
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[4]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">dates</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">tsa</span><span class="o">.</span><span class="n">datetools</span><span class="o">.</span><span class="n">dates_from_range</span><span class="p">(</span><span class="s">&#39;1980m1&#39;</span><span class="p">,</span> <span class="n">length</span><span class="o">=</span><span class="n">nobs</span><span class="p">)</span>
   <span class="n">y</span> <span class="o">=</span> <span class="n">pd</span><span class="o">.</span><span class="n">TimeSeries</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">index</span><span class="o">=</span><span class="n">dates</span><span class="p">)</span>
   <span class="n">arma_mod</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">tsa</span><span class="o">.</span><span class="n">ARMA</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">order</span><span class="o">=</span><span class="p">(</span><span class="mi">2</span><span class="p">,</span><span class="mi">2</span><span class="p">))</span>
   <span class="n">arma_res</span> <span class="o">=</span> <span class="n">arma_mod</span><span class="o">.</span><span class="n">fit</span><span class="p">(</span><span class="n">trend</span><span class="o">=</span><span class="s">&#39;nc&#39;</span><span class="p">,</span> <span class="n">disp</span><span class="o">=-</span><span class="mi">1</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[5]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="k">print</span><span class="p">(</span><span class="n">arma_res</span><span class="o">.</span><span class="n">summary</span><span class="p">())</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stdout output_text">
   <pre>
                                 ARMA Model Results                              
   ==============================================================================
   Dep. Variable:                      y   No. Observations:                  250
   Model:                     ARMA(2, 2)   Log Likelihood                -245.887
   Method:                       css-mle   S.D. of innovations              0.645
   Date:                Sun, 22 Mar 2015   AIC                            501.773
   Time:                        14:58:18   BIC                            519.381
   Sample:                    01-31-1980   HQIC                           508.860
                            - 10-31-2000                                         
   ==============================================================================
                    coef    std err          z      P&gt;|z|      [95.0% Conf. Int.]
   ------------------------------------------------------------------------------
   ar.L1.y        0.8411      0.403      2.089      0.038         0.052     1.630
   ar.L2.y       -0.2693      0.247     -1.092      0.276        -0.753     0.214
   ma.L1.y        0.5352      0.412      1.299      0.195        -0.273     1.343
   ma.L2.y        0.0157      0.306      0.051      0.959        -0.585     0.616
                                       Roots                                    
   =============================================================================
                    Real           Imaginary           Modulus         Frequency
   -----------------------------------------------------------------------------
   AR.1            1.5618           -1.1289j            1.9271           -0.0996
   AR.2            1.5618           +1.1289j            1.9271            0.0996
   MA.1           -1.9835           +0.0000j            1.9835            0.5000
   MA.2          -32.1811           +0.0000j           32.1811            0.5000
   -----------------------------------------------------------------------------
   
   </pre>
   </div>
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[6]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="n">y</span><span class="o">.</span><span class="n">tail</span><span class="p">()</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt output_prompt">
       Out[6]:</div>
   
   
   <div class="output_text output_subarea output_pyout">
   <pre>
   2000-06-30    0.050999
   2000-07-31   -0.206404
   2000-08-31   -0.170874
   2000-09-30    0.257949
   2000-10-31    0.245237
   dtype: float64
   </pre>
   </div>
   
   </div>
   
   </div>
   </div>
   
   </div>
   <div class="cell border-box-sizing code_cell rendered">
   <div class="input">
   <div class="prompt input_prompt">
   In&nbsp;[7]:
   </div>
   <div class="inner_cell">
       <div class="input_area">
   <div class="highlight"><pre><span class="kn">import</span> <span class="nn">matplotlib.pyplot</span> <span class="kn">as</span> <span class="nn">plt</span>
   <span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">(</span><span class="n">figsize</span><span class="o">=</span><span class="p">(</span><span class="mi">10</span><span class="p">,</span><span class="mi">8</span><span class="p">))</span>
   <span class="n">fig</span> <span class="o">=</span> <span class="n">arma_res</span><span class="o">.</span><span class="n">plot_predict</span><span class="p">(</span><span class="n">start</span><span class="o">=</span><span class="s">&#39;1999m6&#39;</span><span class="p">,</span> <span class="n">end</span><span class="o">=</span><span class="s">&#39;2001m5&#39;</span><span class="p">,</span> <span class="n">ax</span><span class="o">=</span><span class="n">ax</span><span class="p">)</span>
   <span class="n">legend</span> <span class="o">=</span> <span class="n">ax</span><span class="o">.</span><span class="n">legend</span><span class="p">(</span><span class="n">loc</span><span class="o">=</span><span class="s">&#39;upper left&#39;</span><span class="p">)</span>
   </pre></div>
   
   </div>
   </div>
   </div>
   
   <div class="output_wrapper">
   <div class="output">
   
   
   <div class="output_area"><div class="prompt"></div>
   
   
   <div class="output_png output_subarea ">
   <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAk0AAAHrCAYAAADSV6k2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd4VEXbBvD77KYnpNAJLfSmFCH0EqRIlSpIEeEFFKUp
   gqhgR4UPxI4ILy8qKiChSq8h0kMvgdACgUBIAmmkJ3u+P4YNCaRsObt7Nrl/15UrbDlzhgDhycw9
   M5IsyyAiIiKiwmls3QEiIiIie8CiiYiIiMgALJqIiIiIDMCiiYiIiMgALJqIiIiIDMCiiYiIiMgA
   ZhVNkiS5SJJ0VJKk05IkhUqS9JVSHSMiIiJSE8ncfZokSXKTZTlFkiQHAAcATJdl+YAivSMiIiJS
   CbOn52RZTnn0SycAWgAPzG2TiIiISG3MLpokSdJIknQawD0A+2RZDjW/W0RERETqosRIk06W5aYA
   qgDoKElSgNm9IiIiIlIZB6UakmU5QZKkLQBaAAjSPy9JEg+3IyIiIrshy7KU3/Pmrp4rK0mS96Nf
   uwLoBuBUPje32MfHH39s1+0Xh99DcfgaNWnSxK77Xxz+DNi+7e/B9m1/D34vsv09CmPuSFMlAL9J
   kqSBKMBWyLK8x8w2jRIQEGDX7VvjHvbevjXuUbFiRYu2zz8Dtm8P92D7tr8Hvxep4x4FMXvLgSJv
   IEmype9BZK5PPvkEn3zyia27QUQlHL8X2Z4kSZAtMT1HVFzY8icXIiI9fi9SNxZNRERERAZQbPWc
   sSQp35EvshBOkRIREZnHZpmmR3OGFr03CfxaExERGYaZJiIiIiIzsWgiAhAUFGTrLhAR8XuRyrFo
   IiIiIjIAi6Z8hIWFoWnTpvD09MSPP/5o6+6QFXCZLxGpAb8XqRuLpnz83//9H7p06YLExERMmjTJ
   1t0plJ+fH/bu3WvrbhARERV7LJrycfPmTTRs2NDo67Kzsy3Qm8JxZZwymCMgIjXg9yJ1Y9H0hOef
   fx5BQUGYNGkSPD09cfbsWYwaNQrly5eHn58fvvjii5wi5ddff0W7du0wbdo0lC1bFp9++ikyMjIw
   ffp0VK9eHRUrVsQbb7yBtLS0nPY3btyIpk2bwsvLC7Vr18aOHTsAAMuXL0fDhg3h6emJWrVqYcmS
   JTnXxMbGok+fPvDx8UGZMmXQsWNHyLKMV155BREREejbty9KlSqFBQsWWPeLRUREVJJY+jRicYun
   FfS8GgQEBMjLli2TZVmWX3nlFbl///7yw4cP5Rs3bsh169bNeW358uWyg4OD/OOPP8rZ2dlyamqq
   /NZbb8n9+vWT4+Li5KSkJLlv377y+++/L8uyLB89elT28vKSd+/eLcuyLEdGRsqXLl2SZVmWt2zZ
   Il+/fl2WZVnev3+/7ObmJp86dUqWZVl+77335AkTJshZWVlyVlaWfODAgZy++vn5yXv27Cn096Pm
   rzUREZGaPPo/M9+aRrUjTZKkzIc5srOzsXr1anz11Vdwd3dH9erV8c4772DFihU57/H19cXEiROh
   0Wjg7OyMpUuXYuHChfD29oaHhwfef/99rFq1CgCwbNkyjB07Fl26dMm5tl69egCAXr16oUaNGgCA
   jh07onv37ggODgYAODk54e7du7hx4wa0Wi3atWtn3m+MiIiIjKbaokmWlfkwR2xsLDIzM1G9evWc
   56pVq4bIyMicx1WrVs35dUxMDFJSUtC8eXP4+PjAx8cHPXv2RGxsLADg9u3bqFWrVr732rZtG1q3
   bo0yZcrAx8cHW7duxf379wEAM2bMQO3atdG9e3fUqlUL8+bNM+83Rk9hjoCI1IDfi9RNtUWTGpQt
   WxaOjo64ceNGznMRERGoUqVKzuPcZ+iVLVsWrq6uCA0NRVxcHOLi4hAfH4/ExEQAosC6evXqU/dJ
   T0/HoEGD8O677yI6OhpxcXHo1atXTnbKw8MDCxYswLVr17Bp0yYsXLgQ+/bte+r+REREZDksmgqh
   1WoxZMgQzJo1Cw8fPsTNmzfxzTffYOTIkfm+X6PRYPz48XjrrbcQExMDAIiMjMTOnTsBAGPHjsXy
   5cuxd+9e6HQ6REZGIiwsDBkZGcjIyEDZsmWh0Wiwbdu2nGsAYPPmzbh69SpkWYanpye0Wi00GvFH
   V6FCBVy7ds3CX4nij3ujEJEa8HuRurFoKsIPP/wAd3d31KxZEx06dMCIESMwZswYAGKU58mRnnnz
   5qF27dpo3bo1vLy80K1bN1y+fBkA4O/vj+XLl+Ptt9+Gt7c3AgICEBERgVKlSuH777/HkCFDULp0
   aaxcuRL9+vXLafPq1avo1q0bSpUqhbZt22LixIno1KkTAOD999/HnDlz4OPjg4ULF1rpq0JERFTy
   SLK5wZ+ibiBJcn734P5C1sOvddGCgoL4Ex4R2Ry/F9neo/8z882+cKSJiIiIyAAcaSoB+LUmIiIy
   DEeaiIiIiMzEookI3BuFiNSB34vUjUUTERERkQGYaSoB+LUmIiIyDDNNRERERGZi0UQE5giISB34
   vUjdWDQRERERGYCZphKAX2siIiLDMNNEREREZCYWTU+YP38+Bg8enOe5KVOm4K233rJRj8gamCMg
   IjXg9yJ1Y9H0hFdeeQXbt29HQkICACArKwurV6/Gq6++auOeERERkS2pNtMkfZrvdKLR5I+N//31
   7NkTgwYNwrhx47B582a89957OH/+vCL9sQVmmoiIiAxTWKZJtUWTLa1atQqLFy9GUFAQXn75ZTRr
   1gwzZ860dbdMpuavNRERkZqwaDJSamoqKleujODgYLRp0wYXL15ElSpVbN0tk6n5a60WQUFBCAgI
   sHU3iKiEK0nfi2RZhk6ny/mc30d+r8mynPN/Wn6/1j/O/d7cv879+Mn3ybKM5557rsCiycEaXxh7
   4+rqikGDBmH48OFo1aqVXRdMRERE5pBlGdnZ2cjKysr3IzMzExkZGUhPT0dGRgays7NzPnQ6XZ7P
   +l/rH0uSBEnKP46jf/7JH/plWS7ytdzvyd1+Ub9OTU0t9GvBkaYCHDhwAB07dsTy5cvtPgSu9q81
   ERFZXnZ2NlJTU/MtfPQFj/5zRkZGTjGUlZWVpwjJTT86o9Focj70hVBRHwAKLJhsJT4+HgMGDOBI
   k7GqV6+eM+JERERkT2RZRlpaGhITExEXF4fo6GjEx8fnOzIDiOIld+Gj/3BxcckphIhFU750Oh2+
   /vprDBs2DB4eHrbuDllBScoREJF6mfq9KDs7G0lJSUhISEBsbCxiYmKQlpaWU+w4OzvDw8MDGg13
   GjIHi6YnJCcno0KFCqhRowa2b99u6+4QERHlkXsU6cGDBzmjSPrXHB0d4ezsDC8vL44QKYyZphKA
   X2siIvuVexQpOjoasbGxSE9PzzOK5OzszAJJAcw0ERER2Qn9KFJCQkKeLJL+Nf0okouLC4skG2DR
   RARmmojItmRZxpUrV7B+/XrUr18/Z7Was7MzPD09WSCpBIsmIiIiG8rOzsbZs2cRHh4OBwcHeHl5
   2bpLVACbFk2snEktOMpERLaQmZmJ48ePIyoqCj4+PihdurStu0SFsFnRxGAyERGVZKmpqTh69CgS
   EhLg7e3NgQQ7wA0biCAyTURE1pKUlITg4GAkJSXl2Rrg9OnTNu4ZFYaZJiIiIit68OABDh8+DAAo
   VaqUjXtDxrDZPk1EREQlzZ07dxASEpKzbQCpC/dpIiIisjFZlhEeHo4zZ87Aw8MDjo6Otu4SmYCZ
   JiIw00REliPLMi5evIjTp0/D09Oz0IKJmSZ140gTERGRhWRnZ+P06dOIiIiAt7c3D8y1c8w0ERER
   WUBGRgZCQkIQHR3NLQXsBDNNREREVpaSkoIjR47g4cOHLJiKEY4TEoGZJiJSTmJiIoKDg5GSkmL0
   uXHMNKkbR5qIiIgUEhsbi8OHD0Or1cLDw8PW3SGFMdNERESkgMjISISEhMDV1RXOzs627g6ZgJkm
   IiIiC5JlGdeuXcO5c+e4B1Mxx0wTEZhpIiLT6HQ6nD9/HmfPni1yDyZDMNOkbhxpIiIiMkFWVhZO
   nTqFW7ducQ+mEoKZJiIiIiOlp6fj2LFjuH//Pry8vLilQDHBTBMREZGCkpOTcfjwYaSkpLBgKmE4
   lkgEZpqIyDDx8fEIDg5GWlqa0XswGYKZJnXjSBMREZEBoqOjceTIETg6OnIPphKKRRMRgICAAFt3
   gYhULCsrCyEhIXB2drboHkxNmza1WNtkPk7PERERFeHevXvIzMzkppUlHIsmIjDTRESFu3btmlUK
   Jmaa1M2sokmSpKqSJO2TJOmCJEnnJUmaolTHiIiI1CApKQkPHjyAi4uLrbtCNmZupikTwNuyLJ+W
   JMkDwAlJknbJsnxRgb4RWQ0zTURUkIiICEiSZJWtBZhpUjezRppkWY6SZfn0o18/BHARgK8SHSMi
   IrK1rKwshIeHw93d3dZdIRVQLNMkSZIfgGYAjirVJpG1MNNERPm5d+8esrKy4OBgncXmzDSpmyJ/
   Cx5NzQUCmPpoxCmP0aNHw8/PDwDg7e2Npk2b5kyH6P+z4mM+tuVjPbX0h4/5mI9t/3jfvn04e/Ys
   6tatC+BxQaOfQrPE46tXr1q0fT5++rH+11FRUcjIyEBhzD57TpIkRwCbAWyTZfnbfF7n2XNERGR3
   EhMTsXfvXh6VUoIUdfacuavnJADLAITmVzARERHZq1u3blktAE72wdxMUzsAIwF0liTp1KOPHgr0
   i8iq9MPyRETA4wC4tY9LYaZJ3czKNMmyfADcIJOIiIoZfQBcq9XauiukIix4iMB9mojoMVmWrbYD
   +JO4T5O6sWgiIiLKhTuAU0FYNBGBmSYiesyWAXBmmtSNRRMREdEjtgqAk31g0UQEZpqISIiKirJp
   AJyZJnVj0URERATbBsDJPrBoIgIzTUQkAuBxcXE2DYAz06RuLJqIiIgAREREcAdwKhSLJiIw00RU
   0qklAM5Mk7qxaCIiohIvKioK2dnZ3AGcCsWiiQjMNBGVZGoKgDPTpG4smoiIqERTQwCc7AOLJiIw
   00RUkqkpAM5Mk7qxaCIiohJLLQFwsg8smojATBNRSaW2ADgzTerGoomIiEokNQXAyT6waCICM01E
   JVFiYqLqAuDMNKkbiyYiIiqRIiIioNFoVBEAJ/vAookIzDQRlTRZWVm4ceMG3N3dbd2VPJhpUjcW
   TUREVOKoLQBO9oFFExGYaSIqSWRZxtWrV1UZAGemSd1YNBERUYmSmJiI+Ph4VQXAyT6waCICM01E
   JYmaA+DMNKkbiyYiIioxMjMzVRkAJ/vAookIzDQRlRRqD4Az06RuLJqIiKhE4A7gZC4WTURgpomo
   JLCHADgzTerGoomIiEqEmzdvqjYATvaBRRMRmGkiKu4yMzNx8+ZN1QfAmWlSNxZNRERU7Kk9AE72
   gUUTEZhpIirO9DuAqznLpMdMk7qxaCIiomItMTERCQkJXDVHZmPRRARmmoiKM3sKgDPTpG4smoiI
   qNiylwA42QcWTURgpomouLK3ADgzTerGoomIiIolewqAk31g0UQEZpqIiqOEhAS7C4Az06RuLJqI
   iKhYioiIsJsAONkHFk1EYKaJqLix1wA4M03qxqKJiIiKnbt379pVAJzsA4smIjDTRFSc2HMAnJkm
   dWPRRERExUpCQgISExPtKgBO9oFFExGYaSIqTuxpB/AnMdOkbiyaiIio2MjMzERERITdBcDJPrBo
   IgIzTUTFhb0HwJlpUjcWTUREVCzYcwCc7AOLJiIw00RUHBSHADgzTerGoomIiIoFew6Ak31g0UQE
   ZpqI7F12dnaxCIAz06RuLJqIiMjuZWVlQafT2W0AnOwDiyYiMNNEROrATJO6sWgiIiIiMgCLJiIw
   00RE6sBMk7qxaCIiIiIyAIsmIjDTRETqwEyTurFoIiIiIjIAiyYiMNNEROrATJO6sWgiIiIiMgCL
   JiIw00RE6sBMk7qxaCIiIiIyAIsmIjDTRETqwEyTurFoIiIiIjIAiyYiMNNEROrATJO6sWgiIiIi
   MgCLJiIw00RE6sBMk7qZXTRJkvQ/SZLuSZJ0TokOEREREamREiNNywH0UKAdIpthpomI1ICZJnUz
   u2iSZflfAHEK9IWIiIhItZhpIgIzTUSkDsw0qRuLJiIiIiIDOFjjJqNHj4afnx8AwNvbG02bNs35
   yV6fJeFjPrblY/1zaukPH/MxHxv/+Ny5c/Dw8MgZrdHng+zp8dWrVzF48GDV9KckPNb/OioqChkZ
   GSiMJMtyoW8whCRJfgD+kWX52Xxek5W4B5ElBQUF5XzjJSL7k56eju3bt8PLy8vWXTHL6dOnOUVn
   Q/Hx8RgwYABkWZbye12JLQdWAjgEoK4kSbckSRpjbptE1saCiYjUgAWTupk9PSfL8jAlOkJERESk
   ZgyCE4H7NBGROnCfJnVj0URERERkABZNRGCmiYjUgZkmdWPRRERERGQAFk1EYKaJiNSBmSZ1Y9FE
   REREZAAWTURgpomI1IGZJnVj0URERERkABZNRGCmiYjUgZkmdWPRRERERGQAFk1EYKapuDt/Hti5
   09a9ICoaM03qZvbZc0REavfZZ8CdO0D37rbuCRHZM440EYGZpuIsPh7YsQM4fRpITrZ1b4gKx0yT
   urFoIqJibc0aoFs3oGlT4NAhW/eGiOwZiyYiMNNUnP3+OzBqFNC5M8ABRVI7ZprUjUUTERVb168D
   ly4BPXoAAQHAvn227hER2TMWTURgpqm4+uMP4OWXAScnoE0b4OxZ4OFDW/eKqGDMNKkbiyYiKpZk
   +fHUHAC4uQHPPQccPGjbfhGR/WLRRARmmoqjI0cABwegRYvHzzHXRGrHTJO6sWgiomJJP8okSY+f
   Y66JiMzBookIzDQVN+npwN9/AyNH5n2+TRuxO3hSkm36RVQUZprUjUUTERUpM9PWPTDOli1A48ZA
   tWp5n3dxEdN1Bw7Ypl9EZN9YNBGBmabCZGYC1asD587ZuieGyx0AfxJzTaRmzDSpG4smIirU8ePA
   3bvAt9/auieGiY0VRdGgQfm/zlwTEZmKRRMRmGkqTFCQyAatWwdER9u6N0VbvRro1Qvw9Mz/9Vat
   gNBQIDHRuv0iMgQzTerGoomIChUUBAweDAwZAvz8s617U7TCpuYAkWtq2RL491/r9YmIigcWTURg
   pqkgmZnA4cNAhw7AW2+JoiktTfn7bArbhPC4cLPbCQsDIiKArl0Lfx9zTaRWzDSpG4smIirQ8eNA
   rVpA6dJAgwZAs2bAX38pe4+94XsxNHAoRm8cDVmWzWprxQpg+HCxqWVhmGsiIlOwaCICM00FCQoS
   BYbetGkiEG5mbZPjRvwNDF87HP8M+wfJGcn47cxvJrel04mz5l55pej3tmwpRqXi402+HZFFMNOk
   biyaiKhATxZNXbuKgmnPHvPbTslMwYDVA/Be+/fQtWZX/NLnF8zcPRP3U+6b1N6//4rwd5MmRb/X
   2VkEwplrIiJjsGgiAjNN+cmdZ9KTJJFtWrjQvLZlWca4TePwTPlnMLXVVABAc9/mGNpoKGbunmlS
   m/kdm1KYgADmmkh9mGlSN7sumhISgB07bN0LouIpJORxnim3ESOAEyeAixdNb3vh4YUIux+GJX2W
   QMpV5cx5fg62X92OAxHGbdmdmgqsXy/yTIbq3Jm5JiIyjl0XTdu2AT17Aps22bonZO+YaXrak1Nz
   ei4uwBtvAN99Z1q7u67twoLDC7B+6Hq4Orrmec3T2RPfvPANJmyegIzsDIPb3LgR8PcHfH0N74e/
   P3DlChAXZ/g1RJbGTJO62XXRdP060K0bMHYscPKkrXtDVLwEBYnRmPy88YbYRDI21rg2w+PC8cr6
   V7By0EpU86qW73sGNxyMql5V8c3hbwxud8UKwwLguTk5iQN8mWuyjrQ0sbM8kT2z+6Jp4EBg8WKg
   Xz/g9m1b94jsFTNNeWVkPJ1nyq1CBWDAAOCXXwxvMzkjGf1X98esDrMQ4BdQ4PskScJPvX7C/EPz
   cSP+RpHtRkUBBw+K/hiLWw9Yz3vvAa1bi1gFFYyZJnWz+6KpZk1xxtTkyUDfvsDDh7buFZH9O34c
   qFMH8PEp+D1vvw389JMosIoiyzL+s+k/aFaxGSa1nFTk+2v61MS0NtMwaeukIvduWrkS6N8fcHcv
   uh9P4iaX1nHrlgjqt2sHTJxo694Qmc6ui6bwcFE0AcCMGUDz5sCwYUB2tm37RfaHmaa8Csoz5fbs
   s0CjRmKarijzD81HeFw4FvdZnCf4XZjpbafjetx1rL+0vtD3mTI1p9eiBXDtGvDggWnXk2G++AIY
   Px74739FQb5ypa17pF7MNKmb3RZNmZnAnTtAtUexCEkSRzykpADvvGPbvhHZO0OKJkCMNi1cWPhm
   lzuu7sC3R77FuqHr4OLgYnAfnLRO+Ln3z5i6fSqS0pPyfc/580BMjGF9zY+jI9C2LRAcbNr1VLTw
   cGDNGuDddwE3N+DPP4GpU8VxN0T2xm6LpogIsVLG0fHxc46OQGCg2Ibgp59s1zeyP8w0PVZUnim3
   Hj3Ecv/9+/N//dqDaxi1YRRWD16NKp5VjO5LJ79O6FKjCz4O+jjf11esAEaOBLRao5vOwVyTZX32
   GTBpElCmjHjcvLnYWf7VV8Uu7pQXM03qZrdF0/XrQI0aTz/v4wNs2QLMmSO2JCAi4xiSZ9LTaMRm
   l9/ks9DtYcZD9F/dHx93+hgdqhtQgRVgfrf5+PPcnzh191Se57OzDT82pTDMNVlOWBiwebMYkcxt
   xgzx52fuJqlE1ma3RVPuPNOTatYE1q4VP8mcPWvdfpF9YqbpsX37jJvuGjVKjExdufL4OVmWMWbj
   GLT0bYk3WrxhVn/KuZfDl89/idc3v45s3ePA4t69QKVKQMOGZjWP554DbtwwfvsEKtonn4hRJW/v
   vM9rtSIY/n//B5w5Y5OuqRYzTepmt0WTfuVcQdq2FZvv9e3LvUGIjFHY/kz5cXMTId/cm13OPTAX
   txJu4afePxkc/C7MmGZj4KR1wpITS3KeMycAnpujo1jVxVyTss6dEwX45Mn5v+7nJ0aaRowQU7xE
   9qDYFk2AWEk3bhzw4osiIE5UEGaahIwM4MgRw/JMuU2cKAK+cXHA1itb8WPIj1g7ZK1Rwe/CaCQN
   FvdZjI+CPkLUwyg8fChOAhg2TJHmmWuygI8/FuFvD4+C3zNiBPDMM8D771uvX2rHTJO62XXRlF+m
   6UmzZwMNGoifSBk6JCpcSAhQt+7T0ylF8fUF+vQBvvrlCsZsHIO/B/+Nyp6VFe3bM+WfwdhmYzFt
   xzSsXy8Ku/LllWmbuSZlnTgBHD0qdo4vjH7V87p1wM6d1ukbkTnsumgqaqQJEP8oly4Vy5L50wwV
   hJkmwdCtBvLz2uQkfBPVHx93+AztqrVTsls5Puz4IQ7dOoRvNu1SZGpOr1kzsSI3Jka5NkuyDz8E
   Zs0CXF2Lfq+PD/Drr8B//sNcGcBMk9rZZdEUHy/2aSpb1rD3OzuLE9DXrRMFFBHlz9SiSSfr8E34
   qyiX2g4+119Xuls53J3c8UnLH3Gmypvo1jNNsXYdHID27QveOoEMd+gQcOGCOBPUUM8/D7z8MvD6
   64Xv+UVka3ZZNOlXzhmTLy1TRmxFMHs2sHu35fpG9omZJtPzTADw5b9fIuphFL7v8QO++cay//FF
   H+iDak6N8d2JuYq227kzc01K+PBD4KOPxA+rxvjiC7EC89dfLdItu8FMk7rZZdFk6NTck+rWBf7+
   Gxg+HAgNVb5fRPbM1DzT5subsfj4YqwdshYDXnTGgwditMESZFksVV/Q9Tv8eOxHhMWGKdZ2QABz
   Tebau1dMc44aZfy1zs5iMcG774qjbYjUyG6LJkNC4Pnp1AmYP1+EVqOjle0X2S9mmkybmguLDcN/
   Nv4HgUMCUalUJWi1BW92qYTTp8VK2AFdqmBWh1l4c+ubRR7oa6hmzYDISODePUWaK3FkWYwyffJJ
   3pMajPHssyIL9corQFaWot2zG8w0qZtdFk2FbWxpiFdfFaNN/fsDacrFIojsmrGbWiamJ6L/6v74
   ssuXaF2ldc7zo0eLAiw8XOkeilGmkSPFTuSTW03Gg9QH+OvcX4q0rdWKqUnmmkyzfbvIm778snnt
   TJkCuLsDX32lTL+IlCQp9VNagTeQJFnpe/ToIf5h9eplehs6nSicAOCvv8Q3YaKSKj1dLKy4dcuw
   6TmdrMOA1QNQuVRlLOq96KnX331XLNZQcsQpKwuoXBk4cEAc8wIAR28fRf/V/RH6Zih8XA0496UI
   CxeKXM3PP5vdlN35OeRn7AnfAzdHN6M/XB3c0L+3G6ZPdcPwl1yh1ZhxGCDEYezNmgH//AO0bGnY
   Nenp6di+fTu8vLzMujeVbPHx8RgwYABkWc43NW2XRVPdumJju/r1zWsnNVWs2ujaFfj8c2X6RmSP
   DhwQ02rHjxv2/k+DPsXu8N3YM2oPnLROT71+6xbQpIk4nsTTU5k+bt0qzpR8Mi/15pY3oZN1WNxn
   sdn3OHlSbLh48aLZTdmVbVe24bXNr2FBtwVIz05HSmaKUR+xCSmIe5gCZw/x2EnrlG9x1btOb8zq
   OMugPgUGAh98IP5MCtsgU49FEymhqKLJwdodMld2tgga+vmZ35arK7BxI9C6tfjJ1ZTwIhUPQUFB
   JXoFnTF5pk1hm/DfU/9FyPiQfAsmAKhaFejeHVi27OnDWk31++/5H5vyZZcv0fCnhni1yatoU7WN
   Wfdo0gSIihIfFSua1ZTdiEiIwJiNY7DmpTUmHays04mv2/++EllRWZYLLLy8XQxfZTB4sFjxPG0a
   sGRJ0e8vLk6fPs0VdCpmd5NSkZFi+wAXZU5nQPny4hTu6dOZZSiJdLrHx3+UZIYWTZdiL2HcpnFY
   O2QtKnoUXlVMmybOo1Mi0JuQAGzbBgwd+vRr3i7e+Lr715iwZQKydObdTKsFOnYsOavoMrIzMDRw
   KKa1mWZSwQSIFcnu7kDv3uKxJElwcXBBadfSqOJZBXXL1EXTik3RtmpbNCxn3OnK330ntojZuNGk
   rhEpzu7pFvzzAAAgAElEQVSKJnND4Plp2FDkmoYMAS5fVrZtUq9bt8RoyNixwI0bAbbujs2kp4sj
   L9q3L/x9CWkJ6LeqH+Z1nYeWlYsOmrRsKTJIGzaY38fAQKBLF6B06fxff/mZl1HevTy+O/Jd/m8w
   Qkk6UmXmrpko51YO09tON+n6rCxxxtycOcbtm2coT09xMPPrr4vRv5KAo0zqZndFk6l7NBWla1fx
   D793b+D+feXbJ/WQZeCPP4DmzUWmbdmykr2pYUgIUK9e4QFwnazDyPUj0b1md4xpNsbgtt9+W5kw
   +O+/Fz59LkkSFvVahK8OfIWIhAiz7lVSDu9dG7oWG8I24Lf+v0EjmfZfwR9/iGnMLl0U7lwu7doB
   48eLY1a4WzjZGoumXMaPF9sQDBggfvqm4ic2FnjpJWDuXGDHDhE07doV2Ls3CNnZtu6dbRQ1NaeT
   dZi1ZxYS0xOx8IWFRrXdv7+YUj961PT+3bghNqMtarVsnTJ1MKXVFEzdPtX0mwFo3Fj8Pblzx6xm
   VO3qg6t4Y8sb+Hvw3yavOszIAD77TCyiscQoU24ffSTOBSwJqxq5T5O62WXRZOrGloaYN09kpsaP
   5081xc2WLeI/RD8/sUqsWTPxfIUKYrn9qVM27Z7NFFQ06WQd1l1ch2a/NMOe8D1Y89IaOGqN27XQ
   wQGYOtW80aY//hBT5075Z87zmNluJkJjQrEpbJPJ99NoineuKTUzFS+teQkfd/oY/pX9TW5n+XKx
   gKZjRwU7VwBHR5E9/PjjkreykdTF7oomS2SactNoxDfp0FBxFhLZv6Qk4LXXgEmTgJUrgQULnl5I
   0KdPQImYknmSPs+U+7w5WZax8dJGNF/SHHOC52BO5zk4Ou4oyruXN+keY8cCO3eKVa/G0h+bYujK
   VmcHZyzqtQhTtk1Bckay8Td8pDjnmqZun4p6ZerhTf83TW4jLU3EGay5VUvduuKeI0eKUa7iipkm
   dbO7osmS03N67u5iU7WlS8V/smS/DhwAmjYVW1WcOSOO0clPST2s9dgxsd+Zl5colv4J+wctlrbA
   x0Ef45NOn+DEayfQt15fSGbMv3h6il34f/zRtP4Bhm9wCABdanZB+2rt8en+T42/4SPFNde04swK
   7L+5H0v7LjXrz/SXX4DnnjPuz0UJr70G+PqKo1qIbMGuiqbkZLH02Br7p1SqJAqn114TOxuTfUlP
   B2bOFPmlhQtF2LuwTRYdHIJw8GDJ+7MOCgI6BcjYemUrWv63JWbtnYVZHWbh5Osn0a9+P7P+Y81t
   yhTgf/8DHj407jr9KJOx3fi6+9f49fSvOHfvnHEXPvLMM2Ibitu3TbpclS5EX8C0ndMQ+FIgSjmX
   Mrmd5GSRCfzsMwU7ZyBJEv+Wly8H/v3X+ve3Bmaa1M2uiqbwcJFH0R95kpKZYvI3RUM0biyyLjdu
   WOwWZAFnzgD+/mL7iDNngH79ir7Gy+tx1qmkkGUZa89sx+YKbfDurncxs91MnJ5wGgMbDDR5NVVB
   atQQozfLlxt+TUaG2ANo5Ejj71fBowI+7/w5JmyZAJ2sM/p6jUaMShaXKbqHGQ/x0pqXML/bfDxb
   4Vmz2vrpJzGd26SJQp0zUvnyYhbglVfED9FE1mRXRdOTU3PLTy3HqA2W3ca7Th1xFhWpX3a2+Am4
   a1fgnXeAdevEN1hDBAQElJgpOlmWsevaLrRd1g5nfadhZse3cfaNsxjccLDixVJub78tNis0dJXi
   1q1iDzVTd/8f33w8dLIOy04uM+n6gIDiUTTJsozXN7+O1lVaY3TT0Wa1lZgoMoGfmj7zqYg+fcRq
   ykmTbNsPS2CmSd3sqmh6MgQeeDEQ5+6dQ2pmqsXuyaLJPly7JkYGduwQo0Wvvmr8lM7zzxfvokmW
   Zey5vgcdlnfA5G2T0d1rMp47fA5jWg61aLGk17at2Jxy82bD3m9MADw/GkmDxb0XY9beWYhOjjb6
   +uJSRC85sQTn7p3Dj71MCJU94dtvxYHpDRoo0DEzLVgg9hhbtcrWPaGSxK6KptwjTfce3sOpu6dQ
   r2w9nI6y3BwwiyZ1k2VxLlWrVsCgQcCePUD16sa3ExQUhI4dgSNHiuceXUE3ghDwWwDe2PIGJrSY
   gAtvXoDDpWHoHGDeafTGkCTDN7t88ED8WQ4ebN49m1RsglFNRmHGrhlGX9uwoRhZMWXVn1qcvHsS
   s/fNRuCQQLg5upnV1oMHwPffi2X/auDmJlY6T5kidvcvLphpUjeziyZJknpIknRJkqQrkiTNVKJT
   BcldNG24tAE96/REh2odEHInxGL3rFMHuHrVYs2TGe7eFcP0v/wCBAeL/5A1ZvyN9vYWO2PrV2wV
   B8E3g9H5t84Yt2kcxjYbi9CJoRjZeCS0Gq1Rh/QqZfBgMSp48mTh71u9GujZU2TNzPVJwCfYF74P
   e8P3GnWdRmPfU3TxafF4ac1L+KnXT6hbpq7Z7X39tdj4t1YtBTqnkBYtgLfeEiPLOuOja0RGM6to
   kiRJC+BHAD0ANAQwTJIkiw3c5t7YMvBiIAY3GIwWvi1w/I7l0rscaVKnwECxOWXz5sDhw2JUwBwB
   j6qH4jIlczDiILr+3hWjN4zGqMajcGnSJYxqMgoOGgcAYjTt2LGiz5tTmqMjMHly0aNNK1aYNzWX
   m4eTB5b2XYqhgUOx5sIao66116JJlmWM2TgGvWr3wpBGQ8xuLzoaWLwY+PBDBTqnsJkzxaKB77+3
   3qipJTHTpG7mjjS1BHBVluUbsixnAlgFwIC1SsaTZbGKrUYN4H7KfRyLPIYetXvA39ffoiNNNWqI
   YyCK82ZqaibLMpadXIZsnUgPx8WJ1VQffCBOPv/sM8N2ijaUvRdNh28dRvcV3TFi3QgMe2YYwiaF
   YUyzMTnFkt7RoyKXosRIjrHGjxe5poKOKblyRfyA1L27cvd8ofYL2D5iO2bsmoGZu2YiS5dl0HX2
   +vfhmyPfIDIxEgu6L1CkvXnzgOHDgWrVFGlOUVqtKLIXLNAiPLyQfUWIFGBu0VQZQO7Z5NuPnlPc
   vXti08lSpYCNYRvRrWY3uDu5o1H5RriVcAuJ6YmWuC0cHYEqVUQInaxvy5UtGPfPOKy7uA67d4tl
   zt7e4siTVq2Uu0/Qo+GEDh1EuDTVcmsLLOJY5DH0/LMnXl77MgY3HIzLky9j7HNjCzz2xBZTc3o+
   PsCIEWLpen5WrACGDRNHsCipuW9zHH/tOE5GnUSPP3ogJjmmyGsaNABSUuxr25FDtw5h3sF5+Pul
   v+Hs4Gx2e5GRYquIDz5QoHMWUqMGMHduFr74ohUWLaqF/fvLISbG/N+7LTDTpG7mfluy2ulsufNM
   gaGBGNVEjN07aBzQpGITnLhzAp1rdLbIvfVTdPXqWaR5KoAsy/hs/2cY1vAVTFo1H46/Dsb/lkmK
   jkA8qVQpsbHh4cNiNZ09OHr7KPqs7IPPO3+ODUM3GPQfZVAQMH265ftWkKlTxWq6WbNEoFdPpxNF
   07p1lrlvWbey2D5iO2btnYUWS1tg7ZC1aOHbosD3S5IoLvfvN33rA2uKTYnFy4EvY9mLy+Dn7adI
   m19+KY7CqVRJkeYsZsQIHaKjT+DatcrYsaMCFi6sC2fnbDRsmIhGjcRH7dpJcHLioaK2JMtAVpaE
   zEwNMjM1yMjQICPj8WPxnISMDM1Tz+l/rdNJ0OkAnU6CLOf9rH9NlvN+NvS96emFTyuZWzRFAqia
   63FViNGmPEaPHg2/R99xvL290bRp05wMif4n/KIe374dgBo1gH92/IOgoCCsHrw65/VKsZUQcicE
   nWt0Nrg9Yx67ugJXrijXHh8b9njntZ24dz4a0Zu+QurwN7FsazCc7suPRkksd/9atYB9+wLw/PPq
   +noU9HjRsUV40/9NTGgxwaD3Z2QAx44FoH172/a/bVtg9uwgvPji49d//FG83rSpZe8/t+tc+Pv6
   o8tnXTCh+QTMGz+vwPf7+oq/D6++qo4/74Ie62Qden3RC+1Lt0efun0UaX/VqiCsWAFcu2b7319R
   jyUJyM7ej2ee8cDIkU0hy8Du3aEID3fH7dtdsXNnBdy8eQyVK6fC378JGjVKhCwHwds7MydHpB/l
   sfVjPVv2R5aBY8fOIC1Nixo1WiI5WYuzZ08hNVWLChXaIDnZAVevhiA9XQt39/ZITnZAVNRhpKVp
   4ejYCZmZGjx8+C+ysjSQpABkZGiQlhaM7GwJWm1HODrK0Gj2wcFBB3f3DnB01CErKxharQ4+Pu3g
   5KRDSsq/cHCQUa5cGzg66pCYeAAODkDFiq0hSTJiYg5Do5FRqVJrSBJw795hSBJQuXIraDTA3bvi
   cZUq4vGdO0eg1cqoVq0VJElGZORRSJJY9BERcRSJiZHIyip86l6SZdOrbkmSHACEAegC4A6AYwCG
   ybJ8Mdd7ZHPuoffZZyK8Wm/I71h3cR02vLwh57U/z/6JDWEbsOYl40KehvrhB3Gy9qJFFmme8iHL
   Mtr9rx3qJ05G9J5hePHTJdgUtgmbhxu4yY8Zdu0Sm/cdOGDxW5lNlmXU+r4W1g1dh6YVDQuQBgeL
   USZbrxIMCgImTBCHY2seBQVeew2oXRt4913r9CE0JhQDVg/A837P49se3+Y7SnfpktibSO1TdHOC
   52DntZ3Y++repzJspho3ThxbNWeOIs1ZVHp6OrZv3w6vQoJ6qalaXLpUChcueCI01BMXLngV29Go
   7GwJSUkOSEx0QFKSY87npCQHpKQ4IDlZi+RkB6SkiM95f61FSooDnJx0cHPLgrt7Ntzds+DmJj7n
   /XV2nve4umbD2VkHR8fHH05O8qPPOjg46KBVcWY/Pj4eAwYMgCzL+e70Z9a/LFmWsyRJmgRgBwAt
   gGW5CyYlhYeLlT6BoYFPrQbxr+yPWXtnWeK2AMQ38U2bLNY85WNv+F7EptzH7u+G4O9VQNMmo/DR
   vo8QGhOKhuXMXCqXj6CgoJyfXtu1A06fFmdsubsrfitFnb13FgDQpILhZ1rYMs+UW6dOgKsrsH27
   2N05NVWsijx71np9aFiuIY6NO4ZXN7yKgN8CEPhSICp75o1l1qsnfmALD3+8etcW9u8XQez8+rA3
   fC8WhSzC8deOK1YwXbkCbNhQvFYPu7pmo1mzeDRrFg9ATBVFRrrmFFE7dlREZKQratV6mFNINWyY
   iHLlrLd52+nTp/OsoMvMlPIUPQV9FgXS48+pqVp4eGTB0zMTnp6ZKFVK/NrDQxQ4Xl6Z8PVNe1QA
   PS56chdBWq39F49KM/tflyzL2wBsU6Avhbp+HRg0PBFBx4Pw+4Df87xWu3RtxKfFIyY5BuXcyyl+
   b247YH2fB3+OjpiNm3W0aN0aAFwwqeUkLDi0AP/r9z+L3tvNTWxncPCgsiu4LGH9pfUYUH+AUQfr
   2jrPpCdJwLRp4kDlXr3EAdnNm4uFF9bk5eKFdUPXYe6BufBf6o9Vg1ehY/WOefqp33rAVkXT1q3A
   6NHi135+wNChwJAhQNWqwJ2kOxi5biRWDFgB31K+it3z00/FHkg+Poo1qTpi6iYVVaqk4oUX7gHI
   Oxr1ZDaqYcNEuLhkIztbypWJkR49Rp7nn3zu8Wv5P6d/HB2dBKA5kpIckZDgiKwsKafgyf25VKlM
   eHpmoUaNZHh6Pn6s/+zmlmXWvnWUP4XXp1jO9evAdYct6FC9A7xdvPO8ppE0YmXMnePoWaen4vf2
   8wOiooC0NMDFRfHm6QnBN4NxO/E2bv84DP/77+Pn32jxBur8UAdznp+j6H8OwOOMhJ5+qbnai6Z1
   F9dhUW/D543T0sTqQGvvz1SQoUPFPjtnzyq7N5OxNJIGH3T4AM0rNcdLa17CrA6zMLnl5JxiNCBA
   /H0YM8b6fYuOFtNka9aIUdB9+8TRIU2bAvUbZuFej2EY0WQCutTsotg9Q0PFNPXPPyvWpN0obDQq
   LKwUMjI00GplaDR49Fl8aLVyzvNOTro8z4vPyPP48TXIeazRyHBzqwpPz7Cc4sfFJdvoI6HIcszK
   NBl0AwUyTenpgKcn0Pv3Qehbtw/GNHv6O9d7u9+Dm6MbPur0kVn3KkjdumKo2txNFKlo3VZ0Q42H
   w3Bp5X8QHJz3tSnbpsDN0Q1zu861aB/27QPef18cq6JW1x5cQ9v/tcWdaXeg1RgWEggOBmbMEPs0
   qcUXX4iv87//ArdvAx4etu3P9bjrGPT3IDQs1xBL+y6Fm6MbLl8WB0HfvGn8mYbmkGWgb1+gcWOx
   ii23jAxgxPL3cfD6SaQs2YZmTTUYOlQcJ1TOzAH3l14CWrYUf1fshSGZJqKiFJVpsovBu5s3gco1
   HmL39V3oVz//vTMtvcklp+is49CtQ7hy/wr+/fGVfHcffrv121h6cimS0pMUva9+FY5emzbA+fPi
   7DG1Wn9pPfrX629wwQSIYlANeabcXn8d2L0bePFF2xdMAFDTpyYO/ucgNJIGbZa1wfW466hTB8jK
   EiPe1vTzz2KPuk8+efq1nTc240jKHzjz0R+IuqvB1Kki91SnjhghXbZMnBdnrFOnxNT0xIlmd59M
   wH2a1M0uiqbr14FSzbahTdU2KO1aOt/3+Ff2R0hkCCw1csaiyTo+D/4cL7i9D08PR3Tt+vTrNXxq
   oGvNrlh6cqlF++HiAvj7i9EPtVp/aT0GNBhg1DVqCYHnVrasOLF+2jRb9+QxN0c3/N7/d4xrNg5t
   lrXB9qvb0LmzdY9UCQ0Vh+P++efTu97fjL+JsZvGYtWgVSjnXg4uLkD//sDKlWKn9fHjRQ6qRg1x
   PuOKFYb/APDRR2KU1c28832JiiW7KZoeVhdnzRWkqmdVyJBxO/GpbaIUwaLJ8kIiQ3A++jwOLx6N
   2bMLngaZ0XYGvj3yLTKzMxW795OZJkBsbqnWIzTuJt1FaEwonq9h+A6c+jxTu3YW7JiJJk4UGR01
   kSQJk1tNxtohazHun3FIbj4He/dZ51TY9HRxbMlXX4loQG4Z2RkYEjgE77Z9F+2qPf2H6eYmptfW
   rhXTncOGiTxU1ariwN1Vq8TK0PwcOQKcOSO2fiDb4Nlz6mYXRVPY9VTccd2B/vX7F/geSZIsOkXH
   osnyPg/+HL29ZkKjc0afPgW/r4VvC9QqXQt/X/jbov1R87ljG8M2omftnnDSGn7w3tGjQKNGIh9I
   hmtfrT1CxofglvM2rHMcgPjUBIvfc9YscQLC2LFPvzZ953RU8qiEaW2KHporVUocWbNpk4g59OsH
   /PYb4OsrVt+tXZv3yKCPPgJmzwac7fMEEiKLs4ui6WjsDtRye67I7QT8fcUUnSWwaLKsU3dP4cTd
   EzixdFyho0x6M9rOwPxD8xWbjn0y0wSIIOzly+KQYLVZf2k9BjYYaNQ1apyasxe+pXxx6PV90CZX
   QbOf/XEh+oLF7rV7txgNWrr06X8Hay6swZYrW/Br/1+N2mYCEGc2jh4NbNsGXLsmgu2LFonjUUaM
   AObPF8/bYoUgPcZMk7rZRdF02SEQfWoWPDWn51/ZciNN1aoBMTHi8E5S3px/56BvmelITnDBQANq
   gZ61eyJLl4Xd13dbrE9OTiIQ/uQKPluLS43D4VuH0aN2D6OuY9FkHmcHJ/R3+gkB2g8Q8FsA1lxQ
   /gSC+/dFYfPrr0CZMnlfu3z/MiZunYi/B//91LYrxipbVkzB7dkDhIWJcwC3bgXmzROHlBNR/lS/
   T1NaZjriym3Bq63mF/neFr4tcOLuCciybPRPYUXRakWo8to14NlnFW26xDsffR4HIw6i1uYVmDUL
   Bm3IJkkSpredjvmH5qNbrW5m9yG/TBPweIquX/6LNm1iy5UtCPALgIeT4UvN1JxnsiedOwN79ozG
   9gXPYtDfgxByJwRfdvlSkV24ZRkY+1oGeg67A5e6t7Hq/G2xX9mjjyO3j+Czzp+huW9zBX4nj1Wo
   IDJlXC2nDsw0qZvqi6Z1Z3ZBG/ssGlQp+ojt8u7l4ensiasPrqJOmTqK90U/RceiSVlzguegf4Vp
   2H3bDUOHGn7d8GeHY/be2Tgdddrgc9eM1bmzWBKvJqZMzTHPpIyAAODDD4E/KzXH8deOY9jaYXjh
   jxdyVrEVJiUzBZGJkXkKoduJt3E7SXy+cu82Hj4Th8reFXFhVxVU8RQf1byqoW3VtpjZbiZa+Law
   zm+UiPKl+qLpr9OBqBhX9NScnj4Mbqmi6epVxZst0S7GXMTe8L14Zs9/8cEHgIMRfyOdtE6Y0moK
   FhxagD8G/mFWP3KfPZdb8+bioNbYWDGlYWspmSnYfX03lvRZYtR1QUGiACTz1Kwp/o5evgzUq1cW
   20dsx6y9s9BiaQv80ucXaCTN00XRo4+UzBRU9qycUwxVKVUF9cvWR9eaXZEdXwWvf1EFZ7aWx7ON
   VHyaKVnck2fPkbqoumjKyM7A/qh/0MHB8CO29WHw4c8OV7w/deoAJ04o3myJ9uWBLzGw8lRsv+SB
   kSONv/715q+j5vc1EZEQgWpe1RTvn6OjmNIKCgIGG167W8zOazvRvFJzlHErU/Sbc9m3D3jvPQt1
   qgTJfQ5dvXqAVqPF3K5z4e/rj/d2v4eybmVziqLnKj2HvnX75jwu61Y239hARobIFH0+E3i2kdV/
   S0RkBFUXTfvC98FHVxfPVDP8BE//yv74aJ9ljlKpU0esaiFlXLl/BduubEOLQz9i5synN/AzhJeL
   F8Y0HYNvj3yLhS8sNLkvBWWagMf7NamhaNIf0GuMtDTg+HHmmZTSuTOwY0feadtBDQdhUMNBJrX3
   ySdAxYrAG28o0z+ybxxlUjdVr54LDA2Eb/xg1Kxp+DXNKzXH6ajTyNJlKd4fbjugrK8OfIXB1Sbj
   XIiXWcucp7aail9P/4r4tHjlOpeLWvZryszOxObLmwvdryw/R44Azzwj9uwh8+lHmpTY7WL/frFS
   7n//s+6ZdkRkGtUWTVm6LGwI2wBt2CDUqGH4dV4uXqjiWQWhMaGK96lKFbFnT0G76ZLhwuPCsTFs
   I+6sm4IZM8SxJaaq6lUVfer2weLji01uI799mvSaNgXu3gWiokxuXhHBN4NRy6cWqnpVNeo6bjWg
   rBo1xN/XsDDz2omLA0aNEmfElS+vTN/I/nGfJnVTbdG0/8Z++Hn74e5FP6NGmoDH59ApTaMRQVCG
   wc0398BcDK7+Bo4F+2D8ePPbm952Or4/+j3Ss9LNb+wJWi3QsaN1zx3Lz7qL64yemgNYNFlCQIB5
   o4+yDEyYILay6NlTsW4RkYWptmgKDA3EwHqDERkJVK9u3LUtKrXA8TvHLdIvTtGZLyIhAmtC1yB2
   81t46y3A3d38NhtXaIzGFRrjz3N/mnR9YZkmwPZTdDpZhw1hG4w+oJd5Jssw9/DeFSvEgbzz5inW
   JSommGlSN1UWTdm6bKy/tB6tvQahYkXjA8KW3Bm8dm0WTeaad2AeBtcYj+DtZfHmm8q1O6PtDCw4
   tAA6WflDVW1dNIVEhsDL2Qv1y9Y36jrmmSzDnFzTtWvAO+8Af/0FuLoq3TMisiRVFk0Hbx1EpVKV
   oIutbfTUHAA0rdgUoTGhFpmq4UiTeSITI7Hy/Eok7XwHkyYpu9ni8zWeh4uDC7Ze2Wr0tYVlmgCx
   oemDB+LUeFswZdUcwP2ZLKV6dTFCevGicddlZQEjR4pDcblJLuWHmSZ1U2XRFBgaiEENBuH6dRgV
   Atdzc3RD3TJ1cebeGcX7xqLJPPMPzcfAmqOxc315TJmibNuSJOUc5Ks0jQbo1Mk2o02yLIs8k5FT
   cwDzTJZkSq5pzhzAywuYPNkiXSIiC1Nd0aSTdVh7cS0GNxyM8HCYNNIEPN7kUmksmkwX9TAKv5/5
   HRlBM/D664CPj/L3eKnRS7gZfxPHIo8ZdV1RmSbg8X5N1hYaE4q0rDQ0r2TcmWOpqcwzWZJ+is5Q
   Bw8CixcDy5cbdr4ilUzMNKmb6v7pHrl9BD4uPqhftj6uXzejaLJQrsnXF0hKAhITFW+62Pv60Nfo
   V3MENq+shLfftsw9HDQOeLv12xYZbbJVrkk/NWfsIdRHjogpIA/Dz/UlI+iLJp0BEbqEBDEtt2QJ
   UKnoYzSJSKVUVzQFhgZicEOx9bJZRZOvZYomjQaoVYvbDhgrJjkGy04tg/bwTIwZA5Qr/GxTs4x9
   biyCbgTh2oNrBl9TVKYJABo0EKM3N26Y3jdTcGpOnapVE5m8UAO2hJs0CejRA3jxRcv3i+wbM03q
   pqqiSZblnKk5ACZnmgDgmfLP4Eb8DSSlJynYQ4FTdMb75sg36FtzCNb9WgXTp1v2Xh5OHnjtudew
   8LDpx6rkR3/umDVHm27E38CtxFtoX6290deyaLI8Q0Yf//oLCAkBvv7aOn0iIstRVdF0/M5xuDq4
   olG5RkhMFD/Vm7pTrqPWEc+Wfxanok4p20mwaDLWg9QH+OXEL3A7+R6GD7fO9MTkVpOx8vxKxKbE
   GvR+QzJNgPWn6DZc2oAX674IB41xx0SmporDpZlnsqyick03bgBvvSUKJzc3K3WK7BozTeqmqqJJ
   PzUnSVJOCNyc85gYBleH7458h55+/fH3Ej+8+6517lnRoyIGNRiEn479pGi7+qJJiXPHDLH+0nqT
   puaYZ7KOgABxflx+uabsbOCVV4AZM4DnnrN614jIAlRTNMmyjMCLYqsBwLw8k56lwuAsmgyXkJaA
   n0J+gs+FDzBggMiBWMs7bd/BouOLkJKZUuR7Dck0AeLPXqcTGxRaWnRyNM5EnUHXml2NvpZTc9ZR
   pYpYBXr+/NOvzZ0rNuZ95x3r94vsFzNN6qaaoul0lPiL0rSiGJpUpGiyUBicRZPhfjj2A7r69cLK
   n2rh/fete+/6ZeujdZXW+O30b4q1KUlitGnvXsWaLNCmsE3oXqs7XByMP82Ym1paT35TtseOAd9/
   D/z2G7cXICpOVPPPOTA0EIMbDM5ZVm1OCFyvXtl6iE2Jxf2U+wr08LGKFcWZXvHxijZb7CSlJ+G7
   o0hih1kAACAASURBVN+h0pVZ6NlTrDq0thltZ+Drw18jW5dd6PsMzTQB1tuvaf2l9RjYYKDR1+nz
   TG3bWqBT9JQnc00PHwIjRgCLFomRKCJjMNOkbqoomvRTc/pVcwDM2thSTyNp8Fyl5xQ/vFeSxBl0
   3HagcItCFqFztW746/t6+OAD2/ShXdV2KOtWFhsubVCsTWvkmhLTE/HvzX/Rq04vo689fBho3Jh5
   Jmt5Mtc0dSrQsSMwaJBNu0VEFqCKoulCzAWkZaWhhW+LnOeUmJ4DOEVnK8kZyVh4ZCGq35yFTp3E
   Hke2kPtoFbmQKsfQTBMgRkBdXIBLlxToYAG2XtmK9tXaw9PZ+MP5mGeyLl9fse/Y2bNAYCAQHAx8
   952te0X2ipkmdVNF0fTk1JxOJ5bq+vmZ3zaLJttYfHwx2lXpiD+/bYRZs2zbl/71++N+6n0ciDig
   WJuW3nrA1Kk5gEWTLXTuDPzxBzBxIvDnnxzlIyquVFM0DWr4eCz7zh2xIkWJfU1a+LbgtgNWlpqZ
   igWHF6Bu1Gz4+wNNmti2P1qNFtNaTyv0aBVjMk2AZYumtKw07Li6Ay/WM3776JQU4ORJ5pmsLSBA
   bF45ZQrQsqWte0P2jJkmdbN50XQx5iLi0uLQukrrnOeUyDPp+Xn7ISM7A3eS7ijT4CMsmgq29ORS
   +Fdqhb8WNsHs2bbujTC66WgcjTyKizEXFWmvc2fDzx0z1p7re9C4QmOUdzd+Z9cjR5hnsoUXXgCm
   Twfee8/WPSEiS7J50bT24loMajAIGulxV5TKMwEi0+JfWflNLlk05S8tKw3/d/D/0CT+QzRsCPj7
   27pHgqujK95s8Sa+Ppz/WRbGZJoAoGpVwMsLuHBBgc49Yd3FdRhQ3/gNLQFOzdmKjw8wfz6g1dq6
   J2TvmGlSN5sXTbkP6NVTsmgCLJNrKlcOyMoCHjxQtFm7t/zUcjQu3wR/fd1cNaNMehNbTsTai2tx
   N+muIu1ZYr+mLF0W/rn8j0m7gAPcn4mIyJJsWjRdfXAVUQ+j0K5q3gOy7KFo0m87wNGmxzKyMzD3
   4Fz4p36IqlWB9safMWtRZd3KYvgzw/HDsR+ees3YTBNgmf2aDkYcRBXPKvDz9jP6WuaZiOwfM03q
   ZtOiaW3oWgxsMBBaTd4xbSU2tszNv7I/jt85XuiSc1Nwii6v38/8jrql62H1163x4Ye27k3+prWZ
   hiUnliApPcnstgICxPLy7ML3zTSKOVNzR46I0L27u3L9ISKix2xaND25oaWekkFwQBze6ubohutx
   15VrFCyacsvMzsSX/36JDrqP4OMjRmHUqFbpWuhcozOWnVqW53ljM00AUKkSUKECcOaMMn2TZRkb
   wjaYvNXAvn3MMxHZO2aa1M1mRdON+Bu4EX8DHat3zPN8SorICfn6Kns/S0zRsWh67K9zf6G6V3UE
   LmyPDz8U05dqNaPtDHxz5BtkZmea3ZaSWw+cvHsSzlpnNCzX0KTrGQInIrIsmxVNa0PXon+9/nDQ
   OOR5Xr+ppdKHXFpivyYWTUK2Lhtf/PsFujh8BAcHoGdPW/eocC0rt4Sftx/WhK7Jec6UTBOgbNG0
   /tJ6DKg/IGeTV2OkpACnTjHPRGTvmGlSN5sVTQVNzSmdZ9Lz9/XH8bvKnkGnL5oseQaZPVh9YTXK
   u5fHhu8CMHu2ukeZ9Aw5WsUQAQHAgQNiJaW51l1cZ/KqucOHmWciIrI0mxRNtxNv4/L9y3i+xtPB
   F6XzTHotfFvg1N1TRZ52b4wyZUSBEBurWJN2RyfrMCd4Dnq6fYS0VAn9+9u6R4bpVacX0rPSsSd8
   DwDTMk2A2HqialWxas0cYbFhSEhPQMvKpm0nzak5ouKBmSZ1s0nRtO7iOrxY70U4ah2fek3p7Qb0
   fFx9UMGjAi7FKnfKqiRxim7P9T1wcXDBlh+6YdYs5adVLUUjaTC97fRCj1Yx1PPPm79f0/pL69G/
   Xv88m7wag0UTEZHl2eS/OP0BvfmxVNEEMAxuCTuu7UBj536IjZEwZIite2OcEc+OwLl753Am6ozJ
   mSZAmVzT+kvrTZ6aY56JqPhgpkndrF403U26i3PR59C1Ztd8X7d40cQwuKJ2X9+N8xu74YMP7O8I
   CWcHZ0xpNQULDi8wq51OnUSmKCPDtOtvJ97G1QdX0al6J5OunzNHjHYxz0REZFlWL5rWX1qP3nV6
   w9nB+anXZFlkmiwRBAfEJpccaVLOvYf3cO3+DUSfaokRI2zdG9NMaDEBWy5vwcp/Vprcho+P2B0+
   xMS/WhsubUCfun3yna4uyt9/A3/9BSxbVvR7iUj9mGlSN6sXTfmdNacXEwO4uACenpa5d7OKzXA+
   +jwysk0cEshHSS6a9oTvgR86oXdPBzga//+9Kni7eOOjTh9h5u6ZiHoYZXI75kzR6bcaMNaZM8DE
   icD69SKQTkRElmXVoikmOQYn757EC7VeyPd1S07NAYC7kztqla6Fs/fOKtZmSd52YNf1XXC63Q2t
   W9u6J+Z5q/VbGD9wPLr83gUxyTEmtWFq0XQ/5T6O3zmO7rW6G3fdfWDAAOD774FmzYy/LxGpEzNN
   6mbVomnDpQ3oUbsHXB1d833d0kUT8Gi/pjvK7dfk4wM4OwPR0Yo1aRdkWcaua7tw90A3tGlj696Y
   78NOH2Jg/YHouqIr7qfcN/r6Dh2AY8eA9HTjrvvn8j/oWrMr3BzdDL4mKwsYMgQYPBgYNszIjhIR
   kcmsWjQVtKGlnrWKJobBzRd2PwyyToO0O3VRp46te2O+oKAgfNb5M/So1QPdVnRDXGqcUdd7eQEN
   GohDc41hytTcu+8CDg7AV18Zdy8iUj9mmtTNakXT/ZT7OHL7CHrWLviMDUuGwPUYBlfGrmu7UN+p
   K9q0luxiB3BDSJKEuV3nIsAvAC/88QIS0hKMut7Y/ZqSM5KxL3wfetfpbfA1v/8O/PMPsGqV/a1W
   JCKyd1YrmjaFbULXml3h7lTwumhrjDQ1rtAY1+KuITkjWbE2S2TRdH0X3O7af55JT79PkyRJ+Lr7
   12hZuSV6/tkTSelJBrdhbK5p+9XtaF2lNXxcfQx6//HjwDvvABs2iGlhIip+mGlSN6sVTYEXC97Q
   Us8aRZOT1gmNyjXCqahTirVZ0oqmzOxM7L+5HzHHuhSboik3SZLwfc/v8Wz5Z9H7r94GF9jt2onj
   VFJSDLvPukvrMLDBQIPee+8eMHAgsGQJ0KiRYe0TEZGyrFI0xafF49+b/6J33YKnITIygKgocY6X
   pSmdayppRdOxyGOo4V0T54+WR0vTjkpTnSfPntNIGvzc52fUKl0LfVf2RUpm0ZWQh4c4NPfQoaLv
   l5GdgW1XtqFfvX5FvzcDGDQIGDNGrJgjouKLmSZ1s0rRtPnyZnSu0RmezgVvwHTzJlC5sgi4WprS
   uabatYGrV0vOtgO7ru9CY/du8PMTAejiSiNp8N++/0Vlz8rov6o/0rLSirzG0Cm6feH7UL9sfVQq
   VanI906dKg6H/vhjQ3pNRESWYpWiqbCz5vTCwy0/NafXwreFokWTlxfg5gbcvatYk6q26/oueMZ2
   LVZTcwWdPafVaLG833KUdi2NQX8PQnpW4XsKGFo0rbu4zqBVc0uWiMN4V6ywn8OQich0zDSpm1W+
   De+7sQ996/Ut9D3WyDPpNSjbAFEPo4xeVl6YkjJFl5CWgDNRZxB7skOx2J/JEA4aB6wYsAIuDi4Y
   Ejik0B3l27YFzp4FHj4suL1sXTY2hm0s8oDegweB2bOBjRstt0s+EREZzipFU/tq7eHt4l3oe6xZ
   NGk1WjSr2Awn7p5QrM2SUjQF3QhC6yqtcfywa7EaaXoy0/QkR60jVg5aCVmWMXztcGTpsvJ9n6sr
   0Lw5cOBAwW0duX0E5d3Lo3bp2gW+5/b/t3fn0VFdd57Av7/SvoFAQhLahcBgbAeBEDYOILCAOD5u
   O0zsjLOA7SwknU7aMz2ZmcTO6aQ7x5k+TvokTnJmYtw5Tuym3cYLdty2CSU2Ox0MRqwGGwxix8II
   0I72O3+8KlESJVHLq/fue/X9nKMj1avSe7eMfPXTvd937xng/vuB3/8euOGGEN4AEbkCM016s6Ro
   +vyNn7/ua6wsmgCGwSNV31iP+fnL0NxsLOYYT5ITkvHi/S+is68TK9evxMDgQNDXXW+9pustaNnd
   bdwp993vAnfdFW2riYjILJYUTaHcIWTFwpaBzA6Dx0vR5G30IqdlGebNc1fGZrRM00gpiSl45Quv
   oLmrGQ+/9nDQwmmsXJNSCq98MPpSA0oB3/oWUF4OfP/7ITaeiFyDmSa9WfJrLyc957qvsWWkiUVT
   WE63nkZzVzM+2V/lqqm5cKUlpeG1B17D6bbTWP36agyqwWHP33or8OGHQGuQBcX9m0V/Kv9TQc/9
   618De/YAzzwD16y0TkTkFlqMFVy+DAwOAhMnWnfNKROmoKuvC00dTaacb+pU4Ngx4324lbfRi7op
   ddjxrsd1IfDrZZpGSk9Kx+tffB1HLh3Bt9/4NlTAehMpKUbh9Pbb136ff2pOglREmzcDP/2pseJ3
   xugL5xORizHTpDctiib/KJOVf1mLiLH0gEm5pqwsY+mBs2dNOZ2W6hvrUVe+DO+9ZxQF8S4zORNv
   fulN7Du/D49seGRY4TTaFN36D9cHvWvuxAngS18C1q61dpqaiIhCp03RZMcvirmTzV2vyc1TdINq
   EPWN9SjtX4bCQmtHBa0QaqZppKyULLz15bew/cx2fG/j94YKp2BF07FLx3C+4zzmFw8fpuvsBD73
   OSPDVFcXUTOIyCWYadKbFkWTlQtbBqopqsGuc7tMO5+bi6b95/cjOzUbJ/eVxXWeKZjs1Gz86St/
   wuYTm/HopkehlEJNjTFde+nS1det/3A97p1+LxI8CUPHlAK+9jVj+5VHHrGh8UREFLKIiyYRuV9E
   DorIgIjMiaYRVofA/fxhcGXS/iduLpq8x7xYOmUp3n0Xriyaws00jTQxbSK8K71446M38A/b/gFJ
   ScZCl9u2XX1NsKm5J54wiqvf/pbBbyJipkl30Yw0HQCwAkCQuGt47CqaisYVIcmThJOtJ005n6uL
   pkYvlk1ZhnffhetC4GbJTc9F/ap6rDu4Do+//fiw9Zo+bv8Yhy4cwh0Vdwy9fsMG4MkngfXrjUUx
   iYhIbxEXTUqpD5VSR8xohF1FE+Bbr8mkMPi0acbGvW7T3d+N7We2oyp7Cc6cAW66ye4WmS/STNNI
   eRl52LRqE57d/yxOlz4xlGt67fBruGvaXUhOSAZgFNerVgHr1gHFxaZcmohcgJkmvdmeaRoYMLaM
   KCuz5/pmrtc0dapRALpt2YH/PPWfuDnvZhzel425c4HERLtbpLfJWZOxedVmvPXJGhyb9Et88snw
   VcDb2oB77wV+8hNgwQKbG0tERCEb89efiHgBFAR56lGl1OuhXuShhx5CeXk5ACA7OxtVVVVDf9m/
   +OJWZGUBKSnGY3+2xP98rB8nn07GxgMb8cSyJ6I+X3o6kJm5FevWAQ88YM/7icXj3zX8DktnL8X2
   bUBh4VZs3apX+8x47D9m5vk3P7gZN+y4FQ/89CB25W/Hy194GZs3b8Xf/z2wcOFifPOb+rx/PuZj
   Nzw+cOAAMjMzh0Zr/PkgJz0+evQo7rvvPm3aEw+P/V83NTWht3f0DdkBQKINQYvIFgD/Qym1e5Tn
   1VjX2LIF+PGPhwdmrXSx6yIqnqxAy/db4JHoB96WLAEeewxYutSExmmiek01fvGZX+Dx1YvwN38D
   3HOP3S0y39atW4c6XjM99vPj+FVbLe6YORuvPfAafvxjoL7eyDolJ5t+OaK41dPTgw0bNmD8+PF2
   NyUqe/fu5RSdjVpaWrBixQoopYLemmPW9FzE9/3YmWcCjC1ectNzcbj5sCnnc1sYvLmrGR9d/Ajz
   Cm/Djh3uvHMOMC/TNNJ9dRXIe307nrzzSbz6KvC73wEvvcSCiYiCY8Gkt2iWHFghIqcB3AbgDRF5
   K5Lz2LWwZSAz12tyW9G0+fhmLCpbhMaPkpGTA+Tl2d0iZ5k1C7h8sghHd5XjG98AXnkFKAg24U1E
   RNqL5u659UqpEqVUmlKqQCn12UjOY9fCloHMDoO7qWjyHru61IBbR5mA6NdpGo3HA9TWAnffDfz8
   50BNTUwuQ0QuwXWa9Gb73XN2T88B5hZNbhppUkrB23h1UUuuzxSZVauAH/wAePBBu1tCRETRYNEE
   YM7kOdh/fj/6BvqiPldlpbH5an9/9O2y29FLR9E32IeZk2Zi+3Z3jzTFKtMEACtWAD/6UcxOT0Qu
   wkyT3mwtmjo6jM1K8/PtbIWx6Wp5djne/+T9qM+Vlmbkfk6dMqFhNqtvrMfSKUvR3i44ftzI5xAR
   EcUrW4um48eB8nI99tziFN21/Fun7NwJzJ4NJCXZ3aLYiVWmiYgoHMw06c3WokmHqTm/mkJzt1Nx
   etHUP9iPLSe2uHqTXiIionCwaPKZWziXI00Bdp3bheJxxSjILIiLEHgsM01ERKFipklvLJp8ZhXM
   wpGLR9DV1xX1udxQNNU31mPZlGVQChxpIiIiggZFk90LW/qlJqbixkk3Yl/TvqjP5YaiyZ9nOnoU
   yMgACgvtblFsMdNERDpgpklvtgfBdRlpAswLg0+ZYtw91xf9Cga26OjtQMO5BiwqW8RRJiIiIh/b
   iqbBQaNo0mWkCTCvaEpJMUZmTpyIvk122HZiG+YWzkVGcobr12fyY6aJiHTATJPebCuampqA8eON
   qR9d1BTxDjrgap4JQFyEwImIiEJhW9GkUwjcb+akmTjTdgat3a1Rn8vJRZO30YtllcvQ2QkcPmys
   0eR2zDQRkQ6YadKbbUWTblNzAJDoSURVQRUaPm6I+lzTpgFHj5rQKIudaz+Hc+3nUD25Grt2Abfc
   Ykw3EhERxTuONI0wt3CuKVN0Th1pqm+sxx0VdyDBkxBXU3PMNBE5m4hAKQWllN1NiQozTXpj0TSC
   WWFwpxZN3kYvlk5ZCgBxEwInIudLSkpCZmYment77W4KuRiLphFqimqw69yuqM9TUQGcPQs46f9f
   pdQ1i1rGy0gTM01EziYiKC8vR3d3t91NiQozTXqztWjSLdMEAFMnTkVLdwsudF6I6jxJSUBxsZHd
   coqDFw4iLTENlRMrcfIk4PEAJSV2t4qIKDQFBQWumKIjfdlSNHV3AxcvAkVFdlx9bB7xmLYPndOm
   6LzHvENLDfin5kRsbpRFmGkicr6MjAyMGzfO0VN0zDTpzZai6cQJoLQUSEiw4+rXV1NoznpNU6c6
   rGjyLTUAxNfUHBG5g1um6EhfthRNuuaZ/GqK4i8M3jvQiz+f+jOWlC8BEH8hcGaaiNwhPz/f0VN0
   zDTpjUVTEP476KL9n85JRdP209sxPXc6ctJzcOUKcPAgUF1td6uIiMKTkZGB7Oxs9PT02N0UciFb
   iiYdF7YMVDyuGABwuu10VOdxUtHkbbyaZ9qzB7jxRiA93eZGWYiZJiL3KCsrc+wUHTNNeuNIUxAi
   Ykquqbzc2GPPCf/vBhZN8TY1R0TuUlBQAACOnaIjfbFoGkVNYfTrNSUmGoH3xkaTGhUjl69cxqEL
   h3B7ye0A4jMEzkwTkXukpaVh4sSJjhxtYqZJb5YXTUo5pGiKozD4lhNb8OmSTyMl0dhkjiNNROR0
   ZWVlzDWR6SwvmpqbjYUfx4+3+srh8Y80DarBqM7jhKIpcH2mM2eMVcx1L2rNxkwTkbvk5+cDcN4U
   HTNNerO8aDp+3Bm/kCdlTEJ2ajaOXjoa1XkcUTSNWJ8pnha1JCJ3Sk1NRW5uriOn6EhflhdNTpia
   86spij4MrnvRdPzycbT3tuOWvFsAxO/UHDNNRO5TWlrquCk6Zpr0xqJpDP71mqKhe9FU31iPpVOW
   QnxDS/EYAicid8rPz4eIOG6KjvTFomkMZuxBV1oKXLgAdHWZ1CiTBS410NsL7N0LzJ1rc6NswEwT
   kfukpKRg0qRJuHLlit1NCRkzTXqzJdOk88KWgaonV2Nf0z70DkS++WNCgvF+jx0zsWEmGRgcwKbj
   m4aKpr17jf3ysrJsbhgRkUlKS0sdvYEv6YUjTWMYnzoeMyfNxPbT26M6z7RpwNHo8uQxsadpD/Iz
   8lE0rghAfE/NMdNE5E55eXkQEQwORncntFWYadKbpUVTXx9w7pwxZeUUy6Ysw8ZjG6M6h665psCl
   BoD4DYETkXslJycjLy/PUVN0pC9Li6ZTp4DCQmOdJqdYXrkcGxvdWTTVHzdC4H7xPNLETBORe5WV
   laGvr8/uZoSEmSa9WVo0NTY6J8/kN79kPg43H8bFrosRn0PHoqmrrws7z+7E4vLFAIw98lpbjbYS
   EbnJpEmTHDVFR/qytGhyysKWgZITklFbXotNxzdFfA4di6Z3Tr6DqoIqZKUYqe933wVuvRXw2LIb
   of2YaSJyr6SkJEyePNkRU3TMNOnN8pEmpxVNALB8yvKock3FxcClS0Bnp4mNilLgUgNAfE/NEZH7
   lZSUOGaKjvTFoikEyyuNoinSBdI8HuN963QHnX9RSz//9inxipkmInfLzc2Fx+PRfoqOmSa9sWgK
   wQ05NwAADl88HPE5dJqiO99xHidaTmBe0TwAQH8/0NBgTM8REblRUlISCgsL0aXrSsPkCJZnmpwW
   BAcAERkabYqUTkXTpuObsLh8MRI9iQCAAweMZSDGj7e5YTZiponI/UpKStDf3293M8bETJPeLCua
   WlqMbTpyc626orncVDSNzDNxfSYiigc5OTlISEjAwMCA3U0hh7KsaPLfOefbF9Zx6irq8M6pdyLe
   UkWXokkphfrGeiyrZAg8EDNNRO6XmJiIoqIirafomGnSm2VFk1PzTH456TmYnjM94i1VdCmaDl88
   DIFg2sSrCzLFewiciOJHcXExR5ooYpYWTU7MMwWKZkuVwkKgvR1oazO5UWHyb50iviG/5mbg/Hlg
   5kx722U3ZpqI4kNOTg4SExO1LZyYadKb5dNzThbNlioeD1BZaf+yA95G77CpuR07gHnz4ndRSyKK
   LwkJCSguLtZ6io70xem5MES7pYrdU3R9A33YdnIb6irqho4xBG5gpokofhQVFWm7XhMzTXpj0RSG
   aLdUsbto2nl2JyonVGJSxqShYwyBE1G88d9Fp/vyA6QfS4qmgQHg1CmgvNyKq8VWNFuq2F00eRu9
   w1YBHxgA3nuPi1oCzDQRxROPx4PS0lItp+iYadKbJUXTuXNATg6QmmrF1WIrmi1VdCiaAtdnOnQI
   yM83/m2IiOJJUVFRxFtjUfyypGhyw9ScXzRbqthZNLV2t2L/+f1YULpg6Bin5q5ipokovkyYMAFJ
   SUnaTdEx06Q3Fk1himZLlYICoLvbWB3dattObsNtxbchLSlt6BjXZyKieOWfouvs7LS7KeQgLJoi
   EGnRJAJMnWrPsgPeY14srVg67Nj27Rxp8mOmiSj+FBYWajdFx0yT3iwrmpy+sGWgaLZUsWuKbuT6
   TC0twOnTwM03W98WIiIdTJgwASkpKejr67O7KeQQlhRNbljYMlA0W6rYUTSdbj2Ni1cuoqrg6lz5
   jh1AdTWQmGhtW3TFTBNR/BERlJWVaXUXHTNNeuP0XIQi3VLFjqLJ2+hFXUUdPHL1n5shcCIiYPLk
   yXY3gRzEkqKptdUIQbtJpFuq2FE01TfWD1tqAGAIfCRmmojiU3Z2tlZTdMw06c2Soqmiwn17m0W6
   pcrUqdYWTYNqEPWN9cMWtRwcNKbnWDQRUbwTEZSXl2s1RUf6sqxocptIt1TJywP6+oBLl2LUsBH2
   n9+P7NRslGWXDR07cgTIzjYWtiQDM01E8Wvy5Mna3EXHTJPeLCma3JZn8otkSxURa6foNhzdwKk5
   IqIxjBs3Dunp6ejtDf+OaIovLJqiEOmWKlYVTUopPLP3GXz5U18edpzrM12LmSai+OW/i+7KlSt2
   N4WZJs2xaIpCpFuqWFU0vX3ybSR5kjC/eHiFxJEmIqLh/FN0ukzTkZ6YaYpCpFuqWFU0rdm9Bqur
   V0NEho61twPHjgGzZsX++k7CTBNRfMvKykJmZqbtd9Ex06Q3Fk1RWl65HN5Gb1jfY0XRdLHrIt44
   8ga+8qmvDDv+3ntAVRWQnBzb6xMROYmOC12SfiwpmrKyrLiKPeoq6vD2ybfD2lLFXzTFchT42X3P
   4p7p92Bi2sRhxzk1FxwzTURU4FtQ0M4pOmaa9Oay1ZOsF8mWKjk5xl10zc2xaZNSamhqbiSGwImI
   gsvMzERWVhbvoqNRRVw0icjPROQDEdknIq+IyHgzG+Yk4W6pEutlB/586s8QCD5d8ulhx5XiSNNo
   mGkiIv9Cl3beRcdMk96iGWnaCOAmpdQsAEcA/MCcJjlPJFuqxLJoChYAB4wAeGoqUFQUm+sSETmd
   DlN0pK+IiyallFcpNeh7uANAsTlNcp5ItlSJVdF06colvH74daz81MprnuMmvaNjpomIACAjIwPZ
   2dno6emx5frMNOnNrEzTVwG8adK5HCeSLVViVTQ9t+853H3D3chJz7nmOU7NERFdX1lZGbq7u+1u
   BmlozKJJRLwiciDIx18FvOYxAL1KqX+LeWs1Fu6WKrEomsYKgAMMgY+FmSYi8sv3bcxpxxQdM016
   SxzrSaXUsrGeF5GHANwFoG6s1z300EMoLy8HAGRnZ6Oqqmrol5R/WsTpj5ffvBw/+8vPsGXLFojI
   dV8/a9ZiHD0KbNmyFSLmtOcvp/+C9sPtGLhxAPDtz+t/ft68xfjwQ6C9fSu2brX/vxcf8zEf87HO
   jydMmIAdO3YgOTl5qJDxT53xsbse+79uamq67p2TEmklLSJ3AvhnALVKqVFvnhcRFQ+BOqUUyn5Z
   ho0rN2JG7oyQvic3Fzh4EPD9URO1VetXoaqgCn83/++uee6dd4DvfQ/YscOca7nN1q1bhzpO0PpD
   wQAAEKZJREFUIqKTJ09iz549yM7OtvS6e/fu5WiTjVpaWrBixQoopSTY89Fkmn4NIBOAV0T2iMj/
   jeJcjhfJlipmTtFdunIJfzz8R6yatSro85yaIyIKXV5eHgDeRUfDRXP33DSlVJlSarbv49tmNsyJ
   wt1Sxcyi6V/3/yvumnYXctNzgz7PEPjYOMpERIHS0tKQk5NjeSCco0x644rgJgp3SxWziialFNY0
   jB4AV4ojTURE4SorK7Nt6QHSE4smE4W7pcrUqeYUTdvPbEffYB9qy2qDPn/qlPG5tDT6a7mVPwBK
   ROSXl5cHEbF0io7rNOmNRZPJwtlSxayRpjUNa7B6zrUrgPv5p+ZGeZqIiIJITU3FpEmTbN1WhfTC
   oslk4WypMm0acPSoMX0WqctXLuPVD1/Fg1UPjvoaTs1dHzNNRBRMSUmJpRv4MtOkNxZNJptfMh9H
   Lh4JaUuV8eOB9HTg448jv97aA2vx2WmfHTUADjAETkQUKTum6EhfYy5uSeFLTkjGorJF2HR8E75w
   0xeu+/pbbgFqaoxNdAsKjDWb8vOvfh34OStr+BSbPwD+5J1Pjnr+nh7gwAFg7lwz3p17cZ0mIgom
   JSUFeXl5uHz5MtLT02N+Pa7TpDcWTTHg31IllKLprbeAc+eApibg/Pmrnw8dArZsGX58YGB4EYXi
   HTibewWH3lyMy5OHF1mZmcb5d+8GZswwRrSIiCh8ZWVlOH/+vCVFE+mNRVMMLKtchp/95WdQSo0a
   zvZLTgbKy42P6+noMAoofxH184/WYE7Xahw4IPB6rx5vagISEq6uNH7nnVG/JdfjKBMRjSY3Nxci
   gsHBQXg8sU21cJRJbyyaYmB6znQAwOGLh0PeUiUUmZnGR2Ul0Nrdiq89uR6H/9c/IS9j+OuUAtrb
   rxZRM8xrAhFR3ElOTkZBQQGam5uRkZFx/W8g12IQPAYi2VIlXGsPrMXyyuXIy8gLcn1g3Djj7ryF
   C4FJk2LWDNfgOk1ENJbS0lL09fXF/Dpcp0lvLJpiJNwtVcKhlMJTDU9h9ZzgK4ATEZG5cnNz4fF4
   MDAwYHdTyEYsmmIk3C1VwvHeuffQ2duJJRVLTD93vGKmiYjGkpSUhJtuugltbW0xHXFipklvLJpi
   JNwtVcKxpmENvjHnG/AI//mIiKxSWVmJefPmoaury/KNfEkP/K0bQ+FsqRKqtp42vPzBy3io6iFT
   zxvvmGkiolAUFRVh4cKFGBwcREdHh+nnZ6ZJbyyaYiicLVVCtXb/WiydshT5mfmmnpeIiEIzceJE
   1NbWIjU1FW1tbVwtPI6waIqhcLZUCQUD4LHDTBMRhSMjIwMLFy7ExIkT0draalrhxEyT3lg0xVDg
   lipm2HVuF9p62lA3pc6U8xERUeRSUlIwf/58lJaW4vLly7yzLg6waIox/5YqZmAAPHaYaSKiSCQk
   JGD27Nmm3VnHTJPe+Ns3xpZVGmHwaIdu23ra8NIHL+Hh2Q+b1DIiIjKDiGD69OmoqalBZ2cnenp6
   7G4SxQiLphgL3FIlGs8feB51FXUoyCwwo1k0AjNNRBSt4uJiLFy4EP39/ejs7IzoHMw06Y1FU4yZ
   taXKmt1rsLqaAXAiIp3l5OSgtrYWycnJvLPOhVg0WSDaLVUazjXg0pVLWDplqYmtokDMNBGRWTIz
   M7Fo0SJMmDAh7DvrmGnSG4smC0S7pQoD4EREzuK/s66kpAQtLS0YHBy0u0lkAv4WtkA0W6q097Rj
   3aF1eLiKAfBYYqaJiMyWmJiIOXPmYMaMGWhtbUV/f/91v4eZJr2xaLJIpFuq/Pv7/44l5UswOWty
   DFpFRESxJCK48cYbUV1djY6ODt5Z53AsmiwS6ZYqDIBbg5kmIoql0tJSLFiwAH19fWPeWcdMk95Y
   NFkkki1Vdn+8Gxc6L2DZlGUxbBkREVkhNzcXtbW1SEpK4p11DsWiySKRbKnydMPT+PqcryPBkxDD
   lhHATBMRWSMrKwuLFi1CdnZ20DvrmGnSG4smC4WzpUpHbwdeOPgCA+BERC6TmpqK22+/HUVFRbyz
   zmFYNFkonC1VXnj/BSwqW4SicUUWtIyYaSIiKyUmJqK6uhrTp09HS0vL0J11zDTpjUWThcLZUuWp
   hqcYACcicjGPx8M76xwm0e4GxBP/lireY17MyJ0x6uv2fLwHTR1N+EzlZyxsXXxjpomI7CAiKCsr
   Q1paGnbs2IGKigq0tLTA4/EgISFh6LP/axGxu8lh88+uKKWGPkY+Dve5UP87+F8X+D1jfX9fX9+Y
   52PRZLHllcvx3P7n8N1bvzvqa57ezQA4EVE8ycvLw5IlS9DS0oLe3l50d3fjypUr6OnpQU9PD7q7
   u9Hb2xv0F/7g4OCYRVaoBYZSCoODg0OfA78O/BxO0eIvdEQEHo9nWDs9Hg8SExOvaXvgewh8L/6v
   R74n/9fBjkX6/Ggk1rc8iojibZVXXey6iCm/moIL//MCkhOSr3m+s7cTJb8owf6/3o/iccU2tDA+
   bd26laNNRGS7sfoipRT6+/vR19c37MNfZAV++IutgYGBoMVA4KiL/7OIIDExEUlJSdd8HvmRmJh4
   TWETWOD4C7bAr51CRKCUCtpgjjRZLHBLldry2muef+HgC1hQuoAFExERDSMiQ0VLqAYGBoYVWAMD
   A0PFzcjCx+NhzPl6ONJkg8c2PQYAeLzu8Wueu+1fbsMPF/0Qd99wt9XNIiIiintjjTSxrLTBaFuq
   7Gvah7PtZ3Hn1DttaBURERGNhUWTDUbbUuXp3U/j67O/jkQPZ02txnWaiEgH7Iv0xqLJBsG2VOnq
   68Lz7z+Pr87+qo0tIyIiotGwaLLJyC1V1h1ch9tLbkfJ+BIbWxW/eOccEemAfZHeWDTZZOSWKmsa
   1mD1HK4ATkREpCsWTTYJ3FLlwPkDONV6Cp+d9lmbWxW/mCMgIh2wL9IbE8c2CdxS5aNLH+Frs7/G
   ADgREZHGuE6TjdYdXIenGp7C3qa92PPNPSgdX2p3k4iIiOIa12nSVF1FHbYc34Lbim9jwURERKQ5
   Fk02yknPwZKKJfhOzXfsbkrcY46AiHTAvkhvDNHYzLvSC4+wdiUiItIdM01EREREPsw0EREREUWJ
   RRMRmCMgIj2wL9IbiyYiIiKiEDDTREREROTDTBMRERFRlFg0EYE5AiLSA/sivbFoIiIiIgoBM01E
   REREPsw0EREREUWJRRMRmCMgIj2wL9IbiyYiIiKiEDDTREREROTDTBMRERFRlFg0EYE5AiLSA/si
   vbFoIiIiIgoBM01EREREPsw0EREREUWJRRMRmCMgIj2wL9JbxEWTiPxERPaJyF4R2SQiJWY2jMhK
   e/futbsJRETsizQXzUjTE0qpWUqpKgCvAviRSW0islxLS4vdTSAiYl+kuYiLJqVUe8DDTADN0Tcn
   fLEeyrRiqNTp78EN/41OnDgR0/Pz34Dnd8I1eH77r8G+SI9rjCaqTJOIPC4ipwA8COCfzGlSePgD
   4P7zW3GNWA+J89+A53fCNXh++6/BvkiPa4xmzCUHRMQLoCDIU48qpV4PeN33AUxXSj0c5Bxcb4CI
   iIgcY7QlB0xZp0lESgG8qZS6OeqTEREREWkomrvnpgU8vBfAnuibQ0RERKSniEeaROQlANMBDAA4
   BuCvlVKfmNg2IiIiIm3EfBsVIiIiIjdw3IrgItJxnee3iki1Ve3RhYgUi8hrInJERI6KyC9FJGmM
   1/83EUmzso06uN7PT7wTkc+JyKCITLe7LU7A/uha7ItCw75obLr2RY4rmgBcb2hMhfAaVxERAfAK
   gFeUUjcAuAHG2lmPj/FtjwBIt6B5uomrn40IfBHAf/g+h0xEnNiXmIH9UQD2RWGJm5+LCGnZFzmy
   oxORWhEJXPLgNyLyoJ1tstkdAK4opf4AAEqpQQD/HcBXRSRdRH4uIgd82958R0S+C6AQwBYR2WRj
   u20hIhkiUi8iDSKyX0Tu8R0vF5EPRGSNiLwvIn8SkVS722sVEckEcCuA7wD4r75ji0XkbRH5DxH5
   UET+n+8XI0Skw/eztRfAbfa13F7sj4ZhXxQG9kXB6dwXObJoCiKu/poL4iYADYEHfCu2nwLwdQBl
   AGYppWYBWKuU+jWAcwAWK6XqrG6sBq4AWKGUqobRyf9zwHNTAfzGt3xGC4DP29A+u9wLYINS6hSA
   CyIyx3e8BkbnNRNAJYD/4jueDuBdpVSVUuovlrdWX/HcH7EvCg/7ouC07YvcUjTFu9E6aAGwGMBv
   fX/xQSl12apGacwD4P+IyD4AXgCFIpLne+64Umq/7+sGAOU2tM8uXwTwou/rF32PFYCdSqkTvp+h
   5wEs8L1mAMDLlreSdMa+KDzsi4LTti9KtOIiMdCP4QVf3IUIRzgE4L7AAyIyDkAJgEYYHRZd9WUA
   uQDmKKUGROQ4AP/Qd0/A6wYQJz9bIjIRwBIAN/tW8U+A0Um9geG/CAXAoO/rbsXbbwH2R4HYF4WH
   fdEIuvdFTh1pOglgpogki0g2jGHNuKWU2gQgXURWAoCIJMAY5n0GwEYA3/Qdg4hM8H1bO4BxNjRX
   B+MBfOLrpJbAmDKId/cBeFYpVa6UqlBKlQI4DmARgHm+jIUHRr7gz3Y2VEPsj3zYF4WNfdG1tO6L
   HFU0iUgigB6l1BkA6wC8D+AFALttbZgeVgC4X0SOADgMoAvAowD+BUaeYL8vJOe/E2ENgA3xFL70
   //wAWAtgrojsB7ASwAcBLxv510q8jKQ8AGD9iGMv+46/B+A3MEYRjiml/K+Ll/82QbE/GhX7outg
   XzQmrfsiRy1uKSKzADyllIrbO3Uocvz5CZ+I1AL4nlLqr+xui27480SR4s9O+HTpixwz0iQi3wLw
   bwB+aHdbyHn48xMV5/xlZRH+PFGk+LMTFdv7IkeNNBERERHZxTEjTUThEJESEdkiIgd9i8P9re/4
   RBHxirHFw0ZfcNf/PT8QkY98C6ctDzhe7VuQ7yMRedKO90NEzmRyX/S4iJwSkXY73gtxpIlcSkQK
   ABQopfb6VpdtAPA5AA8DaFZKPSEi/xvABKXU90VkJowh8xoARQDqAUxTSikR2QngO0qpnSLyJoBf
   KaU22PLGiMhRTO6L5sEI03+klMqy5Q3FOY40kSsppZqUUnt9X3fAuCulCMA9AP7ge9kfYHRegLEC
   7fNKqT6l1AkARwHcKiKTAWQppXb6XvdswPcQEY3JrL7I9/07lVJNFjafRmDRRK4nIuUAZgPYASBf
   KXXe99R5APm+rwsBnAn4tjMwOraRx8/6jhMRhSXKvog0wKKJXM03HP4ygEd8e2AN8a0gy/lpIoq5
   KPsi9lOaYNFEriUiSTA6qeeUUq/6Dp/3ZQzgm3r7xHf8LIytHvyKYfyFd9b3deDxs7FsNxG5iwl9
   EfscTbBoIlcSEQHwOwCHlFK/DHjqjwAe9H39IIBXA44/4NsKowLANBibQzYBaBORW33nXBnwPURE
   YzKrL7KqvTQ23j1HriQiCwC8DWA/rg5t/wBG57MOQCmAEwC+oJRq8X3PowC+CmMD1keUUn/yHa8G
   8HsYG2a+qZT6W8veCBE5msl90RMwtp+ZDOBjAE8rpf7RsjdDLJqIiIiIQsHpOSIiIqIQsGgiIiIi
   CgGLJiIiIqIQsGgiIiIiCgGLJiIiIqIQsGgiIiIiCgGLJiIiIqIQsGgiIiIiCsH/B0y1VWs7mE+f
   AAAAAElFTkSuQmCC
   "
   >
   </div>
   
   </div>
   
   </div>
   </div>
   
   </div>

   <script src="https://c328740.ssl.cf1.rackcdn.com/mathjax/latest/MathJax.js?config=TeX-AMS_HTML"type="text/javascript"></script>
   <script type="text/javascript">
   init_mathjax = function() {
       if (window.MathJax) {
           // MathJax loaded
           MathJax.Hub.Config({
               tex2jax: {
               // I'm not sure about the \( and \[ below. It messes with the
               // prompt, and I think it's an issue with the template. -SS
                   inlineMath: [ ['$','$'], ["\\(","\\)"] ],
                   displayMath: [ ['$$','$$'], ["\\[","\\]"] ]
               },
               displayAlign: 'left', // Change this to 'center' to center equations.
               "HTML-CSS": {
                   styles: {'.MathJax_Display': {"margin": 0}}
               }
           });
           MathJax.Hub.Queue(["Typeset",MathJax.Hub]);
       }
   }
   init_mathjax();

   // since we have to load this in a ..raw:: directive we will add the css
   // after the fact
   function loadcssfile(filename){
       var fileref=document.createElement("link")
       fileref.setAttribute("rel", "stylesheet")
       fileref.setAttribute("type", "text/css")
       fileref.setAttribute("href", filename)

       document.getElementsByTagName("head")[0].appendChild(fileref)
   }
   // loadcssfile({{pathto("_static/nbviewer.pygments.css", 1) }})
   // loadcssfile({{pathto("_static/nbviewer.min.css", 1) }})
   loadcssfile("../../../_static/nbviewer.pygments.css")
   loadcssfile("../../../_static/ipython.min.css")
   </script>