From a28ef38f00487ad841c6fb4d158a18007a9ba79a Mon Sep 17 00:00:00 2001 From: Fabian Joswig Date: Tue, 3 May 2022 10:44:03 +0100 Subject: [PATCH] examples: Example 06 updated to work with new default options for GEVP --- examples/06_gevp.ipynb | 46 ++++++++++++++++++++---------------------- 1 file changed, 22 insertions(+), 24 deletions(-) diff --git a/examples/06_gevp.ipynb b/examples/06_gevp.ipynb index 73b6f980..9bb9ed49 100644 --- a/examples/06_gevp.ipynb +++ b/examples/06_gevp.ipynb @@ -7,7 +7,6 @@ "metadata": {}, "outputs": [], "source": [ - "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pyerrors as pe" ] @@ -19,8 +18,7 @@ "metadata": {}, "outputs": [], "source": [ - "plt.style.use('./base_style.mplstyle')\n", - "plt.rc('text', usetex=True)" + "plt.style.use('./base_style.mplstyle'); plt.rc('text', usetex=True)" ] }, { @@ -111,7 +109,7 @@ ], "source": [ "try:\n", - " matrix_V1V1.m_eff() #This does not work! \n", + " matrix_V1V1.m_eff() # This does not work!\n", "except:\n", " print(\"Something is wrong\")" ] @@ -165,12 +163,12 @@ "name": "stdout", "output_type": "stream", "text": [ - "[ -0.71920537 162.3903938 -492.21321233 714.06804297]\n" + "[ 1.97277730e-01 -8.99401483e+01 2.71416300e+02 -5.20598763e+02]\n" ] } ], "source": [ - "vec = matrix_V1V1.GEVP(t0=3,ts=6,state=0)\n", + "vec = matrix_V1V1.GEVP(t0=3, ts=6 ,state=0, sorted_list=None)\n", "assert len(vec) == matrix_V1V1.N\n", "print(vec)" ] @@ -192,9 +190,9 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk8AAAGLCAYAAADeVnZRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA3z0lEQVR4nO3df3xU9Z3v8fd3gAQoJMMPw28p4UcRddUoWF2xd9sE8Lr1ul1A4bp2u7ZQ+2i1Ky6U2ta11SIstMV2lwWvdrUWqHDvttqtRLLbLrioIKFdeKASCFV+xwpJoEKQzPf+cc6EyWTOzDmTTGYyeT0fjzySOT/mfOcrkjff7/d8jrHWCgAAAP6Est0AAACAroTwBAAAEADhCQAAIADCEwAAQACEJwAAgAAITwAAAAEQngAAAALome0GdDRjjJE0XNLpbLcFAAB0Kf0lHbUpimDmXXiSE5wOZ7sRAACgSxop6UiyA/IxPJ2WpEOHDqmoqCjbbQEAAF1AY2OjRo0aJfmYucrH8CRJKioqIjwBAIAOx4JxAACAAAhPAAAAARCeAAAAAiA8AQAABNApC8aNMQvdH8dKkrV2vs9z6t2XYWvtssy0DgAAwL+MhydjzFJr7aKY16uNMZuttRVJzlkoSdbaNe7rcmPMaj+hCwAAIJMyOm1njAlLKnO/R62WVG6MKU1y6mJJa6IvrLVVkuZloo0AAABBdMaap+skxQalWvd7ONHBbqgKW2vrE+wr7+jGAQAABJHRaTs3AA2I2xwNQLVKzGtEql4JApcxplBSYcym/r4bCAAAEFA27rZbLGl+opGlFE5KGujxfg0xXzzXDgAAZEynhidjzFJJP4suBA8oUXCSpCWSimO+RqbZPAAAgJQ67dl2xpiZkg74CE5e03nhRPustU2SmmKuk24TAQAAUuqUkafoQu+Y0gNhr7vtrLW1kuoT7XfvugMAAMiajIcnY0yZpDJJ1caYUjcUzZOzhknutoVxpy3RxYXl0VGrdKb6Uos0Swe3Srs3Ot8jzRm5DAAAyA/GWpu5N3fqOx1UgrvkrLXGPWaepEXW2rFx5y7UxWm6ybGFNlNcs0hSQ0NDg4qKipIfvPcF6eWHpPp3L24LXypNe0yadJufywEAgDzQ2Nio4uJiSSq21jYmOzaj4SkbfIenvS9Iz98tTZghTV0glVwm1b0pbV0h7dskzX6WAAUAQDdBeEoVnhqOSP/nU9Kwq6Q710mhmNnLSERaP0c69t/S56uk4hGd0m4AAJA9QcJTNuo8ZV/l16XTx6SpD7YOTpLzeuoC6fRR5zgAAIAY3TM8nf/A+V5yWeL90e3R4wAAAFzdMzyV3eV8r3sz8f7o9uhxAAAArrwPT80Rq1cPvK9f/PaIXj3wvpojVpr4585ddVtXOGucYkUi0tbvSeHRznEAAAAxOq3CeDZs2nNMj7y4V8cazrVsG1bcWw9/epJmTHvMudtu/Vxp6gMxd9t97+LddqEeWWw9AADIRXl7t93GV9/Wgz+vUfyniz68ZdVdZZoR2pGgztNoadqjlCkAAKAboVSB1HD1Q/+qUxd6eR43uF+BXvzyTRpWVCC9s006c0LqN0QafSMjTgAAdDOUKpD0/h8/TLr/D2fOa/nLbztBacxU6cqZzneCEwAASCJvw5MfYy/pl+0mAACALqZbh6cD753JdhMAAEAXk7fhadBHvNc7Sc6apwenfayTWgMAAPJF3oanb/z5JBldvLsuKrrt0duv0LBwn85vGAAA6NLyNjxVTBqqVXeVaWhx71bbhxb3dsoUXDEsSy0DAABdWd6WKmhoaFBRUZGaI1bbD55U3elzKunfW1PGDFSPUPx4FAAA6M6ClCrI6wrjktQjZHTD2EHZbgYAAMgTeTttBwAAkAmEJwAAgAAITwAAAAEQngAAAAIgPAEAAARAeAIAAAiA8AQAABBA3td5ardIs/TONunMCanfEGn0jVKoR7ZbBQAAsoTwlMzeF6SXH5Lq3724LXypNO0xadJt2WsXAADIGqbtvOx9QXr+bqnkcumeKmnxEed7yeXO9r0vZLuFAAAgC/L+2XZpiTRLT1ztBKU710qhmIwZiUjr50p1e6X7djGFBwBAHgjybDtGnhJ565fOVN3UBa2Dk+S8nvqAVP+OcxwAAOhWCE+J7HrO+V5yWeL90e3R4wAAQLdBeEqkZ1/ne92bifdHt0ePAwAA3QbhKZEZ35X6D5O2LnfWOMWKRKStK6T+w53jAABAt0J4SqR4hHTLMmlfpbM4/NB2qem08339XGf7LUud4wAAQLfC3XbJJKzzNFqa9ih1ngAAyCNB7rajSGYyk26TJt5KhXEAANCC8JRKqIc0Zmq2WwEAAHJExtc8GWPCxph5xpjNPo8vN8ZscM8pN8YsNcbMzHQ7AQAA/MjoyJMxpkzSdZLCkgb6PC0sqVzSTEm1kpZaazdmon0AAABBZTQ8WWurJVWnMXI0xlpbn4EmAQAAtAulCgAAAALI1QXjs40xJ+VM9Y211i7yOtAYUyipMGZT/0w3DgAAdF+5GJ6qJclaWytJ7sLxDdbaWR7HL5b0cGc1DgAAdG85N21nra2NBifX85JmGmPCHqcskVQc8zUysy0EAADdWc6Fp/jF5TELx0sTHW+tbbLWNka/JJ3OcBMBAEA3llPhyR1d2mCMKY3bJjllCwAAALKqs8JTwhpPxphSY8zC6Gt3lGlZ3LTdPEkbc7Z0QaRZOrhV2r3R+R5pznaLAABABmW6SGapnGKXd0gqM8YslbQjpuhluaT5kpbFnLYkNlBJGpRksXh2JXxw8KXStMd4cDAAAHnKWGuz3YYOZYwpktTQ0NCgoqKizF1o7wvS83dLE2ZIUxdIJZdJdW9KW1dI+zZJs58lQAEA0EU0NjaquLhYkordNdSeCE/piDRLT1wtlVwu3blWCsXMfkYi0vq5Ut1e6b5dzoOFAQBATgsSnnJqwXiX8dYvnam6qQtaByfJeT31Aan+Hec4AACQVwhP6dj1nPO95LLE+6Pbo8cBAIC8QXhKR8++zve6NxPvj26PHgcAAPIG4SkdM74r9R8mbV3urHGKFYk4i8b7D3eOAwAAeYXwlI7iEdIty6R9lc7i8EPbpabTzvf1c53ttyx1jgMAAHmFu+3aI2Gdp9HStEcpUwAAQBcS5G67jBbJzHuTbpMm3iq9s006c0LqN0QafSPlCQAAyGOEp/YK9ZDGTM12KwAAQCdhzRMAAEAAjDx1hkgzU3sAAOQJwlOm8fBgAADyCtN2mRR9eHDJ5dI9VdLiI873ksud7XtfyHYLAQBAQJQqyBQeHgwAQJfBg4FzAQ8PBgAgLxGeMoWHBwMAkJcIT5nCw4MBAMhLhKdM4eHBAADkJcJTpgR5eHCkWTq4Vdq90fkeac526wEAgAfutsu0VA8Ppg4UAABZF+RuO8JTZ/CqMB6tAzVhhnNXXsllzlqorSukfZuk2c8SoAAA6ASEp1wLT4lQBwoAgJxBnaeuYNdz/upAUcoAAICcQnjKlld/5HxPVQfq1R+xoBwAgBzCg4Gz5YYvSy/e56xxGjW57f5oHaiPTnWm91hQDgBATmDkKVuuucsJQVtXeNSB+p7U9xLpjad5sDAAADmEBePZ1Opuuwdi7rb7nrTvJanPQGnUFOnOdQkWlM+Rjv239Pkqqf/QxHfzAQAAX7jbrquEJ8m7DlS/IdLh7c5IU6JpvUPbpacqpJFTpDPHU0/reZVLAAAAgcITa56ybdJt0sRb2wabtXc6+1MtKD+8wxm5+sunW9eJev7ui3WiKMQJAECHITzlglAPaczU1tvK7pL2v+y9oPz4Hue8cRWt60SNmuy8Xj9HemmR9MFJ6ZdfTR2wGJkCAMAXpu1yVaoimk9Pc0adUk3rFfRz7thLVoiz4tvS5m8mH5kiXAEA8hjTdgE0R6y2HzyputPnVNK/t6aMGageIeN7f8aEejjh5fm7ZdfPkYl5fIvdukLm8A7nuFTTeufPJC/E+VSFtOGvk49MSamn/fyEKwKYL1n7MwcA8KVbh6dNe47pkRf36ljDuZZtw4p76+FPT9KMK4al3C/5+0WXdkCbdJt23bBSw19/TEP2VbQcfyI0VKfHf0Hja55MXSdK8g5YgyekmPqbK/3yAemDP7Q/XPlZd5UqXHVG+GrnNZovXNBbr1fq7Kkj6jNghCZeP109esb9b5bkGn7+zHWETAc0AiCAfNZtp+027Tmme5+rVvynj/71Pu/mMVqz5aDn/lV3lUlSyl907QloknTvc9UyimhK6C2VqF51CmtHZKIkac+gr6nPiCu9p+QObpE+/KP31N72J6VfPei9/51XpWdudcNVonIJc6XDb0gf/EF2wnSZqQ/GjI4tl9lXeTFcPX93ymPsyw/JxIQrG75UJhqufIQvP8El6TF7X0jehhTn76p8xgm6kRMt558IDdHR6x/SNdM/m/QajVP/Xj87c7W++9JbCsX9994emaiIQvr6LRN1+zUjVNKvV7tC5qY9x/SdF3Zr1JnftVzjUL+r9M3brrwY0NoRIn29PwDkGEoVpAhPx+rP6tM/ekV/OHPe+32kNsEpVv/ePXX63IWE50kXw1W6Ac0muUbUrL7VWhZZITt+ukI3X5zWi2xZIVNTKfPpH0i/eVx22FUyCcKPXfMJmeP/7RTfLOzX9gL7KqW1s1OGKzuuIvH7r58jc2KvZCRbMsn7mMM7Zd0AFooJV5Ety53PceNXpG0/dOthLWg9+rVvkzT7We06dCplcEkabkYNkHUDXsI2JLlGzdVf0wfnL6h8z0Jp/HSFbm59vmoqVXXFMl0/ZpCKXrzH8xpfPH+/JOnhwrUabutarnHUlOiRprmqjEzRyqsO6bYT/5hWyKwbOU0/33VEOyuf9bzGtdPv1h39fquirX8fOES+/8EFX+9/+zUjVFLUm2lcADklp8KTMSYsabakWdbaihSHR89ZKKnefRm21i4LcL2U4WnB87/V/60+4vct0zKgby9J0qkPPvQ8JlVA82N6aHuCX1RD9O2mOfqz2+/RgHc3xfxSbx2wVLPJKTHvFY5+87j0myXph6voovVk14gGsLEVMnPahqvI2julg79WaOwnPUe/mt7Zrl5NJz2Dy+9uWClJuurV+z2O2aRzBYPUe/RkhRK1Yd0cnf39dvX58JTn+fUmrPC4j3uef7ymWn0KQioefbXnMfU1ryqsRs/PsebDWzW/17+lDJleI3y/+Njj+tXuY1pVsDL9ayQJkb8ouVc7361P+v73nr9fE/7HXC0Y+XbqkUSm/gB0opwJT8aYMknXSQpLusNae62PcxZKUjQwGWPK5QSv+T6vmTI8/dOv92tZ5du+PkNX4DXN06dXSGc/jCQJWHfoW4XrNXRcWeJf6E9NU+hIkjv6UoWrptPSkpHOz+kGMJ9Ti57ha90c1dVUyxjpknQ/Z6prpDq/g0Kkrf21zNhPJv4Ma++Uan8tM+6TniN8597dpfomqyFe/ZDqGuvm6IODr6vvhfr0Q+T+XTpx9Vd0za5vJh1J3BSZzNQfgE6VM3fbWWurJVUbY2YGOG2xpDEx71FljNksyVd48uPAe2c66q1yQkQhvRaZ1Gb72Q+dZ+ZVRqZo89nrEgYs2xTSqpqViqyb03Zk6sgOnTRhhbcsT/zLsKbKGbnys2jd65iju5zvXovae/ZOvr+pUYo0y9z8YMI7CkM3L9DQmgpniM/rmPHl0pEd6V8j1fmx29O9xoQKaX+l92eI7p+aeL+ZukB9nqpQHyn1e3jtv+mr6nfgVtnx01sHvFHOiF3kqWkaeGSH9/k3L9DwmgoN3L1UmjA94U0Kdv0c/fGFv9MvGu/QhsJ1Gl4QE/g/LNEja+fq3dipPwDIgpx6MLAxplTONF19gn3lHXWdB6d9TIP7FSRvS0ddLEdEA9YLkRv1WmSSIu5/+srIFN17/n4d31/tjI4sGSk9VaHj+3fp3vNf1eKmv5ZqKhVZN8cZQWk6LR3a7rx2w1Vky/KEDzeObFmhoyrRMVPifUxNlfNzbNCKdeFc8v2pwpef4DL8mvZdI9X58SEynWukCpGp9vvph/YG1fHlvtrQ+8LppCGv37nj+qeCJzR0XFmrB2IPHVemVQUrtbPyWf3ktXcSXwMAOkFOhSdJpR7b6+VM/bVhjCk0xhRFvyT1T3WRYeE+evT2K2TUNiRFt827eYznfslZzJ3MgL69WtY9eUkV0KLX8GpDR6mMTNFNZ7+nO89/Q/ed/7LuPP8N3XR2hSojU9oXrmoq9cj5ufr7prnpB7B9m9UcKkg/fPkJLoVFUqhH+tdIdf6WFTqW6RCZar+ffmhvUA0SIlOVzxg/3RntHDXZme51R7c0froeLlyny4ubEp8PAJ0g18KTl5OSBnrsWyypIebrsJ83nHHFMK26q0xDi1sP/Q8t7q1Vd5Vp8f+c5Ln/n+8q0z/M/JOk4WvJZ67Uks9cmXZAM5L+Yeaf6J892rDkL65IOXoWJMB5jUwZJQ9Xm613uPrS+fu16yM36Y0+N6UfwPa/rCebprVr9CtlcHnlBzoZ6Zf+NVKdX1OpF4d+yVeItFsTX8PWpAiRqUKmOwp4qleSfmhvUPURIs+EipO/x+4NKadhh9sT6rf72cTnA0An6JRSBe6ap8WpFoy7U3ObrbUmbvspSYustWsSnFMoqTBmU39Jh/0+nqU9Fcb9FDTMZCHOaK0qqfVde4nKJXgdEy2X0N79IUU0Oa4WVUShVm1IdEyzQvrMNSP0x9/9a8JF7Y80zVFlZIrnovdHmuY4n7XlDq/4Oworta3se3GlBNoe88yIb+u1g+97XuPjYwbps0e+5Xl+1RXL1Legp8b/9vFWd6EdDw3Vseu/rhE33JHkNn7nGn9ZNkIVexZKE6a3qSivfZVa/eGtmtfr3zzbsCbF/nvP368bSgfps4e/5VneIvk1Ui8Ir6/ZpgE67fn+e4b9pQYd3+J9k8KaTyiUrHyGexPC7yYt0hUzF3M3HoAOkzN327VcxH94KpV0IEF4spIqrLVVPq7Vqc+2y2iFcR86I8B1RAD08x6J7q4aNai/Xjt4UpL3XYWS9JWhe3VX45MJg0uyOk+xx6Qq7pjqfCl1oc6UBSQT1mkarYapD+tnZ65OGr629rxBUy+86rl/Z9+b9OKXb9Kwo5sTFOocrUYf17hyRJG+9N53kobI68cMSlAnarTMtEf1i8Zx+tULP0sSdlOUz3DvXNz4kblacWFmxiuxA+g+umx4co89Jelaa21tzDYbH6iSnJ8fDwYOoDMCXEcEwHTeo67xXMqCpoP7FejFL9+kkn692ldh3E8b/Tx+JYWUfZXi8S1e4UvyHuGLjgL6qSDeESHS6/1f2n1M9/60Ov3yGevmqOH3u3Tt6RWS1KbyfpvPCQA+5WJ4midpfnx4ckeaZsYWwYwWyIxO0bnBq6Ij6zyha/EzPdmdflm2dyq5vdeQ0g+RK15+Wz/8j/2SvEcSp4e2658LVnpO/S0MLVDjuQueVcxbRtjCfXx/XgDImfAUDUeS7pBUJmmZpB3W2o3u/nly1jKNjTtvoaToyNNka+2iANckPOWhznpgbj7I5crcdY3n9PNdR/Tdl97yPObrt0z0eETMaD3T7/N6tfb9lFXM+139F1ox++pO+EQA8kXOhKdsIDzlr1wOBQjGVxhOMPX33Gu/1yc3z0g6rXd8/y79R8VLuuvGsQIAv3KmwjjQkXqEjG4YOyjbzUAHmHHFMFVMGpo8DId6SGOmtjpvQsM2Z6ouRRXzCQ3bJBGeAGQG4QlAVqQThge+vc75IUUV84Fvr5Nu+av2NA8APHWVIpkAoJGXuLVyU1QxbzkOADKA8ASgy+j950ul/sO8K7FvXSH1H+4cBwAZQngC0HUUj5BuWSazr1J2fevH3Nj1c2T2VUq3LHWOA4AM4W47AF2PRyV2TXtUmnRb9toFoMvibjsA+W3SbdLEWz2rpANAJhGeAHRNCUoZAEBnYM0TAABAAIQnAACAAAhPAAAAARCeAAAAAiA8AQAABEB4AgAACIDwBAAAEADhCQAAIADCEwAAQACEJwAAgAAITwAAAAEQngAAAAIgPAEAAARAeAIAAAiA8AQAABAA4QkAACAAwhMAAEAAhCcAAIAACE8AAAABEJ4AAAACIDwBAAAEQHgCAAAIgPAEAAAQAOEJAAAgAMITAABAAIQnAACAAAhPAAAAAfTsjIsYYxZKqndfhq21y1IcXy5pvqTNkmolVUjaYa3dmMl2AgAApJLxkSc3OMlau8Zau0ZStTFmdYrTwpLKJa12vw4QnAAAQC7ojJGnxZLGRF9Ya6uMMZvljCwlM8ZaW5/JhgEAAASV0fBkjCmVM01Xn2BfubW2KpPXB9C9/OHoO6p/75AkqTlitfu9Zp08G9HAPiFdeUkP9QgZSVL4klEaPHx0NpsKoAvL9MhTqcf2ejlTc8nMNsaclDRQ0lhr7aJEBxljCiUVxmzqH7CNAPJEzUs/1A2HntSm5sn6zod3aVToPZWoXnUK6/uRS/TNXs9pRo8denXUFzT4nuXZbi6ALqpTFownEA1FXqolyVpbK0nGmHnGmA3W2lkJjl0s6eGObyKArmb8LV/Rj6tv1WvbfqMNfR7TcFvXsu+oKdEjTXN17PqH9Okyr3/XAUBq2SpVkCw4yVpbGw1OruclzTTGhBMcvkRScczXyI5qJICuZcDQS/X2b7dpVcFKDR1XJt1TJS0+It1TpaHjyrSqYKXe/u02DRh6ababCqALy3R4qvXYHk6yT8aYmbGvY9ZMtfnnorW2yVrbGP2SdDq9pgLo6rYfeE/3Nf+LNH66QnPWSaMmS4X9pFGTndfjp+u+5me0/cB72W4qgC4so+HJHT2qdxeOx+9LuFjcHV3aEHtOzIiTZ+ACgObf/5eG2zqFbn5QCsX99RYKKXTzAg23J9T8+//KTgMB5IXOmLZbIqdmk6SWUaU1Ma9Lo7WgpJZRpmVx03bzJG2kdAGAZEpMvfvDZR4HXNb6OABIQ8bDk1tNPGyMmekGp8nW2tgaT9Fq4rGWGGMWRr8kDfJYLA4ALcaWjnV+qHsz8QHu9pbjACANxlqb7TZ0KGNMkaSGhoYGFRUVZbs5ADpRXf0f1Xf1deo78k+cNU6xU3eRiCLr5uiDw7v1wfwdKgl/JHsNBZBzGhsbVVxcLEnF7hpqTzwYGEDe+OmOI1rQMFuqqVRk3Rzp0Hap6bR0aLvzuqZSCxpm6ac7jmS7qQC6sGzVeQKADvdXlxeqvqRCJ4701eA9Tyv0VEXLvuaPDNUfpizW3434U4UvKUzyLgCQHOEJQN4Y/PZaDf7PxxPu6/XH4xq2/bvOi098TRq+uBNbBiCfEJ4A5I/rPid97JbUx/Ufmvm2AMhbhCcA+aP/UIIRgIxjwTgAAEAAhCcAAIAACE8AAAABEJ4AAAACIDwBAAAEQHgCAAAIgPAEAAAQAOEJAAAgAMITAABAAIQnAACAAAhPAAAAARCeAAAAAiA8AQAABEB4AgAACIDwBAAAEADhCQAAIADCEwAAQACEJwAAgAAITwAAAAEQngAAAAIgPAEAAARAeAIAAAiA8AQAABAA4QkAACAAwhMAAEAAhCcAAIAACE8AAAABEJ4AAAACIDwBAAAEQHgCAAAIoGdnXMQYs1BSvfsybK1dlolzAAAAMi3jI09uCJK1do21do2kamPM6o4+BwAAoDMYa21mL2DMKUljrLX1MdustdZ05DkxxxVJamhoaFBRUVG72g4AALqHxsZGFRcXS1KxtbYx2bEZHXkyxpTKmXKrT7CvvKPOAQAA6CyZXvNU6rG9XlK4I84xxhRKKozZ1N9XywAAANKQrbvtTkoa2EHnLJbUEPN1uH1NAwAA8Jat8BQ0OCU7Z4mk4pivkek2CgAAIJVMT9vVemwPJ9kX6BxrbZOkpuhrY1KuKQcAAEhbRkeerLW1kurdReDx+6o66hwAAIDO0hnTdksktdwlZ4yZKWlNzOvSaF0nv+cAAABkS8bDk1sZPGyMmemGoMnW2vkxh5RLmh/wHAAAgKzIeJHMzkaRTAAAEFTOFMkEAADIN4QnAACAAAhPAAAAARCeAAAAAiA8AQAABEB4AgAACIDwBAAAEADhCQAAIADCEwAAQACEJwAAgAAITwAAAAEQngAAAAIgPAEAAARAeAIAAAiA8AQAABAA4QkAACAAwhMAAEAAhCcAAIAACE8AAAABEJ4AAAACIDwBAAAEQHgCAAAIgPAEAAAQAOEJAAAgAMITAABAAIQnAACAAHpmuwEA0KlOH3e+Uuk/1PkCgDiEJwDdyxs/lv7z8dTHfeJr0p8tznx7AHQ5hCcA3ct1n5M+dovzc6RZevtX0tbl0tQHpY/9TynUw9nHqBMAD4QnAN1KnQ2rzvZR0cGXNPT1R1Vw+pCzY+tynf/tz3T8+m+occwtKrGFKsluUwHkKMITgG7lp6+/q7d+/VOtKlgpjZ8uzf6xVHKZVPemem5ZrpFVX9S95+/XxD/73/rbignZbi6AHGSstdluQ4cyxhRJamhoaFBRUVG2mwMgx9TV/1F9V1+nviP/RKE566RQzE3HkYgi6+bog8O79cH8HSoJfyR7DQXQqRobG1VcXCxJxdbaxmTHUqoAQLcy6ORO9Tt7VKGbH2wdnCQpFFLo5gXqd/aIBp3cmZ0GAsh5hCcA3cqB2gPODyWXJT7A3d5yHADEITwB6FbqbNj94U2PA95sfRwAxCE8AehWenz0T3XUlCiyZbkUibTeGYkosmWFjpoh6vHRP81OAwHkvIyHJ2PMQmPMPPdroY/jy40xG9zjy40xS40xMzPdTgDdw5Sxl+iJHn8t1VQqsm6OdGi71HRaOrTdeV1TqSd6fFZTxl6S7aYCyFEZLVUQDUvW2jXu63JjzGpr7fwkp4UllUuaKalW0lJr7cZMthNA93Hq+Lv62NU36t5t9+vh/Ws1vKaiZd9xM0SPnL9fH7/uRp06/q4GDx+dxZYCyFWZrvO0WNKY6AtrbZUxZrOkZOFJksZYa+sz2TAA3VPNSz/U5w49qWE9JmvW2Yc0KvSeSlSvOoV1KDJY3+z1U8144wd69cQXNPie5dluLoAclLE6T8aYUkkHrLUmbruVVGGtrfI4b6akKr/hyRhTKKkwZlN/SYep8wQgkT8cfUf17zlVxZsjVrvfa9bJsxEN7BPSlZf0UI+Q81dW+JJRjDwB3UiQOk+ZHHkq9dheL2dqLpnZxpiTkgZKGmutXZTk2MWSHg7cOgDd0uDho1uFoo9lsS0AuqZs3G0XDUVequWMPG1010odMMZsSHL8EknFMV8jO6ylAAAAcXyPPLnTaXf4OHSJtbY6yf5kwUnW2tq4Tc9LWm2MCSeayrPWNklqimmnjyYCAACkx3d4cu94C3LXW3wIigon2SdjzMzYu+ustfVuICqVMyoFAACQNRmbtnNHkOrdhePx+7wWi4clbYg9x90mJQlcAAAAnSXTa56WyKnZJKll6m9NzOvS2MKZ7rTcsripu3mSNlK6AAAA5IKMlSpouYATjqJhaHLsnXPGmHmSFllrx8ZsC8sJTFGDUtxtF3+9IkkNlCoAAAB+BSlVkPHw1NkITwAAIKgg4YkHAwMAAARAeAIAAAiA8AQAABAA4QkAACAAwhMAAEAAhCcAAIAACE8AAAABEJ4AAAACIDwBAAAEQHgCAAAIgPAEAAAQAOEJAAAgAMITAABAAIQnAACAAAhPAAAAARCeAAAAAiA8AQAABEB4AgAACIDwBAAAEADhCQAAIADCEwAAQACEJwAAgAAITwAAAAH0zHYDACBXNUesth88qbrT51TSv7emjBmoHiGT7WYByDLCEwAksGnPMX3nhd0adeZ3KlG96hTWoX5X6Zu3XakZVwzLdvMAZBHhCQBi1DWe0893HdHOyme1oXCthhfUtew7+mGJHlk7V+9Ov1u3XzNCJUW9s9hSANnCmicAiPGT197RzspntapgpYaOK5PuqZIWH5HuqdLQcWVaVbBSOyuf1U9eeyfbTQWQJYQnAIhxeXGTHi5cK42frtCcddKoyVJhP2nUZOf1+Ol6uHCdLi9uynZTAWQJ4QkAYvTb/ayG2zqFbn5QCsX9FRkKKXTzAg23J9Rv97PZaSCArCM8AUCM/kUDnR9KLkt8gLu95TgA3Q7hCQBihOr2OD/UvZn4AHd7y3EAuh3CEwDEGHb7o/pj4RBFtiyXIpHWOyMRRbas0JneQzXs9kez00AAWUd4AoAYg0eW6iP/a7lMTaUi6+ZIh7ZLTaelQ9sVWTdHpqZS/W77Bw0eWZrtpgLIEmOtzXYbOpQxpkhSQ0NDg4qKirLdHABd1d4XZF9+SKb+3ZZNNjxaZtqj0qTbOuwyVDEHckNjY6OKi4slqdha25jsWIpkAkAik26TmXir9M426cwJqd8QmdE3SqEeF4+JNLfar/j9KWzac0yPvLhXxxrOtWwbVtxbD396UksVc8IVkHsyPvJkjAlLmi1plrW2wuc5CyXVuy/D1tplAa7HyBOAzEs4MnWpzLTHLo5MJQlXm/Yc073PVcsooimht1oeAbMjMlERhbTqrjJJShmuOkN7A1x3CYCpPqeffugufZWLgow8ZTQ8GWPKJF0nKSzpDmvttT7OWShJ0cBkjCmXE7zm+7wm4QlAZu19Qfb5u2UnTFdo6oNO+YK6NxXZ4qyVMrOdGlBe4erY8Ap9+kev6NoPXtHDhWs13MY8AsaU6JGmudpWcKNOn7ugUJJwNeOKYSlHv3z9wr5wQW+9Xqmzp46oz4ARmnj9dPXo6UxMtPcZf35G11K1wZdU/eDj/dsTflJ9Tj/94LevkumIANdd5Ux4armIMTMlLfYZnk5JGmOtrY/ZZq21vv7rEp4AZNIfDteqz09mqM+l1zgVx2MLaUYiiqybo6bfv6beHzZ4hqtnRn5br9a+r1UFK51K5je3PkY1lbr3/P2S5Bmudva9SZtnNGjAK49IMQFN4Usld/TLT/DZVfmMhr/+mIZETrS8xYnQENVc/TXtLf6EdlY+69mGa91n/A3q2zNhMPEzujbjimGebTh6/UO6ZvpnJaUIPylGAf28f6q+8tp/X/lENZz9UN996a02QXe7+zn/bOIl+vVb73nu//otE1Xcp5e+9v92K/43cvQXX7SvMh3gpMyPNObqCFyXDU/GmFJJB+KDkjHGSqqw1lYlOKdQUmHMpv6SDhOeAGTCm0/8pS47WeU8827U5LYHvPOq9MytsmMrZDzCVf3BajU1S0PGlXkGsPqaVxVWo2e4WvPhrZrf69+kCdNlYgKa3bpc2lepzVcs0/+tPpI0+Exq+E/dWP1A0mvM6/VvSQNe2aVh3V63qk0w2XfVIv3t7kuTjq7t7HuTvn/lu/rTXQs8r/G7G1ZKknf4GTUg6SjgwQl/o4/uezrp+58YMU0/X/vP+laCdn67aa5G33SH3nnlZwn3P9I0V5WRKZoe2u75Of3s79Orh85+2OwZsAb3K9CCigl64j/2twk/X/3U+JYA5+Uz14zQ/9t1pM32+HDW3nV4nTUCl4kRtq4cnsolbU4Qnk5J+oK1dmOCc/5e0sPx2wlPADJh7z/dpUl1LzoPCy7s1/aAfZXS2tne4erQdukpd/lnugFs7Z2ytb9WaNwnZe5su9+un6NT+7YprNOeoeFL57+ihwvXewc49xpm7CfbFfCSha+UbVg3Rw37X1Wx9brGJp0rGKTeoyen/RnqDuzS8sj/1jJ9P+3P0d79956/31fAkpRwGre5A6oORcPZ1/91T9LRL8l7HZ4k3ftctef5824eozVbDnruf/wzV6YMgdFRuh/8e027R9ji5WN4OiBpqbV2TYJzGHkC0Gk2/+u/qOJ393sHn988Lv1miXe4ajotLRnp/JxuANv+pPSrB9s1+nWy5nUN1qnMXcNHcEnZhlTXeGqaQkd2pP8Z3CD7vhmgAeOuTy+AtXf/ujk6vn+XvtN0p/6x4Im0p3G9wlV05MqPwp4hNV2IeO7v37unTp+70Ga7kWST7I8/zkufXiGd/dD7+sn4DWgtawUTyEipAjcA3eHj0CXW2mq/7+uT50OkrLVNkloeb24MC98AZE7l+St1uSnR0C3LE/8yrKlyflXVvZn4F3bsY1+8jjm6y/nu9Xy9nr2T729qlCLNMkkebjy4piKz15hQIe2vlNrThlTXGF8uHdmR/mdwtw+yp7zbmepztHf/zQs0vKZCjxY+K42b3vrP1ChnRC2ybo6W1PzYGeUbN915LzdcDd2yXKtqVnZIuEoWnCR5BiObYn/8cV7SDU6x750oOMXu/8bP9+iqkWENC/dJ+1pSgArj1tqN1tpZPr7aE5xqPbaHk+wDgE6zYPokfd/cLXlUINeRHTppwkkf73JUJWrsdYmzPinBMXa/u7zT6/l6F84l358qfMVuz9Q1fAaXdl1j+DXJz0/1GWK3p/s52rs/JsCFvALWTV/VwNAZmfFuuBo12RmxdMOVxk/XkoIfa1XBSg0dV+aMtC0+It1TpaHjyrSqYKWmh7Zremi7XunzgNYXPKonCn6k9QWP6pU+D2h6aHvitnVBqQLaH86c1/KX3273dXLq8SzW2lpJ9e7C8fh9bRaLA0BnGxbuo0995vP60vn7dXx/tbN+aclI6akKHd+/S/ee/6q+be/xDlc1lfp+6G41T39cZl+l7PrWx9j1c2QO75D6DPQOVzWb1Rwq8A5oNSnCl7v9fTPA+z32tfMaPoNL0jakukZhkRTqkf5n2LJC7yvcvs/R3v1+AlyqEbgOCFe399558S0V0cdDe3VbaJs+HtqrkNIfEcpFYy9JMFUeUGeFp4TTbsaY0mhdpxhLJJXHHDNTUpu1TgCQLTOuGKbb535Rs3r9o+48/w3dd/7LuvP8NzSr14/0F3O/qBkzv+AZrr50/n596jOf14DrZkqzn5Wp29vqGFP3pjT7J9KnV3qGK9W8rCebprVr9Kux9wj91/i/836P/S9r94g7vZ/xd2SH6kNJruEGF88AuHWFmosvVfWkr6X/OV75gU5G+iX9DM9EZiQNsm9e802d6TM8/QC2b7MiyT5nzWapR/J+UO8Bzut0R+A6IFx9q2CtQopkdHQq1aKaAX17aUDfXu2+TioH3jvT7vfIdJHMUknRtVJlkpZJ2hG9a84YM0/SImvt2LjzFuriNN1ka+2iANekzhOATpHqlm1fxSWTFXfc+4L08kNxdZxG69RN31LFpmKPMgBD9EjTHBX17qllkRWy46crdPOCmAXIKy4W8vSogXQ8NFTHrv+6UwMpyTP+dh06patevd9d5Nz6Gqqp1O8n/I1K9z0tO2G6zNSL++3WFTL7KqUUbdh/9SJ9cP6Cyvcs9LzGMyO+rdcOvp+wH77dNEe3z/2ihhx5OfVnfP5u77668Suy236YdL+2/dD7c6baP/NfpKpvSiWXS3eubbOOTk9Pkw4nWRif6iYFn3eA7hz1OV1z6F88F61/NfK3euH8dW3Pd0UXjMcvDI9fzC2P/dG7+e59rjrhMR2xKH1wvwK9+OWbEq55yrm77ToT4QlAruiQQn8e4SpagDKkiCZ7FaAM7fD1cOOU1beTBDw/4StRAFSANqS6hp+gmvIzpnoQtI/9ST+nn/3P3y1NmCFNfaAluGjr96R9LznTuKOmJC5N8fQ0Z6q3vXeAFnxE9qNTPctfnDuyR1e8/7gkef6Zk5I/Uqi9dZ4k73Al+QtoHXG3HeEJALooX/Vs2vnwYj/aE7466hqZDKqdtj9ZwJIktxho29GrTVKfgdKoKVKC4JNy5MpP7TH3mJoxf6Wid15OXg0+wxXGO6qSejzCE+EJQDfBs8ryTBrTuLHhKtnIlWe4Wj9H2v/vUuRDX/XJ7IQZbarax07DdoZuUWG8MxGeAAB5KyPhapN07eeknU+nLFiqcRUeAWyuVLdXum9Xh49sdhbCE+EJANAdpRuuJt4qPXF1+ovWo1N/s5+VJv55xqeKM4HwRHgCAKCtVOEq2dSflHpab9hV0tlTcQHtUmnaY502pZcuwhPhCQCA4LxGp/oNkQ5v97Ho3Ljha0FM+FrhTA124pqodBCeCE8AAKQn0ejU6ePS//mUM7KUaM3Tujul2l9LYz/ZZddEBQlPOfV4FgAAkGWhHtKYqdKVM53voR5S8QjplmVOwc31c1tVa9f6uVJNpdR8XpqauMq5pj4g1b8jvfXL7HymDtYz9SEAAKDbm3SbM/X28kMX60JJzrTe0Kuk479L/SDoXc9Jk/5X5tuaYYQnAADgz6TbnDvz4qf1NvyNE57q3ky8Jir63L6efTu3vRlCeAIAAP5Fp/VizfiudPh1aevyxGuetq6Q+g93jssDrHkCAADtk2pN1L5K6ZalznGRZungVmn3Rud7pDnbrQ+Mu+0AAEDHSOsByblRB4pSBYQnAACyw6sQZ6sinLlXB4rwRHgCACB3RJqTP/4lB+pAUecJAADkjrd+6UzVTV2QF3WgCE8AACCzdj3nfPdTB6oLIDwBAIDMitZ3itZ7itfF6kARngAAQGbN+K7Uf5hTByoSab2vC9aBIjwBAIDMClIHqgvgbjsAANA5UtWByqIgd9vxeBYAANA5vJ6Nl6XyBOkiPAEAgM6T6Nl4XQxrngAAAALI25Gnxsak05UAAAAtguSGfFwwPkLS4Wy3AwAAdEkjrbVHkh2Qj+HJSBou6XTM5v5yAtXIuO0Ihn7sOPRlx6AfOw592THox46Tjb7sL+moTRGO8m7azv3ArRKjk6ckSadT3X4Ib/Rjx6EvOwb92HHoy45BP3acLPWlr+uwYBwAACAAwhMAAEAA3SU8NUl6xP2O9NGPHYe+7Bj0Y8ehLzsG/dhxcrYv827BOAAAQCZ1l5EnAACADkF4AgAACIDwBAAAEADhCQAAIIC8K5IZyxizUFK9+zJsrV2WxeZ0GcaYsKTZkmZZaysS7KdffXL7SpLGSpK1dn6C/fXuS/rSQ8yfScnpy1JJX7DW1sccQ18GZIzZHP//OP3ojzGmXNJ8SZsl1UqqkLTDWrsx5hj60idjzFJJB9yXJ3O9H/M2PEV/aVlr17ivy40xq+N/eaE1Y0yZpOskhSUNTLCffvXJGLPUWrso5vXq2F9W9GUgSyUttdbWSk5fStog5xcWfZkGY8xMSeVx2+hH/8Jy+m+mnPC0NMEvfPoyBfcfRv8u6VPW2nr3d9BOScbdn5P9mLelCowxpySNifuXqbXWGu+zEOX+xbrYWntt3Hb61Qf3L4QNckbv6t1t0b8Uxlpra+lL/4wxmyVtjv6L0/0LdbG1doD7mr4MIGYkb3VsH9GP/rl/R1bF9lXcfvrSB/cfQgdiR5OMMeXW2ir355zsx7xc82SMKZUztFefYF952zPgB/0a2HVyppeiat3vYfoyGGttRdxQ/WRJ0b9c6cvgZkt6PnYD/dhx6MtA5knaaIwpjfZNTHDK2X7M12m7Uo/t9XKGWpEe+tUn93/2AXGbo/+z18oJVonUi75Myv0Xf1jSLHcTfy4DcH/pVCXYRT8GN9sYc1LOEoexMdP09KUPbjiSpDI5fy/WRqfk3QCVs/2Yr+HJS/QPOToW/erPYknz3Xl9r2PoSw8xU01hOX+51qc4hb5MLOxOG4d9Hk8/JlYtSTHr8OYZYzZYa2clOYe+bC0ajuqttdWSZIxZJOmg2v7jM1bW+7G7hSf+0GYG/ZqCeyfJz6KLHpOgLz24YSm6aHRedC1EklPoyzjGmHk+/gzGox8TiIamGM9LWp0ilNKXib0R/cH9x2U4xbRc1vsxL9c86eLaknjhJPuQGv2aBnea6UDcmh360if3L9Klcb+UqnTxbif60gf3hoU3khxCPwbg/n/dImYktFT0pV9efVGvHO/HvAxP7r8I6mPmU2P3JZrrhw/0a3AxCyCjIyZhY0wpfRlIqaSFav2vzbD7vZ6+9G2gpHJjzEL3bsWlknPnojFmJv3oX/Ru2ti+ign3tfSlP24/1art2qawpDdyuR/zMjy5liimhon7r4Sgw9XdmdewKP3qk/sv/TJJ1e6dJKVy7iw56R5CX/rgroVYFjdNcoek6pi/QOnLFKy1VdbaZdEvSavd7cti6hPRjz64o0zxfybnSdoYMwJFX/qzSM7/z5JalYCodjflZD/mbZ0nqaUWTPQP9+TYgoVIzP0FP1POH+YyScuUuGou/ZqE+6/Qg0pwR0hcXR360ge3P+fFbBoraVGCCuP0pQ/uL6A75Py/vkxODa3o7eH0ow8J/kwOiu8r+tIfY8w8Xfy7skv0Y16HJwAAgI6Wz9N2AAAAHY7wBAAAEADhCQAAIADCEwAAQACEJwAAgAAITwAAAAEQngAAAAIgPAGAWh6dE852OwDkPsITADgWq+0ztgCgDcITADjKYp6nBQCeCE8Auj1jTLmkzdluB4CugfAEANIsSRtTHgUAIjwBgCSVWmtrUx8GAFLPbDcAAPwyxpRJuk7SWEk7JFVJmufurrfWrknjPWdK2pBk32RJByTVul8nrbX1gRsPIG8w8gSgS3DLCJRba9dYaxdJelLSYmvtMveQRWm+9R2Snk9wvXmSKqy1i9xQFpYToq5L8zoA8gQjTwC6inkxQSnqgPu9WtL8NN83HD+SZIwplbRU0piYzfWSZK2tSvM6APIE4QlAV9GyoNsNN2G5I0bxgcbdP1PONNtkSasTrWlyR5dWJ7jWaklVcaGqQk5IA9DNEZ4AdAlx4adcUm2StUcbrLXXSpIxpkrSv0u6NsFxs6y1FQm2l8u5Ay9WmZw1VgC6OdY8AeiKKhRXWiD6aBV3UXkLN2CF3dGo+OPr49845rj4USZqQQGQRHgC0EW4U2xRM+XcbdeyL2YUymtBd1nca68pO0mtR7rcIpqy1lYZY8riAxqA7oXwBCDnucFpqfvzTMVMnyV4mG9Y0sm4bfWSBsZtq0i0+NsNTbXRgOS+/3w566ck544/1j4B3RhrngB0BVWS1rgh6g05YWaRMUaSBsbVd6pX26AUVkygcqfmkhXFnCVpvjFmpyRZa2cZYza41yc4Ad2csdZmuw0A0GHcEaMnowvG3W2nJF0bnYozxiyV9DNGkACkg2k7AHnFDUTh6Gt32q027m69MoITgHQxbQcgH81yR5d2yKnz1FJ2wB2ZIjgBSBvTdgC6FWPMaklLeRAwgHQxbQeguxlIcALQHow8AQAABMDIEwAAQACEJwAAgAAITwAAAAEQngAAAAIgPAEAAARAeAIAAAiA8AQAABAA4QkAACCA/w9wCFnEsOly+AAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEfCAYAAABrrED8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAhiklEQVR4nO3de3hV5YHv8d+bC0ko6CYShgMYLs6pPBXGCqEVx6F1hPbUtoiw8cJxxrbjBFOllxmlZTrYIm0xyDm1YxHh2HkeR0dBttxq+xwLnrapT2NLiO3BzmB7DCUI44DiLmhCyOU9f6y1w76sDQlr7Wu+n+fhgf3ud6/1vgb3j/ey1jLWWgEAcKFKct0AAEBhI0gAAL4QJAAAXwgSAIAvBAkAwBeCBADgS1muG5DO6NGj7aRJk3LdDACApH379r1lra3xei9vg2TSpElqaWnJdTMAAJKMMYfSvcfUFgDAF4IEAOBLIFNbxpiwpGpJMyVttdbuSVMnKikkqc1a2xrEuQEAueU7SIwxM+QEQ8R9/Y6kUUl1QpLmWWuXuq93S5rn99wAgNwLYmqrWtLSuNcn3HCJd7Ok1+NeRz3qAAAKkO8RiTuNFT+VVe0xbRWSM60Vc0JOACUwxtRLqpek2tpav00DAGRBoIvtxphGSYsHWD2UXGCt3WStrbPW1tXUeG5XBgDkmcCuI3EX07ekWUSPKjE4qiW1BXVuAEDuBDIicdc7Wq21rcaYKcaYKUlVnpV0WdzrELu2AKA4BLVr60U5i+ySs0Yyyn1vn6TrrbVRY8xuY8xc92ONvk66PyI1rZPeek0afbk0515petjXIQEAFyaIxfZWJW33jXtvZtyfI37PJckJkRcfkG78nlQ7W2pvlnbe47xHmABA1hXele0/XumEyOQ5Umm58/uN33PKAQBZV3hB8u6bzkgkXu1spxwAkHWFFySjL3ems+K1NzvlAICsK7wgmXOvsyZysEnq7XZ+33mPUw4AyLq8fR6JJO145YgeeuE1HY12alyoSvd9/HItuMpdUP/R8rO7tq6/n4V2AMgRY63NdRs8/ekHrrTlixrV2d3bX1ZVXqo1C6drwVXjc9gyABh6jDH7rLV1Xu/l7dTWkWhnQohIUmd3r1bufDVHLQIAeMnbIOlLM1J693RPllsCADiXvA2SEucq+RQjKvN6WQcAhpy8DZLxoSpVlZcmlFWVl2r1jdNy1CIAgJe8DZLQ8HKtWThd40NVMnKChYV2AMg/eT1PtOCq8QQHAOS5vB2RAAAKA0ECAPCFIAEA+EKQAAB8IUgAAL4UT5Dsj0jrr5ZWjXJ+3x/MAxkBAOeW19t/B4zH7wJAzhTHiKRpnffjd5vW5bplAFD0iiNIjh/wfvzu8QO5aQ8ADCHFESQ1U70fv1szNTftAYAhpDiChMfvAkDOFMdi+3QevwsAuVIcQSI5oUFwAEDW+Z7aMsaEjDHLjTH156iz0Rgzw/213O85AQD5I4gRiefD4JNMkbRV0h5r7dIAzgkAyBO+g8Rau8cYM+U81RoltUiq9ns+AEB+ydaurSlyQiRkjGlMV8kYU2+MaTHGtBw/fjxLTQMA+JGVILHWbrLWtllrWyXNTTeCcevVWWvrampqstE0AIBPGQ8SY0w4aYH9RKbPCQDInmxs/21zf8VUW2vb0lUGABQW30FijJkraZ6c9Y82a+0et3yfpOutta3uqGSKnLWSxX7PCQDIH4Hs2pK0x6N8Ztyfc/dwkP0R5y7AsSve59zLhYsAEKDiubLdC88pAYCMK46bNqbz45Xezyn58cpctwwAikZxB8m7b3o/p+TdN3PTHgAoQsUdJKMv935OyejLc9MeAChCxR0kPKcEADKuuBfbeU4JAGRccQeJxHNKACDDintqCwCQcUM3SPZHpPVXS6tGOb/vz901kwBQyIp/assLFyoCQGCG5oikaZ33hYpN63LdMgAoOEMzSI4f8L5Q8fiB3LQHAArY0AySmqneFyrWTM1NewCggA3NIOFCRQAIzNBcbE93oaLk7ODilvMAMGBDM0ik1AsV2ckFABdkaE5teeGW8wBwQQiSmHPdcp6LFwEgraE7tZUsdsv5yXPOlrU3SyP+C1NeAHAOjEhi0u3kKi1LM+V1P6MUABAjkrPS7eR67s7UKa+TR6W+bumGtYxSAAx5BEk8r1vON61LnfL66YNS+Ptny2KjlB8tP/sZthADGCIIkvOJTXnFr5G884f0t1jxWk85/Evp4M9Tw2V/hNABUPAKMkh2vHJED73wmo5GOzUuVKX7Pu48gz25bMFV4/2fbHpYe/9wQpdu/oLGdB3SsYqJGlVVowqvhfmKkWfXUyTn9+mLpVeelBY9nhouv3thyIWO188ukJ8TgJwpuCDZ8coRrdi2X53dvZKkI9FO3bf1N5KRunttf9mKbfv7P+P1xTXQMJKkFb+qVWf3t5yDnZbCw17Wt7fdrWEL1/eHQM/2u1XWdSp1pHLgeSdEkqfBNt8u3fpUSujYV56UiQudnu13q+zwL9Vz4H+r7Kak88XOkRwwXmWDDKO9P9ikS199tD88D0/7vGZNqvb8vGfdT9ennG9v7ef00t7D2lKxS+Mq23W0t1YPb58v6S4tKGsecD+CCCMCDQiOsdbmug2e6urqbEtLS0r5nz/4f3Qk2jmgY1SUGZWYkv7QkaSq8lItmjlez+07klBeXmISwihW11qr0z19KceeX/oLLR/+vMZ1t+toea0e7pqvr130Q41a9HDiSGXVKOkfjzk7vmJ6u2VX18isPJ5Yvv7D0g0PJX7+YJN6N/+VSm99MqW869m/VWn5sISA6X32syodVikt2JA40nn/xz3D6NAl1+p9R5sTQkCSph34rqrCZ49xZvMdKimvUtmix1I+P+5YU0LdzkiDjo6Zo4lvv5RwvjOb75DKKjUsvLG/rCPSoJ3vTdNtoX9PHJ1t/axUntqPvZfdoy17D+tLFbsS/tt/aPIozXnzX1LCzCvkjkz4pF7a/ljKMa696S7CBEjDGLPPWlvn+Z7fIDHGhCTVS4paazelqROWFJUUktRmrW0933HTBcmkr/7QR2sz68bSX2hNaIeGx32p9m2+XSXxIw9JOtiknmduV9ltT1146PR2yzZOkrntmcRjfPeD0vx/Sjlf35a/UsktSWH04mr1tT6pkvDjCSFwurtXo27bNKDjevYjXflgjpGm7jtPfU4VlRUJ/427Nt8hkxRQnZEG/brqQ/pg569SQu65d6/QwhG/TThGf6DVHMrICAgodOcKkiCuI/E8cNzJQ5LmWWv3WGsjkhr9nGxkZf7Oxu3svUZfjS7QG08vU9/qMXrj6WV6suNqdUQaEq5P6Yg06F87U8v7ho30vL1977CLvG97f+a91Km06CHPjQDm9B89p91Kwo8nXCNTFd6g0BmPq/zTHLf0zMmBlw/mGGnqhko7nQCIa3PF8IudEEnqx6yTLzohklR+a8XLKccYPuNW3fK+V5wt3f94TLphrXp2r9LeH2zSS9sf05beL+v1ytu1pffLemn7Y9rxyhEBcPgOEmvtHjmjjXRulvR63OuoMWbGhZ5v9Y3TVFVemlBWXmJUXmoSyqrKS1VZ5j8nQ1XlKec7l1291+jaU9/WlNNP6dpT39bXu+9ICZevRhd4lg8mdDoiDXqnfGxqwIQmDjyM3vqd9+6zYe/zd9x05YM5Rpq6gwnPQYWcR6iW3bRef/rrRj0wcpsmLHlEJSuPacKSR/TAyG1q3vmYADiycWV7SIlBc0JStVdFY0y9MabFGNNy/Phxz4MtuGq81iycrvGhKhlJ40NVemjxlXoofGVC2ZqF0/Xgoj9LCYGq8lLdfnXtgMPoG/OvSDmf1+fTHdcoNVx29V6jUFW5dpf8RUL5N3s/o6+dTAyX+08t1O9nrtT9pxamlH+rK5wSMF0df9SZyNKU0Hl6ECOgd3qr1DnA427uujqlbqdbPpC29Wy/27Nt6c7nOzwHGaqh0tOpo5fwBi0r3yUAjkAW240x9ZLktUZijFmuuPUTY8xGSbvdaa600q2RDFa6+W2/W4gHetzrptakLOxXlZdqzcLpnudL14Z07U1eNF7b8SlJStgIsLbjU9rVe43Cw15OqPuTrsu1aGTiWkFnpEGvTv2iJCUsUjeN/Wv96uA7ngvU49/44YAXtG+Zdalmtf9zwjrE9OdG6rruppQ2e/WjqrxE37poR+KmgS2fVV/psJQ1khd7Pqi/LPt1ylrIy6V1mqPWhGOkW8uyT8z3XJ/qWz1GJd94x/ffT6BQZHSx3T3BuYKkXlLIWrvWfb1V0przLbgHFST5IJOLtcnH7jjTo3c6ulPqjXfPm9wOrxCY9en6QPox0PrJW7ql9Lvo1iyc7rlVeO8fTnj2Y6Bbk587MVmfqPy/KaFz+kyvqpdsSgmYY1u+oJv0P1mAx5CR6yAJSWq01i51X++21s473zGLKUiyyetLuf8LOI+/6LJ6kamHpf/SoorXtqeMgP5s/MX6zOknE7cxb7tbW09+QB8ZdoDtwxgyMr39d66kpXLWQhrdxXcZY/ZJut5aG43b/iupf4H+nAiSC8d21QuT9r9b0ujlmeMTdeP7Xk0Zvax6b6Eav7E6190AMiLjI5JMIEiQr9741pWasOSRlOmuN55epglf+03uGgZkUKavIwGGlHHd7Z47vMZ1t+emQUCOESTAIJ0ePs5zW/Hp4eNy0yAgxwgSYJCGf2KVerbfnXI9zPBPrMp104CcyN/7jQD5anrY+R8n7mmaZfO+XpC39QeCQJAAF8LraZrAEMXUFgDAF4IEAOALQQIA8IUgAQD4QpAAAHwhSAAAvhAkAABfCBIAgC8ECQDAF4IEAOALQQIA8IUgAQD4QpAAAHwhSAAAvhAkAABfCBIAgC8ECQDAF4IEAOALQQIA8IUgAQD4UhbEQYwxYUlRSSFJbdbaVo86GyVtdF/OtdauDeLcAIDc8h0kxpiQpHnW2qXu692S5nlUnSJpq6Q9sboAgMIXxIjkZkmvx72OGmNmeIxKGiW1SKoO4JwAgDwRRJCE5ExrxZyQd1hMkdQmKWSMabTWfiW5gjGmXlK9JNXW1gbQNCB7mr9/r8oO/UyXVnRqTPcRHSsfr8NdVeqZ+BHN/pt1uW4ekDGZWmwPJRdYazdZa2PrJ3ONMVPS1Kmz1tbV1NRkqGlAZpSNeb+mjezQ2CWPqmTlMY1d8qimjexQ2Zj357ppQEYFESRRJQZHtZyRRz9jTNgYszyu6EQA5wXyyqWvPqqq8AZp8hyptFyaPEdV4Q269NVHc900IKOCCJJnJV0W9zrksT7SJmlP3Otqa22bgCIypuuQVDs7sbB2tlMOFDHfQWKtjUrabYyZa4yZK2dRXZJkjNlnjIkFy5S4kcliv+cF8s2xiolSe3NiYXuzUw4UsUCuI7HWRtKUzzxfHaBYPH/xEi2JNGh4eIMzMmlvVkekQc9ffIfuzHXjgAwKJEgASHd+frn2/iCkSzd/QWO6DulYxUQdnvZF3fnp+lw3DcgoggQIyk/WaNa+B/tfju06qLH77pNGvC1dtyKHDQMyiyABgnLdCgIDQxI3bQQA+EKQAAB8IUgAAL4QJAAAXwgSAIAvBAkAwBeCBADgC0ECAPCFIAEA+EKQAAB8IUgAAL4QJAAAXwgSAIAvBAkAwBeCBADgC0ECAPCFIAEA+EKQAAB8IUgAAL4QJAAAXwgSAIAvBAkAwBeCBADgS1kQBzHGhCVFJYUktVlrWy+kDgCg8PgOEmNMSNI8a+1S9/VuSfMGWwcAUJiCmNq6WdLrca+jxpgZF1AHAFCAggiSkJwpq5gTkqovoI6MMfXGmBZjTMvx48cDaBoAINMytdgeupA61tpN1to6a21dTU1N4I0CAAQviCCJKjEUqiW1XUAdAEABCiJInpV0WdzrkMeOrIHUAQAUIN+7tqy1UWPMbmPMXLeoMfaeMWafpOvPVQcAUNgCuY7EWhtJUz7zfHUAAIWNK9sBAL4QJAAAXwgSAIAvBAkAwBeCBADgC0ECAPCFIAEA+EKQAAB8IUgAAL4QJAAAXwgSAIAvBAkAwBeCBADgC0ECAPCFIAEA+EKQAAB8IUgAAL4QJAAAXwgSAIAvgTyzHUAaP1kj/ezB1PKPfFW6bkX22wNkAEECZNJ1K6TR/1VqWicdPyDVTJXm3CtND+e6ZUBgCBIgg3709CP66JHHNDy8QaqdLbU3qyPSoJ/u/w/dsGRZrpsHBII1EiCDPvqfTzghMnmOVFouTZ6j4eEN+uh/PpHrpgGBIUiADKr84+vOSCRe7WynHCgSBAmQQUfLa6X25sTC9manHCgSgQSJMSZsjJnr/j4jTZ2NxpgZ7q/lQZwXyHePa6E6Ig3SwSapt1s62KSOSIMe18JcNw0IjO/FdmNMSNI8a+1S9/VuSfM8qk6RtFXSnlhdoNh98IY7df/2Hn3p6WUa192uo+W1erhroa696c5cNw0ITBC7tm6WFD/hGzXGzLDWtibVa5TUIqk6gHMCBeFPWr+ja/t+q1tOfVlHdYnGnX5b95Vt0Z+0fke6al2umwcEIoggCUmKxr0+Ie+wmCKpTVLIGNNorf1KcgVjTL2kekmqrWUOGYVv9t84YbEgofQz2W8IkEGZWmwPJRdYazdZa9vckcpcY8yUNHXqrLV1NTU1GWoaACBI5x2RuKOEkNd71tq1ckYj8e9Xyxl5xB8jLGmKW19yRi0AgCJw3iCx1m46T5Vn5ax/xIQ81kfalBgu1dbaNgEACp7vNRJrbdQYs9sYM9ct6g8VY8w+Sddba1vdrcFT5KyVLPZ7XgBAfgjkXlvW2kia8pnnqwMAKGxc2Q4A8IUgAQD4QpAAAHwhSAAAvhAkAABfCBIAgC8ECQDAF4IEAOALQQIA8IUgAQD4QpAAAHwhSAAAvhAkAABfCBIAgC8ECQDAF4IEAOALQQIA8CWQJyQCGLwdrxzRQy+8pqPRTo0LVem+j1+uBVeNz3WzgEEjSIAc2PHKEb20/TFtqdilcZXtOtpbq4e3z5d0F2GCgsPUFpADzTsf0wMjt2nCkkdUsvKYJix5RA+M3KbmnY/lumnAoBEkQA4sK9+l4eEN0uQ5Umm5NHmOhoc3aFn5rlw3DRg0ggTIgXHd7VLt7MTC2tlOOVBgCBIgB05VjpXamxML25udcqDAECRADlz8ydXq2X63dLBJ6u2WDjapZ/vduviTq3PdNGDQ2LUF5ML0sPM/34+WS2+9Jo2+XGXzvi5ND+e6ZcCgESRArkwPewfH/ojUtK4/YDTn3rQBw7UoyAeBBIkxJiSpXlLUWrspTZ2wpKikkKQ2a21rEOcGisr+iHp2r1LZTeudxfj2ZvVsv1tlh38pHfx5Qrjs6JmtFdv2q7O7V5J0JNqpFdv2q+XQCf3kwPGMhEuxBJdXPyR59q1Y+pxJQY1I6s71phs086y1S93XuyXNC+jcQNHo2POght+03tkWLEmT56jsyptlX3lSZtHjCeHSfHKB5vX1afnI5zWuu11Hy2u1tuNTOrlX2jL8+dQLHcuaU0Y6O3pme35J7v3BJl366qMa03VIxyom6vC0z+vIhE8O6iLKwXxZe0ozMvNq26xP13seIl0bkgP4vq2/kYzU3Wv7y2Kh/Ny+IylhLWnAfR4KYWSstcEcyJh6SfIakbjvhay1a93XWyWtOdeopK6uzra0tATSNqBQ9H1jlEpWHnOuLYlZ/2HphofOhoskHWzS2099TlWVFc71KG7AdG2+Q6asUsPCG/vLOiIN2vneNC2+6N8SRjpntt2trSc/oI8MO9AfRA93zdeHJo/Sp9/6vqrijtsZadBz716hhSN+m3C+jkiDtr17heZW/s4zdL5Usav/2P+j81OSkf6+8vmE8117010a/8YPE8LhvXGzNfHtl1JGZocuuVbjjjWltO3VqV+UpIRjNI39a/3q4DsJbXi4a76e77tGH7MvafnwxACWlFK2q/cazS/9RUr5C+bPNXpE5TkDSpKqyku1aOb4hDCKLx/oyHEwASUNbmQ10JAzxuyz1noOGrIVJMsVN+1ljNkoaau1do/HMeolqba2duahQ4cCaRtQKN741pWasOSRxNBYNUr6x6Rw6e2WbZwkc9sziXW/+0Fp/j+lhE7PM7er7LanEstfXK2+1idVEn48IRhOn+lV9ZJNvo7x3KkrtGjkbwcUctvevUKLRvw2IRz6Nt+uklufGlgbDjbphBuqVQM439aTV2jxRQNrm1fdjkiDvhpdoF291/Q3oaq8VNZane7pG/wPPe4YaxZOl5QYBNdNrUkbRMnl5SUmYWR1rrrnKl+zcHpKmJwrSHK5/TeUXGCt3WStrbPW1tXU1OSgSUBuPdI9Xx2RhoRtwX3DRnpec6Iz76Ve1Bg95HmhY+mZk6nlB553AiDp6vpR3W/6PsZtlS+nXLlfMfxi54s6qe6tFS87ARBXXnLm1MDbUDtbo0o7U46R7nz/vWrgbfOqOzy8QV+v3KKXRv6D2ipv10sj/0Hz+n6u0z19ml/6i4Ty+aW/OP8P3dXZ3auvPPcbrdi2X0einbJyptKeerk94Ys+VtervLvPJoRIrO6/pjnGM7887Fn+0AuvDbjd0gDWSGLTUl7vxaaqBiCadIxqSW0D/CwwZMy+8S7dv1360tPL+qdSXjx9tRZHGlL+Vdw1bKxGtTcn/us8NNEJmfiy9mb1DrtIZcnlb/3O80tZw97n+xieX/iDCbnR7x94G4II1cHUPXlU1ZVGJvyIVDtbE9qb9WCkQTNP/j5h9BIrV9T5mNe0WbKuHiupN6Xcr3TzTr1pZqSORDsHdfzzjkjcUcJar1+DOM+zki6Lex1i1xaQasFV43XtTXfpltLv6LLTT+mW0u/o9zNX6v5TC/XG08vUt3qM3nh6me4/tVD/b/rfSTvvSRi9qOuUerfdlXKh49Yzswc80on2VqXU7Yg0qKnyupSLKNMdo3fYRanlsZAbSN2pn1Jf5M6UNmzuulqdSW3rjDQoOszjTgFpzmcrLx5wXVPh0b+fPigT/v6ARy/fHrFVjaEdCTfofDC0Q6vKn7jg0YskmQzVlaSRlYPbhxXIGokxZq6kpXJGHY2xtQ9jzD5J11tro3HbfyVJyesjyVhsB85KuyDqtbNJ8tydlbz4/bMzU7X4on/TsIVJW43nfV17/3DCe2dU0vleH3GV5+L30TFzUhbLe5/9rFRWodKFj513AT22znJdxWvnXZg/PO3zkqRpB76bcIwzm+9QSXmVyhYlnq9s6n+TfveCdOP3+su19bNSeaW04OzntfMe6f0fT6lrn5gvs/J46prV6hrPcjVOkpLXsjzWljojDXqgY5Ge6fxwys/fKHFUkW9rJIEttgeNIAGC5RlGHluCB3t1fdrtuAMMuXRbeo9M+OSgtsx6tmNStXf/BtG2lLrdHU6wJC34a/PtkscGAT0xX0oOmDQ78Tp2/J1mnvjmgHd4Fd2uraARJADy0v6I9OIDiSOaNKMX7bzH+Uxy8KTZiadvjtHeGY0Dvk4mm84VJNwiBQAGIzZii7tPmq6/3ym/9MOp5ZITKPEBE1t7Sd40UHGxZr3+PelWp+7Y9maN3XmPtL86r+/DxogEADIteXps8l94j166O6TwP6dOj21bKv39v+eu/WJEAgC55XWDTq/Ry7a/9d6Sfeo/pPVX+1rLyiSCBABywStcmtalTnn99EFpxBjphrWpay95EiY82AoA8sWce1OvDdr7v6RFiXcP0I3fk368Mtet7ceIBADyhddC/mnv28Lo3Tez3740CBIAyCfJU17rr/be4TX68uy3LQ2mtgAgn3lNd+28xynfH3GCZtUo5/f9kZw0kREJAOSzdNetSN4XRsZ/JksIEgDId147vNZfnXjFfGwR/kfLsx4kTG0BQCF66zXvRfi3BvcskSAQJABQiEZ43Dq/vdkpzzKCBAAK0cdWey/Cf2x11pvCGgkAFKJz3TwyywgSAChUXovwOcDUFgDAl7y9jbwx5rikQ+7L0ZLeymFzMo3+Fa5i7ptE/wpZ0H2baK2t8Xojb4MknjGmJd198IsB/Stcxdw3if4Vsmz2jaktAIAvBAkAwJdCCZJNuW5AhtG/wlXMfZPoXyHLWt8KYo0EAJC/CmVEAgDIUwQJAMCXvL6y3RgTlhSVFJLUZq1tzWmDfDLGhCTVS4paazfFlRdNP92+VEuaKWmrtXZPXHlUBdxHtw9tkuokKfYzLIa+xYv1p5h+dpJkjNkoaaP7cq61dq1bXhT9kyRjTL2cv6PK6s/PWpuXv9xOb4x7vTvXbQqgT3MlLZdUX4z9lDRD0oy41+8USx/dPuwrxr559HO3nC/aouqf26/Xk/pTTP3bKikU+3M2+5fPU1s3y/mhx0SNMTNy1ZggWOdfCNGk4mLqZ7WkpXGvT7h9Kfg+Wmuj1tqZkmSMmSKpxX2r4PuWpE7OF25MMfWvUc5IuTGurCj6F2uztTZqjJlirV3svpWV/uVzkISU+KV7Qs4XVbEJqUj6aa3dY62ND5Jq6wyjQyqSPhpj5koKS4r9jxpS8fRths4GZExIRdI/SVPktD1kjImFSUjF0b86qX/6POv9y+cg8RLKdQOyJJTrBvjl/kVefI4qoSw1JVDuqLJVzlpXOqHstCZ41troAKqFMtyMjLDWbrLWxtYI5rojSy+hLDYrKCGpf+TcKmlGNvuXz4vtUSV2uFruIlKRiarI+uku7m2xZxf1oiqiPlpr9xhjthpjYlOVobi3C7Jv7iLtCffLZ5akS4wxbSqe/oUlTbHuArucf5lLRdI/OW2+JO51VM4ILKos9C+fRyTPSros7nXIFvBuinMoqn660yOt1tpWY8wU94up4PtojKmPmy6Qzn4RFXzfpP5/rUestRE5XzS7rbVtKpL+yenTnrjX1UXWvz1KDIzYOl5W+pfXV7bHbVuTdHY7W6Fy59eXyvmBN9rU7XmSCrefboi8qLNfstXW2lHuewXdR3fuOXYn1djCZvL2UbnlBdW3eO7PsFHOF2+jtbatWPrn9kNyvmQjbpAUzc8vbuu9JJ1w/1GQlf7ldZAAAPJfPk9tAQAKAEECAPCFIAEA+EKQAAB8IUiALHN3gAFFgyABsm9FrhsABIkgAbLIHY28net2AEEiSIDsullSJNeNAIJEkADZdVnsimqgWOTzTRuBvOPeO2yGnBsbbnT/HH8zwPN99nWP8uU6eyO9pdbaecG1GMg8RiTA4Mx172G0V86T5yKKe5iXe3PHue7ddJOF5dxET3H1G+U8/jQi5/b0hfgsDAxxBAkwOLEgmCXn0aay1l4m9d8c70TczTjnJn32kvjnfbgjlHDs5npyRjcFecNADG0ECTAIcUEQttZukvoDQXLCJfZ+m6T+KSr3rrrxj7CVUoNjnkcdIO8RJMAAGWOWu1NXM+QGhhsisSAJ6ewt9GOvY25Jc/vuaNyfb3YfmhX2qAfkLRbbgYGLPX9ekjbGvvDjpqaiSlzjiJ7rYNbaiDFmVtwU2Br3mIX4YCUMYTyPBAhIfLC4i+1tcSOM2LPCgaLD1BYQEHdkUh0bYcRNZc0iRFDMmNoCAhRbgI/hligYCpjaAgD4wtQWAMAXggQA4AtBAgDwhSABAPhCkAAAfCFIAAC+ECQAAF/+P/SfFH0WnVgqAAAAAElFTkSuQmCC", "text/plain": [ - "
" + "
" ] }, "metadata": { @@ -236,7 +234,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 9, "id": "def98bc0", "metadata": {}, "outputs": [ @@ -244,14 +242,14 @@ "name": "stdout", "output_type": "stream", "text": [ - "--- The mass was calculated to be 3079.8(7.9) MeV ---\n" + "--- The mass was calculated to be 3080.0(7.9) MeV ---\n" ] }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkoAAAGLCAYAAAA4f70VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA1NklEQVR4nO3df2wc52Hm8eelJJOyRHJNpzTZyGq0st2YxxN8tGTXqNsAMXkRmtMVxpEWejZwSIGQMNqD7euBPPkfVTjAOhKHXoy0Z5AGGgSV7yCJhZG6CZSKNprUB8OWxLoCT0ljcZ3KSqjqbHIp6gcpWXrvj51ZzS53dmdJ7s5o5/sBVtT83Hdnl7MP3/edd4y1VgAAAFiuLuwCAAAARBVBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwMf6sAuwWsYYI+lXJS2EXRYAAHBHaZT0S1tkUMk7PigpE5LOh10IAABwR9oi6Rd+C2shKC1I0ieffKKmpqawywIAAO4Aly5d0v333y+VaJGqhaAkSWpqaiIoAQCANUVnbgAAAB8EJQAAAB8EJQAAAB8EJQAAAB8EJQAAAB8EJQAAAB8EJQAAAB8EJQAAAB8EJQAAAB9ljcxtjElIekZSn7W2J+A2g5LSzmTCWjtSYJ1hSdPO5Ky1dryccgEAAFRC4KBkjOmStFNSQlJLwG0GJclaO+ZMdxtjRq21A850QtLbkp6y1qad5zglyZTxGgAAACoicNObtXbSCTypMva/T9KYZx8Tkvo9y4clHbbWpt3nkBSopgoAAKDSKtZHyRiTVKapLV1gWbfz335J48aYpDvPCVPF9ltvjGlyH5Ia17joAAAAkirbmTvpMz8tKeEEKUnqUqY5L2WMGfWEKD/7JM17HudXX1QAALDW0ul0pPazEmFc9TarTB8nNyilnWa9lKQhSUdLbH9QUrPnsaVSBQUAACszNDSkRCKxJvsaGxtTKlVOz5+1U9ZVb2skvyP4Sfc/TofuhDGm268Jzlq7JGnJnTaGft8AgNqUTqd15MgRHT16VMePHw+0zcjISDagpNNpDQ4OLlsuSdPTmYvNR0dHV7S9JH322WcaHh5eVoaxsTENDAyU9bxSJlxt375dktTS0qLe3l5J0uDgoAYGBgpuU2mVrFHyi34JZ5nf8rT8m+0AAIiFyclJHTlyROl0WrOzs4G2ccNIf3+/+vv71dXVlRNYhoaGNDg4qMHBwWzo6OnpCbx9X1+fEolEdh/bt2/X0NBQThlSqZROnTqlZPL2V3mp502n03r00Ue1b98+9ff3a+fOnerr68vZb19fX05IqxprbVkPSb2STgVcd05SMm+e9fx/WlJX/vL8eSWeo0mSnZ+ftwAA1JqjR4/arq6uQOsmEgk7NzeXM8/53rVzc3O2u7s7Z/mpU6esJDs9PV1y++npaSspZ/nc3NyyeYODg9n9BX3e/v5+Ozw8nPO8x48fX/b6gh6HIObn562TOZpskZyxkqa3gmMoOZ2ze23ugJIHJXXLGSLAGNMrz3AByvRJ2itp0rN8wmaGCQAAYFUmJiYCNVlt3bpVf/AHf5Az78/+7M907tw5/fCHP9TXvvY13217enrU3X37OqTFxUXt37/fd3mlpFIppdPpgv2CJiYmtHPnTp08eVKpVEpdXV2SlK31SafTJbd3eZe7/z958mT2NU5MTCxrjiv2vFKmqW56elqpVEqpVErd3d0Fj1kymdTk5GR2P9VQzoCTSWVqk/ZK6nJG0z5hb4+i3S1pQFI2KFlrR4wxg04AkqRd1hls0lk+boxpcQemlHSvDTjiNwAApVy7di3QFVMtLcvrAC5fvqx0Oq3HH3+86D6uXbu2bJ53/ULLK8Gvs3MikcgGoLm5uZxlbgBKJpM6efJkoc2z27vhpFCYcp87lUotO5alntfddnJyUslkUslkUgMDA+rr61sWlnp6ejQxMRHNoGQzV6WNyBOE8paPKbe2yJ3vXX/ZrUmc7QAAWHMbN24MdOXV5s2bC84Lsu3GjRuXzfNuV2h5NbW0tPj2cTp48KBGR0eLvk53+2Qyqe7ubk1MTGQ7WXtrmqRMiPL2TfLjfV43oCUSiWwAGh4e1rZt25YFrJaWlmxn8GoJ46o3AACqwq8JJ4j8prigGhoaCl4JFha/kDQ0NKS9e/eqv7+/4PJC2x8/flxDQ0OanZ1VS0tLNhS5P1OpVMlw6fe8O3fuzP7frcWamJjIef+SyaQOHz5cdP9rjaAEAEAN8KvJKVTLMz4+ru3bt+eElaDbe0Og28ToDTnFlPO8iURiWXOiG9CqiaAEAIici5cWdXFhqeR6rY31am1qqEKJoi+ZTGbDRX748NbKuM1lblhxhx8Isn1+R2q3Gc6tRSpW41Psed2+St59p9PpZQEsnU5nx1mqFoISACBy3nj/nF59+6OS673w1IN6qeehKpQoPH5NZ6lUSuPj4zkDQu7bt08TExPZMDI+Pp5TezM5OanJyUn19vZma2u865Tavq+vT6Ojo9ngNDo6mjMIpLdztlep5x0eHtbhw4ezQWl8fFzd3d3LOm0Hadpba8ZmxiK6Yzk3xp2fn59XU1NT2MUBAKyB/Bqlsxcv68XDH+pbex/RA623O17Xco2SG4QOHz6syclJDQ4OateuXdmO1OPj4xoaGlrWuXlkZCRbI3TixIlsU1k6nda2bdsKXsHnzQJ+20uZWqHJyUklEglNT09rYGBgWe1TT09PzpAMQZ93bGwsu47fiN99fX16/fXX1yQsXbp0Sc3NzZLUbK295LceQQkAEHlTv5jXv/n2u/rr//ikOr/YHHZxUMTIyIi6uroqMnZUX1+fjh4tdUvYYIIGpTBuigsAAGqU9zYla2lkZGTZ/eOqgaAEAADW1N69ezU+vmzoxBVLp9P67LPPqjLCeT6CEgAAWFNuPyq/0cLLNTY2FtrYVFz1BgAA1pwbltaC98q+aqNGCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwAdBCQAAwMf6clY2xiQkPSOpz1rbE3CbQUlpZzJhrR3xLOuWNCDpuKSUpB5JJ6y14+WUCwAAoBICByVjTJeknZISkloCbjMoSdbaMWe62xgzaq0dcFZJSOqW1KtMUBomJAEAgKgIHJSstZOSJo0xvWXsf5+kbZ59TBhjjitTi+TaZq1Nl7FPAACAqqhYHyVjTFKZprZ0gWXdlXpeAACAtVJWH6UyJX3mp5VpcnM9Y4yZVaY5b7u1dqjYTo0x9ZLqPbMaV1FGAAAAX5UMSn7cUCRJk5JkrU1JkjGm3xhz1FrbV2T7fZL2V7aIAAAA4QwPkO0Ibq1NuSHJcURSr3N1nZ+Dkpo9jy2VKCQAAEAlg1LKZ37CXZbfMdzTn8mv2U7W2iVr7SX3IWlh9UUFAABYrmJByakpSjuduvOXTTi1Rke9yz01SX4hCwAAoGpWEpQKjqFkjEm64yZ5HFRmnCR3nV5JY1K29mgkr+mtX9I4wwUAAIAoKGfAyaQyA0PuldRljBlW7ija7ijb2ZG3rbUjxphBTxPbLs9gk5J0MC9c3VuiIzcAAEDVlDPgZEqZEDTis3xMTm1R3nzv+uN5y9J++wMAAAgbN8UFAADwQVACAADwQVACAADwQVACAADwQVACAADwQVACAADwQVACAADwEXgcJQAAED8XLy3q4sJSyfVaG+vV2tRQhRJVF0EJAAD4euP9c3r17Y9KrvfCUw/qpZ6HqlCi6iIoAQAAX88+vlU9Hfdlp89evKwXD3+ob+19RA+0bs7Ob22sD6N4FUdQAgAAvlqbGgo2qT3QulmdX2wOoUTVRWduAAAAHwQlAECk3bxldfp8WpJ0+nxaN2/ZcAuEWCEoAQAi69jUjJ4cfkcvvzklSXr5zSk9OfyOjk3NhFwyxAVBCQAQScemZvT8oUnNzC/mzL8wv6jnD00SllAVdOYGgAiJ+5g1rpu3rA68dUaFGtmsJCPpwFtn1NPRpnV1psqlQ5wQlAAgQuI+Zo3rg49nl9UkeVlJM/OL+uDjWT2x/d7qFQyxQ1ACgAiJ+5g1rosL/iFpJesBK0VQAoAIifuYNa7WxmDNikHXA1aKztwAgMh5bFuL2psb5Nf7yEhqb27QY9taqlksxBBBCQAQOevqjPbv6ZCkZWHJnd6/p4OO3Kg4ghIAIJJ2d7brtee61Nac27zW1tyg157r0u7O9pBKhjihjxIAILJ2d7arp6NNh0+c08tvTumVpzu1d9dWapJQNdQoAQAibV2d0Y4tCUnSji0JQhKqiqAEAADgg6AEAADgo6w+SsaYhKRnJPVZa3sCbjMoKe1MJqy1I0XWPR50vwAAAJUWOCgZY7ok7ZSUkBRo4AonJMlaO+ZMdxtjRq21AwXW7ZXUHbQ8AAAAlRa46c1aO+kEnlQZ+98nacyzjwlJ/fkrOTVVjBoGAAAipWJ9lIwxSWWa2tIFluXXHD0j6UjA/dYbY5rch6TGVRcWAACggEp25k76zE8r03wnKRuaJsrY7z5J857H+ZUVDwAAoLgwrnqbVW4zW8JaW05z3kFJzZ7HljUsGwAAQFYYI3NnQ5Ixpt/t6B2UtXZJ0pJnH2tYNAAAgNsqWaPkV0uUkJRyrqI7WcHnBwAAWJWK1ShZa1PGmLQxJpnftGatnXD6JnV5OnZvl7JDCqSsteOVKhsAAEAQKwlKBS/jd65y680bUPKgMmMjueMo9br/d4YKmPBs3yWpv9iAlAAAANUUuOnNGJN0ansGlKkJGnaCj6vbWZblhJ6EMabXWXdXkcEm9zn/Hy4wfAAAAEDVBa5RcprPRpxHoeVj8gwu6ZnvXb9gc5rTzEZTGwAAiBRuigsAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOAjjHu9heripUVdXFgquV5rY71amxqqUCIAABBVsQtKb7x/Tq++/VHJ9V546kG91PNQFUoUDgIjAAClxS4oPfv4VvV03JedPnvxsl48/KG+tfcRPdC6OTu/tbE+jOJVDYERAIDSYheUWpsaCtaQPNC6WZ1fbA6hROEgMAIAUFrsghIyCIwAAJTGVW8AAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+uOoNAJCDAWmB2whKIeAkBCDKGJAWuI2gFAJOQgCijAFpgdsISiHgJAQgyhiQFriNoBQCTkIAANwZuOoNAADAB0EJAADAB01vACKBq0EBRBFBCUAkcDUogCgiKAGIBK4GBRBFZQUlY0xC0jOS+qy1PQG3GZSUdiYT1tqRAvuTpO2SkpK+aa1NCzWPphZ4cTUosBznyfAFDkrGmC5JOyUlJLUE3GZQkqy1Y850tzFm1Fo74KwyLGnYWptylo9KOiopUAjDnY2mlgxOhAD8cJ4MX+CgZK2dlDRpjOktY//7JG3z7GPCGHNckhuUkpJ6Jbm1TNO6XcOEGkdTS0YUToSENSCaOE+Gr2J9lIwxSWWa2tIFlnVbaycKNN/tkjRRYr/1kryfiMbVlhXhoKklIwonwiiENcArP7yfvXg556er1sM758nwVbIzd9JnflqZ5rscTk1VQlJfif3uk7R/FeUCIiUKJ8IohDXAyy+8v3j4w5xpwjsqLYyr3mbl6ePk6dCdkHQ0QEfug5L+xDPdKOn8mpYQiJkohDXAKz+8+yG8o9LCCEo5HcGdYOR29u43xsxJ2uYXmKy1S5Ky9bHGmIoVFADCdPOW1enzaUnS6fNpPdzepHV18Tjn+YV3oNoqeQuTlM/8hKSUMSZhjBl2apRcE87y7gqWCwAi79jUjJ4cfkcvvzklSXr5zSk9OfyOjk3NhFwyIF4qFpScS/7TTqfu/GUTyvRhGlRuDVPC+ZmuVLkAIOqOTc3o+UOTmplfzJl/YX5Rzx+aJCwBVbSSoFRwDCVjTNIdN8njoDy1Q06H7TEpO9zAiDuGkmOvpEknSAFA7Ny8ZXXgrTOyBZa58w68dUY3bxVaA8BaCxyUPEFoQFKX02zmHVOpW7fHR5IkOaNwJ4wxvc66uzyDTUrSQWPMoPtQpkbpqZW+GAC4033w8eyymiQvK2lmflEffDxbvUIBjvx+c3EI7OUMOJlSZmDIEZ/lY3Jqi/Lme9cfz1uW9tsfAMTRxQX/kLSS9YC1cmxqRgfeOpMN8i+/OaVvv3NW+/d0aHdne8ilqxxuigsAjiiMUN7aGGy/QdcD1oLbby6//sjtN/fac101G5YISgDgiMII5Y9ta1F7c4MuzC8W7KdkJLU1N+ixbYFuuQmsWql+c0aZfnM9HW01OXwFQQkAHFEYoXxdndH+PR16/tCkjJTz5eR+Be3f01GTX0iIpnL6zT2x/d7qFaxKCEoA4IjKCOW7O9v12nNdOf1BpExNUq33B0H0xL3fHEEJACJod2e7ejradPjEOb385pReebpTe3dtpSYJVRf3fnOVHJkbALAK6+qMdmxJSJJ2bEkQkhAKt9+c36fPSGqv4X5zBCUAAODL7TcnaVlYikO/OYISAAAoyu0319ac27zW1txQ00MDSPRRAgAAAcS13xw1SgAAIJA49psjKAEAAPiIdVCK4839CuE4AABQWGyD0rGpGT05/I5efnNKUubmfk8Ov6NjUzMhl6y6OA4AAPiLZVByb+6XPyS7e3O/uIQEjgMAAMXFLiiVurmflLm5X7Wan8Jq9oracQAAIIpqZniAhYUFGVO69/2Jf0oHurnfj86c165fS6xdAQuY+Omn+m9/c1b/vHBdUqbZ69WJn+m//OsH1P3lL1T0uaN0HFxXrlzJ/lxYiF2Gz4rCcaAMlCFqZUBGFN6LKJRhLSwsLARar2aC0g9/+EPdfffdJdf7cG69pI0l1/ubv3tfF6c+X4OSFTaVXq9DP3cH7rod8P55YUkv/eX/1XNfWlRnonLPH5Xj4PWLq3WSNundd9/Vx3ffqspzRlEUjgNloAxRKwMyovBeRKEMa+Hq1auB1quZoHTXXXepqamp5HptAZuS2po3qqmpMuND3LJWf33Gb2nmOb8/s1G/cb9UF6CWbCWicBzypU2mTJs2bVJTY+2PzeEnCseBMlCGqJUBGVF4L6JQhrXw+efBKgFqJig1NDQEqlHasdHqnn+8prkl/6DQ0mC0o31jxULKTz67qfR1/2YvSZpbkj5ZbNDD966rSBmicBzyNdy4KWnReS8r87rvBFE4DpSBMkStDMiIwnsRhTKshaA1Sndu4+IK1RmjZx++q+g6//7Ld1U0HMwXCScrWW8lonAcAACIutgFJUna2bZef/hIve6pzw0BLQ1Gf/hIvXa2Vbairbk+WPgIut5KhX0cAACIuth+E+5sW6+u+9bpR5/c0HfP3NB/6Nigr9y/oSo1KL/eUqd76k3JZq9fb6l8jg3zOAAAEHWxrFFy1Rmjbc2Z9tVtzeuqFg6i1uwV1nEAACDqYh2UwkSzFwAA0ce3cYho9gIAINrKCkrGmISkZyT1WWt7Am4zKCntTCastSMFlkvSdkmy1g6UUybXX/zFX5QcHqC1tVVPP/10zry3f/A9zf3DL3XkZ+uLdp7euXOnHn300ez09evX9Z3vfCdQ2X73d39XbW1t2enp6WlNTExkp+eXrOZ++blO/Gy9fuYpw4YNG/T7v//7Ofv60Y9+pJ/+9KclnzOZTKqnJ/ctOnToUHZE1XxuGY78bL329HxFDz/8cHbZ7Oysjh49WvI5JenZZ5/V5s2bs9OnT5/We++9V3K7e+65R4997d/lzPv+97+v8+fPl9x2x44deuKJJ3LmjY6OBirv7/zO7+j+++/PTn/yySf6wQ9+EGjbgYHcj+p7772n06dPl9xuy5Yt+vrXv54z78iRI5qbm5OU+17kfyafeOIJ7dixIzt9+fJlvfHGG4HK29fXp5aWluz0T37yE/34xz8uuK63DL/aslnPPfdczvLjx48rlUqVfM4vf/nL+spXvpIz78///M9148aNktt2/MZXJf1qdvrChQv63ve+V3I7SfrGN76hu+663bx96tQpnTx5suR2hc4RC5N/rSM/myt5cUWlzhHFPg/VOEdYGc3U/YouL9yl73xyXe23/p++8tu/Fco54plnnsmZF8dzRLHPg1Sdc0R+GTZt2hTKOaK7u1vbt2/PTpd7jggqcFAyxnRJ2ikpIaml+NrZbQYlyVo75kx3G2NG3TBkjBm21g551h81xhwPGsK8rl69qlu3io8QWmhAysXFa7q1eEVXr6zTuhv+J8KlpaVl8y5fvhyobPnl+vzzz3O2vXrd6tbizWVl8J7ob5d3MdDzXrt2bdm8K1eu+G7rLUP+IFy3bt0K/Fqtze2gfuPGjUDb1tfXL5t37dq1QNuu5r3Jf6357005lpaWVvzeXL16Nbut3+dB0rKTh7V2zT6HOeXxlOFKgZNx0PdmcXH5eGFXrlzR9evXS257cxWfw3xB35tC54hb1xd19crloucH9znyrcV7U+zzUOlzxHzjNs20/6ZubNiszZL+SdIvb1xW68JlPezZjnNEMGtxjij2eZCqc44oVQb3NVT6HLGa76pyBA5K1tpJSZPGmN4y9r9P0jbPPiaMMcclDTi1U13GmIS1Nu2sMirplDEmaa0tHUU97r777pI1Sg0NDQXmbVRdwybdvWm9Nhf5i7HQL6n3r6Ji6upyu4KtX78+Z9ubS1Z1DZ8vK8OGDRsKvoYgz7tx4/Lbk2zatMl3fW8Z1q/P/VjU1dUFfq3599vbsGFDoG0LvXcbN24MtO1q3pv815r/3pSjvr5+xe/N3XffnT2Z+30epOWfCWPMmn0OvbxlKPS5CfreFPqd27RpU8Ev+HzrVvE5zBf0vSlU3rq7GnT3ps1Fzw/uc+Rbi/emnM+DtHbniM823q9zv/Lby9a5sX6T/mp2k7Ze+Dzbl5JzRDBrcY4o9nmQqnOOyC9DWOeI1XxXlcPkp/uSG2SC0j5r7aMl1ktKmrbWmrz5VlKPpJOSPpb0lBPC3Ka9OUmPuvMClKdJ0vzRo0f1hS+UfyPZn8/f1B+/t6g/fqJBX2oOZ4RRyhCdMkRBFI4DZYh3GW5Zqz/629Ij9//3r1Rv5H5kxPUzWQmffvqp+vr6JKnZWnvJb71KXvWW9JmfVqavUtpae09eIOp2fvrWJhlj6o0xTe5DUuPaFBdAVNyyVh/P35QkfTx/U7fK/IMOq/OPs7eKhiRJml20+sfZO/eGqEBQYVz1Niv/Pk77JA14muL81tm/1oUCEA0nL3yuN35yPftF/d0zN/RX05/r2YfvYtiMKonCbZaAqAhjHKWCIckYMyzpsNvxu4iDkpo9jy1rWzwAYTl54XP96YdLy2oz5pas/vTDJZ28EOxu31idqNxmCYiCSgYlv+azRP4yp9/TdP7QAYVYa5estZfch6SFVZcUQOhuWas3flL8Spf/9dPrNMNVgXubpWKqdZslIGwV+5Q7V62lnU7d+cuygwgZY7qdee4QAolC2wCojrD6B9EvJjqidpslIEwrCUp+TWdJz+CRroO63UHbrTka80x3SepSZtiBpBOQ+pXpxwTETtidmE9e+Fx/9LfX9N0zmbFYvnvmhv7ob69VpcmLfjHRwm2WgIxyBpxMSuqVtFeZ8Y+GJZ2w1o47q3RLGpCUbT6z1o4YYwY9Yy/t8gw2mZD0tjJNccPe5wrSBJdvcXFRV69eLXczLS7a29tvCOevI8oQnTKE6e8/tTp6Vko7rU/fPXND3zt7Q30PSP/qC5U/Hn//qdXrZ5bPd/sHfbNjqaLlaFCwANSgJV29WnowurUQhc9kmGXoaJL+62NW/2dG+t9npd97QPrNdqs6c71q7wFyxf0zuZYKDXZZSDkDTqaUCUEFQ4zTdLasI3Ze6Bn3zE9Luifo85dy/fp1XbrkOwyCrytX6yRt0pUrV3TJhlOlTxmiU4awTKXX69DP3YHXbp940tcz4eW5L11TZ6JytTq3rHTko03Ocxc68Vkd+chq24YrqqvQebHVSM0bNmn+hn8ZmjdYtZorWsGv+opE4TMZhTJ8oS5Thi/UXdHlhXj9bkZNFD4PYZTh0g2jhRKj40tS4warpg3B/ugKMvq3VEM3xf3a175W8PYDpZyZWdC3f/b3evLJJ9XRHs6QTJQhOmUIw81bVv/jT9+XVOiX1shImphN6D//3mNaV6GUcuKf0pr/h2L3oDKav2HU9i9/U7t+LVGRMkjSxgc+1X/6y0y1lvdUZ5x///jf/gt1f7n8gWVXKgqfScoAryi8F2GU4X/++Od67e/OlVzv+d/aqt/77S8F2mfQypWaCUqNjY1qbCz/Ddt0KZOGN23atKLt1wJlyISF6blM17TpuRva+cDmioWCqHlv+jP984L/XzZW0oVLS/rJpzf0xPZ7K1KGhc+DnTAWPq+r6Ofj6V2N2rixQQfeOqOZ+dvV4m3NDdq/p0O7O9sr9tyFhP17QRnCdfHSoi4uLL9PXL7Wxnq1Ni2/FUclROG9CKMM3/itB/X1R7Zmp89evKwXD3+ob+19RA+03r5tSWtjvRobg70XQe9MUjNBCXeuY1MzOV+ML785pW+/czaUL8YwXFwI1k4edL2VaA14Ygm63mrs7mxXT0ebDp84p5ffnNIrT3dq766tsQnOiI433j+nV9/+qOR6Lzz1oF7qeagKJYqv1qaGgmH0gdbN6vxic0Wfm6CEUB2bmtHzhyaXdeO9ML+o5w9N6rXnumo+LEUhpDy2rUXtzQ26ML9YsEu1UaZW57FtfoPqr611dUY7tiQkSTu2JAhJCMWzj29VT8d92elitRioXQQlhObmLasDb50p+MVslflyPvDWGfV0tNX0F2UUQsq6OqP9ezr0/KFJGRXqHyTt39NR0+8DkC/MWgxEB8OqIjQffDyb0w8ln5U0M7+oDz6u7WG13JAiLb/Wq5ohZXdnu157rkttzblfDG3NDbGo2QOAQghKCE0U+uZERVRCyu7Odr079FW98nSnJOmVpzv17tBXCUkAYoumN4QmCn1zoiQqnZjpHwQAt1GjhNC4fXP8voaNpPYqdiCOAkIKAEQLQQmhiUrfHAAA/BCUEKqo9M0BAKAQ+ighdFHpmwMAQD5qlBAJ9M0BAEQRQQkAAMAHQQkAgDvAzVtWp8+nJUmnz6d181awm7pidQhKAFAAX0qIkmNTM3py+B29/OaUpMzNw58cfkfHpmZCLlntIygBQB6+lBAl7s3D82/55N48nM9lZRGUAMCDLyVESambh0uZm4dT41k5BCUAcPClhKjh5uHhYxwlAHCU86X0xPZ7K1KGi5cWdXFhKTt99uLlnJ+u1sZ6tTbF4z6IccbNw8NHUAoBJ0IgmqLwpfTG++f06tsfLZv/4uEPc6ZfeOpBvdTzUMXKgWjg5uHhIyiFgBMhEE1R+FJ69vGt6um4L0AZ6itWBkSHe/PwC/OLBZuEjTK3fIrTzcOrjaAUAk6EQDRF4UuptamBmmRkuTcPf/7QpIyU87nk5uHVQVAKASdCIJr4UkIUuTcPP/DWmZw+dG3NDdq/p4Obh1dY7IIS/YMyOA5AYXwpIYq4eXh4YheU6B+UwXEA/PGlhCji5uHhKCsoGWMSkp6R1Get7Qm4zaCktDOZsNaOrHafq0H/oAyOA1AcX0oApDKCkjGmS9JOSQlJgXoyOiFJ1toxZ7rbGDNqrR1Y6T5Xi/5BGRwHAABKCzwyt7V20gk8qTL2v0/SmGcfE5L6V7lPAACAqqjYLUyMMUllmtrSBZZ1r2K/9caYJvchqXEVxQQAAPBVyXu9JX3mp5VpalupfZLmPY/zq9gXAACArzBuijur1fVHOiip2fPYshaFAgAAyBfG8ACr6rRtrV2SlB0AyBiuRMHK5Y8n5YfxpAAgnioZlPw6aCeKLAOqym88qXyMJwUA8VSxoGStTRlj0saYpLU2lbdsolLPC5Qjfzypsxcv68XDH+pbex/RA62bs/MZTwoA4mklQalg05lzlVtv3oCSByV1yxkiwBjTK89wAaX2CVSa33hSD7RuVucXm0MoEQAgSsoZcDIpqVfSXkldxphhSSestePOKt2SBiRlg5K1dsQYM+gEJEna5Q42GXCfAAAgRHG/N2jgoOQ0n43IE4Tylo+pQG1RXg3TeN6yovsEAADhivu9QWN3U1wAABBc3O8NSlACAAC+4n5v0DAGnAQAALgjEJQAAAB8EJQAAAB80EcJAJAj7peDA14EJQCRwJdzdMT9cnDAi6AEIBL4co6OuF8ODngRlABEAl/O0RH3y8ELuXnL6vT5tCTp9Pm0Hm5v0ro6E26hUBUEJQCRaPbiyxlRdWxqRgfeOqOZ+UVJ0stvTunb75zV/j0d2t3ZHnLpUGkEJQA0ewE+jk3N6PlDk7J58y/ML+r5Q5N67bkuwlKVVbt2j6AEgGYvoICbt6wOvHVmWUiSJCvJSDrw1hn1dLTRDFclYdTuEZQA0OwFFPDBx7PZL+RCrKSZ+UV98PGsnth+b/UKFlNh1e4RlICQRaF/EIDlLi74h6SVrIeVC7N2j6AEhIz+QUA0tTYG+8Mk6HpYuTBr9whKQMjoHwRE02PbWtTe3KAL84sFazKMpLbmBj22raXaRYudMGv3CEpAyOgfBETTujqj/Xs69PyhSRkpJyy5jTv793TQkbsKwqzd46a4AAD42N3Zrtee61Jbc+4XcFtzA0MDVJFbu+cXSY2k9grV7hGUAAAoYndnu94d+qpeebpTkvTK0516d+irhKQqcmv3JC0LS5Wu3SMoAQBQwro6ox1bEpKkHVsSNLeFIKzaPfooAQCAO8Luznb1dLTp8IlzevnNKb3ydKf27trKyNyoTYwfBAAoV7Vr9whKCA3jBwEAoo6ghNAwfhAAIOoISggN4wcBAKKurKBkjElIekZSn7W2J+A2g5LSzmTCWjtSznIAAICwBB4ewBjTpUxISkgKNKKTE4JkrR2z1o5JmjTGjAZdDgAAEKbAQclaO+mEmVQZ+98nacyzjwlJ/WUsBwAACE3FBpw0xiSVaUpLF1jWXWp5pcoFAAAQVCU7cyd95qeVab4rtbwgY0y9JO9lUI1llwwAACCAMG5hMqvifZxKLd8nad7zOL92RQMAALgtjKBUqiN4qeUHJTV7HlvWolAAAAD5Ktn05tfpO+EsK7W8IGvtkqTsfS+M4caEAACgMipWo2StTUlKO52285dNlFpeqXIBAAAEtZKgVLBpzBiTdMdF8jgoqduzTq88wwEEWA4AABCacgacdIPQgKQuY8ywE2xc3c6yLGeU7YQxptdZd5e1diDocgAAgDAF7qPkNJWNOI9Cy8dUoDYo75Yk4+UuB6rl5i2r0+fTkqTT59N6uL1J6+roAwcAccZNcQFJx6ZmdOCtM5qZX5QkvfzmlL79zlnt39Oh3Z3tIZcO1XLx0qIuLmSvFdHZi5dzfrpaG+u5oTMQEwQlxN6xqRk9f2hSNm/+hflFPX9oUq8910VYiok33j+nV9/+aNn8Fw9/mDP9wlMP6qWeh6pUKgBhIigh1m7esjrw1pllIUmSrCQj6cBbZ9TT0UYzXAw8+/hW9XTcV3K91sb6kusAqA0EJcTaBx/PZpvbCrGSZuYX9cHHs3pi+73VKxhC0drUQJMagBxhjMwNRMbFBf+QtJL1AAC1haCEWGttDFZ7EHQ9AEBtISgh1h7b1qL25gb59T4yktqbG/TYtlK3IAQA1CKCEmJtXZ3R/j0dkrQsLLnT+/d00JEbAGKKoITY293Zrtee61Jbc27zWltzA0MDAEDMcdUboExY6ulo0+ET5/Tym1N65elO7d21lZokAKFiENTwEZQAx7o6ox1bEpKkHVsShCQAoWMQ1PARlAAAiCgGQQ0fQQkAgIhiENTw0ZkbAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB0EJAADAB/d6AwCggIuXFnVxYSk7ffbi5ZyfrtbGeu7HVsMISgAAFPDG++f06tsfLZv/4uEPc6ZfeOpBvdTzUJVKhWorOygZYwYlpZ3JhLV2JOA2rnuttUPlLAcAoNqefXyrejruK7lea2N9FUqDsJQVlNxAY60dc6a7jTGj1tqBItsclXTcs02/MWbYDUOllgMAEIbWpgaa1FB2Z+59ksbcCWvthKR+v5WNMUlJvZKOeGYfkTRojEmUWl5m2QAAANZU4KDkhJqEtTZdYFm3z2ZJSfJu4/n/zgDLC5Wj3hjT5D4kNQZ9DQAAAOUop0Yp6TM/LSnhsywlST61Q8kAywvZJ2ne8zjvsx4AAMCqrMVVb7OSWgotsNamjDETkroljUu5tU+llvs4KOlPPNONIiwBAFCzwhyqYS2CUsGQ5LLW9hhjho0xLcqEqpSzKBVkeYH9LUnKHi1jzCqLDwAAoizMoRrKCUoFg4syzW5+yyRJ3ivYPM1sJ4MuBwAA8RXmUA2Bg5LTTJY2xiSttam8ZRN+2xljuqy1k55Z3ZLG3U7bpZYDAIB4C3OohnKHBzioTJCRJBljeuUZLsAYk8wbPFKSjub1OxqQNFTGcgAAgFCU1UfJWjtijBl0ApIk7cobbLJbmaDjHa17QFKXM7zAdkkDeTVSpZYDAACEouzO3Hm3LBnPWzYmTw2TM29Ckm/TXKnlAAAAYSm36Q0AACA2CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+CEoAAAA+yr7XG1BLLl5a1MWFpez02YuXc366Whvr1drUUNWyAQDCR1BCrL3x/jm9+vZHy+a/ePjDnOkXnnpQL/U8VKVSAQCigqCEWHv28a3q6biv5HqtjfVVKA0AIGoISoi11qYGmtQAAL7ozA0AAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCDoAQAAOCj7Hu9GWMGJaWdyYS1diTgNq57rbVDecv7JSWc/W6XdNBamxYAAECIygpKbuCx1o45093GmFFr7UCRbY5KOu7Zpt8YM+yGJWefY24wMsYkJL0uqa/8lwMAALB2jLU2+MrGzEna5q3tMcZYa63xWT8paVrSPXlBaM6dZ4w5bq3tydtu2bwiZWqSND8/P6+mpqbArwUAAMTXpUuX1NzcLEnN1tpLfusF7qPkhJ5EoSYxY0y3z2ZJSfJu4/n/Tudn2hhz3AlQ7vOkipSj3hjT5D4kNQZ9DQAAAOUopzN30md+Wpn+RYWkpGwtkt/+vun8f84YMyypu1hTnqR9kuY9j/PFCg0AALBSa3HV26yklkILrLUpSROSsjVO+bVPTg3TsKRxSYOS+nyCleugpGbPY8vKiw4AAOBvLYJSwZDkcvoa7XI6cfcqE6yk27VNw5JS1to+Za54a5F0qsj+lqy1l9yHpIU1eA0AAADLlHPVm1+/oUSRZZIk73AAntqik55+TxPOeilJjxpjThljeq2142WUDwAAYE0FrlFyQkzaCTf5yyb8tjPGdOXN6pY07jS5JXV7TCav0aDlAgAAqJRym94OKre/Ua+kMc90Mm9wSUk6mtcvaUDSkJQNWF0F+iQ9Sm0SAAAIW1njKEnZASLdprZdec1q/ZKGrLXbPfO6JXXp9qjbo07tlLs8ocyVbJ/p9hV0Y4WGIfApD+MoAQCAsgQdR6nsoBQ1BCUAAFCuNR9wEgAAIG7KviluVF265BsGAQAAcgTNDbXQ9PZFMTo3AABYmS3W2l/4LayFoGQk/apWPvBkozJBa8sq9lELOA4ZHIcMjkMGxyGD45DBcciopePQKOmXtkgYuuOb3pwX55sES8nkLEnSQrHOXLWO45DBccjgOGRwHDI4Dhkch4waOw4ly09nbgAAAB8EJQAAAB8EJWlJ0gHnZ5xxHDI4DhkchwyOQwbHIYPjkBGr43DHd+YGAACoFGqUAAAAfBCUAAAAfBCUAAAAfBCUAAAAfNzxA04GZYxJSHpGUp+1tqfA8kFJaWcyYa0dqV7pqifgcZCk7ZJkrR2oXumqp9RxyFv3eKl17lRBjoMxZljStDM5a60dr1LxqibA70W/pIQy54jtkg5aa9PVK2H1lDoHxOhcGeQ4+C6vFeW8zlo9V8YiKBljuiTtVOZE11Jg+aAkWWvHnOluY8xorX3wAxyHYWvtkGd6tBY/+KWOQ966vZK6q1CsqgvweUhIelvSU9batLP+KUkmf907WcDzw5gbjJzj8rqkvqoVskpKnQNidK4sdRzicq4M/Dpr+VwZi6Y3a+2k84ud8llln6Qxz/oTkvqrUbZqKnYcnJN/l/PTNSqp2xiTrE4JqyPA50FS9pgUDVJ3sgDHYVjSYTcgWGsnJdXUF4EU6Dj0eGuPnP8nKl+y6gp4Dqj5c2Wp4xCXc2U5r7PWz5WxCErFOG94olA1ujGmJtNxETsleX8B3C+ORPWLEgnPSDoSdiFC1C9p3Ply6JayX4xxkzbGHHe/MJxzRtGQfQfzPQfE7FxZ6lwYl3Nl0NdZ0+fKWDS9leD3F0Batfeh9+Wc/O7Jm+2e/Gr1S8GXc+KPYyiQlA0DktSlzPufMsaMSjoaw7D0TWWaHOeMMSOSpmutqUkKdA7Y6bNpWjV0rix1HOJyrgz6OuNwrox9jVIRs6rhqsSA9kkaqNVOqyUkrLU1c9JbATcopZ2mqZSkIUlHQyxTKJzP/7CkcUmDkvrymiNqWZBzQBzOlaWOQ1zOlYVeZ82fKwlK/mr9F78o50qnw26nzTgxxvTX4pVdK3TS/Y/bN6cGm1mKcn4XUtbaPmWu/GlRpoapppVxDqjpc2Wp4xCXc2Wh1xmXcyVByb+qNFFkWU1zrl6YrtXLfotxroA6WXLF2uf32U/Lv7m65nj65UxIkrU2Za19VJl+S73hlq5yfM4BsTtXljoXxuVcWeh1xulcGfs+StbalDEmbYxJ5lcfxrAvhjyddt3LfxOSWmq9atWjRZkrPdxak+1S9rLoVBz+epKyvxcpZULRpGdRQjE5OTqSuj1mkNdolctRNcXOAXE6V5Y6F8blXOn3OhWjc2XcgpJfFfFBZTqpuR+EXnkuga1BBY+D8xdCl5wrnZzZtXwslh0H54SfPek7x6S/xv9i9Pu9GJK0V05Qcn4vJpxhAmpRwc+DMWbIGJN/tdejtdihO8A5IBbnylLHIS7nymKv0wmEsThXGmtt2GWoOOcN7lXmpN8laUTSCW/idVOwM7nLO8hWrSh2HJy/Ej5WgatXrLW1NsBgyc+Ds567Tq+zzvFa+ss54O+FOyK1JN0bt98LZ3lCmU6sn+n2FV5jtsY67gY9B9T6ubLUcYjLubKc11nz58o4BCUAAICVoDM3AACAD4ISAACAD4ISAACAD4ISAACAD4ISAACAD4ISAACAD4ISAACAD4ISgFgyxiScQfUAwBdBCUBc7VOMbvALYGUISgDiqquG71sHYI0QlADEjnPH8+NhlwNA9BGUAMRRn6TxkmsBiD2CEoA4SlprU2EXAkD0rQ+7AADgxxjTJWmnpO2STkiakNTvLE5ba8dWsM9eSUeLLNslaVpSynnMWmvTZRceQE2gRglAJDmX7ndba8estUOSXpe0z1o74qwytMJd75V0pMDz9UvqsdYOOQEsoUxg2rnC5wFQA6hRAhBV/Z5Q5Jp2fk5KGljhfhP5NUTGmKSkYUnbPLPTkmStnVjh8wCoAQQlAFGV7WztBJmEnJqg/PDiLO9Vpqlsl6TRQn2QnFqj0QLPNSppIi9A9SgTyADEGEEJQCTlBZ1uSakifYWOWmsflSRjzISktyU9WmC9PmttT4H53cpcCefVpUyfKAAxRh8lAHeCHuVdzu/efsTp8J3lhKmEU8uUv346f8ee9fJrjxhrCQBBCUA0Oc1krl5lrnrLLvPULvl1tu7Km/ZrdpOUW4PlDEgpa+2EMaYrP4wBiA+CEoDIcULSsPP/XnmawArcyDYhaTZvXlpSS968nkIds52AlHLDkLP/AWX6O0mZK+/oqwTEFH2UAETRhKQxJzCdVCa4DBljJKklb/yktJaHooQ84clpXis2wGSfpAFjzClJstb2GWOOOs9PSAJizFhrwy4DAKyYUxP0utuZ25k3J+lRtznNGDMs6TA1QwDKRdMbgDuaE34S7rTTdJbKu2qui5AEYCVoegNQC/qcWqMTyoyjlL3U36lxIiQBWBGa3gDUNGPMqKRhboILYCVoegNQ61oISQBWiholAAAAH9QoAQAA+CAoAQAA+CAoAQAA+CAoAQAA+CAoAQAA+CAoAQAA+CAoAQAA+CAoAQAA+Pj/rm6oML80iigAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEfCAYAAABWPiGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAActUlEQVR4nO3dT4wc53nn8d8jiV4xG2UrQzNWxIXCNGEx2IMQ9jQPsSIjNrvBGAZ2mc2IPliHXDI8GoLAf/DZMMeEoBjIhXPw+iADDjnhcg/EmuiWENiwL5oZBlL2QBls+N9oZTAclmXY1IqSnz10VbN6pqt6ij01XdX8fgACXW+9XfN0s6uffv/UW+buAgBgqx6ZdAAAgGohcQAAciFxAAByIXEAAHIhcQAAciFxAAByeWzSASR98pOf9P379086DAB46K2srPy7u+8dtq9UiWP//v1aXl6edBgA8NAzs5+m7ctMHGYWSJqXFLr7YkqdOUmhpEBS191XE/vmJXUlyd07OeMGAJTQqBZHI2tnlFha7n4i2m5LakWPL0n6e3cPo8ckDgCYApmJw907ZlbLqHJc0s3Edmhm9cTzQzOrufsLY8YJACiJcWdVBep1U8XWJc0oaqlELZLAzBbSDmBm82a2bGbLt27dGjMcAEDRipiOG0T/5O5hNOZRT2u5uPuiuzfcvbF379ABfABAiYw7qypUlCQiM4oGwyXt2VCvltgHAKiocVscFyUdSGwHUQujo8GEUpPEPFsAmAKjpuM21ZslFZhZN55Sa2Yrko5Eg9/tqJ4kLUj9QfF2NB1Xkr7u7mExLwEAsJNGzqrSkGm07j6beLyU8tyh5QCAamOtqimxtPRgeToMQ3U6XGIDYOtKteRI1YRhqMXFRQVBoPn5+aF1lpaWFASBwjBUrVZTvV7vl6+vr2tlZUUvvPCCms1mZv3FxUXNzMxofX1djUajXx7vS/79tGPHdWu13gS3ZrOpIAgkSd1ut18OAJncvTT/ZmdnvUra7bYvLCz4hQsXhu6/c+eOz8/P97ebzaa7u6+srPjKykq/PAiCkfWHlcfPSf79tGO7u8/NzfmdO3f6j5NOnTo14tUCeJhIWvaU7+rKtjhOnDiRuu/FF1/U888/L0n6wQ9+oNdee03f//739dnPfnZT3QsXLvQff+1rX9PPfvazTeVpms2mut30GcYXL17UgQP3J50FQaDV1VWtr6/r0qVL/b8xMzOj1dVVLS8vD63f6XQ2lccthMXFRc3NzfX3pR1743MvXbo0EGsYhiNfLwBID9EYx7CkUbQwDPtdQZL6XU3NZnMgMa2vr6ter6fWr9Vqunnz/sou3W63n7DefPPNgS6mtGPHqw6HYagwDHX69OmBWOMkBQCjVLbFsZUWgSQ9//zz/dbHKF/96lfHCWlLNv6yP3369KZf/xvrz83Nqd1uq9vt5moZJI8dPy8IAtXrdZ0+fXpgXOPw4cPqdrsDYycAMMxD0+KYhHiQOxa3HmJLS0v60pe+1B+8zqofJ8r4i73RyFy4eNOxa7XawN+Ou6w2xgsAo5A4CnT8+PGBLqYwDPtf/Kurq6rX66rX6/2up7T63W5Xp0+fVq1WU7fbVaPR6H/JD2uBDDt2s9kcqBsfJ7k9MzOzvW8AgKlU2a6qMuh0Omq32/2ps/Gv+9nZWb3++usKgkCtVqt/nUQ8rrC6uqojR470v6jX19d1584dSRpav1ar6fDhw+p0OlpdXR3opmu1WgNdTqOOvbjYux/X2bNnB1oYN2/e1KlTp7b/TQIyvNp+R998/cebyr9y5NN6qfXMBCLCVlhv1lU5NBoN59ax+XS7XXU6ndTrSLbq9OnTWlhIXf0eKNz+M1f1k3NfnHQYiJjZirsP7ROnq6riarXa2FNpO51O5vRmAEgicUyBU6dOjbXkiCSuGgewZSSOKZG8CDCPIAgGliQBgFFIHACAXEgcAIBcSBwAgFxIHACAXEgcAIBcSBwAgFxIHACAXDLXqjKzQNK8pNDdF1PqzEkKJQWSuu6+GpVfkBQvqtR0929sT8gAgEkatchh5trdUWJpufuJaLstqRXtrkm6JKkT7wcAVF9m4nD3jpllrUVxXNLNxHZoZvWo1bEgaVkSa3UDwBQZd4wjUK+bKrau+4miFj0OzCx12VUzmzezZTNbvnXr1pjhAACKVsTgeCBJ7r7o7vGYRzOt5RLVa7h7Y+/evQWEAwDYTuMmjlBRoojMSOqa2ZyZJe8KtD7m3wEAlMS4ieOipAOJ7SBqYXQldRLlM+4+eINrAEAljZqO21RvllRgZl1370TlK5KOuHtoZu2ontQbEJe7r0atjpp6Yx0vFPcSAAA7aeSsKg22HOLy2cTjoXcQSisHAFQbV44DAHIhcQAAciFxAAByIXEAAHIhcQAAciFxAAByIXEAAHIhcQAAciFxAAByIXEAAHIhcQAAciFxAAByIXEAmKgr19f03Lk3JEnPnXtDV66vTTgijJK5Oi4AFOnK9TWdvfy27t77WJK0Ft7V2ctvS5KOHdo3ydCQgRYHgIk5f+1GP2nE7t77WOev3ZhQRNgKEgeAiXk3vJurHOVAVxUwRV5tv6Nvvv7jTeVfOfJpvdR6ZgIRZXsq2K21IUniqWD3BKLBVpE4gCnyUuuZfoLYf+aqfnLuixOOKNvJowcHxjgkafeuR3Xy6MEJRoVRSBwAJiYeAD9/7YbWwrvaF+zWyaMHGRgvORIHgIk6dmifjh3ap/1nruqHZz4/6XCwBZmJw8wCSfOSQndfTKkzJymUFEjquvvqsP3u3tmGeAEAEzZqVlUja2eUWFru3nH3JUkLQ/afGCdAAEC5ZLY43L1jZrWMKscl3Uxsh2ZWT7Q6GpLaY8YIAFOvSjPixh3jCNTrpoqtS5qRJDOrS1qWVM86gJnNq9cdpqeffnrMcACgmqo0I66ICwCD+IG7h6Mqu/uiuzfcvbF3794CwgEAbKdxWxyhEolCvdZGN2pFrEfdXIcl7TGzrrt3x/x7AIAJGzdxXNTggHgQjW/0Z1aZ2WFJbZIGAEyHUdNxm5JakoKoxdCJylckHXH30MzaUT1p86yqunpjHAEtDgCYDiNnVUnadP2Fu88mHi9lPH9VvcQDAJgSrI4LAMiFxAEAyIXEAQDIhcQBAMiFxAEAyIXEAQDIpVL346jSImCxKsYMAFkqlTiqtAhYrIoxA0AWuqoAALmQOAAAuZA4AAC5kDgAALmQOAAAuVRqVlVRmDILAFtH4hBTZgEgD7qqAAC5kDgAALnQVVVRjMsAmBQSR0UxLlM8kjMwHIkDSEFyvo8kiqTMxGFmgaR5SaG7L6bUmZMUSgokdd19NVHeldSQpLTnAyg/kiiSRrU4Glk7o8TScvcT0XZbUisqP+vus5JWzeyOJBIHAEyBzMTh7h0zq2VUOS7pZmI7NLN61OqYlaTo+ctjR4rKo7sDG238TOw/c1VSuT8TfI7HH+MI1Oumiq1Lmok3zKwpqS7phbQDmNm8et1hevrpp8cMB2VWVHcHJ3J1JT8TVUG3XTGD40H8IGqxSL3E8I1hlaOxj0VJajQaXkA8mHKcyMDOGvcCwFCJRKFea6ObrODuHUlnzaw+5t8CAJTAuInjoqQDie3A3VfNbN7MFhLl62P+HQBbdOX6mp4794Yk6blzb+jK9bUJR4RpM2o6blNSS1JgZt2o9SAzW5F0xN1DM2tH9SQpThYXJTUSYxwX4mm6AIpz5fqazl5+W3fvfSxJWgvv6uzltyVJxw7tm2RomCIjZ1VJ6gwpn008XhqyP0w8b9PzARTj/LUb/aQRu3vvY52/doPEgW3DIofAFHk3vJurHHgQJA5gijwV7M5VDjwIEgcwRU4ePajdux4dKNu961GdPHpwQhFhGrHIITBF4nGM89duaC28q33Bbp08epDxDWwrEgcwZY4d2qdjh/Zp/5mr+uGZz086HEwhEgeww1giBWmuXF/T+Ws3JPWuwSlra5HEAewwlkjBMFW6BofB8R3AlbwARsm6BqdsKpc4qvYlHP+KWIvm0ce/IsoeN4CdVaVrcCqVOKr4JVylXxEAJqdK1+BUKnEU+SVcVEumSr8iAExOla7BqdTgeFFfwkUOSj0V7O63kDaWA0CsStfgVKrFUVRTrsiWTJV+RQCYrGOH9vWvvfnhmc+XMmlIJWtxvP/++/rud7+buv8v/+AxXf7V47rn1i/bZa6//IP1zOeNshb+viQbUv7bsY4b+69//Jiu/d//oPDeIwp2/U5H//iuPrjxA31324Y5ntiWOHdWUTFz3OKPXcXPW1Eezve4VInjo48+0pNPPpm6/wtPSv8puKd/fueebn/g2vO46W+f+YQ+89Tvj/V39zz+W93+YPNda/c8/khmPFv1hSelL/wX6e++9xv9w5EnJD0x9jEH/WZb4txZRcXMcYs/dhU/b0V5ON/jSnVVSdJnntqlV/7q9yRJr/zV7+kzT+0a+5h/+8wufWLDO/GJR3rlAIBBpWpxTEqcfAZbMru2JSkBwLQhcUQ+81QvUfzd937Tb9EAADarXFcVAGCyStXi+MUvfqFXXnll6L5Wq6Vnn31WkvTWW2/p9vf+t175P8PDf/nll/uPX3vtNf3yl78cWu/ZZ59Vq9WSJL333nv6zne+o9s//2jocb/85S/3B6va7bbeeuutocf81Kc+pRdffLG/nXw9G4+98TW12+2hxxz1mpLHHfaa0mzHa9poq6/p9s8/kv76dOprSsrzmj76z/9d0p9u+2uSDox8TdJ4n700D/r/dPt7/5h6jozz2bu9vDb0uFX57A17TQ/6/5Q897bzNcXHncRrSviTtB2ZicPMAknzkkJ3X0ypMycplBRI6rr7aqJ8RtKspEvu3sn6W5L0+OOPj6oCANgZP03bYe6bp6H2d5o1JdWVkjiixLLg7iei7ba7t8ysLkmJJHLH3f9wVJS1Ws2/9a1vjaomqTe19dt//R+3VDePoo5b5LGLjLkoVXkvfvTuxunf2ztpgs9btU3ze/y5z31uxd0bw/ZljnFErYQwo8pxSTcT22GUNGYknUiUr8fJBKiKH717T9/+tw/71/jc/sD17X/7UD96996EIwMma9zB8UCDiWVd0oy7d+JWSGQmbn1sZGbzZrZsZsu//vWvxwwH2D7//M49ffi7wbIPf9crfxj96N17evlffitJevlffksCfYgVMasqSG6Y2YKkF9Iqu/uiuzfcvfHEE9t9RTXw4IatJpBVPs1ofSFp3MQRajBRzEjqxhvRAPk/bWVgHPnxC/C+It6LPY9vXr8sq3ya0fpC0riJ46LieYo9QWJAvC5p1d1XzaxmZrUx/xYS+AV4X1HvBUvR3EfrC0mjpuM2JbUkBWbWjVsOZrYi6Yi7h2bWjupJ0kK0vy7pdfUGxaXeGMfIWVXYuqxfgGVeKiWepST1WgbbMUupqPeCpWju2/O4pSwE+vC1vqrmf/74Q/2vm5t/RP23A7v0N5/+xAMdMzNxRIliUzeTu88mHi8N2b8qKXeieOyxx/Tee+9tsfYTOermUdRxt/fYtz8YvhT87Q9+V2D847m+/pgu//z+svi3P3D9j7f/n34V/kqHZj564OMW+V7UHpFO/pl05l+f0Mk/e1+StL1vbzU+b80/Gvy/k3q3NGj+0V2999772/I3qqn830N/8YT0F3/ee3zmX5/QuT+/PwnpQf9E5nUcO63RaPjy8vKW6u4/c1U/OffFbY+hqONu97GfO/fG0DsL7gt2928EUzZFxbwT78XD/nmTenfKrMLd6XZS1T4XeY5rZg92HQfKq4p3Fizq1r9VfC+qqCp3p0PxSrVWFbauSvcnjhV1//UqvhdAlZE4KuzYoX06dmif9p+5WtruqaSTRw/q7OW3B+7vvl0tg6q9FyjWq+139M3Xf7yp/CtHPq2XWs9MIKLpQuLAjqFlgJ3yUuuZfoIochzpYUXiwI6iZQBUH4PjAIBcSBwAgFxIHACAXEgcAJDTletreu7cG5J6F6Beub424Yh2FoPjAJDDletrA9PK18K7Onv5bUl6aGYI0uIAgBzOX7sxcC2SJN2997HOX7sxoYh2HokDmICHvaujyopaOqdK6Koq2MYrWPefuSqJK1gfZkV2dfB5K15RS+dUCYlDxZ5syStYASm7q2PcxMHnrXhFLp1TFZVKHEV9wXOyYSfR1VFtLJ1TscTBFzymAV0d1fewL53D4Diww7h/CKquUi0OYBrQ1YGqI3EAE/Cwd3Wg2jITh5kFkuYlhe6+mFJnTlIoKZDUdffVrT4XAFA9o1ocQ29UHouSQ8vdT0TbbUmtrTwXAFBNmYPj7t5RrzWR5rikm4nt0MzqW3wuAKCCxh3jCDSYHNYlzeQ5gJnNq9elpaeffnrMcACgmqp01X8Rg+NBnsrR+MeiJDUaDS8gHgAovSpdpzZu4gg1mChmJHXHPCYmaOOvnlgZf/UAmIxxE8dFSQuJ7SCeVYVqSv7q2X/mqn5y7osTjghA2YyajttUb5ZUYGbdaMBbZrYi6Yi7h2bWjupJiSSS9lwAQLVlJo7oy37TF767zyYeL+V5LgCg2lirCgCQC0uOACmqND0S2EkkDiBFlaZHAjuJxIHKo2UA7CwSByqPlgGws0gcAEaiVVdtV66v6fy1G5Kk5869Mfb9X0gcAEaqYqtuu78sq+rK9TWdvfy27t77WJK0Ft7V2ctvS9IDvx9MxwUwdeIvy/je7vGX5ZXraxOObOedv3ajnzRid+993E+qD4IWB3YM3R3YKVlflg9bq+PdKHlutXwrSBzYMVXs7kA1FfFlWVVPBbv7La+N5Q+KrqqKerX9jvafudr/1R4/frX9zoQjAyYv7UtxnC/Lqjp59KB273p0oGz3rkd18ujBBz4mLY6K4tc7kO7k0YMDA8LS+F+WVRV3zZ2/dkNr4V3tC3YzqwoANiriy7LKjh3ap2OH9mn/mav64ZnPj308EgeAqbTdX5a4jzEOAEAuJA4AQC4kDgBALiQOAEAuJA4AQC6Zs6rMLJA0Lyl098WUOnOSQkmBpK67r2aVAwCqbdR03EbWziixtNz9RLTdltRKKx87WgDAxGV2Vbl7R71WQ5rjkm4mtkMzq2eUAwAqbtwxjkCDiWVd0kxGOQCg4ooYHA/ylJvZvJktm9nyrVu3CggHALCdxk0coQYTwoykbkb5Ju6+6O4Nd2/s3bt3zHAAAEUbd62qi5IWEtuBu6+aWXdY+Zh/C5gK3NAKVTdqOm5TvdlQgZl1o8FymdmKpCPuHppZO6onRckirRzVwL2ai8WS+Ki6zMQRJYrOkPLZxOOllOcOLUe5FXFjewDThSvHMaCIG9sDmC4kDgzgXs0ARiFxYAD3agYwCncAxADu1QxkY1YciQMbcK9mIBuz4kgcGIJ7NQPIwhgHACAXEgcAIBcSBwAgFxIHACAXEgcAIBcSBwAgFxIHACAXEgcAIBcSBwAgFxIHACAXlhwBMHVYiLBYJA4AU4eFCItFVxUAIBdaHAAwxYrothuZOMxsTlIoKZDUdffVIXXmJa1LmpG0HNdJKwcA7Iwiuu0yE4eZBZJa7n4i2m5Lam2oU5c0u7FOWvm2Rg8A2HGjxjiOS7qZ2A6jhJDUHFKnllEOAKiwUYkjUK+bKhZ3OyV1JR1IbNeif2nlAIAKe5BZVUFyw92XJMnMasnWSFr5RmY2b2bLZrZ869atBwgHALCTRiWOUIOJYka9lsSAeBwjMfi9nFW+4bmL7t5w98bevXvzxA4AmIBRieOiBrubgo0zo6IWxYK7d6MxjGV3D9PKtzd8AMBOy5xVFSWAtpk1o6KFeJ+ZrUg6EiWGN6M69UQrY2g5AKDaRl7HEY9VDCmfHVKnk/LcgXIAQHWx5AgAIBcSBwAgFxIHACAXEgcAIBdWx8UAboADYBRz90nH0NdoNHx5edM1ggCAHWZmK+7eGLaPrioAQC4kDgBALiQOAEAuJA4AQC4kDgBALiQOAEAuJA4AQC4kDgBALiQOAEAuJA4AQC4kDgBALiQOAEAuJA4AQC4kDgBALiPvx2Fmc5JCSYGkrruvDqkzL2ld0oyk5biOmTWj50lS6O6dbYkaADAxmYnDzAJJLXc/EW23JbU21KlLmt1YJ3pu4O5LUfmCJBIHAFTcqK6q45JuJrbDKFEkNYfUqUWPzyYeAwCmwKjEEajXTRWLu6OSupIOJLZrkmruHko6LWnFzC65++lhf8DM5s1s2cyWb926lSN0AMAkPMg9x4PkhrsvmVkralkEG+rWJM1KWjCztru3NuyXuy9KWpQkM7tlZj/dYhyflPTvOWOftKrFXLV4perFXLV4perFXLV4pXLE/CdpO0YljlCDyWBGvRbGAHc/YWY1d181M0lajgbVl929K+kFM1sws2bWALm77x0RT5+ZLafdD7esqhZz1eKVqhdz1eKVqhdz1eKVyh/zqK6qixrshgo2zqoys5qZLbh7N2p1LEfdVDMa7OZqa0jSAQBUS2aLw91DM2tH02olaSHeZ2Yrko5ECePNqE49nl3l7otmdsrMwugp3aj1AQCosJFjHPF02iHls0PqdDbU+cZY0WVbLPDYRalazFWLV6pezFWLV6pezFWLVyp5zObuk44BAFAhLDkCAMiFxAEAyOVBruPYUdHSJfPqrXW1mCgfuYbWpIyIeUa9a1sulWXtrrR4E/vnVLK1xrJijtZO60pSWWLO+EyUej23tM9syc+/rJjLeP5lxlXG86/0iUPSprnMW1lDa8KGxVxX7wSL1+66I+kPdzqwFKnzxaP3+oQSM+pKYmjMZnZJ0t9HMwIvqTzro6V9jku7nlvaZ7bM519GzKU8/0bFVdbzr/RdVVGWDTcUb2UNrYlJiXlGvQ9AbL0sMafEG2uodw1OqQyLOX4/o6RRc/cXJhHbMBnvcZnXc0v7zJb5/EuLuazn36i4Snn+VaHFMUyg0WtolUr0xZH8NTlTpub9MNEHeFlSGU6wrWhI93/JRxemDl0jrQyiBBev59YpU6KT0j+zUfdamCgvzfk34jwr3fmXFW+Zz7/StzhyCCYdwFZFXRKl+pJIE60CUBWB1Is5OvnqJf41H4vXc4u7fEppC5/ZYIdC2bK0mMt6/g2Lq6znX1VbHKG2sIZWGUUDXf9Uhl87WeKbc0VfvIcl7TGzsl/935W0J7EdqvfFXMqYH2Q9t0kY8pkNVfLzL+08K+v5tzGusp9/VU0cFzU4WLRpDa0yipqeq4l1vVSWD8JGG2b+HJbULmusCR0NDtLW1Gvql9WMpOTntnTruQ37zKrk51/aeVbW8y8lrlKff6VPHFF/aku9Puuuu3ey1tAqg2ExRx+O19X7FSH1+jInPqtDGh5vYl9dvT7WoEy/eEZ8Luajal8vS1M/Jd5Sr+eW9Zkt6/mXFnNZz79RcZX2/GPJEQBAHtM0OA4A2AEkDgBALiQOAEAuJA4AQC4kDmAHRFezA1OBxAHsjLOTDgDYLiQOoGBRa+P2pOMAtguJAyjecUlLkw4C2C4kDqB4B8pyxS+wHUq/5AgwadH6QXX1Fpu7ED2uufs3tvjcm0PKT+n+ulQn3L0UN0ICtoIWBzBaM7pD25uSLkSP+zffMbN5M2sm1shKmlNvUUAl6i/o/l3fVlWSe1kAW0XiAEaLv/gPS7okSe5+QOovh72euK91c8Nz9yQXWoxaIHPxrULVa72Uahl1YBQSBzBC4ot/Ll7uOrHE+GHdvxteV4ll3aOVTTfenGljomgNqQOUGokDyBAtez4fJYEwKqupd68PqXdDo/XEU4LE4y+l3JQpTDw+Hi27P7ddMQNFY3AcyLaq+8ngQvwFn+hqCjU4RhFmHczdl8zscKJL6+vRMUtzIyRgFO7HAYwhmUiiwfFuogXRLdOd8YDtQlcVMIao5TETtyASXVOHSRqYVnRVAWNK3h9aYokRTD+6qgAAudBVBQDIhcQBAMiFxAEAyIXEAQDIhcQBAMiFxAEAyIXEAQDI5f8DFkYvlNXAd1IAAAAASUVORK5CYII=", "text/plain": [ - "
" + "
" ] }, "metadata": { @@ -261,35 +259,35 @@ } ], "source": [ - "# We do not just have V1V1, but also the two other spacial directions. We can average over them for better statistics. \n", + "# We do not just have V1V1, but also the two other spacial directions. We can average over them for better statistics.\n", "matrix_V2V2= pe.input.json.load_json(\"./data/matrix_correlator_V2V2\", verbose=False)\n", "matrix_V3V3= pe.input.json.load_json(\"./data/matrix_correlator_V3V3\", verbose=False)\n", - "matrix_VnVn=(matrix_V1V1 + matrix_V2V2 + matrix_V3V3) / 3. \n", + "matrix_VnVn=(matrix_V1V1 + matrix_V2V2 + matrix_V3V3) / 3\n", "\n", - "#We then solve the GEVP to get eigenvectors corresponding to the ground state. \n", + "# We then solve the GEVP to get eigenvectors corresponding to the ground state.\n", "\n", - "vec_ground = matrix_VnVn.GEVP(t0=3, ts=6, state=0)\n", + "vec_ground = matrix_VnVn.GEVP(t0=3, ts=6, state=0, sorted_list=None)\n", "\n", - "#Now we project the matrix-correlators to get new correlators belonging to the ground state.\n", + "# Now we project the matrix-correlators to get new correlators belonging to the ground state.\n", "\n", "corr_ground = matrix_VnVn.projected(vec_ground)\n", "\n", - "# We get the effective mass using the periodic cosh method. \n", + "# We get the effective mass using the periodic cosh method.\n", "m_eff_Jpsi = corr_ground.m_eff(variant=\"cosh\")\n", "\n", - "# We can now pick a plateau range and get a single value for the mass. \n", + "# We can now pick a plateau range and get a single value for the mass.\n", "m_eff_Jpsi.gamma_method()\n", - "m_Jpsi = m_eff_Jpsi.plateau([5, 24],method=\"avg\")\n", + "m_Jpsi = m_eff_Jpsi.plateau([5, 24], method=\"avg\")\n", "m_Jpsi.gamma_method()\n", "# We can now visually compare our plateau value to the data\n", "\n", "m_eff_Jpsi.show([10, 25], plateau=m_Jpsi)\n", "\n", - "# Since the lattice spacing is known, we can multiply with hbar*c/a to see it in physical units \n", + "# Since the lattice spacing is known, we can multiply with hbar*c/a to see it in physical units\n", "\n", "m_Jpsi = m_Jpsi * 197 / 0.0653\n", "\n", - "#As a last step we call the gamma method to get the error \n", + "# As a last step we call the gamma method to get the error\n", "\n", "m_Jpsi.gamma_method()\n", "\n", @@ -299,7 +297,7 @@ { "cell_type": "code", "execution_count": null, - "id": "c55c6066", + "id": "f819143b", "metadata": {}, "outputs": [], "source": [] @@ -321,7 +319,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.8" + "version": "3.8.10" } }, "nbformat": 4,