You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
311 lines
90 KiB
311 lines
90 KiB
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 11,
|
|
"metadata": {
|
|
"scrolled": true
|
|
},
|
|
"outputs": [
|
|
{
|
|
"name": "stdout",
|
|
"output_type": "stream",
|
|
"text": [
|
|
"N = 10100\n",
|
|
"ion has left the boundary!\n",
|
|
"n = 1\n",
|
|
"simulation time = 9.000900090009e-09\n",
|
|
"number of positiions = 2\n"
|
|
]
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAQYAAADtCAYAAABdwdvJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABowElEQVR4nO29d3wc5bk9fmZ3pZW0klZdVi+WXGRLsiwXHEoCmBJiXCjGJJdyDZcSTCBwk0tCCNxvSAK5hMAnQHIT8kvITUgAY4gBBwKkECC4q8vqvW3X9v7+/pDeYbbP7M5KMtrz+egDXq1mZ3fnPfM+z3Oe8zCEECSQQAIJcCFZ7BNIIIEElh4SxJBAAgkEIEEMCSSQQAASxJBAAgkEIEEMCSSQQABkEX6fKFkkkED8wSz2CfgjsWNIIIEEApAghgQSSCAACWJIIIEEApAghgQSSCAACWJIIIEEApAghgQSSCAACWJIIIEEApAghgQSSCAACWJIIIEEApAghgQSSCAACWJIIIEEApAghgQSSCAACWJIIIEEApAghgQSSCAACWJIIIEEApAghkUAIQQulwsejwcJl+4EliIiGbUkIDK8Xi9cLhdsNhsAgGEYyGQyJCUlQSqVQiKRgGGWnG9HAssMTIQ7VuJ2JhIIIfB4PHC5XAAAt9vNPu71etnnUVKQy+VISkpKEMXywJL7ghM7hgUAN3TwX+QMw0Aqlfo8t729HVVVVUhLSwPDMEhKSoJMJoNMJgPDMAmiSCDuSBBDnOH1euF0OkEIYRf11NQU+vv7kZqaiuzsbGRnZ0OhULC/l0gkkEqlkEqlIITA6XTC4XAAACQSCZKSktjQI0EUCcQDCWKIEwghcLvdcLvd7GL3eDw4c+YMHA4Hmpqa4Ha7odfrMTw8DIvFAoVCgezsbLjdbjYpyd1R0MecTiecTieAOaLwz1EkkECsSOQY4gB6l/d6vewd3Ww2o729HSUlJSgtLYXL5fK50xNCYLFYYDAYMDIyAolEgszMTHZHkZKSErAzoN8dN0eRIIqzEktuy5cgBpFBE4w0dACAiYkJjI6Oor6+HhkZGSxxAAgaBnR2dqKiogKEEOj1euj1ejgcDqSnp/sQhT/8iWJsbAwVFRUJolj6WHLEkAglREKw0MHlcqGzsxMymQxbtmyBTDb3cdNdRChSpmSRkZGBjIwMlJeXw+v1wmw2Q6/X48yZM3A6ncjMzERWVhays7Mhl8vZv6Ohh1qtRnl5eUDoQZOZCaJIIBQSxCACvF4vJicnoVAokJqaCoZhMDs7i46ODlRVVaG4uDjm16ChRWZmJioqKuD1emEymaDX69HV1QWXywWlUons7GxkZWUhOTmZ/TsKQggIIXA4HGwyUyqVsrsJWvVIIIEEMcQArjZBrVZDIpEgNTUVQ0NDmJmZQVNTE9LS0gQfN9xugkIikUCpVEKpVKKyshIejwdGoxF6vR7j4+PweDyw2+1Qq9XIyspCUlJSQAWDaijsdjv7GCUKuqNIEMXyRIIYooS/NoGGDqdOnYJCocCWLVsWdJsulUrZ/AMwl+s4duwYjEYjRkdHQQhBVlYW+xNME5EgigQoEsQQBYJpE+x2O3p6elBXV4f8/PyYjs9nxxAJNDRYuXIlgDml5ezsLFseZRiGzU8olcqgmghKFKdPn0ZlZSVSU1MTRLFMkCAGAeCGDnSX4PV60dfXh9nZWaxcuTJmUgCCVypihUwmQ25uLnJzcwEALpcLBoMBGo0GAwMDkEqlLFFkZmb6EAUhhE1Uer1e2Gw2n0Rngig+e0gQA08E0ybYbDa0tbUhLy8PxcXFSEpKEvX14omkpCTk5+ezROZ0OmEwGKBSqdDf3w+ZTIasrCzk5OT4iK0oIdJzTBDFZxMJYuCBYKHDzMwM+vv7UVdXh+zsbAwMDIi2mMUIJYQiOTkZBQUFKCgoAAA4HA7o9XpMTU1Br9eju7sbubm5yM7ORnp6OtvcFYkoaI9HgijOLiSIIQxCyZp7enpgt9uxefNmn7LgZ8lbQS6XY8WKFVixYgWcTicqKythsVgwPj4Os9kMuVzOJjvT09NZkvAnCo/Hw3aS0mRtdnY2ZDJZonN0CSNBDCHg9XoxPT0NhULBioeorLm4uBhr1671uagZhvGRJkeC2+0OuTAWY8cQDoQQyOVyZGZmoqioCABgs9mg1+sxOjoKs9mMtLQ0VkPBbQjjvj+j0QiVSsWqNqkXBf1JEMXSQYIY/MBNME5OTqK8vBxyuRzj4+MYHR3F+vXrkZmZGfB3fBczIQQTExMYGhoCAB+Zc2pqqujvRyz4L9jU1FSkpqaiuLgYhBCWKPwbwuj7op8P7RwFPt2RUY+KBFEsHSSIgYNgoYPb7UZbWxskEomPrNkffIjB7Xajq6sLDMNg8+bNAACr1Qq9Xo/e3l44HA5kZGT4SJiXAiK9L4ZhkJaWhrS0NJSUlLANYXq9HgMDA7BarUhPT0dSUlKAKY1/eTRBFEsDCWKYB7Vc41YdaK9DTU1NRFlzJGIwmUxob29HRUUFSkpK2NdKT09Heno6ysrKWJlzf38/RkZGMDo6CqVSiZycHFaUtBjgNoTxAcMwPu+LEAKz2Yzx8XEYDAYcO3YsaENYMKJwuVwBRJGwwYs/lj0xBNMmEEIwPDwMvV6PVatW8ep1oMrHYKBhCO2uDHcMpVLpo1CkoqSRkREACBAlnQ1gGAYZGRnIy8tDSkoKKioqgjaE0RyFXC5n/87f3YoShVarRWpqKkuYCXcrcbGsiSGYNsHpdKK9vR1paWkoLi5mqw6REGzHQEMHAGHDkFCQSqXIyclBTk4OgEBRkkwmY++6GRkZcZNgC90xRDpOsIYw2ucxOTkJt9vtQxT0O+AShdFohFQqTbhbxQnLlhiCaRO0Wi3OnDmD2tpaFBQUoK+vj3d1wJ8YaOhQXl6O0tJSQecWKizxFyVRrcHk5CRMJhNSUlJ8rOLEgljEQAnYHxKJhN0hAQjaEMbdRdFcBTeRSc+TEgUloIQXRXRYdsQQLMHo9XrR398PvV6P5uZmNualv+MD7mLmGzrECq7WgBACu90OnU7HVgYcDgcmJyeXTMWDL8EEawijIRVtCPN6vUhJSUFGRoaPz0XCBk8cLCtiIITAaDQiOTnZR9bc3t6OnJwcbN682efCFSJaouKn9vZ2EEKiCh0ootExMAyD1NRUlJSUsJWBo0ePwuPx+FQ8cnJykJ2dzTtEAsQNJaJZjP4hldvtRnt7O8xmM1paWkI2hNG/pa8NJIiCL5YNMdAy2MmTJ/G5z30ODMNApVKhr68Pa9euZS86LoSIlux2OyYmJlBbW4uSkpJFj23p3bOsrMyn4sGN4xe64iEWwchkMiQnJ6OiogIKhSJsQ5hSqfSpXnCJgoYeo6OjSEtLY8OUhLvVMiAG/9CB3o3PnDkDm83mI2v2B98dw8TEBIaHh5GXlyc4nxAM8VA+BjN24VvxWOwdQzBw8xV8GsL8k7TcxKTNZkNKSkqAuxXXBm+5uVt9pokhmDbB4/Hg6NGjKCoqwpo1a8J+2ZFyDB6PB11dXSCEYM2aNdDpdKKc90JIooVUPMQ6F7EIhh4rFMmEawjr6elBcnKyD1F4vd6AHcJyt8H7TBKD/zg4+oVPTEzAZrPhnHPOCSpr9ke4UIL2TZSVlaGkpAR6vX5J9TcIRbiKh9VqRXt7e8BwHKEQkxhoVYIPuElaYC7soxUPs9nMalikUmlAQxj33JeTu9Vnjhj8LdcYhvHREygUCl6kAIQOJSYmJjAyMuJTdRDaXRluN7IUmqi4i8lkMqGmpsanFyKaHg8hizmex0pJSUFRURGKiopACEFbWxtkMllAQ1h2djY7JnC5EcVnihiCaROMRiM6OjpYKfLHH3/M+3j+OwYaOni93oCqg5BEJSEE09PTIIQgOztbVIOXeCBcLwS34pGdnY2cnJyQOZvF2jGEAy1Zr1ixAqmpqSCEwGq1wmAwYGhoKGRDWCSi+M53voN7770XNTU1MZ/jYuAzQQzBtAlU1jw1NYXGxkYfwQ/fC5R7V6ehQ2lpKUpLSwP+nu9dnpKLx+OBXC7H2NgYALBuSUqlcknsGMLBvxeCVjx0Oh06OztDVjzEzjGIdSyPx8MmWxmGgUKhgEKhCNoQZrPZAoiC/h33fIaGhhatt0UMnL1nPo9QsuaOjg6kpqZi69atPncWuuj4XFT0uZOTkxgeHg4rWOKzmC0WC9ra2lBWVobCwkL2nOkMS7Vajf7+fng8Hnbh0Zh3KYNb8QAQsuJhs9lEFXwtxO4jVEOY/26JVnSoOI52lJ6tOKuJwePxYHBwEPn5+exsR51Oh+7ublbW7A+6C+C7DVWpVLDZbBEFS5FCienpaQwODrJ+DnRwLe0Y5Cb+hoaGYLPZ2JhXoVCwFYRgo+mWGkJVPLRaLfr7+zE+Pr4gPR58IeR6oA1h4SaEvfHGG7DZbLBYLMjLy+NzzBQAHwCQY25NHiSEPMwwTBWAPwLIBXASwA2EECfDMHIAvwXQDEAL4DpCyHAUbz0kzkpi4IYOBoOBlc729/dDp9Nh48aNIRNifGXOZrMZZ86cgVwuR0NDQ8S7U6jko9frRW9vL6xWKzZv3swrn5CUlAS5XM6aoFgsFuh0OvT09MDhcPhs05d6fgL4tOJhMBiQl5eHtLS0sD0eC71DiiVfEawhzGw24+2338b+/fthMpnw4YcfRlKaOgBcRAgxMwyTBOBDhmH+DOA+AD8hhPyRYZifA7gFwM/m/6snhNQwDLMPwOMArovqDYTAWUcM/toEiUQCm82G7u5uZGdnY9OmTWG/ZD7EQEOH6upqGAwGQWEHF3a7nXWRXr16Ne8Lnnss7laW3qH8+wZo0o+q/JYqqPbAv8fD3/2JT8VjoXIwHi/BgMYKvdWJrNQkrMxXQCYJr3256KKLkJKSgvfeew8ulyui/JzMvRnz/D+T5n8IgIsAfHn+8RcAPII5Ytg1//8AcBDAMwzDMETED+WsIQZ/bQLNJzgcDpw5cwbr1q1jZyaEQzhi8Hg86O7uhtvtxpYtW2C1WnmLlvyJgXZqhpJbRwuJROLTYOSfn0hKSmL7IZZaAjNYd6VYFY9YEIqwCSE4PmLAsNaKtGQpBjVWaCxOfK4qOyzJc3NYAtr2pZgLF2oAPAtgAICBEOKef8o4gJL5/y8BMDb/Wm6GYWYxF25oeL0YD5wVxBBMm+D1enHmzBmYzWasWrWKFykAoYkhWNVBaBMVVcsNDg5Cq9X6dGoKgZDSp39+gop3RkdHYbVa0dnZyS6qxc5P8En68q14UPfpeGb+bS4vhrVWFGbMmQFnyGUY19lhKfEgXR7+dYVWTQghHgAbGIbJAvAagDWxnHusWPLEEEybYLFY0N7ejhUrVqC4uFjQ9jkYMdDQwd/oVcgCpXLrU6dOIT09PWJIEy9wxTvHjh1DRUXFkslPRFNiDFXxUKlUITsr4wkCINK9wuPxRP3dE0IMDMP8DcA2AFkMw8jmdw2lACbmnzYBoAzAOMMwMgBKzCUhRcOSJYZglmvAp4t43bp1UCqVGBwcFGTbziUG/9DB/+4jZMdgMplgNBrR0NCAwsJC3ucTDGLpGCLlJwBf/US8iUwM7QHtnJTL5di0aVPMrlbhPufUJAkqclIxrLNBkSyD1elGaVYq0uXhyUdoqZJhmHwArnlSSAVwCeYSin8DcA3mKhM3AfjT/J8cnv/3v+Z//1cx8wvAEiWGYNoEt9uN7u7uANUh9UHgC0oMVFNQUlKCsrKyoBcsnx0DIQRjY2OYmJhAWloab1JYDG2CkPxEPPQT8ejSFOpq5f/6kTQMmyuzkZ2WBJ3VjZy0NKzMj1w1obJqASgC8MJ8nkEC4GVCyJsMw3QB+CPDMI8COA3gV/PP/xWA/2MYph+ADsA+IS/GB0uOGMLJmsvLywO8DqRSqeAdg0ajgUajCTkjgvvccMemPRgSiQSbNm3CiRMneJ9HOCyU8jFcfsJsNrPVASGfbziEsnaL5jihFnOkigfVhNCKR6RSpUzCYM0KYaIsWlnhC0JIG4CmII8PAtgS5HE7gGsFnZRALBlioNqEzs5OlJeXQ6FQgBCCkZERTE5OoqGhIeiHHc6d2R8ejwdarTbijAjusUMtUG53ZWlpKZt4FAOLJYn2by6i+gm73Y5jx47FnJ8Qy4+Br+6AT8UjLS0NbrcbTqdTtIoHJdWzGUuCGLjaBLfbDY/HA6fTic7OTsjlcmzZsiVkUkkqlfp0uIUCDR2ocIhPNjtUKOGvYqTPFQIxiSQe4OYnZmZm0NzcHHN+QqxQItqdR7CKh1qtxvDwsE/Fg4Ye0VY86M7kbMaiEkMwbYJUKoXBYMDY2Bhqamoixux8BEtTU1MYGhrCunXroNPpBJcgKbxer89A22iz+jabDa2trXC5XD4+jNx5CkuNNPzzEzTpR+3xkpOTI+YnlpoTlEQigUKhQEZGBurq6gJ6PKKteCSIIQaEslwzGo3Q6XRobm7m1ecvlUpDJh89Hg/OnDkDl8vFLmSDwRCV87PdbkdraysKCgoiOj+FAxU+rVmzBmlpaex2vaurC263G9nZ2UtavUjhn/QLlZ/g6ifE3DGI9RlxOyvFmuNhsVji6g6+EFgUYghmuWa329He3g6GYVBbW8vb/CPUjiFU1UGIJTyFRqNBT09PTCpGQgiGhoag0WiwadMmyGQyuFwuVmdPfRj1ej0mJibY8ie9UBejh0AIQuUnuPoJp9MJj8cTsyhJrCQmPVaoBR5txSOxYxCIUNoEtVqN3t5erFmzBrOzs1HrEihoDoBqHfyfzzdZST3/hoaGsGnTJnarLxTU7jwlJYUVPgXb5UilUuTl5bGtvUVFRdDpdBgZGYHZbPaRBkd7LguBUPoJtVrNkn8s+gmxTWWjsYgLVfGwWq3slHQ+GBsbQ3l5+d8AFGJOP/ULQsjTDMPkAHgJQCWAYQB7CSF6Zo59ngZwBQArgJsJIacEvu2IWDBiCKZNoJ2HFouFXXhms1mQLoEbSng8HvYOFSoHwHfHQEfVEULQ3Nwc9YVoNpvR1taGyspKXjMwgU9DmJSUFBQXF7NdlmazGVqt1ifsoGPclvIcS5qfkMvlaG5uZrtiheQnuIhXKCEEoSoe77zzDt58803o9XocP34c9957L9atWxfyOPO7p/sJIacYhskAcJJhmHcB3AzgfULIYwzDPADgAQD/BeCLAGrnf7Zirqlqq+A3EAELQgzBtAlWqxVtbW0oLCz06TwMlzMIBrrQ6fGKioqwdu3akBcXH2IwGAzo7OxEbW0tHA5H1Bch3bmIMZGK6wNAww6DwQCdTofBwUHIZDLk5ubGTZwkFhiGCZqf0Ol0YfMTXCwV70gu6C7p6quvRmdnJ7Zt24aioqKI/qLzodcpACCEmBiG6cZck9QuAF+Yf9oLAP6OOWLYBeC380rHTxiGyWIYpogQMhXzm+AgrsQQzHIN8K0S+G/16aBSvpBKpbBarWhpaQl6PH+EIwaqYpycnERTUxPS0tLQ39/P+1wovF4v+vr6YDabQ+5cwi1cPlUJqVSK3NxctnnM4XAELC6XywWHwxFT2LEQ1RH/nVEk/wmxdwxiJ3utViuysrKwbds2QX/HMEwl5oRORwEUchb7NOZCDYDTWTkP2nV5dhADjc9Pnz6NpqYmVtZ85syZkL0JgLAdA12AdrsdF1xwAa/yYShioOIqmUyGzZs3R701dzqdaG1tRXZ2NjZu3Lhgd265XO6T/DObzWhtbWXDDhrTCw07xPRW5AM+/R10orUYBOH1ekXv0Iwm+cgwTDqAVwHcSwgxcj9zQghhGGZB69dxIQY6Dg6Y2yYyDOMz/TncCDe+EmcaOqxYsQJGo5G3piAYMdA8AHWSjhYejwfHjx/HqlWr2G1yNIhVx0DDDrlcjqampqBhB612LOWwAwiunxgaGoLJZMKJEycE5yf8Iebug0Ko8pGZc216FcDvCSGH5h+eoSECwzBFAFTzj9POSgpu16VoEJUYQmkTRkdHMTExEVLWzAWfpqjp6WkMDAywysPJyUne5+hPDDSsiTUPMDY2BrvdjnPPPVdoA03cwSfsoEThH3Ys9I4hEpKSkpCeno60tDSUlpYKzk/4Ix6hhBAdw/wN4FcAugkhT3J+RTsoH0NgZ+UBhmH+iLmk46zY+QVARGIIVnVwuVyw2WwwmUxhZc1chAslqDmLw+HAli1bolIeUmLwP1ao7SQlt1CLg7Zue71epKWliTJuPt7Kx2BhB1dkxQ07liK4d3mh+YlgxxK7oiOk7fqjjz4CgBsAtDMM0zL/8LcxRwgvMwxzC4ARAHvnf3cEc6XKfsyVK/9dtBPnQDRi4AqIAECv16OrqwtJSUmoq6vjfdcJRQzc0CFc1SESqI7h+PHjKCwsjHgsSiTBLh4qbS4uLkZZWRk++eQTUe6wCx3T02pHRUVFQNghlUrhcDhgMpmWTNgR6vuIxn8iHqGEkO7K8847D4SQUB/qxf4PzFcj7orh9HhB1FCChgGDg4PQaDRoampCW1uboMUSjBho6LBu3bqY72IGgwF6vR7Nzc1s3BoOVG/hfyFSNWRdXR17HDEvtMXqlfAPO6hbFp+wY6HA9zPm09/hcrmgVCpFDZncbvdZ4d4dDqISg8PhQEtLC5RKJTZv3gyJRAKpVAq32827pZVLDLRpic51iOXDJoSw9vKZmZm8SAEIbKTylzZzFwffEICWRTUaDXJzcwPi4aXURJWcnAy5XI5169ZFDDsWSmQVrfIxmH6io6MD09PTGB4eFpyfCIWlsKuKFaISw+DgIKqqqnyGbMhksqiUjNzQIZamJWCuhNjW1galUommpia0tLTw/ltu6zWVNqempgb1dOSzoL1eLzv/sry8HLOzs2yjF11kS6mJinsnjRR2LFS1Q6xdWUpKCuRyOWpqapCSkiLK/A7aTn+2k4OoxEBbV7kQqmRkmLkRc6dPnxYUOoT6MqiKkZYQPR5PVL0YtKRZVVWFoqKikM8NRwy0Q3PFihUoLS1l+//Ly8t9FplWq4XH40Fqaipyc3PZictLDUKqHWLavovZREUl0dHkJz7LiLskmoYSfEBDB5fLhW3btgkOP7iVBVomnZqa8plMJbS7kmEYqFQqTE5ORixphtsx0GQs7dD0PwfuIissLMTo6CikUikGBwdhtVqRmZnJSp4XMn4VcvcLVu2g79vlcsHpdEKr1cYcdixEE1W0/hPxtrRfKCwIMfDZMXBDh9TUVEF3GP/F7na70dHRgeTkZGzZsiVgqC1f0JkGTqeTlzFLKMen0dFRTE5Ohh2d538cqVSKkpISlJSUsOeh1WrZ6dj0TpyZmRn33US0bkncGY9OpxOnTp2CXq+POexYjF4JPv0dSqUSGo1GkOpx//79+PWvf60CoCKErAeAxe6sBEQmhmBfLh9imJmZQX9/Pxs6CBEs+b8GVVgK6WYMBofDgba2NkilUqxevZrXXdp/x8DNJ8Qis/afreByuaDT6TA5OYkzZ86wBqe5ubmiVwrESoJKJBIkJyejpqYGQPCwg2/ibyk0UQXTTwwNDeH+++/H4OAg7rrrLlx11VW4+OKAiqMPbr75Zvz617++HHNDaikewCJ2VgILsGOQyWQhQwlu1WHz5s1Rx6F0x8BnXD0fcPMSKpWK9+II5vhUWFiIiooKQVZnkZKYSUlJKCwsRGFhoY+oh1spyM3NFWVRi2nHxj2Of9hB34N/IjZY2CFmjkGM49D8RH19PZ577jk89thjuO6662A0GiP+7QUXXADMWcBzsaidlcAihhI2mw1tbW1BrdLollzIaPK+vj4A4OX+HArcGRG0u1Kj0Qi2gqNx9Zo1a8KOzhPzoqRJM5rE1Gg0bNcpLYkuZhIzHMEEew+zs7PQ6XQYGhqCVCpFdnY2cnNzkZ6eLmqOQWxQOfT8go8Wi9pZCSwQMfg7JtEEDlccxIUQoZDNZoNWq0VxcbGgidL+8Hg86OrqAgAf+bbQ+ZVTU1PQarUR8wmxtl2HAjeJOTs7i9WrV7PlRJvNhszMTDZpxic8WgwDV3/vRWqpNjY2BpPJBJfLBbVajYKCgkWfx+kPsa3jF6OzEliAHINMJmPt3bmOTeFCh2BVhmCglnA5OTkoKCiI+gKmiU8qbQ62e4kEr9cLvV6P1NRU3n0hC4HU1FSfJCY12x0dHQXDMGxuIiMjI667iVi2//6WaidPnvQx+qV6g+zs7EX/3EXye1zUzkpgAcuV3NAh0p09UsKSEIK+vj4YjUZs3rwZw8PDgrQS9BgMw7DS5lCaCT47Bqr4TEpKQnV1dcwXZ7yUjxKJBFlZWcjKykJ1dTWcTid0Oh3Gx8dhMpmgUCjYsIMmMeOVY4gWtEGvvLycvU5o2DE8PAyJRMLuNiKRXTw+Y6FTqEJgUTsrgQUiBpPJhFOnToUMHYL9TaiFTqsFWVlZaG5uZp2hhGoTvF4vhoeHodPpwhq9Rtox0ETlmjVroFKpQj5PCBYqD5CcnOxzJ7ZYLAGekmJ0iwLitm9zj+UfdgQjO/p7/7AjHg1UQgfaXn/99cDccNo8hmHGATyMRe6sBOIcSni9XoyNjcFsNuNzn/tcVP0SXNCknr8RSjTqypaWFigUiohGr+FIZ3x8HGNjY1ElKiNhoXsluAlAKnfW6/WYmZmBVqtFa2sru8CiSWLGixj8EYzsuNUOrsw5HklMs9kcNuHsjz/84Q/4wx/+EExKK3pnJcMwdwC4Y/6fSgDDhJALgz03bjsGGjpkZWVBqVQKKkX6L3QyP8NyZmYmaFJPyI7BZDLBZDJh7dq1KC0tjfj8YNt66uXgcrl88glCQoDZ2VmMj48HHYe2FOTP1MpeLpdDIpGgsrISOp0OAwMDsNvtPkpMPlUgsfsH+BwrXLVjeHgYwJwYzmg0ipZjsVgsqKioiPk48QAh5OcAfj7vGPVXAE+Gem5ciIFbdUhLS0N7e7ugv+fau9HGJblcznZs+oMvMVC3pszMTJ9Gr3DwP7bD4UBrayvy8/MDvBz4EgPVW1RUVMBoNGJkZAQSiYStJkil0iXTXUnhn8SkC4yeOy0nhlpgS6HE6B92GAwG9Pf38wo7+EKkHEO88TSAvxJC3gj1BNGt3Xp6elh35OTkZHZIrRBQXweqYgzXuASANRMJBe7MyS1btqC9vZ33DoNLDDSfsHr16qDEEokYCCE+VRmv1+sz5YjeyehsDZVKhZycnEXV3ge70/v3EfjH9cE8G8QUJYlFmjKZDGlpaairq4sYdvD9DoTmGBYaDMPcDKACwIFwzxP9iktPT8eqVat8kkNCiUEqlUKj0WBwcJC3T2SohU7v8Hl5eayQSkjoQZOP/vmEUOcR6qJ1uVxoa2tDZmYmmpqaAMDHJp+rBKQ+AWazGaOjo2ymnQp8FjrUiPR6/nG9v2dDdna2qN2VYoE7bCZS2MG32rGUx9MxDNMM4D8BnE8ICbsARHdwKi0t9VkcQi9ij8eDmZkZeDwe3irGUM7Soe7wQolhenoaKSkp2Lx5c9jzCVXBsFgsaG1tRXV1NVasWAEg/F1PIpFAJpOhurrap6xI+woyMjLYsmK8Oy2F3p39PRvo1KnJyUnMzs7CZrOx556amhp1g5YYCFeViLbaIbbASWQcAJAD4G/zn+EJQsitwZ64pPpDqdAoLS0N6enpvLdv/s7S3MExsSQrHQ4HBgcHkZKSgsbGxogXZLBQguok6uvrI04lCgX/OzLttBwfHweAuIqUYk0aymQy5OXlgWEYKBQKdh5nf38/7Ha7jzhpoUMmIeXKYNUOvV7PmrpkZWVBo9HAZrOJMumaYZjLMZcLkAJ4nhDyWKzHJITwLm8uGWKgCct169ax5qN8wV3oHo8HnZ2dkEgkITsa+RDD7OwsOjo62AuBbxacEgOtpKhUKsEDccPlKhiGYSdkV1VVweVysSRB43t6R15KoJ8hnfdYWloaNIkZabsuZlI2lrmVNOwoKyuD1+uFwWDASy+9hNbWVuzbtw9f/OIXcf/990eVxGQYRgrgWQCXYK4f4jjDMIcJIV2CDxYlRCeGUBd1qMXlP84tOTkZarU6aju41tZWlJaWoqysLOTzIxHDxMQERkdH0dTUBKvVCq1Wy+s8aCjh9XrR2dkJhmGCWsCJiaSkpIDdhE6nQ0dHBywWCwYHB5Gbmxu1b0M8lY+hkphU+0KTmLm5uWyOQsyyp1gCJ0poDz30EN577z28/PLL+PDDD2Npgd8CoJ8QMggA80rHXQDOXmIIhlAW7FTFmJOT4zPOje80Ku7xbTYbbzu4UMRAqxd0WrZMJoPNZhPURGW323H8+HEUFRUF9F3wRbSSaO5uorKyEseOHUN6erqPbwMtiQpJBoqxEPlUJYIlMbVaLTo6OuD1epGdnQ2lUimqrZvYpO31elFUVIS9e/dGfnJoBOuijIvvQigsCDHQOzqXGHQ6Hbq7u4OW/oRUMgghGB8fh9lsxnnnncfrgg9GDE6nEy0tLT7Vi1DPDQWbzYbx8XE0NDQIUr/FCwzDoKCgAAUFBSEXGt1NhFogYm3dheoYuEnMyspKuN1u6PV6qFQq1n+Tq8SMBmLPrfwsmMBSxCWUCHgRjlM0IQTDw8NQqVRobm4OGoPxJQaXy4X29nakpKQgIyOD913Qf7HTfEKwmZN8264nJycxOTmJ4uLimEkhHk1UwRaaTqfD9PQ0ent7kZaWFtQFaqk0UclkMuTn5yMzMxMulws1NTXQ6XTo6+vzcXQWksSMR6+ESJ/XgnVRhsKC7RjooNuOjg629BeuVBSJGKj4qbq6Gnl5eTh1ir/1HZcYJicnMTIyElKfEKmJioqWrFYrVq5c6aNNiBYLcdeRyWQ+uwmuC5TH42Ft1sTcMYgVkkgkkqBJTK1W65PEjKT7iDb5uAA4DqCWYZgqzBHCPgBfXsgTWDBiMBqN7EKmtfxwzw9HDFRSTMVPXq9XULKSlje7u7tht9vD6hPChRJc0dKGDRugUqnCKjCFYCEl0f7iHu62XavVQiKRICUlBbm5uVFLhcUmBi78k5jB/CRppYa7qxR7x+B0OkURchFC3AzDHADwDubKlf8fIaQz5gMLwIKEElarFTqdDk1NTbzEH6GSj/7SZrqYhW69vV4vRkZGUFJSEnGYTahjBxMt8TV1oX/f398PpVIZsOAWO06l2/b8/Hyo1Wp2zgWVCtPdRFZWFu+FJdbwWD6L2d9PklupobmVnJwc0ZOPYqoeCSFHMNdmvSiI646B3pWdTidqamp4K8L8BUtAcGkzhZCFZDQaMTQ0hNzcXKxcuZLXufgv9lCiJb4ERf++qqoKNpvNZ8HRZOBSaqKSy+UoLy9npcJ6vR4ajQb9/f3sToIqGUNBrCYqoTsP/0oN3Q3NzMxApVLBYrGgoKAAubm5MXtPmM3mJSuHFoq4EQPVFJSUlCAtLU3Qhe7/xfM1V40Emk+gMmM+4CYfI4mW+BDDyMgIpqenWZMZhmFQVlbGSoep0Mtut2NqampRh8dScL8P2o5NK0lU59Hb2wun0xnS3TmeoYQQcHdDbrcbK1asgM1mQ29vL6tgFNo4RXGWdFbyQlxCCTonYv369VAqlRgdHRXcSAWEniYlFNRrktrU6/V61ocyEvxFSxKJJKRoKVwoQT0c3G43Nm/eDIZhYLPZAHzaG0EXHCEER48ehdPpRGdnJ7xeL5tMW4gBM1xEIjqaBCwrK2MdqrVaLQYHB5GcnMzuJpYKMXBBCIFCoUBeXp6PgpG6UwsdimO1WhM7hlCg8lyu2atMJhOcrSeEoL29Pay0mQ+cTifrPkS9JoVoEyQSCdxuNytaKi8vD/vcYAvJ5XKxGonKykoAn4praOKUEqdEIgHDzE2iqqioYBuRuANmaDItNzd3QZqo+C5o/1mW1MG7v7+fNUORSqUxjagTkxj8qxJcSTbwaRJzZGSEtYUPN4szEUqEgb+KERDeem21WmG1WlFRUcHLZSkUaCWktrYWBQUF7OPROD41NzdHDGOChRI0SVlTU4OCggJ2qC7DMOyi9nq9IISwv+P+l2GYgNKi2WyGRqNBW1sbALCLMV4t2dEeMzU1FaWlpSgtLcXAwAAkEglrZZ+UlBTVvAsxfR0ikUywJKZWq0V7ezsIIQECsSXeWSkIcQkl/L84IYNtqSV8SkoKSkpKBL0u94umbk0bNmwIYHG+xEBzEmlpabxyG/7EQJOMtKzKJQXuZ0TPWSqVsl4M5eXlIIT4lGLp3Y0KlaqqqgJasqnlmlhNVGLqGOi5AXO7Ca5VHK3ORLKAF9MJSuhQI//mNb1ej+npafT09GBiYgL/+Mc/eDuD+R37WgCPAFgLYAsh5ATnd98CcAsAD4CvEULemX9c9O5LLhakXMlnx0AIwcDAAPR6PTZv3ozTp08LKnHRxc4wDCs4CuXnEIkYuKKlzZs349ixY7zPgS4kmhvZtGkTkpOT2TxFMOKkMBqNrOM0rcnTv+MSBA1DGIYJ6C8wGo3QarUYHR2F1WrFyMgIcnNzoVAoFr2JirsIg1nFabVaNranuyD/3YTYOYZoj5WUlOSzixscHMQ777yDN998E++//z4ef/xxbN++ne/hOgBcBeB/uQ8yDFOHOXHTOgDFAN5jGGbV/K/j2n25IAInriQ6GKhQKD09HZs2bWLzAEKUadTerbu7G1lZWdiwYUPIC1qIaEloaYxOtHK73WyS0uPxsAss1PFUKhUGBwfR2Njoo8CUSCTsxUtDDn+yoJ8Xd/htdXU1jh49iqSkJAwNDcFms/G+K8cD4QgmmECJJjDp9Cx63vGQMccKhmGwcuVK1NfXY/v27fjyl78cMH0tHAgh3fQ4ftgF4I+EEAeAIYZh+jHXeQnEuftyQSXRwUClzStXrkRhYaHP3wjJS3i9Xpw+fRqrVq3yyScEQyhiCCZaEgKaJKysrERVVRUARCQFWgLV6XRobm4Om0zkhhwAAvIR/glMiUTCTmT2vyvTGJ/elUNhMXol5HK5z3nTXdDIyAhcLhc7cj7aXVC8QMuVKSkpYo3OKwHwCeffdFYlEOfuy0UNJfylzXz+JhimpqZgNBrR2NgY0AQVDEJES3xhsVjQ0dGB1NRUVFdX+2z9Q93haAkTADZs2CD4ThhsN0GJghACt9sNt9vNPo97V7bb7WzFwG63sxOyY6kYhEO0BMOdngUAg4ODcDqdGB4ehsVi8dlNxLtCEwnhjGC3b9+O6enpgMe///3vY9euXfE+NcGIy47BPwnnH0rQBeF0OkPmAfjmJajrcn5+Pm8hEJcYaLenWq0W7LREodVqcebMGaxevZrVbETKJ9CQJS8vD+Xl5THf+bi7Cbfbjba2NhQVFbHhDd25SKVSMAzDJndLSkpY41OtVouBgQEfNeNS6a6koCVFGttzcyqLbZobTuD03nvvRXPIcF2Wce2+XDCjFrrI7XY7WltbUVBQEDCXgYtIZi1Op5MdaNPU1ISuri7BlvDUBk4qlUbttMRNMtI7Njc5GAzU27K6ujpi2CMUdrsdbW1tKC8vZ8Mh/90EAJ+dhL/xKVUz9vT0sPr/tLQ0KJXKqON7scqM3BwDwzA+ORUhprnxkJzHwSH6MIAXGYZ5EnPJx1oAxwAwiHP35YIQA70ghEibqbAoGGhegmoDAGGhBz32iRMnIoqWQoGrZKSk4na7IZfLcezYMWRmZiI/Pz9gLoRer8eZM2ewfv16UUxDuTCbzejo6MDq1at9ZoQKzU1w1Yyjo6Ow2+2sVDs1NZXNTQjZXYlVZgxHMEJMc+ORxIxWEs0wzB4APwWQD+AthmFaCCGXEUI6GYZ5GXNJRTeAuwghnvm/iWv35YIQAyEEDocDvb29AeYsTrcXozorvAQoz0lFStLcxRtqxzA9PR103oQQ0ZLRaITVauUlWgoGl8vFqim5SUYAWLduHYA58xe1Ws0m+vLy8uD1eqFSqdDU1CRWcoqFTqdDb28v6uvrI961wuUm6HuhJEENXqjIh+4muL4NeXl5EaXa8Sp7hkIw3YG//bvL5RKtVRoAq44UCkLIawBeC/G77wP4fpDH49p9Gfccg9vtZvX+/uYsNqcHP//nEEa0VjAMg7z0ZBz4QjWUqUlB51f29fXBZDJh8+bNAVtDvj6RQkVL3NdnGIatXNAqCk3yAb5JRm7CzGazobu7GyaTCXK5HKOjo8jPz49pa+7/niYmJtDU1CQ4R8LdTXCJgZZDXS6Xj3hMoVBAoVCwvg3BpNrBJMOL3SuRlJSEwsJC9jujideOjo4AFWO055loouIJi8WCtrY2lJWVwWKxBHzgHw1oMaK1oiRrrjlq2mjH250zuG5TqQ8xcLUF/nJrimCt2lxEK1oCPiU6Or6svr4eGRkZvJKMHo8HfX19yMjIQFNTE7xeL3Q6HaamptjFRJunhGbVCSEYGhqC0WjExo0bY64m0AVH/0snXa9fvz5oyBFMqk0lw4Dv1n2xiYELGi6lp6dj/fr1rIoxVtNcq9Uatf/kUkPciIHGpLTDcnx8PEDJqDY7IZd9+u+0ZCk05rlmK6lUCpfLFVLn4I94iZboscfGxjA9PS1IyUhdsIuLi1l5t1QqZdt+6WJSq9VoaWkBAOTl5SE/Pz9ijd7r9aK7uxtSqZTXMByhGB8fx8zMDKut4IqruApMrriK6ylJ511QK3j67/z8/JjKivFIYvqrGC0Wi2DTXABs1eezgLgQAzV75XZYBnOKrslX4F8DWri9SZAwwKzNjQtq09nnazQaTE5O8ppfGSr5GKtoyev1wmq1Qq/Xs0lGPqRgMpnQ2dmJVatWhexboPF7RkYGm1XXaDQYGBiA1Wpl43d/pSItR+bm5opS6uSCynvNZjM2bNjAvi435KBEQUugwXYT/vMujh8/znp0SCSfTvYWKlISK4kZSlXLcGzuaGcr7YlYCNPcpYK4EENhYWHATAWZTAa32+2zNdtUkQWVyYH3z6jhBXDuyhx8YdWcH8H09DSMRiO2bdvG6w4jkUgCZKixipZoklEqlWLNmjW85c1qtRoDAwO8EoFcJCcn+yj+DAYD1Go165SUl5eHjIwM9PT0oKKiIuwOKhrQSotEIkFDQ0PYC51PApO7m5BKpaiqqoJUKoXT6YRWq8Xw8DCsVquPSCmSOYpY1QS+x+Eau/ib5tKBvVQ3EQ05fOMb38ATTzxxBoATwACAfyeEGIDFa6AC4kQMaWlpAXfvYHd0hmHwpfoVuKyuAF4CJMvmFvepljb2C+G77eQmz7iiJe6uxR/hvkhuknFiYiJkktH/eGNjY1Cr1di4cWNM2W5/bwCLxYKJiQm289RsNiMlJUU04xaPx4P29nYolUpUVlYKOmaocih3N8HNxyQnJ7PtzFzJ8/DwMGtYQ63W/M9joYmBC+5ugmuaOz09jZtuuglmsxm//OUvcfnll4edhMbFJZdcgieeeGL9vAHs4wC+BeC/FrOBCohjVcIf4XQGMuncF0QHiVRXVyMlJQUTE/zFXPRuzle0RBOKwc6VKhlpknFiYgJGoxE5OTlh5c09PT3wer1oamoSvUZut9uh0+mwdetWJCcns6PcTCYTMjMz2cUUzQAVujMqKioS1OoeCv67ieHhYcjlcpYouOIqruR55cqVEaXaYhGDGNbx3N3E3//+d1x++eWwWCz4xS9+ge9973u8jnHppZeCEEIFO58AuGb+/xetgQpYwKG2kTosZ2Zm2O13RkYGTCaTYEt4p9PJW7REcwX+Fxmdkt3c3Ay5XA6v14vS0lKMj4+jr68P2dnZyM/PR3Z2Nvu3dPBNdna24LstH0xMTLCTu+kuhFt6m52dhUajYe+4+fn5yMvL45Uhp0rU6upqXn0mQkBb6R0OB9sLEklcxZVq03CKSrXlcjnsdjscDkfME6TEFjjRnNC9994by2H2A3hp/v8XrYEKWEBiCLVjIISw1l9cfYLQ7kqbzYaZmRk0NTXxMinxr2LQO77T6QxIMtJkk9frhV6vZ81kFAoFlEolJicnUV1dLXrMz00EhipHMgzD3nFrampgt9vZ3IrD4UBOTg7y8vKCWr2bzWa0t7dj7dq1Eed9RnPu3d3dkEgkWLduHUuW/rsJ7g/dwXG7Q/2l2qdPn0ZfXx9cLpfPbkLoIhfbOj6cexOfBiqGYR7EnLrx96KdVAxY0FDCX+JM77Tp6ekx2cHRLk3uRRQJXGKg5cysrCysWbMGQPB2aW42nRCCqakp9PX1ITk5GWNjY7DZbGypMVZ4vV50dXUhKSkpYiKQi5SUFNZOzePxQKfTYWZmBj09PVAoFMjPz0dubi4sFgvrLiW2T6HX60VHRwcUCgWqq6vD+jCE85pwu91g5v0vJZK56VPJycnYsGEDa2NPk7NCpdper7hzK2NpoGIY5mYAOwBcTD5t4li0BipgEXcM3HxCsDKikO5Kq9WKhoYGDA4O8j4fSgxWqxUtLS3sefBJMgJzoc/4+Di2bt2KlJQUOBwOaDQa1vqdhhzR3M0oUeXn50fVx0ERTDOh0Whw4sQJ2O121hlZzFIbLaXSrlG+CJbA5IYb3NmnVA/DddUOJtUOpz2IRygRDcG+/fbbAPBNAJ8nhFg5v1q0BipggXMM1ClapVKhv7+fzScEQ6TeB3/Rkt1u590rQY+v1+sxPDzMirD4KBnp9p6qDeldRy6X+7Qxc+/UQtSNNpsNbW1tqKqqErXzkmomDAYD5HI5GhsbMTs7i6GhIVgsFmRlZbG5k2iTctQNu7S0FEVFRTGdL91NyGQylrxmZmbYG4Z/OdRfqq3X6zE1NcXulPyVjGLPrYzWIfrAgQMAkAHg3flr7hNCyB2L2UAFLHAo4XK50N/fD4PBELTfgQuJJPSU6WCiJb69EhQ2mw3Dw8NskpGvvLmrq4vdzoZrGefeqU0mE9RqNU6fPs3e6fLz8wOSg9Tzsa6uDkqlkvd74QNKaBaLha2aKBQKH80EnS4ll8vZBCbfZq94JjElEgmmpqYwOTmJpqYmlhy45VCu14S/9oBKtamSMTc3l7W6EwvROkT39/cDvqEBi8VqoAIWcMcAzDktFRYWslOYogFN/PmLliL1SlAQQliTmPXr17OVh0geCtT/gYq3+ILb5UfLcf7JQToVaWBgIMDzUQxwhUv19fUB7zGYZkKj0aCzsxMej4c9x1CaCeov4d/uLRbGx8fZrlR/JSYfcZW/VJvu5mZnZ9kEZjDPBiEI5950NmJBiMFsNqO3txepqalYvXp1VMeIJFri03btdrvR2toKpVLJxqb0bhOOFKjPQW1tbUwj8oDgycH+/n6YTCbk5ubCaDQiKSlJNJuyaIRLdFtOJcHU08BoNAZoJkwmEzo6OrBu3bqo1KWRMDo6Cq1Wi8bGxqBbf6FeE7TLUqfTsZoN+v4YhmGrOEKl2p+lmRLAAhADzSfU1tZCpVJFdQw+oqVwoQfw6SzNqqoqrFixgp2aTe+CoS4CrVbLNoOJ/cVLJBLMzs5CLpejqamJvVOPjIz4bIejHc1HY/6SkhIUFxdHdQyZTOajmTAajVCr1RgZGQEhBHa7PW6kMDQ0BJPJhMbGRt6JwkhSbbfb7dO3o1AoWM8GrlTbYrH4uGpHqmBYLBZRxGFLBXH1YxgYGIDBYMCmTZvgdrsxNTUl+Fg0do3WaQmYMzHp7u72STLm5uZicHAQo6OjyM3NRX5+PjIyMnwIgnZUxipvDgavd24WplwuZ7f31KaMhhxqtZqdFk7PUalU8rqT0SSmmDE/9xw1Gg16e3tRXl6O8fFx9Pf3x1SJ4YJeO3a7HevXr4/6WKF2Ey6XCxaLBcAcedLfh5Jqj4yM+IzfCzY567M0txKIEzF4vV60tLQgLS0NGzduZL8goYNtPR4PTp48ibVr10Y9WWl8fBzj4+MBScasrCw0Nzez7cAjIyMwm81sR6NGo4HL5RLF58AftKJSUFAQMl+RkpKCsrIydhK2TqfDxMQEuru7I0qgqXApHklMYM5Fa3R0lG1BB8DqCriaCVqJEUKqtATt9Xp9hFFigIabXV1dqKysRFpamk9fh7+4imu24z/rwn9GR7TuTQ899BAeffTRNgBeACoANxNCJpm5N/40gCsAWOcfPwUADMPcBOA784d4lBDyQqyfjT+YCKaYUTtmqtVqn+2lx+PBiRMnsHUrP/Xm5OQkOjs7sXXrVt7b1I8//hif+9znAMxdYDRcqK+v59Uu7fV62T4J2odPs/Nixfyx3sm5EmitVoukpCT2HFNTU6HX69mO0njcwahPQ2NjY8jtNVczodFoAIAliXDuzVQtKZVKsWrVKtGl5W63G6dPn0ZZWZmPdsZfXEXXBDeByQVXqm0wGPCPf/wDJ06cwI033ojrrrtO0DnN522Y+df7GoA6QsgdDMNcAeBuzBHDVgBPE0K2MgyTA+AEgE2YW58nATQTQvTRfSrBEbccg1Kp9In5hVQNqGgpOzs7qi08N8nY2NgIIPLgF2Cu8jA0NITa2loUFhayJiq0zBhrzD87O4uurq6Y7uT+EmibzcaGHFarFV6vF2vXrhW9skGTv0aj0cenIdQ5+s/XpINuqGYiLy8POTk5Ptv8rq4upKSkYOXKlaKTAs23lJeXB0jXhSYw/as4WVlZ+PDDD/HMM8/g8ccfx1tvvcVbx+F301Pg05vxLgC/nVdCfsIwTBbDMEUAvgDgXUKIDgAYhnkXwOUA/iDwIwmLuBED1/eR/jsS6BZbqVRiw4YNaG1tFRx+UDs5mmTkVh7CkQJdtNy+Aa6JCjfmd7lcbMzPt+2ZejRs2LAhamIJhtTUVNasZWZmBsXFxZiZmUF/fz/rVJ2bmxtTOET9Nl0uF7v7EgL/2J1qJmhjVG5uLjQaDbKzs1lzXTFBSaGiooKXaEyo10RtbS3S09Px4x//GCUlJYKNfhmG+T6AGwHMArhw/uESBDZLlYR5XFQsqI4hHEKJloQQg8fjQUtLS1AlYyR58/DwcNhF6x/zU9syk8kEpVKJgoKCkKrBsbExH5s0MUETdVarlc3nFBcXsyGHWq3G4OAg5HI5K6wScuFSCzmZTIa6urqY7+T+d1tq3UcIgcvlgsfjQV5eHu8kayQIJYVg5wuE9pqg8vnZ2Vmkp6cHvX4iNVERQh4E8OC8McsBAA8LPlGRsSSIIZRoSQgxjI+Pw+Fw4JxzzoFCoeAtbx4eHoZer8fGjRt5L1puCc/faSk1NdUnL9HX1weHw+GThBULXN9Hf+ESN+Sora2F1Wr1ES2FqsRw4fF40NHRgczMzLi0k7vdbvT29qKiogIlJSUs4dIka0ZGBjubIxpCdblcOH36NKqqqkSrzATbTXR3d6OzszPk5yNgCtXvMadofBihm6gmMBdOcB//O/93wA8LTgzchp1IoiU+MmduklGpVLLa+kikQBeVRCKJam4kBfcOSK2/qLkrVcOtXbtWdFLweDxsRyifRZuWloby8nKUl5cHVGKUSiW7AOmdkeZpwlVOYkGwvopgmgmq6+A2TfFJqsaDFPwhkUjQ29uLW2+9FX/961+jarvv6+tDbW0t/ecuAGfm//8wgAPzRixbAcwSQqYYhnkHwA8YhqES00sx5/okKuKaY/AHrQzQnQAf0VK4HYN/kvH06dNsXTpSkrG9vR35+fkB3pSxgFp/JScnQ6vVorKyEjKZDL29vXA4HIK1CKHgdDrR2toatXCJa9Tq9XrZkIPOrczOzsb09DQqKiqiMtDlc/4tLS2orKwMub0Ppuvg271Kj19dXY28vDzRz59iaGgIN910E1544QXU19dHdYwHHngAhw4d6sBcuXIEwB3zvzqCuYpEP+bKlf8OAIQQHcMw3wNwfP55/48mIsVE3MqVbrc7YFGfOHECDQ0N8Hq9vERLg4ODSE1NDZrhpUrGyspKdkpST08PDAYDCgsLUVBQEDSWtlgsrB19PO4ktG/A//gejwdarRZqtRpGozHoXZoPaLlz5cqVcbno9Xo92tvbkZSU5NPwJdaQWIfDgZaWlpjOn2omNBoN9Ho90tLS2PANwIKQwujoKK677jo8//zz2Lx5c6yHW3L20gtKDKdPn0ZxcTH6+/t5iZboFrK0tNTncToDkyYZqYcCwzBwOBxQqVRQq9VsIqugoAAKhQJ6vR69vb1Yt26d6HMjgU8rG5EkwjQxqFKpoNPpkJKSwpZCw5VnaV9CvIRLlDTXrFmDrKwsuFwuaDQaqNVqUVqz7XY7WlpawlrqCwU3fFOr1TCbzSgoKEBFRUXcJl5PTEzg2muvxc9+9jNs27ZNjEMuH2LgZmwpPvnkE7jdbmzcuJFXnX18fBwejwcVFRU+j42Pj2PDhg1sZ2SofAK9sFUqFYxGIwghWLNmDfLz80W/YFQqFYaGhtDQ0CC4HMm9sAGwJMGNpYXMpowGlHRCDdvlJln1ej1SU1PZ3QQfrQndSVHSERt0J1JZWQmv1wuNRhMyfxILpqencc011+Cpp57CBRdcIMKZA1iuxEBFS9PT01i/fj3vDsWpqSnY7XZUVVUFuDXxHfxCPSWtVitWrFgBrVaL2dlZZGZmoqCgQJQLZnR0FGq1Gg0NDTGXI51OJ0sSdrudnZatVqtZMhQbVC3Z0NDAi7DpXZruJoDwykaqLYlXsxUlhdraWp+dCM2fUJUoLdlSlahQqFQqXH311fjRj36Eiy++WMy3sPyIgStacrvdrOCGD1QqFWZnZ1FVVYW2tjZkZGSgpqYGAHiRAi23KRQKHzWd/1Y+NTUVBQUFgqXPlKxcLhfq6uriUnno7e2FWq2GTCZj736xCpa4oDqHxsbGqCdw0wlaGo0GFovFx0nbYrGgo6MD9fX1cWlLDkUKwUBLtrQPhvpM8EkGazQaXHXVVXj00Udx+eWXi/kWgOVEDHTryRUt9ff3IyMjg3dZR6PRYHp6GiaTySfJyEfJaLfb0dbWhtLS0rCZe3r3U6lU0Gg0rPQ5VPKSIhTpiAWucGn9+vVgGIatHtC7Hw05ot1FTE1NsWGZWMIrrpO2RqOB0+lEVVUVSkpKRO9QpTmLaAxiaGMaTQZnZGSwjWn+n4Ver8dVV12Fhx56CDt27BDzLVAsH2LQ6/U4ffq0j2hpaGgIcrmcd4ltdHQUfX19aG5uRlZWFm+jVhovr1mzRvAFQ6XPKpUqIHlJFz8tF3KH1YoJrtowVDOR1WplQw5CCO9huBR0YlZDQ4OobskUer0eZ86cQW1tLUwmE9tMRasHQo1Q/EFJQYycBdd+T6vVQiKRsF2hqampuPbaa/GNb3wDe/bsiel1wmD5EIPL5YLdbve5S4yNzUm8+QhmJiYmMDw8DIVCgcbGRl6hA/BpTwLfeDnSe6BxtNVqRU5ODjIzM9lGq3iUw6hwKTs7GxUVFbwWD93Kc88zVI2fEMIaoETT98AHdJKUf3jCPU+bzRa1f4OYpBAM1PH74Ycfxvvvv4/m5mbcd999uPDCC+M1zXr5EAMhhHWFppicnITD4QjbKMNNMq5cuRKDg4Oor6/nlWQcHR2FRqMRJQnoD4/Hg7GxMQwNDSEpKQnZ2dmiJS8pYhUu0fPU6/VsfobKimleore3F263W5S+h2BQqVRs30m40IE7F8JgMLAzLyLleWw2G1pbW+NGChQWiwV79+7FjTfeiJKSErz33nv44Q9/mCCGeYhKDDMzMzCZTGwC0R90JgFNMlL3pg0bNkAmk4WVN585M6ckpVOpxQZttGpsbIRcLg9IXtJ4P1pCohd8TU2NaDsRrhWbVquFw+GAQqFAXV2dqB2eFNPT0xgbGxOcs6D+DTQvQbfy/iVb+hmtXbs2LjoO7utcd911+PKXv4z9+/fH7XU4WN7EoNVqodFoghrC2mw2tguOdge63W4MDg5Co9EgJSUFBQUFAYuPTrPKycnhvfUW+j6oIWl9fX3ABe8vsInGtyHewiVqCJuWloaUlBQf8ZdYqsbJyUlMTU2FNXDhC4fDwZIELdnSEC5enxGF3W7Hl7/8ZezZswe33XZbXHZVQbB8iAGY+4K5MBgMmJycRF1dnc/jVMm4bt26kElGekehi6+goAAZGRno6ekRfTgLBZVZ0603n51IsORluMUXb+ES7ScpLCz0UZD6qxqDDevlC5rIDOXkHAs8Hg+mpqbQ398fULIVO1x0OBy44YYbcOmll+Luu+9eKFIAljsxmEwmVh1IMTExgdHRUdYLgU+7tM1mw+joKCYmJpCWloYVK1agoKBAVNciepelRi3RXCT+iy8nJwcFBQXIyspijVVGRkbY8ERs0GaiioqKsCVibolRr9cjPT2d9+IbHh6GwWBgRWdigyom6+rq2CnoNDQKN7xHKFwuF26++Wace+65uP/++xeSFIDlRgxOp9PHxclqtaKnpwdNTU0+Scb6+nq2xZpP5WFqagpjY2PsxUjv0E6nky0vxrI9FiMJ6A+v1wudTscmBalyc+PGjXEhBZqfqampETQLw790Fyo04k62isXJORyoYjKUTJt2XKrV6pi6V91uN2655RY0NTXhW9/61kKTArDciYFOc9qwYQPa2tqQnp7O9qLzlTcPDAzAbDZj/fr1AbGs2+1meyMsFgtyc3NRUFAg6EKhjURiDJcJ9R76+/tZWbZYyUsu/JuhYoHNZmMXn8vlYiXF09PTca1uRCIFf9DuVY1GE1CNCZfz8Hg8uOOOO1BbW4uHH354MUgBWO7E4PF4cPToUTAM45Nk5KNkpP4NKSkpqK2tjfgF0ilPtIFKqVSisLAwbAxNRTl8L0ahCCZcEiN5yQWdfxkPCTINjQYHB+Fyudj29mjyEuFASSHa9+BfjUlKSmJDDu5n6vF4cPfdd6OoqAg/+MEPFosUgOVGDC6Xy8eBSa/X4/jx49iyZYsgJaPD4UBbWxuKiooCWrD5gMqzVSoVG0PT3giaLKOzEhoaGqLuGQgHvsIlmrxUq9Vwu93srodPaEQTmWKIu4KBEMIO9a2urmYl2jqdjrcOIRLoTAwxic1/18MwDIxGIw4fPgylUon/+Z//iUsoJADLlxgmJycxMjICj8eD8847j1eSERB3biTgezehZVCJRAKn08nqJcRGtDkLasFGQyOavFQqlQEXMm2GilcHptfrRUdHB9LT01FVVeXznfnrEKLd9VBSaGhoiNtUJ7fbjVOnTuFb3/oW+vv7sWPHDtx5553YsmVLXF6PJ5YfMXg8HvT19cFisaC+vh7Hjh3D1q1beZECHcser1IevdhtNhsYhmHLoEKdlMNBLOGSf/KSaw2vUqlEb4biglZo6G4nEri7Hr5W+1TLEU9SAOY+x+9+97uwWCx45plncOLECaSkpKCpqSlur8kDy4sYqKZdoVCwScaPPvoIq1evjqiPp5brDQ0NonflAZ9u7ZVKJXsHpMNb1Go1vF4v22UZ7bY8XsIlrjX81NQUvF4vampqUFBQIPpn5fF40NrayvpjCgV1flar1azVvr9xykKRAiEE3/ve9zAzM4Pnn38+XvLmaLC8iIG2JXOTjHq9HpOTkzCZTGy/AZckvF6vj54/HrGfw+FAa2srysrKQk4MooYp0ZZB4y1c4pYLq6qqWJ8BhmHYXU+ssme3242WlhYUFxeLUrYN5gKVkZGB6elpbNiwIS55EQpCCB5//HEMDAzgt7/97VIiBWC5EQPNMQSrPFBRjUqlgsFgYIeETkxMICsrKyCOFQs0ZyHEd9C/DOovVPJHvIVLVJFJx9Fxz8F/Gx+triPcSDex3gMdgCuXyyGTydhzjcd4vaeeegqtra148cUXY84jVVZWIiMjA1KpFDKZDCdOnIBOp8N1112H4eFhVFZW4uWXX0Z2djYIIbjnnntw5MgRpKWl4Te/+Q02btzof8jlRQxutxtOp5OXPoFeJBKJhN1JiOlUBHx6F1+/fn3UGW86+JZbBqVdlhKJxKfDMx6JTCEzHoMlL/m0OVPFZDxnMszOzqK7u5utoND+iFjFSv4ghOC5557Dxx9/jJdeekmUUKuyshInTpzwyRl985vfRE5ODh544AE89thj0Ov1ePzxx3HkyBH89Kc/xZEjR3D06FHcc889OHr0qP8hlxcxHDlyBDU1NSguLg57IXKHvWZmZsJoNEKlUkGr1SItLY0tLcay0KhaUsy7OHdrrNPpQAhBUlJSXJOAtORZWVkp+FypY5HBYAg51zJaxaQQUFJobGwMGu74W+1H689JCMHzzz+Pd999F6+++qpo33swYli9ejX+/ve/o6ioCFNTU/jCF76Anp4e3H777fjCF76A66+/PuB5HCw5YojrJKrJyUk8/vjj8Hg82LFjB3bv3h0w4GV6ehojIyM+cyPpoJGamhqYzWa25TlUh2U40GlXBoMBGzduFPUuTqdQZWVlobOzE8DcjMuTJ09Gda7h4HK52Fkc0bhG0VbmvLw8tmSrUqkwODjInmt6ejq6urqiskrji0ikAICtDhUUFPj4c9KBOHys9gkheOGFF/DnP/8Zr7/+uqghHcMwuPTSS8EwDG6//XbcdtttmJmZYRf7ihUrMDMzA2CuF4ibtC0tLcXExATvadiLhbgSw6233opbbrkF09PTOHToEO666y5YLBbs2LEDV155JQ4fPoxzzz0XmzZtCrpguePUa2pqWG/G06dPQyaTsRdPqAuE69PQ2NgYl0RmsLt4bW2tz7nGWgbl2wzFF9wpT/Rcx8fH0d3dDYVCAaPRiJSUFNE9GwwGA86cOYMNGzbw/hy4MziBT632W1tbAQS32geA3//+9zh06BAOHz4sumDtww8/RElJCVQqFS655BKsWbMm4JwXUUUpCuI+u5JhGBQVFeGuu+7CXXfdBbVajZdffhk7duxAXl4enE4nsrKysHr16ogfpkKhQFVVFaqqqmCz2TAzM4PW1lZIJJKAhed2u9He3s57tmM0oMIl7vzFYOdqt9uhUqnQ2dnJlkGDXczBQHUQ8erdAOburjqdDlu2bEFycjLUajW6u7tjSl76IxpSCAaFQgGFQoHKykq2ckTH1mVlZWF0dBQ6nQ4vvvgi3nrrrbhUOuiOraCgAHv27MGxY8dQWFiIqakpNpSgNgAlJSWspSEwNxclHj6hYiOuOYZQeOCBB1BWVobrr78ef/rTn3Do0CFMTEzgsssuw549ewSXKenCo/qD7OxsaDQaVFRUxG3LFu2C9S+D5ubmorCwMOjCo0rAeJqT0N6KYBoCWo2hE574Ji/9QedWxEoK4eDxeDA6Oop77rkHJ06cwI4dO3DjjTfi0ksvFfV1LBYLvF4vMjIyYLFYcMkll+C73/0u3n//feTm5rLJR51Ohx/96Ed466238Mwzz7DJx6997Ws4duyY/2GX3PZiUYjB6/UGXFizs7N44403cOjQIQwODuKSSy7B7t27BYcAdPaiXC4HwzCsSElMLYFYwiUq/pmZmQkogxqNRnR1dcVtHgPw6V2cT2+Ff3k5IyODV+WIjgWMl1Sbi7feegs/+clP8MYbb6Cvrw9TU1OiOzsPDg6yx3S73fjyl7+MBx98EFqtFnv37sXo6CgqKirw8ssvsxPQDxw4gLfffhtpaWn49a9/jU2bNvkfNkEMfGA2m3HkyBEcPHgQZ86cwUUXXYRdu3Zh8+bNYUnCX1TkcrnYu7PD4RBlWxyvRiVaBqUVDrfbjVWrVmHFihVxyY3Q9xHNXZybvNRqtSETgjqdDn19fQtCCn/5y1/w2GOP4a233opbyBVHJIhBKGw2G9555x0cPHgQra2tuOCCC7Br1y5s27bN5041OTmJ8fHxkOVIrkjJarWyXYvh9Pv+oB2Y8RIuAZ/OwKysrITBYIBOp2O7QSN5C/CFRqPBwMCAaAuW2zZOd2lJSUmsMWy8SeGvf/0r/vu//xtHjhyJm+4izkgQQyxwOBx477338Morr+DEiRP43Oc+h507d+If//gHPv/5z+Pzn/88rzo3rZOrVCqYTKaISkYAcRcuAXPkNjk5icbGRrbESR2V6KQsWlqkA1GEgqoyI9m7RwuHw4Hh4WFMTk76GNBkZGTEJQH8wQcf4MEHH8Rbb72FFStWxHw8j8eDTZs2oaSkBG+++SaGhoawb98+aLVaNDc34//+7/+QnJwMh8OBG2+8ESdPnkRubi5eeuklwdoSDhLEIBZcLhfeffdd3HfffZBKpdi0aROuuuoqfP7znxd0wVPhz8zMDIxGI7KysnzMR6jjkt1ux7p16+LWtz8yMgKdToeGhoaw5EbLoFxTXL5l0KmpKUxMTPgQj9jQarXsbkQikbAEbDabg/bGxIJ//etf+M///E+8+eabomX6n3zySZw4cQJGoxFvvvkm9u7di6uuugr79u3DHXfcgcbGRtx555147rnn0NbWhp///Of44x//iNdeew0vvfRStC+bIAYx8dxzz8Fut+Puu+/Ghx9+iIMHD+KDDz5AY2Mjdu/ejYsuukhQ/Oxv6JKRkQG73Y6MjIyQo+JiBbWrs9lsgomHW43xeDxhE63j4+OYmZnBhg0b4tZARN2dgu1Gok1ehsLx48dxzz334PDhwygvLxfl/MfHx3HTTTfhwQcfxJNPPok33ngD+fn5mJ6ehkwmw7/+9S888sgjeOedd3DZZZfhkUcewbZt2+B2u7FixQo2lIoCCWIQE7QpiwuPx4OPP/4Yr776Kt5//33U1dVh9+7duOSSSwQlC10uF06fPg2GYeB2u4O6Polx/mfOnAHDMLx0HOFAx7/NzMz4JFozMjIwNjYGrVYbcTcSC9RqNYaGhniFKP5mOXK5nN358NntnT59Gl/96lfx+uuvh51qJhTXXHMNvvWtb8FkMuGJJ57Ab37zG5xzzjno7+8HMGcF8MUvfhEdHR1Yv3493n77bdZRbOXKlTh69Gi0nhtLjhjiLnCKJ4ItJKlUivPPPx/nn38+vF4vTpw4gVdeeQWPP/44Vq5ciZ07d+Lyyy8P6+lIhUu0LZvG+TMzMxgaGkJqaip7IUebb6AmMWlpaaJMy05OTmbbo2kZlIYnUqkUa9eujVsYREmhqamJV4jCVV5SRStVM9LkZShL+Pb2dtx555149dVXRSWFN998EwUFBWhubsbf//530Y57tuKsJoZIkEgk2LJlC7Zs2QKv14vW1la88sorePrpp1FaWoqdO3fiiiuu8HFSDiZcYhgGmZmZyMzMZC/kmZkZnDx5EsnJySgsLBTUE0Fl1HR6ltigcnGaWC0sLMTMzAx6e3sDukFjhUqlwsjICG9SCAaumpF2Wfb09LACMJq8PHPmDP7jP/4DL730Emv8IxY++ugjHD58GEeOHIHdbofRaMQ999wDg8EAt9sNmUzmo1qkisbS0lK43W7Mzs6ejWXSkDirQ4loQQhBR0cHDh48iCNHjiAvLw+7du1CYWEhWlpacOeddyIzM5PXsbjJQLog8/PzQ5boaDOUWOYnwUBndng8Hh+/BkIIm0MRowxKSSFe3aR059Pe3o57770XbrcbDz30EG6++ea4JU8B4O9//zueeOIJvPnmm7j22mtx9dVXs8nHhoYGfPWrX8Wzzz6L9vZ2Nvl46NAhvPzyy9G+5JILJZYlMXBBTU8ef/xxHD58GJs2bWKbvPLz8wVt8W02G1QqFVQqFeukVFBQwCZAqXNUZWVlXEbq0ffT3d0NqVQaNmHqXwYVGufPzMywE8TiuUiBObXhv/3bv+GOO+5gDWkfe+yxuL0elxgGBwexb98+6HQ6NDU14Xe/+x3kcjnsdjtuuOEGnD59Gjk5OfjjH/+I6urqaF8yQQxLEXq9Htdccw1efPFFmM1mvPrqq2yr7pVXXoldu3ZhxYoVgkjC4XCwJOH1epGVlQW1Wo01a9bwdo4SCiEmLv6gOx86bZpWOIJVdSgpNDU1xU3TQTEyMoJ9+/bhV7/6VTAp8WcFCWJYqvCvcNAp14cOHcJrr70Gr9eLK6+8Ert370ZpaamgRUf7N1JTU0EI8ZFmiwWv14v29nZkZmbGnJTzH8zLLYNOT0+zjtTxJoWJiQlce+21+NnPfoZt27bF9bUWGQliOBtBCMHU1BRLElarFV/60pewa9euiANvqTsV7V6k05xUKhVsNptPWTHaygRNZubm5opW06egZVCqEgWA9evXh1WJioGpqSlcc801ePrpp3HBBRfEdCy73Y4LLrgADocDbrcb11xzDf77v/97oVSNfJAghrMdhBCo1Wq89tprOHToEPR6Pb74xS9i9+7dATG9VqtFX19fWAszqj2IdtZmqDH3YoOqJouLi6HVatk27EhS8mgwMzODq6++Gk888QQuuuiimI9HxwCmp6fD5XLhvPPOw9NPP40nn3xyIVSNfLB8iOHtt9/GPffcA4/Hg1tvvRUPPPBAtIda0tBqtaynxPT0NC699FLs2bMHLS0tyM7Oxvbt23kl8+iszZmZGR9r/ezs7JCLjjo5BzOKEROTk5OYmpryUU36D8Ch8yJyc3NjKoNqNBpcddVV+P73v4/LLrtMrLfAwmq14rzzzsPPfvYzfOlLX1oIVSMfLDliiEuQ6PF4cNddd+Hdd99FaWkpNm/ejJ07d6Kuri4eL7eoyM3Nxf79+7F//34YDAa88cYbuO2226DVanH11VejuLgYDQ0NERcLd6wblQ9PT0+jp6cn6EBeavcWzwoHMBfnB5NS+3tI0jJof38/FAoFCgsLBZdBdTodrr32WjzyyCOik4LH40FzczP6+/tx1113YeXKlcjKymLPj3ox0vdMfRplMhmUSiW0Wm1Mk8TONsSFGI4dO4aamhq2fLNv3z786U9/+kwSAxe0Aau8vBzvvPMO/va3v+EnP/kJent7WU+JTZs2RSQJiUSC3Nxc5ObmghDC9hj09vYiIyMDOTk5GB0djXnsXSRQUmhsbAwrpWYYBtnZ2ewcBVoGHR4eRnJyMq8y6OzsLPbu3YsHHngAO3bsEP29SKVStLS0wGAwYM+ePawXaALBERdiCOaMG8RL/zOJ7du346KLLkJSUhKuu+46XHfddbDZbHj77bfxy1/+EnfffTc+//nPY9euXTjnnHMi9i4wDIOcnBzWDUitVqOrqwsymQyTk5PweDyi9m9QjI+PQ6VSRSSFYOfrrxKlcudQZVCTyYS9e/fi3nvvFd1xyR9ZWVm48MIL8a9//WvZqhr54DMtiV4MSKXSgIWUmpqKPXv2YM+ePXA4HHj33Xfxu9/9Dl//+tdx7rnnYvfu3Tj33HMjbrttNhvbvahUKtk7M7d/I9Yx9MBcs5BGoxFMCsHAlTvTMmhnZyc8Hg873+J73/sebr/9duzduzem1woFtVqNpKQkZGVlwWaz4d1338V//dd/4cILL8TBgwexb98+vPDCC9i1axcAYOfOnXjhhRewbds2HDx4EBdddNFZ7/osFHEhhrPVGXchIJfLsWPHDuzYsQNOpxN/+9vfcPDgQXzjG9/Ali1bsHv3blxwwQUB226LxYK2tjasW7eOlWtz78xmsxkqlQqnTp3ivX0PBi4pBAt5To0a8MR7AzDZ3bhifQFuP78SEp6LJiUlBWVlZSgrK4PT6cTrr7+Ob3/727Bareju7sbIyEhcekempqZw0003wePxwOv1Yu/evdixYwfq6uqwb98+fOc730FTUxNuueUWAMAtt9yCG264ATU1NayqcbkhLlUJ6lf4/vvvo6SkBJs3b8aLL76IdevW8fp7IbMBPytwu9344IMPWE+JpqYm7N69GxdeeCG6urpgNBqxadMmXqIoq9XKqi65w1siWaxR6/VQydI+lRnX/eok7C4vACA1SYKvbC7FfdtXCn6/drsd119/Pa666ip85StfwTvvvIP6+nqsWrVK8LE+A1hy25G4lSuPHDmCe++9Fx6PB/v378eDDz7I+2+FzAb8LMLj8eCjjz7Cq6++ijfffBMOhwMPPPAA9u3bJ9iAlpq5qFQqEEJYkvDXVYyMjECv14etoPz0b4P4+T9HfC6K7LQkfPSf5wk6J4fDgRtuuAGXXXYZDhw4sOy26UGw5D6AJSlwEjIb8LOMEydO4Pbbb8dDDz2EDz/8EO+++y5qamqwe/duXHbZZYIl1bSlWaVSwe12s4lAjUYDg8GA+vr6sBWTn38wjGc/GIbH++llkZ+ejH/cdy7vc3C5XLjppptw/vnn47777kuQwhyW3IewJImhqqqKFfbQ2YBZWVkwGAxzJ0UIsrOz2X9/VjE7Owur1cqKl7xeL1paWvDKK6/gnXfeQVlZGespIXS+BbXWHx4ehsPhQFlZGVasWAGFQhFysU4b7dj18+OwONzwEiAlSYL/uqQG123ilz9yu93Yv38/mpub8cADD8RMCmNjY7jxxhsxMzMDhmFw22234Z577ol1JP1iIEEMfDAxMeEzG/CnP/0pdu7c6UME2dnZ0Ov1i3F6SwLUU+KVV15hbdN37dqFHTt28O7eHBoagslkwpo1a1jTVpvNxk7HCta/MWGw4fmPRjFrc+GKdQXYvpafuMrj8eD222/HqlWr8PDDD4uyU5iamsLU1BQ2btwIk8mE5uZmvP766/jNb34Ty0j6xUCCGITikUceQXp6On75y18uu1CCL6h35MGDB/Hmm29CqVRi586duPLKK5GXlxd0EQ4ODsJisQQY0NL+DersTB2ghPRv+MPj8eDuu+9GcXExvv/978ctfNi1axcOHDiAAwcOxDKSfjGw5IghPiaAMcBisbBdfBaLBX/5y1+wfv16trYMwKfm7I/9+/ejoKAA69evZx/T6XS45JJLUFtbi0suuYTdaRBC8LWvfQ01NTVoaGjAqVOn4vzu4gOGYbB27Vo89NBD+OSTT/Dcc8/BYrHg+uuvx44dO/C///u/mJ6eBr0J9PT0BCUFYE6HUVhYiPr6emzduhU5OTmYmJjAJ598gu7ubuh0Oni9Xt7n5vV6cd999yE3NxePPvpo3EhheHgYp0+fxtatWwWPpBeC7373u3jqqafYfz/44IN4+umnY38DSw2EkHA/C46BgQHS0NBAGhoaSF1dHXn00UcJIYRoNBpy0UUXkZqaGnLxxRcTrVYb9O//8Y9/kJMnT5J169axj33jG98gP/zhDwkhhPzwhz8k3/zmNwkhhLz11lvk8ssvJ16vl/zrX/8iW7ZsifO7W1h4vV4yPDxMfvzjH5Pzzz+fnHfeeWTHjh3kzjvvJCaTiVgsFt4/JpOJjI6OkhMnTpD33nuPHDt2jIyMjIQ9jslkIl/96lfJgQMHiMfjidv7NJlMZOPGjeTVV18lhBCiVCp9fp+VlUUIIeRLX/oS+ec//8k+ftFFF5Hjx48Leq2hoSHS1NRECCHE4/GQ6upqotFoYjh7Qkj4NbgoP0uOGMTA0NCQDzGsWrWKTE5OEkIImZycJKtWrSKEEHLbbbeRF198MejzPmvweDzk3nvvJVu3biUXXnghOeecc8ijjz5K2tvbidlsFkQSZrOZjI+Pk1OnTpH33nuPHD16lAwNDRGj0ehDCl//+tfJ7bffHldScDqd5NJLLyU//vGP2cfi/X1v376dnDp1ivz5z38mV199dYzvgBCyBIjA/2dZSKKFbi2XQMwpOlwuF4qKivDxxx+DYRioVCq89tpruO+++2AwGHDFFVdg9+7dqK2tjbjd92+aogNuBwcH4fF40NXVhdHRUej1ejz//PNxs60nhOCWW27B2rVrcd9997GP07DzgQceCJA6P/PMM9i3bx+OHj0KpVIZ1Xd966234je/+Q2mp6exf/9+0d7PkkIE5jgr4b9jiOfW8rMAjUZDnn/+eXLFFVeQ5uZm8u1vf5scP348qp1EZ2cnueKKK0h+fj7ZsWMHee211+J23v/85z8JAFJfX08aGxtJY2Mjeeutt0KGnV6vl3z1q18l1dXVZP369VF/1w6Hg6xatYpUVVURt9stxltZ9B2C/8+y2DEUFhZiamqKzVJT/4JET8cccnNzccstt+CWW26BwWDA4cOH8eijj2JkZASXXHIJ9uzZE1H8RHHo0CEoFApMTk5iYGDA5/MVG+edd95cPBwE77//fsBjDMPg2Wefjfl1k5OTceGFFyIrKytuk70WG0uuKhEPhKpo7Ny5E7/97W9BCMEnn3wS9dbys4SsrCzceOONeP311/HBBx9g48aNePLJJ3HuuefiO9/5Do4fPx60KkEIwbPPPouTJ0/id7/7HWQyGVavXo3t27cvwruIL7xeLz755BO26eoziQhbirMO+/btIytWrCAymYyUlJSQ559/Pqqt5b//+7+T/Px8n5Dk4YcfJsXFxT7bVoof/OAHZOXKlWTVqlXk7bffXrg3vECwWCzk4MGD5Prrryf19fXkwIED5N133yVGo5GYzWbyk5/8hHzpS18idrt9sU81rujs7CRVVVXkvvvuE/Owix46+P985ohBLAQrez788MPkf/7nfwKe29nZSRoaGojdbieDg4OkurparNhzScJms5HDhw+TG2+8kaxfv55ccMEF5MILLyQ2m02U4wcjZa1WS7Zv305qamrI9u3biU6nI4TMkfvdd99NVq5cSerr68nJkydFOYcFxqITgf/PsgglosEFF1zAW1r8pz/9Cfv27YNcLkdVVRVqampw7NixOJ/h4iElJQVXXnklXnjhBZw8eRK7d+/Gyy+/HHQ4TTS4+eab8fbbb/s89thjj+Hiiy9GX18fLr74YnYS1Z///Gf09fWhr68Pv/jFL3DnnXeKcg7LHQliEIhnnnkGDQ0N2L9/P6ugFENRd7YiOTkZX//610X1ngxGyn/6059w0003AQBuuukmvP766+zjN954IxiGwTnnnAODwYCpqSnRzmW5IkEMAnDnnXdiYGAALS0tKCoqwv3337/Yp7RsEE+ZcwKBSBCDABQWFkIqlUIikeA//uM/2HAhUfZcWDAMk/BxiDMSxCAA3C3qa6+9xjZq7dy5E3/84x/hcDgwNDSEvr4+bNmyhX3u2NgYLrzwQtTV1WHdunVs081nvblLTFAtCoCEFmUBkCCGELj++uuxbds29PT0oLS0FL/61a/wzW9+E/X19WhoaGBnRgDAunXrsHfvXtTV1eHyyy/Hs88+6yN8kclk+PGPf4yuri588sknePbZZ9HV1ZVIqAlAQouywIhQtkggDti5cyf5y1/+kmjuCgGxtChnERa9POn/sywk0UsJsfgGLJc74R/+8Iegj8dT5pyALxKhxALCbDbj6quvxlNPPcXOhqBIJNQSWEpIEMMCweVy4eqrr8ZXvvIVXHXVVQASCbUEli4SxLAAICS8bwCQSKi9/fbbWL16NWpqatgkbAKLiAhJiAREgJi+AaOjo+QLX/gCWbt2LamrqyNPPfUUIeTsbvByu92kurqaDAwMEIfDQRoaGkhnZ+din9ZCYtGTjf4/CWI4yzA5Ock2ChmNRlJbW0s6OzvP6gavjz/+mFx66aXsv3/wgx+QH/zgB4t4RguORScC/59lH0ocP34cDQ0NsNvtrHNyR0fHYp9WSBQVFbFDUjIyMrB27dqwEuCzocErIWteeog0V2JZgGGYRwGkAEgFME4I+eEinxIvMAxTCeADAOsB3AfgZgBGACcA3E8I0TMM8wyATwghv5v/m18B+DMh5OCinHQQMAxzDYDLCSG3zv/7BgBbCSEHFvfMli+W/Y5hHv8PwCUANgH40SKfCy8wDJMO4FUA9xJCjAB+BmAlgA0ApgD8ePHOTjAmAJRx/l06/1gCi4QEMcwhF0A6gAzM7RyWNBiGScIcKfyeEHIIAAghM4QQDyHEC+CXAGizxtmw6I4DqGUYpophmGQA+wAcXuRzWtZIEMMc/hfAQwB+D+DxRT6XsGDmVFC/AtBNCHmS8zi3nrkHAE2UHAawj2EYOcMwVQBqASypJAMhxA3gAIB3AHQDeJkQ0rm4Z7W8sewl0QzD3AjARQh5kWEYKYCPGYa5iBDy18U+txA4F8ANANoZhmmZf+zbAK5nGGYD5uaNDgO4HQAIIZ0Mw7wMoAuAG8BdhBDPAp9zRBBCjgA4stjnkcAcEsnHBBJIIACJUCKBBBIIQIIYEkgggQAkiCGBBBIIQIIYEkgggQAkiCGBBBIIQIIYEkgggQAkiCGBBBIIwP8PQB5n2tACWUYAAAAASUVORK5CYII=\n",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {
|
|
"needs_background": "light"
|
|
},
|
|
"output_type": "display_data"
|
|
},
|
|
{
|
|
"name": "stdout",
|
|
"output_type": "stream",
|
|
"text": [
|
|
"done!\n"
|
|
]
|
|
}
|
|
],
|
|
"source": [
|
|
"import numpy as np\n",
|
|
"import scipy.constants as sc\n",
|
|
"import matplotlib.pyplot as plt\n",
|
|
"\n",
|
|
"class Ion:\n",
|
|
" def __init__(self,position, charge, mass, v_z):\n",
|
|
" self.position = position\n",
|
|
" self.charge = charge\n",
|
|
" self.m = mass\n",
|
|
" self.v_z = v_z\n",
|
|
" self.F = np.array([[0],[0],[0]])\n",
|
|
"\n",
|
|
"class infiniteLengthQuadrupole:\n",
|
|
" def __init__(self,R,r_0):\n",
|
|
" # rod radius\n",
|
|
" self.R = R\n",
|
|
" # rod placement radius\n",
|
|
" self.r_0 = r_0\n",
|
|
" # rod positions\n",
|
|
" self.pPole = np.array([[[ (R + r_0)/np.sqrt(2)],\n",
|
|
" [(R + r_0)/np.sqrt(2)],\n",
|
|
" [0]],\n",
|
|
" [[-(R + r_0)/np.sqrt(2)],\n",
|
|
" [(R + r_0)/np.sqrt(2)],\n",
|
|
" [0]],\n",
|
|
" [[-(R + r_0)/np.sqrt(2)],\n",
|
|
" [-(R + r_0)/np.sqrt(2)],\n",
|
|
" [0]],\n",
|
|
" [[(R + r_0)/np.sqrt(2)],\n",
|
|
" [-(R + r_0)/np.sqrt(2)],\n",
|
|
" [0]]])\n",
|
|
" # rod pseudo charges\n",
|
|
" self.rodsPseudoQ = np.zeros(4)\n",
|
|
" # K factor for simplification ofcalculating the equations\n",
|
|
" self.K = 1 / (4 * np.pi * sc.epsilon_0)\n",
|
|
" \n",
|
|
" self.forces= []\n",
|
|
" \n",
|
|
" def check(self):\n",
|
|
" print('pole positions')\n",
|
|
" print(self.pPole)\n",
|
|
" print('pseudo charges')\n",
|
|
" print(self.rodsPseudoQ)\n",
|
|
" \n",
|
|
" def phi_0(self,U, V, f, t):\n",
|
|
" return U + V * np.sin(2 * np.pi * f * t)\n",
|
|
" \n",
|
|
" def coulombForceOnQ1(self,Q1,Q2,r1,r2):\n",
|
|
" a = self.K * Q1 * Q2\n",
|
|
" r12 = r1 - r2\n",
|
|
" mag_r12 = np.linalg.norm(r12)\n",
|
|
" # calculate force\n",
|
|
" return a * r12 / (mag_r12**3)\n",
|
|
" \n",
|
|
" \n",
|
|
" def calcNewIonPos(self,U,V,f,ion,t_s,t):\n",
|
|
" \n",
|
|
" # check if ion has left the r_0 boundary\n",
|
|
" if np.linalg.norm(ion.position) >= self.r_0:\n",
|
|
" return None\n",
|
|
" \n",
|
|
" signLUT = [1,-1,1,-1]\n",
|
|
" index = 0;\n",
|
|
" \n",
|
|
" # reset force\n",
|
|
" F_c = np.array([[0],[0],[0]])\n",
|
|
" \n",
|
|
" for Q in self.rodsPseudoQ:\n",
|
|
"\n",
|
|
" \n",
|
|
" # calculate the pseudo charge points of the quadrupole\n",
|
|
" Q = signLUT[index] * self.phi_0(U,V,f,t)/2 * self.K * self.R\n",
|
|
" \n",
|
|
" # match pseudo charge with ion in z direction\n",
|
|
" self.pPole[index][2] = ion.position[2]\n",
|
|
" \n",
|
|
" #print('ion position', ion.position)\n",
|
|
" #print(self.pPole[index])\n",
|
|
" \n",
|
|
" # calculate force on ion\n",
|
|
" F_q = self.coulombForceOnQ1(ion.charge,\n",
|
|
" Q,\n",
|
|
" self.pPole[index],\n",
|
|
" ion.position)\n",
|
|
" #print(F_q)\n",
|
|
" \n",
|
|
" F_c = F_c + F_q\n",
|
|
" self.forces.append(F_c)\n",
|
|
" \n",
|
|
" # update index\n",
|
|
" index = index + 1\n",
|
|
" \n",
|
|
" ion.F = F_c\n",
|
|
" # http://lampx.tugraz.at/~hadley/physikm/problems/rvaF/index.en.php\n",
|
|
" newPosition = 0.5 * t_s**2 * 1 / ion.m * ion.F + ion.position\n",
|
|
" \n",
|
|
" # velocity component in z\n",
|
|
" newPosition[2] = newPosition[2] + ion.v_z * t_s\n",
|
|
"\n",
|
|
" #print(newPosition)\n",
|
|
" \n",
|
|
" return newPosition\n",
|
|
" \n",
|
|
"N_periodes = 100\n",
|
|
"N_points_per_period = 101\n",
|
|
"f = 1.1e6 \n",
|
|
"t_step = 1 / f / N_points_per_period\n",
|
|
"t_sim = 0\n",
|
|
"N = N_periodes * N_points_per_period\n",
|
|
"\n",
|
|
"# coffein atom 524.50002164 dalton (8.709527e-25 kg)\n",
|
|
"\n",
|
|
"\n",
|
|
"ion = Ion(np.array([[1e-10],[1e-10],[0]]), sc.e, 1.661e-27,1)\n",
|
|
"\n",
|
|
"# quadrupole\n",
|
|
"quad = infiniteLengthQuadrupole(0.0101, 0.01)\n",
|
|
"\n",
|
|
"#quad.check()\n",
|
|
"\n",
|
|
"positions = []\n",
|
|
"\n",
|
|
"\n",
|
|
"positions.append(ion.position)\n",
|
|
"print('N =',N)\n",
|
|
"for n in range(0,N-1,1):\n",
|
|
" \n",
|
|
" newPos = quad.calcNewIonPos(2.4764,17.335,f,ion,t_step,t_sim)\n",
|
|
" \n",
|
|
" if(newPos is None):\n",
|
|
" print('ion has left the boundary!')\n",
|
|
" break\n",
|
|
" \n",
|
|
" \n",
|
|
" ion.position = newPos\n",
|
|
" \n",
|
|
" if n % 100 is 0 and 0:\n",
|
|
" positions.append(newPos)\n",
|
|
" \n",
|
|
" if 1:\n",
|
|
" positions.append(newPos)\n",
|
|
" \n",
|
|
" t_sim += t_step\n",
|
|
"\n",
|
|
"print('n =',n)\n",
|
|
"print('simulation time = ', t_sim)\n",
|
|
"print('number of positiions = ', len(positions))\n",
|
|
"\n",
|
|
"\n",
|
|
"if 0:\n",
|
|
" print(positions[:10])\n",
|
|
"\n",
|
|
"if 1:\n",
|
|
" ax = plt.axes(projection='3d')\n",
|
|
" #ax.scatter3D(positions[:][0], positions[:][1], positions[:][2], c=positions[:][2], cmap='Greens')\n",
|
|
" ax.scatter3D(quad.forces[:][0], quad.forces[:][1], quad.forces[:][2], cmap='Greens')\n",
|
|
" ax.set_xlabel('x')\n",
|
|
" ax.set_ylabel('y')\n",
|
|
" ax.set_zlabel('z')\n",
|
|
" plt.show()\n",
|
|
"\n",
|
|
"print('done!')\n",
|
|
"# dummy check\n",
|
|
"if 0:\n",
|
|
" ion_1 = Ion(np.array([[0.0000005],[0],[0]]), 1000 * sc.e, 1,0)\n",
|
|
" ion_2 = Ion(np.array([[-0.0000005],[0],[0]]), 1000 * sc.e, 1,0)\n",
|
|
" print(quad.coulombForceOnQ1(ion_1.charge,ion_2.charge,ion_1.position,ion_2.position))\n",
|
|
" "
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 9,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"name": "stdout",
|
|
"output_type": "stream",
|
|
"text": [
|
|
"[array([[1.e-100],\n",
|
|
" [1.e-100],\n",
|
|
" [0.e+000]]), array([[1.00000000e-100],\n",
|
|
" [1.00000000e-100],\n",
|
|
" [2.27272727e-007]])]\n"
|
|
]
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAQcAAADyCAYAAABAgwC5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABjTklEQVR4nO19eXgb5bn9+bR439ckdrzFMdm8JoYA/QVKSUmhWUqhhLbsFG7LWsq9vVBKuTyX0hZKS6GF27KWJaFJSEMJpKVlLRCSQLzEduJ4i3dJtmxLtvbR9/vD/iYjWcuMpJHlZM7z5IE40mgka868y3nPSyilUKBAgQJvqOb6BBQoUBCbUMhBgQIFPqGQgwIFCnxCIQcFChT4hEIOChQo8AmFHBQoUOATmiD/rvQ5FSiQH2SuT8AXlMhBgQIFPqGQgwIFCnxCIQcFChT4hEIOChQo8AmFHBQoUOATCjkoUKDAJxRyUKBAgU8o5KBAgQKfUMhBgQIFPqGQgwIFCnxCIQcFChT4hEIOChQo8AmFHBQoUOATCjkoUKDAJxRyUKBAgU8o5DAHoJTC4XDA5XJBWQ2gIFYRzOxFQYThdrvhcDhgs9n4n6nVami1Wmg0GqjVahASk94fCk4zkCB3LuW2FiFQSuFyueByuUAIgdPp5H9OKYXb7eZJwW63IzU1FXFxcQpZnB6IyV+wEjlEASyNEBIAAyEEhBCoVCr+sZ2dnSgpKUFSUhIAJbJQMDdQyEFmuFwuPkpgRMCiBV8XOXuMWq2GWq3mowqr1co/XqPR8H8UslAgFxRykAneaUSoF7CvyILjOLhcLv4xGo2GjyxUKpVCFgoiAoUcZIDb7YbT6eTTCKkXK4su/P2b8HjeZEEI8YgsFLJQECoUcogg2IXK0gh2tw/0+HAvXF9k4Z3KqFQqaDQaxMXFKWShQDQUcogQKKVwOp3gOC5otCCsPfj791D1D77Ior+/HwCwcOFCJbJQIBoKOUQAbrcbJ06cQG5uLjQajeiLbdLuQkO/CRNWF/LT4lFVkIZ4TWR1aUKyYAVOp9PpEVmwmoVarVbIQgEPhRzCgLDoODg4iJycHNEXlpNz44PjRjg4N1LiNegetcDi4LCuPCusyCEYWCdE+B68yYIVNxnRKWRxekIhhxDhrV0ghMDtdot+vsnmwpTDhfzUBABAXko8hk12ODjKHz+S8HeB+yILh8MBu90OYLpuotVq+chCIYvTBwo5hABWdGQFxWA1BF/QqAjcFHBTChUhcLndIARQEf8XcjQQiCyExU1hGqLg1IRCDhLgrV0QXhgqlUoSOaQlaLAkJwkdBgs0KgKOUqwpzoBWrYp4WhFOV0RIFuycHA4HGhoasGzZMj79UMji1INCDiIRTLsgNa0ghGBNUToWZybC6uCQmqBBbmp8pE87ohAWNm02G08aDocDDocDAJTI4hSCQg5B4K1d8JdzS73bT0xMID4+HgvTE8I+1lxBKPUGTg6ReZOFcC5EIYv5A4UcAkCqdkHMBc1xHNra2mC1WuF2u+FyuZCWlobMzExkZmZCq9VG8i3ICn9DZAyMLOx2+6wCp0IWsQ+FHPyA+S4Ii46BIIYcJicn0dzcjMLCQixduhTA9AU0MTGBsbEx9Pf3w+12g1KK+Ph4pKSkQKOZv7+iYGRBKfVIQaRoRBTIj/n7zZMJwjTCu+gYCMEKkoODg+jp6cGqVauQlpbGE49KpeKjBmB6irO1tRVmsxkNDQ1QqVTIyMhAZmYm0tPT5/Wd1hdZuN1u3vhmcHAQhYWFiIuLUyZOYwAKOQhAKYVer0d8fDwSEhIkfTH9RQ4sjeA4DmeeeSYfCfh7vEajQVJSErKyspCVlQWn04mxsTHo9Xp0dHRAq9XyZJKamjqvLx5vshgeHsaiRYsUl6wYgUIOM2BpxNDQEHJzc5GYmCjp+b66FcI0orCwMKQvtlarRV5eHvLy8gAANpuNT0EmJyeRkJCAzMxMZGVlISkpad5fPN4aC28vC4UsoofTnhy8tQts/kAqvCOBgYEBnDhxApWVlUhNTQ3rWEIkJCRg4cKFWLhwISilsFqtGBsbQ3d3N6amppCSksJHFlIJLtbgy8tCIYvo4bQmB1/aBaliJgb2PI7j0NraCkqpRxohBwghSEpKQlJSEgoKCkApxeTkJMbGxtDe3s57UVJKJRNULEIMWSguWZHDaUkOgbQLUsVMDIQQWCwWdHR0YPHixSgoKAhLlRhq9JKamorU1FQUFRXB7XbDbDajp6cHAwMDGB4eRnp6Oh9ZzOdOCODfJYsVkxni4uIQHx+vTJxKxPz+doSAYPZtKpUqJHIwm80wmUyoq6uLmbu0SqVCeno6srKyoNFokJeXx7dNe3t7AcCjEyLM9+cjfJFFb28vX7cRjqcrXhbBcVqRgxj7NqlphcvlQltbGywWC5YsWRIRYpBLIalWq/kuCDB97uPj4xgZGUFnZyfUajVf3ExNTZ3XbVPg5OfI0gzFUk8aTgtykKJdkJJWTE5OoqmpCUVFRUhOTo7YxRQt+bRGo0FOTg5ycnIATM9IjI2NYWhoCMeOHUN8fDyfgqSkpMzLC8ftdvO/FzGWegpZnMQpTw5SJNDAdOTAcVzQ4/b396O3t5fvRpw4cSKkdCSWEBcXh/z8fOTn5wMA3wnp7e3F5OQkkpOTebKYD7MfgCc5eMMXWSguWSdxSpODVAk0ME0O7MvhC0zBCECUqCkUxMrgVWJiIhITE7Fo0SJQSmGxWDA2NoaOjg5YLBa0trbyaUh8fHgTpXK930Dk4A0pLlnzmSwIIYSK+MBPSXII5LsQDIFqDmazGc3NzSgqKkJhYaHHv8XKBS0XCCFITk5GcnIyCgsLceDAARQWFmJsbAxtbW1wOBwenRCpA2SRcOL2BSnk4I1AxjeDg4PIz89HUlLSvLHUI4SkUkrNYogBOAXJgf0CP/vsM5x55pmSf1m+uhWUUgwMDHikEd6Q2gIN1BWZD0RDCEFaWhrS0tJQXFwMt9s9a4CMdUrS09ODtk1jkRy8ISQLo9GI/Px8D5csFlnEsJfFfYSQYwAOAOiklFoDPfiUIgdhcYnVGKTC+yJnaQQhJKCoSaVSeWyhCoZYv/ilwtcAGSOL7u5uEEL4f/c1QBbJizgax+U4jl9ZCHi6ZDkcDtx666247777sGzZsoi/dhi4GUAHgEsAfEoI+SeAHkrpuK8HnxLk4E+7EMrdSJhWsDSiuLgYBQUFAZ8n9W4/OTkJt9uN5ORkn74I0TKYlQsajQbZ2dnIzs4GgFkDZBqNxqNtOh8iB+/jClMOoZwbmI4sYlC+3kIpPZcQcj6AawH8AMDHhJBdAP7uHUnMe3Lwp11gYbtUYQ/rVrBuRFVVFVJSUoI+T2xaQSnFiRMnMDQ0hLi4OFgsFqSmpiIrKwuZmZlhF/aigVCIy3uAzG63w2g08gNkcXFxsNvtfFckUkQhFzkAgQmXzbnECgghGgA5AEApfR/A+2T6DVwD4C4AjwMoFj5n3pJDsNVzoZKD2+3GyMjIrBHrYBBzt3e5XGhubkZ8fDxWr17NE5rZbIbRaERraytcLhfUajVSUlKQl5cXEdVipKOQSNzl4+PjPQbIWBekp6cnogNkcpFDsM/UYrHEFDkA0ALYDwCEEPaBUErpCwBeIIQs8H7CvCQHMdoFtVotWXdgNpvR2tqK+Ph4VFZWSnpuMGUlS1FKS0uxcOFCuFwunhxYYa+kpAQcx/GtwsOHD/O5/KmiWvQFQggSEhKQlJSEVatW+R0gY2QhJbqSM3IIBKfTibi4uKi/bgDYANwDAJTSWRcGpXTY+2fzjhzEahfEipmAk92Ivr4+nHHGGRgenvU5BUWgtIKNb4tJUVjUkJKSgoKCAl61ODg4CLPZjISEBF4CnZiYOCetMznqA95KRl8DZGNjY3x0xdqmGRkZAdumzG0r0ojllqUvzLQvB9nfZ9KMLEqp3t9z5g05SNUuiB2gcrlcaGlpgVqtxplnngmbzRbyVKZ35MBxHI4ePQqXyxXy+LZQtcj8G4xGIzo7O2G1Wvk7alZWVtTuVHKQQ6BjsgGy9PR0ProSDpBRSj06IXIPkAVLKeQqrkYYZwD4NYANhBCVr2hiXpCDr9VzwSCGHFioX1JSgkWLFol+ni94k4PFYkFTUxMWLlyIoqIiSV8Wf1GI0L+hsLDQ44565MgRcByHjIwMZGVlISMjQ7aLRI42rJQ7vJQBMjnOlbUxAyFWCUKgjhSKdXyeaMyTg6/Vc2IQ6CKnM2vp+/v7Z4X6oZq9CMnBYDCgvb0dK1euREZGRkjHEgNfd9Tx8XEYjUZ0dXXx7UKHwxFxy3s50opQjxlogMxisaCxsTGiA2TByIEVlWMcqQDmpwiKUoqpqSlelSY1b/RXc/BOI7x/iaGavTAyam9vh8lkQn19fcAwP9gXNBSCUqvVHtoCu92OsbExGAwGGAwGjIyM8C3TcOoV0U4rpEKYipnNZlRUVHgMkCUlJfGRRSifA8dxAb+PU1NTSE5ODvdtyAUCgAJIAWAJ9MCYJAemXThy5AiWLl0aUkvIV7fCZDLhyJEjHmmEN0JNK1wuFwwGAwoLC7F69eqwvuiREkHFx8djwYIFfOSQlpbGtwxtNpuHvkJKvULugmSk4W+ArLOzk9eZsMgiIWH2BjJf5xooMog1jYMfZGC6g+EXMUUO3tqFUNqRDMKLPFAaEeh5YsGq6MnJyfyymliD9+AUq1cYjUYMDAzA7Xbz9QoxRb1YjhyEx/SG9+dAKeXrNkePHvUYIMvIyPBJmsHSihiPHBjUAHSBHhAz5OBLu6DRaMImB5fLhSNHjkCj0fhMI7wh5a7N1I46nQ7Lly/HwMBASOcazjlIOaYQwnpFaWkpX9QbHR1FZ2fnLHmzt+9BpDFX0UiwATJW5GVkodFogpLD5ORkzEYOgq7EbgBpMz/z2fOPCXLwp12QolXwhkqlwtTUFLq7uwOmEd4Q+wVlpBMXF4f6+nrYbLaYHaYSc17eRT1Wr+jv74fZbPZYtAPIU5CMdFoRyjEDDZD19PSAEIK4uDj+u+mLJJgEPFZBCFkIIBNA18zfNZTSWVODc0oOwezb1Gp1SORAZ/ZPjo+PY82aNRH/RflqgZ5qZi+sXrFgwQI+TzcajTh+/DgsFgvcbjf0en3Elv/Gah3D1wBZd3c3TCYTPv/8c34DmTDCCkU6TQhZDODPAPIxXTD8I6X0ca/HEEzPQFyM6WLitZTSL0Qen1BKKSFkOYCrAPwQwE0AXgJwNyHkI0rpxx7vXdI7iCDEaBdCqTmwO7rD4cDixYsjTgxs56V37SLULsd8gDBPX7x4McxmM44fP46pqSn09fVFRIQUK5FDMGi1WqSkpPBaE+8I68SJE/jwww+RlpYmlfBcAH5EKf2CEJIK4HNCyDuU0lbBY74GYOnMn7MAPDXzXzFQAeAA/BeAvQD+D9N1BwBYAWAIwNyTgxwSaMCzG0EIgdUasI0rCW63G21tbXA6nT7Vjqda5BAIKpUK8fHxKC0t9ahXMBGSRqPhUxCxuoJYjRx8geM4/vfvHWHl5+fjo48+wr///W9UV1fjpz/9KS6//PKgx6SUDmH6AgWl1EwIaQNQAEBIDpsB/HlGxLSfEJJBCFk481yxSATQDuBcACMzP4sHMO79wKiSg1QJtNi0glKKvr4+DAwMoLq6GsnJydDpdBG7k1utVjQ2NgZUO4YqnpqP8L6QfdUrjEajhzGtUF/h75jzIXIApsnB1/AXIQQlJSVYsWIFvvSlL+H666+H3W6XfHxCSAmAWgCfef1TAYA+wd/7Z34mhhzYxfAZgCoAZ2KaYM7FNGH0eT8hauQgZmeEN9RqddAP1+l0oqWlBVqt1qMbEapewRtM7bhixQq+SOULUtMKnU4Hk8mE7OxspKWleXwesR45BLvLe49jT01NeUxY+vKadLvdEd/AFcrIvhiIEUEtXrwYKpUqlIXMKQB2AbiTUmoK70xnHVtNKf0NIeR/MC2CehiACcAdvmoXspODt3ZBigQ6WM1hYmICLS0t/Bi093ND7XSw8+7o6MD4+HhQtSMg/oJmKkqLxYLs7GwMDAzg6NGj/N2VdQNOFRBC+CnTxYsXw+12w2QywWg0etQrnE5nxCXewS7iUCGXCIoQosU0MbxCKX3dx0MGACwW/L1w5mdBMZOKcDP//zMAPxO87rWEkBZK6YjwObKSgz/7NrHwV3OgM2vOBgcH+TTC13NDjRzcbjcOHTqEjIwMrFmzRtR5iyEHh8OBxsZGZGVlobq6Gk6n0+PuajQacfToUVitVmg0GqSlpfG99XARyXw+nPqASqVCRkYGP3PicrkwNjaGEydOwGg0QqfT8SlIuHMQcvtH+kMo5DDTiXgWQBul9DE/D3sDwK2EkO2YLkROiK03kGlruFQAU5ieqZic+W8PgPsBfIiTNQgAMpIDKzoeOHAgJBdowPfdn8mq4+PjA4qaQiWH8fFxWCwWLF26lF/uIgbByGFiYgJHjhxBRUUFcnNzPc5NeHctKirCyMgIhoeHMT4+ju7ubl6QlJ2dHdIFE4tOUAwajQa5ubkwm81IS0tDSkqKxxwEc4TKysoSJW0WYq7IIUQXqHMx3WJsJoQ0zPzsXgBFAEApfRrAW5huY3ZgupV5nYTjvwrgKAAHgCQAcZh2h3IDKMF0euGBiJODt3YhVBdoYDY5sAusrKxsVhrhDankwKKRoaEhpKSk8MU1sQj0Hvv7+9HX14fa2lokJSUFPZZKpUJCQgLKy8sBzC7wpaSk8CnIXHhOyjl4lZCQMKteYTQacezYMb5ewUbSg6UhcpJDoOOGIoKilP4bfkanBY+hAG6RdOCTMGK6jtFEpm3iNJgmCAeATviYs4goOYixb5MCVnMQphE1NTWiPngp5MC0EayoeejQoYgUs1j7k5m9iD2edxTiXeCbnJz08JyMhoeDEHIUS31dyN4RFZM2M6IEMGd298HSiljZtC7ADZhpV85IqB0zf0AIeR3T6YYHIkYOlFLY7XbJvguBwFbTNTQ0BE0jvCFWQMWW4UbC8EUIm82GxsZG5Ofno7i4OGJ3WiKwUCsuLvbp4ZCdnY2srKyIujj7Oo9IQkw04i1tdjqdGB8fh16vx/Hjx/nFv+y9z2VaEWvkQCn1bosK/+0OXz+PGDkwQvD1Cw41DGUhZWVlJRYsmGWOGxBiBFRM7ei9xSpczYLRaERbWxuWL18eUvdBSivT28PBZrPBaDTyLs6pqalwu90RtZCLFcGSVqtFbm4ucnNzAZxc/MveOxuqCqVeEc65zpOpzKCIaFrh644bikW8MI1ITk6WTAz+zoXB7Xbz47m+1I6hRg6UUvT09ECn02H16tUR/UKKRUJCAhYtWsR7F5jNZnR2dvKfZ0ZGBrKzs32G4WIRq2Yv3r4N7e3tvI+n0+n00FeE0wEKdq5CBeV8huzvgBUVxZKDsBtRX1+PgwcPhvS6/i5wq9WKpqYm5OfnY/ny5T5/yaHMSXAcB5vNhqmpKdTX14cVzkZKBMXunJmZmfxUJds6xcJwVthMSkoSfXFGq+YQDgiZ3luZnZ2NnJwccBzH6ytOnDgBQgifgqSlpUl67UCfUyyL16RCdnLQaDRwuVyiwlrWjViyZAkfLYT6Yfv6BYpVO0pNK5hPoVqtxsqVK0M632iAtQ2FYTirVVgsFj4Ez8rKCtoJiMXIwRtCwmGms8J6xdjYGHQ6Hdrb23mizMzMDLtWE6ma21wjouTgb7IyWO7PTFOGh4dFt/ukgFKKzs5OjI2NYc2aNUHbf1LSCkY4q1atQktLi6Tz8vcFipZ8OjExEQUFBSgoKPBwhurv7wellCcK7ztrrKYV3ggUjXiv52NEKazVsMjC+/sS7HdzqkQPUUsr/MFb1BTp6rLD4UBTUxPS0tIkqR2DkQOlFF1dXTAajaIIZ64hphMgdIZid9bh4WG0t7d7LNOJlYJkJI8pJEpf7WKmr/Ceg/GGHE7fc4WopBX+yGF8fBwtLS0eaUQkwXEcDh48iKVLl/J3CDEIFjm4XC40NTUhKSkJq1evjviXOhYGr4R3VipYptPR0YHJyUlotVokJyeHXdxjkItwQtF9+GoXM31Fd3c3LBYLuru7fa4oDKVTQQh5DsDXAegppat8/Pv5APYA6J750euU0gclvzGJmJO0IhppRG9vL2w2G84991zJxw9EDkwX4WvYK1YRLtEQ4rlMZ2hoCBMTEzCZTHxxj0UVoe7zlCNyiNTglXCJjt1u582Eh4aGcOzYMSQkJPD1DLPZHIp0+gUAT2LaCcofPqKUfj3U9xAKopJWuFwn7emcTieam5uRmJgoKo1gIb7YXzLbS6HRaPgvtFT4u3PrdDp0dnbO0kUIIeUO6O+xsRA5BAIji6KiIgDw2OdpMpk8/CbFjizPl1SF4zi/9YpPP/0UP/rRjxAXF4eXX34ZF154oaiImFL64YyHQ0xB9vXDwshhfHwcBw4cQEFBAZYvXy7qFyelODg5OYkDBw4gJycHK1euDFmv4P081jPv7+9HfX29X2KQclGzidX5CO8LmS2RWb58Oc4880yUlZXxo+kHDhxAe3s7RkZGAr7f+UQO3sdk9Yr169fjT3/6E5YtW4a+vj7s2rUrki99NiGkkRDyNiEkKi0x2dMKjUYDu92O7u5u6HQ6yWkEI5dgee3Q0BC6u7s97ursIg/FgZiRAytopqeno66uLuAXWGwL1Ol0orGxEQ6HAwD4iUvWFYj1yCHQuREvv0lhvt7T0wOVSuWRgrDPc76QQ7A6ht1uR3FxMe65555IvuwXAIoppZOEkIsB/BXTPpKyIipmL/39/cjJyQmpGxGs28HUjna7fZbaMdyluCaTCc3NzaILmiwFCrbToKmpCWVlZcjMzATHcR799sTERKSkpIRlVBMNSDHsEZrYeBuyJicnIzs7O6xdmf4gV+QQ7Pcbaem00BGKUvoWIeQPhJAcb3OWSENWchgfH0dHRwdSU1OxfPnykI4RiByCqR3DSSvY5J/YKVAgeFqh1+vR0dGByspKpKSkwOFweAiT6IwF/PDwMEwmEw4ePOghd46V5azh3OW9DVlZy9Bms/HvN5ITppEmHDmMXoKBELIAgI5SSgkhZ2K6HDAa0RfxAVnIQThjcMYZZ8BoNIZ8LH8X+MjICI4dOxZQ7RiKtb3b7cbw8DCcTifOOussSW06f+TANBFMhBUXF+fzcSwkX7RoESwWC1asWOHh6hwXF8dPXEa6wyMFkTKDFbYM9Xo96urqZk2YZmVlITs7W9YJUymQyQVqG4DzAeQQQvoxbeGmBXiTl8sAfJ8Q4sK0e9NWKjHvJIQ8CMBIKf3tzN8fwnTr9HF/z4l4zcHhcKC5uRnJyck488wzMTU1Bb1eH/IxvSMHKWpHqdb2drsdjY2NiI+PR25uruT+va+aA/OKSEhIQF1dnaSLynviUqg1sNlsSE9PR3Z2NjIzM2MmqggHvraEe6sWWYoSySlTKRBjLitFUwMAlNIrg/z7k5hudYaD5wC8DuC3M2YvWzHtQO0XESUHu90+S3QUrtGr8PmMeFJTU0WJj6SkFWwZ7rJly+B0OmGxBNxO7hPeyko2c1FUVISCggJJx/F1Y/CWO09MTGB0dBQ9PT38hSW3jwMgT/HQF7xNbpi8+8iRI+A4jpc3RzPlCuaSHavj2pTSHkLIKCGkFtNbtQ5TSgOmJhElB7Y3UsjqgRSSYsDIgQ1lSVE7iiEHOrPzYnBwEHV1dUhMTAx554Xwoh4dHcXRo0exatUqpKenSz5WMHgbn3jfZdPS0viogp1bpBAtchCCTZimpaWhpKTEY5FOR0dHyBOmUuFvZwWDHDWHCOIZANcCWIDpSCIgIu7n4B3uhRs5EEJgMBgwNTUluQ0ajBw4jkNr6/RCofr6+rB3XrDIgak/Q525CKWVKbzLsiGq0dFR9PX1wWaz8evkvTdmh4JYaLN6L9LxNWHqdDojbnkvk7lstLAbwIOYrmd8O9iDZW9lhrMp2+VyQafTQaVSeVy8YhGs09HY2IhFixZh8eLFHhdMOE5Qx48f5yMoOSzKxEA4RMXOCQDfPkxJSeFTkFBy97mIHILBX8rV1NTE78ZgKUg4vxcxNYdYTCsAgFLqIIS8B2CcUhr0oow4OXjf9UL9Ek1OTvL1hZSUlJBySn8RAAv5/XU6QjF7sdlsGB0dRUFBASoqKsL2A4jk3VmtViMtLQ05OTl8+3B0dBRHjhyB2+32O5od7BwjhUhHIiqVCmlpaUhMTMTq1atneU0KJ0wTExMlvZe5aGVGCjOFyLUAgi/vxBxu2Q6E4eFhdHV1obKyElNTU5iammWMKwq+ZNA9PT0wGAwBQ36paQWbLk1PT0d+fn7M3VWFELYPWe4uHM1OTEzkowp/NneRjhzkVkd6e01aLBaPrg8zuRGu5wt03GAiqFgzlwUAQsgKAG8C2E0pPS7mOTFFDm63G8eOHYPNZkN9fT20Wi2sVmvITtDCi5y1FOPi4rBmzZqAd0gpacXAwAB6e3tRW1uLnp6eiNm7RSuv9yXCYpGV0+nkw/GMjAz+M5sv5ODvIhZOmArX8/X29vL2cdnZ2T4nTOej8zQAUEpbAZRJeY7saQVDsC8As3LPy8vDsmXL+MeGU9BkIqipqSk0NTWJbimKSSsYkdntdtTX10Oj0cT8TEQwCOciioqKeNt7YUcgOzub33saKchBDmLHtYXr+crKyuB0OmE0Gv1OmIohh7kUqEUSUYkcgjlQszuVLyv3cMhBpVJhbGwMfX19qKysRFpamqTz9QfhzkshkYVrac8QKyTjLUpi4bjRaMTo6CiMRiMfjoejM5hrFyghtFot8vPzkZ+fz0dSRqOR3xDOWsbZ2dk+9Q5ybfaeC0SFHJjWwftDE2O1Fio5UEqh0+lgNpuxdu1aSVX5QBe52WxGc3MzysvLZ+ktQilk+kKs1ixYOM5ao2wGhe3zZEQiVWcQqxOZ3hOmbrcb+/fvh9lsRl9f36wJ01Bw/fXX4/nnn9fDvwsUAfA4pndkWgBcSyn9Ipz3JRaypBXeYIYvwguUqR1TUlICqh1D0Rw4nU40NTVBpVJhwYIFktt1/l6Tmb1UVVX5rEhH8o4fC5GDP7DZCu9py9HRUd5GTehkHUyGHkuRQyCoVCpotVp+h6nD4eANeU0mE37/+98DmK5DFRYWijrmtddei+eff34D/LtAfQ3T49lLMb1Z+6mZ/8qOqEQO/hbiilE7So0c2J29rKwMGo0Go6PSh9e8IwBKKTo6OmAymfhCqb/nSSlkdnZ2IjU1lb/jMhKL1chBCO9zjI+P55fpCEVYvb29PJH42xIeq5FDMMTFxfETpuz1br75Ztxwww1IT0/HX/7yl6DHWLduHTC95NYfNgP488yg1X5CSAYhZCGldCgy78I/oppWMKnywMCAaLWjFHJghi/szj42Nha2E5TL5UJzczOSkpKCmr2IIQfmKmWxWFBfXw+r1cprDpgdfGZmZkTSE+9zixSCXcxCEVZZWRl/h2Vbwr0HqOaLlwPgP6JTqVSoqqpCdnY2/v73v0fy91cAoE/w9/6Zn80/cvCXVrA0QqVSSV6IG4wcmCWZxWLxMHwJVZ3JyIENThUXF/NLdsU8zx+Ya3Vqaiqqq6vhdDo9NAfCSvnU1BRaW1t5zUE4EuBIpyhS7/TCO6z3AJXb7UZKSgpcLldEL+hoRA7eEKoj50odG0lEJXJg+wrLyspE52IMwe7Gws5BbW3tLBl0qDMSTqcThw8fljQ4FehcGdGUlJTwU4beYJXyvLw8TE5OorCwEKOjo+jv7wcAPv3wFZpHG6G+vvcAldPpxODgIMbHx3Hw4EEkJyfzKUg4u0Dmoo4hhwsUgAEAiwV/L5z5meyQnRyGh4eh0+lQVFQkmRiAwF9CVruoqKjg1W9ChEIOzDbfbrdj3bp1kr6g/siBjYOvXLkSGRkZoo4jvIhKS0tnheZs6lJMwS/SiGQkotVqkZ6eDrvdjqVLl/Kb1dkyGSZIkjoTMRcWcTJJp98AcCshZDumC5ET0ag3ADKmFUwkZLVaUVJSEvEvcH9/P/r6+gLWLqSSA8dxaGlpgVqtRlJSkuQ7l69WZn9/P/r7+31u3WZuzMJpUH/wDs1NJpNHwY9FFdFwTJJLIUkIQUpKClJSUngR1tjYGAwGg8dMRHZ2dlDL+7kwlw1l6OrKK68EgE/h3wXqLUy3MTsw3cq8LpRzDwWy3HKY2jE3NxfLli3D0NAQ7HZ7RI7tdrvR1tYGl8sVtHYhpZhps9nQ0NCARYsWoaioCJ988onkcxPWOCilHlJw4XkyW3qNRgO3280XawkhvPt0IBBCPAp+Qi8H1kZkXg5yRBWRJgd/F7JarfYYy/YWJDF/TV9+k8FMWUKBmIlMqZHDtm3bsG3bNr/bkWa6FLdIOmiEEPFvjtlsxueff+6hdgzX04GBkU5+fj6Ki4uDfkHFRg4s7A+2fTsYWFrhcrnQ2NiI9PR0VFdXe9ivu91uvjrPvrxut9vj32w2G9xuN1wuF1QqVdA7oLeXA4sqenp6+NdISEiI2EU9VwVO4UwEMwAaHR1FV1cXtFotH1UkJSUFTQFCwVw4T88lIk4OKSkps9SOkSCHQBJrfxBDDizsZy5Q4YAQwlvllZaWemw78iYG78IpMP05TU5OorW1FUuWLAEw/YVknx2LKoINjbFZgSVLlsBms+Ho0aMYHh7GwMAA7+4cruR5rgevvC3vbTYbb/ZitVr540YyeprP49qhIOLkoFKpZuXq4VjFUUrhdDpx/Phxn3l7IATqHrB9F06nMyQjGV+Ympri7eaEHY5AxCDE6Ogojh8/jlWrVvFfMmFUwT5DjuN4kggWVSQkJCA1NRUZGRnIzMzE+Pg4r2RkbtZicnghopVWSEFCQoKHCKu1tRVWqxUNDQ1BRVhiMc9doCQjagrJUFa/sQKh2+1GTU2NJGIA/N/dHA4HGhoakJub63PfRSjo6+vD8PAwFi1aNIsYhDUFf6/Fnl9XV+ch9xZGFVqt1qNGwXGcR1FTTFQhvNsyARbL4VlnQDie7QuxPrLNZM4LFixAenq6TxFWKPoRMa1MMXqY+YI5kU+LAdMFLF68GC6XK2J5Ltti5a/9KRXC0e2lS5fCbDbz/8bqD4D/TgRTTDocDtTV1QWNYITRAosqOI7j0yd/UYWviy8xMRGFhYV8Di8cz05ISOCjCqmkLBVyaxJ8ibCE+hHhOsJAJCUmclBqDgHg68OVmlYYDAa0t7fzAqTR0dGIFDSZw5SYLVYsJQn0ZWEDXunp6Vi2bBlGRkZAKRWdRjBpdlpaWkjWcsKoAgAfVbD/ss+MnUsg+BrP9jZ9YXqDWI8cAP+E460fcTqd/Ibwo0eP8uv5srKyZqXHYgqSSloRBN65vti0wt8Idyibq7yPe/z4cUxOTs7ap+kPwTwoLBYLGhoaUFZWxhcemc5BDDHYbDY0NTVh8eLFWLjQbydLEnxFFXa7HRMTE8jJyYHT6RRV1AROdgbYMtyxsTHeg9HhcECv1yM/Pz8sFSNDpDZoCSE2GtFqtcjLy0NeXh4opZiamsLo6KhPERbHcQHTEKXmEALEdA3YXZSZggp/seE4WFNK8cUXXyAtLW2WvDoQ2IXuixz87aQghPB37kDEMDExgdbWVixfvlyUYjIUqFQqWK1WNDc3o7S0FFlZWfy5CYubYgRYQr0B+zxdLhdaW1v55TLCLeFSIYcmIZRURSjCKi4u5ndjMFJ0u918J8hXATeWnadDQVTIIdgFyTZPl5aW+ryLhtoKZea0ZWVlku/O/giNLcDx7pxQSpGYmAir1YoDBw4gKysLOTk5s4p7Op0OPT09qKmpCbt1GggmkwktLS1YsWIFT2DCqMKbJJjAJ1hUwf69sLAQZWVlHga1x44d48Ny4Qh6MMTqyLZwNwYTtQHga0TeS3/Zyj4p2LdvH772ta8dA6AG8Ayl9BfCfyeEXAvgEZycp3iSUvpMWG9MJKKSVgQCM1CprKz0+8GGklawukVycjLy8/MlPReY7QYlbH2uWbPGp+JRq9WitrYWHMfBaDRCp9PxF0xOTg4sFgtMJhPq6uoiumjFG6yoWF1d7VNa7iv98CaMYK1SdjF7G9ROTU1hZGTEw/Y+Jycn4DKd+WD2wj6PvLw8ZGRk8AVcpq147bXXYDAYMDw8jPLyclFkx3EcbrnlFmDa0KUfwEFCyBszZrBCvEYpvTVib0Yk5sx9mtUBzGZzQAMVQFrkQClFd3c3RkdHUV9fj4aGBtFmo0II5yScTicaGxuRmZk5q/XpK41Qq9UeF4zJZEJbWxvsdjuSkpLQ39+P3NxcWeYg+vv7fbZE/UFsUVMYVfi70wvDcjZxOTY2NmuZTnZ2tsfvO1YjB28IC5LeBdyMjAzceOONeOSRR3DnnXfik08+CVqPOXDgAMrLy9HZ2dkFADPDVZsBeJPDnCCq5MC+BA6Hg6/yBzNQAcTXHDiOw5EjR6DVavm6Rahj2+x5U1NTaGxsxJIlSzwiELEdCafTiY6ODn6zltPpxMjICLq6ujA1NYXMzEzk5OSErVikM9vH2drAUI8lplXKfh4M3sU+1kJsamoCcHIEXQ6zFyDyjlqBbjJLly4FpRS7d+8GAFGf/8DAABYvFk5jox++LeC+SQhZB6AdwA8ppX0+HhNxRI0chDbxzc3NkhbiiokcmBpu8eLFHqPh4ZDD2NgYent7ZzlXCxWLgQal2HtdsmQJr6mIi4vzUPKNj4/DYDDw2gKW40rRFjBFoFarRVVVVcQuCl9RxdDQ9LQwU65SSkUJsAKNoBuNRkxOToLjuDkZQReLYFOZrOsSYVL6G4BtlFI7IeRmAC8CuCCSL+APstUcvKFWqzEwMICBgQFROgPv5wbalWA0GtHW1uZzcCrUTsfU1BRMJtOsORGh4jHQF8FoNOLYsWNYtWqV31qKt2KR5estLS38hZKbmxtQnMO0Frm5uSgqKpL8PqWgv78fBoOBj8q8owrhoFiwkF4oTDp27BhSUlIwOTk5JyPoYhFI5xBKalRQUIC+Po8gYJaRC6VUaIL6DIBfSXqRMBAVimYRg8Fg4BfASEGgyKG3t9dn90D4XCnqSlZ4dDgcWLlypV9iCJRKDAwM8DMWUnQAzAa9uLiYt4zr6+uD2Wzmd10K9yWwKdXS0lLRUVgoYCmLxWJBbW2tx4UvpVYR6PipqakoKCjgPSfZVOnU1BTS09NlHUEXi2AiKKkEUV9fj+PHj4MQUoppUtgKr+3XXmaymwC0ST7xECH7J22329HY2AitVoulS5eG9Mv1dfdnvg4cxwUcnJISOTidTjQ0NCA7O3uWtFpYeAwkhe7o6IDFYhElhQ4E7+UqExMTGBkZ4cewU1NTMTIyghUrVsimlQBOfs5qtRqVlZV+v/xiahX+iML7ooqLi/M5gn7ixAmPQqDU/RjhItDFH4q8X6PR4Mknn8Qll1zyd0y3Mp+jlLYQQh4EcIhS+gaA2wkhmwC4MO1SfW2o5y/5/OQ4KPsAmU/CsmXLoNPpQlY5ekcOjHDy8vKC+jqIrTkwrQUrPLa3t0uSQrNiaHJyckTzfmD682Tim/LycgwNDeH48eNISkrC0aNH/WoqwgXHcWhubkZ6ejpKSkpEvydhrcK7kOktwArU/WDHYu8dmD2aLYwqGBnLufMjkOI1FN3KxRdfDEpphfBnlNL7Bf9/D4B7JB84ApCFHCilfLjPfBJGRkZCmswEPL9kzDfyjDPO4B2CAkEMOYyMjODYsWOoqqri6wNCtaNYKXRhYaHsU3ksZWFbvHxpKnJzcyWJkHyBtW8XLlwoar+oPzCi8BZgsdSDTZYyEhYzgi4s6DLDl+7ubmi1Wn4uItp1ilNNHQnIRA7Dw8OYmJjwCPfD3XnJcRwGBwdx4sQJ0Tsv2HMDkcOJEycwPDw8q/CoUql4u/RAxGA2m3HkyBEsW7YsLBepYGBzJ5OTkx4pi7emYnJyEiMjI2hsbAQAvvshxcdAzloGSz+YRd7g4CC/DY2RhVivCpVKhczMTP5zZyPox48fx9TUFI4fPy5qBD0SONVcoACZyGHhwoXIycnx+DKGSw4TExNwu92SC5r+Xte7ZiH88rACWXt7O3Q6HXJzc5GTkzPrTmwwGPj1eHJ+Mdi5ajSagCkLIYTfg8HahSMjI+ju7sbU1BQyMjKQk5ODrKwsv/UQtpFcbrIDprsfo6OjPNl5e1Ww/xfTKgVOjqDn5eWhpaUFWVlZs0bQ2bZsqQiWqiiRg0j4utNqNJqQ0gqn04mWlhZQSlFTUxPSWLN35MB2XeTk5MzKpVkqkZWVhbVr1/JdFnYnZkQxMjKC0dFRrF69WlYpNFuEk52djeLiYknP9aWpGBkZQWdnJ+Lj4/n3wro8bCCssrJS1ulCpmKdnJxEdXW1R+oRileFN5gewdcI+rFjx2aNoIuJKoIJtU41izggyiIoqQ7UwoGs3t7ekPJIb3Jgx/Teku2v8MjkwOxOrNfrcfjwYbhcLixcuBBms1m2sJXVMoqLi0OaDxEikKbC5XIhMTERZrNZUsoWCphs3ul0Bu1+ALNbpWKiCl+1i0Aj6ImJiTyR+Gs9BxNAKeQgEv5EUFLSCr1ej46ODv4u1t3dHdK5CFuZbBhLWHgExEuhCSHQ6/W8CnN8fJwvBKakpPCFwEhEEqyWIddYt1BT0dfXh97eXqSmpqKxsdGnpiISoJTybdEVK1ZIIvtgUYVQgBWssOk9gs6iCubhwPwmhSPoYsxllbQi1BcS6QbFCm9jY2NYs2ZNWBV34OSezp6eHuj1etTX13scU6ywyWKx8N4ILOIQfsEmJyeh1+vR29vLf/lyc3NDuguzoprctQwA6OnpwdjYGNauXcsLxkwmEwwGA68rCOe9MLjdbhw5coSPwsLpJgQbFrPZbAA80xB/IITwRFlUVOR3BD0xMfG0coECopxWBKs5CA1f6urqIhKqE0IwPDyM1NRUrFmzZlbhUQwxjI+Po62tDStXrvSYsRC+BisEMjt41h612+28qCo9PT3oRTE4OIiBgQHRU5WhgoX3DofDI+8n5OTCHACz3ktmZiZyc3MlpVIcx/F1Ezkk3sKowmq1oru7G0uWLOEjDEYQYoqavkbQR0dH0dfXB7vdjq6uLp9+k6GkFfv27cMdd9yB9vb2Dvj2cogH8GcAqwGMAriCUtoj6UXCQFTJIVDkwAxli4qKwuqrC+FwONDZ2QmtVotVq1Z5/DLFmL8CwNDQEL92T+wwVEJCgodx6+joKAYGBtDW1oa0tDTk5ubOGjBiRTrm9xDphSxCsEGtuLg4rFy5MiBheb+XsbGxWT4Vvjo5DGzBz8KFC2XXgLAW7LJly5CRkeHTq8LlcoEQwguwAkE4gp6eno7h4WGkpKTwfpNsBD0rKwsWi0X0PhXgpJfDO++8gyVLlqyAby+HGwCMUUrLCSFbAfwSwBUhfDQhISZqDv5s18IBKzwuWLAATqeTPyex9QU2T8B0BaHm3mq12mNsmYXsTLTD6hTd3d1Qq9UeG7LkALuLZ2ZmoqSkRNJzvXP1YJoKtgIgEgXVYLBarWhsbMTy5ctnOV8J0w9vARZ7XLCogvlHCn+X7P3/8pe/xJ49e1BXV4fKykrU1NQEJR7m5VBWVgZKqcOPl8NmAA/M/P9OAE8SQgiVUwIqgGyRg7cblK9WJlNS+hIh+YLYwRZh4dHlcnmMGYuVQre0tCAhISGiF6swZC8vL4fVaoVOp8OBAwegVqv57kcg16RwwFq4hYWFYZva+tJUMKUis0sbHx9HRUWFrENhwHTU2dTUhBUrVvhM+xi8BVjC0XsgcKvUuyApfP8//elPMTExgczMTPzmN7/B3Xffjerq6oDnLNLLoQBAHwBQSl2EkAkA2QBGAn8ikcGcpRVutxstLS0ghMwSIflCMA0+MH3xnzhxwqPwaDKZZoWXgYjBbrejqakJCxcu9PCFkAOs+7Fs2TJkZ2fzW7PNZjPS09P59CMSKQa7s5aXl4uSnUuFcFhqcnISDQ0NSE9PR1dXFwYGBmZpKiIF5pmxcuVKSf6NwqjC17Ig71ZpsG4Fx3HYuHEjzj333LDfU6wgauQg1Buw3HDBggUoKioSdZdkLUl/JMLyaAAehUf2PDGFx8nJSRw5cgRLly7lxTNygb3WGWecwSsRmb8BmxlgCsyEhAT+4grFCn5ychLNzc0eZrNygb1WVVUVfxe3WCwemors7Gzk5OSIKtCKea1IiLaCtUodDkfANqnUVqYYL4eZvy8G0E8I0QBIx3RhMiqIWlrBvgTj4+NoaWnh75ZiwSIPXxoCltv6mtJUq9WYmJjgZdD+yIXJbCsrK2VvH7I18v5ey3tmgKk0m5ubQSnlW4tijFDGxsb4oTK53xdTWHq/VlJSEoqKivhWobBAm5qaykdIUvQhTAcix/vyrlWYTCYMDw9jxYoVfmsVUp2nmZdDd3c3ysrK4uDDywHAGwCuAfApgMsAvButegMQZQ9Jh8OBtra2kFR4/gqaZrOZt50TejCwNEKj0WDVqlUwGAzo6enhZcO5ubn8Xbivrw86nU729iEw3f3o7+9HbW2t6CiA9eFLSkr4eQmhB2Vubi4yMzNnEZ9er0d3d3dIe0alghFeMMt9jUbj4VPhS1ORk5MT8II3mUxobW31664dSUxNTaGlpQVVVVVISUnxK8AyGo2Sohfm5XDRRRcB0wYuvrwcngXwEiGkA9NeDlsj/f4CgQQhopBZyul08h+e2z29T3JgYADnn39+SJV/JqMWsjNTUbJfHH/SAeoLFosFer0eIyMj/PnFxcWhqqpK1sk91qqcmJhAVVVVROoIbrcbY2NjMBgMGBsb48e1c3JyoNfrMTQ0hOrqallnP4DpAjBbMxjOBiymqRgZGYHNZvOpqRgfH8fRo0dRXV0t694P4OQQWqC0xe1245NPPsE111zD2/WFgNjwwfOCbOTgcrnAcRxfIc/KysLw8DDOOeeckPLMlpYWFBQUICMjA5RS9PT0YGRkBNXV1bMUj2IKj6z/zi5Sm80mSawkBcx6jhCCM844QxYSYq01g8GAgYEBcByH4uJi5OXlyZpOsEiopqYmoiTENBUjIyM88SUmJsJgMKCurk72SIgRQyAPUAA4dOgQbr/9duzZs0fyYJwAMUkOsqYVLORnQ07sbh3KXZOlFazLoVKpZq3NE6t4tFqt/EAT23PpLVaKVLeAqT4zMzODulaFAybYGRgYQFZWFsrKyngZts1m481qxU4hikF/fz90Oh1qa2sj7u3orakYHBzkp0mbm5t5IpfiUyEWrDUajBgaGhpw22234fXXXw+HGGIWskUOAwMDvNaAhWSHDh1CZWVlSKEns0Xr7+/32eUQq3hkRbNAlXtKKW8ZbzQakZiY6NfTIRCYnV0kl+X6A5tdSEpKwpIlS2aNoRuNRhgMBkxMTPBFwHAGq3p6ejA+Po7KykpZ1ZzAtEiuo6ODT1uYpsJgMPAGtJFq+zJiCNYaPXLkCG688Ubs3LkTFRUVfh8nEjEZOchGDqOjo4iPj/cINQ8fPowzzjgjpCJSa2srDAYDVq5cOatPL2Z5LTDtUHXixAlUVVWJzleZvt5gMGBkZASEEL6gGeh9sFZlRUWFJFltKGCeDzk5OUFnF4RFwNHRUWi1Wr77IeYzYcpRm82GFStWyO6wxNSkNTU1PolZ6FNhNBoRHx/Pvx+pqQfTggQTU7W1teG6667D9u3bsWLFCsnvyQdOL3Jg3oBC+CoqioFer0dLSwuKi4tRVlZ28uQkSKFZMbCysjKsEJgVzfR6PZxOJx/eCgdxWOV+1apVsk/qsTZuUVERnyJJgdVqxcjICAwGg9/3w0BnlslSSrFs2TLZfRr1ej1OnDghqZ7BNBUGg0GSpsKX/NoX2tvbcfXVV+OVV15BZWWl5PfkBwo5tLa2YuHChaLtx9hFPTo6ipycHKjVav7OKJYYmDhKo9GgoqIionc61rM3GAy86YtarcbY2FjYlXsxYCFwpERb3u9HGK4TQtDa2oqEhIRZaYsc0Ol06O3tDavQyVqMBoMBJpMJqampvE+F8Jg2mw0NDQ1BiaG7uxtXXnklXnzxRdTW1oZ0Tn6gkMOxY8d4Jg8GlkNrNBosW7YMw8PDsNvtKC0tFV14ZDs58/LyZN8GxXEcjh49CqPRCLVajZSUFOTl5UXM/MUbTATkb4w8XLBdGSz9sNvtyMjIwLJly2QnvaGhIX4zWqQKnSydYvZ+rODJvELZJKc/9Pb24oorrsAzzzyD+vr6iJyTADFJDrIqJL0h1g3KbrejoaEBCxcu5C9qJl31XkfnD0xzL9csgRButxvt7e0ghODcc88FIcTD/EXoERCJFhxLW+QUAbFdGSkpKTCbzcjLy4NarcaRI0fAcRyf10e6WzA4OIihoaGwlgH7gnDojXluDA8Po7GxEXFxcdDpdOA4zqeYbGBgAFu3bsVTTz0lBzHELKKqkBTjBmUymdDc3DxrLwUzi2Hz+IGIgbXwopHzsxHojIwMD7NaofmL1WqFwWDgd2CGc2HpdDrenl/uOzjbACac4mSr+oSu1mxTeFZWVlhpW39/P/R6PWpqamTvgDAToJqaGqSlpfFisvb2diQlJfEFZ5vNhiuuuAKPP/44zjnnHFnPKdYQVXII5gal0+nQ2dk5a9EupRRarRY6nY6X3/qrrPf39/N3HrkvHrGtysTERH62wPvCYvoDMc5KfX190Ov1qK2tlV31yKK3srKyWao/rVbrsa6OtX1ZuzmUtm9vby9GR0dRXV0tOzHY7Xa+c8bqX0JNBetOXX/99Whvb8fXvvY12QfWYhGy1RwopXA4HB4/Gxwc5OsG3o9lvpHecl/hzD27sFinICcnB3l5eXx0wAQ/K1eulP0LxtKWcFqVbrebL5iNj4/71R8IF9qsWrVK9vfGKvdS35t32xeAR9vXX5TU09PDd5Lkbo06HA4cPnwYS5cuDfjeRkZG8M1vfhN33303vxbw4Ycfluu0YrLmEFVy0Ov1MJlMKC8v53/GjFVY4VGs4lFIFBaLBRzHISMjIyq9dzbpGMm0xVt/EBcXx9+Bu7q6QAiJSvuQkV6wyr0YOBwOniisVqvPWYnu7m6YzWasWrUqasRQXl4esLszNjaGSy+9FPfddx82btwo6znNQCEH1iZbtmwZAN+FR+HzxXQkmDdEeno6OI6DyWRCRkYG8vLyfBaXwsXw8DB6e3tRVVUlq77fYrHw9QW1Ws1vcpJzToJ1QORYasNmJViUxCYcCSExRQwTExN8xHDppZfKek4CnF7kAGDWEpvx8XEMDAxg5cqVfOHRl6+DWMWjyWTivSFY7shyYL1ej7GxMb6lyHQSoYK5TBmNRlRVVUV8lsAbTqcTTU1NyM/P5+dS9Hq9bANibNqxqqpK9jFopj2ZnJyESqXi933m5OTI8tpMKFZWVhawc2U2m3HZZZfh1ltvxRVXRM3HFVDIYfrD7+7uRl5eHrq6ulBdXT2r8ChG2ARMpyhdXV0Bv8zCUH1kZAQJCQk8UUgpljFlIMdxWL58uex3OVbo9GXMygbEmLAnLS0NeXl5Yc0VsNmF6upq2acdKaVob2+H2+3m0ySmOjUYDHA4HB5DYuGSn9PpxOHDh4MSw9TUFL71rW/h+uuvx1VXXRXWa4aA048cHA6HhxuUxWLB559/jsTERL+FRzFSaHYHr6yslFS1n5qagl6vh8Fg4O9WeXl5AS8IVoxKS0sLexmL2HMUW+iMxIAYS138zS5EEoxkAeCMM87w+Vn6Ij9fVv5iwIihtLQ0oM+C1WrFFVdcgSuvvBI33HCDtDcVGZze5MD0ABMTE1i3bl1Io9Zs2zQrzoVzB7fZbDAYDNDr9eA4jq+qC3Nt5kVRUFAg+84F4OTEaLBRYV8IZUBscHAQg4ODUTGEYavwNBoNli5dKopkhSpNo9HIW/mLEZMxjUZJSUlAYrDZbPj2t7+NLVu24Oabb5ad/P3g9CUHoaHs0NAQ1q5de/IFRBIDy8HZ5GEkf4lOp5MnCpbTp6amoru7GxUVFbKbzQInhVuRcjiy2+0wGAx8qO49UMV0BZFypQoESilaW1sRHx8f1lwGE5MZDAZwHMe/J28rf5fLhcOHD/NmN/7gcDhw1VVXYf369bjtttvmihiA05EcnE4nxsbG+IWwmZmZ+PTTT3mlmXBnQKAogA0YlZWVyb4DgeM49Pb2oqenB3FxccjOzkZeXp5sm7SBk1u15ArtvQeqmNOyHCYt3mDmPMnJyR4TteHC6XTy72lychIZGRnIyclBWloampqaUFRUFHCRjtPpxHXXXYdzzjkHP/rRj+aSGIDTkRx6e3t5xSMLbT/55BOcffbZoguPY2NjOHr0qGwDRt5gOXhVVRXi4uL4Ve1MpMSGqSJ1t+3t7cXIyEhUOiCsGGixWJCYmDjLdzLSqQUbnmOLb+SCsEM1ODiI5ORkFBYWIjc31yfZulwu3Hjjjaiursa9994718QAnI7k0NPTg9zcXI8v/SeffIKzzjpLFDEMDg6iv79fdk0BcHL7Fgu1vS9U7ynFxMREvvMRykVFKUVHRwev6JS7A8JyfrVajYqKCn51ABsQY5OKeXl5ERkQc7vd/Nq9aFiouVwuNDQ0oKCgAGlpaXz6AcDDyt/tduP73/8+lixZggceeCAWiAE4HcmBmczyB6MUn3zyCYqLi2eRhseLzrgNTU1NRUUuLLVVyYp/zMWaTV3m5eWJmudghVXmMSH3F5TdwVNSUgJ2XHzl9EyeLuUcWfE5JyfHe+WbLOA4Dg0NDVi0aNGsGRdm5W8wGHD//ffDZDKhtLQUf/7zn2XvzkjA6U0OrPA4OTmJ4eFhv3dfJqdOTExEeXm57BcOx3H8hVNWVhbS67GLSq/Xg1LKV9R9qRn9TXHKBfZ62dnZkjwtmDyd+TSKHRDjOA6NjY3Iy8uTfZ0gez1/xCCE2+3GD3/4Q5hMJuTm5uLw4cP48MMPI37j4TgOa9asQUFBAd58802xTzt9ycGX+Su7++p0OoyMjCAuLo63ry8sLERBQUE4Ly0KrFW5aNGiiL0emyfQ6/Ww2+38cFhqaqrHOvpovD8Wai9atCisVqzYATHh+4tG65cR0YIFCwK+ntvtxn//938DAH73u9/JmsI99thjOHToEEwmk0IOgeByueBwOETVFwwGA1pbW6HVahEfH8+vOpdr7JqJjeTci8m6BHq9HmazGU6nE0VFRVGJGJhk2JfKMhx4q07ZBrHMzEy0tbWhoKBAdqdt4CQx5OfnByRat9uNn/3sZzCbzXj66adlJYb+/n5cc801+MlPfoLHHnts3pODrOXxV155BWVlZaipqQn4S2ELY9esWYPk5GQ+TG9ubgYAPp+P1Iaj8fFxtLW1hSQ2kgLmPZGcnIympiaUlJRgamoK+/fvj4js2R+YriTYkFEoEDoqlZeXw2KxYHh4GAcOHEBcXBxsNpvkpbJSIZYYKKV46KGHMDo6imeffVb2ou+dd96JX/3qVzCbzRE5HiHkGgBvUkpHCSEkmnsyAZnJIS0tDb/73e9w7NgxXHDBBdi8eTPq6+s9Uou+vj5+ixErEAnNUex2O/R6PVpbW3klY35+fsgDOmx/ZG1trewdEOAkEQn3d7DOB1vnl5yczNdewm1nMk1IME/ESEGj0WBkZAQrV65ERkYGRkZGPBbp5OXlRXRAjHVB8vLyghLDr371K/T29uLFF1+Uvaj95ptvIi8vD6tXr8b7778f9vEIIX8CcAOA7wF4NtrEAMicVjBYrVbs27cPO3fuRGNjI8477zxccsklePPNN7F161bU1dWJYnVhPu9wOJCTk8PfmcV8+Xp7e2EwGFBVVSW7XBg4uUMy0ECTsJ04MjICrVbLtxOlplRsJb3cERFDoGlHOQbE3G43Ghsbg3ZBKKV4/PHHcfjwYbz66qtR+V3fc889eOmll6DRaGCz2WAymXDppZfi5ZdfFvN0jy8vIeT3AOIA7ABwN4BbKaXtkT/rICcVDXIQwm63469//Svuvvtu5OXloba2FpdeeinOPfdcSb9EodmL1WpFdnY28vPzZ0lpgZPiH6fTGRUzGGBaozEwMDBrl2cwsEW/rEfPai/BUqqJiQm0tbWhsrJS1pCegXlxiEldhJGS0WgMaTqWRQzZ2dlBieGpp57Cv//9b/zlL3+Zk3bl+++/j0cffTSkmgMhZCuAqwFcAkAL4P8AvE4p/RshRE0pDe7QHCFEnRwA4P7770d1dTU2btyI9957D7t27cLHH3+MM888E1u2bMF5550n6ZfKcRxPFJOTk8jKykJ+fj7S09P5Hn9ycnJU9i0A0+KvsbGxsOcW2HyEL1s84fsQulHLvXkaOLnnIRSLvFAGxNxuN79vNFA7llKKZ599Fv/4xz+wa9cu2T1E/SFMcigAMMjSCELI7QBuAVBPKTVF/GQDndRckIMvuFwufPTRR9ixYwc++OAD1NbWYsuWLbjgggsk1QbYXki9Xo+JiQm4XC7k5+dj6dKlUVEhyhWheOsOmEDJ4XCgp6cH1dXVUbkYmL9kpGoawQbEGDFkZGQEVVq+8MIL2LNnD/bs2ROVelIEMeuOJYwSCCF/BNAI4A8AEK36Q8yQgxAcx+GTTz7Bzp078e6772LFihXYsmUL1q9fL7oQabFY+Iq2zWbDxMQE0tPT+bw30kTBBowSEhJkF28xAuzp6YHJZEJ+fj4WLFggy/sSghU7I+Ev6Qu+Nm5NTU0hJycn6GzGK6+8gu3bt+Nvf/ub7E5WMkAYOagope6Z/9dQSl2EkJsArKOUfjeqJxWL5CCE2+3GwYMHsWPHDrzzzjsoLy/Hpk2bsGHDBr9FN+aLIBzWYsYoOp0OY2NjER2iYotss7Ozo7aKne14WLVqFSYnJ3nPA+GmrUgOck1NTfHbp6MxAMcEXGxyN9CA2I4dO/Dcc89h7969su8pkQkEOBktEEISALgopa6Zn2sAHAPwIKX0xaidVKyTgxButxsNDQ3YuXMn3n77bSxevBibNm3CxRdfzIe4Op0OPT09ATdpMyEPGzhKSkoKuZXIVJbCxS9yo6enB+Pj46isrPQgNkopzGYz3/lgYjJ/04liEe0uCKXUYxbE14BYZmYmKKVobm7GU089hTfffHM+75YgAmJIBLAPwFOU0u2Cny+ilA5G9aTmEzkIwb5AO3fuxN69e/nJu7i4ODz22GOiOx/si8dk3KySnpubG/QYLP+Oxso9dq4dHR2w2+2iahrC4TCVSsUX/qQULRkxyOFI7QuUUrS0tCApKcmv/4PVakVLSwtuv/129PX14aabbsJVV12FFStWyH5+cmAmfeAIIXEA/gbgIKX0vjk/r/lKDkJwHIebbroJX3zxBeLj45GamopNmzZh48aNyM3NlZT/C+c9NBoN30r0vvMyG/cVK1ZE5Y5FKcXRo0dBCPHrvxgI3rZ43guBfIG5e1dVVUWlPcoco9gm70D4xz/+gYcffhjbtm3DZ599BofDgWuuuSZi52Kz2bBu3TrY7Xa4XC5cdtll+J//+Z+IHZ/B7XZDpVIRQogK08TQQin9L8Cz/jAXOCXIwWw249lnn8Xtt98OQgg6Ozuxa9cu7NmzB3Fxcdi0aRM2b96MBQsWSLqohJoDdufNy8uD1WrFsWPHoqYpYFbu7KKJhCMz6xAwjUheXh7fIQBO6iaiYVUPSLOSe++99/DAAw9g7969sjmDsZZrSkoKnE4nvvSlL+Hxxx/3sDgMBTNk4PEzQogawDsAmiilP5z52ZwSA3CKkIM/MAOXXbt24a9//Svcbjc2btyILVu2oLCwUNJFZrPZoNfrMTAwAKvViqKiIhQUFMiuK2Du12zEW47jC4fDMjMzkZSUhIGBAdTU1ERFN8GMaLRabdBOz0cffYR7770Xe/fuxYIFC2Q/N2D6JvGlL30JTz31FM4666yQjyMkhoGBAVBK2fewDsCFlNJHgNggBuAUJwchKKUYGhrCrl27sHv3blitVlxyySXYvHmzaB+H/v5+DA8PY/ny5bx9nMvl4iOKSEcRbAQ6Pz8/Kt4IbrcbfX196O7uhlar5Vu/kbTF8wZLlzQaTVBi+PTTT3H33XfjzTffjMrIO8dxWL16NTo6OnDLLbfgl7/8ZUSO+5//+Z/o6+uDxWJBTU0NHnzwQWErM6oqyEA4bcjBG3q9Hrt378brr78Oo9GIiy++GFu2bPHpzEQpRXd3N0wm06wOAQvRdTod7HY7TxRS3ZO8wazVo9kFYcttmNGtsKPDdmKIKdSKBSMGtVod1K7+0KFDuP322/HGG29IMq2JBMbHx/GNb3wDTzzxBFatWhXSMTiOg1qtxhNPPIGPPvoIL7/8MrZu3Yq8vDw8/fTTMTmyfdqSgxCjo6PYs2cPdu3aheHhYVx00UX4xje+geXLl/MKvfj4+KC7Mlwul8dyX1+5vBiwuYWysrKAOxciiZGREXR2dqK2ttZn25NpKdhCoHC9Jpk1HyEkqFVeQ0MDvv/972P37t0RdbCWggcffBBJSUm4++67JT2vq6vL45y3b9+O5cuX49VXX0VfXx9effVVEEKSAdRSSj+O8GmHBYUcvDA+Po6//e1veP3119HZ2QmNRoOvfOUruP/++yWF1t65PBtfzsjICHghsPZoKHMLocJgMKC7u1u0Nb7Qa9LtdvOdD7FpFZOZU0qDdl6OHDmCG2+8ETt37kRFRYXo9xQuDAYDtFotMjIyYLVa8dWvfhU//vGP8fWvf130MXbt2oVrrrkGN998M379618DAJ544gncc8892LJlCz+xSQh5EEAqgLtjodbAoJCDH5jNZmzevBklJSUwm804duwYvvKVr2Dz5s1Ys2aNJJkys1nT6XQBt4Azdyq55Mm+oNPp0Nvbi5qampDShUC2eL4uekopjh8/DrfbHZQY2tracN1112H79u1R1zA0NTXhmmuu4Zc6f+tb38L9998v+vlGoxFbtmzB1q1b8fbbb6OsrAyPP/44AOCOO+7Ahx9+iMcffxx79+7Fr371q08BXEIpHZPp7YQEhRz8YGpqCh9//DG++tWvAvD0pGhqasJ5552HzZs3Y+3atZIiCrZjQafTYXx8nPc50Gq1/Mh1tCTAw8PD/DKdSNQRhLZ4bDqWWcgxK/yOjg64XC5+ia4/tLe34+qrr8Yrr7yCysrKsM9tLnDs2DGUlJRgaGgIN998M0pLS/H0008DAB566CGYTCZMTU3h97//fS6ldGSOT3cWIkIO+/btwx133AGO43DjjTfyZp4MdrsdV199NT7//HNkZ2fjtddek6UtFy3YbDa888472LlzJz7//HOcc845+MY3voFzzz1Xkvya+RwwE5qsrCwsWrQIOTk5sjsXsT2ZNTU1sizTYdESm45NS0uDy+WCRqPBihUrAhJDd3c3vv3tb+OFF15AbW1txM9tLtDf34/vfe97KCwsxJ/+9Cc0NjYKJ01PzYIkx3GoqKjAO++8g8LCQtTX12Pbtm0eYeAf/vAHNDU14emnn8b27duxe/duvPbaaxE4/bmHw+HAe++9h507d+LTTz/lPSnWrVsnKn9nhcCqqio4HA6P7gAr+kX64h0YGMDw8DBqampkJyHgpIjLZDKBEBJwiKq3txdXXHEFnnnmGdTX18t+btHE4OAgfvSjH+HQoUPQarX46KOPmFHOqUkOn376KR544AH8/e9/BwA8/PDDAKZtsxguuugiPPDAAzj77LPhcrmwYMECGAyGWNk2FDG4XC58+OGH2LFjBz766CPek+LLX/6yz6o+W73nXQj0Zdufn58fkTYi8+ysrq6OCjEAQGdnJ2w2G3/DEA5RMYl6fHw8HA4HvvWtb+H3v/89v0/1VMMvfvELPPnkk/jggw+EEvGYvBDCviUNDAx42HYVFhbis88+8/sYjUaD9PR0jI6ORmVYKZrQaDS44IILcMEFF4DjOHz88cfYtWsXHnjgAaxcuRJbtmzBhRdeiKSkJLS0tMBms6Gurm5WZEAIQUpKClJSUrBkyRJ+gOrw4cP8xRSKx2Rvby+MRmNQN/BIoqurC1arFStXruRvBqmpqUhNTcWSJUtgtVqh0+nw3e9+F319fbj00kujpnyMNkZHR/HPf/4T//znP4POjsQC5N3cehpDrVZj3bp1WLduHdxuNw4cOICdO3fi4YcfRkpKCtRqNV577TVRKUNycjJKS0tRWloKq9UKvV6PpqYmEEL4wbBgeoOenh5MTEygqqoqasTQ3d3NrzT0FyUmJiYiKSkJHMfhd7/7HcbHx/HQQw/hmWeeiWhk2dfXh6uvvho6nQ6EENx000244447InZ8McjOzsabb745b1yqwiaHgoIC9PX18X/v7++fJW1ljyksLITL5cLExIRsi2RiESqVCmvXrsXatWvx4IMP4t///jdqa2tx8cUXo6ioiPekENO+TExMRHFxMYqLi3nb/paWFnAcxxOF96BUV1cXJicnUVlZGVViMJvNAYkBmK65XH755XjooYewYcMGAMD3vve9iJ+PRqPBr3/9a9TV1cFsNmP16tVYv3591Fuk84UYgAjUHFwuFyoqKvCvf/0LBQUFqK+vx6uvvoqVK1fyj/n973+P5uZmviD5+uuv4y9/+UsETn/+Yd++fVi/fj3UajXvSbFjxw689dZbyM3NxebNm/H1r39dsgCK6Q10Oh2cTicvddbpdPwm70jciR/7Vyee/6QXbgqUZCfitRvWICXB8x7D7OtWrVoVkIzGxsZw6aWX4r777sPGjRvDPjcp2Lx5M2699VasX78+qq/rBzFZc4hIK/Ott97CnXfeCY7jcP311+MnP/kJ7r//fqxZswabNm2CzWbDVVddhcOHDyMrKwvbt28PKoMN1h597LHH8Mwzz/Abrp977rmoWbTJATZnsHPnTt7VaNOmTfj6178u2ZOCzXt0d3fD6XSioKDAr22/FOz8YhD3v3nM42el2YnYe8vJMeYTJ07wLlWBiGFiYgLf/OY3cffdd+PSSy8N+ZxCQU9PD9atW4cjR45ExfJOBE5dcog0xLRH33vvPZx11llISkrCU089hffff/+UaY9SSj08KeLj47Fx40bRnhRMnux2u1FeXs6rM4Wu1aFsofr2c5+jod/THZ0AaLn/ywCmC55jY2NBicFsNuOyyy7DrbfeiiuuuELSOYSLyclJnHfeefjJT34SdVIKgJgkh+gkoBJx4MABlJeXo6ysDHFxcdi6dSv27Nnj8Zgvf/nLfG69du1a9Pf3z8WpygJCCMrLy/HjH/8YH3/8MV544QUAwLXXXosNGzbgiSeeQF9fH3wRO4tAAGDZsmXQarXIz89HVVUVzjzzTGRmZqK/vx/79+/H0aNHYTQafR7HF1ITZpeoGL+wTkgwYpiamsLWrVtx0003RZ0YnE4nvvnNb+I73/lOLBFDzCImycFXe3RgYMDv45999ll87Wtfi8apRR2EEBQXF+Ouu+7Chx9+iNdeew1JSUn4j//4D1x44YX4zW9+g66uLlBKwXEcGhoaoNFofE46qtVq5ObmYtWqVTjrrLOQm5uL4eFh7N+/H62trRgdHYXb7X/u596LlkLldY+7Zu1i9PX1YXR0NGgnxGq1YuvWrbjqqqtw1VVXhfW5SAWlFDfccAOWL1+Ou+66K6qvPV8hqlvhcrnw8ccfS5YHRwMvv/wyDh06hA8++GCuT0V2EEKwaNEi3Hbbbbj11lt5T4q77roLY2Nj0Gg0OP/88/GTn/wkaMqgUqmQnZ2N7OxsUEp585r29nakpqYiPz9/1l7L4uwk7P3BWfjft9sxaeewpXoBzl1AeVFVIGKw2Wz4zne+g8svvxzXXXddxD4Tsfj444/x0ksvobKyEjU1NQCAn//857j44oujfi7zBaJqDh0dHbj55puxf/9+vPjii7jssstkPSkxqksA+Oc//4nbbrsNH3zwgWxegvMBLpcLV1xxBdxuNxwOB3Q6nYcnhZTagnCv5ejoKJKTk5Gfn+9z3oPtzgimtnQ4HPjud7+Lr371q7jtttvmtTL26aef5oenJiYmUFJSgvfeey/cw8bkByKKHMxmM7RaLc477zw88sgjvLCHTdpFuncupj16+PBhXHbZZdi3bx+WLl0a0defbzCbzdi9ezeuvvpqANOeFG+88QZef/11nDhxAuvXr8eWLVskC6C892AkJCTwRKHX66HT6YISg9PpxLXXXotzzz0XP/rRj+Y1MQjhdDpxwQUX4L/+678i0YaNyQ9FdLdi//79+MEPfoAvvvhC/rNC8PbohRdeiObmZt5CraioCG+88UbAYwZrjzLs2rULl112GQ4ePIg1a9ZE/L1FE2azGXv37sWuXbvQ3t7Oe1KsXr1aMqmzmYjBwUG4XC4sWbIE+fn5fgfMXC4XbrjhBtTW1uKee+45ZYgBAH7wgx8gNzc3Unb1MfnBiCaHn/3sZxgaGsIf//hHTE5O8uKdjRs34vzzz5f9RMOFmPYoMH0xXXLJJXA4HHjyySfnPTkIYbFY8Pbbb2PXrl04cuQI70lx1llniR7CGhoawuDgICoqKnjvBmYbxwaogOnP+z/+4z9QXl6OBx544JQihhdeeAE7duzA3/72t0hFzTH54QQlB7fbDbvdjg0bNuCee+7Bhg0bcMcdd6CgoAClpaV48cUXsXr1alkWfkQSYusYd955J9avX49HHnkEjz766ClFDkIIPSm++OIL3pPinHPO8Vt0ZsTgPerNbPv1ej3GxsbwySefoK+vD8XFxfj5z39+ShHD559/jmuuuQYfffQRMjMzI3XYmPyAgtKeSqXC8PAwTCYTvvrVr8LhcGD79u3YsWMH7HY7Hn30UXR3d6O9vZ1/Dsdxs3rn7O9WqxXDw8MRfhvBIaY9+sUXX6Cvrw+XXHJJtE8v6khISMDGjRvx4osv4vPPP8c3vvEN7Nq1C+eccw5uu+02/Otf/4LD4eAfPzw8jIGBAZ81hoSEBBQVFWHNmjWorq7G0aNHsX//fnzwwQd4/vnnZTn/66+/Hnl5eSG7QYeKJ598EkajEV/+8pdRU1ODG2+8MaqvH00E7EsePXoUv/3tbzEyMoLa2lqoVCq89957KCwsxI4dO/Doo4/ygpznn38eFosFSUlJHl8ejuNACOHDrw8++AA33ngjCgsLYypsd7vduOuuu3jB0emEuLg4bNiwARs2bPDwpLjnnntQV1eH/Px8mM1m/OpXvwrYyna73Xj00UexePFi7N69G+Pj4+jp6ZHlnK+99lrceuutfBE2WpCL7GIRASOHgoICVFRUYGhoCNu2bcMzzzyDgYEBLF++HCUlJXjyySfx2Wef4fDhw1Cr1fjjH/+I2tpa3HfffWhsbAQwLbwR5mV1dXXYsWMHEhMTYTJNS3E5bnqHx+TkpFzvM+j0KNt9ef7556OkpAT79+/Hpk2bcOjQIdnOKRbBPCmeeuop3gV7+/bt+Oyzz3DzzTfjjTfegMVimfU8t9uN+++/Hw6HA7/73e+gUqmQlZWFuro6Wc5z3bp1UXPnPm1BKQ30h4dOp6MHDx6kTqeTfve736WXX3453b17N21ubqY2m41/3OHDh+lPf/pTevbZZ9MzzzyT3nXXXbSjo0N4KPqPf/yDbt68mdrtdv5n7777Lr388stpZWUl/e1vf+vxb5GA0+mkpaWltKuri9rtdlpVVUWPHDni9/HnnXcePXjwYETPYb7B6XTS6667jo6Pj1OO4+inn35K77rrLlpVVUW/+c1v0pdeeonqdDo6OTlJf/zjH9NrrrmGulyuqJ1fd3c3XblyZdReT0YEuw7n5I9ochDCZDLR5557jm7evJk++OCDlOM4n49ra2ujP/vZz+i+ffv4n5nNZvrzn/+c3n333fzPurq6aHV1NX333XfpwYMH6YUXXkgHBgb4fzcYDP5ORRL27t1Lly5dSsvKyuj//u//Ukop/elPf0r37Nkz67FiyeHtt9+mFRUVdMmSJfThhx/2+ZjXXnuNLl++nK5YsYJeeeWV4b2JGADHcfTQoUP0xz/+Ma2pqaErVqygW7ZsiSoxUKqQg9x/QiIHIby/EG63O+CXpKuri1555ZX0rbfe4n/G7jqUTt+t7rzzTrpr1y7qdrvpe++9Ry+44AJaWVlJH3roIep2u8WcVlTgcrloWVkZ7ezs5KORlpYWj8e0t7fTmpoaajQaKaXTEdipBI7j6J49e6jZbI76ayvkIO+fsJu03pVrQgj/M0rprEGegYEBmEwmfOlLX+J/9uGHH+L666/nj2c0GpGQkIB//vOfeP755/HDH/4QL730Etra2jA4OBjuKUcMYqZH//SnP+GWW27h216nmsxbpVJh06ZNUdu1oSB6kHUqU9jffvfdd/HII49g3759WL16NVJTUwFMW5j19fVh3bp1AKbdgT7//HOsW7cOf/jDH7BhwwZceOGFqK6uhsVi4QesKJ37fTti2qPt7e1ob2/Hueeei7Vr12Lfvn3RPs1TEldeeSXOPvtsHDt2DIWFhXj22Wfn+pROOcg+Ysk6FWlpadi/fz/+/e9/g+M4nH/++fjyl7+Mjz/+mN9PYLfb8c4776C8vBwAcPz4cWzYsIH33fvss8/w6KOPyn3KEYXL5cLx48fx/vvvo7+/H+vWrUNzczMyMjLm+tTmNbZt2zbXp3DKI2p+DmvWrMGuXbug0+nw/PPPw+12w+Vy4f/9v/+HlJQUdHV14b333sNbb72FH/zgB/jggw9QV1fHh+MffvghFi9ejNLSUgCICdWdGHPdwsJCbNq0CVqtFqWlpaioqMDx48ejfaoKFEjGnJi9bNy4EV/5yleg0WiwePFiFBYW4pJLLsEvfvELfOc738H69esxOTmJvLw8OJ1OANPik/POO28uTtcv6uvrcfz4cXR3d/PK0U2bNnk8ZsuWLXj//fcBTDstt7e3z9kaeQUKJCFIxTKq6O/v5/+/p6eHrl+/nra1tdGPP/6Yrlq1inZ2dlJKaUx1LIK1R91uN/3hD39Ily9fTletWkW3bdsm6rjBWqQnTpyg559/Pq2pqaGVlZV07969kXtTCqKNOe9M+PoTU+QghNVqpQ899BCtqKigl156KX377bfn8nSiCjEt0u9973v0D3/4A6WU0paWFlpcXDwHZ6ogQphzIvD1J7Y83wRISEjAvffei3vvvRfj4+OnVQFP2CIFwLdIhePlhBBefj4xMYFFixbNybkqOHURkwaz3jidiAEQ1yJ94IEH8PLLL6OwsBAXX3wxnnjiiWifZkSxb98+nHHGGSgvL8cvfvGLuT4dBZgn5KBgNrZt24Zrr70W/f39eOutt3DVVVcFdI6OZXAch1tuuQVvv/02WltbsW3bNrS2ts71aZ32UMghBiGmRfrss8/iW9/6FgDg7LPPhs1mw8jISFTPM1IQozRVEH2c9uRw8OBBVFVVwWazYWpqCitXrsSRI0fm9JzEtEiLiorwr3/9CwDQ1tYGm82G3NzcuTjdsCF1T4mC6CBmC5LRQn19PTZt2oT77rsPVqsV3/3ud6PuLuQNjUaDJ598EhdddBFvsLty5UoPg91f//rX+N73voff/OY3IITghRdeiAlhmIJTB8E8JE8LEELiABwEYANwDqWUm+NTkgWEkOcAfB2AnlI6iwHJNLs8DuBiABYA11JKZbcbJ4ScDeABSulFM3+/BwAopQ/L/doK/OO0TytmkA0gBUAqgIQ5Phc58QKADQH+/WsAls78uQnAU1E4J2CamJcSQkpniHorgMB7BhTIDoUcpvF/AH4K4BUAv5zjc5ENlNIPARgDPGQzgD/PCHP2A8gghCyMwnm5ANwK4O8A2gD8hVLaIvfrKgiM077mQAi5GoCTUvoqIUQN4BNCyAWU0nfn+tzmAAUA+gR/75/52ZDcL0wpfQvAW3K/jgLxOO3JgVL6ZwB/nvl/DsBZc3tGChTEBpS0QoEQAwAWC/5eOPMzBachFHJQIMQbAK4m01gLYIJSKntKoSA2cdqnFacTCCHbAJwPIIcQ0g/gZwC0AEApfRrTOf/FADow3cq8bm7OVEEsQNE5KFCgwCeUtEKBAgU+oZCDAgUKfEIhBwUKFPiEQg4KFCjwCYUcFChQ4BMKOShQoMAnFHJQoECBTyjkoECBAp/4/2RO87/mDLGXAAAAAElFTkSuQmCC\n",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {
|
|
"needs_background": "light"
|
|
},
|
|
"output_type": "display_data"
|
|
},
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"[<matplotlib.lines.Line2D at 0x7f8f553ca198>]"
|
|
]
|
|
},
|
|
"execution_count": 9,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEQCAYAAACgBo8fAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAkQUlEQVR4nO3dd3gVdf7+/+crjdA70nvvJXRILHQVEEGxYAMpgpToruvuuuuqP9fy2VCUakOwoYICSrckdAm9I70Kofcm7+8fiftjEUyAJHPK/biuc3nOmWHOa0i4ncw5ucecc4iIiP8L8XoAERFJHwp0EZEAoUAXEQkQCnQRkQChQBcRCRAKdBGRAOFpoJvZ+2Z2wMzWpMO2bjOzFZfdzppZx3QYU0TEL5iXn0M3s2jgJDDOOVc9HbebD9gMFHfOnU6v7YqI+DJPj9CdcwnA4cufM7NyZjbDzJaa2Vwzq3wDm+4MTFeYi0gw8cVz6GOAp51z9YBngRE3sI2uwKfpOpWIiI8L83qAy5lZDqAJ8IWZ/fZ0lpRlnYCXrvLH9jjnWl+2jSJADWBmxk4rIuJbfCrQSf6J4ahzrvaVC5xzk4BJadjGfcBXzrkL6TybiIhP86lTLs6548A2M+sCYMlqXedmHkCnW0QkCHn9scVPgYVAJTPbbWbdgYeA7ma2ElgLdLiO7ZUGSgDxGTCuiIhP8/RjiyIikn586pSLiIjcOM/eFC1QoIArXbq0Vy8vIuKXli5detA5V/BqyzwL9NKlS5OYmOjVy4uI+CUz23GtZTrlIiISIBToIiIBQoEuIhIgFOgiIgFCgS4iEiAU6CIiAUKBLiISIPwu0A+dPMdLU9dx/KzKFEVELud3gT5/yyHGLthGy7h45qzb7/U4IiI+w+8CvX2tonz1VFPyZougx7hE+n+6nEMnz3k9loiI5/wu0AFqlcjDlH7NGNSiItPX7KNFXDyTV+xBzZEiEsz8MtABIsJCGNCiAt/2b06p/NkZ8NkKun+YyN6jZ7weTUTEE34b6L+peEtOJvZpwt/vrMKCLQdpNTiBjxfv4NIlHa2LSHDx+0AHCA0xejQvy6yBMdQsnpu/fbWGB95ZxLaDp7weTUQk0wREoP+mZP5sfNyjIa/fW4N1+47TZkgCYxK2cPHXS16PJiKS4QIq0AHMjPvrl2RObAzRFQvy6rQNdBq5gPX7jns9mohIhgq4QP/NLbkiGdOtHm8/WIc9R85w91vziJu1kXMXf/V6NBGRDBGwgQ7JR+t31SzKnNgY7q5VlGHfb+auYfNYtvOI16OJiKS7gA703+TNHsHg+2vzwWP1OXnuIveOXMBLU9dx+vxFr0cTEUk3QRHov7mtciFmDYrmoYYleX/+NloPSWD+5oNejyUiki6CKtABckaG80rHGkzo2YiwkBAeencxz325imNnVPYlIv4t6AL9Nw3L5mf6gOb0jinHl8t20zIunplrf/F6LBGRGxa0gQ4QGR7KX9pW5uunmpI/RxZ6jV9K34+XkXRCZV8i4n+COtB/U6N4bqb0a8qzrSoye91+Wg6OZ9Ky3Sr7EhG/okBPER4aQr/bKzBtQDPKFshO7OcreXzsEvao7EtE/IQC/QrlC+Xki95N+OfdVVm89TCt4uIZv3C7yr5ExOcp0K8iNMR4vGkZZg2Kpm6pvLwweS1dxyxia9JJr0cTEbkmBfofKJEvG+OeaMCbnWuy4ZfjtBk6l5E/quxLRHyTAj0VZkaXqBLMiY3htkoFeX3GBjqOmM/avce8Hk1E5H+kGuhmVsLMfjCzdWa21swGXGUdM7NhZrbZzFaZWd2MGdc7hXJFMrpbFCMfqssvx87R/u35vDlzA2cvqOxLRHxDWo7QLwLPOOeqAo2AvmZW9Yp12gIVUm49gZHpOqUPaVujCHNio+lYuxjDf9jCncPmsnTHYa/HEhFJPdCdc/ucc8tS7p8A1gPFrlitAzDOJVsE5DGzIuk+rY/Iky2C/9xXiw+faMDZC5foPGohL05Zy6lzKvsSEe9c1zl0MysN1AEWX7GoGLDrsse7+X3oY2Y9zSzRzBKTkpKuc1TfE1OxIDMHRfNIo1J8uHA7rQYnkLDJ//dLRPxTmgPdzHIAE4GBzrkbuvyPc26Mcy7KORdVsGDBG9mEz8mRJYx/dajO570akyU8hEfe/4lnv1jJ0dPnvR5NRIJMmgLdzMJJDvOPnXOTrrLKHqDEZY+LpzwXNOqXzse0/s156tZyfLV8Dy3iEpi+ep/XY4lIEEnLp1wMeA9Y75yLu8ZqU4BHUj7t0gg45pwLujSLDA/lz20qM7lvUwrlzEKfj5fR56OlHDhx1uvRRCQIWGoFVGbWDJgLrAZ++42avwIlAZxzo1JC/22gDXAaeNw5l/hH242KinKJiX+4il+78OslxiRsZeh3P5M1PJS/31mFzvWKk/xXJSJyY8xsqXMu6qrLvGoUDPRA/83mAyf5y8RVJO44QvMKBXj1nhqUyJfN67FExE/9UaDrN0UzWPlCOfi8V2Ne6lCNZTuO0HpIAmPnb1PZl4ikOwV6JggJMR5pXJqZg6KJKp2PF6eu477RC9l8QGVfIpJ+FOiZqHjebHz4eH3+06UWPx84Sbuhcxn+w2YuqOxLRNKBAj2TmRn31ivOnNgYWlQtxJszN9Lh7fms2aOyLxG5OQp0jxTMmYURD9Vj1MP1SDp5jg7D5/P6DJV9iciNU6B7rE31wswZFMO9dYsx8scttBs6lyXbVfYlItdPge4DcmcL543Otfioe0PO/3qJLqMW8o/Jazipsi8RuQ4KdB/SrEIBZg6M5vGmpRm/aAet4uL5YeMBr8cSET+hQPcx2bOE8c+7q/Fl7yZkyxLG4x8sIXbCCo6cUtmXiPwxBbqPqlcqL9/2b8bTt5dnysq9tBwcz7er9uHVb/aKiO9ToPuwLGGhPNOqElP6NaNI7qz0/WQZvcYv5cBxlX2JyO8p0P1A1aK5+OqpJjzftjLxm5K4Iy6ez5fs0tG6iPwPBbqfCAsNoVdMOaYPaE6VIrn488RVdHvvJ3YeOu31aCLiIxTofqZswRx89mQjXulYnRW7jtJ6SALvzdvGryr7Egl6CnQ/FBJiPNyoFLMGRdOwbD5e/mYdnUct4Of9J7weTUQ8pED3Y0XzZOWDx+oz5P7abD94ijuHzWPYdz9z/qLKvkSCkQLdz5kZHesUY3ZsDK2rFyZu9ibavz2PVbuPej2aiGQyBXqAKJAjC289UId3HoniyOnzdBw+n39PW6+yL5EgokAPMC2r3sKsQTHcX78EoxO20mZIAou2HvJ6LBHJBAr0AJQ7azj/7lSTT3o05JKDrmMW8bevVnPi7AWvRxORDKRAD2BNyhdgxsDm9GhWhk9/2kmrwQl8v2G/12OJSAZRoAe4bBFh/P2uqkzs04SckWE8MTaRgZ8t57DKvkQCjgI9SNQpmZdvnm7OgDsq8O3qfbSIi2fKyr2qDxAJIAr0IBIRFsKglhWZ+nQzSuTNSv9Pl/PkuKX8ckxlXyKBQIEehCoXzsWkp5ryt3ZVmLc5iZZx8Xz6004drYv4OQV6kAoNMZ6MLsuMAdFUK5aL5yet5sF3FrPj0CmvRxORG6RAD3KlC2Tnkx6NePWeGqzZc4zWQxJ4d+5WlX2J+CEFuhASYjzYsCSzYqNpWq4Ar3y7nk4jF7DxF5V9ifgTBbr8V5HcWXn30SiGPVCHXYdPc9dbcxkyZ5PKvkT8hAJd/oeZ0b5WUebExtCuRhGGzPmZu9+ax4pdR70eTURSoUCXq8qXPYKhXevw3qNRHDtzgU4j5vPKN+s4c15lXyK+SoEuf+iOKrcwKzaarg1K8u68bbQeksCCLQe9HktErkKBLqnKFRnOq/fU4NMnGxFi8OA7i3l+0iqOq+xLxKco0CXNGpfLz/QB0fSKLsuEJbtoGRfPnHUq+xLxFakGupm9b2YHzGzNNZbfambHzGxFyu0f6T+m+IqsEaE8364KX/dtSt5sEfQYl8jTny7n0MlzXo8mEvTScoQ+FmiTyjpznXO1U24v3fxY4utqFs/DlH7NiG1ZkRlrksu+vl6+R/UBIh5KNdCdcwnA4UyYRfxMRFgI/e+owLf9m1Mqf3YGTlhB9w8T2Xv0jNejiQSl9DqH3tjMVprZdDOrdq2VzKynmSWaWWJSUlI6vbR4reItOZnYpwkv3FWVhVsO0WpwAh8t2sEl1QeIZCpLy4/IZlYa+MY5V/0qy3IBl5xzJ82sHTDUOVchtW1GRUW5xMTEGxhZfNnOQ6d5/qtVzN98iIZl8vHavTUpUyC712OJBAwzW+qci7rasps+QnfOHXfOnUy5Pw0IN7MCN7td8U8l82fjo+4NeePemqzbd5w2QxIYHb+Fi7+qPkAko910oJtZYTOzlPsNUrapy8wHMTPjvvolmBMbQ3TFgvx7+gY6jVzA+n3HvR5NJKCl5WOLnwILgUpmttvMuptZbzPrnbJKZ2CNma0EhgFdnT7qIMAtuSIZ060ewx+sy96jZ7j7rXnEzdrIuYuqDxDJCGk6h54RdA49uBw5dZ6Xv1nHpOV7KF8oB6/fW5N6pfJ6PZaI38nQc+giaZE3ewRx99fmg8frc/rcRTqPWsC/pq7l9PmLXo8mEjAU6JKpbqtUiFmxMXRrVIoP5m+n1eAE5v2ssi+R9KBAl0yXI0sYL3Wozue9GhMeGsLD7y3mz1+u5NgZlX2J3AwFunimQZl8TB/QnD63lmPisj20jItn5tpfvB5LxG8p0MVTkeGhPNemMl8/1ZT8ObLQa/xS+n68jKQTKvsSuV4KdPEJNYrnZkq/pvypdSVmr9tPi7h4Ji7drbIvkeugQBefER4aQt/byjNtQDPKF8rBM1+s5LEPlrBHZV8iaaJAF59TvlBOvujVmBfvrsqS7YdpFRfPuIXbVfYlkgoFuvikkBDjsaZlmDkwmrql8vKPyWu5f8xCtiSd9Ho0EZ+lQBefViJfNsY90YA3O9dk4y8naDt0LiN+3MwFlX2J/I4CXXyemdElqgRznonh9kqFeGPGRjoOn8+aPce8Hk3EpyjQxW8UyhnJqG71GPlQXfYfP0eH4fN5c+YGzl5Q2ZcIKNDFD7WtUYQ5sdHcU6cYw3/YQrthc0ncrqskiijQxS/lyRbB/3WpxbgnGnDuwiW6jF7Ii1PWcuqcyr4keCnQxa9FVyzIrEHRPNq4NB8uTC77Stik69VKcFKgi9/LniWMF9tX44tejckSHsIj7//Es1+s5Ojp816PJpKpFOgSMKJK52Na/+b0va0cXy3fQ4u4BKav3uf1WCKZRoEuASUyPJQ/ta7MlH5NuSVXFvp8vIze45dy4PhZr0cTyXAKdAlI1YrmZnLfpjzXpjLfbzxAi7h4vkjcpbIvCWgKdAlYYaEh9Lm1HNMHNKdS4Zz86ctVPPL+T+w6fNrr0UQyhAJdAl65gjmY0LMxL3eoxrIdR2g9JIGx87ep7EsCjgJdgkJIiNGtcWlmDoqmful8vDh1HV1GL2TzgRNejyaSbhToElSK583G2MfrE3dfLbYknaTd0Hm8/f3PKvuSgKBAl6BjZnSqW5zZg2JoWe0W/m/WJtq/rbIv8X8KdAlaBXNmYfiDdRndrR4HTyaXfb02XWVf4r8U6BL0WlcrzJxBMXSuW5xR8VtoN3QuP21T2Zf4HwW6CJA7Wzivd67JR90bcv7XS9w3eiEvfL2GE2cveD2aSJop0EUu06xCAWYNiuaJpmX4aPEOWg9O4IeNB7weSyRNFOgiV8gWEcY/7q7Kl72bkD1LGI9/sITYCSs4ckplX+LbFOgi11CvVF6+6d+M/reXZ8rKvbSIi+ebVXtVHyA+S4Eu8geyhIUS26oSU59uRtE8Wen3yXJ6jV/KfpV9iQ9SoIukQZUiufjqqSY837Yy8ZuSaBEXz4QlO3W0Lj5FgS6SRmGhIfSKKceMgdFUKZKL5yau5uH3FrPzkMq+xDco0EWuU5kC2fnsyUa80rE6K3cdo/WQBN6bt41fVfYlHlOgi9yAkBDj4UalmDUomsbl8vPyN+u4d+QCNu1X2Zd4J9VAN7P3zeyAma25xnIzs2FmttnMVplZ3fQfU8Q3Fc2TlfcejWJo19rsOHSKO4fNZdh3P3P+osq+JPOl5Qh9LNDmD5a3BSqk3HoCI29+LBH/YWZ0qF2MObExtKlehLjZm2j/9jxW7jrq9WgSZFINdOdcAvBHxRYdgHEu2SIgj5kVSa8BRfxF/hxZeOuBOrzzSBRHTp/nnhHz+fe09Zw5r7IvyRzpcQ69GLDrsse7U577HTPraWaJZpaYlJSUDi8t4ntaVr2F2bEx3F+/BKMTttJ2aAILtxzyeiwJApn6pqhzboxzLso5F1WwYMHMfGmRTJUrMpx/d6rJJz0acsnBA+8s4q9frea4yr4kA6VHoO8BSlz2uHjKcyJBr0n5AswcGM2Tzcvw2U87aRWXwPcb9ns9lgSo9Aj0KcAjKZ92aQQcc87tS4ftigSErBGh/O3Oqkx6qim5s4bzxNhEBny2nEMnz3k9mgSYtHxs8VNgIVDJzHabWXcz621mvVNWmQZsBTYD7wBPZdi0In6sdok8TH26GQNbVGDa6n20HJzAlJUq+5L0Y159M0VFRbnExERPXlvEaxt/OcGfJ65i5a6jtKhSiFc61qBw7kivxxI/YGZLnXNRV1um3xQV8UClwjmZ1KcJf7+zCvM2H6RlXDyfLN7JJdUHyE1QoIt4JDTE6NG8LDMHRlO9WG7++tVqHnx3EdsPnvJ6NPFTCnQRj5XKn51PnmzIa51qsHbPcdoMTeCdhK0q+5LrpkAX8QFmRtcGJZkdG0Oz8gX4/6atp9OI+Wz8RWVfknYKdBEfUjh3JO88EsVbD9Rh95Ez3PXWXAbP3qSyL0kTBbqIjzEz7q5VlNmxMdxZowhDv/uZu96ay/KdR7weTXycAl3ER+XLHsGQrnV4/7EoTpy9SKeRC3j5m3WcPn/R69HERynQRXzc7ZVvYdagaB5qWJL35m2jzZC5LNh80OuxxAcp0EX8QM7IcF7pWIPPejYixODBdxfzl4mrOHZGZV/y/1Ogi/iRRmXzM2NgNL1iyvJ54i5aDY5n9jqVfUkyBbqIn4kMD+X5tlX4um9T8maL4MlxifT7ZBkHVfYV9BToIn6qZvE8TOnXjGdaVmTW2v20jIvn6+V7VPYVxBToIn4sIiyEp++owLf9m1G6QHYGTljBE2OXsPfoGa9HEw8o0EUCQIVbcvJl7yb8466qLNp6mFaDExi/aIfKvoKMAl0kQISGGE80K8OsQdHULpGHF75eQ9d3FrFNZV9BQ4EuEmBK5MvG+O4NeOPemqzfd5w2QxIYFb+Fi7+qPiDQKdBFApCZcV/9EsyJjSGmYkFem76Be0YsYN3e416PJhlIgS4SwG7JFcnobvUY/mBd9h07Q/u35/GfWRs5d/FXr0eTDKBAFwlwZsadNYswe1AM7WsX5a3vN3PnsHks3aGyr0CjQBcJEnmzRxB3X23GPl6fM+d/pfOoBfxr6lpOnVPZV6BQoIsEmVsrFWLmoGi6NSrFB/O303pIAnN/TvJ6LEkHCnSRIJQjSxgvdajO570aExEaQrf3fuLPX67k2GmVffkzBbpIEGtQJh/TBjSnz63lmLhsDy0GxzNjzS9ejyU3SIEuEuQiw0N5rk1lJvdtSsEcWej90VL6fryMpBMq+/I3CnQRAaB6sdxM7teUP7WuxOz1+2kRF8/EpbtV9uVHFOgi8l/hoSH0va080/o3p3yhHDzzxUoe/WAJu4+c9no0SQMFuoj8TvlCOfiiV2P+1b4aidsP03pwAuMWblfZl49ToIvIVYWEGI82Kc3MgdHULZWXf0xey/1jFrIl6aTXo8k1KNBF5A+VyJeNcU804P+61GLT/pO0HTqXET9u5oLKvnyOAl1EUmVmdK5XnNmx0bSoUog3Zmyk4/D5rNlzzOvR5DIKdBFJs0I5IxnxUD1GPVyX/cfP0WH4fN6YsYGzF1T25QsU6CJy3dpUL8J3sTF0qlOMET9uod2wuSRuP+z1WEFPgS4iNyR3tnDe7FKLcU804NyFS3QZvZB/Tl7DSZV9eUaBLiI3JbpiQWYNiubRxqUZt2gHrQcnEL9JZV9eSFOgm1kbM9toZpvN7C9XWf6YmSWZ2YqUW4/0H1VEfFX2LGG82L4aX/ZuTGR4CI++/xPPfL6So6fPez1aUEk10M0sFBgOtAWqAg+YWdWrrDrBOVc75fZuOs8pIn6gXql8fNu/Of1uK8/kFXtoERfPtNX7vB4raKTlCL0BsNk5t9U5dx74DOiQsWOJiL+KDA/l2daVmNyvKYVzR/LUx8voPX4pB46f9Xq0gJeWQC8G7Lrs8e6U5650r5mtMrMvzazE1TZkZj3NLNHMEpOSdI5NJJBVK5qbr59qynNtKvP9xgO0iIvn88RdKvvKQOn1puhUoLRzriYwG/jwais558Y456Kcc1EFCxZMp5cWEV8VFhpCn1vLMWNAcyoXzsWfv1zFI+//xK7DKvvKCGkJ9D3A5UfcxVOe+y/n3CHn3G/lye8C9dJnPBEJBGUL5uCzno14uUM1lu04QushCXwwfxu/quwrXaUl0JcAFcysjJlFAF2BKZevYGZFLnvYHliffiOKSCAICTG6NS7NrNgYGpTJx7+mrqPLqAVsPnDC69ECRqqB7py7CPQDZpIc1J8759aa2Utm1j5ltf5mttbMVgL9gccyamAR8W/F8mTlg8fqM/j+Wmw9eIp2Q+fx9vc/q+wrHZhXb1BERUW5xMRET15bRHzDwZPn+OeUtXy7ah+VC+fkzc61qFE8t9dj+TQzW+qci7raMv2mqIh4pkCOLAx/sC6ju9Xj8KnzdBwxn9emq+zrRinQRcRzrasVZnZsDJ3rFmdU/BbaDp3L4q2HvB7L7yjQRcQn5M4azuuda/Jxj4ZcvHSJ+8cs4oWv13Di7AWvR/MbCnQR8SlNyxdg5sBoujcrw0eLk8u+fthwwOux/IICXUR8TraIMF64qyoT+zQhe5YwHh+7hEETVnD4lMq+/ogCXUR8Vt2SefmmfzP631GBqSv30jIunm9W7VV9wDUo0EXEp2UJCyW2ZUWmPt2MYnmz0u+T5fQcv5T9Kvv6HQW6iPiFKkVyMalPE/7arjIJm5JoERfPhCU7dbR+GQW6iPiNsNAQekaXY+bAaKoWycVzE1fz0LuL2XlIZV+gQBcRP1S6QHY+fbIRr95Tg1W7j9FqSDzvzt0a9GVfCnQR8UshIcaDDUsyOzaaJuUK8Mq367l35AI27Q/esi8Fuoj4tSK5s/Leo1EM7VqbnYdPc+ewuQyd8zPnLwZf2ZcCXUT8npnRoXYxZg+Kpm31Igyes4n2b89j5a6jXo+WqRToIhIw8ufIwrAH6vDuI1EcPX2Be0bM59Vp6zlzPjjKvhToIhJwWlS9hVmx0XRtUJIxCVtpMzSBhVsCv+xLgS4iASlXZDiv3lODT55sCMAD7yzi+UmrOR7AZV8KdBEJaE3KFWDGgGh6RpdlwpKdtIpL4Lv1+70eK0Mo0EUk4GWNCOWv7aow6amm5M4aTvcPE+n/6XIOnTyX+h/2Iwp0EQkatUvkYerTzRjUoiLT1+yj5eAEJq/YEzD1AQp0EQkqEWEhDGhRgW/7N6dkvmwM+GwFPT5MZN+xM16PdtMU6CISlCrekpOJfZrw9zurMH/LQVrFJfDJ4p1c8uP6AAW6iASt0BCjR/OyzBoYQ43iufnrV6t58N1FbD94yuvRbogCXUSCXsn82fi4R0Ne61SDtXuO03pIAmMStnDxV/+qD1Cgi4iQXB/QtUFJZsfG0LxCQV6dtoF7Ry5gwy/HvR4tzRToIiKXKZw7knceqcdbD9Rh95Ez3DVsHnGzN3Huou/XByjQRUSuYGbcXasos2NjuLtWUYZ99zN3vzWP5TuPeD3aH1Kgi4hcQ77sEQy+vzYfPFafE2cv0mnkAl7+Zh2nz1/0erSrUqCLiKTitsqFmDUomocaluS9edtoPSSB+ZsPej3W7yjQRUTSIGdkOK90rMGEno0ICwnhoXcX85eJqzh2xnfKvhToIiLXoWHZ/Ewf0JxeMWX5PHEXLePimbX2F6/HAhToIiLXLTI8lOfbVuHrvk3Jlz2CnuOX0u+TZRz0uOxLgS4icoNqFk8u+3q2VUVmrd1Pi7h4vlq+27OyLwW6iMhNCA8Nod/tFZg2oBllC2Rn0ISVPD52CXuOZn7ZlwJdRCQdlC+Uky96N+Gfd1dl8dbDtIqLZ/yiHZla9qVAFxFJJ6EhxuNNyzBrUDR1Sublha/X0HXMIrYmncyU109ToJtZGzPbaGabzewvV1mexcwmpCxfbGal031SERE/USJfNsZ3b8AbnWuy4ZfjtB06l1HxGV/2lWqgm1koMBxoC1QFHjCzqles1h044pwrDwwGXk/vQUVE/ImZcV9UCebExnBrpYK8Nn0DHUfMZ93ejCv7SssRegNgs3Nuq3PuPPAZ0OGKdToAH6bc/xK4w8ws/cYUEfFPhXJFMrpbFCMfqssvx87R/u15vDdvW4a8VloCvRiw67LHu1Oeu+o6zrmLwDEg/5UbMrOeZpZoZolJSUk3NrGIiB9qW6MIc2Kj6VC7GKXyZcuQ1wjLkK1eg3NuDDAGICoqyn+v8yQicgPyZIvgP/fVyrDtp+UIfQ9Q4rLHxVOeu+o6ZhYG5AYOpceAIiKSNmkJ9CVABTMrY2YRQFdgyhXrTAEeTbnfGfjeefWrUiIiQSrVUy7OuYtm1g+YCYQC7zvn1prZS0Cic24K8B4w3sw2A4dJDn0REclEaTqH7pybBky74rl/XHb/LNAlfUcTEZHrod8UFREJEAp0EZEAoUAXEQkQCnQRkQBhXn260MySgB03+McLAL53hdaMpX0ODtrn4HAz+1zKOVfwags8C/SbYWaJzrkor+fITNrn4KB9Dg4Ztc865SIiEiAU6CIiAcJfA32M1wN4QPscHLTPwSFD9tkvz6GLiMjv+esRuoiIXEGBLiISIHw60IPx4tRp2OdYM1tnZqvM7DszK+XFnOkptX2+bL17zcyZmd9/xC0t+2xm96V8rdea2SeZPWN6S8P3dkkz+8HMlqd8f7fzYs70Ymbvm9kBM1tzjeVmZsNS/j5WmVndm35R55xP3kiu6t0ClAUigJVA1SvWeQoYlXK/KzDB67kzYZ9vA7Kl3O8TDPucsl5OIAFYBER5PXcmfJ0rAMuBvCmPC3k9dybs8xigT8r9qsB2r+e+yX2OBuoCa66xvB0wHTCgEbD4Zl/Tl4/Qg/Hi1Knus3PuB+fc6ZSHi0i+gpQ/S8vXGeBl4HXgbGYOl0HSss9PAsOdc0cAnHMHMnnG9JaWfXZArpT7uYG9mThfunPOJZB8fYhr6QCMc8kWAXnMrMjNvKYvB3q6XZzaj6Rlny/XneT/w/uzVPc55UfREs65bzNzsAyUlq9zRaCimc03s0Vm1ibTpssYadnnF4GHzWw3yddfeDpzRvPM9f57T1WmXiRa0o+ZPQxEATFez5KRzCwEiAMe83iUzBZG8mmXW0n+KSzBzGo45456OVQGewAY65z7j5k1JvkqaNWdc5e8Hsxf+PIRejBenDot+4yZtQD+BrR3zp3LpNkySmr7nBOoDvxoZttJPtc4xc/fGE3L13k3MMU5d8E5tw3YRHLA+6u07HN34HMA59xCIJLkEqtAlaZ/79fDlwM9GC9Oneo+m1kdYDTJYe7v51UhlX12zh1zzhVwzpV2zpUm+X2D9s65RG/GTRdp+d7+muSjc8ysAMmnYLZm4ozpLS37vBO4A8DMqpAc6EmZOmXmmgI8kvJpl0bAMefcvpvaotfvBKfyLnE7ko9MtgB/S3nuJZL/QUPyF/wLYDPwE1DW65kzYZ/nAPuBFSm3KV7PnNH7fMW6P+Lnn3JJ49fZSD7VtA5YDXT1euZM2OeqwHySPwGzAmjl9cw3ub+fAvuACyT/xNUd6A30vuxrPDzl72N1enxf61f/RUQChC+fchERkeugQBcRCRAKdBGRAKFAFxEJEAp0EZErpFasdQPbm2FmR83smyueL5NSLLg5pWgwIuX5GyoeVKCLiPzeWCA96xbeBLpd5fnXgcHOufLAEZI/2kjKf4+kPD84Zb1UKdBFRK7grlKsZWblUo60l5rZXDOrfB3b+w44ccX2DLid5GJBSC4a7Jhy/4aKBxXoIiJpMwZ42jlXD3gWGHGT28sPHHXJxYLwv+VcN1Q8qHIuEZFUmFkOoAnwxWUHyllSlnUi+Tder7THOdc6cyZMpkAXEUldCMlH07WvXOCcmwRMuoFtHiK5Az0s5Sj88nKu34q7dl9P8aBOuYiIpMI5dxzYZmZd4L+Xj6t1k9t0wA8kFwtCctHg5JT7N1Q8qC4XEZErmNmnJLddFiC5DO+fwPfASKAIEA585py72qmWq21vLlAZyEHykXZ359xMMytL8tWb8pF8ycGHnXPnzCwSGA/UIfnN2a7OuVTbNhXoIiIBQqdcREQChAJdRCRAKNBFRAKEAl1EJEAo0EVEAoQCXUQkQCjQRUQCxP8DNGgnpenFdnwAAAAASUVORK5CYII=\n",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {
|
|
"needs_background": "light"
|
|
},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"print(positions[:2])\n",
|
|
"\n",
|
|
"\n",
|
|
"# https://jakevdp.github.io/PythonDataScienceHandbook/04.12-three-dimensional-plotting.html\n",
|
|
"ax = plt.axes(projection='3d')\n",
|
|
"#ax.scatter3D(positions[:][0], positions[:][1], positions[:][2], c=positions[:][2], cmap='Greens')\n",
|
|
"ax.scatter3D(positions[:][0], positions[:][2], positions[:][1], cmap='Greens')\n",
|
|
"ax.set_xlabel('x')\n",
|
|
"ax.set_ylabel('z')\n",
|
|
"ax.set_zlabel('y')\n",
|
|
"plt.show()\n",
|
|
"\n",
|
|
"plt.plot(positions[:][0],positions[:][1])\n",
|
|
"\n",
|
|
"\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.7.3"
|
|
}
|
|
},
|
|
"nbformat": 4,
|
|
"nbformat_minor": 2
|
|
}
|