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="s1">&#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:                Mon, 04 Apr 2016   No. Observations:                  235
   Time:                        22:56:03   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="s1">&#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="s1">&#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="s1">&#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="s1">&#39;q&#39;</span><span class="p">,</span> <span class="s1">&#39;a&#39;</span><span class="p">,</span> <span class="s1">&#39;b&#39;</span><span class="p">,</span><span class="s1">&#39;lb&#39;</span><span class="p">,</span><span class="s1">&#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="s1">&#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="s1">&#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="s1">&#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="s1">&#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_stderr output_text">
   <pre>
   /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)
   /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 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.880099  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.033434  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>
   </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="s1">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#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="s1">&#39;a&#39;</span><span class="p">],</span> <span class="n">ols</span><span class="p">[</span><span class="s1">&#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="s1">&#39;red&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#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="s1">&#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="s1">&#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,iVBORw0KGgoAAAANSUhEUgAAAgoAAAFWCAYAAAAIZbVEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz
   AAALEgAACxIB0t1+/AAAIABJREFUeJzsvXl8FFW+9/8+2fewJIFASCAssoMoCqgIwyKKG+ooqKjc
   GfS5PjrLvdefM+N91HFmvM6i3lnueGcURFxwGQRBQRAEBNmEIIQlCRCyQ9ZO0umkk17O74/q1HRC
   d9JLZUHO+/XqF12nqk6dqgbOp77nuwgpJQqFQqFQKBSeCOnpASgUCoVCoei9KKGgUCgUCoXCK0oo
   KBQKhUKh8IoSCgqFQqFQKLyihIJCoVAoFAqvKKGgUCgUCoXCK90uFIQQaUKIL4UQJ4UQ2UKIH7na
   +wohtgohcoUQW4QQiW7n/EkIcVoI8a0QYrJb+8NCiDzXOQ91970oFAqFQvFdR3R3HgUhxEBgoJTy
   WyFEHHAYuANYBlRLKX8nhHga6Cul/JkQ4mbgCSnlQiHEtcAfpZTThBB9gUPAFEC4+pkipazr1htS
   KBQKheI7TLdbFKSUF6SU37q+NwCngDQ0sfCW67C3XNu4/lztOv4AkCiEGADcBGyVUtZJKWuBrcCC
   brsRhUKhUCguA3rUR0EIMRSYDOwHBkgpy0ETE0CK67DBQLHbaSWutvbtpa42hUKhUCgUBhHWUxd2
   LTv8A/ixlLJBCOFtDUR42JYe2nG1t7+OylGtUCgUissOKaWnedJvesSiIIQIQxMJb0spP3E1l7uW
   FFr9GCpc7SXAELfT04AyV3u6h/aLkFJeEp/nnnuux8dwKX/U81PPsDd81DP87j3D4uJiHA6H4f1+
   9tlnZGdnd8mYjaSnlh5WAiellH90a9sAPOL6/gjwiVv7QwBCiGlArdSWKLYA84QQiS7HxnmuNoVC
   oVAoAsZms3HixAl9Oy0tjZAQ46fLuXPnMm7cOEP6klLy/vvvU11dbUh/7vREeOR1wAPA94QQR4QQ
   WUKIBcBv0Sb+XGAO8BKAlHITcE4IcQb4G/C4q90E/Aot8uEA8EupOTUqFAqFQhEwDoeD4uJinE6n
   of06nU4++OADLBYLABEREQhhyOoAQgjmzJlDv379DOnPnW73UZBSfg2Eetk918s5T3hpXwWsMmRg
   vYBZs2b19BAuadTzCx71DINHPcPg6YlnKKXEarUSHR1NVFQUCxYYH0QXEhLC9OnTiYmJMaS/uro6
   Dh06xJw5cwBITk42pN/2dHsehe5GCCG/6/eoUCgUiuDIyckhLy+P22+/3dB+GxoaOHv2LJMmTTK0
   XwC73U5ubi5jx469yDIhhEAa5MzYY1EPCoVCofjuMXToUAoLC3t6GJcNGRkZFBQUdOk1lEVBoVAo
   FIbhepPt6WFcNnh73kZaFFRRKIVCoVAoFF5RQkGhUCgUCoVXlFBQKBQKhULhFSUUFAqFQqFQeEUJ
   BYVCoVAoFF5RQkGhUCgUlxWrVq1i4sSJxMbGMmjQIB5//HHq6uoA+OUvf8nSpUs9nrdnzx6uu+46
   +vTpQ1JSEjfccAOHDx/uzqH3CEooKBQKheKy4eWXX+bnP/85L7/8MvX19ezfv5/CwkLmz5+P3W4H
   8JhW2Ww2c9ttt/HjH/8Yk8lEaWkpzz33HJGRkd19C92OyqOgUCgUCsPozXkUzGYzgwYNYtWqVdx9
   9916u8ViYfjw4bz00ksUFhZy9uxZVq9e3ebcw4cPM2/ePGpqarp72B2i8igoFAqFQmEQe/fupbm5
   mUWLFrVpj42NZcGCBXzxxRdezx01ahShoaE88sgjfP7559TWXj41CJVQUCgUCkX3IkTwnwCoqqoi
   KSnJY8no1NRUqqqqvJ4bHx/Pnj17CAkJ4dFHHyUlJYU77riDysrKgMZyKaGEgkKhUCi6FymD/wRA
   UlISVVVVHstHnz9/nqSkpA7Pv+KKK1i5ciVFRUUcP36csrIyfvKTnwQ0lq7C6NLYoISCQqFQKC4T
   pk+fTmRkJB9//HGbdovFwubNm/Vyzb4watQoHnnkEY4fP270MAPG6XTy+uuvY7FYDO1XCQWFQqFQ
   XBYkJCTw7LPP8uSTT7JlyxbsdjsFBQXce++9pKen62GRDoeD5uZm/dPS0kJubi6vvPIKpaWlABQX
   F7NmzRqmT5/ek7fUhpCQEB544AFiY2ON7dfQ3hQKhUKh6MU89dRTvPjii/zHf/wHiYmJTJ8+nYyM
   DLZt20Z4eDgA77//PjExMcTExBAdHc2IESNISEjgwIEDXHvttcTHxzNjxgwmTpzIH/7whx6+I/j0
   00/173FxcYb3r8IjFQqFQmEYvTk88ruIEAKTyUSfPn0uajcqPFIJBYVCoVAYhhIK3YvKo6BQKBQK
   haJHUUJBoVAoFAqFV5RQUCgUCoVC4RUlFBQKhUKhUHhFCQWFQqFQKBRe6XahIIRYIYQoF0Icc2ub
   JITYJ4Q4IoQ4KISY6rbvT0KI00KIb4UQk93aHxZC5AkhcoUQD3X3fSgUCoVCcTnQExaFN4Gb2rX9
   DnhOSnkl8JxrGyHELcBwKeVI4DHgf13tfYFnganAtcBzQojE7hm+QqFQKBSXD90uFKSUewBTu2Yn
   0DrR9wFKXd9vB1a7zjsAJAohBqAJja1SyjopZS2wFVjQ1WNXKBQKxaVNXl4eU6ZMITExkb/85S+G
   9Tts2DC+/PJLw/rrTYT19ABc/BTYIoR4GRDADFf7YKDY7bgSV1v79lJXm0KhUCgUXvnd737H7Nmz
   ycrK6umhXDL0FqHwr8CPpZTrhRD3ACuBeWiiwR0BSA/tuNo98vzzz+vfZ82axaxZs4IcrkKhUCiM
   pqamhuLiakJCBBkZKSQkJBh+jcLCQpYsWWJ4vz3Nzp072blzZ9d0LqXs9g+QARxz265tt7/W9ef/
   Ave5tecAA4DFwP+6tbc5rl1fUqFQKBTdg7f/cx0OhywqKpXHjp2W584VSbvd3mZ/dXW1/OyzY3LH
   DpPcvr1afv75UVlfX39RPzabTTY2NkqHw+H32L73ve/J0NBQGRUVJePj4+XRo0fl0qVLZXJyshw6
   dKj89a9/rR/rdDrlr371K5mRkSEHDBggH374YVlXV6fvX716tczIyJBJSUnyN7/5jRw6dKjcvn27
   fu5//dd/yeHDh8ukpCR53333SZPJJKWU8oMPPpCZmZnSbDZLKaXctGmTHDhwoKyqqpJSSimEkH/6
   059kZmamTE5Olk899VSH9+TtebvaDZmzeyo8UtDWKlAqhLgRQAgxBzjtat8APORqn4YmIMqBLcA8
   IUSiy7FxnqtNoVAoFL2QnJx8srOtVFcnceqUk2PHTrepUVBcXE1MzBDi4/uQmNiP0NBUzp+vbtPH
   hQsVfPnlcXbtyufrr49jsVj8GsP27du54YYb+Otf/0p9fT1/+MMfMJvNFBQUsHPnTlavXs2bb74J
   wJtvvsnq1avZtWsX+fn5mM1mnnjiCQBOnjzJ448/zrvvvktZWRnV1dV6+WmAP/7xj2zYsIHdu3dT
   VlZG3759efzxxwG49957mTFjBj/60Y+oqanhhz/8IStXrqR///76+evXrycrK4usrCw++eQTVq5c
   6fM9FhQUcOzYsc4P9AejFIevH+A9oAxoBoqAZWg+CYeAI8A+4Eq34/8CnAGOAlPc2h9BExR5wEMd
   XK9DNaZQKBQK4/D0f25zc7PcvPmo/OYbpzx0SMpDh6TctOmEbGho0I85ejRPfvlljb5/69ZymZd3
   Tt9vsVjk5s1H5b59zfLQISl37DDJPXuO+T2+WbNmyRUrVkiHwyEjIyNlTk6Ovu9vf/ubnD17tpRS
   yjlz5sjXXntN35ebmysjIiKkw+GQL7zwglyyZEmbsUVEROgWhTFjxsgvv/xS319WVibDw8N1K0ht
   ba1MT0+XEyZMkP/6r//aZnxCCLl161Z9+69//aucO3eu1/tp/7wrKirkuXPnDLUodLuPgpTyfi+7
   rvZy/BNe2lcBq4wZlUKhUCi6G+lmUcjISKGsrJDaWgdOpxM4z6BBI/T9VqsViCc8PAKA+Pg+VFae
   w+FwEBoa6ve1q6qqsNlspKenu40hQ7cMlJWVkZGR0Waf3W6nvLycsrIyhgwZou+LiYlpYxEoLCxk
   0aJFhISE6PcZHh5OeXk5qampJCYm8v3vf59XX32Vjz/++KKxpaWltbluWVlZh/eydetWbrzxRiIj
   I0lOTiY5OdnPp9ExKjOjQqFQKLqUiIgI0tKiqag4R0NDHRUVRaSkCGJjY/VjEhMTmTFjKIMHm8nI
   aGTGjBFt9kdGRiJlA3a7HQCLxUxsbGhAIgEgKSmJsLAwCgsL9bbCwkIGD9YC6AYNGnTRvrCwMAYM
   GEBqairFxf8MvGtsbKS6+p/LJOnp6WzevJmamhpqamowmUxYLBZSU1MB+Pbbb1m5ciVLlizhySef
   vGhs7n0XFRUxaNCgDu+ltd+uQgkFhUKhUARN6wTujbFjhzN+fCR9+lRwxRUwadJIhGgbwJaQkMDo
   0cMYNWpoG5EAEBsby4QJSZhMJ6iszMHhyGfy5GEBjzckJIR7772XZ555hoaGBgoLC3n11VdZunQp
   AEuWLOHVV1+loKCAhoYGnnnmGRYvXkxISAj33HMPn376KXv37sVms/Hss8+2sY489thj/OIXv6Co
   qAiAyspKNmzYAGiWkaVLl/LSSy+xcuVKysrKeO2119qM7fe//z21tbUUFxfzxz/+kcWLF3d4LxMm
   TCAyMjLgZ9EpRq1h9NYPykdBoVAouhSr1Sr/53/+R7a0tHj1wjfyWmazWdpstoDOnz17tlyxYoWU
   UkqTySQffPBBmZycLNPT0z1GPQwZMkSmpKTIhx56SNbW1ur7V69eLdPT02VSUpJ88cUX5bBhw9pE
   Pbz66qvyiiuukAkJCXLEiBHymWeekVJK+dOf/lTecsstej9Hjx6V/fv3l2fOnJFSaj4Kf/7zn2Vm
   ZqZMSkqSTz31lHQ6nV7vx9vzxkAfBSGl1/QD3wmEEPK7fo8KhULR01itVqKiohBCoP7PDZyQkBDO
   nDlDZmamT8d7e96udk85h/wfkxGdKBQKheLyoqKiok2Cn6ioqJ4bjKJLUUJBoVAoFH6TkJDQqZOd
   wn/a+234g8PhoLGx0cDRaKilB4VCoVD4RF1dHQ6Hg379+nk9Ri09dC/uz/vAgQO0tLRwww03GLr0
   0FtqPSgUCoWil3P27FlCQ0M7FAqKnmPq1Kl67gYjURYFhUKhUHjF34RGyqLQvQghaGxsJDo6+qJ2
   5cyoUCgUii7nww8/pKSkpKeHoeiA9iLBaJRFQaFQKBReaWpq8msiGjp0aJuMhoquJSMjg4KCgova
   lUVBoVAoFF2CzWZj8+bNeqZFf99WCwoKPCbt2bJlC1VVVUEn/3E6nZw/f77LkvSdPXvWkHF6++zZ
   s8fv8ZeVlbW5f/d9nkSC0SihoFAoFAqdsLAw0tLSggrTa8U9VG/+/PltCif5Q35+PtnZ2YD2pjxw
   4MCgx9aKzWZj3759ul9FZmZmwOP0hdTU1IvSU3eE0+lk+/bt+rM04nfxFyUUFAqF4jLHbrfrFQqF
   EEyYMCHgYkutWCwWVq9e7aoEGRxxcXH06dMn6H48ERoais1mw+FwdEn/dXV1fPHFF/p2ZmYm8fHx
   HZ7T2NhIZWUloGVqfPDBB4mJiemS8fmCEgoKhUJxmVNdXc3hw4cN7TM2Npbly5cHHK53/PhxWlpa
   AEhJSWlT1jlYqqqq9AqNISEhzJw5k7CwrskWEBcXx6BBg/yKBCkoKODMmTNdMp5AUEJBoVAoLkOk
   lPpb9IABA7jtttuC7jMvL4+tW7fq28FYJaqrq7FYLEGPyRN1dXWYTKYu6Rvg4MGDukNnaGgo48aN
   63DJoNU3olVMjB07lunTp3fZ+PxFRT0oFArFZciBAwdobm5m5syZhvVptVppbm4mMTHR73OllFRU
   VDBgwADDxuPOuXPnyMjI6JKERO0pKioiMTHR5+cgpWT9+vXMmzePuLg4Q8agoh4UCoVCERSTJ09m
   xowZQfeTlZVFVVUVoBWGCkQkANTX17N9+/YuSdYkpeTYsWPU19cb3jdo/hifffaZPvb09PROn0Nd
   XR2lpaWANqkvWrTIMJFgNEooKBQKxWXCnj17qK6uBiAyMtKQdfno6OiAPfGdTqfuh5CYmMj9999v
   mFe/0+nUHQKFENxxxx1d5hAZExPDsGHD/DqnurpaFwq9HSUUFAqF4jIhOTmZyMjIoPqQUuqOgABj
   xowJOJzwwIEDHDhwIKjxeOPChQvs2bOnS/oG+Pbbb8nLywM0ITJ27NhO/RCys7N1v5DMzEyuueaa
   LhufkSgfBYVCofgOU1JSQlpammH92e12/vGPf3DnnXcSFRXl9/lSSn1CtdvthIaGGmZFaGpqIjw8
   vMsiGNw5f/48kZGRPhfIklLyxRdfMH369E7DI41A+SgoFAqFolPsdjtfffUVzc3NQfdls9kALSHT
   4sWLAxYJ7777rr78ERYWZmgCoS+++IJz584Z1p87zc3NbNiwQbcIpKamdioSqqurOXv2LKBN3PPn
   z+8WkWA0yqKgUCgU3zGcTqeh3v3FxcXs2bOHJUuWBN1XdXU1/fr1M0wgtLS0EBERAbS1VhiNlJLc
   3FxGjRrl87MtKyvjwoULTJkypUvG1BHKoqBQKBQKj1y4cIF3333X0D7T0tK4++67Azq3traWLVu2
   6BEB/fv3N3Sp4Y033tDf8o0WCSdOnGiTOnr06NEdigQpJQcOHNCtL4MGDeoRkWA03S4UhBArhBDl
   Qohj7dqfFELkCCGyhRAvubX/XAhxWghxSggx3619gev4PCHE0915DwqFQtFbGTBgAIsWLQq6n2++
   +abNJNn61u4vcXFxZGRkBD0ed1rTQkdHR7N8+fKg0017IyUlhdTUVJ+PF0LgcDgMWerpTXT70oMQ
   4nqgAVgtpZzoapsF/AK4RUppF0IkSSmrhBBjgPeAqUAasA0YCQggD5gDlAHfAIullDkerqeWHhQK
   xXea/Px8rFYrY8eONazPqqoqoqKiAortLyoqIjQ0lMGDBxs2nla++eYbLBYLs2bNMrxvu93Op59+
   ysKFCwkPD/fpnIqKCioqKhg/frzh4wkGI5ceut41tB1Syj1CiPby8l+Bl6SUdtcxVa72O4D3Xe0F
   QojTwDVoQuG0lLIQQAjxvuvYi4SCQqFQfNeJiYkJ2tNfSsnhw4eZNGkS4eHhJCUlBdyX1Wo19C3f
   3fdg4sSJXWZBCAsLY+zYsX75d4SGhnZJkqjeRG/xURgFzBRC7BdC7BBCXOVqHwwUux1X6mpr317i
   alMoFIrLgsrKSn0tfODAgaSnpwfVnxBCT8EcCCUlJfqEOWrUKIYPHx7UeFqRUrJ69Wo9q6JRiaJa
   OX36NIcOHdK3R40a1aEQkVKya9cumpqaAM3nYsKECYaNJ2h27YI5cwztststCl4IA/pIKacJIaYC
   HwGZaJaD9kg8Cxyvku7555/Xv8+aNatLTFYKhULRnRw8eJCJEycGVVVRSkl1dbVuPbj++usD7mff
   vn3Mnz8/4BTO3hBCcOutt3ZZWGFSUhIJCQl+jadPnz6GlM82DCnZ+cor7PzTn8BsBgPrd0APhUe6
   lh42uvkobEJbevjKtX0amAYsB5BSvuRq/xx4Dk1APC+lXOBq/5l2mPyth2spHwWFQvGdwOiwx9ra
   WjZu3MiDDz4YUMRAY2MjMTExho2nldLSUrKzs1mwYIHhfUsp+eSTT7jpppuIjo726ZwLFy5QVFTU
   +zIpSglbt8ILL0BVFfznf8KSJeDKT3Gph0cK2loL1qM5JiKEGAVESCmrgQ3AfUKICCHEMGAEcBDN
   eXGEECJDCBEBLHYdq1AoFN9JHA4Hr7/+um7yDobWt+E+ffq0EQl2ux2LxYLdbu+0j/Lycv7xj38E
   PRZPJCcnM2nSpC7pWwjBhAkTfHZWBIiNjTXcUhIUUsKmTTB9Ovz0p/DEE3DyJCxdCl2QlbInoh7e
   A2YB/YFyNAvB28CbwGSgGfh3KeUu1/E/B34A2IAfSym3utoXAH9EEzsrWq0OHq6nLAoKheI7gdls
   DtoEf+TIEaqrq5k7d26b9poaE0eOFGGzRRAe3sKUKen07du3zTGt/5e2CguHw2GYY+Hu3bsZOXIk
   AwcONKQ/dwoKCigpKfFraWXr1q1MmzbNr2WJLkdK+PRTzYJgtcKzz8Ldd4MHK5ORFgWVmVGhUCh6
   KXV1dRw/fpzrrrvOsD5bqzW650Ww2+3s2nWCqKhRREVFY7U2YbXmceON49o4Du7atYu4uDiuuuqq
   i/oNlnPnzpGcnNwlpZbNZjNms5lBgwb5fM7p06dJS0vzeXmiS3E6YcMGTSA4nZpAuPNOjwKhqKiI
   QYMGER4efskvPSgUCoWiE6KiooiOjg46/O7LL7+koqIC0ARC++RJzc3N2GwRREVFu64bjc0WcVEE
   xNVXX23YkkBjYyO7d+/W723YsGGGiQQpJRs2bNAjJeLj4zsVCefPn2fHjh369siRI3teJDidsHYt
   XHkl/OpX8NxzcOQI3HWXR5EAkJ2djclkMnQYSigoFApFL8JqtVJbWwtooYBTpkwJOjXxsGHDOlyy
   iIyMJDy8Bau1yTWGJsLDWwgPD2fdunU0NjYC2lq9UaGJkZGRREREdEkOAiEEEydO9Gui79u3L5mZ
   mYaPJSAcDvjgA5g4EX77W3jxRTh0CO64A9r9XSgvL+fkyZP69sKFC0lOTjZ0OEooKBQKRS8iJyeH
   EydOBNWHw+Fo08ewYcM6nDTDwsKYMiUdqzWPyspTWK15TJmSTkREBBMnTiQyMjKo8bRSXFxMSUkJ
   oCUquvbaaw2L4igtLeXLL7/Ut4cOHdqpw+Knn35KZWUloFlvjE417TcOB7z3HkyYAP/93/CHP8CB
   A7Bw4UUCoRWXL0KXDkv5KCgUCkUPY3TVQ5vNxtatW5k/f75f3v12u52Kioouq3iYl5dHSEgII0aM
   MLxvq9VKVVUVaWlpPp9TUlLCgAED/HpGXYLdrgmE3/wGkpO1JYa5cz2KA6fTyUcffcSdd97ZoYD7
   LoRHKhQKhcLFhg0bOHv2bFB9SClpaGgAIDw83K96Ba2EhYURHx+vV2MMFiklp06dapOx0UiRsGnT
   pjYWgc5EQllZGZ9++qm+nZaW1rMiwWaDN9+E0aNhxQp47TXYvRvmzWsjEqSUekhrSEgI06ZN69Zx
   K4uCQqFQ9DD19fXExcUFZYY/e/Ys2dnZ3HnnnX6fW1FRQXR0tOHZD51OJ5s2bWL27NnExsYa2jdo
   Hv4DBw70ubKlzWajurq6S0Iw/aKlBd56S/M9yMzUohhuvNHr4bt37yY0NJQZM2b4fIkeCY8UQgwG
   /h2YCfQDbpdSHhdC/ATYJ6U8YMSAjEYJBYVC0dtwOBzs3r2b6667Lqg3Q/e8BlJKpJQBiY29e/eS
   kpJiyNt+S0sLJpOJAQMGBN1XeyoqKjh8+DA333yzz+esX7+eqVOndkklS79pboaVK+GllzQrwv/7
   f+Aht4OUEpPJRL9+/QBtWSUiIsKv37bblx6EEOOAbGApWlnnDKBVwmUAPzZiMAqFQvFdoLMMhyEh
   IcTGxhoS9tjqtCiE8GsiMZvN+vcZM2YYtiRQVlZGdna2IX21p2/fvowbN67T49yf63XXXdfzFgSr
   Ff7yFxgxAjZu1CIatmzxKBIALBYL69ev1+8jKirK0NTd/uKTRcFVYyEeuAmwAi3A1VLKLCHE94Hf
   Sil7SVxJW5RFQaFQdCfeMhxKKamqqjI0dK2uri6gkEWHw8Ebb7zBQw895HMIod1up7m52WP1xtra
   WuLj47uk/PO2bdsYPXq0z06K58+fZ/fu3dx7772Gj8VvGhvh73+H3/0OrrpKW2KYOtXjoWfPniUl
   JUVf/gnWwbUnnBmvRyva1MDFVRrLgR6WawqFQtHz2O12jhwpIipqFMnJY4iKGkVWVhF2u52amhq2
   bdsWlBXB6XSybds2PbtiYmKiXyKh1SEuNDSURx991GeRUFNjYteuE+zeXcSuXScuSuizfft2zp8/
   7/M4/GHMmDGdiiv3ZzpgwADmzZvXJWPxGYsFXn4Zhg+Hr77S0i5v3OhVJIBWNrzVGRXwWyRIKVmz
   Zg11dXUBD9sbvgqFjuppJgHBVylRKBSKSxxPGQ5bWsJpbm6mf//+LF68OKi3xJCQEL0ktL8UFBSw
   du1afdvXcXgTP61ZDwHuuusuv8ISO8JkMvHJJ5/ok//gwYM7zePw8ccfU1hYCGjPqH2Nim6joUGz
   HgwfDvv3w+efw8cfg4dQ08rKSvbu3atvT5s2jdTUVL8uJ6XUs2cKIZg9e3aXlOP2VSgcBJZ52Xcv
   8LUxw1EoFIpLl/YZDgsKDlNW9q0+0QUiEux2uz4JAkyePNlnL3930tPTWbhwod/neRI/DQ021qxZ
   c1GRKCNITExk8uTJnfbZah0BmDdvHunp6YaNwW/q6+G//ksTCEeOwLZt8NFH0EG669jYWN1ZMVAO
   HDjA/v379e2BAwd2iS+Drz3+CrhNCLEVzaFRAnOFEG8Bi4DfGD4yhUKhuMRon+EwLk5yyy2zg0p7
   XF9fz9GjRwM6d9euXZw5cwbQ3rRjYmL87qNV/DQ1NeJ0OrBam4iLC+eRRx4xTCB89dVXbcbZWYbE
   8vJy3nnnHX07ISHBULHiM7W1Wg2G4cPh+HHYuRPWrIHx4y86VErJO++8o6fnjomJYfTo0X5dzuFw
   kJ+fr29PnTqVmTNnBnULvuBPeORC4L+B4W7NBcD/lVJuNn5oxqCcGRUKRXfy7bffkpqaSlxcnEfH
   P19wOp3Y7faALAfulJeXk5iYSFRUVFD9mEwm1q3bQkhIDEOGZHgsQR0MFy5cIDExsUOfCYfDQUhI
   iB4K2tjY2CW5GXzCZII//lGLZLj1VvjFL2DUqIsOk1LS0tKiW5Sqqqro379/wKKmubmZjRs3ctdd
   d3VqOejRMtNCiBFAClAtpcw1YhBdiRIKCoWiOzl58iQpKSkB+xKAltfA6XRyvZfwOW/YbDa+/vpr
   Zs6caYiRtYfIAAAgAElEQVQJ2t3z3mKxIKUkKioq6MJQDQ0NbN68mXvuucfnSXPt2rVMnjyZ4cOH
   d35wV1FdDa++qmVQvPNO+PnPtZBHLxw6dIi6ujrmzJkT8CUPHz5MWlqa33kpulUoCCEigP3Az6SU
   W424aHeihIJCoehqqqqqghIG7bHb7YSGhgbk+f7NN99w5ZVXBp3itzWE8uGHHw7aItEeKSVFRUWd
   LjHYbDb9PqxWq+Hj8JnKSnjlFS3U8e67NYEwbNhFh0kpKSkpYciQIUBbK0ig5OXlkZSU5Lc/Q7eG
   R0opW4BhgOfMIQqFQnEZ43Q6+fTTT9uEtgXCunXrqKioADRfB18nF4vFQlFREaBNDtdcc40hdQBC
   Q0NZvHixYZPz/v37OX78OKCNszORUFlZydtvv61v94hIKC+Hp56CK67Q/BGysjSx4EEkgCbwdu/e
   rSfaCkTslZSUsG7dOn171KhRQTs9BouvCZc+BPKllD/r+iEZi7IoKBSKrsDdLG9E9ccLFy6QnJzs
   d9KikpISzp07xw033BDU9QHOnDnDuXPnuiQPQVVVFdHR0R36FbS0tBAaGqo/A3eLQrdy/jz8/vew
   ahXcfz88/TS4rATtycvLIyEhIajsjxaLRX8udrudhoYG+vTpE3B/0DMJl/4MLBFC/EEIcb0QYrgQ
   ItP9Y8RgFAqF4lLAZDLx9ttvBxUe2NzczP79+/U+Bg4c6LNIMJvN+ltrWlqaISKhta+pHSQF8ger
   1cqaNWv0cSYlJXXqfPjZZ5+18ervdpFQWgo/+hGMGwdOJ2Rnaw6LXkQCaMsLwVTblFLy7rvv6nkp
   wsLCghYJRuOrRcE94ZLHE6SUxufuNABlUVAoFEbTvmhPINjtdvbv38+MGTP8djzcsGED48ePJzMz
   +He0zZs3M3XqVEN9LFopLi4mLS2tQyHl/jbtcDi6JA10pxQXa4Wa1qyBZcu05QYvFoKamhr27t3L
   rbfeGvDlqqqqsNvtuhXCCItUe7o96kEI8XBnx0gp3zJiQEajhIJCoTCC8+fPU1dX53fsuzstLS2Y
   zWb69+/v97l2u12PNjByYiksLCQ1NTXoUEyAI0eO4HA4uPrqq306vqamhvXr17Ns2bKeyYNQUKAl
   SvroI1i+HP793yEl5aLDnE4nQgiEEDgcDgoLC4MSaSdPnsTpdDLeQ74Fo+jR8MhLDSUUFAqFEVy4
   cAGTycSYMWMC7iM3N5fS0lK+973v+XVec3Mzb7zxBo8++mjQ5vi6ujqOHDnCrFmzgurHEyaTifDw
   cOLi4tq0uxeUstlshIWF6ffhdDq7vzJifj68+CKsWwePPQb/9m/QgUVlzZo1zJw5M+BS1Y2Njeze
   vZv58+d3myBSQsEPlFBQKBSBYjabiYqKCmpyNiJEDowLD7TZbJw6dYoJEyYEPSa73c4HH3zAPffc
   47UeQ/tqmrW1Z5gwYQKjPCQo6nJOn9YEwsaN8Pjj8JOfgIflI6fTSUNDAwkJCYC2PBITE+PX83L3
   X5FSkp2dzfjx47tNFPXE0sPKTg6RUsofGDEgo1FCQaFQBMqmTZsYMWJEUJPaJ598wpgxY/zu48yZ
   MxQUFDB37tyAr91Kbm4uCQkJfhcd8oWSkhIGDx7scRK12+3s3HkcKVPo128QVmsTTU25zJo1Puik
   TX6Rmwu/+Q1s3gxPPqk5LHbgMJiTk0N+fj633HJLwJfcsGEDY8aMYeTIkQH3EQw9IRQKuNiJsR8Q
   D9QCtVJKnxZshBArgFuBcinlxHb7/gP4HZAkpaxxtf0JuBmwAI9IKb91tT8MPOMa12+klKu9XE8J
   BYVC4TNGhz1arVYiIyP97qepqQmr1WpIquTc3Fzi4+MZNGhQ0H2dOHGC2tparrvuuk6PtVgsbNt2
   iuLiE0yb9hBCCCorT3HDDendk3755En49a+1Ik0//jE88QQkJl50mNPpJCcnhzFjxugWAH9/L6fT
   idlsJtHVv9lsJjY2tvuXVVx0e3iklHKolHJYu08iMAu4ANztxzXfBG5q3yiESAPmAoVubTcDw6WU
   I4HHgP91tfcFngWmAtcCzwkhLv71FQqFwg+klKxevVoPVQu02uPHH3+sl/+NioryuZ+9e/dSXV0N
   QHR0dMAiweFwcPToUd38fcUVVxgiEkCrQjlhwoQOj6mvr6epqYnIyEji4sKZPPn7CCGwWpsID2/p
   tGx00GRnw333wezZWgXHs2fhmWc8igTQfuf8/Hyampr0bX8pLCxkx44d+nZ8fHyPiQSjCeoupJRf
   Aa+i5Vnw9Zw9gMnDrleBp9q13QGsdp13AEgUQgxAExpbpZR1UspaYCuwwP87UCgUin8ihGDhwoXE
   x8cH3EdYWBiTJk0KyK+hf//+huUOOH/+PDabLeh+nE4n7733HhaLBdAmwNa1e28cOHCA4uJivZpm
   c/NpKitPYbXmMWVKetctO3z7rZZied48mDpVEwhPPw0efs+8vDxOnz4NaL/7rbfe6ld1TSklR44c
   0ctdDx06lDvuuMOY++hlGPFr5QNXBtOBEOI2oFhKmd1OyQ0Git22S1xt7dtLXW0KhULhF01NTRw/
   flxPNBRIPoGGhgZKSkr00ElfCxfZ7XZyc3MZN24coL35B0pDQwMNDQ164qYFC4x5dwoJCWHWrFkd
   TqJOp5Pz58/rUQHumR379u3LjTfG61EPXSISDh/Wyj0fPKjlQHj7behk0vfXObE9Qghqampoamoi
   Nja2Z8I7u4mgfjEhRBjwCNoEHmgf0Wi+Bp5yhrZ/8gLNJ8HTL+LVEeH555/Xv8+aNatLwoIUCsWl
   SUhICFarNSh/BJvNRmVlpd85FpxOJ4WFhYwePTroRENlZWXU1NQElUq4ldOnT1NUVKRXPexs2aI1
   /O++++7z+AzDwsK6RiAcPAgvvKBZEp5+WkuY5KVUtcViYe3atTz44IOEhISQlpbm9+VOnDiBw+Fg
   4kTNvS6YqpBGs3PnTnbu3NklffvqzPilh+YIYBTQH/g/UsrXfb6oEBnARinlRCHEeGAb0IgmANLQ
   LATXAC8AO6SUH7jOywFuBGYDs6SU/8fV/r/ux7W7lnJmVCgUbXA4HFgslk5N6B3RKi6ivUxM3mhp
   aaGxsdGQNL0VFRUkJSUZvhbe2NhIS0tLh2OsqqoiIiIiqGcYMPv2aQLhxAn42c/gX/4FPISO2u12
   hBC6CCsvL/e7XLN7nofKykqEEF2SxdJoeqLWQwjaJO7+MQMfA3P8EQkuWvtASnlcSjlQSpkppRyG
   Zp24UkpZAWwAHgIQQkxDi64oB7YA84QQiS7HxnmuNoVCoeiU3Nxc9u3bF1Qf+/fvJycnJ6BrHz16
   NKhrt7Jr1y5qamqC7kdKyfvvv09tbS2gmeU7EzK5ubmUl5cHfW2/2LNH8z9YsgQWLdLyIjz+uEeR
   ALBx40bOnj2rb/srEhoaGvjb3/6mO4UmJyf3epHQ2NhoeJ/dnnBJCPEeWrREf6AceE5K+abb/nzg
   arfwyL+gOSpagGVSyixX+yP8Mzzy1yo8UqFQ+EqgxZwCDZ1saWkhPDw86HVsKSVms7lL3uIvXLhA
   SkqKV+uE3W6noKCAESNGGH7tTtm5U7MgFBRo0QtLl4KHlNN2u53q6mpdELinvfaV2tpaYmJi9JTW
   jY2Nfjk59iQOh4PXXnuN5cuXExUV1b0WBSHEs0IIj4tUQohUIcSzvl5QSnm/lHKQlDJSSpnuLhJc
   +zNbRYJr+wkp5Qgp5aRWkeBqXyWlHCmlHOVNJCgUCkUr27dv59SpUwB63n5/kFLy1ltv6W/w/pz/
   0UcfUVZW5tf1PFFeXs7GjRuD7ge0cL5Nmzbp2wMHDuxwCcNut3P8+HHdy7/LkRK2b4cbb4Qf/hAe
   ekhLnPSDH3gUCaAtDRw8eFDfDsQvYvfu3Vy4cEHf7u0i4eTJk1RUVAAQGhrK448/bnj4qa8+Cg5g
   upTyoId9VwEHVfVIhULRmzGZTMTFxQUVflhXV6cn1OkMd4uDzWYL+Lrt19mNKgjlS4Gq0tJSoqKi
   AipiFTBSwtatmgWhuhr+8z9h8WLwMOlLKTl8+DATJ04MuKhVTU0N5eXlQdXw6G7c/w6cPHmSfv36
   XeTE2hM+Ch1drC/QbMBYFAqFwjCklBw8eFDPJdC3b1+/J2uz2cy2bdv0pQpfRYLZbGb16tX623cw
   4mTLli1tfCGCEQlr167V/QoiIiI6FQAVFRW630KXIyV89hlMmwY//amWRfHECXjwQY8iAbRn0dzc
   jNVqDeKyUs8RcSmQn5/Pxx9/rG+PHTvWkEiXjvBqURBCzAJaS5z9J1pGxdJ2h0UDCwGLlHJqF40x
   KJRFQaG4PJFS8vXXX3PllVcGnC7Y4XBw8uRJxo8f7/cEXVlZSXJyckDXdfe0D2Sd3RudRUm0Foxq
   Df/rFqTUijS98AI0N8Ozz2pJk7yM8cyZM9TU1HDNNdcEeDnJRx99xG233eZ3xEpP4HQ6OX36tJ5j
   w26343A4Ol1eMNKi0NHfvhvRBAJoDoPLPBzTApwEfmTEYBQKhSJY6uvrSUhIQAjB9ddf7/f5dXV1
   NDY2kpqaSmhoaKfpils5c+YMFouFSZMmAQQsEpqbm1mxYgWPPvpo0PkHzp8/z/79+1m0aBEAKSkp
   HR4vhKCsrIxx48YFndehU5xOWL9eS5QkpSYQ7rzTo0BwN7X379//ojLWnSGlxOFwEBYWhhCCa6+9
   NuClip7g5MmTDBs2jIiIiK7LSdEBvvooOIFpnnwUejvKoqBQXD6YzWY+/PBDli1bFnBugdzcXBoa
   Grjqqqv8Oq+yshKbzRZwTQX3ydAoT3uHw4HJZOowpK+goICIiAjDakF0itMJa9dqAiE8XBMIt98O
   Xiw2NpuN119/neXLlwe8hLN7925CQkJ8KmTVG9ixYwfp6ek+Z/j0RLdXj7yUUUJBobi8CMTZr6mp
   ya/iTa0cPXqU0aNHB+1lfuzYMaqrq5k9e3ZQ/YCWO2D8+PGkpKT4lDI5JyeHqKgohg4dGvS1O8Th
   gI8+0gRCbCw89xzccotHgdDS0oLT6STKlR/BbDb7VX+jfUppm82mWxN6Iw6HA7PZrOeuqKioICEh
   Qb//QOgWZ0YhRLoQItzte4cfIwajUCgU/nL27Fm2bPlnvrVAJoONGzdSXFzc+YHtaGho0CsOBsOI
   ESOYNm1a0P0AjB07jjNnatm9u4hdu05gMrWtwdfS0sLevXt1B83Ro0d3rUhwOODdd2H8ePjTn+CV
   V+DAAVi40KsVYefOneTl5enb/hbpak0p3XqPRuSw6ErOnj3bJgFYSkpKUCLBaDpyZtRDIl1LDx2+
   lqvwSIVC0RM0NzfT2NgYcElm0N7ofFmTdzqdlJeXk5qaGvC1QLN6rF27lrlz5wadyrm6upovvviC
   ++67D4fDwa5dJ4iKGkVUVDRWaxNWax433jhOtyw4nU6+/vprZsyY0bV+CHY7vPce/PrXMGCAZkGY
   M8erBaGkpITMzEwgMKtQVlYWmZmZhqTG7mpsNhtffvkl8+fP7zIB013OjP8CnHX7rmZbhULRKzhz
   5gyJiYkkJycTGRnpt+m/ubmZNWvW8MADDxAeHu7ThGm326msrGTXrl1eix/5ihCCadOmGZJhsV+/
   fsyZM0cPFbTZIujTR/Pmj4qKxmyO4NSpU8TFxTFs2DBCQkK44YYbgr6uV2w2rXrjiy/CkCHwt7/B
   rFlerQegCYUTJ04wbNgwHA5HQJUmw8LCui8ZVAA4HA5AS4oUFhZGcnKy7mDZ21E+CgqF4pIjOzub
   Pn36MGTIkID7qKio6DQKALQ38IqKSk6cuIDNFkF4eAtTpqT7bcGoqKjgxIkThvghbN26lYyMjIvK
   Utvtdo8WhaFD44iMjAyoYqLPtLTAW29pAmH4cM1JceZMr4fv37+fMWPGtMlNUVNj4siRIp+ec0FB
   ATk5OYaV0+5qPv74Y8aPH8+oUaO65XrKmdEPlFBQKL4b+OvQ1p7KykpKSkq48sor/Tpvz549nDlz
   npEjb/VqzveF5uZmiouLDamVYDKZiI+P93h9k8nEN9/kU1BwmszM0Vx1VUZQyzKd0twMK1fCSy/B
   6NGaQPAhuuDEiROkpaXpQsGbyHF/zq2WBtCqdzY2NtKvX7+uu7cgqKmpwWQy6ZELwWTnDIRuWXoQ
   Qqz0ox8ppfyBAeNRKBSKi5BS8o9//IM777wz4EkvIiLC5zV59zXyiRMnYjb3IyqqrTm/ubm5U6GQ
   lZVFWlqaHoEQqEior69n/fr1LF26FCFEh8+gb9++zJ49kUOHWrjyytFd5xRntcIbb8BvfwsTJsAH
   H2hZFb1QUFBAbm4uN910EwDjxo1rs9/bsknrc3Y6nbzxxhssW7aMmJgYoqKiepXDX3usVmubrJbd
   KRKMpiNnxgLa+iX0ARIBO1CNVv0xDKgDTFLKzC4daYAoi4JC8d3APVuhr9TU1BAZGelXZkYpJW+/
   /Ta33nor/fr18+lN1xu5ubkkJSUFXCuhVbBIKamuru4wH8KpU6eQUjJ27NiAruUzjY3w97/D738P
   V12lWRCuvtrjoe5ZJa1WKxaLxeuz8PScTaYsrrlmpL5E5KvTaU/Q0tLCW2+9xbJly3qF30G3hEdK
   KYdKKYdJKYcBS4EGYDEQLaVMRUvfvMTV/qARg1EoFIpWLBYL7733nu6gFkgCpezsbEpL22ee7xgh
   BLfffrv+1h4WFsaUKelYrXlUVp7Cas1jypR0j5NBS0sLhw4d0revuOKKgEXCV199xbFjx/QxdSQS
   QLMkdKkZ3mKBl1/W/A+++go+/RQ2bPAqEqSU/P3vf6ehoQGg0+JSnp5zYqKTyspK/ZjeJhIqKir0
   OhMRERHcddddvUIkGI2vmRkPAquklH/1sO//Ag9LKQNLvN3FKIuCQnFpIqXk/PnzfmcMDGQt2GQy
   sX//fhYsWOA1msFut3fqjW+32/nqq6+48cYbA5rU3Jc86urqiImJ8XovNpuNTZs2sXDhwq6dnBoa
   4K9/1fIf3HCDVs3Rlaa6PU1NTbS0tOh+By0tLX6lSm5oaODQoUNMnTrV76iHnmDLli2MHTs2KKfa
   rqInqkdOAM542XcaGG/EYBQKxeVNbW0tp0+fBrT/6PwVCU6nkxUrVmA2m/06LyEhgZEjR3YY8hgW
   FkZsbOxFk1dVVRUXLlzQj/ne974XkEhoampi5cqVehhdYmJih4InPDyc0aNHd10iofp6LYIhMxOy
   smDbNi2zoheRAJoFp/X3A/yupxAVFUVcXBwxMTG9UiQUFBSwf/9+ffumm27qlSLBaHy1KJwDtksp
   f+hh30pgtmuJotehLAoKxaVDeXk5hYWFAVcGBN/fYgsKCggNDQ36P/qTJ0/idDoZPz6w9yV3K0J1
   dXWH5vkTJ05gtVr9rkPhF7W18Oc/a1kUb7oJnnkGxozxeGhLSwsnT55k8uTJAV9u/fr1XHXVVb12
   wrVYLLqPS11dHQ0NDXpq6N5MT1gU/hv4FyHEZ0KIR4QQN7v+3AQ8DLxixGAUCsXlR0tLCy0tLQAM
   GDDAb5FQXV3NJ598oqfr9fUt1maz6W/v/lJUVKRfb+zYsbpIsNvtWCwW7Ha7T/0cOHCAAwcO6Nud
   +TOkpqYybFgXvZPV1GjZE0eMgLNn4euv4Z13vIoE0PxGKioq/Ep0JKXEYrHo2zNnzuy1E6/VamXV
   qlVtrDy9daxdic95FIQQPwCeA9wzdhQDv5RS+hNK2a0oi4JC0bv54osvSElJ0csz+4vD4aC0tJT0
   9M5LzhQWFpKenh6Uub41VHPBggVt8jr4mizI3YLQ0NBARESEV3HjcDj48MMPueuuu4IuPOWV6mp4
   9VV47TWtzPPPf66JBS/s27ePtLS0gC0AOTk5nD59mttuuy3QEXcpR44cYdiwYXoq6EDSSfcGeizh
   ktCeVhqQCpwHSnr7LKyEgkLRuwkk7LGkpASn0+mTOGhFSsm6deuYO3eu36mTnU4ndXV1XvMX+BpC
   abPZWLVqFQ8//LDPlo/i4mLS0tKMn6wqKzUHxb//He6+WxMIXqwV7pNlUVERffr08fkZOp1OcnJy
   GDNmjB7qCYEV7+oq3P8OHjt2jLS0tF6byMlXemLpAdCyKkkpi6WUB11/qhlYoVD4RWtBpJqaGiCw
   sMfm5mY9LK2za7WG5wkhuOuuuwKqr1BUVMRXX33V4Xhstog2SZlsNi1ZUOs4QHNAvPvuuzsUCdnZ
   2ezZs0ffHjJkiLGTank5PPUUXHEF1NVpjop//7tXkXD+/Hk++OADfTs9Pd2vZyiE4OzZs/rvJYTo
   VSLh+PHjfP755/r2xIkTL3mRYDT+LD0kALcA6UD7dFhSSvkrg8dmCMqioFD0PsrKyhg4cKBfIqG6
   upp+/fr5NclcuHCBbdu28eCD/qd6aWpqIjIyUh9jRybojiwKx48fp7a2llmzZnm9lnvfZrMZIQRx
   cXF+j7lDzp/XkiStWgUPPABPPw1eaj9YLBZiYmIQQuB0OmloaPBLHBw/fpzIyEhGjhxp0OCNw263
   k5+fr9dcaGlpQQhxSWdO9ES3WxSEENcBBcB7wEvA8x4+CoVC4RG73a4nDwIYNGiQ35aE7du3U1FR
   0elxUkrduW7gwIHcf//9/g3WxWeffUZhYaG+3VnopHuyoKamXD0p06hRo5jWQWpjh8PBm2++SWNj
   IwDx8fHGioTSUvjRj2DcOHA64fhxLaqhgwJRH330EVVVVYBm8fFFJLi/kPXt27dNsafehJSSEydO
   6H9HIiIivnMiwWh8DY/8BggFlgPZUsqWrh6YUSiLgkLR87S0tLBjxw7mzp3rV44B97VjX53KduzY
   QUJCQkAhhO7JmgLxnbDb7TQ1NfH+++9z//33d5g62v1+qqqq6N+/v7Em+aIirQ7DmjWwbJm23DBw
   oMdDzWYzFouFga79/t67yWRi3bp1LFu2rFctK7SyadMmJk+e7HdejkuZnvBRGAP8p5TycLAiQQix
   QghRLoQ45tb2OyHEKSHEt0KIta5ljtZ9PxdCnHbtn+/WvkAIkSOEyBNCPB3MmBQKhfFIKfW35IiI
   CG666Sa/RILFYuH111/X3/x8nYCuvfbagOL6zWYzK1as0N+M/RUJUkrCwsKIj4/nnnvu6VAkZGdn
   s337dn07KSnJuAm2oAAeewwmT4a4OMjJ0VIvexEJoPkhFBQU6Nu+3LvZbNbDBvv06cM999zTa0SC
   zWajrq5O354yZUqnKbAV3vH1X0IRYFRszpvATe3atgLjpJST0TI9/hxACDEWuBdNqNwM/FVohAB/
   cfUzDlgihBht0PgUCoUB5Ofnt3ES85fY2FgeeOCBTictp9PJRx99pIuSmJgYvwRJqzCIj48P+I34
   1KlTbNmyRd/2FB3hnmtg5MiRXH/99X5fp0Py8+GHP9QKNSUlQV6eZlFwFVRyR/Op2KWPqbPlEU9s
   3rxZz0gphAjISbSrOHnyJEePHtW3Bw4c6HeWSMU/8VUo/BL4mfubfqBIKfcApnZt26SUrf+K9vPP
   XA23A+9LKe1SygI0EXGN63NaSlkopbQB7wN3BDs2hUJhHJmZmdxxh3//LAsLC9ukyPVlrT4kJISp
   U6cGVHJ47969ba7nT64C9yXNzMxMbrzxRq/HOp1OXn/9dT3RkKElkk+f1pYWrrkGBg3Stn/zG00s
   eCE0NJTIyEi/Ek6ZTCby8/P17e9///u9JvlQY2Mj27Zt07cnTZrEzJkze3BE3y18TaZ9KzAAOCeE
   2AfUtNsvpZQPGzSmfwHWuL4PBva57St1tQm0ZE+tlKCJB4VC0YMcPHiQ6OhoJkyYgBDC75oHffv2
   9cns3TpptfohDB06NJDhMnny5IDeNKWUrFq1ikWLFtGnTx+vAqO1LHJISAgPPPCAX+WuOyUnRxME
   mzfDk0/CmTPgShLkiQMHDhATE6P/Nv5aEBobG6mpqSEzMxPo+TwIdrud0NBQhBBERUXRr1+/SzY5
   Um/HV6FwPSCBejRTf3sM8RYUQjwD2KSUrULB0y8u8WwJ8TqG559/Xv8+a9asDsOUFApF4AwfPtzv
   N+Vz586RnJxMXFwcCQkJPpmww8PDA5oQnE4n77zzDvfccw8xMTHExMT4dX7rRCSEYNGiRR169mdn
   Z1NSUsLNN98M+GYd8YmTJ+HXv9aKNP34x/CXv4CXcbQKFYARI0b49ds4HA42bNjArbfeSnh4OIMH
   D+41FgSA999/n1mzZpGWlkZISAhTpkzp6SH1KDt37mTnzp1d0rdfmRkNu6gQGcBGKeVEt7aHgUeB
   70kpm11tP0OzVvzWtf05WhppATwvpVzg6bh211JRDwpFF3Ly5ElGjBgR8Brwvn37GDJkCGkdhOuB
   lhMhOjo66LC78vJyUlJS/BYa+fn5ZGdnd7ic4h41YbPZEEIYVwUxO1sTCDt3wr/9Gzz+OLilkG5P
   bW0tH374IcuXL/f5XqWUSCl1q05eXh6ZmZm9opJjeXk5VquVjIwMILBy4pcTPZaZ0UAEbtYCIcQC
   4P8Dbm8VCS42AIuFEBFCiGHACOAg8A0wQgiRIYSIABa7jlUoFN1MeXl5myI/vuBeBnr69OmdigTQ
   LA+tsf3+UFRUxI4dO/TtAQMG+DVxtpKens7cuXM7PPaNN96gvr4e0KwehkywR49qKZbnzYOpU7WC
   TU8/7VEk1NbW6gWp+vTpw8MPP+yXINq5cyeHDx/Wt0eNGtUrRAJoSx+tWTYBJRK6k1YF2dkHiAV+
   BPwD2AGMdLUvBkb70c97QBnQjBZNsQzNSbEQyHJ9/up2/M+BM8ApYL5b+wIg13Xuzzq4nlQoFP5h
   s9lkQ0ODtNlsHvc3NjYG1fdrr70mm5qaOj3WZDIFfJ1WLBaLLC0tDejcd955R164cMHrfqfTKa1W
   q77d3Nwc0HU8cuiQlHfcIWVqqpSvvCKlxdLpKevWrZPFxcU+X8LpdMqKigp9u7GxUTocjoCGazT1
   9V14B3gAACAASURBVPVyxYoV0ul09vRQLklcc5/Pc3xHH18TLg0BdqJFI+QA44GpUsosIcTfgFAp
   5Q+D0Ctdhlp6UCj8o7MqiFarlTfffJPly5f79bYp3RzNfEno43A4WLFiBUuXLiU6Otqve/jqq68Y
   P358QDn73cdZX19PfHy817fy48ePU1hYyMKFC/2+jlcOHoQXXoBvv9UsBz/8IXi5/7q6OiorKxnR
   QbXHjjCZTGzevJklS5b0CifA0tJSBgwYQFhYGFJKTCaTqrsQID2x9PAymgVgJHAVbZ0MdwEqDkWh
   +A5gt9s5cqSIqKhRJCePISpqFFlZRbo5G7TQvkcffdQvkVBWVsbatWv17Y5EQmvIXmhoKMuXL/db
   JEDgcfMlJSV8+OGH+nZCQsJFE2hrvgaAsWPH6s6KQbNvH9x8M9xzD9xyixbF8OSTXkUCaKKtsrLS
   r8t8++23+j307duX+++/v1eIBICsrCxMJi16XgihREIvwVehMA94TkpZxMXRBa0hiwqF4hLHWxXE
   goICtm7dqh/nb9jjwIEDmTNnTqfHnTt3jnXr1unbHU1gdrsdi8Wi/3ngwAF936hRo3yOMnC3OA4a
   NKhD64CUknfffZfa2lpAEzyBVL9sw549mv/BkiVw551aHoTHHwcPEQpSSj777DO9KuWAAQOYPn26
   X5draWnRz+9pcnJyyMrK0rdvu+02kpOTe3BE3w2MtqL7+koQAZi97EsEbMYMR6FQ9CSRkZGEh7dg
   tTbpVRDDw1sYNCij01BCu91Oc3MzkZGRhIWFcerUKSIjI8nMzCQkJMRjtsL2ZGRkkJqa2ulx7ZdH
   JkxIxeFwBBRHv3btWqZOnUpGRgYhISEXCQyn00ljYyNxcXEIIfjBD34QvDiQEnbtgl/+EgoL4Re/
   gIceAi9WkNb7EkKQmZnp1z3m5+dTUlKiJyC65pqeTTljsVj0fBLJycl+JX1SdM6RI0f0Eu5G4auP
   wj7guJRyuRAiFE0YXO3yUXgNGCWl7Px1oQdQPgoKhX+YTCaysoqorW0kIsLO9deP73SS9+TXUF9f
   T2RkpF5oyBtffvklQ4YM8bkkcWtJZ6s1lpiYeCIiEvSSzr4uh7gLCrPZrIsAT5w6dYr8/Hxj/BCk
   hO3bNR+E8+fhmWe0ks8dePBnZWXR0NDgV6ZBu92uP4uGhgYaGxtJ8ZDKubupq6vjgw8+8CtkU9Ex
   NTU15Obm6pYlq9VKaGgoERER3e6j8HvgB0KI1/mnP8JYIcQvgR+49isUiu8Affv25cYbx5GeHsbw
   4f07FQnufg0hIdFERo4gK6uIwYMHdyoSACZMmEB6errP42tdHnE6rTQ3N+jLI76a0ysrK3n77bf1
   bU/OiiaTSTffjh49mltuucXn8XlEStiyBa6/Hp54Ah59FE6dgkce8SgSmpqa9O9XXHEF1157rc+X
   stvtvPbaa/rziIuL61GRsG/fPt0nIjExUYmEIHE6nRQVFenbUVFRbaxgUVFRhoeO+iQUpJQfA48D
   3wdaE2qvBn4CPCGlDLzyi0Kh6DVYrVa9CuLUqVMZO3Zsp+e0TtyRkVGUlBxFypYOJ26bzca2bdt0
   k3NycrJPNRaklJw5c4aICM1q0a/fcJKTh+vLIx310RrmBVqlxkWLFnV47MaNG9s41QU8sUkJmzbB
   9Onw059qIuHECXjwQfBi/bBaraxcuVJ/PrGxsZ0+n7q6Oj03RVhYGI899phfdSuMxr0AVmRkZBtn
   WCUS/Mf976+Ukt27d+vPtDUtd1fiV2ZGIUQsMANIBqqBvVJKb74LvQK19KBQ+M7q1auZP3++T5aA
   ViwWC/v3nyEqapTu19DRUoCUkqysLCZNmuRX5ERrSuGbb76ZpqYmsrK8h3C2Z+PGjYwcOZLRoz0X
   mbXb7dTX1+te9oH4OrRBSti4UVtiaG6GZ5/VkiZ58W2oqKggJiZGfzP0JXzUnT179tC/f3/GjBkT
   +JgN4uDBgzQ2NqpU+QbyzjvvMHv2bL9SaBsZHtkjKZy7EyUUFArfcV/b9oXWNed77rmHI0eKvU7c
   ZrOZ6upqv4s3tbS0YDab6d+/v8exujtPtsd9sm1oaCA2Ntbr5H/27Fny8vKCD3V0OuGTTzSBIKUm
   EO6806tAaGXXrl0MGTJEL7jUGWazmZycHKZOnRrceA2gubmZoqIi3cekubmZsLAwvyNjFP8kOzub
   6OhoPT9GU1OT32HCPSIUhBB9gJ8C09HCIUuBvcB/SylrjRhMV6CEgkLhnebmZjZu3Mgdd9wR8Lpm
   a879jibu0tJSCgsLmTFjhl99///snXd4VNed9z93+qgXVKhCYEBYNGNjMBiLalzjbmM7bmuT5u15
   d5NsNnac4s3meZI3b3aTuCU27jau2MaAMYgOookiBBIIoYJQm1GbmTv1vH9czfWMNDMaCdHv53n8
   mFvm3DMjac7v/Mr3d/ToUerr61mwYEG/Xtfe3s7777/PU089FdU4aGhoICcn58wrGEAxED76CH75
   SyWk8Mwz8K1vQZRnt7e3c/z48QE3MnK73ezdu7ffpZFnA4fDwYYNG7j11lu1sMIAcTgctLe3M2zY
   MEDR8zCbzWdUKnrODQVJkqai5CakAjuARpS207OANmChEOLgYExosNEMBQ2N2FRVVZGfnx/3l/zh
   w4dpbW1l7ty5Me9ra2sjMTGx3wZIX2qI0QjGcUO9CLG0FD788EOKiooYMmRIv54Tht8PH3ygGAgJ
   CfDss4pYUh9zdzqdHDx4sF9Jih9//DFz5sy5IKoXPvjgAxYuXBhXyatGZEI7e9bU1FBdXd2vypa+
   OB+GwkYUw2CJEOJkyPnRwGrgtBBi3mBMaLDRDAUNjXBcLhcNDQ1xu7l70tXVhc/nIy0tLeZ9n332
   GZMmTSI/P79f47/99tssWLCgX3kSAF999RVZWVlMmzYt4nW3243NZotLp6FPfD547z2lm2NammIg
   LFkS00BYuXIl119/fdxqg0IIPJ5vkjQbGxsZMmTIeXHpy7KM1+slubsRVUNDA1lZWRdMw6iLjc7O
   TpYvX87TTz991rww58NQcAKPCSFWRLj2APCqEKJ/jd3PEZqhoKERTnNzMwcPHuyXO3/fvn2MHTuW
   lJSUmPeFtv7tT0JgMGTR39eF5iE4nU4sFkvUUEJdXR3l5eUsXrw4rrEj4vPB228rBkJOjmIgLFwY
   1UAInV9tbS25ublxe1gOHDhAXV3dmZdmDgI7duzAYDBwzTXXnO+pXJQIIVi5ciVLlizB0q246fF4
   BtyaPR7OR6+HVpReD5GQu69raGhcoAQCAbxeRUA1Kyur3zF/IKzELRJut5uXXnpJvS/exb6lpYW3
   3npLLf+K93VOp5MXX3xRLSNMSEjoZSScOHFCnc+IESMGbiR4vfC3v0FBAbz6Krz0EmzaBIsWRTUS
   ysrKWLVqlXo8cuTImEZCIBCgqqpKPS4sLBy8PhL9pL29neLiYvV41qxZmpHQT2pqatSW45IkceWV
   V4b9fp5NI2Gwidej8B/AzcBiIYQcct4KrAW+EEL85qzN8gzQPAoaGrB582YMBkO/kt9Onz5Nbm5u
   n9UFocS7S+qZTxBvtUXP14XKAUdi9erVXHPNNQPPQ/B4YPlyeP55uOIK+NnPIEYc2WazqaGFoGEW
   rwfB7/fz0Ucfcccdd5yXRST0Z+f1eikvL2fKlCnnfB4XK0IIvF6v+hlu376dESNGMHLkyPMyn/MR
   evgV8ChgAVbxTTLjLYALRXwpOJAQQjw7GJMbDDRDQUNDWYj1en3cu3WPx8O7777LkiU3cehQQ9Sy
   x4qKCk6ePNnvnXpf+QTR2Lx5M0ajkVmzZkW87nA4aGpq6ndeRC/cbsWD8JvfKF6EZ5+FPio2/H4/
   f/vb33jkkUdU93JPehpdhw4dIjU19bwtJqG89NJL3HvvvVrHxgGya9cu2tvbWbRo0fmeCnB+DIVA
   nzd9gxBCXDAFtJqhoHG5smbNGqZPn96vEqvQTOxgT4VYQkqyLON2u0lNTe1z7NDcA5fLhcViictw
   CY3zu1wuTCZT1IS+5uZmDh06xPz58+N6v72QZXjlFfjv/4YpU5QyxxiVCadOncJgMKiVCLHyKyL1
   w7Db7SQmJpKTkzOw+Z4BtbW16HQ6VcSnvxoalzstLS3s2rVLDQ/5/X50Ot0FUyJ6znMUhBC6fvx3
   wRgJGhqXMwUFBXEt4EGqqqr45JNP1ONoLac3btxIS0tL9zlLXM/weDy88soreDweAKxWa1xfqF6v
   l7/85S+qG99qtfYyEg4fPowsKxHRrKysgRkJTif84Q8wdiysXQsffwxffBHTSABobW2lvb1dPY72
   noL9MPz+TKqr92OxjGfv3hpGjRp1To2E0E2T0+kM6ymhGQmx8fv97N+/Xz1OSUkJkzjvj8fuYmMQ
   lEY0NDQuBIQQHD16VD3Oy8uLGuv2+Xw4HI6wBMXRo0eHdUgMbTkNqD0Vhg8fHtW1Hg2TycQ999wT
   d+w9mKBoNBp54oknYsb5W1tb6erq6td8VBwO+N3vFANh0yb4/HNYuRKiJO51dXXx9ddfq8eTJ0/u
   s+ul0+nsLi80kZExjIkTF/e7kdVg0NzczNtvv60eT5gwQVX+04iM2+1Wfxd1Oh319fWq0WoymcjL
   yzuf0ztnxGUoSJL0/RjXzJIk/WnwpqShoTEQfD4fR44cUXft0bDZ7GzcWMbmzTW8/fZK9u3bByhf
   hKEGgMFgYPr0UTgchzl8+EtkuYLp00dRUFAQU8goyNGjR9m0aZN6HG/se9euXWGvS0gIr7zu6Ojg
   8OHD6vHcuXP7n6zY2amEF8aMge3bYfVqRVnxqqtivsxqtZKZmUl/wpkrVqygo6MDo9GD2+0mKSkz
   rkZWZ4oQgpMnT6pzzczMHJxW2ZcRK1asoLGxEVC8Rbfccsugd2a8GIg3R8EPrASeFELYQs5PAt4F
   8oQQyWdtlmeAlqOgcanTn3rsnnkHTU0ngCYWL746quvZ5XKxfv16Fi1a1K+FraurC6/XG5d6X2hu
   RF+9Atrb2zl06BBz5syJey4qHR3wv/+rhBkWLoSf/hQmTVIvR6rw+PzzzyksLIw7QbKtrQ2Hw6HG
   /oM5Fna7vV+NrM4UIQQff/wxCxcu7FcI6nJm9+7dGI1Gpk6dCgxCc7DzyPnQUbgZRa55vyRJ87on
   8Y9ACYq+wtWDMRkNDY3+UV9fz4oVvXTQoiLLMqdOnVINi+zsfCQppZcLPKhiCMou+tZbb43LSPjy
   yy/VmH1SUlLcRsKLL76I0+kElJBHTyNh3759ahvl1NTU/hsJbW1Ko6axY+HwYdi4Ed55J8xICPW0
   fP11qdpmevbs2YwaNSruR7W2ttLQ0KAeBxMx09PTKSoqZO7cURQVFZ4VI6G0tJSysjJAWSjuvvtu
   zUiIgc1mo6KiQj0eM2ZMmGLpxWokDDbxJjOuBaYCZcA6SZL2AL8H/gLMEkJUxHq9hobGwImUTxBk
   2LBh3HfffXGPZTab8XjsdHUpfdyiucArKyvD3PvxMm7cuLjzF4KxX71ezxNPPNErxBBKsEa939jt
   SmnjFVdAVRVs3Qpvvgk92jEHkw0tlvFIkpmqqkr27q3B5/ORkZERUzbZ6/WyatUqAgGlOGzs2LFR
   xYkMBgOJiYmDljgohAjLzxg2bJjaWEijN0II1eAE5eceepyRkaHKVGt8Q7/aTEuStBD4HDADe1AE
   mC7YzpGghR40Lm4ildSdOnUKn8+nukf7QghBe3u72pshmgvc7XZjMpn6tYuy2+0cOXKk310MDxw4
   QH19fVTlwba2No4dOzZwNcDWVvj97+GFF5Q2z//xH4o3IQoVFRUcP+4jO/vKblGnAK2tFcydOyqq
   oFPQLS2E4ODBgxQWFp7zPgynT59m3bp1fPvb3z6nz71Y6erq4q233uI73/nOJe8tOOehB0mS9JIk
   PQ+sAdYDDwEjgVJJkq4fjIloaGiEE7rLzcqaqJbUZWZmqvHveGhpaeHzzz9Xk9qiucBXrFgR5jKP
   B6vVitVqjeveoAcBYOLEiTFFmgasTNjcDD/+MYwfDy0tsHs3/PWvMY0EIQSlpaUI0YEsu5AkCY/H
   EzPZcN26dRw4cABQvpCnTJlyTowEIQQbN25UvSu5ubk8/PDDZ/25FytCCN58803Va5CUlHRZGAmD
   Tbw5CtuAfwX+jxDiViHEu8A04CiwQZKk5+J9oCRJf5UkqVGSpAMh59IlSVorSdJRSZLWSJKUGnLt
   j5IkVUqSVCpJ0rSQ849JklTR/ZpH432+hsbFQqiOQWvrSQwGPV6vieTk5Liy/IOGQVZWFg8//HDY
   l2PQBR66uC1dujQut/WRI0fCdBTiUVcUQvDXv/5VzV8wGo293O87d+6ktVVpG5OQkNA/b0JjI/zb
   v8GECUrC4r598OKLECUBsa6uTu2rIEkS999/P9ddNwFZrqC5uVyt8AjOMeiVCTJr1iwmheQ3nE2E
   EGpYQ5IkEhISwsJQ2qIXTnl5Oc3NzYDy2dx4441hXiHt8+o/8RoKKSi5CH8InhBCnBZCLAF+DPx7
   P575KrCkx7kfA+uEEBNQPBY/AZAk6WZgrBBiHPBd4IXu8+nAM8AMYCbwbKhxoaFxKRCqY9DYWIHd
   3hB3Sd2uXbvYsmWLehzpy7Gzs5PXXntNXYTijZt7vd648wWCXgRJknj00UdjJtalpKT0f1fe0AD/
   +q9KzoHbDQcOwJ//DH0kH/r9/l45H7GSDZuamvjiiy/U46SkpHMWZiguLmb37t3q8YwZM+L24lwO
   eL1etfkSKD/bUO9VdnZ21I6iGvERb3lkghDCGeP6NCFEadwPlaQ84DMhxJTu4yNAkRCiUZKkXGCD
   EGKiJEkvdP/7ve77yoF5wPzu+7/fff4vQHHwvh7P0nIUNC5KvF4vXV1dAyqpc7lc6PX6Pl34oU2M
   ouH3+zl69GiYCl08VFRUcPjwYe68886I1+12OwcOHKCoqKhf4wJQX6/oILz5Jjz6KPz7v0O3NyRS
   iaPb7ebzzz/nrrvuinvROHjwIAUFBQNqm30mOJ1OGhoaGNsdLpFlGZPJpC12Udi3bx9dXV3MnTv3
   fE/lguJ8SDirRoIkSUmSJOVJkmQMuR63kRCFbCFEY/dYp4Hs7vPDgdqQ++q6z/U8X999TkPjksDr
   9fLSSy+RkJAQd0ndhg0bVPe41WqNaCRUVFSoAksQnwhSIBDg+PHjfbaZhvBW1GPHjuWWW26Jem9i
   YiKZmZl9jhlGTQ08/TRMngwmk1Lq+Ic/qEZCaIljcfEhtcTTbDbHnfwZxGaz4XA41ONz5bKWZZkT
   J06oxxaLRTMSQmhqauKdd95Rj6+66irNSDjLxF2jI0nSbcAvUMokBXAtsFeSpFeA9UKIt2O9foD0
   /MuUup8d6S82qtvg5z//ufrvefPmMW/evEGYmobG2cNoNPLUU0+pYYZ4wgJDhw7tUzUuMzMzrrr6
   rq4uZFlmyJAhGI1Gbr/99rjm/eabb3LLLbeQnZ2NXq/v5Z7fsmULY8aMYdiwYZhMpvjj/NXV8F//
   BR98AMuWwZEjkJ0ddkto8mdampUDB76gsfEkt902j8TExD7lio8dO4bdbmfGjBkAA/N0DAAhBG+/
   /TZ33XUXCQkJZGRkXDAdCC8E/H4/mzdvpqioCEmSyMzMjFotczlTXFxMcXHxWRk7LkNBkqQ7gQ+B
   r4EfAb8NuXwCeAw4E0OhUZKknJDQQ1P3+TqU6oogI4BT3efn9Ti/IdrgoYaChsaFSmtrK/v27VMX
   ib5yEYI7/WCvgYKCgoj37dmzh8LCQiwWS9w7+BMnTuBwOLBarWEu/EiEdh1cunRpTB2FESNGkJKS
   EtccAEX74PnnlSZN3/seHD0KURI5ZVmmq8tHWpoSvx827Br27i1hw4ZjJCcbIoZtQjtTZmZmxiVN
   PRgEPRWJiYlIksSCBQv63T/jUqalpYXU1FSMRiN6vR6z2YzP51OPg6W+Gt/QcxP83HNx1xj0Sbz+
   rGeBV4UQNwJ/6HHtENDf9F+JcK/ASuDx7n8/Dnwacv5RAEmSZgFt3SGKNcBiSZJSuxMbF3ef09C4
   aElOTu5XkxmPx0NZWVlY4la0++JpPtTS0qJWSgwfPgKHw8rmzTVs3FimqhT2pKamhg8++EA97rnY
   2e32sCTA0aNHx7cYV1bCE0/AtdfC0KHK8a9/HdVIACU5s7q6BFl2dfe9aCQhIZuhQ6eqpaWhoRGP
   x8Of//xn9Vx6ejq5ubkRx44lejUQSkpKqK6uVo+HDh162YcXQnPJtm7dqlbWAFx33XWXZY+FC4V4
   kxll4HYhxFeSJOkBL3CNEGKvJEk3AGuFEHGZw5IkvY3iDcgEGlGMkE+AFSjegxrgvqCQkyRJ/wvc
   BDiAJ4QQe7vPPw78FCXk8CshxOtRnqclM2pcsHR0dODxeOJuahRU4tPpdFF3+oFAgFOnTjFixIi4
   5yGE4K233uK2224jKSkprB+ELLuQ5QqKigoxGAz4fD61pa4QAo8neiWG3++nurqaMWPGxBfjP3JE
   MQhWr4Z/+Af4x3+EGLvHiooKRo8ereZj2Gw29u2rpbPTT3l5HbNmzSUlRfEiNDeXc9VVGaSkpKhV
   A06nM6YipDJmb9Gr/sovNzc3U1FRMbD+FJcBW7duRafT9Vu4SyM6g5nMGG+OQgcQ7ZtsNNAc7wOF
   EA9FuRQxKCeE+Pso518DXov3uRoaFyInT57E5/PFbSjs2rWLffuOkJc3I+qi1dnZyY4dO7jnnnti
   Ls5CCJxOp+r+DmotOBwOvF6T6sK3WKx0dprUZk0rVqzg+uuvZ+TIkUiS1MtI2LBhA6NHjyY/Px+9
   Xq9m78ekrAx+9Sv4+mv4p39SGjfFkUtRW1tLZmamGlLJyMigqCgFh8NBQoLAZFL2L0Gp6tLSUvLz
   8xk/fjzQuzNlT3rmPciyi717KygqSu4zbyS0WdeAEjcvYRobG6mrq+Pqq5U2QVddddXARbY0zjrx
   +rq+An4iSVKoaS8kSTIDfw98Oegz09C4RPH5fKqbdfLkyVzVR2vj0Nd1dJi48sp7wpQafT4ffr8f
   WZYBpWnSvffe2+cOvq6ujs8//1w9Dt4fqt8AyiKr17tUg+Duu+9m5MiRvQfsZuLEiQwdOjSu98TB
   g3D//bBgAUybBsePKx0doxgJ9fX1lJZ+U2S1cOHCXguwwWAgNTWVGTPG0N6+n/LyNaqA0pIlS1Qj
   IR5CRa9AMZq8XlOfoZxAIBDW6CohISFqDsnlgN/vp66uTj02m81hIkgJCQmD1v9CY/CJ11D4KZCL
   osT4Coq7/8dAKUoi4c/PxuQ0NC5FPv74Y06ePBn3/WvWrOHYsWO43W58PjMJCcoXbOiitWPHDvbu
   3dvnWD6fTxVYGjlyZMSGUgaDkvgXVClsbd3FyZN71S/ynh6EtrY23nvvPdX4yc3N7Tsxr7QU7r4b
   Fi+GGTMUA+FHP4I+GvJYrda4Ew7T09O57roJDB9uGnC3xkhGUzTRq+PHj9PUpORh63Q6vv/97/fp
   sbiUCQ35er1eNm3apJ5LS0u7rA2ni424m0JJkjQCeA5FVTEbaAVWA88IIWpjvfZ8ouUoaFxoyLLc
   rwz31tZWtaQxNHfA5XLidldSVFSITqdDkqQ+vQiffPIJV155ZZ+76qC6ndfrxWQyEQgEouYhCCFo
   aGiIr2vhnj1Ku+dduxTJ5e9+F2Ispn6/nzfffJMHHngg7s/sk08+Yf78+YPWXjlaEy0IF2E6dOgQ
   qampMb0tlxMvvfQS999/v1ahcJ4YzByFfnWPvBjRDAWN843P52PNmjUsXrwYk8kUUTkwFL/fz5Yt
   W7j++ut76RAEFy2Px8iRI+u5885bye/uZxBt3NDFLDRuHouPP/6YwsLCqAbFunXrGD58OBN7tGuO
   SkmJYiCUliqeg6eegigyxEIItRQOlHh2dnZ2TCMotESzpqYmLk2J/hDps62traWkpIR77rln0J5z
   MbNnzx6ys7NVQ8nlcmlS0+cRzVDoB5qhoHE+6LmwlJWVUVBQQHt7R58Z9EIISkpKoiZ4Bcd2u92q
   smK0zHyn08kbb7zBsmXL+iy/83q96uIayaAINTjsdjtJSUl9L8bbt8NzzykKij/5iVLy2IdnYOPG
   jZjNZmbNmhV77G727t1La2trzG6Ug4EQgpMnTzJ69GhA+Tm4XC6S+wiXXKp0dHTgcrnIyckBoLq6
   mpSUlLjUPjXOPpqh0A80Q0HjXGOz2dmz5wRtbQ4yMlLURdvn80UtOwSlhC5WEqDNZmPbtm3ceuut
   YbvrWOMaDAY6Ozv7XMza2tp4//33WbZsWcSde0dHBx9++CGPP/54fGWOW7YoBkJlpWIgPP44RAld
   BAIBGhsb1fceWn4ZCSEEp06dUltte71e9Hr9WdchCAQCrFixgjvuuOOyFUfy+/2ql+vo0aO0t7dz
   7bXXnudZaUTinPd60NDQiI9gOR3k0NlpC6tMiJVB39LSEtYhMBKpqakUFhb2WkB7jtvUVE5NTZWa
   mR/NSAjtspeWltbLCBBCqMlnKSkp3HnnnX0bCcXFSgXDo4/C0qWKofDd70Y1EkDRMtiwYYP6LIPB
   EPM5Xq+XDRs2qHM3Go39NhLiFVDavn07x48fB5QExf7kSlxqNDU18dprr6nHEyZM0IyEywTNo6Ch
   MUgIIejs7GTbtnqysr6J3Tc3l3PddUqi386dx0hMvLI7GdGBy3WU+fOnRC0NO378OEajkVEx2ib3
   9Ci0tzfj8VSzaNFVMUvOvvzyS4YPH86UKVMiXl+7di25ublRr4e8cVi/XslBOHUK/vM/4aGHIEZY
   4vDhwwwbNizuRLeysjKysrLI7tHfYSDEElASQuBwONTKioaGBpKSki7L8ILf72flypV861vfQq/X
   98od0biwOR+CSxoaGn2wd+9e7HY7RmMWsuxSwwAuVws7dzoJBBKQZSdO524sliGcOrWfMWNyn2BM
   zQAAIABJREFUz7h+XK/X09R0mLQ0N52dyRiNHmbOvCLiuKG5B4sWLep1T2jvg9mzZ8cu7xMC1q5V
   DITWVsVAWLoU4ng/sizHJSsdRKfTMRgGf18CSlVVVRw4cIC77roLIH49iEuEqqoqhg4ditVqRa/X
   M3HiRPVzlyRJMxIuUzSPgobGIOH1epEkic7OTrWcTqdz4nQ6ycq6RjUcHI7DzJx5BUajEavV2svN
   fvz48fglj7uprq5m2LBheL3eqNUUTqeTV199le9///sRXfVdXV28+eabfOc734ntyhcCvvxSMRA6
   OuBnP1NEk3pUaIRy+vRpysvLmT9/flzv5/Tp0+zcuZM77rgjrvvjxeFwsHlzjerxCQT87N//EQ88
   cBPJyclhi+LlQE8vQXFxMRMnTlQTFDUuXs5JMqMkSc/0YxwhhPjlYExosNEMBY2zydatW7niiit6
   fbEGcxJ8Ph/bt58iK0vZmZWWfkJW1kgWLZoYpkwXRAjBypUrWbBgQUx39+nTp6mqqmL27NkRrwef
   bzQaw3aCoZUNoHgQhBBqglrM3gdCwGefKQaC2w3PPAP33ANRjIrQKglZlmloaFBLOSMRrBIJzr+1
   tXXQFyyfz0dx8SEslnFYrYnIsovjx7/k3ntvuizFkTZu3IjRaIz6e6Rx8XKuDIVAj1OC8I6PoecR
   QkTfTpxHNENB42xy7NgxcnJyoi7qPfMHTp+uxGTqYv78yequP5jb0J/2yw6Hg/r6+og6B6Ex+GB4
   I7T9bChr164lKysrqoy0z+fD7XJhWbMG/a9/rRgLzzwDd94Z1UAIvqdXX32Vu+66Ky5FRCEEr7zy
   Cvfcc89ZL6/78MMP8XqtpKXlD7jJ08VKfX09R48eZcGCBYCShxAU69K4tDjn5ZGSJF2J0vL5JeBd
   lK6POcCDwDLgNiHEkcGY0GCjGQoag01NTY3aEKkvZFlmz549yHJSVO2E06dPs379eh56KFq/NIU9
   e/YwduzYmAmAPp+P9ev3qwmTDkcHHk8V8+ZNUg2T0DwFt9uNyWSK+F5sLa3U/8+r5L3xChh06J79
   KUkPPQQ97g16L4KdJYOekr7KMjs6OnC73WRlZanjnA29/46ODmw2m6p/IMsyer1e7Xp5KfcY8Hg8
   VFRUMGnSJEDxGNntdrW0VOPS5XwkM/4v8IoQ4rch52qA/5YkSQf8CVg4GBPS0LiQ8fv9bN++nSFD
   hsTtqvZ4PNxww5VhC1Ow9FCn05Gbm8uDDz7Y5zhBKeVYdHZ2sm/fRubPVyoVEhNTcDrN6kLeM08h
   oiyz34//3Xcx/fQZRiekc/qffkfT1fOR3ZUU+f1hC2uo96Kl5ShDhli46aabgOhlmUFqa2uRZVk1
   FM7Wgt3V1UV9fb1qKATLGy/VxDyXy4XFYkGSJHQ6HdXV1WpZbUJCwmUZYtE4M+L1KDiAO4QQ6yJc
   WwR8KoToHXC9ANA8ChqDQajQTF94vV6cTmfMXgPr168nLS2N6dOnR73H7XZz9OjRPssThRBqTwbF
   o3CAxMSJavKk01nOnDkF6gLRM09BxeeD996DX/0Kf0oK+25/EummZaoHobm5nLlzR6keA7fbzcqV
   xQwbdkOv3hPB66E7dlmW2b59O/PmzTurrm6fz8fbb7/Ngw8+eMkaA7F4+eWXuffeey/ZcEpfEuga
   CudDcKkdiKaPemP3dQ2NS5KWlhaWL19OV1dXnwI9oNT879+/P+Y91113HdOmTetzrMbGxj7LAnfs
   2MGOHTsAZVc+Y0a+2vlRlitwu09x5Mg3kcFei6fPB8uXw5VXwgsvwP/8D2LLFtqvm4nsVlpXR+qa
   KMsyjY2n1PGs1gS8XhMNDafZuLGMzZtr2LixDLvdDiidGJOSkvr0igyE9vZ2tc22wWBg8eLFl80i
   smXLFioqKtTjJ5988pI1Emw2e8TfLY2zS7wehV8BPwFeAFbwTY7C/cB3gOeFEP2pkjhnaB4FjTPF
   ZrNTUlIJJEdNfgvVHwjN9g+9vmLFCm6//faYrt+WlhYAhgwZEnNOPSsEekoet7W1YTQaMZvNSJIU
   2Rvi9cIbb8Cvfw2jRilJivPmqR6ESF0TT506RWpqKiNGjIgoHe1wHAZQcyQOHlxFQoKOu+/urdkw
   mHz55ZdMmDCBMWPGnLVnXCg0NTXhcDjUCpLm5mYSExMv+ZBCX1LlGuGcD4/CM8DzwKPA18AhYD3w
   SPf5nw/GZDQ0LhROnDjBgQMHVIGe5OTJZGVNDJNkDuWjjz6iuroaiFyDr9PpmD17dp/d9Orq6mhq
   aop5j9/v55VXXsHlcgG9JY9lWeadd97BYrFgMBh6GwkeD7z0EowfD2+/Da++Chs2wPz5YYmK6enp
   FBUVMnfuKIqKCklPT8disageBIPBwPTpo8K8FxMnZuN261Q56bFjryMpaUS/xJXiob6+npKSEvX4
   5ptvvmSNhGBVTBBZlunq6lKPs7KyLnkjAXpLlYdKoGucXfoluCRJUhowBcgFGoADQogLOuygeRQ0
   BkJzczMul4vMzMwwgR7lWnisHpSEucTExLAFu7W1laqqKmbMmBHzWadOnWLYsGEx7xFCqMmQ0LtC
   wO124/f71QUj1MMRchP87W/wm9/AxImKUNKcOb2eFRoDbm9vp7i4OGYr5dD7a2tr+fzzr5g+/ZFB
   3/WFelHa29tpbW29ZI2DUJqbm1m1ahWPPfbY+Z7KeUXzKPSP89YUSgjRJoTYJIR4Xwix+UI3EjQ0
   +kNzc7O6O8nKymLUqFGYzWaMRg+yrOzeg7F6SZL47LPPVM9CUlJSL0+CxWKJ2CY6lEAgwIYNG3A6
   nTHvKy0tZePGjepxzy/GHTt2hOUhhBkJsgz/8z8wdix8/jm8/z6sXh3RSLDZ7GzYcIBNm06ycWMZ
   QEwxHiEElZWVWK1WDAYDo0eP5qGH7gnzMkyfPuqMv8i9Xi8vvvii+nmnpqZeskZCIBDgtddew+Px
   AMrv4qOPPnqeZ3X+ieTBGozfLY2+idujIEnSUOCHQBGQAbQCxcDvhRCnz9YEzxTNo6ARL2vXrmX8
   +PFqGV2QSLH6tLQ0ysrKyMvLUxdJUFziSUlJMSseAoEAXV1dEQWWQnfnPp9PLeWLJIzT0tKi5jJE
   yovA6VRCDL/9LVxzjZKDcM01UecV3LGVl++noGAhVmtGXDu2L7/8kjlz5oS9n8HITC8vLyc3N1fN
   B+mr8uRizoY/cOAAeXl56u/N6dOnyc7OPuutsy9GLuaf87nknOsoSJI0HtgMpANbgWMo4Yd/Ah6V
   JGmuEKJyMCakoREPg/VlEeqiv/HGGyPek5yczDXXjMbpdOJwOEhPT8dms9PY6Keuri4swbGuro7s
   7OyYhsKxY8eorKzk1ltvDTsfqklgMLg5enQrjzzyCCkpKb0WSLfbzaeffspjjz3Wuy2zwwF/+Qv8
   7ncwaxZ88QVEUV4EZbfe0dGBxWLB6zUxY8ZSDAbFE9LZaVI1GIJUVlbidrtVEZ+bb76515gGg2FA
   P5eess/BXTUQ00iI1RHyQkSJt3vVLpU+ny8s7yU3N/d8Te2CZ6C/WxoDJ96qh4+BScBiIUR1yPk8
   YC1QJoS4+2xN8kzQPAqXHoO1KAQCAV5++WW+/e1vR+y70PNZbnczVqvMggULwmKlNtspdLrmmDtv
   WZbVCoRIjYeCPQgMhjySk9NV/YPQFtQOhwO/36/u3Ht5ETo74c9/ht//Hm64QenmOHVqn59DZWUl
   J06cYObMmWFtsENjwKFVFY2Njfj9/j7zKvpLRUUFFRUV3Hbbbf163cUYu966dSsmk6nP/BUNjYFy
   PnIU5gM/CzUSAIQQJ1EqHuJrCaehcYaEtgmOVYUQDzqdjocffpjExER8Ph8OhyNsnK6uLtat243f
   P5SsrIkMGTIDvT4Hh8OhZl/7/T4OH16Fy6Us5A6HA1mWe4318ccfU19fj8/nw+l04vf7w+bidrtp
   amqmqmozQPfY1rCM7oMHD3L8+HH1WDUSOjrg+eeVHITSUvj6a1ixIqqREAgE2LNnj6pnkJk5BL0+
   h+3bT+F0Omlu3h0WA/Z4PLz00kvq/Tk5OYNiJAQCAbVSBCAvL4+FC/sv8HoxZMPX1NTwySefqMdz
   5szRjASNi4Z4zW0T0BnlWmf39TNGkqR/AZ4EAsBB4AlgGEp/iXRgL/CIEMInSZIJeB24GmgBHhBC
   1AzGPDQuXIKLQlraN4tCJPd4NJQW0Hu54YYbkCSJpKSkiB4Km62Nd9/dQWvrafT6BubOLWL48Dw6
   O5Vfdb3ehSy7sFisXH31QzQ372HnzmO0tXk5ebKe/PyRJCVJTJo0lNzcXO6//37a2zvYuLEs7DkW
   iwWLxYLZbCY7OwuzeRygJE0aDG61Rh5g1qxZ4W+mrQ3++EclUXHJEti4Ualm6ANJkmhra1O7SwYN
   r7S0b7QQpkxJJyUlRa2iWLp06aDHy/1+Pzt37mTkyJHo9frIctJxEJpwGvQo9BSHOte43W527drF
   9ddfDyihhGAjJg2Ni414//JLgX/o7uugIilbmh90Xz8jJEkaBvwDMF0IMQXFiHkQ+G/gd0KICUAb
   iiFB9/9tQohxwB+A3/YeVeNSI1oVQryLgsVi6ZV019NDsX37UT74YBepqYsxm29GiCV8/fVB2trs
   GI0empqaaG4+omZfO53lSBKYzeNpbbWSlLSAqioPu3Zt4/XX97B+/X5sNltET8g777xDa2srBoOB
   q6/Ow+M5pu7mp04dzvbt2/F4POEeD5sNnn0WrrgCjh+HrVvhzTdjGgllZWWUlSlVDJIksXDhQqxW
   a8TdeCCQwLZt22hoaFBfHyvnoj+sX7+e+vp6QFGIfOCBB+KWxo7GhZIN39TUpHpdgloTwWOTydSv
   7qAaGhcS8f4l/QL4HCiXJOk9FA2FXOA+YBxwa4zX9gc9kNjd4toKnEIJawQ75iwHngVeBO7o/jfA
   ByiNqzQucYKLwt69FXR2frMzj7UouN1uHA4HGRkZGI3GsJbKPT0UZrOZQ4d24nCMYPjwIRiNJqqr
   T9PcbKOtrZQFC6aRlpbGqFGj0Ov1uN1ufD4f27efQq/X4fUaSU1NpqzMz8SJSwkEbOh0WZSUHAMS
   SU214PG4VE/Ivffeqy4g6enpTJsWwOl0MnToUAwGAw8++KDq8aDVyZhP/0belx+iu+su2LFDMRai
   EJrDMGTIkIgegaDh1d7ehCy3k5o6AqPRw+233z4ofRICAeX9BJP2JkyYcFaSDBVxqORzng0f+hmv
   X7+eG2+8kYyMDHQ6nepN0NC42InLoyCEWA3chhJm+ClKt8j/BLpQWkyvPdOJCCFOAb9D6UpZj9I/
   Yi/QJoQIisPXAcH+qMOB2u7X+oE2SZLObiN7jQuCSIqBsTh+/DiHDh2KeK2nh8LtdjNu3GQSEnx0
   dtpJTExh+PAMxo5tY9y4TNLT05EkSV2MEhMTSUxMxGj0UFW1HZ+vEru9CUnSYTIlYDB4SE5OAxIR
   wkFDQwUHD36uekJ6KurV1dVhs9nUhc7n83FoQykTXn2Dou/chrXTz/Y/vonvxRdjGgkOh4NXXnkl
   LK8g2KUxlKDh5XRWcvr0QXU3PljNlMrKyti2bZt6PHz48LOmIhj8eZwrI2HdunUcOHBAPV66dCkZ
   GdpXkMalR7+UGQEkSUpAyRewCyFiq8T0b9w04EMUL0U7Sk+Jj4BnhBDju+8ZAXwhhJgqSdIh4MZu
   AwNJko4BM4QQ9h7jimeffVY9njdvHvPmzRusaWtcoAQCASRJiqtLYU1NDRs27CQnZ5LqoWhra+OD
   D/bididiNjtYvHgshYWFURdQu93Ojh0VdHT4qK6ux+GAzMw8pkwZjySBz3eSqVOHs39/PR6PEbPZ
   x/Tpo0hNTaWiooIJEyb0nmtjI57/+i/466u03fQwjY/9CM/QvIjKkKAYOXq9Xl0o29raSEtLizhf
   IQSfffYZixcvxmq1Dlq5qdvtZseOHWoOSER9h4uUuro67HY7kydPBpRkV6vVesahEw2NwaC4uJji
   4mL1+Lnnnju3Ogo90HX/N9h//YuAKiGEDdSSzNlAmiRJum6vwgiUcAQo3oWRwClJkvRASk8jIcjP
   f/7zQZ6qxoXOF198wahRoygs7LtELjc3l8LCfAoKRqr9EZKSkigsLOeqq6aqPQ564vP5eO+997j3
   3ntJTk5m9uwCAMzmq2lpaeHw4UaczjoqKzezaNENZGVlcf31yXR0dJCSkoLFYsHv93P06FHy8/O/
   ybNoaFBEkpYvx/Dgg2x7eQXSyLl9JuqtWrWKyZMnc0W3pyGSkRAIBAgEArjdbsaPH68ucmdSmx5q
   lJlMJkwmk2ogXMxGgs/no6mpSa3wMJvNYd6QYDhFQ+NCoOcm+Lnnnhu0sfujzLgE+DUwDcVIECih
   gZ8KIb4644lI0rXAX4EZgBt4FdgF3AB8JIR4T5KkvwD7hRAvSJL0A2CSEOIHkiQtBe4UQiyNMK6m
   o3CZYbPZ2bmzAkjCZPJF1Fno6OjA6/WSmZkZcQwhBPv372fy5MlRd4w+n4/q6mqSk1M4dKghrJrB
   bDbT1eVg167jBAJWEhIgPz+FEyc6aG93odPJFBVNDZ9XXZ1iILz5Jjz2GPz7v8PQoRGVIdPT05Fl
   mebmZkaOHAlE6e8QQklJCTabDaNx6KAKE7333ntce+21ajfDi5lQD0h7eztfffUV995773melYZG
   /xlMHYV4BZeWAF+gKDK+A5wGhgIPAFcAtwySsfAssBTwAvuAp1C8CMHyyH3At4UQXkmSzMAbwFUo
   ctJLe+o8dI+pGQqXAYFAgOLiYmbNmsX27ZV9iu/s37+fQCAQltgYT6Oh6upqGhoamDChgH37apBl
   HYcPVzBp0nVkZw9Dll20tJRgt58gK6uApKRCkpKS6erqpKRkLddeeyMORwMuVxdpaXplXqdOKY2a
   3n0XnnwSfvhD6KHMFyk00NTUxL59+1iyZEnEuQohaGlpUXMTHA4H27ZVkJBQcEbCRDabDYfDoRoo
   sixH9LhcbAQCAf70pz+xbNmyS+L9aFzenHMJZxRRpbUoiYvBxEIkSQpWQzwHnLGhIIR4rnusUE4A
   MyPc6wbuP9Nnapw9zqUmu06nIy0tLabOgt/vV132UyOIETkcDpqbm2MaCsGM9mCpo9kMOp2BqqpT
   pKVlYrFYcbvNJCWNYt++NjIzGxg/3kdLy1FkOQG9XkdOzngAug6tQ3znO/Dpp7BsGRw5AtnZEZ8b
   bBe9ZcsWZsyYgcViITs7O6qRAOByufjiiy949NFHVU+D328JK4XsjwZFkPb2dux2u2ooXMyL6rZt
   2xg3bhxZWVnodDqefPLJi/r9aGicDeLVUZgK/CnUSADoPv4zSjhCQ0PFZrOzcWMZmzfXsHFjGXZ7
   xPSRM0IIQWNjo3o8ffp0kpOTo+osvP/++9TV1YWpJtpsNrxeLwDDhg1j5sxeNimbNm2is1PRG0tJ
   SSEzM1PVHjAajSQk6Glr209T0zG6ujo5ebKe9PSryMwcCmRx9GgTDkcrZnMnfn8AU91xRjz7GLP/
   4V70OTlw9KjiUYhiJASRJInExMReqo6hHD58mLa2NgASEhJ4/PHHVSNhoBoUTqeT119/Xa2gyM/P
   Z/r06TFfc6Fit9tpbW1Vj7Ozs8Pe/9mqyNDQuJiJ11BwA9HUQpK7r2toAIMrsxyLjo4Ovv76a0JD
   S7HEdxYuXMT+/XUUF1epxsu2bds4depUjKdAZmZmWFJecMHt6FCEkvLzUzEaUzEYAnR2lpGfPxJw
   kZMDUIPNdoKcnEwevGYI455/mAmPXoMjw4KjdC+BX/4SR3fVQSSOHj3Kli1b1OPp06dH7UsBShgg
   mnRxf4SJWltbVQMqISGBxYsXX7SJiaGGVU1NjSr4BHDFFVdoQkgaGn3Qn6ZQk1GaQp0IOT8KJeSg
   NYXSUHE4HGzeXENW1jdKgdFK+gZCPCV3Pp8Pu93O2rVrWbp0Ka2tNpYv34ROdyUWS4CxY4dgNJ6O
   GJ/v6uqisrIyLH+hJ9XV1Xz00UomTlyEyeRlypRhJCQkIISguPgQgUAaej0kJQ0ncPgLbti8Ev1X
   X+F/+mnkZcsw5+TQ0dEZsblVaDvljo4O/H5/1ITD+vp6ysvLWbRoUdyfXzwhoU8//ZQZM2YMeuOn
   c83JkyfZtm0bDz74YN83a2hcQpyPZMbxKO2lU4EdfKPMOAtFVvn6C7XNtGYonHsGo5tftMXs8OHD
   nDx5MmJr456vEUJQX19Pbm4uq1fvpqLCQmZmAXV1L5Oevoi8PB/z5o3pZby4XC727NnTS1nP7Xaj
   0+lULQWPx4PX68VsNqPT6fjqq3XodFm0t3s5caKWyQYP0z57neyDu9H9y7/A3/89dO9eo31Gs2eP
   55VXXuF73/teVM0Gn8+nfiayLGOz2c54Qa+qqqKjo4Np0y7uKKLb7WbNmjXcfvvtSJKE3+9HCHHB
   dpHU0DhbnPPukUKICmAK8EfADEwHLMD/A6ZdqEaCxvnhTLX3Y+U3jBs3jqKioqivWbPmIB9++BV2
   ux1JkhgxYkS3Kz4Bq1XC7xeMHPk9XC4Ahxqfrq2txWazAWC1WiPK765bt47Kym9+1U0mk6oEGAgE
   aGzsxGTKZ6LfymOr3mL+L54mc/5sdFVV8B//oRoJEN7x0ONxotMJvF4TgUCAZcuWRTUShBC89NJL
   as6ExWIZsJEQGqJISUlhyJAhAxonUufNc0lFRYX6bJPJxNixY9VwVKgAlYaGxsDotzLjxYbmUTh/
   DKTqIdJOu6bma+bOncaIESP6fI3X20lLSw3p6UbVg3H48GG2bt3LiBELqK5uxuUS+P3HefzxG9TS
   wT179pCZmcno0aPDxu7q6lKFdXrqFNTW1tLR0UFhYSEOh4PS19Zw5YdvkXRgG43f/iHl8+Yze3FB
   xHBL6Jyrq7djtaaTmqqLGgrx+XyqgJLH48FkOrOGrV1dXSxfvpwf/OAHZ5R7EKnz5tno5RBKUDQq
   +DmtWbOGa6+99qw/V0PjYuKceBQkSVogSZImPaYxYAaivR+pmyFY1Yz7SDQ2NiLLEhaLleTkbPLz
   r8HrNak75rFjx3LzzfMwmZoZNUrH+PEyjz02F4fDoY5x9dVX9zISHA4Hb731lvrsnmJGJpNJKaXb
   swfrgw9yzbPfwz55Fgc/Pc7J+55Gl2yIWFHQ2dnJwYMHmTJlGLJcQWbmUFJTdVG9LuXl5Zw4oaYG
   DdhIKC0txelUVNeTkpL43ve+d0ZGwrlKWu3J6tWr1U6YAEuWLNGMBA2Ns0isb/CvgOuAEoDuFtPF
   wJNaqEHjbBGsKGhpqSUzcwRut0x29pCY7vUDBw7gchmQZaUro1KiuIuxYw2qNPKIESPIzc1VPRxK
   0mExw4cPR5KksPOBQACz2UxiYiLLli1TDQQhBGvXrmX+/PmYTCZyTp4k5xe/gNJSdD/6Ec4XXuBE
   eTPezpNRu1rabHZKSiqpqamitVVi6lSlSVKo18XlcrF//35mzZoFwIwZMwb8eYYmfnq9Xtxut1oC
   eKY9CmJpVgymu7+qqorTp08ze/ZsAG688UYtnKChcQ6JGnrobvU8SwgRNBT0KIqJ1wgh9p67KZ4Z
   Wujh4sNms/HGG+8xevS1JCUZIrqzg0lqwQW+s7OTvXtrsNncHD58FL9fxmrNYuLETIqKJpGeno7H
   48HpdIb1QOjpOpflOoYNGxa14uHgwYNMsNkw/eY3UFYGP/4x/N3fQbdIT7RwS3FxMVdeeSUHDzb0
   meTp8/nYsWMHc+bMiWvHH+2Z+/fvp6WlhYULF8b3wfeTwUhajYQsy1RXV1NQoPTOaG9vx+VykdtD
   rVJDQyM650OZUUPjrCHLMh0dHSQnJ2O1WsnIyOAHP3gKj8fTa/Hz+Xw4nU6WL1/O+PGzECIJo9HD
   pEm5JCXJOBxekpPHkJY2HSEkTp6sYNeuKhYsmMqRI0ew2+1qMmTQdS5Jw8jKykKWXQghmDx5kvq8
   6upqWlpauOaaa2DzZib/4heIykrcP/wh+g8/xNBDoCe0uVLobj6YXxFtB75u3TomTpxIXl4eBoMh
   YjJlJEINHb3eRW6uUe1uOH78eCZOnNjHCPERyRgJJq3u3VtBZ6cpqhclHkJloAOBQJihkJqaSmpq
   6qC8Dw0Njf6jGQoa55Xjx0/w0Ud7kWUjcJSHH76HMWPGYDQae2X+BxdFWdbhcmUTCAxX+yvs31+O
   z9eK35+F0ZiEyaQsxl1dBhwOxeU+ZcqUsPHcbjeyLHH06Adcf/2y7oXbHOY6T01JIaGkBP7t36Cm
   hq5//mdKxl+HRyRi3F4ZNXmvqqqKsrIybr/9dkAR9vH5fBiNDSEhkjZVGXH69On9jrOH5gikpVnp
   6LCxadNHFBQUYDQasVqt/RovGrESFtPT0ykqSj4jqW4hBC+//DJPPfUUVquVhIQEbrrppkGZu4aG
   xpnTV3nkcEmSxkiSNAYY0/Nc6H9neZ4alyCyLPPRR3tJSlpAfv63SE+/nw8/3Icsy73ubWho4KOP
   VqHT5eH3ZyDEBKqq6mlvb8RgMOJ0CmbOnElioo5AoAOPx4Usu+jsXI9O1xWWVBj0SpjNZiwWwbXX
   PoZer+Q4GAxu1q5di9PhgHXrSL/zTrJ/+lN4/HF8hw6x88obMCcVRkzecyk1l4DiQZg/f37Yewgt
   G62s3MChQx+oO/Ds7OyoJZHRcLvdVFQcxO9XnpuSksHYsXPweDz9GicW8SQsDiRpdd26ddTW1gKK
   i/Tpp58eNMNGQ0NjcOnrL/uDCOc+iXLvmWVGaVw2BN3YVVVVuFwyQ4cqu9O0tDza2srp6Ojo1ZjH
   YrHg8Vj55JM9uN1mGhvrKCz0YbG00NWVgcvVTEICjBuXQWfnSSoqvkKSdBQWTuGGGyaHLWI7d+5E
   p9Nx3XXXdbvOK+nqUnbLV08fhevTfVgWLwabDf7zP2HpUjAYcDscUUMHer2eV199lUf3oJIdAAAg
   AElEQVQffZSkpCRMJlNYdYIQghMnTpCfn09RUTIulwuzeU6/KxiCn11iYiJms5nc3GEIocPn89HZ
   2YZO5+yzd0N/GKyExVOnThEIBNQQzKRJk8JyRWK1x9bQ0Di/xPpLf+KczULjsiHUje1wNCNEO52d
   dpKT0+nstGM2O1Tt/fr6egwGAzk5OUiSxO7d1aSl3cmwYfkYjSfZu/dv6HSFZGYOZ+bMGzCZjBw4
   UEJnZzU/+MFdAOpOt6OjQx139uzZau5Aeno6eXmtVJ84wTyHC/3N3yG9qwt+9jO47z4IqQwIbapk
   sVhpa2vC57NhNo9DkiS++93vxqwkKC0tJScnh8TERJKTkwf0+e3Zswefz8ecOXMwGAzMm3cVxcUH
   OXLEjiTpKChIpbOzc9DKBXu+53gbSQkhcDgcqgaFw+EIK3HVEhM1NC4eNMEljXOGz+fjq6/2kJAw
   gaSkNGTZxYkTG6mr68DnS8FsdnDvvdPJz88H4NChQ5jNZjIzh7B+fSkffHCEkSPd+HxTMJsz8XoP
   MXlyDgUFRRgMBiRJoqnpMLNmDQ1bKP1+Py+//DKPP/547xbCQiCvWIH++ecx+v3wzDNwzz34AoGI
   cXe73c7evYqh09paSX7+ELVsryeVlZXodDrGjh074M+sq6uL/fv3M2fOHHw+n5r0F5yTz+dj/fr9
   6PX5JCWl4PN5B6XyIJTQ9xyvqNLJkycpKSnhvvvuG5Q5aGho9A+t6kHjoiE0W97tdlNbe4Jhw1JI
   SkrDYrEyZEgeN9+cjd/vJyUlBbfbrVYLTJo0CZ/PR3HxIdLTp5GTY8Pvn4FO52Do0HRk+Sg5OcmU
   lX1JdvYVZGTkYzJ5SU5Oxu/3I8syiYmJ6PV6vvvd76peBCEEK957j1t9PhJ/9zssQigGwp13gk4X
   NXlPlmWOHDlCUdGM7vc0LuZibLFYBiRo5Pf70el0SJKExWLBaDTS2mqjtLS215zcbjeBQAKZmRmA
   ki8w2FoG8SQsut1u3n33XR555BF0Oh2jRo1i1KhRg/J8DQ2N84sWGNQ4awT7L2zadJKNG8twOp2M
   GVNAevpoALq6OvH720lMTCQ7Oxuz2cynn36q9lwAaGlpYf/+LaSmplFUNAVZ3kZTUwldXRu4//4Z
   zJo1nuHDczEYCOspsX//fnbt2qWOoy7YgQDSBx9w1y9+QcIf/gA//zns2wd33w06XczkPaPRiM/n
   Q6fTRUze6+joYPny5WqfgZEjR0aVnY7F66+/TnNzMxBMgJxOaWltxDmFhgaAuEMD/SVSwuLu3bvV
   xFOz2cyNN96ofs6SJF20bak1NDTC0UIPGmeFoBiP2TyO/fs/YMKExUhSE1OmDOPAgVPYbG6qq2sZ
   PXo46ekWrr46j/T09F4tpH0+H2vW7EKnG0N1dTttbQ5keTdDh/pZtmwZer1e9Vq43W4yMjIiuuir
   Kivp+OtfmfbZZ5CUpHgQbrkFeixmPVtkV1RsRKcT3HbbjIg9G2RZxmg0otfrEULQ0tKi9o+Il6am
   Jvx+P0OHDgVQd+7R5gThbbsHEhoYCLIsI4RQqxN27txJQUGBpnGgoXEBooUeNC54vsmWT2Dy5NtI
   SEinpaWNhIQE5syZwJo1u5k8eQ4Wi4nS0o8AmDcvGYPBQHl5OZ2dnRQWFmI2m7n22it47bVN6PVj
   SUvTMXr0rfh8x9Wdu8FgQK/X884777BkyU1UVLTi9ZowGNxcPXU46WvWMPqXvySQmgq//z3ceGMv
   AyGI4lqX1eS9tLQ8JKm5V3llcDFfuXIlM2fOJC8vD0mS+m0kgKJEGQgEVEOhpzegr4TCwdAyiIct
   W7aQm5vLpEmKINXMmTPPynM0NDQuLDRDQWPQkWWZ0tJSDIZEZNlFYmJG2OLW0HCagwebMRqNGAxe
   JkxYgM/nVuPqZrOFY8dasdlqMBo9TJiQSU6OBWhl7Ng5WCwWmpvbVd2C5ORkJEniscceY9Omw1gs
   40lNNND2p6cxrv8SkZ+H7k9/QrdwYVQDIUhTUxO1taWMGqWjs9OEyeRh+vQr1MW3traOXbsqSEgY
   itHoYeHChWRmZvbr87Hb7axevZqlS5ciSZKqQBiNeBQQQxUhB4uKigqqqqpU8aNFixYN6vgaGhoX
   B1roQWPQ8Xq9lJSUMHr0aHbsqMBoTCMxUemOmJyczFtvfcbu3ac5dSoXCJCWZmfmTAsPPXQnFoul
   V/+AhoYdHDt2Cr1+PElJqeTnp6LT1ZGRIWhqauKWW24Bul30649RULKboa8+j3NILofufpQp//gA
   qSE1+z2x2+2kpaUFXXVqOKHnDt3n8/H551twOAQFBfP71dugsbGRrKwsdDodQijzzsnJ6dfnOpC2
   3f3B4XBw6NAh1VPgcrnw+/1qiaOGhsbFwzlpM62h0R/8fj/t7e0AGI1GcnOH8cYbW9m3z05paRmj
   RiWQnp6Ow+GgtlZGp8vDaJyEyXQVTU06vF4fHo9HDVmYzRYOHfqSQMDLsWMOJkyYiV7voa7uEJs2
   fURnZyfNzQno9dnU1dXhczpx/N//y+wnFpG2+m3K/s//8voTf2GLcSg7S45jt9ujzn3VqlW0trYC
   yh+X1WpVk/d0Oh2rVq1SOy9aLDkUFCiKixaLNayddSw2btxIW1ub+oz+GgkwMAXEWAghaGxsVEM4
   RqMRv9+vXrdarZqRoKGhoRkKGoNDVVUVW7duBZTQwxtvbKK+PpuurtE0N4/g7bfXq4ZEV5eP+vpE
   kpM7AUF2dia5uRMxm81qPN7tlsnMHI3L5cTn89LUdJqNG0twOivYuvU4LS0SZvNQKg46Kfv75/Hm
   jyVt/Xrcf32Fkl//lrWyH3Axa9ZcEhOvDJMdttlsqnwwwEMPPcSQIUPC3k9w8QyW+gkh+lVhUF5e
   zuHDh9Xj+++/n4yMjEH7vM+EUA/b6tWr1RCOyWSKqgmhoaFx+XJB5ShIkpQKvAJMAgLA3wEVwHtA
   HlAN3C+EaO++/4/AzYADeFwIUXoepn1ZE6xSGDduHFdccQWgLMRHj7ah11+FTpdMIGBAkuxUV1cz
   YcIEzGY/en0zyck1dHVNRZIcJCQoru8jR44wffoY9u6twGAw0d5+GKfzNH/7234SE++gri6dQOAk
   a1fuYJx5JY+tf5WW4YUc+tlyOibmUFRUyEyHA6cThg6dqu6+Q7UF2tvbsdvtjBw5EqBXGd+2bdvQ
   6XTMmjULgIKCAtXlHytXILQDYnp6+gVZHrhy5UomTpzIuHHj1LwODQ0NjVhcUIYC8P+AVUKI+yRJ
   MgCJwH8A64QQv5Uk6UfAT4AfS5J0MzBWCDFOkqSZwAvArPM288uQ4uJiMjIy1K6MwYVRr9fT1uYi
   OzuThIThdHW1UlkJQ4YMweFwMGlSAYmJMsePt5KefpwhQ2Dy5OEkJCSg0+nULH6Hw8HatTVkZBjx
   eq10drbjaj3Bw44veaDkfWyjC1n3g/dpHDmK6dPH4LVXqn0QkpMN+HxeDAYDDkcntbW7MRiU3mX5
   +fmq+iMoxk5HR4da5jd16lTVSxBJfKmoqLBXrkBLSwuffvopTz75JHDhSBRXVVXh9XqZMGECAPPn
   z49Y5qmhoaERjQsmmVGSpGSgVAgxtsf5I0CREKJRkqRcYIMQYqIkSS90//u97vvKgXlCiMYer9eS
   Gc8S7e3tJCQk9Op62N7ezq9//SF6fSvHji0AHOTnV7FkyUgsFguynITBMAYhAjQ0lGO1ysyff5Ua
   f29rayMlJQWXy8WmTScxGofy/E//mwdsOu6ve5na4SP5rbmAMffcRmbmJKZMGY3JZAxLLOypLWC1
   Orj22mvD4vvB5EBZlvn000954oknemk49EysDD5Dr9dTUlLC9OnT1ffv9/tj9no4F3g8Hmw2m2qo
   1NfX4/P5yMvLO6/z0tDQOLdcqjoKY4AWSZJeBaYCu4F/BnKCi78Q4rQkSdnd9w8HakNeX999LsxQ
   0Bg8hBCUlJQwbdo0zGZzVKEds9nM1VeP4PjxUUyc2IbBEKCwMI+5c5VuicoiXtXdGKoa/n97dx4e
   1XEnev9bvS9qSa19RxsSWkAgmcXYGDAOcQyxHWfBubbjOM5kJpPMzJNk7kzm5r03fjJv5k3uk9yJ
   k8zYM5kk1+OJHa8k4DheYpt9sUGAEBJaENrQvrV6V5/uev/oVlsCCQPGgE19nkcPR6Wj1uniSFWn
   6lf1I5H9+/vQ6XxUVmayZ88eVq1ahclkQucdpOCPv+H3jY9zyFbO10u/SaCsgHXrEvj4xyvo7HQT
   DJ4hEpk9DdDU1ERSkpmKitI5VwkcOnSIkZEI4MBonOLuu+8+Z6pgrsyJk5PG+KhFJJYPYrqjcLU6
   CTM3qRodHeXw4cNs3rwZgNzc3KtyTYqifHRcSx0FA1ALfE1KeUgI8c/At4H5hgPm6impoYMP0PTy
   wVAoNGcAX3d3N++88w7JyaV4vdDcfJicHI3+/gQWLChiaGiI9PR0RkZGWLu2Cq/Xy8GDPkymhYyN
   9dPZKTh8uJ7MTCtP/9trrDq0j5v3vUjXwhKeevg7uPLWsCnPSUZGGkL0UlxcTHHxuzsZzozYr66u
   xmq1zpnGWdM0TpzoJCtrJWlp+QQCfo4caWXt2kQMBkN8pEGv18/a6KilZSeh0DBmcxkAN9544wdX
   2RdoamqKxx9/nK997Wvo9Xqys7PjnQRFUZTL4VrqKPQCPVLKQ7HPXyDaURgUQmTOmHoYmnF+/ozv
   zwP65nrhRx55JH68bt061q1bd3mv/EPsQtbmT0xMkBzbh2A6wG+u14lOHaRhMBTj82kkJ2/E5/NQ
   UpLN/v1H6Ol5h9rabDyebkpKSmLTDCH6+g4xNHSCqamVZJh0pL7wKz7X9BZDNbez7Vuv0Wnt5HOf
   W86pUxNEIhGE6J01ejA9XfHss8/yZ3/2Zwghzhnt6Onpoaenh9WrVxMMBsnMrCItLXr7WCzWeLDj
   5KQ7HpMghJekJA2vt5WJCR02m2DVqps/sJ0PL9Sbb77J8uXLcTgcmEwmvvzlL1/1KQ9FUa6uHTt2
   sGPHjg/kta+ZGAUAIcRO4M+klK1CiO8CttiXxqSUPxRCfBtIllJ+WwhxB9HRh01CiFXAT6SU57Ri
   KkZhfvNlSZzJ5/Px1FNP8dBDD83ZGDU1NaHT6enu9jA2NkV7exdOp8TvX0RXl42BgT6ECKPTjVJV
   lcGNNy5GiDMsW5aPyWTiscdexWS6idG2YdYefZIVhx7jROEyDtz6CRx138BsttDZ+Qpf+UotKSkp
   s/Ig9Pb2kpGRgc0WvU00TZvViEciEXS66Apgt9vNxMQE+fn588Ye3HRTOXv3tsTLx8cHOHXqVT75
   yY0cOdIF2LFYIh9YLoX5jIyMYDAY4p215uZmCgoKVFCioijz+qjGKAD8NfAbIYQR6AAeAvTAs0KI
   LwHdwGcBpJQvCyHuEEK0E10e+dBVuuYPpZlZEpOTo41lfX0ra9c6ZjW2NpuNhx9+eN6lfkII6utP
   MTycByTS3z+KlFPY7UHOnBlAr88lJSU60DM0NIzVmoDbbeLNN9+kurqa8lQHOc98l4pdT9NYuoEn
   /uLfafAFsRqXUWsw4HaPYzZ7SUxMjG9TPN3BaWtrJCsrl7Vra3A6ned0Eh5//HEeeughrFYrDocD
   h8MBzL8lsqZpNDcfprY2uszT6cwiL6+W+vounM6aeKdirnq63GZ2ctrb23E6nfGOQkVFxfm+VVEU
   5bK6pkYUPghqROFcmqYxNjbGO+8MkpW1OF4+nZFweHiYxsbGeee6z5w5Q05ODkIIXC4XP//5Szid
   y0lMLGNsbJiGhq3U1qawf/8ZTKYFjI66ycsrIjc3xPLl0RGFWyrS0f/kJ4Qf/3eG193Dq0srGHHU
   oWkecnKgs9OFXp+O2ezlM5+ppaioiNHRUfr6+hgaisy5EmFqaopIJBLfTdDr9Z73qVvTNNxuN3q9
   nsTERDRNY9u2naSlrcBudxAI+JmYOAbY56ynD+qJvrm5mdbWVu66664P5PUVRfno+yiPKCgfsOmn
   8UBAR1NTK3p9BunpmbN2GczNzZ13RYOUkj179nDHHXfEn9ABhIhOSyQkJFBYmIYQg3zmM0sxm8sJ
   hcK0tnbg9R5BPzLFjXu2Y/jsU4TuuQf3zrc40hckXVrJNWlUVpaRlZWFpmlMTk6SmJgY38QIwOPx
   EAolz1qJMB1fcPjwYZKSkqipqQF4z4bcYDBw9OhRkpOTWbZsGQaDgfXra6mv72B4ODrSsGJFCQ0N
   ffNmbrwc3G43e/fujSdfWrhwIWVlZZft9RVFUd4PNaJwHTl7bn5oqI/Gxv1UVpYRCAxQXV0Q35jn
   bKFQKL4MUEpJW1sbpaWlRCIRXnppP6dO6ZHSgE4XprhYY9OmVQwODnLixABS2rGMd1P3xnMkvPgi
   4v77GfjCF3itqYnc3BoCAR1Selm5smRWmmYpJVu3buVjH/sYOp0u3jhPvwedTtLbe4KUFNMFJWaC
   aGBmS0tLPPHRzKWFM+tpZoDn2XsyvN8YBSklLS0tlJeXI4QgHA5z8uRJKisrr8ndHBVF+fBRIwrK
   JTl7X4CMjBwqK8tYvjyToaFoMqS5jIyMsG3btlkbErW2tpKdnY3D4WDJkjyamt4iGHRjtWZTU7Mc
   o9FIY2MjJWYzhc/+EsMzzzB1771EGhrQ5+WRpmnkTBhmxUg0NLSydq0TvV4fb0Bzc/PYu7eZ6f0O
   amsL4vEFfr/A7T7Fhg0bgehUw1yrN2ZuhGQ2m2ctmZyrYT47ZfP0TpHvJ3NjJBJBShl/b62treTn
   52O329Hr9VRVVV30ayqKolwJakThOnL2iMLExBBS9rFuXfU5jd/ZT9WBQIBAIBAPqINoAxyJRNi1
   qwm9vgghJD7fBCbTBGuLHOh/9CPEM8/Al74E3/oWL+7bx/Lly8nPz4+mhN7dTXr6u4F5w8PNOBxj
   hEIhliyp4dChDg4fPoPNVhjfffHEiefZvHkDGRkZ8eubuaTx7Cd+KSWPPfYYDzzwwKypkvdbjxfb
   aXjhhRdYsmQJCxcuvCzXoCiKcj6Xc0RBdRSuM9PD6FNTRlpbd7B588cpKZm1azZjY+O89to+gsEI
   eXnRJ/hAIMCuXbvYsmVL/LznnnuO6upqurogPb2CQMBD2+uPc+vbb5N/8A14+GFGv/hF0iorgej0
   xdTU1DlTCFJOIYSJQKCVlSuj+yvs3duClLk0Nw9jteYTCvVRW1tMT88Bbr21kpSUFDRNi23a1I7d
   XhmPIRgePsDKlYvIzs4GmLWk8v26kCWlEA1I9Pv91NbWAucu3VQURfkgqY7CRVAdhXNNPxGbTKZz
   8jRMjzpImYFOZyAUCuPzneSWWxbHlyhOD9cHAgGCwSD797eQPmlnwVP/TNJbW+nd/Cnyf/z/MqbT
   sWvXLj7zmc/M2cACvP12Ow0N+1i69BZuuKEQp9MZH21wOheyf//raJofs3kZlZVJQA833VTO8PAI
   zc2D+HyC5uZeVq1aQ2JitMFubn6FZcvyL/tw/vlyPwSDQfr7++MZNEdHR4lEIrNiLhRFUa4U1VG4
   CKqj8C6/389zzz3HfffdN+fmSceOHSMvL4933hnCZMrlwIHfMzDgp6VlnNraXJYty2R8vIOvfOUr
   8c7Cwf/6L4qeeorkvfvo3HQvpzZ9ktrbbiAzMzP+umc3sAMD7QgxzMaNy9E0DZfLRVJSEhaLBSkl
   PT09tLWNY7GU4fG4OHKkiUBglBUrFrBoUTrt7WPU1/djtaZTWVlMQ0MzweBJNmz4EpoWmpUc6nI6
   e7okFAowMXGaNWsKCAaDHDt2jPXr11/Wn6koinIpVDCjckmsVisbN26cd7tfn88Xmx4Yo7HRw8hI
   AcPDgpQUJ8GgkTNn/JSUmBkYGCDb5YLvf5/lf/wjnZs/x5F//w4hu56pqT6mpqZmve7ZQZQ+3whS
   hujvH6C1dTQ2yjBAbW0BiYmJ7N69mw0bNrB//yGOHDlDT88ImZmJeDyTNDQESUioxGyWJCRU0t7e
   xuLFFRw6NER//zEcDsOs7Z2nXUpcwdnMZnM894PRaGDnzseprb0Fs9mM3W5XnQRFUT6S1IjCR5zL
   5aKvr2/Wbn4zG82xsTHsdjtms5mJCRdPPfUcSUnFvPlmM5FIFnb7FGlpJQSDC8jKmmSx/iglT/+E
   4tOniXz96+xashafMZe0tELc7gkikS5uvbVmVmM8PDzMyy/vpLR0U3zI3uttAsBur8Tl6gV0mM2e
   +EiApmm8/vphjh8P4HTWotMZGR09Sig0xLp1G9m16z9JTf00weAQFRXpRCJdrFxZGk9VPdOFxhW8
   l+3bt1NVVcXp025CIRMGQ4C6usIrup2zoijKhVAjCsoFm964aNrMRtPt7qezs4Gqqo9hMATw+XwU
   F9+G3Z5GRYWR9vYzWK1ZaFoO2cP72Pz6L1jQdRjD3/0t4utfJ6DT4Xq9iYGeerq7vWiaCY/nDJWV
   meTl5cV/psPhYOHCHPz+d7dMrqzMpLHRjcFgZGpqCqs1gVAounGSwWAgGAwSDBqwWjOwWh1MTh5B
   iCk0LUgg4GPVqntpaOjE5+skEvGxfHnxnJtEXehW1XPp7u7GbDbHp1Fqa2tJT0+noED3vkcnFEVR
   PizUX7mPoKmpKaSUmM1mUlNTSU1NBd5tNM3mUiwWPcePjwDLcLtdpKUt4vDhV8nIsODx7MFoXIDT
   6SdvZB/rdn+T4oFu2u/ewuTWfyNgMJAS2wDJZtMRDC7A4SjDZBJEIiGamgbZt28f69evJz09HYvF
   wurVq+MjGT6fj7feeguDIYuWlhA6nYNIxMWCBS7M5ujyweh+ByEikUmmpvxAEnq9n+rqNMLh04RC
   NioqQlRW1pKVlTVvg332tMfMnRzP/h4pJT6fL76jo9vtnpW6Ojc3N36sOgiKolwv1F+7j6B9+/bh
   cDioq6ubVR4MBhkaGmZysh2/v4SOjhBGoyQry4DVaqe/301qajapqU4SW49y18v/TvFYD2Nfephd
   y9aDLZOxjknC4UHKyspIT0+nrCyV+vpGhOjHYAiwZEkhweAZVqxYEe+gQLTzYjQasdvtWK1WVq1a
   RX19F+AHjITDXoJBf3wZYV9fHyMjrSxYUERz8+sIoWPRoiTWrl2Kw3Hhmx/NjCt4ry2Ym5ub6ejo
   iOe4UJsgKYqiqBiFj6RIJEIkEjmnMdU0jT/9qZ7GxnFstmo6OnoJhXQsWOAnMdFHa2s/tSEjK1/7
   D1J7j9H+qQcp+cG3OXCsF7/fTH7+YgIBP0NDh7DbbUQiNnQ6H273JKGQFb9/hLKyDXOuOnjyySfZ
   sGEDOTk5ALOWQI6MDNPRMcb4+NssW1bOzTdXkJycjJSSSCSC1+sFmDP+4ELMtwWz2+3mpZde4t57
   752ez1NbKCuK8pGgYhSuY/NF77/44ovcfPPNZGRkMDHhmhW8ByOUlZWRn59PdXUOhw4dR6ezkp0t
   8Hq9HDnSxmJ3Cw+88zoF3nF67/8GPXc+R0jfB1YrU1N6XK5e8vKqMRiMnDw5zvLllej1QSyWQny+
   Q9jtepKSiggEWqmtjS4XHBkZISsrC4DPf/7zs653+knf7/fR3e1GiExMJjMmUzH19d3xGAKdTjdv
   gqoLNb0FcyAQoKGhgcTERCCawGrt2rXx81QnQVEU5Vyqo/AhMjY2zqFDHfh8YLPBsmUF2Gw2zGYz
   q1evJi0tbc7gvcHBkfiSxaysLFasWMrQkIecnCoGnn+az7z5C3L9w+xft4VfZNxIRpaJsonDrFxZ
   gd1ux2bTkZKyDiEEbvcEQuiw2xOor3+ZurrPYbGkceONORgMhngHpr29nYGBgXhH4eyRAIPBgN/f
   w+hoOz094PVCTk4xHR2DpKX554whuBQ+nw+9Xo/ZbCYhISE+0mKz2RBCxEc4FEVRlLmpqYcPCU3T
   2L59P93dJnS6RCYnu9G0durq1mC1yvhwutfrZdeuLjQNsrMrEEIwONhIf/9B7rvvPiwWC+NjYxz6
   8b+y+MVnMQ4MsH3xAzTf8FfozILExEGKikLk5UmKiorIzc2ND92Pjk5gMGjodCbS02+YtTvhTTeV
   s3v3blatWoXVap2zkff7/YyMjJCfnw9Epx/C4TA/+9l2kpI2kpycjts9jsfzJn/1V5tmpZe+VNu2
   baOqquqcbaoVRVE+ytTUw3XI6/Vy8uQ4mZkfQ6830tU1gRAmbLZ8zGZLfMlfdEg/yJkz3aSk5AFG
   DIYgmzdvpr2tjerBQZzf+x639fWxZ+1a9ixYw8BoFUkJSYTDTfT2DrJ0aS6rVtXFG/vpofvGxkYS
   EhJITU2lvv7dpY61tQV4vT4GB73s2tWJxRKZc6+CiYmJeNZEiMYceL1eyssXMjw8zvi4B4MhRFFR
   /qzVBhejoaGB8fHx+JTCnXfeeemVriiKoqiOwodJOCzx+ycwGm1omo1weBF6vQGLxcrkpJGxsTEy
   MjKoqytkYmKEI0eeJzd3IQLJ6cefp/DJR5mSYUzf+x5iyxZqg0EGXj6OqTtAZ+curNYBUlMTWLw4
   F4PBgNvtZufOnWzatAmDwcDSpUvj17J2rYNdu3aRlJSEw+Fg584TFBV9PD7KUF/fys0329i6dSv3
   3HMPJpOJ7OzseKKmaWazGafTTHp6Nnq9gXBYi+V2uLAkThMTE3R0dMSTL5WUlKhYA0VRlMtIdRQ+
   JILBKcLhIXp6nsftriYY7Ke0tCjWSZhgdPQUr7xyiE9/+tM4HA42bbqZttZWPM/+kcVbX8AQ9HH8
   7of5Y0U+D2z5HAaDgdMnT2I2h7jhhqXU1UUIBv2Ew50kJCSgaRp2u/2ctMg+nw+bzYbBYKCurg6r
   1Tprr4KpKT96vSAUMqFpGjfeeOO8W0ZDNFahtraA+voO/P53Ryjmi0+QUjI8PExGRkb8+2eOPkzv
   gaAoiqJcHipG4Ro3M5VyOJxNe/sZPB6NQKCZJUuKmZgIc/BgE3Z7Anr9ICUlpbCXtXEAAB1oSURB
   VGRlZpB9aC8lT/0HkUCQoT//J6Y23w86HcPDzaxZU4DdbqepqQmDwUBnpxufDzRtkr6+FtLTF5KY
   aGfFipJZ2Q9HRkbYvn07Dz300DnXOJ30qb19B4mJuSQmyotKzPReuRimly5qmsYTTzzBgw8+qDY9
   UhRFmYfKHnkRPswdhbGxcV5+eSfBoGRoSM+qVWuw2RyEQkFcrg6Mxn4OHmwnMfFT9PcPEdHMLDr5
   azYd+R1CwOhXH+HFyCh6YyVr1qxmaspLd/ebbNnyScLhMGazGZdrkv37mwkGjbS3d5KXV4HLZcLr
   nWJqqpHbbitm2bJlmEwmILpHg06ni1/j0NAQPT09FBdHlzVOTRkxmULz5lO4lORMzzzzDGvWrFEr
   FBRFUS6QCma8Dkwvc8zKWonJZMPj6aKhoZVVq+oQQofFEiEnJ5/9+9uxWaxUNb/JxoO/IiQm2HHr
   nzO08h4WFNoQjW00N3fg9w+zaFESKSkW9uw5iaaZ0el8DA724vMFqaraxKlTOt5+e4iiokrGxtyM
   jibw6qsHcDqdFBcXxxv4SCQSb+TNZjNWqzUe8Hi+TsCFJmc6efIkFouFwsJCAG6//fb43geKoijK
   laVGFK4h00/bHR0d5Obm8vbbg6SnR7M+Tk5OcuDATpxOD4WFpaxcWYpBCF798rdYt3snPr2F11b9
   AycKTKSl+4BsBgYmMJtt+P1u0tIyMBjGMRjCFBUVUVp6AxMTY7zzzn7Wr78dIQT79x+mtfUMiYlG
   hLgZKbtJSdFRUNBPSkoqkYgNvT5AY+MOvva1r8VHGS70vU1PT8xcVrl2bRWapuF2u+PTHJ2dnZhM
   JjWCoCiKcoku54iC7r1PUa6EsbFxdu48we7d3Rw+fJLR0dF4jgIAk8lIXV0OZWWZ5GToMT7zLFTX
   sPzAWzyz8ka+tfoBdiX0Ub3YSWFhFl7vKQKBRIzGFIzGfLzeDMzmQqTMp6urFa93AocjGZB4PJMY
   DAZKS7PRtA7c7kmk7CYnJwObTU9b2wChUApO50KEyGfRotWzph8uxHTAo8USTc5kNlvi2SL7+vpo
   aGiIn1tYWKg6CYqiKNeIa27qQQihAw4BvVLKO4UQhcBvASdQDzwgpdSEECbgP4E6YATYIqXsvjpX
   /f4EAgF27TqO07mM5GQHDkch3d2tLFmSw7FjJzlzxk1qahLLlxUgnn4a849+RCAxlbfu/A7DS+5l
   YnQ71ZqZkpLoLop9fb2YzZkUFppJSiqjt/c04fAgBoMVm81CMFiHXm9F00KUlydy5MgTVFZ+DKs1
   wje+sYkXX6zHaNRhNLrJzrbgcoUJBLzU15+IpZIeZPHigVmppN/LzORMoHHgwJPU1KzGbDZTXFxM
   cXHxB1fBiqIoyiW75joKwN8ATcD0pPQPgR9LKZ8TQjwGPAz8W+zfMSnlQiHEFuB/A/dejQt+P8bG
   xtmz5wTHjh0gPd1ORUUxNpsDtzuMyWSirq6Ql7Zu5ebRIOav/AgKCvjDli0MV65hYLCAic7TRCK1
   eDwnCQQGufXWTeTmpqBpAfx+jVCoHZerA73ew8KFH0cIQVPTAUZHfSQkGFi/voa1axdjsVjisQWb
   NgXZt+8dCgqWYzS6WbSokN5eLw7HklmppKe3Z36v4EQpJW+88QZLly6lsbGVUMhEVdUN1NUtUCsX
   FEVRrnHX1F9pIUQecAfwfeCbseJbgc/Hjp8Avku0o3BX7BjgeeDnV+5KL4/pgMWkpBrS09MQIpPD
   h5sRwkww2I3dOMXq1v18+oc/JFJaCr/+NdxyCxUdHbRvf432dhcpKZ/Aaj1FKKRnYMCDzebA4aig
   p+cp/H4rfr8Ok2mM0VE3Bw7Uk5wcYcECSV4e1NREly9KKenv7ycnJwdN0ygoKCAnJwe73Y7ZbGZg
   YIC2tqPodAMYDKF4Kun+/gFaW0fnDE4cHBzEbreTkJCAEILc3FySkpJYuzb1olc9KIqiKFfPtfaX
   +p+B/w4kAQghUoFxKWUk9vVeIDd2nAv0AEgpw0KICSFEipRy7Apf8yXp7+/n4MGDWCzlJCc7qKjI
   oanpDCdOtLCoyMCak7u54dGXcefmIh97jG0jIzy4ahWesXGOHTtDf7+N06d76Ot7gQULkikoyMTh
   SMfv9+B0plFUdBNe7xGWLFnNqVM9LFigR6ezUlOTTzjcSU1NdbyhDofDvPHGG2zc+HEaG/tnNfx2
   u52srCxqa7PR65NISEhE00KEQj6am33Y7ZXx5FOHD59k3ToHEN1Kubi4mISEBACqq6vj7111EBRF
   UT48rpm/2EKITcCglPKoEGLddHHsYyY542uzXmLG12Z55JFH4sfr1q1j3bp1c532gTp7/4D09HRq
   a2tpaxsnEPCTmJhI+YIJiv+wi3XbX6EnJZmfr/lbmhx2NgZsJCcn09/fT0vLCJ2dkJJioqbmbvz+
   MYqLJUlJIxiNOkKhU/T0dKHT+SktXUBSkhMph9Drx9Dp0jEazQSDZvbs2UNpaSlZWVkEg0E++9nP
   sndvSzzrpMfjZt++E2zYsBSLxcLy5dF9EsbHo52IyspMGhvd8eDEwcEmhodbKS9Pp7V1lHA4j85O
   Dykp43MugVQURVEunx07drBjx44P5LWvmeWRQoh/Au4HNMAKOIDfARuBLCllRAixCviulPITQohX
   YscHhRB6oF9KmTHH61715ZHT+we4XD70+iluuWVJvPEcHx/n6L4WTP/3X6h944+Ml5Tz5s1fpslq
   Q6dzEghkkZPjpbo6Gq/w+usn6e93YLNZCYWMtLT04HSOs3RpIgsXptHb6yUYNDIwcJyUlBwyM1fR
   1NSBlPlI2Ud19QI0rYOiIgcGg5G2trHYdssuAgEdRUUrcLkmaWkZYHS0g6VLk7nppnKcTueszs7Y
   2BgvvfQm5eV3YbFY8Xon8ftbYimoK89ZAqlGERRFUa6cj+TySCnl/5BSFkgpi4kGJb4ppbwfeAv4
   bOy0B4Hfx463xT4n9vU3r+T1XqjpOASLpQyTKYlIJJl33unA5XKhuVw4f/lL1n35U5QNdXPokf9F
   ZOszdDgT8Xjs6HS9lJYuQKdLoKnpOL///Z9obx+ivb0dvz+Cz7ef4uIEysuNmEwWXnrpNK2tkJRU
   RE7OCvR6PcFgK2lpflyu36NpB9G0DmprCygoKKCtbQyLpYz09AocjipOn+5hcnKClpYBIJ3U1Gwc
   jirq67vj+zvY7XYMBgNJSUnU1CwkEGhleLiZUOg0VVXZRCK2+CiDxWKNL4FUFEVRPpw+DI953wZ+
   K4T4R+AI8MtY+S+BJ4UQbcAo1+iKB5fLxdSUkeRkK9nZlbhc4xzZ9So5T/6aku2/Ibx+HebXXsNS
   WEj7889TZbNxww35NDYKkpPXEwi4OHnyOOGwh5ISG1VVlej13TQ2vkVRUQ4lJRPo9Ql4PCCED4dj
   CSdPtlJTU0ZLSwsbNtRiNBoRYgmRSASLxYLBYMDr9cYTOQEkJDgoLMxlbOwoo6MBUlP9lJXlYrfb
   GRkxMTU1xfHjxykqKsJoNGI2m6mrq5s1ygDQ0jJKIOCPjygYjVMXnAlSURRFufZckx0FKeVOYGfs
   +DSwco5zgsDnrvClXbRXX30VqzWbQKAIc9CP5f98h6/ufAbfjRt5/i++gWPVLWyuqMBhMJCYmEh/
   fz/Llxfh9zdy4sTjdHVZyctbhMFQjtnsord3lFtuWUNZWSIVFQ4yMzPZvbuLnJyFjI46iERA00x4
   vR4giMViwWKxnHNd0/saeDzueHrnlBQrK1dWYzY34nAU09j4O4qKVmM2T2G1Wvn0pz99zusYDIZZ
   0wrRTJCtuN3vnQlSURRFufZdMzEKH5SrHaMQDodx9/Qw9sj/Jnfr07QWVTP21e9jWbqazs6DjI93
   8cADd+H1+njnnVOAHYslwpIlOTQ2NtLY2EckUkN7+zhmczLp6X6WLi0gGDzFmTPHKSys5ejRAWy2
   QnS6EwSDqYRCAerqclm+vPi8gYQdHad54YV6gkE7ZrOXZcuSyMnJITs7OzbdoMdsDs+bk2E+l5L4
   SVEURbl8VFKoa1goFGLr1q3cfffdmDwe9I8+SvK//AuJd9zBkX97nD+0dpHgTcRaf4K0NBt+/zgT
   Ey6eemo/mmbEZBJUVa2goaGPiopyDh7sISkpl+rqQlpa9tLYeBhNK6CkpAS3OxmPx8yqVWtoaGjF
   47GwZEkSNTWVZGVlnbeR1jSNkycHWbiwgrS0BYTDEVyuI6Snp19QgqfzOXuUQVEURfnwumaCGT8q
   jEYjqxYuxPjII7BwIbK3l8j+/UR+9StGU0qx2xPR6wP4fCY6OzvYsuVTHDnSjV5fQlbWLQQCPZw6
   NUwgoMNgMFBSUgZ0Mz7eBHSi0xkxGEpJTCxGymy6u9uw2RysWlVHRUUBN9+8iLy8vHkb6unAwmAw
   iN8fJhRyYbXaSUhwYDKlY7PZgGhjPx24qCiKoly/VEfhMggEArS1tcHwMHz72+Tfdhtafz/BvXvZ
   dued7O7rIxgMEonYcDhSycszsXJlIeXlxbz44ot4PFMYDCECAUlm5ucIBPSAl1AoxPh4K2VlBbG0
   0uvJzTWTnFzF6dNjJCRkIEQpoVAQTQvhcEQb9/l4vV5+8YtfIKXEbDaTnJxATs7S2HtQgYeKoijK
   udTj4mUQ6OrC9D//J/zpTwTu/hSHfvYb+o0pdG4/QUJCiNOnz1BZWYnROEVV1WbMZjPBYJBw2IPB
   kEBTUyNdXQno9V2YTGaysydYseJ2nE4ny5ZV0N/fhMfjJjk5keTkEgwGPVNTFvLyBO3tTbhceiyW
   yDmBg6FQiBdeeIFPfOITJCUlYbfb+cu//EuEEBgMBhV4qCiKorwn1SpcokgkgtbdjenRR0l+4gmS
   77sPrb6ePe0TtLWdJBzOJiHhVkKhE4RCIxw92svSpXm8/PIf8Hp9lJZWIQTk5KzjyJGDTEwIIpFB
   6ur8lJYuwul0YjAY6O3tpbq6Gr1eT0JCPpFINg0Nx/H7h7Hbs/niF2/BZrPFYwlOnToVCzwUsSmN
   XA4ePMXy5UU4nc5Z6aHfbyyCoiiK8tGnph4uhZQM3nsvYsmS6OeNjfCznxFMT0fTzExODjA52YPN
   5kCIVFau/ALhsAWbzcYnPnELTqeRFStKMBqTOXWqF7t9IbW195CXtwKDoZBQKBKPJbj77rtZuHAh
   q1eXo2kdBIMDVFQYuP/+pdx6aw2pqanxEQKI7tvgdrvjmzwVF6/Bbq+gvr4bTdPOeSsqFkFRFEU5
   H9U6XAohyLzjDsSjj0J2NoFAgFMnTlBeXo7fP0IoZMXtHmVw8Aj5+V727v0DlZXLMZsXYrfb+fzn
   P4+maUjpxe8fw2abYHw8k4mJSSKRMFIO4vP5sNvthMNhgsEgDoeDtWurznn6P3ToEJOTk9x6660A
   1NbW4vV6aW72xTdTslisuN3RHRJVh0BRFEW5GGpE4RLpvvhFRHY2AKOjo2zbto1IJEIkomE0lpGY
   WIqm9TAx0UNmZj6VlTkYDAampqbYu3cver2elStLMBgCGI3NdHXtIiEhGYvFTm3tHTQ09DE0NMzO
   nSfYvbubnTtP4Ha78fv9bN++PX4ddXV18U7CtOnNlAIBP6ACFRVFUZRLpzoKlygUCuH1egFISkrC
   YDDQ1dVFa6ubSKQPKd/hllsqMBotJCWV0N3tY3x8HKPRyMDAAAMDAzgcDrZsWUlJSTUVFUYWLQqw
   efNKsrJyCQR0vP32KXS6BYyNDaPXF1Ff301CQgI33nhj/DqEOHc/jelAxek8DIFAqwpUVBRFUS6J
   2pnxEj377LN4PB6+9KUvATA5OcmBA200NxuxWHLRtDGamg5TWlpOZqYLqzWFxETJ4sU5vPLKLoRI
   ZmBggtRUBwMDY7jdIdLTS1m2bBFSani9bXi9krExO253Kw5HCenpIe64Y/F5l0DOpHZIVBRFuT59
   JLNHfthIKTGZTEC0Qf7pT3+Ky+Vj4cJMNG0In8+Hz+envDwHk8nC5GRffJRgwYKNeDz5WK1rOXbM
   TWLiBtLSqhHCwoEDuzlx4kUKCuz09AwgxAJycz+HEMV0dp5Br9df8DWqQEVFURTl/VItyCW68847
   mZycBKIN8n333c/zzx9Gp0tDrx9FyqMsWpSB3W7D6VxCIODH52tDCDt6vQFNM2Kz2ZmasqNpJ4lE
   PNxww82Mj0vWr1+LwWCgrS3I8PAw4+MTGAwhioryCYfDV/mdK4qiKNcT1VG4RC+99BIej4cHHngA
   gPr6dqqqVtLVNYrXayQQCLBlSw29vR34/dENjVasKKGhoY/JySE0rQ2/34HJ5EWIcmy2aLyDw6GP
   Ty04nWbS07Pj2R01za8CEhVFUZQrSsUoXIKxsXHq6zsJhcyYTCHKy1N5+unnWbz4kzidCwiFgrhc
   HaxbV4zZbCYQCODxeMjKymJ8fJy9e5sYGBhifDxCWpqDkRE3RUX5OJ3mWZkax8fHqa/vJhR6d+fE
   i8niqCiKolyfLmeMguooXCRN09i58wQ+n4VwOEB6ehlebxMu1wjJycuwWh2xp/8O1q6twmAw4Ha7
   eeaZZ3jwwQcxGo3xIEO9Xk84HI7/O1fQoQpIVBRFUS6WSjN9FQWDQUIhE5OTvfh8E+Tn1+B221i0
   KJdXXtlNMGjHaJwkIWGYmppcdDo9R450k5d3E7t2NcVHBS600VcpmxVFUZSrSY0oXKTpEQWzuRSL
   xUYwGMDrbWJ4uJekpCqEsNPS0k8g0MHSpXlMTYVIT78Bi8WKx+PG7T7BbbctxWKxXLZrUhRFUZSZ
   1PLIq2h6MyOPp4nTp/fQ33+AkpJk9HoHFoud7m43SUmVOBxVhMPZnDzpwmAw4nKN09TUwdGjE7zx
   xlHGx8fPee0dO3Zc+Tf0EaLq7/1Tdfj+qTp8/1QdXltUR+ESOJ1O9HoXJ0820NTk57nnDqNpUwhh
   RtOMCCExGKZISclEyggu1xgtLd1AAampxSQkVM2ZpEn9crw/qv7eP1WH75+qw/dP1eG1RU1+X4JA
   IMDRoy7S0v4bDocTt3ucU6e2kZzcitc7QiQyxJIlZYCkosJJINDK6KiH1FQrixZlkZDgYHhYJWlS
   FEVRrn2qlboEk5OTBIN2srOjSxUdDiejo5nU1hawdGkeTU2DBIMDRCJTrF1bjdVqRa8/SkJCNgkJ
   DpWkSVEURfnQuC6CGa/2NSiKoijKlab2UVAURVEU5QOnghkVRVEURZmX6igoiqIoijIv1VFQFEVR
   FGVeqqNwBQkhOoUQx4QQR4QQb8fKnEKI14QQLUKIV4UQSTPO/6kQok0IcVQIsfTqXfnVI4T4pRBi
   UAjRMKPsoutMCPGgEKI19j1fuNLv42qapw6/K4ToFULUxz5un/G1f4jVYbMQYuOM8tuFECdj9fj3
   V/p9XC1CiDwhxJtCiCYhxHEhxF/HytV9eIHmqMO/ipWr+/ACCSHMQoiDsfbjuBDiu7HyQiHEgdg9
   9bQQwhArNwkhfhurw/1CiIIZrzVn3c5LSqk+rtAH0AE4zyr7IfB3seO/B34QO/4E8IfY8UrgwNW+
   /qtUZzcDS4GGS60zwAmcApKA5Onjq/3ernIdfhf45hznVgBHiC6dLgTaAUH0oaIdWAAYgaPAoqv9
   3q5Q/WUBS2PHCUALsEjdh5elDtV9eHH1aIv9qwcOxO6vZ4DPxsofA/48dvxV4F9jx1uA38aOK+eq
   2/P9XDWicGVN3+gz3QU8ETt+Ivb5dPl/AkgpDwJJQojMK3GR1xIp5R7g7P2uL7bOPg68JqV0SSkn
   gNeA27lOzFOHEL0fz3YX0T8ompSyE2gDVsQ+2qSUXVLKEPBb3q33jzQp5YCU8mjs2AM0A3mo+/CC
   zVOHubEvq/vwAkkpfbFDM9GGXgLrgRdi5U8Ad8eOZ96fzwO3xo7vZO66nZfqKFxZEnhVCPGOEOLL
   sbJMKeUgRH+ZgIxYeS7QM+N7z/DuL9b1LuMC66w3Vqbqcm5fiw2N/8eMYfP56mq+ur2uCCEKiY7O
   HODCf3fVfTjDjDo8GCtS9+EFEkLohBBHgAHgdaKjUhNSykjslJn1Ea8rKWUYcAkhUriE+1B1FK6s
   1VLKG4A7iP5yrCHaeZjLXL1stenF+Z1dZ4Jonam6PNe/AiVSyqVE/+j8OFY+X11d93UohEgg+mT2
   N7Gn4gv93VX3Ycwcdajuw4sgpYxIKZcRHdFaQXSK5pzTYv9etjpUHYUrKPbUgZRyGPgd0f/owekp
   BSFEFjAUO70XyJ/x7XlA35W72mvaxdZZL1AwR/l1S0o5LGMTlsAveHfoUdXhHGIBYs8DT0opfx8r
   VvfhRZirDtV9eGmklJPATmAVkCyEmG7LZ9ZHvA6FEHqi8TDjXELbojoKV4gQwhbrTSOEsAMbgePA
   NuCLsdO+CEz/EdoGfCF2/iqiw0uDV/CSryWC2b3gi62zV4GPCSGShBBO4GOxsuvJrDqMNWzT7gEa
   Y8fbgHtjEdNFQCnwNvAOUCqEWCCEMAH3xs69XvwKaJJSPjqjTN2HF+ecOlT34YUTQqRNT80IIazA
   bUAT8Bbw2dhpDzL7PnwwdvxZ4M0Z5XPV7fyudhTn9fIBFBGN0D1CtIPw7Vh5CvAnolHArwPJM77n
   50QjUo8BtVf7PVylenuKaG83CHQDDxGNHr+oOiP6h7wNaAW+cLXf1zVQh/8JNMTuyd8RnW+fPv8f
   YnXYDGycUX57rM7bpu/f6+EDuAkIz/j9rY/VxUX/7l6v9+F56lDdhxdeh4tj9XY0VmffiZUXEY33
   aCW6AsIYKzcDz8bq6QBQ+F51O9+HyvWgKIqiKMq81NSDoiiKoijzUh0FRVEURVHmpToKiqIoiqLM
   S3UUFEVRFEWZl+ooKIqiKIoyL9VRUBRFURRlXqqjoCjXkVia44gQovhqX4uiKB8OqqOgKNcftXmK
   oigXTHUUFEVRFEWZl+ooKMp1TAixQwixWwixQQhxWAjhFUIcF0LcNce5NUKIrUKIESGETwhxUgjx
   92ed841YeVAI0SeE+JkQwnHWOREhxD8KIb4phOgUQniEEC/F9rJPF0I8I4SYEEJ0CyH+bo7rKBRC
   /EYIMSSECAghjggh7r78taMoCqiOgqJc7yRQAvwE+BHwKaAfeG5mHIMQYgWwj+i+8n9DNFX6j4lm
   nps+559iZa8Cm4EfEs1t8NIcP/d+YD3wVeDrwBrgSWAr0fwI9wB/AH4ghLh9xs/II5rAZnHsOj4J
   HAZeEEJsfj8VoSjK3FSuB0W5jgghHiSaxW+hlLJDCPEWsBqokFJ2xM5JJ9pZ+H+klD+Ile0CFgBl
   UsrgHK/rJJp46ikp5cMzyu8j2gG4U0r5UqwsQjSBTaWUMhIr+zHwDaKJbv6/WJk+9povTb+mEOKX
   RDsh5VLKiRk/5zUgTUpZe5mqSlGUGDWioChK23QnAUBKOQwMAQUQT2m7GvivuToJMasAE/Cbs8p/
   C2jA2rPKX5/uJMScJDq68dqM6wgTzXCXP+O8jwMvA24hhD72YYh9X810KndFUS4fw9W+AEVRrrqx
   OcqCgCV27CT6UHHmPK+REvu3f2ahlDIshBid8fVp42d9PnWecsuMzzOALwAPznENESAV8JznOhVF
   uUiqo6AoynsZJ9oI557nnDFAAFlEc9wD8emDVGD0Ml3LKLAL+EHs552t7zL9HEVRYlRHQVGU85JS
   +oUQe4D7hRDfm2f64QDRUYh7gbdmlN8L6IGdl+lyXiE6zdF0nmkQRVEuI9VRUJTrz1xP4u/lb4Ed
   wIFY4GEvUAwslVL+tZRyXAjxf4BvCyF8ROMIKoF/BHZLKf9weS6d/wUcBHYLIX4OdBKdGqkGiqSU
   X75MP0dRlBjVUVCU68/ZS53mWvokZ5ZLKQ8JIW4Cvgf8FDADXcCvZ5zzHSHEEPAXRJc9jgL/F/gf
   53vti7leKWWPEOIG4BHg+0B67Oc0Ak9cxGsqinKB1PJIRVEURVHmpZZHKoqiKIoyL9VRUBRFURRl
   XqqjoCiKoijKvFRHQVEURVGUeamOgqIoiqIo81IdBUVRFEVR5qU6CoqiKIqizEt1FBRFURRFmdf/
   D8IZtGFcztllAAAAAElFTkSuQmCC
   "
   >
   </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="s1">&#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="s1">&#39;black&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#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="s1">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#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="s1">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#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="s1">&#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="s1">&#39;red&#39;</span><span class="p">,</span> <span class="n">label</span><span class="o">=</span><span class="s1">&#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="s1">&#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="s1">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#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="s1">&#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="s1">&#39;dotted&#39;</span><span class="p">,</span> <span class="n">color</span><span class="o">=</span><span class="s1">&#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="s1">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="s1">&#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 0xa9f6adb0&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>