{
"cells": [
{
"cell_type": "markdown",
"id": "e1493f15",
"metadata": {},
"source": [
"# Bessel Functions"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "70de4483-0fb1-4fb2-b5d6-0ef7f983ce43",
"metadata": {},
"outputs": [],
"source": [
"using SymPy"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "6e80a623-96ce-4fc2-b28b-da4b0c017a5a",
"metadata": {},
"outputs": [],
"source": [
"using Plots"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "48720af0-919b-4984-986f-aa66ea09b0f1",
"metadata": {},
"outputs": [],
"source": [
"using SpecialFunctions"
]
},
{
"cell_type": "markdown",
"id": "fa5425ce",
"metadata": {},
"source": [
"## 1. Symbolic Manipulations"
]
},
{
"cell_type": "markdown",
"id": "fbac2af2",
"metadata": {},
"source": [
"We first verify that the Bessel functions satisfy the ordinary differential equation\n",
"\n",
"$$\n",
" z^2 \\frac{\\partial^2 R}{\\partial z^2}\n",
" + z \\frac{\\partial R}{\\partial z} + z^2 R = 0.\n",
"$$"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "c318d6e2",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$z$"
],
"text/plain": [
"z"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"z = Sym(\"z\")"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "fe30749b",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$J_{0}\\left(z\\right)$"
],
"text/plain": [
"besselj(0, z)"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"J0 = sympy.besselj(0, z)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "17eb2747",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$- J_{1}\\left(z\\right)$"
],
"text/plain": [
"-besselj(1, z)"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d1J0 = J0.diff(z)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "8843ed34",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$- \\frac{J_{0}\\left(z\\right)}{2} + \\frac{J_{2}\\left(z\\right)}{2}$"
],
"text/plain": [
" besselj(0, z) besselj(2, z)\n",
"- ------------- + -------------\n",
" 2 2 "
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d2J0 = d1J0.diff(z) "
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "4055407f",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$z^{2} \\left(- \\frac{J_{0}\\left(z\\right)}{2} + \\frac{J_{2}\\left(z\\right)}{2}\\right) + z^{2} J_{0}\\left(z\\right) - z J_{1}\\left(z\\right)$"
],
"text/plain": [
" 2 / besselj(0, z) besselj(2, z)\\ 2 \n",
"z *|- ------------- + -------------| + z *besselj(0, z) - z*besselj(1, z)\n",
" \\ 2 2 / "
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"odeJ0 = z^2*d2J0 + z*d1J0 + z^2*J0"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "663f793b",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$0$"
],
"text/plain": [
"0"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"sympy.simplify(odeJ0)"
]
},
{
"cell_type": "markdown",
"id": "8a4729f1",
"metadata": {},
"source": [
"The above calculation verifies that the Bessel function of the first kind satisfies the ordinary differential equation. Let us do the same calculation for the second Bessel function."
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "22496bcf",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$Y_{0}\\left(z\\right)$"
],
"text/plain": [
"bessely(0, z)"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"Y0 = sympy.bessely(0, z)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "b64e87a5",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$- Y_{1}\\left(z\\right)$"
],
"text/plain": [
"-bessely(1, z)"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d1Y0 = Y0.diff(z)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "4aa5114d",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$- \\frac{Y_{0}\\left(z\\right)}{2} + \\frac{Y_{2}\\left(z\\right)}{2}$"
],
"text/plain": [
" bessely(0, z) bessely(2, z)\n",
"- ------------- + -------------\n",
" 2 2 "
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"d2Y0 = d1Y0.diff(z)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "9d7a5842",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$z^{2} \\left(- \\frac{Y_{0}\\left(z\\right)}{2} + \\frac{Y_{2}\\left(z\\right)}{2}\\right) + z^{2} Y_{0}\\left(z\\right) - z Y_{1}\\left(z\\right)$"
],
"text/plain": [
" 2 / bessely(0, z) bessely(2, z)\\ 2 \n",
"z *|- ------------- + -------------| + z *bessely(0, z) - z*bessely(1, z)\n",
" \\ 2 2 / "
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"odeY0 = z^2*d2Y0 + z*d1Y0 + z^2*Y0"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "1528c962",
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$0$"
],
"text/plain": [
"0"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"sympy.simplify(odeY0)"
]
},