{ "cells": [ { "cell_type": "code", "execution_count": 5, "id": "1d08e1d3-0f62-467e-8217-9eaf458afeb4", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA7gUlEQVR4nO3dZ3hUBeK28fvMpBKS0AOBAKFJr6mABUEREcWCgnQBFUFA1oaKropi17WB9LIINkQs6CIqiEISCKH3IhFIQs2kkDZz3g/7X95lRWmTOZnJ87uufMgwOM+VBebec+bMGKZpmoiIiIh4iM3qASIiIlK+KD5ERETEoxQfIiIi4lGKDxEREfEoxYeIiIh4lOJDREREPErxISIiIh6l+BARERGP8rN6wP9yuVwcPnyY0NBQDMOweo6IiIhcANM0ycnJITIyEpvtr49tlLn4OHz4MFFRUVbPEBERkUuQnp5OnTp1/vI+ZS4+QkNDgX+PDwsLs3iNiIiIXAiHw0FUVNSZ5/G/Uubi4z+nWsLCwhQfIiIiXuZCXjKhF5yKiIiIRyk+RERExKMUHyIiIuJRig8RERHxKMWHiIiIeJTiQ0RERDxK8SEiIiIepfgQERERj1J8iIiIiEdddHysWrWKXr16ERkZiWEYLFmy5KxfN02Tp59+mlq1ahEcHEy3bt3YvXu3u/aKiIiIl7vo+MjLy6NNmza899575/z1V155hbfffpupU6eSlJRESEgI3bt3p6Cg4LLHioiIiPe76M926dGjBz169Djnr5mmyVtvvcVTTz3FLbfcAsC8efOIiIhgyZIl9O3b9/LWioiIiNdz62s+9u/fT0ZGBt26dTtzW3h4OPHx8axZs+acv6ewsBCHw3HWl4iIiLhfQbGTCYs38cm6dEt3uDU+MjIyAIiIiDjr9oiIiDO/9r8mT55MeHj4ma+oqCh3ThIRERFgT1Yuvd/7hYXJ6fx96VZO5RdZtsXyq10mTJhAdnb2ma/0dGtrTERExNekLHmfWe8+z46MHKpVDOSDgTFUqhBg2Z6Lfs3HX6lZsyYAmZmZ1KpV68ztmZmZtG3b9py/JzAwkMDAQHfOEBERESA/N5stM+4n7tQ3tDb8Kawbw2MDu1IjNMjSXW498hEdHU3NmjVZsWLFmdscDgdJSUkkJia686FERETkLxzYvo6sNzoRd+obnKZBav1hvHLvrZaHB1zCkY/c3Fz27Nlz5vv9+/eTlpZGlSpVqFu3LuPGjWPSpEk0btyY6OhoJk6cSGRkJL1793bnbhERETkH0+UiZck7tNo4iWCjiKNUJuv690js1NPqaWdcdHysW7eOLl26nPl+/PjxAAwePJg5c+bw6KOPkpeXx7333supU6fo3Lkz3377LUFB1peWiIiIL8stKGbb+/2Jc3wHBmwK6kDtofNoEVHH6mlnMUzTNK0e8d8cDgfh4eFkZ2cTFhZm9RwRERGvsO2wg9EfpnLTyfmM8VvMugYjiRvwPDa73SOPfzHP3259wamIiIh4luly8cmvW3nq298pKnHxWdhddO8+hIQOV1o97U8pPkRERLxUTvYJdk0fSkvHfoySZ7m2aR1e69OGKiHWXUZ7IRQfIiIiXmjPxtUELRlGBzODYsPOGwmn6XFzDDabYfW081J8iIiIeBHT5SL545dpt/01AowSjlCd7F4f0DOmq9XTLpjiQ0RExEtknzzG3hlDiM/7GQzYUKETDYbPpVaV6lZPuyiKDxERES+Qln6KY3MG0825miLTTmrTvxF/1wQMm+WflHLRFB8iIiJlmGmazFy9n5e/3UF1553UrpCJX683SGh3ldXTLpniQ0REpIzKPp7JpwtnMOn3tgC0bdWS2rf9Slhw2b6a5XwUHyIiImXQjuTlVPrmfoZxjCT/R7my5wAGxNfFMMr+1Szno/gQEREpQ1xOJ8kL/k7M3nfxM1ykG5E81qcLDVvXs3qa2yg+REREyogTWYdInzWYhIIUMGBdaFeajphJxbDKVk9zK8WHiIhIGbBtzTKqffcAbThBgenPxtZPEnfrWK+8muV8FB8iIiIWcrlM3v9pD7tXrOYf/if4zVYH1+2ziG8Rb/W0UqP4EBERschRRwHjP9nIz7uPAR2JrxfKLf3uJyS0ktXTSpXiQ0RExAJbVi/F7/un2V7wKMH+VXjulhb0ielp9SyPUHyIiIh4kLOkhOS5jxN/cAY2w+SZsC9oOmwGjSNCrZ7mMYoPERERDzl6+Dey5gwgsWgTGJBcuSfdhn9AcEj5CQ9QfIiIiHjE5pWLqf3jWFrgIN8MZFuH54i7+X6rZ1lC8SEiIlKKSpwuvv5oCr12PonNMNlrjyag71xiGrexepplFB8iIiKl5Ej2acYs3MC2AxG0DKjJsWrxtBn+PkHBIVZPs5TiQ0REpBSs+2U5I753cfJ0CRUDQ9l981Ju6NDE6lllguJDRETEjYqLClk/6yESMhZwa/FAkmvfxbv92lO/Wvk+2vHfFB8iIiJucuS3XWTPH0RCyXYAutUu4bF7OxLoZ7d4Wdmi+BAREXGDDcsX0OCXR6hFHg4qsDfxJTp2H2z1rDJJ8SEiInIZigoLSJ05hoSsjwDY5deEiv3n0y66qcXLyi7f+6g8ERERDzl4PJ9HpyyiQ+anAKyN6Ef9R34mUuHxl3TkQ0RE5BIs23yERz/dRE5hdaoFDeWmzjEkdOtn9SyvoPgQERG5CAWn80idNY63fm9PjlmXDvUqM7TfJGpXCrZ6mtdQfIiIiFyg9D1bKFw4iI7Ovbzr/wuLEz5mfPfm+Nv1KoaLofgQERG5AOu+nk7T5IlUNE5zkjDyuzzLY11aWj3LKyk+RERE/kJBfi4bZ4wk/sRSMGCbf0uqDfknbWpHWz3Nayk+RERE/sT+A/sx5/Um3nUAl2mQFDWU2MEv4+cfYPU0r6aTVCIiIuewOPV3es3aTkZxBY4Tztauc0gc/qbCww105ENEROS/5Odm8/zXO1m44SgA8+tP5PmbW9Aqsp7Fy3yH4kNEROT/HNi+Dj4ZQrOiptiMoYzt2oTR1zbCbjOsnuZTFB8iIlLumS4XKV+8S6u05wk2igjxy6XZ3S8R27yR1dN8kuJDRETKtbycU2yfPpw4x3IwYFNQB2oPnUdsRB2rp/ksxYeIiJRb+7Yk4ffZUGLMQ5SYNlIajCR+wPPY7Harp/k0xYeIiJQ7pmmyaM0erv3uLiKMk2RRheM3TiUxvrvV08oFxYeIiJQrOQXFTFi8ma82HeE621BGhf9CvXvm0qx6LaunlRuKDxERKTf2bPyF95el8NWpxvjZDGJvGEDrThOx6bNZPErxISIiPs90uUj6+BXab3+VpwhiX/ibTLz7OjrUq2z1tHJJ8SEiIj4t++Qx9s4YSkLeKjDgtwqtmTvsasKrKjysovgQERGftSt1JRW/HEF7M5Mi007qFeOJ7/sEhk2nWayk+BAREZ9julwkLXqB9jvfJMBwctiIIPfm6SS0v9rqaYI+WE5ERHzMqfwiRsxP5cC2dQQYTlJDriJkzK80UXiUGTryISIiPmP9gRM8uHADh7MLSLYPpU6rbnS6dZROs5Qxig8REfF6LqeTpA+f5fTuVRwp+hv1q1bk3bvb07J2uNXT5BwUHyIi4tVOHj3CwVmDSDydDDZ4osF++g4aSWiQv9XT5E8oPkRExGttW/st1b4dSRtOUGD6s6nVEwy/bbROs5Rxig8REfE6LqeTpPlPEbt/Kn6Gi4O22pTcNpu4lvFWT5MLoPgQERGvcjSnkE1Th9I172swICX8epoPn05IaCWrp8kFUnyIiIjX+HXvMcYuSqNGbic6BPzErraPE3uLTrN4G8WHiIiUec6SEj75cilPJAfiMqFyREuO3bGeuCh9Eq03UnyIiEiZduzwb2TMGcRthZtZwLM0i7mKZ29uSXCA3eppcokUHyIiUmZtXvU5kT+MpSXZ5BPIk1eGkXBjG6tnyWVSfIiISJlTUlxEypxHif99DjbDZJ+tPn5955LQpK3V08QN3P4KHafTycSJE4mOjiY4OJiGDRvy/PPPY5qmux9KRER8UObve9n1ShcSD83GZpgkVb2FyId/oa7Cw2e4/cjHyy+/zJQpU5g7dy4tWrRg3bp1DB06lPDwcMaMGePuhxMRER/y444s1n30No+YW8g1g9kZN4n4nsOtniVu5vb4+PXXX7nlllvo2bMnAPXr12fhwoUkJye7+6FERMRHFDtdvPbdTj5YtQ+DbjSslE3cHQ/ToVFLq6dJKXD7aZeOHTuyYsUKdu3aBcDGjRtZvXo1PXr0OOf9CwsLcTgcZ32JiEj5kXFwN6tfvYN5q7YBMCgxmp4Pz6SOwsNnuf3Ix+OPP47D4aBp06bY7XacTicvvPAC/fv3P+f9J0+ezLPPPuvuGSIi4gXSln9I9C8P04U8ngmyE3772/Ropffu8HWG6eZXgi5atIhHHnmEV199lRYtWpCWlsa4ceN44403GDx48B/uX1hYSGFh4ZnvHQ4HUVFRZGdnExYW5s5pIiJSRhQVFpA6aywJmYsA2OXXhIr95xMZ3dTiZXKpHA4H4eHhF/T87fb4iIqK4vHHH2fUqFFnbps0aRL//Oc/2bFjx3l//8WMFxER73N4/w5yFwykScm/T8+vjehL+3v+QUBgkMXL5HJczPO320+75OfnY/uf99i32+24XC53P5SIiHiZ5B+W0HTVSCLJJ5sQ9nd6jYTr7rZ6lniY2+OjV69evPDCC9StW5cWLVqwYcMG3njjDe655x53P5SIiHiJgmInk7/ZzrI1uXwT6McO/2ZUGjSftnUbWz1NLOD20y45OTlMnDiRzz//nKysLCIjI+nXrx9PP/00AQEB5/39Ou0iIuJbDqYfZOTnv7H18L+vZnwizmDoTV3xDwi0eJm4k6Wv+bhcig8REd+x/uvpXJE8kUeK72VtUGfeuLMtXZrWsHqWlAJLX/MhIiJSkJ/LxpkPEH/8CzBgWOhann5gArXCg62eJmWA4kNERNzq4K40ShYNJt51AJdpkFRnCLFDXsHP//yn3qV8UHyIiIjbrFs6hebrn6GCUchxwjnc9R8kXnWr1bOkjFF8iIjIZTtd5GTaR58xdu/jYMDWgDZEDJlPq8h6Vk+TMkjxISIil2V3Zg4PLEhld1YIIX49aNEgirhBk7H76SlGzk1/MkRE5JKYLhfJS9/n0fWV+K24MtVDA2l+1zskNqpu9TQp4xQfIiJy0fJyTrFtxgjis//FK0ZT3m/0Jq/37UC1inrvDjk/xYeIiFyUfVuS8Fs8lFjXIZymgdmgK7MHxmOz262eJl5C8SEiIhfEdLlI/uxN2myZTJBRTBZVONZjCgkJN1g9TbyM4kNERM4rJ/sEu2YMIz7nBzBgY3Acde+ZR/PqtayeJl7Idv67iIhIebblUDZ3fJBMxexdFJt21jYcS6uHv6WywkMukY58iIjIOZkuF/PXHmDS1zspcrp4LuxRJl5Xm4TYblZPEy+n+BARkT/IPnmMvTOHcuRkJEXOm+nWLILX+lxHpQp6i3S5fIoPERE5y67UlVT8cgTtzUya+QVQ59oR3H1tBwzDsHqa+AjFh4iIAP8+zZK06AXa73yTAMPJYaMGuTfPoH/7GKuniY9RfIiICNnHM9k/cwgJ+b+CARtCrqTB8DlEVq5m9TTxQYoPEZFybsP+TKrOu5q25hGKTD82NH+EuD6PYth0QaSUDsWHiEg55XKZTP95H69+t5N+RnfuC/iWgltnEt+ms9XTxMcpPkREyqGTR4/w2uLVLNhf4d/ftxpEeK9nqRNWydphUi4oPkREyplta7+l2rcjudflxzK/l3i4Vyz94qJ0NYt4jOJDRKSccDmdJM+fSMz+KfgZLg761ebjOxvRqGldq6dJOaP4EBEpB45n/s7hWQNJKEwFA9aFXUezETMICa1k9TQphxQfIiI+busvXxGxfBStOMVpM4DNbScSe8toXc0illF8iIj4KKfL5J0fdtNq5Wu0sJ/igC0K+swhrpneNEyspfgQEfFBWTkFjFuUxq97j1OF+3irxg/EDn2d4IphVk8TUXyIiPiazau+YMOPn/Lr6b5UCLDzVO8ruap9P6tniZyh+BAR8RElxUWkzHmU+N/n0Mow+a1KC/oNGU2jGhWtniZyFsWHiIgPyDq0n2NzB5BYtAUMSK5yM48MH01QBYWHlD2KDxERL7fxx0+ou3I8zXGQZwaxPW4ScT1HWD1L5E8pPkREvFSx08Uvs5/gmt+nALDH3pDAfnOJadTK4mUif03xISLihQ6dOs2YhRsISQ/jKn+D5Oq30XbYOwQFh1g9TeS8FB8iIl5mZeo2xnz5O9mniwkN6sAv1y3lyk5XWT1L5IIpPkREvERRYQGps8bSNmMpoUUvUr9OY97p1566VStYPU3koig+RES8wOH9O8hdMJCEkl1gwJONDtJ18GAC/PQW6eJ9FB8iImXchu/m0nDN40SSj4MQ9nV6lR7X9bd6lsglU3yIiJRRhQX5pM0YTfyxzwDY6deUsIHzaFvvCouXiVwexYeISBl04FgeP856gqH5/w6PtbX602Hom/gHBFq8TOTyKT5ERMqYrzYd5vHPNlNU2IUrgjZQ8apRJFx7l9WzRNxG8SEiUkYU5Ofyw/wXGbs/ASd2YuvXILrfd9QKD7Z6mohbKT5ERMqAg7vSKFk0mBtdB9jrd4iCKyfwULcm+Nl1NYv4HsWHiIjF1i2dSvP1T1PBKOQEYVxz3S20urqp1bNESo3iQ0TEIqfzctg84z7iTn4NBmwNaE3EkPm0iqxv9TSRUqX4EBGxwIEdG+DjQcS5DuIyDZLqDidu8EvY/fTPsvg+/SkXEfGwT9alM++LdXxsy+CYUYmM694hsfPNVs8S8RjFh4iIh+QVFDFx6TYWpx4CInmnzjPcc8cttKwZZfU0EY9SfIiIeMD+rUmYn43gt9NDsBlXMP66JjxwzY3YbIbV00Q8TvEhIlKKTJeLlMVv0nrzZIKMYv4etJD8AcuIb1jN6mkillF8iIiUklzHCXZOH0Zczg9gwKagWOrcM5cqNRQeUr4pPkRESsGejb8QtGQYHcwjlJg21jUcTVz/v2Oz262eJmI5xYeIiBuZpslXy7+n+y99CTBKyKAap3pOJSHuOquniZQZig8RETdxFBTz+Geb+GZzIVP821Er1I/oYXOoWbWm1dNEyhTFh4iIG+xO+5nR32Wz86SBv91GRte3ueGqphg2fTaLyP/S3woRkctgulys/XAS9T6/hdG571CnUhCf3N+Rodc0V3iI/Akd+RARuUSO41nsmzmYhPxfwYDaYf58fX8c4aGhVk8TKdMUHyIil2DHuhVU+uo+2nKUItOPDc0eJu7Ox3S0Q+QCKD5ERC6Cy+kkeeFzdNj9Dv6Gk9+NmhT0nkF82yutnibiNUol0Q8dOsSAAQOoWrUqwcHBtGrVinXr1pXGQ4mIeMzJvCLGzvmJ6N1z8DecrK94DeHj1tBI4SFyUdx+5OPkyZN06tSJLl26sGzZMqpXr87u3bupXLmyux9KRMRjUg6cYMzCDRzJLuCE3xjGtnERe/vfdJpF5BK4PT5efvlloqKimD179pnboqOj3f0wIiIe4XI6SfrnRD7abXCkpCMNqoXw5N0jaB4ZZvU0Ea/l9vhYunQp3bt3p0+fPqxcuZLatWvzwAMPMGLEiHPev7CwkMLCwjPfOxwOd08SEbkkxzN/59DsQSQWrKe1PZDwpl14tE9nQgL1cjmRy+H244X79u1jypQpNG7cmO+++46RI0cyZswY5s6de877T548mfDw8DNfUVFR7p4kInLRtv7yNa4pnWldsJ7TZgDb2jzJ3+++VuEh4gaGaZqmO/+DAQEBxMTE8Ouvv565bcyYMaSkpLBmzZo/3P9cRz6ioqLIzs4mLEyHNUXEs5wlJSTPe4K436ZhN0x+s0Vh3jGb+s1jrZ4mUqY5HA7Cw8Mv6Pnb7Qlfq1YtmjdvftZtzZo147PPPjvn/QMDAwkMDHT3DBGRi5aVnUvm+zeRWLgBDEiudCMth0+lQsVwq6eJ+BS3n3bp1KkTO3fuPOu2Xbt2Ua9ePXc/lIiI26zefYwb31nL6rw65JuBpLSbTNy4hQoPkVLg9iMfDz30EB07duTFF1/kzjvvJDk5mWnTpjFt2jR3P5SIyGUrKS5i2r828OrqY5gmfBVxDzfe9BixjVtZPU3EZ7n9NR8AX331FRMmTGD37t1ER0czfvz4P73a5X9dzDkjEZHLkXVoP8fmDqSoIJ8+Rc9wR1wDnunVnCB/u9XTRLzOxTx/l0p8XA7Fh4h4wqYfPyVq5UNUxkGeGcS6Lv/k6muus3qWiNey9AWnIiJlWXFRIevnPEzC4XkA7LU3IKDfPK5upNMsIp6i+BCRciMjfQ+n5g0koXgbAEnVbqPNsHcJCg6xeJlI+aL4EJFy4fttmVT8eAgJbCPHDGZ3wmTiewy1epZIuaT4EBGfVlTi4pVvdzBj9X7qGUN5u+JcqvebSvsGzc//m0WkVCg+RMRnHT6wk08//ZAZx+IA6NoxkaY9hhDop6tZRKyk+BARn5T63XwarXmM0WY+m4Oeok+fAVzfoqbVs0QExYeI+JjCgnzSZj5I/NFPAdjpfwXPDexFrXoKD5GyQvEhIj7j0L6t5C8YRLxzDwBra/anwz1v4h+gz48SKUsUHyLiE9Yvm0OTtY9T2zjNSUI5eNVrJFzb1+pZInIOig8R8WoFxU4mfb2NwpStvOp/mh3+zak8+J+0qdPQ6mki8icUHyLitfZlZjNq0Sa2H3EAV3Nls7rceNd9+PkHWD1NRP6C4kNEvNK6pVMJX/8OhwufpmpIFd64qy1XN7nJ6lkicgEUHyLiVU7n5bB5xn3EnfwaDHi62k90vvdNIsKCrJ4mIhdI8SEiXuO3HesxPx5CnOsgLtMgqe4weg9+Gbuf/ikT8Sb6GysiXiFlyTu02PA8FYxCjlOJjG5vk3jlLVbPEpFLoPgQkTItv6iEZbOe5/aMt8CALYHtqDl0Hi1q1rV6mohcIsWHiJRZOzIcjFqQytGjzYkJrMGR+rcRO/AFnWYR8XL6GywiZY7pcvHTt59y/68VKSwxiQirSuYdP5HQpLbV00TEDRQfIlKm5DpOsmPGcLo4vudWczhHmtzFG3e2oWpFvUW6iK9QfIhImbF3068EfD6MGPMwJaaNXk3DSOwfi81mWD1NRNxI8SEiljNdLpI/fZW2W18l0Cgmk6qc7DmVTnHXWz1NREqB4kNELOU4dZw9M4YSn7sSDNgYnEC9YXNpWq2m1dNEpJQoPkTEMpt+P8X0+Qt5q2AVxdhZ32Qs8f0mYthsVk8TkVKk+BARjzNNkzm/HuDFb7ZT7IymfsXh3HTDjSTEXGv1NBHxAMWHiHhU9omjbJ91H/OO30ixWYvuLSIYfvtkwiv4Wz1NRDxE8SEiHrNz3Q+EfXUvCRzlLf8DpHX/jEEd62MYuppFpDxRfIhIqTNdLpI+fI4Ou9/G33ByyIigwm3/YHDbaKuniYgFFB8iUqpOHcvgwMzBJJxeCwakVryaRsNnU7tSVauniYhFFB8iUmo2b9pAjcW305bjFJr+pLV4jLg7/qarWUTKOcWHiLidy2UyddVe/vGvdD70q0KRPZCi22YR3zrR6mkiUgYoPkTErY5nHebhLw/w4+6TgJ0vmrzEozd3oGJYZauniUgZofgQEbfZ+us3VP/XKBJKElnjP5Bnb27BnTFRuppFRM6i+BCRy+YsKSF5/pPEHfgAu2HSPWATXYZ2oElUDauniUgZpPgQkctyLOMgGbMHkliYBgakVLqBFsOnUaFiuNXTRKSMUnyIyCXb8vMX1FwxhpacIt8MZGu7p4ntPdrqWSJSxik+ROSiOV0mU79dz8Ck+wkz8tlvq4ftzrnENm1n9TQR8QKKDxG5KJmOAsYs3EDS/hPstN3DgJoHaD18KkEVQq2eJiJeQvEhIhds00+f8faqdJLyGxMSYKfrbQ8Q17a21bNExMsoPkTkvEqKi0iZ/TcSD8/jRbMS+RHvMGnA1TSoXtHqaSLihRQfIvKXMtL3cHLeQBKLtwGwv1oXZg3vSlBwiMXLRMRbKT5E5E9t/GER9Vb9jWbkkmMGszv+ReJvvMfqWSLi5RQfIvIHRUXFpM58kITMhQDstjeiQv95tG/QwuJlIuIL9NGSInKW9BP53DktiYzDvwGQVKMPdR/5mdoKDxFxEx35EJEz/rX5dx7+bCuOghJeCrqPqM6Die/W1+pZIuJjFB8iQmFBPhtmjsHM2I+j+CHaRlXm3bvbUadyBauniYgPUnyIlHOH9m0lf8EgEpx7wA7Pt8yh7x034m/XWVkRKR2KD5FybP03s2mSNIHaxmlOUZHfrnydgV11mkVESpfiQ6QcKjidx8YZDxB/fAkYsN2/OZUHzadNVCOrp4lIOaD4ECln9h/LI+OD20gsXgvA2shBxAx5Db+AQIuXiUh5ofgQKUe+SDvEE4s306S4B9MCd3LkmtdIuOYOq2eJSDmj+BApB07n5TLvs8VM3lYVgIDoeFx90mhdpZK1w0SkXFJ8iPi433ak4vp4CIOdh/nc9jzXd+nKmGsb4aerWUTEIooPER+WvOQ9Wm54lgpGIceMSrx6U31adWxi9SwRKecUHyI+KD83my3T7yMuexkYsCWwLTWHzqdVzbpWTxMRUXyI+Jr921KwfTqUOFc6TtMgpf59xA58Abuf/rqLSNmgf41EfIRpmny8Lp2MpR8w1p7OUSqT1f19EjreaPU0EZGzlPorzl566SUMw2DcuHGl/VAi5VZuYQkPfZTGY59t5h/Ft7A0rB/2katpofAQkTKoVI98pKSk8MEHH9C6devSfBiRcm3v5rWkf/Ecy3JHYLcF8vD1zbnpqpuw2Qyrp4mInFOpHfnIzc2lf//+TJ8+ncqVK5fWw4iUW6bLRdLHr1Ln05u4puQXJoR8yUf3JjDymoYKDxEp00otPkaNGkXPnj3p1q3bX96vsLAQh8Nx1peI/DXHqeOkvnEr8dsmEWgUszE4nt4jXySmfhWrp4mInFepnHZZtGgRqamppKSknPe+kydP5tlnny2NGSI+aXfazwR/MZwOZgbFpp3UxmOIu3sihs1u9TQRkQvi9iMf6enpjB07lgULFhAUFHTe+0+YMIHs7OwzX+np6e6eJOITTNPkx8XTqPd5b+qYGRyhOvt6fUL8gL8rPETEq7j9yMf69evJysqiffv2Z25zOp2sWrWKd999l8LCQuz2//8PZWBgIIGB+jRNkb+SnV/Mo59tZOu2AL4OCGBrhXgaDJ9LrSrVrZ4mInLR3B4fXbt2ZfPmzWfdNnToUJo2bcpjjz12VniIyPlt2b6d+5ce4feTpwmwR/B954+4rWtnDJs+m0VEvJPb4yM0NJSWLVuedVtISAhVq1b9w+0i8udMl4ukhc/TftfbNCp+CFuVTrx3d3ta1Qm3epqIyGXRO5yKlEGnjmVwYOZgEk6vBQOGVdtBmwceIyzI3+ppIiKXzSPx8dNPP3niYUR8wo7k5VT65n7acoxC05+0Fo/S+Y6HdZpFRHyGjnyIlBEup5OkBX8ndu+7+Bku0o1Iim6bSXzrjlZPExFxK8WHSBlwPLeQGfNm81jW22DA+rCuXDF8JhXD9O7AIuJ7FB8iFkvad5wxizaQ6YgiMqA7TdskEHPrOJ1mERGfpfgQsYizpIQ1C57noR1XcNQMp2H1EOL6z+SKmqFWTxMRKVWKDxELHMtI58jsgXQu3MDrfq1Y2vIdnru1FRUC9FdSRHyf/qUT8bAtq5dS8/sHacUp8s1AgjvcxWu921k9S0TEYxQfIh7iLCkhee5jxB+cic0wOWCri3HnHGKbdrB6moiIRyk+RDwgKyOdo7PuJrFoExiQXPkmWg2fSnCIXt8hIuWPXk4vUspW7jrKHdNSqVCQSZ4ZxLoOrxA3doHCQ0TKLR35ECklJcVFvPH9Xt5fuQ/w45XqT/HEzW2IadzG6mkiIpZSfIiUgszf93By7kDy8mOB7gxIqMtTPW8gyF+f6iwiovgQcbONPyyi3qq/0ZRcxvqlk3D7g/Ro39jqWSIiZYbiQ8RNiosKWT/rIRIyFgCwx96Q4P7z6dFA4SEi8t8UHyJucOS3nTjmDyShZCcASdXvoO2wdwgMqmDxMhGRskfxIXKZvk/bQ/sl3bgCBw4qsDfxZeK7D7J6lohImaX4ELlERSUuJi/bzuxfDnCfvSe3B6+nYv/5tItuavU0EZEyTfEhcgkO7dvOC0s38E1GOACuxNHUv64RAYFBFi8TESn7FB8iFyl12Wwar53AQ2YVUoNf4oU74+naLMLqWSIiXkPxIXKBCk7nsXHmaOKPLQYDDgWEs2RIG2rWVniIiFwMxYfIBUjfs5nChYOJd+4FYE3kIGKGvIZ/QKDFy0REvI/iQ+Q81n01nWYpTxFiFHCSMNKvfpPELndYPUtExGvpg+VE/kRBsZMnPkujOGkWIUYB2wJaUTxiFa0VHiIil0VHPkTOYU9WLqM/TGVHRg4rjFG8Hr2ZhIHP4+cfYPU0ERGvp/gQ+R8pS95nY+oadhT1pVrFQF6/6wY6Nx5g9SwREZ+h+BD5P/m52WydcR+xp5YRa4OMOp24d/AQaoTqvTtERNxJ8SECHNi+Dj4ZQqwrHadpkFz/PiYMHI7dT39FRETcTf+ySrlmulykfP42rTa9QLBRxFEqc7T7eyR27Gn1NBERn6X4kHIrt7CEDVOGceWpJWDA5qAORA6dR/OIOlZPExHxaYoPKZe2HXYw+sNUGp+IJtHfxroGI4kb8Dw2u93qaSIiPk/xIeWK6XKxZFUKj31/gqISF6fDr2TbTb1JaNXW6mkiIuWG4kPKjZzsE+yccQ9XOdZRuWQyLZo25fU+bagcovfuEBHxJMWHlAt7Nq4maMkwYswMirHzYoc8utweg81mWD1NRKTcUXyITzNdLpI+epn2O14jwCjhCNXJ7vUBXWO6Wj1NRKTcUnyIz8o+eYy9M4aQkPczGLChQicaDJ9LrSrVrZ4mIlKuKT7EJ6Wln2Lr3IfpX/IzRaad1KZ/I/6uCRg2fZaiiIjVFB/iU0zTZObq/by0bAd+rpupXyGdGr2eJqHd1VZPExGR/6P4EJ+RfTyTH+a/yAsZ3TCxcX2rerS6/VvCgvytniYiIv9F8SE+YUfycip9cz+3cow9/oXU7DmBAfF1MQxdzSIiUtYoPsSruZxOkhf8nZi97+JnuEg3IrmtzyAatq5n9TQREfkTig/xWieyDpE+azAJBSlgwLrQrjQdMZOKYZWtniYiIn9B8SFeaWvScqovu5c2nKDA9GdT6yeJvXWsrmYREfECig/xKi6Xyfs/7eG773fyqb+D3+x1cN0+i7gW8VZPExGRC6T4EK9x9FQu4z/bys+7jwHRzKr3MoP63EFIaCWrp4mIyEVQfIhX2LJ6KVW+H8/JwrEE+zfmuVta0Cemp9WzRETkEig+pExzlpSQPPdx4g/OwGaYPBGylOojFtM4ItTqaSIicokUH1JmHT38G1lzBpBYtAkMSK58E+2GTyU4ROEhIuLNFB9SJm1euZjaP46lBQ7yzCC2xzxHXK/7rJ4lIiJuoPiQMqXE6eKTxR9x15aR2AyTvfZoAvrOJaZxG6uniYiImyg+pMw4kn2aMQs3sO5ARWr5tya4RjRthr9PUHCI1dNERMSNFB9SJqT9tISRK20cOW2nYmAAebfO45q29a2eJSIipUDxIZYqLipk/ayHSMhYwEMlVzOv9iO8d3d76lXV0Q4REV+l+BDLHPltF9nzB5FQsh2AerWq89mIeAIDAixeJiIipUnxIZbYsHwBDX55hFrk4aACexNfIr77YKtniYiIByg+xKOKCgtInTmGhKyPANjl14SK/efTLrqpxctERMRT9BGg4jEHj+cz4oN/0TjzGwDWRvSj/iM/E6nwEBEpV3TkQzxi2eYjPPrpJnIK7TwZPJb7O9UhoVs/q2eJiIgF3H7kY/LkycTGxhIaGkqNGjXo3bs3O3fudPfDiJcoOJ1H0rtDWbJwKjmFJXSoV5mJY0fTVuEhIlJuuf3Ix8qVKxk1ahSxsbGUlJTwxBNPcP3117Nt2zZCQnT5ZHmSvmczhQsHE+/cyxX+ITSN6cXoG9vjb9fZPhGR8swwTdMszQc4evQoNWrUYOXKlVx11VXnvb/D4SA8PJzs7GzCwsJKc5qUonVfT6dp8kQqGqc5SRjpV79J6y53WD1LRERKycU8f5f6az6ys7MBqFKlyjl/vbCwkMLCwjPfOxyO0p4kpaggP5eNM0YSf2IpGLAtoBXVBs+nde1oq6eJiEgZUarHv10uF+PGjaNTp060bNnynPeZPHky4eHhZ76ioqJKc5KUor2HjnL4tU7En1iKyzRYU+cemjzyAzUUHiIi8l9KNT5GjRrFli1bWLRo0Z/eZ8KECWRnZ5/5Sk9PL81JUkoWp/5Orw/Ws6qoCccJZ2vXuSQOfxM/f71bqYiInK3UTruMHj2ar776ilWrVlGnTp0/vV9gYCCBgYGlNUNKWX6eg1eWbmDOxlwAfmjwIDf1ep1WtepavExERMoqt8eHaZo8+OCDfP755/z0009ER+uQu686sH0dfDKEG0oqssB4gtFdmzH62kbYbYbV00REpAxze3yMGjWKDz/8kC+++ILQ0FAyMjIACA8PJzg42N0PJxYwXS5SlrxLq43PE2wUEWKrzMd31qZdu8ZWTxMRES/g9kttDePc/6939uzZDBky5Ly/X5falm15OafYPn04MY7lAGwK6kDtofOoGvHnp9ZERMT3WXqpbSm/bYhYaN+WJPw+G0qMeYgS00ZKg5HED3gem91u9TQREfEi+mwXOS/TNPkw6TdaL7ufVsYhsqjC8Runkhjf3eppIiLihRQf8pdyCoqZsHgzX206QmPjfl6s/BWNhk6jWfVaVk8TEREvpfiQP7Vn4y8s/vprvnJ0ws9m0OeGbnTofC82Xc0iIiKXQfEhf2C6XCR9/Artt7/KeFzsCavD/QP60b5uZauniYiID1B8yFmyTx5j74yhJOStAgPSKnTk1WG3E15V4SEiIu6h+JAzdqauJPTLEbQ3Myky7aReMZ74vk9g2Er1XfhFRKScUXwIpmny68KXiN35KgGGk8NGBLk3Tyeh/dVWTxMRER+k+CjnTuUX8fAnm6i96wid/J2kVryKRsNnE1mpmtXTRETERyk+yrHUfZmM/mgLh7MLCLDfwNXtY+hy0wCdZhERkVKl+CiHXE4nSR8+S7U9n+IofJb6Vavx7t3taVk73OppIiJSDig+ypmTR49wcNYgEk8ngwHP1N1Mj2FPUzFQfxRERMQz9IxTjmxb+y3Vvh1JG05QYPqzqdUT3HHbOJ1mERERj1J8lAMup5Ok+U8Ru38qfoaLg7balNw2m7iW8VZPExGRckjx4eOO5hSyevp4bnX8EwxICb+e5sOnExJayeppIiJSTik+fNive48xdlEarpxOtAv8nqy2o4m79UGrZ4mISDmn+PBBzpISvv5sNuPSInGZ0CSiNsV3rSEusorV00RERBQfvubY4d/ImDOIm4vSWGE8QFCHfvz95hYEB9itniYiIgIoPnzK5lWfE/nDWFqSTb4ZyIDEBsT2am31LBERkbMoPnxASXERKXMeJf73OdgMk322+vj1nUtsk7ZWTxMREfkDxYeXyzy0j+NzBpJYvAUMSK56C62HvU9QhYpWTxMRETknxYcX+3FHFp9+9AnvmVvINYPZGTeJuJ7DrZ4lIiLylxQfXqjY6eK173bywap9QEtaVB3JzbcOoEOjllZPExEROS/Fh5fJOLibQ/8cydc5A4DqDOlYn2E33kCgn65mERER76D48CJpyz8k+peH6UAeLwUWkXvHR9zQspbVs0RERC6K4sMLFBUWkDprLAmZiwDY5deEBv2nERmt8BAREe+j+CjjDu/fQe6CgSSU7AJgbURf2t/zDwICgyxeJiIicmkUH2XYrz9/T8sVA4kkn2xC2N/pNRKuu9vqWSIiIpdF8VEGFRQ7mfzNdhatyWFxQHXsARWoNGg+bes2tnqaiIjIZVN8lDHp+3bywJdH2HwkDwhgRYf3GdkjFv+AQKuniYiIuIXiowxZ//V0rkieSDdnDw6F3M3rd7ahyxU1rJ4lIiLiVoqPMqAgP5eNMx8g/vgXYEC34N3c9UAiNSvrLdJFRMT3KD4sdnBXGiWLBhPvOoDLNEiOGkLM4Ffw8w+wepqIiEipUHxYKGXpVFqsf5oKRiHHCedw13+QcNWtVs8SEREpVTarB5RHp4ucTFr4PS3XT6SCUcjWgDaY9/5MK4WHiIiUAzry4WG7M3N4YEEqu7MKcdiHcGtDiBs0Gbuf/qcQEZHyQc94HmK6XKR88S5vpJrsLm5E9dBAevd9jMSG1ayeJiIi4lGKDw/IyznF9unDiXMs53VbVZ5rOJ0X+nWmWkW9d4eIiJQ/io9Stm9LEn6LhxLjOoTTNEiP7suUgddgs9utniYiImIJxUcpMV0ukj97kzZbJhNkFJNFFY71mEJCwg1WTxMREbGU4qMU5OTmsOuDQcTn/AAGbAqKJWrYfJpXr2X1NBEREcvpUls323Iom17vp3DslIMS08bahmNp+ch3VFZ4iIiIADry4Tamy8WCNXt47pu9FDldvBk+lnrdQkiI7Wr1NBERkTJF8eEG2SePsXfmUMKyiylyjqZbs5q81qc1lSroLdJFRET+l+LjMu1KXUnFL0fQ3sykpc3O69cEclv3DhiGYfU0ERGRMknxcYlMl4ukRS/QfuebBBhODhsR5N48ndvbX231NBERkTJN8XEJso9nsm/mEBLyfwUDUkOuouHw2URW1ruVioiInI/i4yKtP3CCwLk30c7cRZHpx4bmjxDX51EMmy4cEhERuRCKjwvkcplM/3kfr363k1j68HLgHIpunUZ8m85WTxMREfEqio8LcPLoEaZ/spT3D9YBoFqb66h8y4OEVgi2eJmIiIj3UXycx7a131Lt25E8YOazwu9FBve6jn5xUbqaRURE5BIpPv6Ey+kkef5EYvZPwc9wcdBem/fvbE3DFnWtniYiIuLVFB/ncDzzdw7NHkRCwXowYF3YdTQbMYOQ0EpWTxMREfF6io//sfWXr6mxfBStOclpM4DNbScSe8toXc0iIiLiJoqP/+N0mbzzw26Cf/qQ+/xOcsAWBX3mENcsxuppIiIiPkXxAWTlFDBuURq/7j2OH3dSN7IW1wx+huCQUKuniYiI+JxyHx+bV33BiR/fJvn0GCoEBDKpdxt6tL/F6lkiIiI+q9ReyPDee+9Rv359goKCiI+PJzk5ubQe6pKUFBexZvo4WqwYzNXmOh6r9CNLR3fmtvZ1rJ4mIiLi00olPj766CPGjx/PM888Q2pqKm3atKF79+5kZWWVxsNdtKzf97HrlS4kHpqNzTBJqnoLA8dMolGNilZPExER8XmGaZqmu/+j8fHxxMbG8u677wLgcrmIioriwQcf5PHHH//L3+twOAgPDyc7O5uwsDB3T2Pjj59Qd+VDVCaHXDOYHXHPE9NzhNsfR0REpDy5mOdvt7/mo6ioiPXr1zNhwoQzt9lsNrp168aaNWv+cP/CwkIKCwvPfO9wONw9CYBip4sf5r1A999eA2CPvSFB/eYS06hVqTyeiIiInJvbT7scO3YMp9NJRETEWbdHRESQkZHxh/tPnjyZ8PDwM19RUVHungTAiu2ZTN5VixwzmKRqt1Hn4Z+po/AQERHxOMvfOWvChAlkZ2ef+UpPTy+Vx+neoiad4+NJ6fkt8aNnExQcUiqPIyIiIn/N7addqlWrht1uJzMz86zbMzMzqVmz5h/uHxgYSGBgoLtn/IFhGEzqrSMdIiIiVnP7kY+AgAA6dOjAihUrztzmcrlYsWIFiYmJ7n44ERER8TKl8iZj48ePZ/DgwcTExBAXF8dbb71FXl4eQ4cOLY2HExERES9SKvFx1113cfToUZ5++mkyMjJo27Yt33777R9ehCoiIiLlT6m8z8flKO33+RARERH3u5jnb8uvdhEREZHyRfEhIiIiHqX4EBEREY9SfIiIiIhHKT5ERETEoxQfIiIi4lGKDxEREfEoxYeIiIh4lOJDREREPKpU3l79cvznDVcdDofFS0RERORC/ed5+0LeOL3MxUdOTg4AUVFRFi8RERGRi5WTk0N4ePhf3qfMfbaLy+Xi8OHDhIaGYhiGW//bDoeDqKgo0tPT9bkxpUg/Z8/Qz9kz9HP2HP2sPaO0fs6maZKTk0NkZCQ221+/qqPMHfmw2WzUqVOnVB8jLCxMf7A9QD9nz9DP2TP0c/Yc/aw9ozR+zuc74vEfesGpiIiIeJTiQ0RERDyqXMVHYGAgzzzzDIGBgVZP8Wn6OXuGfs6eoZ+z5+hn7Rll4edc5l5wKiIiIr6tXB35EBEREespPkRERMSjFB8iIiLiUYoPERER8ahyEx/vvfce9evXJygoiPj4eJKTk62e5HMmT55MbGwsoaGh1KhRg969e7Nz506rZ/m8l156CcMwGDdunNVTfM6hQ4cYMGAAVatWJTg4mFatWrFu3TqrZ/kUp9PJxIkTiY6OJjg4mIYNG/L8889f0OeDyF9btWoVvXr1IjIyEsMwWLJkyVm/bpomTz/9NLVq1SI4OJhu3bqxe/duj2wrF/Hx0UcfMX78eJ555hlSU1Np06YN3bt3Jysry+ppPmXlypWMGjWKtWvXsnz5coqLi7n++uvJy8uzeprPSklJ4YMPPqB169ZWT/E5J0+epFOnTvj7+7Ns2TK2bdvG66+/TuXKla2e5lNefvllpkyZwrvvvsv27dt5+eWXeeWVV3jnnXesnub18vLyaNOmDe+99945f/2VV17h7bffZurUqSQlJRESEkL37t0pKCgo/XFmORAXF2eOGjXqzPdOp9OMjIw0J0+ebOEq35eVlWUC5sqVK62e4pNycnLMxo0bm8uXLzevvvpqc+zYsVZP8imPPfaY2blzZ6tn+LyePXua99xzz1m33XbbbWb//v0tWuSbAPPzzz8/873L5TJr1qxpvvrqq2duO3XqlBkYGGguXLiw1Pf4/JGPoqIi1q9fT7du3c7cZrPZ6NatG2vWrLFwme/Lzs4GoEqVKhYv8U2jRo2iZ8+eZ/3ZFvdZunQpMTEx9OnThxo1atCuXTumT59u9Syf07FjR1asWMGuXbsA2LhxI6tXr6ZHjx4WL/Nt+/fvJyMj46x/P8LDw4mPj/fIc2OZ+2A5dzt27BhOp5OIiIizbo+IiGDHjh0WrfJ9LpeLcePG0alTJ1q2bGn1HJ+zaNEiUlNTSUlJsXqKz9q3bx9Tpkxh/PjxPPHEE6SkpDBmzBgCAgIYPHiw1fN8xuOPP47D4aBp06bY7XacTicvvPAC/fv3t3qaT8vIyAA453Pjf36tNPl8fIg1Ro0axZYtW1i9erXVU3xOeno6Y8eOZfny5QQFBVk9x2e5XC5iYmJ48cUXAWjXrh1btmxh6tSpig83+vjjj1mwYAEffvghLVq0IC0tjXHjxhEZGamfsw/z+dMu1apVw263k5mZedbtmZmZ1KxZ06JVvm306NF89dVX/Pjjj9SpU8fqOT5n/fr1ZGVl0b59e/z8/PDz82PlypW8/fbb+Pn54XQ6rZ7oE2rVqkXz5s3Puq1Zs2YcPHjQokW+6ZFHHuHxxx+nb9++tGrVioEDB/LQQw8xefJkq6f5tP88/1n13Ojz8REQEECHDh1YsWLFmdtcLhcrVqwgMTHRwmW+xzRNRo8ezeeff84PP/xAdHS01ZN8UteuXdm8eTNpaWlnvmJiYujfvz9paWnY7XarJ/qETp06/eFS8V27dlGvXj2LFvmm/Px8bLazn4rsdjsul8uiReVDdHQ0NWvWPOu50eFwkJSU5JHnxnJx2mX8+PEMHjyYmJgY4uLieOutt8jLy2Po0KFWT/Mpo0aN4sMPP+SLL74gNDT0zHnD8PBwgoODLV7nO0JDQ//wOpqQkBCqVq2q19e40UMPPUTHjh158cUXufPOO0lOTmbatGlMmzbN6mk+pVevXrzwwgvUrVuXFi1asGHDBt544w3uueceq6d5vdzcXPbs2XPm+/3795OWlkaVKlWoW7cu48aNY9KkSTRu3Jjo6GgmTpxIZGQkvXv3Lv1xpX49TRnxzjvvmHXr1jUDAgLMuLg4c+3atVZP8jnAOb9mz55t9TSfp0ttS8eXX35ptmzZ0gwMDDSbNm1qTps2zepJPsfhcJhjx44169atawYFBZkNGjQwn3zySbOwsNDqaV7vxx9/POe/yYMHDzZN89+X206cONGMiIgwAwMDza5du5o7d+70yDbDNPU2ciIiIuI5Pv+aDxERESlbFB8iIiLiUYoPERER8SjFh4iIiHiU4kNEREQ8SvEhIiIiHqX4EBEREY9SfIiIiIhHKT5ERETEoxQfIiIi4lGKDxEREfEoxYeIiIh41P8DZ6h3xXW2gXQAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "from matplotlib import pyplot as plt\n", "\n", "x = np.linspace(0, 2*2.0*np.pi, 1000)\n", "y = np.sin(x)\n", "\n", "\n", "grid_bin = np.fromfile(\"../out/grid.bin\", dtype=np.float64)\n", "plt.plot(grid_bin, grid_bin)\n", "grid = np.loadtxt(\"../out/grid.dat\")\n", "plt.plot(grid, grid, linestyle='--')" ] }, { "cell_type": "code", "execution_count": null, "id": "8cda5f82-13ac-4b95-aed4-f999a8dbd185", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "ce5642b3-6387-45f7-aaaf-dfb9085b3d49", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "c4866585-9cb8-4eb8-a8c6-0489f831c55d", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.3" } }, "nbformat": 4, "nbformat_minor": 5 }