{ "cells": [ { "cell_type": "markdown", "id": "21d6145a", "metadata": {}, "source": [ "# GBM Example Notebook\n" ] }, { "cell_type": "markdown", "id": "0d771f71-4337-4215-bed9-a2f1a79e2b94", "metadata": {}, "source": [ "An example use in the Python Interface is:" ] }, { "cell_type": "code", "execution_count": 1, "id": "66fc821f-2c86-4748-87c6-dd837757bbeb", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHHCAYAAABXx+fLAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAkUBJREFUeJzt3Xl4U1X6B/DvTdp03/fS0pa1gFDKYgVFAUEGkXEZRdERUNyRGWXcGOcnOs64jY44ijrqII4zuO8D7oLIoqxl37tQSlu6722a5Pz+SO5t0iZt0mZp0u/nefpob25uTtLQvH3P+54jCSEEiIiIiPoRlacHQERERORuDICIiIio32EARERERP0OAyAiIiLqdxgAERERUb/DAIiIiIj6HQZARERE1O8wACIiIqJ+hwEQERER9TsMgMjp0tPTsWjRIo889qOPPgpJkjzy2LYUFBRAkiSsWbPG00PxenwtuyZJEh599FFPD8Ot+uK/eWeS3/PPPvusp4ficxgAkd3279+Pq6++GmlpaQgMDMSAAQMwc+ZMvPjii54emlO8/PLLfeKD9ezZs3jooYcwevRohIaGIjAwEEOGDMFNN92EzZs3e3p4TrN+/XqPfljLHyzmX+Hh4Rg7dixeeukl6PV6j43N28hBiEqlQlFRUafb6+rqEBQUBEmScPfddzt8/aamJjz66KPYuHGjE0brPObvHZVKheTkZFxyySU9Gqen/z30R36eHgB5h61bt2LatGkYOHAgbr31ViQmJqKoqAg///wzXnjhBSxdulQ59+jRo1CpvC+2fvnllxEbG+ux7BUAbN++HXPmzEF9fT2uu+463HHHHQgICEB+fj4+/fRTrFmzBj/++CMuvPBCj43RWdavX49Vq1Y59Es/LS0Nzc3N8Pf3d9o45s+fj0svvRQAUFtbi/Xr12Pp0qUoLCzE3/72N6c9jjs0NzfDz89zv9YDAgLwzjvv4IEHHrA4/vHHH/fquk1NTXjssccAAFOnTrW47U9/+hMeeuihXl2/N2bOnIkFCxZACIH8/Hy8/PLLmD59OtatW4fZs2fbfZ2e/Hug3mEARHb561//ioiICOzYsQORkZEWt509e9bi+4CAADeOzHdUV1fjiiuugJ+fH3Jzc5GZmWlx+1/+8he8++67CAoK8tAIPUen08FgMECj0SAwMNCp1x43bhx++9vfKt/fddddyMnJwdq1a7sMgIQQaGlp6VM/D2e/No669NJLrQZAa9euxZw5c/DRRx85/TH9/Pw8GvQNGzbM4v1z5ZVXYsyYMVi5cqVDARC5n/f9mU4ecfLkSYwaNapT8AMA8fHxFt93rAFas2YNJEnC5s2b8bvf/Q5xcXGIjIzE7bffDq1Wi5qaGixYsABRUVGIiorCAw88ACGEcv+NGzdCkqROaWV760HefPNNTJ8+HfHx8QgICMDIkSPxyiuvdBrzwYMH8eOPPyopbfO/NGtqanDPPfcgNTUVAQEBGDJkCJ5++mkYDAaL69TU1GDRokWIiIhAZGQkFi5ciJqami7HJ3v11VdRUlKClStXdgp+AGO6ff78+Zg4caLF8eLiYtx8881ISEhAQEAARo0ahdWrV3e6/9mzZ7F48WIkJCQgMDAQWVlZeOuttyzOMa83WLVqFQYNGoTg4GBccsklKCoqghACjz/+OFJSUhAUFITLL78cVVVVnR7ryy+/xJQpUxASEoKwsDDMmTMHBw8eVG5ftGgRVq1apTwv+avjGFauXInBgwcjICAAhw4dsvkzP3LkCObNm4e4uDgEBQVh+PDhePjhh+163a29zgkJCZ0+VNPT03HZZZfh66+/xoQJExAUFIR//vOfAIC8vDxcc801iI6ORnBwMM477zysW7dOua8QArGxsVi2bJlyzGAwIDIyEmq12uI98vTTT8PPzw8NDQ3KaxUaGori4mJcccUVCA0NRVxcHO67775O03Qda4AKCwtx1113Yfjw4QgKCkJMTAyuueYaFBQUWNxP/je6ZcsWLFu2DHFxcQgJCcGVV16J8vJyu1+766+/Hrm5uThy5IhyrLS0FD/88AOuv/56q/fp7n1ZUFCAuLg4AMBjjz2mvFfk52mtBkin0+Hxxx9X3jvp6en44x//iNbWVovz5J/p5s2bce655yIwMBCDBg3Cv//9b7ufc0ejR49GbGws8vPzAQA//fQTrrnmGgwcOBABAQFITU3Fvffei+bmZuU+Xf17MPfaa68pz2nixInYsWNHj8dJzACRndLS0rBt2zYcOHAA55xzTo+usXTpUiQmJuKxxx7Dzz//jNdeew2RkZHYunUrBg4ciCeeeALr16/H3/72N5xzzjlYsGCBU8b+yiuvYNSoUfj1r38NPz8/fPHFF7jrrrtgMBiwZMkSAMDKlSuxdOlShIaGKh+cCQkJAIzp94suugjFxcW4/fbbMXDgQGzduhXLly9XAhbA+CF3+eWXY/PmzbjjjjswYsQIfPLJJ1i4cKFd4/ziiy8QFBSEq666yu7nVlZWhvPOO0+prYiLi8OXX36JxYsXo66uDvfccw8A49TI1KlTceLECdx9993IyMjABx98gEWLFqGmpga///3vLa773//+F1qtFkuXLkVVVRWeeeYZzJs3D9OnT8fGjRvx4IMP4sSJE3jxxRdx3333WQRcb7/9NhYuXIhZs2bh6aefRlNTE1555RVccMEF2LNnD9LT03H77bfjzJkz+Pbbb/H2229bfW5vvvkmWlpacNtttyEgIADR0dGdAk4A2LdvH6ZMmQJ/f3/cdtttSE9Px8mTJ/HFF1/gr3/9a7evYVNTEyoqKgAYa1W+/PJLfPXVV1i+fHmnc48ePYr58+fj9ttvx6233orhw4ejrKwMkydPRlNTE373u98hJiYGb731Fn7961/jww8/xJVXXglJknD++edj06ZNFuOura2FSqXCli1bMGfOHADGD8zs7GyEhoYq5+r1esyaNQs5OTl49tln8d133+G5557D4MGDceedd9p8bjt27MDWrVtx3XXXISUlBQUFBXjllVcwdepUHDp0CMHBwRbnL126FFFRUVixYgUKCgqwcuVK3H333Xjvvfe6fR0B4MILL0RKSgrWrl2LP//5zwCA9957D6GhocrzM2fP+zIuLg6vvPIK7rzzTlx55ZXKv48xY8bYHMctt9yCt956C1dffTX+8Ic/4JdffsGTTz6Jw4cP45NPPrE498SJE7j66quxePFiLFy4EKtXr8aiRYswfvx4jBo1yq7nba66uhrV1dUYMmQIAOCDDz5AU1MT7rzzTsTExGD79u148cUXcfr0aXzwwQcAYNe/h7Vr16K+vh633347JEnCM888g6uuugp5eXlOnRLuVwSRHb755huhVquFWq0WkyZNEg888ID4+uuvhVar7XRuWlqaWLhwofL9m2++KQCIWbNmCYPBoByfNGmSkCRJ3HHHHcoxnU4nUlJSxEUXXaQc27BhgwAgNmzYYPE4+fn5AoB48803lWMrVqwQHd/WTU1NncY4a9YsMWjQIItjo0aNsnhc2eOPPy5CQkLEsWPHLI4/9NBDQq1Wi1OnTgkhhPj0008FAPHMM89YPJ8pU6Z0Gqc1UVFRYuzYsZ2O19XVifLycuWroaFBuW3x4sUiKSlJVFRUWNznuuuuExEREcpzX7lypQAg/vOf/yjnaLVaMWnSJBEaGirq6uqEEO2vaVxcnKipqVHOXb58uQAgsrKyRFtbm3J8/vz5QqPRiJaWFiGEEPX19SIyMlLceuutFuMpLS0VERERFseXLFnS6WdlPobw8HBx9uxZq7eZv5YXXnihCAsLE4WFhRbnmr/XrJGvZe3rzjvv7HT/tLQ0AUB89dVXFsfvueceAUD89NNPyrH6+nqRkZEh0tPThV6vF0II8be//U2o1Wrltf7HP/4h0tLSxLnnnisefPBBIYQQer1eREZGinvvvVe51sKFCwUA8ec//9nicbOzs8X48eMtjgEQK1asUL639t7ftm2bACD+/e9/K8fkf6MzZsyweN733nuvUKvVFu8Fa+R/d+Xl5eK+++4TQ4YMUW6bOHGiuOmmm5TxLVmyRLnN3vdleXl5p+fW8bFlubm5AoC45ZZbLM677777BADxww8/KMfkn+mmTZuUY2fPnhUBAQHiD3/4Q5fPWX4+ixcvFuXl5eLs2bPil19+ERdffLEAIJ577jkhhPWfwZNPPikkSbJ4z3b37yEmJkZUVVUpxz/77DMBQHzxxRfdjpOs4xQY2WXmzJnYtm0bfv3rX2Pv3r145plnMGvWLAwYMACff/65XddYvHixRVo3JycHQggsXrxYOaZWqzFhwgTk5eU5bezmNRq1tbWoqKjARRddhLy8PNTW1nZ7/w8++ABTpkxBVFQUKioqlK8ZM2ZAr9crf9WvX78efn5+Fn+Rq9VqiwLxrtTV1Vn81S+78cYbERcXp3w9+OCDAIwZp48++ghz586FEMJibLNmzUJtbS12796tjC0xMRHz589Xruvv74/f/e53aGhowI8//mjxmNdccw0iIiKU73NycgAAv/3tby2mhnJycqDValFcXAwA+Pbbb1FTU4P58+dbjEetViMnJwcbNmyw67UAgN/85jfK1Ict5eXl2LRpE26++WYMHDjQ4jZ7W6Nvu+02fPvtt/j222/x0UcfYcmSJfjnP/9pMV0ly8jIwKxZsyyOrV+/Hueeey4uuOAC5VhoaChuu+02FBQU4NChQwCAKVOmQK/XY+vWrQCMmZ4pU6ZgypQp+OmnnwAABw4cQE1NDaZMmdLpse+44w6L76dMmdLtvxPz935bWxsqKysxZMgQREZGKu+Njq+F+esmj7mwsLDLxzF3/fXX48SJE9ixY4fyX1vTX46+L+2xfv16AOj08/vDH/4AABZTkwAwcuRIi9c7Li4Ow4cPt/t30L/+9S/ExcUhPj4eOTk5yjSinH01/xk0NjaioqICkydPhhACe/bssft5XXvttYiKilK+l8fszN+V/Q2nwMhuEydOxMcffwytVou9e/fik08+wfPPP4+rr74aubm5GDlyZJf37/gBJX/ApqamdjpeXV3ttHFv2bIFK1aswLZt29DU1GRxW21trcUHvTXHjx/Hvn37bH4Yy0XghYWFSEpK6hTEDB8+3K5xhoWFKXUf5v785z8rrcMzZ85UjpeXl6OmpgavvfYaXnvttW7HNnTo0E7deSNGjFBuN+fIzwqA8vM6fvw4AGD69OlWxxMeHm71uDUZGRndniP/8u/ptCwADB06FDNmzFC+v+qqqyBJElauXImbb74Zo0eP7nJMhYWFSoBozvy1PeecczBu3DgEBwfjp59+wqxZs/DTTz/hscceQ2JiIl588UW0tLQogZB5MAUYi5s7vv+ioqK6/XfS3NyMJ598Em+++SaKi4stauusBf8df+7yB64j/x6zs7ORmZmJtWvXIjIyEomJiTbfD46+L+1RWFgIlUqlTEHJEhMTERkZ2e17HbDvtZVdfvnluPvuuyFJEsLCwjBq1CiEhIQot586dQqPPPIIPv/8807XtOcPMFvj7MnPhiwxACKHaTQaTJw4ERMnTsSwYcNw00034YMPPsCKFSu6vJ9arbb7uPkvalt/yduzTsvJkydx8cUXIzMzE3//+9+RmpoKjUaD9evX4/nnn7daU9KRwWDAzJkzO3W2yIYNG9btNeyRmZmJvXv3oq2tzWJO31atgzz23/72tzbrjLqqk+iKIz8roP3nJY/p7bffRmJiYqfzHOnW8WR31cUXX4yXXnoJmzZtsgiAejMmf39/5OTkYNOmTThx4gRKS0sxZcoUJCQkoK2tDb/88gt++uknZGZmdgp2bL3u3Vm6dCnefPNN3HPPPZg0aRIiIiIgSRKuu+46q+/97n6+9rr++uvxyiuvICwsDNdee61HlsWwNwPY2+eckpJiEUCb0+v1mDlzJqqqqvDggw8iMzMTISEhKC4uxqJFi+z6/eOscVJnDICoVyZMmAAAKCkpcdljyH/pdOymsuevwy+++AKtra34/PPPLf6CsjYVY+sX5uDBg9HQ0GDzl5wsLS0N33//PRoaGiyyQEePHu12nABw2WWX4eeff8Ynn3yCefPmdXt+XFwcwsLCoNfr7Rrbvn37YDAYLD6M5G6dtLQ0u8bYncGDBwMwdgZ2NyZnrN47aNAgAMapI2fS6XQAYDUj11FaWprVn7G113bKlCl4+umn8d133yE2NhaZmZmQJAmjRo3CTz/9hJ9++gmXXXaZk54F8OGHH2LhwoV47rnnlGMtLS12dyb21PXXX49HHnkEJSUlNot6Afvfl468V9LS0mAwGHD8+HElkwQYGwZqamqc9l63x/79+3Hs2DG89dZbFk0d3377badzfXk1676KNUBklw0bNlj9S0Oeb7d3mqcn0tLSoFarLTpoAOPChd2R/2rqmPp/8803O50bEhJi9YNh3rx52LZtG77++utOt9XU1Cgflpdeeil0Op1Fi71er7d7pew777wTCQkJuPfee3Hs2LFOt3d8/dVqNX7zm9/go48+shoAmLcvX3rppSgtLbXo5tHpdHjxxRcRGhqKiy66yK4xdmfWrFkIDw/HE088gba2ti7HJE8T9ObDOC4uDhdeeCFWr16NU6dOWdzWm7+Mv/jiCwBAVlZWt+deeuml2L59O7Zt26Yca2xsxGuvvYb09HSLqeEpU6agtbUVK1euxAUXXKB86E2ZMgVvv/02zpw5Y7X+p6fUanWn1+HFF190+SrXgwcPxsqVK/Hkk0/i3HPPtXmeve9LuVvNnveKvKil3J0p+/vf/w4AVrvRXMXa7x8hBF544YVO5zrj3wM5hhkgssvSpUvR1NSEK6+8EpmZmdBqtdi6dSvee+89pKen46abbnLZY0dEROCaa67Biy++CEmSMHjwYPzvf//rtACjNZdccgk0Gg3mzp2L22+/HQ0NDXj99dcRHx/fKWs1fvx4vPLKK/jLX/6CIUOGID4+HtOnT8f999+Pzz//HJdddpnSHtvY2Ij9+/fjww8/REFBAWJjYzF37lycf/75eOihh1BQUICRI0fi448/tnuePzo6Gp988gnmzp2LrKwsXHfddZg4cSL8/f1RVFSktMyaZ7KeeuopbNiwATk5Obj11lsxcuRIVFVVYffu3fjuu++UNXpuu+02/POf/8SiRYuwa9cupKen48MPP8SWLVuwcuVKhIWF2fvj6FJ4eDheeeUV3HjjjRg3bhyuu+46xMXF4dSpU1i3bh3OP/98vPTSS8rrDQC/+93vMGvWLKjValx33XUOP+Y//vEPXHDBBRg3bhxuu+02ZGRkoKCgAOvWrUNubm6399+9ezf+85//AADq6+vx/fff46OPPsLkyZNxySWXdHv/hx56CO+88w5mz56N3/3ud4iOjsZbb72F/Px8fPTRRxaZjUmTJsHPzw9Hjx7Fbbfdphy/8MILlcDZmQHQZZddhrfffhsREREYOXIktm3bhu+++w4xMTFOewxbOi6tYI2978ugoCCMHDkS7733HoYNG4bo6Gicc845Vmu/srKysHDhQrz22muoqanBRRddhO3bt+Ott97CFVdcgWnTpjn9udqSmZmJwYMH47777kNxcTHCw8Px0UcfWa3bcda/B3KA2/vOyCt9+eWX4uabbxaZmZkiNDRUaDQaMWTIELF06VJRVlZmca6tNvgdO3ZYnGfeOmtu4cKFIiQkxOJYeXm5+M1vfiOCg4NFVFSUuP3228WBAwfsaoP//PPPxZgxY0RgYKBIT08XTz/9tFi9erUAIPLz85XzSktLxZw5c0RYWJgAYNESX19fL5YvXy6GDBkiNBqNiI2NFZMnTxbPPvusxVIAlZWV4sYbbxTh4eEiIiJC3HjjjWLPnj12tcHLSkpKxP333y9GjhwpgoKCREBAgBg0aJBYsGCBRbuurKysTCxZskSkpqYKf39/kZiYKC6++GLx2muvdTrvpptuErGxsUKj0YjRo0d3GpPccvu3v/3N4ri8FMEHH3xgcdzWz3bDhg1i1qxZIiIiQgQGBorBgweLRYsWiZ07dyrn6HQ6sXTpUhEXFyckSVJ+brbGYH5bx3EfOHBAXHnllSIyMlIEBgaK4cOHi//7v/+z/gJ3uJb5l5+fnxg0aJC4//77RX19vcX5aWlpYs6cOVavdfLkSXH11Vcrj3/uueeK//3vf1bPnThxogAgfvnlF+XY6dOnBQCRmpra6Xxr/x6EsP5eR4dW8erqauVnHhoaKmbNmiWOHDli979RW0tQ2BpLx3/LHaFDG7wQ9r0vhRBi69atYvz48UKj0Vg8T2uvQ1tbm3jsscdERkaG8Pf3F6mpqWL58uXKcg0yWz/Tiy66yOqSGPY8n44OHTokZsyYIUJDQ0VsbKy49dZbxd69ezu9j3vy76Hjz5scIwnBCioiIiLqX1gDRERERP0OAyAiIiLqdxgAERERUb/DAIiIiIj6HQZARERE1O8wACIiIqJ+hwshWmEwGHDmzBmEhYVxeXIiIiIvIYRAfX09kpOTu92DjgGQFWfOnOm06zURERF5h6KiIqSkpHR5DgMgK+Tl14uKihAeHu7h0RAREZE96urqkJqaatf2PgyArJCnvcLDwxkAEREReRl7yldYBE1ERET9DgMgIiIi6ncYABEREVG/wxqgXtDr9Whra/P0MKiXNBpNt+2SRETkWxgA9YAQAqWlpaipqfH0UMgJVCoVMjIyoNFoPD0UIiJyEwZAPSAHP/Hx8QgODuZiiV5MXvSypKQEAwcO5M+SiKifYADkIL1erwQ/MTExnh4OOUFcXBzOnDkDnU4Hf39/Tw+HiIjcgIUPDpJrfoKDgz08EnIWeepLr9d7eCREROQuDIB6iFMlvoM/SyKi/ocBEBEREfU7DIDI5xQUFECSJOTm5np6KERE1EcxAOpnSktL8fvf/x5DhgxBYGAgEhIScP755+OVV15BU1OTcl56ejokSYIkSVCr1UhOTsbixYtRXV2tnLNx40ZIkoSoqCi0tLRYPM6OHTuU+3fF/HFCQkIwbtw4fPDBB3Y/n0WLFuGKK66w+3wiIiKAAVC/kpeXh+zsbHzzzTd44oknsGfPHmzbtg0PPPAA/ve//+G7776zOP/Pf/4zSkpKcOrUKfz3v//Fpk2b8Lvf/a7TdcPCwvDJJ59YHPvXv/6FgQMH2jUu+XH27NmDiRMn4tprr8XWrVt7/kSJiKjXdHoDWtp8tzmEAVA/ctddd8HPzw87d+7EvHnzMGLECAwaNAiXX3451q1bh7lz51qcHxYWhsTERAwYMADTpk3DwoULsXv37k7XXbhwIVavXq1839zcjHfffRcLFy60a1zy4wwbNgyrVq1CUFAQvvjiC+j1eixevBgZGRkICgrC8OHD8cILLyj3e/TRR/HWW2/hs88+U7JIGzduVG7Py8vDtGnTEBwcjKysLGzbts3BV4yIqH9q0xtw2YubccHTP+BAca2nh+MSXAfICYQQaPZAlBzkr7a7g6myslLJ/ISEhFg9p6trFRcX44svvkBOTk6n22688Ub87W9/w6lTpzBw4EB89NFHSE9Px7hx4+x7Imb8/Pzg7+8PrVYLg8GAlJQUfPDBB4iJicHWrVtx2223ISkpCfPmzcN9992Hw4cPo66uDm+++SYAIDo6GmfOnAEAPPzww3j22WcxdOhQPPzww5g/fz5OnDgBPz++7YmIuvLF3jM4UloPALj+9Z/x9uIcZKVGenZQTsZPAidobtNj5CNfu/1xD/15FoI19v0IT5w4ASEEhg8fbnE8NjZWqd9ZsmQJnn76aeW2Bx98EH/605+g1+vR0tKCnJwc/P3vf+907fj4eMyePRtr1qzBI488gtWrV+Pmm292+PlotVo899xzqK2txfTp0+Hv74/HHntMuT0jIwPbtm3D+++/j3nz5iE0NBRBQUFobW1FYmJip+vdd999mDNnDgDgsccew6hRo3DixAlkZmY6PDYiov7CYBB4ZeNJAEB4oB/qWnT47Ru/4K3F52LcwCgPj855OAXWz23fvh25ubkYNWoUWltbLW67//77kZubi3379uH7778HAMyZM8fqgoE333wz1qxZg7y8PGzbtg033HCD3WN48MEHERoaiuDgYDz99NN46qmnlMBl1apVGD9+POLi4hAaGorXXnsNp06dsuu6Y8aMUf4/KSkJAHD27Fm7x0VE1B99f+Qsjp9tQFiAH7659yKcmx6N+lYdFvxrO3YWVHl6eE7DDJATBPmrcejPszzyuPYaMmQIJEnC0aNHLY4PGjTIeK2goE73iY2NxZAhQwAAQ4cOxcqVKzFp0iRs2LABM2bMsDh39uzZuO2227B48WLMnTvXoW1C7r//fixatAihoaFISEhQpuLeffdd3HfffXjuuecwadIkhIWF4W9/+xt++eUXu65rvq2FfE2DwWD3uIiI+hshBF7eeAIA8NtJaUiMCMSamydi8Zqd2JZXiQWrt+OdW8/ziekwBkBOIEmS3VNRnhITE4OZM2fipZdewtKlS23WAXVFrTYGXM3NzZ1u8/Pzw4IFC/DMM8/gyy+/dOi65oGWuS1btmDy5Mm46667lGMnT560OEej0XALCyIiJ/klvwp7TtVA46fCTeenAwCCNX5YvWgibvn3Dmw5UYlXfzyJV3473rMDdQJOgfUjL7/8MnQ6HSZMmID33nsPhw8fxtGjR/Gf//wHR44cUQIcWX19PUpLS1FSUoLt27fj/vvvR1xcHCZPnmz1+o8//jjKy8sxa5ZzsmFDhw7Fzp078fXXX+PYsWP4v//7P+zYscPinPT0dOzbtw9Hjx5FRUWFslcbERE57mVT7c+8CSmIDwtUjgdp1Pjd9KEAgD2najwxNKdjANSPDB48GHv27MGMGTOwfPlyZGVlYcKECXjxxRdx33334fHHH7c4/5FHHkFSUhKSk5Nx2WWXISQkBN98843N6S2NRoPY2Fin7a11++2346qrrsK1116LnJwcVFZWWmSDAODWW2/F8OHDMWHCBMTFxWHLli1OeWwiIl/UpNXh492nsejN7Vi8Zgd+yatUbjtQXItNx8qhkoDbpgzudN/RKRFQqySU1rWgpLbzTIC3kYQQwtOD6Gvq6uoQERGB2tpahIeHW9zW0tKC/Px8ZGRkIDAw0MYVyJvwZ0pEvkwIge35Vfhw12ms31+CRq1l2cAFQ2Kx7JJh+NfmfKzbV4LLxybjheuyrV5rzj9+wsEzdXj5hnG4dHSSO4bvkK4+vzvq24UrRERE1GPVjVo8+NE+fHOoTDk2MDoYvxmXgvKGFry3owibT1Rg84kK5fY7Luqc/ZFlD4zEwTN12F1Y3ScDIEcwACIiIvJBW09W4N73clFW1wp/tYSrslNw9YQUTEiLUkoVbr9wMF784Tg+2l0MvUFgemY8RiTZzpxkp0bhPz+fwp6iGjc9C9dhAERERORD2vQGPP/tMbzy40kIAQyKC8E/rsvGOQMiOp2bGh2MZ67Owp1Th+C7Q2W4ctyALq89Ls24EOL+4lpodQZo/Ly3lJgBEBERkY8QQuDO/+zCd4eNi75eOyEVK349stulWjJiQ3DrhYO6vX56TDAig/1R09SGwyV1Xr0ekEdDt02bNmHu3LlITk6GJEn49NNPuzz/448/xsyZMxEXF4fw8HBMmjQJX39tuQXFo48+qmyMKX+5YusD1o77Dv4sichX1Da3KcHPquvH4emrxzh1nTpJkpBtCnr2nKp22nU9waMBUGNjI7KysrBq1Sq7zt+0aRNmzpyJ9evXY9euXZg2bRrmzp2LPXv2WJw3atQolJSUKF+bN2922pjl1YWbmpqcdk3yLK1WCwCd1kEiIvI2p6uN7emxoRrMGeOaIuVs035gu718PSCPToHNnj0bs2fPtvv8lStXWnz/xBNP4LPPPsMXX3yB7Oz2lj0/Pz+rm2M6g1qtRmRkpLKnVHBwsNPWvSH3MxgMKC8vR3BwMHeJJyKvd6bGGAAlR3be3shZ5A1R9xR5dwbIq3/jGwwG1NfXIzo62uL48ePHkZycjMDAQEyaNAlPPvkkBg4caPM6ra2tFhuB1tXVdfm4cnDFjTV9g0qlwsCBAxnIEpHXkwOgAS4MgMakRkCSgKKqZpTXtyIuLMBlj+VKXh0APfvss2hoaMC8efOUYzk5OVizZg2GDx+OkpISPPbYY5gyZQoOHDiAsLAwq9d58skn8dhjj9n9uJIkISkpCfHx8dx6wQdoNBqoVN7byUBEJDtT2wLAtRmg8EB/DI0PxbGyBuQW1WDmyASXPZYreW0AtHbtWjz22GP47LPPEB8frxw3n1IbM2YMcnJykJaWhvfffx+LFy+2eq3ly5dj2bJlyvd1dXVITU3tdgxqtZp1I0RE1GcUV7t+Cgwwrgd0rKwBu09Ve20A5JV/9r777ru45ZZb8P7772PGjBldnhsZGYlhw4bhxIkTNs8JCAhAeHi4xRcREZG3KXbDFBgAjEuLBODdnWBeFwC98847uOmmm/DOO+9gzpw53Z7f0NCAkydPIinJu5fsJiIi6o47aoCA9k6wfadrodMbXPpYruLRAKihoQG5ubnIzc0FAOTn5yM3NxenTp0CYJyaWrBggXL+2rVrsWDBAjz33HPIyclBaWkpSktLUVtbq5xz33334ccff0RBQQG2bt2KK6+8Emq1GvPnz3frcyMiInKnVp0eZ+uNDT3Jka7d2HlIXCjCAvzQpNXjWFmDSx/LVTwaAO3cuRPZ2dlKC/uyZcuQnZ2NRx55BABQUlKiBEMA8Nprr0Gn02HJkiVISkpSvn7/+98r55w+fRrz58/H8OHDMW/ePMTExODnn39GXFyce58cERGRG5WaCqAD/VWIDtG49LFUKklZBdpb2+E9WgQ9derULlfhXbNmjcX3Gzdu7Paa7777bi9HRURE5H2KzdYAcseyHuMGRmLziQrsLqzBDTlpLn88Z/O6GiAiIiLqTO4Ac3X9jyzbyxdEZABERETkA87UGKfA3BUAjTVNgeWVN6K2yfvWxGMARERE5APcsQ2GuagQDWJDjatAF1V73/6YDICIiIh8QLGbAyAASIwwBkBldS1ue0xnYQBERETkA9y1BpC5xHBju31ZXWs3Z/Y9DICIiIi8nBDCbatAm4s3BUClzAARERGRu1U2atGqM0CSgMQI1y6CaE7JANUyACIiIiI3k6e/4sMCoPFz30d7QripBqieARARERG5mbs7wGQJ8hQYM0BERETkbqerPRMAydNt8h5k3oQBEBERkZeTF0FMcXcGKMwYAFU1atGq07v1sXuLARAREZGX89QUWGSwv1JzdNbLWuEZABEREXk5TyyCCACSJLUXQntZKzwDICIiIi/niUUQZYleuhYQAyAiIiIv1tKmR2WjFoBnAqAEL10NmgEQERGRF5Onv0I0aoQH+bn98dsDIGaAiIionzlWVo9qUxaC3EuZ/ooKgiRJbn/8RC9dC4gBEBER9UpRVRN+tXITbv33Tk8PpV/yVAeYLJ5F0ERE1B8dK6uHQQDHzzZ4eij9UrGHFkGUJXIKjIiI+qNy0yrAtc1t0OkNHh6N7zIYBG5442dc+fIWNGl1yvFi0yKIniiABiyLoIUQHhlDTzAAIiKiXjHfBqGmuc2DI/FtZfUt2HKiEntO1eCF744rxz3ZAg+0b4fR3KZHXYuum7P7DgZARETUK+VmARALoV2nqKpZ+f83Nufj4JlaAJ5bBFEW6K9GRJA/AOCsF02DMQAiIqJeOVvf/qFXyQDIZU5VNSn/rzcI/PHj/WjTG1BSKwdAgZ4amrIatDcthsgAiIiIeoUZIPcoMgVA0zPjERbgh72na/H8t8fQphdQSe3FyJ6Q4IWt8AyAiIioV8xrgKqaGAC5SlG1MQAanxaFB2dnAgBe3ngSgDH48VN77iNdDoDM3wt9HQMgIiLqMSEEM0BuctpUA5QSFYTrzx2I8WlRym2eqv+ReeNiiAyAiIiox+padGjVtbe+swbIdeQMUGp0MFQqCU9cORp+KuPKzwOiPBsAJUR431pADICIiKjHyjtMeTAD5BqtOr1SYJwaFQwAGJ4Yht9dPBQALLJBnpAQ5n2rQbt/1zQiIvIZ5h1gAFDVxHWAXOFMTQuEAIL81YgN1SjHf3fxUFw1bgCSIzw8BWbKALELjIiI+gVmgNxDboFPje684WlKlHFKzJPkIujy+lboDd6xGjQDICIi6jE5ABoYbZyWqWIA5BJyC7w8/dXXxIYGQCUBBgFUNFgGxT/nVeJ/+854aGS2cQqMiIh6TA6AhieG4VRVEwMgFzEvgO6L1CoJcWEBKKtrRVldi5IRatLqcPOaHWjS6jE0PgzDE8M8PNJ2zAAREVGPyeu+DE8wfrA1t+nRrNV7ckg+ybwFvq+y1gq/8Wg5mkzvh03Hyj0yLlsYABERUY/JGaCM2BD4q411KNVcDNHp+noGCDDbFd6sLmz9/hLl/386UeH2MXWFARAREfWYHADFhwcgKtjYncRpMOfr6zVAgFkAZMoAtbTpseHIWeX27fmVaGnrO9lBBkBERNRjcht8fFggokMYALlCfUsbqk3LC6RG9+EpsA6t8JuOlaNRq0dyRCDiwwLQ0mbA7sJqTw7RAgMgIiLqEa3OoHwwx4UFKAEQp8Ccq8hU/xMZ7I+wQH8Pj8a2+A6LIX51oBQA8KtzknDBkFgAwKbjfWcajAEQERH1iNzu7K+WEBnkj6g+kgFq0xuw6M3tePx/hzw6DmeR638G9uH6H6A9A1RW14JWnR7fHi4DAFw6OhEXDDUGQJtP9J1CaAZARETUI3L9T2xoAFQqCdF9pAbocEkdNh4tx5qtBWjTG7q/Qx/nDfU/QHsXWFldK7aeqER9iw7xYQEYNzBKyQAdPFPn8feHjAEQERH1iNwCL0999JUaoIJKY8CgNwicrm726FicQX4OKX24/gcA4k0BUG1zGz7ZUwwA+NU5iVCpJMSHByIzMQxCAFv6SDcYAyAiIuoROQMU1yEA8nQNUGFFo/L/BZWNXZzpHbwlAxQe6IcgfzWA9vb32eckKbfLWaDNfaQOiAEQERH1iNwBFhdm/Mu/r9QAyRkgACio8IEAyAvWAAIASZKQEG4MhnUGgZgQDc7NiFZub68DqoAQnt8vjAEQERH1SKcMUB+pASo0y/oUmgVD3kgIoXSBpfbhVaBl8lpAAHDJqESozTZpzcmIgUatQnFNM/L7QGDKAIiIiHrEdg1Qm8fGBFhmgPrCB21vVDRo0dymhyQBA7wgAJI7wQBj95e5II0a49OiABizQJ7GAIiIiHqkqxogT01xNLbqLHYj9/YaIHn6KyEsEAF+ag+PpntyBigy2B/nDYrpdPuUYcZpsJ/6QB0QAyAiIuqRjgFQZLBxkT69QaCuReeRMclTXvK+ZKerm726FV4ugO7rawDJRiWHAwAuz0qGv7pziDFlSBwA4OeTldB5+Ofi59FHJyIirySEaN8HzBQABfqrEaJRo1GrR3WjFhFB7l+1WK7/GZkcgaOldWhpM+B0dTMyYkPcPhZn8JYWeNncMckYEBmEMSmRVm8flRyOqGB/VDe1Ye/pGoxPi7Z6njswA0RERA6rbW6D1vQXfGxogHI8OtQ4DVbpoUJouf4nIyYY6TEhpmPeOw3mLS3wMpVKwoT0aGj8rIcXKpWEyUP6xjQYAyAiInKYnP2JCPJHoH97bYrcCVbtoQBIzgClxYS0B0BeXAjtLS3wjphiCoC2nqj06Dg4BUZERA4726H+R6asBeShxRDlbE96bDBadHoA3t0K700t8Pa6eEQC/rUwADlWiqTdiQEQERE5rGP9j8zTGaBTlXLRcAha24xTdN7aCq/TG1BcYwqAfCgDFBcWgItHJHh6GJ6dAtu0aRPmzp2L5ORkSJKETz/9tMvzP/74Y8ycORNxcXEIDw/HpEmT8PXXX3c6b9WqVUhPT0dgYCBycnKwfft2Fz0DIqL+qX0V6A4BkAdXg25p0+NMrXFc6THBSDNNgRV6aQ1QSW0L9AYBf7VkscAgOYdHA6DGxkZkZWVh1apVdp2/adMmzJw5E+vXr8euXbswbdo0zJ07F3v27FHOee+997Bs2TKsWLECu3fvRlZWFmbNmoWzZ8+66mkQEfU7tjJAntwOQy4YDgvwQ3SIRun8KvLSVni5/iclKthiRWVyDo9Ogc2ePRuzZ8+2+/yVK1dafP/EE0/gs88+wxdffIHs7GwAwN///nfceuutuOmmmwAAr776KtatW4fVq1fjoYcectrYiYj6M1s1QJ7cEFXuAEuLDVb2pQr0V6GlzYDi6make1kr/GlT/U+KD9X/9CVe3QVmMBhQX1+P6GjjOgJarRa7du3CjBkzlHNUKhVmzJiBbdu22bxOa2sr6urqLL6IiMi29gyQ5dRMlAf3AzPvAAOMm3PKnWD5XjgNdtoHO8D6Eq8OgJ599lk0NDRg3rx5AICKigro9XokJFgWVyUkJKC0tNTmdZ588klEREQoX6mpqS4dNxGRt+u4CrQsJtR2AFRU1YT6FtftEyZ3e6XHtAcM3twKX1ZnfI0TWf/jEl4bAK1duxaPPfYY3n//fcTHx/fqWsuXL0dtba3yVVRU5KRREhH5Jptt8DYyQEVVTZj+3EZc/NyPOHTGNVl2uQU+Lbp9qist1hgMeWMrfGWj8TU2X2iSnMcrA6B3330Xt9xyC95//32L6a7Y2Fio1WqUlZVZnF9WVobExMSOl1EEBAQgPDzc4ouIiKxr1elR22zM5HRqgzfVANW16CwKjzcePYs2vcDZ+lZc+89t2HrS+asAy0FOmlkGKEOeAvPCDFBFgzGIlLNq5FxeFwC98847uOmmm/DOO+9gzpw5FrdpNBqMHz8e33//vXLMYDDg+++/x6RJk9w9VCIinyRPf2nUqk77fUUE+UMyNSzVNLVPd/2cVwUACNGoUd+qw6LVO/DF3jNOG5NWZ1BqZsyLnb25Fb49A8QAyBU8GgA1NDQgNzcXubm5AID8/Hzk5ubi1KlTAIxTUwsWLFDOX7t2LRYsWIDnnnsOOTk5KC0tRWlpKWpra5Vzli1bhtdffx1vvfUWDh8+jDvvvBONjY1KVxgREfWOef2PJFm2Z6tVUqdpMCEEfs4zbnvwzxsn4NLRidDqDVj6zh78a3O+U8ZUXNMMgwAC/VUWWSlvboWvqDdlgEI4BeYKHg2Adu7ciezsbKWFfdmyZcjOzsYjjzwCACgpKVGCIQB47bXXoNPpsGTJEiQlJSlfv//975Vzrr32Wjz77LN45JFHMHbsWOTm5uKrr77qVBhNREQ9I9f/xIZZ/2COCjZmheQA6MTZBlQ2ahHor8LEjCi8OH8cFk5KAwA8/r9DTpkOU7bAiAmxCMriw4yt8HqDQLFpZ3Vv0KTVobnNuJWHrdeZesej6wBNnToVQgibt69Zs8bi+40bN9p13bvvvht33313L0ZGRES22FoEURYdosHJ8kZlLSA5+zM+LQoBfsaNUx/99SicqmrChqPlOFxSj8mDY3s1pvYtMCxbxlUqYyv8kdJ65Fc2es1aQJWm+p8APxVCNOpuzqae8LoaICIi8ixbHWCyjlNgcv3PeRntm19KkoQBpgX+5ILq3mjfBLVzgCMXRRd6USF0eUN7B1jHaUZyDgZARETkkO4yQOZrAZnX/5w32HL3b7mAus4JAZC1DjCZHBQVeFErfCU7wFyOARARETmk3MZGqDLzDNBxs/qfrJRIi/PkAMipGaCYzhkgb2yFr2zgGkCuxgCIiIgccsq06WhypPU9qsz3A5OzPxPSoqHxs/zIcVYApDcIZSNUaxkgb2yFr2yUO8CYAXIVBkBERGQ3g0EoU0mDY0OtnmOeAVKmvwZFdzovPNA5AVBJbTPa9AIatQpJEZ2DMm9shZenGWOYAXIZBkBERGS3M7XN0OoM8FdLSI60vkdVtKlupbJB214APSim03nOygDJ9T8p0UFQqzoXDHtjK7ycAeIiiK7DAIiIiOwm19EMjA6Gn9r6R0i0KQN0tKweVY1aBPmrMaZD/Q8AhDupCLqr+h+gvRUe8J5d4VkD5HoMgIiIyG5yAJRhY/oLaK8B0huM67xNSI/qVP8DOD8DZK3+RybfllfuLQEQu8BcjQEQERHZTQ4gBsXZXlAwqkPhrrXpLwCIMK0Y3aozoMW06nFP/JJvnGYbHGc7KBuZFAEAOFBca/OcvqTClAHiNhiuwwCIiIjsJk83ZXSxonKIRm2R8bFWAA0AoRo/yCU7PZ0GO3SmDnuLauCvlvCrcxJtnpeVagyAcotqevQ47qQ3CFSZVtGODWMGyFUYABERkd3kKTBb9TaAcZVnuQ7IVv0PYKzNCe/lNNi7O4z7Rc4cmdBlvczYVOMY8isaUWMKLvqq6iYt5F2i5NeRnI8BEBER2UWrMyjr7XQ1BQa0T4NNSI+Cv41iaaB3rfDNWj0+2VMMALhu4sAuz40M1iDdVAe093TfngaT63+igv1tFppT7/GVJSIiu5yqaoJBAMEatc1tMGTyAn626n9kvSmEXr+/BPUtOqREBeGCId1vpipngfb28Wkwpf6HHWAuxQCIiIjsUlDRXv/T3QadiyanY+rwOFwzIaXL85T9wFocD4De2W6c/rpuYipUVtb/6SjLFAD19TqgCqUFntNfruTn6QEQEZF3yK/ovgBaNmNkAmaMTOj2PCUD1ORYAHS8rB47C6uhVkm4ZkKqXfcZaxYACSH67C7r7S3wzAC5EjNARERklzwHAiB7tRdB6xy637s7igAA0zPjkRBufUXqjkYmh8NfLaGqUYvTfXhF6MpGUwaI+4C5FAMgIiKyS35FAwDnBkA9qQFqadPjo92nAQDzz7Uv+wMAAX5qjEwKBwDs6cPTYBX1zAC5AwMgIiKyiyNTYPbqSQD09cFS1DS1ISkiEBcNi3fo8byhEFrJADEAcikGQERE1K3GVh3K6owfzJ4OgN7dbpz+mjch1ermp13xhkLoCm6D4RYMgIiIqFvyCtBRwf6IdOLifOFBxl4ce1eCLqltxra8SkgSMG+i/dNfMjkDdKC4Fm16g8P3dwd2gbkHAyAiIuqWK6a/AMfb4I+U1AMAhieEYUBkkMOPlx4TgvBAP7TqDDhaWu/w/d1B7gLjFJhrMQAiIqJu5Zd3vwt8Tzg6BXay3FiI3d1K1LaoVFKfngZr0urQbNoYlkXQrsUAiIiIuiVngHoaeNjieABkHEdXO793J7sPB0By9ifAT4UQjdrDo/FtDICIiKhb+XbsAt8TcgDUpNXbVZOT18sMENC3C6HLG9o7wPrqQo2+ggEQEZEPatbqsXjNDvzn50KnXM+eXeB7Isy0GSpgXxZIXoxxUC+m4uQA6GR5Q4+24HCl9vofFkC7GgMgIiIftPlEBb4/chaP/+8Qapq0vbpWdaMWNaatKtJjg50xPIVaJSEs0NgJ1l0AVNfShvJ6Y4akNxmg2NAApEQFQQhgfx/bGb6SG6G6DQMgIiIfJAcKrToD3jNtG9FTctYlKSIQwRrnbyFpbx1Qnqn+Jz4swCJz1BNj++g0WGWjaQ0gboPhcgyAiIh8kBwAAcDbPxdCbxA9vlaBi1rgZeGB9gZAva//kfXVAEj+uTED5HoMgIiIfFB5Q4vy/6erm7HhyNkeX8tVawDJlLWAugmA5Bb43nSAyeQAaM8p487wfYWcAWINkOsxACIi8kFyJkH+IH1rW0GPr9VXAiB5CmyQEwKgcwZEIMhfjYqGVhwu6TsLIlY2cB8wd2EARETkg+QA6I6LBkOSgJ+OVygZFEfluSkAsrcGyBlTYIH+apw/JAYA8MORsl5fz1kquQ+Y2zAAIiLyQfJ6MtkDI3FxpnHH9Le3Od4SL4RweQ1QRHD3AZDeIJS1iIY4IQMEANMzEwAAP/RietDZ5H3AYkKYAXI1BkBERD5GCIGKemMmIS40EAsmpQMAPtx1Gg2tOoeuVVbXiuY2PdQqCanRzm2Bl9mTASquboZWZ4DGT4XkHuwBZs20zDgAwJ6iGlQ19m6pAGfQGwSqTEsWxIYxA+RqDICIiHxMo1av7CcVG6bBBUNiMSguBA2tOnyy+7RD18qrME6bpUYFwV/tmo+McDsCIHn6LiMmBGqVc1ZITooIwoikcAgB/HjM81mg6iYt5Hrs6GAGQK7GAIiIyMfI9T+hAX4I1vhBpZKw4Lw0AMBb2wod6noqrGwCAKS7aPoLAMLtWAhR6QCLd+44ppuyQN8f9nwAJNf/RAX7w89FwSa14ytMRORj5AAoLqy9juQ341MQolHjxNkG7CiotvtaBS7aAsNcexeY7ek5Z2yBYY1cB7TpWDl0duxF5koVXAXarRgAERH5GCUAMvsgDQv0x/QRxg/7X/Iq7b5WQaUcALmm/gewrwbo5FnnLYJobmxqJKJDNKhr0WFXof2BYW/pDQLHy+otsnEVDZZLF5BrMQAiIvIx5fXGRRA7FtLKi//tPV1j97XkKbA0F06B2bMOkJwBcsYiiObUKgkXDTNOg7mzG+yfm05i5vOb8Nw3x5Rj7S3wzAC5AwMgIiIfI7fAx3X4IB2bGgEAyC2qtasOSAhhlgFyfQBU36qzumWHszZBtWWaaZkAdwVAQgh8uNNYjP7KjyeVDVkrG00ZIO4D5hYMgIiIfIzSAh9mGQCNSo6An0pCRUMrztS2WLurhbK6VrS0GaBWSUiJck7ruTVyFxhgPQskL4AY54RNUK25aGgc1CoJx882oKiqyenX7+hYWYOS0dIbBO7/cC+0OoPyc+Mq0O7BAIiIyMcoGaAOAVCgvxqZSWEAgNxTNd1eR87+pLiwBR4A/NUqhGjUAKzXAeUpe4C5biHG8WlRANyTBfryQAkAYGJ6FKJDNDhSWo9XfzypZIA4BeYeDICIiHyMtS4wWVZKJAD76oAKTQFQmgunv2RdFUI7cw8wW6a7cRrsqwOlAIB5E1KxYu5IAMCLPxzH/mLjVBi3wXAPBkBERD6mvQsssNNtWaZC6Nyimm6vk19hnA7KcGEHmEyeBqtr6RwAyWsADXJhIbYcAG3Lq0ST1rHVsh1RUNGII6X1UKskzBiRgF9nJWPGiHi06QXK6tgF5k4MgIiIfIjBIJR2amsZoGxTALT/dG236964MwPU1WrQcgbI2R1g5obGh2JAZBC0OgO2nLB/mQBHfWnK/kwaFIOoEA0kScJfrhiNsAA/5RzWALkHAyAiIh9S09wGnamTytpUyqC4UIQG+KG5TY/jZ7veHb5AWQXa9RkgW1Ng5pugujIAkiQJU4bGAgD22pEd66mvDhoDoF+dk6gcS4wIxMNzRijfswbIPfy6P4WIiLyFPP0VFexvtXBZrZIwekAEtuVVYm9RDUYkhVu9jhBCyQC5sgVeZisAMt8EdYALO9GA9t3uT7moE+xMTTP2FtVAkoBLRiVY3HbtxFScqW1BiEaN0AB+NLsDM0BERD6kqwJomT11QOX1rWjS6qGSgJQoz2WATpo2Y02PCXbaJqi2pJlqnVwVAMnFzxPSohAfZlmfJUkSls0chtsvGuySx6bOGAAREfmQrup/ZGPtCIDk6a8BUUHQ+Ln+o8LWatDyFhiunP6SpUa7OABSpr+SXHJ9cgwDICIiH2JtH7CO5ADoWFm9zY4nd2yCas5WBkjZBNVFawCZG2gKgKoatai30o3WG+X1rdhRUAXAsv6HPIcBEBGRD7G1CKK5xIhAJIQHwCCgbMPQkTu2wDBna0d4uSB5eKL1WiVnCgv0R7RpGwpnZ4G+OVQKIYCslAgMiHRtLRPZhwEQEZEPkTNA3bVSd7cxqrIJqhvWAAKA8CBj4a95Bqi2uQ2HSuoAADkZ0W4ZhzwN5uwtMeT6n1nM/vQZHg2ANm3ahLlz5yI5ORmSJOHTTz/t8vySkhJcf/31GDZsGFQqFe65555O56xZswaSJFl8BQZ2XgyMiMgX2VMEDbQXQu8tsp4ByjdNPWW4cPFBc9amwHYWVEEI4xgSwt3zezzNBXVA9S1t2HbSuLbQbNb/9BkeDYAaGxuRlZWFVatW2XV+a2sr4uLi8Kc//QlZWVk2zwsPD0dJSYnyVVhY6KwhExH1afYGQGNNW2JYK4Q2b4F3xyKIgPUA6Oc8Y9Bw3iD3ZH+A9jogOQPmDPtP10JnEBgQGeS2gJK659HFBmbPno3Zs2fbfX56ejpeeOEFAMDq1attnidJEhITmWYkov7HnhogABidEgFJAoprmlFe32pxfkWDFo2mFvjUaPfUq5hvhWEwCKhUEn7JNxYN52TEuGUMADDQBa3we011VmNSIpx2Teo9n6wBamhoQFpaGlJTU3H55Zfj4MGDXZ7f2tqKuro6iy8iIm/TpjegukkLoOsuMMBY8DvE1FreceVjOfuTHBmEAD+18wdqhZwBEgKob9WhrqUNB0ybg+Z4IAPkzABon6nOSp52pL7B5wKg4cOHY/Xq1fjss8/wn//8BwaDAZMnT8bp06dt3ufJJ59ERESE8pWamurGERMROUdVoxZCGFd7jgrufkPNLBuF0PluboEHgAA/NQL9jR9Jdc1t2FVQDYMwFmEnRbiva0ou+i6ubu52rzR77WMGqE/yuQBo0qRJWLBgAcaOHYuLLroIH3/8MeLi4vDPf/7T5n2WL1+O2tpa5auoqMiNIyYico72DjANVHasmix3gu0qrLY47u4OMJl5HdDP+cb6H3d1f8kSwgKhUaugMwiU1Lb0+noVDa0ormmGJAGjBzAA6kt8LgDqyN/fH9nZ2Thx4oTNcwICAhAeHm7xRUTkbextgZdNGmysrdl6slKZpgHa1wByd8FueGD7atA/5xnrf84b5L76HwBQqSSkmOqenDENJr+ug2JDEGZ6ftQ3+HwApNfrsX//fiQlsfWQiHybvR1gssFxobgqewAA4K/rDkMI4y7yBW7uAJPJGaAztS1m9T/uDYAA57bCy8sMZJm67qjv8GgA1NDQgNzcXOTm5gIA8vPzkZubi1OnTgEwTk0tWLDA4j7y+Q0NDSgvL0dubi4OHTqk3P7nP/8Z33zzDfLy8rB792789re/RWFhIW655Ra3PS8iIk9QOsDszAABwB9mDUeAnwq/5Ffhu8NnjS3wFcYP/nQPTYH9cKQMeoNAanSQR1ZNdmYrvJwBYv1P3+PRNvidO3di2rRpyvfLli0DACxcuBBr1qxBSUmJEgzJsrOzlf/ftWsX1q5di7S0NBQUFAAAqqurceutt6K0tBRRUVEYP348tm7dipEjR7r+CREReZCjGSAAGBAZhMUXZODljSfx5JeHMSYlAvWtOkhS+6rI7iIHQBuPlgNwb/u7uYGmzFdvV4MWQigt8OwA63s8GgBNnTpVSblas2bNmk7HujofAJ5//nk8//zzvR0aEZHX6UkABAB3Th2M93YUIa+8EU99eQQAkBwRhEB/97TAy+S1gJq0egDur/+ROasV/nR1M6oatfBTSRiRxNrSvsbna4CIiPoLexdB7Cgs0B+/nzEUAPDJnmIA7u8AA9ozQDJ3d4DJ5Ocur4fUU3L7e2ZSmNuDSeoeAyAiIh9R4WAXmLn55w7EILOur3QPbNlgHgANiAxy+xScLDXK+Lh1LTrUNrV1c7Zt7fU/kU4YFTkbAyAiIh/R0ykwAPBXq/DQ7Ezle3cXQAOWAZA7V3/uKEijVl7DwqqeZ4HkBSazWADdJzEAIiLyAc1aPepbdQB6FgABwMyRCbhgSCwAYHxalNPGZq9wswDIU/U/st62whsMAgeKjdsqMQPUN3m0CJqIiJyjwlT/E+CnQlhAz361S5KENxZOQGFlE4YnhjlzeHYxzwCd56EOMNnA6GDsLKzucSt8XkUDGlp1CPRXYWh8qJNHR87AAIiIyAecNZv+kqTut8GwJdBf7ZHgBzDuPC9JxhWo3bULvS3yrvA9bYWXF0AcPSACfmpOtvRFDICIiLxMVaMWN7zxCyYPjsHDl46ASiX1qv6nr0iKCMLnSy5ATKimV0GcM/S2FZ4F0H0fAyAiIi/zS14lDpfU4XBJHWqa2vDM1WOUKTBHVoHui0b3kYLh9lb4ngVAudwBvs9jAERE5GWqzVqzP9p9GnqDASmm1u1YL84A9SVyC35JbTO0OgM0fvZPY2l1Bhw+YyyA5h5gfRcDICIiL1PdpAUADI4LQWFlEz7NPYNAf+MHtLdngPqKuNAABPmr0dymR3FNMzIcWBfpaGk9tHoDIoL8PbKgJNmHlVlERF6mutEYAM0YkYBVN4yDv1pCS5sBgHfXAPUlkiT1uA5or9kGqJ6uZSLbHM4A5efn46effkJhYSGampoQFxeH7OxsTJo0CYGBga4YIxERmZGnwCKDNZg1KhGv3DAed/13N7R6g0d2T/dVqdHBOFpW73AAdPBMewcY9V12B0D//e9/8cILL2Dnzp1ISEhAcnIygoKCUFVVhZMnTyIwMBA33HADHnzwQaSlpblyzERE/VqNaQosKti4bs6MkQl457Yc/JJfhSlDYz05NJ8iT1+dcnBPsKOl9QDADVD7OLsCoOzsbGg0GixatAgfffQRUlNTLW5vbW3Ftm3b8O6772LChAl4+eWXcc0117hkwERE/V2VHACFaJRj49OiMT7Nc9tH+KKeTIEJIXC8rAEAMCzBM+spkX3sCoCeeuopzJo1y+btAQEBmDp1KqZOnYq//vWvKCgocNb4iIiogxrTFFhUsKabM6k3BvagFb6ktgX1rTr4qSSHCqfJ/ewKgLoKfjqKiYlBTIxnlzAnIvJl1R2mwMg1zDNAQgi7CpqPlRmnvzJiQxxqnSf3c/ino1arcfbs2U7HKysroVarnTIoIiJXeWtrAf7w/l7o9AZPD6VH9AaB2ub2ImhynZQo49YcTVo9Khq0dt1HDoA4/dX3ORwACSGsHm9tbYVGw3+MRNR36Q0CT315BB/tPo0dBdWeHk6P1Da3Qf41HMkMkEsF+KmRHGHsqiu0sxD6GOt/vIbdXWD/+Mc/AJh2C37jDYSGtu9uq9frsWnTJmRmZjp/hERETnKyvAHNbXoAQH5FIyYNdu10fWVDK8obWpGZ6LxuIHn6KyzQD/7cZNPl0mODUVzTjILKJkxI777IvD0DxB3g+zq7A6Dnn38egDED9Oqrr1pMd2k0GqSnp+PVV191/giJiJxkv2l/JgDIr2hw+eMtWbsb2/Or8OmS8522KWZ7Czwz7u6QFhOCLScq7WqFNxjMOsASmQHq6+wOgPLz8wEA06ZNw8cff4yoqCiXDYqIyBX2F5sHQI6t7eKoljY9dhRUwyCAT/YUOy0Aqm6UO8A4/eUO6aZOsAI7OsFOVzejuU0PjVqFtGhugdHXOZw/3bBhA4MfIvJKB8wCoLxy1wZAh0vqoDcYi3W+OlAKg8F6/aSjrK0BRK6TFmNsZbenBkie/hocHwo/Tk/2eXb9hJ566ik0Ndm3DsIvv/yCdevW9WpQROSdDp6pxdwXN2Pj0c6dop6mNwgcNO3QDRhbm9tc2AlmHmyV1LYg17Q/VEf/+P44rli1RdnfqzucAnOvdFMAZE8G6Cjrf7yKXQHQoUOHkJaWhrvuugtffvklysvLldt0Oh327duHl19+GZMnT8a1116LsDDOfRL1Rx/vLsb+4lq8v7PI00PpRC6ADtGoEeivgs4gcLq62WWPd6DYGGzJS8d8ub+k0zkltc34x/fHkVtUg28Oldp13fZ9wDgF5g7yWkC1zW1K8GnLcbbAexW7AqB///vf+O6779DW1obrr78eiYmJ0Gg0CAsLQ0BAALKzs7F69WosWLAAR44cwYUXXujqcRNRH3TIlGEpqnJdYNFTcgH0qOQI5a/6vHLXFULL9UZXjB0AAFi/v7TTMiJvbyuEzjQ1Zm9bPjNA7hWkUSMx3LjRd3dZoKNsgfcqdhdBx8XF4fXXX8c///lP7Nu3D4WFhWhubkZsbCzGjh2L2FhuwEfUnwkhcKjEGAA5unu2O8gByTkDIlBW14IjpfUuK4RuadMr9SB3Tx+Crw6UorimGfuLa5Vi6GatHu9sP6XcZ2dBlV3XrmpkDZC7pcUEo7SuBYWVjRibGmn1HL1B4KQpoB7OAMgr2F2lNWrUKKxduxYqlQpjx47F5Zdfjuuuuw4zZsxg8ENEOFPboqxQXNvcpvx/XyHX5IxOCVf2aMpzUQB0tLQeOoNAdIgGg2JDMD0zHoAxCyT7NLcY1U1tSAwPhCQZswtn61u6vXZ1E7vA3E2pA6qwHdgXVjZCqzMgyF+NlKggdw2NesHuAOivf/0rbr/9dlxzzTWoqrLvLxUi6j8OmRUYA0BRH8oCmRdAjx4QoQRA+S7qBDPPNkmShNmjEwEAXx4ogRACQgi8ucW4tMgtUzKUhRJ32TENxikw90uLlTdFtf1+kTN+Q+JDoVJ1v2cYeZ7dAdBdd92Fffv2obKyEiNHjsQXX3zhynERkZfpGACdru47AVCeqQA6WKNGRmwoBsXJGSDX1AAp2aYBxsBm2vB4BPipUFjZhEMlddhyohLHyhoQolFj3sRUTEw3Li1iTx1QVSOLoN2tvROsqwCI9T/exu4aIADIyMjADz/8gJdeeglXXXUVRowYAT8/y0vs3r3bqQMkIu9wqMT4oS9JgBB9qw5IzsiMSg6HWiUpGaCyulY0tuoQEuDQr0K7H2/0gAgAQEiAH6YOj8PXB8vw5f5SpVbq6vEpCA/0x4T0aPx7WyF2dFMHJIRQMkDRrAFyG7kTrLCLImi2wHsfh//VFxYWKitBX3755Z0CICLqn+QP9XEDo7CrsLpPdYKZT0kBxl3Uo0M0qGrUIr+iUTnuDK269gJo8+teOjoJXx8sw7s7ilDR0AoAWHR+BgAoGaCDZ2rR0KpDqI2ArKFVp3SNcQrMfdJMq0FXNmpR19KG8MDO2TelBZ5bYHgNh6KX119/HX/4wx8wY8YMHDx4EHFxca4aFxF5kdrmNiXguWRkAnYVVvetDNBpy4wMAAyKDUFVoxZ5Tg6AjpbWo00vEBXsjwGR7cWw0zPjofFTKcHPxZnxSiYqKSIIAyKDUFzTjNxTNbhgqPXGkhpTAXSgvwqB/mqr55DzhQX6IzZUg4oGLU5VNnV6v2h1BmVlcU6BeQ+7a4B+9atf4cEHH8RLL72Ejz/+mMEPESkOm7I/AyKDlCCjqI/UAHUsgJa5qhC6YwG0LCzQHxcObf+9eZMp+yM7N8O403hX02BKCzyzP26X1kUdUEFlI3QGgdAAPyRHBLp7aNRDdgdAer0e+/btw4IFC1w5HiLyQnIB9MjkcKSa6iVOVzU7bf+r3jAvgB4U116fkWEqhLZnV/i88gb8/t09+Ol4ebfnHugw3WZublYSACAzMQznD4mxuG2CaRpsZ6HtAKiaHWAeI0+DWasDOlpqnP4amhBqEfRS32b3FNi3337rynEQkReT639GJoUjKSIQapUErd6As/WtSPTwX8RyRmZkkrEAWjZIzgB1sxbQlhMVuPM/u1DXokNpbQumDO06+92xANrcr7OSYRACE9KiO31QTkw3ZoD2nKpBm94AfyubacpTYFEh7ABzt/QuNkWV63+4AKJ34Xa1RNRr5hkgP7UKyZHGoKcv1AF1LICWydmgvPLGTltUyP77SyEWrN6OuhYdgO6fj1ZnULIB1gIgSZJwZXaKkiUzNyQuFBFB/mjS6jstKSCTM0CRzAC5nZwBsrYdhtwBNpQBkFdhAEREvaLVGXD8rPEDYGSScd0buW24LyyGeMBGRmZgdDAkCahv1aGiwXKTS53egMe+OIiHPzkAvUFg5sgEAEBpXQta2vQ2H+tYmbEAOiLI3+HVgFUqCRPS5PWArE+DVSs1QMwAuVvXGSB5DSC2wHsTBkBE1CsnzjagTS8QFuinfOinRhkDIFdngAoqGtGmN9i83aIAOsUyAAr0VytdWh2nwe77YC/e3FIAALh/1nC8duN4hGjUEAIorrHd3m8+/dWTWpAJpmmwnTYWRJS3wYhmBsjt5ACorK4VTVqdcrxJq1MKozkF5l0YABFRr5jX/8gf+vIUjys7wX44Uoapz27E018esXlOfkUDmrR6BPmrMTiu81/n7dNg7YXQ+0/X4tPcM1BJwCs3jMOSaUMgSVL7c+oiqLM13WaviWaF0Nam5TgF5jkRwf7K6tvmgf1Hu07DIIwZxbiwAE8Nj3qAARAR9Yp5/Y/MnmChtw6XGKfdtp6stHmOUgCdbFkALbNWCP3ShuMAgMvHDsDs0UnKcXuek63pNnuNTokwrRWktVprwiJoz0rrsCmq3iDwxmbjnm6LL8hgB5iXYQBERL0ib4Eh1/8A5jVArlsNWq6HMU7BWZ8Gk4Ozc8yCM3Mdd4U/UlqHrw+WQZKAJdMGW5wrPydb03panQFHSmwXQNsjwE+NsSmRAKzXAcnrADED5BnpMZabon59sBSFlU2IDPbHNRNSPDk06gEGQOQWOr0BL3x3HBuOnPX0UMiJhBDWM0CmWqDuioZ7Q66H0eoNKLDRyq5Mz9kIgAbFWWaAVm04CQC49JwkDIm3rOfoLgA6VlYPrd6A8EA/pEY7VgBtTl4PaEd+5wBI2QeMAZBHtC+G2AQhBP65KQ8AsOC8NARruC2Ut2EARBbyyhvwwc4im23BPfXxnmI8/90xPPrFQadelzyruKYZdS06+KslDDULGKJDNAjRqJVzXEEOBgDgsKn13Jx5cDYiqesMUGFlI46X1eN/+84AAJZMG9LpXDmoOWUjqyUHWx1XgHbUeFMnWG5RTafb5KCPCyF6Rlp0ewZoR0E19hbVIMBPhQWT0z07MOoRBkBk4Y+f7Mf9H+7DVwdKnXZNIQTe+Mn4l1JFfavTrkueJwcYQ+LDoPFr/3ViXjTsqk6wKrMA6EhJ53VzyupaUd3UBrVKsrk/U3JEEAL8VGjTCzz86QEIAcwYkWA1YzRQWeG6yeofCPL6P5mJ1oMte8mPnVfRaJE9a2nTo9n0fSRrgDwiPbZ9NejXNhmzhb8Zn4LYUBY/eyMGQKQQQuBgsfGDZFue7cJSR/14rBzHTOtkNGr10Opsty2TdzHvAOso1SxgcAW5IBgAjljJAMm1SYPjQmxuHKpSSUoWaLtpyunu6Z2zPwCQYmrtr2/VWTy2TN4Bfnhi79aCSQwPRGSwP/QGgRNn27vT5A4wP5WEMBu7xZNryVNgxTXN+O7wWUgScMsFGd3ci/oqBkCkOFvfivpW4/oWttYh6Yk3fsq3+N586oK8m7X6H1l3NTO9VW32PjpqJQCSu8RsTX/J5AAIAKYMjcXY1Eir5wX6q5EQbvxL39pzksfQ293AJUnCCFMW6bBZZqu60Rh0RQZr2G3kITEhGoSaBZ8zRyRY7C9H3oUBECnk1UwBYzdMXUvnv3IddfBMLTafqIBaJSHANEVS09z761Lf0GUGyFQI7YpOML1BoNbsfVRc02zxPWAWnDkQAC2dPrTLc20FddWNWpw1Te86YzuEzCTjNeQgDmj/w4GrQHuOJEnKlhgAcPtFgzw4GuotBkCkkLczAACDMG7K2Fv/MmV/Zp+TqKy6K7cvk3czGIRS4Dw4PqTT7QNjXJcBqmtug1yGIy8+1zEL1F0HmGyMqe180qAYnJsR3eW5tla4lveCSo0OssgQ9JSctTpS2p4BquJO8H2CvCL0+LQojE/r+v1CfRsDIFKY1xsAwE4b+xHZq6S2GZ/vNXbV3HbhIGUVVWaAfEODVqcEIeGBnbMScrBQZKNouDfkYCAswE9Zc8c8WGhsbd+eoLspsFmjErB60QS8euP4bh9XqWvqsML1MSfvBi5nrQ6X1CmvndwBFskMkEddPSEFwxPC8PCcEZ4eCvUSAyBSHDcFQPJfwbY2ZLTXmq0F0BkEcjKiMSYlUlm8jTVAvqHBtEO6Rq2yWmRsXjTccXqqt+T3UGSIPzITjUGHeSH0kdJ6CAHEhwV026EjSRKmZyYgIqj7wMLWFJiz6n9kQ+JDoVZJqG5qQ1mdcWqtxpQ5jQ5hBsiTpg2Px9f3XohxA6M8PRTqJQZApJAzQNefOxCAcR2Srjaa7EpDqw5rfzkFALh1inGeXP7LtdpKBw15n3pTABQWaH3KJ0ijVqannD0NJhcERwdrMFwOgMwKhu2d/nKUrWk9OQCSx9Jbgf5qZZuOw6bMVnsGiAEQkTN4NADatGkT5s6di+TkZEiShE8//bTL80tKSnD99ddj2LBhUKlUuOeee6ye98EHHyAzMxOBgYEYPXo01q9f7/zB+5jKhlZUNWohScAloxIQGeyPljaDspO2o97bUYT6Fh0GxYVgemY8ACAySM4AMQDyBQ2txp9jqI0ACHDdlhjmm4LKU1xHS+thMBini+Tuqe6mvxwlP58zNS3KHwdCCKUGyFkBEABkJll2glWzCJrIqTwaADU2NiIrKwurVq2y6/zW1lbExcXhT3/6E7Kysqyes3XrVsyfPx+LFy/Gnj17cMUVV+CKK67AgQMHnDl0nyNnfwZEBiFY44cJptVoe1oH9FluMQDg5vMzoDJtQin/4uYUmG+o6yYDBLR3gjk7A6RsChrsj4zYEGjUKjRq9ThdbQy07O0Ac1RcaAA0firoDQIlNS0AjNt91Lfo4KeSMCjWeS3RI5LkzJYxuKpmETSRU3k0AJo9ezb+8pe/4Morr7Tr/PT0dLzwwgtYsGABIiKsbzb4wgsv4Fe/+hXuv/9+jBgxAo8//jjGjRuHl156yZlD9zly/c/QeOMv8AnpPa8Dqm1qU3bhnjkyQTkeaapdqGYA5BPkGqCuup6UDFC1cwOgKrMMkL9ahSGm9+2R0jroDUIpiHZ2BkilkjoFdfL0V0ZsiMVq2L01olMGSN4JngEQkTP4XA3Qtm3bMGPGDItjs2bNwrZt22zep7W1FXV1dRZf/Y2cAZLXMJmYLmeAqh3u4Pk5vxJCGFfgTQgPVI5HBskZIE6B+YL2GiDbUzIp0e2dYM5U0yEbIq+bc6S0HgWVjWhpMyDQX2Wxxo+zdAzqnF3/I5MXQ5S3xOA6QETO5XMBUGlpKRISEiyOJSQkoLTU9t5WTz75JCIiIpSv1NRUVw+zz5EDoCGmVU3PGRABjZ8KlY1aZadse207adxGY/LgWIvjUcGsAfIlcg1QV9syDHRRACQXQUeZ9sRq7wSrU6a/MhPDoVY5f8Xkjp1gR53cAi9LCA9AlNmWGFWN7VkvIuo9nwuAemL58uWora1VvoqKijw9JLeTF0EckmAMgAL81MhKMU4zOrotxtaTFQCA84fEWBxv7wLjFJgv6K4LDGgPFoprmqE3OG8toI71MPIGpEdK6l1WAC3ruMmrvAbQMCdngCRJUp7D/uJa5fVmBojIOXwuAEpMTERZWZnFsbKyMiQmJtq8T0BAAMLDwy2++pO6lva1RuRaCqC9Dmhnof11QOX1rThW1gBJAnIyrAdANc1tTl8Yj9xP/kDuqgssITwQGrVxt/XCSscyiV1pL4K2nAIrqGzErkJjwO7sFnhZqllWS28QyhYyzs4AAe2BnZxVlSTYtV4REXXP5wKgSZMm4fvvv7c49u2332LSpEkeGlHfJ09/JYQHWKzoa14HZC95F/mRSeGdijXlDyutzoDmNn2vxkyeZ08NkFolIXtgJABgy4kKpz12exG08bHjQgMQE6KBQQDbTYX7zu4Ak5lPgRVWNqJVZ6w3ko87k9wJJv+7Cg/0h5/a535tE3mER/8lNTQ0IDc3F7m5uQCA/Px85Obm4tQp4wJ6y5cvx4IFCyzuI5/f0NCA8vJy5Obm4tChQ8rtv//97/HVV1/hueeew5EjR/Doo49i586duPvuu932vLzNiTK5A8zyL9jxA40ZoLyKRlQ0tNp1rW2m6a/Jg2M63RasUcNfbazJYB2Q96s3bZbb3d5XFw6LAwD8eMw5AZAQor0g2BRkS5KkFCELYcyUZDp5SkomZ4Bqmtqw05RtGpYQpiz34EzyFFi5aaNVTn8ROY9HA6CdO3ciOzsb2dnZAIBly5YhOzsbjzzyCADjwodyMCSTz9+1axfWrl2L7OxsXHrppcrtkydPxtq1a/Haa68hKysLH374IT799FOcc8457ntiXkap/4m3XMMkIthfSevbmwXaaqMAGjB+SMkFnKwD8n4Nrd3XAAHAhUONAdC2kxU9XlncXKNWjza9cQrVPCCQp4sA44aVIU7YlNSa0AA/xJgCr+8PG6fbnbUFRkfylhgyFkATOY9rfkPYaerUqV3WgqxZs6bTMXtqR6655hpcc801vRlavyKvAdQxAAKA8elROFpWj50FVfjVObbrqADjBpGFlU1QqyRMtLGrdlSwP8rrW5kB8gH2FEEDwKjkcMSEaFDZqMXuwmrkDOqcHXREtakbKsBPhSCzPcjkOiDAddNfstToYFQ2arHJlNVyRf0PYNwSY3BcCI6ZsrTcB4zIeTiZTO1rAFkJgOQ6IHsWRJSzP1kpETanRbgdhu9ozwB1PS2jUkm4YKgxI7jpeHmvH9e8AFqS2rMjI8wyQCOSXBOQyORpMLmWzdlrAJkzz2xxJ3gi52EA1M81aXXK9gFDrfwVK3dy7S+u7XZHb1vr/5hjK7zvsLcGCGifBtvkhDqgjgXQsqEJoZBni1zVASYbGB1k8b0rAyDzdn5ug0HkPAyA+rmTZ42tydEhGqvp9eTIIAyKC4FBtAc41gghlPV/rBVAyyK5H5jPsHcKDACmmDJAB87UotLOgnpbOq4CLQv0V2P26CQMjA5WlnBwFfOOr4ggf8Sbdr13BfNsFougiZyHAVA/d6LcegG0uQuGGD+8umpjzqtoRFldKzR+KowzbaRqDVeD9g1anQGtOmNBc3dTYAAQHx6IzMQwCAFs7mU7vFwDJK8CbW7V9eOw6YFpFss5uEKqWQA0PCHMYirO2SwyQKwBInIaBkD93PEy2/U/MnsCILn+Z/zAKASaFaZ21N4FxgDIm8n1P4B9U2AAcNEw50yDVXdYBNETzDNArpz+AoD4sAAlO8spMCLnYQDUz3XVASY7b3AMVJIxy1Nc02z1nK7W/zEnT4HVNnMKzJvJ9T8hGrXd+23J6wH9dLy8VyuB25oCc6ekiCD4mZ63s7fA6EiSJPzqnESEaNQYPSDCpY9F1J8wAOrnTp61vgiiufBAf2SlRgIAthzv/Ne7wSDaC6CHdB0ARSlF0MwAeTN7tsHoaEJ6FAL9VThb34ojph3Ue6LK9N7xZEeUWiVhUJxxp/lzXFxwDQBPXDkaux+ZaTH1RkS9wwCoH2vV6VFg2p9paILtDBAATDFNg1mr39hXXIvqpjYEa9QYkxLZ5XW4EKJvsGcbjI4C/NQ4z7QG0E+9aIfvCxkgAHj+2rH4+7wsjDX9ceBqAX62p5aJyHEMgPqxnQXVMAggLMCv2y6W883qgAwddvVevTkfADBzZAL8u9mnSJkCYwbIqznSAm/OGe3wyk7wVoqg3WlUcgSuGpfi0gJoInIdBkD9kBAC//2lEDet2QEAyBkU0+0v8eyBUQjWqFHZqLWYvjhd3YR1+0sAALdOGdTtYytdYNwR3qvZuw1GR3Id0PaCKjRre7YhbnWj54ugicj7MQDqZxpbdbjnvVw8/MkBaHUGTM+Mx9+uHtPt/TR+KuSYtrcw7wZbvbkAeoPA+UNicI4dBZoRQca/2vUGgboWXTdnU1/lyBpA5gbHhWBAZBC0OgN+zre9rlRX+soUGBF5NwZA/YDeIJBf0Yh1+0ow96XN+Cz3DNQqCQ/NzsQbCybYvbbI+R3qgGqb2vDuDuNmtbddONiuawT6q5X9mzgN5r2UDFCAY9NQkiThwmGm95GVgvrutOr0aDRljhgAEVFveHQzVHKdVp0ef/vqKHYUVuNYab2yZxEAJIYH4sXrszHRwdVypwyNA3AY2/Or0KrT4z+/FKJJq0dmYhguHGp7+4uOooL90VyrR3WTFgNj2NXijerkGiAHM0AAMDY1Eu9sL8LJ8gaH7ysvoKmSHM8+ERGZ428QH/X94bN4w1ScDBh3zh6eGIbs1Ej87uKhiAl1fOn+YQmhiA0NQEVDK37Oq8KarQUAgNsuHORQIWhEsAZnalvYCebFGno4BQYAqVHGoLeoqsnh+1Yr+4BpoLJz/SEiImsYAPmo3YXVAIBZoxLwwK8ykR4TYveCdbZIkoQLhsTg09wzWPHZAZTXtyIxPBCXjUl26DpRymKInALzVso6QA52gQFAiikAOl3dDCGEQ8GzXADNXdGJqLdYA+Sj9hTVAAAuGZmIwXGhvQ5+ZBeY2pgLKo1/vd98QTo0fo69jeTaDXlPJ/I+cg1QT/bcSooMhEoCWnUGlDu4MapcAB3N+h8i6iUGQD6oTW/AgeJaAMDYgZFOvfb5Zis9hwX4Yf65Ax2+RoS8IzwzQF6rvhc1QP5qFZIiggAYs0COqDKbAiMi6g0GQD7oSEk9WnUGhAf6ISMmxKnXTooIwmDTFgDX5wx0aCVgmTwFxh3hvVdP2+BlA6J6FgDVKBuhcgqMiHqHAZAPyi0y1v9kpUa6pFD0T5eNxG/GpeDOqfa1vncUGcTtMLxdT7bCMNfTQmh52tTepRuIiGxhEbQP2nOqBoBx9WZXmDY8HtOGx/f4/pHMAHk9uQaoJ0XQAJDSwwxQdR/YCJWIfAMzQD4o11QAne2mTRodpWyH4cIM0N6iGpyqdLzNujt6g8ChM3XQ6gxOv7a3EEKYFUH3NgBy7GfEImgichZmgHxMbVMb8iqMO7xn9dEASP7rvdpFGaCiqiZc/epWpEYH44c/THXqtVd8fgD/+fkUIoL8MWtUAi4bk4zJg2Pg180msL6kSauH3rQhbk+KoAEgNbq9Fd4RLIImImdhAORjck/XAADSYoIR3UfrJCJdnAH6Jb8KbXqBvPJGNLbqENLDaZqODp2pw39/MW79Udvchvd3nsb7O08jOkSD2y8chNsv6llNlLeRsz9qlaRsa+IoOQNUXN0Mg0HYXavGImgicpb+82drP5Er1//00ewP0J4BqmvRQad3/lSSXAQOAMU1jmUYbBFC4C/rDkEIYM7oJLxz63m4IWcgokM0qGrU4umvjqClrWe7m3sbpQU+wM+hRQzNJYYHwk8lQas34Gy9/WsByYXzLIImot5iAORj9pg+/Mf25QAoqP2vd1fsCC/XQAE9227Bmh+OnMXWk5XQ+Knw0OxMTBocg79eORrb/3gxIoP9YRDAibOO723ljXrbAg8AfmoVkiIDAdhfB6Q3CGX1cBZBE1FvMQDyIUII7DV9+I91UQeYM/ipVcqHp7Nb4Vva9DhSUq9874wAqE1vwF/XHwYA3Hx+hlK/Ahify7CEMADA8bP1Vu/va3qzDYa5lEhTK7ydAVBdcxuEsfSIO8ETUa8xAPIhhZVNqG5qg0atwoikME8Pp0vtrfDODYAOFNdCZyrQBYAiB4tsrVn7yynklTciJkSDu6Z1rvMZbgqAjpb2jwxQb7bBMJcabeoEq7LvZyQXQIcF+MG/HxWdE5Fr8LeID5GnfkYmhyPAr2fFqe7S3grv3E4w8+kvoPcZoNqmNqz87hgA4N6Zw6x+6A9LCAUAHC/rLxmgnm+DYc58U1R7yMFyZAinv4io9xgA+ZA9p4z1P9lO3v/LFeROMGe3wssB0MR04xRgbzNAL204juqmNgyND8V1E1OtniNPgR3tNwFQ72uAgPYMkL1TYPJO8Jz+IiJnYADkQ+QP/75cAC2TC6GdPQUmvwaXjUkGAJyuaoIQoot72FbVqMWarQUAgIfnjLC51o8cAJ2ubkZjq/OLuvsap9UAOZgBquYaQETkRAyAfERLmx6HSuoAANmpfbcAWuaKDVErGlpxuroZkgTMHp0IAKhv1SmdQ476Ja8SbXqBYQmhmNrF1h9RIRrEhQUAAI73g04wuQaop/uAyeS1gM7UNCsLK3ZFfq9EswOMiJyAAZCPOFRShza9QHSIRpla6Msigp2/Iaq8BtKQuFDEhwUiNtT4GEV2Ftl2tL2gCgCQkxHT7blyHdCxUt+fBpNrgHo7BZYQFgh/tQSdQaC0rqXb85kBIiJnYgDkI/aYLYDY08Xp3MkVGaCOU4DyFIu9NSYd7Sww1lRNzIju9lx5GuxYP6gDas8A9S4AUqkkDIiUO8G6/xkpiyAyACIiJ2AA5CO8qf4HMOsCa3ZiBkhZAykSgPl+U44HQA2tOhw8UwugvaC6K/2pENpZRdBA+8/InmJ1pQiaXWBE5AQMgHyAEAI78o3TNdl9eAFEcxHyhqiNzskAGQxmi0CagsBUU41JT6bAdhdWwyCMdSpJEd1PKSqLIZb5fg1QnVIE3ftAxJFd4TkFRkTOxADIB5w424DSuhZo/FSYYEe2oi+QM0A9LVDuKK+iAfWtOgT6q5SFCduzC45ngHaY6n/OTe9++gsAhppqgErrWpz2nPqqBifVAAFm05R2BKncCJWInIkBkA/46XgFAOOHdWAPd+d2N7kN3llF0LlFxumq0QMilHb1VOXDtecB0AQ7A6DwQH8kRxj3tvL1BRGd1QYP2J8BatXpUVZvLJRmDRAROQMDIB+w+YQxAJoyNNbDI7Gf/CHWpNWjVdf7XdRzrWwCq2y1UN3s0FpAWp1BKSo/N8P+jNqwxP5RB+SsrTAA+9cCemXjSdQ0tSE2VIMh8aG9flwiIgZAXk6rM+DnvEoAwAVeFACFBfpBZWpWq3VCJ1h7EXh7wJIcGQSVBLTqDCivb7X7WvuLa9GqMyA6RIPBcfZ/2PaHOiCd3oAmrTFg7e1WGEB7kFpS24w2vcHqOcfL6rFqwwkAwKO/HuU1WU4i6tsYAHm53aeq0aTVIzZUgxGJ4Z4ejt1UKgkRyjRY1wGQzsYHo8x8B/ixZtuA+KtVSgGzI3VAyvRXWpRDSwoonWA+vBZQY2t7ts4ZU2BxoQEI8FPBIIDS2s5rARkMAg9+tA9teoEZI+IxZ3RSrx+TiAhgAOT1fjpeDgA4f0gsVKq+v/6PuSg7FkP8548nkfXYN/jxWLnNc+Qd4OPCApQ6HFlKDzrB5I66iXbW/8iUxRBdNAWm0xtwx9u7sOz93B5v79FbdaYC6AA/FTR+vf/1IUkSBig/o85B6n9+KcTuUzUI0ajx58vP8Yo1rojIOzAA8nKbj8v1P3EeHonjIu1YDPHHY+Vo1OrxwId7bXZXydNfWSmdF4FMcbAQ2mAQ2Flo/wKI5obEh0KSgMpGLSoa7J9ys9fXB8vw1cFSfLy7WGlFdzdnbYNhLtVGHdCZmmY8/eURAMCDszORHNn3VzgnIu/BAMiL1TRpsa/Y2P10wRDvqf+R2ZMBqmo03lZW14on1h3udHtDqw4f7DwNAMg2m/6SObrj+PGzDahtbkOQvxqjkh2bUgzW+Ckf5s7OAgkh8PpPecr3dR5qtXfmIogya51gQgj836cH0KjVY9zASPw2J81pj0dEBDAA8mpbTlRCCGBofCgSO0z9eIOoEGMAJAc51pjf9t7OImXKDwDa9Abc+Z9dOFpWj9hQDX4zLqXT/W1lF2yR9//KHhgJfxu7v3elN4XQr/54Eq9tOmn1tl2F1UqmC3DuHmqOaGh13hpAso6rQRsMAs9/ewzfHzkLf7WEp38zxuumd4mo72MA5MU2nzAGA944/QUAMabNSisbrH+YCyGUD/pZoxIAAA99tB8NrToIIfDwJ/vx0/EKBPmrsXrRRKtBoKOLIfa0/kcm1wE52gpfWtuCp748gifWH8G6fSWdbjfP/gDO3UPNEc5cA0hmngGqbW7Drf/eiX/8YOz6+sMlwzHUFFQSETmT836LkVsJIbDpmPet/2MuRskAWa+XqW/VoU1vLPZ98qoxOFSyGUVVxrqQ2NAAvL/zNFQS8NL12RiTEmn1GvIU2JmaFuj0BmWRRFt2yitAO1j/IxueKGeAHAuATpxtzxg98tkBnDcoGjGhAQCAgopGfHOoDACQHBGIM7UtqPGpKTB52rABl7+0GQWVTdD4qfDElaNx9fjOWT0iImdgBshLFVQ2obimGf5qCTmDevZh7WnRIcYP+EobU2BVpsxQsEaN6BANnrpqDADg7Z8L8fx3xwAAj19xDi4ekWDzMRLCAqFRq6A3CJRYabM2d7q6CWdqW6BWSVbriewxNL69Fd6RTq38ivYAqLJRi0e/OKR8/+aWfAgBTB0epwR6tR6aAqt34j5gMnnPttrmNhRUNmFAZBA+umMygx8icikGQF5qs6kWZnxaFII13pnIi+mmBqjK9CEfbTrv/CGxmH/uQOX2O6cOxg3dFMeqVGZt1t1Mg8nr/5yTHN7j13RQXAjUKgl1LTqcdWDxxZPljQCMxexqlYQv9p7B1wdLUdOkxfumIu9bpwyyq3POleqduA+YLDpEg3DT9c4fEoMvll6A0SkRTrs+EZE13vnJSdjkxe3vsujuAiBTBkgOlADgj5dmorZZi4zYEPxh5nC7HiclKgj5FY04XdUMDLZ93o9HjUFlT+t/ACDQX420mGDklTfiaGk9EsLtK07PqzAGQJeNScLolAi8svEkHv7kAK4aNwDNbXqMSArH5MEx2GQKfD01Bda+DYbzfnVIkoQXrx+HU1VNmD8xtdtpSiIiZ2AA5IV0egN+Pmnc/sJb63+A9gCoskELIUSnNXzkwCjKLAAKC/THyzeMd+hx7CmEPlnegM/3ngEAXDqmd6sND08IQ155I46V1ePCYfYFqHnlximwQXGhuCJ7AL49VIYTZxvw2iZj8fOtUzIgSZJdSwe4kjIF5sQACAAusvN1IiJyFv6p5YX2nq5FfasOUcH+GJXsvVMFcheYVm9QMgvmOk6B9ZQ9u8I//+0xGARwcWY8xg20fwNUa9JjQ7p9PHMtbXoU1xhbwAfFhSDQX41nrh6j7JWWEB6Ay8YkAwAiTduHOGP/tJ5oL4J2Xg0QEZEneDQA2rRpE+bOnYvk5GRIkoRPP/202/ts3LgR48aNQ0BAAIYMGYI1a9ZY3P7oo49CkiSLr8zMTNc8AQ+RswXnDIiA2ovXRwnW+CHItLGltWkw+VhMbwMgZTFE62sBHTxTi/+ZWs//cIl902pdaW/rtm/toYLKRghhnFaSn+u4gVG4c6pxvu6uqUOUbSeUGiCPdYEZH9eZbfBERJ7g0QCosbERWVlZWLVqlV3n5+fnY86cOZg2bRpyc3Nxzz334JZbbsHXX39tcd6oUaNQUlKifG3evNkVw/eYsjpjN1OinfUlfZkyDWYlAJLXB4rqZQDU3XYYz31j7Cibm5WMkQ6u/mzNANOWDXJWpzv5pgLojLhQi2nA+y4Zji0PTceCSe2F3hFBxteixmMLITq/DZ6IyBM8+lts9uzZmD17tt3nv/rqq8jIyMBzzz0HABgxYgQ2b96M559/HrNmzVLO8/PzQ2JiotPH21fI7dxJXrj6c0cxoRoU1zQrBc/m5DqXXmeATBmZs/WtaGnTI9CUdQKAXYVV+OHIWahVEu6dMbRXjyNLMVt92lptU0dyAfRg09SZTJIkJZiSyRkgW/uiuZor1gEiIvIEr6oB2rZtG2bMmGFxbNasWdi2bZvFsePHjyM5ORmDBg3CDTfcgFOnTnV53dbWVtTV1Vl89WVyBijBBwKgrjrB5KyQvF5Qbx4jWGMMesyzMkIIPPPVUQDANeNTMCgutFePI5ODloZWHeqau9+09KRSAB3SzZnt+6fVNLV5ZEd4V2yGSkTkCV4VAJWWliIhwXLRu4SEBNTV1aG52fjBlpOTgzVr1uCrr77CK6+8gvz8fEyZMgX19bZX5n3yyScRERGhfKWmprr0efRWqQ9OgVVYWQ1aXiE6OqR3H7aSJCmF0HuLatCmNwAANp+owC/5VdCoVVh6sXOyPwAQpFErWavTNd0XQueZpsDsCcDkDJDOIKwWjruSEII1QETkM3zut5j5lNqYMWOQk5ODtLQ0vP/++1i8eLHV+yxfvhzLli1Tvq+rq+vTQVBprTEwsHeNmb5MWQzR2hRYo/HDtrcZIMDYCn+0rB7L3t+LBz/ah7SYEOXD/IbzBnaaauqtlKggVDZqcbq6uctOPSGEWQt89xmgQH81AvxUaNUZUNPU5tZMTKvOoGxNwikwIvJ2XvVbLDExEWVlZRbHysrKEB4ejqAg6x9gkZGRGDZsGE6cOGHzugEBAQgI6P2HrDtodQZUmjIj3rgDfEfyflcdp8BadXolwxEd3LsaIABYfEEGzta34HhZA5rb9MreW8EaNe6aOqTX1+9oQFQQ9p6uRXE3nWCVjVrUteggSUB6TPcBEGDMApXVtaK2uQ3uDNPl+h9JAkK8dPVxIiKZV/0WmzRpEtavX29x7Ntvv8WkSZNs3qehoQEnT57EjTfe6OrhucXZ+hYIAfirJacEBp5mqwtMDoj8VBLCg3r/Np00OAaf330BDAaB0roWnDjbgPyKRoxJiUBcmPODX/NC6K7kmwqgkyOCLIqzuxIZpEFZXavbt8OQA9JQjR9UXrz8AhER4OEAqKGhwSIzk5+fj9zcXERHR2PgwIFYvnw5iouL8e9//xsAcMcdd+Cll17CAw88gJtvvhk//PAD3n//faxbt065xn333Ye5c+ciLS0NZ86cwYoVK6BWqzF//ny3Pz9XUAqgwwN94kPI1n5g5qtAd9dF5QiVSkJyZBCSI4PsXqW5J9pb4buuAXJk+ksm1wG5ezVopf6H019E5AM8+pts586dmDZtmvK9XIezcOFCrFmzBiUlJRYdXBkZGVi3bh3uvfdevPDCC0hJScEbb7xh0QJ/+vRpzJ8/H5WVlYiLi8MFF1yAn3/+GXFxvrHUvlz/4wsF0IDtLjD5e2/Nctm7GKJcAD3YgQ40Ty2GyBZ4IvIlHv1NNnXq1C5beTuu8izfZ8+ePTbv8+677zpjaH1WqQ+1wANAjKnAuaKh1WLNHCUA6uUaQJ4i70Df3WKIJ5UOMAcyQKbFEGvdngEyTYGxA4yIfIBXtcGTb60CDQDRpv3AWnUGNGn1ynElAAr10gDINAVW09TWZbt6XoVpCiy2BxkgN9cAyVNgXAOIiHwBAyAvU1rrWwFQiMbY1g1YToN5+xRYWKA/Ikwbl9rqBNPpDThVaawRciQDFOGhKTBug0FEvoQBkJdRAiAfmQKTJEkphK60FgB56RQY0H0hdFF1M3QGgUB/lUMBbftq0O6dApPrmcKDmAEiIu/HAMjLKKtA+0gABLRPc1WZrQat7ATvpVNgQPeF0HIHWEZsqEMdfZFB7p8C0+kN+HzvGQDAhUNj3fa4RESuwgDIiwghfGobDJm80nOl2WrQcjYoykunwACzQmibAZDjBdCAZ6bAfjxWjvL6VkSHaDA9M6H7OxAR9XEMgLxITVMbtDrjPlbx4d6xcrU9rE2BVTc6Zyd4T1IWQ7TRCSYXQHfcBb47cheYOzNA7+8sAgBcmT0AGj/+2iAi78ffZF5Ezv5Eh2gQ4GffqsHewNpaQOYLIXoruQbI1hTYSQc2QTUnd4HVNmvdsiN8RUMrvj98FgAwb0Lf3SOPiMgRDIC8iK91gMnkOh95CsxgEMoqx96dAep6CkzeBiPDwQyQPC3YphcWSwe4yqd7iqEzCGSlRGB4YpjLH4+IyB0YAHkRXyyABsy3wzAWQdc2t8FgSmx4cwZIDoAqGlrR0mYZqNS3tKG83vh8Ha0BCvRXKdNQrt4OQwiB93YYp7+uYfaHiHwIAyAvImeAEnwsAyQXQcvTXnItUFigH/zV3vsWjQjyR4jGOFXZcUVouQA6LizA4YUFJUlyWyfY3tO1OH62AQF+KszNSnbpYxERuZP3frr0Q762CrSs447wvjD9BRgDFVu7wrevAO1Y9kfWXgfk2gBILn6efU6isrAjEZEvYADkRdqnwHynAwww6wIz1QDJ//XmRRBltlrhT57tWQG0LDLY9Z1gzVo9vsg1rv3D4mci8jUMgLyIz06BmYqgm9v0aNbqfWIVaJlSCN1hNegfjhi7qsakRPTousoUWLPraoC+OliC+lYdUqKCcN6gGJc9DhGRJzAA8iJyBigpIsjDI3GusAA/aEy1PpWNrcoUmC8EQNZa4Y+W1uNQSR381RJmn5PYo+u6Y0PU93ecBgBcMz7VoZWqiYi8AQMgL9HSplc+7HytBkiSJIu1gNqnwLx/qk+uATKfAvtkTzEAYNrweGUqy1GRLt4PrK6lDT/nVwIAfjN+gEseg4jIkxgAeQm5ADrQX4XwIN/bjdu8EFpuh48O8f6i2wEd9gMzGAQ+yzUGQFeN63lgEeHiLrD9p2shhHEKTw7iiIh8ie99kvoo80UQJcn3piPkxRCrGrSoMn2o+0IGSJ4CK6tvgVZnwM7CKpTUtiA80A9Th8f3+LrKjvAu6gLLLaoBAGSlRrrk+kREnsYMkJeQ6398rQBa1p4BavWpDFBsqAYBfioIAZTUNuOT3cbsz5wxSQj07/l2JkobvIsyQPtO1wAAxqZEuuT6RESexgDIS8gZoCQfWwVaZjEF5kM1QJIkKdNgJ8sb8OWBUgDAldkpvbqu3AXWcSVog0Hgxn/9ghve+Bl6Q8/3CdtbVAug511qRER9HQMgL6FkgHw0AIoNNa0G3aBFlY8shCiTa2jWbC1EQ6sOAyKDMCEtqlfXjJC7wDpMgRVUNuKn4xXYcqISx8rqe3TtsroWlNa1QCUB5wxgAEREvokBkJfw1VWgZXIGqLimGS1tBgDevQ+YObkOaNOxcgDAFdnJvW4rl7vAapvaLHaEP3imTvn/XYXVPbr2XlP9z7CEMIQEsEyQiHwTAyAv4as7wcvkAOj4WeMWERo/lbKPlreTF0OUXZnd+7byKFMGSKs3oNlso1XzAGh3TwMgU/0Pp7+IyJcxAPISZXXGwmBfnQKTp7vkHdJjQjQ+0+1mHgCNHhCBIfFhvb5mkL9aWTzSvBX+4Jla5f93neppBsh4DXaAEZEvYwDkBQwG0W+mwGRRPVwgsC+Sp8AA4AonZH8AY3G1XAckF0ILIXDILANUWNmEioZWh65rMAilAyyLHWBE5MMYAHmBisZW6AwCKgmIC/P+zihrYjp0fMnrAvmCtJgQqCRArZIwNyvJadeVO8HkVviyulZUNmqhVklIizEWXjs6DVZQ2Yi6Fh0C/FQYntj7TBURUV/FCkcvUFZr/Cs+NjQA/mrfjFnDg/zgp5KgM7Vu+8I+YLK4sAC8cF02QgLUiA9zXgYvskMnmDz9NTguBOPTolBY2YRdp6pxySj79xuT639GJYf77HuNiAhgBsgryC3wiT5a/wNY7gcG+NYUGADMzUrG9MwEp16zfT8wOQAyTn+NSo7AuIHGNntHM0Cs/yGi/oIBkBfw9VWgZeYBkK+sAeRK8hRYTbOxBkjOAI1KDsc40zpD+07XQqsz2H1NOQM0lgEQEfk4BkBeoMzHW+Bl5nU/vrIGkCspU2BWMkCDYkMQGeyPVp0Bh0rqbF7DnFZnUK4xhgXQROTjGAB5gZJa358CAyy3vmAGqHvtU2Ba1Da1KTvOj0wOhyRJGG+aBrN3QcSjpfXQ6gwID/RDegx3gCci38YAyAv4egu8zDzo8aUiaFeJCGrPAB0sMU5/pUYHKcflaTB764Dk6a+s1EifWYOJiMgWBkB9nN4gcKTUuKfTgA4rCvuaaAZADpELxWua25T1f0Ylta/erBRC27kgorwFBtf/IaL+gAFQH/dzXiUqGloRGeyvfKD5KvMaIAZA3ZNrgGqb2szqf8KV27NSI6BWSSipbcGZmuZur2eeASIi8nUMgPq4z3PPAABmn5MEjZ9v/7jkKTBJaq9vIdsizLrAlA6wAe0BULDGDyOTjN93VwfU0KpT9mHL4h5gRNQP+PYnqpfT6gz48kAJADh1BeG+Si6Cjgzyh7qXu6X3B3IGqKpRi5PljQCMHWDmxqfZNw12oLgWQgBJEYGI9/FaMyIigAFQn7bpWDnqWnSIDwtATkaMp4fjcqMHRGBsaiSunTjQ00PxCnKWrE0voDcIxIZqEN9hq5TsgZEAui+E/u5QGQDW/xBR/8GtMPqwz/cap78uG5PcLzIiQRo1Pl1yvqeH4TVCNGr4qyW06Y3bh4xMjujUvSVngA6eqUOzVo8gjbrTdX48Vo5/bckHAFw+NtnFoyYi6huYAeqjmrQ6fGv6q7w/TH+R4yRJQkRQe62UeQG0bEBkEBLCA6Az2+Xd3JmaZtzz7h4IAVyfMxCzR/O9RkT9AwOgPur7w2fR3KbHwOhgbktANsl1QID1AEiSJCUL9PpPeahq1Cq3aXUGLFm7G9VNbThnQDgeuWyk6wdMRNRHMADqo+Tpr7lZSVyUjmyS9wMDOhdAy67KTgEAfHf4LKY/txFrfzkFg0HgyS8PY8+pGoQH+uGVG8Yj0L/z9BgRka9iANQH1Ta34cej5QCMu4gT2SJngEID/JAWbX37ihkjE/DRnZMwIikcNU1t+OMn+3HJyk14c0sBAODv88Yi1cZ9iYh8FQOgPujrg6XQ6g0YlhCKzMTO0xpEMrkTbERSGFRdFMqPT4vGF3efj0cuG4nQAD+cMK35c+fUwZgxMsEtYyUi6kvYBdYHfWGa/vo1sz/UjThT2/voAZHdnuunVuHmCzIwZ0wSXvj+OPxVEv4wc5iLR0hE1DcxAOpjyutbseVEBQBj+ztRVxZMSoNaknDjpDS775MQHognrhztwlEREfV9DID6mB+PlcMgjIsCpseGeHo41MclRQThvlnDPT0MIiKvwxqgPmaPacuCSYN9f+VnIiIiT2EA1MfsOVUDAMjm2j9EREQuwwCoD2ls1eFIaR0AYJxp8ToiIiJyPgZAfci+07UwCCA5IhAJ3JGbiIjIZRgA9SF7ioz1P9kDmf0hIiJyJY8GQJs2bcLcuXORnJwMSZLw6aefdnufjRs3Yty4cQgICMCQIUOwZs2aTuesWrUK6enpCAwMRE5ODrZv3+78wbvA7sIaAED2wEiPjoOIiMjXeTQAamxsRFZWFlatWmXX+fn5+ZgzZw6mTZuG3Nxc3HPPPbjlllvw9ddfK+e89957WLZsGVasWIHdu3cjKysLs2bNwtmzZ131NJxCCIFcZoCIiIjcQhJCCE8PAjDuWv3JJ5/giiuusHnOgw8+iHXr1uHAgQPKseuuuw41NTX46quvAAA5OTmYOHEiXnrpJQCAwWBAamoqli5dioceesiusdTV1SEiIgK1tbUID3fPVhRFVU2Y8swG+Ksl7H90FjemJCIicpAjn99eVQO0bds2zJgxw+LYrFmzsG3bNgCAVqvFrl27LM5RqVSYMWOGco4nGQwCG46exd6imk637Tat/zMyOYLBDxERkYt5VQBUWlqKhATLjRsTEhJQV1eH5uZmVFRUQK/XWz2ntLTU5nVbW1tRV1dn8eUKqzacwE1v7sDfvz3W6TZ5/Z9xrP8hIiJyOa8KgFzlySefREREhPKVmprqkse5fOwASJJxu4tjZfUWt8krQLP+h4iIyPW8KgBKTExEWVmZxbGysjKEh4cjKCgIsbGxUKvVVs9JTEy0ed3ly5ejtrZW+SoqKnLJ+AfGBGPWSOM4Vm/OV463tOlx8Iwx68QVoImIiFzPqwKgSZMm4fvvv7c49u2332LSpEkAAI1Gg/Hjx1ucYzAY8P333yvnWBMQEIDw8HCLL1e5ZUoGAODjPcWobGgFABworoXOIBAbGoCUqCCXPTYREREZeTQAamhoQG5uLnJzcwEY29xzc3Nx6tQpAMbMzIIFC5Tz77jjDuTl5eGBBx7AkSNH8PLLL+P999/Hvffeq5yzbNkyvP7663jrrbdw+PBh3HnnnWhsbMRNN93k1udmy/i0KGSlRECrM+A/Pxufp3n9jyRJHhwdERFR/+DRAGjnzp3Izs5GdnY2AGPwkp2djUceeQQAUFJSogRDAJCRkYF169bh22+/RVZWFp577jm88cYbmDVrlnLOtddei2effRaPPPIIxo4di9zcXHz11VedCqM9RZIkLJ4yCADw9s8FaGnTcwVoIiIiN+sz6wD1Ja5eB6hNb8BFz2zAmdoWPHP1GPz9m2MorWvBu7edh/MGxTj98YiIiPoDn10HyFf4q1VYODkdALDyW2Pwo1ZJGJMS4dmBERER9RMMgDzkunMHIlijxpnaFgBAZmIYgjV+Hh4VERFR/8AAyEMigvwxb0L7ekPcAJWIiMh9GAB50E3np0Nu+spOZQE0ERGRuzAA8qC0mBDccdFgjEoOx8Uj4j09HCIion6DXWBWeGI3eCIiIuoddoERERERdYEBEBEREfU7DICIiIio32EARERERP0OAyAiIiLqdxgAERERUb/DAIiIiIj6HQZARERE1O8wACIiIqJ+hwEQERER9TsMgIiIiKjfYQBERERE/Q4DICIiIup3GAARERFRv+Pn6QH0RUIIAEBdXZ2HR0JERET2kj+35c/xrjAAsqK+vh4AkJqa6uGREBERkaPq6+sRERHR5TmSsCdM6mcMBgPOnDmDsLAwSJLk1GvX1dUhNTUVRUVFCA8Pd+q1qR1fZ/fg6+wefJ3dg6+z+7jqtRZCoL6+HsnJyVCpuq7yYQbICpVKhZSUFJc+Rnh4OP+BuQFfZ/fg6+wefJ3dg6+z+7jite4u8yNjETQRERH1OwyAiIiIqN9hAORmAQEBWLFiBQICAjw9FJ/G19k9+Dq7B19n9+Dr7D594bVmETQRERH1O8wAERERUb/DAIiIiIj6HQZARERE1O8wACIiIqJ+hwGQC6xatQrp6ekIDAxETk4Otm/f3uX5H3zwATIzMxEYGIjRo0dj/fr1bhqpd3PkdX799dcxZcoUREVFISoqCjNmzOj250JGjr6fZe+++y4kScIVV1zh2gH6CEdf55qaGixZsgRJSUkICAjAsGHD+LvDDo6+zitXrsTw4cMRFBSE1NRU3HvvvWhpaXHTaL3Tpk2bMHfuXCQnJ0OSJHz66afd3mfjxo0YN24cAgICMGTIEKxZs8bl44Qgp3r33XeFRqMRq1evFgcPHhS33nqriIyMFGVlZVbP37Jli1Cr1eKZZ54Rhw4dEn/605+Ev7+/2L9/v5tH7l0cfZ2vv/56sWrVKrFnzx5x+PBhsWjRIhERESFOnz7t5pF7F0dfZ1l+fr4YMGCAmDJlirj88svdM1gv5ujr3NraKiZMmCAuvfRSsXnzZpGfny82btwocnNz3Txy7+Lo6/zf//5XBAQEiP/+978iPz9ffP311yIpKUnce++9bh65d1m/fr14+OGHxccffywAiE8++aTL8/Py8kRwcLBYtmyZOHTokHjxxReFWq0WX331lUvHyQDIyc4991yxZMkS5Xu9Xi+Sk5PFk08+afX8efPmiTlz5lgcy8nJEbfffrtLx+ntHH2dO9LpdCIsLEy89dZbrhqiT+jJ66zT6cTkyZPFG2+8IRYuXMgAyA6Ovs6vvPKKGDRokNBqte4aok9w9HVesmSJmD59usWxZcuWifPPP9+l4/Ql9gRADzzwgBg1apTFsWuvvVbMmjXLhSMTglNgTqTVarFr1y7MmDFDOaZSqTBjxgxs27bN6n22bdtmcT4AzJo1y+b51LPXuaOmpia0tbUhOjraVcP0ej19nf/85z8jPj4eixcvdscwvV5PXufPP/8ckyZNwpIlS5CQkIBzzjkHTzzxBPR6vbuG7XV68jpPnjwZu3btUqbJ8vLysH79elx66aVuGXN/4anPQW6G6kQVFRXQ6/VISEiwOJ6QkIAjR45YvU9paanV80tLS102Tm/Xk9e5owcffBDJycmd/tFRu568zps3b8a//vUv5ObmumGEvqEnr3NeXh5++OEH3HDDDVi/fj1OnDiBu+66C21tbVixYoU7hu11evI6X3/99aioqMAFF1wAIQR0Oh3uuOMO/PGPf3THkPsNW5+DdXV1aG5uRlBQkEselxkg6neeeuopvPvuu/jkk08QGBjo6eH4jPr6etx44414/fXXERsb6+nh+DSDwYD4+Hi89tprGD9+PK699lo8/PDDePXVVz09NJ+yceNGPPHEE3j55Zexe/dufPzxx1i3bh0ef/xxTw+NnIAZICeKjY2FWq1GWVmZxfGysjIkJiZavU9iYqJD51PPXmfZs88+i6eeegrfffcdxowZ48phej1HX+eTJ0+ioKAAc+fOVY4ZDAYAgJ+fH44ePYrBgwe7dtBeqCfv56SkJPj7+0OtVivHRowYgdLSUmi1Wmg0GpeO2Rv15HX+v//7P9x444245ZZbAACjR49GY2MjbrvtNjz88MNQqZhDcAZbn4Ph4eEuy/4AzAA5lUajwfjx4/H9998rxwwGA77//ntMmjTJ6n0mTZpkcT4AfPvttzbPp569zgDwzDPP4PHHH8dXX32FCRMmuGOoXs3R1zkzMxP79+9Hbm6u8vXrX/8a06ZNQ25uLlJTU905fK/Rk/fz+eefjxMnTigBJgAcO3YMSUlJDH5s6Mnr3NTU1CnIkYNOwW00ncZjn4MuLbHuh959910REBAg1qxZIw4dOiRuu+02ERkZKUpLS4UQQtx4443ioYceUs7fsmWL8PPzE88++6w4fPiwWLFiBdvg7eDo6/zUU08JjUYjPvzwQ1FSUqJ81dfXe+opeAVHX+eO2AVmH0df51OnTomwsDBx9913i6NHj4r//e9/Ij4+XvzlL3/x1FPwCo6+zitWrBBhYWHinXfeEXl5eeKbb74RgwcPFvPmzfPUU/AK9fX1Ys+ePWLPnj0CgPj73/8u9uzZIwoLC4UQQjz00EPixhtvVM6X2+Dvv/9+cfjwYbFq1Sq2wXurF198UQwcOFBoNBpx7rnnip9//lm57aKLLhILFy60OP/9998Xw4YNExqNRowaNUqsW7fOzSP2To68zmlpaQJAp68VK1a4f+BextH3szkGQPZz9HXeunWryMnJEQEBAWLQoEHir3/9q9DpdG4etfdx5HVua2sTjz76qBg8eLAIDAwUqamp4q677hLV1dXuH7gX2bBhg9Xft/Jru3DhQnHRRRd1us/YsWOFRqMRgwYNEm+++abLxykJwTweERER9S+sASIiIqJ+hwEQERER9TsMgIiIiKjfYQBERERE/Q4DICIiIup3GAARERFRv8MAiIiIiPodBkBE5HMWLVqEK664wtPDIKI+jJuhEpFXkSSpy9tXrFiBF154gXs1EVGXGAARkVcpKSlR/v+9997DI488gqNHjyrHQkNDERoa6omhEZEX4RQYEXmVxMRE5SsiIgKSJFkcCw0N7TQFNnXqVCxduhT33HMPoqKikJCQgNdffx2NjY246aabEBYWhiFDhuDLL7+0eKwDBw5g9uzZCA0NRUJCAm688UZUVFS4+RkTkSswACKifuGtt95CbGwstm/fjqVLl+LOO+/ENddcg8mTJ2P37t245JJLcOONN6KpqQkAUFNTg+nTpyM7Oxs7d+7EV199hbKyMsybN8/Dz4SInIEBEBH1C1lZWfjTn/6EoUOHYvny5QgMDERsbCxuvfVWDB06FI888ggqKyuxb98+AMBLL72E7OxsPPHEE8jMzER2djZWr16NDRs24NixYx5+NkTUW6wBIqJ+YcyYMcr/q9VqxMTEYPTo0cqxhIQEAMDZs2cBAHv37sWGDRus1hOdPHkSw4YNc/GIiciVGAARUb/g7+9v8b0kSRbH5O4yg8EAAGhoaMDcuXPx9NNPd7pWUlKSC0dKRO7AAIiIyIpx48bho48+Qnp6Ovz8+KuSyNewBoiIyIolS5agqqoK8+fPx44dO3Dy5El8/fXXuOmmm6DX6z09PCLqJQZARERWJCcnY8uWLdDr9bjkkkswevRo3HPPPYiMjIRKxV+dRN5OElwulYiIiPoZ/hlDRERE/Q4DICIiIup3GAARERFRv8MAiIiIiPodBkBERETU7zAAIiIion6HARARERH1OwyAiIiIqN9hAERERET9DgMgIiIi6ncYABEREVG/wwCIiIiI+p3/B3Z8bZgZxSdLAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from pygbm import GBMSimulator\n", "simulator = GBMSimulator(y0=1.0, mu=0.05, sigma=0.2)\n", "t_values, y_values = simulator.simulate_path(T=1.0, N=100)\n", "simulator.plot_path(t_values, y_values, output=\"gbm_plot.png\")\n" ] }, { "cell_type": "markdown", "id": "041a3084-6e57-4dca-9a5a-84bfe4f80247", "metadata": {}, "source": [ "Can also implement in the command line:" ] }, { "cell_type": "code", "execution_count": 3, "id": "5f567e6a-cfb8-4844-b8e3-193dfaefbb47", "metadata": {}, "outputs": [], "source": [ "!pygbm simulate --y0 1.0 --mu 0.05 --sigma 0.2 --T 1.0 --N 100 --output gbm_plot.png\n" ] }, { "cell_type": "markdown", "id": "4127d117-7d18-49b9-a914-9f56492f7bcb", "metadata": {}, "source": [ "One can check usage of this by typing:" ] }, { "cell_type": "code", "execution_count": 4, "id": "82179453-97d1-4804-a90f-7a5150a98632", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "usage: pygbm [-h] {simulate} ...\n", "\n", "Simulate Geometric Brownian Motion\n", "\n", "positional arguments:\n", " {simulate}\n", " simulate Simulate a GBM path and plot it\n", "\n", "options:\n", " -h, --help show this help message and exit\n" ] } ], "source": [ "!pygbm --help\n" ] }, { "cell_type": "code", "execution_count": null, "id": "fe59b7cd-9c01-4248-a3d9-c6d72c145a87", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python (pygbm)", "language": "python", "name": "pygbm" }, "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.12.3" } }, "nbformat": 4, "nbformat_minor": 5 }