Quantile regression
===================


.. _quantile_regression_notebook:

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

.. raw:: html

   
   <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>This example page shows how to use <code>statsmodels</code>&#39; <code>QuantReg</code> class to replicate parts of the analysis published in </p>
   <ul>
   <li>Koenker, Roger and Kevin F. Hallock. &quot;Quantile Regressioin&quot;. Journal of Economic Perspectives, Volume 15, Number 4, Fall 2001, Pages 143–156</li>
   </ul>
   <p>We are interested in the relationship between income and expenditures on food for a sample of working class Belgian households in 1857 (the Engel data). </p>
   <h2 id="setup">Setup</h2>
   <p>We first need to load some modules and to retrieve the data. Conveniently, the Engel dataset is shipped with <code>statsmodels</code>.</p>
   </div>
   </div>
   </div>
   <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">patsy</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">pandas</span> <span class="kn">as</span> <span class="nn">pd</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">statsmodels.formula.api</span> <span class="kn">as</span> <span class="nn">smf</span>
   <span class="kn">import</span> <span class="nn">matplotlib.pyplot</span> <span class="kn">as</span> <span class="nn">plt</span>
   <span class="kn">from</span> <span class="nn">statsmodels.regression.quantile_regression</span> <span class="kn">import</span> <span class="n">QuantReg</span>
   
   <span class="n">data</span> <span class="o">=</span> <span class="n">sm</span><span class="o">.</span><span class="n">datasets</span><span class="o">.</span><span class="n">engel</span><span class="o">.</span><span class="n">load_pandas</span><span class="p">()</span><span class="o">.</span><span class="n">data</span>
   <span class="n">data</span><span class="o">.</span><span class="n">head</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[1]:</div>
   
   <div class="output_html rendered_html output_subarea output_pyout">
   <div style="max-height:1000px;max-width:1500px;overflow:auto;">
   <table border="1" class="dataframe">
     <thead>
       <tr style="text-align: right;">
         <th></th>
         <th>income</th>
         <th>foodexp</th>
       </tr>
     </thead>
     <tbody>
       <tr>
         <th>0</th>
         <td> 420.157651</td>
         <td> 255.839425</td>
       </tr>
       <tr>
         <th>1</th>
         <td> 541.411707</td>
         <td> 310.958667</td>
       </tr>
       <tr>
         <th>2</th>
         <td> 901.157457</td>
         <td> 485.680014</td>
       </tr>
       <tr>
         <th>3</th>
         <td> 639.080229</td>
         <td> 402.997356</td>
       </tr>
       <tr>
         <th>4</th>
         <td> 750.875606</td>
         <td> 495.560775</td>
       </tr>
     </tbody>
   </table>
   </div>
   </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">
   <h2 id="least-absolute-deviation">Least Absolute Deviation</h2>
   <p>The LAD model is a special case of quantile regression where q=0.5</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">mod</span> <span class="o">=</span> <span class="n">smf</span><span class="o">.</span><span class="n">quantreg</span><span class="p">(</span><span class="s">&#39;foodexp ~ income&#39;</span><span class="p">,</span> <span class="n">data</span><span class="p">)</span>
   <span class="n">res</span> <span class="o">=</span> <span class="n">mod</span><span class="o">.</span><span class="n">fit</span><span class="p">(</span><span class="n">q</span><span class="o">=.</span><span class="mi">5</span><span class="p">)</span>
   <span class="k">print</span><span class="p">(</span><span class="n">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>
                            QuantReg Regression Results                          
   ==============================================================================
   Dep. Variable:                foodexp   Pseudo R-squared:               0.6206
   Model:                       QuantReg   Bandwidth:                       64.51
   Method:                 Least Squares   Sparsity:                        209.3
   Date:                Thu, 21 May 2015   No. Observations:                  235
   Time:                        05:58:25   Df Residuals:                      233
                                           Df Model:                            1
   ==============================================================================
                    coef    std err          t      P&gt;|t|      [95.0% Conf. Int.]
   ------------------------------------------------------------------------------
   Intercept     81.4823     14.634      5.568      0.000        52.649   110.315
   income         0.5602      0.013     42.516      0.000         0.534     0.586
   ==============================================================================
   
   The condition number is large, 2.38e+03. This might indicate that there are
   strong multicollinearity or other numerical problems.
   
   </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">
   <h2 id="visualizing-the-results">Visualizing the results</h2>
   <p>We estimate the quantile regression model for many quantiles between .05 and .95, and compare best fit line from each of these models to Ordinary Least Squares results. </p>
   </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">
   <h3 id="prepare-data-for-plotting">Prepare data for plotting</h3>
   <p>For convenience, we place the quantile regression results in a Pandas DataFrame, and the OLS results in a dictionary.</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">quantiles</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="o">.</span><span class="mo">05</span><span class="p">,</span> <span class="o">.</span><span class="mi">96</span><span class="p">,</span> <span class="o">.</span><span class="mi">1</span><span class="p">)</span>
   <span class="k">def</span> <span class="nf">fit_model</span><span class="p">(</span><span class="n">q</span><span class="p">):</span>
       <span class="n">res</span> <span class="o">=</span> <span class="n">mod</span><span class="o">.</span><span class="n">fit</span><span class="p">(</span><span class="n">q</span><span class="o">=</span><span class="n">q</span><span class="p">)</span>
       <span class="k">return</span> <span class="p">[</span><span class="n">q</span><span class="p">,</span> <span class="n">res</span><span class="o">.</span><span class="n">params</span><span class="p">[</span><span class="s">&#39;Intercept&#39;</span><span class="p">],</span> <span class="n">res</span><span class="o">.</span><span class="n">params</span><span class="p">[</span><span class="s">&#39;income&#39;</span><span class="p">]]</span> <span class="o">+</span> \
               <span class="n">res</span><span class="o">.</span><span class="n">conf_int</span><span class="p">()</span><span class="o">.</span><span class="n">ix</span><span class="p">[</span><span class="s">&#39;income&#39;</span><span class="p">]</span><span class="o">.</span><span class="n">tolist</span><span class="p">()</span>
       
   <span class="n">models</span> <span class="o">=</span> <span class="p">[</span><span class="n">fit_model</span><span class="p">(</span><span class="n">x</span><span class="p">)</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">quantiles</span><span class="p">]</span>
   <span class="n">models</span> <span class="o">=</span> <span class="n">pd</span><span class="o">.</span><span class="n">DataFrame</span><span class="p">(</span><span class="n">models</span><span class="p">,</span> <span class="n">columns</span><span class="o">=</span><span class="p">[</span><span class="s">&#39;q&#39;</span><span class="p">,</span> <span class="s">&#39;a&#39;</span><span class="p">,</span> <span class="s">&#39;b&#39;</span><span class="p">,</span><span class="s">&#39;lb&#39;</span><span class="p">,</span><span class="s">&#39;ub&#39;</span><span class="p">])</span>
   
   <span class="n">ols</span> <span class="o">=</span> <span class="n">smf</span><span class="o">.</span><span class="n">ols</span><span class="p">(</span><span class="s">&#39;foodexp ~ income&#39;</span><span class="p">,</span> <span class="n">data</span><span class="p">)</span><span class="o">.</span><span class="n">fit</span><span class="p">()</span>
   <span class="n">ols_ci</span> <span class="o">=</span> <span class="n">ols</span><span class="o">.</span><span class="n">conf_int</span><span class="p">()</span><span class="o">.</span><span class="n">ix</span><span class="p">[</span><span class="s">&#39;income&#39;</span><span class="p">]</span><span class="o">.</span><span class="n">tolist</span><span class="p">()</span>
   <span class="n">ols</span> <span class="o">=</span> <span class="nb">dict</span><span class="p">(</span><span class="n">a</span> <span class="o">=</span> <span class="n">ols</span><span class="o">.</span><span class="n">params</span><span class="p">[</span><span class="s">&#39;Intercept&#39;</span><span class="p">],</span>
              <span class="n">b</span> <span class="o">=</span> <span class="n">ols</span><span class="o">.</span><span class="n">params</span><span class="p">[</span><span class="s">&#39;income&#39;</span><span class="p">],</span>
              <span class="n">lb</span> <span class="o">=</span> <span class="n">ols_ci</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span>
              <span class="n">ub</span> <span class="o">=</span> <span class="n">ols_ci</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
   
   <span class="k">print</span><span class="p">(</span><span class="n">models</span><span class="p">)</span>
   <span class="k">print</span><span class="p">(</span><span class="n">ols</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>
         q           a         b        lb        ub
   0  0.05  124.880100  0.343361  0.268632  0.418090
   1  0.15  111.693660  0.423708  0.382780  0.464636
   2  0.25   95.483539  0.474103  0.439900  0.508306
   3  0.35  105.841294  0.488901  0.457759  0.520043
   4  0.45   81.083647  0.552428  0.525021  0.579835
   5  0.55   89.661370  0.565601  0.540955  0.590247
   6  0.65   74.033435  0.604576  0.582169  0.626982
   7  0.75   62.396584  0.644014  0.622411  0.665617
   8  0.85   52.272216  0.677603  0.657383  0.697823
   9  0.95   64.103964  0.709069  0.687831  0.730306
   {&apos;a&apos;: 147.47538852370633, &apos;b&apos;: 0.48517842367692343, &apos;lb&apos;: 0.45687381301842317, &apos;ub&apos;: 0.51348303433542364}
   
   </pre>
   </div>
   </div>
   
   <div class="output_area"><div class="prompt"></div>
   <div class="output_subarea output_stream output_stderr output_text">
   <pre>
   /build/buildd/statsmodels-0.6.1/debian/python-statsmodels/usr/lib/python2.7/dist-packages/statsmodels/regression/quantile_regression.py:189: ConvergenceWarning: Convergence cycle detected
     warnings.warn(&quot;Convergence cycle detected&quot;, ConvergenceWarning)
   /build/buildd/statsmodels-0.6.1/debian/python-statsmodels/usr/lib/python2.7/dist-packages/statsmodels/regression/quantile_regression.py:189: ConvergenceWarning: Convergence cycle detected
     warnings.warn(&quot;Convergence cycle detected&quot;, ConvergenceWarning)
   
   </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">
   <h3 id="first-plot">First plot</h3>
   <p>This plot compares best fit lines for 10 quantile regression models to the least squares fit. As Koenker and Hallock (2001) point out, we see that:</p>
   <ol>
   <li>Food expenditure increases with income</li>
   <li>The <em>dispersion</em> of food expenditure increases with income</li>
   <li>The least squares estimates fit low income observations quite poorly (i.e. the OLS line passes over most low income households)</li>
   </ol>
   </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">x</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">income</span><span class="o">.</span><span class="n">min</span><span class="p">(),</span> <span class="n">data</span><span class="o">.</span><span class="n">income</span><span class="o">.</span><span class="n">max</span><span class="p">(),</span> <span class="mi">50</span><span class="p">)</span>
   <span class="n">get_y</span> <span class="o">=</span> <span class="k">lambda</span> <span class="n">a</span><span class="p">,</span> <span class="n">b</span><span class="p">:</span> <span class="n">a</span> <span class="o">+</span> <span class="n">b</span> <span class="o">*</span> <span class="n">x</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">8</span><span class="p">,</span> <span class="mi">6</span><span class="p">))</span>
   
   <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">shape</span><span class="p">[</span><span class="mi">0</span><span class="p">]):</span>
       <span class="n">y</span> <span class="o">=</span> <span class="n">get_y</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">a</span><span class="p">[</span><span class="n">i</span><span class="p">],</span> <span class="n">models</span><span class="o">.</span><span class="n">b</span><span class="p">[</span><span class="n">i</span><span class="p">])</span>
       <span class="n">ax</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">x</span><span class="p">,</span> <span class="n">y</span><span class="p">,</span> <span class="n">linestyle</span><span class="o">=</span><span class="s">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;grey&#39;</span><span class="p">)</span>
       
   <span class="n">y</span> <span class="o">=</span> <span class="n">get_y</span><span class="p">(</span><span class="n">ols</span><span class="p">[</span><span class="s">&#39;a&#39;</span><span class="p">],</span> <span class="n">ols</span><span class="p">[</span><span class="s">&#39;b&#39;</span><span class="p">])</span>
   
   <span class="n">ax</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">x</span><span class="p">,</span> <span class="n">y</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;red&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s">&#39;OLS&#39;</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">scatter</span><span class="p">(</span><span class="n">data</span><span class="o">.</span><span class="n">income</span><span class="p">,</span> <span class="n">data</span><span class="o">.</span><span class="n">foodexp</span><span class="p">,</span> <span class="n">alpha</span><span class="o">=.</span><span class="mi">2</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_xlim</span><span class="p">((</span><span class="mi">240</span><span class="p">,</span> <span class="mi">3000</span><span class="p">))</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_ylim</span><span class="p">((</span><span class="mi">240</span><span class="p">,</span> <span class="mi">2000</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">ax</span><span class="o">.</span><span class="n">set_xlabel</span><span class="p">(</span><span class="s">&#39;Income&#39;</span><span class="p">,</span> <span class="n">fontsize</span><span class="o">=</span><span class="mi">16</span><span class="p">)</span>
   <span class="n">ax</span><span class="o">.</span><span class="n">set_ylabel</span><span class="p">(</span><span class="s">&#39;Food expenditure&#39;</span><span class="p">,</span> <span class="n">fontsize</span><span class="o">=</span><span class="mi">16</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,iVBORw0KGgoAAAANSUhEUgAAAggAAAFWCAYAAAAMkGV5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd4VdeVt98tCUkIVZBABTCiVxmDMWBjWtwgcXAh4PrZ
   sRNPJpmU+ZLMxE6xJ/OlTeLMJJlxih0b3GBsjDG2MdWWAdOLAYEoEhJqSEK93qa7vz/O1fGVULnl
   qIDX+zz38d3nnrP3vkLWXnvttX5Laa0RBEEQBEHwJqSvJyAIgiAIQv9DDARBEARBEC5DDARBEARB
   EC5DDARBEARBEC5DDARBEARBEC5DDARBEARBEC6jVw0EpdQIpdRHSqmTSqkspdR3PNcHK6W2KaXO
   KqW2KqXivZ55Uil1Til1Wil1m9f1mUqpE57P/tCb30MQBEEQrnZ624PgBP5Zaz0FmAN8Syk1CfgR
   sE1rPR7Y4WmjlJoMrAQmA3cAzymllKevPwOPa63HAeOUUnf07lcRBEEQhKuXXjUQtNalWutPPe8b
   gGwgDfgysNpz22rgLs/7ZcAarbVTa50P5ACzlVIpQIzW+oDnvpe9nhEEQRAEIUj6LAZBKTUKuA7Y
   DwzTWpd5PioDhnnepwJFXo8VYRgU7a8Xe64LgiAIgmABYX0xqFIqGngL+K7Wuv6zUwPQWmullCX6
   z1b1IwiCIAhXElpr1f1dXdPrHgSl1AAM4+AVrfUGz+UypVSy5/MUoNxzvRgY4fX4cAzPQbHnvff1
   4o7G01r3+9fTTz/d53O40l/yM5SfYX94yc/w6voZut1uLly40CP9/vWvf6WmpqZH5m0VvZ3FoIC/
   A6e01v/l9dFG4BHP+0eADV7X71NKhSul0oFxwAGtdSlQp5Sa7enzYa9nBEEQBCEgamtrycvLA0Ap
   xciRIy0fQynFo48+SlxcnCX9NTU1sXr1alpaWizpr5Xe9iDcBDwELFJKHfW87gB+DdyqlDoLLPa0
   0VqfAt4ATgEfAN/Un5lH3wReAM4BOVrrzb37VQRBEISrjYaGBsrLy7u/0U+qq6tZv369ucMPDw+3
   rO+oqCiWLl1KaGioZX1CL8cgaK1307lRcksnz/wS+GUH1w8D06ybXd+xcOHCvp7CFY/8DINHfobB
   Iz/D4OmLn6HD4SAkJISwsDDS0tJIS7M+5j0uLo5Zs2bhHXMXDOfPn6e2tpbrrrsOgKSkJEv69UZZ
   eV7R31BK6av5+wmCIAjBs3nzZlJTU8nIyLC035KSEux2O+np6Zb2C1BZWUlzczPDhw+/7DOlFNqC
   IMU+yWIQBEEQhECxahd+NdCTm2AxEARBEIQrDvEO97yhJMWaBEEQBEG4DDEQBEEQBEG4DDEQBEEQ
   BEG4DDEQBEEQBEG4DDEQBEEQBEG4DDEQBEEQBMFiVq1axbRp0xg0aBApKSl885vfpLa2FoBnnnmG
   hx9+uMPndu/ezY033kh8fDxDhgxh3rx5HDp0qDenbiIGgiAIgiBYyLPPPsuPfvQjnn32Werq6ti3
   bx8XLlzg1ltvxel0dvpcXV0dX/rSl/jud79LdXU1xcXFPP3000RERPTi7D9DlBQFQRCEKwqPUmBf
   T6ND6urqSEtL46WXXmL58uXm9cbGRtLT0/nNb35DQUEBOTk5vPLKK22ePXToELfeeivV1dU+jdXZ
   z8EqJUXxIAiCIAiCRezZswebzcY999zT5vqgQYNYunQp27Zt6/TZCRMmEBoayqOPPsrmzZt9NhR6
   CjEQBEEQhKsPpYJ/BUBFRQWJiYmEhFy+vKakpFBRUdHpszExMezevRulFF//+tcZOnQoy5Yt65Hq
   kr4gBoIgCIJw9aF18K8ASExMpKKiArfbfdlnJSUlJCYmdvn8xIkTeemllygsLCQrK4uSkhK+973v
   +TR2R2MGgxgIgiAIgmARc+fOJSIigrfeeqvN9YaGBjZv3swtt9zic18TJkzgkUceISsrq9t7q6qq
   WLVqlaWxGWIgCIIgCIJFxMXF8fTTT/Ptb3+bLVu24HQ6yc/PZ8WKFYwYMYKHHnoIrTVutxu73Y7N
   ZsNms2G32zlz5gy///3vKS4uBqCwsJA1a9Ywd+7cbscdPHgwK1assLSAkxgIgiAIQr+nuLiYjz76
   qK+n4RM//OEP+eUvf8kPfvAD4uLimDNnDtdccw07duwgPDwcpRRr1qxh4MCBREVFERUVxbhx44iJ
   iWH//v3Mnj2b6Oho5s6dS0ZGBs8++2ynYx09etR8Hx0dben3kDRHQRAEod9js9m4dOkSI0aM6Ndp
   jr2JUgq73U54ePhl161IcxQDQRAEQeiXlJeXExYWxuDBg9tcFwPBQHQQBEEQhM8lhYWFXLp0qa+n
   8blFPAiCIAhCv8HhcFzmMm+PeBAMxIMgCIIgfG54/fXX+0wYSGiLeBAEQRCEfoN4EHxHPAiCIAjC
   VYvNZuOdd97B5XIBdGscCL1HrxsISqkXlVJlSqkTXtduUEodUEodVUodVErN8vrsSaXUOaXUaaXU
   bV7XZyqlTng++0Nvfw9BEAQheCIiIpg4cSKhoaF+PaeU+ty/epq+8CC8BNzR7tp/AD/VWl8H/MzT
   Rik1GVgJTPY885z67KfyZ+BxrfU4YJxSqn2fgiAIQj/EZrNRUFAAGAv9hAkT/FrwtNYdvo4ePUpV
   VVWnn/vzKisro7q62pK+unvV19eTmZkZ0LM9Sa8bCFrrXUD7GpYXgTjP+3ig2PN+GbBGa+3UWucD
   OcBspVQKEKO1PuC572Xgrh6duCAIgmAJ1dXVnD171vJ+Q0JCaGlpsaSvvLy8Hg2WLCkpoba2FjAU
   EBcsWNBjYwVKWF9PwMOPgN1Kqd9hGC2twtOpwD6v+4qANMDped9Ksee6IAiC0A9xu9243W7CwsJI
   SUkhJSXFkn4vXrxo9pWRkRFUX7W1tcTFGXvV2bNnBz23rigsLCQpKckcrz/SXwyEvwPf0Vq/rZT6
   CvAicKsVHT/zzDPm+4ULF7Jw4UIruhUEQRD8YNeuXURGRlq68DqdTrZv386KFSuIiIgIqi+3280b
   b7zBAw88wKBBgyya4WdorSkpKSEtzdjLWvlzyMzMJDMz07L+WumTNEel1CjgXa31NE+7Tmsd63mv
   gBqtdZxS6kcAWutfez7bDDwNXAA+0lpP8ly/H1igtf5Gu3EkzVEQBKEf4HA4CAsLIyQk+JNtl8tF
   WFiYeQYfTMBeS0uLGSCpte6x4D+Hw8Gbb77JihUrGDBgQI+M0crVluaYo5RqPYBZDLQeTm0E7lNK
   hSul0oFxwAGtdSlQp5Sa7TEoHgY29PqsBUEQhE7ZunUrVVVVgJG+aIVxUFBQwFtvvQUQdDR/YWEh
   b775ptm22jhwuVzU19cDxvd/8MEHe9Q4OHjwIKWlpZb11+seBKXUGmABkAiUYWQtnAD+B4gAmoFv
   aq2Peu5/CngMcAHf1Vpv8VyfCawCBgKbtNbf6WAs8SAIgiD0Ebm5uQwfPjxo9783WmvsdjuRkZGW
   9GWz2Rg4cKAFM7uco0ePUldX12sBiDk5OQwZMoTBgwdLNcfuEANBEASh99Bak5eXx+jRoy3t9+DB
   g0RHRzNp0qSg+9q1axdDhw5lwoQJFszscpqbm02DoyePLADq6+s5dOgQixYtanP9ajtiEARBEK5w
   nE4nR44cwel0WtrvyJEjGT58uCV9TZw4kWuuucaSvtqjtWbNmjVUVlYC1h9ZtCcqKoqEhIQe00MQ
   D4IgCIIQFK1Bg1ahtebYsWNMnTo16H5dLhcff/wxCxcu9Fut0Z8xWufpHfTYE3z66ackJCR0aeSI
   B0EQBEHocwoLC1m3bp3l/VZXV2Oz2YLuJzQ0lLi4ONxutwWzupySkhLWrl3bZryeJD4+nujo6B4d
   oxXxIAiCIAgBY1XQoNaampoaEhISgp5TY2Mj5eXlpKenB91XR3inV1oZNNkRzc3N7Nq1i1tvvdXn
   IwvxIAiCIAh9wsmTJzl16hRgLEZWLI41NTW8++67lpyn19fXU1hYGHQ/nbF582ZOnz4NWPf9OyMy
   MpJhw4b1SXlr8SAIgiAIflFaWkpoaChJSUlB9+Ud6R9M1H9dXR2RkZG9Ui66rq6O6OhoS3QdOuLk
   yZMMGDCA8ePH+/xMU1MTERERhIaGigdBEARB6D2KiopwuVwAJCcnW2IcHD16tI1EcDBR/7t27TIr
   RFpNc3Mzr776qlkIKjY2tseMAzDiDOLj4/16ZtOmTRQVFXV/ox+IB0EQBEHolk2bNjFz5kyGDRtm
   WZ82mw2lVMBCSk6n01Qm7GnNgZKSElJTU3ukb4fDwfbt27n99tt9DnLUWlNbW2saEm632zRaxIMg
   CIIg9CjepZOXLl1qiXGQmZlpyi9HRkYGbBw4HA7+9re/4XA4AOs1B7Kzszl8+LDZ7injAGDAgAGM
   GDHCr2cuXbrE+++/b7Z7wqMhHgRBEAThMhwOBy+88AJf+9rXLD3XP3fuHKmpqZZUTLTb7ZbKOHtT
   VVWF2+0mMTGxR/o/e/YsDoeDqVOn+vxMbW0tUVFR3XpNxIMgCIIg9Bjh4eF89atfDdo4cLvdnD17
   1myPGzcuYOPg/PnzbN682WxbXeNh27ZtpvbC4MGDe8w4AIiLi2PIkCF+PbNz506Ki4vNdk8rNYqB
   IAiCIABQUVHBzp07zbYVRYxcLhfZ2dlmgGMwpKWlMXv27KD76QilFMnJyT3SNxjHNe+9954pQz1s
   2DBSUlK6fMbtdlNWVma2v/SlLzFq1Kgem2N7xEAQBEEQAIiJibEsdbGpqQkwPBHLli0LWDJ5z549
   lJeXA4bHwAohpVaqq6s5fvy42Z42bVqPaRqEhoYyevRov3b9NTU1ZGZmthFm6k3EQBAEQfgcU19f
   T3V1NWAswFZUTDx37hzbtm0Luh+ApKSkHlu0Q0JCsNvtPdI3QF5eXptAx8mTJ3drKF26dInm5mbA
   OOZYuXJlrxsGrUiQoiAIwueYI0eOADBjxgzL+tRao7UOKLLe5XJx6tQpMjIyLJuPNzk5OSQnJ/dK
   PYOqqiqampr8qkS5Y8cOxowZE9RRggQpCoIgCAHhnb44Y8YMS4yDHTt2kJ2dDRgLVDBpdyUlJZbE
   LHRERUUF9fX1PdK31pr333+/jQegO+PA5XKRn59vtr/whS/0apxBV4iBIAiC8Dlj3bp1lqsOzpgx
   g7Fjxwb0rMPhMLURwsLCuOOOOywrH+12u9sswHPmzOk2ODBQlFKMHTvWr4qOzc3NfPrpp31Sa6E7
   xEAQBEH4nHHnnXf6LczTHrfbTWZmphmVn5CQYObn+0tubi7Hjh0Laj6dYbPZOHToUI+Vey4sLGTv
   3r1me8KECd2mhhYVFVFXVwcYgaF33XVXn8UZdIUYCIIgCFc5LpeLLVu2mG77qKiooBckpRSDBg0K
   eOG12+3mrnnSpEksWrQoqPl409zcbC7AUVFRLF++vMdqJ8THx/utslhYWGgGhvZnxEAQBEG4ygkN
   DSU5OTloo6ClpYWSkhLAMBBmzZoVsFjR22+/3WMlmY8fP95GnMlqNm/e3MYDcM0113R5v91uN+Mz
   AObOndvtM/0BMRAEQRCuQlpaWrh48SJgLObXXnutX2fjHVFVVcXBgwcDft77nP3ee+9l5MiRQc3H
   G+/Aw9mzZ3P99ddb1nd7xo4d65dhpLUmPz+/X8YZdIUYCIIgCFchlZWV7N+/P+h+3G63eTSRlJTE
   smXLAuqnsbGRl156yTySCDReobM5rlmzhoaGBsv69Ka0tLRNWWpfDITz589TUVEBGEWplixZ0i/j
   DLpCdBAEQRCuEoLRH+iM3bt3ExISwo033hh0X5WVlX7XH+gKh8NhBgT2ZLlnm81GSUkJo0eP9vmZ
   Y8eOMXjw4KCDQQPBKh0EMRAEQRCuEg4ePEhjYyMLFy60rE+n00loaGhARkdubi41NTXMnDnTsvm0
   kpOTw7Fjx7j33nst7xtg+/btTJ8+3eeCTTabjaysrB492vCVK1YoSSn1olKqTCl1ot31byulspVS
   WUqp33hdf1IpdU4pdVopdZvX9ZlKqROez/7Qm99BEAShP5KRkcFNN90UdD/vvPMOlZWVgHEUEKhH
   IiEhwVLNAbfbbZ7jjxkzJuDjDl8YM2aMX1UnQ0NDaWhouOLiDLqiL2IQXgLu8L6glFoEfBnI0FpP
   BX7nuT4ZWAlM9jzznPrMh/Rn4HGt9ThgnFKqTZ+CIAifB/bu3Wsu5hEREZac7c+cOZP4+PiAnv30
   00/bKAn6mwLYFW+99RYXLlwAjF2yVWJKYBx/bN261Wynp6d3W83y9OnTZlbHgAEDWLhw4RUXZ9AV
   vW4gaK13Ae0TQP8R+JXW2um555Ln+jJgjdbaqbXOB3KA2UqpFCBGa33Ac9/LwF09PnlBEIR+RjAC
   Ra3Y7fY2RYWGDx8ecMaDzWbDZrMFNR9vvHfkS5cu7bH0wNjYWL9iDIAe01YICgslqvvLtxsHzFdK
   7VNKZSqlWg9xUoEir/uKgLQOrhd7rguCIFz1tKYvAkycOJHY2Nig+lNKUVdXF5Dokdvtpqjosz/H
   c+bMsawkc21tLa+88oppJAwaNMjSHfrOnTvNuQ8YMKBbqejm5uY25ZfHjx9vqYckKFwuWL0aJk+2
   rEvr/DPBEQYkaK3nKKVmAW8A/plynfDMM8+Y7xcuXGhp8I4gCEJv43K52LFjB/fee2+3LvCucDqd
   NDY2Eh8fT3h4eMBKho2NjXzyySesWLHCcvd6bGwsX/ziF3vMbZ+enu6XMRMREUFUVFSPZkz4jdNJ
   5lNPkfnCCxAXBwsWwLlzlnTdJ1kMSqlRwLta62me9gfAr7XWH3vaOcAc4GsAWutfe65vBp4GLgAf
   aa0nea7fDyzQWn+j3TiSxSAIwlWB2+221KV98uRJysrKWLx4cUBzcTgcREZGWjafVg4fPkxYWBjX
   Xnut5X3X1taya9cuv4yOrKwsoqKi/D5+6HEcDli1Cn71KxgzBn72M5g/H7iCsxg6YQOwGEApNR4I
   11pXABuB+5RS4UqpdIyjiANa61KgTik12xO0+LCnD0EQhKuO8vJyXnvttaD78c4CmDx5MosXL6au
   ro7s7PNkZ5835YO74+DBg5aIMHXEqFGjGDNmTI/0HRMTw4QJE/x6JjY2lujo6B6ZT0DY7fDcczB2
   LKxfD6+9Btu3m8aBlfS6B0EptQZYAAwByoGfAa8CLwLTAQfwfa11puf+p4DHABfwXa31Fs/1mcAq
   YCCwSWv9nQ7GEg+CIAhXPFprGhoaiImJCaqfDRs2MHXqVPOsva6ujr178wkLM0K4XK5i5s4d1WFM
   g7db3e12o5SyxM3e0tLCpk2buP3227utghgI+/btIzEx0edS1DabjR07drBkyZL+FYRos8Hzz8Nv
   fgPTp8NPfwqzZ3d4qwgl+YAYCIIgXKnk5+fT1NTEZAuDzhobG9tUcszOPk9JSRzx8Ya6YU1NJamp
   tUyadLk7/ZVXXuGOO+4gKSnJsvm0kp2dzbhx4yxNW2zl4sWLREdH+2xcaa05deoUkyZN6h8GQlMT
   /O1v8NvfwvXXG0cJHQhPud1uCgoKGDVq1FV3xCAIgiB4ERkZSVRUVFB9OJ1ONm7ciNPpBILLAli2
   bJnPqoLdUVpayvHjx832pEmTLDMOmpqaePvtt82MjJSUlG6Ng2PHjnHy5EnA2H1PmTKl742DxkZ4
   9lkjvmDnTnjvPXjnnQ6NAzCCV/fv32/WzbACMRAEQRD6CZWVleZinpyczKhRo4Lqb8CAAYwfP75T
   TYO0tERcrmJqaiqpqanE5SomLc0wAmpqanj33XfNmIXY2FjLIvfDw8MtLdbkzcCBA5k2bZpfc01O
   Tu4/6YoNDcYxwpgxsH8/bNlixBpcd91lt2ZnZ1NaWgoYP9OVK1da6oURA0EQBKGfsG/fPlOZL1Ca
   mprIzc012xMnTux0NxwbG8vcuaNITa0lNbW2TfxBbGwskydPtswoOHbsGE1NTYChsDhp0iRL+gU4
   cuRIGw/A2LFju5y33W7njTfeMHfbw4YNs0y7IWDq6uCXv4TRo+HTT2HHDnjjDcjI6PQRz1FCj01J
   YhAEQRD6EKtz6isqKsjKygpI8yU/Px8gaM9FR+zbt4/x48czePBgy/suLy8nMjLSL8GovLw887y+
   T6mpgT/9Cf74R7jjDvjxj2HixA5vrays5JNPPuHLX/5yl11KDIIgCMIVTktLC88//7y5sw4Uu92O
   3W4HIDExMWBBuNZy0VbgcDjIyckx23PmzLHMOHA4HKxbt870AAwdOrRb4+DIkSMcOnTIbKenp/et
   cVBVBU8/baQr5uTAJ5/AK69cZhx4p6bGx8czffr0XpuiGAiCIAh9RGhoKCtXrgw6GHH37t1kZ2cH
   9GxeXp4Z0Jeenk56enpQc2nFbrdz5syZHnGBh4eHM2PGDL/qRaSnpzOxk515r1JRYXgJxo2DoiLY
   t8+QSB4/vsPbX3/9dVNaOzQ0lJEjR/baVH0+YlBKhQB3AvOBwcAzWusLSqmFwFmtdXAHZz2AHDEI
   gtDfqK+vJysri7lz5wbVT3tdgkCi7rXWbNy4kQULFgRcvdGbmpoaQkNDg9Zr6IgTJ07gcDiY2UkU
   f3ucTievvfYaDzzwQI/oK/jNpUtGVsLf/gbLl8OTT0IHxlhLSwsNDQ3ExcUBRv0HfyW1e/WIQSmV
   AOwB3ga+DvwfDKEjMOSQfxTsRARBEK5kfFUkbI3gD2bzorXmlVdeoaamBvCvqqDWmvr6esBYSJYt
   W2aJcQCGLLF34SYrSUtL80nuuPXnOmDAAJYsWdL3xkFZGfzgBzBhghGIePSoYSR04qk5e/Ys+/bt
   M9vB1NsIFp88CEqpF4A7gBXAAQy1w+u11keUUo8C/6K1tk7NwyLEgyAIQm/QnSKhw+GgqanJsoUY
   oKqqioSEBL/P0S9dusTmzZt5+OGHfbq/rq6O4uIKwEiLbH/WX1VV1SOBhy0tLbz11lssW7aMiIgI
   n545fPgwzc3NzJs3z/L5+E1JCfzHfxjHBw89BP/6rzB8+GW3ud1usrKyzNRMK4JWeztIcRnwE631
   ng4+KwRGBDsRQRCEK5Xi4grCwtKIjx9CfPwQwsLSzEUVjHx1b2GgQKivr+fjjz8224MHD/Z5IdFa
   m3EGSUlJPPTQQz4912r4lJTEUVISx969+W28I06nk/Xr15sBklYSGhrKDTfc0K1egvcmcOLEicya
   NcvyufhFURF8+9swdSooBSdPGlkKHRgHYCzmxcXF2Gw2s+0v9fX1vP7665bHe/hqIEQDnfmNIoF+
   UvdSEAShf+D9xzojI4P5QRbTaU3jC2QR2LVrFwcPHjTbvi5CHRk+RUWXaG5uBgw3/uOPP+7zDr87
   zpw5w549n+1DR40a1eXxicvl4q9//atpoAwaNMiyufhNQQF885uGbkFkJJw6Bf/5n9CBAFN2djZn
   z54FjH+LJUuW+H2U4HA4TKMvOjqahQsXWp6V4auBcBa4vZPP5gMnrJmOIAjClUdHioSnTx/j3Llz
   QGC7QjB2hpcuXQKMxfi6664LqK9Zs2Zx/fXXBzSH9hQU5LXxZFi5KKWkpPhUbbF1YQwLC+OBBx7o
   O6MAID8f/uEfDKXD2Fg4fdqom5Cc3OkjsbGxfmk2dMSGDRsoKCgAjH+DnlCC9DUG4Qngv4GfA68D
   OcCtwDWe609orV+1fHZBIjEIgiD0Fu3P6sHY2QWj6Z+VlYXNZvN7cddas27dOpYsWRJUqeLWI4aQ
   kGGEhIThchUze/ZI4uLiLDEMApnnoUOHqK+vZ9GiRUGPHxS5uYby4YYN8I//CN/7HnRSq6KhoYF1
   69bxyCOPBPxza2pqorq6mrQ0I86lpaWl0zTPXq/mqJT6NfAD2nod3MBvtNY/DnYiPYEYCIIg9BZu
   t5tPPvmEOXPmBFVnwG63Ex4eHvQCXFhYSFpaWtBFh+rq6nj11deZMWMWkyePC3rn256CggKGDx/e
   5Ty9F0O73U5oaGiPVH70iXPn4Be/MIonfetb8N3vQgdBmk6nk5CQEHPeFRUVQRW7KiwsJDc31ycR
   rF5XUtRa/wgYC/wD8FPgm8D4/mocCIIg9CZKKcLDw033d6CsX78+oFTB6upqDhw4YLZHjBgRlHHg
   XaTp8ce/ypw5My0xDvLy8vjoo4/M9siRI7s1Dv7yl7+YcQ8RERF9YxycPg0PPww33mjUS8jJgX/7
   tw6NA4CNGzeSl5dntv01DtxuN9u2bTPVIkeMGBGwQmagdOtBUEpFAKXAI1rrjb0yK4sQD4IgCD2J
   1pqqqiqGDBnS/c0+4nK5AloAm5qaOH36NDNmzAh6DuXl5WzdutXnbAd/aGpqorGxkaSkpE7v0Vrj
   crlMT4zNZiMyMtLyufjEqVPw7/9uFE/63vfgn/7JiDVoh8PhoLKykpSUFKDrIwBfOXr0KJMnT/Y7
   xqLXPAhaazvgAmzBDiYIgnA1UVNTwwcffBBUepnNZuO1114zd4r+GAelpaVUV1cDEBUVZYlxAEYq
   ZHcFgXxFa8369evbzLMr4wCMugnegZB9YhwcPw4rVsCiRTB9uhFz8NRTHRoHYBwhHD161GwHYhzs
   37+f/fv3m+3rrruuTwMwfQ1SfB7QWusnen5K1iEeBEEQegJvMRsrhG1KSkoCikI/ePAgCQkJjB07
   NqjxATIzMxk6dCiTJ1uveVdUVERKSkqXi6a3l8DlchESEhJ0/ERAfPop/PznsGePoYD4jW9ABwGU
   Wmv279/PjBkzglJrbGxsZNCgQYCRtRIeHh60UdDbMQibgKVKqbeUUg8ppb6glFrs/Qp2IoIgCFcC
   p06dYvv27WY7EOOgtSRzK/4YB607cTDSF60wDsDQarCqr6KiIjZv3my2hw8f3qVx4Ha7efHFF82q
   lmFhYb1vHBw6BMuWwdKlcPPNcP68YSB0kl3R+u/ucDgCHrK2tpbXXnvN9EDFxMT0bcpmO3z1IHQX
   daO11sEdtvQA4kEQBMFq7HY7LpfL3PUFQmVlJYWFhX6X7nW73bz00kvcf//9QVeAtNvtvP/++3z5
   y1+2POjVujbGAAAgAElEQVTP4XBQVVVFchdaAG63G5vNZn4PK87sA+LAASPY8NgxQw75a1+DTkSL
   zpw5Q01NDbNnzw54uNzcXFJTU01hJCs8UO3p1TRHT8XGLtFaZwY7GasRA0EQBCs4ceIEKSkpQaWp
   1dTUEBEREVDxHe/ARSsXlJycHMaMGWNJf++++y6zZs3q0ijw5ujRo1RWVnLLLbcEPXZA7NljHCWc
   OgU/+hE89pihgNgO70qZNTU1OJ3ObmMouuLjjz9m4sSJDBs2LOA+uqPXdRCuRMRAEATBCk6ePElS
   UhJDhw4NuI/t27czcuRIxo8f79dzBQUF7Nmzh/vuuy/gsVvJy8ujvr6ejIyMoPtqz8WLF0lKSmrj
   jehIPKo1VdLtdqOUsnz33C27dhkeg5wcI+jwkUegE7e+w+Hg+eef54knnghY26K4uJiioqKgvA7+
   IgaCD4iBIAhCoFhRpTDQlEVvtNY4HA5LzqYvXTLqKIwcOTLovsrLy9m7dy/Lli3r8PP2FS6dziLy
   8/fy4IMPEhMTE/T4fqE1fPyxYRgUFMCPf2xoGnSw6LcWTWoNmPQOIvR9uM+8PLW1tVRWVvpUqtoq
   ejVIUSn1kVLqw05eHymlPgx2IoIgCP0FrTXvvPMO9fX1Affhdrt5/vnnaWxs9PvZzMxMTp8+DRh/
   7AM1DrTW7N271yxmlJSUZIlxADBkyBBuuOGGTj8vLq4gJCSZiIgw4uOHMGDAcObPv7V3jQOtYft2
   WLAAvv51ePRRQ/Doscc6NA4Adu7cSW5urtn21zhwu9387W9/M4Wd4uLietU4sBJfw0SV5xXi9UoC
   5gHj8aOao1LqRaVUmVLqsgJPSqnvK6XcSqnBXteeVEqdU0qdVkrd5nV9plLqhOezP/g6viAIQme0
   ehyVUjz66KNBLWYhISE8+uijAQUzTps2jfT09IDHbkUpRWhoKE6nM+i+ALZs2UJ+fj5g5Pm3igJ1
   RmVlLhcuHGozn15Ba9iyBebNM4SNnngCsrON44R2hoHNZiMnJ8ds33rrrUyZMsWv4ex2u5mBERIS
   wv333x9QrEl/wycDQWu9UGu9yPPf1tc0YBJQBfzCjzFfAu5of1EpNQKjANQFr2uTgZXAZM8zz6nP
   fsP+DDyutR4HjFNKXdanIAiCr9TW1vLqq6+2MRL8paKigg8++MBs+7pIuFwutm7daoolDRkyJGCv
   QWNjI2fOnDHbN9xwQ1AFm7y57rrrzGJBnVFWVobWmrS0ROLjo0hOnm5WuGyNQ+gxtIb334c5c+D/
   /l/49rfh5El46CHo5KjH6XSaVTchsH/3vXv3tvmZW12voq8IKtFUa50L/Br4rR/P7AKqO/jo98C/
   tLu2DFijtXZqrfMxqkjOVkqlADFa61bh8ZeBu/ycviAIgklsbCxLliwJapcbHx/PxIkT/X4uNDQ0
   qABIb5xOJ6WlpZb0VV1dzZtvvmkaTUOHDu0yWE9rzfbt26mrqyM2NpYbb0wnNbWW1NRa5s4d1XML
   p9awcSPMmgVPPgk//CGcOAH33QcdpE7u3r2b2tpawNAeWLJkiV/D2Ww2Tp48abYXLFjAddddF9x3
   6IdYkfxaAXRfwLsLlFLLgCKt9fF2/3OmAvu82kVAGuD0vG+l2HNdEATBZ8rKyqiurmbixIkopQJK
   YywpKQEMsaOwsDCfjwbq6+spLy830wz91UTw5uLFi8TGxjJo0CDi4+NZsGBBwH15Ex8fz0033dSl
   0WS326mrqyMpKQmlFA8++KD5WWxsbM/upt1ueOcdI10R4Gc/M8SOuhFZGjJkSNBCTMXFxUyePLlv
   MjF6iaAMBKVUIvDPQG5393bRRxTwFMbxgnk5mHl588wzz5jvFy5c2OvVsARB6L+0FgUKhrq6uoAW
   m6amJkpLSxkzZkxQ4wOcO3eOa665JijxplYyMzNJTk42jabuVB4LCgooLi7u3b+tbje89ZZRRCk8
   3MhOuPNO6GShPn/+PGfPnuWOO4yT6EmTJvk95Pbt27n22mtJSkoiMjKS2267rfuHeonMzEwyMzMt
   79dXoaQ8QNN24Q4HhnmuL9dav+PzoEqNAt7VWk9TSk0DtgNNno+HY3gEZgNfBdBa/9rz3GbgaYw4
   hY+01pM81+8HFmitv9FuHElzFAShDY2NjYSHhwec1w6GYE5cXJzfO8eamhqioqKC0u4Hw7ApLy/v
   EbGdioqKbiV/8/Pzuy3T3CO0tMCbbxqGQXQ0PP00LFnSoWHgcDjMn7PNZsNmsxEfH+/XcN4iSfn5
   +QwbNuyKCD7s7VoMHwM7Pf9tfb0L/BSY6I9x0B6t9Qmt9TCtdbrWOh3j6GCG1roM2Ajcp5QKV0ql
   A+OAA1rrUqBOKTXbE7T4MLAh0DkIgvD5oX0aWyBs2rSJS5cu+f3c3r17KSwsDGpsMBa8LVu20NLS
   EnRfDQ0NvPbaa7jdhqJ+YmJil8aB1ppPP/2Uurq6oMf2GZcLXn0VpkyBP/4Rfv972LfPqJvQgXHQ
   mmLamlkQGRnpt3Fw5swZ3n33XbM9atSofm8ctH5fq+h1oSSl1BpgATAEKAd+prV+yevz88D1Wusq
   T/sp4DGMktPf1Vpv8VyfCawCBgKbtNbf6WAs8SAIgtCGQKWKA63g6L2TDQaXy4XD4Qi6BkN7tNaU
   lpZ2mbLY2NhIdXU1w4cPt3TsbnG54LXX4Be/gORkw2OweHGHRkFjYyNOp9M0BAIRqSovLzeDRZ1O
   J1prS/7teoPS0lK2bdvGww8/3OtCSS96dvAdfXaNUupFXwfUWt+vtU7VWkdorUd4Gweez0e3Ggee
   9i+11mO11hNbjQPP9cNa62mezy4zDgRBEMBYAF999VUzaj0Q46ChoYEXXnjB3GX7Yxw8//zzlugQ
   HD58mCNHjgTdDxiejE8//RQwvkv3egaV5OXlWTK2Tzid8Pe/w4QJsGoV/O1vhhLiF77QaZzByZMn
   28zRX+OgpaWFTZs2mUqKAwYM6NfGgdaa3bt3mzE0w4YNaxMgagX+VHOc45VW6P3Z9Rhu/z4o3N01
   4kEQBAEMieHExMSgos1ra2uJi4vz6V5vD4PT6Qw43sFut5vufiuLNFVXVzNw4EBTTrg9Wmuys7MZ
   P3685ZUeu8Ruh9Wr4Ve/grFj4ac/hfnzO7nVTlZWFjNnzgx4uLy8PCIiIvwqt93XeP8e7Nu3j4yM
   jMu8Sr0dg9AVw4BmC/oRBEGwBJvNxqFDnyn4tabg+cOFCxc4fPiw2fbVOMjJyWHz5s1mO1DjQGvN
   qlWrgvJ8tGK321m9erW520xISOjUOGgdq6SkhIaGhoDH9AubDZ57DsaNg/XrjWOFbds6NQ7A0I6o
   qakxvTqBYLfbLVOZ7A0yMzM5ePCg2Z4zZ47lR07edOpBUErdDdyNkbnwILAZQ/PAmyjgZuCM1rrz
   f8k+QjwIgvD5xOFwsHfvXubPnx/wwlpdXU1tbS2jRo3ye+zm5mafDYr2eEfOW1HsqZXS0tIuSzHX
   1tZSXl7OuHHjLBnPJ5qb4fnn4T/+A6ZPNzwGXVQ93L17NyNHjgy4nkRNTQ2ZmZksW7bsitAuaG5u
   pqyszPwdbG5uJiIiotvskd7wIFwDzMcwAACme957v6YCnwBfD3YigiAIweB2u83iSuHh4SxYsMDv
   RaCgoMCMBE9ISPDZONi1axdlZWXm2IEaBxcuXGDdunVmOxjj4PDhw+zb95nOXFfGARi76crKyoDH
   84umJvjP/4QxY4xiShs2wHvvdWgceG/yrrnmGr8rbLpcLrOP2NjYoASpepumpibOnj1rtgcOHNir
   qaW+xiDkA3dprT/t8RlZiHgQBOHzw+nTpzl//jxLly4NuI+PP/6YMWPG+B2tf/78eYYOHRpQzQPv
   M2W3243dbrckna6uro6wsLAuXdBHjhxh8uTJXR43WEpjI/z5z/C738FNN8FPfgJdSBQXFxezc+dO
   7r///oCHfPXVV1m0aFG3NST6A1pr1q5dy1133RXU74BVHoReT3PsTcRAEITPD4EWWWpqavL7HNfp
   dHLixAmuu+66oF3V7733HpMnTw66JLDL5WLVqlXMmjWP8PBw0tISu5U53rt3L5MnTw7Y4+Ez9fXw
   P/9jeA0WLDAMg4yMTm6tJzo6GqUUbreb5uZmvxQibTYb9fX1JCUlAdalmfYUzc3NaK3N38GioiJS
   UlII7aCGhK/0+BGDUmq+UirG632Xr2AnIgiC4C+ZmZlm0ZxANPFbF1W73e7Xc0opKisrLREquumm
   m/yOc+iIpqYmhgyZyKVLiZSUxLF3b/5lYkZVVVUcP37cbM+dO7dnjYO6OkPDYMwYOH4cPvwQ3nij
   U+MA4K233qKqysh0DwkJ8Vs+uqCggFOnTpnt/mwcAOzZs4fz58+b7eHDhwdlHFhJV0GKZmqj531X
   aK11//hGXogHQRCubqqqqoiOjg5qEWhpafHpD7LdbqehoYEhQ4YEPBYYmgrvv/8+X/nKV4I+T87K
   yqKiooKFCxeSnX2ekpI44uON+dXUVJKaWsukSZ95Jmpra8nNzWXGjBlBjdstNTWG4uGf/gR33AE/
   /jF0UuWytraWpqYmU4vB33TOlpYWPvnkE+bNm9f70s8BUFFRQU5ODnPmzOmxMazyIHQVAbMYyPZ6
   LwiC0KdorTly5AgZGRkMGDDA74A1MM61Dx48yF13GRXiuzMO6urqKC6uoKioAHBy6623dnl/dwwa
   NIgbb7zRksUsPT292+qRe/bsISMjg+joaOLi4nrWOKiqgj/8wThOuPNO2LPHSF3sgvLyckpLS6mp
   MbLlfTka8SYkJISwsDBcLle/9RZ4a2EMHDiQmJiYPp6Rb3T6G6q1ztRa13u97/LVazMWBOFzTVNT
   k99HAt4kJyczb948n+4tLy9nz548SkriCAmZhts9LKAaBCdPniQrKwswdncjRozwuw8wghhXr15t
   6hMMGjTIdMGnpSXichVTU1NJTU0lLlcxaWmJREdHB6UV4BMVFYaXYNw4KC6G/fvhpZc6NA6cTicf
   fvihOadhw4bR1BRHSUlcp0cj7dm7d695VKKU4sYbb+y3xoHb7eYvf/mLmR0zaNAgpkyZ0sez8g0J
   UhQEod9TX18f1K7r9OnTREVF+Z0///zzfycxcTYjR04FOnbb+0JZWRmhoaEkJib69VxHdKUKWVdX
   R3Z2LqWlJSxadLNfO/EAJwPPPmtoGSxfDk8+Cd3EU2itOXToENOnT2fAgAE+HY1AW1XJqqoqBg4c
   2G+LJ+Xl5REdHW0GSgajphkIvRGk+JFS6kMfXh8ppT4MdiKCIAgd0dDQwNq1a4PaBQ8YMMBnTQHv
   TcXcufOJje26TkFHuN1uPvzwQ1Olb9iwYQEbB+fOnWujzNiVKmRsbCzTpk1gzJhretY4KCuDH/zA
   qJVQVwdHj8Jf/9qpcdB+xz9r1iy/FsyKigpeeeUVsz148OB+axyAYag1N38mMNybxoGVdBWkmAlo
   DCVFgPFAMpCPUYVxGIaYUimGkuKiHp6r34gHQRCuDrzVBX2lqKiItLQ0vwLeGhsbWbNmDY899hgh
   ISHU1dWxd28+YWFGDr3LVczcuaN8WnwPHz7M1KlTuyyd3BnegXo2mw2n09mlByUzM5OMjIyAYjL8
   oqQEfvtbo17CQw/Bv/wLdKIZ4a0C2dWOv7OfcXV1Nampqebi6mswaV9QWFjIoUOHuPvuu/t6KkAv
   eBC01gu11ou01guBPwBOYK6n2uIcrXU6MBdwAP8V7EQEQRBaycvLY8sWs3ir38aB1pp9+/b5HS8w
   aNAgli9fbo4XGxvL3LmjSE2tJTW1tkvjoKqqitOnT5vtmTNnBmwcvP7662aqX2RkZLfHK2lpaQGN
   5TNFRfDtb8NU46iFkyeNLIVOjIOqqipeeukl0xvT1Y6/s5/xqVOnqKmpMe/rT8aB1prCwkKznZyc
   zOLFV18sv69KilnAb7XWqzv47FHg+1rradZPLzjEgyAIVyaBphQGctZ77tw5ampqmDVrll/Ptaei
   ooL8/Hyuv/76gJ739hpUVlYyePDgTr0fly5d4vjx43zhC18IeL4+UVAAv/41rF0Ljz8O3/8+dCLZ
   XFFRQXx8vOk18FegqKioiKqqKjK60EjoL7jdbtauXcs999zTeyqUftDb1RzHYhwrdMQloBerewiC
   cDWSn5/PpUuXAIiIiPDbOKiurmbVqlX4uylITEwMuPhPTk4ONpvN7CdQ4+DChQts3LjRbA8ZMqTL
   o5G4uDiuueaagMbyifx8eOIJo4BSbCycPm0cLXRRz+GTTz4x61GA/wJFkZGRPVqZMFj2799PXl4e
   YHi0HnjggX5pHFiJrwZCPvCN9heV8Rv8hOdzQRCEgKmvrzdTwQIhISGBRx55xKeYg0OHDrUpyjRs
   2LCAxiwuLjbLMfuLtyGTmprKwoULu7x/27ZtXLx4ETAW37FjxwY0bpfk5hqegpkzYehQOHvW8CAM
   HXrZrdXV1W0KCS1btsyvegcul4sXXngBh8MBGAZWj3ynANFa09jYaLaHDx8etEjWlYavpcKeAV73
   HDWsA8owghS/AkzEKActCILgF42NjWYe/7Rp/p9Snjp1itraWubOnQv4vmt1u904HA6/d6xut5uS
   khKzmNOCBQuAz8SUwHehn3Xr1jFv3jxSUlIYMGBAt5LHkyZN6rkgxLNn4Ze/NCoqfutbcO4cdDOW
   0+n02zhyOp243W4iIiIICwvjrrvu6rf6BefPn+fYsWPcc889AFdEsSer8VkHQSl1C/BvwPXAAIyg
   xYPA01rrHT02wyCQGARB6L9orVm9ejV33nlnwDuz2tpa3G43CQkJXd7ndrspKioK+Cihlbq6OrZs
   2cLy5ctNT4U/mQ7ecQZVVVUkJCR0GWewZ88eli1bFtScu+T0afh//w+2bIHvfMcIRIyP7/BWrTXv
   vvsut99+e8ABkZs3byY1NbVfxhm0tLSwZ88e5s2b13qGD/hf/Ks/0GfVHJVSoUAiUKG1Dr5SSQ8i
   BoIg9G8CSV88fvw4Y8aM8auIT0NDA5s3b+bee+/1+w9++2p77fFV6OfixYvs2rWLFStW+DRuS0sL
   paWlPbNzPXnSMAx27IDvfQ/+6Z+MWIMO8DZqzpw5w+jRo30OBLXZbBQVFZlHB4H8e/ck3kaA1po9
   e/Ywa9asfuvV8JXeDlI00Vq3aK3L+rtxIAhC/6OpqYk1a9aYVRADWSxsNpsZGNgVLS0tplhNdHR0
   m12/Pxw4cKDNWbs/eG9Qhg0bxm233dbl/R988AH5+fmAkdZnuXFw/DisWAGLFxsBiLm58NRTnRoH
   hw4dYvfu3WZ7woQJfmWJOBwOzp07Z7b7k3EAsGnTJrPyo1KKm2666Yo3DqzEnyOGMcAKYARwWeim
   1voxa6cWPOJBEIT+hdaa4uJi8wzfVyoqKvxWIty3bx9Op5Obb77Zr+fAMGRaPQbdVRfs6ojhnXfe
   ISMjo8uCSt79X7p0iYSEBJ9VH33m00/h5z+HvXsNBcRvfAM68cB4x4U0NTURFhbm16K5Y8cOZs2a
   1fMyzwHQ0NBAVVWVedTU1NTEwIEDr8hjhK7oVQ+CUuou4DTw78BdwCKv12LPfwVBEC6jrq7O3EUq
   pfw2Dmw2Gxs3bsTlcnV7r7cc8w033OBzUSZvWlpaWL16tel96G7xaC/0M2fOZzLHCxYs6DIdsaKi
   grVr15rtpKQka42DQ4dg2TJYuhTmzzc8Bt//fqfGQXNzM6tXrzZ/jlFRUT4ZB94bseHDh1tv4FhE
   bW0tFy5cMNtRUVFXnXFgJb4KJZ0ALgIPaq0v9fisLEI8CILQ95SXl5Obm2tmGviK93l1d7v41nte
   fvllli5dahbJ8QdvkaVAz8orKyvZtGkTDz30UKfzbX/uXVVVZX363P79hsfg2DH413+Fr30NOlEy
   vHjxIjExMURHRwP+f/eTJ09SWFjIHXfcYcnUrcTlcrF27VpWrlx5xdZDCIReDVJUSjUC92itt3R7
   c/d9vQh8EShvVV9USv0W+BKGbHMu8FWtda3nsyeBx4AW4Dta662e6zOBVRjHHZu01t/tYCwxEASh
   D3A6nWitAz7PzcnJISsri7vuusuv5wKt+pidnc25c+f48pe/7Pez7Rf82tpa4jvJBAB47733mDBh
   AuM6KIUcNHv2GIbBqVNGZcWvfhW6EfPZtWsXI0aMYFQ3VRi9qampMb9ja+ntHpV69oP6+nrCw8PN
   +QRSk+NKp7eDFM8AVpm4LwHtTc2twBSt9bXAWeBJAKXUZGAlMNnzzHPqs3/lPwOPa63HAeOUUv3P
   fBWEzyk7d+4kKysr4OfT09O7DegDY6HasGGDuUj7Yxx4bx7Gjx/P0qVL/Z8osHXrVrMGg1KqQ+PA
   ++hj/vz5jBkzJqCxOmXXLrjlFnjgAbjnHsjJgX/8xw6Ng5qaGg4dOmS2b775Zr+MA6fTyf/+7/+a
   lSojIiL6jXEAxu9eUVGR2R4+fPjnyjiwEl8NhH8BnvIEKgaF1noXUN3u2jatdev/QfuB1kPKZcAa
   rbVTa50P5ACzlVIpQIzW+oDnvpcxYiMEQegHLFy4kBkzZvj1zIEDB9pE8PsiYhQbG8v06dMDWgDW
   rl1LSUmJOZ4/5+ZtS0LPZfz48Z3eW1VV1UYCOjY21ppofq3ho49g0SJ49FHDODh3zpBI7sJzEx4e
   7vfPKz8/n4oKQwhqwIABPPHEE/3GZV9SUtLG4PniF79ovQH2OcXX39KngcHAKaVUllJqp9drl1Jq
   p4VzegzY5HmfChR5fVYEpHVwvdhzXRCEPmLDhg1UVlYCgVXeS0lJ8Ukp8Pz582008f3Z/XrzxS9+
   kZSUFL+fq6+v58UXX2yz4Lf/vm632/QaJCQksGLFCut2sVrD9u2wYIFhDDz6qCF49Nhj0Mmi/fbb
   b5vVIaOiopg5c6ZfQ9bU1LSRHe7rHXmr9wKM79PVkY4QOL6azC0Yxwyd/VZYctCvlPox4NBav25F
   fwDPPPOM+X7hwoXd6p0LghAYs2bN8usPtdaaU6dOMWnSJEJCQhgxYoRPz4WEhAS0A6+qqmLbtm3m
   Yu1vGl5roGRMTEy3gkvvv/8+Y8eOZdKkSSilzADAoNAatm6Ff/s3qKqCn/4UVq6ETjwfLS0tpuEy
   e/Zsv75vZWUlhw8fNo95pk+fHvz8LcLpdPLcc8/xrW99i7CwMOLj4z/3BkJmZiaZmZmW9+u3kqIl
   gyo1CnjXu0S0p2z014EvaK1tnms/AtBa/9rT3ozhzbgAfKS1nuS5fj+wQGvdpqCUBCkKQs/R0tJC
   dnY2U6dODeh5t9vNli1bmD9/fpeqiFprcnJyGDNmTFCuea015eXlARVmyszMJDY2tstjE+8siObm
   ZiIjI63ZaWsNmzYZwYcNDYZh8JWvQBdemhMnTlBYWOhXXIV39oLT6SQ/P79nAikD4MyZMyQnJ5v1
   KgIp6/15os+UFHsCT4DhD4FlrcaBh43AfUqpcKVUOkZZ6QNa61KgTik12xO0+DCwodcnLgifY1pa
   WigsLDRVEX1Ba019fT1geAKWLFnik2TymTNnzOf84cCBA20CCP0xDrw3F9dffz3XXnttp/fW1NS0
   OXawRHxHa9i4EWbNMjISfvhDOHEC7ruvQ+OgNUYAYOLEidx+++1+Dffyyy9TXl4OGHEG/cU4AONY
   p1WXAhDjoJfwR0lxOPB9YD5GPMKdWusspdQ/A3u01vt97GcNsACjnkMZhkfgSSAcqPLctldr/U3P
   /U9hxCW4gO+2plp6pTkOxEhz/E4HY4kHQRD8oLuqhFprbDYbAzvJqe+OsrIytm3bxkMPPdTlfVpr
   6urquq1w6Mt4UVFRfqc+2mw2Vq9ezde+9rVO4ylahZtagxsdDoc1Mr1uN2zYAP/+70b7Zz8zxI66
   8J64XC5efvllHnzwQZ8zChwOB01NTaZ73ls9sq85d+4c586dCziz5PNOb+sgTAF2YcQi7MPQMbhe
   a31EKfVfwFCt9QPBTsZqxEAQBN/xpSphXl4ehw8fZvny5X717S105IsQT3l5Odu2bePBB/2rJO90
   Otm+fTu33XZbQIGS3vOsra3t0kB57733SE9PZ8qUKX6P0yFuN7z1lmEYhIcbhsGdd0InnojCwkIi
   IiIYOnRoQMMdP36cmpoa5s+fH8ysLcHtdlNcXGzGodhsNlpaWvwqyCV8Rm8bCJuBGAwtgmYMQaNW
   A2EF8Butdedi432EGAiC4Du+VCXUWuN2u/1afA8cOEBLS0u3Sort+/ZFPbGjPk6cOMGUKVP8NhD2
   7t2L2+3mpptu6vQe7122y+WyRlK4pQXefNMwDKKj4emnYcmSTg2DVk6cOEFUVJTPKX1Op5OjR48y
   a9asPs9CaI/T6eTNN99k5cqVARl2Qlt6OwZhHoYR0NEhYBmQHOxEBEHonxw+fJjjx48Dxh8ef/+A
   T5482SdNhI8//pjDhw+b7Y4Wsbq6OrKzz5OdfZ66ujrAyIPPzs42n8nIyPB5jt4biGuvvZbZs2d3
   em99fT2vvPKK+UzQxoHLBa++ClOmwB//CL//PezbZ9RN6OC719fXs23bNrM9bdo0v/L9Q0NDaWho
   8CtmpCfJzMw0dSgGDBjAAw88IMaBBVi5Kfb1N9xN56mMiRheBUEQrmDS0hIpKMinpsZou1zFpKWN
   wuGI7fRcu7OYhZ07dzJjxgyio6N9TvGbM2dOl2f47Y9ACgrymTt3FCEhIQHtiF0uFy+88AJf/epX
   iYiI6PD83VtGOCYmhq9//evB776dTnjtNfjFLyA5Gf77v+ELX+jQKPCWcY6KiiI5Odkvz8onn3xC
   UlIS48ePJyQkhMWLFwc39yBwu93YbDbz5zxmzJigY0yEtrzzzjtMnjzZsv58PWLYAdRpre9WSoXR
   9ohhLTBIa32nZbOyCDliEAT/aF3wi4sLue66qV0WEeoqZiErK4v09PQuz5Ddbjf/+7//y5133umT
   EUqbZhEAACAASURBVOF9BFJcfILIyERGjGhucwTiC94LbHe1GzZv3szw4cMDTuVsg8MBr7wCv/wl
   jBxpxBgsXNjlUcKGDRuYPn26X2JQ3kcfZWVlREdH94uz/OPHj1NSUtIvizpdqeTm5uJ0Opk4cSJg
   lLMeNGgQISEhvXrE8HPgTqXUNoyUQoBblFIvA/cAvwh2IoIg9D2xsbFMmjSakJCWbnepxcUVhIWl
   ER8/hLi4BBoatOlNmDp1areLUkhISLcaCJ3R1FSN02nr/sZ2HD16lA8//NBsd2QcVFd/pgR/++23
   B28c2O3w17/C+PGwdi289NJnEsntfsZa6zbpfAsWLGDkyJE+D3Xx4kVef/0znblhw4b1mXFgt9vZ
   vXu32Z42bZoYB0Fit9u5ePGi2R44cGCbf9/o6GhL40t8MhC01h9j1EVIB/7uufxrjNiEZVrrfZbN
   SBCEPsFm+2zBXbx4sU+yx604nTZKS0+2KUrUEVVVVezdu9ds+1plz+FwoLUNl6uYmppKkpKmEBZW
   T1paYrfPes9p0qRJ3HzzzZ3e29jYyPr1681ngvpja7PBc8/BuHHw9tvw+uuwbRt0kTVw/vx5Nm3a
   ZLYTEhK6zfgoLi4255ucnMx9990X+JyDRGttHot4l86GvpdnvlLx/v2trq7m6NGjZjs1NdVnBdJA
   8FkoSWv9vtZ6LDABuBmjwuIYrfUHPTU5QRB6B7vdzt///nczt98XUlOHYLfnU1NTSVNTM2PHTmXE
   iK5T7gYOHBhQrn1DQwP5+UbMQWpqLamptZelYHaE2+3m+eefp6GhAYDIyMjL4hwaGxvNOgODBg3i
   scceC66YUnMz/OlPMHYsvP++kaGweTPceONlt2qtyc/PNxfV0aNHc/fdd/s13L59+6jxBI4opazR
   YgiQN998kwsXLgCGh2jevHnWFKb6nOJwOPjTn/5kBpYmJyf3qjZEn0gt9xYSgyAIvuNv2t6JEyfI
   zc1lwoQMoGNhJYCioiKioqL88kiAoU4YERERkCiTt9ZCd3EGmZmZJCYmBn+U0NRkHCX89reG+uHP
   fgbdFEXSWrN+/Xpuu+02n8WcioqKaGpq6rKCZG/RWsQpLc2IQ2lqarJGRfJzzI4dO5g1a5b5/1Jz
   c7Pf/w/0utSyUmq8UuplpdQ5pVSTUuqsUmq1UmpssJMQBKH3uXjxIlu3bjXb/qbtTZkyhTvvvJNJ
   k0YzadLoTnfzZWVl5g7XHw4fPkxhYaHfz506dYoPPvjMsdl+4dVam+l1YBRxC8o4aGyE3/0OxoyB
   3bsNr8E773RqHBQUFJCbmwsYf8jvvfdev5QelVL9ZgGurKxscyYeFRXVb+Z2pVBeXt7m/4/hw4e3
   +X8xUNVSK/A1i2Eh8AHQBLwPlAPDMBQVBwJLtNaZPTbLABEPgiB0jsPhoKyszK8zzK1btzJu3DjS
   07vWRauoqCAxsfv4gPZ0p17YGe0LDbnd7k5TM5ubm3nrrbe4//77g8u7b2iA//kfQ79gwQL4yU8g
   I6PbxwoKCnA6nT5rGNjtdt54441+oRPQ1NTEhg0buO++++ToIAi8K20eOHCA+Ph4Sz1Cva2keBiw
   A7dprRu8rscAW4FwrbV/BcZ7ATEQBKEtly5dAiApKSmg58vLy0lISOiyWE5LSwsvv/wyK1eu9Cve
   wG63s2rVKh5//HG/vBlaa1588UXuueceEhISOryndYdmSVngujpDu+C//svQL/jJTwyxo06w2Wxs
   3LiR5cuX+7yotq/zUFJSQkpKSp/szqurq4mJiSEsLAytNcXFxT4HlwqXc+rUKc6dO8eyZct6bIze
   NhCagfu01u908NldwFqtdWSwk7EaMRAEoS3Hjx9nwIABTJo0yaf7XS4XBw4cYM6cOd0ubt4leH0V
   89Fa43A4zN2+PyJA3l6DhoaGLrUU9u/fT0xMTHAiMjU1huLhn/5kSCE/9RR48s/b0xrN3zq/vLw8
   Ro0a5fN327hxI2PHjrVU9CZQ1q9fz5w5c0hNTe3rqVyR1NfXs2vXLjO40OVyERIS0qMemN6OQSjG
   qLbYEeFAUbATEQShZ7Db7WaUfEZGhs/GARjyvC0tLd1mN1y4cIH169ebbV8XwiNHjrTJlff1udzc
   XDZs+KzCe3vjwO12k5OTY7Znz54d+GJbVWXURxg7FvLyYM8eePnlTo0DgA8++IBTp06Z7fT09C6/
   m91ubxNvsXTp0j4zDvLy8kxpbYB77rlHjAM/0Fpz6tQp8/+5QYMGMXr06DYS3VfK8YyvfrzfAM8o
   pfZorYtbL3pKQD8D/KoH5iYIggW8+eabLF682Oc/8i0tLVRVVREREUFxcQWJiWnYbLYu0+dGjhxJ
   SkqKz/23nr9Onz7d5z+W3h6DUaNGmZHzHeFyuTh27BijRo0KvGZCZaURX/CXv8Ddd8OBAzC6Y9VG
   rTU1NTXmEceiRYuIjPTdqVpXV0dWVpYZD2JJESg/8C5VHR0d3eexDlcarWmIoaGhKKU4f/48o0aN
   IioqipCQEFPp8ErD1yOGV4CFQBJGuefWAk1zPO8zve/XWv8fi+cZEHLEIAhtXf++UFxczN69e4mJ
   Gd9l6ecdO3aQmprql0eiNV7g7rvv9jvt8dVXX+WWW24hObnj2nDl5eUopQKOrzC5dAmefRaefx6W
   L4cnn4RupI7Lysr46KOPOhUpal+zIiYmhszMTObOneuXIdETNDQ0sGrVKr71rW9JXEGArFu3junT
   pzN2bP9I6uvtGIR8jGJNnQ3Y2okC9P9n783Do7jOvO37dEuttRftQguSkMSOBQZjZDDYeIsd4zXg
   LTF2SJxlksnMO5NJnJkk9oyzzTeZSWbe+PXEidfEdrzFCxgM2GAwmwwYJECAdqF9bXW31FJv5/uj
   pKJbUkstIYSAuq/LF6qu6qrTbdB5znN+z++ZKq2ftQBB43LE5XKxadMmbr/99pADg8F75qG0fu7o
   6MBoNIb0DH9twVjquv2zBt3d3SPaBhcXF2MwGMa/WmtuVjwMnnsOHngAfvADpWdCEEpLS8nNzVVX
   3sH0E8F6VlRVVZGfnz8u46hzpbi4mLy8vIlvXX2ZcOLECZxOJ4v7S1n9s2JTgUnVIEgps6WUOf1/
   Dvdfjt/5KREcaGhcrhgMBhYsWDCmX/g7d+7k8OHDI17j9XrZvHmzqkeIj48PKTior6/nzTffVI9D
   DQ7q6+t57bXX1OPBwYHH46G4uDhAXzGu4KChAf7+72HOHKWhUnGxUr44Sg+EpqYm1aERgusnBnpW
   OBz1tLSUEBaWTn19GwUFBZMaHPgvlnp6egKstbXgYGQcDgenT59Wj1NTUwMaaE2l4GAiuTiUEhoa
   GiPS19dHVVWVepyXlxd0NVtaWklpaSU2m019vbCwkCuvvFI9Tk9PVPseWK3teDz1TJ+eQkZGxpjH
   Nm3aNG666aaQrvX5fOpElpaWNqLtsBCC+vr6MdlDB1BXB9/9LgyYJB0/rlQpBPmMZ86cCQiirr/+
   +lG3SQYsnAGSknLJzl46vrGeI0ePHmXHjh3q8bJly8a8xXO54d80y+12BxhCxcfHj9jp9FIhpABB
   CPGr/jbPw51LEkJsnNhhaWhojAWbzRag2g92zb591TQ0mKmvN/GnP72tOgpGRkYGiAVNJhOFhdmY
   zQ0IcULVH4Samfj4449Vt0CdThey/8Bbb71FdXU1oAQAg7MNdXV11NUpRVN6vZ5bb711TPoKAGpr
   4VvfUkyNIiLgxAn4r/+CUUSWMTExQX0WhqO3t5cXXniB1NQ4PJ56nM5eenqceDz1ITWZOhe8Xm9A
   VUR+fj7XDNMLQmN43G43zz77rCo+jIuLY9WqVRd4VJPPWHwQSoAHpJQVfq/fDLyIojuYcnUwmgZB
   41LG5/Ph9XpDniAH6wqqq48wa1YM8+fnB31PS0sLlZWVLFu2bExja2lpwWw2B3Uz9Md//3Y0L/+y
   sjJAmfDGTFUV/OIX8Oab8Nhj8H/+DyQnDxEQDggxvV4vL774Ig8++GDIQsIzZ85gNpvVewxoKII9
   43zR09PDxo0bWbt2rSY8DJF3332X5cuXqw6gY/HkmGpMtg/CUiAGOCyEWC+EMAgh/hPFfvkgMLq/
   qIaGxoSyb98+Dhw4EPL1breL1tYy9dhiyRx277Sjo4O+vj4AkpOTQwoOnE4n77//vtqaNjk5OaTg
   oLW1lZdeekk9Huzl7/F42Ldvn7rtkJ+fP/bgoKICNmyAJUsgORlOn4Zf/lINDgayKg0NZvburaS9
   vR1QMhRr1qwJ6XMMUFdXF+CrP5CVMZlMo/asOFc2b95MZ2cnoHyP69atu2gnuMmgurqapqYm9Xjp
   0qUB/2+07y50kWIJsAT4C/A8UAN8C/g7KeUaKWXb+RuihsblSTC9wABLly4dU9o4JSUOm62Szs42
   VVcwXKq7qKgooJlRKERGRpKfnx/SL1Wfz6cGEomJiUFLA0GZpN1uN263e0zjAZRAYP16uPpqRVdQ
   VgZPPQV+PSIGBIQWSwIWSwINDQ3s2XM26EpKShrxM7W2tgYYPRUWFjJ9FHHjROH1eunp6VGP58yZ
   c0Eb+0x1fD5fgKh0sFBz2rRpF7RV9lRkTO2ehRCPAM+guCceA+6QUlafl5FNANoWg8bFSrDSuKKi
   Iq644gqSk5NDuo/H46G3t1d1GgyW6u7r6xvTShkUN0OPx8OsWbPG9L7333+fvLy8oP4J1dXVeDye
   8deUl5bCz34GH36oiBD/9m8hiAbi+PFyysr6yMxUeil0dLSQnm5n7tzgjZQGl2zW1NRcECOc/fv3
   4/F4WLFixaQ/+2LkxIkTVFdXq5bHlzKTusUghDAJIV4FnkPJICwHwoEjQogHz3UQGhoagQxe2Q6U
   xuXl5Y0pRX306FEOHTqkHg+X6na5XPzhD38Y8yo9Kioq5BXrgNgL4MYbbxxxQtXr9eMrGzt+XPEv
   WLVKaZ5UUQE/+UnQ4AAgPj6G2to9albF52siIyO40ZKUkueff15N5UdFRU1acNDd3c2+ffvU46uv
   vloLDkbAZrPx5z//Wd2emjNnzmURHEwkoWoQjgI3AfdIKb8lpdwHLAbeAP4khHhpxHf7IYR4TgjR
   LIQo8XstXgixTQhxWgixVQhh8Tv3uBCiTAhxsl8UOfD6YiFESf+534b6fA2NiwkpJS0tZ3UD2dnZ
   owrm/LNmV155JStXrhzxOoPBwDe+8Y1RxY5SSg4cOKCWFaalpYWUTrfZbPzhD39QnzdYhOjxeNi+
   fbsaRGRmZo7aTjqA4mJYu1bprLhokRIYPP44BAmkiouLVZ3AtGnTeOSR+0lPt5GW1jXELRIURftA
   aloIwb333jsxXSFDYGArBggwZBoYi8ZZpJTs2bNH/XtkNBq54YYb1PPa9zV2Qg0QaoACKaXaHUVK
   2SOl/DqwFvjiGJ75PPCFQa/9ENgmpZwJfNR/jBBiLnAfMLf/PU+Ls/+X/x+wQUqZD+QLIQbfU0Pj
   omXAh6Cjo4WqqoP09VWHXBr31ltvqSVuQohhfzGWlZWxefNm9TiU0kUhBG63G5fLNeq1Ukr1F7XJ
   ZGL9+vVBf0Hr9Xri4uICJsOQOHIE7rkHbrkFli1TAoN/+icwGkd8m8/nC/BOGE1A+Pnnn3P06FH1
   2Gw2T9pk8+KLL9LS0gJAeHg4hYWF2kTnh8PhUAW1A9/LwN9PIQSpqana93UOhFrmOOJmvhAiU0p5
   Jtj5Ya7PBt6XUi7oPz4JrJJSNgshUoGdUsrZQojHAZ+U8lf9121BaQ5VA3wspZzT//r9wHVSym8O
   eo6mQdC4KHG73TidznGVxnV1dWEymUb8xehyuejr68M4ymTa3d1Nc3MzM4I0KQrG1q1bSUpKYtGi
   RcOer6iooKenhwULFozpvgAcPAj/9m/Kn9//vlKyGB0dVF/R0NBAaWlpwGpyJFwuF6dOnVLHNpnl
   bq2trXi9XrXfxFhsqS9HNm7cyNy5c8f89/NSZ7KtlmX/QwuEEN8VQvxUCDGt/7V8oOscx5EipWzu
   /7kZSOn/OY3AVtJ1QPowr9f3v66hcdHT2NjIa6+9FnJpXF9fH1u2bFFX7MFWuJ988ola1mUwGEYN
   DkAJEPwNd0bCf1W+atUqFi5cGPRao9E4JtMhAA4cgC9+Ee66C266SckY/N3fqcFBYLlilVr5ER8f
   z8yZM0N+jE6no66uTs1oTOYKtLW1lba2s0VhWnAQyJEjR9i5c6d6fPvtt2vBwXkkJANuIUQE8Gfg
   nv6XJPA+0IjSCvo0/dsC54qUUgohJmzZ/8QTT6g/X3fddVx33XUTdWsNjfNCamoq69atC/l6g8HA
   tGnTRp3IsrKyQspCNDU1ERcXR0REBMnJySFVTDidTp577jm++c1votfrh1REeL1eNm3axBe+8AUM
   BkPIVRgA7N0L//qviuPh44/D228rDoh++Is6pZTs2PEqPT1NzJmTT3p6otpGORh79uxh+vTpZGZm
   EhYWxq233hr6+M6Bzs5Otm3bphoazZ07d1Kee7FgtVopLy9nyZIlgOKDcan2PTgXdu7cGRA4TRSh
   duj4GXAD8GVgG8oqf4DNwN9wbgFCsxAiVUrZ1J+ZaOl/vR7w/5edgZI5qO//2f/1+uFu7B8gaGhM
   VU6ePElfXx8FBQUIIUYtOezp6aG9vZ3MzEyEEBQUFAy5xuVycfToUZYsWYIQIqC5zEgcPHgQszmJ
   xMSkEbc2BnQGYWFhREVFsWHDhqC/vPV6PTNnzgywcx6V3bvhySehvBx+9CN47z0IUqfu9XpxuXqA
   BLq7u/F6C7Bao2loMFNbWz2s+NC/U+T06dMnTXjY1tZGQkICQggsFgsrVqzQ9sn7kVLS3NysbrGE
   hYUF/J0ZqZvn5czgxe+TTz45IfcN9V/rA8CPpZSvAJ2DzlUD2ec4jveA9f0/rwfe8Xv9/n7nxhwg
   HyiSUjYBNiHE1f2ixa/4vUdD46IjISGBaaP0AvCno6ND7XUQDL1ej8PhGFX8J6VU09o2m43w8Om4
   XHk0NJjZt696WJMmgE8//ZSioiL1eHB1RVlZWcD52bNnjy6GlBJ27IDrr4dHHlHKFsvKFJ3BCCY2
   nZ1NNDbux2ptp7KyChDk5uYHlIj6U1tby+uvv64eZ2ZmDrvlMppZ1XjYtm2bWiYphCAtbcq51E8q
   A63GYWjH0NjY2IAmYhqTS6gixV7gNinlx/1Nm1zAEinlYSHELcA7UsqQNsv6/RRWAYkomYifAO8C
   rwPTUQKOdVJKa//1PwK+CniA70kpP+x/fTHwAhAFfCCl/NthnqWJFDWmLHV1dSQnJ4fs3tbW1kZT
   Uyd6vT7oyt7pdOJwOEhKCl7LPxiHw8Fbb73FV77yFU6dqg7o12C1tpOW1sWcOco+r8fjUSd5l8s1
   ZIXnj9VqxeVyhbadICV89JGyldDYCP/yL/DggxCk9NLj8VBeXq56EEgpsdvt1Ne3UVFxBp8vi/T0
   7IDPkJhoJDExESEEPp8Pl8s1YsloMLOqsVolnzx5EuCCmCldDLz66qusWrXqsg+UJpLJ7sVQDQTz
   dL0KOBXqA6WUD0gp06SUBillppTyeSllh5TyRinlTCnlzQPBQf/1P5dS5kkpZw8EB/2vH5JSLug/
   NyQ40NCY6pSWlgYI0kbCZrPx5pubOXasfcSVfU1NDadOjf7PccBhEZRV2sMPPzxq+t/tdvPMM8+o
   ZWQGgyHgPT6fj9dff121/7VYLKMHB1LCli2wfDl85zvwjW8oTojr1wcNDgY4deqUKswUQqiizpUr
   FxERYQtoVZ2WlsCWLVvU70yn043qJxHMrCoU/EtBTSYTZrM5pPddDpSWlga0Jr/99tvHlD3TmDxC
   zSA8Dvwz8A3gbaAbpTeDBXgTeEJK+d/ncZzjQssgaEw1/FfgY6G0tJK6ulji45UJ139lP1oHxOHY
   vXs34eHhQxoxDV41u911LF2aSXx8PDC6JXNtbS0ZGRmjaw2khE2blIxBdzf8+MeK2dEIArSjR4+S
   kJBARkZG0Gv8P0dx8Ul8Pi8LF84bV4Okwd0vB2dTgtHW1sa7777Lhg0bxvzMSxGn04ndbleDxZqa
   GlVYq3F+mOwMwv8HbARe5qwG4VNgO4pI8X/OdSAaGpc6Ho+H3//+9wENYkbC5/Px8ssvqy5+Ot3w
   k+d7771Hff2wGt0A/K2Uly9fPmyXRpPJRGFhNmlpXaSldRETY+PIkSPq+cHBwalTpwLU09OnTx85
   OJAS3n1X6az4ox8pPgYlJXD//SMGB6CURobaL8JkMpGWlkhqasK4uycOmFX5ZyKGM6uSUnLw4EH1
   +01MTGT9+vVDrruc8F+YNTY2cvz4cfU4KytLCw4uEsbarOlaFEfDZKAd2Cyl/OQ8je2c0TIIGlON
   3t7eUVPb/jQ3N5OcnIzdbh+ysr/mmhxMJlOAGj8YPp+PZ555hvXr14+qBHe73artstvtRq/XB71/
   d3c3brd79AoAnw/eeUfJGAih9Ei4804YYdxtbW3s2LGDtWvXjnzvfnp6evjggw+49957J6wqIJj5
   EgQaKO3evZuFCxeG5C1xqWO323n55Zf51re+pVVnXCAmKoMwpgDhYkMLEDQuNB0dHXz++efccMMN
   I042A9hsNk6ePMnSpUuHPVdf30Zvr5MDB3bz2GOPodPpQp7EXC7XqILI0QIJKSUvv/wyd999d2iT
   oc8Hb72lOB9GRCiBwe23K0HCMLhcLsLDw1UhYVtb24g6hoEKDZ1Oh5SS6upqsrOzz/vEtG/fPoQQ
   w2ZhLjeklGzZsoXVq1erGR7NAfLCogUIIaAFCBqTzeDJOioqioqKCtLS0kJSxDudTo4dO8ZVV101
   6nNMJtOISvuKigpKSkq46667RryXlBKPx6NmDYYLJPwDjdbWVrUaICheL7z+Ojz1lNIb4Sc/gVtv
   DRoYDPDSSy9x0003hZyCfvvtt1mwYAH5+fkhXT9e3G43zc3Nqv7B4XAQERExaoOrS5XGxkaMRqPa
   RvzEiRPk5uaOuWW4xvlBCxBCQAsQNCaTgcna57Pg8bgwGHrUyXokwZvD4cDj8YyYpj916hRWq5Wr
   r7464PWR7uvxeHC5XERHR4847uLiYurq6oK2wj116hSVlZWhuQt6PPDaa0pgkJAAP/2pYos8QsbA
   ZrORmKjs7ftvbwyH2+2mo6ODlJQU9f2hlomeC52dnezatYs777zzvD9rquL1elUjrF27dpGVlUVW
   VtYFHpXGcEy2SFFDQ2MUBsri3O4u3G5ryGVxJ0+eDCj7Go6UlBRyc3NHvdfp09tpb28FFBe6YMGB
   v2Bx3rx53HLLLQHn/QPr7OxsVq1aNfKDPR548UWYOxd+/3v43e/g00/h5ptHzBpUV1dTUqJ2fh91
   Rd7U1MThw4fV47EGB6EaH0kp+etf/6qWbMbFxV3WwcHBgwf5+OOP1eOVK1dqwcFlgJZB0NCYALxe
   L6dOVdPYaAlYzRuNdRiNRhwOO2fO9GA2K02DXK5ali/PHVFhf+DAARYsWDBiBmDwFoPdfpJVq+aO
   2AhJSsmzzz7LunXrhs1aSCl57rnnuPfee0cXH7rd8NJL8POfQ1aWspUwQr8Tn8/HwYMHWbJkSUi2
   y1JKPv30UwoLC8dVHurPaMZHbrcbj8ej7p1XVVWRkZFxWW4jtLe3c+TIEbUD5mhiVY2phZZB0NCY
   Qrzzzjt4vT0BZXFdXaepq+umocGMzZYBSIzGOtLSumhoOIzT6RzxnmFhYQEdEodDSklX10mmTbOS
   ltbF6tULggYHA+Y9QggeffTRIZO/f/fCYMGD382UTMHMmcqWwgsvwMcfjxgcDNy7p6cnwEhotOsN
   BkPI14/EaMZHu3fvprS0VD3Oycm5bIIDr9erOj6CYp41ffp09Tg8PFwLDi5DtAyChsYE4HQ6iYyM
   VO1+ARwOOzZbxrD6gOEMh7xeL7W1teTk5IT8XJ/PR01NzajvKSsr49ixY9x9993Dnj916hQnT54c
   PY3e1wd//CP88pcwZ46iMbgmmMmqQnFxMTqdjvnz5498736OHDmCy+UatpLjXBis12hqqsLnO8ma
   NYq2wl+IeTng8XjQ6XTodDp8Ph/vvvsut99++2UTFF3KTFQGIWjOTghRhdLWeeAhwWZagdKlWWvK
   rXHZ4PV62bp1KzfccAMGg0FNS5tMJjVlrexzK9f39XVz8uQWUlMVkeFwam+n08mRI0dGLdMrKSkh
   MjKS/Px8dDqdGhwMrqCIjIxU9+hzc3OHBBH+ro7DnR80OPjDH+BXv4KFC+GNN2CQYNIf/8k2NTV1
   1Ba9/gFTTk7OeWnpm56eSE1NJVbrwBjbiIs7W6p5OQUHAK+88go33HAD6enp6HS6oMGjxuXLSDmj
   T4Bd/X9+ghJMpKP0ZSgCalDaLOv7z2toXDbo9fpRLYX9nfh6epwkJ08jMzOwpt/j8ahCuNjYWO6+
   ++5RJ6r4+Pgh6f+B/fWGBrPaq+Hll1+muVnpzK7T6QL28KWUvPDCC3R0dADKdoa/4G9AzHfy8HF6
   f/ELyM2F7dsVF8SNG0cMDhwOB88++6wqdExOTiYhIWHE6//4xz+q15vNZrV8biIxGo2Ul39KXFwj
   aWldrFw5ixUrVkz4c6Yqhw8fDnA0fOihh0hPT7+AI9KY6oTai+Ex4O+Am6WUdX6vZwIfAv8lpXz2
   vI1ynGhbDBoTidfrpampKeRfqm1tbTQ0NBAerkx2w5kj7d+/H6/Xy/Lly0d87p49e1i+fHnQlfVA
   +jw21khYmAGrtZ3k5HYWLJgZcJ1/WWCwvgo2m42iHSfI+uBDpr/5NNY5C4j5xU+JvfbaIdcNZCwS
   EmKJj49XgxC73T6ikVJ9fT3x8fFq5mW8PSpGo6GhAYPBoJZRXk4GPlarlY6ODmbMUJK7LS0tGAyG
   0YWnGhc9ky1S/CeUhkx1/i9KKc8ATwA/ONeBaGhMdbq6uigqKiLUoNPlcuHz+ZgzZwZz5sxQuXqq
   CwAAIABJREFUg4MBMSDA1VdfPWJwAMrq32AwqJ0Lg2GzNXL48Fvq8eAJt6ysjI0bN6rHw5ra2O30
   PvEUKx+5g8SKY1T8bhvHnvwLZxIDg6LBGYs33tgYIPAbzWXx1KlTtLe3Bx3rRNHS0oLVqjaHveSD
   A3/hq9PppKWlRT1OTk7WggONMRFqBsEJ3CelfG+Yc3cCf5FShm4wP0loGQSNicDfIGY0urq6iI6O
   Dir0Gkjtr1mzRl3VDkdHRwednZ2jeh8MWBPb7Xb27q1Cr09FpwtTS/jCw8OJjIxUrYullMN/lq4u
   +L//F377W7qWLqP07h8QvlAJXIbrYnjkSCmVlR6yshYA0NnZhtlcT2ysEhgMzpa0tLRw5swZFi9e
   POLnOVeampo4cODAZelZ4HA4ePHFF/n2t799SeopQrEq11CY7AxCKfB9IURA+C2EiAa+D5w414Fo
   aExF9u3bx5tv/nVUY50Bdu/ezZkzZ4KeF0Jw3333jRgcgNLUqaura9Tnvf3229TW1mIymbjmmhwy
   MrpJS+tS6/tfe+012tqUX6o6nW5ocGC1Kg2U8vLg5EnYtQvxyp/oyjCO2MWwt7eHzs5q9binx8GR
   I40BGgj/7ysiIuK8rd79V8mJiYkUFhael+dMNaSUvP7663R3dwOKhuVSbZA0nMYmlH+PGudGqBmE
   G4APgK7+P5uBVOA2wATcJqX86DyOc1xoGQSNc8Fms/Hpp+WEhaWh14cH7Z/g301xuFK5zs5OduzY
   MaoAsbKykszMzFHLzPy1A4P37qWUdHd3qyK/oNmPjg74zW/g6adhzRql9bJfP4PBqzWj0cgnn3zC
   smXLiIyMHGI6VFlZRErKXNLTs/o/cxunT7/Bww9/5bwIDv0/78svv8y99947apfKS4GysjISEhKI
   j48HoK6ujtTU1PO2RTNVGMlSXGMok5pB6J/8FwLbgJXAd4EVKALFgqkYHGhojJf9+/fT3NxMfX0b
   ERFZJCSkDmusA8oE/fvf/57e3l5g+FI5s9nM0qVLR13ZlZeXB+yXD0dLSwuvvPKKejx4YqiqqmLr
   1q3q8ZDgoK0N/vmflWCgvh4OHIDnnw8IDkAp1/TXTgghsFgsqn7CZDJRWJhNWloXaWldFBSkERUV
   hdt99ntYvHjZeZm0jx49qlpTCyF4+OGHL9ngwOv1qhkCULYRBv6uAWRkZFzywYHGhUMzStLQGER5
   eTlJSUk0NLSHtGpxOBxDVsllZWXodLoRNQQ+n4+WlhZSU1NHHI/b7Q7YHhicNejq6lIn8YG/70OC
   kZYW+PWv4dlnYe1aePxxyM4OuGRw1qCuro729nauHVS9MBw2m4033/yQsLAE0tMLgmZbxot/1uTM
   mTNERUWNuk1zKXDw4EEcDgfXjeJQeakzmk22RiAXxGpZKMwTQlzb/+elt9mlcVlSX1+vTq55eXmY
   zeYAHwP/vfiWlhZ27typvne4FHpERMSorW+tViu7d+8etSpi06ZNVFRUqMeDV4ybNm1S9+GFEIHB
   QVMT/OM/wuzZ4HDAkSPwv/87bHCwb181dXUx6h6vxWKhoKAg6LhcLhenT58GlIzCPffcxFVXZQdo
   ICaC2tpa3n77bfU4MzPzkg0OWltbeeuts5UoixcvvuyDAxiasdKCg8kh5AyCEOLrwFNAkt/LzcCP
   pZR/OA9jO2e0DIJGKPh8Pv7yl79wxx13DElVD6ecbm1tpaTkJNOmpauvSSk5deoUM2fOHNE8yel0
   IoQgMjKw6GfwcwwGg3rN4IyBz+fDarWq+9DDWgQ3NMC//7vSSOnLX4Z/+ifIyAg6rtLSSs6ciebo
   0TdYteqb2O22Ufd4e3t72bp1K2vWrJlQYZzP5+Ozzz7jqquuQqfTIaVEShn0e72Y1e1er5cDBw5Q
   WFiIEAKPx0NHRwfJycmjv1lDIwiTmkEQQjwE/C9QDHwVRZz4VaAE+L0Q4sFzHYiGRiiE2q43FAb2
   03U6HQ888MCIwUF4uA+Xy4XNZuPw4UaknDtETV1VVYXD4RjxmZ9++inl5eVDnuOv0N658wTPP/+8
   mlkYnDGoq6tj165d6nHA5FxXB9/9Lsyfr7RZPn4c/vu/gwYHdrtd3ePW68NZufIb6PXBRZJ79uyh
   tVVpJx0ZGckdd9wxYcHBwOfV6XQ4nc6A5lIjBQcXm7rdZrOpTbh0Op3aRRKU/9dacKAxVQi1iuEo
   UCKl/PIw514GFkgpF56H8Z0TWgbh0mIi9yE7Ojp45513ePTRR4ed4AY/q6FhL1dcMR2DwajqEqSU
   NDVVkZvLqCvtgWzAcKv90tJK6upiMJvj0evDsVrbSUnpYP78s8LBlpYWEhISVB3CkPvU1sIvfgF/
   +Qts2AD/8A8wirYBYMuWLej1kVgscQHtqP2/W/9nlZeXk5KSMqoR0lj56KOPSExMHHFLYzguRnX7
   G2+8wfLly0lLS7vQQ9G4RDnvzZoGMQvFTXE4/gy8c64D0dAYDf92vaCU8NfXt40rQIiPj2fdunUB
   3RcHav1Pnaris89K8XoTueKKbGJjTcA1GAyBvgR2eytVVbtJTl5IaWklDocdKRUXwYFUt9vt5o9/
   /CP3338/LS1d6nMGj7mqai8pKblkZFwBDK0+2Lt3L8uWLVMFjWpwUFWlBAZvvgmPPaZ4GYywAu3p
   6aGyspL58+djs9kQIg2dLr2/qdQpjMa6/vErwUFlZSXFxcXcddddgKLPmAh6e3tpb29XbauXLl1K
   dHT0hNx7qrF//350Op3anXLt2rUXeEQaGqERaoBgBzKDnEvvP39OCCEeB74M+FC2Lh4FYoC/AFko
   TaLWSSmtftd/FfACfyul3DrMbTU0VGpra7FarVxxhTIJ+3y+gCxBaekxenp6qKuLpK1tOgZDET09
   epYtO+v+N21aPNXV1f0dAfWkpeVTV9dNWJiZEyfsCGEgK8vA7t2fsmhROrNm5fDggw9SVHRGfU5t
   bbXqchgVFUV6eiKZmTkYDOmqGDIlJYPGxkamTZsGoE7QKhUV8POfwzvvwLe+BadPQ4jCvQFB4+CA
   C2ZhNHaRlpagBjCZmZmjVlmMh46ODk6cOKEGCOPNSKSnJ1JbW612aFSEpNkTNMrx0draSmNjo/r3
   bN68eQGNsDQ0LhZCrWLYDPxMCLHS/0UhxDXAz/rPjxshRDbwdeBKKeUClA6R9wM/BLZJKWcCH/Uf
   I4SYC9wHzAW+ADwthBhTRYbGxUewqoJQiYqKCpiI/CdIiyUBqzWC6moXZvMcZs68Art9JW1tMVRU
   lKnP+vzzz4mNtatq6unTYzCbZ9HX58NsnkN4eAbHj5+gtdVOSYmHffuqKS+vDXhOWFg6lZX1PP/8
   83i9XkwmE8uXzwhQaLtcLg4dOjRUc3H6NKxfr3RTTE+HsjJ46qkRg4OPPvpI7doYHR3N6tWrg17r
   8/kC3PnCw8MnZGXv9Xp54403cLvdAKSlpXHTTTed832ngrp9oFx1AL1eH7D9YzQaR61o0dCYioSa
   QfgBsAzYKYSoAxqBaSjtnssIvv0QKjbADUQLIbxANNAAPA6s6r/mRWAnSpBwJ/CqlNINVAshyoGl
   wP5zHIfGFGZgMji7JTD6ZNDe3k5sbCwREREkJSWRlJQU9FqPpwef7xRSLiU6OpacnFSczlKSk8Mo
   LFzUb2d8DQaDQZ0AlIn7rMDOau1Cp5tLXFw8RqOSWWhuPkZUVAY+nxefT2m4FBERwTe+8Q11K8Fk
   MtHV1UVKSgoGgwGTycTKlSvVDEd0zWnMf/ouscUH0H3ve1BeDiM03vHXDeTk5AypmlC+P2X1feZM
   M0LoiIx0kZmZzde+9rUJER729vYihCAiIgK9Xs/ixYtHrPAYLyaTadKDAv/vt6+vj82bN/Pwww8j
   hCA+Pl6tMNHQuJgJ1UmxEVgE/C3KJGzv//M7wCIpZdO5DEJK2QH8GqhFCQysUsptQIqUsrn/smYg
   pf/nNMC/s2QdylaHxiXOYIe/0Th06BD19fXDnhuckUhI0JGRkYfNdpLGxkrc7jpmz46ira1Kndgi
   IiICJs+Be9TV7aCz8yAORxVhYT58vmaSk5VVfXKyBY+nnpKSDzl9eo+ajRisMygtLVVX+qBkOMx1
   dhb9+3dY8o930pOziLLNH8GPfzxicFBeXh7QtXHGjBnDZgEGAq6YmGaMxlZ19T1RVQkfffQRlZWV
   AeMItenVVOe5555TXS+joqJYv379JdkDQePyZko4KQohcoH3gWtR+j28AbwF/I+UMs7vug4pZbwQ
   4n+A/VLKP/e//gfgAynl24PuK3/605+qx9ddd51mOnIZ4N8bYTROnDhBXV0zmZlZASLF5mYryckW
   Zs/Ooaenh5SUlKATgM1m4/Tpanw+Hw5HNxUVNlJS5hATY6S3t4prr1UqA86caUGn06kixd7eXhoa
   GpgxYxjFfXExtu//kMiig7R8+R9oXfttOtyuoAr97u5utUzT5XLh8XiCbg3Y7XZ27drFbbfdNqGT
   WmtrK3V1dSxatAgI4s9wkXL06FGSkpLUygObzYbRaLxkPp/Gxc3OnTsDzNuefPLJSa1iAEAIsQCl
   F0M80AHslFIeP9dBAEuAvVLK9v7nvA0UAk1CiFQpZZMQYhowsNFXT6BoMqP/tSE88cQTEzA8jYsB
   m83GmTMtfPTRFr70pXtCKiNLSEggKiqKnJwc9bUZMzJwuRwsXaqU3A2XqWhpaWHv3r2sXr2a+vo2
   YmJi1Yl/5sw6jhypwOHwUlVVzM03F+B0OocELU6nk7KyssAA4fPP4d/+Dfbtw/A3f8Oub/8EEZsP
   bldQAZ7P5+Oll15i/fr1REdHYzAYhojifD4fdrudhoZ2fD6fKn48V/ybQYWHhwc0mrqYJ8/u7m56
   e3tJSFAEnDExMQHf6cVkxqRx6TN48fvkk09OyH1D9UEIQ9EAPDDM6VeA9VJK77gHIUQBSrnkVUAv
   8AJQhFK90C6l/JUQ4oeARUr5w36R4isouoN0YDuQN9j0QPNBuHzw9y3o63Og1wcXrNXV1ZGcnBxU
   We50OqmsrGTevHlBn3XmTAsNDfXY7eGYzbMARUE/Z46FkpJWIiOzAXC5zrBgQSLHj3cQFpZOR0c1
   UVE+Vq6cFTi2gweVtsuHDimuh1//OkRHB3UJ7OzsxOfzqRPYaFmTP//5zxgM04iPL1DHeq6CPq/X
   y+9+9zu++c1vXhIqff+Mx/Hjx7HZbJdN62iNS4vJ9kH4KbAW+DHwJ862e36o/1wl8JPxDkJKeVQI
   8RJwEKXM8TDwe8AIvC6E2EB/mWP/9SeEEK8DJwAP8G0tErh8cTgcbNu2A7N5eX/JXgJWa3tQj4Tj
   x4+j1+sDVtG1tbXExMSoGYXhgoP9+/cTGxtLfb2PsLB06uqcNDd3UVgYTmysCasVPvjgXXp7jWRk
   pJCUZAEyOXLkGFFR87FYEmhvP4XXm3B2bAcOKIFBcTH84Afw+uvgJygMJsCrqalBr9erAcLg4MDj
   8WCz2VSx3Lx5V9LWlnTOHhI1NTVYLBbMZjN6vZ7HHnvskggOWlpa2LRpE48++ihA0OBQQ+NyIlRJ
   8ZeBn0kpfyalrJFS9kopq6WUP0Ppz/CVcx2IlPLfpZTzpJQLpJTrpZRuKWWHlPJGKeVMKeXNAx4I
   /df/XEqZJ6WcLaX88FyfrzFxTKQdciiE0hipr69P/fmWW24ZkmJvb2/Hbh/ZziMzMxOXC7Vk0WiM
   Q69Po6FB0ct2d9tpbtbhcl1Ja2ssJSWNWK1tdHScLYHLzb2G2NhEoj4/BLfcAuvWwZo1SlXCd74T
   EBz443Q62bNnj3q8cOFCFixYEHSsVVVVfPbZZ+rxYHHleGlubg74fzpcdcTFgNfrZdOmTarddlJS
   Evfff/8FHpWGxtQi1AAhDdgT5Nw+tAoCjX4myxvf5XKpiv/w8HBWrrwGr7dhWI+E1tZWnn76/1FS
   cjpgLM3NzerPixYtIntQh8O+vj62b9+uljCmp6cTFXVW+JeUZMHrbaGqagutrXU0N5cya9YSoqO7
   AS9Op57W1tNYLBFqtQS73mfhP36R6Y//I3zpS4qPwTe/CaMEOAOllQMT2mB8Ph9FRUXq+fz8fG65
   5Rb1/Hg9JKqrq9m8+azNydKlS8nMDOaZNrWpq6ujt7cXULwKMjMz8XqVnVEhBFFRURdyeBoaU45Q
   A4RGYEWQc4UopYkaGkPMh8LC0tU99ImkoqKCo0ePqsfBDHNsNhu7d5djMl3HwYN9bN+u7C17vV62
   bNlCT09P0GcYDAbi4uICJuX09ETc7jpaW+vwePqYOVNyww3Xk5XloqAgjaSkaaSk9GEy1ZOQ0Mni
   xdO5Y80aVrhrWfb4F1nwX98h+usP4zh8mNIVN1BaURc0gPr444/VMkG9Xs8111wTVGcghKCnp0dt
   cDSYUA2FpJQBgVNqairLli0L+h1NZaSUagAASsXKQGkiwBVXXBEgqtTQ0AgkVJHiUyhmSf+GokEY
   MEq6H3gC+JWU8sfnb5jjQxMpTj7ns3nOwP/LUFLlNTU1nDx5kri4VHbs6MNkUkoNu7pOsHp1jFqh
   MJiGhgacTie5ublB771v3z5qauooKFisCgdtNhsnT1Zx9GgDYWHhJCRMJ8LgZEVvFTG//jU0N8O/
   /As8+CC2np6gTaf8qwKam5sxm81B0/iHDx8mPDx8xK2GseLxePjTn/7EQw89dNFPnrt27UKv17N8
   +fILPRQNjUllokSKoQYI4ShVDMNt0r0KPNLvajil0AKEyedcOy4GU+0DbN68maysLObOnTvqe1wu
   F1arlaKiUkpLM5k2bQZOZzXNzTtYuHAGt9++iuGora2lp6eH2bNnB7zu7zPg8/kQQqiBSmNjI1u3
   HiAj41q6ux00N51gZW8VuX9+kTCbTQkM7rsP+ts2BwuijEYDn3zyCQ899FDQ78fj8ajtn9va2lTX
   xXNh//79TJ8+/aLvLtjY2EhFRQUrVijJTpfLRXh4+EVdbqmhMR4mtYqhf/J/UAjxcwJ9ED6ZIB8E
   jUuE8dghDzA4uBhoajTw/hUrVqiT9HDvaW+v5OTJBm644QpMJhPJyckkJzdQUlKP3W5Gylj0+hmk
   pJx1IZRScuLECWbPno1er2f69OlDxiWl5JVXXmHdunWYzeYhaf7W1i6EiMVsimP6sf2seuYnRLrt
   hD31r7B2LYzgHuhwtAdoHL70pS8FvdZqtfKXv/yFxx57DCEEiSE2Zxru87hcLlXYmZKSMuR7DYWR
   grnJwO12U11dTX6+0hbbaDSSkZGhnr8Uqis0NC4kQQMEIcRq4DMppSrtllIeA45NxsA0Ll7G640/
   uLtgZ6dk69aPWbPmViIiIobt+Of/np6eZny+2IDyvdbWelJTXUhZDUBqajizZuUE3KOhoYHMzMyA
   MUsp6enpISYmBiEEGzZsCAgMPv/8c9LT00lOTiY8LJyC6kbyfnEVwuOm4sG/J2zdtcyZN3xrZP8O
   hMePv09W1kzS0+erfQv8aW5uJj4+nvDwcCwWC48++ug5r4hPnjzJ6dOnufPOOwECTKJCZbRg7nzh
   drsJCwtTv4Njx46Rl5eHEILY2FhiY2PP6/M1NC4nRsogbEdp0FQEIITQAzuADVLKskkYm8ZljhCC
   yMhIRtomslo7ATMAaWlzlUoButTzy5cvx+fz0dysiNPS0xOJjIykvr6e9PR0hBDDdhU8ffo0p0+f
   Zs2aNcBQn4Ho6Gj0QsDbbzPziSfo7nFR8ZUf0XrNrXh8jRRmJg+5p1L+WUFfXy/z5mVit3eRlrZy
   xNV3UVERS5YsUcsyx7MqdrvdfP7551x11VUIIZg1a9aQLZSxMjiYG6+vwlh57rnnWLduHXFxcYSH
   h3P33Xef1+dpaFzOjMVqWaBUMoyvcbuGxiikpydSU1NFQ0MH0dHxeDz1rFx5TVCRns/n49SpEhIS
   vHR3O2hvt2Gznaa8vBOPx0NWViqW/qZG/t31mpubOXjwIDZbn/pck8mEw+FQMwYzZ85k5syZ6nvs
   djsHDx7k+uuvB5+PWcXFiiVyRAT6n/0MVq4kvKGdNOzDbqsMrLgdDh19fT6czo5hV9zNzc20tbWp
   Rj0DAcp4GHAGDAsLw+Fw4PV6CQsLOy8dFc8X+/btIzk5WRWNbtiwQdVgaGhonF+0f2kaUwaTycSC
   BYm8//4mVq68gYyM4VPWnZ2dNDV1AnDPPfdgs9nYvLkElyuWzs5YDIYwKivDaWg4m/bu6uoiKioK
   g8FAVFQUkZE5NDQomYeB9Pg777zDrbfeSlJS0pA0flRUFHEmE/LVVxFPPQWxsfDLX8Ktt4IQmACT
   2TxkrG63m61bt5KdPYuwsHSys88KE4dbcet0uhEzJv6MpAF47733mDNnDvn5+QghWL16dUj3DBX/
   bRIgaJ+IsdLe3k5PT4/qtZCdnR2wbaAFBxoak0fQKgYhhA9YJqUc2GIIA1zAEinl4ckb4vjRqhim
   LgNlgS0tVlJSLMycmY25f4Id3AXQfyL0eLrZvn0XV1xxf/9xPRERDlpaDAgRT1ubGdCTlOQgLi5K
   LbH84IMPmDVrFrm5uWoVgdkcT19fN729faSldTFrVnbA6rqoqIjU1FSmp6XhfP55dL/4JV6LBfnj
   fyHmrrtgBC2A/2c4duwYEEFzc/yQyoWZM7N49dVXWbt27ahukIO/P38NgNNZyezZceTl5annY2Nj
   zzlbMFIQMlEixd7eXjVLVFVVhdVqVTtCamhojJ3JqmLIEEIMuNyE+b1mHXyhlLJy8GsaGsNhs9nY
   vv04ZWV6dLpMjh49yoEDh3j44S9hMpmGBAf+E2Fl5Wmysm4M2PtubT2AyxVNRER8wHNcrrMmSLfd
   dtuQcXR01FJd/Rm5udcDQ3UGKfHxJG7ciPfXv8YZbaH6W/9B56Jr8XgbKLTbg06Iu3btwmg0qpPc
   /Pnz+ydTZcUtpQ+X6wzp6Xno9XpWr149Zm3BUEFnLcXFJ9QAYSK0AKMJEccrRvWntbWV999/n69+
   9avA+ASTGhoa54fRAoQ3h3ntnWFek0DwWi4NDT/q69vo6orFZMrEaLRgs0UTHl45bMq9vr6NpqYO
   HI4G3O4k2tu9uFwtREYasNn66OnpJifHjMtlpK9PR1dXKUIYSEkxc+zYZgoLHwm4n8PhYNq0eGpr
   a9Hr05gx4zo1PW6329m5cye333wz4k9/IuvnP4esLGp++hTl6Tcr7pCA1SqGjNXpdKpWvQsXLhxi
   2+tf/llcfJi0tBT1/ePxH3C7XZw8uYWrr34AIXTExU0nLW3oFse5cD6EiD6fj9dee421a9cSHh5O
   YmIijzzyyASNWENDYyIZKUD46qSNQuOyYCAlXVFxhs7ODsLCIgELQoSj1w+fXrfb7VRV9dDZaQRS
   6Ozsor6+iNbWUiIjFyCEHotFx9Kl8dhsfSQnR6PT6TEaw1m9+hFVpDjApk2bWL58+bBeDT6nk6Wf
   fw7f/S7k58MLL8C119JTWjmimXhHRwdvv/02GzZsQAgxZAJ1u900NjYyffp0TCYT+fnTx7WX3t3d
   TUREBGFhYWRnp3HyZB319dW0t3fi9bYwa9bEOSpOJCdOnCAzMxOj0YhOp2PFihVqtsbfcEpDQ2Nq
   EfS3lJTyhUkch8Yljn+62ufTY7WWodefpK9Pj9fbwMyZuoDmQcePHycvLw+Hw0F1dQNe75UkJfVh
   NkcQFTWd8PBY8vIySUqy4PH0Ybd30dXVgtvtZuXKlep9pJTY/bYD1q1bp05IJpOJvXv3EubpxrR7
   N7pf/pKUefPg1VehsFC9x3CCPL0+EpfLhcFgID4+fkR/AqfTyeHDh8nMzFSrCsbDli1bWLJkCVlZ
   WZhMJq65ZjabNx8jLCydlJS5HD/eMSFp/wHGK0T0eDy43W41i2K32+nt7VV9LIYzo9LQ0Jh6hGS1
   fLGiiRSnDqWllVRU+EhNzUUIQX19DU7ncWJiYkhJsTBrVk7AxPbJJ5+Qk5PDhx8ep7hY0tUVRWbm
   IaKjbwMamTbNQkHBYqT00dXVSVpaF/n509HpdAFagubmZrZt28aXv/zloYNyOrH9538S+/TT6BYt
   gp/8BNvs2cMK7wYL8vbu3UtBQcGQttED7NmzhwULFpzTZF1fX09DQwOxsUqKPy0tQRVyDnyn56vv
   xQDjESLu2rWLiIgIrr766gkbh4aGRuhMqtWyhsZYGTyxSCmpqNiF2ZxIdLSFmJhY8vNnB0xmXV1d
   6gS4atUqSksrSU6ezfTpjdTW+mhtvYnY2BquvNJAdHQfnZ1tHD36OjNnLiY9fY66MvdPxaekpAT0
   NnA4HGx+6y3ubW9H9x//gWnpUnjvPVi8eERRXl9fH15vD/PnzwfgC1/4woiff/DWRqj4N2vq6+uj
   oqKTpKQZ/eOpmRS3Qn9CyUicOXOG4uJivvjFLwJw7bXXatsGGhqXABePY4rGRcPARNvQYKa+3sS+
   fdWYTBHMmrUYl8tLfX0NlZWf4XDY1VbH3d3dvPHGGwGtlRsaztDc/Dnz5mWTnt5DVNTnLFjg4I47
   ruLGG+eTnm5j9epVrFo1J2AS2759OzU1NeqxOlk5HMQ8/TT3fP/76D79FD74AN55BxYvBkZuVa3T
   6QJaBw+mrKyMLVu2qMfz5s0b80TudDp5+umn1e+grw+SkhYHbZ2dnp6Ix1OP1dqO1drevwUwvv4M
   Y6Gvr4+ioiL1ODk5OaAltBYcaGhcGmgBgsaEMzDRRkYaKC19H71+Gna7i2uuycForKO5+QQpKbOx
   2TLYu7cKm81GTEzMkH4HCxfOx2JJoKSkDr0+j2nT8ggLi2br1q2YTCbmzJnBokXzVCMkxcq4kvz8
   +SQlJan32fvhhzR+73uQm4s4dAj9Rx/B22/DwoVBP4OUktLSzfT1KW6LCQkJFBQEtoju6TlbRpmZ
   mTmutsJlZWXqfaKionjsscdC9i4YqIxIS+siLa3rvGYXrFarauA04Mw4EMhERESQkJC3FxoLAAAg
   AElEQVRwXp6roaFx4dC2GDTOG5GRRgoK7sTtViYWxeMA9Ppk+vp8uFw1dHTUkJVlVP0PNm/ejNmc
   SFxcAunpieTkWHA4PBiNkJx8BS6XE4ulKeA5VquVN954g8zMpYSHK938amuruGZePKaXXmLZb36D
   XL0aPv4Y+i2Mh0Oxeq5URXkpKWlkZaWq5/23TVJT43jllVf42te+RkREBJGRkUEtoUeiubkZs9lM
   dHQ0QIBZUigiwYkUJY7Ee++9x+233058fLzq3aChoXFpowUIGhNKa2srDkc7Ho9jyMSmbD2UU1YW
   hcdTT2xsHwUFZ41xbDYbDkcEOl06TmcEpaXH8Xg6sVq7SExcTGysCavVTWysEafTiU6nIyIiAovF
   wvLlq2lstGCxJOBtb8DzX/9C9Gd74Pbb0e3eDSE0JyotLUWn6yItTZlwr7pqSYBIcfv2IxgM6URF
   WaitPcNDDz00JvdDUDopNjQ0qBPsihUrgl57Lq2zz5UdO3aQlJSkai4efvjhSXmuhobG1EELEDQm
   FJ1OR3d3N11drbS1nWD27EwKCxU1/2efHaWlpZOWFh9u9xyqqytwOPaSmDiT2bNzqK9vIz6+AIsl
   AYfDQVmZntjYGHp74cSJNiACj6ceozGGv/71M2bNymNxv35ACIG+q4Npr/2WpNd/R+W8Zez+5f+w
   +P41I06qnZ2dxMXFAYrBUVhYmCoS9Ke+vg2nE6KifFgsCVit0NjYMaoYUUpJS0sLKSkpgLIVkZ6e
   HvL3OVkZgoaGBrq6upgzZw4AV155pZrV0NDQuDzRNAga50x3d7e6V+90OvnoowaKimKorJzGjh01
   qhCxudmK0xlHTEwu0dExmM3TSEhIp7u7L0Cc2NfnoKRkC0IkkJycSUZGEseO7eb995+lrq4Guz2D
   mJilHD7cQlHRUexVVUT+6z9T+PBi3NWn+eM3/sBf73gce+py9u2rVp8/GKfTyVtvvRWwl+4fHHR1
   dbFjxw71ODV1HsnJ+WP6bnp7e/nwww/VZ8TExKh+ABcSr9dLa2urejy4PNRsNhMeHn4hhqahoTFF
   0AIEjXNm3759lJeXA7B3bzE1NVFYrfFYrfHU1ESydesupJRER4dRVeXAZgsnLMxJT08XFouJuXMX
   oNfrVVV+T48TiKanp4zy8qP88Y8f09m5FIvFwsaNTdTXV3LmjJ2+M/FEPfFbIgoWkRIejvPTXZR8
   71sYZk+nsPAK0tOzhyj/q6ur6erqAhRR4GBhpL9vRnR0NImJSlXAWCoGdu7cSXt7u/qMhx9+eEq0
   WPb/bMqWyXb1ODU1lVmzZl2IYWloaExRLvxvLT+EEBYhxJtCiFIhxAkhxNVCiHghxDYhxGkhxFYh
   hMXv+seFEGVCiJNCiJsv5NgvN/wnmxtuuIF5/eK/urpm2toi8Xoz8XozaW+PoLGxgb6+PozGWGbN
   shAXV4LRWILF4sRgsJGenkh5eTnt7e0UFmaTnm5n7lwzXm8b+/efprMzh/Z2Fw0NNxMevpzPt2xl
   +V+fZsN/rCbK66Hof3dR8/gTWBYuJDc3k6ysbGJjh0/Lt7S0BGQUBpfkvfbaazQ2NgJKhiEsLIbS
   UqUPWbCKASmlmkEByMjIGJdg8Xzi8/l45pln6O3tBSAuLo4HHnjgAo9KQ0NjKjPVNAi/BT6QUn6p
   v710DPDPwDYp5b8LIX4A/BD4oRBiLnAfMBdIB7YLIWZKKX3Bbq4xMUgpeeWVV7j11luJj48PmGTj
   4mLx+fpwubrw+SLxet2kpWXh8XiIiYmlsLCAzs4uOjqsGAywdKkFk8lER0cHer0+YM89JaUek6mP
   sLBUoqOnIxqL2FDxJ25u2s6pq9byyg/fJW/VMjyePkDJCgxW/tvtp2hocKqGTEuXLg34LF6vl56e
   HjXtf9ttt2EymYKaJg3nUnjo0CGsVis33ngjgNpR8UJTVFREbm4uCQkJ6HQ6HnzwwSkXuGhoaExd
   pkyAIIQwA9dKKdcDSCk9QJcQ4g5gVf9lLwI7UYKEO4FXpZRuoFoIUQ4sBfZP9tgvN4QQfOELX1DF
   ff7Mnz+TsrJSYBNW6yKys/XMn5/Dxo0bufLKK4mI6CYjIwOLxURV1W58PmUyzs7OBgJFg1FRM8jJ
   yaP7VBEPnnyGL7S/x578BXz8o/+hWZdOSko+Hk9fQPnfYOV/QkIeLS0tAWP0L1e029vo6uripptu
   AlCdHEfqZNjb28uxY8dYsmQJoAj6poI5kN1ux+v1qsJJo9EYoKnwt2nW0NDQGI0pEyAAOUCrEOJ5
   oAA4BPwdkCKlbO6/phlI6f85jcBgoA4lk6BxHujr6+PIkSMsXboUIcSwxjgDqXuLxU1rayzh4SeJ
   iwtj2rRFLF68Dp1OR2qqjZMnq6ioqCcyMpbOzmS2bz9OZmY0Ho+HvXt3kZ9fgNEYgThTwt3bPuGH
   JR+wffpy/vWmX3Dtl2Zxx5qb/Sb5riHlfx988AE33nijqh9ITk5Wz3V2drJp0x5SU5chhMDttnPN
   NaP3DJBSBhgF2e12pJQIIS6ovmBgDACnTp1Cr9ezaNEiALUiQUNDQ2M8TKUAIQy4EviOlPIzIcRv
   UDIFKlJKKYQYqfuS1pnpPKHX6+nr6wuYkPzZsWMHdXV2zOYr6Ow0U1PTxqxZnXR25vH22wcoLMwj
   MzMFt9uN0WgkN/dqwsLCqays4NixOqZNi8HhCKO2NhtDvZtrP/0Nhcd3cOjKK3nnl38lNieXGyKU
   vX8ILP8b0AAMeBLceOONxMfHD/s5Ghra6elxYzSaCAszBGQG/DMLRqOhX5SovO/EibeYPv16QAkQ
   rr/++on8esdFTU0N+/fv57777gNQMxoaGhoaE8FUChDqgDop5Wf9x28CjwNNQohUKWWTEGIaMJAv
   rgcy/d6f0f9aAE888YT683XXXcd111038SO/CAmlS5/P58Nut2M2mwkLCwtoozz4Ph0dLmJiZtPX
   5yM8PJ+EhNl0dtbR2enF4XATG+uhpqYUi8VLevoMurvt1NbaaG2txecLo6yslRkylscOvE5B1VZO
   rHyUt37+PWYstbDIEgd4gxoFFRcX09jYqDZQGsgcDHDo0CHMZjN5eXnodDpyc1cSFmYY8jn8NQfd
   3WXk5MTQ21tHc7OVhQsLVC+DC0VfXx87d+7k5ptvRghBeno6d9111wUdk4aGxoVn586d7Ny5c8Lv
   O6XaPQshdgFfk1KeFkI8AQw4tbRLKX8lhPghYJFSDogUX0HRHaQD24E8//7OWrvn4Rk8GXo89cP6
   +JeXl1NaWsqaNWuG3ENKyccff0xHhwGHw0JDQwtudyvx8eHU1+dw6lQv4eEGwsP70OnOcNNNc8jI
   yCAtrYuYmDA+/PA4VussHA4P+ooT3HbkGWZX7GP/lbdx9Lp7ic64noiIagoL4wKEgTabjbq6Vtrb
   2ygoUJo0eb1edDpdQGbD5/Opqf+Ghgaio6OxWCxBP3t9fVtA6+Ty8gOEhTURG5sz6vd0PqmuriYz
   MxO9Xo+UkiNHjlBQUDAlyiY1NDSmJpdqu+fvAn8WQhiACuBRQA+8LoTYAFQD6wCklCeEEK8DJwAP
   8G0tGgiNkQR4/uTl5ZGbmzvsPYQQdHbaqahIxGzOJCoqgerqdqKj9TgcpfT2xhAZmUZCwjGkPNvP
   wOfzsW3bNubOvYLmnSUUbPsj04p3U3T1gzyx7B46fIlkR8zE3VXKzJm6AK+Bgcldp0ultHQ3dnsY
   K1bkDxl3a2srmzZt4pFHHgEgLS1NPRfMvri8vJZTp4pYuvQ+hBAkJubhdPaG9D1NJANah4EAoKSk
   BLPZTFxcHEIIVV+goaGhcb6ZUhmEiUbLIAQysB1QUXEGny+b9PQsAKzWdtLSupgzZwaffvopRqNx
   SOdCUCb3pqYmdcLduPETiottpKWtJiwshqqqY8THl2A2R9LUFIvLpaOz04vRmEJqqp2sLGVyNtbU
   4HniCXwf7aD2nm+yMctMdOp8Zs60ICX09HhISbEwa1aOOhlXVlbS2NhOb+8Mvwn77Ljb29uJi4tD
   p9MhpcTpdI5qFWyz2YiNjUWn09HV1cWHH35GfPxChBB4PPXEx0tstoxhn3e+2LJlCykpKVogoKGh
   MW4u1QyCxnnCP7Xu8+koKfkcgJiY2IAywQULFgStlbfb7ezevZt169YhhCAlxYKUzTgc7ej1vYSH
   W5k/P4vy8lLmz59BeHgG3d0Oqqp2ERtrYLlRR+yGDfh27aL1/odofGcj3ugYbo6JITMzecSVufv/
   b+/O4+I6z0OP/96ZYYaB2Rm2GUAIjBDarMWyjC3Zsh05VpzGWbq4dbM1nzg3bZruTdLVTZuk7b03
   bZLWt2lu2qStr7M6sR3FbizHsmIbyZGFJCQhIXYxoAEGZoYBNOt7/5hhDALJaJel5/v58OFw5sxw
   5tUR55l3eZ5kknQ6fdbHd+7cyd13301ZWRlKqUXVEXjmmWfYtm0bZWVlOJ1O7r//1jk9CwAtLeeu
   pnixent7GRkZYePGjQDcfffdmM3mN3mWEEJcftKDcINob++eM8YeCPRiMPRRV1fF2NgQGzduXLAy
   odaadDqNyZSNJVOpFL29vdx0000MDAzw5JN7iUYLcbmK8fnMvO1tq7Db7fT39zM5mUIpRdXIAPZ/
   +AfYs4fTn/gE3y2pRRdV4PevPeu4fiAQ4IUXdnHLLbdRVVU6L3nR5GSI6ek+tm/ftOgu/56eHhKJ
   RD6l8NlWZMy2mMmc5yMejxMIBKiry/ZCjI+PMzmZzQ0hhBCXgvQgiItSXGzH56umqamOlpbhOcsE
   Z2ttbWVsbCyfJTAej3P06FFKS0s5cmSM2toNdHa+Qio1yq233ovD4UBrza5du3jfkiXYvvhFeP11
   Jj/xCYoff5yevlNUBhy43dm5BbPH9bXWRKNROjp6aW0dpKhoKYODTk6e7M0HETPzB0KhcczmbOKf
   mVTIC93A0+l0PlmQ2WyeM7lvMcmNLkU1xWQymS98lEwmOXLkSD5AcLvdCyacEkKIq016EG4Q85fx
   dbJ1a9O8m9+Zn5iLiopIJpNkMhmKi4vzxx0+fIJg0IPLVUImkyYcDuFyDbFx482wdy/6s59FHToE
   n/oUqQ99iG9861u8//3vp7s7MKcnY/a4/hNPPIHRWEYs5iIUcmO1plm9upJUKk5ZWYj29gO8973v
   zfdmvNlqjHA4zLe//W0eeeSRS5Lp8EJ6E7TWfPnLX+aRRx7BarVe9DkIIcSbuVQ9CBIg3EBmbnDx
   eJyWll088sgjc1Lxztxw+/u7KStrpLAwSXNzLUePHkVrTXNzM5AdZvjKV/6JFSt+Ga83O2Gxq2sf
   7vbn2NbyMrS3k/yDP2DqoYdwlpXllyYqpbDbzRw5MobJ5EfrDFNT3flApbX1CKFQOePjUUZHnWht
   oLR0Ao+nGJ8vQmGhwu12MzQ0BkAsNjFvEmEqdZB77tmS/8R+tp6RC2m7xSwNheya5GXLluUnc6ZS
   qXxQI4QQl5sMMYjzNru7/Oabl8/7VJ0NHpxYLDWMjcUJhU7R1XWUDRuWs2JFXX7M3mQy8fDDv8ar
   r55gfNyMu62Fe/7jc3jGhuDP/gw++EEOtbWR7O5mRWHhvBvrypUeJiYinDo1SCJxCocjm+q4sNCK
   UoqyMi/BYC+RSDvpdDEORw1+fy3RaJTvfrcFo7EMr7eEYHCI8vI3eiMAotEIk5OT+XoElyI4mGmb
   sy15DAaDZDIZKisrAaivr59T90CCAyHEW5H85bpB9Pf3c+jQId75zncCc8ffU6kUBw8eJJUy0t6e
   wGQqZ2zs+xw54qGmpp6JiQTt7S/icul85r5CiwW16/vc/PxO7JFxMp/+FIMPPEDV0qUAbNiwAcjO
   D5i5sWqt6erqIhI5TXV1GVpryssriUajOBwOystdtLcfxmRagsEwyeTkKJWVKrfcMMqzzx4mHG7E
   ZnPS3x/A662iu3s34XAnfv86UqkAb3/7vVckkVG2gmTW+Pj4nAChurr6bE8TQoi3DAkQbhCVlZVn
   XfqnlCIYDBKNZggG05hMCYLBNUAxJSU+HA4n6XQBNTUFhEZHKXzlFdJ/+ldsC4fpefgPaL2rmY3N
   S9n97LO8p6LinGPtkcgA4+OagYGpfK/CzCTETCbD1FQ/09OTaF2E2dzAyIiLoSEXBw60kUh4KC52
   YjYnSST8TE31sn79MhwOJxbL/KJNl3IFwuwy0pFIgIGBFpqbfxWA5cuXX/DrCiHEtUoChOtYb28v
   FouFyspKCgoKMJvN+Rn/xcUmgsFxbDY7druZZNJFf/9hXK5bGRw8islURGNjALPZBTgxKCMFLzyP
   +vpXMSWS9P7KHxK+/1eJTI5z6tQYiVeyvRNnBgep1CTDwx1Adp3/kiU3YTZr4vFsr0J3dwsORw2B
   wChNTXVs2HA7g4NOxsYiGAxOwEgiEcNoLEOpFMlkP5HIEUym2zGZhlm79o4Fb/xnzhno7++94DTJ
   8XicHTt2sH37doaGxqisLObeex+6oimXhRDiSpMA4To2PT2dL1E8N4fABAcPvoTbbaKp6T66u3+O
   zebD5Wqkr+8UJSV1hMM9RCKKIusY3j07uPfVf8Np0hgffZT2VRs4FfQwNtxJX98AWq9BKTMtLfNv
   wiUlJaxbV4DWESCbgOjYsR76+wcZH58mlcqQTqfO+T6Ghlrw+VaRTvdRVdVEKHQzqVQH27evPutN
   erHppM/m4MGDLF++HIvFgsViobm5GYfDkZ/bIIQQ1zsJEK4zExMT2Gw2lFI0NTXl9wcCoxiNPoxG
   M93dAaLRRkpKJkgkxkkm3bz22jHq6m6huPjnjIxE2HiLifUDfZT/659RgGLst34b48c+QnBkBF95
   Gf2BPiYmFFqvwWpNU1e3lFQqTk/PIAcP/pCHH34Yo9GIz+ebUwth//797N9/lHC4ltHRYjKZAhoa
   BvH7VwLZrvy+vm5MpmIikXaUMuPzVVJYOMn27auYmEjQ0FCM33/7Jf0EP7OUc2ZSYzQaZXp6Ov9z
   bW3tJftdQgjxViABwnXmqaeeYtu2bQuWJj5xYhdTUw4ikTIiESNDQyPU1Ci0BqXMGI1OigpvYsvw
   07z9C9+isNDCoQ//HmNbfh0MBjr39jM2doQVK9bidmsikSBKmamtraGoyEo0GsdsNnPffffNWT6Z
   SCRmpQ+2UF29maYmD8PDowwPx0mlJvJzBRwOB0bjCEoNcO+91WgNdrv3vOcQzJ4zAG+eJvmFF17A
   7/ezevVqALZs2bLo3yWEENcjyYNwnQmHw/k8AbNvqtFolMcf38X4eDEWy1K6ug5RVORBqeOkUlEa
   6lazvqeLxu9+kYRKkfjUZ4jdfS/dPQqLxYzbXUUg0EsweIy6uux8gkikA9CMjAziclXjclnnDTEk
   Egm++tWv8vGPfxyTyTQn5XMsFmPPnqN4PEmKizOUlDhpbq7NF1C6WOeapNjZ2Ul/fz/33HMPsLi0
   y0II8VYgiZIW4XoJEM51ozt9+jTf+973eOihh5iampozMS+ZHGB6upv7778fq9XKa68d4OWX26io
   WEcqNc0PfvAqlWVu7goeZvOu/8RSWsrAh/+E4Pommm9fSiAwytGjUQyGDDU16zl48HWUMrNq1Qpi
   sVEyGRN2+wBWqxWTyZQ/t+HhYaxWK3a7HZibKGhmLoTR6KOjo4vOzpOsWrWEgoIoPt/ay1YtcWpq
   iuPHj+erJE5OTpJMJmVOgRDiuiOJkm4Q0WiUnTsPEw5nx8LN5i6WLSvDbrfnb8j33HMPJpNpwYl5
   RuM4kUgEq9XK8uV1hEKakZFTdB4bZmvgJA//7MvEbWZefNfvE91Uxda7N+KIDGOz2fD7ob8/hslU
   TTgcIp0eprx8BZOTY3R1vcLSpXdht9vn3dA7OjqorKzMBwizEwU5HA5uu20J3/jGNwmHi3C7NxKP
   VzE1FcDpnLikbRcOh/MBgNFoJBqN5h+bnTZaCCHEfBffjysuq+PHe+joyBCP1xIOl/DTn46wd2+Y
   trYQLS29RKPROZMAU6kEw8Mn8j/b7XZ2794NZG/Ot2+spuGVJ/ntxz7J/f0v8MNt/8iXfuk5ehre
   Qc0SH01NdUxOTjI9PZ0vjuTzRTh9eh/33tuIxRIhnTawdOlduXF9L/39/Tz11A7a27uJRqNs3ryZ
   +vr6Oe8jFAoRCoUAcDqdbNlyD+vWbcPtLgKMTE/bCAbb8fu9l6TdtNZ85zvfYWpqCshmVLzrrrsu
   yWsLIcSNQHoQrnHBYBijsRq7vYRoNILJtJpEYhxIYjLVzVm65/d7OXHiKMHgKeJxAyMjHdx8s49N
   mzbR39lJzUsv4fz851npKeHZX/0AL6ktmM0rSUwOkNEHWLv2QZRS+cqN8EZ65nA4iMfjobnZO2u4
   oxaAtrZhJicdDA46z5pvoK+vj8LCQkpKsr0bVmsRxcU2Vq+2MDISprBwkpUrfRe1MuHHP/4xjY2N
   1NfXo5TikUceueDXEkKIG50ECNe4oiITQ0PHicXiTE+HyGTA662gvn4Z4XD2E3kqlSKVSuFwONi6
   dQVPP91Db+/L1NVtJRYqwPC1/8nKZ54gvXIlxm98g4DXj3fQyTvj0xw/3kUsNkhDQzlVVVX539vT
   00MwGOS2224DyBdqgmzQ8MQTT1BWdh/DwxEKC5dSUTE330AikWDPnj351M7r16+f875mVhmAH7fb
   it0+xvLlS8+rbfr6+kilUvneiubm5vywhhBCiIsjAcI1LBqNEg4rHA4jicR+MpkYHo8Ftzs7J2Bo
   6AChUJxdu3ZRXFzIxo23YbebcTjKmQzFafzp91j6rS8RrVnGK7/1SZY89Es0NDTgi0TYs+d5/P4t
   rFmzilTKzcqVHo4e7UIphd/vxePxzCl0lMlkiMfj+UyJ27Ztw+PxMDwcyR8zNRXOJ2ZyuVysWbPm
   rO9tZvhidm/Em/UepFIpIpFIvhcC5taUcLvd59G6QgghzkVWMVyjotEou3e3Mjxso6KinFjsNBMT
   Yerq4lRUVDI4OMiuXf2YzSvp6wvg9R7lpptuJzRwkuYjr7D2J//FiL+RyCe/wHjD6nmrA3bt2kVB
   QRHj49MUFZkYG4OurtfxejeTyQyzffuqOT0K+/btIxqN5pcFzj7PmVUJbW1PUl9/c75886Uwe/lh
   IBBg3759PPjgg5fktYUQ4nokyxwX4a0aIMzcdNva9pNOr6a4uIDVq2tJpZL5G/2XvvTPnD59G+Am
   HHaSnujkzvavsvXnOxi7aQUvbn6ITmcCi6WeFStcrFlTxtjYGHZ7dhKg3W6mrW0Eg6GcQGCI/v5J
   rFYzBQU1DA72YbcfY9u2Rm6//XYg24OglJrzif3YsWNYLBZKSkoIBEbRWlNVVXrW2gjnWzgpkUjw
   9a9/nUceeWRO4iUhhBBnJ8scr2MzyxWXLy+mpyfJ9LSZrq4T+P0F+YmBXm8lvb0nsabNbNn7H2ze
   +0+crKnie7/xz5xeXkln5zChASNe7zQDA3EslhjDw9PU1GTH67u7XyOTiVNamsFudzM5mWZ6ugCt
   J0mlypmYmOall36OUoW4XK78TX12TgOr1YrZbM5PZDyb8ymc9PLLL7Nu3TqKi4sxm8382q/9mgQH
   QghxFUiAcA2Y+XSdTqeJxUI4HKUAVFbWY7fH6O7uwe0eB9LY7asZGBggHY6y4plv8ba2vRwrX8P/
   3PbneO+uJBQaY/RQjFRqPSUlbYTDThyOWlpbuzGZynC5hqipWYXRWIbRWMCSJWuIxaKk022cPh3C
   bF6J2WzF5SpnZGSCY8cKqKnJrk6oq7PS0tLCBz7wAQCWLFmyqPd3rsJJkUgEg8GQn1zocDjIZDL5
   5zqdzkvY0kIIIRZLAoSrbPan60wmTU/Pa9x3XwWpVCBfR8Dn09x66wb27t3LkVdfZfTRx3hfy3N0
   1TbwpXe9n5I7N7Cp2geYcToNZDJGDIYCjEYfqVQNyeQUbrczV9K5lerqlXi9JQSDxxgf96GUYv36
   Ytrbe0kkxqmoqGRsrBOXq550egSDoZpAIEkiMcYDDzxw0e959rBPW1sbJSUl+cJS55rYKIQQ4sq5
   phIlKaWMSqlWpdQzuZ89SqnnlVIdSqmfKKVcs479jFLqhFLqmFLqvqt31hcuGo3y0kv76e+fxGQq
   wOMpo77+fqanMzQ311JSMkxRUQ/NzbW4lWLtU09x0/3bKRkO8/8+9i12PPwN0subcLlKKCgoJ5Uq
   JRKBggIndnuUdLoCpQ6TTnfR2FhPJhPC4VhNJDKGxRJl48YKenqewueL8K533crHP/4+li0DtztC
   ZWUBFssU6XSYtrYThEJuRkdL2LcvMCcj4WL4/d5cwBOiq+vndHT8KJ8QafPmzXOqTgohhLg2XGs9
   CL8DHAVmFrN/Gnhea/33SqlP5X7+tFJqBfArwArAD+xUSi3TWmcWetFr0UzPQSAAo6NHiceLqavz
   EgyGmJ4ew+/3YrMVMHysE8ezT8Njj+F9xzv42kcf4cDkHXjM6yEKIyNLsdkGKStrxustQOs009MR
   lJpE60Gi0TY8ngZOnQpSVTVNdbWR0tIwVVXZokh1dXX5uQAOh4O6uuPE40OUlpZz8uQU4XAN0ah7
   TknnmcmGi5l0GIvFeP3112lu3kAgMEpZmQOf7xcuaalmIYQQl941EyAopaqAdwCfA34/t/tdwEx+
   3G8Cu8gGCQ8CT2itk0CvUqoTuBXYcyXP+WLMjMs3NdVx6FAl4+Mpdu58jdJSNx5PNfuea6V579NU
   /+u/knzf+yh47TWmvF4m//UbxFqPMj5uw+Px43R24nR6ALBaHdjtbqzWA0QiaSoqlpJK3cz4eAEe
   Twqns5jOzqOsXt2Yv0EXFRURDAaxWq0EAqM4naXU1vrwer00NmaXWiplpq5uKTabjXA4zsTEBC0t
   sQUnHWqt6e3tpba2FqUUVqsVl8uVq9kgQYEQQrxVXDMBAvAPwB8Bs+8i5VrrYPxHSMcAABUGSURB
   VG47CJTntn3MDQYGyPYkvCW0trYSCAxjtd6CzWZjzRo/L7+8F4/HQBlHqf7nf2L5y//N4B2bGf/3
   f2fAZOIer5edOw8zOLgEKGNqagibbT9r1tTR2Ohlaio7ZyGTMVJQUMzateuJxzNkMj60TlFebsfl
   stLUVI7H48mfy8jICK+++io2W0Puhu/k9dcDNDdnVyfceec6Wlp6SaXihMNxUqkASjFn0uH4uKa9
   vTsfIBw8+BplZWUUFxdjNBrzFRSFEEK8dVwTAYJS6p3AsNa6VSm1daFjtNZaKXWupAYLPvboo4/m
   t7du3crWrQu+/GWz0Pr/2tpaysvLaWt7YyLiEkuQW3b9N0t++iQHl2/gs+/+HomKKJvSXqZjvfT0
   DBIOW6isrMJqPUksdisFBScxmcZpa9vHPffcQ2fnYQCamsrJZOzE4xHS6RhaJ3hj1AaeeuopNm/e
   nM+CWF19ExMT2Rt+LBajv1+xe3crd965bsGMh4HAKLOnIRw9+ixKGVmx4u0AeDwrSafTl7llhRBC
   QDbx3a5duy75614TiZKUUp8H3g+kgEKyvQhPAhuBrVrrU0qpSuBFrfVypdSnAbTWf5t7/nPAX2qt
   957xulc1UdLsFQqhUDfFxQa2bGnId+9Ho1EG9x0i+blHWfH6Pto33MUPGj5C3DtIKrWSsrIaXK5R
   Vq82EQyGOXhwisLCBszmJKOjcczmEM3NhRQWFnLkyDDl5U1YLGaOH/8xfn8TBQV+jh4dAIw0NpZS
   VDRFc3MtXV1dDA1pLJbsMsXu7tcoL1+B01lCW9sQ09NGSkrG8fsLFsxXcODAAY4cCVBRka3T0NnZ
   QkXFSqqqsrUUwuHQvMyNQgghrozrKlGS1vpPgD8BUErdBfyh1vr9Sqm/Bz4I/F3u+w9zT3ka+H9K
   qS+SHVpoAF674if+Jmav/x8ZOcLkpJPdu1upr6+mijiOxx7D8fjjdNx+O11PPcWU1YXh+VGMiQps
   tknC4VHGxroIBoeor99KMnma7u79LFlShNk8yZIlJYTDBsLhAsbHlxGPR1m1agnl5WuoqjLjcGSo
   rLRz5MghCgsnaW6+FYfDQWGhE4vFmR8iKC9vIhhsJxTyMz1djNUao76+gVQqSSAwSkFBAcPDw/m8
   Bz6fD5vNRjKZ7YFwOKqJRmV+gRBCXE+uiQBhATMf+/8W+I5S6iNAL/DLAFrro0qp75Bd8ZACfvNa
   zKmcSMTz25WVa2lpOUSdKU7R//lLzLufIv4bH8bS3o6KRtm/fz/337+WgYEp2toUnZ0JIEphoQmT
   yU8qFePuu5s5fLgVpYa45Zbs0sCJiSoGBnZSUFCP0biUkZEQNlsp4+PdbNq0FoANG1ZjMJx9RWtx
   sZ21aysJBsfQOkF9fQNWaxETE9kAYHJyko6OjnyAUFZWRllZWf75Mz0lM8MlqVQgn/FRCCHEW9M1
   McRwuVzNIYZkMsm//Mu/UF9/BxbLEkL7XmHZk//GymO7GfqFD/CjBh9b3vtuliypoLd3iFdeeZF3
   v/tBiouLefLJ5+jpGcFqvRmtIZl04HYfZ9Om9zAxEaGiYhybrYC2tl5GRjx4PHb6+uLE4wWUlIzj
   9SZxu5Ns2bJlwXOLRqPs3HmESMQGgNMZ421vWwlAS0svSpWxf//j3HzzPdxxR/2iliReSK0FIYQQ
   l54Ua1qEqz0HIZ1OM9XWxum/+CtsP91J2+YHSf7m31Dor2XPnifw+Wxo7Z13oz5y5AQvvHAMu30N
   6bSXw4f3smFDCdXV1QSD7ZSWavr6Bli69L7cHAMTWrfg8TSwfn01jY1L37Q2ws6dhwmHs+WcXa44
   RuMId955J0ajMZ/2uaamXG70QgjxFnNdzUG4XgwPD/Pzn/88m464owPj5z5H8Y9+RHD7L/P0o19h
   6DSYjg6wArBaCwiFRgmFynE4qolEfszQUClVVT2YTCYMBhdGow2z2Ux5uWJsbC+joz00Nm5kfDzF
   5KTCYjHR3LyGrq4TuN0ruffe2xZ1Qw8ERjEYXNTXl2GzeQmHQ5hM2ZwFhYWFEhQIIYS4tlItv9V5
   PB7WW63w678Od9xBpq6OEz9+gb4P/w1L1m7DZktRWmrDYOjjF3/xPsrLl2IweLHbXbjdG0mlehge
   DmOz2Vm1aiXl5ZqiogEslkLAjdW6gdFRE0ajBUgzONiHzeZgyZJali+/6Zw3dq018fgbcyKmpsaY
   mgrnfy4rK6ewsPAyto4QQoi3EgkQLtLIyAinTp0itmcPk+96N6W/9EtE/H6O7XieLztd9EdnZvqX
   A5qyMhtlZS4OHWqno+MgY2OH6ehoJxBIksn4KC/PllaOxdqJx0+STqewWt2Ul98EdGI0lqI1GI1l
   aF1IOBzKTQr0nvM829vbee6554DsHIGSEgdms3vRzxdCCHFjkTkIF6n7hz/E9ZWvUNx6kP5f/C06
   7vllDnSdYPXqtfT1/QyDIUNl5TKczmVonSEa7QQ0RmMFhw79kEOHxvF41mG3l2C1nuRjH9tKVVUV
   /f39dHefJBZLkcnUUlxcRFtbO8lkFV7vJE5njOrqImw2+7xJgdFolM7Ofg4efJ33ve89+RLKSimU
   UvljZFKhEEJcfy7VHATpQbgA8XgcvX8/vOc91H384yQ3bOLV/9hP8AN/wIlTBygoWE08nqGiogGt
   DVRXF+HzRejq+iFKjWAy+Th5MkYy+XZstncCU6xYUURx8QiRyGkA/H4/HR1H2LRpBRZLBDBQV1eD
   y3WclSszvO1tK9m48WaamuryKY73799POBympaWXsbEKHI61vPpqD9FoFIPBkA8OIFuYqampLv98
   IYQQYjaZpHi+Tp1i9IEHKB8YwPSZz8DjjzPWd4rMYDHT03EikV4mJvyUl7tZteo2XK46bLYYTU11
   ZDLTHDvWyfh4CKOxFqt1FKu1GK0riUSmKC1tIhweB8BoNPLRj34UpRTNzZb8p/2tW+/I39ATiQRK
   KQoKClBKMT4+Tk/PYD45k9udnYAYCIxKECCEEOK8SIBwvtxuKj7yEYwf/jBYrYyOjqL1aSKRIB0d
   BsDK1FQPPT1xzOYEfX0vccstHwFg5cqVVFdX8+1vv0Is5iaVOkg4HMNiWU0gMEUsNsymTVXA/CGA
   hdIW79ixg1WrVtHQ0ADAvffeS3t795VpByGEENc1mYNwkX7yk59w4sQJVqxYyyuvhCgqslFd7WNk
   ZByte/F4TGzbdjder5exsTFOnjxJSUkJzz57mFOnYqRS46RSdTQ2+nC5imhoSOP3e/M1HCCbmbC5
   uZb+/n6i0Si33347kF2ZMHvYAObWf5j9XOlBEEKIG4MkSlqEyxUgTExMUFxcjMFgoLW1lZdffpnC
   wgbC4WUkEs8CxSSTPpYsqaSmxpe/SYfDYV58cTfLl69iamqS118/Rjhso6FhObW1S0mlkvh82VUP
   J04YiMUmmJoapqysnoaGTH6lwZvd7GUCohBC3LgkUdJV9LWvfY3m5maam5tZt24dShXS319EPJ7A
   YrmPkyd7SKdTLF1ay7FjO6ipaebYsR7GxmBy0sKLL05x+vQk8XgJwSBYLEaGhg6xbJmBxsYahoZC
   HD06gc1WTTweZHR0gMpKOw7H4qojOhwOCQqEEEJcFFnFcAEMBkO++FFbWxvPPPM9Rkcj1NU5qKws
   xOlMUF9fgt1ux2QyMzERZHg4TEFBFWVlzTidK8hkKikoqGbVqk3Y7XHKy51UVhbyzDPPkEqlUcqM
   0eikqOgWlDJzHXf0CCGEuAbJEMMFCAaDOBwOkskkLS29xGIFnDgRQikzbneQRGIEn68Jp3MZp09H
   MRjCeL0GotEqxsenGRmxMTR0EkhQVNSB378Cn68Kny9CU1Md7e3dnDhhJJEwAmA2p2loSC84UVEI
   IYSYTYYYroKZsf0dO37A+vXrqKysxWCowG6fytdEMJsT+P3lbNmyctY8gHogWynx9OkJwuEYBoMC
   NFpX4HA45pRI9vu99Pf3Ulz8xkRDKZ8shBDiSpIehEWavTpA6wzp9BBut2Z01M3hw0+ydesnmJiI
   5HsBZiSTSWKxGG63m2g0yr59B0kmU3g8HgC0Brt94WyIMtFQCCHE+ZJVDItwKQOE9vZuBgeduFwl
   9PTspaDAic8XZ2xMEQyGSaWcpNNDbN++iqqqN3IZ7Nt3kGBwiAceuF9u8kIIIS47SbV8FQ0NtTM5
   OYLNZmflSg/x+CTT01GSSSc/+tF+HnvsMSKRCC0tvWi9grKye2lp6SUajV7tUxdCCCEWRQKERfL7
   vaRSAcLhEE1ND1BS4sTv9xKNxkmnT5PJeFGqiWh0JUr5OH68N5/y2GQqIBBIsnt3qwQJQggh3hIk
   QFgkh8NBc3MtPl+EROIAoVAvu3e3curUKZJJIwaDG7vdhd3uwmKpZng4DEAsFqWtrZdQyM3IiGfB
   noRdu3Zd+Td0nZE2vHjShhdP2vDiSRteOyRAOA8OhwO73cy+fSfo7y+hvb2aF1/sY2oqxtTUFBMT
   IdLpAF5vCeXlLlKpAF1dJ5ietmG1pqmrW4rJ5M9PPpwh/yEunrThxZM2vHjShhdP2vDaIcscz9OB
   A11YrfdTWfnGSgWvdz+p1HGMxjK83hIsliiNjUsB2L27FaXM1NUtxWazEQ7Hr9apCyGEEIsmAcIl
   UFpayp13rsv1DGTw+98ojnTnnetoaekllYoTDsclp4EQQoi3hOt+mePVPgchhBDiSpM8CEIIIYS4
   LGSSohBCCCHmkQBBCCGEEPNIgCCEEEKIeSRAuAKUUr1KqUNKqVal1Gu5fR6l1PNKqQ6l1E+UUq5Z
   x39GKXVCKXVMKXXf1Tvzq0cp9W9KqaBSqm3WvvNuM6XUBqVUW+6xL13p93E1naUNH1VKDeSuxVal
   1PZZj0kbnkEpVa2UelEpdUQpdVgp9cncfrkWF+kcbSjX4iIppQqVUnuVUgeUUkeVUl/I7b+816HW
   Wr4u8xfQA3jO2Pf3wB/ntj8F/G1uewVwACgAaoFOwHC138NVaLMtwDqg7QLbbGYC7mvArbntHwP3
   X+33dpXb8C+B31/gWGnDhduwAlib27YBx4EmuRYvSRvKtXh+7ViU+24C9gCbL/d1KD0IV86ZS07e
   BXwzt/1N4N257QeBJ7TWSa11L9l/2FuvyBleQ7TWPwPGz9h9Pm22SSlVCdi11q/ljvuPWc+57p2l
   DWH+tQjShgvSWp/SWh/IbceAdsCPXIuLdo42BLkWF01rPZXbNANGsv+3L+t1KAHClaGBnUqpfUqp
   j+b2lWutg7ntIFCe2/YBA7OeO8Ab/5ludOfbZmfuDyBtCfDbSqmDSqmvz+qSlDZ8E0qpWrI9MnuR
   a/GCzGrDPbldci0uklLKoJQ6QPZ6e1FrfYTLfB1KgHBl3KG1XgdsB35LKbVl9oM629dzroQUkqzi
   DItoM7Gw/wMsBdYCQ8D/vrqn89aglLIB3wd+R2s9MfsxuRYXJ9eG3yPbhjHkWjwvWuuM1notUAXc
   qZS6+4zHL/l1KAHCFaC1Hsp9HwF+QHbIIKiUqgDIdfsM5w4PANWznl6V2yfOr80Gcvurzth/Q7el
   1npY5wD/lzeGr6QNz0IpVUA2OPhPrfUPc7vlWjwPs9rwv2baUK7FC6O1jgA7gA1c5utQAoTLTClV
   pJSy57aLgfuANuBp4IO5wz4IzPzheRp4SCllVkotBRrITioR59lmWutTQFQptUkppYD3z3rODSn3
   R2TGe8heiyBtuKDce/46cFRr/Y+zHpJrcZHO1oZyLS6eUso7MwSjlLIC24BWLvd1eLVnZl7vX2S7
   0A7kvg4Dn8nt9wA7gQ7gJ4Br1nP+hOykkmPA26/2e7hK7fYEMAgkgJPAhy+kzchG2W25x758td/X
   VW7D3yA7KekQcDD3h6Fc2vCcbbgZyOT+/7bmvu6Xa/Gi23C7XIvn1Yargf25NjwE/FFu/2W9DqUW
   gxBCCCHmkSEGIYQQQswjAYIQQggh5pEAQQghhBDzSIAghBBCiHkkQBBCCCHEPBIgCCGEEGIeCRCE
   uAEopT6klMoopequ9rkIId4aJEAQQgghxDwSIAghhBBiHgkQhLgBKaV2KaV+ppR6m1Jqv1JqUinV
   ppSaVxteKXWzUuoHSqlRpdSUUuqYUurTZxzze0qp40qpuFJqUCn1lZkaJLOOySil/lop9UdKqX6l
   VEwp9SOlVKlSqlIp9X2lVEQp1aeU+uMFzmOpUupxpdSwUuq0Uqp1ofMVQlwaEiAIcWPSQD3wj8D/
   At5LtuTud5VS9TMHKaVuBVrI1hT5XeAdwBeZVUNeKfV5sqV6/xt4J/D3wIeAHbmCMLN9ALgT+Bjw
   28AW4L+Ap4DXyRbteRb4W6XU9lm/oxrYSzYn/e8Cv0A2N/33lVK/cLGNIYSYT2oxCHEDUEp9CPg3
   4CatdbdSahfQDKzQWnfljiklGyT8udb6C7l9u4ElQKPW+vQCr+vJPedxrfVvzNr/MPCfwINa62dy
   +zJki8qs0Fpncvv+N/B7wJ9prT+f22ckW2Rqx8xrKqW+Tjb4WK61Hp/1e34ClGqt112ShhJC5EkP
   ghA3rhMzwQGA1nqEbD35asiWKgduJ3vznxcc5NwGFJDtBZjt20CKbG/BbM/PBAc5x3Pf/3vWeaTJ
   VpqbXbf+fuDHZEvVmma+yFawu1kpZXuzNyuEOD+mq30CQoirZmyBfXGgMLftJvshYuAcr+HJfR+a
   vVNrnVJKhWY9PmP8jJ8TZ9mfnHUeAGVk691/cIFz0EAJEDvHeQohzpMECEKIsxkHMsz9JH+mmSCj
   Emif2Zn7dF/CwkHIhRgFdgN/d5bHh86yXwhxgSRAEEIsSGs9pZR6Gfh1pdRnzzLM0EK2F+Ah4Kez
   9v8K2b8vuy7R6TxHds7E0XMMdwghLiEJEIS4cZ25wmChfX8IvAS05CYUBoA64Gat9Se11uO5/Z9R
   Sk2SXYHQBPw18DOt9Y5LdH5/AbwG7FZK/RPQR3YIZBWwVGv9kYv4PUKIBUiAIMSNQ5+xvdASpjn7
   tNb7lFJ3AJ8FvgJYgF7g32cd86dKqRHgfwC/SXY44JvAZy7gvBY8P631SaXULcCjwOeBUiAEtOV+
   lxDiEpNljkIIIYSYR5Y5CiGEEGIeCRCEEEIIMY8ECEIIIYSYRwIEIYQQQswjAYIQQggh5pEAQQgh
   hBDzSIAghBBCiHkkQBBCCCHEPP8faQngz0IWyFYAAAAASUVORK5CYII=
   "
   >
   </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">
   <h3 id="second-plot">Second plot</h3>
   <p>The dotted black lines form 95% point-wise confidence band around 10 quantile regression estimates (solid black line). The red lines represent OLS regression results along with their 95% confindence interval.</p>
   <p>In most cases, the quantile regression point estimates lie outside the OLS confidence interval, which suggests that the effect of income on food expenditure may not be constant across the distribution.</p>
   </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="kn">from</span> <span class="nn">matplotlib</span> <span class="kn">import</span> <span class="n">rc</span>
   <span class="n">rc</span><span class="p">(</span><span class="s">&#39;text&#39;</span><span class="p">,</span> <span class="n">usetex</span><span class="o">=</span><span class="bp">True</span><span class="p">)</span>
   <span class="n">n</span> <span class="o">=</span> <span class="n">models</span><span class="o">.</span><span class="n">shape</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
   <span class="n">p1</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">q</span><span class="p">,</span> <span class="n">models</span><span class="o">.</span><span class="n">b</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;black&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s">&#39;Quantile Reg.&#39;</span><span class="p">)</span>
   <span class="n">p2</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">q</span><span class="p">,</span> <span class="n">models</span><span class="o">.</span><span class="n">ub</span><span class="p">,</span> <span class="n">linestyle</span><span class="o">=</span><span class="s">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;black&#39;</span><span class="p">)</span>
   <span class="n">p3</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">q</span><span class="p">,</span> <span class="n">models</span><span class="o">.</span><span class="n">lb</span><span class="p">,</span> <span class="n">linestyle</span><span class="o">=</span><span class="s">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;black&#39;</span><span class="p">)</span>
   <span class="n">p4</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">q</span><span class="p">,</span> <span class="p">[</span><span class="n">ols</span><span class="p">[</span><span class="s">&#39;b&#39;</span><span class="p">]]</span> <span class="o">*</span> <span class="n">n</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;red&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s">&#39;OLS&#39;</span><span class="p">)</span>
   <span class="n">p5</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">q</span><span class="p">,</span> <span class="p">[</span><span class="n">ols</span><span class="p">[</span><span class="s">&#39;lb&#39;</span><span class="p">]]</span> <span class="o">*</span> <span class="n">n</span><span class="p">,</span> <span class="n">linestyle</span><span class="o">=</span><span class="s">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;red&#39;</span><span class="p">)</span>
   <span class="n">p6</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">plot</span><span class="p">(</span><span class="n">models</span><span class="o">.</span><span class="n">q</span><span class="p">,</span> <span class="p">[</span><span class="n">ols</span><span class="p">[</span><span class="s">&#39;ub&#39;</span><span class="p">]]</span> <span class="o">*</span> <span class="n">n</span><span class="p">,</span> <span class="n">linestyle</span><span class="o">=</span><span class="s">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s">&#39;red&#39;</span><span class="p">)</span>
   <span class="n">plt</span><span class="o">.</span><span class="n">ylabel</span><span class="p">(</span><span class="s">r&#39;\beta_\mbox{income}&#39;</span><span class="p">)</span>
   <span class="n">plt</span><span class="o">.</span><span class="n">xlabel</span><span class="p">(</span><span class="s">&#39;Quantiles of the conditional food expenditure distribution&#39;</span><span class="p">)</span>
   <span class="n">plt</span><span class="o">.</span><span class="n">legend</span><span class="p">()</span>
   <span class="n">plt</span><span class="o">.</span><span class="n">show</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_stderr output_text">
   <pre>
   /usr/lib/python2.7/dist-packages/IPython/core/formatters.py:239: FormatterWarning: Exception in image/png formatter: LaTeX was not able to process the following string:
   &apos;lp&apos;
   Here is the full report generated by LaTeX: 
   
   
     FormatterWarning,
   
   </pre>
   </div>
   </div>
   
   <div class="output_area"><div class="prompt"></div>
   
   
   <div class="output_text output_subarea ">
   <pre>
   &lt;matplotlib.figure.Figure at 0x7fb596d15250&gt;
   </pre>
   </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>