
{
"cell_type": "markdown",
"id": "15186d4d",
"metadata": {},
"source": [
"## 2. Plots of the Bessel Functions"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "5c9fa483",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.4400505857449335"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"besselj1(1)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "d1a1444c",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdeWBU1dkw8Ofce2fPPpNlsu8kYd/FgAJqAQEXELeiFfVVa30/q7WLbdXWVqttra3UaitCX7GuVVQQJKDIIjuEBBIghCxk35NJZp97z/fHDUMy904yEwIJM8/vr8lhcuYkTPLkOctzCKUUEEIIoWDFjPQAEEIIoZGEgRAhhFBQw0CIEEIoqGEgRAghFNQwECKEEApqGAgRQggFNQyECCGEghoGQoQQQkENAyFCCKGghoEQIYRQULsyAuHq1asrKip8fLLL5cK6cSPC6XSO9BCCkSAIgiCM9CiCEb7hRwSl1OVyDW+fV0Yg/PLLL8+cOePjkx0OBwbCEWGz2UZ6CMHI6XQO++8F5At8w48IQRCG/U+QKyMQIoQQQpcIBkKEEEJBDQMhQgihoIaBECGEUFDjhqWXmpqaoqKiiIiI2bNnyz7BarX+97//bWlpuf766ydMmOBuLykpKSgoiIqKWrFihVarHZbBIIQQQr4bhozwueeeGzdu3KOPPvqHP/xB9gkul2vevHnvvPNOc3Pz/Pnzv/jiC7F969atc+bMaWxs/PDDD2fPnu1wOC5+MAghhJBfhiEQ/uQnP+ns7Pzxj3/s7QmbNm3q6OjYvHnzSy+99Morr/z2t78V259//vkXX3zx5Zdf3rhxI8/zn3zyycUPBiGEEPLLMATCsLAwQsgATygoKFi4cKFCoQCApUuXHj16tKWlpbu7e+/evUuXLgUAlmVvvPHGgoKCix8MQggh5JfhWSMcWF1d3axZs8THUVFRarW6vr5eo9EQQmJjY8V2o9F49OhRbz00NTWtWbNm+/bt4ocGg2GABHRbjSsu1DFJP1BsRpeC3W5XKpUjPYqgY7fbCSFYROLywzf8iOB53uFwMIyvWZxCoRj0yZdj16jHTymllBDikUSKjd56YBhGp9NFnKfRaAZ4uQ4nWfUd7oZFCCHkk8uREcbHxzc1NYmP29vb7Xa70WhUq9WU0qampoSEBABobGw0Go3eeoiOjr7rrrsWLFjgy8vdns4/d4KpsrFjwjEpvKwcDodKpRrpUQQjQgimJpff0N7wW7ZsWbdu3aUYT8CLjY1dvXo1z/OEkOH9VXMJA+GpU6dSU1PVavX3vve9n//853/+8585jtu0adOUKVOio6MB4Oqrr964ceMjjzzC8/zmzZt/8YtfDNdL58fAviaKgRAhNNrs37+fEHLbbbeN9ECuMBaL5fHHH1+9evWl6HwYAuGePXtee+21srKy5ubm22+/fe7cuY8++qjT6czNzT106NC0adOWLFny0ksvLVq0aOrUqWvWrFm7dq34ic8+++xdd91VWVl5/PhxlmWXLVt28YMRTTaQonZcMkEIjUZjx45dsWLFSI/iCtPV1fX4449fos6HIRAmJyf3/U9NTU0FAI7jNmzYkJmZKT7esWPHp59+2tzcvGPHjvHjx4vPXLBgwZ49e7Zu3Tpu3Ljly5cPY6qbF0G+rseLaRBCCA1ueAJhcnKyRyMh5JZbbnF/qNFovv/970s/Ny8vLy8v7+LH4CErjJZ1DXuvCCGEAlBg7q5MCSH1FurCnBAhhNBgAjMQKhiI1ZA6Cy4TIoQQGkRgBkIASA6Bcz0jPQiEEEKjXsAGwgQtqTNjRogQQmgQl+NA/YhI0EGdZaQHgRBCo9v69etLSkoAgOO4rKysO+64Q61WX4oXevDBB6+55pp77723b+OmTZvsdvvy5csBwGw2/+lPfyouLh4zZswvfvGL8PDwSzEMWQGbERq1pAHXCBFCaEAbNmw4depUenp6dHT0G2+8IV6EcCmYTCabzebR+M0333z11Vfi41WrVh08ePChhx6qqKgQQ+NlE7AZoVEDha0jPQiEEBr1pk2b9tBDDwHAokWLcnJyzGazTqcDgLNnz37++edWq3Xx4sWTJk0CAEEQNmzYcOzYMYVCMXv27Pnz5wPAt99+++2331JK8/Ly7rjjDrHPrVu37t+/32AwrFy50pfcrrq6+vPPP6+vr9fr9XPnzo2Ojj527Jj4opdBwGaEsRrSbMOMECGEfGK32zdt2pSRkaHVagFgx44d1113ncPhCA8Pv/XWW8W87aWXXnr11VczMjKio6O3bt0KABs3bvyf//mfuLi45ORksQUAHn/88RdeeCE+Pr6iomLGjBk9PYNvXDxy5MiYMWP0ej0AqNXq6dOnHzx48BJ+tf0FbEYYo4Em60gPAiGEBrO+XHjnzGU69ZwZRt7IZz0aX3zxxb/85S9Wq5XjuI8//li8COipp57661//KtZFSUtL+/3vf79w4cKDBw+uXLnyvvvuc3/uwYMHFyxY8MgjjwDAAw88AAAlJSUffvhhRUWFGFAbGhreffdd8QkDaGxsjIqKcn+o1+sbGxuH5Uv2RQAHQtJkxYwQITTazTcSo8YzOF0iBrl9ML/85S9//etfA8COHTtuueWWI0eOpKSkFBcXv/LKK6+//joAWCyWsrIyAHj44YfvvvvudevW3XjjjQ8++GBCQsJdd921aNGir7/+evHixffdd9+4ceMKCwvtdvvNN98sdl5RUeG+d3YAOp3Obre7P7TZbOL07OURsIHQoIIOOwgUGLyCAiE0iiXoSMLl+50/kHnz5iUlJX3zzTcPPvigUqn829/+lpaWJv6TmCYuWrSooaFh586d//nPf6ZNm1ZRUZGXl1dRUXHgwIFPPvlk1qxZhYWFGo0mMzPzo48+cnfrSx3pxMTE6upq98W0VVVVSUlJl+arlBGwa4QcAyEK6HCM9DgQQugKUVxcfPbs2aysLIZhFixY8O6770ae53Q6AaClpUWtVi9YsOCtt97q6Ohoa2trbW1lWfbqq69+5ZVX0tLSysvL8/PzKyoqysrKxE8MCQlxOAb/RTxnzhxKqbjKeOjQoerq6oULF17yL/i8gM0IAcCgJq02qldhSogQQl69+uqr69ats9lsdrv9V7/61bx58wDg73//+5133jlhwoSMjIyKioqrrrrqn//85913393Z2ZmamlpSUvKDH/wgMTHxiSee2Lp1a05OTn19fUxMzLXXXqvRaNatW7d8+fK8vDyO40pKStasWXPDDTfIvjTP8yzLAoBSqXzttddWrlw5ceLEoqKiV199NSws7LJ9BwI6EKqgzQZw+Q5lIoTQFeadd94Rsz2WZfvGnvj4+J07d1ZXVzc2NiYnJ8fHxwNAQUFBeXl5a2trUlJSYmIiALz66qs1NTW1tbUGgyErK0v83JtuumnBggWnTp0SBCErKyskJAQA1qxZo1QqPV69srJy+vTp4uMVK1bMnz+/rKwsIyMjJibm0n/pFwRyINSroc1OATAjRAgheWKUkkUISU1NFa+YdbdkZWW5A54oKSlJup6nUqkmTpzYt8UjwysuLv7zn/+8c+fOV155xd2o1+tnzZrl/xdxsQJ2jRAA9CrS5lnHACGE0MhTKpVz5swpKiryCKsjIpAzwkgVtNsHfxpCCKHLLCcnJycnZ6RH0SuQA2GUirTb8SghQgh59cYbbxw7dgwA1Gp1bm7uypUrB5gsvRj33HPPvHnz7r///r6NH3/8sc1mu+eeewDg22+/PXDgQEVFxaOPPuoxrXqpBfLUaJQKj08ghNBAtm3bZjKZrr/++vHjx3/wwQeLFy++RC9kt9tdLpdH4759+3bt2iU+/v3vf19RUbFhw4aqqqpLNAZvAjkjxKlRhBAa1NixY1esWAEA1157bXZ2trvodklJyWeffWY2m2+++eaZM2cCgMvl+uCDDw4fPqxWq2fPnr1kyRIA2LJlyzfffCMIwtixY90J34YNG8Si26tWrTIYDL4MY/v27QCwc+fOS/RlDiCQM8JIJenEqVGEEPJBZ2fnRx99lJubK9YI3bp165IlS8LDwzMyMr7//e9//vnnAPDCCy+sXbt21qxZubm5hw8fBoBPP/30Jz/5yYQJE6ZOnVpYWCh29dBDD73xxhvjx483mUwzZswwmUwj+HX5IpAzwggVdOLUKEJodDPv32r+btPleS2FMTXy7p94NIpFt10uF6X0/fffF4uc/exnP/v73/8uzpQmJCQ8//zzN998c2Fh4S233OK+awkACgsL58yZs3LlSkLI3XffDQDFxcUbN26srKwUL/itrq5ev379j370o8vzBQ5NIAfCSCV04NQoQmh004yfpUhIvzyvxWhkNsK4i24fOnTouuuuO3jwYEZGRmlp6dNPP/3cc88BgN1uF++C+H//7//deeedb7311o033vjwww+np6ffc889S5cuTUtLu/HGG++///5p06YVFRX19PTMnj1b7LyhoaHvtRKjUyAHwggV6XTg1ChCaFRjdGFK3eUrJzaA6dOnJycn79q1a8yYMSqVau3atRkZGeI/iWni/Pnz6+vr9+3b95///GfmzJnV1dXZ2dmnTp0qKir65JNP5s6de/jwYZ1Ol52dvW3bNne3vhTdHlmBvEYYocSpUYQQ8tXBgwfLy8tzc3MJIYsXL167dm14eLhYO9tsNgNAXV0dx3Fz5sx59dVXu7u7Ozs76+vrAWDSpEnPP/98ampqdXX17Nmzq6uri4uLxU/UaDRW62i/GzaQA6GaBQJg9dyvixBC6IIXX3wxKioqLCxs+fLlL7744pw5cwBg9erVFRUVY8aMWbBgQXZ29ssvvwwADz300Lhx4xYvXjxhwoQf/vCH8fHxr732Wnp6+uLFiydOnJienj537tyYmJj33nvv3nvvzc/Pnzt37pgxY4qKiry9tMvl4rjeicklS5ZERUWVl5evXLkyKipK3IxzeQTy1CgAhCvB5ARNgH+VCCE0RO6i2xqNRtzeIoqJifnqq6/a2trq6uoSExPFdb4vv/yyvr6+ubk5KSlJr9cDwEsvvfTUU0/V1tbGxsYajUbxc7/3ve9VVlZWVlY6nc709HSx1vb69evFiyb6qqiocBcXff/99/seNAwNDb2EX3Z/AR4iwpWky0FjNVh3GyGEZAxcR0av14sBzy0+Pl68icLNYDBITwoyDONeXxR5rBQWFxf/9re//e6771avXi22XM7I5yGQp0YBIByXCRFCaPQJCQlZvnz5iRMn0tLSRnosgZ8RQhcGQoQQGmXS09PT0y/ToZFBBXhGGKYgJjxBgRBCyLsAD4ThSuhyjvQgEEIIjWIBHgjDlNCNU6MIIYS8C/BAGKoAE2aECCGEvAvwzTJhCtJsxTVChNBowXHcm2+++cUXX4z0QK4wPM+7j94Pu0APhEooH+0XgCCEgsgTTzxx4403jvQorkjR0dGXqOdAD4Q4NYoQGk1CQkKmTp060qNA/QT4GmGYEo9PIIQQGkiAB8JQBXRjRogQQsi7wA+EPRgIEUIIeRf4gRDXCBFCCA0g4AMh6XbiGiFCCCGvAj4Q4hohQgihgQR4INRy4BTAJYz0OBBCCI1WAR4IASCEgx7X4E9DCCEUnIIgEOIyIUIIIe+CIRDiCQqEEEJeBX4gxKOECCGEBhD4gTCEw42jCCGEvAqCQKggPbhGiBBCyIvhCYRvv/12SkpKRETEypUrzWazx7/u3r17Wn9FRUUA8O677/ZtrKmpGZbBeAhV4K5RhBBCXg3DNUwnTpx48sknv/7665ycnNtuu+35559/+eWX+z5h6tSpH330kfh4+/btzzzzTF5eHgA0NTWlpaW5nxwXF3fxg5HCzTIIIYQGMAyBcN26dcuWLZs2bRoA/PKXv1y+fPlLL71ECHE/QavVpqeni4+3bdt2zz33KBQK8cPQ0FD3P10iGAgRQggNYBimRk+fPj1+/Hjx8fjx41tbW9vb22Wf2dbWtnHjxvvvv9/d8uWXX2ZkZOTn569fv/7iRyJLx4EZp0YRQgh5MQwZYUdHR2hoqPg4LCwMANra2vR6vfSZ69evnzx5sjgvCgDz58+fM2eO0Wjcv3//Qw89pFKpbr/9dtmXKCkpWbhwofvDjIyMY8eOeRuPxWJxuVwM0xvjlZRrM9OeHn5IXxzyg9ls7jsTgC4Pu91OCFEqlSM9kKCDb/gRwfO8w+HgeV9/pavVao4bJNINQyDU6/Umk0l83NXVBQAGg0H2mf/+978fe+wx94eTJ08WHyQlJZ04ceKDDz7wFgjHjh27Zs2aBQsW+DIehmHUarU7EEZqhQYHDQlhfftq0NBRSkNCQkZ6FEFHoVBgIBwR+IYfEWIg1Gg0w9jnMEyNZmdnHz9+XHx8/Pjx6OjoyMhI6dMOHjxYXl7uLdQpFArfI7xfcI0QIYTQAIYhEK5atWrDhg0HDx40mUwvvPDC/fffL04XPPvss5s2bXI/be3atStWrBDnTkUff/xxZWWlxWL5+uuvV69efeutt178YKRCcI0QIYSQd8MQCMeOHfvaa6/deeed6enp8fHxzz77rNheU1PT0dEhPuZ5/syZMw899FDfT9y9e/e1115rNBqfeOKJZ5555r777rv4wUjhgXqEEEIDIJReAUFi4cKFTzzxhI9rhBaLpe8a4Z5G+otD/J6lw7AaigbW3d3t3jaFLhvcLDNS8A0/IkbpGuEoF6LAqVGEEEJeBX4g1OIaIUIIIe8CPxDqODDjrlGEEEJeBEEgVBCz6wpYB0UIITQiAj8Q4vEJhBBCAwj8QMgxwBKwYYU1hBBCcgI/EALW3UYIIeRdkARCYsYz9QghhOQERyDEo4QIIYS8CI5AyIEFAyFCCCE5QREINRgIEUIIeREUgRBPUCCEEPImKAKhFjfLIIQQ8iIoAiHW3UYIIeRNUARCrLuNEELIm6AIhLhrFCGEkDdBEQgxI0QIIeRNUARCHUcseAEFQgghOUERCLU4NYoQQsiLYAmEeDcvQgghWUERCHUcWPAaJoQQQnKCIhDi1ChCCCFvgiIQ6hRYWQYhhJC8oAiEmBEihBDyJlgCIZ4jRAghJCsoAiFWlkEIIeRNUARCDQtW3DWKEEJITlAEQi1WlkEIIeRFkARCPFCPEEJIXlAEQhULAoBLGOlxIIQQGn2CIhACLhMihBDyIlgCoQ5PUCCEEJITLIEQ98sghBCSFTyBEI8SIoQQkhFcgZDvaHbWV470WBBCCI0iwRIINRzYHELLP55uWf1T3tQx0sNBCCE0WgRLINRyANUlRKXRTMi3Hts10sNBCCE0WgRPICSq6mPq3Gnq3Gm204UjPRyEEEKjRbAEQg0L6sYyZUqOMjXHee70SA8HIYTQaBEsgVDLga61UpmQzkZEU4HHZUKEEEKiYAmEEdTMOSxsRDQAKOJSXE3nRnpECCGERoVgCYRGa70pLB4IAQAuJtHVXDvSI0IIITQqBEsgjDY3doQYxcecId7V1jCy40EIITRKBEsgjLQ0dWhjxcecIc7VioEQIYQQQPAEwnBzc6smRnzMRsXyHc0jOx6EEEKjRLAEQp2lpVUVLT5mI6JdGAgRQggBQPAEQo25rVl5PhCGRFC7lTodIzskhBBCo0GwBEJlT1ujIrL3A0LYMD3f1TaiI0IIITQqBEUgpLyLc5ibSLi7hQ3X8yYMhAghhIIjEArdnYIu3MITdwsTHoUZIUIIIQDghqWX6urqzz77TKFQrFixIjo62uNfW1tbd+zY4f5w5syZycnJ4uPjx48XFBRERkbecccdOp1uWAYjxXe3U12Utc/FvGxYlGBqv0QvhxBC6AoyDBlhSUnJpEmTzp49e+jQoUmTJjU2Nno84eTJkw8++ODH550711vebMuWLddee21ra+uGDRvy8/PtdvvFD0aWYOogYZEW/kILGxqJ5UYRQgjBsGSEf/7zn1etWvWXv/wFAG677bY333zzN7/5jcdzkpKSPvroI4/G3/3ud3/4wx8efvhhQRCmTJnyySef3H333Rc/Him+u5MLibD0yQiZ0EissoYQQgiGJSMsKChYunSp+HjJkiXbtm2TPqenp+edd9757LPP2tt7JyS7u7v37dsnfiLDMIsWLSooKLj4wcgSejq50HCLi7pb2NBIvrvzEr0cQgihK8jFZoQ8zzc1NcXFxYkfGo3G+vp6j+cwDJOcnLx79+6KiooHH3xw48aNs2bNamhoIITExMS4P/Ho0aPeXqWpqWnNmjXbt28XP4yOjn788ce9PdlutxNCGOZCjHd0tTOhUdYOcM++ulQal6n90k3GBie73a5UKkd6FEFHfMNTSgd/KhpW+IYfETzPOxyOvr/hB6ZQKAZ98sUGQkJI3x9CSikhxOM5+fn5u3btEh8/88wzTz755L59+zyeJvuJbizLhoSEREb2HgTU6XQDfGHMeRc6N3epEjIcNQCEYQgAgCIsSujp9P1biXzh8W1HlwfDMB5/+aHLA9/wI4JSOuzf+YsNhAzDxMbGNjY25uXlAUBDQ4PRaBzg+QsXLly9ejUAGI1GSmlTU1NCQsKgn2gwGO68884FCxb4MiSFQuH5J4ClWxkeqWaBZxQqDgCARhiopVuhUPjSIfKR+J0f6VEEHUEQCCH4nb/88A0/IhiGoZQO73d+GILqggULvvjiC/Hxxo0b3eGqpKTEZrMBgCAI7ifv2rUrMzMTAEJCQvLz88VP5Hl+8+bNPsa5IRDMJkYXrmHBen7jKFEogWGp3XqJXhEhhNCVYhh2jT711FOzZ8+22+3d3d2HDx9+8803AcDpdI4bN+7QoUPTpk378Y9/XFdXl56eXl5evmPHjg0bNoif+Nxzz91xxx1nz549ceKEUqlctmzZxQ9GlmA2MbpQLUcsLqpX9U7AsrowvqeTU2ku0YsihBC6IgxDRpibm1tcXDxhwoR58+YdO3ZM3P/CcdwXX3yRlZUFAD//+c+XL19uNBpXrFhRVlY2b9488RNvuOGGvXv3Jicn33vvvbt37750y868mBFy0PdMPRMSLphNl+gVEUIIXSmGp7JMQkLCD3/4w74thBD3mYqEhARvBwRzcnJycnKGZQzeUN4FLiej0mg5vt9RQl0YBkKEEEKBv+VJnBcFQvquEQIAow0TejAQIoRQsAuCQGjpJppQANBw0C8jDAkTLBgIEUIo2AVDIOxhtKEAoPVYI9Ti1ChCCKGgCITdYiDUsKRvlTVGFyqYu0duXAghhEaFIAmEIQCg4cDWb40wVLBiIEQIoWAX+IGQXsgIPTbLhAqWnhEbFkIIodEh8AOhYDUzGh0AaD02y2hxahQhhFBQBMIeRtM7NWr1CIQWDIQIIRTsgiMQimuELLHyfTbLaENwjRAhhFAQBEKLmchmhGqdYLMC3uKGEELBLfADIbWe3zXK9lsjBIZhVBrBah6pgSGEEBoNAj8QClazuEbosVkGcHYUIYRQUARCm4VRaUGcGuX7/RPR6ChmhAghFNwCPxBSaw/R9h6fsHpkhBo8SogQQsEu0AMhpYLDJmaE6v67RgGA0egEKwZChBAKagEeCAWrmVFpgRAQK8vIrBHi1ChCCAW1QA+ENjPRaMXH0jVCRo0ZIUIIBbsAD4TUambUOvGxdI2QaEIoBkKEEApuAR4IBbuFUZ/PCFlJRqjR4dQoQggFuQAPhNRmIeczQg3X7z5CEAOhDQMhQggFtQAPhILV7M4IZadGBQsGQoQQCmoBHgipzUI05zNC6dSoWkftGAgRQiioBXggFGwXMkKOAQLgEC78K6PR4hohQggFuYAPhBZyPhACgIYFW5+kkKixxBpCCAW7AA+E1G4Vy8qIPG9iws0yCCEU9AI8EHpmhByx9tk4yqi1eCUhQggFuQAPhLTPGiFI98swLFEoqcN2+QeGEEJolAjwQCjYrP0CofQCCrVWsFku97AQQgiNGgEeCKnNQlQa94fSu3mJGovLIIRQUAvwQCjYPXeNSqqsaSnul0EIoSAW4IGQnr+eXuSxWQbECyhwahQhhIJYoAdCu42oL0yNaljp1KgWjxIihFAwC+hASKngsDLKfmuEkiprWsGOGSFCCAWvQA6E1GEjnBKYC1+jun9lGRCLy+DUKEIIBbFADoSC3cr0mRcFPD6BEEJIIpADIbVbSZ+dMuBljRADIUIIBbNADoRC/0OEIO4a5T12jeLxCYQQCmqBHAipdGpUskbIqDAjRAihoBbggZAo+wVCucoyWtwsgxBCwSyQA6Gvm2Xs1ss6LIQQQqNJIAdCaUYoLbGGGSFCCAW5QA6Egt3at9AoAGg4sHiUWFMN/UC91QVdjqEPDyGE0GgQyIGQ2q2MymONkFg91wg1Q8sInQLM+Nw15mNnI06sIoTQlSzAA6HH8QlpZZneS+r9t6FKiNHA7enMP0r5wZ+NEEJotArkQCjYrUSl7tsiXSMEhiUcN4RL6jdU07szmJWZzKdVdPBnI4QQGq0CORBSe7+K29C7Ruj5NEY1lKRwZ4NwfQKZZiANFtqEs6MIIXTFCuxAaPOsLMN6Hp8AAKLSUD/3y1T3UABICSEMgatiyP5m4eJGihBCaMQEciCUTo1qJSXWQNwv4+dRwmNtdIqeiI+nGsiRVpwdRQihK1UgB0LqsDIqz+MT0oxwCBdQnOiA8VG9gXCSnhS393lRu7Xz47/bTh8dwoARQghdfgEdCO22wTfLABCV1t+p0dIOmhfZGwjHRpKSjgsZoWnrf1xtDe3rX8YSpgghdEUYnkD41ltvJSYmhoeH33XXXT09PR7/WlJSsnTp0ujo6PDw8BtvvLGiokJsf/vttzP6qK6uHpbBuEmPT3AMEABH/xW9IZygON1Fs8N6A2FGKKmzULsYXwXefHBbxIr/VWWMtxbtvoixI4QQukyGIRAeP378pz/96caNGxsaGrq6up5//nmPJ3R2di5btqy0tLS2ttZoNN51113u9ilTpmw7Lz4+/uIH05cgOVAPcvtliMrvM/XlJpoV3hsIOQaSdeRsNwUAR/VpNsLA6eM0E/JtJw4MfegIIYQul2EIhOvWrVu2bNnkyZO1Wu3TTz+9du1aSvttHsnPz1+1alV0dHRoaOhjjz1WWFgoCL1JWWhoaPp5CoXi4gfTF3XYiFLt0QT8IVMAACAASURBVKjlPGdH/a273WYHhkCU6kJLVjic6aIAYD97XJU5EQBUWRPtFSeA4iYahBAa7YYhEJaVlY0bN058PG7cuLa2tvb2dm9P/uqrr2bMmMEwva+7adOmlJSUq6666t///vfFj6Qfgae8iyiUHs1qjtj6lxv19/hEZTdNCyV9WzLCSEU3AICj+rQqJQcA2LAoRq1ztdQNbewIIYQuG+7iu+jo6AgNDRUfh4WFAUBbW5ter5c+c+/evX/84x+/+eYb8cPrrrtu7ty5RqPxwIEDDzzwgFqtvvPOO2VforS0dOHChe4P09LSioqKvI3HarU6nU7itBGFuru72+NfVUTVYrLp4UIsdALDm9qJ5JnenGpmEtRcd/eFJDJewZ1ug+5ul62mjLv+Lld3NwAQY6rpzHGFJtzHbgOAdHkYXQZ2u50QolR6/s2HLjV8w48InucdDofLJTkA4IVarR50unEYAqFerzeZTOLjrq4uADAYDNKnHTly5NZbb33vvfcmT54stkyaNEl8cOuttxYVFX300UfeAmFeXt5bb721YMECX8bDsqxarabdLrNa447Qbjqli1HpQvukdObwSEdHo/SZ3jTzQkYEDQ29MDeaYxB2twk6hWCymsOTM4EQAICUMUJ7ve/dBoZg+3pHA6VSiYFwpOAb/vITA6FG47n/42IMw9TomDFjiouLxcfFxcUxMTGRkZEezykuLl6yZMm//vWvRYsWyXbCsqx74XBYULtVukAI4hqhZLOMX2uENWaapOs3NZoSSs71gKuxRhGb1BsFAThjqrNhmPfBIoQQGnbDEAhXrVr12Wef7d27t7Oz84UXXrj//vsJIQDw61//euPGjQBw6tSp6667btWqVYmJiUeOHDly5IiY1X7wwQdnz541mUwFBQWrV69etmzZxQ/GTZCcnRBJjxIyKi315/hEnRkSdf1aknXkXA91NtdwsUnuRkVcsqvpnF9jRgghdPkNw9RoXl7e66+//oMf/MBkMi1duvTZZ58V2xsbG8WZ0tLS0pSUlIKCgoKCAvGfCgoKoqKi9u/f//TTT5tMpqSkpN/+9rf33nvvxQ/GjTps0rMTAKDhiNVFAS6kdETt34H6WjNN1PX7AyJSBTwFc1Od2pDgbuSiYnlTO3U6pBt2EEIIjR7DEAgBYOXKlStXrvRoXLNmjfhg2bJlstneX//617/+9a/DMgApapc5OwEAGtbzAgrGz6nROgsk6DwbE3XE3FQfMuXqPv2ybFSsq61BEZfie+cIIYQus4AtsUYdnvXVRNJzhESl8X1qlAI0WalRSzzaE3XAtzRyBmPfRs4Q72qp92PQCCGELruADYSCw8bIZoSSzTKM2o+MsNUGoQpQSr5tiTrCdjVyUXF9GzmDkW9r8H3MHr6qped68Eg+QghdWsMzNToKeds1Kt0s41fR7QaLTDoIAGkqK3E5mJB+pwY5vXHIZ+r/USq8clxwCnDiNi5smEvuIIQQuiBgM0Lq8LyVVyTNCIlCSXme8j4dz2y0Qpzc8ZVUV0unNsajkY2KdbU3+TjgvngKLxYJG25grzGSt07hrb8IIXQJBW4g9JoRytzNy6g11G7zpdsmK43TyGSEia6WFpVnGQEuKpYfUiDc2UDjtTAhijw4hnnvbL9A6Gqq6fzkH66mmiF0ixBCSCpgA6HgsDNym2U0nOeuUfBndrTJCrFyGWGMraVWEe3RyEbF8B3NvnTrYUuNsDSZAYA5caSqmza4hybwrWufp3Zr69u/BUFysyJCCCH/BWwgHGiNUHpJvUrt436ZRguNlVsjDLe3nyOe5VUZtQ4ABJvZl5772tlI5xkJALAErjUyOxt6k0JryQFWFxZ590/YkAjrif3+dosQQkgqcAOhlzVC6fEJ6M0IfQqEzVaIkQmvoOlpLWcM0i2ebEQ039HiS89uVheUdtBp0b3hNj+W7G3u7dhauEs77ToA0M68wVK4069uEUIIyQroQKhUSdu9TI2qfQ2ENhqtlskIqanVpIpqk6wzspEx/gbConaaE0HUbO+HM6LJoRYKAECp7fRR9diZAKDOm2E/fRSGtTorQggFp4ANhILXyjJiibV+GJVW8O1MfbMVYuTWCPmuNhpmaLR69syG6/muVp9GfN6xNjpJfyHWTtKT4+2Up+BsqGJ0YWy4HgDY0Eg23OCoK/erZ4QQQlIBGwi91RqVnxpV+3o3b7PNSyA0tbNhkU2SYMpGGPhO/wLh8Q46IepCIAxVgFFLyrqoo+qkKi3P3a5MH+uoPOlXzx66HBfz2QghFCACNxB62ywjNzXqe7nRVrmpUcq7qN2qCwtvkmaEEQZ/M8KSDjoust9LjI8iJzqo49xpRXK2u1GZkuOoPu1XzxcGDLDiaz72P87XS3FyFSEU7AI3EDrs3naN2oa6WcbkBAUD7tU7N8HUzoRGxmiJTEYYpue72n0eNQDAyU6aG9EvEI6NhNIOcNSeVSZmuhuViZmO2iFOjX5RLZSb6Inl3HNHeI8x28sKzfu/wtVHhFDwCOBA6GWNUFJZBsSpUR/WCGXTQRDnRUMjYjVEJiMM1/OmNl8HDdBhB5sLjNp+jbkRpKzD5WquURhT3Y1cbBLf0UQddt87d3vjpPCzCUxmGFmexqwruxDzHJUl7e/+yXygwLRl/RC6RQihK1HABkLBYZM/UC+5hgkAGKVGcAweCFusYJA7O8F3d7JhUdEakMkIw6P4Lj8CYVkXzQ73jLXZ4cTUWM+G6fuGdsJyXHSi0/+7fzvssL+Z3pzCAMD3M5mPKy8Ewq7N68NvelC/6tc9320SzCZ/e0YIoStRgAZCSqlTfmpUy8mUWPPxJqZWu3wg7J0aVZNmSUbI6MKpzeJjIVMAKDfRLLlAyLXUKOKSPdoVcSmuxmofe3bbXi/MiSNaDgDg6hhS1U0brQAArvYmZ0OVZvI1bFiUOm+65ei3/vaMEEJXosAMhNTpICwHRGYa0/tmmcF3jbbZqEElNzXa3cmGRcZooFkaTAlhQiKE7g7fBg5nuyE91LMxTAE5rhpblDQQJjsb/c4Id9TT6+J7/985Bq41MjvqBQCwHd+nGTeTsBwAaCdfYy3+zt+e+/qkUthSg3dIIYSuAAEaCL0sEML5zTIev6GJb0W32+ygl80IezqZkIgYNTTL9cGGRfEmX/fLVJhoRphMrB3H1zeHxHs0crFJruZaH3t2+66J5sdeeIk5cWRPEwUA2+kj6tzpYqMqa7LjXJmPRQakXisRfnNU+NFe/sOK/ptuBMFavNfVipcVI4RGkcAMhIKX6+kBgCGgZMDmcROTb0W3W21UL58RdrAhETEa0mKTyYHY8CjfN45WdtO0UJmXSHXUVaskgTA6wdns3zUUPU44a+p3YH9WDNnfTEEQHBWlqswJYiNRqpRJmfbKUr86F3U74feF/KfXs/+Zx/78oODqEwo7Pnm9e9v7zX99cgjxGyGELpHADITUy/X0IumZeh+LbrfZvKwRdncyoRFaDhiAbqfnvzKhkb5PjVb1QGqITHuMueE0a/Ro5AzxfHuTX0cdjrXRcVFE0ee/fZKenO6k5rpKJlzP6MLc7aq0cY7KEt97dvuwQpgTx2SFk1kxJFEHW2p7h+dsqLKd2B/92B/Drr+9a+PaIfSMEEKXQsAGQm9TowCg4TyrrBGV1pfNMt6mRvmeTiYkAgCi5ZJCNjTSx6lRpwBNVpqo88wIBZtF6bKVOCM92olCyYRGuDr8uPKwsI1O1vfrX81CZhipPnlKlZrTt12ZOsQD+x9XCHdn9L7Eykzmo4reb4h5zyZd/hKi0ujyl9grTgztgipRWRd95gh/vB3XIBFCwyAwAyF4OU0v0kr2yxDfMsJ2O42SmxoVerrYkHAAiFHL7JdhwiJ53zLCWjON0xBO8n/CtzXwkXGVcrc5cfp4V2uDL52LPOq3iSbpianqjCIpu2+jIjnbca4MqH/BxuyCfc10QWLv17A0mWypEXgKIAiWoj26adcBAFEoNRNnWwp3+dWzW6sN5n3Jd9rhhi2uOnO/4fXs+rz1X884qi6q8hxCKNgEZiAcJCNkpVOjWurDOcI2G0RJL7QQeGqziJOK0WpolV5AERrl49TouR5IlpsXdbU3KfVxVd0y/8QZjLxfgbCdjo/0DIQTogjXeFaZkN63kQ2NJAqly8+8bU8jnWIgIYreDxN0JE5LCtuoo/okG6Fno2LEds34WbaSA3717PZyEb8slay+mn1gDPN84YVpYWvxdz17NmomzG57+3mhp3NonSOEglBgBsIBNsuAXHEZolBSl3PQxbZ2O+glgZDvMRFNiHhUw6CWnRoN57t9+r18zkxTQuQ247Q16mLjWmzUKRkgp49ztfkaCCnAyU46VhIIx0fQyK4aRXyaR7siMdNZd9bHzkW7GoW5xn79zzOSnQ3UdrpQPWaqu1GVOdFRWz6EXal2HtaVCU9NYADgx+PYjyqE3tLhlHZ9+X+Rtz2mu2qBZvI13Ts+8bdnN5cALxwTntzPd0qLkgu8vykyQmj0C8xASB32ATbLaFiweJQbJYRRDn6UsE1ualQw986LAkC0GlokGSHjc0ZYa4ZEnUy7q71JERVr1JJasyTK6uNcbb6uEdabqZaDSEksz6MNrVyE9B5jRXyas77Kx85Fe5vo1bH93lT5seS7JmovL1ZlTXQ3EoVSmZRlr/B7M85XtcKEKCL+uRCthnnxzOfVAgDYz54gLKvKngQAIdfeYj5Q4HsRAw+/Osx/Uy+YnHDH1/16sJ85Vv/MXY2/u28I1XwQQqNZoAZCG1HI3Mor0norNzrgUUI7DzwF96Sfm9DTxbgDoYa0ymSEET5mhDU9NEmyUwYA+I5mNjI2WQfVPZ7/xOnj+HZfA+GpLhgjKVsDAJGd58rUydIESGFMdTZU+tg5AAgUjrbSGdH9XuKqGHK4yemoOaNMze3brsoY5/A/EH5RTW9NvfCmvTWVfFZNAcBauFM7db7YyOmNiphE++lCfzsHgHM9dG2Z8MF87p+z2UYrbD5fE0CwWdrX/1G/6tehN9zZvv7li8kLLS54/6xQ0Y2ZJUKjReAGwoGmRolFcjcvUWkHzgjb7RCplGkXzCZG1xsIDXIZIVFpgFLqGPzAvveMsJmLikkOIed6PIfNRcW52hsH7Vl0upPKBkK+6VxnePKpTs/OFXEpLn8q15R1k1gtiej/XUoNJemWKhoRy6j7lRJXpubaq/3e1bK1ji5MvPAlLExkdtQLLp5aS/ZrJlztbldPyLee2Odv5wDw5knh3iwmWg0sgZ9NYF4r6Z06MH+3STVmiipzgm7WIsIprMeH0jkA2HmY96VrbZkw6wtXUf9dr6at7zW+8ICp4P2h9YwQGrLADYQD7BqVbJYBAEalHnjJqt1OI2VP0/d0MSG9x+8MKpBmhOBzUlhrljk7Ab0ZYUySDmokG0cZXRh1OnxcbCuTK2QKAM6mc7whSRoIuegEV1sDCJJvlhfFHYzH2QzRQihvM2R7NCqTxzjPnfErtTrTRVkCfb+EaDWkhpJjp6sJq+CiE9zt6txp9lNHfO9ZRAHeO0tXZff+UCxPYw610EYrAKXmfVtCrrlZbA+55hbz3i/97Vz0txLBqCXbFnF/nME+vOfCN9Zy+BvrsV1R9/7CevRby5Fvhta56ONK4XeFgvQPMsFmESxyG64QCnoBGgjtAwZCuXKjg15J2GGX2zLaf43QoCbSXaMAwIRECD4EwjoLTZAEQuqwUZeD0YUlhcisEQIhbGSMj3s7y7toVphMu6u5VhOXWNYlyZIVSjbc4PvxjOIOMkkuEE6xl5eFZno0MrowRhfqV7m1XY30mjjP/q81ktrior4LkACgiE2mvMv3bUSiI61Uw4H7VmQ1C4uSmM+rBUdNGbAKZVKW2K4Zf7Wj+rTQ0+VX5wDgEODV4/yL0xkAuDeLMTvhm3oKAJR3dW3+d+QdP1YmZUXe+eOuL/895AXOv5cKzx0Rasx07iZX3ze57dSRxt/d1/j7VeYDBUPrGaEAFpiBUHAOuFlGvu62WhjwTL3XQ4R9pkZlj08AABMSMehRQocAXQ6Iloya72hhI6IBIElHaiRTowDARcXw7T4FwjMmkM0IXc11hsTEM3LXLnGxSU6fy6GVdBHpIUUASO6uOMBlSNsViZmOmjM+dg4A+5rp1bGe/efHEnLuhCpjvEe7KmO8/ewJ3zsHgC01dHFSv/6XJJHNNdR6fF/feVeiUKpyplj9P/6x6ZyQG0HyIggAEIAf5jJrTgsAYDuxn4uKVabmAIAyNZeLirWd2O9v5wDQZIXfHOE3fo/912x2gp78qbh3k7Fg6e5478/6B56LefK1ro1vX0x9u3+UCrp/Oyd86qrqv8ZprzjR/s5L5r2bh9wzQiMoMAPhIFOjsptlBjtK2OGQ2W8J4mYZXe+FEbLHJwCADR38Aop6M43TEEYSR1yd5wNhiMzUKACwkTG+VGnhKZzrkSlkynd3AMumxoSVm2RGrohO8P33ZmkXO96z+g2AIOjaz213eV6dAQCKhAxnvR+bcQ4005nRnuO/KoYktZYq0/M82pXpYx1+1krdXifckNDvJ+KGRGZng2A9edhdjlykyZ1uO3nIr84B4MMK+v3MC/3fkcFsqREsLrAc/lo3c4G7XTdzoeXQdn87B4C/l/J3ZTBi0fbfT2VeL+XFP/h6dm5Qj5ulSh/LGeJD568wffXuEDoHgN2N9MUiofQ27r4s5ravefdtZs7a8ra1v1dlTujZs7Fn52dD6xwABAp/LBaWbHVtrfV8K9rLi837tviy0I7QEARoILQPuFmGlbuScLBL6jvsXgJhn4wwQgUWF0hP+zEh4bx5kJm0BivEa2Xa+c5WLsIAAAlyxycAgI2M9mVqtNZMDWqiZiX9t9Zz0fGZYaTcJLNex8Uk+hgIOx3Q46RJknOQrrYGLjzyjE0t/eNDEZ/mrK/wpXMAMLugqoeOl2Sc8a52BXXVq+I82pUpOX6VmLG6oLCNzu4/9apXwQRNj7OlQZkypm+7aswU+5kivxY4HQIU1Ao3pVz4idOrYKqB7Kg0288UqcfPcrerx8+ylxf7e8iSp7D2NP1hXm//GWFkRjT5pEqgvMu8d0vovOViu+7qG22nDvMmX4vf9vWzg/wrM5mUEPLkeEbDwftne9/oHf99PeKmB3RX36h/4FlTwXtD6xwAnjvKf1Et3J3J3LfT9V3The9tz87POt5/1VZ6sGX1T6nDPrTOAaDJCj8/yL9wTPB4K1Knw1q816+/yVCACdBAOEitUW9XEg4cCGmkUm4ni9nkrlVNAKJU0Cb5UWVDB18jrDPTeNmdMp2tbIQBAPRqsPIyI+ciYviOloE7B4AKuZsOAcDVUs8Z4sOVoGZl6sNx0fGulrpBOweA0g46JpxKvwBnQ5XSmJoZRk5J1iAVxlTfzykea6O5Ef3KhYsc58rqIrIPt0o6j09ztTf5nkMcaqV5kUTHebavICWN0bniNY1ubLie0YU5G3wdPAB810RzIojH1PeNSUzZsUJlSg6jvrBdmFFrlWl5ttNHfe8cAHY3UqMWxHlX0T1ZzPtnBfupI1x0gnsnEaPWasbnW4/u8KtzANjfTNvssCKt9z/gV5PYv54QAMBeflywmrXTrwcATm/UTp3Xs2soSeGZLvrPk8Kn13N3ZzCv57OPfscLFADA1VJn2vZ+9GMv6x94jotOMG39zxA6B4BuJ1y7yWXjobCNLv/aJZx/v1CHrWX1T3t2bWh981fmg9uG1rnok0rhBzv5LyXXcDpry80HtgpmubUHNDoEZiAEp50o/TxHqNIMfBOT16nRPoEQAPQqmaOETEjEoEW/6i1eMsKuVnFqlAAkaEmd9Ex9ZDTfOXgg9HbBk6utgTPEA0BGGDkrmR3lDPE+7mc52UnHhMlkSM6Gai4uJS+SlHZIz37ECtYewSY34StR2EanyO3EcdaUOY2ZhW2SFJ/lFHEpjlpfK+PsbaJXx8j0P9NSeliXK21XpY+zV/ixBrm9Trg+wbP/7yUS5dmjqtxpHu3qnGk2P3e9bqgSlqX2+3FenMR810i7Cvdoplzbt10z+RpL0R6/OgeAdWXCA9mMe+r+ewmkzQ5F7dS8b3PI7CXuS7BD5txkOVjg+05jtz8dFx4by8ZoAACWpTJqFjadEwDAVPBe6NxlbGQMAITf8pB5/1dD2KYEAL89yufHkr/NYj+cz3bY4Z0zvels15f/x8UkRD/2p+jH/tj1+VuuNl8PI3lYc1r4xSHhqhjy6Hf9ruG0HPmm9a3n7KeONv35sYspNF9uosu280sLXCf6/xzxpo6uTeu6v/kvdUorISFfBWYgFBz2gQ/Uy2eEQ50aZfsEQoN8udEIvnuwqVELNWq9ZIThevFxgg7qJMHa90CYLnflr6u1gdMbASAjlEhPebMR0YKl25e86nQXzQ6VCYSupnOKuOTcCHJScjwDCFHEpbgaqgftHACK2uhEuUDoqKsITckslGSEAKBMzHT6vBnnYAu9Si4QGltPb4Ycabu/a5A7G+g8o+eP29hIMqnjWGfiRI921ZjJ/hYE2FxDlyT3G3+IAvJjwVJyUDPuqn6dZ010NZ7zqxyrS4ANVcJdGRf6ZwjclU4+PW2xlR7STpnrbueiE1i90ebn4Lud8N9K4ZGcC9+f/x3LvHFSEHo6bSUHdPmLxUY2LEoz/mrz/q/86hwAmq2wrkx4cToLACyBV2ayvysUeAp8Z4vl8NcRtz4CAFxMYsg1twwt42ywwNOH+E0L2B/mMp/dwD6+j2+3AwDwHc2dn75p+OEfon7wdMjsJe3v/2UInQNAqw2u28zPjiWLk5gbNrtqzv81LPR0tfztCep0OCpLW996zq8b2fqy8bDyW55923ntJldjn9+ClHd1ffavppcfMX317kUWFyxso1trqUsyQL67w/d7yy+dwAyEdMBdo16OT2gGOz4hMzVKnQ4Q+L7FyQzqIWaEDRYwetY4AwDgu9rYcIP4WD4jDDfwXW2Dvk2ruuVvOnS1NXCGOABIC4UK6TEzQjh9nC9/Jp/qhKwwmZ9DZ9M5RWxyTjickvtLQBGb5Gzy6W7h4x10otyWVGfd2ZSs9GK5HyVFYobvtVIPtnjWxAEAyruYpoqTIVnSUgaqlBxH1SkfO7e6oKhdJtAKpnYD7dlNUz1HHptMHTZf/r4RVfWA2SWzgHqXqrxDESHOKLgRllNlT7Kd9CPj3NNEU0NJcv8F4GVpTGPRYWXKmL4zIgCgnTTHWrTb984BYEOVcE0cE9Pn/b88lTnQQpv2f6seN4vRXHjj6mYtHMIE5ttlwvI0JvZ8/1fHklgNbK4RenZv1E2/3j3+kGtusp3YP4Q1zj8V8/dmMWK1isl6clMK87cTPACYtr4XMnuJIi4ZAELn3SaY2v1N9EW/OMTflkaeHM88ksv8KI99fF/vD1rnhjc1k+ZE3Pqw/v5ngNKe3Z8PoXMAeGwvb+fB9APFXCO5ZZvLHa46Pvybq7U+8q4n7eVFXRvfHlrnAPCbo/wt2/jnjvLzNrvMfX73dm//sOkPDzW9/Ihpy/ohdz4sAjQQDrZZRq6yjGbgXaOdDoiQVJYRLN2Mtt/Km14ts0bIhIQPOp/TYPGyRtjVykb0ZoSJOqiVzCMSTsFoQgY9sF/VQ1Plpkb5tiY2Kg4A0sNIpVzdL1bv0+zo6S6aKc0IKXW11HExibkRRHpgHwC4uBSXD6U7BQolHTLlwgVLN7VZU5PiTE7aLvm2KxIyHHU+bcZptILNJfP9cTVUc1GxebE66RokF5MoWLp9XPg53ErzIohWsgDpqCgxGXP3SCfMCFGmj/X9+Mc3DTA/XrrjGGZ1H9sZMln6fPWYKbYyP9YgN9cIi5M8f1dMNZDJrQd7Mq/yaFePu8pWetCvBOK/lfT29H7D13CwOIlpPbSzb7oJAMqUHBB4p88z3qL/KxMeyO43/gdzmHdO85ZD29zpJgAwmhDNhHx/t+yaXfB/Z4Qnxl3o/6nxzD9PCTZTl7V4T8i1t57vnQm97o6ebz/1q3MAqOymn1cLz07u3ef21ATmUAs91kad9ZX28uKwBSsBAAiJuO3R7u0fDWFj7b5muq2OrruG1XHwm6lsCAdviad6Th91VJZG3fu0Mjlbv+oZy+FvnLXl/nYOAF/W0PVn6JFbuH03cemh5In9vdPm1mO7zQcKYp/+V9wv11iL9gxtp/RwCdBAONjxieGaGvVYIITe4jKeT2NDwgVL98C/GhqsECfJCKnLSW0W967UeC2pt8jFqojBZ0dlM0LqsAl2CxsWBQBpofKBkDMYBz1T7xLgXA9N03l+uqujmdGGEaU6M4xUdFPJXl1QxCb6UsO6qofqVSRc8oeIs75SEZ9KCBkbSU5I1iAVcSmullpfDqcfbaVTDHLzrjVlyqSsKQY4Ip16JUSRlOk4VzZo5wCwv5nOkpyABAB71UldRt6+Zplvuyotz1Hpay3WPU3E49IPUUTd8b268dLNxqrsSfYzRT52DgBba/tVthMRSuf2HN0ZOtWjndMbiSbU91zc4oKdDTKBdoWhQ9le61EqAQA0E2dbj+/1eexwrI26KMzsn44vT2V6ThdCRGzfgkQAoJ1+neXot753DgCfVAr5sUzfmlDZ4SQvghzZtk0zblbfP5Q1k69x1Ja7fC4OLHq9VLg/m3G/+dUs/O9Y5rUSoWfnhpBrbnZvhlDEJivTx5oP+h1OflfIPzuZEasoE4CXZ7B/OCY4BTBtfid8yX1i/4wuLHTB3V1b/D544xLgif38P2ezBjUQgL9fzX55jh5ro9Tp6Pzsn1Hf/ykbGsnowqLu+VnXxrcF2yDXHlw6gRsIBy66LVnLH7yyjIN6CYT9MkKDmrRJjxIyLFFrBy5w1WihcRrp1FkHExbp3omQoIN6ubcKG2HgO1sH6NwpQIuNxkvWIF3tTWxkrNh/WghUyl55qI/jxYrSsAAAIABJREFUB5sareqhRi1RSc5muJpruZgEANBwEKMmVZJAy8Uk+XI8o7QD8qQnFAGcjdWcMRUAxkeSE5IL64lSxUZE+7LrtbCNytbEcdSWKxIzp+jJUdk1yKQsR61Pa5AHW2ROQAKAo+pkck7emS5qlgRrZWqe3eep1++aYY6k5g4IvKP6lCp93O5GmT1QQKmPlXearXCuh06TjN/ZWM2qNF90x0o/RT1msu+7Xr9toFMMMn/lXG06sitkkpV65tHqsTP9qmawoUpYluo5+HAl3GPfezYp36NdlT5O6O7wq+DR+2eFlZme/d+TxTDFO7QzbujbSDiFdtI11iN+bNl1CfBuufA/Of1+Ua/KZgrO9liK9+quWti3PSR/iXnfFt87B4CyLnq0la7sc7x1qoFkhsH2faWCpVszYba7XTfje85zp33cQ+7230ohXgvubWKhCvjpBOalIsG8/ytlcrZYRAIAFAkZqqzJ5u+GWLnw4gViIBR4KghEIVchGwC8ZYTqQa5h6rBDhGSNULB0M9r+GaHX4jLhvPdlQqcAXQ4wSMvKdLWyYXr3h/Fya4QAwIbrB84Ia800TkM4yf82397ERfX+IkvUkSar3JWHBuOga4RnusBb8TZFdKL4ODscyiTTw5w+ju9qG3TDW0kHHRsht0DYUKUwpgJAXiQpkZt6VRjTfDnkUOQlEDprzyoSMycbyDHJrlQQCwLU+DRZdLiVTpdbgHTWV+pSMsdFEelmH0Vihqu51petgA1W6HFBtqRmkKP2LBsVOyUxtO+ZPDffb//Y1SjMiWNYybfHfqZIO2bSrkZB2rsqe7LvGWdBnbAgUeYXETlzpCZ+6s4GSTqbmsu3Nw1aqsltUw1dmizpn9KZ7Yf+q5kpeVWiHneV1efKPp0O2NdMF0v6vymkSWNp51M8Cx5pplxrOebHAur2epoRRjL7b3OLVsOPYH9L/HiP6ShV1kTB2u3XqZ5/lwn3ZjEef8I+lMN07d2im7XI/Sc4ABCFUjvjBn93Kv29VPjxuH7fnAfGMNtr+c6dX7jPtopC5y/v2bNxyPt9LlIABkLqdDDez07AAJtlvE+N2nigANI1Hrk1QtJml4tVunDB+8bRJiuNkSsrw5va3FtGYYCMMDKa7xooI6zugRTZnTJ9AiHHgFFLaqSbcaIGv/u33EQzZYu3tdaLGSEAZIWTM9LiNQzLRsUO+gf4yU6aK1kgBABX4zlFbDIAjI0kJZKpUQBQGFNcDVUDdw4Axe1yO3EEwdlYrYhPS9IRpwBNkneHMjHT6cMaZIsNOu00Q7Jl19VQzenjiEoz3UAOSQIh4RSKuBSnDxnngRaYYZA5wemoOqlMy7s6lshOvSpT8+y+7Xrd3UivkZt3tZ89HpkzIUol851XZYxzVJ30sWLq13X0unhJ/5Tay45F5E39ul5an4JVZY63l/m0MbXJClXddJZkm5KzvkKlVr/XaZSGcc3YGbbSg750DgBbaoRrjYz0+KmybN/xuJnb6iXfmbSxgqnd99nRjyuE29Nkfksv7tr9efg1nq2EaCddYz22y8fOKcD7FfSeTM/+b05wTWja1zNuvke7dsb1liM7fF/9LeuiFd10Sf9J71AF/DjsVCfPKdPG9m1XJGSwYVG200PZTHTxAjEQDrhACAAaudsnBp4ald0pA3JrhHq5NUIQ98uYvWaEjVaIld8y2t43EMZpSKNV5m3Ihhv4zjZvnQNAjZkmS2q+AADf3sRGxbg/TA2BKsnsKBcVy3c0D/zuP9tNM2QPKbbUiYcUASA7jJyRnKkHAIUPxWtOddFc2YywsZqLSwaAPNnjGQAKY4qzcZDjGVYX1JhlLqhytdazoRHi7VETokiRZOqVi07gezoHPQd5tJVO1hOZQFVTJhbynmogMmuQAMrkLIcPGeehVjJDdoGz6qQqJWeynpzupNK//JRpeT4e/9jTRPOlC5yU2s+eUKWPzY8l0qlXRhPC6o2+7K1osUGtWWaB1llfwejCrso07JBkhCBmnGU+ZZzb64S5RkY6F2I7dTRk7LRoNZGeQFVlTnRUn/Zx18lmSX3a3v5LDihzr5IergdC1HnTbb5N7fIUNtUIN6dI5hLsVn1j6T/5qdIZdc2E2dZiXxdQD7VQDQvSzcZw5nBLRMZnbREezYrYZEYX5uPfTwDw/ll6Z7rMN39557efRnpGWRAD7eGLunplyAIyEA50mh683EfIqNQDVJbptFPpvCjIZYQGuV2jAMAMeJSwSW6nDADwpnZxJ4tIxUKYQm4zTrh+kIywG5JlM8KOZvGosig1lFRJzgkQpYrR6AY+61PeRTNlp0Zb6t2bETLDiWw5Uy46YdCFh1NyNymKJ1LY0EgAMGrBKch8Z7jYFOdgVyqWdtLscJl5Y2dDpcKYJj6eoCfHJYEQCFEYUwYtjlPY5m0nzhlFYiYATPESCBVJWY6awTfjHG6F6dEy7Y5zp5UpOSoWxkbKrHH2Vt4ZrJCb2QWnO2VKGbha64lSyUZEz44je2WnXtPH2n2Yet3dKOTHEpl517MnVJkTpkeT8i7aJZkeVmVOtJcXD9o5AHxTL5duAtjPHFNlT74+gXwtSdqISqNMyvRlyy4F2FYnSLcRUbvVca5s8oyJ0qKpAKDOnW47ediXwR9qobEaIt3MbDt1RJWWlxur/brOM11WpowRzF0+rv5+ViXcKlk9BQDr8X3s+PwNVTKzlJoJs63F3/nSOQB8XCncni75uRL48PLvNobNPi35s1g7cY6t9NCIVAYIwEAIzoFO04O3qVGlmroc3maoOx0QIX+avluSEcqcIwRx46j3cqONFhone5re1M70CYTgZZlw0M0yXjPCtiYu6kKVzpQQqJY9QREVxw94/W9FN0in/kDg+c4WVt/bf1YYyF9wMVggbLCAmpW5A8vZVMPFJLk/zAmXOaHBRcfzHc3U5Ryg/xNyBzMAwFlfpYg/Hwij5AIhgCI+fdByqccGXIAEgNwIcq5Hbr+MDwUBBArHOsiUKM+xCWaTYO7mYhIBYEa03NQryyni0we9/eNQC52ol9kG5agsVaXmAsCsGPlA6GPGuaeRzo6T+S1kP3tclTFewcD0aCJd41TEJgl2K9810CyI6NsGOk8SCCnvsleWqjLHz4snO6RTrwCqrMn2M8cG7fxEJwlXeh6vBAB7ebEyeUy2QQMA0l/3quzJjooTvswbf1UrLJJEWQCwnTykzp1+YxLzlTTQEuJ7oJVfPRUEW8mBifmzDjRTk+TnRjN+ls23K6/PdNFOh+dmXQCwnz3BRcVelR27oUp63jpckZjpb3HBYRGAgZA67AMcIgQANQsuATy38hPCKNWCl6OEXqdGJRmh17rbuoGOEjZ6yQiF/hkhABi10CAtLhNu4LtaB5i9rDHTJPlCps1s5IVsIiWEVPfIfDoXFTvAqoZAoUruXgtXexMbrndX6UwNJbVmuc04gwXC011UuhMEAFxNNYrYC4FQtngNYTlOHzdw/yUddJxsIGyoFHfiAMA4ueMZ0Fs3fJBizUXeFyCVCekAoGAgL4IUS6de41Jc7U0Dl5k+Y6JRSpm/EnrTTUIAYHo0OdwiO/WaPejxD9kbPwDAXn1KmZoHADkRpMNBpVVqlam5vhQ9/0523hXAUVmiSh8LAPmxZE+j5E1DiCo9z372+MCd15ip2UVzJJPqznOnFTGJjFp3TRyzr1mm3Ikqa6L9zOAZ565mVjbdtJUdU4+ZDADXxZOv6yS/7rUhXEySo3Lwb05BLZXdRmQ7dUSdO21hIpHPOHOm+VItocZMGy0yRSQc506zEYZQQ/SsWLJdknEq4tMoz7t8KILxxTm6NFlmRcB6Yr96/KybUpgvquUyzvGzbMd9CrTDKwADITjtMGBGCF7qbg+wX6bDy/X00kBIACJVID3cPfCu0UYLjZWcnQCxrIwkI5QeJSRKFVEoBatcEAMAgHM9kKTzbKROh2A1/3/23jSwjeu8Gn7ubNhJgAQJ7jspSqSofbcsW7ItW45jx04c20kTZ3PaNOkX57X7dUvauMmbuq6bNmkaN0uTNvvmxE7sKF4VWZKtXZREUtz3BSBBEAuxzGDmvj+GAoG5d0BQkrOYc35Rl8PLKwCcZ57lnJO+f40DDdMsD9nCkqRfNxBORLFLoLDFkzOTqQYhAAgMlFkp+3NF5dmHZbqDlLooAEi+MTXjUdHspOh6g8rQyPp32xHALVRuxuRwKhC2uNClORoPsrQm+5BeNAkjEdq92DfG5hemNInWF1IGUxHLccWV2QPt6Rm8oZCyLo72CFULTsKb3IgUBAA1EC5Vej0xjcmHegAQhy+pjhwIYGsROkEEWnUIK7tTZiwJHQG8magbJ2cmgWHVov2uEuYYddintmXJqdfXpvB1JRSdgUT/goFlgQmq7ZQ2oVDVJPlGl6wbv+ZjyHQTABJ9C07Re8uQTo9zfaJviR5nUISLAYoBp+QdAYbhisrXuFBcBlIW0bRqg9h/YcmM84UxfHM5Q7468Uunzc2bAeC2CkantLs5fmnpjPP5UeUAtXvaecKyZuueUtQ1h8lehrllW7xreWoM1wRvwUCoiHHEZ8sIAcBGDYRmqx6DIltGaNN6OlCro6zdmSUj9OoNyxAZYZne4Gh+turoaIRSGpXnptl8d/qEdJUN9DLCLKXRgTDUUn0tZiZTdVEVDXnQR1RH2bwCJR7LwqXtntPJCH0ZGWGzE7qp8zIlVdk5+51zGaYNKrAkysEZrmghkNs48FjQANHj5Etrk5PDWf5uOwJ4VT7FNEOaGODL6lL/XF+I2mkMDWEplbizM3g9URcFAOlyAxIAmp1oIkrptAlVTdLI0qVRCvFDEpO+cb58wWx5axE6MU15uheqm8WR7iybn57BLS5kIQV3hrqE2gWDyW1F6MwMLWmrbUkMLVF6PTKFr6PqGPRfFOpb1a93l1CGfRDHC5WN2adCFAyvT6PribquEg3Lfq86BnVDKaLTSxralqSXvDaFtxVTfNMSve3mpg0AgAD20jNOB1dUJg5ne+UB4MVxfDOhAg8A8UunzM0bAeCWCvQCsTnkpggfkeDkNN5bpn1xkjMTWIzzZXUCAzeUMi8QGSdXWIIsDvGKJGyuBm/BQIjFOGNaIiO00udldOVGg1mmRq3aIECdl8musuaNUXqEWE7iREzTgyyzokm6uIxbb14mKAIGIAnLcmA6fVIGACrtaCJKyXvYrKXRgRCuo4q3zUykZ4SgY3ABCHHuUllfvKYniJvyKevJzIxwVT6iyplynsosKm4RCXwxWl3XO8q5y4BZvA9RxWsYqx2ZzFlssM7PUiRAQdXEKV8MhG0FlNIoAPBldWLWHuRZP15fQFmXxvqFy4GQRdBWQMk4OXeZMh/MovPgjUE0SXlzpbE+vqQKcbz6zy20jBAAhOpV2QPhcR2hc3Go01S9QLXOF6DaTnlx+MqG5PRE9tnOo15KRgUYi8NdpsuBdpeH0oMElQGSdV7mYgAXmSnPr4mBDqF2tfrJqbAhB0+p2JvqWsSxvuxTIa9OKqRKOwAkettTaju6GWfjuuw9TgzwygTFDkWJz0uTwyqxYbUTSQqQA26mxjZxsDN73/3QJN5aRDE1i186Y27epD5831qpU9pdvTlxRYqsV4O3YCAEScxOnwBdKqFVrzQ6J+L83KZGQdeJaYkeIfkXpYQDjN0JmWX2UqteRlioNzswptMgTAZ8nCtj3FBgwG2mlF65Ao+sHwgHw7iOOjLqn+Tcpekrqv0veWX2NmFPiMIWx5IoB/1cWsZZ50ATUZwgiDFccaWkT8+4FMSr8ikMTsk7rGolp9Digk5abZsvrcnC0LgYwG30QDiQmsQBgLYCdDGAydyBr6jPTlU8RwuESjSsRCPpTyEbCtEZMuNEiK+oz0JyODVDEZQBAHG0V6hsSv1zSxFzapqSFAtVq7LnJcd9lB4VAIjD3ULNovXV9mL0BlEdRSzHl9Vm6XGGJBgI4w3EmJI0NczY8hn7Ajdglwcd9dLS2Zo1iawSd69N4R1uyg+K/RdMda2pf9IzTsHMl1Znf0r47SSmyOZhnOg7b2pYCIQ3lKJDk9Rhn3XZM87zs9hlQhXEbSHRd8FUuyb1iLOvjDJVy1jsnKcqewP4xXHlpnJaFO8+bVq1Uf36lnL04jglXTav2nhl0uRXg7dgIFySRwjZxGV0AmGCkhFiMQEIkRI2bjP4iedUxpaXRW7UF6P0COXQLJunbV7pyo3qB8LReaikcSfkTO6Eimo7jBDVUdZVnMXgYjAMdKfDmUkuszRan0czuMjqeigpMBqhJCVJ/yTrKk7P2DgGqu2UQMsXVyR943qH75rTYShODnOe6vQVXc5+SXWWNuGFWZ1JnLSRVADIF6DARLHB4stqpckhvcOPRLDAQrFZ+11pvJ8vq01/hNpQSBGvAQC+vCELVfH0DN0DUhzt4S83IAGgxAIWjqKfx1c2SOP9WbRCqA1ILImSdzRVdwWAbcXoOC3jNFU3Z7kdn/DhjYU0J+fBzlTdFQCq7IhnKIUKoXa1ONqbpdN2xIu30wJhYqBDqFukilMDIah+lvoZZ0iCniBd1o6xOlLc4loHEhjUQw6m1rWKo71ZMs5DE/hGqkhC7zlTQ1vqn3vL0KtEIIQFrdpsGefLE5hMN0GRE30XzE3rU4e3sFQ1hrXiWF+WBu1MnG5lczW4NoFwcHDwySef/PKXv+z10vOGaDT6jW984/HHHz99OiPUnzt37vHHH//6178eiegOeiwbUmLJQEjn1Au6pdEcrSdUUE3qEcshk4VahorLEJcp+8tBbYMQAMp0M0LdHuHYPCYf/YBWGgWAKjtlngVxPGN16AXagTCtNIpx0j+lKY3WOWil0ay63kNhXEZXMR3n0+qiKlblI3JaHZksjNWuV73sCtADYdI7qskI1zgp3sIAwC2VEa6lZGwRHJ/nMl98KkODMdtYW54eLYw+jwogjvXzFfXpKxvdlJEQWBDHyRYIN9EYkFJmRggAm2nzOIzFzuQVSD76pJI3BmGJIrgjjfXxnqr058ttReg4dV6mujlLxvmGT6fuOnzJVJPhtLyjmCK+w5htXGFJFm7M0Sm8o0gbCLGYkCaHhKpVqZVdHnSEGgjrW0V9Y+djXrylCAnE7TnRf9FU35q+cn0JOkxmnCYLX5It4zxETTcBEn3n01XObyxFhyZpSVvjuri+oMFUDCajlFxcHOllC4pTuTgA3FRO6XEiwSRUNWV5SvjqJfy9wWucwl2D7S5evLhhw4bx8fH29vb169dPTmr/aJPJ5O7du59++un5+flbb7315z//ubr+/PPP7927NxwOP/fcczt37kwkso2J5w68FI8Q9Eqj+nKjcyJ2ElOjeoGQakkI+oOjvhgupo2MKuEAQwRCjwVNxyg1tCw9wtGIHndimnW6NYvVdvq8TJY24VAEaoiXQY7MIcGU7tQIAPV5aDBMyW6ylEZ7Q9BIbRBOZzQIVTTlQzfd9bBK73bcNQertQIaAGmaNSk0O1FviDo4Wp3UCYS+GEgKkFrn0uQgV1qjKXq3FQDdVbGsVq862u4HOkNxvD89owKANU40EMZx4uGPr6hfIiMkAiFOxJKzPs1TwuYiHXGcykZJp3p5agZvdtMEd0a6heqMKLvGhSaieI4c9qluFod1dcnf8CnUQJgY6hI0gdBDKb2CygDRITmMzmNRwXV27U+JI918WW16FG/KR9EkJg1AhJo14tAlvXT5tSmFoqIOIA5cUOddU7i+FB2mtgkb1urFEgXDa1PKDcQkizIfSo35qKiyIxtH6XEKtS3SxIBeg/bVCeV6ujjtOVPj+vSVfWXoFb2pWn0JvZcn8B7PNR4rvQaB8Iknnvjwhz/8r//6r9/4xjd279791FNPaS549tlno9Hos88++9hjj33xi1987LHH1PXHHnvs8ccf/9znPvf000+zLPvTn/706g8DALAUjxAArByal8hnQN0eYVCkDJtQR0Yhi7iMLZ9qX6err0aMjAIAz4DLBNOkuExeoRyki7/olUaTAR9HZoQ2RJrQAgDnKpYDlEAYl8Efx+XEvV72T6nG9+mwcZDHwwRxU+DcpXpJT28QN5BUfYCkb1yTbgLAqnxKmQgAOH0Vt0tBCrcBJyV5blqzv5WDEgslo+U9VZJ3lFq97JjTYyguEjNSWFuALlBLr+V1egwKvYxQGh8QMgOhiYWmfJpTladSDvqphRBvDGI0j0ZxfIAvrUkvSkM2lbgmPc7+yWmFnEcFAHGkJz2jAgAWwUY3OklUR9mCYgAsBygMDQxwnFZ3VaIRJejnSzKK3tSMEFQnLJ3B1GNevKOYcvMUBztNaXVXAEAAuzwMOY/D2PJYZ6HeO/vaFL7OQ9UZuCjUaTPC1+jDPmtFHZ7lxQB2mxFJXE4MXBRq12je2RtK0SFSfIcXhIoGvalaqogBAMR7282Zplo3lDKHJxXy4dLcuCGu0+OMSHBuFnYW/eEFwpdeeun22xfMLW+//fYXX9T6R7/44ov79+/nOE694Ny5c9PT0+Fw+Pjx4+oPMgxz2223kT94ZVDEOLNURmjjIEqXG6VnhLqB0EIJhIUmoIvLOPIVmn2uLwbFtMAth2ZV/TANSmltQtapq7KmWxqdm2ZdWm2uaj0qoU5GOBTGlXbKsEnSr20QqqjPQ/1EeZjNdyvRMPUBszekw52YHkuXlVHRlI+oDAquuIJKJRQVGI7gRlJldHqcKyxJSQGk0OJCnSRn32Rh7flUg44Ls7iVGqjSqPoprNUVr9Hl7LfP4nVERoglMemf1KSzALCO5nEBDMuX1lDti8/66RKp0mhvah41hU1udGaG8iyQRSWOSswAAHGkh88MhKCK49A1AeiDqQMhbGERmYuLw5f4yiZgMu57qhwrzQlrdUKnB/m6jzaPCpAY7NRoSQPATo+O+E5da4JWHU3IcNZPqeuqD4uaP6vGfBRLYlIrX6hdIw5fovY4D+upqPddMDVo7TL2lKLfUku7DWv1VO5enaQ0IHFSEocvCZl13WILVNopj1BCVaPsn6KmDa9N4U2FFP+Dq8TV7ifLstfrLSlZeG9KSkomJrRTDxMTEzt27FC/drlcZrN5fHzcYrEghIqLi1M/eOqULklzamrqG9/4xksvLXhO5ufnP/roo7pHSsSSiM1eaDUhFIpDItMmQmF5Oein/mAggawga74jhuawyUpen8ei6SiQ69jiEOf8DLE+HkZFJsr1yaCfbVhPrpeY0XBQWaMpy3BmnIjF5yOpia8URiNMMS9qtsHRMHCCiBFkfqNUgOEwQ3kR8grF0V5yvTeAqq0Lh08kEoKw8LwQnxoDZxF5fY0N9cwq21zk031JdGLBXDAdPXPo5hLtOwUAkm9cyXdr9q+1oO4gohy+oCRx7gi5fikIFVYGS9pvJEb7GXc5eX2TA52fhttKiKf74sroaJ/JoU3fz8+gNhckEtqbkTg+wLfu0uxfbYKxeSYwn9D8kWN3hTg+QB4mIsHEPFNtFhOJBEIIXw5D0lgfW1gqygrIGT+yNh+dmYZELXH4kurYcDeUa2PbiSnU5qR8LGPD3UJdq2Y9D4GVZS7NJOoynwyxu0KaGErEopo8AwBOTqP/2Kr9m1KiESUcUJzFmv3XO+GHg0wioR3ZR+X1sYFOpnmrZv3IBGxxUz7G0f6LbEUDud7iRG9MitcVZ744jkIlHovNTDHE8+jRKeaJzTj9Aw8AgLE41GW7+881+28tQD/sR+THgK1sinef5rfdqll/fRqtykO8oj1movscX7uGPPzOYuaVUfG+2sxVhmcKSuYHu/jMbi4AHBpn3laJyfPE+8477vqoZv8dhfCXx1E8kdRENqZ6deQ33zPfRNzNohCIo0ab9p2VBju44koJcZobzvXF6MURZV0eUSWqXRPpOmVaqzWMfHEUXV+ME4kEw+SaxXEcx7LElIHmmhz30gNCKP2PUFEU8lcyDIPTHhYxxizLav4b1B9cPCXH2e12l2vh4+h0OrO8CsKmfaaymuwvk42HaBI0sgqMxSqLceoPBiVwmZDmOzgWYW0O8voiC8yK2s0BgHM4lfkQeb0vDiVWyvVKeI7LLyCvL7WCN065nrE7ITLHFGiNUsejUGnXbpMMzrDOInLzageMzgO5zhZ4EheOkevD81DrWDgMwzCpC3DAy9e1kNfXOWBwnnJ4trAEB3xMGrVORW8ImvK11yuxCJYSXH6hps3msYKCYVZkNM6OvKdSnh4jD9MbhuZ8yn9Wnh7jS6rI9TUueHmCcni+pErxjTKt2zXrnUF4oJ64HmPZO8qXVmv2FxhoyoNLIWZzZt+WcZcpkTkkJTQN184QNOcDzzIKwyCEUrsp3mGurJY8/PpC+NkI5fBCRYM00k1efy4Ad1dTrpcnBoTr7ySv3+SG9jmmQdPQtdhYV7HiG+PKMu7To/PAMVBBiDxIE/1ceT1D3Aq2FcMjpyjvlKm6ef7lH5Hrp/ywtYj2zo71WnYcINe3F8PxabieaMsJNavl0R6udUf6YiwJXXOwsRAlY0z6VknfKDLbuHzt89AmN3QHISovuMAvbl7bEjn4HfIwx3yw00M5vDTUxdesIdd3e+DYNPNAvWYZhLqW5OAiI1MFBjjqg3/egrR/U/F5xT/JVzZqbnM1DrBw0BdGqzLfWaFmjTw1jJJartoRH+wuAZb42Ej9F4SGteTh95bBf3XDX5Ify4a2ZP9Fy7rdmvVDXvjXzTj9VrMkEKWuocXVBkKGYTwez+Tk5Jo1awBgamqqtFTbGSorK5uaWigc+f3+RCJRWlpqNpsxxlNTUxUVFeoPptJKEm63+7777tu/f38uRzI3tJnM5uwvk12QExjxmbPVSatDEhM8r82oEjLIipRv0a6jRJS155PXe2zgT0jkOudwJf2T5PqMKNc5tIcBACUSMLmKWeL6CrvsTVCuZ51FaD7IezJGSIIiIJBqOoLRAAAgAElEQVQKrdpNkuEA5yomD1PIA8dIYYXXyFei4vLI3DR5/UhUrstfOAzP86kLlDmfqegW8voml/L8KOZ57Z1OKK7As1Oa60UFpmJSg0v7fxUnp7miMl6gaBysciYHomypprNVWIITMVaWVE+lFHojypoCymHwzISpdTt5+DY3/o8umee1fzWmstp4z1nN9Riga05aX6TdRp6bRoLJ5KQIo7UVyl1htIOgUfMl1TA9xmc2nzpDyno35nlWURSEUOrXzE8NmyobycNv8sDFgMRyvDYuVzfFjh8kr2+fTf7fLSzPZ1yNxYQ8O2WubCDrxluKlXMBfF8D8c5WNSlTg3zm/MvZgLKliPLKxycGTNWryMPUO0HBklfkNEV+tnb13Fg/z7KaaudJf/Leeu3hAWNppKfgPY+Qf1M7S5Qf9lPOY6pdLY/28BsynP/e8OPWAjnPwoWTfPpRxdFeU30LeXieh/WFybNz7N7MzhlfUomwwoQD6VZoAPDGjPy+RsofuDTclbfnLnL/PWX4m72Uj6WlsS168mXN9d1BbObkOqf24nhvD1+1SjBTphVuKJOPzqBWd+Z5eF4or1fGes2XeYEqDvvkfeWUwycHOxz77iUPf2M5fOA1CbO8ZkTW1rzR/7//pLnen4CBsLS9hMNJhdzqanANeoS33nrrs88+q3797LPP3nrrQqZ//vz5WCymXvD8889LkqResHnzZrfbbbfbr7vuumeeeQYAZFl+7rnnbrvttqs/TI6w8RRlGWSyUHuE1AYhZKVPhCRC1HuBU0/rEcahmPz4YayE5xiCRwgqlZDqU09rE2ZpEHJEg1BFlY3SJmRdxfLcNDkSMqSrr0YZlgE9cRkAtqiMnJcZCOEKG4UKlpxeNL7XgD4vgxBXVE7Oy3QFcDNVxdS74PerQXM+6glSRnY5muvh2Dy2cDTTjIkhskGoYlltwna//qRMmnhbCk4B3GaKNzJfWpv0jWmaSYEETMcp3VlpYoD3VJFREFRFU2obr7KRpL2rI6PkxeSkTApbixhSv4ax2FmnW/PiJ2S4OEshfiSnxxmLjdp331GM3vDRmOk1qxPE4OgxL95JJWYMdmjmUVPQ068Rals0tH0McMyrkI4cSiSohAN8ZmKtYl0hGp/H5ICeUNeaGOjQ/M0ensRk4gsqcSKNQZiOPSXoEF2/po1sE1KJGVgSxdFeU522ewoALhM05aMTxLASX1anRIIaytahCeU6D+WGcPW4Bls++uij3/3udx966KH777//zJkzH/3oRwFAkqR169Z1dHQAwIEDBzwezy233PKpT33q0Ucf/exnP6v+4D/8wz98+tOffvjhh2+77TaLxXLXXXdd/WFyhJWlDMvoEeqDNO4EACjRMLJSxjEZBE6BorvN6lgSemkeTEosgngT2fAD1YCCNtzK0aiEumx62qSMimoHZXAU8QJjsZOuhINhij4ZTkrKfJDkZgBAvS6VsCw5re0u94awDndi0fheAyqVEAB4D0V6m+73qyjJ6XHOo53EAQA7D0UWNEgyxz3VSd+YZhS+IwA6I6NDmqnFFPSF1miBcJamWYOxRrwtHRuo0t68wBWWJiczYsm5WbyukDIDJY718Wnj9enY5EZn/FR9GYqiqX4g7OartD0tFfpCbloSxTk/bsqnqMCLQ11CZp0whSo7YhDlnRWqVkkTA5qnhGM+vENnUsZETMqo0NOvMdVppcM7A9hloo50dgg1q4FW5WMR7PCg1wiPDtbhYuz5GrWH307hPTlPyqi4oRT9lqpfQyimjkRwRMJriI+9ONjJl9Vqavsp7KWSKBAyNa7TkChemqDol14TXINNV61adfHixc2bN+/fv//cuXNFRUUAwHHcc88919jYCAAsy7788ssf+9jH6urqjhw5cuDAAfUH9+3b98Ybb9TX13/4wx8+fPiwQKt0vUnQ4RHaME36OShBPi0LV6IRakYIuiprTmpGSFXclkMBUlZGhV5GyOQXKgTnXZ9NT5GVUVFN86kHAJYmvU0NhPKsl813A606XWwBSQGSE0alEvaFoJHOnRjj3HqBELppQmhccaUmI8QA3XMU7kTSP8nmFZKCQSrWOIEyOCqYWEeBJqPV1ZSZGuJoz/UAsNaVq+KoguECbWQ0OetFJrNGnzaFDYV0Wj1f0aCROdbTlKGOjKooMkMeT6OXlNdL3pF0aUoMcIo2MioH/ViWOaLJrWJbsR6tfpUmEOpR6RPDl/QyNgDY6aGQKJDJwhWVSWlUSwzwulchM0JlPqSE5/Ry/V0lzBs+CgmVHBw94sVUBmFi4KIpkziRjutLGB02oTZp+y0tI1TiUck7opeL1ziQiaWxCWvWSJND6dybVyfxDaUUu4947zkNcSIde8sY0mEYAMxNG+I9Gfo1dMGaa4FrE11LS0sfeuihBx98MDXPghA6cOBAfv7C87zJZHrXu9718Y9/vLk544msqanp4x//+L333ms2L8H8u7awcjBPaMYiE51QP5egl0ZxLMLQMkLQU1nTkRv1xnCxmZyUmWWIKUQV5XoGFE53kl4apVycDMywTp2M0I5GaIGWdCUMSZBQoIh460iV0XTU0aqjnLNImQ9qRKF6ghRuA6j0Bp2MsNmpmxFqOPVj89jOI/KdlaZGqOmgihYX6ghQ1vnSaikzqboY0OFOTFC4EypKrcAgILkxfHltcjJDaG0gjAtMiFQjksYHNFT6dGxw04XWhIoGaSwjaTtDo9JD1owQADbTjA8RL3BF5ekZbU8QF5iQm/jYiMPdgk46CABb3OjMDCWWmKpXi0OZgVDPOmqoS6ihZ4Sg7zBsyhQd7QliG4/KSZXOgQ6hppmasQFAoQnKrZSnHKG8Tg5MpwtOHZ6kB0KxX0ulT4cuyaE+g+SgChCSxFxx4KJQvYpaf1JxYylFaw3xgkYC5uVxvJfqS9VzztS0QW/z3SXorB9HiBuyqWlDouds6mM/HMEhkS5hf/V4C2qN5gJ1alQDxmylZoRUWRnQ5xGCKi5DTPxT5UaTCoQlSicpS0bosaCZOM0jIt8tz5EZIehlhJxuIMxVXIaaDgJAksamT4Euvc0wbIFHozjaF8L0jHB6nNPpETbkocEwxbWH9KC4NAfNNE2ZpHeEJ0h4KegpjnKlNdJkRvXyIs3vF8vJ5PS4XmkUFqqj2kXGYkdWe3rGqet6P94nVOgHwkKgSG8D8JVaxdEzNHE1LInJ6Qlqm0rFFjc6qUerT2P7ndBlEHbrlS4BwGWCMhvlxefLapIBnxKfT6284cM7SAnTRCw5M5nlKUGf7dcipgXCLE7CAq0HlsL1pTTRUYYVapoTadXRw1PUjG1e8o3pFY0BYLMb9QVxgGjHmBraEv0XUvecVyd0lNV621NC3lTQq5eEBMzLE5h0KlZikaR3JEsubuVgs5siFMe5SxEnpBrAqoHimxIGV2wgtHJonjYso9CUZYLiskujbjNMEztR5UZ9cVxoBvLt1WPTAwCLoNAM3hjJqXcrREZI11dTZCUcYGmDiwBQbacIKAMAV1giZ9LGh8K4hpiABwCZcCJMB9WVEAA4t9ahtydI0VeTQ7PACXq5uImFChtFvZorKk/6p0BZ7Ax3BjBpQwgAkneE0w9Uulb1mQ69MobuOdxCBMLk9DjrLNKruwLAugIdY8LyDGPCc368nurHO6YVV0tHmRVxCMgGsFBeL00Np16ckAQTUUr3VBrv15uUUbFFj/ZenWFDccKnn7FV06tzKugSMAwrVDamksKJKI5IlEKCOHxJKK/PcviNbtQbwiEyL6lrTQx0pmKJrsdh1tIlAOymiYICgKmhTbyctPWHMAYg9VfF/otCdXOWw/MM7PCgw2SbML+QseenFFNfmaRnbPHecyktbCr2lTGvTtAkYJo2xC8Hws45bGIph0/0nhfSHC2o2F/B/GaMVh1t3pRyovjNON5f8SbFwZUaCOnGvCyHWJaUbJ8TwUkq1SiKkogxFlrZMYvKmsMpZ4rLeGNA9aanCo2mUG5F4/PaRTa/UA7NajLOsXkgS6Ny0M84nCTHWUU1TXcbaBmhriWvfypLaVTfjCljXiaWBF8MVxOBNukb43XqoiqanXCJbONxPOssSt9fz3ciOZUtI1ztRD1BSsapCYS9QVxqpfixSRODWTIqAFhXiM5R24Tl9dLYYptQlX0hL5PG+nidHp4KqigoMlk4V1GqtHtmBrcVIFIrUhzpyVK6VDc/66e8OJp5luPTNPclRZHG+pYMhPSkrXZNyoZCFT8jXxqNKQRlEwY2u9HrxP6s082YzKkG9ms0WRYsxqXJYb0em4o9JegwTcDaVL82fjkQ6mphE+JkJG4sZV6hOUWYG9epnTYM8OqEQgZCJRKU/d7s72ypFcptNAmYyiY56Fdn6F4Yw7foOf1mUixI7K9AB2nehKbVm+NdJwFAUuCVCeVmmrXTNcEKDYRWDkhFJdBRWQuJOI+n5PuM2arfEqDrbrPEvEw2b3qdjBAAym0UlTXEcozFrgm01GGZZGCaddFHEgCgyAKxJJAle66wRPZrS6NUS97kzARVX02FXiDkM+dl+kK41kG5HWv8eEk06zj08p5KKW1wtHOOMt4GGEveUSp3QoWVg3IbhYTAFVfIAV/qKepCQMePd3xAb6RTxXqagy4A8OX1YppNxFlaaVQOB7Ak6g2bqNATBeWrVqVEQfXFz3RHOlXk8VBtpyim8p4qJTSr1kJiSegMUCZxpMlB1ulmLPREX4Ve9dJU25KqLuqVLhMDHdTx/XTsLqHMXgKAcLnTNj6PgyLl+Skx0CFUNGRJ9AGg3IZcJgo9hq9qkmcm1BfnFR13pHiaGa8ebipHL9EM5U1NGxPdZwDg/CzOFxD5ZBnvOWtqWKv3WJzCLeU0H12GMTWtV5O2g2MKNWOLXzplWr05++brClFEgl6iu29uWi+N9Cjx6BEvbsxD1FvlNcEKDYQ2nUDImC0K0SYMSjShUf1JGQAoMsMMTZmdnJfx0SZlQO0REhIVKeiaMTndcnDRbygoAlC96eemOReF26ACqdVRKpUwOJNeXRwIYUpGiHFy1stm6xFCH1Uau6g8Ob042NkT1FEZ9VFURtOx2qXjl1RSLXkX51moGWHSP8Xa8/XmvFVQq6OI5Th3WaqfcX4Wt1GN4ycGhKwZYbMTjUZoupdpDrqTURBlIO9o0ugS6SAAbHYz9OplZWOqjacbCIe7s2dsALCtiObkwDB8VZNavTw9g9e4kIXIlRODnVnaSCpWO9FsAk8RTQehdrU00q1+Mo9M4euIHhuWk9JIt5CpSEBiTylDJ8xdDoSHJvH1tKlIjXuRHm4spfjcIpYT6loSvecwwCsTCqXHFgnKs0tkbACwoRBNxymio6amdeJQF5bE3+hmbKdNzZuWPPxtlczzo7Tq5eot8c4T80l43Yv3ERmbNDmEEJvlyVIFAnhbFfrVKDmPbRZqWxKXTv9yWLmj6k2MVis3EJKlUQBAtHkZqiuvHptehduMpnUyQjm3jFAJzWYpjZZZ0TidU59BJdTlTsx69bgTKmocFAYFYjnG4UqmGfsNhKGOaAnI4QDiTRoNl4zD21BYArIZw6kOupfRE4ImGolQ8o1p1HM0WO2kjHoDAO+pSk4tzMuoN1PylZcmh0i9Uw3aCoBOckiziWj3A92YnjBI0m7CwBoXpU3Iuoqxoqj84jN+uk2gONq75O1yaxE6RVPHTvf2o5YuF+gBS93Rdngo1UVQnRwGOwHgmI5ctUiTq9aAQbDLwxwhkjbGYmcLS8XR3pAE3UFKFJdG+9jC0uzpJgDsKEbnZylPIabGdYm+dsCYKicN6lRk1mETFTeXo5eoPIFVG+PdZy/MYoeASLuPePcZU2Pbkhkbg+DmcuYgEUsYs40vr0/0tf96VLm1krjhY5zoOmVevWXJw+8uQZfmsI94CrGs2ZLoOfvSSGJrMcoj+oDxi2+YW7RKsFS8vZr5xRDlxbGs3RE7f/QXw/jO6jerQQgrNhBaOYqyDAAwZislI6RaT8xnD4Q6GaHDqWh7hBQ2PQDI4UDW0iiM69nzpgWqkYieNz3FkjcdNQ5KRghqdfQylRADDEdoJEL/JOfWrYsCAAJozEdkGYTNL1QSsRQtqYdmkAQASd/oEhmhE12aozkhlC5ayXcGKJMsACBN6bLdU6AOdgIAX1aXmmehTnXKoQCWZT3WSgq6lkaVTSrJ4fQMxbscAKTRHkGf26Ci2AJ5AuXF58vrkjMTOBGbiOJokuKXKw5f4qua9HoBKezMJqFyEbINm3RoLGepuL4U/VbPfq/v/JEpHT/bvnY92ZR0WDnY6KbY+3EFHiRYpKnhl8ZpU5HRcNI3JtQukc4CwL5y5sgUxRXS3Lw53nXq4BjeT83Yuk6am5coLao4UImeIwIhAJhbtgUvHD8zg/cSUVwc7mbs+dnL6SoEBvZXML8cIZ5C7E6+pLrj1Lk7qynRJHb+qGXtzlwOf1MZuhjAk8RtzdK2K9J1ygHim0ScULFCA6GNp/AIAQCZbGQgpLvyZi2NZgmEcjiDhuaNUfTVcFLCiViWQFtu08kIXUXpKmuj83RL3mTAp6evpqLWQVHZAACucNFKfmIe5/FAzoMkpydIp0ANGvN0hdBSbL/uObyKKI1iSZSD/iwNSABwCuAQ0CgRyDlPVXJmQi2g6bHdkxNDfFlN9sOv02njCeX1qoOuPwEhifKIII33ZU8HVeh7+y1olZ2cVvRkWZYsXQLA9mJK9RKxHF9WJ470qGx0SvWPcNqjYrUThSRMNrBNtWvE0T5Zkqj6YepA75IfG1AJbfRAuC7R2/7KhHIjodQKAPGec+ampTM2ALilnHmRmrQ1bxw9exIDkA9nie6zQn1rlpHOFJwCtBVSAjnnqUQMe75z8ACZsSlKvOuUeU1OSdWBSubQpELWuiyt2yPtr+8tpVSkYxeOWtbu0K7q4O4a9NNBWg+1bXd+z2t31xCNhplJOTiThf6YDhMLb69mfjxABtr8kfym/8Mfz/GQV4YVGggFBjCARLynVCohlT6RhUQIAEX6pVGyR0hOjSqhAONwZXn6LrcCOTUKRGl0dB5XUukNAV+WYRkAqLXDIFVcxl2SMt4bCFNGpUGdlFnqjrbKCT3UeZY0B93uICUQJqfHucLSJctELU64SNDeES+w+W7JNw46JD9YcArM1sMDgLo8NJfApIQeX14njfcDxmdn8LoCmvf6aJ9QuUQPD7KQEKqa1OoltYcnz/oA0JLpJmRhjteuTgx2HvXinTRLWHGwY8keGwAggN0lDHmvRyYL76m6dLGLbgnbdz7H2+X6QjQZpeQNpsa2xGDn4dHEPiKpwmJCGunOpXQJ+uOL5tWb5y6eupU2DBLrPJFLaVHFHVXMs0RSBQC4eXv5yOukn21i4CJX4KGqFZJwmWBrEfo10cnjiiv8jP2DFsJbEePYuSOkw4Mebq9ijk9TqqNHiq+7KXi8lNfmFtHTr1jWX09VmKLiTxqYb/VoDy8q8A3r3hsmXspxkyvDCg2EoDM4iswUcRkdV94IVWhURZ4ACRnIGgjjcGmmRqn29HKY4k2fjgobGqP3CIvSDbtHI1BJ43fIAZ9G8F6DWgcaoDIcCktTVL/+kP7IaNESgbBJz0r+siKoNwYMAlJ8RPKOcEu1qQCgtUCH7VdWq9LeqVOdWBKTs94ssjIqEMA6mlYZY8tjbHnJmYmzfrosizTWKxDmcCRanGh0nkJoE6qbxZHuwZDCMYhs/SaGOrPIpqRDTwDaVN8qDlw86qWULnFSEkf7cskIQfU0p2s0rxu/0E71Lk/0nDNl5bGlwCLYV8a8QCRtjNkGpXXF3gtkdzPR185XNGQfgEphoxvNiZicajY1bnD4+t5eRDx+KnK884SldVsumwPAO6rRL4Youu3Hina+M3LETDzgxc69Zll3XY6bA8B99cz3+7W7zybgJ3l7tk4c0qyLQ52I57OPMafDxsGd1cx3+rQv/tfGXbHSVbFzr2WsYhw9/oJt6y05nx1uLENRGTRPaT8dVLy1O3nvACmafw2xcgOhjUPzEtlYpmaEOF8gbg1Zh2UQQJGFwqBgKTxCXEJkhHJoVk9WRkW+ABgo8yasqzg5t9gjHJ3HVURGqEQjAIgx0xmQKurydEqj7rJUabQ/TOkkgVoaLcrG84MFhgMtUBVXSlMjAHCJpgIKS5H8UljrymbjoGDooAVCaWqIKyrPpca10Y3OUEkIFY3iaG8Wa4Us+mQpcAxsKKTo8TN2J2N1tHePUIU0xaEuoSanQLW+EI1SzQpqWxNDXd2zEhlLxKFLvKcqx1iyr4w+x29etdE0eOYmMhBinOg5uyTVLIXbKtHztE5YT8mW98qnSb5NvONEjqVFAEAAd1YzPx/S7u9XTMfsa3fNagt0ib4LXGFpLom4isZ8VGIFUg7tq9HmAhzVCGRjORk995plw54cNweAd9Yyr04q05l9me/1KZGWffL517CY8Y3I6wetywlUAPDRZuaprgxm/WAYH/MqtTfdETn8TPqVsQuvM3kFvL7OEQkE8MkW5vHzi4FWwfDP7con1lts190ReulHyzrqsrCCAyFPzQi1wzIYICgte2oUANwmmCblRh0uJa1HKGMIJCh5j5x1ZFQFtU3I5hcqoUDKBmEkAlUkmz7gzZ4OAoBTAI6htDk5d6l8ORD2haCBpu2cS49wlRP1BCnzLFxJlUr165zTkX2ZGl5ymAUA1uraONRJ4/19IVxkpky4SWP9Qg49PADY7EandLXEek5NU4ZZ5KAfJ6Xs3c0UdnjQMarAdF3rTHcH1QMoMdCRXdkkBRbBTg/FT4Cx2ufzyt7N95qIvCTRm2vGBgAtLpRQKJwwuarFExzZ69TK94mjvYzNkX16Kx23VzIvjisJotzyA2HbFt/rWqcwjGMXX7e05TSvoeJdtcyPiE7VTweVqYY9yfZDmvXo6VcsG5cRqEAtAHZn7D8QxudnoWD7vvk3Dqavx88f40urc/zMqMjj4Z4a5r+6MmLJVzqVBzYWm+rXzh9/IbUuh2bjF9+wbc/J5DWF7cWoxAo/7F/c/wvtykebmfy1W3FSjHe8cfm3KqHffM+x795lbQ4AH1zFnJ/FKWWAb/UoNh5uq0T2PXclus+kZBOuOVZwIKSVRhmTNiOcl8DEAEf2sJcKhEUWWiA0W7GcTNGup+PgMgH5DKuEZtk8uv5ZCpU2GCXqNIjlGFueHPIDANZj0896cxkSq6fR3hmrA1hWre72BTGp3iuHAsCy2V8ZAMjjwSlQzJ644gp5dgrLyc4AXfZFmhrmSpcOhC0u1BfC5L1SqKiXxvp1hTqz6pOlQ19LrHl+oCsgUhiQ4vClLEKaGlznYY7SmN2mulbTyEVSlFmJReSZCT6HBqSKfWUMSWgDgE73+jukdnI93n0694wNdMYXD01zl9xtpr6T2s0vvm5u2Z775sUWaCtAmupoUIQfRyqsdps41Jm+nug7zzhcuYzhpLCnFHljoKkofLtXWbdrpzjaK88uth5wIha7cMy66cbcNweA9zUyz40q6WzI/+hQHmxi8nbdHj31ihJdfFAIH3ravvvty9ocAB5pY77cKacMXn7QrxSaYXcJctx8X/jlH6emskO/+Z5t2y1L/qmS+Kct7F+fVNQe+TEv/tWI8kgbCwjl3/mRuZ99VVUGCL/yE8Zqz30MJwUzC0/tYt/3W/nkNP7NGP6bU/JTu1gEwJitrns/4f/W5+KdJyK/+V7stWeW3ms5WNmBkCgtkvQJfcXtbFOjoM7LEHKgkDk4Sq2LAoAcnM3CpldRaaMMRoJqlhTwAcBUFPIFIOfE5Fkfm1sg1DcOHIcFs0ByTmw8u+xLCqud0EX4JSGWYwtKkr6xDqpitSTKAd+SdVcAMLNQ50CkXxLrKsZy8tLYLFWfTBztFaqWLl0CQFM+CoqUqQGhqlGeHNxZmKRMygx2mpYijKdwXQl6w4fJYa5Y1dqWQDt5eGmgQ6hZnUtRVwVdJQTgx8zGNT5toFIiwaR3LLs+mQZ3VjNPE5ywp4cUZfWu2PkjmvVo+xHLul25bw4AD9Qz3+nNOP+PBpSbyxn7ln3pSQ8AzJ94wbZl37I2ZxF8pJn5j87F85+Yxr4Y3Fxjsm29OXL4F6l18eSL5ubNWWhOVBSY4H2NzOfPLjymjc/j/+1VPtnKsE63Zd114Rd/oK7Hzh/FUiJH7kE6VuWjd9YynzgmY4DxefyXJ5R/2cYCgFDZaG7eHPjxl0BRYuePxTuOO265f7mbA8AuD3pPA7rp+eSTF5R7Xkp+YzenFszMzZusm/b6vvjJ2f/9p8ix5wre88gVbA4A+yvQv2xl7n1F/uQb8g/3cqkWhrllu/Ndnwi//JPkzKRp/fVXtrkeVnYgJEujFq0loa7idkxXcVtFsQV8NAYF6yhIVUenovr6akuVRitsMEYbHOVcRSrnfSRCEeoEgOTsVE4ZoQP6aYOjXHG5ND3hiwGLKKYZSd84n1sg1LNx4MtqpMmhDppitTQ1zLnLcrzd67XxhMrGyFAPxVpBTkpTw3x5TkkVAthWjI4RnuZIMM84Kt/G9pE/khjoEHLgyalwCtCYT7HfO5QoYTkO+0Y062LP2Sw2NyTWFiBJ0Ror9oXwEcsaU3BKYwse63jDtGpD7lEWAG4uR11zGSonCRl+MaRsvWFnou+CMr+ouS6OdIMs5zJDlI531zEvTyjpJI2vdikfaWZs226OnT+WGsyWQ7PxjhPWLTcta3MA+Nhq5ukhJSXd/pnT8iNrGRaB/Ya750++pA5OK/Mh8dgvHTfft9zNAeDvNrA/G1IOjuG4DO//rfwXrWyZFQFA/oH3R88cip5+RRy6NPfT/3C96+NLEjep+Oet7GAY73o2ue1Z+ZE2JuXF4bz7z5RYZOIzD8w9/Z+FH/r7JRUG9PCFLeyn1jIDYfyzm7gDlYsnzLv9/a57P2FqaPM88uSJwUIAACAASURBVJXc+6Yk7qtnBt/Ndb2T08gXWFq3F33iCed7HmGW+fCxJFZwIOQpnHrGbE33cwGdkVHIpTSqkxGyDqccuhwI9dj0If/SpVE7IuWUAIB1eeRZLwAMRyiTMgAg+6eyWEOkQOW8AwBXVJH0jenrny3Bdk+hRc/GoaRmbmRQwVBKSNNk8V4noc/GazJN9pLDLNL4AFdUhgRSXp2O6zzMUZqZwEl7y47IRc0iTsSkqSVEmTW4uVxb/QOAg2M4WLMgQ5wOsfu0eSk5Rw3uqkZPD2ac/+kh/LYa3tK6PXbucPp67Oxh6zIfwAUG7q1j/rt7cf+fDSmb3KjCZbWs3ZHeCZs/8ivbjluXe7vPF+D9jczj7Quvz7PDioLh5nLE2J3WTTeGXvi+uh46+N0rq/65zfBX69j3HZLDEvzbRWUqCh9pZgCAzS/Mu+V+/7f+MdF/wf+tz/Hr9+hZS2ZHoQl+uo/7wG+TlT+Qyqzob9cv3IcZe777oX+MHPr57Hf+yXnPny8ptaMHGwev3s59egP7/H724dbFmzwSTO6H/tHzyJdL/u5bS2ovZMd7G5iv7GRJkSBT43rbzgNX8Jr/frGCAyEHEdqwDM50YrqKQEjpEULmvIyuvlpw6WGZKjulxwYAbEFxciEQQjXtgS/HHmFDHkVaGgB4T0XSN9YTopD8AEDyji5JP1ChN9gpVNSFhgfWUkl+4/18Wa6BcEsROkFr4824V22MdZMiBuJQV+49PAC4voRCEphPwvP8Ws+Uts2W6Duf3fiUxIFK5lcjGfsrGJ4fVSo3bY9feD19PTkxCICWe0e+v575Xn/GGP93epX76hnr5r0ZIxWzPnG019ySKz0ghY+vYZ7qktWiCwb4l/PKJ1pYALDfeE/kt79QO1XJmYlYx3HbjgPL3RwA/mY9++MB5aVxPD6PP/G68uS2hVZ73q3vjZ0/Gjn6q8jhZ+KXTjn2P3AFmwPAw63MhkLk+Z70zW7l5zez/OU7pX3PO6yb9wWf+bqpvtV80xVuDgA7PWjg3Xz73dz/3sCmTwnw5XXF/+fLJZ/+9rJYEyR4Bm6rRFSdP9ZZtKyP4krAyg2Edpq4DGO2aTJCao8QiwkAyC42r1sazSuQF0ujuJTICLGcVKJh1kEzjU1DlQ1GaPa5XMFCRjgUwaRuIai60vqK2Cms0qP6FVcmfaM9QbxKj96QWyBsLUBdczQ/o7I6ztu3TsdjSFhKVDqFDYXo0hwmVTaOW5rXRrpTg7UpJAY7TMt5AN9ajHqCWlr94UksVbcpoz2aOfX4pdPmVUvrGqdjpwdNRHE6ieWYFxeZUc2GjZJ3JF02IXH+iKlteT02ANhejDgEKdvxw1M4iWF3CTI1rANFTvScU9fDh39h23pz9o86Faud6KZy5tOnZAD4z07FzMLtVQgA+NIac8u2wI/+XQ4HZr/7RN7N92XvtevBbYYf7ePeeyi5+qfJT7YyN13m0TO2PPef/t/4hdfjXaeK/uwL2WlCWcAg+PJOdua9/IV7OI1IkOOGu4s/9aW8296XO1WcCgsHZbSCkIHfPVZuIKT3CM1WJaYtjVK4E0s1CAGg2Ix81GGZPJeyWBoFUmVDCQcYh3PJYlGVHY3M0xgIBQviL0NhTGaESmQOcXwWRewU3GZggJLUqg633bPJZlIRWxLl8GyO0942DipsqJtUHHUVJ2Vlu5lQ81QUcXwg98FIMwtrCyizna8GHZLDnbIqXQDGYv+FHMVNVAgM7CllXsh0E/31mHJ9tY2vWhXvXrTtBozjHcdzlB5OgUXwrlrmO32L5/9On3J/PYNYzrJ+dyppw3IyfvoV86a9y9pcxaNtzGdOyxgAA3zmtPxoG4MAAKG8/e+de+ZrWBKlyaHoyZfsN77zCjYHgH/bwR4cw9ueSX6hXfmfPYtpj/Mdf4pleeofHxRqVtv3vOPKNgeA60vQyP38+AN8evUPAPiSKvefft790X/MZa4qO6zLaIwa+CPGig6EpOUeY7Zp/Ajn6LIyS9RFAcBjAS/F7h7YvALVxxIAJqOUHqEc9LP5SysqWTmwc0AOLrIFHnluGjAeCgMpd5mcmeRySAdVNDuRnsNtxDtB0hvk6TGuqGJJ/bMUNrrRGZpoZ6e1cX1cO28ieUfY/MJlPeBfT/MEPzyJ+fq2+OWMJ7U54s1L0is1uKMKPZNWvcQAzwzjt1cjy7pd6Sob4kg38EIu9EcNPtzMfP2SIioAAP4E/HRQeX8jAwD23XfOH/2lWl2MnnyJ81TlWI7W4L0NDIPgz47If35UlhRQNwcAy4br+bJa3xf/v5mv/rXz7j/Lru2QBYUmOHUX94Ut7IW7ufQBYySYCj/wt+X//IzzroeubBgkBYEBh1HkM3DVWLmB0M6jeWJYBvECKApOLkZIqqyMMh9ecuCq2KIjN+ookMMLgdBLywhz4U6ooFrJI15grI5k0D9Es4ZIzkxm8Y7XoFnHz4gpqc6bHSZJhIp3hM+B5JfCZjc6RWRscyKcNDd5/N2a9eX28ADghlLmlYmMjG10Hs8mcFnbRtVKNIV41ylT8zJ4ciruqmEOjiqpx6kjUziPhzVOZF1/fbzrhBJbqFzPH3/BuvlKMrZ1BWh9IXzxggIAf39afncdow4Q8SVV5tVbAz/5sjjaE3ruf2y3ve8KNgcABsEzN3MCCyyCX97CpXeqCh54xHnXR90fe3y5JDkNrBzsLUOuXCeQDBj4/WAlB0JKRggLDIrF6ugctTQaDTO2JTJCBw8KprgepmeEE7QeoTw3nUtGCKpZEt0josQ3NmmnWkPkIASawhonhYoHAHOu2m3KMCkykJwaXlKxOh3UeZYT0zhS2iwNdmrWxdzcD9KxpxSdnskQ7Tw4im+pYCyr1ksj3enN4HjHG5blz4MUmeHGMuZ7l1U2nupSPriKAQDGlmdu2R557ZcAIIdmY+1HbNtvXe7mKr68g/1Sh3LT88kXx/HnNy+m2s57/gwnJf9/fy7/roe43EQAqHCZ4Es72C/vZLVMGIRMTetzUbMzYOAtgJUbCPVN6jM49TqK22HGSpMXy4THgqaINuFCjxDjaBIkhRJl5aCfzV+CO6Gi1qHnEVE2PTleT52UyUH/LIUWF7pIG+wctFWvlwbJdWVqOHd6AwBsdqMLs1p7tmNepbBxjTjam56XA0Ci/0LuPDwVVg52l2SI8f9iWLmjCiGTxdS4Pnb+mLooz81Ik8PL4uGl8Ggb84VzSkSCc3788oTyoVULf1B5t74ncvgX8a6Tge8/ad91+3I51ynUONCZd3AfW8OcuJNLz6uQYC588G9L//5/rzJjM2DAABiBkARjsafPy9BLozlkhABQQmsTIo5HJosSDU9G9WRlZnLMCHVdA91l4ckJ0jseAJLTY7lPELS64AKN6neSq6sODWhXMZYnB4TlaOxaOWgt0CaFhyfxtgobX1KdriuYnJnESYnPjaGYjnvrFsX4fTE45sVvq2IAwLr15vljz6nr86//2rphz5UNlO8oRgcq0Y3PJe94Qf7idjalX8oVlha859Hgr77NuUsd+99zBTun4LHA3TUMlcNjwICBa4KVGwjtPIoQ7hMAgMw2HMuhNJqDKEOxBXmpnPq8Ajnon4xSOOMAIM/N5Gg/VudA/dRAWFSGZybqaZE6F2uIFMptSMFAer8dTnjMcjRdHwRUb1XBwtiXYH1ocGMpSm/jRZNwagZfV4LMzZvSaePxrpPm5k1XMFjxzlrmmFdRteK+2qW8q45Ry8WWlm04Ho1dOCYH/ZGjv7LfcPdyd07hK7vYv93A/GQfe399xl+TefVmz6Nfcb7z48vSZDFgwMDvHis3EFKnRgGAsWSIy8wlrsR6QkWJBaaIKAKqfW5wZiKKqSyi3ANhfR70hyjrfFGFeW6sgSC8Lyhi25Yu6qawnjDewwBn/MBVNKomsSmIw5e4imVrVdxSwfw6TZ351Um8yY0cPFjW7kiVLgEgdv6ouXXZAr4AYOPg4Vb248fkE9P4K53y/992+QPPMK77Hg786N99X/ykY++7ch8gIoEA7qpmqNZIBgwY+KPAyg2EesMyjNmmLJkRzodzCSelVkqPEADY/EI56J+IQhnJBcBYDs7kKNNXbUeTUSwSnHSuuLxgfrKRiNRJ3+hyq4ubCaGygRC2cshR16yR+RcHO9mq5SlGAsB1HtQfwimT4Z8NKndWMwDAVzQAwyQGOgAg6Z+UJoeWKyGWwqNtjNuM3vGi/OWdbLqBolCz2vPof7ofesyx9wp5cgYMGHhrYGUHQlqPEFnsOLYo2aLTIwzllBFadTLCfLc8R88IlfkQEsw5CnnwDFTYKFbySDDPsnmNik+zLk2NcMucA9xahI5nSku/4cPbipFQ16pGqRQS/Re46lzdFVLgGXhnLfM/vRgAwhI8M6y8u27hM2nffaeqxB9+8Yf2HQeuWBSKZ+A7N7DjD3CpnVNg8wv5smWMuRowYOAtiZUbCB06PULGYktNjUoKJBSwEXfg3IdlJqmceqdbnpsZn4dyokeYDPhydygFgKZ8IBVBp2IwYKm2zQ5rN58aWi6te4eHecOH0yUpj3rxLg8y1a4RR3tVqTkAkMMBOTjDLoc7kcInWpivdMqBBHzxonJrBZPqm9p23KpEQtNfeiTRd96+711XsLMBAwYM5IKVGwj1S6OLKmuqBxPZ/FHmc+oRllrRZJRWGnW65eDMeBSXE665cmCaW04gbHai7qB2sTOAg64qaUobCKWJoeVKM5dYwG1G6bOjr07iPaUImSxC1apE74KQWLzzhKlpw5VJL7a40HvqmS3PJP+rS/mnrWlK+SxX9PHHbbvvKH7433PRhDNgwICBK8PKDYRWDuIyKEScSi+N6rvy5lQaLbPCBH1YpkgvI5QDPta1DB+v5nzURTAcuuawXFwjTWRS/TCWJgavoBJ4Uzl64bKJ63AEzybwugIEAJZ1u6JnF/x6YmcPX41Y/hPb2O/ewJ59B1eZ+WSATBbrhj3Lmu4xYMCAgeVi5QZCBGChUQnTeYTUSRksiaAoyETzT8pEiQVNxzAZa1mnWw74qBlhctabi318CmtcqIMQf7kYwI7KOnGsP2Nn/ySy2q4gqNxRxfz8stv4z4fw2yoZBgEAWDfsiXeelEOBpHdUHOu3XNFUZwrbixFpjWTAgAEDvwOs3EAIAA5adZSx2FIqkbrcidzCCcdAgYnCqWcsdgzIAxFS214OeJdVGm11oc6ANtaen8VVddVywKfqMqsQR3qE5dMbAGBfGRqK4K45jAH+u0d5T8NlE1Fbnm3HbYHv/8vsD57Mu+neK3DqMWDAgIE/BKzoQGjnUZiYl8nMCCmlUWU+p7qoinIbGqe1CZN5ReuZGcr6MjPCfAEKTGggjVavYLgwi9uKeL68ThztTa2Lw5eEmmVPdQIAx8AnWti/PCF/7ZIiMHBj2eILknfgfULtGkvbdVdjpmPAgAEDv1+s6EBIzQiRxZbqEarDMhoo86FcRkZVlNvQ+DwlEEbsxWuQlt4AAPLMVI5+filsKERn0qh+3UHssSCnAKbaNekMB3GgY7mi1Sl8qpURGPTvFzNc5QAAsVze/vc49r7zKs10DBgwYOD3iBUdCO0chCmlUbsSTQ3LgJNwkFGioVwUt1VU2GBsnrLut5Y0yl7tzvMhQCj3dFPFlqIM+9njPry1CAGAqaEt0XsutXNyeoJfPuFdhYmFn93Edr6Ta3EZAc+AAQNvNazoQOjgqYHQpohxUBQACCSwi2pGaCfd2emotKFRWkY4LhRXSNpAmPQvwywwhR0e9Lpv8Vcc8+EdnoVAKI32qnJx8a5TpsZ1huilAQMGDJBY2YFQoPQIASHGbFHjBzUjlOeDuSdtlTYYjVDWe9nS4uikZnFZHkkpbCtC7bM4Nf766gTeU4IAAAlmU9MG1So9euply4brl7uzAQMGDKwErOhAaNfV3Xao1VG60GgklHtGWGVHI7SM8DwqdYTJQDh+BYHQysG2IvTyuAIAPUEck6G1YCGLte++I/zyT+KdJ6XJIUvbruXubMCAAQMrASs6EFJLowCALDYlFgaAuQR20oRG2Zwzwmo7DNG8c0/IHj7oVQuwKSSnx7niihx3TsfdNcyPBjAA/KAf312zOLhialxvXrM18MMvuu57+Iq1Og0YMGDgrQ0jEFLWGetCRhgQwUUOyywnIyy3oem41iAiIsGsLHAOV3I2o00oeUc5z7K9ZwHg/nrmhXHldR/+r0vyR5sz3lPnOz5a+tj3r9i6wYABAwbe8ljRgTCP2iMEYKx2JapmhFQPpmDu+iwsgnIrGolk/Jb+MK51IM5TlfSOLK5inPSN8VcUCF0meGIru+/55EPNjDHYacCAAQPLwooOhA4eQroZYRgAAiKmZITzy8gIAaCOsM/tD+H6PMSXZOhiJ2cmWIczF+U2Kh5sYqIP8v+wkb2yHzdgwICBFYsVHQjzspVGwwAQSABJn5AjQca2jEDYkIf6Mp2SeoLQmAd8aY00MZRalCYGDG88AwYMGPjdY0UHQoeAQiK1NOrA0UhcBgxgyaTeYUkEObksV6CGPNSrDYR4lRPx5fXSeF9qURzt48vrl/cfMGDAgAEDV42VHQj1MkKbQ5kPBRLguroGoYpmJ+oJZgTCrjm82on40upkYDplAiwOdwvVq5a1swEDBgwYuHqs6ECYx0NIpKwz1jwlGg6I2EUqbkeCjN25rN/SnA9dc4v/xABdc3iNEwHDCpUN4lAXAGA5KY32XJkotgEDBgwYuBqs7EAo6AzL2PLk+XAgQeFOyJEgu5xJGQCodaCZOE6lnsNhbOeRurOpYZ0qByoOdXFFFYzFvtz/ggEDBgwYuEqs7EDI6/QIbQ4lGppNYEppNBJc1sgoADAI1jhRu3/hF52bxesLFr5lbtkWaz8KGMfaj5hbty/3/AYMGDBg4OqxogNhvk5GyNry1B5hAVEalSNzyy2NAsAmNzpzORCenMabixa2FSobkckcOfpc9PSrti03LXdbAwYMGDBw9bg2gfCpp54qKSnJy8u79957w2GtpFhHR8eBAwcKCgrsdvutt97a17cwKvnNb36zPg3Dw8PExm8ueAY4BqJJ7bqqLDMXp5EIw3OsY9mBcFvxokHEMS/eUbz4srve9YnIqz913HwfW7AMY3oDBgwYMHCtcA0C4fnz5//qr/7qhRde8Hq98/Pzn/3sZzUXBIPB++67r7e31+v1VlZWPvDAA+r63Nzcxo0bX7yMsrJl601fPfKonHqGZczWaCREC4QBZvmBcE8JOjShYIBoEk7P4J2exURTqFld8ulvO264e7l7GjBgwICBa4JrYFD37W9/+5577mlrawOAv/7rv77zzjufeOIJlGZZvnPnzp07d6pff+xjH9u2bZuiKAzDAIDD4airq7v6M1wx8gUUFHGJRVsCZWx5iVCwoELbDpQjQdbhWu5vqXGgAhM67sPj83hbMXIY8tcGDBgw8AeDa5AR9vT0tLa2ql+3trbOzs76/X69i3/9619v27ZNjYIA8Oyzz5aXl2/ZsuWb3/xmll+BMY5EIoHLiMfjV39sFfkCBKkMCrtTDIcLaKXRK8gIAeD9Tcy/XFC+eFH5QNOK7ssaMGDAwB8acsoIx8fHf/7zn5PrDzzwQEFBQSAQsNsX5v4dDgcAzM7Out1u8vojR4488cQThw4dUv9500037d27t6ys7Pjx4x/4wAcsFkuqaqrBpUuXHnzwQY5bOG1paenx48f1ThuLxSRJSsXa7LCzwlRQDJtlzTo225Sw36LEw+GMbyVD/hgSEkQfdEk8WIV+M8yVmtGBotjyf/qPA5EIzYPYwJuMRCKBEBIEYsTZwJsM4wP/e4Esy6IoJpPEcIcOzGYzzy9RhcspEMZisdSESzpEUQSAwsLCUGhBVToYDAJAUVERefHJkyfvueeeH/3oR+vWrVNXUl+8/e1v/4u/+Iuf/OQneoFw9erVDz/88P79+3M5LcuyZrM5x0BYYJFFlnc4tBcnnYV8IFLutDgcaVVTjEORYF5pBWKXXVJ2ALx8h/olkWa+haA+CRn4XUIQBCMQ/r5gfOB/91ADocVyhf4EVOR0Q29oaPi3f/s3ve82Nze3t7erX7e3t3s8HqdTWzxsb29/+9vf/vWvf/2WW26hbsIwDMYUSt+bDacAc7TSKOtwcmNzmtKoEg0js/UKoqABAwYMGPiDxTXoV33wgx985plnXnvttdnZ2c997nMf+tCH1EmZv/mbv3nmmWcAoKura9++fR/60IfKy8tPnz59+vRpNav9/ve/393dPTs7e/DgwS996Uv33HPP1R9mucjSI7Qm5grNGUM0cmiWzSugXG3AgAEDBv5ocQ2Sm+bm5qeeeuojH/nI3NzcnXfe+Xd/93fq+szMjFpD7+rqqqmpOXjw4MGDB9VvvfDCCwUFBadPn/7MZz4TDAarqqo+//nP/8mf/MnVH2a5cApojiYug+wuR+K8RllGDvqNQGjAgAEDbzFcmyrf/ffff//992sWv/a1r6lf3H333XffTeHJPfnkk08++eQ1OcAVwylAT5CyPm9xlclzTCarQgnNMkYgNGDAgIG3Flb6KL/TRO8RBk2uIjmgWZSDfja/8HdxLAMGDBgw8LvCSg+ELgEFEpTS6Axf4BbJQDjDOim0EAMGDBgw8MeLFR8ITRCgZYTTihkAcCKWvigH/Wy+EQgNGDBg4C2FlR4InQIEEpT1mQRErIXy3Ez6YjIwzTqN0qgBAwYMvKWw0gNhgUmnNBqHmLVQDmZoxckBH+cyPCIMGDBg4C0FIxDCLC0jnI5hyeFOzk2nVrCYwGKcsS3PldeAAQMGDPyBY6UHQhMLPAMRwolpOg7IWSQHfKmVZMDLuooBaX0qDBgwYMDAHzVWeiAEgAIT8hPV0ek45lzF6YFQ9k9xhSW/26MZMGDAgIE3HUYghEIzpTo6HQeL25P0e1MryZlJIxAaMGDAwFsPRiCEQhP4CX9DXwycJaWyfzK1kpyZ4Nxlv9OTGTBgwICBNx9GIAS3Gc3EtaVRXxwXl3rkcAAnF/qH/6+9+49p4u7jAP45uPKjtViYtLYFpOOHZg+P5hnTsT7iqIKbIktI3ELM/GMs2/40bMuT7Q+2P5YYt+nQLC4zUZMleyLJhsxQZsLArcTVjKLOuTBFKDIptBQoFGyxv+7547obuuUZEspld+/XX3dfGvpOc7l3v3fXu4jXxWpzVjwdAAAkFoqQ1qTRxP2HRmfDxDKkTElOztJFJkb5wbDnV1aXK0I+AABIJBQhrUkjb/C+GaE7wGnTGSJS6PIi7l+JKDYfiN2dZbNwjhAAQGpQhJSdxnjvP0foCdLadCIihcEUdjmJKOwaVOjz8dsJAADpQRGSNp3G7y9Cd5DTKxkiSsktDN3pJ6LQ8M2UvGJR4gEAQEKhCEmXzjxwaHQsQHolEVFK/mOh2ze4aOTewLXUgn+Kkw8AABIJRUjaNPLc95AJGgvEZ4RJqgxWvy54xRYa+iV1/b/EyQcAAImEIqS1SsZ9/4zQFSCDMr6srqj1NTcpSy1JaSoRwgEAQIKxYgcQnyaFQjEKREj524fhussZlfGvCOmbynX/yWez8VN6AABpwoyQiMigZMYCv08KR+5SzoLpH6vLpaRkEWIBAEDioQiJiAxKcgXiyxzRyF0udxV+KQEAIAsoQiKiHBUzcjc+I5yYp7RkUuGYMQCAPKAIiYjyVtGvc/Hl27NcvhrTQQAAuUAREhHlqZjhufiM0DnLmVCEAACygSIkIjKpmaHZeBEO+KkwQ9w4AACwclCEREQFGTTgjy/fmuGKVmNGCAAgFyhCIiKTmhkNcPeiRER909xjGhQhAIBcoAiJiBRJVKBmbsxwUY5+meb+kYkiBACQC/xKIG7TI8yVCS6JyKhk1Aqx0wAAwEpBEcaZtcxFN3c3Qv/WYToIACAjODQatyuXsd6J/XcgVrMORQgAICMowjiTmtlrSlKytCcXnwkAgIzg0OjvjptxZ20AANnB7AcAAGQNRQgAALKGIgQAAFmTYBG2tLT09/eLnUJ2wuHw0aNHxU4hRzabzW63i51Cjo4cORKNRsVOITt9fX1tbW3L+z8lWIRtbW1Xr14VO4Xs+Hy+jz/+WOwUcvTtt992d3eLnUKOjhw5EggE/vp1sKx6enra29uX939KsAgBAAAWD0UIAACyhiIEAABZYziOEzvDXyssLAwEAunp6Yt58fj4uEqlUqlUiU4FC8VisZGRkby8PLGDyI7P52MYRqPRiB1EdoaHh/Py8hgGN2VcUXNzc8FgMDs7e5Gv37dv33vvvff/X/P3KEKv1zs7O7vIF4dCIYVCga1z5d27dy81NVXsFLLDX7iYnIz7Iq00bPCiiMVi0WhUoVjsQ4L0ev1fTqL+HkUIAACQIDhHCAAAsoYiBAAAWUMRAgCArKEIAQBA1qT2PMLvv/++tbVVo9HU19cbDAax40gWx3EOh6Ozs3NqaqqkpGTfvn0pKSlE5PF4zp07J7zMYrEUFRWJF1OCBgcHu7q6hNWamhq9Xs8vnzt3zmazGY3GV155JSMjQ6SAknXy5MlYLCaslpSUmM1mv9/f3NwsDJaVlW3cuFGMdJLCcdzg4ODly5dnZmb279+/8JrPS5cunT17NiMj46WXXsrJyeEHg8HgyZMnh4aGNm/eXFdXt4SfDEhqRnj+/Pmamhqj0Tg2NrZlyxafzyd2Ism6c+fOCy+8MD09nZube+LEicrKSv4ifqfT+fbbbzt/Mzc3J3ZSqXE4HIcOHRI+4fn5eX786NGjDQ0NBQUFdrt9x44dC3fZsCycCzQ0NFy/fp2IvF7vgQMHhPHp6WmxY0rBwMBAWVnZp59++tprr/n9fmG8o6Nj9+7der3e6/Vu2bJlcnKSH6+trbVarUVFRQcPHmxsbFzKW3ISUl5efvz4cX65qqqqqalJ3DwSFg6Hw+EwicQhhAAABPNJREFUvzwzM5OSknL16lWO4+x2+/r160WNJnFnzpzZuXPnA4OhUMhgMHR1dXEcF4lE8vPzz58/L0Y6Wejt7U1PT/f5fBzHDQwMZGZmip1IaqLRKMdxU1NTROR2u4Vxi8Vy7NgxfnnXrl0ffPABx3EOh0Oj0QQCAY7jrl+/rlar/X7/w76jdGaEkUjEbrdXVVXxq1VVVTabTdxIEsayLMvGj6tHIpFoNLpq1Sp+dXZ29qOPPjpx4sTQ0JB4AaVsdHT08OHDp06d8ng8/MitW7cmJyeffvppIkpOTt6+fTs2/sQ5derU888/L9zKJxQKHTt27JNPPrlx44a4wSQjKelPiikWi128eLGyspJfFfbwNptt69at/OHTkpIStVq9hKcPSacIx8fHo9GoVqvlV3U63djYmLiRZOLAgQO1tbWFhYVElJqaWlZWNjU11d3dvXHjxq+//lrsdFKjVqs3bdo0PT3d1ta2YcOGy5cvE5Hb7c7KyhJuLqPT6UZHR0WNKVnz8/PNzc319fX8KsuyFRUVXq+3t7d38+bNn3/+ubjxJGxiYiIcDgt7eK1Wy+/h3W73wtutabXaJWz80rlYhr/jTiQS4VfD4TB/+QYkVGNj47Vr17777jt+9fHHH29paeGXm5qa3nrrrd27d4sWToqqq6urq6v55YaGhnfffddqtbIsK2z5RBQOh3HrrwRpaWnJzMzctm0bv7pu3Tqr1cov79ix44033njxxRfFSydlD+zhI5EIv4dnWXbh45GXtueXzoxwzZo1qampLpeLX3W5XMLVdJAgBw8ePHv2bGdnZ1ZW1h//ajabnU7nyqeSD+ETNhgMPp8vGAzy49j4E+f06dMvv/zyn16XaDabPR4PLhBLkMzMTKVS+cc9vNFoFAZjsdjY2NgSfi8gnSJkGKa6uvqLL74gokgk0tra+txzz4kdSsqampo+++yzzs5O4WAFEQn7YiJqa2srKSkRI5qUCZeJchxntVr5T7iwsLC4uLi1tZWIZmZmvvnmG2z8iTA0NNTd3b1//35h5IEN3mQyCSfLYdnt2bPnyy+/JKJoNCrs4aurqy9dusQfJr1w4UJaWtoTTzzxsP9ZOodGieidd96prKwcHBwcHh5WqVR79+4VO5Fk3bx58/XXXzeZTDU1NfzI4cOHKyoq3nzzzR9//PHRRx/lryb/6quvxM0pPXV1dX6/32g0/vzzz36/v6Ojgx8/dOhQfX19R0eHw+F49tlnS0tLxc0pSadPn37mmWdyc3OFkffff7+9vX3Dhg0jIyM//fTTmTNnRIwnJU899RT/nW/nzp0KhaK3t5eIGhsbt2/ffvv27ZGREZZl6+rqiMhkMr366qvl5eXbtm1rb2//8MMPhev4Fk9qT5+YmJi4cOGCRqOxWCyLf04HPKxgMNjX17dwpKCgQKPRzM3N9fT0eDye7OzssrIyfDtedj6fr6enZ2pqSq/Xm83mhadDhoeH7XZ7Tk7O1q1b8RiyROjv78/IyFi7dq0wEgwGHQ6Hy+XKysp68skn8VTI5XLlypWF3SR8sZucnOzq6lq9erXFYlm48f/www9Op7O0tLS4uHgJbye1IgQAAHgo0jlHCAAAsAQoQgAAkDUUIQAAyBqKEAAAZA1FCAAAsoYiBAAAWUMRAgCArKEIAQBA1lCEAAAgayhCAACQNRQhAADI2v8A4RWpSFM8trIAAAAASUVORK5CYII=",
"image/svg+xml": [
"\n",
"\n"
],
"text/html": [
"
"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"plot(0:0.01:100, besselj0, label=\"Bessel J0\")\n",
"plot!(0:0.01:100, besselj1, label=\"Bessel J1\")"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "f71394c9",
"metadata": {},
"outputs": [],
"source": [
"# savefig(\"figbesselJ0.png\")"
]
},