Crystal_structure.ipynb

{
 "metadata": {
  "name": ""
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Setup of crystal structure, energy and stress calculations\n",
      "==========================================================\n",
      "\n",
      "*This example uses a local installation of Quantum Espresso. Without such installation you will not be able to execute this tutorial.     However, you can still read it and lern from it. The \"Remote execution\" tutorial contains similar material and is configured for remote execution of Quantum Espresso.*\n",
      "\n",
      "To use the `qeutil` system you need to do three things first.\n",
      "\n",
      " * Read a basic introduction to the iPython interpreter and the notebook environment (The links to the documentation are at the top of the page under menu \"Help\"). Note particularly \"Keyboard shortcuts\" document.\n",
      " * Import a set of libraries into your notebook environment (see below).\n",
      " * Configure the qeutil for your particular computing setup\n",
      " \n",
      "The `qeutil` library is built as an extension to the [ASE](https://wiki.fysik.dtu.dk/ase/index.html) library. Many of the presented functions depend on the functionality of the `ASE` library. Thus, it is highly recommended to familiarize yourself with its structure and function. The full documentation is available on the web.\n",
      " \n",
      "During working out the examples pay attention to the comments in the code, which start with the \"#\" character and continue to the end of the line. They are added for your convenience - to provide detailed descriptions of the procedure. \n",
      "\n",
      "The whole notebook system is written in the python language. You do not need high-level knowledge of the languade to understand the tutorials. However, some familiarity with its syntax is still required. Introduction to the python language may be found in the Help menu, if you need it. Please note that the indentation in the python code *is important*. The blocks of code are defined by the common indentation. \n",
      " \n",
      "Importing\n",
      "---------\n",
      "This notebook is presented as a read-only document as linked from the main site. In the top-right corner of the page there is a \"Download\" link. To work on this exercise download the file from the link into your notebook directory and open it as the iPython notebook.\n",
      "\n",
      "Configuring\n",
      "-----------\n",
      "This part is a little bit more complicated. This notebook is configured for a local execution of the Quantum Espresso programs. The qe-util has some support for the remote execution with a help of the queue management system. \n",
      "This topic is covered by the second notebook named \"Remote calculation\". \n",
      "\n",
      "For this exercise just review the configuration in the second cell and verify that the commands actually execute the Quantum-Espresso programs.\n",
      "\n",
      "Additional libraries\n",
      "--------------------\n",
      "The following notebook includes a minimal set of libraries required for it to run (see the first cell of this notebook). You may need additional libraries if you extend the analysis presented below. It is as easy as adding additional 'import' clauses at the beginning. \n"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Import the basic libraries\n",
      "\n",
      "# ASE system\n",
      "import ase\n",
      "from ase import Atom, Atoms\n",
      "from ase import io\n",
      "from ase.lattice.spacegroup import crystal\n",
      "\n",
      "# Spacegroup/symmetry library\n",
      "from pyspglib import spglib\n",
      "\n",
      "# The qe-util package\n",
      "from qeutil import QuantumEspresso\n",
      "\n",
      "# iPython utility function\n",
      "from IPython.core.display import Image"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 1
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Configuration\n",
      "-------------\n",
      "\n",
      "Configure `qe-util` library for local execution of the Quantum Espresso on four processors. The commands should execute the following programs from the Quantum Espresso package:\n",
      "\n",
      "* `pw.x` the basic program for calculating electronic structure and basic static properties of the crystal such as:\n",
      "    * Electronic structure\n",
      "    * Energy\n",
      "    * Stress tensor\n",
      "    * Charge distribution\n",
      "* `ph.x` the program for calculating the a second derivative of the energy with respect to atomic displacements or unit cell deformations.\n",
      "* `matdyn.x` the program for processing dynamical matrix of the crystal.\n",
      "* `q2r.x` the program for the transformation of the dynamical matrix into a matrix of force constants in real space.\n",
      "\n",
      "The commands should execute the `pw.x` program of the Quantum Espresso suite using `pw.in` and `pw.out` files as input and output respectively. The example below is a fairly standard command for running the `pw.x` on the four-core PC."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Configure qe-util for local execution of the Quantum Espresso on four processors\n",
      "QuantumEspresso.pw_cmd='mpiexec -n 4 pw.x < pw.in > pw.out'"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 2
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Crystal definition\n",
      "------------------\n",
      "\n",
      "Here we define a $\\beta$-SiC crystal: a cubic zincblende crystal with a spacegroup F-43m (space group number 216) and an experimental lattice constant (here, A=4.3596 A). The atomic positions are specified in the fractional (crystallographic) coordinates; i.e. coordinates measured in units of lattice constants in the coordinate system defined by lattice vectors. The unit cell is specified as three lengths (in angstrom) of the lattice vectors and three angles (in degrees) between these vectors. The crystal axes are oriented in the conventional way in the Cartesian (X,Y,Z) coordinate system:\n",
      "\n",
      "* Vector $\\vec{A}$ is along the X axis\n",
      "* Vector $\\vec{B}$ is in the XY plane"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "a=4.3596                                         # Lattice constant in Angstrom\n",
      "cryst = crystal(['Si', 'C'],                     # Atoms in the crystal\n",
      "                [(0, 0, 0), (0.25, 0.25, 0.25)], # Atomic positions (fractional coordinates)\n",
      "                spacegroup=216,                  # International number of the spacegroup of the crystal\n",
      "                cellpar=[a, a, a, 90, 90, 90])   # Unit cell (a, b, c, alpha, beta, gamma) in Angstrom, Degrees"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 3
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "We can display the picture of the crystal. The `ase.io.write` procedure, used here, is a very flexible tool. It can write the crystal in numerous formats, it can even run external tools for rendering a 3D scene. Here, we just use a simple renderer build into the write procedure to store a picture of the crystal into a disk file and then to display it."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Write the image to disk file\n",
      "ase.io.write('crystal.png',       # The file where the picture get stored\n",
      "             cryst,               # The object holding the crystal definition\n",
      "             format='png',        # Format of the file\n",
      "             show_unit_cell=2,    # Draw the unit cell boundaries\n",
      "             rotation='115y,15x', # Rotate the scene by 115deg around Y axis and 15deg around X axis\n",
      "             scale=30)            # Scale of the picture\n",
      "\n",
      "# Display the image\n",
      "Image(filename='crystal.png')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "png": "iVBORw0KGgoAAAANSUhEUgAAALwAAADWCAYAAABrJuiXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXd4FFX3x7+zO2VnN4UQEgIkpFAMoROBUAURkSIgxQZi\nAaniq6IoiBLfV0VsIF1AiiIgiFjAAAqCIF0IPdQ0SighdXdmZ3dnfn9M9semkLq7M5vs53l4zJjJ\n3LPJd++ee+655wBevHjx4sWLFy8ej1ZpA7wohgZAXQAhAHIBiMqa48VL1aAA1C/h/38G4CIAHsAt\nAAkAbgMY4D7TvHipOD6QZ2lHngawD0A6ADOA4yX8XAsA0QD0BdcEAAnyG8GLF0UgUNzdbAVgM4Bj\nADIBmAA0LHJPUwDdAYQDICswnuTwdXsA3wIIq8DPe/FSIQwAvgawDcA5AEYATxW5JwTAMMiCDIb8\npnAWjoL3BfA/AIko/gnixUu5eAPyrLkbQDKKuw9aABMB9Ifscvi50zgUFrwd+yeEL4CxkNcEHo8z\nZ4maSl8AD0N2I8IB3AQwsMg9LwOwAUgt+Gf3sdVCfMG/kmgIYAWAegAeBMC5xyTX4BV86TSDHL2w\ni7kBZJfC5nDPcACNcE/MyQBuuNdMl0MAaAngZMF1UwAXlDOn8tRkwdcCMAj3xNwQwOsATjvc0xXA\nE7gn5lTIf/SaHLMOgRz9+RPAZADZyppTMaqr4AkAj+GemMMBbADws8M9oQBmQRZxWsF/D0DehPFS\nOna/fh4AC+TFrUdMAp4q+BjIsWS7mC8DmF/knk0A7uLezPwPgCtutLGmoIc8438FYBnkN4BqUaPg\n66CwmBkAM4vcMx+yC2IX81EAe9xoY3UjHvdftJaH1pA/LSdBXsOoFncLngYQgcJ+82zIcWc7r0GO\nQdvFfAHAN261suYhwXlaIADMBfAd5ImoWuMPeUfwcQCvQF7NO/IoZPdjF4CVkGcVfzfa56VkSorD\nVxYNgDEArgPo4sTnuh0N5FhsHOQZOK7I92sDyIMc5dgKYBFk8XtRP84UvB3HXJ/H4GETmwZAVsG/\nIwB+hLzN7QgBda4JvJRNvIufvxDyhlxPF4/jNAjIUY9kePjHlBfFaAU5/x4AguABk6MEOWJiN3QY\n5CQmL14qygbIUbU27h64otlwZtzz9zoAOAvgEada5KUm8AzkZLpApQ0pjZIWNo0hL1YBOX/am07q\npTJ8Djk8XUtpQxyJL+P7ayAfTigavfGifuIVHr8+gOUqsKNCEJAzBzs7XHvxDFwRlqwMds30h5y0\n51EaWgBgMe6tzL2oF7UI3s5DkLNS5yltSEWoDeBLAG8rbYiXMlGb4AH5FFhIwdehACIVtKVSDAfw\nPLwLWzWiRsE7MgTywfX33TlofBV/vgOA/ZBTDryoi3ilDSgH9SAfpwTk87a0qwd0xixAAAgo+DoK\ncl67Fy8VZSDkYlJD4MKFrbM/9gZCrnw1y8nPJSGXvaDhYat8LxXiUchpyC7DFX5ebciFgwA5bUFf\nyr1FYSDH/F/R6XQ/sCybQhCEFYCo1WotBEHYAIgMw2SwLLsVwFQAveD+EhheXE9LAOtQjoVtRWZA\nZx4SKImBkLPqZkDedr7fGyyKoqjJAEb7+/tLUVFRVHh4OBsWFoaQkBDQ9D3Xzmaz4fbt20hPT0dq\naqqQnJzM3bx5kyFJ8hee57+AnPnppeonnpTGAGAK5FKDrQBY73djRQQcD9f/UjpBzp0ueqQPAHqy\nLPuBzWZr36lTJ03Xrl3poKCgCg+Qn5+PgwcP2nbv3s1brdarPM9/COB7qD9S4UpcPZm5CxKy2FnI\nEcEVAATHG9T8ImMB/BfAfxmGmURR1NDHH39c37Zt20KzeGURRRFJSUn49ddfjdnZ2Sd5nn8G8pHC\nmkh1ETwD+dhoW8ieQj3Ibk6+/QY1v0gawFyCIF6OjY21Dh06VMeyrNMHsdls2LVrl3XHjh1mm832\npiiKX6PmzfaeInh/yIK2n4leDzk+b2cigDdxr/SKEXJgJN5+g1pfJEHT9CyapiePGDFC36xZM/A8\nj3379qFbt25gGMbpA964cQOrV682Zmdn7+N5fjDk+uk1BTUInoB8vsIu5pMAzjt8vyPks9CORbE+\nBHC1jOcWem1q8+EBQEPT9NcBAQHPTJ482eDj4wMAyMvLw08//YSUlBS88847LhG91WrF6tWruQsX\nLpwwm8294fBRWM2Jh+v/thTkUoXhkH+v/zp8j4RcwYzHPTHPg1x81o4Gsngr+ulbacG7YxYgaJpe\nEBQU9Pwrr7xiKMmFyczMRGCgfG4gPT0dYWHOLWMuiiLWrl3Lnzp16rjZbO4JdRU9VTN6yGL2BXC4\nyPe2AOgN+VxrKuQTT0ULZ+kh17x3NuoVPEVR79SqVWvG66+/btDrSw/JG41GzJkzB3Xq1MHIkSNh\n/yRwBqIoYsWKFdylS5f+KHBvappPD8j+cjsAsSzLPgQgRhRFAwBoNBoLQRA8gCSTyaSBfM6ZgVwV\n+SCAUUWeFQR5BleiKplqBd+KpumD06ZNYwMCAsq+G7ILcvToUbRv3x5arRZWqxUkWZHGF6U/e9as\nWca7d++OlSRprVMeql60kCMaZgCDdDrdWxaLJdLf31/gOE5vsVgIrVYLPz8/PP/889BqtRAEATdv\n3sTFixeFtLQ0/vbt2yxN0yk8z38GYC0KF9dSknhUctHqSsFTDMOcGjx4cNNOnTpVaowCgaJ9+/bo\n2bOnU3z8tLQ0zJ8/P89isTQFkFHlByoHAznqlefw/0Ig707ay4BbtVqtpnHjxtYePXr4NG3aFBzH\nITc3FwEBASgrQiaKIs6fP4/du3fnX758WaPRaL4VBOFLyDkvqqEibSsJFF5EOA2SJN9v2LDhY0OG\nDKEJonLvKY1GgxYtWuDEiRNo2LAhynKJyoO/vz/MZrPmxo0bbaxW6/dVfqDroFG4Zj0gH5mbDrl9\nzf8gV0Xe7/B9C4B0jUYjkiTZsWfPnpoRI0Yw3bt3p4OCgqDRaEDTNHx9fUFRZTf/IAgCQUFBaN++\nPd2+fXtKo9G0vnr16jiCIEhRFPdDJdWFlQ5FAUAoRVEXpk+fXm5XpjxIkoT169ejbdu2iI6OrvRz\nrFYrPvroI2NWVtZwyC0e3Y29y54jkwD0wb36nPsglzd0ZADkwlmpkBs0FH1DtGAYZmODBg3Cnnvu\nOYMzf/d2srKysGbNGuPVq1evms3m4QBOOX2QCqJ4Y2KSJKd26NCh84MPPugc59sBSZKwefNm1K1b\nF3Xq1KnUMzQaDViWpS9fvhxqtVpXOdfCEukB4F3IzQZmQJ6lvyhyTy3IJVJWQ569V5TwnAuQF5G5\nKPKG0Wq1EymK2jh48OCQoUOH0q7Y0AMAlmXRoUMH2sfHJ/DixYsvEgSRJ4riIZcMdn/i4eCZKD3D\nUxRF3XrjjTdq1atXzyUDWK1WaDQaaDQanDhxAhEREfD3r1iZQ0EQMGPGDF4QhBhUrRx0OIBxKNyo\n4WEU9nM7QN4at8ej7TuGzoCgKOp9vV4/dfLkyfrKTgKVITMzE/PmzTOZTKYvLRbL+3Bf5KvQ2lPp\n43aDQ0JCtK4SOwCQJAmNRn6Z165dw+zZs3HmzJkKPYOmacTFxWkoinqltNsgl/r+EnJZwqOQDygU\nhQOwHfIs3g3FmzQcRvEWlk6Boqj3fH19p77xxhtuFTsABAYGYsqUKXo/P7/XKYoqKTnQLSi608qy\n7L5hw4Z1iY2NdeZjSyUzMxMMw8DHxwdZWVnw9/f//zdEafzzzz/YtGmTWRTFlZBn5r0ofHhFA2AO\nZDfCPjOfg0pa6JAkOUGv138+ZcoUfUU/4ZxJbm4uvvjiC1N+fv5Um8220A1DqiYOT5AkaZw5cybr\n6+vrxMeWn/Xr1yM1NRU9evSAJEnIysrC3bt3QdM0hg8fXuje33//Hbt27bJYrdaPIW+Ln4bKu104\nEE3T9LG33nqLrUxKtbO5c+cOPv30U04QhAchr0VciWpcmkYMw4hKiR0AnnrqKfTt2xf79+/H5cuX\nIUkSGjVqhHbt2hW7t1+/fmjUqBEHedb+DZ4jdi3DMBsef/xxRg1iB4A6depg4MCBDMMwG3CvAbKr\n+MDxQknBx4aGhhYNlbkVgiDQqlUrpKamYsSIEejXrx/i4uLQqFGjEu+PioryIUmyo5vNrBJarfbN\nkJCQqC5duii9XitE586dNfXq1YsgSfItFw8V73hRkV/CB2XfUn5IkuwYFRXlvAQYNxAWFqahabp7\n2XeqhoYajWbmqFGjDOVZp7gTjUaD5557zqDRaN6D3PfLPeNW4N54Zw5M03SjwMBAdf0VyqB27doQ\nRdF1ISUnQ1HUpLi4OI09u1RtBAYGomPHjtoyol9ORUnB6cuzZa0mKIqCJEnOT8R3DQyA8V27dlW1\nvd26daMhNznWuWM8xQRPEARZ2bwZZ9OnT59y3afVaiFJkqd8Kg1t0KABUbeuumvZBgcHIzQ0FCje\nM8xZxDteKPbHkySJs1rvW03BrfTt27fsmwBYLBYQBCGUfafysCw7pWfPnsqFwCpAz549fVmWfdNF\njy+0yaWYDy+KYjbHcc58pMvhOA4ajcYVp3KcDS0IQstmzZopbUe5iI6OhtlsjoHshrmUigjeqdvB\nPM8fTk9P96iD0tevXweARKXtKAct/P39OWeUM3EHNE0jICCAg1xBzKUo6Y/+m5yc7FHnRVNSUswc\nx/2ttB3l4MHw8HDFM2ErQkREBAngQRc8WjUbT8fv3Lmjt9kU3XsCACQklC/NPTk5mUfh0/aqhGGY\nLpGRkQal7agIERERep1O19UFj453vFBS8PkURd26ceOGgibIbN++vcx7rFYr7ty5o4cHuDQkSTYL\nDvasFrrBwcHQaDQuX3QottMKADabbcORI0c8Iupx8uRJUBSViMLnQlWJJEmsp/jvdgr2ZFwei1cs\nSgMAFotlwcGDB0VBUL/md+3alcfz/Gyl7SgntNpSCcqioNqEy9+lrs5UK4sMSZLO//DDD631OgaC\nYIbFYoHVagFBaEBRFCiKgt7HFw0ahCIsLAx16tSBszesytp4un79Om7dumUD8KtTB3YdZjWsjSpC\nwZ6MK4IY8XCYrN0t+EYA+vn76LtJkDqaeKF+w5DaQtO6OsTGRMHPh4WOocHQJGw2CWbBAl6w4OrN\nTBw9cxo7ErYgz8ghvGEo6jUIQ8OG4WjWrFmVqwmXtfG0Z88eXhTFBVCmkFBlMJnNHhUAQ4G9rtiY\nmYlKCj4elXNrtAD6+/vop9pEsd3Ang+iY8vGbOvoCMQ0CgVNVaxy0p2sXCQmpeB4Ugp2Hz2ETT9u\nRIcOHRDXqXOlD2qXRlpaGo4dO2axWq1FS8OpFqvVeuLmzZsdY2I8p4VWRkYGrFbrCVeP48oTT3Vp\nihxLktpXI+oHMa+O7Os7+OH20DHOddOSr97Csk27sHbrPoSFhaFjpy6IiYkp17G9svDg6mPPt2zZ\ncsHo0aM9Jv16xYoV+SdPnnwVwEonP7qQbl3h0ugNLPO5TRRffKJXB2nck73ZNtERLhhGJjI0GB//\n52m8N24INu88jCUbduLXXzbjiSHDUNWt9YSEBIvRaPxHkqR1TjLXXfyblpbmUfUwU1NTJbhhj8PZ\ngu9mYJn1j8S1DPjy7ed1gbXcl7vE6mg8278rnu3fFbsOncakj1bi1MkmeHzgoDLLxCUkJBTz4y9e\nvIi///6bt1gsz8Pziqkm5efnMxzHlfna1QDP88jLy9PBNedbXbLTqjewzGJ/X/32r+PH1l896xXW\nnWIvysMdW+DQug/RKIjBF599inPnzpV6f9GNp7S0NCxbtsxksVgGwTNrSlppmj546pTihb7KxalT\np8AwzGGU0oysCsQ7Xjhj46mrgWUuPhLX8vnjP37KDnjIfSU3SsPPwGL+9BexPH4Mft38Izb8sB48\nX3auWkpKChYuXMgJgvAsgL9cb6lr4Djus127dql+kwyQ9zg4jvvUHWNVNEpTCIIgnjGwzPIlM1/W\nq0XoRenRoTkOrfsQU7/4HosXLsCYseNwv0oJJ0+exJo1a0yCIAwH8Lt7LXU6CZmZmcK+ffsgSRLM\nZnmPw2azQqPRgqIokCSJWrVqISwsDIGBgU7f3ygP6enpyMzMNAPY6o7xKp1RR5HaCX4+7IKEJdP1\nXdtVvlipO2BoCv27t0VmVg5WrP0JzZu3KOTbbtu2DRkZGdzOnTvvCILQH545sxsAdCeA4QF+humk\nVvu5KIq1IeRpGgYyqOunRVgAjYggA0L8KPjRNuhgRlrqFSRs/wPbd2zH5YsXcOPGDRhNJuh0Orf4\n/1u2bOGvX7/+uSRJu100RDyqWluSpsjJtXz1n2xfOkMfGepZSUoL1m7DvLV/YMKkVxAQEIDTp09j\n1apVFo1Gs0IQhClQTyH/8tLMwDKv2WziyMbhIZZOrZuysc2j6DbREWjSsB602vJ5rbcyc3DifCqO\nnUvG0bPJOHLqEho0aICOnbqgefPm0Gqdn2187do1zJ0712ixWKIgd9tzBZWuPAYA0Go0z/n76r/+\na+VMNry+Ogr7VJSv1iRg3trt0Gip/JycnJwCf90T8tztUAAG+/vopwJoPnrow9RLQ3qSoXWdV53A\nLFjwy64jWLJhJ9IyMtExrjM6xsXBz8/PKc+32Wz45JNPjHfu3HlVkqSSqh87iyp18Tvsa2B//GP5\nDDY6soGzDXMrMxdskJZv2nnNyJmbwXO69Wlpipys1WpmRkfU104e0dd3QI9Y0JRrM0ROXkjF0h93\n4eedR9A8phn6DRhY4QrMRUlISLDs3r17f0HjOFeGfStfW5LV0Tk/zX3Tv1Obps43y81IkoRR0xZw\nuw6dXmXkzBOVtqccNPE16DY0blivyfzpLxlaNHFu98LykJNvwtzvfsfKzXswYOAgxMbGVmqhW9BK\nKN9isTRD2X1Wq0rlBM8ytPTC4B7CrNef9axE61K4m5OPdsOmmrLzTD8AeElpe+6DliK1r5Gk9n/v\njR9Kj3+yt1bp1N/EpBS8HL8MBr8APDF0eIVm+5s3b2LOnDkcz/PPAvjZdVb+P/GoRFOzYfWCAjYe\n2zgbrK7a6B0AsPXvYxgxdR4gRznUVpGgia9et6FxeEiTb/47wRAVpp4aM2bBgtnf/IoVm3eXe7a/\nefMm5s2bx5lMpkmSJDk7Z6ZclEfwQayOvvjzvLf8O7Zq4nKDlKBW3AswsMxilbk2j7AM/fN7E4bq\n1DCr34/EpBSMmbkU9cIi8cSQofdN2ktJScGSJUs4s9k8QZKk1W428/8p87foa9B988KgHmx1Fbsd\nitQ+D7kjhxp4wsAyv276aoph4tN9VCt2AGgTHYFdK96DkHcH69d+j6IHT6xWK7Zs2WJZuHBhHs/z\nTyspdqBswXdhGfrh9ycMq15+TAl8Ne1Fva+BXQ6F+14RBPG0nw/7/dbF09jObR5Q0pRy42dg8fNX\nU+BHWfHdt6v+X/Tp6emYNWuWce/evbsLet0qfmKs1N0EPx923lsvDmrRqU1TdRSBdBGSJOGlIT2x\ndOOfOiNn3gngmkKmDPA1sGu2L32Xbdm0oUImVA6K1OKJR9pj61+Hse/AEZw4edK4detWs9FoHG+1\nWt+GcqHfeJRzpzVYx1CpSVvm6mr5elSJk0oz//sE8bMVv/6Ua+SGl3230+luYJmEXxe8rY9tHqXA\n8M6B4wX0nzBLOnfl2mHOLPSF3CtWScrX8oamyLFP9Oog1RSxA8CIAd00gsU6AIB7W9wBAXodvfnb\nWa94tNgB+VzCzwumEkG1/VoAeFRpe4pyP8FrSVL7n3FP9lb/6QEnUtvfBwN6xEqkVjPGneP66nVf\nP9W3i6FXnMtLK7oFPwOLlR9ONLAMvRSAemKpuL/g+0fUD2JceTRPrUx8+lGWoanX4b4u5Y/rWab/\n/yY/perGBRUltnkUXh7WS+er162CsoGAsk88+fvqx4x/srdH1BZ3BrOWbf7/r9vFRCGoth8LwB3N\nywL0Onr18v+O1/vo3dIAw61MH/sEHeDv0w3AkwqaEe94UaLgbTaxQ3WPuzsy+5tfCl13a9eMgmsq\n2RaiwJXRd4v1jDruFUXH0KpzbUoSfIBgsdZu3DDE7caohfYtGun8ffQPuXiYbjqGrnauTFFim0dh\n9NCHdb4Gdq7StgAlC77dA5H1TeU9OFAdaR0dAQlSe1eO4e+jf+etlway1dGVKcobzw+gLRbrYLg/\n+lWMYqomgNhOrZvWqOjM26MHFbpuFtUAHC+EQE4ocwX1zRZLr6f7danWG3p2avv7oP9D7URSqxmt\nwPDxjhfFBF/Lz/BQbExUtU8lcGTay08UuqYpEpGhwSYAbVwxHkOTE4b36ST5GWrOvDLx6T56hqbe\ngPt7EpTe1MxqtbVrHR3uPnNUSseWTWgA7VzwaEqj0UwaN7x39fdlHGgXE4kGwbVZAOXrEeoiigle\nsFr96tau2vGt6kD94AAWgCtaWA9sEl6PVOLEkpIQBIHJI/v6FpzDVYxigrfZRIphPKtDtivQMRRo\ninS6D1/LV//qpKf71Jg9DkeG9u4IwWKNg4IhyqKCJ6w2kdTRNUvwjhtPdnQMDZLUOrv6LsGZLbHV\nNe5eFnodg5ZNG5oBuDQCVoQyd1prROTAkaIbTwBAgADh/PSChgxNauoF1XLyYz2Hzm2aGrRajTsF\nH+94UVTwklajsZoFdbSEVxLeLMBqtTk7hzu2VdNwUYmSdmqhbbNI0g2bevel2AxPajUWs+ApnV1c\nBy9YYLZYnVqMlCK1HeJaN6k5+dYl0CY6ArxZcEm4tzwUFzypNd7N8ZS6RK7jTlauACDXmc80sLr+\n7Zp5dr57VWlYrw4AQgegnpuGjHe8KCZ4miJPn7qY5iZb1EHRnVYAOHzqEgfnNiEmeEGIrokp144Q\nBIEWTcJ4AO4qN136xlOukdtz7GxyjXLii+602mwizqdcZwEcc+Iw9cyClazJC1Y7UaF1KQD1lRi7\npDj8kf2J5z2tgq5TuZSWAYoks+Hc85jP0xSpSA12taFnaS3c0HW7JEqqwvnvmUvpOkmSauwfJzEp\nBRSpdebsDp1O15uh3HWISt0YWEYLQJFEopLi8BkAuNQbd9xti2IU3Xj69+wVITvXuNuZY4ii2MoV\nNdY9EYokNZBLfruDso/4MTR1/N8zl91jjgoouvH0z7EkXnJuC8VAm83ma7HWqKXRfTHxZhtc03W7\nJOIdL0oUfFaucfXKzbs9oiGWs7mUloErV28RAPY68bHt6taty3n3N2RMnNkK9wm+EPfLTd5w9Mxl\nIvmqq7qQqJdlG/80S5CWAjA78bENgoKCtDabDV7RA9l5JhsUai10v9YRHEEQy5dv2jnxo/88U2MO\ng5h4M9Zs2SvxZsv8Um6jADSHHEeur9FoDARBkDabzQi5nNwpyO7QbYef0el0Om2DeiE4d+Uaanos\n/sT5VBHAGSXGvm+vFI4XFqz6eff4GeOGVrua8EWxbzz9uOMgSK32AIDUIre0pml6PEmSPXmej/Lz\n8+MjIiKIoKAgA0VRhEajgdVqhclksqSkpJhu3LihIwjCSJLkMZPJtBYAq9FoiAahoUhMSqnRgud4\nAddu3mUhTwzuIB7lbYjg76v/e9Zrz3Z7tn9XF9ukPJIkof1T0/IvpWUMB7ANAANgGMuybxME0ahb\nt270Aw88QDZo0AAMU3qhAUmScOfOHaSlpeHgwYP5ycnJVJ06dbStWrUig3VmLJj+ojtekio5cvoS\nhr3+5eWcPFNjNw1ZqLZkqd2wcvJMn365ekvbpx7r7FPdqxjsT7yAG7ezjAB2AHiKoqivQ0NDNT17\n9vStaNtGgiAQFBSEoKAgxMbG+mRmZmLfvn3Yu3cvQoNr9k5rYlIKJFE6oNT4ZbV/+/1mZvaFxT/s\naP3Ks49V2yAybxYw/oOlRiNnnsEwzBaDwdB91KhRhoiICKc8PzAwEIMGDUKvXr3wwQfxMAsWMDXs\nkI2dw6cumXKN3D6lxi9r2hbzjPxTHy39yXwpLcMtBinB/5ZsEm7dzUmiKOqLzp0795o2bZrTxO6I\nj48PQoKDcPayqxvXqZfDpy5b4dw9jrIoe+OpCJcsVtuMF99dCJtNdJFNynHk9CUs/mGHRkvSMZMm\nTfIbNGgQTVGum32bPhCNjTsOuez5aubclWu4dTdHBHDCjcPGO16UyzG3Wm1fXU6/hcU/7LCVfbfn\nwJsFPPnGXEEUJXLKlCmsK2b1onTq3AVrt+4DxwsuH0ttfL3hD14SpUUAFNuMKO9KVDTx5vkfLf3J\nfDH1hksNcidPv/mV1SLK51YDA11RkaM4derUQVhYGDbvPOyW8dRCnpHDDwn7wQuWxUraUZHQy6tW\nq+3NQZM/Nd3KzHGZQe5i1rLN4oGTl8hXX33V7Yvxjp26YMmGne4eVlE2bj8g0RS5B67vvF0qFYo1\nWqy2xVk5xrn9JswyZud5bsr85j8PS3PXbMPQoUOl4OBg9Onj3mJYMTExuHEnG8fPJbt1XKWQJAnz\nvk8w5uSbPlNg+HjHiwoH1zmzMOPG7azvBkz8xJid63miT9h7HOM+WGYJDw8X2rdvTwBA37593WqD\nRqNBx06dsfiHP906rlLsT7yAO1l5OQB2KTB86Uf8yoFk5MyTrqTfWtVr9H+NNzOznWSX6/khYb/0\nwvSF+YRGI40YMUKn5AGXTp0648+Dp3Eg8YJiNrgDwWLFKx99YzRy5ncg73oqSkUEH+/wtWjizZOv\n38r6sscLH5hOX0x3slnOxWYTMefbrbbXZq/KgkZ7tF+/frS/v7L1Mw0GAwYPGYrx/1sOE+/MxEx1\nMfubny137uYelCTpe6VtASom+JlFriXOLLx/KzNnYu8x/zN+svxnixoPOFxKy8DDL31g/HL1byc4\nXnhMFMW4Dh06qOLsYqtWrVC3Xig+WLxJaVNcQmJSChat22HOM/HPQQWzO+CEWt02UVzNmYVmC9dt\nO9B5xHtGtcz2NpuIBWu32bqNet907sq1GXlGvr1Wq+3Xtm1bsGzh45QJCQkKWQkMHDwEG7cfqnau\njWCx4sUiWrjnAAAW7klEQVQZi4y8YJkIQMlYdoV3WstDep6R73E5PeM/apjt7bP67G9+PsHxQmvB\nYp0LQKPVal/t3r17sdPy27dvV8BKGR8fn2rp2ji4MmsUNiXe8aIigv+gjO9Loih9w5mF6IXrtu9v\nMWiKcd6a30V3VTGTJAlHTl/CSzMWcV1HvsfZZ3UAlwpuaWMwGKjQ0FC32FMRWrVqhZD6DTHug+XV\nIn3j97+PY9G6HSY1uTJ2KrtoLY2reUau583MnF6fr/ztp2aPv24e8/4S7t8zVyBJzn/tJt6Mb3/Z\ng/ZPTcsbPPmzm7/+dTSeFyxhBbO6o3piIyMjVZvxOezJp3Dxejb+88kql/ye3MWeo2cx+v3FRs4s\nPAJlXZkSKSs9uLJIAA7lGrnhAOr8vPPwmIS9x19vEFybHTXoId+2zSLRqmlD+Faix5EkSUi9cQcn\nklKw5+hZYf3v/4gUSR7IyTd9CjmXvcQpkmGYLpGRkfoqvSoXQlEUXnhxNJYuWYR35q7DJ68943F1\ngQ6euIhn3vzKxPHCAABHlbanJNz5G9UC6GNgmSEURXbJN/GNgmv78Q82b0R0bNXEp0WTMPgZWOgY\nGixDwWoTwZkFmM0WXLt1F8fOJlsPnLhgPH0xXSdB4nQ0dTLXyO202cRvAaSUNTjLshfHjh3bODIy\nstj3EhIS3L75dD+MRiO+Wb4UnVuE46t3XoCnHLzZefAUnntnvsnEC0MAKLcoKk48ynvEz8VQAGIA\nxBpYpjNNke0lCQZREhlRlHQEAZtGo+EJgjATQEaukdtts4lHIOdSV/SjkiAIwjJ79mwtTav/fC7P\n81i9cgUigg1Y/t/xMLDq7l38446D0uSPV+RzvNAXwD9K21OEQkf8KiL4eJTfj1cbNEEQ3Jw5czxj\nugRgsViwcP48WMxGfD1zLLq2i1bapGJkZufhtU9WcbsOnb5r5Mz94d489/JSSPBV2XjyJHRardaj\ncvkpikLzlq3QrlkkXp65BFO/WAMjp56w5W+7/0W74W9zOw+d/sbImZtCnWIvhsfMeFXEKknSfT/N\nlNx4Kg1RFNE0oj72r/kQ2XlGdBk5A/uOJSlqU2Z2Hp57Zz43/oOl13LyTL1NnHkyAJOiRpVOoXC6\nq6I0aoO32Wzk/Soib9++XTWLVkcsFgtYhkWAvw+Wxo/D738fx8szl6BL22i8PKwXOrRs7LZIzs3M\nbKz+eY9twbptZqtNXGHizG9D3UK3E+94UVMEL1IUlZOVleVfu3ZtpW0pN3nZWagXdK8rer/ubdG5\nbVOs+W0vJvx3OfQsjTFDe2HYo3Hw0Tu/3LokSdifeAHzvk+Qdh8+Y6NI7YZ8Ez8bwEmnD+YmKiL4\nsnZaVQ1FUSfT09O7eZLg06+mo0304EL/r5avAa88+xgmPv0odh85i+WbdiJ+0UY82acTBvdqj1ZN\nw6skfptNxKW0DOw5elZauHabMTMnP8fEmZeKkqQ3C5bpkPc5CKhsB7W8VETw8a4ywh3wPL8nPT29\nc+vWrVW72+oIz/O4m5WNByJK7gyj0WjwcMcWeLhjC6RnZGL1z7sxY956nLt8DQ3r1UHr6HC0iY5A\n2+hINGvUAL56tlBMX5IkmAULUq/fQWJSCv49e0U4eOICl5R8XU9TZJZGQ+zPzefmA/gLhcUdAiAB\nchDjN3iY8KuzS0MAqA0gHEC4KIoZly9fzgdQLBHe3Uf8ysPVq1fxQFQoSLLs92dYSCBmjB+KGeOH\nwmK14tyVa0hMSsGJpFRs2nEIScnXYOIFkFoNdDQFS0EVYw2hkVgdfZsitceyc427C2riHxMs1rul\nDHcTwHTIgt8PQO2dM+Khko2nqqKF3PowHHJpa8etbALyH4aGXBg1FcA6iqJWfPjhh7qyakOqgW0J\nvyO8FjD79Wed8jxJkiBYrOAFCyhSi6Qr1zBg4uxME28OAVCV1FYtgI8ALAKgxvaPlY7DxzvdlNLR\nAWgCoGMJ31sHOUJwBMAXAHoV+b4EIApALQCtAQwEsI4kyQPHjx93mcHOwmaz4fChQxj1eDenPZMg\nCDA0BX8fPfQ6Bu1iohDZIIgG0L+Kj6YgTzjHATStsqEupiIzfKF3ihPwB9AQwDkUnmEmQ/7IrA3g\nGuRf5NAiPxsMuWkwX8ExHwsKCtowdepU37y8PFgsFlgsFthsNpAkCYqiQFEU/P39odEot0WRmJiI\n0//+gz+XvevScX5I2I+pX6zZn5Nv6uKExwVBdm8kAD0B7IOCBZccqHRqQUUET0AWpb1JgJ3OABZD\ndkNIyK7Gw5DdDzshBd+7AcAZu6M0gBYAYn30uq4ajWYEb7Zoa/n7gGVoMDQFitTKH/dmCzjeDCPH\no2FoA9QPDUX9BmEICwtDcHCw294EXy9eiDdH9sLQ3nEuHYc3C2jc91U+38S3xL1zA1VFC3kx2wjA\nMLivDvz9KH+57FKgINdPdxRzLQAbIYs5rOB7I1E4c+4cgBcg+3p3UfIK3xlVW3UAhtfy1U/JN5lj\n6gXV4tu3aKTp0LKJoU10BFo2bVgoIWvWss2FmhNn5xpx4nwqjicl4+iZC9izawcys3LRpHEjdIjr\njJiYGJeJPzk5Gbdv38LjPR50yfMd0TE0Rgzopl2+aecIm010VtjZBqAfgD6459P7AlCqZ1ih11XR\nGf4AZDHXBfAxCvv1WgC9Ic/aaVCmh08Uy9CvSJBebtcsEpOefcynR/vmZWYb1op7AdkHV5V6T3au\nEb/vPY4lG3YiIzMHHeI6oWPHOPj6+jrNeEEQMOfLz/HJq8Mx6OH2Tntuafyy6whe+2TVnqxcYw8X\nDrMPsi6mQeGFbXkFrwVwBfLM+SKAP6AO/8xOL39f/ftWq63DcwO7a14e2otu1DCk3D9cHsE7cvxc\nMpZu3Infdv+LmOYx6P5QTzjj6OBvv/4CA4z49uOJVX5WeUm9fhudnn0328QLAS4cxgfAm5AP6Ox3\n4ThlUtFFaBfIWXH5AOpA+RhsgK9e9zWrY/q/P3GYfsgjHaDXVTzkWFHB28nKycearfvw5eqteLB9\nB/R5rC9IsnJeYnJyMtZ8uxKH132EwFrO+9QoC0mS0KDneLOJN0fAOe5keZgP4AKAJXDzxFlRR/Qf\n3PPbNwP4FYC7evUUZYBeR18e1qfTwGM/ztaPHNCtUmKvCgH+Ppj87GM4tPZDaM1ZmPvl50hLq/gn\n9u3bt/Hdt6swf9qLbhU7IIcrWzQJ4yF3JXQXywAMADDODWPFO15UZZt9LeRIzA2497BugK9e922A\nn2H6d59M9hs7/BGSpqq2YSxJErrFNqv0zxv0Ogzt3QHBAT74bNFqGI1GREU1KtfC9u7du1j69WJM\nHzMQTz7WqdI2VIVLaRnk0TOXr0iStNtNQ94EsAbyzq4EeZ/EH3IY2tn8BYeFa1VCDWYAXwI4VnC9\nEsC7AFx5UDqO1dGXhvXpNPDfjbP1VRGpI44RmspCEASG9+l0b7af8wUyMzNL/ZmbN29i0YL5eG1E\nH7z0RM8q21BZ2sVEkn4GtrsCQ9vDznrIHsObrh7QmbG1DyHvao5y4jMdeYTV0TtXfjix9py3n2dc\nkQ7rDIID/bHu08l45amHsWjhfGRkFHeLJUnCoUOHsGjBfMRPHIKJTz+qgKX3CK7tD4Ig6ihownrI\nu7TrCq7DIIe5nY4zk8cuA3gS9xbCz0D+6HJGieRBBpZZt+HLN9gubR9wwuNcC0EQGP9kb9Ty1WPa\n3EV4aczLCAsLAwBkZ2dj048bYDHl4beFb6FV0/AynuZ6dAwFUZQqXjPFuRhxL5Q9CMB7AN6B7Dk4\nDVdkS9o3k4yQFyfLAcyqwvMeNbDMui2L3mHbNiteYkPNPN23C3z0Okz+eBlGvfAiUlJSsHvXTox7\nshfefOFxVHXt4Sx0DA1REtWUUbcAwB7IMz0g69SGyqUiu+2I36+Qd1n9Cq6bQM5/uXnfnyhOZ72O\n3vzjnCkeJ3Y7Ax6KRZ6Rx5RPF6NrbAy2LnobLZqElf2DbkRDEICkuvPNp3AvLWEM5F37KQAq2gIx\n3vHC1S/SDOB2wdcPAzgDYEI5fzaEZeiE1R+/ou/URvVJeKXyTL8umDlxONKu30LjhnWVNqcYnFmA\nRkNUNBHPnSwD8A1kN7lKuPNd/TXkVN9zBddMKeMTvgbdt2OHP8L27tzKLca5mrHDH0HjhvUw+5tf\nlDalGLzZAoJQteBtkH351wquOwP4DECFd4fd/TF2GcDugq9fhByHLZrLDoIgngn09+08fewT1aY/\nO0EQ+HLqKHy/ZS+Onb2itDmFyM4zgiAIT2rNeAVy3H4bKpgtoKTf9jXkkzJF478hOppasvKjiQaG\nrjZ6ByCHLGe9/iwm/m85eLN6GhOfOJ8qGjle0RyXCpIBYCyAHpAXsrUADEHJ4o93vFBS8BKAH3Gv\nolkvAEt89Lr14558ROepi9SyGPJIRzQJV5drcyDxfL5ZsFZ0MagGuIL/hkDW0V8oHoipchc/V3EM\nQP3a/j7dpr1cfVyZohAEgS/eGoU1v+1FYlKK0uYAAE6cTyUhu5eeShKAdpBPylkhz/QllntQk+Cz\n/Qxs9BdvjdJUN1emKMGB/nhtVH8sXLtNaVNwKzPH3mrH07sk23Av9TgSckhzRdGb1CT4DqyOrt8r\nroXSdriFEf27YseBk7h9N1dROxLPp0CvY87As+rL2I+QFp0ZhwFIhPxpxUI+efVfxxtUI3g/A/vm\npGf6sEoennYntfwMGNjjQXz329+K2nH0zBXRxJuVNaI4JOSwtSPRkA+QJEHexT+H4qnphwCMhrzJ\nacA93/7/UYu66ggW64CRj3dXiz1uYfTQh7Fi81+KNTITRRGrNv/FmQXrBjcPXfTvTEJOF94L+Qig\nCcU3mTIAzIEcjQkGEIh7ezp20iHP7vbqCcVQRTIHqdWMHtAjVqrt76O0KW6lTXQEQgL9sWP/CfTt\n1tbt4+88eBq8YLkGub6PK3kTQBwKqsBBrmow2uH7Vsgx9WuQz75eBVA0bpsNucRflVCD4AmGpl6f\n8NSjSmfrKcLoob2w4qe/FBH8grXb8nPzudmomv/eD/LpJbuYjShePMt+sN9eBe5WCc9xSz9XNQg+\nnNRqfdvFVM+4e1n0694Wb372LWw20a0NzFKv38ahkxcJyLno9yMawNO4J+ZwyGceHMuzWAGchTxD\n2wVdlI3OsNkZqEHwsa0fCLd6WotGZ+Hvo0fdwFq4kHoDzaIauG3cReu3W6w221HIizq7mKegcFUB\ne77TPshHOlNxb7PHzo6Cfx6B4oKnKbJDXJumNct5L0Kb6HAkJqU4VfBb/z6GK+k3kZ6RifSMTPR/\nqB1GDpBrVeYaOXy/Za/VahMzAORA9qlTUbxK2Al4SO+m8qK44A0s071ds8gaFZ0pSttmkUhMSsEz\n/cpX4vHkhVScvXwV6TcycfVmJuoH18bbowcVuufPAyfB0BTC6weha7totH7g3smq6XPX8gDxE+Qc\n8xqF0oInOF5o2SY6QmEzlKV1dAS27pHPwt/MzMbpi+nyzHzjDnLyTfj8rcLHhH//+zgupWUgLCQQ\nrR8IR0yj4gdK5rz9Qolj7T58Bj/9ccho4oVJTn8hHoDSgg9naIoIqeOS87oeQ+sHwnHqYhpsNhE7\nD57GDwn7EVYvEGEhdRDbvFGx+98ZM7iEp5RNrpHDmJlLTCZeGAHZlalxKC34pk3CQ9RUsk8R/H30\n8PfR4/rtLDzbvyue7d/VJeNMn7uW582WzVBXa3i3orTgWYNKy224Gz3LuDRHvqa7MnaUXiyyLEMr\nbYMqYGgKnIsEf+7KNTw3bQFn4oWnUUNdGTtKi43R0VTNDMAXgaFJCEJVWi2VTMr12+g/YZbJxJnH\nAfjT6QN4GEoLXhCsVk9KS3UZFqsNlJPr1KRcv40+L39oyjfx79hE8TunPtxDUdqH502c4BU85PYz\nOicefElKvoZ+42dx+Sb+bcFiXeC0B3s4Sgue43izV/AAOLMFDOMcwe85ehYj357PmTjzOO/MXhil\nBZ+SfO2WR3TGdiUcLyAzKw/1qrgfkW/i8e5X6/gN2w+YOF54Bh6U4+IulBb8hew8E5Wda0QtP4PC\npijH6UtpaBpRD1U5y/v3v+cw5r0lJhNv/pXjhYkAspxnYfVB6UWraGCZ82o5va8Ux8+loLLpFfkm\nHv+ZtZJ/asqcu7fu5jyZb+KfgVfs90XpGR682bL3eFJKyx4dmtfY8GRiUgratyieQlAa127dxcrN\nf1mXbdwp2ETxN44XJsAr9DJRXPCcWThw8MSFUUB/9zY3UhGJSSkYO/yRMu+TJAl/Hz2HBeu2Gf8+\nek5DajVrjJx5LuQDGF7KgeKCB/Dv8XPJNXZ253gByVdv3TcXPjtPbpJ89PRlccVPf5lyjdztfBM3\nW5LwvbnwySMv5UANgr+Qa+TES2kZaFyB3qrVhd1HzqBZowZIuXYbvGBBVk4+TlxIlQ4kXsg/fi6Z\nuJuTT/nodRfMgmWfiRfWQu6k6A3lVhJVzKw6hvr8+UEPTZ79xkhaaVvcTb/xH0tHz1wx6mgqV6Mh\nzARB5PCC5QDHCwcgl5w4j3vNv7xUEVUIHkCEgWXOXUyYp3N3r1UlSb56C51GvJvPmy3BKH5W1IsL\nUDosaSeF1GoPbvrDEwvYVp7lm3YKBEF8A6/Y3YZaBI+cfNOn89b8nidJNcM95XgBq37eLXK8MF9p\nW2oSqhE8gO3Xb2WZjp319CK25ePnXUdAarVHIXdF8eIm1JTHIkmA5ezlq91GDOhKV+c6NUbOjBFT\nvzLeyc6bBK/g3YqaZnhYrbaFZy9fTVm5ebcy1UXdxMyFG8z5Jn4bvMldbkeN02iMnmWOHvj+Qza8\nfpDStjid/YnnMey1L7JNvNAIwF2l7alpqMmlsXNboyGkw6cux1U318bImdF/wixTVq5xJOTC/V7c\njKpcGjtmwfppdXRtZi7cYDaa+ATIXcq9KIAqBQ/Amm/in3xv/nr+7OWrStviFP7YfxLfb9lryjPx\nY5W2pSajVsEDwFnObBnTb8IsLvlqSeXEPYcDiRfw/PQFJo4XBsDrtyuKGn14R07bbGLWpj8O9RzS\nuyPlZ/C8ngknzqdi8KufmYyceTCAPUrbU9NRu+AhitIRmyjafkjY33nAQ7FUgAcdBTx08iIGTf6U\nM3HmkRKwVWl7vHiA4AHAahP/ESzW7HW//9Pz0c6tqaDafkqbVCZ/HTqNJ6fMNRk581AJ2KK0PV5k\nPELwgDzTmwVr8rqt//St5WvQtm0WQagxZClYrPho6U+WaXPWmjiz0BfALqVt8nIP9SmmbGJ89LqN\nMY1Cw5d9MM6gps2pxKQUvDhjkfHO3dyDeSb+OQA3lLbJS2E8ZoZ34LZgsS69k50rrti8u7MaZnv7\nrP7aJ6u4O9l5E80W61sA8hQzyMt98cQZ3pEYH71uY/PGYeEf/+cZQ2zzKLcOLooi/jp8Bm9+/p13\nVvcQPF3wAECSpHYyQ5HvhNYNZCeP7Os75JEOcOXJqaycfHy/dZ+4YO02zsiZM/JN/HuSJK2H96yp\n6qkOgrejBdDH30c/1WqzdXzu8e6al4f1ohs58WD48XPJWLx+B/frX0cJmiJ/zzVynwM4CK/QPYbq\nJHhHInUM9QqAl/199ETbZpFSXOsmPm2jI4nWD4SXq6zfzcxsJCal4Pi5FPHgiQv5iedTSYvFahQs\n1jkWq205gNsufxVenE51FbwdLYCmAGJZHd1JR1Pd8k180wA/gyUqtK6N1dGEXscQNEUSvGCROF6Q\njBwvXUzLoDjeDL2OOW3kzHsEi/Uw5AoCKfDO5h5NdRd8SdjfBBEA2IJ/NAC+4J8JwAV4xe3Fixcv\nXrx48eLFixcvivN/DbH0fkGcnvoAAAAASUVORK5CYII=\n",
       "prompt_number": 4,
       "text": [
        "<IPython.core.display.Image at 0x4341510>"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Check the spacegroup (symmetry) of our creation. The `spglib` library provides various functions dealing with the symmetry of crystals. For example, it has a symmetry finder, which can identify the symmetry group of the crystal. Here, we use this function to check if the structure we have created is indeed a zinc-blende cubic crystal (F-43m group)."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "print 'Space group:', spglib.get_spacegroup(cryst)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Space group: F-43m      (216)\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Defining the `Calculator`\n",
      "-----------------------\n",
      "\n",
      "Now, when we have our crystal build we need to define a `Calculator` for our computations. It is done by creation of `QuantumEspresso` object, specifying various parameters for the calculation. The same calculator may be used for a number of structures and a number of computations. It is a provider of functions which allow the crystal object to respond to questions such as: \"What is your total energy?\". You need to specify a number of parameters for the calculator which are specific to the case. The parameters used here *should not* be used in the production runs. They are defined just for the presentation purposes. For a real calculation you *need* to select the parameters according to the case you are dealing with. Unfortunately there is no simple \"rule of thumb\" for these parameters. The list of possible parameters covers a fairly complete set of parameters used by Quantum Espresso package. For the description of possible parameters you need to consult the documentation on the [Quantum Espresso website](http://www.quantum-espresso.org/), particularly the  [document describing input parameters](http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PW.html)\n",
      "\n",
      "The meaning of the parameters in greater detail:\n",
      "\n",
      "* `label` -- This is a label for the calculator. It is used as a first part of the directory name used by the calculator to store and execute the calculations.\n",
      "* `kpts` -- A list of k-vectors for the sampling of the Brillouin zone. Here it is specified as a grid size (n x m x k) which is a typical approach. It may be also specified as a list of k-vectors\n",
      "* `xc`, `pp_type`, `pp_format` -- These are the parts of the name of the used pseudopotentials. The pseudopotential name is constructed as: \n",
      "\n",
      "    (`Element_symbol`)_(`xc`)_(`pp_type`).(`pp_format`)\n",
      "\n",
      "\n",
      "* `ecutwfc` -- Cut-off energy (in Ry) for the plane waves used in the calculation. This needs to be adjusted according to the pseudopotential used.\n",
      "* `use_symmetry` -- Controls use of symmetry in the calculation. If set to True the calculator will internally extract a primitive unit cell out of the crystal and perform all calculations on the primitive unit cell."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Create a Quantum Espresso calculator for our work. \n",
      "# This object encapsulates all parameters of the calculation, \n",
      "# not the system we are investigating.\n",
      "qe=QuantumEspresso(label='SiC',                    # Label for calculations\n",
      "                       wdir='calc',                # Working directory\n",
      "                       pseudo_dir='../../pspot',   # Directory with pseudopotentials\n",
      "                       kpts=[8,8,8],               # K-space sampling for the SCF calculation\n",
      "                       xc='pz',                    # Exchange functional type in the name of the pseudopotentials\n",
      "                       pp_type='vbc',              # Variant of the pseudopotential\n",
      "                       pp_format='UPF',            # Format of the pseudopotential files\n",
      "                       ecutwfc=70,                 # Energy cut-off (in Rydberg)\n",
      "                       use_symmetry=True)          # Use symmetry in the calculation ?\n",
      "\n",
      "# Check where the calculation files will reside.\n",
      "print qe.directory"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "calc/SiC.mjVUKG\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Assign the calculator to our system\n",
      "cryst.set_calculator(qe)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 7
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Computing - single-point calculation\n",
      "------------------------------------\n",
      "\n",
      "We are ready to perform our first Quantum Espresso calculation. It is as simple as asking the crystal for its energy or stress tensor. \n",
      "\n",
      "You just need to call an appropriate function: e.g. `cryst.get_stress()`\n",
      "\n",
      "The calculation may take some time (5-30s, depending on your system). Be patient.\n",
      "\n",
      "**Note 1**: The default notation for stress tensors used here is a *Voigt* notation. Where the independent components of a symmetric stress tensor $\\sigma$ are collected into a 6-component quantity (note: it is *not* a vector in the tensor analysis sense). This is a common notation in the tensor algebra of symmetric tensors. The components of the stress tensor in the Voigt notation are:\n",
      "$$\n",
      "[\\sigma_{xx}, \\sigma_{yy}, \\sigma_{zz}, \\sigma_{yz}, \\sigma_{xz}, \\sigma_{yz}]\n",
      "$$\n",
      "\n",
      "**Note 2**: The convention for signs of external pressure is opposite to the sign of stress.\n",
      "\n",
      "**Note 3**: In case of the interrupted calculation it is easy to recover from such a crash by adding `recover = .true.` to the input file."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Run the calculation to get stress tensor (Voigt notation, in eV/A^3) and pressure (in kBar)\n",
      "print \"Stress tensor   (Voigt notation eV/A^3):\", cryst.get_stress()\n",
      "print \"External pressure                (kBar):\", cryst.get_isotropic_pressure(cryst.get_stress())*1e-3"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Stress tensor   (Voigt notation eV/A^3): "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "[ 0.02404228  0.02404228  0.02404228 -0.          0.          0.        ]\n",
        "External pressure                (kBar): -38.52\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Do the same but get the results in GPa\n",
      "# Note that this time we get the results immediately. \n",
      "# We did not change the system, so it is not necessary to repeat the calculation.\n",
      "print \"Stress tensor  (Voigt notation, GPa):\", cryst.get_stress()/ase.units.GPa\n",
      "print \n",
      "print \"Stress tensor (Tensor notation, GPa):\"\n",
      "print cryst.get_stress(voigt=False)/ase.units.GPa\n",
      "print\n",
      "print \"External pressure              (GPa):\", cryst.get_isotropic_pressure(cryst.get_stress())*1e-4"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Stress tensor  (Voigt notation, GPa): [ 3.852  3.852  3.852 -0.     0.     0.   ]\n",
        "\n",
        "Stress tensor (Tensor notation, GPa):\n",
        "[[ 3.852  0.     0.   ]\n",
        " [ 0.     3.852 -0.   ]\n",
        " [ 0.    -0.     3.852]]\n",
        "\n",
        "External pressure              (GPa): -3.852\n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Computing - a series of calculations\n",
      "------------------------------------\n",
      "One of the advantages of this system is that you can run a series of calculations automatically, as is illustrated below. \n",
      "For example we can find a minimum of energy of our crystal - which is its equilibrium lattice constant.\n",
      "To do this it is needed to modify the crystal at each turn of \n",
      "the loop and collect the results. \n",
      "Alternatively you can create a whole bunch of systems and calculators\n",
      "and run them all at once. If you have many CPUs this may \n",
      "considerably speed up the calculation. \n",
      "Here we will do a sequential run. The topic of parallel execution will be covered in the \"Remote calculation\" notebook."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# A sequential run for a series of lattice constants\n",
      "\n",
      "# We will store the results in this list\n",
      "result=[]\n",
      "\n",
      "# Our prototype crystal is just a copy of the structure defined above\n",
      "cr=Atoms(cryst)\n",
      "\n",
      "# It needs a calculator as well. This may be the same calculator or we can define a separate one.\n",
      "cr.set_calculator(qe)\n",
      "\n",
      "print \" Scale     A(A)     Energy(eV)   Pressure(GPa) \"\n",
      "print \"===============================================\"\n",
      "\n",
      "# Iterate over scales between 98% and 102% of the starting unit cell size. \n",
      "# We use 11 points in the interval\n",
      "for x in linspace(0.98,1.02,11):\n",
      "    # Modify the crystal by scaling the lattice vectors\n",
      "    cr.set_cell(cryst.get_cell()*x,scale_atoms=True)\n",
      "    \n",
      "    # Calculate energy and stress and store the results in the result list\n",
      "    result.append([x, x*cryst.get_cell()[0,0], cr.get_potential_energy(), 1e-4*cr.get_isotropic_pressure(cr.get_stress())])\n",
      "    \n",
      "    # Print it as well\n",
      "    print \"% 5.03f   % 6.04f    %+6.4f      % +8.3f \" % tuple(result[-1])\n",
      "\n",
      "# Prepare the collected data for plotting\n",
      "# This will make an array (matrix) out of a list and transpose it for easier access later\n",
      "# Transposing the matrix means that we can specify just a column to get the whole column as a vector.\n",
      "result=array(result).T\n",
      "\n",
      "# Let us save our calculated data to a file. \n",
      "# To have data in columns we need to transpose the array again.\n",
      "# This is a consequence of the row-column conventions and has no deeper meaning.\n",
      "savetxt('e+p-vs-a.dat',result.T)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " Scale     A(A)     Energy(eV)   Pressure(GPa) \n",
        "===============================================\n",
        " 0.980    4.2724    -263.2246       +10.575 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 0.984    4.2898    -263.2407        +7.403 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 0.988    4.3073    -263.2522        +4.380 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 0.992    4.3247    -263.2594        +1.501 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 0.996    4.3422    -263.2620        -1.242 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 1.000    4.3596    -263.2607        -3.852 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 1.004    4.3770    -263.2551        -6.335 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 1.008    4.3945    -263.2458        -8.694 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 1.012    4.4119    -263.2327       -10.933 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 1.016    4.4294    -263.2157       -13.058 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        " 1.020    4.4468    -263.1959       -15.062 "
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Analysis\n",
      "--------\n",
      "\n",
      "All plotting in the notebook environment is done by the [matplotlib](http://www.matplotlib.org/) graphics library. It is a very flexible and useful tool for scientific plotting. You will find an extensive documentation on the web. It is linked in the \"Help\" menu above. Below we will plot the energy and stress curves and find the lattice constant corresponding to the zero stress (minimum energy)."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Let us plot the results and save the figure\n",
      "figsize(12,7)\n",
      "\n",
      "# To make the plot nicer we define a shift to energy. \n",
      "# Rounding to second decimal digit in eV\n",
      "E0=round(min(result[2])-(max(result[2])-min(result[2]))/20,2)\n",
      "\n",
      "# Plot the result\n",
      "plot(result[1],     # Arguments (x-axis)\n",
      "     result[2]-E0,  # Values (y-axis)\n",
      "     'o-',          # Symbol and line style\n",
      "     label='Total internal energy')\n",
      "\n",
      "legend()            # Add a legend\n",
      "\n",
      "# Set the axes labels\n",
      "xlabel('Lattice vector length ($\\AA$)')  \n",
      "ylabel('Energy %+8.2f (eV)' % (-E0))\n",
      "\n",
      "# Store the figure\n",
      "savefig('e-vs-a.pdf')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAt4AAAG7CAYAAAD0RZIrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XtcVHX+x/EXCIJX1DRNMFFwFS3xgprWKpWKorJplloZ\n62qapWA3K80VKyvLShHXrEwz81LWpgUSWiFdfkiZZZuaSqKAt7bC8gLIeH5/nGJFQVBm5szA+/l4\n8Ig5nDnnPUXw4Tuf7/frYRiGgYiIiIiIOJSn1QFERERERKoDFd4iIiIiIk6gwltERERExAlUeIuI\niIiIOIEKbxERERERJ1DhLSIiIiLiBA4tvJOTk2nXrh1t2rRhzpw5pZ4TExNDmzZtCA0NZdu2bcXH\nn376aTp06MDVV1/NbbfdRkFBgSOjioiIiIg4lMMKb5vNxqRJk0hOTmbHjh2sWrWKnTt3ljgnKSmJ\nvXv3smfPHl5++WUmTpwIQFZWFq+88gpff/013333HTabjdWrVzsqqoiIiIiIwzms8M7IyCA4OJjA\nwEC8vb0ZOXIk69atK3HO+vXriY6OBqBHjx7k5eVx5MgR6tevj7e3NydPnqSoqIiTJ0/i7+/vqKgi\nIiIiIg7nsMI7NzeXFi1aFD8OCAggNze3Quc0atSIBx54gCuvvJLmzZvToEED+vbt66ioIiIiIiIO\n5+WoC3t4eFTovNJ2rM/MzGTevHlkZWXh5+fHLbfcwptvvsntt99e4rzg4GAyMzPtkldEREREpCxB\nQUHs3bu3UtdwWOHt7+9PdnZ28ePs7GwCAgIueE5OTg7+/v6kpqbSq1cvLrvsMgCGDRvGF198cV7h\nnZmZWWrhLuIK4uLiiIuLszqGyHn0vSmuTN+f4qoqOqh8IQ5rNQkLC2PPnj1kZWVRWFjImjVriIqK\nKnFOVFQUy5cvByA9PZ0GDRrQtGlT2rZtS3p6OqdOncIwDDZt2kT79u0dFVVERERExOEcNuLt5eVF\nQkICERER2Gw2xo4dS0hICIsXLwZgwoQJREZGkpSURHBwMHXq1GHp0qUAdOrUiTvvvJOwsDA8PT3p\n0qUL48ePd1RUERERERGH8zDcuFfDw8NDrSbislJTUwkPD7c6hsh59L0prkzfn+Kq7FF3qvAWERER\nESmHPepOh7WaiIiIiJSlUaNG/Prrr1bHEDlPw4YN+eWXXxxybY14i4iIiNPpd7i4qrK+N+3xPeuw\nVU1EREREROR/VHiLiIiIiDiBCm8RERERESdQ4S0iIiIi4gQqvEVERESczNPTkx9//LHUr0VGRvLG\nG284OdGlCQ8PZ8mSJVbHcBsqvEVERET+ULduXerVq0e9evXw9PSkdu3axY9XrVpV6nNSU1Np0aKF\n3TIkJSUxevToCp1rdeHr4eGBh4eHZfd3N1rHW0RERFxGYmIa8fEpFBR44eNTRExMfwYN6u20axw/\nfrz481atWrFkyRJuuOGGi7q/M1W26D1z5gyenu49DltUVISXl3uUtO79b1pERESqjMTENGJjPyQl\n5Uk2b44jJeVJYmM/JDExzanXKE1BQQFTpkzB398ff39/7rvvPgoLCzlx4gQDBw7k4MGD1KtXj/r1\n63P48GEyMjLo2bMnDRs2pHnz5kyePJnTp09X6F5nj2IvW7aM6667joceeohGjRrRunVrkpOTAZg+\nfTqffvopkyZNol69esTExACwa9cu+vXrx2WXXUa7du14++23i6/997//nYkTJxIZGUndunX55JNP\nCAwM5Pnnnyc0NJQGDRowcuRICgoKAMjLy2Pw4MFcfvnlNGrUiCFDhpCbm1uh12EYBs888wzBwcE0\nbtyYESNGFG+alJWVhaenJ8uXL6dly5Y0adKEp5566qKe+9prr9GyZUv69u3LmTNneOCBB2jSpAmt\nW7cmISEBT09Pzpw5w9tvv01YWFiJbC+88AI33XRThV6HXRluzM3ji4iIVFul/Q7v33+6AcZ5HxER\nj1X4uva4xp8CAwONjz76yDAMw5gxY4bRs2dP46effjJ++ukno1evXsaMGTMMwzCM1NRUIyAgoMRz\nt27damzZssWw2WxGVlaWERISYsybN6/46x4eHkZmZmap9w0PDzeWLFliGIZhLF261PD29jZeffVV\n48yZM8aiRYuM5s2bl3quYRjG8ePHjYCAAGPZsmWGzWYztm3bZjRu3NjYsWOHYRiGER0dbfj5+Rlf\nfPGFYRiGkZ+fbwQGBho9evQwDh06ZPzyyy9GSEiI8dJLLxmGYRg///yz8e677xqnTp0yfv/9d+OW\nW24xbrrppjLvf7Z58+YZPXv2NHJzc43CwkJjwoQJxqhRowzDMIx9+/YZHh4exvjx4438/Hzj22+/\nNXx8fIxdu3ZV+LnR0dHGyZMnjVOnThmLFi0y2rdvb+Tm5hq//vqrceONNxqenp6GzWYz8vPzjUaN\nGhk7d+4sztapUyfj3XffLTV3WfWlPepOjXiLiIiISygoKL1d4MMPa+DhQYU+UlJKv0Z+fo1KZVu5\nciX//Oc/ady4MY0bN2bmzJnFEyCNUnYz7NKlC927d8fT05OWLVsyfvx4Nm/efEn3btmyJWPHjsXD\nw4M777yTQ4cOcfTo0eKvn33/Dz74gFatWhEdHY2npyedOnVi2LBhJUa9b7rpJnr27AmAj48PADEx\nMTRr1oyGDRsyZMgQvvnmGwAaNWrE0KFD8fX1pW7dukybNq3Cr2Px4sU8+eSTNG/eHG9vb2bOnMna\ntWs5c+ZM8TkzZ87Ex8eHjh07EhoayrfffgvASy+9VO5z4+LiqFWrFr6+vrz11ltMmTKF5s2b06BB\nAx599NHify8+Pj7ceuutrFixAoDvv/+e/fv3M3jw4Ir/R7ATFd4iIiLiEnx8iko9HhFhK2UMu/SP\n/v1Lv4avr61S2Q4ePEjLli2LH1955ZUcPHiwzPN3797N4MGDueKKK/Dz82P69On8/PPPl3TvZs2a\nFX9eu3ZtoGQv+tl93vv372fLli00bNiw+GPlypUcOXKk+NzSJoKefY9atWoVX//kyZNMmDCBwMBA\n/Pz86NOnD8eOHavQ1ulZWVkMHTq0OEf79u3x8vIqzlLaa/vzvvv37y/3uWe/jkOHDpV4HBAQUCJL\ndHQ0K1euBOCNN95gxIgReHt7l/sa7E2Ft4iIiLiEmJj+BAVNL3EsKGgakyf3c+o1StO8eXOysrKK\nHx84cIDmzZsDpU9wnDhxIu3bt2fv3r0cO3aM2bNnlxittZdz733llVfSp08ffv311+KP33//nYUL\nF17S9Z9//nl2795NRkYGx44dY/PmzRiGUaHC+8orryQ5OblElpMnT3LFFVfY5blnv/YrrriC7Ozs\n4sdnfw5wzTXXULNmTdLS0li1alWFV42xNxXeIiIi4hIGDerN/PkRRETMoE+fOCIiZjB//oCLWtXE\nHtcozahRo3jyySf573//y3//+18ef/zx4uKtadOm/Pzzz/z222/F5x8/fpx69epRu3Ztdu3axaJF\niyp1/7I0bdqUzMzM4seDBw9m9+7drFixgtOnT3P69Gm+/PJLdu3aBZTeFnMhx48fp1atWvj5+fHL\nL78wa9as884p65p3330306ZN48CBAwD89NNPrF+/vkL3vdjn3nrrrcyfP5+DBw+Sl5fHnDlzzvuj\nZPTo0UyaNImaNWvSq1evCuWwNxXeIiIi4jIGDepNcvITpKbGkZz8xCUVzPa4xrkee+wxwsLC6Nix\nIx07diQsLIzHHnsMgHbt2jFq1Chat25No0aNOHz4MHPnzmXlypXUr1+f8ePHM3LkyBKFYEWXASxt\nneyzH8fGxrJ27VoaNWrElClTqFu3LikpKaxevRp/f3+uuOIKHn30UQoLC8u83oXuOWXKFE6dOkXj\nxo3p1asXAwcOvGCes8XGxhIVFUX//v2pX78+PXv2JCMjo0L/Di72uXfddRf9+/enY8eOdO3alUGD\nBlGjRo0SSyWOHj2a77//njvuuOOCr9+RPIyL/dPHhXh4eFz0X24iIiJiPf0OF0fasGEDEydOLNEe\ndOrUKZo2bcq2bdsICgoq87llfW/a43tWI94iIiIi4tby8/NJSkqiqKiI3NxcZs2axbBhw0qcs2jR\nIrp3737BotvRNOItIiIiTqff4WJPp06dok+fPuzatYtatWoxePBg5s+fT926dQEIDAzEw8OD9957\nj9DQ0Atey5Ej3iq8RURExOn0O1xclVpNRERERETcnApvEREREREnUOEtIiIiIuIEXlYHEBERkeqn\nYcOGFV7LWsSZGjZs6LBra3KliIiIiFQbERGPkZLyZCnHZ5Cc/ESZz9PkShERERGRizBmTH88PaeX\nOBYUNI3Jk/s5/N4a8RYRERGRamPiRDhwIA2bbSP5+TXw9bUxeXI/Bg3qfcHnaR1vFd4iIiIiUkFf\nfQVDhsCOHXCxrdxqNRERERERqYAzZ+Cee+CZZy6+6LYXFd4iIiIiUuUtWQI1a8Lo0dZlUKuJiIiI\niFRpP/8M7dtDSgqEhl7aNdTjrcJbRERERMoxYQL4+sL8+Zd+DXvUndpAR0RERESqrIwMeP99c0Kl\n1dTjLSIiIiJVks1mTqicMwcaNLA6jQpvEREREamiXnkFatWCO+6wOolJPd4iIiIiUuX89BN06ACb\nNkHHjpW/niZXqvAWERERkVKMGwf16sGLL9rneppcKSIiIiJyjvR02LDBNSZUnk093iIiIiJSZfw5\nofK558DPz+o0JTm08E5OTqZdu3a0adOGOXPmlHpOTEwMbdq0ITQ0lG3btgHwww8/0Llz5+IPPz8/\n4uPjHRlVRERERKqAxYuhfn0YNcrqJOdzWI+3zWajbdu2bNq0CX9/f7p168aqVasICQkpPicpKYmE\nhASSkpLYsmULsbGxpKenl7jOmTNn8Pf3JyMjgxYtWpQMrx5vEREREfnD0aNw1VXwySfmxEp7skfd\n6bAR74yMDIKDgwkMDMTb25uRI0eybt26EuesX7+e6OhoAHr06EFeXh5Hjhwpcc6mTZsICgo6r+gW\nERERETnbI4/AnXfav+i2F4dNrszNzS1RLAcEBLBly5Zyz8nJyaFp06bFx1avXs1tt93mqJgiIiIi\nUgV88QWkpMDOnVYnKZvDCm8PD48KnXfukP3ZzyssLOT9998vsz8cIC4urvjz8PBwwsPDLyqniIiI\niLi3oiJzQuXcueYSgvaQmppKamqqfS72B4cV3v7+/mRnZxc/zs7OJiAg4ILn5OTk4O/vX/x4w4YN\ndO3alSZNmpR5n7MLbxERERGpfhYtgkaNYMQI+13z3AHdWbNmVfqaDuvxDgsLY8+ePWRlZVFYWMia\nNWuIiooqcU5UVBTLly8HID09nQYNGpRoM1m1ahWjXHFKqoiIiIi4hCNH4PHHISEBKthwYRmHjXh7\neXmRkJBAREQENpuNsWPHEhISwuLFiwGYMGECkZGRJCUlERwcTJ06dVi6dGnx80+cOMGmTZt45ZVX\nHBVRRERERNzc1KkwZgy0b291kvJpy3gRERERcUuffgq33WZOqKxb17H3cunlBEVEREREHKWoCO69\nF154wfFFt72o8BYRERERt7NwITRtCsOHW52k4tRqIiIiIiJu5dAh6NjRbDVp184597RH3anCW0RE\nRETcyujREBAATz/tvHvao+502KomIiIiIiL2lpYGmze79g6VZVGPt4iIiIi4hdOnzR0qX3wR6tSx\nOs3FU+EtIiIiIm5hwQLw94dhw6xOcmnU4y0iIiIiLu/gQXNC5RdfwF/+4vz7a3KlCm8RERGRauG2\n26BVK5g925r7q/BW4S0iIiJS5X3yibkt/PffW9fbrZ0rRURERKRKO33a3KFy3jz3nFB5NhXeIiIi\nIuKy5s+HwED429+sTlJ5ajUREREREZeUkwOdOkF6OgQHW5tFrSYiIiIiUmU9+KC5brfVRbe9aOdK\nEREREXE5H30EW7bAa69ZncR+NOItIiIiIi6lsBAmTTL7u2vXtjqN/ajwFhERERGX8uKLEBQEQ4ZY\nncS+NLlSRERERFxGdjZ07my2mQQFWZ3mfzS5UkRERESqlPvvN9tMXKnothdNrhQRERERl5CSAlu3\nwvLlVidxDI14i4iIiIjlCgpg8mSIj4dataxO4xgqvEVERETEci+8AO3aweDBVidxHE2uFBERERFL\nHTgAXbrAl19Cq1ZWpymdJleKiIiIiNu77z6IjXXdotteNLlSRERERCyTnAzffgtvvml1EsfTiLeI\niIiIWOLPCZULFoCvr9VpHE+Ft4iIiIhY4rnn4KqrYOBAq5M4hyZXioiIiIjTZWVBWBh89RUEBlqd\npnyaXCkiIiIibmnKFHNSpTsU3faiyZUiIiIi4lSJifD997BmjdVJnEuFt4iIiIg4TX4+xMTAv/4F\nPj5Wp3EutZqIiIiIiNM8+yx07gwREVYncT5NrhQRERERp/jxR+jeHb7+Gq680uo0F0eTK0VERETE\nbUyZAg8+6H5Ft72ox1tEREREHO7992H3bli71uok1lHhLSIiIiIOdeoUxMbC4sVQs6bVaayjVhMR\nERERcahnnjE3y+nXz+ok1tLkShERERFxmL174Zpr4JtvICDA6jSXTpMrRURERMRlGYa5ZvfUqe5d\ndNuLerxFRERExCHWr4d9++C996xO4hocOuKdnJxMu3btaNOmDXPmzCn1nJiYGNq0aUNoaCjbtm0r\nPp6Xl8fw4cMJCQmhffv2pKenOzKqiIiIiNjRyZPmhMqFC6v3hMqzOazwttlsTJo0ieTkZHbs2MGq\nVavYuXNniXOSkpLYu3cve/bs4eWXX2bixInFX4uNjSUyMpKdO3eyfft2QkJCHBVVREREROzs6aeh\nZ0+44Qark7gOh7WaZGRkEBwcTGBgIAAjR45k3bp1JQro9evXEx0dDUCPHj3Iy8vjyJEj+Pr68umn\nn/L666+bIb288PPzc1RUEREREbGjPXtg0SL49lurk7gWh4145+bm0qJFi+LHAQEB5ObmlntOTk4O\n+/bto0mTJowZM4YuXbpw1113cfLkSUdFFRERERE7MQyYPBkefRT8/a1O41ocNuLt4eFRofPOXZbF\nw8ODoqIivv76axISEujWrRtTpkzhmWee4fHHHz/v+XFxccWfh4eHEx4eXpnYIiIiIlIJ770H2dnm\naibuLDU1ldTUVLte02GFt7+/P9nZ2cWPs7OzCThnHZlzz8nJycHf3x/DMAgICKBbt24ADB8+nGee\neabU+5xdeIuIiIiIdU6cgClT4PXXwdvb6jSVc+6A7qxZsyp9TYe1moSFhbFnzx6ysrIoLCxkzZo1\nREVFlTgnKiqK5cuXA5Cenk6DBg1o2rQpzZo1o0WLFuzevRuATZs20aFDB0dFFRERERE7eOopuO46\nUANC6Rw24u3l5UVCQgIRERHYbDbGjh1LSEgIixcvBmDChAlERkaSlJREcHAwderUYenSpcXPX7Bg\nAbfffjuFhYUEBQWV+JqIiIiIuJYffoDFi2H7dquTuC5tGS8iIiIilWIYEBEBAwfCffdZncYxtGW8\niIiIiFjunXfg0CGYNMnqJK5NI94iIiIicsmOH4f27WHFCujd2+o0jmOPulOFt4iIiIhcskcegYMH\n4Y/1MqosFd4qvEVEREQss3OnOcr93XfQrJnVaRxLPd4iIiIiYok/d6icMaPqF932osJbRERERC7a\n22/DTz/BPfdYncR9qNVERERERC7K77+bEypXrTI3zKkO1OOtwltERETE6aZOhaNHYdkyq5M4jwpv\nFd4iIiIiTvX99+aW8P/5DzRtanUa59HkShERERFxGsMwN8mZObN6Fd32osJbRERERCpk9WrIy4O7\n77Y6iXtSq4mIiIiIlOu33yAkBNauhZ49rU7jfOrxVuEtIiIi4hQPPGCOdi9ZYnUSa9ij7vSyUxYR\nERERqaL+8x944w1zYqVcOvV4i4iIiEiZDAPuvRdmzYImTaxO495UeIuIiIhImVauhOPHYfx4q5O4\nP/V4i4iIiEipjh0zJ1T++9/Qo4fVaaylyZUqvEVEREQc5r77zNHuV16xOon1NLlSRERERBxi+3az\nzUQTKu1HPd4iIiIiUoJhwD33wOOPQ+PGVqepOjTiLSIiIiIAJCamER+fwv79Xhw6VMQVV/QHelsd\nq8pQ4S0iIiIiJCamERv7IZmZs4uP3X//dGrUgEGDVHzbg1pNRERERIT4+JQSRTdAZuZsFizYaFGi\nqkeFt4iIiIhQUFB6I0R+fg0nJ6m6VHiLiIiICB4eRaUe9/W1OTlJ1aXCW0RERKSaMwwoKOiPn9/0\nEseDgqYxeXI/i1JVPZpcKSIiIlLNLV0KJ070ZulSWLx4Bvn5NfD1tTF58gBNrLQj7VwpIiIiUo1l\nZsI110BqKnToYHUa12WPulOtJiIiIiLVVFER3HEHzJihotsZVHiLiIiIVFOzZ0P9+jBpktVJqgf1\neIuIiIhUQ+npsGgRfP01eGoo1in0r1lERESkmvn9d7PFZNEiaN7c6jTVhyZXioiIiFQzY8eChwe8\n+qrVSdyHPepOtZqIiIiIVCPvvgubN8M331idpPrRiLeIiIhINXHwIHTuDOvWmUsISsVpOUERERER\nqZAzZ+Dvf4d77lHRbRUV3iIiIiLVwIIF5qTK6dPLP1ccQ60mIiIiIlXcd9/BDTeYSwgGBVmdxj2p\n1URERERELig/H26/HZ59VkW31TTiLSIiIlKFPfAAZGXB2rXmEoJyaRy+nODRo0d5++23SUtLIysr\nCw8PD1q2bEnv3r255ZZbuPzyyy948eTkZKZMmYLNZmPcuHE8/PDD550TExPDhg0bqF27NsuWLaNz\n584ABAYGUr9+fWrUqIG3tzcZGRmVeJkiIiIi1c9HH8GaNfDttyq6XUGZhffYsWPJzMxk4MCB3H33\n3VxxxRUYhsGhQ4fIyMjg1ltvJTg4mFfLWHndZrMxadIkNm3ahL+/P926dSMqKoqQkJDic5KSkti7\ndy979uxhy5YtTJw4kfT0dMD8qyI1NZVGjRrZ+SWLiIiIVH2//GKuYrJ0KVx2mdVpBC5QeMfGxtKx\nY8fzjoeEhHDDDTfwyCOPsH379jIvnJGRQXBwMIGBgQCMHDmSdevWlSi8169fT3R0NAA9evQgLy+P\nI0eO0LRpUwC1kYiIiIhcAsOACRNg+HDo18/qNPKnMidXJicnk52dfcEnl1aY/yk3N5cWLVoUPw4I\nCCA3N7fC53h4eNC3b1/CwsJ45ZVXLvwqRERERKTY8uWwaxc8/bTVSeRsZY54Hzx4kF69etGyZUtu\nu+02brnlFpo0aVLhC3tUsJGorFHtzz77jObNm/PTTz/Rr18/2rVrx1//+tfzzouLiyv+PDw8nPDw\n8ApnFBEREalqfvwRHnzQ7O/29bU6jftKTU0lNTXVrtcss/CeN28eL7zwAmlpaaxevZonnniCjh07\ncttttzFs2DDq1at3wQv7+/uXGDHPzs4mICDggufk5OTg7+8PQPPmzQFo0qQJQ4cOJSMjo9zCW0RE\nRKQ6KyqC0aPh0UfhAo0JUgHnDujOmjWr0te84Drenp6ehIeH89JLL5GTk8P999/PvHnzinuwLyQs\nLIw9e/aQlZVFYWEha9asISoqqsQ5UVFRLF++HID09HQaNGhA06ZNOXnyJL///jsAJ06cICUlhauv\nvvpSX6OIiIhItfDMM1CrFkyZYnUSKc0FlxP80/bt21m9ejVvvfUWjRs35ukKNAx5eXmRkJBAREQE\nNpuNsWPHEhISwuLFiwGYMGECkZGRJCUlERwcTJ06dVi6dCkAhw8fZtiwYQAUFRVx++23079//0t9\njSIiIiJVXkaGuS381q3gqS0SXVKZG+js3r2b1atXs2bNGjw9PRk1ahQjR46kdevWzs5YJm2gIyIi\nIgLHj0OXLvDUU+ZKJmJ/9qg7yyy8g4KCGDlyJKNGjeKqq66q1E0cRYW3iIiICIwfD4WFsGyZ1Umq\nLofuXJmZmVn8+f79+9mzZw99+/bl5MmTFBUVUb9+/UrdWEREREQqb9062LQJvvnG6iRSnnI7gF5+\n+WWGDx/OhAkTAHPlkaFDhzo8mIiIiIhc2OHD5kY5b7wBGhN1feUW3gsXLuSzzz4rHuH+y1/+wtGj\nRx0eTERERETKZhgwZozZZnLttVankYoot/D28fHBx8en+HFRUVGFN8cREREREcdYuBB++QVmzLA6\niVRUucsJ9unTh9mzZ3Py5Ek2btzIv/71L4YMGeKMbCIiIiJSih07YNYs+OIL8Pa2Oo1UVJmrmvzJ\nZrOxZMkSUlJSAIiIiGDcuHEuMeqtVU1ERESkuikogGuugXvugbvusjpN9eHQ5QTdgQpvERERqW4e\nfhh++AH+/W9wgXHQasMedWeZPd6DBg3i7bff5uTJk+d97cSJE6xZs4bIyMhK3VxEREREKi41FVas\ngFdeUdHtjsoc8T569CgJCQmsXbuWGjVqcMUVV2AYBocPH6aoqIgRI0Zw77330qRJE2dnLqYRbxER\nEakufv0VQkPh5ZdhwACr01Q/Tms1OXz4MPv37wegZcuWNGvWrFI3tRcV3iIiIlIdGAbcdhs0aQLx\n8VanqZ4cunPl2Zo1a+YyxbaIiIhIdbNyJWzfDl99ZXUSqQxNrhQRERFxYVlZ0L07pKRAp05Wp6m+\nHDq5UkRERESsZbPBnXfCQw+p6K4KKtRqIu4rMTGN+PgUCgq88PEpIiamP4MG9bY6loiIiFTAs8+C\nlxc88IDVScQeyiy8O3fuzM0338yoUaMICgpyZiaxk8TENGJjPyQzc3bxsczM6QAqvkVERFzc1q3w\n4ovmPz3Vo1AllPmfMS8vj7y8PK6//nq6devGiy++yMGDB52ZTSopPj6lRNENkJk5mwULNlqUSERE\nRCri5Em4/XZYsABatLA6jdhLmYV3gwYNmDt3Lvv37+f5559n9+7ddOnSheuvv56XX37ZmRnlEhUU\nlP6GRn5+DScnERERkYvx4IPQrRuMGGF1ErGnct+48PDwoHfv3ixatIicnBwefvhh/u///s8Z2aSS\nfHyKSj2XWtKRAAAgAElEQVTu62tzchIRERGpqA8+gA0bICHB6iRib2UW3m3btj3vmJeXFwMGDGDp\n0qUODSX2ERPTn6Cg6SWO+flNY/LkfhYlEhERkQs5cgTuuguWLwc/P6vTiL2VObly9erVzswhDvDn\nBMoFC2aQn1+DGjVsfPfdAE6f1sRKERERV2MYMHYs/OMf8Ne/Wp1GHOGCG+js3LmTgwcP0qNHD+rW\nrVt8fMOGDQwcONApAS9EG+hcvPR0+NvfICMDWra0Oo2IiIj86aWX4NVX4YsvoGZNq9PIuRy6gU58\nfDw33XQTCxYsoEOHDrz33nvFX5s2bVqlbirWueYacxH+UaPg9Gmr04iIiAjArl0wYwa8+aaK7qqs\nzBHvq666ivT0dOrWrUtWVhbDhw/njjvuYMqUKXTu3Jlt27Y5O+t5NOJ9ac6cgSFD4Oqr4ZlnrE4j\nIiJSvRUWQs+eZm/33XdbnUbKYo+6s8web8MwittLAgMDSU1N5eabb2b//v0qdt2cpye8/jp06QJ9\n+oALdA2JiIhUW3Fx0Lw5TJhgdRJxtDJbTS6//HK++eab4sd169blgw8+4Oeff2b79u1OCSeO07ix\n+XbWP/4BublWpxEREamePv0Uli2DJUvAw8PqNOJoZbaaZGdn4+3tTbNmzUocNwyDzz//nOuuu84p\nAS9ErSaVN3s2pKTARx+BV5nvf4iIiIi9HTsGoaGwcCEMGmR1GimPPerOC65q4upUeFeezQYDBkCv\nXjBrltVpREREqo/Ro6FePfjXv6xOIhXh0FVNLmSQ/iyrMmrUgDfeMJcv+vhjq9OIiIhUD6tXw5df\nwty5VicRZ7qkEe+DBw/SvHlzR+S5KBrxtp9NmyA6Gr7+Gpo2tTqNiIhI1XXgAISFmdvCd+1qdRqp\nKLWaqPC2qxkzYMsWSE42Vz4RERER+7LZoG9fiIiARx6xOo1cDIe2mmRnZzNu3DgeeeQR8vLyGDNm\nDFdddRWjR4/m6NGjlbqpuKaZMyE/H+bMsTqJiIhI1fTCC+Z+Gg89ZHUSsUKZhfff//53QkND8fPz\n45prrqFt27YkJSXRvXt3Jk6c6MyM4iReXrByJcyfD599ZnUaERGRqmXbNnjuOVi+3JxjJdVPma0m\nnTp1Kl7H+8orr+TAgQOlfs1KajVxjMREmDjR/AFx2WVWpxEREXF/p06Z/dzTp8Ptt1udRi6FQ1tN\nzr7w6NGjS3zNZrNV6qbi2gYNghEjzMmW+rtGRESk8qZOhU6dVHRXd2UW3lFRUfz+++8AzJ49u/j4\nnj17aNu2reOTiaWeegp+/hlefNHqJCIiIu5twwZYv17rdYtWNZELyMqCHj3g/fehe3er04iIiLif\nn34yR7rffBPCw61OI5Xh0FaTLVu2cOzYMQBOnjzJP//5TwYPHszDDz9cfFyqtsBAeOklGDkS8vKs\nTiMiIuJeDAPuugvuuENFt5jKLLzHjBlDnTp1AIiNjeW3337jkUceoVatWowZM8ZpAcVaQ4fC4MEw\nbpz6vUVERC7Gq6/C/v3wxBNWJxFXccHJlV5eXgBs3bqVefPmcd111xEXF0dmZmaFLp6cnEy7du1o\n06YNc8pYHDomJoY2bdoQGhrKtm3bSnzNZrPRuXNnhgwZUtHXIw7w3HOwbx8sWmR1EhEREfewZw9M\nm2Yu01uzptVpxFWUWXh36NCB1157DYDQ0FC+/PJLAHbv3k3NCnwH2Ww2Jk2aRHJyMjt27GDVqlXs\n3LmzxDlJSUns3buXPXv28PLLL5+3Pvj8+fNp3749Hh4eF/3CxH58fGDNGoiLAxdYRVJERMSlnT5t\nrl4SFwchIVanEVdSZuH96quvsnnzZlq3bs2OHTvo1asXrVq1Yty4cbz66qvlXjgjI4Pg4GACAwPx\n9vZm5MiRrFu3rsQ569evJzo6GoAePXqQl5fHkSNHAMjJySEpKYlx48ZpAqULCA6G+Hi49Vb4Y7Eb\nERERKcUTT0DjxnDPPVYnEVfjVdYXGjRowOuvv86xY8fYt28fRUVFBAQE0KxZswpdODc3lxYtWhQ/\nDggIYMuWLeWek5ubS9OmTbnvvvt47rnn+O233y72NYmDjBwJH38Md98NK1aA3ogQEREp6fPP4ZVX\nzE3o9HtSzlXmiPef/Pz86NSpE2FhYRUuuoEKt4ecO5ptGAYffPABl19+OZ07d9Zot4uZNw+2b4c/\nupBERETkD7/9BqNHw+LFcBElk1QjZY54b9++nfHjx5OTk0NkZCRz5syhYcOGAHTv3p2MjIwLXtjf\n35/s7Ozix9nZ2QQEBFzwnJycHPz9/XnnnXdYv349SUlJ5Ofn89tvv3HnnXeyfPny8+4TFxdX/Hl4\neDjhWq/HoWrXhrfegt69zTW+r7rK6kQiIiKuISYG+vWDqCirk4g9pKamkpqaatdrlrmBzrXXXsuM\nGTPo0aMHS5Ys4bXXXmP9+vUEBwfTuXPn81YgOVdRURFt27blo48+onnz5nTv3p1Vq1YRctYsg6Sk\nJBISEkhKSiI9PZ0pU6aQnp5e4jqbN29m7ty5vP/+++eH1wY6lnn9dXj2WcjIgD9WnRQREam23n4b\npk83W0z0e7FqskfdWeaI9++//86AAQMAePDBB+natSsDBgxgxYoVFbuwlxcJCQlERERgs9kYO3Ys\nISEhLF68GIAJEyYQGRlJUlISwcHB1KlTh6VLl5Z6La1q4nqio81+75gYWLLE6jQiIiLWycmBSZPg\ngw9UdMuFlTniHRoaSlpaGn5+fsXHtm/fzrBhw/j111/5+eefnRayLBrxttbx4xAWBo89Zu7KJSIi\nUt2cOQP9+5s7Uz72mNVpxJEcumX81KlT2bFjR4ljHTt25OOPP2bYsGGVuqlUDXXrmv3e990Hu3db\nnUZERMT55s2D/Hx49FGrk4g7KHPE2x1oxNs1LF5s7mqZng6+vlanERERcY7t2+HGG835Tq1aWZ1G\nHM2hI97nevrppyt1I6m6xo+Htm3hgQesTiIiIuIc+fnm7pTPP6+iWyquwoX3W2+95cgc4sY8PODl\nlyE5GdautTqNiIiI4z36qLkd/OjRVicRd1LmqiYiF8PPD9asgchI6NIFWre2OpGIiIhjpKSYA03f\nfqvdKeXiXLDHOzAwsHgpv4MHD9K8eXPzSR4e/Pjjj85JeAHq8XY98+eb28l//jnUrGl1GhEREfv6\n+WcIDTX3s7jxRqvTiDPZo+6s8OTKimya42wqvF2PYcDQoeaI9wsvWJ1GRESk8hIT04iPTyE/34td\nu4q49tr+vPtub6tjiZM5dAMdkUvh4QGvvWa2m4SHa9tcERFxb4mJacTGfkhm5uziY99+O53ERBg0\nSMW3XJwKT6689tprHZlDqpBGjWDVKrjrLjhwwOo0IiIily4+PqVE0Q3w44+zWbBgo0WJxJ1VuPBO\nSEhwZA6pYnr2NJcXHDUKTp+2Oo2IiMilKSgovTkgP7+Gk5NIVVDhwlvkYj34INSvDzNnWp1ERETk\n0nh4FJV63NfX5uQkUhWo8BaH8fSE5cvNjw8/tDqNiIjIxfntNzhwoD8NG04vcTwoaBqTJ/ezKJW4\nM20ZLw63eTOMHAlbt8IfK1KKiIi4tIICc2+KNm1g8OA0EhI2kp9fA19fG5Mn99PEymrIqcsJuiIV\n3u7j8cfhk09g0yaoobY4ERFxYTabOWBkGObmcPq9JWCfulOtJuIU06ebrSdPPml1EhERkbIZBkye\nbG6Us2KFim6xL414i9McOgRdu8Kbb8L111udRkRE5HyzZsG6dZCaai4QIPInp4x433///Xz//feV\nuokIwBVXmFvsjh4NR49anUZERKSkl16CN96ADRtUdItjlFt4h4SEMH78eLp3785LL73EsWPHnJFL\nqqh+/SA6Gu68E86csTqNiIiIae1aeOIJSEmBpk2tTiNVVYVbTXbt2sWyZctYuXIl1113HXfddRfX\nW9wvoFYT91RUZG4nP2QIPPyw1WlERKS6+/hjczJlSgp06mR1GnFVTptcabPZ2LVrFzt37qRJkyaE\nhobywgsvMGLEiErdXKonLy9zS/kXX4QvvrA6jYiIVGdff20W3W+9paJbHK/cEe/77ruP999/nxtu\nuIFx48bRvXv34q+1bduWH374weEhy6IRb/f2/vswaRJs2waNGlmdRkREqpu9e6F3b1i4EIYOtTqN\nuDqnrOO9dOlSbr31VurUqXPe1/Ly8mjQoEGlAlSGCm/398AD5g++994DDw+r04iISHVx6BBcdx08\n8gjcdZfVacQdOKXw3rp1Kx7nVER+fn60bNkSLy+vSt28slR4u7/CQvMH3223wZQpVqcREZHqIC/P\nnGs0fDg89pjVacRdOKXwvuaaa9i6dSsdO3YE4LvvvqNDhw4cO3aMRYsWERERUakAlaHCu2rYtw96\n9IDEROjWzeo0IiJSlZ06BQMGQMeOEB+vd1ul4pwyubJ58+Z88803bN26la1bt/LNN9/QunVrNm7c\nyNSpUyt1cxGAVq1g0SJzcotWqxQREUcpKjLfYW3eHObPV9Etzldu4f3DDz/QoUOH4sft27dn165d\nBAUFndeCInKpbr4ZBg6EcePM7XpFRETsyTBg4kQ4ccLczM2zQuu6idhXud92HTp0YOLEiWzevJnU\n1FTuuece2rdvT0FBAd7e3s7IKNXE3LnmRMvFi61OIiIiVc2MGfDNN/DOO1CzptVppLoqt8f71KlT\nLFy4kM8//xyAa6+9lnvuuQdfX19OnDhBvXr1nBK0NOrxrnp274Zrr4VNmyA01Oo0IiJSFcTHm0sG\nfvYZNGlidRpxVw6fXFlUVES/fv345JNPKnUTR1HhXTWtXAmzZsHWrVC3rtVpRETEna1aBVOnwqef\nQmCg1WnEnTl8cqWXlxeenp7k5eVV6iYiF+O22+CvfzV78fR3lYiIXKqUFHOp2qQkFd3iGspdiLtO\nnTpcffXV9OvXr3gTHQ8PD+Lj4x0eTqqv+HhzacFly2DMGKvTiIiIu8nIgDvugHffhauvtjqNiKnc\nwnvYsGEMGzaseAUTwzC0mok4XO3a8NZb5gYHPXpA+/ZWJxIREXexaxf87W+wZIm5SZuIqyh3ciXA\nyZMnOXDgAO3atXNGpgpTj3fV99pr8MIL5shF7dpWpxEREVeXk2MW2zNn6h1TsS+nbKCzfv16Onfu\nzIABAwDYtm0bUVFRlbqpSEWNGQOdOkFsrNVJRETE1f3yi7kr5T33qOgW11Ru4R0XF8eWLVto2LAh\nAJ07d+bHH390eDARMHcVW7QI0tLM1U5ERERKc/IkDBkCERHw0ENWpxEpXbmFt7e3Nw0aNCj5JG33\nJE5Urx6sWWOOeu/ZY3UaERFxNadPw4gR0Lo1PPectoIX11WhnSvffPNNioqK2LNnD5MnT6ZXr17O\nyCZSrFMnePxxuPVWyM+3Oo2IiLgKw4C77gKbzZwXpLFBcWXlfnsuWLCA77//Hh8fH0aNGkX9+vWZ\nN2+eM7KJlHD33RAcrLcQRUTkfx5+GH74Ad5+G7y9rU4jcmEVWtXEVWlVk+onLw+6dDHfSrz5ZqvT\niIiIlZ5/3lwy8NNP4bLLrE4jVZ3Dt4wH+OGHH5g7dy5ZWVkUFRUV3/jjjz+u1I3tQYV39fTllzBo\nEGzZAq1aWZ1GRESssHw5zJgBn30GLVpYnUaqA6cU3h07dmTixIl06dKFGjVqFN+4a9eu5V48OTmZ\nKVOmYLPZGDduHA8//PB558TExLBhwwZq167NsmXL6Ny5M/n5+fTp04eCggIKCwv529/+xtNPP31+\neBXe1daLL8Lq1eYoR82aVqcRERFnSkqCf/wDPvkEQkKsTiPVhVMK765du7J169aLvrDNZqNt27Zs\n2rQJf39/unXrxqpVqwg56/+QpKQkEhISSEpKYsuWLcTGxpKeng6Ym/bUrl2boqIirrvuOubOnct1\n52w/pcK7+jIMiIqCtm1h7lyr04iIiLP83/+ZP//ffx+uucbqNFKdOGUDnSFDhrBw4UIOHTrEL7/8\nUvxRnoyMDIKDgwkMDMTb25uRI0eybt26EuesX7+e6OhoAHr06EFeXh5HjhwBoPYf2xQWFhZis9lo\n1KjRRb84qbo8PGDZMnNb+Q8+sDqNiIg4w/ffw003mW0mKrrFHXmVd8KyZcvw8PBg7jnDivv27bvg\n83Jzc2lxVtNVQEAAW7ZsKfecnJwcmjZtis1mo2vXrmRmZjJx4kTat29foRck1cdll8GqVeYky6++\ngoAAqxOJiIijHDgAAweaEyoHDrQ6jcilKbfwzsrKuqQLe1Rw9fpzh+z/fF6NGjX45ptvOHbsGBER\nEaSmphIeHn7e8+Pi4oo/Dw8PL/UcqbquvRamTIFRo8xeP69yv6NFRMTd/Pe/5o6U990Hd9xhdRqp\nLlJTU0lNTbXrNcssU5599lmmTp0KwNtvv80tt9xS/LVp06bx1FNPXfDC/v7+ZGdnFz/Ozs4m4Jwh\nyXPPycnJwd/fv8Q5fn5+DBo0iK+++qrcwluqp6lTzaL79tvTyMtLoaDACx+fImJi+jNoUG+r44mI\nSCUcPw6DB8Pf/mYW3iLOcu6A7qxZsyp9zTJ7vFetWlX8+blF9oYNG8q9cFhYGHv27CErK4vCwkLW\nrFlDVFRUiXOioqJYvnw5AOnp6TRo0ICmTZvy3//+l7y8PABOnTrFxo0b6dy5c8VflVQrnp5w551p\nvPPOh6SkPMnmzXGkpDxJbOyHJCamWR1PREQuUWEhDB8O7dtDKYubibgdh70x7+XlRUJCAhEREdhs\nNsaOHUtISAiLFy8GYMKECURGRpKUlERwcDB16tRh6dKlABw6dIjo6GjOnDnDmTNnGD16NDfeeKOj\nokoVsHx5Cjbb7BLHMjNns2DBDI16i4i4oTNnYMwYc8nYl182J9WLuDuHdsQOHDiQgefMgJgwYUKJ\nxwkJCec97+qrr+brr792ZDSpYgoKSv9Wzs+v4eQkIiJSWYYBDzxgTqj88EPN35Gqo8xv5e3bt1Ov\nXj3AbPf48/M/H4u4Eh+folKP+/ranJxEREQqa84c2LQJ0tLgj9WFRaqEMgtvm00Fi7iPmJj+ZGZO\nJzPzf+0mNWpMY9CgARamEhGRi7VkCSxeDJ9/Dg0bWp1GxL7K3bnSlWnnSjlbYmIaCxZsJD+/Br6+\nNjp37sdrr/Xm/fehe3er04mISHnWrYO774bNm+Evf7E6jUhJTtky3pWp8JbyJCaak3PWroXemmMp\nIuKyPv0Uhg2DpCTo1s3qNCLnc8qW8SLubNAgc3fL4cMhJcXqNCIiUprt282f0ytXquiWqk2Ft1R5\nN94I//43jB4N69dbnUZERM62bx9ERkJ8PPTrZ3UaEcfSAj1SLVx7rfn25aBBcOoUjBhhdSIRETl6\n1NwK/pFH9HNZqgcV3lJtdO0KGzfCgAFw8qTZ+y0iItb4/XdzpHvECJg0yeo0Is6hwluqlauvhk8+\ngb59zeL73nutTiQiUv0UFMDQoeaAyOOPW51GxHlUeEu185e/mJsy3HijWXw/9JDViUREqg+bDe68\nE/z84F//0lbwUr2o8JZqKTDQLL779oUTJ2DmTP3wFxFxNMOA2Fizt3vDBqhRw+pEIs6lwluqLX9/\nc5OGfv3M4vvZZ1V8i4g40hNPwGefmT97fX2tTiPifFpOUKq1yy83e743bzYn95w5Y3UiEZGq6aWX\n4PXXITnZbDMRqY60c6UI8NtvMHgwBAXBq6/q7U8REXt65x2YPNncnTIoyOo0IpdGW8ar8BY7OnHC\nnGXfsCGsWAHe3lYnEhFxf598Yi4Z+OGH0Lmz1WlELp22jBexozp1zJ0tT52Cm2+G/HyrE4mIuLdt\n28yie80aFd0ioMJbpARfX/Mt0Vq1ICrKXG5QREQu3t695m7BixbB9ddbnUbENajwFjmHtzesXGmu\nejJggNn/LSIiFXf4sLkV/MyZ5juIImJS4S1Siho1YMkSc6fLvn3hl1+sTiQi4h6OHTMHLaKjYcIE\nq9OIuBZNrhS5AMOAqVMhJQU2bjSXHxQRkdLl55tF91VXwYIF2htBqhataqLCW5zAMGDWLFi9Gj76\nyGxBERGRkmw2uOWW/7XraVlWqWrsUXdq50qRcnh4QFycuepJ796waRO0amV1KhER6yUmphEfn0J+\nvheZmUVcdll/MjJ6q+gWKYMKb5EKeughqF0b+vQxi++//MXqRCIi1klMTCM29kMyM2cXH6tZczqb\nNsGgQb0tTCbiujS5UuQi3Huv2XZy/fXwn/9YnUZExDrx8Sklim6Afftms2DBRosSibg+jXiLXKQx\nY8yR7759ITERuna1OpGIiPOdOFF6CZGfrz4TkbKo8Ba5BCNGmJvsDBwI//43XHut1YlERJzn22/h\n66+LSv2ar6/NyWlE3IdaTUQuUVQUrFgBQ4eaq52IiFQHy5eb7/jde29/goKml/haUNA0Jk/uZ1Ey\nEden5QRFKiktDYYPh6VLze2RRUSqooICuO8+c3L5u++aa3UnJqaxYMFG8vNr4OtrY/LkfppYKVWW\n1vFW4S0uIiMDhgyBhQvNIlxEpCrJzjZ/tvn7m4MMfn5WJxJxPnvUnWo1EbGD7t3N3S1jYuCNN6xO\nIyJiP5s2mT/jbr4Z3nlHRbdIZWhypYidhIaavd79+8OpUzB+vNWJREQu3Zkz8Mwz5tbvb74JN9xg\ndSIR96fCW8SOQkIgNdWceHTihNkPKSLibvLyIDoajh6FL7+EgACrE4lUDWo1EbGzoCDYvBkWLYLZ\ns8s/X0TElWzfDt26wZVXmj/LVHSL2I9GvEUc4MorzdVO/hz5nj0bPDysTiUicmFvvAH33w/z5sHt\nt1udRqTqUeEt4iDNmpltJxERZvH94ovgqfeYRMQFnb1U4Mcfw9VXW51IpGpSGSDiQI0bmxMuv/wS\nJkwAmzZ0ExEXk50NffrAoUPmzyoV3SKOo8JbxMEaNDCXGszMhDvvhNOnrU4kImL66CNzqcChQ81N\ncbRUoIhjaQMdESc5dcrcgMLHB1atMv8pImKFM2dgzhyIj9dSgSIVpZ0rVXiLmykshFGj4ORJc3Sp\nVi2rE4lIdXP2UoFvv61VS0QqSjtXiriZmjVhzRqz9zsyEn7/3epEIlKd/LlUYIsWWipQxAoOL7yT\nk5Np164dbdq0Yc6cOaWeExMTQ5s2bQgNDWXbtm0AZGdnc/3119OhQweuuuoq4uPjHR1VxCm8vOD1\n16FNG3OXy7w8qxOJSHWwYgXceCPExUFCgjkQICLO5dBWE5vNRtu2bdm0aRP+/v5069aNVatWERIS\nUnxOUlISCQkJJCUlsWXLFmJjY0lPT+fw4cMcPnyYTp06cfz4cbp27cp7771X4rlqNRF3Zhjm8l1p\naebky8aNrU4kIlVRYaH5syYlxWxx06olIpfG5VtNMjIyCA4OJjAwEG9vb0aOHMm6detKnLN+/Xqi\no6MB6NGjB3l5eRw5coRmzZrRqVMnAOrWrUtISAgHDx50ZFwRp/LwMNf2joz831JeIiL2lJMDvXtD\nbi589ZWKbhGrObTwzs3NpUWLFsWPAwICyM3NLfecnJycEudkZWWxbds2evTo4ci4Ik7n4QFPPmnu\nENe7Nxw4YHUiEakqPv7Y7Oe+6SYtFSjiKhy6c6VHBffIPnfY/uznHT9+nOHDhzN//nzq1q173nPj\n4uKKPw8PDyc8PPySsopYado0qF3bLL43bYLgYKsTiYi7MgxzqcD58//X1y0iFy81NZXU1FS7XtOh\nhbe/vz/Z2dnFj7Ozswk4Zwr1uefk5OTg7+8PwOnTp7n55pu54447uOmmm0q9x9mFt4g7mzIF6tSB\n8HCzF7N9e6sTiYi7OXbMXCrw8GFzF0qtWiJy6c4d0J01a1alr+nQVpOwsDD27NlDVlYWhYWFrFmz\nhqioqBLnREVFsXz5cgDS09Np0KABTZs2xTAMxo4dS/v27ZkyZYojY4q4jLvugmeeMUeo/ljgR0Sk\nQr77DsLCwN9fSwWKuCqHjnh7eXmRkJBAREQENpuNsWPHEhISwuLFiwGYMGECkZGRJCUlERwcTJ06\ndVi6dCkAn3/+OStWrKBjx4507twZgKeffpoBAwY4MrKI5e64w2w7GTAA1q2Da66xOpGIuLo33zTf\nNXvxRfNniIi4Ju1cKeKikpLg7383d5br08fqNCLiigoL4f774cMP4Z13oGNHqxOJVF0uv5ygiFy6\nyEhzl8tbboHkZKvTiIiryckx/yjPyTH7uVV0i7g+Fd4iLuz66+G99+DOO81/iojA/5YKjIoylwps\n0MDqRCJSEQ7t8RaRyuvVyxzxjoyEU6dg1CirE4mIVQwDnn0W5s3TUoEi7kg93iJu4j//gYgIGD48\njV27Uigo8MLHp4iYmP4MGtTb6ngi4mDHjpnzPg4dMud+nLX3nIg4gT3qTo14i7iJq66CmTPTuOee\nD7HZZhcfz8ycDqDiW6QK++47uPlm6NcPVq8GHx+rE4nIpVCPt4gbeeedlBJFN0Bm5mwWLNhoUSIR\ncbQ334QbboAZM2DhQhXdIu5MI94ibqSgoPT/ZfPzazg5iYg4WmEhPPAAbNgAH32kVUtEqgKNeIu4\nER+folKP79hhY98+J4cREYfJyYHwcDhwAL76SkW3SFWhwlvEjcTE9CcoaHqJY61bT2PgwH506waP\nPw75+RaFExG7+OQTc6nAIUPg3//WUoEiVYlWNRFxM4mJaSxYsJH8/Br4+tqYPLkfgwb15sABcwe7\nbdsgPh4GDbI6qYhcDMOA556DF14wlwrs29fqRCJyNnvUnSq8RaqYDz+EyZMhJMRc67dVK6sTiUh5\njh2DMWMgNxfWrtVSgSKuSFvGi8h5IiLMpceuuQa1n4i4gf/8x/x/tVkzSEtT0S1SlanwFqmCfHzg\n0Ufh66/h22/NNcATE61OJSLnWrkSrr/eXCrwX//SUoEiVZ1aTUSqAbWfiLiWwkJ48EFISoJ33oHQ\nUN+UhU0AABeVSURBVKsTiUh51GoiIhXyZ/tJjx5qPxGxWm6uuVRgVpa5VKCKbpHqQ4W3SDXh4wPT\npsHWrf9rP0lKsjqVSPXy51KBgwfDe+9pqUCR6katJiLVVHIyxMRA+/Zm+0lgoNWJRKouLRUo4v60\nnKAKb5FKKSiA5583i4HYWHjoIfD1tTqViPtLTEwjPj6FggIvatQo4tSp/hQV9WbtWrjySqvTicil\nUOGtwlvELvbvh/vug+3bzc13IiOtTiTivhIT04iN/ZDMzNnFx+rXn86yZREMHdrbwmQiUhmaXCki\ndtGyJbz7LiQkmCPfN91kTvwSkYsXH59SougG+O232SxevNGiRCLiKlR4i0ixAQP+t5lHWNj/t3f/\nQVXV+R/HXwgIu+iIv8ICNhNpBUG4SRq5aJY/CpXCyaKaVSvTZUatrTYzaya1dtZ13fJHtemWm6PD\ntuks6IpoP0StdFFiMwfXBRcU2NjW3yKCguf7x/ly5cIFAbnnXuT5mLkT93DOvW9unzm+OLzP5yMt\nXszsJ0BrVFdLxcU+Tr9XVeVtcTUAPA3BG4ADPz9pwQJz9pO8PHP2k23b3F0V4NnOnZN+9zspLEw6\nebLG6T7+/rUWVwXA0xC8AThV136ycqU5+0lyMu0nQEPl5eY0nQMGmL+sbtkiffTROIWFLXDYLyzs\nFc2ZM9ZNVQLwFM7/HgYA/++BB8zFd5YtM9tPnnvOXHGP2U/QmRUUmFe4//IX6fHHpZwcM3xLks1m\n3kC5cuVrqqrylr9/rebMuV8TJnBjJdDZMasJgBarP/vJypVmKAc6kwMHpCVLpOxsKTVVmjNH6tvX\n3VUBsALTCRK8AbfYts1sP4mKkt56i8V3cGMzDGnHDjNwFxZKzz8vzZghdevm7soAWInpBAG4RV37\nSVyc+XjjDWY/wY2npkZKS5PuuEN64QVp+nTp6FGz3YrQDaAtuOIN4LoUF5vtJ4cOmYvv0H6Cjq6y\nUlq71ryvISREmjfPXFTKy8vdlQFwJ1pNCN6Ax6D9BB3dyZPSu++aC0nFx5uBOz7e3VUB8BS0mgDw\nGHXtJ0OH0n6CjuX4cbN9JDxcKioyb5xMTyd0A2h/BG8A7cbfX3r1VXPmh9xcKTpayspyd1WAc4cO\nSVOnSjab5Otr/uL44YdSRIS7KwNwo6LVBIDLZGaa7SdDhpjtJ7fe6u6K0NkZhvTll+YMJbm55vhM\nTZUCA91dGQBPR6sJAI+WmGheVbzjDrMF5c03pepqd1eFzujKFbN9ZMQI6amnpEmTzLaS+fMJ3QCs\nwxVvAJYoLjb7aPPzzdlP7r/f3RWhM6iuljZskJYulQICzBsmJ0+WvL3dXRmAjoZZTQjeQIdD+wms\ncO6ctHq19Pbb0uDBZuAePZopAQG0Ha0mADqcuvYTm432E7S/8nKzfWTAALOHe8sWaft26d57Cd0A\n3I/gDcBy/v7Sa69J+/ebD2Y/wfUqKJBmzZIiI6Xz56WcHHPVSZvN3ZUBwFUEbwBuc9tt5g1vb70l\nzZ5t9t4eO+buqtCR7N8vTZki3X23FBQkHTliLoAzYIC7KwOAxujxBuARqqrMG+CWL5eef16KiNit\nP/xhh6qrfeTnV6O5c8dpwoSR7i4THsAwpB07zCkBCwvN8TJjhtStm7srA3Aj4+ZKgjdwwykqklJS\nduubb7arpuZN+/awsAVavnw84bsTq6mRPvlE+u1vpcuXpZdekh57zFz8BgBcrUPcXJmVlaVBgwYp\nPDxcS5YscbrP3LlzFR4erpiYGOXl5dm3P/XUUwoKClJ0dLSrywTgIW67TQoM3OEQuiXp6NE3tXLl\np26qCu5UWSm98450++3Se+9Jb7xhrjI5dSqhG0DH4tLgXVtbq9mzZysrK0v5+flKS0vT4cOHHfbJ\nzMxUYWGhCgoKtHr1aqWmptq/9+STTyqLO66ATqe62sfp9s8+89YDD0iLF0uffWZOGYcb18mT0qJF\n5i9jn35qzse9e7c0YQIzlADomFwavHNycjRw4ED1799fvr6+SklJUUZGhsM+mzdv1rRp0yRJw4cP\n15kzZ1ReXi5JSkhIUM+ePV1ZIgAP5OdX43R7QkKtZs40Z61YuFC6+WZzPvBZs6SPPpL+9S+z/xcd\n2/Hj5mJL4eHmzbbZ2eZNuPHx7q4MAK6PS4N3WVmZQkND7c9DQkJUVlbW6n0AdC5z545TWNgCh21h\nYa/oxRfHKjnZ7PHds0c6fVr64x/NKeS2bZPGjpX69pUmTjTnB9+5U6qocNMPgVY7dMhsH7HZzBaS\n776TPvhAiohwd2UA0D6c/z23nXi18G+BDRvVW3qcJL3++uv2r++55x7dc889LT4WgGequ4Fy5crX\nVFXlLX//Ws2Zc3+jGyu7dpWGDTMfzz5rbisrk/buNR8LFkjffmv2BsfHm1POxcebU83RqmC9rVt3\na8UKx5lqEhNHas8e85ep3FxzVdMVK6TAQHdXC6Czy87OVnZ2dru+pkuDd3BwsEpKSuzPS0pKFBIS\n0uw+paWlCg4ObvF71A/eAG4cEyaMbNMMJsHB0sMPmw/JXBUzL0/6+mspI0N6+WXp0iUzgNeF8bg4\n6cc/bucfAA62bt2tZ5/drqNHr940e/DgAvXoIdXWjtSLL0obN5qLKwGAJ2h4QXfhwoXX/ZouDd5x\ncXEqKChQcXGxbrnlFn388cdKS0tz2CcpKUmrVq1SSkqK9u3bp8DAQAUFBbmyLACdiJ+fdNdd5qNO\nSYkZxPfulX71K7PFISLCMYzfeitXxdvTihU7HEK3JJWXv6mgoNeUmztS3t5uKgwALOTS4O3j46NV\nq1Zp/Pjxqq2t1dNPP62IiAi9//77kqRZs2YpMTFRmZmZGjhwoAICArR27Vr78Y899ph27dqlkydP\nKjQ0VIsWLdKTTz7pypIBdAKhodKjj5oPSbp4UfrmGzOMb9wovfCCub1+EB86lKuxzamtlb7/3vyl\n5vhx8791j+PHpW+/df7PTWCgN6EbQKfBAjoA0IBhmLNp1PWKf/21dPiwFBXlGMbr3Rd+QzMM6cSJ\nxmG6/tfl5VKfPuZnUvf4yU+ufj1v3qvateuNRq89fvxryspa7IafCgBah5UrCd4ALFJZKR044BjG\nu3a9GsTj46U77jBbWzqa8+edh+m6r0tLzav9DcN0/a+Dg83PoynOerzDwl7R8uWNb5oFAE9E8CZ4\nA3ATw5D+/e+rQXzvXunIESkmxjGMt+JecZeorjaDc1MtICUl5vLrzsJ03fOQEKlbt+uvZevW3Vq5\n8tN6M9WMJXQD6DAI3gRvAB6kokLav98xjAcEOAbx2NjGV4adTbPXkkBav6+6qRaQ06fN8N9UC0ho\nqNSrFzeSAsC1ELwJ3gA8mGFIhYVXZ1DZu1c6etRcIKYuiFdU7NbChQ1bMBbo7bfH6667RjZ5lbqk\nxAzdvXs33wISFCRuXgSAdkDwJngD6GDOn5dycq6G8U8/fVU1NY1vOuzS5TX16LG42RaQa/VVAwDa\nT3vkTpdOJwgAcNS9u3TffeZDkkaN8tHu3Y33u/tub+3ZY21tAADX6uLuAgCgM/P3r3G6PSCg1uJK\nAACuRvAGADeaO3ecwsIWOGwLC3tFc+aMdVNFAABXoccbANyMafYAwPNxcyXBGwAAABZoj9xJqwkA\nAABgAYI3AAAAYAGCNwAAAGABgjcAAABgAYI3AAAAYAGCNwAAAGABgjcAAABgAYI3AAAAYAGCNwAA\nAGABgjcAAABgAYI3AAAAYAGCNwAAAGABgjcAAABgAYI3AAAAYAGCNwAAAGABgjcAAABgAYI3AAAA\nYAGCNwAAAGABgjcAAABgAYI3AAAAYAGCNwAAAGABgjcAAABgAYI3AAAAYAGCNwAAAGABgjcAAABg\nAYI3AAAAYAGCNwAAAGABgjcAAABgAYI3AAAAYAGCNwAAAGABlwbvrKwsDRo0SOHh4VqyZInTfebO\nnavw8HDFxMQoLy+vVccCniw7O9vdJQBOMTbhyRifuJG5LHjX1tZq9uzZysrKUn5+vtLS0nT48GGH\nfTIzM1VYWKiCggKtXr1aqampLT4W8HT84wFPxdiEJ2N84kbmsuCdk5OjgQMHqn///vL19VVKSooy\nMjIc9tm8ebOmTZsmSRo+fLjOnDmj8vLyFh0LAAAAdCQuC95lZWUKDQ21Pw8JCVFZWVmL9vnPf/5z\nzWMBAACAjsTHVS/s5eXVov0Mw2jze4SFhbX4fQB3WLhwobtLAJxibMKTMT7hicLCwq77NVwWvIOD\ng1VSUmJ/XlJSopCQkGb3KS0tVUhIiC5fvnzNYyWpsLDQBZUDAAAA7c9lrSZxcXEqKChQcXGxLl26\npI8//lhJSUkO+yQlJWndunWSpH379ikwMFBBQUEtOhYAAADoSFx2xdvHx0erVq3S+PHjVVtbq6ef\nfloRERF6//33JUmzZs1SYmKiMjMzNXDgQAUEBGjt2rXNHgsAAAB0VF7G9TRZAwAAAGgRj165sra2\nVjabTZMmTWr0vQ0bNigmJkZDhgzRiBEjdPDgQfv3+vfvryFDhshms2nYsGFWloxOormxmZGRoZiY\nGNlsNg0dOlRffPGF/XssDAUrtHV8cu6EqzU3Nuvs379fPj4+2rRpk30b5064WlvHZmvPmy5rNWkP\ny5cvV2RkpM6fP9/oewMGDNDu3bvVo0cPZWVlaebMmdq3b58kc0aV7Oxs9erVy+qS0Uk0NzbHjBmj\nBx98UJL03XffKTk5WYWFhfaFoT777DMFBwfrzjvvVFJSEm1UaHdtGZ8S5064XnNjUzLDz7x583T/\n/fc7bOPcCVdry9iUWn/e9Ngr3qWlpcrMzNSMGTOcTjkYHx+vHj16SDIX3yktLXX4Ph00cJVrjc2A\ngAD71xUVFerTp4+kli0qBVyvto7POpw74SrXGpuStHLlSj388MPq27evfRvnTrhaW8dmndacNz02\neP/yl7/U0qVL1aXLtUv84IMPlJiYaH/u5eWlMWPGKC4uTmvWrHFlmeiEWjI209PTFRERoQceeEAr\nVqyQ1LJFpYDr1dbxKXHuhGtda2yWlZUpIyNDqampkq6uB8K5E67W1rFZ93VrzpseGbz/9re/6aab\nbpLNZrvmbxE7d+7Uhx9+6NDz9dVXXykvL0/btm3TO++8oz179ri6ZHQSLR2bDz30kA4fPqwtW7bo\n5z//OVcRYYm2js86nDvhKi0Zm88995x+85vfyMvLS4Zh2PdjoTy40vWMTan1502P7PH++uuvtXnz\nZmVmZqqqqkrnzp3T1KlT7XN+1zl48KCeeeYZZWVlqWfPnvbtN998sySpb9++Sk5OVk5OjhISEiz9\nGXBjaunYrJOQkKCamhqdOnVKISEhLVoYCmirto7PkydPqnfv3pw74TItGZu5ublKSUmRJJ04cULb\ntm2Tr69vixbkA9rqesZmUlJS68+bhofLzs42Jk6c2Gj7sWPHjLCwMGPv3r0O2y9cuGCcO3fOMAzD\nqKioMO6++25j+/btltSKzqWpsVlYWGhcuXLFMAzDyM3NNQYMGGAYhmFcvnzZGDBggFFUVGRUV1cb\nMTExRn5+vqU1o/No7fjk3AmrNDU265s+fbqxadMmwzA4d8I6rR2bbTlveuQV74bq/sxUf/GdRYsW\n6fTp0/Z+G19fX+Xk5Ki8vFyTJ0+WJNXU1OiJJ57QuHHj3FM4bnjOxuamTZu0bt06+fr6qlu3bvrz\nn/8siYWhYL3WjE/OnbCSs7HZFM6dsFJrxmZbzpssoAMAAABYwCNvrgQAAABuNARvAAAAwAIEbwAA\nAMACBG8AAADAAgRvAAAAwAIEbwAAAMACBG8AAADAAgRvAAAAwAIEbwBAh1ZUVKQZM2Zox44d7i4F\nAJpF8AbQ6XXr1q1F+509e1bvvfdek88lacSIEe1aW1s5q62tWvr5tFTD2oqLixUdHd2iY6urqzVq\n1CjVX3S5rKxMycnJKi8vd9hv5MiRunLlSvsVDgDXieANoNPz8vJq0X6nT5/Wu+++2+RzSfrqq6/a\ntba2clbbtRiG4RBo67T082mpttRWZ8OGDZo4caJDTT/72c+0ZcsWTZ061b7Nz89PCQkJSk9Pv+56\nAaC9ELwBwInk5GTFxcUpKipKa9askSS9/PLLOnr0qGw2m1566SXNnz/f/nzevHmSHK8Or1u3TjEx\nMYqNjdW0adMkSevXr9fw4cNls9n0i1/8wukV2fnz5zsE09dff13Lli1r9vj67zV16lSntf3+979X\ndHS0oqOjtXz5cknm1eaf/vSnmjZtmqKjo1VaWtrkZ+LsvYuLixUREaGZM2cqKipK48ePV1VVlf2Y\nxYsXa9CgQUpISNDjjz+uZcuWNarNy8tLtbW1Tb5GfWlpaXrwwQcdth07dkzp6enaunWrw/akpCSl\npaU1+fMAgOUMAOjkunXr1mjbqVOnDMMwjMrKSiMqKso4deqUUVxcbERFRdn3afi8/msdOnTIuP32\n242TJ0/aXy8/P9+YNGmSUVNTYxiGYaSmphrr1q1r9N55eXnGqFGj7M8jIyON0tLSJo9v+F6nT59u\nVNuBAweM6Ohoo7Ky0qioqDAGDx5s5OXlGUVFRUaXLl2Mv//9781+PocPH3b63kVFRYaPj4/x7bff\nGoZhGI888oixfv16wzAMIycnx4iNjTWqq6uN8+fPG+Hh4cayZcsa1dbca9RXU1Nj9OvXr9H2pUuX\nGps3bzZGjx7tsL2qqsq45ZZbmvy5AMBqPu4O/gDgiZYvX25vUygtLVVBQYFuuukmh30MJ20Zdb74\n4gs98sgj6tWrlySpZ8+e2rBhg3JzcxUXFydJunjxovr169fo2NjYWP3www/6/vvv9cMPP6hnz54K\nDg7WX//6V4fjq6qq1K9fP509e9bhvQIDA3XmzBmH1/zyyy81efJk/ehHP5IkTZ48WXv27FFSUpJu\nvfVWDRs2rNnP4/PPP3da+8iRI3XbbbdpyJAhkqShQ4equLhYktl289BDD6lr167q2rWrJk2a1ORn\n1tRr1HfixAl1797dYVtFRYV8fX01ceJEvfjii8rLy5PNZpNktptcuXJFVVVV8vf3b/bnAwArELwB\noIHs7Gx9/vnn2rdvn/z9/TV69OgmWx+a4uXl5TRkTps2Tb/+9a+vefyUKVO0ceNGlZeXKyUlpdnj\nV61a1ewvAc7qMQzD3icdEBBwzXoMw3D63sXFxfLz87M/9/b21sWLF5t8z6Y09RrO6qhv7dq12rlz\np/7xj3/I19dXy5Yt0/r16x32b+8edQBoK3q8AaCBs2fPqmfPnvL399c///lP7du3T5LUvXt3nT9/\n3r5fw+f1jR49Wp988olOnTolybyh8L777tPGjRv1v//9T5J06tQpHT9+3Onxjz76qNLS0rRx40ZN\nmTJFkpo8/t577230Xg1rq7vR8OLFi7pw4YLS09OVkJBwzcBepzW11xkxYoS2bNmi6upqVVRUaOvW\nrfLy8mr2c2tOnz59VFFRYX9eU1OjoqIipaena+3atdq+fbsyMjJUUlIiyZzZxNvb2yHUA4A7EbwB\ndHqVlZUKDQ21P44cOaKamhpFRkZq/vz5io+PlyT17t1bI0aMUHR0tObNm9fouXR1BpDBgwdrwYIF\nGjVqlGJjY/XCCy8oIiJCb7zxhsaNG6eYmBiNGzfOYQq8+iIjI1VRUaGQkBAFBQVJUpPHR0ZGNnqv\nhrXZbDZNnz5dw4YN01133aVnnnlGMTExDjU3xcvLq9naGx5f9zwuLk5JSUkaMmSIEhMTFR0drR49\neqhXr14OtXl5eTX5GvV5e3srKipKR44ckWRe/T9w4IDOnj0rSSosLJSfn5+Sk5N14cIF5eXl2f/f\nAYAn8DJaerkDAIBWunDhggICAlRZWalRo0ZpzZo1io2NbfPr/elPf9J///tf+y86zXnllVd05513\nKjk5uc3vBwDtieANAHCZJ554Qvn5+aqqqtL06dNbFJibc+nSJY0ZM0a7du1q9kp9dXW1xo4de839\nAMBKBG8AAADAAvR4AwAAABYgeAMAAAAWIHgDAAAAFiB4AwAAABYgeAMAAAAWIHgDAAAAFiB4AwAA\nABb4P3JUaNZ3yV4gAAAAAElFTkSuQmCC\n",
       "text": [
        "<matplotlib.figure.Figure at 0x4485e50>"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Fitting the Equation of State\n",
      "-----------------------------\n",
      "\n",
      "Below we will fit the Birch-Murnaghan logarithmic equation of state to our lattice constant-pressure data to find:\n",
      "\n",
      "* Equilibrium lattice constant $A_0$\n",
      "* Bulk modulus $B_0$\n",
      "* Derivative of bulk modulus $B'_0$\n",
      "\n",
      "The Birch-Murnaghan equation of state has a following form:\n",
      "\n",
      "$$\n",
      "P(V)=\\frac{B_0}{B'_0} \\left[\\left(\\frac{V_0}{V}\\right)^{B'_0}-1\\right]\n",
      "$$\n",
      "\n",
      "To fit this formula to our data points we use a standard least-squares non-linear optimization procedure `leastsq` from the `optimize` module of the `SciPy` library. The documentation for this library is also included in the \"Help\" menu."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Lets do the same with pressure. \n",
      "# But this time let us fit a Birch-Murnaghan equation of state to the data\n",
      "\n",
      "# We need a fitting package from scipy\n",
      "from scipy import optimize\n",
      "\n",
      "# Define a B-M eos function\n",
      "def BMEOS(v,v0,b0,b0p):\n",
      "        return (b0/b0p)*(pow(v0/v,b0p) - 1)\n",
      "\n",
      "# Define functions for fitting\n",
      "# The B-M EOS is defined as a function of volume.\n",
      "# Our data is a function of lattice parameter A^3=V\n",
      "# We need to convert them on-the-fly\n",
      "fitfunc = lambda p, x: [BMEOS(xv**3,p[0]**3,p[1],p[2]) for xv in x]\n",
      "errfunc = lambda p, x, y: fitfunc(p, x) - y\n",
      "\n",
      "\n",
      "figsize(12,7)\n",
      "# Plot the data\n",
      "plot(result[1],result[3],'+',markersize=10,markeredgewidth=2,label='Pressure')\n",
      "\n",
      "# Fit the EOS\n",
      "\n",
      "# Create a data array: lattice constant vs. isotropic pressure\n",
      "ap=array([result[1],result[3]])\n",
      "\n",
      "# Estimate the initial guess assuming b0p=1\n",
      "# Limiting arguments\n",
      "a1=min(ap[0])\n",
      "a2=max(ap[0])\n",
      "# The pressure is falling with the growing volume\n",
      "p2=min(ap[1])\n",
      "p1=max(ap[1])\n",
      "\n",
      "# Estimate the slope\n",
      "b0=(p1*a1-p2*a2)/(a2-a1)\n",
      "a0=(a1)*(p1+b0)/b0\n",
      "\n",
      "# Set the initial guess\n",
      "p0=[a0,b0,1]\n",
      "\n",
      "# Fitting\n",
      "# fit will receive the fitted parameters, \n",
      "# and value of succ indicates if fitting was successful\n",
      "fit, succ = optimize.leastsq(errfunc, p0[:], args=(ap[0],ap[1]))\n",
      "\n",
      "# Ranges - the ordering in ap is not guaranteed at all!\n",
      "# In fact it may be purely random.\n",
      "x=numpy.array([min(ap[0]),max(ap[0])])\n",
      "y=numpy.array([min(ap[1]),max(ap[1])])\n",
      "\n",
      "# Plot the P(V) curves and points for the crystal\n",
      "\n",
      "# Mark the center at P=0, A=A0 with dashed lines\n",
      "axvline(fit[0],ls='--')\n",
      "axhline(0,ls='--')\n",
      "\n",
      "# Plot the fitted B-M EOS through the points,\n",
      "# and put the fitting results on the figure.\n",
      "xa=numpy.linspace(x[0],x[-1],20)\n",
      "plot(xa,fitfunc(fit,xa),'-', \n",
      "     label=\"\\nB-M fit:\\n$A_0$=%6.4f $\\AA$,\\n$B_0$=%6.1f GPa,\\n$B'_0$=%5.3f  \" % (fit[0], fit[1], fit[2]) )\n",
      "\n",
      "legend()\n",
      "xlabel('Lattice vector length ($\\AA$)')\n",
      "ylabel('Pressure (GPa)')\n",
      "\n",
      "# Save our figure\n",
      "savefig('p-vs-a.pdf')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAtwAAAG7CAYAAADwsEIWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl4TVfj9vHvjsggkgglhtCYxyDGxhhjVA3VVlGlFK02\n0hqLKomaPcXPULRKS6uG9mlLURSNeSiNeVYUJYYYImkiw3n/8MpThITkZCcn9+e6zlU5Z++z7uRK\n9Lay9tqGxWKxICIiIiIiVmFndgAREREREVumwi0iIiIiYkUq3CIiIiIiVqTCLSIiIiJiRSrcIiIi\nIiJWpMItIiIiImJFphbuN998E09PT3x8fJKeCwkJwcvLC19fX3x9fVm9erWJCUVERERE0sbUwt29\ne/eHCrVhGPTv35+wsDDCwsJo0aKFSelERERERNLO1MJdv359PDw8Hnpe9+IREREREVuRKddwT58+\nnSpVqtCjRw9u3LhhdhwRERERkadmmH1r9zNnztC6dWsOHDgAwOXLl8mfPz8Aw4cP5+LFi8ydO/e+\nc0qVKsWpU6cyPKuIiIiIZC8lS5bk5MmTaXoP+3TKkm4KFCiQ9OeePXvSunXrh445deqUlp1IphUS\nEkJISIjZMSSNQkLuPmyJvjclM9P3p2RWhmGk+T0y3ZKSixcvJv35xx9/vG8HExGRjDJypNkJRETE\nVpg6w92pUyc2btzI1atXKVq0KCNHjiQ0NJS9e/diGAbFixfns88+MzOiiIiIiEiamFq4Fy1a9NBz\nb775pglJRNKPv7+/2RFEkqXvTcnM9P0ptsz0iyafhmEYWsMtIlZlGKC/ZkREJD16Z6a7aFJEREQk\nveTNm5fr16+bHUOyAA8PDyIiIqzy3ircIiLJCA42O4GIpIfr16/rt+KSKumxG8kj31tLSkRERMRW\nqTNIaj3qeyU9vocy3baAIiIiIiK2RIVbRERERMSKVLhFRERERKxIhVtERESyLcO4+zDrfMkeVLhF\nRJIREmJ2AhGxdd7e3uTKlQtXV1cKFixI9+7diYqKMjuWWIEKt4hIMkaONDuBiNg6wzBYsWIFkZGR\n/PHHH+zevZvRo0ffd0x8fLxJ6f4nISHB7AhZngq3iIiIiMkKFy7M888/z8GDB7Gzs2PmzJmULl2a\nsmXLArBixQqqVq2Kh4cHdevW5cCBA0nnTpgwAS8vL9zc3ChXrhwbNmwAYNeuXdSoUQN3d3cKFizI\ngAEDAAgNDaVo0aL3je/t7Z10XkhICK+88gpdunTB3d2d+fPnc/PmTXr06EHhwoXx8vJi+PDhJCYm\nZsSXxiaocIuIiEi2cW/N9YNrrx/1fHqf/6B7+zufO3eOVatW4evrC8CyZcv4/fffOXz4MGFhYfTo\n0YM5c+YQERHB22+/TZs2bYiLi+PYsWN8+umn7N69m1u3brF27Vq8vb0BeP/99+nXrx83b97kzz//\npEOHDo/5vO4PvXz5ctq3b8/Nmzd57bXX6NatGw4ODpw6dYqwsDDWrl3LF198kfpPNJtT4RYREREx\ngcVi4cUXX8TDw4P69evj7+/Phx9+CMDQoUPJkycPjo6OfP7557z99tvUrFkTwzDo2rUrjo6ObN++\nHXt7e2JjYzl06BBxcXEUK1aMEiVKAODg4MCJEye4evUquXLlolatWqnOVqdOHdq0aQPAzZs3+eWX\nX5gyZQrOzs7kz5+fvn37snjx4vT/otgoFW4RERHJNiyW+x8pPZ/e5/+bYRgsW7aM69evc+bMGWbM\nmIGTkxPAfUs+zp49y6RJk/Dw8Eh6nD9/nosXL1KyZEn+7//+j5CQEDw9PenUqRMXL14EYO7cuRw/\nfpzy5ctTq1YtVq5cmeqvk5eX133jx8XFUahQoaTxe/fuzZUrV1L9ftmdCreISDKCg81OICLZ2b+X\neBQrVoxhw4Zx/fr1pMft27eTloh06tSJzZs3c/bsWQzDYPDgwQCUKlWKb7/9litXrjB48GBeeeUV\n/vnnH1xcXIiOjk56/4SEhIfK87/HL1q0KI6Ojly7di1p/Js3b963jlweT4VbRCQZ2hZQRDKLXr16\nMXv2bHbt2oXFYiEqKoqVK1dy+/Ztjh8/zoYNG4iNjcXR0REnJydy5MgBwDfffJNUpN3d3TEMAzs7\nO8qUKUNMTAyrVq0iLi6O0aNHExsb+8jxCxUqRPPmzenfvz+RkZEkJiZy6tQpNm3alCGfvy1Q4RYR\nERHJRB68gLF69erMmTOHPn36kDdvXkqXLs2CBQsAiI2NZejQoeTPn59ChQpx9epVxo0bB8CaNWuo\nVKkSrq6u9OvXj8WLF+Po6Ii7uzszZ86kZ8+eeHl5kTt37vuWsBiG8VCGBQsWcOfOHSpUqEDevHlp\n3749ly5dsvJXwnYYFktqVxplHoZhkAVji4iISAZTZ5DUetT3Snp8D2mGW0RERETEilS4RURERESs\nSIVbRCQZumhSRETSi9Zwi4gkwzBSv5euiGRe6gySWlrDLSIiIiKSRalwi4iIiIhYkQq3iIiIiIgV\nqXCLiIiIiFiRCreISDKCg81OICIitkK7lIiIiIjNUmeQ1NIuJSIiIiIiWZQKt4iIiIiIFalwi4iI\niIhYkQq3iIiIiIgVqXCLiCQjJMTsBCIiYiu0S4mISDIMA/TXjEjWp84gqaVdSkREREREsigVbhER\nERERK1LhFhEREckGfvzxR4oWLYqbmxt79+6lUqVKbNq0yexY2YIKt4iIiIhJvL29yZUrF66uruTN\nm5dWrVpx/vz5Rx4fGhqKnZ0dL7300n3P79u3Dzs7Oxo1avTIcwcOHMjMmTO5desWVatW5eDBgzRo\n0ACAkJAQunTpkj6flDxEhVtEJBnBwWYnEJHswDAMVqxYQWRkJBcvXsTT05OgoKDHnpM/f3527NhB\nRERE0nPz58+nTJkyGIaR7DkWi4W//vqLChUqpGt+SR0VbhGRZGhbQBHJaI6Ojrz88sscPnz4scc5\nODjw4osvsnjxYgASEhJYunQpnTt3TnY3jdjYWFxdXUlISKBKlSqULl0auDu7vn79elavXs24ceNY\nsmQJrq6u+Pr6pv8nl82pcNsww7j7EBERkczrXkmOjo5myZIl+Pn5pXhOly5dWLBgAQBr1qyhUqVK\nFC5cONljHR0duX37NgD79+/nxIkTwN3ZdcMwaNGiBR9++CEdO3YkMjKSsLAwAMaPH0/r1q3T/PkJ\n2JsdQERERCS7slgsvPjii9jb2xMVFUWBAgVYvXp1iuf5+fkRERHB8ePHWbBgAW+88QbR0dFpyvHg\n7PiQIUOe+v3kfprhFhERETGJYRgsW7aM69evExsby/Tp02nYsCHh4eH89ddfuLq64urqipub20Pn\ndunShenTpxMaGkq7du10g59MTIXb1jlfMzuBiIiIpIJhGLRr144cOXKwdetWihUrRmRkJJGRkdy6\ndeuh419//XVmzZrFCy+8gJOTU5rHFutR4bYh99ZsJ63ddrgN71TBaDYEwyFKa7pFnoAumhSRjHJv\nZtpisSTNdpcvXz7F84oXL86mTZsYM2ZMmjMULFiQM2fOaJbcSlS4bdmd3PDZHnA7D4EVocwKsxOJ\nZBkjR5qdQESyi9atW+Pq6oq7uzvDhw9nwYIFjy3c/56NrlOnDgULFkx6/nEz1Y97rX379gDky5eP\nGjVqADB27Fhatmz5RJ+LJM+wZMF/yhiGoX+BpcK9nyuLBdb/uZ53V71LxfwVmdpiKkXdi5obTiST\nM4y7PzsikrWpM0hqPep7JT2+h0yd4X7zzTfx9PTEx8cn6bmIiAiaNWtGmTJlaN68OTdu3DAxoe1o\nUqIJ+3vvx7egL76f+TJp2yTiEuLMjiUiIiJi80wt3N27d39o65vx48fTrFkzjh8/TpMmTRg/frxJ\n6WyPo70jwxsOZ0fPHaw5tYYac2qw/dx2s2OJiIiI2DTTl5ScOXOG1q1bc+DAAQDKlSvHxo0b8fT0\n5NKlS/j7+3P06NH7ztGvh1Ln30tKHmSxWFh6aCn91/anVelWjGs6jrzOeTM2oEgmpiUlIrZBnUFS\ny2aXlCQnPDwcT09PADw9PQkPDzc5UdZlsTy6MBiGQYdKHTj87mEccjhQcWZFvt73tf5SEvn/goPN\nTiAiIrYiU99p8nFX24b8a88uf39//P39MyaUjXF3cmd6y+m8UfUN3l7xNvP2zmPWC7Mo90w5s6OJ\nmErbAoqIZE+hoaGEhoam63tmyiUloaGhFCxYkIsXL9KoUSMtKckgCYkJzPx9Jh9v+pje1XvzYf0P\ncc7pbHYsERGRp6bOIKmVrZaUtGnThvnz5wMwf/58XnzxRZMTZR857HIQVDuIfb33cTziOJVmVWL1\nydUpnygiIiJiJRs2bODDDz8kNjbW7ChPzdQZ7k6dOrFx40auXr2Kp6cnH3/8MW3btuXVV1/lr7/+\nwtvbm6VLl5InT577ztO/VjPG6pOrCVwVSPVC1fm/Fv9HYdfCZkcSERF5IuoMWd+IESN4+eWXOXLk\nCB07drTaODY7w71o0SL+/vtv7ty5w7lz5+jevTt58+Zl3bp1HD9+nLVr1z5UtiXjtCjVgoPvHKRs\nvrJUmV2FaTunkZCYYHYsEREReUonTpzAycmJLl26PPKY119/nUKFCuHm5kaJEiUeunV8Sq8/bpwZ\nM2ZQo0YNnJyc6N69e6oy582bl2XLllG9evVUHZ8ZZbolJZK5OOd0ZlTjUWzuvpkfj/5IrS9q8fuF\n382OJWJ1umhSRMwSHR1NvXr1rPLegYGB1KpV67G3eR86dCinT5/m1q1b/PLLL0yfPv2++6ak9Prj\nxilSpAjDhw/nzTffTHXmNm3acOvWLUqXLp3qczIbFW5JlXLPlGND1w30rd2XNovb0GdVH27G3DQ7\nlojVjBxpdgIRya6mT5/O9u3bSUhI398qL168GA8PD5o0afLYJRIVK1bEyckp6WN7e3sKFCiQ6tcf\nN067du1o27Yt+fLlS3Xuzz//nDlz5nDzZtbtHSrckmqGYdClShcOvXuIuIQ4yn9ansUHF2ttnIiI\nSDoJCwujTJkyODg4cPHixRSPb9WqFR4eHsk+2rRpk3TcrVu3CA4OZsqUKan6//a7776Li4sLFStW\n5KOPPqJatWqpej2146S2O9wr2b6+vsyaNStV52RGKtzyxPI65+Wz1p/x/avfM3bzWAK+CeDEtRNm\nxxIREcnS4uPjWbp0Ke3ataNgwYJcuHAhxXNWrFjB9evXk30sX7486bjhw4fTs2dPChcu/NjlJPfM\nnDmT27dvs27dOj766CN27dqVqtdTO05qMgDMnTuXrl270rt3b6ZNm8adO3dSdV5mo8ItT61O0Trs\neWsPASUD8Jvrx8cbPyY2Putu2SMiImKmTz/9lJ49ewKkunCnxt69e1m/fj19+/YFUj+7bBgG/v7+\ntG/fnkWLFqX4+pOMk5oMCQkJ7N+/nwoVKvDSSy+RkJDAN998k6rsmU2mvtOkZH45c+RkQJ0BtK/Y\nnvdXv0/l2ZWZ2XImTUo0MTuaiIhIlnHq1Cl27dpFnjx52LJlC/Hx8fz9998AREREMGfOHAoUKEDl\nypXv263j+eefZ8uWLcm+Z4MGDVi5ciWhoaGcOXOGYsWKAXD79m0SEhI4cuQIu3fvTjFbXFzcY9dc\n33t948aNqR4nNTPcP/zwAxcuXKBTp04AFCpUiEmTJj3RBZeZhQq3pIti7sX4scOPLD+2nB7Le1C3\nWF0mNZ9EwdwFzY4m8lSCg81OICLZhcVi4auvvuLrr7/Gzu7u4oOwsLCkGe6vvvqKRo0aUa1aNd54\n4w0WLlyYdO4vv/yS4vu//fbbvPbaa0ljffLJJ5w5c4bZs2c/dOyVK1dYv349rVu3xsnJiXXr1vHd\nd9+xbt26FF+vVKlSUjl+1DgJCQnExcURHx9PQkICsbGx2NvbkyNHjoe+JkuWLGH16tVJr509e5ZS\npUqxfPnypPXp3bp1wzAMvvzyy9R9sU2iJSWSrtqUbcOhdw/h5eqFzywfpu+cTnxivNmxRJ6YtgUU\nkYywY8cOWrduzcmTJ0lMTARgy5Yt7N+/nw0bNhAaGsrp06cpVKgQ9vb2REREPPEYzs7OFChQgAIF\nCuDp6Unu3LlxdnZOmrVu2bIl48ePB+7OPM+ePRsvLy/y5cvH8OHD+frrr6lZs2aKr6c0DsCoUaPI\nlSsXEyZM4JtvvsHZ2TnZfbx79OjBjh07WL9+fdJzoaGhODg48M4777Bz504Azp07Z7UtFNOTqXea\nfFq6a1TWcPjKYfqs6kPEPxF82vJT6hara3YkERHJZmyhMwQGBjJs2DAKFy5My5YtWbVqldmRMoU7\nd+7g6+vL/v37H5ohfxo2e6dJsW0V8ldgfdf1DKk3hA7fd6DbT90Ivx1udiwREZEspWzZsoSHhxMT\nE4Obm5vZcTINBwcHDh06lC5l29o0wy0ZIjI2kpEbRzJ/33yCGwbTu0Zv7O10CYGIiFiXLXSGa9eu\nMW/ePNzd3fHx8cHPz8/sSDbJmjPcKtySoQ5dPkSfX/pwI+YGn7b8lDpF65gdSUREbJg6g6SWCvcD\n9MOTtVksFhYfXMygXwfRrGQzJjSdQAGXAimfKJKBQkJ04aSILUhNZzBGpu4mLI9jCVYvyepUuB+g\nwm0bbsXe4uONHzN/33xCGobQu0Zvcthl/nVYkj0YBuivGZGsT51BUkuF+wH64bEtBy8fpM+qPtyM\nvallJpJpqHCL2AZ1BkktFe4H6IfH9lgsFhYdXMSgXwcRUDKA8U3Ha5mJmEqFW8Q2qDNIamlbQLF5\nhmHwms9rHAk8goeTB5VmVmLm7zNJSEwwO5qIiIhImqhwS6bi5ujGpIBJbHhjA0sPLaXmnJpsP7fd\n7FgiIiI24c6dO/To0QNvb2/c3Nzw9fVl9erVqX79nsWLF1O+fHly585NqVKl2LJly2PHPXHiBE5O\nTnTp0iXFjIsXL6Z27drkzp0bT09PnnvuOWbNmpX0ure3N7ly5cLV1ZWCBQvSvXt3oqKinuCrkPFU\nuCVTqlSgEr+98RsD6wzkle9eoceyHlyJumJ2LMlGgoPNTiAi2cmePXt44YUXaNy4MXPnzmXu3Ln0\n7duXkHTeLik+Pp5ixYqxadMmbt26xejRo3n11Vc5e/Zsql4H+PXXXxkyZAjz58/n9u3bbN68mRIl\nSjx23MDAQGrVqoVhPH5HmEmTJtG3b18GDx5MeHg44eHhzJ49m61btxIXFwfc/a34ihUriIyM5I8/\n/mD37t2MHj06jV8ZK7NkQVk0tjylmzE3Lf1W97Pkn5jf8umuTy3xCfEPHXN3ta0J4UREJFPLSp2h\nQ4cOlpUrVyZ9HBkZabGzs7NERUVZddzKlStbfvjhh1S/7ufnZ5k3b16q33/RokWWV1991RISEmJ5\n/fXXH3ncjRs3LC4uLo/NYrFYLN7e3pb169cnfTxw4EBLq1atLBaLxTJu3DhLyZIlLa6urpYKFSpY\nfvzxx1TnfNT3Snp8D2mGWzI9N0c3JgdMZsMbG1hyaAm1vqjFjvM7zI4lIiKSrnbs2EHjxo2Bu5sJ\njBs3jsDAQHLlyvXIc1q1aoWHh0eyjzZt2qQ4Znh4OMePH6dixYqpej0hIYE9e/Zw+fJlSpcuTdGi\nRQkKCiImJibZ82/dukVwcDBTpkxJ8cLD7du3ExsbS9u2bVPMfe+9zp07xy+//EK1atUAkpa33Bv3\n9ddf59KlSym+n7WpcEuWUalAJULfCKX/c/15aclL9FzeU8tMRETEJhw6dIh8+fKxceNGVq9eTZ8+\nffD29mbatGmPPW/FihVcv3492cfy5csfe25cXBydO3emW7dulClTJlWvh4eHExcXx3//+1+2bNnC\n3r17CQsLe+SSjuHDh9OzZ08KFy6c4nKSq1ev8swzz2Bn9796WqdOHTw8PMiVK1fSOnGLxcKLL76I\nh4cH9evXx9/fnw8//BCAV155hYIFCwLw6quvUrp0aXbt2vXYcTOCCrdkKYZh0LlyZ44EHsHVwZWK\nMysy6/dZYGg3ExERybp+++03Xn75ZQICAmjRogVTpkxh/PjxnDx50irjJSYm0qVLF5ycnJgxY0aq\nX3d2dgYgKCgIT09P8uXLR//+/Vm1atVD77F3717Wr19P3759AVKc4c6XLx9Xr14lMTEx6blt27Zx\n/fp18uXLl/S8YRgsW7aM69evc+bMGWbMmIGjoyMACxYswNfXN2mW/+DBg1y7du0Jvzrpz97sACJP\nI4+zOzAFCrzJu0cCodcXsHImhlH7vuO09aqIiGQFGzduJCgoKOljBwcHXF1dk2a+P//8cwoUKEDl\nypWpXr160nHPP//8I3cIadCgAStXrnzoeYvFQo8ePbhy5QqrVq0iR44cqX7dw8MDLy+vVH9OZ86c\noVixYgDcvn2bhIQEjhw5wu7dux863s/PD0dHR3766SdeeumlVI3xb2fPnuWtt95iw4YN+Pn5YRgG\nvr6+mWIfdhVuydou+8CXG6HyQujYDk60hHXjIDq/2ckkiwsJufsQEbE2i8XCtm3b+Prrr5OeW7ly\nJTdv3qRp06Z89tlnNGrUiGrVqvHGG2+wcOHCpON++eWXJx7vnXfe4ejRo6xbty5pZvhJXu/evTvT\np0+nRYsW2NvbM2XKFFq3bv3QcW+99RadOnVK+hw/+eQTzpw5w+zZs5PNlSdPHoKDg3n33XexWCw0\nb94cFxcX9u/fn6pt/6KiojAMg2eeeYbExEQWLFjAwYMH7zumW7duGIbBl19+meL7pScVbsmS7v/H\nqoFhvA7HWtP3pxAW1qtIiH8Ib1d/G8jxiHcQebyRI1W4RcT6wsLCWLp0KfHx8cydOxeAa9eucfr0\naTZv3oyLiwunT5+mffv22NvbExERkabxzp49y+eff46Tk1PSWmeAzz//nE6dOqX4Otxdl3316lXK\nlCmDk5MTHTp0YNiwYQC0bNmSBg0aMGTIEJydnZOWoADkzp0bZ2dn8uXL98h8gwYNokiRIkycOJGu\nXbvi4uJCiRIlmDhxInXq1Hns51ahQgUGDBiAn58fdnZ2dO3alXr16t13zPnz55M+j4ykW7uLTbh3\nHYbFAvvD9/P+6veJ+CeCaS2m0dC7obnhJEvSrd1FbIMtdIbAwECGDRtG4cKFadmyZbLrpSVld+7c\nwdfXl/379z+0jAZ0a3eRJ1LZszIbum7go/of0fWnrnT8viPnbp4zO5aIiMhTKVu2LOHh4cTExODm\n5mZ2nCzLwcGBQ4cOJVu2rU0z3GIT/j3D/W/RcdFM3DqR6bum07d2XwbWGYhzTueH30DkAZrhFrEN\nttAZrl27xrx583B3d8fHxwc/Pz+zI9kka85wq3BLtnDmxhkGrh3Inot7mNR8Eu3KtUtxP1DJ3lS4\nRWyDOoOklgr3A/TDI09rw+kNvPfLexTMXZCpLaZSsUDyd9YS0S4lIrZBnUFSS4X7AfrhkbSIT4xn\n1u+z+HjTx7xW6TVC/EPwcPYwO5aIiFiBOoOkli6aFElH9nb2BNUO4vC7h4lNiKX8p+WZs2cOCYm6\nW6WIiIikP81wS7YXdjGM91a/R9SdKKY/P526xeqaHUlERNKJOoOklpaUPEA/PJLeLBYLiw8u5oN1\nH9Dw2YZMaDqBIm5FzI4lIiJppM4gqaUlJSJWZhgGnXw6cSTwCN55vKkyuwrjNo8jJj7G7GgiIiKs\nXLmS6Ohos2PIU1LhFvmX3A65Gd14NDt77mTnhZ1UmlmJ5ceWa3YkG9IOJSKSmQwYMICEhLRfa/T6\n669TqFAh3NzcKFGiBGPGjHnksefPn6d169bky5ePQoUKERQUlJTB398fZ2dnXF1dcXV1pXz58g+d\nv3jxYsqXL0/u3LkpVaoUW7ZsSXP+rEpLSkQeY+2ptby/+n2edX+W/2vxf5R7ppzZkSSDaB9uEduQ\nVTrDnj17GDFiBP/88w+dO3cG4MCBA+TJk4fOnTszbtw45s2bl+ZxDh06RMmSJXFycuLYsWM0bNiQ\nr776ihYtWjx07EsvvYS7uzufffYZ169fp1mzZvTq1YugoCAaNWpEly5dePPNN5Md59dff6VXr14s\nXbqUWrVqcfHiRSwWC4ULF07z52At1lxSYp+ms0VsXPOSzdnfez8zds2g3rx6vFHlDUY0HIG7k7vZ\n0URExIZUr14dV1dXAgMDadmyJQC3b9/G3d2d/Pnz07t373QZp2LF++8/YW9vT4ECBZI99tChQ0yd\nOhUHBwc8PT0JCAjg8OHDSa8/roQGBwcTHBxMrVq1AChUqFA6pM+6tKREJAU5c+Skn18/Dr17iJux\nNyn3aTnmhc0j0ZJodjQREbEhO3bsoHHjxsDdMjtu3DgCAwOJjo5OKq6P0qpVKzw8PJJ9tGnT5r5j\n3333XVxcXKhYsSIfffQR1apVS/Y9AwIC+Pbbb/nnn3+4cOECq1evvm8mfOjQoeTPn5969eqxcePG\npOcTEhLYs2cPly9fpnTp0hQtWpSgoCBiYrLvdVFaUiLyhH6/8DvvrX6PhMQEpj0/jee8njM7kliB\nlpSI2Ias0hkOHTpE165dGTt2LBaLhZ9//pmqVavSq1cvq4xnsVjYuHEjr7zyCqtWrUq20EdERNC0\naVMOHDhAQkIC3bp1S1rWsmvXLipWrIiDgwOLFi2iT58+7N27lxIlSvD333/j5eVFjRo1+Pnnn7G3\nt6dt27b4+/szevRoq3w+6UG7lIhkIjWL1GTrm1sJqhXEy0tfpttP3bgYedHsWCIikoX99ttvvPzy\nywQEBNCiRQumTJnC+PHjOXnypFXGMwwDf39/2rdvz6JFix563WKxEBAQQPv27YmOjubq1atEREQw\nePBgAGrVqoWLiws5c+aka9eu1K1bl1WrVgHg7OwMQFBQEJ6enuTLl4/+/fsnvZ4dqXCLPAU7w44u\nVbpwNPAoBXMXxGeWD//Z+h/uJNwxO5qkk+BgsxOISHayceNG6tWrl/Sxg4MDrq6uHDp0KOm5iIgI\nJkyYwJdffsmePXvuO//5559P2jHkwccLL7zwyHHj4uJwcXF56PmrV6+yZ88e+vTpQ86cOcmbNy/d\nunVLVWlvWnI/AAAgAElEQVT28PDAy8srNZ92tqHCLZIGro6ujG86nu09trPx7EYqzazEiuMrssSv\nL+XxtC2giGQUi8XCtm3b7lvWsXLlSm7evEnTpk2Tnps/f37S7iCTJ0++7z1++eUXIiMjk32sXLkS\ngCtXrrB48WKioqJISEhgzZo1fPfdd7Rt2/ahTM888wyFChVi1qxZJCQkcOPGDebPn0+VKlW4efMm\na9asISYmhvj4eBYuXMjmzZvvW9/dvXt3pk+fzpUrV7h+/TpTpkyhdevW6f2lyzK0S4lIOiidrzQr\nXlvBqhOrGLB2ANN2TmNywGQqFahkdjQREcnEwsLCWLp0KfHx8cydOxeAa9eucfr0aTZv3nzf7POf\nf/7JK6+8gr29PREREU88lmEYzJ49m3feeQeLxUKZMmX4+uuvqVmzJgAtW7akQYMGDBkyBMMw+OGH\nHxg4cCDjxo3D3t6eJk2aMGXKFOLi4hg+fDhHjx4lR44clC9fnmXLllGqVKmksYYPH87Vq1cpU6YM\nTk5OdOjQgWHDhqXxq5V16aJJkXQWlxDH7N2zGbVpFC+Xf5mPG31Mfpf8ZscSEcmWbKkzBAYGMmzY\nMAoXLkzLli2z9Zpoa9BFkyJZSM4cOQmqHcTRPkdxtHekwswKfLLtE2LjY82OJiIiWVjZsmUJDw8n\nJiYGNzc3s+PIE8i0M9ze3t64ubmRI0cOcubMya5du5Jes6V/rYrtO3b1GAN/HciRK0f4pPkntC3b\nFsMwzI4lIpIt2FJnuHbtGvPmzcPd3R0fHx/8/PzMjmRTrDnDnWkLd/HixdmzZw958+Z96DVb+uGR\n7GPtqbX0X9OfAi4FmBIwhSoFq5gdSR4jJEQXTorYAnUGSa1su6REPyBiS5qXbM7e3nt5teKrBHwT\nQK/lvQi/HW52LHmEkSPNTiAiIrYi0xZuwzBo2rQpNWrUYM6cOWbHEUkX9nb29K7Rm6N9juLu5E7F\nmRUZv2U8MfHZ93a3IiIiti7TLim5ePEihQoV4sqVKzRr1ozp06dTv3594G4ZD/7XXSn8/f3x9/c3\nKanI0ztx7QQfrPuAfZf2MbHZRF4u/7LWd2cSurW7iG3QkhJJrXvfK6GhoYSGhiY9P3LkSNtdw/1v\nI0eOJHfu3AwYMADQD4/Yng2nN9BvTT/cHd2ZEjCF6oWrmx0p21PhFrEN6gySWtluDXd0dDSRkZEA\nREVFsXbtWnx8fExOJWI9jYs35o+3/qBrla60WtSK7su683fk3w8dZxh3HyIiIpJ1ZMrCHR4eTv36\n9alatSq1a9emVatWNG/e3OxYIlaVwy4HPav15FifY3i6eOIzy4fRm0bzT9w/ZkfLlv61ak1ERCRN\nssSSkgfp10OSHfx5/U8GrxvMrgu7mNB0Ah0qdsDO7u70tr79RURSR51BUitb7sP9OPrhkexk09lN\n9F3dF+eczmwbMQUu1FLhFhFJJXUGSS0V7gfoh0eyi6T12kYiVFkAjYfB6cawfhzc8ko6Tj8OIiLJ\nU2eQ1FLhfoB+eCS7eOgCSYfbUG881JgFu4Jg6yCIc1HhFhF5hLx583L9+nWzY0gW4OHhQURExEPP\nq3CLZDP3CviZ62cZsn4IW/7awrgm43jN5zXsjEx5DbSIiEiWpsItks3cK9z3vv23/rWVfmv6YRgG\nUwKmUKdoHfPC2ZiQkLsPERHJ3lS4RbKZBws3QKIlkW8PfMvQ9UPx8/JjXJNxlMxb0pyANkQ3vhER\nEbDhG9+ISOrZGXa8Xvl1jvU5RhXPKtT+ojb91/Qn4p+H16GJiIhIxlPhFslCLJZHz7rmypmLYQ2G\ncejdQ8TEx1BuRjkmb59MbHxsxoYUERGR+2hJiYiNOnLlCIPXDebg5YOMbzqe9hXaY+i+8KmmJSUi\nIgJaw212DJEs4bfTvzHw14HktMvJJ80/oV6xemZHyhJUuEVEBLSGW0RSoVHxRvze63f61OpD5x86\n8/LSlzlx7YTZsTK94GCzE4iIiK3QDLdINvJP3D9M2zmN/2z7D519OjO84XCeyfWM2bFEREQyLc1w\ni8gTcc7pzOB6gzkSeIRESyLlPy3Pf7b+h5j4GLOjiYiI2CzNcItkY8euHmPI+iGEXQxjbJOxdKzU\nUXesFBER+RddNCki6WLz2c0MWDsACxY+afYJDb0bmh1JREQkU1DhFpF0k2hJZOmhpQxdP5QqnlWY\n0HQCZZ8pa3YsERERU2kNt4ikGzvDjo6VOnIk8Aj1itWj3pf1CFwZyOWoy2ZHM0VIiNkJRETEVmiG\nW0SSdS36GqM2jeKb/d8wwG8AfZ/ri3NOZ7NjZRjtwy0iIqAlJWbHEMkWTkacZOj6oey6sIvRjUbT\nuXLnbHFhpQq3iIiACrfZMUSylW3ntjFg7QDuJNzhk2af0Kh4I7MjWZUKt4iIgAq32TFEsh2LxcL3\nh79nyPohVMhfgYlNJ1I+f3mzY1mFCreIiIAumhSRDGYYBu0rtufwu4dp7N2Yhl81pPeK3oTfDjc7\nmoiISKalwi0iT8zR3pF+fv041ucYLjldqDizIiGhIUTGRpodLd0EB5udQEREbIWWlIhImp25cYYR\nv43g1z9/ZVj9YbxV/S0ccjiYHUtERCTNtIZbRDKVfZf2MXT9UI5dO8aYxmN4teKr2WJHExERsV0q\n3CKSKf12+jcGrxtMoiWRCU0n0KREE7MjiYiIPBUVbhHJtO7taPLhhg8p4VGC8U3G41vI1+xYIiIi\nT0SFW0QyvbiEOL744ws+3vQxjYs3ZlSjUZTwKGF2LBERkVTRtoAikunlzJGTd2q+w4mgE5TLV45a\nc2rx/i/vcyXqitnRHiskxOwEIiJiKzTDLSIZ6nLUZUZvGs23B77l/drv08+vH7kdcj90nGHc/a9Z\nP+q68Y2IiIBmuEUkCyrgUoBpz09jZ8+dHLl6hDLTyzDr91nEJcSZHU1ERMQqNMMtIqb64+IfDFk3\nhDM3zjC2yVheLv8yhmFohltERDIFXTQpIjbj11O/MnjdYHLmyMmEphNoVNwfUOEWERFzqXCLiM0w\nDMBIhIpLoMkwuFoO1o2D8Cr3HZdRP/oq3CIiAlrDLSK2xmIHBzvBjKNw4nnoEgDtukKeMxkeJTg4\nw4cUEREbpRluEcmUDANwvMWIXyYx4/cZvFHlDT6s/yHP5HrG7GgiIpKNaIZbRGxbrBsjG43k0LuH\niI2PpdyMcozdPJbouGizk4mIiKSaCreIZHoFcxfk0xc+ZXuP7ey9tJfS00szZ88c4hPjzY4mIiKS\nIi0pEZFM6XHbAu66sIvB6wZzMfIiYxqP4aXyL2HcO0FERCQdaZcSEcm2LBYLa06tYej6oeQwcjC2\nyVialWim4i0iIulKhVtEsr1ESyLfH/6e4b8Np7BrYcY2HotfUb80v29IyN2HiIhkbyrcIiL/X3xi\nPPP3zmfkxpFULViV0Y1HU9mz8lO/n/bhFhERUOE2O4aIZEIx8THM3j2b8VvG06REE0b6j6RU3lJP\n/D4q3CIiAtoWUETkIU72TvR9ri8ngk5QLl85nvviOXqv6M2FWxfMjiYiItlUqgp3VFQUR48e5dix\nY0RFRVk7k4hImrk6ujK84XCO9TmGm6MbPrN8GLR2ENeir5kdTUREsplHLimJjIxkzpw5LF68mKtX\nr+Lp6YnFYiE8PJx8+fLRuXNnevXqRe7cuTM6s5aUiMgT+zvyb0ZtHMV3h7/jvdrv0e+5frg6uj7y\neC0pERERsPKSkhdffBFXV1d+/vln/vzzT7Zv386OHTs4ffo0K1aswMXFhbZt26ZpcBGRjFLYtTCz\nWs1iZ8+dHL92nFLTSzFl+xRi4mOSPT44OIMDioiIzcqUF02uXr2avn37kpCQQM+ePRk8ePB9r2uG\nW0TS6kD4AT767SPCLoYxouEIulXthr2dvdmxREQkk8mwXUoiIiI4ceIEsbGxSc81aNAgTQM/SkJC\nAmXLlmXdunUUKVKEmjVrsmjRIsqXL/+/0CrcIpJOtp/bzocbPuTCrQuMajSK9hXbY2foenIREbkr\nQwr3nDlzmDZtGufPn6dq1ars2LEDPz8/NmzYkKaBH2X79u2MHDmS1atXAzB+/HgAhgwZ8r/QKtwi\nko4sFgvrT6/nw/UfEpcYx5jGY3i+1PO6a6WIiGTMtoBTp05l165dPPvss/z222+EhYXh7u6epkEf\n58KFCxQtWjTpYy8vLy5c0HZeImI9hmHQtERTdvbcyYgGIxj06yAafNWAzWc3mx1NRERsQIqF28nJ\nCWdnZwBiYmIoV64cx44ds1qg1M4oGUbIvx6hGMajb8McEnJ3x4EHHzpex+t4Hf/v4+3sDF6q0I7D\nffbzzNledP2pKy0XtiTsYliWyK/jdbyO1/E6Pu3Hd+sWSkhISNIjPaS4pKRdu3bMmzePqVOnsn79\nejw8PIiPj2fVqlXpEuBBO3bsICQkJGlJybhx47Czs7vvwkktKRERawsJgaEfxTLnjzmM2TyGBs82\nYFSjUZTJV8bsaCIikoEyZA33lStXyJ8/PwChoaHcunWLFi1a4ODgkKaBHyU+Pp6yZcuyfv16Chcu\nTK1atXTRpIhkOMP43z7cUXeimL5rOpO2T6Jt2baMaDiCYu7FzA0oIiIZwqpruH/++Wfy589P5cqV\n8fLyYuvWrfj7+9OmTRurlW0Ae3t7ZsyYQUBAABUqVKBDhw73lW0RkYzm4uDCkHpDON7nOJ4unvh+\n5ku/1f24HHXZ7GgiIpIFPHKG28fHh++++45y5cqxc+dOBg0axKZNmzI6X7I0wy0i1vbvGe4HXbp9\nibGbx7LwwEJ6VevFoDqDyJcrX8YGFBGRDGHVGW57e3vKlSsHQO3atYmMjEzTQCIitqJg7oJMe34a\ne9/ey42YG5SdUZYRv43gRswNs6OJiEgm9MgZbi8vL/r375/U6KdMmZL0sWEY9O/fP0OD/ptmuEXE\n2h43w/2g09dPM2rTKH4+/jPv136f92u/j6ujq3UDiohIhrDqDHfPnj2JjIzk9u3b3L59O+njew8R\nEVsWHJz6Y4t7FGde23lse3MbR68epdT0UkzcOpGoO1HWCygiIllGqm7tntlohltEMrPDVw4TEhrC\n5r8280GdD+hdozfOOZ2TPdYw7v5Xf6WJiGROVt0W8ODBg5w6dYq2bdsC0LdvX27evIlhGPTp04dq\n1aqlaeC0UOEWkaxg36V9BIcGs/vv3QytN5Se1XriaO943zEq3CIimZtVl5QMGTKEZ555JunjtWvX\n0qpVK/z9/fn444/TNKiISHZQpWAVfur4Ez91/ImVJ1ZSZkYZ5uyZQ1xCnNnRREQkAz1yhrt69ers\n2bMn6ePatWuzc+dOAOrWrcvWrVszJmEyNMMtIlnR9nPbGRE6gj+v/8mIBiPoXLkzOXPYA5rhFhHJ\nrKy6pKRMmTIcP3482ZNKly7NiRMn0jRwWqhwi0hWZRjAsxuh8XBwCYfQEDjUASz3/8JRf8WJiGQO\nVl1SUrhwYXbs2PHQ89u3b6dIkSJpGlREJLMLCbHim59tCF9uhFUz4Lmp8E5lKP9fMBKtOKiIiJjl\nkTPcu3btokOHDnTr1o1q1aphsVj4448/+Oqrr1iyZAm1a9fO6KxJNMMtItb2JPtwp20cC5ReRbWB\nI0i0JPKx/8e0KtMK497VlCIiYiqrLikBCA8PZ8aMGRw+fBiAihUrEhgYiKenZ5oGTSsVbhGxtowr\n3Hf/m5hoYdmxZYz4bQTOOZ352P9jmpdsruItImIyqxfuzEqFW0SsLaML972xEi2JfH/4e4JDg8nn\nnI9RjUbRqHgj6wcREZFkWXUN9wsvvMB3331HdHT0Q69FRUWxZMkSWrZsmabBRUTkfnaGHa9WfJWD\n7xykd43evLXiLRrPb8yWv7aYHU1ERJ7SI2e4L1++zIwZM/j+++/JkSMHhQoVwmKxcOnSJeLj4+nQ\noQOBgYHkz58/ozNrhltErC6jZrhTEp8Yz4J9Cxi1aRRl8pVhVKNR1CpSy+xYIiLZRoYtKbl06RJn\nz54F4Nlnn6VgwYJpGjStVLhFxNpCQqy8U8kTupNwh3lh8xizeQxVC1YluGEwNQrXMDuWiIjN0xpu\nEZFsJiY+hjl75jBh6wR8C/mqeIuIWJkKt4hINhUTH8PcP+Yybsu4pBnvmkVqmh1LRMTmqHCLiGRz\nMfExzAubx7gt46jsWZnghsFa4y0iko4yrHBHR0dz7tw5ypYtm6bB0osKt4jI/WLjY5kXNo+xW8bi\nU8CH4IbB1PYy7wZlIiK2wqrbAt6zfPlyfH19CQgIACAsLIw2bdqkaVAREUlfjvaOvFPzHU4GnaRN\n2Ta0/649zy98nh3nd5gdTUQk20uxcIeEhLBz5048PDwA8PX15c8//7R6MBERM2WmHUqehKO9I71r\n9OZE0AleLPsiHb7vQItvWrD93Hazo4mIZFspFu6cOXOSJ0+e+0+yS/E0EZEsbeRIsxOkjaO9I2/X\neJvjfY7Trlw7Ov63IwHfBKh4i4iYIMXmXLFiRRYuXEh8fDwnTpwgKCiIOnXqZEQ2ERFJo3vF+0TQ\nCV4u/zKd/tuJ5l83Z9u5bWZHExHJNlK8aDI6OprRo0ezdu1aAAICAhg+fDhOTk4ZEjA5umhSRKwt\ns9xpMr3dSbjD/L3zGbtlLKXyliKkYQh1i9U1O5aISKZl9V1K4uPjadasGb/99luaBklvKtwiYm22\nWrjvuZNwhwX7FjBm8xhK5S1FcMNg6hWrZ3YsEZFMx+q7lNjb22NnZ8eNGzfSNIiIiGQuDjkc6Fmt\nJ8f7HKdjxY50/bErTRc0ZfPZzWZHExGxOfYpHeDi4oKPjw/NmjXDxcUFuNv0p02bZvVwIiJmCQ42\nO0HGyJkjJz2q9aBrla58vf9rui3rhnceb4IbBtPg2QZmxxMRsQkpruH+6quvHj7JMHjjjTeslSlF\nWlIiImIdcQlxfLP/G0ZvHs2z7s8S3DCYht4NzY4lImIa3dpdRESsIi4hjoUHFjJ602iKuhclpGGI\nireIZEsZUriLFy+e7MBm3vxGhVtEJGPEJ8azcP9CRm8eTRHXIgQ3DMbf2x/DMMyOJiKSITKkcF+9\nejXpzzExMXz//fdcu3aNUaNGpWngtFDhFhHJWPGJ8Xx74FvGbB5D/lz5+ajBRwSUDFDxFhGbZ9qS\nkmrVqvHHH3+kaeC0UOEWETFHQmIC3x/+ntGbR+OYw5GPGnxEm7JtsDN0B2IRsU0ZUrj37NmTNIOR\nmJjI7t27mTVrFvv27UvTwGmhwi0i1hYScvchyUu0JLL82HJGbxpNbEIsw+oPo32F9uSwy2F2NBGR\ndJUhhdvf/39r9ezt7fH29mbgwIGULVs2TQOnhQq3iFibrd/4Jr1YLBbWnFrDqE2juBp9laH1htLZ\npzM5c+S877h7K0/0NRWRrEa7lIiIWIkK95OxWCyEngll9ObR/Hn9T4bUHUK3qt1wtHcEVLhFJOuy\n+p0mAaZOncqtW7ewWCz06NGDatWqsWbNmjQNKiIitsUwDBoVb8T6rutZ+NJClh9fTslpJZm6YyrR\ncdFmxxMRMVWKhXvu3Lm4ubmxdu1aIiIiWLBgAUOGDMmIbCIikgXVKVqHla+tZFnHZWw8u5ESU0tA\n3YngEGl2NBERU6R4a/d7U+grV66kS5cuVKpUyeqhREQka7u7hKQ68AMUOAj1x8L7JTD8g2BnEMR4\nAFpiIiLZQ4oz3NWrV6d58+asWrWKgIAAbt26hZ2dtn8SEdsWHGx2AhtyuRL891uYtxXynIb3SkGT\nDyHXFbOTiYhkiBQvmkxMTCQsLIySJUuSJ08erl27xoULF6hcuXJGZXyILpoUEcla/n3R5JkbZ5iw\nZQJLDi2hW9VuDKwzkMKuhc0NKCLyCBly0eT27dspW7YsefLk4euvv2b06NG4u7unaVAREcm+vPN4\nM6vVLA68cwCLxUKlmZUIXBnI2RtnzY4mImIVKRbu3r174+Liwr59+5g8eTKlSpWia9euGZFNRERs\nWBG3IkxpMYWjfY7i5uhGtc+r0WNZD05GnDQ7mohIukqxcNvb22MYBj/99BOBgYEEBgYSGakrzUVE\nJH0UcCnAuKbjOBF0gqLuRfGb68frP7zO4SuHzY4mIpIuUizcrq6ujB07lm+++YZWrVqRkJBAXFxc\nRmQTEREbYbGkvCNJXue8hPiHcOq9U1QqUIlG8xvxytJXCLsYljEhRUSsJMXCvWTJEpycnJg3bx4F\nCxbkwoULDBo0KCOyiYiYJiTE7ATZl5ujG0PqDeHP9/6kbtG6tFrUilbftmLbuW1mRxMReSqpurX7\nmTNnOHnyJE2bNiU6Opr4+Hjc3NwyIl+ytEuJiFibbu2eecTEx/Bl2JdM3DaRYu7FGFpvKAElAzDu\nbX0iImJFGbJLyeeff0779u15++23ATh//jzt2rVL06CPExISgpeXF76+vvj6+rJ69WqrjSUiIpmf\nk70T79R8hxNBJ+hVrRcD1w6k+ufVWXpoKQmJCWbHExFJUYoz3FWqVGHXrl0899xzhIXdXUfn4+PD\ngQMHrBJo5MiRuLq60r9//0ceoxluEbE2zXBnXomWRFYeX8m4LeO4Gn2VD+p+QJfKXXC0dzQ7mojY\noAyZ4XZ0dMTR8X9/icXHx1v913gq0yIi8ih2hh2ty7Zm65tbmdN6Dt8f/p6S00oyeftkbt+5bXY8\nEZGHpFi4GzZsyJgxY4iOjubXX3+lffv2tG7d2qqhpk+fTpUqVejRowc3btyw6lgiIpI1GYZBQ++G\nrH59Ncs6LmPH+R0Un1qckNAQrkVfMzueiEiSVN3a/YsvvmDt2rUABAQE0LNnzzTNcjdr1oxLly49\n9PyYMWN47rnnyJ8/PwDDhw/n4sWLzJ079/7QhkFwcHDSx/7+/vj7+z91HhGRB4WEaKeSrOj4teNM\n3DqRH478QLeq3RjgN4AibkXMjiUiWUhoaCihoaFJH48cOTLNqy8eW7jj4+OpVKkSR48eTdMgT+vM\nmTO0bt36ofXiWsMtIiKPc/7WeSZvn8xXe7/ipfIv8UHdDyiTr4zZsUQkC7L6Gm57e3vKli3L2bNn\n0zTIk7h48WLSn3/88Ud8fHwybGwREbENXm5eTA6YzImgE3i5eVF3Xl1e/e5V3URHREyR4pKS+vXr\nExYWRq1atXBxcbl7kmGwfPlyqwTq2rUre/fuxTAMihcvzmeffYanp+f9oTXDLSIiT+D2ndt8vudz\nJm+fTKUClRhabygNnm2gvbxFJEXp0TtTLNwbN24E7t85xDAMGjZsmKaB00KFW0REnkZsfCxf7/+a\nCVsnkD9XfobWG8oLZV7AzkhxDwERyaasWrj/+ecfZs+ezcmTJ6lcuTJvvvkmOXPmTNNg6UWFW0RE\n0iIhMYH/Hvkv47aMIz4xniF1h9ChUgfs7ezNjiYimYxVC/err76Kg4MD9evXZ9WqVXh7ezN16tQ0\nDZZeVLhFxNq0S0n2YLFYWHNqDeO2jOPczXMMqjOI7r7dcbJ3MjuaiGQSVi3c/76bZHx8PDVr1ky6\n06TZVLhFxNp0p8nsZ9u5bYzbMo7df++mb+2+vFPzHdwc3cyOJSIms+ouJfb29sn+WURExBbVKVqH\nnzv9zNrX17L/8n5KTC3BsPXDuBx12exoIpLFPXKGO0eOHOTKlSvp43/++QdnZ+e7JxkGt27dypiE\nydAMt4hYm2a45c/rf/Kfrf9hyaEldKrUif5+/SmZt6TZsUQkg2XILiWZkQq3iFibCrfcc+n2Jabv\nnM5nez6jcfHGDKoziJpFaj7y+Hs7Der7R8Q2qHCLiFiJCrc8KDI2krlhc5m8fTKl8pbig7ofEFAy\n4KG9vFW4RWyLCreIiJVolxJ5lLiEOJYcWsLErRMBGFRnEB0rdSRnjrtb56pwi9gWFW4RERGTWCwW\n1p5ay8RtEzlx7QT9nutHz2o9cXNy/f+vmxxQRNKFCreIiIjJDAMovBvq/AdKrIc9b8HO9+B2wfuO\n0/+2RLImFW4RERGT3beE2+MU+E0Gn2/hcHvYNgCulQVUuEWyKhVuERGRTMYwgFxXCF75KTN/n0nd\nYnX5oM4H+BX1MzuaiDwFFW4REZFM5t8XTUbdieLLvV8yafskvNy8+KDOB7xQ5gXsjEfed05EMhkV\nbhERK9EuJfK0ktulJD4xnu8Pf8/ErROJiY9hYJ2BdPbpjKO9ozkhRSTVVLhFRKxE+3DL03rctoAW\ni4UNpzcwcdtEDl4+yPu13+ft6m/j7uSesSFFJNVUuEVErESFW6xt76W9/Gfbf1h9cjU9fHvwfu33\nKeJWxOxYIvKA9OidWkQmIiJigqoFq7LwpYXseWsPsfGx+Mzyofuy7hy6fMjsaCKSzjTDLSKSDM1w\nS0a7Fn2NWbtnMX3XdGoVqcUHdT6gXrF6D906XkQylpaUiIhYiQq3mOWfuH+Yv28+k7ZPIp9zPgb4\nDaBd+XbY29mbHU0kW1LhFhGxEu1SImZLSExg2bFlTNo+ib8j/+b92u/Tw7cHro6uZkcTyVZUuEVE\nRLKBHed3MHn7ZNafXk8P3x68V/s9vNy8zI4lki3ookkREZFs4Dmv51jafim7e+3mTsIdKs+qTOcf\nOvPHxT/MjiYiqaAZbhERkSzmRswN5uyZw7Rd0yiVtxQD/AbQsnRL3cFSxAq0pERERCQbi0uIY+mh\npUzaPonouGj6PdePrlW64pzT2exoIjZDhVtERESwWCxsPLuRSdsnsevCLnpX701grUAKuBQwO5pI\nlqc13CIiVqIdSiQrMQwDf29/fu70Mxu7bSQ8KpyyM8rSc3lPDl85bHY8kWxPM9wiIsnQPtyS1V2J\nusKs3bOY+ftMqhWqxgC/ATQu3lg30hF5QlpSIiJiJSrcYiti4mNYuH8hk3dMJqddTvr79adjpY44\n5HAwO5pIlqDCLSJiJSrcYmsSLYmsObmGyTsmc/jKYYJqBfF29bfxcPYwO5pIpqbCLSJiJSrcYsv2\nXXxXeRYAAB1CSURBVNrH5B2T+fnYz3T26Uzf5/pSMm9Js2OJZEoq3CIiVqLCLdnB35F/M2PXDOb8\nMYcGzzag/3P9qVO0jtZ5i/yLCreIiJWEhGinEsk+bt+5zVd7v2LKjinkz5WfAX4DaFe+HfZ29mZH\nEzGdCreIiIikm4TEBJYfW86k7ZM4f+s8QbWC6FGtB3mc8jx07L1JcP3vWGydCreIiIhYxe8Xfmfq\nzqmsOrGK13xe473a71EmX5mk11W4JbtQ4RYRERGr+jvyb2b+PvP/tXfvUVXVifvHn8NFNHUMzbLA\nMhUDBOEgYuaYmpfMC6mVlo2iKF4YreZazfwhzUWb1XK+OTYkkVimaaOm6HjLVEorRYHUJB1NMKHI\nvKGgoMD+/cEPgjggIpvN5f1ai5Vnn73PeaDP2j5sP+ez9VbSWwrxCNELD76gQfcPkpNTcePmr2M0\ndhRuAABQJ65ev6r3D7+v1/e9LsMwdGTJ89KhX8m43sLqaICpKNwAAKBO/LRwiSHdv1N6cKHkuVdK\nipD2R0qXPYqf5a9nNDK10TudaikLADQqrFACVMYmpQ2SVm6QlnwmuV2WIv2lJyZIHolWhwPqJa5w\nA4ADrMMNVK3shyYv5l1UXEqcFiUu0t2t7tbzvZ/XWJ+xcnV2tTYkUAuYUgIAJqFwA1VztEpJybKC\nr+97XScvnNTsXrMV0TNCbVu0tSYkUAso3ABgEgo3ULUbLQuY8n2KFu5bqPhj8Rrffbye7/28fNr7\n1F1AoJZQuAHAJBRuoGrVXYc7KydLiw8s1uIDixXYIVAvPPiChnYZKicbHyNDw0DhBgCTULiB2pVX\nkKdVX63Swn0LdfX6VT3f+3lNCpikls1aWh0NqBKFGwBMEhXFSiWAGQzD0KenPtXr+17X7lO7NdU+\nVb8O+bXubXOv1dEAhyjcAACgwTp54aTeSHxD7x58V4PuH6QXHnxBfTz7yPbTot+A5SjcAACgwbuU\nf0nvfPmO/rXvX2rboq3mhMzRuO7j5ObiZnU0gMINAAAaj8KiQm06vkmLEhfp8A+HFREUoZnBM+Xx\nCw+ro6EJa7B3mly9erW6d+8uZ2dnJScnl3tu/vz58vLykre3tz766CMr4gEAAAs4Ozkr9IFQbZ+4\nXbvCdun81fPye9NP49eM155v93CxDQ2WJVe4jx49KicnJ82YMUMLFixQUFCQJCk1NVUTJkzQ/v37\nlZmZqcGDB+t///ufnJzK/17AFW4AAJqG7LxsvXvwXb2R+IZaNmupOSFz9IzfM2rh2sLqaGgiGuwV\nbm9vb3Xr1q3C9vj4eD3zzDNydXVVp06d1LVrVyUmJlqQEEBTxwolQP3QpnkbPdf7OR2dfVTzB83X\n2q/X6r7X79NLH7+kUxdPWR0PqJZ6ter8d999J09Pz9LHnp6eyszMtDARgKbqlVesTgCgLCebk4Z1\nHaZNEzbp86mfK78gX0FvBWnsB2O1M20n//KNes3FrBceMmSIsrKyKmyfN2+eRo0aVe3XqWxpoKgy\nl58GDBigAQMG3GxEAADQAHVt21X/N+z/9NdH/qr3Dr6nOVvmyCabZofM1sQeE7mZDm5JQkKCEhIS\navU1LV2lZODAgeXmcL/66quSpJdeekmSNGzYML3yyivq3bt3ueOYww3AbNxpEmg4DMPQrvRdWpS4\nSJ+e+lRhAWH6da9fq0vbLlZHQyPQYOdwl1X2GwgNDdWqVat07do1paWl6fjx4woJCbEwHQAAqO9s\nNpseuf8RrRu/TknTk+Tq5KoHlzyoke+P1LYT21RkFFkdEU2cJVe4161bp+eee05nz55VmzZtZLfb\ntWXLFknFU07i4uLk4uKihQsX6tFHH60YmivcAEzGFW6gYbty/YpWHl6pRYmLdLXgqmb3mq2wwDD9\nwu0XVkdDA8ONbwDAJFFRrFQCNAaGYWjPt3u0KHGRPj75sSb4T9DskNnyvsPb6mhoICjcAAAA1ZR5\nKVOLDyxWbHKsetzVQ3NC5mi413A5OzlbHQ31GIUbAADgJuUX5Os/R/6jRYmLdO7qOUUGRyrcHi73\nFu4O9y9ZMI3q0TRRuAEAAG7Bvox9WpS4SJuOb9I433GK7BWpgA4B5fahcDdtFG4AAIBakJWTpdik\nWL2V/Jbua3OfIntF6gmfJ+Tm4kbhbuIo3AAAALWooKhAG49tVPSBaB364ZCm2qdq/hMzpOz7KNxN\nFIUbAEzCKiVA01V6k+t2x6TgxVLAMul0X2l/pPTNUMkovo0JVaRpoHADgElYhxtoukoLdwnXXMlv\nlRTyb8ntknRgppQyRcaVdpbkQ92icAOASSjcAEqUFPCiIkOJmYmKPhCtDcc26PEHHldkr0iFeHBX\n7MaMwg0AJqFwAyjh6EOTZ6+cVVxKnBYfWKx2t7VTZHCkxvuN122ut1kTEqahcAOASSjcAEpUtUpJ\nYVGhtn2zTdH7o7U3Y6/CAsI0M3imvNp51W1ImIbCDQAmoXADKFHdZQHTLqQpJilGcSlxst9tV2Rw\npEZ0GyEXJxfzQ8I0FG4AMAmrlACoqbyCPK1JXaPo/dHKuJShGT1naFrQNN3V6i6ro6EGKNwAAAD1\nWMr3KXrzwJtanbpaw7oOU2RwpH557y9lq7AUCuorCjcAAEADcDHvot798l1FH4iWm7ObIntF6ln/\nZ9XarbXV0XADFG4AAIAGxDAM7UzbqegD0dqVtksT/CdoVvAsdb+zu9XRUAkKNwAAQAOVcSlDsUmx\nik2OlVc7L83oOUNP+DwhNxc3q6OhDAo3AABAA3e98Lrij8UrJilGB7MOKiwgTNN7TmdpwXqiNnqn\nUy1lAYBGhRVKANQVV2dXPen7pLZP3K7Pp34um82mvnF9NXjZYK0+slrXCq9ZHRG3iCvcAOAA63AD\nsFJ+Qb4+/PpDxSTF6Ni5Y5oSOEURQRG63/1+q6M1OUwpAQCTULgB1Bdf//i13kp6S+8dek/B9wRr\nZvBMjew2khvq1BEKNwCYhMINoL65ev2q1qSu0eKkxUq/mK6p9qmaFjRN97a51+pojRqFGwBMQuEG\nUJ99deYrxRyI0YrDK9T33r6a0XOGHuv6mJydnK2O1uhQuAHAJBRuAA1B7rVcfXDkA8UkxSgrJ0vT\n7NM0NWiq7ml9j9XRGg1WKQEAk8yda3UCALixls1aKtwern3T9mn9+PXKvJwpv2g/jf1grLad2KYi\no8jqiBBXuAEAABqVy/mX9f7h9xWTFKOLeRc1ved0TQmcorta3VXlcTZb8X+pWOUxpQQAAAAOGYah\n/d/tV8yBGH149EMN6TxEM4NnamCngbKVtOsyKNyOUbgBAABwQ9l52Vp+aLkWJy1WfkG+ZvScobDA\nMN1x2x2l+1C4HaNwAwAAoNoMw9AXGV8oJilG8UfjNaLbCE0Pmq6H73tYTk62/7+PxSHrGQo3AAAA\nasR223kpYJkUFCs5FUjJ06SDYVLuneX2a+qVi8INACaJiir+AoDG6qdp3IbU8Yvi4u2zTjo5WEqK\nKP6v4UzhpnADgDlYhxtAU2OzSXLL1pt7Vio2OVbnrpxTuD1c4fZwef7C0+p4lqFwA4BJKNwAmpqf\nf2gy+ftkvZ38tlZ9tUp9OvZRRFCERniNkKuzq3UhLUDhBgCTULgBNDWVrVKSey1Xa1LXKDY5Vicv\nnNTkwMmaap+qLm271H1IC1C4AcAkFG4ATU11lgVM/TFVS5KX6L1D78n/Ln9FBEVojPcYubm41U1I\nC1C4AcAkFG4AqFx+Qb7WH12vt1Pe1pdZX+pX/r9SRM8I+bb3tTparauN3ulUS1kAoFGZO9fqBABQ\nf7m5uGm833htn7hd+6bt022ut2nwssHqG9dXS1OWKvdartUR6xWucAMAAOCWFRQVaPPxzYpNjtVn\n336m8d3Ha1rQNPW8p6fV0W4JU0oAAABQ72ReytTSL5dqScoSuTd3V0RQhCb4T1Cb5m2sjnbTKNwA\nAACot4qMIu04uUOxybHafnK7Hn/gcUUEReihjg/J9tOdd+o1CjcAAAAahDO5Z7Ts4DK9nfy2nGxO\nmhY0TRN7TFT7lu2tjlYlCjcAAAAaFMMwtOfbPYpNjtWGYxs0uPNghdvDNbTLULk4uVgdrwIKNwCY\nJCqq+AsAYJ7svGz958h/tCRliU5fOq2wgDBNCZwir3ZeVkcrReEGAJOwDjcA1K0jZ45o6ZdL9d6h\n9/RAuwcUbg/Xk75PqlWzVpbmonADgEko3ABgjeuF17Xp+CbFpcRp97e79YTPEwq3h6uPZx9LPmhJ\n4QYAk1C4AcB631/+Xu8dek9LUpbIyeak8MBwTQyYqA6tOtRZBgo3AJiEwg0A9YdhGPr89OeKS4nT\nh0c/1MP3PazwwHAN9xouV2dXU9+bwg0AJqFwA0D9lHMtR6uPrFbcl3E6fu64JvaYqHB7uHza+5jy\nfrXRO51qKctNWb16tbp37y5nZ2clJyeXbk9PT1eLFi1kt9tlt9sVGRlpRTwA0Ny5VicAADjSqlkr\nTbFP0e4pu/XplE/l4uSiQcsGqc+SPopNitWl/EtVHm+zFX/VJUuucB89elROTk6aMWOGFixYoKCg\nIEnFhXvUqFE6fPhwlcdzhRsAAAAlCooKtO3ENi1JWaKdaTs12nu0wu3h6ndvvwoftCx5WN0qWRu9\n05LVxb29va14WwAAADRCLk4uGtFthEZ0G6EzuWe0/NByzdo0S9cKr2lK4BSFBYTJ4xceluWzZEpJ\nVdLS0mS32zVgwADt2bPH6jgAAABoQO5sead+2+e3+mrWV1oxdoW+zf5W/m/6a/iK4VqTukZyzq/z\nTKZd4R4yZIiysrIqbJ83b55GjRrl8Jh77rlHp0+flru7u5KTkzV69GgdOXJErVu3rrBvVJlbwA0Y\nMEADBgyoregAAABo4Gw2m0I8QtTbM0Ry/ae2+HyoLfZo6bezpOQI2Wzzyu1fMmskISFBCQkJtZvF\nylVKBg4cWG4Od3WfZw43AAAAqqPCByTdT0p3J0upT5bbXFm1bLCrlJRV9hs4e/asCgsLJUknT57U\n8ePH1blzZ6uiAWjCyvwjGgCgATOM8l+60FlKfbLidhNZUrjXrVunjh07au/evRoxYoQee+wxSdIn\nn3yigIAA2e12PfXUU4qJidHtt99uRUQATdwrr1idAADQWHDjGwBwgBvfAEDjZMWygJZPKQEAAAAa\nM0vW4QYAAACsYMW/XnKFGwAAADARhRsAHJg71+oEAIDGgg9NAgAAAJXgQ5MAAABAPUfhBgAAAExE\n4QYAAABMROEGAAAATEThBgAHoqKsTgAAaCxYpQQAHODW7gAAiVVKAAAAgHqPwg0AAACYiMINAAAA\nmIjCDQAAAJiIwg0ADsyda3UCAEBjwSolAAAAQCVYpQQAAACo5yjcAAAAgIko3AAAAICJKNwAAACA\niSjcAOBAVJTVCQAAjQWrlACAAzabxGkGAMAqJQAAAEA9R+EGAAAATEThBgAAAExE4QYAAABMROEG\nAAfmzrU6AQCgsWCVEgAAAKASrFICAAAA1HMUbgAAAMBEFG4AAADARBRuAAAAwEQUbgBwICrK6gQA\ngMaCVUoAwAGbTeI0AwBglRIAAACgnqNwAwAAACaicAMAAAAmonADAAAAJqJwA4ADc+danQAA0Fiw\nSgkAAABQCVYpAQAAAOo5CjcAAABgIgo3AAAAYCIKNwAAAGAiCjcAOBAVZXUCAEBjwSolAOCAzSZx\nmgEANNhVSv7whz/Ix8dHAQEBGjt2rLKzs0ufmz9/vry8vOTt7a2PPvrIinjALUlISLA6AuAQYxP1\nGeMTjZklhXvo0KE6cuSIDh48qG7dumn+/PmSpNTUVH3wwQdKTU3V1q1bFRkZqaKiIisiAjXGXxqo\nrxibqM8Yn2jMLCncQ4YMkZNT8Vv37t1bGRkZkqT4+Hg988wzcnV1VadOndS1a1clJiZaEREAAACo\nFZZ/aDIuLk7Dhw+XJH333Xfy9PQsfc7T01OZmZlWRQMAAABumYtZLzxkyBBlZWVV2D5v3jyNGjVK\nkvT3v/9dzZo104QJEyp9HZvNVmFbly5dHG4H6otXXnnF6gioBY3xNMPYRH3G+ER91KVLl1t+DdMK\n9/bt26t8/p133tHmzZu1Y8eO0m0eHh46ffp06eOMjAx5eHhUOPbEiRO1FxQAAAAwkSVTSrZu3arX\nXntN8fHxat68een20NBQrVq1SteuXVNaWpqOHz+ukJAQKyICAAAAtcK0K9xVmTNnjq5du6YhQ4ZI\nkvr06aPo6Gj5+vpq3Lhx8vX1lYuLi6Kjo5k6AgAAgAatQd74BgAAAGgoLF+lxJHCwkLZ7fbSD1eW\ntWLFCgUEBKhHjx7q27evDh06VPpcp06d1KNHD9ntdqaiwBRVjc34+HgFBATIbrerZ8+e2rlzZ+lz\nW7dulbe3t7y8vPSPf/yjLiOjCanp+OTcCbNVNTZL7N+/Xy4uLlq7dm3pNs6dqAs1HZ83c+60ZErJ\njSxcuFC+vr66fPlyhec6d+6sTz/9VG3atNHWrVs1ffp07d27V1LxiiYJCQlq27ZtXUdGE1HV2Bw8\neLAef/xxSdLhw4c1ZswYnThxQoWFhZo9e7Y+/vhjeXh4qFevXgoNDZWPj09dx0cjV5PxKXHuhPmq\nGptSceF58cUXNWzYsHLbOHeiLtRkfEo3d+6sd1e4MzIytHnzZk2bNs3hfev79OmjNm3aSCp/05wS\nzJCBWW40Nlu2bFn655ycHN1xxx2SpMTERHXt2lWdOnWSq6urnn76acXHx9dZbjQNNR2fJTh3wiw3\nGpuStGjRIj355JNq37596TbOnagLNR2fJap77qx3hfs3v/mNXnvttdI7UVZlyZIlpTfNkYp/0xg8\neLCCg4MVGxtrZkw0QdUZm+vXr5ePj48ee+wx/etf/5IkZWZmqmPHjqX7cEMnmKGm41Pi3Alz3Whs\nZmZmKj4+XrNmzZL00/03OHeiLtR0fJb8ubrnznpVuP/73//qzjvvlN1uv+FvDLt27VJcXFy5OV2f\nffaZUlJStGXLFv373//W7t27zY6MJqK6Y3P06NH6+uuvtXHjRk2cOJGrhqgTNR2fJTh3wizVGZsv\nvPCCXn31VdlsNhmGUbofq5TBbLcyPqWbO3fWqzncn3/+uTZs2KDNmzcrLy9Ply5d0qRJk7Rs2bJy\n+x06dEgRERHaunWr3N3dS7fffffdkqT27dtrzJgxSkxMVL9+/er0e0DjVN2xWaJfv34qKCjQ+fPn\n5enpWe6GTqdPn5anp2ddRUcTUNPxee7cObVr145zJ0xTnbGZlJSkp59+WpJ09uxZbdmyRa6urhVu\nhse5E7XtVsZnaGjozZ07jXoqISHBGDlyZIXtp06dMrp06WJ88cUX5bbn5uYaly5dMgzDMHJycoyH\nHnrI2LZtW51kRdNS2dg8ceKEUVRUZBiGYSQlJRmdO3c2DMMwrl+/bnTu3NlIS0sz8vPzjYCAACM1\nNbVOM6PpuNnxybkTdaWysVnW5MmTjbVr1xqGwbkTdetmx+fNnjvr1RXunyv556SYmBhJ0owZM/SX\nv/xFFy5cKJ1L4+rqqsTERGVlZWns2LGSpIKCAj377LMaOnSoNcHR6Dkam2vXrtWyZcvk6uqqVq1a\nadWqVZIkFxcXvfHGG3r00UdVWFioqVOn8il7mOpmxifnTtQlR2OzMpw7UdduZnze7LmTG98AAAAA\nJqpXH5oEAAAAGhsKNwAAAGAiCjcAAABgIgo3AAAAYCIKNwAAAGAiCjcAAABgIgo3AAAAYCIKNwAA\nAGAiCjcAoMFKS0vTtGnT9NFHH1kdBQAqReEG0KS1atWqWvtlZ2frzTffrPSxJPXt27dWs9WUo2w1\nVd2fT3X9PFt6err8/f2rdWx+fr769++vsjdIzszM1JgxY5SVlVVuv4cfflhFRUW1FxwAbgGFG0CT\nZrPZqrXfhQsXFB0dXeljSfrss89qNVtNOcp2I4ZhlCuyJar786mummQrsWLFCo0cObJcpl/+8pfa\nuHGjJk2aVLrNzc1N/fr10/r16285LwDUBgo3APzMmDFjFBwcLD8/P8XGxkqSXnrpJX3zzTey2+36\n4x//qJdffrn08Ysvviip/NXgZcuWKSAgQIGBgQoLC5MkLV++XL1795bdbtfMmTMdXoF9+eWXyxXS\nqKgoLViwoMrjy77XpEmTHGb75z//KX9/f/n7+2vhwoWSiq8uP/DAAwoLC5O/v78yMjIq/Zk4eu/0\n9HT5+Pho+vTp8vPz06OPPqq8vLzSY/7617/K29tb/fr104QJE7RgwYIK2Ww2mwoLCyt9jbJWrlyp\nxx9/vNy2U6dOaf369dq0aVO57aGhoVq5cmWl3w8A1CkDAJqwVq1aVdh2/vx5wzAM48qVK4afn59x\n/vx5Iz093fDz8yvd5+ePy77WV199ZXTr1s04d+5c6eulpqYao0aNMgoKCgzDMIxZs2YZy5Ytq/De\nKSkpRv/+/Usf+/r6GhkZGZUe//P3unDhQoVsBw4cMPz9/Y0rV64YOTk5Rvfu3Y2UlBQjLS3NcHJy\nMvbt21flz+frr792+N5paWmGi4uLcfDgQcMwDGPcuHHG8uXLDcMwjMTERCMwMNDIz883Ll++bHh5\neRkLFiyokK2q1yiroKDA6NChQ4Xtr732mrFhwwZj4MCB5bbn5eUZ99xzT6XfFwDUJRerCz8A1DcL\nFy4snY6QkZGh48eP68477yy3j+Fg+kWJnTt3aty4cWrbtq0kyd3dXStWrFBSUpKCg4MlSVevXlWH\nDh0qHBsYGKgzZ87o+++/15kzZ+Tu7i4PDw+tW7eu3PF5eXnq0KGDsrOzy73X7bffrosXL5Z7zT17\n9mjs2LFq0aKFJGns2LHavXu3QkNDdd999ykkJKTKn8eOHTscZn/44Yd1//33q0ePHpKknj17Kj09\nXVLx9JrRo0erWbNmatasmUaNGlXpz6yy1yjr7Nmzat26dbltOTk5cnV11ciRI/X73/9eKSkpstvt\nkoqnlRQVFSkvL0/Nmzev8vsDALNRuAGgjISEBO3YsUN79+5V8+bNNXDgwEqnOFTGZrM5LJdhYWGa\nN2/eDY9/6qmntGbNGmVlZenpp5+u8vg33nijyvLvKI9hGKXzoFu2bHnDPIZhOHzv9PR0ubm5lT52\ndnbW1atXK33PylT2Go5ylLV06VLt2rVLX375pVxdXbVgwQItX7683P61PQcdAGqCOdwAUEZ2drbc\n3d3VvHlzHT16VHv37pUktW7dWpcvXy7d7+ePyxo4cKBWr16t8+fPSyr+oOCgQYO0Zs0a/fjjj5Kk\n8+fP69tvv3V4/Pjx47Vy5UqtWbNGTz31lCRVevwjjzxS4b1+nq3kA4RXr15Vbm6u1q9fr379+t2w\nqJe4mewl+vbtq40bNyo/P185OTnatGmTbDZblT+3qtxxxx3KyckpfVxQUKC0tDStX79eS5cu1bZt\n2xQfH6/Tp09LKl6pxNnZuVyZBwCrULgBNGlXrlxRx44dS7+OHTumgoIC+fr66uWXX1afPn0kSe3a\ntVPfvn3l7++vF198scJj6acVPbp3764///nP6t+/vwIDA/W73/1OPj4++tvf/qahQ4cqICBAQ4cO\nLbeUXVm+vr7KycmRp6en7rrrLkmq9HhfX98K7/XzbHa7XZMnT1ZISIgefPBBRUREKCAgoFzmyths\ntiqz//z4ksfBwcEKDQ1Vjx49NHz4cPn7+6tNmzZq27ZtuWw2m63S1yjL2dlZfn5+OnbsmKTiq/0H\nDhxQdna2JOnEiRNyc3PTmDFjlJubq5SUlNL/dwBgNZtR3UscAADchNzcXLVs2VJXrlxR//79FRsb\nq8DAwBq/3jvvvKMffvih9BecqvzpT39Sr169NGbMmBq/HwDUFgo3AMAUzz77rFJTU5WXl6fJkydX\nqyhX5dq1axo8eLA++eSTKq/M5+fna8iQITfcDwDqCoUbAAAAMBFzuAEAAAATUbgBAAAAE1G4AQAA\nABNRuAEAAAATUbgBAAAAE1G4AQAAABNRuAEAAAAT/T/WwNNuUwGhpQAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x45f7190>"
       ]
      }
     ],
     "prompt_number": 12
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 12
    }
   ],
   "metadata": {}
  }
 ]
}