diff --git a/Exercise sheet 6/exercise_sheet_06.ipynb b/Exercise sheet 6/exercise_sheet_06.ipynb index abd31a8..8af69d0 100644 --- a/Exercise sheet 6/exercise_sheet_06.ipynb +++ b/Exercise sheet 6/exercise_sheet_06.ipynb @@ -552,7 +552,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEtCAYAAAAbeVcBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAxTklEQVR4nO3debxVdb3/8ddbRAUnMFARUEzQHME8WV7tOqBpdkut7DqU1rUoh6vZbdBG60ZZ2mC3sEi96s0hf5plpaWZ6HUkcEbtioICogyCE4mgn98f3++WxWaffdY5nH32Gd7Px2M/zprXZ+2z1v6s7xq+X0UEZmZmbVmn2QGYmVnP4IRhZmalOGGYmVkpThhmZlaKE4aZmZXihGFmZqU4YZiZWSlOGGZmVooThq0VSbMlHdjsOMqSNEPSfp28zB0k3SfpJUmnduayC+voUd9zbybpYknfLjltr/q/OWE0iKQpkpZIWr/ZsXR3XXlQRcTOETGlkxf7RWBKRGwcET9Z24X1th+ZMiRtlLf7mMKwjSU9LenDki6TdFHVPPtKWixpWNdH3Dc5YTSApFHAu4EAPtCA5a/b2cvsqbrJd7ENMKMjM3aT+OtSsnuN4WMl9euMdUTEy8AE4DxJQ/Pg7wPTIuJq4FTgUEkH5XVvAPwS+I+ImN8ZMVgJEeFPJ3+ArwN3AD8E/lAYfgZwddW05wE/yd1bAdcAC4FZwKmF6WYDXwIeBJYD6+blPQG8BDwCHFGY/u3AfXnc/wN+DXy7ML7VddXYnpHAb/K0i4GfVsV1YO4OYHRh3MVV6/wSMC/H9HdgPPA/wBvAP4CXgS925LuoEfMa66oR77/mdVY+y0klhdLfD/BX4HXg1byM7YEdgSnAUlIi+UDVPK3GX+f7mA18Ps/zQv5/btDB/+fLwMjc/Zn8f9si938BuKBq+m3zcg8pDHs3sAjYZW2XXzXvxcAVwH6kfW1YYdyReds2BL4L3FBnObPzuh4EXgEuBLYAbsj7xF+AwYXpW/2fAbsD9+b5fg1cScljicL+1hs+TQ+gN36AmcBJwB7AisLBsg2wDNgk9/cD5gPvIpX2ppOSzXrAW4EngYPztLOB+0k/3gPysCPzzroO6cfvFWBYnv8p4DSgP/BB4LXKTt7Wuqq2pR/wAPCjfKBuAOxTGP/mAUGdhAHsAMwBtsr9o4DtqpdRJr5a30VVzKXXVZhnE+BR4NPt+X7yvFOAT+bu/vn//+U87wGkH5odqr6zevGvEWMeNjX/vzfLsX6mvf/PPP0cYCdAwEPA48Dbcv9MYGyNef6Z9AO+P7AnVQlkbZdfmHcw6ZhYBHyixvirgetyLFvXWc5s4G5SkhgOLCD96O8OrE9K9N9o63/GqmPp9Dzdh0nHdKljqbX9rad+mh5Ab/sA++Qdakjufww4vTD+duC43H0Q8ETufifwdNWyzgT+O3fPBv6tjXXfDxyWD+55gKrW++0y66oavlf+cVjjLL4QV5mEMToftAcC/VtbRmd8F+1ZVx62DvAH4Pz2fj953BRWJYx3A88C6xTGXwGcVRVDvfhrxTgb+Gih//vAzzsY70OkH/2DgWvz513Ae4Fb68R1MOmHegGF0mxnLb8w/19IJ1ab1hi3BakEc1oby5gNHFvov6by/839/w78tq3/GelYeobVj6U7KXks1fpf9uSP72F0vuOBGyNiUe6/PA+j0H907j4m90MqfWwlaWnlQzrj2aIw75ziiiQdJ+n+wvS7AENIZ6HzIu+xNeYts66KkcBTEbGy7U1vXUTMBD5LOggXSLpS0latTN7u72It1gUwEdiYdJ287PpbsxUwJyLeKAx7inSWW9Rq/HU8W+heBmzUwXiX5Hk/S7ok+iLpzP4koN5N+6eBlaSSwuw603V0+Uj6KKlE+Bfge9XjI+I5UumjzD2j5wrd/6jRX/n+6v3Pah1LTxW612Zf6XG6/Q23nkTSAOAjQD9JlYN7fWCQpLER8QDpfsIPJI0AjiCdwUP6AZkVEWPqrOLNnVbSNqSbfuOBuyLidUn3kw7m+cBwSSrs6CNJ9zvKrqtiDrC1pHVLJI1lwMBC/5bA3DeDj7gcuFzSJsAvSD8IHytuVzviq55n9ZGtr2s1ko4iJfB3RMSKdqy/Nc8AIyWtU/gB2hr4v3bEX3fbamhvvEuBFtL9gSmSPgSMJZ1w/LbWDJK2A24i3Xt5Cbhe0oERUeuHu93Lz+vYnHTp8yOkkvkMSZdHxG0lt6uj6v3Pah1LW9OxY6nHcwmjcx1OugG6EzAuf3YE/hc4DiAiFpIuYfw3aUd7NM87FXhR0pckDZDUT9Iukt7Ryro2JP2wLASQ9AnSAQlwV47jFEnrSjqMdImgoj3rmko6aM6WtKGkDSTt3UpM9wPH5OUdAuxbGZHfVTggP2b8KukM7/U8+jnStd+OxLeGNtZVnG534L+Aw/P/pTPWfw/pXtIXJfVXeufj/aQbpWVVfx9taW+8S0jX5Ctn+y+S7nf9IiJqfU9bATcDEyPi4oi4hnQD/kZJteJs1/ILfkq6THRLpCefvgj8sgseTa/3P7uLVKo6NR9LH6Tjx1KP54TRuY4nXbt8OiKerXxIB8KxhUcoLyddX69cjiIfSO8nJZlZpGL3BcCmtVYUEY8APyDt0M8Bu5KezCIiXiPd6D6BdLb3UdI1+uXtXVdh2tGkSxJzSTfYazktT7sUOJbVzybXB87O63oW2JxUdIf0xMtXc5H+8+39Lmqot66iw0iXSm6X9HL+3LA268/f/QdI1+sXAZNI96weKxk7VH0fJdbZ3niXkK4uVPa/F4FBpBJrLYtJj6+eX1jnZaRLTAs6YflIOpx0/+8LhXVcQNrfvt7afJ2h3v+scCx9nLRd/0p6YrAy79ruqz2KVr80Z72VpHtIN0n/u9mxmFnP5BJGL6X0FuyWuRh9PLAb8Kdmx2VmPZdvevdeOwBXkZ4EeQL4cPiNWDNbC74kZWZmpfiSlJmZleKEYWZmpThhmJlZKT0mYUj6rqTPNjuO3kbS1vn9g365f4qkT+buYyXdWJg2JI3uqljMrPEkfUBSqRdLe0TCUKof/zhSFQ/NWP8XJD2s1KLaLElfqBo/W9I/Ci9/3djastpYz1WSDpK0fqFqkYbKLxluVOsN3Ii4LCLe06h1q6qhoHqxdFeSNpO0UNLtVcPHSZouaVn+O65q/OmSnpX0gqSLOvI2s6SdJN2Ru7+lBrX21xcVT5y6cJ2nSJomabmki2uMHy/psbxP3aJUPVBlnCR9T6lBqcWSvi9JhfGj8jzL8jKKx911wC6Sdmsrxh6RMEhvWV4fEf9o0vpFSliDgUNIVW4cVTXN+/OP3UZr8SO7B6kK5t2AhzscbTegHtAwUCf5Hqmq8TdJWg/4HfAr0j5zCfC7PBxJB5PaMhlPqmjvrcA3O7DuPUhVa1e67+3AMroNlyx5Bvg2cFH1CElDSG+Yf41Uvf00UtscFRNIVRONJf1+/Aupqv6KK0jt47wF+ApwtVY1VFUZP6HNCJtdXW6ZD6nu+mLVzrcCH8rd+5DqVDo09x8I3N/geH4C/FehfzZrWYUx6Yflydx9IvD9qvGbkhqBmU+quvzbQL88rh9wLqlagieBk/N3sm6t+Ei1uP4qd4+qmnYKq6rq/jhwe2G+INXo+mRe1znkKqHztHeQKo97Pse3Xf7fLc7TXwYMytOv0VBQjVi2IrV98DypvYJPVW3DVcClpMrwZgAthfE1G1Dq5P1gL1LVLJ+o+p7ew5rVyz9Nbj+CVGXGdwrjxgPPdmD95wHH5+5ngI0K4/YjVavxRVL1HfNJPyiHkirVex74cmH6dVjVINfi/N1uVhj//0jVrLwA3AbsXBh3KKkBr5fydn++1v5T2IdG5+6LgfOB60l1OR1I/caIzspx/Cqv6yFSg1Vn5m2cA7yn5DHzcVKV/+eSqvyYBbw3j5vI6o1i/ZR00vijvJ4XSA0zrdF4VCftV98GLq4aNgG4s9C/IenYeVvuvxOYUBh/AnB37t6eVC3QxoXx/0tuTyX3702q265ubD2lhLEr6aCvuJV0QECqr/5JVlV09895/BokHaNCNcQ1Plu3FUgu5r2bNatXvixfmrhR0tiyG5aLmUtJO/uI3H0ecHKOqbJdl5AqQRtNagTmPUClyPwp0hnF7qRaQj9cdv0dcERex9tJdTH9W2HcO0n/i81JB51I9SJtRaqEcSTpoCciPkb6Ea2UzL5fY11XkH70tiJt03ckjS+M/wCpgrhBpMTyU0iVDwKnkGqg3ZjULsPsWhsj6Yx6+0RrX0I+G/5ZXk/1y0w7Aw9GPhKzB/PwyvgHCuMeALaQ9JbW1le17ptybCcD/yXpRVJ12nMl3VCYdEtSg1fDSfUx/ZJUr9gepH3461pVeeCppISyL+n7XpK3r+IGYAzpf3svKflXXAh8On/Xu5BOEso6hlXVy98J/J70fQwnJdLP5hJZxftJJxuDSWfMfyYlu+HAt1j9snW9YwbS/vp3UpMA3wculKSI+ArpB/WUvG+ekuf9Z9KP7yBSnVKLa22QpEl19qkH2/HdFK22z0TEK6TkXm+fKo57MiJeamU8pFLyKKXanVvXiAzZgIy7gpxJY9UZ2YO5+0+knaCSTW8FPtjAWL6Zv+z1q7LzAFLV3meSzsQGtXO5vyL9AG9IOgMcWBi3BekMYUBh2NHALbn7r6x+tvAeGlfCKDbTeRJwc2Hap9vYxsOB+wr91XG9GQspubzO6mdF3yWfeeVt+Eth3E7AP3J3qw0odeJ+cDqrGlyq/p6+BlxZNf1l5EaUSAd68Xvsn7d7VDvWvz2pvWtIFSt+oWr8fqQz0MoZ9cZ5He8sTDOdVFMvpB+M8YVxw0jHXa3mbwflZW2a+58mXf7YpGq61b6Xwj5ULGFcWhjXVmNEZwE3Fca9n1QCqN7GQbR9zHwcmFkYNzDPu2X1cZD7DyAdl++i0NBSg/atWiWMC4Gzq4bdAXw8d7/O6r+RY/L2iFSt/91V804srqOwD7baimFEzylhLCHtDBV3AdtL2oJUS+SlpPrsh5CqHm5I/fmSTiHdy3hfRCyvDI+IOyLiHxGxLCK+S6qt9d0llzk3ny0eTTojWkBqlOUZST/Mk21D+ofOL5z5/oJ0tge5AZjCYosNvHS26vVs1co4JG2u1HjRvHwW/CvS2VwZWwHPx+pnRdUNEVU3KLSBUrsd7W1AqV3ysk4lXQuu5WVSk69Fm5Auo9QaX+l+iTbkG6NLyWeIufs/WVW77eaFyRfHqgcIKvf/WmtEaBvg2sL+9SjpR2gLpSq7z5b0RP4/zs7zVP6XHyJdlnpK0q2S9qK89jbsVR3/ohrbuBFtHzNQ2H8iYllh3jVExF9JJdifAc9Jmtzm2Xjn6sg+9XKkbNDWvLDq93VpvSB6SsJ4kHRGBbz5z51Oqk774UhVEN8JfI7U5OmiWgtRekz05TqfVi9JSfo38o3KiJjb2nSVEEmZvU0RMYJ0I/0vETEImAycHBGDIuJzebI5pLOlIXn4oIjYJCIqRcr5pDPyiurteIU1GzbqqOr1PFPcnKppv5uH7RYRm5Auh6jO9EXPAJtJKp4obE26Ft2miLg8IvYh/XAENVpvA5D05Xr7RCuL35N0Bv6I0tNs5wF7Kj311I90uXK34lMqpBuRlcuYM0g3JyvGAs9FRM1LHFXb9dO8n9xKOuvdhtQi3KZ5v6hV3XgZc0jX8AcVPhtExDzSZaPDSCW2TUklQcj/y4j4W0QcRvox/i3p/gdU7XeSau13xX2g0hhRMYaNI+LQDm5PvWOmLWvsmxHxk4jYg3QpZ3sKVbEXSfp5nX2qTEuBtay2z0jakHSPsN4+VRz31qpjqTge0iXj2RHxYr0gekrCuJ5CYzzZraTrx5X7FVOq+tcQ6THRjep8nq41n6Rjge8AB0XEk1Xjtpa0t6T1lBoX+gLpzKvyuOMopfcXRtXZvuITLm8nPQFRjHs+cCOppb5NJK0jabvC/Y2rSA28jJA0mJTYiu4HjlJqHGZt73F8QdJgSSNJCfvXdabdmHR2s1TScNY8wFptKCgi5pBOAr6bv9fdSDfyLqs1fZFKNqCU1/OdevtEK6u4gfSjOS5/vk66nj4un+1Oyes7VekR6VPyfJVr+5cCJyg9FjsY+Crp8kwl/otV47HKKmNJpYy30zlPR/0cmKj8qKakoUoNb0H6Py4nXbMfSDoWKrGul0/ENo3UYuGLrPquK6WgcZI2IN+/qqPTGiMqccy0ZbV9U9I7JL1TUn9SInyV1vepz9TZp1pNWEo1S29AeoilX97vK08bXkt69PVDeZqvky7LV9pZuRT4nKThuQT8H+R9KiL+j/Qb8I28zCNIJzDXFFa/L2m/rqunJIxLgUOVmkCtuJW0I9/WSn9n+jbpcbS/Fc4Ufp7HbUx60mMJ6ez3ENKZWuVscSTpUkq9M+M9gHvzGenbqN1e8XHAeqSnUZYAV5POciHdzPwz6QC9l0IDL9nXSGcjS0j3YC6n435HKt3dD/yRdG21Nd8k/aC9kKetjquthoKOJv0wP0M6YL4RETeViLFsA0odEhHLY/UGsl4AVuRucon3cNL/bCnpwYDD83Ai4k+km6y3kPaNp4BvFFYxknzCUUsuCT+fS9pvZ9WjtWvjPNKDAzdKegm4m3RPAdLxV9mHH8njij4GzM6Xqz5DKklWfqi+RWqf+3HSU0mtis5vjKjeMdOW84APS1oi6SekSzi/zMt5ipQ8z+1gXK35Kunk5gzSd/iPPIxILUJ+iHTvYQnpf1N8tP8XpAcGHiI9kv9HVn8A4CjSwypLSMfGh2P1ViaPpsR7bj2mtlpJ3wEWRMSPmx1Le0j6KrAwIrrspcNcmplFuuHbVjvc1o0ovavxAOky3oq2pjdbW5LeD3wsIj7S5rQ9JWFYeU4YZtYIPeWSlJmZNZlLGGZmVopLGGZmVkqPriBuyJAhMWrUqGaHYWbWo0yfPn1RRAxte8rV9eiEMWrUKKZNm9b2hGZm9iZJHaoNwpekzMysFCcMMzMrxQnDzMxKccIwM7NS+lzCeHjeC5x02XReXdFjmo02M+sW+lzCWLpsBdc/9CyTbpnZ7FDMzHqUPpcw9hkzhCN2H875tz7BzAVttldjZmZZj34Po6O+8r4d+etjCzjuwqmce+RYiu3cjBg8gJGbDawzt5lZ39QnE8aQjdZnp2GbcNeTiznmgntWG3fa+DGcftD2rcxpZtZ39cmEsejl5Twy/0WGDxrAOUfuhli9hGFmZmvqkwlj4h8fZdlrK7nmxL0YvfnGbc9gZmZ976b37Y8v4tr75nHivts5WZiZtUOfSxiDBvbn0F235KT9Rzc7FDOzHqXPXZLaZfimTDp2j2aHYWbW4/S5EoaZmXWME4aZmZXihGFmZqU4YZiZWSkNSxiSNpA0VdIDkmZI+mYefpakeZLuz59DC/OcKWmmpL9LOrhRsZmZWfs18imp5cABEfGypP7A7ZJuyON+FBHnFieWtBNwFLAzsBXwF0nbR4TrITcz6wYaVsKI5OXc2z9/os4shwFXRsTyiJgFzAT2bFR8ZmbWPg29hyGpn6T7gQXATRFRqenvFEkPSrpI0uA8bDgwpzD73DysepkTJE2TNG3hwoWNDN/MzAoamjAi4vWIGAeMAPaUtAtwPrAdMA6YD/wgT65ai6ixzMkR0RIRLUOHDm1I3GZmtqYueUoqIpYCU4BDIuK5nEjeAH7JqstOc4GRhdlGAM90RXxmZta2Rj4lNVTSoNw9ADgQeEzSsMJkRwAP5+7rgKMkrS9pW2AMMLVR8ZmZWfs08impYcAlkvqREtNVEfEHSf8jaRzpctNs4NMAETFD0lXAI8BK4OTe9ITUw/NeYNKUmfzwI+PYoH+/ZodjZtZuDUsYEfEgsHuN4R+rM89EYGKjYmqmpctWcP1DzzJ66Ew+954dmh2OmVm7+U3vLrLPmCEcsftwzr/1CWYueKnZ4ZiZtVufq968mb7yvh3562MLOO7CqZx75Fik1ZuGHbnZwCZGZ2ZWnxNGFxqy0frsNGwT7npyMcdccM9q404bP4bTD9q+SZGZmbXNCaMLLXp5OY/Mf5HhgwZwzpG7IVYvYZiZdWdOGF1o4h8fZdlrK7nmxL3cnriZ9Ti+6d1Fbn98EdfeN48T993OycLMeiQnjC4yaGB/Dt11S07af3SzQzEz6xBfkuoiuwzflEnH7tHsMMzMOswlDDMzK8UJw8zMSnHCMDOzUpwwzMysFCcMMzMrxQnDzMxKccIwM7NSnDDMzKwUJwwzMyvFCaOXe3jeC5x02XReXdFrWrs1syZpWMKQtIGkqZIekDRD0jfz8M0k3STp8fx3cGGeMyXNlPR3SQc3Kra+pNI07KRbZjY7FDPr4RpZwlgOHBARY4FxwCGS3gWcAdwcEWOAm3M/knYCjgJ2Bg4BJknq18D4+gQ3DWtmnaVhCSOSl3Nv//wJ4DDgkjz8EuDw3H0YcGVELI+IWcBMYM9GxdeXfOV9OzJwvXU57sKp3DlzEXc9sfjNz5znlzU7PDPrIRpaW20uIUwHRgM/i4h7JG0REfMBImK+pM3z5MOBuwuzz83Dqpc5AZgAsPXWWzcy/F7DTcOaWWdoaMKIiNeBcZIGAddK2qXO5KoxLGosczIwGaClpWWN8bYmNw1rZp2hS9rDiIilkqaQ7k08J2lYLl0MAxbkyeYCIwuzjQCe6Yr4ejs3DWtmnaGRT0kNzSULJA0ADgQeA64Djs+THQ/8LndfBxwlaX1J2wJjgKmNiq+vcNOwZtZZGlnCGAZcku9jrANcFRF/kHQXcJWkE4CngSMBImKGpKuAR4CVwMn5kpatBTcNa2adRRE99zZAS0tLTJs2rdlhmJn1KJKmR0RLe+fzm95mZlaKE4aZmZXihGFmZqU4YVjDuOJDs96ldMKQtGEjA7HexxUfmvUubSYMSf8k6RHg0dw/VtKkhkdmPZ4rPjTrXcqUMH4EHAwsBoiIB4B/bmRQ1nu44kOz3qPUi3sRMUdaraonX5S2UlzxoVnvUSZhzJH0T0BIWg84lXx5yqwtrvjQrPcokzA+A5xHqmp8LnAjcHIjg7LewxUfmvUeZRLGgIg4tjhA0pYNisd6kUrFh6ceMNrJwqwXKHPTe5akK3KNsxXXNyog6z1c8aFZ71ImYTwE/C9wu6Tt8rBajR2ZrWaX4Zsy6dg92KC/m2Y36w3KXJKKiJgk6QHg95K+RI2W8MzMrHcrkzAEEBF3SBoP/Bp4W0OjMjOzbqdMwji00pGbVT0A+KfGhWRmZt1RqwlD0kcj4lfA0VUv7VXc1rCozMys26lXwqhUNujnIc3MrPWEERG/yH+/2ZEFSxoJXApsCbwBTI6I8ySdBXwKWJgn/XJEXJ/nORM4gVT1yKkR8eeOrNvMzDpfmdpqvy9pE0n9Jd0saZGkj5ZY9krgPyJiR+BdwMmSdsrjfhQR4/Knkix2Ao4CdgYOASZJ8vOYZmbdRJn3MN4TES8C/0KqGmR74AttzRQR8yPi3tz9Eqn+qeF1ZjkMuDIilkfELGAmsGeJ+MzMrAuUSRj9899DgSsi4vn2rkTSKGB3oFJd6SmSHpR0kaTBedhwYE5htrnUTzBmZtaFyiSM30t6DGgBbpY0FHi17AokbQRcA3w2l1TOB7YDxgHzgR9UJq0x+xovCEqaIGmapGkLFy6sMYv1dW4a1qwx2kwYEXEGsBfQEhErgGWky0dtktSflCwui4jf5OU9FxGvR8QbwC9ZddlpLjCyMPsI4Jka8UyOiJaIaBk6dGiZMKyPcdOwZo1Rqk3viFgSEa/n7lci4tm25lF6eeNC4NGI+GFh+LDCZEcAD+fu64CjJK0vaVtgDDC13GaYreKmYc0ao1SLex20N/Ax4CFJ9+dhXya9CDiOdLlpNvBpgIiYIekq4BHSE1YnV5KUWXt95X078tfHFnDchVM598ixFF8+HTF4ACM3G9jE6Mx6JkW0Xo9gLiWMiIg5rU7URC0tLTFt2rRmh2Hd1NGT7+auJxevMdxNw1pfJ2l6RLS0d766JYyICEm/BfboaGBmzeCmYc06X5lLUndLekdE/K3h0Zh1EjcNa9b5yiSM/YHPSJoNvEJ6/DUiYrdGBmbWUW4a1qwxyiSM9zY8CrNO5KZhzRqjzHsYT5Hejzggdy8rM59Zs7hpWLPGKFP54DeALwFn5kH9gV81MigzM+t+ypQUjgA+QLp/QUQ8g9vIMDPrc8okjNcivawRAJI2bGN6MzPrhcokjKsk/QIYJOlTwF+ACxoblpmZdTdtPiUVEedKOgh4EdgB+HpE3NTwyMzMrFtpM2FI+l5EfAm4qcYwMzPrI8pckjqoxjC/m2Fm1se0WsKQdCJwEvBWSQ8WRm0M3NHowMzMrHupd0nqcuAG4LvAGYXhL3WkmVYzM+vZWr0kFREvRMRs4KvAs/kt722Bj0oa1DXhmZlZd1HmHsY1wOuSRpNa0NuWVPowM7M+pEzCeCMiVgIfBH4cEacDw9qYx8zMepkyCWOFpKOB44A/5GH9GxeSmZl1R2USxieAvYCJETFL0raUqHxQ0khJt0h6VNIMSafl4ZtJuknS4/nv4MI8Z0qaKenvkg7u6EaZmVnnq9um91otWBoGDIuIeyVtDEwHDgc+DjwfEWdLOgMYHBFfkrQTcAWwJ7AVqQqS7SPi9dbW4Ta9zczaryFteucFzyJXPFgUEW+tN19EzAfm5+6XJD0KDAcOA/bLk10CTCFVn34YcGVELAdmSZpJSh53ldwWMzNroDIt7hWz0AbAkcBm7VmJpFHA7sA9wBY5mRAR8yVtnicbDtxdmG1uHmZmZt1AmRb3Fhc+8yLix8ABZVcgaSPSo7mfjYgX601aa/U1ljdB0jRJ0xYuXFg2DDMzW0tlLkm9vdC7DqnEUaoBJUn9Scnisoj4TR78nKRhuXQxDFiQh88lNQVbMQJ4pnqZETEZmAzpHkaZOMzMbO2VuST1g0L3SmA28JG2ZpIk0ot+j0bEDwujrgOOB87Of39XGH65pB+SbnqPAaaWiM/MzLpAmfYw9u/gsvcGPgY8JOn+POzLpERxlaQTgKdJ90SIiBmSrgIeISWmk+s9IWVmZl2rXm21n6s3Y1Wpodb426l9XwJgfCvzTAQm1luumZk1R70SRqn7FGZm1je0mjAi4ptdGYiZmXVvbT5WK+mSYnXmkgZLuqihUZmZWbdTpi6p3SJiaaUnIpaQXsIzM7M+pEzCWKeqgsDNKPc4rpmZ9SJl38O4U9LVpDevP4KfZDIz63PKvIdxqaRppOpABHwwIh5peGRmZtatlLq0lBOEk4SZWR9W5h6GmZX08LwXOOmy6by6wpUUWO/jhGHWiZYuW8H1Dz3LpFtmNjsUs05XKmFI2kbSgbl7QG5Bz8yq7DNmCEfsPpzzb32CmQteanY4Zp2qzIt7nwKuBn6RB40AftvAmMx6tK+8b0cGrrcux104lTtnLuKuJxa/+Znz/LJmh2fWYWVuep9Mair1HoCIeLzQSp6ZVRmy0frsNGwT7npyMcdccM9q404bP4bTD9q+SZGZrZ0yCWN5RLyWmrcASetSoyU8M0sWvbycR+a/yPBBAzjnyN1QodLmEYMHNDEys7VTJmHcKunLwABJBwEnAb9vbFhmPdfEPz7KstdWcs2JezF6c9/us96jzE3vM4CFwEPAp4Hrga82Miiznur2xxdx7X3zOHHf7ZwsrNcp86b3G8Av88fM6hg0sD+H7rolJ+0/utmhmHW6NhOGpIdY857FC8A04NsRsbgRgZn1RLsM35RJx+7R7DDMGqLMJakbgD8Cx+bP74HbgGeBi1ubSdJFkhZIergw7CxJ8yTdnz+HFsadKWmmpL9LOriD22NmZg1S5qb33hGxd6H/IUl3RMTekj5aZ76LgZ8Cl1YN/1FEnFscIGkn4ChgZ2Ar4C+Sto8I169gZtZNlClhbCTpnZUeSXsCG+Xela3NFBG3Ac+XjOMw4MqIWB4Rs4CZpHc/zMysmyhTwvgkcJGkjUjVm78IfFLShsB3O7DOUyQdR7oH8h+5Bb/hwN2FaebmYWZm1k20WcKIiL9FxK7AOGBcROwWEVMj4pWIuKqd6zsf2C4vaz6pcSag8GZTYdW1FiBpgqRpkqYtXLiwnas3M7OOKtUehqT3ke4vbFB54zsivtXelUXEc4Vl/hL4Q+6dC4wsTDoCeKaVZUwGJgO0tLT4jXMzsy5SpvLBnwP/Cvw7qSRwJLBNR1YmaVih9wig8gTVdcBRktaXtC0wBpjakXWYmVljlClh/FNE7CbpwYj4pqQfAL9payZJVwD7AUMkzQW+AewnaRzpctNs0pvjRMQMSVeRWvVbCZzsJ6TMzLqXMgnj1fx3maStgMXAtm3NFBFH1xh8YZ3pJwITS8RjZmZNUCZh/F7SIOAc4F5S6cDVhJiZ9TF1E4akdYCbI2IpcI2kPwAbRMQLXRGcmZl1H3VveueKB39Q6F/uZGFm1jeVedP7RkkfUuV5WjMz65PK3MP4HLAh8Lqkf5AerY2I2KShkZmZWbdSpj0MtwJjZmalXtyTpI9K+lruH5krIDQzsz6kzD2MScBewDG5/2XgZw2LyMzMuqUy9zDeGRFvl3QfQEQskbReg+MyM7NupkwJY4WkfuTaYyUNBd5oaFRmZtbtlEkYPwGuBTaXNBG4HfhOQ6MyM7Nup8xTUpdJmg6MJz1Se3hEPNrwyMzMrFtpM2FIOg/4dUT4RreZWR9W5pLUvcBXJc2UdI6klkYHZWZm3U+ZJloviYhDgT2B/wO+J+nxhkdmZmbdSpkSRsVo4G3AKOCxhkRjZmbdVpk3vSslim8BM4A9IuL9DY/MzMy6lTIv7s0C9oqIRY0OxszMuq8yj9X+XNLgXH/UBoXhtzU0MjMz61bKXJL6JHAb8Gfgm/nvWSXmu0jSAkkPF4ZtJukmSY/nv4ML487MT2L9XdLBHdkYMzNrnDI3vU8D3gE8FRH7A7sDC0vMdzFwSNWwM0hNvo4Bbs79SNoJOArYOc8zKVdHYmZm3USZhPFqRLwKIGn9iHgM2KGtmfIlq+erBh8GXJK7LwEOLwy/MjcBOwuYSXqM18zMuokyN73nShoE/Ba4SdIS4JkOrm+LiJgPEBHzJW2ehw8H7i6uMw9bg6QJwASArbfeuoNhmJlZe5V5ce+IiFgaEWcBXwMuZFXJoLPUai88WolnckS0RETL0KFDOzkMs57p4XkvcNJl03l1xevNDsV6sfa8uEdE3BoR10XEax1c33OShgHkvwvy8LnAyMJ0I+h4Kcasz1m6bAXXP/Qsk26Z2exQrBdrV8LoBNcBx+fu44HfFYYfJWl9SdsCY4CpXRybWY+1z5ghHLH7cM6/9QlmLnip2eFYL1XmHkaHSLoC2A8YImku8A3gbOAqSScATwNHAkTEDElXAY8AK4GTI8Jla7N2+Mr7duSvjy3guAuncu6RY5FWXekdMXgAIzcb2MTorDdQRM1bBT1CS0tLTJs2rdlhmHUbR0++m7ueXLzG8NPGj+H0g7ZvQkTWHUmaHhHtrnm8YSUMM+tai15eziPzX2T4oAGcc+RuiNVLGGZrywnDrJeY+MdHWfbaSq45cS9Gb75xs8OxXqirb3qbWQPc/vgirr1vHifuu52ThTWME4ZZLzBoYH8O3XVLTtp/dLNDsV7Ml6TMeoFdhm/KpGP3aHYY1su5hGFmZqU4YZiZWSlOGGZmVooThpmZleKEYWZmpThhmJlZKU4YZmZWihOGmZmV4oRhZmalOGGYmVkpThhmZlaKE4aZmZXihGFmZqU0pbZaSbOBl4DXgZUR0SJpM+DXwChgNvCRiFjSjPjMzGxNzSxh7B8R4wrtyp4B3BwRY4Cbc7+ZmXUT3emS1GHAJbn7EuDw5oViZmbVmpUwArhR0nRJE/KwLSJiPkD+u3mtGSVNkDRN0rSFCxd2UbhmZtasFvf2johnJG0O3CTpsbIzRsRkYDJAS0tLNCpAMzNbXVNKGBHxTP67ALgW2BN4TtIwgPx3QTNiMzOz2ro8YUjaUNLGlW7gPcDDwHXA8Xmy44HfdXVsZmbWumZcktoCuFZSZf2XR8SfJP0NuErSCcDTwJFNiM3MzFrR5QkjIp4ExtYYvhgY39XxmNnaeXjeC0yaMpMffmQcG/Tv1+xwrIG602O1ZtYDLV22gusfepZJt8xsdijWYE4YZrZW9hkzhCN2H875tz7BzAUvNTsca6BmPVZrZr3IV963I399bAHHXTiVc48cS75HCcCIwQMYudnAJkZnncUJw8zW2pCN1menYZtw15OLOeaCe1Ybd9r4MZx+0PZNisw6kxOGma21RS8v55H5LzJ80ADOOXI3xOolDOsdnDDMbK1N/OOjLHttJdecuBejN9+42eFYg/imt5mtldsfX8S1983jxH2369Jk8fC8Fzjpsum8uuL1LltnX+eEYWZrZdDA/hy665actP/oLl2vH+ftek4YZrZWdhm+KZOO3aPLX9rz47xdz/cwzKzH8uO8XcsJw8x6LD/O27WcMMysx2rG47x9ue4sJwwz67Ga8Thv5Wb76KEz+dx7duiSdVY0O1n5preZ9UjNepy3mTfbm/1kmBOGmfVIzXqcF9LN9oHrrctxF07lzpmLuOuJxW9+5jy/rGHrbfaTYYrouc1it7S0xLRp05odhpn1QUdPvpu7nly8xvBG32xf9PJyxv/gVjZcr1+HnwyTND0iWtq7bt/DMDNrp2bWndXMJ8OcMMzM2qmZdWc1M1l1u4Qh6RDgPKAfcEFEnN3kkMzM3lS52X7qAaObUtFiM5NVt7rpLakf8DPgvcBOwNGSdmpuVGZmqzTzZnuzngyr6FYJA9gTmBkRT0bEa8CVwGFNjsnM7E3NqjsLmpusoPsljOHAnEL/3DzsTZImSJomadrChQu7NDgzs2ZqZrKC7pcwVGPYas/9RsTkiGiJiJahQ4d2UVhmZtbdEsZcYGShfwTwTJNiMTOzgu6WMP4GjJG0raT1gKOA65ock5mZ0c0eq42IlZJOAf5Meqz2ooiY0eSwzMyMHl41iKSFwFPNjqMDhgCLmh1EF/M29w19bZt76vZuExHtvgncoxNGTyVpWkfqcenJvM19Q1/b5r62vd3tHoaZmXVTThhmZlaKE0ZzTG52AE3gbe4b+to296nt9T0MMzMrxSUMMzMrxQnDzMxKccJoIEmHSPq7pJmSzmhlmv0k3S9phqRbuzrGztTW9kraVNLvJT2Qt/cTzYizM0m6SNICSQ+3Ml6SfpK/kwclvb2rY+xsJbb52LytD0q6U9LYro6xs7W1zYXp3iHpdUkf7qrYupITRoOUadtD0iBgEvCBiNgZOLKr4+wsJdsyORl4JCLGAvsBP8hVwPRkFwOH1Bn/XmBM/kwAzu+CmBrtYupv8yxg34jYDfhPeseN4Yupv82VY+B7pJoqeiUnjMYp07bHMcBvIuJpgIhY0MUxdqYy2xvAxkqt1m8EPA+s7NowO1dE3EbajtYcBlwayd3AIEnDuia6xmhrmyPizohYknvvJlUi2qOV+D8D/DtwDdCTj+O6nDAap822PYDtgcGSpkiaLum4Louu85XZ3p8CO5JqIH4IOC0i3uia8JqmzPfSm50A3NDsIBpN0nDgCODnzY6lkbpV5YO9TJtte5C+/z2A8cAA4C5Jd0fE/zU6uAYos70HA/cDBwDbATdJ+t+IeLHBsTVTme+lV5K0Pylh7NPsWLrAj4EvRcTrqQDdOzlhNE6Ztj3mAosi4hXgFUm3AWOBnpgwymzvJ4CzI738M1PSLOBtwNSuCbEp+mQbL5J2Ay4A3hsRi5sdTxdoAa7MyWIIcKiklRHx26ZG1cl8SapxyrTt8Tvg3ZLWlTQQeCfwaBfH2VnKbO/TpNIUkrYAdgCe7NIou951wHH5aal3AS9ExPxmB9VIkrYGfgN8rIeWltstIraNiFERMQq4GjiptyULcAmjYVpr20PSZ/L4n0fEo5L+BDwIvAFcEBF1H9vrrspsL+mJmYslPUS6VPOliOiJVUO/SdIVpCe+hkiaC3wD6A9vbvP1wKHATGAZqZTVo5XY5q8DbwEm5TPulT29RtcS29wnuGoQMzMrxZekzMysFCcMMzMrxQnDzMxKccIwM7NSnDDMzKwUJwwzMyvFCcPMzEpxwrBeQdIgSSc1O462NDpOSZ/K7avcL+mNQvcPG7VO6zv84p71CpJGAX+IiF26QSwiHVtr1MTbkTjrLa/OPMOBOyNim7LzmLXFJQzrLc4Gtstn0+dI+qikqbn/F7lxGySNkvSYpAskPSzpMkkHSrpD0uOS9ixMc0luNe7qXNcXeRlrLDvP86ikScC9wEhJv83V1s+QNKGVOEcVW3GT9HlJZ7WyvJrb1IpdSFXIm3UaJwzrLc4AnoiIccBFwL8Ce+f+14FjC9OOBs4DdiPVlnsMqQruzwNfztPsAEzOrca9CJwEIGnHOsvegdRY0u4R8RTwbxGxB6km01MlvaUYZ0R8oY1tenN5wMA2tqnarkCPrJfMui9XPmi90XhSOyN/y5XfDWD1VtBmRcRDAJJmADdHRORKEUflaeZExB25+1fAqcC5dZZ9G/BUblWv4lRJR+TukaRmWp9tx3YUl9fWNlXbBbipHesya5MThvVGAi6JiDNbGb+80P1Gof8NVh0T1Tf3Kv01l53vTbxS6N8POBDYKyKWSZoCbFAjlpWsXtIvTvNKobutbaq2K/CjktOaleJLUtZbvARsnLtvBj4saXMASZtJau/N360l7ZW7jwZub+eyNwWW5GTxNuBdNeIEeA7YXNJbJK0P/Esr8ZTeJknrkEozj5XZULOynDCsV8itut2RbyB/AvgqcKOkB0mXZoa1c5GPAsfn+TcDzs/reaTksv8ErJun+U/g7uo4JZ0TESuAbwH3AH+glR/5dqwX0j2auRGxvJXxZh3ix2rNqnSnR3TNuhOXMMzMrBSXMMzMrBSXMMzMrBQnDDMzK8UJw8zMSnHCMDOzUpwwzMysFCcMMzMrxQnDzMxK+f9xUyRzxADgswAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEtCAYAAAAbeVcBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABAkElEQVR4nO3dd3wc1bn/8c9XxZLcJBdZ7r1hDBhsWugYgyEQ05IAISEkgZDAhUAS2g/S6KSSm0BCC3BDINyEUBya4WJ6s4lxNxj3glxwL7IsPb8/zgjW8koayVqtVnrer9e+tDv1mdXMPnNmzpwjM8M555yrS1a6A3DOOZcZPGE455yLxROGc865WDxhOOeci8UThnPOuVg8YTjnnIvFE4ZzzrlYPGE455yLxROG2yOSFkk6Lt1xxCVplqSjG3mZwyT9R9ImSZc25rIT1pFR33NLJukBSTfGnLZF/d88YaSIpMmS1knKS3cszV1THlRmtreZTW7kxV4JTDazDmb2+z1dWEv7kYlDUvtou89JGNZB0hJJZ0p6WNL91eY5StJaST2aPuLWyRNGCkjqDxwBGPClFCw/p7GXmamayXfRD5jVkBmbSfy1UrB/kuH7ScpujHWY2WbgQuAOScXR4NuBKWb2D+BS4CRJ46J15wP3AD80s5WNEYOLwcz81cgv4CfAG8BvgIkJw68G/lFt2juA30fvewL/BFYDC4FLE6ZbBFwFTAfKgJxoeR8Dm4DZwGkJ0x8A/Cca97/A34EbE8bXuK4k29MHeDyadi3wh2pxHRe9N2BwwrgHqq3zKmB5FNM8YCzwP0AlsA3YDFzZkO8iScy7rStJvF+N1ln1KiOUFGJ/P8D/ARXA9mgZQ4G9gMnAekIi+VK1eWqMv5bvYxHwo2ieDdH/M7+B/8/NQJ/o/UXR/60k+vxj4N5q0w+Iljs+YdgRwBpg5J4uv9q8DwCPAEcT9rUeCeO+HG1bO+AW4NlalrMoWtd0YAtwH1ACPBvtEy8CnRKmr/F/BuwPvB/N93fgUWIeSyTsby3hlfYAWuILmA98HxgNlCccLP2ArUDH6HM2sBI4hFDam0pINm2AgcAC4IRo2kXANMKPd0E07MvRzppF+PHbAvSI5l8MXAbkAqcDO6p28rrWVW1bsoEPgN9GB2o+cHjC+M8OCGpJGMAwYCnQM/rcHxhUfRlx4kv2XVSLOfa6EubpCMwBvluf7yeadzLwneh9bvT/vzaa91jCD82wat9ZbfHvFmM07N3o/905ivWi+v4/o+mXAiMAATOAj4Dh0ef5wH5J5jmS8AN+DHAQ1RLIni4/Yd5OhGNiDXB+kvH/AJ6KYulby3IWAW8TkkQvYBXhR39/II+Q6H9a1/+Mz4+ly6PpziQc07GOpZr2t0x9pT2AlvYCDo92qK7R57nA5QnjXwe+Eb0fB3wcvT8YWFJtWdcAf4neLwK+Vce6pwETooN7OaBq670xzrqqDT80+nHY7Sw+Ia44CWNwdNAeB+TWtIzG+C7qs65oWBYwEbirvt9PNG4ynyeMI4BPgKyE8Y8AP6sWQ23xJ4txEXBuwufbgT81MN4ZhB/9E4B/Ra9DgBOBV2qJ6wTCD/UqEkqzjbX8hPlfJJxYFSYZV0IowVxWxzIWAV9L+PzPqv9v9Pm/gCfq+p8RjqUV7HosvUnMYynZ/zKTX34Po/GdB7xgZmuiz3+LhpHw+ezo/TnRZwilj56S1le9CGc8JQnzLk1ckaRvSJqWMP1IoCvhLHS5RXtsknnjrKtKH2Cxme2se9NrZmbzgR8QDsJVkh6V1LOGyev9XezBugBuAjoQrpPHXX9NegJLzawyYdhiwlluohrjr8UnCe+3Au0bGO+6aN4fEC6JbiSc2X8fqO2m/RJgJ6GksKiW6Rq6fCSdSygRvgjcVn28mZUSSh9x7hmVJrzfluRz1fdX2/8s2bG0OOH9nuwrGafZ33DLJJIKgK8A2ZKqDu48oEjSfmb2AeF+wq8l9QZOI5zBQ/gBWWhmQ2pZxWc7raR+hJt+Y4G3zKxC0jTCwbwS6CVJCTt6H8L9jrjrqrIU6CspJ0bS2Aq0TfjcHVj2WfBmfwP+Jqkj8GfCD8LXE7erHvFVn2fXkTWvaxeSziIk8APNrLwe66/JCqCPpKyEH6C+wIf1iL/WbUuivvGuB8YQ7g9MlnQGsB/hhOOJZDNIGgRMItx72QQ8I+k4M0v2w13v5Ufr6Ea49PkVQsl8lqS/mdmrMberoWr7nyU7lvrSsGMp43kJo3GdSrgBOgIYFb32Al4DvgFgZqsJlzD+QtjR5kTzvgtslHSVpAJJ2ZJGSjqwhnW1I/ywrAaQdD7hgAR4K4rjEkk5kiYQLhFUqc+63iUcNLdKaicpX9JhNcQ0DTgnWt544KiqEdGzCsdG1Yy3E87wKqLRpYRrvw2Jbzd1rCtxuv2B/wZOjf4vjbH+dwj3kq6UlKvwzMcphBulcVX/PupS33jXEa7JV53tbyTc7/qzmSX7nnoCLwE3mdkDZvZPwg34FyQli7Ney0/wB8Jlopct1Hy6ErinCaqm1/Y/e4tQqro0OpZOp+HHUsbzhNG4ziNcu1xiZp9UvQgHwtcSqlD+jXB9vepyFNGBdAohySwkFLvvBQqTrcjMZgO/JuzQpcA+hJpZmNkOwo3ubxPO9s4lXKMvq++6EqYdTLgksYxwgz2Zy6Jp1wNfY9ezyTzg1mhdnwDdCEV3CDVerouK9D+q73eRRG3rSjSBcKnkdUmbo9eze7L+6Lv/EuF6/RrgTsI9q7kxY4dq30eMddY33nWEqwtV+99GoIhQYk1mLaH66l0J63yYcIlpVSMsH0mnEu7//ThhHfcS9ref1DRfY6jtf5ZwLH2TsF1fJdQYrJp3T/fVjKJdL825lkrSO4SbpH9JdyzOuczkJYwWSuEp2O5RMfo8YF/guXTH5ZzLXH7Tu+UaBjxGqAnyMXCm+ROxzrk94JeknHPOxeKXpJxzzsXiCcM551wsnjCcc87FkjEJQ9Itkn6Q7jhaGkl9o+cPsqPPkyV9J3r/NUkvJExrkgY3VSzOudST9CVJsR4szYiEodA+/jcITTykY/0/ljRToUe1hZJ+XG38IknbEh7+eqGmZdWxnsckjZOUl9C0SEpFDxm2T/YErpk9bGbHp2rdqtZRUG2xNFeSOktaLen1asNHSZoqaWv0d1S18ZdL+kTSBkn3N+RpZkkjJL0Rvf+FUtTbX2uUeOLUhOu8RNIUSWWSHkgyfqykudE+9bJC80BV4yTpNoUOpdZKul2SEsb3j+bZGi0j8bh7Chgpad+6YsyIhEF4yvIZM9uWpvWLkLA6AeMJTW6cVW2aU6Ifu/Z78CM7mtAE877AzAZH2wwoAzoGaiS3EZoa/4ykNsCTwF8J+8yDwJPRcCSdQOjLZCyhob2BwM8bsO7RhKa1q96/34BlNBtesmQFcCNwf/URkroSnjC/ntC8/RRC3xxVLiQ0TbQf4ffjZEJT/VUeIfSP0wX4f8A/9HlHVVXjL6wzwnQ3lxvnRWi7PrFp51eAM6L3hxPaVDop+nwcMC3F8fwe+O+Ez4vYwyaMCT8sC6L33wNurza+kNAJzEpC0+U3AtnRuGzgV4RmCRYAF0ffSU6y+AituP41et+/2rST+byp7m8CryfMZ4QWXRdE6/olUZPQ0bRvEBqP+zSKb1D0v1sbTf8wUBRNv1tHQUli6Uno++BTQn8FF1TbhseAhwiN4c0CxiSMT9qBUiPvB4cSmmY5v9r3dDy7Ny+/hKj/CEKTGTcnjBsLfNKA9d8BnBe9XwG0Txh3NKFZjSsJzXesJPygnERoVO9T4NqE6bP4vEOutdF32zlh/P8SmlnZALwK7J0w7iRCB16bou3+UbL9J2EfGhy9fwC4C3iG0JbTcdTeGdHPojj+Gq1rBqHDqmuibVwKHB/zmPkmocn/XxGa/FgInBiNu4ldO8X6A+Gk8bfRejYQOmbarfOoRtqvbgQeqDbsQuDNhM/tCMfO8Ojzm8CFCeO/DbwdvR9KaBaoQ8L414j6U4k+H0Zo267W2DKlhLEP4aCv8grhgIDQXv0CPm/o7sho/G4knaOEZoiTvPrWFUhUzDuC3ZtXfji6NPGCpP3iblhUzFxP2Nl7R+/vAC6OYqrargcJjaANJnQCczxQVWS+gHBGsT+hldAz466/AU6L1nEAoS2mbyWMO5jwv+hGOOhEaBepJ6ERxj6Egx4z+zrhR7SqZHZ7knU9QvjR60nYppsljU0Y/yVCA3FFhMTyBwiNDwKXEFqg7UDol2FRso2RdHVt+0RNX0J0NvzHaD3VH2baG5hu0ZEYmR4Nrxr/QcK4D4ASSV1qWl+1dU+KYrsY+G9JGwnNaS+T9GzCpN0JHV71IrTHdA+hXbHRhH34J/q88cBLCQnlKML3vS7avirPAkMI/9v3Ccm/yn3Ad6PveiThJCGuc/i8efk3gacJ30cvQiL9QVQiq3IK4WSjE+GM+XlCsusF/IJdL1vXdsxA2F/nEboEuB24T5LM7P8RflAvifbNS6J5jyT8+BYR2pRam2yDJN1Zyz41vR7fTaJd9hkz20JI7rXtU4njFpjZphrGQygl91do3blmqciQKci45USZ1D4/I5sevX+OsBNUZdNXgNNTGMvPoy87r1p2LiA07X0N4UysqJ7L/SvhB7gd4QywbcK4EsIZQkHCsLOBl6P3/8euZwvHk7oSRmI3nd8HXkqYdkkd23gq8J+Ez9Xj+iwWQnKpYNezoluIzryibXgxYdwIYFv0vsYOlBpxP7iczztcqv49XQ88Wm36h4k6USIc6InfY2603f3rsf6hhP6uITSs+ONq448mnIFWnVF3iNZxcMI0Uwkt9UL4wRibMK4H4bhL1v1tUbSswujzEsLlj47Vptvle0nYhxJLGA8ljKurM6KfAZMSxp1CKAFU38Yi6j5mvgnMTxjXNpq3e/XjIPp8LOG4PISEjpZStG8lK2HcB9xabdgbwDej9xXs+hs5JNoeEZr1f7vavDclriNhH6yxF0OzzClhrCPsDFXeAoZKKiG0EvkQoT37roSmh1PSfr6kSwj3Mr5oZmVVw83sDTPbZmZbzewWQmutR8Rc5rLobPFswhnRKkKnLCsk/SaarB/hH7oy4cz3z4SzPYg6gElYbGIHL42t+np61jAOSd0UOi9aHp0F/5VwNhdHT+BT2/WsqHpHRNU7FMpX6Lejvh0o1Uu0rEsJ14KT2Uzo8jVRR8JllGTjq95vog7RjdH1RGeI0fsb+Lx1224Jk6+1zysQVN3/q6kToX7AvxL2rzmEH6EShSa7b5X0cfR/XBTNU/W/PINwWWqxpFckHUp89e3Yq3r8a5JsY3vqPmYgYf8xs60J8+7GzP6PUIL9I1Aq6e46z8YbV0P2qc0WskFd88Lnv6/rawsiUxLGdMIZFfDZP3cqoTntmRaaIH4TuILQ5emaZAtRqCa6uZZXjZekJH2L6EalmS2rabqqEAmZvU5m1ptwI/1FMysC7gYuNrMiM7simmwp4WypazS8yMw6mllVkXIl4Yy8SvXt2MLuHRs1VPX1rEjcnGrT3hIN29fMOhIuh6iW6ROtADpLSjxR6Eu4Fl0nM/ubmR1O+OEwkvTeBiDp2tr2iRoWfxDhDHy2Qm22O4CDFGo9ZRMuV+6bWEuFcCOy6jLmLMLNySr7AaVmlvQSR7Xt+kO0n7xCOOvtR+gRrjDaL5I1Nx7HUsI1/KKEV76ZLSdcNppAKLEVEkqCEP0vzew9M5tA+DF+gnD/A6rtd5KS7XeJ+0BVZ0SJMXQws5MauD21HTN12W3fNLPfm9lowqWcoSQ0xZ5I0p9q2afi9BSYzC77jKR2hHuEte1TieMGVjuWEsdDuGS8yMw21hZEpiSMZ0jojCfyCuH6cdX9isnVPu/GQjXR9rW8liSbT9LXgJuBcWa2oNq4vpIOk9RGoXOhHxPOvKqqO/ZXeH6hfy3bl1jD5QBCDYjEuFcCLxB66usoKUvSoIT7G48ROnjpLakTIbElmgacpdA5zJ7e4/ixpE6S+hAS9t9rmbYD4exmvaRe7H6A1dhRkJktJZwE3BJ9r/sSbuQ9nGz6RIrZgVK0nptr2ydqWMWzhB/NUdHrJ4Tr6aOis93J0fouVagifUk0X9W1/YeAbytUi+0EXEe4PFMV/wNKUq2ymv0IpYwDaJzaUX8CblJUVVNSsULHWxD+j2WEa/ZtCcdCVaxtohOxQgs9Fm7k8++6qhQ0SlI+0f2rWjRaZ0Qxjpm67LJvSjpQ0sGScgmJcDs171MX1bJP1ZiwFFqWzidUYsmO9vuq2ob/IlR9PSOa5ieEy/JV/aw8BFwhqVdUAv4h0T5lZh8SfgN+Gi3zNMIJzD8TVn8UYb+uVaYkjIeAkxS6QK3yCmFHfrWGz43pRkJ1tPcSzhT+FI3rQKjpsY5w9juecKZWdbbYh3AppbYz49HA+9EZ6XCS91f8DaANoTbKOuAfhLNcCDcznyccoO+T0MFL5HrC2cg6wj2Yv9FwTxJKd9OAfxOurdbk54QftA3RtNXjqqujoLMJP8wrCAfMT81sUowY43ag1CBmVma7dpC1ASiP3hOVeE8l/M/WEyoGnBoNx8yeI9xkfZmwbywGfpqwij5EJxzJRCXhT6OS9gF8XrV2T9xBqDjwgqRNwNuEewoQjr+qfXh2NC7R14FF0eWqiwglyaofql8Q+uf+iFArqUbW+J0R1XbM1OUO4ExJ6yT9nnAJ555oOYsJyfNXDYyrJtcRTm6uJnyH26JhWOgR8gzCvYd1hP9NYtX+PxMqDMwgVMn/N7tWADiLUFllHeHYONN27WXybGI855YxrdVKuhlYZWa/S3cs9SHpOmC1mTXZQ4dRaWYh4YZvXf1wu2ZE4VmNDwiX8crrmt65PSXpFODrZvaVOqfNlITh4vOE4ZxLhUy5JOWccy7NvIThnHMuFi9hOOeciyWjG4jr2rWr9e/fP91hOOdcRpk6deoaMyuue8pdZXTC6N+/P1OmTKl7Quecc5+R1KDWIFrNJam7Jn/MkrVb657QOedcUq0mYRjGdU/OxG/yO+dcw7SahHHBEQMp3bCdidNXpjsU55zLSK0mYeRmZ3Hz6SO5YeJsNmzzB2idc66+Wk3CABjdrzPHjSjh9ufm1j2xc865XbSqhAFw1QnDmTS7lKmL16U7FOecyyitLmEUts3lupNHcO3jMyivqEx3OM45lzFaXcIAOGXfHpQU5nPvawvTHYpzzmWMjH5wr6EkcciAztz23Fxui3E/47KxQ7h83NA6p3POuZasVSYMM+PthZ9y3qH9OPvgvrRrk0Ofzm3rntE551qxVnlJ6unpKyndsJ3rTh5B/y7teHLacj4s3eQP9TnnXC1aXcLYsLWcGyfO5ubT9yE3O4v83GwuOXYIvTsV8IuJs1m+flu6Q3TOuWap1SWM256fy7gRJYzu12mX4W3b5HDV+OFs2l7Oy3NXeWnDOeeqaVUJY+riT3lxdilXjh+edHx+bjbDu3ekfX4OD7+zxJOGc84laDU3vcsrKrn28Zlcf/IICgtya532wP6dObB/Zx55dwl5OVmctn8vJDVRpM451zy1mhLGPa8toHthPifv2yP2PGcf1Je9exYyZfE6Vm3cnsLonHOu+Ws1CUOIGyaMrHdJYVj3Dozo0ZG/vr2YDdvK/TKVc67VajUJ43tHD6Jvl4Y9a9EuL4crjh/G2s1l3PTvOazZXNbI0TnnXPOXsoQhKV/Su5I+kDRL0s+j4T+TtFzStOh1UsI810iaL2mepBNSFVtDDSxuzw/GDWXbjgqeneH9ajjnWpdU3vQuA441s82ScoHXJT0bjfutmf0qcWJJI4CzgL2BnsCLkoaaWUUKY6y39nk5tM/LYcGaLbw8dxWHDe5Km5xWU1BzzrViKfuls2Bz9DE3etV2A2AC8KiZlZnZQmA+cFCq4ttTRw0t5pjh3fjjy/OZNLs03eE451zKpfTUWFK2pGnAKmCSmb0TjbpE0nRJ90uqeoKuF7A0YfZl0bDqy7xQ0hRJU1avXp3K8GO5fNxQ9ulVyMTpK1i3Zcdu4++a/DFL1m5NQ2TOOde4UpowzKzCzEYBvYGDJI0E7gIGAaOAlcCvo8mTVV/arURiZneb2RgzG1NcXJySuOure2E+Rwwu5rEpS9levusVNMO47smZXrvKOZfxmuTiu5mtByYD482sNEoklcA9fH7ZaRnQJ2G23sCKpoivMRS2zeW7Rw1i2tL13PrsXLbtCInjgiMGUrphOxOn+01y51xmS2UtqWJJRdH7AuA4YK6kxCfnTgNmRu+fAs6SlCdpADAEeDdV8aXKIQO7cNFRA1m6bisvz1tFbnYWN58+khsmzmbDtvJ0h+eccw2WyhJGD+BlSdOB9wj3MCYCt0uaEQ0/BrgcwMxmAY8Bs4HngIubWw2puIratmFoSQfKyiv5sHQTQ0s6cNyIEm6P0VmTc841V8rka+tjxoyxKVOmpDuMWm3dsZPfTvqQwwZ35cp/TOeuc0fv1lKuc841JUlTzWxMfefzBwhSrG2bHP7fF0cwul8nDh/clav/OZ3yisp0h+Wcc/XmCaOJdMjP5eoTh5OdJX734ofpDsc55+qt1TRv3hx065jPdV8cwbn3vcMfX/64zukvGzuEy8cNbYLInHOubp4wmpCZcfdrC7hq/HAuOGIAUxev4+PVWzjrwD5kZXl/G8655s0TRhN6evpKSjds5ztHDCAnO4uDB3ahXV4OH63azM7KSvbuWZjuEJ1zrkZ+D6OJbNhazo0TZ3Pz6fuQm/351z6yVyGDu7XnP0vW89bHa6mozNxaa865ls1LGE3ktufnMm5ESdIqtdlZ4txD+mFm/PL5eYzp34ljh5ekIUrnnKuZlzCawNTFn/Li7FKuHD+81ukkceX44YzoUciT05azfP22JorQOefq5gkjxcorKrn28Zlcf/IICgtyY83TvTCf40d0Z9KsT1i1aTs7/bkN51wz4Akjxe55bQHdC/M5ed8edU+coKBNNt88bACbtu/kxn/PYdXG7SmK0Dnn4vGEkWJC3DBhJFLDqs0OKm7PT04eAYJfvzDPGzB0zqWNtyWVQdZuLmPlhu2s3bKDI4d0bXAScs61bt6WVCvQpX0eI3sVkpstHnl3qd/bcM41KU8YGegLg7pyzsF9eXLaCu55dUHS3vy8a1jnXGPzhJHBzhjdmwmjejJ18TreXrB2l3HeNaxzrrF5wshw3TrmM7pfJ0o3bmfRmi18umUH4F3DOucanyeMFkASE0b1omuHPP7yxkLeXrDWu4Z1zjU6TxgtSPu8HH54/DAO7N+Z37wwj8KCXO8a1jnXaFKWMCTlS3pX0geSZkn6eTS8s6RJkj6K/nZKmOcaSfMlzZN0Qqpia+mys8R/jR3C2s07OHpoMS/MLmXq4nXpDss5l+FSWcIoA441s/2AUcB4SYcAVwMvmdkQ4KXoM5JGAGcBewPjgTslZacwvhYtN2o+fUTPjhzQt4hrHveuYZ1zeyZlCcOCzdHH3OhlwATgwWj4g8Cp0fsJwKNmVmZmC4H5wEGpiq+16N2pLX86dzTt8nI45563+WSDNzHinGuYlDZvHpUQpgKDgT+a2TuSSsxsJYCZrZTULZq8F/B2wuzLomHVl3khcCFA3759Uxl+iyGJH44bxrn3vcMht7xU5/TeNaxzLpmUJgwzqwBGSSoC/iVpZC2TJ2vnYreHCMzsbuBuCE2DNEacLV1i17DfO3oQZsbqzWW8t3Ada7eU8dUD+5CX41f/nHO1a5JaUma2HphMuDdRKqkHQPR3VTTZMqBPwmy9gRVNEV9Ll9g1LIQSR7cO+Xxx3x4cNbSY1ZvKuOPFj1i5wfvfcM7VLJW1pIqjkgWSCoDjgLnAU8B50WTnAU9G758CzpKUJ2kAMAR4N1XxtRY1dQ1bpV+XdvTu1JZvHtafjdt28vA7i5m+bH3TB+qca/ZSeUmqB/BgdB8jC3jMzCZKegt4TNK3gSXAlwHMbJakx4DZwE7g4uiSltsDtXUNm6iwIJfCglwGFrdj8rzVzFqxgaWfbmPciBKys7xVXOdcChOGmU0H9k8yfC0wtoZ5bgJuSlVMrU1V17CTrjgq9jy52VmMG1GCmbF1RwXzPtnE7JUbOXFkd9rlxdtd7pr8MV/cpwd9u7RtaOjOuWbIn/RuoRrSNWwiSRzYvzMjenbkoP6deWbGSqYtXc+KGP2Me8OHzrVMnjBaqIZ2DZtM3y5t+fKYPgzo2o5nZ37C4rVb+GDp+hqn94YPnWuZPGG0UHvaNWwyhQW5fPvwAfQqKmDJp1t5+oMVLFyzhYrKXUsS3vChcy2TJ4wW6ntHD0rZPYSc7CxO2a8np+zXk/Vbd/DL5+exbssONpft/Gya0f06e8OHzrUwsROGpHapDMRlpv37duLqE4dTacY9ry7gzY/XsCVKHFedMJxJ3vChcy1GnQlD0hckzQbmRJ/3k3RnyiNzGaVL+zwuHzeULwzqypPTVvDfL31E27xsrjt5BNc+PsMbPnSuBYhTwvgtcAKwFsDMPgCOTGVQLrOdc3BfLjl2MHNXbuKDJevokJ/Dva8tTHdYzrk9FKtivZktrXbz1B+oc7WSxD69CxnZqyOT563i/AemcFuM+xne8KFzzVechLFU0hcAk9QGuJTo8pRzcfzlzcX8+IShjN2rhEVrtnDU0G4UtPHGDp3LNHEuSV0EXExoanwZoTOki1MYk2tBqho+vPDIQQzv3pH9+hTx4FuL2LpjZ90zO+ealTgJo8DMvmZmJWbWzczOJXSG5FytkjV82KOwgIuOGsSMZRu46d+zd6mK65xr3uIkjIWSHolanK3yTKoCci1HbQ0fHjywC5ePG8qaTWU8/M7i3R7+c841P3ESxgzgNeB1SYOiYd58qatVVcOHV44fXuM0bdvk0L9rO0b1KWL2io18VLqpCSN0ztVXnIRhZnYn4Wb305JOIUlPeM5VqW/Dh3v3LGTvnh15b9E6npy2vAkidM41RJxaUgIwszckjQX+DtR82uhavYY0fJiVJc45uC9mxv2vL6R3pwKO37t7CqN0ztVXnBLGSVVvzGwlcCyhq1XnktqThg8l8a3DB9CzqIC3Pl4bqzl151zTqLGEIelcM/srcHYNB/6rKYvKZbTvHT2o7onqMLJXIdt2VHD/Gws568A+FBbkkpOki1nnXNOp7QisamywQw0v51KqoE02Fx8zmB0VldwwcTbzV21Od0jOtWpKVa9okvoADwHdgUrgbjO7Q9LPgAuA1dGk15rZM9E81wDfJjQ9cqmZPV/bOsaMGWNTpkxJSfyueamsNDaV7eR/3lrEeV/oT4f83W+me9ewzsUjaaqZjanvfHFaq71dUkdJuZJekrRG0rkxlr0T+KGZ7QUcAlwsaUQ07rdmNip6VSWLEcBZwN6EeyR3SvL2IxwQbooXFuTy1QP78t6iT/mwdNNuXcB617DOpVaci8LHm9lG4GRC0yBDgR/XNZOZrTSz96P3mwjtT/WqZZYJwKNmVmZmC4H5wEEx4nOtSHGHPI4dXsKqjWXc8dJHu4zzrmGdS604CaOq7H8S8IiZfVrflUjqD+wPvBMNukTSdEn3S6p6DLgXsDRhtmXUnmBcK3b4kK5cNnYIL89bxf2vL6S8otK7hnUuxeIkjKclzQXGAC9JKga2x12BpPbAP4EfRCWVu4BBhEYMVwK/rpo0yey7XVuQdKGkKZKmrF69OsksrrWQxDHDunHk0GIWrN7C1MXrvGtY51KozoRhZlcDhwJjzKwc2Eq4fFQnSbmEZPGwmT0eLa/UzCrMrBK4h88vOy0D+iTM3htYkSSeu81sjJmNKS4ujhOGa+EGd2vP0JL2LFyzhenL1nPpsUO8a1jnUiBWxXYzW2dmFdH7LWb2SV3zKDy8cR8wx8x+kzA88fHf04CZ0fungLMk5UkaAAwB3o23Ga61k8SZo3szokdHHnxrEafu38u7hnWukcXqca+BDgO+DsyQNC0adi3hQcBRhMtNi4DvApjZLEmPAbMJNawurkpSzsWVk53FVeOHs3HbDv5v7ip+9fw8rjlpr3SH5VyLUOtzGFEpobeZLa1xojTy5zBcbV6eW8r5D8TbP7xrWNeaNPQ5jFpLGGZmkp4ARjc0MOfSwcz4y5uLuWr8MMp2VnLpsUPIyvJW+Z3bE3HuYbwt6cCUR+JcI6rqGvY7Rwxk3IgS3ltU79rgzrlq4iSMYwhJ4+Po2YkZkqanOjDnGqp617B79yykqG0bXv9oTbpDcy6jxbnpfWLKo3CuESXrGnZY9w68OKeUQwZ29lZvnWugOM9hLCY8H3Fs9H5rnPmcS4fauoa9+JjBTJy+kh07vaqtcw0Rp/HBnwJXAddEg3KBv6YyKOcaIk7XsAf07cSj7y1p4sicaxnilBROA74EbAEwsxV4fxiuGYrTNWzfLm35ypg+TJ63qgkjc65liJMwdlh4WMMAJLWrY3rn0iJu17D5udnMWrGRpZ9ubaLInGsZ4iSMxyT9GSiSdAHwInBvasNyrv6+d/Sg2J0nXXTUINZu2cGWsp0pjsq5liPOTe9fAf8gNCI4DPiJmf0+1YE5l0rZWaJXUQG/r9anhnOuZnFuet9mZpPM7Mdm9iMzmyTptqYIzrlUKu6QxyXHDmbm8g3pDsW5jBDnktS4JMP82QzXInTIz2Xq4nW8v8SbQneuLjUmDEnfkzQDGBY94V31Wgj4k96uxfjGof3omJ/Dui070h2Kc81abSWMvwGnEPqpOCXhNdrMzm2C2JxrEpIo6ZjP7178kIrKmltvdq61qzFhmNkGM1sEXAd8Ej3lPQA4V1JR04TnXNPokJ/Ltw4fwJrNZekOxblmK849jH8CFZIGE3rQG0AofTjXovTr0o63F6z1h/qcq0GchFFpZjuB04HfmdnlQM2P0jqXwSaM6kW7vBzKdnpnj85VFydhlEs6G/gGMDEalryhHudagNF9O3Hzv+d40nCumjgJ43zgUOAmM1soaQAxGh+U1EfSy5LmSJol6bJoeGdJkyR9FP3tlDDPNZLmS5on6YSGbpRzeyIrS1xw5EAWrN6S7lCca1biPOk928wuNbNHos8LzezWGMveCfzQzPYCDgEuljQCuBp4ycyGAC9Fn4nGnQXsDYwH7pSU3ZCNcm5P9e7Ulm3lFTw5bXm6Q3Gu2YjzpPdCSQuqv+qaz8xWmtn70ftNwBygFzABeDCa7EHg1Oj9BOBRMyszs4XAfOCgem+Rc43kgL6dMAv9gzvn4vW4NybhfT7wZaBzfVYiqT+wP/AOUGJmKyEkFUndosl6AW8nzLYsGuZc2py6fy/+/MrHnHNwXzrk+60717rFuSS1NuG13Mx+BxwbdwWS2hOq5v7AzDbWNmmy1SdZ3oWSpkiasnr16rhhONdgZ47uzeR5vq85F+eS1AEJrzGSLiJmB0qScgnJ4mEzezwaXCqpRzS+B1BV6X0ZoSvYKr2BFdWXaWZ3m9kYMxtTXFwcJwzn9kiX9nkcNayYR9/1nvpc6xbnktSvE97vBBYBX6lrJoVebO4D5pjZbxJGPQWcB9wa/X0yYfjfJP0G6AkMAd6NEZ9zKdcxPxcDVm3aTrcO+ekOx7m0qDNhmNkxDVz2YcDXgRmSpkXDriUkisckfRtYQrgngpnNkvQYMJuQmC42M68I75qNsw/qy0tzStm3tyjukJfucJxrcqqpBoikK2qbsVqpIS3GjBljU6ZMSXcYrhXZUraTP7+6gCvGDU13KM41mKSpZjam7il3Vds9jA51vJxrddrl5XD5cUN4bubKdIfiXJOr8ZKUmf28KQNxLlNIorzCeGP+Gg4b3DXd4TjXZOLUknowsTlzSZ0k3Z/SqJxr5k7Zrye9igpYsnZrukNxrsnEaUtqXzNbX/XBzNYRHsJzrlXrWVTAva8vYMfOynSH4lyTiJMwsqo1ENiZeNVxnWvR2uRkcdX44XxYuom7Jn/spQ3X4sV9DuNNSf8gPHn9FeCmlEblXIZol5fDuq07mLF8PW8tWMuD5x9IeATJuZYnTtMgDwFnAKXAauB0M/ufVAfmXKY4YkgxV40fzrJ1W5k43WtPuZYr1qUlM5tNeKDOOZdEn05tGV7SgV88PZsjhxZTWOANFbqWJ849DOdcHbKyxM2n78PBAztz+3Nz0x2OcynhCcO5RlLUtg1fO6gfT05bwdTF69IdjnONLlbCkNRP0nHR+wJJ/qS3c0kcOrgLFxwxgKv/OZ3yCq9u61qWOA/uXQD8A/hzNKg38EQKY3Iuo106dghlOyu448WP0h2Kc40qzk3viwldpb4DYGYfJfSS55yrRhLXn7w3Fzw0hT+8PL/O6S8bO4TLvTFDlwHiJIwyM9tRVbdcUg5JesJzzgVmxv+8vZjvHjmQvl3a8rWD+6U7JOcaRZx7GK9IuhYokDQO+F/g6dSG5Vzmenr6Sko3bOdHJwwjNzvL72W4FiNOwria8MDeDOC7wDPAdakMyrlMtWFrOTdOnM3Np+9DbnYWXxnTh0ffW0rZTu8LzGW+OD3uVQL3RC/nXC1ue34u40aUMLrfZ82vcfTQYp7+YCVnju6dxsic23N1JgxJM9j9nsUGYApwo5mtTUVgzmWaqYs/5cXZpUy64qhdhvfp3JaCNtm8OX8NX/D+M1wGi3NJ6lng38DXotfTwKvAJ8ADNc0k6X5JqyTNTBj2M0nLJU2LXicljLtG0nxJ8ySd0MDtcS4tyisqufbxmVx/8oikzYJ0bZ/HKx+tZuuOnWmIzrnGEaeW1GFmdljC5xmS3jCzwySdW8t8DwB/AB6qNvy3ZvarxAGSRgBnAXsDPYEXJQ01M7/w6zLCPa8toHthPifv26PGaa48YThTF6/jwP6dvEVbl5HilDDaSzq46oOkg4D20ccaT5fM7FXg05hxTAAeNbMyM1sIzCc8++FcRhDihgkja00E2VkiS/DYlKVNGJlzjSdOCeM7wP2S2gMCNgLfkdQOuKUB67xE0jcI90B+GPXg1wt4O2GaZdEw5zLC944eFGu6Mf0707tTW1as30bPooIUR+Vc44rTH8Z7ZrYPMAoYZWb7mtm7ZrbFzB6r5/ruAgZFy1pJ6JwJQiLabdXJFiDpQklTJE1ZvXp1PVfvXPp1bd+GP73ysXft6jJOrP4wJH2RcH8hv6rIbWa/qO/KzKw0YZn3ABOjj8uAPgmT9gZW1LCMu4G7AcaMGeNPnLuMk5OdxUVHDWLtljJ6FHopw2WOOI0P/gn4KvBfhJLAl4EGtXUgKfGO4GlAVQ2qp4CzJOVJGgAMAd5tyDqcywQ9iwp4Z8GnvDl/TbpDcS62ODe9v2Bm3wDWmdnPgUPZtTSQlKRHgLeAYZKWSfo2cLukGZKmA8cAlwOY2SzgMUKvfs8BF3sNKdfSTRjVkwozzLyg7DJDnEtS26O/WyX1BNYCA+qayczOTjL4vlqmvwm4KUY8zrUIkjhkYBdufXYuV5843KvaumYvTgnjaUlFwC+B94FFwCMpjMm5ViM3O4uxe5WwbN22dIfiXJ1qLWFIygJeMrP1wD8lTQTyzWxDUwTnXGtw0IDOPDtjJTsrjQFd26U7HOdqVGsJI2p48NcJn8s8WTjX+I4bUcJLc0rrntC5NIpzSeoFSWfIL7A6lzK52Vl8+/ABPP7+snSH4lyN4tz0vgJoB1RI2kaoWmtm1jGlkTnXykiiotKYtWIDe/csTHc4zu0mzpPeHcwsy8xyzaxj9NmThXMpcObo3pRXGBu2lac7FOd2E+fBPUk6V9L10ec+UQOEzrlGJon+Xdpy58vz0x2Kc7uJcw/jTsLDeudEnzcDf0xZRM61ckVt2/BfY4cw75NN6Q7FuV3ESRgHm9nFRA/wRa3LtklpVM61cu3zcnhh1icsWbs13aE495k4CaNcUjZR67GSigFvZtO5FPvuUYPYsmMnOyv8cHPNQ5yE8XvgX0A3STcBrwM3pzQq5xxtcrLIzc7i3tcXpjsU54AY1WrN7GFJU4GxhCq1p5rZnJRH5pxjcLf2rN1cxs6KSnKy45zfOZc6cWpJ3QF0NrM/mtkfPFk417QOHtiF30z6kE3bvaqtS684pyzvA9dJmi/pl5LGpDoo59yuvnPEQGYs81Z5XHrFeXDvQTM7CTgI+BC4TdJHKY/MOfeZzu3aUFKYz9MfJO2I0rkmUZ+LooOB4UB/YG5KonHO1WhQcXs+3bKDikrvcMmlR5x7GFUlil8As4DRZnZKyiNzzu3mvC/0559Tl3nScGkRp/HBhcChZuadDzvXDIzqW8QT/1nOGaN7pzsU18rEqVb7J0mdovaj8hOGv5rSyJxzSQ0t6UBRQS7Tl61n395F6Q7HtSJxLkl9B3gVeB74efT3ZzHmu1/SKkkzE4Z1ljRJ0kfR304J466JamLNk3RCQzbGudaiuEMej7+/nP9+6SNvPsQ1mTg3vS8DDgQWm9kxwP7A6hjzPQCMrzbsakKXr0OAl6LPSBoBnAXsHc1zZ9QciXMuCUlcfeJwVm8u47onZ2Lm9zRc6sVJGNvNbDuApDwzmwsMq2um6JLVp9UGTwAejN4/CJyaMPzRqAvYhcB8QjVe51wN8nOzOWJIV+aXbmLi9JXpDse1AnFuei+TVAQ8AUyStA5oaGXwEjNbCWBmKyV1i4b3At5OXGc0bDeSLgQuBOjbt28Dw3CuZRg3ojs7dlby06dmceTQYgoLctMdkmvB4jy4d5qZrTeznwHXA/fxecmgsSTrLzxpGdvM7jazMWY2pri4uJHDcC7zjN2rhM7t2nDbs/54lEuterVmZmavmNlTZrajgesrldQDIPq7Khq+DOiTMF1vGl6Kca5Vyc/N5rYz9uWF2Z8wdfG6dIfjWrCmbv7yKeC86P15wJMJw8+SlCdpADAEeLeJY3MuY+3ftxNj9yrhir9Po9z7z3ApkrKEIekR4C1gmKRlkr4N3AqMi54cHxd9xsxmAY8Bs4HngIvNrCJVsTnXEt186kja5+dw72vef4ZLDWVydbwxY8bYlClT0h2Gc83Gy3NLOf+BeMfEZWOHcPm4oSmOyDVHkqaaWb1bHo9TS8o5lwHMjL+8uZiLjhrICXt3Z2SvQnK90yXXiHxvcq6FeHr6Sko3bOeHxw+jX5d23PzMHHbs9PsZrvF4wnCuBdiwtZwbJ87m5tP3ITc7i87t2vCtwwawcM0Wynb67UDXODxhONcC3Pb8XMaNKGF0v8+aZ6NP57Z0yM/hlmfmetJwjcIThnMZburiT3lxdilXjh++27ieRQV896iBLFyzhe3lnjTcnvGE4VwGK6+o5NrHZ3L9ySNqbBakR2EBHfNzucXvabg95AnDuQx2z2sL6F6Yz8n79qh1up5FBVw+bigzlm/wkoZrME8YzmUwIW6YMBIpWXNsuypq24aeRfn88vl53hy6axBPGM5lsO8dPYi+XdrGnr5HYQFXnzicF2aXeknD1ZsnDOdamdzsLPbrXcQ9ry5Idyguw3jCcK4V6l6Yz3+NHcLf3lnCth1e0nDxeMJwrhU7bq9uPDFtebrDcBnCE4ZzrVi3jvmcfVBf7njxIy9puDp5wnDOcfbBfXjto9Vee8rVyhOGc45uHfI5dng3fv70bLbu2JnucFwz5QnDOQdATnYWFx8zmDkrN3mVW5eUJwzn3GeKO+QxpKQ9Nz8zhy1lXtJwu/KE4ZzbRcf8XC4bO4TFa7d60nC7SEvCkLRI0gxJ0yRNiYZ1ljRJ0kfR3051Lcc5lxpd2udR0jGPW5+dy2ZPGi6SzhLGMWY2KqFf2auBl8xsCPBS9Nk5lyZd2udx+bihLF67xZOGA5rXJakJwIPR+weBU9MXinMOoHO7NvQoLOC2Z+f6jXCXtoRhwAuSpkq6MBpWYmYrAaK/3ZLNKOlCSVMkTVm9enUThetc69W5XRuuPnE4M5ZvYNP28nSH49IoXQnjMDM7ADgRuFjSkXFnNLO7zWyMmY0pLi5OXYTOuc+0y8thcHF7fjvpo88e7rtr8scsWbs1zZG5ppSWhGFmK6K/q4B/AQcBpZJ6AER/V6UjNudccp3ateG6L+7FczM/YeP2cgzjuidn+tPhrUiTJwxJ7SR1qHoPHA/MBJ4CzosmOw94sqljc87VLitLHDqoCw+9uYgLjhhI6YbtTJy+Mt1huSaSjhJGCfC6pA+Ad4F/m9lzwK3AOEkfAeOiz865ZqaobRsuOXYIj767hGu/uBc3TJzNhm1+b6M1UCYXJ8eMGWNTpkxJdxjOtUobtpbz7MyVTF++AQE3nbZPukNyMUmamvBIQ2zNqVqtcy6DFLbN5ayD+tIhL4cXZpUydfG6dIfkUswThnNuj3z/mMGctn8vrnl8OuUVlekOx6WQJwzn3B4pLMjlyvHD2FpWwQ8f+4ANW8upqEzdpW6vzps+OekOwDmX+XKys7j1jH059753eOqDFXVOf9nYIVw+bmiD1lVVnffB8w9EUoOW4RrGE4Zzbo+ZGXe/toCrxg/ne0cP+mz4ms1lfFi6iXZtcpg0u5Tj9y6hS/s8ehbmN3hdFxwxkCf/s4KJ01dyyn49GyN8F5MnDOfcHnt6+kpKN2znO0cM2GV41/Z5dG2fB8B+fYqorDQmzSnlf6cs5bT9e/Hqh6s5elg3+nRuG3tdudlZ3Hz6SL731/c5cmgxhQW5jbotrmZ+D8M5t0c2bC3nxomzufn0fcjNrv0nJStLnLB3d35w3FD6dWnHcSNKKNtZyRP/Wc6vnp/H2s1lLFi9uc6nx0f368xxI0q4/bm5jbkprg5ewnDO7ZHbnp/LuBEljO5X/y5sehQWADC4W3sAtpdX8PK81cxYvoGO+bms2rSd8SN7JC1FXHXCcMb99hVOP6B3g9ZdX3dN/pgv7tODvl3il4ZaGi9hOOcabOriT3lxdilXjh/eKMvLz83mzNG9mTCqF0cPK+bQgV0xM34z6UMefHMRy9dvY9m6UEOqsG0u1508gmsfn9Ek1XnT2XZWc6kZ5gnDOdcg5RWVXPv4TK4/eURK7iNIom+XthS1bcMV44Zy3hf6k5slXv1wDTOWbeD+1xeys6KCko553PvawkZff3XpbDuruTT06JeknHMNcs9rC+hemM/J+/ZosnV265jPOQf3BWBkr44s+XQrlQY/+t/p3BbjfsaeVOdN58325lIzzBOGc65BhLhhwsi0PQshib6d23L9k7N2qc67YVs5bbKzeG7WSsp3GkNK2vPy3FUcu1cJqzeV8dbHazmwfyeys1Tv2BNvtjdl21nNpWaYNz7onMtYT32wgj/+33wmXnp4nTW0AEo3bufjVZsZXNKev761mKHdO9A+L4cZyzZw+ujebNxWTu9OBXTIr/kHecPWcsb99hXuOnd0k9xsT3Ttv2Y0SkOPDW180BOGcy4jNeYP946dlRjGax+uoXTTdg4e0JmnPljJwQM6s728AjM4bHBXcrJFbnZWvRNVY2msbW5owvBLUs65jLQn1Xmra5MTfvSPG1Hy2bArxnUAYOP2chau3kLZzgrueW0xJR3zGNi1HdvKK7jlmbl8YVAXynZWcsTQrkxfuoHuhXl0zM9l4/ZyehQWsGNnJbk5WbRrk73Hl+8Sa4Y1dbICTxjOuQxUVZ130hVHpXxdHfNz2a9PEQCXjh3y2fCbThvJ1+97l/vfqLuG1pFDunLLGfvy9/eWsnfPjuRmi2lLN3DGAb3494yVZCk80PjEf5Yzqm8RZeUVzFm5ibMP6suj7y2hqCCXI4YW8/pHaxjVp5C2ednc+9rCXZphaQp+Sco5l1HKKyo5+fevc8mxg9NWY8jMOO8v73HowC5N8qNdUWlUVBqbtpfTJieLVz5czSV/+0+seZPVDPNLUs65ViEd1Xmrq6ntrFTJzhLZWaJL+zzMjMemLNutocem0Owe3JM0XtI8SfMlXZ3ueJxzzUu6q/PWp+2sVGjqZJWoWZUwJGUDfwTGAcuA9yQ9ZWaz0xuZc665aOqz6uoa82Z7fVUlq7vOHZ2WZNXcShgHAfPNbIGZ7QAeBSakOSbnnAMav+2s+kpnsoLmlzB6AUsTPi+Lhn1G0oWSpkiasnr16iYNzjnXeqW67ay6pDtZQfNLGMkuSu5SjcvM7jazMWY2pri4uInCcs61dum82Z7uZFWluSWMZUCfhM+9gbo7CHbOuRRL58325lAzDJrZTW/gPWCIpAHAcuAs4Jz0huScc+m92Z7ummFVmlXCMLOdki4BngeygfvNbFaaw3LOubRKd82wKhn9pLek1cDidMfRAF2BNekOoon5NrcOrW2bM3V7+5lZvW8CZ3TCyFSSpjTksfxM5tvcOrS2bW5t29vcbno755xrpjxhOOeci8UTRnrcne4A0sC3uXVobdvcqrbX72E455yLxUsYzjnnYvGE4ZxzLhZPGCkUp28PSUdLmiZplqRXmjrGxlTX9koqlPS0pA+i7T0/HXE2Jkn3S1olaWYN4yXp99F3Ml3SAU0dY2OLsc1fi7Z1uqQ3Je3X1DE2trq2OWG6AyVVSDqzqWJrSp4wUiShb48TgRHA2ZJGVJumCLgT+JKZ7Q18uanjbCxxthe4GJhtZvsBRwO/ltSmSQNtfA8A42sZfyIwJHpdCNzVBDGl2gPUvs0LgaPMbF/gBlrGjeEHqH2bq46B2wgtVbRInjBSJ07fHucAj5vZEgAzW9XEMTamONtrQAeFBnHaA58CO5s2zMZlZq8StqMmE4CHLHgbKJKU3hbk9lBd22xmb5rZuujj24RGRDNajP8zwH8B/wQy+TiulSeM1Kmzbw9gKNBJ0mRJUyV9o8mia3xxtvcPwF6EFohnAJeZWWXThJc2cb6XluzbwLPpDiLVJPUCTgP+lO5YUqlZNT7YwtTZtwfh+x8NjAUKgLckvW1mH6Y6uBSIs70nANOAY4FBwCRJr5nZxhTHlk5xvpcWSdIxhIRxeLpjaQK/A64ys4p0tyibSp4wUidO3x7LgDVmtgXYIulVYD8gExNGnO09H7jVwsM/8yUtBIYD7zZNiGnRKvt4kbQvcC9wopmtTXc8TWAM8GiULLoCJ0naaWZPpDWqRuaXpFLns749ohu7ZwFPVZvmSeAISTmS2gIHA3OaOM7GEmd7lxBKU0gqAYYBC5o0yqb3FPCNqLbUIcAGM1uZ7qBSSVJf4HHg6xlaWq43MxtgZv3NrD/wD+D7LS1ZgJcwUqamvj0kXRSN/5OZzZH0HDAdqATuNbNaq+01V3G2l1Bj5gFJMwiXaq4ys0xsGvozkh4h1PjqKmkZ8FMgFz7b5meAk4D5wFZCKSujxdjmnwBdgDujM+6dmd6ia4xtbhW8aRDnnHOx+CUp55xzsXjCcM45F4snDOecc7F4wnDOOReLJwznnHOxeMJwzjkXiycM55xzsXjCcC2CpCJJ3093HHVJdZySLoj6V5kmqTLh/W9StU7XeviDe65FkNQfmGhmI5tBLCIcW7u1xNuQOGtbXi3z9ALeNLN+cedxri5ewnAtxa3AoOhs+peSzpX0bvT5z1HnNkjqL2mupHslzZT0sKTjJL0h6SNJByVM82DUa9w/ora+iJax27KjeeZIuhN4H+gj6Ymo2fpZki6sIc7+ib24SfqRpJ/VsLyk21SDkYQm5J1rNJ4wXEtxNfCxmY0C7ge+ChwWfa4AvpYw7WDgDmBfQmu55xCa4P4RcG00zTDg7qjXuI3A9wEk7VXLsocROkva38wWA98ys9GElkwvldQlMU4z+3Ed2/TZ8oC2dWxTdfsAGdkumWu+vPFB1xKNJfQz8l7U+F0Bu/aCttDMZgBImgW8ZGYWNYrYP5pmqZm9Eb3/K3Ap8Ktalv0qsDjqVa/KpZJOi973IXTT+kk9tiNxeXVtU3UjgUn1WJdzdfKE4VoiAQ+a2TU1jC9LeF+Z8LmSz4+J6jf3qj4nXXZ0b2JLwuejgeOAQ81sq6TJQH6SWHaya0k/cZotCe/r2qbq9gF+G3Na52LxS1KupdgEdIjevwScKakbgKTOkup787evpEOj92cDr9dz2YXAuihZDAcOSRInQCnQTVIXSXnAyTXEE3ubJGURSjNz42yoc3F5wnAtQtSr2xvRDeTzgeuAFyRNJ1ya6VHPRc4Bzovm7wzcFa1ndsxlPwfkRNPcALxdPU5JvzSzcuAXwDvARGr4ka/HeiHco1lmZmU1jHeuQbxarXPVNKcqus41J17CcM45F4uXMJxzzsXiJQznnHOxeMJwzjkXiycM55xzsXjCcM45F4snDOecc7F4wnDOOReLJwznnHOx/H8Dct9Lorn/yQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -566,7 +566,9 @@ "source": [ "# Plotting\n", "with h5py.File(\"xy_data.hdf5\", \"r\") as f:\n", - " plt.scatter(temperatures, f[\"cluster-size\"], marker='4', s=100)\n", + " plt.figure()\n", + " # My girlfriend told me to use this marker.\n", + " plt.plot(temperatures, np.array(f[\"cluster-size\"]), marker=\"4\", markersize=20, linestyle=\"dashed\", linewidth=.5)\n", " plt.xlabel(\"temperature $T$\")\n", " plt.ylabel(\"average cluster size\")\n", " plt.title(\"Average cluster size for the $w \\\\times w$ XY model \\n\"\n", @@ -741,7 +743,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "id": "1015f826", "metadata": { "deletable": false, @@ -757,11 +759,45 @@ "task": false } }, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEtCAYAAAAWZydGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABANUlEQVR4nO3deXxU1fnH8c83ISEgS5B9CauIrKKgYOuCRSpordraumstFmnr0lq3/n6ttbXW2vZXq1WLuFRrXWvdanErFlFZBBRRQGQTCCggO4SQkDy/P+5NHMJkkrlMMlme9+s1r8xd5tznTu7MM/ece8+RmeGcc85VJiPdATjnnKvbPFE455xLyBOFc865hDxROOecS8gThXPOuYQ8UTjnnEvIE4VzzrmEPFE455xLyBOFa/QkLZQ0qh6V20/Se5J2SLoy1eWH2/hE0kk1UbZLnqSHJP26muum/H/niSIUvrlFktpVmD9fkknqmabQXIpV/CCZ2UAzm5bKMlNVbiWuA6aZWUszu/NAC2uMSUFSi3C/z4uZ11LSaklnSXpU0oMVXnOCpE2SOtd+xOnliWJfK4FzyyYkDQaapS+c9JHUJN0xuEr1ABZGeWF9+L8qcESc+YdLykzFNsxsJzABuENS+3D274C5ZvY0cCVwiqQx4bZzgPuAn5jZp6mIoV4xM38E/V19AvwMmBMz7w/A/wIG9AzndQH+CWwkSCxXxqx/A7Ac2AEsAs6Ms41rgAXANuBJIKeSeK4H1oZlLQFGh/OPAN4N5z8JPAH8OuZ1BhwSM/1Q2fJqxnd9GN8eoEmi/a3kPbw2fP0u4AGgI/BSuM3/AG2q834BRwLvhcv+Ee7rr6vzXlbxP3oEKAV2AzsJfp1/ApwULj87nF/22EPw673SeOOVGRNjWbn9gWnAVoIv+a9HOTaA14ESoDDc1qHVLHuf/2ui96M68SR5XOwE8sLnEwmO0Y7h9LXA/RXW7xWWOzZm3nHA58CgAy2/wmsfAh4HRgGbgM4xy74V7ttBwK3ASyk89qv6n1X6Oa/qvSfmuEvZ92MqC6vPj7I3l+BLuT+QCawh+PVmQE+CM7B5wI1ANtAbWAGcHHNgdQnXOzs8YDpX2MY74ToHA4uBiXFi6Rduu0s43RPoE25zFfBjIAs4Cyim+omiOvHNB/IIzqQS7m8l7+Gs8APSFdgQHuxHAE0JvuR+UeGDuF88Mft5Vbif3wCK2D9R7PdeVifmih+kyj5YQKuw3Muqev/ilcEXx1QWsAz4nzCmrxB8AfRL9tgI150GXBo+r27Z5f/Xyo79OPPixhPhuFgDDAAEfAAsBQ4Lp5cBh8d5zfEEX9wnAkdTIXEcaPkxr20DfEqQhC6Js/xp4IUwlu5VfH9U69iv6n9Ggs95dd77eP/PA3141dP+HgEuAsYAHxH8qi9zFNDezH5lZkVmtoLgdPQcADP7h5mtM7NSM3uS4IA9ukL5d4brbAb+BQyNE0MJwcE1QFKWmX1iZsuBkQQHzp/MrNiCU+Q51d2xJOJbY2a7q9rfSvzZzNab2VrgTWC2mb1nZnuAZwk+OFXFM5LgbObOcD+fIfjSqijeexkl5v1IygAeIzibuLeKeKsyEmgB/DaM6XXgRWKqORPsTyrLLvu/Vldl8ST7Hm8NY/wqwRfkh0AuMBZYa2bvV3yBmU0HziP4on4RmGBmL6eq/JjtbCH4Rd8ceCbOKj8k+CL/lZmtrqycUHWP/ar+Z4k+5yk5vpNV5+sr0+ARYDrB6e/fKizrAXSRtDVmXibBQYGki4CrCc4AIDgY9mkcBz6LeV5A8IttH2a2TNKPgJuAgZJeCcvtQnDgx/YNv6p6u1Xt+NbEPE+4v5VYH/N8d5zpFtWIpxn772dsXGXivZdRYo7nFqAlQV11VfFWpQuwxsxKY+atIvjlGavKY+MAyo73/lWlsniSfY+3ELxXPwJuAy4h+CX/AyBRY/xqYC/Br+hPEqwXtXwkXUDw//xP+NqJscvNbL2kz6lem1B1j/2q/meJPuepOr6T4omiAjNbJWklcAowvsLiNcBKM+tb8XWSehBk9tHATDMrkTSf4PQ3ShyPAY9JagXcS3AQ3w90laSYg6g7Qb15mQKCX0dlOgH5ScRX8cs57v4eqCri+ZT99zOPffezMtWJOeEgLJLOIfh1d5SZFVcj3qrKXAfkScqI+XLoDnxc1c5UQ3XLThRfsoPSJHtcbAWGE1TTTZP0TeBwYBDwXLwXSOoDvEbQtrIDmCLpJDOL94WddPnhNjoAtwPfJqg9WCjpsfBspiZV9T+Ld/yXfc5r7DOZiFc9xTce+IqZ7aow/x1gu6TrJTWTlClpkKSjCBq8jKAuFUmXEByoSQuvk/+KpKYEjZa7CaqjZhL8wrpSUhNJ32D/qo/5wHlhbGOBE8L5UeJLtL8HKlE8Mwn29/JwP0+nelU81Y15PUHd7n7Cq23+DJxhZhurGW/CMoHZBO0Z10nKCu+tOI2ggfJApaLsRLHHk+xxsYWgvr3s1/12gvane82spOLKkroAU4FbzOwhM/snQcP6q5LixZlU+THuAp4zs/9acCXTdcB94eeuJlX1P0v0Oa/Jz2SlPFHEYWbLzWxunPklBP/QoQRXG3xO8Cu/tZktAv6P4J+8HhgMvB0xhKbAb8PyPwM6AP9jZkUEDbvfIfhwnM3+9apXhTFuBc4n/EUVJb5E+xtxv2LLrjSemP0cH+7HBQR1uHuqUW51Yr4V+JmkrZKuqVDE6QTVFm9J2hk+XqrG+1dpmeH+fB0YF8ZzD3CRmX1U1f5UY39TUXai9yPeNpM9LrYQ1F48Fk5vJ2hDuK+S9TcRXIb6l5htPkpQlbQhBeUj6QzgWIIrlcq2cT+QT9BQXGOq+p8l+pzX5GcyEe1bDebqG0kPAflm9rN0x1KTJM0GJpnZX9Mdi3ONjZ9RuDopvAu2U3jqfTEwBKjsqhfnXA3yxmxXV/UDniK4UmQ5cJY1xjtinasDvOrJOedcQl715JxzLiFPFM455xLyROGccy4hTxSNgKRbwy5BnGvUJF0p6bfpjqO+8UTRwIV97V9E0A1IOrZ/raQPFYzGtlLStRWWfyJpd8zNba9G3M5TksZIairps6pf4apD0ihJ+bW8zVMlvRXeAPiZpPsktYxZ3lTSg5K2h8uvrvD6oZLmSSoI/w6NWTwZuCDsvsNVkyeKhu87wJQkew1NJREkqjYEvXleHvalFOs0M2sRPr4acTvDCLp1HkLQe2i9pXowuFANa03QpXYXgi7/uwG/j1l+E9CXoIO8Ewm6whgLICkbeB74O8Ex9zDwfDgfMyskGCPiotrYkYbCE0XDNw54o2xC0hthp2lIOlbBMK+nhNMnKejoLmXM7Hdm9q6Z7TWzJQQf4i+nchuS2hBc6r2JoHO4dyssnybp15JmhGct/5LUVsFwl9slzVHMULeSDpP0mqTNkpZI+nbMslMVjFe9XdIaSTfFLMuR9HcFw2VuDcvtGC7bZ7hRSTdJ+nv4vGf4fxgvaTXB2AVI+q6kxZK2SHpFQceEZa83ST+QtDQ8W7tZUh9JM8PYnir7cgzX/5qCYX23hu/DkJhln0i6RtICSdskPRnuy0EEX6pdYs74ukg6WtLccDvrJf3xQP+HsczsMTN72cwKwm7A72PfY+Yi4GYz22Jmi8Pl3wmXjSK4P+xPZrbHgqFiRdBVeJlpwKmpjLmh80TR8A0mGIypzBsEHyYIBodZwRcdBx5PTFKJJem88Eumskf3qgKRJIKRyir2APqopI2SXpV0eHV3TNJoBd0trwG6hc/vAH4YxnRCzOrnABcSdOXch6DPpr/yxaA8vwjLPIig19LHCPrYOhe4R9LAsJxdBF9UuQRfNt9X0G8QwMUEv4bzgLYEXVYncyZ3AsEv6JPDMv+HoM+f9gTdSD9eYf2xBGdSIwk6tJtM0L9XHkGHheeG+3Qk8CBwWRjXvcAL2rfzu2+H5fUiOCv7Ttgp5jhgXcwZ3zqC9/gOM2sVvpdPxdsZSd2rOGbOi/e6OI4nPGbCHwVdgNgxJt4Hyv4/A4EFFbroXhCzHIL/d7WPM+eJojHIJeimucwb7JsYbo2ZPoFKEkX4Ky83waOqQV0gqDLIIPiCLnM+wXgAPYD/Aq9Iyq3OjpnZVDPLJej48FsESeATgoFdcs0sdl/+Gnb2uI3gV/JyM/uPme0lGGq1bFCZrwGfmNlfw7OgdwmGnTwr3OY0M/vAgsGLFhB8eZe9f8UEX8SHmFmJmc0zs+3V2ZfQTWa2K6wmvAy41cwWhzH+Bhgae1YB3GZm28Outz8EXjWzFTH7WLZP3yPoSXV2GNfDBB0sjowpK5lBk4qBQyS1M7OdZjYr3kpmtrqKY+axeK+LpWDM6ov5oqO+sjEdtsWsto1g7JCy5bHLKi6H4PNQo53oNTSeKBq+Lez7IZkJHBpWiQwlGJwpT1I7gq6Ma6QvfkmXE/wSP9WCEb8AMLO3zWx3WM1wK0FvscdVs8z88CziXIK66A0ECWddnOqQ6g4q0wMYEfvLlyCZdQq3OULSf8MzoG0EZw1lgxc9ArwCPCFpnaTfScqqzr6EKg4adUdMDJsJqlBiByRKZp9+UmGf8th3YKSKgxS1oHLjCcbq/iisXvtaVTsWhaSRBGd2Z5lZ2VgNO8O/rWJWbcUXP4Z2VlhWcTkEn4eKycQl4Imi4VtA8KEGwMwKCMbcvQr4MOzSeAbByG3LzezzeIVIOj+mnjreo9KqJ0nfBW4ARptZVVfQGNUc7MnMuhFUl/wnPLOYDPww/LV6dcIXV24N8EaFX74tzOz74fLHCMZQzjOz1sCksngtGLbyl2Y2APgSwdlJWaPpLvYfUGq/XaoQx2UV4mhmZjMi7tMtFcpqbmYVq7Li2a+PHzNbambnElTN3QY8HVbZ7SOsekp0zJxf2UYVjAvyAvBdM5sas+0tBAP7xFYdHc4X1ZkLgSFhNWeZIexb3dmffauuXBU8UTR8U/iiaqTMG8DlfFHNNK3C9H7M7NGYeup4j7hVT+GXwW+AMRaM7xu7rLukL0vKDhtPryX4df52uLyskbdngv0ru9oJ4Ehgv3FEkvQiwRnXhQoGlcmSdJSk/uHylsBmMyuUdDTBuM5l+3OipMGSMgnGRCgmGIAJggGlzgnLG05YlZXAJOCnZW0jklpL+lbEfboPmBieDUnSQQoa5VtW+crgLKWtpPKqGkkXSGpvwehsW8PZ+w0QFFY9JTpmHo23QUmDCHoKvsLM/hVnlb8RjJ/RRtJhBFVrD4XLpoWxXKngMtrLw/mvx7z+BIKqOVdNnigavr8Bp0hqFjPvDYIvvOmVTKfSrwnq7efE/JKcFC5rCfyFoHpsLcHZwTgLrl6CoHpkVbisMsOAd8NfkIdRvbGNK2VmO4CvEjR+ryOokrmNYDApCAbP+ZWkHQT15rENuZ2ApwmSxGKC9/Xv4bKfEzT8bgF+yReD7FQWx7Phdp+QtJ2gDWJcxH2aS/Blele4/WV8cZVQVa/9iKAdZkVYbdWF4P+0UNJOgobtcyy47DRVfkLQgP9AzDET+3/9BUGPwqsI3uPfm9nLYbxFwBkEZ3Jbge8SjFZYBMGVaQTDHD+cwngbPO89thGQ9Btgg5n9Kd2xJEPSz4CNZpaWmwVdwyPpCoJqw+vSHUt94onCOedcQl715JxzLiFPFM455xLyROGccy6hBtf5WLt27axnz57pDsM55+qVefPmfW5m7eMta3CJomfPnsyde6CX0jvnXOMiaVVly7zqyTnnXEKeKJxzziXkicI551xCniicc84l5InCOedcQp4onHPOJeSJwjnnXEKeKJxzziXU4G64S9btr33MHVOXVrr8qtF9+fGYQytd7pxzDV2D62Z8+PDhFvXO7LPvnQnAk5cdk8qQnHOuzpM0z8yGx1vmVU/OOecS8kThnHMuIU8UzjnnEvJEESopNbYUFLF2y26mLl5PSWnDartxzrmoPFEQJIkLH5jNsg07yd+6mysef48LH5jtycI55/BEAcC0JRuYv2YrZXmhoKiE+Wu2Mm3JhvQG5pxzdYAnCmDhuu3sLirZZ97uohIWrduepoicc67u8EQBDOzSimbZmfvMMyB/624Ki0viv8g55xoJTxTAqH4dGJqXS4aC6aZNMmjdrAlPzlnDuffNSm9wzjmXZo2+Cw+AzAzxyPgRjLtjOgV7Svjl6QMZ1a8Dcz/ZTEFYJVVYXMI/5q7hW8PzyMnKrKJE55xrODxRhDIzRJvm2bRpDqP7dwRgRO+25ctf/2gDP39+IXe+voyJJ/Th/BHdPWE45xqFtFY9SRoraYmkZZJuiLP8Wknzw8eHkkokHZyOWE8Z3JknJ4ykb4cW3PziIo773X954K2Vfgmtc67BS9sZhaRM4G5gDJAPzJH0gpktKlvHzH4P/D5c/zTgx2a2OR3xQnCG8VjvtsxasYk7/rOUf72/ju9+uScApaVGRlkjRzV4r7XOufoinVVPRwPLzGwFgKQngNOBRZWsfy7weC3FltDI3m0ZOaEtOwqLkcSmnXv45l9mcOExPatdJfXjMYeWJwLvtdY5V5els+qpK7AmZjo/nLcfSc2BscA/K1k+QdJcSXM3btyY8kAr0zInC4DthXvp3LrZPlVSflmtc66hSOcZRbx6msoq/E8D3q6s2snMJgOTIRiPIpkg4lUB9bzh3+XPq1MF1KvdQTw+YWR5ldTNLy5i8vTlvHb1CbQKk4lzztVX6UwU+UBezHQ3YF0l655DDVU7xVYBHaiyKqnZKzbxzsrN5Univ0s2cEzvtn6VlHOuXkqYKCSVmFlNfbvNAfpK6gWsJUgG58WJoTVwAnBBDcWRciN6ty2/tHbN5gK++9Ac2rVoyvdP6MN5flmtc66eqaqNovqX8STJzPYClwOvAIuBp8xsoaSJkibGrHom8KqZ7aqpWGpS3sHNefx7IzmkfQt+5W0Yzrl6qKqqpxq9ScDMpgBTKsybVGH6IeChmoyjppVVSZW1Yfzu5Y84bUhncrIyy8fBKNhTwtTF6xnVrwOZSVxm65xzNU1mleeC2KonSWOAbwN3m9l8SRPCRuQ6Zfjw4TZ37tx0h5HQms0F5B3cnL0lpYz4zVQ27yrCgObZmQzNy+WR8SM8WTjnapWkeWY2PN6yZC6P/QFwLXCBpK8AQ1MQW6OUd3BzAF7+8DO2FBSVn7b5OBjOuboomUSx0cy2mtk1wFeBo2oopkZjxee7qHhCV1BUwozln6cnIOeciyOZRFF+c4GZ3QD8LfXhNC7xxsHIEAzr0QaA4pLSdITlnHP7qHaiMLPnK8yaFHdFV20Vx8Fonp3JyN5tOXlgZ4r2lnLKHW/y6xcXsa2gOL2BOucatUhdeEi6H1gvaY2k2ZLuk3RFimNr8MrGwTikQwu65Tbjz+ceUd6QXbi3hGE92vDA2ys54Q//5aG3V/oZhnMuLaL29XQc0NHM8oBvAM8CB6UsqkakbByMrm2aMbp/x/KrnVrlZPHbbw7h31ccx4DOrbjpX4s4+fbpfLptd5ojds41NlG78JgFtAE2mNlagjurpyR+iYtiQJdWPHrpCF7/aAPPvreWji1zANheWOz9SDnnakXUM4rJwBuSrpF0XNjNhqshkhjdvyN3nXckGRliy64iTvjdf7n+6QVs2FGY7vCccw1c1DOKvxNc9dSE4P6KIZJyzKxPyiJzlcrMFN88shsPz/yEfy1Yxw9G9eHS43pXuw8pHzTJOZeMat+ZXWH+m2Z2XIV5Tc1sTw3EmJT6cGd2RVEHLlr5+S5++9JiXlm4ni6tc5hy1XHkNs+ulW075xqWRHdmRz2jmC/pKjO7o2xGXUgS9UmqxsG498LhzFqxiTeXbixPEvlbCujWpnnqg3bONUpRE0VH4CRJ1wPvAu8D883sHymLrIFL+TgYYbfmyzbsZOyfpnPywE5cP/Ywurf1hOGcOzCRuhk3s2+bWX+gF3Aj8DHBGNguzbrk5nD5Vw7h9Y82cNIf3+DWKYvZXug37Dnnokt4RmFmCRNJWN30bvhwdUDz7Cb86KRDOeeo7vzh1SVMfnMFz89fx7RrR/mASc65SNI5FKqrQZ1a5/CHbx3Od77UkwX528qTxPtrtnJ4Xm56g3PO1StR76Nw9cSgrq05b0R3AGYu38Tpd7/NRQ++w5LPdpQPmrR2y26mLl5PSWmNjlPlnKunEl4eW+Mbl8YCdwCZwP1m9ts464wC/gRkAZ+b2QmJyqyPl8fWlqK9pfxt5ifcOXUpOwr30qFVUzbu2EOp+aBJzjV2Kb88VlJT4JtAz9gyzOxXSZSRCdwNjAHygTmSXjCzRTHr5AL3AGPNbLWkDlHidYHsJhlcelxvvnlkN675x/tM/eiLAZJiB00a3b9jGqN0ztU1UauengdOB/YCu2IeyTgaWGZmK8ysCHgiLDPWecAzZrYawMx86LcUaHNQNofn5e53SVtBUQkfrt2Wlpicc3VX1MbsbmY29gC33RVYEzOdD4yosM6hQJakaUBL4A4z22/AJEkTgAkA3bt3P8CwGoeyQZMKikr2mf/EnDUM6tqarxzWAcmroJxz0c8oZkgafIDbjvctVLHBpAkwDDgVOBn4uaT97lIzs8lmNtzMhrdv3/4Aw2oc4g2a1K9jS5o2yWD8w3M5e/Is3lu9Jb1BOufqhKiJ4lhgnqQlkhZI+kDSgiTLyAfyYqa7AevirPOyme0ys8+B6cDhEWN2MeINmjTlquN47eoTuPmMQazYuJPrnl5AqV8J5VyjF7XqaVwKtj0H6CupF8F4FucQtEnEeh64S1ITIJugaur2FGzb8cWgSW2aU96AnYm4cGQPzjyiK59t201Ghti5Zy9/nrqU7x3fm3YtmqY5audcbYuUKMxslaTDCUa6A3jTzN5Psoy9ki4HXiG4PPZBM1soaWK4fJKZLZb0MrAAKCW4hPbDKDG75LRo2oRDOrQEYNbyTdz/1kr+PmsVE0/ow/jjetE8O/lDx7s3d65+inQfhaSrgO8Bz4SzzgQmm9mfUxhbJH4fRXKq28348o07+d3LH/HKwvV0aNmUH510KOcenRe5wdu7N3eubkl0H0XUNorxwAgzu9HMbgRGEiQO10D1ad+Cey8cztMTjyHv4Oa89OGnflWUc41E1DYKAbHXVZZQSU+zrmEZ3vNgnp54DDv27AWCsS+u/+cCrh7Tj2E92qQ5OudcTYiaKP4KzJb0bDh9BvBASiJydZ4kWuVkAbB6UwEfr9/JN/8yg3GDOnHtyf3o3b5FmiN0zqVS1MbsP0p6A/gywZnEJWb2XkojczUmFaPrlfnSIe2Yds0o7n9zJZOnL+fVReu5cGQPfnHaAK+acq6BiNzNuJnNA+alMBZXS1I5uh7AQU2bcNVJfTlvRHfunLqUUrPyJLFnbwlNm/g4GM7VZ0klCklvmdmxknaw713UAszMWqU0OlevtG/ZlJvPGETZlXTzVm1h4t/ncdXovpx9VB5Zmd6rvXP1UVKfXDM7Nvzb0sxaxTxaepJwZcrOJpplZdKr7UH87LkPOfn26bz84WeYmY+D4Vw9E/U+itvM7Pqq5qWD30dRt5gZ/1m8gdte/ohlG3ZywqHtKC4xZq3Y5ONgOFeH1MR9FGPizEtFtx6ugZHEmAEdefmq4/jtNwaTd3Bz5q/ZStlJROw4GM65uimpRCHp+5I+APqFnQGWPVYCH9RMiK4haJKZwTlHd6dDyxx2V+jafHdRCYvWbU9TZM65qiR71dNjwEvArcANMfN3mNnmlEXlGqx442BI0KGVdzboXF2VbGP2NjP7xMzOBbYDHYEewCBJx9dEgK5hqTgORnZmBgJ+9a9FvPB+xV7mnXN1QaQ2CkmXEowN8Qrwy/DvTakLyzVUFcfB+MsFRzL1J6MY3K01OU388lnn6qKoN9xdBRwFzDKzEyUdRpAwnKtSvHEwHv/eyPLLah+e8QkdWjZl3ODO6QzTOReK+hOu0MwKASQ1NbOPgH6pC8s1NmVJoqTUeOH9dXz/0Xf50RPvsa2gOM2ROeeiJop8SbnAc8Brkp5n/2FMnUtaZoZ4YsJIfnRSX/614FNO/tN0pn+8Md1hOdeoRUoUZnammW01s5uAnxP0HHt6KgNzjVdWZgY/OulQnv3Bl2iR04RLHppD/paCdIflXKMVqY1C0lTg/8xsipm9Ec6bDExIspyxwB0EQ6Heb2a/rbB8FMG42SvDWc+Y2a+ixOzqnyHdcnnximOZuWIT3do0B2DN5gLyDm6eVDk+BKtzByZqY3Yv4HpJR5lZWSN23Fu/KyMpE7ib4C7vfGCOpBfMbFGFVd80s69FjNPVczlZmZzYrwMAM5Z9zoUPvsNlx/fmqpP6VrtX2tjecn0IVueSF7WNYiswGugo6V+SWkco42hgmZmtMLMi4Am8+solMLhba846shv3TFvO6Xe9zeJP/W5u52pD5KFQzWwv8ANJ3wHeApIdB7MrsCZmOh8YEWe9YyS9T9BYfo2ZLdwvGGkCYbVX9+7dkwzD1ZYDHTCpZU4Wt501hK8O7Mj1//yAr9/1FjeM68/4Y3vVWMzOueiJYlLZEzN7KOz/6YdJlhGvq9CKXdm+C/Qws52STiG4yqrvfi8ymwxMhqD32CTjcLUkVQMmje7fkVd/3IafPfcB2X6TnnM1LupQqPdWmJ4HfDfJYvKBvJjpblS4xNbMtsc8nyLpHkntzOzzJLflGpiDD8rm7vOOLJ9+fv5athfu5YIR3X0IVudSLNneY98K/+6QtD3msUNSshXGc4C+knpJygbOAV6osL1OCj/1ko4O492U5HZcAyWpPCm8umg9P3/uQy568B0+3bY7zZE517CkbYS7sI3jcoJ+ohYDT5nZQkkTJU0MVzsL+DBso7gTOMeijLTkGry7zj2Cm88YxNxPtnDy7dN57r21+KHiXGpE7RTwturMq0p4H8ahZtbHzG4J500ys0nh87vMbKCZHW5mI81sRpR4XcMniQtH9uClq46jb8eW/OjJ+cxfszXdYTnXIPgId65B6dnuIJ667BgeuHg4R3QPLsRbvcnv6nbuQPgId67BycxQea+0H6/fweg/TuO6p99na0ERWwqKWLtlN1MXr6ek1KumnKsOJVOPG95Y14Y6PMLd8OHDbe7cuekOw9URe/aWcMd/lvKXacvJysyguKQUA5pnZzI0L5dHxo8gM8OvknJO0jwzi9vDRqQR7oDzgeOAi81sFdAivCrJuTqlaZNMrht7GDeMO6w8SQAUFJUwf81Wpi3ZkNb4nKsPorZR3A0cA5wbTu8I5zlXJ+3ZW7rfvN1FJSxa592AOFeVqIlihJn9ECgEMLMtQHbKonIuxQZ2aUWz7H07EczJyqR/55Zpisi5+iNqoigOe381AEntgf1/sjlXR4zq14GhebmUNUc0y8qgqKSURd6xoHNVipoo7gSeJeg99haCTgF/k7KonEuxzAzxyPgRHNKhBd1ym3HnuUfw9SGd+eNrS3lk1qp0h+dcnRa1r6dHJc0j6Goc4AwzW5y6sJxLvcwM0aZ5Nm2aw5gBnRjVrwM79uzlxuc/pHWzLL5+eJd0h+hcnRT1zuymwJFAa6At8C1JN6YyMOdqWlZmBneddyRH9TyYq5+cz7xVdeIKb+fqnKjdjD8PbAPmAXtSF45ztSsnK5P7Lx7O3f9dxsAuUcbfcq7hi5ooupnZ2JRG4lyatMrJ4qfj+gOwraCYjTv3cEiHFmmOyrm6I2pj9gxJg1MaiXN1wJVPvMd5983y/qGcixE1URwLzJO0JOzr6QNJC1IZmHPp8LNT+1NUUsqFD85mw47CdIfjXJ0QNVGMIxiS9KvAacDXwr/O1Wt9O7bkoUuOZuOOPVz0wDts212c7pCcS7tIicLMVsV7pDo459JhaF4uky8czvKNO/n5cx+mOxzn0i5SY7akq+PM3gbMM7P5BxSRc3XAsX3bMfmi4fTvlNTAjc41SFGrnoYDE4Gu4WMCMAq4T9J11S1E0tiwnWOZpBsSrHeUpBJJZ0WM17mkndivA51a51BSajw5ZzWlPn6Fa6SiXh7bFjjSzHYCSPoF8DRwPMG9Fb+rqoCwr6i7CUbLywfmSHrBzBbFWe82grG1nUva7a99zB1Tl+4zr+cN/y5/ftXovvx4zKGVvv61RZ9x/T8/YOG67fzy6wORfPwK17hETRTdgaKY6WKgh5ntllTdG/COBpaZ2QoASU8ApwOLKqx3BfBP4KiIsbpG7sdjDk2YCKpy8sBOTDi+N5OnryC3eTZXH0BZztVHURPFY8AsSc+H06cBj0s6iP2/6CvTFVgTM50PjIhdQVJX4EzgKyRIFJImEFR/0b1792pu3rnqkcRPxx3G1oIi7py6lDbNs7jky72q9dp4ZzOxqjqbca4uiNop4M2SphDcTyFgopmVjT96fjWLiXf+XrES+E/A9WZWkuh038wmA5MhGAq1mtt3rtok8ZszB7O1oJjfv7KEU4d0pkPLnCpfF3s2c/a9MwF48rJjajRW51It6hkFwAogE8gBmks63symJ/H6fCAvZrobsK7COsOBJ8Ik0Q44RdJeM3suctTORdQkM4M7zz2C5Rt3VitJONdQRL089lLgKoIv9/nASGAmQRVRdc0B+krqBawFzgHOi13BzMrP7yU9BLzoScKlU05WZnnngY/NXk2f9gcxonfbNEflXM2KennsVQRtBqvM7ETgCGBjMgWY2V7gcoKrmRYDT5nZQkkTJU2MGJdztaKwuIQH3lrBpQ/P5cO129IdjnM1KmqiKDSzQgjGpjCzj4B+yRZiZlPM7FAz62Nmt4TzJpnZpDjrfsfMno4Yr3MplZOVySPjR9Aypwnf+es7rPx8V7pDcq7GRE0U+ZJygeeA18Krnyq2LzjXoHXJbcYjl46g1OCC+2fz2TbvRNA1TFH7ejrTzLaa2U3Az4H7Ce6BcK5R6dO+BQ9fcjTbdhcz/eOkal+dqzeiNmYPB/4X6BGWIeAWYEjqQnOufhjcrTX/vWYU7Vs2TXcoztWIqJfHPgpcC3wAlKYuHOfqp7Ik8c7KzTw0YyW3nz2Upk0y0xyVc6kRNVFsNLMXUhqJcw3Aqk27mPLBZ4j3ufPcI8jM8H6hXP0XNVH8QtL9wFSgvG8nM3smJVE5V099a3geWwuKuWXKYlo1y+I3Zw7yTgRdvRc1UVwCHAZk8UXVkwGeKFyj973je7OloIh7pi2nTfMsrht7WLpDcu6ARE0Uh5vZ4JRG4lwDcu3J/di6u5jVmwsoKTWvgnL1WtREMUvSgIpjRzjnApK4+fRBAGRmiD3FJWwpKKJgTwlTF69nVL8OnjxcvRE1URwLXCxpJUEbhQAzM7881rlQWSJYu2U3J/3xDQqLSzDgisffY2heLo+MH+HJwtULURPF2JRG4VwDNn/NFvbsLSnvQ7+gqIT5a7YybckGRvfvmNbYnKuOqONRrEp1IM41VMs37sIqjJKyu6iEReu2e6Jw9ULUvp6cc9U0sEsrmmXve/Nds+xMBnRplaaInEuOJwrnatiofh0YmpdLWXNEhmBIt9aM6tchvYE5V02eKJyrYZkZ4pHxIzikQwu65uYw6YJhPHrpSG/IdvVGpEShwAWSbgynu0s6OrWhOddwZGaINs2z6damOV8d2Iltu4t5as6adIflXLVEPaO4BzgGODec3gHcnZKInGsEHp7xCdf9cwFvLf083aE4V6WoiWKEmf0QKAQwsy1AdrKFSBoraYmkZZJuiLP8dEkLJM2XNFfSsRHjda5O+f6oPvRufxDX/3MBO/fsTXc4ziUUNVEUS8ok6N8JSe1Jsrvx8PV3A+OAAcC5kgZUWG0qQXchQ4HvEgyQ5Fy9l5OVye/POpx123Zz65TF6Q7HuYSiJoo7gWeBDpJuAd4CfpNkGUcDy8xshZkVAU9QYZQ8M9tpVn4F+kFAhavRnau/hvVow/gv9+LR2auZscyroFzdlfQNdwr6TJ4OzANGE3TfcYaZJfuzqCsQ25qXD4yIs70zgVuBDsCplcQ0AZgA0L179yTDcC59fvLVfhSVlHJIhxbpDsW5SiWdKMzMJD1nZsOAjw5g2/GuDdzvjMHMngWelXQ8cDNwUpx1JgOTAYYPH+5nHa7eaJadya/CzgPNzMeucHVS1KqnWZKOOsBt5wN5MdPdgHWVrWxm04E+ktod4Hadq3PWby/knMmzmLViU7pDcW4/URPFicBMScvDq5I+kLQgyTLmAH0l9ZKUDZwD7DO8qqRDwqouJB1JcGWVf5Jcg9MypwmfbivkuqcXUFDkV0G5uiVq77HjDnTDZrZX0uXAK0Am8KCZLZQ0MVw+CfgmcJGkYmA3cHZM47ZzDUbz7Cb87qwhnDN5Fr9/ZQm/OG1gukNyrlzk3mMltQH6Ajkxi5LqVdbMpgBTKsybFPP8NuC2KDE6V9+M7N2Wi47pwUMzPmHcoM4c3evgdIfkHBC9C49LCa58egX4Zfj3ptSF5VzjdP3Yw+ia24w/v7403aE4Vy5q1dNVwFHALDM7UdJhBAnDOXcADmrahL9+5yi65DZLdyjOlYvamF1oZoUAkpqa2UdAv9SF5Vzj1bdjSw5q2oTC4hJWbdqV7nCci3xGkS8pF3gOeE3SFhJc2upcY3X7ax9zx9R9q5F63vDv8udXje7Lj8ccGve1Ex6ZR/6WAqZceRw5WZlx13GuNuhALyKSdALQGnjJzIpTEtUBGD58uM2dOzfdYTh3wN5cupELH3iHy07ozU/H9U93OK6BkzTPzIbHWxbpjKJsHIoKhgK/ilKec25/x/Vtz7lH53Hf9BWMHdiJI7q3SXdIrpGK2kaxK+ZRQnBfRc8UxeScC/30lP50bJXDdU8voLC4JN3huEYq6n0U/xc7LekPVLir2jl34FrlZHHrNwbzp/8sZWtBMZ1ae1uFq31RG7Mrag70TlFZzrkYo/p14Pi+7cnwMbZdmkS94e6DsI+nBZIWAkuAO1IbmnOuTEaG2LyriFtfWsyevV4F5WpX1DOKr8U83wusNzPvycy5GvT+mq3c+8YKmmZmcPVX/bYlV3sinVGY2aqYx1pPEs7VvBMP68A3juzK3dOW8+HabekOxzUiUS+PvTrRcjP7Y7RwnHOJ3Pi1Aby59HOu+cf7vHD5sWQ3iXrhonPVF7XqaThBX09lVzqdRtBJ4JpKX+GcO2C5zbP5zZmD+d7f5nLfmyv44YmHxF0v3h3hsRLdEe5cRZHuzJb0KvBNM9sRTrcE/mFmY1McX9L8zmzXGPz17ZWcOqQzHVrmVLnu2ffOBODJy46p6bBcPZbyO7OB7kBRzHQRfsOdc7Xmki/3AqC01Cg1o0mmV0G5mhM1UTwCvCPp2XD6TOBvqQnJOVcdu4tKuPjBdziubzuuGN033eG4BizqVU+3AJcAW4DNwMVm9ptky5E0VtISScsk3RBn+fkx92vMkHR4lHida4iaZWfSoVVT7nx9KUs+25HucFwDFvWGu28BS83sDoKeY2+UdESSZWQCdxP0EzUAOFfSgAqrrQROMLMhwM3A5CjxOtdQ/fLrA2mVk8U1/3ifvSWl6Q7HNVBRKzZ/bmY7JB0LjAEeBiZV8ZqKjgaWmdkKMysCngBOj13BzGaY2ZZwchbQLWK8zjVIbVs05eYzBvHB2m1MfnNFusNxDVTURFHWh8CpwCQzex7ITrKMrux7OW1+OK8y44GX4i2QNEHSXElzN27cmGQYztVvpwzuzCmDO/HMu2sp9rMKVwOiNmavlXQvcBJwm6SmJJ904vVwFvdaXUknEiSKY+MtN7PJhNVSw4cPP7CRmJyrh35z5mCym2SQ5Vc/uRoQ9aj6NvAKMNbMtgIHA9cmWUY+kBcz3Y04w6lKGgLcD5xuZpsiRetcA5fbPJvm2cE42zOX+8fEpVbUq54KzOwZM1saTn9qZq8mWcwcoK+kXpKygXOoMKaFpO7AM8CFZvZxlFida0x++9JHXPzXd1i2YWe6Q3ENSNrOU8OOBC8nODNZDDxlZgslTZQ0MVztRqAtcI+k+ZL8lmvnEvjBiX1onp3JdU+/T0mp18K61EhrhaaZTTGzQ82sT3hvBmY2ycwmhc8vNbM2ZjY0fMS9vdw5F+jQModfnDaAd1dv5a9vr0x3OK6BiNp7bFPgmwTddpSXYWa/Sk1YzrmozhjalX8v+JTfv7KE0f07pjsc1wBEPaN4nuCeh73ArpiHcy7NJHHLmYMZ0bstUTr9dK6iqJfHdqsLPcU65+Lr2CqHv333aEpKjS0FRRTsKWHq4vWM6teBTB972yUp6hnFDEmDUxqJcy6lSkqNcyfP5OP1O8nfupsrHn+PCx+Y7Y3cLmlRE8WxwLywQ78Fkj6QtCCVgTnnDsy0JRv4YO328umCohLmr9nKtCUb0hiVq4+iVj2NS2kUzrmUW7huO4XFJfvM211UwqJ1272R2yUlUqIws1WpDsQ5l1oDu7SiWXYmBUVfJAsDmmZ5Nx8uOVHPKJDUBugLlI/FaGbTUxGUc+7AjerXgaF5ucxasYlSg5ysoC+osYM6pzs0V89EHY/iUmA6wV3Vvwz/3pS6sJxzByozQzwyfgSHdGhBt9xm3H3ekbz38zF0P7g5paXGQ2+v3K9qyrl4op6DXgUcBawysxOBIwDv39u5OiYzQ7Rpnk3XNs0Y3b9j+djas1du5qZ/LeLb985k3dbdaY7S1XVRE0WhmRVCcJe2mX0E9EtdWM65mnRMn7bce+EwVmzcxdfveovZK7zHWVe5qIkiX1Iu8BzwmqTnidNFuHOu7jp5YCee++GXaJWTxfn3z+apuWuqfpFrlKJe9XRm+PQmSf8lGDf75ZRF5ZyrFYd0aMlzl3+Z6/6xgF7tDkp3OK6OitqYLUkXSLrRzN4A5gNDUxmYc652tMrJYtKFwziq58EAPDZ7tbdbuH1ErXq6BzgGODec3gHcnZKInHNp8/nOPdw6ZbG3W7h9RE0UI8zsh0AhgJltAbJTFpVzLi3atWjKszHtFg/P+MR7oHWRE0WxpEyCGz2R1B4oTVlUzrm0KWu3OOHQ9vzihYX873Mfpjskl2ZRE8WdwLNAR0m3AG8DtyZbiKSxYceCyyTdEGf5YZJmStoj6ZqIsTrnktQqJ4v7LhrOlaP7Mrhr63SH49Is6lVPj0qaB4wOZ309vJei2sIzkruBMUA+MEfSC2a2KGa1zcCVwBlR4nTORZeRIa4ec2j59MsffsbBB2VzdK+D0xiVS4eoQ6EOB/6XL4ZCvUwSZjYkiWKOBpaZ2YqwzCcIRs0rTxRmtgHYIOnUKHE651KjtNT48+tLWfLZDm48bQAXjuyBFH8ApNtf+5g7pi6ttKyrRvflxzEJyNV9UTsFfBS4FviA6G0TXYHYO3zygRFRCpI0AZgA0L1794jhOOcqk5EhHp8wkh8/MZ8bn1/IB/nbuPmMQeRkZe637o/HHFqeCM6+dyYAT152TK3G61IrahvFRjN7wcxWmtmqskeSZcT7ORLp8gozm2xmw81sePv27aMU4ZyrQnm7xVcO4R/z8jn73pnsLvJOBRuDqGcUv5B0PzAV2FM208yeSaKMfCAvZrob3g2Ic3VaRoa4+qv9GNi1Ne+t3kqz7P3PKFzDEzVRXAIcBmTxRdWTAckkijlAX0m9gLXAOcB5EeNxzsWI107Q84Z/lz8/0HaCkwd24uSBnQD4cO023lu9hQsStFu4+i1qojjczAYfyIbNbK+kywnGssgEHjSzhZImhssnSeoEzAVaAaWSfgQMMLPtlZXrnNu3naCmPf7Oah6dvZoFCdotXP0WNVHMkjSgwqWsSTOzKcCUCvMmxTz/jKBKyjlXR918+iDaHpTNna8v4+MNO5l0wZF0bt0s3WG5FIramH0sMD+8WW6BpA8kLUhlYM65+qGs3eLeC4exbP0OTvvzW6z8fFe6w3IpFPWMYmxKo3DO1XsnD+zE85d/mfvfXElem2aUlBpbCooo2FPC1MXrGdWvA5kZ3oZRH6mhdfg1fPhwmzt3brrDcK5RKyk1zpk8k7mfbMGA5tmZDM3L5ZHxIzxZ1FGS5pnZ8HjLolY9OedcpaYt2cD7a7aV3xhVUFTCOys3c9frS9lb4v2H1jeeKJxzKbdw3XaKKySEvaXGn19fRmmYPaYt2cCMZZ9TWOw37dV1UdsonHOuUgO7tKJZdiYFMXduN8/O5CdjDiW7SfD79I+vfcyC/G1kN8ngyO65HNO7Hccf2o4jureJtE3vY6rmeBuFcy7lSkqNCx+YzawVmyi1+G0U2wuLmbNyMzOXb2Lmik0s+nQ7Ywd24i8XDAPgbzM/YVDX1gzp2pommclVfngfU8lL1EbhZxTOuZTLzBCPjB/BuDumU7CnhF+ePnC/q55a5WQxun9HRvfvCMDWgiJ2FO4FYMOOQm58fiEAB2VnclSvgzmmd1vGDupEj7YH1f4ONXKeKJxzNSIzQ7Rpnk2b5pQng0Rym2eT2zwYUblDyxzm/uwkZq/YzIzlnzNzxSamLdlIuxZN6dH2INZsLuDVRev5Up+29OvYkoyYBOSX5aaeJwrnXJ3UrkVTTh3SmVOHdAZg/fZCmoedEM5euZmbXww6hmjTPIuRvdtyTJ+2nHZ4F3746Lss27CTUoMrHn/PL8tNAU8Uzrl6oWOrnPLnZw3rxjF92gbtG8s3MWvFJl5e+Bm5zbKYv2Zr+ZVVBUUlvLt6C8/NX8uZQ7vuc+bhqs8ThXOuXuqa24yzhnXjrGHdMDM+3VbI0/Py9xsjo7C4lJ889T63vfQR7/zvSQDc/+YK1m8vpGOrHDq3bkan1jl0zQ3+RpGuK65qa7ueKJxz9Z4kuuQ2i3tZbtMmGXxrWDf6dGhRPu/d1VuYungDe/Z+ca/HwC6t+PeVxwFw3dPvs333Xjq1zqFz6xw6tc6hT/sWDOraOu72y3rrLSm1hA34qVZbowl6onDONRij+nVgaF7ufpfl/vL0Qft8Yd9z/jDMjG27i/l0WyGfbSskdiiNPXtLWbZxJ28u3ciuMOmc1L8j918cXD067o43ycyATq2alSeSod1yuXvasi/aRx57jyF5rfn7+BFJX96brJpuwPf7KJxzKZXuG99S/at+R2FxeSI5pENLzIyfPfcha7fu5rNthXy2vZCtBcWc2K89s1du3udspkx2kwwmntCHq8ccyvbCYr5xzwyaNskgJyuTpk0yaNokg7OG5XHqkM5sLSji9tc+/mJZ+PdLfdoxoEsrthcWM3P5pvLXN8kQv/73Ihbkb6v0npXq8PsonHO1pjYHTYon2ctyq9IyJ4uWOVnl05K45cx9x23bXVTCn19fyrQlG/d7/TG923J4Xi5HdM8FwAwO7diCPcWl7NlbSmFxCTsK97KrKLiHZEfhXp5/fx2FxSXs2VtK2W/5m88YxIAurVizuYDLHplXabwFRSXMX7OVaUs2pGT/wROFc84dsGbZmQzr0SZutyWXHtdrny/s1s2yuOf8YZWWlXdwc+bf+FUAzIyikiChZIfVV33at+DFK45lz95S9hSX8I95+Tz73tp9ythdVMKiddtTlijS2imgpLHh4EfLJN0QZ7kk3RkuXyDpyHTE6ZxzVSlrHymr7SmrAhrVr0PkMiXRtEkmrXKyyoeYzcnKZFDX1gzr0YYvHdKOrw3pXH5/SZlm2ZkM6NIq8nYrStsZhaRM4G5gDJAPzJH0QoXhVccBfcPHCOAv4V/nnNtHvLaRnjf8u/x5TbaNxNt2QVEJM5Zvos//TKnRbVfWgH8gCaqitDVmSzoGuMnMTg6nfwpgZrfGrHMvMM3MHg+nlwCjzOzTysr1xmznXGOTigb8utqY3RVYEzOdz/5nC/HW6QrskygkTQAmAHTv3j3lgTrnXF0U70xm/MNf/FBuCDfcxUt3FU9vqrMOZjYZmAzBGcWBh+acc3VfbV1hls7G7HwgL2a6G7AuwjrOOedqUDoTxRygr6RekrKBc4AXKqzzAnBRePXTSGBbovYJ55xzqZe2qicz2yvpcuAVIBN40MwWSpoYLp8ETAFOAZYBBcAl6YrXOecaq7TecGdmUwiSQey8STHPDfhhbcflnHPuC2m94c4551zd54nCOedcQp4onHPOJeSJwjnnXEINbjwKSRuBVemOI4J2wOfpDqKW+T43Dr7P9UMPM2sfb0GDSxT1laS5lfWz0lD5PjcOvs/1n1c9OeecS8gThXPOuYQ8UdQdk9MdQBr4PjcOvs/1nLdROOecS8jPKJxzziXkicI551xCnihqmaSxkpZIWibphkrWGSVpvqSFkt6o7RhTrap9ltRa0r8kvR/uc73uJVjSg5I2SPqwkuWSdGf4fiyQdGRtx5hq1djn88N9XSBphqTDazvGVKtqn2PWO0pSiaSzaiu2VPNEUYskZQJ3A+OAAcC5kgZUWCcXuAf4upkNBL5V23GmUnX2maCH4EVmdjgwCvi/cIyS+uohYGyC5eOAvuFjAvCXWoippj1E4n1eCZxgZkOAm2kYjb0PkXify47/2wiGU6i3PFHUrqOBZWa2wsyKgCeA0yuscx7wjJmtBjCzDbUcY6pVZ58NaClJQAtgM7C3dsNMHTObTrAPlTkd+JsFZgG5kjrXTnQ1o6p9NrMZZrYlnJxFMFplvVaN/zPAFcA/gXr9OfZEUbu6AmtipvPDebEOBdpImiZpnqSLai26mlGdfb4L6E8wzO0HwFVmVlo74aVFdd6Thmw88FK6g6hpkroCZwKTqlq3rkvrwEWNkOLMq3h9chNgGDAaaAbMlDTLzD6u6eBqSHX2+WRgPvAVoA/wmqQ3zWx7DceWLtV5TxokSScSJIpj0x1LLfgTcL2ZlQQny/WXJ4ralQ/kxUx3I/gVXXGdz81sF7BL0nTgcKC+Jorq7PMlwG/DEQ2XSVoJHAa8Uzsh1rrqvCcNjqQhwP3AODPblO54asFw4IkwSbQDTpG018yeS2tUEXjVU+2aA/SV1CtsrD0HeKHCOs8Dx0lqIqk5MAJYXMtxplJ19nk1wRkUkjoC/YAVtRpl7XoBuCi8+mkksM3MPk13UDVJUnfgGeDCenx2nBQz62VmPc2sJ/A08IP6mCTAzyhqlZntlXQ5wRUQmcCDZrZQ0sRw+SQzWyzpZWABUArcb2YJL7+ry6qzzwRXwTwk6QOCapnrzay+ddFcTtLjBFdvtZOUD/wCyILy/Z0CnAIsAwoIzqjqtWrs841AW+Ce8Bf23vreu2o19rnB8C48nHPOJeRVT8455xLyROGccy4hTxTOOecS8kThnHMuIU8UzjnnEvJE4ZxzLiFPFM455xLyROEaLEm5kn6Q7jiqUtNxSvpeOL7JfEmlMc//WFPbdA2L33DnGixJPYEXzWxQHYhFBJ+3/XrFjRJnovISvKYrMMPMelT3Nc6Bn1G4hu23QJ/w1/PvJV0g6Z1w+t5wUBkk9ZT0kaT7JX0o6VFJJ0l6W9JSSUfHrPNwOErb02FfXIRl7Fd2+JrFku4B3gXyJD0Xdh+/UNKESuLsGTtqmqRrJN1USXlx96kSgwi6cXcuKZ4oXEN2A7DczIYCDwJnA18Op0uA82PWPQS4AxhC0HPteQRdYV8D/E+4Tj9gcjhK23bgBwCS+icoux/BIEVHmNkq4LtmNoygZ9ErJbWNjdPMrq1in8rLA5pXsU8VDQbqbb9hLn28U0DXWIwmGOdjTtgpXTP2HXVspZl9ACBpITDVzCzsqLBnuM4aM3s7fP534ErgDwnKng6sCkexK3OlpDPD53kEw6F+lsR+xJZX1T5VNAh4LYltOQd4onCNh4CHzeynlSzfE/O8NGa6lC8+JxUb9Mqm45Ydtj3sipkeBZwEHGNmBZKmATlxYtnLvmf7sevsinle1T5VNBi4vZrrOlfOq55cQ7YDaBk+nwqcJakDgKSDJSXbqNtd0jHh83OBt5IsuzWwJUwShwEj48QJsB7oIKmtpKbA1yqJp9r7JCmD4Ozlo+rsqHOxPFG4BiscRe3tsGH4EuBnwKuSFhBUwXROssjFwMXh6w8G/hJuZ1E1y34ZaBKuczMwq2Kckn5vZsXAr4DZwItU8uWexHYhaIPJN7M9lSx3rlJ+eaxz1VCXLrV1rrb5GYVzzrmE/IzCOedcQn5G4ZxzLiFPFM455xLyROGccy4hTxTOOecS8kThnHMuIU8UzjnnEvJE4ZxzLqH/B2SlrH3TNIpzAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "# Plotting\n", - "# YOUR CODE HERE\n", - "raise NotImplementedError()" + "\n", + "# As copied from last week:\n", + "def jackknife_batch_estimate(data,observable,k):\n", + " '''Divide data into k batches and apply the function observable to each \n", + " collection of all but one batches. Returns the mean and corrected \n", + " standard error.'''\n", + " batches = np.reshape(data,(k,-1))\n", + " values = [observable(np.delete(batches,i,0).flatten()) for i in range(k)]\n", + " return np.mean(values), np.sqrt(k-1)*np.std(values)\n", + "\n", + "with h5py.File('xy_data.hdf5','r') as f:\n", + " k = 100\n", + " mean_sq, std_sq = np.array([jackknife_batch_estimate(data, np.mean, k) for data in f[\"square-magn\"]]).T\n", + " plt.figure()\n", + " # I do plot a line to visualize the trend.\n", + " plt.errorbar(temperatures, mean_sq, yerr=10*std_sq,\n", + " fmt='.--', markersize=10, capsize=4)\n", + " plt.xlabel(\"temperature $T$\")\n", + " plt.ylabel(\"mean square magnetization $\\\\overline{m^2}$\")\n", + " plt.title(\"Mean square magnetization for the $w \\\\times w$ XY model \\n\"\n", + " \"(w = {}, #measurements = {})\"\n", + " .format(width, measurements))\n", + " plt.show()" ] }, {