From 68a7015ab1996c994736b5fa0a25764ef070b0a7 Mon Sep 17 00:00:00 2001 From: samwaseda Date: Thu, 14 Mar 2024 21:03:28 +0000 Subject: [PATCH] update in the text --- notebooks/first_steps.ipynb | 60 ++++++++++++++++++++++++++++++------- 1 file changed, 50 insertions(+), 10 deletions(-) diff --git a/notebooks/first_steps.ipynb b/notebooks/first_steps.ipynb index 171060f..539e7c5 100644 --- a/notebooks/first_steps.ipynb +++ b/notebooks/first_steps.ipynb @@ -182,7 +182,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "id": "12103d66", "metadata": {}, "outputs": [], @@ -201,10 +201,21 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 11, "id": "d8f7f333", "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAGwCAYAAACq12GxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAABCLElEQVR4nO3deXxU9b3/8fdsmSwkgwSyAGFTCCCIGIqA5WpdokCl3qsF5SpiReWnVpQuSmndb9FaqXUBtQpai4i4PbyPS9VYLaAISgxuICAgYUkICWayLzNzfn8kGRISIBNm5mQyr+fjMY9kvvM9k8+cWuft93y/32MxDMMQAABAlLKaXQAAAICZCEMAACCqEYYAAEBUIwwBAICoRhgCAABRjTAEAACiGmEIAABENbvZBXR2Pp9PBw4cUGJioiwWi9nlAACAdjAMQ+Xl5erdu7es1uOP/RCGTuDAgQPKyMgwuwwAANABe/fuVd++fY/bhzB0AomJiZIaTmZSUpLJ1QAAgPYoKytTRkaG/3v8eAhDJ9B0aSwpKYkwBABAhGnPFBcmUAMAgKhGGAIAAFGNMAQAAKIaYQgAAEQ1whAAAIhqhCEAABDVCEMAACCqEYYAAEBUIwwBAICoRhgCAABRjTAEAACiGmEIAABENW7UapKqOo9KKupkt1lks1hks7bxaGxvz03mAABAxxCGTLJm2yH9v+Wft6uv1aJWAanhYZXNKtmtVlmbflqanltkt1r8P48buBoDWfP+/p+N7Ucf57BZZbdaZPf/tMhhtTaEu2avO2xW2Zpeb+rb2M9ha/jd399mafGeDlvD5yEMAgBCiTBkoliHVT6f5PH55DOO3c9nSD6voXrvcTp1YU2hqXlYchwVxOxWqxyNQcxusx45plkQi7FbFdP48+jnMY3HxNhtjT+tcjbr12b/pudNbfaGYAcAiCyEIZNMGpmuSSPT/c8Nw5DXZ8jb9PPoh2HI4zXkMwx5fIZ8voafzV9v+t3/mmHI6z3Bezbr7//Z+LcajvPJ61OLnx5fw+v1Pp+8Tb97G9t9hjxeX4vX672NbT5DHl/ja96G372N/TzehmPbUu81VO/1SvXh+l+n46wW+cOT86jQ5GgWmo6EsiMBzHlUP6fdpvgYm+JiGn86mn63Kz7GplhHQ3tTnxiblVE0AOgAwlAnYbE0jnCYXYiJmgKhx9cYrhoD0pEA1Sxc+Xyq9zb29/pU72sWuFq87msMUz7VeXz+n3Ve46jnjY/mbY2/1za1+V83jhzj8bX4DD5Dqqn3qabep/Iwnz+b1dIsMNmO+t3eqq0hRNkV5zh+6IprbHPYWG8BoGuK5u9edDL+QGiTYh02s8tpF6NxdO344amprWW/5oGqrVBWU+9TdZ1H1fVeVdV5VV3X8LOm8XlV42tNl0+9PkMVtR5V1HpC8lkdNktjcLL7A1LzsBQfY1O3WLtOiY+RK86hU+JjdEqCQ93jY9S98XlSnINLiQA6HcIQcBIsloaJ4A6bVQlOc2qo9/pUXX8kLFXXeVVd72kMTE3Pm35vaG/e/0jAavladZ1XVfVeeRsvXzaMsHlUVtPxsGWxSEmxDp0S3xiS4htCUvOf3eNjGl6Pa2xPiFFCjI1LgABChjAERDhH4zyjpFhH0N/bMAzVeX2qqfOpqjFg+UNX/ZFw1RSoymo8Kq2qU2lVvX5o9tNdVa/yWo8MQ3JX18tdXS+VVAXwGS1yxTWEpFPiY+SKd/h/PxKqWgcspz0yRhgBmIswBOCYLBaLnHabnHabXDq5sFXv9am0ql7u6jr9UFWvHyobwlJp4/PSqjr9UNnwvClE/VBV75+nVVxRq+KK2oD+ZnyMTd3jGkebEpqNNjUbhUpNciozNVG9Ep2MPgFRijAEICwcNqt6JTrVKzGw64nVdd7GYNQwwvSDf9SpKTQ1BKnS6iOjUaVVdfIZ8o9aHXDXnPDvdI93KDM1UZlpjY/URA1JSwzJiBuAzoUwBKBTi4uxKS4mTr27x7X7GJ/PUHltwyU7f1hqNtrkbvz5Q1Wd9pdW6/viSpVW1Wvj7sPauPtwi/fq0z1OQ1K7KTMtSZlp3ZSZmqRTUxK4BAd0IYQhAF2O1WqRK84hV5xD/ZNP3L+m3qvviiq0/WC5thWWa1vjzwJ3jfaXVmt/abU+3HbI399mtWhgz4RWI0n9esTLymo5IOJYDMOIzm2N26msrEwul0tut1tJSUlmlwMgjNzV9UcCUuPj28KyY66oi3PYNDi1W6uQxHwkIPwC+f4mDJ0AYQhAc4Zh6GBZbePoUZm2FVZo28Ey7ThYodqjNuFsckq8Q0NSEzU0LdF/uW1IaqISmY8EhAxhKIgIQwDaw+sztKekssVltm0Hy/V9ceUx7z3Yp3ucMtMSmwWlRA3qxXwkIBgIQ0FEGAJwMprmI7UISYXlKixre4WbvXE+0pC0RA1tXNE2NC1RGacwHwkIBGEoiAhDAELBXVXfEI78l9saQtLx5iM1rGpLVFb/U3T5WX1l535xwDERhoKIMAQgXAzDUGFZzZEJ240jSTuKKlrdFPjSUb21aNoobqALHEMg398RtbR+7dq1euSRR5Sbm6uCggK9+eabuuyyy47Z/4033tCSJUu0efNm1dbW6vTTT9e9996riy++OHxFA0A7WSwWpbvilO6K03mZKf52j9enPYertK2wXF/td+u5dbv0v18ckMfr01+vHK0YO4EIOBkR9f+gyspKjRo1Sk8++WS7+q9du1YXXXSRVq9erdzcXP3kJz/RpZdeqry8vBBXCgDBY7dZdWqvbpo8Ml13XjJUT1+dpRibVf/8ulC3vPy5aj1es0sEIlrEXiazWCwnHBlqy+mnn67p06fr7rvvbld/LpMB6Iz+va1IN76UqzqPTz/J7KUlV2cp1sEqNKBJIN/fETUydLJ8Pp/Ky8vVo0ePY/apra1VWVlZiwcAdDbnZaZo6bU/UqzDqg+3HdINf9+kmnpGiICOiKow9Oijj6qyslLTpk07Zp+FCxfK5XL5HxkZGWGsEADa78eDe2rZrLGKj7Fp3Y5i/eKFz1RV1/ZqNADHFjVhaMWKFbr33nu1cuVKpaSkHLPf/Pnz5Xa7/Y+9e/eGsUoACMz4U5P14i/GKiHGpvU7SzRr2WeqqCUQAYGIijC0cuVKXX/99Xr11Vd14YUXHrev0+lUUlJSiwcAdGY/GtBDL80+W4lOuz7dfVjXLv1U5TX1ZpcFRIwuH4ZWrFihWbNm6eWXX9aUKVPMLgcAQuKsfqfoH7PPVlKsXbl7ftDVz38qdzWBCGiPiApDFRUV2rx5szZv3ixJ2r17tzZv3qz8/HxJDZe4Zs6c6e+/YsUKzZw5U48++qjGjRunwsJCFRYWyu12m1E+AITUqIzuevmGceoe79AXe0t19XMbVVpVZ3ZZQKcXUWFo06ZNGj16tEaPHi1JmjdvnkaPHu1fJl9QUOAPRpL0zDPPyOPx6JZbblF6err/MXfuXFPqB4BQG9HHpRU3jFOPhBh9td+tq/62UYcrCUTA8UTsPkPhwj5DACLR9oPlmvG3jSquqFVmaqKW33C2enZzml0WEDbsMwQAUW5IaqJeuXGcUhKd2nawXFc+u0FFZTVmlwV0SoQhAOiiTkvpppU3jVe6K1bfFVXoymc3qNBNIAKORhgCgC5sYM8ErbxxvPp0j9Ou4kpNf/YT7S+tNrssoFMhDAFAF9cvOV4rbxqnjB5x2lNSpenPfKK9h6vMLgvoNAhDABAF+p4Sr1dvGq8ByfHa90O1pj/zifaUVJpdFtApEIYAIEqku+K08qbxGtQrQQfcNZr+zAbtOlRhdlmA6QhDABBFUpNitfLG8Rqc0k2FZTWa/uwGfVdUbnZZgKkIQwAQZXolOvXKjeM0NC1Rh8prdeWzG7StkECE6EUYAoAolNzNqRU3jNPpvZNUXFGnK5/9RN8c4FZFiE6EIQCIUqckxOjl2eM0qq9LP1TVa8bfNuqrfQQiRB/CEABEMVe8Qy/NPluj+3WXu7peM57boLz8H8wuCwgrwhAARLmkWIf+/oux+tGAU1Re49E1z3+q3D2HzS4LCBvCEABAibEOvXDdWI0b1EMVtR7NfP5TbdxVYnZZQFgQhgAAkqQEp13LZo3Vj0/rqco6r2Yt+0zrvys2uywg5AhDAAC/uBibnrt2jM4d0kvV9V5d98JnWrv9kNllASFFGAIAtBDrsOnZmVm6YGiKaj0+zf77Jn34bZHZZQEhQxgCALTitNu05OosXXx6quo8Pt30Uq5ythw0uywgJAhDAIA2xditenLGWZoyMl11Xp/+3z9y9c7XBWaXBQQdYQgAcEwOm1V/vfJM/ezM3vL4DN3ycp7+94sDZpcFBBVhCABwXHabVYumnan/OquPvD5Dc1/J01t5+80uCwgawhAA4IRsVoseuWKUpo/JkM+Q7nh1s17L3Wd2WUBQEIYAAO1is1q08L9G6r/P7ifDkH7z2hd65dN8s8sCThphCADQblarRQ9eNkKzJgyQYUh3vfGVXtqwx+yygJNCGAIABMRiseieS4dr9o8HSpL+8NbXWvbxbpOrAjqOMAQACJjFYtGCKcM059xTJUn3/e8W/W3tLpOrAjqGMAQA6BCLxaI7L8nUbeefJkn6n9Vb9dSH35lcFRA4whAAoMMsFovmZWdq3kVDJEmPvLtNf31/h8lVAYEhDAEATtptFwzWby/JlCT95f3tevS9bTIMw+SqgPYhDAEAguLm807T76cMkyQ98cF3euidbwlEiAiEIQBA0MyeOEj3XjpckvTMml168P+2EojQ6RGGAABBNeucgXrwshGSpOc/2q173/6GQIROjTAEAAi6q8f118OXj5TFIr34yR4teOtr+XwEInROhCEAQEhM/1E/PXLFKFks0ssb8/XelkKzSwLaRBgCAITMFVl9deWP+kmS8vJLzS0GOAbCEAAgpEb0SZIkbTtYbnIlQNsIQwCAkMpMTZQkbS8kDKFziqgwtHbtWl166aXq3bu3LBaL3nrrrRMes2bNGmVlZSk2NlaDBg3S008/HfpCAQB+gxvD0AF3jcpq6k2uBmgtosJQZWWlRo0apSeffLJd/Xfv3q3Jkydr4sSJysvL0+9+9zvddtttev3110NcKQCgiSvOoXRXrCRpB5fK0AnZzS4gEJMmTdKkSZPa3f/pp59Wv3799Nhjj0mShg0bpk2bNunPf/6zLr/88hBVCQA42pDURBW4a7StsEJZ/XuYXQ7QQkSNDAXqk08+UXZ2dou2iy++WJs2bVJ9fdtDtbW1tSorK2vxAACcnMy0xnlDjAyhE+rSYaiwsFCpqakt2lJTU+XxeFRcXNzmMQsXLpTL5fI/MjIywlEqAHRpQxrnDW1jEjU6oS4dhiTJYrG0eN60JfzR7U3mz58vt9vtf+zduzfkNQJAV+dfUcbIEDqhiJozFKi0tDQVFrbc8bSoqEh2u13JycltHuN0OuV0OsNRHgBEjdNSuslikUoq61RcUaue3fj3LDqPLj0yNH78eOXk5LRoe++99zRmzBg5HA6TqgKA6BMXY1P/HvGS2G8InU9EhaGKigpt3rxZmzdvltSwdH7z5s3Kz8+X1HCJa+bMmf7+c+bM0Z49ezRv3jxt3bpVS5cu1fPPP69f//rXZpQPAFHNP2+IS2XoZCIqDG3atEmjR4/W6NGjJUnz5s3T6NGjdffdd0uSCgoK/MFIkgYOHKjVq1fr3//+t84880w98MADevzxx1lWDwAmYEUZOquImjN03nnn+SdAt+WFF15o1Xbuuefq888/D2FVAID2YEUZOquIGhkCAESuIyNDFcf9D1sg3AhDAICwGJCcIIfNoopajw64a8wuB/AjDAEAwiLGbtWgnt0ksaIMnQthCAAQNkPSWFGGzocwBAAIm8xURobQ+RCGAABhw15D6IwIQwCAsGlaUbajqEJeHyvK0DkQhgAAYZNxSrxiHVbVeXzaU1JpdjmAJMIQACCMrFaL/1IZO1GjsyAMAQDC6shO1BUmVwI0IAwBAMIqk5EhdDKEIQBAWLHXEDobwhAAIKyaRoZ2F1eq1uM1uRqAMAQACLPUJKeSYu3y+gztOsSKMpiPMAQACCuLxdLsDvZcKoP5CEMAgLA7sqKMMATzEYYAAGHHyBA6E8IQACDsuEcZOhPCEAAg7JrC0N7D1aqs9ZhcDaIdYQgAEHY9EmLUK9EpqeGmrYCZCEMAAFP4d6JmEjVMRhgCAJiCeUPoLAhDAABTZKZ1k8SKMpiPMAQAMAV7DaGzIAwBAEwxuDEMFZXX6ofKOpOrQTQjDAEATNHNaVffU+IkcakM5iIMAQBM419RRhiCiQhDAADTDEljRRnMRxgCAJgmk0nU6AQIQwAA0zRfUWYYhsnVIFoRhgAAphnUK0E2q0VlNR4dLKs1uxxEKcIQAMA0sQ6bBiTHS2LeEMxDGAIAmCozjXuUwVyEIQCAqbhHGcxGGAIAmIq9hmA2whAAwFRNew1tP1gun48VZQg/whAAwFT9e8Qrxm5VTb1Pe3+oMrscRKGIC0OLFy/WwIEDFRsbq6ysLK1bt+64/ZcvX65Ro0YpPj5e6enpuu6661RSUhKmagEAJ2K3WXVar26S2HwR5oioMLRy5UrdfvvtWrBggfLy8jRx4kRNmjRJ+fn5bfb/6KOPNHPmTF1//fX65ptvtGrVKn322WeaPXt2mCsHABxPZhrzhmCeiApDixYt0vXXX6/Zs2dr2LBheuyxx5SRkaElS5a02X/Dhg0aMGCAbrvtNg0cOFA//vGPddNNN2nTpk1hrhwAcDxHVpRVmFwJolHEhKG6ujrl5uYqOzu7RXt2drbWr1/f5jETJkzQvn37tHr1ahmGoYMHD+q1117TlClTjvl3amtrVVZW1uIBAAitzLSGy2TsNQQzREwYKi4ultfrVWpqaov21NRUFRYWtnnMhAkTtHz5ck2fPl0xMTFKS0tT9+7d9cQTTxzz7yxcuFAul8v/yMjICOrnAAC01jQytPNQheo8PpOrQbSJmDDUxGKxtHhuGEartiZbtmzRbbfdprvvvlu5ubl65513tHv3bs2ZM+eY7z9//ny53W7/Y+/evUGtHwDQWp/ucermtMvjM/R9SaXZ5SDK2M0uoL169uwpm83WahSoqKio1WhRk4ULF+qcc87Rb37zG0nSGWecoYSEBE2cOFEPPvig0tPTWx3jdDrldDqD/wEAAMdksVg0JLWbPs8v1bbCcv9IERAOETMyFBMTo6ysLOXk5LRoz8nJ0YQJE9o8pqqqSlZry49os9kkNYwoAQA6D1aUwSwRE4Ykad68eXruuee0dOlSbd26VXfccYfy8/P9l73mz5+vmTNn+vtfeumleuONN7RkyRLt2rVLH3/8sW677TaNHTtWvXv3NutjAADa4F9RxiRqhFnEXCaTpOnTp6ukpET333+/CgoKNGLECK1evVr9+/eXJBUUFLTYc2jWrFkqLy/Xk08+qV/96lfq3r27zj//fD388MNmfQQAwDFwjzKYxWJwvei4ysrK5HK55Ha7lZSUZHY5ANBlFVfUasyD78tikbbcd4niYmxml4QIFsj3d0RdJgMAdF09uzmVnBAjw5C+K2LzRYQPYQgA0Gkc2YmaS2UIH8IQAKDTYEUZzEAYAgB0GqwogxkCDkPLli3TqlWrWrWvWrVKL774YlCKAgBEJ/89yhgZQhgFHIYeeugh9ezZs1V7SkqK/vjHPwalKABAdBrcODJU4K6Ru7re5GoQLQIOQ3v27NHAgQNbtffv37/FHj8AAAQqKdah3q5YSdIORocQJgGHoZSUFH355Zet2r/44gslJycHpSgAQPQaksaKMoRXwGHoyiuv1G233aYPP/xQXq9XXq9XH3zwgebOnasrr7wyFDUCAKKIfydqJlEjTAK+HceDDz6oPXv26IILLpDd3nC4z+fTzJkzmTMEADhp7DWEcAs4DMXExGjlypV68MEHtXnzZsXFxWnkyJH++4MBAHAymvYa2lZYLsMwZLFYTK4IXV2Hb9Q6ePBgDR48OJi1AACg01K6yWKRfqiqV3FFnXolOs0uCV1cwHOGrrjiCj300EOt2h955BH9/Oc/D0pRAIDoFeuwaUBygiT2G0J4BByG1qxZoylTprRqv+SSS7R27dqgFAUAiG5DUhs2X2QnaoRDwGGooqJCMTExrdodDofKysqCUhQAILr5V5QxMoQwCDgMjRgxQitXrmzV/sorr2j48OFBKQoAEN3YawjhFPAE6j/84Q+6/PLLtXPnTp1//vmSpH/9619asWJFm/csAwAgUM33GmJFGUIt4DA0depUvfXWW/rjH/+o1157TXFxcTrjjDP0/vvv69xzzw1FjQCAKDOgZ4IcNosq67zaX1qtvqfEm10SurB2XSZ7/PHHVVNTI0nKz8/X5MmT9fHHH6uyslLFxcX64IMPCEIAgKBx2Kw6tRd3sEd4tCsMzZs3zz85euDAgTp06FBIiwIAwL8TdWGFyZWgq2vXZbLevXvr9ddf1+TJk2UYhvbt2+cfKTpav379glogACA6ZaYlSl8wMoTQa1cY+v3vf69f/vKXuvXWW2WxWPSjH/2oVZ+mCW5erzfoRQIAos+RkSHCEEKrXWHoxhtv1FVXXaU9e/b4J0snJyeHujYAQBRrWlH23aEKebw+2W0B7wYDtEu7V5MlJiZqxIgRWrZsmc455xw5nce/V8yKFSs0depUJSQknHSRAIDo0/eUOMU5bKqu92rP4Sr/hGog2AKO2ddee+0Jg5Ak3XTTTTp48GCHigIAwGq1+G/LsZ1LZQihkI05GoYRqrcGAEQJ/7whJlEjhLgACwDotDLTuEcZQo8wBADotFhRhnAgDAEAOq2mkaHvS6pUU8/WLQgNwhAAoNNKSXTKFeeQ12do16FKs8tBFxWyMNS/f385HI5QvT0AIApYLJYjd7Bn3hBCpN1h6NNPP22xu/TRq8Vqa2v16quv+p9//fXXysjICEKJAIBoNiStYXk9K8oQKu0OQ+PHj1dJSYn/ucvl0q5du/zPS0tLddVVVwW3OgBA1POPDDGJGiHS7jB09EhQW/sIsbcQACDY2GsIoRbUOUMWiyWYbwcAgD8M7fuhWhW1HpOrQVfEajIAQKd2SkKMUhIbbgO1g9EhhEBAYWjLli368ssv9eWXX8owDH377bf+5998802oamxh8eLFGjhwoGJjY5WVlaV169Ydt39tba0WLFig/v37y+l06tRTT9XSpUvDUisAIDjYiRqh1O671kvSBRdc0GJe0E9/+lNJDZfHDMMI+WWylStX6vbbb9fixYt1zjnn6JlnntGkSZO0ZcsW9evXr81jpk2bpoMHD+r555/XaaedpqKiInk8DLMCQCQZkpqodTuKta2wwuxS0AW1Owzt3r07lHW0y6JFi3T99ddr9uzZkqTHHntM7777rpYsWaKFCxe26v/OO+9ozZo12rVrl3r06CFJGjBgQDhLBgAEAXsNIZTaHYb69+8fyjpOqK6uTrm5ubrrrrtatGdnZ2v9+vVtHvP2229rzJgx+tOf/qSXXnpJCQkJmjp1qh544AHFxcW1eUxtba1qa2v9z8vKyoL3IQAAHTIkjRVlCJ2ALpOZqbi4WF6vV6mpqS3aU1NTVVhY2OYxu3bt0kcffaTY2Fi9+eabKi4u1s0336zDhw8fc97QwoULdd999wW9fgBAxw1Oadh48VB5rQ5X1qlHQozJFaEribjVZEfPSzreXCWfzyeLxaLly5dr7Nixmjx5shYtWqQXXnhB1dXVbR4zf/58ud1u/2Pv3r1B/wwAgMAkOO3K6NEwos+lMgRbxIShnj17ymaztRoFKioqajVa1CQ9PV19+vSRy+Xytw0bNkyGYWjfvn1tHuN0OpWUlNTiAQAwH/OGECoRE4ZiYmKUlZWlnJycFu05OTmaMGFCm8ecc845OnDggCoqjqw+2L59u6xWq/r27RvSegEAweXfiZrbciDITioMPfTQQyotLQ1SKSc2b948Pffcc1q6dKm2bt2qO+64Q/n5+ZozZ46khktcM2fO9PefMWOGkpOTdd1112nLli1au3atfvOb3+gXv/jFMSdQAwA6J/YaQqic1ATqP/7xj5o2bZq6d+8epHKOb/r06SopKdH999+vgoICjRgxQqtXr/avdCsoKFB+fr6/f7du3ZSTk6Nf/vKXGjNmjJKTkzVt2jQ9+OCDYakXABA8zUeGwrG3HaKHxTiJu6smJibqiy++0KBBg4JZU6dSVlYml8slt9vN/CEAMFGtx6vhd78rr8/QhvkXKM0Va3ZJ6MQC+f6OmDlDAIDo5rTbNLBngiT2G0JwnVQY2rJli+mbMQIAood/RRmTqBFEJxWGMjIyZLPZglULAADH5Z83xMgQgojLZACAiJGZ1rATNSvKEEyEIQBAxBjSbONFn6/D63+AFghDAICI0T85QTF2q2rqfdr7Q5XZ5aCLIAwBACKGzWrx37SVnagRLAGHoSuuuEIPPfRQq/ZHHnlEP//5z4NSFAAAx8I9yhBsAYehNWvWaMqUKa3aL7nkEq1duzYoRQEAcCxD0ppWlFWcoCfQPgGHoYqKCsXExLRqdzgcKisrC0pRAAAcC3sNIdgCDkMjRozQypUrW7W/8sorGj58eFCKAgDgWJpGhnYeqlCdx2dyNegKAr5R6x/+8Addfvnl2rlzp84//3xJ0r/+9S+tWLFCq1atCnqBAAA019sVq25OuypqPfq+pNK/3B7oqIBHhqZOnaq33npL3333nW6++Wb96le/0r59+/T+++/rsssuC0GJAAAcYbFYNCS1YUXZt1wqQxAEPDIkSVOmTGlzEjUAAOGQmZaoz/NLG+YNjTK7GkS6gEeGPvvsM23cuLFV+8aNG7Vp06agFAUAwPFwjzIEU8Bh6JZbbtHevXtbte/fv1+33HJLUIoCAOB42GsIwRRwGNqyZYvOOuusVu2jR4/Wli1bglIUAADH07SiLP9wlarqPCZXg0gXcBhyOp06ePBgq/aCggLZ7R2aggQAQEB6dnMqOSFGhiF9V8Tmizg5AYehiy66SPPnz5fb7fa3lZaW6ne/+50uuuiioBYHAMCx+OcNsaIMJyngoZxHH31U//Ef/6H+/ftr9OjRkqTNmzcrNTVVL730UtALBACgLZlpifpkVwnzhnDSAg5Dffr00Zdffqnly5friy++UFxcnK677jpdddVVcjgcoagRAIBWjqwo4zIZTk6HJvkkJCToxhtvDHYtAAC0W2Zaw8aL3KMMJ6tDYWjbtm164okntHXrVlksFg0dOlS33nqrhg4dGuz6AABo0+DGkaHCshq5q+rliufqBDom4AnUr732mkaMGKHc3FyNGjVKZ5xxhj7//HONHDmSe5MBAMImKdah3q5YSdL2IkaH0HEBjwz99re/1fz583X//fe3aL/nnnt055136uc//3nQigMA4HiGpCXqgLtG2wrL9aMBPcwuBxEq4JGhwsJCzZw5s1X71VdfrcLCwqAUBQBAe7ATNYIh4DB03nnnad26da3aP/roI02cODEoRQEA0B7sNYRgaNdlsrffftv/+9SpU3XnnXcqNzdX48aNkyRt2LBBq1at0n333ReaKgEAaENm2pGRIcMwZLFYTK4IkchiGIZxok5Wa/sGkCwWi7xe70kX1ZmUlZXJ5XLJ7XYrKSnJ7HIAAM3U1Hs17O53ZBjSpwsuUEpirNkloZMI5Pu7XSnH5/O169HVghAAoHOLddg0IDlBkrS9kM0X0TEBzxlqr5EjR2rv3r2hensAACRJQ1IbNl/cxiRqdFDIwtD333+v+vr6UL09AACSmq0oYxI1OihkYQgAgHAYktZ0jzLCEDqGMAQAiGhNI0M7DpbL5zvhmiCgFcIQACCiDeiZIIfNoso6r/aXVptdDiIQYQgAENEcNqtO7dV4B3sulaEDCEMAgIjn34maMIQOCFkYeuaZZ5Samhr09128eLEGDhyo2NhYZWVltXlrkLZ8/PHHstvtOvPMM4NeEwDAXP6dqFlRhg5odxiaPHmy3G63//n//M//qLS01P+8pKREw4cP9z+fMWOGEhISglNlo5UrV+r222/XggULlJeXp4kTJ2rSpEnKz88/7nFut1szZ87UBRdcENR6AACdw5GRITZeRODaHYbeffdd1dbW+p8//PDDOnz4sP+5x+PRtm3bglvdURYtWqTrr79es2fP1rBhw/TYY48pIyNDS5YsOe5xN910k2bMmKHx48eHtD4AgDmGNo4M7SyqkMfrM7kaRJp2h6Gjb2HWjluaBVVdXZ1yc3OVnZ3doj07O1vr168/5nHLli3Tzp07dc8997Tr79TW1qqsrKzFAwDQufXpHqf4GJvqvD59X1JldjmIMBEzgbq4uFher7fVPKTU1FQVFha2ecyOHTt01113afny5bLb7e36OwsXLpTL5fI/MjIyTrp2AEBoWa0WDU49cgd7IBDtDkMWi0UWi6VVW7gd/TcNw2izDq/XqxkzZui+++7TkCFD2v3+8+fPl9vt9j+4vxoARIbMpnuUMYkaAWrfcIkaQsesWbPkdDolSTU1NZozZ45/knTz+USh0LNnT9lstlajQEVFRW2uWisvL9emTZuUl5enW2+9VZLk8/lkGIbsdrvee+89nX/++a2Oczqd/s8IAIgcQxgZQge1Owxde+21LZ5fffXVrfrMnDnz5Cs6hpiYGGVlZSknJ0f/+Z//6W/PycnRz372s1b9k5KS9NVXX7VoW7x4sT744AO99tprGjhwYMhqBQCEXyb3KEMHtTsMLVu2LJR1tMu8efN0zTXXaMyYMRo/fryeffZZ5efna86cOZIaLnHt379ff//732W1WjVixIgWx6ekpCg2NrZVOwAg8jXdo+z74krV1HsV67CZXBEiRbvDUGcwffp0lZSU6P7771dBQYFGjBih1atXq3///pKkgoKCE+45BADomnolOtU93qHSqnrtPFSh03u7zC4JEcJihHuNfIQpKyuTy+WS2+1WUlKS2eUAAI5j2jOf6NPdh/WX6aP0n6P7ml0OTBTI93fELK0HAOBEmi6VbStkJ2q0H2EIANBlDEljRRkCRxgCAHQZR0aGCENoP8IQAKDLGNK48eL+0mqV19SbXA0iBWEIANBldI+PUWpSw8a5O4qYN4T2IQwBALoU/07UXCpDOxGGAABdin/eEJOo0U6EIQBAl8KKMgSKMAQA6FLYawiBIgwBALqUwY0ryooralVSUWtyNYgEhCEAQJcSH2NXvx7xkqTtBxkdwokRhgAAXY5/RRnzhtAOhCEAQJeTmdZwqYwVZWgPwhAAoMthryEEgjAEAOhyMtOO7DVkGIbJ1aCzIwwBALqcQT27yW61qLzGo8KyGrPLQSdHGAIAdDkxdqsG9kyQxB3scWKEIQBAl8RO1GgvwhAAoEtiJ2q0F2EIANAlsdcQ2oswBADokppWlO0oKpfXx4oyHBthCADQJfXrES+n3aqaep/2Hq4yuxx0YoQhAECXZLNa/DdtZSdqHA9hCADQZbETNdqDMAQA6LL8K8oYGcJxEIYAAF0Wew2hPQhDAIAuq2lkaNehStV5fCZXg86KMAQA6LLSXbFKdNrl8RnaXVxpdjnopAhDAIAuy2Kx+C+VMW8Ix0IYAgB0aawow4kQhgAAXVomew3hBAhDAIAujRVlOBHCEACgS2taUZZ/uEpVdR6Tq0FnRBgCAHRpyd2c6tktRoYhfVdUYXY56IQIQwCALq9pEvU2JlGjDYQhAECX519RxrwhtCHiwtDixYs1cOBAxcbGKisrS+vWrTtm3zfeeEMXXXSRevXqpaSkJI0fP17vvvtuGKsFAHQGmf69hrhMhtYiKgytXLlSt99+uxYsWKC8vDxNnDhRkyZNUn5+fpv9165dq4suukirV69Wbm6ufvKTn+jSSy9VXl5emCsHAJiJvYZwPBbDMAyzi2ivs88+W2eddZaWLFnibxs2bJguu+wyLVy4sF3vcfrpp2v69Om6++6729W/rKxMLpdLbrdbSUlJHaobAGCu8pp6jbz3PUnSF3dnyxXvMLkihFog398RMzJUV1en3NxcZWdnt2jPzs7W+vXr2/UePp9P5eXl6tGjxzH71NbWqqysrMUDABDZEmMd6tM9TpK0vYjRIbQUMWGouLhYXq9XqampLdpTU1NVWFjYrvd49NFHVVlZqWnTph2zz8KFC+VyufyPjIyMk6obANA5DGnaiZpLZThKxIShJhaLpcVzwzBatbVlxYoVuvfee7Vy5UqlpKQcs9/8+fPldrv9j7179550zQAA87ETNY7FbnYB7dWzZ0/ZbLZWo0BFRUWtRouOtnLlSl1//fVatWqVLrzwwuP2dTqdcjqdJ10vAKBzyWSvIRxDxIwMxcTEKCsrSzk5OS3ac3JyNGHChGMet2LFCs2aNUsvv/yypkyZEuoyAQCdlH/jxYPliqC1QwiDiBkZkqR58+bpmmuu0ZgxYzR+/Hg9++yzys/P15w5cyQ1XOLav3+//v73v0tqCEIzZ87UX//6V40bN84/qhQXFyeXy2Xa5wAAhN9pKd1ktUilVfU6VF6rlKRYs0tCJxExI0OSNH36dD322GO6//77deaZZ2rt2rVavXq1+vfvL0kqKChosefQM888I4/Ho1tuuUXp6en+x9y5c836CAAAk8Q6bBqQnCCpYXQIaBJR+wyZgX2GAKDrmPNSrt75plC/nzJMsycOMrschFCX3GcIAICTxYoytIUwBACIGv4VZdyjDM0QhgAAUSMzrWHjxR0Hy+XzMUsEDQhDAICo0T85QTE2q6rqvNpfWm12OegkCEMAgKjhsFk1qFfjijI2X0QjwhAAIKpkph3ZfBGQCEMAgCjTtBM1K8rQhDAEAIgq3KMMRyMMAQCiStNlsl2HKlXv9ZlcDToDwhAAIKr06R6n+Bib6rw+7SmpNLscdAKEIQBAVLFaLRrsv1TG5osgDAEAolBmasPmi6wog0QYAgBEIf+KMiZRQ4QhAEAUyuSGrWiGMAQAiDpNy+u/L6lUTb3X5GpgNsIQACDq9Ep0qnu8Qz5D+q6ISdTRjjAEAIg6FouFnajhRxgCAEQl/07UhKGoRxgCAESlIWmsKEMDwhAAICpl+i+TMWco2hGGAABRaUjjxov7S6tVXlNvcjUwE2EIABCVusfHKDXJKYnRoWhHGAIARC1WlEEiDAEAoph/RRmTqKMaYQgAELWGcFsOiDAEAIhimVwmgwhDAIAoNrhxRVlxRZ2KK2pNrgZmIQwBAKJWfIxd/XrES2J0KJoRhgAAUc2/ooxJ1FGLMAQAiGqZaQ2Xyrax11DUIgwBAKIaew2BMAQAiGqZzW7YahiGydXADIQhAEBUG9Szm+xWi8prPSpw15hdDkxAGAIARLUYu1UDeyZIkrZxqSwqEYYAAFFvSBoryqIZYQgAEPWGNt2jjJGhqEQYAgBEPe5RFt0iLgwtXrxYAwcOVGxsrLKysrRu3brj9l+zZo2ysrIUGxurQYMG6emnnw5TpQCASNF0j7IdByvk9bGiLNpEVBhauXKlbr/9di1YsEB5eXmaOHGiJk2apPz8/Db77969W5MnT9bEiROVl5en3/3ud7rtttv0+uuvh7lyAEBnltEjXrEOq2o9PuUfrjK7HISZxYigTRXOPvtsnXXWWVqyZIm/bdiwYbrsssu0cOHCVv3vvPNOvf3229q6dau/bc6cOfriiy/0ySeftPk3amtrVVt75GZ9ZWVlysjIkNvtVlJSUhA/DQCgM7n0iY/01X63nr46S5eMSDO7HJyksrIyuVyudn1/R8zIUF1dnXJzc5Wdnd2iPTs7W+vXr2/zmE8++aRV/4svvlibNm1SfX19m8csXLhQLpfL/8jIyAjOBwAAdGrsRB29IiYMFRcXy+v1KjU1tUV7amqqCgsL2zymsLCwzf4ej0fFxcVtHjN//ny53W7/Y+/evcH5AACATu3IPcoIQ9HGbnYBgbJYLC2eG4bRqu1E/dtqb+J0OuV0Ok+ySgBApOHu9dErYkaGevbsKZvN1moUqKioqNXoT5O0tLQ2+9vtdiUnJ4esVgBA5Gm6R9nu4krVerwmV4NwipgwFBMTo6ysLOXk5LRoz8nJ0YQJE9o8Zvz48a36v/feexozZowcDkfIagUARJ60pFglxtrl8RnaXVxpdjkIo4gJQ5I0b948Pffcc1q6dKm2bt2qO+64Q/n5+ZozZ46khvk+M2fO9PefM2eO9uzZo3nz5mnr1q1aunSpnn/+ef3617826yMAADopi8Xi329oG5fKokpEzRmaPn26SkpKdP/996ugoEAjRozQ6tWr1b9/f0lSQUFBiz2HBg4cqNWrV+uOO+7QU089pd69e+vxxx/X5ZdfbtZHAAB0YkPSErVpzw+sKIsyEbXPkBkC2acAABDZXlz/ve55+xtdOCxVz107xuxycBK65D5DAACEGnsNRSfCEAAAjYakNuw1lH+4SlV1HpOrQbgQhgAAaJTczame3Rr2mttxsMLkahAuhCEAAJphJ+roQxgCAKAZdqKOPoQhAACa8e81xMhQ1CAMAQDQzJA0VpRFG8IQAADNDE5pmDN0sKxWpVV1JleDcCAMAQDQTGKsQ326x0mStrOiLCpE1O04AAAIh8y0RO0vrdb1L3ymPqfEKd0Vq/TucertilW6K07p3WPV2xWnNFesYh02s8vFSSIMAQBwlEkj0vTvbUUqr/Xo28JyfXuclWXJCTFK794Qkno3hqZ0V6x6N/5MTYqVw8aFmM6Me5OdAPcmA4DoVFHr0YHSah0orVaBu0YFpdU64K5RgbtaBaU1OuCuVk2974TvY7VIvRKdDWGpMTQ1D0u9u8epVzenrFZLGD5V9Ajk+5uRIQAA2tDNadeQ1ET/vkNHMwxDpVX1OtAYjgrcjWGpWWgqdNeo3mvoYFmtDpbVavPetv+W3WpRalLskbDUeBmueWjqkRAji4XAFAqEIQAAOsBiseiUhBidkhCj03u72uzj8xkqrqw9EpaOCk0F7hodLKuRx2dof2m19pdWS/qhzfdy2q0Nc5eah6XmP7vHKSnWEcJP3HURhgAACBGr1aKUxFilJMZqVEb3Nvt4vD4Vlde2DEulNUcuz7mrVVxRp1qPT9+XVOn7kqpj/r1Ep12ZaYka0cel03snaUQfl05L6cacpRNgztAJMGcIAGC2Wo9Xhe4af1gqcB8JS00/3dX1bR4bY7dqWFqiTu/j0ojeLo3ok6QhqYldfhUcc4YAAOhCnHab+icnqH9ywjH7VNV5tPdwtbYUuPX1/jJ9vd+tLQfKVF7r0Rf73Ppin9vf1261aEhqokb0SWocRXJpWHqi4mOiMxYwMnQCjAwBACKVz2co/3CVvj7QEJC+OeDWV/vdKq1qPYpktUin9urW4hLb8N5JETsPKZDvb8LQCRCGAABdiWEYOuCu0df73fpmv1tfH2gYRSoqr22z/4Dk+BaX2Eb0dumUhJgwVx04wlAQEYYAANGgqKxG3zQGo6aRpIbVba316R7nHz0a2cel0/skKSUxNswVHx9hKIgIQwCAaPVDZZ2+OVCmrxoD0jf73cdczZaS6NSIPi6N6J3UMJLUx6XerljT9kYiDAURYQgAgCPKauq1pXEEqWkkaeehCvnaSBOnxDv8E7SbLrH16xEflt22CUNBRBgCAOD4quo82lpQrm8OuBsus+0v0/aD5fK0kZASnXYNb7zE1hSQBvXqJluQAxJhKIgIQwAABK7W49X2worG+UcNE7W3FpSpztP6fm7jByVrxY3jgvr32WcIAACYymm3aWRfl0b2PXKrknqvTzsPVfj3QfrmQMOltsy0tu//Fi6EIQAAEBYOm1VD05I0NC1JV2T1lSR5fYZq6r2m1sXNSgAAgGlsVosSnOaOzRCGAABAVCMMAQCAqEYYAgAAUY0wBAAAohphCAAARDXCEAAAiGqEIQAAENUIQwAAIKoRhgAAQFSLmDD0ww8/6JprrpHL5ZLL5dI111yj0tLSY/avr6/XnXfeqZEjRyohIUG9e/fWzJkzdeDAgfAVDQAAOr2ICUMzZszQ5s2b9c477+idd97R5s2bdc011xyzf1VVlT7//HP94Q9/0Oeff6433nhD27dv19SpU8NYNQAA6OwshmEYZhdxIlu3btXw4cO1YcMGnX322ZKkDRs2aPz48fr222+VmZnZrvf57LPPNHbsWO3Zs0f9+vVr1zFlZWVyuVxyu91KSkrq8GcAAADhE8j3d0SMDH3yySdyuVz+ICRJ48aNk8vl0vr169v9Pm63WxaLRd27dz9mn9raWpWVlbV4AACArsvc28S2U2FhoVJSUlq1p6SkqLCwsF3vUVNTo7vuukszZsw4bkJcuHCh7rvvvlbthCIAACJH0/d2ey6AmRqG7r333jaDR3OfffaZJMlisbR6zTCMNtuPVl9fryuvvFI+n0+LFy8+bt/58+dr3rx5/uf79+/X8OHDlZGRccK/AwAAOpfy8nK5XK7j9jE1DN1666268sorj9tnwIAB+vLLL3Xw4MFWrx06dEipqanHPb6+vl7Tpk3T7t279cEHH5zwuqHT6ZTT6fQ/79atm/bu3avExMR2Ba9AlJWVKSMjQ3v37mU+UghxnsOD8xwenOfw4DyHT6jOtWEYKi8vV+/evU/Y19Qw1LNnT/Xs2fOE/caPHy+3261PP/1UY8eOlSRt3LhRbrdbEyZMOOZxTUFox44d+vDDD5WcnBxwjVarVX379g34uEAkJSXxf7Yw4DyHB+c5PDjP4cF5Dp9QnOsTjQg1iYgJ1MOGDdMll1yiG264QRs2bNCGDRt0ww036Kc//WmLlWRDhw7Vm2++KUnyeDy64oortGnTJi1fvlxer1eFhYUqLCxUXV2dWR8FAAB0MhERhiRp+fLlGjlypLKzs5Wdna0zzjhDL730Uos+27Ztk9vtliTt27dPb7/9tvbt26czzzxT6enp/kcgK9AAAEDXFhGrySSpR48e+sc//nHcPs1njA8YMKBdM8jN5HQ6dc8997SYo4Tg4zyHB+c5PDjP4cF5Dp/OcK4jYtNFAACAUImYy2QAAAChQBgCAABRjTAEAACiGmEIAABENcJQiC1evFgDBw5UbGyssrKytG7duuP2X7NmjbKyshQbG6tBgwbp6aefDlOlkS2Q8/zGG2/ooosuUq9evZSUlKTx48fr3XffDWO1kSvQf56bfPzxx7Lb7TrzzDNDW2AXEeh5rq2t1YIFC9S/f385nU6deuqpWrp0aZiqjVyBnufly5dr1KhRio+PV3p6uq677jqVlJSEqdrItHbtWl166aXq3bu3LBaL3nrrrRMeY8r3oIGQeeWVVwyHw2H87W9/M7Zs2WLMnTvXSEhIMPbs2dNm/127dhnx8fHG3LlzjS1bthh/+9vfDIfDYbz22mthrjyyBHqe586dazz88MPGp59+amzfvt2YP3++4XA4jM8//zzMlUeWQM9zk9LSUmPQoEFGdna2MWrUqPAUG8E6cp6nTp1qnH322UZOTo6xe/duY+PGjcbHH38cxqojT6Dned26dYbVajX++te/Grt27TLWrVtnnH766cZll10W5sojy+rVq40FCxYYr7/+uiHJePPNN4/b36zvQcJQCI0dO9aYM2dOi7ahQ4cad911V5v9f/vb3xpDhw5t0XbTTTcZ48aNC1mNXUGg57ktw4cPN+67775gl9aldPQ8T58+3fj9739v3HPPPYShdgj0PP/zn/80XC6XUVJSEo7yuoxAz/MjjzxiDBo0qEXb448/bvTt2zdkNXY17QlDZn0PcpksROrq6pSbm6vs7OwW7dnZ2cfcAfuTTz5p1f/iiy/Wpk2bVF9fH7JaI1lHzvPRfD6fysvL1aNHj1CU2CV09DwvW7ZMO3fu1D333BPqEruEjpznt99+W2PGjNGf/vQn9enTR0OGDNGvf/1rVVdXh6PkiNSR8zxhwgTt27dPq1evlmEYOnjwoF577TVNmTIlHCVHDbO+ByNmB+pIU1xcLK/Xq9TU1BbtqampKiwsbPOYwsLCNvt7PB4VFxcrPT09ZPVGqo6c56M9+uijqqys1LRp00JRYpfQkfO8Y8cO3XXXXVq3bp3sdv5V0x4dOc+7du3SRx99pNjYWL355psqLi7WzTffrMOHDzNv6Bg6cp4nTJig5cuXa/r06aqpqZHH49HUqVP1xBNPhKPkqGHW9yAjQyFmsVhaPDcMo1Xbifq31Y6WAj3PTVasWKF7771XK1euVEpKSqjK6zLae569Xq9mzJih++67T0OGDAlXeV1GIP88+3w+WSwWLV++XGPHjtXkyZO1aNEivfDCC4wOnUAg53nLli267bbbdPfddys3N1fvvPOOdu/erTlz5oSj1Khixvcg/7kWIj179pTNZmv1XxlFRUWtUm+TtLS0Nvvb7XYlJyeHrNZI1pHz3GTlypW6/vrrtWrVKl144YWhLDPiBXqey8vLtWnTJuXl5enWW2+V1PClbRiG7Ha73nvvPZ1//vlhqT2SdOSf5/T0dPXp00cul8vfNmzYMBmGoX379mnw4MEhrTkSdeQ8L1y4UOecc45+85vfSJLOOOMMJSQkaOLEiXrwwQcZuQ8Ss74HGRkKkZiYGGVlZSknJ6dFe05OjiZMmNDmMePHj2/V/7333tOYMWPkcDhCVmsk68h5lhpGhGbNmqWXX36Za/7tEOh5TkpK0ldffaXNmzf7H3PmzFFmZqY2b96ss88+O1ylR5SO/PN8zjnn6MCBA6qoqPC3bd++XVarVX379g1pvZGqI+e5qqpKVmvLr0ybzSZJnf6m4JHEtO/BkE7PjnJNSzeff/55Y8uWLcbtt99uJCQkGN9//71hGIZx1113Gddcc42/f9OSwjvuuMPYsmWL8fzzz7O0vh0CPc8vv/yyYbfbjaeeesooKCjwP0pLS836CBEh0PN8NFaTtU+g57m8vNzo27evccUVVxjffPONsWbNGmPw4MHG7NmzzfoIESHQ87xs2TLDbrcbixcvNnbu3Gl89NFHxpgxY4yxY8ea9REiQnl5uZGXl2fk5eUZkoxFixYZeXl5/i0MOsv3IGEoxJ566imjf//+RkxMjHHWWWcZa9as8b927bXXGueee26L/v/+97+N0aNHGzExMcaAAQOMJUuWhLniyBTIeT733HMNSa0e1157bfgLjzCB/vPcHGGo/QI9z1u3bjUuvPBCIy4uzujbt68xb948o6qqKsxVR55Az/Pjjz9uDB8+3IiLizPS09ON//7v/zb27dsX5qojy4cffnjcf992lu9Bi2EwvgcAAKIXc4YAAEBUIwwBAICoRhgCAABRjTAEAACiGmEIAABENcIQAACIaoQhAAAQ1QhDAAAgqhGGAHRa5513nm6//fZO+TcGDBigxx57LOj1AAg/whAAAIhqhCEAABDVCEMAIsI//vEPjRkzRomJiUpLS9OMGTNUVFTkf/3f//63LBaL3n33XY0ePVpxcXE6//zzVVRUpH/+858aNmyYkpKSdNVVV6mqqqrFe3s8Ht16663q3r27kpOT9fvf/17Nb9tYVFSkSy+9VHFxcRo4cKCWL1/eqr5FixZp5MiRSkhIUEZGhm6++WZVVFSE7oQACBrCEICIUFdXpwceeEBffPGF3nrrLe3evVuzZs1q1e/ee+/Vk08+qfXr12vv3r2aNm2aHnvsMb388sv6v//7P+Xk5OiJJ55occyLL74ou92ujRs36vHHH9df/vIXPffcc/7XZ82ape+//14ffPCBXnvtNS1evLhFEJMkq9Wqxx9/XF9//bVefPFFffDBB/rtb38bknMBIMhO/sb3ABAa5557rjF37tw2X/v0008NSUZ5eblhGIbx4YcfGpKM999/399n4cKFhiRj586d/rabbrrJuPjii1v8jWHDhhk+n8/fdueddxrDhg0zDMMwtm3bZkgyNmzY4H9969athiTjL3/5yzFrf/XVV43k5OSAPi8AczAyBCAi5OXl6Wc/+5n69++vxMREnXfeeZKk/Pz8Fv3OOOMM/++pqamKj4/XoEGDWrQdPaozbtw4WSwW//Px48drx44d8nq92rp1q+x2u8aMGeN/fejQoerevXuL9/jwww910UUXqU+fPkpMTNTMmTNVUlKiysrKk/3oAEKMMASg06usrFR2dra6deumf/zjH/rss8/05ptvSmq4fNacw+Hw/26xWFo8b2rz+Xzt/ttG49yh5mHpaHv27NHkyZM1YsQIvf7668rNzdVTTz0lSaqvr2/33wJgDrvZBQDAiXz77bcqLi7WQw89pIyMDEnSpk2bgvb+GzZsaPV88ODBstlsGjZsmDwejzZt2qSxY8dKkrZt26bS0lJ//02bNsnj8ejRRx+V1drw35ivvvpq0OoDEFqMDAHo9Pr166eYmBg98cQT2rVrl95++2098MADQXv/vXv3at68edq2bZtWrFihJ554QnPnzpUkZWZm6pJLLtENN9ygjRs3Kjc3V7Nnz1ZcXJz/+FNPPVUej8df30svvaSnn346aPUBCC3CEIBOr1evXnrhhRe0atUqDR8+XA899JD+/Oc/B+39Z86cqerqao0dO1a33HKLfvnLX+rGG2/0v75s2TJlZGTo3HPP1X/913/pxhtvVEpKiv/1M888U4sWLdLDDz+sESNGaPny5Vq4cGHQ6gMQWhbDaLaZBgAAQJRhZAgAAEQ1whAAAIhqhCEAABDVCEMAACCqEYYAAEBUIwwBAICoRhgCAABRjTAEAACiGmEIAABENcIQAACIaoQhAAAQ1f4/37nsPFHZv6oAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "E_diff = np.diff(E_lst, axis=-1).flatten() / len(fcc)\n", "plt.xlabel(\"lambda\")\n", @@ -214,10 +225,18 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 12, "id": "b70a0bb0", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "The free energy difference between bcc and fcc at 300 K is 0.957321594716624 eV\n" + ] + } + ], "source": [ "print(\"The free energy difference between bcc and fcc at 300 K is\", E_diff.sum() * np.diff(ti_lambda).mean(), \"eV\")" ] @@ -240,7 +259,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 13, "id": "0dee5e1a", "metadata": {}, "outputs": [], @@ -253,7 +272,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 14, "id": "aa94c95b", "metadata": {}, "outputs": [], @@ -271,10 +290,31 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 15, "id": "7abfdf91", "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "[]" + ] + }, + "execution_count": 15, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAGwCAYAAACpYG+ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAABTPUlEQVR4nO3dd3hT9f4H8HeSpune0D0pe9MyWmZFAQVUuCiKIiDg5SIXBFFAvCxFrgLOn4IiAioKKuBWyoayCqXMUkZpaUsXbaG7SZOc3x+ludYympLkZLxfz5PnoScnyaeH0vPmOyWCIAggIiIisgFSsQsgIiIiMhUGHyIiIrIZDD5ERERkMxh8iIiIyGYw+BAREZHNYPAhIiIim8HgQ0RERDbDTuwCzI1Wq0VOTg5cXV0hkUjELoeIiIgaQRAElJWVISAgAFLpndt1GHz+JicnB8HBwWKXQURERE2QlZWFoKCgOz7P4PM3rq6uAGovnJubm8jVEBERUWOUlpYiODhYdx+/Ewafv6nr3nJzc2PwISIisjD3GqbCwc1ERERkMxh8iIiIyGZYZfD55JNPEB4eDgcHB0RFReHAgQNil0RERERmwOqCz+bNm/HSSy9h/vz5SE5ORt++ffHwww8jMzNT7NKIiIhIZBJBEASxizCknj17olu3bli1apXuWNu2bfH4449j2bJl93x9aWkp3N3dUVJSwsHNREREFqKx92+ravFRqVRISkrCoEGD6h0fNGgQDh06dNvXKJVKlJaW1nsQERGRdbKq4FNYWAiNRgNfX996x319fZGXl3fb1yxbtgzu7u66BxcvJCIisl5WFXzq/H0OvyAId5zXP2/ePJSUlOgeWVlZpiiRiIiIRGBVCxj6+PhAJpM1aN0pKCho0ApUR6FQQKFQmKI8IiIiEplVtfjY29sjKioKO3bsqHd8x44diI2NFakqIiIiMhdW1eIDALNmzcLYsWMRHR2NmJgYfPbZZ8jMzMSUKVPELo2IiIhEZnXBZ/To0SgqKsKSJUuQm5uLDh064Pfff0doaKjYpREREZHIrG4dn/vFdXyIiIiMo7pGg9PZJegR7mXw927s/dvqWnyIiIjIfKg1WhxKK8JPJ3MQfy4PlTUaJL42EN4u4kwsYvAhIiIig0vNK8WmxCz8ejoHheUq3fEAdwdkFFUy+BAREZFlq9Fosf1cHr48fBWJ6cW6455Ocgzt5I9HOwciOtQTUunt19YzBQYfIiIiui9l1TXYcCgDXx25ivxSJQBAJpVgcHtfPBEdjD6RPpDLzGMFHQYfIiIiapK6wLPmQDpKqmoAAD4uCozpEYwxPUPh5+4gcoUNMfgQERGRXipVaqw7mIE1B67gZmVt4GnRzBnTHojE0I4BsLczj9ad22HwISIiokYRBAF/ns3DG7+mIKekGgAQ0cwZMwa2xLBOAZCJOHansRh8iIiI6J7SCyuw8Odz2H/xOgAg0MMRswe3wqOdAy0i8NRh8CEiIqI7Uqm1+Gj3JXy67wpUGi3sZVJM6R+Bfw2IhKO9TOzy9MbgQ0RERLeVdr0cL206iTPXSgAA/Vs1w+JH2yPMx1nkypqOwYeIiIjqEQQBm49lYfEvKaiq0cDDSY63RnTEwx38IJFYTrfW7TD4EBERkc7NShXmbT2DP87mAQBiW3jj3Se7mOXU9KZg8CEiIiIAtV1bz68/hqtFlbCTSjB7cGu80DdC1JWWDY3Bh4iIiHAorRBTvkpCabUagR6OWPVsN3QK8hC7LINj8CEiIrJxm49lYv62s1BrBXQL8cBnz0XDR6RNRI2NwYeIiMhGabUC3t6eik/3XQEAPNo5AO+M6gQHueVNU28sBh8iIiIbpNUKmP/jWXybmAkAmDGwJV56sKXFz9q6FwYfIiIiG/PX0COVAO+M6oxRUUFil2USDD5EREQ2pDb0nMG3iVmQSoCVT3bGiK62EXoABh8iIiKb8ffQ8+6TXfB410CxyzIp8903noiIiAxGEAQs+PmsTYcegMGHiIjIJny2/wq+PpIJiQ2HHoDBh4iIyOr9cSYXy/5IBQAsGNbOZkMPwOBDRERk1U5m3cRLm08CAMbFhGJC73BxCxIZgw8REZGVyiquxKQNx6BUaxHXuhn+M6yd2CWJjsGHiIjICpVW1+D59cdQWK5CW383fDSmG+xkvO3zChAREVkZQRAwb+sZXCooh6+bAl+Mj4aLgivYAAw+REREVuebxEz8djoXdlIJVj8bBX93R7FLMhsMPkRERFbkfG4plvySAgB4dUhrdA3xFLki88LgQ0REZCUqVWpM++aEbjDzpD4RYpdkdhh8iIiIrMSCn84h7XoFfN0UWPlkF0il1r3TelMw+BAREVmBrSey8UNSNqQS4IOnusLL2V7skswSgw8REZGFyy2pwoKfzgEAZgxshV4R3iJXZL4YfIiIiCyYIAh4fdtZlCvV6BbigWkPRIpdkllj8CEiIrJgv57Oxa7UAshlErz9j06QcVzPXTH4EBERWagbFSos+rm2i+vFuEi09HUVuSLzx+BDRERkod74LQVFFSq08nXB1AHs4moMBh8iIiILtO/idWw9cQ0SCfDff3SCvR1v6Y3Bq0RERGRhKpRqvLb1DABgXEwYunF15kZj8CEiIrIwH+2+jGs3qxDo4YhXBrcWuxyLwuBDRERkQbKKK/FFQjoAYNGj7eHMXdf1wuBDRERkQf77ZypUGi16R3rjwbbNxS7H4jD4EBERWYikq8X47XQuJBJg/iPtIJFwzR59MfgQERFZAK1WwJJfzwMARkcHo12Am8gVWSYGHyIiIgvwy+kcnMq6CWd7GWYNaiV2ORbLqoJPWFgYJBJJvcfcuXPFLouIiOi+VKk0ePuPVADA1LhINHd1ELkiy2V1Q8GXLFmCyZMn6752cXERsRoiIqL7tzbhCnJKqhHo4YiJfcLFLseiWV3wcXV1hZ+fn9hlEBERGURxhQqr9qYBAF4d0hoOcpnIFVk2q+rqAoC3334b3t7e6NKlC5YuXQqVSnXX85VKJUpLS+s9iIiIzMWaA1dQodKgQ6AbHu0cIHY5Fs+qWnxmzJiBbt26wdPTE4mJiZg3bx7S09Px+eef3/E1y5Ytw+LFi01YJRERUeMUV6iw4VAGAOClga04fd0AJIIgCGIXcTeLFi26ZzA5duwYoqOjGxzfsmULRo0ahcLCQnh7e9/2tUqlEkqlUvd1aWkpgoODUVJSAjc3ThUkIiLxvP1nKlbtTUOHQDf8Mq0Pg89dlJaWwt3d/Z73b7Nv8Zk2bRqeeuqpu54TFhZ22+O9evUCAFy+fPmOwUehUEChUNxXjURERIbG1h7jMPvg4+PjAx8fnya9Njk5GQDg7+9vyJKIiIiMbs2BK6i8NbZnILemMBizDz6NdfjwYRw5cgRxcXFwd3fHsWPHMHPmTDz66KMICQkRuzwiIqJGY2uP8VhN8FEoFNi8eTMWL14MpVKJ0NBQTJ48Ga+++qrYpREREemFrT3GYzXBp1u3bjhy5IjYZRAREd0XtvYYl9Wt40NERGTJ1iawtceYGHyIiIjMRKVKja+PZAIApsW1ZGuPETD4EBERmYkfkrJRUlWDMG8nPNTOV+xyrBKDDxERkRnQaAWsTUgHADzfJxwyKVt7jIHBh4iIyAzsPJ+Pq0WVcHeUY1RUkNjlWC0GHyIiIjPw+YErAIBneobAyd5qJl2bHQYfIiIikZ3MuoljGTcgl0kwLjZM7HKsGoMPERGRyOpaex7tHAhfNweRq7FuDD5EREQiyr5RiT/O5gEAJvYJF7ka68fgQ0REJKL1BzOg0QroE+mDdgFuYpdj9Rh8iIiIRFJWXYNNx7IAABP7srXHFBh8iIiIRPLjyRyUK9Vo0cwZA1o1E7scm8DgQ0REJAJBEPDN0drtKZ7pGcrtKUyEwYeIiEgEyVk3cT63FAo7Kf7RjQsWmgqDDxERkQjqWnuGdQqAu5Nc5GpsB4MPERGRiZVU1uCXUzkAgDE9Q0SuxrYw+BAREZnY1uRsKNVatPFzRbcQD7HLsSkMPkRERCYkCAI26gY1h3BQs4kx+BAREZnQsYwbuFxQDid7GR7vGih2OTaHwYeIiMiEvjl6FQDwaOcAuDpwULOpMfgQERGZSHGFCr+fqd2X65meoSJXY5sYfIiIiExkS1I2VBotOga6o2OQu9jl2CQGHyIiIhMQBAGbj9fuy8Up7OJh8CEiIjKBM9dKcLmgHAo7KYZ18he7HJvF4ENERGQCW5KyAQCD2/txULOIGHyIiIiMTKXW4udbKzX/I4r7comJwYeIiMjI9lwowI3KGjR3VaBPpI/Y5dg0Bh8iIiIj23qitptrRNdAyKRcqVlMDD5ERERGdKNChd2pBQCAkd3YzSU2Bh8iIiIj+uV0Dmo0AjoEuqG1n6vY5dg8Bh8iIiIjqpvNNbIrW3vMAYMPERGRkVwuKMOp7BLYSSV4tEuA2OUQGHyIiIiMZsuJawCAAa2bwcdFIXI1BDD4EBERGYVGK2DbreDDQc3mg8GHiIjICI5eKUJeaTXcHOwwsG1zscuhWxh8iIiIjOCX07UrNT/S0R8KO5nI1VAdBh8iIiIDq9Fo8cfZPADA8M4c1GxOGHyIiIgMLOFyIW5W1sDHxR49w73ELof+gsGHiIjIwH49lQsAeLiDP+xkvNWaE/5tEBERGZBSrUF8Cru5zBWDDxERkQHtv1iIsmo1/NwcEB3qKXY59DcMPkRERAb0y6n/zeaScid2s8PgQ0REZCBVKg12ns8HAAzr7C9yNXQ7DD5EREQGsudCASpVGgR6OKJrsIfY5dBtMPgQEREZyK+3Fi0c1tkfEgm7ucwRgw8REZEBlCvV2HW+AAAwvBNnc5kriwk+S5cuRWxsLJycnODh4XHbczIzMzF8+HA4OzvDx8cH06dPh0qlMm2hRERkk3adz4dSrUWYtxPaB7iJXQ7dgZ3YBTSWSqXCE088gZiYGKxdu7bB8xqNBkOHDkWzZs2QkJCAoqIijBs3DoIg4KOPPhKhYiIisiW/3Fq0cFinAHZzmTGLCT6LFy8GAKxfv/62z8fHxyMlJQVZWVkICKhtYly5ciXGjx+PpUuXws3t9ulbqVRCqVTqvi4tLTVs4UREZPXKlWrsv3QdADC0E2dzmTOL6eq6l8OHD6NDhw660AMAgwcPhlKpRFJS0h1ft2zZMri7u+sewcHBpiiXiIisyL4L16G61c3Vxs9V7HLoLqwm+OTl5cHX17feMU9PT9jb2yMvL++Or5s3bx5KSkp0j6ysLGOXSkREVubPc7X3mcHt/djNZeZEDT6LFi2CRCK56+P48eONfr/b/bAJgnDXH0KFQgE3N7d6DyIiosZSqjXYk1o7m2twBz+Rq6F7EXWMz7Rp0/DUU0/d9ZywsLBGvZefnx+OHj1a79iNGzdQU1PToCWIiIjIUA6lFaFcqUZzVwW6BHmIXQ7dg6jBx8fHBz4+PgZ5r5iYGCxduhS5ubnw968dWBYfHw+FQoGoqCiDfAYREdHfbT/7v24u7s1l/ixmVldmZiaKi4uRmZkJjUaDkydPAgAiIyPh4uKCQYMGoV27dhg7diyWL1+O4uJizJ49G5MnT2b3FRERGYVGK2BHSu3eXIPbs5vLElhM8FmwYAE2bNig+7pr164AgD179mDAgAGQyWT47bffMHXqVPTu3RuOjo4YM2YMVqxYIVbJRERk5Y5nFKOoQgV3Rzl6RniJXQ41gsUEn/Xr199xDZ86ISEh+PXXX01TEBER2bzt52pbewa2bQ65zGomSls1/i0RERE1gSAI2H5rGvsQdnNZDAYfIiKiJjiXU4prN6vgKJehX6tmYpdDjcTgQ0RE1AR/3prNNaB1MzjIZSJXQ43F4ENERNQE2/+yWjNZDgYfIiIiPaVdL8elgnLIZRLEtWkudjmkh0YHnx9//BEajcaYtRAREVmEurV7Ylr4wN1RLnI1pI9GB59Ro0YhMDAQc+bMQWpqqjFrIiIiMmu7ztcGn4fasrXH0jQ6+GRmZuLf//43tm3bhvbt26NPnz5Yt24dKioqjFkfERGRWblRoULS1RsAwG4uC9To4BMQEID58+fj4sWL2L17N1q0aIHp06fD398fkyZNwuHDh41ZJxERkVnYc6EAWgFo4+eKIE8nscshPTVpcHP//v2xYcMG5Obm4t1338X58+fRp08ftG/f3tD1ERERmZVd5wsAAA+29RW5EmqK+9qywsXFBXFxccjIyEBqaiouXrxoqLqIiIjMjkqtxb6L1wHUblNBlqdJLT6VlZXYsGED+vfvj1atWmHz5s2YNWsWMjIyDFweERGR+UhML0a5Ug0fFwU6B3mIXQ41gV4tPgcPHsQXX3yB77//Hmq1GiNHjsTOnTsRFxdnrPqIiIjMxs5bs7keaNMMUqlE5GqoKRodfFq1aoW0tDR07doVb7/9NsaMGQN3d3dj1kZERGQ2BEHArtS63dg5vsdSNTr4DBkyBBMnTkTnzp2NWQ8REZFZulRQjqziKtjbSdG3pY/Y5VATNTr4fPjhh7o/q9Vq7N27F2lpaRgzZgxcXV2Rk5MDNzc3uLi4GKVQIiIiMdV1c8W28IaT/X3NDSIR6f03d/XqVQwZMgSZmZlQKpV46KGH4OrqinfeeQfV1dVYvXq1MeokIiISVd00dnZzWTa9Z3XNmDED0dHRuHHjBhwdHXXHR4wYgV27dhm0OCIiInNQVK7Eicza1ZoHcrVmi6Z3i09CQgIOHjwIe3v7esdDQ0Nx7do1gxVGRERkLvZcuA5BANr5uyHAw/HeLyCzpXeLj1arve0u7dnZ2XB1dTVIUUREROZkt242F1t7LJ3eweehhx7C+++/r/taIpGgvLwcCxcuxCOPPGLI2oiIiESnUmux/2IhAI7vsQZ6d3W99957iIuLQ7t27VBdXY0xY8bg0qVL8PHxwbfffmuMGomIiERz/Grdas326BTI9essnd7BJyAgACdPnsSmTZuQlJQErVaLiRMn4plnnqk32JmIiMga7LtQuzdXv1ZcrdkaNGkhAkdHR0yYMAETJkwwdD1ERERmZc+F2mnsca05vscaNGmTUiIiIltw7WYVLuaXQyoBV2u2Egw+REREd7D3VmtPtxBPeDjZ3+NssgQMPkRERHewJ7V2fM+A1s1EroQMhcGHiIjoNpRqDQ6l1U5jH8DxPVZD7+Azfvx47N+/3xi1EBERmY1j6TdQqdKguasC7QPcxC6HDETv4FNWVoZBgwahZcuWeOutt7hNBRERWaW62Vz9WzWDRMJp7NZC7+CzZcsWXLt2DdOmTcP333+PsLAwPPzww/jhhx9QU1NjjBqJiIhMrm5gcxw3JbUqTRrj4+3tjRkzZiA5ORmJiYmIjIzE2LFjERAQgJkzZ+LSpUuGrpOIiMhkMosqkXa9AjKpBH04jd2q3Nfg5tzcXMTHxyM+Ph4ymQyPPPIIzp07h3bt2uG9994zVI1EREQmtfdibWtPVKgn3BzkIldDhqR38KmpqcGWLVswbNgwhIaG4vvvv8fMmTORm5uLDRs2ID4+Hl999RWWLFlijHqJiIiMbu+tbSq4WrP10XvLCn9/f2i1Wjz99NNITExEly5dGpwzePBgeHh4GKA8IiIi06qu+d809rg2XL/H2jRpd/YnnngCDg4OdzzH09MT6enp91UYERGRGI6mF6O6Rgs/Nwe09nUVuxwyML2Dz9ixY41RBxERkVmom801oDWnsVsjvYPPyJEjb3tcIpHAwcEBkZGRGDNmDFq3bn3fxREREZna/ou143v6t2I3lzXSe3Czm5sbdu/ejRMnTuiScHJyMnbv3g21Wo3Nmzejc+fOOHjwoMGLJSIiMqbsG/+bxh4byWns1kjvFh8/Pz+MGTMG//d//weptDY3abVazJgxA66urti0aROmTJmCOXPmICEhweAFExERGcv+i7WDmrsGe8DdkdPYrZHeLT5r167FSy+9pAs9ACCVSvHvf/8bn332GSQSCaZNm4azZ88atFAiIiJjq+vm6sduLquld/BRq9VITU1tcDw1NRUajQYA4ODgwAFhRERkUWo0Why8XNviw+BjvZo0q2vixIl47bXX0L17d0gkEiQmJuKtt97Cc889BwDYt28f2rdvb/BiiYiIjOVk1k2UKdXwcJKjY6C72OWQkTRpHR9fX1+88847yM/PBwD4+vpi5syZmDNnDgBg0KBBGDJkiGErJSIiMqK6bq6+LZtBJmWvhbXSK/io1Wps3LgREydOxPz581FaWgqgdqbXX4WEhBiuQiIiIhPQje/hpqRWTa8xPnZ2dvjXv/4FpVIJoDbw/D30GMvSpUsRGxsLJyenO26HIZFIGjxWr15tkvqIiMhyFVeocPpaCQCO77F2end19ezZE8nJyQgNDTVGPXekUqnwxBNPICYmBmvXrr3jeevWravXzebuzn5aIiK6uwOXrkMQgDZ+rvB1u/OWTGT59A4+U6dOxcsvv4zs7GxERUXB2dm53vOdOnUyWHF/tXjxYgDA+vXr73qeh4cH/Pz8jFIDERFZp7r1e7has/XTO/iMHj0aADB9+nTdMYlEAkEQIJFIdFPaxTJt2jRMmjQJ4eHhmDhxIl544YV6aw79nVKp1HXdAdCNWyIiItsgCAIOXOL6PbZC7+Bjzruuv/HGGxg4cCAcHR2xa9cuvPzyyygsLMTrr79+x9csW7ZM15pERES2JzWvDAVlSjjKZYgO8xS7HDIyvYOPIcf2LFq06J6h49ixY4iOjm7U+/014HTp0gUAsGTJkrsGn3nz5mHWrFm6r0tLSxEcHNyozyMiIsu379Zsrl4RXlDYyUSuhoxN7+ADAF999RVWr16N9PR0HD58GKGhoXj//fcRHh6Oxx57rNHvM23aNDz11FN3PScsLKwpJQIAevXqhdLSUuTn58PX1/e25ygUCigUiiZ/BhERWTZuU2Fb9A4+q1atwoIFC/DSSy9h6dKlujE9Hh4eeP/99/UKPj4+PvDxMd56CcnJyXBwcLjj9HciIrJtlSo1jmfcAMDgYyv0Dj4fffQR1qxZg8cffxz//e9/dcejo6Mxe/Zsgxb3V5mZmSguLkZmZiY0Gg1OnjwJAIiMjISLiwt++eUX5OXlISYmBo6OjtizZw/mz5+PF154gS06RER0W0fTi6HSaBHk6YgIH+d7v4AsXpMGN3ft2rXBcYVCgYqKCoMUdTsLFizAhg0bdF/X1bBnzx4MGDAAcrkcn3zyCWbNmgWtVouIiAgsWbIEL774otFqIiIiy3bg1jT2vi19uLm2jdA7+ISHh+PkyZMNBjn/8ccfaNeuncEK+7v169ffdQ2fIUOGcH8wIiLSS9009r4t2c1lK/QOPq+88gpefPFFVFdXQxAEJCYm4ttvv8WyZcvw+eefG6NGIiIig8stqcKlgnJIJUBsC2+xyyET0Tv4TJgwAWq1Gq+++ioqKysxZswYBAYG4oMPPrjnDC0iIiJzceBSbTdXpyAPeDjZi1wNmUqTprNPnjwZkydPRmFhIbRaLZo3b27ouoiIiIyqLvhwN3bb0qTgU8eYU9GJiIiMRasVkFA3vofT2G3KnTexuoP8/HyMHTsWAQEBsLOzg0wmq/cgIiIyd+dySnGjsgYuCjt0CfYQuxwyIb1bfMaPH4/MzEz85z//gb+/P6f/ERGRxdl/q7UnpoU35DK92wDIgukdfBISEnDgwAHdXlhERESWRrcbO8f32By9Y25wcDAEQTBGLUREREZXoVQj6WrtNhV9uH6PzdE7+Lz//vuYO3cuMjIyjFAOERGRcR1NL0KNRkCQpyPCvJ3ELodMTO+urtGjR6OyshItWrSAk5MT5HJ5veeLi4sNVhwREZGh1U1j79uyGcep2iC9g8/7779vhDKIiIhMg+v32Da9g8+4ceOMUQcREZHR5dyswmXdNhUMPraoSXP40tLS8Prrr+Ppp59GQUEBAODPP//EuXPnDFocERGRISXcau3pHOwBdyf5Pc4ma6R38Nm3bx86duyIo0ePYuvWrSgvLwcAnD59GgsXLjR4gURERIZSt35P30i29tgqvYPP3Llz8eabb2LHjh2wt//fpm5xcXE4fPiwQYsjIiIyFK1WwMHLtwY2c5sKm6V38Dlz5gxGjBjR4HizZs1QVFRkkKKIiIgMjdtUENCE4OPh4YHc3NwGx5OTkxEYGGiQooiIiAytrpurVwS3qbBlev/NjxkzBnPmzEFeXh4kEgm0Wi0OHjyI2bNn47nnnjNGjURERPetbmBzv1Yc32PL9A4+S5cuRUhICAIDA1FeXo527dqhX79+iI2Nxeuvv26MGomIiO5LpUqN41drF9jtw4HNNk3vdXzkcjk2btyIN954AydOnIBWq0XXrl3RsmVLY9RHRER0346mF6NGIyDQwxHhPs5il0Mi0jv41ImIiEBERIQhayEiIjKKBN02FT7cpsLGcXQXERFZvQO3Bjb34TYVNo/Bh4iIrFp+aTUu5pdDIgF6c5sKm8fgQ0REVq1uU9KOge7wdLa/x9lk7Rh8iIjIqiXUbVPBbi5CE4PPgQMH8OyzzyImJgbXrl0DAHz11VdISEgwaHFERET3Q6sVkHBrm4o+kdymgpoQfLZs2YLBgwfD0dERycnJUCqVAICysjK89dZbBi+QiIioqVLzylBYroKTvQzdQj3ELofMgN7B580338Tq1auxZs0ayOVy3fHY2FicOHHCoMURERHdj7rZXD3DvaCwk4lcDZkDvYPPhQsX0K9fvwbH3dzccPPmTUPUREREZBB13Vx9W7Kbi2rpHXz8/f1x+fLlBscTEhK4oCEREZmN6hoNjqbXblPBgc1UR+/g889//hMzZszA0aNHIZFIkJOTg40bN2L27NmYOnWqMWokIiLS27GMYqjUWvi5OSCyuYvY5ZCZ0HvLildffRUlJSWIi4tDdXU1+vXrB4VCgdmzZ2PatGnGqJGIiEhvB7hNBd1Gk/bqWrp0KebPn4+UlBRotVq0a9cOLi5M00REZD72X7y1fk8rju+h/2nyAoY5OTkoKipCx44d4eLiAkEQDFkXERFRkxWUViM1rwwSCdAnkuN76H/0Dj5FRUUYOHAgWrVqhUceeQS5ubkAgEmTJuHll182eIFERET6qpvN1SHAHV7cpoL+Qu/gM3PmTMjlcmRmZsLJyUl3fPTo0fjzzz8NWhwREVFT/HV8D9Ff6T3GJz4+Htu3b0dQUFC94y1btsTVq1cNVhgREVFTaLXCX4IPx/dQfXq3+FRUVNRr6alTWFgIhUJhkKKIiIiaqnabCiW3qaDb0jv49OvXD19++aXua4lEAq1Wi+XLlyMuLs6gxREREemrbpuKXhHe3KaCGtC7q2v58uUYMGAAjh8/DpVKhVdffRXnzp1DcXExDh48aIwaiYiIGo3je+hu9G7xadeuHU6fPo0ePXrgoYceQkVFBUaOHInk5GS0aNHCGDUSERE1SpVKg8SMum0qOL6HGtKrxaempgaDBg3Cp59+isWLFxurJiIioiZJvLVNRYC7A1o0cxa7HDJDerX4yOVynD17lkt/ExGRWTpQt1pzy2a8V9Ft6d3V9dxzz2Ht2rXGqIWIiOi+6Mb3tOL4Hro9vQc3q1QqfP7559ixYweio6Ph7Fy/KfHdd981WHFERESNlV9ajQv5tdtU9G7B4EO3p3eLz9mzZ9GtWze4ubnh4sWLSE5O1j1OnjxphBKBjIwMTJw4EeHh4XB0dESLFi2wcOFCqFSqeudlZmZi+PDhcHZ2ho+PD6ZPn97gHCIisk51rT2dAt3hyW0q6A4a3eJz5coVhIeHY8+ePcas57ZSU1Oh1Wrx6aefIjIyEmfPnsXkyZNRUVGBFStWAAA0Gg2GDh2KZs2aISEhAUVFRRg3bhwEQcBHH31k8pqJiMi06tbv4WwuuhuJ0Mht1WUyGXJzc9G8eXMAtXtzffjhh/D19TVqgXeyfPlyrFq1CleuXAEA/PHHHxg2bBiysrIQEBAAANi0aRPGjx+PgoICuLm5Nep9S0tL4e7ujpKSkka/hoiIxKXVCoheuhPFFSpsfqEXekZ4i10SmVhj79+N7ur6ez76/fffUVFR0fQK71NJSQm8vLx0Xx8+fBgdOnTQhR4AGDx4MJRKJZKSku74PkqlEqWlpfUeRERkWc7mlKC4QgUXhR26hXqKXQ6ZMb3H+JiDtLQ0fPTRR5gyZYruWF5eXoPWJ09PT9jb2yMvL++O77Vs2TK4u7vrHsHBwUarm4iIjGPfhdpurtgW3pDLLPLWRibS6J8OiUTSYE2E+10jYdGiRbr3vdPj+PHj9V6Tk5ODIUOG4IknnsCkSZPuWY8gCHetc968eSgpKdE9srKy7ut7IiIi09t/a3xP/9Yc30N31+jBzYIgYPz48bod2KurqzFlypQG09m3bt3a6A+fNm0annrqqbueExYWpvtzTk4O4uLiEBMTg88++6zeeX5+fjh69Gi9Yzdu3EBNTc1dxyEpFAruKk9EZMFKqmpwIvMmAKAfBzbTPTQ6+IwbN67e188+++x9f7iPjw98fBq31sK1a9cQFxeHqKgorFu3DlJp/caqmJgYLF26FLm5ufD39wcAxMfHQ6FQICoq6r5rJSIi83TociE0WgERzZwR7OUkdjlk5hodfNatW2fMOu4qJycHAwYMQEhICFasWIHr16/rnvPz8wMADBo0CO3atcPYsWOxfPlyFBcXY/bs2Zg8eTJnZxERWTFdN1crtvbQvem9crMY4uPjcfnyZVy+fBlBQUH1nqubbSaTyfDbb79h6tSp6N27NxwdHTFmzBjdOj9ERGR9BEHQDWxm8KHGaPQ6PraC6/gQEVmOS/lleOi9/bC3k+LUgkFwtJeJXRKJxODr+BAREZmbfbd2Y+8Z7sXQQ43C4ENERBarLviwm4sai8GHiIgsUpVKg6PpxQAYfKjxGHyIiMgiHU0vgkqtRYC7AyKbu4hdDlkIBh8iIrJIum6u1s3ueycBsh0MPkREZJHqgg9XayZ9MPgQEZHFySquxJXrFZBJJYiNbNwOAEQAgw8REVmgvbdae7oGe8DdUS5yNWRJGHyIiMji7E0tAADEtWkuciVkaRh8iIjIolTXaHAwrRAA8ACDD+mJwYeIiCzKkStFqK7Rws/NAW38XMUuhywMgw8REVmUvbc2JY1rw2nspD8GHyIishiCIGB33fie1uzmIv0x+BARkcW4UliBzOJK2Muk6M1p7NQEDD5ERGQx9txq7ekZ4QVnhZ3I1ZAlYvAhIiKLsedCbfAZwG4uaiIGHyIisgjlSjUSb+3GHtea21RQ0zD4EBGRRUi4VIgajYAwbydENONu7NQ0DD5ERGQR9rKbiwyAwYeIiMyeIAi68T1crZnuB4MPERGZvZTcUuSXKuEol6FHuJfY5ZAFY/AhIiKzV7dac+9IbzjIZSJXQ5aMwYeIiMzebu7GTgbC4ENERGatsFyJE5k3AHB8D90/Bh8iIjJru88XQBCAjoHu8Hd3FLscsnAMPkREZNbiU/IBAA+29RW5ErIGDD5ERGS2qlQaJFyuHdj8UDsGH7p/DD5ERGS2Ei4XorpGi0APR7T1dxW7HLICDD5ERGS2dt7q5nqonS8kEonI1ZA1YPAhIiKzpNEK2JXK8T1kWAw+RERklk5m3URhuQquDnboGcHVmskwGHyIiMgs7bjVzRXXujnkMt6uyDD4k0RERGZpR0oeAOBBzuYiA2LwISIis3PlejnSrlfATirBgNbNxC6HrAiDDxERmZ2d52u7uXpFeMPNQS5yNWRNGHyIiMjs7Eyp3ZSUixaSoTH4EBGRWSmuUOH41WIAwMC23JSUDIvBh4iIzMrO8/nQCkA7fzcEeTqJXQ5ZGQYfIiIyK3+cyQUADOngJ3IlZI0YfIiIyGyUVNUg4XIhAOCRjv4iV0PWiMGHiIjMxs6UfNRoBLTydUFkcxexyyErZCd2AWTbajRaXMgrw5lrJTidXYK8kiqoNFqo1FqoNALUGi2auyoQ4uWEYC8nhHg5IbK5C8J9nLlhIZEV+v1WNxdbe8hYGHzI5ArLldh6Ihu/n8lDSm4pVGrtXc8/d5tj/u4O6BPpg76tmqF3C294uyiMUywRmUxpdQ0OXGI3FxkXgw+ZhFqjxf5L17H5WBZ2nS+AWivonnN1sEOnIHd0DPRAuI8TFHYy2NtJYS+TQiaVIL+0GleLK5FZXInMokpcyC9Dbkk1vk/KxvdJ2QCAHuFeeCIqCI909Iezgj/WRJZo1/l8qDRaRDZ3QStfV7HLISvFOwQZlSAI+O1MLt7+MxVZxVW6452D3PFEdDD6RPog1NtJr26r6hoNEtOLkXC5EPsvXkdqXhkS04uRmF6MRT+fw9BO/ngyOhhRoZ7sDiOyIL+drt2b6xHO5iIjkgiCINz7NNtRWloKd3d3lJSUwM3NTexyLNrJrJt489cUHL96AwDg4STHiK6BGN09GG38DHdtc0uqsPXENXx/PAsZRZW6452D3PHP/i0wuL0fZFIGICJzVlZdg6g3d0Kl1uLPl/oa9HcE2YbG3r8tYlZXRkYGJk6ciPDwcDg6OqJFixZYuHAhVCpVvfMkEkmDx+rVq0Wq2nYVlFbjpU3JePzjgzh+9QYc5TK89GBLHJr7ABYOb2/wX2j+7o54MS4Se2YPwHf/jMGoqCAo7KQ4lV2CqRtP4IGVe/H1kauortEY9HOJyHB2pxZApdYiopkzWrObi4zIIrq6UlNTodVq8emnnyIyMhJnz57F5MmTUVFRgRUrVtQ7d926dRgyZIjua3d3d1OXa9N2p+Zj9venUVxRG0r/0S0IrwxuDT93B6N/tkQiQY9wL/QI98Lch9vgy0MZ+PLIVVwtqsTrP57FB7su4V/9W2BMzxA4yGVGr4eIGk83m6uDP7uoyagstqtr+fLlWLVqFa5cuaI7JpFIsG3bNjz++OONfh+lUgmlUqn7urS0FMHBwezq0pNSrcF//0jFuoMZAGqXmv/vPzqiU5CHqHVVqtT47lgW1hxIx7WbtWOMfN0UeDEuEqO7B0NhxwBEJLYKpRrd3tgBpVqL36f3RbsA/u4l/VlVV9ftlJSUwMvLq8HxadOmwcfHB927d8fq1auh1d59qvSyZcvg7u6uewQHBxurZKuVdr0cIz4+pAs9E3qHYduLsaKHHgBwsrfD+N7h2DN7AJaN7IhAD0fklyqx4KdzGLB8L745mokazd1/RojIuHanFkCp1iLM2wlt/dnNRcZlkS0+aWlp6NatG1auXIlJkybpjr/55psYOHAgHB0dsWvXLixYsADz5s3D66+/fsf3YovP/dmdmo9p3ySjUqWBl7M9lo/qhIFtfcUu646Uag2+P56Nj/dcRm5JNQAg1NsJMx9sheGdAzgImkgE//o6CX+czcPUAS3w6pA2YpdDFqqxLT6iBp9FixZh8eLFdz3n2LFjiI6O1n2dk5OD/v37o3///vj888/v+tqVK1diyZIlKCkpaXRNnNXVeF8ezsCin89BKwA9w73w4dNd4etm/LE8hlBdo8GmxEz8357LKCyvHY/U2tcVswa1wqB2vhxjQGQiZdU1iH5zJ5RqLX79dx90COS4TGoaiwg+hYWFKCwsvOs5YWFhcHCovZnm5OQgLi4OPXv2xPr16yGV3r2n7uDBg+jTpw/y8vLg69u4VggGn3vTaAW89ft5rE1IBwCMjg7GmyM6QC6zvJ7TCqUa6w9l4NN9aSitVgMAOgS6YdZDrRDXujkDEJGRfXc8C6/+cBotmjlj56z+/DdHTdbY+7eos7p8fHzg4+PTqHOvXbuGuLg4REVFYd26dfcMPQCQnJwMBwcHeHh43GelVKdKpcFLm5Ox/Vw+AOCVwa0xdUALi/1l5ayww4txkXi2Zyg+O5CGdQczcPZaKZ5ffxxdgj0w66FW6NvSx2K/PyJzt+3ENQDAyG5B/HdGJmER09lzcnIwYMAAhISEYMWKFbh+/bruOT+/2hU+f/nlF+Tl5SEmJgaOjo7Ys2cP5s+fjxdeeAEKBfdxMoTS6hpMWHcMSVdvwN5OihVPdMajnQPELssg3J3keGVwGzzfOxyf7b+CDYczcDLrJp77IhFdgj3wYlwkBrZpDinHABEZTM7NKhxJLwIAPNbFOn6XkPmziOATHx+Py5cv4/LlywgKCqr3XF1PnVwuxyeffIJZs2ZBq9UiIiICS5YswYsvvihGyVbnRoUKz32RiDPXSuDmYIe147uje1jDWXWWzttFgXmPtMXEvuFYvfcKNh69ipNZNzH5y+No5euCqQMiMayTP+wssFuPyNz8dDIHglC7116Qp5PY5ZCNsMhZXcbEMT4NFZRVY+znibiQXwZvZ3t8NbGnzayzcb1MiS8OpuOrw1dRrqwdAxTo4YixMaEYHR0MT2d7kSskskyCIGDw+/txMb8c/x3ZEU/1CBG7JLJwFjG42Rwx+NSXW1KFZ9YcxZXCCvi6KbBxUk9ENre9dTZKqmrw1eEMfHEwQ7cqtcJOise6BOC5mDDORCHS09lrJRj2UQLs7aQ4Nv9BuDvKxS6JLJxFDG4m85ZVXIkxnx9BVnEVAj0c8c3kngj1dha7LFG4O8ox7YGWmNgnAr+cysH6QxlIyS3Fd8ez8d3xbLQPcMOIroEY3jnAYqb0E4lpW3LtoOaH2voy9JBJscXnb9jiUyv7RiWe+uwIsm9UIdTbCd9M7oVAD0exyzIbgiDgROYNbDh0FX+czUWNpvafkVQCxLbwwaOdAzCgTTM0d2UIIvo7tUaLmP/uxvUyJT5/LhoPtjPfRU/JcrDFh5rs2s0qPL2mNvSE+zjj28m9TLLJqCWRSCSICvVCVKgXiiva47czufgx+RqSrt5AwuVCJFyuXZ+qfYAbBrRuhgGtm6NLsIdFrnVEZGgH04pwvUwJTyc5+rVqJnY5ZGMYfKienJtVeOqzw8gqrkKYtxNDTyN4OdtjbK9QjO0ViqziSvx08hriU/JxOrsE53JKcS6nFB/vSYPCTor2AW7oEuyJLiEe6BTojmAvJ26TQTZn24lsAMDwzgGwt+N/Bsi02NX1N7bc1VUbeo4gs7gSod5O2PRCL/i7s3urqQrLldh/8Tr2XLiOA5eu42ZlTYNz7GVSBHs5ItzHBeE+TvB3d4SPqwI+zvbwdlHA01kOZ3s7OMhlegckrVaASqNFjUaLGo0AtUaLGq0AjUZAjVYLrbbhP32ZVAJ7OynsZVLY20mhsJPBQS7lwnJkMBVKNaLf3ImqGg22TY1F1xBPsUsiK8GuLtJLzq3urcziSoR41bb0MPTcHx8XBUZ2C8LIbkEQBAHphRU4mXUTp7Ju4mTWTZzPK4NKrUXa9QqkXa+45/vZ20nhKJdBLqsLIRLU5RGtVoBaK0CjFaDWaqHW1H5tCHZSCVwd7ODqIIergx28nO3R3NUBzVwVaO6qgK+bA0K9nRDi7QQ3Bw5Spbv782weqmo0CPdxRpdgD7HLIRvE4EO1Y3putfQEezni2xd6IYADmQ1KIpEgopkLIpq5YGS32kU4NVoBuSVVSC+sQEZhBa4UVqCgTImiciWKylUoqlDhRqUKdW2yKrUWKrW2yTVIJYCdTAo7qQR2UglkUkm9lhxBqA1LKrUWKo1W97lqrYAblTW4cZsWq7/zcrZHiJcTWjRzQVt/V7T1d0MbP1d4u3D1dKr1bWImAOAf3QLZkkiiYPCxcdf+MqYnxKu2e4uhxzRkUgmCPJ0Q5OmEvi1vP8BTEAQo1VpUqTSoqtGgUqWBWlsbSgQBECBAEAC5TArZrUBjJ5PATlrbVSWXSSCX1XZd6bvdhlqjRVWNBuVKNcqq1SirrkFplRqF5UpcL1eioFSJ62VK5JZUIbO4EoXlKhRX1D5OZt2s916+bgp0DvJAVKgnokI90SHQHQ5yWVMvHVmo1LxSHL96A3ZSCZ6MDha7HLJRDD4mJAiCWf0PJ/tGJZ5eU7tOT+itgcwMPeZFIpHAQS6Dg1wGU4+EsJNJ4SqTwtVBDv9GrM9YrlTjalEFrhZV4mJ+GVJzy3A+rxRXiyqRX6pEfEo+4lNqN7eVyyToFOSB3pE+6NvShzPebMTGI7WtPYPa+6I517sikXBw898Ya3DzlqRsfJ+UhQ3P94DCTvz/6WYVV+qmrId5O+FbDmQmI6lQqpGSW4oTV2/gROYNJF29icJyZb1zXBR26BXhhQfa+OLBds25/pEVqlCq0fOtXShXqrFxUk/0jvQRuySyMtyyoomMEXxuVKgwYMVelFTV4OkeIVg2sqNB3repUvNK8dzaRBSUKRHm7YRNL8RwyjqZjCAIyCyuxJErRThwqRAHLxfWGz8kkQBdgz0wqL0fBrf3Q7iPba4Wbm2+TczEvK1nEO7jjF2z+uvd9Up0Lww+TWSsFp+9FwowYf0xCAJE3ZAv6WoxJqw7htJqNVr7uuLLiT24xQKJSqsVkJJbin0XryM+JR+n/jY+qHOQOx7rUrsdSDNXDpK2RIIgYNhHCTiXU4r5j7TF5H4RYpdEVojBp4mMuY7P/+2+hBXxF2Evk+K7KTEmn8q5J7UA/9qYhOoaLaJCPfHFuO5wd+L0YzIveSXV2HE+H/Hn8nAorQga7f+2A+kd6YMnooMxuL2vWXQZU+OcyrqJxz4+CHs7KY7MGwgvZ3uxSyIrxODTRMYMPlqtgClfJyE+JR/+7g745d994GOiab5bT2Tj1R9OQ60VENe6GT55JgqO9rxxkHkrLFfi11M5+PFkTr2ZYp5OcoyKCsLTPUIQ0cxFvAKpUV794RS+O56NEV0D8d7oLmKXQ1aKwaeJjL1yc1l1DR7/+CDSrlegZ7gXvp7U06izWWo0Wiz7PRVfHEwHAIzoGoh3RnXiDBqyOBmFFdiafA3fHctCXmm17nivCC9M6B2OB9v6cvsPM1RSWYOey3aiukaLLf+KQVSol9glkZVi8GkiU2xZcbmgHI9/fBDlSjWejA7CWyM6ws4IQSS/tBovbjyB41dvAACmDmiB2YNac1AhWTS1Rot9F6/jm6OZ2HOhAHULVId6O2FCbBieiA6Gs4IrdZiLdQfTsfiXFLTxc8UfM/qa1ZIeZF0YfJrIVHt1bT+XhylfJ0EQgD6RPvh4TDeDjrc5cqUI075JRmG5Eq4KO6x8sjMGtfcz2PsTmYOcm1X4+shVbDyaiZKq2plhbg52eLZXKCb2CeeK0SITBAEPvbcflwvK8cZj7TE2JkzsksiKMfg0kSk3KY0/l4eXNp9Epap235o1z0Ujsvn9jVcoqarB/+2+hC8OZkCjFdDGzxWrno3ilGCyapUqNbYkZWNtQjoyiioBAI5yGZ7pGYIX+kVwsTyR7Emtnc3qbC/D4dcGci83MioGnyYy9e7s53NLMWnDcVy7WQVXBzv835hu6N/q9tsX3I1ao8W3x7Lw3o6LKK5QAagdz7N0RAc42bPZn2yDVitgx/l8/N/uyzhzrQRA7eauY3qEYOqAFgxAJvbk6sNIzCjG5L7hmD+0ndjlkJVj8GkiUwcfoHbmyr++TsKxjBuQSIC41s0xNiYU/Vs2u+d4nCqVBvsuFuDdHRdxMb8cABDZ3AXzh7ZFXOvmpiifyOwIgoC9F6/jo12XcCLzJoDaFqDn+4Thn/1bsOXBBJKuFuMfqw5DLpPgwKsPcJFUMjoGnyYSI/gAgFKtwaKfz+HbxCzdsRAvJzzTMwSxLXzgaC+Fwk4GR3sZKpW1YWd3agEOpRVBeWvHbk8nOWY+1ApP9wjhrC0i1AagQ2lFWBF/Acm3ApCHkxxTB7TAczFh3CjViCZtOI6d5/MxOjoYb4/qJHY5ZAMYfJpIrOBT58r1cmw8monvj2ehtFrdqNcEejhiWGd/TO0fyQUJiW5DEATsSMnH8u0XcKmgtmU0wN0B8x5pi2Gd/DnTyMAu5pdh0Hv7IZEAO2f1RwuutUQmwODTRGIHnzpVKg1+PnUN3x3PRs7NKlTVaFCl0kCp1kIqAaJCPfFAG1880KY5Wvm68Bc3USNotAK2JV/Du/EXkFNSuxZQ9zBPLBjWHh2DGrEFPTXKrO9OYuuJa3i4gx9WPRsldjlkIxh8mshcgs+daLUCtIJglHV/iGxFdY0Gn+2/gk/2XkZ1jRYSCfBkVDBeGdLaZKupW6trN6vQ/509UGsF/PRib3Q28dY8ZLsae//m3dPCSKUShh6i++Qgl2H6wJbY/fIAPNYlAIIAbD6ehYEr9+HbxExotfz/YFN9fuAK1FoBsS28GXrILPEOSkQ2K8DDER881RVb/hWDdv5uKKmqwbytZ/DEp4eRmlcqdnkWp7hChU23Jmj8a0ALkashuj0GHyKyeVGhXvh5Wm+8PrQtnOxlSLp6A8M+TMCyP86jukYjdnkW44uEdFTVaNAh0A19In3ELofothh8iIgA2MmkmNQ3Ajtn9ceQ9n5QawV8uu8KHv7gAI5lFItdntnLK6nG5wlXAADT4lpywgWZLQYfIqK/CPBwxOqxUVjzXDR83RRIL6zAk58exqKfz6FS1bglJmzRyvgLqK7RIjrUE4Pb+4pdDtEdMfgQEd3GQ+18ET+zP56MDoIgAOsPZWDw+/txOK1I7NLMTkpOKX44kQ0AmD+0LVt7yKwx+BAR3YG7oxzvjOqML5/vgUAPR2QVV+HpNUfwxq8pHPtziyAIeOv38xAEYFgnf3QN8RS7JKK7YvAhIrqHfq2a4c+X+uLpHsEAgLUJ6Rj2UQLOZJeIXJn49l28joTLhbCXSTFnSBuxyyG6JwYfIqJGcHWQY9nITvhifDR8XBS4XFCOEZ8cxAc7L0Gt0YpdnijUGi3e+v08AGBcbCiCvZxErojo3hh8iIj08EAbX8TP7IeHO9TO/Hpv50U8+elhZBZVil2ayf2QlI2L+eVwd5RjWlxLscshahQGHyIiPXk52+OTZ7rhvdGd4aqww4nMm3jkwwPYkpQNW9kFqLS6Bit3XAQATB/Ykhskk8Vg8CEiagKJRIIRXYPw+4y+6B7miXKlGi9/fwrTvk1GSWWN2OUZ3eKfU3C9TIlwH2eM7RUqdjlEjcbgQ0R0H4K9nLDphRi8Mrg17KQS/HY6F0M+2I8jV6x32vv2c3nYciIbUgmwfFQn2NvxVkKWgz+tRET3SSaV4MW4SGz5VyzCvJ2QW1KNp9ccwTt/pqLGygY+F5Yr8drWMwCAF/q1QHSYl8gVEemHwYeIyEA6B3vgt+l9dYsefrI3DaNWHUJ6YYXYpRmEIAiYv+0MiipUaOPnipkPcUAzWR4GHyIiA3JW2OGdUZ3xyTPd4OZgh1PZJRj64QFsSsy0+IHPW09cw/Zz+ZDLJFj5ZGco7GRil0SkNwYfIiIjeKSjP/58qR96hnuhUqXB3K1n8MJXSSgqV4pdWpPk3KzCop/PAQBeerAV2ge4i1wRUdMw+BARGUmAhyO+mdwLcx9uA7lMgh0p+Rj8/gHsSS0QuzS9VKk0mLrxBMqUanQN8cA/+0WIXRJRkzH4EBEZkUwqwZT+LfDji73RytcFheVKTFh/DPO3nUG50vx3e1drtPj3t8k4mXUT7o5yvPtkF9jJeOsgy8WfXiIiE2gf4I6fp/XBhN5hAICNRzMx+L39OHi5UNzC7kIQBCz4+Rx2ns+Hwk6KteOiEe7jLHZZRPfFYoLPo48+ipCQEDg4OMDf3x9jx45FTk5OvXMyMzMxfPhwODs7w8fHB9OnT4dKpRKpYiKi+hzkMiwc3h7fTOqJIE9HXLtZhWc+P4rXzLT15+M9l/HN0UxIJMAHT3Xl1HWyChYTfOLi4vDdd9/hwoUL2LJlC9LS0jBq1Cjd8xqNBkOHDkVFRQUSEhKwadMmbNmyBS+//LKIVRMRNRQb6YPtL/XTrXj8za3Wn50p+SJX9j/fH8/CivjaLSkWDW+PIR38RK6IyDAkgoXOr/z555/x+OOPQ6lUQi6X448//sCwYcOQlZWFgIAAAMCmTZswfvx4FBQUwM3NrVHvW1paCnd3d5SUlDT6NURETXUorRBztpxGVnEVAGBgm+ZYOLw9QrzF2elcEAR8feQqFv2SAo1WwJT+LTD34Tai1EKkj8bevy2mxeeviouLsXHjRsTGxkIur90Y7/Dhw+jQoYMu9ADA4MGDoVQqkZSUdMf3UiqVKC0trfcgIjKV2Ba1rT//GtACcpkEu1IL8OB7+/DejouortGYtBalWoO5W87gPz+dg0Yr4ImoILw6uLVJayAyNosKPnPmzIGzszO8vb2RmZmJn376SfdcXl4efH19653v6ekJe3t75OXl3fE9ly1bBnd3d90jODjYaPUTEd2Ok70d5gxpgz9m9EOfSB+o1Fp8sOsS4lbsxddHrkKlNv62F/ml1XjqsyPYfDwLEgkw9+E2eGdUJ0ilEqN/NpEpiRp8Fi1aBIlEctfH8ePHdee/8sorSE5ORnx8PGQyGZ577rl6K6FKJA3/gQqCcNvjdebNm4eSkhLdIysry7DfJBFRI0U2d8FXE3vgk2e6IcDdAbkl1Xj9x7OIW7EX3xzNNFoAOnS5EMM/SkBy5k24Odhh3fjumNK/xV1/dxJZKlHH+BQWFqKw8O5TOcPCwuDg4NDgeHZ2NoKDg3Ho0CHExMRgwYIF+Omnn3Dq1CndOTdu3ICXlxd2796NuLi4RtXEMT5EZA6qazTYlJiJT/amoaCsdrXnIE9HjOkZgpFdg+Dn3vD3or6OZRTjvR0XcSitdif5Vr4u+GxsNMI4ZZ0sUGPv33YmrKkBHx8f+Pj4NOm1dXlNqaz9hRATE4OlS5ciNzcX/v7+AID4+HgoFApERUUZpmAiIhNxkMswvnc4nuoRgm+OZmLVvjRk36jCO39ewIrtF9A70gejooIwqJ0fHO0bv2eWVisgKfMGPtx1CQcu1f7HUy6T4KnuIZjzcBu4KES9LRAZnUXM6kpMTERiYiL69OkDT09PXLlyBQsWLEBubi7OnTsHhUIBjUaDLl26wNfXF8uXL0dxcTHGjx+Pxx9/HB999FGjP4stPkRkjqprNPj5ZA5+OJGNxPRi3XG5TILWfq7oFOSBzkHu6BDoDheFHSSQoK6nKq+0GscyinEsvRhJV2+gtLp2zSA7qQRPRAfjxbgWCPIUZxYZkaE09v5tEcHnzJkzmDFjBk6dOoWKigr4+/tjyJAheP311xEYGKg7LzMzE1OnTsXu3bvh6OiIMWPGYMWKFVAoFI3+LAYfIjJ3V4sqsOXENWxJysa1m1V6v97JXobhnQIw7YFIBHsx8JB1sKrgY0oMPkRkKQRBQPaNKpzOLsHp7Js4lX0TqXllUKm1EARAQO2vdzcHOaJCPREd5oXuYZ5o6+8GOffbIitjEWN8iIio6SQSCYK9nBDs5YShnfzFLofIIjDyExERkc1g8CEiIiKbweBDRERENoPBh4iIiGwGgw8RERHZDAYfIiIishkMPkRERGQzGHyIiIjIZjD4EBERkc1g8CEiIiKbweBDRERENoPBh4iIiGwGgw8RERHZDAYfIiIishl2YhdgbgRBAACUlpaKXAkRERE1Vt19u+4+ficMPn9TVlYGAAgODha5EiIiItJXWVkZ3N3d7/i8RLhXNLIxWq0WOTk5cHV1hUQiue/3Ky0tRXBwMLKysuDm5maACulOeK1Nh9fadHitTYfX2nSMca0FQUBZWRkCAgIgld55JA9bfP5GKpUiKCjI4O/r5ubGf0gmwmttOrzWpsNrbTq81qZj6Gt9t5aeOhzcTERERDaDwYeIiIhsBoOPkSkUCixcuBAKhULsUqwer7Xp8FqbDq+16fBam46Y15qDm4mIiMhmsMWHiIiIbAaDDxEREdkMBh8iIiKyGQw+REREZDMYfAzgk08+QXh4OBwcHBAVFYUDBw7c9fx9+/YhKioKDg4OiIiIwOrVq01UqeXT51pv3boVDz30EJo1awY3NzfExMRg+/btJqzWsun7c13n4MGDsLOzQ5cuXYxboBXR91orlUrMnz8foaGhUCgUaNGiBb744gsTVWvZ9L3WGzduROfOneHk5AR/f39MmDABRUVFJqrWcu3fvx/Dhw9HQEAAJBIJfvzxx3u+xmT3RoHuy6ZNmwS5XC6sWbNGSElJEWbMmCE4OzsLV69eve35V65cEZycnIQZM2YIKSkpwpo1awS5XC788MMPJq7c8uh7rWfMmCG8/fbbQmJionDx4kVh3rx5glwuF06cOGHiyi2Pvte6zs2bN4WIiAhh0KBBQufOnU1TrIVryrV+9NFHhZ49ewo7duwQ0tPThaNHjwoHDx40YdWWSd9rfeDAAUEqlQoffPCBcOXKFeHAgQNC+/bthccff9zElVue33//XZg/f76wZcsWAYCwbdu2u55vynsjg8996tGjhzBlypR6x9q0aSPMnTv3tue/+uqrQps2beod++c//yn06tXLaDVaC32v9e20a9dOWLx4saFLszpNvdajR48WXn/9dWHhwoUMPo2k77X+448/BHd3d6GoqMgU5VkVfa/18uXLhYiIiHrHPvzwQyEoKMhoNVqjxgQfU94b2dV1H1QqFZKSkjBo0KB6xwcNGoRDhw7d9jWHDx9ucP7gwYNx/Phx1NTUGK1WS9eUa/13Wq0WZWVl8PLyMkaJVqOp13rdunVIS0vDwoULjV2i1WjKtf75558RHR2Nd955B4GBgWjVqhVmz56NqqoqU5RssZpyrWNjY5GdnY3ff/8dgiAgPz8fP/zwA4YOHWqKkm2KKe+N3KT0PhQWFkKj0cDX17fecV9fX+Tl5d32NXl5ebc9X61Wo7CwEP7+/kar15I15Vr/3cqVK1FRUYEnn3zSGCVajaZc60uXLmHu3Lk4cOAA7Oz4a6WxmnKtr1y5goSEBDg4OGDbtm0oLCzE1KlTUVxczHE+d9GUax0bG4uNGzdi9OjRqK6uhlqtxqOPPoqPPvrIFCXbFFPeG9niYwASiaTe14IgNDh2r/Nvd5wa0vda1/n222+xaNEibN68Gc2bNzdWeValsddao9FgzJgxWLx4MVq1amWq8qyKPj/XWq0WEokEGzduRI8ePfDII4/g3Xffxfr169nq0wj6XOuUlBRMnz4dCxYsQFJSEv7880+kp6djypQppijV5pjq3sj/mt0HHx8fyGSyBv9bKCgoaJBc6/j5+d32fDs7O3h7exutVkvXlGtdZ/PmzZg4cSK+//57PPjgg8Ys0yroe63Lyspw/PhxJCcnY9q0aQBqb86CIMDOzg7x8fF44IEHTFK7pWnKz7W/vz8CAwPh7u6uO9a2bVsIgoDs7Gy0bNnSqDVbqqZc62XLlqF379545ZVXAACdOnWCs7Mz+vbtizfffJMt9AZkynsjW3zug729PaKiorBjx456x3fs2IHY2NjbviYmJqbB+fHx8YiOjoZcLjdarZauKdcaqG3pGT9+PL755hv2yzeSvtfazc0NZ86cwcmTJ3WPKVOmoHXr1jh58iR69uxpqtItTlN+rnv37o2cnByUl5frjl28eBFSqRRBQUFGrdeSNeVaV1ZWQiqtf5uUyWQA/tcaQYZh0nujwYdL25i66ZFr164VUlJShJdeeklwdnYWMjIyBEEQhLlz5wpjx47VnV83ZW/mzJlCSkqKsHbtWk5nbyR9r/U333wj2NnZCR9//LGQm5ure9y8eVOsb8Fi6Hut/46zuhpP32tdVlYmBAUFCaNGjRLOnTsn7Nu3T2jZsqUwadIksb4Fi6HvtV63bp1gZ2cnfPLJJ0JaWpqQkJAgREdHCz169BDrW7AYZWVlQnJyspCcnCwAEN59910hOTlZt3SAmPdGBh8D+Pjjj4XQ0FDB3t5e6Natm7Bv3z7dc+PGjRP69+9f7/y9e/cKXbt2Fezt7YWwsDBh1apVJq7Yculzrfv37y8AaPAYN26c6Qu3QPr+XP8Vg49+9L3W58+fFx588EHB0dFRCAoKEmbNmiVUVlaauGrLpO+1/vDDD4V27doJjo6Ogr+/v/DMM88I2dnZJq7a8uzZs+euv3/FvDdKBIHtdURERGQbOMaHiIiIbAaDDxEREdkMBh8iIiKyGQw+REREZDMYfIiIiMhmMPgQERGRzWDwISIiIpvB4ENEREQ2g8GHiGzWokWL0KVLF6N+RkZGBiQSCU6ePGnUzyGixmHwISK9jB8/HhKJBFOmTGnw3NSpUyGRSDB+/HjTF3YPEokEP/74Y71js2fPxq5duwz2GePHj8fjjz9e71hwcDByc3PRoUMHg30OETUdgw8R6S04OBibNm1CVVWV7lh1dTW+/fZbhISEiFiZflxcXODt7W3Uz5DJZPDz84OdnZ1RP4eIGofBh4j01q1bN4SEhGDr1q26Y1u3bkVwcDC6du2qO/bnn3+iT58+8PDwgLe3N4YNG4a0tLR673Xo0CF06dIFDg4OiI6Oxo8//liva2jv3r2QSCTYtWsXoqOj4eTkhNjYWFy4cKHe+/zyyy+IioqCg4MDIiIisHjxYqjVagBAWFgYAGDEiBGQSCS6r//e1SWRSBo86s7VaDSYOHEiwsPD4ejoiNatW+ODDz7QvXbRokXYsGEDfvrpJ91r9+7de9uurn379qFHjx5QKBTw9/fH3LlzdbUCwIABAzB9+nS8+uqr8PLygp+fHxYtWqTPXxER3QGDDxE1yYQJE7Bu3Trd11988QWef/75eudUVFRg1qxZOHbsGHbt2gWpVIoRI0ZAq9UCAMrKyjB8+HB07NgRJ06cwBtvvIE5c+bc9vPmz5+PlStX4vjx47Czs6v3Wdu3b8ezzz6L6dOnIyUlBZ9++inWr1+PpUuXAgCOHTsGAFi3bh1yc3N1X/9dbm6u7nH58mVERkaiX79+AACtVougoCB89913SElJwYIFC/Daa6/hu+++A1Dbbfbkk09iyJAhuveIjY1t8BnXrl3DI488gu7du+PUqVNYtWoV1q5dizfffLPeeRs2bICzszOOHj2Kd955B0uWLMGOHTvu/BdCRI1jlD3fichqjRs3TnjssceE69evCwqFQkhPTxcyMjIEBwcH4fr168Jjjz0mjBs37ravLSgoEAAIZ86cEQRBEFatWiV4e3sLVVVVunPWrFkjABCSk5MFQRCEPXv2CACEnTt36s757bffBAC61/Xt21d466236n3WV199Jfj7++u+BiBs27at3jkLFy4UOnfu3KBOrVYrjBgxQoiKihIqKyvveC2mTp0q/OMf/2hwbf4qPT293vfz2muvCa1btxa0Wq3unI8//lhwcXERNBqNIAiC0L9/f6FPnz713qd79+7CnDlz7lgLETUOO52JqEl8fHwwdOhQbNiwAYIgYOjQofDx8al3TlpaGv7zn//gyJEjKCws1LX0ZGZmokOHDrhw4QI6deoEBwcH3Wt69Ohx28/r1KmT7s/+/v4AgIKCAoSEhCApKQnHjh3TtfAAtV1T1dXVqKyshJOTk17f22uvvYbDhw/j2LFjcHR01B1fvXo1Pv/8c1y9ehVVVVVQqVR6zwo7f/48YmJiIJFIdMd69+6N8vJyZGdn68ZI/fX7rfueCwoK9PosImqIwYeImuz555/HtGnTAAAff/xxg+eHDx+O4OBgrFmzBgEBAdBqtejQoQNUKhUAQBCEegGg7tjtyOVy3Z/rXlMXpLRaLRYvXoyRI0c2eN1fQ1VjfP3113jvvfewd+9eBAUF6Y5/9913mDlzJlauXImYmBi4urpi+fLlOHr0qF7vf7fv+a/H//r91j1X9/0SUdMx+BBRkw0ZMkQXYgYPHlzvuaKiIpw/fx6ffvop+vbtCwBISEiod06bNm2wceNGKJVKKBQKAMDx48f1rqNbt264cOECIiMj73iOXC6HRqO56/scPnwYkyZNwqeffopevXrVe+7AgQOIjY3F1KlTdcf+PlDb3t7+np/Rrl07bNmypV4AOnToEFxdXREYGHjX1xLR/ePgZiJqMplMhvPnz+P8+fOQyWT1nvP09IS3tzc+++wzXL58Gbt378asWbPqnTNmzBhotVq88MILOH/+PLZv344VK1YAQINWkbtZsGABvvzySyxatAjnzp3D+fPnsXnzZrz++uu6c8LCwrBr1y7k5eXhxo0bDd4jLy8PI0aMwFNPPYXBgwcjLy8PeXl5uH79OgAgMjISx48fx/bt23Hx4kX85z//aTBIOiwsDKdPn8aFCxdQWFiImpqaBp8zdepUZGVl4d///jdSU1Px008/YeHChZg1axakUv5KJjI2/isjovvi5uYGNze3BselUik2bdqEpKQkdOjQATNnzsTy5csbvPaXX37ByZMn0aVLF8yfPx8LFiwAoF8X1eDBg/Hrr79ix44d6N69O3r16oV3330XoaGhunNWrlyJHTt2NJhyXyc1NRX5+fnYsGED/P39dY/u3bsDAKZMmYKRI0di9OjR6NmzJ4qKiuq1/gDA5MmT0bp1a0RHR6NZs2Y4ePBgg88JDAzE77//jsTERHTu3BlTpkzBxIkT64U0IjIeiXCnDnUiIhFs3LgREyZMQElJSb2BxUREhsAxPkQkqi+//BIREREIDAzEqVOnMGfOHDz55JMMPURkFAw+RCSqvLw8LFiwAHl5efD398cTTzxRb1o6EZEhsauLiIiIbAYHNxMREZHNYPAhIiIim8HgQ0RERDaDwYeIiIhsBoMPERER2QwGHyIiIrIZDD5ERERkMxh8iIiIyGb8PyxmWq/H+ROtAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "plt.xlabel(\"Magnetization\")\n", "plt.ylabel(\"Free energy eV\")\n", @@ -286,7 +326,7 @@ "id": "b8bf7959-f89f-43a1-ba9f-ca4e30678c2d", "metadata": {}, "source": [ - "The free energy minimum shows the most stable state. In this case it is around 0.5. Metadynamics, however, often requires a meticulous sampling, meaning with the standard energy increment of 0.001, you might overshoot and miss the global minimum. For a real measurement, it is recommended to make it as small as possible, meaning you should estimate the amount of time needed to run your calculation with Metadynamics with a small energy increment value and see what would be an acceptable total computation time." + "The free energy minimum shows the most stable state. In this case it is around 0.6. Metadynamics, however, often requires a meticulous sampling, meaning with the standard energy increment of 0.001, you might overshoot and miss the global minimum. For a real measurement, it is recommended to make it as small as possible, meaning you should estimate the amount of time needed to run your calculation with Metadynamics with a small energy increment value and see what would be an acceptable total computation time." ] }, {