Skip to content

Instantly share code, notes, and snippets.

@brandondube
Created January 23, 2019 01:57
Show Gist options
  • Save brandondube/ae6fbe729fb6ecfdcc6a8dd1097e8051 to your computer and use it in GitHub Desktop.
Save brandondube/ae6fbe729fb6ecfdcc6a8dd1097e8051 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"\n",
"from prysm import FringeZernike, PSF\n",
"from prysm.thinlens import (\n",
" object_to_image_dist,\n",
" mag_to_object_dist,\n",
" image_displacement_to_defocus\n",
")\n",
"\n",
"from matplotlib import pyplot as plt\n",
"\n",
"plt.style.use('bmh')\n",
"\n",
"\n",
"def one_over_e(psf):\n",
" \"\"\"1/e radius of a PSF, simplified in 1D and with low accuracy.\"\"\"\n",
" u, d = psf.slice_x\n",
" l = len(u) // 2\n",
" u, d = u[l:], d[l:] # positive side only\n",
" diff = abs(d - 1/np.e)\n",
" idx = np.argmin(diff)\n",
" return u[idx]"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array([ 3981.07170553, 4085.30503983, 4192.26743524, 4302.03034466,\n",
" 4414.66709181, 4530.25292015, 4648.86504321, 4770.58269614,\n",
" 4895.48718863, 5023.66195925, 5155.19263117, 5290.16706937,\n",
" 5428.67543932, 5570.81026725, 5716.66650191, 5866.34157803,\n",
" 6019.93548138, 6177.55081561, 6339.29287075, 6505.26969354,\n",
" 6675.59215965, 6850.37404771, 7029.73211533, 7213.78617709,\n",
" 7402.6591846 , 7596.47730859, 7795.37002325, 7999.47019267,\n",
" 8208.91415964, 8423.84183668, 8644.39679955, 8870.72638315,\n",
" 9102.98177992, 9341.31814086, 9585.89467917, 9836.87477662,\n",
" 10094.42609266, 10358.72067645, 10629.93508176, 10908.25048496,\n",
" 11193.85280601, 11486.93283267, 11787.68634794, 12096.31426086,\n",
" 12413.02274072, 12738.02335479, 13071.53320962, 13413.77509612,\n",
" 13764.97763834, 14125.37544623])"
]
},
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"distances = np.logspace(3.6, 4.15, 50) # 5.6 to 21 meters, log spacing\n",
"lens1 = dict(efl=400, fno=2.8, name='400/2.8')\n",
"lens2 = dict(efl=600, fno=4, name='600/4')\n",
"lenses = (lens1, lens2)\n",
"distances"
]
},
{
"cell_type": "code",
"execution_count": 28,
"metadata": {},
"outputs": [],
"source": [
"curves = {}\n",
"for lens in lenses:\n",
" efl, fno = lens['efl'], lens['fno']\n",
" focusing_distance = mag_to_object_dist(efl, 1/15) # 1:15 is about headshot magnification\n",
" nominal_displacement = object_to_image_dist(efl, focusing_distance)\n",
" displacements = object_to_image_dist(efl, distances)\n",
" deltas = (displacements - nominal_displacement) * 1e3 # mm to um\n",
" wavefront_errors = image_displacement_to_defocus(deltas, fno, wavelength=0.55, zernike=True)\n",
" psf_sizes = []\n",
" for error in wavefront_errors:\n",
" p = FringeZernike(Z4=error, dia=efl/fno, samples=1024, mask_target='fcn', opd_unit='um') # mask target='fcn' will just accelerate computation\n",
" ps = PSF.from_pupil(p, efl=efl, Q=2)\n",
" size = one_over_e(ps)\n",
" psf_sizes.append(size)\n",
" \n",
" curves[lens['name']] = psf_sizes"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEXCAYAAACzhgONAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsnXd4HNW1wH93m7SqVrNkyw2DC2DAYIMBG2wwmBoIhARCIBQnEAIhhCSPEF5CCSWEl0ASSgqEEkhMCAm9mWpjwOBecbfc1MuuymrrfX/MzGolbZmVtmk1v+/bT7szszP3aGbnzCn3HCGlxMDAwMDAIBKmdA/AwMDAwCCzMRSFgYGBgUFUDEVhYGBgYBAVQ1EYGBgYGETFUBQGBgYGBlExFIWBgYGBQVQMRWFgkMEIIaQQ4rKQz0VCiP8KIZzquglCiA+FEI+nYCxXCiF8yT6OQeZhKAqDlCCEOFkI8bIQoka9wf1vmG2uVNf1fZ022H0PYUYB/w75fB1wAjBbXbcXuBC4OVEHFEKMUf+P8/qseh6oTtRxDIYOlnQPwGDYUABsAv4BPBRlOz8wps+ylgTte8ghpazrs2gSsFFKuT5kWaz/T6LG4gJcqTiWQWZhWBQGKUFK+YaU8lYp5fOAO8a2dX1enkTtOxT1qfkHQojnhRCdQog9QoiLhBDFQojnhBDtQoidQoiv9fnePUKIzUKILiHEXiHEn4QQxX22+aYQYocQolsI8YkQ4lz1eHPU9fPUz6cLIZao+9okhDgjzBgvU9/vBhYCp6rLP1SX93M9CSGuV/fnFkI0CCH+HbLuUiHEciGEQwjRJIR4XQgxOeTre9W/H6jH2a1+r5/rSQhxthBiZchxHhVC5Iesf0oI8a4Q4hrV4nOq1l+F3vNkkH4MRWGQaZjVm3OtegM8N3Sl6pOXQogrE3S824A3gKOA14BngEXAYuBo4HXgGSFEWch3XMA1wGHAlcA84A8hY5wBPAf8U93vb4hs6fwfcK+63QrgeSHEiAjbHgv8C1iK4na6MNxGQog7gfuBR4EjgDOBNSGb5AC/Ao4BTkex4l4XQtjU9ceof7+mHufYCMc5EngFWAJMB64AzgX+FGbcpwDnqGOZrsptMFSQUhov45XSF7Ab+N8wy08Avo1yIzkB5eYqgYUh21QDXwIXxLPvCNtK4KGQzxXqsj+GLCtRl50bZT8XoFgyJvXzc8DSPtt8T93PHPXzPPXzhSHbVKnLzugzxstCPj8FvNtn3x8Cj6vv81EU2U/iOB+l6nFmq5/HqJ/n9dnuSsAX8vnvwOd9tjkfCADjQ8bbCOSEbPMzoDbd16Hx0v8yYhQGGYOU8lPg05BFnwohSoFbgCfUbfYDUxN42LUhx28UQviBdSHLWoUQHmCktkwIcSFwE3AIUIRimdtQbvQHUCyNd/sc51PCE3zSl1LWqcevHIQ8hwO5wDuRNhBCTAduR1HI5YBQV40HlsV5rPf7LPtI3d9hQI26bLOUMtQluJ/ByWiQYgzXk0Gm8wkwIYn79+pYJlF/K0KIWcALKO6WC1DcNN9Tt7P1+Y4ewsVfEvG7DHt8IUQeihKRwNXAcSiuIUnv8Q/qOH2W95VR0qOcDIYAhkVhkOkcTU9wNROYAzRJKYMpuEKIi/psswnFdRbK8ckeWMixu4EzgPVh1h+K4mK7TUq5GUAIcSK9b9zajd0c41gbgbl9ls1FUQSb4hu2QSZjKAqDlCCEKEBx1YDqplFdIB1Syu3qNncAnwNbUQKuFwHfAW4M2U818B5wq5Tyv3r3nUC2ABVCiIXAByiK4/t9tvkd8IUQ4i7gWRRX2Y/VdUltACOl7BBC/Ba4QwjhQgnK24GzpZT3obiD3MAP1O0mAL/uM64moANYIITYCLillK1hDvcAsEoI8TvgL+q+/gg8J6Xckwz5DNKD4XoySBUzgdXqaxRwvfo+NK2zCHgE5Ul4KcpT8TeklI+EbGMFpgCh6ah69p0QpJSvAfegZCqtBy4Bftpnm5XAt9TXeuBWQLNAuhM9pjD8AiWb60ZgA4qr6Rh1bE3AZSjZThtRso9+ghKA1sYfQPkffgPFmlsd7iBSynXAeShWxFqU4Pbr9LjiDLIEIaXR4c7AINkIIb4NPAmUSSnb0j0eA4N4MFxPBgZJQAjxExTXVAtKsPh+4AVDSRgMRQxFYWCQHI5EiUuUorhvnkVJSTUwGHIYricDAwMDg6gYwWwDAwMDg6gYisLAwMDAICpZGaP48MMPZU5OTrqHYWBgYDCk6Orqapo/f36/yr5ZqShycnKYOnXg5YBqamoYP358AkeU2RjyZjeGvNlNIuVdtWpVTbjlhuspDFarNd1DSCmGvNmNIW92kwp5DUURhuLi4tgbZRGGvNmNIW92kwp5DUURhqampnQPIaUY8mY3hrzZTSrkzcoYxWAxnkiyG0Pe7KaoqIj29naGyxyx/Px8nE5nXN8RQlBQUIAQ+qq9G4oiDB5P1BbNWYchb3Yz3OTt7u6mrKwMm20g7TWGHl6vN+44hcfjoaOjg8LCQl3bG66nMLhcrnQPIaUY8mY3w01eKeWwURIAgUAg9kZ9sNlscVlchqIIQ1VVVbqHkFL6ytvS5cXZ7UvTaJLPcD+/2Y7dbk/3EFKKkfWUJurq6tI9hJQSKq/L6+fa/3zJDS9vweuP/0llKDCcz+9woLs7FS0/ouNwOLjiiiuYNWsWs2bN4vPPP6e1tZULLriAmTNncsEFF9DWphQSllLys5/9jBkzZjBnzhzWrl3ba18XXXQRBw4cCH6+5ZZbGDt2bPCz1xuum29iMRRFGIaT2Qq95d3a2IWj20ddu4ePd2dnRezhfH6HA3oDtMnk1ltvZf78+SxfvpylS5cyZcoUHnroIebOncuKFSuYO3cuDz30EADvvvsuO3bsYMWKFTz44IP8+Mc/Du7H5XLR1tbG6NGjAVi9ejUOh6PXsVIhr6EowqA3wJMthMq7qaEz+P7ljdmZZjicz+9wIN0T7pxOJ5988gmXX345oCjq4uJi3nzzTS655BIALrnkEt544w0A3njjDS655BKEEBx77LE4nc6gFfjxxx8ze/ZsAPx+P7fffjt33HFHr+OZzbFamw8eI+spDM3NzRQUFKR7GCkjVN7NIYpiU0Mn25u6OKQ8L11DSwrD+fwOB0KzvBY8HraL66B55ztHR1xXU1NDeXk5N9xwAxs2bOCoo47ivvvuo6GhIRgvqqqqorGxEYDa2lqqq6uD3x89ejS1tbVUVVXx3nvvcfbZZwPw17/+lTPPPLNfzMnn8yVdWRgWRRhKSkrSPYSUoskrpWRzQxcAx48rAuCVTdlnVQzX8ztcSLerzefzsXbtWq666io++ugj8vLygm6mcITLPtLcScuXL+f444+ntraWl19+mWuuuabfthZL8p/3DYsiDC6Xi6KionQPI2Vo8h5wenB0+yjOtXDNrGo+2+Pk/R0tfOe40RTlZs+lMlzP73DB7/cH30d78k8Wo0ePZvTo0cycOROA888/n4ceeoiRI0dSV1dHVVUVdXV1VFRUBLffv39/8PsHDhygqqqK3bt3U11djc1mY/369ezatYsZM2YA0NXVxYwZM1i5ciWBQMCwKNJBJmRNpBJNXs3tdNjIfMYU5zJzTCEev+Sdrc3pHF7CGa7nd7gQqijSQWVlJdXV1Wzbtg2Ajz76iClTpnDmmWeyaNEiABYtWsRZZ50FwFlnncWiRYuQUvLFF19QVFREVVUV7777LvPnzwdgwYIFfPnll6xdu5a1a9eSl5fHypUrgYHNo4iX7HlMTCDDLe9ck1cLZB9aqcQkzjusghX72nl1cxMXTBuJ2ZT+bJJEMFzP73AhNzc33UPg/vvv59prr8Xj8TBhwgQefvhhAoEAV199Nc8++yxjxozhySefBOD0009n8eLFzJgxA7vdzsMPPwzAe++9x/333x/zWKkI3huKIgx1dXXDqp69Jm+oRQFw7JgiKgts1LZ7WLHPyaxx2VEzaLie3+FCJlhQRxxxBO+//36/5S+99FK/ZUIIHnjggV7L3G439fX1jBs3Luz+9+7dG3zv9XpJdqM2w/UUhkx4Ikklubm5uLx+drW4MAmYpGY5mU2CrxxWDmRXUHs4nt/hRCrSRZNNTk5OWEUTDpMp+bdxQ1GEYbiVALDb7Wxt7CIgYWKpHbu154d25uQybGbBF/uc7He40zjKxDEcz+9wIhsURTwYiiJNtLa2pnsIKaW1tTUYnzisMr/XuqJcC6ccrKRXvrq5MeVjSwbD8fwOJ4ZbtVyfL/l12QxFEYaysrJ0DyGllJWVBeMTh47M77f+9EmlAKyr7UjpuJLFcDy/w4l0z6NINamYR2EoijC0t7enewgxWXOgnYUvbGJTfWfsjWPgdDqDE+0OC6MoCnOUC9Hrz45GMEPh/CaS4SZvKorkZRKpSAdOiaIQQvxNCNEghNgQsqxUCLFYCLFN/VuiLhdCiD8IIbYLIdYJIY4J+c4V6vbbhBBXJGu8Q8F0Xb7HwV6Hmw92DN6tsN/RHZxoV1XY/2nMZlbSYr2B7FAUQ+H8JpLhJu9w6WynkQp5U2VRPAWc2WfZz4D3pJSTgPfUzwBnAZPU1zXAY6AoFuB2YBZwHHC7plwSzVDIO3f5lEk2O1sG35SmCcWKOGxkfthKlFazcplkS9nxoXB+E8lwkzdTsrz8fj9z584NFgKsqanhtNNOY+bMmVx99dVBBe52u7n66quZMWMGp512Gnv27Om1n1NOOQWn08nFF1/MrFmzOOGEE7jzzjuD60PnUXi9Xr7//e8ze/ZsZs2axYMPPpgQWVKiKKSUS4CWPovPB55W3z8NfDVk+TNS4TNghBBiFHAGsFhK2SKlbAUW01/5JIShUL/f5VVu2rtaXIN+olhdo8y81iba9cWqTrTLFtfTUDi/iWS4yZsJ8ygA/vSnPzF58uTg5zvuuIPrrruOFStWMGLECJ599lkAnn32WUaMGMHKlSu57rrrelWH3bNnD6NGjcJisXDDDTewfPlyPvroI5YvX87ixYuB3q62l19+GbfbzbJly/jggw946qmn+imegZDOCXeVUspaACllrRBipLq8Gtgbst0+dVmk5f1oaGhg4cKFWCwW/H4/F154Iddffz11dXXk5+djNptxOp1UVFTQ0tKClJKKigrq6+spKCjA7/dTU1NDZWUljY2NCCEoLS2lsbGRoqIi/H4/nZ2dwZotVquV4uJimpqaKC4uxuPx4HK5guttNhuFhYU0NzdTUlKCy+Wiu7s7uD43Nxe73U5raytlZWW0t7fj8XiC6+12OzabDYfDQXl5OQ6HgxaHElju8PjZsreeysKcqDIBdHR0hJVpR5uSNVEa6MDtdveTyWxXvu/2+YPrkyGT1+sNrtdznqLJFO08dXZ20tbWlpLzlCqZol17QghqamqySqZo58lsNuP3+/H5fMF7gJQSq9WK1+vFZDIhhMDv9+teD0rQ2Ov1BtNv/X4/Vqs1mHVkNpuDlVz379/P22+/zc0338yjjz6K2+1m6dKlPPzww/h8Pi666CIeeOABrrzySl577TV++tOf4vf7OfPMM7nlllvweDxIKXnnnXeYO3cuVquVE088EbfbjcViYdq0aezZs4dAIEAgEMDr9SKEwOfz0dnZSVdXFw6HA6vVis1mC64PlamzsxObzdbrPEVCpMqfJ4SYALwmpZymfm6TUo4IWd8qpSwRQrwO3Cel/Fhd/h7wP8CpQI6U8m51+S+ALinlb/se69NPP5VTp04d8Fjb2toYMWJE7A3TyE9f38ZaNQvprgUTOX6As6ZdXj9ffXodQsB/v31krzkUGt2+AOc9tRabWfDaVdMHNe5MYCic30Qy3ORtamqivFyZKPpW1YlJOcaZdZ9EXX/FFVfwox/9iI6ODh5++GEeeeQRFixYEKzPtG/fPr7xjW/wySefcOKJJ/LCCy8ES40fc8wxLF68mLKyMr71rW9xzz33MGHChOC+HQ4H8+bN47///S8TJkwIKkRQrIvvfe97LFmyBJfLxd13382VV14ZdoxOp7NfschVq1atnD9//sy+26Yz66ledSmh/m1Ql+8DxoZsNwY4EGV5wunbQSoT0VxPADubBx6n2NrYhaT/RLtQQl1P2RAoHArnN5EMN3nTnfX09ttvU1FRwfTpPQ9V0UqJR1rn8Xg4cOBALyXh8/n4zne+wzXXXBNcHpr1tHLlSsxmM5s2bWL16tU8+uij7N69e9AypdP19ApwBfBr9e/LIctvEEIsQglcO1TX1NvAvSEB7AXArckYmPY0ksm4vD0Xxy6dAW2X18+L6xuwmAVVBTlUFdpYuV9Jnew70S4Us0lgEhCQ4JdgGeK1AYfC+U0kw03e0LpHsZ78k8Hy5ct58803Wbx4MW63m/b2dn7+85/jcDiCT/9aKXHoKTNeXV2Nz+fD6XRSUlLCkiVLmDVrVq9933TTTRx88MFcd911wWWh8yhefPFF5s+fj9VqpaKiguOOO47Vq1f3UjYDIVXpsf8EPgWmCCH2CSEWoiiI04UQ24DT1c8AbwA7ge3AX4HvA0gpW4BfAV+or7vUZQlnKDyBaVlPADt0Kor3trfyzKo6/vZFLfd+sJsbX9nKorX1QPiJdqFkU+bTUDi/iWS4yZvudOBf/vKXbNy4kbVr1/L4449z0kkn8Ze//IU5c+bw8svK8/CiRYuCneu0MuOgBKNPOukkhBC89957nHbaacH93nPPPTidTu69995exwu1KMaMGcOSJUuQUtLZ2cmKFSt6BdQHSqqynr4ppRwlpbRKKcdIKZ+QUjZLKedLKSepf1vUbaWU8nop5cFSyiOklCtC9vM3KeUh6uvJZI033aarHkJdTwecbrp9sW/gNa2KQjmyqoDZ44s5uMxOntVEaa6JGdXR+yoH51JkQebTUDi/iWS4yZup7tE77riDxx57jBkzZtDS0sJll10GwGWXXUZLSwszZszgscce4/bbbwdg2bJlwX7Z+/fv57e//S1btmxh3rx5nHzyyTzzzDMAvPXWW0HlsXDhQjo7OznxxBOZP38+l156KYcffvigx26UGQ9DpuedSymDrqcxxTnsc7ipaXUxpSK6VbCnTSnq9/UjRwZLhkspcbvd5OZGr2kfjFNkwaS7TD+/iWa4yZsp8ygA5syZw5w5cwCYMGEC7777br9tcnNzeeqpp3ot279/P6WlpcGCjtXV1bS0hHegnHvuuZx33nkAFBQU9NtXIjBKeIQh0/POPX5JQCo3b60k+M6W2Lnjex3KNmNH9PyQhBDU19fH/K5e19PWpi5q2zO7ymymn99EM9zkzZR5FIOhurqaF154Qde2qbAYDUURhvz86E/m6UazJnKtJiaWKk8csTKfujx+mjq9WE2CyoLeZTr0yGvRMemu3e3jR69u5X9e344/gy2PTD+/iWa4yWuUGU/CMZJ+hCFIpl9oWiA7z2oOKopYmU/71F4S1cU5/Vqa6pHXqiNG4ez24fVL6js8bKzP3EqzmX5+E81wkzdcGZpsJhXyGooiDE6nM91DiIrLoyiKXhZFjFIee9r6u5009MgbVBRRGrl7QpTIhzvbYu4zXWT6+U00w03e4Ra8z5rqsUONioqKdA8hKi6fcmHYLSZK8ywU51ro8Php7Iz8A9HiE+PCKAo98lpNWowisjIKXffxrraMdT9l+vlNNMNN3tzc3LSnyKaSgfSj8Hg8cVkiRtZTGFpaWsjLC18gLxPQUmPtVjNCCA4qzWXNgQ52trgYWRC+acteNeNpbHH/Jux65NXjevKEBLrbun2sq+3g6Bhpt+kg089vohlu8nZ2dmK1WrMiqK2Hrq6uuM+vECJmfadQDEURhkzNw9boURTKU/7EUjtrDnSwq8UVseZTuIwnDT3y6nE99VUiH+5szUhFkennN9EMN3kBCgsz77pLFg6Ho1/NpkRjuJ7CkOmmupb1lBeiKCBy5pM/IDmgBrPHhLEodLmezLFdT5pFUZGvzMlYtrsNXwa6nzL9/CYaQ97sJhXyGooiDNq8gkCGPolpFkWuWsRPUxSRSnnUtbvxBiQV+dawhf/0zKOw6UiP1dZNKs9jbHEOTrefNQcyrw2nHnmzCUPe7CYV8hqKIgwFBQX84eO9XL5oI12e5GcUxEtoMBuUALVJRC7loc3IDud2AnT5Ki26sp6UdTazYO5EpXbjRzsH36o10cTjm80GDHmzm1TIayiKCKypbaex00tNW+YFxIIxCptiHdgsJsaOyCUge+o5hRKMTxQPvLSBnqwnLT3WajYxd6LS/2DZbkdWFBI0MBjOGIoiDB0dHUHfurPbl+bR9CeoKCw9p69nPkV/xba3TUuN7R+fAEXeWOjJevKGWBTjS+xMKMmlw+NndYa5n/TIm00Y8mY3qZDXUBRhqKysxKfeEJ3uTFQUvYPZED2gvTeG66mysjLmMXsURewJdzY18K25nzJt8p0eebMJQ97sJhXyGooiDI2NjcEqqY7uDIxR9AlmQ4+i2NTQ++lCShk1NRYUeWMRzHqKksWkxS80paK5nz7Z3dZrjkW60SNvNmHIm92kQl5DUYRBCBF8cs5o11OIRXHkqALsVhPbmlzsd/RUb23r9tHu9it9J+zhp83omaFp1ZH15PH1tijGFOdycJmdLm+Alfsyx/1k1ALKbgx5E4+hKMJQWloajFE4MlFR+Pq7nnIsJmaPVybbfRCSaRTqdop0QZWWlsY8Zs+Eu9gxCm1bgJMPUqyKDzMo+0mPvNmEIW92kwp5DUURhsbGxiERzM7tMydi3sFqTGBHa3A2biy3E8TpeooWowjIXttCT5zisz0OXV34UoHhmshuDHkTj6EowlBQWIj24Ox0Z26MIjTrCeCY6iKKcszsaetmpzr5Tst4ClfjSUPP9H89ridv0PXUY1GMLsphSkUeLm+AJRliVSS73EGmYcib3aRCXkNRhKHb02NFZKJF0R3MeuptUVhMgpPVJ/gPdig3Zc31FK5qrIaeMsW60mODwezel9VXDi0H4JVNTTGPkwpSUZY5kzDkzW6MMuNpwtnRGXyfiTGKLm9PP4q+nKK5n3a2EtCR8QRKtc1YxNOPItSiAMX9VJhjZmtTF182xD5WstEjbzZhyJvdpEJeQ1GEobR8ZPC90+3LqOqbASmDvv5cS//Td3hlPhX5Vho6vKw50E59uwezUFxAkaiqqop5XF0zs33ahLve48qxmDhjchkAr25Ov1WhR95swpA3u0mFvIaiCMOB+obg+4CEzgyq99StWhM5FlO/lqYAJiGYp7qf/r6qDgmMKsoJ9rwOR11dXczj6nM9yV7bhnLuoeUIFEsn3VaaHnmzCUPe7CYV8hqKIgyiT4/hTJp0p/XL7hvIDkVzP22sV0zSaG4nAKvVGvO4+lxPPSU8+jK6KIeZY4rw+iVvb22OebxkokfebMKQN7tJhbyGogiDPb9305NMKuMRDGTbIp+6g8vsvbKcxkXJeAIoLg7f7CgUPf0ovP7+6bGhfOUwJaj92uamtLZJ1SNvNmHIm92kQl5DUYShsbml1+d0u0pCCQayLf37SmgIIYJWBcS2KJqaYscN9PSjiGZRABw7pojKAht17R5W7nfGPGay0CNvNmHIm92kQl5DUYQhN693ffdMSpHV5lDkhcl4CiUeRRGXRaEr6yn82MwmkRGpssYTZ3ZjyJt40q4ohBA/EkJsFEJsEEL8UwiRK4Q4SAixXAixTQjxvBDCpm6bo37erq6fkIwxdbs9vT5nlqJQXE/hUmNDqS7OZc6EEYwtzgkWDIyEx+OJuh4IBsP1uZ4iB87PmFKG1Sz4Yq+TWqc74nbJRI+82YQhb3aTCnnTqiiEENXAjcBMKeU0wAxcAtwPPCilnAS0AgvVrywEWqWUhwAPqtslnM7u3jcwRwbNzu4pCBjZ9aTxy9MO4omvH0ZOlMA3gMsVvoVqKHqynnpcT5GPV5xrYe7EEiRKrCId6JE3mzDkzW5SIW/aLQrAAtiFEBYgD6gFTgX+ra5/Gviq+v589TPq+vkiCaUTi0aU9PqcURaFT5/rKR50zaPQkfUUtCiipOJCz0ztt7Y2405D/Scjzz67MeRNPLruNkKIYiHEL4QQ/xFCvBP6GszBpZT7gf8D9qAoCAewEmiTUmp3531Atfq+Gtirftenbl82mDGEI5OD2UHXU5Rgdrzom0ehJ+tJtShiWDBTK/KYVG6n3e1PS09tI88+uzHkTTzhGxT05wUUt9B/gYTZOUKIEhQr4SCgTT3OWWE21e5O4R5V+925GhoaWLhwIRaLBb/fz4UXXsj1119PXV0d+fn5mM1mnE4nFRUVtLS0IKWkoqKC+vp6CgoKgrWeCnPMtLv9NDm76OrqorGxkaKiIvx+P52dnVRVVVFXV4fVaqW4uJimpiaKi4vxeDy4XK7gepvNRmFhIc3NzZSUlOByueju7g6uz83NxW6309raSllZGe3t7Xg8nuB6u92OzWbD4XDg7FT+BZ6udtxut26ZQGmZWFlZSWNjI0IISktLlUq5Ph9NTU1RZXJ0KKVAPP4ANTU1YWXSrIPa/XspzLNHlWnemBy2Nbl4cc1+Zlfn4nA48Hq9wfWDlSnaeero6KCtrS2p56m8vDylMkW79qSU1NTUZJVM0c5Td3c3Tqczq2SKdp46OjpoaGhIiEwR79V6ylMIIZxAmZTSG3PjOBBCfB04U0q5UP38beAE4OtAlZTSJ4Q4AbhDSnmGEOJt9f2nqquqDqiQfYT49NNP5dSpUwc8rjc2HOChz+o5qCSXXa3djBuRy+MXHTrg/SWSxz/fz7/WNXDVzFF8c3piTM6Ojo6YF4rL6+f8p9eRYzHx6pVH9VsfkJIzn1gDwNsLp8dsptLtC/Ctf26g3e3n4fOnMLkib+ACxIkeebMJQ97sJpHyrlq1auX8+fNn9l2u19H9MZCMO+Ue4HghRJ4aa5gPbAI+AC5St7kCeFl9/4r6GXX9+32VRCJoa1faiZblKzMeM8v1pD+YrZfm5tgzpW0x+lGEZjzpCRvlWkwsmKQ0XHl1c2r7B+iRN5sYjLwBjxfX3lpaV6yn8d1PcDe2xP5SmolH3s23/54NP74vo+q5xUsqrme9rqcrgTeO2AnPAAAgAElEQVSEEMuB+tAVUsq7BnpwKeVyIcS/gVWAD1gN/AV4HVgkhLhbXfaE+pUngL8LIbYDLSgZUgknJzcP6KTUriiKdrePgJSYwtwAu/bUEvB4KDhkfDKG0o9kBLNLSkpibmM2CUxCqX3lD8h+daaC3e1iBLJDOffQcl7c0MgHO1q5ZlY1hTl6L8fBoUfebEKPvJ7mNvY+819c++rorm3CXd+Eu64RT3Nb7w2FoPiYwxi5YA4jF8yhYOrEjGs9qvf8eprbqPnz8wCM/+7FFE6dmMxhJY1UXM96f5n3AGOB3UBol4xBq2Ep5e3A7X0W7wSOC7NtN4pbKqm41HkUuVYT+TYznR4/HW4/Rbm9/10Bj5flX7kWX2cX81a9hLUo+eauy6NvHkVc+3S5dDU/sZgEHr/EG0ZRxJpsF47q4lxmVBeycn87b29t4aIjRsb+UgLQK2+2EEvegM/Hiktvxrn2y37rhNmMbWQpuZXlmHJzaFu1EcdK5bXtvj+TO6aKkafPpmLBbMpOPAZTji2ZouhC7/l1hMjbvPSLIasoUnE961UUlwCTpZS1yRxMpqApCotJUJyrKIp2t6+fomhethJ3vTIXoPWzNYxcMCf5YwsWBUyc66m7u1vXdlazCY/fj9cf6FfiPDiHwhLf0+VXDitn5f52XtvcxIXTKsJabYlGr7zZQix5dz36D5xrvyR39Egm/vAKckdVkFNZTk5VOTnlJb2KZPo6u2j+6AsaFi+jcfEyuvfVsefJF9nz5IuY8/MYd9WFTL7turRaGXrPr2PN5uD75o++YMJ3L07WkJJKKq5nvYpiJ5DQQHYmoxQFdGE1CdUd4sHR7ae6z0z5+tc/DL5vXrYyJYqiW2cJj3jQm4cdrR1qzxyK+MY1a2wxIwusHHC6WbW/nZljkv+kb+TZ99C+eQfbH3gcgGkP/pzyuf0M+V5Y8vOoPHsulWfPRQYCONZspuGdj2l8Zxntm7az6+FnEWYTk2/9XkJliAe95zdUUbR8spqAx4vJNvQqz2bMPArg78ArQohvCiFODX0lc3DpotXZDoDFbKJYtSL6BrSl30/Dm0uCn1s+XpWSsXXpLOERD3rzsKPNzo5VEDASZpPgnKnKBLxUNTUy8uwVAl4f6394N9LrY8zl58dUEn0RJhMjjjmcyT+7ltnvP8Mxz/wGYTaz8/fPUPP4C4kY+oDQc36llDhVRWEdUYi/y0Xbyg3JHlpSyKR+FNcDo4B7UQLK2uvxJI0rrQiTYmpbTSLobupbarx1+To8zW3Yx47ClGOjfeO2/oG/JNBTFDBxrqfc3OhFAzWiFQaMVWI8GmdOLsNiEizf46ChI/l1a/TKmy1EknfXw3/HuW4LuWOqmHr7DYM+zsgFczj8tz8DYPMvHqL2pXcHvc+BoOf8uuuacDc0YykuZNTXzgCgecmKZA8tKaTietb1q5ZSHhThNTSjP7FQfbIWs6A4R3nft4xH3esfAFB13qmMOPYIAFo+Sb5VobcoYDzY7dGLBmpEtyjC98vWQ0melZMOGkFAwutfJt+q0CtvthBO3vZN29n+uycBOOLBn2MpyE/IscZccg6Tb/seSMm6H9xF89LU33z1nF/Hmk0AFB81lfK5swBoWvJ5UseVLFJxPWdCraeMo7NLCQ71sihCFIUMBKh/4yMAKs85hbI5MwBoWZYCReFL/DyK1lZ9ZTSCMYowTYc019NALAroqf/05pfNEedqJAq98mYLfeUNeH2su/FXSK+PsVdcQNlJ/eZXDYqDbric8d/9BtLrY9VVP8OxbktC9x8LPedXi08UTz+U0hOnIyxmHKs343W0J3t4CScV17PeWk9HCSHeF0K0CCE86ssrhMjKer5mm9IRzhKiKELboTpWb8Jd20hudSXFRx9KqaoompetTOq4fAGJ1y8xCcgZwJN7JMrK9JXLijbpzjsIiwLg8Mp8DirJpa3bx8e7HQPah170ypst9JV35++fpn3DNuxjRzHll9cn/HhCCKbeeSOjLjgdf0cXKy+9ma7d+2J+z9/VzYH/vkPNE/9G+gdesVnP+Q1VFJaCfEbMmAaBQEq8AokmFdez3se/fwLLgJNRZmgfCkwlObO1047LrSR4WcwmitVJYI6QGEX966o1cfZchBAUH3Uo5jw7ndtq6K5L3izjnoKApoSmH7a363uKitaTwjtIi0IIwVcOqwDg1U3JnamtV95sIVRe54at7HjoKQCmPXQblvzklE4RJhNH/P5/KTv5WDxNrXxx8U1hZ3XLQIDmj1ew/od38/6R57LuujvYfNvv2PvMSwM+dqzzK6UMzhkpOkop9VN28rGAkiY71EjF9az3V10F/FJKuUFKuSP0lczBpQuPT7khh7qe2lXXk5QyGJ+oPGceACarhZLjpwPJdT8lo3wH6G98oidGEa1pUSxOPbgEq0mwob4zqBSTwXBtbBPweFl/491In59xV19E2exjknpck83K0X+7l6Ijp+KqOcDKS2/G194JQMeWXWy55zE+OvZrfHHRjex//g38HV3kT5oAwPb/ewJfR+eAjhvr/Hbt3o+3rR1bRSm5o5VJnmVzFUXRtGToKYpMalz0NHBpMgeSSVhylCwCq1lQlGsm39lG9Ysvsvb6O/js3Gtw1RzAVl5CiRrEBoJxiuaPw7ufAh4v2+7/66BS8LSbpz2BgWyIYx5FlJ4U3gGmx4aSZzMH05Hbk9gsarjOo9jx4FO0b9qOffxoJt92XUqObSnIZ8Zz/0feQWNwrt/KF9/4IZ8suIqP536LXX/8O93767GPHcXBP7qSOR//kzlLnmPEsUfgaW5j1yPPDeiYsc5vMJA9/dCgZV48/VAshfl07dyLa+/QmlecSfMofg38Sm1Z+n7oK5mDSxedXUoldYtJUJxjYdaHb3H4W69S++I7OFZuBGDMZef1mrFaGiOg3fr5OnY8+CSbf/H7AY+rx6JIrKLQP48ick+KgZTwCEeBmmXW7k5eIcbhOI/CuWErO//wDABHPHQblvzUZX7lVJQyc9GD2CpKcazehHPdFixFBYz51lc47qVHOXn5C0y65RoKDhmPEIIpv1RSdXf96Z8DcuXGOr+h8QkNk8VCqWphpSNTazBkUj+KfwO7SHA/ikxFCuVmZzULCnMtFLQrwdVx11xM5RknkTd+NPYxvbV40eGHYCkuxLXnAF01B8gbP7rXen9nFwDOdV/id7kx23PiHldPQcDEup50p8fqilEMLnaiFQbsSKJFMRzTY2uffg3p9zPm8vMpPeHolI8hb3w1x77wB/Y8+SJls2dQsWA25tzwv4GSY4+g8px51L/+Idt/8zjTfndrXMeKdX6da5T4RKiiACg7+Tga3lpK00efM+bSr8R1zHSSiutZr6KYjtKPYlg4d/1SudlZTAKLSZDnVnRj4bwTIvp1hdlM6YlH0/DmElqWreqvKLqVf530+XGs3UypGtOIh9BgdiKx2fQVcuuJUfR3PSXeokieotArb7Zgs9nwtCoPO8VHDbxPy2ApnDqRw+//qa5tJ992HQ1vL2XfotcZ/91vUHjowbqPE+38Bnw+nGq6bl9FUa7GKZqXrkQGAog4y9Gki1Rcz3r/E0uBw5I5kEyi26u4PSzqhWJXFUV3DM1dfOQUADp37um3LuDu0bFtKwYWp0iW68nh0JeO2jMzO4xFERh8MBug0KYqCk/yFIVeebMFh8OBz6n0WLEWFaZ5NPrInziWsd++AAIBtt79aFzfjXZ+O7fV4Hd1Yx87ClvZiF7r8iaOJbe6Em9LG+0btw1o3OkgFdez3jvOLuAdIcSfhRB3hb6SObh0IcyKoaXd9HK6FUXRlRNdUZjsShA84OlfPzHgCVUU6wc0rmRlPZWXl+vaLmrWk08LZg9OiRWmIEahV95soby8PKgoLEWJmYGdCg65+SrMBXk0vvcpzR/rjxtEO7+O1f3jExpCiGCabFOGpMkGpGRXi4tWV+SarKm4nvX+qvNQmgnZUPpSaK8xSRpXWnGrFoXmk7d1KfGFTlt0RWFWa/EHuvt76EKXtX2xfkAdtZKV9aTbojBFznoaTAmPUApSEKMYjhaFN6gohoZFAWArL2HiDy4HYMtdjyDDXHfhiHZ+wwWyQwm6nzIgTbbd7eOOxTu59j9fcvFzG7j2xc089uk+PqlpoyPkQSoV17OuGIWU8qpkDyST0J6YLWaB9PuxuFxIIXCao/sCTTlKcC4QJq851MrwNLfRtXs/+QfFp2eTURAQwOvVV0E+WtaTpjwGOuFOQ7Mokqko9MqbLXi93iFpUQBM+O7F7HnqPzjXbaH2v4sZrRbwi0a086spiqIIiqJsjlLOpHX52gEnnSSCHc1d3PXuLmrbPditJgIBya7Wbna1dvPfjY2YBBxSlsf00QWMMndRUuFPuKchlNT0nhxqmEyAH6tJBCcIuXNy8XqjP9GYcpRa9qHxCI2A293rc9sX6wegKBJfEBAS048iURZFKlxPw3EexQ4tRlE8dCwKAHNeLpP+57ts+NG9bL33T1SeMy9itpRGpPMbcHto37xdaed61JSw29jKSyicNon2Ddto/WId5aorKpW8vbWZPy7bi8cvOaTMzi9OO4iyPCtfNnSx5kA7a2rb+bKhi61NygvgkbXrmVqRx1GjCzlpwggmliU2E2pohPVTjFu9IVtMJrwO5Qfmzs2LmYmjtYEMF6Pwa13z1HapA4lT9HS3S9c8ishZT4lKjy2wqRPukhjMHm7zKGpra3tcT4VDy6IAqP7GWRQcejDd++vZ87cXY24f6fw6N25Hen3kHzI+arXc8pOVvhzNH6W2mqzHF+DBpXv47ZI9ePySs6aU8dBXJjOqMAeb2cSRowr49oxR/O7cybx4+RHce+bBXHzkSCYUmQlIyYb6Tp5bXcea2sSX9DAURRi0B2aLWeBTmxi57faYT7kmmxajcPdbF1DrR5WeoKTFtn4xAEWRpGB2fr6+m0e0rCfPIPpRhFKQAteTXnmzhTyLFen1YcqxxXwaz0SE2RwsXrjj90/jaXVG3T7S+Y0Vn9AoS0OcorbdzU2vbuXNLc3YzIIfnzyOH500DluEh0K71czMMUUsPK6ae08dzb8vO4I7T5/IBdMqktIl0lAUYVAf3LGaRLDssDvXHtuiyFUUhT+s60lZNuLYIxFWCx1f7gw+5eklWcFss1mf4omW9ZSIEh4ARSlwPemVN1uQnUrZ/KFoTWiUz5tF2cnH4nO0s1MtahiJSOfXuVafoig57ihMOTac67fiaUp+Ce/P9zq44aUtbG92MarQxu/Pm8wZk/VXhDWbzRTkWDhhfDHXHT+GcSMS38howHccIcTJQoiSRA4mU/AFetwovV1PsSwKNUYRJT3WWlxA0ZFTQEocqzbGNa5kzaNwOqM/oWlE60cxmA53oQSznpLoetIrb7bgqKsHwDLE4hOhKKU9rgchqHnyRbpqDkTcNtL5DVoUR0dXFGZ7DiXHHQlErt2WCPwBydMra/nft3fS7vZz/LgiHvnqFA4ui6+ibyqu58H8qj8Edgsh7knQWDKGoOvJJPCpiqJbh0WhJz3WZLMxYuY0IH73U7JcTxUVFbq2s0bpRzHQntl9KbD1uJ4CA0gh1oNeebOFIotyXVqHsEUBUDRtMqMvOhPp8bL1vj9F3C7c+fV1dtGxrQZhMVN42CExjxUsO54k95Oj28dtb+/gudV1mARcfewo7jh9YvBBKR5ScT0PWFFIKU3A4cCmxA0n/QSkDLqeLCaBNxij0BHMzo2SHqu6nkw5NkpmKlVn4w1oJ8v11NLSv09AOKLXekpMCQ+zSZBnNSGBziRZFXrlzRZaDyjBXUtxQZpHMngm3fJdTDk26l56F+f68J3zwp1f57otEAhQeNghuuI0ZXOVgHbTR58PaM5TLG5/Zyer9rdTnGvhvrMO4ZKjqjANsMdMKq7nQf2qpZT7pJQDqwWcofhUt4rFJBBC4OsVo9DpegoXo/BoisIa7LHtWLUprk5eySoKqPeHkOx+FBrJLgyYjB9+JuN3KinelsKhryjsY6oYd+WFAOyMUIY83PnVG8jWKJo2CWtpMd376+naFbs7Xzz4ApLNDZ0I4NELpnD06MG5BFNxPetthXqzEGK6+v54IcQeIcROIcQJyR1e6vH5exQFEIxReOx2uryBoCIJRzA9Noyi8ActihxyqyrIHVOFr72Tjq27dY9Ncz0luiigftdTtJnZiSnhAaFzKZKjKIab6ylPKP9PaxZYFADjr7kYYTFT/+oHYWMV4c5vvIpCmEzByXeJTpN1uHxIoCjXQkX+4Av6ZZLr6Uco9Z4A7gN+B9wDPJSMQaUTX5/idlp6LAXKjyyaVRFNUWjpsdo2WtOjlk9W6xqXlDJprqf6+npd20WdmZ1AiyLZPSn0ypsttB1QGvFkg0UBYK+uZNQFC5B+P7v/sqjf+nDnN15FAcnreqfVbSqxJ2a+cyquZ713nGIppUMIUQgcBfxRSvkEEH56YxwIIUYIIf4thPhSCLFZCHGCEKJUCLFYCLFN/VuibiuEEH8QQmwXQqwTQiS8l2OwCmofi8JcpCmKyE+5mhKIlh6rzd4uP/V4APY994ou09HjlwSkMq7BZhb1paBA3w0k+szsxFkU2qS7ZGU+6ZU3WzB71NplWWJRABz0faXh5v5/vIanpXeto77n19PiwFVzAJM9h/zJE3Qfo+wkRVG0fLySgC9xDy2tLmVfJXZrQvaXiutZ7696rxDiROASYImU0i+EKAIS8Uv+PfCWlHIqihLaDPwMeE9KOQl4T/0McBYwSX1dAzyWgOP3whdS5wkIzqPQUgujWhTR0mNVRaFlRo06bz62ilLaN23X1Wc7WeU74iH6zOzElPCA5LuehhvZFKPQKDz0YMpPPQG/q5s9T/0n6rYOdf5E0bTJmCz6n+Lzxo0i76Ax+No7ca79clDjDSXRFkUq0HvX+SlKl7vbgF+py84FBuW8U5XNycATAFJKj5SyDTgfpU836t+vqu/PB56RCp8BI4QQowYzhr5o/netF4UWzLYFFUVsi0J6vP0qXQaD2ersbVOOLRiU2/2X52OOK1mBbICODn0T/6LPzE5MUUBIfr0nvfIOBm+bU3e102TTrTYtyoasp1A0q2LPEy/gd/VUQ+h7fp0DcDtpJKPseFvQokiMokjF9azrVy2lfENKOVpKOUFKqc1AeQE4b5DHnwg0Ak8KIVYLIR4XQuQDlVLKWvXYtcBIdftqYG/I9/epyxJG0NeuuVnU2dO5JbEtCiFExHpPwXkUuT3Bq7Hf/irCZqVx8TI6d+4lGi6PGshOgkVRWVmpa7tIWU/+gOIWMwlIgEGR9DIeeuUdKF17ann/iHPZ+D+/Sepx9GJRrdGhVhAwFqWzj6HoqKl4mtvY/683gsv7nt+BxCc0ytU02eYliQto91gUiXE9Jft6Bp3VY4UQE6Os3jnI4x8D/EBKuVwI8Xt63ExhhxJmWb/H24aGBhYuXIjFYsHv93PhhRdy/fXXU1dXR35+PmazGafTSUVFBS0tLUgpqaiooL6+nmafahX4vXR3d+NtU2Y9Wgvt0OamvrWdpiZJZ2cnVVVV1NXVYbVaKS4upqmpCWG1gNtDzfYdVE86mLq6Omw2G361emyXx0N7fT3d3d1UVVVRtOBEHK99xLbHnqPkhm9SVlZGe3s7Ho8nuH+73U6jU1EUOSY4cOAAXq83uD6WTJoPs6Ojg8rKShobGxFCUFpaSmNjI263m9LS0ogyFRcX4/F4qG9VXBjdHi+1tbUUFhbS3NyMvbBYOZkCPB4PdXV15ObmYrfbaW1tjSiTzWbD4XBQXl6u9ExQZfJ2KE/ATe1d1NTUDEimoqIi/H5/WJl27tzJ2LFj8Xg8uFyu4HqbzRaUqaSkBJfLFTxP8ci0b+lnSK+Ppi/WUVNTk7DzFE2m0PPUV6auxmYA3GZBfci1N9jzlE6ZtPNUddWFOG+6l21/fIaRXz+D+sZGurq6qK6uDsrUslKpgGA/7GBqamrikqnkmEPBJGj9YgMdza00dzgHLVO9U7kX+LscuN0jBn3tbdu2jcrKyoScp4g3Xj2BVCFEQL0hazfq4JeklAP2hQghqoDPpJQT1M8noSiKQ4B5Uspa1bX0oZRyihDiz+r7f6rbb9G2C93vp59+KqdOHVhv4M0Nnfzwla1MrcjjobMP5p2xJyPMZupffp5n19TzraOruGJGZG/X+9POwdPUyinrXiVnZE+9lvemnoG3rZ1TN72JrbQ4uLx903aWnfptzHl25q1+KeJT34p9Tn7+1g6OHl3I/WfHnlkaD/v27WPMmNglzx3dPr7+7HoKc8y8ePmRweXObh8XhVk+UJbsbOXu93czZ0Ixvzwt2jPKwNAr70A58O+3WHfDXeRUlXPKmleSdhy9fDD7Ytw79jL7/Wd0zUoeSgR8PpbOvgRXzQGmP34PVeee0uv8dtc28uHR52MpKmD+l28NqA/2p2d/F8eqjRzz9wcYefrsQY/5lje2s/pAO/eccTDHjh18Ab9EXs+rVq1aOX/+/Jl9l+t1PZmklGb1rwnF3fMX4PLBDEpKWYcSKNeyp+ajzPR+BbhCXXYF8LL6/hXg22r20/GAo6+SGCyhTYt8wUB2AUWqmRhz0l2EFNm+6bEahYcdQumcGfi7XOz7x6sR95usOk8ApaWluraLlPWUyNRYCE2PTY7rSa+8AyU496apNSPiFLJr6BcFjITJYmHCtd8EYNcjzyGl7HV+HWuUwhHFR00dkJKAxHe9a0twMDvZ1zMMcGa2eoO/CWVOxWD5AfCcEGIdMB24F/g1cLoQYhtwuvoZ4A0UV9d24K/A9xNw/F4E+yqYRDA+YS0q0J2JEy5FVkoZMuGuv19ywjUXA1DzxL8jpuElaw4FQGNjo67tImU9BQPZA/wh9kWbmZ0sRaFX3oGiZcpJnz+oNNKJVq9sKBcFjMaYS87BWlqMY/UmWj9b0+v8xupopwetnEftf96Ju+JzOBKdHpvs6xkGV8JjCkov7UEhpVwjpZwppTxSSvlVKWWrlLJZSjlfSjlJ/duibiullNdLKQ+WUh4hpdTfcV0nPSU8TD0WRVGh7kwcLUVWhgSzpc8PgQDCbA6bnldx2onkHTSG7n11NL3/Wdj9JqsgIEBRkT7zV5ut7pf0KtiXqIKAGsFgtic5WU965R0o2nUD4GlKb10p6fcT6HKBEFgKBv1zzUjMebmMu+prgGJVhJ7fwQSyNUpmHcWI447E09zGjgefGtRY/QGJU72HFCfIokj29Qz6S3gsFUIsCXmtAJajzNDOKrTUT4s5xKIoLtD9lBu0KEIqyPbUeQo/XV+YTFRffDYA9W98FHYbl0+1KBJcvgPAr7PelBAibOZTokqMaxQl2aLQK+9ACX3qTEU/g2horXwthfkDdr0MBcZf9TVMuTYa3/2Eji1KEQkpZXD+w2AUhRCCQ391k1Li/PF/0bljz4D35ez2EZBK3xXtwWuwJPt6Bv0WxeMocx2016+BI6WUiXA9ZRTahDubSeBt0yyKAt0NdcxhKsgGU2PDuJ00Rp5xEgCNi5eFLRQYtChsibcoOjs7dW/bE6focT8lql+2ht1qwiQUmaPV1hoo8cg7EHwZpCg011c2xidCsZWXMOaScwHY99cXAOjavR9vWzu2ilJyR4+M9vWYFB05hTHfPBfp9fHl7X8Y8H4S7XaC5F/PoD+Y/XSf17+llNuSPbh04AuxKLQ6T9YRhfotijAVZLU5FZEsCoCCqROxjx+Np7mNtpX9GxoFFUUSLIpIzejDEW7SnTeB5TsATEKQH+xLEVkx72nt5ra3drClMb4fSjzyDgRvL9dTui0KzSrOzvhEKBO+dwmYTDjfWUZ3bWOvjnZigCW8AXa3uvj285tYedb5mAvyaHz3Exrf+3RA+9LmUIxI4KzsZF/PEEVRCCFuC3l/V6RX0keYYrwhZcaDT2NFhSE3Lj/+OCvIBtQ5FNqs7HAIIYJWRcNbS/utT2YwO1Iz+nCEy3xKZIlxDT2KefH2Fr7Y5+S+D2rw+PRnF8Uj70AItSjcjRliURRlt0UBkDdhDFXnzEN6fdT89V84Vg8+PuHs9nHH4p3Ud3hY6hAccvPVAHx5++/DluqJRWuCZ2VD8q9niG5RhCbmjo3yyiq8IRk8vpAYhdkkKLCZYzbUCaco/N3RYxQamqKofzucokheMNtq1W8Gh4tRJLIgoEZhMKAd+X/d1Kn8Xw843Ty/Tn8FzXjkHQiZaFFYirLfooCesh57//4SzcuUIhIDVRT+gOSe93dxwKlcZ81dXsZ/5+vkTRxL5/Y97Hnyxbj32ZbgWdmQ/OsZoigKKeV1Ie+vivRK+ghTTKjryRuS9QT6itWFtSjUJw9zbnRFUTLrSKwlRXTt2EPHtt291u1uVXLhKwsGX7++L8XFxbE3UulxPfU8wSd6HgX0tESNFhNq6ux5olu0tp79DnfEbUOJR96B0CtG0ZxmReHQUryz36IAKD76MIpnHYWvvZP2DYp3fKCK4s/L97P6QAcjci1YTYJOjx+3MDH1zhsB2P7bv+FujC+rrUW1KBLpekr29QzRXU8T9bySPsIUE9q4SFMU1hGaotDcIfFVkNWUhrBF1/wmi4WK+ScCvd1Pjm4fe9q6sZkFk8rtccmjh6amJt3bRrcoEul6iq2UG1VFcdjIfLx+ySOf7tVVsj0eeeNFStlr7kS6LQrvMLMoAIouPiP43j52FLayEXHv460tzby0sRGLSXD7aQdRlq/8dlu6fFScdiLlpxyPz9nBtvv/Etd+NYuiNC9xVkAyr2eNaL6C7cC2kL/bInzOKrwhjYt6nsaUOijxWBS90mODJcZj9+odeaYap3jn4+CyjfXKOKZW5Ce8FwXEaVGEiVH0/M8SN7aCGO1QpZRB19Mt88ZTYDOzYl87S3e1xdx3Mp/A/J1dEGJtudPtehpGMQqN0WfNpWCq8u0avnEAACAASURBVAw7EGtiY30Hf1imFOn8weyxHF5VQJl6Y2/u8iKEYOpdNyIsZvY992rE3t3hSEaMIq0WRWjZDuA7wCJgKpCr/v0HsDDpI0wxvpDGRV5n714UeibdadVhtQC28j52eqxG+SmzMOXYaFuxIWjWbqhTsnqmVSXnx+7x9G+0FImeLnch6bG+JFgUMVxPTrcfj1+SbzMzqiiHq48dDcBjn+2PGkOC+OSNF82aMNmVh4J0WxTBONswsii8Xi+Tf34dZnsuVV89La7vNnR4uOvdXfgCkvMPq+CsKUq9ttIQRQFQMGkC4xZeBFKy+RcP6e5b3Rp0PSXOokjm9ayh9xHwV8B3pJTb1J4R24BrgbuTN7T04FNvgEp6bF+LInYmjpbZFM71FCuYDWDJz6NszgyQksbFywDYUKeMY1pVcvoJuFwu3dv29M1O3oQ7CFHKEW76mjVRoboEzp5axpSKPJq7vDyzKnr5r3jk3f/8G3x+4Q24G5p1ba/Nys4bOxphNuNztIdtjZsqtMl/2daLIhoul4uRC2Zz2s73qDpnnu7vuX0B7nx3J60uH9NHF3Dt8T0dDMr6KAqAQ26+GlvZCFo/W0vdy+/pOkai6zxBfNfzQNH7yzYBE/osGw8kPgUnzXhCSngEJ9wV93U9RbYotA52gThmZvelQst+emspLq+fbU1dmAQcOjI5FkVc8yjCpccGkhDMjuF60uIT5aqiMAnBjbPHYhLw8sZGdjR3Rdy3XnnrXvuA9TfdQ8snq2j6UF8/gmBcq6Qo6Bv3NMd2hyWL4MNOFnW3i4V2fuOZOyGl5HdL97CtycWoQhv/e+pBvWZOa4qiJURRWIsLmXTrtQBs+dUj+NXii5EISElbt2pR5GbJPIo+PAi8L4S4VwhxnRDiXpQWpQ8mb2jpwRfSuKgnPVappaLLosjRgtlh0mOjzKMIZeQCpZRx85LP2bynBb+EiaX24FyORBPXPIqgRRGS9eRLXnpsJKWsZTxV5Pf8TyeV5/GVQysISPjjsn296lGFokfels/WsO76O0Hdh7fVEeMbCto1YykqwFZeAqQ3TuEbhhbFQOYVPL+ung92tGK3mrhzwUSK+tzIw1kUAGO+eS5FR0yme389Ox95NuoxtPIdhTnmhFrf6Z5HEURK+QBwFVCJ0tWuCrhaSpkZLbwSSDBG4fMScHsQNmsw7qDdvJzRYhRqwDrszOwY6bEauVUVFB99GIFuD9tfV7KfkuV2ArDpVGAQGqPoH8xORtZTRIuiQ/n/ahaFxpUzR1Fqt7CpoZOPdoZ/ko8lb8eWXay64hYCbg/WUsUq0BpYxUKLUViLC7BVKIoinXGKnhIew0dRxHM9A3y2x8GTXyjuylvmjWdCSf/MwqCi6OytKITZzNRf3QTArkeexbUv8k07GJ9IoDUB8cs7EHSrNSnlW1LKhVLKs6SUV0sp30rmwNKFdtMzdymuC2tRQdCE1RejCJce6+61Tg9V588HwP3WB0DyAtkAhYX6A53hZ2Ynrl+2RoFN/V9HiFE0dvW3KADybWbOPbQcIKL7KZq83bWNrLj0ZnyOdkaedTIH36S0RfG26lQUzp65N5pFkU5F0VPCY/goiniu5z2t3fz6g91I4IoZozhxfPhU2kgWBUDp8dOpOn8+gW4PW371SMRjtSWhzhPEJ+9A0f3LFkJMF0L8QAhxZzaX8AgGs11q1c2QGjl6CgOGTY/t1p8eqzHqq6eBEJSuXUuOq4tplcn7oTc36wvUQvieFB5fEiyK3Fiup/AWBfT8qB3d4b8bSV6vo50Vl95M9/56Rhx7BEc9emdPnEGnReELtSg0RRHnpKxE0lPCY/goCj3Xs8vr54V19fzk9W10eQOcfNAILp0eufe0No8inKIAmPKL6zHZc6h7+T1aPl0ddpvWJASyIb7f70DRW2b8GmAZcCpwC3AE8GOUlqVZhfakbO7osSg0dFkU4dJjdRQF7EtuVQX246Zj9vs4dsf6hE7Q6UtJSYnubcMWBQxoFkXiZ2ZHcj31xCj6/1+01EPtCa4v4eQNuD2svupWOjbvIP+QcRzz9G8w23Owlig56npjFD0WRQE5aY5RSCl7SnhkefXYUKJdz+1uH8+uquWyRRv56+cHaOv2cWRVAT8+eVzU4Hee1USOxUS3L0BXGCvXPqaKiddfBsDKy37K9t/+DV9nb4s2GamxEN/vd6DotSj+BzhTSnkB4FL/XgTEXxUrw9FiFCa1dK82Kxt0TriLmh4b3wXSNkcJah++PuH9mXoRV3pslKKAiQxm51pMWEwCj1/i7lPwT0oZEqPor3y18ghtESyKvvJKKVn3w7tp+WQVOSPLmPGPB4N9za0jlEQGb2t7v/2Eo8eiKMRWprSoTJTrafdfnmf5BdfHzK7RCLjcSK8PYbMGy98PB8Jdz60uL098vp/LF23kmVV1tLv9HF6Zz91nTOSBcw6JWUNNCBHV/QRw0PWXUXnOPPydXWx/4HGWHP8N9jz1HwJe5TpMRmospCY9Vu+IR0optZoSASGESUr5phDiuWQNLF0EYxSqogg12QtCXE8BKTGFeQIJmx4bxzyKUDZMOZKjLBbyN22m+0DDoGvqR6K7W9+NB8K7nhJdZhyUH2ZhjplWl48Ot5+ckPLq7W4/br8kz2oKmwlWrAYLI1kUfeVteu9T6l56F3N+HjP+8Vvyxo0KrrOWqIpCr+spJOtJe2hIlKLY+/eX6dy2G8eazZSeeHTM7bXyHaYs7WwXidDz29Dh4YV1Dby5pSn4QHNMdSGXTq/kiKqCuFJoy/KsHHC6ae7yMnZEbr/1ZnsORz9xLy2frmbLrx7FsWojm372f+z+67+YfOu1tBYcBCReUcTz+x0oeke8TwgxQUq5G9gKnC+EaALSN5MoSWjpsUKzKEJiFFazCbvVhMsbwOUNhL1JBYsCesIoijiyE6SUrHFKCqdMY/LGNdS+9G6wMmaiGVA/iiSXGQfF/dTq8tHu8QV9xBA+NTYULaskkkURKq+Ukm2/eRyAQ35yNUXTJvfa1hanogjOoyguDLYeTVRhQK2tql7Fo1k3thHJb5WZSVRVVdHh9vHXzw+weFtL0Etwwrhivjm9kqkDnI9UlqdcV5EsCo3SE47m+Nf/Qv3rH7L13j/RtWMPa75zG+MmHkz1KV9hxOmJLZGXSfMofgNoRVPuAp4F3gfuTMag0ol2UYmO/hYFxE6RjVY9Vm96LMB+pxtHt4+9xx4PwIH/vK37u/EyoH4Ugf5ZT4kMZkPkmFBTV+RANig9O2xmgdsXCPbxCCVU3sZ3Psa57ktsFaWMu+LCfttaigpAKHNqAr7YPbzDzaNIhEUR8PqCmVd6K5Zq8QlpT376ZCZRV1fHG1828+aWZvwBybyJI/jzhVO5c8HEASsJiJ751BchBFXnnsKcj57jsF//BFtFKUU7d3DxEw/h+emdtG/eMeBx9CUj5lEIxTZbAiwGkFK+CZQAJVLKx5I7vNQTnEjWoVkUfRVF9IB2uPRYvxrYNsfhelqv1ncqmns8lqIC2jdsC/YCTjS5uf3N6Eikome2RqS5FA0dvWdl90UIEXQ/hct80uSVUrLtAcWamHjj5Zjz+v8fhMkUjFP52mLHKbTZ/NbiAmxlPcFsvbWAIhGqbPQqnuCcjmE0hwKU83ugXfnNXTOrmp+fehAHlQ6+6nLfek96MFktjLvyQk7+7F+sOeM8PLYcXEs+Y9n8K9hw830JmbUfz+93oMT8ZUvlCl8PBEKWeaSUHZG/NXQJ3gC1pvTFvc32YEA7glvDFC5GoXNmttsXYPWBdp5eWcvza5WnhGljR1B17ilA8qwKu13/jyh8P4rkWBQFEdKRe+o8Rf5/BgPaYeIUmrwNby6hfcM2cqrKGXvZVyPuSwto60mR7amtVIQ5Lxdzfh7S4+3Vo2IghFoR7iadFoV6TM19Nlyw2+00qg8To4sSF8QPlvHojD+Hx5Rn56OTz+BvN99B9RUXIkyCff94laUnfZN9/3xtUA8S/9/ee4c5cpVp37+j1FJLrc5hunumZzwe54ADToAxDoBJBrMmL8lkE5ewwL4v3+6y7MILuyy7sESDgQUMCw4smGBsDDhh43HAYewZj6cndU5qdSvrfH9UnVJJKklVakk9M9R9XX21pErnqVNVTz3pfkKhkMam/Lt7uOfl7yY5MVP3virB7ivg/cAxNdc6AqBcTzJeTAioEK1lUejZJTmr9NgKrqd0Ns9X7t7P5d99iL+9aRffu3+Sg7E0bT4PZ23sZMPlzwVg4rqb1/xmaoWFBfuuEet+FM2xKIyiu1LXU5XUWIWuoJ4ia6HQFxYWkPl8wZp4z+vxhio/UAopstUVhcznC66nDi0+EejT6zDW6H4y12LYtij0seQCjQ2eHupYWFhgRn+ZGIg0LhXVieupFPFUjpwE0dPNyZ/5EM/83ffofdaZZOaXePgD/8w9L7uK+BN7HO9XSsn+X/2eP172Tv70yvczf8d29nz1Wsf7qQW7V9BtwC+FENcA+wDjKSGl/GbDR7WOUA/AfAWOnFocRNUaF1lVZu9dSPLPv32K3fNJBLC1N8SpGyKcuqGDk4fCRNp8yPNOo21DP4l9Eyzd/yhdp5+4NiFL0Nvba3vdQnpsczvcQeV2qAVCwMoWRWeosuupt7eXqZ/dRvyxJwkODzD62hdXHUchRba6osjGV0FKvJF2PD7t+G39PSTGD5KeXSC8dVPV7auhyPVkN0ahX8PhfvvzeySgt7eX6bjm0qlmdTreb42iu2ooLbYLb93EmT/6AhPX/Zod/99/sHD3A9xx0evZctVr2fq+N1Z9cQFNQczffh+7Pnc1C398ENAy9La86zVsetPLHY+vFuwqimcATwHPLvldAkeUolAWRX65kA9vRq0YhdcymF2eHiul5KbH5/jKXftJ5STD0TY+/pzNHNNfnsooPB66zzqFyRtvYXXPgYYriuXlZSIRe35s66ynxqfHgjlGUfywn6lSla3QVSVFNra4yO7PXQ3AUe97Q80aA3+3dg3UynzKmjKeFBpFDJiqw6LIGhZF83sqH0qYnl9iNZOnzecxrqFGoCdUUBRSSkeptQsWLVCFEAy//Hn0X3Quj3/qy+z/7o3s/vdvM3nDbzjh0x+i74Kzy/YjpWT+Dl1B3K0pCG9nhKOueh1jb345vkhzCittKQop5XOacvRDENm8BClJ6+RebQPFb2M1LQpDUZgsChWjMFF4/PbJBb5wu9ZF6+JtPbz73FHaq7DDqsrarB47aSScNS6q3I+iWVlPMZNS1jrbaed2oEr/8EKMovztb/YXfyD+xFMERwYZffWLao7DsChqKIpMrJwuo1GZT2YrIjXjTFHI9r+cYjuAyWWtrqA/7Hf0MK+F9oCXdr+H1UyeeDpnXJ92sFCF58nfFeWkz/4tI1dcyiMf+X/Ed+zmT6/6AEMvvZjj//F9xjNo7o7t7PrsN1i4+wF9uw42v/M1iIvP4qgT6+sLbhd/Wc5LG8jk8nQuzJFdWCLQ20VwpJj/pWbWU1s5hUfOojL70WntgX/FyQO89ewRasEX1iyNXLxyn4V6UV8/ChPXUxNIAaEQzDZnPa2kcySzeUJ+D+3+ysfrqpD1JHM5Fq+5AYCtf/MmW0SNAZs0HmaeJ2PbBvE9mQPYuZVVcqtJyywtM5Ti6hkZXtOxDzfIUCcQa6jbSaGn3c/qklZ050RR2KnK7j7rFM67+Rr2fPVadv3r1Uze8Btmb72bLe96DbO/u5cFnUPK19nBlne8irG3vAJfR5iU6VnTLDS+AXMdEEJ4hRD3CyF+pn/fIoT4oxBipxDih0KIgP57m/59l758c6PHkslLBg+MA9B52gllbyQqoLV/yXpyhF+7EGQ2h8xpDzgrrqeY/gA7qtdexpFhUcQbb1HU149CD/pLWYhReBpsUVi0QzXiE+3V3xYr0XhMXH8zq0/uI7RpmJFXvMDWOOzGKMzMsQqNsyiKt7fjylIWxXLuiKuLrYrdExpJXiMD2QpWDYzswC7Pk8fv46h3v45n/u779F98HtlYnJ2f/hoLd92Pr7ODoz/yVp5970/Y+oE3Gc+EQ6KOokV4H/CY6ftngM9LKbcBCxR6c18JLEgpj0ZrmvSZRg5CSkk2JxnavwfQFEUpThgM4xGwY2bFsjezEMJEDKhdTPmkohkvKIqlpLZtp01ueq9e5ZttgkVRV3qsrhyyeYkEvAK8jVYUqsud6TzP2ghkgynryRSjyGez7PpXLaS29W/ehMdv79wrGo9a6bGGRWFyPTWKGFDFKIRuAdlRPMqiaO9tPmncoYRYVnvBaIZFUW/mk1Pm2PZNGzj9u5/lad/4FN1nn8rRH34Lz773Jxz9N28qy8R0cv/WCzsFdx4hxIXqrb7REEKMAi8EvqF/F2gstT/WV/k2oJLcL9O/oy+/SDTQCZmXWnR+w37donhaud8vHPByXH+YvCz0si6FikUol5OyKLym9FhV2R21ab42M0bhqHFRCSmgQQjoa/w7R8SChNFu2mOnBY3HzM13sPrUfoJjwwz/1fNsj8NIj60ZoyhunQsQ6GsMMaByXUWO2azvr7YrS1kUbT3WPRaOVCyktWuyv0oMq17UqygKvSjsu6tUdffZN36Zoz/45jIFoXBINC6SUuaBG6WUzbJf/x2NnVY5vXuBRSmlusP3A8qJP4KWnou+fElfvyHI5CUil2Pg4F7AWlEAPG1Ym7D7D1pX6hZSZHVFYUEKqHzn0aC9rAwVo2iGRbG0ZI9CG8pJAZtBCKhgLrhT9SN2LQojPTZR2HZ1934A2s97mpG+ageG62nerkVh4XpaA9+TzOVIz2tz1HGcxhNkR/EoRZEgX2PNIwtTMS2YPVAlK65eGCmyK7XpXMyoFsxeK5zcv/XC7t3yeyHEOVLKuxt5cCHEi4BpKeV9QogL1M8Wq0obywxMT09z5ZVX4vP5yOVyXH755Vx11VVMTk4SDofxer3EYjH6+/uZn59HSkl/fz9P7Z+gd2YSfyaDf2SAfDjIvn37EELQ09PDzMwM0WiULWHtxtt+IMb4+Dh+v5/Ozk5mZ2fp7OwEnbI4GYszlVwhp7ue9k9O0iskiUSCmG6KLk1PIFZDhEIhFhYW6O3tZXl5mXQ6zdDQEJOTk4RCIZIeTczUUoyDBw+SyWSM5dVkmpqaMlJf4/E4g4ODzMzMFMkUCASYnZ1lZWXF2GepTOl0mkQigQxrb8jJdJZ4PM5TB6YB8JBnfHzc2D4YDNaUKRAIsLS0RF9fH0tLS5YyBTyQzsOu8X0EvYLJmHbeQjLN3NxcRZmi0ShBnyCZlSzEEyzPzzA/rikKEQ2zuLhoyKSOGQgE6OjoYG5uju7ubhKJBMlkkmhYCxon5xaIxWIVZVqamAJgOZ9hZWWFpaUlEknNAkxMzTI7O1vXPEVyAvJ5fF1R8rq1sjIxbXntmWVSrrJIfy/j4+NFMjV6nuq99qLRKLlczta1V2ue1PL5pP4SE5slHvc0VKb8ilafMbEYZ3x83LZMs3FNea3MT7EY6HEsU7V5ymQyTE9PN2SeKsGuohgHfiGEuJHygrtP2NyHFZ4BvEQI8QIgCETRLIwuIYRPtxpGgYP6+vuBjWhstj6gEyizwQcGBrjjjjvKDjY2NmZ8Vs0+2tsLdQtDw6MM7f85AL1PP4VgMMjGjRvLtj+3o5PA3bPsWUgRHdhmvCWEw5p7yBcKkUY7uZs2beIx3fW0edtWhMdDoD1COn8Qv1dwzFFjRlA2GtXeXM2Tpo65MKhTjCdSDA8Ply2vJJN5uSqsK5Xp4MGD9PX10dfXV7aNkklBuX7ywkMkEmFgyA8sEQr4jO3M21eTCaCrq6vsOEXbBw8wu5qhs2+IgUiAhYd3AbB5sIve3s6KMgF0hQ4yuZwmkfcyNjbGUlJXMoN9xnGtjqnGqsae1oPY+eUVotFoRZmWMtoDamDzJsLhMOFwGDk4yE6Ph3xshZ5oJ56A3/E8LT+qyRwc6KVz4zCTQH5xueo85bNZ8qtJEIJcwMfY6GiRTGXneY3zVO+1p2Dn2itdXjpPoN1vCyntkXDqtjFCfm9DZTp2bBj+vJN41mOsW0umTZs2sZzWLMCTjt5MsMRNW0umWvPU3t7OwMBA3TKZl4+Pj2MFu/6CEHADmoIYRXtYq7+6IaX8mJRyVEq5GXgVcKuU8rXAb9EaIwG8AbhR//xT/Tv68ltlAzktMnnJUJX4hELA5+FEvTXpgwfL4xTmojupKwnh9yE82ulWGU/RNp/tPO9mxigyGfv+1tL0WCM11tOcvIhSvqeZGhTjZhToxrVtMsr9E3VWlOTvtMcgm42VF9wJr9dogpSer48ATgXCA33dtPVrFl0tBtnssuai9HWEyeYqN9o60rCYzJLNa/VOtZoR1YN6YhQr6RyZvCTk95QpiUbAyf1bL+wW3L2p2QMpwd8C1woh/gmNZ+pq/ferge8KIXahWRKvauRBs3nJkJ4a22WR8WTGaSMR7j+4zP0Hl7lga3FWiRGjSKWrpsZ22oxPAHibGKOoqx+Fnh7bLPoOhdK6lWq9sktRSgyomDoHtm5xNAbFIJtZiJFdXDbiDqWo1J860N9DenaB9OwCwaF+R8eGQiA70N9tO93WTHfein4FhwoUGWAzMp6gmEHWbnX2Qh2BbCdoxfzaGrkQomKnDSnl7kYMREp5GxqnlNrnWRbrJIErGnE8K6TjCfqmDpL3eOg4qToH4mnDHcAED1gEtM3psSo+4bUMZNu/cAp1FI1XFJOTk0XmZzUoi0I1eGpmMBuKi+5W0jmDmiFSpYpdoZRqXCmKpXwGp70C/V1RMgsx0ouxiopCPZxLqenb+rqJU3/mk7Ie2vp7DIuiVgGfYd1EI47m93DHdBPIAM0I6tdePJ0jlsrZSm83aiiCzRlTK+bX7pNqF5rbyaw+lcun8fbdOiH+yBN48nkWR0bxhavnJh/d204k4GViOc3kcoqhjgJNgrIecqmUtUXhMDUWMLql5eKrjnlmasHKF1xxHKaCOyll07rbKZiL7gr04vaoGVRxk2FR6A/q6IjzNzA7NB6qu5254A7WXnRXsCh6bO8rEys03nIyv+uN+dUMT80nOGO0Pmp01Uu9WRYFaO6neDrH3ErGpqJoTq9shVbMr63XQCmlR0rp1f97gGHga8BfN3V0LcbqQzsAWByr7ZrwegSnbFBpssVxClVYJ9MZy9TYmF5s58Si8AT8eNoCyFyOfKKxJfter31d7xECn6nLXbMIARUMbq10zhSfsPdmZm6JmltNkksk8bQFKuajV4NRS1ElRbaSRWEQA9ZJ46G4ndr6e/D3dIIQpOeXbMVLfNGIo/ldT2Tzko/ctIuP/fJJds7WZzkb10iTLApw3sComamx4Oz+rRd13d1Sykng/cC/NHY464vEnx8HYHGTPTNOcz9R5n4yLIpk2pJifMmIUTh7wzDiFCuNdT/FYvb6QSuYe1I0O0YRMcUo7NCLm2F2Pak6hkBvF8vLtTvVlaIWg6zM541EA+UmVGiYRdHXg8fn05SWlFUpRbK6ReGPRhzP73rhph2z7F3U0kh3zSXq2oeyKAaaaVE4pBs3LIr25lgUrZjftbwGHguUc2Ifxkg/rFkUq9uOtrW+WVGYk68U+V8+ZVIUJubYguvJ2ZuAaobT6Myn/n5nAVZz5lOz+mUrmKnGnQSyoZhBVj2kA71djuWF2q6n7PIKSImvI4woecNbs6KYVTGK7qL/1eIUGZNFUY+8rcZKOsd3txc4iw4uJevaT8GiaK7rCezzPS022aJoxfzaDWb/geLCtnbgROAfmzGo9UA2vkJuZo6M3098cLD2BsDGrjbCAS8LiWxRYMurK4V82qQoguVZT05cT4DBNd/ogPb8/HxRXnUtaJlPWsqfyn5qtuspnsqhmM3t+p+7TVlPKpAd6Ot2LC+YGWStFUWljCegEIBeYzA7oO8n0NcNjz9FamaejgrbGBZFZ6QueVuNHzwwyVIyS5tXkMpJDsTqI4KYjhfiWM2C0xRZZVF0NSlG0Yr5tTvyb5R8XwEelFLubPB41g2+SJjorT/hCz+6l002yeKEEHQGvaykc8RTWUNRGOmxSVN6bJHryRkhoDHGJlkUTktRzK6ndDZf9FujodqhxlI5VjLaebP7EChyPWV1RdHbVVc72QKDrDVdglUNhcJamhfJfL5gDen7CdhQPGaLohntcxuJyeUU1z+s9Xl+69kjfPHO/RyMObcosnnJfCKDwL57sh44VxTNTY9txfzaraP4du21Dn9kPR4W+gc5vkYbQjO0B1maWCpnEFIZPSnS1umxhuvJqaJQPSkaHKNYm+upRRZFOmscy67ryUwMmEppD1V/va6nGgyymcXCg7kUgd76XU+ZxWVkNqcFpfVOfG02guNZU9bToe56+ua9B8nkJc/Z2s2FW7v54p37ORBLO87um1vJkJfQE/IZCRfNgFNFcUS7noQQttxKa6TwOKSgagOyaftZRR0WzXUKdRSpCllP9cUovE2qzp6amnKUh62sh2xeksk316Iwn1/V6c6u68nvLeS8x6d0901vt2N5wUaMokLGExTHKJw+/NIlbifz52qKxxhPNFKXvK3CY9Mr3LZ7kYBXcOXTh4m0+egM+lhKZplbzTiyDBS9TGeguW/YPXpQ2o6ikFI2PT22FfNbbeR26DkObZvWIZS/Peigx3ApxQQU0mPz6UyNymynMYrmVGfb7ZetoKqz0zlJOttci0JlPc2uZsjkJG1e4agPcmfQRzydY2VaDwj3dRNyKC+Y0mMrxigUxXh5/r8vHMIbCpJLJMnFV8uyoqpBuatUABvsdc0zV2b765C3FZBS8tW7DwBw+UkDRmvbkWgbS8ksB5ZSdSmKgUhzW7/2mILZeSnxVFH8qxnN6m7zeZpCKQLO7996UPFJtQ60HesOlerphI7FqrlOUXpsstiiSGbzpHISv1c45n0xbbKa3AAAIABJREFUgtlN4HtyAnNPCqVcmxej8BrHAs337OSNvCvk40AsRXK2EKOoB8r1VElRVLMoQLMCEnsPkppdcKQoDIvC1HzITnA8Y2qi5IwQu3X4w1OLPDq9QlfQxytPLSSQDHe28ej0CgdiKU4drhSuL8e0Tt/RG2puXUHA6yHa5iWWyrGUyNLdXvnF0k4L1MMBjl8DhRADQojLhRDHNWNA64msau+Zs39rqbfbmNn1ZE6PTas6Ck1R1EMIqGBUZ1vEKPLZLI9/8kvM6311nSAet27AVAnmnhTN6pet4PWIot7YduMTCqroLjNbiFE4lRdqu56qZT1B/SmyqdkCfUfpvlIzVVxPy4Xx1CNvs5HO5fnGvRop9OvP2EDYRMkyEtUsggMV2g1XgrIoOjzNJ0G0G6dodiAbnN+/9aDq3S2EGBVCXC+E2CGE+JYQ4kS0lqVfAR4SQjSUlG+9kdUfetGI/Tc+q77OnirpsfUQAipUY5Cdv/N+nvrS93jo3f9o9Oq2i0Gb6cAKZmLAjOF6al7w0NzE3mnao2pglFvQLIq2vm5H8qZzeZ6aT2gZZ1UYZKtlPYFZUTirzraMUaiueVWaIZkVl9P5bQVufGSGyeU0Y11BLj22uPfYsFIUMYeKQrcotgw1v/Wr3aI7u72y14JWzG+t18Avo7G0fgCN5+lXwFuklANo5Hwfb+7wWguDETVpvyq0Q39jNQezDZrxZNpoh6rSY+shBDT2WyVGoR4oyQNTzP3hT472OzMz42j9YtdTcy0KoCgm4ZTDR1kU6JZAoLfbkbzfuOcgb79uBw/PJPB3aUpAdbIzo5ZF0WbDCrBCWl/frCjMWU9WqZFSSsOi8Ecjjue32VhKZvn+A1qTp7eePVzWa32kU1MUBx0qCkUI6Es3/w3bvkXRfNdTK+a31t19HvBOKeUvgHcBg2h9KZBS3ggcmqkUdaIQo7D/dhyxtChUK9Ryrifloup0QAioUC1GYe51sP/7P3O0X6cusGLXU/MtiohJUTh2PYX8+NJpRDKJCPjxRtodybtL5xzaOZsw3E9pi1qKgkXRYNfTTHFVNoC3PYg33I5MZyyvhXwihcxk8QQDeNqcxXRagf/ePslKOsfpIx083YL8T7meDsZS5B3UCCj6jt4WxAPs8j01OzUWnN+/9aCWovCrXtlSylVguaRR0KF1Ba4RKkbREbZf5WgEs1PlwWxzeqy3xPXUUYdFYWQ9WcQoMvOFh9fUL39vVCLbQU9PT+2VTChyPTWZFBCgK7HCC374TTY/8YjjQqrOoI/QivYQb+vrNlpW2sWU/vCZjqerxikKFkUt15NTi6Lc9aTtr6toedFYVHyiQ1NaTue3mdi3mORnj80ggLedNWL5kAsHvHQFfaRz0uiRXgvJbJ5YKoffIxgb6q29wRpxKFkUrZjfWne3TwjxHCHEhUKICy2+Hx60lDahFEUqYT+ryGA3NSsKq/RY/bd6CQGheowiPacrCo8Hmc5w8Ce/sr1fp6arz+R6ajbNOMDWH3yP4/58H2fe/hvHMYqukI/2Ve3BqTKe7MqbzUvjQTAZT1dNka2d9bRWi6L4YVAt80m5xtRYDiXX0zfuOUhOwvOP7eWo3spU/sr9ZDdOoayJvrCfudnZtQ+0BuzyPRViFEe262ka+CZaZ7mrgbmS79NNHV2LoVxPkfbqvSjMKNBgm11PldNjl+skBITinhSlUNQSQy+8AID93/9f26X95v68dhCwyHpqlutp7vb76Lxd63/esbTomOytK+gjtKI/OHVFYVfe2ZW0wS81HU8XGGQtFEWlXhQKbX3O+Z6klCb6jlKLonJ1tkHf0VHef3k98dDEMnftXSLo8/CGMzZUXddp5lOhhiLQEnmdWxTNcz21Qt6qak7vZf0Xg6wemPU4qCM0t+pUVbceU8/sfEq70EstivqC2ZVJAVWMYuRVL2T+zu3Ed+xm6f5H6Tr9xJr7zTnMkip2Pcmi3xqJfDrDox/7nPE9ElugI+DsOF1BH+1x7cGpLAq78qq8fICp5Vqup8bHKLLLK+RTabztobJGWtWqsw36Dn0sTue3WfjNTm2sl5/Ub/j4K8FpQHvaaIHqb4m8TrOemul6aoW8zXMsH4ZQrqdcxj5zZZvPg98ryOQkKf2hWeB6MpECBpWiqI8QEKpXZqf1GEXbYC/DV1wKaFaFHaysOCvg8xe5nppnUez56rWs7BxHjg6Tagviz2TILjrrJRENFlxPfr1oza68iokU9ILKDmuLIp/NalaeEBWL6QrEgPbTY829sktRje/JCKzrFoXT+W0Wds5p1+2ZNrrXDddpUfSHAy2RtzvkR6AFq3P5yi+WzW5aBK2ZX1dRmKD87b26L9ouzD0TwJQeW9SPosT1tBZFsbxS5lZSrqdATxejr3kxABPX/8ZWkyOnzdnNWU+ZJpECJvZN8OS/fQuA4AffyXKn9mBMTjjzdno9gq6kdiPldRPdrrxT8eIXhkS7pgRKLYrssnaOfdEIwmN9HvzdURCCTI3OdGakDfqO8mBloIorK1NiUTid32Ygnc2zZz6BALZWiU0ojDispVA1FP2RQEvk9XkEnUEfeVnIbCpFIpMjlc0T8BYXjTYarZDXVRQmKIsiHrOmkq6EjkDB/QTmyuxMWYc7I5hdR3qsaodKPl/UDlVKaVgU/u5OIsdspuvMk8itrDL501tr7ndycrLmOmYYrqcmBrMf+8QXyCWSDL3kIk576fnI/j4Akgech8U6dEWR0eMHduWdLlEU8TbtAVeaHmu8wVdpserx+Qj0qHaq9q6v0j4UZlRzZWWXimMUTue3GdizkCQnYbSzzRbnkbIoJmKpqm/sCtNGjMLfMnlruZ/M1kQzU1hbIa+rKExQ7LFtNvtRKJRmPhmV2alCwZ2iiC40LaovYazgfiqYm7mVVWQ6gzcUxKtTpI++5iUA7P9B7ZoKv9+ZWWxYFE3qmT198x1M/+L3eMPtHPcP76Ur5Oe0p20GnFsUABHd9ZTs0CwKu/JOx4s76sX8mqIotSiM1NgK8QkFp3GKQgtUC9dTtRjFcqFpETif32bgCb0eZVufvdTz9oCXnpCPTN5eiuxMvOB6apW8tQLazW5YpNAKeV1FYYKqMu7scMbGWMogazQuSqWKGhethRBQwao6Oz2vPbj8PQWX2dBlF+INt7N4z0Ms79hddZ+dnc5cbeZ+FI3umZ1LpHjs7z4PwNEfvpLgBo1rPzis0RQkD0453qfKekq0a/NqV17lejplSNtuzhcEymMUBYuiOoGd0wZGqoq7rc/CouivkvVUkoHldH6bgZ0OFQVo5IAAB2o0MZJSGi1QByKBlslbW1E0P5ANrZlfV1GYoFxPq8vOmpWXMsgW+lFkitJjYya3U72mqFUtRUbPeAqYFIUv3M7wy58LwL7v3FB1n7MO8859JtdTowvudv/nd0nsPUjk+K2MXXmF8XtwwwAAyYPOc8bblrUH53JIe+DbkVdKabylnqQrimmv9uBqmUUxW831ZMOiiGrXitP5bQZ2zSlFYT/13G6K7HIqRzKbp93vIRzwtkzeWoqiFVXZ0Jr5dRWFCcr15NSiMFxPSWVR6HUU5vRYk6Ko1+0E1plPRnyip/jNYuMbXgbAgR/dVDWoXbdFUZQeu3aLIr5zD7u/+F0ATvz0h/CYXIDBEaUonFsUvpj2YF/UYwx25F1KZknlJJGAly09miUxgTavZRaFidK7Guz0kTDDir5Dwd/VgfB5ycbiRhfFwngOLYsik8uzZ16zCrb22rco7BbdmTOeoHXyGjQeFVxjrXI9uRZFi6FIAWXeGYN/pK04mC18XvB4IJ8nt6rdIJ62QN0tUM1QfE85U4zCnPFkRvTEbXSddQq5+CoHf/LrivtMp501slepsOlsvtCPYo2tJ2U+zyMf+gwynWH0NS+m++xTi5YbFsWEM4sil0ghkkmyXh+LHu1BYkdelZc/EAkwpDfCOZDzWTLIGgVuFZhjFdocup4q0XeAxu9jKJ4SupbMsmqDql0rTue30RhfSJLJS0aibUV04rUwEtUUdC2LQrmd+iPag7tV8tp3PTXXomiFvK6iMEG9HecctEKFQpX1su56EkIYKbKKxdPTFijUUNSR8aRgHaPQHhSlFgXApjdeDsC+a66rWKmdSNhny4VC1tNqptAGda1ZHfu+eyMLf3yQtoFejv3EVWXLg8NarCJ5cMpRM3lFxZ0IR1jUz78deVV8YiDip7vdh98jWEznDWVgZpC1bVHYaDhkRiX6DgUjoF1ioWSXiinPnc5vo1GIT9h3OwEMR7V7qJZFMR0vtihaJW+trKdWNS1qhbyuojBBxSiGBpw1K6/GIKty2j0Bv8n1tAaLwjJGoVsUFvUfQy+8gEBvF8uP7mLx3j9b7rPeOorVjPbgXWt8Ijkxw+Of/BIAx3/qA0YFtBm+SBhfNEI+ma7YZc4K6m1bUxTa+bcjr1IUg5EAHiEM6hCP7s4xxykKFkWjs57KKcbt7E9dG+paWe86ip1z2oPMSSAbCimyk8vpqimy5kA2tE7eWnxPrQpmH/F1FEKIjUKI3wohHhNCPCKEeJ/+e48Q4mYhxE79f7f+uxBC/IcQYpcQ4iEhxOmNHI+i8FicdxYcsmaQ1fv26vv0BtvWRN+h4AtXi1GUt/n0tAUYfa2WKrv329dZ7tNxHYXuZlrRLai1uJ2klDz60c+Si68y8PxnMfii51RcNzjsPE6hHqKr4Q7j/NuRdyZe4A4CGNTdGjk9fmWupTDe4G1mPdlRFNmVBLnVBJ62QJVqb02BlGY+GW1QO53VjTQLyqI42qGiCPm99Lb7yeZlWU2LGQWLQpujVsnbFfThEbCYzBovmWa0omkR/GXUUWSBD0opjwfOAa4SQpwAfBS4RUq5DbhF/w5wKbBN/3sbWmOlhkG5nkJtzojnIpYMssUXhyfgXxMhoIJhUcQtLAoL1xPA6OsuAyGY/N/fWqZTBgLO5DVcT7qiCPjqVxRTP/st07+6HV9HmBP+5UNVXViFzCf7tRRFFoV+49qR12xRaP/1jKewXp1tsmrsZj21OVAURsaTTo1uBSvFk89mtVa5QuDV+aGczq8Vbntygddd+zB/3OusGDWbl+ye1y0KGxXZpbBToW0Es/W5aoS8duD1CCP+YGVVtMr11Ap511VRSCknpJTb9c/LaG1WR4DLgG/rq30beKn++TLgO1LD3UCXEKI6DaUDZOrMelIP/riZQTZYPHki4F8TxbiCFYNstRgFQPumDfRffB4yneHAteUFeB0d9hvYQ8H1VLAo6ruMMosxo2bimP/zLqNmohIKmU9OFEUhRhFP58jk8rbkNd5S9YfPQIf2Pxkqp/EoUIzbraOw7kxXNO4q9B0KapmZP8qKTsTp/Jbi7r1LfPq2PUzHM/x8hzNre+9CkkxOsqEjYCR9OIEdckBF3zGgxyjWKq8T9LRrMpXGKVLZPKuZPH6PMFzTzUIr5G1+KyibEEJsBk4D/ggMSiknQFMmQogBfbURYJ9ps/36bxPmfU1PT3PllVfi8/nI5XJcfvnlXHXVVUxOThIOh/F6vcRiMfr7+5mf127a/v5+0lntwTc9OUF+dYnBwUFmZmaMZjczMzNEo1FyuRwrKysMDQ0xOTlJSmoXwuJqhsXFRdLpNOa8KdGmEZVNLWoPlFximfHxZWP7YDBIKBRiYWGB3t5elpeXSafTxvJQKEQgEGBpaYmcTheyODFNKpVicnKSxPQcAHOpVdpXV4tkmpqaIhKJ0PXyi5m5+Q7Gr7ke34ufjcfnNWRKJBL09fUVyeT3++ns7GR2dpbOzk7S6TSJRIKhoSFmprTTrUgQvUIyNTVFMpl0JNP+T36Z1PQc0TNPwnvR0xkfHzeWW81Tql17aEw9tpPwnCZzPB6vOk+zu8e1c67oO+aXWZgYZ9OmTUUyTU5OEggE6OjoYG5ujsll7cGUWZwm1TmML6kphlggSDsQn5xhYmKCdDpNWlcasVwGz+IiS0tL9PX1sbS0RCaTMfbf3t6OJ9hGPpEiPjvPUipRNk9KJu8Bzb2W7wixurpqee3FhHa9rk7MMD4+TmdnJ8tPabeHr6Od8fFxAoEA8XicQCBAd3c3iUTC0TyNr3r4zF1zBt36AweWefKpPYwOb6h5P01NTXG/Nk0Mt0MymbR1P5mvvZ6AduDH9s/wvK3RsnmKdnYxq1sUUX+e8fFxVlZWGB0dtXU/Wc1TLZnM89StvyQ+sW+SscgGQ6ZJXbF1Br3s3bu36v1UKpPTedq9ezdDQ0MNkakSDglFIYSIAD8B3i+ljFVxP1gtKHs1GxgY4I477ihbcWys0Lm1u1t7u2tvL/hNcwhAMjq8gUGdknrjxo2W2/f19Rm/aYG2GRJZSbSzE48QPNkRQb0DedsCRCIRMsIPZBnb0M+Y7q8171Pxypsnzby8q6uLiZ5dAAQRtLW1MTY2xm49eDl24vEE29uLZFLb97zkEvb/89dJ7D1I8In9DDz3GcbyWCxGNBotkkkhHC73j4+NjgCFN+pQwFfU4N2OTAt/fJDp//kVIuDnlM9/nMjoqOX2RfN0/DZmgLaVFL29Whcz9b/SPE3omVneHt1N4/GzefPmMg5/8zbethDx9Dx+r+CkbZvxCMEJm4fhkZ3E27Rz61lNsWGDZsw+pZ//wc2bCHV10dXVVXbu1P4Dfd0k90/iXU0yOmYtc29vL/t+fz8AnRtHaG9vt7z2ho/bxgSQX1w2lnv8moLxd0aN39T8QnHvglrzlA7382+37SSTl7zguF4enlxh72KS1VCfce0pWN1PY2Nj/OyAprietqmPYDBo635SCIfDbJWL8OcFlmVb2TEjkQhzKxlyUrPSO9pDdJiuZzv3kzqO1fJKMin09vbSv28vEEeEokXzNJlZASbobvfXvp9KZAJn82S+ntcq0/j4eNn4YP1jFAgh/GhK4ntSShVtnVIuJf2/8jXsBzaaNh8FDjZqLJds6+G523rIpqtTBpTC6xGEA14kBXeMxxTnUJ/XQgioUBqjKCUErATh8bBJL8Abv/pHRcvqTY81vtcRzFZ1HZvf9koi2zbb2qYQzHYeo1BNixYS2ZryKrfTQFjLeIJCrMKg8Vgsj1HUcj2B/ThFysh4Ki+2K92XOe5kMMeaAuD1pE/uW0zy8V8+yWomz7O3dPGe8zZyxogm33377Wed7ZzVjn10HfEJqF2dbSYDVGhlOnBvhaK7VjQsUjji02OFZjpcDTwmpfw306KfAm/QP78BuNH0++v17KdzgCXlomoE3vuMjXzo2WNkHdZRgDlFtphBFgqB7WZUZlsRAlbCyGtejDcUZO5397L82JPG78mkM8VYqhjqaVq0cO9DAAxc8gzb2xiKwgExoFIUqmnRUjJbU15zDYVCX9iPR5QrCiN47PEYweO8lFz74CSPTZf3CbCb+VSNENDYl0VdhsE7ZQqsO53f6Xiaj/5iF0vJLGeOdvCRC8bwegRnjOqK4oC9niC5vOTJOeccT2ZsMFJkrVlkZ0pqKMC5vGtBpaK7VqXGQmvkXW+L4hnAXwMXCiEe0P9eAHwauEQIsRO4RP8OcBOwG9gFfB14VzMGVU9ecqEnhaYovGaLItjWEEJAKLwp5vSgpRUhYCUEuqOMvPpFAOz5yg+M3+utozD265C+IxOLE9+xG+H3ET31ONvbmRWF3aI79RANDWguqsVEtqa85qpsBa9H0B8OkNR7Uqj0WNVNzh8NG8HjP+2P8c17J/i/v3rSyHxRqNbC1IxaxXYAgd5CZbbUu5wZGVimhAwn87uQyPDRX+xiZiXDCQNh/u9FW4wXgZOHIvg8gidmVo2XnmrYt5QklZMMRgJ1p4QHfR76wn5yUqunKMW0qso2KYpW1o1UKrprVVU2/AXUUUgpb5dSCinlKVLKp+l/N0kp56SUF0kpt+n/5/X1pZTyKinlVinlyVLKPzVjXPXkJStFEUsV8z2BpjQaQQgI4FV1FDp3k0EI2FteQ2GFzW97BQjBwet+TXJKy2Cptx9Fpe+1sHjfwyAl0VOONejX7cAXbsfX2aEV3dns6aAsio4B7aG6mMzWlHe6JDVWYTASIBnSzr9Kj1Vv8OZe2Y9Macojlsrx1T8eKNqHbYvCSI+trCg8fp/WECmfL4xnuTxV1+78rqRz/N0vn2T/UoqjeoJ88nlHFfWOCPm9nDgYRgIPHKxtVdRbkV2KQops+ZvzjIXrqZV1I5UsilalxsJfRh3FIYlgMOh4m9Kiu6IYRYOqsqG8MrsQn7DXYL198yiDl56PzGTZ+62fAM7lXatFsXjvwwB0P/0UR9sBRgqtnaK7XCJFbmVVs1x6tfOzlMjWlHeqpNhOYaCjYFEo11MhPlF4MO+YLqQu37JrgT+ZfPrV+kiYYVCMV7EooNxCKdCJFBSXnflNZfN84te72TWXYDga4J+ff7RxTZtx+oh999Ou2foqskthkANaxCmsXE/13L/1oqeG66nZxXbQGnldRWGBUMj5G1BHaU8Kc4wiGGDPgvY2tFYmSXOMQkpZkRCwGja/49UA7Pv29eRWk47l9QiBWTc4ZY5d1OMTXU8/2dF2YOpLYYMcUNVQBPq6jRt2MZmpKe90BUUxFAmY6ii0B6WqofDpPE+5vOTxGU2Jv+xETal94fZ9JBTdiU1iwGqEgGYYdOO6rBllUUQLwWw78/uVu/fz58k4ve1+/uXSo40HYCnO0Ptdbz+wXNP9Z1RkO2CMtYKyKA7Gyl1PpYSAUN/9Wy86gz68QotNprN54/eFFloUrZDXVRQWWFiwx8VjRqS0J0VbwaUiAn5+9JD2BvysLfYf6Fbw+H1aMZ/OTFur2M4KXU8/mc7TTiCzEOPAj26qS16fyd3khOspn82yuP1RALrPqsOiUOSAB2pbFOZAtlLQi4lsTXlLq7IVBiIBUsEQUgiyS8vks1nDslAZT/uWkqxm8gxE/Lz17BGO6gkxFU/z3e2T+vi1OEvclExQilwyRTYWR/i8+LtsMtKWWhSmDKxa8u6aXeWmHXN4BXzqeVvZ0FHZHXh0b4hom5epeLpqtXQuL9mlczwdvVbXU5UGRlZKvZ7ruV54hDCU6rwpHrXYwmB2K+R1FYUFVG6+E5S3Qy1cuLG8hz0LSfra/Vyyrfoboh0ovqfcympVQsBKEEIYVsWer/2Qnu7KmTWVYHY3ObEo4o89SW5llfbNIzXdKlZwZlGYFIXu8ltKZqvObzYvmV/NICi0QFUYjASQHg9ZPU6RXYqTKbEolNvpuP4wPo/gb561CY+A6x6e5omZVbrOOAlvpJ3440+xutc6YU+5pQJ93UaAvBJKM59KCQGh+vUspeS/7t6PBC47sZ+jaqSxeoQouJ/2V3Y/HYilSGbz9IX9aw7oFiyKYsWUzuVZSGTxCOgxHaOe+3ctsIpTtDKY3Qp5XUVhgeVle+l/ZnQEKrue9q5qJukVpww0pBOcOU5RjRCwGgZf+GyCo0Os7t7HgZ/f5ngM5hRZJzIt6PGJrjOdu53AWYzCcD31ml1P2arzO7uSJi8133NpkH5Qp/FImOIU2RKeJ5USe9yAts4x/e289MR+8hI+f/te8j4ffRecDcDMr2+3HrcN+g6FQjMk3fW0VN4bo5q8v39qkYcnV+gM+njdafayZ5T76b4Dlesp6ml9WgkbOtoQaFlPZvI9VbvQ2+7Ha7oe67l/14JSRZHO5llJ5/CKAg9cM9EKeV1FYYF6GoF0lDQvMqfHLuUEXUEflx7X15DxmXtS1CIErASPz8fmt74CgOnv3lhj7XL4i1xP9i0KIz5Rh9sJIDiiemfbsChMb+btfg9+jyCRyRNPVHaZVMp4Ao2dVACrQe2tO7MYMyjG/Z3aw3OHriiOHyg8IN9wxgYGIwGenEtw3Z+nGXjuM7Vj3WytKFIztTOeFFT3O8OiUOm6Joui0vWczOb5+j1aVtYbz9xgm4tJWRQPTsQtWVNBc2dBYxRFwOehP+InLzGoVaByLKnVjZpKO90pOvuukN8o2Gwm3MZF64R68pJLGWTN6bFZn5+XnzywpvoJM1SXO82icB6jUBh9zYvxdYRZ3f4YSw897mhbf5HryYFFcY+mKLrrCGSDU4ui4HoSQhhkjMGuygrbqthOwe/10NvuJ6EHtNPzS0VNi1bTOfYsJPGK4gBuyO/lfc/UCAW+s32C3FlngMfD/J33F/UVMcZdow+FGWVZTxbd9ipdzz9+aIrpeIajekI8/xj77ov+cIBNXUESmbxlUSEUKrLrYYy1glWFthHILnERtrr/RmlfilYGsuEvoI7iUMVa6ijihuupOD32Rcc3xpoAE4PsyqqRP+/UogDNhTXyGq0Ab+/V/+No22LXk723puTBaZIHpvB1hIkcu8XR8RQKRXcztRlYS6qyVUB79/7KSkYV21lZFKC9vRq1FCaLwtcZ4YnZVSRaX+i2kpeCM0ejXLi1m3RO8qVHl+g64yRkJsvs7+4pO4Zig7XqlV2K0hiFFYWH1fU8HU/zwwe18/Cuc0eKXDd2UI3OIy8lu9ZYkV0Koy1qzKworC2KVvffKC26K6TGtkZRuHUU64T60mOLXU/mGMWWoaijXsG1UByj0B+GDmMUCpve+HIQgokbflPWe7kazFaEXYtiQe+w13XmSTWDtJVgFN2l0mRqjNfseoLCjZsSlQOMpfTipRgsqaUwZxk9ZuF2MuMd54wQbfNy/8E4S6efph3v1+XklXZTY6E460lKWaDwMLVltbqer773IKmc5PwtXZyywTlNdTU6j4lYitVMnp52X8U0W6cYtqilKG2BqtDK9Fgoj1G0MpANbnrsuqGeRiAdJX2zF7KFU7tteG0psaUwYhTmYLaDrCczwltG6Tr/TPKpNPu//1Pb25ldT3YtisU/6YqijkI7M5RVkahBDliwKLSHqXI9reYqj7dSaqzCoLmWYj5WlPWkMp6O7bfuSNcV8vP2c0YAuDaiMXYs8eVlAAAby0lEQVTO/OZOg35DwaDvqMLzpGBYFHML5FaTyGwOTzBQZNGWXs8PT8b57ZMLBLyCt5w1XPMYVqhG5/GE4XZqjDUB1plPVjUU0LrGRQoFRaGdh1ZWZcNfQOOiQxVLS866eIHGSePzCFLZPOlsnu2zhZzvULt9mgo7UDGK1PScbULAauh42YUA7L3mevLZ2hw+AD6P8xjFWuMTCqrTXaoGOWAh60l3PemKYnrJ2q8OlQOkCgORAMn2gutJtUH1RSPsmFEWhbWiALj46B7O3hjlQNcAif5+MvOLRl2JMW4HFoW3PYQnpPW4SE1qAf7Slqzm6zkvJf91134ArjhlkKEqNRPVUETnMVFsVTQykK1QqKUwKQoTy68Z9dy/a0GpRTG/2rqqbGiNvK6isIDixncCIQqdrGZW0jwwU1AUHoetVWtBxShW92oM6/UEss0Ye+GFtB+1keSBKWYsXCFW8Duso8iuJFh+eCfC66Xz9BPqHisUOt0lDti0KAzXk3bjZrzWlAdSypqKYqgjYASztRiFZlEs+tpYSGSJtnkZjlaebyEE73/WJjqCPh47+kQApkvSZNM26TvU/tr07KiVJ/WmRdFiRWW+nn/9xDy75hL0tft5xSkDrAWV6il2Njg+Adp59whNkWdyWrq5QQhYMlf13L9rQUebF79HsJLOkcjkWm5RtEJeV1FYoF4NrdxPNzwyy4ooxCTM8YpGQMUoEuOaorBLCFgJseVlNr3pcgDGbQa1nabHLj3wGDKXo+PEo42CwXph9M6uYlHkkilycY3nSZ0vlcn0xEzccpvFZJZ0ThIJeCvGlDSLwhSj0Kk8diW183HcQLgm6WNvu593nzfK7mM1y2ril38oWq6C2XYsCigowpXde4FigkIoXM8r6RzfvFe7Zt5y1nAR4V89sKLzkFKaOJ4a5zsPeD30hwPkJUwsp1lN51hJ5wh4RVkP+lZbFMJcnb2abSnFOLgWxbohkylvlG4HKqD9i8dnyXkLF4mZzqMRUBZFQq/stUsIWAmZTIaRV74Qb3uI+Tu2F/WqqISAQ9eTEZ+os9DOjEIDo8rZS6WpsQBnbezE7xXsmMsYloMZtawJKGaQTU7NkUskwePh8RXtLfe4Km4nMy44qpvNzzmTVFuQ5M49xMe1eoZ8JqvVxng8tjPZlEJZ3W1tUajr+fv3T7KYzHLCQJjnbHVejV8KKzqPyeU08XSO7pDPcMk0CmZyQHPGU6lirvf+XQvM7qdWNi2C1sjrKgoL1JuXrCyKdE4SChfcG96Gu54KMQqoP+NJYWhoCH80wsgrLgUwWGWrwWkwe1HFJ85qpKKoXHRXGsgGCAe8nDfWiQRu2VXeD6JWIBugzech0KUpZsP11xlhx4zmbjm+3561JITg3edv5uCxmvvp1u/8Sh+3Hlfp6UR47b3xqzTalSc1i6I0RjE0NMSBpSTXP6Kdr3edO7omqnsFKzoPMxFgI45hRoFuPGWkMZfWUEDr6yigkCI7v5oxCu7cOoojHPXmJXeYTOBzji4UMDXa9eTtKH4YrTVGoeTd9KaXA3Dwf35pUEFUgt8BKaDM57UeFLTQojBSY4uVqOLa+s3O+bI6DKuGRVaIDmoP5rzeZdAbjRjulmNtKgrQYibHv+wCACZ/dTv7FpOFcTvgwSq4nqwtij37J/jkLXvI5iXP3dbDMQ7GWAtm9xPAzrnGu50UlEVx0GRRlKbGQuvrKKBgUUzF0yyncnjE2lsK2IVbR7FOsGqAbgeKAsEj4PxjB43fGx/MLh6fE0JAKyh5I8duofdZZ5JLJDlw7c+rbuMkmL2yc5zM4jLB4QFCo2t/+ynEKCoX3Zl5nsw4YyRKtM3DvqUUT8yuFi2brlKVbUZPfxfS9LacD4fJ5CWbuoK2aTAUznvlxUghGHlqJ//+qx0kpmt3tiuFUiqpifKsp1xe8o1HEuyeTzASbeNtZ484Gl8tKIvigYllsnlZsCgaGMhWMDcwquYmrPf+XQt62rV5V61fu4K+ltB3QGvkdRWFBbw2Tf5SqDz9czZ10t9dmDxPg/OcVYxCwSkhYCnM8m668q8AGP/mj8vy+4uO6YAUUPXHboQ1AeALh/B3dSDTmYoNgEqrshW8HsGzNmrFaL/ZWex+suN6AhiIBkkFC2/MiTbtc6VCu2oI9HTSeebJeHM5knfdxx3b92i/26jKViitt1AEhVJKvnTnfh6cThJt8/JPz9va8LdcM53Ho1MrRmrsMc1QFJ2FvhSV6Dug/vt3LVAWhaJWb1VqLLRGXldRWCAWq8yKWQ3PP6aXl5zQxzvOGSluXNSk9FiFeug7zDDLO3DJMwht3EBi/CDTN1dOlfUVVWZXf3NaVBXZDYhPKLSZrAorlKbGmnG6nk342ycXjFRLsBfMBi1VUwW0AWJ+LR5VqdCuFoae+wwAjtrxMNsf1txHbTYIARVKyQNVv+wf/3man+2YxeeBf7jkKONB22goOo9fPj5LLJWjM+izfICvFUMdbUaKrKrQtqqgr/f+XQuUolDjalV8Alojr6soLNDf31/Xdr1hP+8+b6N2QZsynTzBBiuKjhLX0xrTY83yCq+XsbdcAcCer/6w4jZ2LYp8OsPs7+8FoLtBFgVAqEacwvD1W5yb07cMsaU7SCyV414TV5FdRTEQCRhU4wDzHm39eiwKwGCTPXb3o4RimoLz2qjKVihVhv7OCH94apGv36MF299/7gZOHIpYbdoQKDqPW5/UzvnRvaGGB7JBK/Ic6gggwXAblhbbQf3371qgFIVyhLZSUbRCXldRWGB+vjwjxik8AZNF0WDXk7fU9bTG9NhSeUdf82K8kXYW7rqf2J+tWWXtZj3t+++fkpqYIXLMFqInH7OmcZrRViPzqZLrCbSOYBebgtoAiUyO5VQOv1fUJHMz03gALHjbCPo8bO6uL4AbPmYzobFhfLFljnlCC/rfE7d/a5bGMybyPj5z2x4A3vz0DZzQYa/avl4oOg/FON7IQrtSDOtxCkVvXkrfAY25f52iNBW4la6nVsjrKgoL1GIltQOz66nR6bEenw+PibJjremxpfL6OsKMvubFQGWrwg4pYHYlwZOf/xYA2z76NtvpnnZQ06KoEMwGTd4Lt/bgEXD33hixZLZALx4O1AxCmmspAFLBEMf0tTtmYFUQQhhWRcespvh+tyC5/6C9hjT+7mjRuf3WjhjpnOTSY3t55SmDDbmeq0HReSistfVpNaiANmhZhlZFg82W1wrhgJc20wtTKy2KVsjrKgoLNMKU8/h8xs3b6PRYoKi6uV5CQAUreceuvAI8HiZu/A3Jqdmy5cqK8AioZFDs/eaPSc/MEz31OAYuPX9NYyxFW43q7FLmWDP6+/vpDfs5faSDbF7yu90Ltt1OAO0BL7mOQmZRKhSq2+2koBSFwmqkg8/eNl5GuGcF4fEUWU7zws8ZIx285xkbEUK0xDWhsp+gNRYFWAeyYX1cT0IIo5YCWldsB67rad0wNVW7KY4dqCB2oyuzoRCnWCshIFjL2z42zOCl5yMzWcsCPBWj8Hs9lv7ozNIyT33pvwE45uPvaLjPOqQ63VXge6rmelLyXny07n7aNW+qobB3g/u7TIoi2G67IrsSus8+tSj2NDw2yOxqhi/csc/WG6PPJGf/UDf/56ItBnFjo67nalD1FB1tXoZsKNt6YQ7IW9VQQGvktYKZUr1VvSigNfK6isICkUhjAn+B/m5EwF9WANUIqMyntRbbQWV5N7/tlQDs+84N5Erahyp3U6X4xJ6v/IDM4jI9551O7/lPX/MYS9GmOt1ZWBT5VLrA8xQtl03Je97mLkJ+D49Nr7Jd7/9cKzXWOL7JiksGQ2tWFJ6An77nnGN8f/+LTibk9/CHpxa5eWd1H/S+xSS7ZeEB+tEXnlDEVdWo67katvWGeOc5I3z42WNNCWQrqAZGULlnSCvktYI5TtFK11Mr5HUVRRNx+jWf4ek//Pc1k+BZwasX3a0146kaus46heipx5GZX+Lgj39RtEwFs61SY1Mz80ZsY9vH3t6UB0dR0V0+X7TMiufJch8+D+dv0c7f7Xs0YjU7rieA9r6Comjv6WgIr9GAnibr744y3BvmqnNHAfjSXfuL+jCY8bvdC7z7xseZb9MVlcfD0ODaXx6cQgjBy04a4JxNzT32UEfAcHU2IwV3LShWFIfW2NYKV1FYIB63Zhd1io7jt9Jz7mkN2VcpDItijRlPUFleIQSb3/4qAPZ87UdFLhCfR+BLp7G6HXb/53fJrSbov+QZa+49UQlFRXclne5Ss5UD2VAsr3I/Kdi1KKIDhX2PDtvvN10N/RefR3Bk0LDALtnWw7O3dJHI5Pn0b/cYmT4AmVye/7prP5+6dQ+JTJ6eEa04xNcRLuse2Kjr+VCA1yOMHhqVXE/rJa9SFIJC8W0r0Ap5WydNAyGEeD7wBcALfENK+elG7n9wcLD2SusM5c9ea8YTVJd36MUX8vgnv8TKzj3s+a/vI/M5Yg89wcr9j/He/ROsdnWzY+elDL/8uXScuI3kgSn2XnMdoGU6NRPB4UEyi8skD04XpYiWNiwqhVnekzdEGIj4bfM8KfQMdKMIQLZsakw/AH9XlGff82MjCUIIwXufuZFHplfYMbPK9++f5PVnbGA6nuafbnmKHTOr+DyCt589wkmJzey8qbgFqsLhcD07wfGDYQ7EUhX5pNZLXqUookFf3Rlw9aAV8h52FoUQwgt8CbgUOAF4tRBibZ1wSjAzU5mV9FBBI2MU1eT1+H2MvVkjC3z8k1/iiU99hcn/vZX8/gnyQtC+uMCeL3+fOy9+I7ef/1oefMcnkOkMQy+9mOiJ29Y8tmoIqjhFSYpspkpVNhTL6xGCi3SrQgB9Nt0ZfRsKimnblsZlnZSmEHe0+fjbZ48hgO8/MMkPH5ziXdfvYMfMKv1hP//6om1cdmK/oSitYjKHw/XsBO9/5ka+88oTGKtQt7Je8ipF0cr4BLRG3sPRojgL2CWl3A0ghLgWuAx4tOpWDtDMYFyjoB4IldwrTlBL3o2vfxmzv72HfDZL9ORjiJ58LAujY3z44QTnJqe5YvIxJn56Kys792j783rZ9uG3rHlctRDUM5/Gv/4/zP3uXuP35R1aP41KFkWpvJds6+GHD04xEm2z3dZ1w2gfe4Gcx8O2Uft0G/Xg1OEOXnHKAD98aJqr9cZDZ4528NELNhvcTW0DmvvLnI2lcDhcz04Q8HqqtnBdL3mPHwhz7lgn5zY5TlOKVsh7OCqKEWCf6ft+4GzzCtPT01x55ZX4fD5yuRyXX345V111FZOTk4TDYbxeL7FYjP7+fubnNbrp/v5+pqamiEQieDwexsfHGRwcZGZmRutg1dPDzMwM0WiUXC7HysoKQ0NDTE5O4vf76ezsZHZ2ls7OTtLpNIlEwlgeCATo6Ohgbm6O7u5uEokEyWTSWB4MBgmFQiwsLNDb28vy8jLpdNpYHgqFCAQCLC0t0dfXx9LSEuL80xiamkWeezKzs7M1ZQLNl2klU1tbG7Ozs1Vl2vyVTxTJ1OP380p/lk3Bbka3PIued7+aud/fS/aPDyOPGiYeaSMXizmSKZPJGMvtzBMbNJfP/J3bmb9ze9mFku/tZHx8vEymZDLJ4uJi0Tx95Mwo3e0B4vG4rXnq6ekhe/H50BFG5LOMj+9vjEwV5umC/gz39QR4cj7NS7aGePO5o0xPHSCuyxTf2MfAKy4leuFZRTIHAgGCwSDj4+MNu/YaJVOz7iePx0Os2ddeBZnefIyXaNRT835q5DMimUwyPT3dEJkqQaxHFeNaIIS4AnielPIt+ve/Bs6SUr5HrXPXXXfJ4447ru5jjI+PMzY2tuaxHi44XOXNJVJM/u+tZOOrZct8kXaGXvQcvO3l/bEPV3nTuTzLqZzjDKvDVd564cpbP7Zv337fRRdddGbp74ejRbEf2Gj6PgocbOQBotG1ZxIdTjhc5fWG2oyufE5wuMob8HrobXceVjxc5a0XrryNx2EXzAbuBbYJIbYIIQLAq4CfNvIAuSp9GI5EuPIe2XDlPbLRCnkPO0UhpcwC7wZ+BTwG/EhK+Ugjj7GystLI3R3ycOU9suHKe2SjFfIejq4npJQ3ATc1a//r0Zx9PeHKe2TDlffIRivkPewsilZgPZqzrydceY9suPIe2WiFvK6isMANN9yw3kNoKVx5j2y48h7ZaIW8rqKwwHXXXbfeQ2gpXHmPbLjyHtlohbyuorBANtvc1pGHGlx5j2y48h7ZaIW8h13BnR3ccsstM8B4vdvPz8/39fT0lLd1O0Lhyntkw5X3yEaD5R276KKLysjLjkhF4cKFCxcuGgfX9eTChQsXLqrCVRQuXLhw4aIqXEVRAiGEVwhxvxDiZ+s9lmZDCNElhPixEGKHEOIxIcS56z2mZkII8QEhxCNCiIeFED8QQpQzBh7GEEJ8UwgxLYR42PRbjxDiZiHETv3/2nnpDxFUkPez+vX8kBDieiFE83oFtxhW8pqWfUgIIYUQjemiVQJXUZTjfWjUIH8J+ALwSynlccCpHMFyCyFGgPcCZ0opT0Lrjviq9R1Vw3EN8PyS3z4K3CKl3Abcon8/UnAN5fLeDJwkpTwFeAL4WKsH1URcQ7m8CCE2ApcAe5t1YFdRmCCEGAVeCHxjvcfSbAghosD5wNUAUsq0lHKx+laHPXxASAjhA9ppMOvwekNK+XtgvuTny4Bv65+/Dby0pYNqIqzklVL+WueDA7gbjV36iECF+QX4PPARoGmZSa6iKMa/o53w/HoPpAU4CpgBvqW72r4hhAiv96CaBSnlAeBzaG9dE8CSlPLX6zuqlmBQSjkBoP8fWOfxtBJvBn6x3oNoJoQQLwEOSCkfbOZxXEWhQwjxImBaSnnfeo+lRfABpwNfllKeBqxwZLkliqD75i8DtgDDQFgI8br1HZWLZkEI8XdAFvjeeo+lWRBCtAN/B3yi2cdyFUUBzwBeIoTYA1wLXCiE+O/1HVJTsR/YL6X8o/79x2iK40jFxcBTUsoZKWUGuA44b53H1ApMCSE2AOj/p9d5PE2HEOINwIuA18oju1BsK9qLz4P6c2sU2C6EaDidrKsodEgpPyalHJVSbkYLct4qpTxi3zillJPAPiHEsfpPFwGPruOQmo29wDlCiHahdaO/iCM4eG/CT4E36J/fANy4jmNpOoQQzwf+FniJlLK8R+4RBCnln6WUA1LKzfpzaz9wun5vNxSuovjLxnuA7wkhHgKeBvzzOo+nadAtpx8D24E/o137X1vXQTUYQogfAHcBxwoh9gshrgQ+DVwihNiJlhnz6fUcYyNRQd4vAh3AzUKIB4QQX1nXQTYQFeRtzbGPbMvMhQsXLlysFa5F4cKFCxcuqsJVFC5cuHDhoipcReHChQsXLqrCVRQuXLhw4aIqXEXhwoULFy6qwlUULtYFQoi/r1bQqLO8XtDCIanj3iaEeIv++bVCiL8Emg8XLqrCVRQumgIhxBuFEH8WQqwKISaFEF92QvkspTxRSnnbGsdwjRDin+rdXkr5PSnlc5t9nMMNQojNOqW1b73H4qI1cBWFi4ZDCPFB4DPAh4FO4BxgDK0IKrCeY3Ox/nAVzOEHV1G4aCh0+vJ/AN4jpfyllDIjpdwDvAJNWZhpUYJCiB8KIZaFENuFEKea9rNHCHGx/tkjhPioEOJJIcScEOJHQoge07rPFELcKYRYFELs062ZtwGvBT4ihIgLIf63wngv0RvdLAkhvggI07I3CiFu1z8LIcTn9cYxS3pjnJMqHcc03mUhxKNCiJeV7lcI8TkhxIIQ4ikhxKWm5T1CiG8JIQ7qy28wLXuRXnG8qMt8is15uUCv5v2ILsOEEOKlQogXCCGeEELMCyE+blq/2jn/vf5/UZf5XH2bNwutAdaCEOJXQogx0/6kEOIqvUJ8Z6XzaUcWF+sAKaX75/417A+tsUoW8Fks+zbwA/3z3wMZ4K8AP/Ah4CnAry/fA1ysf34/hd4CbcBXTfvZBCwDr9b30ws8TV92DfBPVcbaB8RMY/iAPva36MvfCNyuf34ecB/QhaZMjgc2VDoOcAUaS60HeCUaO+8G034zwFvRGii9E603hmJK+DnwQ6BbH9ez9d9PRyP1O1vf7g36eWqzMS8X6LJ9Qt/nW9Fo5r+PRnlxIpAEjrJxzjej9T7wmfb/UmCXfl58wP8B7jQtl2hNhXqAULXz6f4den/rPgD378j6Q7MYJiss+zRws/7574G7Tcs8aH0inqV/30NBUTwGXGRad4P+oPWhdTC7vsLxaimK15eMQaARq1kpigvROqadA3icHEdf5wHgMtN+d5mWtesP0iFdtjzQbbGPLwOfLPntcaVIahz/AiABePXvHfoxzzatcx/wUhvn3EpR/AK4smQ+V4Ex/bsELjQtr3g+3b9D7891PbloNGaBvgp+6A36coV96oOUMo/2kB622G4MuF53tyyiPcRywCCwEXiyzrEOl4xBmr+bIaW8FY1w7kto1N1f091slhBCvN7kIloETkKzYBQMhk9ZYDmN6PLMSykXLHY7BnxQ7VPf70asz5kV5qSUOf1zQv8/ZVqe0MegjlXpnFthDPiCaf15NMU7YlrHfK4dnU8X6wtXUbhoNO4CUsDl5h+F1j3vUrS+zQobTcs9aG4Oq/ak+4BLpZRdpr+g1LrW7UPj5bdCLcbLiZIxCPP3sp1J+R9SyjPQ3DTHoAXry46j++a/Drwb6JVSdgEPY4p/VME+oEdYZ4jtAz5Vch7apZQ/sLFfp6h2zq3O6z7g7SXrh6SUd5rWKdquyvl0cYjBVRQuGgop5RJaMPs/hRDPF0L4hRCbgf9Bsxi+a1r9DCHE5br18X40BXO3xW6/AnxKBUeFEP1CiMv0Zd8DLhZCvEII4RNC9AohnqYvm0Jr+VoJPwdONI3hvWjunzIIIZ4uhDhbCOFHizck0d6wrY4TRnsozujbvgnNoqgJqbUr/QXwX0KIbv38na8v/jrwDn0cQggRFkK8UAjRoR/nGiHENXaOYwPVzvkMmnvsqJL1PyaEOFFfv1MIcUWlndc4ny4OMbiKwkXDIaX8f8DH0XpUx4A/or1xXiSlTJlWvREt0LsA/DVwudS6z5XiC2gNeH4thFhGUyZn68faC7wA+CCau+MBQGVPXQ2coLtDbijdqZRyFi3o/GlgDtgG3FFBrCjag3oBGNfX/5zVcaSUjwL/imZdTQEnV9mvFf4aLR6wAy14/X59vH9CC0J/UR/HLrR4h8JGh8ephmrnfBX4FHCHLvM5Usrr0VKirxVCxNAsqEutdw1UP58uDjG4/ShcHJIQQuwFXiel/H3NlV0gtPqUB4FTKihbFy7qhlv44uKQgxCiH+hHy3xyYQNSyjRaiqkLFw2H63pycUhBCPF0YCfwn7pbyYULF+sM1/XkwoULFy6qwrUoXLhw4cJFVbiKwoULFy5cVIWrKFy4cOHCRVW4isKFCxcuXFSFqyhcuHDhwkVVuIrChQsXLlxUxf8PqUemz8MLeYwAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots()\n",
"\n",
"ax.plot(distances/1e3, curves['600/4'], label='600/4')\n",
"ax.plot(distances/1e3, curves['400/2.8'], label='400/2.8')\n",
"ax.legend()\n",
"ax.set(xlabel='Object distance, meters', ylabel='Blur radius, mm', title='15:1 magnification')\n",
"plt.savefig('blur-vs-object-distance.png', dpi=200, bbox_inches='tight')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment