{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Bilayer Sonophore model: computation of balance quasi-static deflection" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Imports" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from PySONIC.core import BilayerSonophore" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Functions" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "def plotZeq(bls, ng_range, Q_range, fs=15):\n", " fig, ax = plt.subplots(figsize=(15, 4))\n", " ax.set_xlabel('$Q_m\\ (nC/cm^2)$', fontsize=fs)\n", " ax.set_ylabel('$Z_{eq}\\ (nm)$', fontsize=fs)\n", " for ng in ng_range:\n", " ZeqQS = np.array([bls.balancedefQS(ng, Q) for Q in Q_range])\n", " ax.plot(Q_range * 1e5, ZeqQS * 1e9, label=f'ng = {(ng * 1e22):.2f}e-22 mole')\n", " for item in ax.get_xticklabels() + ax.get_yticklabels():\n", " item.set_fontsize(fs)\n", " for key in ['top', 'right']:\n", " ax.spines[key].set_visible(False)\n", " ax.legend(fontsize=fs, loc='center right', bbox_to_anchor=(1.8, 0.5), frameon=False)\n", " fig.tight_layout()\n", " return fig" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Parameters" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "a = 32e-9 # in-plane radius (m)\n", "Cm0 = 1e-2\n", "Qm0 = -71.9e-5\n", "bls = BilayerSonophore(a, Cm0, Qm0)\n", "charges = np.linspace(-80, 40, 200) * 1e-5\n", "gas = np.linspace(0.5 * bls.ng0, 2.0 * bls.ng0, 5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Balance deflections" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAqYAAAEYCAYAAACOUzQHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xd81fX1+PHXSXKzBwkr7LBBpogCKiKgiKtaa+totcu2ztrW7vrt+Lb9Wdva2tbaYce3Ym3VqqgVQVCWDJW9kZUEGSEJ2Tv3nt8f70/wGjLI4t4k5/nwPm7ymeeTAJ573ktUFWOMMcYYY0ItItQBGGOMMcYYA5aYGmOMMcaYMGGJqTHGGGOMCQuWmBpjjDHGmLBgiakxxhhjjAkLlpgaY4wxxpiwYImpMcYYY4wJC5aYGmOMMcaYsGCJqTHGGGOMCQtRoQ7gbJs/f74uXrw41GEYY0xXI6EOwBjT+XW7imleXl6oQzDGGGOMMQ3odompMcYYY4wJT5aYGmOMMcaYsGCJqTHGGGOMCQuWmBpjjDHGmLBgiakxxhhjjAkLlpgaY4wxxpiwYImpMcYYY4wJC91ugn1jTMvVVPkpLaiktKCKitJqKktrqSytprKs7r2Gmio/tTUBaqsD1FbXfe0nEFAABMH7z30fKURFR+KLjsAXE+m+jnGv2AQfcUk+YhOjiUvyEZcUTVyij4QeMcQnRSMRNpe7McZ0RZaYGmPQgFJyspKCnHIKj5dTdKKckoIqSgsqKTlZSVVZbYPnxcRHEZvgIzbR5xLKxGiioiOI8kUQFR1JlC+CiMgIQFEF1LsfoH6lpsZPbZXfJbXVfmqqAlSU1HAiq4SKkmoCfj3tnpFRESSmxZCUFktyz1iSesaS3CuO1PQEeqTH44uO7LCfkzHGmI5liakx3YiqUlZYRe7hUnKzSyg4VkbB8XIKT5TjrwmcOi46LoqktFgS02JIH5pyKhFMTI1xVcxEHzHxUV7S2XGxVlf6qSipprK0hvLiasoKqyjJd8lycX4lh7bnU1Fc/cFJAklpsaSmJ5DaL560fgn0HpREWv8EIqOs55IxxoQ7S0yN6cJKC6o4frCI3OwScg+XkHe4hIqSGrdT8CqN8Qwam+oqjn3jSU2PJzbRh0hom8tFhJi4KGLioqBP48fVVvspyqug4Fg5BcfLKDhWxsnj5Rx5r+BUsh0RJfQakEjvwUn0HpxEnyHJpA1IILIDE2tjjDEtJ6qnN5V1ZVOnTtUNGzaEOgxj2l0goJw8Wsqx/UUcO1DE8QNFlJysBCAiUkjr76qHvQa55KzngASiY7vuZ9NAQCnOrSD3cAm5WSWcyC4hN7uE6grXLSEqJpL0ocn0G55Cv5E9SB+agi/GugG0gXX8Nca0mSWmxnRSqkrB8XLe33OSw7sLOPpeAdWVfgDiU6JdwjW8B+nDU+g1IJFIn1UHVZWi3Apys0o4drCIY/sLyXu/FBQkQug9KJH+o1IZPDaNfiNTiPJZotoClpgaY9rMElNjOpHy4moO7z7J4d0neX/3ScqKXP/K5N5xDByTSv8RPeg3PIWknrEhb4rvLKoqajl+wCWpR/cXknOomIBfifRFMGBkDwadk8agsWmk9U+wn2nT7IdjjGkzS0yNCWOqysmjZWRuz+PQ1jxyMotBITbBx8AxqQwck8qgsWkk94oLdahdRnVlLUf3FboPALtOUnC8HICElGgyJvZi6KTeDBydahXo01liaoxpM0tMjQkzAX+Ao/sKObQ1j8zteRTnuX6ifYYkkTGxF0PG96T3oCSby/MsKTlZyeHdJ8nekU/WrpPUVvnxxUQyeFwaQyf1Zsj4nsQm+EIdZjiwP5DGmDazxNSYMBAIKMf2FbJ/4wkObD5BRUkNkb4IBo5JZejEXmRM6EVCj5hQh9nt1db4eX9PAYe25ZG5NY/y4mokQhg0JpURU/sybHIvYuK7bZJqiakxps3CLjEVkT8Bkap6RxPHTAV+A5wLHAF+rKpPnsn1LTE14UIDyrEDRS4Z3XSC8uJqonwRDJnQixHn9WHI+J42SjyMaUDJySrm0JZc9m88QXFeJRFRwpBxPRkxtQ8ZE3p16VkPGmCJqTGmzcLmX01xowp+BHwR+GsTx/UGlgBPA58HLgf+KiLHVfX1sxGrMW1RlFvOnvXH2bv+OCX5lUT6IhgyvicjznPJjCWjnYNECOlDU0gfmsL064dzIrOEfRty2L8hh0Nb84jyRTB0cm/GzEhn4Jg0IqzrhTHGNCssKqYiMgyXjI4HyoGljVVMReQ7wBeAEaoa8Lb9HRigqvOau5dVTE0oVFXUsn9DDnvfPs6x/UUgMHB0KmOmpzN0cu/uVlnr0lwlvJD33j3B/g05VJXXktAjhtHT0xkzPZ3U9IRQh9hRLPM2xrRZuPzfcAZwELgF+Hczx84EVtUlpZ4VwOMiElFvuzEho+qa6neuPsKBTbn4awKkpscz/fphjLognaS02FCHaDqARAj9R6bSf2QqF398BJnb8tmz7hibl2SxaXEW6cOSGTOjHyPP72sfSIwxpp6w+FdRVf8J/BM4k3kCBwKb6207CsQDaUBee8dnTEtUltWwd/1xdr51lIJjZUTHRjJ2Rj/GzOhHn4wkmwuzG4nyRTLivD6MOK8PZUVV7H37OHvWHmPFP/ey5vn9jJ6WzvhLBtBzQGKoQzXGmLAQFolpC8UDlfW2VXnvDZagROSLuL6rDB48uOMiM92WqpJzqJgdq46wf+MJ/DUB+mQkM/u2MYyc2tf6jRoSUmKYMm8I514+mOMHi9m56gi71xxjx8oj9BuewrhLBjB8Sm9bbcoY062FRR/TYCKyAtjfRB/T7cBLqvpg0LbLgdeBNFUtaOr61sfUtCd/bYADm06w9c33OZFZjC82ktEXpHPOzP70HpQU6vBMmKssrWHP+mPsWHWEohMVxCb4GDezP+NnDSQxtdNND2ZNAcaYNuuMFdPDQL962/oDpUDR2Q/HdEcVJdXsXH2U7Svfp7yomh5947nk5lGMnp5u/QbNGYtN9DH5ssFMmjOI998rYMeKI2xcksXm17MZMbUPk+YOos+Q5FCHaYwxZ01n/D/oW8BnRUT0g3LvbGCNDXwyHS3/SClb3zzMe2/n4K8NMOicNObcNojB56TZSkym1dwk/WkMGpNGUW4F25e/z661R3nvnRz6jUhh0txBDJ3U26acMsZ0eWHflC8i0bhBTSdVtVpE+gJ7gWeAR4HLgEeA+ar6ZnPXt6Z801KqypG9BWxaksXh3QVE+SIYPT2dibMHkda/y079Y0KsuqKWXWuOsm35+5TkV5LUM5ZJcwYx9qJ+4VqVt6zZGNNmnSExvRRYDsxW1RXetunAb4GJQBbwA1VtbpopwBJTc+Y0oBzamsfGJVmcyCwmPjmaiXMGMm7mAFsb3Zw1gYByaGsuW984zLH9RcQkRDFx9iAmzh4Ybn8OLTE1xrRZ2CWmHc0SU9Mcf22A9945zqYl2RTmlJPcO44p8wYzenq6jZg2IXX8YBEbF2eRuS2PqJhIxs/sz+TLBpPQIywGSlliaoxpM0tMjfHUVPnZ9dZRtizLprSgil6DEplyxRCGn9ubiMiIUIdnzCn5R0rZtCSLfe/mIJHCmBn9mDJvMCm940MZliWmxpg2s8TUdHs1VX52rDzC5qVZVJTU0H9kD6bMH+IGNNlk+CaMFeVWsHlpNrvXHkX9yoipfTn/6oxQLXtqf1mMMW1mianptmqqvYT0dZeQDhqbyvlXD6XfiB6hDs2YFikrqmLLssPsWHUEf7WfURekM/WqDHr0PasVVEtMjTFtZomp6XZqqv3sXHWETa9nU1FczcAxqVxwjSWkpvMrL65m89Jsdqx4H79fGT2tL1OvGkpK77izcXtLTI0xbWaJqek2amv87Fx1lE1Lsij3EtLzrxlKf0tITRdTVlTF5tez2bHqCOpXxsxI57wrM0ju1aEJqiWmxpg2s8TUdHkBf4A964/z7n8PUVpQxYDRqVxwTQb9R6aGOjRjOlRZURWbFmexc/VRVJWxF/Xn/KszSEjpkFH8lpgaY9rMElPTZakqBzfnsv6lgxTmlNN3aDLTrx/OwNGWkJrupbSgio2LM9m1+igRUcKkOYM494ohxMS160T9lpgaY9rMElPTJR3ec5L1Lx7gRFYJqenxTL9+OEMn9bJR9qZbK8ot5+2XD7Hv3RxiEqI4b34GEy4d0F7z89pfLmNMm1liarqUE1nFrF94gMO7C0hMjeGCa4cxenq6rTFuTJDc7BLWLzxA9q6T3t+ToYye3q+tf0/sL5kxps0sMTVdQnF+BesXHmTfuznEJviYelUG4y7pbys1GdOE9/cWsO7FA5zILCa1XwIXfnQ4GRN7tfZylpgaY9qsXTsYGXO2VVfUsnFJFluXHQaB864cwpR5Q4hu375zxnRJA0encuO3zjvVF/vw7pNtSUyNMabN7P/eplMK+APsWnOMd145SEVJDaOm9WX6dcNJSosNdWjGdCoiwvApfRg6qRe1NYFQh2OM6eYsMTWdTtbOfNb8Zz8Fx8roNyKFa+4dSZ8hyaEOy5hOLSIygujIiFCHYYzp5iwxNZ1G/pFS1j6/n+xdJ0nuHcf8L41n2OTeNtLeGGOM6SIsMTVhr7K0hvUvH2TX6iNEx0Vx8cdHMn7WACKjrLpjjDHGdCWWmJqwFfAH2Ln6KG+/fJDqSj8TLh3I+VcPJTbRF+rQjDHGGNMBLDE1YenovgJW/Xsf+UdKGTA6lZk3jaRn/8RQh2WMMcaYDmSJqQkrpQWVrH1+P/s2nCAxLYb5XxzPsHOtH6kxxhjTHVhiasKCvybAljey2bAoEw3A1KszmHLFEHzRNkG+McYY011YYmpCLnNbHquf20dxbgXDJvfmohtHkNwrLtRhGWOMMeYss8TUhExxfgWrn9lH5rY8UtPj+ciXJzPonLRQh2WMMcaYEAmLxFREIoGfAJ8BkoDFwD2qmtPI8XOAnwHjgOPAn4BfqKqelYBNm/hrA2x94zDvvnoIgBk3DGfS3EFE2uTexhhjTLcWFokp8EPg08DtQD7wOPA8cHH9A0VkBPBfXGJ6MzAF+AdQBvz+7IRrWuvovgJWPP0eBcfKGDqpFzNvGmXLiBpjjDEGCIPEVESigfuBL6vqUm/bzcAhEblQVdfWO2U+UKGq/+t9f1BEPgFcgSWmYauipJq1z+9nz/rjJKXFcvXdE8mY2CvUYRljjDEmjIQ8MQUm45rvV9RtUNVMEckEZgL1E9NcIE1EbgGeAc4BLsFVWU2Y0YCya81R1r14gJoqP1PmD2HqVRk22t4YY4wxpwmHxHSg936k3vajwKAGjn8e+CvwT2ABEAk8i+ujasJI7uESVj69l5xDxfQf2YNZt44mrV9CqMMyxhhjTJgKh8Q0Hgioak297VVAQ50PewBDgJ/jKqYTgEeBH3iv04jIF4EvAgwePLh9ojaNqq6s5Z2XD7Ft+WFiE31c9pmxjJqWbpPkG2OMMaZJ4ZCYVgARIhKlqrVB22NwA5rqexjwq+q3ve83i0gU8EcR+a2q5tc/QVX/DPwZYOrUqTZyvwNlbs9j5dN7KS2sYtzMAUy/bhixCba2vTHGGGOaFw7z8xz23vvV296f05v3AaYDG+ptexvwAVYODZGyoiqWPLGDV3+/jei4KG74+nlceutoS0qNMcaElIjEi8ifRSRPRApF5AkRSWzi+EwR0UZeg71jRorIIu9674vIj7wiWVtjvUlEtohImYjsF5Fve1Nq1u2fIiLLvPseFZG/iMhZnQDc+zl8qqOuHw4V061ACTALeApARDKADGBVA8e/D0yst208EAAOdFCMphGqyu61x1j7/H5qqv1M+8hQzp03hMiocPjMY4wxxvAn4DzgGlwR62/etk82cvz5uPErdRKA5cAqVc0WkVRgNbALmA0k4lplBwKfb22QInIlbvzMV4DXgHOBJ7yYfywi/YFluLE29wFpwB9w42wua+19w03IE1NVrRKRx4FfikgecAI3wn6lqq73ppNKA06qajXwG+C/IvIg8DRuVP6vgMdVtTg0T9E9FeaUs/ypPRzdV0j/kT249JOjSU23wU3GGGPCg4gMAG4F5qrqem/bHcByEfmmqp7WMququfWu8QegFm+sCm7e9XjgRlU9GXTNt0Tkx6qa2cpw7wSeV9XHvO8PiMhY4LPAj4GbgErgTlX1e/e9B1glIoNVNbuV9w0r4VLWehD3KeEp3KeSLOBGb9+FwDHvHVVdBNwAXA9sww18+hPwtbMbcvflrw2wYVEm//7xO+QfKWX2p8Zw/VfPtaTUGGPMKV6T7+dEZJWIVIjIHm8wcvAx3xCRbBEpF5EXReQ3IrKiketd2kQTe2YjYVyEa1FdE7RtDeCngUV8GrjnJFxCeq+qlnubRwI765JSz2bv/RLvvBgReUREjolIsYisFJHpzdzuJ8CP6m0LAKne1y8DN9UlpUH7CTqmfvyZInKPiPzX+x0cFJFrReQGEdknIqVel4ReQeeM97YViMhJEXkyeH8D9/iCiOz1rr9dRD7dzHM2KeQVUwBv0NMD3qv+vhWA1Nu2EFh4VoIzH3L8YBHLn9rDyaNljDivDxd/YiQJKTGhDssYY7q8jG+/+ihu7u9Q2JL5s6u/0orzHgbuwRWSHgD+ICJLVDVLRO4Dvg/cBbyLqxjeT8Pd+MDNa15/PEodfyPbBwIngmf+UdVaETlBw1NS1vdD4C2vKFbnKHCtiESoal1imOG99/HenwSGAZ8AcoBbcFXaSar6XkM3UtV3g78XkWTcz2axt/8Ap3dZ/BZuPM6OJp7hZ7jfwf3AI7hC4C4vpiRc14BvAN/yulKuAV7BzSWfCjwGLBWRqfWSYkTkLlwyfTcuOZ8BPCYiqOo/moipUWGRmJrwV11Ry/qFB9i+6giJPWK46u6JDLWVm4wxxjTtb6r6LLjqKHAHcAGuZfQB4Neq+pR37FdFpNEqpted73gL7x+Pa/6ur7EpKU8RkaHAR4Cr6u16Fvgf4GER+QEuufstrrk/WtzS6Z8AxqvqTu+cH3nP9gDwpeaCFpF4XAEuDvh2I8f8DNdv9vr6CWM9L6vqk945TwDXAd9R1Q3etqW4sTrgEsxC4LN1yby41Th34VbefLXetb8H/EhV/+N9f0BEhgDfxS0X32KWmJpmHdrmpoAqK6pi4qUDmXbdMKJj7Y+OMcacTa2sWIbaqeqgqhaKm886WkR64uYkX1fv+DU0UhUWkZm4QUENyVLVcQ1sr8BNP1lfY1NSBvskbuag14M3quo+Efk4H3QjLMVVficCRbhBSwBvy4fn746pi0VESutd89QsAV6z+cu4MTSXq2pW8LHeKP3HcAnuXar6cjPPsT/o67pnDq68VgA9va/HA+/UqzDv9sYAjScoMRWR3sAA3Bihh4OuFwVEiUi092GiRSy7MI2qKK1m9TP72PduDmn9E5j/pfGkD00JdVjGGGM6j6oGtgmuuggtG+uygca7MtRfpKfOYaCPiEQGDRiKwjW5NzQlZbDrgGdU9bT5z1X1FeAVEekH5OOqr4/iEr66RHgGLukLVvfzaCz5zsAlwknAJaq6rd7+WFzFdj7wKVV9uplngIZ/NoEGtkHD1WVwsxTUv05d0nkfQcvKB6ltYFuzLDE1p1FV9m88wepn3qOqvJbzrxnKefNtCihjjDHtQ1WLRCQbmMaHm4en0XAyi6pW8OHq35lYg8t1ZgBvedsuxiXEaxo7SUQScJXPBxvYdzGuX+U8VT3mbbsZV41cywf9TPuq6rKg834P7AYeU9XTnkNE+uAGgPuBC1X1UL39EcBzwBzgWlVd0tzDt8JO4HYR8QU15Z+D62u6K/hA73d4BMgIfh4RuROYrKp3tiYAS0zNh5QVVrHyX3s5tDWPPkOSuO4rY+k5oNF5iI0xxpjW+jnwkIjswVVDP4NbRGdFe91AVY+IyLPAX0Xkc7hq7RPAgrqpokQkDkgBcoP6ak7EVQm3N3DZPcAUXB/T3+Oqn78D/p83bWWxiDwD/Nmbzuk94HO4wV3zmgj390AvXOJZISLpHzyG5uAGQl2D66e7NWg/QH4DS7u3xmO4CujfReQhXEL6O9yc8280cPxPgF95HzLewH2w+BXud9sqlpgawFVJ96w7xlvP7cdfG+DCG0Ywae5AIiKtSmqMMaZDPI5LxB4BeuAqpwtxA5ba0x245GoRrnn5P7hJ7OvcBPwdGApketvqRv83tMx5nohc68V9N26U/g9U9dF69/yZd90UXKX0BlVtKLmrS45vwFVy36m324/L1+oWBPhLA5eYyQcV4VZT1RwRuRz4Be7DQhmuv+s3G0p8VfWPIhKDG9X/O1z3iB/jnr1VpIGuE13a1KlTdcOG+iuadm/F+RWs+OdeDu86Sb8RKcy5bSw9+rb3vwvGmC5Omj/EmA+IyHxge/Ak9yKyGDiiqq1eQcl0blYx7cY0oOxYdYR1Lx5AgUtuHsX4SwYgEfb/F2OMMR3u08BQr7k7H9dMfRlNN3ebLs4S026qMKecNxfs5tj+Igadk8alnxxNcs+4UIdljDGm+7gXN5L9Ndwo9D3ALar6ZkijMiFliWk3E/AH2PLGYd555RBRvgjm3D6GMTP6UW+uNWOMMaZDqWo+cFuo4zDhxRLTbiT/SClvPrmbE1klDJ3Ui1m3jCahhy0naowxxpjwYIlpN+CvDbBxcRYbX8skJj6KeXeMY8R5faxKaowxxpiwYolpF3ciq5g3n9xN/pEyRp7fl5k3jSQuMTrUYRljjDHGnMYS0y6qtsbPu69msvn1bOKTfFx190SGTuwV6rCMMcYYYxrV4sRURCYAFwDpuLVhT+JWNVirqgXtG55pjeOHinjzH7spOF7O2Av7cdGNI4iJ94U6LGOMMcaYJp1RYioiw3BLYX0S6AsEgELcerY9cKs0BERkJW5FgmdUNdAhEZtG1Vb7eeeVQ2xZlk1CjxiuuW8SQ8b1DHVYxhhjjDFnpNnEVET+gktI3wL+F1gL7AxaTxYR6QWcD1yBWx/1hyLyeVVt8/JY5swcO1DEm0/upjCnnHNm9ueiG0YQHWc9NYwxxhjTeZzJQuiVwBhVvVxV/6iq24KTUnDrxqrqa6r6FWAI8H1gQAfEa+qpqfbz1nP7eOGXG/HXBPjI/ZOZ/ckxlpQaY4wxQcR5TUQePIPjvi0ih0SkQkQ2isjVQfv/T0S0kdftbYzxXhHZIyJlIrJLRO6ot/8yEVknIiUikiUivxCRs7Y6johkeM95cUfdo9nsRVXvbckFvSb8Z1odkTljR/cV8OaTeyjKrWD8rAHM+OhwomMtITXGGGOCiUg08AdgPrCmmcO/Bnwbt2TqduAW4CURuUBVNwH3e/uDPQGMAF5sQ4x3AT8D7gTWAbOBx0WkSlUXiMgk4FVcy/TtwDDvvmnA51t733BjWUwnVFPlZ93CA2xf/j7JvWK57qvnMnB0aqjDMl2c1tbiLynBX1joXkVFaEUFgfIKAhUVaKX3dWUlWl0NAT/qD3zw7q91F4qMQiIjEV+U+zrKvSLi44iIj0fi44k49UogMjmJyNRUInv0ICI2NrQ/BGNMpyMiU4C/Aim48THNiQceUNWXvO9/KiJfB2YBm1S1CCgKuv5NuIR3iqqWtCHUO4Hfq+pT3vcHRGQG8FlgAfA5YKuq/o+3f5+IfA/4i4jcqao1bbh32GhRYioiMbgf0GjcaPwdwDZVPdABsZkGvL+3gOULdlOcV8nE2QOZfv1wfDGRoQ7LdGKBsjJqTpygNucEtbknqM3Jcd+fyKU2J4favDz8hYUESs7s31uJjUWio5GICIiMdO9RUe4dUL8f9ddCrR+trXXf19RATfP/pkps7KkkNSotjai+fYnq2wdf375E9elDVJ+++NL7Etmzpy0gYUyIiYjiKnmfwY1DyQJ+pap/DjrmG8B9QC9gCZANTFLVSxu43qXA8kZul6WqGY3smwu8gRsns625uFX1x0H3jMUlhPHAigZiigN+AfxaVbcHbU8FHgGuAwRYD3xVVfc2cesv454/WACoqzw9gUtQ6++PBhJoIOn2fge3AvcAU4GD3vOcC3wXSAYWAZ9R1SrvnIuBnwJTgHJcK/i3VbW8getH4KrHX8L9DncBP1DVRU08Z5NaWjF9Grgel5AmABkuLikDduIy+TtbG4xpXHVlLeteOMCOVUdI6R3HRx+YQv+RPUIdlukEVBV/QQE12dlUZ2dTnX2Y6uwsarIPU52djf/kydPOiUhIcIle377ETZrkksGUFPfq4b0nJxORkIDExRMRF0tEXJxLSiPOpOt6A3FWVxOoqCBQXv7Bq6wMf1HxB1XagoJT77X5+VTt20dtXh4EPjwJiMTHEz1wIL7Bg4geNJjowYPwDRzk3vv3R3w2fZrphH6Y8igwOUR338IPi77SivMexiVF24AHgD+IyBJVzRKR+3BjUu4C3sVVDO8HVjVyrbVAv0b2+RvZjqr+ou7rlnxgFZEbgP/gEsv/UdXNDRz2RVwl9qdB5wku2SvCDQovxyXfb4nIGFXNbyTOlfXuPxjXjeB33v4d9fb7cN0O1qtqU5XgXwNfwE3t+Q8vtreBK3GFxqdxP/M/iMg04E3gt7hkcyjwR+/92gau/RBwg/dzOICrHL8gIvNVdUUTMTWqpYnpPOA+VX0cTn1SmABMDHq1mIhEAj/BfapKAhYD96hqTiPHDwQexf3CK3B/cL7eUDbfFRzefZLlC/ZQUlDJpLmDmHbdMHzRViU1p6stKKDqvX1UvfceVfv2udf+/R+udooQ1S+d6EGDSZo7B9+gwfjS+xLVp67q2IfIxISzHrtERxMZHU1kSkqLztPaWmrz86k94VV7jx2n+nC2S7wzMylb/RZaVfXBCT4fMRkZxIwcQczIke41YgS+QYOQSPt7ZUw7+5uqPgunqqN34OZCz8Ilqr8Oarr+alODalS1GjjewfEGextXWZwDPCwiOar6RN1Or1r4ZeBxr3m/zlxchThNVYu9bXeJyFxcAvdQczcWkd64/qTHcf1O6++PBP4PGAc0NxDpr6r6infeAuAx4G5VPQTsEJEtwHjv2AeADar6de/7PSJyJ7BIRMYBZUExJOIwcdvUAAAgAElEQVQ+SHxMVZd4mx/z+sJ+hwYqzGeipYlpNnCo7htVrQDe8V5t8UNcJ+PbgXzgceB5Gvhhe90JlgLHgIuAnrhPAAGgRQO1wl11RS1rnt/PrreO0qNvPDd8/Tz6DW/Z/7RN16SBADXZ2VRs30Hljh1UvreXqn378eflnTomMiWFmFGjSLn2GqIzMvANHkz04MH4BgwgIiYmhNG3L4mKwte3L76+fWHChNP2ayBAbW6uVzE+THXmIar27adi23aKF732wXViYogZPpyY0aOJHT+OuPHjiRkzpkv9rEwn17qKZai9V/eFqhZ6FctoEemJm8VnXb3j19BIVVhEZgKvNbQP15Q/ru3hfkBVjwBHgK0iMgr4Bq45vc6FfDAAKdi5QCRwtF6FNhYY61VCdzUWuzd3/Gu47gOz6iW9iEg88C9cce5GVd3QzKPsD/q6DJcvZQZtqwDq/qEbj0uIg60O2vd20Pax3nnPiUhws5UPaLCweCZampj+DLibxv9gtJg3Uu5+4MuqutTbdjNwSEQuVNW19U65FVfKv7BupSkR+SGuCaDLyN6Zz/Kn9lBWWMW5lw/mgmuHEmVV0m5JVak9epSKHTup3LGdih07qNyx81QVVGJiiBk1isRLLiFm1MhTVcCo3r2tnyUgERGnEtf488//0L5AWRlVBw+6KrNXYS5dtYqiF72BtVFRxIwcSdz4ccSOG0/suHHEjh6FREeH4EmM6ZSqGtgmgDca8oymrayzgca7MrTbwB9vaqi9qhqc0G3HFc+CXQe8o6oH622vxo3DmdbA5UuBPD78HKdiF5FzcTlWAS7POVwvtp64xPEc4GpVfeMMHqn+z0ZVVRs5trKBbXX/I6l/nWrv/QY+nPxCE10rmtOixNSbruACEVmKK0WvbodRYJNxzfcrgu6TKSKZwExcn5JgVwBLg5c/VdW/AX9rYxxhoaq8hrf+s589a4+Rmh7PDd88j/ShViXtTrS2lsq9e6nYuJHyjZso37QRf65XCfX5iB01iuSrriJuwnhiJ0wgZvhwJMom2GiNiIQE4iZMIC6o0qqq1ObkULljx6kPASVLl1H43H8A90EgbuJE4s6bQvx55xE3eTKRSUmhegRjOiVVLRKRbFzyFlyhm0bDyWxdK239BKgj/ALXzzK4FfYCPlzlBJejvNnA+TtxUzhRl9x6Te//BF7wujac9hwiMgZY5u27qn5fVK/75BJgMHCpN3VVe9uJa40ONtN7311v+z5csjpQVRcHxfkDXMX4+60JoKWj8h/AdWIG14eiRkT2AFu917a6qmcLDPTej9TbfhQY1MDxo4A3ReTHwKcABV4AHlTVhjL9TiNzex4rntpDeUkNU+YP4fyrM4jyWZW0qwtUVlKxdRsVmzZSvmEjFVu2EChz3Xh8/fuTMGMGcZMnEzdhAjGjRlnTcgcTEXzp6fjS00m67DLAJas1R466ivXmLZRv3Ej+E38h/49/AhFiRo8mfsoU4s6bQsIFFxDVu3eIn8KYTuHnwENeHrEBN85kOq3sm9haXsKXAuR6Cwj9Cjd/6Abcqpcfxa2AeV3QOQJMwg0Squ8N3Cj8Z0Xkflyz9rdxg4f+t4lQnsRVLG8DfCKS7m2vVdU879zJwEdw3QTSg8490U5LwT8MbBaRX+K6KGQAvwcWqepuEcmoO1BVy0XkV7jfYTHud3gNLiFt9byqLS2zfA94CngQNyp/Iu4XMxHXHD8AlyW3RDwQaKDyWoXrj1FfMu6BXwM+7t3zMaA3rp/qaUTki7gOxwwePLiF4XW8yrIa3np2H3vfPk5a/wSuunsifYYkhzos00HU76dy127K1q2jbN1aKjZucvN+ihAzciQp132EuCnnEX/eFHz9GhuAas4mESF64ACiBw4gef58wHUDqNi2jfKNm6jYtJHChQspePppAGJGjiB+xgwSps8g/oLziUxMDGX4xoSrx3FTDD0C9MBVThfi8oKz6Sbg77iR55mq+hcRicIN4BkC7MUN8AmeAikVl6OcNsJeVVVErgd+CbyE64e5GZivqvWrrgB4fVjr+hrVn1LqAG7y/k/icqz6fUDBFfLeb/5Rm6aqO0TkGtyA9Ptwz/dvXN7XkAdxTfq/BPripqP6kqr+X2tjkMa7GTRwsEgOcIuqNlS6RkR6NDNlQUPnfAw3qt6nqrVB29fgRobdX+/43bg5u0bVLY0qIjcCzwG9GpuGoc7UqVN1w4bm+gmfPQe35LLy6b1Ulroq6dQrM4j0tW66HROeVJWaw4cpW7uOsnXrKF+/Hn+R68seM3o0CTNmED/tAuKnTGnxiHQTPrS2lsrduylbv57ydesp37jRzQYQGUncxIkkzJh+qvrdRaersg7NpkVEZD6w3RtkVLdtMXBEVbvMSkamZVpaMX0K14TfYGLa0qTUU9ext1/Q1wD9Ob15H29bZV1S6qn7BJJBA59ewlFFaTWrn9nHvndz6DkwkWvunUTvwdZPrasIVFdT/u67lK5YSenKldRkuzmTo/r1I3HuXBJmzCBhxnSievUKcaSmvUhU1Af9Vb/wBQJVVVRs3uIq4+vXkffHP5H3+B+ISEwk4aKLSJw1i8RLZtqfAdOdfRoYKiL34P7ffQ1wGW5qStNNtTQxzQK+IiJHgT/WSw5baytQglvq6ykArw9DBg1Psrsa+IKI+IKa/8fjRoBltkM8He7AphOs/NdeqspqueDaoUy5YgiRUVYl7exqck5QusolomVr16Hl5UhMDPHTp5F2++0kXHQh0RkZNlK+m4iIiSFh+jQSpk8DvoK/uJiyt9+mbNUqSleuomSJm/Yvdvx4l6TOuoTY8eNbvUCBMZ3Qvbg5yV/DDYLeQxOtsqZ7aGlTfikf9P0owHUK3oI3+Km1S5OKyM9wnZ4/A5zA9TupVNVLvemk0oCTqlotIn1xo8Zex3UEHohbA/cNVf1cc/cKZVN+eXE1q/79Hgc2naD34CTm3D6WXgOt71lnpapU799P8dKllC57g8pdrnAf1a8fibMuIXHWLBKmTyciLi7EkZpwo6pU7d5N6cqVlK5cRcXWraBKZM+eJM2ZQ9K8y0mYNq2zTUtln7iMMW3W0sRUcJPJTsSt+FS36tNw3FxkZara4vZor5Pxw7iyvo8PVn7KC1obd3bd8lYicg5uia2ZuDnBngK+U7fOa1NCkZiqKvs3nmDVv9+jurKW868eyrnzBhMZaZWRzkZV3fRBr79OydKlVGdmgghxkyeTOHs2ibNmETNqpFVFTYvUFhRQtno1JcuXU7ZyFYHyciISE0mcPZukyy8j8eKLiYg/2+NBWsz+0Btj2qxFiWmjFxGJxTWnj2/LSKyz4WwnpmWFVax4ei+Z2/LoMySJOZ8eS8/+ViXtTNTvp2LTJoqXLqVk6TJqjx2DyEgSpl1A0rx5JM6Zg69Pn1CHabqIQFUVZWvXUrJ0GaVvvIG/qAiJjSXh4otIvvxyEi+9NFwHyVliaoxps3ZJTDuTs5WYqiq71xxjzfP78dcGmHbtMCbNHUiEVUk7BQ0EKN+wgeJXF1GybBn+/HwkOpqEiy8m6fLLSZp9KZE9eoQ6TNPFaW0t5Rs2UrJ0KSXLllGbkwNRUSTMmEHy1VeRdNll4TQVlSWmxpg2a3FiKiLXASNxS2rtBHZ4qzF0CmcjMS3KrWD5U3s4sreA/iN7MPu2MfToE/bNcN2eqlK5cxfFr75K8aJF1ObkIHFxJM2+1FVGZ84kIiEh1GGabkoDASq3b6f49dcpeW0xNUePIjExJM6aRfLVV5M46xIiYhua+vmsscTUGNNmLe1j+mfc5PbHcZPhxuFGwx8EtuEGQP2kA+JsNx2ZmAYCyrY3D/P2SweRSOGij43gnIv6IxH273U4qzp4kOL/vkrxq69SnZUFPh+JM2e6itTs2Z2hb5/pZlSVis1b3IeoxYvx5+cTkZBA0mVzSb76ahJmzAjFXKn2D50xps1ampgWAr9Q1Z963w/ng5WfJgITVHVkRwTaXjoqMc0/WsryBXvIOVTMkAk9ufTW0SSmhrR6YZpQk3OC4v/+l6JX/0vVrt0gQvy0aSRffRXJl19uzfSm09DaWsreftt1O1m6lEBJCZGpqSTNv4Ie111H7KRJZ2swniWmxpg2a2liehj4jKq+0XEhdaz2Tkz9tQE2Lcliw6JMouOimHnTSEZO7WujssNQoKqK0jfeoHDhQsreWgOBALETJ5Jy9VUkzb8SX18bwGQ6t0B1NWWrV1P86quUvLkcrawkOiODlOuvJ+Uj1+Lr378jb2//6Jmw5s0stAhY01zrrojcBnwLt0zpDuBBVV0atD8aeAi3TGg8bo71e1X1UBtjvAm3FOpI4BjwF1xBsG6lyynAz4GpQLn3PN9U1ZNtuW8LY1TgNlV9qiOu39IJ9v8PmA902sS0PeVkFrN8wW7yj5Qx8vy+zPzESOKSOtW8g12eqlK5dSuFCxdSvOg1AsXFRPXrR88vfoGU664jZujQUIdoTLuJiI4mae5ckubOxV9aSsmSJRS9uJDcRx8l9ze/IX7aNFKuv47kefOsi4rpVrxE8g+4HGZNM8feCvwD+B7wPHA58LKIXFk3bSXwJ9wqVbfi5l//rXfMRG3lqHIRuRL4J/AV3KID5wJP4KbR/LGI9AeWeTHdh5vj/Q/As14sXUJLK6ZfAh4AHgMeD17bvrNoj4ppTbWfd145xNZl2cSnxDDr1tEMnWjLCoaTmpwcil56maIXX6T60CEkNpakeZfT46MfJX7aNFtdx3Qr1YcPu78PL71EzeHDSHw8yfPmkXL99cRfcH57/X2wiqkJS16V8a9ACpAKPNJUxVREtuAGdn8qaNsTwAhVnS0iw4ADwNy6VapEZDRu4Z+5qrq/lXG+hFtc6Kagbf8DfFZVh4nIV4FvAIOCKqgzcatkDlHV7NbctxVxhlXF9BFcyfpR4Ecishq36tMW3MCnVv0yOpMjewtY/tQeinIrOGdmfy68YQQxcS39MZqOEKiqomTZMopeXEjZ2rUQCBA39Tz63fF5kq64Ipym1THmrIoeNIje995Dr3vupmLjRgoXLqTktcUULVyIb+BAUm74KD0++lF8/fqFOlTThXgJzOdxqzqej1vW/Feq+uegY76Bq/71ApYA2cAkVb20getdiltwpyFZqprRyL65uJbe/8UN1G7OSFwlMthm4HZvQaB5QG7w0qmquhcYEhRrDPD/cBXVBO/8b6nq+ibu+xOgrN62AC6ZBngZ2FBvOfiA956K+9l9iIhkAr8ArsT9HI4B9+OqsA8D/XCJ7e2qmuedMx7XXWAGoMB/ga/V7W/gHl8Avg4MBvYDv1TVfzTxnE0Ki5WfzqbWVkyrKmpZ98J+dq4+SnLvOGZ/agwDR6c2f6LpcFX791P43HMULXwJf1ERvv79XZ+6668jevDgUIdnTFgKVFRQsuwNCl94nvJ160GE1Ns+Rfp3v9vaS1rFtINN+MeER4HJIbr9lu2f3v6VlpzgJaZ5wD24hPAB4HPAMFXNEpH7cMnbXcC7wJ24pGlVI4lp3RLlDfGrau4ZxJQJ/KWZiuk+4BlVfTBo289x1crewHeBC3HN99/ytq0Bvqqq73vHP4PLl74G5AC3AN/GJd3vNRend41kYBewWlVvaeSYl4EpuIqpv4H9mUBP3O9gDa7AOMe77r1AEq5rwBOq+i0RycAVHF8BfoZLeB/DJcBTVdUfXDEVkbuAHwF345LvGd7x97c2OW1Rqc/rN3HAe70Y9OCnVn5qTRCdwdZl2ex66yiTLx/MBdcOxRcdGeqQurVARQXFi5dQ+OyzVGzeDD4fSZfNJfXjHyd++nRrqjemGRFxcaRcew0p115D9fvvU/TCC0RnZIQ6LNP1/E1Vn4VT1dE7gAtw1dMHgF8HNQl/VUQubuxCqlqNm66yoy0AviYiy4EVwCW4yi9ANJAMjMHF/1WgCjcQ6g0RmQQMBD6BWw1zp3fej7xnewD4UnMBiEg8sBA3Lee3GznmZ8A1wPUNJaVBXlbVJ71zngCuwy3jvsHbtpQP8re7gUJc94Eab//NuER2PvBqvWt/D/iRqv7H+/6AiAzBJe8dk5h6I9OebuqhVbUS2OC9EJERQD9VXd2aoMLRuVcMIWNiL/oMSQ51KN1a5Z49FD77HEWvvEKgpITojAz6fPObpFx/HVFpjX2QNsY0JXrgQHp/+cuhDsM0o6UVyzBxqjqoqoXejDXRItIT1/S9rt7xa2ikKuz1p3ytkftkqeq4tocLuEphH+9ekbjFhH6BSz6LgBpcf9Ub60bhi8iNuGbyq7xzAN6uN0NPjPdCREqDd6jqqb5mItIL12x/DnC5qmYFHysikbiq5JeAu1T15WaeJ7ibZV1XgQNB2ypwVVVwCeo7dUmpF9tuEcnz9p1KTEWkNzAA+KWIPBx0vSggSkSivQ8TLXImFdMHcKPBFgD/UdWtDR3k/SGbD9wMXMoHny66BF90pCWlIRIoK6No0SIKn/sPldu2IdHRJF1xBamf+DhxU6fa1FzGGBO+qhrYJkDd4OmWNG9toPGuDDWNbG8xL5m6V0QeAFJV9biIfBnIUdUyETmC67p4KOicEyKSj5teqi4RnIFL+oLV/TwaS74zcIOokoBLVHVbvf2xuFH484FPqerTZ/BIDf1sAg1sA6hsZHtkA9epSzrvw1WW62vVAPlmE1NVnezNq3Uf8D0vy9+N6zdShVsBaiiu02sB8BRwp6oeaU1AxtSp2rePgn/9i6KXXiZQVkb0iOH0/e53SPnIR2wCfGOM6cRUtUhEsoFpfLh5eBoNJ7N4y593+CBrEfkJUKKqD/NB14HrcQkjuDlLfywiY1V1t3dOOm4A1wFchRWgr6ouC7ru73H502MNDRYXkT64wV1+4ML6c6KKSATwHK6P6LWquqQ9nreenbhBXr6gpvxzcH1NdwUf6P0OjwAZwc8jIncCk1X1ztYEcEZ9TFX1GeAZb6Wny3AdbdNxI81ycCO61gArgsu/xrSU1tRQ8sYbFDz9L8rfeQfx+Ui+6kp63HQzcedOtuqoMcZ0HT8HHhKRPbhq6GeA6TRcfeswIhKHa5rP9botZgKPiMh2YA9uXtHzcYO0wOU8q4GnReRuXPP4o8BeYJGqVnuDn/4sIvfgujN8Dje4a14Tofwel9zOASq8ZBfcEJ8c7/7X4Prpbg3aD5DfTvnXY7hC5N9F5CFcQvo73ICohuaw/wnwK+9Dxhu4Dxa/wv1uW6Wlg5/qBj4Z065qck5Q+OyzFD77LLW5ufj696f3A1+jx8c+Zn1HjTGma3ocl4g9gmt9fRU34Odsr/5wE/B3XOtvpqr+xUv6/oRLzN4F5nhTQqGqKiIfAX7pxezDTXx/e1CfyjtwfVX/jkt6dwM3NLZyppcc34Dr2vBOvd1+XL72Se/7vzRwiZnAWy187tOoao6IXI7rU7sBl3S/jFtd6rTEV1X/6E2N9Q1cAnsE+DHu2VulRdNFdQXtvSSpaT1Vpfyddyl4+mlKli0Dv5+EmTNJvfUWEi+5BIm0mQ+M6USsOcO0iIjMB7YHd/0TkcXAEVXtUuNUzJmzmeHNWecvLaXopZco+Ne/qN5/gIiUFNJuv53UW262eUeNMab7+DQw1Gvuzsc1U19G083dpouzxNScNacGMy18iUB5ObHjxtHvpz8l+aoriYiLC3V4xhhjzq57cX0zX8ONQt8D3BK8opLpfiwxNR1K/X5KV67k5JMLKF+/HomOJvnKK0n95K3ETphgg5mMMaabUtV84LZQx2HCiyWmpkP4S0spev55Tj71T2oOHyYqPZ3eX/saPT5+I1GptpSrMcYYY07X5sRURMYAe7W7jaIyDarOzOTkU/+k6IUXCJSXE3fuufR54GskXXYZEmWfg4wxxhjTuPbIFP4KTBKRbcDmoNf2M12Kylte6ye4OcySgMXAPd68Xc2d+18gUVUvbVX0ps1UlfJ16zj55AJKV66EqCiSr5xP2m23EzdhfPMXMMYYY4yhfRLTY8AfgR24Jba+h5t0f6iI5KnqxDO4xg9xo/Nux43Mexx4Hri4qZNE5EvA1cDK1gZvWi9QUUHRy69Q8NQCqvbtJ7JnT3rddRc9br4JX58+oQ7PGGOMMZ1MeySmw1R1gff1ZhFZCvxUVS8SkQHNnSwi0cD9wJdVdam37WbgkIhcqKprGzlvBPD/gHXt8AymBWqOHaPg6X9R+Oyz+IuKiBk7ln4PPeRG18fEhDo8Y4wxxnRS7VIxFZEr6tZsVdX3RWSq9/WRpk8FXJU1iaAlyFQ1U0QycSsZnJaYek3/TwIPA6OAEW18BtMMVaVi8xZOLniSkteXgipJc+eSdvttxE2daqPrjTHGGNNm7ZGYfhG3Xux9uOWrxgJ5LTh/oPdeP4k9Cgxq5JzvAIpbDuzPLbiXaSGtrqZ48WJOPrmAyh07iEhOJu3Tnyb11luJHthsQdwYY4wx5oxFtPUCqnpEVWfh+oXWAMuBa1twiXgg0MAarFVAbP2DRWQK8ADwaVUNnMkNROSLIrJBRDbk5ua2ILTuqzY/n9zHH2ff3Lkc/ea3CJSVkf6D7zNy+Zv0/eY3LCk1xhjTKYjICBFZKCJ5IpIrIs+JSJPLDIrICRHReq8Hg/b/ooH9+9sh1ntFZI+IlInILhG5o97+y0RknYiUiEiWF8dZW6FGRDK8Z21yDFBbtNv8Paq6CFjUilMrgAgRiVLV2qDtMUBZ8IEiEgs8BTyoqmf8B0BV/4xXWZ06dapNa9WEyt27OfnkAopffRWtriZh5kzS/t9tJFx0ERLR5s8xxhhjzFkjIgnAEmAXMAeX9zwCvCYiU1S1qoFz+gK9gUuAfUG7SoK+Hg/8HjejUJ3gHKY1sd4F/Ay4Ezd+ZjbwuIhUqeoCEZkEvAr8HDdYfBjwBJAGfL4t9w4n7TGP6Q3Ab4BSVR0rIn0AVdUzLU0e9t77BX0N0J/Tm/en4boKPCwiD3vbYnCJbSlwjqpmt+Y5ujP1+yl5800K/vEk5Rs2IHFx9LjxY6R+6lPEDBsW6vCMMcaY1poHDAbOVdViABG5HcjG5RSrGjhnPC7JfLuJaS/HA8+q6vF2jPVO4Peq+pT3/QERmQF8FlgAfA7Yqqr/4+3fJyLfA/4iInc20PLcKbVHCex/gSlAZdA1n27B+Vtxn0Jm1W0QkQwgg9P/wLwDjMQNmKp7vYjr2zoZ1y/VnCF/cTH5f/s7B+ZdwZH7vkzN0aP0+eY3GblyBenf/74lpcYYY1rNa/L9nIisEpEKr4n6i/WO+YaIZItIuYi8KCK/EZEVjVzv0gaaz+temY2E8Q5wVV1S6qnrBtjYMoTjgQONJaUikoIbH7O7kfMRkVQR+ZuI5IvISRFZJCKjGzve82Xc9JvBAkFxPgHc3cD+aCChkThURG4RkbdEpNLrHjBdRO4SkcMiUiQi/xKRmKBzLhaRlV53gRwR+a2IxDdy/QgR+a7XraBMRN4Vkauaec4mtUdTfrGq5taNylbV4yLS80xPVtUqEXkc+KWI5AEncP1VV6rqem86qTTgpKpWAB9qwheRYqCiJU373V3VoUMULHiKwoUL0fJy4qdOpc+3vknSnDm2OpMxxoSp3WPGPoorwoTClrF7dn+lFec9DNwDbMOND/mDiCxR1Sxv0PT3gbuAd3EVw/tpuIoJbpaefo3s8ze00ZsdqH7r67eBcuCtRq41HqgVt4DPVO/8R4OmxqxbOeazIlJXiHsN+K6qFolLiBYBRcAV3r3uA94SkTGqmt9IrB+ak93rB3sL8Dtv/456+33A14D1qlrYyLMA/Br4AvAe8A8vtreBK4HRuGLiKtzvZhrwJvBb4EvAUFyyPJSGxw89BNyAGwh/AJgPvCAi81V1RRMxNao9spBDXr8HhVM/qMQWXuNBwIfrP+rDW/nJ23chbkDVbIKmlDIto6qUrVnLyQVPUrZyFeLzkXz11aTdfhux55wT6vCMMcZ0TX9T1WfBVUeBO4ALgCxcovrroKbrrzY1qMarYLap6dzrx3kvcF9jCSIwDugJ/A9u0aArgb97Y2H+7u0HtyDQdbik7RHgHBGZA8wFzgfSgiq1d4nIXFwC99AZxNkb15/0OK7faf39kcD/ebE0NxDpr6r6infeAuAx4G5VPQTsEJEtfJBsPwBsUNWve9/vEZE7gUUiMo6gsT8ikoj7IPGxuilDgce8nPA7tDJna4/E9Gu4hHKAiDyAy5YbnBS/Md6gpwe8V/19K4BGJ8lU1Tsa22e81ZleepmTCxZQfeAAkb160eu+e0m96SaievUKdXjGGGPOUCsrlqH2Xt0Xqlrota5Gey2rQzh9kZw1NFIVFpGZuMpkQ7JUdVwj++rO/x5usNJDqvpYE4fOBqJVtW6w01YRGYLLd/6Oa1J/QVXrpsbcLiI5wHpc18ZzgUjgqHx4ju9YYKxXCd3VWOwiMsx7znhglqoW1XuOeOBfuGrsjaq6oann5sMtzWW45v/MoG0VuPE64BLUV+udvzpo39tB28d65z0nIsGzJPlwK4C2SrOJqZds/v/27jxOrqrO+/jnV1vvS3oLWTqkkxCyJ0AQBmQRkUUQ0XHEYXd5FGXUYdyYGR51HOdxQx0dcdzHBZBFHWFGJKAIigw+DzOQhJCNbJC9O73vtfyeP86p6urq6i3ppKq7f++87uveOvdW1Tnp7VvnnnvuL1V1R7b9qnpIRC4F3gisAu7BBVWTQ+7uTPfS8sCDJNraKFy2jNlf+Dxll19OIBLJdfWMMcZMD0Ouesd1NiWvYB/PtS7J60myGfbCHxEJ4IYIvg/4hKp+caQ38VfqZ9Z7I+60OqqqDJ2vfaNf1wP9QDPu4qpMnf656e1I1V1ETsOF0hbgHFVNvygcH+h/BSwDrlDV347UlszX99S3IZveLGXJdJ35OskxuG8lY5glwwytGIux9Jh+CfhbEblaVbOOx/Dzif6nX0yOZL0708UXUzsR9fsAACAASURBVHXTjRSdfrrdnckYY0xe8GMxk1fGp/fQnUX2MEu260zG6Bu4IQTvVNUfjnSgiISAXcBXVPWrabvWApv8MXcCr1PVMzL2g+sJ7cRdG0Py+hd/6v0eXE/rA9naISJLgN/4fW/MHGogbr7SdbhZBi5U1f8ZteXjtwk4N6PsPL/OvNhrOy6szlXVR9Pq+Slcj/Enj6YC4zmV/5iIXK+qv8i2U0RWAJ9R1bceTUXM0dP+ftrXrXN3Z9q4MXV3pqrrriU8xybCN8YYk5e+CHxORLbgekNvBs5mAq8n8VeIvx/4B+BRETkpbXerqvb6sZKlqnpQVWMi8h/AHSKyAxc0rwZuAK7wz/sF8Nci8kXcHOkLcD2y96jqNhHZjjut/4CIfBh3Wvt23MVDnxmhuj/G9VjeAITT6hrzwwY+g+tpvQo3TCC9LYfHetOhUXwBeN6H7+/iZki6C3hEVTf7WZMAUNVuEfkK7mvYjvsaXokLpEc9r+pYg+l1uIGsD4jIbar6L1mOKcENAjYnSKy5mdb776fl3p8Sa2wk0tDAzE/+byrf/GYCJVlnjjDGGGPyxTeBGtyFQ5W4ntNf4sZWTpTr/fpTfkl3A27o4Uf9vuRpxdtwp9K/jpsFYAvwdlV9DEBVnxGRq4BP4y7U7sBd2f53fr+KyNW426Y/hBuH+Txwmaqmjy1NEZHFuAumALZm7N4BLMJlsSBDx4CCG0Kwd5j/gzFT1RdF5ErcWNwP4i7wug93kXo2d+BO6d8JzAR2Au8brWd6JDL8MAN/gBvQejbwAm6C17cBd6rqJzKOOwt4RlWDR1uZE2Ht2rX63HOjjRPOb71bt9L84x/T/h//6e7O9NrXUnXTjXZ3JmNMLtlYITMuInIZsNFP6ZQsexTYp6pT5k5GZnzGfCrfT9NwjYh8FfiYiMwFbp4qdxrIdxqP0/nUUzT/6Md0/+lPSGEhFW99C1U33EDBwoW5rp4xxhgzXjcBDSJyK65n7krgYtzdmsw0Ne7polT1NhHZhxuHMFNE3pI2pYKZYPHOTtp+8Quaf3I30VdfJTRrFnUf/QiVb3sbwcrKXFfPGGOMOVp/Bfwz7ir0Mtwp879U1SdyWiuTU0c1j6mq3unD6b8BfxCRyye2WqZvxw5a7rmXtl/+kkR3N0WnnUbdR/6GsosvtrszGWOMmfT8Vec35LoeJr8cdcJR1Z+KyEHcver/i5GvNDNjoPE4nU8+SfPdd9P9X8+m7s4047rrKFq5YvQXMMYYY4yZxMYSTBcwdCJZAFT1dyJyPq4b/tsTWbHpJNbSQtvPf07LvT8lun8/oZNOova226j8i7cRqqrKdfWMMcYYY06IUYOpqu4eZf8GEfkzXDhdMkH1mhZ6X3qJ5rvvof1Xv0L7+ih+zWuou/0TlF10kZ2uN8YYY8y0MyHpR1VfEZEzGbjzgRmG9vfT/vjjtNx9Dz3PP48UFVFx9dXMuO5aChcvznX1jDHGGGNyZsK65VS1G/j9RL3eVBM9fJjWBx6k5f77iDc2EZ43j5l/ezsVb3kLwfLyXFfPmGFFE1G6o910RbvojnbTl+gjGo/SH++nP9FPf7yfaCI6aB3XgdskqyqKDtkGCAfChAKhgXUwTFjChINhQuIeFwYLKQoVURwuTq0jgYjdYtcYY6YgO198HKkqPf/937T89D7a162DWIyS88+j6p+up+S1r7XJ8M0Joap0Rjtp7Wulra+N1r7Wwdu9brs92p4KoMkQ2hXtoj/Rn+smDBGUoAupoWKKw24pj5RTUVBBRaSC8oLyoeuCCqoKq5hRMINgIK/vA2KMMdOWBdPjIN7RQdtDD9N6/330bX+ZQFkZVdddy4xrryVy8sm5rp6ZIqKJKI3djTT1NNHY08iRniM09rjHTd1Nqe0jvUeIJWJZX0MQyiJlVBZUUh4ppyRcQmVpJSXhEkrCJRSHiykJlaQeF4WLKAgUEAlGiAQjhAOudzMScI8jgUiqtxMY1KuZ3Bb/T1FiiRjRRHTYdTQepTfeS3esm55oj1vHeuiO+rUv74x20t7fzsGug7T3t9Pe105Ms7c5IAEqCyqpKaqhurCa6qLqQdu1xbWcVHwSM0tmUhQqmuCvmjHmRBORZcBXgHOAPuDnwCdUtW0Mz/1L4B9VddEIx9wFXKGq8yegrn+Fm9+1HtgDfEVVv5e2/2LgH4EVQDPwAPBJVe051vceY/3mA7uA81T16ePxHhZMJ1DPi5toue+ntP/qEbSnh8IVK5j1T5+l/PLLCRRP5K1/zXTQ0d/Bga4DHOg84NZp2/u79tPY3TjotDi40DejcAa1RbXUFNewqHIRNUU1zCicQWVBJZUFlVQUVKS2yyJlU7L3UFXpjnW7nuD+9lTvcHNvM0d6jqQCe3NPM3va93Ck9wh98b4hr1NZUMnM4pmcVHJSakk+ri+rp7aodkr+/xkzVYhIKfAb4EngNcAM4Hu4edjfOspzrwS+D+wf4ZhLgQ/gQuSx1vX9wOeBW3DTcL4O+KaI9KnqT0RkNfAr4IvAjbhZk74LVAFT5hauFkyPUaK7m/ZHHqHlvvvpffFFdzHTlVdQ+fZrbO5RM6J4Is6h7kO80vEKr7T7peMV9nbu5WDnQTqig2+oFgqEmFUyi1klszh71tnMLp3NzOKZ1BXXuZ6+olqqCqsIBezHWkRSvbyzmT3q8cnhDsle54NdB1PLoe5DHOw6yAuNL9DWN7iDJRwIM6d0DnPK5lBfWs/csrluKZ1LfVk9xWH7QGpMjp0MPA38L1XtAhCR7zDC3OsiUgR8DbgZ2AyUDHNcFS64PgXMn4C63gLcpap3+8c7/KxH7wR+ArwLWK+q/9vv3y4ifw98T0RumSq3iLe/YEep7+WXabnvftoeeohERweRRQuZeccdVFz1JruYyaQkNMGBrgPsbtudCqCvdrzqAmjHXqKJgd8jBcEC6svqmVs6l7Uz17oQWuqC6OyS2VQXVRMQG5d8PIi4IQ1lkTLmV8wf9rjuaDeHug9xoOsA+zr3sbdjL692vMrejr1sOLxhyIeJqsIq5pbNZX75fBoqGmgob6ChooH6snrCwfBxbpUxuSUiiuvJuxk4k4FT099JO+ZjwAeBGmAd8AqwWlUvzPJ6FwK/G+bt9mQ7la6qm4C3p73GYtzdph4boep1wELgz4A3AdcPc9y3gIdxParvyajrDODLwJsBAZ4FblPVrSO874dw7U+XwPXygusd/UmW/RFceG7NfEH/NbgWuBU3c9JOXMA9Dfg7oBx4BLhZVfv8c14L/BNwOtAN3A/c7i9yz3z9AHA78D7c1/Al4FOq+sgI7RyRBdNxSPT30/HY47Tedx/dzz2HhMOUXXopM95xDUVnnGFXCU9j0USUV9tfZWfbzoGldSe723fTExsY+lMUKqK+rJ6FFQu5sP5C5pXN4+Tyk6kvq6euuM6CZ54rDhe7gFnRkHV/W1+bC6udLqwmg+uz+5/l4R0Pp44LSpC5ZXNpKG9gfsX81Gs2lDdQWVh5oppjJpm7bnnin4E1OXr7F2791kV/fRTP+wIuFG0APgL8q4isU9U9IvJB4JPA+4H/h+sx/DDDz/DzDDBrmH3xYcpTROQFYDUuIF893HGqugd4vX/Om4Z5retxwW01cFvGPsGFvTbgUly4+yDwtIgs8bdizfa+T2W8zjzgL4F/8ftfzNgfBv4GeFZVh4TSNF8F/hewDfiRr9ufgMuBU4F7cf/n/yoiZwFPAF/Hhc0GXABvwIX0TJ/DDYl4L7ADuAz4hYhcpqpPjlCnYVkwHaO2hx7i0Be+SLy5mXB9PXUf/QgVb32r3ZlpmoklYrzS/grbWraxrWUbu9p2sbNtJ6+0vzLoYptZJbNYULGAM2aewYLKBTSUN3By+cnUFNXYB5gprKKggoqCCpbXLB+yr7O/kz3te9jZtpNdbbvY3b6bXW27eGb/M4NmPqgurOaUGaewqHIRi2csZlHlIhZWLrRhAWay+oGqPgCp3tH34MZ67sEF1a+mnbq+zffWZaWq/cDBY6jLu3A9i18Aficiq7P1Ao5GROpxp/rfrKpdWX6nvx7XQ1ylqu2+7P0i8npcgPvcGN6jFjee9CBu3Gnm/iDwQ2A5MOz/mfd9Vf0P/7yfAN8APqCqu4AXfWBPjj38CPCcqn7UP94iIrcAj4jIcqArrQ6luA8Sf66q63zxN/xY2L/FjesdNwumYxSsqqLo9NOYcc07KDn3HJvqaRo40nMkFUC3tWxje8t2drTuSIWIoASpL6tnQcUCLpp3EQsqFrCgYgENFQ0WIswQpZFSltcsHxJa44k4+7v2s6ttF7vadrG9ZTsvt77Mz7b9jN54L+AuaptTOmdQYD1lxinMK59HOGBDAqaLo+yxzLVtyQ1VbfUhLiIi1bjxn/+VcfwfGaZXWETOw91lMps9qjr0E2EaVf0f/zp/DuzF9ZreO4Y2pNdBcIHwByNclX4aEAT2Z4TWQmCp7wl9abi6i8gCXDuLgQsyZw8QkWLgp7je2Lep6nOjVPvltO0u3On/3WllPUCB316BC8Tp/pC2709p5Uv98x4UkURaeRg4NEqdhmXBdIxKzzuP0vPOy3U1zHEQjUfZ0baDrc1bUwF0W8s2jvQOnG2pLapl8YzFnL30bE6ZcQqLZyymoaKBSDCSw5qbqSAYcB9w6svqOX/u+anyeCLOvs59bG/dngqr21u28/u9v0/dwCAUCLGwYiFLqpawtHopS6qWsKRqCSXhrNdqGJMLQ6e7cGMuk6eYxtPL8xzDD2XIeuGPn95otao+lCxT1QMicgSYM473TpoHXAT8mb+KHlwQC4tIJ+70eD9uKqezsjy/E2hicDtSdReR03ChtAU4R1VfzWhPNS44LsNNUfXbMdQ58/9GVVWzHgm9WcqS6TrzdZKnet7K4PALYxhaMZy8Caa+W/qzuEHSZcCjwK2qmjV1i8g1uK7iU4ADuOkfvqSqR/2fYaa+vngf21u289KRl3jpyEtsbt7M9pbtqYuQCoIFLKpcxPlzz08F0FNmnEJVoQ3ZMCdWMBBkXvk85pXP4/XzXp8q74/3u55VH1i3tmzl6X1P89CO1N9d5pXNSwXVpVVuXV1UnYtmGJOVqraJyCu48JbeQ3cW2cMsfq7OzAA0mtcA94rInGSeEJEGoJbBvZZjtQ+XO9J9AHeB1YV+fwFuCidU9WX/nkHgHuAXfmjDkHaIyBLc1FYvA2/MHIvqZwtYhwvHFyZ7gCfYJuDcjLJkr9zmjPLtuLA6V1UfTavnp3A9xp88mgrkTTAFPg3chJub6wjwTdwkuEPGTojI5bgv8F/jPlmchrtaLYybeNYYuqPdbGvZNiiE7mjdkeptKo+Us6x6Gdcvu55lVcs4tepU5pXNs3kpTV6LBCOcWnUqp1adOqi8sbuRzc2b2XxkM1uat/Bi04us270utb+uqI4l1a5HdVnVMpZWL2VWySwb82xy6YvA50RkC6439GbgbI5ybOIw/hN3Jfo9InIbruPrX3BDCH4NqbGSpao66vhVVY2RESpFpBmIpYXQ3+Kuwn9ARD6MO619O+7ioWGnqQJ+jOuxvAHXA3uSL4+papN/7hrgKtwwgZPSnntYVRMcuy8Az4vInbhcNR+4C3hEVTf7HmjA3YpeRL6C+xq2476GV+IC6VHPq5oXwVREIrgBtB9S1cd92TuAXSJyjqo+k/GUW4Cfq+o3/OMdIrIUN9eXBdNpqCfWw5bmLWxs3Mjm5s28dOQldrfvJuF/TqsKq1havZQL5l7Asmr3R3l2yWz7o2ymjNriWmqLawcNB2jra2Nr81Y2N7uwuqV5C0/vezr1czGjYAbLqpexrHoZy2uWs7JmJXXFdblqgpl+vombYujLQCWu5/SXuLGVE8KHp0twV6b/HlDg34G/SQtyHwU+xcAp62N9TxWRq4E7gYdwPajPA5epatZeWj+N1Zn+YeaUUjuARcB1uJ7IzDGg4O4UtXcC6v6iv7HAZ3EzCRwB7gPuGOYpd+BO6d8JzMR9CHifqv7waOsgww8zOHFE5DW4AbUNqro7rXwX8C1V/ULG8WcCXelfYD/J7EdVdQYjWLt2rT733GjjhE0+iyfi7GjbwYtNL7KxaSMvNr3I9pbtqZ7QuqI6llYvdQG0ailLq5cys3imhVBjcB/i0oezbDqyKXUm4YoFV/D584ZcADxW9gNmxkVELgM2quq+tLJHgX2qOmXuZGTGJy96TIG5fr0vo3w/7lPAIKr6/9Ifi0g5bh60RzOP9fvfi5uigXnz5h1rXc0JpKoc7DqYCqAbmjbw0pGXUnODloXLWFGzgneteBcra1ayomYFtcW1Oa61MfmrKFTEqtpVrKpdlSrrifWwtXkrhaHCHNbMTEM3AQ0iciuuZ+5K4GLgkpzWyuRUvgTTYiCR5XZafbjpFYblp034JVCEG8MxhL/LxHfA9Zgec23NcdPW18ampk2pILqxaWPq6vhwIMzSqqW8ZdFbWFGzgpU1K5lXPs8mpTfmGBWFilhTl6t528009lfAP+PGepYBW4C/VNUnclork1P5Ekx7gICIhPzA4qQC0iZzzSQiNbjbgS0D3uDv1mAmiXgizsutL7O+cT3rG9ezoXEDu9t3p/Y3VDRw7pxzUyF08YzFNj2TMcZMEf6q8xtyXQ+TX/IlmCbn6ZqVtg0wm6Gn94HU3GSP4T5lna+qG45j/cwEaOtrY0PjBl5ofIH1jevZ2LiR7pi76UZVYRWra1fz5kVvZkXNCpZXL6csUpbjGhtjjDHmRMqXYLoe6AAuAO6GVPCcT5Z75opIHfA73ASu5/jbapk8ktAEu9p28cLhF1I9ojvbdgIQkACnzjiVNy18E2vq1rC6djVzS+faxUnGGGPMNJcXwVRV+0Tkm8CdItIEHMZNI/GUqj7rp5OqApr9vXLvwk0xcRHQkzaXlw43Ib85vjr7O9nQtMGF0MPr2dC0gY7+DsDdP3x17WquXHAlq2tXs6Jmhd2y0xhjjDFD5EUw9e7ATZB/t18/Ctzq952D6yF9nYj8CXf7qwDwfzNeI05+tWlKUlX2tO9hfeP61Gn5l1teRlEEYWHlQi6dfymra1ezunY188vnW2+oMcYYY0aVF/OYnkg2j+n4dUe72XRk06DT8q19rYCbrmlV7apUCF1Zu9LGhhozPdmnT2PMMbPeRTOIqrK3c2/qlPz6xvVsa9mWmry+oaKBC+svZE2tGxu6oHKBTddkjDHGmAlhwXSa64318tKRl9xped8jmpw3tDhUzMqalbx75btTPaIVBRU5rrExxhhjpioLptPMwa6Dblyo7w3d3LyZWMJNHVtfVs85s89xIbRuNYsqFxEK2LeIMcYYY04MSx1TWDQeZUvzFl5ofCHVG3qo201aUBAsYHn1cm5cdmOqN7S6qDrHNTbGGGPMdGbBdApp6mkaNDZ005FN9MX7AJhVMovT605ndd1q1tSuYXHVYsKBcI5rbIwxxhgzwILpJBVLxNjesn1gyqbD69nbuRfw95SvXso1p16T6g2dWTIzxzU2xhhjjBmZBdNJoq2vbdAFShubNtIT6wGgpqiGNbVreMeSd7C6djVLq5dSECzIcY2NMcYYY8bHgmkeSmiCna07B01gv6vN3XU1KEEWz1jM1YuuZnXtatbUrWF2yWybwN4YY4wxk54F0zww5HaejRvoiLrbeVYWVLK6djVXLbyK1bWrWV693G7naYwxxpgpyYLpCRZPxNnRtoONjRvZ0LSBDY0b2NG6I3U7z0UzFnFpw6WpCexPLj/ZekONMcYYMy1YMD3OGrsb2dC0IRVENzVtojvWDUBFQQUralZwycmXsLpuNStr7HaexhhjjJm+LJhOoORdlDY2bWRD4wY2NG3gYNdBAEIS4tSqU7lq4VWsql3FqtpVzCubZ72hxhhjjDGeBdOjlNAEe9r3sKFxQyqIbm/ZTkzdXZRml8xmTe0aVi5dyaraVSypWkJhqDDHtTbGGGOMyV8WTMcoOV1TMoRubNpIR7+7QKkkXMKK6hXcvOJmVtWsYmXtSmqKanJcY2MyqEI8CtFuiPZArAdifRDvd+Xx/rFvJ+KgiSyLZt9GB8ok4JZA0G/LQJkE07aTx/nyYASCYb/226GCLOVpSygC4WIIF0GoyL2WMcaYvGXBdIwe3PYgX/ufrxGQAAsrF3LJyZewqnYVK2tWsqBiAcFAMNdVNFNJIg79ndDXCX0dfmn3ZR1Dl2j3QOActHQP3tb4xNc1GSCRoaFSAiDpx+ADazwjxCYGh1104usJECp0ITVcPBBY09eR9O1SKCjzS3nadkZZKHJ86mqMMdOQBdMxurzhcpuuyYydKvR3QW8r9LRAj1+nPx603eYCZjJ49neO7X1CRVBQCpGSwQGruHpo6AoXQbhwcA9iuHBwD+RYtyXoezyP0xhp1bQAm4BEzPfYJntt+wb34May9fD2ux7hWK8L5f2Z4b17YOnvhK7Gwcf1d/qQPIpgwdDAWlgOhZVQNMMvlYPXyX2FFe7/0RhjDGDBdMzmlM5hTumcXFfD5IKq663saoLuI27paoLupsFlmeEzERv+NQOhtOBSCcVVMOPkgWAzqLdumCVSBsEp+iMs4kNv8tR7Du5kpuoCanqP9ZDe6mxlHdD6KvRudN8HI37IEBdikwE2+T1RXAXFNVBS4z5klNQMPC6qmrpfd2PMtGe/3cz01NcBnYeh85BbjxQ4u4+43rdsQkU+LPgwUT47e89YUeXg7Ujp8ettNBNDxPVER0qg7KSjf51Yv/+wkuwhT+85bxlc3tMCra9AT7Pbzl4x/2EmI7iW1A4tK53ptq1X1hgzSVgwNVNHrN+djk2GzUHrjLJoV/bXKCgf+KNeUQ+z1wz+Y19cAyXVA2WRkhPbRjP5hCJQWueW8YjHXEDN/LCU+fjIDnj1T24729ADCbjv19KZUFrr1iV+Paiszn1wsgvEjDE5ZMHU5L++Tug4CB37of1A9qDZecj9Ec+msNL/Aa6DOWcMbKfWde4Pd3G1Xchi8kcwNL5Am0i4nthUcG30Pxv+5yP5oa3pZbeO9w19jZVvhz//7sS2wxhjxsGCqcmdRNz9gew44AJnxwFo3z84hHYccOP4MoUKB3p8qhfCyedkBM600BnKwfhEY060QMCPTa0CFo98rKq74G7QGYbDULXghFTVGGOGkxfBVESCwGeBm4Ey4FHgVlU9NMzxa4GvAacB+4B/VNUfn5jamlElLxZKBstU4DwwuKzz0NBTj4EQlJ7kxvTVLoYFF0L5LCib7crKZ7vQWVBmYzSNOVrix6kWVULNKbmujTHGpORFMAU+DdwE3AgcAb4J/Bx4beaBIlILrAPuBd4NvAH4vogcVNXHTlSFp7X+bmjfB2173dK+D9pehTZf1r4/+xjOwkoXLMtOgrplPnD6JRk+S2ptjJsxxhgzTeU8mIpIBPgw8CFVfdyXvQPYJSLnqOozGU95D9AGfFhVE8AWETkd+ChgwfRYxWPQeXAgdKaCZ9rjbGM5S2dCxVyoWwqnvMGHTR9Ck+EzYvO/GmOMMWZ4OQ+mwBrc6fsnkwWqultEdgPnAZnB9Dzg9z6UJj0JfFNEAhnlJp2qm4JmUO9mMnD6xx0Hht4dqKACKua44DnnDLdOLuVzXAC1cZzGGGOMOUb5EEzn+vW+jPL9QP0wxz+f5dhioApomtDaTSapU+wjBM9Yz+DnBCMuXFbMhYbzBrbTg2dheW7aY4wxxphpJR+CaTGQUNVoRnkfUDjM8b1ZjmWY4xGR9wLvBZg3b97R1zSX4jF/8dC+LKfZfRAdcopd/Cn2OTBzGZxyiQ+cyfBZ76ZJsjGdxhhjjMkD+RBMe4CAiIRUNf0ejgVAtlnQexh6f8Lk46yzpqvqd4DvAKxdu1aPrbrHgSp0N7uAmTV4Jk+xZ4xSKKgYCJpzz/S9nfUDwbNsts3LaYwxxphJIx+C6at+PSttG2A2Q0/vJ4+flVE2G+jEXRSVf/q7XI9m5mn19McjnmI/f+C0ejJ42il2Y4wxxkwx+RBM1wMdwAXA3QAiMh+YD/w+y/FPA+8UEVHVZO/n64A/5uTCp+Qp9qzTJiWvYs+857UMXMU+czksviwtePqxnXaK3RhjjDHTTM6Dqar2icg3gTtFpAk4jJvH9ClVfdZPJ1UFNKtqP/B94OPAt0Tkn4GLgWuBy45rRQ9tgt1PZwTPfdlPsRdWQLkPmHPP9KFz7sBpdzvFbowxxhgzRM6DqXcHEMb1mIbxd37y+84BfofrFX1SVQ+JyGXA13FX5+8BblTVJ45rDXc+Bev+FoIFA6fSGy4YGM+ZHjwLyo5rVYxJUlWicaU/niAaS9AfT9Dv19F4glhciSWUeCJBPAGxRIJ4QlNLLKEkEsljBpfHVYnHE+4Y9WVxRVPvDYqSPG+hyUK/ndw/cCxpx/oCAIGACIHUWgYeB2TQPhEIppWJiH/stgMiBAMQCgQIB4VQIEAoKESCAUJBtx32Zcn94VCAcECG7A8FBLG7ixljzAklA2fDp4e1a9fqc889N/4n9rRCPAolNXYrTDNINJ6gNxqnN+rWfbE4Pf0JemPxVHlP1G33pR3Xk9yOJfcl6EuGyrRw2Z8WOKOptabK803yx0NwYVHSygWBQft9aFVIqAvD+fQrKRleI6EABaEABeEABaGg2w757XDadih57NDyzOcWR0IURQIUhUMURYIUR4IUht06HJyUw3jsF6Mx5pjlS49p/iuqzHUNzDFQVXqjCbr6Y3T3xemOxujqi9Pdn7buj9OT8bi7L1keTz23qz82KIjGEkeXpIIBoSgcpNAHluQ6HApQEAxQGA5QXhgiHHRhJ+LX4Yx1JChDygv8OhQQwsEAgYDrAQymLSHf85jsHQyIC2HBoBCUtGPSnpveUwkD4XIiexbVh9P0oJpQJaEQTyiaZTuh6h9DXJVY3IX3WMKtk73H0YTvRY4niCaSx/ljwegY0QAAEWBJREFUfe9wNG2/e97Ah4Pkh4e+WNytowlauvp9mTuuLzb4Q8bRCPnvjaKIX/x2cWo7RFHYhdtkmE0eU1IQpCQSoqTAL5GgX4coKQgSmpyh1xgzTVgwNXkrGSY7+qJ09sbo7IvR2RujI7nujdLZN/A4c39PNE5XX4xuHyrH0xNXFHZ/4IsjIYr9H/bSghB1ZQWpMJAMlYV+XRQOUhB2vV6FoQBFkeR2+nED25O0V+y4E3/KPoBM+l9QiYSmBdp4KsAme8x7+uN097se8+7+ZFksbdut049p6YqmPTdGb3R8ATgSClBa4L6vk+tkcD33lBpuOPvk4/g/YowxI5vsv/dNnoonlM7eGO29Udp6orT3RGnvjdKRFiDTQ2UqZCb3+/Kx9EYWhAKUFbrgWOrXsysLKfY9RMUR12tU7P8Ipz8u8T1SJZEQxb6nqSgcJBCws5Lm2AUCQmHAfSBxw+ePj1g8kQqwyQ9jnX0xuvtjdPa5nv/OtA9pXX2u9z9Z1tEb41B7Lw21JcetjsYYMxYWTE1Wqkp3f9yFyt4o7T0x2nuigx4PbEd9AI2ltjv7Ru+hLAwHKC0ID4TKghDzqoopLQxRlgqZYUoLQ5SnHeP2h1MhNBKynkczvYWCAcqCAcoKj1/4NcaYE8GC6RQWiydo64nS2hOltTstRPZEae+NpQVKHzz9/ja/Pz5Kb2VpgQuM5UVhyovCzKksYumsMsoLw1T4svLCUGq7rDBEeWE4FS7tVLYxxhhj0lkwnQR6o3Ha0wJma3c/rT1R2rqjtPb0u7KMx23dUTr6YiO+bkEoMChAVpVEaKgpobwwTHmRD5SFyf3JsOnCZVlhyC6iMMYYY8yEsmB6giRPjbtw2e9DpA+aPf5xerD0PZet3e5Ch+EEA0JlUZiK4jCVRWHqygpZXFfmH0eoLA5TWezCZUVawCwrDPlxb8YYY4wx+cGC6TglEkpHXyxLb+XAtguWQx9H48OfGo+EAszwYbKiOMy8qmIfKiNUFLlwWVE0EDaTZaUFIZsE3BhjjDFTggXTMbr72T18+bGttPVEGWnoZUkkOChMLp5ZSkWy5zIVMAd6MpNB03ovjTHGGDPdWTAdo4aaEq5cNTuttzKSCprJsFlRFLYrxI0xxhhjjpIF0zE6d1EN5y6qyXU1jDHGGGOmLOveM8YYY4wxecGCqTHGGGOMyQsWTI0xxhhjTF6wYGqMMcYYY/KCBVNjjDHGGJMXLJgaY4wxxpi8YMHUGGOMMcbkBQumxhhjjDEmL4jqCPfXnIJEpBHYc5RPrwGaJrA6uWbtyX9TrU3Wnvx2LO1pUtXLJrIyxpjpZ9oF02MhIs+p6tpc12OiWHvy31Rrk7Unv0219hhjJh87lW+MMcYYY/KCBVNjjDHGGJMXLJiOz3dyXYEJZu3Jf1OtTdae/DbV2mOMmWRsjKkxxhhjjMkL1mNqjDHGGGPyggVTY4wxxhiTFyyYjkJEGkTkYRFpFZEDIvJdEanMOOY2EdkjIt0i8riInJKr+o6FiERE5MsiclBE2kXkVyLSkHHMpGpTkoh8TESGjE+ZTO0RkdNF5Df+e26/iHxPRKrS9gdF5HP++7FTRH4mIjNzWeeRTLb6ZhKRmSLyI1//VhFZJyIr0vZfJyJbRaRHRJ4VkTNzWd/xEJGzRSQmIhemlV0iIi/49mwQkctzWEVjzDRjwXQEIhICHgHiwJ8Bfw68Fvhu2jHvBv4B+AhwFtADPCoiBSe8wmP3beDtwLXAOUAR8LCICEzaNiEiK4F/zFI+adojIrOB3wC7cN9zfwG8Bngg7bBPAzcBNwLnA3OBn5/Qio7Pp5lc9U0RkQDw78Bi4M24n5c24LciUi0iFwM/AL4MnA5sBB4TkdocVXnMRKQE+AkQTCtbBjwMPAicBjwE/FJEluekksaY6UdVbRlmAVYACqxMK/sA0J72eCvw6bTHpUAHcG2u6z9Mmxb4Nl2UVnYq7m5YiyZjm3wdI8B64Hfu23rQvknTHuA2YD8QTCs7z3/N5vl2tgM3p+2f7/efk+v6D/N1mTT1zVL/03xdl6aVFQBduKC9Dvhh2r4AsAP4u1zXfQxt+3by5wW4MK3syYzjfgd8J9f1tcUWW6bHYj2mI2sGEsB7RaRQRGpwPY3PAYhIHa4n5cnkE1S10+8/74TXdmwuARpV9YlkgapuVdWTVfXlSdomgM8C+4DvpxdOwvY8DFyjqvG0soRfzwDWAGUMbs9uYDf52Z7JVt9MrwBX4j7cJCUAwX09zmVw2xLA78nztonIG4ErgA9l7DqPtPZ4T5Ln7THGTB0WTEegqvuBDwLvxPWQNAJ1wDX+kLl+vS/jqfuB+hNRx6OwGNgpIteKyHo/hvFBEUm2ZdK1SUTOw32N3pNl96Rqj6ruUNU/ZBR/Alf/F5lk7WHy1XcQVT2iqr/ygTPpQ0Ah7sNNCZOsbf4D9vdwPy8tGbvnMsnaY4yZWqZ1MBWR+SKiwyy9fnzZEtyYv3OBS4EYcL+IBIFi/1K9GS/dh/vDdcKN1iagHNemj+BOG/8FMBM3Zq6QPGvTGL5GZcCPgQ/5DxKZJlV7shz/eVyP3Qd8L2oxkFDVaMahOfueG8Vkq++IROQq4HPAV3DDXyBPvrfG4dvAf6jqo1n2FTP52mOMmUJCua5Aju0Dlg6zLwFch7tA6GRV7QIQkatxY8jeiOtJADfmLF1yDFoujNam24AK4G2qugtARN4GHMC1KfnHNl/aNFp7vg48p6o/HeaYHr+eLO0B3JXswDeA9wHvV9WH/a4eICAiIVWNpT03l99zI5ls9R2WiNyMu/DxPuDjuFP5kD/fW6MSkZtw42ZXDXNID5OoPcaYqWdaB1Pfi7NluP0i8mFgSzKU+ufsFJEmYBHwJ188C3g57amzgc0TX+PRjaFN+4CuZCj1zzksIkeABuBpX5wXbRpDe24GekWk0xeFfHknLtQ97ssnRXsAfM/1A8BlwPWqem/a7lf9elbaNrj2ZJ6CzQeTrb5Zicjf48YxfwPXO68i0owLbLMyDs/ntt2MO11/MDkJhy//tYj8CPc1mkztMcZMMdP6VP4Y7AUWp08rJCKzgGpgu6oeBrYDF6TtLwXW4i6AyEd/AEpEJNVrJyInATXAjknYplOAlbiLbNYAf+/L1wAPT7b2+OEjDwKvB96UEUrBzTzQweD2zMdd6Z537WHy1XcIEfk4LpR+UlU/qKoKfuoHeIbBbQvgpsTK17ZdDyxj4OflUl/+HuCTuA+mF2Q853Xkb3uMMVNNrqcFyOcFmIO7OOABYDlwJu4X9/NAyB9zC9AJvAM3vdTDwDYgkuv6D9Mmwf2ReR43T+Yq4Alc72FkMrYpo33XM3S6qEnTHuBW3PQ97wZOyljC/pjPAwdxPaqnA8+SMcVPPi2Trb4ZdV+FG1f+/SxfjxLfpqj/ui3FnepvAmpyXfcxtm8ug6eLWgn04+b9XQJ8Bnd6f2ku62mLLbZMnyXnFcj3xQeZX+OmjjoA/AiozTjmdtx4007gUWBBrus9SpsqcVflNuN6s/4dmDuZ25RW7yHBdDK1B9cDp8Msr/XHhHATujfhJnu/P5+D0GSrb0bd/88IX487/DHvBHb6APdH4PRc13sc7RsUTH3ZFcAm3EVQLwAX57qetthiy/RZRHXI3RuNMcYYY4w54WyMqTHGGGOMyQsWTI0xxhhjTF6wYGqMMcYYY/KCBVNjjDHGGJMXLJgaY4wxxpi8YMHUGGOMMcbkBQumxhhjjDEmL1gwNcaMSkTqReS3IrJZRDaJyBfF32zdGGOMmSgWTI0xYxEDPqGqS4HTgLOAt+a2SsYYY6YaC6ZmWhGRsIjcJiL/LSKdItIuIutE5KwJfh8RkfUictMx1vP/ikibiPT4Ot8mIpGMY98uIgePZw+mqh5Q1ef8dj+wAahPq8NdIvL94/X+xhhjpge7JamZNkSkEngMWAB8DXgWqAZuA9YAb1DV30/Qe10DfAlYqKrRcT53BvAbYCHwL7j7rwNcDtwC3KCqD6Qd/yMgrqrvmoi6j6F+1bh7qF+iqpt92XxgC7BCVV8+EfUwxhgz9VgwNdOC7018Ahf2LlDVXWn7CoFtwGFVXTtB7/dH4ElV/fujrOcS4HWquiVj/1rgSLL+IhIADgAfUNWfT0TdR6lfAfAo8J+q+uWMfb8B1qvqR453PYwxxkxNdirfTBfvBi4E3p0eSgFUtRe4GzhDRGqP9Y1EZBFwDvCzjPKAiHSJyC0i8lkR2etP03/XB0yAm3w9b8kMpb6uz2XU/0xgBvB42vsE/Sn/DSLSKyIHROQeP7wgWYcPi8hXROSwiLSIyMf8c28QkZf8MIdfiEhR+usC9wDPZ4ZS7+fAdWltMcYYY8bF/oCY6eKjwP+o6uPD7H/Vr2dNwHu9HugC1meULwCKfV2KgHcCdwLvAd7sj/kbYLOqPjTG97oC+IOqtkOqB/VnwCdxIfJK4O+AsLrTI8k63Ab0AdcCvwK+KCLfAN4BfAy4Hbja1zHp20AHMFyP6DPATGDlGOtujDHGDBLKdQWMOd5EZDFwKsMHKoAav+6egLc8AxcuExnlq/z6LlX9qt9+XETeD5wiIifjQt0d43ivK3ABNOk24A3Amcnxn96/ZdThK6r6dQAR2Q5cBywFLvYBFhF5L+7/DRE5F9fr/CLwvL/O6gfJ1/A2AXHgNQwN5cYYY8yoLJia6WCFX780wjFLgF5g9wS830lAU5bylUA7cFeywI8prfTHJ3saXxzLm4jILNzUTdf6xwHg48A3MkJpZh1agX9NKyvx68/r4EHnJUAzgKr+ERjxqn9VjYlIK679xhhjzLjZqXwzHZT5dbawiIiUAFcBj6pqbALerxB3mjzTStxp9/60soW40/qbgApfdmiM7/NGYKeqbvWPVwF1wL+P8JyVwNMZMwWsws1TmpqRQESKgfmMMSSn6cO13xhjjBk3C6ZmOjjg1/XD7P8AUApku6DnaDTjekEzrWLoKe7VQAIXAA/7stljfJ8rcONDk5LjYw9kOTa9Di9kqcMWVU0P0ytxvx82jLEuSZX4XlZjjDFmvCyYmungGdwp9CHzfIrIBcBngO+q6tO+7Nt+eUJE9ovIzSLycRF5VkR2iMiSUd5vK9CQ8T5FuN7RzFC4Ctfr2QX8l6/nO8lCRF6bth0BLmZwMD3o10uHeX6yDpnhOFtgXoW7gGtHttca5vVrcRdWbRvrc4wxxph0FkzNlKeqncAngCtF5H4ReZOIXCEiXwPWAQ8Dt6Y9ZQ3Qg7u6/lrcJPcbVfVs4EFGvxXnH4F5GVNPLcf9vGULgBuy1PMhEXmbiFwgIu8SkV8DX0h73vn+9Z5KK9uIGxLwryJyvYhcJCJ/JSLJ541Uh2xlm7JcwDWStYDiPggYY4wx42YXP5lpQVW/JSKHcVMh3Yfr2esH3qmq9yaP8xcQLQJer6oqIgo8q6q/9ocEGP1U9ZP+mMuAn/iylbgr/jPvirQK+HFGPQ/gppT6oS9+BTf+8wdpz7sC+E366Xd/8dGbcFNQfdW3cTvw+Yw6pHpB/V2m5jL0lP2qLGWjuQx4SlWPjPN5xhhjDGB3fjLTlIhcBvwauFFVf5JWvhT4oaqe5R9/CKhW1U/5x+uAf1DVEXsFfW/sIlW94jjVfxvwJVX97vF4/fHyk+/vAW5X1btzXR9jjDGTk53KN9OSqj4K3At8NeOU+xrg+bTHp2U8XsnYehK/BFzo51CdcKq6OF9CqfcXuOEP9+W6IsYYYyYvC6Zm2lLV61S1RlUb04pXM/gCpdOSj0XkJKDLjwUd7bX34iakn4g7SU0Ggrvd60RMt2WMMWaaslP5xhhjjDEmL1iPqTHGGGOMyQsWTI0xxhhjTF6wYGqMMcYYY/KCBVNjjDHGGJMXLJgaY4wxxpi8YMHUGGOMMcbkBQumxhhjjDEmL1gwNcYYY4wxeeH/A0/spSTEjL7mAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plotZeq(bls, gas, charges)" ] } ], "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.7.3" } }, "nbformat": 4, "nbformat_minor": 4 }