From 4aa24c3343c285ef9b6854f6dbc9c2fa6076808c Mon Sep 17 00:00:00 2001 From: felixhenneke <felixhenneke@mi.fu-berlin.de> Date: Mon, 10 May 2021 22:49:41 +0000 Subject: [PATCH] fixed some errors - wrong sign in f - wrong index shift for right boundary in B --- five_point_stencil.ipynb | 47 ++++++++++++++++------------------------ 1 file changed, 19 insertions(+), 28 deletions(-) diff --git a/five_point_stencil.ipynb b/five_point_stencil.ipynb index 5eaf5c8..dc1c4b3 100644 --- a/five_point_stencil.ipynb +++ b/five_point_stencil.ipynb @@ -40,7 +40,7 @@ "source": [ "#right-hand-side of the poisson equation\n", "def fun(x,y):\n", - " return (17*(x**2 + y**2)*np.sin(x*y) + 4*x**2*y**3+(5*x**2 + 3*y**2))" + " return -(17*(x**2 + y**2)*np.sin(x*y) + 4*x**2*y**3*(5*x**2 + 3*y**2))" ] }, { @@ -129,7 +129,9 @@ { "cell_type": "code", "execution_count": 7, - "metadata": {}, + "metadata": { + "scrolled": true + }, "outputs": [ { "data": { @@ -182,7 +184,7 @@ " # and upper boundary nodes are needed\n", " for i in range(N-1):\n", " #right boundary:\n", - " B[(i+1)*(N-1)-1][2*(N+i-1)] = pow(h,-2)\n", + " B[(i+1)*(N-1)-1][N + 2 * (i + 1)] = pow(h,-2)\n", " #upper boundary:\n", " B[-(N-1)+i][-N+i] = pow(h,-2)\n", " B = sp.csr_matrix(B)\n", @@ -318,13 +320,6 @@ "appr_u = sp.linalg.spsolve(A,RHS)" ] }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - " " - ] - }, { "cell_type": "markdown", "metadata": {}, @@ -357,13 +352,13 @@ "name": "stdout", "output_type": "stream", "text": [ - "-8.180678655519896 -8.991060763833046\n", - "-8.270439559649837 -11.480133561857317\n", - "-10.760438683213074 -12.53128157061408\n", - "-12.195353800544211 -12.961737935032156\n", - "-12.821966137009229 -13.145532611303409\n", - "-13.087791494720165 -13.228477586622512\n", - "-13.203877060573376 -13.267565784364171\n" + "-8.987207596607227 -8.991060763833046\n", + "-11.479386388270894 -11.480133561857317\n", + "-12.531189312354407 -12.53128157061408\n", + "-12.96172877159427 -12.961737935032156\n", + "-13.145531806826055 -13.145532611303409\n", + "-13.228477521210696 -13.228477586622512\n", + "-13.267565779309312 -13.267565784364171\n" ] } ], @@ -406,7 +401,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAha0lEQVR4nO3de3RU53nv8e8zuqILEtJIgEbchY2xg5FiG9tAEieu69QmNnKzcnFX0ySN18lybyerPW26mp5eVk7S1SanTW+p2+Y0p3Wck8TgW3xJ0ji1jQ3YSMbgS4zAXCQwIAkE6H55zh8zkgWW0Ahpa8+Mfp+19hpmz8zeP5kxj979vvt9zd0REREBiIQdQEREUoeKgoiIjFBREBGRESoKIiIyQkVBRERGqCiIiMiI7LADTEU0GvWlS5eGHUNEJK3s2rWr1d0rxnotrYvC0qVLeemll8KOISKSVszs0Hiv6fKRiIiMUFEQEZERaVkUzGyTmd3X0dERdhQRkYySlkXB3R9193tKSkrCjiIiklHSsiiIiEgwVBRERGREWg9JvVRNJ85ypL2beYW5lBXkUlaUS2FuFmYWdjQRkVDNyqLw8MtH+dufNp23Lzc7QllBLvMKcykvHPWYKBplBbmUFb6zzSvIITtLDS0RySyzsih86salfHBVJe2dfe9sXX20n+vjVFcfbZ19NJ/qor2zjzM9A+MeZ25+NuVFecwryKGsMI+ywnEe1RoRkTQxK4tCtCiPaFFeUu/tHxziVFffeQXkVGe8cJzq7KO9q5/2zl5aTnezp+U07Z199A+OvZpdblYk3soY1Ropu6CgzCvMoTzxOK8glxy1RkRkBs3KojAZOVkRKovzqSzOT+r97s653gFOdfbT1tmbKCjxwnHhY8vpbtrO9U66NTL60lZ5US7z5+azpLyQojz9dYrI1OhfkWlmZhTn51Ccn8Pi8oKkPtM/OMTprv6LXs461dnH0dPd7G3poL2zj77BoXcdJ1qUx9LyApaUF8Yfo4Ujz0vm5Ez3jyoiGUhFIQXkZEWoKM6joji5S1ruTmff4MhlrGOnuznY1sWhtk4OtnWyramVBxt6zvvMvIKcd4pFeSFLo/HHZeWFlBbkqK9DRAAVhbRkZhTlZVOUl82isgLWLip913u6+wY53N7FwbbORLGIF40XD57i4d1H8VHdHnPzs1lSXsiS8gKWDj9G448VRXkqGCKzSFoWBTPbBGyqqakJO0rKmpObxeULirl8QfG7XusdGORIe/d5xeJgWxd7Wjp4Yu/bDA69UzEKcrPOb2GMamnML84nElHBEMkk5j72SJl0cM0117jWU5he/YNDtJzqTrQwus57PNLedd7IqrzsCEsuLBaJlkZV6RyyVDBEUpKZ7XL3a8Z6LS1bChKcnKwIS6OFLI0Wvuu1wSHn6OnuUcXinZbGM2+epHfgnc7v3KwI1WVz3rkcNeoxNm+OhtqKpCgVBUlaVsRYVFbAorICNqyMnvfa0JBz/GwPB1u73nVZavuBNrr6Bs87TvW8OWNelloWLVQLQyREKgoyLSIRY2HJHBaWzOGGFeXnvebunDzXG29htMZHSA0XjcZDpzjb+859GiVzcrhheTnrV0ZZv6KcZdFCdXSLzCAVBQmcmY3cAHjt0rLzXnN32jv7ONjWxYGT59j5Vjvbmlp58tW3AVhYks/6mijra8pZvyJK5dzkbiIUkUujjmZJOe7OwbYunmtq5fmmVl440Mbprn4AVlYWJYpElHXLy5ibr5vyRCbrYh3NKgqS8gaHnNeOnmHb/la2NbXy4sF2evqHiBisqS5lQ02UG2vKqVs8j/ycrLDjiqQ8FQXJKL0DgzQcOs3ziSKxu7mDwSEnLzvCdcvKuHFF/HLTlVUl6rQWGYOKgmS0sz397DjQPtKSePP4OWBUp3VNOetrouq0FknQfQqS0Yrzc7h59XxuXj0fgBNne3hhfxvP7Wvl+f1t53Va37giyoaV6rQWGY9aCpLR3J1Dw53W++NFYrjTuqayKN4fsaKc61eUq9NaZg1dPhJJGBpyXjt2hm1NrWzb38bOt9rO67QeHvpat0Sd1pK5Mq4ojJoQ73P79u0LO46ksd6BQRoPn+b5plaeu6DT+tqlZdxYU86Gmqg6rSWjZFxRGKaWgky3sz39iRvo2tjW1MrPj58F4tOL37CiPDH8NcpydVpLGlNHs0iSivNz+NAV8/nQFfFO65Nne0eGvm5rauOpV48DsGDuqDuta6LMV6e1ZAi1FESS5O4cbh++07qN5/e3cmpUp/X6FeWJO63LtfyppDRdPhIJwHCndbwl0cbOt9rp7h8kYvCe6lI2JqbjeO+SeeRma6pwSR0qCiIzoG9giMbDp9i2P94f8fKR0wwOOXNysli3vIwNNVE2rqzgsvlF6o+QUKkoiITgTE8/2/e38VxiZNOBk50AVBTnsaEmGt9Wqj9CZp6KgkgKaDndzbZ9rTzbFO+4bu/sA+Cy+UVsqKlg48oo1y0rozBP4z8kWCoKIilm9E10zzW1svOtdnoHhsjJMmoXz2NjohWxprpU90fItFNREElxPf2DvHTwFM82nWRbUyt7W84A8fsj4vM1xS83LSkvUH+ETJnuUxBJcfk5WfF/+BNrX7d39sVbEfviLYnhSf2q581h48ooG2oquHFFOfMKc8OMLRlILQWRFOfuvNXaybamVp7d18oL+9s42zuAGbwnVsL6migba6K8d+k88rI1X5NMTJePRDLIwOAQu5s7Eq2IkzQePs3AkJOfE+G6ZeUj/RGrFhTrUpOMSUVBJIOd6x1gx4E2nk1camo6EV9kKFqUy/rE0NeNKytYUKKhrxKnPgWRDFaUl33efE3HOrrZ1tTGc/tO8lxTGw+/fBR4Z/2IjSvjU3EUaeirjEEtBZEM5u688fZZnkvcHzG8fkR2xKhdXMqGmgo2rIxydXUJ2VmaimO20OUjEQHi60fsOnRqZFTTnpYO3KE4P5sblpePDH3VetaZLeOKghbZEZkepzr7eH5kKo6THGnvBiBWOocNNVHWr4yPbNLQ18yScUVhmFoKItPrUFtnvMN6X3xN6zM9A+RkGbdetZC71y1m3bIytSAygIqCiEza4JCzp6WDR14+yg92HeFMzwArKgq5e90S7qqrpqRAa0akKxUFEZmS7r5BfrjnGPfvOETj4dPkZUfYdHUVd69bzNpFpWo9pBkVBRGZNq8e7eA7Ow7zUGMLnX2DrF44l7uvX8wda2Ma5pomplQUzKwU+FVgKaPua3D335q+iJdGRUEkPOd6B3j45Rb+Y/thXj92hsLcLO6sjXH3uiWsrpobdjy5iKkWheeB7cAeYGh4v7t/ezpDXgoVBZHwuTsvHznN/TsO8+juo/QODFG7uJS71y3h9jULyc/RfEypZqpFocHd6wJJNkUqCiKppaOrnwcbmrl/xyH2n+ykZE4Od9VV88l1i6mpLAo7niRMtSj8d+Ac8BjQO7zf3dunM+SlUFEQSU3uzo632rl/x2Ge3HuM/kHn+uVl3L1uCb945QJys3X3dJimWhTuBb4MnAaG3+zuvnw6Q14KFQWR1Nd6rpfvv9TMd3Ye4kh7N9GiXD56zSI+ed1iFpUVhB1vVppqUdgPrHP31iDCTYWKgkj6GBpynm1q5f7th/jJ68dx4H0rK7h73WI+uKpScy/NoKnOkvoq0DW9kURktolEjPdfVsH7L6vgWEc3/+/FI3x35xHu+fddLJibz8evW8THr12sKb5DlkxLYStwJfA05/cpaEiqiEzJwOAQP33jBPfvOMwz+04SMeNDqyq5+/olbKyJEonoprggTLWl8FBiExGZVtlZEW65cgG3XLmAw21dPPDiYb734hF+9NpxFpcV8InrFvPRa6qJFuWFHXXWuGhLwcyygKfc/eaZi5Q8tRREMk/vwCA/evU49+84xPYD7ZqQLwCX3FJw90Ez6zKzEnfvCCaeiMg78rKz2HR1FZuurqLpxFm+s+MIP9h1hEd3H9WEfDMgmT6F7wHXAz8GOof3q09BRGaKJuSbXlMdkvqpsfZrmgsRCYMm5Ju6Kc+Sama5wGWJpz939/5pzHfJVBREZi9NyHfpptpS+ADwbeAgYMAi4FPu/sy0ppwELccpIsM0Id/kTbUo7AI+6e4/Tzy/DHjA3d877UknSS0FERlNE/IlZ6pF4RV3XzPRvjCoKIjIWMabkO+PblvNVbGSsOOFbqpF4VvEJ8L798Suu4Fsd//0tKa8BCoKIjKR4Qn5/uXZA1SXFfDwvevDjhS6qd7R/HngXuC3iPcpPAP8w/TFExEJTrQoj89/YAXZEePLj7/O/pPnWFGhS0njmXBaQnfvdfevu3u9u2929//t7r0TfU5EJJXcsbaKiMHWhpawo6S0CYuCma03sx+b2ZtmdmB4m4lwIiLTpXJuPhtWVrC1sYWhoYmH4s9WyUxg/q/A14ENwLWjNhGRtFJfG6PldDc7D4a+cGTKSqYodLj7E+5+wt3bhrfAk4mITLNbrpxPYW6WLiFdRDJF4Wkz+0szu8HM6oa3wJOJiEyzgtxsbr1qIY/vOUZP/2DYcVJSMqOP1iUeRw9fcuCD0x9HRCRYd9XFeLChmR+/dpxNV1eFHSflTFgU3P2mmQgiIjITrl9ezsKSfLY0NKsojGFSK2Wb2WNBBRERmQmRiHFnbYxn9rVy8qxG119oUkUBiAWSQkRkBtXXxhgcch7ZfTTsKClnskWhMZAUIiIzaOX8Yt4TK2FrY3PYUVLOpIqCu38mqCAiIjNpc22MvS1nePP42bCjpJRk7mi+3cwazazdzM6Y2VkzOzMT4UREgvKRtVVkRYwtumfhPMm0FP4a+BRQ7u5z3b3Y3bWskYiktWhRHu+/rIKHGlsY1LQXI5IpCkeAvZ7Mup0iImmkvi7G22d62H5AkzQMS+bmtf8BPG5m/wWMjN9y968HlkpEZAbcfMV8ivOyebChmfU10bDjpIRkWgpfBrqAfKB41CYiktbyc7K4bc1Cntz7Nl19A2HHSQnJtBTK3P2WwJOIiIRgc22M7754hKdefZvNtdVhxwldMi2Fn5iZioKIZKRrl5ZRPW+ORiElJFMU7gWeNLNuDUkVkUwTiRiba2Nsa2rl+JmesOOELpnlOIvdPeLuczQkVUQy0ebaGEMOD7+s1sKEfQpm9r6x9rv7M9MfR0Rk5i2vKGLtolK2NLRwz/tWhB0nVMl0NP/eqD/nA9cBuwhxPQUz2wRsqqmpCSuCiGSYu+pifOnhV3nt6BlWV83eiyHJXD7aNGr7BeAq4Hjw0S6a6VF3v6ekpCTMGCKSQW5fU0VOlrGlYXZPkjfZWVIBmokXBhGRjDGvMJebLq/k4d1HGRgcCjtOaJLpU/hb4stvQryIrAV2B5hJRCQU9XUxfvTacZ5rauUDl1eGHScUyfQpvDTqzwPAA+6+LaA8IiKhuWlVJSVzctja2KKiMB53//ZMBBERCVtedha3r1nIgw3NnOsdoCgvmd+bM4vWUxARGaW+rpqe/iGe2HMs7Cih0HoKIiKj1C0uZWl5wayd9kLrKYiIjGJmbK6tZvtbbbSc7g47zoxLpigMr6fwRTP7wvAWdDARkbBsro3hDg81zr7WgtZTEBG5wOLyAq5dOo+tjS3MtoskWk9BRGQMm2ur+cOte9jT0sGa6tKw48wYracgIjKG296zkNzsyKzrcNZ6CiIiYygpyOHmKyp5dPdR+mfRtBdaT0FEZBz1tdW0dfbxzJsnw44yY8btUzCzVe7+hpnVjfW6uzcEF0tEJHzvv7yCssJctjS08KEr5ocdZ0ZcrKP5C8A9wNfGeM0JcT0FEZGZkJMV4SNXV/GdnYfp6O6nZE5O2JECN+7lI3e/J/F40xibCoKIzAqba2P0DQzx+CyZ9iKZuY/+3MyyRj2fa2b/J9hYIiKpYU11CSsqCtk6S0YhJTP6KBvYaWZrEkNTXyS+HKeISMYzM+rrqtl5sJ0j7V1hxwlcMqOPvgj8PrAD+DfgNnf/u4BziYikjDtrYwBsnQXTXiRz+eh9wN8Afwb8DPg7M6sKOJeISMqIlc7h+uVlbGlozvhpL5K5fPRXwEfd/Svu/kngPuCnwcYSEUkt9XXVHGzrovHI6bCjBCqZonCDu782/MTdtwDrg4skIpJ6PnzVAvKyI2xpaA47SqCSWY5z0MxuA64kPlPqsD8LLJWISIopzs/hF69cwGOvHONLt68mLztr4g+loWT6FL4JfAz4TcCAjwJLAs4lIpJyNtfFON3Vz9NvZO60F8lcPrrR3X8VOOXufwrcACwKNpaISOrZWBMlWpTH1sbMvYSUTFEYXo+uKzHqqB9YFlwkEZHUlJ0V4Y61Vfz0jROc6uwLO04gkikKj5lZKfCXQANwEPhugJlERFJWfV2M/kHnsQyd9iKZm9f+3N1Pu/uDxPsSVrn7l4KPJiKSelYvnMvl84szdhTShKOPEvMe3QYsHX6/meHuXw82mohI6olPexHjK0+8wVutnSyLFoYdaVolc/noUeDXgHKgeNQmIjIr3bE2hhlszcDWwoQtBaDa3dcEnkREJE0sKMlnQ02ULY0t/M7NlxGJWNiRpk0yLYUnErOjiohIwubaGM2nunnp0Kmwo0yrZIrCdmCrmXWb2RkzO2tmZ4IOJiKSyn7xygUU5GZl3D0LyRSFrxG/Ya3A3ee6e7G7zw04l4hISivMy+bWxLQXPf2DYceZNskUhX3AXk+h+WLNbJOZ3dfR0RF2FBGZxTbXxTjbM8B/vn4i7CjTJpmicAz4mZl90cy+MLwFHexi3P1Rd7+npKQkzBgiMsvduCLK/Ll5GXXPQjJF4S3gP4FcNCRVRGREVsS4c22M/3rzJK3nesOOMy2SmTr7T2ciiIhIOqqvq+afnjnAo7uP8un16T8tXDItBRERGcflC4pZvXBuxqzfrKIgIjJF9XUxXmnuoOnE2bCjTJmKgojIFH1kbRURgy0N6d9aSGZCvArgc4yaEA/A3T8TXCwRkfRRWZzP+y6r4KHGFn73lsvTetqLZFoKDwMlwE+AH47aREQkYXNtjKMdPWx/qy3sKFOSzIR4Be7++4EnERFJY7esXkBRXjZbGlq4cUU07DiXLNmV134p8CQiImlsTm4WH75qAU/sOUZ3X/pOe5FMUfht4oVBE+KJiFxEfV01nX2D/Oi1t8OOcsmSWY6z2N0j7j5HE+KJiIxv3bIyYqVz0noU0rh9Cma2yt3fMLO6sV5394bgYomIpJ9IxLiztop//Nl+TpzpoXJuftiRJu1iHc1fAO4hPnX2hRz4YCCJRETS2Obaav7+6f08svsov75xedhxJm3couDu9yQeb5q5OCIi6a2msoirq0t4sKElLYuC7mgWEZlmm2tjvH7sDK8fS78xOSoKIiLTbNPVVWRHLC0nybtoUbC4RTMVRkQkE5QX5fGBy+PTXgwOpcyilUm5aFFILMH50MxEERHJHPV11Zw428u2ptawo0xKMpePtpvZtYEnERHJIB9cVUlxfnbaXUJKpijcBLxgZvvN7BUz22NmrwQdTEQkneXnZHH7miqe3Ps2nb0DYcdJWjIT4n048BQiIhmovi7GAzsP8+Tet7nrvdVhx0lKMtNcHAJKgU2JrTSxT0RELuKaJfNYVDaHLY3NYUdJ2oRFwcx+G7gfqExs/2Fmvxl0MBGRdGdmbK6t5vn9bRzr6A47TlKS6VP4LLDO3f/Y3f8YuJ74SmwiIjKB+toY7vBQ49GwoyQlmaJgwOjJwQcT+0REZAJLo4XULS5lS0Mz8VH+qS2ZovAtYIeZ/YmZ/QmwHfjXQFOJiGSQ+rpq9p04x6tHU3/ai4nuaI4AO4BPA+3AKeDT7v7XwUcTEckMt69ZSG5WJC3WWZjojuYh4Gvu3uDu33D3v3H3xhnKJiKSEUoLcvngqkoe2d3CwOBQ2HEuKpnLRz8ys7vMTP0IIiKXaHNdjNZzfTy7L7WnvUimKHwB+D7QqzWaRUQuzU2XV1JakMODDal9z0IyfQq3JtZoztUazSIilyY3O8KmNVX8+LXjnOnpDzvOuJLpU/irGcoiIpLR6uti9A4M8cSeY2FHGZf6FEREZsjaRaUsixam9CikyfQp9KlPQUTk0pkZ9bUxdrzVzpH2rrDjjCmZCfGKE30KOepTEBGZmjtrYwA8/HJqthaSmRDPzOxXzOxLieeLzOy64KOJiGSeRWUFXLesjC0NLSk57UUyl4/+AbgB+GTi+Tng7wNLJCKS4eprYxxo7WR3c0fYUd4lmaKwzt3vBXoA3P0UkBtoKhGRDPZLaxaSmx1hSwres5BMUeg3syzAAcysAkjt+7RFRFLY3PwcfmH1fB7dfZS+gdT65zSZovANYCtQaWZfBp4D/legqUREMtxddTFOdfXzs5+fCDvKeSZco9nd7zezXcCHiK+jcKe7vx54MhGRDLZxZQXlhblsbWzhlisXhB1nxIRFAcDd3wDeCDiLiMiskZMV4SNrq7h/+2E6uvopKcgJOxKQ3OUjEREJQH1tNX2DQzy2J3WW6lRREBEJyVWxuaysLEqpaS9UFEREQmJmbK6LsevQKQ61dYYdB1BREBEJ1Z1rY5jB1sbUaC2oKIiIhKiqdA43LC9na2NqTHuhoiAiErL6umoOtXXRcPhU2FFUFEREwnbrVQvIz4nwYAp0OKsoiIiErCgvm1uvXMBju4/SOzAYahYVBRGRFLC5rpozPQP89PVwp71QURARSQHrV5RTWZzHlpBHIakoiIikgOysCHesreLpN07Q3tkXWg4VBRGRFFFfV83AkPPYK+FNe6GiICKSIq5YOJdVC4pDHYWkoiAikkLuqqtm95HT7D95LpTzp0xRMLNCM/u2mf2zmd0ddh4RkTDcsbaKiMHWkFoLgRYFM/uWmZ0ws70X7L/VzH5uZk1m9geJ3fXAD9z9c8BHgswlIpKqKufms2FlBVsbWxgamvlpL4JuKfwbcOvoHYn1nv8e+DCwGviEma0GqoEjibeFe/eGiEiI6mtjtJzuZufB9hk/d6BFwd2fAS78qa4Dmtz9gLv3Ad8F7gCaiReGi+Yys3vM7CUze+nkyZNBxBYRCdUtV86nMDcrlEtIYfQpxHinRQDxYhADtgB3mdk/Ao+O92F3v8/dr3H3ayoqKoJNKiISgoLcbG69aiGP7zlGT//MXjgJoyjYGPvc3Tvd/dPu/nl3v3/GU4mIpJC76mKc7R3gx68dn9HzhlEUmoFFo55XA6mzQKmISAq4fnk5C0vy2dLQPKPnDaMovAisNLNlZpYLfBx4JIQcIiIpKxIx7qyN8cy+Vk6e7Z258wZ5cDN7AHgBuNzMms3ss+4+APwG8BTwOvA9d381yBwiIumovjbG4JDzyO6Zu5iSHeTB3f0T4+x/HHg8yHOLiKS7lfOLeU+shK2NzXx2w7IZOWfK3NEsIiLvtrk2xt6WM7x5/OyMnE9FQUQkhX1kbRVZEWPLDN2zoKIgIpLCokV5vP+yCh5qbGFwBqa9SMuiYGabzOy+jo6OsKOIiASuvi7G22d62H6gLfBzpWVRcPdH3f2ekpKSsKOIiATu5ivmU5yXzYMzcM9CWhYFEZHZJD8ni9vWLOTJvW/T1TcQ6LlUFERE0sDm2hhdfYM89erbgZ5HRUFEJA1cu7SM6nlzAh+FpKIgIpIGIhFjc22MbU2tHD/TE9x5AjuyiIhMq821MYYcHn45uNaCioKISJpYXlHE2kWlgV5CUlEQEUkjd9XFeOPts7x29Ewgx1dREBFJI7evqSJalMehts5Ajh/oLKlBMbNNwKaampqwo4iIzKh5hbns+MMPkRUZaxHLqUvLloLuaBaR2SyoggBpWhRERCQYKgoiIjJCRUFEREaoKIiIyAgVBRERGaGiICIiI1QURERkhLkHv+ZnUMzsJHBoiocpAaayruelfn4yn0v2vRO9b6LXo0BrkplS3VT/XlPlnNNxzEs5Rhjfz2Tekynf0TC+n6PPu8TdK8Z8h7vP6g24L4zPT+Zzyb53ovcl8fpLYf99pMrfa6qcczqOeSnHCOP7mcx7MuU7Gsb3M9nz6vIRPBrS5yfzuWTfO9H7pvqzppMwftYgzjkdx7yUY4Tx/ZzsedNZWD/nhOdN68tHMr3M7CV3vybsHCLj0Xc0eGopyGj3hR1AZAL6jgZMLQURERmhloKIiIxQURARkREqCiIiMkJFQZJiZnea2T+b2cNmdkvYeURGM7PlZvavZvaDsLOkOxWFWcDMvmVmJ8xs7wX7bzWzn5tZk5n9wcWO4e4PufvngF8DPhZgXJllpun7ecDdPxts0tlBo49mATN7H3AO+L/uflViXxbwJvALQDPwIvAJIAv4ygWH+Iy7n0h87mvA/e7eMEPxJcNN8/fzB+7+yzOVPRNlhx1Agufuz5jZ0gt2Xwc0ufsBADP7LnCHu38FuP3CY5iZAV8FnlBBkOk0Hd9PmT66fDR7xYAjo543J/aN5zeBm4FfNrP/FmQwESb5/TSzcjP7JlBrZl8MOlwmU0th9rIx9o17LdHdvwF8I7g4IueZ7PezDdAvK9NALYXZqxlYNOp5NXA0pCwiF9L3MyQqCrPXi8BKM1tmZrnAx4FHQs4kMkzfz5CoKMwCZvYA8AJwuZk1m9ln3X0A+A3gKeB14Hvu/mqYOWV20vcztWhIqoiIjFBLQURERqgoiIjICBUFEREZoaIgIiIjVBRERGSEioKIiIxQURABzOx3zKwgwONXTcdc/2Z2bjryiIxH9ymIAGZ2ELjG3VvDznIxZnbO3YvCziGZSy0FmVXMrNDMfmhmu81sr5l9zMx+C6gCnjazpxPvu8XMXjCzBjP7vpkVJfYfNLO/MLOdia1mjHO838xeTmyNZlZsZkuHF5Exs38Z9fpJM/ufif2/Z2YvmtkrZvanF/kZvpzIv93M5gfx30lmLxUFmW1uBY66+9WJBV2eTMwAexS4yd1vMrMo8EfAze5eB7wEfGHUMc64+3XA3wF/PcY5fhe4193XAhuB7tEvuvuvJ167A2gD/i2xxOlK4usIrAXem1h85kKFwHZ3vxp4Bvjc5P8TiIxPRUFmmz3AzYnf9je6e8cY77keWA1sM7OXgU8BS0a9/sCoxxvG+Pw24OuJFkhpYh6f85hZPvB94Dfc/RBwS2JrBBqAVcSLxIX6gMcSf94FLB3/RxWZPK2nILOKu79pZu8Ffgn4ipn9yN3/7IK3GfBjd//EeIcZ58/D5/iqmf0wcY7tZnYz0HPB274JbHH3n4w651fc/Z8m+BH6/Z2OwEH0/7BMM7UUZFYxsyqgy93/A/groC7x0lmgOPHn7cD64f4CMysws8tGHeZjox5fGOMcK9x9j7v/BfFLT6sueP1eoNjdvzpq91PAZ0b1XcTMrHIKP6rIJdFvGTLbvAf4SzMbAvqBzyf23wc8YWbHEv0KvwY8YGZ5idf/iPhC8gB5ZraD+C9VY7UmfsfMbiL+m/xrwBPAwlGv/y7Qn7g0BfBNd/+mmV0BvBBfDptzwK8AJ6b6A4tMhoakikxCugxdFblUunwkIiIj1FIQEZERaimIiMgIFQURERmhoiAiIiNUFEREZISKgoiIjFBREBGREf8f6yXLYPU7gx4AAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEKCAYAAAAMzhLIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deXxU9fX/8ddh34OKqKwBAwIiKEbUuhStCy6otVjr0lZRqf5qrbV1Qy2KbV2odtOq1L22+lUBJYhaFxS/VpRFgYRNxIWAEEANa0JIzu+PDN+GNJncJHPnzkzez8djHsncuXPvGRk5nPv5fO4xd0dERKQ2zaIOQEREUpsShYiIxKVEISIicSlRiIhIXEoUIiISlxKFiIjE1SLqAMLQpUsXz87OjjoMEZG0Mm/evA3uvnf17RmZKLKzs5k7d27UYYiIpBUz+7ym7br0JCIicWVUojCzUWY2qbi4OOpQREQyRkYlCnfPc/exWVlZUYciIpIxMipRiIhI4ilRiIhIXEoUIiISlxKFiEgGKCuvYPK8QsorEt86IqPWUZjZKGBUTk5O1KGIiCTNvz/ZwK3TCli+bgtZbVtywqB9Enr8jKooNOtJRJqSNd9s56f/nM/5f3ufbTvKmfTDQ/nOwK4JP09GVRQiIk1B6c5yHn7nU+57cwUV7lx9Qj8u//b+tGnZPJTzKVGIiKSRmUuLuC2vgM82buPkA/fh5tMG0XPPdqGeU4lCRCQNfLFxGxOmF/D6kiL6dmnPE2OG8+3+/3X/vlAoUYiIpLDtO8p54K0VPDhrJS2aGTecMoAxR/WhVYvkDTErUYiIpCB359WCtdw+fQmrv9nOGUO7Me7Ugeyb1SbpsWRUotD0WBHJBCuKtnBbXgHvfLyBAft25JmxR3BE370ii8fcE784I2q5ubmufhQikm62lO7kz298zKP/+yltWzXnlyf258IjetOieXIuM5nZPHfPrb49oyoKEZF05O68+NEafjdjCUWbS/l+bg+uGzmALh1aRx0aoEQhIhKpxWs2MX5aPnM++5ohPbJ46IeHckivPaIOazdKFCIiESjeVsY9ry3jqdmfk9W2JXeefRDfz+1Js2YWdWj/RYlCRCSJKiqcZ+eu4u5Xl/HNth1ceERvrjmxP53btYo6tFopUYiIJMlHq75h/Iv5LCgs5rDsPbjtjMMZ1K1T1GHVSYlCRCRkG7aUMvGVZfzP3FXs3bE1fzh3KGcd3B2z1LvMVBMlChGRkOwsr+Cp2Z9zz2vL2b6jnLHH9uVnx+fQsU3LqEOrl4xKFFpwJyKp4v2VGxk/rYClazdzdE4Xbj1jEDldO0YdVoNkVKJw9zwgLzc397KoYxGRpmndphJ+N2MJL360hu6d2/LABcMYOXjftLnMVJOMShQiIlHZsbOCx979lD+/8TFlFc5Vx+dwxYgc2rYKp0dEMilRiIg00qzl67k1r4CV67dywsCu3HL6IHrv1T7qsBJGiUJEpIFWfbWN37y0mFcL1pG9Vzseu+gwjhuQ+FakUVOiEBGpp5Kych56eyV/fWsFzcy49uQDuPSYPrRukf6XmWqiRCEiEpC78/qSIiZML2DVV9s5bch+3HTqQLp1bht1aKFSohARCWDl+i3clreYt5evp1/XDvzz0sP5Vk6XqMNKCiUKEZE4tpbu5L6ZK3j4nZW0adGcW04fxI+O7E3LJPWISAVKFCIiNXB3pi/8kt++tIS1m0r43rAeXH/KAXTtmPxWpFFTohARqWbZ2s2Mn5bP7JVfcWC3Ttx/wSEc2nvPqMOKjBKFiEhM8fYy/vj6cp5873M6tmnBb84azHnDe9E8BXtEJJMShYg0eRUVzuT5hdz1ylI2bt3B+cN78auTDmCP9qnbIyKZlChEpElbVFjMr6fl8+EX33BIr848fvFwBnfPijqslJLyicLMBgI/B7oAb7j7AxGHJCIZYNuOndzzr+U89u6n7Nm+Fb8/ZyhnH9I9JVuRRi3URGFmjwKnA0XuPrjK9pHAn4DmwMPufmdtx3D3JcDlZtYMeBJQohCRRpm1fD3jpi6i8OvtXHB4L64bOYCstunVIyKZwq4oHgfuo/IveADMrDlwP3AiUAjMMbNpVCaNO6q9f4y7F5nZGcAVwN9DjldEMtjXW3dw+0uLmTJ/NX27tOfZnxzJ8D5NdzZTUKEmCnefZWbZ1TYPB1a4+0oAM3sGONPd76Cy+qjpONOAaWb2EvDPmvYxs7HAWIBevXolJH4RyQzuzrQFa5iQt5ji7WVceVwOVx6fQ5uWmXlvpkSLYoyiO7CqyvNC4PDadjazEcDZQGtgRm37ufskYBJAbm6uJyJQEUl/q7/Zzi0v5PPm0iKG9sjiqUsPZ+B+naIOK62k/GC2u78FvBVxGCKSZioqnL/P/py7X1lKhcPNpw3k4qP6NPk1EQ0RRaJYDfSs8rxHbFujqWe2iAB8vG4z109eyPwvvuGYfl343XcPouee7aIOK21FkSjmAP3MrA+VCeIHwPmJOLB6Zos0baU7y3ngrU+4f+YK2rduwb3fH8p3D+me1v2qU0HY02OfBkYAXcysEBjv7o+Y2ZXAq1TOdHrU3QvCjENEMt+8z7/mhskL+bhoC2cM7cavRw2iS4fWUYeVEcKe9XReLdtnEGdguqF06Umk6dlSupPfv7qMJ977jP06teHRi3I5fsA+UYeVUVJ+MLs+dOlJpGmZuayIm6fms6Z4Oz86ojfXjhxAh9YZ9ddaStB/URFJOxu3lDJh+mJe/GgNOV078PzlRzbp24CHrc5EYWadgR8B2VX3d/erwgurYXTpSSSzuTsvfLSaCXmL2VK6k6tP6McVI/andQstnAtTkIpiBjAbWARUhBtO4+jSk0jmWvXVNm56IZ9Zy9dzSK/O3PW9IfTfp2PUYTUJQRJFG3e/JvRIRERqUF7hPP7vz/j9q8toZnDbGQdy4RG9tXAuiYIkir+b2WXAdKB010Z3/yq0qEREgKVrN3H95EUsWPUNxx2wN7/57kF079w26rCanCCJYgcwEbgJ2HUPJQf6hhVUQ2mMQiQzlJSVc//MFTzw1id0atuSP/3gYM4Y2k0L5yJi7vHvn2dmK4Hh7r4hOSE1Xm5urs+dOzfqMESkAeZ89hU3TF7IJ+u3cvaw7tx82iD2VEvSpDCzee6eW317kIpiBbAt8SGJiPzH5pIy7nplKU/N/oLundvyxJjhfLv/3lGHJQRLFFuBj8xsJruPUaTc9FgRSU+vLV7HLS/kU7S5hEuO7sM1J/anvRbOpYwgfxIvxB4pT2MUIull/eZSbs0r4KWFXzJg3448+MNDObhn56jDkmrijlHE2pa+7u7HJS+kxtMYhUhqc3eem1fIb19awvYd5Vz1nRzGHrs/rVo0izq0Jq1BYxTuXm5mFWaW5e7F4YUnIk3FFxu3cePUhby7YiPDs/fkd2cfRE7XDlGHJXEEufS0BVhkZq9ROV4BaIxCROpnZ3kFj737Gfe8towWzZrxm7MGc/7wXjTTwrmUFyRRTIk9REQapGBNMTdMXsSi1cWcMHAffnPWYPbNahN1WBJQnYnC3Z8ws1ZA/9imZe5eFm5YIpIJSsrK+dMbHzNp1kr2aNeKv14wjFMG76uFc2kmyN1jRwBPAJ8BBvQ0sx+7+6xwQ6s/zXoSSR2zV27kximL+HTDVr6f24Nxpw6kczstnEtHQVZmzwPOd/dlsef9gafd/dAkxNcgmvUkEp3i7WXc+fISnv5gFb32bMcdZx/EUTldog5LAmjMyuyWu5IEgLsvN7OWCY1ORDLCK/lr+fWL+WzYUspPju3L1Sf0p20r9YpId0ESxVwzexh4Kvb8AkD/XBeR/7NuUwnjXyzglYK1DNqvE49edBiDu2dFHZYkSJBEcQXwU2DXdNh3gL+GFpGIpA1355k5q/jdjCXs2FnB9SMHcOkxfWjZXAvnMkmQWU+lwL2xh4gIAJ9u2MqNUxYye+VXHNF3T+44ewh9urSPOiwJQZBZT0cBtwK92b1ndsr1oxCR8O0sr+Bv73zKH15fTusWzbjz7IM497CemvKawYJcenoE+AUwDygPN5zG0fRYkXAtX7eZa59bwILCYkYeuC8TzjyQrp20cC7TBUkUxe7+cuiRJIC75wF5ubm5l0Udi0gmKSuv4KG3P+HPb6ygQ5sW3H/+ME4bsl/UYUmSBEkUM81sIpW38ajaj2J+aFGJSMpY8uUmrn1+AfmrN3HakP2YcMaB7NWhddRhSRIFSRSHx35WXYThwPGJD0dEUkVZeQUPvPUJf3nzY7LatuSBC4ZxykGqIpqiILOe0qoXhYg0XsGaYq59biGLv9zEmQd3Y/yoA9W3ugmrV69BM5vu7qeHFYyIRGvHzgrum7mCv85cQed2rXjoh4dy8oH7Rh2WRKy+TWm7hxKFiEQuf3Uxv3puAUvXbua7h3Rn/KhBuomfAPVPFB+GEoWIRKZ0Zzl/eWMFD7z9CXu1b8XDP8rlhEH7RB2WpJB6JQp3HxNWICKSfAsLv+FXzy1g+botjD60B7ecNoisdrrnp+wuyMrs04Hb+c/KbAPc3TuFHJuIhKRqQ6G9O7TmsYsO47gBXaMOS1JUkIrij8DZwCKvq3lFxLQyW6RuH37xNdc+v5AVRVs4N7cnN50+kE5tVEVI7YIkilVAfqonCdDKbJF4SsrK+cNry/nbOyvZt1MbnhgznG/33zvqsCQNBEkU1wEzzOxtdl+ZrbvJiqSJeZ9/xbXPLWTlhq2cN7wX404dQEdVERJQkETxW2AL0AbQXDmRNLJ9Rzn3/GsZj7z7Kd2y2vLUJYdzdD+1JZX6CZIourn74NAjEZGE+uDTr7ju+QV8tnEbFx7RixtOGUiH1vWdES8SLFHMMLOT3P1foUcjIo22bcdO7n5lGU+89xk99mjLPy87nG/trypCGi5oK9RfmVkpUIamx4qkrNkrN3Ld8wv54qtt/PjI3lw3cgDtVUVIIwW5KWDHZAQiIg23tXQnd72ylCff+5zee7XjmbFHcETfvaIOSzJEkAV3x9a03d1nJT4cEamvf6/YwHWTF7L6m+1cfFQ21558AO1aqYqQxAnybbq2yu9tgOFUtkVVPwqRCG0p3ckdM5bwj/e/oE+X9jz7kyM5LHvPqMOSDBTk0tOoqs/NrCeVq7VFJCL/+/EGrp+8kDXF27n06D788qQDaNuqedRhSYZqSH1aCAxMdCAiUrfNJWX8bsYSnv5gFX33bs/zlx/Job1VRUi4goxR/IXK1qcAzYCDAfXLFkmyt5ev58bJC1m7qYSfHNuXX5zYnzYtVUVI+IJUFHOr/L4TeNrd3w0pnhqZWXvgbeBWd5+ezHOLRK14exm/fWkxz84tJKdrByZf8S0O6bVH1GFJExJkjOKJhh7czB4FTgeKqq7uNrORwJ+A5sDD7n5nHYe6Hni2oXGIpKs3l65j3JR8ijaX8P9G7M9V3+mnKkKSLux+FI8D9wFPVjlec+B+4EQqxzvmmNk0KpPGHdXePwYYCiymcsaVSJNQvK2MCdMXM3l+If336cBDPzyKoT07Rx2WNFGh9qNw91lmll1t83BghbuvBDCzZ4Az3f0OKquP3ZjZCKA9MAjYbmYz3L2iPnGIpJPXF69j3NRFbNy6g58dn8OVx+fQuoWqCIlOFP0ouseOuUshcHhtO7v7TQBmdhGwobYkYWZjgbEAvXr1SlCoIsnz9dYd3JZXwAsfrWHAvh159KLDGNw9K+qwRNKnH4W7P17H65OASQC5ubkp32RJpKpXC9Zy09R8vtm2g59/px8/PS6HVi2aRR2WCBBNP4rVQM8qz3vEtjWaWqFKuvlq6w7GTysgb8EaBu3XiSfGHMaB3VRFSGqJoh/FHKCfmfWhMkH8ADg/EQdWK1RJJzMWfcktL+SzqaSMa07szxUj9qdlc1URknpC7UdhZk8DI4AuZlYIjHf3R8zsSuBVKmc6PeruBfU9tki62rCllPEvFvDSoi8Z3L0T/zjncAbsq7v2S+qyusaozWwzlbOOUr4fRZVLT5d9/PHHUYcjsht3Z/rCLxk/rYAtJTv5+Qn9GHtsX1URkjLMbJ675/7X9sRNZkodubm5Pnfu3Lp3FEmS9ZtL+fWL+bycv5ahPbKYeM5Q+u+jVi+SWmpLFLVeejKzAe6+1MyG1fS6u+t+TyJ1cHfyFn7J+Bfz2VpaznUjD2DsMX1poSpC0ki8MYprqFyXcE8Nrzkp2I9Cs54klazfXMotL+TzSsFahvbszO9HD6GfqghJQ7r0JJJg1auIX5zYn8uO6aMqQlJebZee6vzmmtntsfsz7XreycweS3SAIplg/eZSLn9qHlc9/SG99mrPS1cdzRUj9leSkLQWZHpsC+ADM7sY2IfKm/z9JdSoRNKMuzNtwRrGTytg245ybjhlAJcerSpCMkOQ24zfaGavA+8DXwPHuvuK0CNrAI1RSBSKNpdwywv5vFqwjoN7dub35wwhp6vGIiRzBFlHcSzwAPAUcBCwB3CJu68JP7yG0RiFJEP1KuKXJ/bn0mP60ryZRR2aSIPUe3psFb8HznH3xbEDnQ28CQxIbIgi6aNocwk3T83nX4tVRUjmC5IojnT38l1P3H1K7E6yIk1O9SrixlMGqIqQjBdkjKLczE4DDmT3LnMTQouqgTRGIWGqWkUc0qszE0cPJadrh6jDEgldkFaoDwLtgOOAh4HRwAchx9UgunushMHdefGjyipie1k5404dwCVHq4qQpiPIpadvufsQM1vo7reZ2T3Ay2EHJpIKijaXcNPUfF5TFSFNWJBEsT32c5uZdQM2AvuFF5JI9KpWESVl5dx06kDGHN1HVYQ0SUESxXQz6wxMBOZTeZ+nh0ONSiRCRZtKGDc1n9eXrGNYr85MPGco+++tKkKariCD2bfHfp1sZtOBNu5eHG5YIsnn7rzw0WpunbZYVYRIFUEGs5sDpwHZu/Y3M9z93nBDqz/NepKGUhUhUrsgl57ygBJgEVARbjiNo1lPUl/Vq4ibTxvIxUepihCpKkii6OHuQ0KPRCTJqlYRh/beg7tHD1EVIVKDIIniZTM7yd3/FXo0Ikng7kz9cDW3TiugdGeFqgiROgRJFLOBqWbWDCgDDHB37xRqZCIhWLephJumLuL1JUUc2nsPJo4eQl9VESJxBUkU9wJHAos8E9vhSZOgKkKk4YIkilVAvpKEpKt1m0oYN2URbywtIjc2FqEqQiS4IIliJfCWmb0MlO7aqOmxkurcnSnzV3NbnqoIkcYIkig+jT1axR4pS9NjZZfqVcTEc4bSp0v7qMMSSUtBVmbfloxARBKhahWxo7yCW04fxEXfylYVIdIIQSoKkbSwblMJN05ZxJtLizgsew/uHq0qQiQRlCgk7bk7k+evZkKsivh1rIpopipCJCGUKCStrS0uYdxUVREiYQpyU8C9gcuoclNAAHcfE15YIvHtqiJuyyugTFWESKiCVBQvAu8ArwPl4YYjUre1xSXcOGUhM5etZ3j2ntw9egjZqiJEQhMkUbRz9+tDj0SkDu7O8/MKmTB9MWXlFYwfNYgfH6kqQiRsQTvcneruM0KPRqQWu1URffbk7u+pihBJliCJ4ufAODMrJcVvCqiV2ZnH3XluXiG3T1/MznLn1lGD+JGqCJGksky8hVNubq7PnTs36jCkkapXERNHD6H3XqoiRMJiZvPcPbf69lorCjMb4O5LzWxYTa+7+/xEBiiyi6oIkdQS79LTNcBY4J4aXnPg+FAikiZNVYRI6qk1Ubj72NjP45IXjjRVqiJEUpdWZkvkVEWIpDYlComMqgiR9BA3UZiZAT3cfVWS4pEmQlWESPqImyjc3c1sBnBQkuKRDKcqQiT9BLn0NN/MDnP3OaFHIxnty+Lt3DhlEW+pihBJK0ESxeHABWb2ObCV/6zMHhJqZJIxVEWIpLcgieLk0KOQjKUqQiT9BemZ/bmZDQWOiW16x90XhBuWpDtVESKZo1ldO5jZz4F/AF1jj6fM7GdhB1bl/CPM7B0ze9DMRiTrvNJwXxZv5+LH53Dd8wsZuF8nXrn6GC46qo+ShEiaCnLp6RLgcHffCmBmdwHvAX+p641m9ihwOlDk7oOrbB8J/AloDjzs7nfGOYwDW4A2QGGAeCUiqiJEMlOQRGHs3tmuPLYtiMeB+4An/+9gZs2B+4ETqfyLf46ZTaMyadxR7f1jqLzU9baZ7QPcC1wQ8NySRBqLEMlcQRLFY8D7ZjY19vws4JEgB3f3WWaWXW3zcGCFu68EMLNngDPd/Q4qq4/afA20DnJeSR5VESKZr66V2c2A2cBbwNGxzRe7+4eNOGd3oOpK70Iqp+DWFsPZVM686kxldVLbfmOpvNstvXr1akR4EpSqCJGmoa6V2RVmdr+7HwJE0n/C3acAUwLsNwmYBJWNi8KOqylTFSHStAS59PSGmX0PmOKJaYe3GuhZ5XmP2LZGUyvU8KmKEGl66myFamabgfbATqCEevbMjo1RTN8168nMWgDLge9QmSDmAOe7e0HDPsJ/UyvUxKteRVw/8gBVESIZpt6tUGNvagaMdPd3G3jSp4ERQBczKwTGu/sjZnYl8CqVM50eTWSSkMRTFSHStAUZo7gPOKQhB3f382rZPgOY0ZBjxqNLT4mlsQgRgQArs4mNUcR6U6Q0d89z97FZWVlRh5L2tLpaRHYJMpj9E+AaoNzMtlPPMQpJL6oiRKS6IDcF7JiMQBJBl54aR2MRIlKTIDcFNDO70MxuiT3vaWbDww+t/nTpqWHcnWfnrOKke2fx/sqvuHXUIJ657AglCREBgl16+itQARwP3E7lDfruBw4LMS5Jki+Lt3PD5EW8vVxVhIjULFCHO3cfZmYfArj712bWKuS4JGTuznNzY2MRFRqLEJHaBUkUZbE7vjqAme1NZYWRcjRGEYyqCBGpjyCJ4s/AVKCrmf0WGA3cHGpUDeTueUBebm7uZVHHkopURYhIQwSZ9fQPM5tH5S03DDjL3ZeEHpkklKoIEWmoIBUF7r4UWBpyLBICVREi0liBEoWkJ1URIpIIGZUoNJhdSVWEiCRSRiUKDWarihCRxMuoRNGUqYoQkbAoUWQAVREiEiYlijSmO72KSDJkVKJoSoPZutOriCRLnT2z01Em98xW72oRCUuDemZLalEVISJRUKJIAxqLEJEoKVGkOFURIhI1JYoUpSpCRFKFEkUKWltcwo1TFjJTVYSIpICMShTpPj1WVYSIpKKMShTpfK+n3aqI7D2ZeI6qCBFJDRmVKNJR1SqirLyC8aMG8WNVESKSQpQoIlS9irh79BCyu6iKEJHUokQRAXfn+XmFTFAVISJpQIkiyVRFiEi6UaJIElURIpKulCiSQFWEiKQzJYoQqYoQkUygRBGStcUljJu6iDeXFqmKEJG0llGJIhVWZrs7k+ev5ra8AsrKK/j16YO46FuqIkQkfWVUooh6ZbaqCBHJRBmVKKKiKkJEMpkSRSOpihCRTKdE0UCqIkSkqVCiaIB1m0q4cYqqCBFpGpQo6mFXFTEhr4AdqiJEpIlQogioahVxWPYe3D16KH1URYhIE6BEUQdVESLS1ClRxKEqQkREiaJGqiJERP5DiaKadZtKGDdlEW+oihARAZQodjP1w0LGv6gqQkSkqpRPFGbWDLgd6ATMdfcnwjrXiqItHLBvR1URIiJVhJoozOxR4HSgyN0HV9k+EvgT0Bx42N3vjHOYM4EewEagMMRwufqE/jQ3UxUhIlJF2BXF48B9wJO7NphZc+B+4EQq/+KfY2bTqEwad1R7/xjgAODf7v6QmT0PvBFWsC2bNwvr0CIiaSvUROHus8wsu9rm4cAKd18JYGbPAGe6+x1UVh+7MbNCYEfsaXl40YqISE2i+Cd0d2BVleeFsW21mQKcbGZ/AWbVtpOZjTWzuWY2d/369YmJVEREUn8w2923AZcE2G8SMAkgNzfXw45LRKSpiKKiWA30rPK8R2ybiIikoCgSxRygn5n1MbNWwA+AaYk4sJmNMrNJxcXFiTiciIgQcqIws6eB94ADzKzQzC5x953AlcCrwBLgWXcvSMT53D3P3cdmZWUl4nAiIkL4s57Oq2X7DGBGos9nZqOAUTk5OYk+tIhIk2XumTfua2brgc8beZgsoDHXsBr6/vq8L+i+de1X1+tdgA0BY0p1jf1zTZVzJuKYDTlGFN/PuvbR9zNx5+3t7nv/16vurkcND2BSFO+vz/uC7lvXfgFenxv1n0eq/LmmyjkTccyGHCOK72dd++j7Gf55tRS5dnkRvb8+7wu6b137NfazppMoPmsY50zEMRtyjCi+n/U9bzqL6nPGPW9GXnqSxDKzue6eG3UcIjXR9zN8qigkiElRByASh76fIVNFISIicamiEBGRuJQoREQkLiUKERGJS4lCGsXMzjKzv5nZ/5jZSVHHI1KVmfU1s0diTc+kgZQomjAze9TMiswsv9r2kWa2zMxWmNkN8Y7h7i+4+2XA5cC5YcYrTUuCvp8r3b3ONgUSn2Y9NWFmdiywBXjSYz3NY61ql1OlVS1wHrW0qnX3otj77gH+4e7zkxS+ZLgEfz+fd/fRyYo906R84yIJjyemVa0BdwIvK0lIIiXi+ymJoUtPUl19W9X+DDgBGG1ml4cZmAj1/H6a2V5m9iBwiJndGHZwmUoVhTSKu/8Z+HPUcYjUxN03Ujl+Jo2gikKqU6taSWX6fkZAiUKqC61VrUgC6PsZASWKJizZrWpF6kPfz9Sh6bEiIhKXKgoREYlLiUJEROJSohARkbiUKEREJC4lChERiUuJQkRE4lKiEInDzK42s3YhHr9bInolmNmWRMQjUhOtoxCJw8w+A3LdfUPUscRjZlvcvUPUcUhmUkUhAphZezN7ycwWmFm+mZ1rZlcB3YCZZjYztt9JZvaemc03s+fMrENs+2dmdreZLTKzD8wsp4ZzfNvMPoo9PjSzjmaWvasxj5k9XOX19WY2Prb9WjObY2YLzey2OJ/ht7H4Z5vZPmH8d5KmSYlCpNJIYI27D401yXkldmfcNcBx7n6cmXUBbgZOcPdhwFzgmirHKHb3g4D7gD/WcI5fAT9194OBY4DtVV9090tjr50JbAAej7WX7UdlH4aDgV8NQgoAAAGhSURBVENjDX2qaw/MdvehwCzgsob9ZxD5b0oUIpUWASea2V1mdoy7F9ewzxHAIOBdM/sI+DHQu8rrT1f5eWQN738XuDdWqXSO3bdoN2bWBngO+Jm7fw6cFHt8CMwHBlCZOKrbAUyP/T4PyI7zWUXqRf0oRAB3X25mw4BTgd+Y2RvuPqHabga85u7n1XaYWn7fdY47zeyl2DneNbOTgZJquz0ITHH316uc8w53f6iOj1Dm/xlwLEf/b0sCqaIQoXL2EbDN3Z8CJgLDYi9tBjrGfp8NHLVr/CE2rtG/ymHOrfLzvRrOsb+7L3L3u6i8XfaAaq//FOjo7ndW2fwqMKbKWEh3M+vaiI8qUm/6V4dIpYOAiWZWAZQBV8S2TwJeMbM1sXGKi4Cnzax17PWbgeWx3/cws4VAKVBT1XG1mR0HVAAFwMvAflVe/xVQFrusBfCguz9oZgOB9yrbk7MFuBAoavQnFglI02NFEiBdptGKNIQuPYmISFyqKEREJC5VFCIiEpcShYiIxKVEISIicSlRiIhIXEoUIiISlxKFiIjE9f8BYftgemwsmVQAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] @@ -440,7 +435,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEKCAYAAAASByJ7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAlFUlEQVR4nO3dd3xUdb7G8c83DQgldJBekiAoUqQpSFexsLq7ruBaaKKyKmBh7+7qrvdu0ysW7Ig0QcRV71rWXQUL0qugoggJhCKi0gMB0n/3j0zcyIZkkszkzGSe9+uVVzJnZs55oqNPzjm/8zvmnENERCTK6wAiIhIaVAgiIgKoEERExEeFICIigApBRER8VAgiIgJAjNcBysPMhgPDa9euPT45OdnrOCIiYeWTTz456JxrdPpyC+frEHr06OE2bNjgdQwRkbBiZp8453qcvlyHjEREBFAhiIiIjwpBREQAFYKIiPioEEREBFAhiIiIT1gWgpkNN7MZ6enp5Xr/J7sPs/W7YwFOJSIS3sKyEJxz/3DO3ZKQkFCu9z/4r60Mm7acCS99wlffqhhERCBMC6GiZo7qwcTBiSxPPchlTyzntvmfsGWfikFEIltEX6l89GQ2s1fsZM7KXRzPymXYOU2ZOCSJTs3qBDCliEhoOdOVyhFdCIXST+Ywa+VO5qzYyfGsXC49pwkThyRxTrPyHZISEQllKgQ/pJ/MYfbKncxeuZPjmblc0qkJk4aqGESkalEhlEH6qRzmrNzJrBUFxXBxpyZMGpLEuc1VDCIS/lQI5ZB+Koe5K3cxc0UaxzNzGdqxCZOHqhhEJLypECqgsBhmrUjjmIpBRMKcCiEAjmX69hiWFxZDYyYNSaZzCxWDiIQPFUIAHcvM4cWVu5i5Yifpp3IYcnZjJg1N4rwWdSs9i4hIWVWpQii8hWZiYuL41NRUz3Icz8zhxVW7eGF5QTEMPrsxk4Yk0aVlXc8yiYiUpkoVQqFQuYXm8cwc5q3ezQvL0zh6ModBHRoxaWgyXVUMIhKCVAiVQMUgIuFAhVCJMrJyeXFVwcnnIydzGNihEZOGJNGtVT2vo4mIqBC8kJGVy7zVu3hhWUExDEhuxKShSXRXMYiIh1QIHsrIymX+6t3MWLaDIydz6J9csMdwfmsVg4hUPhVCCDiRlcv8NbuZsSyNwyeyuSipIZOHJqsYRKRSqRBCSPHFkMT5ret7HU1EIoAKIQSdyMrlJV8xHPIVw6QhSfRoo2IQkeBRIYSwk9kFxfD80oJi6JdYsMegYhCRYFAhhIGT2bksWLOH55ft4GBGNn0TGzB5aDI9VQwiEkAqhDByKjuPBWt3M31pQTFc2L6gGHq1VTGISMWpEMLQv4shjYMZWVzYvgGThiTRu10Dr6OJSBhTIYSx04vhgnYNmDQ0iT4qBhEphypVCKEy22llO5Wdx8vr9jB96Q4OHM+iT7v6TB6arGIQkTKpUoVQKFL2EE6XmZPHy2v38FyRYrj74g46xyAiflEhVEGFxfDsxzs4mJFFv8SG3HWxrnwWkZKpEKqwwnMMz328g0MnshmQ3Ii7Lta02yJSPBVCBDiZncu81bt5fmnBJHpDzm7MXRcnc25z3fNZRP5NhRBBCu/HMGNZGumncrikUxMmD02mU7M6XkcTkRCgQohAxzJzmLNiFzNXpHE8M5fLOzdl0pBkOjSt7XU0EfGQCiGCpZ/MYdaKNGav3MWJ7FyuPK8Zk4YkkthYxSASiVQIwpET2bywPI25q3ZxKiePq7o0Y+KQJNo1quV1NBGpRCoE+cGhjCxmLE9j3qrdZOXm8bPuLZg4OIlWDeK9jiYilUCFIP/hwPEsnl+6g/lrdpOb77imewvuGJxIy/oqBpGqTIUgZ7T/WCbPfryDl9ftIT/fcW3PltwxKJFmdWt4HU1EgkCFIKX6Nv0Uzy7ZwSvr92AYI3u15FcDE2maUN3raCISQCoE8ds3R0/x9EfbeW3D10RFGdf3bsWEge1pXFvFIFIVqBCkzL4+fJKnP9rO6xv3Ehtt3NinNbcOaE/DWtW8jiYiFaBCkHLbdfAET36UypubvqFaTDSjLmzDLf3bUb9mnNfRRKQcVAhSYTsOZPDUh6m89dk+4mOjGd23DeMvakfdeBWDSDhRIUjApH5/nCc+TOWdz7+ldrUYxvRry7h+bUmoEet1NBHxQ4ULwczqADGFj51zhwMXr3xUCN7a+t0xnvgglXe/+I7a1WMYf1E7xvRtQ+3qKgaRUFbuQjCzW4E/AqeAwhc751y7gKcsIxVCaPhyXzrTPkjl/S3fUzc+lvEXtWP0hW2oWS2m9DeLSKWrSCGkAhc45w4GK1x5qRBCy+a96Tz+QQofbd1P/Zpx3Nq/HTde0Jr4OBWDSCg5UyFE+fHeHcDJwEeSqqZziwRmj+7Jm7f3pXPzBB58dyv9H17CzOVpZObkeR1PRErhzx5CN2AOsBbIKlzunJsY0CBmVwNXAI2BZ5xzi0t7j/YQQtsnuw/z+PuprNh+kEa1q3H7wPaM7NWK6rHRXkcTiWgVOWS0DlgBbAbyC5c75170Y6OzgSuB/c65c4ssHwY8AUQDM51zDxV5rh7wiHNuXGnrVyGEh7Vph3js/RTW7jxM0zrVuX1Qe67t2ZJqMSoGES9UpBBWOecuLOdG+wMZwLzCQjCzaCAFuBjYC6wHrnPObfE9/yiwwDm3sbT1qxDCy6odB3lscQobdh+hWUJ17hicxDXntyAuxp8jlyISKBU5h7DEzG4xs7PMrH7hlz8bdc4tA04fntoL2O6cS3POZQOvAFdZgf8F3i2pDHxZNpjZhgMHDvgTQ0LEhe0b8tptFzB/XC+aJFTnd29sZvCjH/Pq+q/JycsvfQUiElT+7CHsLGax38NOzawN8E6RPYRrgGHOuZt9j28EelOw1zCKgj2GT51z00tbt/YQwpdzjo9TDvD4+yl8vjed1g3imTg4iau7NSc6yryOJ1KlnWkPocTxgGYWBfzGOfe3QGYpZplzzj0JPBnA7UgIMzMGdWjMwORGfPjVfh7/IIV7XvuM55bu4J6Lkxl2blPMVAwilanEQ0bOuXzg9gBvcy/QssjjFsC+AG9DwoSZMbRTE965sx/PXd8dgAkLNnLVMytZnnqAcJ5aRSTc+HMO4X0zu9fMWpb1HMIZrAeSzKytmcUBI4G3y7ICMxtuZjPS09MrEENCiZlxWeezWDS5P1OvOY9DGdncOGsd172who17jngdTyQiBPUcgpktBAYCDYHvgQecc7PM7HJgGgXDTmc75/5SxtyAziFUZVm5eSxcu4enl2znYEY2Qzs2YcqlHejQtLbX0UTCnmY7lbB0IiuXOSt38vyyNDKycrm6a3PuGppMqwbxXkcTCVsVuQ4hFpgA9Pct+hh43jmXE+iQZaVCiBxHT2YzfWkac1ftJDfPMbJXSyYOTqJxHd3WU6SsKlIIM4FYoPDK5BuBvMJho14ws+HA8MTExPGpqalexRAPfH8sk6c+SuWVdV8TE22MvrAttw3QTXpEyqIihfCZc65Lacu8oD2EyLXn0Eke/yCFNz/9hlrVYri1fzvG9G2rKbdF/FCRK5XzzKx9kRW1AzR1pXiqVYN4Hh/RlXcnXUSfdg14ZHEKA6YuYc7KnWTl6uMpUh7+7CEMoWC20zQKLiprDYxxzi0JfrySaQ9BCm3cc4Sp721jddohmtetweShSfysewtd9SxSjAqNMjKzakAHCgphq3Muq5S3VAoVghTlnGPF9oNMXbSNz/emk9i4lq56FilGRQvhQqANP76n8rxABiwLnVSWkjjnWPTldzyyOIXt+zM4r0UCUy7tQL/EhioGESp2Unk+0B74lH+fO3CBvkFOeWgPQUqSl+/4+8a9TPsglW+OnqJPu/r8etjZdG9Vz+toIp6qSCF8BXRyIXgFmwpB/FHcVc/3XprM2U3reB1NxBMVGWX0BdA08JFEKke1mGhG923L0imDuPeSZNbuPMRlTyxn8iub2H3ohNfxREKGP3sIS4CuwDp+fE/lnwQ1mR+0hyDlcfpVzyN6tmTikCSa6KpniRAVOWQ0oLjlzrmlAcpWZjqpLIGw/1gmTxa56nnUhW2YMKC9rnqWKk+T24mcga56lkijQhApxdbvjvHo4hTe3/I9DWvFcfugRH7ZuxXVYqK9jiYSUCoEET8Vd9XzT7s1JybanzEYIqGvIqOMRCJK91b1eHl8b14a15sGteKY8vrnXDptGe9u/la39JQqrdRCMLO+Zva+maWYWZqZ7TSztMoIJ+IVM6NfUkPeur0v02/ojpkxYcFGfvL0Spal6F7PUjX5M8poK3AX8AlFZjl1zh0KbrTS6ZCRVJa8fMcbm77h8fdTdNWzhL2KDDtd65zrHbRk5aBhp+IVXfUsVUFFCuEhIBr4Oz++MG1joEOWlfYQxCsnsnKZu2oX05fuICMrl6u6NOOui5Np3aCm19FESlWRQijuvgfOOTc4UOHKS4UgXtNVzxKONOxUJIj2H8vkqY+2s3Ddnh/u9fyrQe2pUz3W62gi/6Gi90O4AjgH+OHPHufcHwOasBxUCBJqil71XD8+jrsvSWZkz1a6c5uElHJfh2Bm04ERwJ0U3DHtFxTcRlNETlN4r+e3b+9H+0a1uO+NL7jiyeWs3H7Q62gipfLnwrQLnXM3AUecc/8DXAC0DG4skfDWuUUCf7u1D89e352MrFyun7mWm1/cwM6Dmm5bQpc/hXDK9/2kmTUDcoC2wYskUjWYGZd3PosP7h7Afw07m9U7DnLJ40v58ztbSD+V43U8kf/gTyG8Y2Z1ganARmAX8EoQM5XKzIab2Yz09HQvY4j4pXpsNBMGtmfJlIH8vHsLZq3cycCpS5i/ehe5eflexxP5QZlGGZlZNaC6cy4k/k+sk8oSjr7cl86f3tnCmrTDJDWuxe+v7ET/5EZex5IIUpGTyvFm9nsze8E5lwU0NrMrg5JSJAKc0yyBheP78PyN55Odl89Ns9cxdu56tu/P8DqaRDh/DhnNoeAK5Qt8j/cCfw5aIpEIYGZcek5TFt/Vn99dfjbrdx5m2LRl/PfbX3L0ZLbX8SRC+VMI7Z1zD1NwMhnn3CkKhp+KSAVVi4nmlv4F5xeu7dmSeat3MWDqx8xduZMcnV+QSuZPIWSbWQ3AAZhZe4rMaSQiFdewVjX++tPO/GvSRXRunsB//2MLw6YtY8m2/V5HkwjiTyE8ALwHtDSzBcCHwK+DmkokQp3dtA7zx/Vi1qge5DsYM2c9N81eR8r3x72OJhHA36krGgB9KDhUtMY5FxKXXWqUkVRl2bn5zF+zmyc+SOFEdh7X927F5KHJ1K8Z53U0CXNlnsvIzLqXtEJNfy1SOQ6fyGbaByksWLuHmnHRTBySxE0XtCEuRnfAlfIpTyEUTntdHegBfEbBHsJ5wFrnXL8gZfWbCkEiSer3x/nTP79iWcoB2jasyX2Xd2RIx8aYaYyHlE2Zr0Nwzg1yzg0CdgPdnXM9nHPnA92A7cGLWjpdqSyRKKlJbeaN7cWcMT2JMrh53gZunLWOrd8d8zqaVBH+3CDnU+dc19KWeUF7CBKpcvLyWbBmN49/kMrxzBxG9mrF3Rcn07BWNa+jSRioyB3TFgIngJcoGHp6A1DLOXddMIKWhQpBIt3Rk9k88WEq81fvpkZsNHcOSWTUhW2oFhPtdTQJYRUphOrABKC/b9Ey4DnnXGbAU5aRCkGkwPb9Gfz1X1/x0db9tG4Qz28v68il5zTR+QUplm6hKRIBlqYc4M/vbCF1fwZ92tXn91d24pxmCV7HkhBTnlFGrzrnrjWzzfiuUi7KOXde4GOWjQpB5D/l5uWzcP3XPLZ4G0dP5TCiR0vuviSZxrWrl/5miQjlKYSznHPfmlmxt8t0zu0OcMYyUyGInFn6qRye+jCVuat2US0mitsHJzK2b1uqx+r8QqTTISORCJV2IIO//msrH3z1PS3q1eB3l3fksnOb6vxCBCvzdQhmdtzMjhX5fqzo4+DGFZFAadeoFjNH9WDBzb2pVS2GXy3YyIjn17B5r67jkR/THoJIBMnLd/xt/dc8ungbh09m8/PuLZhyaQea1NH5hUhSkTumjStm2UOBCiYilSc6yvhl71YsmTKQW/q34+1P9zHokY95+qNUMnPyvI4nHvNndqxrzOz6wgdm9izQOHiRRCTY6lSP5beXdeT9u/vTP6kRjyxOYcijS3n7s32E81EDqRh/CuFnwGgzu87M5gHZzrmxQc4lIpWgdYOaTL/xfBaO70NCjVgmLtzENdNX8+nXR72OJh4oadhp/SIPawNvAiuBPwA45w4HO1xpdA5BJHDy8h3/98leHl60jYMZWfysW3OmDOvAWQk1vI4mAVae6xB2UnBBmhX5Xsg559oFI2hZqBBEAi8jK5dnl2xn5oqdRBlMHJLErf3bEx2lYapVRZW6DsHMhgPDExMTx6empnodR6RK+vrwSf7yz69478vv6NmmHo9d25WW9eO9jiUBUJHJ7W4qbrlzbl6AspWb9hBEgss5x5uffsMf3vwSBzwwvBPXnN9CF7WFuXIPOwV6Fvm6CPhv4CcBTSciIcnM+Gm3Frw7+SLOaVaHKa9/zoSXNnL4RLbX0SQIynzIyMwSgPnOOc9LQXsIIpUnL98xc3kajyzeRt34OKZecx4DO2gEejiqyB7C6U4CSRWPJCLhJDrKuHVAe966vR/14mMZPWc9f3jrC05l64K2qiKmtBeY2T/49/TXUUAn4LVghhKR0NWpWR3evqMfUxdtY9aKnazcfpBpI7rRuYXuuxDu/DmpPKDIw1xgt3Nub1BT+UmHjES8tXL7Qe559TMOZmQxeWgSEwYmanhqGCjPbKdnAzjnlgJrnHNLnXMrnXN7zaxPELOKSJjom9iQRZP7c1nns3hkcQrXPr+aPYdOeh1LyqmkcwgvF/l59WnPPRuELCIShhLiY3nqum48MbIrKd8f57InlvHqhq81J1IYKqkQ7Aw/F/dYRCLcVV2b897k/nRukcCvX/+c2176RMNTw0xJheDO8HNxj0VEaF63Bi/f3If7Lu/Ikq0HuHTaMpZs2+91LPFTSaOMWpjZkxTsDRT+jO9x86AnE5GwFBVljO/fjn5JDZn8yqeMmbOeG/u05neXd6RGnO7nHMpKKoQpRX4+fSiPhvaISIk6nlWHt+7oyyOLtjFzxU5W7jjItBFdOa9FXa+jyRmE5eR2hTTsVCQ8rNp+kHte+4wDxwuGp942oD0x0eW5LlYCIZBXKouIlMmFiQ15b1J/LvcNTx0xY42Gp4YgFYKIVIqE+FiePH146noNTw0lKgQRqVSFw1PPa1GXX//f59w6/xMOZWR5HUvwb+qKRsB4oA1FTkKHwn2VdQ5BJHzl5ztmrdjJ1EXbqFMjlqnXnMegszV7amWoyDmEt4AE4APgn0W+RETKrXB46lt39KVhrTjGzF3P/W9u1uypHip1tlMg3jn3X0FPIiIRqeNZdXjz9r48urhgeOqq7Yd4fERXurSs63W0iOPPHsI7ZnZ50JOISMSqHhvNfVd0YsHNvcnMyePnz63iyQ9Tyc3L9zpaRPGnECZRUAqnzOyYmR03s2OBDmJm7cxslpm9Huh1i0h4uLB9Q96d3J8rzjuLx94vmD1196ETXseKGKUWgnOutnMuyjlXwzlXx/e4jj8rN7PZZrbfzL44bfkwM9tmZtvN7De+7aQ558aV79cQkaoioUYsT4zsxpPXdWP7/gwue2I5r6zbo+GplaDU+yGYWffivvxc/1xg2GnrjQaeAS6j4O5r15lZp3KlF5Eq6yddmvHe5P50bVmX3/x9M7doeGrQlXRS+W7gFuDRYp5zwODSVu6cW2ZmbU5b3AvY7pxLAzCzV4CrgC3+BDazW3y5aNWqlT9vEZEw1axuDV4a15vZK3fy8KJtXDptOQ9f05nBZzfxOlqVdMY9BOfcLb7vg4r5KrUMStAc+LrI471AczNrYGbTgW5m9tsScs1wzvVwzvVo1KhRBWKISDiIijJuvqgd/7ijHw1rxTF27gbue2MzJ7NzvY5W5fgz7DTQiru5jnPOHQJuq+wwIhIeOjStzVt39OWxxSnMWJ7G6h0anhpoXkxdsRdoWeRxC2CfBzlEJMxUi4nmt5d35OWb+5CZk8fPNDw1oEosBCvQsqTXlMN6IMnM2ppZHDASeLssKzCz4WY2Iz09PcDRRCQcXNC+Ae9O7s9w3/DUX2h4akCUWAiuYJzXm+VduZktBFYDHcxsr5mNc87lAncAi4CvgFedc1+WZb3OuX84525JSEgobzQRCXMJNWKZNrIbT13XjR0anhoQ/kxu9www1zm3vnIi+U+T24kIwLfpp7jn1c9YteMQQzs24aGfd6ZhrWpexwpZFZncbhCw2sx2mNnnZrbZzD4PfEQRkfI5K6FgeOr9V3RkWeoBhk1bxkdbv/c6VtjxZw+hdXHLnXO7g5LID2Y2HBiemJg4PjU11asYIhKCtn13nEmvbGLrd8f5Ze9W3H9FR+LjvBhQGbrOtIfg1z2VzawLcJHv4XLn3GcBzlcuOmQkIsXJys37YXhqmwY1eXxEV7pqeOoPyn3IyMwmAQuAxr6vl8zszsBHFBEJjMLhqQvH9yE7N5+fP7eKaR+kaHhqKfw5ZPQ5cIFz7oTvcU1gtXPuvErIVyLtIYhIaY5l5vDAW1/yxqZv+EmXZkwb0ZWoqOKuj40cZ9pD8OfAmgFFb2GUR/FXG4uIhJw61WN5fERXkprU4uH3tlG/ZhwPDO+Emf43djp/CmE2sNbM3vA9vhqYFbREfihyUtnLGCISRiYMaM+hjGxmrdhJw1px3DE4yetIIae0K5WjgLXAGOAwcAQY45ybFvxoZ6YL00SkrMyM+y7vyE+7NeeRxSksXLfH60ghp8Q9BOdcvpk96py7ANhYSZlERIIiKsp4+JrzOHIym/ve2Ey9+FiGnXuW17FChj8Xpi02s5+bDriJSBUQGx3Fs9d3p0vLukxc+Cmrdhz0OlLI8KcQ7gZeA7KCeU9lEZHKEh8Xw5zRPWndIJ5b5n3CF99ookzw7xzCMN89lePKek9lEZFQVTc+jnnjelGnegyj56xj10HNllrabKf5wCOVlMVvmv5aRALhrIQazBvXm7x8x02z17H/WKbXkTwVlucQNMpIRAIlsXEt5ozpxcGMLEbNWc+xzByvI3mmLOcQsnUOQUSqoq4t6zL9hvPZvv84N7+4gcycvNLfVAWVWgi+cwZRzrlYnUMQkaqqf3IjHvlFF9bvOszEhZsict4jfya3MzO7wcx+73vc0sx6BT+aiEjluqprcx64shOLt3zPfW98EXF3X/PnkNGzwAXAL32PM4BngpZIRMRDo/u25c7Bifxtw9c8snib13EqlT9zGfV2znU3s00AzrkjZhYX5Fwl0lxGIhJMd1+czMGMbJ5ZsoP6Nasxrl9bryNVCn/2EHLMLBpwAGbWCPD04JpGGYlIMJkZf776XIad05Q/vbOFNzd943WkSuFPITwJvAE0NrO/ACuAvwY1lYiIx6KjjGkju9KnXX3ufe0zPt623+tIQefPKKMFwK+BB4Fvgaudc68FO5iIiNeqx0bzwk096NC0NhNe2sjGPUe8jhRU/uwh4Jzb6px7xjn3tHPuq2CHEhEJFbWrxzJ3TC8a16nG2LnrSf3+uNeRgsavQhARiWSNaldj/tjexEZHcdPsdew7esrrSEGhQhAR8UOrBvG8OKYXGZm53DhrLUdOZHsdKeBUCCIifurUrA4zR/Xg6yOnGDN3PSeycr2OFFBhWQia7VREvNK7XQOeuq4bn+89yoQFG8nOrTpTXIRlIeg6BBHx0qXnNOXBn3VmWcoBprz+Gfn5VWOKC3+uVBYRkdOM6NmKgxnZTF20jXrxcTwwvBMhdJeAclEhiIiU068GtudQRjazV+6kUe1q3D4ovKfTUSGIiJSTmXH/FR05fCKLqYu20aBmHCN7tfI6VrmpEEREKiAqypj6iy4cPZXD797YTN34OIad29TrWOUSlieVRURCSWx0FM9e350uLesy8ZVNrEk75HWkclEhiIgEQHxcDLNH9aRV/XjGv7iBL/eF37B4FYKISIDUqxnHvLG9qF09hlGz17P70AmvI5WJCkFEJICa1a3BvHG9yM3P58ZZ69h/PNPrSH5TIYiIBFhi49rMGd2TA8ezGDV7Pccyc7yO5JewLARNXSEioa5bq3pMv/F8Ur8/zvgXN5CZk+d1pFKFZSFo6goRCQcDkhvx6LVdWLvzMBMXbiI3L7TnPQrLQhARCRdXdW3OA8M7sXjL99z/5hc4F7rzHunCNBGRIBvTty2HMrJ5esl2Gtaqxr2XdvA6UrFUCCIileCeS5I5dCKLp5dsp37NOMb2a+t1pP+gQhARqQRmxp+uOpfDJ7L54ztbaFArjqu6Nvc61o/oHIKISCWJiY7iiZHd6N22Pve8+hkfb9vvdaQfUSGIiFSi6rHRvDCqB0lNajPhpY1s2nPE60g/UCGIiFSyOtVjeXFsTxrVrsaYuevZvv+415EAFYKIiCca167O/HG9iImK4qZZ69h39JTXkVQIIiJead2gJi+O7cnxzFxumr2OIyeyPc2jQhAR8dA5zRJ4YVQP9hw+yZi56zmZnetZFhWCiIjH+rRrwJMju/H53qNMeGkjOR5NcaFCEBEJAcPObcpff9qZpSkHmPLaZ+TnV/4UF2F5YZqZDQeGJyYmeh1FRCRgRvZqxaET2UxdtI16NeP4w5WdMLNK235Y7iFotlMRqap+NbA9Y/q2Yc7KXTz78Y5K3XZY7iGIiFRVZsbvr+jEYd+eQoOacYzs1apStq1CEBEJMVFRxtRrunDkZA6/e2Mz9WrGcek5TYO/3aBvQUREyiwuJorpN3TnvBZ1uXPhJtakHQr6NlUIIiIhKj4uhjmje9KqfjzjX9zAln3Hgro9FYKISAirVzOOeWN7Uat6DDfNXseeQyeDti0VgohIiGtWtwbzxvYiNz+fG2ev5cDxrKBsR4UgIhIGkprUZvbonuw/lsWo2es4lpkT8G2oEEREwkT3VvV47obuHD6RzXfpmQFfv4adioiEkYEdGvPxlIFUj40O+Lq1hyAiEmaCUQagQhARER8VgoiIACoEERHxUSGIiAigQhARER8VgoiIACoEERHxMecq/76dgWJmB4DdFVxNApBeye8v63v8fX1pryvt+YbAwTLkCmUV/fcaKtsMxDqD/RkN1OeztNfo8xm47bZ2zjX6j2edcxH9Bcyo7PeX9T3+vr601/nx/Aav/32Eyr/XUNlmINYZ7M9ooD6fpb1Gn8/gb1eHjOAfHry/rO/x9/Wlva6iv2s48eJ3DcY2A7HOYH9GA/X5LOt2w5lXv2eJ2w3rQ0YSWGa2wTnXw+scIsXR5zP4tIcgRc3wOoBICfT5DDLtIYiICKA9BBER8VEhiIgIoEIQEREfFYL4xcyuNrMXzOwtM7vE6zwiRZlZOzObZWave50lnKkQIoCZzTaz/Wb2xWnLh5nZNjPbbma/KWkdzrk3nXPjgdHAiCDGlQgToM9nmnNuXHCTVn0aZRQBzKw/kAHMc86d61sWDaQAFwN7gfXAdUA08OBpqxjrnNvve9+jwALn3MZKii9VXIA/n687566prOxVTYzXAST4nHPLzKzNaYt7Adudc2kAZvYKcJVz7kHgytPXYWYGPAS8qzKQQArE51MCQ4eMIldz4Osij/f6lp3JncBQ4Bozuy2YwUQo4+fTzBqY2XSgm5n9NtjhqirtIUQuK2bZGY8fOueeBJ4MXhyRHynr5/MQoD9UKkh7CJFrL9CyyOMWwD6PsoicTp9PD6gQItd6IMnM2ppZHDASeNvjTCKF9Pn0gAohApjZQmA10MHM9prZOOdcLnAHsAj4CnjVOfellzklMunzGTo07FRERADtIYiIiI8KQUREABWCiIj4qBBERARQIYiIiI8KQUREABWCCGY22czig7j+ZoGYp9/MMgKRR+RMdB2CRDwz2wX0cM4d9DpLScwswzlXy+scUnVpD0EihpnVNLN/mtlnZvaFmY0ws4lAM2CJmS3xve4SM1ttZhvN7DUzq+VbvsvM/tfM1vm+EovZxgAz+9T3tcnMaptZm8Kbv5jZzCLPHzCzB3zLp5jZejP73Mz+p4Tf4S++/GvMrEkw/jlJ5FIhSCQZBuxzznXx3YjlPd8srvuAQc65QWbWELgfGOqc6w5sAO4uso5jzrlewNPAtGK2cS9wu3OuK3ARcKrok865m33PXQUcAub6bkmaRME9ALoC5/tuGnO6msAa51wXYBkwvuz/CETOTIUgkWQzMNT3V/5Fzrn0Yl7TB+gErDSzT4FRQOsizy8s8v2CYt6/EnjMt+dR1zcnz4+YWXXgNeAO59xu4BLf1yZgI3A2BQVxumzgHd/PnwBtzvyripSd7ocgEcM5l2Jm5wOXAw+a2WLn3B9Pe5kB7zvnrjvTas7wc+E2HjKzf/q2scbMhgKZp71sOvB359wHRbb5oHPu+VJ+hRz375N+eei/Xwkw7SFIxDCzZsBJ59xLwCNAd99Tx4Havp/XAH0Lzw+YWbyZJRdZzYgi31cXs432zrnNzrn/peBw09mnPX87UNs591CRxYuAsUXOVTQ3s8YV+FVFykV/YUgk6QxMNbN8IAeY4Fs+A3jXzL71nUcYDSw0s2q+5++n4IbvANXMbC0Ff0wVtxcx2cwGUfAX/BbgXeCsIs/fC+T4DkcBTHfOTTezjsDqgltXkwHcAOyv6C8sUhYadirip3AZnipSXjpkJCIigPYQRETER3sIIiICqBBERMRHhSAiIoAKQUREfFQIIiICqBBERMTn/wH5IaufEml2JgAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEKCAYAAAAMzhLIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3dd3xV9f3H8deHsFfYewQIQ1YQAihWrbhwL9TaWltxtv6sFVeptta9tdU6q9ZaR2WJII4iSrGIIiAJYe+9RyCQkPX5/ZFLGylcLiQ3596b9/Px4JF7T+54B655e+459/M1d0dERORQqgQdQEREYpuKQkREwlJRiIhIWCoKEREJS0UhIiJhqShERCSsqkEHiIYmTZp4SkpK0DFEROLKrFmztrp70wO3J2RRpKSkMHPmzKBjiIjEFTNbdbDteutJRETCSqiiMLPzzOyV7OzsoKOIiCSMhCoKd5/g7tcnJycHHUVEJGEkVFGIiEj5U1GIiEhYKgoREQlLRSEikgBy84t4/7u1RGPpiIT8HIWISGWRV1DEO9+s5oUpy9ias49OTevSu02Dcn2OhCoKMzsPOC81NTXoKCIiUZVfWMx7M9fw/OdL2bgrj+M7NubFK/uWe0kAWCKucJeenu76ZLaIJKKComLGzFrLc58vZd3OXNLbN2T4GV0Y1KlJmR/bzGa5e/qB2xNqj0JEJFEVFTvjvlvHs58vYdW2vaS1bcAjF/fixM5NMLOoPreKQkQkhhUXOx/O3cAfP1vM8i176NGqPq/9LJ3B3ZpFvSD2U1GIiMSg4mLn03kbeeazxSzelEPX5vV46cq+nNG9BVWqVExB7KeiEBGJIe7O5AWbeXrSYuZv2EXHpnV49opjObdXywoviP1UFCIiMcDdmbpkK09PWkzGmp20b1ybpy9L4/y0VlRNCvYjbyoKEZGAfbVsK0//czEzV+2gdYNaPHZJLy7u24ZqARfEfglVFPochYjEk29Xbufpfy5m+vJttKhfkwcu7Mnl6W2pXjU2CmI/fY5CRKSCzVmzk6f+uYgvl2ylSd0a3HRKJ64Y0I6a1ZICzaXPUYiIBCxrXTbPTFrM5IWbaVSnOr89uxs/PS6FWtWDLYjDUVGIiETZwo27eGbSYj6dt4nkWtW448yu/GxQCnVrxMev4PhIKSISh5ZuzuGPny1m4twN1K1elVtO7cw1J3agfs1qQUc7IioKEZFytnLrHp6dvIRxc9ZRs1oSv/xhJ647sSMNalcPOtpRUVGIiJSTtTv28tzkpYyevZZqSca1J3bkhpM60rhujaCjlYmKQkSkjDZk5/Lnz5cycuYazIyfHteeX57SiWb1agYdrVyoKEREjtLm3Xm88MUy3pmxGnfn8v5tuemUVFom1wo6WrlKqKLQB+5EpCJsy9nHy1OX8+b0lRQUOUP7tuH/BqfStlHtoKNFRUIVhbtPACakp6dfF3QWEUk8O/fm88rU5bzx1UryCoq48NjW/GpwZ1Ka1Ak6WlQlVFGIiETDrrwCXvtyBa//ewU5+YWc27sVt5zamdRmdYOOViFUFCIih5Czr5A3pq3glanL2ZVXyJAeLfj16Z3p1qJ+0NEqlIpCROQAuflFvDl9JS9PXc72Pfmc2q0Zt57ehZ6tk4OOFggVhYhISF5BEe98s5oXpixja84+TurSlFtP68yx7RoGHS1QKgoRqfTyC4t5b+Yanv98KRt35XF8x8a8eGVf+qc0CjpaTFBRiEilVVzsjJ69lj99toR1O3NJb9+Qpy9PY1CnJkFHiykqChGplJZtyeE3YzL5duUO0to24JGLe3Fi5yaYBbMudSxTUYhIpVJQVMxfvlzOHz9bQq1qSTx5aRqX9G2tgghDRSEilUbWumzuGpPJvPW7OLtXC/5wfo+EmccUTQlVFBrhISIHk1dQxLOTl/Dy1OU0rF2dl67sy5CeLYOOFTcSqig0wkNEDvTtyu3cNSaT5Vv2cGm/NtxzTneSa8fXwkFBS6iiEBHZL2dfIY9/spA3p6+iTcNa/P2aAZzYuWnQseKSikJEEs6URZu5+/0s1mfncvUJKdx+RlfqxMn61LFIf3MikjB27MnngYnzGTt7HanN6jL6xkH0a1+5P1VdHlQUIhL33J2P5m7k3vFZ7NxbwK8Gp3LT4FRqVE0KOlpCUFGISFzbvCuPe8Zl8c/5m+jVOpk3hw2ke6vKNd012lQUIhKX3J1RM9fywMT55BcWM+Ksblzzgw5UTaoSdLSEo6IQkbizetteRryfybSl2xjYoRGPXtKbDgm+ylyQVBQiEjeKip03vlrJk58uIqmK8dBFPbmifzuqVNH4jWhSUYhIXFi8aTd3js5kzpqdDO7WjIcu6knL5FpBx6oUVBQiEtPyC4t5ccoy/vzFEurVrMafftSH89NaaYhfBVJRiEjMylizkztHZ7Jo024u6NOK35/bncZ1awQdq9KJuCjMrH7p27v79qgkEpFKLze/iKcnLeK1f6+gWb2avHpVOqd1bx50rErrsEVhZjcA9wF5gIc2O9AxirlEpJL6atlWRoydy6pte/nxwHb85qxu1K+pIX5BimSP4nagp7tvjXaYstKYcZH4tSuvgEc+Wsi7M1bTvnFt3r3uOI7v1DjoWEJkRbEM2BvtIOVBY8ZF4tNn8zdx97i5bNm9jxtO6sivT+tCreoavxErIimKEcBXZvYNsG//Rnf/VdRSiUilsC1nH3+YMJ8JGevp1qIef7kqnd5tGgQdSw4QSVG8DHwOzAWKoxtHRCoDd+eDOeu5b8I8cvYVMvz0Ltx4cieqV9X4jVgUSVFUc/fhUU8iIpXC+p253DMui88XbubYdg14/JLedG5eL+hYEkYkRfGxmV0PTOD7bz3p9FgRiVhxsfPOjNU8+vFCioqd35/bnZ8NSiFJ4zdiXiRFcUXo64hS23R6rIhEbMXWPdw1JpMZK7bzg9QmPHJxL9o2qh10LIlQ2KIwsyrAb9z9vQrKIyIJpLComFf/vYJnJi2mRtUqPD60N5f2a6PxG3EmbFG4e7GZ3QGoKETkiMxfv4s7x2SQtW4XZ/ZozgMX9KRZ/ZpBx5KjEMlbT5+Z2e2UlMWe/Rt1jEJEDmZfYRF//nwpL05ZRoPa1XjhJ305q2cL7UXEsUiK4vLQ15tKbdMxChH5H7NWbefO0Zks27KHS/q24XfnHkOD2tWDjiVldNiicPcOFRFEROLXnn2FPPHpIv42fSWtkmvxt2EDOLlL06BjSTmJZChgNeAXwEmhTVOAl929IIq5RCROTF28hRFj57I+O5erjmvPHUO6UbeGVjBIJJH8a74IVANeCF3/aWjbtdEKJSKxb+fefB6cuIDRs9bSsWkdRt1wPOkpjYKOJVEQSVH0d/e0Utc/N7OMaAUSkdj3SdZG7hmXxY69+dx0SiduHtyZmtU0xC9RRVIURWbWyd2XAZhZR6AourFEJBbt3JvPvePn8cGc9fRoVZ+/DetPj1bJQceSKIukKO4AvjCz5YAB7YGro5pKRGLOFws3c9eYTLbvyefW07rwy1M6US1JQ/wqg0jOeppsZp2BrqFNi9x9X7j7iEji2J1XwIMfLuC9mWvo2rwer/+8Pz1bay+iMon01IR+QEro9n3MDHd/M2qpRCQmTFu6lTtHZ7IhO5df/rATt5zWmRpVdSyisonk9Ni/A52AOfz32IQDKgqRBLU3v5BHP17Im9NX0bFJHUb/YhB92zUMOpYEJJI9inSgu7t7tMOUldbMFim7mSu3c9uoDFZt28uwEzpwx5ldtSxpJRfJkagsoEW0g5QHd5/g7tcnJ+v9U5EjlVdQxMMfLeDSl6dT7M4/rj+O35/XXSUhEe1RNAHmm9kMvr9w0flRSyUiFSpz7U6Gj8xg6eYcfjywHb89+xh9ulr+I5JXwh+iHUJEgpFfWMyfP1/C81OW0bRuDc1okoOK5PTYf1VEEBGpWAs37mL4exnM37CLi/u25t7zepBcq1rQsSQGad9SpJIpLCrm5anL+eNni0muVY1XftqPM3rExWFICYiKQqQSWbo5h9tGZZCxZifn9GrJAxf2pFEdrRch4akoRCqB4mLn9WkreOLTRdSqnsRzVxzLeWmtgo4lcSKSD9ydQMkB7fah2xvg7q4V7kTiwOpte7l9dAYzVmzntGOa8fDFvWhWT2tXS+Qi2aN4DbgVmIWmxorEDXfn7W9W8/BHC0gy48lL07ikb2utXS1HLJKiyHb3j6OeRETKzfqdudw1JpMvl2zlxM5NeOyS3rRqUCvoWBKnIimKL8zsCWAs3//A3eyopRKRo+LujJm9jvvGz6PInQcv7MlPBrbTXoSUSSRFMTD0Nb3UNgcGl38cETlam3fn8duxc/lswWYGpDTiiUt7075xnaBjSQKI5AN3p1REEBE5eh9mrueecVnszS/innOOYdgJHahSRXsRUj4iOj3WzM4BegD/OVXC3e+PVigRicz2Pfn87oMsJmZuIK1tA566NI3UZnWDjiUJJpLTY18CagOnAK8CQ4EZUc4lIocxaf4mRoydS3ZuPnec2ZUbTupIVS1NKlEQyR7FIHfvbWaZ7n6fmT0F6CwokYBk5xZw/4T5jJm9lmNa1ufNYQPo3qp+0LEkgUVSFLmhr3vNrBWwDWgZvUgicihTF2/hrjGZbN69j5sHp3Lz4M5Ur6q9CImuSIriQzNrADwBzKbkjKdXo5pKRL5nz75CHv5oAW9/s5rUZnUZe2U/0to2CDqWVBKRnPX0QOjiGDP7EKjp7tnRjSUi+329fBt3jM5g7Y5crj+pI8NP70LNalp1TipOJAezawO3Ae3c/Toza2dmJ7r7h9GPJ1J55RUU8fgni/jrVyto16g2I284nv4pjYKOJZVQJG89/ZWSOU/Hh66vA0YBKgqRKPlu9Q5uG5XB8i17uOr49vzmrG7Urq5hzxKMSF55ndz9cjO7AsDd95rmAYhExb7CIv702RJe+tcyWtSvydvXDuSE1CZBx5JKLpKiyDezWpQcxMbMOlFq5pOIlI+sddncPiqDhRt3c1l6G+45tzv1a2ppUgleJEVxL/AJ0NbM3gZOAH4ezVAilUlBUTEvTlnGs5OX0LBOdV7/eTqDuzUPOpbIf0Ry1tMkM5sNHEfJokW3uPvWqCcTqQSWbNrN8JEZzF2Xzflprbjv/B401NKkEmMOWRRm1veATRtCX9uZWTuNGRc5ekXFzqtfLuepSYupW6MqL/ykL2f30udYJTaF26N4KvS1JiUjxjMo2aPoDczkv2dBicgRWLl1D7eNymDWqh2c0b05D13Ui6b1agQdS+SQDlkU+8eLm9lYoK+7zw1d70nJGtoicgSKi52/f72KRz9eSNUk45nL07iwj5YmldgXycHsrvtLAsDds8zsmChmEkk463bmcufoDKYt3cbJXZry2CW9aZFc8/B3FIkBkRRFppm9CrwVuv4TIDN6kUQSx4FLkz58US+uGNBWexESVyIpiquBXwC3hK5PBV6MWqIDhPZebgGaAJPdvcKeW6QstubsY8TYuUyav4kBKY148tI02jWuHXQskSMWyemxecAzoT9HxMxeB84FNrt7z1LbhwB/ApKAV9390TDPvwC40cyqAG9SgSUlcrQ+ydrI3e/PZXdeIXeffQzDftCBJC1NKnEq3OmxI939MjObS+hT2aW5e+8IHv8N4M+U/ILf/7hJwPPA6cBa4FszG09JaTxywP2HuftmMzufkr2av0fwnCKByc4t4L7x8xj73Tp6tq7Pu5f1oUvzekHHEimTcHsU+99qOvdoH9zdp5pZygGbBwBL3X05gJn9A7jA3R851HO5+3hgvJlNBN452jwi0fTlki3cObpkUaFfndqZmwenUk1Lk0oCCHd67IbQ11Xl/JytgTWlrq8FBh7qxmb2Q+BioAbwUZjbXQ9cD9CuXbvyyCkSkb35hTz68ULenL6KTk3rMPYXg7SokCSUcG897abkLSfj+289GeDuXiGL9Lr7FGBKBLd7BXgFID09/X/eKhOJhtmrd3DbyAxWbN3DsBM6cOeQrlpUSBJOuD2KaL2xug5oW+p6m9A2kbiRX1jMHz9bzEv/WkbL5Fq8c91ABnXSOHBJTId9A9XMrjnItkOepRSBb4HOZtbBzKoDPwLGl+HxRCrUgg27uOD5abwwZRlD+7Xhk1+fqJKQhBbJ5yguMbM8d38bwMyeB2pF8uBm9i7wQ6CJma0F7nX318zs/4BPKTnT6XV3n3dU6UUqUFGx88rU5Tw9aRHJtarx6lXpnNZd48Al8UVUFJSccVQMDAF2uvuwSB7c3a84xPaPCHNg+miZ2XnAeampqeX90FLJlR7kd1bPFjx4YU8a19UgP6kczP3gx33NrPQq7vWAccA04PcA7r496umOUnp6us+cOTPoGJIA3J23vl7Fwx8tpFqScf8FPbmgTyuN4JCEZGaz3D39wO3h9ihm8f2zngw4J/THgY5RyCkSMzZk53Ln6Ey+XLKVEzs34fGhvWmZHNG7riIJJdxZTx0qMohIrHB3Ppiznt9/kEVBkfPAhT25cmA77UVIpXXYYxRmdtXBtrv7mwfbLhLPtuXs455xWXyctZF+7Rvy1KVppDSpE3QskUBFcjC7f6nLNYFTgdmUmt8kkggmzd/EiLGZ7Mot5K4h3bj+pI4a5CdCZNNjby593cwaAP+IWqIy0FlPcjR25xVw/4T5jJq1lmNa1ueta9Po1qJCBg+IxIVI9igOtAeIyeMX7j4BmJCenn5d0FkkPny1bCt3jMpkQ3YuN53SiVtO7UL1qhrkJ1JaJMcoJvDfWU9VgO7AqGiGEom2vIIiHvtkIX+dtpIOTeow6sZB9GvfMOhYIjEpkj2KJ0tdLgRWufvaKOURibo5a3YyfOQclm/Zw8+Ob89dZ3WjdvWj2bkWqRzCTY/t5u4L3f1fZlbD3feV+t5x7v51xUQUKR8FRcU8N3kJz09ZRrN6NXjrmoH8oLNmNIkcTrj/jXoH6Bu6PL3UZYAXDrguEtMWb9rN8JFzyFq3i4v7tube83qQXKta0LFE4kK4orBDXD7YdZGYVFTsvPbv5Tz5z8XUq1GVl67sx5CeLYKOJRJXwhWFH+Lywa7HBJ0eK6Wt3raX20dlMGPlds7o3pyHL+5FEw3yEzli4YqijZk9S8new/7LhK63jnqyo6DTYwVKRnC8O2MND06cT5IZT12axsV9W2sEh8hRClcUd5S6fOAoVo1mlZi0aVced43JZMqiLZyQ2pjHh6bRuoEG+YmURbihgH+ryCAiZTU+Yz2/G5fFvsIi/nBed646PoUqGsEhUmY6eVzi3o49+dzzQRYTMzfQp20Dnr4sjY5N6wYdSyRhqCgkrn2xcDN3jslk59587jizKzec1JGqSRrBIVKeVBQSl3L2FfLQxPm8O2MNXZvX442r+9OjVXLQsUQSUiSznpoC1wEppW8f6brZIuXtm+XbuH10Bmt35HLDyR0ZfnoXalRNCjqWSMKKZI/iA+BL4DOgKLpxykafo0hseQVFPPnpIl6btoJ2jWoz6objSU9pdPg7ikiZmHv4z86Z2Rx371NBecpFenq6z5ypM3gTydy12QwfOYclm3O48rh2jDjrGOrU0DunIuXJzGa5e/qB2yP5L+1DMzvb3T+KQi6RsAqKinnhi2U89/kSGtetzt+GDeDkLk2DjiVSqURSFLcAvzWzfUABJZ/MdnfXEmASVUs35zB85Bwy12ZzQZ9W3H9+T5Jra5CfSEWLZCnUehURRGS/4mLnja9W8tgnC6ldPYnnf9yXc3q3DDqWSKV12PUozOyg48TdfXb0YklltW5nLrePzGD68m2c2q0Zj1zSi2b1agYdS6RSC7dHMRy4HnjqIN9zYHBUEkml5O6Mmb2O+8bPo9idxy7pxWXpbTXITyQGhJv1dH3o6ykVF0cqo605+xgxdi6T5m9iQEojnrosjbaNagcdS0RCdH6hBOqTrI3c/f5cducVcvfZxzDsBx1I0iA/kZiSUEWhD9zFj115Bdw3fj5jZq+lR6v6vHt9H7o013kTIrEo7PQ0K9G2osKUlbtPcPfrk5M18yeWfbV0K0Oemcq4Oev41eBU3v/lCSoJkRgWdo/C3d3MPgJ6VVAeSWC5+UU89slC3vhqJR2b1GHMLwbRp22DoGOJyGFE8tbTbDPr7+7fRj2NJKw5a3YyfOQclm/Zw88HpXDXkG7Uqq5BfiLxIJKiGAj8xMxWAXv47yeze0c1mSSEgqJinpu8hOenLKNZvRq8fe1ATkhtEnQsETkCkRTFmVFPIQlp8abdDB85h6x1u7i4b2vuPa8HybU0gkMk3kQywmOVmaUBJ4Y2fenuGdGNJfGsuNh5fdoKHv90EXVrVOWlK/sxpGeLoGOJyFGKZOGiWyhZuGhsaNNbZvaKuz8X1WQSl9Zs38vtozL4ZsV2Tu/enIcv6kXTejWCjiUiZRDJW0/XAAPdfQ+AmT0GTAdUFPIf7s7ImWu4f8J8zIwnhvZmaL82GsEhkgAiKQrj+yvbFYW2iQCweXceI8bMZfLCzRzXsRFPXppGm4YawSGSKCIpir8C35jZ+6HrFwKvRS+SxJOP5m7g7vfnsje/iN+d252rB6VQRSM4RBJK2KIwsyrA18AU4AehzVe7+3dRznVUNMKj4mTvLeDe8VmMm7Oe3m2SefqyNFKb6dPVIokokjWzv3P3YysoT7nQmtnR9eWSLdwxKpMtOfu4eXAqN52SSrWksNNgRCQOlGXN7Mlmdgkw1g/XKpLQ9uYX8ujHC3lz+ipSm9Xllav60buNRnCIJLpIiuIGShYxKjSzPLRmdqU0a9UObhs5h5Xb9nLNDzpwx5ldqVlNIzhEKoNIjlEMcfdpFZRHYkx+YTF/mryYF6cso2VyLd65biCDOmkEh0hlcrjpscVm9mcgro5RSPlYuHEXt76XwYINu7i0Xxt+f1536tXUCA6RykbHKOR/FBU7r365nKf+uZj6taryl6vSOb1786BjiUhAjuQYRZGZ5aJjFAlt1bY93D4qg29X7mBIjxY8dFFPGtfVCA6RyiySoYA6Ob4ScHfenbGGByfOJ6mK8czlaVzYp7VGcIhIREMBDfgJ0MHdHwgtjdrS3WdEPZ1UiE278rhrTCZTFm3hhNTGPDE0jVYNagUdS0RiRCRvPb0AFAODgQeAHOB5oH8Uc0kFmZCxnnvGZbGvsIj7zu/BT49rrxEcIvI9Ea1w5+59zew7AHffYWbVo5xLomzn3nx+98E8JmSsp0/bBjx9WRodm9YNOpaIxKBIiqLAzJIABzCzppTsYUicmrJoM3eOzmT7nnxuP6MLN57ciaoawSEihxBJUTwLvA80M7OHgKHAPVFNJVGxZ18hD320gHe+WU2X5nV5/ef96dk6OehYIhLjIjnr6W0zmwWcSsmpsRe6+4KoJzsKmh57aDNXbmf4yAzW7NjL9Sd1ZPjpXTSCQ0QictjpsfFI02P/a19hEc9MWsLLU5fRpmEtnhyaxsCOjYOOJSIxqCzTYyVOzV+/i+Ej57Bw426uGNCWu8/pTt0a+icXkSOj3xoJqLComJenLuePny2mQe3qvP7zdAZ30wgOETk6KooEs2LrHm4bOYfZq3dyTq+WPHhhTxrW0dnMInL0VBQJwt1565vVPDxxAdWSjD/9qA/np7XSCA4RKTMVRQLYmJ3HnWMymbp4Cyd2bsITQ9NokVwz6FgikiBUFHHM3RmfsZ7fjcuioMh54MKeXDmwnfYiRKRcqSji1I49+dzzQRYTMzfQt10DnrqsDx2a1Ak6logkIBVFHPp84SbuGjOXnXvzuXNIV244qRNJGuQnIlGioogjOfsKeWjifN6dsYZuLerxt6sH0L2V1o8SkehSUcSJGSu2c9uoOazdkcuNJ3fi1tM7U6OqRnCISPSpKGJcXkERz0xazCtfLqdtw9qMvOF4+qc0CjqWiFQiKooYNm99NsPfy2DRpt38eGA77j77GOpoBIeIVDD91olBpUdwNKxdnb9e3Z9TujYLOpaIVFIqihizYuseho+cw3erd3Ju75Y8cIFGcIhIsFQUMaK42Hnrm1U8/NECalRN4tkrjuX8tFZBxxIRUVHEgg3Zudw5OpMvl2zl5C5NeXxob5rX1wgOEYkNKooAuTsfzFnP7z7IorDIeeiinvx4gEZwiEhsUVEEZPuefO4ZN5eP5m6kX/uGPHVpGikawSEiMUhFEYDJC0pGcGTn5nPXkG5cf1JHjeAQkZiVUEVhZucB56WmpgYd5aB25xXw4IcLeG9myQiOv18zgGNaagSHiMS2hCoKd58ATEhPT78u6CwH+mb5Nm4blcH6nbn84oed+PVpGsEhIvEhoYoiFuUVFPHUPxfx6r9X0K5RbUbdeDz92msEh4jEDxVFFGWty+bW9+awZHMOVx7XjhFnaQSHiMQf/daKgsKiYl6csow/TV5C47rVeePq/vxQIzhEJE6pKMrZsi05DB+ZQcaanZyf1or7L+hBg9oawSEi8UtFUU6Ki503p6/k0U8WUrNaEs9dcSznaQSHiCQAFUU5WL8zlztGZzBt6TZ+2LUpj12iERwikjhUFGXg7rz/3TruHT+PomLn4Yt6ccWAthrBISIJRUVxlLbl7OPu97P4ZN5G+qc05MlL02jfWCM4RCTxqCiOwqT5mxgxNpNduYWMOKsb156oERwikrhUFEdgd14BD3w4n5Ez13JMy/q8dW0a3VpoBIeIJDYVRYSmL9vG7aMy2JCdy02ndOKWU7tQvWqVoGOJiESdiuIw8gqKeOLTRbz27xWkNK7NqBsH0a99w6BjiYhUGBVFGHPXZnPryDks3ZzDT49rz4izu1G7uv7KRKRy0W+9gygoKuaFL5bx3OdLaFK3Bm8OG8BJXZoGHUtEJBAqigMs3ZzDbSPnkLE2mwv7tOK+83uSXLta0LFERAKjoghxd974aiWPfryQWtWTeP7HfTmnd8ugY4mIBE5FEWJmzF+/ixNSm/Doxb1ophEcIiKAiuJ7HryoJ9WTqmgEh4hIKSqKUrQ0qYjI/9InxkREJCwVhYiIhKWiEBGRsFQUIiISlopCRETCUlGIiEhYKgoREQnL3D3oDOXOzLYAq8r4MMlAdgXf/0jvE8nty+M2TYCtR8Z+hmcAAATCSURBVJArlpX13zVWnrM8HjPar9FIb1vW16hen+X3vO3d/X8noLq7/hzkD/BKRd//SO8Tye3L4zbAzKD/PWLl3zVWnrM8HjPar9FIb1vW16hen9F/Xr31dGgTArj/kd4nktuX120SRRA/azSeszweM9qv0Uhvq9fofwX1c4Z93oR860nKl5nNdPf0oHOIHIxen9GnPQqJxCtBBxAJQ6/PKNMehYiIhKU9ChERCUtFISIiYakoREQkLBWFlImZXWhmfzGz98zsjKDziJRmZh3N7DUzGx10lnimoqjEzOx1M9tsZlkHbB9iZovMbKmZ/SbcY7j7OHe/DrgRuDyaeaVyKafX53J3vya6SROfznqqxMzsJCAHeNPde4a2JQGLgdOBtcC3wBVAEvDIAQ8xzN03h+73FPC2u8+uoPiS4Mr59Tna3YdWVPZEozWzKzF3n2pmKQdsHgAsdfflAGb2D+ACd38EOPfAxzAzAx4FPlZJSHkqj9enlA+99SQHag2sKXV9bWjbodwMnAYMNbMboxlMhCN8fZpZYzN7CTjWzEZEO1yi0h6FlIm7Pws8G3QOkYNx922UHD+TMtAehRxoHdC21PU2oW0isUCvzwCoKORA3wKdzayDmVUHfgSMDziTyH56fQZARVGJmdm7wHSgq5mtNbNr3L0Q+D/gU2ABMNLd5wWZUyonvT5jh06PFRGRsLRHISIiYakoREQkLBWFiIiEpaIQEZGwVBQiIhKWikJERMJSUYgcgpn92sxqR/HxW5XHOglmllMeeUQORZ+jEDkEM1sJpLv71qCzhGNmOe5eN+gckri0RyGVnpnVMbOJZpZhZllmdrmZ/QpoBXxhZl+EbneGmU03s9lmNsrM6oa2rzSzx81srpnNMLPUgzzHyWY2J/TnOzOrZ2Yp+xflMbNXS31/i5ndG9p+h5l9a2aZZnZfmJ/hoVD+r82seTT+nqTyUlGIwBBgvbunhRbI+SQ0FXc9cIq7n2JmTYB7gNPcvS8wExhe6jGy3b0X8Gfgjwd5jtuBm9y9D3AikFv6m+5+beh7FwBbgTdCS8t2pmQNhj5Av9BiPgeqA3zt7mnAVOC6o/trEDk4FYUIzAVON7PHzOxEd88+yG2OA7oD08xsDvAzoH2p779b6uvxB7n/NODp0J5Kg9DMou8xs5rAKOBmd18FnBH68x0wG+hGSXEcKB/4MHR5FpAS5mcVOWJaj0IqPXdfbGZ9gbOBB81ssrvff8DNDJjk7lcc6mEOcXn/czxqZhNDzzHNzM4E8g642UvAWHf/rNRzPuLuLx/mRyjw/x5sLEL/XUs50x6FVHpm1grY6+5vAU8AfUPf2g3UC13+Gjhh//GH0HGNLqUe5vJSX6cf5Dk6uftcd3+MklHZ3Q74/k1APXd/tNTmT4FhpY6FtDazZmX4UUWOiv7PQwR6AU+YWTFQAPwitP0V4BMzWx86TvFz4F0zqxH6/j3A4tDlhmaWCewDDrbX8WszOwUoBuYBHwMtS33/dqAg9LYWwEvu/pKZHQNML1manBzgSmBzmX9ikSOg02NFyiheTqMVOVp660lERMLSHoWIiISlPQoREQlLRSEiImGpKEREJCwVhYiIhKWiEBGRsFQUIiIS1v8DJ6Rlp/j1r6IAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] @@ -456,17 +451,13 @@ "plt.ylabel('error in Euklidean norm')\n", "plt.xlabel('step size h')" ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Obviously the algorithm is not consistent, although it is supposed to be of order 2.\n", - "Even after several hours of debugging, I couldn't figure out the problem." - ] } ], "metadata": { + "@webio": { + "lastCommId": null, + "lastKernelId": null + }, "kernelspec": { "display_name": "Python 3", "language": "python", @@ -482,7 +473,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.7.8" + "version": "3.8.5" } }, "nbformat": 4, -- GitLab