Skip to content

Instantly share code, notes, and snippets.

@ahartikainen
Last active January 28, 2019 10:00
Show Gist options
  • Save ahartikainen/089ce30d95fd13c53cbb137a9038f9ed to your computer and use it in GitHub Desktop.
Save ahartikainen/089ce30d95fd13c53cbb137a9038f9ed to your computer and use it in GitHub Desktop.
Adhoc edgecorrections for Scipy gaussian_kde
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import numpy as np\n",
"from scipy.stats import gaussian_kde\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"plt.style.use(\"ggplot\")\n",
"plt.rcParams[\"figure.dpi\"] = 150"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"def kde(arr, x, lower_edge=None, upper_edge=None, bw_method=None):\n",
" arr = np.asarray(arr)\n",
" arrs = [arr]\n",
" \n",
" # check edges\n",
" if lower_edge is not None and upper_edge is not None:\n",
" assert lower_edge < upper_edge\n",
" if lower_edge is not None:\n",
" assert lower_edge <= arr.min()\n",
" if upper_edge is not None:\n",
" assert upper_edge >= arr.max()\n",
"\n",
" if lower_edge is not None:\n",
" arrs.append(lower_edge-arr)\n",
" if upper_edge is not None:\n",
" arrs.append(upper_edge*2-arr)\n",
" arr_ = np.concatenate(arrs)\n",
"\n",
" if bw_method is None:\n",
" bw_method = len(arr)**(-1./(1+4)) #Scott\n",
" bw_method /= len(arrs)\n",
" \n",
" kde_func = gaussian_kde(arr_, bw_method=bw_method)\n",
" y = kde_func.pdf(x) \n",
" # normalize by integration from lower to upper\n",
" # adhoc\n",
" y *= len(arrs)\n",
" return y"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"rng = np.random.RandomState(seed=19)\n",
"data = rng.rand(10000)\n",
"x = np.linspace(0,1,100)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"y = kde(data, x)\n",
"y_le = kde(data, x, 0)\n",
"y_ue = kde(data, x, None, 1)\n",
"y_lue = kde(data, x, 0, 1)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"''"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAwoAAAIJCAYAAAAF28o7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3Xd8U/X+P/BXRtM0Helu6aQtHZQCLaMtW1ARFypcQXHc6w83bvR6vVcu1+t1fMWLXr1X9LqugiiIiqKIsqGsUsrsooPuvdOVNOP3R2namlJG056M1/Px4ME5OSfJu00zXvkskcFgMICIiIiIiKgXsdAFEBERERGR5WFQICIiIiIiEwwKRERERERkgkGBiIiIiIhMMCgQEREREZEJBgUiIiIiIjLBoEBERERERCYYFIiIiIiIyASDAhERERERmWBQICIiIiIiEwwKRERERERkgkGBiIiIiIhMMCgQEREREZEJBgUiIiIiIjLBoEBERERERCYYFIiIiIiIyASDAhERERERmWBQICIiIiIiE1KhCxBaZWUlDAbDsN+vt7c3AKC2tnbY75uEw8fdfvGxt1987O0XH3v7JdRjLxKJ4O/vb7bbs/ugYDAYBAkKve+f7A8fd/vFx95+8bG3X3zs7Ze1P/bsekRERERERCYYFIiIiIiIyASDAhERERERmWBQICIiIiIiEwwKRERERERkgkGBiIiIiIhMMCgQEREREZEJBgUiIiIiIjLBoEBERERERCYYFIiIiIiIyASDAhERERERmWBQICIiIiIiEwwKRERERERkgkGBiIiIiIhMMCgQEREREZEJBgUiIiIiIjLBoEBERERERCYYFIiIiIiIyASDAhERERERmZAKXQCRvUkpqENTk0qw+08MchXsvomIiMh6sEWBiIiIiIhMMCgQEREREZEJBgUiIiIiIjLBoEBERERERCY4mJmIiGiIcRIDIrJGbFEgIiIiIiITDApERERERGSCQYGIiIiIiEwwKBARERERkQkGBSIiIiIiMsGgQEREREREJhgUiIiIiIjIBIMCERERERGZYFAgIiIiIiITDApERERERGSCQYGIiIiIiEwwKBARERERkQkGBSIiIiIiMsGgQEREREREJqRCF0BERERkq1IK6tDUpBLs/hODXAW7b7J+bFEgIiIiIiITDApERERERGSCQYGIiIiIiExwjAIREZGVEWk10Dc3ob2uCtDrIFN6Qqb0gN7RCRBLhC6PiGwEgwKRBdE3NUBVmIOWynKoqsuhqq1BU30tmhsb0NLWCrnMEQonBZycneGkcIaTiyucXFzhNXIUfCbNAFyVQv8IRDRI+k4NVLkZqM48gcq8bJQVF6KsuhrNbe1QqTVQdXZC1amFWm/ocz0RAGepGM5SCZwdHOAikyHY2wshwSGojR+H6ImTETpxMqTOLsL8YERkdRgUiIRgMKCtpAAVaSkoOnMCBQV5yC2vQEFzC7SGi1+9P04SMWI9lYgJDUXUmLEInTwdyvHJgMzRvLUTkdnoOlpRdWQfClJTkJN1BvnlFShsbkG7Tn/Zt2UA0KLVo0WrBzo6AbThVF0jkJMP7NgNAHAUixCudMHYkCBMnzYN029dCN/YsRCJROb9wYjIJjAoEA0DnbYTZ/fvwfo9u7Ar5SCOF5Wgsl1t1vto1+lxrKYBx2oagLQTwGdr4SmTYnZ0JGbetBDBN9wOEUMDkaBaivKQv+sn5J48hqyCfGTU1Ju0DAwltd6ArAYVshqysPFkFvDeR4h2d8H0MTGYMXsOkm+9Ha4jAoatHiKybCKDwTB8r1AWqKKiAkL8Cnx8fAAANTU1w37fNPS0nRpk7NqOQ7/+jCPHjiO1qBSNGu0lX18pc4Cnkxzuzs7wcHOF0k0JZ1c3aNQdaG9rQ2t7G9o7OtDa0QFVeweKm5qhu8ifcaDCEdclJOCeB+9H5JzrIBJzLoPhxOe8fSrPysDJHb9g2/adOJabi3PNrZd8XV8nR4R6eSI4IAABIWFQ+vrDSekBJw9vOHn5wNHTGwp3L4jEErSrmqBpaoC6qR7q5kaom5ugqqlEef5ZFJcUoai6GkVNLbjYu52DSIQZoQG46bq5uO7/PQj3oJDB/QLsnI+Pz/l1FJoEq4HrKAhDqNd8kUiEESNGmO/2GBQYFGjwdDodsvbtxoGtP+LQ0aM4UliK5s6LBwMPRwdE+/kiYmQYQmLiEDhxKnxjx8PhMr/5V7e3ozT9EApT9yMv8zSyi4pQOMACP3Geblhy/VwsemElnDw8L+u+6MrwOW8fyjJP4+APm3H40EEcyslDkartotdxEIkQ7eOJ2PAIRIyNh39sPHziJsDJzd1sdSUGuaKttRUF6WnISj2EQykp2J+RjfLW9gteRyYWYWZYMG6edx3mLn0Ibn7m+/BhLxgU7BeDgo1gUBBGaqlwq1QCg3/hbG9rxZldO3B0zy6kHT+OIwVFF20xEAMY4+uJ8dGjETkhCSHTr4EyOGzI+ga3NNYj47v12LvtBxwpLEZ/PZ595TI8fONc3L3i73D28RuSOqiLvT/nbZHBYEBpxikc+vF7HDrYFQxKWi78wbubv0KO8aEhiI4di7DE6QhInAmZwnlIa+3vNc9gMKDg1Ans/24TDhw4gAN559B0gdcxR7EIcyJH4vbbb8fsPzwAmZNiSOu1FQwK9otBwUYwKAjDmoKCwWBASW4OTu/ZhbQDKTiWmYXTlTXQXKRfsVQkwng/LySPHYOkWVdh0k23YFRsnCBvGo1lxTj65YfYvWcnsmrqTY57yqS4f+4c/OGvf4cyMHhYaxtOQv7dKZVKTA/3stvnPCDs798cH5a0nZ3IObAfqdu3ITUtDUcLilDR1nHR6wW7KDAhIhwx8ZMQOft6eMaMG/bBw5fy82u1WhzevAlbNnyFn4+dQJ26s9/zvBwdcGvSRCy+/0HEzpnLgdADYFCwXwwKNoJBQRiWGBQMBgMa6uuRm3YEWalHkJVxBlkFhciprkFLp+6itykGMM7XE1PHjsGU2XOQOP82uHj59DnHEt40lJVZ+N8/38CGlMMmgyjdpBLcN3s6Hn7tTbjZ4IBGBgVhWVtQaK6pxolff8axAylIP3MGacVlaL6E14JQVwWSo0YhecoUTJ1/KxKuulrw5/3l/vydGg0Of7sRP2zcgJ+Pn0aDpv/QEOulxO1zr8Vty56AT1iEOUq1KZbwms+gIAwGBRvBoCAMIT4wtLe1orG+Ds01VfBsKUfp2bMoKTyHkrIylFbXoLSxGS2XMK6gm6NYhHh/b0yKiUbi9OlIvGUh3PwH/nBtCW8a3ZrKCvHL269g84EDJlMxeskc8P9uuQXJj78IsaPcrPcr5JsWg4KwLDko6LRa5B05iGO7duB4+jEcyy3A2Ybmiw7+BYCRrgpMiR6FKVOnIfnmWxEYG9fnuCU87wfzvNOoO7Bv/ef4+ssvsT0rt99ZmqQiEWaPCsXtt9+Oa+57AI5D3JXKWlj7Y09XjkHBRjAoCMNcHxja21pRW1mBxvpaNDXUo6mhDs211VBVlkFVU4XG+jo0NDWivqUF7ZcRAvoTpHBEfHAgJo0fh0mzrsKYa+ZB5nJ5L8CW8KbxWy1V5djxziv4Zs8etGj7flsaqXTB/Q88iqjf3QeYqXsBg4J9PucBywoKVXm5OL79Z6SnHsGJ7BycKq+GSntpLYdjvD2QODoKk6dMReINN8MvMnrA61jC895cz7uGinL88N472LT1Z6RX1vZ7jrtMilsTJ+D2pQ9i/LXz7Lprki099nR5GBRsBIOCMC7nA0NDXQ0KcrJQei4PNZXlqKmsQE1lOWory6FqNv+Lr5tUgmhPJWKCAzE6OhqjEyYgeuoMuIUOfuCxJbxpXEhrfQ12vf0yvtyxw6SFYXpIAO5+biX8EmcN+n4YFOzzOQ8I9/vXtrdBnn0Ix/btxrGTp5BeWIyy1ouPLQAANwcpEgL9MCF2NCZPn4EJ198MV9/LG/hvCc/7oXje5R5KwdcfvIdvUg5fcF2YCKULbp4+FfPvW4roKdPNXoOls9XHni6OQaGXgoICnDp1Cnl5ecjNzUVDQwMcHBzwxRdfXNHttba24uuvv0ZqaioaGxvh7u6OyZMnY9GiRXB2Nm9zJoOCMC70gaGmshx5madRcDYLBdmZOHc2Cw115vsdidA1cNfb0QGBbi4I8vZCcMAIBIeORHBkNIJjx8AjMgbiIZrRwxLeNC6mqSgf37z6J/x44nSfbhdiAAsmxuPWFavgPOLKBzwzKNjncx4Yvt9/e3U5zu3Yguy0w8jIz0VGdd0lLWomBhDp4YaJo8IxYeIETJxzLSKSpkIiHdzapNbwvB8MXWcn8n/cgD3fb8TunDx06PtfVTraU4n5s2Zg/tIHEZ4wcZirFIYlPPYMCsJgUOjljTfeQFpaWp/LrjQoqFQqvPjii6ioqICfnx/Cw8NRWlqKkpIS+Pv745VXXoGrq/n+6BkUhNH9gUGr7UT2yXSkHdiL9IP7UHIu77Jvy0Ekgo+jFF4yB+P/nnIZ3N09oPTyhpuPP1z9A+A6IgSKwFDoPH0xPiYEokG++V8JS3jTuFSlh3bj81Uv4WhZZZ/LlQ5S/L9bbsW0J1dAcgUrPTMo2OdzHhi6339bWTHyd25BVtohnMjNRXZ90yWNLfBzckRCSCAS4sYgYdp0jLt6Lly9fc1enzU97wdLXVeF41/8F9t3/IL0qv67JgFArJc7Zk+eiDk334KJN9wMB5lsGKscPsbHvrER2qY6NJ/LRVPJOXQ0N6Gzox2ajjZ0dnSgU92BTo0aOq0Ozq5ucHZ3h4uHNxSe3nD28oWTjx9kI0IAqcNl18CgIAwGhV42b94MjUaDiIgIRERE4MEHH7zioPDvf/8b+/btQ2JiIp5++mlIJBIAwCeffIJt27Zh5syZeOyxxwZbshGDwvBramrCmi8349jBfThx5ADaWi7+4cHHUYoxbgqEKhwRIJch0Knr3wilG+QRo6EZMRJqT3+oPf3Q4emHTldPYICVh4V64bS2DwwGvR5nNn6CTz5ag0JV3xVlI93d8OATz2LUjYsu6zYZFOzvOd/NXL9/bWsLCrdvxql9O3A0MxPZDc0XvY5MLMJYXy9MiInGxOQkTLhmHgJHjzFLPRdjbc97c2nMOoG0rz/D3iOHcKa28YLnuTlIMTMqAlfNmY3Zi++Cf1j4MFZpPqrqShQeT0dBxmkU5p5FcWkpqhoaUdbQhKrWNjRdwoxZA3F3kGC0tydGhYQgLCoWwQmJ8EpIhkgx8Gsqg4IwGBQGsGjRoisKCo2NjXjooYcgFouxZs0auLv3rErZ2dmJRx55BC0tLXj//ff7HBsMBoXhYTAYcPz4caxbtw7ff/89Ojou3D842FmOsa5yxLkpMMZNgTFKBXwdu75FafcOQGtQJFqCI9EaFIkOnwBAdOFAcCEMCpens6Md+1avxP9++slkwPM1URG4c8X/wTMq7gLX7otBwT6e8/254t+/wYC6k0dw5pfvkZ5+FEdLyk3G0fyWp6MDEkKCMXrMWIxKmombF94EubPLld3/IFnr895cEoNckX/oALZ89jF+2H8QOY0D/x0EuyqQEBaKCfHxSJh5FeKumgO5hSzw1qZqRmF6GgpOncC5szkoKCzCucpKFNY3oaZDM+z1yMUijPZyx5TxCZhw8+/glXSVyXsig4IwbCUoDH/fiwEcP34cBoMBsbGxJkHAwcEBEydOxO7du3HixAlcddVVwhRJl0WlUuHbb7/FunXrkJmZ2e85DjIZJkRGYp5Ujavd5Rjp3Lc7i1rpjfL4maiLnwmNu0+/t0FDy0HuhKv//AYS734Im17+I344lWE8tuNsPvb9/nYsmjEd1//xFSiGoOsG2R9dpwYlO7fg2C9bcOD0KZxTtQ14vrdchgkjQxE3fgIir7oOfvHJEPVqVRQqJFCXiCnT8NSUaXjSYEDe/j3Y+e3X2H34CFLLqkwWryxRtaHkVBZ+OJUFfP4lHEQixPp6ISY0CCNDRyIsMgqhsWMQOj4BSk8vs9apUatRkZuN0qxMlObnoayoEKXlFSiprkFhbT3KL2GBvYGIAHjLHeEmd4TcQQpHBwfIpFLIHBzgKJMBIhFa2trQ3NYGVXsHmtQaNGs6L9iVrkNvwPGaBhzfsQvYsQsRrgrMGDsWk+fdghFzboLIwTa7dNHwsaigUFRUBAAICwvr93hYWBh2796NwsLCYayKrkRZWRn+/e9/Y9OmTWhrM32Dd1W6Y+rV8zAlLg43lp+BX3kugJ75+vVSBzSMTkRdwiyoRo6+olYDMj/XkAjc9+E3uHr3T/j4n6/g1PlVnjV6A9bt3Y8tB+fg7ptuxqwnV8DBQr4BJOvR2dyAnB++Qtre7UjJyb3gysAAIJeIMTkkEBMmJmH03PnwHTupTzAgyyQSiRA5czYiZ87GwwBaqiuR8vWX2L19O/Zk5KC0nw/inQYDTlbV4mRVLZB6os8xT0cHBChdoVQ4QalQwNXFGW6urnBzU8LV1RUGgwE6nQ4GvR46nQ46nQ56vR4tLSo0NDahoVmFhpYWNLa1o6FdjQa1BgO3VQ3MUSzGSHcXhPn4IDQoEOEREdC7+8HJewTcRkbALSAUEofLG2eg0+nQUl+LyuOHUHwiDQVns5BXWoqChiaTWvNVbcg/eAT/O3gE/q/+DddNiEfAn15A0MTEQfxUZM8sKijU1nYNfPL09Oz3uJeXV5/zyPJUVlbi3//+N7744gtoNKbNsMnJybj77rvhOzoJoak/w+/ADxDre7qy6BwcUTHrNtROnAOdnAv2WKqQ2Tdi5czrcOyDN/HRhvWoPt/k3tSpxX+++w6bf92G39/9B0z8/WMQnx9nRNSfjvpqZHy7Fod278CBc0XoGKBLUaSHGxLj4jBu5lyEXXMzHLiol9Vz8fXHvGVPY96yp6HX6VB2Mh3pO7fj+LE0HM/Nw+nq+gFnrKpXd6K+un4YK+6aQCNU6YxwH2+EBQUhLCIcYaPHICw+Af6RMcaxlYB5up1JJBIoffygnHsroufearxc3d6G8vTDyNixBQdSDyO7tqHP9So7NPjsYCrWzr8N10WOxNL7lyLpznshlljURz+ycBb119Ldb93Rsf+ZVLovV6v7n6+5P88884zJZTKZDK+//joAwNvb+3LLNAvp+Rl3uvuwWbvq6mqsWrUKH3zwgcn4A3d3d9xzzz1YunQpYmNjock4gcrVL0FWV9HnvJa4ZNTe9iC0Hr4Y6k4CPj7mba6+VFKptOtFX6kU5P7N7ZoXXsP0h5fjp5efw2e/bjeOXyhrbcerH6xB7LcbcM/9j2LSXQ8Zp5gU6ncPAEqVcN84SyQSSKVSm3nOX4nu339LeQnSvvwQ+3b+ioOFxSZdT7pJRSJMDg7AtOmzMPn238Nn9Lgrvm8h/+5s7Xl/uS7nd+/nfwMmXHeDcV+tasbJX37C0T27kZeXh4LSUhTW1OFcU+sFp2E1B1+5DEFKVwT5eCE0IAAhoaGIGB2LmMlJGJkwEdJLnPFtSB97pRK+N9yG+Btuw10Aqs5mIHXD/7B/7y6kl5YbWxv0AH7OLcTPz69A7P+twiOLb8fdL74EZx92Ex1KtvI5z6KCQvegYntexdHaNDY24o033sB7771n0sUoMDAQf/rTn3DvvffCyckJANBxZB8aV70ImbanS0Gnuw9qFzyE1rjkYa2dzEPu4Y2Fqz/F3MJcbFj5NDampkN7/rmcWVOPF177B0L/8zYWL1yEq5a9AEC4D2wknNLsTPy66p/Yt3cnUksqjH8jv+UilWBG1ChMvfo6xN9+H5x9/c1y/ykFdWa5nSshYavaFXN0dUPi7+5E4u/u7HO5XqNBWeZp5B47isqSYjTW1aGxsQFNjU1obFahuaUFqvY2iCCCRCyGWCKGWCSGRCKGRCyBwkkOL3d3eHp5wcvHB16+fvAOCIBPUDBC4sZDYeaxD8PBL2oMbl6xCjevAJoqyrDvw9XYtOV7lPSasS6zvhmPr/kYKz/+DE8vvAWPrX4Xzl7CfGFK1sGigoLxw+QFZsTpbkm4UItDf1avXj3g8draWs56dAX0ej02btyIV199FXV1fd+AfX198fjjj2PJkiWQy+VoaWlBS0sL9If3wPDp28D5b4EMIjGqkueh4qqF0MvkwDDOCFJTM3TfRA3Ex8cHOp3ONmc/8fDF7e98gatPHMbGN1/CL7nnjIeKmlvwxqef4MMv1uLe66/BH/78V/gGhQx7iU1NQzDrkV4PkV4Hg0gEdP9D9/89lEoltFqt1T7nr0ThyePY+tnH2LYvBekVNRcckOkuk2LW6Bgkzb0JkTcugvT8+BYtYBPPle5vk23hZ7kSQ/V6Kw8MwdjAEIw18+226vRoNdPzVLDXfIULZj75V0xb9me07P0WH3/4IfafKzEertdoseLLb/CfzVvw5O0LcMeLL0HGAf9mxVmPhkB3N6D6+v77G3Z/IBWqu5AtGcw0kXlZp/HfVS8jN+NUn8vdPDyx4J77MW/hnXCUO+FUbSeArpYD76M7ELL1fxCd/6igdVQgb8mzaA2JvuI6yDJ5xifj4XU/4+adP+CHT97DjrxCYxN4g6YT//r+Z7y/ZRtuiI/DDbfehqvuuAsKC3yD0jY3ojE3A7V52agtLkBNRTlUqiaoWlvR3NaGpvYONHWo0ajRotOghxiirpyArn9ikQhSsQjujjJ4OMnh6eKCEV4ecHd1hZe3NwJDQhEWOwah4yfCxd/fJlpSNeoOHP3hO+z6cQt2px9HTv2F1zfwdpRhVlwsEq+bj/DrF0Aqk1/wXCK6MhKpFNf+fimu/f1S5KQexqf/fAObDh01Ti9c2a7BC59/hfe//R7L770btzz7AqSX8WUs2T6LCgqhoaEAgHPnzvV7vPvy7vNoeDU3NmDde6ux44dNfVphHOVOWPiHh3DT4nvg1M/gQr+ULQja+ZVxv1Phhtx7nke7/8jhKJsEEnj1fDxy9XwsOXkEv3zwFr4/cQpt59+c1HoDvks/je/ST8PppX9gTmwUrr/pZlx91+/h5uExbDVq1e2ozzyJ6swTqMjLQWlJIUpralDepEJl+6WPhQIAHQzo+5W5AdABzZ3tKG5pB2oagF7f6PXm7eiAUA83hPn5YmRICKLGxCF6UiJCJ0yGg4W/aZedzcaejeuxe89e7M89Z7LORm+BznLMGjcOk6+/FcFzbr7s2V+I6MpFJybj9Q3fYnlJMd7983NYu+eAcXxQUUs7nnjvQ/xn/Vd48cnHMfuBR23iywsaPIsKCvHx8RCJRMjKykJTU1OfwT+dnZ04duwYRCIREhISBKzS/uj1emzfvBHr1ryFlua+zafTrrkef3jij/D266eZy2BAwK6NGJHyg/EijZsnzt7zAtTeAUNdNlkI5fgkLHrvK9xSdg7716zC1/tTjLMkAUC7To+fTmfjp9PZkP3fm0gOC0FcdBRixo1DTNJUjBqfcFndDX9Lo1ajNCsDRadPIv34aVSWFqG8sgJFNbUoUbVdsK/8cKpVd6K2sg7HKuuAk1nAll8AAI5iEcLd3RAV4I+oURGIGhOHyAmTERo/ATL58H8D39nZicz9e5C2czvS0o/hWH4RylrbB7xOpLsrrk+ahBvuWIK2mGmcwtQOCbnQIcAFx37LJzgEf1+7AQ/lZOHtPz+PDUeOQXf+ZTCnUYV7XnoVc9etw99W/wuhkzitqr0TJChs27YN27ZtQ2JiIpYsWWK83MPDA9OmTUNKSgo++ugjPPXUU8ZBYOvWrUNzczNmzJhhtlWZ6eJKCwvw3qsrkHXyWJ/Lg8NG4f5nX8S4SRcYgGwwIPjnz+F79FfjRR2efsi95wUummanHAPDcM0/3sMzzlps//A/+PmXX7AzrxAt2p7+yxq9Afvyi7AvvwjYuh3APyEViRDhqURUUCDc3dzg5CSHk5PT+X8KODkroO7oQGN9fddgxqZmNKpUaGhpRbWqBeUtbcY3wcslBjDCxQmB7u7w8/aGu7sHXJUecHb3gLOnN5w9faDw9oXUyRkw6KHX6wGDHgadHgaDHp0d7WitrkRLTSXamxuB1kZUVlejtqERJXX1KGxUQT3AzC1qvQFZ9U3Iqm8CzuQAm7cC6JoNaKTSBaP8fRERGoqwUaMQGBaBgMhoBMSMhsJlcF25NBoNSrIzce7EcZzLykBBQQFyi0twoqzyoisiO4pFmBISgNlTp2D27xYhPHGq8ZtJoT8wElGPwOjRWPXND3j4eBr++ZcX8P3JnkVRf80vxt7bFmDZ3Nl45J/vQOE+fC29ZFlEBjOM5E1PT8c333xj3M/NzYVIJMKoUaOMly1cuBATJkwAAGzcuBGbNm3CrFmzsGzZsj631dzcjL/85S+oqqqCn58fIiIiUFJSgpKSEvj5+eGVV16Bm5vbYEs2qqiosMvBzBd7w+7s1GDz2o+x8ZP3oO3smaHISeGMxQ88hhsX3Q2p9MLdBnwPbUXwr18Y99t9g3D2nhegdbGMkCfUN0zmmFPb2vX+3bfX1yLli8+w9aef8Gt2Pho7tcNejxhAsKsCod7eCA4MQsDIcHiHRcIrMhbK8ChIHczT9UepVGJ6uFef57xOp0PVuXwUnjyOwqxMFObnIa+oGGcrqlDU3HrFCz95ODogwNUFfu5KODvJ4STv9c/JCXInJ2g0GuNEAy1tbWhpa0dLRweqm1UoaW69rJaWUBcnzI4bjdnXzcPU391xwRlj7DUo2PtgZqEJ2aJgCa/5l/rzZ+zbhRefew6ppZV9Lg92dsLfnlyGuY88CTFbBC8ZBzP30tzcjNzc3D6XGQyGPpc1N194UFtvbm5ueO2117Bx40YcPXoUqampUCqVmDdvHhYtWgSXQX5TRhd3NuMU3nv1RRTlne1z+dSr52HpM3+Gp/fAcy8ryvIRtKNnTEJrQDhy734eOifLeeyE+sAi5Dz+lsjJ0xvXPr4c1z6+HJ2qZhz5ej2OHz6E7PwCZFdUI6+5BVoz5HgRgBEKRwS5uyPQ1wf+I4LgHx4Jv9h4eMXGQzpMC3d1fWD4zd+efASkSSMwKukGjAJwzfmLNa3NqD2Zhoozx1Gcm4Xi0mIU1tSiWNV20QDRoO5Eg7oBGb9ZgMkcHEQijPWPTXctAAAgAElEQVTxwMSoSExKTsLEa6+H/5ix7M9MZOXGzJyDbw4dxXdvr8LL//kANR1d47RKWtux9NU3MeerDXh9zQcIjBsvcKU0nMzSomDN2KLQo6O9Des/+Bd+2rC2qwvFeZ4+vnjwuZVImnX1RW9X0tGK0R/8BY6NXT+XxtUDmQ+/Bp2CfUQBfrMIXN63e+q6GuQfPoCs9DScy89HW3s72jvUaFer0a7WoF3T9c9BKoWHizPcXV3hrlTC3cMd7h5e8PTzQ0hMLILHx0Ou9BD0G21zPfba1hbUZRxDZcZJlOfnoLSkCBX1DahUtaC6Q33F3ax+SwQgyFmOkZ7uCA/wR1joSIxLSsbYuddD4XllM8+xRcF+n/dCYovC5f/8zbW1eOuZx/DJrpQ+rYuuUgn++vsluGPlKxBzfZABsUWBbMqJwylY8/pKVFeU9bn8ugV34J5ly+HscgkvNAYDQrd8ZAwJBpEI5xY+xpBAV8zRywexN96K2BtvFboUiyF1doFf4iz4Jc7Cb7/X07e3obU4D43nclFXcg7NNdVQd7RB3dEBtboDIl3n+ZClhkzqABeFE1wUCri4OMPFxQXOLq7w8PFG+JhxCEmYCLmH9S06RUSD5+btjZWff4U7DqXgxeVP42BROQBApdXhuY/XYsv2HVi15r8Iip8gcKU01BgU7FxzUwM+fft17Nn6fZ/LA0JGYtmf/4HYhEmXfFveaTvhkZlq3C+/6ndoCY0xW61ENDCxkwKu0ePgGj0Owf0c5+wvZI8EbUm0gO6mg/r5g8fj2a92IPWj1fjP559C1dk1/fG+4gpcfcutWLZwIZKeXAnRBVoX+Jpj/YT/CyZBGAwG7P/1Jzy++MY+IUEikWLhHx7CW+u+v6yQ4FRZhOBf1hn3m8PGoHL6fLPWTERERMNLJBYj6cFn8e667zEjvGcdqxatDv+3YSPevH0OVDmnBrgFsmYMCnaorKwMrz77CFavWI7mhp5VsEfFjsWbn32Dux95GrLLmLderOlA+KZ3IdZ1zY7U6eyGcwseBTg7AhERkU1QjhyFp9dvwwsPPAQ3h54OKQfLqvDk/Utw9sv/ClgdDRV+krMjarUa7777LmbNmoW0lD3Gy2WOcvzhyefx+kdfYWRk9GXfbshPn0JeVwEAMECEcwuWWcw0qERERGQeIpEIifc/jX+t/wEzR4UZL2/QaPH826ux+el7oWtvFbBCMjcGBTuxd+9eXHPNNXj99dfR3t6zkur4xKn415dbcMuS+4yL210Oz5P74XUqxbhfOWM+VOFxZqmZiIiILI97SDieWrcVz923FHJJz0fJzw6m4pVF16I564SA1ZE5MSjYuLKyMjzwwANYsmQJCgoKjJe7uXvg8RWvYuU7H8M/sL9hjxcnVrcj6Nf1xn1VSDTKr1o46JqJiIjIsolEIkx9+Dms/s+niHDvGbR8vLoeTz5wN85+8Z6A1ZG5cNYjG9Xa2ooPP/wQ7777Ljo6OoyXi8Vi3Hvvvbj6rkfg4qYc1H34HfwJDm1dC+npHBzPj0vgvMpERET2YkRCEl7ZtANrn30AP53oGtTc2KnF8++8g+y0I7jpjQ8hdpAJUhtnXRo8tijYmPb2drz//vuYMmUKVq1a1SckTJo0CT///DNeeeWVQYcEqaoBfoe2GverptyATuWVLcJERERE1svRVYn7P9iIP97/YJ+uSJ8ePIKP7rsFOlWjgNXRYDAo2Ai1Wo1PPvkEU6dOxcsvv4y6ujrjMW9vb7z11lv47rvvEBdnnvEDAXu+haSza3n3ToUbqqbeaJbbJSIiIus05YFnsPq9/yHCw8142c+55/DG3TdBXVEsYGV0pRgUrFxHRwfWrl2LadOmYcWKFaiurjYec3JywrJly7Bv3z4sWrQIYjNNVyqvKYP38T3G/YpZt0Hv6GSW2yYiIiLrNSI+ES9v2Ibk0CDjZamVtfjr7xdAlXlcwMroSnCMgpXKyMjAV199hW+//RaNjX2b9ORyOe655x4sW7YMPj4+Zr/vwJ0bIDLoAQAdnn6omTjH7PdBROYn5Aq1RGQ/nJSeeGbdj/j0sXvx88mucQtnm1rw/KO/x0t//z/4zLxe4ArpUjEoWBGVSoXNmzfjyy+/xMmTJ02OOzg44K677sLjjz8Of3//IanBuTgH7jnHjPtlcxYDEv4ZERERUQ8HmRwPfLAB3iuewNrt2wEAFe0aLH/hWfz9qXKMvH2pwBXSpeAnPAum0+mQlZWFw4cP4/Dhw9izZ0+fNRC6KRQKLFiwAE888QQCAwOHriCDAUHbe6ZDbQ2MQGNs4tDdHxEREVktkUiEBf94F14jXsM7n38GPQCVVofnV7+Jf3S0I/Kex4QukS6CQcEC6HQ6NDY2oq6uDjU1NTh16hQOHTqEo0ePorm5+YLXS0hIwJIlSzB//ny4uLgMeZ3u2WlwKc0z7pdeeycgEg35/RIREZH1mrXsBbiPCMJrq16DWq+HRm/AyjXv4U0vXwTcsEjo8mgADArDrL6+Hn/605+gUqlQU1OD6upqNDQ0QK/XX9L13d3dsXDhQtx5550YPXr0EFfbi06LwJ0bjLuNUQloCR3G+yciIiKrNX7BPXjFZwRW/ulJtGp1aNfpseK1v+NNT294JHOso6ViUBhmUqkUP/300yWfL5FIMG7cOCQlJSE5ORkzZsyAXC4fwgr7552+B/K6CgCAQSRC2dV3DHsNREREZL0iZlyDv6z4B1b+7c/oNBhQr9Fi5QvP4PX3P4ciepzQ5VE/GBSGmaurKxwcHNDZ2WlyzMHBAV5eXggJCUFSUhKmTJmCiRMnDku3ooGIOtUI2Putcb8ufhY6fIMGuAYRERGRqdHzbsOztdV47d23AAAlbR145Yml+Nunm+AQECpwdfRbDArDTCQSYcWKFfD19YWPjw+kUim8vLzg5eUFV1dXiCywz7/XqQNwaG0CAOilMpRftVDgioiIiMhaJd79EB6pqcKar7omSMlsVOGtR+7Cs5//ALHSU+DqqDcGBQEsXbrUuL5BTU2NwNVchMEAn6Pbjbu18TPR6cYnMREREV25uU//FQ01Vfhq504AwKHKWnzyyB24/5PNgFxhlvsQcu0YpUqM6eFegt2/uXBlZhqQc8lZKKp6ll2vmXytgNUQERGRrbj9H+9ibvx44/5P+cXYsvw+wGAQsCrqjUGBBuTbqzVBNXI0xyYQERGRWYjFYtz/7udIDOsZm/BJ2kkUffm+gFVRb+x6JKCUgjo0NQnXLHYx0pYmuGemGver2ZpAREREZuQgc8RTH27CXxddi7z6RgDAG/9dg7dnzIVDcITA1RFbFOiCvNN3QazXAQA0rh5ojJ4ocEVERERka5xcXfHEq/+Cw/kJXcrbNVj/52XAJa4xRUOHQYH6p9fB59gu427txDmAhA1QREREZH7BCUm497aeWRU3ny1E9vuvCVgRAex6RBfgnpMOWXM9AEAvlqBmIldNJPMQchYKIiKyXNcv/xtSD6fgdHklAODNL7/Ef2bfCMfR8QJXZr/YokD96j0lamNsIrQu7gJWQ0RERLZOIpVi2ao1cJJ0fTyt02jx3xVPAZ1qgSuzXwwKZEJeUwa3cxnGfQ5iJiIiouHgN2o07r/3/xn3d5VUIv3NvwhYkX1jUCATPkd3GLfb/ELQGhwlYDVERERkT2Y/tBxJEWHG/bd+3IqWI7sFrMh+MShQH2J1O7xO7jPu10y+Fjg/CwERERHRUBOJRHhg1QdwkzkAAFq0erz78p+B9haBK7M/DArUh+epA5BoOgAAWkcF6sdOFbgiIiIisjdegSF4ZNnTxv3UmgYc/w9nQRpuDArUw2DosxJzXfxM6GVyAQsiIiIie5W0+D7Mjos17n/6048wNNYKWJH9YVAgI5eibDjVlBr3ayZfI2A1REREZM9EIhEW/+U1SM93gS5pUyP1368IXJV9YVAgI+8Te43bzeFxUHuNELAaIiIisnd+4dG4MTnJuP/p9h3QV5cLWJF9YVAgAIBI2wn37DTjfu2E2QJWQ0RERNTl5uf+Dpm46yNrVUcn9v/rJYErsh8MCgQAcCs4DYm6HQCgc3BEUyRXQSQiIiLheQWG4JZZs4z7n+87AF1JgYAV2Q8GBQIAeGQcMW43RcZzEDMRERFZjBuX/w0KqQQAUK/RYtfbbFUYDgwK1NXtKOeYcb9hTNIAZxMRERENL6WPHxbOnWfcX3s4Ddq8MwJWZB8YFAhu+ex2RERERJbtuqdWGBdhU2l1+PmtlwWuyPYxKBA8Mnt1O4qKh8HBUcBqiIiIiEw5K91x+/xbjfvr009BfTpVwIpsH4OCnTPpdhTLbkdERERkma5+9Hl4ybu+0OzQG/DDv14VuCLbxqBg59jtiIiIiKyFk7MLFi9aYtzflJGDtiO7BazItjEo2DmPzMPGbXY7IiIiIks3a+mT8HNWAAA0egM2v/9PgSuyXQwKdkyk1cA9J924z25HREREZOlkcjnuvPs+4/6POQXQFp0VsCLbxaBgx9jtiIiIiKzRlCUPQOnYNQNSu06Po2vXCFyRbWJQsGOc7YiIiIiskUwux7yp04z7W/btAzrVAlZkmxgU7JRIq4F7du/ZjpIFrIaIiIjo8sz8w2PG7aymVlRu3ShgNbaJQcFOueWfhkTTAaC729F4gSsiIiIiunQBMXGYEBxo3N/x7VcCVmObGBTslEdG725HCex2RERERFbn2gV3GLe35p7joGYzY1CwQ12zHXGRNSIiIrJuExbc02dQc9rnHNRsTgwKdsgt7xS7HREREZHV6xrUPN24/8N+Dmo2JwYFO9R3tiN2OyIiIiLrNfO+voOaqzio2WwYFOyMSKflImtERERkMwKix2BCSJBxfzsHNZsNg4KdcS7JNXY70ksd2O2IiIiIrN61t/Ud1KwrzBGwGtvBoGBn3PJPGbdVoTHsdkRERERWb8KCu6F0lAE4v1IzBzWbBYOCnekdFJojxglYCREREZF5yORyXDetZ6XmH1L2AxoOah4sBgU7Im1tgnNFoXGfQYGIiIhsxazfDmr+eYOA1dgGBgU74pZ/xritcfVAh0/gAGcTERERWY+AqN8OamZQGCwGBTti0u1IJBKwGiIiIiLzunbhncbtrbnnoCspELAa68egYC8MerjlnzbuNkeMFbAYIiIiIvObcOtdUMp6VmrO/P4LgSuybgwKdsKpqhgOrU0AAANEaA6PE7giIiIiIvOSyeWYMa7ny9CDB/YJWI31k5rzxjQaDTZv3owDBw6gtrYWLi4uGD9+PBYvXgwvL6/Luq0TJ05g69atyM/PR1tbG5ydnTFq1CjceOONGDuW34ZfLre8ntaEtsBw6BSuAlZDRERENDQmXXcLfkzrWlx2f1EZHqytgMh7hMBVWSeztShoNBq8/PLL2LRpEzo6OjBp0iR4eXlhz549eP7551FZWXnJt/Xjjz/i1VdfxcmTJxEYGIikpCT4+PggPT0dL7/8Mn799VdzlW03eo9PaOJsR0RERGSjRs+9BQqpBADQqtMjbwtXar5SZmtR+O6775CTk4OoqCi8+OKLkMvlALo+9H/++edYs2YNXnrppYveTnNzM9avXw+pVIq//vWviImJMR47fPgw3nrrLaxduxYzZ8403gcNTKzpgEtxzwqFHJ9AREREtkomlyM5Ogq7MrIAAIf27kTkfU8LXJV1MkuLglarxbZt2wAAS5cu7fMB/qabbkJoaCiysrJQUHDxkee5ubnQarWIi4vrExIAIDk5GSEhIVCr1SgtLTVH6XbBtTATYr0OAKB1VKA1aJTAFRERERENnaRrbjBu780vBFSNwhVjxcwSFLKzs9Ha2go/Pz+EhYWZHE9KSgIApKWlXfS2HBwcLuk+XVxcLq9IO9Z7fIIqfAwglghYDREREdHQGnPzIjiIu6aBr9doUfzzJoErsk5mCQpFRUUA0G9IAIDw8PA+5w0kIiICCoUCZ86cQXZ2dp9jR44cQXFxMaKjo+Hv7z/Iqu2HW0HvaVE5PoGIiIhsm7OrEpPDRhr3j+zcJlgt1swsYxRqa2sB4IIzG3l6evY5byDOzs54+OGH8c4772DlypWIjo6Gp6cnqqurkZ+fj/j4eDz66KPmKNsuyBprIK+rMO5zfAIRERHZg6SrrsXB/P8CAPbmnMUidTsMjk4CV2VdzBIUOjo6AACOjo79Hu8es9B93sUkJyfDxcUFb731Vp9WBaVSiTFjxsDV9dKn9nzmmWdMLpPJZHj99dcBAN7e3pd8W+YklUohkUigVCqH9H7cMg4atzW+QXAKHQU+RYQjkXR1+xrqx50sDx97+8XH3n7xsRfWtHsewtsf/xcGAGXtGjQe3I6QW+8ZlvuWSCSQSqXw8fEZlvsbKmYJCgaDYVDHf2vLli1Yt24dJk+ejEWLFsHX1xfV1dXYsGED1q1bh9zcXCxfvnwwJdsNRU66cbstOkHASoiIiIiGj6d/ABKCA5FeUgYAOLjlm2ELCrbCLEHByanrO2q1Wt3v8e7LL2U608zMTKxduxZhYWF45plnIBZ3DaMICQnB8uXL8cILL+DIkSM4efIkxo8ff9HbW7169YDHa2trLzvImIOPjw90Oh2ampqG7k50WoSdPWHcrQmOQfNQ3h9dVPe3SkP6uJNF4mNvv/jY2y8+9sJLmjYT6V99CQDYdeoUbqivg0Fi1vWG+6VUKqHValFTUzPk99WbSCTCiBHmW1zOLIOZu7vv1NXV9Xu8vr6+z3kD2bt3L4CumZK6Q0I3sViMxMREAEBGRsYV12svnMvyIVG3AwD0EilaQmMucg0iIiIi25Fw293G7TxVOxoObBewGutjlqAQGhoKADh37ly/x7vXT+g+byDdoaK7leK3ui9vaWm57DrtjTKvZzXmlpBo6GVcoI6IiIjsh8/ICMT49Ey2k7b1WwGrsT5mCQoxMTFQKBSoqqrqNywcOXIEADBhwoSL3lZ3M11+fn6/x7sv9/X1vdJy7YZbfk9Q4LSoREREZI+mJE8zbu8/cQLQ6wWsxrqYJShIpVLMmzcPAPDJJ5/0md3oxx9/RFFREWJiYjBqVM+KwNu2bcNTTz2F9evX97mt7q5FKSkpJgu0HT16FCkpKRCJRMbzqH+SNhUU5T2hrXkUgwIRERHZn4Tblhi3TzWo0Hri4ABnU29mG82xYMECnD59Gjk5OXjyyScRExOD2tpa5ObmwtXV1WTtg+bmZpSXl6OhoaHP5ZMnT0ZycjIOHz6MN954AxEREfDx8UFNTY2xNeGOO+5AQECAuUq3SW4FZyBC1yDtThd3tPsGC1wRERER0fALih2PEKUriptUAIDjP36N6ROmC1yVdTBbUJDJZFi5ciW+++47pKSk4OjRo3B2dsasWbOwePHiS16vQCQS4emnn8bu3buxd+9eFBcXo7CwEAqFAgkJCbj++usRHx9vrrJtlktRz/oTzeFjAJFIwGqIiIiIhCESiTB9UiLW79wJADiQdhTTDQZ+NroEZp0fSiaTYfHixVi8ePFFz120aBEWLVrU7zGRSIQ5c+Zgzpw55izPrrgW9wSFlhDOdkRERET2K2H+ImNQSK2qhyb3NGRR7JZ9MWYZo0CWRdLeAqfqUuO+itOiEhERkR0LT5wBH0XX7I96AKe//0rYgqwEg4INcik+a9zuVLhB7WW+hTeIiIiIrI1YLMb08T1d14+mHxWwGuvBoGCDXPp0O4piHzwiIiKye+NmXWPcTi2pADraBKzGOjAo2CCX4hzjNldjJiIiIgIi59xk/ODb0KlF1f5fBK3HGjAo2BhRpxrOvdZP4EBmIiIiIsBZ6Y5Y/54FezP37RCwGuvAoGBjnEvzINLrAAA6mRxt/iECV0RERERkGRLG9YxTOJ5xRsBKrAODgo1x7bV+QktwJCCWCFgNERERkeWInTXXuJ1WWQNDY52A1Vg+BgUb02d8Qki0gJUQERERWZbw6dfAUdL18bddp0fJzi0CV2TZGBRsiU4L59I84y7HJxARERH1kMnlGB8UaNw/c2ivgNVYPgYFG6KoLIKkUw0A0EukaA2MELgiIiIiIssSP2GScTs9O3uAM4lBwYa49Bqf0BYQDoODTMBqiIiIiCxPzJwbjdsnaxugLS8SsBrLxqBgQ1x7jU9QcXwCERERkYmQiVPgJnMAAGgNQMH2HwSuyHIxKNgKg54LrRERERFdhEQiwYSwkcb906kpgtVi6RgUbIS8thzS9hYAgAEitAZHClwRERERkWUaN3mKcTstNw8wGASsxnIxKNiI3uMT2v2CoZM7C1gNERERkeWKuXa+cTu7qRUdZ08JWI3lYlCwEX3XT2C3IyIiIqIL8Y8eA1+F3Lh/duePAlZjuRgUbIRrUa+BzByfQERERHRBIpEIE6OijPsnj6UKWI3lYlCwAbLGGsiae5YgbwnljEdEREREAxmXPMO4fbSgENDrhCvGQjEo2IDe3Y46PP2gdXEXsBoiIiIiyxc9t2ecQkmbGs3pBwWsxjIxKNgAlyKOTyAiIiK6HB6BoQhzdzPuZ+3eKmA1lolBwQa4FPfMeNTChdaIiIiILknC6Fjj9skT6QJWYpkYFKyctLUZTrXlxn2OTyAiIiK6NHHTrjJupxaXwaBuF64YC8SgYOV6j0/QuLhD7eEnYDVERERE1iPympshOr9dp9Gi5sB2QeuxNAwKVq7P+gmh0YBINMDZRERERNTNxcMLsb7exv3MfTsErMbyMChYOS60RkRERHTlEsaOM24fP80VmntjULBiIq0GTpVFxv2W4EgBqyEiIiKyPrEzrzVuH6uohqG1WcBqLAuDghVTVBZBfH5xEL3UAe2+wQJXRERERGRdwmfNhfR81+1WnR5Vh/cIW5AFYVCwYs6l+cbtthFhgEQqYDVERERE1sfRyRkxvl7G/fwj+wWsxrIwKFgx57I843ZrYISAlRARERFZr9jIKON2RtYZASuxLAwKVsy5rKdFoTVolICVEBEREVmvqInJxu1TpeWAXi9gNZaDQcFKSdpUcGyoNu6zRYGIiIjoyoTNvM64XdamRkv2CQGrsRwMClaqd2tCp7MbNErvAc4mIiIiogvxDApFgLPCuH8uZaeA1VgOBgUr5Vzae3zCKC60RkRERDQIcWGhxu2cU+kCVmI5GBSsVJ/xCex2RERERDQoMWMTjNunCwoErMRyMChYI4OBA5mJiIiIzChi2hzjdmZ9E3T1VQJWYxkYFKyQY30lpB2txv3WgHABqyEiIiKyfoHxSXCSdH001hqAkn3bBa5IeAwKVqh3a0K7dwD0csUAZxMRERHRxUgcHBAX4G/cz0s7JGA1loFBwQr1XpGZ4xOIiIiIzCN29BjjdubZbAErsQwMClao74rMHJ9AREREZA6Rk6cZt0+UVwGdGgGrER6DgpURaTVwqiwy7rcGsUWBiIiIyBxGzrjWuN3QqUV9+gEBqxEeg4KVUVQWQazXAQD0Uge0+wYLXBERERGRbXD28MIoD6VxP//QXgGrER6DgpXpPT6hbcRIQCIVrhgiIiIiGxMX0dNbI/vMCQErER6DgpVRlPceyMzxCURERETmFBU/ybh9uqgEMBgErEZYDApWhjMeEREREQ2d8F7jFHKbW6Eusd9VmhkUrIikTQV5Q88qgVyRmYiIiMi8fKPj4OEoM+4X7v9VwGqExaBgRXovtNbp7AaN0lvAaoiIiIhsj0gkwriQIOP+2eNHBaxGWAwKVsS5tPf6CRGASCRgNURERES2afSYccbtjLxcASsRFoOCFendosDxCURERERDIyJ5lnH7VHUd9G0qAasRDoOCtTAYfhMUOD6BiIiIaCgEJ8+Cw/meG+06PaoO7hK4ImEwKFgJx/oqSDtajfttgeECVkNERERkuxydFBjt2zMWNC81RcBqhMOgYCWcy3rGJ7R7B0AndxawGiIiIiLbFhsVZdzOzDwjYCXCYVCwEn1WZOb4BCIiIqIhFTVxinH7ZGk5YNALWI0wGBSsRO8WBQ5kJiIiIhpaI2f2LLxW0a6GKvOEgNUIg0HBCoi0GjhVFhn3OZCZiIiIaGh5BIYi0EVh3C+wwwHNDApWwKmyGGK9DgCglzqgzS9Y4IqIiIiIbF9caKhxO+8MWxTIAjlXnDNut/mHAhKpgNUQERER2YfI0WOM29mFhcIVIhAGBSugKO8VFEaECVgJERERkf0Y2WtAc2ZNPfTqdgGrGX4MClZAUcGgQERERDTcApNmGj8st+n0qD12UNB6hhuDgoUTdWrgVFNm3G8LYFAgIiIiGg5yZ1eM8lQa94vSGBTIgjhVFUPUayBzu0+gwBURERER2Y+Y0JHG7bys08IVIgCzjYrVaDTYvHkzDhw4gNraWri4uGD8+PFYvHgxvLy8Lvv2KisrsXnzZpw+fRqNjY2Qy+UYMWIEEhMTMX/+fHOVbfFMBjKLJQJWQ0RERGRfImLjgOMnAQBZRcUCVzO8zNKioNFo8PLLL2PTpk3o6OjApEmT4OXlhT179uD5559HZWXlZd1eamoqli9fjt27d8PFxQWTJ09GWFgYqqursX37dnOUbDU4kJmIiIhIOKGTpxu3cxqaoGtpFrCa4WWWFoXvvvsOOTk5iIqKwosvvgi5XA4A+PHHH/H5559jzZo1eOmlly7ptgoLC/H222/DyckJzz33HGJiYozH9Ho9zp07N8C1bQ8HMhMREREJJ2DiFMjEImj0Bmj0BlSl7kPAnJuELmtYDLpFQavVYtu2bQCApUuXGkMCANx0000IDQ1FVlYWCgoKLun2Pv30U2i1Wjz66KN9QgIAiMViREREDLZkqyHSciAzERERkZAcZHJEeXsa9wuPHRKwmuE16KCQnZ2N1tZW+Pn5ISzM9INsUlISACAtLe2it1VaWoqsrCyMGDECEydOHGxpVs+pqoQDmYmIiIgEFtPrM25uTqaAlQyvQXc9KioqAoB+QwIAhIeH9zlvIGfOnAEAjKge/C8AACAASURBVBs3DhqNBgcPHjS2RISGhmLKlClQKBSDLdlqOJf3tMK0+YVwIDMRERGRACLi4oEjXV96Z5WUClzN8Bl0UKitrQWAC85s5Onp2ee8gZSUlAAAZDIZ/vjHP6K8vLzP8fXr12P58uWIjY0dTMlWo89A5oBwASshIiIisl+hiTOBjz8CAOQ2qtBZVw0HL1+Bqxp6gw4KHR0dAABHR8d+j3ePWeg+byCtra0AgK1bt8LZ2RnPPvss4uLi0NjYiE2bNiElJQWrVq3C6tWr4eHhcUn1PfPMMyaXyWQyvP766wAAb2/vS7odc5NKpZBIJFAqlRc8x7W6ZwouQ0TsgOeSdZBIulqF+FjaHz729ouPvf3iY287XKbNhkIiQZtOBz2A+hOHELXg3gueL5FIIJVK4ePjM3xFDoFBj1EwGAyDOt6bXq8HAOh0Ojz++ONITEyEQqFAQEAAnnjiCURERKC1tRW//PLLoGq2BqJODWSVPUFBHTRKwGqIiIiI7JdEKkWsf08LQu5R+xjQPOgWBScnJwCAWq3u93j35b1nQ7qQ7nM8PT0xfvx4k+OzZ89Gfn4+MjIyLrm+1atXD3i8trb2ssKMufj4+ECn06Gpqanf44qy/D4DmWuclMAFziXr0f2t0oUed7JdfOztFx97+8XH3rZEhYcjrawCAJCZeRrTB3hclUoltFotampqhqs8AIBIJMKIESPMdnuDblHo7rpTV1fX7/H6+vo+5w3E17crqV2omab78uZm21/oos/4BA5kJiIiIhJUxNgJxu3MsvIBzrQdgw4KoaGhAHDBhdB6z1p0MSNHjgQAtLS09HtcpVIBuLTWCWvnzIXWiIiIiCxGyJSrjNtFLe1Ql198Rk9rN+igEBMTA4VCgaqqqn7DwpEjRwAAEyZM+P/t3XtwVPX9//HXsrmHECBACKIhJMKWgkGuWviB/dbRoEyHixB0xksH7SjactFWWxyRIkXr/LDtdIaxRRShaBGFUX40Oq1CCX7lFtLG3Ay5ESQJWSAEg0nYZH9/2F1ZTwJJzglnkzwff53dc/LZ9847gfPas5/PMez7rnHjxik8PFxVVVWtrpKUl/fNurVtLcXak0RdvjQqKx4BAADYKu7GMRoQFup/fOJ/P7GxmmvDdFAICQlRWlqaJGnTpk0Bqxvt3r1b5eXlcrlcSkn5djJuRkaGli1bpm3btgWMFR4erlmzZqm5uVkbN24MGCs7O1v79u2Tw+HQ7bffbrbsoMYdmQEAAIKLw+HQ9xLi/Y9Ls69+M+HuzvRkZkmaN2+ecnJyVFhYqKVLl8rlcsntdquoqEgxMTFasmRJwPF1dXU6deqUzp07ZxjrnnvuUX5+vrKysrR06VKlpKSorq5OX3zxhbxerxYtWhQQOnoi7sgMAAAQfEbdOEqfln9zw7XC41/oDpvr6WqWBIWwsDCtWrVKO3fuVGZmpg4fPqzo6GjNnDlT6enpHbpXgW+s999/X/v371d2drZCQ0M1duxY3X333e36ClN3x0RmAACA4DMydZL0j48lSXmVVZLXKzkcNlfVdSwJCtI3J/jp6elKT0+/6rELFy7UwoUL2y4qJETz5s3TvHnzrCqvW2EiMwAAQPC54Qf/I/3f30mSqr5u0lfF+eqbMsbmqrqO6TkKsF7U5UGB+QkAAABBIXb4CMVHfbv65onP9tlYTdcjKAQZh6dJkadP+h9zRQEAACB4jLlumH+7JCfLxkq6HkEhyDCRGQAAIHjdOMrl3y4sKbaxkq5HUAgyhonMTsumkQAAAMCkkTdP8W/nVp2W978f8PZEBIUgw0RmAACA4HX9D/7Hv32uyaPzucdsrKZrERSCDBOZAQAAgldU3BAlxkT7H5cd+peN1XQtgkIQYSIzAABA8Pve9cP92yWf/9vGSroWQSGIBExkdjKRGQAAIBjd6Pr23gkFpSU2VtK1CApBJOBrR0OZyAwAABCMEifc4t/OrzmrFs8lG6vpOgSFIBKw4hFfOwIAAAhK102ZIcd/t7/yNOvc50dtraerEBSCSFRVmX+boAAAABCcImIHaES/vv7HJw4fsLGarkNQCBKOZo8iqyv8jy8mjLCvGAAAAFzR6OHfziUtzc+xsZKuQ1AIEhE1J9XHN5G5j1MNQ4Zf5ScAAABgl+TL7tD8RVnpFY7svggKQSKqssy//XX89fIykRkAACBojRg/2b+df9rdI+/QTFAIEpcHhYtDR9hVBgAAANrhuqkz/dt1l5p1rgfeoZmgECQCrigwPwEAACCoRQ4cHHCH5oojmTZW0zUICsGgpUWR1Sf8D5nIDAAAEPxGXzfMv12a+x8bK+kaBIUgEHGmUs5LjZIkr8Ohi/E32FwRAAAAriblxtH+7S9Ke96EZoJCELj8a0cNg66TNzTctloAAADQPok3TfJvF5yukbelxcZqrEdQCAJRlZffkXmEfYUAAACg3a675dsJzeeaPDpf2LPup0BQCAIBd2RmxSMAAIBuIXpIgob3jfQ/PnG4Z01oJijYzetVZGW5/yFXFAAAALoP17BvJzSX5WbbWIn1CAo2C6utUUjjRf/ji0MTbawGAAAAHZGSMsq/XVRSYmMl1iMo2Ozy+QkNA+PVEhFlYzUAAADoiMSbJvq386tP96gJzQQFmwXckZmvHQEAAHQrwy+b0Hym8ZIuFOfbWI21CAo2CwgKTGQGAADoVvomXK+EqAj/44oeNKGZoGAnr5crCgAAAN2ca9hQ/3ZpzjEbK7EWQcFGzvNnFHqxzv/4a64oAAAAdDs3BkxoLraxEmsRFGwU/uW3v0hN/eLkie5nYzUAAADojBHjJvi3C6qqJa/XxmqsQ1CwUfjJb4MCXzsCAADonoZPneHfPt3QpNriAhursQ5BwUbhX3671i73TwAAAOieYq4fqfjIcP/j0k8/trEa6xAUbBR4RSHJxkoAAABghish3r99POuQjZVYh6Bgk5bz5xRaW+N/zFePAAAAuq+UkSn+7cLjX9hYiXUICja5VPLtL9Cl6H66FDPAxmoAAABgxuUTmvNPVdpYiXUICja5VFLo3744dITkcNhXDAAAAEy5/rIJzZUXG1VZVHiFo7sHgoJNPJddUeBrRwAAAN1bvxE3alBEmP/x0Q//bmM11iAo2CTgigJBAQAAoFtzOBz63tAh/sdH//eAjdVYg6BgA+/FejVXnvQ/ZsUjAACA7u/yCc3Zubk2VmINgoIdKkr9m57wKDX1H2xjMQAAALBC0tjx/u1/l39pYyXWICjYwFvx7f0Tvk5IZCIzAABADzB8yv/xb5+s/1pnK8ptrMY8goIdyi+/I/MI++oAAACAZfonf08Dw0L9j3P2/tPGaswLsbuA3sh7gjsyAwAA9DSOPn30vaGDdeLsOY1PvF7RMbF2l2QKQeEa8zY2SgETmUfYVwwAAAAstXzr/9PAIUM1fWScampq7C7HFL56dK1VnZS8LZKklrBwNcQl2FwQAAAArOIMj7S7BMtwReEacyQmq88ftim27pzyC0qlPmQ1AAAABB+Cgg0cUX0Vnpikr6ITpfPn7S4HAAAAMODjbAAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYhVg3U1NSkXbt26cCBA3K73erbt69SU1OVnp6uuLi4To9bWVmpp556SpcuXVJqaqpWrlxpVckAAAAA2mDJFYWmpiatWbNGO3bsUENDgyZNmqS4uDjt3btXTz/9tKqqqjo99p///Gd5PB4rygQAAADQTpYEhZ07d6qwsFCjRo3SH/7wBy1fvly//e1v9cADD6iurk4bNmzo1Lgff/yxcnNz9aMf/ciKMgEAAAC0k+mg4PF4lJGRIUlavHixIiIi/Ptmz56txMRE5efnq6SkpEPjnj9/Xlu2bNG4ceM0bdo0s2UCAAAA6ADTQaGgoED19fWKj49XUlKSYf/UqVMlSUeOHOnQuK+//rqampr0yCOPmC0RAAAAQAeZDgrl5eWS1GpIkKSRI0cGHNceWVlZ+vTTTzV37lwNHTrUbIkAAAAAOsh0UHC73ZLU5spGAwcODDjuahoaGvTaa69p2LBhmjNnjtnyAAAAAHSC6eVRGxoaJEnh4eGt7vfNWfAddzVvv/22ampq9NxzzykkxPzqrStWrDA8FxYWphdffFGSNGjQINOv0RkhISFyOp2KjY215fVhD6fTKUn0vRei970Xve+96H3v5XQ6FRISosGDB9tdiimmryh4vV5T+y9XXFysjIwMzZgxQ2PHjjVbGgAAAIBOMv2RfWRkpCSpsbGx1f2+5y9fDak1zc3NevXVVxUVFaUHHnjAbFl+69evv+J+t9vdoTBjlcGDB6u5uVnnz5+/5q8N+/g+VaLvvQ+9773ofe9F73uv2NhYeTwe1dTUXNPXdTgcSkhIsGw800HB99WdM2fOtLr/7NmzAce15cyZMyorK1P//v0NJ/f19fWSpOPHj+v5559XRESEnnnmGbOlAwAAAGiD6aCQmJgoSSotLW11v+/+Cb7jrqa2tla1tbWt7quvr1deXp6ioqI6USkAAACA9jIdFFwul6KiolRdXa3S0lLDMqkHDx6UJE2YMOGK4wwZMkTbt29vdV9ubq5Wr16t1NRUrVy50mzJAAAAAK7C9GTmkJAQpaWlSZI2bdoUsLrR7t27VV5eLpfLpZSUFP/zGRkZWrZsmbZt22b25QEAAAB0AfPrj0qaN2+ecnJyVFhYqKVLl8rlcsntdquoqEgxMTFasmRJwPF1dXU6deqUzp07Z8XLAwAAALCY6SsK0jf3JVi1apXmz5+vsLAwHT58WKdPn9bMmTP10ksvcXdlAAAAoJtxeO1YGzSIVFZW2rY8ambJGZZM62VYKq/3ove9F73vveh97xUbG6vpI+O6/fKollxRAAAAANCzEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABiEWDVQU1OTdu3apQMHDsjtdqtv375KTU1Venq64uLi2jVGfX29jh07pqNHj6qsrExut1sOh0PDhw/X9OnTdccddygkxLKSAQAAALTBkrPupqYmrVmzRoWFhRowYIAmTZqkmpoa7d27V1lZWXrhhRc0dOjQq47zwQcf6L333pPD4VBSUpImTpyouro6FRYW6vjx4/rss8+0cuVKhYeHW1E2AAAAgDZYEhR27typwsJCjRo1Ss8++6wiIiIkSbt379abb76pDRs2aPXq1VcdJyIiQnPnztWdd96pgQMH+p+vrKzUmjVrVFBQoHfffVf33XefFWUDAAAAaIPpOQoej0cZGRmSpMWLF/tDgiTNnj1biYmJys/PV0lJyVXHmjNnju69996AkCBJCQkJ/nBw4MABsyUDAAAAuArTQaGgoED19fWKj49XUlKSYf/UqVMlSUeOHDH1OiNGjJAknTt3ztQ4AAAAAK7OdFAoLy+XpFZDgiSNHDky4LjOqq6uliT179/f1DgAAAAArs50UHC73ZLU5spGvq8R+Y7rrD179kiSJk2aZGocAAAAAFdnejJzQ0ODJLW5EpFvzoLvuM746KOPlJOTo+joaM2ZM6dDP7tixQrDc2FhYXrxxRclSYMGDep0XWaEhITI6XQqNjbWlteHPZxOpyTR916I3vde9L73ove9l9PpVEhIiAYPHmx3KaaYvqLg9XpN7b+avLw8vfHGG3I4HHrssccME50BAAAAWM/0FYXIyEhJUmNjY6v7fc9fvhpSe5WXl+vll1+Wx+PRT37yE02ZMqXDY6xfv/6K+91ut+kw0xmDBw9Wc3Ozzp8/f81fG/bxfapE33sfet970fvei973XrGxsfJ4PKqpqbmmr+twOJSQkGDZeKavKPi+unPmzJlW9589ezbguPaqqqrS2rVrVV9frwULFmjWrFnmCgUAAADQbqaDQmJioiSptLS01f2++yf4jmuPs2fP6oUXXlBtba3uuusuLViwwGyZAAAAADrAdFBwuVyKiopSdXV1q2Hh4MGDkqQJEya0a7yvvvpKa9eu1enTp3XbbbfpwQcfNFsiAAAAgA4yHRRCQkKUlpYmSdq0aVPA6ka7d+9WeXm5XC6XUlJS/M9nZGRo2bJl2rZtW8BYjY2NWrdunSoqKnTrrbfq0UcflcPhMFsiAAAAgA4yPZlZkubNm6ecnBwVFhZq6dKlcrlccrvdKioqUkxMjJYsWRJwfF1dnU6dOmW4y/Jbb72loqIi9enTR06nUxs2bGj19R5//HErygYAAADQBkuCQlhYmFatWqWdO3cqMzNThw8fVnR0tGbOnKn09PR2T2Sur6+XJLW0tCgzM7PN4wgKAAAAQNdyeO1YGzSIVFZW2rY8ambJGZZM62VYKq/3ove9F73vveh97xUbG6vpI+NYHhUAAABAz0NQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYBBi1UBNTU3atWuXDhw4ILfbrb59+yo1NVXp6emKi4vr0Fj19fV65513dOjQIdXW1qp///6aPHmyFi5cqOjoaKtKBgAAANAGS64oNDU1ac2aNdqxY4caGho0adIkxcXFae/evXr66adVVVXV7rEuXLigX//619qzZ4+cTqcmT56syMhI/f3vf9evfvUrXbhwwYqSAQAAAFyBJVcUdu7cqcLCQo0aNUrPPvusIiIiJEm7d+/Wm2++qQ0bNmj16tXtGmvz5s2qrKzUlClTtHz5cjmdTknSpk2blJGRoc2bN+uJJ56womwAAAAAbTB9RcHj8SgjI0OStHjxYn9IkKTZs2crMTFR+fn5KikpuepYtbW12r9/v5xOpx5++GF/SJCk+++/X/369VNmZqZqa2vNlg0AAADgCkwHhYKCAtXX1ys+Pl5JSUmG/VOnTpUkHTly5KpjHTt2TF6vV2PGjFH//v0D9oWGhmrixIlqaWlRdna22bIBAAAAXIHprx6Vl5dLUqshQZJGjhwZcJyZsZKSkvTJJ5+orKysE5W2zuFwWDZWR/VxOOTsY9/r49rr89/fN/re+9D73ove9170vvfy9f5an2da/Xqmg4Lb7ZakNlc2GjhwYMBx7RnL9zPf5XuN9ozls2LFCsNz0dHRWrNmjSRp6NCh7R7LajNHxUuKt+31YSf63nvR+96L3vde9L63SkhIsLsEU0x/9aihoUGSFB4e3up+35wF33FmxvI939jY2OE6AQAAALSf6SsKXq/X1P7WjrXyssn69estG8tKzzzzjCTpxRdftLkSXEv0vfei970Xve+96H3v1VN6bzooREZGSmr7U37f85evhnS1sdq6+uAbq60rDt1JU1OT3SXABvS996L3vRe9773ofe/VU3pv+qtHgwYNkiSdOXOm1f1nz54NOK49Y/l+5rt8r9GesQAAAAB0numgkJiYKEkqLS1tdb/v/gm+48yM5Xu+PWMBAAAA6DzTQcHlcikqKkrV1dWtnuAfPHhQkjRhwoSrjjV+/Hg5HA7l5+fr/PnzAfsuXbqko0ePyuFw6OabbzZbNgAAAIArMB0UQkJClJaWJknatGlTwPyC3bt3q7y8XC6XSykpKf7nMzIytGzZMm3bti1grAEDBmjatGnyeDzauHGjmpub/fu2bt2quro6TZ8+3XAzNgAAAADWcng7sixRG5qamrR69WoVFRVpwIABcrlccrvdKioqUkxMjNauXRtwv4Lt27drx44dmjlzph5//PGAserq6rRy5UpVV1crPj5eycnJqqioUEVFheLj47V27Vr169fPbMkAAAAArsD0FQVJCgsL06pVqzR//nyFhYXp8OHDOn36tGbOnKmXXnqpQzc169evn9atW6e0tDR5PB4dOnRIFy9eVFpamtatW0dIAAAAAK4BS64oAAAAAOhZLLmiAAAAAKBnISgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMAixu4CeoKmpSbt27dKBAwfkdrvVt29fpaamKj09XXFxcR0aq76+Xu+8844OHTqk2tpa9e/fX5MnT9bChQsVHR3dRe8AnWVF7+vr63Xs2DEdPXpUZWVlcrvdcjgcGj58uKZPn6477rhDISH8qQYbK//uL1dZWamnnnpKly5dUmpqqlauXGlh1bCC1b2vqqrSrl27lJOTo9raWkVERCghIUFTpkzRj3/84y54B+gsK3ufnZ2tPXv2qLi4WBcvXlR0dLRSUlJ09913a9y4cV30DtBRJSUl+s9//qPjx4+rqKhI586dU2hoqP761792arzudp7nfP7555+3u4jurKmpSWvWrNH+/fsVGhqqm266SV6vV1lZWdq/f7+mTJmivn37tmusCxcuaOXKlTp27JhiYmI0duxYff311zp69KgOHjyo6dOnKzw8vIvfEdrLqt6/9957euONN3Ty5EkNGjRIo0aNUnR0tIqLi5WVlaXc3Fz94Ac/ICwEESv/7r9r/fr1qq6uliQNHTpUM2bMsLJ0mGR17w8dOqTf/OY3Kikp0ZAhQzR69GhFR0fr1KlTKi8v11133dWF7wYdYWXvd+/erT/96U+qrq5WUlKSbrzxRnm9XuXl5elf//qXYmNjlZyc3MXvCO3xl7/8RR9++KG+/PJLNTQ0SJKcTqfmz5/f4bG643keZx4m7dy5U4WFhRo1apSeffZZRURESPrmH4E333xTGzZs0OrVq9s11ubNm1VZWakpU6Zo+fLlcjqdkqRNmzYpIyNDmzdv1hNPPNFl7wUdY1XvIyIiNHfuXN15550aOHCg//nKykqtWbNGBQUFevfdd3Xfffd12XtBx1j5d3+5jz/+WLm5ubr99tv1j3/8w+qyYQEre19WVqbf//73ioyM1C9+8Qu5XC7/vpaWFpWWlnbJe0DnWNX7uro6bdu2TSEhIXruuecC+v7ZZ5/plVde0ZYtWzRjxgz/a8A+o0aN0ogRI5ScnKzk5GT99Kc/7fRY3fE8jzkKJng8HmVkZEiSFi9eHPAHPXv2bCUmJio/P18lJSVXHau2tlb79++X0+nUww8/7P/lkaT7779f/fr1U2Zmpmpra61/I+gwK3s/Z84c3XvvvQEhQZISEhL84eDAgQMWVg8zrOz95c6fP68tW7Zo3LhxmjZtmqU1wxpW9/7111+Xx+PRkiVLAk4WJalPnz58ohxErOx9UVGRPB6Pxo4da+j7LbfcohtuuEGNjY06efKktW8CnTJnzhwtXLhQEydOVP/+/Ts9Tnc9zyMomFBQUKD6+nrFx8crKSnJsH/q1KmSpCNHjlx1rGPHjsnr9WrMmDGGX8TQ0FBNnDhRLS0tys7OtqZ4mGJl769kxIgRkqRz586ZGgfW6arev/7662pqatIjjzxiSZ2wnpW9P3nypPLz85WQkKCJEydaXiusZWXvQ0ND2/Wanf36IoJTdz3PIyiYUF5eLkmt/qMhSSNHjgw4zsxYvufLyso6Wia6gJW9vxLfd9XNfIoBa3VF77OysvTpp59q7ty5Gjp0qPki0SWs7P3nn38uSbrpppvU1NSkvXv3atOmTdq0aZP++c9/6uLFixZVDStY2fvk5GRFRUXp888/V0FBQcC+gwcP6sSJExo9ejT/FvQw3fU8jzkKJrjdbklqc6UD31dJfMe1Z6zvfv3Ex/ca7RkLXc/K3l/Jnj17JEmTJk0yNQ6sY3XvGxoa9Nprr2nYsGGaM2eONUWiS1jZ+4qKCklSWFiYfvnLX+rUqVMB+7dt26Ynn3xSY8aMMVMyLGJl76Ojo/Xoo4/qj3/8o1atWqXRo0dr4MCBOn36tIqLizV+/HgtWbLEuuIRFLrreR5BwQTf7Pe2Zqj7vsPoO87MWL7nGxsbO1wnrGdl79vy0UcfKScnR9HR0ZxABhGre//222+rpbwmNwAABj1JREFUpqZGzz33HCtbBTkre19fXy/pmw8DoqOj9dRTT2ns2LGqra3Vjh07lJmZqZdfflnr16/XgAEDLHoH6Cyr/+5vueUW9e3bV6+88krAVYXY2Fh9//vfV0xMjMmKEWy663ke/yuZ4PV6Te1v7ViHw2GqJlwbVva+NXl5eXrjjTfkcDj02GOPtfkJBK49K3tfXFysjIwMzZgxQ2PHjjVbGrqYlb1vaWmRJDU3N+tnP/uZUlNTJUlRUVH6+c9/rsrKShUXF+vDDz/UokWLOl80LGH1v/kffPCBtm7d6l8/f8iQITp9+rT+9re/aevWrSoqKtKTTz5ppmQEme56nsccBRMiIyMltZ3+fM+3Z3kz31htfRrhGyvY1tftrazs/XeVl5fr5Zdflsfj0UMPPaQpU6Z0vlBYzqreNzc369VXX1VUVJQeeOABa4tEl7Dy7953zMCBA/0h4XI//OEPJUm5ubmdqhXWsrL3eXl52rJli0aMGKEVK1bohhtuUEREhG644QY9+eSTSkpK0sGDB/Xvf//bujcA23XX8zyuKJgwaNAgSdKZM2da3X/27NmA49ozlu9nvsv3Gu0ZC13Pyt5frqqqSmvXrlV9fb0WLFigWbNmmSsUlrOq92fOnFFZWZn69++v9evXB+zzfS3l+PHjev755xUREaFnnnnGbOkwycq/+yFDhkiSBg8e3Op+3/N1dXUdrhPWs7L3+/btk/TNSkl9+gR+XtunTx9NmTJFpaWlys3NbTVEonvqrud5BAUTEhMTJanNm+L41lP2HWdmLN/z7RkLXc/K3vucPXtWL7zwgmpra3XXXXdpwYIF5guF5azufW1tbZvrZtfX1ysvL09RUVGdqBRWs7L3vqWPv/rqq1b3X7hwQVLnrkrCelb23nei6PuE+bt8z7f1u4Huqbue5xEUTHC5XIqKilJ1dbVKS0sNS14dPHhQkjRhwoSrjjV+/Hg5HA7l5+fr/Pnzio2N9e+7dOmSjh49KofDoZtvvtnaN4FOsbL30jf/Iaxdu1anT5/WbbfdpgcffNDymmENq3o/ZMgQbd++vdV9ubm5Wr16tVJTU7Vy5UprCodpVv7djxs3TuHh4aqqqpLb7TZ8ipiXlyep7aUUcW1Z2Xvf/+/FxcWt7vc977vqhJ6hu57nMUfBhJCQEKWlpUn65vbbl3/vbPfu3SovL5fL5VJKSor/+YyMDC1btkzbtm0LGGvAgAGaNm2aPB6PNm7cqObmZv++rVu3qq6uTtOnT2c9/SBhZe8bGxu1bt06VVRU6NZbb9Wjjz7a7SY79SZW9h7di5W9Dw8P16xZs9Tc3KyNGzcGjJWdna19+/bJ4XDo9ttv7+J3hfawsve+eWeZmZmGG7QdPnxYmZmZcjgczE/rpnraeR5XFEyaN2+ecnJyVFhYqKVLl8rlcsntdquoqEgxMTGGtZDr6up06tSpVu+0+9BDD6moqEgHDx7UsmXLlJycrIqKClVUVCg+Pp5PmYOMVb1/6623VFRUpD59+sjpdGrDhg2tvt7jjz/eZe8FHWPl3z26Fyt7f8899yg/P19ZWVlaunSpUlJSVFdXpy+++EJer1eLFi0KOPGEvazq/eTJk3XLLbfos88+0+9+9zslJydr8ODBqqmp8V9NWLRokYYNG3bN3hvalpWVpXfffTfgOY/HE3C1d/78+f6rST3tPI+gYFJYWJhWrVqlnTt3KjMzU4cPH1Z0dLRmzpyp9PT0Dk1K6devn9atW6ft27fr8OHDOnTokGJjY5WWlqaFCxdyO/cgY1XvfRNXW1palJmZ2eZxBIXgYeXfPboXK3vvG+v999/X/v37lZ2drdDQUI0dO1Z33313u7+6iGvDqt47HA4tX75cn3zyifbt26cTJ06orKxMUVFRuvnmmzVr1iyNHz++i98N2quurk5FRUUBz3m93oDn2rvoQHc8z3N4zS74DgAAAKDHYY4CAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAg/8PmmakW6BifI8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 900x600 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(x,y)\n",
"plt.plot(x,y_le, c='k')\n",
"plt.hist(data, 20, (0,1), density=True, alpha=0.3)\n",
";"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"''"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAwoAAAIJCAYAAAAF28o7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3Xl4VPW9P/D3mS2ZzCSTfSNkJzuEJYQdAUVRUawLuPXWpdUW6t7q7a9erdcuXtfW2xb1Wm1dsCoqCmLEJYDsSQgQyEL2hex7MslMMpn5/RFzSJwQApnkzPJ+PQ8P5ztz5swnOZnJvHO+i2CxWCwgIiIiIiIaRiZ1AUREREREZH8YFIiIiIiIyAqDAhERERERWWFQICIiIiIiKwwKRERERERkhUGBiIiIiIisMCgQEREREZEVBgUiIiIiIrLCoEBERERERFYYFIiIiIiIyAqDAhERERERWWFQICIiIiIiKwwKRERERERkhUGBiIiIiIisMCgQEREREZEVBgUiIiIiIrLCoEBERERERFYYFIiIiIiIyIpC6gKkVl9fD4vFMuXP6+/vDwBobm6e8ucm6fC8uy6ee9fFc++6eO5dl1TnXhAEBAcH2+x4Lh8ULBaLJEFh+POT6+F5d108966L59518dy7Lkc/9+x6REREREREVhgUiIiIiIjICoMCERERERFZYVAgIiIiIiIrDApERERERGSFQYGIiIiIiKwwKBARERERkRUGBSIiIiIissKgQEREREREVhgUiIiIiIjICoMCERERERFZYVAgIiIiIiIrDApERERERGSFQYGIiIiIiKwwKBARERERkRUGBSIiIiIissKgQEREREREVhgUiIiIiIjICoMCERERERFZUUhdAJGr2VfWgo6OLsmePz3MU7LnJiIiIsfBKwpERERERGSFQYGIiIiIiKwwKBARERERkRUGBSIiIiIissLBzERERJOMkxgQkSPiFQUiIiIiIrLCoEBERERERFYYFIiIiIiIyAqDAhERERERWWFQICIiIiIiKwwKRERERERkhUGBiIiIiIisMCgQEREREZEVBgUiIiIiIrLCoEBERERERFYYFIiIiIiIyAqDAhERERERWWFQICIiIiIiKwwKRERERERkRSF1AURERETOal9ZCzo6uiR7/vQwT8memxwfrygQEREREZEVBgUiIiIiIrLCoEBERERERFY4RoGIiMjBCKY+yHp7YOpoBcwDUHr5wOKuxoDKHZDJpS6PiJwEgwKRHRH0XTCeKUN3/RnoG+rQ2dyIrtZmdLa1Qt+jh7ubGzw8NPDQaKDWaKHWekHt6QXfyBnwSEmD2UMr9ZdARBNkMfXDWFmMplPHUF9cgDNVZahraEBnTy+6jH3o6utDV78JXaYBmCyDjxEAaBQyaORyaJRyaJRKeLq5ISzAH9PDI2GYm4q4eekITJkFmdpD0q+PiBwHgwKRFCwW9NdVoz5nH6pO5qK8tAQltbU43d6F3gHzRR0ywE2Bmf6+SIiORmxyKsIXLINixixYlCobF09EtmI29qI5ex/Kj+zD6fw8lNbWoqKjCx39Axd0HAuAbpMZ3SYzYOwHYADQBdQ3A3mFwOcZAAAvhRwzfDwxM2I6lixbhsXX3QhdTBwEQbD510ZEjk+wWCwWqYuQUl1dHaT4FgQEBAAAmpqapvy5aeqZB0wo3b8Xx3Z/i2++24/jlVWo0hsm/XmjNO5YkZSAZevWw3flNbAolJwqTyJ8zbuugICA76fI7EBPdRnKvv0cJcdzUFBWilONLdBf5B8HbEEGYJafDktSkrDs0kuRtu4GqP0DJavH2Qw/91Lhe740pHrPFwQBISEhtjsegwKDAtnegKkfBZlf4+CXX+BwzlEcqaxBi7F/3I9XygR4u7vBx8MD3lotdF46aLSeMBoN6OnpQU9vL/SGXuiNRnQbjGjpPX/oSNRpcHlaGm67916ELloGQca5DKYSX/Ouqa4wH8e//hK7vv4GOaeLUdLRPe7HeqmUiPD1wfTgYIRMj4QuMAgeOh+463yh9vWHm68/PLx9Icjk6O3qQF9HG4wdbTB0tcPQ0YHulkbUlhWjqqYKFY3NaBrH+4S7TIaV0WG4es0VWH3Hz6ANmTaRL9/lMSi4LgYFJ8GgQLYwMDCAwr27sX/ndhzMzsbh8mp09JvO+ziNQo64QH9ER0QiIj4ZoXPSEZQ4Gxov3QV1BWhvbkJF9n5UZB9AceEpFFTXoM1gHHVfAcCiIF/cvPYqXPur/well27cz0MXj69511Bz6gQObv8Uhw4ewKGiElR09Zz3MXIBmOHng8TIKMSkpCI4aTYCkufA0y/AZl2C0sM80drSguKcIyg4chgHDx7AgcJitBr6zvkYN5mAVTERWHvlGlx258+gDQy2SS2uhEHBdTEoOAkGBWkcqZFulUpg4m+cht4enPz2axzJ/AbZubk4XFaJ9r7zB4NEf2/MmhGP2LkLELF4Jfyi4yGbhL/sWywWNJ6pwolt72HP11/gVF3DqPtFaNyxcd3VuPH/PQl3Hz+b10Fnufpr3hlZLBZUnzyOgzs+xcEDB3HodAmqu3vP+7gAdzekRkxHXGIKoucvxbSFy+Gu9ZrUWkd7zzObzSg4chh7P/sY+w8cxKGyynOOkXKXyXBZfBRuuukmXPKTn0Lprp7Uep0Fg4LrYlBwEgwK0nCkoGCxWFBbVooTmV8he98+ZOcXIK+uEUbz2D83cgGYGeiHBSlJWLh8BdLWXov4mamS/NKoLS7A4ff+gW+/242aTuuuD8HuKvz8qtW47fGn4BFkuzcYeyPlz51Op8PSaD+Xfc0D0n7/bfFhydTfj8J9e3Dkqy9xJCcHWWUVqO8Z/crdcGFaNebGxCAhdS5iV6yBf9KcKR88PJ6v32g0Yt9H72P7Bx/gy2N56DzHVdEAdyWuW5CGm356D5JWruZA6DEwKLguBgUnwaAgDXsNCobeXhRnHUL+oQM4lXcCBaVlyK9rQnvf+ccXyADMDPDF4pnJWLRyFeZfsw5eAUEj9pH6l4bFYoG8NBtv/eUlbMs6ioEf/Oj7qRT42eoVuOvp/4EmyPm6GTAoSMvRgkJHUyOOfbkTOfu/Q87JU8ipqkWX6fyzEUV4emBRfCwWLlqMhWuvxdyVlznch0WjwYDvPnwPO7Z+iC+PnzpnaEjy0+HGyy/Djzbej8DoWFuU6lSkfs8HGBSkwqDgJBgUpCHlBwZjVwcC2stRefIkyk4XoqKiEhW1dahobsGZ7l6M96fBTSZgdrA/0hLikb50KdLX3QCv4NAxH2MPvzSGNJUVYedf/oCdR7LQ94OrI0HuKtz9o+sxb+NvIFO52fR5pfylxaAgLXsOCgMmE4oP7cfRzG9wNCcHOSVlKG7rHNf7QbSXBgvjYrFwyRIsumYdQhNTRtxvD6/7ibzujIZe7H7nX/jw3+/h68JS9I/yO1MuACtiInDTjTdh9V0/g7uGa7oAjn/u6eIxKDgJBgVpTOYHBktPN9qK8tBYdBINlWVoqq9FQ1MT6tvaUNvZjdZxXB0YTYhahTnh05A2axbmr1iJ5EuvgJvnhfUrtodfGj/UVlOJL1/+PT7bt9+qf3KCjyd++rNNiLn+J4ANB1VKhUFBWvYUFBpLi5H7VQaOHj6M3MJCnKhtHNfVArkApPj7YH5iPNIXLUb6VdcgIDZuzMfYw+veVq+71jM12Lb5ZXy0MwPHGlpG3UenUuCatNm46a6fYt6atS7dNcmZzj1dGAYFJ8GgIA1bfGAQjL1oydqLkkN7UVNVgZqGelS1tKGqu3fUv3iNl0omYIa3JxLDQpEcH4+kufOQtGQ5fKJjJ/wLzx5+aZxLZ2Mdvnrpv/HB7t1WVxiWR0zD7b/+HQLmL5vw8zAouOZrHpDu+z9g6IHH6cPI2ZOJ7GPHcLS8GtX68w86BgCdUoG5YcGYm5SI+cuWY+6V10DjH3BBz28Pr/vJeN2dPvgdtr66GR/tO4T63tHHakR4abB28UKsu+MuJC29xOVCg7Oeezo/BoVhysrKcOLECZSUlKC4uBhtbW1QKpV49913L+p4er0eH374IY4cOYL29nZ4e3tj/vz5WL9+PTQazUTLHYFBQRoX/IHBYoGssQZnDnyDoqyDyDtdhNz6ZrSNYwrS0SgEAdO1akT6+SAqNASRERGIjItHZHIKps+aA9UFXikYL3v4pXE+raVF+PCPv0HGyfwRt8sF4Pq5s7Huv56DJmT6RR+fQcE1X/PA1H3/jU11KP9mBwqzDyK/pBh5DS0wmM+/qJkMQJyPF+bNiMbcufOQdulqRKUvglyhmFA9jvC6n4iB/n6U7ngfez/7EN8UFp/zex3j7Ylrli/FNXfcjYQFi6a4SmnYw7lnUJAGg8Iwzz77LLKzs0fcdrFBoaurC48//jjq6uoQFBSE6Oho1NTUoLq6GsHBwfjDH/4AT0/b/dAzKEhjvB8YZO3NKPnwDWR8uRN7zzScd6ahIQKAUI0a0319EOzvj8DgEASGhsMvIgq+0fFYsXgOFG627Xs/HvbwS2O8Kvd9hX89/zRy6xpH3K5TKnDXuuuw5IH/gvwixi8wKLjmax6YvO+/ob4GZd/sQEHWAeSeLkJBSwfGs9ZxkNoNc8KnYU5KMuYsWYrUy66A1u/CrhaMhyO97ifK2NqIY+++hm+++RJZdef+WY/18cLKeXOwau01SF97HdzVzjndqnju29th6e5AV2UJOqvLYehoR7+hF/3GXpgMBvQZDOjvM8A8MACNpxc8vH2g9fGH2scPGv8guPsHQR4YCov8wkMrg4I0GBSG2bZtG/r6+hATE4OYmBjcc889Fx0U/vrXv2Lv3r1IT0/HQw89BLlcDgB44403kJGRgeXLl+OXv/zlREsWMShIY6wPDMKACYYDu7Dnw7fx6fFTqB1jQSABgwsVxYWHY1p4JIJj4hAQPxP+cSlw8/A45+OkeuN0tA8MFrMZx9/7P7z55quo+sHCUTO8PfHT+36FuLUbLuiYDAqu+ZoHbPf9N/XoUf3NZzi++ytk559CfmvHeQcdKwUBKYF+mJcQh7SFCzH3sjWYlpRsk3rOx9Fe97bSVXQC2R/+C98dPoDcxrZz7uehkGNxTCRWrrgEK9ffivD4BIfsoqRvbkRF7lFU5Oehovg0qqprUN/ahtr2DjR296Klr3/ck2WMxt9NiQR/X8RMn46IuESEzU6H/+wFEDRjXwFnUJAGg8IY1q9ff1FBob29Hffeey9kMhk2b94Mb29v8b7+/n784he/QHd3N1555ZUR900Eg4I0RvvAINN3ovj157H926+xt7511DdUlUyGlGkhSE6eiRkLliN66SpovC78Z4FB4cL09/Zg94tP4q2dO9H9gwGfK2dE49bH/wT/hNRxHYtBwTVf88AEvv8WC9pOZuFUxqfIyTmCrKoz0J9jYbAhOpUCc8KnIyl5JmLTl+HaG6+BWsvXvRTSwzxRlXUY2//5OrZ/tx95LWN/HwI93DEncjrmpKZiztLlSF11GTztZAV5fUc7Ko5mo+zkcZSfPo3yikpU1DegorUdjWP8UWuyuMkExPvqkJ4yC7PXrEPI8isg+8FVBwYFaThLUJhYx0sby83NhcViQVJSklUQUCqVmDdvHjIzM3Hs2DGsWLFCmiLJ9iwW9H61Ff/78gs42NRudbdCELAkKRGrbrkTiZdcAaVKJUGRrk2p9sDq3z6HhT/+OT76/aP47PgpMchlFpdh/5034/pFi3D1Y7+HNmjsKWKJxsNs6kfNtzuQ8+Vn2H/iBEo79WPu7+OmxJyIcMxMnYu4FVcgeM4iyL6/Ig1AspBAg8LnL8Cm+Quw0WJB5aH9+PajD5B58BAOVNVZjWlo7DHgy/xifJlfDLy3FTIAM/y8kRAehsjwcETNmIGIxGREps5BQHCITa8+9Pb0oPZ0AWoKCnCmrAQ1VZU4U1uHmsYmVLS2jWuBvfPxVinh5e4GN4UCKoUCKqUCKoUSKqUSgiCgu7cHXT296DQY0WE0wjBGKDaaLTjR3I4Tu/cCu/ci0P03WJKYgLmr1iDu6pug0vDnnibGroJCZWUlACAqKmrU+6OiopCZmYmKiooprIomk6KuEt/86TFszjpuNf5gurcOl1+5Fktuvwc+/kHnOAJNJc/wGNzx2ke4bPcXeOOF3yO3cXB6xD6zBf/efwA7rl+NW6+8CqseegJu/AVFF6i/qx3F299H1u5d+K7wNFqM557KWCUTkDZ9GubMnY+ky69FcGr6iGBA9kkQBEQuWoq7Fi3FXQB6W5tx+KP3kbnrS+zOK0DJD7o4AoAZQFFLO4pa2oHckyPu0yjkCPHSQqdWQ6dRw0ujhZenFl5eXuJ4xoGBAQwMDMA8MACz2YyBgQHo9Xq0tXegrasLbd16tPX0osNgPOfCcuOlkgkI99IiKtAfEdOmISo6GmZdENQBwfAKj4YuLAoqd/cLOmaf0Yiu1ibU5x5CzYmjKC8pQmlNDcrbO6wWzWw09OGT3BP4JPcE3F56HssTZsC4cROWrPsRZDLZhL42ck12FRSam5sBAL6+vqPe7+fnN2I/clxCfx9a3v1fPPfOOyjsGjlN4aKEBFx974NIXLicb2x2KmzFlfiv5Zcj97UX8H//flecGrHbNIDXtm/Htq934cc3345Fdz8IuVIpcbVkzwytjcj/5B0cyvwa+8oqrNbyGC7G2xPzk5Mxa/nliLn8Wqg8uKiXo1P7+mPFzzZhxc82wWI2o7ngJI5+vQu52UeQW1SM4w0tY65voTcNoKS1A8DUdetyGwoDAX6ICgtDZHQUor+/whEclwDFsFmybNHtTOXmBr+QMPiF3Ijkq24Ub+8z9KIu9zAKv/0ch7MO4Xhd04hB/EazGV/lF+GrX96PiMd/i1uvW4f19z+EwKDgi66FXI9dBQWDwQAAcDvHbDRDtxuN47/09/DDD1vdplKp8MwzzwAA/P39L7RMmxh6Ixnqw+ZKWo7sx8d334m3iypGjEPw06hx/2/+G0tvuG3SB7IFBPhN6vHPRaFQQC6XQ6ezj/62E7XysT9gyb2PIOMPj+KNjAx09g/+Qm/sNeKFN/+BmG1bcfsd92DJnb+E4vsuY1J97wFA1yVd8JTL5VAoFC75mh8y9P3vOlOFnH+/jr3f7MKBiiqrdTuGyAUgLSwUS5cuR9qN/4Gg5DkX/dxS/tw52+v+Ql3I9z4wKAhJKy7F7d+3TT165H/zJbIyv0VJ8WmUVdegorEZ5R1daO8//wJ5F0MGINjDDWE6L4QF+CE8NBThERGITkxCfFo6ImbPhWKcM75N6rnX6RCwZh1mrVmH9QA6aypw7MN/4sDur3GgtGLEeLLK9i786Z/v4Ll/vYurF6Thnkf/E5dd7dqL4U02Z/mcZ1dBYWhQMX9wnVflzm1Y++OfoOgHVxGuW3057vz9X+Cps80gdZo6Km9fXPvc67j0vjJ8/NTD2HLwiNiNrLStA0+99ByCXvsrbrr2R7j8/t8CkO4DG0mnKv8kMp57Ad/tzcSRmlqrLhNDNAo5ls6IwZJLr8Ds9XdCG2ibQXn7ykZfRXgqyNkl6qIpPDSYdc31mHXN9SNut/T3o7mkCCU5WairqkR7azPaW9vQ0dGB9o5OdHR1o6unB4IAyGVyyGUyyGQCZPLBbY1aDV8fb/j6+sEvMAB+QcHwDQ5BQNh0hCXPgpsD/i7yCovE8od+h+UP/Q6m7k4cf/vv2PHx+9hbXS/+Uc5kseDTQ1n49PobkBIehl/91xNYf9vtI66CEA1nV7MePf/88zhy5AjuuOMOXHXVVVb3Z2Vl4bnnnkN6ejp+9atf2aRWzno0dWp2fISbH/41yvVnrwhF+vvhnt89h8T5i6e0Fs56NHm6TuXgo+d+h88Kiq1mrtIo5Lh51XL89LdPIDw2bsprm5RZj8xmCOYBWAQBGPqHof/PcsVZj0pzsvHFO28iY+8+5Nafu8uot0qB5QnxSF99NeLWboDSw7YLa0pt6K/Jzvy6H4srz7oj+Xu+xYKQ8my89/eX8X7WMdQbrMf9hPn54p5f3odbbr8dHmNMK04XhrMeTYKhbkCtra2j3t/S0jJiP7p4Uz1NpHHne/jVs39CTe/Z6eNuWLsOGx57mrMYORnP5Hm445/bsfa7DOz4x1+xs7AU/d+Hcb1pAP/YlYk3d2ViVVIcrrrmWqy+7T/g62d/VxnM+k50lhSgpaQAzVVlaK6rRVdXBzq79eju6UFH7+Dgx46+fvSbLZAJgABhMCsAkAkCFIIAb3cVvNVq+Hpq8bafL3SenvDz80NIeASiklMQnjoXav8Ap7iS2tujx+FPP8a3n+/A7tzjKG0/9/tMgLsKlyQnYf4V1yL6yuuhUF3YAE8iGgdBwPTlq/Do8lV4qKEOma+8jLc+247M+rPrWtS0tOKJp57Ciy88jzvu/inu/ulPzzlWlFyPXQWFiIgIAEB5efmo9w/dPrQfOYbej17Hg3/5CxqHzWByx4/vwLpf/qeEVdFk81+2BncsW4P1BbnIfO1FbM06KvYpNgP4Ov80vs5/HvJnX8DimEhcedWVuOL2OxA8bdqU1Wgy9qI1/ziaCo6jrqQIZ6orcKapCbUdXajrMYxrdd/zaTT2Ax16oL4ZKK4YdZ9gtQoR3jpEBAUiKiIccckpmJE2HxFz5otjO+yRxWJB+ak87Nn6b2Tu3osDZZVjDkaerlVj+axZSFtzHaavWsuB7kRTSBkUgsuf/BNW//px5P9zM1596x18WtME0/eXftu79fjzX/6CV195Bf9xxx249957ERTEGQddnV11PWpra8PPf/5zyOVybN68ecTgn6EF17q6uvDqq69ywbUJmpIrChYLut/5Cx587XW09J2dcu7ee3+JNXfZbnXti8GuR1NvoKEGR157AR9kZqJCbxh1HwHArGlBSImNRUJKChIXLELC3DT4+Phc9PP26PWoKTiFqpMncPTYSTTUVKG2oQ5VzS2o7tKfs6+8PXCTCYjx8UJcaAhmxERjRlIyYuemIXL2PLip1VNeT09PD/L2fIuczG+QfTQXOeWVaD7PIlOJPl5Ys3A+rtxwC7rjF0NwwZnMXL3rkdSk7PpkD+/55/r6LZ3tOLPlDbz+7hZsqWywWsTQTaXCLbfeio0bN2LaFP4Bx1k4S9cjSYJCRkYGMjIykJ6ejltvvXXEfS+//DL27duHBQsW4MEHHxQHgb355pv44osvsGzZMtx33302q5VBYZJYLOj4vz/h/n+9i87vZ14QANz34KNYectdk/vc48CgIJ25Ogv2/+t1fLHzc2QUlqJpjLnyhwRrPRAfGgJvL094qN2hdldD7aGGh4cGHh5qGA0GdLS1o6OjHR2dXejo7ka7vgdN3Xo09IweSsYrUO2Oad46BPn7wVvnA61OB43OBxpfP2h8A6D2C4BSrYHFbIbFYobFbAbMg//3G3qhb25Ad1MDejvbYenpQENjE5rb21Hd0oqaTr3YLetCyAUg3FOLmOAAxEZEIDImBiGRUQiNjUdofAJ0Pj4T6sqk1+tRceoEyo4fQ1lBAcorylFcVYNT9U3nrVctl2FJxDSsXLIYK2/cgPB56WItUq6MLSUGBWkxKIz99Vtam9D24T/x9vsf4vXyBrT2jVxLQqFQ4MYbb8SmTZsQHR09maU6FQaFYY4ePYqPPvpIbBcXF0MQBMTGxoq33XDDDZg7dy4A4IMPPsDWrVtxySWXYNOmTSOO1dnZid/+9rdoaGhAUFAQYmJiUF1djerqagQFBeEPf/gDvLy8JlqyiEFhcnR//A/8/PkXxL9QyAUBD/7md1i6bsOkPu94MShIZ/j3fqCrA9kfbMEX2z9DxqkiVNtg1dMLJQMQ5umBSH9/hIVOQ0hEFPwiZ8BvRiK8o+OhcrfNX+5HG8zc39+P2pLTqDxxDBWFBagoLUVpVTWK6xtR1aW3Ggw+Xh4KOUI9NQjSeUHt7g61mxvU7u5wd3eDWq2Gu7safX1GdHfr0a0fHHPR3duLboMRjV3daND3nv9Jhpmh02DFzCSsuuJKpF+/Hu7eo18BYlBw3de9lBgUxvf1W+pr0P3W37Fl1zd4paxhRHdhYPAD6FVXXYVNmzYhNTV1Mkp1Ks4SFGwyRqGzsxPFxcUjbrNYLCNu6+zsHNexvLy88Kc//QkffPABsrKycOTIEeh0OqxZswbr16+HVssFduydrOwU/vtvfxVDglImw6+fehbzL18rcWVnSfWBRcp5/O2R3FOHBXf/Agvu/gWe6NHj1PaPkXtgPwpLSlBwpg6FbV3iFamJUAjANI0a07x1CA0IRFDINARFxyIwcTb8kmZDMUWz7Ax+YPjBz55nONyXhCNhybVIGHZzX3cHmk5ko+5kLmpKClFVXY2K5mZUjaO7VI9pACVtnShpG9/77oVwl8kwJ8gX8xLikLZwEeZdfiV8ZiQ4xWBsIlcmBIdB++s/4J7LrsHtW17FhyeKsLmsXpyExGKx4PPPP8fnn3+OJUuWYNOmTVi+fDlf+05uUroeORJeUbAtuUGPd+9ah/dLa8Tbfvv755G22n5CgpT4l8UL++ueubMddbnZKMg+grKSYuj1Peg1GtFjMKDXaESvsQ89fX1QyuXw1mrh7eUJnZcOOh9vePv6wScwCOFJyQhJngml1kvSv2jb6tybevRozz+K+vzjqCstQk11Jepb21Df2Y2GXoO4hsVEuckEhGvViPLzRVRoMKIiozBr0SIkrrocKt3FjRnhFQXXfd1LiVcULvzrtxiNsHzxIfq+2IptVY34W2k9yvTWV3yTk5OxceNGXH311VBycoIRnOWKAoMCg4LtWCw4/fv78diOr8Sb1l97HW757TO2fy4HxQ8M0v7SdoagMKY+Iww15WgvP4226gp0NtXDaOiF0WBAn8EAmbnv+4BlhEqhhNZDDU+NBhqNFp6eWmg8veAT4I+olFkISUmFwsu2i04xKLju615KDAoX//VbGmth/vfrMJ/Iwq6Gdmwua0Buu95qv6CgINxyyy247bbbEBoaOpFynQaDgpNgULCdvi/+jZ/+/r/RZRrscpQ4PQxPvbcTSqX9Tu841fiBgUHBUT8w2AKDguu+7l2VU5x7iwU++Ycxfec/odB34nBrNzZe2mZXAAAgAElEQVSX1SOzybpro0wux/ylK7Hm+psxK30xFobrRjmga3CWoGBX6yiQ41LWlOJ3f35ODAlalRL3v/QPhgQiIiJHJghoS16IrsgkTP/in1h46jAW+nmioLMHr5Y3YHtdO/rNg7/7zQMDOLznaxze8zVCwiLw41s34Nprr0VUVJTEXwRdLI6spAmT9Rnw6RP34/iwy5EbH30CwdO5MB4REZEzMGm8UH7j/Si96QH0a7yQ6OWBP6dG4cjKFDyWFIEQ/4AR+9fVVOLZZ5/F0qVLsWbNGvztb39DVVWVRNXTxeIVBZqwmv/9Hf7vVKnYvnLlpVh0zU0SVkRERESToT0pHV2RiZie8Rb88g7Az02JTZH++EWEBZ95L8dblU3IOfjdiG7deXl5yMvLwx//+EfMnj0bV111FZYtW4bk5GRxvSyyTwwKNCGWPdvxxLYdYjsyMAA/+d3zElZEREREk2nAwxMV129Ce3waIrf/H+TGXsgEAdd1VODy+HAc/Nm72HXwEHL37kJhYeGIxx47dgzHjh0DAHh7e2Px4sVYsmQJli5dipiYGE63amcYFOiiyYy9+Ouf/wct36/i6CaX48EXXoObjRaoIiIiIvvVnrwA+aFRiN76v9DUlgEAPBqqsGLbXxG99i7EPvkYTp8+je3bt+Ozzz5DSUnJyMe3t2Pnzp3YuXMngMHZk1JTU5GSkiL+Cw0NZXiQEIMCXbS6LX/DV7XNYvvuTQ8hIi5RwoqIiIhoKvX5BKLozicw7et/I+hwBgBA3m9E1CebUV9ZiPar7sCyDfdg6fqfoaqsGAe+ycDxIwdQnJ8H88DIBTUbGhqwa9cu7Nq1S7xN66VDVFwiwqNjERwWjuBp4QgOC0dQaBiUqrEnTJF6pjdnwKBAF0Xe2YrX3n9fbCdPD8Nlt94tYUVEREQkBYtCiZo1P0ZXZCIiP30VCkMPACDgaCZUHc0ou+l+mN08EBETh4iYONxyz/3o1XfjVG428nIO4UTWIVQUF4567O7ODuRlH0Je9qERtwuCAP+gEASFhsE3MAi+/oHwDQiCb0Ag/AIGt3v9IqBWs5fDRDAo0EU5/dqzONp2dk70W3/9JC8NEhERubCOhDQUBEcg6qO/Qlsz2M1IV5qH+DefRsltj6Lf8+yq7mqNFmlLVyBt6QoAQGd7G4ryclFWVIDy0wUoO52Pprracz6XxWJBU30tmurPvQ8AeHh4wM/PD35+fvD19RX/1+l00Ol08PHxgU6ng7e3N3Q6HTw9PaHVauHm5sbPNWBQoIugrK/C5p0ZYntxcjKSFiyTsCIiIiKyB33eATj9k98ictsr8D11GMDguIWE159E8W2PwhAYNurjvLx9MH/ZKsxftkq8rbuzA+WnC1BeXIi66krU11ShrqYKTfW1Vt2WzqWnpwc9PT2orq6+oK9DqVRCo9GIwcHDwwPu7u6j/pPJZJDL5ZDL5ZDJZFAoFPD09IRcLsc111yDoKCgC3pue8KgQBfsyP8+jVK9AQAgE4D1//l7iSsiovFw1ZWRiWhqWRQqlN/wS/R5+SP44OcAAFVnC+LfeAqlGx5Ed1TyuI6j9dJhZtpCzExbOOJ2k6kfTfV1qK+pRFNdLVqbG9HS1Ii2pka0NDWgtbkRnW2tE/oa+vv70d7ejvb29gkdJz09nUGBXIfs9HG88t3ZfoKXL1nGAcxEREQ0kiDDmctvRZ+3P6Z/8RYEWKAw9mDGO/+DiuvuRdvMJRd9aIVCiZCwcISEhZ9zn/7+PsSo+9HS0oKWlha0traiubkZLS0tYgDo6Oiw+n/4+g+2IJM59trGDAo0fhYLvn3592g09gMYnA71hl8/JXFRREREZK+a0i9Hv6cPoj7+G2SmfsjMA4j++O+oMPWjZc6KSXtepVKF4GA/BAcHj/sxFosFPT096OrqQnd3N7q7u8Xt3t5eGAwGq39GoxFmsxkDAwMwmUzitlKphMlkgo+Pz/mf2I4xKNC4CVmZeCP3lNhed/W18A8OlbAiIiIisnftifNx+ie/Rcx7L0DZM9gFMmL76zB5eKIjfp7E1Z0lCAI0Gg00Gs2EjxUQEAAAaGpqmvCxpOTY10No6gyYsO1vz6HLZAYAeLm54Zr7/lPiooiIiMgR6MNmoOiu36FP6w0AECyWwYXaqookrozGwqBA49K3ayveLaoQ2zfd+hNovXTSFUREREQOxegXjJLbH8OA2+DaBjJTP2Lfex7ujTUSV0bnwqBA5yX0G7Hljc3o/36AT7CXJy6/c5PEVREREZGj6Q0KR8nNj8AsVwIAFIYezHjnGSg7miWujEbDoEDn1b3rI3xa1SC2b7n3Aajc3CSsiIiIiBxVd2Qiym/YBMv3C5qputoQ9/YzkPdwCmd7w8HMNDaLBRkfviM2owP8sORHt0hYEBERETm69sT5qLrqTkR8/gYAwL2lDrFbnkPxf/w/mFXuNnkOKdeO0XXJsDTaT7LntxVeUaAxyU4fx6fFFWL76pt+DLlcLl1BRERE5BSa0y5F7YobxLb2TCmiPv47YOO1DOjiMSjQmLLeeUWc6UirUmLRTT+WuCIiIiJyFnXLf4TGtMvEtndRDvyPZkpYEQ3HrkcS2lfWgo4O++2PJ+9qx9aDh8X26uUroPaY+NzCRERERAAAQUD1lT+BW1sjdKUnAABhX76Dzqhk9PkGSVwc8YoCndOZj95AYVev2F51130SVkNEREROSSZD5bp7YFJrAQDyfiOitm0GzGaJCyMGBRqdeQCff/6p2EyLjUZYTJyEBREREZGz6vf0QeXVd4ltbXUxgvdvl7AiAtj1iM7BfPhbfFndKLZX3/4zCashZyLlLBRERGS/2pMXoKVoCfzy9gMAQnZ/hI7YVPSGREpbmAvjFQUaVeaWf4gLrAVqNZi3+hqJKyIiIiJnV33VT9Dn5QsAkJkHEPXJ3yGY+iSuynUxKJAVZX0Vtp44JbbXXL0OcgUvPhEREdHkGnDXoGLdvWJb3XQG0775QMKKXBuDAlnJf/cV1Bv6AQBKmYAV//FziSsiIiIiV9EVnYKGBWvEduChDGjLT43xCJosDAo0gszYi23ffiO2l8+dBx//QAkrIiIiIldz5tIN6PUPBQAIsCBq26uQGXskrsr1MCjQCJ27PsKB5g6xfekdGyWshoiIiFyRRalC+fUbYZHJAQCqzhYEHdgpcVWuh0GBzrJYsGvru2IzNigQCWmLJCyIiIiIXFVvSBTql5ydTCXo4E4o9B1jPIJsjUGBRPKiXGwrqRTba266DYIgSFgRERERubL6JVePWIgteB/XVphKDAokOvze6+gyDa6CqFUpsfiG2yWuiIiIiFyZ2c0D9UuvFdsBWV9D2dEiYUWuhUGBAACCqR+7srLF9uply6H20EhYERERERHQOH81+jx9AACygX6E7P1E4opcB4MCAQD6j+5Fdkun2F6y4U4JqyEiIiIaZFGqULf8OrHtn7sHbi31ElbkOhgUCACQ/dmHsHy/HerliehZ8ySth4iIiGhIy5wVMPoMTtcuWMwI3b1V4opcA4MCQTD149ujR8X2siXLOIiZiIiI7IZFrkDtihvEts/JQ1A3VElYkWtgUCD0ZY/sdpR2/W0SVkNERERkrTVlMXoDwgAMLsIWmvmhxBU5PwYFQvb2D852O9J5ImbmXEnrISIiIrIik+HMqpvEpnfRUWhqiiUsyPkxKLg4wdSPzNxcsb1sMbsdERERkX3qiJ8HfWi02A79llcVJhODgotjtyMiIiJyGIKAM5duEJte5afgWXZSwoKcG4OCi2O3IyIiInIkXVHJ6IxMEtsh322TsBrnxqDgwgRTH75ltyMiIiJyJIKA2pU3ik3PigK4NddKWJDzYlBwYcasPcgZ1u1o/g23S1gNERER0fjop8ehJ3C62PY/milhNc6LQcGFZe/4cES3o+iUOZLWQ0RERDQugoDmeavEpv+xvRBMfRIW5JwYFFyUYOpDZu4xsb188XJ2OyIiIiKH0TprCcwKFQBA0dsN74IsiStyPgwKLuqH3Y7SbuBsR0REROQ4Btw1aE1ZJLYDctj9yNYYFFxU9vaz3Y6m6bzY7YiIiIgcTtOw7keelRzUbGsMCi5IMPUh89jZbkec7YiIiIgcUc+0GPQEhYvtgJxvJazG+TAouCDjkd0jF1ljtyMiIiJyRD8Y1Ox3/DsOarYhBgUXlL1jq7jNbkdERETkyFpmLsGA0g3A4KBmn3wOarYVBgUXIwyYkMlF1oiIiMhJmN090DZsULN/zjcSVuNcGBRcjOlUDnJau8Q2ux0RERGRoxsxqLmqCO5NZySsxnkwKLiY/K+3i7MdBWo82O2IiIiIHF5PaDR6giPEtj8HNdsEg4KLyck5Im6nJSez2xERERE5PkEYcVXB7/h3EPo5qHmiGBRciLy7HQeqz84vPPOSyyWshoiIiMh2WmcuPjuo2aCHT8GR8zyCzodBwYU07clAk9EEABAAJF12tbQFEREREdmI2c1jxErN7H40cQwKLuTE3q/E7aSQIHh5+0pYDREREZFtNaddKm57VhXBraVewmocH4OCq7CYcSQ/X2zOnZcuYTFEREREttcTEoWewOli27uQaypMBIOCixAqTyOnuV1sJ19+rYTVEBEREU0CQUB74nyx6V3AoDARClserK+vD9u2bcP+/fvR3NwMrVaL1NRUbNiwAX5+fhd0rGPHjmHnzp0oLS1FT08PNBoNYmNjcfXVV2PmzJm2LNslFO/6DH3mwYlRtUolYuYtlLgiIiIiIttrT0hD6J6PAQDaM6VQdrai34vdrS+Gza4o9PX14emnn8bWrVthMBiQlpYGPz8/7N69G4899hjq68ffR2zHjh344x//iOPHj2PatGlYsGABAgICcPToUTz99NPYtWuXrcp2GblH9ovb8+JmQKFQSlgNERER0eToDQqH0SdQbHsX5UhYjWOz2RWFTz75BEVFRYiLi8Pjjz8Od3d3AIMf+t966y1s3rwZTz311HmP09nZiS1btkChUOCJJ55AQkKCeN+hQ4fw0ksv4e2338by5cvF56CxyfoMOFRWKbZTl1wiYTVEREREk0gQ0JaQhuCDOwEMdj9qmr9a4qIck02uKJhMJmRkZAAA7r777hEf4NeuXYuIiAgUFBSgrKzsvMcqLi6GyWRCSkrKiJAAAAsXLkR4eDiMRiNqampsUbpL0GftQZneILaT1/xIwmqIiIiIJld7Qpq47VlRAHlvt4TVOC6bBIXCwkLo9XoEBQUhKirK6v4FCxYAALKzs897LKVyfF1itFrthRXpwvK++ULcjvDRIXBauITVEBEREU0u/fQZ6Nd6AwAEixm607kSV+SYbBIUKisHu7WMFhIAIDo6esR+Y4mJiYGHhwdOnjyJwsLCEfcdPnwYVVVViI+PR3Bw8ASrdh3ZJ46J22mpcySshIiIiGgKCDK0x88Vmz6c/eii2GSMQnNzMwCcc2YjX1/fEfuNRaPR4Oc//zlefvllPPnkk4iPj4evry8aGxtRWlqK2bNnY+PGjbYo2yXImutwqK5RbKesulLCaoiIiIimRnvCfAR8vzqzV+kJyPoMMKs4vvVC2CQoGAyD/d/d3NxGvX9ozMLQfuezcOFCaLVavPTSSyOuKuh0OiQnJ8PT03PctT388MNWt6lUKjzzzDMAAH9//3Efy5YUCgXkcjl0Ot2kPk/ZJ/9El8kMAFDKZEhfeyPUHh6T+px0bnK5HAAm/byT/eG5d108966L515iqQsx8JEGcoMeMlM/gutKoZ+1eEqeWi6XQ6FQICAgYEqeb7LYJChYLJYJ3f9D27dvxzvvvIP58+dj/fr1CAwMRGNjI95//3288847KC4uxiOPPDKRkl1G9p6vxe25UREMCUREROQaFEr0JM2H59HdAADtiQNTFhSchU2CglqtBgAYjcZR7x+6fTzTmebn5+Ptt99GVFQUHn74Ychkg8MowsPD8cgjj+A3v/kNDh8+jOPHjyM1NfW8x3vxxRfHvL+5ufmCg4wtBAQEYGBgAB0dHZP3JAMmHBh2RSZ1/uLJfT46r6G/KvE8uB6ee9fFc++6eO6lJ8SkikFBfeoIOltbYJHbdL3hUel0OphMJjQ1NU36cw0nCAJCQkJsdjybDGYe6r7T0tIy6v2tra0j9hvLnj17AAzOlDQUEobIZDKkp6cDAE6dOnXR9boKU0Eu8trOTgeWfMU6CashIiIimlqdsbNg/n6RWYWxB57l+RJX5FhsEhQiIiIAAOXl5aPeP7R+wtB+YxkKFUNXKX5o6Pbubs6Hez6FX32GoWslAR5qTE+cKWk9RERERFPJrHJHZ8wsse1dyNmPLoRNgkJCQgI8PDzQ0NAwalg4fPgwAGDu3LlW9/3Q0GW60tLSUe8fuj0wMHDU++msnJwj4vb85GQIgiBhNURERERTr23Y4mvehTmA2SxhNY7FJkFBoVBgzZo1AIA33nhjxOxGO3bsQGVlJRISEhAbGyvenpGRgQcffBBbtmwZcayhrkX79u2zWqAtKysL+/btgyAI4n40Opm+EweqasX2zOVcupyIiIhcT0f8XFiEwY+8Sn0HNDXFElfkOGw2muP6669HXl4eioqK8MADDyAhIQHNzc0oLi6Gp6en1doHnZ2dqK2tRVtb24jb58+fj4ULF+LQoUN49tlnERMTg4CAADQ1NYlXE26++WaEhobaqnSn1Lw3A43GfgCAACBp9TXSFkREREQkgQG1Fl2RifAqHxzf6l2YDX14vMRVOQabBQWVSoUnn3wSn3zyCfbt24esrCxoNBpccskl2LBhw7jXKxAEAQ899BAyMzOxZ88eVFVVoaKiAh4eHpgzZw6uvPJKzJ4921ZlO62TB3aL24nBgfDy8ZWsFiIiIiIptSekiUHBpzAbZ1bfCrBL9nnZdH4olUqFDRs2YMOGDefdd/369Vi/fv2o9wmCgFWrVmHVqlW2LM+l5BWeHdU/e+b5p5ElIiIiclbtCWkI/+JfAAC3tkaoG6rQG3z+SXZcnU3GKJB9kfV04Wh9s9iesXiFdMUQERERSazfyxf6aTFi27swe4y9aQiDghNqPbQbrX0mAIPjE2KXXCptQUREREQSa4+fJ257lZ2UsBLHwaDghIoO7ha3Z/j7Qqvzlq4YIiIiIjvQGZ0ibmvOlELWZxhjbwIYFJzSyfyzKTklMUnCSoiIiIjsQ09IFExuHgAAwTwAbWWhxBXZPwYFJyP0G5F7pk5sJyxYJmE1RERERHZCJkNX1Nk/oLL70fkxKDiZ7twDONPbJ7ZjLrlCwmqIiIiI7EdXVLK47fn9dKl0bgwKTqZ437fi9nSdJ3wDgyWshoiIiMh+DA8KHg1VUOg7JazG/jEoOJn8vGPi9qy4OAkrISIiIrIvBv9Q9GnPTvLCqwpjY1BwJgMmHK2qEZsJ8xZJWAwRERGRnREEdEWz+9F4MSg4kb6iEyjt7hXbsSvWSFgNERERkf0Z3v3Ii0FhTAwKTqR07y5xO8BDjaDImDH2JiIiInI9nVFn11Nwa2uEqr1JwmrsG4OCE8k/dnY58tSYaAiCIGE1RERERPanX+cHg1+I2PYs41WFc2FQcBYWM46VV4jNxLnp0tVCREREZMc6R3Q/4noK58Kg4CTMVSXIb+8W27FcP4GIiIhoVFbrKVgsElZjvxgUnET5ngyYv9/2UikRljhL0nqIiIiI7FVXVBIsGOyirdR3wr2xWuKK7BODgpMoPHpE3E6NjIBMxlNLRERENJoBtRY9IRFim7MfjY6fJp3E8ZIScTsxda6ElRARERHZP6vuR2SFQcEJCI1ncKKlQ2zPWLZawmqIiIiI7F/XsGlSPSsKAPOAhNXYJwYFJ1D93S4YzYODcNzlMkTMXSBxRURERET2rTs8Dma5AgAg7zNAc6ZM4orsD4OCEyg8ckDcTgmbBqVSJWE1RERERPbPrHKHPmyG2PbkNKlWGBScwInTheJ2SkqqhJUQEREROY6R6ylwnMIPMSg4OKGzDbmNLWI7bslKCashIiIichxd0WeDgqa6GEK/UcJq7A+DgoNr2P81uk2DKygoBAFRi1ZIWxARERGRg9CHRmNA5Q4AkA2YoK0qkrgi+8Kg4OCKDn8nbicGB8LdQyNhNUREREQORK5AV0SC2PQqY/ej4RgUHNzJgrM/0ClJyWPsSUREREQ/NGKaVI5TGIFBwZH1G3GsrlFsxi1YJmExRERERI5n+DgFj7oKyIw90hVjZxgUHJj+VA4ajf1iO3LxKgmrISIiInI8vYFhMLl5AAAEWKCpKZW4IvvBoODAKg7tFbeneWrgHRAkYTVEREREDkiQQT89Vmxqq09LWIx9YVBwYKdPHhO3kyIiJKyEiIiIyHF1T48TtzU1xRJWYl8YFBxYfkWluB3HhdaIiIiILop+WFDQ1pQAZrOE1dgPBgUHZelsw8mWdrEdmc6BzEREREQXQz8tBhZBAADIjb1wb6qRuCL7wKDgoBoP74bBbAEAKAUB0+cukrgiIiIiIsdkVrmjN+hsN25tNbsfAQwKDqss+6C4HRfoBze1WsJqiIiIiBxb9/QZ4raW4xQAMCg4rKKifHE7IWbGGHsSERER0fkMDwoaznwEgEHBMVksOFVTKzZjZ6dJWAwRERGR4xs+oNm9tQEKfYeE1dgHBgUH1F9VgpKus6sGRnChNSIiIqIJ6dP5o0/rLbY1HKfAoOCIqg5+K27rVEoExyZIWA0RERGRExAE6MOHTZPKoMCg4IhKjuWI20nTQiF8P50XEREREV287rDhQYHjFBgUHFBhydmEG5+QJGElRERERM5j+IBmj9pyCKZ+CauRHoOCo+k3Iq+hSWxGz1soYTFEREREzqM3JBJmuRIAIBvoh0d9hbQFSYxBwcF0njiClj6T2OZAZiIiIiLbsMgV0E+LFtuuPqCZQcHBVBz+TtwO99LCyy9AwmqIiIiInIs+bNjCay4+ToFBwcGcPnVC3E6KjJSuECIiIiInNGKF5upiwGKRsBppMSg4mPyKCnE7LiVVukKIiIiInJB+WFBQdrdD1dEsYTXSYlBwIOaOVuS3nl0lMHLBcgmrISIiInI+Jo0OBt8gsa1x4e5HDAoOpP5gJozmwctfKpmAsDmc8YiIiIjI1rqnc+E1gEHBoZTlHBS34wP9oXRzk7AaIiIiIuekZ1AAwKDgUIqKCsTtxNgZY+xJRERERBdr+IBmdUMlZH0GCauRDoOCo7BYcPJMrdiMSU2TsBgiIiIi52UImIYBNzUAQLBYoKkpkbgiaTAoOIi+ymJUdPeK7cglXGiNiIiIaFIIMnQPW0/BVRdeY1BwEJUHvhW3fdxUCIiOl7AaIiIiIuc2fJpUbY1rznzEoOAgSo5li9vJYaEQBEHCaoiIiIic2/CZjzTVJYDFLGE10mBQcBAFpWf7xiUkJEtYCREREZHz00+LgeX7P8wqjD1wb6o9zyOcD4OCI+g3Iq+hSWxGpy2SsBgiIiIi52d2U6M3KFxsa2pcb5wCg4IDaMs9hPb+AbEdvpgDmYmIiIgmmz40RtzWnCmTsBJpMCg4gMqcA+J2pM4TGm9fCashIiIicg36adHitqaWQYHsUFnBSXE7IXy6hJUQERERuY6e0LNBQd1YDcHUJ2E1U49BwQEUVlaK2zHxSRJWQkREROQ6egPDYFYoAQCCeQDq+iqJK5paDAp2zmI0oKC5VWxHzl0oYTVERERELkQmR09wpNh0te5HDAp2rvXYIXSbBuftFQBMS18ubUFERERELkQf6rrjFBS2OlBfXx+2bduG/fv3o7m5GVqtFqmpqdiwYQP8/Pwu+Hj19fXYtm0b8vLy0N7eDnd3d4SEhCA9PR3XXnutrcq2e1U5B8XtKG9PqD29JKyGiIiIyLX0DBvQ7OFiQcEmVxT6+vrw9NNPY+vWrTAYDEhLS4Ofnx92796Nxx57DPX19Rd0vCNHjuCRRx5BZmYmtFot5s+fj6ioKDQ2NuKrr76yRckOo7QgT9yOnx4+xp5EREREZGvDryi4N9VCZuyVsJqpZZMrCp988gmKiooQFxeHxx9/HO7u7gCAHTt24K233sLmzZvx1FNPjetYFRUV+POf/wy1Wo1f//rXSEhIEO8zm80oLy+3RckOo2j4QOYEDmQmIiIimkpGv2AMuKkhN/ZCgAUedRXojkyUuqwpMeErCiaTCRkZGQCAu+++WwwJALB27VpERESgoKAAZWXju1Tz5ptvwmQyYePGjSNCAgDIZDLExMSc45HOx9JnQEFzm9jmQGYiIiKiKSbIoA+JEpuuNE5hwkGhsLAQer0eQUFBiIqKsrp/wYIFAIDs7OzzHqumpgYFBQUICQnBvHnzJlqaw2vJPQT9wNmBzKHpy6QtiIiIiMgFDV9PwZXGKUy461Hl911jRgsJABAdHT1iv7GcPDm4sNisWbPQ19eHAwcOiFciIiIisGjRInh4eEy0ZIdRPWxF5ihvL6i1HMhMRERENNVcdYXmCQeF5uZmADjnzEa+vr4j9htLdXU1AEClUuHRRx9FbW3tiPu3bNmCRx55BElJrtFXvySfKzITERERSW34FQW3tkbIe7ow4OEpYUVTY8JBwWAwAADc3NxGvX9ozMLQfmPR6/UAgJ07d0Kj0eBXv/oVUlJS0N7ejq1bt2Lfvn147rnn8OKLL8LHx2dc9T388MNWt6lUKjzzzDMAAH9//3Edx9YUCgXkcjl0Ot059zldfXb1v4SZs8fclxyDXC4HAJ5LF8Rz77p47l0Xz70T8fLCgMYLcn0nACCwowE9IWHn3F0ul0OhUCAgIGCqKpwUEx6jYLFYJnT/cGbzYH/8gYEB3HfffUhPT4eHhwdCQ0Nx//33IyYmBnq9Hl9++eWEanYEZkMv8pvOrsgcu4DjE4iIiIgkIQgwTJ8hNt2qSyQsZupM+IqCWq0GABiNxlHvH7p9+GxI5zK0j6+vL1JTU63uX7lyJQO6Oy4AACAASURBVEpLS3Hq1Klx1/fiiy+OeX9zc/MFhRlbCQgIwMDAADo6Oka9v/lwJnq+H8gsA+CTMu+c+5LjGPqrEs+l6+G5d108966L5965eASGQ1OYAwCQl+WPeV51Oh1MJhOampqmqjwAgCAICAkJsdnxJnxFYajrTktLy6j3t7a2jthvLIGBgQBwzss0Q7d3dnZecJ2OpiL77IrM0T5ecNc4fz84IiIiInvlijMfTTgoREREAMA5F0IbPmvR+URGRgIAuru7R72/q6sLwPiuTji6soKzA5m5IjMRERGRtIbPfKTqaoOyq22MvZ3DhINCQkICPDw80NDQMGpYOHz4MABg7ty55z3WzJkz4ebmhvr6+lFnScrPzwdw7qlYnUlRZYW4HZuYLF0hRERERAST1ht9Xr5i2xWuKkw4KCgUCqxZswYA8MYbb4yY3WjHjh2orKxEQkICYmNjxdszMjLw4IMPYsuWLSOO5ebmhiuvvBIDAwN4/fXXRxzr2LFj2LNnDwRBwGWXXTbRsu2a2diLgpZ2sR3BFZmJiIiIJKcf1v1Ic8b5g8KEBzMDwPXXX4+8vDwUFRXhgQceQEJCApqbm1FcXAxPT09s3LhxxP6dnZ2ora1FW5v1JZsbb7wRBQUFOHr0KB544AHExsais7MTp0+fhsViwc033zwidDij5qMH0TtsIDNXZCYiIiKSXk9oNHwKswG4xhUFmwQFlUqFJ598Ep988gn27duHrKwsaDQaXHLJ/2/v3oOjqu//j7+W3dwTEsIlCeVrEhIxXxsMctXCD+zUn4Iy/XEpRJ3x0kE7irZctNUWR6RI0TqDbaczjK2iCEWLKIwyNPprFUvwZ7gr5mbIjViSkA0kwcUkbLK/P/zumvUk5HJO2E3yfPy1e87hs+/1zeJ57dnP58xRdnZ2r+5V4B3rnXfe0YEDB3TixAmFhIQoMzNTt99+e49+wjTQnT767UTmtBGxCouMDmA1AAAAkL5zReFMmeTxSDZbACvqX5YEBembE/zs7GxlZ2d3e+zSpUu1dOnSrotyOLRo0SItWrTIqvIGlNKOE5m5IzMAAEBQ6LjykePrrxTaUKfWEWMCWFH/Mj1HAdYrPl3pe5yWwURmAACAYNAWEaXm+ATf86hB/vMjgkKQaW++qCImMgMAAAQlv/spDPIJzQSFINNxIrPdJo2dNivAFQEAAMDLb55CNUEBV1Dldycyc0dmAACAoOF/h+ZyydMewGr6F0EhyJQWdZzIzB2ZAQAAgsnFpBR5/melI3trs8Kd1QGuqP8QFIJM8enTvsdMZAYAAAgu7aHhah71Pd/zwXw/BYJCEGlvvqgi57cTmVMm3xjAagAAANAZ1/e+cz+FQYqgEETqjv0/Nbd/O5E5cdrMAFcEAACA7/Kfp0BQwBVQcfRj3+P0EXHckRkAACAIdVz5KLLmtNTeFsBq+g9BIYiUFeb7HnNHZgAAgOD09Zj/ksf2zWn0MHfroJ3QTFAIIl9UMZEZAAAg2HlCQvX16A4TmqvLA1hN/yEoBIn21ha/iczckRkAACB4XRyb6nscWV0RsDr6E0EhSDhPfOI3kTmJOzIDAAAErYuJKb7HXFFAvzrd4Y7M4+OGK5SJzAAAAEHrYlKK73FkTeWgvEMzQSFIlBV+e0fmCePGBbASAAAAdOfrxGR59O0dmsPqawJckfUICkGi5HSl73HahP8OYCUAAADozjd3aE7yPR+MPz8iKASBdrdbRXXnfM+TJ00LYDUAAADoiYtJg3tCM0EhCJz//Ki+cn9zow6bpLEzZge2IAAAAHSLoIB+V9Xhjswpw6MVHhsfwGoAAADQE34TmqsrJI8nUKX0C4JCECgv+Mz3+JpxYwNYCQAAAHrqYmKy77Gj5aJCz58NYDXWIygEgS8qvp38kpaeEcBKAAAA0FPt4ZFqjk/0PR9sPz8iKASYp71dhbVO3/PkrCkBrAYAAAC90fHnR1GDbOUjgkKANRSfVOMlt+/592bMCWA1AAAA6A3/Cc0EBVio6nCu7/G46AhFjk68zNEAAAAIJoN5QjNBIcAq8j/1Pc4Yy0RmAACAgaRjUHB8/ZVCG51dHzzAEBQCrKS81Pc4Pf3qAFYCAACA3mqLiFZL3Gjf88E0oZmgEECe9nYV1ny7jFbyxMkBrAYAAAB94f/zo8EzT4GgEEDnTxWqvuWS7/m46dyRGQAAYKAZrHdoJigEUPnHH/geJ0aEKXpcSuCKAQAAQJ8YVj4aJBOaCQoBVHriiO9xRhKrHQEAAAxEF5O+vUNziKtJ9sb6AFZjHYJCABWf+sL3+Oq09ABWAgAAgL5yR8WqdXi873nYf0ovc/TAQVAIoMIzNb7HKZnXB7ASAAAAmNHx50fhVacCWIl1CAoBUnPqC9VcbPY9Hzf9fwWwGgAAAJjRceUjrijAlKP/N8f3eFRYqGJSuIcCAADAQOXqcEUh7EuCAkw4+vFB3+OMxNGyDaMVAAAAA5XfHZob69XWcC5wxViEs9MAOfF5vu/x1alpAawEAAAAZrljRuhSdNy3z0uLA1iNNQgKAfLZ6S99j1O+nxXASgAAAGAFV4erCpfKCArog/PVZ3T6wkXf86uYyAwAADDgdVz56FJpUQArsQZBIQDy9397R+a4EIdiJ2QGsBoAAABYoeM8hUv89Ah98dnhT3yP/zthFBOZAQAABoGOVxTa62rk+aopgNWYxxlqAHxeUOh7fHVK6mWOBAAAwEBxaXi8LkXGfLvh9MBeJtUR6AKGopMdJjKnXntdACsBAACAZWw2XUxKUeS5WkVfc60uhYYHuiJTCApXWJPTqbLGr3zP/2vKDwJYDQAAAKxUeudjGh4/UrPGj1RdXV2gyzGFnx5dYaePH1GE/Zv/7DEOu0ZkTglwRQAAALCKxz54vocfPO9kgMj833NVlP+5zp78VB8f+VTDHLQAAAAAwYez1ABwxMQq6/8s1oWJN6mxsTHQ5QAAAAAG/PQIAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAgcOqgVpbW7Vnzx4dPHhQTqdT0dHRysrKUnZ2tkaOHNnncaurq/XYY4/p0qVLysrK0po1a6wqGQAAAEAXLLmi0NraqvXr12vXrl1qbm7W1KlTNXLkSO3fv1+PP/64ampq+jz2X/7yF7ndbivKBAAAANBDlgSF3bt3q7i4WBMmTNAf//hHrVq1Sr/73e90zz33qKmpSZs3b+7TuB988IHy8/P1ox/9yIoyAQAAAPSQ6aDgdruVk5MjSVq2bJnCw8N9++bPn6/k5GQVFhaqrKysV+M2NjZq27ZtmjhxombOnGm2TAAAAAC9YDooFBUVyeVyKSEhQampqYb9M2bMkCQdOXKkV+O+8soram1t1QMPPGC2RAAAAAC9ZDooVFZWSlKnIUGSxo8f73dcTxw7dkwff/yxFi5cqMTERLMlAgAAAOgl00HB6XRKUpcrG8XHx/sd153m5ma9/PLLGjt2rBYsWGC2PAAAAAB9YHp51ObmZklSWFhYp/u9cxa8x3XnjTfeUF1dnZ566ik5HOZXb129erVhW2hoqJ599llJ0qhRo0y/Rl84HA7Z7XbFxsYG5PURGHa7XZLo+xBE74cuej900fuhy263y+FwaPTo0YEuxRTTVxQ8Ho+p/R2VlpYqJydHs2fPVmZmptnSAAAAAPSR6a/sIyIiJEktLS2d7vdu77gaUmfa2tr04osvKjIyUvfcc4/Zsnw2bdp02f1Op7NXYcYqo0ePVltbmxobG6/4ayNwvN8q0fehh94PXfR+6KL3Q1dsbKzcbrfq6uqu6OvabDYlJSVZNp7poOD96U59fX2n+8+dO+d3XFfq6+tVUVGhuLg4w8m9y+WSJJ06dUpPP/20wsPD9cQTT5gtHQAAAEAXTAeF5ORkSVJ5eXmn+733T/Ae152GhgY1NDR0us/lcqmgoECRkZF9qBQAAABAT5kOChkZGYqMjFRtba3Ky8sNy6Tm5eVJkiZPnnzZccaMGaOdO3d2ui8/P1/r1q1TVlaW1qxZY7ZkAAAAAN0wPZnZ4XBo7ty5kqQtW7b4rW60d+9eVVZWKiMjQ+np6b7tOTk5WrlypXbs2GH25QEAAAD0A/Prj0patGiRTp48qeLiYq1YsUIZGRlyOp0qKSlRTEyMli9f7nd8U1OTzpw5o/Pnz1vx8gAAAAAsZvqKgvTNfQnWrl2rxYsXKzQ0VIcPH9bZs2c1Z84cPffcc9xdGQAAABhgbJ5ArA0aRKqrqwO2PGpuWT1Lpg0xLJU3dNH7oYveD130fuiKjY3VrPEjB/zyqJZcUQAAAAAwuBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYEBQAAAAAGBAUAAAAABgQFAAAAAAYOKwaqLW1VXv27NHBgwfldDoVHR2trKwsZWdna+TIkT0aw+Vy6fjx4zp69KgqKirkdDpls9k0btw4zZo1S7fccoscDstKBgAAANAFS866W1tbtX79ehUXF2vEiBGaOnWq6urqtH//fh07dkzPPPOMEhMTux3n3Xff1dtvvy2bzabU1FRNmTJFTU1NKi4u1qlTp/TJJ59ozZo1CgsLs6JsAAAAAF2wJCjs3r1bxcXFmjBhgp588kmFh4dLkvbu3avXXntNmzdv1rp167odJzw8XAsXLtStt96q+Ph43/bq6mqtX79eRUVFeuutt3TXXXdZUTYAAACALpieo+B2u5WTkyNJWrZsmS8kSNL8+fOVnJyswsJClZWVdTvWggULdOedd/qFBElKSkryhYODBw+aLRkAAABAN0wHhaKiIrlcLiUkJCg1NdWwf8aMGZKkI0eOmHqdlJQUSdL58+dNjQMAAACge6aDQmVlpSR1GhIkafz48X7H9VVtba0kKS4uztQ4AAAAALpnOig4nU5J6nJlI+/PiLzH9dW+ffskSVOnTjU1DgAAAIDumZ7M3NzcLEldrkTknbPgPa4v3n//fZ08eVJRUVFasGBBr/7s6tWrDdtCQ0P17LPPSpJGjRrV57rMcDgcstvtio2NDcjrIzDsdrsk0fchiN4PXfR+6KL3Q5fdbpfD4dDo0aMDXYoppq8oeDweU/u7U1BQoFdffVU2m00PPfSQYaIzAAAAAOuZvqIQEREhSWppael0v3d7x9WQeqqyslLPP/+83G63fvrTn2r69Om9HmPTpk2X3e90Ok2Hmb4YPXq02tra1NjYeMVfG4Hj/VaJvg899H7oovdDF70fumJjY+V2u1VXV3dFX9dmsykpKcmy8UxfUfD+dKe+vr7T/efOnfM7rqdqamq0YcMGuVwuLVmyRPPmzTNXKAAAAIAeMx0UkpOTJUnl5eWd7vfeP8F7XE+cO3dOzzzzjBoaGnTbbbdpyZIlZssEAAAA0Aumg0JGRoYiIyNVW1vbaVjIy8uTJE2ePLlH43311VfasGGDzp49q5tuukn33nuv2RIBAAAA9JLpoOBwODR37lxJ0pYtW/xWN9q7d68qKyuVkZGh9PR03/acnBytXLlSO3bs8BurpaVFGzduVFVVlW688UY9+OCDstlsZksEAAAA0EumJzNL0qJFi3Ty5EkVFxdrxYoVysjIkNPpVElJiWJiYrR8+XK/45uamnTmzBnDXZZff/11lZSUaNiwYbLb7dq8eXOnr/fwww9bUTYAAACALlgSFEJDQ7V27Vrt3r1bubm5Onz4sKKiojRnzhxlZ2f3eCKzy+WSJLW3tys3N7fL4wgKAAAAQP+yeQKxNmgQqa6uDtjyqLll9SyZNsSwVN7QRe+HLno/dNH7oSs2Nlazxo9keVQAAAAAgw9BAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgAFBAQAAAIABQQEAAACAAUEBAAAAgIHDqoFaW1u1Z88eHTx4UE6nU9HR0crKylJ2drZGjhzZq7FcLpfefPNNHTp0SA0NDYqLi9O0adO0dOlSRUVFWVUyAAAAgC5YckWhtbVV69ev165du9Tc3KypU6dq5MiR2r9/vx5//HHV1NT0eKwLFy7oN7/5jfbt2ye73a5p06YpIiJC//jHP/TrX/9aFy5csKJkAAAAAJdhyRWF3bt3q7i4WBMmTNCTTz6p8PBwSdLevXv12muvafPmzVq3bl2Pxtq6dauqq6s1ffp0rVq1Sna7XZK0ZcsW5eTkaOvWrXrkkUesKBsAAABAF0xfUXC73crJyZEkLVu2zBcSJGn+/PlKTk5WYWGhysrKuh2roaFBBw4ckN1u1/333+8LCZJ09913a/jw4crNzVVDQ4PZsgEAAABchumgUFRUJJfLpYSEBKWmphr2z5gxQ5J05MiRbsc6fvy4PB6Prr32WsXFxfntCwkJ0ZQpU9Te3q4TJ06YLRsAAADAZZj+6VFlZaUkdRoSJGn8+PF+x5kZKzU1VR9++KEqKir6UGnnbDabZWP11jCbTfZhgXt9XHnD/ufvG30feuj90EXvhy56P3R5e3+lzzOtfj3TQcHpdEpSlysbxcfH+x3Xk7G8f+a7vK/Rk7G8Vq9ebdgWFRWl9evXS5ISExN7PJbV5kxIkJQQsNdHINH3oYveD130fuii90NVUlJSoEswxfRPj5qbmyVJYWFhne73zlnwHmdmLO/2lpaWXtcJAAAAoOdMX1HweDym9nd2rJWXTTZt2mTZWFZ64oknJEnPPvtsgCvBlUTfhy56P3TR+6GL3g9dg6X3poNCRESEpK6/5fdu77gaUndjdXX1wTtWV1ccBpLW1tZAl4AAoO9DF70fuuj90EXvh67B0nvTPz0aNWqUJKm+vr7T/efOnfM7ridjef/Md3lfoydjAQAAAOg700EhOTlZklReXt7pfu/9E7zHmRnLu70nYwEAAADoO9NBISMjQ5GRkaqtre30BD8vL0+SNHny5G7HmjRpkmw2mwoLC9XY2Oi379KlSzp69KhsNpuuv/56s2UDAAAAuAzTQcHhcGju3LmSpC1btvjNL9i7d68qKyuVkZGh9PR03/acnBytXLlSO3bs8BtrxIgRmjlzptxut1566SW1tbX59m3fvl1NTU2aNWuW4WZsAAAAAKxl8/RmWaIutLa2at26dSopKdGIESOUkZEhp9OpkpISxcTEaMOGDX73K9i5c6d27dqlOXPm6OGHH/Ybq6mpSWvWrFFtba0SEhKUlpamqqoqVVVVKSEhQRs2bNDw4cPNlgwAAADgMkxfUZCk0NBQrV27VosXL1ZoaKgOHz6ss2fPas6cOXruued6dVOz4cOHa+PGjZo7d67cbrcOHTqkixcvau7cudq4cSMhAQAAALgCLLmiAAAAAGBwseSKAgAAAIDBhaAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwMAR6AIGg9bWVu3Zs0cHDx6U0+lUdHS0srKylJ2drZEjR/ZqLJfLpTfffFOHDh1SQ0OD4uLiNG3aNC1dulRRUVH99A7QV1b03uVy6fjx4zp69KgqKirkdDpls9k0btw4zZo1S7fccoscDj6qwcbKz31H1dXVeuyxx3Tp0iVlZWVpzZo1FlYNK1jd+5qaGu3Zs0cnT55UQ0ODwsPDlZSUpOnTp+vHP/5xP7wD9JWVvT9x4oT27dun0tJSXbx4UVFRUUpPT9ftt9+uiRMn9tM7QG+VlZXps88+06lTp1RSUqLz588rJCREf/vb3/o03kA7z7M//fTTTwe6iIGstbVV69ev14EDBxQSEqLrrrtOHo9Hx44d04EDBzR9+nRFR0f3aKwLFy5ozZo1On78uGJiYpSZmamvv/5aR48eVV5enmbNmqWwsLB+fkfoKat6//bbb+vVV1/Vl19+qVGjRmnChAmKiopSaWmpjh07pvz8fP3gBz8gLAQRKz/337Vp0ybV1tZKkhITEzV79mwrS4dJVvf+0KFD+u1vf6uysjKNGTNG11xzjaKionTmzBlVVlbqtttu68d3g96wsvd79+7Vn//8Z9XW1io1NVVXX321PB6PCgoK9O9//1uxsbFKS0vr53eEnvjrX/+q9957T//5z3/U3NwsSbLb7Vq8eHGvxxqI53mceZi0e/duFRcXa8KECXryyScVHh4u6Zt/BF577TVt3rxZ69at69FYW7duVXV1taZPn65Vq1bJbrdLkrZs2aKcnBxt3bpVjzzySL+9F/SOVb0PDw/XwoULdeuttyo+Pt63vbq6WuvXr1dRUZHeeust3XXXXf32XtA7Vn7uO/rggw+Un5+vm2++Wf/85z+tLhsWsLL3FRUV+sMf/qCIiAj98pe/VEZGhm9fe3u7ysvL++U9oG+s6n1TU5N27Nghh8Ohp556yq/vn3zyiV544QVt27ZNs2fP9r0GAmfChAlKSUlRWlqa0tLS9LOf/azPYw3E8zzmKJjgdruVk5MjSVq2bJnfB3r+/PlKTk5WYWGhysrKuh2roaFBBw4ckN1u1/333+/7yyNJd999t4YPH67c3Fw1NDRY/0bQa1b2fsGCBbrzzjv9QoIkJSUl+cLBwYMHLaweZljZ+44aGxu1bds2TZw4UTNnzrS0ZljD6t6/8sorcrvdWr58ud/JoiQNGzaMb5SDiJW9LykpkdvtVmZmpqHvN9xwg6666iq1tLToyy+/tPZNoE8WLFigpUuXasqUKYqLi+vzOAP1PI+gYEJRUZFcLpcSEhKUmppq2D9jxgxJ0pEjR7od6/jx4/J4PLr22msNfxFDQkI0ZcoUtbe368SJE9YUD1Os7P3lpKSkSJLOnz9vahxYp796/8orr6i1tVUPPPCAJXXCelb2/ssvv1RhYaGSkpI0ZcoUy2uFtazsfUhISI9es68/X0RwGqjneQQFEyorKyWp0380JGn8+PF+x5kZy7u9oqKit2WiH1jZ+8vx/lbdzLcYsFZ/9P7YsWP6+OOPtXDhQiUmJpovEv3Cyt5//vnnkqTrrrtOra2t2r9/v7Zs2aItW7boX//6ly5evGhR1bCClb1PS0tTZGSkPv/8cxUVFfnty8vL0+nTp3XNNdfwb8EgM1DP85ijYILT6ZSkLlc68P6UxHtcT8b67s9PvLyv0ZOx0P+s7P3l7Nu3T5I0depUU+PAOlb3vrm5WS+//LLGjh2rBQsWWFMk+oWVva+qqpIkhYaG6le/+pXOnDnjt3/Hjh169NFHde2115opGRaxsvdRUVF68MEH9ac//Ulr167VNddco/j4eJ09e1alpaWaNGmSli9fbl3xCAoD9TyPoGCCd/Z7VzPUvb9h9B5nZizv9paWll7XCetZ2fuuvP/++zp58qSioqI4gQwiVvf+jTfeUF1dnZ566ilWtgpyVvbe5XJJ+ubLgKioKD322GPKzMxUQ0ODdu3apdzcXD3//PPatGmTRowYYdE7QF9Z/bm/4YYbFB0drRdeeMHvqkJsbKy+//3vKyYmxmTFCDYD9TyP/yuZ4PF4TO3v7FibzWaqJlwZVva+MwUFBXr11Vdls9n00EMPdfkNBK48K3tfWlqqnJwczZ49W5mZmWZLQz+zsvft7e2SpLa2Nv385z9XVlaWJCkyMlK/+MUvVF1drdLSUr333nu64447+l40LGH1v/nvvvuutm/f7ls/f8yYMTp79qz+/ve/a/v27SopKdGjjz5qpmQEmYF6nsccBRMiIiIkdZ3+vNt7sryZd6yuvo3wjhVs6+sOVVb2/rsqKyv1/PPPy+1267777tP06dP7XigsZ1Xv29ra9OKLLyoyMlL33HOPtUWiX1j5ufceEx8f7wsJHf3whz+UJOXn5/epVljLyt4XFBRo27ZtSklJ0erVq3XVVVcpPDxcV111lR599FGlpqYqLy9Pn376qXVvAAE3UM/zuKJgwqhRoyRJ9fX1ne4/d+6c33E9Gcv7Z77L+xo9GQv9z8red1RTT27tSQAABIFJREFUU6MNGzbI5XJpyZIlmjdvnrlCYTmrel9fX6+KigrFxcVp06ZNfvu8P0s5deqUnn76aYWHh+uJJ54wWzpMsvJzP2bMGEnS6NGjO93v3d7U1NTrOmE9K3v/0UcfSfpmpaRhw/y/rx02bJimT5+u8vJy5efndxoiMTAN1PM8goIJycnJktTlTXG86yl7jzMzlnd7T8ZC/7Oy917nzp3TM888o4aGBt12221asmSJ+UJhOat739DQ0OW62S6XSwUFBYqMjOxDpbCalb33Ln381Vdfdbr/woULkvp2VRLWs7L33hNF7zfM3+Xd3tXfDQxMA/U8j6BgQkZGhiIjI1VbW6vy8nLDkld5eXmSpMmTJ3c71qRJk2Sz2VRYWKjGxkbFxsb69l26dElHjx6VzWbT9ddfb+2bQJ9Y2Xvpm/8hbNiwQWfPntVNN92ke++91/KaYQ2rej9mzBjt3Lmz0335+flat26dsrKytGbNGmsKh2lWfu4nTpyosLAw1dTUyOl0Gr5FLCgokNT1Uoq4sqzsvff/76WlpZ3u9273XnXC4DBQz/OYo2CCw+HQ3LlzJX1z++2Ovzvbu3evKisrlZGRofT0dN/2nJwcrVy5Ujt27PAba8SIEZo5c6bcbrdeeukltbW1+fZt375dTU1NmjVrFuvpBwkre9/S0qKNGzeqqqpKN954ox588MEBN9lpKLGy9xhYrOx9WFiY5s2bp7a2Nr300kt+Y504cUIfffSRbDabbr755n5+V+gJK3vvnXeWm5truEHb4cOHlZubK5vNxvy0AWqwnedxRcGkRYsW6eTJkyouLtaKFSuUkZEhp9OpkpISxcTEGNZCbmpq0pkzZzq90+59992nkpIS5eXlaeXKlUpLS1NVVZWqqqqUkJDAt8xBxqrev/766yopKdGwYcNkt9u1efPmTl/v4Ycf7rf3gt6x8nOPgcXK3v/kJz9RYWGhjh07phUrVig9PV1NTU364osv5PF4dMcdd/ideCKwrOr9tGnTdMMNN+iTTz7R73//e6WlpWn06NGqq6vzXU244447NHbs2Cv23tC1Y8eO6a233vLb5na7/a72Ll682Hc1abCd5xEUTAoNDdXatWu1e/du5ebm6vDhw4qKitKcOXOUnZ3dq0kpw4cP18aNG7Vz504dPnxYhw4dUmxsrObOnaulS5dyO/cgY1XvvRNX29vblZub2+VxBIXgYeXnHgOLlb33jvXOO+/owIEDOnHihEJCQpSZmanbb7+9xz9dxJVhVe9tNptWrVqlDz/8UB999JFOnz6tiooKRUZG6vrrr9e8efM0adKkfn436KmmpiaVlJT4bfN4PH7berrowEA8z7N5zC74DgAAAGDQYY4CAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAg/8Pgcm5kIcVSV0AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 900x600 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(x,y)\n",
"plt.plot(x,y_ue, c='k')\n",
"plt.hist(data, 20, (0,1), density=True, alpha=0.3)\n",
";"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"''"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAwoAAAIJCAYAAAAF28o7AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAXEQAAFxEByibzPwAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3Xd4W/W9P/D3kWTJWl7ynvFI7OztOCEDwmggJFBGQinclj3CbSnQcnuhQC+0pVDKLbcQ4EehBRpWgNCkSQiEJJC9t+N4721rWMOyxu8PE8WKHMeOZR9Ler+eJ0/01TjnYx8fW2+d7xDcbrcbREREREREPUjELoCIiIiIiEYeBgUiIiIiIvLBoEBERERERD4YFIiIiIiIyAeDAhERERER+WBQICIiIiIiHwwKRERERETkg0GBiIiIiIh8MCgQEREREZEPBgUiIiIiIvLBoEBERERERD4YFIiIiIiIyAeDAhERERER+WBQICIiIiIiHwwKRERERETkg0GBiIiIiIh8MCgQEREREZEPBgUiIiIiIvIhE7sAsTU0NMDtdg/7fmNjYwEALS0tw75vEg+Pe+jisQ9dPPahi8c+dIl17AVBQGJiot+2F/JBwe12ixIUeu6fQg+Pe+jisQ9dPPahi8c+dAX6sWfXIyIiIiIi8sGgQEREREREPhgUiIiIiIjIB4MCERERERH5YFAgIiIiIiIfDApEREREROSDQYGIiIiIiHwwKBARERERkQ8GBSIiIiIi8sGgQEREREREPhgUiIiIiIjIB4MCERERERH5YFAgIiIiIiIfDApEREREROSDQYGIiIiIiHwwKBARERERkQ8GBSIiIiIi8sGgQEREREREPhgUiIiIiIjIh0zsAohCzfayVhgMJtH2n5+qFW3fREREFDh4RYGIiIiIiHwwKBARERERkQ8GBSIiIiIi8sGgQEREREREPjiYmYiIaIhxEgMiCkS8okBERERERD4YFIiIiIiIyAeDAhERERER+WBQICIiIiIiHwwKRERERETkg0GBiIiIiIh8MCgQEREREZEPBgUiIiIiIvLBoEBERERERD4YFIiIiIiIyAeDAhERERER+WBQICIiIiIiHwwKRERERETkg0GBiIiIiIh8yMQugIiIiChYbS9rhcFgEm3/+ala0fZNgY9XFIiIiIiIyAeDAhERERER+WBQICIiIiIiHxyjQEREFGAEhx0SqwVd+la4XQ7II2IApQpOeTggkYpdHhEFCQYFopGkw4jOmjJ0NNaio7EexuZGmNpaYNC3wWqxIFyhgFKpglKtgUqtgVKrRbhaC13mGGgm5cOl0oj9FRDRILkdXbBVnEbTiUNoKC5EbVUF6hsbYbRa0NFph7GzCyaHAx0OJ5zus69TSSVQyyTQyGRQhcmgVSiQGheL9PRRsEybjDHTZyJx4hRIlCrxvjgiCigMCkRicLthr69C/f7tqD5+COWlJSipq8NpvQk2l/vCr++FTi7DhNho5GVmImf8ZGTMmg957mS4w+R+Lp6I/MXVaUXLvu9QtncHThceRWltPSqMJhi7nAPelsXpgsXpQnOn4/t7OoDGVuB4EbD+SwCAVibB6KgITMhIw9x5czHn+psQlZMLQRD8+FURUbAQ3G73xb0rCRL19fUQ41sQFxcHAGhubh72fdPwczkdKNmxDYe2fINvtu/EkcoqVJttQ77fNJUCC8fmYu51yxC7cCkQJudUeSLhOR+64uLivp8i0wBLVSnKtqxH8ZEDKCwtxYnmVlicLtFqEwBM1EVi7vixmHv55Zh53Y1QxSWIVk+w6XnsxcLf+eIQ63e+IAhISkry3/YYFBgUyP8cXXac3PI1dn25AXsOHMS+ylq02bv6/fowiYAYZTii1SpEabSIioiASqNFZ6cNVosFZqsVFpsVZlsnTLZONFusF9zmmAgVrpo+A7fffx9S5syHIOFcBsOJ53xoqj91Eke+/hJffrUZB4qLUWro6Pdro+RhyNBFIy0xCUlpoxAZnwBlZDSUUTEIj46FIiYWqigdJFIprCYDOg3tsBvaYDMa0Gk0wNTahNrS06iqqUJFUzOaLBf+cEIhEXBpZiquXfQDXHnHvdAmpQzmyw95DAqhi0EhSDAokD84HA6c3LYFOzesw679+7G3orpfXQc0YTKMjotFdkYG0vMmIGXqLCSMmwKVRjugrgDG9jZU7N+Bin07UFx4AoXV1Wi1nv9Nwaz4aPxo8dW47pdPQB4Z1e/90MXjOR8aak4cxa61X2D3rp3YXVSCCpPlgq+RCsAYXTTyMrOQPWEyksZNQdyEadDGxPqtrvxULdrb2lByYB8K9+3Grp07saPwNFpt9vO+RiERcFl2Oq69+mpcecc90MQn+q2eUMGgELoYFIIEg4I49taIt0olMPhfnBZzB458vQkHtm3BvkOHsbe8CsYuR5+vEQCMjY3GpNxcjJ46C2lzLkNs1tD0DXa73Wipr8GxNR9g66Z/41h9Y6/PS1Mp8OB1i7HsiWcQHq3zex10Vqif88HI7Xaj+vgR7Fr3BXbt3IXdp0tQ3XHhq3vxSgUmZ6RjzNgJyJx5CVILFkChHto3c739znO5XCjctwfb//U5duzchd1llTA7ev+AQyERcMWYLNx888249Kd3IyxcOaT1BgsGhdDFoBAkGBTEEUhBwel0oqKoEMe3foP9O3fiQGEhTjS2wnGBnxuZIGBSgg6zJoxDwfxLMf3apRgzYZIofzQaS4uw54O38M2321BpMPo8nhAux31XX4Hbnvwt1InJw1rbcBLz5y4yMhJzs3Qhe84D4n7//fFmydHVhaLt32LPVxux98AB7CurQIOl84KvS9UoMS07G3mTp2P0pYugGzdl2AcP9+fr7+zsxM7PPsHajz/CxsNHYbD3/uFHrCIM18+ajmV334txC6/iQOg+MCiELgaFIMGgII6RGBTcbjeaGhtRfGAvCvfuwamTJ1BYXonTTa2wOi/cjUgqAJPiYjBn4gTMvmwhZi69HhpdnNdzxP6j4Xa7EVZ+EP/4y5/x+Z4DXlMrAkCMXIa7r1iAu597AeqE4OtmwKAgrkALCobmJhz+cj0O7NyOg8dO4EB1bb+6FGZoVZidm4OC2XNQcO1STLvsioB7s9hps2H76g+xbvUn2Hj4+HmvmI7TReKmq67ADx/8GeKzcvxRalAR+3c+wKAgFgaFIMGgII7hfsPgdrthNZuhb2uBvqkBMR11qCkuRlV5GWpq61Dd3IxqvRGdA5h9RCWVYEpiLKaPzcPMS+Yif+kN0Cb2fXKOhD8aZ7SUF2PDX57Duj17YT9nStb4cDnuuv6HmLHivyGRK/y6XzH/aDEoiGskBwWnw4Hi3TtwcMtmHDxwAAdKylDcbkR//jpkRahRMCYHBZdcgtlLrkPy2Alej4+E834w512nzYpt//wHPvngQ3xdVOLz+wLo/qDk0uwM3HzTzbjyznsQruaaLkDgH3u6eAwKQYJBQRz+esPQabOhpbEO7a0tMLS1wtDeCkNzI0wNtTA2NUDf1op2owFtpg7YHH2PIeiLACBDE46p6amYNmkSZl56GfIWXoWwAf4xHAl/NM6lr63Cpleew5rvtsN6TlDKi9birntWIOeGnwB+6l7AoBCa5zwwsoJCY0kxDn21AQf37sGhU0U4WteEjvP0z+9JKgDjY6ORPzYX+bPnIP+aJYjLGdPna0bCee+v866trhZfvPYKPl2/AYcaW3t9TqRchqUzpuLmO+/GtEWLQ7prUjAdexoYBoUgwaAgjoG8YegwGlB+uhDV5SVoqqtFU0Mdmhvq0FxfC0N7m99riw6TYawuErlpqRiXl4u8adORO3s+VKlpg/6DNxL+aJyPqbkBX738P/hoyxafTwznZ6Tgtl8+g7iZ8wa9HwaF0DznAfG+/06bBcqi3TiwbQv2Hz6CgxVVqOnnOiYRYTJMS0nEtHF5mDlvPqZdvQSauPgB7X8knPdDcd4V79qOT954DZ/t2I3684zVGBWhxrWXzMbSn9yJcXPnh1xoCNZjTxfGoNBDWVkZjh49ipKSEhQXF6O9vR1hYWH45z//eVHbM5vN+OSTT7B3717o9XpERUVh5syZWLZsGdRq9WDL9cKgII7zvWHQt7Wi5OQxlJ8uRFnRSZQVnURTfa3f9isBoFPIEKcIQ0qEFqlxOqQnJyN1VCbSR49B6rgJiMweA2GIZvQYCX80LqSt7DRW//6/sOHYSa/7pQJww7QpuO43L0KdlHbR22dQCM1zHhi+77+tqR7lm9fh1P5dOFlajOONrbC5LtytUAJgTHQEpo/OwrRp0zH98iuQlT8HUplsUPUEwnk/GM6uLpT9+yNs++ITbD5VfN7vdU6UFkvmz8WSO+5Gbn7BMFcpjpFw7BkUxMGg0MMLL7yA/fv3e913sUHBZDLhySefRH19PRISEpCVlYWamhpUV1cjMTERv/vd76DV+u+HnkFBHGfeMLhcLpQUHsfBndtwYMc2lBQeH/C2pAIQKw9DnEIGnTwMsQoZYsMViIyKRpQuFhFxiYhITEFEcjrCU9LhjI7H5Lx0CIP8438xRsIfjf6q3P413n3pWRys855aNTJMhjuvux6X/Pw3kF7E+AUGhdA854Gh+/7b6qtRunkdCvftxOHi0yhsNaA/o40SlApMTU/B1AnjMfWSuZh8xQ98JiDwh0A67wers7URh//5Jr7evAn7G87/s54THYGFM6Zi4bVLMXPxdQhXBud0q55jr9fDZdLDVFkCY3U5bEY9uqxWdNms6Oq0wm7rRJfdBrfTCZU2AuqoaKijdVBFx0IVmwBlbDyk8SmALGzANTAoiINBoYc1a9bAbrcjOzsb2dnZuPfeey86KPz1r3/Ft99+i/z8fPziF7+AVCoFALz99tvYuHEj5s+fj4ceemiwJXswKAw/q9WK1z9ehwM7tuLgzu9gaO+9n2tPkWFSjI9QIV2lQKpSjhSlHKlKOZIiI6DOGYeupFHojElEZ0wCbDEJ6NLGAH2sPCzWL85Ae8Pgdrlw9MO38M7br6PynIWjRkdF4J6fPYrRi5cPaJsMCqF3zp/hr++/w2JG1ddf4Oi2r7Hv5AkUthkuOOg4TBAwMUGHabljMKOgANOuWISUceP9Us+FBNp57y/GoqM48PE/8O2enTjc3H7e56lkUlySk4nLLl2Ay5bdivTcvGGs0n86mhtReeggyk8eR0XxaVTV1KChrR117QY0mS1oPc90s/2lk8uQGxuD7LQ0ZIwei7QpMxE7rQASdWSfr2NQEAeDQh+WLVt2UUFBr9fjvvvug0QiwcqVKxEVdXbF2K6uLjzwwAPo6OjA66+/7vXYYDAoDJ9Tp07h/fffx6effgqj0Xcu/zPiwuWYqFViYqQK4yOUmBCpQkq4HIIgwBqbDHPqaHSkjYY5dTRsccmAcP5AcD4MCgPTZbNi20tP4x/r/+0z4HPh6Czc+pvnocud1K9tMSiEzjl/rov+/rvdaD++Dyc2foEDB/ZiX1UtzBeYoSxKLsPU9DSMHT8RObPmYemNS6DU8LwXQ36qFpV7d2PtP97C2u924nhr39+HBFU4po5Kx7QpkzF13gJMuuxyaLQRw1Rt38wGPSoO7kfZ8SMoP30a5RWVqGhoREWbHk19rHI9VOQSAWOiI5A/YRKmLLoOKQsWQSL1vlrOoCCOYAkKw9/3og+HDh2C2+3GuHHjfIJAWFgYpk+fji1btuDw4cO49NJLxSmSBsRqtWLdunV4//33fbqnnSGRSjEhMxOLwp24MlqJMZpwrwFvnZGxqJ8yH61T5sMe5f9uAXRhYeFKXPHEC5h1+3349Llf4V9HTng+wf2muAzbf7ocN86Zg2sefw6aeP/9gqLQ5XJ0oeabdTjw5b+w4+hRlBrNfT4/WhGGaRnpmDhlOkYv+AESpxZA8v0VaQCihQTqlpFfgIfyC7DC7Ubl7h34ZvVH2LJ7N3ZWNfiMaWi02LDx5GlsPHkaWPUJJAByY6ORl56KUenpyMgZjczxE5AxaQpiExL9OkDaarGgrvgUagtPoqa0BDVVlaitq0dNUzMq2tr7tcDehUQpwhChUCBcJoM8TAaFLAzysO5/EASYrRYYLVaYbDYYO+3nXS0bAOwuN463GnB823fAtu8QF/7fmJOXi+mXLcKYa2+GQjMyAhYFrhEVFCorKwEAmZmZvT6emZmJLVu2oKKiYhiroovR3NyM119/HR9++CH0er3P40qVGrMWXIGCiROxtKUIqQ1lXo+7ZGFoH5uP1qkLYBo19qKuGpD/adOz8dM3P8UVWzfgby89h8NN3d3G7C43Pti+A+v2XIlbr74Glz38GyjUfGNGA9Nl0qN47UfYt3UTvjt1Gq2dXed9rkIiwfS0ZEybno9xVy1F4uR8CH10N6SRQRAEjJo9F3fOnos7AVham7H704+wddMmbDleiLJzujgCgAtAYUs7ClvagYPHvB5Ty6RIjtAgUqVEpEqFCI0aEVotIrQR0EZoATfgdDrhdDnhcjrhcrngdDph7uhAu8GAdqMJ7WYz2i1W6G2dMA6ye5BCIiA9UotRcTpkpKQgKzsbrsgEKOMSoU3PQmTKKMjDwwe0za4uOzraWtFwaBdqjh5EefEplNbUoKzdAMc5PSKabXZ8cfgYvjh8DPK//AnzckfjrgdXYN71N0DC84MuwogKCi0tLQCAmJiYXh/X6XRez6ORp62tDa+99hreeecd2Gy+0w9OmjQJt912G9KmzEfWgU1I2LEOEtfZT0ucYQrUL/ghWqYvhDPcvzNckf+kXno1npp/FQ6++Se89eEqNFi7P2UzdTnwxr/+hTVfb8Ltt9yGgjt/DmnYwAffUeiwtTfj5GfvYfeWr7G9rMJnLY+esqO0yB8/ARPnX4nsq5ZCruKiXoFOpYvDwnsfwsJ7H4Lb5ULTiWM4uPlLHNq/D4eKinG0qRUdjvP/TJgdThS3GYC24evWJZcISI9QIzMuFlmpqRiVlYmsseMxavJUJOWO9YytBPzT7SwsTI7ohCREL7oBYxfd4Lm/q9OGhsN7cWrzv7Fn7y4cqm/yGsRvd7mxufA0Nv/nz5H25BP48fXXYdnPfoGECywMStTTiAoKZ95YKhS9z6Ry5v7Ozv5f+nvkkUd87pPL5Xj++ecBALGxsQMt0y9k38+4c6YPW6Bra2vDyy+/jFdffRUdHR1ej6nVaixfvhz33HMPpk+fDvuJw2j48/9A3lrv9byOCQVo+eG9cETHY6j//MfF6YZ4D72TyWSQSqWIjOx78FmgWPj47zH3vsew4blf4u2NX8L0/SXyRosNf3r7LeSs+QS33XEf5vxkBWRyOQDxvvcAEGkS7xM1qVQKmUwWNOf8xTjz/e+oq8b+D/4fvt28CbsqqtDZy0q/QPeMZjNSkzF37nzMuOk/kDB+6kXvW8yfu2A77wdqIN/7+IQETFh4Bf7j+3aXuQMnNn+J/Vu+QXHxaZRV16CiuQXlhg4Yuy68QN7FkKB7RqzUqAikxsYgPTkZ6RkZyB43Hrkz8pExZRpk/ZzxbWiPfSRir1qCCVctwU0ATLWVOPzJ37Fry9fYUVru+X0MANWGDjz/j3/iT++uwtX503HvL/8LVy1ZEnLrWgynYHmfN6KCwplBxfzBDRwmkwkvv/wyXnnlFZ8BynFxcXjsscdw1113ISKiu5+kbc+30L/4JOSOs10KuqLi0HLDfTBPCI15tYONPCoG1/3pb7jiZ2X49JlH8MHuvZ43fiVtBjzz0gtIfP3/sOy6G3DFz54AIN4bNhJPdeEJfPniS/ju22+wt7rep8vEGWqZFHNHZ+OSy3+AKcvu8NuYl+1lF55dbaj0/ISZBiZMrcGUpTdiytIbve53d3WhufgUivfvQ311BQytrWhvb4fBYIDeYISxowMmiwWCIEAqkUAiSCCVSiCRdP9TK5WIiY6CTqdDTGwcdAkJ0CUlIy41HWkTJkER6Z8JU4aTNiUD8x5+GvMefhqODiOOvr8S//7sI2yrqvdcaXC43Vi7Zz/W3nQTxqUm47HfPIXlt92OMF75pfMYUbMe/elPf8LevXvx05/+FNdcc43P4/v27cOLL76I/Px8PPbYY36plbMeXRy3243PP/8czz33HBobvefZj4qKwoMPPog77rgDKpXKc79r91a43/lf4PuBa25BgsaCRai/9Ea45APrszlYnPVo6BiPH8Cnf3oGawuLfaasVMuk+NHlC3D3E08hLXv0sNc2JLMeuVwQXE64BQE48w9n/j8rFGc9Kj24Hxvfewcbv92Ogw3n7zIaJZdhfl4u8q9cjDHXLkeYKri6HZ75NDmYz/u+hPKsO6L/zne7kVJ5AB+8+go+3HsI9TbfcT8pumjcs+Ih3Hrb7X5f1DaUcdajIXCmG1BbW1uvj7e2tno9jy7eYN4wlZ8+hbdeeg4nD3vPYqTWRuC6W+/A4mW3Q6XR4HibE2jr3k/svq+Rvv7vEL5/6+hQqFBy62Mwp+de/BdBI1LEhOm44+9rce23G7Dub69iQ1Epur4P42aHE299+Q3e2fQNFo7LxeKl1+GKW29DdMzIu8rgMhthKD6J1tJTaKksQ0t9LUwmI0zm7k8qDVYbjLZO6O1d6HK5IRGA7oggQBC6uy/IJBJEKeSIUikRo9HgXV00IiO00Ol0SErPQOa4CciYPA3KuPiguJJqs1qwZ81n+Obfa7H10BGU6M//eyYuXI4F48dh5g+WIuvqGyAb5g8LiEKCICBl7mV4bO5leLipAVte/wve/WIdvmk4+z6rtrUdz/zPs3j5pZfw07vuxl133+0ZE0o0ooJCRkYGAKC8vLzXx8/cf+Z5NLxMBj0+ePMVfPnZh3D1mM4uTC7HdT++E9f/+E6oe5nrOmH7WqRu/tDT7lJFoPj2x2FNHDUcZZNI4uZfjTvmX43lJw9h85svYfX+Q54+xU438NWJInx14gXInn8Rc3Iycc3ia/CD236KeD9+EnIhjk4r2k4eQdPJw2goLUJNdSVqm5tRZzCiwdLZr9V9++ZEc2cXYDQDaAFKKnp9Vny4HKOiI5GREIfM9HSMmTABY2bkI2PqTM/YjpHI7Xaj4uRxbPvkQ2zZtg07Siv7HIycqlZiweRJmLHoeqQtvJYD3YmGkSw+EVc+9Qdc8diTKPzH63jj3fexproJju8v/RrMFvzllVfw5htv4Paf/AT3338/EhISxC2aRDeiuh61t7fj/vvvh1QqxcqVK70G/5xZcM1kMuGNN97ggmuDNJArCm63G5vXfor3Xn0JRr336pr58y/HHQ//FxJT0np7IZK/+RhJ2//lucseEYPTt/8anbHJF127P7Dr0fBzNFZj9xsv4eMtW1B9nrnIBQBTU5MwfnQOxk6ciLH5BcidNmNQAwGtFguqTx5H1fGjOHj4OBprqlDXWI+q5lZUd5jhHP7Tv98UEgHZ0REYnZSIMTnZGD1uPHKmzcCoKdOhUCqHvR6r1Yrj27Zg3zdf4cChQzhQXoVma9+TS+RFa7Fo1kxcc8ut6MidE5JTmIZ61yOxidn1aST8zj/f1+826lG76m289c9VWFXZ6LOIoUIuxy0/+hEefPBBpKamDkepQSVYuh6JEhQ2btyIjRs3Ij8/H7feeqvXY6+88gq2b9+OWbNm4eGHH/YMAnvnnXewYcMGzJs3D//5n//pt1oZFPpWW1WOlX94CicO7vO6PyktA3c98gSmz5nf+wvdbqRteBfx+zZ57rLFJKD49l+PiEXTGBTEMy3SjW/feRMbNqzHplOlaOnHvOUpWjVyU5IQFREBZbgCSqUKSmU4VCo1lCoV7J02GNraYTDooTcaYegwo73DjJYOMxosvtP09peA7lVik6MikRAbi6jIaGgiI6GJjIEqOgZqXRyUMXEIU6rgdrvhdrvgdrkAV/f/XZ02mJsb0NHcCJtRD5fFgMamJrTo9ahuaUe1scPTLWsgpAKQrtUgJzEeOaPSMSo7B8mjspCUMwbJuXmIjI6+6K8ZACwWCypOHkPZoUMoLypEeXk5iquqcbyhGfbzzFB0RrhEgktGpeCyS2Zj4U23IH16vqdblZgrY4uJQUFcDAp9f/3utma0f/J3vPfRJ3irvBFt5/xOlslkuOmmm7BixQpkZWUNZalBhUGhh4MHD+LTTz/1tIuLiyEIAnJycjz33XjjjZg2bRoA4OOPP8bq1auxYMECrFixwmtbRqMRTzzxBBobG5GQkIDs7GxUV1ejuroaCQkJ+N3vfueZQccfGBR653B0Yc37f8PHf3sNXfazy9KHK1W4+c4HsOSWnyCsjy4R8bvWI23T2aBojU/F6dt/DYdmZMwkwaAgnp7fe4dRj70f/RMb1q3FxhNFqLPa+3jl0JAASNOqkBEbi9TkFCSNykLsqBzE5IxDdHYuwhT+6Tvf22Bmh8OBuuLTqDx2GBWFJ1FRVobSqiqcrm9ClcnsMxi8vzRhUiRq1EiKioQyPBxKhQLK8HCEhyugVCoRHq6E3d6Jjg4zOsxmdFgs6LBa0WGzocnYgXqzdUD7y4lQY8HEsVh41SLMuukWKKN6DyoMCqF73ouJQaF/X7+7oQYd776GVZs24/WyRjSds+ChIAhYvHgxVqxYgUmTJg1FqUElWIKCX8YoGI1GFBcXe93ndru97jt36szziYiIwB/+8Ad8/PHH2LdvH/bu3YvIyEgsWrQIy5Ytg0bDBXaG2ukTR7Hy979BRUmR1/2zFlyBux99ErEJiX2+XlVbitSvz45JMCdnofi2x+FUjpxjJ9YbFjHn8R+JZBFRmHPPCsy5ZwV+azHj6BercWjXTpwqKUFhbT2K9KY+F1vqL6kApKiUSImKRHJ8HBKTUpGQlYOEcVMQM24qZErVhTfiB91vGM752dOmQTEnDblzlqDn0H672YjmI/vQcPwwqksKUVVdjYqWFlSZLtxdqqPLiZJ2I0ra+/d7dyDCJRJMjo/BjLwxmF5QgOlXXg1d7tigGIxNFMqExFRofvk73HvFEty26g18crQIK8saUPP9Bzhutxvr1q3DunXrMHfuXKxYsQLz5s3juR/khqTrUSDhFYWzrBYzPnjjFfz74/e8BitH6+Jwzy9/g9mXXXXB7UptZox94wko9N1fl10bjZO8ld2YAAAgAElEQVT3/wFOVehOj9cTP1kc2Kd7LqMeNQf2onD/XpSXlMBiscDa2QmLzdb9f6cdFrsdYVIpojQaRGo1iIqMQmR0FCKjY6BLSETa2PFImTgJYZoIUT/R9texd1jM0J88iIYTh1FXeho11ZWob2tDk6kDDdbOC3YN6i+5REC6RolMXQwykxOROWoUJhXMxrjLfwB55MV1beIVhdA978XEKwoD//rdnZ1wb/gE9g2rsaaqCa+WNqDM7DseaeLEiXjggQewePFizwJj1I1XFCioHNz1HV7/49Norq/zuv+q65fhPx56rNfZjHy43chY+5YnJLgFAeU3PsSQQBdNEhGF9MuuQno/QmqokKnUiJ0xD7Ez5mHCOY+5O22w1ZRDX34a7dUVMDQ3wG6zotNmg91mg+DqgrWzO2TJZWHQqJTQqtVQqzXQajVQa7WIjotD5vhJSJk0BbKIkdFVkIiGl6BQQLj+NijmLMTNH76FG47uw6ZGPVaWNuCwweJ53rFjx/Dggw8iMTERt956K2699Va/vkkl8fGKQohfUTDq2/H2y3/Ato3/8no8JSMTD/z6fzB+6sx+bzN239fIWP+Op1172c1omH+9fwoOEvxkUdxP94LhisJgiL3wFa8ohO55H6qC4ti73Yg+uQdp6/8OmdmI3W0dWFnWgK3Nvl0bJVIp8uctxKIbf4SJMwpQkH7xM9YFOl5RoIDmdrvx7Zfr8PbLv/ea8lQqleGH/3E3br7jAcgVin5vT9lQibQv3/e0jZnj0TB3qV9rJiIiomEmCGgfXwDTqHFI2/B3zD6xB7N1Wpw0WvBGWSPWNujh+L67ssvpxO6tX2H31q+QlJaBn/z4R1i6dCnXvwpgvKIQglcUqqur8cAvfolDu77zun/0+ElY8d/PIiNnYKslS+w2jH3zSYS31gMAutQROHn/H0bMDEcjSVB8ujRIvKLAKwqhZiQcexJHMB77qJN7kb7+HYSZu68oNHd24cN6A95t7EBja2uvr5k8eTKWLFmCJUuWhMyaDLyiQAHHZrNh5cqV+Otf/wqb7ezc8uFKFW69/2Fcc/OPPetWDET6v9/xhAQ3BJTfsIIhgYiIKAjpx+XDNGos0ja+C92xnYhThOE/R8XiwQwdvohcgPerWnBg13deH8IeOXIER44cwXPPPYepU6di8eLFmDdvHsaNGwdJCC7CGEgYFELEpk2b8Mwzz6CystLr/qkFc3H/488gPvniEn7Mke+gO7rd026YtxSmrHOHWBIREVGwcKq0qLhhBfS5MzBq7f+DtNMKqSDgBmMlFuWmY9fd/8SXu3bh0LebUFTkPdX6oUOHcOjQIQBAVFQU5syZg7lz52Lu3LnIysridKsjDINCkKuoqMBTTz2FzZs3e90fFROLn/78ccz/wbUXfVJKOq1I3bTK0zal56Lu0hsHVS8REREFBv34WTiZnIms1f8HdV0ZAEDVWIVLv/grsq69EznP/BeKioqwdu1a/Otf/0Jpaan36/V6rF+/HuvXrwcAJCYmYsqUKRg/fjwmTJiA8ePHIzk5meFBRByjEKRjFAwGA1auXIk33ngD9h4rK0ulUtx555249Jb7oBrk4nVJW1Yj+dvPAQDOMAVOrHgBXZGxg9pmsAvG/qoDxTEKHKMQakbCsSdxhMqxFxxdSPn6QyTs2eh1f/O0y1B1zU8BqQxutxuVJaexc/NGHNm7EyWFx7zWbDofbWQUMseMRVpmDpLS0pGYko7E1DTEJ6UiTC7v87Vir6EBcIwCjTAmkwl/+9vf8Oabb/r8Ypo9ezZ+97vfITc3d9B/sGWmdiTsWu9pN86+hiGBiIgoBLllYahZdDtMo8Zi1BdvQGbrXmsh7uAWyA0tKLv5Z3ApVBg1OhejRufi1vt/DnOHCScP7cOx/XtwdP8uVJac7nXbJoMeR/ftwtF9u7zul0gk0CUkIiE5Dbq4BMTExSM6Nh66uHjExMUjJi4BVl0GlErlkH/9wYxBIUiYzWb8/e9/x2uvvQa9Xu/1WGJiIp566iksXbrUb5fvkrd+BmlX9yqNXaoINM5Z7JftEhERUWAy5M1AYWIGMj/9KzQ1JQCAyNJjyH3nWZT8+Ffo0p5d1V2t0WLmvIWYOW8hAEDf1oqiY4dRfvokyk+fQvnpQjQ31PW6HwBwuVxorq/zWSj2XCqVCjqdDjqdDjExMZ7/IyMjERUV5fUvMjISWq0WGo0GigFMER/MGBQCnMlkwqpVq/Dqq6+i9ZxpybRaLe655x7cf//9UKvVfttneHMtYg9t9bTrF/wQLgUTOxERUaizR8Xh9E+ewKg1ryPmxB4A3eMW8t56GsU//hVs8b1PnhIVo8OsBZdj1oLLPfcZDe2oKC5C+elC1FdXoqGmCg01VWhuqOtXtyUAsFgssFgsqK6uHtDXERYWBo1G4/mnUqkQHh7e6z+JRAKpVAqpVAqJRAKZTAatVgupVIolS5YgISFhQPseSRgUApDb7cb+/fuxatUqrF27Flar1etxtVqNu+66C/feey+io6PPs5WLl7L5Iwju7hPUFpOA5ukL/b4PIvK/UB0jQETDyy2To/zGh2CPiEXirn8DAOTGVuS+/VuULn8YHZnj+7WdiMhoTJpRgEkzCrzudzi60Fxfh/qaKjTX16KtpQltzWf+NaKtpclrMdmL0dXVhfb2drS3D247+fn5DAo0PFpaWrB69WqsWrXKZ+YAAFAqlbjjjjvwwAMPICYmZkhqUFcVIarogKddu3A5IOWPEREREfUgSFB71a2wR8UibcO7EOCGrNOC0e//ERXX34f2iZdc9KZlsjAkpWUgKe38Kz532e3IUtrR2tqKtrY2tLa2ev7p9Xro9XoYDAaf//09wU2grxPBd3gjmNVqxeHDh7F7927PP4fD4fO8qKgoLF++HA888IBnlP2QcLuR+tXZ6VDNKdnQj8sfuv0RERFRQGvOvwpd2mhkfvYqJI4uSFxOZH32GiocXWideumQ7TdMLkdSkm5AMwC53W5YLBaYTCZ0dHSgo6MDJpMJJpMJVqsVNpvN519nZydcLhecTiecTidcLhccDgfCwsLgdDqHpGfHcGJQEJnb7YbJZEJLSwtaW1vR0tKCI0eOYM+ePTh8+LDX1Kbnmjt3Ln70ox9h0aJFCA8PH/Jao07t9wxOAoCaK38EcG5jIiIi6oN+7Eyc/skTyP7gJYRZurtAZqx9Cw6VFobc6SJXd5YgCFCr1X4Z1ynW9Kj+xqAwzFpaWvDQQw/BYDCgpaUFzc3N6Orq6vfrExMTsXz5cixfvhwZGee/5OZ3TgdSNn/kaerHTEVHxtjh2z8REREFLHPqaBTd+QzG/P1ZyDv0ENxuZK3+P5y+/dcwp+eKXR6dB4PCMFMoFPjuu+/6/XyVSoUZM2Zg1qxZKCgowMyZMyGVSoewwt7FHtyK8NZ6AIBbEFB7+S3DXgMREREFrk5dIkpuexy57/wPpJ1WSBxdyPngTyi64+nzzoZE4mJQGGYajQZyudynS5FEIkF0dDR0Oh3S09NRUFCAWbNmYeLEiQgLCxOp2m5CVyeSt33mabdOWcATmoiIiAbMmpCOklsexej3/wiJswsymwWj338ep+56hgu3jkAMCsNMEAT88Y9/RGJiIuLj4yGVSqHT6RAZGSnKlYL+0B3dgTBz9yrPLpkcdZfeKHJFREREFKg6Ro1F+Y0rkPXJXyC43ZCb2jHmvedx6s6n4VRpxS6PemBQEMGyZcsCZ5CL2424fV95mi1T5qMrYmimXiUiIqLQoB87E1XX3IGMf78NAAhvrUfOqhdR/B//DZfcPxO0iLl2TKRJgrlZOtH27y+BPbkrDTl19WmoGqs87eaZV4pYDREREQWLlhmXe/VS0NSWIvOz1wA/r2VAF49BgfoU3+NqgmnUWI5NICIiIr+pn/9DNM24wtOOKjqA2INbRKyIemLXIxFtL2uFwSDeZbELkXUYEHVyr6fdxKsJRERE5E+CgOqrfwJFexMiS48CAFK/fB/GzPGwxySIXBzxigKdV+zBbyBxOQEAdm009CNoURQiIiIKEhIJKq+7Fw6lBgAg7epE5pqVgMslcmHEoEC9czkRd+AbT7Nl+kJAygtQRERE5H9d2mhULr7T09ZUFyNxx1oRKyKAXY/oPKKKDkJubAMAuCRSNE9fKHJFFCzEnIWCiIhGLv34WWgtugS6YzsAAElbP4UhZzKsSaPELSyE8YoC9arnlKj6cflwaKJErIaIiIhCQfU1P4H9+2nYJS4nMj9/DYLDfoFX0VBhUCAf4c21iCg/4WlzEDMRERENB2e4GhXX3edpK5trkbL5YxErCm0MCuQjbt/XntuWhHSY08aIWA0RERGFElPWBDTOWuRpx+/eCE2PDzBp+DAokBdJpxW6I9962s0zrwQEQcSKiIiIKNTUXr4c1thkAIAANzLXvAFJp0XkqkIPgwJ5iTm6A1K7DQDgUKjQNnGOyBURERFRqHGHyVF+w4NwS6QAALmxFQk714tcVehhUKCz3G6vlZhbp8yHSx4uYkFEREQUqqxJmWi4ZImnnbBrPWRmg4gVhR4GBfLQVJ6CsrnG026eeUUfzyYiIiIaWg2XLPZaiC1xO9dWGE4MCuQRe3ib57YxawI6dUkiVkNEREShzqVQoWHuUk87bt/XCDO0ilhRaGFQIACA4OhC1Kn9nnbLtMtErIaIiIioW9PMK2HXRgMAJM4uJH37ucgVhQ4GBQIARJQdg7TTCgBwhilgGD1F5IqIiIiIugc218+/3tOOPbQNitYGESsKHQwKBACIPrHHc9swegoHMRMREdGI0Tr1UnRGxwMABLcLyVtXi1xRaGBQoO5uR0UHPO328bNErIaIiIjIm1sqQ92lN3ra0cd3Q9lYJWJFoYFBgRBRym5HRERENLK1TZgDa1wqgO5F2JK3fCJyRcGPQYEQfbJHt6MxU+AOU4hYDREREVEvJBLULrzZ04wqOgh1TbGIBQU/BoUQ59PtaBy7HREREdHIZMidDnNylqed/A2vKgwlBoUQx25HREREFDAEAbWXL/c0I8pPQFt2XMSCghuDQoiLPrnbc5vdjoiIiGikM2WOh3HUOE876bs1IlYT3BgUQpjgsCOq6KCnzW5HRERENOIJAuouu8nT1FYUQtFSJ2JBwYtBIYSx2xEREREFInPaGFji0zzt2INbRKwmeDEohDDOdkREREQBSRDQMn2hpxl7+FsIDruIBQUnBoUQJTjsiDrVc7ajAhGrISIiIhqYtkmXwCWTAwBk1g5EFe4TuaLgw6AQoiJKj0FqtwE40+1ossgVEREREfWfM1yNtgmzPe24A+x+5G8MCiEq+kTPbkdT2e2IiIiIAk5zj+5H2koOavY3BoUQ1D3bERdZIyIiosBmScmGJSHd04478I2I1QQfBoUQFFFylN2OiIiIKPCdM6hZd+Q7Dmr2IwaFEOQ92xG7HREREVHgap14CZzfv5eRWTsQfZKDmv2FQSHECE4HF1kjIiKioOEKV6G9x6Dm2AObRawmuDAohBh1dbGn25FLFsZuR0RERBTwvAY1VxUhvLlWxGqCB4NCiIkoPeq5bcrIY7cjIiIiCniW5CxYEjM87VgOavYLBoUQ0zMoGLMniVgJERERkZ8IgtdVBd2R7yB0cVDzYDEohBCZ2QB1fYWnzaBAREREwaJt4pyzg5ptZkQX7hW5osDHoBBCIkqPe27btdGwxaWIWA0RERGR/7gUKq+Vmtn9aPAYFEKIT7cjQRCxGiIiIiL/aplxuee2tqoIitYGEasJfAwKocLtQkTpMU/TmD1RxGKIiIiI/M+SlAlLfJqnHXWKayoMBoNCiFA2ViHMbAAAuCHAmDVB5IqIiIiI/EwQoB8709OMKmRQGAyZPzdmt9uxZs0a7NixAy0tLdBoNJg8eTKWL18OnU43oG0dPnwY69evR2lpKSwWC9RqNXJycrB48WJMnMhPwwcqouTs1QRLShacKq2I1RARERENDX3eDCRv+wwAoKktRZixDV0RMSJXFZj8dkXBbrfj2WefxerVq2Gz2TBjxgzodDps3boVjz/+OBoa+t9HbN26dfj973+PI0eOICUlBbNmzUJcXBwOHjyIZ599Fps2bfJX2SGj5/gEA2c7IiIioiBlTUhHZ3S8px1VdEDEagKb364ofP755ygqKsKYMWPw5JNPIjw8HED3m/53330XK1euxG9/+9sLbsdoNGLVqlWQyWR46qmnkJeX53ls9+7dePnll/Hee+9h/vz5nn1Q3yR2GzRVRZ42xycQERFR0BIEtOfNQOKu9QC6ux81z7xS5KICk1+uKDgcDmzcuBEAcNddd3m9gb/22muRkZGBwsJClJWVXXBbxcXFcDgcmDBhgldIAICCggKkp6ejs7MTNTU1/ig9JGgrTkLicgIAHAoVzKk5IldERERENHT0eTM8t7UVhZBaO0SsJnD5JSicOnUKZrMZCQkJyMzM9Hl81qxZAID9+/dfcFthYWH92qdGoxlYkSGs5/gEU9Z4QCIVsRoiIiKioWVOG40uTRQAQHC7EHn6kMgVBSa/BIXKykoA6DUkAEBWVpbX8/qSnZ0NlUqF48eP49SpU16P7dmzB1VVVcjNzUViYuIgqw4dEWU9p0Xl+AQiIiIKcoIE+txpnmY0Zz+6KH4Zo9DS0gIA553ZKCYmxut5fVGr1bj//vvxyiuv4Omnn0Zubi5iYmLQ1NSE0tJSTJkyBQ8++KA/yg4Jcn0zwlvrPW2OTyAiIqJQoM+bibjvV2eOKD0Kid0Gl5zjWwfCL0HBZrMBABQKRa+PnxmzcOZ5F1JQUACNRoOXX37Z66pCZGQkxo8fD622/1N7PvLIIz73yeVyPP/88wCA2NjYfm/Ln2QyGaRSKSIjI4d0PxEndnpu2+NToczIgXJI90h9kUq7u30N9XGnkYfHPnTx2IcuHnuRTS6A81M1pDYzJI4uJNaXwjxpzrDsWiqVQiaTIS4ublj2N1T8EhTcbvegHj/X2rVr8f7772PmzJlYtmwZ4uPj0dTUhI8++gjvv/8+iouL8eijjw6m5JChKjrouW3JnSpiJURERETDSBYGy7iZ0B7cCgDQHN05bEEhWPglKCiV3Z9Rd3Z29vr4mfv7M53pyZMn8d577yEzMxOPPPIIJJLuYRTp6el49NFH8etf/xp79uzBkSNHMHny5Atu789//nOfj7e0tAw4yPhDXFwcnE4nDAbD0O3E6UDm6cOeZnNaHoxDuT+6oDOfKg3pcacRicc+dPHYhy4ee/EJ2ZM9QUF5Yi+Mba1wS/263nCvIiMj4XA40NzcPOT76kkQBCQlJflte34ZzHym+05ra2uvj7e1tXk9ry/btm0D0D1T0pmQcIZEIkF+fj4A4MSJExddb6hQ15ZC2mkFALikMnRk5F3gFURERETBw5gzCS5Z94yask4LtOUnRa4osPglKGRkZAAAysvLe338zPoJZ57XlzOh4sxVinOdub+jg/PhXkhkydnVmDvSczmAh4iIiEKKSx7uNeNj1CnOfjQQfgkKeXl5UKlUaGxs7DUs7NmzBwAwbdo0n8fOdeYyXWlpaa+Pn7k/Pj6+18fprIjSs0GB06ISERFRKGrvsfha1KkDgMslYjWBxS9BQSaTYdGiRQCAt99+22t2o3Xr1qGyshJ5eXnIyTm7IvDGjRvx8MMPY9WqVV7bOtO1aPv27T4LtO3btw/bt2+HIAie51HvpBYTVHVnQ5sxh0GBiIiIQo8hdxrcQvdb3jCzAeqaYpErChx+G81xww034NixYygqKsLPf/5z5OXloaWlBcXFxdBqtT5rHxiNRtTV1aG9vd3r/pkzZ6KgoAC7d+/GCy+8gOzsbMTFxaG5udlzNeGWW25BcnKyv0oPShFlxyGge5B2lyYK1vg0kSsiIiIiGn5OpQamUWMRUd49vjXq1H6Y03NFriow+C0oyOVyPP300/j888+xfft27Nu3D2q1GgsWLMDy5cv7vV6BIAj4xS9+gS1btmDbtm2oqqpCRUUFVCoVpk6diquvvhpTpkzxV9lBS1N5dv0JY9Z4QBBErIaIiIhIPPq8GZ6gEH1qP2qvvJXvjfrBr/NDyeVyLF++HMuXL7/gc5ctW4Zly5b1+pggCFi4cCEWLlzoz/JCirbqbFDoSOdsR0RERBS69HkzkL7hHwAARXsTlI1VsCZeeJKdUOeXMQo0skitHVA21XjaJk6LSkRERCGsKyIG5pRsTzvq1P4+nk1nMCgEIU3Vac/tLlUEOnX+W3iDiIiIKBDpc6d7bkeUHRexksDBoBCENF7djsawDx4RERGFPGPWBM9tdW0pJHZbH88mgEEhKGmqijy3uRozEREREWBJyoRDoQIACC6n18Qv1DsGhSAjdHVC3WP9BA5kJiIiIgIgkcCUOc7TZPejC2NQCDLqmhIILicAwCkPhyUxXeSKiIiIiEYGU+Z4z23t99Ol0vkxKAQZbY/LaB1powGJVMRqiIiIiEaOnkFB1VgFmdkoYjUjH4NCkPEan8BVB4mIiIg8bLHJsGuiPG1eVegbg0IwcTqgrinxNDk+gYiIiKgHQYApi92P+otBIYioGioh7eoEALikMq+FRYiIiIjIu/tRBINCnxgUgkjPab4syVlwh8lFrIaIiIho5DFmnl1PQdHeBLm+WcRqRjYGhSCi7TE+wcTxCUREREQ+uiJ1sOmSPG1tGa8qnA+DQrBwu7jQGhEREVE/GL26H3E9hfNhUAgS4S11kFk7AABuCDCnjRa5IiIiIqKRyWc9BbdbxGpGLgaFINFzfII1IQ3OcLWI1RARERGNXKbMcXBDAACEmY0Ib6oWuaKRiUEhSHivn8BuR0RERETn41RqYEnK8LQ5+1HvGBSChLayx0Bmjk8gIiIi6pNP9yPywaAQBOT6ZsiNrZ52RwZnPCIiIiLqi6nHNKnaikLA5RSxmpGJQSEI9Ox2ZItJgKPH0uRERERE5KsjfQxcUhkAQGq3QV1bJnJFIw+DQhDQVHJ8AhEREdFAuOThMKeenSVSy2lSfTAoBAFN1dkZjzq40BoRERFRv3ivp8BxCudiUAhwMrMRypY6T5vjE4iIiIj6x5R1Niioq4shdHWKWM3Iw6AQ4HqOT7BrotAZnSBiNURERESBw5ycBac8HAAgcTq83lcRg0LA81o/ISMXEAQRqyEiIiIKIFKZ17TyEWXsftQTg0KA40JrRERERBfPa5pUjlPwwqAQwASHHcqGSk+7I210H88mIiIionP1HKegqq+ApNMiXjEjDINCAFM1VELy/eIgLlkYrPFpIldEREREFFis8alwKFQAAAFuqGtKRa5o5GBQCGA9f5AtSZnA94uGEBEREVE/CRKY03I8TU31aRGLGVkYFAKYurbEc9ucki1iJURERESBqyNtjOe2uqZYxEpGFgaFAKauPXtFwZya08cziYiIiOh8zD2CgqamBHC5RKxm5GBQCFBSiwmK9iZPm1cUiIiIiC6OOSUb7u+nmJd2WhHeXCNyRSMDg0KA6nk1oUsdAXtkrIjVEBEREQUulzwc1oQMT1tTze5HAINCwFLX9ByfkMOF1oiIiIgGoec08xqOUwDAoBCwvMYnsNsRERER0aD0DApqznwEgEEhMLndHMhMRERE5Ec9BzSHtzVCZjaIWM3IwKAQgBRtDZDZzJ62OTlLxGqIiIiIAp89MhZ2TZSnreY4BQaFQNTzaoI1NhmucJWI1RAREREFAUGAOb3HNKkMCgwKgajniswcn0BERETkHx2pPYMCxykwKAQg7xWZOT6BiIiIyB96DmhW1ZVDcHSJWI34GBQCjOCwQ9lQ6WmbU3lFgYiIiMgfrEmj4JKGAQAkzi6oGirELUhkDAoBRtVQCYnLCQBwycJgjU8TuSIiIiKi4OCWymBOOTtJTKgPaGZQCDA9xydYkkYBUpl4xRAREREFGXNqj4XXQnycAoNCgFHV9RzIzPEJRERERP7ktUJzdTHgdotYjbgYFAIMZzwiIiIiGjrmHkEhrEMPuaFFxGrExaAQQKQWE8LbGz1trshMRERE5F8OdSRsMQmetjqEux8xKASQngutdakjYI+MFbEaIiIiouDUkcaF1wAGhYCirum5fkI2IAgiVkNEREQUnMwMCgAYFAJKzysKHJ9ARERENDR6DmhWNlZCYreJWI14GBQChdt9TlDg+AQiIiKioWCLS4FToQQACG63V6+OUMKgECAUbY2Q2cyetqXHYiBERERE5EeCBB091lMI1YXXGBQChLr2bJK1xibDGa4WsRoiIiKi4NZzmlRNTWjOfMSgECC8VmTm+AQiIiKiIdVz5iN1dQngdolYjTgYFAJEzysKHMhMRERENLTMKdlwfz/DpKzTgvDmOpErGn4MCgFAcNihbKj0tDmQmYiIiGhouRRKWBPSPW11TeiNU2BQCADKhipIXE4AgEsWBktCmsgVEREREQU/c/LZXhzq2jIRKxEHg0IAUNeXe25bEjMAqUzEaoiIiIhCg7nHLJPqOgYFGoFUdT2CQlKmiJUQERERhQ5L8tmgoGyqhuCwi1jN8GNQCACqegYFIiIiouFmjU+FSxYGABBcTigbqkSuaHgxKIxwQpcdyuZaT9uSzKBARERENCwkUlgSR3maodb9iEFhhFM2VkHoMZDZGpcickVEREREocOcHLrjFPw2KtZut2PNmjXYsWMHWlpaoNFoMHnyZCxfvhw6nW7A22toaMCaNWtw7Ngx6PV6hIeHIykpCfn5+Vi6dKm/yh7xfAYyS6QiVkNEREQUWiw9BjSrQiwo+OWKgt1ux7PPPovVq1fDZrNhxowZ0Ol02Lp1Kx5//HE0NDQMaHt79+7Fo48+ii1btkCj0WDmzJnIzMxEU1MTvvrqK3+UHDA4kJmIiIhIPD2vKIQ310HSaRWxmuHllysKn3/+OYqKijBmzBg8+eSTCA8PBwCsW7cO7777LlauXGLIj24AACAASURBVInf/va3/dpWRUUF/vd//xdKpRK//OUvkZeX53nM5XKhvLy8j1cHHw5kJiIiIhJPpy4RToUS0k4rBLihqq9Ax6ixYpc1LAZ9RcHhcGDjxo0AgLvuussTEgDg2muvRUZGBgoLC1FW1r9LNe+88w4cDgcefPBBr5AAABKJBNnZ2ed5ZfARHBzITERERCQqQQJzjw9rQ2mcwqCDwqlTp2A2m5GQkIDMTN83srNmzQIA7N+//4LbqqmpQWFhIZKSkjB9+vTBlhbwlI3VHMhMREREJLKe6ymE0jiFQXc9qqysBIBeQwIAZGVleT2vL8ePHwcATJo0CXa7HTt37vRcicjIyMDs2bOhUqkGW3LA6JlYLQnpHMhMREREJIJQXaF50EGhpaUFAM47s1FMTIzX8/pSXV0NAJDL5fjVr36Furo6r8dXrVqFRx99FOPGjRtMyQHDayBzjyRLRERERMOn5/swRXsTpBYTnCqtiBUNj0EHBZvNBgBQKBS9Pn5mzMKZ5/XFbDYDANavXw+1Wo3HHnsMEyZMgF6vx+rVq7F9+3a8+OKL+POf/4zo6Oh+1ffII4/43CeXy/H8888DAGJjY/u1HX+TyWSQSqWIjIw873O0TWdX/3Nnj+vzuRQYpNLuq0I8lqGHxz508diHLh77IBIRAac6AlKzEQAQb2iEJSn1vE+XSqWQyWSIi4sbrgqHxKDHKLjd7kE93pPL5QLw/9u79+CoyvuP459lc90lJCRACD9rCKGYsWCQqxZG6tRRUKbDRQg646WDdhQvXLTVFkekSNE6g22nM4wtoghFiyiMMjQ6rUIJjuHeIrkYEhKDJCEbSBYWk7DJ/v6wG7OeBJKcE3Y3+379lT3n8Ox350vgfPac5zlSS0uLnnjiCU2aNEkOh0PDhg3Tk08+qczMTHk8Hn300Uemag4HtkvNimn3mPCma0YGsRoAAIAIZrOp8Qc/bHsZW3kiiMVcPaavKMTHx0uSmpqaOtzv395+NaTO+I9JTk5Wdna2Yf+tt96q0tJSHT9+vMv1rV279rL7XS5Xt8KMVQYPHqyWlhY1NDR0uN/xdWnAROba+ESpk2MRPvzfKnXWd/Rd9D5y0fvIRe/7FseQa+UsOiRJspcVXLaviYmJ8nq9qq2tvVrlSZJsNpvS0tIsG8/0FQX/rTt1dXUd7j979mzAcZczZMgQSer0Mo1/u9vt7nad4SZgfgITmQEAAIIqElc+Mh0U0tPTJanTB6G1X7XoSoYPHy5JunDhQof7z58/L6lrVyfCnZMHrQEAAISM9isfxZw/p+jz54JYzdVhOihkZWXJ4XCopqamw7CQn58vSRo3btwVxxozZoxiY2NVXV3d4SpJBQUFkjpfirUvaZ9UWfEIAAAguLz9k9Q8ILntdSRcVTAdFKKiojR9+nRJ0oYNGwJWN9q5c6cqKiqUlZWlkSO/m4ybm5urJUuWaMuWLQFjxcbGasaMGWppadH69esDxjp69Kj27Nkjm82m2267zWzZIY0nMgMAAIQeT7svb51f9/2gYHoysyTNmTNHx44dU3FxsRYvXqysrCy5XC6VlJQoISFBixYtCjje7Xbr9OnTOnfOeMnm7rvvVmFhoQ4fPqzFixdr5MiRcrvd+vLLL+Xz+bRgwYKA0NEX8URmAACA0HNx2AgNLDooKTKuKFgSFGJiYrRixQpt375deXl5OnDggJxOp6ZNm6acnJxuPavAP9YHH3ygvXv36ujRo4qOjtbo0aN11113dekWpnDHRGYAAIDQE3BF4XSZ5PNJNlsQK+pdlgQF6dsT/JycHOXk5Fzx2Pnz52v+/PmdFxUVpTlz5mjOnDlWlRdWmMgMAAAQetrPG4365oJi6mvVPHBIECvqXabnKMB6jvZBgfkJAAAAIaEl3qnG5NS2184+fvsRQSHE2LzNij9zqu01VxQAAABCR8DzFPr4hGaCQohhIjMAAEDoCpinUEVQwFVkmMhst2waCQAAAEwKfELzScnXGsRqehdBIcQwkRkAACB0XUwbLt//VjqyNzcqzlUV5Ip6D0EhxDCRGQAAIHS1xsSpcdB3t4b35ecpEBRCCBOZAQAAQp/n/773PIU+iqAQQgImMtuZyAwAABCKAucpEBRwFQTcdjSUicwAAAChqP3KR47qr6T/fdHb1xAUQkjAikfcdgQAABCSvhnyA/ls355G9/M299kJzQSFEOKoLm/7maAAAAAQmnzRMQG3iLe/K6QvISiECFuLV/E1lW2vL6YND14xAAAAuKz2q1M6qsqDVkdvIiiEiLjaU+rnn8jcz67GIdcEuSIAAAB05uLQ4W0/c0UBvap9Ev0m9QfyMZEZAAAgZLW/+8NRXdEnn9BMUAgR7YNC+4QKAACA0PPN0HT59N0TmmPrqoNckfUICiEi4IoC8xMAAABC2rdPaE5re90Xbz8iKISC1lbF13zV9pKJzAAAAKGv/SqVfXFCM0EhBMTVVcl+qUmS5LPZdDH12iBXBAAAgCshKKDXtf+L1Tjo/+SLjg1aLQAAAOiagAnNVeWSzxesUnoFQSEEtL+njduOAAAAwsPFoeltP0c1XVTMuTNBrMZ6BIUQEPBEZlY8AgAACAutcQ41Jg9te93Xbj8iKASbz6f4qoq2l1xRAAAACB/tz92cfWzlI4JCkMXU1yqq6WLb6/aXsAAAABDaAic0ExRgofZ/oRqTU9Ua5whiNQAAAOiOvjyhmaAQZAFPZOa2IwAAgLDS/vwt6psLimlwBa8YixEUgiwgKDCRGQAAIKy0xPdXU9Lgttd9aUIzQSGYfD6uKAAAAIS5wNuP+s48BYJCENkb6hR90d32+huuKAAAAISdvvqEZoJCEMV+Xdr2c/OAFHmdA4JYDQAAAHrCsPJRH5nQTFAIothT3wUFbjsCAAAITxfTvlvePtrjlr2hLojVWIegEESxX5e1/czzEwAAAMKT15mo5gHJba/b3zUSzggKQRR4RSHjMkcCAAAglLU/l4urPBHESqxDUAiS1oZziq6vbXvNrUcAAADhq/25HFcUYMqlsi+/+9k5QJcSBgaxGgAAAJjhaXdFof1dI+GMoBAkl8qK236+OHS4ZLMFrxgAAACYEvCE5oY6tdSfDV4xFiEoBIm33RUFbjsCAAAIb96EgbrUP+m716XFlzk6PBAUgiTgigJBAQAAIOx52p3TtT/XC1cEhSDwXfSopepU22tWPAIAAAh/7c/pLpUWBbESaxAUgqHyZNuP3liHmpMGB7EYAAAAWKH9XSKXuPUIPeGr/G4m/Ddp6UxkBgAA6APaX1Fora2W74I7iNWYR1AIhor2T2QeHrw6AAAAYJlLA5J1yZHw3YavwnuZ1KhgFxCJfF/xRGYAAIA+x2bTxbThcpytUf/rrtelmLhgV2QKQeEq8zU1SQETmYcHrxgAAABYqvSepzUgOUVTR6SotrY22OWYwq1HV1v1KcnXKklqjYlVY0pakAsCAACAVXz2vvM9fN/5JGHClp6pfn/cokT3ORUWnZT6kdUAAAAQeggKQWBz9FdseoYuONOlhoZglwMAAAAY8HU2AAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAgyirBmpubtaOHTu0b98+uVwu9e/fX9nZ2crJyVFKSkqPx62qqtLTTz+tS5cuKTs7W8uXL7eqZAAAAACdsOSKQnNzs1atWqVt27apsbFREyZMUEpKinbv3q1nnnlG1dXVPR77L3/5i7xerxVlAgAAAOgiS4LC9u3bVVxcrFGjRumPf/yjli5dqt/97ne6//775Xa7tW7duh6N+8knn+j48eP66U9/akWZAAAAALrIdFDwer3Kzc2VJC1cuFBxcXFt+2bOnKn09HQVFhaqrKysW+M2NDRo06ZNGjNmjKZMmWK2TAAAAADdYDooFBUVyePxKDU1VRkZGYb9kydPliQdPHiwW+O+8cYbam5u1sMPP2y2RAAAAADdZDooVFRUSFKHIUGSRowYEXBcVxw+fFifffaZZs+eraFDh5otEQAAAEA3mQ4KLpdLkjpd2Sg5OTnguCtpbGzU66+/rmHDhmnWrFlmywMAAADQA6aXR21sbJQkxcbGdrjfP2fBf9yVvPPOO6qtrdXzzz+vqCjzq7cuW7bMsC0mJkYvvfSSJGnQoEGm36MnoqKiZLfblZiYGJT3R3DY7XZJou8RiN5HLnofueh95LLb7YqKitLgwYODXYoppq8o+Hw+U/vbKy0tVW5urm655RaNHj3abGkAAAAAesj0V/bx8fGSpKampg73+7e3Xw2pIy0tLXrttdfkcDh0//33my2rzdq1ay+73+VydSvMWGXw4MFqaWlRQ0PDVX9vBI//WyX6HnnofeSi95GL3keuxMREeb1e1dbWXtX3tdlsSktLs2w800HBf+tOXV1dh/vPnj0bcFxn6urqVF5erqSkJMPJvcfjkSSdOHFCL7zwguLi4vTss8+aLR0AAABAJ0wHhfT0dEnSyZMnO9zvf36C/7grqa+vV319fYf7PB6PCgoK5HA4elApAAAAgK4yHRSysrLkcDhUU1OjkydPGpZJzc/PlySNGzfusuMMGTJEW7du7XDf8ePHtXLlSmVnZ2v58uVmSwYAAABwBaYnM0dFRWn69OmSpA0bNgSsbrRz505VVFQoKytLI0eObNuem5urJUuWaMuWLWbfHgAAAEAvML/+qKQ5c+bo2LFjKi4u1uLFi5WVlSWXy6WSkhIlJCRo0aJFAce73W6dPn1a586ds+LtAQAAAFjM9BUF6dvnEqxYsUJz585VTEyMDhw4oDNnzmjatGl6+eWXeboyAAAAEGZsvmCsDRpCqqqqgrY8al5ZHUumRRiWyotc9D5y0fvIRe8jV2JioqaOSAn75VEtuaIAAAAAoG8hKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADAgKAAAAAAwICgAAAAAMIiyaqDm5mbt2LFD+/btk8vlUv/+/ZWdna2cnBylpKR0aQyPx6MjR47o0KFDKi8vl8vlks1m0zXXXKOpU6fq9ttvV1SUZSUDAAAA6IQlZ93Nzc1atWqViouLNXDgQE2YMEG1tbXavXu3Dh8+rBdffFFDhw694jgffvih3n//fdlsNmVkZGj8+PFyu90qLi7WiRMn9Pnnn2v58uWKjY21omwAAAAAnbAkKGzfvl3FxcUaNWqUnnvuOcXFxUmSdu7cqbfeekvr1q3TypUrrzhOXFycZs+erTvuuEPJyclt26uqqrRq1SoVFRXpvffe07333mtF2QAAAAA6YXqOgtfrVW5uriRp4cKFbSFBkmbOnKn09HQVFhaqrKzsimPNmjVL99xzT0BIkKS0tLS2cLBv3z6zJQMAAAC4AtNBoaioSB6PR6mpqcrIyDDsnzx5siTp4MGDpt5n+PDhkqRz586ZGgcAAADAlZkOChUVFZLUYUiQpBEjRgQc11M1NTWSpKSkJFPjAAAAALgy00HB5XJJUqcrG/lvI/If11O7du2SJE2YMMHUOAAAAACuzPRk5sbGRknqdCUi/5wF/3E98fHHH+vYsWNyOp2aNWtWt/7ssmXLDNtiYmL00ksvSZIGDRrU47rMiIqKkt1uV2JiYlDeH8Fht9slib5HIHofueh95KL3kctutysqKkqDBw8OdimmmL6i4PP5TO2/koKCAr355puy2Wx69NFHDROdAQAAAFjP9BWF+Ph4SVJTU1OH+/3b26+G1FUVFRV65ZVX5PV69fOf/1yTJk3q9hhr16697H6Xy2U6zPTE4MGD1dLSooaGhqv+3gge/7dK9D3y0PvIRe8jF72PXImJifJ6vaqtrb2q72uz2ZSWlmbZeKavKPhv3amrq+tw/9mzZwOO66rq6mqtXr1aHo9H8+bN04wZM8wVCgAAAKDLTAeF9PR0SdLJkyc73O9/foL/uK44e/asXnzxRdXX1+vOO+/UvHnzzJYJAAAAoBtMB4WsrCw5HA7V1NR0GBby8/MlSePGjevSeBcuXNDq1at15swZ/eQnP9EDDzxgtkQAAAAA3WQ6KERFRWn69OmSpA0bNgSsbrRz505VVFQoKytLI0eObNuem5urJUuWaMuWLQFjNTU1ac2aNaqsrNTNN9+sRx55RDabzWyJAAAAALrJ9GRmSZozZ46OHTum4uJiLV68WFlZWXK5XCopKVFCQoIWLVoUcLzb7dbp06cNT1l+++23VVJSon79+slut2vdunUdvt9jjz1mRdkAAAAAOmFJUIiJidGKFSu0fft25eXl6cCBA3I6nZo2bZpycnK6PJHZ4/FIklpbW5WXl9fpcQQFAAAAoHfZfMFYGzSEVFVVBW151LyyOpZMizAslRe56H3koveRi95HrsTERE0dkcLyqAAAAAD6HoICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAgyirBmpubtaOHTu0b98+uVwu9e/fX9nZ2crJyVFKSkq3xvJ4PHr33Xe1f/9+1dfXKykpSRMnTtT8+fPldDqtKhkAAABAJyy5otDc3KxVq1Zp27Ztamxs1IQJE5SSkqLdu3frmWeeUXV1dZfHOn/+vH7zm99o165dstvtmjhxouLj4/WPf/xDv/71r3X+/HkrSgYAAABwGZZcUdi+fbuKi4s1atQoPffcc4qLi5Mk7dy5U2+99ZbWrVunlStXdmmsjRs3qqqqSpMmTdLSpUtlt9slSRs2bFBubq42btyoxx9/3IqyAQAAAHTC9BUFr9er3NxcSdLChQvbQoIkzZw5U+np6SosLFRZWdkVx6qvr9fevXtlt9v10EMPtYUESbrvvvs0YMAA5eXlqb6+3mzZAAAAAC7DdFAoKiqSx+NRamqqMjIyDPsnT54sSTp48OAVxzpy5Ih8Pp+uv/56JSUlBeyLjo7W+PHj1draqqNHj5otGwAAAMBlmL71qKKiQpI6DAmSNGLEiIDjzIyVkZGhTz/9VOXl5T2otGM2m82ysbqrn80me7/gvT+uvn7/+/tG3yMPvY9c9D5y0fvI5e/91T7PtPr9TAcFl8slSZ2ubJScnBxwXFfG8v+Z7/O/R1fG8lu2bJlhm9Pp1KpVqyRJQ4cO7fJYVps2KlVSatDeH8FE3yMXvY9c9D5y0ftIlZaWFuwSTDF961FjY6MkKTY2tsP9/jkL/uPMjOXf3tTU1O06AQAAAHSd6SsKPp/P1P6OjrXyssnatWstG8tKzz77rCTppZdeCnIluJroe+Si95GL3kcueh+5+krvTQeF+Ph4SZ1/y+/f3n41pCuN1dnVB/9YnV1xCCfNzc3BLgFBQN8jF72PXPQ+ctH7yNVXem/61qNBgwZJkurq6jrcf/bs2YDjujKW/898n/89ujIWAAAAgJ4zHRTS09MlSSdPnuxwv//5Cf7jzIzl396VsQAAAAD0nOmgkJWVJYfDoZqamg5P8PPz8yVJ48aNu+JYY8eOlc1mU2FhoRoaGgL2Xbp0SYcOHZLNZtONN95otmwAAAAAl2E6KERFRWn69OmSpA0bNgTML9i5c6cqKiqUlZWlkSNHtm3Pzc3VkiVLtGXLloCxBg4cqClTpsjr9Wr9+vVqaWlp27d582a53W5NnTrV8DA2AAAAANay+bqzLFEnmpubtXLlSpWUlGjgwIHKysqSy+VSSUmJEhIStHr16oDnFWzdulXbtm3TtGnT9NhjjwWM5Xa7tXz5ctXU1Cg1NVWZmZmqrKxUZWWlUlNTtXr1ag0YMMBsyQAAAAAuw/QVBUmKiYnRihUrNHfuXMXExOjAgQM6c+aMpk2bppdffrlbDzUbMGCA1qxZo+nTp8vr9Wr//v26ePGipk+frjVr1hASAAAAgKvAkisKAAAAAPoWS64oAAAAAOhbCAoAAAAADAgKAAAAAAwICgAAAAAMCAoAAAAADKKCXUBf0NzcrB07dmjfvn1yuVzq37+/srOzlZOTo5SUlG6N5fF49O6772r//v2qr69XUlKSJk6cqPnz58vpdPbSJ0BPWdF7j8ejI0eO6NChQyovL5fL5ZLNZtM111yjqVOn6vbbb1dUFL+qocbK3/v2qqqq9PTTT+vSpUvKzs7W8uXLLawaVrC699XV1dqxY4eOHTum+vp6xcXFKS0tTZMmTdLPfvazXvgE6Ckre3/06FHt2rVLpaWlunjxopxOp0aOHKm77rpLY8aM6aVPgO4qKyvTf//7X504cUIlJSU6d+6coqOj9be//a1H44XbeZ79hRdeeCHYRYSz5uZmrVq1Snv37lV0dLRuuOEG+Xw+HT58WHv37tWkSZPUv3//Lo11/vx5LV++XEeOHFFCQoJGjx6tb775RocOHVJ+fr6mTp2q2NjYXv5E6Cqrev/+++/rzTff1KlTpzRo0CCNGjVKTqdTpaWlOnz4sI4fP64f//jHhIUQYuXv/fetXbtWNTU1kqShQ4fqlltusbJ0mGR17/fv36/f/va3Kisr05AhQ3TdddfJ6XTq9OnTqqio0J133tmLnwbdYWXvd+7cqT//+c+qqalRRkaGfvjDH8rn86mgoED//ve/lZiYqMzMzF7+ROiKv/71r/roo4/09ddfq7GxUZJkt9s1d+7cbo8Vjud5nHmYtH37dhUXF2vUqFF67rnnFBcXJ+nbfwTeeustrVu3TitXruzSWBs3blRVVZUmTZqkpUuXym63S5I2bNig3Nxcbdy4UY8//nivfRZ0j1W9j4uL0+zZs3XHHXcoOTm5bXtVVZVWrVqloqIivffee7r33nt77bOge6z8vW/vk08+0fHjx3Xbbbfpn//8p9VlwwJW9r68vFx/+MMfFB8fr1/+8pfKyspq29fa2qqTJ0/2ymdAz1jVe7fbrS1btigqKkrPP/98QN8///xzvfrqq9q0aZNuueWWtvdA8IwaNUrDhw9XZmamMjMz9Ytf/KLHY4XjeR5zFEzwer3Kzc2VJC1cuDDgF3rmzJlKT09XYWGhysrKrjhWfX299u7dK7vdroceeqjtL48k3XfffRowYIDy8vJUX19v/QdBt1nZ+1mzZumee+4JCAmSlJaW1hYO9u3bZ2H1MMPK3rfX0NCgTZs2acyYMZoyZYqlNcMaVvf+jTfekNfr1aJFiwJOFiWpX79+fKMcQqzsfUlJibxer0aPHm3o+0033aRrr71WTU1NOnXqlLUfAj0ya9YszZ8/X+PHj1dSUlKPxwnX8zyCgglFRUXyeDxKTU1VRkaGYf/kyZMlSQcPHrziWEeOHJHP59P1119v+IsYHR2t8ePHq7W1VUePHrWmeJhiZe8vZ/jw4ZKkc+fOmRoH1umt3r/xxhtqbm7Www8/bEmdsJ6VvT916pQKCwuVlpam8ePHW14rrGVl76Ojo7v0nj29fRGhKVzP8wgKJlRUVEhSh/9oSNKIESMCjjMzln97eXl5d8tEL7Cy95fjv1fdzLcYsFZv9P7w4cP67LPPNHv2bA0dOtR8kegVVvb+iy++kCTdcMMNam5u1u7du7VhwwZt2LBB//rXv3Tx4kWLqoYVrOx9ZmamHA6HvvjiCxUVFQXsy8/P11dffaXrrruOfwv6mHA9z2OOggkul0uSOl3pwH8rif+4roz1/dtP/Pzv0ZWx0Pus7P3l7Nq1S5I0YcIEU+PAOlb3vrGxUa+//rqGDRumWbNmWVMkeoWVva+srJQkxcTE6Fe/+pVOnz4dsH/Lli166qmndP3115spGRaxsvdOp1OPPPKI/vSnP2nFihW67rrrlJycrDNnzqi0tFRjx47VokWLrCseISFcz/MICib4Z793NkPdfw+j/zgzY/m3NzU1dbtOWM/K3nfm448/1rFjx+R0OjmBDCFW9/6dd95RbW2tnn/+eVa2YhwGFwAABi9JREFUCnFW9t7j8Uj69ssAp9Opp59+WqNHj1Z9fb22bdumvLw8vfLKK1q7dq0GDhxo0SdAT1n9e3/TTTepf//+evXVVwOuKiQmJupHP/qREhISTFaMUBOu53n8r2SCz+cztb+jY202m6macHVY2fuOFBQU6M0335TNZtOjjz7a6TcQuPqs7H1paalyc3N1yy23aPTo0WZLQy+zsvetra2SpJaWFj3xxBPKzs6WJDkcDj355JOqqqpSaWmpPvroIy1YsKDnRcMSVv+b/+GHH2rz5s1t6+cPGTJEZ86c0d///ndt3rxZJSUleuqpp8yUjBATrud5zFEwIT4+XlLn6c+/vSvLm/nH6uzbCP9Yoba+bqSysvffV1FRoVdeeUVer1cPPvigJk2a1PNCYTmret/S0qLXXntNDodD999/v7VFoldY+XvvPyY5ObktJLR36623SpKOHz/eo1phLSt7X1BQoE2bNmn48OFatmyZrr32WsXFxenaa6/VU089pYyMDOXn5+s///mPdR8AQReu53lcUTBh0KBBkqS6uroO9589ezbguK6M5f8z3+d/j66Mhd5nZe/bq66u1urVq+XxeDRv3jzNmDHDXKGwnFW9r6urU3l5uZKSkrR27dqAff7bUk6cOKEXXnhBcXFxevbZZ82WDpOs/L0fMmSIJGnw4MEd7vdvd7vd3a4T1rOy93v27JH07UpJ/foFfl/br18/TZo0SSdPntTx48c7DJEIT+F6nkdQMCE9PV2SOn0ojn89Zf9xZsbyb+/KWOh9Vvbe7+zZs3rxxRdVX1+vO++8U/PmzTNfKCxnde/r6+s7XTfb4/GooKBADoejB5XCalb23r/08YULFzrcf/78eUk9uyoJ61nZe/+Jov8b5u/zb+/s7wbCU7ie5xEUTMjKypLD4VBNTY1OnjxpWPIqPz9fkjRu3LgrjjV27FjZbDYVFhaqoaFBiYmJbfsuXbqkQ4cOyWaz6cYbb7T2Q6BHrOy99O1/CKtXr9aZM2f0k5/8RA888IDlNcMaVvV+yJAh2rp1a4f7jh8/rpUrVyo7O1vLly+3pnCYZuXv/ZgxYxQbG6vq6mq5XC7Dt4gFBQWSOl9KEVeXlb33//9eWlra4X7/dv9VJ/QN4XqexxwFE6KiojR9+nRJ3z5+u/19Zzt37lRFRYWysrI0cuTItu25ublasmSJtmzZEjDWwIEDNWXKFHm9Xq1fv14tLS1t+zZv3iy3262pU6eynn6IsLL3TU1NWrNmjSorK3XzzTfrkUceCbvJTpHEyt4jvFjZ+9jYWM2YMUMtLS1av359wFhHjx7Vnj17ZLPZdNttt/Xyp0JXWNl7/7yzvLw8wwPaDhw4oLy8PNlsNuanham+dp7HFQWT5syZo2PHjqm4uFiLFy9WVlaWXC6XSkpKlJCQYFgL2e126/Tp0x0+affBBx9USUmJ8vPztWTJEmVmZqqyslKVlZVKTU3lW+YQY1Xv3377bZWUlKhfv36y2+1at25dh+/32GOP9dpnQfdY+XuP8GJl7++++24VFhbq8OHDWrx4sUaOHCm3260vv/xSPp9PCxYsCDjxRHBZ1fuJEyfqpptu0ueff67f//73yszM1ODBg1VbW9t2NWHBggUaNmzYVfts6Nzhw4f13nvvBWzzer0BV3vnzp3bdjWpr53nERRMiomJ0YoVK7R9+3bl5eXpwIEDcjqdmjZtmnJycro1KWXAgAFas2aNtm7dqgMHDmj//v1KTEzU9OnTNX/+fB7nHmKs6r1/4mpra6vy8vI6PY6gEDqs/L1HeLGy9/6xPvjgA+3du1dHjx5VdHS0Ro8erbvuuqvLty7i6rCq9zabTUuXLtWnn36qPXv26KuvvlJ5ebkcDoduvPFGzZgxQ2PHju3lT4OucrvdKikpCdjm8/kCtnV10YFwPM+z+cwu+A4AAACgz2GOAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAAOCAgAAAAADggIAAAAAA4ICAAAAAIP/B8yyQIg4qhlbAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 900x600 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(x,y)\n",
"plt.plot(x,y_lue, c='k')\n",
"plt.hist(data, 20, (0,1), density=True, alpha=0.3)\n",
";"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.5"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment