From 6d5ce32ca3955adf674ea1e294242cd475d0cf9a Mon Sep 17 00:00:00 2001 From: ParkSuMin Date: Tue, 23 Sep 2025 15:05:28 +0200 Subject: [PATCH] Add jupiter notebook for results' visualization --- graphics.ipynb | 89 ++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 89 insertions(+) create mode 100644 graphics.ipynb diff --git a/graphics.ipynb b/graphics.ipynb new file mode 100644 index 0000000..e8e008a --- /dev/null +++ b/graphics.ipynb @@ -0,0 +1,89 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 22, + "id": "18bd929b", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGzCAYAAAD9pBdvAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAAUFBJREFUeJzt3Qd4VMX6BvA3BULvSO+9NxFBERQiTQWvF6+AF0TgguJfFEXETkcRkQsWsIANUbyCiiAgUqRJV3qRLk0UqRJCcv7PO4cN6WST3T1b3t/zHLJ7sjmZTMLutzPffBNmWZYFEREREYeEO/WNRUREREjBiIiIiDhKwYiIiIg4SsGIiIiIOErBiIiIiDhKwYiIiIg4SsGIiIiIOErBiIiIiDhKwYiIiIg4SsGISIDZv38/wsLCMG3aNASzJUuWmJ+TH0UkuCkYEfEjDDD4Arxu3TqnmyIetnLlSrz00kv466+/nG6KiN9RMCISYMqVK4e///4b//73v51uirgZjAwdOlTBiEgqFIyIBBiOnOTIkQMRERHwVxcuXHC6CSISQBSMiARBzsgDDzyAPHny4LfffkOnTp3M7aJFi+LJJ59EXFxckq+Pj4/H66+/jlq1apmgplixYujbty9OnTqV5HFfffUVOnTogJIlSyIqKgqVKlXC8OHDU1yvZcuWqF27NtavX49bbrkFuXLlwjPPPJPuz7Bjxw7885//RKFChUwbrr/+enz99dcZ+vl/+ukntG3bFvnz5zffq0WLFlixYkWSx3A6hH20a9cu3H///eax7I/nn38e3Kj80KFD6NixI/Lly4fixYtj3LhxKb5PTEwMXnzxRVSuXNn8/GXKlMFTTz1lzifG7/PII49g9uzZph/4WPbtd999l6Q9gwYNMrcrVKhgvoYHf5ciomBEJGgwSGjTpg0KFy6MV1991bxI80V2ypQpSR7HwIMvjDfddBMmTJiAnj174pNPPjFfGxsbm/A4BjsMagYOHGge16hRI7zwwgt4+umnU3zvP/74A+3atUP9+vVNoHPrrbem2c6tW7fixhtvxPbt28212MbcuXObIGrWrFnp/ow//PCDCXjOnDljAoVRo0aZaY/bbrsNa9asSfH4f/3rXyb4GjNmDJo0aYIRI0aY9kVHR6NUqVJ4+eWXTbDBoG3ZsmUJX8evueuuu0w/3nnnnZg4caJp3/jx4801k1u+fDkefvhh3HfffXjllVdw8eJF3HPPPaZf6B//+Ae6dOlibvMaH330kTkYIIkIwHcJIuInpk6davG/5dq1a9N8zL59+8xj+FiXHj16mHPDhg1L8tgGDRpYjRo1Srj/448/msd98sknSR733XffpTh/4cKFFN+7b9++Vq5cuayLFy8mnGvRooX52rfffjtDP2OrVq2sOnXqJLlGfHy81axZM6tKlSoJ5xYvXmyuy4+ux/Dzbdq0MbcTt7NChQpWdHR0wrkXX3zRfO1//vOfhHOXL1+2SpcubYWFhVljxoxJOH/q1CkrZ86cpg9dPvroIys8PNz0V2L8GXndFStWJJzj/ezZs1t79uxJOPfzzz+b8xMnTkw4N3bsWHOOvz8RSUojIyJBpF+/fknuN2/eHHv37k24P3PmTDNlwZGBkydPJhwc9eAoyOLFixMemzNnzoTbZ8+eNY/j9ZgPwmmWxDg1wRGWa/nzzz/N6Ma9996bcE0eHEHgyMzu3bvNVFNqNm3aZD7ftWtX83jX154/fx6tWrUyIxsc0Uisd+/eCbeZY8PpIMYPvXr1SjhfoEABVKtWLUU/1ahRA9WrV0/STxyBocT9RK1btzbTWC5169Y1U0CJrykiaYtM53MiEkCYe5F82L9gwYJJckH4Yn769Glcd911qV7jxIkTSaZTnnvuORM8cFokMV4jMU55ZM+e/Zpt3LNnjwkGmLvBI6028HrJse3Uo0ePNK/PdvFndilbtmySzzMQYz8VKVIkxXnXlIrre3EaKa1plMT9lNr3Sa3vRSRtCkZEgkRGVtdw5ICBCHNEUuN68WUeBnNO+O5+2LBh5l0/X8Q3bNiAwYMHpxiBSDyKcq3vT8zR4EhIapjDkd7Xjh071uSmpIajO9fqk7T6yZ5xufq96tSpg9deey3VxzKZ1d1rikjaFIyIhBAGFd9//71JXk0vgGDVU44UfPnllyZh1GXfvn1Z+v4VK1Y0H7Nly2amNtzhmgZhgOTu17qL3+vnn3820z9c9eIJnrqOSDBSzohICGGuBlfdcIlucpcvX04oyOV6p5/4nf2lS5fw5ptvZun7c1SGS4EnT56Mo0ePpvj877//nubXMq+FQQJXuJw7d86tr81MPzF35Z133knxORacY56Ku7hiiFT0TCQljYyI+KH3338/SZ0KlwEDBmTpupx64dLe0aNHm4TQ22+/3YxSMEeCSZtcwsv6H82aNTM5D8zPePTRR827ei5F9cS0wxtvvIGbb77ZTIP06dPHjJYcP34cq1atwuHDh82IRGrCw8Px7rvvmiXErOPBhFnmljBoYEIpR0y++eYbeAKr237++ecmIZjX5kgSgzgm7vL8/PnzTTKsOxhM0bPPPmuWALPfuWzYFaSIhDIFIyJ+6K233kr1PIubZdXbb79tXhg5OsHiZJGRkShfvrwpDsYXXWKtkjlz5uCJJ54wSawMTPh5TlukleuRUTVr1jR777A0OmuZcDqIIyYNGjQwdUzSw1EVBi0c2Zk0aZIZIWHRMtYQYZDlKQx8WMSMNUE+/PBDU/+EBdYYODEgrFq1qtvXbNy4sWk3+5+BJvNSOO2lYEQECOP6XqcbISIiIqFLOSMiIiLiKAUjIiIi4igFIyIiIuIoBSMiIiLiKAUjIiIi4igFIyIiIhJYdUZYYIh7U8ybN8/s3sl9JKZOnZpmASCWlb711ltTnGf1RdYHyAiuxz9y5Ajy5s2rksoiIiIBgtVDuEN3yZIlTf0ejwQj3IGSRZEYXDAY4aZarNyYeJfMtOzcudNUSHRJa9fQ1DAQSb4xlYiIiASGQ4cOoXTp0p4JRl5++WUTFHAkxKVChQoZ+loGHwUKFMjQY2NiYszh4qrLxmqFHB3xlNjYWFPqmcEVSzOLd6iffUd97RvqZ99QPwd+P3NUhHHCtV673arAyjLOLAXN/SOWLl1q9oV4+OGHzf4SaXFN05QrV84EGLVr18ZLL72UUHY6Nfw8S0UnN336dFOSWURERPwf0zm6du2K06dPJ5kdyVIwkiNHDvNx4MCB6Ny5M9auXWv2aeBeC9xQK63pGQYkzClhMMKNrrjh1k8//YSGDRtmaGTkzJkzZkTm5MmT6f4wmYkGFy5ciOjoaEXdXqR+9h31tW+on31D/Rz4/czX7yJFilwzGHFrmoaJpAwqRo0aZe5zY6stW7akG4xUq1bNHC7cDfTXX381G1AxKElNVFSUOZJjJ3njD9Jb15Wk1M++o772DfWzb6ifA7efM3o9t5b2lihRwkzVJFajRg0cPHjQrcbdcMMN2LNnj1tfIyIiIsHJrZER5nlw2iWxXbt2mXwQd2zatMkENiIiTuDs9OXLlxEXF5fpYe3IyEhcvHgx09eQa1M/+38/R0REmK/NatkNt4KRxx9/3EyzcJrm3nvvxZo1azBlyhRzuAwZMsTUIvnwww/N/ddff91k0taqVcv8oMwZ+eGHH7BgwYIsNVxEJDMuXbpk6hwxsS4rwQzrJHG5omofeY/6OTD6mQtLOMCQPXt23wQjjRs3xqxZs0zAMWzYMBNkMNjo1q1bwmP4nzzxtA3/4z/xxBMmQGGD69ati++//z7VQmgiIt7EvDeWCOC7ORZh4pNnZp58eZ1z584hT5486RZykqxRP/t3PzOI4Wv877//bv5fValSJdO/J7crsN5xxx3mSMu0adOS3H/qqafMISLiND5x8omXq/OyUiaA1+C1uMJQL5Leo372/37OmTOnSVI9cOBAwjUyQ79dEQk5emET8a//T/ofKSIiIo5SMCIiIiKOUjAiIhLgmIQ7e/ZsBKoHHngAnTp1croZQfk7LV++vFlo4u8UjIiIBPiLNVcxtmvXzqdtClXcO61+/frwN9OmTUt1M1pu2/Kf//wH/s7t1TQikorYWODcuZTH2bOpn+d6/MqVgapVgSpVgCJF+FbI6Z9CAhRrRDiNyzxZMIsFsMR/FC1aFIFAIyMiyXHvyDVrgJEjgSefBPr1A+6/H+A709atgRtvBGrXBlh5uHBhbqZkBxeFCgFly3J7a+55ANx2G9CxI8A6PH37Ak88Abz4IjB2rH3tnj1Z1hi47jr7a/k1fCx3rJ4+nW9pgL/+cro3gv93ff68M0fG9yh1a0h///795v6XX35p6jlxCXO9evWwatWqJF+zfPlyNG/e3CzN5FLnRx99FOfZriu4dxj3IuPW7wx2uPPqiRMnEj7PDVD5febNm4dGjRqZ/cR4zdSwmBYLZfKde6FChdCxY0fTzvSWmo4ePRqVKlUyxbS4D9oXX3yR4nvPnz/ffI4/w2233Wbax/ZwmxJuysY2Jy5u57oua2Txa9gvqV130aJF5mdn37HQp6vyOEcfuKP8zz//bB7HI3k5i8TX4tYnuXPnNj83K5hz+avLW2+9ZX4+1rqpVq1amnu1JW7XX4meD1jJnOfYj/x8z549zWZ0rnZxBCe1aRrWAWP/s6YI+4i/l+PHj6cY+WF7+LX58+fHfffdh7N8Y+VFCmFFiC8MGzYAn39uH+k8UaaLQUmePEmPvHlTnuPoyO7d9nHokB10MPjgkRzf2XD0xDWK4rrNkZXcubP8o4c0vlDx95GJd3EpB8TdxL8BL/7+nn32Wbz66qumEBVvd+nSxewJxpELblbatm1bjBgxAu+//74pWvXII4+YY+rUqQklwocPH25eKPkiz93aOV00d+7cJN/n6aefNt+nYsWKKFiwYIp28Dpt2rRB06ZN8eOPP5rvz+/L7//LL7+kWrWTAcPHH3+MN9980wQjGzZswP3332/e5bdo0SLJC+ekSZNM0MAXVR4MiqZPn26KeN19992YOHEiBg8enOS63NyV/bJs2bJUr8v+GjdunDnfr18/PPjgg1ixYgX+9a9/mc1hv/vuO1O8k/hinRy3GuC0Wp8+ffDpp5+a+husWO4qsMfiodzxnkFC69atMWfOHBNMlC5dOlMFQRkw8VovvPBCQuDEYCM5BmOuQGTp0qWmnf379zd/G4nzU/j3wfts16lTp0y/jhkzBiP5JspbrABw+vRpvoUwHz3p0qVL1uzZs81H8R6/7ef4eMvauNGynn7asipWZDhy9ciVy7L++U/LGjTIsoYOtaxx4yxr8mTL+uQTy/rqK8tatMiyfvrJsrZts6yDBy3rzz8tKyYmc+24cMGyNm+2rP/9z7LGjLGsXr0sq3lzyypePGmbUjtKlrSsli0tq08fy3rlFSv2iy+sudOm+V9f+4m///7b2rZtm/lonDt37T721sHvnUE9evSwOnbsmObn+fw4a9Ysc3vfvn3m/rvvvpvw+a1bt5pz27dvN/d79epl/ec//0lyjR9//NEKDw+/2jfJrF271lzj7Nmz5v7ixYvNff7fTs9HH31kVatWzYrn/7crYmJirJw5c1rz589P8fNdvHjRypUrl7Vy5UorLi7OOnXqlPnINnfp0iXJ9/7+++8Trjl69Ghz7tdff00417dvX6tNmzYprpvYta777bffmnOufnnxxRetevXqpfsz//HHH+ZrlixZkurnmzVrZvXh/9lEOnfubLVv3z7V36mrXewLl40bN5pz/H3T1KlTrfz586f4XuXKlbPGjx9vbi9YsMCKiIiwDvI5K9nfxqJFi0w/8+djP505cybhMYMGDbKaNGmS8f9XmXj91siIhBb+H9+8+eoICEcmXHLmZIlh4N57gfbtueGCb9rE78tpHx7JcWjUNYLCY9euq7f/+AM4csQ+liwxD+d/aKYxWhyWbdMGuP12eyqIU0mSEn/HHKFwE99hnjlzxgxzZ7rgk5f/vrj1hotrY1KOcFSvXt1MM3BU4pNPPkl4DF//XOXyOc2xfv16M/LAx/LdMT/nGuZPvHs7pzPSw6/niAynexLjXmV8B54cH8uplejo6CTnObrAKZm0fsZixYqZERKO0CQ+xxGJrFw3cd+V5TRsBnAqiqNIHBHi9+PoB0cXXNfavn17iqTSm266CRMmTIA38ftySo6HC3+XnEbiprctW7Y05zg9k/j3xXYnnqLzBgUjEhq2br0agOzYcfU8Sxcz8GAAwkDE36Y9+ITQsKF9JPfnnykCFWvLFoTx2LSJk8rAyy/bL3p8kmFgwqN6dSXLurAfMvM75wszdzfl1/ppNVeW6HZxTQ+4AgpOYfTt29fkiSTHF1zmjvCFlAcDFk5XMAjhfb54J8aciPTwezGnJHHgk15yJR9P3377rXkRTLxnCqdg0vsZE993nUv8M7uuW6pUqSSPu9Z1yXWdjOJ0F/uXUzqfffYZnnvuOSxcuBA3MufMTeFX/sbsAZOr01/ekl4/eouCEQleDDo++8wOQLZtu3qeTzxcBukKQJK9YwsYTHpt0sQ+rrgcG4tF06ejNf9z//ADwN2xjx0DOM/vmuvnuyJXYNKqlZ2EKyGlYcOG2LZtGyoz7ygVmzdvxh9//GHyBFzvotetW5fp78UX4+uuu86MJF0L36kzOGDwwwTbLI9ApXLdxPkh7mKOC1cNZQRHXHhwc1nmzDCXhcEIR56Yg9KjR4+Ex65YsSLJiFNqQRuXcLvycpjA6m67+H2ZTMzD9Xvl3wETY5kb5CQFIxJcOI3hGgHhdIwLk+Q4bcEA5K67gAw8KQaqmAIFYHG0h090rmkpBiU8li2zE2bfe88++K6vceOrwQnftSV7VyT+gSslkr8AFS5cOMmQe0YxoZMvikxY7d27txnd4IsS37kzIZSjI3xxY/InEziZtMlk1szgru5jx441iZPc7Z1JmlxVwtU+3ESV9xPj9MCTTz6Jxx9/3CRYcmUHX2S5GohBSeIXcHckvi7f5d98882mTxkEuHNdTmFwKou/C7ad100+ssLPT5kyBXfddZfZHZpJpbt370b37t3N5wcNGmSmbRiocArnm2++Mf3hSopNjkEjf8+cNmMSKadUmGCbvF0c/eFKIK4S4pRV8s0g+b3q1KljfidMeGX/PvzwwyY4Sz5V5XNWAFACa2Dzej/v3m1Zo0ZZVv36SRMEIyMtiwlhH3xgWYkSv0K6r8+ft6zvvrOsgQMtq3btlEmVefNa1l13WdYbb9j9mijpMBikl2jnjsSJlb7ABE8+ByY/mHyZVgIrExxd2FaeYyKky5o1a6zo6GgrT548Vu7cua26detaI0eOTPj89OnTrfLly1tRUVFW06ZNra+//jrJdVNLqkzL0aNHre7du1tFihQx16tYsaJJ4HQ9pydP0GWy6+uvv24SX7Nly2YVLVrUJKIuXbo0ze+dWgJn8mTTzFw3eaIoE2Hvueceq0CBAuY8v29yx44dszp16mSVKFHCyp49u0kifeGFF5L8vbz55pumH9iOqlWrWh9++GGSayT+ndLy5cutOnXqWDly5LCaN29uzZw5M0m7qF+/flbhwoXNef7syRNY6cCBA9Zdd91lfud58+Y1ibNHjhxJ+HtOLUGXX8/reDOBNezKD+3XOEzH5VOMYjMyzJdRnHPjMrX27dunmCMT+H8/8x0+380sXnz1HAsusRYIR0BYFySVpYbBzO2+/u03YOFCe9SEH0+eTPr5ChXsEROOKrFuSirLGAMJkyb5rpV1JjK71bnHEljlmtTPgdHP6f2/yujrt6ZpJDDxxbNrV3tFSUSE/UL5r3/ZAYhyIDKOiXwPPGAfTFDbuPHqlM6KFRxvBiZPtg9OdXXoYBeA40et0BERD1EwIoGFL5gjRrDakT2xwFUmM2cCiZbzSSbxHVGjRvYxZIi95JVLhhmYzJ9v5+PMmmUf3AOjc2c7MLn5Zr9dUSIigUHPIBI4OArCd+Qsqc5AhOv0+e5dgYh3sIIjVxv9978Aqzr+8gvw1FP2aAorxr7zDsBVCez/Z59lEQOnWywiAUrBiAQGlknnKMh339m1QbgfBKcOsjDvL26qU8euW8L9Nbhs+MEH7VVJvD9qlL0nD0dVxo/nGkSnWysiAUTBiPg3joC89ZY9FXDwoL0fy08/2Ymr4gzm6HD/DC4NZg0T1nK58047eZj7+wwcCHC5JpNeP/44UxVOvS0A8vZFAoYn/j8pGBH/xR1EuS7/4YdZs9lOTmXhpUSlmsVhLGXPlUtff22PhkyaZNcqYW4Pc03+/W/W5LZzSziqdfmyo811rTBKvJOriGSN6/9TVlZLKoFV/BNzFO65xy7jznfiY8YATzyhMub+rEgRoH9/+9izB2D5b46MuG7zYGDSpYsdnHDazce/z4iICLMPh2ufDRaFcpX7dncpJMuic0mjlpx6j/rZv/uZIyIMRPj/if+v+P8rsxSMiP/54gs7H4GbxBUvbk8D3HKL060Sd3A6jYnGL7wAcKOyjz4CZswAjh8HuIkfD+6Rw6CkWzeWj/RZ04rzb+rKxmeZxSfhv//+Gzlz5sxUMCMZo34OjH5mIOL6f5VZCkbEf3Djp8GD7QRIYgDCF7ArO11KAOITm2v/HP5euUSYoyVffWXvHfTcc/bBOjEcUWGpfuaeeLVJYWYDNu6VktnNxvh1y5Ytwy233KKCiV6kfvb/fubjszIi4qJgRPwDK4GyaBmX6tKgQfYKDS+/MIkP8UmOS4V5nD4NfPmlHZiwgi5X5/Bg4mvfvkCfPvaUjhfxCTSzT6L8Ou7rwWqTepH0HvVz6PSzJuHEeXwxYv4AAxEuFeWL1CuvKBAJZiwr37MnsGiRXeWVRdaYc3L4MPD88/bOwqywy78JrXwRCXoKRsQ5XHExerS9lwzn77lKhqtl7r7b6ZaJL5UrZ4+CMRBhbglX43D65NNP7SXd3E10yhR7dZWIBCUFI+KMU6eAjh2BZ56xgxLujbJqFVClitMtE6dwrxsmtPLvYP16O4mZRe1+/tmeumHl18ces8vSi0hQUTAivsfN2Fipc84c+wWIZcXff5/rLJ1umfgLTtuxqBpzicaNAypVsvNMJkwAqlWzdxJmEqzDdUtExDMUjIhv8QWmaVM7T4DLOVeuBHr3Vv0QSV2hQnZFV46GzJtnJ7/yb2XhQrsIHoMUTvFkYZmuiDhPwYj4RHhMDCK4QoKBR0yM/aLC0uF8ByxyLSzE1LYt8M03wN699hLwwoXtLQK4SR8TXl1TPEp4FQk4CkbE+w4cwC2DByP8gw/sF5WRI+0h9oIFnW6ZBCKOqLEiLxNe+Td1ww32dgGs8NqsmT0FyBE4lXwXCRgKRsS7Ll1C5L33Iv/+/bCKFrX3K2HSqko7S1YxuZV7F3HjRO7qzKXCPMecJI7AsWYJtxDgrsIi4tf0iiDeNXQowjZuxKW8eXGZ+SGtWjndIglG119vJ0FztGTsWKBiRXvF1muv2XklrFnCIEVE/JKCEfGeH3+064gA2PTQQ3Y9CRFvYh7Jk08Cu3cD335r17CJi7NrljA/ife5e7DySkT8ioIR8Q4uw+T28ZaF+O7dcZRz+SK+wmnA9u3tVTccEeFmfCz9zoqv7drZBfaYb8JcExFxnIIR8Y5HH7Xn6itUQByHykWcUr++vQcOV+FwmXCePMCWLXahvQoV7K0HGDyLSOAEI7/99hvuv/9+FC5c2Gw3XKdOHaxjCe90LFmyBA0bNkRUVBQqV66MadOmZaXN4u8+/xz48EP73SnLe3O/GRGnlS1rF1A7dAh4+WWgZEngyBF7mTCXBjPZlZ8TEf8ORk6dOoWbbrrJ7Oo3b948bNu2DePGjUPBdJZo7tu3Dx06dMCtt96KTZs24bHHHkPv3r0xn1uJS/BhAmG/fvZtrpq56SanWySSVIECwFNP2YX3+Maodm3g7Fk72ZWJr5xeZAl6EfEZt7ZFffnll1GmTBlMnTo14VwFDnOm4+233zaPYdBCNWrUwPLlyzF+/Hi0adMm1a+JiYkxh8uZM2fMx9jYWHN4iutanrxmSIuPR0SPHgg/dQrxjRohjjuxJvqdqZ+9T33tBlZy5SqbLl0QNn8+wsePRzh3kOaUzscfI751a8Q//jgsJr0mqxCsfvYN9XPg93NGrxlmWRlPK69Zs6YJIA4fPoylS5eiVKlSePjhh9GHlTXTcMstt5gpmtdffz3hHIMZjpCcTmOe9qWXXsLQoUNTnJ8+fTpyaf8Sv1Xx669R5/33cTkqCkteew3nubGZSADJv2cPKn/1FUquWIFwbuDIXOzy5bGnY0f81rw5rEi33r+JhLwLFy6ga9eu5vU+XzpT9m4FIzlYUAjMARuIzp07Y+3atRgwYIAZ/ejRo0eqX1O1alX07NkTQ/gu+Yq5c+eaqRs2knknGRkZ4YjMyZMn0/1hMhOxLVy4ENHR0WbqSbJg82ZENm2KsEuXEPfGG4hPFKCqn31Hfe0h+/cjfOJEhL//PsLOnzenrNKlEf/II4jv3RuxOXOqn31Af8+B3898/S5SpMg1gxG3wvz4+Hhcf/31GMWNqQA0aNAAW7ZsSTcYyQwmuvJIjp3kjT9Ib103ZFy8aK9M4DLJO+5AxEMPISKVje/Uz76jvs6iKlWA//7XFO3D22+b22GHDyPi6acRMWoUwnv3Ro5atdTPPqJ+Dtx+zuj13EpgLVGihJmqSYw5IAe5WVUaihcvjuPHjyc5x/uMkFIbFZEAxI3KNm8GrrvO3hNEO/BKsGByPkd19++3/7Zr1OBbPUS89hpa9+2L8IcftpcMi0iWuBWMcCXNzp07k5zbtWsXyqVTWbNp06ZYxEJDiXA4iOclCPB366ojwidrBiQiwYYjtQ8+aNcnmTMH8TffjIjLlxHx7ruci7ZX4Gzb5nQrRUIjGHn88cexevVqM02zZ88ek1A6ZcoU9O/fP+ExzA3pzs2rrujXrx/27t2Lp556Cjt27MCbb76Jzz//3FxLAtyffwKu6Tku573jDqdbJOJdrJ3ToQPifvgBP44ahXiuCGS5ea7A4RLhe+4BNmxwupUiwR2MNG7cGLNmzcKnn36K2rVrY/jw4WaVTDeWWr7i6NGjSaZtuKz322+/NaMh9erVM0t833333TSX9UqAYN4zA5DffrPfGb76qtMtEvGpP2vWRNw33wAs+nj33fb/iS+/BBo1skvOL1/udBNFAobb69TuuOMOc6QlteqqLVu2xEbtmBlc+E5w5kyASx0/+QTIndvpFok4g8EHg5CtW+2NIbkpHzfj49GihZ1TlUqtEhG5SnvTiPuYzOeamnvpJXv7dpFQV6uWHaTv2gVwaTtXESxdCtx+O9CkCfDVV6YwoIikpGBE3MP5cSbrsXw2S70//bTTLRLxL5UqAVOm2KtsBgwAuGpw7VqgUyegXj1gxgz7/5GIJFAwIu7hDqecC8+b194Ej9uyi0hKpUsDrDzNkUQuD+b/Ga7G6dLFXiL8/vt2bR4RUTAibli/HnjhBfv2xIn29usikj4ud2ehyAMHgGHDgEKFgN27gV69gMqVgUmTgL//drqVIo5SMCIZc+ECwFVTly8D//wnkGj5tohksIDa88/bQQlXnxUvDhw6BPzf/9mBPUcdOf0pEoIUjEjGDBoEsOBdyZJ2eWytDBDJnDx5gCeeAPbtA958E2DRSFapHjzYvs2k8FOnnG6liE8pGJFr+/Zb+0mTuHS7cGGnWyQS+Ljx6EMP2VM2U6fa9XoYhHA/nPLl7SlRFhYUCQEKRiR9J07YZbDpsceA6GinWyQSXLgEmBtNspz8Z58BdeqY/W8wfLgdlHBqR0GJBDkFI5I2VpTs3dsOSFjqmgWdRMQ7uDLt3nuBTZuA//0PqFvXziEZMcIOSp57DvjjD6dbKeIVCkYkbe+8A7DcdfbsdpVVDiuLiPf3v/nHPwBWrWZlV9YmYVAycqQdlDzzDHDypNOtFPEoBSOSOlaRdG1myBERvksTEd8GJdzzhhvvzZoF1K8PnDtn/3/k6hvWLlFQIkFCwYikFBsL3H+/vZz3ttvsXBERcS4oYfVWBiWzZ18NSsaMsUdKWAX599+dbqVIligYkZSYOMfy1QUKAB98YD8ZioizuJy+Y0c7KOE+Nw0aAOfPAy+/bI+UcGmwghIJUHqVkaRWrrTnpmnyZLuktYj4V1By1112ReSvv7Z3DWZQwqJpHCl56ik76VwkgCgYkauYJMfpGe4sys3wmNkvIv6JQcmdd9qjmEw05+7ZnFodO9YeKWGhQgUlEiAUjMhV3GGUVSFZBZJ7z4hIYAQld9wBrFkDzJkDNG5sByUsOc+RkieftCu8ivgxBSNi4zw0q0DyiY278ebP73SLRMQd/L/boQPw00921eQbbrA34Bs3zh4pGTgQOHbM6VaKpErBiNjeeMP+eN99QPPmTrdGRLISlLRvD6xeDcybBzRpYgcl48cDFSva0zdKdBU/o2BE7FLT06fbt/v3d7o1IuKpoKRtW2DVKuC7764GJZy+4UgJi6epoqv4CQUjYk/PXLxoV3ps1szp1oiIp4OSNm3soGTuXDvRlatvXMXTuCHfX3853UoJcQpGQh1Xzrz11tVRET5xiUjw4f/tdu3sRFfWKXGVmXdtyMeP3KBPxAEKRkLdggXAr7/aCatduzrdGhHxVZ0SJq1/8QVQqxZw+rQ9QsKRElZ2ZYVXER9SMBLqXImr3MI8d26nWyMivsLKyvfcA/zyCzBjBlC9up0/xj1vGJQwt4RLhEV8QMFIKNu/314CSA8/7HRrRMSpoORf/wK2bLGX9VeubG/Ax1U3XH0zYYKdUybiRQpGQtnbbwOWBURHA1WrOt0aEXFSRIRdgXn7duD99+08EhZL40aZlSrZo6gxMU63UoKUgpFQxXc6775r39aoiIi4REYCPXsCO3cCU6YAZcoAR44AjzwCVKlin+PO3iIepGAkVM2cadcY4BMNS0mLiCSWPTvQpw+we7c9KlKyJHDoENC3L1Ctml0S4PJlp1spQULBSKgnrvbrZ78TEhFJTVSUPXrKVXevvw4UK2bvYfXgg0CNGsDHHwNxcU63UgKcgpFQxK3HuX9FtmxAr15Ot0ZEAkGOHPZmmnv32ittihQB9uyxd/iuXRv47DO7bpFIJigYCUVvvml/7NzZfpcjIpJRuXIBTzxhj46wimuhQsCOHfa+Vg0aAF9/bSfGi7hBwUio0T40IuIJefIATz9tByXDhgH58tk1Szp2BG68EVi4UEGJZJiCkVDdh6Z+faBpU6dbIyKBjkHI88/bQQkLpnHkhCXnb78daNkS+PFHp1soAUDBSKjuQ8OENO1DIyKewumaUaPsnBLWJmHi67JlwC232Bv1MUARSYOCkVCifWhExNuYhzZ+vJ3c6lqtx+eeJk3sKRxO5Ygko2AkFJfzsqCR9qEREW8qXdoeiWXxtB497LLzTG7lbsFMdmXSq8gVCkZCBedzXfvQPPSQ060RkVDB/W2mTQO2brX3wCEuA+ZuwXxjxOcmCXkKRkLF5Mnah0ZEnMNdgbk78M8/A3fdZeewMUjh8xHfIB0+7HQLJVCCkZdeeglhYWFJjur8A0vDtGnTUjw+BwvniHP70Gg5r4g4qW5d4Kuv7MKLXHHDkvLctJO7BQ8cCJw44XQLJRBGRmrVqoWjR48mHMuXL0/38fny5Uvy+AMHDmSlvZIZn39+dR+aDh2cbo2ICHDDDcD8+cDSpUDz5vaOwEx8rVABeOYZuyaShAy3NyWJjIxE8eLFM/x4joa483iKiYkxh8uZM2fMx9jYWHN4iutanrymP4p44w0Tdcb16YN4TtX4+OcNlX72B+pr31A/exDrHX3/PcK+/x7hL76I8HXrTGVX6403YA0YgMgaNdTPAfz3nNFruh2M7N69GyVLljTTLU2bNsXo0aNRtmzZNB9/7tw5lCtXDvHx8WjYsCFGjRplRlfSw2sOHTo0xfkFCxYgFwvqeNhCVgoMUvn37EHLNWsQFxmJBWXL4tLcuY61JZj72d+or31D/exhzz6L4mvXovr06ci/fz+yDx+O1nnzYveCBdjXrh3iWbtEAurv+cKFCxl6XJhlZbxe77x580xwUa1aNTPlwoDht99+w5YtW5A3b94Uj1+1apUJXurWrYvTp0/j1VdfxbJly7B161aU5rIvN0ZGypQpg5MnT5ppH09GbOz86OhoZOOmcUEook8fhH/wAeK7dEHcBx840oZQ6Gd/ob72DfWzl8XHI+x//0P40KEI37XLnLJKlkT8kCGI5wqc7NmdbmFQifXi3zNfv4sUKWJigPRev90aGWnXrl3CbQYYTZo0MaMen3/+OXqlsvsrR054uDRr1gw1atTA5MmTMXz48DS/T1RUlDmSYyd54z++t67rOM65cgkdk4MeeQThDv+MQdvPfkh97RvqZy/q2hWxd9+NTYMHo8HXXyPswAFE/N//IeK117iaAujWDYiIcLqVQSWbF/6eM3q9LC3tLVCgAKpWrYo9rLSXwUY1aNAgw4+XLNI+NCISyCIjcahVK1zesgWYONGu7sq6JCyixlU5//ufNuMLElkKRjhl8+uvv6JEiRIZenxcXBw2b96c4ceLh/ah4XJe7UMjIoGKI+WPPGJvZzFmDFCwILBtG/DPfwKNGwPffaegJJSCkSeffBJLly7F/v37sXLlStx9992IiIhAly5dzOe7d++OIdy18Yphw4aZpNO9e/diw4YNuP/++83S3t69e3v+J5GkuGTOtQ/Nld+PiEhA4zYWgwfboyPcKThPHmD9euYQAC1aaIfgUAlGDh8+bAIPJrDee++9KFy4MFavXo2iRYuazx88eNAktrqcOnUKffr0MXki7du3N4ksDGJq1qzp+Z9EknrzTfuj9qERkWDDN1nDhtk7BLNQGkdOGIhwh2AGJgxQJKC4lcA6g6V807FkyZIk98ePH28O8THtQyMioYBvhMeNAx5/HBgxAnjvPXvKhsc999gBi978BgTtTROMWFpZ+9CISKhgqQg+73En4Pvvt3PkmNxap46d7KrN+PyegpFgw9UzfHdA2odGREJJpUrARx8Bv/wC3H23ncj/4YdAtWrAww8DR4443UJJg4KRYN2HhlVx77jD6daIiPhe7drAl18Ca9bYm/GxJDlXFzJYeeop+zlS/IqCkWDzxhv2x379VBBIREIbl/1yZeHixay6aY8cjx1rb8bHwmlX9j0T5ykYCSbcYIrvBFgqOZWKuCIiIallS4A7zDOxn0Ugz54FuP9ZxYp2AuzffzvdwpCnYCQYl/N27gxcd53TrRER8R9Mam3f3l72y+lsJvdzuubJJ4EqVYApU3y+o7lcpWAkWPA/1aef2reVuCoikrrwcPsN29atdrJ/mTLAb78Bffvay4D5PMrEV/EpBSPBYtq0q/vQ3Hij060REfFvkZHAgw8Cu3cDEybYNUu4b1rXrkCDBsCcOSox70MKRoKB9qEREckcVm999FG7misLp3Gbey4NvvNO4KabWM3T6RaGBAUjwbYPDaN6ERFxD/e5efZZu0Aa97/JmRNYtQq49VagTRt7gYB4jYKRYFrOy31ocuVyujUiIoGrUCF7Z2C+wWOhNE7nLFhgLxPmLsHbtzvdwqCkYCTQMYqfO9e+zf84IiKSdSVK2G/0du4E/v3vqyXmWVDtgQeA/fudbmFQUTASLPvQsMogl6eJiIjnsBYJS8pv3gx06mTn6H3wgb00mLkmx4453cKgoGAkkGkfGhER36hVC5g1C/jpJ6B1a7smycSJdon5Z54BTp1yuoUBTcFIsOxD06GD060REQl+N9wALFwILFoENGkCXLgAjB5tj6Dw4/nzTrcwICkYCWTah0ZExBm33Wavtpk92x41+esve4SEIyWTJgGXLjndwoCiYCRQaR8aERFnMam1Y0fg55+Bjz+2R0eOHwf+7/+AatXsXJO4OKdbGRAUjAQq7UMjIuIfODLdrZu97JfPzVyJw9U2PXoA9erZoyeq5pouBSOBSPvQiIj4H45UP/SQXVb+5ZeBggXtPXDuvtvepuOHH5xuod9SMBKIpk61V9Jw/wTtQyMi4l9YfPKpp+wS86zqyvucVm/VCoiOBtaudbqFfkfBSCDvQ8MiZ9qHRkTEPxUoYO93w2quzCPJlg34/nt7Rc4//gFs2+Z0C/2GgpFA3IeG0Tb/yLUPjYiI/yteHPjvf4Fdu+w8Er6JZM2SOnXsbTz2q5qrgpFA8/779kftQyMiEljKlwemTbOruTKPhCPdvF/1SjVXrsQJUQpGAgn/cFloh+691+nWiIhIZrAuyZdf2tVcmUcSm6ia63PPAadPI9QoGAkkXMvOksN58wLXX+90a0REJCuYO8Icku+/t3cFZvXWkSOBChWAV16xq7uGCAUjgcQ1KtKihb2ttYiIBL5WrexREo6W1Khhv+kcPBioXNneDJUjJ0FOwUggca1RZxliEREJHmFhdh4J80mYR1KuHHD0qF23hAHK9On2VH2QUjASKBgZL1tm31YwIiISvNVce/QAdu60V+CwwjaXBrPCK2tLzZkTlNVcFYwEChbJ4XxikSL2cjAREQleUVF2bRIGIqxVki8f8MsvwJ13As2bAz/+iGCiYCTQ8kVuvRUI169NRCQk5MljV3Hdtw8YNAjIkQNYsQK45RagfXtg0yYEA72qBQrli4iIhK5ChewVNtz3pm9fezpn3jx76qZLF2D3bgQyBSOB4O+/gZUr7dsKRkREQlepUvYKmx077CCEZsywk1wZpPz2GwKRgpFAwEDk0iWgdGmgShWnWyMiIk6rXNleYbNxoz1dExcHTJlin+cmfX/+iUCiYCSQ8kU4KqKN8URExKV+feDbb+3VljfdZO/oPnYsULGiXUDt3DkEAgUjgUD5IiIikh7XChsu/a1b1y4pz9LyLDHPUvMxMfBnCkb8Hf+guKyXFIyIiEhaOHLeoYM9dfPJJ/boyIkT9iZ81asDH35oT+f4IQUj/o6RLqvuMVekTBmnWyMiIv4uPBzo2hXYvh14802geHFg/367mFq9esBXX/ld4TS3gpGXXnoJYWFhSY7qjLbSMXPmTPOYHDlyoE6dOpg7d25W2xy6+SIiIiIZlT27XU6ehdPGjAEKFAC2bgU6dQKaNQOWLEHAjozUqlULR48eTTiWL1+e5mNXrlyJLl26oFevXti4cSM6depkji1btmS13aFD+SIiIpIVuXLZG+/t3QsMGQLkzAmsXm0X0WzTBtiwAU5ze+vXyMhIFOeQTwZMmDABbdu2xSBWjQMwfPhwLFy4EJMmTcLbXCedhpiYGHO4nDlzxnyMjY01h6e4ruXJa3rU778jG8v/so033xywOzf6fT8HEfW1b6iffUP97IVqrkOHAv36IXz0aIS/+y7CFixAtgULcH2zZrhcoYJdr8SDMvq7czsY2b17N0qWLGmmXZo2bYrRo0ejbNmyqT521apVGDhwYJJzbdq0wezZs9P9HrzmUHZYMgsWLEAuRngexgDJH5VcsQKNmcNavjyWuJJYA5i/9nMwUl/7hvrZN9TPXtC2LXLVq4fqM2ag9LJlKLVyJVZ++SV+Z0VXD7pw4YLng5EmTZpg2rRpqFatmpmiYcDQvHlzM+2SN2/eFI8/duwYihUrluQc7/N8eoYMGZIkiOHISJkyZXD77bcjHzcL8mDExj/y6OhoZMuWDf4mnGvHGczedRfas6hNgPL3fg4m6mvfUD/7hvrZB3r1wsUNG3Bw3DjUHzQI2Zhn4kGumQ2PBiPt2rVLuF23bl0TnJQrVw6ff/65yQvxlKioKHMkxz9Gb/xBeuu6WXYluSiidWtE+GP7gqWfg5D62jfUz76hfvayhg2xo1s3VMye3eP9nNHrZWlpb4ECBVC1alXs4cY9qWBuyfHjx5Oc4/2M5pyEtEOH7I2PuESLuzOKiIgEqSwFI+fOncOvv/6KEiVKpPp55pQsci1NvYJDbjwv17B4sf2xcWMgf36nWyMiIuIfwciTTz6JpUuXYv/+/WbZ7t13342IiAizfJe6d+9u8j1cBgwYgO+++w7jxo3Djh07TJ2SdevW4ZFHHvH8TxJsVF9ERERChFs5I4cPHzaBxx9//IGiRYvi5ptvxurVq81tOnjwIMI5rXBFs2bNMH36dDz33HN45plnUKVKFbOSpnbt2p7/SYIJK+OpvoiIiIQIt4KRGTNmpPv5JalUc+vcubM5xA3MwTl82K6exyp5IiIiQUx70/gj16gIAxEv1FURERHxJwpG/JHyRUREJIQoGPE33KHXtZJGwYiIiIQABSP+hpsInjwJ5M5tL+sVEREJcgpG/DVfhIXOPFyWV0RExB8pGPE3yhcREZEQo2DEn1y+DCxdat9WMCIiIiFCwYg/Wb8eOHsWKFgQqFfP6daIiIj4hIIRf8wXufVWICLC6daIiIj4hIIRf6J8ERERCUEKRvzFxYvAihX2bQUjIiISQhSM+IvVq+2ApEQJoHp1p1sjIiLiMwpG/EXiXXrDwpxujYiIiM8oGPEXyhcREZEQpWDEH3A575o19m0FIyIiEmIUjPiD5cvtgmcVKgDlyzvdGhEREZ9SMOJP+SKtWjndEhEREZ9TMOIPlC8iIiIhTMGI0/74A9i06WrlVRERkRCjYMRp3BjPsoCaNYHixZ1ujYiIiM8pGHGa8kVERCTEKRhxmvJFREQkxCkYcdKRI8COHXbF1RYtnG6NiIiIIxSMOGnxYvtjw4ZAwYJOt0ZERMQRCkacpHwRERERBSOO4Qoa5YuIiIgoGHHMvn3AgQNAZCRw881Ot0ZERMQxCkacnqK58UYgd26nWyMiIuIYBSNOUb6IiIiIoWDEqXwRVzCifBEREQlxCkacsG0bcPw4kDMn0KSJ060RERFxlIIRJ7hGRZi4GhXldGtEREQcpWDECcoXERERSaBgxNfi4oAlS+zbyhcRERFRMOJzGzcCf/0F5M8PNGjgdGtEREQcp2DEqSkabozHgmciIiIhLkvByJgxYxAWFobHHnsszcdMmzbNPCbxkSNHDoQs5YuIiIgkkem35mvXrsXkyZNRt27daz42X7582LlzZ8J9BiQh6dIl4Mcf7dvKFxEREcn8yMi5c+fQrVs3vPPOOyhYsOA1H8/go3jx4glHsWLFEJJ++gm4cAEoWhSoVcvp1oiIiATuyEj//v3RoUMHtG7dGiNGjMhQ8FKuXDnEx8ejYcOGGDVqFGql82IcExNjDpczZ86Yj7GxsebwFNe1PHnN9IQvXIgIAPEtWyLu8mWECl/3cyhTX/uG+tk31M+B388ZvabbwciMGTOwYcMGM02TEdWqVcP7779vpnNOnz6NV199Fc2aNcPWrVtRunTpVL9m9OjRGDp0aIrzCxYsQK5cueBpCxcuhC/c9OWXKALgl6JFcWDuXIQaX/WzqK99Rf3sG+rnwO3nC5wNyIAwy+JGKRlz6NAhXH/99abBrlyRli1bon79+nj99dczHCXVqFEDXbp0wfDhwzM8MlKmTBmcPHnS5J94CtvCnyU6OhrZsmWDV50/j8jrrkMYR3e2bwcqVUKo8Gk/hzj1tW+on31D/Rz4/czX7yJFipjBiPRev90aGVm/fj1OnDhhplpc4uLisGzZMkyaNMkEEBERnIhIG3/QBg0aYM+ePWk+Jioqyhypfa03/iC9dd0k1qzhbxwoWxbZqlVjIg1CjU/6WQz1tW+on31D/Ry4/ZzR67kVjLRq1QqbN29Ocq5nz56oXr06Bg8efM1AxBW88Brt27dHSEm8S28IBiIiIiIeCUby5s2L2rVrJzmXO3duFC5cOOF89+7dUapUKZP3QcOGDcONN96IypUr46+//sLYsWNx4MAB9O7dGyFF9UVERERS5fESoAcPHkR4+NUVw6dOnUKfPn1w7Ngxswy4UaNGWLlyJWrWrImQwfLv69fbt2+91enWiIiIBFcwssS16Vsa98ePH2+OkLZ0KRAfz6VFQKlSTrdGRETEr2hvGl/ni4iIiEgSCkZ8QfkiIiIiaVIw4m3HjwNbtti3W7Z0ujUiIiJ+R8GIty1ebH+sXx8oXNjp1oiIiPgdBSPepnwRERGRdCkY8Tbli4iIiKRLwYg3HTgA/PorwMq0zZs73RoRERG/pGDEF6MiN9zA8rVOt0ZERMQvKRjxJuWLiIiIXJOCEW+xLOWLiIiIZICCEW/ZuRM4cgSIigKaNnW6NSIiIn5LwYi3uEZFbroJyJHD6daIiIj4LQUj3qJ8ERERkQxRMOIN3KHXVXlVwYiIiEi6FIx4A/ei+fNPezlv48ZOt0ZERMSvKRjxhp9/tj82aABERjrdGhEREb+mYMQbtm2zP9aq5XRLRERE/J6CEW9QMCIiIpJhCka8GYzUrOl0S0RERPyeghFP+/tvYO9e+7aCERERkWtSMOKNyqtc2luoEHDddU63RkRExO8pGPFmvkhYmNOtERER8XsKRjxN+SIiIiJuUTDiaQpGRERE3KJgxNO2brU/KhgRERHJEAUjnhQTA+zZY99WMCIiIpIhCkY8adcueyVNgQJAiRJOt0ZERCQgKBjxVr6IVtKIiIhkiIIRT1K+iIiIiNsUjHiSVtKIiIi4TcGIJ2mDPBEREbcpGPGUS5eA3bvt2xoZERERyTAFI57CJb2XLwN58wKlSjndGhERkYChYMQbyataSSMiIpJhCkY8RcmrIiIimaJgxFOUvCoiIpIpCkY8RSMjIiIivg9GxowZg7CwMDz22GPpPm7mzJmoXr06cuTIgTp16mDu3LkIKrGxwM6d9m0FIyIiIr4JRtauXYvJkyejbt266T5u5cqV6NKlC3r16oWNGzeiU6dO5tiyZQuCxq+/2gFJ7txAmTJOt0ZERCT4g5Fz586hW7dueOedd1CwYMF0HzthwgS0bdsWgwYNQo0aNTB8+HA0bNgQkyZNQlBO0YRr5ktERMQdkciE/v37o0OHDmjdujVGjBiR7mNXrVqFgQMHJjnXpk0bzJ49O82viYmJMYfLmTNnzMfY2FhzeIrrWlm9ZvjmzYgAEF+9OuI82L5g4al+lmtTX/uG+tk31M+B388ZvabbwciMGTOwYcMGM02TEceOHUOxYsWSnON9nk/L6NGjMXTo0BTnFyxYgFy5csHTFi5cmKWvb7RoEUoD2B4Whj3Blg/jR/0sGae+9g31s2+onwO3ny9cuOD5YOTQoUMYMGCAaTCTUb1lyJAhSUZTODJSpkwZ3H777ciXL59HIzb+LNHR0ciWLVumrxP5/PPmY7W770bV9u091r5g4al+lmtTX/uG+tk31M+B38+umQ2PBiPr16/HiRMnTM6HS1xcHJYtW2ZyQDi1EhHBCYurihcvjuPHjyc5x/s8n5aoqChzJMdO8sYfZJauyxLwu3aZm5FM5tV/mDR56/cnKamvfUP97Bvq58Dt54xez61sy1atWmHz5s3YtGlTwnH99debZFbeTh6IUNOmTbFo0aIk5xiB8XxQ2LePSS5AzpxA+fJOt0ZERCTguDUykjdvXtSuXTvJudy5c6Nw4cIJ57t3745SpUqZvA/itE6LFi0wbtw4k/TKnJN169ZhypQpCKqVNDVqaCWNiIhIJnj81fPgwYM4evRowv1mzZph+vTpJvioV68evvjiC7OSJnlQExQb5ImIiIhvlvYmtmTJknTvU+fOnc0RlFQGXkREJEs0r5BV2iBPREQkSxSMZEVcHLB9u31bIyMiIiKZomAkKw4cAC5e5FpkoEIFp1sjIiISkBSMeCJ5tXp1IJVlzSIiInJtCkayQsmrIiIiWaZgJCuUvCoiIpJlCkayQiMjIiIiWaZgJLPi4xWMiIiIeICCkcw6eJB7IwPZswOVKjndGhERkYClYCSzXKMi1aoBkVkuZCsiIhKyFIxklqZoREREPELBSGZpgzwRERGPUDCSWRoZERER8QgFI5lhWaoxIiIi4iEKRjLj8GHg3Dk7cbVyZadbIyIiEtAUjGSGa1SkalUgWzanWyMiIhLQFIxkhpJXRUREPEbBSGYoeVVERMRjFIxkhpJXRUREPEbBSFZW0mhkREREJMsUjLjryBHg9GkgIgKoUsXp1oiIiAQ8BSPuco2KcElvVJTTrREREQl4CkbcpXwRERERj1Iw4i7li4iIiHiUghF3KRgRERHxKAUj7q6kUcEzERERj1Iw4o7jx4FTp4DwcKBaNadbIyIiEhQUjGRmiqZSJSBHDqdbIyIiEhQUjLhD+SIiIiIep2DEHcoXERER8TgFI+7QyIiIiIjHKRhxhwqeiYiIeJyCkYz6/Xfg5EkgLEwraURERDxIwYi7oyIVKgC5cjndGhERkaChYCSjlLwqIiLiFQpGMkrJqyIiIl6hYCSjlLwqIiLifDDy1ltvoW7dusiXL585mjZtinnz5qX5+GnTpiEsLCzJkSNQK5dqZERERMQrIt15cOnSpTFmzBhUqVIFlmXhgw8+QMeOHbFx40bUSmPEgEHLzp07E+4zIAk4f/xh70tD1as73RoREZHQDUbuvPPOJPdHjhxpRktWr16dZjDC4KN48eIIaK5RkXLlgDx5nG6NiIhI6AYjicXFxWHmzJk4f/68ma5Jy7lz51CuXDnEx8ejYcOGGDVqVJqBi0tMTIw5XM6cOWM+xsbGmsNTXNe61jXDN29GBID4mjUR58HvHyoy2s+Sdepr31A/+4b6OfD7OaPXDLM43+KGzZs3m+Dj4sWLyJMnD6ZPn4727dun+thVq1Zh9+7dJs/k9OnTePXVV7Fs2TJs3brVTPmk5aWXXsLQoUNTnOf3yuVAjY/a776LSnPmYHenTtj2wAM+//4iIiKB6MKFC+jatauJAZi24bFg5NKlSzh48KC58BdffIF3330XS5cuRc0MJHYyQqpRowa6dOmC4cOHuzUyUqZMGZw8eTLdH8ZdbM/ChQsRHR2NbNmypfm4iLZtEf7DD7j8zjuwevTw2PcPFRntZ8k69bVvqJ99Q/0c+P3M1+8iRYpcMxhxe5ome/bsqFy5srndqFEjrF27FhMmTMDkyZOv+bX8IRs0aIA9e/ak+7ioqChzpPb13viDvOZ1t283HyLr1OGDPf79Q4W3fn+SkvraN9TPvqF+Dtx+zuj1slxnhLkgiUcxrpVnwmmeEiVKIGCcOgUcPWrfrlHD6daIiIgEHbdGRoYMGYJ27dqhbNmyOHv2rMnhWLJkCebPn28+3717d5QqVQqjR48294cNG4Ybb7zRjKT89ddfGDt2LA4cOIDevXsjYFwZFUGZMlyn7HRrREREQjsYOXHihAk4jh49ivz585vEVAYinGci5pKEh18dbDl16hT69OmDY8eOoWDBgmZaZ+XKlRnKL/EbKnYmIiLiP8HIe++9l+7nOUqS2Pjx480R0LRBnoiIiFdpb5pr0ciIiIiIVykYuRZtkCciIuJVCkbSw8qvhw/bt7WSRkRExCsUjGRkVKRkSaBAAadbIyIiEpQUjKRH+SIiIiJep2AkPcoXERER8ToFI+nRyIiIiIjXKRhJj4IRERERr1MwkpZz54ADB+zbWkkjIiLiNQpGrrUnTbFiQOHCTrdGREQkaCkYSYuSV0VERHxCwUhalC8iIiLiEwpG0qIN8kRERHxCwUhaNDIiIiLiEwpGUnP+PLB/v31bOSMiIiJepWAkNTt3ApYFFC0KFCnidGtERESCmoKR1GiKRkRExGcUjKRGyasiIiI+o2AkNRoZERER8RkFI6lRwTMRERGfUTCS3N9/A3v32rc1MiIiIuJ1CkZSW0kTHw8UKgRcd53TrREREQl6CkbSyxcJC3O6NSIiIkFPwUhyyhcRERHxKQUjyWkljYiIiE8pGElONUZERER8SsFIYjExwJ499m0FIyIiIj6hYCSxXbvslTT58wMlSjjdGhERkZCgYCSt5FWtpBEREfEJBSOJKXlVRETE5xSMJKbkVREREZ9TMJKYRkZERER8TsGIy6VLwO7d9m0VPBMREfEZBSMuXNJ7+TKQNy9QqpTTrREREQkZCkZSyxfRShoRERGfUTDionwRERERRygYcVEwIiIi4v/ByFtvvYW6desiX7585mjatCnmzZuX7tfMnDkT1atXR44cOVCnTh3MnTsXfkm79YqIiPh/MFK6dGmMGTMG69evx7p163DbbbehY8eO2OrKt0hm5cqV6NKlC3r16oWNGzeiU6dO5tiyZQv8ChNXd+60b2tkRERExH+DkTvvvBPt27dHlSpVULVqVYwcORJ58uTB6tWrU338hAkT0LZtWwwaNAg1atTA8OHD0bBhQ0yaNAl+t5ImNhbInRsoU8bp1oiIiISUyMx+YVxcnJmCOX/+vJmuSc2qVaswcODAJOfatGmD2bNnp3vtmJgYc7icOXPGfIyNjTWHp7iuFX9lpCa+enXzc4GHeLyfPfm7k9Spr31D/ewb6ufA7+eMXtPtYGTz5s0m+Lh48aIZFZk1axZqpjG1cezYMRQrVizJOd7n+fSMHj0aQ4cOTXF+wYIFyJUrFzxt75w5qAHgcL582OivOS1BYOHChU43IWSor31D/ewb6ufA7ecLFy54JxipVq0aNm3ahNOnT+OLL75Ajx49sHTp0jQDkswYMmRIkhEVjoyUKVMGt99+u0mc9WTExs6vciVyKxUdjRLt23vs+pK0n6Ojo5EtWzanmxPU1Ne+oX72DfVz4Peza2bD48FI9uzZUblyZXO7UaNGWLt2rckNmTx5corHFi9eHMePH09yjvd5Pj1RUVHmSI6d5I0/yIgdO+yPdeogQn/wXuOt35+kpL72DfWzb6ifA7efM3q9LNcZiY+PT5LfkRincxYtWpTkHKOvtHJMnBDG/JBdu+w7WkkjIiLic5HuTp+0a9cOZcuWxdmzZzF9+nQsWbIE8+fPN5/v3r07SpUqZXI+aMCAAWjRogXGjRuHDh06YMaMGWZJ8JQpU+Avch0/jjAGUzlzAuXLO90cERGRkONWMHLixAkTcBw9ehT58+c3BdAYiHCeiQ4ePIjw8KuDLc2aNTMBy3PPPYdnnnnGLAnmSpratWvDX+Q9dMi+UaMGkKjtIiIi4ofByHvvvZfu5zlKklznzp3N4a/yHjxo39AUjYiIiCNCfigg7+HD9g0FIyIiIo5QMOKaplEwIiIi4ojQDkbi4q6OjGiDPBEREUeEdjBy4AAiLl2CxZomFSo43RoREZGQFNLBSNi2bfaNatWAiAinmyMiIhKSQjsY2b7dfLS4rFdEREQcoWCEwYiSV0VERBwT0sEINDIiIiLiuJAORuIfegh7O3SA1aCB000REREJWW7v2htMrO7dsblIEZQpV87ppoiIiISskB4ZEREREecpGBERERFHKRgRERERRykYEREREUcpGBERERFHKRgRERERRykYEREREUcpGBERERFHKRgRERERRykYEREREUcpGBERERFHKRgRERERRykYEREREUcFxK69lmWZj2fOnPHodWNjY3HhwgVz3WzZsnn02nKV+tl31Ne+oX72DfVz4Pez63Xb9Toe0MHI2bNnzccyZco43RQRERHJxOt4/vz50/x8mHWtcMUPxMfH48iRI8ibNy/CwsI8GrExwDl06BDy5cvnsetKUupn31Ff+4b62TfUz4HfzwwxGIiULFkS4eHhgT0ywh+gdOnSXrs+O19/6N6nfvYd9bVvqJ99Q/0c2P2c3oiIixJYRURExFEKRkRERMRRIR2MREVF4cUXXzQfxXvUz76jvvYN9bNvqJ9Dp58DIoFVREREgldIj4yIiIiI8xSMiIiIiKMUjIiIiIijFIyIiIiIoxSMiIiIiKNCOhh54403UL58eeTIkQNNmjTBmjVrnG5SUBk9ejQaN25syvhfd9116NSpE3bu3Ol0s4LemDFjzLYJjz32mNNNCTq//fYb7r//fhQuXBg5c+ZEnTp1sG7dOqebFXTi4uLw/PPPo0KFCqafK1WqhOHDh19zszVJ37Jly3DnnXea0ux8jpg9e3aSz7N/X3jhBZQoUcL0e+vWrbF79274QsgGI5999hkGDhxo1lZv2LAB9erVQ5s2bXDixAmnmxY0li5div79+2P16tVYuHCh2Rny9ttvx/nz551uWtBau3YtJk+ejLp16zrdlKBz6tQp3HTTTWZX03nz5mHbtm0YN24cChYs6HTTgs7LL7+Mt956C5MmTcL27dvN/VdeeQUTJ050umkB7fz58+a1jm/EU8M+/u9//4u3334bP/30E3Lnzm1eFy9evOj9xlkh6oYbbrD69++fcD8uLs4qWbKkNXr0aEfbFcxOnDjBtzXW0qVLnW5KUDp79qxVpUoVa+HChVaLFi2sAQMGON2koDJ48GDr5ptvdroZIaFDhw7Wgw8+mOTcP/7xD6tbt26OtSnYALBmzZqVcD8+Pt4qXry4NXbs2IRzf/31lxUVFWV9+umnXm9PSI6MXLp0CevXrzdDUIk34+P9VatWOdq2YHb69GnzsVChQk43JShxFKpDhw5J/q7Fc77++mtcf/316Ny5s5l2bNCgAd555x2nmxWUmjVrhkWLFmHXrl3m/s8//4zly5ejXbt2TjctaO3btw/Hjh1L8vzBDe6YwuCL18WA2LXX006ePGnmJIsVK5bkPO/v2LHDsXYFs/j4eJPDwGHu2rVrO92coDNjxgwz3chpGvGOvXv3mqkDTu8+88wzpq8fffRRZM+eHT169HC6eUHl6aefNtvaV69eHREREeb5euTIkejWrZvTTQtax44dMx9Te110fc6bQjIYEWfetW/ZssW8uxHPOnToEAYMGGDycpiMLd4LqDkyMmrUKHOfIyP8m+b8uoIRz/r888/xySefYPr06ahVqxY2bdpk3sww8VJ9HZxCcpqmSJEiJto+fvx4kvO8X7x4ccfaFaweeeQRzJkzB4sXL0bp0qWdbk7Q4ZQjE68bNmyIyMhIczB5mIlovM13lZJ1XGFQs2bNJOdq1KiBgwcPOtamYDVo0CAzOnLfffeZFUv//ve/8fjjj5sVeuIdrtc+p14XQzIY4bBqo0aNzJxk4nc9vN+0aVNH2xZMmCPFQGTWrFn44YcfzDI98bxWrVph8+bN5t2j6+A7eA5p8zYDb8k6TjEmX5rOnIZy5co51qZgdeHCBZPHlxj/jvk8Ld7B52cGHYlfFzlVxlU1vnhdDNlpGs77criPT9o33HADXn/9dbPsqWfPnk43LaimZjjM+tVXX5laI655RyZFcQ27eAb7NnkeDpfksRaG8nM8h+/MmVjJaZp7773X1CWaMmWKOcSzWAuDOSJly5Y10zQbN27Ea6+9hgcffNDppgW0c+fOYc+ePUmSVvmGhYsK2NecChsxYgSqVKlighPWeuHUGGtEeZ0VwiZOnGiVLVvWyp49u1nqu3r1aqebFFT455XaMXXqVKebFvS0tNc7vvnmG6t27dpmuWP16tWtKVOmON2koHTmzBnz98vn5xw5clgVK1a0nn32WSsmJsbppgW0xYsXp/qc3KNHj4Tlvc8//7xVrFgx8zfeqlUra+fOnT5pWxj/8X7IIyIiIpK6kMwZEREREf+hYEREREQcpWBEREREHKVgRERERBylYEREREQcpWBEREREHKVgRERERBylYEREREQcpWBEREREHKVgRERERBylYERERETgpP8Hi9B8u3nfcxsAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "\n", + "N = 20\n", + "x = np.linspace(0, 10, N + 1)\n", + "with open(\"matrix_linear.txt\", 'r') as file:\n", + " data_str = file.read()\n", + "y = np.fromstring(data_str, sep=' ')\n", + "plt.plot(x, y, label=\"Linear element solution\", color='red')\n", + "plt.title(\"Linear element\")\n", + "plt.grid(True)\n", + "plt.legend()\n", + "plt.show()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "73525bc7", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGzCAYAAAD9pBdvAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAATgdJREFUeJzt3Qdc1VX/B/APe8kSlA0uFAE1cYJ777JhpZZW6vP0T8uR+mRPPWmWWmpm+pTa8mmYZaUtTXGAGrj3QnEAKkNQtuz7f53DSGQIeu/93fF5v16/uPd67+VwQvlwzvecY6JSqVQgIiIiUoipUp+YiIiISGAYISIiIkUxjBAREZGiGEaIiIhIUQwjREREpCiGESIiIlIUwwgREREpimGEiIiIFMUwQkRERIpiGCEyMnPnzoWJiQlSU1Pv+dwmTZrgueeegzaJzyc+LxEZD4YRIh138eJF/POf/0SzZs1gbW0NBwcHdOvWDcuXL8ft27eVbh7dw8cff4y1a9cq3QwinWaudAOIqGZ//PEHRo0aBSsrK4wbNw7BwcEoKCjA3r17MWvWLJw+fRpr1qzR2OePiYmBqSl/Z3nQMOLq6qr1ESYifcIwQqSjLl++jKeffhp+fn7YuXMnPDw8Kv5s8uTJiI2NlWFFk0QIIiLSNP7KQ6Sj3n//fWRnZ+Pzzz+vFETKtWjRAlOnTpW3r1y5IutAqpsOEI+LOpG7iZqRJ598Uk77uLi4yPfKy8u7Z81Ieno6pk+fLv9MhBVvb285alOXGpRvvvkGHTp0gI2NDRo2bCjDVkJCwj1fV1JSgg8//BBBQUFyqsrNzU1OXd26datKe4cPH46IiAh07NhRfp42bdrI+8LPP/8s74v3EO04evRolc917tw5PPHEE7J94nnifX799ddKzxH9LPr1r7/+wowZM9CoUSPY2dnh0UcfxY0bNyq1R4xeRUZGyueLq3fv3vf8eomMDcMIkY767bffZJ1IWFiYRt5fBBERPhYuXIihQ4fio48+wj/+8Y9aXyPCUY8ePbBixQoMHDhQ1q28+OKL8gf41atXa33tu+++K0OLv78/PvjgA0ybNg07duxAz549ZcCpjQgeYlqqvFbm+eefx7fffotBgwahsLCw0nPFiNGYMWMwYsQI+bWJwCJui+eLEPXMM89g3rx5shZH9IEIOuVEcOjatSvOnj2L1157DUuXLpUhY+TIkdi4cWOVdr388ss4fvw43nrrLfzf//2f/H82ZcqUij8XAUqEtYCAAHz99dfy+ve//13r10pklFREpHMyMjJU4q/nI488UqfnX758WT7/yy+/rPJn4vG33nqr4r64LR57+OGHKz3vpZdeko8fP3684jE/Pz/V+PHjK+7/5z//kc/5+eefq3yekpKSGtt35coVlZmZmerdd9+t9PjJkydV5ubmlR4Xn0983nJ79uyRn/Pbb7+t9No///yzyuPideKxqKioise2bt0qH7OxsVHFxcVVPL569Wr5+K5duyoe69evn6pNmzaqvLy8Sl9XWFiYyt/fv+Ix0c/itf3796/0dU+fPl1+nenp6RWPBQUFqXr16lVj3xCRSsWRESIdlJmZKT/a29tr7HOIupO7f8sXNm/eXONrfvrpJ7Rr105OR9xNTEHUREyPiBEIMRIhpnPKL3d3dzlSsmvXrhpfu2HDBjg6OmLAgAGVXiumWRo0aFDltYGBgQgNDa2436VLF/mxb9++8PX1rfL4pUuX5MebN2/K2hzRxqysrIrPk5aWJkdgLly4gGvXrlX6XGIk6c6vW4waFRcXIy4ursavh4iqYgErkQ4SdRyC+KGoKSIE3Kl58+Zy5YyoP6mJmNp4/PHH6/25xA9yMUhz9+csZ2FhUetrMzIy0Lhx42r/PCUlpdL9OwOHIIKM4OPjU+3j5XUnYnpHtPHNN9+UV02fy8vLq8bP5ezsXOk9iahuGEaIdDSMeHp64tSpU3V6fk2jEuK39LqqbWTjQYlREfH+W7ZsgZmZWZU/FyMctb1WBBFR81EdUTx6p+rev7bHS2eySj+PMHPmTDkSUh1RNFyf9ySiumEYIdJRYlWI2EMkOjq60rRDdcp/I7+7ELS26QIx4tC0adOK+2JkQPxArm33UzF6UteAdPfrxA9o8flatmxZ79du375dFq+K1TGaIoqFy0dp+vfvr7b31WTIIzIUrBkh0lGzZ8+WKzkmTpyI5OTkaqdMxMqS8pEUsbHW7t27q2y4VZP//ve/le6LFTLCkCFDanyNmKIRq0eqW1lS22jAY489JkcRxCqWu58n7ou6jJqIGg4xwjN//vwqf1ZUVHTPlTh1JUZfxLLb1atXIzExscqf37lktz7E/0N1tZHIUHFkhEhHiRGBdevW4amnnkLr1q0r7cAaFRUlCzvv3ANEhJZFixbJj2JvDBFMzp8/X+umag8//DAGDx4sR1/EHiBiSawoUK2JWF77448/yl1hX3jhBVlEKgo/xT4cq1atqvG14mt55513MGfOHFmTIpbKiuJc0QYRbEQhqJgeqU6vXr3k0l6xTPfYsWNySbEYvRAjO6IPRCAT+4Kogwho3bt3l3uRTJo0SY6WiCAo+kcsXRZBrL5EH33yySfy6xfTPCL0iGJaIrqD0st5iKh258+fV02aNEnVpEkTlaWlpcre3l7VrVs31YoVKyotQc3NzVVNmDBB5ejoKJ/z5JNPqlJSUmpc2nvmzBnVE088IZ/r7OysmjJliur27duVPvfdS3uFtLQ0+VwvLy/ZHm9vb/mc1NTUe34tP/30k6p79+4qOzs7eQUEBKgmT56siomJqXFpb7k1a9aoOnToIJfoijaLJbizZ89WXb9+vVJ7hw0bVuW14usVn6e65dCLFy+u9PjFixdV48aNU7m7u6ssLCzk1zl8+HDVjz/+WGVp78GDByu9ViwTvnu5cFJSkmyTaLP4My7zJarKRPznznBCREREpE2sGSEiIiJFMYwQERGRohhGiIiISFEMI0RERKQohhEiIiLSrzAiDooSR3C7uLjI3RDFevxDhw7V+PyIiAi5A+HdV1JS0oO2nYiIiIxt0zNx+JPYkrlPnz7yjAlxJoTYeKh8K+raxMTEVBz+JdR06FV1xBbV169fl5skcWtlIiIi/SB2DxEHfoqztsRBnGoJI++99548+fLLL7+seOzOsy1qI8KHk5MT7ocIInefuElERET6ISEhAd7e3uoJI2LLZ3GapdgKOjIyUh6l/dJLL8ltk+/loYceQn5+vtzOeu7cuXKEpSbieeIqV74vm9g6WoyOqEthYSF27dolR3pqO8KcHgz7WXvY19rBftYO9rP+97MYFRGDFvf62V2vHVitra3lxxkzZshAcvDgQUydOlWeSTF+/Pgap2dE3Yg4K0MEjM8++wxff/019u/fj5CQkGpfI8KKOFDrbuKcDltb27o2l4iIiBSUm5srz7zKyMioVKrxQGHE0tJShgpxSFe5V155RYYScZBUXYmDr3x9fWUoqcvISGZmppymSU1NrfWLuZ80GB4ejgEDBjB1axD7WXvY19rBftYO9rP+97P4+S1OFL9XGKnXNI2HhwcCAwMrPSZOE/3pp5/q1bjOnTtj7969Nf65lZWVvO4mOkkT35Cael+qjP2sPexr7WA/awf7WX/7ua7vV6+lvaLOQ0y73EkcUe7n51evxoljwEWwISIiIqrXyMj06dMRFhaGBQsW4Mknn8SBAwewZs0aeZWbM2eO3Ivkq6++kvc//PBDWbwSFBSEvLw8WTOyc+dObNu2Tf1fDRERERl2GOnUqRM2btwoA8fbb78tQ4YIG2PHjq14TmJiIuLj4yvuFxQU4NVXX5UBRRSftm3bFtu3b5dVu0RERET1CiPC8OHD5VWTtWvXVro/e/ZseRERERFVh2fTEBERkaIYRoiIiEhRDCNERESkKIYRIiIiUhTDCBERESmKYYSIiIj0a2kvkTHLzi9CYvptXM/IQ3puAfILS5BfXIL8wmLkF5WUXcUoLFLBysIUthZmsLE0g52VOWwtzWBjUXrb0cYC7o7WaGhrCVNTE6W/LCIiRTGMEN0l43YhTlxNx+nrmbh6KxeJ6Xm4JgJI+m1k5hWp9XNZmpmisYMVPByt4e5oU/rRwRpNG9mhpZs9PB2tYWLCsEJEho1hhIxaQVEJziVl4lhCesV16UZOra+xtzaHl5MNnG0tYW1hCitzMzkKYmVedtvcFBbmpnLUJLegCLkFxWVX6e3bBcW4mVuA1Ox8FBSX4Oqt2/ICblX5XHaWZmjhZo+WjRvA301c9ghwt4eHo40Ge4WISLsYRsjoJGXk4c9TifjzdBKOxKfLQHI3n4Y2aOvthKYudvBwsoank40MIGLkwt5aPadaFhaXICUrH0kZt5GYkSfbJT6KEZjYlGxcTs1BTkExjieky+tObg5WeMjHCe19neXHtt6OsLXkX2ci0k/814uMQsLNXGw5lYgtp5JwNL7yD3ZRv9HOx0n+UH/IxxHtvJ3g0sBK422yMDOVAUdcNYWVK6k5uJCSjfPJWfLjheQsXLyRg+TMfGw9nSwvQZSdiGmddt6OsEw3QcesfHg15JHrRKQfGEbIYKVk5mHD4avYfDJR1n/cqYOfM4YEu6NvQGM0dbXTyboMEVbEtIy4hrbxqHhcTPecupaJo/G3KqaWxIjKuaQseQFm+Pr9SLRo3ADdmrsgrIUrujZ1gaMtwwkR6SaGETI4MUlZ+HTPJfxy7BoKi1UVIwddmrpgSBt3DApyh5uDNfSVmI7p3LShvMqJKZ5jCbdw8HIath67gmu5JnKqR1z/i46TX3+wlyO6tXCVASzE1xlmXMVDRDqCYYQMgkqlwt7YVHy65zJ2n79R8XhHP2c80cEbAwLdtDL1ohSxTHiwowf6tXJF25KLCOs9AIcTMhB1MU1eIpScuJohr08iLsLJ1gJ9WjWWwaRny0ZyqoqISCkMI6TXRPHpb8evy5GQ0imK0lGQIcEemNijqSzwNEYibAwO9pCXkJyZh6iLqYiMuYFdMTeQnluIjUevyUuMkHRq4oz+rd0wMNAdvi62SjefiIwMwwjp7UiIKN6c//sZuQeIIDYVe7KjDyZ0bwqfhvyBeicxLfVoe295FRWXyFVEO84lY+fZFFkYu+/STXm988dZBHs5yBqVYW084Odip3TTicgIMIyQ3rl4Ixtzfz2NPRdS5f1G9lZ4vlsTjO3sxyLNOjA3M62oOZkzpDXi0nKw81wKws8kY9+lNFkcK673/4xBoIcDhrX1kOFEFPoSEWkCwwjpjZz8IqzYGYvP916Shali99J/9mqGl3q3kFuu0/0Rox/Pd2sqr7Ts0iXDYhm0qDU5k5gpr8VbS4PJyPaeeOQhL70uACYi3cMwQnoxJfPHyUS8+8dZuYRV6NOqEd4aEYQm/G1drUSR75guvvK6mVOAbaeTZN/fGUwWbjmHbs1dMbK9FwYHu6OBFf8ZIaIHw39FSKeJKYQ5P5+UPwzLd0Z9a3gQ+rVurJN7gxiShnaWeLqzr7xu5RRg86lEbDp6DQev3JIrl8T1xqaTsuj10fZe6OHvKqeAiIjqi2GEdJbYsn3WhhPIyi+S572I6RgxLWNtwSkZbXO2s8TYLn7yik/LlXu4iJU4l1Jz8Ovx6/IStTuPh3jjyY7eaNaogdJNJiI9wjBCOrlcd9GWc/jir8sVe4Use+ohrpDREWLp78v9/DGlbwu5b4kIJWJ59Y2sfKyKvCivzk0a4slOPhjaxp1n5hDRPfFfCdIp4pC4KeuOyKWnwj96NsOsQa3k1uikW8Q0mTjTR1yvD22NneeS8f3BBESev4EDV27KS6x6GtHOE0938pGH+XFqjYiqwzBCOiMiJgXTvz+GW7mFsLc2x9JR7TAwyF3pZlEdWJqbVmyylphxGz8dvoofDl1F/M1cfHcgXl5Bng54pqsfHm7nCTsWvRLRHfgvAimuuESFD7efx8pdsVCpxBkqDvh4TAfuBKqnPBxtMKWvv6zx2X/5Jn44lCBX5IjDCkUx8oI/zuLREC9Zf9LK3V7p5hKRDmAYIUVl5RXi/745IldmCGO7+OLN4YEsUjUApqYmCG3uIq//DA/Ej4ev4tv9cbiSlouvouPkJbahF6MlYomwlTn/nxMZK4YRUkx6bgHGf3EAx69myK3cFz7WRm6oRYa5GmdSz2Zyq36xTPubfXEIP5sslwmLy7WBJcZ09sXYrn7cUI3ICDGMkCLEyotnP98vD7cTh7p99UJntPV2UrpZpIXRku7+rvISh/etP5Ag60mSMvPw0c5YfBxxUW49Pz6sCUJ8nVjwSmQkGEZIkRUzz3y2X+5RIfam+GZCF9YOGCExAjK1vz9e6tMc204nY23UZTlKUr5viVh981xYE3k2DqdwiAwb10uS1ndUHbUqWgYRT0dr/PDPUAYRIyeWbYvAseHFMPz+cnc80cFbrs4Re5jM+OE4ui3aieXbL8hzc4jIMDGMkNaIo+pFELmWfhtNXGyx4f/CeBIsVRLs5Yglo9oh+rW+cn8ZdwdrpGYXYNn28whdtBNzfj6B2JQspZtJRGrGMEJakZANjP38IFKy8tHSrYEcEfFyslG6WaTDB/ZN7tMCe/7VBx+Nbi+nbMTOvN8dSED/D3bjuS8P4K/YVHmIIhHpP9aMkMaJ1TL/PWOG28WFaOPlKItVxeoKorpM4YhN0ka09ZD1JJ/vvYRtZ5IREXNDXgHu9pjUo5nc5VVM7RCRfmIYIY1KuJmLf3xzBLeLTdDRzwlfPN8ZDtYWSjeL9IxYVdO5aUN5XUnNwZd/XZY7vIrVWK9uOI4l22LksmFxwnAD7u5KpHf4qwRpTMbtQjy/9iBu5hTC206Fz54NYRChB9bE1Q7zHgnGvjn9MHtwK7kiKzEjD+/8cRahC3fgvT/PISUzT+lmElE9MIyQRhQWl+Clbw8jNiUbbg5WmNSqmOeRkFo52lrILef3/qsP3nu8DZo1skNWXhE+ibiI7u/twms/ncDFG9lKN5OI6oBhhNROFBW+sfEU/opNkzurrnmmPZyslG4VGSqxB8lTnXyxfXovrHm2g9wsraC4BOsPimLXSPzfN4dx8mqG0s0kolrwV1VSu1WRl/D9oQSYmgArx7RHoIcDrhxVulVkDLu7ilOexXXoyk2siryI7WdTsOVUkrx6+LvKkZSuzRpyZ1ciHcMwQmq1+WSinLMX3hoRhL4BbigsLFS6WWRkOjZpiM+aNERMUpYMJWJH1z0XUuXV3tdJhpJ+AY1lgCEi5XGahtTmaPwtTP/+mLwttvEW54sQKUns7rvsqYcQMbM3nu3qJ5f/Ho1Px6SvDmHw8t3YdPQaiopLlG4mkdGrdxi5du0annnmGbi4uMDGxgZt2rTBoUOHan1NREQEQkJCYGVlhRYtWmDt2rUP0mbS0SW84h/4/KIS+Rvnm8MDlW4SUQWfhraYPzJYFru+2Ku5XP57Pjkb074/JutKfjiYIIuuiUgPwsitW7fQrVs3WFhYYMuWLThz5gyWLl0KZ2fnGl9z+fJlDBs2DH369MGxY8cwbdo0TJw4EVu3blVH+0kH5BYUYcL/Dsptu0V9iNgx04zD36SDGttb47UhAfjrtb6YObAlnG0tcCUtF7N/OoHeiyPw9b445BUWK91MIqNTr5qR9957Dz4+Pvjyyy8rHmvatGmtr1m1apV8jggtQuvWrbF3714sW7YMgwYNut92kw6Z//sZ+VtmY3srfP5cRy7hJZ3naGOBKX398Xy3pli3Px6rd1+SZya9uekUVuy4gH/0bIYxXXxha8nvZSJtqNfftF9//VUGiFGjRiEyMhJeXl546aWXMGnSpBpfEx0djf79+1d6TLyHGCGpSX5+vrzKZWZmyo+iEFKdxZDl78UCy/sXfiZFnhciFicseSIYrrbmVfqT/aw97Ov6sTQFngv1wdMdPbHh8DWs2XMZSZn5cgO1jyNi8UJYEzzTxadKwGY/awf7Wf/7ua7vaaKqx0lT1tbW8uOMGTNkIDl48CCmTp0qRz/Gjx9f7WtatmyJ559/HnPmzKl4bPPmzXLqJjc3V9ad3G3u3LmYN29elcfXrVsHW1vbujaXNCyjAHjvuBlyikzQ16MEjzThnDvpt6IS4OANE4RfM0VafulUo525Cn08S9DDXQVrM6VbSKRfxM/5MWPGICMjAw4ODuoJI5aWlujYsSOioqIqHnvllVdkKBEjIOoKI9WNjIjpodTU1Fq/mPtJbOHh4RgwYICsg6G6E982E78+gt0X0tDa3R4b/tkFVjUcVMZ+1h72tXqIFTa/nkjExxGXEXczVz7mZGOB58P88GxXX1ibqdjPWsDvZ/3vZ/Hz29XV9Z5hpF7TNB4eHggMrLxKQtSA/PTTTzW+xt3dHcnJyZUeE/dFo6oLIoJYdSOuu4lO0sQ3pKbe15Ct/euyDCIigIiC1QY2995ilf2sPezrByO67qnOTfB4B1/8duI6VuyIxaXUHCzbEYsvouLwXKgvPIrYz9rCftbffq7r+9UrjIiVNDExMZUeO3/+PPz8/Gp8TWhoqBwJuZNIYOJx0k/nk7OwYEvpxmavD20Nfzd7pZtEpBHmZqZ4tL03Hm7nhd9PXMdHOy7g4o0cLN95ETZmZki0v4iJPZvDngdAEmlvae/06dOxb98+LFiwALGxsbKGY82aNZg8eXLFc8R0zLhx4yruv/jii7h06RJmz56Nc+fO4eOPP8YPP/wg34v0T35RMV757igKikrQu1UjjAutOYgSGQqxVP2Rh7ywbXovrBjdHi0a2eF2sYkMJeJQvpU7LyA7v0jpZhIZRxjp1KkTNm7ciO+++w7BwcGYP38+PvzwQ4wdO7biOYmJiYiPj6+4L5b1/vHHH3I0pF27dnKJ72effcZlvXpqydYYnEvKQkM7S7z/RFue8UFGF0pGtPPE71PCMN6/GM1c7ZBxuxBLtp1H9/d24r+7YpHDUEJUb/VeRD98+HB51aS63VV79+6No0d5Upq++ys2FZ/uuSxvv/d4W7mBFJGxhpIQVxXmPBOGrWdvYPmOC7h0IweLt8bgsz2X8I+ezTE+zI/7lBDVEc+moTpJzy3Aqz8cl7dHd/bFgEA3pZtEpDPTN+HTe2HZU+3Q1NUOt3IL5WGRPd/fJYMJd3QlujeGEaqT+b+fRVJmnhyWfnN4a6WbQ6RzoUQUuoZP74mlo9rBz8VWHo8gNk/rtXgXvoq+IuutiKh6DCN0T4fjbuKnI1fl7SVPtuPQM1Etq28e7+CN7TN6YdFjbeDlZIPkzHz855fT6LM4At8diOeBfETVYBihWhWXqPDWr6fl7Sc7eiPEt+ZDEYmolIWZKZ7u7IudM3th/iNBcHOwwvWMPMz5+ST6LY3Ej4evyr9bRFSKYYRq9f3BBJy6lgl7a3PMHhygdHOI9IqVuRmeDW2CyFl98ObwQLg2sET8zVzM3HAcA5dF4o8TiShhKCFiGKHai1YXby3d3Gx6/5ZwbXDvXVaJqCprCzNM6N4Uu2f3wWtDAuBkayE3T5u87giGr9iLneeS5RELRMaKYYRqtCz8vFwZ0NKtAZ7l5mZED0zUW73Yq7kMJVP7+aOBlTnOJGbihbWH8PgnUYiKTVW6iUSKYBihap1NzMTX++Lk7bkjguQcOBGph4O1BaYPaIk9s/vgnz2bwdrCFEfi0zHms/0Y8+k+HIm/pXQTibSKP2GoCjFcLIpWxVT20DbuCGvhqnSTiAySs50l5gxtjd2z+mB8qB8szEwQdTENj30chUlfHcK5pEylm0ikFQwjVMVvJxJx4PJN+dvav4dVPqWZiNSvsYM15j0SjF0ze2NUB2+YmgDhZ5IxZPkeTFt/FHFpOUo3kUijGEaoktyCIiz446y8/VLvFnKfBCLSDm9nWywe1U4eyCdGJUVN66Zj1+Vy4H9vPImkjDylm0ikEQwjVIk46EvstOrT0Ab/6NlM6eYQGaUWjRvg47Ed8NuU7ujVshGKSlT4dn+83M11weazuJVToHQTidSKYYQqXEnNwae7Sw/Ce2NYoFyOSETKaePtiP+90Bnf/6MrOvo5I7+oBGt2X5Ln3qzceYEnBJPBYBihCvN/P4OC4hL08HfFQB6ER6QzujRzwYYXQ/HFcx0R4G6PrPwiLNl2Hr0WR+B/UVdQUMQt5km/MYyQtP9SGnacS4G5qQneGhEEExMTpZtERHcQfyf7Brhh8ys9sPzph+DbUBzGly9XvvX7IAIbj3KLedJfDCMkLd9xQX58qpOPnK8mIt1kamqCRx7ykofxzR8ZjEb2Vki4eRvTvz+Oocv3YPsZ7uZK+odhhHDwyk25t4HY4+ClPi2Ubg4R1YGluSme7eqHyFm9MXtwKzhYmyMmOQsTvzqEUauicejKTaWbSFRnDCOEj8pGRZ7o4MOlvER6uMW8WIa/Z3ZfudW8lbkpDsXdwhOrojHxfwcRk5SldBOJ7olhxMgdjruFPRdSZa3IS72bK90cIrpPjrYW8hA+cULw6M4+MDM1wfazKRi8fDdm/HAMV2/lKt1EohoxjBi58lqRx0O84dPQVunmENEDcne0xsLH2mLb9J4VG6f9fOQa+i6JxNu/ncFN7lFCOohhxIgdjb+F3edvyN+gJrNWhMigNG9UunHaL5O7Iay5i1y2/8Vfl9GrbI8Ssdsyka5gGDFi5bUij7X3gq8LR0WIDFE7Hyd8O7ELvnqhM4I8HSrtUfLNvjgUFnOPElIew4iROp6Qjl0xpaMiU/pyVITI0Pco6dmykdxeXuxRIo57uJGVjzc2ncLAZbux+WQilwOTohhGjNSKnaWjIo885Ak/Fzulm0NEWtyjZMeM3pj3cBBc7CxxOTUHL317BCP/+xeiLqYq3UQyUgwjRujUtQxZZS+OKX+5r7/SzSEiBfYoGR/WBJGz+2BqP3/YWprh+NUMjPl0P5778gDOJmYq3UQyMgwjRryCRvyG1NSVoyJExqqBlTmmD2gplwOPC/WTS/wjYm5g6Ed78OoPx3Et/bbSTSQjwTBiZE5fz0D4mWSIo2e4goaIBLGl/NuPBCN8Ri8Ma+MhlwP/dOQq+iyJwILNZ5Gey+XApFkMI0ZmxY5Y+XFEW0+eQUNElYiR0v+ODcGmyd3QpWlDeRrwmt2X0PP9XVgdeRF5hcVKN5EMFMOIERHzwH+eTpKjIi9zBQ0R1eAhHyes/0dXfPlcJ7Rys0dmXhEWbjmHvksi8NNhng5M6scwYkQ+ibgoPw5t4wF/N3ulm0NEOr4cuE9AY2ye2gOLn2gLD0drXM/Iw6sbjmPYR3sQEZPC5cCkNgwjRiIlKw9bTiXK2//Xi2fQEFHdiL2IRnX0wa6ZvfGvwQGwtzbHuaQsPPflQTzz+X65Oo/oQTGMGIn1BxJQWKxCiK8Tgr0clW4OEekZawsz/F/v5tg9qw8mdm8KSzNT/BWbhuEr9mLq+qNIuMmD+Oj+MYwYgaLiEqzbHy9vi70FiIjul7OdJd4YHogdr/bCyIc85WO/HLuOfksj8c7vZ7jyhu4Lw4gREEt5kzLz4NrAEoOD3ZVuDhEZAHHK94dPt8fvL3dHtxalB/F9tveyXHmzZjdX3lD9MIwYga+i4+THpzv5wsrcTOnmEJEBEdO+30zogv+90BkB7qUrbxZsPidHSjYevYoSrryhOmAYMXAXkrMQfSlNbv0+pouv0s0hIgNdedOrZSP88UoPLBnVTq68Ebu3Tv/+OEas3Iu9F3jmDdWOYcTAfb2vdFRkQKAbPJ1slG4OERn4ypsnOnj/vfLGyhynr2fKVTfjvziAc0k884aqxzBiwLLyCuUGRcK4UBauEpF2V96Ig/ie79YEFmYmiDx/A0OX78HsH48jKSNP6SaSjmEYMWAbj15DTkExmjeyQ1hzF6WbQ0RGpqGdJd4aEYTtZWfeiPKRHw5dRe8lu7Bka4z8hYlIYBgxUGJnxPLC1We7+sk5XSIiJfi5lJ558/NLYejo54y8whKs3BWL3osj8HX0FRQWlyjdRNKnMDJ37lz5Q+3OKyAgoMbnr127tsrzra2t1dFuugdRtBqbkg1bSzM81sFb6eYQESHE1xkbXgzF6mc7oJmrHdJyCvDmL6cxaNlubD2dxO3ljZh5fV8QFBSE7du3//0G5rW/hYODA2JiYiru8zd07fi6bFTksRAvOFhbKN0cIqKKnwGDgtzRN6Ax1h+Ix4fbL+BSag7++fVhdG7SEHOGBqC9r7PSzSRdDyMifLi7u9frG68+z6cHl5hxG9vOJMvbLFwlIl1kYWaKZ0ObYGR7L6yOvIRP91zCgSs38ejHURje1gOzBwXAw4G/SBmLeoeRCxcuwNPTU063hIaGYuHChfD1rXn/iuzsbPj5+aGkpAQhISFYsGCBHF2pTX5+vrzKZWaWLgcrLCyUl7qUv5c631MXfBN9RR7x3bmJM5o2tFb86zPUftZF7GvtYD+rj7UZMLVvMzzZwRMf7ojFxmPX8fuJRDltM6ajF1oVs5/1+fu5ru9poqrHJN2WLVtkuGjVqhUSExMxb948XLt2DadOnYK9fdUj6aOjo2V4adu2LTIyMrBkyRLs3r0bp0+fhre3d621KeK977Zu3TrY2trWtblGqagEmHvEDFmFJniuZTHau3AOloj0x7Uc4Jc4U8RklJY02pipMNC7BD3dVTDnkgu9k5ubizFjxsgMIMo21BJG7paeni5HPT744ANMmDChTgmpdevWGD16NObPn1+vkREfHx+kpqbW+sXUl2hPeHg4BgwYAAsLwxgO/O1EImZsOAk3eyvserWHHApVmiH2s65iX2sH+1nz9lxIxaI/Y3A+JUfe93a2wcwB/hga7MbaQz36fhY/v11dXe8ZRuo9TXMnJycntGzZErGxsXV6vvgi27dvf8/nW1lZyau612viL76m3lcJ6w6UbnI2uosvbK2r9qGSDKmfdR37WjvYz5rTN9BD7o8076s/sTPFFldv3ca0H05gbbQT/j2sNTo1aah0Ew2OhQa+n+v6fg/0a7OYsrl48SI8PDzq9Pzi4mKcPHmyzs+n+jmfnIVDcbdgbmqCMZ15Dg0R6f/28l0bq7BtWjfMGNBSblVwLCEdo1ZF459fH8KlG9lKN5HUpF5hZObMmYiMjMSVK1cQFRWFRx99FGZmZnLaRRg3bhzmzJlT8fy3334b27Ztw6VLl3DkyBE888wziIuLw8SJE9XVfrprx1WhT0BjNHbgfi5EZBhsLc3xSj9/RMzqLQ/8FAd/bj2djIHLduOtX07hZk6B0k2kB1SvaZqrV6/K4JGWloZGjRqhe/fu2Ldvn7wtxMfHw9T073xz69YtTJo0CUlJSXB2dkaHDh1kiAkMDHzQdtNdxDHdvx67Lm+PfMhL6eYQEaldY3trLHi0DZ4Pa4JFW85hx7kU/C86Dj8fuYbJfVvgubAm8lwcMvAwsn79+lr/PCIiotL9ZcuWyYs0T0zPiCO7G1iZo1/rxko3h4hIY/zd7PH5c50QFZuKdzeflScDi3AiNnucPbgVRrT1hKkYPiG9ofxSC1KLTcdKp2gGB7vzNwMiMgphLVzx25Tu+ODJdvBwtJa/kE1dfwwjP/4L+y+lKd08qgeGEQNQUFSCzScT5W1O0RCRMREjII+FeGPnq70xa1Ar2Fma4cTVDDy1Zh8mfXUIF1nkqhcYRgxA5PkbSM8tRGN7K4Q2d1G6OUREWmdjaYbJfVogYlYfPNPVV67ECT+TLA/hE0Wuadl/711FuodhxICmaB5u5yn/AhIRGatG9lZ4Z2QbbJ3WA/0CGqOoRCWLXHsvjsCqyIvIKyxWuolUDYYRPZeVV4jtZYfiiQOniIgIaNG4tMh13cQuCPJ0QFZ+kSxy7bc0Er8cu4YH2HycNIBhRM+Jtfb5RSVo3shO/oUjIqKqRa5LR7WDu8OdRa5ROHjlptLNozIMI3puU9lGZ6Jwlec1EBFVX+T6eAdv7JrZGzMHtpRFrsfLdnJ98evDuJJaev4NKYdhRI+lZOYh6mKqvP0IV9EQEd2zyHVKX3/smtUbozuX7uT65+kkDFgWiXm/nUZ6LndyVQrDiB779fh1lKiAEF8n+LrYKt0cIiK92cl14WNt8Oe0nujdqhEKi1X48q8r6Pn+Lny25xLyi1jkqm0MI3rsl/Lt31m4SkRUby3d7LH2+c74ekJnBLjbIzOvCO/8cRYDPtiNP04ksshVixhG9JTYyOfktQy5lHdYG56CTER0v3r4N8Ifr/TAe4+3kUuD42/mYvK6I3hiVTSOxN9SunlGgWFET/1SVrja098VLg2slG4OEZFeE7/YPdXJFxEze2NqP3/YWJjhcNwtPPZxFKasO4KEm7lKN9GgMYzoITF0uIlTNEREamdnZY7pA1rKlTdPdPCGWKT4+4lE9PsgEgu3nEVmXqHSTTRIDCN66GhCuhxGtLU0w4BAN6WbQ0RkcNwdrbFkVDv8/nJ3hDV3kWeArY68JHdy/Sr6CgqLS5RuokFhGNHjKZqBgW6wtTRXujlERAYryNMR307sgs/Hd5SbS97MKcB/fjmNQR/ulrtfs8hVPRhG9IxI42LIUHiEUzRERBonNpTs19pNLgWe/0gQGtpZ4tKNHEz86hDGfLofp65lKN1Evccwomf2xqYiLacALnaW6NHCVenmEBEZDQszUzwb2gQRs3rjxV7NYWluiuhLaRixci9mbjiOpIw8pZuotxhG9MxvZYWrw9t6wNyM//uIiLTNwdoCrw0JwI4ZvTCinSfETM2Ph6+iz5IILAs/j9yCIqWbqHf400yPFBWXYMe5FHl7WFtPpZtDRGTUfBraYsXo9tj4Uhg6+DnjdmExlu+4IItcfziUgGKxRTbVCcOIHjkUdwsZtwvhbGshv/GJiEh57X2d8eOLofh4bAh8G9oiJSsfs388geEr9uKv2NLzw6h2DCN6RFRuC30CGssNeoiISHeKXIe28UD4jJ54Y1hr2Fub42xiJsZ+th8T1h5EbEq20k3UaQwjekIsHws/WxpGBrTm3iJERLrIytwME3s0Q+SsPngurAnMTU3k9LpYCvyfX07JpcFUFcOIHp1FE5eWC0szU/Ro2Ujp5hARUS3E8t+5Dwdh6/Se6N/aTdaPfBUdh16Ld2HN7os8GfguDCN6YvvZ0sLV0OYuaGDFjc6IiPRB80YN8Nn4jlg3qQsCPRyQlVeEBZvPof8HkTwZ+A4MI3pWL9K/dWOlm0JERPUU1twVv73cHYufaAs3Bysk3LxdcTLwUZ4MzDCiD9Ky83G47JtV7AJIRET6Ryw8GNXRRx7CN63/3ycDP/pxFF7+7qhRnwzMMKIHdsXckJvqBHk6wNPJRunmEBHRAxBnik3rX3oy8Kiyk4F/O35dngz83p/nkGWEJwMzjOjRFA1HRYiIDOtk4MVlJwOHNis9GfiTiIty07Rv98fJjS6NBcOIjssrLMbuCzfkbS7pJSIyzJOB103qgs/GdUSzRnby/LF/bzyFIcv3ICKmdPGCoWMY0XHiEKbcgmJZ8BTs5aB0c4iISEObpvUPdMPWaT0xd0QgnGwtcCElG899eRDjvjiAmKQsGDKGER23o2yjM7FOXXyzEhGRYZ8M/Fy3poic2QeTejSFhZkJdp+/gSHLd2POzydxIysfhohhRIeJ9efbz6RUhBEiIjIOjrYW+PewQGyf0QtDgt0hztz77kC8PBn444hYOYVvSBhGdNjp65lIysyTy7/EZmdERGRc/Fzs8MkzHfDDP0PR1tsR2flFeP/PGPRbGolfjl0zmE3TGEZ0WHjZKpqeLV1hbWGmdHOIiEghnZs2xKaXuuHDpx6Ch6M1rqXfxtT1x+QeJYfjbkLfMYzosB3nuKSXiIhKmZqaYGR7L+x8tTdmDmwJW0szHEtIx+OfRMvdXPV50zSGER2VmHEbp65lys1w+gZwC3giIiplY2mGKX39ETGzN57u5CN/TohzbsTUzcItZ5Gph5umMYzo+MF4Ib7OcG1gpXRziIhIxzR2sMaix9vij5d7oFsLFxQUl2B15CW5adrX+/Rr0zSGEZ0/GI9TNEREVLNATwd8M6ELvniudNO0mzkFeHOTfm2axjCig3LyixB9MU3e5im9RER0L2Ifqr4BpZumzXs4CM56tmkaw4gO2nPhhhxu83OxRYvGDZRuDhER6dGmaePDmiCimk3TXt+ou5umMYzooPA7NjrjrqtERKSOTdPW7dfdTdPqFUbmzp0rfzjeeQUEBNT6mg0bNsjnWFtbo02bNti8efODttmgFZeosKtsjo/1IkREpKlN0347fl1nNk2r98hIUFAQEhMTK669e/fW+NyoqCiMHj0aEyZMwNGjRzFy5Eh5nTp16kHbbbCOX02XxUcO1ubo2MRZ6eYQEZEBbZq27Kl2FZumvfzdUTz+SZTcq0Rp5vV+gbk53N3d6/Tc5cuXY/DgwZg1a5a8P3/+fISHh2PlypVYtWpVja/Lz8+XV7nMzEz5sbCwUF7qUv5e6nzPB7WnbFSka7OGQEkxCkt0ayjNUPrZULGvtYP9rB3sZ/UbHuyGfi1d8UVUHNbsuYwj8ekYteYAQlxMEdgpC00a2av189X1/129w8iFCxfg6ekpp11CQ0OxcOFC+Pr6Vvvc6OhozJgxo9JjgwYNwqZNm2r9HOI9582bV+Xxbdu2wdbWFuomApKu+P20GKwyhcPtRGzefB2GRJf62dCxr7WD/awd7Gf1awrgtTbAH/GmOHDDBEfSTLFh214EOat32iY3N1f9YaRLly5Yu3YtWrVqJadoRGDo0aOHnHaxt6+appKSkuDmVrnuQdwXj9dmzpw5lUKMGBnx8fHBwIED4eDgAHUmNvFNPmDAAFhYWEBpoqBo1sFdAEowaURPuV7cEOhaPxsy9rV2sJ+1g/2seaMBnEi4iU9+P4Cpo/rB0tJSre9fPrOh1jAyZMiQittt27aV4cTPzw8//PCDrAtRFysrK3ndTXwzauIbUlPvW18H4jJQUFQCNwcrtPRwNLiVNLrSz8aAfa0d7GftYD9rVlufhhjhVyKDiLr7ua7v90BLe52cnNCyZUvExsZW++eitiQ5uXQn0XLifl1rTozNX7Gp8mO35q4GF0SIiIg0Ekays7Nx8eJFeHh4VPvnoqZkx44dlR4TQ27icaoqqmzX1bAWrko3hYiISDfDyMyZMxEZGYkrV67IZbuPPvoozMzM5PJdYdy4cbLeo9zUqVPx559/YunSpTh37pzcp+TQoUOYMmWK+r8SPSdOWTxxtXR5VVhzF6WbQ0REpDX1qhm5evWqDB5paWlo1KgRunfvjn379snbQnx8PExN/843YWFhWLduHd544w28/vrr8Pf3lytpgoOD1f+V6Ln9l27KHfKautrB08lG6eYQERHpZhhZv359rX8eERFR5bFRo0bJi2oXdbG0XiSUoyJERGRkeDaNjoiKTasoXiUiIjImDCM6QJyiGJNcerwzR0aIiMjYMIzogOhLpaMirT0c0NBOvRvOEBER6TqGER0QVbG/CEdFiIjI+DCM6ND+It24vwgRERkhhhGFJdzMRfzNXJibmqBT04ZKN4eIiEjrGEYUFl02KtLOxwkNrOp9iDIREZHeYxhR2F9l+4tw11UiIjJWDCMKUqlUf59Hw/1FiIjISDGMKCg2JVvuMWJlbooQPyelm0NERKQIhhEF/VW2pLdTk4awMjdTujlERESKYBhRUMUUTQvWixARkfFiGFFIcYkK+8p2XuV5NEREZMwYRhRy6loGMvOKYG9tjmAvR6WbQ0REpBiGEYWnaLo2c4GZqYnSzSEiIlIMw4hCori/CBERkcQwooD8omIcvHJT3uZ5NEREZOwYRhRwND4deYUlcG1gBf/GDZRuDhERkaIYRhTw966rLjAxYb0IEREZN4YRBZQv6WW9CBEREcOI1hUWl+DE1XR5u2OThko3h4iISHEMI1p2NjFT1os42ligmaud0s0hIiJSHMOIlh2JuyU/tvd1gin3FyEiImIY0bYj8aVTNB18nZVuChERkU5gGNGyw2UjIyF+DCNEREQCw4gWpWTm4Vr6bYjZmXY+Tko3h4iISCcwjGjRkfjSUZGWbvZoYGWudHOIiIh0AsOIEvUinKIhIiKqwDCiRL0Ii1eJiIgqMIxoSUFRCU5ey5C3WbxKRET0N4YRLTl9PUMGkoZ2lmjiYqt0c4iIiHQGw4iW60Xa+zjxcDwiIqI7MIxoeedVTtEQERFVxjCi5WW9LF4lIiKqjGFECxIzbiMxIw9mpiZo5+OodHOIiIh0CsOIFhyJK60XCXC3h60lNzsjIiK6E8OIFqdouNkZERFRVQwjWsDNzoiIiGrGMKJheYXFco8RgWGEiIioKoYRDRNBpLBYBdcGlvBpaKN0c4iIiAwrjCxatEhu4DVt2rQan7N27Vr5nDsva2trGFvxqhgV4WZnREREVd330o6DBw9i9erVaNu27T2f6+DggJiYmIr7xvRDuaJehMWrRERE6gsj2dnZGDt2LD799FO8884793y+CB/u7u51fv/8/Hx5lcvMzJQfCwsL5aUu5e+lzve8k0qlqlhJ09bTXmOfR9dpup/pb+xr7WA/awf7Wf/7ua7veV9hZPLkyRg2bBj69+9fpzAiwoufnx9KSkoQEhKCBQsWICgoqMbnL1y4EPPmzavy+LZt22Brq/5D5sLDw6EJN/OBlCxzmJqocO1kNG6cgVHTVD9TVexr7WA/awf7WX/7OTc3VzNhZP369Thy5IicpqmLVq1a4YsvvpDTORkZGViyZAnCwsJw+vRpeHt7V/uaOXPmYMaMGZVGRnx8fDBw4EA55aPOxCY6f8CAAbCwsIC6/X4iEThyEsGejhg5oiuMlab7mf7GvtYO9rN2sJ/1v5/LZzbUGkYSEhIwdepU2ei6FqGGhobKq5wIIq1bt5b1JvPnz6/2NVZWVvK6m+gkTXxDaup9j1/Lkh9D/BryL5IG+5mqYl9rB/tZO9jP+tvPdX2/eoWRw4cPIyUlRU61lCsuLsbu3buxcuVKWedhZmZ2z4a1b98esbGxMJrD8Vi8SkREpJ4w0q9fP5w8ebLSY88//zwCAgLwr3/9655BpDy8iPcYOnQoDH2zszPXS4enQnydlG4OERGRYYQRe3t7BAcHV3rMzs4OLi4uFY+PGzcOXl5esghVePvtt9G1a1e0aNEC6enpWLx4MeLi4jBx4kQYshNXM1BUooKbgxW8nLjZGRERUU3UfoRsfHw8TE3/3kvt1q1bmDRpEpKSkuDs7IwOHTogKioKgYGBMJbzaIxpXxUiIiKth5GIiIha7y9btkxexqaiXoTn0RAREdWKZ9NoaLOzoxXFq6wXISIiqg3DiAYk3LyN1OwCWJqZIsjTUenmEBER6TSGEQ04ca30cLzWHvawtrj3CiMiIiJjxjCiAaeulS7pDfLiqAgREdG9MIxowOnrGfKj2AaeiIiIascwooHi1fLNzoI81XeODhERkaFiGFGzpMw8pOUUwMzUBK3c7ZVuDhERkc5jGFGz02X1Ii0aNWDxKhERUR0wjKjZ6fIpGi9O0RAREdUFw4ianSorXuX+IkRERHXDMKJmLF4lIiKqH4YRNbqVU4Br6bfl7UCGESIiojphGNFAvYifiy0crC2Ubg4REZFeYBjRwGZnnKIhIiKqO4YRTaykYfEqERFRnTGMqBFHRoiIiOqPYURNcvKLcCk1R97myAgREVHdMYyoybmkTKhUQGN7KzSyt1K6OURERHqDYUTN9SLBXhwVISIiqg+GETU5dY31IkRERPeDYUTtK2kYRoiIiOqDYUQNCopKcD45S95m8SoREVH9MIyogQgihcUqONpYwNvZRunmEBER6RWGETUejhfo4QATExOlm0NERKRXGEbUuNlZsBfrRYiIiOqLYUQNTnEbeCIiovvGMPKAiktUOJvIlTRERET3i2HkAV1Jy0FuQTGsLUzRrFEDpZtDRESkdxhG1LTZWWsPB5iZsniViIiovhhG1LSShlM0RERE94dhRG07r7J4lYiI6H4wjDwAlUqFU+XLehlGiIiI7gvDyAO4npGH9NxCmJuaoKU7i1eJiIjuB8PIAzhdVrzaonEDWJmbKd0cIiIivcQwooZ6kWAvTtEQERHdL4YRNWwDz5U0RERE949h5AFwJQ0REdGDYxi5T2nZ+UjMyJO3W3vYK90cIiIivcUw8oCjIk1d7WBvbaF0c4iIiPQWw8gDhpFA1osQEREpF0YWLVoEExMTTJs2rdbnbdiwAQEBAbC2tkabNm2wefNm6DsWrxIRESkcRg4ePIjVq1ejbdu2tT4vKioKo0ePxoQJE3D06FGMHDlSXqdOnYI+i0nKqjggj4iIiO6f+f28KDs7G2PHjsWnn36Kd955p9bnLl++HIMHD8asWbPk/fnz5yM8PBwrV67EqlWrqn1Nfn6+vMplZpZOiRQWFspLXcrfq77vWVBUgsupOfJ2cxcbtbbJEN1vP1P9sa+1g/2sHexn/e/nur7nfYWRyZMnY9iwYejfv/89w0h0dDRmzJhR6bFBgwZh06ZNNb5m4cKFmDdvXpXHt23bBltbW6ibCEf1kZgLFJWYw8pMhSN7d8LERO1NMkj17We6f+xr7WA/awf7WX/7OTc3VzNhZP369Thy5IicpqmLpKQkuLm5VXpM3BeP12TOnDmVAowYGfHx8cHAgQPh4OCg1sQmOn/AgAGwsKj7ipgtp5KA4ycQ4OGEYcO6qK09hup++5nqj32tHexn7WA/638/l89sqDWMJCQkYOrUqbLRohhVU6ysrOR1N9FJmviGrO/7Xky9LT+2dLPnX5B60NT/P6qKfa0d7GftYD/rbz/X9f3qFUYOHz6MlJQUhISEVDxWXFyM3bt3yxoQUedhZlb5wDh3d3ckJydXekzcF4/rq9iU7IowQkRERFpcTdOvXz+cPHkSx44dq7g6duwoi1nF7buDiBAaGoodO3ZUekyMrIjH9dX55NKVNP5uDZRuChERkd6r18iIvb09goODKz1mZ2cHFxeXisfHjRsHLy8vWYQqiGmdXr16YenSpbLoVdScHDp0CGvWrIE+Kiz+eyWNP0dGiIiIdG8H1vj4eCQmJlbcDwsLw7p162T4aNeuHX788Ue5kubuUKMvrqTmoKhEBTtLM3g6aq5uhoiIyFjc19LeO0VERNR6Xxg1apS8DMH55NJ6kRZu9nL3WSIiInowPJumni6klNaLtGzMehEiIiJ1YBippwtlIyMsXiUiIlIPhpH7HBlh8SoREZF6MIzc50oa7jFCRESkHgwj9VxJU1jMlTRERETqxDBSDxfKdl7lShoiIiL1YRi5n51XuZKGiIhIbRhG7mNkpCVX0hAREakNw0g9XKg4k4bFq0REROrCMHI/Z9JwmoaIiEhtGEbqKC7t75U0Xk42SjeHiIjIYDCM1PdMmsYNuJKGiIhIjRhG6r0NPOtFiIiI1IlhpI7Ol28Dz3oRIiIitWIYqaPYspERbgNPRESkXgwjdVxJcymVp/USERFpAsNIPVbS2MozabiShoiISJ0YRupTvNq4AUxNuZKGiIhInRhG6rWsl/UiRERE6sYwUgcXylbS8EwaIiIi9WMYqcc0DVfSEBERqR/DyD0U3bGSRuy+SkREROrFMHIPV9JyK1bS8EwaIiIi9WMYuYcLyVkVoyJcSUNERKR+DCP3cCGlfFkv60WIiIg0gWHkHs6XjYxwJQ0REZFmMIzcQ2z5yAjDCBERkUYwjNxrJc2NHHmb0zRERESawTByj5U0BcUlsLHgShoiIiJNYRipRWzZzqtiioYraYiIiDSDYaROZ9KwXoSIiEhTGEbqsKyX28ATERFpDsNIHTY847JeIiIizWEYqQFX0hAREWkHw0gN4m5yJQ0REZE2MIzU4GJZvUjzxnZcSUNERKRBDCM1uJJWOkXT1JX1IkRERJrEMFKDy6m58mMTF1ulm0JERGTQGEZqcCW1dGSkiYud0k0hIiIyaAwj95imaeLKMEJERKQzYeSTTz5B27Zt4eDgIK/Q0FBs2bKlxuevXbsWJiYmlS5ra2vourzCYiRm5MnbTRlGiIiINMq8Pk/29vbGokWL4O/vD5VKhf/973945JFHcPToUQQFBVX7GhFaYmJiKu6LQKLr4tJK60UcrM3hbGuhdHOIiIgMWr3CyIgRIyrdf/fdd+Voyb59+2oMIyJ8uLu7Q59cLqsXEaMi+hCeiIiIjCaM3Km4uBgbNmxATk6OnK6pSXZ2Nvz8/FBSUoKQkBAsWLCgxuBSLj8/X17lMjMz5cfCwkJ5qUv5e939nhdTSj+fj7ONWj+fsaqpn0n92NfawX7WDvaz/vdzXd/TRCXmW+rh5MmTMnzk5eWhQYMGWLduHYYOHVrtc6Ojo3HhwgVZZ5KRkYElS5Zg9+7dOH36tJzyqcncuXMxb968Ko+Lz2Vrq/mltusvmiI6xRSDvEsw1KdE45+PiIjIEOXm5mLMmDEyA4iyDbWFkYKCAsTHx8s3/vHHH/HZZ58hMjISgYGBdUpIrVu3xujRozF//vx6jYz4+PggNTW11i+mvkR7wsPDMWDAAFhY/F0bMvbzgzhw5RaWPB6MRx7yVNvnM1Y19TOpH/taO9jP2sF+1v9+Fj+/XV1d7xlG6j1NY2lpiRYtWsjbHTp0wMGDB7F8+XKsXr36nq8VX2T79u0RGxtb6/OsrKzkVd3rNfENeff7xt+8LT82d3PgXwA10tT/P6qKfa0d7GftYD/rbz/X9f0eeJ8RUQty5yjGvepMxDSPh4cHdNXtgmIkZXJZLxERkbbUa2Rkzpw5GDJkCHx9fZGVlSVrOCIiIrB161b55+PGjYOXlxcWLlwo77/99tvo2rWrHElJT0/H4sWLERcXh4kTJ0LXNztzsrWAk62l0s0hIiIyePUKIykpKTJwJCYmwtHRURamiiAi5pkEUUtiavr3YMutW7cwadIkJCUlwdnZWU7rREVF1am+ROlt4P24DTwREZHuhZHPP/+81j8XoyR3WrZsmbz0yeXy03p5QB4REZFW8Gyamg7IY70IERGRVjCM3OVK2VbwLF4lIiLSDoaRmkZGWDNCRESkFQwjd8jJL0JKVukyZU7TEBERaQfDSDXLehvaWcLRhhvsEBERaQPDyB2upJbWi/hxJQ0REZHWMIxUMzLSlPUiREREWsMwcgcu6yUiItI+hpFqRkYYRoiIiLSHYeQOl8tqRjhNQ0REpD0MI2Wy8gqRml26rNfPlQWsRERE2sIwUiaubOdVFztLOFhzWS8REZG2MIyUYb0IERGRMhhGynAbeCIiImUwjNxdvMp6ESIiIq1iGCnDaRoiIiJlMIyU4TQNERGRMhhGypb1puUUyNscGSEiItIuhhG5rPe2/OjawAoNrMyVbg4REZFRYRi584A8Fq8SERFpHcOIDCOlK2lYL0JERKR9DCN37L7KehEiIiLtYxgRIyM3y/cYYRghIiLSNoaRO0dGOE1DRESkdUYfRnKLgFu5hfK2nwsLWImIiLTN6MPIjbzSj43trWDHZb1ERERaxzBy20R+ZPEqERGRMhhGykZGmrJehIiISBEMI3kcGSEiIlKS0YeR1LIwwt1XiYiIlGH0YaR8msaP0zRERESKMOowkp5biNyismkahhEiIiJFGHUYKT8gz83BCjaWZko3h4iIyCgZdRj5e+dV1osQEREpxajDyN+n9TKMEBERKcWow0hc2QF53AaeiIhIOcYdRspGRvwaMowQEREpxajDyOjOPujlXoLWHvZKN4WIiMhoGfXJcE+EeME26Th8nDkyQkREpBSjHhkhIiIiPQsjn3zyCdq2bQsHBwd5hYaGYsuWLbW+ZsOGDQgICIC1tTXatGmDzZs3P2ibiYiIyFjDiLe3NxYtWoTDhw/j0KFD6Nu3Lx555BGcPn262udHRUVh9OjRmDBhAo4ePYqRI0fK69SpU+pqPxERERlTGBkxYgSGDh0Kf39/tGzZEu+++y4aNGiAffv2Vfv85cuXY/DgwZg1axZat26N+fPnIyQkBCtXrlRX+4mIiMhYC1iLi4vlFExOTo6crqlOdHQ0ZsyYUemxQYMGYdOmTbW+d35+vrzKZWZmyo+FhYXyUpfy91Lne1JV7GftYV9rB/tZO9jP+t/PdX3PeoeRkydPyvCRl5cnR0U2btyIwMDAap+blJQENze3So+J++Lx2ixcuBDz5s2r8vi2bdtga6v+lS/h4eFqf0+qiv2sPexr7WA/awf7WX/7OTe3dD8vtYeRVq1a4dixY8jIyMCPP/6I8ePHIzIyssZAcj/mzJlTaURFjIz4+Phg4MCBsnBWnYlNdP6AAQNgYWGhtvelytjP2sO+1g72s3awn/W/n8tnNtQeRiwtLdGiRQt5u0OHDjh48KCsDVm9enWV57q7uyM5ObnSY+K+eLw2VlZW8rqb6CRNfENq6n2pMvaz9rCvtYP9rB3sZ/3t57q+3wPvM1JSUlKpvuNOYjpnx44dlR4T6aumGhMiIiIyPub1nT4ZMmQIfH19kZWVhXXr1iEiIgJbt26Vfz5u3Dh4eXnJmg9h6tSp6NWrF5YuXYphw4Zh/fr1cknwmjVrNPPVEBERkWGHkZSUFBk4EhMT4ejoKDdAE0FEzDMJ8fHxMDX9e7AlLCxMBpY33ngDr7/+ulwSLFbSBAcHq/8rISIiIsMPI59//nmtfy5GSe42atQoeRERERFVh2fTEBERkaL04tRelUpVryVC9VnOJNZAi/dlpbbmsJ+1h32tHexn7WA/638/l//cLv85rtdhRBTLCmKvESIiItIv4ue4qDWtiYnqXnFFB4jlw9evX4e9vT1MTEzU9r7lm6klJCSodTM1qoz9rD3sa+1gP2sH+1n/+1lEDBFEPD09Ky1w0cuREfEFiBODNUV0Pr/RNY/9rD3sa+1gP2sH+1m/+7m2EZFyLGAlIiIiRTGMEBERkaKMOoyI82/eeuutas/BIfVhP2sP+1o72M/awX42nn7WiwJWIiIiMlxGPTJCREREymMYISIiIkUxjBAREZGiGEaIiIhIUQwjREREpCijDiP//e9/0aRJE1hbW6NLly44cOCA0k0yKAsXLkSnTp3kNv6NGzfGyJEjERMTo3SzDN6iRYvksQnTpk1TuikG59q1a3jmmWfg4uICGxsbtGnTBocOHVK6WQanuLgYb775Jpo2bSr7uXnz5pg/f/49D1uj2u3evRsjRoyQW7OLfyM2bdpU6c9F//7nP/+Bh4eH7Pf+/fvjwoUL0AajDSPff/89ZsyYIddWHzlyBO3atcOgQYOQkpKidNMMRmRkJCZPnox9+/YhPDxcngw5cOBA5OTkKN00g3Xw4EGsXr0abdu2VbopBufWrVvo1q2bPNV0y5YtOHPmDJYuXQpnZ2elm2Zw3nvvPXzyySdYuXIlzp49K++///77WLFihdJN02s5OTnyZ534Rbw6oo8/+ugjrFq1Cvv374ednZ38uZiXl6f5xqmMVOfOnVWTJ0+uuF9cXKzy9PRULVy4UNF2GbKUlBTxa40qMjJS6aYYpKysLJW/v78qPDxc1atXL9XUqVOVbpJB+de//qXq3r270s0wCsOGDVO98MILlR577LHHVGPHjlWsTYYGgGrjxo0V90tKSlTu7u6qxYsXVzyWnp6usrKyUn333Xcab49RjowUFBTg8OHDcgjqzsP4xP3o6GhF22bIMjIy5MeGDRsq3RSDJEahhg0bVun7mtTn119/RceOHTFq1Cg57di+fXt8+umnSjfLIIWFhWHHjh04f/68vH/8+HHs3bsXQ4YMUbppBuvy5ctISkqq9O+HOOBOlDBo4+eiXpzaq26pqalyTtLNza3S4+L+uXPnFGuXISspKZE1DGKYOzg4WOnmGJz169fL6UYxTUOacenSJTl1IKZ3X3/9ddnXr7zyCiwtLTF+/Hilm2dQXnvtNXmsfUBAAMzMzOS/1++++y7Gjh2rdNMMVlJSkvxY3c/F8j/TJKMMI6TMb+2nTp2Sv92QeiUkJGDq1KmyLkcUY5PmArUYGVmwYIG8L0ZGxPe0mF9nGFGvH374Ad9++y3WrVuHoKAgHDt2TP4yIwov2deGySinaVxdXWXaTk5OrvS4uO/u7q5YuwzVlClT8Pvvv2PXrl3w9vZWujkGR0w5isLrkJAQmJuby0sUD4tCNHFb/FZJD06sMAgMDKz0WOvWrREfH69YmwzVrFmz5OjI008/LVcsPfvss5g+fbpcoUeaUf6zT6mfi0YZRsSwaocOHeSc5J2/9Yj7oaGhirbNkIgaKRFENm7ciJ07d8pleqR+/fr1w8mTJ+Vvj+WX+A1eDGmL2yJ404MTU4x3L00XNQ1+fn6KtclQ5ebmyjq+O4nvY/HvNGmG+PdZhI47fy6KqTKxqkYbPxeNdppGzPuK4T7xj3bnzp3x4YcfymVPzz//vNJNM6ipGTHM+ssvv8i9RsrnHUVRlFjDTuoh+vbuOhyxJE/shcH6HPURv5mLwkoxTfPkk0/KfYnWrFkjL1IvsReGqBHx9fWV0zRHjx7FBx98gBdeeEHppum17OxsxMbGVipaFb+wiEUFoq/FVNg777wDf39/GU7EXi9iakzsEaVxKiO2YsUKla+vr8rS0lIu9d23b5/STTIo4turuuvLL79UumkGj0t7NeO3335TBQcHy+WOAQEBqjVr1ijdJIOUmZkpv3/Fv8/W1taqZs2aqf7973+r8vPzlW6aXtu1a1e1/yaPHz++Ynnvm2++qXJzc5Pf4/369VPFxMRopW0m4j+ajzxERERE1TPKmhEiIiLSHQwjREREpCiGESIiIlIUwwgREREpimGEiIiIFMUwQkRERIpiGCEiIiJFMYwQERGRohhGiIiISFEMI0RERKQohhEiIiKCkv4fSR2d5UiQO74AAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "N = 20\n", + "x = np.linspace(0, 10, 3 * N + 1)\n", + "with open(\"matrix_cubic.txt\", 'r') as file:\n", + " data_str = file.read()\n", + "y = np.fromstring(data_str, sep=' ')\n", + "plt.plot(x, y, label=\"Cubic element solution\")\n", + "plt.title(\"Cubic element\")\n", + "plt.grid(True)\n", + "\n", + "\n", + "plt.show()" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": ".venv", + "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 +}