202 lines
110 KiB
Plaintext
202 lines
110 KiB
Plaintext
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 1,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"import numpy as np\n",
|
|
"import matplotlib.pyplot as plt"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 45,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUIAAAHHCAYAAADK0FH3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAABZcElEQVR4nO3deXhc5Xnw/++Zfd+l0TZavW+AF4QXwNjGxgYHSNI0KUkMyRvepIZAnLZv/F4BAk1iSNuEEEhIeVMgvcKPhLSQJg0QMNiUYHmRrMWWZcm2ZMna932bmfP7Q5awsB1sazyj8dyf6zoX1pmZM89zkG8/y/08R1FVVUUIIRKYJtYFEEKIWJNAKIRIeBIIhRAJTwKhECLhSSAUQiQ8CYRCiIQngVAIkfAkEAohEp4EQiFEwpNAKIRIeBIIxWWnKMoFHbt27Yp1UUWC0sW6AOLK9+///u+Tfv7lL3/JW2+9ddb5uXPnRrNYQkxQZNMFEW333XcfzzzzDPKrJ6YL6RoLIRKeBEIhRMKTQCiESHgSCIUQCU8CoRAi4UkgFEIkPAmEQoiEJ4FQCJHwJBAKIRKeBEIhRMKTJXZCiIQnLUIhRMKTQCiESHgSCIUQCe+KCYTPPPMM2dnZmEwm8vPz2bdvX6yLJISIE1dEIPz1r3/Ntm3beOSRRygqKuKqq65iw4YNtLS0xLpoQog4cEXMGufn57Ns2TKefvppAMLhMIFAgPvvv59vfetbMS6dEGK6i/ut+kdGRigsLGT79u0T5zQaDevWrWPPnj3n/Mzw8DDDw8MTP4fDYTo6OvB6vSiKctnLLEQ0qapKb28vaWlpaDSXrxM4NDTEyMhIRK5lMBgwmUwRudaFiPtA2NbWRigUwu/3Tzrv9/upqKg452d27NjBo48+Go3iCTFt1NXVkZGRcVmuPTQ0RE6WjaaWUESul5KSQnV1ddSCYdwHwkuxfft2tm3bNvFzd3c3mZmZrGITOvQxLJkQkRdklPf5I3a7/bJ9x8jICE0tIU4WZuOwT63V2dMbJmtJDSMjIxcUCL/zne+c1bCZPXv2eRtC5xL3gdDn86HVamlubp50vrm5mZSUlHN+xmg0YjQazzqvQ49OkUAorjCnZwGiMexjsyvY7FP7njAX//n58+fz9ttvT/ys011caIv7WWODwcCSJUvYuXPnxLlwOMzOnTtZvnx5DEsmROIJqeGIHBdLp9ORkpIycfh8vov7/EV/4zS0bds2tmzZwtKlS7n22mt58skn6e/v55577ol10eKOotOhcTnpvy6PuvUKGbNbyHa0E1Y1NA/aaeqx09duQdeux9SuYGlScR4bQPmgJNZFnz40WrQ2K8H5OdTeYsVwVSczPG0YtCE6hyw09trp7bSgbddjbNdgblFxVo+ge6cw1iWfsjAqYaaWiHIpn6+qqiItLQ2TycTy5cvZsWMHmZmZF/z5KyIQ/vVf/zWtra08/PDDNDU1cfXVV/PGG2+cNYEihIgfPT09k34+35BWfn4+L7zwArNnz6axsZFHH32U66+/nkOHDl3wuOgVEQhh7KHh9913X6yLIURCCxPm4ju2Z18DIBAITDr/yCOP8J3vfOes92/cuHHiz4sWLSI/P5+srCx+85vf8OUvf/mCvvOKCYRiihQFxWBAE0ij8ZZURm7q5kszC1lnO4xHO0RTyErJYBZ7TLmUkUrfiI1wtw7tiIquqoHIJE3EP0VvQOtPou2mTFrXjnDzvINsdJeRpuukK2ShfDidA+ZsykilZ9iO2q1BOwKm6jaCsS58BIRUldAU12iMf76urg6HwzFx/lytwXNxuVzMmjWLY8eOXfB3xv1kiRDiyuRwOCYdFxoI+/r6OH78OKmpqRf8XRIIhRARMz5ZMtXjYvzd3/0du3fvpqamhg8++IA777wTrVbL5z73uQu+hnSNhRARE0YlFOVZ41OnTvG5z32O9vZ2kpKSWLVqFQUFBSQlJV3wNSQQig/TPRbkcGKzhdz8k9yd/gH5pjqsGoXWkIa6US9HB1I40e2lr82KqUmHoyaM+39qCba2xroGsacoaGw21DnZnLjVjue6Jr4eKOI683GcmmE6wibqRr0c6U+lsjOJnlYbpkYdjuow3j83Eqw+GesaxK2XX355yteQQCiEiJhY5RFOlQRCIUTERHLWOJpkskQIkfCkRZjgFJ0Orc9L5405NG4IsnnRAT7t2c9MXR8ANUEDRYPZfNCdR0lzOj2Ndiy1OjxHQ9jePESwvz/GNYg9RW9Am+yjfXUmzWuDrJtfwq2eEnL1bQCcCHo4NBhgf1cWFa1+BhpsWOu0eCqCWN86RHBgIMY1iJzw6WOq14g2CYRCiIgJRWDWeKqfvxQSCIUQERNSx46pXiPaJBAmKkVBYzSi5ASovTUJ/Q3tPJC3l5utR/Brw/SGoWrUzb6BPAo6cqhsSiJcb8F5UoOvbAjtu0Ux6cJMK2emzGyy47qumfsDB8m3HMOjGaI7bKRiJJWiviyK2zM41exGV2/EXQO+kj4oKJV7OE1IIBRCRIyMEQohEl4YhdAl7DD90WtEm6TPCCESnrQIE5FGi9ZhY3jxDE7eqmfRsiq+mLqHxcYmTIpCc0hD2XAaH/TOoLAtQEOjG8MpA+4TkLS3jVB5ZaxrEHMfTZlZM6+Uzd7i86fMNNown9Lirgzh2FVFqL0jxjW4PMLq2DHVa0SbBEIhRMSEItA1nurnL4V0jYUQCU9ahEKIiInXFqEEwgQzvpV869pMOtYP8al5BdzpLCRXP8SIClWjZooGc9jTlUtZSyr9DXasdVrcR0PY/nSIUKIvqfsLuYMuzTC9YcPH5g5eyY81CKsKYXWKs8ZT/PylkK6xECLhSYtQCBEx0jUWQiS8EBpCU+xoxmLoQAJholAUNGYzzMyi+jY37uub+IesAlZbqvBooDusUjHqY19/HgVtORxv9qGcMuM8qeArGUDzfnHCr4v9uNzBmqCbkoEsiroDk3IHPUdD2HdfubmDZ1IjMEaoyhihEEJEn7QIhRARI2OEYvrSaNG6nQzk51G7ScOKxeV8PnkPCwztGBSFU0EdpcPpFPTOoKg1g+YGF6ZTBpzHwng/aCR4oibWNYitjzyhzpU/OWWmK2ykaiQlYVNmzhRSNYTUKY4RxmCJnXSNhRAJT1qEQoiICaMQnmL7Sh7nKYSIazJGKKYlRW9Am5FK0/o0+tf2sWV2Ibc5ignoRhlRoWLUyoGBXD7ozOVIi5/Behu2Oi3e8lFMbx4kGAzGugoxNb4ksX11gOa1o6yb9+ET6sKqcs6UGUvd2HZb9ncrCXV2xroK4gJIIBRCRExkJkukayyEiGNjY4SyVb8QQsQdaRFeqRQFjcWCOjeHY7fbyVhezwMZH7DSXINTo9ARhvKRFPb0zeBAeyY1jV609Sbc1eA72Av7ymIwdzeNnM4dDM8d227Lc10T92cUn5U7eKA3h+L2dBqaXegajGP3rzSxcgfPFI7AWmOZNRZCxDUZIxRCJLwwGskjFNODotOhcbvpW5lD3S2wfnExn/HuY56+G4CTQT2FQ9ns7c6luDWNjgYn5jo9rmNh3LtrCDY2xbgGsXWhKTMHujI52prMYKMNyylJmYlnEgiFEBETUhVCU9xGa6qfvxQSCIUQEROZjVmj3zWW9BkhRMKTFuGVRFFQDAY0gTQab0lldE03X5lxgNscJaRpQwyoUDXqZG//DPZ1ZlPRnMxovRVHrQbv4RH0bxUSjMGM3bQxnjIzP4fjG214889OmakYTuNgX9aHKTPj220Vj6UcJWLKzJnCqobwFGeNwzJrLISIZ9I1FkKIOCUtQiFExISZ+qxvLB4SJoHwSqHRorVZCS7I4cRmCzOuq+GLaXvIN9Vh1yi0hhTKRtLY25vH/rYsTjW50Z8y4j4BSQc6CZcciXUNYkrRG9CmJNN2U4DWNSOsnTs5d/D4qJfSwUyKugIf5g5ObLdVQairO9ZVmBYik1Ad/Y6qdI2FEAlPWoRCiIiJzFrj6LfPJBBeARSdDq3PS+eNOTTeMsrmhQf4tGc/M3V9ABwbNZ1eUpdDSXM6PY12LLU63JUhHO8k+JKwM3aZOX6bDd+1zXw9UMQy84lJKTOFfdmUtqfR0ORG12DAUz2WMqPul5SZM8XrfoQSCIUQEROvLUIZIxRCJDxpEQohIiYyCdUyRiguhqKgMRpRcgLU3pqE4cY2Hsjdx83WI/i1YXrDUDXqZk//TPZ2ZHOsKYlwvRnnSQ2+siG0u4oSenzro9ttrZ9fzEZ32TlTZipa/Aw1WrGc0uI5GsT2TgWhnp5YV2HaCasK4anmEcZg9xnpGgshEp60CIUQEROZZ5ZI11gIEccis/uMBEJxoTRatA4bw4tnUHObnmuWVfL5lAIWG5swKQrNIQ1lw2l80DuD/a2ZNDW6MZ7S4zoOvoIWQkePxboGsaPRorFazrvdVnvYTMVwGsV9mZS0p9HY5EZfP5472IN64FBCj61eiSQQCiEiJoRCaIoJ0VP9/KWQQCiEiJh47RrLrLEQIuFJizAOjee/ta7NpGP9EJ+eV8CdzkJy9UOMqFA1aqZoMIc9XbmUtaTSX2/HWqfFUxHE+mYpoaGhWFchZj5uu62q0SRKBjI52BWgsjWJoQYr1lNa3EeDWN8+TLi/P9ZVmNZCTL1rG4vxVwmEQoiIideusQRCIUTExOumC9M6EO7YsYP//M//pKKiArPZzIoVK3jiiSeYPXv2xHuGhob45je/ycsvv8zw8DAbNmzgpz/9KX6/P4Ylv0wUBY3ZDDOzqL7Njfv6Jv4hq4DVlio8GugOq1SM+tjTN5P97Vkcb/ahnDLjqgZfST/KByUx2QZ9WhjfwXt+DtWbLCQtm7zd1njKTGFvFqXtaTQ1uTCMp8wUdaMePJy49y4BTOvJkt27d7N161YKCgp46623GB0dZf369fSfMU7zjW98g9///ve88sor7N69m4aGBj75yU/GsNRCJC719H6EUzlUSZ+Z7I033pj08wsvvEBycjKFhYXccMMNdHd384tf/IKXXnqJNWvWAPD8888zd+5cCgoKuO6662JRbCESVrx2jad1i/CjurvHHpDj8XgAKCwsZHR0lHXr1k28Z86cOWRmZrJnz57zXmd4eJienp5JhxAicU3rFuGZwuEwDz74ICtXrmTBggUANDU1YTAYcLlck97r9/tpamo677V27NjBo48+ejmLG3kaLVq3k4H8PGpv1bBq8SH+JmkvCwztGBSFU0EdpcPpFPTOoKg1g+YGF6ZTBpzHwnj/p57gybpY1yBmFL0BbZqf1psyaFszzM1zPtxua1TVTKTMFHZlUtWSxHCjBdvpdCPLzkOEBwZiXYW4EettuB5//HG2b9/OAw88wJNPPnnBn4ubFuHWrVs5dOgQL7/88pSvtX37drq7uyeOurrEDRJCRNL4xqxTPS7F/v37+fnPf86iRYsu+rNxEQjvu+8+/vCHP/Duu++SkZExcT4lJYWRkRG6uromvb+5uZmUlJTzXs9oNOJwOCYdQoj41dfXx1133cVzzz2H2+2+6M9P60Coqir33Xcfr776Ku+88w45OTmTXl+yZAl6vZ6dO3dOnDt69Ci1tbUsX7482sUVIuGNd42negBnjeMPDw+f93u3bt3KrbfeOmm+4GJM6zHCrVu38tJLL/G73/0Ou90+Me7ndDoxm804nU6+/OUvs23bNjweDw6Hg/vvv5/ly5dfOTPGioJiMKBNS6FxQxoDa/rYMruQ2xzFBHSjjKhQMWrlwEAuH3TmcqTFz2C9DVutFu+RUUxvHiQYDMa6FrExnju4MJfqTWZSlzXyjfSDLDZXY9eM0BU2Tc4dbHRjqNfjrQZfURfh4nLJHbxIYTRT3lh1/POBQGDS+UceeYTvfOc7Z73/5ZdfpqioiP3791/yd07rQPizn/0MgNWrV086//zzz3P33XcD8KMf/QiNRsOnPvWpSQnVQoj4VldXN2nYymg0nvM9DzzwAG+99RYmk+mSv2taB0JVVT/2PSaTiWeeeYZnnnkmCiUSQvwlIVUhNMVZ4/HPX8j4fWFhIS0tLSxevPjDz4dCvPfeezz99NMMDw+j1Wo/9jundSBMeIqCxmYjPDebY5ttZK44xYPpe1hprsGpUegIQ/lICnv6ZnCgPZOaRi/aehPuakgq6kXdX8bH/1NyBRofTkj107w2nc6bhtg4u4j1rjKy9R2MqhpOjPrOTpkZH054u5TwXxiPEucX7fSZtWvXUlZWNuncPffcw5w5c/g//+f/XFAQBAmEQogIUiOw+4x6EZ+32+0TecXjrFYrXq/3rPN/ybSeNRZCiGiQFqEQImKmwzNLdu3addGfkUA4TSk6HRqvh55VOdTfEmb9VcV8xruPefqx9dYng3oKh7LZ251LcWsaHQ1OzHV6XMfCuHfXEGw8/xLDK9rpp/uNXpXLiY0m0pc0sCWtmGXmE1g0o3SELFQMp1LYm01ZeyrNTa4PU2YOdBAurUjMcdUICatTWyI3fo1ok66xECLhSYtQCBExslW/ECLhjW+uOtVrRJsEwulGUdAYjSiZ6dRv8qOu7uRvZ+xjg+0wadoQAypUjTrZ2z+DfZ3ZVDQnM1pvxVGrwXt4BP1bhQQvIBH9inPGUsSmm9PoXj3IplkHWOc8PJE7WDWSTOlggMLOTI61+BhptGKr00juoJBAKISInEiuLIkmCYRCiIiRMUIxdad3Sxm9Ko+a20zMve4En0/dwzJjPXaNQmtIoWwkjb29eexvy+JUowd9vQH3CUg60Em45EisaxAbGi1ap4ORq3M4vtFI9pJTfCm1mCWmmomUmcPD6RT1ZFHWnkprkxNjvR7fCfAWScqMkEAohIigMBFYayyTJUKIeKZGYNZYHucphIhrsX5406WSQDhNKDod2iQf7WuyaVk/wp0L9vJJ9wFydWNPUDs2aqJwKJs9XXmUtaTS02jHUqvDXRnCsbOCUFd3jGsQA6dTZjSBNJpuTqHvhgFunbWfm52HCOi6GFK1VIz4KR3IpKgzwPEWH8EGC/Y6DZ4jo5jeOkg4UXfvFpNIIBRCRIzMGgshEl68do1l0wUhRMKTFmGsKQoasxklJ8DJzV4s17fyzZwC1liOkqRV6Q3D0VEve/vz2NuRzbGmJML1ZpwnNfjKhtDuKiIU6zrEgkaL1u1k6JocajfqmXFNLV9NPcgiYx0WzSitISvlQxkU9WZyqC2VtiYHxnoDvuMq3gNthMorEyN38PTjHuiNztfJWmMhRMKTrrEQQsQpaREKISImXluEEghj6fQa2aFledTcquPapRX8TXIBVxtbMCkKDSEtZcPp7OmZQWFbBk2Nboyn9LiOg29PM6HK47GuQfSdsU1ZwwY/Q9f3cueMItY6yknTdTOkaikfTuXgQBYlnemcaPYRajTjqNXgOTKC4a2DhMKJMaqq6A1o/Um0rEiCX0fnO+M1EErXWAiR8KRFKISImHhtEUogjBFFb0Cb5qd5XQbdawf5m3l7+YSziGzdCCMqVI2a2T+YS0FXLoeaUxlosGGr1eI+GsT6ZimhoaFYVyH6NFq0HhcDy3Kp26Bl7tXV3Ok/yEJTHQbCNAXtlA+nU9STyeG2FDoanZjq9XiOh/HsS6ChhNPbuQUX5lK90Yx7bm3UusYqU09/iUVakwRCIUTExGuLUMYIhRAJT1qEQoiIidcWoQTCaFMUNBYL6qxsTmx2krSyka9l7uF6y3E8GugOqxweSaagbwb7O7KobvKhnDLhqgFfST/KByWEY12HaBtPmckJUH9LEsFV3fx1bhmr7UdI0fUyENZzeDSF4v5MijszqGn2ojaYcNZq8JQPo99ZRCgRnux3xpP8mtel0XnjEBtmH2S5poy7o1SEeA2E0jUWQiQ8aREKISImXluEEgiFEBGjqgrqFAPZVD9/KSQQRpGi06FxOelfnkftrbD66jI+5ytgnqETg6JwKqijeDjAnp48iloCtDU6MdXpcR0L4/lzPcGTdbGuQtQpOh0at5v+63Ko26Bw9aIqPpFcwjxjPQbCNASdlA+lc6A7iyNtfroaHZjrdTiPh/HsqSdYfTLWVYiO08s1hxfncuwWAzmLT7ElpZSrzCdR+xIw5/QiSSAUQkSM7EcohEh4MkYozm88rSE9lYaNaYys6eYrMw+wyV5Kmi7IkArlI3b2D+awpyOXI80pDDdYsddq8JSPYnyziGCC7Jgy4fTO3eRmUrfRg7Kyky/klLLSWkmyto9+VU/FSColp1NmTjZ7UepNuGoVvIeH0O46SDBBUmbO3I1ncFUft884yI2OCtJ1XQypOo6M+oFTsS7ptCaBUAgRMTJZIoRIeNI1FkIkPGkRinPTaNFYLYQW5HJis4Xc607yhfQCVphOYtcotIYUykdSKejL40B7JicbvejqjXhOgK+oG7XwcKxrEHWKTofW56VnRTan1qssW3CU23wlzDE2okWlPuTk0GCAgz0Bylv99DbZMdfpcB0L4/5zLcFT9bGuQlRMpBbl51C3XmHBVSe4LbmU+cZT6JUQTSEH5UPp7GtPAQpjXdxpTQKhECJi1Ah0jaVFKISIayow1cn6WMz1y6YLQoiEJy3Cy2h8rKvrhhwaNgS5dVERn/bsZ7a+B4ATQQNFg9ns7cmlpCWNzgYn5jod7qowznePE2ptjXENouz0FmXMzKJmowvj8na+nF3CcmsVXs0APaqREyPJFPVlUdqRTl2zG83p3EFf2SCa94sTJ3dwPMdykweWd/E3OWWsslWSpO2lXzVQNTy2LVlJRzo1NZaoFS2MgiIrS4QQiSxeZ42layyESHjSIrwcxpc9ZWVQd1sy2hs6eGBGATdbj+DXhukNQ9Wom30DeezrzOZoUzLBeguOkxp8h4bRvZMgOyqfQdEb0Po8dN6YTePaECvnl3Obt4SZhmYAaoNuDg0GKOoOUNHqp7/JiqVOh7sqhON/qgk1t8S4BtGh6HRok3x0rcqiYe2HqUWzDM1olDB1QQ/lg+kUdQc42pZMX6MN4/Fg1MoXVhUUSagWQiQyVY3ArHEM2gDSNRZCJDxpEQohIiZeJ0skEEaaRovWYWPkmjxqbjWw6Noq7krZyzJTAyZFoTmkoWw4bWxJXWsm9Y1uDPUGPMchaV87ocNHY12D6FIUNDYb6pxsTmy047yuha8GSsi3HMelGaRHNVIxnEZpX4Di9nQaml3oGox4asBb2o+yp4SE2KBs/D7NzqJ6owNrfhtfzixlubUKl2aQ3rCJyuEUivsyKW1Pm3SfnAe7qYpSMSUQCiESXrxOlsgYoRAi4UmLUAgRMfE6ayyBMIIUnQ5tip+2NZm0rR/ik/P28inXAXL1Q4RUODFqYv9gLnu6cjnUmkJvgx1LnQ5PRQjb2+WEentjXYWoUvQGtCnJtK0O0LJmhDVzS9ngKSNP30oIhZqgl5KBTA52BahsTWKo0YrllBZ3ZQj7uxWEurpjXYWo+Oh9Wj23jA3uQ+QaWgirGmpGfZQOBijpyjjnfRrubItaWccC4VTHCCNUmIsgXWMhRMKTFqEQImJk1lgIkfBUpr6fYCwWl0ogjITT2yIpOQFqbvfiuL6Zv88q4EZLFUlale6wSsWoj339eextz+Z4UxJqvRnXSQVfySCa/zlIONZ1iCaNFq3NSnBhLtUbzSQtbeb+jGKWWk7gUIbpCpupGE7l4OmcuKYmF4YGA55q8BX3oB44lBi5g6dzUkcX5IzdpyXN/G1GCcvM1dg1Q7SHrFSNpFDcm0lZeyrNTS4MDfrT92nsMQ8JcZ8iQAKhECJipGsshBBx2jeOq0D4+OOPs337dh544AGefPJJAIaGhvjmN7/Jyy+/zPDwMBs2bOCnP/0pfr8/OoXSaNE6HQwty6PmVh0rlpXzueQCrjK0YVAUGkJaiocyKOidQVFrBs0NLkynDDiOq3j3NBI6Vh2dck4Tit6ANs1P87oMOm8aYsPsg6x3HSKg6yCEQuVoMqUDmRR3ZXCsxcdIoxVbnQZPRRDLO4cIDwzEugqXn6KgGAxoU/00rU+n+8Yh1s86yFpnOQF9O6OqlqoRP6UDAYq7Mjje4iPYaMFep8FTMYrp7VLU4eHYlD0CLUJkZcn57d+/n5///OcsWrRo0vlvfOMb/P73v+eVV15h9+7dNDQ08MlPfjJGpRRCxKO4CIR9fX3cddddPPfcc7jd7onz3d3d/OIXv+CHP/wha9asYcmSJTz//PN88MEHFBQUxLDEQiSm8ZUlUz2iLS4C4datW7n11ltZt27dpPOFhYWMjo5OOj9nzhwyMzPZs2fPea83PDxMT0/PpEMIMXXjkyVTPS7Gz372MxYtWoTD4cDhcLB8+XJef/31i7rGtB8jfPnllykqKmL//v1nvdbU1ITBYMDlck067/f7aWpqOu81d+zYwaOPPjrlsp053tWzboC/mbuXO52FBHSjjKhQMWrlwEAuBV05HG5OYbDBhq1Wi+dIEPObxYRGR6ZchrgxngpyVS4nbjERWFbPltQSFpursSqjtIRsHB1Oo7g3QFl7Kq1NToz1enwnwFvUQbi0IjFSjE6POQ9fk8vxDQZyl9bxpdQSFprqsCojtIbsVAynUtSTyeG2FDqanJjq9fhOqHgPtBEqr4xJHl4sZWRk8PjjjzNz5kxUVeXFF1/k9ttv5+DBg8yfP/+CrjGtW4R1dXU88MAD/OpXv8JkMkXsutu3b6e7u3viqKuri9i1hUhoqhKZ4yJs3ryZTZs2MXPmTGbNmsX3vvc9bDbbRQ2PTesWYWFhIS0tLSxevHjiXCgU4r333uPpp5/mzTffZGRkhK6urkmtwubmZlJSUs57XaPRiNFovJxFFyIhRXL3mY8OWV3I39tQKMQrr7xCf38/y5cvv+DvnNYtwrVr11JWVkZxcfHEsXTpUu66666JP+v1enbu3DnxmaNHj1JbW3tRN0EIMf0EAgGcTufEsWPHjvO+t6ysDJvNhtFo5Ktf/Sqvvvoq8+bNu+DvmtYtQrvdzoIFCyads1qteL3eifNf/vKX2bZtGx6PB4fDwf3338/y5cu57rrrLl/BFAWNxTK2vfxtTvwrG/haoIDrLcfxaKAjDOUjKezpm8GB9kxqmrxoTplw1UBScR8UlCbOOM54Tlx6Kk03p9J94yC3zCxivauMdF0Xo6qWipEUDg5kUdKZzokWL6FGC45aDe6KUUxvHSQcjN7jKGNm/BGwmek0rvczcH0ft84oYY2jnHRdF0OqjsPDaZQNBCjuzKCm2YvaaMJRq8F7ZAT92wcJhafBgroIJlTX1dXhcDgmTv+l1uDs2bMpLi6mu7ub3/72t2zZsoXdu3dfcDCc1oHwQvzoRz9Co9HwqU99alJCtRAi+iK5xG58FvhCGAwGZsyYAcCSJUvYv38/P/7xj/n5z39+QZ+Pu0C4a9euST+bTCaeeeYZnnnmmdgUSAgx7YTDYYYvYnVN3AXCWFN0OjQuJ30r86jbCGuuLuWvfXuZZ+jEoCicDOopHsqkoCePopZ0OhqdmOv0OI+F8fxPHcFT9bGuQvRotGjdToaW5FK7XkfeNWOpIFebTmJSgrSG7JQPpVPUk0l5m5/OJgemej3uEyrefS2Ejh5LiCGE8d+pgfxcatdrmXdVNZ/wlzDfeAq9EqIp5ODQYICDPQGOtPrpabJjrtfhPB7GU9BA8ERNrKswWZT/p23fvp2NGzeSmZlJb28vL730Ert27eLNN9+84GtIIBRCREwsdp9paWnhi1/8Io2NjTidThYtWsSbb77JzTfffMHXkEAohIicGOw+84tf/GKKXzjN02eEECIapEV4oc5IA2nYmMbImm6+MvMAm+ylpOmCDKlQPmJn70AeezpyqWhOZrTBir1Wg7d8FMObBwjGYjV5LIyngmRlUH9LMiOrerkzr5DVjiOkaHsmUkFK+jMp7UqnutkLDSactRo85cPodxYRSoR7dXpnc3IzqdvogRVdfDb7EDfYj5Kk7aVfNVAxlEZJf4DSjnTqWtxo6k24ahW8h4fQ7j44DX+nlNPHVK8RXRIIhRCRE6cbs0rXWAiR8C6pRdjf38/jjz/Ozp07aWlpIRyevC/IiRMnIlI4IUScidMW4SUFwv/1v/4Xu3fv5gtf+AKpqakoSvT79FGl0aKxWggtzOXEbRZyrzvJF9ILWGE6iV2j0BpSKB9JpaAvj/1tWdQ2etDVG3FXg+9AN+rBw7GuQdQoOh0at5v+/BzqNigsWHSC25JLWWSsm8iJKx9Kp7A7Kz5y4i4TRadD6/PSvSqb+nUqyxYcZZOvjDmGRjRKmPqgi0ODAYq6AxxtS6av0YbllA7XsTCuP58kWN8Q6yqc2yXsHnPOa0TZJQXC119/nf/+7/9m5cqVkS6PEEJE3SUFQrfbjcfjiXRZhBBxLpLbcEXTJQXCf/zHf+Thhx/mxRdfxGKxRLpM08p4F6brhhwaNgS5dVERn/bsZ7Z+bK+0E0EDxUNZ7OnOo6Qljc4GJ+bTXRj3rmqCTc0xrkGUnN6Rh7wAdbd4UFZ08jfZh1hlq/wwFWQ4leL+TEo60qlr8ZxOmVHwHh5G927RNEwFuQwUBY3Nhjo7i+pbHJjz27g78xDLrVV4tf30hk1UDqdQ2hegtCONU81utA1GXLUKvrJBNO8XT+/7lEhjhP/yL//C8ePH8fv9ZGdno9frJ71eVFQUkcIJIUQ0XFIgvOOOOyJcDCHEFSGRJkseeeSRSJdDCHEFUNSxY6rXiLYprSwpLCzkyJEjAMyfP59rrrkmIoWaFsaXieUEqNuUhPaGDu7P28cG22H82jADKlSNOjkwkMuejlyOtiQzfHpJnadiFOPrRQSnw47Bl5uioOj0aH0euq7PpmFN+JypIEeG0idtI2Vq0OGoVvEWNBGqSoC80/ElmslJtK0O0LJ6lOvnHWKD5zB5+hYAakZ9HBlKo6grQGVbMv3NVsyndDhPhPG8X0ew7lSMK3EBEmmMsKWlhc9+9rPs2rVr4qFJXV1d3HTTTbz88sskJSVFsoxCCHFZXdISu/vvv5/e3l4OHz5MR0cHHR0dHDp0iJ6eHr7+9a9HuoxCiHgRg8d5RsIltQjfeOMN3n77bebOnTtxbt68eTzzzDOsX78+YoUTQsSZROoah8Phs1JmAPR6/VnrjuOOoqBotWjsdkauzuHkRiOzr63mi2kfsMxYj0ujoSOscngkmYK+GezvyKK6yYdSb8JeP5YTp3+7MNa1uPzG75PNSmhOFtXrbBiXdfC5zHJW2sdyBwfCRmpGfBT3Z1LWmUZti4dwswlbgwbXsRC2Px0i1N8f65pcXuP3yelgdEEW1etMuJe0cnd6OfnW43i0ffSGTRwf8VPaF+BQZyp1LW6UJiP2Bg3uyiDmN4oIJsKT/GLokrrGa9as4YEHHqCh4cP1jvX19XzjG99g7dq1ESucECLOqBE6ouySAuHTTz9NT08P2dnZ5OXlkZeXR05ODj09PfzkJz+JdBmFEPEiTgPhJXWNA4EARUVFvP3221RUVAAwd+5c1q1bF9HCCSFENFxyHqGiKNx8880X9aSoaU+jRWM2o0n20b4qjeYbQ6xeeIi/9u3jKkM7Fo2WtlCIwyPJ7OvPo7Ajc2yb+UYjlkYFZ00Qw66S+H0E5fm2U1M0p/+jgKJB0WpQzGZI8dGW76M1P8SsWbWs8x9hnqkeu2aI1pCdE8N+SvsyKO/009TqRGkyYW1UcJwMYd95JH7HBy/mPlktqKnJNK9w07EsyMJZ1azxVTDH2IhFM0xT0Mmx4RTK+tI50uGnpcWJtskw8ftkfrMYNZ7GB6/0lSVPPfUU9957LyaTiaeeeuovvldSaIRITFf8ypIf/ehH3HXXXZhMJn70ox+d932KokggFELElQsOhNXV1ef885VEl+wjnJNB02IbXYuCXDO7ho2eUq4ytJOqszGqhoBBRlUdfSEjg0E94VENuqCCZhS0w+qFd2MuZldv5fxzWopG+eiJyT+f8fqkncQ1msnlOP3zxHtOd+/G/quAVoui04HRQMhppT/TSleejr7cEP6sDha4GgjoO7AqI4yqWrpCVtqCNtqHLfQOGQkP6DAMgW5QRd8fJtTbe4F1n/p9Ouseneu9l+M+eWz0ZlnpztXSlxskPdDOAkcDafpOTMro6ftkoW3URvuwlb4hI+qQFt2QMnafekOooyMXUO+Pu0dK9CYg4jSP8JJmjR977DEGBgbOOj84OMhjjz025UIJIUQ0XVIgfPTRR+nr6zvr/MDAAI8++uiUCyWEiE8KH44TXvIRg3JfUiBUVfWcD2wqKSmRLfyFEHHnotJn3G43iqKgKAqzZs2aFAxDoRB9fX189atfjXghoyXY0oahd4i05iQMfWkc1GfhNQ7g8g4AnVg0WgD0ShCLZgSjLohGpxLWQlgHYd3pcaIL2Ur9YrZbV8+/nZd6ESsaP/YbLzQtpEGDrcaMrSqJ1uu8NOs8HDYNkmLsxqUdwKIZxq4ZxK3rx20YxGococ8cImzUETIpBM0aTBYL4QtJn4nAfbqYewSRvU+Oagv2qrH0mXqDh3JzKqmGblyaAUyaURzaIdz6flyGASxGO/1GEyGTduw+WbToNFr4uO3cPu4eRXNr/ys9fQbgySefRFVVvvSlL/Hoo4/idDonXjMYDGRnZ7N8+fKIF1IIESfidLLkogLhli1bAMjJyWHFihXn3HhBCCHizQUHwp6eHhwOBwDXXHMNg4ODDA4OnvO94+8TQiSYK71F6Ha7aWxsJDk5GZfLdc7JkvFJlFAoTreoD4cIDw6i1tXj/WMfluYc3ulZSOdSM59PKWCxsQmnRmGOoZVRq45RVUswrKE27GEgZEIT1KG/8Rq0u+L0KX7nG0s6PfY2PtamjgLDwyi9vSQ1tuKqzKLm5ix+s8TKQMBIvvU4KbpurJph7JohrLphynRp1Gk99GuNqBotmjXzMb9RfGF5ctPNxd6nnj5SGuy4q7IpvzmXtqutbEozsdRygjRdJxZl7D7ZdCOUatNo1Dnp15hQNTo86xdjePNAbB72ewmu+JUl77zzzsSM8LvvvnvZCiSEENF2wYHwxhtvPOefhRBiwpXeNT7TG2+8gc1mY9WqVQA888wzPPfccxPb9bvd7ogWMqrUsWVyofYOTO8PMaspkyOts3jieidbsgq40VJFklZlsbEBq2YYm3aYvbpsjul99BgshHUmvLol6HcWxU135pKM36euLrT7B8irS6bteAYvrF7B8bm+iaezzTfW49IOkGzopcgYoNKcTK/ZStCkw2VdjPv9OoKn6mNdm8tHVVFHRwi1d2D4cz8zq1NoXpvKL250cnKWh7WucrL1bVxtqsWj68Oj76fUnE6VKYlei4WQSY/Ldi2OXccItbXHujYfL04D4SUlVP/93/89PT09AJSVlbFt2zY2bdpEdXU127Zti2gBhRDicrukFmF1dTXz5s0D4D/+4z/YvHkz3//+9ykqKmLTpk0RLaAQIn7E62TJJbUIDQbDxKYLb7/99sST6zwez0RLUQiRgBLpcZ6rVq1i27ZtrFy5kn379vHrX/8agMrKSjIyMiJawJhRVcIDAyjlx8nuSKalKcC/rLuZ2rle7nAVkq0bYYGhHbtmEI+unz2GXMqMqfQa7IRMRtz2a7H96dCFLSOLZ6qKOjxMsK4B72vduCqz2X/LAiqvTeYTGWXkW4+RrW/Dq+0jWd9DirmXUksaDVYXbTYjQ65MfKU+lA9KYl2Ty0tVCQ8NEa6pI/mVTtwVuby78WqOLEnhjvQSFptryNO34rINkGrowm/KpMyaSrPFxahNz5BrFr7iXtQDh2Jdk78skcYIn376aXQ6Hb/97W/52c9+Rnp6OgCvv/46t9xyS0QLKIQQl9sltQgzMzP5wx/+cNb5v7RztRDiyhevY4SX/PCmUCjEa6+9xpEjRwCYP38+n/jEJ9BqtRErnBAizsRp1/iSAuGxY8fYtGkT9fX1zJ49G4AdO3YQCAT47//+b/Ly8iJayJg5nQMWrG8k6XcD2BryeHnTcmqWeLnLv4erDG3M1A9i1xzFqe3HYxjggDFAs8lFyGRgxLYI754mQseuzEcbTBIOEerpQSmsIPdUEu1VAZ5bcz0V8/xs8pSRa2jhatNJknU9+A09HDQHqLQm0WW3MuKy4k7Jx767ilB7R6xrcnmFQ4S6e9DuK2dmnZ+WNek8vdrN2tlHudl9iGx9G0tMNSTpevAbeyi2ZHDM5qPDbmXY7cCTei2Wdw5f+WPPUXZJY4Rf//rXycvLo66ujqKiIoqKiqitrSUnJ0ce3CREIpvq7tTx9ID33bt3U1BQMGk3aq/Xy+OPP87KlSsjVjghRJxJpK6x0Wik9xxPIevr68NgMEy5UNNOOESooxPTe4eZ3ZhFcfNcTl7vZktWAastVaRpQ1jNNXh1fbj0A+w1ZlNt9tFhNjFiTyXJ70T5c3GsaxEVY0MJDbj/q2cslWbjAiry/dwRKCXfeoyArgOXbQC/vps0S4ASWxqNdhcjDiMDvtn4SvuhoDTW1bi8xlOOautJ+m0X7ooc3tu4iPJlfu7IKGGZuZo8fSteex+phi6KTFkctqXQbHcx4tRjT16Er7gHtfBwrGtyxbikrvFtt93Gvffey969e1FVFVVVKSgo4Ktf/Sqf+MQnIl1GIUS8UCN0RNklBcKnnnqKGTNmsGLFCkwmEyaTiZUrVzJjxgx+/OMfR7qMQog4MeUn2EUg/eZSXFTXOBwO80//9E/813/9FyMjI9xxxx1s2bIFRVGYO3cuM2bMuFzlFEKIy+aiAuH3vvc9vvOd77Bu3TrMZjN//OMfcTqd/Nu//dvlKt/0ccaSu5wuP81NGfzT2vXUzvPyCWcR2boRrjK0YXcM4dH1U2DM4bAphW6LjaDVgsd9LeY343RH5oulqoR7e1GKj5Lb4KOjKpPn1txA+fxUNnrLmGloYrG5hhR990QqTZXNR6fDxrDHhjstgVJpenrQ7D/CjLpk2qoC/GzNairmpHKz+zC5hpaJVJpUUzfFlgyO2310OCwMe5y4M67F+s4RwucYrxcX56K6xr/85S/56U9/yptvvslrr73G73//e371q18RDl/k8xKFEFemRBgjrK2tnbTN1rp161AUhYaGhogXTAgRf+J1jPCiAmEwGMRkMk06p9frGR0djWihhBAimi5qjFBVVe6++26MRuPEuaGhIb761a9itVonzv3nf/5n5Eo4zaijIwTrGkh+rR97fR4v3bqCE4t9E0vu5uj7cViP4NH14TYMcNCUQavZSdBiwOlajPfPjQSrT8a6GlGhjo4QbGjE9YdenJWZFG2cx9Hrkrk9s5RV1kqyde0T204VWzIptabRYHcx4jQykDQbb+kAyp4rfHsumFjG6X1tLPfyzxsXcXhZCncGSsi3HJ+UU1hszuSQLYUmu4sRlwF70gJ8JX2o+8tiXY0PxeETKi7pAe9n+vznPx+xwggh4lwirCx5/vnnL1c5hBAiZi55Gy4hhPiohNuPMKGdufa4KZvi5rnUrnLzxdNrjzN0QSyaajzaPrz6fvabsqg2+2i3mBixp+ErdiXE2BfwYU5haRU5DR46j2XzizXXU7kw+aztudKMXRRZAhy1JY/lFLqtuNLyceyOk0dZTsUZ25jlnUqivTLAs2tupGJ+ChvchyZyClN03aSZujhoCXDM7qPDYWXIa8edkY/93QpCXd2xrUecdo0vaYmdEEJcSaRFKISImHjtGk/7FmF9fT2f//zn8Xq9mM1mFi5cyIEDByZeV1WVhx9+mNTUVMxmM+vWraOqquryF+z0kjv18DFy/r2OoVf8/HPxzfx/Xcs4GdTj0sBiYwu3OEvZ6D/M4qxaDDN66JoXpnGVlaHbrkXRX4Fblp2HOjpCsKkZ5+/LmP3zAfa/voAfHlvH270LGAgbyda1s8payS1Jh7kp8xhpea0MzB2iZamG5k/NQl1+VayrEBUfbmN2mNnPDlPwxkKePL6Wt3sX0K8ayNa3cYOtgk1JZVyfeZzUma0MzB2meZmG5s/Mg+sWxbgCETqibFoHws7OTlauXIler+f111+nvLycf/mXf8Htdk+85wc/+AFPPfUUzz77LHv37sVqtbJhwwaGhoZiWHIhRDyZ1oHwiSeeIBAI8Pzzz3PttdeSk5PD+vXrJ56JoqoqTz75JN/+9re5/fbbWbRoEb/85S9paGjgtddei23hhUhEUW4R7tixg2XLlmG320lOTuaOO+7g6NGjF13saR0I/+u//oulS5fyV3/1VyQnJ3PNNdfw3HPPTbxeXV1NU1MT69atmzjndDrJz89nz549573u8PAwPT09kw4hxNRFe63x7t272bp1KwUFBbz11luMjo6yfv16+i/y4VbTerLkxIkT/OxnP2Pbtm383//7f9m/fz9f//rXMRgMbNmyhaamJgD8fv+kz/n9/onXzmXHjh08+uijESnjmUvubA15vLB5Fcev8fE3SXtZYGg/95I7k4ugRY/TuRjvnxsI1tRGpCzTnqoS7u9HKTlKblPy2PZca6+nYv7ZT7pLMXZTZAlQaR9LpRnyWHEnSirNmduYNU5+IuAGz2FmGpq42lRLkq6XNGP32DZm9mmyjVmU02feeOONST+/8MILJCcnU1hYyA033HDB15nWgTAcDrN06VK+//3vA3DNNddw6NAhnn322XMu97tQ27dvZ9u2bRM/9/T0EAgEplxeIUTkfLSnZjQaJ+1zcC7d3WN5lGc+WO5CTOuucWpqKvPmzZt0bu7cudTWjrWgUlJSAGhubp70nubm5onXzsVoNOJwOCYdQogIiOAYYSAQwOl0Thw7duz4i18dDod58MEHWblyJQsWLLioYk/rFuHKlSvPGvisrKwkKysLgJycHFJSUti5cydXX301MPavyN69e/na174W7eIKkfAimUdYV1c3qZHyca3BrVu3cujQId5///2L/s5pHQi/8Y1vsGLFCr7//e/zmc98hn379vGv//qv/Ou//isAiqLw4IMP8t3vfpeZM2eSk5PDQw89RFpaGnfccUf0Cnp6yZ35vXJmt2Szf/MCald4uDvzA1aYq8+95M7io91sYsSRjq/YnThL7gA1GPxwe66qLA5sXEDFdWOP/FxhrSJb147H1keqvos0cxbFlnQa7C5GnUYGfbPwlg2gfHDl369Jj0Y9OnafjuSncHtmKcutVWTr2/Bq+86/jVnZAMqe0rjcFgu4qN7afffdxx/+8Afee+89MjIyLvq7pnUgXLZsGa+++irbt2/nscceIycnhyeffJK77rpr4j3/8A//QH9/P/feey9dXV2sWrWKN95446wNZIUQURDlyRJVVbn//vt59dVX2bVrFzk5OZf0ldM6EMLYM5Rvu+22876uKAqPPfYYjz32WBRLJYQ4l2gvsdu6dSsvvfQSv/vd77Db7RPZIk6nE7PZfBHfqapx2nCOnJ6eHpxOJ6u5HZ2in9K1FKMRbaqfpg3p9K/t43OzC7nNUUxAN8qIqnIiaOPAQC4fdOZypMXPYL0Na50Wb/kopjcPogaDEapVfFD0BrTJPjpuyKRpbZDr51eyyVPGTMPYBFh90EXpYCbF3RlUtPoZaLRhPqXFXRnC8d4JQq2tMa5BdIzfp87rM2lcE2L5/GPc6islT9+CRgnTEHRzaDCDoq4AlW3JH96nqhCW3UfZ2fL/6O7uvmwTg+N/h+be9320xqn1xkLDQxx5+v9eUHkVRTnn+eeff5677777gr9z2rcIhRBxJAZd40iQQCiEiBzZj1AIIeKTtAgjTB0eJnSqAf9/9NNfn8cLt66i5hovf+3bd9aSO68x7/RT7k4vuXMtTawld3zkSXfnSKXJ1HWOPelO30WKqZdSaxqn7G5GnUYGkmbgK01LnFSahkacv+/GcTSbko1zqbgumduzylhurSKg68BlHSBZ30OauWfiPo04jVicOfD/olNO5fQx1WtEmwRCIUTkxGnXWAKhECJiZIdqIYSIU9IivAzUYJBQeweWd4eY3ZRDwW0LqVnpYUtgD9ebT5xzyd0J8+kld/Z0kopdUFAa62pEz/i2UyWV5DT66Dg2eXuumYbm09tOffikuwqbny772JPu3Kn5iZFTOL6NWdlRshu9dFVl88LaVRxekMptvhJmGZrPfiKgPZkufRT/mkvXWAghiMu1zdI1FkIkPGkRXi6nn3KnlFWR255CY1Maj6/ZwMkzltwtNrbg0g7g0/fxgSGXI2Y/XRYbozYbXu+yhFtyd6GpNH5991mpNIO+GXjLEiSVJhgk2NyC4499OCozKd84m6PXJXNHTinLrccmUmn8+m7SzAGK8HEqSmWL18kSCYRCiMiJ0zFC6RoLIRKetAiFEBEjXWNxNlUdW3JXV0/Kb/vGltxtGlty9zlfAfMMnX9xyZ3LuRTPn+sJnqyLdU2i56OpNMczx57gNt/Pbd4S8vStE6k0KcZuii0ZVNj8dNpPP+kuNYGedNffj3K4iswmD72VOfxy3UoOL0zltqRS5hgaWWiqI0nXg3c0mfM/3DbS5UK6xkIIEY+kRSiEiBjpGgshRJx2jSUQRsHEkrt3hpjdmMOezQupWXH2kjuXdgCvvp+9xmxqzF7aLCaGHRkkFbsTa8kdZ+QU/r4XZ1U2+zcuoPK6ZO4MlJBvOU62vgOX9sPtuSaedOc4/aS70oGEeDKgGgwSamnF9kYfcyozObpxFkdXJLM5+xA32I+SrutCb+mLYoGIy0AoY4RCiIQnLUIhRMTIGKEQQsRp11gCYbScsfY4rz2FxsY0nli7gVOz97PRXnp67XETLs0Abn0/e405lJtS6LJYx9Yee5Zh+lNirT2eyCksPkpuYxLtxwM8u+ZGKuansNFTelZOYZElQKUtmU7H2PZcrrR8HLuqCLV3xLoml9f471bFMQLNbvoqc3j55uWUL0xlc3IJ2SRQHuolkkAohIgYRVVRpviIzal+/lJIIBRCRI50jcXH+siSu77GPH6x8Qaqr/bx1769zDN0Ms/Qi0tTTpKuF7dhkGJTOm0m59iSO1cCLrnjdCpNfQPu/+rBWZVNwcaFHMn386nMg2el0qSZs8ZSaWwuRpxGBr2z8ZX2J0T6kRoMEmprw/pWP3OqMjixMY8frvBxc3IxcDLWxZvWJBAKISJGZo2FECJOu8aSUC2ESHjSIoyB8SV31reHmN2Qy59PL7n7QkbBWUvuPIZ+DpgyqbGcseSuyAX7ymJdjegaT6U5eJTchrFUmp+uuYmKealsOk8qzVFbMp12G8MeG+60fOzvVhLq7Ix1TS6v8VSayhNktHYyUJHNazcsBX4Xla+XrrEQQsRp11gCoRAiYuK1RShjhEKIhCctwlgZH8sprSSvLYXGxnSeWHPuJXc+fe9HltzZ8fqWYXyjCMKhWNckqsa252rC/bsenMdy+PPGRVTkJ3NHxuTtuZL1PaSYxx752WB3MeIwMuCbg6+4LyHGV8dzCs07e5lVlRy9RXbSNRZCiNh0badKusZCiIQnLcJYmrTkboC+hpxzLrmza47g0fafY8ndMjz/c4pg3alY1yS6wqGxVJrCCvLqk2k7HuBna1ZTOTeFW9xl5OlbWWKqIVnXQ5qxiwPmLKpsvrFUGrcdd3o+9ncrCHV1x7oml5eqEh4aIlRVE9XvZKqbJsimC0KIeCazxkIIEaekRSiEiByZNRaXamL7pHeGmN2QPbHkbkvGB6w015ClG8VuOY5X1zdpyV27xcSIPYDvoDshUkI+amx7rka8r/XgrMrhvU2LOLLMzyczism3HCNX34ZH20eyvodiSyYl1jQa7W5GHQYGffPwFfeh7k+A+xbFFCslPHZM9RrRJl1jIUTCkxahECJypGsshEh08TprLIFwulBVwn19Y0vu2lNpaMrg8TW38DezDnCbvYQ0XfCcS+46zTZGbXY83mUY30y8JXeEQ4R6etAcOMKMBj+txzN4es1qjs7xs9FdRraufSKn0G/o4aA5QKUtiQ6HdSynMO1arO8cIdzbG+uaXBniNI9QxgiFEAlPWoRCiIiRrrGYuvEld7WnSH2lf2zJ3S0fLrlbeL4ld2Yno9bxJXd1BE/Vx7omUaeOjhCsa8D3n924qnLZtelqji7zc0d6McvMJ8ZSaex9pBq68JuzKLWm0WRzMeIwYPctwFfcg1p4ONbViH9xOlkiXWMhRMKTFqEQImKkayyEEHE6ayyBcBqaWHK3c4hZDdm8v3khtSvcfCF9z/mX3JlPL7lzZJJU5EmMpWMfFQ4R6u0dS6Wp99N8Ip2f3ORm42w/611lZOs7ztieq5tCSyZVtiQ6HBaGPU68KcswvV2KOjwc65qIKJNAKISIGOkaCyGEzBoLIUR8khbhdDW+5K6kkhltqTQ0Bnh8jZPPz9rPJnvpWUvu9hhyqTD7P1xy50vQJXenczGDtfUk/0c37mO5/GnjYo4s9fPJtIMsMdWQq2/DZRvAr++m0JRNmTWVZpuLUbseR9JivEWdhEsrYl2TuCRdYyGECKtjx1SvEWUSCIUQkSNjhEIIEZ+kRTidfWTtcW9jLs+dXnv8WV8B8/Tdk9Yee40DE2uPgxY9zgRee0w4RKi7B21BOTNrU2g6kc6PVnvYNOswG1xlBHRdLDNXk6TrJc3UNZZTaE2i3W5hyOPBk7EM058OogaDsa5JXFGIwBhhREpycSQQCiEiJ05XlkjXWAiR8KZ1IAyFQjz00EPk5ORgNpvJy8vjH//xH1HP+BdDVVUefvhhUlNTMZvNrFu3jqqqqhiWOvLGl9zZ3i5n1i+GeO+tRTxes4m3BrLpDatk6Ua53nKcWz0lrEmvJDunheFZg7RfpdB0aybKkvmxrkJsjKfSnDyF/5UKsp/T8Po7S3nq5Do+GMhjVNUy09DCTbYjbEw+xKqsE3hmdtA7b4SmfD0ddy1DO29WrGsRV8bTZ6Z6RNu0DoRPPPEEP/vZz3j66ac5cuQITzzxBD/4wQ/4yU9+MvGeH/zgBzz11FM8++yz7N27F6vVyoYNGxgaGophyYVIUGqEjovw3nvvsXnzZtLS0lAUhddee+2iiz2tA+EHH3zA7bffzq233kp2djaf/vSnWb9+Pfv27QPGWoNPPvkk3/72t7n99ttZtGgRv/zlL2loaLikmyGEiD/9/f1cddVVPPPMM5d8jWkdCFesWMHOnTuprKwEoKSkhPfff5+NGzcCUF1dTVNTE+vWrZv4jNPpJD8/nz179pz3usPDw/T09Ew6hBBTp6hqRI6LsXHjRr773e9y5513XnK5p/Ws8be+9S16enqYM2cOWq2WUCjE9773Pe666y4AmpqaAPD7/ZM+5/f7J147lx07dvDoo49evoJfDudYcvfEGgd3zTxw7iV3xlyOmFLoNFsJWpx4fEsx/KkwJjNyMRcOEerqQvfBYWbWpdF0c4Af3uhh88xD3Ow8TJqum+vMx0nS9eA39lBkCXDC5qXdbmHY48OT6cTwpwRcrngpwqePqV4DzmqgGI1GjEbjFC9+btO6Rfib3/yGX/3qV7z00ksUFRXx4osv8s///M+8+OKLU7ru9u3b6e7unjjq6uoiVGIhRKQEAgGcTufEsWPHjsv2XdO6Rfj3f//3fOtb3+Kzn/0sAAsXLuTkyZPs2LGDLVu2kJKSAkBzczOpqakTn2tububqq68+73Uv578sQiSyS+nanusaAHV1dTgcjonzl/Pv7LRuEQ4MDKDRTC6iVqslHB5rO+fk5JCSksLOnTsnXu/p6WHv3r0sX748qmUVQhDRWWOHwzHpuJyBcFq3CDdv3sz3vvc9MjMzmT9/PgcPHuSHP/whX/rSlwBQFIUHH3yQ7373u8ycOZOcnBweeugh0tLSuOOOO2Jb+MvhzCV3v+mnt+EvL7kbe9xnGh1GF6NWAy5XPu7d1QSbmmNdk+gbv3cnakn5TTdDVTn8buO1HFmcwh0pxVxtOskcQzNebR+phm6KLJmUWlPpsDsZcehxJF2Ld28Locrjsa7J9BanK0umdSD8yU9+wkMPPcTf/u3f0tLSQlpaGv/7f/9vHn744Yn3/MM//AP9/f3ce++9dHV1sWrVKt544w1MJlMMSy6EiJa+vj6OHTs28XN1dTXFxcV4PB4yMzMv6BrTOhDa7XaefPJJnnzyyfO+R1EUHnvsMR577LHoFUwIcU6x2Jj1wIED3HTTTRM/b9u2DYAtW7bwwgsvXNA1pnUgFOf24ZK7IWY15vDebYuoWe7h7owPuM50ctJT7nzGPvYas6m1eGizGBlx5OI74EM9eDjW1YiNcIhQRyfG9weYVZtOw4ZsfnSDl0/MKGOt/TDpuh7slmMk6XpINvZSbMmg2ual3WFmyOPHm+lCv7MoMdOQLkQMusarV6+etOz2UkzryRIhhIgGaREKISJGCY8dU71GtEkgFEJEjswai6gaX3J38Ch5bWk0NgX4/k23nHPJnVM7yD5jNhWmZDrNVkZtTrxJS9G/laBL7lSV8NAQHKshraOLgWO5/HbDdVQtTuKO5IPMNzYwz9iIR9tHqqGLA+Zsyq1+Om0ORpxGnP58vHuaCJ6oiXVNRIRIIBRCRE6cPrxJAqEQImIiucQummTWWAiR8KRFGM9OLxsL19SNPeWuIYfnbrmBk1d7+Ix3H/P03czW92CxH8Kn78VtzKPElEan2cmoxYDLnY9r1wlCzS2xrklshEOE2jswv9vP7JoAVbfO5J9WJvHJ3BJusFWQruvBpR0gRddNqimbYmsG1XYvbXYzQ540vIe96N6VnMJJZLJECJHwVKa+H6GMEQoh4lm8jhFKILwCqMEgoZZWbG8NMqshh12br6bmOi9fzNjDdaaT5OpGsJuP4dX24dH3U2jKHFtyZzUyYs/DV5iUuEvuTqfSKFXVBP69m97jOfzy5pUcWZTC7UnFzDecmUrTzQFTFkesfrrsY6k0ruR83HvqCZ6UzX3jmQRCIUTkqERgjDAiJbkoEgiFEJETp5Mlkj4jhEh40iK8UnxkyV1DUyY7bnLyhZn72GgvI00bYpmpAZe2H6++f/KSO3uCL7njjHHWN/uYczyT8k2zObbCxyezS1hlrSRd14PDUoVf302qKYsiS4A6u4dWh4lBbwDfIR+a94oT9v5NCANKBK4RZRIIhRARE6+zxtI1FkIkPGkRCiEiJ04nSyQQXknOWHKX9pt+ehpyePaW1dRd7eHTnv3M1vcwT9+N/TxL7tyufJy7E3jJnaoSHhhAKT9GVouXruPZ/GLd9VQsTOETvoPM0rew0FiPV9tHirGbQksW5VY/XXY7I04zLn8+rv+pIdjYFOuaxE6cBkLpGgshEp60CIUQkROnLUIJhEKIyJH0GTFdjOfE2f80yMyGXP50+2JO5Hv5Ytoe8k11E2uPXdr+D9cemz20Wo0MO/JIOuAjXFwe62rEjBoMEmxuwfHHPuwnsinaOI/K5Ul8KrOY5dYqArouXJYB/LpuUkzZFFsyqLO5aXGYGPTm4CtNQflzcayrEROSPiOEEHFKWoRCiMiRMUIxrZxecqcpqmBGewb1TdnsWO3k8zP2scF2mDRtiOtMY6kgXn0/BcYcKs1Jp5fcufAmJ/aSO1SVcH8/SmkVOU0+Ok5k8uzaG6mc72eTp5SZhhauNp3Cqxvbnmu/OYsKq58uu41hlwV3Sj6O3ccItbXHuibRFVZBmeLvTFi6xkIIEXXSIhRCRI50jYUQIgKBMAY7s0ogvJKNL7k7cZKMl/voasjhpxtu4uQi31lL7jy6PrzG3LOX3L17nFBra6xrEjPq6AjBhkZcv+/FcSybP9+6iIr8ZD6VUcwKSxXZuk5c1gH8+m5STDkUW9JpsLkYdRgZ9M7CV9wH+8piXQ3xMSQQCiEiR7rGQoiEF1aZctdWZo2FECL6pEWYANRgkGBLG843BrE05PLmbYupuc7DF9M+YJmxfmzJnaUKj27scZ/7jVmcMrtpsxgZdswg6UBSQi+5Q1UJ9/aOPQahMZnWEwGeXrOaY/OSucVVSra+g6uNdXhOb89VZAlw1JZMp83GsNuOJ+1arDuPEO7tjXVNLj81PHZM9RpRJoFQCBE5MkYohEh4cTpGKIEwUYRDhHp70R2oYGZLOnXNOXzvRhdbZuxlna187Cl3xnpcmoGxJXemHI6akum0WBi1ufD6lqDfWZS4S+44nUpT34jv1R5cx3J5+9ZrOLLMz6fTi1hmPkGevh2vto9kfQ8p5l5KrGk02NyMOgzYvQvwFfegFh6OdTXEOUggFEJEjnSNhRAJTyUCgTAiJbkokj4jhEh40iJMJKpKeGgI5cRJMv6/PjobcvjJhjXULPLyac9+Zur6WGjoxKX9cMldsSmdLrODUasRt/tanO8m4NZSZzo91qo5cISZ9X6aq9P50U1uNs1KYYOrjGxdJ9eaqknW9eA39FBozqTKlkSHzcqw24kn9Vosb5cSHhqKdU0uD+kaCyESXjjMlB86Eo5+HqF0jYUQCU9ahEKIyJGusYgX40vuXK8PYq3P5fXNSzmR7zv/kjvTmUvuZpG0v5NwaUWsqxE7p7c3C9bWk/zbXlzHcnh941KOLvHzqbQilphqmKlvxWPvw6/vptCUTZk1lWabi1GHHofvGnz7OwgdPhrrmkRenAZC6RoLIRKetAiFEJEjS+xEXDlzyV1rBrXNOey40ckX8vadteTOrR8Ye8qdMYkuk4Wg1YPPsxjtrqJY1yK2wiFCXV3oPjjMzLo0Gjdk8i83ePjEzDJudhwiXdfDdeYTJOl68Rt7KDIHOGb10W6zMOzy4k298pYtqmoYdYq7x0z185dCAqEQInJUdeotOhkjFEKI6JMWoRAictQIjBFK+oyIqvEld8drCLzUe3rJ3U3ULPxwyd08QydWTTk27RAW3QiHdKn06GyE9Sa8pmUY3jxwRY1xXbTTqTShE7WkvtzN4PE8/vOWfI5e4+fT/gMsMtYzx9CMS9uPT9+Ly5DDIVMKvUYHIZMRt/1abH86RLi/P9Y1iYxwGJT426FausZCiIQnLUIhRORI11gIkejUcBh1il1jSZ8RMTGx5O6Pg1jr83j9trEld59PLWChsQGXZphsQxstVgedbgvVw3oGh010jurxjixG905hrKsQe+EQoY5OTLvKmH0yQM3GXP7peg+fzC3hBlsFds0QuYYWum0W+oJGKkd0DAStKGEtsADb6yVX7tZccUACoRAicqRrLIRIeGEVlPgLhDJrLIRIeNIiFGPG1x7vP8LM5nTqGnJ4eHkS81KbcRsH6B01Ut/npK3TTrjTgGFQQQnDiEuHKT2NYH1DrGsQe+N5mVXVZHR001udw7+vWcXOmbPJsHcRDGtoGbDT2mNjuNOEfkBBCUHQrMCcXCguj3UNpk5VmfIO1dI1FkLEMzWsok6xa6wmWtf4vffeY/PmzaSlpaEoCq+99tqk11VV5eGHHyY1NRWz2cy6deuoqqqa9J6Ojg7uuusuHA4HLpeLL3/5y/T19UWxFkKICWo4MsdFeuaZZ8jOzsZkMpGfn8++ffsu6vMxbRH29/dz1VVX8aUvfYlPfvKTZ73+gx/8gKeeeooXX3yRnJwcHnroITZs2EB5eTkmkwmAu+66i8bGRt566y1GR0e55557uPfee3nppZeiXZ34d7prR+VxUiqPw49hkLEDenHQhuMcHwtGtZDTnxoMEmpuwfJqCzNfHTvXefo1M21knuMz0c+cu3L8+te/Ztu2bTz77LPk5+fz5JNPsmHDBo4ePUpycvIFXUNRY9EOPQdFUXj11Ve54447gLHWYFpaGt/85jf5u7/7OwC6u7vx+/288MILfPazn+XIkSPMmzeP/fv3s3TpUgDeeOMNNm3axKlTp0hLS7ug7+7p6cHpdLKa29Ep+stSPyFiJaiOsovf0d3djcNxrn/Kpm7i75By55T/DgXVUXapr15wefPz81m2bBlPP/00AOFwmEAgwP3338+3vvWtC/rOaTtrXF1dTVNTE+vWrZs453Q6yc/PZ8+ePQDs2bMHl8s1EQQB1q1bh0ajYe/evVEvsxAJL8pd45GREQoLCyfFCY1Gw7p16ybixIWYtpMlTU1NAPj9/knn/X7/xGtNTU1nNX11Oh0ej2fiPecyPDzM8PDwxM/d3d0ABBmdci6oENNNkFEgOpMQkfg7NF7enp6eSeeNRiNGo3HSuba2NkKh0DnjREXFhT9gbNoGwstpx44dPProo2edf58/xqA0QkRHb28vTqfzslzbYDCQkpLC+02R+Ttks9kIBAKTzj3yyCN85zvficj1P2raBsKUlBQAmpubSU1NnTjf3NzM1VdfPfGelpaWSZ8LBoN0dHRMfP5ctm/fzrZt2yZ+7urqIisri9ra2sv2i3Il6OnpIRAIUFdXd9nGmuLddLxHqqrS29t7wWPml8JkMlFdXc3IyEhErqeqKoqiTDr30dYggM/nQ6vV0tzcPOl8c3PzX4wBHzVtA2FOTg4pKSns3LlzIvD19PSwd+9evva1rwGwfPlyurq6KCwsZMmSJQC88847hMNh8vPzz3vtczWxYWwMcrr88k5nDodD7tPHmG73KBr/wJtMpolsjmgxGAwsWbKEnTt3Tky0hsNhdu7cyX333XfB14lpIOzr6+PYsWMTP1dXV1NcXIzH4yEzM5MHH3yQ7373u8ycOXMifSYtLW2iwnPnzuWWW27hK1/5Cs8++yyjo6Pcd999fPazn72s//oJIaaPbdu2sWXLFpYuXcq1117Lk08+SX9/P/fcc8+FX0SNoXfffXd8q4pJx5YtW1RVVdVwOKw+9NBDqt/vV41Go7p27Vr16NGjk67R3t6ufu5zn1NtNpvqcDjUe+65R+3t7b2ocnR3d6uA2t3dHamqXZHkPn08uUex8ZOf/ETNzMxUDQaDeu2116oFBQUX9flpk0cYS8PDw+zYsYPt27efs8ssxsh9+nhyj+KTBEIhRMKbtgnVQggRLRIIhRAJTwKhECLhSSAUQiS8hA+EU93H7EqzY8cOli1bht1uJzk5mTvuuIOjR49Oes/Q0BBbt27F6/Vis9n41Kc+dVZmfyJ5/PHHURSFBx98cOKc3KP4ktCBcHwfs0ceeYSioiKuuuoqNmzYcNayvUSye/dutm7dSkFBwcQej+vXr6e/v3/iPd/4xjf4/e9/zyuvvMLu3btpaGg4536SiWD//v38/Oc/Z9GiRZPOyz2KM5chtzFuXHvtterWrVsnfg6FQmpaWpq6Y8eOGJZqemlpaVEBdffu3aqqqmpXV5eq1+vVV155ZeI9R44cUQF1z549sSpmTPT29qozZ85U33rrLfXGG29UH3jgAVVV5R7Fo4RtEUZqH7Mr3fgWZR6PB4DCwkJGR0cn3bc5c+aQmZmZcPdt69at3HrrrZPuBcg9ikfTdtOFyy1S+5hdycLhMA8++CArV65kwYIFwNgekAaDAZfLNem9Z+4TmQhefvllioqK2L9//1mvyT2KPwkbCMXH27p1K4cOHeL999+PdVGmlbq6Oh544AHeeuutqO+2Ii6PhO0aR2ofsyvVfffdxx/+8AfeffddMjIyJs6npKQwMjJCV1fXpPcn0n0rLCykpaWFxYsXo9Pp0Ol07N69m6eeegqdToff70/4exRvEjYQnrmP2bjxfcyWL18ew5LFlqqq3Hfffbz66qu888475OTkTHp9yZIl6PX6Sfft6NGj1NbWJsx9W7t2LWVlZRQXF08cS5cu5a677pr4c6Lfo7gT69maWHr55ZdVo9GovvDCC2p5ebl67733qi6XS21qaop10WLma1/7mup0OtVdu3apjY2NE8fAwMDEe7761a+qmZmZ6jvvvKMeOHBAXb58ubp8+fIYljr2zpw1VlW5R/EmoQOhqk59H7MrDefYHxJQn3/++Yn3DA4Oqn/7t3+rut1u1WKxqHfeeafa2NgYu0JPAx8NhHKP4otswyWESHgJO0YohBDjJBAKIRKeBEIhRMKTQCiESHgSCIUQCU8CoRAi4UkgFEIkPAmE4rK4++67ueOOO2JdDCEuiOw+Iy6aoih/8fVHHnmEH//4x0iuvogXEgjFRWtsbJz4869//WsefvjhSc81sdls2Gy2WBRNiEsiXWNx0VJSUiYOp9OJoiiTztlstrO6xqtXr+b+++/nwQcfxO124/f7ee655+jv7+eee+7BbrczY8YMXn/99UnfdejQITZu3IjNZsPv9/OFL3yBtra2KNdYXOkkEIqoefHFF/H5fOzbt4/777+fr33ta/zVX/0VK1asoKioiPXr1/OFL3yBgYEBALq6ulizZg3XXHMNBw4c4I033qC5uZnPfOYzMa6JuNJIIBRRc9VVV/Htb3+bmTNnsn37dkwmEz6fj6985SvMnDmThx9+mPb2dkpLSwF4+umnueaaa/j+97/PnDlzuOaaa/i3f/s33n33XSorK2NcG3ElkTFCETVnPvJSq9Xi9XpZuHDhxLnx58eMP061pKSEd99995zjjcePH2fWrFmXucQiUUggFFGj1+sn/awoyqRz47PR4XAYgL6+PjZv3swTTzxx1rVSU1MvY0lFopFAKKatxYsX8x//8R9kZ2ej08mvqrh8ZIxQTFtbt26lo6ODz33uc+zfv5/jx4/z5ptvcs899xAKhWJdPHEFkUAopq20tDT+/Oc/EwqFWL9+PQsXLuTBBx/E5XKh0civrogc2apfCJHw5J9VIUTCk0AohEh4EgiFEAlPAqEQIuFJIBRCJDwJhEKIhCeBUAiR8CQQCiESngRCIUTCk0AohEh4EgiFEAlPAqEQIuH9/+sfmFbXIofYAAAAAElFTkSuQmCC",
|
|
"text/plain": [
|
|
"<Figure size 640x480 with 2 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUIAAAHHCAYAAADK0FH3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAABqkUlEQVR4nO29eZAkV3Xv/7mZWfve+z7Ts2gBsQgEsjT8bAxjyWYxemD74Sc7hLCRgZFAyMZGLwAhYRjE84OxkCyMwpbkMLIwD4MBGxFYgAge2pAMRg+QWAQazdJ7d62ZVZWZvz+yqrqql5leqqu7ps4nIqOrs6uybt+pPnPPye/3XOW6rosgCEIHo233AARBELYbCYSCIHQ8EggFQeh4JBAKgtDxSCAUBKHjkUAoCELHI4FQEISORwKhIAgdjwRCQRA6HgmEgiB0PBIIhZZw1113oZSqHYZhMDw8zJve9CaOHTu23cMTOhxjuwcgdBY33XQT4+PjmKbJQw89xF133cW3v/1tnnjiCYLB4HYPT+hQJBAKLeW3fuu3uOCCCwD44z/+Y3p6erj55pv54he/yO/93u9t8+iETkVSY2Fb+f/+v/8PgJ/97GfbPBKhk5FAKGwrv/jFLwBIpVLbOxCho5HUWGgpCwsLTE9PY5omDz/8MDfeeCOBQIDXvOY12z00oYORQCi0lIMHDzZ8v3v3bv7xH/+RkZGRbRqRIEggFFrMbbfdxllnncXCwgJ///d/z7e+9S0CgcB2D0vocCQQCi3lpS99ae2u8WWXXcbLXvYy/sf/+B88+eSTRKPRbR6d0KnIzRJh29B1ncOHD3P8+HFuvfXW7R6O0MFIIBS2lZe//OW89KUv5ciRI5imud3DEToUCYTCtvPud7+biYkJ7rrrru0eitChSCAUtp3Xv/717N27l7/6q7/Ctu3tHo7QgSjZ11gQhE5HVoSCIHQ8EggFQeh4JBAKgtDxnDGB8LbbbmP37t0Eg0EuvPBCHnnkke0ekiAIbcIZEQg/85nPcN1113HDDTfw+OOP84IXvIBLL72UycnJ7R6aIAhtwBlx1/jCCy/kJS95Sc2d4DgOo6OjXHPNNbznPe/Z5tEJgrDTaXuvcbFY5LHHHuP666+vndM0jYMHD/Lggw+u+BrLsrAsq/a94zjMzs7S3d2NUmrLxywIrcR1XTKZDENDQ2ja1iWBpmlSLBabci2/39/SrRvaPhBOT09j2zb9/f0N5/v7+/nxj3+84msOHz7MjTfe2IrhCcKO4ejRo1vW7sw0TcZ3RTk52RxB/MDAAE8//XTLgmHbB8KNcP3113PdddfVvl9YWGBsbIyX8SoMfNs4MkFoPmVKfJt/JxaLbdl7FItFTk7a/PKx3cRjm1t1pjMOu178C4rFogTCtdLT04Ou60xMTDScn5iYYGBgYMXXBAKBFXvgGfgwlARC4QyjchegFWWfaEwRjW3ufRxaX55q+7vGfr+fF7/4xdx///21c47jcP/993PRRRdt48gEofOwXacpR6tp+xUhwHXXXccVV1zBBRdcUGvplMvluPLKK7d7aO2JUii/H70rhT3UTX4kQq5fx+xRFJMu5biNFilh+G18Pq8mVCzqlC0D19RRloZWVGhFhV5UaEXQSvWH630tg1Z2vcP2zivbO7Ty4mNlOyjHe0zZQTkOuJXvHQds73vv3OLj2vnK813XBaeyPHKrP6v7vvK4JqRwKn+Q9cIKZ/GPtEFw4SwXXyifgRaP4Qx0UxiKkBswKPQqiimXUtxGhW2MQBnDZ6OUS6lUN4dFDc06zRyWK4+r81d2V5xDxyrAd/61OZ+N0+Dg4rA5IcpmX78RzohA+N//+39namqK97///Zw8eZIXvvCF3HfffctuoAiCIKzEGREIAa6++mquvvrq7R6GIHQ0Dg6bTWw3f4X1c8YEQqFJaDpaJIzq78EaTZEeC5AbVhSGbPx9eUZSaUajc/QHMiSMAj5lk7GDTFhxpswos2aYjBkgb/opWQalJWleLWWupnhFVUnp6lNlammdVnZRlfRZOd5jL012aqlx/VfXqaS4tfS5Lm12qimzDa5XHq+ly3WpsgLve10Dx21MlXW9liqr6jkAXdWep5RChULQk8IaTpId8ZMd0SgMOGh9Bfq60ozG5hkIpmtzaDq+lefQNChZlXKDpdBKCr06l+Wlc6hWnEOn0Lo/c9t1sTfp0djs6zdC298sEQRB2CyyIhQEoWnIzRJBEDoeBxdbAqHQzijDQEsmcId6yY3FSY8Z5EZcSoMWvX1p9iWn2ROZZsQ/S1LPoymHjB3CdHxoyqkcLqrucBWgqH1FuaCU9z2Vc7D4PYvf104phavcBpmtqxRqpT8YpaD+vFKN8hcATavVEJehqRWlMMveY4XXKqXA50OLRrAHusmPRcgMG+SHXYoDRbr60ownZxmPzDAcmKNLz6Irl5wT4Fgxtcp71X2tm6uG+RJ7/KaRQCgIQtOQ1FgQhI5H7hoLgiC0KbIiFEAptEAArbeH0kg32bEQmTGN/JCDPlBgT+8cZycmGQ9NMeSbJ6YVAJi1o5woJZkoxpmxIsxZYTKWn4Llp1Q0cCxPQ6isOptYUdVsYqqMpwssexpCzxrmfdVqFjFQZQfNrugHHU/rp5WdZRrCBltdVTfo1Fnpqta7pZY7t65eWP2++rhKnbWOpRY8AE1DBQLLtIP5QQejP8+ennn2xafYFZxlyD9HTDOxUczbEY4VUxw3kyvqMGuWxQYNIXUaQk6rw3Tt1q2wnMqx2Wu0GgmEgiA0DbsJd403+/qNIIFQEISmYbvesdlrtBoJhJ2OpqNHIzDQS2EsSWbMT3ZEYQ6ViQ1k2Nc1zf7YFGOBGXqNNBHNouQaTJbjPGN1c8JMMGVFmTdDLBSCmJZv1ZROlaiz1p0mpSvVWeyq1rr6LjRL0+JTWOtq5+skMytZ67wfLKbFa+1Cg8+HFo9h96WWSWZSvRn2pGbYG51mLDBDt54lqJUwHR/HSymOWUkmrDjTZoR5M0TWDGCZPsqWDhVrXS0drisrVA+91smnfv5YtCbaLm6p7bcl2nIkEAqC0DSkRigIQsfjoLA3qfCWDtWCIAjbgKwIOxjl86N1JXGGesnvipAeNcgNu9hDJkN98+xNTHNWZLJmqfOpMjknwPFSiqNm19pqWyUqdcLltS1VrtS4KtKZWo2rUtuqtdwqOw1fsV2vS3X1a51sZsUO1atJZmBzXao1DRUKQldyzZIZgJlylGeLXRy3EhXJTIR0nWTGMY3F+mqtdZlXS9Vrnaqr87dYH6zWVJfKjtxS65JNx+W0DsW1XKPVSCAUBKFp2E1IjTf7+o0gqbEgCB2PrAgFQWga7boilEDYiSiFFgqh9fVQHOkiOxasWOps/AN59nTPcXZigt3BGYZ8c4Q1C/BqW88UezhuJZk0o8zUalsBipaBY+ooS0fVdl9bSfu2snZQq2kHl1jqypX2+bbT0J4fh0qdcHGnumZZ6latD25CO+hTNiVX3xLt4Gnrq3Yra4QKZ2lPtQ1co9VIaiwIQscjK0JBEJqGpMaCIHQ8Nhr2JhNNu0ljWQ8SCDsMZRhosRgM9JDflSQzapAdUVhDJVL9afakZjgrOlnxFmcIakVKrsHJUoJnrG4mrDgTZox5M0TGDFAwfZSr3uK6bTtrbbeW1gfL1RZSLppd9Ru7lRZcbs1r7B3N8Ra7S+uB26AdnLfDa9cOllRlHuu82evUDtZve+qWW1cjdJtQI3SlRigIgtB6ZEUoCELTkBqhsLNRCuX3o3elsId7yO6KkBnVyQ27uIMFRnvn2Z+YYm94ihH/DEk9j17Zpe5oqYtjVopJM8aMGWHODJEz/ViWD9s0wNLQamndouSj3g5W66i8tF1Utd2WXZV+VKxzdanxqpa6qmTGrUuV12qpW6tkpnpO00DXV5XMVHeoa5VkZqV2W6taEsut86zZrobtbrJGuA0WO0mNBUHoeGRFKAhC03BQOJtcX8l2noIgtDVSIxR2LpqOFgqi9fdSHE2RHguSHVUUhm2C/Tl2d89yVnyS3cFpBowFYnoB29WYKsf5pdVzasnHKXap04veLnUrWupOIfk47S51VclMMy11p9qlrsWSmY2226qXzNTPm7K3Q5nXXkggFAShaTTnZomkxoIgtDFejVBa9QuCILQdsiI8w1GGgZaI4w72kRuLkR4zyI1AcahIT1+avalpzo5OMOyfq2w1WcR0/BwvpXi2mGLSijFRiLFgBcmYAcyCv85Spzz9YIlFDVwZVm+5VbHVnc5SV3YatYPVGt5K7fjrW/IvrQ3C5ix19drBgRT5keXawb2pGcYjM63VDi7RXp52O4MWWuycJniN5a6xIAhtjdQIBUHoeBw00REKO4iqpa6nm/JoD9mxMJlRjfyQgxo02d07x1mJSfaEphn0zdUsdfN2hGeLXRyzkkxanqVuvhAib/mwTD+2qUNRQzO12u5q+lq7UJc9Kci6LXVLu1C3yFKnQkHoTmENJ8iMBMiOKAqDNr7+Anu6PcnMntA0/b6FFSUzE4U481aIhUKQguWjaPpwKynxVktmlkmOhFMigVAQhKZhuwp7k220Nvv6jSCBUBCEptGcxqytT41FPiMIQscjK8IzEU1Hi4RR/T1YoynSuyr1reEy0YEs412znB2bYHdwml4jTUSzsF2NY6UujhVTHDeTFUtYmEydJWxpF+pajbAJXZQ31IV6A5a6ZbVBaLTUadriDnVVycyIQX7IpThYpKt3ZcmM6fqYKCU4ZiU5YSaYMSMsWEGyZgDTrLQrK2qtk8wsaVfWKhxXw9nkXWNH7hoLgtDOSGosCILQpsiKUBCEpuGw+bu+2yH2kUB4hqEMAy2ZwB3qJTcW9yx1oy7lQYv+vgX2JafZH5lkxD9Lt5HFp8rknEDFEpZiwooxWbHUZWuWuiWWsGo7/uLy+pYq02ipq9W5lljqTlXfOpWlrqFO2ERLHYCuo4KBVbWD+3rm2BObXqYdnLGjHCumOGnFT68d3JbaauvacDVHUN36RFVSY0EQOh5ZEQqC0DSa4zVu/fpMAuGZglJogQBabw+lkW6yu0I1S50+mGdv7xxnxyfZE5qqpHUFAGbKUZ4tdnPcSjJpRZkxI6TNoLdLnenDqVrq6napq1nqTrVLnb26pU6VvdTVS+3q0zlWttQ1qws1nDIt1hJx7P4k+VGvy0xu2KU0sCiZ2ReZYjgwt27JTL0dcXkH7w1IZtayw199SaGV8pk27UcogVAQhKbRritCqREKgtDxyIpQEISm0RxBtdQIhY2g6ejRCAz0UhhLkhnzkx1RmENlYgMZ9nVNsz82xXhgil4jQ1ArUnINTpYSPFvs4oSZYMqKMlvwLHWm5aNkGstaRjV0oS7W1QRXaLmlbHexxlWrb23SUuc4jbVBaKwP1ktmKj87bW1QqZpkxhyp7FA3rCgM2fj6Ti+ZaRs7YotwXIWzWR3hNnSfkdRYEISOR1aEgiA0jebsWSI3SwRBaGOq3Wc2e6wV27Z53/vex/j4OKFQiL179/LBD35wsSSyRmRF2OYonx+tK4kz1Et+V4T0qKd/s4dMhvrm2ZuY5qyKpS6p52uWumeL3Svq3yzTR3mDLaO8+lalRljeYMuo9Vjqlrbjh3VZ6pRhoOKxRe3gSFU7aNHTm2E8ubJ28HgpyTErxaQZW1k7aDXqLpdqB1fUXa6iHVzTVgaraQfrtzU4Q7n55pu5/fbbufvuu3nuc5/Ld7/7Xa688koSiQTveMc71nwdCYSCIDQNG4W9SUH0el7/ne98h9e97nW8+tWvBmD37t380z/9E4888si63lNSY0EQmkYzU+N0Ot1wWJa17P0uvvhi7r//fp566ikAvv/97/Ptb3+b3/qt31rXuGVFKAjCjmR0dLTh+xtuuIEPfOADDefe8573kE6nOeecc9B1Hdu2+dCHPsTll1++rveSQNiuKIUWCqH19VAc6SKzK0h2VCM/ZOMfyLOne46zExPsDs4w5JsjphewXY2pcryy3WSSSbPqLQ6QNwMULQPH1FGWjrKU5ylesT64ynadNe1gnf6toh/EcdHKzurbTbbSW6yUt1VnV9JrtzUaqGkH/f159nXPsz8+xa7gTE07aKOaph1cVXe5FVsZ2HZLt/O0WV9qu9o1AI4ePUo8Hq+dDwQCy577z//8z3z605/mnnvu4bnPfS7f+973uPbaaxkaGuKKK65Y83tKIBQEoWk0Z88S7/XxeLwhEK7Eu9/9bt7znvfwxje+EYDnPe95/PKXv+Tw4cMSCAVB2B5a3XQhn8+jaY3P13UdZ513ynd0IDx8+DD/8i//wo9//GNCoRAXX3wxN998M2effXbtOaZp8qd/+qfce++9WJbFpZdeyt/8zd/Q39+/jSPfWpRhoMViMNhLfixBZtQgO6Kwhkqk+tPs65pmX2SKscBMg6XuWKUL9QkrwZQZZc4MeZY600e5mtpZWi2t01foQr1c+uGi1ck+GrpQ10k+Vk3t6tK5FWUza7XUrVcyM5AiPxL2JDNDXrutnr40e1PT7AlPn1YyM2eGvFZl1qklM/qpWpU1QzJTLScslcysNHdnIK997Wv50Ic+xNjYGM997nP5z//8Tz72sY/x5je/eV3X2dF3jR944AEOHTrEQw89xNe+9jVKpRKXXHIJuVyu9px3vetdfOlLX+Kzn/0sDzzwAMePH+f1r3/9No5aEDoXt9KPcDOHu44a4yc+8Ql+53d+h7e//e2ce+65/Nmf/Rl/8id/wgc/+MF1jXtHrwjvu+++hu/vuusu+vr6eOyxx/jVX/1VFhYW+Lu/+zvuueceXvGKVwBw5513cu655/LQQw/xK7/yK9sxbEHoWFqdGsdiMY4cOcKRI0c29Z47ekW4lIWFBQC6uroAeOyxxyiVShw8eLD2nHPOOYexsTEefPDBVa9jWdYyjZIgCJ3Ljl4R1uM4Dtdeey0HDhzgvPPOA+DkyZP4/X6SyWTDc/v7+zl58uSq1zp8+DA33njjVg53a1AK5fejd6Wwh3vI7oqQGdXJDbu4gwVGe+fZn5hib3iKEf8MST2PrhwydoijpS6OWSlOmnHPFlZtx7+eGteKso8W1bhWstRtRDLTncIcTjS02zqdZOYZq5sJK86EGWO+UldtkMzUtSrbMZIZt25OpA3XaWmbFeGhQ4d44oknuPfeezd9reuvv56FhYXacfTo0SaMUBCEamPWzR6tpi1WhFdffTVf/vKX+da3vsXIyEjt/MDAAMVikfn5+YZV4cTEBAMDA6teLxAIrCjOFAShM9nRK0LXdbn66qv5/Oc/z9e//nXGx8cbfv7iF78Yn8/H/fffXzv35JNP8swzz3DRRRe1eriC0PFUU+PNHq1mR68IDx06xD333MO//uu/EovFanW/RCJBKBQikUjwR3/0R1x33XV0dXURj8e55ppruOiii868O8aajhYKovX3Yo11ebawUUVh2CbYn2N39yxnxSfZG5yi10g3WOp+afVwvKIdnK1Z6rwal2MajTWuunb89fXB9baUP62lrlrjaoKlblltEBrqg8owUIk49kCK3EiYzIhOfmix3VZVOzgWmKm0KmvUDq61rqpbypunraqrLrUirrWu6tbNyxbjoG26sep2NGbd0YHw9ttvB+DlL395w/k777yTN73pTQB8/OMfR9M03vCGNzQIqgVBENbKjg6Ea1HEB4NBbrvtNm677bYWjEgQhFNhuwp7k6ntZl+/EXZ0IBQ8lGGgJeK4g33kdsW8LtQjUBxatIWdHZ1g2D9Hn5HGp8qYjp/jpRTPFlOcMONMm1EWrCDpQtDrQl3rlqLQLM1L54pr2aVujZa6stOY2lXTtpW6UK9mqau30G3EUleVzPR2YQ7FyYz4yVUkM4H+PONdc+yNTy+TzEyVY2fY7n4tTI3bVD4jgVAQhKbhNqH7jLvJ12+EHX3XWBAEoRXIilAQhKbR6j1LmoUEwp2MUmiBAFpPN+WRbrJjYTKjGrlhB23AZHfvHGclJtkTmmbQN1ez1M3bEZ4tdnHMSjJpxZgqRElXpB9Fy4dt6sulH6fqQl2q1AWrcpmqfKYV0g/YkKWuKpkpD1bbbXmSmfKARU9vmj2VHeoG/fO1dls5x8+JUqo2bw2SGdPvzVuxbt6Kp5m3+u7dO2HeWoCnnNpsjbBJg1kHkhoLgtDxyIpQEISm0cxW/a1EAqEgCE2j2lx1s9doNRIIdyqajh6NQH8P5miKzJi/YqkrE+3PMt41y9mxCXYHp+k10kQ0C9vVOFbqas+d1pa2jdqopU4pVDgMvSkKQ7Fauy1zyCbQl2dP9yz7YtPsCk3Ta2ROqx0sNGguF+dNX0lzWW3Hv1bt4Om2MdiUdnC55lJYHQmEgiA0DXGWCILQ8UiNUGgayudHSyZwh3rIj8VIjxnkRlzKgxb9fQvsS06zPzLJiH+WbiOLT5XJOQGOV3epMxNMmxEWrCBZM4BZ8FO2dKhYw/T6LtRLdqnTK51mGix1y1LjylFetNMtk36cylLXkOadRvoBa7PUKYXyLUpmciMhMsM6+UEXe9CkrzfNnsQMeyLTjPhna11m6iUzJ804c1aY+ULojJMaCadGAqEgCE3DoQleY7lZIghCO+M24a7xerbzbBYSCAVBaBrSfUbYPEqhhUJovd2UhrvI7AqRHdXIDzkYg3n29sxxdnySPaGpStuoAgAz5SjPFrs5biWZtKLMmJGapc4yfThVa5ipVaQflfrgStKPUqP0w6ttVWqEZQflgGY7qLJXi/JkM/WSD1auc22yCzWcoj6o66hwCLqTFEbiZEZ85IYUhaEyob48e7vmapKZfmOBiFasSWaqO9RVJTNpM4B5OsnM0vpg/bxV64Ona1O2XsnMSm3K6ufpFDXVtfT17HQkEAqC0DTkrrEgCB1Pu6bG0nRBEISOR1aEOwVNR49Hob+Xwq4k6VEfuVGFOVgi3p9lb9c0+2NTjAem6DUyBLUiJdfgZCmx9rbybWoNO612MJmgNJgiNxIkO+y127IHTAZ6F9ibmGFvZIpB3/yq2sFlNVXrFNrBU9VU7dW1gw011abuULeGmmoLEa+xIAgdj6TGgiAIbYqsCAVBaBrtuiKUQLjdKIUyfOjdKezhHnJjUTKjOrlhF2eowHDvAvuTU+wNT9V5ZD1v8bPFbo5ZSU6YCa+tfNVbbPqwTQOKmuctNuu8xVVPcXF17aCq266zodXWqbzFVZ1gq7zFFe2g29tFYSRGZtggN6wwB8tE+nOMpebYF5tid3CGHiNd0w6eLCc8P7aVYMqMMmeGvJqq6aNs+nAtrVJT1bx5spb6ildqU7bGLU5tF8pOba5qNdW6+TmtdnA9W5wCSimUPwiZLfsEN9CugVBSY0EQOh5ZEQqC0DTadUUogXA7UQotHEbr66E40kVmV4DMqEZhyCYwkGdX9yxnVSx1A8YCMb2A7WpMleM8W+zyLHVm1VIXIG8GKJoGjqWjLB1l1cll1rrbWk0yU9eFuiKbwXHRys5yycySNK6plrqlXaiVQvn9qGSc0mDSa7c1opEfdHEHCgz1LrA3Mc14eKZBMpNxghwvpjrHhkhlJ794FLs7CNNb+kle/Cdi8/KX7TAESiAUBKFptOuKUGqEgiB0PLIiFAShabTrilAC4TahDAMtFoPBXnK7EmRGDXIjYA0V6epPszc1w1nRScYCM3TrWYJaEdPxn1r+UW0dZWm1HepWbR3VDvKP+voggK6jRcI4fSnywzGyIwa5IU8yE+3PsqsimRkLzNYkMyVX53gpxbPFFJNWjIlCjAUr6M1ZwV/XbkuhWdpibXCpDbG8wTkrO2vfwqCJkhlCQehJURiOMd9rw2Nb8jFeRrsGQkmNBUHoeGRFKAhC02jXFaEEQkEQmobrKtxNBrLNvn4jSCBsNRUdnN6Vwh7pJbsrXLHUOTBosat3jv2JKfaFJxn0zZHU8+jKIWOHOFrq4piVqrWOWqjq4KyKpW5J66iqdlAvnkY7uNFtJ6s1rSZuO7lafVD5/ahUguJQstJuy9vCgAGT4W7Phrg7NMOgf56knkPHJeMEV9dbWgaO6ektF1uUbZ3esqXawWSc0kCy1pYsmypszWf5DEICoSAITUP6EQqC0PFIjVA4PZon/1D9PVijKdJjAbIjisJwmXB/jvHuWc6OTbA7OEO/b56IZmG7GidLyUqK50lmZmspnp+SZeCYxmIX6pKqpHpeOlefFqsy6HX2sKrko9pxppbiVbvN2Mu7TzdY6qoymiakeKumxLqOqkpmRqK1LjPWQIlYf5bdqTn2RycZDc7SpWdrkpmJcpxnrC4mrBjTFZlR1gxgNexQt4Jkpi4t1ktV61ydZKaue3eDZKbamWe9kplTdeaB9XWZiYRxepOetGi4Ii0aKhOMtqj1TBsjgVAQhKYhN0sEQeh4JDUWBKHjkRWhcEqUYaAlE7hDveTG4p5kZgSKQxa9fWn2JafZH51kzD9Dt5HFp8qYjr9mDzthxtdX66qzh7Wk1rUFljoVDEAqQXEwSXY04ElmBh3UQIHR7gX2J6YYD0/T50vXJDPzdpgTpSTHrSQnCzHmrDALhSAFy++1KDMNVNGrpy5tUVatp6pldrr6edtgi7Im1lMb5gxQfp83TwOJRWnRoIM7YDLUs8CYcZyntuZjfcYggVAQhKbhNiE1lhWhIAhtjcviQnYz12g10nRBEISOR1aEW41SaIEAWk83pdEesmMhMqMa+WEHfaDAnt45zk5MMh6aYsg3T0wroCuHeTvCs8UujlnJmqWu2lq+aPmwTX2ZpU5fzR5W31reXt1S19BafiVL3VJ72FotdWupddXVB5Wuo2JR7L4UheEImRGd3DAUB0ok+zKMp2bYG51mODBHr5EhqEpeu61qizIzwXT9rn4FP2VLB8trUabXWxBX0Q7qpXqtZbUN/w7a1Q9QuoaKRDzt4EiU7LBObkhhDZaI9Xkayz3RafpKk/zz1n3CG3BQKHGWCILQybTrXWNJjQVB6HhkRbiVaDp6NAIDvZijSdK7/GRHPNtTbCDDntQs+2OT7A5O02ukiWgWJdfgZCnJsWKK42ayYqkLk6mz1HmSGa3WNaWa4ulr6JpSk8pUH6/FUldNf50V0rxqWmzbtRR4w5Y6TUMF/NCVxBpKkB3xkx3WKAw66P0FdvfMsy8+za6Qt0NdTPe6qszb4c6xIFLpQO33Lc7TsJ9cRTKj9RfY1TPP3vh0TVpk5LNb9hFfiuMqlAiqBUHoZKoxfrPXaDWSGguC0PHIilAQhKbRrjdLJBBuEcrnR+tK4g72kNsdJTNikBtxKQ+aDPR5HZXPikwy4p8lqefxqTI5J8DxUmp9EhBrqWTmNJa6nSQBqZfM+AxUNIrdnyQ/EiEzbJAfdikOFEn1ZtiTmmE8MsNYYJYuI0tQlTBdHxOlBMesJCfMhNe1uzpfZqVrd1FbnK8VO1BXpDKlpbXUxR3qdsx8UZHMxGKVeQqTGarMU3+RVF+G3cnZBmmRT5VJt1COIoFQEISOp11vlkiNUBCEjkdWhIIgNI12vWssgbDZKIUWCqH1dlMa6SazK1ix1Nn4B/Ls755jf3yKPaEphnxzhDULgJlylGeL3V77KDO24o5rFDU0U2vcca2u3ZZWsYUtrXedbsc1Tw/X4h3XqvVBTUMFg9CVwByuagcVhUEbX3+BPd3z7ItPsSs4S79voaYdnLGjW6+1XI92cC1ay81qBwMB6EliDsZXnKc9sRnGQ9MN8zRvR/i5Gd2iD/tyvF9zszXCJg1mHUhqLAhCW3Ps2DH+4A/+gO7ubkKhEM973vP47ne/u65ryIpQEISm0eq7xnNzcxw4cIBf//Vf5ytf+Qq9vb385Cc/IZVKres9JRAKgtA0XDbfT3A9r7/55psZHR3lzjvvrJ0bHx9f93tKIGwmmo6eiMNAD/mxBJkxH9kRry1Ssj/D3q5pzopOMhaY8dpHacWKtzjBs8UuTpgJpqwoswWv3lVoaMe/WO/SV2jHX6sPrlTvKtXXvapHRQ+3kXqX4zTWuqCx3rUe7WAshj2QIjcaJjOskx9yKQ1Y9PRmGE9WtYMzJPV8TTt4vJTkmJXipBlnzgwzZ4bImX4sq6IdXNKebGkr/tVqg6u1J2PJnK2qHVxvezJvYk7fit8wUPEo5f4kudEw2SGd/LBLqb9xnkb8c3TVtnnwcaKU4piV5FgusHWf+S0knU43fB8IBAgEGn+XL37xi1x66aX87u/+Lg888ADDw8O8/e1v5y1vecu63ktqhIIgNI1qarzZA2B0dJREIlE7Dh8+vOz9fv7zn3P77bezf/9+vvrVr/K2t72Nd7zjHdx9993rGresCAVBaB5NzI2PHj1KPB6vnV66GgRwHIcLLriAD3/4wwCcf/75PPHEE3zyk5/kiiuuWPNbtlUg/MhHPsL111/PO9/5To4cOQKAaZr86Z/+Kffeey+WZXHppZfyN3/zN/T397duYEqhDB96dwp7pJfcWIT0qJe+OEMFhnsWODs5yd7wFEP+Obr0LLpyyNihiqXu1BYxzVye5i2meitLZlayiK0rzatPideb5sGp02JNQ4VD0J2kMBSrSUHMQZtAf57xrjn2xqfZFZyh37dQkxhNlWMrlhBMy0fJNHAr9sNquy3dWqGEsMoOdctKCE61VVmLSgiVny2VzKhwCLc7SWE4RmbYR646T30rz5PjakyV4xyzUkxYcSbMGNOZFiZ+TbhZQuX18Xi8IRCuxODgIM95znMazp177rl87nOfW9dbtk1q/Oijj/K3f/u3PP/5z284/653vYsvfelLfPazn+WBBx7g+PHjvP71r9+mUQqC0EoOHDjAk08+2XDuqaeeYteuXeu6TlsEwmw2y+WXX84dd9zRcFt8YWGBv/u7v+NjH/sYr3jFK3jxi1/MnXfeyXe+8x0eeuihbRyxIHQmrtucY628613v4qGHHuLDH/4wP/3pT7nnnnv41Kc+xaFDh9Y17rYIhIcOHeLVr341Bw8ebDj/2GOPUSqVGs6fc845jI2N8eCDD656PcuySKfTDYcgCJunmTdL1sJLXvISPv/5z/NP//RPnHfeeXzwgx/kyJEjXH755esa946vEd577708/vjjPProo8t+dvLkSfx+P8lksuF8f38/J0+eXPWahw8f5sYbb2zOAJVCC4fR+noojnSR3h0kO6IoDNsEB3Ls7p7lrLjXjn/A8KxPdqWO80urh4linJOFGHNWmIVCsLG9fFFDWXVymdV2qVsqAzmNpU4rOy1pLw9L6oOatigFGUyRGwl5kplBF3vQpK83zZ7EDHsi0wz65+nWs/iUTc7x16Qgk5ZnP5wvhFbe0a9qpVvrXG1UMrPRHf28SVmbZCYZpzSYJDcSIjusefM0YNLXs3yeNOWQdwIcL6Y4biWZtKK1nQ/TubrtEM5AXvOa1/Ca17xmU9fY0SvCo0eP8s53vpNPf/rTBIPBpl33+uuvZ2FhoXYcPXq0adcWhI7GVc05WsyOXhE+9thjTE5O8qIXvah2zrZtvvWtb3Hrrbfy1a9+lWKxyPz8fMOqcGJigoGBgVWvu5IwUxCEzSPdZ7aAV77ylfzgBz9oOHfllVdyzjnn8Bd/8ReMjo7i8/m4//77ecMb3gDAk08+yTPPPMNFF120HUMWBKEN2dGBMBaLcd555zWci0QidHd3187/0R/9Eddddx1dXV3E43GuueYaLrroIn7lV35ly8enDAMtEYeBXnJjcTJjBtlRKA4W6elPszc1zb7IFGOBGbr1LEGtiOn4OVlO8IzVzaQVY8qMMmeGPD1cwe9Z6iytoofTPL3gCu34l9e8XLQ67eApLXVlp6YZrOnhVmovX99ifmm9C9ZlqVO6DqEg9KQoDMdIj/rIDykKQ2XCfTlGU/Psi02zKzRNr5GpaeJOlhPe1gWVrTprc3U6++HSuarXWtpVPeEa5mo92sGV5grWbj9kUTvo9KUoDMXIjBrkhhTmYJlIX46x1Bz7YlOMBmcb5mmilOTZYooTZpxpM8qCFax9poq58tb+IdTTREF1K9nRgXAtfPzjH0fTNN7whjc0CKoFQWg9smdJi/jmN7/Z8H0wGOS2227jtttu254BCYLQ9rRdINwRKIXy+9F7urGHusnuCpMZ1ckNOzBksrt3jrMSk+wJTTPomyOp52uWuqOlrlrnlBkzwoIZ9DqnmH5PBlKs65xSkYGo8hZ2TtlKGUi9ZGaJFCQ35OIOFBjoSbM3McPeyBSDvvnKjn42GSe4TAoyXwhSsPxYpm+xY3ddl5maZOZUHbvtxbnyvjooBzTb2d6O3YDy+1CJOMWhJLmRoCeZGXJgwGS4e4E9iWnGwzO1eapKZp4tdnnzZEZX7mye2/K/iEa2IbXdLBIIBUFoGpIaC4IgtOnNkh0tqBYEQWgFsiJcL5qOFgmj+nuwRlOkxwIVS12ZyECOPd0z7I9Osjs4Q79vnohmYbsaJ0vJSi0nwUQhzrwVqtRy6ix1dS2ktLoWUvX1QVWu7MRWqXlVpR/r3nVtqaWuybuugSeZUZEwTm+S/HCM7LBRaSNVJtqfZVdFCjIWmKXHSBPRipRcneOl1KpSkEXJjEKztMXd/KrttjYrLyo7jXXU6hxslbwIULqGikQq8xQhM+JJZqyBEvH+LOOpWfZGPclMl54lqJVwXK3Swi3FCTPBdH0Lt4KfsqVDtS1Z3t7SP4lGVOXY7DVaiwRCQRCah6TGgiAI7cmGVoS5XI6PfOQj3H///UxOTuLUbcgD3j4CgiB0IG26ItxQIPzjP/5jHnjgAf7wD/+QwcFBlGp9Tr8dKMNASyZwhzxLXXrMIDfiUhq06O1Lsy85zdnRCUb8s3TXdhPzb67mtWqL+Z1b8wJPE0cqQXEgsaImbn9yit2hGQb98yT1HDouGSd4Wk2csvQ6O93GW5PtBO2gUgr8PuhKUhyMkxkJkBvWyA86aP0FdvXMszc+zXh4mj5fmqTuCQIzdqhWb54yo8zW5mnlFm6uuYV/FEtpRveYdpHPfOUrX+Hf/u3fOHDgQLPHIwiC0HI2FAhTqRRdXV3NHosgCG1OR7Xh+uAHP8j73/9+7r77bsLhcLPHtPNQCi0QQOvppjTaQ3YsRGZUIz/soA8U2NM7x9mJScZDUwz55htSmFUtdatsRq5ba7PUbXiXuq3urFwvBRmJkh3WyQ2DNVAi1pdld2qO/dHJWveUoCpho2rlg0krxkQhdsryQa0jT7MkM/UdeVopmYnFsPuT5IfDZIYN8kMuxYEiyd4s46kZ9kanGQvM0mVkCaoSJVdf7MZzGsmMXrfroWu2MNXspBrh//7f/5uf/exn9Pf3s3v3bnw+X8PPH3/88aYMThAEoRVsKBBedtllTR6GIAhnBJ10s+SGG25o9jgEQTgDUK53bPYarWZTzpLHHnuMH/3oRwA897nP5fzzz2/KoHYUmo4ejcBAL+ZokvQuP9kRhTlUJjaQYV/XNPtjU4wHpug1MsS0Ag4ak+UYzxa7OW4lOWnGVm6PVNTQTK1RClIvmVmpjVS5WhPcWW2klFIQCEAqjjWUIDvsr5OCmA1SkGobKV055CptpJbuUJe3fI2tycwzY4c6pZTXqbsrgTkUJzvsJzusKAzaGH0F9vTMsy8+xa7gLEP+OWKaia4c0naw9nmq36HOa+F26rZkrtWaPxXvF6VzaoSTk5O88Y1v5Jvf/GZt06T5+Xl+/dd/nXvvvZfe3t5mjlEQBGFL2ZDF7pprriGTyfD//t//Y3Z2ltnZWZ544gnS6TTveMc7mj1GQRDahU7azvO+++7jP/7jPzj33HNr557znOdw2223cckllzRtcIIgtBmdlBo7jrNMMgPg8/mW+Y7bFqVQhg+tK4k72ENuV9TbUWzEpTxoMtQ/z1nJKfaHJ9kVmPZ2qVMliq7OjB1tqOfMWWGyVoCC5W+sD1panUVsDfVBe7E+qOxKm62V6oOOCw6VetUq9cG1ttyq4qz86VRKgc9ARSPYfUnywyEyQ4uauFRfhj01TdwMvUaaiOYVrWbK0YY66pwVZr4QXF5HrdS9tqyO6rqr1wfXWkeFxeetMGe1eYrHKPcnyA2HyA7p5IdcSgMWPb0ZxpMztV0Pe40MQa2I7WrM2tGa9fBkwZunhcq2BUVzuaVOW2I9dEvN/uM489hQavyKV7yCd77znRw/frx27tixY7zrXe/ila98ZdMGJwhCm+E26WgxGwqEt956K+l0mt27d7N371727t3L+Pg46XSaT3ziE80eoyAI7UKbBsINpcajo6M8/vjj/Md//Ac//vGPATj33HM5ePBgUwcnCILQCjasI1RK8Ru/8Rv8xm/8RjPHszPQdLRQEK23m9JQatFbPOSgD+bZ3zvLecnj7A9NsNs/7bVPVzYlV2PGjpBzAuQdP5ZjYNo+SrZO2dawbQ3XUeAolK28llAuKEdVanmLglTv/PLvcUG5ric6daq+WG/YtfNu5Xz163pQCgW4GuBUEwYHdFXRE+reKUOB0lDhIE4ySmEgQm7QR26wURN3dmKSfaFJRvwztTqqjca8HWbBjrBgh0iXA+TLfsyyQcnWscsablmDsoYqK1RZodkszln94dR9dfDqpm71XHUOaKib1uapSt18Vb9f21SpU8+T7s0nuo4KBbFTMcyBMNkhg/ygwhywCfTnGe+a4+zEBHtCU4z6ZunWvRZuNl59cN4Os1AOkS4FG+apXNJxyhrYCsoKrVydn8Z5opVl+zPdWXLLLbdw1VVXEQwGueWWW075XJHQCEJncsY7Sz7+8Y9z+eWXEwwG+fjHP77q85RSEggFQWgr1hwIn3766RUfn1EohQoE0KIR6E7VrGLZUY38sE1oKMvZvZNckHyGF4R/yR5jliFDEVUBdOXDcktM2HP4lY1flfEpG59mYygbXXPQNIe8FqCouTiajq1ruJbmZRMaoHmplncoHB30koujKzQd77AVrl5pKaVrOI4nq/HkNFpt9zpcd7HFVDVNrjyutYNyWEydT5UWVjuQK4WrFBgark+nHPFRTBgUunTMXkWhz8HtNxnqXeC5XSd4fvRZzgkcZ7exQI+uE1UBACy3yISdJaJZhDWLsFYkpJcI6lF8mkNac8jrfkq6gaMbODq4lTnQdXBLClcH16g+dnHL4FTSQ63sohmq1nLLLXtpsTdP2mIbrrXMU13ZYSPz5Ph17IgPK+mj0KXV5kkNFBjrmee81AnOizzLc4LHGNWz9OoGIeUHoOCaHLdzhJU3RwGt3DBPmdo8ubh69fNUmZvKnLk6uOVN/E2slzbVEW7orvFNN91EPp9fdr5QKHDTTTdtelCCIAitZEOB8MYbbySbzS47n8/nufHGGzc9KEEQ2hPFkht8Gzm2YdwbCoSu6664YdP3v/99aeEvCELbsS75TCqVQimFUoqzzjqrIRjatk02m+Wtb31r0wfZMlwX17KwS2W0YolAsYRW6kI5QZSjk3Ni/L+yTr7sZyYR4Vj4JGO+mYp8pkjJ1Ziyu/lFsYdni12cMBNMWVFmC2Gv7bzlo2QauJV26l67JNXYdn7VXeuo1bw821hjm35Vs41V7GH1X5e0lz9la6nqPCxtN185X/0XV0rh13UCwQCRVAxzIEJ20CA/G+L4vJ+pvig/7+nh7PgQe0JTS+QzfmbtLo4Wu3nG6mLCijFtRpkzQ2TNAJbp29zufnVzdcr2/EvnqTofFQvd6Sx23pScep50pdB1nUAoSKQ7htlfkc/MhXlmIMBkf5SfdvXww/gQ+0OTDPnm6uQzQabKcX5Z7OGYlWTCijNtRphvmKcl7flX2NmvtW24znD5DMCRI0dwXZc3v/nN3HjjjSQSidrP/H4/u3fv5qKLLmr6IAVBaBPa9GbJugLhFVdcAcD4+DgXX3zxio0XBEEQ2o01B8J0Ok08Hgfg/PPPp1AoUCgUVnxu9XmCIHQYZ/qKMJVKceLECfr6+kgmkyveLKneRLFtu6mDbDmOjZPN4haL+EyLhNmDYUbRTZ2cGeEnRZ1sMcB8MsxcJMKuis3Or2w0HMKaRVArEdDL+LUyPt3G0G00zUDpLq7mQlUvqMDVXFxN1fSEtfMrHGqxSEftE6MBDp52TbkNd91cpVa22a1mIztde6m67SmhhFswUdkc4fko/tkEwYUw2bROLhPhpzk/s70hJlNR9kWijAVm6DI8/aCGQ0wvkDAK5OwAeaNIwfBRNHTKho5ddnANBbbu6QR1tfjVqejknMqh0aC/dDUXV69YEpXCVa43XRrgePOslvy1uZXat+dhVKvPT8NUrW+etGyOyFwM/3yCwEKQbFonn4nx1KCf2Z4w08koeyJeu7JuPUtQK+JTZRJ6nowRJGsHyBs+zNo8aTi28mx2ulvRVipcG+/QAZsN3hLdGGe8s+TrX/967Y7wN77xjS0bkCAIQqtZcyD8tV/7tRUfC4Ig1DjTU+N67rvvPqLRKC972csAuO2227jjjjtq7fpTqVRTB7ktVKU00zNoVpFYvgcjn8Ao+MiaIY6ZBvmij0xXgEw0yFhghj4jTVhZdBtZHDRsV8N2FY6rcGsHlFxwXC9Xc4BKvxfqpaTKVTi1T8SSn7ug4VZe73iZjwGqDC4aygUXB7WYM4OueRKaSscUhebJQCpPAUDzzjV0n6mmfZpaTPXqHruuC8USzM6jF0ximQT+hTj+tJ9sxs9cNsX380FmuiPMxCMNu7Ml9Tya30HD9Q7lHQBKuViq0jBGaaA0XLw0F1TdVFXT2SXztOh2Q6vMk1IKDQdVBsfQvM4syutkjXK9uQNvnjRtcQ6bNU9WEaZn8eULJBaSBBZi+NM+cpkg01kfmb4gM10RZuMRdgVn6PctENNMuo0sunIqhzdHqnLkFRQBR+k4Sl9UnihqdRS3eLoPexNp00C4oerBu9/9btLpNAA/+MEPuO6663jVq17F008/zXXXXdfUAQqCIGw1G1oRPv300zznOc8B4HOf+xyvfe1r+fCHP8zjjz/Oq171qqYOUBCE9qFdb5ZsaEXo9/trTRf+4z/+o7ZzXVdXV22lKAhCB9JJ23m+7GUv47rrruPAgQM88sgjfOYznwHgqaeeYmRkpKkD3G7cchl7fh6tWCRcsDDy3Rj5ABnTx7yZ5AeWn3R3kHQ8yJ6QnyHfHBHNotdIo+F4rbiUg6E56JqDrrnktEr9S+nL619qsb6l1eQcsKz+pdZe/3JxvK7YtV/KXV/9i4qlTFOL9bBK/UsphVuV3OQLKMsimC3gS6fwZypSmmyEn2f9zPcFmUlGWIiGGPHPVXa0K7InMElYtwjpRfya17psXnfIai6W5qOs697bapr3vprrtZuqyY0qnaKVW3vslr2/J606V5qGVnZxlObV18oOKA3KDhiat4sdzmKdsGqhq8qR3MpSx3ZOPU9Qmc/KGkNzls2Tm8ujrCKhXAFfOkkgHSKb0clnovwk52euN8RsMsyecIjhgGe5i2km44EpAloJv1bGr9n4tQgLukNWD2DqPuzaHGm4Si3Kslq5i10n1QhvvfVWDMPg//yf/8Ptt9/O8PAwAF/5ylf4zd/8zaYOUBAEYavZ0IpwbGyML3/5y8vOn6pztSAIZz7tWiPc8OZNtm3zhS98gR/96EcAPPe5z+W3f/u30XW9aYMTBKHNaNPUeEOB8Kc//SmvetWrOHbsGGeffTYAhw8fZnR0lH/7t39j7969TR3ktuO6OPk87okyftMime/BV4iSzuvkC1F+YhlkiwFyyQD5cIAR/wxJPU+3kSWolTy7nVbCUA6G8mqFWc3F1HzYmtFY/1KuV99Zof5FxUJX1YepWq1QQ7NdHFarfymU2vr6F4BrO7CQRjMtYukE/oVEnaawi0y/pymcT4TYFQx6Wjm9wJBvjrBWJKwXiRgWQT3OvFFmQQ961jvd19iOXsOz22mVMVW2OtA0tzJfi7VUV/PadHm1Qm8uNKVBuTKXtoNSLq7D4le78tfouN7fpeu18ELXFq10FSuecrztBl23qht0FotOjubNE57lrzZP5XJNexlPJ/EvxPCn/eTSAaazXaT7KvMUD9c0hWHNYpd/mphmEtUtQnqJgBHFr9ssaEFM3aFkGDi6jlvXtp9W6gjblA3VCN/xjnewd+9ejh49yuOPP87jjz/OM888w/j4uGzcJAidzGY6U1fT6nZZET7wwAM89NBDDd2ou7u7+chHPsKBAweaNjhBENqMTkqNA4EAmUxm2flsNovf79/0oHYkSy13hV6MfAJf3iBbCHGssGi5y0f9tQ4iEc1iT2CSgFYioJUJ6JVuNJVdyAqaQ1k3cJWOrWmLqbC2mNpV7VKuqutSo7wMFlVrquLJacoAnhwEqKSu60z7KhazhvQYvTHtc1wv/V4xPbbBNFETRYI5T0oTyITJLBjksxF+kfOT7gswk4ywtyKl6TKy3nwFLRJ6gahuccJIENDLLPiCZA0H0/Bh6wa2oXmpX91ubdVd/vRKekxdJxqtrCopskIruSilgVaxqTleZxqlLVrtsJWXx1ZKCmhe+uwq5c2TW5mH00mOKv8U9aUEVWtUUyc5KpbqpDRhMhUpzU9zfuZ7Q8ymFqU0ST1Pr5EhrFkkjDwxX4qTRpyA4c1TzvBj6X5sXffmydKguA2Rpc3YUGr8mte8hquuuoqHH3649g/60EMP8da3vpXf/u3fbvYYBUFoF9wmHS1mQ4HwlltuYd++fVx88cUEg0GCwSAHDhxg3759/PVf/3WzxygIQpuw6R3smiC/2QjrSo0dx+F//a//xRe/+EWKxSKXXXYZV1xxBUopzj33XPbt27dV4xQEQdgy1hUIP/ShD/GBD3yAgwcPEgqF+Pd//3cSiQR///d/v1Xj23Ess9zlujAKi5a7J4o+st0BsvEAu4PTDBiePGTUN0NEswjrRUJ6qda5ekELUtBdipqLoxk4mlaThNTLZxq6L1fqXW55UWqj2d6OZY6qykPwaldlF1d5O7tVbXU1e131v17HwdW1xV3bdG25PMR1G+tfq8hDoCIFAtxSGXchjWaaRNNJ/AsJAlUpTS5Fpi/IXE+Y2ViE8dB0TUoz4p8hqptEDYuIkWTKiDKr22T0AHnd70lENB3X0CqdqxulNK7u2cpczcXRlVf6rEhoalKackVKU96glAa8+WuWlGZuAT1vEsska23McukAMzmDTD7gtTGLeW3MqvNUL6WJ6EWmjCjzRogF3aunlsyKlKZYLWYKq7Gu1Pgf/uEf+Ju/+Ru++tWv8oUvfIEvfelLfPrTn8ZxZKIFQaAzaoTPPPNMQ5utgwcPopTi+PHjTR+YIAjtR7vWCNcVCMvlMsFgsOGcz+ejVGplewtBEITmsq4aoeu6vOlNbyIQCNTOmabJW9/6ViKRSO3cv/zLvzRvhDuRquXu+En8BZNkoRdfPkI6r5MzYzxl+cgWA2QSwQbLXa+RJqhKtVZKQd1rpzSnO+Q0P5buYuvGyq2UKprBms1Oq3tcsZO5aolOrux6mwHYePa7SmuuRZ1cpdZV1clVa4cb0cmt0Jqr+jq3YKImpwnkTYx0En8mgi9jkM9E+Hml7dRcKsR4JMxYYHaZpjBuJAkbCWaMiKcp9AU8TaG5RZpC5a4wVy3QFJom6uR0rY1ZIB0mk13UFM72hJhLhRmPLGovu4zsqTWFPh+lUot3lWxD2eKGNniv5w/+4A+aNhhBENqcTnCW3HnnnVs1DkEQhG1jQ4JqQRCEldjumyUf+chHUEpx7bXXrut1G+5H2PEs2+7T8x4bpkHWDHK00E2u6COTavQex/QCe9SU13JKK9a8x3N6iIzuYGoOZc2Ho7m4mlY53Jr3uEFbuNR7XKkXenrCapsurzWX9wlj0XtcxmtNr2mLOjm8ay7Tya3kPXZq4rhTt+aq1gnrNYWZJP50An86QC7jZz6b4r8qWrn5eHjNmsKC4aNo+pqvKVRurRa4LZrC+jZm6QSBBT/ZdIC5jFHbGnU+HlqmKax6tGOG5c2TESZjBMhaLZS3bWNq/Oijj/K3f/u3PP/5z1/3a2VFKAhC25PNZrn88su54447NrSvugRCQRCaRjNT43Q63XBYlrXq+x46dIhXv/rVHDx4cEPj3vGp8bFjx/iLv/gLvvKVr5DP59m3bx933nknF1xwAeClFDfccAN33HEH8/PzHDhwgNtvv539+/e3ZHzLLHeFLox8kGzez5yZ4gnLf0rLXUArL1ruNIe05pDXXEq64VnudHAttaLlDrXEcle/e9tKlruKrc5VLhrUUr1Fy52XDq7Jcqcqko91WO5WktIE0hF86cX2XAu9wfVJaartuSpSGsfQ0FeQ0mi6VzKgJkdaIqUpN6k91wY7fa8opZkoEsybp9wRsH6eknqOSMiT0kSMFGGjyIwRQSu2cL3TxNR4dHS04fQNN9zABz7wgWVPv/fee3n88cd59NFHN/yWOzoQzs3NceDAAX7913+dr3zlK/T29vKTn/ykYen70Y9+lFtuuYW7776b8fFx3ve+93HppZfywx/+cJn4WxCE9uHo0aPE4/Ha9/X65frnvPOd7+RrX/vapv7ed3QgvPnmmxkdHW2Q7YyPj9ceu67LkSNHeO9738vrXvc6wPND9/f384UvfIE3vvGNLR+zIHQ0TVwRxuPxhkC4Eo899hiTk5O86EUvqp2zbZtvfetb3HrrrViWtaYN5XZ0jfCLX/wiF1xwAb/7u79LX18f559/PnfccUft508//TQnT55sqAskEgkuvPBCHnzwwVWva1nWsvqDIAibp9XymVe+8pX84Ac/4Hvf+17tuOCCC7j88sv53ve+t+ZdNXf0ivDnP/85t99+O9dddx3/83/+Tx599FHe8Y534Pf7ueKKKzh58iQA/f39Da/r7++v/WwlDh8+zI033ti8gdZb7kyLVL4XXyGMYerkCouWu1wyQD603HIX1q2a5S6gRzzLnb42y52mrcNyp1w0ewdY7laQ0lTbc61VShPzxZkwlux0Z9btdLeClKa6092qUhqNzbfnauZOd2vdETAeZldohkHffIOUJmGkiBhJDKuF22e0WD4Ti8U477zzGs5FIhG6u7uXnT8VO3pF6DgOL3rRi/jwhz/M+eefz1VXXcVb3vIWPvnJT27qutdffz0LCwu14+jRo00asSAI7ciOXhEODg7ynOc8p+Hcueeey+c+9zkABgYGAJiYmGBwcLD2nImJCV74wheuet1AILBi4VUQhE2yA7zG3/zmN9f9mh29Ijxw4ABPPvlkw7mnnnqKXbt2Ad6Nk4GBAe6///7az9PpNA8//DAXXXRRS8cqCML2W+w2yo5eEb7rXe/i4osv5sMf/jC/93u/xyOPPMKnPvUpPvWpTwHUPIV/+Zd/yf79+2vymaGhIS677LLWDrZquZuaRisWied7MXJxjPwSy11XgHykvSx3qmqxa6LlbtX2XGmvPdcvsn4W+oLMJMPsjU6vqClcactPy/BR1jepKWzWlp9boilMNm6N2jBP3pafvUbG0xRW5ilgLrbIE1ZmRwfCl7zkJXz+85/n+uuv56abbmJ8fJwjR45w+eWX157z53/+5+RyOa666irm5+d52ctexn333ScaQkHYDnZAarwRdnQgBG8P5de85jWr/lwpxU033cRNN93UwlEJgrASzUhtJTU+A3DLZey5BTTTIlLoxZdPLVrurBT/ZflJdwdPabnzV46mWu7KFTXLBix3NWlIMy13lQ2/lklpqjvdZf3M51I8UdnpbiUpTULPEzdMjvsSp5fSGN68uZqXKp9WSrORne6WztcWSWmi6fjiPKX9LGRTPNFftyNgeJo+X5qknmPEP4Mezm7hJ/7MQAKhIAjNQ1JjQRA6njYNhDtaPiMIgtAKZEW4FTi2Z7k7dgJ/wVxiuYs2WO6ssI9B39wyy11IL+HXbQJ6uGmWOyoSmh1huatrzVWT0pyc8qQ0mXopTZhf5HwrSmmSep5gsERMN4nqFhO+OEEjwrwvtCilMXRsQ8exNK9cp4Nb9Oar+VKaJfO1VVKayeJiG7OFCJm0Ti4T5hdZP/N91R0BZxgNztKlZ0lo+eZ/xleh8jHb9DVajQRCQRCaR5umxhIIBUFoGu0qn5EaoSAIHY+sCLeKlSx3+YrlrhDkqLlouctEgqta7k7oiUbLnXIpa8aGLHeu8spTO9JyB7i2vWlNYcRIENQ9TWHaCJA3/JQMA0c3zixNYTqDZppEMkl8Ga89Vy7jI51N8oP+ILM9XnuusdAsSWd6yz7my5DUWBAEgW0JZJtFUmNBEDoeWRFuMUstd0a+C6MQIFvwdrmrt9ztDU7Ra6RPb7nT/eu23NVLaHaS5W5perxWKc3Sne5WktJMGxHmjRBZI4Clu6eV0uiaQtPcyly2gZTGKnrde3IFjIUUgYUIvrROPhPil1kfc30hZlIRRlhbu/pm0K43SyQQCoLQPNq0RiipsSAIHY+sCAVBaBqSGgurU2e5CxRMjFwvvnwYo7BoucuX/OQSAfaFJ09pufPpEYyNWO7U8rrXTrDcNbTmgjVLaU61091KUppaey599fZcaJVd70q0lZTGTWfQLItINoE/kyCQDpDN+MhmEzwxEGAy2LoaoaTGgiAIbYqsCAVBaBqSGguCILRpaiyBsFWcxnL3S7OHjOU/peUuoJUJ6SUCenRjljvN9eqIbWK5O6WmMBvhFzk/C72epnBvNMSIf25FTeEpd7qr1xTW7XTXTppCANcqoqZm8OdNkgteey5/2iCXCXMs0d3kD/MpaNNAKDVCQRA6HlkRCoLQNKRGKAiC0KapsQTCFlPzHhdLRArWMu/xD4o+st0BsrHl3uOgKhHWihvyHlPzG7s1H3JbeI8rc+bOzS9qCtMJ/JkAuYyf+aynKaxtZblEUxjVTaKGRcwXP/2Wn5U64ZZqCgEcd0s1hco0CWcT+NJJ/OkAWo9vSz7LZxISCAVBaBrKdVHu5pZ0m339RpBAKAhC85DUWFgzjo2TzeIWiwSsIkauu85yF+PHlo9sT2CZ5W7AN09Es9ZtuaumwdVUWdOabLkDTx60VZa7SrrYIKVJJwmkI/gynpTm5zk/c72hZVKabj1LJGiR0AsrSmlMw4ddkdK4ulbrWr11Uhqv9LClUppiCTU9iz9XIJlOoffJn/npkBkSBKFpyF1jQRCENk2NRVAtCELHIyvC7aJquZuYRDPNmuVOLxjkzBC/LBjkin4yqQ1a7mo1wIoEprqb3VZY7pRabJ+1VZa7ap1wJSlNOlBpz5WsSWnm4yF2BWdXlNJEjCRTRpRZ3SajB2pSGkfTcQ3Na8t1CilN/XYIO1pKk82hikXC6eCWfpTrkdRYEAShTVNjCYSCIDSNdl0RSo1QEISOR1aE28xKljtfPkA272fGTJG3PMtdIeZjd3CG/oqWcE2WO92z3KFVal6awi2yacsdFb1gyy13qqJaXKopzKQIZMJk0ga5bISf5/3M9waZS4UZj3iawl4j3aApjBtJwkaCGSPSqCk0jYrlbrmmUK/TEdbmc62aQuWuoMM8haZwszrMOk0h0/kt+eyuiKTGgiAI25PabhZJjQVB6HhkRbgTWGK50/PdGIUwuqmTr7PcpRMh9oX967bc2ZqGa2k1+QzFzVnu3IqkZtstdw1SmhS+dAJ/xk8242cum+L7+WBlp7sQe0J++n0LhDXrlFKavO6nZBg4+spSGr36eBUpDQpvXpZKaZRbS4HXJKWBRflRtaxg2955Xa0gpdG971eS0rSyiYFbGfNmr9FiJBAKgtA05K6xIAhCmyIrQkEQmofcNRY2RZ3lTjdN4vk+fLkYRn5tlrugKq/dcqctsdxVbXYbtdzh4KJti+UO18XNF1DFEsFcAV86hT8TIZvRPSlNzk+6L8hcMsy+yBTDgTm69cX2XFHdJGKkCBtFZowIc0aInOHH2qCUxtVUxXa3XEpz6vrqGqQ0ul6pGa6vvtrKwKIcFt93E9doNZIaC4LQ8ciKUBCE5iGpsSAInU673jWWQLjDaLDc5bsbLHfTVnJVy90u/xRhzVrVclfUvB3bHA3PcletcVUtd5XzG7PcNWritsVyV6cpjGVS+BcS+NOepnAmlyLTH2CuO8x8PMyu4IzXnksz2eWfJqEXSBgpIkaSgFFm3giRMZZoCvVNagqVp4X0NJtbrCmsr68Cyq2pO7eeNtURSo1QEISOR1aEgiA0DUmNheaxquXO2JDlTtdc8rqDpfuxtcqG5ppWS+POFMtdg5Qmb+LLpPBnwmSyOvlslJ/m/Mz3hphNhdkTDq0ipbE4acS9rjSGvSilMSo73VmLUho0hbNUStPwuE5Ko1yU5m5MSgNeWWGtUhq9Mh9189Yy2vRmiaTGgiB0PLIiFAShaUhqLAiC0KZ3jSUQ7lRWstzlK5a7QohfmouWu3zUz7B/jj4jvaLlzqfZLBhB0pqLpfko6xXLnd4Cyx0tttzhSZDc2Tm0QoFYtgt/Oo4/7SeXCTCTM8gWAsx1rSyliWkmUd0iohcr7bnCnlVRdypdvxulNLpVqa/WbHbg6qeQ0pTq5m2tUhqlcD3t0HIpDTTOW7UdV8O8SQXsdEggFAShaUhqLAiCIHeNBUEQ2hNZEe5w6i13YbMHI5/CKATIFvxMm4uWu2wsgBX0ncJyF8enOSxoQQq6c0rLHYrmWe6Ui+usYrmzK6LALbDc1TSFxyYIZguLO91VNIU/yfmZ6w0xmwyzLxJkxD9LUs/Ta2QIaxaxiqYwbBSZ8Xk73XmaQj+2rjdoCjUd9KKq6Qvd0uqaQhSockWfaa9DU1j5vZZpCqvzWp03112uKbRbpyOU1FgQBMGpu4mzmWu0GAmEgiA0D6kRCoIgtCeyImwHKt5jVSwSKFikct348mGMwtq8x0GtREArEzRi+HSbeT1ETnMp6r6Vvcd1Grit9R6v4KGF5nmPWdQU6qZFLJPEV9EUZjMBpnNdZPqCzHeHmI+F2RWapt9YIKIVGQ9MktDzxI2uZZrCgu6jbBg4WkVTqFc0hXpdnbA2h42awup8rqoprGoxl2oKq1rMjWgKVQtrhDShRtiUkawPWREKgtA8qs6SzR5r5PDhw7zkJS8hFovR19fHZZddxpNPPrnuYUsgFAShbXnggQc4dOgQDz30EF/72tcolUpccskl5HK5dV1nR6fGtm3zgQ98gH/8x3/k5MmTDA0N8aY3vYn3vve9XqqAJ5+44YYbuOOOO5ifn+fAgQPcfvvt7N+/f5tH32TqLXeWRTzfu6LlLpfyL7Pc7QtMENaKXmsurYxPszH0EFnd2bzlTlVlNEssdzigtPVb7qqSmiZa7nBd3FwOVSoSypv40snF9lyZKE/lfcz3hphPhtgbCTLom6fLyK4qpZmv7HRXNHzYpr64052+gpRGV7ildUhpNHArj5srpWkNrZbP3HfffQ3f33XXXfT19fHYY4/xq7/6q2u+zo5eEd58883cfvvt3HrrrfzoRz/i5ptv5qMf/Sif+MQnas/56Ec/yi233MInP/lJHn74YSKRCJdeeimmaW7jyAWhQ3GbdADpdLrhsCzrtG+/sLAAQFdX17qGvaMD4Xe+8x1e97rX8epXv5rdu3fzO7/zO1xyySU88sgjgLcaPHLkCO9973t53etex/Of/3z+4R/+gePHj/OFL3xhewcvCMKmGB0dJZFI1I7Dhw+f8vmO43Dttddy4MABzjvvvHW9145OjS+++GI+9alP8dRTT3HWWWfx/e9/n29/+9t87GMfA+Dpp5/m5MmTHDx4sPaaRCLBhRdeyIMPPsgb3/jGFa9rWVbD/y7pdHprfxFB6BCU66I22Uar+vqjR48Sj8dr5wOBwClfd+jQIZ544gm+/e1vr/s9d3QgfM973kM6neacc85B13Vs2+ZDH/oQl19+OQAnT54EoL+/v+F1/f39tZ+txOHDh7nxxhu3buBbiFsuY8/Oo1nFFS13haKPdCm4quXOk9LY+DWbgG43Wu60DVrutJUsdxqUFwtGK1ruoHG3Nruyw139bm3NstwBbrGEOz2LnjeJZ1P4M3H8GR+5bJCJnEGmP8BsV5h9sTC7gzP0GOmalCaqe+25QnqJgL64011NSqPruMYKUppabbBOSlOtsa4mpVna1my9Upolbc1a2qrfYVEGtZlrAPF4vCEQnoqrr76aL3/5y3zrW99iZGRk3W+5o1Pjf/7nf+bTn/4099xzD48//jh33303f/VXf8Xdd9+9qetef/31LCws1I6jR482acSCILQS13W5+uqr+fznP8/Xv/51xsfHN3SdHb0ifPe738173vOeWor7vOc9j1/+8pccPnyYK664goGBAQAmJiYYHBysvW5iYoIXvvCFq143EAicdpktCML6aWZqvBYOHTrEPffcw7/+678Si8VqmWAikSAUCq35Ojt6RZjP59G0xiHquo5TkV6Mj48zMDDA/fffX/t5Op3m4Ycf5qKLLmrpWAVBoKl3jdfC7bffzsLCAi9/+csZHBysHZ/5zGfWNewdvSJ87Wtfy4c+9CHGxsZ47nOfy3/+53/ysY99jDe/+c2AVyu59tpr+cu//Ev279/P+Pg473vf+xgaGuKyyy7b3sFvJRXLnVYqrWi5e9L0kev1k00EyIf8jPhniWkFBnzzBLVSRVNYXNlyp+vYWqssd6fRxW2B5a5BU5gr4Eun8GdC+LIGuVycJ/N+FnqDzCdCjIdnGPbPkdRzDBgLxLRCrY3/pBFlxoiQNoIrawqr7bkqNcOa7U6v0xSWPV2mprFcU1j2NIWU3YrW0ps7Vfa0l6eeOxprrG4L1zst3rPEbdL+Jjs6EH7iE5/gfe97H29/+9uZnJxkaGiIP/mTP+H9739/7Tl//ud/Ti6X46qrrmJ+fp6Xvexl3HfffQSDwW0cuSAI7cSODoSxWIwjR45w5MiRVZ+jlOKmm27ipptuat3ABEFYEWnMKrQW18UxTdyllruCQc4M8QvTIGsFyKSC5KN+xvwzdBtZknqOYKB0esuddgZb7qiT0pgWiUwKfzqGL+sjlw1xPOcj3R9kNhVhIRZiLDBbk9LsC54kYeSJGykihteVZq4ipTF1H2Xdh6O7lZ3uNFzdbZTSlOqkNNXzZU65Q6CyF6U0KAVlx5s7ba07BLZQPtOm23nu6JslgiAIrUBWhIIgNA3lULthtZlrtBoJhIIgNI82TY0lELY5K1ru8gGy+UXLXbYcIB/zszs4Ta+RXtVy59McMprTMZY7XBc3k0VZFuFcAV/Gk9Jkswb5bJwfDQSY7wkxmwizNzzNoH9+RSlNuCKlmdeDFAwHy/Dh6Dq24TbudKdVpDQNLbnqLHjVWqHtSWlQLko5KE01SmkUYIOyXVwcULpXK7Qd3ErNUDmu9zvrVQ2UcCokEAqC0DzWKYhe9RotRgKhIAhNo9UWu2Yha2ZBEDoeWRGeCSy13BV68BXCGGZ1lzuDXK+fdDzIvvAkQ775FS13fr1MwAivbrlTa7DcKRdNayPLHRVN4dQMPtMikfU0hf60QS4b4ljOIN0fYD4VZn80yGhwll4jU9MUem38UwT1EgE9yoIvWNEU+is73VX0mFZFj1ms6C5rbfzr2nNV57TMou1O0yo6wsqclSt6TK1xzlbUY1bmEf3Mtdg1CwmEgiA0j+r2KZu9RouRQCgIQtNo1xqhBMIzhZUsdzlvl7usFebnpkGmL0AuFSAfCTDin11muQtoZYJ66dSWu2qKdyrLXbn9LHfYNk46g7KKRLImvkwSfyZIJusjn03ww8EA890hZhMR9oanGPTNE9MLDPnmiOkFYrpJ1PC60vj1CGndJm8EKBoGjqnjVFJhryMNXorc0MF6dSmNZ7urkyE5bkVKo9CU0yBDWklKI/KZ0yOBUBCE5uHShBphU0ayLiQQCoLQPNr0ZomsmQVB6HhkRXiG0WC5y9fvchdgykp5u9x1ebvcLbXcBbQSYb24zHKX11xKmntmW+4q51zLwp2awZ8vYGS68GeiZDIGuVyIo1kfC31B5lJh9lfac3UZWWKayb7ABAm9QETvImwUmTaizBk2WSOApfsoG25lpzuFq2voel0bLr2uNqgvl9K4SqGVXZSm4dqu97gipXGUhnJOI6XRWrjecfD+3Td7jRYjgVAQhKbRrneNJTUWBKHjkRWhIAjNo01vlkggPBOpt9yZRZJmD0YhhF4wyBfi/Liyy10uEWA85K9Z7oZ9s0Q0q8Fy59MjGLpDTvdvmeUOtIad2rbNcletHTpOTVMYziXxpZP4MwGyWR/ZXIL/l6+054pHGA9PL9MUJowuokaRoB5jzgizoAcpGH6KuoGjGzVNoV7d3a5YVyssraAp1E5RZz2FplCheb+b08LEr00DoaTGgiB0PLIiFAShebTpilACoSAIzUPkM8KOouo9PjmBYRVJ5Hu87T7zBrlCmKctg1zRTzoZPK33OKBHG7zHJe0U3mO1Me9xrb3UdnuP62qHrmXBxDT+XIFkpgt/Ooovo5PLhvhl3mChN8hcKsTeqKcp7DXSNU1hVDeJ6BYnjAQBvcyCL0jWCGDqLmVDx9V1T0+oKU87qC/1HlObV2+z4EXtpqcjrNRZ7SWawjJL6qytbcMl8hlBEIQ2RVaEgiA0D6kRCjsRt1zGnplFtywiNcudn1whwITZRc7yk+0KLLPcjfpnTm2501ex3GlU7HbULHeedKZNLHfVP8J6Kc1CGlUsEskl8WWT+DN+sjk/C7kkTxT8LPSESMeD7AoGGfLPEdNMRn2zxDTT2+nOSDBlRJmttefyUzI8KY1taLimWmKzW5w/rZIe65U5q82fBpRBaUvnz2vB1TB/rWzD5Sz+G27qGi1GUmNBEDoeWREKgtA8JDUWBEFoQiDchs6sEgg7AcfGzmTQikWCBQs934Mvv2i5+6HpI9Mb2DLLXYNlTLlodfXCllruXLdSmDyN5a5y/fo6oVsw4eQUgbyJkUnhz0TwZQzy2Qg/z/mZr7TnWoiEGAvM0K1nSep5IkHLa+OvW4SNEtNGhHkjVGnP5UlpbF3HsTSv3FlnudN00DWF0lmUKGkuWlnVLIxayW2YP+W43pxqbk1K4zr61n/G2hwJhIIgNA9JjQVB6Hgcl02ntnLXWBAEofXIirBTqG33OdVouSsY5MzTW+4CWun0ljtNw9XVqpY7VMVyV79d5eksd+VKra8ZljtFpTbIhix3OA7O3DyaaRLNpPCnE5X2XH7m8im+nw8y1xMmEw+yKzhDv2+BmGayyz9NVDeJGhYhvURQLzOr22T0AAXDR1H34Wo6tu7Nn1tp01XVZGolKnXXxjb+VV2hslnUFJaXawpxNmv+XQeu4x2bvUaLkUAoCELzkBqhIAgdT5vWCCUQdhhuqbjccpf3k8uf2nK32z9NWCue3nKna6e03FW7V9db7lBeFryy5W4Fy9g2W+7cXB6KpQYpTTajk8tG+Gnex3xviNlUmD1hT0qT1PN061kiQcvrSmOkCBtFZowIc0aIvOFgGj5s06ikx9pix2pdVWQ0dVIkve6xptDKK0tpNNvBVRqOLbcCTocEQkEQmoekxoIgdDwuTQiETRnJupA1syAIHY+sCDuResudWUTPd3uWO7MFljtVZ7lTnuXOVcpTqbSL5Q6gVMKZnUMzTWIVKY0v6yeXCTCTM8gWAsx1hUnHQ+wKTdNvLBDRiuzxT5HQC0T1biJGkYBRZt4IkdED5PVqey4d19BwLVVpyVWx1NW36jqdlKYqRdJc3FbWCCU1FgSh43EcNr3piNN6HaGkxoIgdDyyIhQEoXlIaiy0FbVd7iYxTKtxl7vTWO58gXLNcufpCW0M3SarO5jKT7nSwt/VKuLApZY7rc5yV79LW73ljooTbi2Wu8rvA6233LnZHKpYIljTFIYxcjr5bJSf5H0s9AaZTYTZE5lmxD9Ll5H1NIUhrz1XSC8yacSYMSIsGDY5w4+l+71aq16ptVZbchUXW/m7pSXtzeo0hV6ttdKyy3ZbqyNs00AoqbEgCB2PrAgFQWgeYrET2pEGy12hFyOfPK3lLqhKjPpmCaoSPmUT0MoYmoOuXDTlUtB8lDUXV1WlNNVUuNFWh1JoS9Li6s+8tNh7iZe2aN5m5VBJXestd6Dqfraq5Q4a02NotNw5bsXrd5r0GBqkNG6xiDs7h25axLMp/Jk4/oyPXC7IRN4g0xdgvitEOhZkd3CGHiNNRCsyHpisSJFKhCpdfeb0EBndwdT8lHUDR1U70WiNJYaa/KgipVH13X684Wklbw4dX+u6z7iug7vJ7jGbff1GkEAoCELzqP9PZzPXaDFSIxQEoeORFaEgCM3DbUKNUOQzwrawbJe7bnyF5Za7dDzIntA0/b4FwppFUs9T8hvYaJRcDcdVVP8MTFdRdhWuqyr1c4Vy3UoRq/K+lfPeFpBQ0bPUztfKeNSlLu5ivUu5Xl1MVepjSqta5yrPUXW1MY3F4lm1Tlj9g1ObrxNi27iZDMqyCOUKGJku/JkgvqxBPhPnybzfk9LEI+yNTNHnSxPTCnQZ2VX/WUygDLjoOLVZcCvzVE9l3mqPqT1Pw5PXtAzHoeZT3CjbUCOU1FgQhI5HVoSCIDQPSY0FQeh0XMfB3WRqLPIZYftYyXKXq1juCmF+nvcx2xPmRDLBSHieLl+OmG5iuxplR6Ps6NiOhu1oOI6G6yrc6u5prkJVFwp19UHlgnLdytfGx9XnVJ+3bJGxkkTjVCuJU0k6NvqHV9clpVpvdC0LNTWDkTdJpLvwZ+L4Mga5bIjjWR9zvWFOJOMMR+bpD2SI6ha6cnBche0qnMrhVuqctV+pOocN4158uOxnS34unBoJhIIgNA9JjQVB6Hgcd5Xl6ToQQbUgCELrkRWh0MDS7T59mSTBWT/5CT+F7hRPJJL8IFrGCJXx+ctomovjKIqWD9vUwdLQLA2tqPAXvdZRWom6w/W+liuPy6DZLlrZRZU9z7Bmu6hypU2/7aJszzesHNf73nHqvlZa7jvOYpt+x9MA1lpCOd5z3Kr9q1oTXPp9pY5Y0xfWt+mv/77+OVXqapCuU8ZNp1GmSTibxz+fIjQbIjepY/bE+WUyys9j/WjhMobfxjBsXFdRKhqUTQOKGsrS0C1vDrW6OdSLdfNWAq1ceVz2tjLQSm7DHCrHRTNLTf+crEp1m4RNX6O1SCAUBKFpuI6Lu8nUeNl/Mi1gW1Pjb33rW7z2ta9laGgIpRRf+MIXGn7uui7vf//7GRwcJBQKcfDgQX7yk580PGd2dpbLL7+ceDxOMpnkj/7oj8hmV1frC4KwhbhOc451ctttt7F7926CwSAXXnghjzzyyLpev60rwlwuxwte8ALe/OY38/rXv37Zzz/60Y9yyy23cPfddzM+Ps773vc+Lr30Un74wx8SDAYBuPzyyzlx4gRf+9rXKJVKXHnllVx11VXcc889rf51zhwcGzudhnQa/ScQxzu2G3fJ152Oa1k4mQz88ihhILxd43BbmBpvA5/5zGe47rrr+OQnP8mFF17IkSNHuPTSS3nyySfp6+tb0zWUux3r0BVQSvH5z3+eyy67DPBWg0NDQ/zpn/4pf/ZnfwbAwsIC/f393HXXXbzxjW/kRz/6Ec95znN49NFHueCCCwC47777eNWrXsWzzz7L0NDQmt47nU6TSCR4Oa/DUL4t+f0EYbsouyW+yb+ysLBAPL41/6XV/obUf9v031DZLfFN9/NrHu+FF17IS17yEm699VYAHMdhdHSUa665hve85z1res8de9f46aef5uTJkxw8eLB2LpFIcOGFF/Lggw8C8OCDD5JMJmtBEODgwYNomsbDDz/c8jELQsfT4tS4WCzy2GOPNcQJTdM4ePBgLU6shR17s+TkyZMA9Pf3N5zv7++v/ezkyZPLlr6GYdDV1VV7zkpYloVlWbXvFxYWAChTap+8SxDWSBkvNW5F8teMv6HqeNPpdMP5QCBAIBBoODc9PY1t2yvGiR//+Mdrfs8dGwi3ksOHD3PjjTcuO/9t/n0bRiMIrSGTyZBIJLbk2n6/n4GBAb59sjl/Q9FolNHR0YZzN9xwAx/4wAeacv2l7NhAODAwAMDExASDg4O18xMTE7zwhS+sPWdycrLhdeVymdnZ2drrV+L666/nuuuuq30/Pz/Prl27eOaZZ7bsg3ImkE6nGR0d5ejRo1tWa2p3duIcua5LJpNZc818IwSDQZ5++mmKxWJTrue6Xi/IepauBgF6enrQdZ2JiYmG8xMTE6eMAUvZsYFwfHycgYEB7r///lrgS6fTPPzww7ztbW8D4KKLLmJ+fp7HHnuMF7/4xQB8/etfx3EcLrzwwlWvvdISG7wa5E758O5k4vG4zNNp2Glz1Ir/4IPBYE3N0Sr8fj8vfvGLuf/++2s3Wh3H4f777+fqq69e83W2NRBms1l++tOf1r5/+umn+d73vkdXVxdjY2Nce+21/OVf/iX79++vyWeGhoZqv/C5557Lb/7mb/KWt7yFT37yk5RKJa6++mre+MY3bun/foIg7Byuu+46rrjiCi644AJe+tKXcuTIEXK5HFdeeeXaL+JuI9/4xjfqmzPVjiuuuMJ1Xdd1HMd93/ve5/b397uBQMB95Stf6T755JMN15iZmXF///d/341Go248HnevvPJKN5PJrGscCwsLLuAuLCw061c7I5F5Oj0yR9vDJz7xCXdsbMz1+/3uS1/6Uvehhx5a1+t3jI5wO7Esi8OHD3P99devmDILHjJPp0fmqD2RQCgIQsezYwXVgiAIrUICoSAIHY8EQkEQOh4JhIIgdDwdHwg328fsTOPw4cO85CUvIRaL0dfXx2WXXcaTTz7Z8BzTNDl06BDd3d1Eo1He8IY3LFP2dxIf+chHUEpx7bXX1s7JHLUXHR0Iq33MbrjhBh5//HFe8IIXcOmlly6z7XUSDzzwAIcOHeKhhx6q9Xi85JJLyOVytee8613v4ktf+hKf/exneeCBBzh+/PiK/SQ7gUcffZS//du/5fnPf37DeZmjNmMLtI1tw0tf+lL30KFDte9t23aHhobcw4cPb+OodhaTk5Mu4D7wwAOu67ru/Py86/P53M9+9rO15/zoRz9yAffBBx/crmFuC5lMxt2/f7/7ta99zf21X/s1953vfKfrujJH7UjHrgib1cfsTKfaoqyrqwuAxx57jFKp1DBv55xzDmNjYx03b4cOHeLVr351w1yAzFE7smObLmw1zepjdibjOA7XXnstBw4c4LzzzgO8HpB+v59kMtnw3Po+kZ3Avffey+OPP86jjz667GcyR+1HxwZC4fQcOnSIJ554gm9/+9vbPZQdxdGjR3nnO9/J1772tZZ3WxG2ho5NjZvVx+xM5eqrr+bLX/4y3/jGNxgZGamdHxgYoFgsMj8/3/D8Tpq3xx57jMnJSV70ohdhGAaGYfDAAw9wyy23YBgG/f39HT9H7UbHBsL6PmZVqn3MLrroom0c2fbiui5XX301n//85/n617/O+Ph4w89f/OIX4/P5GubtySef5JlnnumYeXvlK1/JD37wA773ve/VjgsuuIDLL7+89rjT56jt2O67NdvJvffe6wYCAfeuu+5yf/jDH7pXXXWVm0wm3ZMnT2730LaNt73tbW4ikXC/+c1vuidOnKgd+Xy+9py3vvWt7tjYmPv1r3/d/e53v+tedNFF7kUXXbSNo95+6u8au67MUbvR0YHQdTffx+xMgxX6QwLunXfeWXtOoVBw3/72t7upVMoNh8Puf/tv/809ceLE9g16B7A0EMoctRfShksQhI6nY2uEgiAIVSQQCoLQ8UggFASh45FAKAhCxyOBUBCEjkcCoSAIHY8EQkEQOh4JhMKW8KY3vYnLLrtsu4chCGtCus8I60Ypdcqf33DDDfz1X/81otUX2gUJhMK6OXHiRO3xZz7zGd7//vc37GsSjUaJRqPbMTRB2BCSGgvrZmBgoHYkEgmUUg3notHostT45S9/Oddccw3XXnstqVSK/v5+7rjjDnK5HFdeeSWxWIx9+/bxla98peG9nnjiCX7rt36LaDRKf38/f/iHf8j09HSLf2PhTEcCodAy7r77bnp6enjkkUe45ppreNvb3sbv/u7vcvHFF/P4449zySWX8Id/+Ifk83kA5ufnecUrXsH555/Pd7/7Xe677z4mJib4vd/7vW3+TYQzDQmEQst4wQtewHvf+17279/P9ddfTzAYpKenh7e85S3s37+f97///czMzPBf//VfANx6662cf/75fPjDH+acc87h/PPP5+///u/5xje+wVNPPbXNv41wJiE1QqFl1G95qes63d3dPO95z6udq+4fU91O9fvf/z7f+MY3Vqw3/uxnP+Oss87a4hELnYIEQqFl+Hy+hu+VUg3nqnejHccBIJvN8trXvpabb7552bUGBwe3cKRCpyGBUNixvOhFL+Jzn/scu3fvxjDkoypsHVIjFHYshw4dYnZ2lt///d/n0Ucf5Wc/+xlf/epXufLKK7Fte7uHJ5xBSCAUdixDQ0P83//7f7Ftm0suuYTnPe95XHvttSSTSTRNPrpC85BW/YIgdDzy36ogCB2PBEJBEDoeCYSCIHQ8EggFQeh4JBAKgtDxSCAUBKHjkUAoCELHI4FQEISORwKhIAgdjwRCQRA6HgmEgiB0PBIIBUHoeP5/5BQdJh8bVkoAAAAASUVORK5CYII=",
|
|
"text/plain": [
|
|
"<Figure size 640x480 with 2 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"#8.5.1\n",
|
|
"\n",
|
|
"dt = .01\n",
|
|
"p = 10\n",
|
|
"D = 1\n",
|
|
"tau = 5\n",
|
|
"xmax = 100\n",
|
|
"tmaxs = 50\n",
|
|
"tmax = int(tmaxs / dt)\n",
|
|
"\n",
|
|
"deltaT = np.zeros(xmax)\n",
|
|
"deltaR = np.zeros(xmax)\n",
|
|
"T = np.zeros((xmax, tmax))\n",
|
|
"R = np.zeros((xmax, tmax))\n",
|
|
"\n",
|
|
"\n",
|
|
"S = np.zeros(xmax)\n",
|
|
"S[int(xmax / 2) - 1] = 1\n",
|
|
"\n",
|
|
"for t in range(0, tmax):\n",
|
|
"\n",
|
|
" if t > 0:\n",
|
|
" for x in range(0,xmax):\n",
|
|
" T[x,t] = T[x, t-1]\n",
|
|
" R[x,t] = R[x, t-1]\n",
|
|
"\n",
|
|
" for x in range(0,xmax):\n",
|
|
" deltaT = S[x] + p*T[x,t]**2 / (1 + T[x,t]**2) - R[x,t] * T[x,t] / (0.1 + T[x,t]) - T[x,t] \n",
|
|
" if x > 0 and x < xmax - 1:\n",
|
|
" deltaT += D*(T[x-1,t] + T[x+1,t] - 2*T[x,t])\n",
|
|
"\n",
|
|
" deltaR= T[x,t] - R[x,t] / tau\n",
|
|
"\n",
|
|
" T[x,t] += deltaT* dt\n",
|
|
" R[x,t] += deltaR* dt\n",
|
|
"\n",
|
|
"plt.title(\"T\")\n",
|
|
"plt.imshow(T, cmap='viridis', interpolation='nearest', extent=[0, tmaxs, xmax, 0])\n",
|
|
"plt.xlabel(\"Time\")\n",
|
|
"plt.ylabel(\"Position\")\n",
|
|
"plt.colorbar()\n",
|
|
"plt.show()\n",
|
|
"\n",
|
|
"plt.title(\"R\")\n",
|
|
"plt.imshow(R, cmap='viridis', interpolation='nearest', extent=[0, tmaxs, xmax, 0])\n",
|
|
"plt.xlabel(\"Time\")\n",
|
|
"plt.ylabel(\"Position\")\n",
|
|
"plt.colorbar()\n",
|
|
"plt.show()\n"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 49,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhoAAAGTCAYAAABu7rurAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA+rUlEQVR4nO3de3RU1d3/8c+ZSTITIBluJeESID+h3C8REAOuihqhSKnpY636s4IUtbYBxfRpl/hT8B6ti4IIgtQibS0LqxZoqWIxCj6WqFzEB6xgUWpQSVAwCQkkmZxzfn8kGRlIMJPJ5GQy79daZ2Vmz9ln7zmCfLP39+xt2LZtCwAAIAJcTncAAAC0XwQaAAAgYgg0AABAxBBoAACAiCHQAAAAEUOgAQAAIoZAAwAARAyBBgAAiBgCDQAAEDEEGgAAIGIINAAAiDGPPPKIDMPQvHnzznne888/r8GDB8vr9WrEiBF66aWXQm6LQAMAgBiyY8cOPfXUUxo5cuQ5z9u+fbuuu+46zZ49W++++66ys7OVnZ2tffv2hdSewaZqAADEhvLycp1//vl68skn9eCDD2r06NFasmRJg+dec801qqio0KZNmwJlF154oUaPHq2VK1c2uc24cDsNAABCU1lZqerq6rCvY9u2DMMIKvN4PPJ4PA2en5OTo2nTpikrK0sPPvjgOa9dUFCg3NzcoLIpU6Zow4YNIfWRQAMAgFZUWVmp9H6dVHTUDPtanTp1Unl5eVDZwoULde+995517rp167R7927t2LGjSdcuKipSSkpKUFlKSoqKiopC6iOBBgAArai6ulpFR00d2tVPyUnNT5UsO2EpfcwnOnz4sJKTkwPlDY1mHD58WLfffru2bNkir9fb7Dabg0ADAAAHJCe5wgo0AtdJTg4KNBqya9cuHT16VOeff36gzDRNvfHGG1q2bJmqqqrkdruD6qSmpqq4uDiorLi4WKmpqSH1j6dOAABwgGlbYR9Nddlll2nv3r3as2dP4Bg7dqyuv/567dmz56wgQ5IyMzOVn58fVLZlyxZlZmaG9D0Z0QAAwAGWbFlq/oOfodRNSkrS8OHDg8o6duyobt26BcpnzJih3r17Ky8vT5J0++236+KLL9aiRYs0bdo0rVu3Tjt37tSqVatC6icjGgAAQIWFhTpy5Ejg/YQJE7R27VqtWrVKo0aN0gsvvKANGzacFbB8E9bRAACgFZWVlcnn8+nzA33CTgbtNehTlZaWfmOOhpOYOgEAwAGmbcsM43f9cOq2JqZOAABAxDCiAQCAA1ozGdRJBBoAADjAki2TQAMAAERCrIxokKMBAAAihhENAAAcECtPnRBoAADgAKvuCKd+NGDqBAAARAwjGgAAOMAM86mTcOq2JgINAAAcYNq1Rzj1owFTJwAAIGIY0QAAwAGxkgxKoAEAgAMsGTJlhFU/GjB1AgAAIoYRDQAAHGDZtUc49aMBgQYAAA4ww5w6CaduayLQAADAAbESaJCjAQAAIoYRDQAAHGDZhiw7jKdOwqjbmgg0AABwAFMnAAAAYWJEAwAAB5hyyQzj932zBfsSSQQaAAA4wA4zR8OOkhwNpk4AAEDEMKIBAIADYiUZlEADAAAHmLZLph1GjkaULEHO1AkAAIgYRjQAAHCAJUNWGL/vW4qOIQ0CDQAAHECOBgAAiJjwczSiY0SDHA0AABAxjGgAAOCA2hyNMDZVY+oEAAA0xgpzCfJoSQZl6gQAAEQMgQYAAA6oTwYN5wjFihUrNHLkSCUnJys5OVmZmZl6+eWXGz1/zZo1Mgwj6PB6vSF/T6ZOAABwgCVXq66j0adPHz3yyCMaOHCgbNvW73//e1155ZV69913NWzYsAbrJCcn68CBA4H3hhF6XgiBBgAAMWD69OlB7x966CGtWLFCb731VqOBhmEYSk1NDatdpk4AAHCAaRthH5JUVlYWdFRVVX1z26apdevWqaKiQpmZmY2eV15ern79+iktLU1XXnml3n///ZC/J4EGAAAOMOueOgnnkKS0tDT5fL7AkZeX12ibe/fuVadOneTxeHTrrbdq/fr1Gjp0aIPnDho0SKtXr9bGjRv17LPPyrIsTZgwQZ9++mlI39Ow7ShZWgwAgHagrKxMPp9Pa94dpQ5J7mZf5+QJUzdmvKfDhw8rOTk5UO7xeOTxeBqsU11drcLCQpWWluqFF17Q008/rW3btjUabJzO7/dryJAhuu666/TAAw80uZ/kaAAA4ADLdskKYwlyq26coP4pkqZISEjQgAEDJEljxozRjh079Pjjj+upp576xrrx8fHKyMjQwYMHQ+onUycAADigpaZOwmFZVpNyOqTavI69e/eqZ8+eIbXBiAYAAA6wpEBCZ3Prh2L+/PmaOnWq+vbtqxMnTmjt2rXaunWrXnnlFUnSjBkz1Lt370COx/33368LL7xQAwYMUElJiR577DF98sknuummm0Jql0ADAIAYcPToUc2YMUNHjhyRz+fTyJEj9corr+jyyy+XJBUWFsrl+nqU5KuvvtLNN9+soqIidenSRWPGjNH27dublM9xOpJBAQBoRfXJoCt2j1Nip+b/vn+qvEY/O3+HSktLm5yj4QRGNAAAcEBzlhE/s340iI5eAgCAqMSIBgAADrBkyFI4yaDNr9uaCDQAAHAAUycAAABhYkQDAAAHhLvoVkss2NUaCDQAAHCAZRuywlmwK4y6rSk6wiEAABCVGNEAAMABVphTJ1aUjBUQaAAA4IDwd28l0AAAAI0wZcgMYy2McOq2pugIhwAAQFRiRAMAAAcwdQIAACLGVHjTH2bLdSWioiMcAgAAUYkRDQAAHMDUCQAAiBg2VQMAAAgTIxoAADjAliErjGRQO0rW0SDQAADAAUydAAAAhIkRDQAAHBAr28QTaAAA4AAzzN1bw6nbmgg0AABwQKyMaERHOAQAAKISIxoAADjAkktWGL/vh1O3NRFoAADgANM2ZIYx/RFO3dYUHeEQAACISoxoAADggFhJBiXQAADAAXaYu7farAwKAABiHSMaAAA4wJQhM4yN0cKp25oINAAAcIBlh5dnYdkt2JkIYuoEAABEDCMaAAA4wAozGTScuq0pOnoJAEA7Y8kI+wjFihUrNHLkSCUnJys5OVmZmZl6+eWXz1nn+eef1+DBg+X1ejVixAi99NJLIX9PAg0AABxQvzJoOEco+vTpo0ceeUS7du3Szp07demll+rKK6/U+++/3+D527dv13XXXafZs2fr3XffVXZ2trKzs7Vv376Q2jVs246SdBIAAKJfWVmZfD6f/u9r/1cJnRKafZ3q8mqtvXStSktLlZyc3KxrdO3aVY899phmz5591mfXXHONKioqtGnTpkDZhRdeqNGjR2vlypVNboMRDQAAHFCfoxHOIdUGLqcfVVVV39i2aZpat26dKioqlJmZ2eA5BQUFysrKCiqbMmWKCgoKQvqeBBoAADjAkhFYhrxZR12ORlpamnw+X+DIy8trtM29e/eqU6dO8ng8uvXWW7V+/XoNHTq0wXOLioqUkpISVJaSkqKioqKQvidPnQAAEMUOHz4cNHXi8XgaPXfQoEHas2ePSktL9cILL2jmzJnatm1bo8FGSyDQAADAAXYznhw5s76kwFMkTZGQkKABAwZIksaMGaMdO3bo8ccf11NPPXXWuampqSouLg4qKy4uVmpqakj9ZOoEAAAHhDVtEubOr4E+WFajOR2ZmZnKz88PKtuyZUujOR2NYUQDAIAYMH/+fE2dOlV9+/bViRMntHbtWm3dulWvvPKKJGnGjBnq3bt3IMfj9ttv18UXX6xFixZp2rRpWrdunXbu3KlVq1aF1C6BBgAADmjtlUGPHj2qGTNm6MiRI/L5fBo5cqReeeUVXX755ZKkwsJCuVxfX3PChAlau3at7r77bt11110aOHCgNmzYoOHDh4fULutoAADQiurX0bjyHz9RfMfmr6Phr6jWxsmrw1pHozWQowEAACKGqRMAABzQnP1KzqwfDQg0AABwQLhPjrTEUyetgUADAAAHxEqgQY4GAACIGEY0AABwQKyMaBBoAADggFgJNJg6AQAAEcOIBgAADrAV3iOq0bLaJoEGAAAOYOoEAAAgTIxoAADggFgZ0SDQAADAAbESaDB1AgAAIoYRDQAAHBArIxoEGgAAOMC2DdlhBAvh1G1NBBoAADggVraJJ0cDAABEDCMaAAA4gBwNAAAQMbGSo8HUCQAAiBhGNAAAcABTJwAAIGKYOgEAAAgTIxoAADjADnPqJFpGNAg0AABwgC3JtsOrHw2YOgEAABHDiAYAAA6wZMhgCXIA7cUnn3yin//85xo0aJASExPVrVs3XX311frPf/7jdNeAmFT/1Ek4RzRgRAOIETt27ND27dt17bXXqk+fPvrPf/6jFStWaNKkSfrXv/6lDh06ON1FIKZYtiGDdTQAtBfTpk3TD3/4w6Cy6dOnKzMzUy+++KJuuOEGh3oGoD1j6gSIEYmJiYHXfr9fx44d04ABA9S5c2ft3r3bwZ4Bscm2wz+iAYEGECNOnTqlBQsWKC0tTR6PR927d9e3vvUtlZSUqLS01OnuATGHHA0A7crcuXP1zDPPaN68ecrMzJTP55NhGLr22mtlWZbT3QPQThFoADHihRde0MyZM7Vo0aJAWWVlpUpKSpzrFBDDYmWvEwINIEa43W7ZZ0zqPvHEEzJN06EeAbEtVp46IUcDiBHf+9739Mc//lHz5s3TqlWrNGvWLC1dulTdunVzumsAWkFeXp7GjRunpKQk9ejRQ9nZ2Tpw4MA566xZs0aGYQQdXq83pHYZ0QBixOOPPy63260//elPqqys1MSJE/Xqq69qypQpTncNiEnhPjkSat1t27YpJydH48aNU01Nje666y5NnjxZ//rXv9SxY8dG6yUnJwcFJIYR2kgKgQYQIzp37qzVq1efVc7KoIAzagONcHI0Qjt/8+bNQe/XrFmjHj16aNeuXfrOd77TaD3DMJSamtqcLkpi6gQAgKhWVlYWdFRVVTWpXv1j7V27dj3neeXl5erXr5/S0tJ05ZVX6v333w+pfwQaAAA4oKXW0UhLS5PP5wsceXl539i2ZVmaN2+eJk6cqOHDhzd63qBBg7R69Wpt3LhRzz77rCzL0oQJE/Tpp582+XsydQIAgAPsuiOc+pJ0+PBhJScnB8o9Hs831s3JydG+ffv05ptvnvO8zMxMZWZmBt5PmDBBQ4YM0VNPPaUHHnigSf2MihGN5cuXq3///vJ6vRo/frzeeecdp7sEAEBYWmpEIzk5Oej4pkBjzpw52rRpk15//XX16dMnpD7Hx8crIyNDBw8ebHKdNh9oPPfcc8rNzdXChQu1e/dujRo1SlOmTNHRo0ed7hoAAFHDtm3NmTNH69ev12uvvab09PSQr2Gapvbu3auePXs2uU6bDzR+85vf6Oabb9asWbM0dOhQrVy5Uh06dGgwex4AgKhht8ARgpycHD377LNau3atkpKSVFRUpKKiIp06dSpwzowZMzR//vzA+/vvv1//+Mc/9PHHH2v37t368Y9/rE8++UQ33XRTk9tt0zka1dXV2rVrV9CXdrlcysrKUkFBQYN1qqqqgjJuLcvS8ePH1a1bt5Cf/QUAxBbbtnXixAn16tVLLleEfxcPd2O0EOuuWLFCkjRp0qSg8meeeUY33nijJKmwsDDoe3/11Ve6+eabVVRUpC5dumjMmDHavn27hg4d2uR223Sg8eWXX8o0TaWkpASVp6SkaP/+/Q3WycvL03333dca3QMAtFOHDx8OOX+hrTtzC4KGbN26Nej94sWLtXjx4rDabdOBRnPMnz9fubm5gfelpaXq27evvuP5geLivDK8XpVefJ5uuGuTxnk/0fvVPbWtZJD2fNFbJcVJij8WJ89xQwmlthJOWEo4YSnuZI3clX65TtVI1X4ZNTWSv0Z2jSmZpmRZkmnW/ke0rNpVVM7cDdPtluHxqHJUP13+yJua1mmfPqzurq1lQ7TrizQdPepT/BfxSvjKkKfEVvwJWwknTMWfrJH7VI2MUzUyqv0y/P7Q23a5ZHg8qhnYR1lPFOjmzv/Wf/x+5Z8crG3Hvq2Dxd3l/yJRCcfc8pRInlJb8RWW4itMxVXUyFVVI1flad+9xpRqTNmmJVlm3aozlmTZQd9X8fGyU7rq4PXJ+umlr+rbCUX6vKaz9lX00fslqTpyPFn+Eq/iS92KO2EooVyKL7flrrQVV2nJXVV7uKpNuWosqcaWYVmSacs44y+MXT9adeYvIJZqz7UsybQky6q9Rv1705Jddw+D7mXde9tW7feyLdn1388+7f6Gs6wfgDanRn69qZeUlJQU8bZae2VQp7TpQKN79+5yu90qLi4OKi8uLm50lTKPx9Ngxm2cEa84I0GGK0Fx8V4ldopTp0SXOlS5lVCTIHeFR65Er1zeOLk9htwJtuLirdojrkZut1sut19yu2RYbsnllu2qkWyz9h8ew5QtWzIsKfDzNIa7tu04r7yd4tQpyaUO1W4lWPFyn6xvO/60tm3FxZt1bdfIcPtluF0yTFcz2nbJMBKkuraTk1zq5HfJ64pTfGWCXB28ciV65fa65faotv1qS3Fxte27avxyud2nfXdTctXIti1JZl2AYUmGHfR95YqX7fbI5a1tt4PHrUR/nBKMeMX5PXJVeuWq8spV5Za72pC7uq5t01acacldd7gsUy7bkmxLRt13bHKgYdQFGoZV10dLhupe25Zkm7JlNXAvTcmw6l7Xfj878P1Ov79R8jcdQNPU/ZVujal2dm9tAxISEjRmzBjl5+crOztbUm3ORX5+vubMmdMibZitnA9r2JLfipNpG/LbcbJsV+0OfPX/XtUl+Bj1iT6Wvg5bwwlfTVNxx8q1bNckfTC4p76qTtRHx7ur5GiS4r+MU8djhjxf2fKWWoo/YSq+vEauk365KqtlVPmlar9UUyO7pqb2p2mdMZpiB37Tt08f2Tj6hQbkSi+ri162O9cV+pWgTxRKvrN9xk8AQHRo04GGJOXm5mrmzJkaO3asLrjgAi1ZskQVFRWaNWuW010DAKD5bCPkhM6z6keBNh9oXHPNNfriiy+0YMECFRUVafTo0dq8efNZCaIAAEQTcjTakDlz5rTYVAkAAG1CS61B3sZFRaARadZpw09GhP/D2YYU76qR27DlPjNp80yGapMb65OSwklOcrtV062T7hibr7ldPtEhf7m2dP+2Xu02RPuKeqoiuaP8SW75v3LLU+JSwgm34iviFVeRIPcpv4zK+qdeap86MWpqAk9rGGfkaMhw1SaOxsdLiV6Z3ZJUdl5HHRtmyN+3Sh2SquRyWfL741RdGSfrZJxcp9xyVxpyVxpyVUkuf91h6uvcS5dkuSUrQTITJNNry/JIpteSnWBLcZbktgO3ybaM2rqWUXdIRn2ZbXydB3Pa/Q78UTDs2vtv1Jef+w+GceYQ5rlOj47RTiAmWacqpV9sdLob7QqBBgAADuCpEwAAEFlRMv0Rjja/1wkAAIhejGi0Npchd1ND2DNPCys92ZLLb2pfRS+933G/PvL30IGTqfqs3KdTJzyKK3cp7qQUd9JWXP3KnJWmXFU1MqpPy8/w13y9noZZu+hVYE2N01YlPX1tDfs/tpJ2SUmsqAmgjaux/TrcSm0xdQIAACInRp46YeoEAABEDCMaAAA4ov4Z+nDqt30EGjHEqLFk2S59y23Jb3+l3p4SdUs8qWMdO6q6Y5zMU27VVNZubuaqcdXuK2baX68tYhi1Gw25XTLc7sA+J4ZZt4NrfU6GadWWGZZs25AhU5JkW67T1towyNMAENuYOgEAAAhPs0Y0Kioq9Mgjjyg/P19Hjx6VZQWvcPnxxx+3SOcAAGi3YmREo1mBxk033aRt27bphhtuUM+ePWuH0wEAQNOxe2vjXn75Zf3973/XxIkTW7o/7V79nwvTNmTarTtzZZSf0qt7h2hS5/1Kcp2S33bL6/bLk1CjKq+pmg4uGX5DhmXINlyy3IZMjyF3pVvuqni5/JYM06rNxZBku12yPG75k+JUkepW2f+R1P+k+vU4rm95y+UybFWacSqtTtSJKo/KKz2qqoyXWe2WXe2S4XdJNUZdLoghw9LXEXr9HiMuW7ZLkqv2tdx1Pw0F70ei0+ra+vpG23X7mwR91tDNafjnufY4OWt/k9P70Gilb/gcgKNac6+TWNm9tVn/0nXp0kVdu3Zt6b4AAIB2plmBxgMPPKAFCxbo5MmTLd0fAABig90CRxRo1tTJokWL9NFHHyklJUX9+/dXfHx80Oe7d+9ukc4BANBukaPRuOzs7BbuBiLN9tfIOHpMQ/5fhda6xsm2LKmmRjJr1MtfqF51a2IE1sKoy8MIrHvRwGSgIcldd3gldTvts2OnvXZJ8tUdANCWteZeJ7GiWYHGwoULW7ofAADEFMOWzpFv3qT60SCslUF37dqlDz74QJI0bNgwZWRktEinAABo91hHo3FHjx7Vtddeq61bt6pz586SpJKSEl1yySVat26dvvWtb7VkHwEAQJRq1lMnc+fO1YkTJ/T+++/r+PHjOn78uPbt26eysjLddtttLd3H9iWcxc2aUde2bdmmKbvaL+tUpaxjx2V+8aXML47J/KpUZmmZrJMnZVVWyq6qku2vll1TI1lm7RHug94AgIbVJ4OGc0SBZo1obN68Wa+++qqGDBkSKBs6dKiWL1+uyZMnt1jnAABot2Jk6qRZIxqWZZ31SKskxcfHn7XvCQAAiF3NCjQuvfRS3X777fr8888DZZ999pnuuOMOXXbZZS3WOQAA2q0YWbCrWYHGsmXLVFZWpv79++u8887Teeedp/T0dJWVlemJJ55o6T62Oy7Dktuw5TYsuQxLrtOfUarf4+P0PTykr/9LGYZkGLUb2Rmu+gueVt+QXLXlgXUxTLN2bQzTDByyLXIwAMBJMRJoNCtHIy0tTbt379arr76q/fv3S5KGDBmirKysFu0cAADtFiuDnpthGLr88st1+eWXt2R/AABAO9LkQGPp0qW65ZZb5PV6tXTp0nOeyyOuAACcGyuDnmHx4sW6/vrr5fV6tXjx4kbPMwyDQOMcjLo9RJIMW13d5eoSf1IdE6rl8pgyE+JkxhsyEwxZ8basBENmtUtGglsuy5Jh2rJVl7phGDKMGtkuQzJNGe7T9ikxrdrPTVO2S5Ll+nrPEuO01wAA5/B4a7BDhw6pW7dugdeNHR9//HGTG3/jjTc0ffp09erVS4ZhaMOGDUGf27atBQsWqGfPnkpMTFRWVpb+/e9/N/n6AACgVl5ensaNG6ekpCT16NFD2dnZOnDgwDfWe/755zV48GB5vV6NGDFCL730UkjtNuupk/vvv18nT548q/zUqVO6//77m3ydiooKjRo1SsuXL2/w81//+tdaunSpVq5cqbffflsdO3bUlClTVFlZ2ZxuAwAQs7Zt26acnBy99dZb2rJli/x+vyZPnqyKiopG62zfvl3XXXedZs+erXfffVfZ2dnKzs7Wvn37mtxuswKN++67T+Xl5WeVnzx5Uvfdd1+TrzN16lQ9+OCD+sEPfnDWZ7Zta8mSJbr77rt15ZVXauTIkfrDH/6gzz///KyRDwAAoo2hr/M0mnWE2N7mzZt14403atiwYRo1apTWrFmjwsJC7dq1q9E6jz/+uL773e/ql7/8pYYMGaIHHnhA559/vpYtW9bkdpv11Ilt27XrOJzhvffeU9euXZtzybMcOnRIRUVFQY/M+nw+jR8/XgUFBbr22msbrFdVVaWqqqrA+7KyshbpT4uwbCUcr9Tm4mEal/ix/LZbHVzV6uI5qaMdO6k8KU5+vyHDql0jw4ozFJ/gknnKLXdVnFzVllx+S4Zp1eZjSLLjXDK9caruHKfynm6d+D+WPP1PqF/Xr9Q54ZQk6WRNgk74PSqr9OpkVbyqquJlVrmlapeMGpdkqrbN+vm++ieuXLbkkmyXXfu6vqyxP9319U9/5CpK5hABQJKsU5XSbRud7kZIzvx3zuPxyOPxfGO90tJSSTrnv9sFBQXKzc0NKpsyZUpIv/CHNKLRpUsXde3aVYZh6Nvf/ra6du0aOHw+ny6//HL96Ec/CuWSjSoqKpIkpaSkBJWnpKQEPmtIXl6efD5f4EhLS2uR/gAA0KJaaFO1tLS0oH/38vLyvrFpy7I0b948TZw4UcOHD2/0vKKiopD/HT5TSCMaS5YskW3b+slPfqL77rtPPp8v8FlCQoL69++vzMzMUC7Z4ubPnx8UfZWVlRFsAADanhZ66uTw4cNKTk4OFDdlNCMnJ0f79u3Tm2++GUYHmiakQGPmzJmSpPT0dE2YMKHBjdVaSmpqqiSpuLhYPXv2DJQXFxdr9OjRjdZr6pARAADtQXJyclCg8U3mzJmjTZs26Y033lCfPn3OeW5qaqqKi4uDyoqLiwP/RjdFkwONsrKywBfJyMjQqVOndOrUqQbPDeULNyY9PV2pqanKz88PBBZlZWV6++239bOf/Szs67c625ZdUyPXf47IuClJD3mukUxLRrVf8tcorfoLqeaI7Jqa2r1JrLr1MOrWvLDr1sg4/XpSbbpEXN3RQVKPuo9NScdOaz5OUte6AwDQsBrbr09bq7FWXkfDtm3NnTtX69ev19atW5Wenv6NdTIzM5Wfn6958+YFyrZs2RLS7EWTA40uXbroyJEj6tGjhzp37txgMmh9kqhpmk26Znl5uQ4ePBh4f+jQIe3Zs0ddu3ZV3759NW/ePD344IMaOHCg0tPTdc8996hXr17Kzs5uarcBAGiTWntl0JycHK1du1YbN25UUlJSIM/C5/MpMTFRkjRjxgz17t07kOdx++236+KLL9aiRYs0bdo0rVu3Tjt37tSqVaua3G6TA43XXnstkJn6+uuvN7mBc9m5c6cuueSSwPv63IqZM2dqzZo1+tWvfqWKigrdcsstKikp0UUXXaTNmzfL6/W2SPsAADimlUc0VqxYIUmaNGlSUPkzzzyjG2+8UZJUWFgol+vr50QmTJigtWvX6u6779Zdd92lgQMHasOGDedMID1TkwONiy++uMHX4Zg0aVLtVuaNMAxD999/f0iLgAEAgLOd69/belu3bj2r7Oqrr9bVV1/d7HabtY7G5s2b1alTJ1100UWSpOXLl+u3v/2thg4dquXLl6tLly7N7lBrc8uS67Txp3PtumsbRtMXSHG5JOu0nArLkm2asqurpdIT35h7AQBo59jrpHG//OUvAwuE7N27V7m5ubriiit06NChsxb2AAAAZwtrVdAw8ztaU7NGNA4dOqShQ4dKkl588UVNnz5dDz/8sHbv3q0rrriiRTsIAACiV7NGNBISEgKbqr366quaPHmypNplTNvUkt8AALRVLbQyaFvXrBGNiy66SLm5uZo4caLeeecdPffcc5KkDz/88BsX/2gTbFtuw5JLtlyGJZdhyX3afh62S3W73Ri1r121rwM/6w7DMGQbrtp6tiHDkuz63AzbDuRp2LZduzcJeRkAgHrkaDRu2bJliouL0wsvvKAVK1aod+/ekqSXX35Z3/3ud1u0gwAAIHo1a0Sjb9++2rRp01nlixcvDrtDAADEgtZesMspzQo0JMk0TW3YsEEffPCBJGnYsGH6/ve/L7fb3WKdAwCg3YqRqZNmBRoHDx7UFVdcoc8++0yDBg2SVLs9e1pamv7+97/rvPPOa9FOtjRXTe3Pbm5bqXGl6uUp1ceJ3VXSoYNqOsTJfcqQu9KQyy8ZpkuGFSdZksuQXIYhuV2S2yWjxi3DtGpzMurzMk7PtzDN2v1MTFO2aUqmWbdMe21Oh23Vn1uXr0GuBgCgnWlWjsZtt92m8847T4cPH9bu3bu1e/duFRYWKj09XbfddltL9xEAgPYn3DU0ouR302aNaGzbtk1vvfVWYO8TSerWrZseeeQRTZw4scU6BwBAu8XUSeM8Ho9OnDhxVnl5ebkSEhLC7hQAAO0egUbjvve97+mWW27R7373O11wwQWSpLffflu33nqrvv/977doB1ucbavjpyf1wM5p2jvkPbkMW4cquulYRQfVVLnlrqkdkrLdkhUnWfGGTI8hw3JLLsl2u+SKd0tmgoz6vUwMQ3a8WzWd4nUyJUElA12qHnJS3+55VL07lMplWKqo8aikOllfVSbqRKVHlacSVFMVJ7vKJcPvklFjyDBV9wfPkAw7sIaH7bIlt2S768bLXKobO5MCm6/Yp/20jeD3AIAmsU5VSrdvdLob7UqzcjSWLl2qAQMGaMKECfJ6vfJ6vZo4caIGDBigxx9/vKX7CABAu8NeJw2wLEuPPfaY/vrXv6q6ulrZ2dmaOXOmDMPQkCFDNGDAgEj1EwAARKGQAo2HHnpI9957r7KyspSYmKiXXnpJPp9Pq1evjlT/AABAFAsp0PjDH/6gJ598Uj/96U8l1W6oNm3aND399NNyuZo1C9Pq7Gq/XP8u1KA50gdWp7r1LcrVwz6hHnV7kATWt7AbXt+isfQHt6SkukOSTEmFgU9PSSpRJ0mdWvILAQBaTI3t16et1ViMJIOGFB0UFhYGbQOflZUlwzD0+eeft3jHAABoz2IlRyOkQKOmpkZerzeoLD4+Xn6/v0U7BQAA2oeQpk5s29aNN94oj8cTKKusrNStt96qjh07Bsr+8pe/tFwPAQBor6JkVCIcIQUaM2fOPKvsxz/+cYt1JuJMU3ZNTW0ORkO5GBL7jQAAWkeM5GiEFGg888wzkeoHAABoh5q9TTwAAGi+cBM6oyUZlEADAAAnMHXSvti2ZFsWuRkAgDYhVkY0omOVLQAAEJViZkQDAIA2hakTAAAQMQQakZeXl6e//OUv2r9/vxITEzVhwgQ9+uijGjRoUOCcyspK/eIXv9C6detUVVWlKVOm6Mknn1RKSkpIba3e84ZSfB69W+XStorB+ufx8/TRl9108ssOiiuJU0KpobgKKb7CVtwpW3FVttxVtowaWy7TlmHW/he1DcmOc8n0GvInulSdbKjaZ6g62VZNB1t2vF07IWVJht+Qy2/I5Zdc1XU//ZKrpvYwTFuGJRmnp4oYku2SbJch2y3ZcZLlVu1rlySj9hzptPm5+uVorbrX1tdlDf5BrL+G8XWb9a8DPwMdavh1U+YG7TOvBQBtnFlVKS3d6HQ32hVHczS2bdumnJwcvfXWW9qyZYv8fr8mT56sioqKwDl33HGH/va3v+n555/Xtm3b9Pnnn+u//uu/HOw1AADhi5W9Thwd0di8eXPQ+zVr1qhHjx7atWuXvvOd76i0tFS/+93vtHbtWl166aWSahcNGzJkiN566y1deOGFTnQbAIDwMXXS+kpLSyVJXbt2lSTt2rVLfr9fWVlZgXMGDx6svn37qqCgoMFAo6qqSlVVVWdd80S5pUSXpYoqqfKkX/6Kapknq2SdcsmqjJNZZcioklzVtoxqW/Lbsv11UyfWGVMntkum25AZ55JZbcisMmRV2rJctuya06ZOagzJb0h+yfYbsv2qPWok22zC1Imr9jymTgCgdZhVlZJq9/ZCy2gzgYZlWZo3b54mTpyo4cOHS5KKioqUkJCgzp07B52bkpKioqKiBq+Tl5en++6776zykeOOnvbuk5bqNgCgHTp27Jh8Pl9kG2FEo3Xl5ORo3759evPNN8O6zvz585Wbmxt4X1JSon79+qmwsDDyf2jambKyMqWlpenw4cNKTk52ujtRh/vXfNy75uPehae0tFR9+/YNjKxHUqws2NUmAo05c+Zo06ZNeuONN9SnT59AeWpqqqqrq1VSUhI0qlFcXKzU1NQGr+XxeIK2sa/n8/n4S9dMycnJ3LswcP+aj3vXfNy78LhcrGfZUhy9k7Zta86cOVq/fr1ee+01paenB30+ZswYxcfHKz8/P1B24MABFRYWKjMzs7W7CwBAy7Fb4AjBG2+8oenTp6tXr14yDEMbNmw45/lbt26VYRhnHY2lLjTG0RGNnJwcrV27Vhs3blRSUlKg8z6fT4mJifL5fJo9e7Zyc3PVtWtXJScna+7cucrMzOSJEwBAVGvtqZOKigqNGjVKP/nJT0JaJuLAgQNBo2M9evQIqV1HA40VK1ZIkiZNmhRU/swzz+jGG2+UJC1evFgul0tXXXVV0IJdTeXxeLRw4cIGp1Nwbty78HD/mo9713zcu/C06v1r5WTQqVOnaurUqSE306NHj7MeygiFYfMMDwAAraasrEw+n09Dch6W2+Nt9nXMqkp9sPyusxJ/G8tVPJ1hGFq/fr2ys7MbPWfr1q265JJL1K9fP1VVVWn48OG69957NXHixJD6SbYLAABOaKEcjbS0NPl8vsCRl5fXIt3r2bOnVq5cqRdffFEvvvii0tLSNGnSJO3evTuk67SJp04AAIg1p6+V2Nz6khoc0WgJgwYNCtp7bMKECfroo4+0ePFi/fGPf2zydQg0AACIYq35KPMFF1wQ8npXBBoAADghClcG3bNnj3r27BlSHQINAAAc0NqPt5aXl+vgwYOB94cOHdKePXvUtWtX9e3bV/Pnz9dnn32mP/zhD5KkJUuWKD09XcOGDVNlZaWefvppvfbaa/rHP/4RUrvtOhl0+fLl6t+/v7xer8aPH6933nnH6S61SXl5eRo3bpySkpLUo0cPZWdn68CBA0HnVFZWKicnR926dVOnTp101VVXqbi42KEet12PPPKIDMPQvHnzAmXcu8Z99tln+vGPf6xu3bopMTFRI0aM0M6dOwOf27atBQsWqGfPnkpMTFRWVpb+/e9/O9jjtsM0Td1zzz1KT09XYmKizjvvPD3wwANBm4Fx/2p900JVTblPx48f1/XXX6/k5GR17txZs2fPVnl5eSt+i/Dt3LlTGRkZysjIkCTl5uYqIyNDCxYskCQdOXJEhYWFgfOrq6v1i1/8QiNGjNDFF1+s9957T6+++qouu+yykNptt4HGc889p9zcXC1cuFC7d+/WqFGjNGXKFB09evSbK8eYbdu2KScnR2+99Za2bNkiv9+vyZMnq6KiInDOHXfcob/97W96/vnntW3bNn3++echLfgSC3bs2KGnnnpKI0eODCrn3jXsq6++0sSJExUfH6+XX35Z//rXv7Ro0SJ16dIlcM6vf/1rLV26VCtXrtTbb7+tjh07asqUKaqsrHSw523Do48+qhUrVmjZsmX64IMP9Oijj+rXv/61nnjiicA53L9a9QtVLV++vMHPm3Kfrr/+er3//vvasmVLYMuMW265JbyOtfLKoJMmTZJt22cda9askSStWbNGW7duDZz/q1/9SgcPHtSpU6d07Ngxvf7667rkkktC/prtdh2N8ePHa9y4cVq2bJmk2t1h09LSNHfuXN15550O965t++KLL9SjRw9t27ZN3/nOd1RaWqpvfetbWrt2rX74wx9Kkvbv368hQ4aooKCAVVpVOyR5/vnn68knn9SDDz6o0aNHa8mSJdy7c7jzzjv1z3/+U//zP//T4Oe2batXr176xS9+of/+7/+WVLvhVUpKitasWaNrr722Nbvb5nzve99TSkqKfve73wXKrrrqKiUmJurZZ5/l/jXizPUjmnKfPvjgAw0dOlQ7duzQ2LFjJUmbN2/WFVdcoU8//VS9evUKqQ/162gM++nDcieEsY5GdaXef+oulZaWtul9bdrliEZ1dbV27dqlrKysQJnL5VJWVpYKCgoc7Fl0KC0tlaTA7oW7du2S3+8Pup+DBw9W3759uZ91cnJyNG3atKB7JHHvzuWvf/2rxo4dq6uvvlo9evRQRkaGfvvb3wY+P3TokIqKioLunc/n0/jx42P+3km1jxrm5+frww8/lCS99957evPNNwMrP3L/mqYp96mgoECdO3cOBBmSlJWVJZfLpbfffrvV+xxt2mUy6JdffinTNJWSkhJUnpKSov379zvUq+hgWZbmzZuniRMnavjw4ZKkoqIiJSQknLUEbUpKSsib67RH69at0+7du7Vjx46zPuPeNe7jjz/WihUrlJubq7vuuks7duzQbbfdpoSEBM2cOTNwfxr6exzr906qHREqKyvT4MGD5Xa7ZZqmHnroIV1//fWSxP1roqbcp6KiorP294iLi1PXrl3DupdsE4+YlJOTo3379oX8nHSsOnz4sG6//XZt2bJFXm/zh0BjkWVZGjt2rB5++GFJUkZGhvbt26eVK1dq5syZDveu7fvzn/+sP/3pT1q7dq2GDRumPXv2aN68eerVqxf3L1pE4eOtzdEup066d+8ut9t9VmZ/cXGxUlNTHepV2zdnzhxt2rRJr7/+uvr06RMoT01NVXV1tUpKSoLO537WTo0cPXpU559/vuLi4hQXF6dt27Zp6dKliouLU0pKCveuET179tTQoUODyoYMGRLIeq+/P/w9btgvf/lL3Xnnnbr22ms1YsQI3XDDDbrjjjsCy09z/5qmKfcpNTX1rAcJampqdPz48bDuZf2IRjhHNGiXgUZCQoLGjBmj/Pz8QJllWcrPz1dmZqaDPWubbNvWnDlztH79er322mtKT08P+nzMmDGKj48Pup8HDhxQYWFhzN/Pyy67THv37tWePXsCx9ixY3X99dcHXnPvGjZx4sSzHqP+8MMP1a9fP0lSenq6UlNTg+5dWVmZ3n777Zi/d5J08uRJuVzB/wt3u92yLEsS96+pmnKfMjMzVVJSol27dgXOee2112RZlsaPH9/qfY427XbqJDc3VzNnztTYsWN1wQUXaMmSJaqoqNCsWbOc7lqbk5OTo7Vr12rjxo1KSkoKzDn6fD4lJibK5/Np9uzZys3NVdeuXZWcnKy5c+cqMzMzpp+akKSkpKRALku9jh07qlu3boFy7l3D7rjjDk2YMEEPP/ywfvSjH+mdd97RqlWrtGrVKkkKrEfy4IMPauDAgUpPT9c999yjXr16nXPHyVgxffp0PfTQQ+rbt6+GDRumd999V7/5zW/0k5/8RBL373TftFDVN92nIUOG6Lvf/a5uvvlmrVy5Un6/X3PmzNG1114b8hMnQWJk6qTdBhrXXHONvvjiCy1YsEBFRUUaPXq0Nm/efFbCD6QVK1ZIqn3G+nTPPPOMbrzxRknS4sWL5XK5dNVVV6mqqkpTpkzRk08+2co9jU7cu4aNGzdO69ev1/z583X//fcrPT1dS5YsCSQzSrXP8VdUVOiWW25RSUmJLrroIm3evJl8GElPPPGE7rnnHv385z/X0aNH1atXL/30pz8NLL4kcf/q7dy5M2j9h9zcXEnSzJkztWbNmibdpz/96U+aM2eOLrvsssDf56VLl4bVr1hJBm2362gAANAW1a+jMfIn4a+j8b+r2/46Gu12RAMAgDaNqRMAABAxMRJotMunTgAAQNvAiAYAAA6IlWRQAg0AAJzA1AkAAEB4GNEAAMABhm3LCGOFiXDqtiZGNIB26sYbb4y5FSCBqGK3wBEFGNEAopBhGOf8fOHChXr88cfFenxA20UyKIA268iRI4HXzz33nBYsWBC0QVmnTp3UqVMnJ7oGAEGYOgGiUGpqauDw+XwyDCOorFOnTmdNnUyaNElz587VvHnz1KVLF6WkpOi3v/1tYLPBpKQkDRgwQC+//HJQW/v27dPUqVPVqVMnpaSk6IYbbtCXX37Zyt8YaIdiZOqEQAOIIb///e/VvXt3vfPOO5o7d65+9rOf6eqrr9aECRO0e/duTZ48WTfccINOnjwpSSopKdGll16qjIwM7dy5U5s3b1ZxcbF+9KMfOfxNgOhXP3USzhENCDSAGDJq1CjdfffdGjhwoObPny+v16vu3bvr5ptv1sCBA7VgwQIdO3ZM//u//ytJWrZsmTIyMvTwww9r8ODBysjI0OrVq/X666/rww8/dPjbAIgG5GgAMWTkyJGB1263W926ddOIESMCZSkpKZKko0ePSpLee+89vf766w3me3z00Uf69re/HeEeA+1YjCzYRaABxJD4+Pig94ZhBJXVP81iWZYkqby8XNOnT9ejjz561rV69uwZwZ4C7R9PnQCIeeeff75efPFF9e/fX3Fx/O8CQOjI0QDQqJycHB0/flzXXXedduzYoY8++kivvPKKZs2aJdM0ne4eEN146gRArOvVq5f++c9/yjRNTZ48WSNGjNC8efPUuXNnuVz87wMIV3t/4kSSDJulAwEAaDVlZWXy+Xwa86OHFBfvbfZ1avyV2vXn/6fS0lIlJye3YA9bFpOuAAA4wbZrj3DqRwECDQAAHMBTJwAAIHJiZB0NsrkAAEDEMKIBAIADDKv2CKd+NCDQAADACUydAAAAhIdAAwAAB7T2NvFvvPGGpk+frl69eskwDG3YsOEb62zdulXnn3++PB6PBgwYoDVr1oT8PQk0AABwQv06GuEcIaioqNCoUaO0fPnyJp1/6NAhTZs2TZdccon27NmjefPm6aabbtIrr7wSUrvkaAAAEAOmTp2qqVOnNvn8lStXKj09XYsWLZIkDRkyRG+++aYWL16sKVOmNPk6jGgAAOCAlpo6KSsrCzqqqqpapH8FBQXKysoKKpsyZYoKCgpCug6BBgAATmih3VvT0tLk8/kCR15eXot0r6ioSCkpKUFlKSkpKisr06lTp5p8HaZOAACIYocPHw7aVM3j8TjYm7MRaAAA4ICW2uskOTk5Iru3pqamqri4OKisuLhYycnJSkxMbPJ1CDQAAHBCG9+9NTMzUy+99FJQ2ZYtW5SZmRnSdcjRAADAAa29jkZ5ebn27NmjPXv2SKp9fHXPnj0qLCyUJM2fP18zZswInH/rrbfq448/1q9+9Svt379fTz75pP785z/rjjvuCKldAg0AAGLAzp07lZGRoYyMDElSbm6uMjIytGDBAknSkSNHAkGHJKWnp+vvf/+7tmzZolGjRmnRokV6+umnQ3q0VZIM247w2AsAAAgoKyuTz+dT5nfvV1y8t9nXqfFXqmDzApWWlkYkR6OlkKMBAIADWioZtK1j6gQAAEQMIxoAADjBsmuPcOpHAQINAACccNrqns2uHwWYOgEAABHDiAYAAA4wFGYyaIv1JLIINAAAcEIbXxm0pTB1AgAAIoYRDQAAHBAr62gQaAAA4IQYeeqEQAMAAAcYti0jjDyLcOq2JnI0AABAxDCiAQCAE6y6I5z6UYBAAwAABzB1AgAAECZGNAAAcAJPnQAAgIhhZVAAAIDwMKIBAIADWBkUAABEDlMnAAAA4WFEAwAABxhW7RFO/WhAoAEAgBNiZOqEQAMAACfEyDoa5GgAAICIYUQDAAAHxMpeJwQaAAA4IUZyNJg6AQAAEcOIBgAATrAlhfOIanQMaBBoAADghFjJ0WDqBAAARAwjGgAAOMFWmMmgLdaTiCLQAADACTx1AgAAEB5GNAAAcIIlyQizfhRgRAMAAAfUP3USztEcy5cvV//+/eX1ejV+/Hi98847jZ67Zs0aGYYRdHi93pDaI9AAAMAJ9Tka4Rwheu6555Sbm6uFCxdq9+7dGjVqlKZMmaKjR482Wic5OVlHjhwJHJ988klIbRJoAAAQI37zm9/o5ptv1qxZszR06FCtXLlSHTp00OrVqxutYxiGUlNTA0dKSkpIbRJoAADghBYa0SgrKws6qqqqGmyuurpau3btUlZWVqDM5XIpKytLBQUFjXazvLxc/fr1U1pamq688kq9//77IX1NAg0AAJzQQoFGWlqafD5f4MjLy2uwuS+//FKmaZ41IpGSkqKioqIG6wwaNEirV6/Wxo0b9eyzz8qyLE2YMEGffvppk78mT50AABDFDh8+rOTk5MB7j8fTYtfOzMxUZmZm4P2ECRM0ZMgQPfXUU3rggQeadA0CDQAAnNBCj7cmJycHBRqN6d69u9xut4qLi4PKi4uLlZqa2qQm4+PjlZGRoYMHDza5m0ydAADggNZ+vDUhIUFjxoxRfn5+oMyyLOXn5weNWpyLaZrau3evevbs2eR2GdEAACBG5ObmaubMmRo7dqwuuOACLVmyRBUVFZo1a5YkacaMGerdu3cgz+P+++/XhRdeqAEDBqikpESPPfaYPvnkE910001NbpNAAwAAJziw18k111yjL774QgsWLFBRUZFGjx6tzZs3BxJECwsL5XJ9Pdnx1Vdf6eabb1ZRUZG6dOmiMWPGaPv27Ro6dGiT2zRsO0p2ZQEAoB0oKyuTz+dT1nnzFOdufuJmjVmlVz9aotLS0iblaDiFHA0AABAxTJ0AAOCEGNkmnkADAABHhBloiEADAAA0JkZGNMjRAAAAEcOIBgAATrBshTX9YUXHiAaBBgAATrCt2iOc+lGAqRMAABAxjGgAAOCEGEkGJdAAAMAJMZKjwdQJAACIGEY0AABwAlMnAAAgYmyFGWi0WE8iiqkTAAAQMYxoAADgBKZOAABAxFiWpDAW3bKiY8EuAg0AAJwQIyMa5GgAAICIYUQDAAAnxMiIBoEGAABOYGVQAACA8DCiAQCAA2zbkh3GVu/h1G1NBBoAADjBtsOb/oiSHA2mTgAAQMQwogEAgBPsMJNBo2REg0ADAAAnWJZkhJFnESU5GkydAACAiGFEAwAAJzB1AgAAIsW2LNlhTJ3weCsAAGhcjIxokKMBAAAihhENAACcYNmS0f5HNAg0AABwgm1LCufx1ugINJg6AQAAEcOIBgAADrAtW3YYUyc2IxoAAKBRthX+0QzLly9X//795fV6NX78eL3zzjvnPP/555/X4MGD5fV6NWLECL300kshtUegAQBAjHjuueeUm5urhQsXavfu3Ro1apSmTJmio0ePNnj+9u3bdd1112n27Nl69913lZ2drezsbO3bt6/JbRp2tIy9AADQDpSVlcnn82mS8QPFGfHNvk6N7ddWe71KS0uVnJzcpDrjx4/XuHHjtGzZMkmSZVlKS0vT3Llzdeedd551/jXXXKOKigpt2rQpUHbhhRdq9OjRWrlyZZPaZEQDAAAntPLUSXV1tXbt2qWsrKxAmcvlUlZWlgoKChqsU1BQEHS+JE2ZMqXR8xtCMigAAA6okT+shUFr5JdUO0JyOo/HI4/Hc9b5X375pUzTVEpKSlB5SkqK9u/f32AbRUVFDZ5fVFTU5H4SaAAA0IoSEhKUmpqqN4tCS6psSKdOnZSWlhZUtnDhQt17771hX7ulEGgAANCKvF6vDh06pOrq6rCvZdu2DMMIKmtoNEOSunfvLrfbreLi4qDy4uJipaamNlgnNTU1pPMbQqABAEAr83q98nq9rdpmQkKCxowZo/z8fGVnZ0uqTQbNz8/XnDlzGqyTmZmp/Px8zZs3L1C2ZcsWZWZmNrldAg0AAGJEbm6uZs6cqbFjx+qCCy7QkiVLVFFRoVmzZkmSZsyYod69eysvL0+SdPvtt+viiy/WokWLNG3aNK1bt047d+7UqlWrmtwmgQYAADHimmuu0RdffKEFCxaoqKhIo0eP1ubNmwMJn4WFhXK5vn4gdcKECVq7dq3uvvtu3XXXXRo4cKA2bNig4cOHN7lN1tEAAAARwzoaAAAgYgg0AABAxBBoAACAiCHQAAAAEUOgAQAAIoZAAwAARAyBBgAAiBgCDQAAEDEEGgAAIGIINAAAQMQQaAAAgIgh0AAAABHz/wEyzDKOjmfuwAAAAABJRU5ErkJggg==",
|
|
"text/plain": [
|
|
"<Figure size 640x480 with 2 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAg0AAAGFCAYAAAB+E8C8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA16klEQVR4nO3de3hU1b3/8c/eezKTAEm4CAmpQXIUFVQQxXIQT+sFsXirbY/WHvRE7amtJ6CYXn5yngpqlWj7tFIqovWo2Md6tDdsj630Uap4bEURqvWuKBaqAiImIYHMZe/1+2MmQ4YE3JkJszOZ9+t59pPMntl7rWxuX9b3u9ayjDFGAAAAn8AOugMAAKAwEDQAAABfCBoAAIAvBA0AAMAXggYAAOALQQMAAPCFoAEAAPgSCroDAAAUm46ODsVisZzvEw6HVVpa2gc98oegAQCAPOro6FDdIUO0ZZub872qq6u1cePGvAUOBA0AAORRLBbTlm2u/r5urCrKs68SaN3p6ZDj31UsFiNoAABgIBtSbmlIuZX19Z6yvzZbBA0AAATANZ7cHHZ/co3Xd53xidkTAADAF0YaAAAIgCcjT9kPNeRybbYIGgAACIAnT7kkGHK7OjukJwAAgC+MNAAAEADXGLkm+xRDLtdmi6ABAIAAFGJNA+kJAADgCyMNAAAEwJORW2AjDQQNAAAEoBDTEwQNAAAEoBALIalpAAAAvjDSAABAALzUkcv1+UbQAABAANwcCyFzuTZbpCcAAIAvjDQAABAA1yjHrbH7ri9+ETQAABCAQqxpID0BAAB8YaQBAIAAeLLkysrp+nxjpAEAgAB4JvejN8aOHSvLsrodDQ0Nvu/BSAMAAEVg7dq1cl03/frll1/W6aefrvPPP9/3PQgaAAAIgJtjeqK3144cOTLj9c0336xDDz1Un/3sZ33fg6ABAIAA9FXQ0NramnE+EokoEons99pYLKb7779fjY2Nsiz/faCmAQCAAHjGyvmQpNraWlVWVqaPpqamT2z74YcfVnNzsy655JJe9ZmRBgAACtjmzZtVUVGRfv1JowySdPfdd2vWrFmqqanpVVsEDQAABKCv0hMVFRUZQcMn+fvf/67HH39cv/nNb3rdJkEDAAABcGXLzaFKwP3kj/To3nvv1ahRo3TWWWf1+lpqGgAAKBKe5+nee+9VfX29QqHejxsw0gAAQABMl2LGbK/vrccff1ybNm3SZZddllWbBA0AAAQg3+s0SNLMmTNlTPbbY5KeAAAAvjDSAABAAFxjyzU5FEJmP2CQNYIGAAAC4MmSl8OAv6f8Rw0EDQAABCCImoZcUdMAAAB8YaQBAIAA5F7TQHoCAICikKxpyD7FkMu12SI9AQAAfGGkAQCAAHg57j3B7AkAAIpEIdY0kJ4AAAC+MNIAAEAAPNks7gQAAD6Zayy5Oexymcu12SI9AQAAfGGkAQCAALg5zp5wSU8AAFAcPGPLy2H2hMeKkAAAFIdCHGmgpgEAAPjCSAMAAAHwlNsMCK/vuuIbQQMAAAHIfZ2G/CcLSE8AAABfGGkAACAAue89kf//9xM0AAAQAE+WPOVS08CKkAAAoJ9ipAEAgACQngAAAL7kvrgTsycAAEA/xUgDAAAB8IwlL5fFnQLYGpugAQCAAHg5pieCWNyJoAEAgADkvsslNQ0AAKCfYqQBAIAAuLLk5rBAUy7XZougAQCAAJCeAAAAAxYjDQAABMBVbikGt++64htBAwAAASA9AQAABixGGgAACEAhbljFSAMAAAEwsuTlcJgs6iHee+89XXTRRRoxYoTKysp0zDHH6Pnnn/d9PSMNAAAUgY8//ljTp0/XKaecokcffVQjR47UW2+9pWHDhvm+B0EDAAAByHd64pZbblFtba3uvffe9Lm6urpe3YP0BAAAAejc5TKXQ5JaW1szjmg02mN7v/vd7zRlyhSdf/75GjVqlCZPnqy77rqrV30maAAAIABuapfLXA5Jqq2tVWVlZfpoamrqsb133nlHy5Yt07hx4/THP/5RV1xxha688krdd999vvtMegIAgAK2efNmVVRUpF9HIpEeP+d5nqZMmaJFixZJkiZPnqyXX35Zd9xxh+rr6321RdAAAEAAuqYYsr1ekioqKjKChn0ZPXq0JkyYkHFu/Pjx+vWvf+27TYIGAAAC4MmWl0OVQG+vnT59ut54442Mc2+++aYOOeQQ3/egpgEAgCJw9dVXa82aNVq0aJE2bNigBx54QD/96U/V0NDg+x6MNAAAEADXWHJzSE/09toTTjhBK1as0Pz583XDDTeorq5Oixcv1uzZs33fg6ABAIAA9FVNQ2+cffbZOvvss7Nuk/QEAADwhZEGAAACYHLcGtsEsGEVQQMAAAFwZcnNYtOprtfnG+kJAADgCyMNAAAEwDPZFTN2vT7fCBoAAAiAl2NNQy7XZougAQCAAHiy5OVQl5DLtdmipgEAAPjCSAMAAAHI94qQfYGgAQCAABRiTQPpCQAA4AsjDQAABMBTjntPBFAISdAAAEAATI6zJwyzJwAAQH/FSAMAAAEIYmvsXBE0AAAQAGZPAACAAYuRBgAAAkB6AgAA+FKIe08QNAAAEIBCHGmgpgEAAPjCSAMAAAEoxJEGggYAAAJQiEED6QkAAOALIw0AAASgEEcaCBoAAAiAUW7TJk3fdcU30hMAAMAXRhoAAAgA6QkAAOBLIQYNpCcAAIAvjDQAABCAQhxpIGgAACAABA0AAMAXYyyZHP7hz+XabFHTAAAAfGGkAQCAAHiyclrcKZdrs0XQAABAAAqxpoH0BAAA8IWRBgAAAlCIhZAEDQAABID0BAAA6Jeuu+46WZaVcRx55JG9ugcjDQAABCCI9MRRRx2lxx9/PP06FOpdGEDQAABAAEyO6YlsgoZQKKTq6uqs2yQ9ARSRzuHJ7du3B90VAH2ktbU144hGo/v87FtvvaWamhr90z/9k2bPnq1Nmzb1qi2CBgAAAmAkGZPDkbpPbW2tKisr00dTU1OP7U2dOlXLly/XypUrtWzZMm3cuFH/8i//op07d/ruM+kJAAAC4MmS1QcrQm7evFkVFRXp85FIpMfPz5o1K/39xIkTNXXqVB1yyCH6xS9+oa9+9au+2iRoAAAgAH1VCFlRUZERNPg1dOhQHX744dqwYYPva0hPAEVo+/btuuCCC1RRUaERI0boqquuUkdHR9DdApBHbW1tevvttzV69Gjf1xA0AEXoggsuUEdHh5qamnTmmWdqyZIluvzyy4PuFlBUOhd3yuXojW9961tavXq13n33Xf3lL3/RF77wBTmOo6985Su+70F6AihCdXV1+u1vfytJamhoUEVFhW6//XZ961vf0sSJEwPuHVAcOgsac7m+N/7xj3/oK1/5ij766CONHDlSJ510ktasWaORI0f6vgdBA1CEGhoaMl7PnTtXt99+u/7whz8QNAAD1IMPPpjzPQgagCI0bty4jNeHHnqobNvWu+++G0yHgCLEhlUACpJl5f8vH6DYFWLQQCEkUITeeuutjNcbNmyQ53kaO3ZsMB0CUBAIGoAitHTp0ozXP/nJTyRlLv4C4MDK9+yJvkB6AihCGzdu1LnnnqvPfe5zeuaZZ3T//ffr3/7t3zRp0qSguwYUjXzPnugLjDQAReihhx5SJBLRNddco9///veaM2eO7r777qC7BaCfs4wJIlYBAKA4tba2qrKyUuPuv0bOoNKs7+Pu6tBbF92slpaWrJaRzgbpCQAAAlCIsycIGgAACIDRnu2ts70+36hpAAAAvhRE0LB06VKNHTtWpaWlmjp1qp577rmguwQAQE460xO5HPnW74OGhx56SI2NjVq4cKHWr1+vSZMm6YwzztC2bduC7hoAANkzfXDkWb8PGn70ox/pa1/7mi699FJNmDBBd9xxhwYNGqR77rkn6K4BAFBU+nUhZCwW07p16zR//vz0Odu2NWPGDD3zzDM9XhONRhWNRtOvPc/Tjh07NGLECNbXBwDslzFGO3fuVE1NjWz7AP+/OtcUA7MnMm3fvl2u66qqqirjfFVVlV5//fUer2lqatL111+fj+4BAAaozZs36+CDDz6gbRTiipD9OmjIxvz589XY2Jh+3dLSojFjxugknamQSpInLUuJz0zSu+dbumzK0/rXir+pwnb097itNR3/pGeaD9Wr26q1a+sgRbaFVPqxFPnYU6TFVcnOhJz2uOxdUVmxmBSNy8TjUjwu47qS68q4RjKejJf8mvEra1nSlAl6+8Iy/du0v2j20HUaYYf1vhvX2o5aPdV8hF76sEYtW8oV3u4o8pGl0o9Nsu22hELtMdntMVnRmNQRS7btJmTiiU9uW5I98UhtmF2uL//LX1Q/9HmNcsq0xd2ttR01Wt1yhP764cH6eGuFwh+GFPnIUiTVdnhnQqH2eLLtjqgUi8vEElIitqdto+TXzrZTnKGVapt+qLb8s6NBh7aopqJFnrHU0lGmj1sHKd4SUUmzo5JWS5EWo3CbUajdU0m7K2dXQk5HQlYsIcXishIJKeFKridjjNS1LSv1vwLbSo4qdY4s2amvXpdnYYzS65p1vYfJ/AyA4pIwMa1u/h+Vl5cH3ZV+qV8HDQcddJAcx9HWrVszzm/dulXV1dU9XhOJRBSJRLqdD6lEIWtP0KBQqewyS6VDSlRebqvCtjUkbqssFFJJPCxnUER2Wamc0pCcsOSEPYVKXIVCCTmOLduRLNuSbEvGTv4DZSxXslwZy5PkyVhGkqeMapV026WKdGm71U22HU6E5bQn27ZLHTkRS07YpNsOObZsx0q1rWTbni1j2Z/ctiTbicgu7dK2Y6vdtTWoxFHY7dp2qIe2nVTbkmxbxrYlS3valpEsu0vbSY4VVqgk9fMM6lBocESeseQ4EdmJUtmxiOwOR0402Z5TYhQq8RQKuXJCCTlOXJbjSI4ty3Mk25WMm/pH30/Q0DnEuKdPGUGD9hE0eAQNQLHKRzqbxZ36WDgc1vHHH69Vq1bpvPPOk5SsUVi1apXmzJmT/Y2t7nkqb6//VR7Q/2R6RvKkuHHkGsmTJ9dY8lJ1qZ07l1k99SHHjtlbd6jyjUr94lPHadehYR1Z9oHixtFbu6v07s4R2tleKqvDlh21ZMclO2FkuZLlGcnrMnLheRmjCemfqwduc7MGP/aKjnjxIMVHD1XbyIMVK7eVKLVU7kiWJ9lxyYl6CkWNnA5Pod2u7KgrO5aQEp7kdgYGVjIIsCyl/7hYyeBNti3j2Mmvnee6Bg+dfe8coej8vut5P7L9y6Sn+1NnA/QvblT6OE9tGSu3ugSChu4aGxtVX1+vKVOm6NOf/rQWL16s9vZ2XXrppUF3DQCAotLvg4Yvf/nL+vDDD7VgwQJt2bJFxx57rFauXNmtOBIAgEJCIeQBMmfOnNzSEQAA9DcFuPlEQQQNA06qiLHEcuVYki1bjuXmpWmvariajzSqP3ydrhi+VqOcwfog0aa/hHZqtxfWjt2DtH1nRF7ElheyZBxLxlG62LN7ceFeP5ebLBA1XvJ9y7Ykx0kWMkoyjq34IFvRobYSgyTPSdZu2DHJ6bAU6jAKdVhyS205HY6caEh2wpMV91J1FSZ9H+NYMiW23IgjN2LJK7HklljJfttKFWmmakNMsnbC8kzqa+d5k36/m850Y+pnTqcP95dG7KlsoYdzAaQiAfiQiHdIb+enrUIshOz3K0ICAID+gZEGAACCUmAzuwkaAAAIQCGmJwgaAmB6MTe/x7UacmDFEgrttvRhrFwfurYi1m7tNJbavbCibkgJ15ZcK7U2Q9caAJNZ6ut5+2jAlmxLdii1noLjyHKSX82u3Qpv3KZhWyIyZWGZSIm8EjtZL6HOtoyshCcr4UkJT5bnJVeA9LzufZCSi2ql1mZQ5xoNzp76i72ftWV6WJ/B2+u9Lrr9WpHQAwa0hBv95A8VMYIGAACCwOwJAADgj6X9T8fyc31+MdgKAAB8YaShn3H3U9jSU8691xKurITU7oa10wtrp7dbzV5ELe5gtcRLFY2HZCWs1KFkbYNrlNokI7Uj5F79SK3ZYEmySkJSJCKrNCJTGpGJhOSVOD2v6yDJjnfZz8KYZFupeoP0z+sk6xXM3ntbdO5e2bV+wVbGfg7dntl+9pn4xFoTQmwAfYn0BAAA8KVYgob29nbdfPPNWrVqlbZt2yZvr0r6d955p086BwAA+o+sgob/+I//0OrVq3XxxRdr9OjRedl3HACAAaVYtsZ+9NFH9fvf/17Tp0/v6/4UnWSpgNf7WoZs6xs++liVb43SkzVH6MOOIfrUoGbtjJfqH21D9cFHlXI/iiiy3VGkWQrvNCrZ5cnp8GTHXFnx5HoJ6bYtWwqFZIUkKxSSGTJIiVEVav9UqXaNdBSrlNxSydhGtmvJiktOXLKjkhMzyf0m4pKdMLITqXUa3NTPm9orIv3j2p37X+z53qS+V9c/d51fO7to9rxO70GRPrefZ7j3fhNd7723/d1m75IK4mugX0vEO6SX89NW0exyOWzYMA0fPryv+wIAQPEowJqGrOrBv/e972nBggXatWtXX/cHAAD0U1mNNPzwhz/U22+/raqqKo0dO1YlJSUZ769fv75POgcAwIBVLDUN5513Xh93A5LkmuTAT48jTvsbhtp7/YKuLCsj8eXu+FgHrXhFB/15hGKjRumNyoPlltoytlTtSU7Mkx1NyOlw5URdWdG4rI64rHhCisVlEgmZREJyXRnXk1w32T3Pk1p3KrTtI1W+Xaqhnes0lJbIhFJ7QliddQXJdR8sz+vy1duT4DNmzz4T+/qZpPSeFel1GfwU5Ha9p9+EIIW+QNHI594Tlule99Tb6/Mtq6Bh4cKFfd0PAACQJzfffLPmz5+vq666SosXL/Z9XU6LO61bt06vvfaaJOmoo47S5MmTc7kdAADFI6BCyLVr1+rOO+/UxIkTe31tVkHDtm3bdOGFF+rJJ5/U0KFDJUnNzc065ZRT9OCDD2rkyJHZ3BYAgOIRQE1DW1ubZs+erbvuuks33nhjr6/PKmiYO3eudu7cqVdeeUXjx4+XJL366quqr6/XlVdeqf/5n//J5rZFxewjGeUZu+uHcotCrdS9LMmy9/zm8nZ3yHp/q0IfNaskEpZKSqSQsyd331lb4HpSIiHjeTKum6o78JK1DJ3rNaRWAzXGpOocXFmxmEybI8vpPOxkX2wrYyEw01lT4BnJeHvalrrvM2EyVx1N/2xd2b38A7S/WhAAxcnEgu5Br7W2tma8jkQiikQiPX62oaFBZ511lmbMmJG/oGHlypV6/PHH0wGDJE2YMEFLly7VzJkzs7klAADFpY/SE7W1tRmnFy5cqOuuu67bxx988EGtX79ea9euzbrJrIIGz/O6TbOUpJKSkm77UAAAgB70UdCwefNmVVRUpE/3NMqwefNmXXXVVXrsscdUWlqadZNZLe506qmn6qqrrtL777+fPvfee+/p6quv1mmnnZZ1ZwAAQO9UVFRkHD0FDevWrdO2bdt03HHHKRQKKRQKafXq1VqyZIlCoZDc1PT5T5LVSMNtt92mc889V2PHjk0Pi2zevFlHH3207r///mxuWVx6SL17suTKlmcseZ4tecm9FyzT+dVInlLHfkJT25I8W1ZJqn7AcSTbzqgrSPYhtdaB6yVzeDF9cr2BbUlyZHXWE3SpM7A6P9u1H51cT7KT/Td+1zzYuz7B209829tahlyvAzBweXn8eyGPsydOO+00vfTSSxnnLr30Uh155JH6f//v/8lxHF/3ySpoqK2t1fr16/X444/r9ddflySNHz9eM2bMyOZ2AAAUnzzOnigvL9fRRx+dcW7w4MEaMWJEt/P7k/U6DZZl6fTTT9fpp5+e7S0AAChaA3pFyCVLlujyyy9XaWmplixZst/PXnnllTl3DAAAHDhPPvlkr6/xHTTceuutmj17tkpLS3Xrrbfu83OWZfkOGp566in94Ac/0Lp16/TBBx9oxYoVGftaGGO0cOFC3XXXXWpubtb06dO1bNkyjRs3zm+3+yXLM7JcS7u8sHaakEq8hD7yhujjxGC1xksVi4Zkx2zZccmOG1mukoeX2qehyx4NkpK5+dT6DpZtywqFZEXCUmlEpiwiEw7JC9nJOobOvSi67v2QWm/BeHvdd++9GfauR2BPBgADjRuVduSprQLcGtt30LBx48Yev89Fe3u7Jk2apMsuu0xf/OIXu73//e9/X0uWLNF9992nuro6XXvttTrjjDP06quv5jRlBAAA9F5WUy5vuOEG7dq1q9v53bt364YbbvB9n1mzZunGG2/UF77whW7vGWO0ePFiffe739XnP/95TZw4UT/72c/0/vvv6+GHH86m2wAAIAdZBQ3XX3+92traup3ftWuXrr/++pw7JSVHM7Zs2ZIxI6OyslJTp07VM888s8/rotGoWltbMw4AAPobS3uKIbM6AuhzVrMnjDEZc/o7vfjiixo+fHjOnZKkLVu2SJKqqqoyzldVVaXf60lTU1POgYuXmsZi5TIVZl+MUWhbqwZvGqI/vDtBO2KDNTzcrm3Rcr3dcpA++KhS7vaIIjtshVulknajkl2enA5XdsyVlUjt++CZVI2CLYVCskKSVVIiUzFEsdEVah8d1u6RtuJDJDciyUrWUdhxyY5JTlRyYib5fdzITkiWa5JrQnQuv9A5Gyj11dhWsnSiy3lpz9oLVtcaiC7fBlHhCwDZcOMd0oY8NRbAhlW56lXQMGzYMFlWctGgww8/PCNwcF1XbW1t+sY3vtHnneyN+fPnq7GxMf26tbW127rcAACg93oVNCxevFjGGF122WW6/vrrVVlZmX4vHA5r7NixmjZtWp90rLq6WpK0detWjR49On1+69atOvbYY/d53f529wIAoN8YyLMnJKm+vl6SVFdXpxNPPLHHTav6Sl1dnaqrq7Vq1ap0kNDa2qpnn31WV1xxxQFrFwCAvBjIQUNra2t6F63Jkydr9+7d2r17d4+f7brb1v60tbVpw4Y9yaONGzfqhRde0PDhwzVmzBjNmzdPN954o8aNG5eecllTU5OxlkO+7Dcvv/d6Bt0utjI+4276h2qXt8s8OkLvjDpSb1Q4ckuSayhUJ4zsmJHTEZXT4cqJurKicVkdcVmxuBRPyCSSh1w3uXdEaqMR47VLrTsV3vqhIq+XSmWlMqVhmXBIpsSRLEvGSv0snifLNRlf5XoZ6zRYXfaSMF33aehMS1k9nOvNcwGAfibhRoPuQr/mO2gYNmyYPvjgA40aNUpDhw7tsRCys0DS725Zzz//vE455ZT0685ahPr6ei1fvlzf+c531N7erssvv1zNzc066aSTtHLlStZoAAAUvAG9jPSf/vSn9MyIJ554ok8aP/nkk/fsptgDy7J0ww039GrtBwAACsJATk989rOf7fF7AACQhYEcNHS1cuVKDRkyRCeddJIkaenSpbrrrrs0YcIELV26VMOGDevTTuaTl9rDwXTOf+36i2qUmafvuk/Dvlip9bM6ywBStQFea6usjg5FtoRVGolIJSGZkJOuf7C8LnUGnpesY3DdVN2Bl6xl6NyHwvNS3TFSZ61DLC7tbJPlOKnDllL7pXdNLRljkms+SJLx9vxcUnIviv2x8z9HGAAOKC8WdA/6taxWhPz2t7+dXmnxpZdeUmNjo84880xt3LgxY40EAADQs5xWg8yxHiJbWY00bNy4URMmTJAk/frXv9Y555yjRYsWaf369TrzzDP7tIMAAAxIBbgiZFYjDeFwOL1h1eOPP66ZM2dKkoYPH85eDwAADFBZjTScdNJJamxs1PTp0/Xcc8/poYcekiS9+eabOvjgg/u0gwdEOnefjNJcI7mWUUy24sZRzA3J8yzJVWovhuSeDEoPCX1CLYNtSZ4ty0l+bzmO1FlXYNnJ99NrHdgyiYQs102WPXSe76wr6Kw5MKn6hc57h+w9P0eKJe2pT+jsR6fONRi0Z6+Inp/NPq4HAPStAiyEzGqk4bbbblMoFNKvfvUrLVu2TJ/61KckSY8++qg+97nP9WkHAQAYiIqmpmHMmDF65JFHup2/9dZbc+4QAADon7IKGqTkrpYPP/ywXnvtNUnSUUcdpXPPPVdOalofAADYjwJMT2QVNGzYsEFnnnmm3nvvPR1xxBGSpKamJtXW1ur3v/+9Dj300D7t5IFgxz0p7mhbrFzvu4M02Itrc3yEPohVqjlapng0JDtqy44reSSMbDdVX5BeH6GHXzHbTtYWhB1ZkbCswYNkSiMypSXySpz0Hg5W55oPrpHleTKprxn7P3SpL7A80/P+D3t/DwDInhuVduSprVxTDIVS03DllVfq0EMP1ebNm7V+/XqtX79emzZtUl1dna688sq+7iMAAOgHshppWL16tdasWZPei0KSRowYoZtvvlnTp0/vs84BADBgFUt6IhKJaOfOnd3Ot7W1KRwO59wpAAAGvGIJGs4++2xdfvnluvvuu/XpT39akvTss8/qG9/4hs4999w+7eCBEn73Qx20tlYrvBP0+MFHqLw0qmgipJadZYq3RFSyw1HkY0uRj40irUYlbZ6cXQnZHQlZcVdWwk3uCdFZd2BbsuxQcj2GsjK51cPUNmaw2mocRYdKiUFGsiXLley4JTsqhTokOyo5USMnJjlxIzthZLmda0GkOmtJxraSi4fZqe/t5Peyuq+7YHX2yXS+7vLmJ/0mozwCQBFz4x3Shvy0VYhbY2dV07BkyRIddthhOvHEE1VaWqrS0lJNnz5dhx12mH784x/3dR8BAEA/0KuRBs/z9IMf/EC/+93vFIvFdN5556m+vl6WZWn8+PE67LDDDlQ/AQBAwHoVNNx000267rrrNGPGDJWVlekPf/iDKisrdc899xyo/gEAMDAN9JqGn/3sZ7r99tv19a9/XVJys6qzzjpL//3f/y3bzirTkX+WJVm23G0fauQjUY18dpjiwwcpXj5EdsRSqWXJ8ozsWEJO1JMTdWV3xGV3JKRYXFYsLsXiMglXSiQk15VxXcl1JSX3irDad8lqaVX530tVMbhMZlCpvEiJTMTJWGvB8oyshJdcq8F1ZbmpNSC6rNVg7bUWhMnYt6LL9/vbJ6Kn9SQAAN0k3GjQXejXevUv/aZNmzK2vp4xY4Ysy9L777/f5x0DAGAgG/B7TyQSCZWWlmacKykpUTwe79NOAQBQFApsILhXQYMxRpdccokikUj6XEdHh77xjW9o8ODB6XO/+c1v+q6HAACgX+hV0FBfX9/t3EUXXdRnnTmgLEtWSVhWaSS5J0QkIpWEZOIJlezYpVCrLXXWZXTuCeG6UiK1JkOqhsEkXMlzpXgiuUaDm1yvYU8zVvI+nidFozKeK7XvlhNyJMeRHDtzr4jOPSw672VM8r6pOgRjPHVjpfrZtT7iE/afMKbAwlkACIIXy19bA70Q8t577z1Q/QAAoKgUzeJOAACg+GS1jDQAAMjRQE9PFDKnfIjs0Z9SfHSFokNL5EZsGSu5VoITN7KjRk7Ukx1z5URdWbGE5KXqBFJrJqhrXYBtyfIkhULJ2gLLlkIhWSEneS7kyDipOgnLyvy17XI/y/WSY0zGJGsdjEnuHZGuk3BS7XWphbCsVJtd6hisLn3duy2p23oP+7W/NR/2Z39tZHtPAMgjq4cysgPWVgGmJ4omaAAAoF8pwJEGahoAAIAvBA0AAATB9MHRC8uWLdPEiRNVUVGhiooKTZs2TY8++miv7hFoeqKpqUm/+c1v9Prrr6usrEwnnniibrnlFh1xxBHpz3R0dOib3/ymHnzwQUWjUZ1xxhm6/fbbVVVV1au2tlw4Qc7Zu3Te2Gc1OtysD2JD9crO0Xp9+yjt3DpE4e0hRT6yFGk2irR6KmnzFGpPyNkdl9WRkBWLy0q4ybUbOtdVcGzJcWTKInIHh5UoDyte7ihRassNWzKpkMxyJds1shOSHTfJI2Fkxz3Z6b0nOmsZTLLEwZJkWzKp+gVjWzKh1Fcn9XWvMgHLKFUToeR9PCXvKUmmy/cpJl0j0fPr5Ll95M32Orf3vTPuBwAFIpHokDblp6181zQcfPDBuvnmmzVu3DgZY3Tffffp85//vP7617/qqKOO8nWPQEcaVq9erYaGBq1Zs0aPPfaY4vG4Zs6cqfb29vRnrr76av3v//6vfvnLX2r16tV6//339cUvfjHAXgMAUHjOOeccnXnmmRo3bpwOP/xw3XTTTRoyZIjWrFnj+x6BjjSsXLky4/Xy5cs1atQorVu3Tp/5zGfU0tKiu+++Ww888IBOPfVUSckFpsaPH681a9bon//5n4PoNgAAueujQsjW1taM05FIJGO7h564rqtf/vKXam9v17Rp03w32a9mT7S0tEiShg8fLklat26d4vG4ZsyYkf7MkUceqTFjxuiZZ57pMWiIRqOKRvdsbdp5TzfWIe2KKtoW1+5wQtFYXPH2mNxdUXm7Q/I6QnKjltyYUSLuyYp7yWWj3Xhq2+q4LC+1hHRnesKyJcuRcSU34SmR8JSIO0o4tjxlpieMZ+TtnZ5I7CM9odTvhb3TE0p9NcnzPaYnRHoCALKVSCT//cjL0vt9FDTU1tZmnF64cKGuu+66Hi956aWXNG3aNHV0dGjIkCFasWKFJkyY4LvJfhM0eJ6nefPmafr06Tr66KMlSVu2bFE4HNbQoUMzPltVVaUtW7b0eJ+mpiZdf/313c6/vvwGabn0Sl93HAAw4Hz00UeqrKwMuhu+bN68WRUVFenX+xtlOOKII/TCCy+opaVFv/rVr1RfX6/Vq1f7Dhz6TdDQ0NCgl19+WU8//XRO95k/f74aGxvTr5ubm3XIIYdo06ZNBfMboL9obW1VbW1tt9+Q8Ifnlz2eXfZ4drlpaWnRmDFj0iPeB1JfFUJ2zobwIxwO67DDDpMkHX/88Vq7dq1+/OMf68477/R1fb8IGubMmaNHHnlETz31lA4++OD0+erqasViMTU3N2eMNmzdulXV1dU93mtfuZzKykr+AGWpN78h0R3PL3s8u+zx7HJj23mYJ9APFnfyPC8jpf9JAp09YYzRnDlztGLFCv3pT39SXV1dxvvHH3+8SkpKtGrVqvS5N954Q5s2bepV4QYAAMVu/vz5euqpp/Tuu+/qpZde0vz58/Xkk09q9uzZvu8R6EhDQ0ODHnjgAf32t79VeXl5uk6hsrJSZWVlqqys1Fe/+lU1NjZq+PDhqqio0Ny5czVt2jRmTgAAClq+12nYtm2b/v3f/10ffPCBKisrNXHiRP3xj3/U6aef7vsegQYNy5YtkySdfPLJGefvvfdeXXLJJZKkW2+9VbZt60tf+lLG4k5+RSIRLVy48BOnn6A7nl1ueH7Z49llj2eXm7w+vzynJ+6+++4cGkuyTF7mlQAAAClZrFpZWanx/7lITqQ06/u40Q69dvt/qaWlJW/1K+w9AQAAfOkXsycAACg2ljLW0svq+nwjaAAAIAj9YMplb5GeAAAAvgzooGHp0qUaO3asSktLNXXqVD333HNBd6lfampq0gknnKDy8nKNGjVK5513nt54442Mz3R0dKihoUEjRozQkCFD9KUvfUlbt24NqMf918033yzLsjRv3rz0OZ7dvr333nu66KKLNGLECJWVlemYY47R888/n37fGKMFCxZo9OjRKisr04wZM/TWW28F2OP+w3VdXXvttaqrq1NZWZkOPfRQfe9738vYM4Hnl/TUU0/pnHPOUU1NjSzL0sMPP5zxvp/ntGPHDs2ePVsVFRUaOnSovvrVr6qtrS2nfnVOuczlyLcBGzQ89NBDamxs1MKFC7V+/XpNmjRJZ5xxhrZt2xZ01/odtijvG2vXrtWdd96piRMnZpzn2fXs448/1vTp01VSUqJHH31Ur776qn74wx9q2LBh6c98//vf15IlS3THHXfo2Wef1eDBg3XGGWeoo6MjwJ73D7fccouWLVum2267Ta+99ppuueUWff/739dPfvKT9Gd4fknt7e2aNGmSli5d2uP7fp7T7Nmz9corr+ixxx5Lr2B8+eWX59Yx0wdHng3YKZdTp07VCSecoNtuu01ScqnM2tpazZ07V9dcc03AvevfPvzwQ40aNUqrV69Ob1E+cuRIPfDAA/rXf/1XSdLrr7+u8ePH73O30WLT1tam4447TrfffrtuvPFGHXvssVq8eDHPbj+uueYa/fnPf9b//d//9fi+MUY1NTX65je/qW9961uSkvsCVFVVafny5brwwgvz2d1+5+yzz1ZVVVXG3PsvfelLKisr0/3338/z2wfLsrRixQqdd955kvz9Pnvttdc0YcIErV27VlOmTJEkrVy5Umeeeab+8Y9/qKampld96JxyedTXc59y+cqdTLnMWSwW07p16zK21LZtWzNmzNAzzzwTYM8KQ2+3KEdyddOzzjor4xlJPLv9+d3vfqcpU6bo/PPP16hRozR58mTddddd6fc3btyoLVu2ZDy7yspKTZ06teifnSSdeOKJWrVqld58801J0osvvqinn35as2bNksTz88vPc3rmmWc0dOjQdMAgSTNmzJBt23r22Wdz60ABjTJIA3T2xPbt2+W6rqqqqjLOV1VV6fXXXw+oV4Whr7YoLyYPPvig1q9fr7Vr13Z7j2e3b++8846WLVumxsZG/dd//ZfWrl2rK6+8UuFwWPX19enn09Of42J/dlJypKa1tVVHHnmkHMeR67q66aab0vsI8Pz88fOctmzZolGjRmW8HwqFNHz48JyeZb6Xke4LAzJoQPb6aovyYrF582ZdddVVeuyxx1Ramv0wYzHyPE9TpkzRokWLJEmTJ0/Wyy+/rDvuuEP19fUB967/+8UvfqGf//zneuCBB3TUUUfphRde0Lx581RTU8PzwwEzINMTBx10kBzH6Vahvr8ttbFni/Innnhin1uUd8XzTKYftm3bpuOOO06hUEihUEirV6/WkiVLFAqFVFVVxbPbh9GjR2vChAkZ58aPH69NmzZJUvr58Oe4Z9/+9rd1zTXX6MILL9Qxxxyjiy++WFdffbWampok8fz88vOcqquruxXRJxIJ7dixI7dnWYCFkAMyaAiHwzr++OMzttT2PE+rVq1iS+0esEV59k477TS99NJLeuGFF9LHlClTNHv27PT3PLueTZ8+vdvU3jfffFOHHHKIJKmurk7V1dUZz661tVXPPvts0T87Sdq1a5dsO/OvcMdx5HmeJJ6fX36e07Rp09Tc3Kx169alP/OnP/1Jnudp6tSpWbddiFMuB2x6orGxUfX19ZoyZYo+/elPa/HixWpvb9ell14adNf6HbYoz155eXm69qPT4MGDNWLEiPR5nl3Prr76ap144olatGiRLrjgAj333HP66U9/qp/+9KeSlF7v4sYbb9S4ceNUV1ena6+9VjU1NenK92J2zjnn6KabbtKYMWN01FFH6a9//at+9KMf6bLLLpPE8+uqra1NGzZsSL/euHGjXnjhBQ0fPlxjxoz5xOc0fvx4fe5zn9PXvvY13XHHHYrH45ozZ44uvPDCXs+cyFCAK0IO2KDhy1/+sj788EMtWLBAW7Zs0bHHHquVK1d2K3ZBfrYoL2Y8u56dcMIJWrFihebPn68bbrhBdXV1Wrx4cbqQT5K+853vqL29XZdffrmam5t10kknaeXKldSPSPrJT36ia6+9Vv/5n/+pbdu2qaamRl//+te1YMGC9Gd4fknPP/+8TjnllPTrxsZGSVJ9fb2WL1/u6zn9/Oc/15w5c3Taaael/zwvWbIk7z9L0AbsOg0AAPRHnes0TLxskZxwDus0xDr0t3vyu07DgB1pAACgXyvA9MSALIQEAAB9j5EGAACCUIAjDQQNAAAEoBBXhCQ9AQAAfGGkAQCAIJCeAAAAfljGyMph1YNcrs0W6QkAAOALQQMwQF1yySVFt1wwUFAKcMMq0hNAAbIsa7/vL1y4UD/+8Y/Fgq9A/1WIsycIGoAC9MEHH6S/f+ihh7RgwYKMHSOHDBmiIUOGBNE1AH4VYCEk6QmgAFVXV6ePyspKWZaVcW7IkCHd0hMnn3yy5s6dq3nz5mnYsGGqqqrSXXfdld79tby8XIcddpgeffTRjLZefvllzZo1S0OGDFFVVZUuvvhibd++Pc8/MYD+gKABKCL33XefDjroID333HOaO3eurrjiCp1//vk68cQTtX79es2cOVMXX3yxdu3aJUlqbm7WqaeeqsmTJ+v555/XypUrtXXrVl1wwQUB/yRA4etMT+Ry5BtBA1BEJk2apO9+97saN26c5s+fr9LSUh100EH62te+pnHjxmnBggX66KOP9Le//U2SdNttt2ny5MlatGiRjjzySE2ePFn33HOPnnjiCb355psB/zRAgaMQEkB/NnHixPT3juNoxIgROuaYY9LnqqqqJEnbtm2TJL344ot64okneqyPePvtt3X44Ycf4B4D6E8IGoAiUlJSkvHasqyMc52zMjzPkyS1tbXpnHPO0S233NLtXqNHjz6APQUGPmZPABhQjjvuOP3617/W2LFjFQrx1wXQp5g9AWAgaWho0I4dO/SVr3xFa9eu1dtvv60//vGPuvTSS+W6btDdA5BnBA0A9qmmpkZ//vOf5bquZs6cqWOOOUbz5s3T0KFDZdv89QHkqpBmTiT7y5JxAADkTWtrqyorK3X8+TcqVFKa9X0S8Q6t++V31dLSooqKij7s4b6RpAQAIACFWAjJ+CIAAPCFkQYAAILA7AkAAOCH5eV+9EZTU5NOOOEElZeXa9SoUTrvvPMyNrrzg6ABAIAisHr1ajU0NGjNmjV67LHHFI/HNXPmTLW3t/u+B+kJAACC0EfpidbW1ozTkUhEkUik28dXrlyZ8Xr58uUaNWqU1q1bp8985jO+mmSkAQCAAPTVLpe1tbWqrKxMH01NTb7ab2lpkSQNHz7cd58ZaQAAoIBt3rw5Y52GnkYZ9uZ5nubNm6fp06fr6KOP9t0WQQMAAEEwJnnkcr2kioqKXi/u1NDQoJdffllPP/10r64jaAAAIABBLe40Z84cPfLII3rqqad08MEH9+paggYAAIqAMUZz587VihUr9OSTT6qurq7X9yBoAAAgCHle3KmhoUEPPPCAfvvb36q8vFxbtmyRJFVWVqqsrMzXPZg9AQBAAPpq9oRfy5YtU0tLi04++WSNHj06fTz00EO+78FIAwAAQeijQkj/H8993WlGGgAAgC+MNAAAEIBC3BqboAEAgCCwyyUAABioGGkAACAApCcAAIA/nkkeuVyfZ6QnAACAL4w0AAAQhAIshCRoAAAgAJZyrGnos574R3oCAAD4wkgDAABByPMy0n2BoAEAgAAw5RIAAPhTgIWQ1DQAAABfGGkAACAAljGycqhLyOXabBE0AAAQBC915HJ9npGeAAAAvjDSAABAAEhPAAAAf5g9AQAABipGGgAACAIrQgIAAD8KcUVI0hMAAMAXRhoAAAgC6QkAAOCH5SWPXK7PN4IGAACCUIAjDdQ0AAAAXxhpAAAgCAW4uBNBAwAAASjEZaRJTwAAAF8YaQAAIAgFWAhJ0AAAQBCMpFymTbIiJAAA6K8YaQAAIACFWAhJ0AAAQBCMcqxp6LOe+EZ6AgAA+MJIAwAAQWD2BAAA8MWTZOV4fZ6RngAAIACdhZC5HL3x1FNP6ZxzzlFNTY0sy9LDDz/c6z4TNAAAUATa29s1adIkLV26NOt7kJ4AACAIea5pmDVrlmbNmpV9eyJoAAAgGH0UNLS2tmacjkQiikQiufRsn0hPAABQwGpra1VZWZk+mpqaDlhbjDQAABCEPhpp2Lx5syoqKtKnD9Qog0TQAABAMPpoymVFRUVG0HAgkZ4AAAC+MNIAAEAA8r1hVVtbmzZs2JB+vXHjRr3wwgsaPny4xowZ4+seBA0AAAQhz1Mun3/+eZ1yyinp142NjZKk+vp6LV++3Nc9CBoAACgCJ598skyO+1UQNAAAEATPSFYO/4h7bFgFAEBxYJdLAADgT45Bg/IfNDDlEgAA+MJIAwAAQSA9AQAAfPGMckoxBFAISXoCAAD4wkgDAABBMF7yyOX6PCNoAAAgCAVY00B6AgAA+MJIAwAAQSjAQkiCBgAAgkB6AgAADFSMNAAAEASjHEca+qwnvhE0AAAQhAJMTxA0AAAQBM+TlMNaC17+12mgpgEAAPjCSAMAAEEgPQEAAHwpwKCB9AQAAPCFkQYAAILAipAAAMAPYzyZHHaqzOXabJGeAAAAvjDSAABAEIzJLcXA7AkAAIqEybGmgdkTAACgv2KkAQCAIHieZOVQzBhAISRBAwAAQSjA9ARBAwAAATCeJ5PDSANTLgEAQL/FSAMAAEEgPQEAAHzxjGQVVtBAegIAAPjCSAMAAEEwRlIuUy5JTwAAUBSMZ2RySE8Y0hMAAKC/YqQBAIAgGE+5pSdYERIAgKJAegIAAAxYjDQAABCAhInmlGJIKN6HvfGHoAEAgDwKh8Oqrq7W01v+kPO9qqurFQ6H+6BX/lgmiKQIAABFrKOjQ7FYLOf7hMNhlZaW9kGP/CFoAAAAvlAICQAAfCFoAAAAvhA0AAAAXwgaAACALwQNAADAF4IGAADgC0EDAADw5f8D+fOj+P8pl9AAAAAASUVORK5CYII=",
|
|
"text/plain": [
|
|
"<Figure size 640x480 with 2 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"#8.5.6\n",
|
|
"\n",
|
|
"dt = .001\n",
|
|
"D = 1\n",
|
|
"xmax = 20\n",
|
|
"tmaxs = 100\n",
|
|
"tmax = int(tmaxs / dt)\n",
|
|
"\n",
|
|
"a = np.zeros((xmax, tmax))\n",
|
|
"b = np.zeros((xmax, tmax))\n",
|
|
"\n",
|
|
"a[:,0] = .1 # influences number of stripes\n",
|
|
"b[:,0] = .1\n",
|
|
"a[xmax - 1,0] = 1\n",
|
|
"\n",
|
|
"for t in range(0, tmax):\n",
|
|
"\n",
|
|
" if t > 0:\n",
|
|
" for x in range(0,xmax):\n",
|
|
" a[x,t] = a[x, t-1]\n",
|
|
" b[x,t] = b[x, t-1]\n",
|
|
"\n",
|
|
" for x in range(0,xmax):\n",
|
|
"\n",
|
|
" delta_a = a[x,t] ** 2 / b[x,t] - a[x,t]\n",
|
|
" delta_b = a[x,t] ** 2 - b[x,t]\n",
|
|
" \n",
|
|
" if x > 0 and x < xmax - 1:\n",
|
|
" delta_b += D*(b[x-1,t] + b[x+1,t] - 2*b[x,t])\n",
|
|
"\n",
|
|
" if x == 0:\n",
|
|
" delta_b += D*(b[x+1,t] - b[x,t])\n",
|
|
"\n",
|
|
" if x == xmax - 1:\n",
|
|
" delta_b += D*(b[x-1,t] - b[x,t])\n",
|
|
" \n",
|
|
" a[x,t] += delta_a * dt\n",
|
|
" b[x,t] += delta_b * dt\n",
|
|
"\n",
|
|
"plt.title(\"a\")\n",
|
|
"plt.imshow(a, cmap='viridis', interpolation='nearest', extent=[0, tmaxs, xmax, 0])\n",
|
|
"plt.xlabel(\"Time\")\n",
|
|
"plt.ylabel(\"Position\")\n",
|
|
"plt.colorbar()\n",
|
|
"plt.show()\n",
|
|
"\n",
|
|
"plt.title(\"b\")\n",
|
|
"plt.imshow(b, cmap='viridis', interpolation='nearest', extent=[0, tmaxs, xmax, 0])\n",
|
|
"plt.xlabel(\"Time\")\n",
|
|
"plt.ylabel(\"Position\")\n",
|
|
"plt.colorbar()\n",
|
|
"plt.show()\n"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
}
|
|
],
|
|
"metadata": {
|
|
"kernelspec": {
|
|
"display_name": "venv",
|
|
"language": "python",
|
|
"name": "python3"
|
|
},
|
|
"language_info": {
|
|
"codemirror_mode": {
|
|
"name": "ipython",
|
|
"version": 3
|
|
},
|
|
"file_extension": ".py",
|
|
"mimetype": "text/x-python",
|
|
"name": "python",
|
|
"nbconvert_exporter": "python",
|
|
"pygments_lexer": "ipython3",
|
|
"version": "3.8.10"
|
|
}
|
|
},
|
|
"nbformat": 4,
|
|
"nbformat_minor": 2
|
|
}
|