aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers/fm/Python animation/Bessel-FM.ipynb
blob: f01db05522e34c1a554ec32851bdc62e96c4bb66 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from scipy import signal\n",
    "from scipy.fft import fft, ifft, fftfreq\n",
    "import scipy.special as sc\n",
    "import scipy.fftpack\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib as mpl\n",
    "# Use the pgf backend (must be set before pyplot imported)\n",
    "# mpl.use('pgf')\n",
    "\n",
    "\n",
    "def fm(beta):\n",
    "    # Number of samplepoints\n",
    "    N = 600\n",
    "    # sample spacing\n",
    "    T = 1.0 / 1000.0\n",
    "    fc = 100.0\n",
    "    fm = 30.0\n",
    "    x = np.linspace(0.01, N*T, N)\n",
    "    #beta = 1.0\n",
    "    y_old = np.sin(fc * 2.0*np.pi*x+beta*np.sin(fm * 2.0*np.pi*x))\n",
    "    y = 0*x;\n",
    "    xf = fftfreq(N, 1 / N)\n",
    "    for k in range (-4, 4):\n",
    "        y = sc.jv(k,beta)*np.sin((fc+k*fm) * 2.0*np.pi*x)\n",
    "        yf = fft(y)/(fc*np.pi)\n",
    "        plt.plot(xf, np.abs(yf))\n",
    "    plt.xlim(-150, 150)\n",
    "    plt.show()\n",
    "    #yf_old = fft(y_old)\n",
    "    #plt.plot(xf, np.abs(yf_old))\n",
    "    #plt.show()\n",
    "    \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAeIUlEQVR4nO3dfZAc9X3n8fdXu3pCqwfQrh7QAxKRMAiDDbXGduGzsU3OQlDS5ezYECd2HCe6OhtffHYlh487nOBLlR+uEh8XYoJ9rsQ+wqPPjuLIkTGGghgjkBACJCG0EhIroefn593Vfu+P6dntGU3P9O7O9Mx0f15VW5rp6Z35Tu/qM7/99q+7zd0REZHmN6reBYiISHUo0EVEUkKBLiKSEgp0EZGUUKCLiKSEAl1EJCUqBrqZfd/M9pnZqxGPm5ndY2ZdZvaymV1b/TJFRKSSOCP0vwMWl3n8JmBh8LUc+M7IyxIRkaGqGOju/jRwqMwqy4AfeM5zwBQzm1mtAkVEJJ7WKjzHLKA7dH9nsGx3uW9qb2/3efPmVeHlRUSyY+3atQfcvaPUY9UI9NjMbDm5tgxz585lzZo1Sb68iEjTM7MdUY9VY5bLLmBO6P7sYNl53P1+d+90986OjpIfMCIiMkzVCPQVwKeC2S7vAY66e9l2i4iIVF/FlouZPQjcALSb2U7gq8BoAHe/D1gJLAG6gFPAZ2pVrIiIRKsY6O5+W4XHHfh81SoSEZFh0ZGiIiIpoUAXEUkJBbqISEoo0KVu1u44zMa3jtW7jKZw4mwfP1lXcjawyIBEDywSCfvod54FYPvXb65zJY3vv//kVX68bhfz2ifwzjlT6l2ONCiN0EWawO6jpwE41dNX50qkkSnQRURSQoEuIpISCnQRkZRQoIuIpIQCXUQkJRToIiIpoUAXEUkJBbqISEoo0EVEUkKBLiKSEgp0EZGUUKCLiKSEAl1EJCUU6CIiKaFAFxFJCQW6SDPxehcgjUyBLiKSEgp0kWZi9S5AGpkCXUQkJRToIiIpoUAXEUkJBbqISEoo0EWaiaYtShkKdBGRlFCgizQTTVuUMhToIiIpoUAXEUmJWIFuZovNbLOZdZnZHSUen2tmT5rZOjN72cyWVL9UEREpp2Kgm1kLcC9wE7AIuM3MFhWt9t+AR9z9GuBW4G+qXaiIiJQXZ4R+HdDl7tvcvQd4CFhWtI4Dk4Lbk4G3qleiiAzQtEUpozXGOrOA7tD9ncC7i9b5M+DnZvYFYAJwY1WqExGR2Kq1U/Q24O/cfTawBPihmZ333Ga23MzWmNma/fv3V+mlRUQE4gX6LmBO6P7sYFnYZ4FHANz918A4oL34idz9fnfvdPfOjo6O4VUsIiIlxQn0F4CFZjbfzMaQ2+m5omidN4EPA5jZFeQCXUNwkWrTgUVSRsVAd/c+4HZgFbCJ3GyWDWZ2t5ktDVb7MvBHZrYeeBD4fXfX7hsRkQTF2SmKu68EVhYtuyt0eyNwfXVLExGRodCRoiIiKaFAFxFJCQW6SDPRnikpQ4EuIpISCnQRkZRQoIuIpIQCXaSZ6MAiKUOBLiKSEgp0EZGUUKCLNBNNW5QyFOgiIimhQBcRSQkFuohISijQRZqJpi1KGQp0EZGUUKCLiKSEAl2kmWjaopShQBcRSQkFuohISijQRURSQoEuIpISCnQRkZRQoIs0EU1ykXIU6CJNxJXoUoYCXUQkJRToIk3E1XSRMhToIiIpoUAXaSLqoUs5CnQRkZRQoIs0EQ3QpRwFuohISijQRZqIq4kuZcQKdDNbbGabzazLzO6IWOfjZrbRzDaY2T9Ut0wRAbVcpLzWSiuYWQtwL/CbwE7gBTNb4e4bQ+ssBL4CXO/uh81sWq0KFhGR0uKM0K8Dutx9m7v3AA8By4rW+SPgXnc/DODu+6pbpogAGqJLWXECfRbQHbq/M1gWdhlwmZn9ysyeM7PF1SpQRETiqdhyGcLzLARuAGYDT5vZVe5+JLySmS0HlgPMnTu3Si8tkh069F/KiTNC3wXMCd2fHSwL2wmscPded38DeJ1cwBdw9/vdvdPdOzs6OoZbs4iIlBAn0F8AFprZfDMbA9wKrCha5yfkRueYWTu5Fsy26pUpIqBD/6W8ioHu7n3A7cAqYBPwiLtvMLO7zWxpsNoq4KCZbQSeBP7E3Q/WqmiRrFKgSzmxeujuvhJYWbTsrtBtB74UfImISB3oSFGRJqIBupSjQBcRSQkFukgT0blcpBwFutTc4ZM9rN1xuN5lpNKuI6d5bc+xepchDUKBLjV323ef46PfebbeZaRC8fj8+q//ksXffqYutUjjUaBLzb2253i9S0gNdVykHAW6iEhKKNBFmoqG6BJNgS6J0QwNkdpSoIs0EX0mSjkKdEmMwkikthTokhjl+chpG0o5CnRJTLiHrn760BgG6K8cKU+BLolRFg2frlQkcSjQJTHh0aVGmsOjYJdyFOiSmDhhdK7f+S+PvcyWvdk8ujT//rv2ZfP9y8go0CUxBSP0iHW27j/Bw2u6+dwDLyZSU6PZsu84D6/p5vMPrCv5uP6ykXIU6CIiKaFAl7qImuWiEWhOVHtKm0fKUaBLYoYS1ma1q6OR5acnigyHAl0SEx51RmW7ZnHkRH34af6+lKNAl8Ro2uLwaeQucSjQJTHK8OHTXy4ShwJdElNw6H+FgMrqiLTSvgP9ZSPlKNAlMXGySIGVo80gw6FAl7qoFNzZneVSnlovUo4CXRITZ/StEbrI8CnQJTkK6xHTB56Uo0CXxBTMQ1cwlVU831znQ5c4FOiSGIVRZVH7DtQ7lzgU6JIYL7gdda4SBVc52jpSjgJdEjOUw9Yto9Nc9FeMjESsQDezxWa22cy6zOyOMut91MzczDqrV6KkUfS5SpKto9FUevs6l4uUUzHQzawFuBe4CVgE3GZmi0qsNxH4Y2B1tYuUdFAUVZbPa20rGY44I/TrgC533+buPcBDwLIS630N+AZwpor1SYrEuWJR1lXah6DtJuXECfRZQHfo/s5g2QAzuxaY4+7/XMXaJGUKpy2Wv8BFNjvooQ+9os0zcG4bJbqUMeKdomY2CvhL4Msx1l1uZmvMbM3+/ftH+tLSbOIcKZrxxFKLXEYiTqDvAuaE7s8OluVNBN4OPGVm24H3ACtK7Rh19/vdvdPdOzs6OoZftTQlj7gd1p8foWd0iF655aLEl2hxAv0FYKGZzTezMcCtwIr8g+5+1N3b3X2eu88DngOWuvuamlQsTSveuVw89rppFDn7R0EuMVQMdHfvA24HVgGbgEfcfYOZ3W1mS2tdoKRHnEP/+zOeW5VmuWT1g07iaY2zkruvBFYWLbsrYt0bRl6WZFcusdRyERk6HSkqifEYTfSsj0AHRuhRs4ASrEWajwJdEhPnXC7aKVqazrYocSjQJTFxDlvP+qHt/Z5vOWX0E01GRIEuiSk4UlQ7RUuq3HLJ+AaSshTo0lAUWApyGT4FuiQmzrlcBg/9z2bLQdMWZSQU6JKYOKPMrAdWxt++jJACXeoiqkc8uFMwyWoaR3+FnQgKfClHgS6JidVySaSSxlXx/Wf9TxgpS4EuiYkTRf0ZD6zBWS71rUOakwJdEhNus0QGVsaDLGo/w8CBRUkWI01HgS6JGcoIPaMt9IGNpGmKMhwKdElMYQ+9/BWLsqryRaKjlmd8wwmgQJdEVQ6dgR56Rqe5RO1DqHjhC+W5oECXBMU622IilTSuSsGsszBKOQp0aShZbx0M991nfbtJjgJdEhPnmqJZz6VKl+DT5CApR4EuiYk622J4dDlwPvSEamo0Oh+6jIQCXRIT61wuGR9rDrd1oqAXUKBLgqKmLYaXZ/6KRZV2ikYuV6KLAl0SFKcvXKmHnHZRh/5r2qLEoUCXxESNygvWyXgwVT6wKOMbSMpSoEti4hzlmA/97LZc1EOX4VOgS11o2mJpw72mqnroAgp0qZPCUfmgrE9bjPqo07RFiUOBLomJE0ZZ7xEP9+1nfLNJQIEuiYnaKRp1O4vybz/6nC06l4tEU6BLYmKN0DMeTcMfoWd7u0mOAl0SE+egmMEDi7LZRR/ufHPFuYACXRIU2UZQy2VA1CwXHVgkcSjQpS6iAijrl6Cr1DqJfFSBLijQJUGFp8/Vzr1yireDpi1KHAp0SYymLVY27J2i+igUYga6mS02s81m1mVmd5R4/EtmttHMXjazJ8zskuqXKs1P0xYrqdgr18W1pYyKgW5mLcC9wE3AIuA2M1tUtNo6oNPdrwYeA75Z7UKl+cUJnYEeekab6P39w/s+5blAvBH6dUCXu29z9x7gIWBZeAV3f9LdTwV3nwNmV7dMSYOoS9DFOQtjVgweWBTxeIwTnEl2xQn0WUB36P7OYFmUzwI/G0lRkk5DGqFndJ7LsM+2WOU6pDm1VvPJzOx3gU7gAxGPLweWA8ydO7eaLy1NoOCEXK5ReSkDI/QhRrS2oUC8EfouYE7o/uxgWQEzuxG4E1jq7mdLPZG73+/une7e2dHRMZx6JSWi2i/9WU+mSpegy/r2kbLiBPoLwEIzm29mY4BbgRXhFczsGuBvyYX5vuqXKWkQJ4qynlfD/UDTtEWBGIHu7n3A7cAqYBPwiLtvMLO7zWxpsNq3gDbgUTN7ycxWRDydZFj0VMXzz+WS0RZ6jEvQDfMbJRNi9dDdfSWwsmjZXaHbN1a5LkmhOKPIrLdc8m9/qFcuyvZWkzwdKSrJ8dJ3wov7zuXujW7J5hA9/6HXH5HounSflKNAl8TEyZxzwZE12Z22mPu3b4hDdPXQBRTokqA4h/v3BkGW1dZL/l1HjtCHeMCRZIsCXRITOYoMLe47lxuhZzWg8juIzxVtAI3AJQ4FutRFVDz1Bj30rAZYpZaLTjss5SjQJTGRLZdQHPX1a4QOcK4o0CudD10HHAko0KWKTp7t408eXc/RU70lH4+3UzQ/Qs+m/PsuDvSK31e0+j1PbOH5Nw5VpyhpGgp0qZoHn3+TR9fu5N6nuko+XnD+logzLOZbLllN9OiLQA/tg+4vH3+dj//tr6tTlDQNBbpUTT6MIi8GHeM58jtFszrLZdiH/mdzc0kRBbpUTcWLUoRCJ3whh3AW9Wa85RJl4AjSIe4slWxRoEtiwqHTc670pXkGpy1mM6CiznGT1xux3cKiQl/ST4EuVRfn4JdwMIWDSztFB9/5wP4EBrdH1AdheNsWz2GX7FCgS1309IUCPbR8YB56USa5OwdOlDzNftM6cbaP0z3nCpaF33dfib5UeLsVfF/o9lBnyEh6KNAlMeGwCgdTeLQ+OA+9MJS+98wbdP6PX9B96BRp8favruK9X3+iYFmpDzcYHJlHtVzC22uo54GR9FCgS9XFOWV3uHUQHqUOHilaaNWGPQDsOXamChU2jiNFc/bDs1z6QtvoTG9uG8UaoZ9ToGeVAl2qLs7RjOFgOtMbGqFHTFvMfwCMbkn3r2xBLzw00h4I9Bg99IJWjWRKuv93SEMpGKGHAv107+AIPR9i4ZAPrz8q5WfVPRPaFr2hQD89MEKv/PePeujZpUCXxBT00Eu0E2Cw5VK8szAf6L0pbyecCr3vvhJtqagReph66NmlQJcEhafklQ70fLsgPGqHwSDrixFozSbcijpd4sMN4EzwgdbTV7hdBp9j8HZfyj/0JJoCXarGgkNFI0/xGjHLpVSInerpK/jeNI/Qw+8//JdJ/j2f6/eK77/gMn7qoWeWAl3q4mxf6VkuR0/nZn2c6e0vOOKxt8K0vWZ24uzgh1d4Wxw/k98W54d8saidqZItCnSpmkqH6xfOsT6/5dLf7xw51TNwgegzfecHWZwecrMJh/ip3nOMCWbyHAk+3Ar/gomatqh56KJAlyoaOGw/xqH/pVouJ3r66HeYOXl8bnm4/ZCyEXp4FB3eEXqm5xwzp4wD4MipntyyiBF6wemINUIXFOhSRZVGhgUn5+rrZ9zoUYxpGcXBk7ngyl8YY8bkXKCdKnHAUVoCPfw+ThWM0PuYMSkf6LntcfBEbvtMHj+6oFUVdQUojdCzS4EuVVPpLH/h0Dlxto+2saOZPnkse4/mjv48HIxILy4R6Hlp2SnaV+KgIci95/a2sYxuMQ4Hgb472D4LprUV9NujtsQ57RTNLAW6VE3fQMul8kyMAyd6aG8bw8xJ4wcC683gPC1Xz54SrHP+ybhSM0LvKz1CP3D8LO1tY7h4yni6D+e2x56jpwF4+8WTOHq6d6DtEnUFKE1bzC4FulRNvncbdfrWcAAdPHmW9raxzLpwPNsPngSga98JzOD6Be3A4Mg03BPujZjl0WwKWy65UffJs30cO9PHjMnjWdDRxtZ9JwDYfvAU40aPYsH0iQAcClpU4a2sHrqAAl2qKD9CjxMoB0/0MLVtDO+YPZm9x86y68hpXuo+wvypE7hk6gXA4Mg0PCc9LS2XUicny594bObkcVw+cyJd+05w7Ewv6948zNWzp9DRNhYY/MslauezeujZpUCXqsn3biv9yX+u39l3/AwdbWN538IOAP7nqs0823WQD14+jXGjW5g2cSzb9udG7ifPnj/bpdmFP5jyLZf8+51z0Xg+dPl0+vqdv/jpJtbvPMq/WdBOx8RcoO85ev4ZJzVCF1CgSxVVGqHnQ2fHwZOc6e3nshkTWTCtjSVXzeDH63YxtnUUn7l+HpDro6/rPgLA7mCkDunpD4dPYbA3GJm/1H2YllHGopmTuXbuFN63oJ2H13QzdcIYbnv3XBZObwNg897jQGHfXPPQBaC13gVIc/pV1wHmtU9g1pTxA8vys1wiT/EahM7LO48CsGjmJAD+6hPv5Jar93HlxZOYfWGu3XL9gqn8YtNeXttzjA1vHRt4jmNnekmD8Pt49a2j9Pc7//LqHq6dO4XxY1oA+O6nOnnitb1cO/dC2oN2y5yLxrM++KCLmrYYteO4v99Zsf4tbrl6Jq0pPw1xVumnKkPm7nzye6u55Z5nCpbnR4a7S7QEct+X+3flK7tpbxvLFUGgj21tYclVM7lk6oSBdZe+42LaxrZyx49e4f8+t4NZU8Zz2fQ2dgQ7UPPWbD80cIh8o3J3nt16oGB64vYDuRksH7p8Gqu3HeJPf/QyW/ef5JPvvmRgnfFjWrjl6ou5OPSh+YHLOnhmy4GBHaOlvHXkdMnlj63dyRcffokfPrdjpG9JGlSsQDezxWa22cy6zOyOEo+PNbOHg8dXm9m8qlcqDSN/wMvhoqvt5A96eePAyfO+B3KB/sSmvTy+aS+3vmsOLWVObj61bSzf/NjVvLbnGFv3n+DOm69gfvuEgT4z5A5E+th9v+aW//2vI31LNfXwC938zndX8+ia7oFlbxw4ySiDP196JTMmj+OxtTv5eOdslr7j4rLP9an3zqPnXD93/9OGgtZWuMmS3/5tY1sLZhbtCoJ+//F0XZtVBlVsuZhZC3Av8JvATuAFM1vh7htDq30WOOzuC8zsVuAbwCdqUbDUX9QIfGPQGjl0soenX9/P+y/rKHj8a/+8kWOne1k0cxKf++BvVHydJVfN5P2XdXCu35k8fjS7Dp9m1Ya9bNp9jCtmThroJe84eIpz/V72A6Kent16EICtwYdR37l+Vm3YwzvmTGHORRfwyy/fwNHTvQM7Pcu5bPpEPv/BBdzzxJaBfQwwOCX00Mkent6yH8gdvNV96DRzg1lDadmhLNHi9NCvA7rcfRuAmT0ELAPCgb4M+LPg9mPAX5uZeaWzNUlDc3de23OcE2f7eNe8i+jvdw6e7OHJzfsG1tl77AxTLhjNyld2s3H3Mf7DBy7lZ6/s4VPff545F41n4bSJA+vOmjKeT7xrDl/40EIuGBNv903b2MH1/v21s/ibp7r43e+t5reumcXG3YO99cc37uH6Be1cMKaVw6d6WL3tEB+5cnriveLn3zjEpPGtvG36RM729bP/+FnW7jgM5FpNbWNbeWbLfrbsO8Ff/841AIxpHRUrzPP+840Lmd9+AQ+u7mbHwVzrZvkP1zJ5/GjePHiKfnfuXnYld/3jBr79xOvcdcsiJo8fzYtBHS++eZjuQ6eYMXkco1tGsefoGV7ddZQPXj6tYT8UJR6rlLlm9jFgsbv/YXD/94B3u/vtoXVeDdbZGdzfGqxzIOp5Ozs7fc2aNUMueH33EZbd+6shf58ko2PiWNydAyfO7/Fe2jGhxHcMTbjl0uxquT0u7ZiQqm2VNqu++H7eNmNi5RVLMLO17t5Z6rFEZ7mY2XJgOcDcuXOH9RxTLhhdzZJkGCaNa+W9vzGVbftPsiU4mnHc6FHccNk0WlvOH+Ft2n2MaRPHMbVtzIhfOz8zJj8MMeDYmT6efn3/iJ+7Vi4PpmdCru5qjoEvnzGRF7Yf5l3zLmSUFT7zopmT2PDWsYJ9GjdfNTPYL6Gwr6eJ42oTvXGedRcwJ3R/drCs1Do7zawVmAwcLH4id78fuB9yI/ThFHzJ1Als//rNw/lWEZFUi9NgfAFYaGbzzWwMcCuwomidFcCng9sfA36p/rmISLIqjtDdvc/MbgdWAS3A9919g5ndDaxx9xXA/wF+aGZdwCFyoS8iIgmK1chx95XAyqJld4VunwF+u7qliYjIUOhIURGRlFCgi4ikhAJdRCQlFOgiIimhQBcRSYmKh/7X7IXN9gPDPY9nOxB5WoE6Ul1D16i1qa6hUV1DM5K6LnH3jlIP1C3QR8LM1kSdy6CeVNfQNWptqmtoVNfQ1KoutVxERFJCgS4ikhLNGuj317uACKpr6Bq1NtU1NKpraGpSV1P20EVE5HzNOkIXEZEiTRPoZvZlM3Mzaw/um5ndE1yY+mUzuza07qfNbEvw9enoZx1RPV8LXvclM/u5mV3cIHV9y8xeC177x2Y2JfTYV4K6NpvZR0LLy14EvEp1/baZbTCzfjPrLHqsbnWVqDPx1wy99vfNbF9wBbD8sovM7PHgd+ZxM7swWB75e1aDuuaY2ZNmtjH4Gf5xI9RmZuPM7HkzWx/U9efB8vmWu1h9l+UuXj8mWJ7oxezNrMXM1pnZTxOry90b/ovcxTNWkZu33h4sWwL8jNwFYN4DrA6WXwRsC/69MLh9YQ1qmhS6/Z+A+xqkrn8LtAa3vwF8I7i9CFgPjAXmA1vJnQ65Jbh9KTAmWGdRDeq6Angb8BTQGVpe17qKakz8NYte//3AtcCroWXfBO4Ibt8R+nmW/D2rUV0zgWuD2xOB14OfW11rC56/Lbg9GlgdvN4jwK3B8vuA/xjc/lzo/+mtwMM1/nl+CfgH4KfB/ZrX1Swj9L8C/pTBK49B7sLUP/Cc54ApZjYT+AjwuLsfcvfDwOPA4moX5O7HQncnhGqrd10/d/e+4O5z5K4wla/rIXc/6+5vAF3kLgA+cBFwd+8B8hcBr3Zdm9x9c4mH6lpXkXq85gB3f5rc9QTClgF/H9z+e+DfhZaX+j2rRV273f3F4PZxYBMwq961Bc9/Irg7Ovhy4EPkLlZfqq58vY8BHzazmlwV28xmAzcD3wvuWxJ1NXygm9kyYJe7ry96aBbQHbq/M1gWtbwWtf2FmXUDnwTy54eve10hf0BupNRodYU1Ul313halTHf33cHtPcD04HZdag3aAdeQGw3XvbagrfESsI/cIGkrcCQ0qAm/9kBdweNHgam1qAv4NrlBaH9wf2oSdSV6kegoZvYLYEaJh+4E/iu5NkLiytXl7v/o7ncCd5rZV4Dbga82Ql3BOncCfcADSdQUty4ZPnd3M6vbtDQzawN+BHzR3Y+FB5H1qs3dzwHvDPYV/Ri4POkaipnZLcA+d19rZjck+doNEejufmOp5WZ2Fbm+6vrgl2c28KKZXUf0xat3ATcULX+qmnWV8AC5Kzp9tRHqMrPfB24BPuxBY65MXZRZXtW6ItS8rirVUi97zWymu+8O2hb7guWJ1mpmo8mF+QPu/v8aqTYAdz9iZk8C7yXX4mkNRrvh1451MfsquB5YamZLgHHAJOB/JVJXLXcKVPsL2M7gTtGbKdzx8nyw/CLgDXI7Hi8Mbl9Ug1oWhm5/AXisQepaDGwEOoqWX0nhzsdt5HYCtga35zO4I/DKGv4Mn6Jwp2hD1BXUkvhrlqhhHoU7Rb9F4Y7Hb5b7PatRTQb8APh20fK61gZ0AFOC2+OBZ8gNZB6lcOfj54Lbn6dw5+MjCfw8b2Bwp2jN60rsF7VKG2c7g4FuwL3kemavFIXEH5DbudYFfKZGtfwIeBV4GfgnYFaD1NVFrh/3UvB1X+ixO4O6NgM3hZYvITdzYSu59kgt6votcn3Ds8BeYFUj1FWizsRfM/TaDwK7gd5gW32WXC/1CWAL8AuCQUC537Ma1PU+cjsbXw79Xi2pd23A1cC6oK5XgbuC5ZcCzwf/Fx4FxgbLxwX3u4LHL03gZ3oDg4Fe87p0pKiISEo0/CwXERGJR4EuIpISCnQRkZRQoIuIpIQCXUQkJRToIiIpoUAXEUkJBbqISEr8f7DBJ3GnZJ8oAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABGbUlEQVR4nO29eZwjV3nv/X2qtPS+zfTsM+7xigdjwJ54AcJusCGxAwRiEt/ADQkhCTcJITeBS14gLG+C74VwSeDNJUCIkwtO4hBiiIkxBC8sxp4xxtt4Fo/Hnp6lZ+t9kVRV5/2jFlWru6cl1Wl1qft859OfkUpS6ZR06qenfuc5zxGlFAaDwWBYHVjL3QCDwWAwNA4j+gaDwbCKMKJvMBgMqwgj+gaDwbCKMKJvMBgMq4jMcr3x2rVr1cDAwHK9vcFgMDQlu3fvPqWU6q/39csm+gMDA+zatWu53t5gMBiaEhF5Jsnrjb1jMBgMqwgj+gaDwbCKMKJvMBgMqwgj+gaDwbCKMKJvMBgMqwgj+gaDwbCKMKJvMBgMqwgj+k2K8jwe+95duE5puZtiMCxKqTDD4/d8F1PKffkxot+knDh0kDv/+n/z7GOPLHdTDIZFOfjQLv7jc3/ByPGjy92UVY8R/SbFdRwAPNdZ5pYYDIvjBVeknusuc0sMRvSbFv8yWXnmctmQfkJbx9g7y48R/SYlFHulvGVuicGwOEb004MR/SYlEntzDhmagEj0PROkLDdG9JuUcuRkTiJD+gn7qYn0l5+qRF9ErhWRvSJyQETeN8/jbxeRkyLycPD36/qbapiFuVw2NBHR2JPpr8vOovX0RcQGPgtcAwwCD4rI7UqpJyqe+o9KqXcvQRsN82Aulw1NhQlSUkM1kf4VwAGl1EGlVBG4FbhhaZtlWIzyQK45iQzpp2zvmCBlualG9DcDh2P3B4NtlbxJRB4RkdtEZOt8OxKRd4rILhHZdfLkyTqaawhRGNE3NA9RNzXdddnRNZD7DWBAKXUpcBfwd/M9SSn1eaXUTqXUzv7+upd4NACEto4RfUMTYCL99FCN6B8B4pH7lmBbhFLqtFKqENz9AnC5nuYZFsLkPRuaiciONJMJl51qRP9B4AIR2S4iOeBG4Pb4E0RkY+zu9cAefU00zEdk75iBXEMzEEb6xt9ZdhbN3lFKOSLybuBOwAa+pJR6XEQ+AuxSSt0O/K6IXA84wBng7UvYZgNmINfQXET91ET6y86iog+glLoDuKNi2wdjt98PvF9v0wxnw0zOMjQT4RWp6a/Lj5mR26SYMgyGZiJK3jFXpsuOEf1mxUT6hiaiHOkb0V9ujOg3KeUZueYkMqQfk22WHozoNymmtLKhqYgGck1/XW6M6DctJnIyNA+RvbPM7TAY0W9aopPI2DuGJqBcNsRE+suNEf0mRZliJoYmwszITQ9G9JsVU1rZ0ESYICU9GNFvUkw2hKGpUMaOTAtG9JsUI/qGZsLMIE8PRvSbFLNylqGZiPqpiVGWHSP6TYpZaNrQTJTLMJggZbkxot+kmIWmDc2EKcOQHozoNyvG0zc0EcaOTA9G9JsUMzBmaCqUmZGbFozoNykme8fQTER2pIn0lx0j+k2KMnnPhibCBCnpwYh+s9IE9k7xmWfwCoXlbsaqwJuepjg4uNzNWBAj+unBiH6TEp08KT2JVLHIwV94AyP/8i/L3ZRVwZl/+AcOvekXl7sZC2JEPz0Y0W9S0r4wulcsoaan8cbGl7spqwJvbAx3dDS1/SEayE3xlelqwYh+05Ju0cd1AFCeu8wNWR0oNxDTlA6UlueVLG87DEb0m5ZyPf2UnuRuIPZmoLkxhP3BTeePrDKRfmowot+kpN0jVY4f6ac18lxxhGIafu4pIxqCMkHAsmNEv0lJ/eQsx9g7jSS0d1RqRd+UYUgLRvSblHL2zvK2YyGMvdNgUm7vkPJss9WEEf1mJeWRftneSakIrTCiK6q0Rvqe8fTTghH9JiX1C6MHEWdq27fSCFN4Uxrpp30MajVhRL9JUaQ80o/snXS2b8UReuZOukXf2DvLjxH9JiXtk7OUGchtKFGevptSe8dE+qmhKtEXkWtFZK+IHBCR953leW8SESUiO/U10TAfqY+cIk8/pe1baXhpz94x9fTTwqKiLyI28FngOmAH8FYR2THP8zqB3wN+rLuRhvlId+Rk7J0Gk3J7J/pRWuZmGKqL9K8ADiilDiqlisCtwA3zPO+jwCeAGY3tMyxAZO+kNJI29k5jaRp7xwQBy041or8ZOBy7PxhsixCRy4CtSql/P9uOROSdIrJLRHadPHmy5sYayqR+WrvJ028sKc/TN5Oz0kPigVwRsYBPAe9d7LlKqc8rpXYqpXb29/cnfevVTcoHxkKbwUT6jUGl3tOvvGFYLqoR/SPA1tj9LcG2kE7gEuBuETkEXAXcbgZzl5a0Xy4r1wzkNpSwH5hI37AI1Yj+g8AFIrJdRHLAjcDt4YNKqVGl1Fql1IBSagC4H7heKbVrSVpsAJrg5DEF1xqLSnekn/YZ5KuJRUVfKeUA7wbuBPYA/6SUelxEPiIi1y91Aw3zk/rSysbeaSjlgmvp/LxTP4N8FZGp5klKqTuAOyq2fXCB5748ebMMi5NyT9/YO43Fa47sHZO0ufyYGblNSnlGbjoj/XL2Tjojz5VGeEWV3uyddAcpqwkj+k1KeSA3nSdR2d5JZ/tWHGEQkFJPP+39dTVhRL9JSXvkVLZ3UnolstJIefZOtLKXsXeWHSP6TUrqPVIzI7ehpD5P30t3ivFqwoh+s6KaI3vHDOQ2CC/l2TspvzJdTRjRb1LSfhIZe6fBpD57x0zOSgtG9JuUtIt+eeWsdEaeKw2V+kg//D+l/XUVYUS/SUm7R2rsnQYTFVxLd6RPWlOMVxFG9JuUtF8ul+2ddEaeKw2V9uydlJcCX00Y0W9W0l61MMreSWn7VhphpF9Kd6Sv0ppttoowot+kpD7Sd8zKWQ0l9fX0TaSfFozoNykq5VULQ3vHDOQ2BpX27B0v5YkHqwgj+s1K2iMnM5DbWNKevUOUvrO8DTEY0W9W0p6yGdkMKbUbVhrlgmvpjPSjH6WUXpmuJozoNynR5XxqRT+wd1LavhVHeEWV0h/ZtAcpqwkj+k1KOXknpZGTWTmrsaQ+e8fYO2nBiH6z4jVH9o4ZyG0MKu3ZOynvr6sJI/pNStpT4MzKWQ0m7dk7Ke+vqwkj+k1K2lM2MXn6jSXl2TsxQ3JZW2Ewot+0pH1gTJmCaw1Fpb32TsprRa0mjOg3KeUCVmkVfWPvNJRQTFMa6ad9Bvlqwoh+s5LySJ8wiySlA4srjdSvnJX2/rqKMKLfpJQHxtJ5uRzZO+Ykbwypz94xop8WjOg3KWmPnMzKWQ0m7dk7pLu/riaM6Dcpqc97Nnn6DSXtK2dFP0omCFh2jOg3OylN2YxsBjOQ2xjSnr0TXpmalM1lx4h+k5L2SD8Sn5R6zCsNlfrsHTM5Ky0Y0W9S0u7ph9k7qW3fSiPt2TspD1JWE0b0m5WUR05leyed9tOKI+3ZO5HYp7O/riaqEn0RuVZE9orIARF53zyPv0tEHhWRh0Xk+yKyQ39TDXHSXobB2DuNpWzvpDPSJ+UpxquJRUVfRGzgs8B1wA7grfOI+leUUs9TSr0AuBn4lO6GGmaT+hmOjsnTbyhNEumb/rD8VBPpXwEcUEodVEoVgVuBG+JPUEqNxe62Y67hlpzI1klp5GTsnQaTetFPeZCyishU8ZzNwOHY/UHgysonicjvAH8A5IBXzrcjEXkn8E6Abdu21dpWQ4xosssyt2MhygujG9FfamZ9xim1d6IgxYj+sqNtIFcp9Vml1HnAHwN/ssBzPq+U2qmU2tnf36/rrVcnaa9aWDIzchtG7DNObfaOsXdSQzWifwTYGru/Jdi2ELcCv5CgTYYqSPvlslkYvXHEf/jTau+Ygdz0UI3oPwhcICLbRSQH3AjcHn+CiFwQu/t6YL++Jhrmo7zkaLpFP63tW1E0g70TZZmZ/rDcLOrpK6UcEXk3cCdgA19SSj0uIh8BdimlbgfeLSKvBkrAMPC2pWy0IXYSpTVyMgujN44miPTNjNz0UM1ALkqpO4A7KrZ9MHb79zS3y7AIafdI49k7SilEZHkbtIJpBnvHlFZOD2ZGbrPSLKIPJmNjqWkieyet/XU1YUS/SUlzpK+UglIJwug+pdHnSiH6gRVJZfZOvI+msb+uNozoNynlAlYp9MyDtkkuB5gTfckJPl/J5dJp78wS/RT211WGEf1mJ4UDY6HwSDbrbzCDuUtL+CObzabS3pn1o28CgGXHiH6TkupIPxCeMNI39s7SotzylVUaI/14HzVXfcuPEf0mpezpL3ND5iGK9DXZO57y+NAPP8STZ55M3LY08fCJh/n4/R9PLoQqZqd5XuomQMXTNE3K5vJjRL9JKWdDpOsEh3IpAF32znhxnK/t/xo/OvqjpE1LFfcduY9b996KqxJG53F7B1J3ZTU70k9ff11tGNFvVtI8I7dC9JNaDo7n76/klZK1K2WEx5NU9FWF6Kctg2e2p7987TD4GNFvUqKISanUCX+lvZPUg/KCYy26xUT7SRslNxB9T1OkH9ppKYv0TfZOujCi36SkOSNCORWin1CEwkjYRPrzM+dHNsWRftoClNWIEf0mJdUnkhtm7wR2Q8LBu9DeWWmRfng8iSP9KE9fj52mm1nff9r66irEiH6TkmbRnzOQm/CSfqVG+kXPF31HJYzMU+/pm5TNNGFEv1nx0nsiKc15+mEkvNJEX5enH+bpW7m8f7+UNtGPp2waT3+5MaLfpMRlPm2DY6HoW5ry9MNIeMXZO0GknzhlU80eyA3ttdQQF/1lbIbBx4h+kzIrYkrbhJfI3tET6a/47B1dA7mptXfinn66ApTViBH9JkWlOA0usnfygd2Q8JJ+xdo7nq6UzWAgN/y80yb6cSsybQHKKsSIfoNQJZfx+47o6/RpHsgtVXj6uuwdb2VF+lH2jmZ7J9WeviaDRzke4/cNRuMZhuoxot8gxu4eZPTfDzL10JCW/TVF9k5e80CuuzIj/TAltV4ieyf8vJ10fU6zsnc0BT3j9x1h9N+fZvKB41r2t5owot8glON3fHdCzwmZ5jQ4FYhONJCb8ERf6SmbySN9//ONPu/U2Tv6rUh3rDBn34bqMKLfIKys/1Grop6JM8pTiBXsM21pcJUpmwlP9JU6OWvJyjCkTPTDnB2xLG0BiiqGaaq2lv2tJozoNwgJOmfYWZOjsALRT9ssx8oyDElniEbZOyvM09dXhqFS9NM5I9eyLG19VZXCPmYkrFbMJ9YgItEvaYz07WCfqRP9ipRNXfbOCvP0tQ/kZkPRT9fnFFo6YtvaI30sI2G1Yj6xRhEG5QVNoq/KkX7a7J1QdMIUQhLaFyu1tHLk6SeekRsO5Aafd8rsnVDoLcvS1lcjm9R4+jVjRL9RBJ3TK2kSaKWwrCDST9s8R6ey4Jqe2jsr1dNPXnunouBa2kQ/snf0+e9eIPpmILd2jOg3CjdY3lDXQK7yyvZOyjr+nNo7CS/pV+rkLF2Rftrz9KOBXNvWGOkH+3HT1febASP6DSIUZl0DuXF7J3UDuaXZtXeS5umvxNo7SqnIttJVhsFK+YxcS2v2Thjpp8vabAaM6DcIpT3Sj6Vsasp99jw9q3Apd3akn/RKZCVm78SvWnTl6essuKaUwtN0BRn2KdGZvWM8/boxot8oPP2ib2m0dwrTDl94z70M7hlOvK85efoJ7YvQ/nA8J3WZSvUSv2rRNpCrMU//4MMn+dJ/vw9HQ3+NBnI1Zu94wRWzMvZOzVQl+iJyrYjsFZEDIvK+eR7/AxF5QkQeEZHvisg5+pva3ISd09N1EmmO9GcmSpQKLmOnpxPvq+zp6ym4Fh/oXCm+vtZIPxzIzerz9MdPz1CYdChpyDaLFm63LH3FAR3j6dfLoqIvIjbwWeA6YAfwVhHZUfG0nwA7lVKXArcBN+tuaNMTRvozGiJ9NTsbQkfwFI05aLhqmFNwLWmefiwSXimiH4/0k9bemTOQqyHS98IgRUd/iPVXLf0rVmQtbUkMzUA1kf4VwAGl1EGlVBG4Fbgh/gSl1PeUUlPB3fuBLXqb2fxEolryojo8de+L8uUy6In0tZ7kjgO2jdhWuPNE+4tHwitlMDc+PqGtnn5e3+SssB94OiLpuL2jIb3YiwdOJtKvmWpEfzNwOHZ/MNi2EO8AvjXfAyLyThHZJSK7Tp48WX0rVwKxzunNJKyqGF7Oa8ze0XmSK6eEZDIQXokktXdikfBKEf1Z9o6mevpRtpSGSF/rlV98IFfH/qbLx2ci/drROpArIjcBO4H/Od/jSqnPK6V2KqV29vf363zr9OPFRT9hZKdmi76OtLXw5NGSseG4/hwCS4Kd68negZVj78RLSugrwxBOzkpuIUZXfjqCAKXX058VNBnRr5lMFc85AmyN3d8SbJuFiLwa+ADwMqVUQU/zVg6zfMjppB5upaevIdJ3NUZ2jgPZbPlHKWk9/bi9s0LSNnUO5EZ9y7Yhk9Hj6WsMAuIzcnWMP8VF30T6tVNNpP8gcIGIbBeRHHAjcHv8CSLyQuD/ANcrpU7ob+YKIL6kbVJ7J4iWLDvM3kmbvePMsneSRmNxe2elFF3TmbJJLDtGMhktnr7SGQTE+quWSH/aePpJWFT0lVIO8G7gTmAP8E9KqcdF5CMicn3wtP8JdAD/LCIPi8jtC+xu1RKP9PXZOzojfb99egZyfU9fQntH40DuSrF3tA7khp+vZfm2ms5IX4u9U+6vWib/xdJIzXKJtVONvYNS6g7gjoptH4zdfrXmdq04ZkVMSTtqlA2hfyBXy2SXMNK39QzkxiPhFTOQG7tiSZyyGfYty/YjfQ15+jrtnVn9VcuM7/hls4n0a8XMyG0UnoKMH/kmFdZ4qVrQNJCrM2WzFIi+hJF+QnsnNjlrpXj6OiP9aCDXEshmtXj6eu2div6aVPjDlOeMmBm5dWBEv0EoV0VLuyW9JC2nbGq0d3QO3DkOZDLllNKE9o4X+1FbKZ6+zuydaKA88vR12jv6MsPK/TVh/3fDFFXbRPp1YES/UXgKCdbJxUkY6VMZOembnKVnxqRbYe8knJG70rN3NOXph6Kvo+CaXk+/XGXT35Bwf+HykFnLZO/UgRH9BqFcVV4yMXGkX15+DjSXYdA4OUtEz0Duis/eSVx7J5iRG0b6Gjx9rXZfOJCraQa5csKqoibSrwcj+o3Ci4t+0lBnCSJ9nQN3oaevacxhJWbv6Iz0VTzSz+q2dzQO5Eb9IeE+XQW2ILbx9OvBiH6DUG7c3kkY6VRMztIR7WivvZPVl6e/ImvvxAuuJV4usZynj61J9JdkIFeXp+8htvgzvo3o14wR/UbhKSRjgaUve0fX5TLotnccJJPVlqfveA552y/TvFI8/fA48nZeQ6RfOZCrYXKW0mnvzLYjE9dccxXYFtjG068HI/oNQnkKLEFsK7mnXzEwlsbsHZ15+p7yaMm0ACvH0w/tnbyd11ZPH9vP09cyOUtn7R2vwt5J7On7kb5YYjz9OjCi3yjCS1JbkmfvVFTZ1BHt6LR3cBw/qtOUp+96Li22L/orJdIvuSVsscnZOX15+iJBpK9j4RP9ZRhEU5CiXOWX7baMp18PRvQbRRjpZ5JH+lTW09dQo1y3vUNWX56+oxxaM63AyhnILbpFcnYOW2xtyyViWaB7IFenp2/rmVeiXA/J+AO5JtKvHSP6DcKPToKMA02RfnkgN32LqEgmqy9P33MjgVxJ9k7GypCxMvrsHctCMnpm5Oq0d7RXhQ09fUuMp18HRvQbROjpY1uJa+/MqaefUk9fV56+q9zIClkx2TtekZzl/5Alrb2jPBdEYvaOhoHcpVpEheRjPJGnb0vyOlarECP6jcJViKUpt1jz5TIszeQsXXn6jnLIWBmyVnbFePqRvWPZeiL94LOWjA06Cq5pHcgNSyvbifcFQR/NGE+/XozoNwg1y9NPmrJZOTCWLnuHkoNky9k7SX1Xz/OwxKLFblk5kb5bJG/nscWetTJYXXheefwkjYuo6F7pLR7pG3unZqoqrWzQgKei7J3EC6PPOYk02ju6BnJte0nsnYK7MhZlK7iFaO6BFnsnivQ1VdlcSntHR/ZO1nj69WIi/QYRppmJDk/fq7B3NGTvhNUU9RVcy2rL03c8397J2/kVJ/ra7J3gs5ZMBqWj4NpSDOTqqr3jxvL0jb1TM0b0G4UbTs7S4UPqL8Ogc2H0aHJWGOkn/JFbqZF+zs6RkYyGKptudFUlmYweT38pUjYtTTNyHT97R8yM3Lowot8glOcXiUKHpx8NjC2Bp68jciqVkGyQvSOSuH2u52Jb9oqK9CNP37IT195R8YFcTXn6Wu0dzf01zNPHwkT6dWBEv1F4XpS9o6vgWppTNskEw0V28vK3YaSft/MrZiA3tHcssbQsjK59IDdcM1nrGrl6xqDiVqmJ9GvHiH4DUEqBh7YZuXOqFmpcLjH5CemCUr7NgF8aQNdAbj6z8iL9jGQSZ+/MGsjVVWVTq70T1ooKUzaTTs7y/Ktm4+nXhRH9RhDVygmyd3RV2VyKSD9p2wLBEbsc6etYGD20d1ZKpD/jzkR5+slLK1cM5OpcI1djwTVtkb4Ty4QzkX7NmJTNRhBVQfSrbCYtuFa5KIWOpbO0RXah6GeDrmVZWuydjGSwLIsZZyZZ+1JCPE9fi70TDuRmNVXZ1LmojublPX1P3wqqbJoZubViIv0GEEZLYgmSEX32js4ZubrsnVD0Ndo7juesuEg/zN7RkbIZt3fIZECpchG2OtFq71SmGCfdZaz2Dp6e/r+aMKLfAKJL5CDST1xwrWJRCp32ji7Rjw/k6lgY3RabnJWj4Gny9L9yIzzxb7W95uGvwm2/puXtC26BlkwLGckknpzl2zvlyVlAYotH6SzDMKe/asrTt/WU7l5tGNFvBKGnGdbT1zU5S9e0dpbA0w/ER0QgYdSp3dNXCvZ9Cw4/UNvrnv0R7L8r8du7novjOeTsnJ+9o2FhdJFQ9P0f26SLo0dBgIaCZkqjHak8PykiWi4RPeMOqwkj+g2gvHC1nto7oUeqcyBXaaq9E4pNKD46Uggrs3cSH28YWdf6A+IWa3/NPIQZSGGefuLsnZIDmfJALgAJK20uhb2jpb+GP0IZy0T6dWJEvxG45U4fFolKYnlU1tNPk70Tio0EImR1tONNTCTaZTxP31Ne8mwXJ7CI6hF9p5DYlA6vVsKUzaT2jjs5gd3e4d8JPndt9o7Wgdzk/TUaH7OtaA1mE+nXhhH9BqBi2Tuh95okv7iyyqbW7B1N9k7o6dudXbgT44n26XpuVHsHoOAk9PVDsXdqFH2nAKjylUKdhJG+roFcb3wCq6sLiNk7CURfeSrqUnpSNiuqwiawI8NihZIJrFIwkX6NGNFfhOnp6eQ7ccNOXx58SuKVVmbveBonZyW2dwL/PvT0rc4OvPFkkb6jnKj2DpB8glYU6de4n+jHItn7xyN9HSmb3vg4dqcf6ZcHcuvfpxcLInRE+mH/tLTYO7GkCI1VZrWc501CVaIvIteKyF4ROSAi75vn8ZeKyEMi4ojIL+pv5vJw7Ngxbr75Zk6ePJloP5FlG3j6kDCCCk6abN5fLNwpJs9oKQ/cafL0gzx9u6MTd3ws0T5dz8WyrCjSTzyYG4p9reJdry1UwSxPX5JH+u74OFZHJxCbH5HA04/3AR3ZO07R/7wyLS3BGyS4yg0j/TBlExLPyj169Cg333wzp0+fTrSfZmFR0RcRG/gscB2wA3iriOyoeNqzwNuBr+hu4HIyMjKCUoqxsWSiFUX6duySNEH9nTBSyrW2AVDUEKXEF1FJFIlFnr4vPlZnZ+JIP5ycpS3SD9fZrcfTr+d1FVQO5Cb19L3xcawo0k9u78Sjex1RdHFmmmy+Ra+nn4ldNSdsY3iej46OJtpPs1DNjNwrgANKqYMAInIrcAPwRPgEpdSh4LEVNT2uGEQopVKyTIi4py928ki/PNnFIpPPU5xJLvrxE0epclXkmvdTMTnL7uzEG6/f01dK+QO5QcomaLR36o30E9o7szz9hCtnKdfFm5zE7vQ9fXSIvqvX3ilNT5NtaYlmDSfJ048P5EbrSCRMK3WCzyo831c61dg7m4HDsfuDwbaaEZF3isguEdmV1DJpBGEncJJOa4/V3pGMBk8/TNkUi1xLKyUNkb4bO9ET/SBVDORanZ14U1N1i1AoiGH2Dmi0d9waf8zrvUKoIB7pZ6xMInvHm5wEmBvpJ8jTjwcAOuyd4sw0udbW8kpqSXYZnjd2PNJP1r4wqDOivwQopT6vlNqplNrZ39/fyLeuC22RfnjiWDF7J8nJ5JUHhnMtrdoj/STRXSRCgfUUDjDWm7YZCmI8e2fGTVh/J8zaqXkgV0+kXzmQm8Teccf8qyi70/f0rbZ2oPw91MNS2Du5ljY9kf68nr6J9GuhGtE/AmyN3d8SbFvxaI/04/aOBk8fhGyrHtHXdUnvDg8DkOnrBcAKbAe3TtEPBdGSpRjIrSdlk9p/LCqIR/pJZ+R6QTqsFYi+3dvjN3FkuP59LpG9E3qGieaoxDx9NHn6JtKfy4PABSKyXURywI3A7UvbrHSwFJG+juydyN6xhFxLixZ7Z1Z0l6BtTiD6dl8fULYdvDoHw0NB1JqyGdk0daZs1moLVVBp73jKq9vXD8dLwkg/0+v/2IY/vvWwdPZO0PeT+DtRymZ5clbSPH0T6VeglHKAdwN3AnuAf1JKPS4iHxGR6wFE5GdEZBB4M/B/ROTxpWx0owjFXqenH0UnSSL9cH+Bp1+cSV5uWJe9454ZRlpasFpbgbIYuXVm8IQ57HF7J3Gkn7KBXKDuaD/8XMOUTbunx2/imTN1t0+/vTNDrqVVi0iX7R19tXfC8zxpcNcsVFVPXyl1B3BHxbYPxm4/iG/7rCi0R/q2hdiB2CeakRuKvpBtbaN4YihR+6C8PJ5/O5m9YwfRJpTFyKtzVu7SRPoJyjDU87oKKlM2gQSRvn8FFU3Oymaxurpwh0fqbp+uACCkND1FrrUV0Ji9k7GQMHAykX5NmBm5Z2FJPP2MDk8/eK1IkL0zlax9zD65k8zw9UW/J7pvdwWR/lhC0bdsWmx/ck/ylM06xXsJZuRmxI+76p2VG0X6QRkG8H39JPaOp3lyVnFmmmws0tdScE1T+jMYT98QQ3uevlWuvaNjRq6I7+lryd5xFVYm+cCYMzJMprcvuh8OMNabqx/ZO1onZ9VRe0cp7ZF+WHsHqLuIXHgFZXd0RNsyPb3JRD+cB5JJvhyhUsq3d2Ipm4kmZznlPH00efpG9A0R2iL9eWrvJEkzi9s7uVbf009aadPzFJngKiSRvXNmtr0TilG9RddCMYxn72gT/Vr2Exd6DaKfkQwZK4MVDG7WH+mPI/k8kstF2+y+vmhAvR5Coc9krMSRvlPwq5LmWlpj2TtJ7J1ywbVyHStj79SCEf2zsDQzcpN31LDsQjbfQralFZTyT64EeJ7CzlrR7Xqp9PQlm0VaWuouxRCKoW2VPX1tA7lusfo6MHFLR8NAbngskb1T50CuNzYeXU2F2L0JI/2gb9pZK7GnH16FZltafeEHSkkSD2LjYybSrw8j+mdBW6QfXZLq8fSnx/waIW3d3eX6OwktHuWWRb/eS3pVLOJNTEQ5+iF2Z/1F16LJWeJHxVkrq28gF6qP2jVG+uGi6EDZ3qlzgpY7MR5lSIWEnn69V3+h0NtZK7G9UwzGm3KtrbQG4w5TY/XXuImXVhYNmXBgIn1DDF2RvjtRBAGrLatlRu7U2Ch2NhtET/7gZlLR9zxFJhsUxKrzHHKCjJF4pA/Jiq6FYhiKo5YlE+NefrVR+1JF+pYf6Ze8+vqYNz4xJ9LP9PaiCgXUVH0D/OFVaCZrJ7Z3wnTiXEsrrZ3dAEyNjdS9v3jtHavDLyPtTSQ7P02kb4jQJvpjRayOHGIJVj5I0Zup/+phemyUtq6eIGXTv2ROWmnT8xR2Qk8/nAVq98wW/cyaNThD9aWVnp7xy9325v195u28Pk8fqp9opdnTDyP98LiGZ+qzY5wTJ8isWTNrW/ij69SZthlF+pnk9k44cTDb0kImmyXf1s50gqq1asYJFiMSrHwGydm4Y8n6g4n0DQC4rht1hqT2jjdexO7yIzvJ2lhtGdzR+jvq1NhodKkc+aQJRT9u79R7orvRbNzZop8bGKD4zDN17fPU9CkA+tv8Wk36Rb/KfekUfadAPuOLfnhcJ6drL0CoPI/iM8+QGxiYtd0Osqfq9fXj9k7ySN/vl7kgOGnt6kpk77ijBezufJQJZHflcMeTfR8m0jcAs6N7HZG+3RnLrujO447W38GmRkdp6/IvlUPRT2LvhMvjZSJPvz5/J6q701sp+ufgDg/j1lGv/OSUL4b9rb445uwcM07Sgmt1WDU67R2vQM7y+0N4XKemTtW8H2doCFUokDvnnFnbk9bfKds7GgZyQ0+/xR97au3qZjqBveOMFrC7y+eS1ZlNLPphUFcqlbSsQpd2jOgvQPxXP2mk78YifQhFv37hmB6PiX5rctEPl8eLIv06o7vi4CAAmfXrZ20PI9F6ov2T0ydpy7TRlvVFo6+lL4r+62aZB3JPTZ2ir8WPxvta+rDE4sT0iZr3Uzx0CGBOpJ/dsAGAUvB91MqsgVxNnn42GHtq6+phKoG9444Wsbvz0X27K483pifSr7y9UjGivwCh6OdyuUQdQbkKb7KEFY/0e5KJvm/v+KKfbUnu6YcnduTp1xndFfbtJ7NhA3ZsdigkFP2pk6xrWxfd39K5hSMTCYu8LuNArlKKIxNH2NrpF661LbvuH7JI9LcPzNqe2bABq7OTmX376mqj52r09Cvsnbaurij7rFaUp3BHC2Tiot+px97JBfMcVoPFY0R/AcIvv62tLVGk700UQVER6efwphy8Yu252aWZGZxCgbbuHqB8MpWSRPqxy/n4/Vop7NtH/sIL5mzPbt0KlkXx6UM17/PU9CnWtq6N7m/p3MKJqRPJfP16ovZ6rg7mYbQwykRpgs0d5XWI+lv7IxurFoqHDiGtrWTWrZu1XUTIX3ABhX3762pjNDkrSNlMMvEvDEZCG7Ktu4epsdG6LERvsgSumh3pd+ZQRQ+vUN85qpTCcRza2oLUZyP6q5fwy29vb08U6bvBpWelpw/UFe2Hg2CVA7lJIv345TzUl6evSiUKBw/ScuGFcx6zcjmymzbVbe+EvjfAlo4tKBRHJ47WvK+IukQ/1gcSiP7ghG+5bOks1ydc27q2zkj/GXLnnFNekSpG/sILKOzbV5dg6+gPURtnpslkc1i2n7XW2tmN8jxmpmpf5CU8X+KibwXBlFunxRMGdO3t/uIzRvRXMfFI3/O8ugd4wktPfaI/4rcrsHcs2ybf1s6khkUzktg7haefhlKJ/EUXzft4bmCAwlNP1bRPpZQf6bfNjvQBBsfr86uBZR3IDdsdF/3+tv66sncKTz01x88PabnoIrzxcZxjx2rer47+EDI1MkxLzO5rCydojY7UvC93JBD9ntmRPiQXfRPpG2ZF+lD/AE8k+jF7JxOJfu0dLMxxbuvqibb1bdnK6SPP1tU+KEdySQZyQyshP0+kD9B62Qsp7N1bU02YydIk0870nEgfSObruwWwssHtGgdyrayeSL9jdqR/ZuZMTfV3SkeOUBocpO2yF877ePg91OPr6+gPIacGn2Xtlm3R/dbAlqzH1y9H+rEAKjivvDp9/fC8NpG+gUJQyybsDPX6+u5owZ+N2zFPpD+SxN7pjrat3bKNU4efTTztPpPgcn7m0UeQbJbc9u3zPt7x4heDUkz96EdV7zOMfuOe/trWteTtfLJI3y1BPpjFWmukn+9MJvrjg/S19EXZSADrWtfhKY8zM9UvfDLxwx8C0P6iF837eP7CC0GEmUcerbmNOvqDvx+XM4OHWRMT/fAKtZ4JWs5oETKC1Z6NtkWRfp0p0JX2TiFhDatmwIj+AowFnbInWImo3ki/dGSCTH9bubomIFkLuydP6XjtvuboiSEQoT1oF8DarecwMz5W1yUzzC6wBbVfziulGP/P79F29VVYsWqPcVouuQSrq4uJH/yg6v1WTswCf5ByS8eWKGKuC6dQFv1aB3LzncnsnYnBWVE+ENlXtVg8kz/4IZn168mdd968j9udnbRedhnj3/tezW1M2h9CRk8M4ZSKrN1ankfQ0efPHh45cbzm/ZWOT5JZ0zprDENabOzuHMUj9VVxDc/r3mBuyViCdNJmwYj+AoyOjtLZ2Uk+70fl9UT6ylMUnhkjP9A157H8ud0UDo7UHEUd2/8k/VvPIZtvibat2eKfVKcO1zfrVVVEdrVezhf27ad0+DCdr3r1gs8R26b9qquYvO/7KLc6G+PJM08CMNA1MGv7QPcAT555sv6sErdYR6Qf/Djku+qO9D3lsX94PwPdA7O2b+/yr472ntlb3X6KRSZ/9CPaX/SieQdxQzpf+UoKe/ZQHKzNCkvaH0LC/rhm6+xIv6t/Pcf2P1lbmxyP4tOj5M/tnrVdRMgNdFM4NFZXfwjP687OTnK5HKN1TCBsNozoL8Do6Cjd3d1ks/6lZD2Rfun4JGrGJbe9e85j+XN78CYdnBPVF8XyPJdj+/ey6aKLZ21fu80X/dOD9fn6lQN3tf4QjX/3OwB0vOLlZ31e1+tfj3PiBBNVRp+7h3azpWMLG9o3zNp+9carOTJxhKfHnq6pnRFOAXLBoiO1evr5jroj/SdOP8GZmTNcvenqWdu3d2+nN9/LrqFdVe1n7I478EZH6Xr968/6vM5XvRKAieD7qZak/SHk9GG/P8btHYBNFz6Ho3v31CTSxcFxVMmj5byeOY/lB7rwxoq4w7V/L+F5nc1m6e7uNqK/mglFP5MJqiDWIfrFZ/xLxfw580T65/k/BDNPjVS9v9ODhylOT7Hpwtmi39bdQ0tnFycOHay5jTDXw60lsvOKRUZu/UfarrySbEW+eCWdr3olmU0bOXPL3y+6X6UUu4d2c/n6y+c89tItLwXgvsH7zrqPgycnePjwyNwH3AT2Tq5j3tfsOnSGw2fO/gN+z+A9WGLxkk0vmbVdRLh8/eXsHtq9aDOUUgzf8vfkzjuP9hfP7+eH5AYGaLnkEs585SuoGq5Uw0y1pJH+yWeepqt/fZRWHLLpoouZHBlm7GT1s5ALT42CMG8AlRvwtxUO1S7Y4XmdyWSM6K9mlFJzIv167J2ZvcPY3Tns3vycxzK9LdhrWph54nTV+zu6dw/AHNEXEbY991IO7n4At04bCurzcEe/9q84J06w9jffuehzJZOh76b/wtQDDzBx3yKCPXqQkcLIvKK/sWMj5/ecz92H7z7rPv7sW0/y+7f+ZO4DdQ3kBkKfa59X9N/1Dw/xqbvOnilz7+C9XLr2UnpaeuY8dvn6yzkycYTjk2f3usfvvJOZJ56g71d/9azWTsiad/4GpWeeZew/7lz0uSFJ+kNIqTDD0w/vZtslz5/zWNh/j+7bU117lGJ6z2myG9qxY4O4Idn1bUhrhpl9tacth+e1ifRXOZOTk7ium8jecc7MMLP3DG2XrV/w5Oy4YgOFp0YpHq2u1vz+B35I55p+utdvmPPYxS95OdPjYzz76MM1tRPmDtxVeznvnD7Nyc98htbnP5+2q69e/AVA702/Qm77do5/+E9xJxY+7jCK37l+57yPv2bgNewa2hX5/vNx+MwUg8PTuJXHM2sgt9oqmwWw85DJz/mhmCo6nJooMDi8cKT/0NBDPHH6Ca4555p5H9+5wT/O+44s/GPoDA9z/OMfp2XHDnre9Maqmt356leTv+ACTnzyk7hVDlLqsHee2v0ApZlpLn7Jy+c81r9tgJaOTvb/+IdV7at4aIzS4ATtV87t9+AvQ9r+wnVMP3qq5pIMlfbO1NTUiq+/Y0R/HsJf+yT2zsQPj4JAx1UbF3xO+xUbkZzN+N2HF93fqcPP8MwjP+H511w374/I9hdeTkt7B4/d892a2gkxeycT1Pqv4nJeOQ7H/scH8MbH2fDRj1QVdYI/O3fjxz9GaWiII7/7e6h58qIdz+GrT36Vy9ZdxtaurfPu51cu/hU6c5385U/+csH3OjY6g+MpToxXVOV0i37EDtUvju4UfcG3c3Mi/SPD/mzooyPzV/9USvGZn3yGta1refNFb573ORf1XsSFvRfyf5/4v/N63d7MDIO/82680TE2fuyjSNAvF0Msi40f/xjOiRMc+9CHqip/4HkKEbDs+u2dJ+75Lh19a9i645I5j1m2zfNe9VoOPHi/n412FpRSjH3vMFZ7hrbL1i/4vParN4KrmLi/tslolfYOsOKjfSP68xAX/Z6eHizL4lgNMxtLxyeZ+NFR2p6/btaU8Uqs1gwdP7uZ6UdOMfXo2dP1Hvi328hkczzvVa+d93E7k+XSV1/Lvh/dV/Vlc8jEsC9WreFEl0UiO1UqcexP/h8m7rmHde9/37ylF85G22WXsfFPP8zkD3/I4Xe9a07J5W8f+jZHJ4/ytue+bcF9dOW6eMcl7+DewXu5bd9tcx6fKjqMTvsndCjKEU7BF287X2Okn5s30h8c8fc/NDYz96oCuOWJW9g9tJt3XfouWjOtcx4H36J7+3PfzlOjT82J9p0zZ3j2197B9EMPsekTf07Ljh3VtTmg9dJL6f/932P8W//B0Mc+tmj2lPIUYgutwcpU42dqK2V96KcP8fTDu3n+Na9DrPkl5gWveT0IPPiNr511X1O7hyjsG6bzZVuxcvaCz8v2t9GyYw0T9w7inKq+JMmxY8fIZrN0dHQY0V/NxEU/n8+zbds2Dhw4UNVr3ckSp7/6JFZrhu6fO3fR53e9civZLR0M37Z/wYGo/Q/8kD33fY/LXn9DNLllPq584y/RsWYtd/71Z5iZrH55wqP7RsjmbdZt8y0P9yxrjpaOHePZX/8NRr/+ddb+t3fT98u/XPX7xOl505vY+Gd/xuQDD3Lw+hsYu+MOlOdxavoUn3jwE1zYeyEv3/rys+7j7c99Oy/e9GI+/uOP882D35x9TLGoezAu+p4HXikWtdewclb0QzE70g/373iKUxOzfxBu23cbn9r9Ka455xrectFbzvoW1w5cy7bObXz0/o8yWvCLko1+8985+PPXM/PYY2z+9F/Qdd111bW3gjW//uv0/dqvMfyVr3L4t36L0omFB1HdksKyhHXbu7Bs4diBkarfZ2p0hO988XP0btzEzp97w4LP61rbz/OveR0//fa/c+inD837nJn9w4z821Pkz+2m4yWb531OnN4bzgPb4vStT1a1Mp1SigMHDrB9+3YT6a929uzZQ3d3N61BBcvzzz+foaGhRSduOKenOfWFR3FOT9N340XzDjpVIrbF2v+yA7szx6kvPsbE/cdmeaiHHvkJd/zVJ1l/7vm86M1nF9hcSyvX/tbvM3L8GF/78w9XXY/n6IERNp7XjZ216N3QNu9J7pw+zYlPfoqDr/85ph95hI1/9mf0/87vVLX/heh5wy8w8NWvYnd1ceQP3suBn/95vnjzTbgTE3ziZz+BJWfvnrZlc/PLbuYF/S/g/fe9nz/5/p9wbMK/Ijs2Whb6IyNx0Q9E3s5BJlfbQG4meM0C9g7A0eC9jk4c5Y/u+SP+9Ed/ytWbruZjL/7YohZY1s7yiZd+gsnhk9zy0bfy1PXXc/QP/5Ds+vUM/PM/03XttdW1dR5EhPV/9N/Z8KEPMnX/jzn4+p/j5Gf+ct6yGEcPjNC3qYNszmbdOV0c2TdS1XuMnz7Fbf/vB5kcHuba334PmQUm6oW89Ffezpot2/jGX/w5g08+Hm1XnmL8viOc+vLjZNa00vfLz0Gsxe1DuztP31supHR0klNffAxnkdpWp0+fZmRkhAsu8CvDdnV10dbWxp49tV0pNxvVGYOriCNHjvDss8/ymte8JjpJzz//fL7zne+wZ88errzyyjmv8Youk/cfY+y7z4LA2l99Li3n98553kLY3Xn6f/NSzvzTXka+foDJB4/TesVaHjnwPR74xm2s2bqNN/zxh7Azi/+InPO8F/Bzv/dH3PGX/4u//+Pf5Wd/+e0858Uvw17AA56eKHLm6CQXXuH7peddto7d3zrE1FiRlpzH5A9+wNi3/oPx73wHVSjQee1rWffe95LbsmXe/dVK6/MuYfvX/5Un/vkLHPnc57j+q0V+Ppcl/5NPc+aqq2n7mZ8hf965C3rYXbkuPn/N5/ncTz/Hlx/7Mt84+A1esfUV9HpXgZXBUq2zB1hDkc/k67B38uVIXykI+sfg8BSWgCczfOvpb/HlAz/m3sF7ERF++wW/zW887zeiBdDnwysWKezdy/TDP6Xr3nv5m/tdrNLTHNrcwpYP/j4DN/7GgjZJrfS+9a20XXkVJz71SU597nOc/vKX6XrNa+i67lrarr6aiTGXk8+O86I3ng/Apgt6ePiuZykVXLL5+e0V1ynx+D3f5b6v3oJbKnH9e//HnAyz+cjmW3jj+/+U2z72Af75Ix/g6l/4JXZsfynTPzhO6fgULRf30ffmC7HaFu/3Ia071rDmV57DmX/cy9CnH6LrVdtov2LDvNZQKO7nn+8fq23bXHHFFdx9992cPHmS/v7+Oa9ZCUiSWtlJ2Llzp9q1q7rJKI2iUChwyy23cOrUKd7znvfQEqz2o5Tib//2bzl+/Djvete76OvrQzkehadHmX70FFOPnkJNO7Rc1EvPDeeT6WtZ5J3mx3Ncjt31GIUfnSJXzFF0Z5jqnGTrqy+j/cJ+7L6WqiIe8HOk7/zrzzB0cD/tPb089+Wv5oIrXsS67ediWeUT4NG7B7n31n284T3PZ01ujOO7D/DN79hcMvV91v3kNiiVsLu76bz2Wvre9jby585fW6ceXM9l99Bubt17K3c9cxfrWvv5RO872PLjQ0z8539SOuqXT5Z8nvz555O/6CJy55xDdtMmsps2kt20iUx/f/SDcGziGP+49x/5l/3/wkhhBKWEnLuJDmsLN11+Bdu7t7POhbVf+SXWvOpjtDz4BVh7IfzyrYs39m9fj1ccY+Ki6xj9/v/izNtvZ7A4yuHxw9zy4G6mrWcoWccRUaxrXcd126/jph03RRPLvGIR98wZSkePUjp8mOKzhykNHqbw1EFm9u6FYEAxd845dLz85RzcuZH3nf4iw4URrh24ljde+EZ+Zv3PYFsL+9q1MrNvH2f+7u8Y//ZdeOPjSC7HkcveypMtV/CG18HaF5zPsZFWvvm5R3nV2y/mObGkBM91GTp4gH0//gGP3/NdpsdG2XTRDl7zm/+NNZvnH3yvRHkK58QUE/tPMPjdh+mY6iJr5SjmCrS+dAPrX7kDq84fu9LJKUa+foDCU6NI3qbtsnW0XrKW/DldSMZiaGiIv/mbv2FgYICbbropet3ExASf/vSn2bx5MzfddFOUvZcmRGS3Umr+tLZqXm9E3xf1gwcPcuedd3Ly5Ene8pa3cPHFF0ePqWmHU08P8cV/vQUbi5/tuJQtp7oQRyE5m5aL++i4eiP5gYX99kpcx2H0xBAjQ0c5+cwhhp7az5G9TzA1OoJYFs9/3mt4zrorkeMu3oQvCJKzyG5oJ7u+Hbs3j92dx+7Jk+nJY3XkkJw1y0JQnsdTDz3II3fdwaGfPoRSily+hf7eNXRncjDTylOF59FbmuCFD/01UiqigAd3vp+Z9n5ese0gm176fNqvuhLR0PmnnWn2De9jz+k9/PTkT/n+ke8zUhihM9fJjRfdyH+95L/SmeuMnl8cHGRq1y4Ke/dR2LuXmX37cE9V1J0Xwerqwu7pxu7uwe7pxuru5r7jQxwunULapxnNjDOTm6aYhUIWXMv/ExEyeGTausnm28hkcpw4/SwdHX205TvwHBfPdaBUoDh1hhmxEU9hKbA9yJegpQj5Ygu9bhfuaAvntq7lynXrUBOTOCPDuGeGcc+cwZucnNPuzIYN5M45h9ZLnkvLJc+j9XmXkN1c9q5HZkb40mNf4rZ9tzFeGmdd6zqu2nQVO9fv5Dl9z2Gge2DBgeFa8IpFf1D9nse59+h5dIw9y+U/+QsA3JZWdr/gt5nIZDiv9XFUS5HRUoGTp09SLBawLJtzL9vJpde8joHnX1bR/xTetIM3WcIdL+KNF3FOz+CcnsY5PUPp2ASq6I8fWR1ZvE3CvuMP8tPHvoPnOrR197Dl4ktYf+75rNmylTWbt9G1bt2soOVsKKUoPjPG5I+P+4kSjkLl4Jm+Ub4//giWbfPrb3ob3Zv7sNqzUdsfeeQRvva1r7Fx40Ze+9rXcs4CaxYsFw0RfRG5FvjfgA18QSn15xWP54FbgMuB08AvKaUOnW2fjRZ95XqokkdhcoaJ0XFGR0YZOj7E0MkhDh0/zNj0OJ25dl657Uq2Zdbhjhbxxgq4Y0VUye+YwzLBf+YeY1gm6ci2sX3zOQzsOI/evm662tuxlEdxepri9BSFqUkKU1MUp6eYmRhnanSEyZHh4G+E8dMnZ6XPdfevY8PWAbZdeDHbtp9PPpuDUgmvVMI5XcI57eAOu7hjCnfKAmeejq9clCqAO4MqTqIKE3gzk6jiFK5XZNJWjGVtJjMwLQUcVcT1Snh4eMoll8vS0tmB3dnFqaNFHNeiZ1MPa7b20buph7a+TjL5HJl8nkw+Tzafh4yNaysKUmLcm2TUnWC0OMpwYZjhmWGOTx7n2OQxjkwc4djkMTzlH3NvvpcXb34xL9vyMl629WVVi5c3NUXp2DFKR49ROnYUZ+gE7shI+W90FHdkhJHTI2RLRXLVpmQmRInFTCaHm8vTu6YLq72DTG8vdl8fdl8vmb4+7N4+/wplyxaymzcvWJyukhlnhrsH7+bbh77N7qHdUTVOQdjUsYltndvob+tnfdt6+tv66W/pp8Nup03ytOL/5VQG5bg4xQJuoYhTLOIUCkwNT3BmcJjTz55h5NgouZzH2g05nLExZiYmcYoOIhaWZLDFxiJDq8rR4Vp0l1zaXJuMnUfybUhLB5JvR+w8ZFoRaYF5xmUk62K1KuxOhd1rkVmbwe7NYufzSC7HTKnI0/v2cOTQUxw9uJ/xM+UJjJZt097TR0dvHx19a2jv7aO1s5N8W7v/195Ovq2DfHs7uZZWHKWYnJrmxNApjjx1mKcOP81EaYq1qpNXFC+hWwXVTjMWmR4/iLI6sjxdOMb3Bh9gsjRNd1snW9dvYcuGTfT1rqGnr4eeNT1kWnNI1p5VTLERLLnoi4gN7AOuAQaBB4G3KqWeiD3nt4FLlVLvEpEbgTcopX7pbPu9cP1m9dkbfwfEt0dBAEGJABb+Jv9/JVL+AzzKt5X4j3uAJ+AKuKJwgv9d8XBEUaDEtBQpyex0tbzKsMHrYcBdx3ZvHcotUXKnmXEmmHHGmQ7+xp1hRounGHPHKHb34nT14rR3gR3zal0HcRzEdRDPQVwXcV1QHuJ5ZFyXjOOScTxyJZe845EveuQdRcYDUQoQJDwwiG4LsW2AWFmsXCeZfBd2vgs704qVzWPbeexMC3Y2T8bOY0kGy7KwxSKDfyXgf9KxSoWxe+Gt+OPz4SkPhUIpD/xb/lURCoUXu+0/R/nfVvC/f2DR7XD7rLes7JcqtiX2OuXvFxW8g1L+D4vnYVsKWxTKcxClUAo8K4tnt4R2PJYzA+6ML05WlrzdT8mbwKMAygPP8f/PtODZef+dFVhuAUsV/W9DBBEbFws/pd3CssLPOngcgv8FYp/27P9V+Xb88FX5W4q/Ptoi1jy3y+8vWNF7n21wXM35zJnzmce/RU95OErhKBdXObieg+eWcL0CbrGA50zhlKZxSlO4zjRuaQqnMIF4RcDzP1f/bEaiPhC8jwieJSgBxxZmshYzOYtixqKUsXFswcnYOLaFsmywbJRto6zgz86gslmIXxW4DpnJcTKjp8iPjdKd6aE7u4a2TBdtmS5as920ZbrI2W1k7VawszxtD/GMdYoha4QZKWd6iYIcWfIqQ05lyGGRURa2Ev9KUAkWYCuwlGAp5WfNKOWrm1KIUtGxh/u0UNE3LMqfLxFuD7+HV3z6fYlEv5qB3CuAA0qpgwAicitwA/BE7Dk3AB8Obt8G/JWIiDrLL8qErfhB7/IvWFDC5bB1msPWae5joVF7C1gT/JWRytPEzvidbYG9LN3RFoI/Q3IWnyhnqId88Fd9goNuxMridfZR7OyjCIzjR7E+U8FfdfNxlECBEgVpvtm71Yj+ZmafCYNAZQpL9ByllCMio/gKOcuAFZF3Au8E2LJhM52l2ROXhHJkISoWyQSRrijBonx71r8wGo721FiqGRlRonwrRTw8UXgolKgocvJvR3GzfxUTRljxx6LDmxv9zgoQZYHt4Y302JSGZiK8WCP+v4RBesU1ZHArdpU6e3v8fnm7pcrntVXxvyDR45aKXzlVvv9yUj7jvOi8D8/F2ed3tDV2jhN/bJ6LvyQ0NGVTKfV54PPge/rv/fj7G/n2BoPB0PT80cc/kOj11eRDHQHiOVhbgm3zPkdEMkA3/oCuwWAwGFJENaL/IHCBiGwXkRxwI3B7xXNuB8JCKb8I/OfZ/HyDwWAwLA+L2juBR/9u4E78lM0vKaUeF5GPALuUUrcDXwT+XkQOAGfwfxgMBoPBkDKq8vSVUncAd1Rs+2Ds9gwwf81Yg8FgMKQGU3DNYDAYVhFG9A0Gg2EVYUTfYDAYVhFG9A0Gg2EVsWxVNkVkHNi7LG/eGNZSMSN5hbGSj28lHxuY42t2LlJKdS7+tPlZzkVU9iYpGpR2RGSXOb7mZCUfG5jja3ZEJFF5YmPvGAwGwyrCiL7BYDCsIpZT9D+/jO/dCMzxNS8r+djAHF+zk+j4lm0g12AwGAyNx9g7BoPBsIowom8wGAyriCUXfRF5s4g8LiKeiOyMbR8QkWkReTj4++vYY5eLyKMickBEPiNpWoq+goWOL3js/cEx7BWR18a2XxtsOyAi72t8q+tDRD4sIkdi39nrYo/Ne6zNRrN+N2dDRA4F59PDYbqfiPSJyF0isj/4f/nWMawREfmSiJwQkcdi2+Y9HvH5TPB9PiIily1fyxdngWPTe94ppZb0D7gYuAi4G9gZ2z4APLbAax4ArsJfKOxbwHVL3c4lOL4dwE/xFwbdDjyFX5raDm6fC+SC5+xY7uOo8lg/DPzhPNvnPdblbm8dx9e0380ix3UIWFux7WbgfcHt9wGfWO521nA8LwUui+vHQscDvC7QEAk05cfL3f46jk3rebfkkb5Sao9SquqZtyKyEehSSt2v/CO7BfiFpWpfUs5yfDcAtyqlCkqpp4ED+IvMRwvNK6WKQLjQfDOz0LE2Gyvxu1mIG4C/C27/HSk+xypRSt2Lv25HnIWO5wbgFuVzP9ATaEwqWeDYFqKu8265Pf3tIvITEblHRH422LaZ+CL1/u3NjW9aYuZbUH7zWbY3C+8OLpO/FLMEmv2YQlbKcVSigG+LyG4ReWewbb1S6lhw+ziwfnmapo2FjmelfKfazjstZRhE5DvAhnke+oBS6t8WeNkxYJtS6rSIXA58XUSeq6M9uqnz+JqSsx0r8P8BH8UXkY8CnwR+rXGtM9TJS5RSR0RkHXCXiDwZf1AppURkxeRur7TjQfN5p0X0lVKvruM1BaAQ3N4tIk8BF+Ivsr4l9tT5FmJvKPUcH2dfUH6xheaXjWqPVUT+BvhmcPdsx9pMrJTjmIVS6kjw/wkR+Vd8C2BIRDYqpY4FdseJZW1kchY6nqb/TpVSQ+FtHefdstk7ItIvInZw+1zgAuBgcIk2JiJXBVk7vwo0YzR9O3CjiORFZDv+8T1AdQvNp5IKL/QNQJhhsNCxNhtN+90shIi0i0hneBt4Df73djvwtuBpb6M5z7E4Cx3P7cCvBlk8VwGjMRuoKdB+3jVgNPoN+F5TARgC7gy2vwl4HHgYeAj4+dhrdgYH9hTwVwQzh9P4t9DxBY99IDiGvcQykPAzCvYFj31guY+hhmP9e+BR4JGgw21c7Fib7a9Zv5uzHM+5+BkePw3Otw8E29cA3wX2A98B+pa7rTUc01fx7eFScO69Y6Hjwc/a+WzwfT5KLMMujX8LHJvW886UYTAYDIZVxHJn7xgMBoOhgRjRNxgMhlWEEX2DwWBYRRjRNxgMhlWEEX2DwWBYRRjRNxgMhlWEEX2DwWBYRfz/1xncqxmjNuwAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Number of samplepoints\n",
    "N = 800\n",
    "# sample spacing\n",
    "T = 1.0 / N\n",
    "x = np.linspace(0.01, N*T, N)\n",
    "\n",
    "y_old = np.sin(100* 2.0*np.pi*x+1*np.sin(15* 2.0*np.pi*x))\n",
    "yf_old = fft(y_old)/(100*np.pi)\n",
    "xf = fftfreq(N, 1 / N)\n",
    "plt.plot(xf, np.abs(yf_old))\n",
    "#plt.xlim(-150, 150)\n",
    "plt.show()\n",
    "\n",
    "fm(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.7651976865579666\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEHCAIAAAACoPcnAAAAAXNSR0IB2cksfwAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeNrsXWdAFFfXPrO9L0tZei/SpKuIWFHsGqNJLNGoSSzE8iVRozGJJcWYqEneJPbYNTYiYgGNIiqCClKkCdI7S9vKsnW+H2tQcVl2V8A2zy+YmVvmzsyz55577nMQFEUBAwYMGF5N4LAhwIABA0ZhGDBgwIBRGAYMGDBgFIYBAwaMwjBgwIABozAMGDBg6DEQXtqeIQiCPR4MGN4QGB3dhXvJ7+pFISYmBmsda/3Naf3F3jg2kcSAAQPmC8OAAQMGzBeGAQOGXkZbW5tKpXqeGkgkkkQi6SXSIRDIZDJGYRgwYAAURSsqKtRqNR6Pf556HB0dm5qaeqfPCoWCQqHY2tq+7BQml8tzcnIUCsWAAQO0XpCYmNjU1DRy5Eg2m429ixgwGAGxWAwAzs7Oz1kPn883MTHptW4XFxfLZLJuscV6kMIGDRpEo9EEAkFmZuazZ5cvX15bW+vn57dx48akpCQmk4m9jhgwGAqVSkUikV65bpNIJIVC8bJTWGpq6oMHD6ZPn/7sqfr6+piYmNLSUhwOV19ff/To0UWLFmGvIwYMRkwkX8UIym7s9ovxheXm5vr5+eFwOADo379/UlIS9i6+4VCj6qTqpOSa5JzWnPL08hEOI3zNfbFheaVRUFCQlpbGZDJHjx7djf77l4LCBAIBnU7X/M1gMAQCQfup9evXb9iwQfP32bNnX+ADwFrvNVSpqs60nqlX1ZMRMhWhZmdn78ne40n0nESdxMKxsJHXASqV6uLiwufzn7/pbqmkHenp6WvXrh0xYkRNTc3q1auvXr1KpVLbz4rF4ry8PJlM1j0WXc8hPz/f39//2eNJSUkjR47U/L1nz56lS5dqDdjForTfhNbjS+ODDwdHno6ML41XqBQxMTEiuWh/zv5+R/qNODniYctDbOR1oLm5uba29vmbbmlp6eyURCI5ePBgcnLyl19+qX8P1Wp1+9+DBw9OSEh48mx5eblYLO6Wj71XrbDbt28HBwcTicTAwMCCgoKWlhYOh3P+/Pl58+ZhVvebicTKxNU3V/tb+P86/FcT8qMVMQaRMddnbrht+KJ/F310+aMjY4/YMe2wseoS/+bVFzWIjSsrlUqtTYUz+js8e0ooFC5duvTzzz8fOXLk8uXLLS0tQ0NDN2zYIJVKn7xs7ty5np6e7f+2u7qkUmlpaamrq+urN5GcM2fO7du3KysrfX19165dO2PGjDFjxjx48MDKyopGo/34449DhgyxsLCwsbGZMGEC9v69gSgXlq++udrL1Gv7yO00Aq3DWTcTt72Re2fHzV6SsOT4+OMUAgUbMd2Izao5l1VjdHEXC7pWCgMAExOTb775RvNR37lzJzQ0NDQ0VKFQPHmNqampZjaqmZDa29vj8Xi1Wj137tyoqCgHB4dXj8IOHTqkY7I9c+bMadOmtba29mY0CoaXB3KVfMX1FSQc6ZfhvzzLXxo4sZ1+Hvrzwn8Xbk3bujZ0LTZourF5qt93bxm5BiIQCEw5nX6JHA5H8weFQtFEoiUlJXWwwlxdXblcbnx8/OnTpwFg9+7dHA7nww8/dHJyWrNmTc/d9YuMzieRSK9iSAuGbsG+nH0Pmh/8GfGnJc1Sx2Wh1qGzvWcfyj00wmHEQJuB2LjpAI2EBzAyRh+VERhkA9hg2LBhHawwc3NzAJg+fbomjgpF0YULF3I4nM2bN/foXWMbjDC8AFSKKvdm7x3rPHaI3ZAuL14WuCyxMnHT3U3/TPqHgMPe2JcCERERui+4fv363r17vb29fX19AWDLli1jxozBKAzDa4ItaVuIOOKKkBV6Wet40qp+qz65+smxB8fmeM/BRq+XYWVl1b7BZsmSJfqbaWq1uhe6h4ntYOht5DTmXKu4Ntd3LpfG1bPIELsh4bbhe+7vkSgk2ABiwCgMw4vEHxl/mJBN3vd636BSywKXCWSC4w+OYwOIAaMwDC8M2Y3Zt2puzfOdRyfSDSroZeY1xG7IwdyDrcpWbBgxYBSG4cXgSN4RBonxbp93jSi7wH9Bi6wl5mEMNowYMArD8ALAa+VdLr88zX2aoSaYBn7mfgHcgKP5R9WoGhtMDBiFYeh1Eyz/iBpVT/ecbnQNMz1nVogqkmuSscF8+SGXy9PT0+/cudOjrWBBFRh6CTKVLLoweqTjSFuG8YrDoxxHcWnco/lHw23DsSF9yaFb9BSzwjC8YrhaflUoF07vM/15KiHgCO/2efdW9a1KUSU2pL2D1tbWQ4cOpaSkrF271iAtoNTU1F27dvV09zArDEMv4Z+if2wZtiFWIc9ZzxS3Kdszt8cWx34S8Ak2qr0A45Qqeg0YhWHoDdSIa1LrUj8J+ASB55Ub5tK4YTZhscWxi/0X4xBsGvGkobsBihOMK8pUqcDUCd47ovWs0UoVGIVheE1wpugMAExyndQttU1ynbTqxqrUutQB1gOwsX0MIh0oRuq+oEolkDsVyDVOqUJDahiFYXjlgQIaWxwbZhNmRbfqlgpHOIxgkphni89iFPYUhqyAISuMKyo2MAlbl0oV2EQSw+uD7IbsGnFNN7quyHjyGKcx50vOfxP6DSaF+ELQpVIFaBM9ffUoTEeyW7lcHh8fLxKJhg4dameHyQq/zogviyfhSSPsR3RjneNcxp0qPHWj+kakYyQ2wj0K45QqQJvoaU+gB72hy5cv3759e35+/pAhQ0Qi0ZOn1Gp1WFhYQkICj8cLDw/PycnBXpTXFWpUfansUrhtOIPE6MZqg7hBFlSLy2WXsRF+w9FTVpjuZLfV1dVNTU2//vqr5sorV65odNEwvH7I5GXyWnljnLpZ7g6H4CIcImKKY6RKKZVAxcb5jUVPWWEdkt1mZGQ8edbW1tbV1XX37t3x8fFJSUk9JOeI4SWZRVIIlKH2Q7u95kinyDZlW1I1lkcZs8J6ADqS3QIADocLCgo6d+6cpaUlg8F40lOGpcJ9nVpHAT0vOO9GcPv3wr/d3joKKAPH2H9rf2tm6xs78i9tKlzdeAVS4epOdnvjxo3Q0FDN39u2bVu+fDmWCve1bD2Tl+l7wPdiycUeav3blG/7H+kvVUixVLg9lwq3J9CNqXB7aiLZnuwWAM6fP69Zgr19+7YmlgSHw7W1tWm6LpVKNfNNDK8fEisTCTjCYNvBPVT/KMdRrcrWlNoUbKixiWQ3Q2uy2/ZUuGFhYX379u3Xr5+VlVV5efm5c+ewJ/G6UliwZXD3rkU+iWDLYCaJeaPqxnD74dhov0TPvfNoqleGwkBbstv2yTaCIIcOHRKJRK2trZaWltgjfy1RJaoq4hdN9Zjag68vjhBmE3a98jo6EH3+3ZcYugXLly+vra318/PbuHFjUlISk8l8VSkMukp2y2Qye/r2MLzIn+KqRADQJ1Pk82CY/bBLZZfymvJ8zHywMe8JtLa2nj592t3d/fz58/379588ebKOi3VHU716FIbhDZ9Fupm4OTAderSVwbaD8Qj+euV1jMKkSqlCrTCurEghIigIDKKWKb9BYjsdoqmSkno85AWjMAw9ArFcnF6fPsenxzPXsslsfwv/xKrEqICoN3zM1yWviyuNM7q4E9vp3FvavdL6i+3ojqbCKAzDK4OkmiSFWjHMblgvtDXUfuiv936tb623pL3RftWp7lP7W/U3erZoadLp6OkvtsPlcpuamjRHmpqauFxuT981RmEYegQ3q25yKBw/C7/eoDC7ob/c++Vm1c1pHtPe5DEfYD3AaPUhfjeJ7bRHU3E4nPPnz8+bNw+jMAyvHlBAU2pSwmzCekdV1dXE1ZpunVKb8oZTWG+iM7EdrdFUGIVheMVQ2FzYIG0YaDOw11ocaDPwSvkVNarGpKi7HYaK7TwbTdWjwJ43hu5Hck0yAshA616lMKFcmNOEqTa9FCCRSL3DXxiFYegRpNSmuHHcuDRur7WombSmVGM7jd44YBSGoZshU8nS69PDbMJ6s1EWieVt5p1ci2X5xigMA4bnQ1p9mkwl601HWLshlsXLEsvF2CPAKAwDhueYRdakkPHkYG5w71OYClXdrbuLPQKMwjBgMB7JNclBlkG9n1gogBvAIDJSajB3GEZhGDAYiwZpQ1FLUW+uRbYDj+D7WfVLrsHcYS8Ycrk8PT39zp07vdMcFheGoTuRWpuKAhpqE/pCWg+1Dr1Wea1GXGPDsMGexYvCoEGDaDSaQCBojybDrDAMrwzu1t1lk9l9OH1eSOv9rfsDQGp9KvYguhGtra2HDh1KSUlZu3atPtL+qampu3bt6rXuvbBUuABQVFSUlJRkbm4eERFBpWJ5tF4TCgu2DH5RIfKuJq5mVLPUutTJrpPfwMGXpNxWVFYYzVM4S0vW2LHPnjJIbKf377oHKUy3eOPJkyc3bNgwc+bM3Nxce3t7f39/7Pt/1VHfWl8pqpzlNetFdQABJNgy+G7tG7ooyT99SnjhotHFpc7OWikMDBHbeX0oTLd4o1wuX7Nmzc2bN21sMJ/Fa2SC1d4FgH5W/V5gH/pZ9btcdrlaXG3LsH3Txt/6u++s1q0zrqxAIDD5T1HnWegvtvP6UJhu8cb8/HwzM7Nt27ZlZGQMHz589erVBAK2sPA6zCJNyCZuJm4vsA8away7tXenuE9508YfR6WCsQ4ZnFqNYxiQpaUzsZ3Xh8J0izc2NDSkpaVt3Ljxu+++mzVr1p9//rl8+XLNKSwV7qvb+nXhdVu87bnYcy/23lk4VnRqNC4P9yaMfC+kwhUKhSqVSnNWKpW2tbXx+fzg4ODOali0aFFaWlpVVZWXl9eKFSumTtWS/+WVT4WblZVlY2Oj+Ts6Onrq1KlYKtxXo/XWFvTBRTTtAJp3FhXVP3mmRlzje8D3WP6xF37vK6+vHHFyRIeDspISwblzLdHR4pQUtUyGpcJ9bVLh9pQVplW88fbt28HBwUQi0cvLi0QiNTY2mpub5+fnOzo6YrOwlx2CKkj4FrJPg1r56AiCA48xMHI9WPR5SRxh7XPJuNK4cmG5I8sRACRJSbxtv7Tl5T2eNDGZpu+/b7bgYxy2Dv7q48WkwiUSib///ntERISDg0N9fT2WCvdlR14sxH4CaiUMWATek4BlBxIeFMTB3V2weyiM2QzBH6TWpZpSTF1NXF94ZzU0mlqX6kC3q//++5Zjf5McHKy++YbWvx+OzpAVFgpiYhp37RJeuGD3x+9kDw/s8WIUph06UuECwIQJE8aOHdvY2Iilwn3ZkboXLq4E+/7w9h4w+S+jGtsWbAKh/8cQsxjOLQNxfWpjaj+rfi9DPlpHlqMlzTKt6nbo79dF/14xmz/f4tP/Q4hEzVmitRVj6BBO2r3qFSvK33/fbscOWnAw9pBfXfRsCKJu8UY8Ho/x18uOrL/hwgroMw7mnH3MX+2gW8CMExAwq+rm5hpxzcswi3w0l7Ts570zQXTlquVXa7mrVrbz1+NZQkiw099/Eyy4lQsXyQoLseeMURiG1xHlyRC7DFyGwDsHoDPlCRweJv2e6hQCAP2UyEvS8dFJrf3uS0lR80zff7+za4jWVg779+GZjMqFC5WNjdjTxigMw+sFaQtEfwgcR3j3MOCJOl8i/D2HQDMUcb60DlqbX3zHMzO5R68meePuj+7CMUfgcu137lS18GtWrwEUxZ45RmEYXiOcWwaSRpi6DyjsLq+915gVZNUPkTTChU9fbK/VbW01q9cQraxOv2WW3pDR5fXkPn24a1ZLkpKaDx3Gnnm3oKCg4OjRo7Gxsd0T9oVRGAZj8OAC5MXCiK/Auutcto3SxipRVaDDMBi2GnJjoOjKC+x44/bt8ooKm00/eNoHpten61OE8957jOHDG375RVFdjT3550RqaupHH31UUlJy4cKFkJCQ1tZWjMIw9DbwajlcWgNcbwiN0uf6e/X3ACDIMgjCloK5B1xcCcq2F9JzRVVV88FD7MmTaf37B1sGV4gqeK08fQpaffMN4HB1GzZiT/9ZGCS2ExIScvPmza+//nrXrl0cDqcXhA+xnYkYOsKjPhZaKmB+XBcusP+QzkunEqiepp6A4GHsT3D4LUjdCwOX9H7PeVu2Ijgc9/+WA0AQNwgAMnmZkU6RXRYkWltZLF1Sv/knupvrq/vg0uLKagqN3GmkVCrZ5vSID7yePWWQ2A6CPFrSkUqlpaWlrq49Pp4YhWHo8MLWuPLiwf89cNBXPDqjPsPfwh+P4AEAXIeD63BI+hWC5wGJ3psdl2ZmCi9dsli6lGBpCQBeZl50Iv0e754+FAYAnNmzW46fMI+LhxUrAIfNTp6CoWI7arV67ty5UVFRDg4OGIVh6E60CNpuJFXxalsVCiWDTfbyNg/xs3xKoPDGzwiqhmFr9KxQrBAXtBQs8n+spAQR62D3cLizEwZ/3pu31vjndoKpqen8eY+mwwi+r3lfPd1hAIAQCBbLl8s/+0x48SJrwoT2420y5a2UmvIyQZtUQaESHV3ZgwfYkUgvI8eFjHWCsUaW5fP5OkI4DRLbQVH0ww8/dHJyWrNmTS/cNUZhbwruptcnnHlIb5C3f3wigLs36hMJuQxf05mzvNhMMrSUQcbhcvNhzhwnPavNashSo2rNrO0RbAKhzxhI/h/0XwBkZu/cXVtOjvjmTe7nn+Moj+PXgi2Dd2TtEMlFTJJe3WCNHVOydSvvt/8xx4xBCISiMv7Z4wW4cgnpv3ALKcCD2w33jxWhLoxpM70dbBlv7OvUmdgOiqILFy7kcDibN2/unZ5gFPb6Q9yq2PVHOqlEQsaBwpXhP8DKx9vMhEkurxSmptYKMhrVmc27sm+5j3OYotgGOEKh5SRnvSvPqM/AI/i+Fn2fOjp0NeweCvcOQNjSXjLBduzEs9mcmTOePBhkGaRG1VkNWeG24fpZYkhjZCR5//6Ws7HHG5wUWXwSgMSS5BRiFRxoxTWnNDS1pWXUFdytoxeLo7+7yw61mP9B3zfzpYqIiNB6/Pr163v37vX29vb19QWALVu2jBkzBqMwDMajpEJ4fNs9VhsqdaR9FBXAYT82Uvq4mfZxM4UZcON2za3jBTXnKm7SKQMj5rUpTPSvP52X7mXmRSPQnjpqEwDOg+HOTghdDLgef8dkRcWihASLJUtw9Ke8b/4W/kQc8V79PX0pDEDSxwO8/M//04wyTSSmxOmL+ro4PB4NexuGvY0bjHfLe9gcsyeblNLwY+mtZV8MoFFf5+/IysqqPRfRkiVLurTO1Gp1b3YPc1u+zsh72Bz9UxpVhjpMdly5JvRJ/noSQ0JtPv9pqIXFw/uSCb/eGKFQ6lu/Uq3Mbsx+ahbZjoFLQFAFeb0hH9hy5AhCJJpMf6/DcTKe7GXmpb87DACahLjLlh9I6A4WdqIvfxj8JH89CW930zU/DgU/E0ad7Nevb4kkcuxle1HAKOy1RVEZ//yvmTgUHbzIZ9LYLta2SWrxu7RvfR2SKHx8bhJdz28ytym3TdkWyA3Ucs59NJh7QPLvPX2bapFIEBvLnjCBYGb27Nlgy+CcphyZSq8w8YelfN4dugohBVT/3T/raBeTThx8EhVkPtKGKVb9b32KtE2JvXIYhWHoNghEslO/ZhBQdOgC3xB/PeRA7u0DmWjovLc4w61Npbg/vrutVHY9HdAYOAHcAG2fOAKhi6EmA6rTe/RO+f+cUbe2cjrZzh3EDZKr5LmNuV3WU1Mn+WdrOglFwj709pw9WpqdLdUjk+v0aZ7mETZskeq3zXewt+6FQNccPiUlJTo6ur6+nsFgSKVSS0vLyZMnh4WFYaP2kkOlVv+56Q6rDXV9x0Uv/lIr4c4ucBkG1n6z3oNvi8pNK+GXn+6s/LKL0LAMXoYjy9Gc2kneh77vwuWv4N4BsA3qqVtF0ZZjx6hBQRRvL+0UZhmEQ3D3ePeCLHX1QdyqOPDTXboKcP6S0BBrtfck3tZtzUeO2gYEdNmFGe94/lHfyszh79iZuXjRo+t56YUNNwqBp6QqqXiEpAaVHNoULCXNh+swth+BTMTe0p61wmQy2bZt2wQCwaZNmw4fPrxjx44DBw78+OOPKIoeOXKkrU3f7SOJiYnR0dEdcn88icLCwnv37mGPoXuxd28Wu1lJGWA+NkK/pcWCiyCsgQGPYrv8fGVydwatQrr/cI4u9gA0syEz2LJzvUAyA/pOg5zTIBP10J22pqXJy8tNn16IfBIsEsuF7dKlO+yPzXdYraj9OHsHazUA4Gg0k7eniC5d0lOE55OoIKEFSZXZfP1mZWlsSv7qWPnJemYtE0FxrWypmCuSmEkUZAVNwCDdVlV8nZC/+5JSKsNe1B6ksMbGxqVLl44ZM4b4hFYcgiCDBg2aMWNGQ0ODPlUvX758+/bt+fn5Q4YMEYm0vMRNTU3jx4///PPPscfQjbibXi9LbxGYET+a56dvmbR9wLYF98dR7Ms+DRGYEITJvNvpdZ0VKhWUtrS1aHeEPbaC5oJcAtmnemoWGR2NYzIZI0fquCbEMiSTl6lGO50X7z+UzayXI34mb090bz/ImTULVan4J07q0w0EB1Er+uHxIsY/d4jJSjxKEPdpM//cz3vTJN+1b/l8Psl31WTfjVOcfxiuHEqWktqYJbSi9VdqUnJeyzcwNTX10KFDp06d6pa8SkZSmK2tbVxc3GeffZaZmSkQCBITExv/+znC4/H29vZd1qtJhXv8+PGvvvpqyJAhR49qcY5++umnS5cuxUinG9EmUyYczJPjYf6nes/dmkugOBGC5gIO334Mj8PN+yxYhoMb+/OaW9p0OMK6oDDbILD2h3sHeuJm1WKx6NJl9vjxT4azPotAbqBmC4F29s6qF6Y0CNiERYuemjOSHBwYgwa1nDyJKvXy09fFJo1jsEyJJqnKarfvIj3njaJxO2aWRQh4p7H9fb6brBpGwQNOGdP0YP+/r5sTQ6X67bff6urqUlNT/f396+rqXgyFoShqamq6devW1NTU7du3V1dX//jjj1euGCCi0iEVbkZGR+Wmc+fOWVlZBQUFYbzTvVNIpgztM9HJylzv/Ylp+wGHh6DZHQ5bc+mhsz2pCtj9273OHGFmVDNNliBdCJgJtVnAy+v2mxXGxaulUvbbXaS8DbYKBoCM+gytLrArf+UqcDD3s2D8M/siObNmKuvrxdev665fKZVl/xBDzyWJCYJLprhqMfdqcheiPY5j+tmvGSwitzAKKDk/nQX1S622aJBSBR6PP3LkyKpVq3766acBAwZ0yIHdE9DuzkcQRCKRIAji6elJJpP79+8/a9asY8eO6V+v7lS4fD7/559/jo+PT0/v6KTAUuEa3XojH6fKpjXR1Xay7LNns/X6BVMrRufub2QFpibc1dp6oyWFWwebt130dFV0KHtTeNMWb9tlP8lK+mgEXxT9fZ7Nu9177/Z79+ItLS+XlkJpqe4rTXAm5zPO0ws70vq9uxSunCjzkabevqKldbXahcUq2L6jWizudADlKrtUlQXOsoRYwg9im6nFRZW01FMPBY2ZxK7CXZFgnOndUqdm5+S1Rxv6MwCHwPOlwlXKZSql8bEdDW1SIkVLVjoej7dkyZIlS5YMHDhw9erVdDo9JCRk8+bNHXziM2fOdHd/NBNPSEiorq4uKSnx9/fXOp3sxlS4nQ5z//79L168aGlp6ePjozni5ubW1tZG0Wm0t4PL5TY1NbX7vDT7P9uxfft2Mpm8cePGqqqqkpKS7du3R0VFtVPY+vXrNTQ6efLkF8ggr1zrP6y+QUUU8/9vgJM9S29rOQayRDYT1052Ha619bFj1VtXXScVU/q/P9Sa+5gCeK28llMtHwd/PNlbj362xbrzMt0nHQYE1133Li8rKy4vt/ziC389Lk6+mZxWl9ah2pu3qy2aCqT21JVLR3TWOq+0tPmvfeMHDSKYa1l1lTbwy7beYODM2gLQIdM/0ByMY5WWRJfW1dkt+USPxEiTIfeXcw71jqyHzb6rJhvx3FtaWmQymWaH9oXffnqQfMPot87Uxm7eLzu1PL22Ng6H88MPPwBAQUFBbm7uyJEjn90j6ejo2L5RPDMzs6ysjMvlmpmZad09LhQKBw8eTKfTe5DCOBzOuHHjamtrb9++jSAIj8draGjo37+/nvXqToX7zjvvaKaQubm59+/fx6aTz4/YuGI2X0kMNjOAvwAg6xiwbMB5SGfnSSTcqHneyTtyD/6ZuXrDoMeOMF46AOhajnwSfu9B9IdQngxO4d11v8LzFwCHY40fp9fbyA28WHKxWlxty7B9xD5tylvHCnEEWLhE17tn8vbbTXv2Cs6eNfvww4781SIs23qTjrJU4WS3iaHtx8eOct6UWEnJ4VfViu2su94H7vPpxJyfz5o0mT7Ycxm4zzUmnoOGWji5GFdWKpVyLLg62EDzhz5KFQCwceNGAFixYsXu3btXr179AiaSjx0i1tbW1taavysqKhITE8lk8sCBXStJ6U6F6+7urrE5GQzGuXPnQkNDMQ56HiiV6rz4CoQEn37gZ0AxSQMUJ0DY0icd+c8ixN8y2buSmSe8dK189PBHnq/0+nQagdbHtI9+39Z4oLDg/snupLD4OFq/fgSuXh+9ZgtUen16O4Xt3pXBkKOuU53ZbLKOgiQnJ1pQEP90tNn8+YA8zs8kF0tLf77GQNnocLrTmI6p5ybP9b26LfPvfdkr1+qluebz+aTc9THsYjOyuPZ5xsQ1ZIBryADjyuoW23kWnSlVyGQyEomkET7k8XjOzs49/fITOrNO23m3HQ4ODhoBM61nn4XuVLgahIeHJyYmYhz0nDh+6gFThnLH2BqmY3X/JKgU4D+jyws/+jjg95U3Mv4pHjLQlkohAEAGL8Of+5/MYZcgUsFzAuTFwLifgUB+/vuVFRTIioqtZs/W83o3Ezc2mZ3By5joOhEA8h42I/kisRVp7KiuPzD221Nq134lzcqi/hfmqpIrCn+8xFabKgbgXMZoSZ3p7W560YlKLZXyGqsSAAAgAElEQVRm5TX4e1t02QSCQzxWRhZ/n+BRayMoq2U7Wb/8r1xnShX5+fkzZsxwdHSsqalxdnaeP39+T/dE+xuPIMjBgwcrKys7HK+srDxw4ID+tetOhYuhWyBuVdQk1QloyLRJfQwrmXUM7ELA3KNrm5pK8J3szFTAX3vvA4BYIS5sKewinKIDfKdBmwBKrnWPCRYXj+DxrFGj9H3LEZy/hb9m8gsAMX/lqBCYtSBAn7KssWNxNBr/n3/aj+RtuWiiNG31Vri83alROWuenwIHF47k6/ulMOncD/wQBFe95w6qVL1UL1gHpQrdE8OAgICsrKw9e/akpKScPXuWSqW+GCvMxMRkzpw5J06cSE9Pb2hooFKpUqmUy+X27dv3gw8+QBAEI45umP01NgpiY8UJCfLSMlSpJNra0vr3M5k6lezublA9R4/k0FXg+ZYrYtCG17ocqMuB8Vv1vHz8KOdNN6poOS0PS/g8Sm5HmcMu4TwEqBzIiwWPblCPEl68SAsNxT+hdNwlArgBN6tu8mX8a1ca2XwlMdjU3kYvwUIcjcaMjBTFx1t99RVCIj346zJHaCqw5PvMmajLA8OlE33YxGzB7bTa0BC9rCozL8d0k2teAte8bcd9PJKg8i60NgKRBlZ9wXMC9H0HiNRX4sUmkUj6hI72uC8MQZDp06dPnz4d45puByqXN+7a1bxvv1oqpfbtyxw1ChBEXlHRcuzv5kOH2ZMm4gL0NXBapUphVrOCRYgYYqBIedYxIJDBd6r+JWYs8D/7Q9rp/dnMqVlaZA51A0+EPuOg4AKoFHpmFekMbbm58ooKs4ULDCoVxA1CAb1dmV4YB2oSfD7PAKche+IEQUyM+PqNZpRNKyTzKY0+y9/qstSs9332r0m+Fv1QTwoDAJxLW8u9HJNmL15GMdfHG5g20MaHqlSIXQqJm2DkevB7D/t89KWwZ6FWq/Py8jRijBiMhqK6umr5/7Xl5rLGjbVYtozk+Dg0VCUUNv+1r+mvv5wSE9v8/Drbuvwkjh3Po6kQj/FOhnVCrYTsU+AxGqgc/Qs52bNw3ixmrjAzs8zT3LOjzGGX8J4EmUeh9Aa4RTyfCRaHEIlMvWeRGvQ170vGkxNjSu2Uzn1muJIIBpis9NBQgplZ7bkrROpIGSL1+DwSwXU9EeGwKSRvE1IO/3Z6XWiQVdfNZJ8aVvC1nGNX07JJIltiPmUcrj20rCwJrqyHfxZA0RWY8BuQaNh31IUvrOOUR6lcv379F198cfjw4V6WZHzNICsqKps+Q1FRYffnH7Zbtz7JXwCAZ7EsPv0/pxPHUQRX/sEHralpXdgjMmXLvUYBAz9qqKNh/Si7BWKeET/pc+f7teLBLscvwCLQ4Jt3GQ4U1vOLIAovX6YPCsOzWAaVIuFJ3tRAy1JbkTnJ4OHC4xljxyuQQBRQ64+DSSx9o5lmzvKW49CE04VdX3pnJ0R/3EJ3pS6NRwYyGMB58NcTe4+cwuHDyzDia8iJhsNvQdvjQHEcDvcqfpIoinaXP4qg3xPER0ZGhoWFKZVKHJafymj+Ki4unz0HIZGcjv9Ncuk0fofi41O5eJHXyVOVixc7HjlM8fTs7MoTpx/QlOAUabjfIf8skOjgarA1xKST6AOBnmQnyWOBocv3BDJ4jIGCC6DeZrQataywUFFZab5ggRFl7XIHE1D82NluRpRtRF04VGuRWSnb1QDrz5RDIXqxSbnCjNyGQJ/OlyYzj0HcavCemEKcPJFh6TLFMiczhlnKEhTXsF1t/vPs4GDICuB6wun5cGgyfHAeyAwAoFKpPB5PKBQSicTnIQWFQqG/As1zoq2tTS6X6xkk3z0UplarDx48eOrUKU9PzwULFryx7nyhXNgkbVKhKnumPRlvWHCAsqmpcuEiBI93OnKY2JWzU8lmO+z7q+y99yoXLnI+fYpgYaHNNFbXpfBUNNz4SAOjGdUqyD8PHqONcw+zAmseZsmZaQ6CaXI2k2TgXHIy3D8J5ck6gml1Q3T1KuBwjBHDDS2YlddgU2ObY30jzMQUwMqgsuXxdzkSa15LhnljDsAHBpV9Z4bXia/vXI4u7JTCypLg3DJwHQZT96nPX9Qcs583oGXHg4pDd/pueHoHqOcEeO8IHJ8Jp+bAjBOAJ5JIJBsbGz6fr1AoULSLvZYKtUKhUiAIQsFTOnzFPB6Py+X2zndEIpHs7Oy6yxjqlMKeDHXD4/Fz584NDQ2tqKh4nfhLIWlTtsmpZqwumSu6MPpcyblifrFGswWH4Pws/Ca6TJziPoWI69o5jSoUVVGfqJqaHA8fIuq3WEPgcu137y57972alasc9v31bHLWsxeKGUqwHGVr8G1XpIC4HrzfMm7QMhszyt2bh6XPPrj//rJlIYYVdo0AEh3yzhpNYeKEa1R/f60a07px8XA+BQfptpfTeXb6birQkGZ1g/IaX4rIbDyBvzdF2dRkUOtW5nS5A41a3lpRLdaStE3SANEfgqkLvHv4yVUOtpN1lX06p8q8KjHTbtjTwR/ukTDhFzi7FK6sg9E/AACdTte9U6elreV4wfG40rhSwaPNpDgE18+q3yTXSeNdxmuC++7fvz9gwIBX8SvulAi/++67GzdutO/DHDhwIIIgjo6OrzxtSduKjl/P/SameNXl+m9Tm37OqvgiMX9NbO5v5wUlNR2pB9DThafH/TNu271tLBLrk4BPtg7dunnw5oV+C6VK6be3v50UM+lWza0uG+Vt+0V6/77NT5sphiyGkN3drdZ9I7l9u3H37mfPFtysERPgrfGGJ3zPjwUiFdxGGUlhvEx3LyuxBQnNE5RUGCgIRaSC6wgojAPUGG0GZX29NCeHGWHw/PfqjQpWi5Lsx7G1sMrgZRjgslGjFdtv4YFoPsuHM3ECqlIJ4+IMbX3SOx44gH9O5j/rEIIzi6BNANP2P5tw02NehFQtFsVXagkTC5wNoYsgZTs8uNCFwY2q/37w95joMTuzdlrTrT8L/mzr0K2bBm+a7zuf18pbm7T2rbNv3a27+0p/0Z1aYQQCITo6esWKFQwGY+jQoWPHjtV/g+RL6kFUqgqPXCPkqykIFUXJUpZEZqpECHiVSIZrQpk1TMGuh5Wsuy4fDqZZmQGARCFZm7T2asXV/lb9v+j/hQfnqRDQqIColJqUn1J/iroS9VHfj5YELkFAu30qvnmz+cABzsyZhi6iAQD7rbfEN5Mat+9gRUY+6T5LSatli1XEfqYEgoHWOKqG/HPgPsq4Va0qURWvlRfIDfT/wPvSlszTB/NWfW2gEHmfsZB/DupzwMrgDIyihARAUWbECEPv+E5MCZEAi2b7Nt4PjC+LV6NqnH5BdPk74tgqs9a+SvO+LgBA8fQUxcWbdqLT3+kdu5n+Y0EkPxQKRDI28wn/w/3jUHQFxm0BSx8tbE+n4Acwqalo4dHEPh88w9qjNkLlXTgbBda3gG2ntd1WZeuq66uuV10Ptw3/ot8XTmynJ88uD1p+rfLa1rStH1/++KO+H9mD/etGYfPnz/fw8ACAxsbG69ev7969+8KFC+0yOK+EuQUFcVB0BZqKQCYUKGyqq99l4WyF+GbSCDP3EeHw9Lq4oKy28tRdVqNJ3bYM3BCGBJXMjptdKihd3X/1TK+ZWulpoM3A4xOOb7qzaff93bxW3oawDc9+GCqhqPbLtZQ+fSxXrTTuPiy/XCNJSqpdt87x0KH2bXrXz5WQEZg5zdNwEkoFYQ14GSnCoTFhArmB7hzTC64MarE4LateL3n+J+dBCA4K4oyhsKsJJGdnkoHb7k7884DdqjYZbs2kkwK5gacKTxXxizr8IGlFTVI2o4LOpzX6vv/IIcUcM7rhf78reTyCgW6j8PEu9w8UnDxR8PFH/8WjtTbD5a/AfgD0+7CzUi5Tw/PunSXlkuViKYnxtOMST4Jp+2FnGJxbBu//82xZgUzw8eWPC1sKvxzw5QxP7XvIhtsPH2g98Me7P+6+v9uP5DdePV4fr4hKIBDGx7feTVVUVKBqNY7BoPTpw4wYQQsJATy+9z/0Tn+LNPwFAObm5lOnTt27d29ERES7fs5LDWUb3NgCv/jA6Xnw8BIQyFXyQU3Vi2hgLlX+5W26wp51D6Cjcc52svZdOZk+20GGkxKS5Kq04ipR1c6RO2d5zerMvAIAMp68Pmx9VEBUTFHM+pT1KHScH/G2bFE1N1tv2oSQjdwbSDAz465a2Zqa1r7HpaxSSK+XqZ1oneWF1IW8s0CgGB0in8HLYJKYriauADD7w74KHFw+9sCwKugWYBcChQZPx9RiceudO4aaYNI2ZcX1WiEFmf5OHwDQJAFp32mkq2CTUBJb1QYSt+WPbWfW6NGgVov+NVhqdXCorYCOa77f+Fj++t9voE0IE3/VLUDEHutMwdGLDidqOcdxhJEboOgqZHaURBYrxIuuLCoRlPwR8Udn/KUBhUBZH7b+s+DPsuXZq2+s1iHPrSGv+s2bHw4ZWrduvTQ9HW9mSrK3R+Vy/qlT5R/MLR4/QXjhgnEugh6hsA7YsmXLgwcPejmi4tdtd3fszDx7oaiWJ9G3TPU92DEIEr4D+wEw9wKseFjC+kJdO1GByFnz3d0XzQauF1xcCXtGQLMWnTwzXxfbteEPIG+4ctgfZZ8NsNRr7rzYf/Fi/8VnHp7ZkbnjKUv+3j3+qVOmc+dSvDyfZxxM3n6bGhjY8Nv/1FIpAJw59QAAJkz1MLgiFIX8WHAdoVmPN47CArgBGmPTwpTKCDRlC1QX/i01rBaPsVCTCSLDJInFN26iCoWhjrBDh7LpCjRworNGlNWWYWtJs9THHVb8+1UyQmVOcaKYPB4rkrMz2cNDGH/JiKFzGWjFkMPFq6UAAHU5kHkUQhcD11t3KdvBfnxSI6WMIG0RajkdMh8cB8GlL0HS+KT/a9WNVQXNBb8M+0XPNObzfOeNpY69XH5Z68+wBpLk5JIJE5sPHWZPmOB85ozbtQT7nTttf/3F6e9jHinJttu24SiU6s9XVC5YqKyv13NMCoqa/z714M/t6b1BYfPnz4+MjNRHoKIbQSwUqzObq85VRH9z54c1NxKSKrpyNR+D/eNArYA5Z2HG3+AUXhKTQryjFuMFDivCTT0dwS0C5sTCu4dBUAm7h0BJ4rPuz69S163y3FWAe2Ajss3ZGqtnV6MCoqa4T9mZtfN61fX/XG/KunXriLa25kuXPO9AIIjlqpVKHq/5wIFWqRItEovNCH3cTA2upyYd+JXgbeQsUigXlghKAiweL5DNmeMrJsD982X6JJ18gsJGA6qGh4YRgSghgWBmRvX3179IQ7NUktkiMCGMiXjsBgqyDNIqQv3Up3XwqkmbucRZah3akWJYo0e3pqcrDZ+OvDXRTYqDrIRKAIAr3wDVBAbrlfiG+5YPESGXHNSmZYjgYOJvIG+FhI3tx/7I+ONm1c0vB3w52G6w/t0LI4d9EvDJmYdn9ufsf/Zs84GDFR99jGeznaNPW3//XYefZIRCYY0b6/xPtNW6b1rv3SudOk2alaXLmlOrj57M3/Rp4pUtmc1Xa+A+vzcozNTU1MnJqZdNxE92jpizNdx3tofCg0EUKvKPFP3w5c3Kmk5UgG9sgZjF4BAKC66Dy1AAKL+cSritECEtbmtGUi2eEMzwngQLbwLbHo69C/nnnqxje9b2a5XXVoauEg/g8k1bTJpM83fp+6V9NeArTzPPr5K+qpPUAUDL8eOyomKrr77CdUcIHzUwkDlqZNNf+86duk9WQ+AIB2NqyYsFPAn6jDWuD1m8LDWqfjIhI4VMcBlpy5KhR08YIo1v6QMmDlBoCIWp1ZJbt+iDB4Mh84DD++6T1BA54ykBj0BuYK2ktlbSqTJXQ1YxJQ8nIDR5LdAy3WaOHg0qlejyZUNHj0ImEPswmS3K8luXoegqDP4cKGy9KCzIQ0BtotfQJDXaMsKZu8OAhZB+GGrvA8CNqht7s/dO85g2zWOaoT1c5L9onPO439J/S6lJecoZsu2X+h9/ZEVGOp8+pSPQGnA4zowZTqdO4mi08jkfSDpRzU9Jq/3ps+v8hFpEjRJDzAYs9vn4j2G9QWEvCkw6aeggu08/67/45yGIvwm9WX7yu7txV5+Zudz4GRK+Bf8Z8H60Zt9fY3YJekUoQUTOK4aTmM9EzbDtYO5FsPaH6A+h9NHvW1p92p77e6a4T9G4D3xWTGyhNDFKqWXnb+vTVRKetGXoFqVaufrmaoVA0PjndnrYQMawod01FNxPP1VLpc03iyVEGD3cqJ+T/LPgMkzPL0frLJKII/qaPRUXMm1SHwEd15DME4jkhhhiY6DkGiikel7elpenamlhhA/Sv4WiMj6+WCK2JndYbdBoBHVmiCEKdfOxfCUonKI6Lvg88n66uZLdXEWXLhsxgJPf7oMAFJ67AyYO0O9j/QvavBuIxxFK/07WfnrIKqCbQ/wX/LaW9cnrPUw9vhzwpXGPeEPYBneO+6obq3itPM2Rpt17mnbv5rz3nu22rYgeP8ZkV1enkyfIri6Vnyx5lsV2781K25tPUKIWo22/+GX4go/8Q/wtSQTcy0thOlLhFhQUHD16NDY2Vs8UAAwaMWpx0IjlAXIiUnSq9ODRJ3Lwpe2HhO/Afwa8tV2zc0XaJGw+8kCFqGyj+ncauUo1gVmnwcwNjs+EumyBTLD6xmoHlsOa/mv+M9IRzy/GSICP3hC3PKzSp5MOTIcv+n9xr/5eyo+fqwQC7qpV3TiYJBcX2Ygparwp042GGPHcau9Dc6nRs0gNhXmZelEIlA5TmUHT3CgqOLT/vgF19RkL8lYou6mvI+zWLcDh6IZkkj99MBcA3p3XMV7Bg+PBIDE6c4dZ3BPTEBZxpBndxrzTn9XRo1tTU1XNzYYOoJM9C+W0Vgj8ZAM+NUj60dTLSUBrpvPo0gZtcy4KC0Z8DeXJP1xZwpfxvx/0vT4Li9pNRQJl69CtMpVsXfI6FFBBTAzvl19YEyZYrftGf/sXb2LisH8/2dWlaslSafajNDRypfrnH1IUaU1ic9K8jQPfndIH6Sbu6UEK05EKNzU19aOPPiopKblw4UJISEhra6uedfb1MovaGCZiE8Q3eXv23gcAKEmEuJXgMRom/6lZ3EHVaMkvVylAo060YjnoXO+nsGHWaaCw4fjMb2993SJr+WnwT1TC49VrIpVi+WEQiqD1f6UrZQq9fmndJkeSA0zOJVMmjNFldRuFZIthCIqOaDMq/3lezCPFG6OgUCuyG7MDuFpkAocOtBNZEFX5wopqsb7VOQ4CMgMK4/W8XJJ0i+Lpqb9AWFpWPaNWpnZluDh0VNzEIbgAiwCti5Il0Un2qJPQUuAwSlf4Pmv0aFSlEl25asQwjrY414YyT5cGG1rQcnJfAkIsOdZJTrPA95Os+8Q13V/sv1BfNfDOngzL8bOQz5Kqk85f/K123Xr6wFCbHzeBget4eDbbYe9evIVF1eIoRW2tUqneujGZViFVeDBWbww3N+1O4TPtPfvhhx98tUH/enWnwg0JCbl58+bXX3+9a9cuDodz584d/Wtms8krvg8XmBFlaY0n9l2HUx+AuQdM3dcuAJ+/5xJbad7mpbAJ10MTimUD7x5KVPIvVV5b7LfIy6yjvg3H3Q4XzmQAp2CHXr4bBJDFd0zUCLp/UDdvmpXL1dJ6hCqrQqIPqcVig8vnx4LTYKCZGtd6flO+TCXrTKl1ygc+eIDj+/Q2xAhkcB4KRXplJlVLJNLMTLohs8jLxwrkOHh/vvbQswBuQBG/SCh/ao2PX1yF3G1rVjV4LenCV0j28CA5OwsvGb4uWZ3uITyOEIUVqQZbcBb+rgJyM62GonVpUoGqfjShu8gV81TdkBPovT7vRbD7m2zYjZia2G7dihCM2ZOPNzW137kDlcuqoj7Z9t0tFk+OBJj832f9ke62mrTX9+mnnyZpg/716k6F277RUiqVlpaWuroatkuGRMB99s1AkQmu6a6ypNUPZpxojxLgpRcySih8cpPHXH1D4SVcz+9t7Dzkirli7a4Zp4mhfEYjs5ZZd6drHWF5aani8jXe2OAY/k2D9rJ0iXNxxRQVWA51VgmFT+og6/eTkguND59zFgmd5+72dDNVODOo1W3p2Tx9a3QbCS3l0FTU5YWtd++iCgUjXN/UIZeulbMFSnoAx6KTX/sgbpAaVWfyMtuPqOSKmr33AKDcR44jdv3FsiIjW+/cUfENXEpL2gZUDjXAjiVFb6cbnObaYqInESGVHtXyGR7OO1ze1rgaMSdc3wwq+XO+aQggyxIoJhL4e449/jmCEMiurjY/b8mXOtPrFCpvVtSiHklUpv1pUalUjeh1RkbGyZMnw8PDIyMj4+PjJ06cqGe9ulPhPvp1Vavnzp0bFRWlySqigf6pcKfYnsiWDLvAW0aMzrI2TwcAUKo8bhOICKmkj6JY75Si56Xn6xWStTIL5Mq3VytJYor1s60j7irPtFbJ6Yd3KvOBqCsE2fr4CQYeL/EbyULLvrz85ULmQh2RsTrw7L1nJ9JoeJzMhC91dqrasTOJyUT1Nu896854ILhLZXhZ1VnjWr8ouWiON0+63OnPmJUjUltKj92fVTlUL7cATaYeBZB95pcSi0jdrXPPxrJJpCuVlWht1wl+UDVkJ9BJOMTVrOLsWe1ROApUgQf8yVsnWygtmiNmdwVOqHOeebGaw9EnGS2ZTHZUKm9u2Srop+9Gd5qsYeSDCw8tJ1KZtQKEfuFkVn3lHX1G/knYqaUmFWZnT0UD6fGXK1QL/xT+6U30RnF+UPXT/f2fl1qMfJ5XjpmZZX05IWukVzQuk3TqBx+Sj9H88qCYyHAab1V319RadvZsdY94rNDOkZqa2q9fv3Xr1n377bcoioaHh8vlclQ/JCUljRw5UvP3nj17li5d2uECDX+tWrWqsxp09w0tvISuY9f9vXZL1NVNSxOEYhmKotk/xVR+cb3s0l1UbxQ0F/gd9Nt0ZxMqqkM3O6F7IlCVEkXRmJiYDleWxd2t/OJGzv/O66hNVlKS5+1T//PPKIqeKz7ne8D3fPF51HA823phScvvC6/+9lsqiqLCS5fy+ngKr141oMY/B6AHxhvduhpVDzk+ZG3SWt0Fd+3J/GPh1Yv/lujbq9+D0SPTumy9aPSYioWL9Kzy71P5fyy8euR4nu7LZl6YOSdujubv0nMplV9cz/4xRmvrnaFoVGTFxwsMeASX1qIbzVBBNYqim9YlbVt0VfPS6r73DqhOyq784kbuzrgnD25I3hB0OKhKVIWiKLp/HLqlD6poM/qVUzQ0FAwILX33PaVCNi122shTIyUKCWoU0nN42xZd+f7/Eh6+O+tBSD95ba2RH7tO6PoZP3PmzLZt28aOfeQasLGx4fH0nSa0p8IFgPPnz2tSNt2+fVuTew5F0YULF3I4nM2bNxvDu20CiF0GVj6W0772m+bKkKM7fkmru5vPbuLw2c2Okf30r+mn1J/YJHZUQBQwLGHsz1CVCql7tbs5x/TjUxoZVTQdq5ONO3biSCSz+fMBYLzLeB8zn/9l/E+hVjz/L01cbBEKMG6SGwAwIiKINjYthw7rW7ixEHj5RqvrAECFsKK5rbnLlEVzZvtIiJB5vkylp46oawSU3QSlLqehorpaXlampyOsVfrUdiIdCOIG5TbmylVyQVmt+oZYBALPZYZtumJGjpKkpKiEIv1WQ6SQcRi8JgLLBgD6RdiTUDgbW2Tog7AZ5CvAN5NKQCV/9FJViarOFJ2Z5jHtUYrMoatBVPvsliMDXA4//KBua7P5cROeQPo69GteK2/3/d1G1CNtU8bvyVEhyMzPgh23/ggqVd26dT1hhOmiMDMzs/Lyco0/S6lUFhQUWFrqu6G3PRXuiBEjGAxGeypczS7L69ev79279/Lly5pVgvj4eMN6ffkrkDTApD+BQB413FHhwWBUSRtPFcrRNtcFBuyhS6hIuFN7JyowikViAQD0nQbuoyDhWxBr3yFhPzcURdGaw2mdecEEFy5wZs7ULJwhgCwNWlojrokpinnOh6RUqmUPhSITgpuTCQAgeDxn1kzJ7duywkK9yufGAIIDz/FGdyC9Ph3+SyirAxQywWmELasNPXpcv+RjbhGgkEJ5ii5P5a1bAMAYpBeF7fvrPl2BhkxxwXc1xQ7kBspUsry6nOpdd3GA4872JVDJBlJYJKpQiBMS9Lo66zhI+TBgkea/iMEOIhJSkcYz4lnQw60pCL0k+pHE0/as7XgE/6Hvf3vFnQeDQyjc+hXUSiMql9y+LbwYZ75ggWYvvZ+F33iX8UfyjtSIawytasfv6cw21GOio5M9i2hnx/38M/H1G4KYmF6lsPnz52/ZsuX777+/fPnywIED33//fYIhaxMzZ868d+/eP//8c+TIETweDwB8Pt/KygoAhg0bplarc/7DmDGG/ACWXIf0wxC2FGweLfB/siTEjlRrgjcV+xC71C98/LuoVmy7t82F7fKOxzuPj47eBEoZXP1WaxG2k7XUWcaWm5XFaZFYaty9B0cimc6f135kkM2gIG7Qnvt7ntMQu/hvKU0JnoMeJ8IxeecdhEJp+fu4fquJseAQCkzjE6xmNGRwyJwOai1a8c7kPgIarv5WvV6Rrk6DgUDRvS4pTrpFtLHRR52irFKozOUL9ZPGD7QMRADh78tkoaaq/kQzH4OTTlP79iVaW+sbpn93N1j7g/1jTUGOrwlbos7IbTC0XafR/SSoQJUlAjVaJii7WHJxltcsLu0J5YzwT6GlHLJPG+xTUirrv/ue5OBg9uHj/LX/F/R/CIL8lvGbQVVdulZOKhZLHWmTxj5aqePMmEHrF1L/wyZlQ0PvUZiJiUlKSsrs2bNnzpy5c+fOFStWGFp796fCVUjh3DIwd4Nhj/NxSqvq/Cg2NW01cQ8NiOg7ln+sXFi+qv+qp1JSm7vDwMpNzHoAACAASURBVCjIPMqRaDfy3eePlKBCWSKvgxCdsqFBeOGCydSpHSQ9F/ovrJXUPqchlnOzWoqDCaMfr9viWSzW6EjB+fPqLvXOm0ugLhu8Jj1PBzJ5mQHcAH3WJRAcDJrqRlXB4QPZXddLpIJjmA4KQ1Wq1tspes4ij/+VjaDotPl6hf5wyJx5dWPcJG58TrPrtMHGDAqCMCNHiW/dUku60iCovAO8POj/VDj+xMkeagQSzhcb/skiiB+dASbll9N2ZO2gEqjzfOc9/Y6OBqu+kLQNUMPSgrQcPiIrKrL8cs2TqipcGneO95y4krjsxmw96xG3KjKji8UkWLTsCbMdh7P+7jtULq///ofeo7DCwsKcnJypU6c6ODhs2bLlypUr8MJxcyu0lMPE3+CJGPGKg3cAINuWy2yQn7mgl39BLBfvvr97sN3gQTbPfB5DVgHD0rdGu4FDIBOJoWw6wn54/Kltt82Hj6AqlekHczpcH2YTFsAN2HN/j9Iowx4AKmvEzGYl3o1BIj31sEymvaMWiURdqibknQVAniecormtuUxQpjWoVSuGDrITWRCVeQK9Il3dRkLDAxBoX6tqu39fJRTR9ZhFXk+pYtTJwIvl7qLXT2Z9WsGU5jG1qirPT8caPTLMyEhUJhMnXu/KiD0CZAb4PJWy08aSLjYnIuWt0jaDXwyXaeFtaonwRuWlskvv9nnXhGzSgVth8OfQUAD55w0gAqm0cccOxtAhjGHDOs7GfOebUc22pG7Rs6o9OzPoSgie4sqkP5VageToaLZ4kTA+XpKc0ksUduTIEU1agaVLl7733ntr1qzROONfGFrKIPl38HsXHB+/07UpuRyZmcReOnv5UBEZeXixQiDoesfS0fyjIrloaeBSbaYjHYatMRUXQoF2TSuXKeFCpBl3Xyb/L45M3dbGP3mSNWqUVl38hX4LayW1caVxxt107JkCBGD0xI55d2j9Qsiurvzo011TmF2IxotsHNLq01BAQywNkMmfPMcHD3Bcny1HmpySnRhi4lu3AI+nh4bqrkOlVt86VSQlwNz5eulYSBv4wlPFSlT+tfPustYKo0eGFhhIsLDoQj5MLoGcaPCZ+qxMbsBQO7IaYi8Y7NQnkIkKZ5Sjtghpcn3fW5uErNckMHeHpK3612l2LVElkXBXahHmpBPpUQFR6bz0G1U3unabZvMIhWKJDXnUcC3TebP580nOznUbNqByeW9QWGtrq4mJyd27d8PCwt566y1PT896vZWAegTxawBPgFFPKce2nCuSqaXuc4eRSLhBMzxoKti3N6tLE+xQ3qHhDsO9TDvJNRv4vohiA/9+05lPlDXGkYKjFx26pvmXf+qUis83nTdX68XhtuF9TPvsz93fmRKT7o9TnC8QMvHe7lqi6tlTprSmpsmKO5+M8CugJvN5TDAAuFd/j0qgept561/E291U7kynVrWlZnb1wlh4AtsOirVTmCTpFrVvXzy7i33pBw/lslvVdkNt9MmopJIrSn69TgaabDy9mt70XOHHOBxz1EhxYqJGx60TL0A0yCUQNPvZM5EjHCVEpPiOMd8UZ5KPQi2b2zjRgqotMRIOD2FLoSYTyvQKR1fW15skJ7MnTSK7ac9QN8Vtij3T/o/MP3S/w6ga4g/kyfHoB1HabXaESLT6aq28vLz5yJHeoLDQ0NCvv/565cqVU6ZMAYDy8nIzwzPHdBuKr0HBRRiy6km3dMnZFJbaVOlL0GhRDA61ldiQ8Q/FeQ91beA4nH9YJBct8Os8FyEO/8B6GjQWwv2TWs/bDQ0QkJpo5WRpkwBUqpbDR2jBwTqkrGZ7zX7Y8rCDhok+uHytnKEElwHaF4LZU95CCATBP2c6N8FiAFDwmvg8A59WlxbADSAYmPnx/fl95Xj491B+11JibhFQfA1UHQ18lVDUlp3dpSOsrFIouNMgYOFnvKPHjlQ1mrf5PFvFkQeCz9DBFlSL59xBwYyMVLe1SZJu6ZpFmnuAnZZAHzwOR/NksYQq3a+rVpysj0nDpzmqnJ/NWfMIftOBYQnJv+tTW8P/fgcAi2VLO7X7cIRF/ovym/KvlOvyJh09kcuWqC3DrazMO93nRB80iDliROOf25XdZA/porBp06YtW7Zs9erVb7/9tlKpnDdvniZk/wVAJYe4lWDuDqGLH7+NCqUiuVmCCjxmPp69z/zIX4VAzP6czmoSyUWH8w6PcBjhY6Yr5riGHQy2wZD4Ayi1T0u5U3zwCLHkyE3hlSvyigrTuXN11DbOZRyXxj2Qe8DQ+85MrJLhYPIE7b+NBDMzxvDhgpgYVNmJPyUvFmwCgeNk9MAL5cIifpFBWcs0sDSn2UXYsFvVf+3L6todJhNBdcdQldbbKahKxehKneLvnZk4FH3rY72U+HP/d54jMRfYidymDwOAQMtATbyI8XPJfv3wHE6n65KNhVB5B4LmdFZ8wiR3ALh8zrC5pEQhOVFw4kE/PgBURqd2NtuEfh9C4SVo7CLyRl5Swo+J4Q8MJVrrWrOe4DLB1cT1z8w/OxOnFknktUn1Qipu1ntdRPNbfrkGVal427b1OIUVFhZaWFhMmjTp4sWLs2fPtrd/cTlO0vZB40MYvQnwj2cKhUeu0REWebA5Qni8pGhvwyD7m7CblXFXtKshH847LJaLF/sv7qJFBIGIb4BfCZnHtE+AAt0F5GZ6DbX+4DGSg4NuTXcijjjLa9btmtsFzQX633RNvYTRqEAdaVRKpxaQydS3lU1N2uXlhDVQlfacs8j0+nQ1qjbIEdaO9972FHAIbRktBUU6rQznYYAjQHHHACtx0i0ck0nRKdN65nwRq0lJ9OdonWh3QMGBK+w6kxZWs8+SRyFygdzAanF1uzaWEUDweGbECFFionbnTsYRwJPAf3pnxZ3sWUIOQVkslhsiextbHCuUCycPmSmgtTAaGNo1qQGg30dApEDKH7pra9yxE0cmNz/jxe/IFAguKiCqmF98sfSi1gsO7LtPVUHoNLcuN3IT7exM58wRnDvflpfXsxT2srjzZSK48TO4DAX3xzu3ZXwxKV8twDU5Tejo650/z09MgMzY0menMEK58HD+4QjHCL0ESVyGgf0ASNr27BznkSE22YeAEPk4J9N5c7tUI3nH4x0akXY477D+933ubCEOhRETdO2Bp4eH4zkcwTlta095ZwHQ53eEkfFkX3Nf44pPX+ynBojena1rfZ/CAruQZylMkpxMDw1FOk+KU8uTFMWVCynIRx92vVr68O9E2gMSn9jou2pC+0FNsK4+2UC6mEuKRFpW2VA1ZJ8C91FAt9BR3CvchqqCC5dK9GwOBfTvB3/7mfv5mftZjPMk4EhlJzqZxtLMwH8GZJ0ASaehWPKKCuHFi5yZM1T0riUuRjmO8jL12p65/dnl9YpqsSpPKDInDR9kp89dmC9cQOBw6n/c3LMU9rK485O2gaQJRm588ljRgUQSjmo2WQsTUcgE11H2LDkcP/XgWRNMIpd0bYK1Y+gq4FfA/RPaKSzIvamtxMJmMGlI11sCmCTmJNdJ8WXxLbIWvd5UNbTktAjouE4T2WusAAKBNXaM6OpVtUikhcKs/cDU5TkpzNfcl4w3Mv2Si4MJJ4zLFqp279XpcnIZDjUZIH08MvKyMkVVlW5H2P5f08kqZNQ87w7hJlr8qKdvkjMRIaGlz+qxT9rsfUz70In053SH0cPC8CyWlrlk2S34f/auM6CKc2nP7umcfui9V8ECiCgEURR7RaNGjSUmJrZUzU2iiSmmehPjNZYk1xhjr2DvioICUkSkN+n1cHov+/2AIODplNzk4/mlZ3dftrw7OzPvzPMI6yHkRcOHT5vkJcNBYYqp5e8P6h9UCipfCnwJAOzDAwS4NkIF0tlv1BOj14FGqa9nDgC4e/cBHs9ZtswklxOQdSPW1Yhqksp79qIf2/8YxWDWMlPXfFAq1WbNGmlGhvjOnX40Yf8T6XxRA6TtgaEvdtbiA4CwppnaQOVTuA6jdC8pJszwFVih9amNXYtuRErR4cLDE90nmqIh+CxN4xwK97brXJpU1dRoH53Go8SnpzNNckkCFio1yjMlJvHk3EqtpivBLdy4XiFzxgxMLhfd7M7AJ26CmnSL9SI7vmFqaQG3INwhvDeDLF8aLOAQFNm8jGz93z/v8aDVQEVylygyBQAMJMIOHslntqlIoWyjKpalR+8QMzERnuf3wSQCtRvlLA7BhdiEGFUDMfJi4/G0cbGiGzewHjFK3gkg0cFvkuHDiXgU70Nn8NRPa4Sm/LnDhYc5ZM5E946IhBblQEapFWf1cFJbe4P/FEjfB0od9CGqhgbB+fPsBQvwtrYmXmyMS8xQ26E/P/65a8NJdl4zpU6u9qaZEs4/y4EsXEDy8W769jtMo+kvE/Y/kc6/tQ20GhjXjQu8+uB9FBDnl/S2cyMojJjmTlXDH134qX/P/12sFK8ettq8E4jZCG2V8OS0ji/Yb79R2qr5hBarWpKCb7yM04vpFeEYcbz4uNaEsumHN2uUCMya6WN0T8qIEUR3d8G58z1dMEzby7XInOYcDaaxLBHWLbR/J1SBg1sHCvSW7DmHAYXdNZaUpKQS3d0JetKvmblNvHtNAgZu1SojhWCFey5TclEhrs33OfvVmQ4r5hWLleLeXCA9Pl4jFEofdsmsqxVQkASBM4Fg/JWJn+4NABfOlRrds0ZUk1KXssB/AfHPpLD75AgJJlDn6Dd/o9eBjKczkuDu+xkA9BUD6cOa4WvqxfVJZc8csasnitUILF5hXrYBweFs33pbWVHBP3Wqv0wYADg4ODx48ODixYsYhtnY2Ay0/WouhNwjEPFa1zW1ltxypogjtBGwfQ1F3VPiPAUMnDCrrf21ESlFR4qOxHvEm+GCtcNvCjiEwN3tPdo1NHy+4Gwic+YMm+lBBJRUftykApxFAYsaJA2dQm163R8ZQmmUq10pPeqb9YExfZokPV3dlUek8BzYBYJtrziIsxqzcAhuqO3QXj5GBxvq8Be9aUps97cZukksUBx4xkB5hyOJqVTSjAx9UWRLm+zWr/lKFJa+FWagnVurUj/5JoleReNbtfl/PK2nGnZnOsw+VItpc1tze3OBtOholErtpglSehXkAgiZb8rhQb4cIR0nKRQY/bQdLTqKImg3dSIUQYIpdGDV3NST0XMfA86hkPZTjwmsbmrinznDmjvX8ELk84hyigq1D/358c9KjbL9c0JvUaEBdHsbK7NN/4Q4q5Hhrf/Z1V8mLDMzc/Xq1SQSKScnh0AgfPvttwOdzr+xFYi0HnJ7Tafy1Fql93LjykBj5/mQtXDw9zwAOJB/QKKSvD7sdfPjBAReeA9aS6DoYtefeceOaeVyzvLlDqMCBTgusRxMIdePdY11pDoeLTpqJPYpI+ABeWGyqb3HzBkzQKMRXvxzqUjSAlX3e5nIb0+EDbEZYoW36v2TnDjWnTiSw+Cqdv1HD/G/dxwIaqGlGABkjx5ppVKdfUVqtfaXbx6SVTBqaYCrk15BX3F9S+HWiyweh2/DC948C0/S2zw7zHYYDsH1NpYkkWhjY0Q3bkBnTJR3EugO4BljatJwlD1NDRW1horvZGpZYllivEd8t6ZuAM95UUqtVJBcpd8RWw+tpT2UCri//QZarfVrr1pwvW8Me6NB0nC27CwAXD9ZqkLgpaUWLvjYbdqk5Av7y4T1hi+sb1ByBV54pyvde83NbJbKWu6l6qYLqe9zEeEsssarCwUltQ1HCo9MdJ/ow/Kx5DQCZwDbA+7v7PwBUyp5hw7TYl5or2amjXUmo9SKk8YdMRyCm+8/P60+7angqaHdmghCChIZ6mDiCRI9PChDhwrO/xlLFl0EraaXrd0KjeIJ90nvo8hOvPbKcIkzGVcoOnS0QHc6DKA9lpSkpCJ4PHXUqB67YFr495cPmAI1M9ouJlJvy1TNjazGHVk0NV02TBP83kydWmrPwnA8Jcg66GHTw15eHT0+Xs3lSrOzAQDkAii5BkPmduo5GMXMaT4KFPhVhpzuq0+vipSihf49SzQIFLLMVc1SWrfm6VnWDJoJLDd48NOzGEIo4p88xZg2leDsbMHFRjpGhtqH/vL4lwePaumtSjSAbmupogclJEQ5bml/mbDe8IX1DZjOENEldaXFxNdrZVqx7zJTScGmLg7EY8j+Q1elaqklLlhnmDN6HdRkQHWHoKTg/Hl1aytneQdDgHv8SBHwtXlSTGu8hSjBN4GAIxwvPq5vh5T0OpYKdRhhXtjOmD5NXlCgrKrqSITZ+IL9kN7c+8ctj5UapQVFrQawYVOEgInjJTeeOPtcfRzLFWx8202YODWVMmI4+twy/w/bM2j1Ck0QY9kS3d98OV+c92UickOiQTSUBa6+i2JNOauRDiPzWvLk6l7JtdDGjkXI5I5YsiAJ1HIY+qLph1tR8Jgb1UaMa2zVy3txquSUN8tbJ/Gk54ujNZi68byeplQUD6Neh6cpUN/hbPJPntBKJJ0T2AKsHb62Sdp05WR6b1wwAJDxhEyrYf1lwnrJF9YHGLe5aza09ORdOrCRUCvT2emGBdkK7VCHavsJttMsdME6sr5LwMq6wxHDsLbfDpD8/amRz9wEQhiThjAqLxhvIeKQOXFucecrzis0unPb969XqRGYPcvXPBM2ZQqgqPDyZZByofJuL10wAMhsykQR1ChTq1kgk/AbPh4toqFNV+sOHsnX4YhVpRBEAnlBQY8oUq3WfvvFfVKFWOZJ2bAhXGfmq+TA9fovU1lCjsCG77U1zi7U1Bs4ynGUSqvqZXUYSqHQXogWXr8OGAb5Z8HaG5zMu3Vjp3jgAM4l6U7ql/HLclty5/rO1W0B7TlCppAhYOlW/AaA0GVAZrSXuWJqNe/QIeqY0eRAy0UCIxwiIgjjbdtYqD/VtheiauW/JRMQUn+ZsN7zhRmQwjW6FQC6ljWrZQrIkoqB5/1ijFnnoI18jNcSbfPjevXyESgwchUUX4bWUvG9e4qyMutXVgLyLELxnhMlw0TyBybRuSX4JggUAp3tZgKRglAna6Wr2UyyWSeIt7W1GjFCdOUqFF4ArRqGzOmluXnY+NCf408n0vv2q0SnEtduGS1i4kR3m3Z8n9Gt/Nh7PCildvkpoNV2NWFNrdJvt6RQa+XqAPq7G0f3GFAlkRcfuFG2+ZpVEVmOlxPm2A55byaBYsbdC7ULJaCE9Ib03l5afLy6qUmWkQJP71nwCQkfZs8larm5ujsZzpSeIeKIM7z1LjE7zxmOIOjTU3qU50k0CFsB+YkgqBVevqxqaDTcEmcK3PMnqVAFNarE4hEETxvoTTS+VVt/mTAAIJPJCQkJCxYsCDBf1dWAFK7RrX+mSZ+dXsnB2xSURp3gjKBmCALxFfyL/CN1NjWkcjDgopv23XkV8CR4sKvtt9/wDg6Mqd00ZRE8TuNPYGCc2uRHxkdyjHBnuJ8u1VGokZhUSsQQG09Llk3okyfJi4qUKcfB2gccQnpzrXK1PLclN9Ixsl/SA0zSe9tekLhQCCXib9+/+4y51OMFwBFpRQU4FosypCMKTrpUfnBrGo2npoy2ffOtkZ0zQqtS1956lPdlYs1n96hFJA2qVo3BB305S1+poKFJjieH2IY8bOx1Oiw2FiEShSf2g0Zl2VoK4qhiKOFuWs8yV6VGeaH8QpxbHJukVxKNE+ghJLaRa4kqiZ6IOGI1IAik72377TeSjw/thRd6c7GZuU0cLrnSteBw5X6LY/CaQ2kIIK6Lw/vFhO3atSsyMvLTTz/98ccf3dzcbGxsDpnDj2FYCtfwVh0BM1dArsDx8VzXCealZn578ptMLYudPZyAIccO5vdqhlJtYdhCxZ0TkrR0zstLn9cH9V7wggKT8W9WGp+pgMzxmZPZmFkl7LmKVJ/VKiIhHk6W8CMyJk8GFBWm5vXeBctuzlZqlKMcR/VThoCIRzdtHm0T70SWaVP+k/f11tTk1FoNngKuEUh1GzVqjFiuOXOudNvG5NpzVRoUGbkq8OUX/XiltU/PpxX8dOnJx2erNifDNRGdz5BQJTCZHvjVTM+Zoy0+n1EOo/K5+SKlqFfvEp1OHTNalJoNTDdwHG7BCH7eKhUCadef9vj9RtUNnoKX4Jtg5NsQ50FCKZUn7+nNLAfNllw4KC8o5KxY3jWGsADXTpSoUJiXMLZV1nqi5IQFIzTnljLF1iIbIcvbpe9NWF1d3d69e3fu3FlaWnr58mUul1tUVPT111+bPq5hKVzDW59Hxf5kAkKym29efpqn4B0vPj7Fc8qk8OEiWyJWKuqtIzZ6PTefjJLw7Pk66n0IVLLcRWVoYagLZvvMxqG4M6XdKvXTshsZMq3NcAtbIPC2tla+zsJqcm/EitqR3pBOQAlG9T4sgFqhaskrr0h6ULT/WnBt2XhnwTBq81heDfF4SubbJ7Jyl0miP27UjCr5KMn3Tv5MDW8sTTiHIrE/kd/4aYbkv5X4VBWjhkZQECVskXIU4rRlVMjHs11ih/fyrCIcI7SYNrMps5fj0MfHqPhKOT3GMgNBJmFKJzKxXt6jBvh06WlnmvNIByPSXC5jh4uAhxUq9S4rjVnf9gSHZ1CY06f3KsnwqInBVeECGDG+I0c7jf718a9StdRsE3YqX61Vei2P6eU9152eLy0tnTBhQkRExJIlS27evEkkEt3d3c3K5RuWwjWw9XkpXJQnG9bqUIdWNz3lwVMzWtuvyq9KVVKfFp+kpCS2B07dYrX7x5SwcFOd3udFSfF8gXc1meEvuXTjigbVsf6NWCuH1VLKjtxPjTRONO6H8ztRcMKj1gP980OS9YDCQfAcylMwJomqDyNsuIpi/OULaUrb8t5MiyuFV5xR52sXr/WR3dJQnoqoPISpoDNQDoqgRAAiUDAgKzUyCkoEIqi1Gg1oEAwwVKsCtQbBtIgSRTENqpEBKPFtKoJWTQYlE6+2pmBkPAACwM271Tdk6BrQEBHisdRjUylTLbvz7XDlZlkhWP4Dbq2thYOQbHj4OsrO3beGBnVQX/C0vIfChxPIE873aMDQBQqHH9Tmde37g3JfHVVHxOZmjwYya5j4/OWLGIIzOuH1IeuuFRtFHezqk5KShmqGPlA8+Pj0xy+QzYhMiVXCEJVHGbXi8X1Jv5gwtVrdbrDweDwOh7NgXDs7u3a9NQDgcrl2dnYmbt26devWrVsBAEGQWbNmAcCTzxIx0AatiY90tTP9BNrkbdtOb5vmNW3VC6vaf/n6aQqnFUZGRjvZU02xX+1/vVv8+/XXPARn68ub7iKAkat0HvikOtGJ5zIkfAjd2UjfmXWd9Rs33qAOp8a7xwOAQKSouJKqcCQtWjhe5183DilXnV5ditqNVKlsZlle13os8VijtvH1Ya/PGtbb4tjaO494d5/SxUw8aqvB1GKiUGgjJDkx6J52NBcbqoNNj8xm9apVqrx73h+Mh9l7YGBx5fqVFlkLYDCrF7cOjp+sdkIItfwwiwZJSkp6841pX711G9dCnjWrox7zh6wf0Hz0g5kf6CZo7Q6tSl3+0XVXHjNI1wk0bNkiJBLsPBtn+mAQPMvohNedZMhrrrn8BAtiLFrYsUSWfyM/g5vx2bTPTK2C1mKFH56XacVRmxaYtfBiXi5sx44dZDJ56tSp27dvJ5PJZDL58ePHpo9rWApX51adaEwvZEk5Ykcpwxz7BQD7n+xXapRda8HiX/TDY3DyiIUURRqhiH/qNGPqVIJ/GDz4SZ88jNPsUARBq0+mGR1wjNMYR6pjZ9f32cRSIgaj4z0tf5iF5/FEpVWwv9CoJohBVKortZi2V4kwLVZxJqXwX+fhiogqoYqZIm2clfPWMUO2zQp+e4bvgrEOEYE0J9se9gtTKKSZWeDBhPJbgGEDbMJGOYwq45VJsF44BSoZlN2gjwlRVlUpSkstHsY6hM2UYpm5TQCg1qrPlZ+LcYkxxX4BAErAq31Qhsa6Mb2nlKeayxWcO8+ck4Bz9jGRzVW3rT9arEJh0ZJnWZ01w9fw5LyjhUdNHKH8TAod2DCM3Hv7pdeExcXFyeVysVgslUrb/yEWi5XmMPYblsLVuVUHMKwtqUSByX2Wx5p1Va2y1uPFx6d7T3dnPNMgCA2xE9kStCWiplapBXeKf+KEVizmrFwJo9dBW0WPfqNOcALcBKQ2SgNZKTLyMqAIOsd3zoOGB+06ow3ZLUISYqDo3IQE5Fmw8aXPnKsoLlZWVFg8TIW6wgpvNdTGwtbI2juPCj86T8zAUAyRBCldP34h+KPZbhPDjFbzSbOyMLlc5BcIokZoLhhgEzbScSQGWKW60vIhSq+DSkafuwxwuG79kmZi5ixfDUDy5acAcKfmTqustZvUqTF4vRilwhStV3oyTfGOHsVUKs7yZTB6LdTnQFWqBeeWkd3EaFPhAxk2XWrBQmxCYlxi2nv4TEmGah6KxMD3WRTbJw9OtwlDEASvC2YNbUAKV+fW51H041GGlqMOREksmll/+pe8X9Ra9etDe5bjT14YQMDg2CGzlyYxlart4MGOasCA6cDxNECGyYn3ISLkihPGp0h7peKZsjMZ2U0MGWY/oheN9FIuPE2BIXMY8fGAor1xxMrV5WEOYeaS5QOAjCfM+zIRrgjxGoJsmMbvq2n+L8fp5IfQCUlKKkIiNQXFACCdLd8DhmDrYDqRXqG23PRDQSJYWeOHTbIKDRVetfz+O9hQJbZEpFoilalPl562t7KPco4y/XASiya1lTKkHEH5s+IMTC7nHTlKj40lenjAsEVAs4P7ljRX3zxdokJg0eKeC2vrhq8TKASHCw8bHaHs0G0rhE6JtTerOsqSQLJPYFgK16hQLlJDlWAC38XmWesmadPpktOzfWa70Hsu1o4YYiuyJWiLzXbEhBcuqJubrV95BQAAxcGoN6A6DWp1VxI5RQUL0TZciaaHXO7zsLeyAi5OmgAAIABJREFUj3aOPlt69s6lcgsq8ru/P+dAq4ag2Xg7O6swy1+hZmlzi6ZllIPZUWTdvce1X91nCtgCG77HJ+N8F8WaO0fFqalWYaFyKxuwD3qexLW/gSJomH1YucrSZRC1HEqvQuB0QPGMyZMVpaW9iSVHxLqQtHA0Kft+/f05vnNQxLz31GXeSASBmjPP5qfg3DlNWxtnxQoAADwJwldCyRXgmsfZn5bdSOeq8EHdXLB2BFoHxrrGHiw4aJi2SM4XE4pBgOO6T4549qseFdH/CRPWS1AJHIzTiBLMcwd+fvwzBpg+gaJJL/oTMDh22BxHDMO4/91P8vendjLwjVgKFLaBhAI50paC0CoS7xsdO8E3gS8SEuolckcSh92L1EDeSbD1b++LpMfHK4qLlZWWxETtRermJcK0WOG+q9qLPAzRdhTHU82+EHVrq6KkpKMo33s8VD0AlWyg02GOo7habntcb77vehsU4naOScaUyQger5sN3DRMGuchJiC1mbUAMMfH7Co/pocjvyutPoa1/X6QHBxsNfLPItKRqwBH7Nr4bQpunSpRothiPQ2qa4avESqEhwoNVY+WH7hNQsk9yZbvfPWPNWFtsircuQNaqRkeU724/mzp2QTfBEeqbhak0BA7kQ1RWyRqaTP1DRHffa6jiGgF4Sug6ALwnuo8xHP6aAkmVGYa55ge6zp2SFsMUYuLmuRh+Z3i10D1g06OY/rEiYAgRiRa9X1pG9KoCNV0VjW1TPHk80R6pZWAxHP9V7RjZJBlVyBJSQUMe2bC1HLLkjW9QbRTNACk1KVYcnDhOaCw2tl1cBwONWqM8MIFixclEBRYw1k2YtZo3CQnmiXpUbsutPri5GRFebn1ii5N3VRbGLYQHh0BUaPpLhijTU0MZOn70AZwAuLc4w4WHNRXISysaaY20PgUrkNElw6KliLIPfKPNWGcWHt1Wxt3/37TD9n7eC+KoKtCVhn6xC3wI2Bw9A9THbHWffsITk49OoogYjUgOEjbrWcKIhBEpgO79paRfiMcgvOvi2mh1viFECy/U09OAYZ1EuwRHBwow4YJzc8oazFtan2qN8HbxMhFxhMWf3GFJbMWuouDP5lFYTMsvgJJairexobs7w8A4DYaCJSBjyU9mB5slJ1ab77p1Kig+BL4TwVcx0Nkzpihqq+XZlneOu4ZLVDi5A6lUZYd3kGrX4moZQru/t8ITk70yd0psKPeAq3aqL5RJ26fLFGhsHCxoU/UG8PekKgkBwsO6tzaTrbssqR7ge71T4BI+8eaMJc5ExmTJ7ft/03dYlL7dBm/LKksaYH/gh6ccDocMWu8pkhoiiMmSUuTZWdbv/Zqz44iugOEzIfsgyDV3aTqNS9Kgcn4t41Ec5dvVLJk9DzH5MSyRMvv1OMT4D4a2M+WXxmTJsnz85XV1WYNU9RWxJVx/fAmuWCiupaqb5LpaqZsBBb0xhToTXYWwyT371Ojozr8XAIF3EZD+e2Bn3K+BN/0hnS11swGr8pkkPG79kXS4uJQKyvhBctjyUv1iRV2jzmN9PpGC+s86LEuZIRa8vMFaUYG5+WXe2pBcTwhZB5k/hckrSZlwXhqQhDLxiAphR/bb4L7hEOFhwSKntwN7WTLPduJnqZAyRWIee8fa8IAwO7ddzCVqvWn3abs/H3W91QC9dWhxokoJ77oT8DgqAlLk6279+AdHFhzdZGcjFkHKjlkHdB5YEe/kcpIv9GjmzVyHLCG4M6UnjGFU18HGh5Dc0EPpRz65MkWxJL36u6hCOpLML6qIKlvrd+ZTtZSkYl03wVje/mU5YWFai6X2lXsw3s8NBf0MtFriQnD+0pUkkctj8w7rOAckGjgFfvsvaJQaOPHCy9fxiwiOubKuMm1yTYjrFCAxMRiy67FbWKYEGkj1JCByWYl6JrAL7wLKrm+SKJbFuxkiRKFxUuNZwneGPaGVCV9XmywnWzZa2WXqYJp4dpHwHLrxgn4zzNhBBcX9oIFvJMnFWVGlorSG9Lv1d57behrLJJxQtfwYfZCa4Km0IgjJsvJkWZkWK96BSHqotO0CwKf8ZC+B/R06nstjFJjmsZzekuCSyv4NJ6a4M+YFzC3TlxnId9L3gnAEXtQIxAcHSjBwXpVpvUgtS41yDqIihjpXpA2cmt+TCODFW4Sy9zGe72JMAShjh7dzYQBQMVAO2JeeC8CSjAvHabVQPEl8JsC+G5JIuaMGRqBQLdKsTEkliVqtJrlsbOEdJzkCV+t1lp2OeRwBgXPUMYuQum6SJNs/CBoFmT8DDK+gUEeZDYweGpSMMsUAigfls8kj0mHCg/xFc/GrLmdw1JZy701FGvms11zj0H9I5iwFfCkf7IJAwCbtWtQCqXlxx2Gkzj/zvq3E81pUcAiE4edvNDfaI1Yy3924W1s2PPm6d1j9HoQN8MT3bpqFFuWmC1gCFmiOt2B8LmTRQAwc45fnFsci8TSSb9jLATTwpPT4BvflZ67wxGLj5c9zlM1NJg4klApzG3JNVqCJOeLq3bcp4AVGs9yjeubPnBxaio5MKCbFJhdENAdB96EkRDScLvhqXXmpMOqUkHSAkE9CcKo0VF4a2uB+bEkBtjp0tPhDuHuDHf/aCcrNZy7ZGGpB7Usgy+ttcJ89ApNjn0flBLI2GcwC1aqQGHJYlNJFl4f9rpMLfv58c+dL6fkao1MK/Z9+ZmXCmo53N4GzmEwZG4vH9nfwITh2GzrVa+Irt+QZetNjl6suFjILXwz9E3TRVtDQ+yEHIK6SNjG0+1DyXJzJffvW7+yEiHr//h4xYL9EHiwS9/ak4F+I4FIgVZJJfZED1dGO5vdrepbPDnPvLtTeReE9To5jumT4gHDRNdMjSXv19/XYtoYZ0PMARqlqnz7TRowYBy9r+yXViaT5eRQe0hGIgh4j4fy24BpB3i+RTlFFbcVt8haTD2g8BwQKOAzscfPCA7HmDJFdPOWRmgeh09GQ0aNqGae7zwAmDHNW4qHonuW1HloxWL+yVMYWkZBaeXH7+rJ1ASC/xRI24PX6A5HklNrmQK11TA2k2nqm+XF9JrrO/dI4ZEyfhkAlJ68SwM2hHUnW36wGwS1EP95Lzl//h4mDAA4y5cTHB0av9gGuiS8pGrpzpydQ6yHTPGcYtawE+b5EbRw+NAT3S7YDzvw1tashQsNDYEgMGY9NOXr8xcM9BudOFlM1ELUlI6myATfBJVWda78nHm35tERIDN16q0S3dzIQYGmx5IpdSksEivYJtiAr1vwzQWmmi0P0bpPHtlXD1eano4pldSo6J4bvMeDlAsNuQM82aKdozHA7tffN2lvTAuFF8B3ok69SObcuZhcLrxoniN2quQUk8SMc48DACIeZQxlM0WaB5kN5l4I/+RJrVjs+1qCCHjIE6UhR0zG92m+rDu3kFQuQ2HJEvN4rjaM2GBFsNr+cLtKIodsqQh4vvO7fBrFTZDyPQRMB/eo3j+vv4cJQ8lku42b5AUFOlUz9zza0yxt/nDUhwiYZ9EjQu1F1gRVoeB5R4xaUiJJS7NZswY1qv4bPA8YTgbKXHX2GynV2rbsVgEVfSGyQ0KmXdnBvFhSLoDCczB0QY8sTNdYUpqTozZBdwoD7H7d/dFOow2UU+TvvMiW2AhdxL6Lx/XhwxXfS0EpFKuw53w673GAoANfWuHH8bOzsrtbe9ekvWvSQdSgj6CNHBRIDgrknzLjmfIUvFs1t2Z4zeiMJxYuCFKikJxkXiyJaTRtB/+gRkaSg4IoYx0oCK308B3duzoOgyGzvZsvg7in3PrNO9VMoYYRxjFRz7QTbDJ7zfA1qfWpGT+doiA0erxbtwXr6x+DRgkTP+0b4wB/EzCmTrEaNar5399r+N1Sj2X8skOFh+b6zrVMsTUuwY+oRf440J3eC8NsrlwluLiwXjShvRZHgIjXoPwWNOrmCHOKChYibbgSddd+o5NnimlqGBLXTa06wTehUlBphg5F7jFQyWCEXg0rxuTJoNWasi5ZxC1qkbVEO0fr26HiTAqzkcFncIPWTe3bJyu5d88qMlLHgomVNTiEDLwJQwAZ6zI2tS61XerVCArOAZ4EvvH6trMSEuT5+fLCIhP/+rmyc0qNMsHvGUErk0lCfGjUFmVRmRkc86LLl1UNDe0dRR5TIoRoG6EIUwr11GdM2IpiGri7vcfP6RcqpDhY8pIlAkWLAhZFYyNcuM48cqvL+C6clNVpkHscot4Ea5//XyYMABw++lArFrf85z9dfYfP0z6nE+hvhr5p2ZgRofZiexKUiJ7WPNPjFF66RKqrs3vrTYRgWrlp+EogUg0sTpMj7SgIvbPfCNNCTWqjiITMnOzddbdJHpNoRNrpEpM/2jl/gNMIcNRru4keHiRfX1NqXG/V3MIhuBgX3Ymw5uxSNF0lRHgB707r22eqrK5WVlfTovWYTp84qEkHhWiAZ9o413ESlSSjMcOYq4NB0XnwjgOSXpEUxowZCIkkOG3SM21P5A+3G95DbSthYZAWgfMnTRbawDDu/v0kb29aTAcNIXuGDwkhlx7Q44ixPaqsYyDrN2h7VsZ48XolU6y1HmlrRbFEtwyH4F4tm45hWFpEF3Z1rQYuvQdMJ4h6u89CtL+RCSP5+bFfeol37Li8oIML6Wzp2eym7LfD3jalkEIf5iwNRDA48XtHRgxTKlt+3KlwcuxZjm8AZCaMWAp5p0CkO2HhNSNS3KXf6PzVcroCc41y6BG0kfHkaZ7TrlVdk2MmUMvW50BjHoS+bMR7nTRJmpWl/pNgUq8Jq74Vah+q8zbKWvj846UaTOmyZrQBWWxLXbAUAKC+oMeEeceBRgVPUwZ4pkU4RljhrW7X3Db+CPg1EDjD0JvMYNAnTBCcP4+ZwFVVra6uFFS2J/K7wtWJpnAmE6ul9U0mlbmK796TFxRav7qqM1nuOHoIn8yl1lL0LY6XOMwGHAFufd7+X41W+/jiUzEeXloUaNk9rLmRZad2zGMU/NR6oKjtTyc087/QmAeTvwai1f9HEwYANuvX49nshs2bMY2mXlz/3cPvRjmOmu3bK6p4fx+OyotqVSvPLWgFgLYDB5TV1S1TpwFqzs2JXAOYBtL1LE6jCBpIpAO79s4jAHhys1aCh/mz/Z/fMcE3Qa6WP1KaUFqZ/TsQrSB4nuG96JPiQaMR3TBE0Fwrqi3hlYx30yEwjKk1FTvukMHKao6LuayTJr1sKSlENzeim5vuza4RQKQOfCxJwpGinKPu1NzBwGCTY+E5wBHA38injjVvnkYgEF03zpH9UPmQRqRN8tCxODN1nj+KwXHTeKJa9+4lODszutPwOS4IxaG4qj9014vICSwYswGenGlvTT16vJAhxzzjnMkkS1wwtUIlvd4gxUTR6xdyyJyP73+swTQgrIdbX4D3+N6LnP6NTRiOQXf45GN5QQH31/9uSd0CCHwe9bm5Wfzn8dLyYBUClw4XqhoaW/fspcfHS33NDNTZ7hA4AzL3g1yoc7vn/A59o8s3K5liDSfMhkjUcfMDrQODrIOylFnGAjAJ5J2CoDlANtKWSPL1JXp5GWbgu1VzCwDGu+owYQV7rjA11vIQzHH0kD5/mphSKU1P1+uCAQCOCB4vDDx3GADEusY2S5sLuAVGTJhnDFCMRADUURFENzfeESPNzEKl8InyyXTP6WRdizPBAdYSBxJaJjbabyTNyJDl5FiveqVHS5x1oLuAzWfw2M3ZegLSqLeA7QaXNorFsvqUJoEVMn+Wv2V3r/jna1SESRprzWbafBT5USG3cE/uHrjwNmhVMO37vn1S/WvCDIjdFhcXHz58+Ny5cwqFwqwx6RMnMiZNatq1s/5Jxr9G/ksfI4VZcLSjEoOZDK7q4cc/YRhm//4mS0aJfhfkQkjXTfre2W9UeS5XiofFi/T2aiT4JjRoGvK5Br+3j46AQgThJgnKMyZNkqana9r0JoNvVt8M4AQ8z4hQm/yYXkvjk1v9lozrj+khzcrWSqVGBA29xwO3HHhVA2zCYlxicAjudrX+WLIpH7jl7ew6xl4ylP3SS9KsLPmTJwb2Ol9+Xg3qron8Hpi+KACHgVEtwda9e/E2Nqw5Oih6fFbHqTB568lC0ClxRKDApK+gKf/qf45RNRA51xuxyDxw8yvptVS+Fdd96igAiHOLm+M755fcnx9W34bxHwPH829jwgyI3T58+HDVqlUVFRUXL14MDw+XSs0jIGx5Y7YUp3n/NmOW5/S+OtuXl4eoQFUm9eKsfo3g7GyRIRwKAdMgbbc+R8xrYZQaU4dopXaRdhSyXv98qtdUIkI8VXJKv/eCQcY+cAkHF5OKs+jx8ZhGI7qt+23kyXmPmh89H0XKuELZxTo5SHzemthPM0SSkoIQiVajRhkxYQADH0uySKxQ+9Cb1fodwIIkQHEQYFLClDl3LmplxTuil10eA+x48XFXnGsAR6/mdLC/tcSJhKvotvTU86k9fiy5/4CzcoXOemwyi6YdRmRgnJKjd3QfHzBN6DqdW20r42jHR7tZ4llrsabDj9WY2nPVM+f6gyGveKg1/3J04g1/sc+fVH+ZMMNit+Hh4ffu3duyZcu+ffvYbHZ6uhm9ga2y1nfzPj83y86+rK11776+OmErjWJodZKI7n6ZMsbyUca+DzIBpO/VuZFszaqUN7hZOc4cayipRCPQhhGHXai4IFTqmallN6C1FEa9YeJJkQMDiG5u+mLJWzW3tJj2+SiyfOctEkJhJHiSWbR+miTilBSrsFAjlXc2vsBy/UtiyXiP+DJ+WXuJue4o0n0MUE1S5cAx6MxZswQXL2p4ursv0urTKgWVo0hGmCZnvxSIAnLykN7wtnX3HhyLxdZfj+27KFaItuFyVTKuQOcO/61eoQXSEpeDlpGdFe+/ztByNCPwVCebTqtGOf/Ot1yBAMW9m7zRbBaQv8qEGRa7Rf5cKJHJZJWVld7e3qZ+t1WSNTfXiJXiZW//lzl7Vsvu3dKsrD454aYvv3SquicjqVozuJbpg/zpiE2FB7tArmN+HDlV+FRBQxFc/WkjJnsUcZRcLe8hlNtlvu8GuqNZkvf0SfGSBw90drpcrrzswfTw53TLepQcuM5SWIu95N3Y6foU6qYmRUkJNTra+K7e46HyLvT11DduwtzjUQS9UnlFl6koheZCCDKDT5W9eDGmVPJPntS59WjRUQ6ZE0wwUoHl78ORu1LI1dInxTqWmGWPHonv3LFeuRK10r/ehyK284OIKKlstw7H9lZKNaEJZTo3OzQnQdZv5t6xlrxySimBT+D6dpX2SP0Rym74j//806jPMhszP0v7rG8fE76fHr9hKdx2aLXa5cuXr1mzxq3LgtTzUrid0IDmoPhgpbpyKW1pfnI+GhbmnpJavnZd9Yb1aiq1N2fLyM11SExsi4ujhShVmYRd2++Fj7ZQjJaJRcbKLxX9/naxQzcTo1BBXTJNQaA2auus623OnzyjJeoV6HTAOXjgPPZn72eVsXosVtDldePLbxc6zSu5cMn0syKRSO4q1d3t24Xdi+BFmChDkDGOPK7rxeKaJKRyegvWUG1HLO2FLqyRe56Z6YBhD9VqxXN/osedd+IzRsoF947taKP6DoDl6vrXPXGep/NPu1a79jQljUn+gFyrwsvrzbg/zj4+9b/8msrhYN3rDfla/h3hnbHksXgEb3TWOXqizbXUk/uyy8f2/NC6/PwLiU5PZbMwY4OwiHxvide1b3+X+T9bizh9NqnoFpWAQ2z9qM01wZxL/0quUIlJDqZenkrjk46SEEqFn6L8zxPgiEuiy76qZ4/KrONAnTqWHHu29GxbTdtEct9lJ7DeQaVSHThw4NcuuHjxIoZhKSkpEyZMaN/nl19+Wb9+fY8D2+3Xpk2b9I3c49ykKulr114LORCSWJbY+aMsv6Bo2PDKRS9pFQqLL0FeVlY0IvTpS4u1ajWGYdu/Sdu5+uaPu89bflOOvoR95YpJuF1/2/FDxn9W30x+UFd373HN+3cLf7lqYIDExMRrT68FHwi+VX2r57Yzr2HbHDBJq3mnpNWWjo+rfv2NHj//UfBH8IHgSn5l5y9KsSz/3aSyTVfFdS1Yf6L2rbdLYsZiWu3z195zVykP+5SN3f4S63/0+OunSk4FHwguaC3oud/u0dj+KeYOLklPL/APaDt6tMfv32d+P+z3YY2SRh3Xrgu792bvWn3z0rWKrj+KU+8X+Ae0HTpkyggaparg/aTyTdeEtc2dF753X/au1TfPnC/FMAwT1GHfemE/jcIUYhOvLu+bxJr3kysvpT/7iV+LbffDdo7A5MKOaYhpP73/afCB4B1ZOwy87Gaht4EkHo9ftmzZK10wdepUMCaFi2HY6tWr2Wz2N998Y5KDKmtZdW1VekP659Gfz/J+5t2QgwIdv/lalpPTsGWLZaG7Riis2/AmSqE4//B9O63lqtXDZXgQ55PFUpWFNyVuCyglkPzs0kor+FAkFDuQYiKdnKJDhEgbrlRtWN9ovNt4B6rDkcLuK/G8p5B3CsJWgpW1eaeEIPT4iZL797XibgIzlysuB1oHejA9On8p+vEKHc/Cj2c/y2X0AzCNRpKaSouONomogMIC57CBz+gDwAS3CXgUf+Vp91iSWwZN+WYF8h351ogIyogR3F//i2mePXqFRnGm9EycW5y9lb2J46xYHiImILnnnyqV2k5PpGXHDoKzM+tFk/LlKAFvtyQYB/jqPffb52FzG06ewxOw8XOm+wAAMJxg3m/QWgpJa0x5s0qP3mG1cfjWPI8pf0oTKcRwZD6opLDwSGf3AgLIltFb5vvN/zXv149SPlJoFL1/Rv2VCzMshZucnPzrr79eu3YtODg4ODj4ypUrBoZKrk2ef35+Ka/0h3E/dLVfHfHIpEm2b24QJJ1r/m672S+SXF77xhplba3zjz/i7TpS7EwmaegcL4Ya/fmnHAsv3jYAQl+GzP8CtxwAMC2c2JerRWDhKyEdljfSjoLQy04bKjrHIbgF/gvSG9IrBF1IX1N+ABQPY9ZbErhNmoQpFOI7z5aiakW1ea15Uz2fralVnE1lC62f4ivdJ43sV9Mgy87WCIW0WJPpXr3HQ12WYWa+/gCTxBztNPrK0yvdalwLzgGCGi7K17s4sfo1VW2t8MIzHeVLFZf4Cv7CgIWmD0Im4QOnutGV2M/7Oqao8OpV2ePHNmvXmNoSB2Ad7KUIAaaaU/DjRbVa25BD1iDw0pouzYyeMTDxU8hPhKsfGh6qMb2Q9AgEuLagt/7sP1PJ4NhCaC2B+QfBttsaa7sVWzdi3YWKC8suL9O7WvKXmzAwKIUbGxur1Wqf/InJkyfrHCG7KXvtzbXrbq7jkDnHph8b5zpOz7RYzV68mLt/f/O335nui2EKRe1bb0lzcpy+/cYqvBv16OQ4j2aWmlAuvnmv2sKLj/0A8CS48QkA/HEsnyXQMCPtPFw7ylA9Z0SKga/NFmtVhlLUCX4JRBzxmbyooBZyj8KIJUB3sOCMKMOG4R0cuvZLXqq8hCDIZI+Omy+oqEceyETA44Ux+9s0iG/fQYjEDr0iE02YVgOVd2HAMd1zer24/mFjF83QgkRwGQkMS4SFaGPHkgMCWvftA40GADDAfi/4PZATONLBvG/GzCneInsiViDIzmvG5PLmb78jBwWyZpnnGPotGcdjcJktrHOfnbVWos5xTp1TtAOj18HotZC2G25/qW8Qfnmt+HSVAmSeb8Z06CUqpXBkPlTdh9l7wFvHO4sAsnro6h/H/dggaZh/fv4XaV/8j5owMEHs1jCWXVn2uOXxe+HvnZh+wovpZSBKctj8UbsVq3tvIyY33mCo4fOrV70qvnvPcetWxiQd/RxDwuVSIuQcK6ttEFty6jQ7iHobCs9Xpd3mpjTx6bgVS0OenS+KkMbYUBFG2fF7htawSOyZ3jOTypK4Mi4AwK0vAEEh+h0L7yaCMCZOlNy7p5XJ2l+exLLEcPtwB6oDAGhV6rpfHyII4rByBNZHMssGILp9mzpqlKGFs56Z8DAgM/+S0oo49zgGkXG29GzH/9sqoSHXgiiy8ynYrF2rrKjgnz4DAMk1yeX88uXByy15O9aOUCBwdX9+48/7VQ0N9h98CDicuYMEvjuNq20NV1iTyY0LEnSVpMVvg9CXIfkbuPQeaHumPkR1Lc0/P0YBz3rJz8qODQAgrIffJkPVfZizr1NSS/dX3jU2aXbSPN95Z8vO/u+asF7i86jPb8y/sWzIMjyKNzozHDZ/ZPfuu8LLl58uWiQvMsRtIs3KqpwzV5ab67x9uz46HTIBi3tlCEELB7/PspC5fPRaNdMj5UglYJCwOqRHobPn9EgRwkPy5GqFoYzb8uDlKq3qUOEhaHwCj49DxGpgOlt8P+nx8Vq5XJycDH9Sgyb4dtSCF+6+wtBy1OEEtp9rfz9W5dOnyspKWmysOfMUD15joey6xcqMFoOEI031nHq96nqHTnVhEgBiuQkDoE+cQAkNbdm1SyuX//bkN2eas86mSKNwtKN6TXaly7CMG82MKZOfydyag8Y2xW0JRaASxJEcam5k6XyzYMZOiHoLMn6B32cAv6Zzi6C8vn5nBgFIlDmOtsO8AQAKz8O+GGirgEXHDNuvdrBIrI8iP7qScOUfa8Jm+8w2nUgaEMT61VWue/eoW1qfzpvf8PEnz8tZywuL6jduqlqyFPB4j6NHGFMNsbyGD7O3GevAFGl2fp9pydkTKAdEn/CV7hHBBQE+nOduPGI1zpGC0MoOG6JDcKO7xbvHHys6JrqxGchMeOGd3txPq7BQvI1Ne43r6dLTndSgtbceMRoYfFqrz/yYAXis4tt3AIA2Lta8w/yngqAOGh8P/Dyc4ztHoVFcfnoZAKAgCVzCgOnSmwHtN76nbm4u+Omr7ObsZUOW4RCcZePMnenLkebXOsWkD19oweFSmfrHJs/BAAAb7ElEQVTgvzMRLUk6I0CqFWHXRWUn7+q2YhM/hbm/QONj+Gkk3PwMhPVNmUUt+x4TMRJhqq1jZBA8TYE/ZsPxJcB0gVU3DRCo6UgdU2z/sSbMkgAuJsbr0kX2woWCpKTyKVMrpk+v37ip8bPP6t9/v3zK1Mo5c0Q3btisXu11Lok8xHjf8uKFQXIvKqlCvO9ns1P7+w/kKepZrta5UbzP2/P6PeAeP1KItuGLtCqJoch3ZchKsUp8siULxm4Ccu+yVChKj58oTk7mC5puVt9spwaVNvPkVxulmMjvzfiBeUai27dJ/v4EJzNzSb6TAMVD0cWBn1RB1kH+HP9TJaeAXw11Ob1xwTrykiNG0OPjsQOnfKWMOT5zLB5HeOHC8If7MExY80Cc/KDWrGPVau2OLx4wxFrHeOdR4/zLQlUiHJ+cBU++PaeW6VooHPoivPEAAqZjKTuLtm5TnGzAMA3N/4ZL03fwQxAcmAaNeTDpS1h1A2z9B/Lp/NNMGADgGAz7zR/53Lpp/+EHBEcn6aNHwgsXpdk5JE9P+80f+STfsX3rTeN00n/irfdGCm2Jqmze/gN5pp/DibPFkrQWIRs/deMiIJAhaa3O4nLGJE8ySi39zVC5QKCV8xgldojNUYQt6/3NocfHa6XS5NM/dlCDarGKXXeIQGbN9yHSqQPwdDRCoSwriz7O/L5xKw64RULxxb9kUr3o92IBtyAnax8A6KOZNguStQtUKPbuHToZT7ZsBDWX27jtS6thIXO3TpATIOuPkuy8ZtPt1/ZtD5hcFSnC5sU5/gCAWRH8t0zhc3isNlbF1ptlp+7pWGtiuFTRVhdJ/0sjLhRqq5zZn9i2HoXmInCLhIRf4a0nMHotoPgBfjR4+IcCb23Nefllzssv99YgoujbmyN/2HofSWvZKcxcvy7caPv+wSP5grtNIhq6fnMknkqEqf+G069A8rcwrufitMvYoU9uJdKqGaLaZrqLnsbJm1tX8bgr7W1OlJ5ZGrS0l5dDHTkSx+GIrl6LWBnhw/Ip+vkKS2kj9pe5hQ/Ql1Ny7x6m0ZgdRXbGklc/BH41sNwGeDrN9J65M2fnoZqrI5yG98lf/6nuGHM85aXLVYLz55kzzK/P0GrrN27CZDKnL7cRHWlz3wpN/Hf27T1PxC8HxEQacW8FIuVPX6Uz21RYMOvVlc8of/EUUvCmWTU3srQ3BFaZUJlxU8aS4R2oBAZZLVGqmsRELsEKoROAJB+uDX5xGaDL/xfedBQGYQxkEv7dz6PFTiRcgfDrj+42NEsMJBe++ypNdLdJxMSv/3RMh2hCyDwYsQTubYdKHeuPzovDUAStOqBHuLD4EjzcP3L4qkjHyF8e/yJRSXp7MTgcf5RfYJFkud/ixvRCq3Iyn8QNWBk/YDdTdOs23tqaEhJiycEB0wDgL4klyXhyguuEGyCv85/Q+9EKuAV3au64LnuVMmJE42efq2przR2hdc9eyf379h9vIXp5AYCvF2vG28PVOCTnQOH+3w2FC5m5TXs2p9LbVKRRNuvW6ZDRc50Q5v/ldHUUQUaRWQmotBIyKROohUQal6EiKOTDNR5bx/ksHAv9v2w9aML6EkQ8umlzFCGMQ+Wpj25N37Mvp7W7DLhMrj56smjXpnuUKqnCh/retuhuoi9TvgMbPzixFFpLe4zM9nMV2gmZYuumjOdWUfnVkLgGnIZD3JY3Q9/kK/jP67xbgCTXZisFBBWrhWcqFZjMe8O4AbuNmEoluXeXFhtrHiPus5vlAfZD/qpYcpECEIBjeE3vh9qZvZNFZi0estR5+3cIAnXvvIuZw5onunmzZfdu5uzZrLnPdGQDfDgvb4mQsAiyBy3bNt65ebe6hwJnaQV/+zcP0vbkE9SY3wLvVSv0Si4gKOIxIzL409leX01grPGjLHNjrPFz++KFIV/M8VkY200O8n8h3ho0TyYCQeG1V4dnRzZfOVREyuEdfvRAQseROGQEB3KegsxXk7SgISP+c70njnN/zgRawaLj8N8JcGQ+rLwGtG4xo88r4+u+ui9KLLYP93/2cZPx4fA8AID5BwBHDLYJnuA+4ff83xcGLOyNUEBGY8ZFdvUiGqX2Sj2H7Imbad1NI76/o8gHDzRCET2+F06f/1RI+QFkPKCwB3gCOBRenkyjnnh6aWXYBjbZ8r/+sPFhan3qu+HvUglUcKY6fvVV7foN9R984LzdpPYS2aNH9e++RwkOdtz6SY9NjnbUD76MOXD4iSytuehIWeaJMi2HSGYQVTK1mqdkSLUkAJkb5eXXhtnbmFaRhyIMN/v/8Rdz0AszD6Ehdh9+EzPytSEaTxqoMaRaApViVKxR2pM8Ejw3fT9Wh/3q8CDcYdExEDfBgWkg7CbOTGbRNEMJDC2n+FBHXp+glsCRF4H3FBYeBrZH+4/rRqyTqWU/5fzUm/PflbPLhmavil5sTfYWuUucooIH8u6Jrl3DMejUMaMtHyJgGmjVUHJ1oB98Yx60FL3hu0Cuke9/st/iYbSY9uuMr51pzosCFrX/Qo+Ls9/4nvDS5cbPPkeMVb1JH2ZWr3oV7+Dgume3TlJDBIUVS4Nf/+4F1jhHLZOAtimhQoxrlAGKQAhzxkfhGz8cbar9+ptg0AuzBBGh9hGh5n+dnMNgyVk4Mh9+jYP5B8D1Gb+d36LY/IIkaj5dWNPMoEmjy74EVTPM299V7tiL6bUwYOHRoqMzfWaG2FiSS0qpS8lpzvmS9q41yaNVUOQTNKCL35hGI7p5izZunOl9fDrgOBwYTlB8EYYtHNBHnncKcET34cuna7lHi44uDVpqZ2WJGMrx4uMlvJId43Z0LXjkrFih4Qta9+1zLCjQTpyI6mGO4p863fjFF0QXZ7f9+3EcjoG/QrMiLF4QCAv+X7yMg17YwMItElZcATwZDkyDqx+CtK3TY3dYMgIFXNWea7A3iqJqg8WnIKAnrfa64es4ZM7XGV9rMbMbBrSYdkf2jlBNwJB0RylISLmHhBcHNKkkTc/Q8Hi9iiIBAEEgcAaUXgelZOBOHdPCk5PgEwdWnNXDVqu16p8f/2zBMFwZ96dHP0U6Rsa5xfXYZPv2W/b/+hftyZOK2XNEt271aEJQVlbWrlnbsHmzVWio+x9/dFISDGLQC/srYD8EVt+Fa5shfS88/BW8xoHTcMARrbllTRo2G5lSpl5SERgY76mjUJ5GpL0T9s6HKR8eKTqyJHCJud//6qbK/aWbUcDbvuynsooTJCZqJR+jVOrAXLfo2jWUSqWZQtNqGEPmQPo+KL5kSgtL36D6AQjqYMJnAOBGd1sUsOhI0ZEE34RAa/MobT998KlcLf9wlG7iB87yZZkiofely7Vr1hLd3anRUQR7B41YJMvOlmbnoFZWdps2Wi9fbuFKyKAXNoi+BIkOM36ENx7AyFeBWwrJ38LNz6D8VkBEkxBrwYknKsV6K29neM8Y5zpuR9YOs1hKWmWtu3N2f1+2nolYY9FW1kM8mbNnaWUy4dWBSippNKLr12mxYxFSr9ezXCOB6QxPzgxoFEmkdupFrh2+lkPmfJH2hVm+8Lnyc7drbm8I3eDJ1CvhI3Nz87pw3um7bwmuroKkc83ff9+2/zetXGG7Yb3PtavWK1cO2q9BL+x/Cbb+MGkbTNrW9XviOKJW8GuxWwEeU2sQvO7WuU/HfDrn3Jx/3fvX4amHTekhxQDbkrrlzcczPMBb7C0NmBkDAJRhw4ienoKziV0X5vsxiszOVnO5vY0iO2PJoDnw8BeQC3rbcWUK1HLIPwsB0zsFqGlE2saRG9+/+/4veb+sHrralDGeCp5+lfFVmH2YUd8ZweGYM2a0F7tiGg2Cww2+KINe2N8JbF8XdTiBg9rl79CbqGKT2duit5XySj9M+dCI3DQAAPxR8EdwMiMSG83n8AJefUaKwJw9S5qZqaqrG4DrEl68hFIotJg+aiMPTgC1AoouDMQjKbwAMh6M6GZ6pnpOneY1bc+jPY+ajeuui1XiDbc3EFHi1y98jZojzThov/56E2ZACrcdJSUlWX0kQfSPgc/8mCq0ktXKLtUn9gcQ5RT1Ttg7155e++7hd4at2K3qW23HH81STeORWoe8162RhTlrFqCoIOlcf18RplIJr1yhT5xoemuqETiHAsdrgGLJR4eA5QYePVN4m0dtdqG7bLi9oUpoSKZXppatvbG2RlTzfez37dRsg/jbmDADUrjt4HK506ZNe/fddwcfQw+0htMFOC7pEeimcAIAgGVDli0NWvpHwR9fpH2hT5svuTY5d//FeYrZfAI36KPpSPemEIKDA3XkSEFiYn+TcInvJGv4fMaM6X056JDZUHEHJC397D3WQ+VdGLEEnvOeaETangl7EEBeu/5aOb9c59Ft8rbXr7+e25L7bcy3YfZhgxP772TCDEvhtuPtt99ev3794DPQARzq9c44KSLSXBfWp+hlyNo0ctOqkFUnik8su7KsG8U+gEKj+OnRT+V7bs6Xz+LhWwM2T8cRdVRjMWfPVlZXS7Oz+/VqBOfP462tqWPG9OWgQxJAq4aCfnYhc/4ADINhL+nc6Ep33Ttxr0qjWnp56aWKnpp4KXUpL55/saCt4Nux3050nzg4qfsJ/ZXO7yGFm5LSU+ri/PnzDg4OoaGhZ86cGXwMz4NizXTdMLr2x3TsvLZK/NB9sm5u9TdD3wy0Dvz0/qdzkuZEO0dHOERQ8JQKQcXN8hvrH00LRSa2UVuHvD9Dp/0CAPqkePSLL/gnT1qF9ZePoBGKxMnJ7AUL+jiz4xAMDiHw6BCMfKXfAmAMco+CVyyw9DLZBnICj04/+s6dd96/9/4veb9McJ9gb2XfKmtNrk1+0vrEk+m5K25XACdgcD73HxCsf4KIs2fPHj9+/NixYz3+3Q4+nz9z5swrV65kZ2dv3rz5ThdZna5SuImJif/fH49Q4f2YyMCxS2iV4hC9fXkSTJKmSHukfMTT8gDAX+S8ueoVG5x9JeFpWyjDMKmAXWIi82Fm5YcfqPunQIyZnmF/5kz1hvVyZ+e+Hdmr5XpI7R+3A7YJKf1Clm0rKhhT9nWmx5o6dqQRWwdYrjI3XZFep6nTghYBxBHnGE4KDyeGo4MrZiZg9uzZlhuiv0QKd9u2bRMmTHj//fcXL17s6ur6008/GZXCHWCYKEo6AH9dzhPlfXS25v27T7ac7RQu1Qe+qC1nV2LlplsVG6+XJ6WapARcXFzgH8Ddv7+frr1y0UvlU6f1y52XcLHPbbErH/TXcz+2GPvGE1PJzJBVVstrhDUycw75H5l1f+2E783L3ttAsl0K9/nfO6Vw2Wz2hQsXVqxYAQBpaWlhYWEEAmH+/PmhoaHt8ebjx4/b/z0InSCxaMGfzSr89RqtnMn9Ma/aSeIyO5zp4dhjNxlPWHUmDVeisUE4AhLXZXkE08skZmeSnx8lNJR39Bhn+XKThGnNgaK8XJaTY79pU7/cGisO+E+BxydgwqeAI/Tx4IJaKL4MURvAHFZVEo7kQncZnLT/hFxYpxSura2tk5NTpxRuUVGRg4ODr6+vr68vANBotPPnz0dGRg4+CUNAkcDXJnELqxqOZzMbrYV7S+uQDC0bwTFJmFaLSVTAw+hqJg2hCNE24ljSkMnmMSOzFy6o3/S+JC2dOrqPHwT/2HGEQGDOmtlfd2bEUshPhJIrlgnTGkLmfgAMwlYOzr7/pyYMAF566aV58+ZJpdJOKUk+v6cac3R0dNdE2CAMwDrQ3XqrO6+0tv5KLtIANC4d34YHAC2mkSJCka2QE+0TFPmCBSMzJk9u+vob3qFDfWvCMLlckJREj483TKvQu3zYOGA6Q+b+PjZhagVk/w5+Uwwk8gfxzzdhAEAkEolE4uBd7kOwfV3Yvh2hilIo0WoxMovWyzERIpG9aGHrnr3Kiop2IuM+gfDKVY1QyF7wYn/6pzgIfwVufg6tJWDj12fD5p8BSStEvDY43/4GIcrgLfj7gsig9t5+tYOzeDFCJHIPHOjD02s7eJDk42MVHt6/dyF0GeBJkPFznw2IYXD/P2AXCJ4xg3Ns0IQN4u8BHIfDmjNHkJikbumbendJWpq8oICzou+XCHqCagPBc+HREZAL+2bA0qvQlA9Rb/X7mQ9i0IQNog/BWbEc02h4hw73jQt24He8tTVz+vSBOPWI1aCUQPbvfTNayg/AcoXghMEpMWjCBvG3Ckvd3BgTJrQdOYKTyXo5lKK8XJyczF68uA/YwUyB03DwioUHu0At7+VI1uJiqE6DMRv6vkpjEIMmbBD9DZv167QSCTs5uZfjtO7eg1Io7EUDSG8fsxFEjZB9sJfDBDacBroDjFg6OBkGTdgg/n4g+fgwp09np95Xt7Za7oKVlQkvX+YsXYpjD6BOmkc0uI+B1B2gUVo+SOk1a3ERjH0fCJTByTBowgbx93TE1q0Ftbp1zx6LR2j5cSdKpXJWrhjoU4/ZBII6ePirhYdjWrj5qYQ86IINmrBB/J1BdHMTREbyjh2XFxZZcLgsN1d044b1smU4JnOgT917HPjEQfI3z3ShzELOH9D4pNBx3mAWbNCEDeLvDe6keDyL1fjZZ2ZTIWo0jZ99jrez46xY/tec+qSvQCmBO1+ZfaC0DW58Cu5j6pgjByfAoAkbxN8bGjLZbuNGWU4O/+Qpsw7kHT8hz8+3f//9ARN26wlbfwhdBpn7of6ReQde/QAUIpj+w2At2KAJG8Q/AcxZM6mjI5u+/lpZVWXiIcqqqubt26lRUYypU/7KU4/7GGh2kPi6GQUW+YmQewxeeAdsB7kJB03YIP4ZQBCnr79GiMS6d9/D5MZtAaZS1W/chBAITl9u+4vPnMyEmbuguQiubTbNdayCC2+CawTEbBx87IMmbBD/HODt7Z2+/kpeUFD/rw9Aa1DzFcMatmyR5eU5fvEF3t7+rz91nzgYsx4yfoEsY/X6Mj4cmQ8ICnN/BXRQU3XQhA3inwVabKz9xo3CK1caPvkENBp99qv53/8WJCbZblhPnzjhf+XUJ2wFnzi49C7kn9Vvv3hwOAF4VbDwCLDdBx/33/VbO3gLBmEAnBXLNSJh6+496lbu/7V3fjFNXXEcP6UUWigjzsK0pBUmhpASSjSIoiINiNURNRKm0pm5EVEJgeiTJMuiPEg2DDHxwQczWEbKEhibk0TNQqBRwp+UpLLIGqIZCUFLS4s2QSptoXu4syFAr1fovb2F7+fp9pxyf+d7fuf+uPf03POT198Qvt/6jWJhdnayrs55789NujLZpUt8+tcsJKW/kF9PkY5y4pwge6uWztPbzKT9a/JmnJT+TJR74WiEsJUxGAwOh6OwsDB+pVVCL1686O3tlclkBQUFEgnWQ/OUhOrqyIQE6436f78o/vTbb+IOHYpKSvJMTc109zju3vXYbAnV1bLKS7xrd7SU6H4jf1SQv74j/9wjuTUkJY+IJGRqlAy3kqFmIv6EfPU72ZYLFyOErUxNTY3FYsnMzKyrq+vt7Y2Li1tc29bWdv369bKyspGREYVCoVar4QzesunMmZidO60//GhruGlruOkvl6jV8psN7CVwWysiCfmyhTxtJYYbpG3RmnthFFGfIQXfk1gZnIsQtjJUKtyxsbGIiAir1arX6y9evOivdbvdtbW1T548kcvl8EFYEJ2Wpmz6yT0+PjswMP/GGSGVxuzaGZ2WFgZNzyojmafIeD+Z/Jt450i8gnyej+CFEPYB6FPhms3mzZs3NzY2mkwmjUZz9erVyEjMyoUBUUpllFIZfu2OEJLk/SR5Pzy4/ghNKtyurq6ioqIHDx7k5+frdLq8vLyamhqqCqlwAdhohF8q3OHhYblcTh13dHSUlJQgFS6swzpS4YZNKtz09PSoqCi73S6Tycxm87ZtWJUDAODTXBh9KlyRSHT79u2CggKlUmm1Wjs7O+EJAACPQhj5UCrc4uLiI0eO2O32z/jwSgoAACFsOfSpcIVCIeIXAGAt4B1JAADuwthBgP3nAADhG8JYWrPGMHrCOqxvHOshF44HSQDARgQhDAAQxgivXbvGz5a9e/eusLAQ1mEd1iGc7iE0hA/AAACAB0kAwMaFR79Iut3uZ8+eeTyenJwcfyH9vq/0tatgZmZmYGDA/zE1NTU5Odn/8fHjx263mxAiEokOHjzIRifQmwi63sWMjo4ODQ3FxcUdPnw4Ojr6oxq2ajj2b8j1Mj8/e9rpBzl72lm6wHn0IJmdnR0TE+N0Op8+/T+PqX/f1/b29uX7vtLXro6XL182NDRQx21tbY2NjadPn/bXbtmyhfookUjq6+vZ6AQaE2zo9WM0Gq9cuVJUVDQxMdHX1zc4OBgTE8OwYauGe/+GVi/z87OqnX6Qs6edrQvcxyfMZrNaraaOJycnlUrl/Py8z+erqqq6c+fO4m/S164d6uVNl8u1uDApKYntHghkgm29CwsL/uMDBw50d3czbNiqCa1/udfL/Pxsa6cf5KxqZ+MC5+9c2JJ9X00mE/PatdPS0lJSUiIWixcXejye/Pz8EydOdHV1saQ6kAm29frXFrpcrrGxse3btzNsWJj6l3u9zM/Ptnb6Qc7NOA/iAODv6nyn0xkbG0sdS6VSp9PJvHbtNDU1NTc3Lyns7OxMT08fGRkpLS199OiRSqUKuupAJtjWS7GwsHDu3LnKykrlst2lg649tP7lXi/z83OjPdAg52acB3EAcB3CvF6vXq/3er3+kq1btx49enT5NxMTEx0OB3XscDgSExOZ166xGUajUSAQ7FqWmGf37t2EkD179pw8edJgMKzatTTWA5lYo14m1n0+X3l5eXJycm1t7fK/CpZ2bvzLBI71Mj8/B9ppBjkH2oM8AHg7F/b27VuFQjE9Pe3z+Y4fP05tjNvf3+92uwPVBosLFy7cunXL/5EyOjc3R02geDye3Nzc+/fvB13+iiYo66zqpeaGzp8/f/ny5SXl7GkPpGi96mXoaLbHNv0gD9Q2Pl/gPAphZ8+e3bFjh1gsVqlUra2tPp9Pr9dnZGRoNBqdTuf1en0+X3x8vMViob6/vDYozM7OJiQk2O12fwlldGhoKC0tTavVpqamVlRULJ4PDhYrmvBLZkkvRU9Pj0AgUL3n4cOHHGhfUdE61svQ0WxrpxnkgdrG5wuc76vz3W734n1fP6qWjcZYrVaZTMZe7nF6Exzr5UA7r/zLpa8/eP51rD24AwAvGAEAwhi8YAQAQAgDAACEMAAAQAgDYQwmZ8FHEYkuADzh9evXg4ODUql0enr62LFj6BCAuzAQTvT19Wm12n379hkMBvQGYAgWVQBeYLFYXr165XK5zGZzdnZ2VlYW+gTgLgyEDaOjoyqVSiQSORyOlJQUdAhACANh9TggEIjF4pycHIVCwd72MgAhDABWsNls1IHJZMrIyECHAIbgF0kQesbHx6empoxG4/Pnz7VarUwmQ58ApvfvmM4HIaenp0ej0bhcLg7eKwZ4kAQgyFB7HyN+AdyFAQBwFwYAAAhhAACAEAYAAAhhAACEMAAAQAgDAACEMAAAQhi6AAAQvvwHJ68WM++hvMYAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "\n",
    "for n in range (-2,4):\n",
    "    x = np.linspace(-11,11,1000)\n",
    "    y = sc.jv(n,x)\n",
    "    plt.plot(x, y, '-',label='n='+str(n))\n",
    "#plt.plot([1,1],[sc.jv(0,1),sc.jv(-1,1)],)\n",
    "plt.xlim(-10,10)\n",
    "plt.grid(True)\n",
    "plt.ylabel('Bessel $J_n(\\\\beta)$')\n",
    "plt.xlabel(' $ \\\\beta $ ')\n",
    "plt.plot(x, y)\n",
    "plt.legend()\n",
    "#plt.show()\n",
    "plt.savefig('bessel.pgf', format='pgf')\n",
    "print(sc.jv(0,1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAoY0lEQVR4nO3dd3hU1drG4d+bRpAuhBp6k1CF0Kt0RIqKKHZFsWFXbOfzWM+xHVGUIiiiWBAVEBXpIIK0gAoHEAy9SuiEAGnr+2ONhxjBhGQma8p7X1cuJ3t2Ms+G8LiyZ++1xBiDUkqpwBfmOoBSSinv0EJXSqkgoYWulFJBQgtdKaWChBa6UkoFiQhXL1ymTBlTrVo1Vy+vlFIBadWqVQeMMTFne85ZoVerVo2EhARXL6+UUgFJRLaf6zk95aKUUkFCC10ppYKEFrpSSgUJLXSllAoSWuhKKRUkcix0ERkvIvtF5L/neF5EZISIJIrIGhFp6v2YSimlcpKbEfoEoOffPN8LqO35GAKMzn8spZRS5yvH69CNMYtEpNrf7NIP+NDYeXiXiUhJEalgjNnrrZBKFZTU06dYNfklMtJOkR5WiLSwaE5GluJEVBmSo2JIjiqDEWe3b6gg0aVeORpXLun17+uNn8xKwM4sn+/ybPtLoYvIEOwonipVqnjhpZXyrtVTXqf15uHnfP60iWCLqchGE8sGU5WVmXX5r6lBmrt79FQAKls82m8LPdeMMWOBsQDx8fG6sobyKydPJFNz41jWRzUkbtg8SDtpP04kwfF9cHwPhQ5upl7Sr9Tbv4H+R3+0XxgRDbHNoXZ3uKg3lK7p9kBUyPJGoe8GKmf5PNazTamAsnraG7TlMAcvGQURhexH4ZJQvAJUaPTXL0hOgh1L7cfWRTDn/+xHmbpQ/3Joci2Uqlrgx6FClzcKfTowVEQmAS2Bo3r+XAWa5BPJ1PntXTYUakS91pfm7ouKxkBcX/sBcHg7bPwOfv0Gvn8Zvn8JqrWHpjdCXH+IiPJZfqUgF4UuIp8CnYAyIrIL+CcQCWCMGQPMAC4FEoEU4BZfhVXKV1ZPGU4HDnO086i8f5NSVaHVnfbjyA745TP4+WOYcjvMeRpaDIFmN8MFF3ott1JZiatFouPj443Otqj8wdHjxzn1n4YcLRRLncd/ABHvfXNjIHEeLH0btiyAyCLQ4nZoe78Wu8oTEVlljIk/23N6p6gKeQlT36Qch4nq+oR3yxzs96vdFW6cBncugbq9YMmb8GZjWPgSnDrm3ddTIU0LXYW0w0ePU3/zeBKjG1AtPpfnzvOqfAMY8B7c9SNU7wAL/w0jLobVH0Jmpm9fW4UELXQV0lZMHUF5OUh01ye9Pzo/l3JxcM3HcPt8e4nj9Hvh3c6wS09BqvzRQlchK+nwMRpufY+thesT28zHo/OzqdQMbp0FV4yDY3vh3S7wzUNw+njBZ1FBQQtdhawVU0dQsaBH59mJQKOBcG8CtLobEsbDqDaweYGbPCqgaaGrkLTv0DGabH+f7YXjqNC0t+s4UKgY9Py3HbFHFIKJ/eHr+yH1hOtkKoBooauQtHTqKCrJAS7o5nB0fjZVWsKdP0Cb+2DVBzC2E+w768zVSv2FFroKObsOHqfxjgnsKVyHmIsvcx3nryILQ/fn4cav4NRRGNcZVoyz17Qr9Te00FXIWThtPDVkL4Uvedi/RufZ1ehor12v3gFmPAKf3wynk12nUn5MC12FlG1JyTTZPp6DhSpTKv4q13FyVjQGrp0MXZ+FDdPhve5waKvrVMpPaaGrkDJz+ic0CNtGZMcHISzcdZzcCQuDdg/AdV/Asd32vPrm+a5TKT+kha5CRuL+4zTZPp5jkTEUb3GD6zjnr1YXGLIAileCj66E5WNdJ1J+RgtdhYxpX0+lVdgGwtvdF7hT2V5YAwbPhjq94LtHYfY/dNoA9T9a6CokrN9zjCbbxnMyogRFWg92HSd/ChWFqyfa6Xh/fAu+HAxpp1ynUn5AC12FhMnffkfX8J+QVndBVBHXcfIvLBx6vQLdnod1U2Di5XDysOtUyjEtdBX01uw6wsU7JpAafgHRbe90Hcd7RKDtfTBgPOxOgAl97LJ4KmRpoaugN/HbhVwWvhSa3QKFS7mO430NroRBk+BgIky4FI7tcZ1IOaKFroLaqu2HaLLzQ5AIotrf5zqO79TqAjdMsbM2ju8Jh7e5TqQc0EJXQe3d75ZxVcT3ZDa5DoqVdx3Ht6q2gZs80wWM7wUHEl0nUgVMC10FrR83H6D+rklEkkFkuyAenWdVqRnc/C1kpMIHffSu0hCjha6CkjGGUTN/4YaIeWRedJldGShUlG9gJ/ZKPwkf9IUjO1wnUgVEC10FpUW/HaDWnmmUIJnwtve7jlPwyjeAG6bZ0y8f9NE3SkOEFroKOsYY3pi1niFRM8ms3AoqN3cdyY2KTewbpScO2lI//rvrRMrHtNBV0Jm7YT+V9s6hotlPWNsQOXd+LrHxcN3ndoT+0ZV2xK6Clha6CiqZmYbXZ2/k3ugZmAtr2TlPQl3V1naqgKQNMOk6nSYgiGmhq6Ayc90+SuxfTt3MzUiboXbqWQW1ukL/0bDtB5g6BDIzXCdSPqA/7SpoZGQaXp+ziQcvmIkpEgONB7mO5F8aDYQe/4L1X8F3w3RJuyAU4TqAUt7y9S97kKRfaVkoAdo/BZHRriP5n9b3QPLvsORNKFoeOj7qOpHyIi10FRTSMzJ5Y+4mnig2G5NZGGl+m+tI/qvrs/aKlwUvwIXVoeEA14mUl+gpFxUUpqzeTcrB3XRL/x65+Hq44ELXkfyXCPQdAVXawLS7YecK14mUl2ihq4CXmp7Jm/N+Y1iphYjJsKcV1N+LKARXfwTFK8Kng3QyryCRq0IXkZ4islFEEkXk8bM8X0VEFojITyKyRkQu9X5Upc5ucsJOjhw5RP/0mUi9vvY0gspZkdL2GvXMNPjkar1GPQjkWOgiEg6MBHoBccAgEYnLtts/gMnGmIuBa4BR3g6q1NmcSsvg7fmJPFJmORFpx6FNiN9IdL7K1IaBE+1c6p/fAhnprhOpfMjNCL0FkGiM2WKMSQUmAf2y7WOA4p7HJQCdOEIViE+W7+DAsWSuNd9A1bYQ28x1pMBToyNcNhw2z4N5z7hOo/IhN4VeCdiZ5fNdnm1ZPQNcLyK7gBnAvWf7RiIyREQSRCQhKUmXylL5k5KazqiFm7mvwgYKndgDbc76Y6dyo+mN0Px2u+j02i9cp1F55K03RQcBE4wxscClwEQR+cv3NsaMNcbEG2PiY2JivPTSKlR9uHQ7B5JPc2v4d3BhTajdw3WkwNbjX1ClNXw1FPatdZ1G5UFuCn03UDnL57GebVkNBiYDGGOWAtFAGW8EVOpskk+n8873m7m1ahJFD/wMre7S2/zzKyIKrvoACpe0c76kHHKdSJ2n3PwLWAnUFpHqIhKFfdNzerZ9dgBdAESkHrbQ9ZyK8pn3F2/lcEoa910wG6JL6G3+3lKsnH2T9Phe+HKwzvkSYHIsdGNMOjAUmAVswF7Nsk5EnhORvp7dHgZuF5FfgE+Bm43RiSKUbxxNSWPsD1u4ujaU3PYdNLsZChV1HSt4VG4Ol74Gm+fD/Bdcp1HnIVe3/htjZmDf7My67eksj9cDbb0bTamze3fxFo6fSmdYqe9hl0CLIa4jBZ9mN8HuVbD4dbv4dO1urhOpXNCTjiqgHDqRyvjFW7mifklKb5wEcf2gRKzrWMGp18tQriFMGQJHd7lOo3JBC10FlHe+38zJtAyeqLgKTh/V2/x9KbIwXDUBMlLhi1shI811IpUDLXQVMPYfP8UHS7fRv3EFYta9D7HN7RJrynfK1II+b8LO5TD/eddpVA600FXAGL1wM2kZhsdqbodDW6DV3a4jhYaGAyD+VjuH+qZZrtOov6GFrgLC3qMn+XjZDgY0jaXcuvegeCzU65vzFyrv6PFvKN8Qpt6h59P9mBa6Cghvz0/EYHio0Sm7LmbLIRCu67MUmMhoe9NRRhpMuUOvT/dTWujK7+08lMLkhJ1c3bwy5dZPgMgidu4RVbBK17TXp29fbE+/KL+jha783oh5vyEi3NeiBKz9HJpcC4VLuY4VmhpfA/UvhwUv2uvUlV/RQld+beuBE0z5aTfXt6xK2Y0f2V/5W93lOlboErFT7RYtB1/eDqeTXSdSWWihK7/25txNRIWHcVe7SrDyPajT0/7qr9wpXAouf8deaTTrCddpVBZa6Mpvbfr9OF/9soeb2lQjZut0SDmgo3N/Ub09tHsAVn8I67PP1adc0UJXfuuNuZsoEhXBHe2rw7LRUK4BVO/gOpb6Q6cnoUIT+Po+OKaLlPkDLXTll9btOcqMtfu4tW01Su1fCvvX2RuJRFxHU3+IiIIr34W0UzD9XtAJVp3TQld+aficTRSPjmBw+xqwdBQUiYEGV7qOpbIrUxu6PQuJc+3pF+WUFrryOz/vPMLcDfsZ0qEGJU5sh99mQfPb7M0tyv80vx2qtYdZT8GRHa7ThDQtdOV3Xp+ziVIXRHJz2+qwfDSER0H8YNex1LmEhUG/kYCBr+6BzEzXiUKWFrryKyu3HWLRpiTu6lSTohnH4OdPoOFAKKqLivu1UlWh+wuwdREkvOc6TcjSQld+5T+zNxJTrBA3tKpmz8mmpUBrnVUxIDS7GWp2hjlP22vUVYHTQld+48fEAyzbcoi7O9WkcHgmrBgL1TtCufquo6ncEIG+b0FYBEzTUy8uaKErv2CM4bXZG6lQIppBLarA+q/g2G6d8zzQlIiFni/Bjh9h+RjXaUKOFrryCws3JbF6xxGGdq5FdGS4vZGodC2o3d11NHW+mlwLtXvYFY4Ob3OdJqRooSvnjDG8PnsTsaUKc1WzyrBzBexOgJZ32isoVGARgcteBwmDrx/QG44KkP5rUc7NXv87a3cf5b4utYmKCINloyC6BDQe5DqayqsSsdD1GdiyAH6Z5DpNyNBCV05lZhqGz9lE9TJFuOLiSnZ5s/XToelNUKio63gqP+IHQ+WWdkbG5CTXaUKCFrpyasZ/9/LrvuM80LU2EeFhsGKcfaLFELfBVP6FhUGfEZB6AmY+5jpNSNBCV85keEbntcsW5bJGFe0//FUToF4fKFnZdTzlDWUvgvaPwH+/hI0zXacJelroypmvft7N5qQTPNStDuFhYs+1njqilyoGm3YPQkw9+PYhOH3cdZqgpoWunEjLyOTNeb8RV6E4PeqXtzehLB8DFZtC5Rau4ylvioiyNxwd2wNzn3WdJqhpoSsnvly1i+0HU3ioWx3CwgQ2z4cDm3TO82BVubm9DHXlu7Bzpes0QUsLXRW40+kZvDU/kcaVS9KlXlm7cdkoKFoe4vq5Dad8p/NTUKwCfPMgZKS7ThOUclXoItJTRDaKSKKIPH6OfQaKyHoRWScin3g3pgomk1fuZPeRkzzcrQ4iAvt/hc3zoMVt9tdzFZwKFYNeL8Hva+08Pcrrcix0EQkHRgK9gDhgkIjEZdunNvAE0NYYUx94wPtRVTA4lWZH582rlaJ97TJ24/IxEBENzW5xG075Xr2+UKsrLHhR1yH1gdyM0FsAicaYLcaYVGASkP334tuBkcaYwwDGmP3ejamCxUfLtrP/+Gke7l7Xjs5TDtmrWxoNhCJlXMdTviYCl74Kmekw60nXaYJObgq9ErAzy+e7PNuyqgPUEZElIrJMRHqe7RuJyBARSRCRhKQkvXMs1Jw4nc6Y7zfTtlZpWtUobTeumgDpJ6HlXU6zqQJ0YQ1o/zCsmwqJ81ynCSreelM0AqgNdAIGAeNEpGT2nYwxY40x8caY+JgYXYEm1HywdBsHklN5qFtduyEjzd4ZWqMTlIv7269VQabt/XY2zRmPQNop12mCRm4KfTeQ9ba9WM+2rHYB040xacaYrcAmbMErBcDxU2mMXbSFS+rG0KxqKbtxw3Q4vkdvJApFEYWg93/sykaLh7tOEzRyU+grgdoiUl1EooBrgOnZ9pmGHZ0jImWwp2B0DSr1P+MXb+NIStqZ0TnYOc8vrAm1urkLptyp0QkaDIDFr8PBza7TBIUcC90Ykw4MBWYBG4DJxph1IvKciPT17DYLOCgi64EFwKPGmIO+Cq0Cy5GUVN79YQvd48rRMLaE3bhzJexaCa3u0jnPQ1mPf9krnL59WOdN94KI3OxkjJkBzMi27eksjw3wkOdDqT8Z98MWklPTeah7nTMbl4+GQjrnecgrVg4uecrOxvjrN3ZiNpVnOjRSPnUw+TTvL9lG74YVuKh8cbvx6G5YNw2a3qBznitofhuUjbOXMaaddJ0moGmhK596Z9EWTqVl8EDXLKPzleMAo3OeKys8Anq9DEd2wI9vuU4T0LTQlc/sP3aKD37cRv+LK1GrrGcknpoCCe/DRZdBqapuAyr/Ub0DxPWHH16HIztz3F2dnRa68plRCzeTnmm4v0uWK1jX6Jzn6hy6v2D/O/sfbnMEMC105RN7jpzkk+U7uKpZLFVLF7EbMzPtpYoVmkCVVk7zKT9UsjK0fwjWT4Oti1ynCUha6Mon3pqfCMC9WUfnW3TOc5WDNvdCySrw3WM6xW4eaKErr9txMIXPE3ZyTYvKVCpZ+MwTy0ZD0XJQ/3J34ZR/iyxsr03fvx4S3nOdJuBooSuvGzH/N8LDhHsuqXVmY9ImSJwLzW/XOc/V37voMnsX6YIX4cQB12kCiha68qrNSclMWb2L61tVpVzx6DNPLB8D4YUgXuc8VzkQgV6vQOoJmPec6zQBRQtdedWbc3+jUEQ4d3WqeWZjyiH45VOd81zlXkxdaHEHrP4Q9q11nSZgaKErr9m47zhfr9nDzW2rUaZooTNPrJoAaSl23halcqvjMChcyt5BqvO85IoWuvKa4XM2USQqgiHta5zZmJ5q14+scQmUq+8unAo8hUtCpyfsJYybZrpOExC00JVX/Hf3UWau28fgdtUpVSTLm57rpsDxvdB6qLtwKnDF3wJl6tibjdJTXafxe7mabdGfzF63j6k/ZV9fQ7m2cd9xShSOZHD76mc2GgM/vg0xF0GtLu7CqcAVHmnvIP1kICSMh1Z3uk7k1wKu0I+cTGNzUrLrGCqbiHDhqUvrUTw68szGrYvg97XQ9y29kUjlXe3u9jLGhf+2b6xfcKHrRH4r4Ap9YHxlBsZXznlH5d7St6FIDDQc6DqJCmQi9majMe1g0avQ89+uE/ktPYeufCNpI/w2295IFBmd8/5K/Z1y9eHiG+wb7AcSXafxW1royjeWjrRLizUf7DqJChad/2F/puY8nfO+IUoLXXlfchL8MgkaX6M3EinvKVrWzsa48VudjfEctNCV9618FzJOQ6t7XCdRwabVPVCiir3ZKDPDdRq/o4WuvCvtpC302j0gpk7O+yt1PiKjoes/7XQAP3/iOo3f0UJX3rXmM0g5AG30RiLlIw2uhErxsOBfdklD9T9a6Mp7MjNh6Sgo3wiqtXedRgUrEej+PBzfA8tHu07jV7TQlfckzoUDG+1t/nojkfKlqm2gTi9Y/AacOOg6jd/QQlfes/QtKFZRVyRSBaPrM5CaDD+85jqJ39BCV96xd429lKzlEF2RSBWMshfBxdfDinFweJvrNH5BC13lX0aaHSVFFoFmN7tOo0JJpychLALmv+A6iV/QQld5l5kBaybD281h/VfQ8g67IIFSBaV4BWh9N6z9HPb85DqNc1roKm82L4B3OsKU2yGqKFw7GbroLdnKgbb3Q+ELYc4/Q35lIy10dX72/wofXwUT+8Ppo3Dle3DHIqjTQ69sUW5El4COj8HW72HzPNdpnNJCV7mTvB++fgBGt4Ydy6Hb83DPSmg4AML0x0g5Fn8rlKoGc54J6SkBcvUvUUR6ishGEUkUkcf/Zr8rRcSISLz3IiqnMtJh2Wh4qxn8NNFOh3vfT9D2Pp0WV/mPiCjo/H92QZU1k12ncSbHQheRcGAk0AuIAwaJSNxZ9isG3A8s93ZI5ciOZTC2I8x8HGKbw93L4NJXoEhp18mU+qv6V0CFJrDgRUg75TqNE7kZobcAEo0xW4wxqcAkoN9Z9nseeBkIzT/JYJKcBNPuhvE94OQRGDgRrv8SytR2nUypcwsLs1MCHN1pF8IIQbkp9ErAziyf7/Js+x8RaQpUNsZ8+3ffSESGiEiCiCQkJSWdd1jlY8bAqgnwdjM7yVbbB2DoCojrq294qsBQvQPU7AKLX4dTx1ynKXD5fjdLRMKA14GHc9rXGDPWGBNvjImPiYnJ70srbzq0BT7oA1/fbyfXuutH6PYsRBVxnUyp89PlaTh52K5pG2JyU+i7gayrMsd6tv2hGNAAWCgi24BWwHR9YzRAZGbY5eJGtYG9v0CfN+GmryGmrutkSuVNxSYQ19/+XJ844DpNgcpNoa8EaotIdRGJAq4Bpv/xpDHmqDGmjDGmmjGmGrAM6GuMSfBJYuU9+3+158lnPWl/Vb17mb11X0+vqEB3yVOQlgI//Md1kgKVY6EbY9KBocAsYAMw2RizTkSeE5G+vg6ofCAzAxYPh3faw8HNcMW7cO1nUKJSzl+rVCCIqQNNrrWrZx3ZmfP+QUKMo1tl4+PjTUKCDuIL3OHtMPVO2PEj1OsDvYdDUX0/QwWhIzvhrabQ6GroFzzn00VklTHmrKe09Ra/UGGMXYNxdFu7HmP/MfZyRC1zFaxKVob4wfDzx3DgN9dpCoQWeihIOQSTb4Rpd0H5hnDXEmgySM+Vq+DX/mGIKGxvNgoBWujBbsv3MKo1bPwOuj4LN38Dpaq6TqVUwSgaA63vgXVTYc/PrtP4nBZ6sMrMsKuif9jPzkZ3+3xo9wCEhbtOplTBajPUztM//3nXSXxOCz0YHdsLH/SF71+GJtfBkAVQoZHrVEq5EV0C2j1oFzHftsR1Gp/SQg82v82FMW3t6i2XvwP9R+rdnkq1GALFKsC854J6EQwt9GCRkWZXbPn4SihaHoYshMbXuE6llH+ILAwdHoWdy+C32a7T+IwWejA4/rs9xbLkDWh2C9w+z95YoZQ6o+mNdhGMec9DZqbrND6hhR7odiyHdzrA3p/tHZ993rCjEaXUn4VH2ikBfl8L66e6TuMTWuiByhhYMQ4m9LYFfttcaHSV61RK+bcGAyCmHix8KSiXqtNCD0RpJ+0CFDMegZqd7fnycvVdp1LK/4WFwSVPwIFNsPZz12m8Tgs90BzeDu91h18+gU5PwKBJULik61RKBY6L+kC5hnaUnpHmOo1XaaEHkq2L7Bqfh7fDoM+g0+N2xKGUyr2wMLjkSTi8FX751HUar9I2CBQJ42Hi5VCkrL1RqG5P14mUClx1e0HFpvD9q5Ce6jqN12ih+7uMdJgxDL55EGpcArfNgdI1XadSKrCJ2Cteju6Anya6TuM1Wuj+7ORh+HgArHgHWg+1i1BEl3CdSqngUKsLVG4Ji16DtFOu03iFFrq/OpAI73aFbYuh71vQ40WdWEspb/pjlH58D6ya4DqNV2ih+6MtC+HdznaEfuNX9g43pZT3Ve8AVdvZtUdTU1ynyTctdH+z+kOYeAUUq2invK3W1nUipYKXCHR+Ck7st+uPBjgtdH9hDMx/AabfCzU6weDZdt4JpZRvVW1jLzhY8gacTnadJl+00P1BeipMvQMWvWpPr1z7GUQXd51KqdDR+R+QctBegBDAtNBdO3kEProC1nxmf6j6jLCTCCmlCk5sPNTuAUtGwKmjrtPkmRa6S0d2wPgesGMZXD7WztesCzcr5cYlT8KpI7BstOskeaaF7sqen+1licf2wg1ToPHVrhMpFdoqNoGLLoOlIyHlkOs0eaKF7sJvc+H9SyE8CgbPspdOKaXc6/QEnD4GS992nSRPtNAL2i+fwadX29v3b5sLZeu5TqSU+kP5BlD/clg2Bk4cdJ3mvGmhF6SlI2HqEKjSGm7+FoqVd51IKZVdpycgLQWWvuU6yXnTQi8IxsDcZ2DWk1CvL1z3hV6WqJS/iqkLDa6A5WMDbpSuhe5rGekwfSgsHm4XcL5qAkRGu06llPo7HYZ5RumBdS5dC92X0k7C5Bvgp4+g42Nw2XCdYEupQFD2InsufUVgjdJzVegi0lNENopIoog8fpbnHxKR9SKyRkTmiUhV70cNMCeP2DlZNn4HvV6117jqNeZKBY6OwyD1RECN0nMsdBEJB0YCvYA4YJCIxGXb7Scg3hjTCPgCeMXbQQPK8X0woTfsWglXvgsth7hOpJQ6X2XrnRmlB8h16bkZobcAEo0xW4wxqcAkoF/WHYwxC4wxf8w9uQyI9W7MAHJ4m73789BWuG4yNBzgOpFSKq8CbJSem0KvBOzM8vkuz7ZzGQx8d7YnRGSIiCSISEJSUlLuUwaKpE0wvpc93XLTdKjZ2XUipVR+lK0H9fvD8ncCYpTu1TdFReR6IB549WzPG2PGGmPijTHxMTEx3nxp9/athfd7QWaavcY8Nt51IqWUN3T4Y5Q+0nWSHOWm0HcDlbN8HuvZ9ici0hV4CuhrjDntnXgBYudKe848ohDcMtPebaaUCg7l4gJmlJ6bQl8J1BaR6iISBVwDTM+6g4hcDLyDLfP93o/px7b+AB/2g8IXwi3fQZlarhMppbytwzBITfb7UXqOhW6MSQeGArOADcBkY8w6EXlORPp6dnsVKAp8LiI/i8j0c3y74LJpNnw8AEpWhltnQim9WlOpoFQuDuL6+f0oPSI3OxljZgAzsm17Osvjrl7O5f/WTYMvb7N/0ddPhSKlXSdSSvlSx8dg/TRYNsouRuOH9E7RvPj5E/jiFqjUDG76WstcqVBQLg7i+tuZGP10lK6Ffr5WjINpd9k5zG+YAtElXCdSShWUjsMg9bgdpfshLfTzsXQkzHgE6l4Kgz6DqCKuEymlClK5+n59Ll0LPbeWvHlm+tuBH+qMiUqFqo6P2VWN/HDtUS303Fj0Gsx5GupfAQPGQ3ik60RKKVf+N0ofAycPu07zJ1roOVn4Msx/HhoOhCvGaZkrpfx2lK6Ffi7GwPwXYeG/oPG1cPkYCM/VVZ5KqWBXrr49/bpstF+N0rXQz8YYmPccLHoFLr4B+o3UhSmUUn/mh6N0LfTsjIE5/weLX4dmN0OfERCmf0xKqWzKN4B6fTyj9COu0wBa6H9mjL2S5ce3oPlt0Hu4lrlS6tw6DLOj9OXvuE4CaKGfYQx8N8zeMNDyLrj0NS1zpdTfq9AI6vaGZSPh1DHXabTQAcjMhG8fsktNtR4KPf+t638qpXKn46Nw6iisHOc6iRY6mZnwzQOQMB7aPgDdX9AyV0rlXsWLoXZ3+PFtOJ3sNEpoF3pmJnz7IKz+ANo/DF2f0TJXSp2/DsPg5CFIeM9pjNAtdGPsvCyrJkC7B6Hz/2mZK6XypnJzu4bwkhGQmuIsRmgWujEw41H7f9M290GXf2qZK6Xyp+NjkHIAVr3vLELoFboxMPMJ+wZG66HQ7Tktc6VU/lVpZafVXvImpJ10EiG0Ct0YmP0PWD7aXpqob4AqpbypwzBI/h1WT3Ty8qFT6MbYGROXvg0t7tBLE5VS3letHVRpA4uHQ/rpAn/50Ch0Y2Des/DjCHsHaK+XtcyVUt4nYlc1Or4HfvqowF8++AvdGJj/gv0/ZrNboNerWuZKKd+p0QliW3hG6akF+tLBX+gLX4IfXoOmN0Lv1/V2fqWUb4nYK16O7oQ1kwr0pYO73b5/Bb5/CZpcD5e9qWWulCoYtbpAxaZ2tbOMtAJ72eBtuEWvwYIXofEg6KtT4CqlCtAf59KPbIe1nxfYywZnyy0ebpeNa3S1Lk6hlHKjTk8o39AzSk8vkJcMvkJfMgLmPgMNBkD/0VrmSik3/jiXfmgzrJtSIC8ZXIW+dKRdbaj+FXD5O1rmSim36vaGsvVh0auQmeHzlwueQl82xq42FNcPrhinCzorpdwLC7PzpR/YBOu/8v3L+fwVCsKKcTDzMbjoMrjyPS1zpZT/qNcPytT1jNIzffpSgV/oK9+z0+DW7Q0D3ofwSNeJlFLqjLAw6PAo7F8Pv37j25fy6Xf3tYT37dJxdXrCVRMgIsp1IqWU+qsGV0DpWvbeGGN89jK5KnQR6SkiG0UkUUQeP8vzhUTkM8/zy0WkmteTZrf6Q7t0XO3uMPBDLXOllP8KC4f2j8Dva2Hjd757mZx2EJFwYCTQC4gDBolIXLbdBgOHjTG1gOHAy94O+ic/fQzT74NaXWHgRIgo5NOXU0qpfGt4FZSqBot8N0rPzQi9BZBojNlijEkFJgH9su3TD/jA8/gLoIuIj2bAWjMZvrrHToBz9ccQGe2Tl1FKKa8Kj7Cj9D0/QeJcn7xEbgq9ErAzy+e7PNvOuo8xJh04CpTO/o1EZIiIJIhIQlJSUt4Sl6gMF/WGQZ9qmSulAkvja+xpYh9dvFGg1/cZY8YCYwHi4+Pz9jtH1db2QymlAk14JFznu7ldcjNC3w1UzvJ5rGfbWfcRkQigBHDQGwGVUkrlTm4KfSVQW0Sqi0gUcA0wPds+04GbPI8HAPON8eG1OUoppf4ix1Muxph0ERkKzALCgfHGmHUi8hyQYIyZDrwHTBSRROAQtvSVUkoVoFydQzfGzABmZNv2dJbHp4CrvBtNKaXU+QjsO0WVUkr9jxa6UkoFCS10pZQKElroSikVJMTV1YUikgRsz+OXlwEOeDFOINBjDg16zKEhP8dc1RgTc7YnnBV6fohIgjEm3nWOgqTHHBr0mEODr45ZT7kopVSQ0EJXSqkgEaiFPtZ1AAf0mEODHnNo8MkxB+Q5dKWUUn8VqCN0pZRS2WihK6VUkPC7Qs/PgtQi8oRn+0YR6VGgwfMhr8csIt1EZJWIrPX8t3OBh8+j/C48LiJVRCRZRB4psND5kM+f60YislRE1nn+rgNiqa58/FxHisgHnmPdICJPFHj4PMrFMXcQkdUiki4iA7I9d5OI/Ob5uCn71+aKMcZvPrDT824GagBRwC9AXLZ97gbGeB5fA3zmeRzn2b8QUN3zfcJdH5OPj/lioKLncQNgt+vj8fUxZ3n+C+Bz4BHXx+Pjv+MIYA3Q2PN56RD4ub4WmOR5fAGwDajm+pi8dMzVgEbAh8CALNsvBLZ4/lvK87jU+WbwtxF6fhak7of9IThtjNkKJHq+n7/L8zEbY34yxuzxbF8HFBaRQgWSOn/ytfC4iPQHtmKPORDk53i7A2uMMb8AGGMOGmMyCih3fuTnmA1QxLP6WWEgFThWMLHzJcdjNsZsM8asATKzfW0PYI4x5pAx5jAwB+h5vgH8rdDzsyB1br7WH3lrEe4rgdXGmNM+yulNeT5mESkKPAY8WwA5vSU/f8d1ACMiszy/qg8rgLzekJ9j/gI4AewFdgCvGWMO+TqwF+Sng7zSXwW6SLTyDRGpD7yMHc0Fu2eA4caYZM+APdhFAO2A5kAKME9EVhlj5rmN5VMtgAygIvb0ww8iMtcYs8VtLP/nbyP0/CxInZuv9Uf5WoRbRGKBqcCNxpjNPk/rHfk55pbAKyKyDXgAeNKzRKI/y8/x7gIWGWMOGGNSsCuHNfV54vzLzzFfC8w0xqQZY/YDS4BAmOslPx3knf5y/UZCtjcMIrBvBlTnzJsK9bPtcw9/fiNlsudxff78pugWAuPNo/wcc0nP/le4Po6COuZs+zxDYLwpmp+/41LAauybgxHAXKC362Py8TE/BrzveVwEWA80cn1M3jjmLPtO4K9vim71/H2X8jy+8LwzuP5DOMuBXgpswr5b/JRn23NAX8/jaOzVDYnACqBGlq99yvN1G4Fero/F18cM/AN7rvHnLB9lXR+Pr/+es3yPgCj0/B4vcD32DeD/Aq+4PhZfHzNQ1LN9nafMH3V9LF485ubY37pOYH8bWZfla2/1/FkkArfk5fX11n+llAoS/nYOXSmlVB5poSulVJDQQldKqSChha6UUkFCC10ppYKEFrpSSgUJLXSllAoS/w/cXHW71ezsNQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "x = np.linspace(0,0.1,2000)\n",
    "ratio = 2\n",
    "first = 1\n",
    "(length,)  = x.shape\n",
    "slop = int(length/6)\n",
    "second = ratio-first\n",
    "odd = ratio % 2\n",
    "\n",
    "first = int(first * length/ratio) \n",
    "second = int( second * length/ratio) + odd\n",
    "slop = np.array(np.append(np.zeros(first-slop) , (np.arange(slop))/slop))\n",
    "#steep = np.ones(int(first * length/ratio)+ odd) - np.exp(-np.arange(int(first * length/ratio) + odd)/200)\n",
    "steep = (np.ones(first) + slop)*0.5\n",
    "\n",
    "step = np.append(steep,  np.ones(second))\n",
    "m = np.sin(5 * 2.0 * np.pi * x) * step \n",
    "plt.plot(x, step, '-')\n",
    "plt.plot(x, m, '-')\n",
    "plt.savefig('m_t.pgf', format='pgf')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'np' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[1;32m/home/joshua/Documents/SeminarSpezielleFunktionen/buch/papers/fm/Python animation/Bessel-FM.ipynb Cell 6\u001b[0m in \u001b[0;36m<cell line: 1>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> <a href='vscode-notebook-cell:/home/joshua/Documents/SeminarSpezielleFunktionen/buch/papers/fm/Python%20animation/Bessel-FM.ipynb#W6sZmlsZQ%3D%3D?line=0'>1</a>\u001b[0m x \u001b[39m=\u001b[39m np\u001b[39m.\u001b[39mlinspace(\u001b[39m0\u001b[39m,\u001b[39m0.1\u001b[39m,\u001b[39m2000\u001b[39m)\n\u001b[1;32m      <a href='vscode-notebook-cell:/home/joshua/Documents/SeminarSpezielleFunktionen/buch/papers/fm/Python%20animation/Bessel-FM.ipynb#W6sZmlsZQ%3D%3D?line=1'>2</a>\u001b[0m y \u001b[39m=\u001b[39m np\u001b[39m.\u001b[39msin(\u001b[39m100\u001b[39m \u001b[39m*\u001b[39m \u001b[39m2.0\u001b[39m\u001b[39m*\u001b[39mnp\u001b[39m.\u001b[39mpi\u001b[39m*\u001b[39mx\u001b[39m+\u001b[39m\u001b[39m1.5\u001b[39m\u001b[39m*\u001b[39mnp\u001b[39m.\u001b[39msin(\u001b[39m30\u001b[39m \u001b[39m*\u001b[39m \u001b[39m2.0\u001b[39m\u001b[39m*\u001b[39mnp\u001b[39m.\u001b[39mpi\u001b[39m*\u001b[39mx))\n\u001b[1;32m      <a href='vscode-notebook-cell:/home/joshua/Documents/SeminarSpezielleFunktionen/buch/papers/fm/Python%20animation/Bessel-FM.ipynb#W6sZmlsZQ%3D%3D?line=2'>3</a>\u001b[0m plt\u001b[39m.\u001b[39mplot(x, y, \u001b[39m'\u001b[39m\u001b[39m-\u001b[39m\u001b[39m'\u001b[39m)\n",
      "\u001b[0;31mNameError\u001b[0m: name 'np' is not defined"
     ]
    }
   ],
   "source": [
    "\n",
    "y = np.sin(100 * 2.0*np.pi*x+1.5*np.sin(30 * 2.0*np.pi*x))\n",
    "plt.plot(x, y, '-')\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "916dbcbb3f70747c44a77c7bcd40155683ae19c65e1c03b4aa3499c5328201f1"
  },
  "kernelspec": {
   "display_name": "Python 3.8.10 64-bit",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.10"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}