
{
"cell_type": "code",
"execution_count": 18,
"id": "172fdcfd",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.08825696421567697"
]
},
"execution_count": 18,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"bessely0(1)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"id": "e5606af8",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dZ2AU1d4G8P/MbHpCQkgCSUhCCwESepGOlNC7IoiIICqKCoKKYrlXQV4VG6KUK4peBQS8IB0podfQIXQILRXSe7I7c94P680NyZbZMmSWfX6fkmWZOTtns8+eM6dwjDECAABwVnx1FwAAAKA6IQgBAMCpIQgBAMCpIQgBAMCpIQgBAMCpIQgBAMCpIQgBAMCpIQgBAMCpIQgBAMCpIQgBAMCpPdQg3LVr19atW+U/X6vVKlcYkAm1oAaoBTVALaiBErXwUIPw6NGjR44ckf/8kpIS5QoDMqEW1AC1oAaoBTVQohbQNQoAAE4NQQgAAE4NQQgAAE4NQQgAAE4NQQgAAE4NQQgAAE4NQQgAAE4NQQgAAE7NsYNw/W1p+XWpuksBAAAOTFPdBbDejmT2xhHJS0Mio+ciHTvRAQCgujhwEH56RpzXgQ/35sbuEZ9txPNcdRcIAAAckKM2pFKL6FwWG16P7xjE+bnSvjRW3SUCAACH5KhBuD1Jig3lXXkiolEN+PW3cKcQAACs4ahBuD+NPR78d2do/7rcjmS0CAEAwBqOGoTH7rGOQX8HYUt/7n4xSy2q3hIBAIBDcsggLNTR7QIWU/PvIOQ56lybP5SO3lEAALCYQwbh+SzWxI/TVCh7xyDu2D30jgIAgMUcMggvZLPm/g/MlmgfyB3PQBACAIDFHDIIL+Wwpn4PBGGbAO5MJkMSAgCApRwyCK/ksijfBx6p5Ua+rlxiHqIQAAAs45BBeC2XIn0rLyTTwp/OZSEIAQDAMo4XhCKj2wWsgU/VIOQSsqulRAAA4MAcLwiTC1ktd85dqPx4dE3uQjZahAAAYBnHC8JbBVTfx8Dj0TW5izkIQgAAsIzjBeHtAhbhbWCniShf7noe02FWPQAAWMLxgvBOAUV4G3jcXaBgD+5mPhqFAABgAccLwrsFLMzL8N6DTfzoci6CEAAALOCAQVjIwgx1jRJRY1/uSu5DLg4AADg2xwvC5EKq62X4nxr7ctfQIgQAAEs4YBAWsRBPoy3CqwhCAACwhIMFYalI+VoKcDf8r5E16Frewy0QAAA4OAcLwrRiVseDM9weJKrrxWWVskLdQy0SAAA4NAcLwtQiCvY0+q88R/W97bz09u4U9tlZKaUIPa4AAI8mBwtCfYvQxBMa+XLX7ReEv12XJu4XbxewjhvE1CJ7HRUAAFTE0YKwiOoYbxESUUMfupFvn3NllNCbR8Wt/YTFXYSJjblpR0T7HBcAANTE0YKwmNX2MPWEhjW4G3ZqES64ID5Rn4+uyRHRuy2FQ+nsPLZ5AgB45DhYEN4rpiB3U12j9gpCnUQ/XpGmRf99fTw0NKUZ//1FrGQKAPCocbAgTC8m0y3CBj6UaI+u0bgUFuHNNfH7X+hObMz956ZUgv5RAIBHi4MF4b0SFmhysEw9Hy6p0A57UKy9JT3V4IGLE+LJtQngtt21f6Pwai47k4lOVwCA6uFoQVhMQUZm0+u58lTbg7tbaFOuMKItd9jQ8MqJO7Ie/+ctOyfWO/Fizy3ik3Hik3GiFj2vAAAPnYMFYUYJCzB5j5CIGvjQTdt6R89lMS8Xalij8omGhHPbkiTRflH4y1VpWxK78KTm0pOaUpHNOo6OVwCAh82RglAnUb6W/N3MPK2BD5do266EccksNtRA3Nb14kI8ueP37ZOEeVp697j42+OCnyu58PRLD83y61JCtiJ9pMU6ktD5CgBgiCMFYWYp+bkSb6ZBSPV9bN2ed0+q1DPY8Gn6hnI7k+0TKYsuSn1D+Zb+f5+olhvNbCHMPmXn7tEiHY3fKwYs14as1P5xE32vAACVOVIQZpSwQHP9okTUoIZNA0clRofSWbc6hq9M71B+d4od4kQn0aKL0lstHjjL5Kb8nlTplm0pXhEjGr1bx4juj3PZ2k8z/ai06Q6yEADgAY4UhJmlVMvkSBk9G1uEF3NYoDtnbJJG19rcyQxm+ySK7cmsrhe18H8g1700NK4R/9NVu2XVj5eljBL6ubvgqaE2Adx/egsvHRAzS+11+L+dzGDt1uu4H7WdNurOYc0BAHA0jhSEGSWslpv5FqGNQXg4nXWubfQs3i4UU5M7es/Wj/sV16XxkQYu/vON+V+vMbvczyvU0T9Oiou7CJr/nqdjEPdkfX7OaXsOyTmbxQZu173VnNc+7/JKUz52m+5ijv2zMLWI5pyWZhwVj9h85QEAKnGkIMwskdUirO1BBVoq0Fp5lmP3WMcgU3HbrQ53IM2mj+NiHW27Kz1Z38DFb+7P1XSlQ+l2+Lj/4bLUPZhvVeuB1/JBa+G3a1Jase2HJyIqEWl0nDi/ozCmIa/haXwk/3l7YXScWGrX0a9H77E2f2rvl7DaHtzoOPGT04r07iZks7lnpEUXpSx7t5gBQOUcKQizSqmWuSGjRMQR1fPmbhVYmSXx91mHQFNB2LUOdzDNps/iHclS2wDO2PbCoxvyqxNt/awXGS24IL3donL91vagsQ35hRftk1RfnZdi/LmnG/7vLBMa8w1qcPMT7JZVdwvZiJ26n7prFnQS3mnJHx+uWX5d+tl+vcd6n52VYrfqcsvYkXssZq32oG1fdAzK09K8c9LYPeL/nZGykbUAauJIQZhZyvxldI0SUX0fsq53tEBLN/NZTE1TZ+lcmz92n9kym3DDbTYswuiVf6Iet/42s/GT+K8kVtuD2gUYeCGvR/NLL0tlNkdJdinNTxDndaj8Qr7swH95XrTXZ/2k/eK0GGFg2N8vpLYH/RkrvBMv2nFI0cKL0m/XpFMjXOZ1EH57XPh3D80TcTr7zmO5nMNardNdyGYDw7ib+azlOt1Ze99M1Un06Vmp4Wpd4HLtqDjRXkvPV5JSxFZclzbclrD9NTxKHCkIs+QNliGi+j5cYp41pzidyZr7cy4mr0otNwr25Kz+oGRE2+5Kg6ssW1OusS/n60I2zlb8+ar0QpThl9HYl2vmx228bWsSLrwoDQ3nG/hUfiGRvtywCP47eyxQvu6WlF5MbzV/4IU09ePeiBHePGafRuG5LDb7tLi5n1C+4XNsKPdFB2HMbrt18KYXU/+/xH+05v/dQxjXiF/aTfiqIz/wL/GOtZ0WVZWINHC7bn+qtKGvkPCEy2NBXJdNOvveT2VEc89ILdbqNt5hiy5KjVZrN9+xf9buTWWDtusartb13qpTaIRzZin9clX6JkGKt9OEYIMKtFSM7wqOw8GCsKarrGfW9+FuWvUpczKDtTXUiqqkUxB3xNrbeKczmJ8bV79KflQ0JILbbMOnQHYpxSVLowzdg9R7Por/xbbexTKJFl0S36zS9ar3dgt+0UXRxrG1EqMPT0jzOvxvsE+5Gc35ExnsmM0f9Izo5YPi/7UTKlXH+Eg+ypf7/JwdPogZ0YR9umcjuQmN//cyRtXnpzfnx+4R7bVK0cT9Yi13bnM/TUxNrrYHvdWc/6WHZuROnY0rS1Q09bC4+Y50/gmX1b2E7QM062I1kw+Kq27YM6vmnpHG7xWfbsjvGCC82ox/+5g07Ygd13EiIvr5qtTkD+1fSexWPntmjzh0h/0HUf+VxNqu19VeoQ1Yrh24XXdBgSUyzmWxKYfE2G26CfvE/Qp04xPRveK/vy4cSLO1d8oYkdHFHGbHfdRt4UhBmF3K/GXMIySi+j50y6qphKcyWBt5QWj1wNEdyayvoWVrKhocxm+5a/3748/bUp9Q3tf4l4aR9fjD99g9G4bMrLspNfPjmvkZfiFRvlzbAM7GT8m1tyQ/N+pX18Ap3AV6vxX/sc3DX9ckSlqJJjY28FcwvxO/IEFMtm3RWiJacV26V0L/bC1UevzN5rwLT4vs0W5edlW6nMN+7i4IFS5V/7rcrFbC2D2i7QvQE9Gii9L+NLZ9gKa83dwpiNs+QJh2VDyRYZ8PsvkJ0u83pPjhmnGN+IY1uJH1+PjhmpMZbOYxu428+vys9NlZ6cBgzapewredhItPahr7Uo/NuowSe52B5p6RXj4ozm4r5D/nkvWsy5BwvucWne29L5VO0W+bLsKbm9lCaB/ITdgnvnLInssUM6JPz0rN/qPdnsxuF7Aph8TOG3XXcu0ZVyKjr89LoSu1I3aKvbeKkWt0a2weFWEr9hDNmTPngw8+kP/8vLy8ir+2XKs9nSHJ+Y9nMqXma7WWFY4xxljMf7SnZJziXKYUtcaa4zPGem3Rbr5j5hRakdX6tSyl0LozsH7btGsSRdPPGbdHt/CCmefoVaoFvZ6bzZxi8x2pw3orL5Fe+/XaDbeMnqJEx0JXas9myno/GKSTWNQa7a5ko0eYeUz38kGd1cdnjBVoWehK7dF0w6e4lC0F/lZ2v1jWoQzWAmMsrYgFLS87n2XgFBJjsVu1X56TVcsm6Mt5LdfAKf5IFBuv0RbaVM+MMbY/VQpeUXanoPIpskpY5Brtqhu2vgTG2JpEsd4qbXJh5VPMitd12agtk3cGY7Wg9815sdkf2tSiBx48cV+qvbwszvjbzCJvHdW1/fOBU+SXsUF/aUfu1OnscQadxMbt0XXdpC2vC4mxRRfF2svL4u/Z5yUUatmAv7SPb9Zeyv77gPtTpaZ/aF87pBPlncF0LVjHoVqEZVRTxqhRsnYqYbGOEvNZtMmRMnrNanKpRcyKcfbFOjp+n/WoY+YUGp76hPLbk6z5lpRdSkfvsYFhZmp2TEN+lbXfwm7ms4RsU+N9iKh/XS6tmKyeX384neWU0eBwo6dwE+j1Zvw3NgxP/c9NKdCDeocYrYuZLYU/EiVbdjJZcEHqWpt7zMhsnCZ+3JiG/CdnbGrxvHdcnBDJGxzexREt6iJ8dlZMLbL++IzopYPiR22ERlXWoCeiJ+vz7QK4j0/Z9BIKtPTcPvHHbpowr8qnqOlGq3oKU4+IKUU2tUgS89mrh8Q/+wghnpVPMbe94O/GvX/C1nbnnlQ275y4rb9Q58G1ONoGcKt7acbu0d22+ZbwkkvStiS2c4Cm4im8XWhdrCZfy962R9P5jSNiahHb0f9/dcERvdKUX9pNGLrDDu1CrUTDd+qC3LmdAzTlu712q8MdHaY5l8VeP1Jtuw7YJwhLSkr++uuvrVu3FhXZ8DdnTnYpqylv1GgNF3LjydKuvws5LLIG5yrjkggctQvkrLjZfvgea+HPebuYf2a/utxfSda87TbflXqF8F4aM0+LDeUuZDPrPl9+u86ebsibvlACRxMiuWXW3olcdEma0pQ3va7si034Dbclqzt4Pz8rvduyco9lRbXc6Pko/uvzVr6EfC3NTxBntzV1md5vJSy/JiVZm7Xns9jWu9L7VfpdyzWqwT3fmP/HSes/X36/IZWI9HJTo6/im47CL9ekSzasovDRKbF7Ha58YHAlbQK4l5vwrx+2/hsPI5q0X3y3pVBpTq0eR7Ssu7DiOrNl8m52KT23V/x3D024t4FT9Ajm3mohPLtXtGWhjFMZ7J+nxPWxQtXGgCtPa3prNt5ha21bTPjnq1JcClsXq/Go8ukxJJyf3VYYvlO0cbTw64dFDw33U/fKN/5ruNDmfppDaexb+828sogdgjArK6tNmzbz5s1bsGBB8+bN09LSbD9mVTqJikXykREheg1qWNwoPJvJDP6pGPRYIGfFYI3dKVIv402QivqGcruSrdnyaf0tNjzC/ClceRoSzq+9ac2f5vLr0rhG5t85Exrzv9+QrLh7kVFCW+9Kzxlaeacifzd6op6VK9LtSWWlIhn78C33Rgz/6zUrp/3pF1Vv7GvqFLU9aFIU//lZK//4/3FSeqelUMPkH8WsVsLGO9IVq77Ll4g067j0TUfBxDeSIA+a1VJ4y9rmyMUc9tt16YvHTH0jea+VcD6bWfe9kIh+uSoVi/RGjNG3U4A7fdeZf+mAaPWcohnHxOH1OINb1vz9hBhe4OjbC1aeoESkZ/eK8zsabpcTkZ8r/d5TePWw9a3/K7nsnXjxP32Mvp1ebMJ3COKmHrb+S9WK69K+NLb88QduZpfzcaENfYVPz1bP6lF2CMLFixfXq1cvLi7ur7/+ateu3bfffmv7MavKKSNfF5IbU1b1jp7NYi3lB2EQd+y+xW/rvans8WBZ1zzUiwv25E5aOBKhRKS4FGmQ8R7Fip6oz627ZfFLiL/PeKL2Jtcc0KvvwzX147betfgUv16ThoTzcrrBX2nK/+uSNV8Xvk2QpsWY3ciEQjy5oRH8D5ctfgklIn17QXy3pfmKeLO5sPKGNWv9nMpgJzKYibaanp8rzYgRPrJqV5PvL0ptA7guxlcc1JvSjL+SS7tTrPn8evOo+H4rIdDktCg3gb7pKEw/as2QkJwyeu+4uLiLqSwnopH1+AY16BurWv97UtmeFPZ/7UxlOc/R0m7C/50RresgnXNabFbzgZUrqmofyE1uwr9yyJqgEhlN2Cd+1EYwNvxN77tOwv40Zt3Yn1v5bPpRcVUvwURjJsKb+6Gr8MweMc/adcGsZocg3LBhw5gxYziOI6IxY8asX7/e9mNWlVPG/OT1i+o18KEbFg4cPZvJWvrLD0I+/p5lA4sLdXQui3UyuX5bRX3rcjss/Ba8O4W19OfkrL9DRH1D+bNZFo8dXXVDMv0HWdH4SP7Xaxb/5f90xegkyEraBHC1PWi7hVfpZj47lC6rUUtEb8Tw31+0uF376zWpXQAv535zkAc904ifn2Dx59fs09LMFry7qY/fv70Wze9NkSwdx5+npS/PiZ+0M3+VXHn6pC0/67jFX0h2JrPEfHrFXJYT0aAwLsKbllyy+CP441Pi8Hp8axlfcL/tJHx53uJxwlqJXjskftuJN3u/o1ENbnqMMPWIxS8hIZv9eEX6rrP5mn6/tXA9j1mx29r8BMlLQ680M1MR3i70c3dhymGLlyFkRM/vF2e2FMx+wA6N4PuGcjOOPuybheZuJcmQlJQUFham/zksLCwpKcnYM3Nzcy9evLhkyRL9r66urk899ZSHh5GNHohEURTFv69IVgn5urDyX82q58WO3WcWTdQ6n82ifSWZZwh0JW8X7kqWLtJk31dFB1OptT+5cnJfRJ9g+r+z4qwWFryEjbfZoDCSeQINUb9Q7s+buheiTL2EirUgMVqdyHb1l3thR4TTm0el+0Xmt1Mud/Q+iYw6BcqtiBejaMklsZ+8Dme9RRfZ+Eacm7yKiPGlRj7sjxu60Q3knkJi9PV5tqSL3Jp+oxnrsJHNjCETM14q1gIRJWTT8ftseXdZde3O0fQYbvYpceXjFlylr8+xvqFclI+siniyHn1+lq1N1I2Q0S2vJzGaeUz6pA3Hy/uz/rwd12+7+EwDZuIqVXIll624zs6P4OWcIMKTJjfh3okX/93d6EuoVAtEND+BRXizwXVlvYbp0dRmA9twSzfYXJ98OUY0+YD0cWsu0NX8GQSiJZ25MXulXnWYn+yrlJjPPjvLDg/m5NR0p0AaGUHTj+iWdbPgvbT4EisR2bSmsq7S5+2ozQa2+bY4oK7hJ1StBdN4nte300ywQxDqdDpB+PvbiouLi1arZYwZPHFOTk5KSsrx48f1v3Ic16tXr6CgIGNHLisrKy39+7tHRiFfQ8OVlsptM4d58CtyuVLZq4PcLSQ3XuPDlZ/QvHa1NIfTpHB3ud+/difxXQJJ/kvo4EdnMjUZhWU+Glmpw4g233HZ0ltXWio3OweHcL/c4J+tZ+oqVayFA/f4IHc+wl0n8yq5EcWGaH6/pn0hUu5V+vGSMK4+lcmuuOEh9E68JjFbF+pp/slEVKyjX65q9vaT+xKIaEpj7osLwvBQueMENidxPhqhg5/cU9RxoX4hwsIL0pvNjL7qirVARJ+c1rzehPGi3OVvJtanr89rztwTm/rKem/klHELLwp7+lpwlf7Zgp91gu9XW2fwDlBVq2/xrjw/MFgr8xSRntQ/RJh7WpzTSu576e1jmhnNmA8n9xTTGlObzS77k7WPBRi+SpVqIa2Y+/ycsDvWgqv0VVv+taNCV/+yqgNSDPr3DV6U+HERcl9CG18aECK8G0/ftpf1zmBELx8UZjSjUFdR5ik+jKEOW1223NL1CZZVEbcL6OPTml19RW2ZrPeeC9Gix/hJh4T4gTo/VwP/pVItmOXu7m42CO0wj7B58+Zr167V/7xjx46IiAhjz7RlHuEfieITuyyY1HUzTwr/3YL5TZtuS/22WTYf6qtz4muHLChS102mZq0Z1Her9k/jc+kqOZUhRVo4u7FAy3z/XZZVYuo5FWth8gHdvLOWTerackfqvFFuqQq0rKblEyhfO6T750m5FbHsijjoL8uukiixRqu1h41MB6yq6ybz8zgrSciSgleUFRkvV8VauJwjBS0vyy+z6Axs3lnxqTi5V+m947oXD1g8h7L7Ju2yK7JeeImO1V+l3Zdq2Z9DcqFU69eyW/my/tfOJKnRam2phVMQl18T2/2pNTahrdIMtnF7dLPiLb5Ko+N0H5yQ9b/uFbPay8vOWDhZNqeUha7UHkyT9b9+viK2/VOrtfAq7UiSIn7X5sp4B0qM9dqitfRDgzE29bDumT2Gr5JK5xF27do1Li5O/3NcXFy3bt1sP2ZVuWXkK3vIKBGFeXP3ipn85SLPZ7MWsm8Q6nW0ZH2ZYh2dybTgBqFev7q8/Btgm++wIcaXMDXIS0O9Q/iN8pZz00q07pYkv4dQr28odzOfXZU3anFNotStDh8sr21XbnJT/scrTOYSKt9flF6LlnFjrQKeo2kxcudRHLnHUgppZD3L/rKia3KPBcodATv3tDQtWpAzCaeiKc34A2mSnJmd6cX0r0vSB60s/nD4rIPw0SlJzhqb31+UYmpy3c1NqK0kxJObFiPMjDd/lXQSvXFU/PIxM5N8qhrbiPfQ0NIr5k+xN5UdSGMm5q4Y801H4V+XJDnbdk4/Kj7biJc/cEHP15W+68S/cMD8GocpReyd42LVyQxmxYZy/erKupP33QWpRKQZzS1+L33aXjiZwey7hp8JdgjCqVOnLl++/JNPPpk3b97ixYtnzJhh+zGrytNSDdm93kQkcBTmbcHA0fNZrLmFb7g2AdylHFYkr8PsyD3W0p/ztLArul9dTn4QbrojDZE3XrSiJ+tzf8ibWb8zmTX25QzOlDJBw9MzDeUubfrjFWlSY8uOT0QxNbmGPrRexmC2g2msQGt42TbTJjbm96VKcrZ0+Pys9FYLXmb3YEXvt+bnnZXMfnW7msu2J0uvRVtc0V4aerelIGfm+JzT4vhI3tKKJqJOQVyHQM7sKgf3S2jeOQP7lsjxVnM+/j7bk2qmIhZckMK8yPSaDwZxRAs7C/84KZoex1si0ssHxQWdzE/YrSrYk+a0EybtN3OffeNt6dg99nFbi4OWiEbU41vV4mYdN1XX+gEsrzYzP4DFoK8eE/almpm5eC6LfXJa/LWH4fkSpnlqaMXjwrSj4sNZjNQOQdikSZPDhw/n5eXdu3dv3759rVu3tv2YVeWWWXCTXK+hJQNHz2VZ3CJ0FyjGn5O51uK+VKlHsMVvh+ianMjosowvjylF7EYe62pupHtVQ8L5g+lMzlS5FdelsbLHi1Y0MYr/9br5FltCNrtdQGbXxDHotWj+OxmTtOZfkKZGm581UZWXhl5uyn9hbhnu81ks/r40wdwMSIPaBXCtanFmp2r846Q0PcbM3EFjJjfhL2TTXpMpcjmHrUk0NUnftHkd+G8SRNNLBMw6Lj7TkG9icqS+MR4amt+Rf+WgqebOrXz26VlRzjBLg5r7cy9E8VNMzkP48ITYshY31PKg1XupCe/tQv93xmhdpxXTy4fEn7sLln51Lrews7DuFjOxQ8i8s1KBlt6TMcPHIG8XWtlTmHJYNDZFNbuUnowT53cSGhqZ+2hWmwBOP4s/t8y6A1jAPivLREdHz5s378svv2zVqpVdDlhVbhn5ulp2QRvV4K7L65ErFSkxnzWRPf6znPxtKOTPIKxkYBi3VcYC3BtvswFhvKVdHETk7UJ9Q3mzEwoLtLT1rvRUA2teQjM/LsKbzE4o/NclaVJjzoqXQEQjIvjbBWT6S8n1PLY/VZpgaIltOaZGC/+5aWbFtdmnpbeaCzLHQVQ1px3/6Vkx3/hoquP32cF09rrlzUE9N4E+78BPP2pqJe5pR8T3WgkyZ+BUVd+He62Z8JrxhWD2prIdSewjqxo6esMi+Bb+3IdGmrYSo4n7xZktjM49l+MfrYWb+exfRr6UbE9iqxLZQmuDlog4on/3EBZfEncmGxoMItHoON3LTYWuFnYdV+TvRqt6CZMOGN7+YtMdacEFaVUviztFK2ofyH3WXhi03cACeIU6GrZTNyycs+6rc7nJTfg+odyg7TqlZxY6zFqjeWVk6bfgRjU4mc3qSzmsoQ/nZvkbu0tt7rCM24RFOjqVyTpb3lwjokFhvJwtmdbfluSPXK/k6YbciutmTvGfm1KPYD5A3n6QVU1uwi82OQksT0srb0gvNbHyDanh6Y0Yfp7JJVo+PytNaWZNX5ZegDu91ISfc9roKU5ksCP3mJxZcca09Of6hvJzjaw+KjGadkT8pK31L4GIRtXnA91pvpHW8/LrUnoxvWZuPplps1rxifnMYGd4ThlN3C8u6Wpli7bc4q7C6kRmcMPCf5wUeY7etPymVEVuAq3uJfzzpLivSuv5Yg57bp9uZU/B6r8FvRBPbnUvzbi9ukorNWolenavGODOWXGPtpJOQdz8jkK/v8RKu0BvuC29cEBcHyvUrbK4q6UmNuZfasJ33yxWvPecVMj6bNVF1uA+72D9d4Vy33QUWtfiumzUnbLTPicGOUwQ5mpNzbIyqJEvd01eEFrRL6rXpTZ/KM38CoKH01mrWpx1n1+9Q7iTGWa6LnPK6Og91t+qTkUiGhjGJ2Qz02te/HxVmmD53btyTzXgT2cyEwt9/XRF6leXD7XhL/PFJvz+NKPTxhPz2frb0lQLhyyvKDMAACAASURBVMlU8nYLYeNt6byh8SaMaPpR8eM2vNXNQb3POgg/XzU8pGXxJYnn6DlrW7Tl/tVVmHdWPJ1Z+RQ38tibx8Sfe9jUSiAiV55+7ym8c1w8XuUjfsxu3fAIo8uKylfLjf7TR5i0v/JufN9flH6/wVb10ljT/f2gxr7c7700o+J0FXsX4zP52K3i1x2Fbja01cp1q8Mt664ZvF239LKkb6NfyWWx23QlIq3saWYpHJmebsjP78jHbtX946R4JpMdTmcvHBCnHpG29NPIWRxKjpkt+Nlt+T5bdeP2il+ck6YcElut0w2L4H/sbp+XwBF911l4uwU/eIeuwwbdkB26NlvdrF4MzxiHCcK8MlbDwq7RyBp0LVfWM89bG4TBnuTnxpldcXh3itTL8huEeh4a6hnCbzHZr7jxttQz2PqGgptAoxvwv1w1+iou57BreWyQtUFLRO4CvdzU6MDLMonmJ0hv2fYt3ktDM1sI750wfIr3T0hTowX58/oNqulGH7cVXj5kYPXkHy9LZaLhrQ0tUseDvuggPL1bLHiwL+h8DvfxKXFZdyuGHVRW34db1EUYsVO8VWEoWXoxDd4hftxGkLMIi1nRNbmfummG7NDF/XfdtexSGrFT5y5wX9ijlUBEHQL/DqpPTku38tm5LDZxv/j9BSluoJkF2+TrGcxt7KuZekQcuF338SnpyTjxmYMuP3QTbOzuq2hQGLd7kOb3G1LQCm29Vbrum3VDwvk/+whW9E4Z82R9/vhwTVYpPbtXnHZErOtFZ0dq2snYdVW+sQ35y6NcutTm0opZwxrcmZGad1vaJQT/Z3wkf3uMy9ePCS9G8RsfL7V0MLB5dp+QYYIt8wgf26A9Insil16ZyNyXlcmZSNR3q3aLuT0CjZm4T7f4oplzdFhv8ZSpiv59VRy+w9TEo4F/aX+/btOebWczpbCVhqcT5eXlvX5Y7swnEzJKmLFJYAsviAMtnNtnUImONVqt3VqlKuOSpYjf7bBtHmNMlFivLdpK0xbPZ0mBv5UlGNoU0DqTD+h6b9Hm/Xee1oVsKXRFmaVzE01beEEMWaFdfk28WyCtvSnWW6Wdc8qex2eM7U6Rwn/Xdt+kHR2nC/itbPoRncxt/+S7mSdN2Keru1LbeI32/eO6PAvnVspRrGMrr4sfnNAtuyImZ+Xb/wSMMcbuFbPEPMnS+XzOSYl5hHZYWebhyCuzbPoEEbnwFObNJeYxs+PTzmaxFv5WFqxHMLftLnu5qdEnZJXS5RzW0cIZhBUNi+CnHdHmaQ3fWblfQofT2ZreNn1HauHP1fehdbcMDIfJKeNWXJfOP2HrW6WWG01pxn9wQvrt8Qe+7uaW0SdnxC397PBWdBPoX12FCfvEE8M1Qf9duS+jhJ7fLy7uYv0AvIp4jlb21HTaqPNxkWY05zmiM5ls6A7x206CnJVFZVrYRZh6RGyxTvdsIy6rlFYnSp+10o6qb9uNtQdNacY39+c+OS2+Hc+ifLklXQQrZpWY1jOYuzpKsyeVZZaweR0Mb1Fko3o+3M/d7dd6MsRdoPLFdfMt3+VUpkB3CnS3//UBmRwmCPO1Fg+WIaIoX7qSayYI04tJJ5HV9417BnPvxIuMjPZZxaVI3YNlbXNojK8r9Qzh196UDPa8rbwuDQm3aQCF3lst+I9OSaMaVO7TWHxNGB7BV93R1ArvtBSi/6Pbmcwqbljz1jFxWISsZZHl6BXCTYriB2zXbeorhHhyqUU0dIduXCNugM33pcrV9qC9g4Sndos/XZFCPOlcFvu+s2DdeFpjBI4WdhaO3mPb7krh3typERo/qcSOx9frVofbPkDZTwA3gfrX5ciCbWMAqoHDBGFeGfNxsfjPKcqXu2LuNuGZTAt2X6oq3Jur4cqdyzK6c8W2u6x/XVs/JcdHcl+dNxyEMlemN2twOP/xKemPxAcahWnF9MM14fhw+3zKe2no5x7CM3t0ewZponw5Ilp0Udqfxk4Mt+f78J9teFeBmq/VxdTkErLZjObCezYPwKsk3Js7MlRz/D7LLKGudTj522RapGMQ1zHo75rNt3ArFQCQzzGCkBEV6izYlbdcEz/usLl5fmezLNiP16C+odyOJMNBKDHadlf6sLWt13lwGP/aYSkhm8U82P+2N5VJjKyYql8VR/RNR+GZPWKfUL58UMlrh8UJDcT6Phb2ShvXM5j7vIPQfbPu6Yb8nQI6n8W2DzC1RZl1ZrXkX4jiE7JZtB8XZHR3E5twRB3sNO4OAKqXY4waLdSSh4asGIfU1I8zuyzLGUu2ITRoQBhvbLZ4/H0W4M7V97H1E1PD05Sm/JdVVjb5/Kw4vbndBmh1q8ONbsCN3KXLKSOJ0azj4u189m603F0XZHq2EX9wsKaOBzcojDszUtPA5otjUKA79QxWKgUB4FHiGEGYr7WmOUhETf3Mz204ncna2DaYuFcwdzqTGdysct0taUQ9+3zQv9qM33r3gbV696SyK7k03qoFvYz5vIPQphbXcLW23ird0Xtsa3+NHUdyl4v05d5tyU+KssOtTQAAGzlGEOZprblBSET+buQukIldpwu0dLfAmsXVKvLQUJ9QfkOVRZ8Z0ZpENqq+fS6yryt91EZ4Yb+on0yaW0aTD4rzO1q8vr5pPEdfdxQuPOkSN1DYM0hjrylZAACq5RhBWKAlSzedKRfjzyVkG/3XM5ksxt/K9S0rGtOA+73KjiH7UpmvK1m6qYUJrzTjQ7y4Idt1K65LfbbqBoVZv+yvaXU8KNK2LwcAAI7CYYLQ6vEUMTW580aW3SKiExnMLossDA7nT2eyWw9OM/rhsjQpyp5XmCNa1VPoW5ffeIe90oz/uqOyM6gAAJyBY9yiydeSt7UlbeHP7UkxFYR9QuwQhO4CTWzML7gglYfTrXy2M1la0tXOAyI1vK0LCgMAQEWO8ZFaoGU+Fi40Wq6lP3fG+K7cx+6xDjas+VLR9Bjh12tSeaPwvRPSq814G1fZBwAApTlIEOqsbxFG1+Ru5LFiQ1MAMkooo8TWkTLlgj3p3VbC03vE+yX0bYJ0KoPNbIGuSwAAtXOMrtECLXlZ27RyE6iJH3cuiz1WpeV35J7UIZCz4zLpbzbnM0pY+O/aFv7c1v72WdwSAAAU5Rgf1QU23CMkovYBXPx9A0F4MI11rWPnwSyftRc+a4+GIACAw3CUrlHmbdU8Qr2OQdxRQ/vI709j3e2xwSYAADguxwjCQhvmERJRl9rcwSorjuZp6UK2gWYiAAA4FccIwgId2bIWV6Qvp5Po5oOT/PakSB2DOHf0YgIAODfHCEIbW4RE1CuE25n8QBBuvcsGhDnGywcAAOU4RhIU6piXxqY+zIFh3OY7/wtCkdHG29LQcPSLAgA4O0cJQpu6RoloYBi/P03KLfv717gUFubNNayBIAQAcHYOEoQ2zCPU83WlfnX5367/vS724ovS84Z2ewcAAGfjGGFge4uQiGbE8F+ckwq0dOQeO57B7LuNHwAAOCjHmFBfqCPbV2l5LIgbFMY9vkWXUsQWd8GyLwAAQOQoQVhk82AZve87C5vuSGFenI1b0gMAwCPDMYLQLi1CIuI5GqbMTrYAAOCgHCAVJEalInk4RmQDAICDcYAgLBbJQyB0ZQIAgBIcIAiLdGgOAgCAUhwiCJmnPUbKAAAAVOUQQWiHSYQAAAAGOUAQFqNrFAAAFOMAQVikIw9slgQAAMpwgCAsFu0ziRAAAKAqBwhCDJYBAADlOEAQ4h4hAAAoxxGCUMQ9QgAAUIojBCFahAAAoBhHCEKR3NEiBAAAZThCENpp6wkAAICqHCEIReYhYNQoAAAowgGCsARdowAAoBgEIQAAODUHCEKMGgUAAOU4QhCiRQgAAIpxgCAsQYsQAAAU4whBKDI3HqNGAQBAEQ4QhKUiuaNFCAAAynCAIMRaowAAoBwHCMJSkdwQhAAAoAwHCEKMGgUAAOU4QBCWoUUIAACKcYAgxMoyAACgHIcIQuaORbcBAEAZDhCEZRK5OkAxAQDAITlAwqBrFAAAlKP2IBQZMUYatRcTAAAcldoTBpMIAQBAUQhCAABwaqoPQonc1F5GAABwYGoPmTKRuWLuBAAAKEbtQYiuUQAAUJTqgxBdowAAoCS1hwxahAAAoCi1B2GZRC5qLyMAADgwtYcMtp4AAABFqT4IsdAoAAAoSe0hg3uEAACgKLUHYZnEXHnMIwQAAKWoPghFckWLEAAAFKP6IMQ9QgAAUJLaQ0aL6RMAAKAktYcMWoQAAKAotYcMWoQAAKAotYcMpk8AAICi1B6EaBECAICi1B4ymEcIAACKUnsQokUIAACKUnvIIAgBAEBRag8ZBCEAAChK7SGDIAQAAEWpPWQQhAAAoCi1hwyCEAAAFGWHkMnKynrmmWeioqL8/f1TU1NtP2BFCEIAAFCUfUKmc+fOCxcuzM7OliTJLgcshyAEAABFaWw/hL+//6uvvpqXl2f7oapCEAIAgKLsEITyiaKYnZ2dmJio/9XNzS00NNT0f0EQAgCAomQFYX5+/pdffln18aeffrpJkybyT3bp0qXt27dv2bJF/yvHcZs3bw4PDzf2/MLCwhKtq65ULCiwc48ryFdYWMhxWOWumqEW1AC1oAaW1oKnpyfPm2lOyQpCjuNcXV2rPm726JXExMRERkbOmTNH5vMZY4wXfDxdvb3x5qs2jDFvb+/qLoWzQy2oAWpBDZSoBVlB6O3t/f7779v3xDLpJNKgaxQAABRjn3uEe/bsyc/PJ6IDBw4EBAT07t3bXh0IuEcIAACKsk8QfvXVV6WlpX369Pnpp5+IqGfPnoJgn+10dYw06BYFAADF2CcIN2/ebJfjVIWuUQAAUJTaQwZdowAAoCi1hwyCEAAAFKX2kBEZCbhHCAAAilF7EKJFCAAAilJ7yGCwDAAAKErtIYOuUQAAUJTag1ArMQ0W9wMAAMWoPQhFhq5RAABQkNpDBl2jAACgKLUHIUaNAgCAotQeMhg1CgAAilJ7yGDRbQAAUJSqg5ARMUY8ghAAABSj6iAUJRJUXUAAAHB4qs4Z9IsCAIDSVB2EEnGYOwEAAIpSdRBiEiEAAChN1UGIuRMAAKA0VecMWoQAAKA0BCEAADg1dQehxATMIgQAACWpOwgxahQAABSm7iDEPEIAAFCYuoMQK8sAAIDCVJ0zGCwDAABKQxACAIBTQxACAIBTQxACAIBTU3UQMuIwjRAAABSl6iBEixAAAJSGIAQAAKem7iDEPEIAAFCYqnNGUnn5AADA8ak6aNA1CgAASlN1EEqMMGoUAAAUpe4gxO4TAACgMFUHoSihRQgAAMpSdRBKhHuEAACgLFUHIXaoBwAApak6CDF9AgAAlKbqoMGEegAAUJqqcwYtQgAAUJqqgwbzCAEAQGmqDkJswwQAAEpTdRBiiTUAAFCaqoOQESEHAQBAUaoOQowaBQAApak6ZzBqFAAAlKbqoMGoUQAAUBqCEAAAnBqCEAAAnBqCEAAAnJqqgxDTJwAAQGkqD0KsLAMAAMpSdRCiaxQAAJSm7iBUefkAAMDxqTpo0CIEAAClIQgBAMCpqToIMWoUAACUpu4gZMQhCQEAQEmqDkLsRwgAAEpTdRCiaxQAAJSm6iDEYBkAAFAaghAAAJyaqoOQOOLQOQoAAEpSdRCiRQgAAEpTdRBisAwAAChN1UGIFiEAACgNQQgAAE5N1UGIrlEAAFCauoMQS6wBAIDCVB2E6BoFAAClqToI0TUKAABKQxACAIBTU3UQomsUAACUpuogRIsQAACUpuogRIsQAACUpuogBAAAUJqqg5ARhxYhAAAoSuVBiHuEAACgLHUHIVaWAQAAhak6CDFYBgAAlKbqIETXKAAAKA1BCAAATk3dQYh7hAAAoDB1ByFahAAAoDCN7YdITExcuHBhfHw8Y6x3794zZ8708vKy/bCEFiEAACjPDkEYHx/v4eExd+5cnufffvvtmzdv/vrrr7YflogklbdYAQDA8dkhCMeMGVP+8+zZs8eNG2f7MQEAAB4OO7e4zpw5ExkZaa+jYR4hAAAoTW6LcP/+/VUfbNy4cZ06dcp/TUhImDt37rZt24wd5MKFC3/99dfy5cv1v7q5ua1bty48PNzY87VaKi2VCgokmYUEJRQWFnK4VVvdUAtqgFpQA0trwdPTk+fNNPlkBaEkSR999FHVx2fMmDF48GD9z9euXevfv//ChQs7depk7DhRUVG1atWaMWOG/leO4+rXr2/ivIKmxMPd1dsbNwqrE2PM29u7ukvh7FALaoBaUAMlakFWEPI8v3v3bhNPSExM7NOnz+zZs5955hlTJ9Noatas2aBBA5mFw/QJAABQmh0aW3fv3o2NjX3nnXeef/55249WEaZPAACA0uwQhGvXrs3Ozv7ggw/8/0un09l+WEKLEAAAlGeHIHzjjTeyHqTR2GFWBhERbkwDAIDCVD0OBV2jAACgNHUHIbpGAQBAYeoOQrQIAQBAYeoOQrQIAQBAYQhCAABwauoOQnSNAgCAwlQdhAAAAEpTdRBi9wkAAFCaqoMQAABAaaoOQgyWAQAApSEIAQDAqak7CDFqFAAAFKbqIAQAAFCaqoOQEePQOQoAAEpSdxCiaxQAABSm7iDEYBkAAFCYqoMQAABAaaoOQrQIAQBAaeoOQtwjBAAAhak7CNEiBAAAhak6CAEAAJSm6iBE1ygAAChN3UGIrlEAAFCYqoMQAABAaQhCAABwaqoOQkYcdqgHAABFqToIAQAAlKbqIMRgGQAAUJqqgxAAAEBpqg5CzCMEAAClqToIAQAAlKbqIMQ9QgAAUJqqgxAAAEBpqg5CtAgBAEBpqg5CAAAApak6CDFqFAAAlKbqIAQAAFCaqoMQ9wgBAEBpqg5CAAAApak6CNEiBAAApak6CAEAAJSm6iDEqFEAAFCaqoMQAABAaaoOQtwjBAAApak6CAEAAJSm6iDEPUIAAFCaqoMQAABAaaoOQtwjBAAApak6CAEAAJSm6iBEixAAAJSm6iAEAABQmqqDEKNGAQBAaaoOQgAAAKUhCAEAwKkhCAEAwKmpOggxahQAAJSm6iAEAABQGoIQAACcmqqDENMnAABAaaoOQgAAAKUhCAEAwKmpOggxahQAAJSm6iAEAABQmqa6CwAA8MiaOXPm7t27q7sUj4ilS5e2bt1aiSOrOgjRNQoADi0+Pn7y5Mlt2rSp7oI4vLfffvvWrVvOGIQAAI6ucePGbdu2re5SODw/Pz/lDo57hAAA4NRUHYSYUA8AAEpTdRACAAAoDUEIAABOTd1ByKFnFAAAlKXuIAQAAMfEGNNqtdVdClkQhAAATiQhIeGpp5566qmnRo8ePW3atIMHDyp0ogMHDlSaN5KRkTF27Njk5OTyR1avXj179mwiKikpmTVrVseOHUeMGHH27FmFimSMqoMQo0YBAOwrPT19586do0aNevLJJ0NDQ/v27Xv48OGHc+qAgIC6detOnDiRMUZEt2/fnjJlysCBA4lo1qxZR44cWbx4cffu3WNjYwsKCh5OkfQwoR4AwLm4u7uPGjVK//OePXsOHDjQuXNnIsrPz1+yZMmlS5eaNGny2muveXp6EtGRI0dWrlyZm5tbr1691157LSgo6OrVq0uXLk1LS6tTp86kSZOaNGlCRGfOnPn9999zcnIGDhw4bNgwY6eeM2dOu3btfvrpp0mTJr300ktTp05t165dcXHxsmXL9u7d27p169atW69Zs2b16tWTJk16KBeDSOUtQgAAsDtJkrKzs9PT03fv3n369Gl9ChYXF3fo0CEzM3PUqFF3797t168fYywxMXHEiBGdOnV69tlnAwICMjMztVptjx49IiIiJkyY0LRp05SUFCKKi4sbMWJEkyZNBg4cOHv27MWLFxs7tZub27Jly955550PPvjg/v377733HhHdvHmzpKSkVatW+uc89thjD7l3VNUtQqw1CgCPmK/PSwsvSg/tdF4udHSoxvPBT/rMzMx27doRUUZGRv/+/du3b09Ey5Yti4qK+uyzz4howIAB0dHR8fHxBQUF/v7+AwcO9PPzi42NJaJ79+7l5+cPHTo0PDy8/IAffvjhvHnz9K3M8PDw0aNHv/LKK8aK1L59+5dffvmLL744efKki4sLEaWnp/v5+ZXPEqhZs+bFixftfSVMUXUQAgA8Yl5swg+v9/C+4bvw5FnlYz4wMPDGjRtEVFZWNnbs2OnTpy9evDghIWHfvn0NGzbUPyctLS0xMfHJJ59s0aJFaGho165dn3rqqYkTJwYFBU2bNq1Zs2YtWrQYPny4vgc1ISFhxowZ7777LhFJknT37l39XUBjHn/88WXLljVv3lz/a40aNYqKisr/tbCwUNGVRatCEAIAPDw+LuTjopauLldX1z59+ixatIiI/Pz8xo0b991331V6zqpVq3JycjZv3jx79uzS0tIpU6bMnTv3gw8+iIuL++yzz27durVo0SI/P781a9Z07NjRumJEREQUFRWlpqYGBwcT0fXr1/UN1odG1fcI0TUKAKCc1NTU1atX62/ODR8+fNWqVVevXtX/0+3bt4uKilJSUnJzc/UZ2aNHD/2vKSkpHh4egwcPHjVqlH4uxMiRIz///POSkhIiYoxdvnzZomIEBAT07dt34cKFRHT16tVdu3aNGTPGzi/VJLQIAQCcS1paGsdxHMcFBwfHxsZ+9dVXRPTYY4999dVXvXv3DggIKC4uFgRh3759586dmzBhQmhoqFar9fDwWL9+/f3793v06OHv7+/h4ZGbm7ty5Uoimjt37uuvv96gQYO6deumpqYOGzbs+++/53lefwtQjvnz5w8bNmzVqlWZmZmzZ89u0KCBgq+/CgQhAIAT6d27t7EbeOPHjx8/fnxqaqqHh4f+Ll3//v2TkpLS09Pd3d1r1aqlf5r+EcaYvieTiLy8vJYtW6bVatPS0gIDA93d3Ymoa9euJ0+eNHii2NjY1NTUio9ERUVdunQpKSmpVq1a+mkbD5OqgxAT6gEAHrLyeNPTaDShoaEVH+E4rk6dOlX/o4uLS1hYmNXn5TjOlv9uC1XfIwQAAFAaghAAAJyaqoMQo0YBAEBpdrhHmJKSsmnTpsTERA8Pj759++pX6wEAAHAIdmgRHj169NSpU8HBwaIoDhkyZOnSpbYfEwAA4OGwQ4tw5MiRI0eO1P8cGBi4fPnyF1980fbDEkaNAgCA8uw5fUKr1Z44cSI6OtqOxwQAcFxubm5jx47Vz6sDW6Snp7/wwgsKHdw+QXj58uVBgwbdu3evZcuWO3bsMPa0O3fuHD9+vHweJc/zH374YWBgoLHnM8aVlpaWlNiljGClkpIS+ctDgEJQC2pgRS38+uuvWVlZCpXHqbi4uISGhpaUlFhaC66urjxv5iagrCDct29f7969qz5+7ty5Zs2aEVFkZOTx48dv3749c+bMqVOn/vjjjwaP4+PjExISUr6aKs/zvr6+Jl5SfW9tLQ+Ni2oWqHVOLi4u+AiudqgFNbCiFgICAgICAhQqj3OytBY4GTfYONObZVjq+PHjPXr0qLihRkWffPJJaWnpnDlzZB4tPz/fx8fHfqUDa6AW1AC1oAaoBTVQohbsMGq0pELf5alTp+rWrWv7MQEAAB4OO9wjfO655+7fvx8REZGcnHzy5MkVK1bYfkwAAICHww5B+PPPPx89ejQ5OTkgIKBz586+vr62HxMAAODhsEPXqKenZ69evZ599tkBAwbYNwUXLFhQgjGj1Q21oAbffvttaWlpdZfC2c2fP7+srKy6S+HslKgFVa81umTJkvv371d3KZzd4sWLMzIyqrsUzm7RokWZmZnVXQpnt3DhQsyFqHbff/99dna2fY+p6iAEAABQGoIQAACcGoIQAACcmp0n1Jv27LPPbtmypWbNmjKff/fu3ZCQEEEQFC0VmIZaUIM7d+6EhoaiFqrXnTt36tata3a9LlCUpbUwduxYs6u4PNQgLC4uTkpKkv/HXFpa6ubmpmiRwCzUghqgFtQAtaAGltZCcHCwh4eH6ec81CAEAABQG7TxAQDAqSEIAQDAqSEIAQDAqSEIAQDAqdlnh3q7Y4ytXr06Pj6+Xr16L774otkxP2AvV69e3bp16927d8PDw8ePH6+f6yKK4k8//VT+nObNm3fq1Kn6yvjo27p1a1JSkv7nGjVqjBkzRv9zbm7u0qVLU1JSevXqNXjw4OoroFP45ZdfKq5pGRUVpd9sdfny5eUPtm3btm3bttVRukdZXl7eqVOnrl27FhMTU/GjJj8/f+nSpUlJST169Bg2bFj540ePHl27dq2Pj8/EiRPDwsKsOKNKW4T//Oc/58yZExkZuW3btoovGJQ2ePDgK1euhIeH79+/v0WLFvpVRrVa7eTJk69cuZKYmJiYmIhFL5W2YMGCLVu26K/2nTt39A+KotizZ89jx441bNjw9ddf//7776u3kI+8W7duJf7XzJkzT548SUR5eXlTpkwpfzwnJ6e6i/kIeumll6ZNm/bpp5/+8ccf5Q9KktS7d+8DBw40atTozTff/Prrr/WP79q1q3///nXq1MnMzOzQoYOVy1Mz9cnPz69Ro8aZM2cYY8XFxf7+/vHx8dVdKGdRXFys/0GSpOjo6GXLlukfJKLCwsJqLZoT6dev38qVKys9uGnTpgYNGuh0OsbYjh076tatq9Vqq6N0Tufy5cuurq7p6emMsdTUVFdX1+ou0SNOFEXG2OTJk6dPn17+4Pbt28PCwvTv+b1799apU6esrIwx1qdPn6+//lr/nMGDB3/22WdWnFGNLcIzZ864u7u3bNmSiNzd3bt167Zv377qLpSzcHd3L/+5pKTE29u7/NclS5YsWLDg1KlT1VEupxMXF/fll19u3LhRkiT9I/v27evdu7d+PYqePXump6ffuHGjWsvoLH788cchQ4YEBQXpf5UksTtyFQAABEJJREFUacGCBQsXLkxISKjegj2qDK4as2/fvp49e2o0GiLq1q1bbm7ulStXGGP79+/v27ev/jmxsbHWhYUagzAtLS0wMLD819q1a6ekpFRjeZzT/PnzXV1dhw4dSkQcx/Xp0ycjI+PSpUs9e/b88ssvq7t0j7jGjRt7eXllZGTMnDkzNjZWp9PRg38XGo3G398/NTW1WovpFHQ63fLly59//nn9r4Ig9OnTJz09/cyZM506dfrhhx+qt3jOo+L7n+f5gICAlJSUzMzMsrKy8seDgoKs+6NQ42AZjUYjimL5r1qtFssaPWTr1q374osvdu3apb/ybm5uO3fu1P/T6NGj+/btO2XKFE9Pz2ot46NswYIF+h8++OCDqKioDRs2PPHEE1X/LlxdXaupgE5k8+bNgiD069dP/2tgYOC2bdv0Pw8aNGj8+PEvvPACVh99CAy+/11cXIhI/01R/4N1fxRqrL+QkJC0tLTy15ycnBwcHFy9RXIqGzdufPXVV7dt29asWbOq/9q5c2edTlc+phEU5e3t3bJly8TERCIKDQ1NTk7WP15YWJiTkxMSElKtpXMKy5YtmzBhgsEVkrt06ZKfn4/Nwx+Oiu//kpKSzMzMkJAQX19fLy+v8setDgs1BmGbNm28vLz0TZD09PRDhw5hpPhDs3379pdeemnjxo36e7R6+sEyeps3b/by8oqIiKiO0jkFURS1Wq3+5/T09Pj4+JiYGCIaMmTI9u3b8/LyiGjdunXR0dH16tWrxnI6g7S0tO3bt0+YMKH8kZKSkvKfN23aFBQUVH7vEBQ1ZMiQXbt26fem37BhQ8OGDSMjI4lo6NCh+sGloiiuW7dOfzfHYraP8FHCb7/9FhgYOHHixEaNGk2bNq26i+NEvLy8ateu3fa/lixZwhj74YcfoqOjn3nmmf79+/v4+Cxfvry6i/koS0tLCwoKGj58+OjRo2vVqjVhwgRJkvT/NGbMmOjo6Oeeey4gIGDLli3VW05n8Omnnz7++OMVH5k3b16rVq3GjRvXu3dvX1/f9evXV1fZHmGLFy9u27ZtYGCg/rPo559/1j8+fvz4pk2bTpgwISAgoPzKX7hwITAwcPTo0d26devQoUNRUZEVZ1Tv7hPXrl07ceJE/fr1O3bsWN1lcSKnTp2q+JYIDg4OCQnRarUnT568efOmr69vu3bt8BVYaVeuXLl48aJOp4uJiWnatGn544yxAwcOpKSkdO7cOTw8vBpL6CSuX7/u6elZsQu6tLT0xIkTd+7cqVmzZvv27WvVqlWNxXtUpaSkVBzzEhoaWqdOHSJijB08eDApKalTp04Vu0OysrJ27dpVo0aNXr16WXePUL1BCAAA8BCo8R4hAADAQ4MgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp4YgBAAAp/b/AK7goDcjE2oAAAAASUVORK5CYII=",
"image/svg+xml": [
"\n",
"\n"
],
"text/html": [
"
"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"plot(0:0.01:100, bessely0, label=\"Bessel Y0\")"
]
},
{
"cell_type": "code",
"execution_count": 20,
"id": "dd61e08d",
"metadata": {},
"outputs": [],
"source": [
"# savefig(\"figbesselY0.png\")"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 1.12",
"language": "julia",
"name": "julia-1.12"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.12.4"
}
},
"nbformat": 4,
"nbformat_minor": 5
}