{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# useful to autoreload the module without restarting the kernel\n", "%load_ext autoreload\n", "%autoreload 2" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from mppi import Utilities as U\n", "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Analysis of the Dos" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "si_scf_file = 'QeCalculator_test/ecut_40.save/data-file-schema.xml'\n", "si_nscf_file = 'QeCalculator_test/bands_12.save/data-file-schema.xml'" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "dos = U.Dos.from_Pw(results=si_scf_file,\n", " label='si_scf',broad_kind=U.lorentzian,eta = 0.1)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[(array([-12.89642137, -12.88642137, -12.87642137, ..., 0.97357863,\n", " 0.98357863, 0.99357863]),\n", " array([0.01111491, 0.01120202, 0.0112905 , ..., 0.01893493, 0.01877357,\n", " 0.01861477]))]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dos.dos" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnXmYHGd17t+v971n18xoJGuxbNmybNmW9xW4wTYYDAEMBLODYy77Jb6sDlvuJZeQmwRI7PgaSDAJDpjFNhizE2LwJsmSZUmWLWubRbP2LL2v3/2j6quu7q6qrq7ume7pPr/n0aOZnp7qb3p569R7zncO45yDIAiCaC9szV4AQRAE0XhI3AmCINoQEneCIIg2hMSdIAiiDSFxJwiCaENI3AmCINoQEneCIIg2hMSdIAiiDSFxJwiCaEMczXrgvr4+vmHDhmY9PEEQxKpk9+7ds5zz/mr3a5q4b9iwAbt27WrWwxMEQaxKGGMnzNyPbBmCIIg2hMSdIAiiDSFxJwiCaEOa5rkTBEGYIZvNYmxsDKlUqtlLWVE8Hg9GRkbgdDot/T6JO0EQLc3Y2BiCwSA2bNgAxlizl7MicM4xNzeHsbExbNy40dIxyJYhCKKlSaVS6O3t7RhhBwDGGHp7e+u6WiFxJwii5ekkYRfU+zeTuBNEGxNP5/DdJ08imy80eynECkPiThBtzP27x/DJH+7Hvz9xstlLIVYYEneCaGOmliTPdi6eafJK2ov3vOc9OHjwYEOOdfvtt2Pbtm24/fbbG3I8AVXLEEQbE0/nAACxVK7JK2kv7rnnnoYd6+6770YkEoHdbm/YMQESd4Joa6KyuC8k2yNy//xDB3BwYqmhxzx7OITPvmqb7s/j8ThuvvlmjI2NIZ/P44477sCdd96Jr3zlK9i5c2fF/fP5PN797ndj165dYIzhXe96Fz760Y/iyJEjuO222zAzMwO73Y7vf//7+OhHP4pYLIYLL7wQn/zkJ/HGN76xYX8XiTtBtDFROWJPZvJNXsnq5ZFHHsHw8DB++tOfAgAWFxdx55136t5/7969GB8fx7PPPgsAWFhYAAC85S1vwSc+8Qm89rWvRSqVQqFQwIMPPohAIIC9e/c2fN0k7gTRxgg7JtEm4m4UYS8X27dvx8c+9jF8/OMfx4033oirrrrK8P6bNm3C0aNH8cEPfhCvfOUr8fKXvxzRaBTj4+N47WtfC0DafbrcUEKVINqYRFYSdYrcrXPGGWdgz5492L59Oz7zmc/gC1/4guH9u7u7sW/fPlx77bW466678J73vGeFVloKiTtBtDGZnFTfnshSQtUqExMT8Pl8uOWWW3D77bdjz549hvefnZ1FoVDA6173OvzVX/0V9uzZg2AwiJGREfz4xz8GAKTTaSQSiWVdN4k7QbQxmZwUsbeLLdMM9u/fj4svvhg7duzA5z//eXzmM58xvP/4+DiuvfZa7NixA7fccgu+9KUvAQDuvfdefPWrX8W5556Lyy+/HJOTk8u6bvLcCaKNScuRO9ky1rnuuutw3XXXldz2u9/9Tvf+5513nmZ0v2XLFvzmN7+puD0Wi9W9Ri0ocieINkaxZUjcOw6K3AmijcnkKXJfTi655BKk0+mS2+69915s3769SSsqQuJOEG1MOiuJeyZfQC5fgMO+Oi/WOect2RnyiSeeWLZjc87r+v3V+UoTBGGKTL4Al0P6mCezqzN693g8mJubq1vsVhNiWEc99fAUuRNEm5LLF5AvcHT7nJiNZRT/fbUxMjKCsbExzMzMNHspK4oYs2cVEneCaFOE3x5wOyRxX6U93Z1Op+VRc51MVVuGMbaOMfZbxthBxtgBxtiHNe7DGGNfZYwdYYw9wxi7YHmWSxCEWUSkHvBIMVw21zm2BmEucs8B+BjnfA9jLAhgN2Psl5xzdTPjGwBskf9dAuBO+X+CIJqEIu5u6WOeya9Oz52wRtXInXN+inO+R/46CuAQgLVld7sJwLe5xOMAuhhjQw1fLUEQphEbmIIeZ8n3RGdQU7UMY2wDgPMBlNf/rAUwqvp+DJUnADDGbmWM7WKM7eq05AhBrDTp8sidxL2jMC3ujLEAgB8A+Ajn3FK3fM753ZzznZzznf39/VYOQRCEScptmWyePPdOwpS4M8ackIT93zjnP9S4yziAdarvR+TbCIJoEmm5aVjQQ5F7J2KmWoYB+AaAQ5zz/6tztwcBvE2umrkUwCLn/FQD10kQRI2UV8tQQrWzMFMtcwWAtwLYzxgTs6A+BWA9AHDO7wLwMIBXADgCIAHgnY1fKkEQtSDq2oPkuXckVcWdc/4oAMOmDlzaF/z+Ri2KIIj6EX1lipE7ee6dBPWWIYg2pbhDVSqFpMi9syBxJ4g2pWITE4l7R0HiThBtSmW1DCVUOwkSd4JoUxpd5/6Gu/6ID9/3dN3rIlYGEneCaFOK7QdEQtW6LZPJFfDU8Xk8sHeiIWsjlh8Sd4JoU4S4++XIvZ7eMtPRlPJ1oUBVN6sBEneCaFOELeN22OCy2+pKqMbTRb8+sUonOnUaJO4E0aZk8gW47DYwxuC0M2TrsGVSKkGPp3ONWB6xzJC4E0Sbks4W56e6HPVF7iTuqw8Sd4JoUzL5PNyNEnfV76otGqJ1IXEniDYlkyuL3Btky8Qocl8VkLgTRJuSVom70944cU9kSNxXAyTuBNGmZHKFoi1TZ7WMaEIGUOS+WiBxJ4g2RW3LuOv23NUJVfLcVwMk7gTRpqRzUikkINkyucLyVsv8+xMn8auDU5Yfg2gsZoZ1EASxCpFsGTsASdyzOes7S1MqWyapsYkpmcnjUz/aDwA4/tevtPw4ROOgyJ0g2pR0XpVQbUC1jMMmbYbSEvfZWLr4uNR9siUgcSeINiWdVdW517lDNZ0rwOO0w+O0I5kxFvfx+aTlxyEaB4k7QbQpmXxpKWS97Qc8Thu8TnuJ/y6Yi2WUr09EEpYfh2gc5LkTRJuSyZWLe32eu9thh0PHlpmLFyP3MYrcWwISd4JoU9JlCdV6SyE9ThucdpuOLVOM3MmWaQ3IliGINqVkE5OjTs89my967jq2jN9lx/oeHyYWzIs75xxfevgQnhlbsLw2QhsSd4JoUyptmXo8dymh6nNpe+6ReBq9ATfWdnkxXoO4H52N459/fxTv+84ey2sjtCFxJ4g2JZ0rVsvU77kXE6rannsGvQEXhru8NdkyJ+Xk63wiU+WeRK2QuBNEG5LLF1DgKNmhWledey4Pt8MOj0uvFDKDXr8La7u9mIqmTPv7Edmrt9uY5bUR2pC4E0QbIoTcVVbnzrm16F2yZUQpZKVwz8XS6PW7MdLlBefA1FJK4yiViIjdxkjcGw2JO9G2jC8k8cTRuWYvoymILo5qz51zIG9xuHUqm4fHYYfXaa9o+VsocERkW2ZttxcAMDpvrtY9IV8F0NDtxkPiTrQtn/7Rfrzx7sdLdk92CmnZFvE45VJIWeSt+u6pbAFupx1eV6XnvpTKIlfg6A24MSKLu1nfXYh7ug7LiNCGxJ1oW353eAYAMLloziJoJ0R/F3VCFYBl3z0tJ1Q9si2jjrRFjXtfwIWhsBeMmd/IJCpvMrmC5asKQhsSd6LtmY52oriXRu4uu+RpWy2HlDYxSbaM+vgAEIlL4t7rd8PlsGFN0GO6HFJt8VDDscZC4k60PbEOHC4hIuLyyN2KuOcLHNk8lz136Thqa0bYXj1+FwBgpNuLMZOee1LdSlijCoewDok70fYkO3Dmp4is1e0HAFjq6S5OFB6nDT6X1LFELe5iR+pQ2AMAWNttfiOT+rXRqp8nrEPiTrQl6kv8TowIRbWM21ns5w5Y89yL4i7VuQOlz+nYfBJ+lx1dPicAyZ6Zj2dNHTuhOo5WiSVhHRJ3oi2JpYoRYaIDI0JxcvM46vfcU4p/b1M8d3ULgrH5JEa6fWByrXqXz4lYOmdqI5M6Wtdqa0BYh8SdaEtKIsIOjNxT5ZF7HZ67OnIX4i5EeWophSeOzWFDn0+5f7ccwS8mq0fvyUwePlfpMYnGQOJOtCXqKDDRgeKuVwpZj7i7HXZ4XXJCVX5O3/mtpxBN5fCaHWuV+4d9UmJ1wUS/mGQ2j275/p1ony0nVcWdMfZNxtg0Y+xZnZ9fyxhbZIztlf/9ZeOXSRC1oY4COzEi1EuoZiwlVIu2jEcVucfTORyaXMJ7r9qIG7YPKffv8kqR+3yieuSeyOSVKhuyZRqLmWEd/wLg6wC+bXCf/+Kc39iQFRFEA1BHgZ0YEaorXACpnztgLXJPa9gyqWwez09FwTmwc0NPyf1FJG6m02Mqk0e3LO6deBJeTqpG7pzz3wOIrMBaCKJhUOSuUwppKaGqEnfZH09k8kq73s39/pL7Bz1SzKhOamvBOUcim0eP7NGnqVqmoTTKc7+MMbaPMfYzxtg2vTsxxm5ljO1ijO2amZlp0EMTRCUicnU5bJ3puWs0DgOseu6V1TLJTF4Zit3rd5fcX4h7NGVsy2TyUsuBLjnSp/4yjaUR4r4HwGmc8/MAfA3Aj/XuyDm/m3O+k3O+s7+/vwEPTRDaCEHv8bk6MnJP5fJw2pnSJ73YW6aOTUyOYuSezOYRiWdgtzGEZY9dEPRI30erRO6pjCTmXUrk3nmv03JSt7hzzpc45zH564cBOBljfXWvjCDqQAh6l89Z0gelU0hnC0qNO1Ac2pG18FwUI3c7XHYbbEwS/Ll4Gt0+F2xlgzZcDhvcDhuiaWNxT2Sln4sEbCe+TstJ3eLOGBtk8u4FxtjF8jE7s4k20TKIJGrY6zQ9FaidSOfySo07ADjrSKiqk7OMMWnUnmzL9MrJ0HKCHkdVW0Z5jeTIvRNfp+WkarUMY+y7AK4F0McYGwPwWQBOAOCc3wXg9QDexxjLAUgCeBO3Ou6FIBpEShW5z3RgP/dUtqAkU4HGJVQBKD3dI/GMUsZYTtDjrLBl9o0u4KnjEbz7yo1gjCnWmc/lgMtuo8i9wVQVd875m6v8/OuQSiUJomVIZvNw2Bj8bkdHVmGoh2MD9XruovJGOobHWRT3s4ZDmr8jRe6l4n7LN55ANJXDJRt7sX0krJyAvU473A4btfxtMLRDlWhLkpmCLBr2ugZDr1bSOWlyksBVR+SezkonCtE7Rpqjmsdc3Lwts5DIKGJ/YGIRAFSRux0uh41smQZD4k60JclsDh6XHBF2YBVGOlcoi9xlz91SQjWvWDKAZMtEUzksJrMVZZCCoLvUltk/vqh8PSFPxhLi7lEidxL3RkLiTrQlyUxedbnfeaKRypbaMnYbA2PW69w9quSsx2lX+rX3BLQj94DHgZiqWubQqSUAAGPAnJwDScrVMj6XHW6nnSL3BkPiTrQlyWypuHdajj+dK5RE24wxOO02a557rjRyD3kcOD4bB4Aqtoxa3KMYDHmwZSCgTG6qTKh23hXWckLiTrQlyWxBsmVkUcpaELXVTLoscgck391qKaS6Zr7X74aYZT0Q1LFlPFJPdzH0+vmpKM4cDKLX71Z2topSSK/LDrezM6+wlhMSd6ItSWZy8MmbboDOG75cnlAFJN/dUkI1V2rL9KqsmDUhj+bvhER/GdmamVhIYqTbi96AS4nck6qEqpsSqg2HxJ1oSxLyEAixkafTokKtyN1ZR+SuPlH0BorR+kBIL3Iv9pdJZfOYT2Qx3OVFj9+ltAJOZKUWCU67Da4OzY0sJ2Za/hLEqiORyUuX+44OFfeyaBuQxN1qP/eQqn9MnypyV2+UUhNwS/ePpXPIyZbYYMiDVDaPpVQWhQJXkt7iOGYmNxHmocidaEsSmZx8uS+JR6dd8qeyebjspcLrctTjuRelYqTbZ3BviWLknsPEolRZMxT2IOx1gnMgms7Jr5F0P7JlGg9F7kRbItkyDqXlbSd57pxzpHIFZSSeoD7PvXii2L42jOu3DeK1F6zV/R21LbMg2zCDYQ/G5BLKpWRWsc4AkC2zDJC4E21JstyW6aAWBNk8R77AFctDUI/nrrZ4XA4b7nrrhYa/o277e0retDQU9iLsjQGQhmeL1wiAvNmsc16jlYBsGaLtyOQKyBU4fE6VLdNBLQiS2dJGXwKrde5iz0AthFS2zORiCl0+J7wuu9L7fbEscu/UNhHLCYk70Xao66ddHRi5i3YLIioWuOw2S+0Hkpk8PK7axD2gEvdTi0kMyiWTJeKezcMre+6uDm0TsZyQuBNtR0LZ1u5QVct0jnAkVd0W1TgdtXvu+QJHOleAz1mbg+t12mG3MURTWZxaTGEoLIm7mLok2TLSXgQAHdsmYjkhcSfagu/vGsX1f/97LCayJd0GO7HOXVfcLXjuSlteV21SwRhTWhBMLqYw1OUFAENbJlfgyo5Won5I3ImqPLRvAp994NmW7s/yhYcO4rnJKB49Mlu6rb0DSyHF319upVjx3PVOFGYIehyYi6cxF89gSLZlvE47nHZWkVAV9lknvU7LDYk7UZW/+P4+/OtjJzA2n2z2UqpyeCqqRO7+Di2F1BNkK71lkhnt5KwZgm4nXpiSqmMGZVuGMWmg9kKiPHLvvNdpuSFxJ6oiLA1R0tZqFAocCVnQDk8uIZGRPPdO3aGa0rVlavfckzrJWTMEPA68MC2J+1DYq9we8joRiaeRzOaLm5icFLk3GhJ3wjSReGvOIl1MZhWv9sRcotRz78DL/WRG+lvLBdlpoVpG3dyrVkQ5JAAMdxUbjIW9TpyMSFeB3XKCtdjgTX99EwtJvPNbT2L3ifma19KJkLgThhRUCa5YujUvmSMJqYVsj9+F0UiiYnwb0FmRu361jHXP3YotI6Y02Vhpy4Kw14mTc1I/+G65H7xoTGZky3xv1yh+e3gGX3joQM1r6URI3AlDYpniwIV4Omdwz+YRiUvift5IGPFMXhkkEfI4ixFhB9VQ6wlyPZ67lYTq6QMBAEDAXcx9AJK4x+XjiuoZM/bZc6eiAIB9Y4tK22BCHxJ3wpAlVae+WKuL+7ouAMAz44tgTPJ2GWMdV0OdUpKg9feWqcdzF+I+3OUtuT2s6jDZ5ZMidzNXWFPRFBw2aRbswYmlmtfTaZC4txmjkQSeGVto2PGWkkVBb1VxF1cU24bDAIBnxhYQ9jphl4Wg08TdqP3ASkbul23uxesuGMHnXr2t5PYe1Wi+rvLI3WAn8Vwsg8tP7wNQnMlK6EPi3mZ86L6n8eqv/6FhvbGXUsXjtKotIy7xz1wTBAAsJLKKaACAy2HvOHEXQzDUSOLOa9qvUE/k7nHa8bc3n4dLN/WW3L6+p+i/iwHbZnoAReIZbO73o9vnxMlIoub1dBok7m3G0yelqP35qWhDjrcabJmEvK6+oAt98pQgcbkPiMi9czx3qYtjpRgL66OWebL1RO56bOjzK1+HPOWRu/brlMrmEUvn0Ot3YSDowXSUPPdqkLi3Eeqt25MNqklfkifYe5y2lo7cGQM8DjvW90j+riixA9Bxw5dTOl0cnXbJpqrFmqmnWkaPzX2SF3/GmoByW7WEqsir9AbcGAi5SdxNQP3c2wi1FTPToDe/iNyHw17EW7QUMpGWGlDZbAzrenzYc3IBa7uLSTy3w95hde55TRtF2DS1iruWxVMPYZ8TD33gypJB29XaD0TlICPkcaI/6MbRmXjD1tOukLi3EepNRo0qFROee3/QrURxrUY8k4fPLb2VRW218N+Bzpvyo9d/XQh0LX3TE+lcQy0ZwfaRcMn3wnPXe53ickmuz21Hf9CNmWganHMwxhq+tnaBxL2NmItllK8bJu7JHPwuO4IeByYWWrP9QCIjrREA/vyaTfC6bHjN+cURcO4O6xWezBa0PXd77Z57NJVTpiotJ9V6yyRV/YK6vC5k8gWksgVLid5OgTz3NkL4kjYGzKqEvh6iqSyCHic8TrvSs6TViKeLPUrWhDy4/bqtJYLkdthW3ZQfzjlGLVaEpDI6kbtD9tyrXMUcnFjCbw9PA5ByLkHP8seA1WwZke/xyYGGtLbGVIS1KyTubcScLO6b+wOYTzRK3KUPt9dpb1lbJpHJwe/Wj+BW43zOnx+YxFVf/i1+dXCq5t9NZuvz3F/x1f/CO7/1FJKZPGLprFLRspxUS6iqW0qEvGI+K4m7ESTubcS8LO6b+v3KxPl6iaazkri7WjhyzxQjdy3cDvuqK4V84lgEALDLQpOsejx3deR8bDaunNyXG4fdBruN6b5OShtnt0NpSLaYbM3qrVaBxL2NmItnEHQ7MBD0NDhyl2yZlo3c0yYi91WWUBU7gycWau+hn8zo1Lmb8NwXksX3zclIYsXEHZDWp2fLqNs4C8uNIndjSNzbiEg8g56AC90+Z0kb3HoQH27Jcy+UdIlsFRLVInenvmi0KnNy5ZPIo9RCKpvXHIunRO4Gz0UsVYyGRyMJJeeyEhjtRxBluD6nHWGv8NwpcjeCxL2NiMQz6PG70OVzgfPS3aVWER9ur9O4VK2ZxDM5w37jLvvqi9yFqM9ZEHc9W6Y4T1b/Cky9C/lEJI5YeuUid7fDIHLP5uBy2OCw25QcQCPe3+1MVXFnjH2TMTbNGHtW5+eMMfZVxtgRxtgzjLELGr9MwgyReAa9fhe6/dKbP9IAa2YplUPI44BXFoZWtGYSae0EosDtXH2euyhrna9R3DnnuuLukWvJUwbJZXXkfnBiCdk8R6/c0mG5MdqPkEjnlXLXYkKVIncjzETu/wLgeoOf3wBgi/zvVgB31r8swiyPH53D6+/8IxYTWUTiGXT7XOiW+6os1Cnu6VwemVxBSagCrSfu2XwBmXwBfsOEqhQRtvKA73IUWyaRqWndmXwBnFcOxwaKkbtRYjwqR+5ru7x4elTqUzQQXBlxN0p8q603t8MGp51RKWQVqoo75/z3ACIGd7kJwLe5xOMAuhhjQ41aIGHMlx95DrtOzOPXz02pPHdJ3Ofjku9uVdREZCQSqkBxM0mrkDAxBs5lt6HAgVwL5gu0yBc4UtkCfC6pbUK8huc8kdZv9OWpsgsUKEbuZw0FId42/Ssm7sYJVfEaM8YQ8jjJlqlCIzz3tQBGVd+PybcRK4C4xH7yWASZfAEDQU9R3BMZfPC7e/Cy//ufNQn8keko/vYXhxVLQNS5S4/XWuJenPFpnFAFWjNfoIW4OhoMSXNHaxEx4Zn73ZXPh8dE5C5+f+tgSLltpSJ3I1tG3WICkKwZSqgas6IJVcbYrYyxXYyxXTMzMyv50G2LaBb2++el53NNyI0u2XOfjWXw8P5JHJ2JYzRivqTuUz96Fl/7zRE89MwpAECXz6nYMq0m7qJEzihyV/qWtNja9RB/k4iYjbpx/uSZCbzzW08qdobowRLQEHe3iRO0EPezh4viPhj26N29oRhtNktmpOZwgqDHQaWQVWiEuI8DWKf6fkS+rQLO+d2c852c8539/f0NeGhCvMEn5Ba/a0IeBN0OOGysZCLTeA310kJcxO7IgaCnaMu0mEAKW8YwoeqovWFWMxFXI0Lcowbi/pH79uK3h2fw7Lg0mShuELmbmVMaS+fgtDNctKFHOo7LbnhV1EjcDjvSOq9RPJ0v2ctAtkx1GiHuDwJ4m1w1cymARc75qQYcl6gC57xigMZgyAPGGHoDLjx2dE65vZZNTaJS46A8ymwg6FZsmVbz3MXJxtBzNzHCrZUQJ6yBoBQxG0XuedluG5uX+tDEZM9dM3J32MCY8RVMLJVDwO1Af9CN77z7Ejz4wSut/REWcBk0eEtmS/cySJE72TJGVD0lM8a+C+BaAH2MsTEAnwXgBADO+V0AHgbwCgBHACQAvHO5FkuUksjkUeDAeSNh7BtbBCBF7gAw0u3DbtXWdbObYbL5AiaXit0fbUwakCD8zVaL3NUNpfSo1k621VDEPSRF7jEdEeOcK0nPKfk1E/fVEncxLDxVJXIPyHXtV27ps/YHWMQooRpPl+5lCHmcVC1Tharizjl/c5WfcwDvb9iKCNOIqP0lWwewb2wRm/v9SpS6vkcSd9H72qy4x1I5cC4NMY7EM/C5HLDbWMt67sUxcMalkIDx5p1WQthiA1VsGfVwlqklqXSyaMton+yqdfeMpnIIuFdmR2o5boNZt+W7kENeityrQTtUm0ShwHHLPU/g3sdPWD6GeHNv7PPjgfdfge+85xLlZ2II8aY+P4Ieh3lxl8Xhog3dAIq+b6vaMqZKIau0k201zNoypeIuR+5p/cgdkE50xgnVLII6v7vc6FXLcM5LSiEBqTw3kcnXNFWq0yBxbxIvzsTw6JFZ3PFjzY2/phAf5KDHgfPWdWEoXBwt9ydnr4HbYcPNO9eh2+cy7bkLYbnhnCG89dLT8JU3nAdAJe4t5lsnTHjuZhKJrUR5QlXPllFHrmKsolFCFZAi92oJ1cAKtRsoR2+QeTpXQIFLU5gEojOk3nND0CSmpqH2ta0iKmW0LqPPWRvGoS9cD5uN4Z5Hj5mefypOGF0+J774mnOU24VAtp4tI8avGdW5C8+9tdauhzjBhrwOuBw2xDLaAiZeq5DqyiyWkXqw6M089TiMbZlYKodNfU0Sd50Gb0pexVkauQPSwI5uv6vidwiK3JvG9FJxDJ5VwYylipG7FjabNHkn4LYbVlyo0Yv8bDZW9ZK+GSilkAZzPt2rzpYRSWIHgm6HbnQqbj+t16+Iezyd07VkAGkjk2FvmWZG7nKDt/INd4r1VraJCSi2RiYqIXFvEtPRorifWrQWxUer+KsCn8uhCEY1FHHXqG32ulqvp3syk4fbIQ160MO1ymwZdR7B73bonpjFhqX1vT7MJzIoFHhFPXg57iqRezSVa5rnLq6wyvcjKIM61AlV+QREG5n0IXFvEjMqcZ+2aNFEq0TugoDbUVEPr4dRQs7rtLdkQtXIbwdUnnuL5Qv0SGTycNoZnHab4WsnXv/1PT4UOLCQzCKWzhk3UTPomZ7NF5DOFXT9+uVG7wpLaxey2pYhtCFxbxLqSodFizvtxGV5tQ+jz2VXop9qGJXSteIc1WqDOgB1nXtrrV2PZCan2EwBt37JnxD90+TKqEg8jcVkVrEstDAqhYybvBJcLvSusLR2IYdoYEdVSNybRDxd/ABbFvd0Fl6nXTd5JvDXELmLDoTajadab45qMpszbD0A1NY4bGIhiWPaKzGpAAAgAElEQVSz8YaszSrqE5bfbVfsl3JiqRxsDFjbLVVJzcUymI9n0OPTTzAaVcuIk0gzq2WAyhyU1kmn6LlT5K4HiXuTiGdyGO6S6piti7u55JffLUXuZjpDxtM5OOTkaTmN9Nw55/jgd5/GA3s12xCZxowtI2aHmhH3V3/9D3jJV36nOaKQc449J+cbMr7QiEQ2r5T9SZ679nMek5OnvX6pZDISz2A+YVw9YpQUV0prmxS5e+UTWoW4ZyqvUAMuBxijyN0IEvcmEUvnMBj2SG9Qi+K+ZDL55Xc7kC9wU+ImtnkzVpmgbKTnfmw2jof2TeDD9+2t6ziJjPbUITVm69w555iNSbmQU4uVjdaeOj6PP/2nP+Lvf/W8xdWaI6k6YVXz3ANuB3oDkpjPxtJYSGTQ7dO3ZfwGFp2Sb2lS5B6QT2ixdHnkLhKqxdfZZmPo9buU14uohMS9SYiStaDbUZfnbipylyMiM+WQsXRe13MVQ7IbQaOsj0QmVzXnwBiDy6BviUD9Omi1SD4yHQMAPPbiXMXPGkk8nYPPKWwZg2oZ+cpN9O8/GUkgV+DoMYjcfW79yimjvjQrgbCiEmV/b0JnL0N/0GO5GKETIHFvElLJmgNhn7MuW8bM8GIhfmY2MsXT+mIp1Ug3JnKfUJV/1pPoTGSM56cKpBpq48dRR4FaIwpPRqTOizaNq5pGkswW/ya/2yE1iNOwgoQt43LYEPQ48IJ88uky8NwDbgeyea75XIjS2pUaiF2OEoSUXVmISN5XdoU2EHSXlBQTpZC4N4l4Rvpghr11iHvKeMOKQFzOmkmqxg0i4UZWy6gHP1v9+wHZwqhiywDGJYACtX+rNVx8VG6rW896zZAosWWk/xMaz3s0nUNALgkcDnux67jUBXTIYLiGOG5C40RfjNyb0zhMVGiVX6kkZKvQVraXYSDoVnrqEJWQuDcJESHXI+7RVNbUB1GItZmNTCIa1KKRCVV1I7PFhHWxNJNQBeSOg1UsJXXJ4YLGmsQJSUv4G0mypFpG31KLpYpNvtb3+pST97pun+6xxfG0TvSxtNzOolmRu/hby96ncZ1y1/U9PkxH0/jP52fIe9eAxL0JpHN5ZPO87sg9atqWMR+5Jwx2ODYyoaoW94U6I3eviUlBboet6iQm9W7HeY0ummLDzHw8Y3nouBnUHRADhmJcPBGLLqCMAUNd+pG7sD60kqoicjdzJbQcKEFIRUI1p1zBqNk+EgbnwNu/+SRuu3f3iqxxNUHi3gTEm9fnsiNscdCvmMJkypZx63+gyzHa4ShqpLX831qZT2TgtEuX2VpRshly+QIy+YKpyN1oyo9ARO6MAfMaaxJ9THIFbjj6rl7UVyNGyfBYqmihbRkIAAAcNma478HoRB+V30/l9sdKIaqeytcmnewq35Pnr+9Wrlx2nZhvyPuynSBxbwLqCfVBj9NSfwypbt1c8quWahlDz10M7GjATs9IPIONfX4A2slLM5hp9ytw6/QKVyNeh3XdPs01LSazymNpRfaNQJSsqhOqQKXgFQoc8UxesVBeetYAhsMe3H7dmYbHDxjaPOaCheXCbmPwOu0V9mFcp4Ir7HXil//jGuVvJmumFBL3JqCeUB90O5DKFmruWFjLbkIj37ZibQbVMo0c2BGJZ7ChVxJ3q7ZU8QrIjC1jr1otE03lJFsj7KnY8l8ocERTWcXPXq5uhOUzYYtirL2xR0SuA0EP/vCJl+LWqzcbHl8pN9TIv8QzzesIKfC7HRXVMolMrqSXu5rBsAdnD4UAFBPehASJexNQt9UV26hrjd6V5JeJSEsIRfmHppxiLkBvRJu8PbzO7oqcc0TiGayTfWKr49LEc2Dm6sVMnbvYFNSlUZ4az+RQ4MWt/svVjVCIrlfVfgCoPDFrbTjS2nhWTtHD16i+aXLkDsjtFjT+VqNmaCPyazI2X7k3oZMhcW8CcWVCvV0RploFzmxHSABKS9xqkbuyE9BgExNQf+SezOaRzhXQF3DD77Kb7ntTzlINVy9mbJmlVBYhj1Nz+LJ4LCEky9WNMFFW062XULW64UhEwFqRu9l9E8uJ31XZbqFaRVRvoNh+gShC4t4ElMkyLofSurRWcS+O2KteCskYM9x2Xr6uarZMvRuZxIew1+9CwKM/jKIayrASEwJnps49mpLELeR1VrSEEN8XxX15bJnymbB6llotJzY1RtU3S8lsS0buS1U6XYa9TjBmLjGfzuXxd798Hs+MLdS91laHxL0JqHumi6EDtUaC0RojNzOdIZUGTTqXwCLJV2+t+3xc+lu7/S4EPU7LkXstvVDcDrsJW0aK3MNeJ+Jlw5eL4u4r+b7RJLPClpGea6nPT6W4C1soVKO4ux022Jj2JqbFZA5hAxFdCcp76WTzBcQzecN12W0MIY/TVGL+u0+cxD/8+gXc8cCBhqy3lSFxbwLxsmoZwILnbkHcq21iMurlDjQuoToXl6oaevxOBNwOyxZHLc+By0T7ASVy17DKhAc/3OWt+FkjUaYOyX+TdNXlqPDIi7ZcbWLMGJN7xGuUeqayTRf3Lp8LC8nK3cvV1tXlc2qWr6rhnOPbj58AADw/GV3WvQqtAIl7Eyj2TC967rVe5itveoMOgGr8LnvV3jIxJRdQxXOvN3KXI6xunwtBj/le8+UovVBM7NJ1V5kdCpTaMkBpdC5en25ffSekamjNhPW5Kq2KWnIu5XT5XBUJ41Q2j0yuYGh/rARhrxML8cpBNtXF3VV1M9yuE/M4OhPHWUMhJLP5qieD1Q6JexMQPdNddhtCHmtDBxaSGdiY1NfaDD6XfndB9boAA8/d1SjPXfpbe/yyuNfpuRvNDBWYaZ0QTWURlG0ZoLREU7w+IY8TQY/+dKR6SWYqa/cDbgdiGT1bpnYx7taIcs2K6HLT7XMhms4plpjZIKbLW92W+fmzk3A5bHj3lRsBaLd1bidI3JuAqCVnjCl+ca1iIcapmd1NqFU/XI7R/FSgcQnV+XhG8UmNxshVQ0yiclSZRAVISddMrqBrzXDOKyP3lDpyL5ZdhjyVCddGUUyoFl8Drba/S6ks7DZmagNXOWGNKLdlxN1femIVfYeqravb56yaUN03toDta8PY1C/tr2j3pmMk7k1A3TPdbmOWBG4xmUNXDR9ErSqEinVVudRvVCnkXFwaKGGzMQTc9SVUzVaL6G0GEqSyBeQKHEG5FBIo3ai0mMzC75JOJMsZuRfr3IuirfXaiZp0M7Xt5WhFua0i7uLxxfpqsWXmq0Tuh05Fcc5wCIMhqffO5GJ772glcW8CUuSunuSuneAyYjFZW/LLTEI1ZrIUMlnnwI75eEYZMCE8dyt9QaImJ1EBUFrj6p3gourIXB6+XD7EXPRJD3kr6+Abha4to5FQtVqTrhXlmo2QlxvxvhC2US0J1Wgqh5xOc7hEJodYOoc1YY8yzKTayWC1Q+LeBEQvd0HQU3uCrtqU+3LMJFTj6RzcDptu4ykxrq7ehGokkVHmfAqBKveUzWAlcteLuJdUVy1hLVtG9Xwva+SezcNpL23+pWUDifyAFcI+F5ZS2ZJZsK0SuXf5ROReo7grr5n26zIblYS8P+CGx2mH22Fr++HaJO5NYClVHLIASB/eWsViqcbI3edyIJnNGw53rtZC2GZjDZnGNB/PoFcWd2VTjQWxrKXRlXIS0Yncl1QJSq/TDoeNlXz4FxJZhOWIXmsHa6NIasyE1apuWaojcu/yOsF5eTWQ/Pe3TOSeUf73u+yGnS6B4vQpvaTqjNxUrC8o7WYNe6t79KsdEvcmoB6yAFiLBGu1ZYQIGlkzZsSyET3d50sid+lvsOK7m215DKh3Zmp/oNWlhYwxhMr67C8ms+jyCltGer2Wo05aayZst0/KS6g3YUVTuZo3MCnHk5OWalsiEpeqr6wes1EIy0R0eJyJpjEQ0u9PL9CqcFIjjtcfKIr7ck/UajYk7k2gXJSCNUaCnHPZA64hclcaUOkLsxmbo95Re4UCx3wiix450ipWC9X+QYuaHBBe+jjVPHfpOS3vs68+mQY9TuQL3FR//FrRmgkrXueSk00iYznK7vJVes7TS2n0BtymKo+WE7+8a/vUQkpZV78cbRshSiX1at2FuPfJ4q7VHK7dIHFvArGyS2oRCZplKZVDvsCVSNIMRj1F1OuqFgl76hy1J7ze7jJbxoqHHUubT6gGq/z95ZuCQh5HqS2jOpmGLPYDMkNSo0lWueXAOcdcPKMIVa1oVYtMR1MYMCGiK8Habh/GF6Qa9KloCmtMRO7Cc9cb2TgjD9LuDUjPZdjrrGsC2GqAxH2FyZcNWQCgDOwwe5k/E5WimoGQ+Q+jKXFP56rOZPU67VUnGhkhmob1+IvJyWrr0qJQ4FiqIakonm89b19dLQOgxJYp371Z3FXceHHQarvbVRaVxjNSV02Rt6gVMUBbvYlnOppuHXHv8mBiIQnOOaaWUlhjJnI3Yct0+5yKd6/VHK7dIHFfYbQ2CgU9DmTzvGrXQsHUkhSFDASrRzSCam9+sbZqSTpPnbaMsAJ6/NIH1nLL43QOnMO0NeV12mFj+ieRxWRW2XMAlCZNyys2rPbgN4PWCUtJMooB3THxHFoT97DXCY/ThsnF4iaemWi6pvfTcrK2y4vx+SRmYmmksgWlh74Rxfp4HXGPll7pdHkrk9TtBon7ClNs1VsauQPmWxCInXVraojczYp7ta389SZU54Qw+eqrlqm1LrvYMEtf3KXWsdKmICmyk+4rBEOcSJTIfRmmMWnVr5cL12y81D+uFcYYhsJeTMrvo3yBYzaWrulKcDkZ7vIims5h70mpLe+m/kDV33HYbQi6HSVNx9TMxtIlz1fYKyWps1WGpq9mTIk7Y+x6xthhxtgRxtgnNH7+DsbYDGNsr/zvPY1fansgoj21/VFs+2tOLE4tClumwZF7qrotI0Xu1j8QStMw2ZbxuxxgDDUPnBYfYuFHm6G8AqbkeInS6qOQV/LcOefKmpVqGXEyXhZbJlvRL6a7bNNNvZE7IPnuInKfi6dR4GgZW+aMNUEAwE/3nwIAbJTHMVYjbJAknYmllTJIAEpZaztbM1XFnTFmB/CPAG4AcDaANzPGzta4639wznfI/+5p8DrbBqVNrTqhWmPb3+enohgOe2oarKDV6VBNMpNHJl9Qdmfq4XXZ66pzVzcNA6Ta+YCr9h265dG0GXr8+lvUy0tL+/xuZPIFLCVzxTI6WRxqPRmbpVDgmtaY32WHx2lTkoKiZbJIDlphqMujjKUbjUizR0Wv+mZz7kgYAPDA3gl0+5zKgJRqhL1O3YTqbDSNPtXzJYKCdrZmzETuFwM4wjk/yjnPALgPwE3Lu6z2Jappy9QmFgcnlnD2cKimx/U47XAZ7MpTt+E1wuu01WXLzCcy8DhtJY2xrExjEsnFWvrrdPtcim9dTrm4D8pJx8mllCKqIrJdLs9dzGktF3fGGIbDXuWKbU6ZZGU90t4yEMTkUgqLiSyOzsQBABv7zEXIy01vwI1z1krv78tP7zPdHK/Lp10Bk8zkEc/kK2wZQL90sh0wI+5rAYyqvh+TbyvndYyxZxhj9zPG1jVkdW2IUnLnrvTczYhFKpvHizMxZeJ7LRht3BCRcHeVSNjncigTm6wwF8sofrugfPqOGWrtZw9IkfucgbirrwLKxd1pZ4ogSC0aWMM9d6MBHENdHkzI1S2TiykEPY6Kevha2DooWR+Hp6I4NhuHw8ZMR8grwSeuPwuXb+7FB15yuunf0UuSll95AcX3TadH7mZ4CMAGzvm5AH4J4F+17sQYu5UxtosxtmtmZqZBD726ELXKaq9YWCFmKkYOT0ZR4Kg5cgeqiXtGvo9x5B6SE1FGbQyMUO9OFVgZ2LGorLcxkXu55y5qwacWJXHvC7iVCJIxJreMaKwwGA3gGAp7MSHXfo/NJ7GuTgvlTCHuk0vYP76ILWuCTd/ApObKLX349/deirNqCGJCOi0FpqOlu1OB6nXx7YCZV3McgDoSH5FvU+Ccz3HOxY6IewBcqHUgzvndnPOdnPOd/f39Vta76hF13uoosZbI/eCpJQDA2UPhmh/bSNxFFz6R6DQ6BufWLYlIPFORCAx4nDX71wsJqZe722E+eu3xS7NRy3MGomZebfGIypFTiymMLySVSF4gNXtrdOReuktWzXCXF9PRNLL5AsbmE6bKA40YCnvQF3DhsaNz2Du6gPPXd9V1vFZA2nWaqdgvUr47FTBXYLDaMSPuTwHYwhjbyBhzAXgTgAfVd2CMDam+fTWAQ41b4uqFc15RajUfzyDkcZQ0QvK7pBpsM5f5ByeWEHQ7LF1CG4u7Oc+9q0o9cTXmE5mKxwi6HYjVmlCtsf0CUKw6KV+71CemtGmW22FHr9+FyaUUTswlsKGsYiPkXdnIfTjsAefAxEISo5H6I3fGGF66dQAP759ENJXDtWes/mCry+tENs8r9mEo4h4svu+q1cW3A1XFnXOeA/ABAD+HJNrf45wfYIx9gTH2avluH2KMHWCM7QPwIQDvWK4Frybe9509uO7vfl8i8JFEtiJyZYwpu1SrcfDUEs4aCplOMqkxtiUqryi0KN8tWStakbslW6bGxmlAsbY+UvYcRHRObOt7fTh0agkTi0ms7ykV0+WYxlTsTFkp7lvWSLXePz8wiWQ2j61Dwbof79arN8HnsuOMNQG8ZOtA3cdrNuK9Wf76ina/6gS0w25DwKAuvh0wVUvHOX8YwMNlt/2l6utPAvhkY5e2ukll83jkwCQA4NCpJZw7Il32zscrPWfAXGfIQoHjuVNLeP2FI5bW1B90YyaWBue8YoLPfCILn6u6zaHVxMosmVwB0VSu0paxMokqUXvkXt5xUDC9pN3OYdtwCN95/CSAokctCHocDR/TNh/Xr90/aygExoDvPinVNlhJqJdz+kAQj33yZfA49Xv4ryZEwnQmmi4p65yJpdDlc8LlKP0b270z5Op/RVuUQ7I3DgDPjhe/jqh6masx0xlycimFeCaPLWusRW39QTeyea75htayS7QQCddqw4i1WFA2MJV77g4kMsa95suZT2RqjtyHuyQrS73tHigm3Mq331+wvlv5+pKNPSU/6/I5Gz7JJxLPgDFta8zncuC8kS4cm42j1+9Sql3qJex11pS3aGXE6ydeT0F56wGBUV18O0DivkyMzhebMqmFXk9EzSTojs/WV4+sjmzKmYmmTW2KqScRJewPrVJIoLbmYbMxc61g1awJecAYlI6DAhGBl+/QfMX2IZyxJoCbd46gt0wc1oQ8mI1lSnqs18ucPH7QrmO5vfOKDQCAt1x6WktVtrQK4spruuyKanIppdmqo93b/ja3M38bMynXJG/q9+OYLMqcc03PGZA83HLRKedoveIeKIp7efQ/tZSqSBpqUU8iSm/bvHqHrploPJsvYD6RRX+gtkZXLocNfQF3STdEQHo+XHZbhc3jcdrxi49eo3ksUSo5HU01bGen3ntDcNOOtbjmjP6mj8JrVXr9bthYZeQ+tZTC5Zv7Ku4f9jrxwnRspZa34tDpf5mYXEzD57LjvJEuHJ2R3kBLyRzSuYJmxNnlq5xIX87x2TjcDpsiLLWiRO6xysj91GKqotxPC5fDBr/LbkncZzQ2kwDFVgxmNwWJ5mO1Ru6AZM2c0rBl+oPuijyEEeK5aqTvPldF3AHJj69lnZ2E3cbQH3SXvCb5Asd0NI3BcOdF7iTuy8TUkiSWm/r8mFhMIZnJK5G58H7V9AfdmJWTnXocn4vjtF6fpUoZ8RiANN1GTSKTQzSVMyXuANATcCn9TWpBPG65KIscRHmVgx7CVrIk7mGPshlIWVc0VXNHxEGlJ7okJIUCxy8OTJr+G7TQy8cQ5hkIekoi99lYGvkCx2C48jMX9rqwmDA/R2G1QeK+TJxaTGIw5MHGfsnqODYbV+yAIQ0R7Q9IyU6jiHg0ksT6Huv9P0IeB/wue4X9IxKMZq8IBoKeihOEGWZiabgdtopSP9Gtr7yKRf840nqtiPu6Hh9G55MlydvRSBJrNU64RhSnGUlr+eHT47j13t2444Fna16TYC5mLu9B6DMU9mBcle8yem+HvU5k8oW65hO0MiTuy8TUUhqDIQ829Un1ycdm41Ujd0DbMgEkv/5kJIF1PdZ3JjLGsL7XjxNz8ZLblQ+Aych9IOjGdLR2O2JGx/4QlQymxT0qdhzWLoRbBgLI5ArKc5DJSTs+zeQb1IS9TvhUJ8pfyGWvfzgyi4KF1gy5fAELyawyxISwxsY+P05EEsrJ+5SBuNdT1rsaIHE3Cecc+8cWTW00KhTk8WBhDzb0Scm2ozMxvDgdg99l1+ybbVTJAkh+bDKbr9hMUyun9fhwQm7xKjgpWr52mTu2JO4WIveodoVLyOOAy27TPbFpHQewNqxC1Ks/PxUFIFXOFDiwocYkNWMMm/sDODIdQzZfwGMvzsHlsGEhkcWLM7Un6aajaXBu/uqJ0GZjnx+ZXEGx3sRJXOtzI3Zb12OltTIk7ib59aFpvOrrj+Lmf368asnebDyNXIFjKOyBz+XAcNiDF6ZjODwVxZY1Qc2EmOKH60TEQoDrFvc+H8YipbbEC9MxeJw20/1KBkIeRFO5mvu6z0TTJc2bBIwx9AVcyk7CaozNJ9EXcMPjrL0++/QB6Urq8GRM/l8qU93cX7vdtWUggBemYtg7uoBoOofbrt4EADg0Ga35WBPKVR2Jez2ISjJRoXZ0RtoXoNU9tFpAtdohcTfJj/ZKvdIOnVrCR+572jAJMyVPlRdT23es78KTxyI4MLGEs3S2jQ9UeaOJgQrr6o7c/cjkCyVJxSPTMWzqC+jWV5ejl5itxnQ0peuT98kJZTOMzSctt6f1uRw4fSCAPSfnAQAHJpZgt7Gaug8Ktq0NY3Iphf94ahQ2Brzt8g1w2BieU+1rMIuwd2r1/olSRI7riFzi+OJMDJt0Ttx6m57aBRJ3kxyaWMIN5wzijhvPxq8OTeNXh6Z17ytmU4pL7Es39WJyKYVoKoertmg3aAq4pWTnxIJ25K6Ie5011eLkcmCiKEBHpmNKRGsGcdKaqsF3z+Sk2nS9Icx9gVrEPVFX7/ErT+/DE8fmkM7l8eSxCLYOBi1dBVxxei8A4P7dY9ixrgt9ATdOHwjgOUuRu/RcDpG418VA0IPhsAe7T84jly/gwMQStg1rd1AVFVIUuXcwuXwBo/MJbOjz4+2XnYZunxMPy/MdtSivinnN+WvRH3RjOOzBNTrd9xhjWNfjU+yXckYjkhVRz4AGQOpR4rQz7BuThg/PxzMYX0hW9E4xQgjrqM5atdAamKBmSKNEUYt8gWN8IVnXFcxVW/qQyhbwwN4J7D4xj6stdkQ8c00QF2+Q2hK8+0rJkjlzMGgpcp9YSCLsddY0OpHQ5uKNPXjyWATPTUaRzOaxY512O2OP046gx1Gxo7VdoHeSCU4tppDNc2zo9cFht+HKLf14/Oic7v0nFlJw2pmS8At5nPjNx66B024zjBBP6/Upl5PlnIjE66qUEXicdmwdDGHPCcmW2C3/v/O0bqNfK2Gk2wvGgBNz5sVdzOvU8/XX9fgwn8gils4ZCtzUkvRa1GNfXH1GP4bCHvzP+58BY8DrLtAaLFYdxhjuecdOHJuJ4zxZQLYOhvDA3omSrpW/PTyNB/dO4IuvOUf3bzsZSZAl0yAu29yLH++dwJd+dgg2BlxxeuXuVMGA3EyvHaHI3QTH5Yz7aXK53LbhEE4tpnR3lE4sJDEU9pZsNgp6nFUv/U/r9WN0PqlZSndkOo4tNVgnRly5pQ+7TsxjLpbGbw5PSztpdaIbLdwOO4bD3poi96KtpCPust1U7ZiiEmVzv/Xnwmm34e/fuAPnrevCHa88G6cPWG/CFfI4S5470Yr3sMqauePHz+JHT4/jB7vHdI9TqzVG6PPKc4cR9jrxhyNzeOnWNYb7Iazu2VgNkLib4LgcoYpaaNFu9aDO5ffEQtJS1cP6Hh8yuUKFlz0XS2M2lsYZFrtBlvPq84aRL3B849Fj+Mm+CbzsrDU1e87rerwVJZVGjM5L99Wq8RfHA6BrSwmEaJ6xpj4hvGRTLx54/xV415Ub6zpOOWcNSu+N5+QqnCPTMeWq5ac6Vl48ncP4QrJhJ+9OJ+B24F/eeRHeecUG/O8/PcfwvmbtwNUIibsJTszG4XHalIoWMb/04ISBuGtsd66GiNzKE3LPT0nRai2+uBFbB4O4+ox+/NPvXkQik8dt12yq+RgbeqWGaGa3bo9GklgT0i9fNBu5Pz8VRV/AXdGlsVVYE3Kjy+fEoVPSa/jIs5Kg37RjGHtPLmiWj4qrEautnIlKzl/fjc++aptuAl+wvteHU0sppHPtt0uVxN0Ex+cSOK3Hr9gsfQE3+gKukktvQTZfwFQ0rRuhGrFNPmnsH1ssuV1cIZzZoA8/Ywxfe/P5+Pj1W3Hvuy/RrSYwYutgEJF4xnQZ2eh8wrDSp8vnRI/fhRemjDcAPTcZrTtqX04YY9g6GFTq5x85MIkL1nfhNTvWIpMvKDkONSJIaOW/q105rdcHzqXgo90gcTfBCblhl5ozB4M4PFUp7iflrc+17ngEJF9+U58f+8dLxf2pYxGMdHsx0MDdi2GvE++7djMu29xr6ffPlk8IetaUGs45XpiKKqPitGCM4eyhkOHxkpk8Dk4s6VY/tArnjXTh2fElHJhYxLPjS7jhnCFctLEHdhvDYy9WJuJ3n5hHt89puZUzYR2xKfBkJF7lnqsPEvcqFAocJyKJCrE+c00Iz09FK6YHvShXu1hNju1Y34VdxyPKcTnnePJ4BJdstCbCy4VIHOpZU2qmo2nMJ7JVrzzOHg7h8FS0Yqi4YO/oAnIFjos29Gj+vFV4+bZBZPIFvO7OP8JuY7jxvCEE3A5sXxvGYxpVVrtPzOPC03qolW8TEI34aqn8Wi2QuFdhcjx81gYAAA2OSURBVCmFTK5QEblvHQwilS1UJACPyP6p3q64arx06wDmE1k8Le+g3HNyAZF4BpdbjLCXi5DHic39fjx1PFL1vmIS1dYqu0C3DYeQyRU07S4AeOzFWTBWOv6uFblgfReuPqMfqWwBf3bxegzJ+ZfLNvdi3+gC4qr2FaORBI7OxivG+BErQ1/AhS6fU+k11E6QuFdBjLYr7xookpvlQnRkOoY1IbcyXahWrj6jHy6HDffLZXP37x6Fx2nDdecMWjrecnL55j48eSxSddTcs7LNVG3u52WbpBPYo0dmNX/+i4NTuOi0Hs0+Ia0EYwz3vG0nfvz+K/C5V29Tbr9sUy9yBY5dKt/953I3yeu2td7r2wkwxnDOcLhkznG7QOJehRd1IvEtawJgrFjyJtg7uoBzLCQoBSGPE2+6aB3u3z2Gbzx6DPfvHsNrzx9pyZ2LV5zei0Qmr+x21ePxo9IW/64qA7gHQh5sHQzid4crWzscmY7huckoXr5tTV1rXilcDht2rOsq6dezc0M3nHaGP74onbwKBY77nhrF9rVhrO9tzKg+ona2DYdweFLfDlytkLhX4YXpGAJuR0UrVp/LgS0DgZLqh7lYGkdn4thZpyf8oZdtwVCXB1/8yUEMhb34i5efUdfxlovLNvXBYWP45cEp3fukc3nsOhHBpZvM2UrXbRvEE8ciFSWR337sOFx2G15zvrXdpK2Az+XApZt68cM940hm8nhg3ziOTMfw7gbX2hO1sX0kjEy+oFxhtgsk7lU4Mh3D5oGAZrLr8s19eOp4RKmRffKY5D/v3FCfJ9wXcONnH74a//aeS/CzD1/VsjXdYZ8TL906gB/uGdeNen7//CxS2YJuT51ybr5oHRiAe/7rqHLbsdk47ntyFK85f9hSD/dW4sMv24KZaBpv/n+P49M/ehYXrO/CjecONXtZHc3lm/vAGPBfL2jbgasVEvcqvDAdw+k6W92vOF1qQCXK2x56ZgJ9ARfOb0CpXsDtwBWn98HfgnaMmpt3rsNsLK3bSO1HT4+hx+/ClVv0+3uoWdvlxZ9dsh73Pn4CD+2bwLHZOG67dzfcThv+4uVnNnLpTWHnhh588TXnYHophZ0benDXLRfCYaePYTPp8buwfW0Yv35Ov9PraqS1laPJnFpMYiaaVjYXlXP1GX3oC7jw7cdO4PSBAH51cBp/dsn6jvqwvmTrALYOBvGVXxzGddsGS3agHpmO4ZFnJ/HeqzfBWcNz8vHrt2L/2CI++N2nAQB+lx3//NadDa3zbyZvvfQ0vPXS05q9DELFq84dxv96+BAOT0Zx5mAQY/MJfPEnB3HB+m78+TWbm708S3SOClng6ZNSovACnY6Jbocd77xiI37z3DRu/NqjsNsY3nt17Vv5VzN2G8MdN56N0UgSH//BM8jJ9kwyk8ft9++D3+3Ae6+q7TkJepz43m2X4StvOA933Hg2fvE/rjEd+ROEFV534QjcDhv+7pfPI57O4T3/ugs/PzCFL/3suYod46sFitwN2HV8Hi6HTWkUpsWfX70Ji8ksnhlbwIdetqUj27ZecXofbr/uTPzNzw/j6Ewc15zRj18dmsLhqSj+8c8usOSTux12vP7CkWVYLUFU0uN34cP/bQu+/MhhPPXlCBaSWdz5lgvw0e/txX1PncT2ke3NXmLNkLjrwDnHb56bwiUbe+By6F/gOOw2fOoVZ63gylqT97/kdIx0e/EPv34BX//tEWwZCOCet+3Ey85aHaWLBHHb1ZtRKHA8cSyCWy49DddtG8TPD0zioX0TuOPGsy1N62omJO46HDoVxfG5RMfZLPVw0461uGnHWuQL3PQ8VoJoFWw2hg+8dAs+oLrtdReO4Md7J/DLg1N41XnDTVubFchz1+FbfzgGj9OGG86hMrVaIWEn2oUrNvdhbZcX//HUaLOXUjMk7ho8eSyCH+wZw5suWo8ev/GuSoIg2hebjeFNF63Do0dmV90mJxJ3mZloGl9+5DncfNdjuOUbT2B9jw8fa9GdoQRBrBxvv2IDun1OfPEnByu6wLYyJO4AZmNpvOGuP+Kff38UuYLUye97t12GoMXmXwRBtA8hjxOfvOEsPHEsgv/98CHT08eaTccnVBeTWbzjW09icimF/7j10rr7whAE0X68YecIDkws4huPHsNoJIHPvXqbpWlrK0lHi/tzk0v4yH178eJMDHe/bScJO0EQmjDG8LlXb8Pabi/+9hfP4+ov/xbXnzOIV24fwpVb+lryKt+UuDPGrgfwDwDsAO7hnP912c/dAL4N4EIAcwDeyDk/3til1g/nHKcWU/jji3N45NlJ/Pq5KfT6XfjG2y/C1SYbWxEE0ZkwxnDr1ZtxwzlD+Jc/Hsf9u8fwk2dOwW5j2DIQwLbhMM4aCmJdjw/run1Y1+Ntquizav4RY8wO4HkAfwJgDMBTAN7MOT+ous9/B3Au5/w2xtibALyWc/5Go+Pu3LmT79q1q971g3OOVLaAWDqHeDqHeCaHpWQOM7E0ZqLSv+loCifmEnh+KopoSpqC0xdw480Xr8M7Lt/Qsl0XCYJoXXL5AvacXMCjL8xg//gi9o8vYTZWOjDe77KjN+BGb8CFXr9bnvzkwhWn9+KqLdYCSsbYbs75zmr3MxO5XwzgCOf8qHzg+wDcBOCg6j43Afic/PX9AL7OGGN8GTIPvzs8jS/85CAS6bwi5kYJbKedoT/gxkiPDzftGMaWgSB2bujGWYMh2KgemyAIizjsNly8sQcXyyMSOedYTGYxGkniZCSB0fkEppfSmIunEYlnML6QxP5xaWymw8Ysi7vp9Zm4z1oA6gr+MQCX6N2Hc55jjC0C6AXQ8AbJIa8TZw2G4Hfb4Xc74Hc54Hc7EJC/97kcCHkc6A+60R90I+x10uBhgiCWHcYYunxSZL59RH8aG+fcMCBtFCuaUGWM3QrgVgBYv369pWNcsL4bF7yltQckEwRB6MEYg30F4k0zde7jANapvh+Rb9O8D2PMASAMKbFaAuf8bs75Ts75zv5+SmASBEEsF2bE/SkAWxhjGxljLgBvAvBg2X0eBPB2+evXA/jNcvjtBEEQhDmq2jKyh/4BAD+HVAr5Tc75AcbYFwDs4pw/COAbAO5ljB0BEIF0AiAIgiCahCnPnXP+MICHy277S9XXKQBvaOzSCIIgCKtQbxmCIIg2hMSdIAiiDSFxJwiCaENI3AmCINqQqr1llu2BGZsBcKIpD15JH5ZhN+0yQutdXlbTelfTWgFabyM4jXNedaNQ08S9lWCM7TLTiKdVoPUuL6tpvatprQCtdyUhW4YgCKINIXEnCIJoQ0jcJe5u9gJqhNa7vKym9a6mtQK03hWDPHeCIIg2hCJ3giCINqSjxZ0x9gbG2AHGWIExtlN1+58wxnYzxvbL/7+0meuU16S5Vvlnn2SMHWGMHWaMXdesNRrBGNvBGHucMbaXMbaLMXZxs9dkBGPsg4yx5+Tn/MvNXo8ZGGMfY4xxxlhfs9diBGPsb+Tn9hnG2I8YY13NXpMWjLHr5c/UEcbYJ5q9nlrpaHEH8CyAPwXw+7LbZwG8inO+HVIr43tXemEaaK6VMXY2pC6c2wBcD+Cf5Lm3rcaXAXyec74DwF/K37ckjLGXQBodeR7nfBuArzR5SVVhjK0D8HIAJ5u9FhP8EsA5nPNzIc1n/mST11OB/Bn6RwA3ADgbwJvlz9qqoaPFnXN+iHN+WOP2pznnE/K3BwB4GWNNnaKtt1ZIInQf5zzNOT8G4AikubetBgcQkr8OA5gwuG+zeR+Av+acpwGAcz7d5PWY4e8A/E9Iz3NLwzn/Bec8J3/7OKQBQK2GMjuac54BIGZHrxo6WtxN8joAe8QHvQXRmnG7tklrMeIjAP6GMTYKKRJuuWhNxRkArmKMPcEY+0/G2EXNXpARjLGbAIxzzvc1ey0WeBeAnzV7ERqsls+VLis6Q7UZMMZ+BWBQ40ef5pw/UOV3twH4P5Aud5edetbaChitH8DLAHyUc/4DxtjNkAa8/LeVXJ+aKmt1AOgBcCmAiwB8jzG2qZnTxaqs91NYofeoWcy8lxljnwaQA/BvK7m2TqHtxZ1zbklAGGMjAH4E4G2c8xcbuyptLK7VzIzbFcFo/YyxbwP4sPzt9wHcsyKL0qHKWt8H4IeymD/JGCtA6jEys1LrK0dvvYyx7QA2AtjHGAOk138PY+xizvnkCi6xhGrvZcbYOwDcCOBlLTqSs2U+V1YhW0YDOXv/UwCf4Jz/odnrqcKDAN7EGHMzxjYC2ALgySavSYsJANfIX78UwAtNXEs1fgzgJQDAGDsDgAut1zwKAMA53885H+Ccb+Ccb4BkH1zQTGGvBmPsekj5gVdzzhPNXo8OZmZHtzQdvYmJMfZaAF8D0A9gAcBezvl1jLHPQPKE1QL08mYm1vTWKv/s05C8yxyAj3DOW87DZIxdCeAfIF0tpgD8d8757uauShv5w/xNADsAZAD8Bef8N81dlTkYY8cB7OSct+TJCADkWctuAHPyTY9zzm9r4pI0YYy9AsDfozg7+n81eUk10dHiThAE0a6QLUMQBNGGkLgTBEG0ISTuBEEQbQiJO0EQRBtC4k4QBNGGkLgTBEG0ISTuBEEQbQiJO0EQRBvy/wG/V+Gg4n3ZdwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dos.plot(plt)\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "dos.append_fromPw(results=si_nscf_file,label='si_nscf',eta = 0.1)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXm8ZFV9Lb72mWq6c/eFpruBbmYZBKUjJEZDxFniSNR8xCkOwZcYY7Qd8lETefr0JT7ynjEBcfgZfcQp8REUkGjUCNFgGgIOgMxCz913rumM+/fH3vucfU6dU1W3hnvq3t7r87mfmk5V7apbZ5111ncilFIoKCgoKGwsaHkvQEFBQUFh8FDkrqCgoLABochdQUFBYQNCkbuCgoLCBoQidwUFBYUNCEXuCgoKChsQitwVFBQUNiAUuSsoKChsQChyV1BQUNiAMPJ6482bN9MdO3bk9fYKCgoK6xJ33nnnUUrpbKftciP3HTt2YM+ePXm9vYKCgsK6BCHkV91sp2wZBQUFhQ0IRe4KCgoKGxCK3BUUFBQ2IHLz3NPgui727t2LZrOZ91JGBsViEdu3b4dpmnkvRUFBYR1hpMh97969GB8fx44dO0AIyXs5uYNSirm5Oezduxc7d+7MezkKCgrrCCNlyzSbTWzatEkROwchBJs2bVJnMgoKCqvGSJE7AEXsCajvQ0FBoReMHLkrKCgcwzh8H/DIv+W9ig2BjuROCCkSQn5CCLmHEPILQsiHU7YpEEK+Sgh5iBByByFkxzAWq6CgsMHxdxcDX3xx3qvYEOhGudsAnkUpPR/ABQCeTwi5OLHNmwAsUEpPA/DXAP7nYJeZL9785jfj3nvvHchr7d69G+eccw527949kNdTUBgYvvVO4DOX5r0KBkrzXsG6R8dsGUopBVDlN03+l/zmXwLgL/j1fwTwKUII4c9d9/jsZz87sNe67rrrMD8/D13XB/aaCgoDwZ7P572CCG4DsMp5r2Jdo6tUSEKIDuBOAKcB+FtK6R2JTbYBeAIAKKUeIWQJwCYAR3td2Ie/+Qvcu3+516en4uytE/jz3zmn7Ta1Wg2vfOUrsXfvXvi+jw9+8IO45ppr8IlPfAK7du1q2d73fbzpTW/Cnj17QAjB7//+7+Od73wnHnroIVx55ZU4cuQIdF3H17/+dbzzne9EtVrFhRdeiPe///141ateNdDPp6AwELgNwCzlvIa6Ivc+0RW5U0p9ABcQQqYA/D9CyLmU0p+v9s0IIW8F8FYAOOmkk1b79DXBt7/9bWzduhU33XQTAGBpaQnXXHNN5vZ333039u3bh5//nH0di4uLAIDXvOY1eN/73oeXvexlaDabCIIAN954I8bGxnD33XcP/4MoKPSKxkL+5O7UgMrmfNewzrGqIiZK6SIh5PsAng9AJvd9AE4EsJcQYgCYBDCX8vzrAFwHALt27Wpr2XRS2MPCeeedh3e9611473vfi8suuwzPeMYz2m5/yimn4JFHHsHb3/52vOhFL8Jzn/tcrKysYN++fXjZy14GgFWZKiisG7iNvFfAlLtCX+gmW2aWK3YQQkoAngPg/sRmNwJ4Pb9+OYDvrVe//YwzzsBdd92F8847Dx/4wAdw1VVXtd1+enoa99xzDy655BJce+21ePOb37xGK1VQGBJ8J+8VjMYa1jm6Ue4nAPh77rtrAL5GKf0WIeQqAHsopTcC+ByALxFCHgIwD+DVQ1vxkLF//37MzMzgiiuuwNTUVMdg6tGjR2FZFl7xilfgzDPPxBVXXIHx8XFs374dN9xwA1760pfCtm34vo9yWXmICiOKIIiue3Z+6xDwvbxXsO7RTbbMTwE8JeX+D0nXmwB+d7BLywc/+9nPsHv3bmiaBtM0cc011+Dd73535vb79u3DG9/4RgR85/jYxz4GAPjSl76EP/iDP8CHPvQhmKaJr3/96zjllFPW5DMoKKwanmTF5KWa5ZP9wM1nDRsII9U4bBTwvOc9D8973vNi9/3gBz/I3P7888/HXXfd1XL/6aefju9973st91er1Zb7FBRyh+yz56Xc5YOKr8i9X6j2AwoKCnFy90eA3JVy7xtKua8CF110EWw7/sP/0pe+hPPOOy+nFSkoDAgx5Z6TLSO/r/Lc+4Yi91XgjjuStVsKCn3iyC+Bf70KePl1gFXJbx1y6mFuyl16X6Xc+4ayZRQU8sQt7wXu/xaw78581zESyl0id+W59w1F7goKeUKo9VrPnToGg5hyz4ncZUIPlC3TLxS5KyjkicI4u6wdyXcdIxFQVcp9kFDkrqCQJ4wCu8yb3D1plOMoBFSV5943FLl3gUH2c8/CbbfdhnPOOQcXXHABGo0R6O2hsDZwOalWD+e8DsmWyYtYVZ77QKGyZbrAIPu5Z+H666/H+9//flxxxRVDfy+FEYJbY5e5e+6SoMjL745lyyjPvV+MLrnf8j7g4M8G+5pbzgNe8PG2m6y2nzsAjI2N4R3veAe+9a1voVQq4Z//+Z9x/PHH4+tf/zo+/OEPQ9d1TE5O4oc//CF838d73/tefPvb34amaXjLW96CUqmEr33ta7j11ltxyy234Prrrx/s51YYXThcMVcP5ruOWEA1J2KN5bmrxmH9YnTJPSestp87wA4IF198MT760Y/iPe95Dz7zmc+EHSVvvfVWbNu2Lezzft111+Gxxx7D3XffDcMwMD8/j5mZGdx+++247LLLcPnllw/9MyqMEIRiPvQLRm6Gle86gNFQ7sqW6RujS+4dFPawsNp+7gBgWRYuu+wyAMCFF16I73znOwCApz/96XjDG96AV77ylXj5y18OAPjud7+LK6+8EobBvvqZmZkhfRKFdQGhUL0msPQEsOnUfNbh1gGjxIh9FDx3Zcv0jdEl95wg+rnffPPN+MAHPoBLL+08MNg0TRBCAAC6rsPz2A/z2muvxR133IGbbroJF154Ie68M+dCFYXRg0xijcX81uE22fQlt5EfsXqK3AcJlS2TwP79+1Eul3HFFVdg9+7dqR0fu8XDDz+Miy66CFdddRVmZ2fxxBNP4DnPeQ4+/elPhweA+fn5QS1dYT0i8IHKLLveyPG3IOam6mZ+nrvKlhkolHJPYLX93Nth9+7dePDBB0EpxaWXXorzzz8f5557Lh544AE8+clPhmmaeMtb3oI/+qM/GvCnUFg3CDygNMPy3O3BDoRfFdwaV+71/FSz/L4qz71vKHJPYLX93IF4j/bLL788DIp+4xvfaNnWMAxcffXVuPrqq2P3f+ELX+htwQrrG4EXVanmVTwEAPYKW0dzOT9iDfzo+jDPHpb3A+MnANxK3ahQtoyCQp4IvKi/jFwlutawV4DCBKAZ+St3og3vAHPw58DVTwL2fH44rz9CUMp9FVD93BUGjsADrDF2Pc/ZpfYK8/51Iz/PXZC7URzeAWbxcXb5wK3Ar71pOO8xIhg5cqeUhpkno4Y8+rlTea6kwsbDqCj35jJQnMxXuVNuyxjF4R1gyLFjVozUJy0Wi5ibm1OExkEpxdzcHIrFYt5LURgWYuSes3IvjAOamaPnzgndLA1vDZQNsoemD+f1Rwgjpdy3b9+OvXv34siRnDvkjRCKxSK2b9+e9zIUhoXAY50hNRPwcmoYRynL1Ak9d7/zc4YB8b5GYXipkCLd8hhQ8CNF7qZpYufOnXkvQ0Fh7RD4jFB1M7/cbqcKgDLlrhv5rSPwGOnq1vCsIWF9jaj1O0hs/MOXgsIoI/CYRaCZ+XndNk/lLYzlny2jGcP9LkSDNLLxbZmO5E4IOZEQ8n1CyL2EkF8QQt6Rss0lhJAlQsjd/O9Dw1mugsIGgyC0PBWzIDyzkr/nPuzvQvTPV547AMAD8C5K6V2EkHEAdxJCvkMpTU6vuI1Setngl6igsIERqtUcFbPDe8pb5fw991C5D4vchXLf+KZFx09IKT1AKb2LX18BcB+AbcNemILChkcQsOyNYVsRnRAq9/KIeO5D7G8jWhsfAxl5qzp8EUJ2AHgKgLSE718nhNxDCLmFEHJOxvPfSgjZQwjZozJiFI55CDLX9HxJNVTulRHx3PXhB1RpTmcna4iuyZ0QMgbgnwD8CaU02eHoLgAnU0rPB/A3AG5Iew1K6XWU0l2U0l2zs7O9rllBYWMgJPcRUu65eu5raMscAy2FuyJ3QogJRuzXU0pbumFRSpcppVV+/WYAJiFk80BXqqCw0RCSu8kVc17KnROeVeGqOWfPfZhpocKWCYLhvP4IoZtsGQLgcwDuo5RenbHNFr4dCCFP4687N8iFKihsOMjKPc+eLg5PhTTL+ebbh2mhQ7SGQnLf+Mq9m2yZpwN4LYCfEULu5vf9GYCTAIBSei2AywG8jRDiAWgAeDVVPQQUFNpDKOS889yFVRFmy+TtuStyHwQ6kjul9HYAbcu5KKWfAvCpQS1KQeGYQMxzHwFbJsxzz5nch2rL8M+qAqoKCgpDQ8yWyXG8nddkpK4POVOlE9akQlUod0XuCgoKw8KoKHfPZm12gZw9d3/4aaEiFVKRu4KCwtDQQu55KfcG68Qo1pKrchfxB5UK2S8UuSso5IVYEVOOillW7qPgua9FQFV57goKCkODIDDdzFm5NyXlnqPnTv3hxx8cpdwVFBSGjWRANTdyt9n0I6D/M4gv/x7w1St6e25YoTqk+IMYSiLea4NjpIZ1KCgcUwjz3Dmh5WbLNOOeO/UZEfYy0OKXN/e+jsBjB5lhWVRiKAlwTJC7Uu4KCnlB9tzz9LqTnru8trWE7LmDDr5FwONSv0NlyygoKAwNyfYDo+K5y2tbS8TIHYO3Zm76U3ZpVlRAVUFBYYhIpkLmastIee5A/2vpxfaQG4cNYg1JTO9gl2e9UCl3BQWFIaKl5W+eqZBCufdhy8gZLqIZ2Wog57kDg/8+3DpwyiWAXlCeu4JCWzSXgfp83qtYv/CTee552jLCc+/DlvEa0nV79c8PPDa4Wij3QROwvQIUJvJta7yGUOSu0DMWPvXbwF/uzHsZ6wMrh4B/+6t4kDCm3HPML5eVu96PcpeUtijzXw3kSUzJ1xsEnDpgjeX7Xa8hFLkr9Izp6kN5L2H94OZ3Ad//CPBESsaGZuZsyzQBg+e5i2BmL8QqE2ZPyl2aJwsM/vsIXN4czVABVQUFhQFB2Fe+E92XLGKiQT4TggblucfIvVflrg8voOq70tQrRe4KCh0RuD2otGMOvCBIJr3YsI4hpf91AqWD89z7Vu7JVMgBWyeByw4cRFO2jIJCN2jUVvJewvqB6EoItCp3YO3TIcXBZhCeu6yG+/Lch0TuvtyYTCl3BYWOaNYVuXeEKOV3pYySZCoksPbK3eYpi4XxaC1Aj557v+Q+5Dx3OWB7DCh31VtGoW84Tg878rGGTuQeEtoak45opFWY4GvpIw0xZss42du1e34sz32Itkw//XPWCZRyV+gbvtvDjnysIlO55+S5h+QulLvw3PvNlunDltH7OHvIfO2ABaw1k+XSA+z2BoYid4W+4auAamdQ3o3Ql76r5LAOYO099yYn9yJX7v147nJ6Yc8BVR3QLX57kOTOX0vT8+2fs4ZQtoxC3/C9nPKz1xOESvRSUiF1M79ujDaPlwzEc+9DuVMqDevg5O4P8IxQ/q7psdH2V5G7Qt/we/FXjzUIsktV7kOyIrpBfY5dlqb5Wvrx3PtQ7uLgJ1tUg/wuxGtpZnSGscGVe0dbhhByIiHk+4SQewkhvyCEvCNlG0II+SQh5CFCyE8JIU8dznIVRgVUCvwp5d4FBNnJpBcb1pFTtszKAXY5toWvJSfPPWZRDVm5DyvVcsTQjefuAXgXpfRsABcD+ENCyNmJbV4A4HT+91YA1wx0lQojB1cKogZKuQP3fIWNmMuCCKR6KcqdaMNRqwKNBeDwfemPLe8DKrOAwQl1UHnuqyXm2FmMIPdhKHepYOxYD6hSSg9QSu/i11cA3AdgW2KzlwD4ImX4DwBThJATBr5ahZGBJwVRA5UtA/y/P2Aj5ppL6Y9n2TKawdLxhtUJEQBu/GPg7y6OZ+oIHL4fmDk1up2X5z7sgq5AsmWOkYDqqrJlCCE7ADwFwB2Jh7YBeEK6vRetBwCFDQRPVu6DPH1e7xDZJ0mEyj0RUBVkOsxUyPtuZJeNhdY17b8LOPnXo/vy8txliyokd2XL9IOuyZ0QMgbgnwD8CaU04xfc8TXeSgjZQwjZc+TIkV5eQmFE4MrKXXnuEeyMal1Bdn7CcxdEsxapkMm1Le9nBDd7VnTfRvXc/bSaAkXuIISYYMR+PaX0Gymb7ANwonR7O78vBkrpdZTSXZTSXbOzs72sV2FE4LvRzh/kNR5uFJFG7pRGgyxkReu7EZmuRUA1ubbaUXZZ3hzdN7A89wHYMoMk39CWkcl9Y6dCdpMtQwB8DsB9lNKrMza7EcDreNbMxQCWKKUHBrhOhRGDbMtQFVCNkEbuvivluad47oCUCjlENZlcm0iDLM9E9+XVzz0MLg9LufPPI9oPyO+5QdFNnvvTAbwWwM8IIXfz+/4MwEkAQCm9FsDNAF4I4CEAdQBvHPxSFUYJcoYMVco9gp3iWMrj55L93IViXwvlniTLulDum6L7BtFbhuiRcvds9icqYDs9V04LHajnLjx9E9Dc+HtuUHQkd0rp7QibUWduQwH84aAWpTD6kLNlqAqoRkhT7q5kUSTz3IUtsxaee/L/JIK81lh0X1+eOydQqxJ9zv/7CuCx24C/yMgiCp8rFTENM1tGNwBf2TIKCpmQq1IDL1sBNX50HRr/fu1aLClfiGZUaemGMeUuk7srZcusQfuBpFUiWxUCfeW58+eYpei9HruNXYqS/07P1XSWGqqZQ8pzP3YCqqr9gEJPCKSAarvT59K/7GZXnn7lkFeUM3QT8Pz4MA6BmHKXvivfjch0LdoPJF9b/N9kch9EP3ez3Pqb8JqM9DOfK9kyAPPdB2rLyHnuSrkrKGTCj3nHG3sn6QptlTsnd2s8MWbPXWPPPUu5W9F9g5ihKtsyAmnfS9pz5QDzQJW7nOeuipgUFDLhO3JAdWPvJF1BEEWqcuf3FScTAdU1znNPvnbg8tYHenSfprH7+smWSVPuHcldKmIChqDcVZ67gkJXiFWldqM2O3mu6x3iO0gjMVEZOn58a567sGPWgnBaPHcnUuoyNLO/gKrsuQt0rdxFgNkaUj93Re4KCm1B5arULnaSDd1cTEz5AdJJrD7PLsdPSCh32ZYZkuferpmX78YtGQHd7C3fnkrZMkkLKO2MJrbOpC0zpICqrjx3BYW2iFWlZpC7bNc4dgfltp4hK8w0Env8x8zqmDwxodw9KaA6pGwZmdBTyT1NuRv9tR8wy60zVDsRdVK5a+YQbZljo4hJkbtCTwjSJgol4EnbOM0NTO4ycaUp94e+C5zxfOa5B26U0y3GygHDS4WU19ZC7k46ufeqmsOAarlVuXc6WKRmywzSllGNwxQUugJNVlqmwHWizBDH5or2p18D9t05zKWtPWTi8hLkTilQPQQcf27UM10QX8yW4SQ/cFumTUuATFumR2KVUyFXq9zlSUzA6g8w1SPtraS0PHeqbBkFhRbIhUskS7nb0Q7u2g2mWL/xFviffc7Q17emkEklqdwFoZpFQC/E75NtEVG4M+hUSD+RVy8jGLQtI+e5pwRv2z43JaDarS1DKfCJ04CvXtHF6yvlrqDQHtKOl0XusnJ37SbsOitB1zeaYgra2DJCyRslwODkLr47uXEYMPggIpCwZbrMlunHliE6+5yBF3+NTkTaT0BV9PN54JbOr6+rrpAKCm0hOkE6VAfJUHmuVJnpOnXUlubWZG1rDr9NQFWodKMQkbu4L0numjnkgGritbNsmZ5TIb34mDx5cEnXAVWZ3LtU7qJ1cTvEbBlVxKSgkAnRCbKJAkiGEvecSCl6dhP15XnpBTZQ3rsgCaPYqtzFbbMU2TKCtJLZKoOuypTXJr9veDvDltGN3lIhRYBYHMSai9Jjncg9pYip2wNMfb7zNqntBxS5Kyi0gu8sNrEyT299qXOk7zTQXInGvPl2Nb6xvQLc+PbuVNioQRByYSJFufOzF6MYBVTFfXLLX2D4nnvytbOyZXpdBw0Syl3qBNnpYJHquXdL7l2cEcoHD0XuCgrZEDnsDilAoxkBVZncXRt2LSL3lYX4mMW52/8/4K4v4sitfzWE1Q4ZggiLE90p91hANem5DzMVMknu7bJlesgxb1HuErmvOhVyFbZMfTW2jK48dwWFtuA7nkMKmQFVX5rW5DsNuLXoNL22GCf3Bw4xf/aR/SM6WzcIAKeW/pis3H0nTtBpyj0roNprlkrbdbexZbKyZXo9yISeewq5d+u5kx6KmLpS7vwsiZBjZhKTIneF3hB4cKkOnxiZnnsg5VUHbhN+PSL3+lJcbRGxw2FEvfjbPgH8j61A0k4CIpIQ04bkXPd2yl3OcweGlC3TJhUys7dMHxWqRI8OYn0pd6v7A4xM7lnPifXOV7aMwqiDUjQe+EE+p5eBCw86KPRMW0ZW7oHbRNCMphTZK8lTaT7sSxSzjBru4ANHlve3PiYrdyBuzYTKvcDUOxD58nL7AWBI2TJ8bUTPsGUGmQoZtFHunTz3tCKmbrNlJHJPFpGFr+9Hn1XZMgqjjkf+/R9R+oeX4P5b/nbN35v4jNyZcs8id2kUn2vH5ou61XiGgwa2c5NRVe5CbaeN0ZM9dyAeVA3JvRQp2i+/Gjj487iaBIaTLSNez6p03zis14NM6Ln3o9x7KGKSlXtW90nflVo9KOWuMOI4/NBdAICFvQ+u+XuTwIVLDAREh5axk8SGaHtNEHsFVcrUa1BLpK8JEhxV5S7OLOyUWaBClYpZpHKZv8j1lytUAWD/f6UUMRVaC436hSBVs9x9tsxqVHPsvQbgufdSxNQNuad14FTKXWFUQUWuOMnh3xh48GEgIAa0DM+dyuTuNqG7K5gnU2hSE7SxkNiWkRoZVXIPM0CWWx8LJHUMxKctxSpUJZW8cpAdyGRyNQqtPVn6Rajcy2tgy7TLlumhQrXrPHfJ4sskd8kC0zQARCl3hREGT6MjdIjTezKgcc89IHomucc6R/oOdLeKplbBEhmDJhe4ACD+iJO7OKVPNt8CJALtoNyLU9H9K9y7l5W7UYgfGAYBsTYzxZYJBmzLUL9NnnsPRUy+012xW32O9coHsg8IvtSBU7yPIneFUQURVkEOp5ck8OBzW0ZHli0jkZzXxJg7h5oxhSoZh+HEyV0o1iz/fmSQZlcIQimMs8ss5V6eAXa9id0WgVmZXPXCYHuYy2tLVe5ZvWV69P4Dv48K1aTn3mULZN9lB5HxLdHtrNePFYwpclcYZYT50nmQuwsfBqjWzpaJdrSp+mM4xX0Qy6UT0dAnYDlx75r4Tf66o7rD8QNpmicuvv9QuUsE7TbZcwVZXXY1cOJFwOIT7LZZirYdinLnazFXYcv001uG6K3KvZuZrOK5hERrkNefBYenppY38e2zyD0RvNYM5bkTQj5PCDlMCPl5xuOXEEKWCCF3878PDX6ZCqkQVkYOGSZa4MIjJlPuWWqb75g2NXDa4r8DAI7u/B00zEkUvbh3Tfi2WcHZkUGaJx6mQgpyl5V7kxG4IC0AqMwCi4+z6yI9Ulwfmue+imyZnvu5c1876blbY9157rJtItbVaR3C9hJnTe1smVjaqa6UO4AvAHh+h21uo5RewP+u6n9ZCl1BrnRcY+jUhU/M9sqdr6+GEnSe6vjMZz4LrjmBShBPKdRCz33t4werQjtbxsogd5nAAUbuDv/8MeVuDV65y6PvZOKjtE2Faq+2TEq2jFHsLkAr/PpwDV0qd2F7iRqDrO2TBw9F7gCl9IcAumi7prDWINzT1gZdst4FtMCFp5mgxICGjNNbvqPVSRkAaw88NTEJmBVYNG5vhMp9VHu9C+WdRh7JIqZkQFUmcICRu0BSufs2cO+NwH3f6n/N8nqTnrs8di6Jnm0ZTtAiK8its8/Uzev5Kb3tgdUr97a2jPLce8GvE0LuIYTcQgg5J2sjQshbCSF7CCF7jhwZ0R4i6wiCEPMYF6YFXLkTPXv4BvfcbY2R2zLGYBo6ArOEEo0rVD3gyn1Udzjhz7bLlkm1ZRqtyn3suOi6KT2mW+z1v/Za4Kuv6X/N8tqssUQrAid6zyR0ni2z2rbMgcf8dTmf3yx116smrbc90PmgIJS7yERqF1DVE+Q+qkJiQBgEud8F4GRK6fkA/gbADVkbUkqvo5TuopTump2dzdpMoUuQgGeY5GTLBJoBqpnQ2yh3mxrwNEYgNY2Tn1mGRbywJzwAaPyzZPr3eSPswZ4WUE3muXdS7puj64ZsyxTjB4ZmSsHUahGmQpbjqYXi86T2lulxWHdyWAcQkXs32TIx5S489w5rEMq92MGWSZ4ZaLoKqHYCpXSZUlrl128GYBJCNnd4msIAoIVByLW3ZXTqISAmqKbDyCL3wIEHAx5hO2pdIncAsBtV4D8/Bzz0r9ADYcuMKrm78cu0x7JSIY1CfHvZlpGVu5FQ0SKjph8EbjyDRRCa38aW6dbvbnkvTqCaFh0gjFJ3HR7TWjF0s4bQcxcB1S4ahwHHhC1jdN6kPQghWwAcopRSQsjTwA4YG3Se2mghtDJyIESDugj4VJss5U581qLA1xi5NQ2mrojFPfilQyje9Kfs9YwzAYwyuXM13ksRk5FU7pItk1TuMpaeALac29t65bXpVpywdUOyZdqR+ypFg2gcBrADmuOyg5dnd2HL+PHe9t3aMu5qAqqK3GMghHwZwCUANhNC9gL4cwAmAFBKrwVwOYC3EUI8AA0Ar6Z0I81QG10IxZ5H+iAjdwuUmJmeO+HK3ee2jGMwdaUXmH2xsPeXmAhfb4RtmSCIiCArW4bo7FRft1qzZcoz8e1lW0b49EDcqwYGo9xFLrsu95Ivt/fc+7JlpMZfADt40WD4tkwYUG1jyyQ992Od3Cmlv9fh8U8B+NTAVqTQNfK0Mgx4oJoB6Nm2DPFduBK5+xajco0r96NPPICTxeuNNLlLxJSl3AVxGMU4wbgpAdXSdHQ9puIT5N7NEIpOEHZEsuIzzJZJC6gKS2S1yt2LK3eAKXff6Wz++duIAAAgAElEQVSv+H3aMqHnnmXLKM9dYR0hJPdclLuHQC+0DaiSwIVPWNAVAAKe0aAXmXJvHnk03NYU5J7l3+eJdnNIAZ6JIZRqosrUrUV2jYBc0BTz3BPknui/0xNE58ekjx4GVFP0XbfVoUnIBCq+D7Pc3eCNtHmyQBe2jFDunWyZFM990O2VRwx9e+4K+UHnBT/5KHcXVDNBiA6DBKBBAKLFtYIWOPCIGSl7Tu4Gt2WslcfDbcdZTH4089xjDdDSAqqOpFiLcXXv1KIsGhmv+cdWZd5C7oPIlvGY3ROzZdDelhFnGmlnKe0Q+CnKvQzoy9ndGmPPTatQXWURU2aFaqJgaxjtlUcMitzXMQz+Q87DyjAp8zAp35kD34OuxYlCtCgwfDa8glSY92yVmJKdaEZTjSpgO+lI2jJp+eGxxyXiSHruWeR++nNa75M9d2ssfTDIahEq94SH7bexZbhtBjdjZmwWYp67IHeRLdPhQJWZLbPaVMisPPeUCtgNrtyVLbOOIXxqkoOVYcADlbIwvJSeKFrgwCcmTI+pcnNqK7ssseDXcf7hlNcdQeXeidxlS0FW7oHPiD5py2RBtmjGjhtMKwLfYQSuJTzsULmn6Dueqgqn3vpYO8Q8d37QsCrd2zLJkYNAuhKnFPjeR4AD9zDlromDF+neljEKqz8zWWdQ5L6OYSIf5U4DHxbxQTQThO8wntu6U+mUKfc7Zy4DAFjHn8UuS0zJTpNWZZrVPjhXyAovy5YJA6qS5+5w5Zum3NMg93uvHBep0n6Qmi2D9raMWK/bB7nLyl1Ovcx8blJZt7FljtwP/PCvgBv+MCoSI6R9w7PkwUMvKOWuMLowqSD3tVW74eBr3QLhys/3WklZD1wExMTTX/HH+Mun/ivOPZt1piiWxzNfezSVu6TwOtkysnIPyb3c3fvIKZNWuXflHvjAbf8LqB6JlHtYxJRI6UymXwIRuTurtWVkgubZ0CKg2rG3jJvw3EVQN+VgX+VnfIu/ird3aGe1tGTjmK2e+03vBn7wP9uvcx1Bee7rGCbySR90HZv9cAwrPH3OUu6+PoEds2N4z4t3hfcXSpFNsZ9swVZ6EABrDWyMpHIXny3jtD9myxSiHuMiICqCfZ1QksjdKAFeq23VFQ7cDfzrVcC+uyRyT9gy4gCUrIoFIlumJ+XOCVqkGZpl7rl3ky0jneHobWwZsS63wZW7RO5ZBxHPjges02yZ//wMu7zkve3Xuk6glPs6RjFMH1x7cmdvXAhJI/Badyo94FWsCZTLFQSUpQPOFbaG99dRGlHlLrUXyGo/oMvZMlxxV9lBK5wS1AkF6YzGLHbOMMlClTflq891YcsMUrlLtoxQxYWxLm2ZjMZhac8T6wpc3nmyFD0nbXtK2Xrkz3oMBFSVcl+HeOCbV0N7/Ec4jTBSX2tbRpA7MUxowpZJ2VEMePBTyN00dNjQUYCHWuUkoHkXAKBOSpjGCtsZ5VzwvBG2zU0ZeCEej+W5821WBLmf0N37EAI8czew+Uzgke/3HvATKZaUsrUVxqP1/fhvgQe/AxzPm7emtR8QB5nVpmLKBG1LE5K6sWVaKlTb2DLyQceTlXuG5y7uk89SkqmQG7CgSSn3NcB9t1yLox85E747mOj8GXd+GKcd+U54e62tDI8H+jTJlvHdFHKnDqiWctoPoMAPTN7kzvC+psbsgLSzgFwhSNaqpJNHMltGKO6VA+xy7Pju3+tZHwCe/Lv8DKBH5S6eF7itvWUe/Bfgx5+SbJkU5W6WWBXt8v7Wx7JAKR+4ocfvL0232jKUtgaLAy+eudONLQOwjB6h3LPUuC+daQrIB+Hka24QKHJfA2y94yPY7B3E0QO/Gsrrr7WV4blCuRegGW2UO/XC6tQsmJt2hNcFubuyYv3PzwIHftrnivtEu1F14nFBTMUJwOYjBKuHAbMS7x/TLcxS79kycipmMs9dQMQF0mwZAJjcDizt7f49hfIV6vu3uG+9+QxOutL3dvvVwMe2x1MtV2XLVKPr1UNxzz1te0HiRtKWSbSJSH6WdQ5F7muAMV592awOoOIwgSotwVjjgKrv8GHWhgXCFZaforZNuCwXvg1K05Ef7ejM6w3PAmpHgZveBffvXzqIZfeO0JYZSycPecpPcZLZGUHALktTrdt3g36GZQuiogH3mq3WdExhGaUFVAFg5lTgyH3dv6ewBoVyP/9VwF8ssQygZKDz+x9jt8WZDZCezQJEiv/og5FNJB8UlvdHAWDdSm8GFir3FFtG9DiMWT0bI/9dkfsawAf7wTdrgyf3Gln7IKTIjNEMK8xzT7NSTHigaZ6uhMJ0FFB19XLs9Y/8ipGL2Zxb/VSgQSLmuXfIcy9OAaBsRqq93H2mTBJ6gRFmpyyTNAhy8iVbJllIVTvK3yeD3LfvYkO8Vw51956CVLN61ciTncQBoC5N7wz81jF4QPTdf2oXcM3T2XWZiN1alAqpZQRu0ywoI5EaKtsyq+2pM6JQ5L4G8PnX7AyB3BukvObkLmIHmqTcgxQSMqkLmnHaf4gyRbt522nhfbbFuiUKcl9cWggfo4Moxe8VWQFVQVaeExFHcZJdNpeA5nJUFr9aCPLppf+J8NzdWnTgMRO59vUO5H7c2exy/pHu3rMduSf7wwsybsjknhhgTUg0e1V8z0u8BXKyLYKYdKVb2QFv8Xi4Jn5dEL9syyhyV3jsJ9/EA9/5bMftCC/o8JweA2QSPCe+s9t6BRqhCFKKiIYF4blrZhGaIWyZxA5BKSziZ5LHV55yPV6Pq7BpIiIdvzwbe61aNfJW64tdKshhIGnLUMqU5v86E7j9r3khDScYmdzt5Xh642ogDoq9WARhEVU9yuRJayesW9lZSWLOa7Vb5Z7w3GUkg6OCxGXvPNkeQDzPd1uDncm2CKEtk5FPn6bcxfcr/rcb0JZRqZB9YPPNb8UY6sBz3tx2Ox0BAICm9F9ZLRy7HvunNYwJwGNByIKxNv/OgCtr3bAQ8ANXkLArPKfJ1plB7u946W8CL/3N2H1kjJG7UO6OHe3Ey0cPoLLl9EEsf/VITloKPNbUq3oI+O5fsGyYpHJvLDLlPr2z5eW6gpEgn9VAqFBBirrZSuL1+exgKhBl+NS6HGQfKne99bFYjn2FDdEGonRJ8fykhSfIPXnW5tSA8a3ACs/mEfEE3Uzvh5OW0y+3QH70tigILm+/zqGUex8YA/sh2Y1q2+1MwlRNMABFYDfj6t8usKk+nrN2P0iffw7dKoDwHTd55hAVOrUPqALA//FejjuCs2BVpvhrsc8SSGc6K8sD6G3eK8JUSK4QfSeeAy4PwRYB1OYSI6WebZkBKHe3zoKyaf+D+lx6jrtAeRMAsgrl3s5zT3R4FOQuq+VkbxkgsmWS5O7Wgclt0W1x0M2yZdKqccX3++htwN9fBnzvo9FjitwVBBaPdJcPPBDl3oz7jaGVkZJnPiwI8jXMAnRDtPyNfzY3zKhpow45nvK6v8Sdz7oehsm2FZ4+lXxQrzkKnvtYdFsmHK+Z4rkv9h9Qld97NZDz42mQTu5ONT3HXUDT2YFKDnq2Qzeeu7Bl0sg92VsGiIqSZFVNKVu73GQtVO5ZnntKnru4vryPXcqZQRvEllHkPgA0VhY6b4TBKHeXK/dFWsEROgFMsGwTZxAdBLsE5eSrGwVofMelCa/TaXDysxLDoVPwzDNm8d8uOY0VRUHqUyOTe2OUyD1BOL4tee6cdKqHGen3Su5GIuC3GiSfI8hVBEnD+zucVZVmgEZ3v+325J5sfcBJvsVzT9oyfFqSbN+4dWa9yKmd4npWG98wz10OqPL3aqQcvJRyVxBoVrvcAQbwo3G5VfGLXR/B0h/8F3SLZR60BDT7hLe4D0dv/khqgEood920QnKXs2Xmb/kIxr7wLLZNoct2t/z1AEm5Szuq32xvfQ0VYtKS7IMnrQJZuROddSwEerdlQuXeC7knDvTCMnr19cDv/0u8yVk7lGfSyS8NbQOqgtw5qYuDtiD3wGdnGMn1pNkydpURvFWJPkdI7sV0cs+qUAWAWsqc2g2i3FVAtUeEAUMATjXbD/ZdB+JkcxC2jNdkPn+lXMFpWzdjvk35fz94+B/ejTMP34wnTrgYJz7l2bHHxBmIYRXg2uxboNKBa+aOvwqv66uozhTKXQRsiXQ2Ejg5krvHm07JJNVcjm8jCJQQZmcscHLvW7n3ElDNIPeZU9hfYYwp8k7kXpqJFxq1Q7uAqkh9dBu89QAPegpbRpBpS0BV2DIyua+wg4JZZgeSwE3YMhIx/+QzrFhr9iy+jkL8tYH0gLFS7sc2asuRonEby5nbNeRg6wC60Il0SqPAdtgwFXFAfWsE7EW2Uy8tHG15TChq0ypC5ySUtGUEjNIqyJ177uFUJ78Bm5oIKAG1V9mhcJAIOytK+dp24n8uE0dpJn/lLlscZuLsyeLpmZ0mRJUHZMuIg4vb4GcVoiqU7xtpyhqQbBmJ3J2VaHShOEMRXTeT/WJufjdw2yeiGIQpWYTielrAWJH7sY2aVGDj17OVe7MWkQAdwI9GkLvOMzcIJ9eB2zI6+/E7Kynkzg8kZrESdoXMInez2D2565wghe1DPBtNYqEBK+7PrjV83gtc9o5byF0ijtL0AJS7yJbp4f/qO0BlNrqdHBZS6JLcSzODCaiGk51q8WIhodzTujYCUYBU/q5rc+y+4gRw3uXsvqmT+fMzWjaI9Ei5kEtcTyP3DWLLdCR3QsjnCSGHCSE/z3icEEI+SQh5iBDyU0LIUwe/zNFDfSXy6oLkKbqEZi1SHWQA5C76uhg8UCkqRNOGZfSDgJ9K+/UU5caVULFUhsF9ctlzr9NIga2K3MVriSIpvwkbBTRQBOm1t/kgEHZWlMg9+T+Xe8iUpqNeK70WMRl9KvcxidyT1anCKus0/q88zQi5G7ITnjtJsWVk5R7r6Ji0ZRLkbpTYZ5GVe5jbPg5c9tfAm78XTbBK9osREC2QZeUuPrvosSPjGFLuXwDw/DaPvwDA6fzvrQCu6X9Zo41Hf3ITZv/5NdEdbUrjnbr02AB+NAFXIWaR/VB1bssMvk0u4e/XSqrUbSKgBJZZCG0hKllOdRINei60GamXhG7Glbvm23CIhSYpgHg5tmT17Hjb3DRbRqRAAvFxef3aMj1lyzhsBqtAktyFYu9GuQPdqfd2nrs8cFs+SIssmKzBIWJgibx/iTbEhXH2t/3C6DFDOvgGQXS/GJpiSAPIw+9EOhCIs6xjRblTSn8IoN1/9yUAvkgZ/gPAFCGky+kE6xOVW/4YkzTauXccvBVP3Pnt1g0pBT1yf3hzEMo94MEyU9gyYRHRYNUGEVW1KYqZ+DZsmNB0DbpodSt1/fOlOH1pYnPX75kkd91vwiMWHFKA3mtv80EgOYc0LVtGJvfSdHS9MIme0E8qpG93Z8t0CnaLz9FNxow4uKelV8pj+4Ryt8Ylzz0lVRGIplrZK5HtJZN7EoK8PTt+hrByiK1BrtJNm2srlP0GUe6DyJbZBuAJ6fZefl+XYfb1B5NGO5xHNWzxDwDffBVwYbwx2H3f+iSedOeHwtuk0zSaLiDI3SqxH2fU22Wwyl33he8dkerh734S3oGfg3g+bGKhBISFR7LnXkC0c0xPSyq2A0wzXrij+zZcUoBPCIy8yd1IsWVKUqqgXFQjz0It9kjufQVUnbhNlFTo4nZS0SchzkC6CaqmtdUVCG0ZSbmPzUbWVpYtY5Yi5T6xFZh/uD25y4VfVFLuKwfilgzQGmQGIktpg5D7mgZUCSFvJYTsIYTsOXKky54VI4hFLVJmh7TZzO3se2+K3R6EcqfcJrEKTKVEanewp5J6wAdySKR63O0fxNaHvwriNeGA7Yji4CJnAhX4we/u4BQYRsppegYMi2feCOUe2HC1AhytCCMYBeUu2zIrwIzUN0ZW6zKx6j3qp34CqnLFLABUEmdP4nc4dhzaYjW2TKjcU1oayMpd+Ozlza0B1RZbhpO7U2VrJVqUmpmq3MV31kwM9DjceiAzrCj4KypmxZnMsWLLdIF9AE6Ubm/n97WAUnodpXQXpXTX7Gw2KY46XC3y7mwtUgR+IqjpFDZF16kObQDKXaQhWkWmPMLc8EEr94D73inZB0V3AS5hO7Fp8mwZHlCjgY8ScfDV0itx6PIbVvWe0VkAe2+T2vA1C65WgunnSO5eii3TWAirgwHET/PL3Z+tZCJ8r1USTRAwi0wmymTGjjj4dDqrCJV7N+SeMu1IwCgwAnUkch8/ngXmfS/6jC22TIltI9o4FMajdgFpWUhyPx65gVj1UKtyByL1vom3nRaBWKXcQ9wI4HU8a+ZiAEuU0vVtyVDadjiEGUQ7XDmI8q9rtXiQjUge32FtdiDkLlSFVUgEVAc8yd0I+JxUv5Xcy81DcEiBv3986IFrMxLevuV4PO/JJ6/uPS1O7qJ3TeDA0wrw9RKsYO3aK7QgzXOvz7HmWq/8IvCq6+Pbyyq+VySVexAAd32xbfA+XJt4/gs/AVz8h60dIX/zT4Gnvh4467L2r9WTck+xZQhhytltROQ+xnPTRc/5tOeaRVaQZa8wYrfGI4so7TuW/z+xNsE0ndzFAfmUS1iM4lkfAECOHXInhHwZwI8BnEkI2UsIeRMh5EpCyJV8k5sBPALgIQCfAfDfhrbaNUDgeQg+PI17r39v5jYWbeIIncRtF10Xm7foNOKFNkRSXU2tMhhy923Y1ISus39dWEQ0YOVuUG6NpHjdU94cPMLe1wgzFNj7N+vsdJh08nNTYPFWCiLzxqQOfL0IzyjBonmSezJbxmFqtjQDnP0S4EkJkhwEuWsGGNHw39AT/wHc+HbgOx9q+7Qwz9soAE97C/D8/9G6zcQJwIs/2TmTxyqzIOVqlHtWp8niJHsdYZeM85bCTi06gLWQe5l9/uZSlB0DsAKttO84DKg2W+si0n6PwpaZ2Arsfgg456XZ/WnWIToagpTS3+vwOAXwhwNbUc5YmDuATaA4+6FPA/jL1G0K1MbDM8/EM17wKtwzuw3T33wJCsRrIXedq97v7fo7bL37k9CCASgCrwkbJsTJr26KVMTBqg2TMoI1gtYf+iYs4LDGEqJ0XYNHNYDPcRXtj0laNkIHiJx5QRQWteHrBVC9FPr4ucBLBFRrR9mZSnlT+vZTJzMCuvjK9Me7ASFxohHqefHx9s/LUsG9ojQDpNU6JJEVFBUY38JyygXpjknkHtoyCUtHkHV9jpM798THjgO0FF0qt2xoGeiRotzFWmWLRy8cO8r9WINd71zmXoANylOzzt/1m/j5RayXSrIdr+7beMA4A8+67DXwNRM67V9di9xvAYNPfh90KqRJhTXCDlA0OYxDWoMPPTyDcTi564UeyF3X4FA9PAuwqAOqF0DNMkoYAVtGHLCW9rLLLG+9PAPsfhB4zn/v730NiWhElkknVRn2Li+2365bdNs8rJ0tAwDjJ3ByrzH/vcyDvPZK9nNlQi6MR1k+ot1AEnKGUXJ6U7v4gvx6hrVhlLsi9wQ8u3OZe4E6oNIPT7QCcBPkbgRNeNybDogJfQC2DPFtuIhOfSPlPlhbxuJK2eTkbifOSnw9Ig8PWpjmKc5e9E7VjykghMBFNOTYggNqFBGYZVjwUjtUrgm8JgvuiQCcUM9Zyh1gtkHWCLtuoUvKXcwN7aQq2wU2e0FpustUyA5nDNM72DzW2lFG0kKFO7Vs1S8foKyxyJbZfGb6e8h57k5CpKWR++Yz2OWklA+iFwbSA2oUoMg9AafRvhLSc2w2G1Qm94Ig9/hzjcCGx9WEr1kDU+6upJqjxl0DJndw9cxtGeGlC3haRB4+jEi58wOcXly9cgcAD0Z4oChQB4FeDP373AZ2ODXuP3NrRgxqLg0gK6YdZOUubIY2gX4Akec+MFtmusuAagdyP/W3maK+94a4Cm9ny8gCoTTNDhAAsPWC9PcIbZlmpNxF4FauQxB44V8Cz/0IcPy50X262VttwQhCkXsCSfWdRFM8LgVoTE5kXsLnMwMbPifBQLNgdEHu7uEHYn3Mk9AT5G6aKeROKR79ym4sPryn4/tlwaJO7DKp3D0zyjP2SKTcPf79GKvo4y7DJQaI74B6DgwSAEYRhO/kjXpO5O7Wo4O5VQEWHmPXx4aczqtbEVkL5S4X56QhN1umQ0B152+xM53GAstMCcm9KjX2SvjisdqBaeDit7FMnwteg1QIZW9Lr9ku7XPqJOA33h737zdQQFWRewJuh9ayUcAw+iEKIvPtOLlb1A67Kwaa2ZHcXbsO8+9+DQ996hWZ22iBA0+TlLsZz1YBgMVDj2Pn/deB/l/WNe/A3f+CIx97MtyV7grHqO+xsxMAFtgP3U4c9AJLIneYIPz9Pf79mcXeyN3nyt0R/wezCMK/X7s+hM6QlAK/vCXei0Tg4M+Bh7/PiEtYMtZYFBSsdCgC6hcy0Qiy6tqWGWBAtbEALD7RfjvfYRWeab1lAEb657ycr60YqXKnGqnsZNVoktwntgLP/vPstgmi1UNziR0MNTOqOu02gylrVN86hCL3BGSCDny/5XG7zh7XJOUuWgEETpLcHQR698p9mfdOP33xtsxtjMCOBTNNSwSRotdeXjgMAJimzCud+97fYNb+Fe6/89/avr+ALdlLRe69O4kh4FSqEHSJGebD+/w7sFbRx12GRzi5i0HgRglakb2X3WYoSs946LvAl18N3H5162PXPh340kvZdRFMFaRkjaf3JxkkZKIRHnKn1sehLTMgz708wzKD/ve5wK9+lL2dCDq3wwU88W7bhRK519ifbrVW8ybJvRNEaqe9zEfxlSPFPt1lzYVS7hsXMkE3U+Z2us3WbBCzIMg9nhNuwQblp8dUt2CifUCwm1mseuDAl3Zc02gldztBxC6f0uR3OfDC4YVIS7SCInFBg6D1uVL6mEMK0LlPGfDtCj0qdw8mtMCJnSGRCTbp3j90b0+v2RaiHazU4I0tJLGDi4O5uBy2JQPEiUbYMp0C/uG80AHZMnJc4cgvs7fz3c5nC9suZGP+Lv1g3HN36+l56DKhJ1sopEE3mfoXyt2ssPc69VL23t1ABVQ3LmTl3qi1knsYMJTI3eKee7KDYoE6ErmbMNH+R9Nc6axMTerAl2wZTeSZS5k4boLcAz6Vx6umzItMgQgqrxC2A7p2oyVYrJcicndJAbpImexTuTvEYnEF/n7EKmPyOKa6Tvrhu1i2xSAhClmSO/S+O+O3Q8XOL0We9jARU+6J0XRZyCrl7xVyd8m0kXTh+3ah3AHgpIuYt64b7AAkPPe07KqsTpvtUJwEmouRcj/5N4DXfqP7Bm4bKKCqZqgmQKX82DSP122IgGFE7oUS+2EG0uxKGgQoEjdsVUp1KywMyoJd60zuBnXCIK2ABz1G7l7CHxcZecEqlXtdHwf8Q2g2qvAS8QR9Iurq7GpW2GhMkHux0hu5u1oRRtAMD6KGVcbmLdvCx1d+dQ/Gz760p9dOhbAxkuR+f7zpW6gsheKsrIVyL0ZnFuJ36dZYfCCtiAdInzrUD7ZJs3eSPexleF2SuwyrwpV7LX29hADP/zizhbpNKy1OMuXuu719B0Yh+s7XOZRyT4BK1oqdkp0hPGVdmjBUFBaERO6hb82LjKAXOtoyThfkblIHNLETuTDCgCYQBTUFCM3uzZ4G8fyGztS506iFn/tOPAkAMHHi2eH2vlYMK1mp20BACYqFlIrALuBoRRh+M1TueqGMTWNFvNBmZfQH9z3a0+tmv6H4rhIpho/+G3DC+dFtUbAkDmprodxjqZDS/9Rtc5AWj/VQZ5CKymZg98Ms06VdXxvfyc6UyYJViTJbsuIXF7+NZbR0i+Ikq6gVc1ZXCxVQ3cCQCNBJ8dyFgrWkVD/DNFllpSc/l9sKIr1LN2GQAIGXTfB+YynzMQGLumGQNlwTMaJJOAD8ZvyMQ/S06ZbchXJ3TU7uzToo/9wrz/4rXLPzUzjl9IjcPa0QkjtxG2jAWlWrXxm+XoIZ2FJKZRmaRvCJ32fDwGrLA1ZVWQHKow8CO54R3Rbes8iXXitbJvTcpTOndr67OAgMSrkDjOCLUx3I3V59ELfAVbZbT++v3gumdwALj7IMn156/Biq/cDGhUTQXqN1JxKBRasU/zHakHKSAdhN0UCLkzsPfDopY+vC1+bk7tJsYjTBSvJjz4MeGwRCJZVHAz9sl0u6HHghhnC7PLXMadQQcHI5/6wz8LbXvxamHv10fL0Q5sMTrwGb9J6p4etFWLQZ1gwY/KxodjMLqAXJQqa7vwzc/r97fr+op7i0Qzt19r+UK1BFifrZLwFmTgWe9Du9v2e3MAqR/ytnYrXLmBHbDUq5CxTGW+fGynAb6f1b2qGyCagfjYrEBoHNp7O2wPOP9lZkphd666E/glDknoA8iDnpMwNRNk2S3B1ixXqfC89YC5W7IPc2wRq+83jIJvcCdUETlXwOsaBJ5CRn/NiNKiyfraVbcvdt9jmCAlOprl0LvfRCuZU0Ar0YTqfSvAaa6IPcjSJrGNYUs2LZTj8+NgabGnGCCQLghiuB7/55ep56NwhTDKX/tSjaKc8AT3oxK08XtszWC4A/vgs47qze3m810K2IaJxaFBRsR+5ujf3WsvLNe0Vxor3n3hO5z7IgrVMd3MHo5Kfz9dTYgO/VQgVUNy40aSiEl5KZEAYMk+QOK9biV3jGWkGMw2OEJ/qdp4LvPFnePA0CFIgL6PE0N5aKKL2utO56dQkF3i5X69KW8RPVfV6TpasxL71VYQVGEYWwRXA9NsBktQh4B0hxpmDx2EbB0DCPctwaOPSz6Hr9aOfJQmkQFofsY4ty+9IMcPnn8ztNNwrxCtWxLczG6GTLDCP/3hoHam0sMbex+mHg5c3sNc3iYNokA8D2X4t8816Uu6GU+4aF5jXRoHy6UYpyF751IZHq52iFsMUvwFEG87oAACAASURBVHxqANB5JSvhqWmend3dUHMYcRkkgJei8F2HP9eMK2OXxN+byBk/tRUUA36gSRm8kQafB4Y1rlZduw64TTRgQdNbfzLUKIa2jOHV+iJ3apZQpHb43YtMJEII6qQM3YnUY/WxKF3RXuxxPoxQwfKBXFbuujl4i6NbiAEXADuzEAevNOW+7y7g734dOHz/4PxrGYUxwGnjubuN1fv8lU3sNWtHBkfuuhH1oJk6qYfnWxtGuStyT8Dwm1giTIHQtNRBng1SKMTVsysV8gCRpSNSJglX7p7bpm+MG+20zZQeNyIDhyQKVFytEKYiAoDmybn6S2G7XL1Lcg/4QcQYY+TuOw1oXh02yUh1M4oowgEohek34PRB7jBLsIgXEphcDNUgJeiSwj7yyE/D6/PzPQZaQ1smQ7nnicIYb1/bZLEgkX6Zlut+9z8Ah+8FfnV770O528Eaa3/G4DVWXzhVlgqTBvldX/I+dimncXYLFVDduNCDJqo62zmS7QQAADwbhCTyjD1ixQg2Se4aV9uek02wlhcpo2SjLkBKr0zsRL6UrQIw31vAqS2hzG0Zo8tRdYE4O+HkTp06NLeOBtKVGTXK0AiF5zRgBQ14eh+2AFd/hLeZLUn58rZehuFH/xNHUuuNaudMo1Skkbus3POE6N9TPcguxczWtKwVOZtmGOsujHfw+nvx3GVyH5ByB1gPmz87AMycsvrn6gXWnC2v9tIDhCL3BEy/iabBd6oUj1rLyAbxEgTr84wToyiUO1O9bltyjwjGabYeWBxeVKUlfG9PL4Z91wHmewvYiwehEZbDnTZVKQ3CeipOzPLPUofuNWBrGcqM95mpryygEDTgGf2Tu1k/iCotomBGYSFXK8GSyN1vLMGmLLfabbQJ9rVDaMtUo3a69TZzOtcSokHWMj+IiYydpHK3q1EbYqB9n/l+1uI72X1XsloItINcCDa5LXu71YKQ3uMO4djI9W/NKHJPwKA2XL3MfHevlWCJ34SDVnvC0+PkHmbVcHLXeTGT38aWKUrDttP6you8ey3RFc/XizClEX6638QKZSpKVrdml8qdcs+9MsmUFXUaMPxapt0iWhHUl+dRoE0EfZA74YRaqe9HDeXYkHHXqMAKogOu5a1gP2Uq1euZ3KVWuiJ42VhgvvWgBl70ClENu7KfXY7zAipZQd/zFeBj24BHpWZzwyB3S2qnmwa3GRXsdYuZU6PrwifPG6JAcAM0D1PknoAV2PD1IpqkEEuLFDC8eqqC9bVCOJoOkMmdecY6t2X8NqmQpaCOeTCidJutO5HNyd1IKPdAL4aTkwDA9BtY0FimC105GPtsXYF/7vEZpqyo24DpN+BmkLtRZu/VWFlACQ0EfQT0zDFGTDPOPtS0+Ot4RiUMDgOA5VUxp7MDkN/rIA9ZBQsbrjGfvyUDtCr34hRTx7It89Ov8StShe0wWiOEk5NSvufAZ0p3tcpdbr422UPwcxjQW1tor1cock+gQJsI9CJsFGLetYDprqChtfZNkXO9gaiPS7HEFE9I7l62eq7QGpY0plyTgz8AaRBGcTx2f6AXUEB0YLGCBlYMRk56/RAAYJ6Odz9k2m2iSU2MlccQUALqMi/dzVDkZoWRu11dQBk2aB/ZJcUJRtZTdKnlew7MMko0+p+UgipWLJZBEqQcDLtCWv54fT5/SwaIOm8uc+VeGI/6sYTbSN/Rzmey+aQ7nzn4tYiziDTlLkRQL50o3/oD4I23tLb7zQthl1Wl3DccCrARGCXYWjGV3C2/BttII/dCmOsNRD1qihVB7nyQtZseiQ88F2Vio8pJ2UvJhxfknuy4SI1SmIoIAFbQRMNiCrjQYJ38lvRpFLscMk28Bmye9tiEBeI123rphTFGhM48HxzdB7mXpiI1Z+tJch9jmT/cGy8HNfiFGTSpGavK7RqUMkIXrQTEa4yKcheEuryPXYrxdLItI883vfANwPseB3Y8ffBrEf3704Kqgtx7aXmw9Smsc+OoQFR/b4Bcd0XuCZSoDWqW4ZAiDL+VYEt+Fa4x3nI/NYowJfUMtwaPamHKpMGHasidI2VU+SCKRpHbDCk59lnkDqOEMrFBeZVmgTbhmxNwqIExl7XIXTE2scEbnWZwgpF7kweNm8SC5tVRbOOll8b52cYca+pF+vB8J6ajQqS6FSdYYlWggTIiD3xU0IBWnkINRZBeyN1tAKCt5F6fH45vvVqIEXGLv2KXhXGm1OUqXXm+6QkXRCQ8aBTaeO6icnW1RUyjCBE36LKae5TRFbkTQp5PCPklIeQhQsj7Uh5/AyHkCCHkbv735sEvdfhgw689wKyw1rMpeeFlWoNvtf6IqVFEkcrk3kADhTBlUrfaB1QbS2wn9UoiQ6X1xxWI3O9ygtz5D1Jk4hTRBDXLqJESpn2W/90oHgeT+G2zdQQ0vwmHk3sdFehOlZF7hpdenmAkbCw+xi7HuxiskIHJySnWhA2AXd4Se4xyC8Kur4AKgiuOo4EiiNuDLSPIPFkcVJ/LP8cd4N45AY4+xG4XJ1jg8cFbo/vqc8BZlwGv+2dg06lZr9Q/rDaeuyD3wgYg93BCVEoa9DpDR3InhOgA/hbACwCcDeD3CCFnp2z6VUrpBfzvswNe55qgzlv8EqsMTy+25IXTIMA4rcJP+RFTgxXfiK6PrOgnyrYQaptm/GgaYlAHV5FpOfZR1WaS3JmibjZrAKUoURuBWUYTRUyCEZZX2cLfp3M+uO5F5F7Tx1F05tiZQQa5T05thk1NTFYfBACUJnon96JlRL1pJuLpcTqPNdSri2issIMhKU6hSUrQ3R52RkHmsnIPfFbiPwq2jG6yMwhnhdkFhQlg+9PYYzdcyc7C6nPAptOAUy4Z7lqEt5+WY9/cSMpdmu26ztGNcn8agIcopY9QSh0AXwHwkuEuKx/YNUHuFXhasSW7ZGX+MCzig4xvaXku4QFTUWikew00SRRgKlfYDz/ISCVr1ph3qo8zoklT7uLAUKrEdyIxrNtp1OA6DRgkADHLaGqSjcILYBq1Lsg9aMLlGUFNYwIzDvd8MwjPMnUcJptwsvcYW99kf+1wf0RYH3X/hHiFocatAbu2gtoyI3ejPAVbK8Vy+7tGSO5cubt1oLEIgI6GLQNEB57x41n+9kVXsuKcI79k6/edtVlru4DqhlTuPdh8I4ZuyH0bAHn0+V5+XxKvIIT8lBDyj4SQEweyujVGkyt3rTgG3yjBonGCnTvIPGVzenvLcwmfuBR2g/QbcCRyF4ScFfgTgzqKM/yrTa2OrcGnBMViPCVR4z9Iu7YcjQa0KmGPF5fqsHhBkl3vnA9u+E14fNqTY05iS8AGbutj2Yp80Ygem9yyo+N7tMO9F16FNzi7cdYF8cCgWRLkvoQmV+5GZQqOVoLp90DuQoXK+eNiCs8o2DIAI3WANQ0DWJHNU1/HCHWBe/FrQe7tAqobSbkfY+TeDb4JYAel9MkAvgPg79M2IoS8lRCyhxCy58iRNvMYc4IoEtILFQRGqSV1sHqY7UyVTa3HLsJPWxs19kM3/AYcKR/eKhRYNWXGj8atM0U9McsOHDRFIRG3zvzlROsDg6ciNlcW0KyJQqcKHN4GoErKMMtsx7NrncndpDY8vnbXmgrvL7SxW2oFRkKH6TQ2T/fX2+QdL/o1fGT3n+KkTfEArsHJ3WmswObDxK3KNFy9vDpyr88zBSpISa78DFsPjEAqJAAcxx3QmZ3RfZP893fgbna5FuSum8z+a6ac+W0o5S5iC8eGLbMPgMxm2/l9ISilc5SGTPhZAKmjximl11FKd1FKd83OrsEMym5BKUBpOFjaKI6BGiUUECf35gL72JPHn9zyEoJghRds+k14ida8DVKIdWyU4XNyH5vajBotgKSMUiNulMUiQ/SAaa7MRWcfhbGwx8uiNg2zGBFjJ5iBHa49KEXkXpzM/p9VN50HAJjTN8eqSnuBrhFsn27NzLHK7KDhNpbhchurND4DzyjFKlc74nPPBT7/vIiUQuVeGz3l/qQXM9I875XRfYLc9/8Xuxxbo31p7DhAKooLIQ6Sw8rUWUscY8r9PwGcTgjZSQixALwawI3yBoSQE6SbLwZw3+CWOHw8+H8uw96P74LDyd0sMXIvJVIHg6UDCCjBpi2ttoxVYUqvKRRl0FrR2UARWsb8S5H9UR6fRoOUQFKUg+FW0SStVaKlSabc7JV5OFyZ68UxNAuMoJbNzRIxdiZ3ixdyAQAdjxy4qRNOy3zO+FNejp8GO7Hn5Dd1fP1eUeTWlt+swuP9X0qTm+AbFRRpl+QeBMDcg8Chn0cqtDjFBpk71Si1cFQ895MuAt77GHD6s6P7RB8WQe6VHvrY94LxE4Dqodb760dZIdhqZ6iOIswSALIhyL1jWRil1COE/BGAWwHoAD5PKf0FIeQqAHsopTcC+GNCyIsBeADmAbxhiGseOE5fvB0AcLjGVXdpDLDKMEgA17Vh8jRGrXoQ82QSmwutlXhFnuvtVBnplIMqjlpxe6Jt4K+5DJfqKJcrWCDF1O0Mr4p6SnWs6AHj1RbQXGF57YXxTahPnwIssAEYFi+mSs5XTUOZNuCb7H0qW88C7mf3H39cayBZ4GkXPBl7pr+H3902hHazHCG5N5ahNRcQUILxyRn4RhlF2l2BVqzoZ4kXXRUnosrPUekIKSM5VWlsC0B0YD+3ZYbRbiANY8ezgyKlLLgrUDuydmsYNgjhhWLrn9y78twppTdTSs+glJ5KKf0ov+9DnNhBKX0/pfQcSun5lNLfppTeP8xFDxI08MPr7sF7AQClylSYXtioR2RYbBzCgp6u6Erj7H5hF4zTFfiFVnI3MsidOMtYIRUQTYNN0g8ClleFrbemI07wqk7aWIBTZbZCeXIW1k4WkDy043dQKgvV28FzpxQVWkfAT7FPOOtiLNEybibPhKZl2y2EEPzajhkUzQGPd5NQHmPfJ7WrQGMByyhjslwANSsowWZpjJ1QPxpdP3wfoBmsbN6qsCB2fQ7QzMh7HUXoBlPR1GeWzWobdvWKia3A3EPAx08CHr8jur96ZO3OHtYCVjnbc7dXgAM/TX9sxDAiDR3yw8rCUYgwUGnuFwCAsenjofGWoU5tBZhiyrjiHMGylZ7mV+HWiN9YhO86qKAJmiB3t01Wh+4so04qmAHrW562XclfQb24o+X+YqnEulg2FuBV2aepTG/GjrMuwH/M/gKXnr4VTXGQ6kDudmMZBRKA8H4rJxy3Gbe88Ds47cQT2j5vLVAul1ClRWjNRWj2EpYxhildi3rZuPXOvm8tQe6FCUmtVYG6yVR7n3GDoWNyG7C8t6UWYKjY8mR2aS8D9/wDs4wAptxnz1y7dQwbyf49Mm54G3DfN4H3PDpaZ3cpOObbDyzNRS1xT6g/CI9qGJ+YglZi5FYTxUUApvw5NIvp5D42yWwZ2lhCdYkRCEn88129nBn4M92V0HJx9XjfcoFSUIdrtZIXIQQrZAyavQTKPeOJqVkQQnDxWdth6hrGxsZRpUWgfhT2w7chWE4JjAGoLbEzDyKltb3gonNxxtb8PeiCoWMOkzAah2E4S6jx74twcnfqXXSGlJX70uNRgzCrzHbo6qEot3yUIdrlriWpnvZsYJpn7QhLC2Df2UaxZQA+dSrjt3TfN9nl/KPx+w/cA9z8nt4HtQ8Bxzy51xYikpvFPJbJGDRdg8UzUBrLjAxcu4EZLCPI2PFN00KVlkCaS6gusufoiXQ63yjByuipXvBWwkZZXsZBoIIaghRyB4AVbQKWPQ/SXECVllAsxk/VCSFYJFMYW34YhS9dhkc/94bU16nzgLDo0T5qWNKnUWgeRcFdQl1n34VIQ+0mzTOm3IEoU6Y0zYh/ef/aquFeccZz2eWO31y79xybBd5xN5t0NPcwu89eAZqLwNS6LG1JR3kmir1kQR6OAgBfeQ3wk08DRx8Y3rpWiWOe3JtL8Xz7KmGEUeJpfzYPUC4cZv9MbXJr5mutkAp0exH1JeZ7m2NxcvfMcVRoupdX9KtwTB70NMooJMg98Fxu9aST7rJ1PCbsgzAa8+EM2Jb1GdM4u/6fAIBTl36cuk1YHFSeSn08b9TMGZSdOUy6R7BsMp9XL/Iagy4KtMJUR0HqIsd9cjuwtI+Te/4WVEec8zLgbT8Cdg0vOykTM6cAi4+zzomikGqqNT143aK8uVUEJFFL1OmILCL5jCZnHFPkfs+n34z7P/G8eBCVk/dBMNthkXvqIgPFXWFksHjwcQBSBWkKlozNKDUPoz7HDgSVxLZBaQYTtBp7f4FyUIPHyT2wxlBGnNzF2UBWn/FmZRs2+4dQaR7Agpke3KpbCWslpUOkzStlrcpoknujMIsZ7zCmgoWwsZjOLTQRzG6L+jybKjTJ01nly8Y8+xvPPoCPFI4/B9By2IU3n86CuQuPRR0rpzcSuW+Kd9sUkPfbJPlrPA1UTM0aARw75E4pzj/wdZxV/Q/cf9s/hXf7PP3xQOl0AEC9zAh5jLeeFY/X5tgReWxz9ulntXA8JtzDsOfYgWBmW7xLHylvhkECrHBlL6NCawi4Kg8KUxijdQReNA1m6SgroDIz+rYEEydiHHVsdx5GrZSuPOuVeH5+Y7n1B+xyG6o0NZrZD9XxU1BCEzoCBDwHX+fTm5yVDmoL4L3apyNS3/oUdjkp/V9nzxjkkjceNvF6h7mHgEMswwybN9B3VtkM2Eut05jkNNra4cSTuFBKKvocccyQe3U5+scEP742eqCxAJuaaIzvAAD4m88CAIxPTMOjGij/h7qLjFynUqpTBZzKCdjsHwGWnkCDWpjZFCdi0Ztl+eiB2P2eY7Ouizy7hoxthkYoluejgpEqD/wWJ9Nzza3tFwAAKmjCnsooNtoWLxw+su/Blk3cJfY+07OthVqjAHL8ueF17zh2vTDBvmdvObnDpaCxwKpPn/9x4NI/B856EbtfkDwQZYUopCMk9weBg/ew4O5GqE4VEIkQ9YQIk9W6TOK+Fw0s6WTnrCGOGXJfOMhOHxfpGE5q3htaEnqTBVFPeP678K9Tl+O0ZzMPU9M1LJMKtCazKejSAdjUwMzm7EIeTJ2MEnGwae5OHNaOg6bHv16L+/jVhXiVnyBxwgOwxgR7j+W56BSvucQCv5VN6ap851N+K7xePiV9ss1Zv/VqfGX89fjWmR9jr3/wETz25Xfhsa+8O9qoehgNamFqajTTvGZOfxoWaQU2NXHcaaxr5MQMO8vwVrpQTWKE3vgW4Bl/GlVVbjoNOPVSVu4v93FRaEVpimUUHfgp8KsfAdt35b2iwUJUJ7eQu/T7qkrX63OIlPvokPsxk+e+fITZKo9UzsdT6/+OpaP7MTm7DYazhKo2gZ2nnIGdf/K52HOqZAJGkwdHa/twVNuEbXr28XD8lKcBDwCnug/iJ2O/jaTGL3NLxV5MkPvR/ZgBYHJVLlrmVuejTB5viT1nYlO6H7x5ZhO+seUdKMw/gEt/4/mp28xOT+DV7/okDh7cC/zy/XAO/hI7HmCt96n7ERCzCL1+BPNkqu3nzBPnnnwCXuJcBQ86bjmJHSw3TYxhiZYR1FrtrhY05tP9YUKA135jwKvdwNj5TOBnX2fXz3hevmsZNEQHzpUDLK4hINJoZ58UJ/p6hqLPGaO5Bw8BzQVG7o0TWOHFoV+xIlqWUpeeXbJoHY+xJlfMjQNYMNuodgA7z7s4vG5vb51jOTbDSNtNKMzaPFPoxWmmyisz7DJ2EFjaizotYHpTdg72y6+8Ci/6s6+gaLXv8TE7uxVVWsTk498J7zv4MKu6KzUPY8nIP6c9C5NlEx983e/go2+8DONF9jnHiwbmMQGt0YVqqs+PTlOw9Yynvp4N4950GpsEtZEwdRK7FJlAAkKVH/ekdIumMDlSyv2YIXdvkZH01Jn/f3v3HhxXdR9w/Pvb1dteSdbqLVmWLNvIMn6rDi52bCe2wTwawOAYDMm0nQBJaNppGwIN07zaIdOGhgkNSaBJk4YmOKVDcQvGDnHAQHmZ4sSy/MQPWcKSLFmWZC167ukf5y56W7Jed3f1+8zsaH337u7veuf+9uy555yfHRd8sdZOQkjqaqI9dvD1UAJJeaR32eelddYSGOJCZcj0adPZPechXo+5itKNfzLgcX9OAV3GQ3djZZ/tbY22n9vntMqT0+37dDX3JPe4lkpqvNl4x6FF7fV6qPNkUtx28KNt9Sf+z8bYcYbmhPDsbw/55Pws1szrmTTj8QjNkkJc2zBjk7s7nSpL4fvlFTGKVsOX9sM9eyFm4EqlEc2XA944Oxqot1DiziixF1xD9ZA/Svol2nJ3xcWztJoEZpbYftrOhlMApAUb6EgcfGRIMHkmaTTRXF+N3zTSmTz8cK+Nd36Zqx/ahT9tYAKJj4unVjKIbe7bIuh2ZoumptvRH6lpmQRMfJ+JEiltVVyIH78henWJdiRPrUml3cTQUV1OV9tFMk0DHamR1+fcEJtNSnv1pXdqOQuYj6pSqTGaMatnidxo4vHY1vuF/i33c/ZXX2huRKg7pneLPlAfNrNUp0xyjwnU0eBJIzl5BufxIU2VtLU2kUwrwSHGNcfl2tEYx3b9CI8YEmYuGXMcDfF5JAf6zm7zNFXSaHwkp9oLql6vh2pvHolNdhZgsLOd3K5qAr7xS7qBWZ8A4JDvaiq9M0lsPEzdiQP2/TMjb1jb+YQC0rrqelpTg2lykn8kzEBV7korhvp+o8lCyyyEllq46IzOCtTbLqr0eRDssjN2w8CUSe5J7edocfqS673ZJLZWU/+B7ZqJmTH42PWZi1YDcMWxJwHIKV055jhapxWQ1VXdZwJR4sXT1Mbk9ily0ZhURHqbbTlUHX2POOnCm790wOuNVtkf3cMv5j1K4V3/TH3SHDI/fJ/6Q3sByChZPW7vM1k6U2fjwUDjyaF3anaSe4omdzWMnMVw7nDfcpeNp+2vlVDN3VCLPdSiDy1NEib97lMmuad21hFwFv1qScghpeMsFz6wLePEjMG7W7JyCjjoLWU6ASq8JeTmj73lHMxaiI8AtacqPtqW3n6G5qS+XzDtafPINnUELtTSUP6Sjad0/JJucmI8d9zxxxRmpdGRXkq6OY//6HbOmEyK50TeCn8xWXZ+QuDMe3D6DTjzzsCdQok/JbyvKagwkLsETBBq7K9ZjLF98DOK7CQn6Olfv3DGNhj6b3fZlEjuJtiNP9hA53Tb/dI+PZ+s7jo+PG0vIubOWzbkc7uve4SXEzfQff2j4xJLutMqPlv+CmCHaGaZetrTr+yzn2++rbxz8o3n8J14nhMyk6LiknGJob9pC+zQybz296nIvAHvJdZtD1e+WYtpMkl0v/c0/Ou18OP10NW3TCK1FbYvNZom3KiJMfMq29Vy3DasuFgHHS0wo7BXt4wz4KHxpJP0ne2a3CfP+bpq4qQbcX6Oe3IWES+d5J96hg/IwJ8+9PDCRcv/kLVfeYaFy8beJQMwe8Fy6k0KnmMvAlD53m4AUub1HTo5f/nHOU0OxW99lTntFZyctWXMtUmHsrxsJT/3/wX/EXMj8zf/zYS8x0S7IieVnd0r8FW90rPx5Kt9d6o5AFl9v0SVGtQ0PxSshEM7bKu9yvklmLfMXkT25dpum+4uu4hamiZ3V5yvshdG4v22+2XmMjvpIidYw+nUj01qLLExMRzyr6ek6XWazr6P53fbOWdSmbdsbZ/94uPiOHX1tzlnUng5dhVlt/714C84DkSEu/7sG9z20FMUZEfmutyF/iT+LWEbbyeu5u8999BOHOZ4zzh+mqrtdPlZg8/eVWqAxbfbBH70RTi229bZzXEGVeQssmUO6yrsRdTMUmf+hPRcaHXZlEjuTad/B0D67MUA5BXMZk/qrVSaLPwbvzzp8WRvup8uvMgTayhtfZPf520lIT5uwH5rNtyE7ysVXP3Af5MyPWnS44wkIsLHFi9gS+PneTKwhje659NxeFfPhevQbMq5G90LUkWWxVvBPxf+83Ow/99h4a09JQ0LVkL9EbsdoOAqW/5wRmHYrOk+JZYfCNaU02oSyCvsuVC47s//ha6gIdaFafZz55bw26t+gO/tRziXPJ9V27425L6pSQOTvhrcfevmUNvcxtxMHztfXsHapifhyE7byvrf70Hh6ugqB6cmljcW7tgOz/+VrbO7/hs9j115C+z5Frz1Q7voXGhWa/bCnouwLpsSyT274S2Ox5ey2NtTvFlEiPW6d+Fw3abNsGmza+8fjfzT43l8m135cvORG6hs3EXB9jvtbEOPF677R5cjVBHHXwyf+a+B21ML4MbvQfkzsOGbPduzF9l++sB512usRne3jDG8/8ovKQhW0TJrg9vRqEm0dn4et7Q+SGDZ3bD0TvjcHjuDUKnxsnQb3PWsba2HFK+zf4/tdiemXqI3uRtDxaM3Uvzbe6kmk4XXf8HtiNQkWl+aRT0pPJv5ebj+O9odoyZH7jJIKYC3nxi00tlkitrkfuD5xyltepU96XfA3S+TkhqeZePUxCjJ9rEgN5nv/vooJ84NXrdWqXHn8cCa+6H6XXj1O64m+BEldxG5VkSOiMhxEXlgkMfjRWS78/hbIlI43oFejlNv7qBo37c46C1l9b2PkZer082nGhHh0U8vwRi44bHXeHjnISo+aCYYdLc1paaAJdvgys2w5+/gqVvg6K6+yxhMEjHDfLOIiBc4CmwAqoB3gNuNMRW99vkCsMgYc6+IbAVuNsZ8+lKvW1ZWZvbt2zf6yLs76W5toKWhhpbzNTRXH6K7poK0ujfJ76rkFLmYzzxH0ezIWwRLjZ+qxgAPv3CYneVnCRpITohhYX4Khf5p5M9IIn9GIiuK0shKTnA7VBVNgkF483F47Z9spSZPrF1YLOMKO1yyeJ0teDIKIvKuMWbY8lcjSe4rga8bY65x/v0ggDHm4V777HL2eUNEYoAaIMNcWNwmFgAABJhJREFU4sVHm9wPvPQUha/dj4/WAY+1mESOx8yhYdYm/uCmL5GSrNPMlVV/sZ09h+vYf+YC5dVNVJ4PcCFgCyD/6K7lXLPg0oVYlBqVzjY4/Tqc3GsnRJ07DE1VsOov4RNfHdVLjjS5j2QoZB7Qe43aKqD/tM6P9jHGdIlIE+AHxn15tFh/Ie+mbqQ7IQ2S/Hin+YlNziQ1/wpmFc1jacKlqxCpqSl9ejxbymaypaxngbaWtk6qL3xITkqii5GpqBabAHM+aW8hwW7o7pjwt57Uce4icjdwN0BBQcGoXqNk6SpKlq4az7DUFOVLiKUkWxsDapJ5vOCZ+AbFSC6oVgO916PNd7YNuo/TLZMCDKhWbIx5whhTZowpy8iIzDVMlFIqEowkub8DzBWRIhGJA7YCO/rtswP4rHP/VmDPpfrblVJKTaxhu2WcPvT7gF2AF/iJMeagiHwT2GeM2QH8GPi5iBwHzmO/AJRSSrlkRH3uxpgXgBf6bfvbXvfbgNvGNzSllFKjFbUzVJVSairT5K6UUlFIk7tSSkUhTe5KKRWFhl1+YMLeWOQccNqVNx+bdCZg5q3L9Jgigx5TZJjoY5pljBl2opBryT1Sici+kazrEEn0mCKDHlNkCJdj0m4ZpZSKQprclVIqCmlyv3xPuB3ABNBjigx6TJEhLI5J+9yVUioKactdKaWikCb3ERKR20TkoIgERaSs32MPOvVjj4jINW7FOBYi8nURqRaR/c7tOrdjGq3hav5GIhE5JSIHnM9mDPUp3SMiPxGROhEp77UtTUR+LSLHnL8z3Izxcg1xTGFxLmlyH7ly4BZgb++NIlKKXQVzAXAt8LhTdzYSfdcYs8S5vTD87uHH+b//PrAJKAVudz6jaLDO+WxcH2Y3Sj/FniO9PQD8xhgzF/iN8+9I8lMGHhOEwbmkyX2EjDGHjDFHBnnoU8DTxph2Y8xJ4DiwYnKjU72sAI4bY04YYzqAp7GfkXKZMWYvdknw3j4F/My5/zPgpkkNaoyGOKawoMl97AarMZvnUixjdZ+I/N75qRlRP497iabPozcD7BaRd51yldEiyxhz1rlfA2S5Gcw4cv1c0uTei4i8JCLlg9yiouU3zPH9ACgGlgBngUdcDVb1t8oYswzb3fRFEfm42wGNN6d6WzQM3wuLc2lSC2SHO2PM+lE8bSQ1ZsPCSI9PRJ4E/meCw5koEfN5XA5jTLXzt05EnsV2P+299LMiQq2I5BhjzopIDlDndkBjZYypDd1381zSlvvY7QC2iki8iBQBc4G3XY7psjknVsjN2AvIkWgkNX8jiohMExFf6D6wkcj9fPrrXX/5s8BzLsYyLsLlXNKW+wiJyM3AY0AG8LyI7DfGXOPUk/0VUAF0AV80xnS7Geso/YOILMH+LD4F3ONuOKMzVM1fl8MaqyzgWREBe87+whjzorshXT4R+SWwFkgXkSrga8C3gV+JyJ9iV4nd4l6El2+IY1obDueSzlBVSqkopN0ySikVhTS5K6VUFNLkrpRSUUiTu1JKRSFN7kopFYU0uSulVBTS5K6UUlFIk7tSSkWh/wcchiU9npUD1wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dos.plot(plt)\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "energies = np.array([-5,0,3])\n", "weights = np.array([1,1,0.5])" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "dos.append(energies,weights=weights,eta=0.1)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsvXncJFV9Lv6c2rvfbXYYGHBAEcJumAu4YDBGIOrVGImaK0YTkZBoFn9K1FyNgY/GuPz43RgTFI3XaLhJxAWJslwSYtQkQmCcYRdGRJhhBmZ937e3Ws/vj3NO1anqqu7q7uq3l7eez2c+3W93ddWpnq7nPPWc70IopShRokSJEtMFZdQDKFGiRIkSxaMk9xIlSpSYQpTkXqJEiRJTiJLcS5QoUWIKUZJ7iRIlSkwhSnIvUaJEiSlESe4lSpQoMYUoyb1EiRIlphAluZcoUaLEFEIb1YE3bNhAt27dOqrDlyhRosRE4t577z1AKd3YbbuRkfvWrVtxzz33jOrwJUqUKDGRIIT8LM92pS1TokSJElOIktxLlChRYgpRknuJEiVKTCFG5rmnwXVd7N69G61Wa9RDGRtYloUtW7ZA1/VRD6VEiRIThLEi9927d2Nubg5bt24FIWTUwxk5KKU4ePAgdu/ejRNOOGHUwylRosQEYaxsmVarhfXr15fEzkEIwfr168s7mRIlSvSMsSJ3ACWxJ1B+HyVKlOgHY0fuJVYWe2p78Njhx0Y9jELxyKFHsHP/zlEPo0SJkWKsPPcSK49Lvn4JAOD+t94/4pEUh1/7p18DMF3nVKJEryiVew5cfvnleOihhwrZ11VXXYXTTjsNV111VSH7K1GiRIk0lMo9B77whS8Utq/rr78ehw4dgqqqhe2zRIkSJZIYW3K/+p8exENPLxW6z1OPmceH//tpHbep1+t4wxvegN27d8P3fXzoQx/Cddddh0996lPYtm1b2/a+7+Ptb3877rnnHhBC8Fu/9Vt497vfjV27duHKK6/E/v37oaoqbrzxRrz73e9GrVbDOeecgw984AN44xvfWOj5lShRooTA2JL7qHDbbbfhmGOOwXe+8x0AwOLiIq677rrM7Xfs2IE9e/bggQceAAAcOXIEAPDmN78Z73//+/G6170OrVYLQRDg5ptvxuzsLHbs2DH8EylRosSqxtiSezeFPSycccYZeM973oP3ve99ePWrX40LLrig4/YnnngiHn/8cfze7/0eXvWqV+Giiy7C8vIy9uzZg9e97nUAWJZpiRIlSqwkygXVBJ7//Odj+/btOOOMM/DBD34Q11xzTcft165di507d+LCCy/EZz/7WVx++eUrNNISJUqUyEZJ7gk8/fTTqFaruOyyy3DVVVdh+/btHbc/cOAAgiDA61//enzkIx/B9u3bMTc3hy1btuCmm24CANi2jUajsRLD7xsBDUY9hBIlAKcONA6NehRTga7kTgixCCF3E0J2EkIeJIRcnbKNSQj5R0LILkLIXYSQrcMY7Erg/vvvx7nnnouzzz4bV199NT74wQ923H7Pnj248MILcfbZZ+Oyyy7Dxz72MQDAV77yFXz605/GmWeeiRe96EXYt2/fSgy/b/jUH/UQSpQAvvI64BNlHaUikMdztwH8IqW0RgjRAfyAEHIrpfSH0jZvB3CYUvo8QsibAHwcwESGglx88cW4+OKLY69997vfzdz+rLPOSlX3J510Eu68886212u12sBjHAb8wIeuTFflyYAGUEh5c5obd38eWHwKeEVnK3KoeOou9kgpUJbeGAhdf/mUQTCSzv/RxGavBfC3/PnXALyclEVRJgrTaMuUdyM94pb3Av/+F6MeBYNnj3oEE49c0TKEEBXAvQCeB+CvKKV3JTY5FsBTAEAp9QghiwDWAzhQ4FhHjvPOOw+2Hf/RfeUrX8EZZ5wxohEVB496ox5C4ZjGu5EVgWcDmjnaMbgNQC+jzAZBLnKnlPoAziaErAHwTULI6ZTSB3o9GCHkCgBXAMDxxx/f68dHjrvuSs5p0wM/mD6VO413IyuCxiFgfvNox+DUgOq60Y5hwtGTIUkpPQLgXwFcknhrD4DjAIAQogFYAHAw5fPXU0q3UUq3bdy4sb8RlxgKptHCmMa7kRWBOwaRXc4YjGHCkSdaZiNX7CCEVAC8AsAjic1uBvBW/vxSAHdSSpO+fIkxxjQq92k8pxWB74x6BOMxhglHHltmM4C/5b67AuCrlNJvE0KuAXAPpfRmAH8D4CuEkF0ADgF409BGXGIomEblPo3nNDQEkoU1DouZvjvqEUw8upI7pfQ+AC9Ief1PpOctAL9W7NBKrCSmReXKN4xeUNoyuSFbMaNSzfIEE5TkPijKIOAcKLKeexa+//3v47TTTsPZZ5+NZrM51GOlYVpUrryIWjqDPcCVfnPeiHr2ypNKqdwHxtgWDhsnFFnPPQs33HADPvCBD+Cyyy4b+rHSMC2RJfJ5TMuEtSKQlbs3IuXuS3ZQ6bkPjPEl91vfD+wruE3a0WcAv/znHTfptZ47AMzOzuIP/uAP8O1vfxuVSgXf+ta3cNRRR+HGG2/E1VdfDVVVsbCwgO9973vwfR/ve9/7cNttt0FRFLzjHe9ApVLBV7/6Vdx+++249dZbccMNNxR73hmQle3UkDui85iYcxqHbExZufsj8tzlSaW01AZGacskIOq579y5Ew888AAuuSQZ9dmOer2O888/Hzt37sRLX/pSfP7znwcAXHPNNbj99tuxc+dO3HzzzQBYJ6YnnngCO3bswH333Yc3v/nNuPzyy/Ga17wGn/zkJ1eM2IHpVLnyOU0EuT/4TeDqNUDt2dGOw61Hz0e1oBpT7qUtMyjGV7l3UdjDQq/13AHAMAy8+tWvBgCcc845uOOOOwAAL37xi/G2t70Nb3jDG/Crv/qrAIB//ud/xpVXXglNY1/9unWjS9SYSJXbBRNH7vd+iT3u3Qmc9IrRjSOm3EdkiciTSrmgOjDGl9xHBFHP/ZZbbsEHP/hBvPzlL+/6GV3XIUrpqKoKz2O3lJ/97Gdx11134Tvf+Q7OOecc3HvvvUMde6+IESEmgAhzYOLuRqw17LE+4kod5YLq1KG0ZRLotZ57J/zkJz/Beeedh2uuuQYbN27EU089hVe84hX43Oc+F04Ahw6NrnZ1jNyD6SP3iVDu5ix7rD0z2nGMw4KqV9oyRaJU7gncf//9uOqqq6AoCnRdx3XXXYf3vve9fe3rqquuwmOPPQZKKV7+8pfjrLPOwumnn45HH30UZ555JnRdxzve8Q68613vKvgs8mHiVG4OTBy5gy+k1vePdhjOGMS5y8ctbZmBUZJ7Ar3WcwfiNdovvfRSXHrppQCAb3zjG23bapqGa6+9Ftdee23s9S996Uv9DXgATB4Rdoc8SU3EOQnFPPIFVYncR0WspXIvFKUts4oxjeQuh3dOxN2IUMz1UZO75LmPKgyxjJYpFKVy7wHTVs99Gsl94s7J4Xd9B3aNdhyxaJkRkbtX2jJFoiT3HjBt9dxLz30MIKyIpd3A8j5g7ujRjMOtA4oO0GB0xFpmqBaK0pZZxZjGOixySOdETFgykY5yUdVtAnoVUPXR2TKych/V3cMUoST3VYypVO7BhE1YgQeYC+x5o62/zcrBbQB6BVC00RFrGS1TKEpyX8WgmO7aMhMxYQU+MLOePW8eGd04HIncR0Ws8nHLBdWBUZJ7Aq1WC+eeey7OOussnHbaafjwhz886iENDRMXNpgDE3dOgQeYc+z5qDJDAcBeZuNQ9dERq9xToCwcNjBKck/ANE3ceeed2LlzJ3bs2IHbbrsNP/zhD0c9rKFg4hYfc2DiKl36LmDwLNWRkvsSYC1w5T4iYg2PS4Y3wRx5EvjTBeDR/zuc/Y8RSnJPgBCC2Vl2sbmuC9d1w7ox04ap9NwnbcIKfIncR9jerrUEmPMsYmbU5K5XhmcN7b2PPd7zxeHsf4wwtqGQH7/743jkULIP92A4Zd0peN+57+u6ne/7OOecc7Br1y68853vxHnnnVfoOMYFZeGwMUDgRfVl3JXvwBXCXgaseUDVRmjLcHLXrOEt6iqc8iZh4h8QpXJPgaqq2LFjB3bv3o27774bDzzwwKiHNBTELIyycNhoELiAMcOej1K524urQ7mLCZ9MP/WNrXLPo7CHjTVr1uBlL3sZbrvtNpx++umjHk7hkJXtRKjcHJg8cvcA1WSk6o1IuVMaKfeReu78N6hZw7t7EPtV1OHsf4ww/dNXj9i/fz+OHGEhac1mE3fccQdOOeWUEY9qOJg4IsyBiTunwGcRKqOMUnFqzKYwx8CWIQqgGsNT7mLRulTuACHkOABfBnAUAArgekrpXyS2uRDAtwD8lL/0DUrpNcUOdWWwd+9evPWtb4Xv+wiCAG94wxvCLkvThokjwhyYuLsR32UqcpR2iM3r25izfBwjJHdF4xPMkL4LUf1yFSj3PLaMB+A9lNLthJA5APcSQu6glD6U2O77lNKJZ8EzzzwTP/rRj0Y9jBXBNLbZm7jErBihjUq58/6pxuzoQyEVbbgTjKjCuQqUe9czpJTupZRu58+XATwM4NhhD6zE8CEvok6Eys0BP5jAJCZFH61iFs2xjRluD43Qc1c0XgJhSN+FiEgqyT0OQshWAC8AkFYe8YWEkJ2EkFsJIacVMLYSQ4as3CeiDksOTJRyDwIAlCv3EZKqUO56dcTlBzxmlwyzeJmwZcb9t1EAckfLEEJmAXwdwB9SSpcSb28H8BxKaY0Q8koANwE4KWUfVwC4AgCOP/741ONQSqc2aagfDJN0J66xRQ5M1DpCIEVujJJUx8WW8d1IuQ8r5l/sN5iO33sn5FLuhBAdjNhvoJS29Y6jlC5RSmv8+S0AdELIhpTtrqeUbqOUbtu4cWPbcSzLwsGDB6dGRQ4KSikOHjwIy7KGsv+Jq8OSAxO1oCpIdOTRMuNiy3jRXcywJjqh3FdB7Zo80TIEwN8AeJhSem3GNkcDeIZSSgkh54JNGj3XL92yZQt2796N/ftH3Cx4jGBZFrZs2TKUfU+Uys2JiaotIwgmXEQcsS1jjNqW8aPvYmjRMly5j/tvowDksWVeDOAtAO4nhOzgr/0xgOMBgFL6WQCXAvgdQogHoAngTbQP+a3rOk444YReP1aiT0wjuU+Ucvclch+XaJlRNuuQI4dK5T4wupI7pfQHADqa4JTSzwD4TFGDKrEymLg6LDkga4qxt/dC5a6OR7SMWFAdZRJTqNyHHC2zCsh9+uOBSmQiRoQYcyLMiYlS7iG566NdyBSx33pl9LVlhu65l+ReYhUgRoRTEj0wWZ67iJYZcSik1wS0CkDIiMsPrITnLmyZMf9tFICS3FcxJiomPCcmKutWTKgi/G9UtoxnA5oZjWWkyl0tPfeCUJL7KkYsm3NK6rlPVHhnGAqpjTYU0m0ySwYYD1tmmGMobZkSqwGxZtJTaMtMjuc+4lBIWbmPuirksC0qERlUknuJacZE+dM5IRP62EfL+LLnPkJS9VrMcxdjGYT4dv4D+9cP5Noyw7BlRN16IGraMcUY22YdJYaPibIwcmKylLvw3EdcOMxrSZ47HwelbIG1V3zzt9njWW/q/bOBx8YxLIvKa0Xf8ZTcqXZCqdxXMaZduY/9Oclx7iONlmlFnruq87GNgPxkz536bIIpEgd/Ej/WlKMk91WMiYoJz4nJUu6SLTPKaBm3FY+WAUZDfnKGKlC8er/pd9gjUUpyLzHdmMbyA/J5jL3nLi+ojjJaJum5A4NPNP3Ekctx7kWMIQm9yh5PeVVpy5SYbkwjuU/U3YhcFXKkce6Sche2TD8TjUyYTq2Pz0v13PsdQyfQADjhFwBzoST3EtONgcm9tQQ0DhU4osExUesI4YKqqC0zIsKRPfdBbBm5BrtoRN0LZM+93zF0gr0EWAvs+y5tmRLTDDlDtR+Ve+gzLwc+MV5VPMd2QXVpL3DnR+ME3hYK6YxmbMkMVaA/8pPH3xe5u8P13J0Gb0hSknuJKUcsJryPwmHrao8WOZxCIM6DgIyXLXPLe4HvfQJ46u7oNblwmGqMMENV8twHsmUkwvTsPj4/ZM/dd9jEoWirIs69JPdVDGFhaIo2NRmq4jw0RRsv5d48zB5ldZvMUAUdjTWTjHOXx9YL5M/00yZv2J574EoVOKfj994JJbmvYghlqyv6QETou32otCFBKPexI3fREkEmvVic+5CsiG4IAsC3pTj3AWyZgZW7Fy0u9zuGTvDF/ktbpsSUQ5CfpmgDFQ5r1pP90keH2DmNFblziKqEQCJaZkhWRDd4fLLReJ9eZYBJJkbuAyyoDlW5awApyb3ElEOQn67oCAaob90aI3KX70bGynMXqfxp5D5MQuuGFv+/s+b5WAaYZGSrYyw9d1cKO51+ci9ry6xihCqXaAMRoeOMKMojBfI6wlgpd0HujkTufiJDFVh50rE5uZuc3AcpP1CIcpc99wK/C0ojz50oLOa93/o5E4JSua9iFEWE4+S5F7WOMDTIpCc36xi5cl/gY1H7H0fhtkyBokH8FmKe/hjd2Q0BJbmvYggiHJjcvVK5d4VIrvKliTBmyxj8/RX+Lu1F9mgWYcsURe5G/2PIQuwuSYmON8UobZlVjKIiSzxn/JS7StTx8twFkXgZ5D6srMxuSHruRZUf6JXcg4Cp69hEVyC5i4lC1SMVP+Xk3lW5E0KOI4T8KyHkIULIg4SQP0jZhhBCPk0I2UUIuY8Q8vPDGW6JIiFiwvtafJTT/EeVfJOCgAZQiAKFKONVOEyEQMbIPZGhCqy8LbO8jz3OHhWNBejTcx9gQZVKpRiGYcuEyn2IoZZjhjy2jAfgPZTSUwGcD+CdhJBTE9v8MoCT+L8rAFxX6ChLDAUi/LGfUEjfi0honDx3ChqS+4oq99qzwBP/nv1+KrnLDbJHFAq5tIdlp1bWRmPpdxyD2DLDtqjkfrXiHMfJthsCupI7pXQvpXQ7f74M4GEAxyY2ey2AL1OGHwJYQwjZXPhoSxQKYcWoitpzKKTrRBdvMEaeu099KFCgEnVlPfcvXgx86ZXZylvEkyc9d6IwD3gYESICD94E/J83pTe/OPwEMH9MFDVSVPkBtwhyH4bnrkeLxqVyj0AI2QrgBQDuSrx1LICnpL93o30CKDFmEBZGP0TouhGhB9742DKUUqiKCoUoK0vuhx5nj42D6e8LspOVu+9GKrKoOuppuOl3gUdvbR8bpcCTPwS2bIteK6r8wEDKfQi2jOy5l7ZMHISQWQBfB/CHlNK+slYIIVcQQu4hhNyzf//+fnZRokBQSqGAWRi92jLyIuq4KXcCsvLkLiAaMCfhpZB74EVkOsxQSLfOHluL8deXngYaB4Djzo1eG0TVDuK5x8JCh2DLiDui0nOPgxCigxH7DZTSb6RssgfAcdLfW/hrMVBKr6eUbqOUbtu4cWM/4y1RIHzqgxBGhL0WDvNjyn2MPHdKoRJ15T13gVaG7sny3EPlPgS1mkRy4qlzgTUnOaiFVYXsV7mrQ46WGWHC2AojT7QMAfA3AB6mlF6bsdnNAH6DR82cD2CRUrq3wHGWGAJkIuxV5XqerNzH5yIJaABCyMp77gL2YvtrvheRi5+IlhFKOYztHuJ3aScmnvoB9ljdEL1WmC3Tq3Ifsi2TGi0zRqGyQ0CeOPcXA3gLgPsJITv4a38M4HgAoJR+FsAtAF4JYBeABoDfLH6oJYqGUO4qUfuwZaILj46RcvepD4UoIISMjy3jyR2KkrYMvwRXIhQy+f8klPuMTO5FVYXsseRvrHHJMKJlJM89nGjHZ61oGOhK7pTSHyCsV5q5DQXwzqIGVWJlENAAKlEZEfYYLRNIFx4do4uEUhouEnujuO1Os2XkyJEscl+JUMgkuYs+pyI7FRhskhFK2JiNjvUv1wAP3Qz83j35PiuHhRYaLSPZPqvElikzVFcxZAujV386ptw7KKzg4VtAaQD11Ff3Pc5eIJT7invuohhVWpMKL6O3qO9FZDrMUMjweAlyF/9vmhG9VkT5Ab0akfv3/1/22K1Il+y5Kwoj4GEod0UHFDd+zClFSe6rGEK59+O5y/VkOil35R9/nT350xQveggQSUwqUVc2Q1XRGXmKyBQZsnKPdWJyIzIdZiikQDKqSYxFlcm9AFvGmGm/S/CdqNtTp8+G0UPGcDx3VQf8ASpfThDKwmGrGGE2J3oPhfRj8drjo4D8IPLcV1S5C1JMU+6ihrs+k1DubqTYV6IqZJIsQ8KTyH2QOwiZ3JN3CXId+06fVaQ7mUKjZeRQyDKJqcSUQ8SE95OhKicu0TEidwoWAbTi0TKCKJw05c4Jv7Imrp7lOPeV8NzbyN3hGbJq9JqiAiAD2jKV9rsEpxu5S547METlXoZCllgFGCQmPFZPJg8RrJBFEtBgNElM4jtIU+7NQ+xxbnOKck967iu4oOo70aQio19iDRdUZ9rj3Ls1zJY990HGkLn/snBYiVUEOYlJlP/NCyqXHMhBSCtVXGwkC6qBHxWhSiMxkfY/f0y25z60vqHSd9C2oOrGLRkBVe/PlhHft55my6Tc0cTGOWRbRvbcV0mce0nuqxjygmrPGarybXcGicp2jWP3GPfcJ0Qo5Ioqd5mE0khs178wEl84LiVaZsi2jDyZJMnSd6Ljy5BjwXtB6LlXUxZvu+yvjdyLVu5y7fzScy8x5RgkmzMWIZNBBK5Uf8ZpdfFcC4JQ7ivqucvnn6bcf/pvwKmvYQ0xAi9SjIEbV6pA8YvTMkGm2TJpyl3R+7Rl5FDIhC3TF7kPW7mX5F5iShELhewxWoZKyoxkEJIjlQUWyr3x7Q/A3vn1Pkabc1xcua9ohqrfgdyDAGgeBtafFIUCCpKVo2UGWcjsODbp/yYtWibVljH6jJaRkpiSx+p2Xm0Lqj1OMI//W9R4JHX/pedeYhUhZmH0HC0jkTtNv0g8OyJ3t9UE9V1U7/lrmN/8rf4GnAOycl85z106/2S0jFCwugWonNyFHy177kDxPjOQQ7mn2TJ9JhDFomUyEqa6fbafBdXAB778GuCLl2RvIyarUrmXWA2IVYXskQjjtkwGuUuLqK7TQHP5cF/j7AWDFEPrG52UuyB3rZKi3KUMVYARfdGEE/Pccyp3ZQDPXdHYeVI/kbzV5bwGsWWa/Hd1+Kcd9i83yBaee7mgWmJKIceE9xotI4iiRXWQDEJyJOXu2U3UFzMaWRSIACOoCimTZnJBVZC9brWTe6pyL7jkr/x/07ag6mYo9z79bpncAaB1RHqvV8+9h+9CVLfMs/9SuZdYDfADptwJes/mFJ57C2Ymucvhj67dQGNJIvc0G6hbLHQOiDZ7K6rcwyiRuc7KXe3guQPDt2WSBJu1oKpq/RcOU7ToPOXmIF0XVAdIYmrkIHe56mRJ7iWmHQHYgmo/GaqUXxg2MTM9d7mhh++00Fo+FP7dWo6reLr7XuCjR8N/9J97GkfbuCiFoqxwnLsgDmuhndzD0gOSLRN67l5cufdrh3QcWydbJovcjQFsGTUqRNaUlHvuaBnhufcw0dVzdHWLlR8oyb3ElCOWzdljtIyIY7ZJB+UuLar5ThNOPfLclw8/G9v2obtuBwD8+Ps39jaO5DEl5b5ihcMEEVrzbEFVPq7wnfU0z92Ne+6qNoRQyD5sGaXPO4jQlrHY3wPZMj0o9zy2jO9GzchLci8x7RANshX0YWEELgJK4BIjU7m7Uiik77Tg1aOLvZ4g9yXOd8vNwTxnOQKoUOX+7CPA9z6ZXkZBEKi1AIDGI0VEuV8th+c+EuWekcQ0CLmrgyj3fmwZ6S4wa0KPVeAsk5hKTDnCyBKl9wxV4rtwoSEgKpSMiySQPPfAbcFvRBd7cyl+K014rW/S68JuAkNLYrrpSuDOjwCLT7W/JwhZNL2QKyDKyj303PlrcoYqMBzPXZ4scse5D5DERNSMBdVeo2V6sWUk5Z7VFSzWr7ZU7iXGHIt7HsUzf3Y6Dj2+va/Px9rs9aHcPagIiAYlQ7nLjbOp2wSVFtic5fittEIYqQ9qpQwtFPLg4+xRVqMCggitBfYox7oLopeV+7feyRYb5QxVoH87pBPChhyV/OUH+g6FTCyoxpT7EOPcZc89q7RwWpG2ktxLjCse++7f4SjnKTx2xxf7+rwgQgLSu+ceeHAJV+4Z9occLRN4Noi9hIAyhe7VDsW2JfxCJgMSsiipUDi5CyJopTXAFsp9jj3GbJkUz315L/DYHexzMrlrBddTkcdmzPSwoNqvLePzBdV+omWkUEWgt3BM2ZbJirhKbYxSxrmXGFOI2i0B6e+/MVTuClPuvahm4jvwoMHvqNzl+O8WFGcZ+7AeHlVAG3Fyj5psD6bcw3WEoj13QVh2So9UQUyc3P9j393YU9vDXgvj3CvRQiPAUuWpH1fOmjVEcq/2aMtMaJx7R+XOv2txvZTKvcTYgl8I/apdubaM+DsvSODBgwZKtEzlLpM79VpQnWU01RksYQakmU7uZEA1Jci9cM9dkE5a04lQuc8CAH773o/hkq/zVPgwzt1iBC+wxMlfXlDVzPaCW4NCEKQx207YQYckpkEyVFMXVPvJUO0hzr2yjh8n4zNy2Ckh7DgluZcYV5ABFYhcFRJAT9aMErjwicptmfTjU7mGu2ej6h5GXZ3HsjIPzY5714STWta+8kKu514ouYvmzsk65UBEhMZc+3tynPvcZgB8P4u72WNSuRdO7nxseppyz6oKOWASUxgKKdsy3Tx3PqkTyXOnQXfrJAiYLTO3ufNxRAy+QEnuJcYa/MeZFYrY9eMIwphwoEflTj340PiCakY9d0m5V1r78VznYSxZx6KmzMNwEguTvlDugy0ohn1heQOSwmPd06IxhCrlyj0GV1LuhAAfOgAc8/PAkSfZ67Ka18zsaI9+Eciee1qce9G2jNK/LSPi0MUYgO6TgrPMJoHZjXz7jOMkY/oVrfTcCSFfJIQ8Swh5IOP9Cwkhi4SQHfzfnxQ/zBKp4ETQL3/J2ZwAegqHVHwHHtFBiQY1S7lzwnapipMP/jMUUOw99hK09AVYXnxhUhG2TBHKXZqwCs9STSMb8Zo5175i4DVZpIpQ/qoGzB4FHPkZ+1uYhTlMAAAgAElEQVQmd3XYtkzeOPdBasvo7baMMZvPlpEXl8U+upG7mDxFpFLWuNuygdVSuQP4EoAOtTQBAN+nlJ7N/10z+LBK5EKodvv7kQoiFLZML8XDFOrCIzoCpQO5c+VeRwUm2PMXXvgqOPoaVP34wiTh55IVM58XlFKoihqdU2HKnZNzmrIOomiZttG7LVY0TMbsxqiSoTZk5Z61oEppZ1tmoKqQUoaqaubz8LPIPdnRKQlhe4kcg6zJIJkNXNoyAKX0ewAOdduuxMpDECKh/VkZcigk0JvKVQMHvqIDigoFGZ/jF1qdVAEwBX/UhvWgxgwsGleoSsDJfUDlLkoqiKSolVHuURKTKxS6gFDuMmY2Rc9jtgz33B+7A3hssBo7bWNL2jLJ0EMZgzTIVrSotoxT48lbOfYX+JHfLo+r26Tg5VXuyWzgktzz4oWEkJ2EkFsJIadlbUQIuYIQcg8h5J79+3MU+ynREcQbXLmLUEgAPRUPYwuqOgKiQc0iUN5Eu6Uwcl/CDHRNBdWrqCTIXR3wLkRAjpYRfxcCsZ9Uz10Q6Cy8BLfDbaYo96Oi5/J7QrnfcClww+sHHjIbW4YtI54raeTO68r3etcTuMzuEElMACf3HA2325R1Ts9dKPeQ3DstqJbKvVdsB/AcSulZAP4SwE1ZG1JKr6eUbqOUbtu4cWMBh17dIAFP/BkgWiYWCtlDtIxKXfiKDqroUDsod5tqcBV2sdcUHk2iV2EQL9bwQ+HnkmXx5IUc5y7+LgSCNDpFy5hz8JBgd7fFIlVkzEq/ffm9ZLRMWjZsr0hGywjCFueTtaAK9E5+yVBIICL3rrZMsra98Ny7jEF47l1tmUSpB0UtF1S7gVK6RCmt8ee3ANAJIRsGHlmJrlBCn7o/WybM5kTvRKhRFwHRQRUVWia52zzRiV2odUHuBiM0u1FD8N1PwH/gm1A5uRcZCin+LgSCNNI8YGlBNd2WSSh32ZbREspdhoioGQSixIGwSgShyQ2jk8irmtuOxQlUUSJy1ir5Gm63ee5avjGIJLFctszqUu5a9006gxByNIBnKKWUEHIu2IQx/JY7JUK1Owi5K0SBovQeLaNSD4FiAB2UOwlcuESDp7AL3daYuiLGDACgdfhpWN/9KABA057PzmVAcperQgIrpNyFuky1ZVpxXx1I2DIJz13GkSeBzWf2NdxobE48gsV34j1SOyn3XiNm5OJcqsmOIfqpdrVlEi0H80bLiKqbFlfuWddCaijkKid3QsjfA7gQwAZCyG4AHwagAwCl9LMALgXwO4QQD0ATwJvoihXSXt1Qw3osg9ky/fjTGnURKCwUUss4PvEduNDhc1vGMdgFqHByP/TUw1gT7s8J9zsIPOrFzqkQ5R4EERFkKXdFAxQFrppQ324DmElYkLItI+wEYDjK3fcYUcbIsiop94xmHUAf5C4lCmkG4ICROw36UO5G9HontCn3HBmqQP+JWhOEruROKf31Lu9/BsBnChtRidwIrYwBFlQ1RQujZXrx3DXqgqoGoGqZtoyoPxNwwvMMRuWqxcj94FM/xol8W50KW2bw8gOqEq0jFKIzYguRGZ47Jw4vSdBus125y4Q+K1s0ic82CwhSE7HsSTXe0ZYRhbX6IXf+WXEXolfYGLqSe0o0ixh/JwhyDz33TnHucoZqGedeYoyh0cF8ai/wWGRJH9EyOlxuy2iZtowSsFh4cVEFFiN3zWKZnO7+n4TbGpSR5qALqn7gxxaJC1HuMqFnhUJypelqCSXs1KNqkQKEADqb4NrKD8hoHsbACMk9YXN0tGVyWiJJyOQu9qFXeZx7tyQmP125d1PXeZV7my0zhPLKY4aBPfcSo4M6YIRJMlqmFyLUwBbPqKJDJRQ08EFkZQS2JuASPRwnqaxln+XK3apFjS8WKEtqUtvTgHqCRz1oilZsKKRMAmmEIIXxycqdUgri1FiMeRK/+x+AXYu/1qbcC4qW6ZXcxZ1Grw3LZYIW52LMsEqa9nKOccrkntP3T3rumbZM4s5AM/uL5Z8glOQ+wdAHVO7ClumncJjOF1QJv5g9z4VuJMmdKXfTYySmzLIgKqPClOxa5+lwWxPsItYGVNp+4Be/oOp1U+5OQrmzbTzqQXfqLMY8ibVb21+TlXtlXXb52l7gO2xxM7Q5ctgyYjJKq4DZCTHPnZO7Xs2nkpOeeLfaMpSyOyAxAYmibZm1ZVIWbIsu9TBmKG2ZCYZYfOzXp/YCtvgosjl7sWUMuIBmgPILxnfbL0KFh0suUaYEK+uOZZ+tMLI7JtjX9plBlXtAA2hEK5bcY7ZMCnlIxORJZOk5DaYY05R7GuSY97nNvSvnNIiJJ1O5dyL3Wvt7HY8lhRuKheUwiSlPnHvagmrK55wGcPUa4K7r2XekWYy4iZo/Q3UYpR7GDCW5TzCEcu/XlpH7jYq/8yDwPKiEAqoBwi8YN4Xc1cCBp+h46qzfx9f8l2LNSS8EECl3k7SPOzNmPic86sUWVIvx3CXCSM1QjQpwuTK527w4WppyT0N1vfR83WDkLicrpdoy/DyS0T1ANMn0eucQ89z592DM5C8/kPTEgXSy3nc/e/zun8UXrDsdJy0aJ7ntgV3A4p7O45wglLbMBMPgVka/atcPmC0TRpbkLBzmuk2YACN3rtwDr/0iVAMHrjqPX7n4Ijz1316CEzYykqvMRAuMNVQxi4hEtAGVu1hQLbRwWFdbJlqsk5W7G1ZFzKncq+ui53olvetTHtT2A596HvD6v2Hj1cx2D1uEdCYXgIFoMpJ7weZB4EekLHIm9Gq+DFXfDZPbAHT23B3u37uteO2eTtUsPSe+ppGm3D9zDnv805RWihOIUrlPMIwwNrx/5a4SNcxQzatyHZvXgeGhkADgue0XlcZLFOiqghM3Ruq1MhOFAu7Vjg2fL9PKQJ57QANQ0JjVVKhyV80OC6opyr3+DHuSjHPPQkUid82KUut7xbMPsce7PhdF8vSi3AXJ9kzukucuomOMak5bpgfPXawF+HZCuXfIhPXt+OKxWi6olhhD1J5+GI1nf4YNGNyWURW151BIj/duhWYwggfgp1y8LIu13dO1TAMOVWEQH4uV44HlxwAATWJhhva/yCWIXC75W6jnbqbURAfice6KBrEu7S7xBeP5zfmOI1T0SRfxzM4+bZk6L8qnqEydWvMRsT37EDuPULmnkXufnrtsfYiJobI2Z/mBDM89ldz5vmkQzwDOIvfAZ2OLKXcjrtynMO+yVO4rgAM/2YEff/PPC9tf8PmLsOmmN0Ih7AfZj08d0KDvwmEeV5SKFkVheCmeO8tibb/tJ4RA48ey57eGrzdIFQqhoN1S1TMgyicUH+cuKivOdI+WkQjKq3HlPpeT3AHgqp8Al/7vwZS7CDv03ShaRijhW/8IuP5CSbmn2DLmAnt9uX3BuyNkchd+fXVDZ7sk/KyfHgqZFh8vTzpuI0HuKccRJC6Tu2rGF8qnMHKmJPcVwMGv/SFO3vkxHHnmqe4b58A8jXux/XjuocqVbJm8Kte1I+WuaDxaJsVz16mLII08gHBiIuufG77WUphi9OQU/8M/yx2SJ84pto5QiOcuyH0ug9yjioOeHp2vV9vHyK7aQx29mQ1MWeuV/glHkBn12+PcBcL2fynKXVGAhS3AYg+/1yAAQCNyP+EC9ji3uZ10n/gBcNM741UZ/YRy71S8TLaLGgejENKsBdU0C0oz4qUk5N/YlKj4ktxXACc1dgAAlheLr6fWpEZfPrUgclVRey4cJit3ES1DUy4qHbxEQQdU1m8Jn9sqJ3dh+zgN0E+/APaNl+cbF1d5hVeFDKs+ZtgycrSMFKvuNg8zH13p4zLTrP6jZYSdE/iR12wkyg6Lu4qs/5+1W4GDu/IfUyhs4blf8ufAZV8HNj6/fUH1W+8CdvwdsLw3/vnUkr/S5wTpyuS+tCeK7snq+JS2eJxU7q60zynx4ktyXwH4/Gtu1QrIOEygTqp9RZgIIteI1nObPUG+qmZAEYo1ZUFVhwua1gxCQmVtZFm4Wly5L+9+CIT6MB/7Tq5xhRNWn8XQMiFIwJjNiHOPPHdXar7hObUoLb5XaBbbbw+5ByGEKg+EcjfawzFr3HJJU+4AcOw5wDMPtWfRZiEkd66+zTngeb/EngtFLch5iYcb1vfHP690yFC99lTg23/InieVu/jOM20ZqUm5gGbygmZ83LJynxKLpiT3FYAgd6defIhVk1T78tw9vgirKr232fM5eRDdBNHYRRikXFQ69UDTojEkzB3zvPB5S2e1ZzyHkfv+/c9EG+bwn9NsmULLD2Qpd88JVaEnKWHXWY7S4nuFIN20QmXdIJS72+ChkCnkXud3kVn/P8e8gNk6+x/Jd8wkuctI+ufCI29Id7JBIoOUEF65kX/fS3uAe7/EnicXerUuce5hwpbsuYvJg3+/snLv1rd1QlCS+wA49NTDeObhf++6HeXk6bZ6jD5IQdLbbioz0EgA2mNXGaHc+ykc5rtcuetmqNz95AVBKUziZt72/88Nn8aH3bfimPVro/1WNsb2VatF9UhaR55BNwhbpvAFVU9S7p7NFCilwD++BbjvxlistVw4zHPqAyh3M37sXiAmQrfBPq8a7ZUpGwfY3UaWZTR/DHusdf/eAXQm92RCkthGrjeT9NyByGZJTuzJEE0R3ZO1cOtFa0TRvhPfr6zc+5lQxxBlKOQAqH/5TTjOfaJr0oOomhgUoAicVgPyZdrS5gAfcB0HhlXJ/FwSocqVUvXzK3feAUozQUXpgkSEi+85UIHM2/4P/fZlONJ4YxiPDiCsc+6J3qvN6CJe3L8b1sbn5DqnwmvLyH1IQZnd4SwDD9/M/s0eHVoDMeXu1oH5E1N2mAODkLtQ7k6dTUKqwZSwjPqBbEsGiBqK5I2YSXruMpK15Pn/TczySXruQFSTJll0zKmzblb1Z9nfgtwVLd1SCcldtmUkT//AY/EJY0rKEpTKfQAc5z4BAHDtzgtfIuwv6De0TYLdjKuWlskJscd9h2GDShQtk9dzD3jYo2qYUYZqwnMPv5MM5W7pKo5eYBfbN/yXYDfdAGOWJfEIT59KC4rLi91rm4frCIVXhZQWVMXfLSliSVbukrXguY3+bRl1EFtGqNE6I7u0/4PGwezFVIC3AiRA7dl8xxTknlarJmnLCHJ3kuSeVO6C3BOZuk4NWIiS38KSyqqRTsxpYZ/i+919N/CZbcC//4W0fWnLlOA4/GznehQi7I8WoNxtOx4W6FVZswfP6a02dSwUssdoGZ9fQJpmQuUKKEhEKTg29+XT0tsTMC+9Hn+37evQDUaQYhIMnIjcnWaXkrFInNMwqkIaErnLhOO2QhUcKz/gNgZbUAX683/DSZEy3zyNxO2lzspd1djYGzkjvPJ47oI0s5S7mkbuie+aUjZpyeGlQrlrGRnEaXHu4vmhx9njk//Rvv2EoyT3AtBYyncBBAX8aNwmI/c6NXGYzoYKxu1x37Jy71XlUm7LaIYVKneaWAsIyThHXZVXnXUs3v/fzw4ngjAhSipc5TW711mJ3Y0MJYlJVu6SFefboaftSraEG7gsIagfiEmxn8iN5GcEuW46NfF6l4m3ur4gck+WPuC/FVm5i16vsc9x5S7fJblNRu7y70om9462jLygysfUSLkjLJV7CYFWLWfHnAJ+NC5X7jvO+Rieefu9UA3u9bo9knuaypUyVO0DT2Df1/4oVQmJtQNVN6Dyi0SOljlw0x9jzd+wCpCqlbMiItgCLRB5+lSymvwci9HyORVqy3g2I60wgiVhywCh0vZIRO4e0L9yH8SWScbHi8XUN/4d8Ju3RSTaSbkDjNzztvoTd30dyZ3/RsSkLcjd91hYYnI8YkFV9tztZU7us9F+xaSrWV1smRTlXj/Qvn2p3Fc3ZJ/d7kDuYUIOIsVbxHFnZudwyvFHg4holZQ4806Q67CkFQ57/O+vwtEPfA5P/OjOts+KOxDNsKDyDFUqkfuGHX8V1ppX85a7BaDw7M6wNrxUWyXo1skHUninXKO+KM9dNeMKNOkDC+UuNfj2CBk8FLKvBdWEehVJPuufCzznhZFH3ZXc1/Wh3FMWVIXF5DZY3L4gd2HLZJVCUFJsGXuZhS0aVVa/HYifjzwZ3vlR4JtXdlbu9ZQ1hSmJlinJvU/UpAU+r5GdnNRsyOFegyt3jyt3jfvTSmhl9PaDFGGDGtFCz10OhXRrTNGkZdWKSUo3rfD4gZeeSKVXeyB3La7ciddCkxoIKAG1u1coFOMvvnAYjxWXIyyylHvgQeHr0i5BvBl2LxiU3GWVmrTGBBkaid6uSVTXp9sWaehky4RVJhvxiUco9zTyBbgt4yWU+yL725iJJn9Ru0dNlPH93ieAnX8fTSZyMxTxnaRFA5Vx7qsb9cXods5rZPvBrXqx5O7zRUbNZD/USLn3tm85bDCtzZ6nsv27S+3KRkSxGNYMVK62aUa9br3ShUAkCFtG3BkQv4UWMVGHBeLkX1DtJ7yzIzy7Xbm3EuGvnCzcwEWF5zV4IIPHuffzm/GcKJQRiJMaEE04ZpeJt7quGHIXzcDderwBSGjLZLT8Uw2mouWJtLafHcucBy54D3tNtCzM8txFmKMc6y++E7kEgsBqUe6EkC8SQp4lhDyQ8T4hhHyaELKLEHIfIeTnix/m+KG5GKVO0w6NjFtyVmoB3da9kNzjyt3v1XOXwgbTCoeJgl9uysRF+QVkWlWoKQuqNRrFExs9kLsWkjvvDevZcGCgBTNXnZUwiUkpuFlHslWd5zAFKcNaE46hwo/tDmLLhEk2fS6ozko15DOVezdyX8/UcZ7CbZ08d1m5y+TeZssklLteYZFIsi2ztJs9WvPAhX8MvHdX9B1rVrzMgYAoc5C2CCuXQBBYRcr9SwAu6fD+LwM4if+7AsB1gw9rvPHQNz6Ok779+uiFDh1z7JgtM7giEOGBOlfuiqin3uMPsnsoJFOfNKXVGnGb8KgC0zRTlXuTRArJrPZD7jwD1m/BISZaxILqdScY+W6k0GYdns0sGTmkL2nLVBi5u76LiiiDPJAtM0AopJ9Q7lnk3k25i+YheXz3MPs0xXOX2/Y5Kco9q7a8XmGfka8v0QbPXGChk/IkFkYY2VHNGABY2svCL2VPPy2Kq8KzpVeLcqeUfg9Ap3uz1wL4MmX4IYA1hJAeClhPHo7f+f/F/j5/79/h/hs+0Lad21iE/2+fDP8mBSj3gF8cusXJXReed5/kLqtcKYlJ4YuTRFbMToPdpnst2DBACIGqiaqQ7RMDAMyszf9TCCcKodz9FlxiwFYsKD2Q+1CSmFRTimBJWVAVyp16sHgNew8kIoxeMWgopNz9qc2W6cFzB/JFzKTVb0ke35Fsmcq67guqeoXdsbWWou9RNEAxU8YeToiteK2Y5b3MGpKzdNPIPVmSYMJRhOd+LAC58PNu/lobCCFXEELuIYTcs39/yu3QhKBJItuhRplKPeOxv27b7sG/fTdOPvxv4d9KAZ67SPAxBLmLBc1ePXepsUVa4TA1EAWVIlI9/Je/AHziBBDfhk0YqWua5ENzWIgIac18/gVVjYd1inPRAhuuYsImFWg5uhKlNesodEFVDulrLQLzUbnimHLnJOESkk5CeTBohqp83OQYRBXFbmMT5J5LuQtyT+vJKil38Xua3RT1Qs1aUNWrzBayl4F5TimiomSa3SUvQsvlBJb3tdfWSU54QGQplXHuvYNSej2ldBuldNvGjTn7So4hjqhRr8v96qbojYTX5ydSt0lQwI8mQe4iQzStzV0nyGGDonCYbMuonFQUiVTXLj/K3vOacFiLbKh6oiEyAJPaOEJn8L+Ut8Rrx3SBZrB9itrwKid3T7WgBznIPSXOvbgFVdmWsRm5S41GhIXhUQ+aakKjlNkyPZx/DCFR9bOgasfrqMgWDRAR31zi9SREw+48i6qCENMyksWCqiPZMrOboto3WRODqGlvL7GJRjUjck+bmOQJsRu5K2pUTVLcwYh9lso9xB4Ax0l/b+GvTS3kL03ugtRqxhNtPDO6JV+kM1AKIHexmGlZ7IIRVkavyl0sPuqKDl30/5RammlcuSt+uy1Qbe2DQ9hxhXIXce6B68AgPrYf+z/wtvde29OYdK7cKb+4jKAFXzHhqhXofndyd7nvr6s6dNE8owArrG1B1XeYmp3dBGw6Ddhwcpg67/oudM2ETimzZfqF1ueCqu+xkgOaFSUrSTXmAQDP/UX2uPaEzvsKlXsOcvcyrBWAfTeqwawSodbnjmFRL57dWbm7XLmbc+zfIl9QTVvLCG0ZO579ai+m2zDijuLYF7BH0eu2JPcQNwP4DR41cz6ARUppSnzR5ODHt38O+396f+b7Om/i/LNgE+aDKGqiVY+Tu0zmh9QNUDLCBXuC14JHFegGu4iUkFx7I3dBeroakbsrjU+Qu5pih8w5++EKchdt5fjEICY4szKHNdXudWVkhPsSYwtseIoFT63CzKHcw3NS0s+pb3g2I54wPNFlhFfdALz9duCKf402DTxomgWNclumX6iJUMjafuB/nQE8/t0uYxWNKQzgnXcBV6Rs/4K3AFf+AHjuyzrvy1oDgPS2oJpV0sCYYR678NnFXYNTy/brZc/dWmBWjNh2dhPaIE+IyQifpHIHIuV+/IuASz4O/Mp1PHFqlZA7IeTvAfwngJMJIbsJIW8nhFxJCLmSb3ILgMcB7ALweQC/O7TRrgAay4dx8n/+EfDl12ZuYwQ2/nPml7Bw1Q7smP/F8HU7odwVfqH9kJwFVzELIXfitWBDD+0OLVxQ7W3focpVIpXrSBOETrnvnaKY1wUH4SnsQtJUHh3BL+4Wjw5Skm3dcsAwuXLn49Cpg0A1EWhVmLS7go2dU5HkLroZCVvGqUVWgTkXU4Vu4ELXq9BA4VX7XEwFWJ11RY/Ies89wJEngf/4TJexCnvEYrbRMS9o34YQ4Ogzuo+hl+JhnRpuA6zKZO2ZSFHPHs0e7eX0eusAu+OgPju+UO4AG1MqWUtWVrLme6rHrkb7O/9KYO5oHis/HZ5713rulNJf7/I+BfDOwkY0YtQWD6MKYCPN/kGbsEEqc1gzN4Pz3/VF/N///T9x0d7PwWnFf1Ca38Jjyok4/X134qlPXQC1AFuGeDYcYkDQidqvcpeIUCMaCEiMCENyDzi5SOsJVbTgcnInigKHqmHvSrvBLt5+yF3XdPiUhARlUBu+VkGgVWChu5rqNmH1DdGHVNz2C2uguq5tUy/woJsL0I1ZuEefO9hxNSsiGkHyQXomcDSADo2v+0He+jLhpJJx3PnNLGqlTbnXOyh3/hvymsyGEVaMmBiSiCn3RC2iVI+eT9bCfgKixKkpQJmhmkCe0rIWtUE19sOzDB0bTjiTfTZR3EoNbHiqhVlTg08MqLQI5c4SewTCRcge1YYgQkNlIY26osMJ0pQ7Iws3UWrYV6Ix+FCZwgLg8LsXxexeDTIJRSFwoYHwi90EU+5Un0EVra79RGXPXSNa7LWB4LUY0YqknyNPsseZDW2buoELTdWhWWvh9ZudKqBJRCNshm4efOh9F0jueZS7+P2l1XMHWGTR4Z/xUsNWGDoas2XalLukzmXlLi9ky4iFQiZsGXE8GSI6ZkGKetLM0nOfViStlSRo4KNCnNhtnigF4LbiPygtaIX2ha/ohZC7ErRCvxuIFjTTeph2gvCnNf4D11U9tvhoghMst0Oajfj34itSI2ioYQy/aCWo9UHuAOAganJsUgdUq4DyaIugS6ak7LkTQmAoRjHk7jTY4pteYckwgtxlxcfhBV545zDwsVUpnV7YDN1CO7MWJ/tFXnLvFAoJAMefz1r7PfmfcaK2a9kTkmylWAvR933UaenHkBe8xfclSD0t3+CXrga2nBu3qLL6sE4gSnJPwG11LlDVarTXqRBE5iXIXSwIAkCgGNAKIBrFt+PkrkuLfBIWH7wDfj27WqVsYQBoI0KTRr430D7puXoUv+4RLbQLhHLXrf7I3SNcuQc+DOIBmhVaPLEibDnOSVf1YmwZt8GIhhCm3g/9lL1ebVfugtw1osWij/qC7P8Km6EruQtbxuq8XV5U1wGNHCWtOyUxAcDJr2RrCHvuZQlW4i7IWY5KSyR9dFltV9cBZ74RWHcie0yDWPtw6tH3JZK5KinK/eRLgMvviMfMl8p9epEk6CRENAiRPGVBZF6icqERtODzi8xXDWi0M9EEnovm1Zvxk5s/nrmN6tvwJHIXceZyyd3m4kEs3HgpnvzrXwEALO95BD/98u/EForaiFCJlCYNAtbcGixmHQCcBLH6euRh+lBB+GdF1Uq1kj95SYYLDQjcyAbSLBCeJt+qd27Y0emcOiKtpjcANA+zjEi3ERGHMROViU2JE3cDF5qiQVO0wZW7XAhL2AzdJoxO8eb9IG/ZX98BQNLLDwDAzHrg+Rez55W1cSIW2aTJcEV5TaOyDjjhAuD3fwRsOCn9GGIyaB5md1tyyYE0WyYNqlkq92mFTNA0xeMVClaRfohGhT33E760QW0EKlMjVDGgofPFvnR4Pyq0gedu/7PMbdTACRczgSg2XP5BHniGJQyfUN8BAPjpTR/BCY//Hzx6963hNrI/LR6FrSHXqg/JPdG7lUpxxg6MMPvW59+fWelTuUOH4ruwxSSrV6BwcnfqnRuRu4HLqlxygjEUo7t63rMd+ORzgfu/1v7edS8Brv059jwkdz5pqWYqYYS2TN6JpRM0iWiEzWB3tg2LV+68eNjHjgMO7OpwXB4u2in88wWXsceZjXFbxqmz2uxJS0e2UlIWr9sg1jhai1FDDzGeNcdlf06GltGHdQJRknsCvkTudopFY3NbRjEl5c6fJz1hAw4C7n0GigG9i+feWOoelaBRO7aYqRvttoydqOQY1Fiph0Ytel0Qj1h4lMmoxc+7Ri1YYFX22r4L6VbWUaBYqloAACAASURBVEyoPGQy4N+f0ady94gGJXBCewd6BWSeRUfYh3d3/KwbuKFqB3LaMvt/zB4fvT3+euBHFQiByP8VJD93VCqRycp9YFtGrk8uyD0ZBZJEWISrIHIXxcPsJeBnP8jeToSLdsJJFwEX/xlwyZ/HbZkkESePnXyeBd1i5906wu4G9Cpwwf/DwjA3n93980Cp3KcZgUTuzVq7DeC2mD2hShX1hHIPEmVpLWqDCuWuGtC7KPfGcg5yDxz4knJXVZWFD0oq0U3444EoZFWLbq9d3w0XHoE4ETpcNS+ROSiEwnWa8BPkrlSiSBCHmGG5Ap9PcGa/tgwxoAatWNTN7CaWSbn19t9sr6Muf9ZPkHse9Zxs3iywd2f876RyT6b0c3iBB03RoKt6QZ57gty7daQSyr0b0eaFHBGUVh5XQISLdoKiAi98JwuLFLV6hHJPC52V1Xpa0lIarAVJuVeB018PvOfHUfZpN5TKfXohq+9mvf1Ccrk9IeqpA6yuOQDAicidBgGPqpHInXa+2Du16xPQqAM/sWjlQouRk5sM5+SCyJdu6ZMqV15QdUQTboWp82ajBi8Z5ildLJ5iQhXx8HxytHoo9SvDUSxofitc2FaNKjZuPj58v/bkjszPuoELQyKYXBEr4kJOEvEj347/LcjdzCb3gAbwqAdDMaCRgjx3P0W5dwoJFd58jsbkuSAnQSXLHMsQJRp6gTHLzifZ8FqAEODFfwi88F35a/RYa4DmER7hxPep9EBzmlXGuU8rqETQaTHvLicvXWr8bPFWclRS7o4decYAANWE0UW5d/OUAeHjx8ndgxarX51c2CU8Bl0+NzdwQ78dQCzO3eMTXFNj5O4066Ei30nZYtb8llOj4ykWdE7u1G0ioASVPj13V6lA85uhx6+ZM9gwX8WbnA8CAPY9me37CtUcO6dut9hZUSi7/gXYFJ1jaAuIlm4p5C6Opat68baMXMLWbbcLo0EIO6v3JLJUzB8D/O5drLhWh74F8JzeF3GN2agMcNZ4X3E1cPFH8++zsoYvqNaigmW9QDWmJkO1JPckpOSHNHIXi6Yyueu6wbI0pTosdoNtRwS5awY0EmT2GgUAv0MvVgGDum3k7hItVtrAT6hsRZCMdG5p/rRYUBW2jK0z68Vp1hHw8z580V/gcyf+JU58/unhZz3FhM5r0RC3gSYMaFpG1EQXsAqQdqjcNasKVSF4/1tYOYhaSk/XzHPKY8uk2RyUAvsfiQpsAVGMdTVB8hLE5FjogmrSlgE6L6o6BSt3ANh0CrNFOip3u/fEKWuBE3G9e1eovFi7FTj0ONtvnkXYJJJNticYJbknIRG0mxJX7YcLhvGLx4YRyx60BcEagtzZD99xsgtg+bxdn0OziVGHA6rGF8s8Hj4oQKWLnwY+DJ9PNPK5+e22jFCaYgLzTIncOWmcdcpJ+O3f+A0YWvTT8VUTOic24jVgk/4TaHy1AoM2o0bgPMz02KPZoiptJu5u/usLwJ0fST8nNUcSkyBNOU/A5Y2cZb95jqe8n/paVgXy1Ne07UpMjoZqFKPc02wZoPOiqtiuSHIH2AJ6J7/fbabXe+mEmQ0sDNWppXvu/WDjyaws8MFd/ZF7ssn2BKMk9wSIrG5TomUCR4T6xZWGTcywUBgA2NxWUMXtJlc1tp2dPk6bTBn5yCZ3izqgiexDhxix0rzyukGrvhyReyflLtkyrsP2FVgsFM1p1UCdbC890CowaFT/3cYA5K5ZMKkdTjAGv0Oam52BTXVQWT0GPvCd9wDf+yQQ+KlWU9eSv4IoZcIUJW6r64GzL2P2jCCKzWcB77qbkUgCYUkHxRjCgmojCg3sRLJOjZf77e/OKRNmF1vGbfZuBc1sZIu0dq24yWjrS9mj10rNIO6KKVpQ7Vo4bLVBbk7h2ynkbqdHgzCClZRxKx4yqXBClmPIkyC86bIBl1kDiUUkKmVtxo9tQpVrjkhE1agtwgwa/NyibTrZMuHkIBpQ2A1A9E012ok7UC0YvLCX6tbRVPpXYawCpB2OQdwhWbqKA6iCyASz597oee3Z/myZtBBDkbRTWQe89jPds0LFrkTBM7WgBdVY+YEaK5jVPJxN7qIaYlF+uwxzHqg/nv2+28ifKCQwu4kpd93qj4jTcOw5rJSv18wXPplEGQo5vVD9FuqUlwxIIXeRKm3NxMndJUYYDgggrBCp8dtNwssEeE62cld5IwOVULhuu3qwW3xi0JKeuxW1xQNA5Iif2hFUaJOfW8KWSapcTkY+n4DUGUHubNGrBROK2v6ToZoVlivQ/AZspcfb89i+KrCoDT/lDqlOqlCdiNwP7/qv6L1DT8Pxnd5tGUGUsu0hyL26nk2wOVVw4QuqxgxT7JSy8ck10JN47A7go0ez2i1F+dcyzCHZMm6dfd/9EHEaFCWqPbP2Ob1/foqUe0nuCah+E0cU5jUHyZrQAOA24FIVhpEkWDNmjYgyBioPkyShcu9A7m508bQa7RewWOgMmwyIYylGWL0RQKyZdKu+KJF7tE3Lb8GUFsBkW0bE6+uc3H27AcVropnlpWsWC/ukFIbfgKsOoByNKnTig/B8Aksi95ZShSpFihx6MmqocvjIIdi+DUtaj9AULX+0jJ1hy/SAZD35gZW7OcfWUpw6895FhE7aguoj32aVOZ99KL2OyqAwZ7vYMh0iXrIgN/EuSrkDLLrm6DOBrRf0/lmVr3MkWmZOIkpyT0DzW6ipjNxpKrk30YLR1huUxWdHM74IHRTZq0oO5W540UVrN9uPLcIDiREnd1e1orrriJO7Uz+CKrHDcxNoeS1Ykr1jqIZky/A+rfPs4vPtOlS3hhZJV2aU+6VOcxlGMCC5c4JQmgcQUAKrGu3LVqrQ/eh7cZaipJpWbRFNrxk7p95sGen7bvZH7mJyNJSCFlRFiv7S0+xRNIl2UhS0rKpTShEPDHOOTSpZpOc2e18UHRa5H38+cOX3+/sewobvBVQTHTFKck9AD1qwtXmW9em2++OK10yNBvEVM2xNB0QFtHTuGSu8j6XfgdxNL7pAk7VcgKjLkZoop+srlTAUEQA0idztw09Hr0sTQFLlpin36gK7OKjThOpl2y2E1/SoLx2GFTThaQMsjnFLwaw/jSVUYUghlY5aDReHAYC0FrFEebnlxhJs347djRiq0V25CxXs1qPkoMZBAKRnBSx77oUsqIbkzssgzB8THzPAbJv7bozKEAPFEqU8FuqnXhMA+ltQlbNO+7FQhgG5yfaEoyT3BIygBU+toAmrveA/WE/RNHL3VBMajX4Q4YKgKRpZc+XuZpN7JaiHYZDJrk5AFHffRu6qCUNqQ6f5TRwEI1z3SETuhkTuTa8ZI0JTM2GHCTPsAp5dy5RV4Dahd/DS1SojwcbiQVi0iaCf5BEOZYYR01zzKdTITOwOydNmwsVhADDcZeyhbALym4ttdyOWasHudpHK/rWwfBoHGbH3GHESC4UsYkE1qdxnN7FKh7JKv+NDwDcuB3ZH6w9ppYgHRljoK8N3dxu9e+4bpIijbs26Vwphw4/JX1QtyT0Bg9rwVQstYkJJJfc67BR7wlctGJJ6Dqsj8tBBlSv3TrZMlTZwkLBwt7QwzJDcK/FwRBaKKHVR8ps4orD9kOV9AACbarE+pLZvx4iwolXgBA78wAd1m/ApwQInd7gNGH4DXobdos+wYzVrh1ClrYHI3ZxnxLTJfRoNJb4w6OszsKRG2aa/jMMat45aS2j5rdjdSEWroOk1QTv5p04NYX0GYc00DvWlfmVbRld0BDSAH/g97ydEktythShlX6D2TPvn1hzf/tqgMDqQu++y8g29krtRBU54KXDstvy1Y4aN0JaZfOVehkImYFIbgVaBTUyoXvstqOEto6mmxHorVhjrDQBUFNDipQmEcg9SomDYByhmaB3P6Jux2TvQ1tUJiAqCGVac9KhWgSn1GDWCJhb1jYD9BLQGu/gPkzXt5C4RYZW3DWx6TRCviRZMVE2T3Uk4TRhBE0fUY1OHbs4w5e4sPguDeKADRGtYC+wiN+GgpSbJfRYVRP8n1aCGRmUTgjoBbdXQ0lowpUiiilaBT/22mjMx2DWWbbr8tETuB/si97DSJi8cBgAe9aB2yFvoCPE9ip6t5jx7TbZlpLITeMm7gaW9rFhW0RATTZrfL763fkIwf+NmNjHkrR0zbIi72SmImCmVewIWbXFyr8Ti1gUq/jIcvZ3cqWbCgHQrx3/wVa7cNV533c8gd7tVg0F8NHRGKn5KSzmPq3k9EWMfaBYsaWIxaRO2uR4BJai2WGOJRXVdbJvk4mOFR+CE5E5MEELQggl4TZi0CT/j4rXmWFSNfeAJAJEH3w9m10SLbLY2H3uPGnOsBLHvAZRiltZBKmtRh4XAWW6bsORzSoXvMoUmQgyFKm0c6is0LxnnDmAw3118j4LcrQVGsjLBNqRGIyf8AvCrn0ttIjIwOtkyg9SzISS77+oooJXkPpWgQYAKbMCY5dEvaeReg6vPt70eaJUw1hsA4NbRojo03ilJkHtWI+vaIqsI6VQYuXl2O7mLqo5mMktUr8AgPgKPKccKt0YasLDgsYiSmrkJVWKD+h4opW2LjzIRKm6UZdokFlS3hiptIshYKK3MMyKkB38CAFDnNqZulwcL66Lb85qVICmuZL3mEuA2ocMDqaxBAxYCpwaf+pkTVioEKYkQQ1m5z/Rhy8jkrhRA7iLa4yAvlmbOsdeWonUU1A+wRhfHv4hFiQwLnchdZA0PIwRzpSEyZTsVZ5sQ5CJ3QsglhJAfE0J2EULen/L+2wgh+wkhO/i/y4sf6vBht+pQCAU1ZlilwxRyn0UdvpmiTLUKLDhh9ybFbaBJpGgUUzR5Tv/RiEYdwcxRmdsFnNyTJQAIV0yiBaAlyJ1UsClgyq5VZYWumvUluIGLgAaZRKj4LTi8ZnxNmYPuLLIJIyPUbY5787NLjISEb94P5meqWKZsLP5s3AYi3I5q1BbhiP6w1gKaxILPcwTSJqyG1z5RAojIfFZKDqJ0YFtGxLnLr/UFc57lNBz5GfvbWsMSdHb/F3DHn7DXGgeBc98B/NatvXvePY2lE7kvRuOddMjt/yYcXcmdEKIC+CsAvwzgVAC/Tgg5NWXTf6SUns3/faHgca4IGsvsR6oYM/DUSljGVsBuLmMWTdBKyoWvW7yxBa+x4jbQghS5McsujrSSBgDQ4o06RNchP6VMAeVq3ppJknuFj6+BwHNhERfQq2iRChTCFhMpr2LYqC2ixePdsywMzW/C4RNTU13Agr0XGgkAM12ZzVarOEAXsLn5GNvvQv+LY4pCsI/wySHRK1O12Hk3a4fR4NUh1eoatEgVLo+iqUgJXl2Vu/CuRVEwp8b++XZf5N7ipQIszSpGuRMSjW1mI1vs2/Z29vf2r7CIDntpONExSXQid5Hc1Gv5gXGEIPdu7QwnAHmU+7kAdlFKH6eUOgD+AcBrhzus0UDEkSvmLDytEgsdBIADe1jXe3XtlrbPCoJt8WqQqt+ErUgLlrPsh08zfjR2jVWENNYytZqWQCVeq87EFZLCFXWrsRw1GDFm0eI1XgJKoC+wGOlWfTEMeUxT7g2vAS1ohX1abX0eR3vM81Vm0+0WQggOqBsxBza++Y05+1Vm4B82/B6+7Z+P+Z97eex1jbf2s+tLqC8xctdn1sJWKgh4clNMufP/k2ZWbLYgKlG+16nHSw/0CDGJVLRKuKDatXBZNwhyF2PcdArw0j9ireRq+/hYC0rd74Q8yt2aBuUuFo5XgXIHcCyAp6S/d/PXkng9IeQ+QsjXCCGpVzch5ApCyD2EkHv27+/QsmtEsEWSkDWDQLVi0SUAsPgMuz2ubkhJuOAzvs3JVfWacKSQScuqwKNKZrlWh9dyn+PEmDoJeA00qAlVjUdfiDjzVu0Qmrx/qmLOwOGhi3VY0HlEi11fCm2KLJWr+S14XNU7xlq2iAnA6GC3LJvM2limFWw66pjM7fLgrb/+Fux9xXXY9rz4fvQqJ/fGImx+p6PProWjVuHxEgtVaVGvq3IXpCQnBw1I7pqiQVf0cEF14Fh3cfcihzeufQ4rZiZaAQ4jIzUJzQIULb0EQUju/S+kjw1CW6ZLO8MJQFELqv8EYCul9EwAdwD427SNKKXXU0q3UUq3bdzY/6LbsODwkruaNYdAr8JCfMW8eYDNcfNHtZO7JiXyAIAeNOFKtgdRFDRgZSoC0ahjbt3RsKkOkjIJKG4DrZQEKnOOEVFz8SDsuiD32bAMwGFlDQxBjPVF1Pi+Z/Uo6kYmQj2w4fGEpUC61e5ktyzNPx8AcFBZB73PRh0Cx6+v4h0vPRGKEg+PE+TuNpbg8JaE1uw6eFoVLlodzymGr74VuPX9AK/CGVfuovRA74TZ8luo8J654q6o5WfnNeTCc17CH18cvbbAtZOoipnRz7VQEMIaTdeebX9vmjx30UZxCmyZPHHuewDISnwLfy0EpVRuj/MFAJ8YfGgrh/u++hHQ5b1Qn38RAECvzKGhVWEllLt3hNkTGzZvbduHnMgDsESimhEniBaxYoWvZAS8CcXMwjrUiQWSsp3i1tBMSaCyeJkAe/kgWtz+Ua1ZOMYC0AQWtU0wBTE2l9ByObkb6UQot/ILZqOOQ/NHnZg6dgBQf+5VcPZ9Gf+17tXYmrnVYLBmmDL0m0th1cjKwkYcUWe4/aF1J3dKgYduYs9FTfaZDay9mlNj0SdAX1ZH02uGx5zhiVz1QaMuznwDsP558V6ma/jluPse9rhSCUDzm+OROgKNQywMUrfa35s0iDu/KbBl8pD7fwE4iRByAhipvwnA/5A3IIRsppTu5X++BsDDhY5yyDjzoU8CAH507H8DABiVWVC9CoP48F0Hqs4SYJTlp3EYc1g72x7nHsZ6c7ugGtRwSH9ubJsmqULNiNygrSX4lGB2ZgF7SRWK264cdHe5LWsTAGbW8BT8+iG0lhi5G3Pr0Vg4EVhkmZ0m9+m95hJqbopyl/zpKm3A47H81uaTAR6Jt/nY7MzHF73k5fh7+m94xQuel7nNoKjwictr1aC0mHKfX7Meu/Uq3MABUMGM1PRBTswK0ZI6OS3zn6xIDnIGtGXcZvg9imM3UrKcewIhwJZz4q+JAmJ7trPHmRUi97nNUVimjPqz45NhOigUlRF8p25XE4Kutgyl1APwLgC3g5H2VymlDxJCriGEiF5jv08IeZAQshPA7wN427AGXDSo1Ene2cfmJHNmIVykbNQj781o7MMhJf2ir8yz1706s1dm6TKCRMikrVSgeumKQLEXsUxmoKgKbGJBS1F8hruMVkp27DxP/PHrh+AsM+VZXbMR6nHbAADPbHpJqHqDVi1UkzK5CzKqOcuYpXUE/BZ70/PPQ0AJ7qI/17EvqqEpeOvLzsAxa4YXjleZ5dU6W8tA4zCWaBULsxVQYxa2wtL856QEM3FnUpMnyrqU9PPMgyxG3JiJmjU3DrLX+vCPm367cs8MwxwEmskad7h11gTa7D8juCfMH8N6y/71C4F9D0Sv155ZuQlmJZAs8SDDs+NF2sYYucoPUEpvAXBL4rU/kZ5/AMAHih3ayqC2dBCCDqz9rD743NpNIGKBtLGMuTWMuGfsZ7FspP+IZ7l69uqHEHgu5tEI29QJuGolVtVQhmovokZmsQaArVahp2xn+TUcttrXqqszs2hRHWgehleb4+PZiONP/nncom/CBeedC8ojRmhzEct8sUi2ZTRFw5w+h8PNQ6yeOg87fM5xx+Efz/sqnnviSRg1ZmZm0KQGlNYhKK0jWMQsjtNUUH0GDe7Pz0h1bUzVhKVaWLQltV6XFvKfeZBFgRDCG2PUWMZndV1f6fBNrxmGl4qF3YFtmSysOY5Fy8wPtnjdE44+gz0++xBw75eAV32K/V3bD6x/bubHJg7GTLbnfst7ge1fBt7/5NgvIK/6DNXF/ZGHeHTjEfiUYG7Neqjco24sHw7fX+sfgF1JX7wSEwCaR1A7wm7tSSVO7o46AyMlMQoAdHcJTW65uOpMKrlXgxpco125E0KwROag2EcQ8AXBhXWboKsKXnnhBZirmJibW0CdmkDjABb3PQQgToQAMG/O42CNkR+RFlLf+MqLsO2U0VftMzUF+7EGWnM/NOcIagr7Log5i5qiQCVKLAIIYOeUSe6Hfxp56yav2bL8TBR+2CNkz33o5L6Or39sOmU4+0/D814R1WAXiVUAm2SmxZYB+G8hI1pm+5fZI8/GDrHv/2/vzIPjqO48/nkzo7kkzanDsq3Dh3xIvvAFztocBmwgyxWHQKiwFEvC5oIlkK0km81mk9TusiRbVHazQCWEgiWVTUKlUlAQsAkJGIfDNg7Gxhe2LNuSNbpGGmmkmdEcb/94PdJIM7JkXTOS+uNyqaf79fRrtfrbv/699/v9DsHObw+kjc4BZr24Bzt8/culsp2AKMBoNJJXoMS6N6DEIBoJ4SVAPGWAMZW8PDNBaYNwgO5ONaPAVDB4UC5usmFOZBZ3S2zA5RI12rFkaFdADwlz5hkJ3QYHlnAbIuQnKG1YrYNFTghBh3BT0PUx4Q9/QZ4kLZmW0+KkUytUYbTnnlUihKDT6MESasXa10HIqH4XBkshQWHAbrSlFVFxWpzDiztAoWb52tzKau9qHPBpXySDxD3pc58MtwzA0uvVz2RB6KnAUQYPn4DaWwd87+EuVdd1MjJRZgu7d2DsZTiGumZ+dSe88xNoOzF5/bpIZpe4JxIwJAVrqGNwytRuoUTNpk37i3QpH63fpy6mwZkewDSwbwGmSAc9AbVP3lBxzyskX2Z+3bPGg/0JyeJ5BVgTg0UhEYuSTxiZKfUBELCUURjxYQq10SkyPwC6TB5qevfTYTTgiadHTjrNTgIRNWaQZ3enbc8Fgnle7H1tOGOtKvMlYLAW0GE04DSl+56dZieBvhRxT960SVFPWunOcug8p5J0jVHce/p6+l1dJoMJq9E6/gHV4ai9Fb70NmyY4kwfBgN4q6GjXvmfO+rV+lzJxz4R2IsGJ2TLxFAjQUutTeBcetssMavE/aPHbuTMv68nER1I3hXrVhfpPErMO83qZs93qc993UoMOrQAJpt3eHHvMBVjCzfT06ra5hcNtmbiNi9O2U0ili6s+YluYppVHs/LH5TaFqDLrx5CYpgpeuH8eZTEfRSEztNhzuxW6DWrtxG/0Ygnnkgrmea0OOnWBh8tBbkp7iFrCcUxH56En4iWL8dkd+E3GnFkKCaSZrn3toPFCU5NwJMVgFzlKuoz3Dmw7SLp6uvCkfJmZc+zT55bBlSeGUMWbuGiJSqIyl+nXFsA7qqp78dkkV8MPRks91SDaGge/WRh9q5GcoVZI+7RvjC13XuojNZxeNfT/evjPcoN0WRXQTi9BWrAsr8KkXaRQ23Kci8oHv71s8c2B1dfM9F21bZo3uBBJpFfjFHIfrdNKgWyh7hZmw1ideOgl3jKQyjQqv5oTM7MbiHprCSfMEuix+i1ZR5kCxUqIfMbDXjicYL+pkHbnRYnQS2M3+aegsCYMdDtWIIZFfUZd6nzMRQU4zcacMj01LEZxd3uhpLl6nO5lkkx1a1Qkil10oWJJ+IEo8HB4m6yT55bJpsko2bbTihfszAowZ8p5HtVhOrQqmmprpq0YC7NUArmTuT9rBH3pFsFwPbBgLiLUDu90kLIrQJaZNkaABwOF1FpRGoXtK9Dias3QwBTkmjBfIoSrdB5hm5pw+UZHMRkcqi3gUDbYFHtC/diFdH+POiiUHuLSBnsDbarZZs7s1VuW7CxfzlSlFmcDFUqyrHNYMKdSNDWeFK5qjSLxGF2EJRhEoC7dHz5YSaLvLkr+5flnFUA2F2lynKPp89wcVqcdEY6B6oxJTM+XvcI3PYMVF+r1pdfOrBTasDQKEnOQCo0D56KGZwB86XTSIp76wk4/xcoXn7xxbFzmWR08lDXTKorJnU5FhkoydmTIYI3S8wacQ+0KF/YaVFORV8dUvO9G8MddIlClu/4Nn9Y8l1WXf8FQKULCIjC/mAZ0dVAj7TgdA8f3CK8izCLOJVtu/EZy9IG9yyaH79niJ8/0JaskaneFsxaZsiutoFXvEineiAUFmV2GSxZOzCw5qrZmrHN2q2f4deLH8VvtuKJx+lqOsXZH1/Lmf9Wg3Ml9hIkcNZgx+nIvQFVgJKll3I0Uc7hRBXzl6rgnkJ3MZ0GAwXR9JJ2JbYSookondpYQn8JPXO+8lsnr5FzPmx+CK76pzHNlunqU2kfHCkh+C6Li45Ix3C7TF/M+crHfmYP1P8ZKj+R7R5NLMkZQUP96knhNlkHu2VS22VKz5AlZo2497QroWxyr8ciorQ2qNH+vEgnQaMTr8fLNXc+REH+wPTAgMGFJaye3uZgI63GEsQFfJxFy5RlXCpbaXOkW8/5mtUd7vQNWt+luVwsLuVysbnVz54Ut0msS+3jKs7s8y+w29m58Wl+W/nPrLjkrzK2sVlM3Hj7XfTJKO54gljjQSoC+6ns3Esk2EGpXbliTua50x5MucLKci9fMD/KPcZ/o2aumq5pskRICIEjml4IpVTLj9/Sq910F6qPes134Yp/GFO/+sU9xS3jsXrwa7OPZhxLtkPdGxALQc0MSxKbjB0IDPGfJwPgSmsHu1+Gs+izzKypoRrtVNaxseoT4P8d7edOUFKxFGusk5BpmNklljLcYbVfYbiJTvOFA0aqll1Cs3RTKjpg2SfTtjuK1P6xwGC3TI9fHcPuVdsLvErcIyntDJ1n8MtC3M7hc2Zvv2Hk2pm+XvWQyI+Zmd/wcv/6xmP7KFmg3grOmHPTagew5hl5/v6rMBoEZpN60Aai6oYqjqSXRks+sJp7m1nqWarcMmMooTcSyTcDZ8psJo/Vgz88Q8X9si/D6bdg/nqo2pzt3kwsycHh5GBxkqRwl9aC79dqQoIQA0Lvz9sg3wAAEg9JREFUmKdb7tkg0dVETBqYU6Os2mCzCkIoiHUSNmeeGRIuKKck7kMmEpTEfYTyLzyLIs9k4vxN/8eu2kfZcM3tadu9xXPpkVakv37Q+oj24HF41fe7StTPeGDAwrcFz9Bsmjtui7qxW3tLSLgpSQz8IQZOH6BUC9Bqsud20YUyp42SwoEkVY1BdU7lvenpaJPi7uvxqSClaE+/+2siae3VHjC2ge/22rz0xnqHTzk8nXFXwpffhpv+K3eKW08Udo+KPk1O80wSbFGzYrzVqqBLMv1xUvRLanSfezYw9fhoF27KKqqJSQMJfz0yEac40UY0P7NFLt0LsIsIZw6+QSG9JLwjzwi4ZN0mtt32dxnzsBiMBpqMZdi66wetT3SeJy4FnlIl6oWFTlpwY+yo629TFDlHt338g5xJITTm1wJwTCykXTrAdwhLKEKelLTbp9fg2PmgejjW9DalTe8sshVhEibVJllo2jnxg8VJt0+xfUDcPVb1htARnoF+95mOe4Ga6plKT5vyxydTLCct9n6LvkYlphumTvJUM2vE3RJuodNUhNlsptlQjKnrLB2tjZhFDOHK7Md2LlbT5AJv/AQAd/WlGdtdDJ22ctyRhkHrzIE6fKKkP6pUCIHPXIkjqN4ugv4mSmUbkaLacR+/MdiI2WCmZMPnOZUo48TSL3LOvBBH13F8R96mKhqlxTa9LLGGYANWmUdpojvN52k0GKl0VFIXqEsR9+FjFcZKa6gVp8U5qBJUUtzbQiMExOjkHiU10PThYGMhcE5F6SZTLSQHVXtaVa3bZCBXjvjdZ424O/pa6TErq8qfV0ZhqAG/VjbP4s1QWQlYtGoT3dLG6sDrdEk7C2rHX10+7KqmLN5EpHdg7rWj9yytlsHWZLBwMXOjZ5CJOOcOvglAwcIN4z5+XaCOCkcFmy7dgvWhv3Dj7ffR5VzG/Gg9oWOvsaAvRovIUG0nhznVeQqvQUt73HIszXpf5FrEqc5TENCmw44xAvVCNPc0U2IfnF9lboF6I0y+WehMI+atVVMhUyNO/afBs2hA3JMumO4mtW7o+iwza8TdnfAT0fyvPfnleKM+us5+CIBrfubkS1arjX1zPwfA+2V3YJsAd4V14SaMQnLm4G4A4n1h5kXP0OMYkge9fAP5hDl3eA/hIy8TlFYWXXLVuI9/zH9MDSwC81xaLpZ567AQZVXjr7HFS/D1Nk1e2PwEI6XkuP84cwtrSEhB97E34Bc74JcDYx6LXIs4132OcPNhlc51EsS9vqueysLBRsL8AvWGcKbrTKZddHKZ+ZohVb9H/YxFlNB7Fqa7Zfyn1fpk2uMcCWSaFeIe7u3GQY/KgQ3E3QvxEsBQ9zo90kpF9aph97383v/g0B3vseXzP5qQvixYcyUxaaD7o9cAOPvRO1hFFFPVpkHtFm66iYg0Ed31PZa17uRgweU4CsaXt7sj3EFLbwvL3IMfZrVbP8sJWU6vtFC88FMAHO84Pq5jTRWNwUa6o90sL13NQbkI197/hFOvw4lXVZANsNS9FInkWPMHkxKyH01EaehuYIFzcH4Ve56dElsJZ7unR/5vnRTmrFJGwNGX1GffIUCqyGabR7lhkrNp/HWa6Cfnx+uW+5TRckYV4TB6lGXl1oJ81nT9iZOWGkym4WeEmkxGVi5bdsFCFReD11vMQcs6yhtfIhGN0LH/efqkkap12wa1m1M6l91Ft7Mo+D5RacR1/bfHfez3m1XNzVXFgx9mXmchifveZP9n9nLX1fcCsN+3f9zHmwr2+fYBcMOST/Ck/DQJDPgsmsieeAWAS0pVxOn+7jqYty7j94yHc93niMkYVc6qtG2VzkpOB06n76ST2xgMKsjt411q1kzdG2p9xSa1bc5KOP8BdDWpfETexSmWuy7uU0b7aVWEw125GoDqNZs5KZTQ96y8a8r7E153HyWynY//Zwc1jc9zwL6Z0jnpM3Y23/dfvLL2p5zcsZPaFWvGfdx3m94lPy+flcUr07Ytm+fl8toqPDYP1e5q3vO9N+7jTQXvNr2L1+qlxruE2KJrWBt+gs2B73EkUUnfUSXuHquHxbYS3rXkweJrJrwPh9tUVaJlnnT33oqiFRz1HyUST5+Dr5PjbPqKKrv32y/Avp8rYS/UXDLz18P5A+oNEaByk0rBYPOkz4/PErNC3PuaDhOXgnnVKwA1H91670u8e/XzbPrre6a8P5uuvY1XXXewtPMtfHgp2vHDjO1sFhPX33Q7a1cN7zYaLbFEjD+d/ROXzrmUPEN6gq1Urph/Bft9+3N+lkcoFmJ3w242z9uMEIIHrq6mz+yi2FnAa4m1mBreU66ZWB9XdQXYZ7XSMufik4KNxIHmAzjMDha70uvHrileQywR40j7kQk/rs4k45gLN/0Emj5Qc9qv/f7AtppbIBaGlx5UPngtzxFlq9QsmxxgVoi7y/cuJ03V2O0DPuv58yu4bMu2rITZGwyCbQ88yaHPfYjj4b+wePHSST/mnsY9tIRauHHRjSO2vXHRjcRlnBdPvTjp/RoPO+t3EowGuXmxCn9fXe7iw+9u4+1vbuX1wpsJCyv84lPw1FZu8p0mIeB3E3xO8USctxreYsOcDRhE+u20pmQNBmFgT+OeCT2uzhSx6jb4+sfw8DEoH0jOR/lGWH0nGExw7Q+UhQ9QtlqVIRyaUTILzHhxbz55gOq+o7TMuSLbXRmEwSBYubgST+HkFZROkpAJHv/gcebmz+WK+SP/HhY6F3JZ2WU8c/iZ/myHuUYkHuHJg0+yzLOMdaUDfnST0YAQgo21S7kn+g3i+aWAoOqmJ9kybwvPHnl2QoOK3ml6h5ZQC9dVXZdxu8fqYeOcjbxc9zIJmTsl2HQuAptL1dpNRQi49Qn4pxZYnRKNXrUF4n1Q96ep7WMGZrS4H/z19yh+bis9WFl43Vey3Z2s8fThpznqP8r9a+8nz3hhl0ySB9c9SKAvwPff+f5AutwcQUrJo3sfpTHYyEPrHspoMW9fMYf3YtW8fOlz8MW3YNVtfG3d1wjFQvzjnn8klkgvmHKxxBNxHv/gcUrtpWytyJyJE2BH9Q4ag428cvqVcR9TJ8cwDJloseAKVbJx38+z058UZqy4n9r7MiuPPMYB60bO7HiJeeUzqAzYKJFS8uxHz/LjAz9me9V2PrkgPZnZcNR6a7n/kvt5tf5VvvPn7+RMfpRIPMK/vvev/ObEb7in9h42zd2Usd3aCjflHhs/3HkMX0C9Ile7q/nWxm+xp3EP9//x/nGNKUgp+dH+H3Go7RAPrnswrR5tKtuqtrHcs5xH9j7Cue7cKcOmMwmYzLD5a3DyNXj/2ax2RYzGKhNCXAf8GDACT0kpHxmy3QL8L7AOaAdul1LWX+g7169fL/fvn5zpduf+8gccL9xNp3DieGAPbvfEZwHMZaSU7Gncw1OHnuJAywGurriaRy9/9IICNNz3PHHwCZ44+ATzCubxNzV/ww0LbsBlnfrEYv6wn531O3nuyHOc6z7H3TV389D6zFZ7kgNnO7jrqffIMxn4279awCdXlbGwKJ/nTzzPI3sfwWK0cMviW9hetZ0VRSswGUZOkhpLxNjn28fPDv2Mfb59fG755/jGxm+MuF99oJ47f38nZoOZ71z2HbZWbM3ZtMo64yTWB7/8jHLN1NwCG+5VxWBMlpH3HQVCiPellOtHbDeSuAshjMAJ4FqgAdgHfFZKeSSlzZeBVVLKLwoh7gBulVKmp0VMYdziHg3R191GsKOFYEcz3Y3HiTUfw936HhXR0zRSTM8dv2XJstVjP8Y0ZVf9Lh5+82FK7CV8afWX2FG9Y1xCss+3j8fef4xDbYcQCGq9tSz3LmeRaxFz8+dSZCui2F5MobkQq9GKceir6ihIyAQ90R46wh20h9vxh/w0BBs41XmKY/5jHPMfQyJZ4V3BA2sfGNZiH8rHzd384OWj7D6hogY9+WZq5zrwuDo5m3iBut63SRDHarSywLmA8sJyXBYXDosDk8FEPBEnlojRGmrF1+PjSPsRemO9uCwu/n7t31/U7/ZU5ym+/ubXOdl5korCCrZWbGV96Xoum3vZoJw0OjOAWB/s/iG8+4Qq2WeyqlKExUtVSceFV8KCy0f6loxMpLhvAv5FSrld+/wtACnlv6e02am1eUcIYQJ8QLG8wJePVdwPvfa/LP7z17GRPm84KK2cNC2mbf42Lrn5frye2WWxJ4nGo+w6s4ttldtG7WMfDUfaj/Bmw5vsbdrLiY4T/QUqhmIxWrCZbFiMFozC2C9+AtFvaUcTUSLxCH3xPiLxCNFENON3eaweql3VbJizgS3zt1DjHdtUxvOdIf54rIUPGzr56HwXDR0hAqEoGHr58vUJpKWe+q56GrobCEQCdPV1EZdxDMKASZgoshVRYi9hmWcZ6+es58ryK8ckyLFEjFdOv8ILJ1/gQMsBookob93+VlbehnSmgEgQTr+pKla1HVd1ZwONsOVh2Dq2wMSJFPdPA9dJKT+vfb4LuFRK+dWUNoe1Ng3a51Nam2GdmmMV9+MH3+H8m0+TsHoQdg+GgiLMhV7c85ZQUVVNvnXixExneKSUtIfbae5ppi3URmuolWBfkFAsRCgWojfWSyQeISETSCnp/6ctmw1mLEYLZqP6aTFasOfZ8Vg9eKwe3FY3ZflluK2Zc+1PBMFIjMaOEGUuK44hfzfJ+2IyXSfhWJjjHcdZXTz73i5nNYk4xKOQZx25bQZGK+5TWolJCHEfcB9ARUXFCK0zs3T1JpauHt0ruc7kIYSgyFZEka1o5MY5SoHFxNI5hRm3TYU/3Gqy6sI+GzEY02fZTMZhRtGmEUjNRztfW5exjeaWcaIGVgchpfyplHK9lHJ9cfHEV8PR0dHR0VGMRtz3AdVCiAVCCDNwBzA0zO9F4G5t+dPAHy/kb9fR0dHRmVxGdMtIKWNCiK8CO1FTIZ+WUn4khPg+sF9K+SLwc+A5IcRJwI96AOjo6OjoZIlR+dyllL8Hfj9k3T+nLIeB2ya2azo6Ojo6Y2XGRqjq6OjozGZ0cdfR0dGZgejirqOjozMD0cVdR0dHZwYyqsRhk3JgIVqB6VgWvgjI7RJFF49+TtMD/ZymB5N9TpVSyhEDhbIm7tMVIcT+0YT+Tif0c5oe6Oc0PciVc9LdMjo6OjozEF3cdXR0dGYgurhfPD/NdgcmAf2cpgf6OU0PcuKcdJ+7jo6OzgxEt9x1dHR0ZiC6uI8SIcRtQoiPhBAJIcT6Idu+JYQ4KYQ4LoTYnq0+jgchxL8IIRqFEB9o/2/Idp/GihDiOu1anBRCfDPb/ZkIhBD1QohD2rWZnOLDk4wQ4mkhRItW3Ce5ziOEeE0I8bH2c/Kqs0wCw5xTTtxLuriPnsPAp4DdqSuFEDWoLJi1wHXA41rd2enIY1LKNdr/34/cPPfQfvf/A1wP1ACf1a7RTOAq7dpkfZrdGHkGdY+k8k3gdSllNfC69nk68Qzp5wQ5cC/p4j5KpJRHpZTHM2y6GfiVlDIipTwNnAQ2Tm3vdFLYCJyUUtZJKfuAX6GukU6WkVLuRqUET+Vm4Flt+Vnglint1DgZ5pxyAl3cx8884FzK5wZt3XTkq0KID7VXzWn1epzCTLoeqUhglxDifa1c5UyhVErZpC37gNJsdmYCyfq9pIt7CkKIPwghDmf4PyMsvxHO7wlgEbAGaAL+M6ud1RnKZinlWpS76StCiMuz3aGJRqveNhOm7+XEvTSlBbJzHSnlNWPYbTQ1ZnOC0Z6fEOJnwEuT3J3JYtpcj4tBStmo/WwRQvwO5X7afeG9pgXNQogyKWWTEKIMaMl2h8aLlLI5uZzNe0m33MfPi8AdQgiLEGIBUA3szXKfLhrtxkpyK2oAeToympq/0wohRL4QojC5DGxj+l6foaTWX74beCGLfZkQcuVe0i33USKEuBX4b6AYeFkI8YGUcrtWT/Y3wBEgBnxFShnPZl/HyKNCiDWo1+J64O+y252xMVzN3yx3a7yUAr8TQoC6Z38ppXw1u126eIQQ/wdcCRQJIRqA7wKPAL8RQtyLyhL7mez18OIZ5pyuzIV7SY9Q1dHR0ZmB6G4ZHR0dnRmILu46Ojo6MxBd3HV0dHRmILq46+jo6MxAdHHX0dHRmYHo4q6jo6MzA9HFXUdHR2cGoou7jo6Ozgzk/wFnGx50zhEOiAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dos.plot(plt)\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }