{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from scipy.linalg import eigh_tridiagonal" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$$ \\begin{bmatrix}\\frac{1}{\\Delta y^2}+mL^2V_1 & -\\frac{1}{2 \\Delta y^2} & 0 & 0...\\\\ -\\frac{1}{2 \\Delta y^2} & \\frac{1}{\\Delta y^2}+mL^2V_2 & -\\frac{1}{2 \\Delta y^2} & 0... \\\\ ...& ... & ... & -\\frac{1}{2 \\Delta y^2}\\\\...0 & 0 & -\\frac{1}{2 \\Delta y^2} & \\frac{1}{\\Delta y^2}+mL^2V_{N-1} \\\\ \\end{bmatrix} \\begin{bmatrix} \\psi_1 \\\\ \\psi_2 \\\\ ... \\\\ \\psi_{N-1} \\end{bmatrix} = mL^2 E \\begin{bmatrix} \\psi_1 \\\\ \\psi_2 \\\\ ... \\\\ \\psi_{N-1} \\end{bmatrix} $$\n", "\n", "$$ \\psi_0 = \\psi_N = 0$$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Define what $N$ and $dy$ is" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "N = 100\n", "dy = 1/N\n", "y = np.linspace(0, 1, N+1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Define potential $mL^2 V$" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "#np.where(((y<0.25) | (y>0.75)), 1e8, 0)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "def mL2V(y):\n", " #return 1000*np.sin(20*y) * y**4 \n", " return np.where(((y<0.25) | (y>0.75)), 1e8, 0)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "d = 1/dy**2 + mL2V(y)[1:-1]\n", "e = -1/(2*dy**2) * np.ones(len(d)-1)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "w, v = eigh_tridiagonal(d, e)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnEAAAFoCAYAAAA4vK32AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAxYklEQVR4nO3de5hlVXng/+/bXd2AgGDAtMit0bQGBQ3aAyjJz0o0EYmCSTQD4w2DdiajuYyX/ExM1JDxiZdxouaHKKMOwUxENIZfJ6JokIoOBgOKcmmE9ABCtyJ3tLl1Vfc7f+xd04fqVZdT55x9TtX+fp6nnjpnn13nrKq1zql3r/WutSIzkSRJ0tKyYtgFkCRJUvcM4iRJkpYggzhJkqQlyCBOkiRpCTKIkyRJWoIM4iRJkpYggzhJWiYiYiIiMiIGvnZURNxSv9Ytg34tSWVjwy6ApOVjnuBhG/Aj4NvA54G/y8zJAZblD4D9gfsy84ODep1BioiXAj9X3/1gZt43tMJIGjnhYr+S+qXLHqBrgJdl5o0DKsstwOHA9zNz7SBeY9Ai4lzgNfXdIzLzlnnOnwCeB5CZMeCy3cIS//tKS509cZIG5ddm3H8c8FzgFcBewNHAVyLimMy8p+nCLUeZOT7sMkhqjkGcpIHIzAsLh/9HRHwQmAAOBA4D/t/6S5LUBSc2SGpUZl4H/HHHoZcNqyyStJQZxEkahi903H5SRDxm5gkRcWREfCgiro2I+yPioYj4fkRcEBEzh2o7f+6WOjfv8PrQ4dMzNmd8nT7Lz6+OiDMiYmNE3BYRD0fEfRFxdUR8ICLWzvWLRcS5Ha+xtj72woi4MCK2RMQjEfGDiPhsRBw313OwKx8O4ObC73DujJ+bd3ZqRDw+Il4XEX9T/21/HBGTEXFXRFweEWdGxBPm+h0ljQaHUyUNw50z7u8PPDh9JyL+DHg7sHLGeYfVXy+vk/h/o5/5dBGxHrgAOGLGQ3tQ5fAdDbwxIn4vMz+2gKdcEREfAX5nxvGDqHogfz0iNmTmJ3os+oJExJOAG9n97wpwQP11HPCmiHh1Zn6+iXJJWhyDOEnD8PgZ9388fSMi/gJ4W313B3A+8FXgIaog6reANcA4cGlEHJeZD3c81wbgMcA59evcWR+b6duddyLiOcA/1T+bwMXAl4GtVBMxngO8qn78oxHxSGaeO8/v+V+A06gCp/OAzcC+wK8DL6IaDflIRFyWmd/r+LkPAxcCvwf8Yn3st4E7Zjz/rfO8/kyrqQK4m4BLgGup/j4rqILjFwC/BOwNnB8RJ2TmFV2+hqSGuMSIpL7pHMaba4mLiHg9VZAFcEtmHlEffw5wGRDAA8BJmfm1GT/7U1QB1vr60H/NzLcWXuMWFrgERkTsC1wHHArcB7w0M/+5cN7PUAU/h9XlW5uZd80451wePQx6HnBGZk7NOO9DVEEawNmZ+Z8Kr9f5XD0vMVL/7X42M78xx3P8IvAPVIHcRGb+4izn3YJLjEhDZU6cpEZFxM8C7+449LmO22+lCuAA3jozgAOoh09fxq7h1/8YEfv3WKzXUwVwAK8uBXD1a28GXlvf3ZtyD1+n7wGvnxnA1f6EqncR4IXdFXdxMvOeuQK4+pxLgQ/Ud8cj4tC5zpc0PAZxkgYiIl464+v0iDgH+Ba7hlN/ALyvPn8P4KT6+N3ArHlimfl94NP13X2AX+mxuK+qv9+Ymf8w14mZ+VWqcrOA1z07M7fP8jw/Aa6s7x4REXsutLAN6Az0jh1aKSTNyZw4SYPy9/M8fj3w8sycnuTwTKoJBFAN4xWDnw5fBs6obx9HNSGhaxGxH/CM+u6P6q2u5rOt/n7kPOddPs/jW6eLQTW54/YFvHbPIuLpVMO0JwDrgP2o8uVKDmmiTJK6ZxAnqSkPUCXmX0UV4H02Mx/pePygjtsL2Yqr85yDZj1rfoeya1TiF+qvhXrcPI/fNc/jnb//wHviIiKA9wBvYeEjMY8dXIkk9cIgTtJALGLvzn07bj+wgPO3ddzed9az5rdfDz+7ap7Hd/bw3IPwx8Af1rd3UM3G/QbVLNcHgMn6saOAP69vl5YjkTQCDOIkjYqfdNzeewHn7zPLz3arMxg8LzNfM+uZS1hE7AX8UX33J8AvZua3Zjl3snRc0mhxYoOkUfHDjtvrFnB+5zk/mPWs+W3tuL2c87+ew67g+GOzBXC1w+d4TNKIMIiTNCq+y64csfGImG+osnNm6L8WHp8eypxzWLde521Tfff4iBilHLDO4dhuh6dnWtNx+3/Pc24jS55I6k0rg7iI+GRE3BER1y7g3MMi4tKIuKreO/Gk+X5GUvfqSQ7Te6oeCJw+27n12mWn1Xe3US3+O9P0MOlChmb/uv7+GHbtFjEKOod6F/J7zOXBjttPnu2kiHgW8OIeX0tSA1oZxAHnAicu8Nw/AS7IzGOAU4GPDKpQkng/u3qfPhARJ8w8ISIeR7VA8HRQ89HMvL/wXDfX3w+IiMPmed2zgO/Xt98WEW+NiFk/HyNiv4j4vYh4wTzP26ubO24/q8fnurLj9usiYub+sETEOuDvaO//BmlJaeXEhsz8WkSs7TwWEU+m+iB/PNUV6+vrvQyTXVPs96O33BtJc8jMyyPivVQJ+PsC/xwRn2bX3qlHAa9j19Dg1cA7Znm6S4CT69ufj4izqfLupoPEazJza/26D9Trw/0z1fv9fcBvR8TfUQ21bquPP4lq8dtxqnXVphcJHpRLOm6/LyIeD9wATO8AsTUzr1nIE2Xm1oj4PNW+rfsD342Ij1H9DVcAzwVeTbXUyXn1bUmjLDNb+QWsBa7tuH8JsK6+fRzw1fr2QcA1wBbgXuDZwy67X36N6hfVRU9WHy09Pc+ZVIFKzvE1ARwwx3PsQxXwzPbzpxd+5qnAt+d53emvh4ETC89xbsc5a+f5Pec9F/jbOcpw7oxzJ+b6+wMHUAVtsz3fDqrRh/GOY++a5bluqR+/Zdjtzi+/2vpllzkQEftQXYV+NiK+A3yMXYuHnkb1QXkI1ZZAn5prmEVS7zLzHVS7KPwVVU/YT6gmPWyhGu77jcwcz8y753iObcDxVPu0fhu4n3nWbcvMG4BnA6dQ5cndCPyYKri5j2ryxXlU+XoHZeaXFvs7duFVwO9QBWh3sasXrmv13+t4qp7Oq6hGHR6kmujwP4DnZuZ/6bG8khoSmTnsMgxFPZz6j5l5VD0b7YbM3G3V94i4jupq+7b6/k3A8Zl5R6MFliRJ6mCPEpCZPwZujoiXQ7U1TUQ8s374VuD59fEjqfJF7iw+kSRJUkNa2RNXJ0qPUy1j8CPgnVSJ02dTDaOuAs7PzDMj4mnAf6fKr0ngDzPzy8MotyRJ0rRWBnGSJElLncOpkiRJS5BBnCRJ0hLUusV+DzzwwFy7du3AX+eBBx5g77173SVH/WSdjCbrZfRYJ6PHOhlNTdTLt771rbsy8/Glx1oXxK1du5Yrr7xy/hN7NDExwfj4+MBfRwtnnYwm62X0WCejxzoZTU3US0R8f7bHHE6VJElaggziJEmSliCDOEmSpCXIIE6SJGkJMoiTJElaggziJEmSliCDOEmSpCVopNeJi4hPAi8G7sjMowqPB/Ah4CTgQeD0zPx2s6Xc5bZ7HmRqZ7UX7e0P7OTmux4YVlFU0O862Xv1Sn76sXv27fkkLQ/bp3ay9b6H+vqc/k8ZTdP18lN7r2a/vVY1/vojHcQB5wL/H3DeLI+/CFhXfx0HnF1/H4pXfPyb3HrPg7sOfH1iWEXRbL4+0dene/oTH8snXvPveMJ+BnOS4H9+8/u8+wvX8+D2Hf1/8q9P9P851buvT/DnLz2KVx1/eOMvPdLDqZn5NeCeOU45BTgvK5cD+0fEQc2UToLrfvBj3n3R9cMuhqQR8MP7H+Ltf3/tYAI4qWDUe+LmczBwW8f9LfWxH3aeFBEbgA0Aa9asYWJiYiCFeeih/nafa2n4xg0/ZGLi/mEXY8natm3bwN6TWhzrZHGuuH1q2EXQkPzbjTcy8fDNjb/uUg/iFiQzzwHOAVi/fn0Oap+zvf71UnjowflP1LISY6vd07AH7gk5eqyTxbn3qi3wne8OuxgagnVPeQrjQxhOXepB3Fbg0I77h9THhuLg/fdiRVS3H3roIfbaa69hFUUF/aiTHZncds+je1yndu7s6TklLQ+TO3K3Y/vsMcaB+6zu+bn9nzKapuvlsXsOJ5xa6kHcRuCNEXE+1YSG+zPzh/P8zMB8esPx//e2V7Kjpx918pOHJzn6XV9+1LHJKYM4STC5Y/fPgpc884n8xa8f3fNz+z9lNA27XkY6iIuITwPjwIERsQV4J7AKIDM/ClxEtbzIZqolRl47nJKqLVat3H0u0OTO3a++JbXPVKEnbtXKGEJJ1BYjHcRl5mnzPJ7AGxoqjsTYit0/kKcKV9+S2qfUEze2YqQXgdASZ+uSurCyEMTtTNhhb5zUelOFzwF74jRIBnFSFyKC1aUhVXvjpNYr5ceWUjCkfrF1SV0aK1xZl67AJbVLKT+29Hkh9YtBnNQl8+IklZQ+B+yJ0yDZuqQuFWeoFmalSWqXUo986aJP6heDOKlL5SDOnjip7babE6eG2bqkLhVz4uyJk1qvtHuLs1M1SAZxUpfKC/7aEye1XelibsyeOA2QrUvqUnligz1xUtuVcmPNidMgGcRJXTInTlJJ6XNg9Zj/ZjU4ti6pS6UcF4M4SaWcOLfd0iDZuqQulXJcXOxXUnE41YkNGiCDOKlLpRyX0nY7ktql1CPv7FQNkkGc1KVSjktpux1J7VKa4OQ6cRokW5fUJbfdklRS6okzJ06DZOuSulTKiXPbLUml3FiHUzVIBnFSl5ydKqmk2BPncKoGyNYldamU41JaWkBSu5R65O2J0yAZxEldKuW4OJwqqZQb68QGDZKtS+pS6crabbcklXLi3HZLg2QQJ3WptHinOXGSthfWi7QnToNk65K65N6pkkpKubEGcRokW5fUpfLEBodTpbYrpVW47ZYGySBO6pKL/UoqKW675WK/GiBbl9Sl0rpP253YILVeaZa6PXEaJIM4qUuri7NT7YmT2s6cODXN1iV1qdQTZ06c1G6Z6WK/apxBnNSlUk6cs1OldttRuJBbuSKIMIjT4BjESV1yiRFJMxXz4VzoVwNmECd1qbjEiBMbpFabLOTDrTYfTgNmC5O6VN6xwSBOajPXiNMwGMRJXSrunVq4CpfUHqUZ6qVJUFI/2cKkLpkTJ2mm7cWFfu2J02AZxEldGiuswO5wqtRupeHUVWP+i9Vg2cKkLhWHU+2Jk1qtlFLh7FQNmkGc1CUX+5U0U3mhX//FarBsYVKXSj1x26fsiZParJQX6+xUDZpBnNSl4jpx9sRJrWZPnIbBFiZ1qZTnYk6c1G6lz4BVhUlQUj/ZwqQulZcYsSdOarNSb7zDqRq0kQ7iIuLEiLghIjZHxNsKjx8WEZdGxFURcXVEnDSMcqpdXCdO0kyldeJc7FeDNrItLCJWAmcBLwKeBpwWEU+bcdqfABdk5jHAqcBHmi2l2qh0dW1OnNRupXXiVtsTpwEb2SAOOBbYnJk3ZeZ24HzglBnnJPDY+vZ+wA8aLJ9aqpTnYk+c1G7FbbfMidOAjQ27AHM4GLit4/4W4LgZ57wL+HJE/C6wN/CCZoqmNiv2xJkTJ7XapDlxGoJRDuIW4jTg3Mz8QEQ8B/hURByVmY+6JIqIDcAGgDVr1jAxMTHwgm3btq2R19HC9atO7n9k9w/rBx5+2PpeJN8ro8c66d41Wyd3O3b3XXf07e9onYymYdfLKAdxW4FDO+4fUh/rdAZwIkBm/ktE7AkcCNzReVJmngOcA7B+/focHx8fUJF3mZiYoInX0cL1q07ue3A7XPqVRx2LFWPW9yL5Xhk91kn3fnTFrXDNNY86dugTD2J8/Jl9eX7rZDQNu15GecD+CmBdRBwREaupJi5snHHOrcDzASLiSGBP4M5GS6nWcdstSTOVlhlydqoGbWRbWGZOAW8ELgaup5qFel1EnBkRJ9envRl4fUR8F/g0cHpm+t9UA1Va7NeJDVK7lT4DVhU+K6R+GuXhVDLzIuCiGcfe0XF7E3BC0+VSu8222G9mEuGHttRGpclN9sRp0GxhUpdWrghKF9g7HFKVWmtyZ6EnziBOA2YLkxbBvDhJnUo9catcYkQDZhAnLUIp16W07Y6kdijlxLnYrwbNFiYtwqqxQk+cC/5KrVWenWpPnAbLIE5ahNIVdmnbHUntUHr/rzYnTgNmC5MWoZTrUtp2R1I7lHJi7YnToBnESYtQ+nCenLInTmqrUk6sS4xo0Gxh0iKUlg6YKiwxIKkdSsOpLvarQTOIkxZhVSEnrpTYLKkdykuM+C9Wg2ULkxahNJzq7FSpvUo5sebEadAM4qRFKOW6uE6c1F6lnFh74jRotjBpEVYXe+IM4qS2KuXEjpkTpwEziJMWobhOnEuMSK1VyoktLQou9ZMtTFqE4hIj9sRJrVXqiStNgJL6yRYmLUIp18XZqVJ7TU45sUHNM4iTFqG0Y4M5cVJ7TZZ64gziNGAGcdIilGanuu2W1F6uE6dhsIVJi1Baid2eOKm9SjmxpQlQUj/ZwqRFKPbEGcRJrVV6/zucqkEziJMWwYkNkjqVlhgqXexJ/WQLkxbBiQ2SOpVz4uyJ02AZxEmL4GK/kjqVh1P9F6vBsoVJi7BqbPcrbPdOldqrPLHBnjgNlkGctAilldhLwymS2qH0/jcnToNmC5MWobQSuzlxUnuVFvtdbRCnAbOFSYtQnJ1qTpzUWuWeOIdTNVgGcdIilGadTU7ZEye1UWaWlxgxJ04DZhAnLYKzUyVNK60RObYiiDCI02AZxEmLUOyJMydOaqWpQj6cy4uoCbYyaRFKs86cnSq1U7Enznw4NcAgTlqE8rZb9sRJbeRCvxoWW5m0CMXhVHPipFYqzkx1UoMaYBAnLUJxYoM9cVIr2ROnYbGVSYtQyncp5cVIWv5KM9NLvfVSvxnESYtQWondnDipnYr7ptoTpwbYyqRFKG67VVhmQNLyVwzizIlTAwzipEUo5cQ5nCq1U2liw+ox/71q8Gxl0iK42K+kafbEaVgM4qRFKM08c7FfqZ3Ki/3671WDZyuTFqE8O9WeOKmNyttu2ROnwRvpIC4iToyIGyJic0S8bZZzfjMiNkXEdRHxt02XUe1U7IlzsV+plUq98K4TpyaMDbsAs4mIlcBZwC8DW4ArImJjZm7qOGcd8EfACZl5b0T89HBKq7Yp5bvYEye10/ZiTpxBnAZvlFvZscDmzLwpM7cD5wOnzDjn9cBZmXkvQGbe0XAZ1VLmxEmaVu6JczhVgzfKQdzBwG0d97fUxzo9BXhKRFwWEZdHxImNlU6tVgri7ImT2qmcEzfK/161XPRtODUingWcBjybKtjaC3gI2ApcCXwmM7/Vr9erjQHrgHHgEOBrEXF0Zt43o2wbgA0Aa9asYWJios/F2N22bdsaeR0tXD/r5MeP7H7l/dDDj1jni+B7ZfRYJ925ZuvkbsfuvvOOvv4NrZPRNOx66TmIq3PXPgr8FhDAD6gCt9uBPal6y8aBN0fEucCGzNyxgKfeChzacf+Q+linLcA3M3MSuDkibqQK6q7oPCkzzwHOAVi/fn2Oj48v/BdcpImJCZp4HS1cP+vk/gcn4dIvP/rgijHrfBF8r4we66Q7P/zXW+Gaax517JAnHsT4+DP69hrWyWgadr30o7/3bcAZwF8Ch2XmIZl5XGY+r/5+CHAY8EHgtfX5C3EFsC4ijoiI1cCpwMYZ51xIFSASEQdSBYw39fbrSPNbNVaY2OC2W1IrTRX3TjUnToPXjyDutcDHM/MtmbmldEJmbsnMNwOfpOqxm1dmTgFvBC4GrgcuyMzrIuLMiDi5Pu1i4O6I2ARcCrw1M+/u8feR5lWaeebEBqmdSov9mhOnJvQjJ+4Q4PIFnvsN4BULfeLMvAi4aMaxd3TcTuBN9ZfUmNLMs6mdSWYS4RW41CYu9qth6celwm3ACQs89xeo8tikJS0iWFlcK87eOKlt3HZLw9KPVvZx4LUR8eGIeHLphIh4ckT8FfDq+nxpySv3xpkXJ7VNaXmhVYWLPKnf+jGc+n7gcOANwBsi4m6qWaSPAHtQLTdyQH3uOfX50pK3asUKHubRH972xEnt47ZbGpY5W1lEPH++J8jMnZn5n4BnAu8Bvg2sBA6sv3+7Pn5MZv5OZtpVoWWhNPusNEtN0vJWmpnucKqaMF9P3Fci4vPAmzPz+3OdmJnXAm/vW8mkEVf6kLYnTmqfySm33dJwLORS4deATRHxrojYc9AFkpaK1W69JYlyLuyYOXFqwHxB3BuAu6m20PpT4HsR8bKBl0paAorDqTvtiZPaprhO3JjDqRq8OVtZZp5NtY3VB4Epqp0XPhMRl0TE0wdfPGl0la60zYmT2qf0vl9VWBBc6rd5W1lm3p+ZbwKOAr5AtT/qOHBVvazI/gMtoTSiSrPPthvESa1TSqNw2y01YcGXCpn5b5n5EuCFVNtgjVENt94YERvCZerVMqUgzq23pPaZLKRRODtVTei6lWXmV4BnAL8L3EO1lMjZwJUR8dz+Fk8aXeWcOHvipLYpDaeutidODVjUpUK9NtxZVPlyH6bKlzsG+HpE/E1EPLGPZZRGUinnxSVGpPYp9cCPmROnBvTUyjLzvsz8A+Bo4ItU+XKnUc1ifVtErOq9iNJoWjVW2jvVnjipbUq5sObEqQl9uVTIzBsz81eBFwHfA/YB3k21vtyL+/Ea0qgpXWmbEye1j9tuaVj62soy82J25cvdCzwZuLCfryGNitKK7PbESe1TyoU1iFMT5tt2a14RsQJ4CvB0qmVIjqpvP3b6lF5fQxpFxZ44F/uVWqeUC+twqprQVRAXEU9i92DtqcDqztNm/Ni9vRRQGlWlFdntiZPap/S+d7FfNWHOIC4i3syuYO1I4DEzT+m4PQncAFwNXDP9PTO39K200ghZVdixwdmpUvsUZ6faE6cGzNcT934g2b13bSu7ArXpoO36zJzqewmlEVVcJ86eOKl1Js2J05DMF8Q9AFzL7r1rDpGq9UorspdWbpe0vJVnp9oTp8GbM4jLzH2bKoi01KwuBXFT9sRJbVPeO9WeOA3enK0sIn6pnn0qaYaxQk6c225J7VPKhS3lzEr9Nl+A9k/AHRHxqYh4WUTs3UShpKWgOJzqxAapdVwnTsMyXyu7Gvgp4BXAZ4C7IuIfI2JDRDxh4KWTRlgp58UdG6T2cXaqhmXOIC4zfw5YC/w+MAGsBE4Czga2RMTl9R6pRw62mNLoKV1pu06c1D6lvVPtiVMT5m1lmXlrZv5VZj4f+GngVcDngQeBY6n2SL02Im6MiPdFxM9HhJcgWvZKV9qlpQYkLW+lpYVKObNSv3V1qZCZ92Xm/8zMlwMHUvXKfRy4HfgZ4C3APwO3R8THI+LkiNiz34WWRkFpRXaHU6V22bEzmbmyUASsNIhTAxbd35uZ2zPzS5n525l5MHA88B7ge8Djgd8C/p4qj+7vI+L0iHhcX0otjYBiT5zDqVKrzLbllgNSasJ8S4w8f6FPlJn/mpl/nJlPB54CvBX4BrAncArwCeB3eyirNFLKOXH2xEltMlVY4NtJDWrKfD1xX4mIz0XE4d08aWZuzswPZOYvAE8AXgf8A1UenbQslGen2hMntYn5cBqmhQyn/hqwKSLetZj8tsy8KzM/mZkvzcz/2n0RpdE0VsqJc9stqVVKve+rx5yZqmbM19LeANwN7AX8KfC9iHjZwEslLQGlIZPSUgOSlq/illuFCzxpEOZbJ+5sYB3wQWAKOAz4TERcEhFPH3zxpNFV2jvV4VSpXVzoV8O0kHXi7s/MNwFHAV8AAhgHroqID0fE/gMtoTSiSttuucSI1C6ltSFd6FdNWXBLy8x/y8yXAC8ErgfGqIZbb6y34fLSQ61SXuzXIE5qk9KFW2nSkzQIXV8uZOZXgGdQLRdyD9Wiv2cDV0bEc/tbPGl0lRb7nZxyOFVqE3PiNEyLammZuTMzz6LKl/swVb7cMcDXI+JvIuKJfSyjNJKKS4y47ZbUKsXFfu2JU0N6ulyot+H6A+Bo4ItU+XKnUc1ifVtErOq9iNJoKuXEudiv1C7lxX7tiVMz+tLSMvPGzPxV4EVU227tA7yban25F/fjNaRRY0+cJHviNEx9vVzIzIvZlS93L/Bk4MJ+voY0Kkp5L5NT9sRJbVLqfXd2qpoy1usTRMQKqr1Sn061DMlR9e3HTp/S62tIo2j1WGl2qj1xUpu47ZaGqasgLiKexO7B2lOB1Z2nzfixexdbuIg4EfgQsBL4eGa+Z5bzfgP4HPDvMvPKxb6e1I3itlvmxEmtUuqJMydOTZkziIuIN7MrWDsSeMzMUzpuTwI3AFcD10x/z8wtiylYRKwEzgJ+GdgCXBERGzNz04zz9gV+H/jmYl5HWqzSOnHu2CC1SykPtrSbizQI8/XEvR9Idu9d28quQG06aLs+M6f6WLZjgc2ZeRNARJwPnAJsmnHenwPvBd7ax9eW5lXKe9luT5zUKsV14pzYoIbMF8Q9AFzL7r1rix4i7cLBwG0d97cAx3WeEBHPAg7NzC9EhEGcGlUK4pydKrVLcTjVxX7VkDmDuMzct6mCdKueUPHfgNMXcO4GYAPAmjVrmJiYGGjZALZt29bI62jh+l0nD07u/uH98COT1nuXfK+MHutk4TbdNrnbsTvvuJ2Jif72dVgno2nY9dLz7NQB2goc2nH/kPrYtH2p8vUm6m1bnwBsjIiTZ05uyMxzgHMA1q9fn+Pj4wMsdmViYoImXkcL1+86eWj7DrjkS486lrHCeu+S75XRY50s3K3/cgtcd92jjh1+yMGMjx/V19exTkbTsOtllPt8rwDWRcQREbEaOBXYOP1gZt6fmQdm5trMXAtcDuwWwEmDUlrQs5QfI2n52l7YL9mcODVlZIO4epLEG4GLgeuBCzLzuog4MyJOHm7pJFhZWAtqZ8LOwjY8kpan0rZbLvarpozycCqZeRFw0Yxj75jl3PEmyiRNiwhWrYzdEpsnd+5kjxUrh1QqSU1ysV8Nk5cLUg9c8FdqN7fd0jDZ0qQemBcntVvp/V76XJAGwSBO6kHpirt0ZS5peSrlxLntlppiS5N6UNx6ywV/pdYo7thgTpwaYhAn9cCcOKndSu/31WP+a1UzbGlSD0of1tvNiZNao9wT579WNcOWJvWgNGxiT5zUHsW9U53YoIYYxEk9KCUwOztVao9SDqyzU9UUgzipB6UP69JsNUnLU6nn3XXi1BRbmtSD8hIj9sRJbVHKgTUnTk2xpUk9KOXEGcRJ7VHadsvhVDXFIE7qQaknzokNUnu42K+GyZYm9aA0C82eOKk9tk/ZE6fhMYiTeuC2W1K7lXrinNigptjSpB6UZ6faEye1RSknzm231BSDOKkHbrsltVup592eODXFlib1oJQT57ZbUnuUcmAN4tQUW5rUg9XOTpVarTw71eFUNcMgTupB6cPanDipPYo9cS72q4bY0qQelHLinJ0qtUep592eODXFIE7qQWl2quvESe1hTpyGyZYm9aC8Y4NBnNQW5SDOnjg1wyBO6kFpex2HU6X2cNstDZMtTerBqsKink5skNqjmBPnYr9qiEGc1AN74qT2ysziupDmxKkptjSpB05skNprR2EodUXASnvi1BCDOKkH5YkN9sRJbWA+nIbN1ib1wMV+pfYqL/RrL5yaYxAn9aC0Mvv2KXvipDYo5b+uGvPfqppja5N6sGrMnjiprUprQpZ2cZEGxdYm9aD0gW1OnNQOk4WcOBf6VZMM4qQeODtVaq9iT5xBnBpkECf1oNQTZxAntYP7pmrYbG1SD0pJzKVlByQtP8WJDebEqUG2NqkHpeUE7ImT2qG45ZbDqWqQQZzUg9LCnk5skNphsjAT3cV+1SRbm9QDJzZI7TU5tft7fbU9cWqQQZzUg1IScylPRtLyU9x2y5w4NcjWJvXAbbek9ir1upsTpyYZxEk9cLFfqb1K7/XV5sSpQSPd2iLixIi4ISI2R8TbCo+/KSI2RcTVEXFJRBw+jHKqvUof2NvNiZNawZ44DdvIBnERsRI4C3gR8DTgtIh42ozTrgLWZ+YzgM8B72u2lGq74nCqPXFSK5S23XJ2qpo0yq3tWGBzZt6UmduB84FTOk/IzEsz88H67uXAIQ2XUS1nTpzUXqVtt0prR0qDMspB3MHAbR33t9THZnMG8MWBlkiaobQ6u7NTpXYo9bq77ZaaNDbsAvRDRLwSWA88b5bHNwAbANasWcPExMTAy7Rt27ZGXkcLN4g6eWhq9w/xh7dPWvdd8L0yeqyThbn21sndjt35o9uZmLi3769lnYymYdfLKAdxW4FDO+4fUh97lIh4AfB24HmZ+UjpiTLzHOAcgPXr1+f4+HjfCzvTxMQETbyOFm4QdfLw5A74py896liywrrvgu+V0WOdLMwtl90MmzY96thhhx7M+PhRfX8t62Q0DbteRrnf9wpgXUQcERGrgVOBjZ0nRMQxwMeAkzPzjiGUUS1XXOzXnDipFVzsV8M2sq0tM6eANwIXA9cDF2TmdRFxZkScXJ/2fmAf4LMR8Z2I2DjL00kDsXJFEDPymDNhR+HDXdLyUsp/XTXmxAY1Z5SHU8nMi4CLZhx7R8ftFzReKGmGVStXsH3GHoqTO3aycsXKIZVIUhNK68SVJjtJg2Jrk3pUWlKg9OEuaXkpLTHiYr9qkkGc1KPS4p4u+Cstf6XFfl1iRE2ytUk9WlW48rYnTlr+JqcKw6n2xKlBBnFSj8ozVO2Jk5Y7Z6dq2GxtUo/K+6faEyctd8WJDfbEqUEGcVKP3HpLaqdS7mspR1YaFFub1KNST5w5cdLyV+6J89+qmmNrk3pU+tB2dqq0/JVnpzqcquYYxEk9Kg2fuPWWtPwV14lzYoMaZGuTelRa7NeeOGn5K+W+utivmmQQJ/XInDipnUrv89XmxKlBtjapR8V14gzipGVvqpA2YU+cmmQQJ/XIiQ1SOxWHU82JU4NsbVKPxko5cU5skJa90sQGZ6eqSQZxUo9KPXHb7YmTlr1ST5zrxKlJtjapR6Urb7fdkpa/Uu6rOXFqkkGc1KPSOnHmxEnL31RxsV//rao5tjapR6WeOBf7lZa/8mK/9sSpOQZxUo9Ks9EmpwzipOXOnDgNm61N6lFxiZHCMIuk5aWUE2cQpybZ2qQeFYdTzYmTlr3SxZoTG9QkgzipR6UPbWenSstfsSfOxX7VIFub1CO33ZLaqRjEjdkTp+YYxEk9KgZx5sRJy15pKSG33VKTbG1Sj4rbbtkTJy1rmTnLOnH2xKk5BnFSj0qL/TqxQVreSgHcyhVBhEGcmmMQJ/VodXF2qj1x0nJWXl7EAE7NMoiTeuS2W1L7FBf6NR9ODbPFST0q5cS57Za0vBW33LInTg0ziJN6VNyxwZ44aVkrL/Trv1Q1yxYn9ch14qT22V7YH3m1QZwaZouTelQaQnF2qrS8ueWWRoFBnNSj0oy0KXPipGWtmBNXyI+VBskgTupRaYV2c+Kk5a04O9XhVDXMFif1qPTBvd2cOGlZK68T579UNcsWJ/WoOJxqECcta6WUCXPi1DSDOKlHxcV+C0nPkpYPF/vVKLDFST0qJTOXlh+QtHyUhlPtiVPTDOKkHq0esydOapvS5CVz4tQ0W5zUo1JPnDlx0vJWnthgT5yaZRAn9ai8Y4M9cdJyVlzs15w4NWxs2AWYS0ScCHwIWAl8PDPfM+PxPYDzgGcDdwP/PjNvabqcardSHsz2HTv5ycOTQyjN0vPQVPq3GjHWyfxKfx9z4tS0kQ3iImIlcBbwy8AW4IqI2JiZmzpOOwO4NzN/JiJOBd4L/PvmS6s2K/XE3fmTRzj6XV8eQmmWqH/ybzVyrJOuuXeqmjbKLe5YYHNm3pSZ24HzgVNmnHMK8Nf17c8Bz48IL4XUKJcVkAT2xKl5I9sTBxwM3NZxfwtw3GznZOZURNwPHADc1XlSRGwANgCsWbOGiYmJARV5l23btjXyOlq4QdXJ9h3JigAnpErtdu+dtzMxce9Antv/KaNp2PUyykFc32TmOcA5AOvXr8/x8fGBv+bExARNvI4WbpB18tybv8n/2nzX/CdKWrZe8UvHMP7Unx7Ic/s/ZTQNu15GOYjbChzacf+Q+ljpnC0RMQbsRzXBQWrUB0/9Of70wmv55s33uNBvl6amphgbG+WPovaxTrrzhP325DXPOXxgAZw0m1F+l14BrIuII6iCtVOB/zDjnI3Aa4B/AV4GfDUzHdRS4w7cZw/OfuWzh12MJWnYV7LanXUiLQ0jG8TVOW5vBC6mWmLkk5l5XUScCVyZmRuBTwCfiojNwD1UgZ4kSdKyN7JBHEBmXgRcNOPYOzpuPwy8vOlySZIkDZtrI0iSJC1BBnGSJElLkEGcJEnSEmQQJ0mStAQZxEmSJC1BBnGSJElLkEGcJEnSEhRt2+AgIu4Evt/ASx0IuJnmaLFORpP1Mnqsk9FjnYymJurl8Mx8fOmB1gVxTYmIKzNz/bDLoV2sk9FkvYwe62T0WCejadj14nCqJEnSEmQQJ0mStAQZxA3OOcMugHZjnYwm62X0WCejxzoZTUOtF3PiJEmSliB74iRJkpYgg7geRMSJEXFDRGyOiLcVHt8jIj5TP/7NiFg7hGK2zgLq5U0RsSkiro6ISyLi8GGUs03mq5OO834jIjIinIXXgIXUS0T8Zv1+uS4i/rbpMrbNAj6/DouISyPiqvoz7KRhlLNNIuKTEXFHRFw7y+MRER+u6+zqiHhWU2UziFukiFgJnAW8CHgacFpEPG3GaWcA92bmzwB/Cby32VK2zwLr5SpgfWY+A/gc8L5mS9kuC6wTImJf4PeBbzZbwnZaSL1ExDrgj4ATMvPpwB80Xc42WeB75U+ACzLzGOBU4CPNlrKVzgVOnOPxFwHr6q8NwNkNlAkwiOvFscDmzLwpM7cD5wOnzDjnFOCv69ufA54fEdFgGdto3nrJzEsz88H67uXAIQ2XsW0W8l4B+HOqC52Hmyxciy2kXl4PnJWZ9wJk5h0Nl7FtFlInCTy2vr0f8IMGy9dKmfk14J45TjkFOC8rlwP7R8RBTZTNIG7xDgZu67i/pT5WPCczp4D7gQMaKV17LaReOp0BfHGgJdK8dVIPPxyamV9osmAtt5D3ylOAp0TEZRFxeUTM1Ruh3i2kTt4FvDIitgAXAb/bTNE0h27/7/TNWBMvIo2iiHglsB543rDL0mYRsQL4b8DpQy6KdjdGNUQ0TtVj/bWIODoz7xtmoVruNODczPxARDwH+FREHJWZO4ddMDXPnrjF2woc2nH/kPpY8ZyIGKPq+r67kdK110LqhYh4AfB24OTMfKShsrXVfHWyL3AUMBERtwDHAxud3DBwC3mvbAE2ZuZkZt4M3EgV1GkwFlInZwAXAGTmvwB7Uu3fqeFZ0P+dQTCIW7wrgHURcURErKZKMN0445yNwGvq2y8DvpouzDdo89ZLRBwDfIwqgDPHZ/DmrJPMvD8zD8zMtZm5lipP8eTMvHI4xW2NhXyGXUjVC0dEHEg1vHpTg2Vsm4XUya3A8wEi4kiqIO7ORkupmTYCr65nqR4P3J+ZP2zihR1OXaTMnIqINwIXAyuBT2bmdRFxJnBlZm4EPkHV1b2ZKiny1OGVuB0WWC/vB/YBPlvPM7k1M08eWqGXuQXWiRq2wHq5GPiViNgE7ADempmOJgzIAuvkzcB/j4j/TDXJ4XQ7BwYrIj5NdTFzYJ2L+E5gFUBmfpQqN/EkYDPwIPDaxspm3UuSJC09DqdKkiQtQQZxkiRJS5BBnCRJ0hJkECdJkrQEGcRJkiQtQQZxkiRJS5BBnCT1WUR8JyKy3hmk8/jW+viHh1U2ScuHQZwk9VFEHA48E7gPmOg4vgZ4Yn33240XTNKyYxAnSf11Sv39osyc6jj+7I7bBnGSemYQJ0n9NR3E/f8zjj+r/v4wsKm54khargziJAmIiBPqfLWMiN+c5ZzjImJbfc77C48/Dvh/gO3AF2c8PB3EXTOjh06SFsUgTpKAzLwM2Fjf/bOIWNn5eEQ8FfgCsDfw18AfFp7mV4Ex4KuZ+ZMZj00PpzqUKqkvDOIkaZc/AnYAPwu8cvpgRDwRuBg4APhH4HWZmYWfLw6lRsQBwGH13W/1ucySWsogTpJqmbkJOLe++86IWBUR+wNfAg4H/hfwm6Xh0IjYAzgRSHb16E17Vsdte+Ik9YVBnCQ92juBh4AjgDdQ9aodDVwDvCQzH5rl534J2Ae4MjN/MOOx6SBusn4eSeqZQZwkdcjMrcD0Yrx/STVR4RbghZl53xw/+tL6+8xZqbAriLsuM7f3XkpJMoiTpJIPAzvr2/cAv5KZP5zt5IgI4CX13QsLpzipQVLfGcRJUoeIGAM+xq7Px8dQDa/O5VjgIOB/Z+Z1M55vP+BJ9V2DOEl9YxAnSbW6R+3jwIuBO4GbgT2BP5vnR2db4BfgGCDq2wZxkvrGIE6Sdnkf8BpgG9Wab2+vj78mIp42x8+9tP4+Vz7cDuC7fSijJAEGcZIEQES8BXgL1QzSX8/MK4DzgauBlcBfzPJz64AjgbuAywqnTOfDfS8zH+x3uSW1l0GcpNaLiFdT9cIlcHpmfgWgXtD3T+vTTo6IEwo/Pj2U+o+ZuaPw+P/dbisi9pnna6yPv5akZS7Ki45LUjtExElUw6BjwH/OzA8WzrkcOA64LDN/fsZjXwd+Hvi1zLxwxmN7Az9m4RfMuz2HJM3GnjhJrRURzwE+SxXAvbcUwNWmc+NOiIjpnjci4vHAc6lmr3658HM/R3efs9/p4lxJLWdPnCQtUkT8FvAJ4B8y8+Rhl0dSu9gTJ0mLN90rd+EwCyGpnUyilaTFuwy4ivLSIpI0UA6nSpIkLUEOp0qSJC1BBnGSJElLkEGcJEnSEmQQJ0mStAQZxEmSJC1BBnGSJElLkEGcJEnSEmQQJ0mStAT9H1+aOTKM3YvIAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(10,5))\n", "plt.plot(y, mL2V(y), lw=5)\n", "plt.title('Potential', fontsize=30)\n", "plt.ylabel('$V/V_0$', fontsize=25)\n", "plt.xlabel('$x/L$', fontsize=25)\n", "plt.grid()\n", "plt.savefig('v3p1.png', dpi=200)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAn4AAAFoCAYAAADJt/Z7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACe9ElEQVR4nOzdd3xT19348c/RsOSBJx4YMGZ4MGwwNisJhOzVZq9mNGmbpvNpn6a7fZ4mv+6dPm3TtGnTNmmSZu/ZhMQEAmba2IAnYMDGA+OBl2RLOr8/rmQMeMj21bLP+/XSC/nec8/9WsLW12cKKSWKoiiKoijK5GcIdACKoiiKoiiKf6jET1EURVEUZYpQiZ+iKIqiKMoUoRI/RVEURVGUKUIlfoqiKIqiKFOESvwURVEURVGmCJX4KYriV0KIdCGEdD/+Geh4FEVRphJToANQFCU0CCHGs+jnHinlMr1jmeqEEOuB9e4v/ymlrPXTfZcB17q/fFlKWeKP+yqKoh+V+CmKooSe9cD97ueFQK2f7rts0H1rgRI/3VdRFJ2oxE9RlPG4zstyHWcecLdOCV2jURRFUbyiEj9FUcZMSvlyoGNQFEVRxk5N7lAURVEURZkiVOKnKIpfjWVWrxAiWQjxSyHEfiFEtxCiVQixXQjxdSFEhLtMrbuuWi/uvVgI8VshRIm7LrsQol4I8aoQ4nYhxLC/E4eKWwgxXQjxgBCiTAjR6X7sFkJ81xPfKPFcI4R4VghxUAjRI4SwuePZI4R4TgjxRSFEwqDyD7gn2dw/qJoPBsXleRQOca98IcT/CiHeFkIccd+rVwhxVAjxshDiDiGEcZg473bf9x+DDv9jiPvWjvC9jvu1H+/rpSjKEKSU6qEe6qEeoz4A6XlMsJ70QXX9c4RyFwCtg+97xqMUmI02yUACtSPUZQL+D3COUJ8EioAUb+IGCoC6EeoqBuKHqSsceH2UWDyP/x503QNeXlN4xv3u9/K6nUDqEPHe7eX1Z70HOr3243q91EM91OPshxrjpyhK0BFCLEL7oPe0mhUDT6AlWsnALcC5wDOMMlZZCCGAZzk1IaUBeBrYA/QAc4BbgXxgFbBBCLFCStkzQrWzgTeAeOBJ4AOgC1gEfAlIQJsB+zvgk0Nc/1PgqkHxPAHsc9cRBSwA1gDrzrjuabSZtLeivQYA/wvsPaNcyxlfhwMOYCvwEVADnHTHPxe4A5iJ9hq8LIQ4V0rZP+j699FevwuB/3If+4P7+GCnvWY6vvbjfb0URTlToDNP9VAP9QiNB35s8QM2DSrzEGAYoswDjNLa5C731UFl/gVEDFFGAD8ZVO7no8QtgTZg1RDl5rrPSbRkK/WM80ag3RMzkDTCa5UIZI/yva/34jVfwTCtae7zYWhJqqfOu4Ypd/egMnd7cd8Jv/Z6vF7qoR7qceqhxvgpijJmQ4ztGu5x9zjqLgDOc3+5F/iKlNJ1Zjkp5QPAxlHqsgLfc3+5Ay2hOaslT2q+j5ZwAnzBfe1IviKl3DZEXYfQklXQkpaLziiSCMS4n78ipWwe7gZSyuNSyopR4hiVlHKHlLJxhPN9wNeBQ+5Dd070njq+9n5/vRRlMlOJn6IoweaaQc8fklI6Ryj7f6PUdRmQ5H7+4FAJ5BmecP8bDaweodxx4KkRzg/uAl10xrneEc4FjPt19iSyK93dtBOh12sflK+XooQqNcZPUZTx8HYB593jqLtg0PMPRilbOMr5tYOexwkhrh2l/MxBzxeOUP/OURLS+sH3HXxCStkhhNgOrAQuFkK8hDZebpM8fVydrtyzZq8FbgDygFRgGkM3AExDS8DOWoB7DHR57QP1einKZKUSP0VRxkz6dgHn1EHPD44SR5sQoh2IHaZI+qDnDw1TZjhxI5w7c/LEmeyDng/VZfwlYANacnWt+9EthNgGbAbeAz7yopXMK0KIWcDLaJMovDXRxC990POJvvZ+fb0UZTJTXb2KogSbSPe/Di9bdLpHOBczwrnRhI1wbkIJhpRyJ9qs38c51ZUZiTZr9gfAh8ABIcTtE7kPgBDCDLzDqaSvBfg7cB9wO1oL4HXux+AW1iHX9BsD3V57f75eijLZqRY/RVGCjSeRMwkhzF4kf5EjnOsa9Hyee+JFUHDHcpcQ4nNoS5Gcgzap5Xy05VfSgSeEEHOklD+dwK0+wamxce8C10kph0yWdU6cdH3t/fh6Kcqkplr8FEUJNscGPZ83UkEhRBzDd/PC6WPtZk0gJp+RUtqklB9IKX8ipbwCbULEt9GWLwH4wQR3o7h40POvDZf0uc2ZwH3O5JPX3g+vl6JMairxUxQl2Owc9PyCUcquH+X84OVeLh1XNH4mpeySUv4SeMF9yIK2Dt9gg7uaR5t9mzzo+YHhCgkhktC6U0cylvv65bX38vVSFMVNJX6KogSbVwY9/9Jw+8e6fXWUut7k1ESMLwohZkwoMv+qHfT8zGE5g7tRR+rqhtN305g/QrnvAuZR6hrLff392tcOeq6GMSnKMFTipyhKUHEP5N/s/nIJ8Hv3UiSnEUI8gDa+a6S6uoH/5/4yHnhbCJEx0jVCiFVCiF+ONW5vCSHyhBD/K4RIHqHMdOAm95cSbV/iwQaPl1s+yi13DHr+o2Fey3uBr4xSz5juq9drr9PrpSiKm5BSjl5KUZQpTwgx+JeFt+v4AbwxeIKGECKdUwnEY1LKu4e41yK0hGW0vXq3Amloa8AdklIOOSZQCPEYp/bMdQCvos0EbUCbvZoI5KDtsjEXOCClXHBGHaPG7U1ZIcR6tNmzTrR9c7cAVUAnWoKUA9zmfg7whJTytJ003InOMbQWui7gl2jJjmcZmVYp5XZ32RlANada6PagbZ/meS2vR0ugG4Ey4BJ3ublSytoz7mtCG7uXBPQDDwJFnJpp2yul3HjGNRN67fV4vRRFGSTQe8aph3qoR2g8OH2f2rE8Ys+oJ33QuX+OcL8LOLXn7VCPMmA2WgIjgT0j1CWA/wFsXsZcOEQdXsU9Wlm0JMvb1+5pIHyYe/zU2/jRdkPpHaF8Hdq4uH8OOpY+zH3vHaGeWr1fe71eL/VQD/XQHmochKIoQUlK+YEQIhv4BvBxtJY9O1CD9gH/MFoy42npaR2hLgn8WAjxKHAPWutSlvtaF9pYtAq0FsQ3pZRFvvie3LFsFELkoG1ptgZYjDbrNQJtPN4RtFa0x6SUH45Qz/eEEKXAXWiTMuIZZu1BKeUrQojlwLfQvvcU4CTauLhXgD9JKU94s0ublPIRIcRh4PNou6wkok2oGK78hF57vV4vRVE0qqtXUZSQ5U4IPOO5fi+lHG2yh6IoypSmJncoihLKvjzo+QcBi0JRFCVEqMRPUZSgJIRYO9QM1EHnv4Q23gy0CQev+yUwRVGUEKa6ehVFCUpCiBrACryFNqv3ONos1vlos4rzBhX/uJRSJX6KoiijUImfoihByZ34jbTgMGiTOz4rpXzSDyEpiqKEPJX4KYoSlIQQq4EbgNVo6/QloM3kbENbx+094GEpZVPAglQURQkxKvHzwvTp02V6erpu9XV3dxMZOdpuR0qgqPcneKn3Jrip9yd4qfcmuOn9/uzatatFSpk41Dm1jp8X0tPT2blz5+gFvVRYWMj69et1q0/Rl3p/gpd6b4Kben+Cl3pvgpve7497rc0hqVm9iqIoiqIoU4RK/BRFURRFUaYIlfgpiqIoiqJMESrxUxRFURRFmSJCMvETQlwuhKgUQtQIIb4zxPn7hBD7hRClQogNQog5g87dJYSodj/u8m/kiqIoiqIogRNyiZ8Qwgg8BFwBLAI+IYRYdEaxYqBASpkLPA/80n1tPHA/sApYCdwvhIjzV+yKoiiKoiiBFHKJH1rCViOlPCil7AOeBq4ZXEBK+YGUssf9ZREwy/38MuBdKWWrlLINeBe43E9xK4qiKIqiBFQoJn4zgaODvq5zHxvOZ9D2+hzPtYqiKIqiKJPGpF7AWQhxB1AAnD+Oa+8F7gVITk6msLBQt7i6urp0rU/Rl3p/gpd6b4Kben+Cl3pvgps/359QTPzqgdmDvp7lPnYaIcTFwPeB86WU9kHXrj/j2sKhbiKlfAR4BKCgoEDquaK2WkE9uKn3xwdsHdBcAW21EDMLkhdB+NiH14bie2Prd1Lf3kunzUFPn4PePicnbTb2t+0iJiyB+TGZhIcZiQgzERFmJD4yjORoK0aDCHToYxaK789Uod6b4ObP9ycUE78dQIYQYi5aIncrcNvgAkKIPOAvwOVSyuZBp94BfjpoQselwHd9H7KiTDHHK6HkKWjeD0374WTd2WWmzYCkRVoSmHMzzMj1WTjS4aD9uedwnGg97bgwGoi5+mrMMyc+4uN4p52So+3sre/gSGsPR1p7ONraQ3OnfaCMIawZc+wOTDG7MZi6AXD2zqS/fQX9J5eBywqA2SiYFRfBrLhwZsdHMG96JHlpsSxOjcFqNk44VuWUvtpaOt54A87Ytt6UnETsDTcgDKE4IkpRhhdyiZ+U0iGE+DJaEmcE/i6l3CeE+CGwU0r5KvArIAp4TggBcERKebWUslUI8SO05BHgh1LK1iFuoyjKeLQego2/gNJnQBghMQvmnANJC7UkL34utB/REsLmcmjaB9segS1/gEXXwgXfh8RM3cM6/sc/cuLPfxnyXMfrbzD3+ecwhId7XZ+UkurmLjZXt1B8tJ3iI23UtfUCYBAwIyac2fHhnJ+ZyOy4cNoM29nd/gYHO/diFEZWJK3lsjkfo6m7kbcOv8zh8JeZNvMtlsatI3faNdh7kjnapiWOe8saaOvpB8BkECxKjWbZ7Fjy58SxLiORuMiwib9AU5Tz5EmOfPoz9B87NuR518lOEj7zaT9HpSi+FXKJH4CU8k3gzTOO/WDQ84tHuPbvwN99F52iTEEnj8HGX0Lxv8BggjVfgnO/BpEJZ5dNzIKMS0593dsOW/8IRQ9D+auQeyus/zbEpesSWnfRNk785RFibrieGT/+8WnnerZu5chn7qHpZz9nxg//34j19DtdbD/UynvlTWwob+ZIq7ZwQGqMlWVpsdy1Jp1labEsSY0hPOxUq9xj+x7jrzt/TXp0Ovfl38fH53+c6eHTB85/Mf+T7D+xnxeqX+DNQ29SebKIf135L+bHLh8o09xpo+RIO8VH2yk50s7zu+p4fOthDAIK5sRz0cIkLl6UzPzEKD1esilBSknDD+6nv7mZ9GeexpqbO/gk9f/9NZoffJCIlSsIz8kJXKCKorOQTPwURQkSUmqtde//GKQLlt8F674J0TO8ryM8Fi78H1j1edj8IGz/K5Q9B2vvg/O/AxPoanO0tXHsW98iLD2dlO9/H3cPwIDIc84h4Z7PcOKvfyPynHOIvvyyM749yY7aNp7deZR39jXSaXNgMRk4d8F0Pnf+PC7ISiI1dviWwncPv8tvdv6GS+dcyq/O/xUGcfb3IoRg8fTFLJ6+mHty7uG2N27jSxu+xBNXPjGQICZNs3Lp4hQuXZwCgNMl2VvfwYbyJt4rb+Znb1Xws7cqmDc9kmvzZnJj/qwR41Kg44UX6Hz7bRLvu4/wpUtPPykEM370Q3rLyqj/+jeY++KLGKMiAxOoouhMJX6KooyPsx9e/5rWypf9MbjsJxNrpYucrtWx5kvw3gNal3FLNVz7JzCPPYmRUtLwve/jbGtj9p8fxhARMWS5xK98he5t22n4wQ8Iz1mCeeZMmjttvLCrnud2HuVgSzdRFhNXLEnhkkXJnJcxnYiw0X91lh4v5bubvktuYi4/Oe8nQyZ9Z0qNSuWPF/2RT739Kb7y/ld49LJHCTed/b0bDYKls2NZOjuW+y7Nor69l/fLm3izrJHfvlvF796rYm1GIreumM1FC5MJM6lxaoPZDxyg8Sc/JWLNahLu+cyQZYwxMcz81S85/Mm7aPrRD0n9xS/8HKWi+IZK/BRFGbvednj2k3BoI6z9hjY2T69B8NGpcN1fIHkxvPsD6KiDT/xbSwzHoO3Jp+j64AOSv/ddrIvO3NznFGE2M/PXv+LQdddT9ZX7+NO1X+ft8hacLsmK9Di+sH4+V+XO8CrZ86jrrOO/3v8vEsMT+f2Fv8dqsnp97ZLpS/jFul/w3x/8N9/b9D1+s/43oyaNM2PDuXNNOneuSefIiR6e23WU53fV8YUnd5MQGcbtq9K465x0EqIsXscxWbnsduq//g0MViupP//FiJM3IgoKmP7FL9Lyxz8See65xFx9tR8jVRTfUH8GKooyNm218PfL4PBHcM1DcNH/6pf0eQgB534VbnoMGkvhbxfB8SqvL7dVVtL8y18Sdf75xN1554hlpZRst0fwwvrbMewrJfnlp7jnvLls+Pr5PPf5c7ipYPaYkr6TfSf50oYv4XA5eOjih4i3xnt9rceFaRfyzRXf5L0j7/HgrgfHdG1aQgRfvzSLzd++kH98agXL58Tx+/drOPcX7/PAq/uob+8dczyTSfNvfoO9ooIZP/sp5uSkUctP//znCC/Ip/GB/0ffkSN+iFBRfEslfoqieK9uF/ztYuhsgDtfgrw7fHu/xdfC3W9AXzc8ejGxbaWjXuLq7aX+vq9jiIlmxs9+eta4Pg8pJW/vbeSahz7ijke38XJ8Dk1rLuKm8nf5amLXuCZK9Lv6ue+D+zjSeYTfXfA75sXMG3MdHncsvINbs27ln/v+ybOVz475eqNBcEFWEn/9ZAHv3beOj+em8kTRYc7/5Qfc92wJB453jTu2UNVZWEjb4/8i7pN3Ms3LNdOEycTMX/4SzGbq7/s6sq/Pt0Eqio+pxE9RFO+01MAT14M5Aj7zHsxd55/7ziqAezZAVAo5ZT+Ghj0jFm9/7jn6Dhwg9ec/xxQ/dGvbrsOt3PDwFj7/xC5O9vbzs+tz2PTtC1j7h19gTk2l6Ve/Hleo/y7/N9sat/HAmgdYkbJiXHV4CCH49spvs3bmWn6+/ecc6xp6yRFvLEiaxq9uWsqH37qAT65J562yRi598EP+9+W9tHTZR69gEpBS0vyLX2LJyCDpG98Y07Xm1FRm/L8HsO3dS8cbb45+gaIEMZX4KYoyup5WeOombamWu171yVp7I4qbA3e9Sr95Gjx1K5xsGLZo+0svY128mKhzzz3r3KGWbr7wxC5ueHgrdW29/OKGHDZ8fT2fWJmGxWTEGBVJ3B23YystxV5TM6YQu/q6+GvZX1k9YzXXLLhmzN/iUEwGE/+7+n8RCB7e8/CE60uNDecHH1/Epm9fwO2r0nhq+xHW/6qQP75fTW+fU4eIg1dvSQl9hw4Rf/fdGMLGvvbhtMsuw5yWRsdLL/kgOkXxH5X4KYoyMkefNpGjow5ufUq39fXGbFoKZTn/o23/9u9boa/nrCK2igrs5eXEXHfdacc7evt54NV9XPLbjWysOs59l2RS+M313LIi7ayt0WI+/nEwmeh4+eUxhff4/sdpt7fz1eVfHfO3NpIZUTO4JfsWXj3wKgfbD+pS5/QoCz+8Zgn/+do6zpmfwK//U8UFvy7kxd11SClHryAEdbz0MiI8nGmXXTZ64SEIIYi97lp6tm+nr+6sXUIVJWSoxE9RlOFJCW98DWo3aRM50lYFNJzuqLlw46Nad+9L94LLddr5jpdeBrOZ6KuuBLTuvTdKG7j4txt5fGstN6+YTeE31/OVizKGnbBhSkggau1aOl59Den0rhWs1dbKY/se45I5l7Bk+pIJfY9DuSfnHsJN4fyx5I+61js/MYpHPlnAs59bQ3K0hfue3cPtf9vGoZZuXe8TaC6bjZNvvUX0pZdMaD2+mGuuASHoeOVl/YJTFD9TiZ+iKMPb8nsofgLWfQtybw50NJqsK7T1/spfg/d/NHBY9vfT8dprTFu/HlNcHPXtvdzz2E6+9NRukqZZeOVL5/HT63JImjb60iox112Lo7mZ7i1bvArpr6V/xea08eVlXx73tzWSeGs8dy26i3cPv8velr26179ybjwvffFcfnztEsrqOrjsdx/y0Ac19Dlco18cAjo3bMDV2XlWS/BYmVNTiVi1io6XXka6Jsdro0w9KvFTFGVo5a/Du/fD4uvhgu8FOprTrf4i5H8KNv8WSp4CoGvTJpytrUy79loe3XyIS367kS0HTvA/Vy3klS+dS86sGK+rn7Z+PcbYWK/GczV0NfBM5TNcPf9q5sWOfxbvaD65+JPEWeL4v93/55P6DQbBHavn8N7Xz+fihUn86p1KPv6Hzew63OaT+/lTx0svY0qdQcTKlROuK/a6a+mvq6N31y4dIlMU/1OJn6IoZ2s7DC99DmYu13bOGGZJlIARAq78FcxbD699FZrLtTFccfF8rsLMj17fz6q58fzna+u4Z+08TMax/aoTYWFEX3UVne9twNnRMWJZz6SLLy794ni/G69EmiP5bO5nKWoooqihyGf3SY628qfb8/nbJwvotPVz05+38Kt3Kuh3hmYLV39TE91bthBzzTUjLtbsrWmXXIIhIoL2MY4BVZRgoRI/RVFO53LBq18GBNz0z3Ftl+YXRjPc8ChYonH8+3OcfP8DXk3MZW9zD7+5aSl/v3sFs+OH3qbNGzHXXYfs6+PkW28NW+Zg+0FeOfAKt2TdwoyoMexPPE43Z91MSmQKv9/9e59Pwrh4UTL/ue98bsyfxUMfHODGh7dwMATX/ut45VVwuYi99lpd6jNERDDtisvpfOttXD1nTzBSlGCnEj9FUU6381E49CFc9mOITQt0NCOLnE73Jb/k5LYahNPB4YL1vPXVtdyQP2vYhZu9ZV28CEtGhjZhZBh/LPkjVqOVz+Z+dkL38pbFaOGLS79IWUsZ7x953+f3i7KY+OWNS3n49uXUnujhqt9v5t/bj4TMzF8pJR0vv0x4fj5hc+boVm/sddfh6umh8913datTUfxFJX6KopzSekgb1zf/Qlh+V6CjGdWuw61c+k4cRw4mYonr5//uWTShVr7BhBDEXHstvXv2YD949jIq+1r28e7hd7lr8V3j2pZtvD4+/+PMjZnL74t/j9Pln7X3rsiZwTv/vY7lc2L57otlfO5fu+jo6ffLvSfCVlpK38GDxFyrz7qKHuH5+Zhnz6Z9hD8KFCVYqcRPURSNywWvfBkMRrj6D8E3rm8QKSX//OgQt/yliNntx7C09xGbZcD06hfBqV9CEv3xj4HROGSr37/K/8U08zQ+ueiTut3PGyaDiS8u/SIHOw6y5Zh3s471kBJj5V+fXsX/XLWQDyqb+fgfN7Pv2MjjHwOt/aWXEFYr0Zdfrmu92h8F19CzbRv99WpNPyW0qMRPURTNjr/C4c1w2U8hZlagoxlWT5+Drz1TwgOv7Wd9ViK/nHZUW7vvCz+ExlLY/KBu9zInJRF13nl0vPLKaWv6ddg7eLf2Xa6cdyVRYWPf03eiLkq7iHhrPC9Wv+jX+xoMgnvWzuOZz62hz+Hi+j9t4cXddX6NwVsuu52Tb7zJtEsuwThtmu71x1xzLUhJx6uv6l63oviSSvwURYETB7Qu3gWXQN4dgY5mWE3dWrLxyp5jfP2STP5y61Jsb73JtPXnY1p1K+TcBBt/AQ2lut3z1Jp+WweOvXHwDfpcfdyQcYNu9xkLs9HM1fOvpvBoIS29LX6///K0OF77r/PIS4vlvmf38L8v7w26Nf+6PGv36dzN6xE2ayYRq1bR/tLLITPmUVFAJX6Korhc8MqXwBgGV/8+aLt4369o4oGtvTR02PjH3Sv4r4sy6NnyEc6WFmI8Mzav+CVEJMDL+nX5Rl1wAYaYmIE1/aSUvFj9IgvjF7IwYaEu9xiP6zKuwyEdvHbgtYDcP3GahSc+s4p7183jX0WHueWRrbTbgif5a3/5ZUwpKUSuXu2ze8Rcey39R47Qu3u3z+6hKHpTiZ+iTHWlz8CRrXD5TyE6NdDRnEVKyd83H+Kex3aSFGHg9f86j/VZSQB0vvcehuhootat0wpHxMNVv4WmMtjxqC73N1gsRF96KV2Fhcj+fva37qeyrTJgrX0e82LmsTxpOS9WvxiwFieT0cD3rlzIn25fTmVjJz8ssrH/2MmAxDKYq7ub7i1bib7qSoTR6LP7RF96CSIsjM533/PZPRRFbyrxU5SprK8bNvwQUpfD0tsCHc1ZHE4X97+6jx++vp+LFybzvZXWgVm7Ukp6thYRuWoVwmw+dVH2VTDvAij8GfS06hJH5Lnn4urpobdsLy9WvYjVaOWKeVfoUvdEXJ9xPbUna9ndHNgWpytzZvDc59cgJdz05y18UNEc0Hh6du4Eh4Ooc8/16X0MkZGEL19Od5HvFtRWFL2pxE9RprItf4TOY9qEDh12NdBTp62fex7fyeNbD3Pvunn8+Y58LKZT3dD9R4/Sf+wYEWvO6MoTAi79Mdg64MNf6xJLxKqVIATtH33Im4fe5NL0S4kOi9al7om4ZM4lRJmj/D7JYyiLU2P4wRor6dMj+cxjO3h8a23AYuneWoQICyN8+XKf3yty9WrsFRU4WvX5I0NRfC24ftMriuI/Jxvgo9/BomtgzppAR3Oa+vZebvrzVjZVt/DT63L43pULMRhOH3vYvVVrZYlcPUTsKUtg+Z2w/RFt4soEmeLisC5cSMPGd+jq7+L6jOsnXKceIswRXDn3Sv5T+x9O9gW+izXOauDZz63hwuxkfvDKPh54dR9Ol/+7obu3biV8+XIMVqvP7xXp/sOjZ9s2n99LUfQQkomfEOJyIUSlEKJGCPGdIc6vE0LsFkI4hBA3nnHOKYQocT/UPHxl6nr/x+BywMUPBDqS01Q2dnL9nz6ivq2Xf35qBbetGnr3kO6irZiSkwmbmz50RRf8jzZh5d0f6BJXxJrVWMprybCmsTzJ9y1J3ro+83psThtvHRx+azl/irSY+Mud+XzmvLn8c0stX3hiF7Z+/yw0DeA4cQJ7ZaVPJ3UMZl28GENU1MAfIooS7EIu8RNCGIGHgCuARcAnhBCLzih2BLgbeGqIKnqllMvcj6t9GqyiBKuGPVDyJKz6HMTPC3Q0A3YdbuPmv2xFSnj+C+ewNiNxyHLS5dLG961ePfzWbNOSYe3XoOJ1qN084dg6c9MxOiV39OdPeDs4PS2KX0R2fDYvVL8Q6FAGGA2C//3YIu7/+CL+s7+JT/1jB502/+z04Wl5izxzCICPCJOJiJUr6d66dfTCihIEQi7xA1YCNVLKg1LKPuBp4LSFmqSUtVLKUiB41hZQlGAhJbzzfQiPg7XfCHQ0AzZWHeeOv20jLsLMC184h6yU4RfdtVdW4mxvP3t835nWfBmiZ8E739OWrZmA1yMP4DDA8qOmCdWjNyEE12dcT3lrOftP7A90OKf51Llz+d0ty9hR28ptf93GiS67z+/ZvbUIQ1QU1sWLfX4vj8jVq+k/epS+OrWLhxL8gus3mHdmAkcHfV0HrBrD9VYhxE7AAfxcSvnyUIWEEPcC9wIkJydTWFg4rmCH0tXVpWt9ir4m+/uT0LKNnNpNVGXcy7FtJYEOB4BtDQ4eKbUzM8rA13LhQOl2hhqZ53lvIt59j2lAqQTXKO9V0sybWFT+IOXPPkBTyoXjis8hHTxT9wrZaVGkfrCZwjUj39PfYpwxmIWZPxb+kZsTbg5YHEP97MQCX14WxkMlHVz14Aa+WWAlIdx3bQ4JH3yAY948Nm6eeCuvt4wmI9OB3Y/9E5uPZxKP12T/vRbq/Pn+hGLiN1FzpJT1Qoh5wPtCiDIp5VmfMVLKR4BHAAoKCuT69et1C6CwsBA961P0NanfH0cf/Ok+mJ5J5q0/JdNoHv0aH3ui6DB/Lt3Lijnx/O3uAqKtw8fkeW+OPPUU/XPnsu66a0e/gWsd/G0jC+ufY+H134awyDHH+O7hd+k60sX086/E/PhrrF22DGNs7Jjr8aVNmzax8ehGfnPebwg3hQckhuF+dtYD565s5dP/3MGvSyT/+kw+C5L030atr66eAy0tzLz3XuL9+DMspaT64YeZ3drGzCD93TGpf69NAv58f0Kxq7cemD3o61nuY16RUta7/z0IFAJ5eganKEGt5EloPQCX/AiCIOn726aD/M/Le7kwK4nHP7NyxKTPQ/b10bNzl/djuAwGbbmazmPaLN9xePXAqyRFJJF9yY0gJd3bto+rHl+6PuN6Ovs7KTxaGOhQhrQiPZ5n7l1Dv1Nyy1+KqGzs1P0ePUXaODt/je/zEEIQuWo13du2qe3blKAXionfDiBDCDFXCBEG3Ap4NTtXCBEnhLC4n08HzgWCa1CMoviKsx82/xZmFkDmZYGOhkc+PMCP3yjnypwU/nxnPlazdzss9JaVIXt6iBjLrM05a2D+Rdq6hX3dY4qzp7+HLfVbuGTOJUQuXYYhIoLurVvGVIc/LE9aTrw1ng1HNgQ6lGEtSo3m2c+txmQUfOKvRZQ36LsETfeWrRgTpxO2YIGu9Xojcs1qnC0t2Kur/X5vRRmLkEv8pJQO4MvAO0A58KyUcp8Q4odCiKsBhBArhBB1wE3AX4QQ+9yXLwR2CiH2AB+gjfFTiZ8yNZQ+A+1H4PxvBXw/3ocLD/DTNyu4KncG/3drHmaj97+KurcWgRBErlw5tpue/y3oaYGd/xjTZZvrN9Pn6uOitIsQZjMRK1bQE4RLdxgNRi5Mu5BNdZuwO30/iWK85iVG8fS9awgzGrjtr0XsO9ahS71SSrq3bSNy9ZqAzLr2LB/To3bxUIJcyCV+AFLKN6WUmVLK+VLKn7iP/UBK+ar7+Q4p5SwpZaSUMkFKudh9fIuUMkdKudT9rz6beSpKsHM6tF0sZiyFjEsDGspDH9Twi7cr+PjSVP7vlmVjSvpAW7/PumjR2MfYpa2Guetgy++hv9fry9478h5xlriBtfsi1qym7/Bh+hsaxnZ/P7go7SJ6HD1sawjuxYTnTo/kmc+tJtxs5Pa/bWNv/cSTP3tVNc4TJ/y2ft+ZzDNnYk5LU+v5KUEvJBM/RVHGaO/z0HYI1gW2te8PG6r51TuVXLMslQdvXoppjEkfdju9e0rHP4Zr3begqwl2P+5V8X5nP5vqNrF+9nqMBq0rOnKNtlNIMH7Ar0pZRZQ5ivcOvxfoUEY1JyGSp+9dQ2SYidv+WkRZ3cSSv0CN7xsscvVqenbsQDocAYtBUUajEj9FmexcTq21L3kJZF0ZsDAeLjzAb96t4vq8mfz25mVjT/qAsJoa6O8nYqht2ryRfh6knQObfweO0btDtzVuo6u/i4vnXDxwzJKRgTE+nu6i4Fuw12w0s27WOgqPFuJwBX/ykZYQwdP3riY63Mztfyti/7Hxj/nr3lqEeU4a5tRUHSMcm8g1q3F1dWHbuzdgMSjKaFTipyiT3b6X4EQ1rPumNsM1AB7bUssv3q7g6qWp/OqmpRgN42t1DKuo0MbZ5Y9zyzQhtLF+nceg+IlRi793+D0iTBGsmnFqqVBhMBC5ehU9W4uCcgbnRWkX0WZvo7i5ONCheGV2fAT//uxqIi0m7nx0GzXNXWOuQzoc9OzYMfS+zX7kmXDUrcb5KUFMJX6KMpm5XFprX2I2LAzMDoXP7TzK/a/u45JFyfzm5vEnfaAlfuF5eRjCJ7BO3bz1MGsFbH5QW9dwGE6Xkw+OfsDaWWuxGC2nnYtYvRrH8eP0HTw4/jh85LyZ5xFmCAvq2b1nmh0fwZP3rEIIwR1/28bR1p4xXd9bVoaruzug3bwAprg4LAsX0r0l+FqDFcVDJX6KMplVvAbHywPW2vd66TG+/UIpazOm88fbxjZ790yOtjbMR+sm/uEuBJz/beg4Cnv+PWyxkuMltNpauTjt4rPOBfM4vwhzBOfMPIcNRzYEZYvkcOYlRvHEPSvp7Xdy29+KaOyweX2tZyZtxKqxbOLkG5GrV9NbXIyr1/sJRIriTyrxU5TJSkrY+CtIWACLr/P77TeUN/HfT5dQMCeeR+4swGLybp2+4fRs02aqjmn9vuEsuBhS82DTb7T1DYew4cgGzAYz580876xzYbNnY545k+6twdmyc1HaRTR2Nwbd3r2jyU6J5vFPr6Stu5/b/1ZEi5d7+3Zv2Ypl4UJMcXE+jnB0kWtWI/v76dm9O9ChKMqQVOKnKJNV1dvQVAZrvwGGiSVdY7XlQAtfeHI3i1KjefTuAsLDJn7/7qIiXFYr4Tk5Ew9QCG2Gb/thKHvurNNSSjYc3sCa1DVEhUUNWUXkOWvo2b4d6XROPB6drZ+1HqMwhlR3r8fS2bH8/e4V1Lf3cuej2+noHTox93D19tJbUjLQChtoEfn5YDKp9fyUoKUSP0WZrLb8AWLSIOcmv952b30H9z6+i/SECB7/9EqmebENmzd6i0vonzcPYdJpi/GsKyBpkbabxxldohWtFRzrPsZFaRcNe3lEQQGuzk7sB87a6jvgYq2xFCQXhGTiB7Bybjx/ubOAmuZO7n18J7b+4ZNr2/79yP5+IgoK/Bjh8AyRkVgXL6K3uCTQoSjKkFTipyiT0bESOPwRrPocGHVKlLxw5EQPd/9jB9FWE499eiWxEWG61Ovq6cFeXU3/3HRd6gO0Vr/VX4TmfXBo42mn3jvyHgZhYP3s9cNebs3NBcBWVqZfTDq6MO1CDnYc5GBH8E1A8cb5mYn8+qalbDvUyteeKcHpGnq8Ym+p9vqH5+rQEqyT8Nyl9O7bp9bzU4KSSvwUZTIq+hOERcHyO/12y5YuO5/8+zYcLhePf2YlM2ImMPP2DLb9+8Hloj89Xbc6Aa01NGI6bP3TaYffP/L+wN63wwmbMwdDdPRA4hFsLky7ENC+l1B1zbKZ/M9VC3lrbyMPvLpvyMkqtrJSzKmpmKZPD0CEQwvPzUH29mI/EJpJtzK5qcRPUSabkw2w9wXIuwOsMX65Zbfdwaf/uYPGkzYevWsFC5Km6Vq/J7nqnzNH13oxW2HFZ6D6HWipAaC2o5aa9poRu3lBW88vfMkSestK9Y1JJymRKeRMzwmJXTxGcs/aeXxu3Tz+VXSYP7xfc9b53tKygdbXYOEZh2oL0v8bytSmEj9FmWx2/E3brWPV5/xyuz6Hi88/sYt9x07y0G3LyZ+j/8zK3rJSzDNnIqOjda+bFfeAMQy2PQwwMC7O02I2EmtuDvbKKlw275ce8acL0y5k34l9NHY3BjqUCfn25dnaji/vVvHv7UcGjjtaW+mvq9Nnwo+OzEHeGqxMbSrxU5TJpL8Xdv4dsq+C+Hk+v52Uku+8UMqm6hZ+dn0OFy1M9sl9bKVlWH01hisqSevyLXkKelp5/+j7LIxfSGrU6Ft/hefmgtOJbX+5b2KbIM8ahKHc3QtgMAh+cWMu67MS+f5LZbxf0QScGl8ZTOP7AIQQhOfk0Buk4z+VqU0lfooymex5GnpbtUkLfvB/G6p5sbie+y7J5OaC2T65h+PECfrr6wnP8WF33uovQH8PHTv+yt6WvZw/+3yvLrMuWQJAb+ke38U2Aekx6aRNS+OjYx8FOpQJMxsN/On25SxOjeHLTxWzt76D3j2lYDBgXbQo0OGdxZqbg72qSi3krAQdlfgpymQhJRQ9DCm5MOccn9/upeI6fvdeNTcsn8V/XbjAZ/fpLdXGSfm0VSclB+auo6jsMVzSxbmp53p1mTkpCdOMGdiCuEvvnNRz2NG4gz7n8NvThYqIMBN/u6uA2HAzn3lsB+3FJVgWLMAQGRno0M4SnuNpDQ6tRbSVyU8lfooyWdRsgJZKWPMlbakSHyo6eIJvPV/KmnkJ/Oz6HIQP72crK/NPq87qL7GFXqYZLSyZvsTry4K9S++c1HPodfRS0lwS6FB0kRxt5e+fWkG3zUHH7j2YFnv/XvlTeI6nNTh4/28oU5NK/BRlsih6CKJSYPH1Pr3NgeNdfO5fu0iLj+DPd+QTZvLtr5He0jIsGRkYIiJ8eh+54BK2REaxql9iEt7vNBKem0P/0aM42tp8GN34rZyxEpMwTYruXo/slGj+fFEKkfZuXuiJxuF0BTqks5gSEzGlzlAze5WgoxI/RZkMmsvhwPuw8h4w6bNo8lBOdNn51D92YDII/nH3SmIi9NmVYzhSSnrLyvwyeP9Q52EaDbCmtRHqdnh9nTUnuBdyjjRHsjRpKVuPBee+wuOV01kPwGv2OB54beg1/gItPCdXtfgpQUclfooyGWz7M5iskP9pn93C7nDyuX/toumkjb/eVUBagm9b4AD6jxzB1dGB1Q/LdWw5tgWAc10mbQFsL1kXLwYhgvoD/tzUcylvLaeltyXQoejGVlaKsFq59Ko1PFF0hL9/VBvokM4SnptLf10djtbWQIeiKANU4qcooc52EkqfgyU3QGSCT24hpeR/XtrLzsNt/ObmpSxP03+tvqGc2o7L9wv0fnTsI9Kj05mZeweUvwZdzV5dZ4yKxLJgftAu5AzaOD+AooaiAEein97SMqyLFvGtq5Zw2eJkfvLGfjZWHQ90WKfxtFQHa2uwMjWpxE9RQl3Zs9DfDQWf8dktHt18iOd21fGVizL4WO7o69vppbesFBEejmWB72YNA9iddnY27mRN6hrIvxtcDij+l9fXW3NysZWWBWV3I8DChIXEWmLZUr8l0KHoQvb3Y9u/n/CcHAwGwW9vXkZm8jS+/NRuDhzvCnR4A6yLFoHBENStwcrUoxI/RQllUsLOf2hLuMxc7pNbfFDZzE/fLOeKJSn890UZPrnHcGylZVgXL0KYTD69T3FzMTanTVvGZXoGzF0Hu/4JLu8mDYTn5uJsa6O/vt6ncY6XQRhYM2MNW45twSWDbyLEWNmrq5F2O+FLtZbgSIu2zEuY0cA9j+2ko6c/wBFqDJGRWBYsCOrWYGXqCcnETwhxuRCiUghRI4T4zhDn1wkhdgshHEKIG884d5cQotr9uMt/USuKD9TthKa9UPBpnyzhUtPcyVeeKiYrJZrf3LwUg8G3y8QMdqpVx/fdvFvqt2AymFiRskI7kP8paD+iTZjxwkCXXmnwfsCfM/McTthOUN1WHehQJszTgjZ4j95ZcRH8+c586tp6+NJTu4Nmpq81NyeoW4OVqSfkEj8hhBF4CLgCWAR8Qghx5gJfR4C7gafOuDYeuB9YBawE7hdC+GewkqL4ws6/Q9g0yLlx9LJj1N7Txz2P7cRiNvC3uwqICPNtq9uZbFVVyL4+v8zo3XJsC3lJeUSY3RNWsj8GkYna6+sFS0YGwmIJ6i69NTPWAEyKZV16y0oxxsVhnjnztOMr0uP5ybU5bK5p4cdvBMc2euE5uTjb2+mvqwt0KIoChGDih5aw1UgpD0op+4CngWsGF5BS1kopS4Ez/+S7DHhXStkqpWwD3gUu90fQiqK7nlbY9yLk3gyWabpW7XC6+PJTxRxrt/GXO/OZGRuua/3e8AyIt/q4xa+lt4XKtsqBCRCAtiRO3p1Q9RZ0jN59K8xmrIsWBfVCzsmRySyIXTAwezmUefZuHmrh8JtXzOYz583ln1tqeXr7kQBEdzrPHy69QdwarEwtoZj4zQSODvq6zn3M19cqSnDZ8zQ4bFDwKd2r/uU7lVqryXVLyJ8Tr3v93ugtLcMYH495pm8nk3jWtzst8QPIv0sbQ7n7ca/qCc/NwbZvH9Lh0DtE3Zybei67m3bT098T6FDGzdnVjb2mZsQhAN+9Ipu1GdP5wSv72H0ksAtrWxYsQFitQb2tnzK1+LfvJoQIIe4F7gVITk6msLBQt7q7urp0rU/RV0i8P1Kycvsf6Y/OorjiBFQU6lZ1UYODR/bYuSjNRFLXAQoLD+hW91gkFG3FOTOVjRs3DhzzxXvzYsuLRBmiaCxtpFmcvoRLTnweUVsfoUiuRBpG3s3DajIRY7Px0dNP45g1S9cY9RLZG0m/q59/vPcPFocv1r1+f/zsmKuqiJeSaiT7RrjXLbMlFXWSTz+6lQfOsRJrCVw7R9zMmTRu3kx5AH+vhMTvtSnMn+9PKCZ+9cDsQV/Pch/z9tr1Z1xbOFRBKeUjwCMABQUFcv369UMVG5fCwkL0rE/RV0i8P7WbYWM9XPYw65et163a8oaT/HPDR6xIj+Phe1b7fDu24Ti7uqlqaCTlhhtYNui90Pu9cUkX9z97P+vmrOPCdReeXSClG56+jfNTbZB91Yh19c2bx4FH/87isDDigvT/z2rHah59+lG6ErpYv3K97vX742fnxIEDNAMrb78dU9zIQ7Tn55zk+j9t4clDVp4M4P/npq1FtD3zDOefey7C7NvdboYTEr/XpjB/vj+h2NW7A8gQQswVQoQBtwKvenntO8ClQog496SOS93HFCW07Pw7WGNg8XW6Vdne08fn/rWLmHAzD92+PGAfkgC2fftASp8v3FzZWkmrrZVzZ547dIGMy2BaqleTPMyzZ2OMiQnqxXqtJiv5yfkhPcGjt7QM8+zZoyZ9AAtnRPOLG3PZUdvGT97Y74fohmbNzUHabNhragIWg6J4hFziJ6V0AF9GS9jKgWellPuEED8UQlwNIIRYIYSoA24C/iKE2Oe+thX4EVryuAP4ofuYooSOruOw/1VYdjuY9Zl04XRJvvJ0CQ0dvTx8Rz5J06y61Dteno3trUuW+PQ+nokOnhmvZzGatLF+NRug9dCIdQkhsOYG/96s56Sew6GOQzR0NQQ6lHHpLSslfAxb+F29NJXPrp3LY1sP89zOo6Nf4AOeP2CC/f+GMjWEXOIHIKV8U0qZKaWcL6X8ifvYD6SUr7qf75BSzpJSRkopE6SUiwdd+3cp5QL34x+B+h4UZdxKngBXv7bWnE5+859KPqw6zg+vWeK37dhG0ltahjktzatWnYnYcmwLmXGZJEYkDl8o705tjcTdj41aX3hODvbqalw9wTt5wjOJJRRn9zqOH8dxrAHrGJf4+fbl2ZwzP4Hvv7yX0rp23wQ3AvOsWRhjY9VCzkpQCMnET1GmLClh12Mw5zxIzNSlynf2NfKnwgN8YuVsPrEyTZc6J6p3bxnhPm7t63X0UtxcfPZs3jPFzITMK6D4CXCOvCOENTcHXC5s+wPXrTiaBbELSApPYmvD1kCHMma9e/cCjKnFD8BkNPDH25aTGGXhC0/spq27zxfhDUtrDc5RM3uVoKASP0UJJYe3QNshWH6nLtXVtnTzjWf3sHRWDA9crf8sz/FwtrdrrTqLz1yXXV97ju+h39XPypSVoxdefid0H4fqd0csZl2oxWyrqNQjRJ8QQrBixgp2NO4Iud0k7BUVIASWrOwxXxsfGcbDdyzneKedrz5TgtPl3+/dunAR9oMHcfX5N+lUlDOpxE9RQknxE9pOHQuvnnBVvX1OPv/ELoxGwUO3L8diGnm5En/xJE3j+XAfix2NOzAKI8uTvdjjeMElEJWsvf4jMCUlYoyLw1YRHLtGDGdF8gpaba0c6hh53GKwsZVXYE6bjTEqclzX586K5YGrF/Nh1XH+8L5/t66zZmeBw0GfmuChBJhK/BQlVNg7Yf/LsOR6CIuYUFVSSv7n5b1UNnXyu1uWMStuYvXpyV5ZAbg/KH1oZ+NOFicsJtLsRRJhNEHuLVD9DnQ1D1tMCIF1YTb2IG7xAwZaOXc07ghwJGNjq6zAmr1wQnV8YuVsrl8+k//bUM3GquM6RTY6S7b2h0wwtwYrU4NK/BQlVOx7Cfp7tMkGE/T0jqO8sLuO/7owg/VZSToEpx9beQXGxOmYpk/32T16Hb2UtpRSkFLg/UV5d4DLAaXPjFjMkpWNvaoqqHfwmDVtFskRyexoCp3Ez9nVTf/hIxP+g0AIwU+uzSEreRpffbqY+vZenSIcWVhaGiI8POhbg5XJTyV+ihIqip+E6ZkwawzJyhDK6jq4/9V9rM2YzlcvytApOP3YKisn3KozmpLmEhwuBytSVnh/UWIWzFqhdfeOMDbOujAb2ddHX23txAP1ESEEK1JCa5yfvaoKONVyNhHhYUYeviMfp1PyxSd3Y3c4J1znaITRiDUzM+hbg5XJTyV+ihIKWqrhaJG2dt8QG9N7q6Onny88uYvpkWH83615GA3jr8sXZF8f9poan3fzesb35SXlje3CZbfD8Qqo3z1sEc/YRFt5xURC9LkVKdo4v4MdBwMdilc8LWVWHRI/gLnTI/nVTUvZc7SdH7/un1Y4y8JsbBUVIZNsK5OTSvwUJRSUPAnCCEtvHXcVUkq+/twemk7aeOj25cRHhukYoD7shw5Bf78urToj2dk0hvF9gy25Hkzh2lqKw7DMm4swmwfGKgarFclaa2eojPOzV1RijInBlJKiW52XL0nh3nXz+FfRYV7bc0y3eodjzc7GdfIkjobQXDxbmRxU4qcowc7pgJJ/Q8YlMG38H3qPbj7Ee+VNfPeKheQFwSLNQ7GVu1t1Fvquq7env4eylrKxje/zsMbAoquh7AXoH3psmDCbCctYEPQtfgPj/EIk8bNVVGDJzkZMoMV7KN+8LIv8OXF898UyDrV061r3mawDEzyC+/+GMrmpxE9Rgt2B96GrUZtcME7FR9r4+VsVXLY4mU+dm65fbDqzl1cgrFbC5szx2T32HN8z9vF9g+XdAfYOKH992CLW7IVB/+HuGee3s2ln0Hc9SqcTe1WVbt28g5mNBv7wiTxMRsEXn9yNrd934/0smZkgxMAfOIoSCCrxU5RgV/IERCRAxmXjury9p48vP1VMSoyVX964VPcWEz3ZKiuxZGYijL5bU3Dc4/s85pwHsWkjdvdas7NxnjiB47j/lgsZj1AZ59d3+DDSZsOy0DdDAFJjw/ntzUspbzjJj1733a4rhogIwubMURM8lIBSiZ+iBLPuE1DxpraGnGnsY/KklHzjuVKaO208dNtyYsLNPghSH1JK7OXlWLN8vH5f004WTx/H+D4Pg0Gb5HFwI7QfGbKIxT05Jdhb/TytnsHe3TswBMCHYz8vzE7mc+fP48ltR3jVh+P9LNnZQf//QpncVOKnKMGs7Dlw9Y+7m9czru97Vy5k6exYfWPTmaOpCWdHh89adeDU+D7PxIZxW3ab9m/Jv4c8HSpjuWZFzSIlMiXoEz97RSWYzVjmzfPpfb5xqXu83wulHDze5ZN7WLOz6T96FGeXb+pXlNGoxE9RglnxEzBjGSSPfR9dz7i+yxencPc56bqHpjd/tOqUHB/H+n1DiU2Dueu07l6X66zTxuhozKmp2IN8gocQghXJwT/Oz1ZRgWX+fESYb2eie8b7hZkMfOmpYp+M9/O0BtsrVXevEhgq8VOUYNW4F5rKtG7FMero7ee//q2N6/vFjblBPa7Pw+5uHbNk+q6rd2fjzomN7xss7w6tq/fI1iFPWxYG/wQPODXO70D7gUCHMixbRblP/yAYTBvvt4zyhpP89E39J2F4ZqwH+6xvZfJSiZ+iBKvSp8FggiU3jOkyKSXfe7GMxg4bf/hEXlCP6xvMVlGJeU4axqhxjr3zwo7GHSyevpgIsw57E2dfBeZI2DNMd29WFn21tbh6/bMl2Hh5lrUJ1u3bHC0tOI+3DLSU+cMF2Ul8du1cHt96mLf3NupatykpCWNcXNCv86hMXirxU5Rg5HJC2fOQcSlEJozp0n9vP8obZQ1847KsoF2vbyi2inKsWb4d37e3Ze/Ex/d5hEVqa/rtf2XINf0sC7PB5cJeXa3P/Xwk2Mf52dwzYH29jd+ZvnlZNrmzYvjW83uoa+vRrV4hBJbsLNXipwSMSvwUJRgd2gidDdps3jGobOzk/72m7cN771rfDoTXk7Orm/4jR7H6cGJHyfESHFKH8X2D5d4C9pNQ+dZZpwa69IK8u3dgnF9jcI7z87SM+XobvzOFmQz8/tY8XBK++nQJDufZYznHy5q9EHt1NdLh0K1ORfGWSvwUJRjteQYsMZB5udeX9PY5+fJTu5lmNfPbm5dhCLJ9eEdir6oCKQf2ufWFnY07MQmTPuP7POaug2kzoPSZs06ZZ87EEBU1MHYxmK1IWUGbvS0ox/nZyiswzZiBMTbW7/dOnx7JT65bwq7DbfzuPf1abq3ZWUi7nb7aWt3qVBRvqcRPUYJNXzeUvwaLrwWz1evLfvj6PmqOd/G7W5aROM3iu/h8wFbh2arNd4mfruP7PAxGyLkJat6D7pbTToVSl97Aen5BOM7PnxM7hnLNspncXDCLhwpr2FLTMvoFXrBkqwkeSuCoxE9Rgk3569DfDUtv9fqS1/Yc49/bj/L58+dzXsZ0HwbnG/aKSgwxMZhSxr8X8UgGxvfp2c3rsfRWcDlg7wtnnbJmZWOvrEQOseRLMJkZNZMZkTOCbpyfy2aj71CtXyd2DOWBqxczb3okX32mhJYu+4Trs8ybizCb1QQPJSBU4qcowab0aW2duNmrvSp+tLWH771URl5aLPddkunj4HzDVlGBNTvbZ8vO7Dm+B4d0UJBcoH/lyYshOQf2PH3WKevCbFw9PfQfPar/fXUkhKAguYBdTbuCapyfvboGnE6/T+w4U0SYiT98Yjkdvf186/nSCb9GwmwmLGOBavFTAkIlfooSTDob4WChNmnAMPqPp8Pp4r+fKUFK+P2teZiNofcjLZ1O7FVVPu3O2928G4MwsDRxqW9usPQWOLYbWk4fB+YZs2gLgb1Zlycvp9XWSu3J2kCHMiBQEzuGsig1mu9dkc37Fc08tqV2wvVZsxdiU4s4KwEQep8SgBDiciFEpRCiRgjxnSHOW4QQz7jPbxNCpLuPpwsheoUQJe7Hn/0evKKMpOw5kC7I9a6b9w/v17DrcBs/uW4Js+N1HLvmR32HDyNtNiw+TPyKm4rJissiKizKNzfIuQmE4axWP0vGAjAaB8YwBrPlScsBKG4uDnAkp9jKKzBERGCePTvQoQBw1znpXJidxE/fqqC84eSE6rJmZ+FsacFx/LhO0SmKd0Iu8RNCGIGHgCuARcAnhBCLzij2GaBNSrkAeBD4xaBzB6SUy9yPz/slaEXx1p5nYGY+TF8watEdta384f1qrs+byTXLZvohON8Y2KrNRxM7+l39lLaUsjx5uU/qB2BaCsxbD6XPnraFm8FqxTJvbtBv3QYwN2YucZY4djXtCnQoA2yVFViysxFetH77gxCCX92YS0y4ma/8e2JbullCZD9nZfIJjp+msVkJ1EgpD0op+4CngWvOKHMN8Jj7+fPARSIU9qxSpramfdoWbV609nX09PPfT5cwKy6CH167xA/B+Y69ohLMZizzfLPuYPmJcnodvQMtWj6Teyt0nL2FmyUrOyS69IQQ5CXlBU2Ln3S5sJdXBEU372AJURZ+c9NSqpu7+PEb+8ddjzVL+75U4qf4mynQAYzDTGDwSOk6YNVwZaSUDiFEB+DZ/mCuEKIYOAn8j5Ry01A3EULcC9wLkJycTGFhoW7fQFdXl671KfoK1Psz78A/mSWMbO1IoX+E+0speXiPncYOJ99fZWXn1s3+C9IHYrd8hCE5mY1btoxadjzvzYaTGwCwHbBRWDu2a8fC4IzhXIOVpncepCqrf+B4RJiZaQ0NbHzjDWSk77aj00N0VzRHO4/yynuvEGOKGfP1ev7sGI8fZ3p3N4eFgYog/H15ebqZJ4qOEG9vYnny+D5Kp8fHc3Tjh+zNyNA5urOpz53g5s/3JxQTv4loANKklCeEEPnAy0KIxVLKswZrSCkfAR4BKCgokOvXr9ctiMLCQvSsT9FXQN4flxN2fR4yLuXcS68eseizO4+yvbGUb16WxacvGL1LONhV/+8PiDznHHK9eM3H8968+P6LpDnSuPqikV9XXXRcS2rlW6Seu3pgDcYuk5mjL75EQcJ0Ilef+TdqcEk4nsDLb76MZYGF9enrx3y9nj87J999l3og95qrCc/N1aVOPZ1znovrH/6Ixyt7ue3y1aTEeL/mpsfRZcvoO3yYPD/8vlGfO8HNn+9PKHb11gODR/rOch8bsowQwgTEACeklHYp5QkAKeUu4AAQmutfKJNL7WZti7alI2/Rdqilmwde3cfqefF8/vz5fgrOdxytrTiOH/fZxA6XdFHcXKzvbh0jWXoL2Dug6u2BQ56uylCY4JGdkE24KZzipsB399orKkEILH5oDRsPz5Zu9n4XX3+uBJdr7Eu8WLOz6Dt0CJfN5oMIFWVooZj47QAyhBBzhRBhwK3Aq2eUeRW4y/38RuB9KaUUQiS6J4cghJgHZAAH/RS3ogyv7FkImzbiFm397qVbzEYDv715GcYQ2pJtOPaqKgAsmb75cK/tqKXd3k5+cr5P6j/L3PMhMkmbne1mmj4dY3w89mr9tvzyFbPBTO70XHY37w50KNirqghLS8MQHh7oUIY1LzGK+z++iI9qTvDo5kNjvt6SmQkuF/YDwbdVnjJ5hVziJ6V0AF8G3gHKgWellPuEED8UQnj6ch4FEoQQNcB9gGfJl3VAqRCiBG3Sx+ellK1+/QYU5Uz9Ntj/Giz8OJiH/5D7/YZq9hxt56fX5ZAaG7wfhmPhSfysmb5peN/VrM1Q9emM3sEMRlhyA1T/B3rbBw5bMjOxVwV/4gfaa1XVVkVnX2dA47BXVWmJUZC7ZcVsLl2UzK/eqWTfsY4xXev5/kLl/4YyOYRc4gcgpXxTSpkppZwvpfyJ+9gPpJSvup/bpJQ3SSkXSClXSikPuo+/IKVc7F7KZbmU8rVAfh+KAkDNu1r3YO5NwxbZUdvKQx/UcGP+LK7KneHH4HzLVlWFMS4O43TfbDNX3FRMvDWetGlpPql/SLk3gbNP22/ZzZKZgb2mJui3bgPIS8rDJV3sOb4nYDG4envpO3IkJBI/IQQ/vyGX2AgzX326ZExLvISlpSHCwgb+AFIUfwjJxE9RJpXSZ7XuwfR1Q54+aTu1dMsDVy/2c3C+Za+qxpKZ6bOt2nY37yY/Od9n9Q8pdTnEzzutu9eamYns7Q36rdsAliYuxSiM7G4KXHev/cBBkDIkEj+A+MgwfnPzUmqau/jZm96P5RQmE2EL5qvET/ErlfgpSiDZOqDqHVhyPRiHnmR//yv7aDxp43e3LiPKMnkm4kuXC3tNjc8+3Bu7G6nvqvffxA4PIbSdPA59CCcbgEFdeiEwzi/CHEF2fHZAx/n5euynL6zNSOQz583lsa2H+aCi2evrrBmZKvFT/EolfooSSOWvgdMOOTcPefqVknpeKq7nKxdmsDwtzs/B+VZ/fT2yp8dnH+6ehYj9Nr5vsJybAAn7XgTAMl+bgW0LkQ/45cnLKTteRp+zLyD3t1dVISwWwtL82EWvg29elkV2yjS++fweWrrsXl1jyczEcfw4jrY2H0enKBqV+ClKIJU9B3FzYebZyUl9ey//8/Je8ufE8aULQn/pljMNTOzw0XIdu5p2EWGKICsuADs/TM+AGcsGunsNkZGYZ88OmUH8+Un59Ln62H9i/DtTTIS9qgrL/PkIozEg9x8vq9nI7z+Rx0mbg289X4qUoy/x4vnDJxRag5XJQSV+ihIonY1ad2DOTVr34CBOl+S+Z0qQEh68eRkm4+T7UfUkfmELfJP47W7ezdLEpZgMAeoez7kJjhVDSw3gmdkbGi1+y5KWAQRs315bdWjM6B1KZvI0vntFNu9XNPPU9iOjllczexV/m3yfJooSKva+CNLl7hY83d82HWTboVbu//gi0hIiAhCc79mqqjDPmoUxSv9tzDrsHdS01QSmm9djyfWAGGj1s2Rm0Hf4MC67d12AgZQQnkB6dHpA9u11tLXhPN4SsokfwF1r0lmbMZ0fv17OweNdI5Y1JSVhiIkJmT8KlNCnEj9FCZSy52DGUkg8/QNu37EOfv2fSq5YksKN+bMCFJzv2aurffbhvuf4HiTSfws3DyU6Feau1d5nKbW1Cp1O+g6Gxprx+cn5FDcX45L+XYLG0/IVyomfwSD49U1LsZgNfO2ZEvqdw7+GQgisGRkq8VP8RiV+ihIIJw7Asd1ntfbZ+p3899MlxEWE8dPrcvy7DIkfufr66DtU67OJHbuadmEymFgyfYlP6vdazk3QegCOFQ/q0guND/i8pDxO9p2kpr3Gr/cNxRm9Q0mOtvLT63LYU9fBHzaM3I1ryczEXl3t1ZhARZkolfgpSiCUPQcIbZeHQX7xdgXVzV386qalxEWGBSY2P+g7eBCcTp/t2FHcXMyihEWEmwK8w8nCq8EYBmXPa4v1ms0hNbMX8Pu+vfaqKowxMZgSE/16X1+4MmcGN+bP4o8f1LDr8PCbRFkyM3F1d+M4dsyP0SlTlUr8FMXfpNQSv/TztO5Atw+rjvOPj2q5+5x0zs8M/Q+9kQy06vhgRq/daWdvy16WJwVwfJ9HeCxkXAp7X0AYDYTNnx8yg/hnRc0iMTxxYNs7f/Fs1TZZWrvv//giUmPD+doze+iyO4Ys42ndDJU/CpTQphI/RfG3hhI4UXNaN29bdx/feG4PC5Ki+M4V2YGLzU/sVVVgNhOWnq573Xtb9tLv6g+OxA+097lLm8FtyQydsVxCCJYnL/frDh7S5fLp2M9AmGY18+Aty6hr6+GHr+0bsoznD6BQ+aNACW0q8VMUfyt7HgxmWHQ1AFJKvv9yGW09ffzulmVYzaG1dtl42KqqsMybhzCbda/bMxPV7zt2DCfzMgibBnufx5qZiaOpCWdHR6Cj8kpeUh5NPU00dDX45X79xxpw9fRMqsQPYEV6PF9YP59nd9bxzr7Gs84bp03DlDojZP4oUEKbSvwUxZ9cLm0Zl4xLIFzbiePlknreLGvka5dksmRmTIAD9A97te+2atvdtJu5MXOJtcb6pP4xM4fDwo/B/tewzJ8LhM5ivZ7k2V/bt02WiR1D+epFmSyZGc13XyyjudN21nm1dZviLyrxUxR/OrIFOo8NTOqob+/lBy/vo2BOHJ9bN/l25xiK8+RJHA0NPvlwd0kXJcdLgqeb12PJjWDvwGLSWntCZSxXZlwmEaYIv63n58uxn4EWZjLw4M3L6LY7+M4LZWfN4LVkZmI/dAjZF5ht8pSpQyV+iuJPZc+DOQKyrsDlknz92RJcUvLgLcswGibHYPbReFq7fPHhfrD9IJ19nQM7TwSNeedDRAKmhncxTJsWMi07JoOJ3MRcSppL/HI/e1UV5tRUjFFRfrmfv2UkT+M77l09/r396GnnLJkZ4HBgr60NTHDKlKESP0XxF2c/7H8Fsq6EsEge3XyIooOt3H/1YmbHT87dOYYysEevD7p6i48H2fg+D6MZFl2LqHoby4LQmdkL2mtZ3V5NZ1+nz+9lD+Gt2rzl2dXjR6/v51BL98BxtXWb4i8q8VMUfznwAfS2Qs6NVDSe5FfvVHLpomRumsS7cwzFVlWFYdo0TDNm6F53cVMx8dZ40qal6V73hOXcCI5eLInmkFqsNy8pD5d0UXq81Kf3kX192A/VTvrEz2AQ/OrGpYSZtF09HO5dPSxz54LJFDKtwUroUomfovjL3ufBGos9fT3//XQJ0eFmfnb95N2dYzj26mosGRk++b6Lm4vJS8oLztd09mqInonVdAxXZyeOxrNndwaj3MRcDMLg83F+9kO14HBM+sQPICXGyk+uW0LJ0Xb+VHgAABEWhmVuukr8FJ9TiZ+i+ENfD1S8AYuu5rfv11LR2Mkvb8whIcoS6Mj8SkqJvaraJxM7WnpbqOuqC75uXg+DAZZcj6VvLxA6W7dFmiPJisvy+Ti/yTyjdygfy03l2mWp/N+GavYcbQfAomb2Kn7gVeInhNgmhHhYCPFZIUS+EGLy7iWlKL5Q/Q70dVGecCmPfHiQ21alcWF2cqCj8jtHUxOukyd90qrjaZEKuokdgy25EUu0tpRHqMzsBe01LW0ppd/V77N72KuqwGTC4oNFvYPV/7tmCUnTLHzt2RJ6+5xYMjPpP3YMZ1dXoENTJjFvW/xWAPcCfwa2A51CiBIhxN+FEP8lhAjSP7EVJUiUPY8rKpl7P7SSFh/B969cGOiIAmJgYocPZvTubtqNxWhhUfwi3evWzYylGGfMxzTNGFKD+JcnLafX0UtVq++SVXtVFZa5cxFhU6ddISbczK9vWsrB4938/K3ygdbOUFnnUQlN3iZ+9wMvA4cBAZiBXOBu4HfATiFEnRDim0KIyb/tgKKMha0Dqt/lo7C11J/s47c3LyPSYgp0VAFxqjtP/xa/kuYSlkxfgtmo/24guhFCa/WL6sZePvT2XcHI04rqy3F+tikwo3co5y6YzqfPnctjWw+zx5QAqJm9im95lfhJKX8kpbxBSjkPiAcuAr4BPAmUu4ulAj8HPhRC+HQRJiHE5UKISiFEjRDiO0OctwghnnGf3yaESB907rvu45VCiMt8GaeiANrYPqed3zTk8sX1C8ifExfoiALGXl2NKTkZY4y+O5T09PdQ0VoRvOP7Bsu5EWtMP321tch+33Wd6iklMoUZkTN8lvg5u7pwHGuYkokfwLcuz2JBUhTf2NSEiIhQ4/wUnxrz5A4pZbuU8gMp5W+llHdKKZcAicBXgVZgNfAXneMc4G5RfAi4AlgEfEIIcWbfzmeANinlAuBB4BfuaxcBtwKLgcuBP6kWSsXX7CXPcpRkHDPy+MpFU2Pg+nBsVdU++XDf27IXh3SERuI3PQPLnBlIh4u+w4cDHY3X8pLyKG4u9skyNJ4WrqkyseNMVrOR392yjJZuBw1xqSrxU3zK28kd0SOdl1K2Sin/AKwEWoBbhRC5OsQ3lJVAjZTyoJSyD3gauOaMMtcAj7mfPw9cJLT1Ha4BnpZS2qWUh4Aad32K4hOy6zim2o287lrDg7fkEWaauhPppcNB34EDPvlw97RELU1cqnvdvmBZfTkA9t2bAxyJ9/KS8jjee5z6rnrd6z61VdvUbPEDWDIzhq9dkkmxKYHO8oqQWedRCT3eDjRqE0IcAEqAYs9DStk0uJCU8pAQ4gG0FrlPAV/TL9QBM4HBe93UAauGKyOldAghOoAE9/GiM66d6YMYx+Tf964jtex4oMNQ3ASwUcf6uuOSSPrybWQkT9Ox1tDTd/gwsq/PJ1u1FR8vZkHsAmIs+nYh+0rYRXfDD5/DVvQ20TffHehwvOJpTS1uLmbWNH0XHbdXVWGIjMQ8M1XXekPN59bN45EnBcbDnWxaswip43KUev9eU8ZPfP1LrLvxywG7v7eJnwAWuB83DBwUogktCdzNqYTwFbTE73J8k/j5hRDiXrSZzCQnJ1NYWKhb3V1dXafVZ4+IoiWxTbf6leAR1eli7kHY13VY1/9DociycxexwN6TnTjG+Vqc+bMD4JIudh/bzfLI5SH1Gs+MMdKztzRkYnZJF1Zh5c2SN5l2dOg/YoZ6f7wRt2MHJCezcaNKTboX11DYIUiyq1FIk1V/S/tZPyfj/dkZD28Tv1wgb9BjGRADpKCNtbt8UNmTgAQyhBB3AruACimlS6eY64HZg76e5T42VJk6IYTJHesJL68FQEr5CPAIQEFBgVy/fr0esQNQWFjIafXpWLcycWe9PxNQU/gq/Z//NtGuNt3qDFXHS0tpMRo555abMVjGt3D1UO9NZWslvUd6uWrZVayfv37I64JRXeYCbPv2s35hIiQvDnQ4Xsl/N5/GnsZh/y+P52dHSkn1t7/DtMsuY+kU/xnB0ccTVS2cvDqaZ+/cpmvVev5eU/Tnz/fH21m9e6WU/5JS3ielvEBKGQfMB24EfgK8BTSitQzGuP81AP8EytDW/dsqhHhIh5h3oCWVc90LSd8KvHpGmVeBu9zPbwTel9qAiVfRxh9ahBBzgQy0dQkVxSfSl60DoGXfrgBHEni2qirC5swZd9I3HM+OEkG9cPMQrPnr6O824dr570CH4rVlScs40H6Ak30ndavT0XwcZ0fHlJ3RO1h/zXuUmg0sSwyBSUpKyBr3SHMp5SEp5YtSyv+VUn5MSjkTrQXwasCF1up3xF08HG0c3ucnGrCU0gF8GXgHbSmZZ6WU+4QQPxRCXO0u9iiQIISoAe4DvuO+dh/wLLAfeBv4kpTSOdGYFGU4pthYumItuGpqAx1KwNl9NKO3+Hgx08OnMytK33FnvmZZsgwA++aXIUQG8ucl5SGR7Gneo1udU22rtpFUlf6LXoOB5Qs+FuhQlEnMq65ed5ftbqB8pC5bKWWzEKIdLaFskVLOdc8IzgOWu/+dMCnlm8CbZxz7waDnNuCmYa79CVorpaL4hSN9BvHHajnRe4KE8IRAhxMQrp4e+o8eJebaMyfgT1xxUzF5SXloE/dDh2eSi63uBOF1O2H2igBHNLqc6TkYhZHi5mLWzlqrS52nZvRO8cSvr4fdDdsgNpJlMwoCHY0yiXnb4vcYUMqpLts/CSHuEUIsF0IMLJMvhEgAfo/W2vcugJTypJRyo5TyQSnlJ/X+BhQl2EVlL2FmC5Q0Tt3uXntNDUip+4d7U3cTx7qPhcb6fWcwz5qFCA/HftICe58PdDheiTBHkB2fretCzvaqKkyJiZjipu7C5gBUvU2xWZBqiSclMiXQ0SiTmLeJXwPauD1Pl+3n0BZp3gF0CSGqhBB70JZHWea+5g/6hqoooSk1dzVhTqgq/TDQoQTMwB69Onf1Fh/XEpBQTPyEwYAlIwO7PRH2vQSu0Bh1kpeUx96WvfS79Nl1ZKpu1XYmWfY8xdZwls04c3UyRdGXt5M7ZgLJwJXA94EXgVpO7du7AMgBLGjj+74vpSwasjJFmWIisxcCcGLf7gBHEjj26mpEeDjm2bNHLzwGJc0lhJvCyYrP0rVef7FkZmBvldDVBLWbAh2OV/KS8rA5bZSfKB+98Cik00lfzQGV+PW2U3fofVqMguXJ+YGORpnkvN4pXkp5HG1CxNueY0KIGLQWvmy0bdvagP9IKdUO04riFjZ/PlIIxMEj2Bw2rCZroEPyO1tVFZYFCxAGfXcu2d20m9zpuZgN5tELByFrZiYdz7+AwzkNU9nzMG99oEMa1eCFnHMTJ7ZBU9/hI9qi3lM98at4neIw7Wcj1GanK6FnQr+FpZQd7vF7f5FS/lhK+ZBK+hTldAaLBeesJGY1u9jbsjfQ4QSENqNX3/F93f3dVLZVkpccet28Hp4xj/Zpa6D8VXDYAxzR6BIjEpk9bbYu4/zUxA63sucpjp7ONPM0FsQuCHQ0yiQ3dTcOVRQ/ispaRFqz1HVQfKhwnDiB88QJ3cf37Tm+B5d0heT4Pg9PS5dNLABbB9RsCHBE3slLyqO4uXjC+8naq6rAYMCyYL5OkYWgrmY4tJHiyGnkJuViNKgdOxTfUomfovhB1MLFJHVA2dGdgQ7F706t06bzxI7mYgzCwNLEpbrW60+mhASMCQnYW4Hw+JCZ3bs8aTmttlYOnzw8oXrs1VWEpaVhsE694Q8D9r1Mh4ADjpMsT1oe6GiUKUAlforiB5bMTAwSWsqLcem2e2FosFdroz/07s4rbiomKy6LSHOkrvX6myUzA3vNAVh0DVS+BX3dgQ5pVJ7u9Ym2YPtqUe+Qsvd5SpK1n41Qbr1WQodK/BTFDzzdnNOPdXOg/UCAo/EvW1UVxvh4TNOn61Znv6uf0pbSSfFBac3MxF5Tg1x8PfT3aMlfkJsbPZdYSyy7m8c/U93V20vfkSNTO/FrPwJHt1GcsgCTMLFk+pJAR6RMASrxUxQ/MM+eDRYLacen3jg/X7TqVLZW0uvoDemJHR6WzExkby/9YjZMS4Wy5wId0qiEEAPj/MbLXnPAJ4t6h5S9LwBQbHCyMGEh4abwAAekTAUq8VMUPxAGA9aMDOafME+pxE+6XNhranSf0bu7SWtpypsEm9kPbN1WUwNLroea96CnNcBRjW550nIOnzxMS2/LuK5Xe/QCZc/TN6uAve01k6L1WgkNKvFTFD+xZGaSdnzi46JCSX9dHbKnR/8dO5qLmRk1k+TIZF3rDQTLggUghJYI5dwELgfsfyXQYY3K09pa0lwyruvtVVUIq5WwtDQdowohTfuhaS/756+jz9WnEj/Fb1Tipyh+YsnMIKKzj86mOpq6mwIdjl/4YkavlJLdzbsnzQxIQ0QE5tmzsVdVw4ylMD0TyoJ/du+i+EVYjJZxj/OzV1dhmT8fYZyiy5eUPQfCSHF0PKAWblb8RyV+iuInnlav2cflwB6zk93AjN75+q3TdqTzCK221kkxvs/DkpmhvVZCaK1+hzdDR12gwxqR2WhmyfQlFDeN7/+yrXoKz+iVUkvu51/A7o4q0qalMT1cv8lPijISlfgpip94PuTmnzCPu3ss1NiqqjDPno0hUr8lVzzj+yZLix9ofxT0HT6My26HJTdoB90D/4PZ8qTllLeW09PfM6brHG1tOI+3TN2JHUe3Q8cR5JIbKWkuUd28il+pxE9R/MSYkIAxLo7ckzEDyctk54sZvSXHS4ixxDA3Zq6u9QaSJSMDnE76DhyAhPkwsyAkZvfmJeXhlM4xb0Vor/TNot4ho+xZMIVzaOYS2u3tLE+ePH/EKMFPJX6K4idCCCyZmcw5DpVtlXT3B/9CvRPh6uujr7bWJzN68xLzMIjJ8+trYOs295hIcm6CxjJorghgVKNbmrQUgRjzOL8pPaPX2Q/7XoKsKyhp14ZCqPF9ij9Nnt+cihICLJmZRNd3IF1O9hzfE+hwfKrvwAFwOnWd0dvp7KT2ZO2kGt8HEDZnDiIsTJvgAbD4OhCGoG/1iw6LJiMuY8wz1e3VVRhjYzElJvoosiB2sBB6TkDOTexu2k2sJZa50ZOn9VoJfirxUxQ/smRmIGx2kjsMk7671xczeg/aDwKTa3wfgDCZCJs/f+A1Y1oyzFuvJX5SBjS20eQl5VHSXILD5fD6GltVFZbMTIQQPowsSJU+C9ZYWHAxu5p2sTxp+dR8HZSAUYmfoviRp/VrTW8qu5p2BTga37JXVyPMZl3XaTtoP0iYIYxFCYt0qzNYWD0zez1yboL2w1C3M3BBeWF50nJ6HD1Ut1WPXhhtUe++6pqpOb6vrxsq3oBF19Bkb6Ouq4785PxAR6VMMSrxUxQ/Cpu/AIC8znjKWsroc/YFOCLfsVVVETZ/PsJs1q3Og7aDLJm+hDBjmG51BgtLRgaOpiacHR3ageyPgcka9N29nokJ3o7z6z92DFdPz9Sc0Vv5FvR3Q+7NA69XfopK/BT/UomfoviRMSoS86xZzDkusDvt7DuxL9Ah+Yw2o1e/D/deRy9H+o5M2qUvPC1gA9291mjIvBz2vQhO77tR/S0lMoUZkTO8Huc3pSd2lD2v7cecdg67mnYRYYogKy4r0FEpU4xK/BTFzyyZmUTVaXuxTtbuXmdHB47GRl0nduxt2YsL16Rd+uKsmb2gdfd2H4dDhYEJykt5SXkUNxUjvRiPOJD4ZUyxrt6eVqh5F3JuAIOBXU27yEvKw2QwBToyZYoJqcRPCBEvhHhXCFHt/jdumHJ3uctUCyHuGnS8UAhRKYQocT+S/Be9omgsmRk4Dx8lIzJ90iZ+Azt26Jj4eSbDLE1cqludwcSUnIwhOvrUzF6AjEvAGhP0W7gtT1pOc28zJxwnRi1rr6rCPHMmxij9FvUOCftf1vZhzrmJDnsHNe01k/aPGCW4hVTiB3wH2CClzAA2uL8+jRAiHrgfWAWsBO4/I0G8XUq5zP1o9kfQijKYNTMTnE7Ody6gpLkEp8sZ6JB0Z/PBjN7dzbuZYZ5BjCVGtzqDibbOY8aprl4AkwUWXQPlr0Hf2HbH8CfP8joH7AdGLWurqpqa4/tKn9X2YU7JHfgjRk3sUAIh1BK/a4DH3M8fA64dosxlwLtSylYpZRvwLnC5f8JTlNFZsrMBWNYRQ1d/F1VtVaNcEXrsFZUYY2IwJSfrUp/D5aC4uZgFlgW61BesrFnZ2CsrkS7XqYO5t0BflzYbNEgtiF1AjCWGGlvNiOVcNht9h2qxLMz2U2RBoq0WjmzV3kuhLXhtNmh7HSuKv4Xa4IJkKWWD+3kjMNSnykzg6KCv69zHPP4hhHACLwA/lsMMShFC3AvcC5CcnExhYeEEQz+lq6tL1/oUffn8/XG5SDKbMRXXwWp45qNnWB+93nf3C4D47dtxpaSwceNGXeqrtdfS6+hlFrMm9c+O1SCI6elh8/PP40xyj0SRLlZbEuku/BNlrcG74HGaIY3q3uoR3x/T4cMkOJ0ccDrZP4nfxzPNqX2GucDWnjTshYV80PABaeY0tm7a6rcY1OdOcPPn+xN0iZ8Q4j0gZYhT3x/8hZRSCiHGurLp7VLKeiHENLTE707g8aEKSikfAR4BKCgokOvXrx/jrYZXWFiInvUp+vLH+3No4UIiexykRqZyMvrkpPr/IJ1OKv/7a8TdcgvLdPq+/rH3H9AIS2KXTKrX6ky9CdOp/dcTLI2JJXrw9+n6JNbND7I+f6G2uHMQOrzvML/e+WsWrVxEUsTQw6fbn3+eBiD/hhsImzPHvwEGipRQeh/MOY81l99MT38P9f+u51NLPsX65ev9Fob63Alu/nx/gq6rV0p5sZRyyRCPV4AmIcQMAPe/Q43RqwdmD/p6lvsYUkrPv53AU2hjABXF76zZ2dgqKshPWs6upl1ezYYMFX2HDyNtNl2783Y27SQ9Op1oY7RudQYjS8YCMBqxVZSffiL3VpAu2Bu8kzwKUgoA2Nk4/ILTtvIKDBERmGfPHrbMpFO/C1oPwNJbANhzfA8O6VDj+5SACbrEbxSvAp5ZuncBrwxR5h3gUiFEnHtSx6XAO0IIkxBiOoAQwgx8DNjrh5gV5SyW7CxcHR2sMs6n1dZK7cnaQIekG1u5lrRYs/VJ/JwuJ8VNxVPig9JgtWKZNxd7ecXpJxIzITUP9jwdmMC8kB2XjVVYR5ypbquowJKVhTCE2kfPBOx5WluIe9E1gDZJySAMk3Z2uhL8Qu2n7+fAJUKIauBi99cIIQqEEH8DkFK2Aj8CdrgfP3Qfs6AlgKVACVor4F/9/h0oCmDNXgjA4rYoYHKt52evqASzGcu8ebrUV9VWRWd/50CL0mRnycrGVll59oncW6GxFJrLzz4XBIwGI/Ms89jZNHSLn3S5sFdUYJ1KEzscfbD3Bci6QluWB+1nPTs+m6iwqAAHp0xVIZX4SSlPSCkvklJmuLuEW93Hd0op7xlU7u9SygXuxz/cx7qllPlSylwp5WIp5VellJNvHQ0lJFgyM0EIYg63Em+NH1jeYTKwVVRgWbAAEabPtmqeRKIgeWokftaF2TgaGnC2t59+YskNIIxB3eq3wLqAgx0HOdF79np+/fX1uLq7B2a1Twk170Fvq5a0A/3OfkqPl7I8Sa3fpwROSCV+ijJZGKMiMafNxl5ZSX5y/iRr8avAmqXfNlS7mnYxM2omKZFDzfmafCxZWmJkqzij1S8qERZcrO3dO3i5lyCywKottzPU/2dbhdZ9rdcQgJBQ+jRETIcFFwGw78Q+7E77lPkjRglOKvFTlACxZi/UJngk53Os+xgNXQ2jXxTkHC0tOI4f1607zyVd7GraNaU+KK3ZWtJ81gQP0CYInKyH2k1+jso7aWFphJvCh0z87OUVYDBMncWbe9uh8m2tpdZoBk4lxJ4FrxUlEFTipygBYs3Oov/IEZZHaknSrubQb/XztFJ5Wq0m6kD7Adrt7VNiYoeHafp0TImJZ0/wAMi6EizRUPqM/wPzglEYWZq4dMhxfraKCsLS0zGEhwcgsgDY/zI47QOzeUFL/ObFzCPeGh+4uJQpTyV+ihIgnrFOs5tdRJmjJsU4P3ulpztPn65eTwvJVJnY4WHJHmaChzkcFl0N+18J2i3cCpILqG6rpsPecdpxe0XF1Orm3fMMJGRAqjaez+lyUtJcMqX+iFGCk0r8FCVArAu1mb19lVUsS1o2Kcb52corMM2YgTE2Vpf6djbtJDkimVlRs3SpL1RYs7OxHziA7Os7+2TurUG9hVt+cj4Sedr/Z2dHB/3Hjk2drdraDsORLVprnxAAVLdX09nfyfJkNbFDCSyV+ClKgJiSkzHGxGCv0CZ4HOw4SKutNdBhTYi9Ur9WHSklOxt3kp+cj3B/eE4Vluws6O/HfvDg2SfnnAsxs7WJA0EoJzGHMEPYaYmfp/VyyrT4lT6r/Ztz88ChgdbrKTReVQlOKvFTlAARQmBZeGqCBxDS3b0umw37wUNa0qKDwycPc8J2Ysp188Kp1mDbUOP8DAbIuQkOvA+djX6ObHQWo4WcxJzTxvnZp9KMXim1pHzOuRB3alu6XU27SI1MnTKz05XgpRI/RQkga1YW9qoqFsdkE24KZ0fjjkCHNG726hpwOgcWp56oqbZ+32Bhc+YgrFbsQ83sBVh2m7aFW5BO8ihILqCitYLOvk5AS2CNCQmYEhMDHJkfHN0OJ2q098jNJV3sbNzJipQVAQxMUTQq8VOUALIszEba7ci6YyxPWs72xu2BDmnc9J7YsbNpJwnWBNKj03WpL5QIoxFLZubZa/l5TM+A2aug+EmthSnIFKQU4JIuSppLALDpOAQg6BX/C8yRsOjagUPVbdW02dtYNWNV4OJSFDeV+ClKAA3u0luRsoKa9hpaelsCHNX42MorMEREYJ49e8J1TeXxfR7W7GxsFRXI4RK7ZbdDSyXUDb1FWiDlTs/FJEzsbNqJ7Oujr7pmamzV1tcN+16CxdeC5dSWbJ4/6FSLnxIMVOKnKAFkmTsXzGbslRUDrQGh2t1rq6jAkp2NMEz810p9Vz1NPU1TcnyfhyU7C1dHB46GYRb2XnwdmCOg5An/BuaFCHMEi6cvZmfTTuyHDiH7+3Vb2zGo7X9Vm3Gdd8dph7c3bGdO9Bw1vk8JCirxU5QAEmFhWBYswFZeQXZ8NtPM09jWsC3QYY2ZdLnc67Tp180LU3N8n4dnrOSw3b3WaFh0DZS9EJRr+hUkF7C/ZT8n95UCTI0Wv+InIH4epK0ZOORwOdjZtJOVKSsDGJiinKISP0UJMGtWFraKCkwGE/kp+SE5zq+/vh5Xd/fAotQTtatpFzGWGObHztelvlBkycwEhtm6zWPZ7dDXCeWv+Skq7+Un5+OQDuqLNyPCwghLTw90SL7VehAOb9YmdQwanlB+opyu/i5WzlCJnxIcVOKnKAFmXZiN073H7aqUVRztPMqxrmOBDmtMbDov17GjcQf5SfkYxNT9FWWMisQ8Jw37cC1+4F4yJD0ou3vzkvIwCiPd+/dhycxEmEyBDsm3Sp4CBCy97bTD2xq1FvwVyWp8nxIcpu5vVUUJEpZBXXqeVoFQa/Wzl1eAwYAlI2PCdR3tPEp9V72aAYnW3etJqodkMGitfoc+hLZav8XljaiwKBYnLCK8tnHyd/O6nFDyb5h/IcTMPO3U9obtZMRlkBCeEKDgFOV0KvFTlADzjIuzV1awIHYB8dZ4tjeEVuJnq6ggLD0dQ3j4hOvyjHFcnbp6wnWFOmt2Fv1HjuDs6hq+0NJPAEJLPILM+ZYcIrqdsCA90KH41sFCOFl31qSOPmcfxc3FrEpRf8QowUMlfooSYMaYGEypM7TlUISBFSkr2Na4bfhlPIKQNrFDn1adooYikiKSmBs9V5f6QplnzKS9qmr4QrGzYd56ravR5fJPYF4q6JwOQE2iM8CR+FjJk2CNhawrTztcerwUm9OmJnYoQUUlfooSBKxZ2QNdeitTVtLc08zhk4cDHJV3nB0d9B87psvEDpd0sa1hG6tnrJ6y6/cN5kmmbeUjTPAAraWp4wjUfuiHqLw3s6EPgK0RoTVmdUx626D8dW0bPbP1tFPbG7djEAbyU/IDFJyinE0lfooSBKwLs+k7dAiXzTYwti1UxvnZKrXJB3qM46psraTd3s7qGaqbF8CUkoIxJmbkCR4A2VeBNUbbySOI9FfV0D7dyub20N2DelRlz4PTflY3L2g/w4viFxEdFh2AwBRlaCrxU5QgYMnOBpcLe3U1adPSSIpICpnEz67jjN6ihiIAlfi5CSGwLBxlggeAORyW3Ajlr4Ktwz/BecFeXo5rfhqHOg7R2N0Y6HB8o+RJSF4CM5aedrjX0cue43tYMUPN5lWCi0r8FCUIDGzdtm8/QghWpaxiR+MOXDK4xmwNxbZvP8bp0zElJk64rqKGIhbELiAxYuJ1TRbW7GzslZXI/v6RC+bdAQ4blD7rn8BG4ezqou/IERJytW7OUFyYfFQNpXCsWHvtzxiaUNxcjMPlUBM7lKCjEj9FCQLmWbMwxsbSu7cMgJUzVtJqa6WmvSbAkY2ut6yM8JycCddjd9rZ3bRbtfadITw3B9nXh726euSCqXlaq9POf0AQTAyy7d0HUjJz5XrirfEDrbmTyq5/gMkKubecdWp7w3ZMwkReUl4AAlOU4YVU4ieEiBdCvCuEqHb/GzdMubeFEO1CiNfPOD5XCLFNCFEjhHhGCBHmn8gVZWRCCKy5Odj2aNtbeWYBBvuyLs7OTvoOHiQ8d+KJX0lzCTanTSV+Z7Dm5gLQW1o6ckEhoODT0LwP6gK/37Mn3oicXFalrKKooSikZqqPyt6pta4uvh4i4s86vb1xO7mJuUSYIwIQnKIML6QSP+A7wAYpZQawwf31UH4F3DnE8V8AD0opFwBtwGd8EqWijEN4Ti72mhqcXd2kRqUye9rsgVX/g5Vt714ArDm5E66rqKEIozBSkDJ19+cdinnmTIxxcfSWlo1eeMmNEDYNdv7d94GNwlZWinlOGsbYWFanrqalt4UD7QcCHZZ+yp6Dvi4o+NRZpzr7Otl3Yp/apk0JSqGW+F0DPOZ+/hhw7VCFpJQbgM7Bx4S2NsSFwPOjXa8ogRCemwNSYtu/D9Ba/XY27sThcgQ4suF5kpHwJYsnXFfRsSJyE3OJNEdOuK7JZKA1uGyUFj8ASxTk3gx7X4SeVt8HN4Le0jLC3X8QrJmxBmDydPdKqSXXyUtg1tmTN3Y17cIlXWr9PiUohVrilyylbHA/bwSSx3BtAtAupfR8itYBM0coryh+ZXWPk7OVacnUqhmr6OrvoqJ1lBmdAdRbVkrYnDkYY2MnVE+HvYN9J/apbt5haK3BB3B2dY9euOBT2vIiewK3k0d/UxOOpqaBIQAzomYwJ3rO5En86ndDY5n2Wg+x3uS2hm1YjBaWJi4d4mJFCayg2zVbCPEekDLEqe8P/kJKKYUQPhswIoS4F7gXIDk5mcLCQt3q7urq0rU+RV+BfH8SpidwZMP7lM2fT79Tm8X51EdPcWnMpQGJZzTTd+ykLytrwq9XSU8JEoml0TJiXVP1ZycMSZyUFD35JP1ZmaOWz4vOwvzhQ2y3LRoyMfEVz/tjKSkhFih3OOh3v1+zXbMpqi9iwwcbMAqj32LyhayK35NksLKlIxXnEP8f3z/2PunmdLZs2uL/4IYxVX92QoU/35+gS/yklBcPd04I0SSEmCGlbBBCzACax1D1CSBWCGFyt/rNAupHiOMR4BGAgoICuX79+jHcamSFhYXoWZ+ir0C+P/UrV9JTUkKe+/6PvfoYDWENQfn/pb+piZqODtIvvpj4Cca3uWgzEW0RfPLST2I2mIctN1V/dhxLl1L9x4fINhpI8Ob7j/0avPx51qebYO5an8fn4Xl/mncXc8JkYs1tt2GwartZOA472Fy4mfjF8aE907W3DTZvgbxbWXvxlWedbupu4tjzx/ha7tdYv2S9/+MbxlT92QkV/nx/Qq2r91XgLvfzu4BXvL1QatPJPgBuHM/1iuIP1pxcHMcacBw/DsDamWspaS6hs69zlCv9zzNrU48ZvUUNRRSkFIyY9E1lprg4zLNnezfBA2DxtdresQGa5NFbVoo1M3Mg6QNYkbICgaDoWIh39+55Bhy9kH/2pA6Aj459BGg/u4oSjEIt8fs5cIkQohq42P01QogCIcTfPIWEEJuA54CLhBB1QojL3Ke+DdwnhKhBG/P3qF+jV5RReJKo3jJttux5M8/DIR1BOTbKVloGJhMW9+LT43Ws6xiHTx5W4/tGEZ6TQ2+Zl4mfORyW3Qblr0HXWDpGJk66XNjK9mI94w+CGEsMixMWB+X/Za9Jqa3dl7ocUpcNWWRT3SaSI5JZELvAv7EpipdCKvGTUp6QUl4kpcyQUl4spWx1H98ppbxnULm1UspEKWW4lHKWlPId9/GDUsqVUsoFUsqbpJT2QH0vijIU68KFYDTS657BuTRpKdPM09hcvznAkZ2tt6wMa1YWBotlQvV4dnRQid/IrLk5OBoa6G/2MpHL/xS4+qH4Cd8Gdoa+2lpcXV0DM3oHW526mtLjpXT3ezFJJRgd2QrHK7T1EofQ7+pna8NW1s5ai/Dj2EpFGYuQSvwUZbIzRERgycjQWtMAs8HM6tTVbK7bHFSL32qtOmWEL534+n1bG7aSYE1QLSSjCHcv5OxZO3FUiZmQvlZroXL5b+u/kYYArJ6xGod0sKtpl9/i0dXOv4MlBpZcP+TpkuYSuvu7OW/meX4OTFG8pxI/RQkyni49T6K3duZamnubqWqrCnBkp/QdOoSru3vCCze7pIttDdtYnbpatZCMYqA1eLQdPAYr+BS0H4EDG3wX2BlspWUYIiIImzfvrHPLkpZhMVrYemyr3+LRTXcL7H8Flt4KYUOvNbmpfhMmg0m1XitBTSV+ihJkwpfm4jp5kv7DhwEGWg821W8KZFinGVi4eYITO8pPlNNqa+Wc1HP0CGtSM4SHY8nKHGgN9kr2xyEyCbb9xXeBnaG3rAxrTg7CePaSLRajhYLkgqAcujCqXf8AZ9+w3bygje/LT85Xi5ArQU0lfooSZDytaJ6B/IkRiSyMX8imumBK/PZgiIwkbO7cCdVTWFeIQRjUDEgvhefk0rt3L9LbrltTGKy4B2reheOVvg0OoL8fW0XFiH8QnD/7fGpP1lLbUev7ePTi6IPtf4P5F0JS9pBFGrsbqWmvUf+XlaCnEj9FCTKWBfMRERGnLd1x3szz2HN8Dx32jgBGdoqt1N2qY5jYr5CNRzeyNHEpcdY4nSKb3MJzc3CdPEmfuzXYKwWfBqMFih72XWBuprp66O8f2IVmKOfPOh+AjXUbfR6Pbva9CF2NsPpLwxbxtMirxE8JdirxU5QgI4xGwhctwjZoLNfaWWtxSidbGwI/Nsplt2OrrCR8hA93bzR2N1LeWj6QCCijO3NbP69EJWr79+552uf795pra4FTE1GGkhqVSmZcJh8c/cCnsehGStj6EEzPggUXDVtsU90mUiNTmRszsVZwRfE1lfgpShCy5uZiKy9H9vUBkDM9h+iwaDbXBX5slL28HByOs9ZpG6uNR7UWn/Wz1+sQ1dRgmX92a7BXVn9RW3R41z98E5ibubYWU2IipuSRt1E/f9b5lDSXBE0L9ogOb4HGUlj9hWG3v+tz9rGtYZtaxkUJCSrxU5QgFJ6bg+zrw1ZVDYDJYOKc1HPYXL8Zl/Tf0hxDOTWxY2IzegvrCpkVNYt5MWfP/lSGJoxGwhcvHljn0WvJi2DeBbD9r9p4NR8x19Zizc0dNflZP3s9TukMqglLwyr6E4THQe4twxbZ3bybHkePWsZFCQkq8VOUIBQ+0KV3enfvCdsJKlorAhUWoE06MSUlYR6lVWckPf09bG/YzvrZ61ULyRhZc3Ow7z/VGuy11V+EzgbY/7JP4nKePImpqcmrIQBLpi8hwZow0OobtFoPQsUb2jjJsIhhi22u24zZYGZlyko/Bqco46MSP0UJQqbUVIwJCad16Z2bei5AwGf32kpLJ9zNu7VhK32uPtXNOw7hObnI/n5slWNc13HBxZCQoY1X88Fi4J6Fpb1Z4scgDKybtY6P6j+i39Wveyy62fYXMJhgxWdHLLapfhMrUlYQYR4+OVSUYKESP0UJQkII90LOp1r8EsITWJywOKDdY872dvoOHx5yO66x2Hh0I9PM01ievFynyKaOU/s5j7G712DQxqk1lGhbj+nM80eKdckSr8qfP/t8Ovs72d20W/dYdGHr0La7W3I9RM8Ytlh9Vz0HOw6qbl4lZKjET1GClDU3h74DB3F2dQ0cWztrLWUtZbTb2gMSU+/efcDEFm52SRcb6zZy7sxzMRvMeoU2ZZhmzMA4ffrYFnL2WPoJbbxa0Z90j6u3rAxHcjLG6Givyq+ZsYYwQxiFRwt1j0UXu/8FfV1aF/kIPBOu1DIuSqhQiZ+iBKnwnFyQEps72QJtPT+XdLHl2JaAxGQrKwUhvG7VGUpZSxmttlbOn62WcRmPU63B40j8wiIg/1PauLW2Wt1iklLSW7qH/vR0r6+JMEewasYqCo8WBtU+1AA4HVo375xzIXXZiEU31W9iVtQs5kTP8U9sijJBKvFTlCAVnqMlV7179gwcW5KwhFhLLB/WfxiQmHr3lBI2dy7GadPGXcfGoxsxCqNqIZmA8Nwc+g4exHny5NgvXvlZEAZdF3R2NDbiPN4ypsQPtNm9dV11HOw4qFssuqh4DTqOaF3jI7A5bGxv3K6WcVFCikr8FCVIGWNjsWQsoGf79lPHDEbOn3U+G49upM/pu2U5hiIdDnp27iQiP39C9RTWFZKXlEeMJUanyKae8Px8kJKenTvHfnF0KuTcDLseg65mXeLp3rYNgP6MBWO6bt2sdQDB1d3rcsGHv4aEBZB15YhFP6r/iF5HLxfMvsBPwSnKxKnET1GCWMTqNfTs2oVr0NIdl8+9nK7+Lj6q/8ivsdj27cPV1UXkmtXjruNY1zGq26rVbN4JCl+2DGG10r21aHwVrP06OO2w5Q+6xNOztQhjXByO1NQxXZcSmcLC+IXBtX1b1VvQtBfWfRMMxhGLvl37NvHWeFakrPBTcIoycSrxU5QgFrlmNdJmo7ekZODYqhmriA6L5p3D7/g1Fk+SEbF6/Imfp2VHbdM2MYawMCLy8+kpGufs3OkLYMmNsONR6D4xoViklHQXFRGxepU2c3iMzp99PnuO76HN1jahOHQhJWz8JcTN1V6fEfQ6etlYt5GL0y7GZDD5KUBFmTiV+ClKEItYsQIMBnqKTrXsmA1mLp5zMR8c+QCbw+a3WLqLirBkZWGKjx93HRvrNpIenU56TLp+gU1RkWtWY6+uwXH8+PgqWPcN6O+BoocmFEffoVocTU1Erl4zruvXz1qPS7qCYxeP6ne15W7Wfh2MIydzH9Z9SK+jl8vSL/NPbIqiE5X4KUoQM06bhjVnyVldepelX0aPo8dv3b0um43e3buJnEBrX1dfF9sbt6vWPp1EuBOt7qJt46sgMQsWXwvbHoGe1nHH0e1udRzvEICFCQtJDE8M/Dg/KWHjLyAmDZbeOmrxd2rfIcGaQH7yxMa8Koq/qcRPUYJc5Oo19JaV4ezqHji2MmUlcZY43qn1T3dvb3Exsq+PiAmM7yusK8ThcnBBmhoIrwfrwmwMMTEDide4rPsm9HVqS5eMU8/WIsypqZhnzx7X9QZhYP3s9Wyu30xPf8+445iwgx9A/U5Y+zUwjry+ZE9/D5vqNnHJnEswjjIOUFGCjUr8FCXIRa5ZDQ4HPTt3DBwzGUxcPOdiCusK6XX0+jyG7q1FYDIRUTD+QeyvH3ydGZEzyEvK0zGyqUsYjUSuXEn31q3jXwcveTFkf0xb2sXWMebLpdNJ97ZtRKxZPaHlTK6adxW9jl42HNkw7jomxDO2L3omLLt91OIb6zZic9q4fO7lfghOUfSlEj9FCXLheXkIi4Werae37FyWfhm9jl6/7N3bXVREeE4OxqjIcV3f0tvC1mNb+di8j2EQ6teOXiLWrMZxrIH+I0fGX8n53wJ7B2x/ZMyX2vaX4zp5ctzj+zzykvKYGTWT1w++PqF6xq12s7aN3XlfA5Nl1OJvH3qbpPAk9UeMEpLUb2BFCXIGi4Xw5XlnjfMrSC4g3hrP27Vv+/T+zpMnse3dO6FlXN469BYu6eJj8z6mY2SKJ+Ea97IuADOWQublsPUhsHeO6dKB8X2rV43//mjdvVfOvZKihiKO94xzsspEbPwFRKVA3p2jFu3q62Jz/WYuTb9U/RGjhKSQ+l8rhIgXQrwrhKh2/xs3TLm3hRDtQojXzzj+TyHEISFEifuxzC+BK8oERa5eg72qCseJU0tvGA1GLplzCZvqNvl0bFTPjh3gck1oGZfXDrzGooRFzIudp2NkStjcdEzJyXQXTSDxA1j3Leht05Z3GYOerUVYMhZgSkyc2P2Bj83/GC7p4s1Db064rjE5vBVqN8G5XwWzddTiHxz9gD5Xn5rNq4SskEr8gO8AG6SUGcAG99dD+RUw3J9u35RSLnM/SnwQo6LoztPaduYH/OXpl2Nz2ny6AG73lq2I8HDCly0b1/UH2g9Q3lquWvt8QAhB5OrV9BQVIV2u8Vc0Kx8WXAwf/U5LAL3g6uujZ/fugdnFEzUvZh6LExb7t7tXStjwQ4hMgvy7vbrkP7X/ISUyhdzEXN/Gpig+EmqJ3zXAY+7njwHXDlVISrkBGFufhaIEMevixRimTTttPT/QxkYlhif6dHZvd1EREfn5GMLCxnX96wdfxyiMXDH3Cp0jUwAiz1mDs70de2XlxCq6+AHobde2K/NCb3EJ0mYjco0+iR/Ax+d/nIrWCqrbqnWrc0Tlr8GRLXDB9yAsYtTiJ/tOsvnYZi6do7p5ldAVasuNJ0spG9zPG4HkcdTxEyHED3C3GEop7UMVEkLcC9wLkJycTGFh4ThuNbSuri5d61P0FazvT8y8efR98AGVZ8S2yLSIjUc28vb7b2M1jN5VNRaG9nYSDxzgxNJcDo7jNXFJFy/Uv0CWJYu92/ZOOJ5gfW8CySAliUDpE0/Sc8nFE6orc8bFpBT9mR3OJfRGzBixbOSrrxJpMLDbbkO635OJvj/RzmgMGPhT4Z+4Ju6acdfjDeHqZ+X2b+KMnMOuk2kD38NItnVtw+FykNiaGHL/D9XPTnDz6/sjpQyqB/AesHeIxzVA+xll20aoZz3w+hnHZgACsKC1GP7Am5jy8/Olnj744ANd61P0Fazvz4l/PSH3Z2VL+9Gjpx3f3bRbLvnnEvnagdd0v2f7K6/I/VnZsmfv3nFdv71hu1zyzyXy9QOv6xJPsL43gVZzxZXy8Gc/O/GKTjZK+eMZUj59+6hFD91yqzx08y2nHdPj/fnCu1+QFz17kXS6nBOua0Qf/UHK+6OlrH7P60s+/+7n5WXPXyZdLpcPA/MN9bMT3PR+f4CdcpicJujaqqWUF0splwzxeAVoEkLMAHD/2zzGuhvcr4kd+AewUv/vQFF8Y2Cc3xnLuixNXEpyRDLvHNK/u7d7axHGmBisCxeO6/rXD75OhCmCC9Mu1DkyZbDI1avp2bET2dc3sYqmJWsLGJe/BrXD7wrj7Oqit6xsQgt6D+fj8z9OU08TOxt36l73gO4T2rp9Cy6BBRd5dUmHvYOiY0Vcmn7phNYsVJRAC7rEbxSvAne5n98FvDKWiwcljQJtfODE+54UxU/C5s3DlJhIzxlLd3iWwthUv4mm7ibd7iel1Mb3rVqFMIz9V4Xdaec/tf/h4jkXE24K1y0u5WwRa1Yje3vpLS2deGVrvgzRs+Cd78EwE0Z6tu8Ap3PC6/cNZf3s9USaI307yWPjL6CvCy79sdeXvFLzCg7p4Mq5V/ouLkXxg1BL/H4OXCKEqAYudn+NEKJACPE3TyEhxCbgOeAiIUSdEMIz7/5JIUQZUAZMB7z/qVeUABNCELFmNd1DzOC8KesmXNLFs1XP6na//sOHcTQ0jHv9vsKjhXT1d6nZvH4QuXIlGAwTW8/PwxwOF98PDSVQNvT/p+6irQiLhfC8ZRO/3xnCTeFcnHYx7x5+F5vDpnv9HK+CHX/TZvEmZXt1iUu6eLryafKS8siO9+4aRQlWIZX4SSlPSCkvklJmuLuEW93Hd0op7xlUbq2UMlFKGS6lnCWlfMd9/EIpZY676/gOKWVXoL4XRRmPyNVrcLa2Yq8+fdbj7GmzOX/W+Txf9Tx9zgl297l5lo4Z7/p9rx98naTwJFamqBEVvmaMicG6aNHE1/PzWHIjpC6H9/4f9J29RmTP1iIi8pdjsIy+y8V4fGz+x+jq76KwrlD/yt/9AYRFwvrven3J5vrNHO08ym3Zt+kfj6L4WUglfooy1Q03zg/gEws/QautVbelXbq3FmFKSSEsPX3M17bZ2thct5kr512pNrH3k8g1q+ndswdXd/fEKzMY4LKfQucx2PrH0045WlqwV1frtn7fUFYkryApIonXD+jc3XuwEKregrVfhyjvF51+qvwpksKTuGiOd+MBFSWYqcRPUUKIecYMwtLT6f7w7P1518xYw9yYuTxZ/qRnFvu4ufr66C4qInL16nENZH/twGs4pEN18/pR5Jo14HAM+UfBuMxZA4uugc0PQuvBgcNdmza776f/xA4Po8HIVfOu4qP6j2jsbtSnUocd3vo2xKbBqs97fdmhjkN8dOwjbsq6CbPBrE8sihJAKvFTlBATfeUVdBcV0d90+kQOIQSfyP4E+07so7RlYoP8uwoLcXV0EH3V2Aey9zv7eXz/4xQkF5AVnzWhOBTvRaxYgTEhgY5XXtWv0st+BgYzvPLlgYkeHa+8gnn2bKxLluh3nyHcknULEsnj+x/Xp8LCn8HxCrjqQa+2ZvP4d8W/MRvM3Jh5oz5xKEqAqcRPUUJMzDXXgMtFx6tnf8BfPf9qIs2RPFX+1ITu0fHSy5gSE4k855wxX/v6wddp6mninpx7Ri+s6EaYzcR87GN0FhbiaPNu27VRxcyEy38Khz+CHX+lv76enqIiYq69xudLmsyMmsmVc6/k+arnabe1T6yyul3w0f9B3p2Q4f0i1119XbxS8wqXpV/G9PDpE4tBUYKESvwUJcSEzZlDeH4+HS+9fFaXbqQ5kusWXMd/Dv+Hlt6WcdXvaGmh68MPibnmaoRxbOPznC4nf9/7dxbGL+Sc1LEnjcrExFx/HfT3c/KNN/WrdNntkHEpvHs/HU//U7vPNdfqV/8IPr3k0/Q6evl3xb/HX0m/DV7+PExLhct+MqZLXznwCj2OHm5fePv4768oQUYlfooSgmKvu/b/t3fn0VFXdx/H39/sQBYgEISwClKQLWgUBBdEUdQKocWNtqLAY59qrVZ5XKAVdwHXYm0tBSu2WhRtkCLFiogWKhFXNhERkS0StiQEAmS5zx8zUQiTMJNMJhPm8zonJzNz7+/+vnPuycw39/5+93J440YO+li37epuV1NaXsqcL+bUqO2C+fOhrIyUrKyAj128ZTGbCjcxptcYLXJbDxJ+8APiT+1OQXZ28Bo1g8t/h4uKI3/ObBr3O5O4tunBa78aXZp1YVC7Qby47kUOlBx7d7Ff3nkIdq2HYdMgIcXvw8pdObPXzaZ3i970bFG309oioaTET6QBSho6FEtIIH/u3GPKOiR34Oz0s3ll/SuUlJUE1K5zjoJ/ZJPQuzfxXboEfOyMVTNon9SeIe2HBHSsBE/TrBEcXLOGg+vXB6/R5DYUn3wTJfmlpPRKDl67fhjXaxwFhwp4df2rgR+85QPPXcmnX+f3Dh0V3t/+PpsKN3FN92sCP69IGFPiJ9IARScmkjRkCIVvLKD80KFjykd1G8Wu4l289c1bAbV76PPPObR+PSlZwwOOaXnuctbuXsv1Pa/XEi71KPmHl0FMDAVzA9rY6LjyV+3DYqNI3vcK7NoQ1Lar06dlHzJbZTJr7azA/pEpKYa5v/DsQhLADh0VXlr3EqkJqVzc4eLjVxZpQJT4iTRQTUdkUV5YSNHixceUDUwfSIfkDry47sWA2szPnuu5SeDSwO/mnblqJmmN0hjWeVjAx0rwxDRvTuKg8yiYNw9XWhqUNsuLi9m3cCHJF19EVKMET0JVFpy2/TGu1zjyDuQFto3b4gdh9wYY/jTEJwV0vi2FW/jP1v94lnCJ1hIucmJR4ifSQDXu14+Y1q3J93E9V5RF8ZPuP2HlzpUs2bLEr/bc4cMUzp9P4gUXEN20aUCxrNq5ipxvc7i2x7XERccFdKwEX9MRIyjbtYuipUuD0t6+RYso37+flCuvgUsfg60fwL8nBqVtfwxoM4Duzbvz3OrnKCsvO/4Bq1/zTPGeMQ5OHhTw+aZ9Mo3YqFiu6HpF4MGKhDklfiINlEVHkzJsGPuXLqMkL++Y8pFdR9I5pTOTP5hMcWnxcdsreu89yvburdE078zVM0mOS9ZaZ2Ei8ZxziG7WLGjTvQXZ2cSmp9M4MxN6XwFn/RJynoUP/hyU9o/HzBjTawybCjexeMuxI9xH2bICsn8B7Qd4dh8J0PLc5SzctJBxvcaR1jithhGLhC8lfiINWEqWZ02/wn/+85iy2KhYJvafyLaibcxcNfO4beVnzyW6RQsSzz47oBg25m/k7c1vc023a2gS2ySgY6VuWFwcyZf/kKK336YsP79WbZVs387+95eTkpWFRXm/MobcD12HenbC2LCo9gH7YUj7IbRPas+MVTOq3pkmfzPMvgaSW8NVf4OYwPYSLikr4eGch2mb2Jbre14fhKhFwo8SP5EGLL5TJxr17Ut+drbPL8MzTjqDy06+jOdWP8c3hd9U2U7pnj0UvfsuKcOGYTExfp/fOccTHz1BQnQCo7prA/tw0nTECFxJCQULaremX8G8eeDc0SPBUdHw4xmQ1h3mXA9562oZ7fFFR0UzpucY1u5ey4Kvfbyng4Xw0lVQehhGzYEmqQGfY9baWXxd8DV397ubhBj/d/cQaUiU+Ik0cCkjsji84SsOrl7js3x85njio+N5JOeRKkdKCufPh9LSgKd556yfw7tb3+WW026heULzgGOXupPQvTvx3bpRkD23xm045yjInkvjM84grl27owvjk2DUyxDbCF66EvbXbMHwQAzvMpzeLXvz0PKH2F60/fuCslJ4dQzs/AKunAUtuwbcdm5RLtNXTmdwu8Gc2/bcIEYtEl6U+Ik0cMmXXILFx1e5aG+LRi34Zd9fsmz7MhZt9j0tl589l4QePUjo6v8X5sb8jTy64lEGthmo0b4wlZI1nIOrVnFoQ82WXyn+5FMOf/NN1Yt5p7SFq/8ORTtg9iiiyg7XPFg/xETFMPmcyZS5MiYsneC50cM5z40mG96CSx+FzufXqO2pK6binOPOM+8MctQi4UWJn0gDF52URNJFF5Gfnc2hjRt91rnqB1fRrXk3pnww5ZgdEAoXLODQ55+T8uMf+X3OkrIS7vrPXTSKacQDAx8gyvRREo5SLr8ci40l74knq74urgqurIydTz5JVOPGJF1czVp2bU+HEc/Clhx6r5wE+3fXMurqtUtqx8T+E/lox0f8ZdUMmP9rz40m/W+EM8bWqM2l25ayaPMift7n57RJbBPkiEXCiz6tRU4AaePHExUfz7bbbve5oHNMVAwT+01kx4EdPLvy2e9eP7x1K7n3TKJRRgbNrvB/6YqnP32az/d8zr0D7qVl45ZBeQ8SfDGpqaSNv52ixYvZ+9JLAR27e/p0DqxYQavf/IboxOPctNNjBIx8juTCL2HGBXW+wPPlJ1/Oxe0G88ynv2f1qhfh7F/DRYHtw1vhUNkhHs55mI7JHRl96uggRyoSfpT4iZwAYlul0fqRhzm0bh15jz/us05GWgYjuozgr2v+yvLc5biSErbfPh6ANo89hsX6t1BtTm4Oz69+niu6XsHg9oOD9h6kbjS79lqanHcueVOmcvCLL/w65sDHn7Dz98+QfNllpIzI8u9EPX/MpxkPwqFCT/K3aVnNgz4OK9jKb9flkFpaxl2dunPgvDsgKvCvs3JXzqMrHmXLvi1M7D9RizVLRFDiJ3KCSDr/fJr97GfsfeGv7HvnHZ91bjv9NjqmdOTGRTeS89DtFH/2Ga3vv4+4tul+naPgUAETlk6gQ3IHxmeOD2b4UkfMjDYPP0xUSrJnRLi4+jUdywoL2T5+PLGtW3PSvZMwM7/PVZjSDca9DYlp8MJw+Gx2bcM/1raPYcYFpBRs55GMW9l8OJ+pK6YG3ExJWQkTlk7g5S9eZvSpo+nfun/wYxUJQ0r8RE4gaeNvJ75bN3InTKRkx7GLOjdNaMrzQ5/nsr0dSJr9Frsv7Euyn9uz5Rblcus7t7KneA+Tz51M49jGwQ5f6khMairpU6ZweONGdjwyucp6zjly75lESV4e6Y8/RnRSYFudAdC8E4z9N7TvD9k/h7k3Qf6WWkTvdWgfvDsV/nKpZ32+sf/mjNNvYEzPMbz25Wv84dM/UFLu316+B0oOcPPim3lj4xvcctot3J55e+3jE2kglPiJnECi4uNJf/wxyouL2X7nnbiyY7e3anKgnJ+9upuCVo25tc9Knv7k6Wov/C935cxeN5us17NYs3sN9w28jx6pPerybUgdaDJgAKnjxpL/yisULnzTZ52C115j38KFtPzVr2jUp0/NT9aoGfz0HzDgZlj1Cjx9Giy4A/btCLytkmJYNg2e6g3vPARdLvCMKqZ1A+CmjJu4pNMl/PGzP3L1/KtZs9v3skYV9hzcw9g3x/J+7vvcP+B+xvUaF9CopkhD5/9KrSLSIMR37kyriRP49rf3sPN300gacuFR5bue+QPle/PpM/slfljwKtNXTmd38W5uzLiRFo1aHHWH7qaCTUz67yQ+zvuYs1qfxaQBk0hP9G9aWMJPy1/9iv05H5B7zz1EN00hqsn3N22U5efz7UMP0/is/qSOq9ndsUeJiYOLHoR+/+sZqVsxAz5+AfrdAJljoGkHqC7h2r8L1mTDe49B0bfQeTAM/g2kn35UtdjoWKaeO5WhHYfy4PIHGfXGKEb3GM2NfW48ahHmw2WH+Sr/K+547w5y9+fy1KCnOL99zZZ+EWnIGlTiZ2bNgZeBjsAm4Ern3N5KdTKAPwLJQBnwkHPuZW9ZJ2A2kAp8BPzMOVe3C0+J1IOmI0eyf9l/2T19OrunTz+mvNWEu2nSoyeTXA9SG6UyfeV0XvvyNeKi4miT2Ia2SW1pntCchV8vJD4mngcGPsDwzsM1MtLAWWws6Y8/xtdZI9h83bFbkkU3a0abyVO+35otGFLawrBpMPAWWPKIZ/Ru2e8gLsmz80dad2jVA2ISIO9zyFvr+dm/03N8+7Ng5HPQcWC1pxncfjCZJ2XyxIdP8JfVf2Hx5sX0admHrfu2sq1oG3kH8nA4kuKS+PNFf6ZvWt/gvUeRBqRBJX7AXcDbzrnJZnaX93nl1TYPANc65740szbAR2b2pnMuH5gCPOmcm21mzwJj8SSJIicUMyN96hT2jxyJKz36uqfolBQaZWR8V+/mvjdzTvo5fLHnC7YWeb4kt+7byupdqxnUbhB397ubFo1a1MO7kLoQ164dJ897nYNffnlMWUL3U4ltlVY3J07t7Nnm7dw74JulsGOtJ9H7fB58PMtTJ7YxtOwGp1zsSQjbZkK7ftWPDB4hOS6Zewfcy9BOQ5mcM5mc3BzSE9Pp17ofbRPbkp6UzpknnclJTU6qm/co0gA0tMRvODDI+3gWsIRKiZ9zbv0Rj7ebWR7Q0swKgMHAqCOOvxclfnKCsrg4Es+ufpSkQkZaBhlpGXUbkISN2PR0YtPracq+Zdejt1RzDvZ9C6XF0LRjjZZlqax/6/7MzZpb63ZETkQNLfFr5ZzL9T7+FmhVXWUzOxOIA77CM72b75wr9RZvBXSxkohIfTKD5Nb1HYVIxLBAt/Gpa2a2CPA1Dj8RmOWca3pE3b3OuWZVtNMaz4jgaOfccjNrASx3znXxlrcD/uWc61nF8TcANwC0atXq9Nmzg7ceVVFREYmJiUFrT4JL/RO+1DfhTf0TvtQ34S3Y/XP++ed/5JzL9FUWdiN+zrkLqyozsx1m1to5l+tN7I5dqMxTLxl4A5jonFvufXk30NTMYryjfm2BbdXEMR2YDpCZmekGDRpUo/fjy5IlSwhmexJc6p/wpb4Jb+qf8KW+CW+h7J+Gto7fPKBiM8XRwOuVK5hZHJANvOCce7XidecZ2nwHGFnd8SIiIiInqoaW+E0GhpjZl8CF3ueYWaaZzfDWuRI4F7jOzD71/mR4y+4EbjOzDXiu+ZsZ0uhFRERE6lHYTfVWxzm3G7jAx+sfAuO8j/8G/K2K4zcCZ9ZljCIiIiLhqqGN+ImIiIhIDSnxExEREYkQSvxEREREIoQSPxEREZEIocRPREREJEKE3c4d4cjMdgLfBLHJFsCuILYnwaX+CV/qm/Cm/glf6pvwFuz+6eCca+mrQIlfPTCzD6vaSkXqn/onfKlvwpv6J3ypb8JbKPtHU70iIiIiEUKJn4iIiEiEUOJXP6bXdwBSLfVP+FLfhDf1T/hS34S3kPWPrvETERERiRAa8RMRERGJEEr86pCZDTWzL8xsg5nd5aM83sxe9pbnmFnHeggzIvnRN7eZ2VozW2lmb5tZh/qIM1Idr3+OqPdjM3NmprsVQ8SfvjGzK71/P2vM7KVQxxjJ/Phsa29m75jZJ97Pt0vrI85IZGbPmVmema2uotzMbJq371aa2Wl1EYcSvzpiZtHAM8AlwKnANWZ2aqVqY4G9zrkuwJPAlNBGGZn87JtPgEznXG/gVWBqaKOMXH72D2aWBNwC5IQ2wsjlT9+Y2SnA3cBA51wP4NZQxxmp/Pzb+Q3winOuL3A18IfQRhnRngeGVlN+CXCK9+cG4I91EYQSv7pzJrDBObfROXcYmA0Mr1RnODDL+/hV4AIzsxDGGKmO2zfOuXeccwe8T5cDbUMcYyTz528H4AE8/ywdDGVwEc6fvvkf4Bnn3F4A51xeiGOMZP70jwOSvY9TgO0hjC+iOefeA/ZUU2U48ILzWA40NbPWwY5DiV/dSQe2HPF8q/c1n3Wcc6VAAZAakugimz99c6SxwL/qNCI50nH7xzsF0s4590YoAxO//na6Al3NbJmZLTez6kY4JLj86Z97gZ+a2VZgAXBzaEITPwT63VQjMcFuUOREYmY/BTKB8+o7FvEwsyjgCeC6eg5FfIvBM1U1CM9I+Xtm1ss5l1+fQcl3rgGed849bmZnAX81s57OufL6DkxCQyN+dWcb0O6I5229r/msY2YxeIbdd4ckusjmT99gZhcCE4FhzrlDIYpNjt8/SUBPYImZbQL6A/N0g0dI+PO3sxWY55wrcc59DazHkwhK3fOnf8YCrwA4594HEvDsEyv1z6/vptpS4ld3VgCnmFknM4vDcxHtvEp15gGjvY9HAoudFlYMheP2jZn1Bf6EJ+nTNUqhVW3/OOcKnHMtnHMdnXMd8VyDOcw592H9hBtR/Plcm4tntA8za4Fn6ndjCGOMZP70z2bgAgAz644n8dsZ0iilKvOAa7139/YHCpxzucE+iaZ664hzrtTMfgm8CUQDzznn1pjZ/cCHzrl5wEw8w+wb8FzweXX9RRw5/OybR4FEYI73fpvNzrlh9RZ0BPGzf6Qe+Nk3bwIXmdlaoAz4P+ecZjJCwM/+uR34s5n9Gs+NHtdpwCE0zOzveP4pauG9xnISEAvgnHsWzzWXlwIbgAPA9XUSh/pbREREJDJoqldEREQkQijxExEREYkQSvxEREREIoQSPxEREZEIocRPREREJEIo8RMRERGJEEr8RETCgJl9ambOu2PMka9v874+rb5iE5EThxI/EZF6ZmYdgD5APrDkiNdbAW28Tz8OeWAicsJR4iciUv+Ge38vcM6VHvH66Uc8VuInIrWmxE9EpP5VJH6vV3r9NO/vg8Da0IUjIicqJX4iIjVkZgO91985M7uyijr9zKzIW+dRH+XNgHOBw8C/KhVXJH6rKo0EiojUiBI/EZEacs4tA+Z5n95nZtFHlpvZD4A3gCbALOAOH81cBsQAi51z+yqVVUz1appXRIJCiZ+ISO3cDZQB3YCfVrxoZm2AN4FUYD4wzjnnfBzvc5rXzFKB9t6nHwU5ZhGJUEr8RERqwTm3Fnje+3SSmcWaWVNgIdABWApc6Wuq1szigaGA4/uRwwqnHfFYI34iEhRK/EREam8SUAx0Am7CM3rXC1gFXO6cK67iuMFAIvChc257pbKKxK/E246ISK0p8RMRqSXn3DagYoHlJ/HcrLEJuNg5l1/NoVne35Xv5oXvE781zrnDtY9SRESJn4hIsEwDyr2P9wAXOedyq6psZgZc7n0610cV3dghIkGnxE9EpJbMLAb4E99/pjbGM/VbnTOB1sBXzrk1ldpLAU72PlXiJyJBo8RPRKQWvCN3M4AfAjuBr4EE4L7jHFrVos0AfQHzPlbiJyJBo8RPRKR2pgKjgSI8a/JN9L4+2sxOrea4LO/v6q7vKwM+C0KMIiKAEj8RkRozs/HAeDx33v7IObcCmA2sBKKBR6o47hSgO7ALWOajSsX1feuccweCHbeIRC4lfiIiNWBm1+IZ7XPAdc65twC8izT/1lttmJkN9HF4xTTvfOdcmY/y77ZqM7PE4/zEBPFticgJznwvJC8iIlUxs0vxTNHGAL92zj3lo85yoB+wzDl3dqWy/wBnAyOcc3MrlTUBCvH/H/Nj2hARqYpG/EREAmBmZwFz8CR9U3wlfV4V1/oNNLOKET7MrCUwAM9dv//2cVwGgX02fxpAXRGJcBrxExEJITMbA8wE/umcG1bf8YhIZNGIn4hIaFWM/s2tzyBEJDLpomARkdBaBnyC72VcRETqlKZ6RURERCKEpnpFREREIoQSPxEREZEIocRPREREJEIo8RMRERGJEEr8RERERCKEEj8RERGRCKHET0RERCRCKPETERERiRD/DyfW/zBdTzXnAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(10,5))\n", "plt.plot(y[1:-1], v.T[0])\n", "plt.plot(y[1:-1], v.T[1])\n", "plt.plot(y[1:-1], v.T[2])\n", "plt.plot(y[1:-1], v.T[3])\n", "plt.title('Eigenstates', fontsize=30)\n", "plt.ylabel('$\\psi$', fontsize=25)\n", "plt.xlabel('$x/L$', fontsize=25)\n", "plt.grid()\n", "plt.savefig('v3p2.png', dpi=200)\n", "\n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, '$mL^2 E/\\\\hbar^2$')" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZkAAAD4CAYAAAA+epuFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAW3UlEQVR4nO3dfbAldX3n8fdHiCaCCoYbCmeAQWtAgYoD3CWsoCUS5SklGneR2fVx3YyusMGHLRfNVjRJUZCsyIZKJDsiq5QIQZDIClHRlaCJKDM44dFZBhhkJiMzPgRctDDgd//ovnBmuDPcO/f0Offc835Vdd0+v+4+/T0FxYdf969/napCkqQuPGPYBUiSFi5DRpLUGUNGktQZQ0aS1BlDRpLUmV2HXcB8stdee9WSJUuGXYYkjZTVq1f/sKompttmyPRYsmQJq1atGnYZkjRSkty/vW1eLpMkdcaQkSR1xpCRJHXGkJEkdcaQkSR1xpCRJHXGkJEkdcaQkSR1Zt6FTJKLk2xOcntP218nWdMu65OsaduXJPl5z7a/6jnmiCS3JVmX5IIkGcLPkaSxNh+f+P8U8BfAJVMNVfXGqfUk5wEP9ex/T1Utm+Z7LgR+D/g2cB1wAvC3/S9Xkrq15KxrOz/H+nNP7uR7511PpqpuBH483ba2N3IqcNmOviPJPsBzq+qmal79eQnwuj6XKkl6GvMuZJ7Gy4EHq+runrYDknw3yd8leXnbtgjY0LPPhrbtKZKsSLIqyaotW7Z0U7UkjalRC5nlbN2L2QTsV1WHAe8DPpvkubP5wqpaWVWTVTU5MTHtJKKSpJ00H+/JTCvJrsDvAkdMtVXVo8Cj7frqJPcABwIbgcU9hy9u2yRJAzRKPZnfBr5XVU9cBksykWSXdv2FwFLg3qraBDyc5Kj2Ps5bgC8Mo2hJGmfzLmSSXAZ8CzgoyYYk72g3ncZTb/i/Ari1HdJ8JfCuqpoaNPBu4CJgHXAPjiyTpIGbd5fLqmr5dtrfNk3bVcBV29l/FXBoX4uTJM3KvOvJSJIWDkNGktQZQ0aS1BlDRpLUGUNGktQZQ0aS1BlDRpLUGUNGktQZQ0aS1BlDRpLUGUNGktQZQ0aS1BlDRpLUGUNGktQZQ0aS1BlDRpLUGUNGktQZQ0aS1BlDRpLUGUNGktSZeRcySS5OsjnJ7T1tH0myMcmadjmpZ9sHk6xLsjbJ8T3tJ7Rt65KcNejfIUmahyEDfAo4YZr286tqWbtcB5DkYOA04JD2mI8n2SXJLsBfAicCBwPL230lSQO067AL2FZV3ZhkyQx3PwW4vKoeBe5Lsg44st22rqruBUhyebvvnf2uV5K0ffOxJ7M9ZyS5tb2ctmfbtgh4oGefDW3b9tqfIsmKJKuSrNqyZUsXdUvS2BqVkLkQeBGwDNgEnNevL66qlVU1WVWTExMT/fpaSRLz8HLZdKrqwan1JJ8Avth+3Ajs27Pr4raNHbRLkgZkJHoySfbp+fh6YGrk2TXAaUmeleQAYCnwHeBmYGmSA5I8k2ZwwDWDrFmSNA97MkkuA14J7JVkA/Bh4JVJlgEFrAfeCVBVdyS5guaG/mPA6VX1ePs9ZwBfBnYBLq6qOwb7SyRJ8y5kqmr5NM2f3MH+ZwNnT9N+HXBdH0uTJM3SSFwukySNpnnXk5Gk+WjJWdd2fo71557c+TkGzZ6MJKkzhowkqTOGjCSpM4aMJKkzhowkqTOGjCSpM4aMJKkzhowkqTOGjCSpM4aMJKkzhowkqTOGjCSpM4aMJKkzhowkqTOGjCSpM4aMJKkzhowkqTOGjCSpM/MuZJJcnGRzktt72v57ku8luTXJ1Un2aNuXJPl5kjXt8lc9xxyR5LYk65JckCRD+DmSNNZmFDJJ9khycpKXbfsf6yS7JfnDPtb0KeCEbdquBw6tqt8E/i/wwZ5t91TVsnZ5V0/7hcDvAUvbZdvvlCR17GlDJskhwF3AF4BvAjcn2b9nl92BD/eroKq6EfjxNm1fqarH2o83AYufpuZ9gOdW1U1VVcAlwOv6VaMkaWZm0pM5B/gW8DxgEXAv8PdJlnZZ2A78B+Bvez4fkOS7Sf4uycvbtkXAhp59NrRtkqQB2nUG+xwFHFtVjwCPAKcm+RhwQ5JjgYe6LLBXkj8AHgMubZs2AftV1Y+SHAH8Tdvzms13rgBWAOy33379LFeSxt5MejLPAqq3oareB1wB3AC8pP9lPVWStwG/A/z79hIYVfVoVf2oXV8N3AMcCGxk60tqi9u2p6iqlVU1WVWTExMTHf4CSRo/MwmZtcDkto1V9V7gczT3ajqV5ATgA8Brq+pnPe0TSXZp119Ic4P/3qraBDyc5Kh2oMJbBlGnJGlrMwmZq4Hl022oqjOBzwB9Gx6c5DKae0AHJdmQ5B3AXwDPAa7fZqjyK4Bbk6wBrgTeVVVTgwbeDVwErKPp4fTex5EkDcDT3pOpqnNobv5vb/vpwOn9Kqiqpgu0T25n36uAq7azbRVwaL/qkiTN3rx7GFOStHAYMpKkzhgykqTOGDKSpM7M5GHM7UoyCRzXflwPrKmqtXMtSpK0MOx0yLRDi1dOfaR9YDPJI8Aa4BZgNXBLVd0xtzIlSaNoLj2Z/9r+PY8mVPYBlgGHAy8Djmm3F7DLHM4jSRpRswqZJB8CllXVqcB+wOeq6gPT7PdsmsA5AjisD3VKEkvOurbzc6w/9+TOzzFOZtuT2Qd4Q7u+GXhgup3aqV/+oV0kSWNqLqPLvoC9FEnSDsw4ZJJ8nHaizCS/CvwR8JIkx3dUmyRpxM2mJ/Mu4Mh2/WHgKzQTT/5Nkvcn+fV+FydJGm2zCZnfAr7crt9M896WY2jeN/NnwKYkq5L8zyTvTDKZ5Jn9LVeSNEpmfOO/qm5OcjmwV1UdneQZNC8sm6QZRXYETw5hhmbo8mNJ7qiqw6f7TknSwjar0WVVdQlwSbv+S+COdvk0QBs8B9MEzmS7/GYf65UkjZDZPifz5zTvb/nG1CuQe7XBc3u79AaPJGkMzTYATge+DjyY5KIkJz/dfZc2eCRJY2i2IfMC4D/RzEn2JuAa4IdJLk/yxiTP6XeBkqTRNauQqarNVbWyqk4EJoA304w4Owm4DNic5ItJ3pFkov/lSpJGyU7fL6mqn1bVZ6vq39IEzinA5TRDnT8B/FOSG5OcmWT//pQrSRolfbkpX1WPVtX/rqq3A3vTvGPmQmB/4Hzg3iSr+3EuSdLo6PvIr6r6ZVV9vap+v6r2p5kl4E+BX5vJ8UkuTrI5ye09bc9Pcn2Su9u/e7btSXJBknVJbk1yeM8xb233vzvJW/v8MyVJMzCrkEnyB0lePJtjqmpVVX2oqg6e4SGfAk7Ypu0s4GtVtRT4WvsZ4ERgabusoOk9keT5wIdpLt0dCXx4KpgkSYMz257MnwCn9jYkmVEPZaaq6kbgx9s0n0L73E3793U97ZdU4yZgjyT7AMcD11fVj6vqJ8D1PDW4JEkd68flsg8k2TzdhiT7JNm9D+fYu6o2tes/oLnvA7CIrd9ps6Ft2177dDWuaOdcW7Vly5Y+lCpJmtKvezLbm4F5BfCTPp0DgHamgafMNjCH71tZVZNVNTkx4ahrSeqnQUz50o9zPNheBqP9O9Vz2gjs27Pf4rZte+2SpAEalXnFrgGmRoi9leatnFPtb2lHmR0FPNReVvsy8Joke7Y3/F/Dk68pkCQNyKwmyByEJJcBrwT2SrKBZpTYucAVSd4B3M+Tgw+uo5ltYB3wM+DtAFX14yR/QvPeG4A/rqptBxNIkjq2MyHTt/sh03551fLtbDpumn2LZtLO6b7nYuDiPpYmSZqlnQmZ/5bkFJpews3Afv0tSZK0UMw2ZL5K8+bLqWXF1IYkNwJrgH9sl9v6U6IkaVTN9s2YrwFIcgBPvvlykiZwjmmXqctpjwOP9K1SSdLI2akb/1V1H3Af8LmptiRL2Tp4DgOeR8f3cCRJ81ffRpdV1d3A3TTvlSFJgBfTBI4kaQx1NoS5Hfl1V7tIksbQqDyMKUkaQTMKmSR7JDk5ycvay2C923ZL8ofdlCdJGmVPe7ksySE0Q5cnaELpliRvqKr72112p3kq/487q1LSvLHkrGs7P8f6c0/u/BwajJn0ZM4BvkUzUmwRcC/w9+1oMkmStmsmN/6PAo6tqkdonns5NcnHgBuSHAs81GWBkqTRNZOQeRbbPOtSVe9r783cAPy7DuqSJC0AMwmZtTTPutzZ21hV703yDJ6cdl+SpK3M5J7M1cC0MyNX1ZnAZ4BMt12SNN6eNmSq6pyqOnEH20+vKp+3kSQ9heEgSeqMISNJ6kynIZPknCTdP7klSZqXuu7J7AWc0PE5JEnzlJfLJEmdGZmQSXJQkjU9y8NJ3pPkI0k29rSf1HPMB5OsS7I2yfHDrF+SxtGs3ieT5Bjglqr6WUf1bFdVrQWWtXXsAmykeYbn7cD5VfXR3v2THAycBhwCvAD4apIDq+rxQdYtSeNsti8tuxF4PMlaYDWwql3WVNXP+13cDhwH3FNV92/z5oFepwCXV9WjwH1J1gFH0kz2KUkagNmGzA3AYcDB7fKmtv3xJN+jCZyp8FnTnxKndRrta55bZyR5S3ve91fVT2hmjL6pZ58NbdtWkqwAVgDst99+nRUsSeNoVvdkqupVVbUncBDNxJjnA98Afg4cCrwNuAD4B+CnwBv7WSxAkmcCrwU+1zZdCLyI5lLaJuC82XxfVa2sqsmqmpyYmOhnqZI09mbbkwGgqu4G7gYun2pLchDNRJqTwBE0PZ7d2WYG5z44kea+0INtLQ/21PAJ4Ivtx43Avj3HLW7bJEkDslMhM532xvxa4FKA9lUAL6EJnH5aTs+lsiT7VNWm9uPrgdvb9WuAz7bvvnkBsBT4Tp9rkSTtwJxDpr18tVt7H+QJVVU0rwe4c9oDd+5cuwGvBt7Z0/xnSZbR9JjWT22rqjuSXNGe/zHgdEeWSdJg7XTIJNkX+DTwiuZjfkpzs/+WnuWuNmz6on07569v0/bmHex/NnB2v84vSZqdufRkPg68EniA5jLZ/sDLaUJnKlh+nuQfq+rouRQpSRpNcwmZlwM3A8dU1b/AE5ezDmuXw2nux/yruRYpSRpNcwmZR4EbpgIGnric9c12AZ64ZyNJGkNzmbvsqzTPy+xQVf1iDueQJI2wuYTM2cBxSY7sVzGSpIVlp0Omqu6keWbl6iTL20krJUl6wlyGMO8NvAv4DeAzwAVJbuDJuctWb/vsjKT+WHJW9y+cXX/uyZ2fQwvfXG78XwicBPwzzUOQ+wNvaJcCSHI/sKqqTp1TlZKkkTSXkHkVcBtwdFX9P4Ak+9MMW+5d3jDXIiVJo2kuIfM48KWpgAGoqvuB+4HPT7Ulcf58SRpTcxld9g3ghU+3U1V9fw7nkCSNsLmEzB8BJyU5pF/FSJIWlrmEzHKaBzKvT/LqPtUjSVpA5nJP5r/QjCIL8KUk99KEztQQ5tur6rG5lyhJGlVzHV12eM9yIM27XKZmYP5FkttohjC/e05VSpJG0k6HTFXdANww9TnJs4GXsnXwvJRmGLMhI0ljqJ+vX/4Z8K12AZ6YgfnQfp1DkjRa+hYy02lnYL6ly3NIkuavuYwukyRphwwZSVJnRipkkqxPcluSNUlWtW3PT3J9krvbv3u27UlyQZJ1SW5Ncvhwq5ek8TNSIdM6tqqWVdVk+/ks4GtVtRT4WvsZ4ERgabusoJk1WpI0QKMYMts6Bfh0u/5p4HU97ZdU4yZgjyT7DKE+SRpboxYyBXwlyeokK9q2vatqU7v+A2Dvdn0R8EDPsRvatq0kWZFkVZJVW7Zs6apuSRpLnQ5h7sAxVbUxyW/QzJn2vd6NVVVJajvHTquqVgIrASYnJ2d1rCRpx0aqJ1NVG9u/m4GrgSOBB6cug7V/N7e7bwT27Tl8cdsmSRqQkQmZJLslec7UOvAa4HbgGuCt7W5vBb7Qrl8DvKUdZXYU8FDPZTVJ0gCM0uWyvYGrk0BT92er6ktJbgauSPIOmrdyntrufx1wErAO+Bnw9sGXLEnjbWRCpqrupZlwc9v2HwHHTdNewOkDKE2StB0jEzLSfLPkrGs7P8f6c0/u/BxSl0bmnowkafQYMpKkzhgykqTOGDKSpM4YMpKkzhgykqTOGDKSpM4YMpKkzhgykqTOGDKSpM4YMpKkzhgykqTOGDKSpM4YMpKkzhgykqTOGDKSpM4YMpKkzhgykqTO+PpljTRfgSzNbyPTk0myb5KvJ7kzyR1JzmzbP5JkY5I17XJSzzEfTLIuydokxw+vekkaT6PUk3kMeH9V3ZLkOcDqJNe3286vqo/27pzkYOA04BDgBcBXkxxYVY8PtGpJGmMj05Opqk1VdUu7/lPgLmDRDg45Bbi8qh6tqvuAdcCR3VcqSZoyMiHTK8kS4DDg223TGUluTXJxkj3btkXAAz2HbWDHoSRJ6rORC5kkuwNXAe+pqoeBC4EXAcuATcB5s/y+FUlWJVm1ZcuWfpcrSWNtpEImya/QBMylVfV5gKp6sKoer6pfAp/gyUtiG4F9ew5f3LZtpapWVtVkVU1OTEx0+wMkacyMTMgkCfBJ4K6q+lhP+z49u70euL1dvwY4LcmzkhwALAW+M6h6JUmjNbrsaODNwG1J1rRtHwKWJ1kGFLAeeCdAVd2R5ArgTpqRaac7skySBmtkQqaqvglkmk3X7eCYs4GzOytKkrRDI3O5TJI0egwZSVJnDBlJUmcMGUlSZwwZSVJnRmZ0meYvp9uXtD32ZCRJnTFkJEmdMWQkSZ0xZCRJnTFkJEmdMWQkSZ0xZCRJnTFkJEmdMWQkSZ3xif8FwqfuJc1H9mQkSZ0xZCRJnTFkJEmdMWQkSZ0xZCRJnVnwIZPkhCRrk6xLctaw65GkcbKghzAn2QX4S+DVwAbg5iTXVNWdXZzPYcSStLWF3pM5ElhXVfdW1S+Ay4FThlyTJI2NVNWwa+hMkn8DnFBV/7H9/Gbgt6rqjJ59VgAr2o8HAWsHWOJewA8HeL75wt89Xsb1d8P4/Pb9q2piug0L+nLZTFTVSmDlMM6dZFVVTQ7j3MPk7x4v4/q7Ybx/+5SFfrlsI7Bvz+fFbZskaQAWesjcDCxNckCSZwKnAdcMuSZJGhsL+nJZVT2W5Azgy8AuwMVVdceQy+o1lMt084C/e7yM6++G8f7twAK/8S9JGq6FfrlMkjREhowkqTOGzJCM43Q3SfZN8vUkdya5I8mZw65pkJLskuS7Sb447FoGJckeSa5M8r0kdyX518OuaRCSvLf9d/z2JJcl+dVh1zQshswQ9Ex3cyJwMLA8ycHDrWogHgPeX1UHA0cBp4/J755yJnDXsIsYsD8HvlRVLwZeyhj8/iSLgN8HJqvqUJpBR6cNt6rhMWSGYyynu6mqTVV1S7v+U5r/4CwablWDkWQxcDJw0bBrGZQkzwNeAXwSoKp+UVX/PNSiBmdX4NeS7Ao8G/inIdczNIbMcCwCHuj5vIEx+Y/tlCRLgMOAbw+5lEH5H8AHgF8OuY5BOgDYAvyv9jLhRUl2G3ZRXauqjcBHge8Dm4CHquorw61qeAwZDVyS3YGrgPdU1cPDrqdrSX4H2FxVq4ddy4DtChwOXFhVhwGPAAv+/mOSPWmuTBwAvADYLcmbhlvV8BgywzG2090k+RWagLm0qj4/7HoG5GjgtUnW01wafVWSzwy3pIHYAGyoqqne6pU0obPQ/TZwX1Vtqap/AT4PvGzINQ2NITMcYzndTZLQXJ+/q6o+Nux6BqWqPlhVi6tqCc0/6/9TVQv+/2yr6gfAA0kOapuOAzp5l9M8833gqCTPbv+dP44xGPCwPQt6Wpn5agSmu+nK0cCbgduSrGnbPlRV1w2vJHXsPwOXtv8zdS/w9iHX07mq+naSK4FbaEZUfpcxnl7GaWUkSZ3xcpkkqTOGjCSpM4aMJKkzhowkqTOGjCSpM4aMJKkzhowkqTP/H7ZQbagBG7NWAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.bar(np.arange(0, 10, 1), w[0:10])\n", "plt.ylabel('$mL^2 E/\\hbar^2$', fontsize=20)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.10" } }, "nbformat": 4, "nbformat_minor": 4 }