1.3.viterbi_lapper.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "3813fddc-14cd-4ec2-a68d-f77246525849",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from hmmlearn import hmm\n",
    "from hmmlearn.base import ConvergenceMonitor\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib import colormaps\n",
    "import matplotlib.colors as colors\n",
    "from static_probabilities import TRANSITION_PROBABILITIES_12UL, EMISSION_PROBABILITIES_12UL, START_PROBABILITIES_12UL\n",
    "from time import time\n",
    "from collections import Counter"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "7e247c0e-9269-4449-aa1e-1fbb4e55f666",
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_matrix(matrix):\n",
    "    plt.imshow(matrix, interpolation='none', cmap=colormaps[\"autumn\"], norm=colors.PowerNorm(0.1))\n",
    "    for (j,i),label in np.ndenumerate(matrix):\n",
    "        plt.text(i,j,round(label, 3),ha='center',va='center')\n",
    "    plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "3c9fbd4e-e76b-434c-bae2-40052ee81153",
   "metadata": {},
   "outputs": [],
   "source": [
    "model = hmm.CategoricalHMM(n_components=7, n_iter=1000, init_params=\"s\")\n",
    "model.transmat_ = TRANSITION_PROBABILITIES_12UL.copy()\n",
    "model.emissionprob_ = EMISSION_PROBABILITIES_12UL.copy()\n",
    "model.startprob_ = START_PROBABILITIES_12UL.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "520d03c9-28e6-4c9b-ae08-997623e21c80",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAGdCAYAAAAv9mXmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABVhklEQVR4nO3dfVhVZb7/8beggIWQT0EoaiXqmAg+gaiTTUOSGaPXNSfRrMifjelPLaOmEVOxScWTWpTPWr/s5OVgzkTHyjTHfDiN5APoHKk85ei50BTQTFAMNeD3x71xs5HN8wL37vO6rnXZWuu71r6/7rXWd6/7XsualZaWliIiItLAPJq6ASIi4p5UYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERSzRv7A8sKSnh9OnTtGrVimbNmjX2x4uISD2UlpZy8eJFgoKC8PCo+h6l0QvM6dOnCQ4ObuyPFRGRBnTy5Ek6duxYZUyjF5hWrVoBcPIk+Pk19qdbKLOpG2CBbU3dgAa2sKkbIDUyo6kbYIGYpm5AwykohOCH7dfyqjR6gSnrFvPzc7MC49vUDbCAT1M3QH6R3PG4c8PrQ02GODTILyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpa4aQvM8uXQpQv4+EBkJOzfX3V8Sgp07w4tW0JwMDz3HBQV2devXAm9e9v/T5pRUfDpp/b158/DtGn2fXTqBM88A/n5DZjT+9Dld+AzGCKfhP1fVR2/6e/Q499MfOgY2PIPx/WXLsPUV6HjCGg5BHqOhlV/c4x5egHcPcqsb/8AjHwejv5vA+a0H7q8Dj6vQORa2H+qmpy+gh5LTXzoCtjyreP60lKY8zncsRhazoPod+G7Hxxj5u+BQW/BLfPgtuSGy6WmlgNdMP/jxUigmkOTTUAPW3wosMXKxtWRq+XkjsedO14fbsoCs3EjJCRAUhJkZkJYGMTEQF5e5fEbNsCMGSb+m2/g7bfNPmbOtMd07AgLF0JGBhw8CPffDyNHwle2L/H0aTMtXgxZWbBuHWzdChMmNFBOn0FCCiQ9BZnvQVgIxEyDvPOVx+/9J4ydBRNGwqH1MGoojHoBso7ZYxJeh63psP7P8M37MH0MTF0Em3fbY/r1gHfmmPXblpoTadhUKC5ugJyyIGEbJN0HmU9DWADErIe8S05yyoaxf4UJfeHQJBjVA0alQlauPebVf8Cb+2DVw7DvKbjVC2Leg6Jr9pirxfDIPTB5QP1zqK2NQAKQBGQCYZj/3bqTQ5O9wFhgAnAIGGWbsixuZ224Wk7ueNy54/UBoFlpaWlpbTdavnw5ixYtIicnh7CwMJYuXUpERESNti0oKMDf35/8fHMnUZnISBgwAJYtM/MlJeauZNo0U0gqmjrVFJYdO+zLnn8e9u2DL75w3pY2bWDRIudFZNMmeOwxKCyE5s2rSexg1asjn4QBPWHZi2a+pASCH4Zpo2HGkzfGxyVCYRF8/Lp92cDxEN4NViWa+V5xEPcAzH7KHtPvcRg+COZNrrwd//0dhD0Kx9Lg7o7V5PRJNTmthQFBsGxEuZxeh2kRMOPXleS0CQqvwsfjyuW0FsIDYVWsObiDlsDzUfDCYLM+vwgCFsG6UTAm1HF/6w7B9K1wIbGaPMrMrWFcFSKBAYDt0KQECAamAZUcmsQBhcDH5ZYNBMKBVfVvToO46XKaW/VqlzvuAEZUk9OTrnN9KLgE/r+B/Px8/JxdxG1qfQezceNGEhISSEpKIjMzk7CwMGJiYshzdntRS1evmruM6OhyjfQw8+nplW8zaJDZpqwb7fhx2LIFHnqo8vjiYkhNNYUjKsp5W8qKYLXFpRpXr0HGUYguV4M9PMx8+pHKt0k/AtEVfinFDHSMH9QbNu+B7/PMSbLzIHybDcMiK99n4U/wzkdwZxAEB9Qzp58h4zRE31Uhp7sg3Ul3RfpJx3iAmK72+BM/Qs4lxxh/H4js6HyfjekqkAGUOzTxsM07OTRJrxAP5u7AWXxjc7Wc3PG4c8frw/U8arvBa6+9xh/+8AfGjx9Pz549WbVqFbfccgv/7//9vwZp0LlzpgAEVEgwIABycirf5tFH4c9/hiFDoEULuPtuuO8+xy4ygCNHwNcXvL1h0iRIS4OePZ2345VXYOLEeqfEuQu2nNo4Lg9oAzk/VLoJOT9AQNuq45f+EXreZfpYvaLgwWdg+Ytwb1/H7VZsAt97zfTpXti+HLxa1DOny1BcCgG+Fdp4qzlZK83pUtXxZX/WZp+N6RxQDFQ89wIAJ4cmObWMb2yulpM7HnfueH0oU6sCc/XqVTIyMogud3vh4eFBdHQ06U5uL65cuUJBQYHD1NB27YIFC2DFCjNm88EH8MknpkCU1707HD5sus4mT4b4ePj66xv3V1AAI0aY4jN3boM3t8Es3QhfHoHNSyDjPVgyHaa8Cn/f5xg3brjpp929Grp1gtGJUHSlSZosIo3kZrg+1Krz59y5cxQXFxNQ4fYiICCAo0ePVrpNcnIyL7/8co0/o1078PSE3FzH5bm5EBhY+TazZ8Pjj8NTtr7G0FDT/TVxIrz0krndBPDygq5dzX/36wcHDsAbb8Dq1fZ9XbwIDz4IrVqZO5wWDVDJ291my6nCgF3ueQhsW+kmBLaF3B+cx/9UBDNXQNoiGDHELOsdAoe/hcXrIbrcbbC/r5lCOsHAUGh9P6TtgrEx9cjpFvBsBrkVfuHlFkKgb+XbBPpWHV/2Z+4luKOVY0y4k+++MbUDPIEKhya5gLPmBdYyvrG5Wk7ueNy54/WhjOVPkSUmJpKfn399OnnyZJXxXl7m4l9+wL6kxMw7Gy+5fNleRMp4epo/q3qEoaQErpSr1AUFMGyYacPmzeYR6Ybg1cI8rbHjgONn7zgAUaGVbxMV6hgPsH2fPf7az2byaOYY4+kBJVXkXFpqpitXa59HeV7NoV8Q7DhhX1ZSAjuOQ5STwcGoYMd4gO3/ssff2dqc7OVjCopg3ynn+2xMXkA/oNyhSYlt3tlQXlSFeIDtVcQ3NlfLyR2PO3e8PpSp1R1Mu3bt8PT0JLfC7UVubi6BTm4vvL298fb2rlWjEhJM91X//hARYd5xKSyE8ePN+ieegA4dINn2LHpsLLz2GvTpY55AO3bM3NXExtoLTWIiDB9u3m+5eNE82rxrF2zbZtaXFZfLl2H9ejNf1pvXvr19P3WV8CjEvwz9fwUR90DKX8yg2vhYW05J0KE9JE8188+OgaFPw5L15hdI6mdw8BtYYxtX8vOFoX3hj29CSx/oHAi7M+E/tsBr003M8VOwcTsMGwjtW8OpXFj4rol/aHD98gFIiIL4NOgfBBEdIOVLKLwG4/vYcvoAOvhBsq1H9dlIGLoOluyFESGQmgUHT8Ma299Bs2YwfSDM2wMhbcyJP/tzCGplHi0tk30Bzv8E2fmmP/7wGbO8axvwrd2hVvucgXigPxABpGCeqLIdmjwBdADKXpN4FhgKLME8SJSKeeBwjbXNrBVXy8kdjzt3vD5ALQuMl5cX/fr1Y8eOHYwaNQqAkpISduzYwdSpUxumRUBcHJw9C3PmmIH98HDzTkpZz1x2tuMdy6xZ5iCZNQu+/94UhNhYmD/fHpOXZwrTmTPg729euty2DR54wKzPzDRjM2DvRitz4oR56bNeOQ2DsxdgzmozEBfeDba+aR+oy85x/LUxKAw2zINZK82tbkgwfLgYepVrW+p8SFwO42bD+QJzEM2fDJN+b9b7eMN/HYaUVPixwAwC3tsH9r4Ft1cYUKxTTr3gbCHM2WkGQ8MDYetj9sHS7PwKOXWCDb+HWZ/DzB3mZP5wDPQq1+P64mDzSOnEj+BCEQzpZPbpU66rcs5OePef9vk+ti7OnfFw3531z6sqccBZYA5mUDsc2Ip90Dsbx26BQcAGYBYwEwgBPgR6WdvMWnG1nNzxuHPH6wPU4T2YjRs3Eh8fz+rVq4mIiCAlJYX333+fo0eP3jA2U5mavAfjkqp5D8YlVfMejMuZ29QNkBqZ29QNsEA178G4ktq8B1PrNzzi4uI4e/Ysc+bMIScnh/DwcLZu3Vqj4iIiIr8cdXqFcOrUqQ3aJSYiIu7npvy3yERExPWpwIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSzZvskzMB3yb79IbXv6kbIL9Ic5u6ARaY29QNkCoV1TxUdzAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQscdMWmOXvQ5ffgc9giHwS9n9Vdfymv0OPfzPxoWNgyz8c11+6DFNfhY4joOUQ6DkaVv3NMebpBXD3KLO+/QMw8nk4+r8NmNNy6NIFfHwgMhL27686PiUFuneHli0hOBieew6KiuzrV66E3r3Bz89MUVHw6af29efPw7Rp9n106gTPPAP5+Q2Ykzt+T/uhy+vg8wpEroX9p6rJ6SvosdTEh66ALd86ri8thTmfwx2LoeU8iH4XvvvBMWb+Hhj0FtwyD25Lbrhcamo50AXwASKBag5NNgE9bPGhwBYrG1dHrpaTOx53N2WB2fgZJKRA0lOQ+R6EhUDMNMg7X3n83n/C2FkwYSQcWg+jhsKoFyDrmD0m4XXYmg7r/wzfvA/Tx8DURbB5tz2mXw94Z45Zv22p+YKGTYXi4gbIaSMkJEBSEmRmQlgYxMRAXl7l8Rs2wIwZJv6bb+Dtt80+Zs60x3TsCAsXQkYGHDwI998PI0fCV7aL/OnTZlq8GLKyYN062LoVJkyofz7gpt9TFiRsg6T7IPNpCAuAmPWQd8lJTtkw9q8woS8cmgSjesCoVMjKtce8+g94cx+sehj2PQW3ekHMe1B0zR5ztRgeuQcmD6h/DrW1EUgAkoBMIAyIAZwcmuwFxgITgEPAKNuUZXE7a8PVcnLX465ZaWlpaW022LNnD4sWLSIjI4MzZ86QlpbGqFGjarx9QUEB/v7+5O8EP9/KYyKfhAE9YdmLZr6kBIIfhmmjYcaTN8bHJUJhEXz8un3ZwPEQ3g1WJZr5XnEQ9wDMfsoe0+9xGD4I5k2uvB3//R2EPQrH0uDujtUk1r/q1ZGRMGAALFtWLqdgc4cxY8aN8VOnmsKyY4d92fPPw7598MUXzj+nTRtYtMh5Edm0CR57DAoLoXnzqtvMwapXu9z39EnV+YD55TggCJaNKJfT6zAtAmb8upKcNkHhVfh4XLmc1kJ4IKyKNcUvaAk8HwUvDDbr84sgYBGsGwVjQh33t+4QTN8KFxKrbysAc2sYV4VIYABgOzQpAYKBaUAlhyZxQCHwcbllA4FwYFX9m9Mgbrqc5la92pWOu4Ii8F8I+fn5+Pn5VRlb6zuYwsJCwsLCWL58eW03rZGr1yDjKERH2Jd5eJj59COVb5N+BKIrVOCYgY7xg3rD5j3wfZ75y995EL7NhmGRle+z8Cd45yO4MwiCA+qZ01VzlxEdXSGnaEhPr3ybQYPMNmXdaMePw5Yt8NBDlccXF0NqqikcUVHO25Kfb7rTqi0u1XDL7+lnyDgN0XdVyOkuSHfSXZF+0jEeIKarPf7Ej5BzyTHG3wciOzrfZ2O6CmQA5Q5NPGzzTg5N0ivEg7k7cBbf2FwtJ3c+7mp9mRk+fDjDhw+3oi0AnLtgLpYBbRyXB7Rx3s+e8wMEtL0xPqdcf+PSP8LEBaZvv7mn+QLXvgT39nXcbsUmeHGpuXB17wzbl4NXi3rmdM6WU4ULYEAAHD1a+TaPPmq2GzLEXGh//hkmTXLsIgM4csQUlKIi8PWFtDTo2dN5O155BSZOrF8+4Kbf02UoLoWACnfWAbfC0XNOcrpUeXzOJft6qDqmKZ0DioGKtTkAcHJokuMkPqdhm1ZnrpaTOx93lo/BXLlyhYKCAoepKSzdCF8egc1LIOM9WDIdprwKf9/nGDduuBkf2L0aunWC0YlQdKXx27trFyxYACtWmDGbDz6ATz4xBaK87t3h8GHTdTZ5MsTHw9df37i/ggIYMcIUn7lzGyGBOnK170lEnKtnR0n1kpOTefnll2sc3+428PSE3AoDxbnnIbBtpZsQ2BZyf3Ae/1MRzFwBaYtgxBCzrHcIHP4WFq+H6HLdL/6+ZgrpBANDofX9kLYLxsbUOIUbc2pnyynXcXluLgQGVr7N7Nnw+OPwlG0sIjTUdH9NnAgvvWR+2QN4eUHXrua/+/WDAwfgjTdg9Wr7vi5ehAcfhFatzB1Oi3r+0gc3/Z5uAc9mkFvhF15uIQQ6GS8M9K06vuzP3EtwRyvHmHAn331jagd4AhUOTXIBZ80LrGV8Y3O1nNz5uLP8DiYxMZH8/Pzr08mTJ6uM92phnhLaccC+rKTEzEeFVr5NVKhjPMD2ffb4az+byaOZY4ynB5RU8YhDaamZrlytssnV8vIyF//yA/YlJWbe2XjJ5cv2InK9vZ72djlTUgJXyv2SLyiAYcNMGzZvNo9INwS3/J6aQ78g2HHCvqykBHYchygnDw9EBTvGA2z/lz3+ztbmZC8fU1AE+04532dj8gL6AeUOTUps886G8qIqxANsryK+sblaTu583Fl+B+Pt7Y23t3ettkl4FOJfhv6/goh7IOUvpq99fKxZ/0QSdGgPyVPN/LNjYOjTsGS9+eWb+hkc/AbW2MYr/HxhaF/445vQ0gc6B8LuTPiPLfDadBNz/BRs3A7DBkL71nAqFxa+a+IfGlz/v4eEBNN91b8/RESYd1wKC2H8eFtOT0CHDpBsexY9NhZeew369DFPoB07Zu5qYmPthSYxEYYPN++3XLxoHm3etQu2bTPry4rL5cuwfr2ZL+uhbN/evp865+SO31MUxKdB/yCI6AApX0LhNRjfx5bTB9DBD5JtI8LPRsLQdbBkL4wIgdQsOHga1tj+Dpo1g+kDYd4eCGljTvzZn0NQK/NoaZnsC3D+J8jON/3xh8+Y5V3bgG/tTp/a5wzEYx6EjABSME9U2Q5NngA6AGWvSTwLDAWWACOAVMwDh2usbWatuFpO7nrcWV5g6iJuGJy9AHNWmwHg8G6w9U37AHF2juOv3EFhsGEezFppulhCguHDxdCrqz0mdT4kLodxs+F8gbl4zZ8Mk35v1vt4w38dhpRU+LHADD7f2wf2vgW3VxjIrlNOcXD2LMyZAzk5EB5u3kkpG/jPzna8Y5k1yxwks2bB99+bghAbC/Pn22Py8kxhOnMG/P3NS5fbtsEDD5j1mZlmbAbs3WhlTpwwL33WKyd3/J56wdlCmLPTDIaGB8LWx+yDpdn5FXLqBBt+D7M+h5k7zMn84RjoVW7E+MXB5pHSiR/BhSIY0sns06dcV+WcnfDuP+3zfWxdnDvj4b47659XVeKAs8AczKB2OLAV+6B3No5dHYOADcAsYCYQAnwI9LK2mbXiajm563FX6/dgLl26xLFj5s24Pn368Nprr/Gb3/yGNm3a0KlTp2q3r8l7MC6pmvdgXFI178G4nBq8B+Ny5jZ1A6RG5jZ1AxpObd6DqfUdzMGDB/nNb35zfT4hIQGA+Ph41q1bV9vdiYiIm6p1gbnvvvuo5U2PiIj8At2U/xaZiIi4PhUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWKJ5k33yNsCnyT5daqJ/UzdAfpHmNnUDLDC3qRvQNHQHIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBI3bYFZvh+6vA4+r0DkWth/qur4TV9Bj6UmPnQFbPnWcX1pKcz5HO5YDC3nQfS78N0PjjHz98Cgt+CWeXBbcsPmA7D8fejyO/AZDJFPwv6vqo7f9Hfo8W8mPnQMbPmH4/pLl2Hqq9BxBLQcAj1Hw6q/OcY8vQDuHmXWt38ARj4PR/+3AXNaDl26gI8PREbC/v1Vx6ekQPfu0LIlBAfDc89BUZF9/cqV0Ls3+PmZKSoKPv3Uvv78eZg2zb6PTp3gmWcgP78Bc3Kz78kdz6WaWA50AXyASKCaQ5NNQA9bfCiwxcrG1ZGr5XRTFpiNWZCwDZLug8ynISwAYtZD3qXK4/dmw9i/woS+cGgSjOoBo1IhK9ce8+o/4M19sOph2PcU3OoFMe9B0TV7zNVieOQemDzAgpw+g4QUSHoKMt+DsBCImQZ5553k9E8YOwsmjIRD62HUUBj1AmQds8ckvA5b02H9n+Gb92H6GJi6CDbvtsf06wHvzDHrty01F4dhU6G4uAFy2ggJCZCUBJmZEBYGMTGQl1d5/IYNMGOGif/mG3j7bbOPmTPtMR07wsKFkJEBBw/C/ffDyJHwle0if/q0mRYvhqwsWLcOtm6FCRPqnw+43/fkjudSTWwEEoAkIBMIA2IAJ4cme4GxwATgEDDKNmVZ3M7acMWcmpWWlpbWNDg5OZkPPviAo0eP0rJlSwYNGsS///u/07179xp/YEFBAf7+/uTPAD+fymMi18KAIFg2wsyXlEDw6zAtAmb8+sb4uE1QeBU+HmdfNnAthAfCqlhzsgYtgeej4IXBZn1+EQQsgnWjYEyo4/7WHYLpW+FCYo3TghFVr458Egb0hGUvlsvpYZg2GmY8WUlOiVBYBB+/Xi6n8RDeDVbZ2tUrDuIegNlP2WP6PQ7DB8G8yZW347+/g7BH4Vga3N2xmpz6V5NTJAwYAMuWlcsp2NxhzJhxY/zUqaaw7NhhX/b887BvH3zxhfPPadMGFi1yXkQ2bYLHHoPCQmjevIoGH6w6H3DB7+mTavJxxXNpbi1inYgEBgC2Q5MSIBiYBlRyaBIHFAIfl1s2EAgHVtW/OQ3iZsmpAPAH8vPz8fPzqzK2Vncwu3fvZsqUKXz55Zds376da9euMWzYMAoLC+vRXEdXf4aM0xB9V7lGepj5dCe39uknHeMBYrra40/8CDmXHGP8fSCyo/N9NqSr1yDjKERH2Jd5eJj59COVb5N+BKIr/PqLGegYP6g3bN4D3+eZE3/nQfg2G4ZFVr7Pwp/gnY/gziAIDqhnTlfNXUZ0dIWcoiE9vfJtBg0y25R1ox0/Dlu2wEMPVR5fXAypqaZwREU5b0t+vulOq7K41IC7fU/ueC7VxFUgAyh3aOJhm3dyaJJeIR7M3YGz+MbmqjnV6pTcunWrw/y6deu4/fbbycjI4N57722QBp27DMWlEODruDzgVjh6rvJtci5VHp9zyb4eqo6x0rkL5mIZ0KbC57dx3s+e8wMEtL0xPqdcX/fSP8LEBaZvv7mnuXisfQnu7eu43YpN8OJSc+Hq3hm2LwevFvXM6ZwtpwoXwIAAOHq08m0efdRsN2SIudD+/DNMmuTYRQZw5IgpKEVF4OsLaWnQs6fzdrzyCkycWL98wP2+J3c8l2riHFAMVKzNAYCTQ5McJ/E5Ddu0OnPVnOo1BpNvG1lt06aN05grV65QUFDgMEnDWLoRvjwCm5dAxnuwZDpMeRX+vs8xbtxwMz6wezV06wSjE6HoSuO3d9cuWLAAVqwwYzYffACffGIKRHndu8Phw6brbPJkiI+Hr7++cX8FBTBihCk+c+c2QgJ15Grfk0hDqXOnQklJCdOnT2fw4MH06tXLaVxycjIvv/xyjffb7hbwbAa5FX4N5RZCoG/l2wT6Vh1f9mfuJbijlWNMeGCNm1Zn7W4DT0/IrTBQnHseAttWugmBbSH3B+fxPxXBzBWQtghGDDHLeofA4W9h8XqILtf94u9rppBOMDAUWt8PabtgbEw9cmpnyynXcXluLgQ6+TudPRsefxyeso1FhIaa7q+JE+Gll8wvewAvL+ja1fx3v35w4AC88QasXm3f18WL8OCD0KqVucNpUc87MnC/78kdz6WaaAd4AhUOTXIBZ00MrGV8Y3PVnOp8BzNlyhSysrJITU2tMi4xMZH8/Pzr08mTJ6uM92oO/YJgxwn7spIS2HEcopwMdkYFO8YDbP+XPf7O1ubEKB9TUAT7TjnfZ0PyamGeEtpxwL6spMTMR4VWvk1UqGM8wPZ99vhrP5vJo5ljjKcHlFTx2EZpqZmuXK19HuV5eZmLf/kB+5ISM+9svOTyZXsRud5eT3u7nCkpgSvlfskXFMCwYaYNmzebR6Qbgrt9T+54LtWEF9APKHdoUmKbdzaUF1UhHmB7FfGNzVVzqtMdzNSpU/n444/Zs2cPHTtWfVR5e3vj7e1dq/0nREF8GvQPgogOkPIlFF6D8X3M+ic+gA5+kGwbwXo2EoaugyV7YUQIpGbBwdOwJtasb9YMpg+EeXsgpI05SWZ/DkGtzGOYZbIvwPmfIDvf9F0fPmOWd20DvrVL4cacHoX4l6H/ryDiHkj5i+lrH29r4xNJ0KE9JE+15TQGhj4NS9abX76pn8HBb2CNbbzCzxeG9oU/vgktfaBzIOzOhP/YAq9NNzHHT8HG7TBsILRvDadyYeG7Jv6hwfXLB8wjyvHx0L8/RESYd1wKC2H8eFtOT0CHDpBsew8iNhZeew369DFPoB07Zu5qYmPthSYxEYYPN++3XLxoHm3etQu2bTPry4rL5cuwfr2ZL+t1bd/evp865+Rm35M7nks1yhuIxzwIGQGkYJ6osh2aPAF0AMpe0XkWGAoswTwQmop56HCN9U2tMVfMqVYFprS0lGnTppGWlsauXbu48847LWlUXC84WwhzdpqBw/BA2PqYfWAxO9/xF+GgTrDh9zDrc5i5wxz4H46BXuVGuF4cbB6/nPgRXCiCIZ3MPn3Kda3M2Qnv/tM+38fWJbMzHu6rZ6pxw+DsBZiz2gwAh3eDrW/aB4izcyrkFAYb5sGslaaLJSQYPlwMvbraY1LnQ+JyGDcbzheYi9f8yTDp92a9jzf812FISYUfC8zg8719YO9bcLvzYbOa5xQHZ8/CnDmQkwPh4eadlLKB/+xsxzuWWbPMBWrWLPj+e1MQYmNh/nx7TF6eKUxnzoC/v3npcts2eOABsz4z04zNgL0brcyJE+alz3rl5GbfkzueSzURB5wF5mAGtcOBrdgHvbNx7L4ZBGwAZgEzgRDgQ8B553/jc8WcavUezP/9v/+XDRs28J//+Z8O7774+/vTsmXLGu2jJu/BuKRq3oNxSdW8B+NyavAejMup5j0YlzS3qRsgVbHsPZiVK1eSn5/Pfffdxx133HF92rhxYz2aKyIi7qjWXWQiIiI1cVP+W2QiIuL6VGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEErX6XyY3qIVN9snyS9W/qRsgv1hzm7oBTUN3MCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZw6QKzHOgC+ACRwP5q4jcBPWzxocAWKxtXieX7ocvr4PMKRK6F/aeqjt/0FfRYauJDV8CWbx3Xl5bCnM/hjsXQch5Evwvf/eAYM38PDHoLbpkHtyU3bD4Ay9+HLr8Dn8EQ+STs/6rq+E1/hx7/ZuJDx8CWfziuv3QZpr4KHUdAyyHQczSs+ptjzNML4O5RZn37B2Dk83D0fxswp+XQpQv4+EBkJOyv5sBKSYHu3aFlSwgOhueeg6Ii+/qVK6F3b/DzM1NUFHz6qX39+fMwbZp9H506wTPPQH5+A+Xjjt+RG55LNeFq17xaFZiVK1fSu3dv/Pz88PPzIyoqik/LnymNaCOQACQBmUAYEAPkOYnfC4wFJgCHgFG2KcvidpbZmAUJ2yDpPsh8GsICIGY95F2qPH5vNoz9K0zoC4cmwageMCoVsnLtMa/+A97cB6sehn1Pwa1eEPMeFF2zx1wthkfugckDLMjpM0hIgaSnIPM9CAuBmGmQd95JTv+EsbNgwkg4tB5GDYVRL0DWMXtMwuuwNR3W/xm+eR+mj4Gpi2DzbntMvx7wzhyzfttSc3EYNhWKixsgp42QkABJSZCZCWFhEBMDeU4OrA0bYMYME//NN/D222YfM2faYzp2hIULISMDDh6E+++HkSPhK9uF/vRpMy1eDFlZsG4dbN0KEyY0QD7u+B254blUE652zQNoVlpaWlrT4I8++ghPT09CQkIoLS3l3XffZdGiRRw6dIh77rmnRvsoKCjA39+ffMCvrq3GVO8BwDLbfAkQDEwDZlQSHwcUAh+XWzYQCAdW1aMd182tenXkWhgQBMtGmPmSEgh+HaZFwIxf3xgftwkKr8LH48q1dy2EB8KqWHPCBi2B56PghcFmfX4RBCyCdaNgTKjj/tYdgulb4UJiLXIaUU1OT8KAnrDsxXI5PQzTRsOMJyvJKREKi+Dj18vlNB7Cu8EqW7t6xUHcAzD7KXtMv8dh+CCYN7nydvz3dxD2KBxLg7s7VtHg/lXnA+aOZcAAWLasXE7B5g5jRiUH1tSpprDs2GFf9vzzsG8ffPGF889p0wYWLXJeRDZtgsceg8JCaN68igYfrCafJ13sOwL4pJqcXPFcmluLWCdulmteAeAP5Ofn4+dX9VW8VncwsbGxPPTQQ4SEhNCtWzfmz5+Pr68vX375ZT2aW3tXgQwgutwyD9t8upNt0ivEg6n+zuIb0tWfIeM0RN9lX+bhYebTndzap590jAeI6WqPP/Ej5FxyjPH3gciOzvfZkK5eg4yjEB1hX+bhYebTj1S+TfoRiK7w6y9moGP8oN6weQ98n2dO/J0H4dtsGBZZ+T4Lf4J3PoI7gyA4oJ45XTV3GdHlDhQPDzOf7uRAGTTIbFPWjXb8OGzZAg89VHl8cTGkpprCERXlvC35+aY7rcriUg23/I7c8FyqCVe75pWp8+FbXFzMpk2bKCwsJKqKM+XKlStcuXLl+nxBQUFdP/K6c0AxUPFYDQCOOtkmx0l8Tr1bU71zl6G4FAJ8K3z+rXD0XOXb5FyqPD7nkn09VB1jpXMXzMUyoE2Fz2/jvK895wcIaHtjfE65vu6lf4SJC0z/fnNPc/FY+xLc29dxuxWb4MWl5uLVvTNsXw5eLeqZ0zlbThUOlIAAOOrkwHr0UbPdkCHmYvvzzzBpkmMXGcCRI6agFBWBry+kpUHPns7b8corMHFiPfO54IbfkRueSzXhate8MrUe5D9y5Ai+vr54e3szadIk0tLS6OnsTAGSk5Px9/e/PgUHB9erweLelm6EL4/A5iWQ8R4smQ5TXoW/73OMGzfcjBHsXg3dOsHoRCi6UukuLbVrFyxYACtWmDGbDz6ATz4xBaK87t3h8GHTdTZ5MsTHw9df37i/ggIYMcIUn7lzGyGBOnC170iaTq3vYLp3787hw4fJz8/nr3/9K/Hx8ezevdtpkUlMTCQhIeH6fEFBQb2LTDvAE8itsDwXCHSyTWAt4xtSu1vAsxnkVvg1lFsIgb6VbxPoW3V82Z+5l+COVo4x4Y2QVLvbwNMTcisMFueeh8C2lW5CYFvI/cF5/E9FMHMFpC2CEUPMst4hcPhbWLweost1wfj7mimkEwwMhdb3Q9ouGBtTj5za2XKqcKDk5kKgk7/T2bPh8cfhKdt4RGio6f6aOBFeesn8ugfw8oKuXc1/9+sHBw7AG2/A6tX2fV28CA8+CK1amTucFvX8te+W35Ebnks14WrXvDK1voPx8vKia9eu9OvXj+TkZMLCwnjjjTecxnt7e19/6qxsqi8voB9QblyVEtu8s866qArxANuriG9IXs2hXxDsOGFfVlICO45DlJMBz6hgx3iA7f+yx9/Z2pwY5WMKimDfKef7bEheLcyTQjsO2JeVlJj5qNDKt4kKdYwH2L7PHn/tZzN5NHOM8fSAkioeRSktNdOVq7XPozwvL3PxLz9gX1Ji5p31Al++bC8i19vraW+XMyUlUK7nmIICGDbMtGHzZvOIdH255XfkhudSTbjaNa9MPYYQjZKSEocxlsaSAMRjHgyKAFIwT0yMt61/AugAlD2u/iwwFFiCeTgqFfMAzprGam8UxKdB/yCI6AApX0LhNRjfx9beD6CDHyTbRuWejYSh62DJXhgRAqlZcPA0rIk165s1g+kDYd4eCGljTpLZn0NQK/MYZpnsC3D+J8jON33Xh8+Y5V3bgK93PXN6FOJfhv6/goh7IOUvpr99vK2NTyRBh/aQPNWW0xgY+jQsWW9+/aZ+Bge/gTW28Qo/XxjaF/74JrT0gc6BsDsT/mMLvDbdxBw/BRu3w7CB0L41nMqFhe+a+IcG1y8fMI8ox8dD//4QEWHecSkshPG2A+uJJ6BDB0i2HVixsfDaa9Cnj3kC7dgxc1cTG2svNImJMHy4eb/l4kXzaPOuXbBtm1lfVlwuX4b168182VBl+/b2/dQpH3f8jtzwXKpR3rjWNQ9qWWASExMZPnw4nTp14uLFi2zYsIFdu3axrexMaURxwFlgDmbQKhzYin1QKxvH27NBwAZgFjATCAE+BHo1SmshrhecLYQ5O83AYXggbH3MPrCYne/4q3BQJ9jwe5j1OczcYQ78D8dAr3Kjdi8ONo9fTvwILhTBkE5mnz7lulbm7IR3/2mf72PrktkZD/fdWc+chsHZCzBntRkEDu8GW9+0DxJn51TIKQw2zINZK003S0gwfLgYenW1x6TOh8TlMG42nC8wF7D5k2HS7816H2/4r8OQkgo/FpgB6Hv7wN634PYKg9l1yikOzp6FOXMgJwfCw807KWUD/9nZjncss2aZC9SsWfD996YgxMbC/Pn2mLw8U5jOnAF/f/PS5bZt8MADZn1mphmbAXs3WpkTJ8xLn3XOxx2/Izc8l2rC1a55UMv3YCZMmMCOHTs4c+YM/v7+9O7dmz/96U88UHam1EBDvQdz05nb1A2wQDXvwbicGrwH43KqeQ/GJVXzHoxLmtvUDWg4tXkPplZ3MG+//XY9miUiIr8kLv1vkYmIyM1LBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWaN7UDXAbc5u6AfKL1L+pGyC/OEXAwpqF6g5GREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEi5dYJYDXQAfIBLYX038JqCHLT4U2GJl4+rI1XJavh+6vA4+r0DkWth/qur4TV9Bj6UmPnQFbPnWcX1pKcz5HO5YDC3nQfS78N0PjjHz98Cgt+CWeXBbcsPmA7D8fejyO/AZDJFPwv6vqo7f9Hfo8W8mPnQMbPmH4/pLl2Hqq9BxBLQcAj1Hw6q/OcY8vQDuHmXWt38ARj4PR/+3gfJZDl26gI8PREbC/moOqpQU6N4dWraE4GB47jkoKrKvX7kSevcGPz8zRUXBp5/a158/D9Om2ffRqRM88wzk5zdMPuB+3xG457lUrwKzcOFCmjVrxvTp0xuoOTW3EUgAkoBMIAyIAfKcxO8FxgITgEPAKNuUZXE7a8PVctqYBQnbIOk+yHwawgIgZj3kXao8fm82jP0rTOgLhybBqB4wKhWycu0xr/4D3twHqx6GfU/BrV4Q8x4UXbPHXC2GR+6ByQMsyOkzSEiBpKcg8z0IC4GYaZB33klO/4Sxs2DCSDi0HkYNhVEvQNYxe0zC67A1Hdb/Gb55H6aPgamLYPNue0y/HvDOHLN+21JzcRg2FYqL65nPRkhIgKQkyMyEsDCIiYE8JwfVhg0wY4aJ/+YbePtts4+ZM+0xHTvCwoWQkQEHD8L998PIkfCV7SJ/+rSZFi+GrCxYtw62boUJE+qXy/Wc3Ow7Avc8lwCalZaWltZlwwMHDjB69Gj8/Pz4zW9+Q0pKSo22KygowN/fn3zAry4fbBMJDACW2eZLgGBgGjCjkvg4oBD4uNyygUA4sKoe7WhIN11Oc6teHbkWBgTBshFmvqQEgl+HaREw49c3xsdtgsKr8PG4cu1dC+GBsCrWnLBBS+D5KHhhsFmfXwQBi2DdKBgT6ri/dYdg+la4kFjDfEZUHxL5JAzoCcteLJfTwzBtNMx4spKcEqGwCD5+vVxO4yG8G6yytatXHMQ9ALOfssf0exyGD4J5kytvx39/B2GPwrE0uLtjFQ3uX00+kTBgACxbVi6fYHOHMaOSg2rqVFNYduywL3v+edi3D774wvnntGkDixY5LyKbNsFjj0FhITRvXnWbOVj1apf7jgA+qSYnFzqXCorAfyHk5+fj51f1VbxOdzCXLl1i3LhxrF27ltatW9dlF/VyFcgAosst87DNpzvZJr1CPJi7A2fxjc3Vcrr6M2Schui77Ms8PMx8upNb+/STjvEAMV3t8Sd+hJxLjjH+PhDZ0fk+G9LVa5BxFKIj7Ms8PMx8+pHKt0k/AtEVfv3FDHSMH9QbNu+B7/PMib/zIHybDcMiK99n4U/wzkdwZxAEB9Qjn6vmLiO63EHi4WHm050cJIMGmW3KutGOH4ctW+ChhyqPLy6G1FRTOKKinLclP990p1VbXKrhbt8RuOe5VKZOBWbKlCmMGDGC6OiKl7cbXblyhYKCAoepvs4BxUDF7zUAyHGyTU4t4xubq+V07jIUl0KAb4XPv9Uc2JXJuVR1fNmftdlnQzp3wVwwA9pU+Pw2kPNDpZuQ8wMEtK06fukfoeddpn/fKwoefAaWvwj39nXcbsUm8L3XTJ/uhe3LwatFPfI5Z8unwkESEAA5Tg6SRx+FP/8ZhgyBFi3g7rvhvvscu8gAjhwBX1/w9oZJkyAtDXr2dN6OV16BiRPrnsv1fV1wr+8I3PNcKlPrApOamkpmZibJyTUbEUpOTsbf3//6FBwcXOtGiriypRvhyyOweQlkvAdLpsOUV+Hv+xzjxg03YwS7V0O3TjA6EYquNG5bd+2CBQtgxQozZvPBB/DJJ6ZAlNe9Oxw+bLrOJk+G+Hj4+usb91dQACNGmOIzd24jJFBHrvQduZJa3bCePHmSZ599lu3bt+Pj41OjbRITE0lISLg+X1BQUO8i0w7wBHIrLM8FAp1sE1jL+Mbmajm1uwU8m0FuhV9DuYUQ6Fv5NoG+VceX/Zl7Ce5o5RgT3ghJtbsNPD0ht8Jgce55CGxb6SYEtoXcH5zH/1QEM1dA2iIYMcQs6x0Ch7+FxeshulwXjL+vmUI6wcBQaH0/pO2CsTF1zKedLZ8KB0luLgQ6+fucPRsefxyeso1FhIaa7q+JE+Gll0zXDYCXF3Ttav67Xz84cADeeANWr7bv6+JFePBBaNXK3OG0qOcvfXC/7wjc81wqU6s7mIyMDPLy8ujbty/NmzenefPm7N69mzfffJPmzZtTXMnjFN7e3vj5+TlM9eUF9APKjUNSYpt31g0cVSEeYHsV8Y3N1XLyag79gmDHCfuykhLYcRyinAx4RgU7xgNs/5c9/s7W5sQoH1NQBPtOOd9nQ/JqYZ4U2nHAvqykxMxHhVa+TVSoYzzA9n32+Gs/m8mjmWOMpweUVPF4TWmpma5crX0eZby8zMW//IB9SYmZdzZecvmyvYhcb6unvU3OlJTAlXK/5AsKYNgw04bNm80j0g3B3b4jcM9zqUyt7mB++9vfcuSI40ja+PHj6dGjB3/605/wLDsSG0ECEI95iCYCSME8UTXetv4JoANQ1pH3LDAUWIJ5mCgV87DKmkZrcfVcLaeEKIhPg/5BENEBUr6Ewmswvo+tvR9ABz9Itg3VPRsJQ9fBkr0wIgRSs+DgaVgTa9Y3awbTB8K8PRDSxpwksz+HoFbmMcwy2Rfg/E+QnW/6rg+fMcu7tgFf73rm9CjEvwz9fwUR90DKX8yA7nhbG59Igg7tIXmqLacxMPRpWLLe/PpN/QwOfgNrbGMWfr4wtC/88U1o6QOdA2F3JvzHFnhtuok5fgo2bodhA6F9aziVCwvfNfEPDa5nPgmm+6p/f4iIMO+4FBbCeNtB9cQT0KEDlPV4x8bCa69Bnz7mCbRjx8xdTWysvdAkJsLw4eb9losXzaPNu3bBtm1mfVlxuXwZ1q8382VDr+3b2/dT55zc7DsC9zyXoJYFplWrVvTq1cth2a233krbtm1vWG61OOAsMAczqB0ObMU+6J2N4+3ZIGADMAuYCYQAHwKN2+qquVpOcb3gbCHM2WkGDsMDYetj9oHF7HzHX4WDOsGG38Osz2HmDnPgfzgGepUbhH5xsHn8cuJHcKEIhnQy+/Qp170yZye8+0/7fB9bt8zOeLjvznrmNAzOXoA5q80gcHg32PqmfZA4O6dCTmGwYR7MWmm6WUKC4cPF0KurPSZ1PiQuh3Gz4XyBuYDNnwyTfm/W+3jDfx2GlFT4scAMQN/bB/a+BbdXGMyudT5xcPYszJljBvbDw807KWUD/9nZjncss2aZi9OsWfD996YgxMbC/Pn2mLw8U5jOnAF/f/PS5bZt8MADZn1mphmbAXs3WpkTJ8xLn/XKyc2+I3DPcwnq8R5Mmfvuu4/w8PBGfw9GGsHcpm5AA6vBezAup5r3YFxSNe/BuKRq3oNxJbV5D6aeT6XDrl276rsLERFxQy79b5GJiMjNSwVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIilmjeZJ88A/Bpsk9veHObugEWmNvUDZBfpP5N3QCp0iVgYc1CdQcjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEjdtgVm+H7q8Dj6vQORa2H+q6vhNX0GPpSY+dAVs+dZxfWkpzPkc7lgMLedB9Lvw3Q+OMfP3wKC34JZ5cFtyw+ZTU8uBLoAPEAnsryZ+E9DDFh8KbLGycXXkajm527G3/H3o8jvwGQyRT8L+r6qO3/R36PFvJj50DGz5h+P6S5dh6qvQcQS0HAI9R8OqvznGPL0A7h5l1rd/AEY+D0f/twFzWg5duoCPD0RGwv5qDqqUFOjeHVq2hOBgeO45KCqyr1+5Enr3Bj8/M0VFwaef2tefPw/Tptn30akTPPMM5Oc3YE5u+D3dlAVmYxYkbIOk+yDzaQgLgJj1kHep8vi92TD2rzChLxyaBKN6wKhUyMq1x7z6D3hzH6x6GPY9Bbd6Qcx7UHTNHnO1GB65ByYPsDQ9pzYCCUASkAmEATFAnpP4vcBYYAJwCBhlm7IsbmdtuFpO7nbsbfwMElIg6SnIfA/CQiBmGuSdd5LPP2HsLJgwEg6th1FDYdQLkHXMHpPwOmxNh/V/hm/eh+ljYOoi2LzbHtOvB7wzx6zfttQU2WFTobi4AXLaCAkJkJQEmZkQFgYxMZDn5KDasAFmzDDx33wDb79t9jFzpj2mY0dYuBAyMuDgQbj/fhg5Er6yXeRPnzbT4sWQlQXr1sHWrTBhQv3zAff8ngCalZaWltY0eO7cubz88ssOy7p3787Ro0dr/IEFBQX4+/uTPwP8fCqPiVwLA4Jg2QgzX1ICwa/DtAiY8esb4+M2QeFV+HicfdnAtRAeCKtizV9a0BJ4PgpeGGzW5xdBwCJYNwrGhDrub90hmL4VLiTWOC2YW4tYJyKBAcAy23wJEAxMA2ZUEh8HFAIfl1s2EAgHVtW/OQ3ipsppbvUhLnfsjagmnydhQE9Y9mK5fB6GaaNhxpOV5JMIhUXw8evl8hkP4d1gla1NveIg7gGY/ZQ9pt/jMHwQzJtceTv++zsIexSOpcHdHavJqX81OUXCgAGwbFm5nILNHcaMSg6qqVNNYdmxw77s+edh3z744gvnn9OmDSxa5LyIbNoEjz0GhYXQvHnVbeZg1atd6XsquAT+v4H8/Hz8/PyqzKvWdzD33HMPZ86cuT59UdU3VAdXf4aM0xB9l32Zh4eZT3fSVZF+0jEeIKarPf7Ej5BzyTHG3wciOzrfZ2O7CmQA0eWWedjm051sk14hHszdgbP4xuZqObnbsXf1GmQchegI+zIPDzOffqTybdKPQHSFu6iYgY7xg3rD5j3wfZ4poDsPwrfZMCyy8n0W/gTvfAR3BkFwQD1zumruMqLLHSQeHmY+3clBMmiQ2aasG+34cdiyBR56qPL44mJITTWFIyrKeVvy8013WrXFpRru+D2VqfVfTfPmzQkMDGyYT6/EuctQXAoBvo7LA26Fo+cq3ybnUuXxOZfs66HqmKZ2DigGKn6vAYCz+8McJ/E5Ddu0OnO1nNzt2Dt3wVwsA9pU+Ow2zvvZc36AgLY3xueUGzNa+keYuMD07Tf3NBfDtS/BvX0dt1uxCV5cai5c3TvD9uXg1aKeOZ2z5VThIAkIAGcdKY8+arYbMsRcaH/+GSZNcuwiAzhyxBSUoiLw9YW0NOjZ03k7XnkFJk6sXz7gnt9TmVrfwXz33XcEBQVx1113MW7cOLKzs6uMv3LlCgUFBQ6TiLiupRvhyyOweQlkvAdLpsOUV+Hv+xzjxg034wO7V0O3TjA6EYquNH57d+2CBQtgxQozZvPBB/DJJ6ZAlNe9Oxw+bLrOJk+G+Hj4+usb91dQACNGmOIzd24jJFBHN8P3VKs7mMjISNatW0f37t05c+YML7/8Mr/+9a/JysqiVatWlW6TnJx8w7hNVdrdAp7NILfCr7vcQgj0rXybQN+q48v+zL0Ed7RyjAm37masVtoBnkBuheW5gLMmBtYyvrG5Wk7uduy1uw08PSG3wkBx7nkIbFvpJgS2hdwfnMf/VAQzV0DaIhgxxCzrHQKHv4XF6yG6XPeLv6+ZQjrBwFBofT+k7YKxMfXIqZ0tpwoHSW4uOOtYmT0bHn8cnrKNRYSGmu6viRPhpZfML3sALy/o2tX8d79+cOAAvPEGrF5t39fFi/Dgg9CqlbnDadEAv/Td8XsqU6s7mOHDh/PII4/Qu3dvYmJi2LJlCxcuXOD99993uk1iYiL5+fnXp5MnT1b5GV7NoV8Q7DhhX1ZSAjuOQ5STQaeoYMd4gO3/ssff2dqc6OVjCopg3ynn+2xsXkA/oNw4JCW2eWfdwFEV4gG2VxHf2FwtJ3c79rxamKeEdhywLyspMfNRoZVvExXqGA+wfZ89/trPZvJo5hjj6QElVTwuVFpqpitXa59HeV5e5uJffsC+pMTMOxsvuXzZXkSut9fT3i5nSkrgSrlf8gUFMGyYacPmzeYR6Ybgjt9TmXoNT912221069aNY8eOOY3x9vbG29u7VvtNiIL4NOgfBBEdIOVLKLwG4/uY9U98AB38INk20PdsJAxdB0v2wogQSM2Cg6dhTaxZ36wZTB8I8/ZASBtz0s/+HIJamcdKy2RfgPM/QXa+6Ys/fMYs79oGfGuXQp0kAPGYh2gigBTME1XjbeufADoAZa9JPAsMBZZgHiZKxTysssb6ptaYq+XkbsdewqMQ/zL0/xVE3AMpfzF97eNt7XsiCTq0h+SptnzGwNCnYcl688s39TM4+A2ssY1X+PnC0L7wxzehpQ90DoTdmfAfW+C16Sbm+CnYuB2GDYT2reFULix818Q/NLjuuVzPKcF0X/XvDxER5h2XwkIYbzuonngCOnSAZNtBFRsLr70GffqYJ9COHTN3NbGx9kKTmAjDh5v3Wy5eNI8279oF27aZ9WXF5fJlWL/ezJf19rdvb99PnXNyw+8J6llgLl26xL/+9S8ef/zxhmmNTVwvOFsIc3aagdDwQNj6mH2gNDvfsTIP6gQbfg+zPoeZO8yJ/OEY6FVuIPDFweZx0okfwYUiGNLJ7NOn3C3unJ3w7j/t831st8Y74+G+Oxs0xUrFAWeBOZhB7XBgK/ZB72wcbzkHARuAWcBMIAT4EOhlfVNrzNVycrdjL24YnL0Ac1abAeDwbrD1TfsAcXZOhXzCYMM8mLXSdLGEBMOHi6FXV3tM6nxIXA7jZsP5AnPxmj8ZJv3erPfxhv86DCmp8GOBGXy+tw/sfQturzCQXaec4uDsWZgzB3JyIDzcvJNSNvCfne14xzJrlin0s2bB99+bghAbC/Pn22Py8kxhOnMG/P3NS5fbtsEDD5j1mZlmbAbs3WhlTpwwL33WKyc3/J6glu/BvPDCC8TGxtK5c2dOnz5NUlIShw8f5uuvv6Z9+/Y12kdN3oNxSXObugFSrblN3QALVPMejEuq5j0Yl1TNezCupDbvwdTqDubUqVOMHTuWH374gfbt2zNkyBC+/PLLGhcXERH55ahVgUlNTbWqHSIi4mZuyn+LTEREXJ8KjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCzRvMk+OQbwbbJPl5qY29QNaGBzm7oB8ovVv6kb0IAKah6qOxgREbGECoyIiFhCBUZERCyhAiMiIpZQgREREUuowIiIiCVUYERExBIqMCIiYgkVGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWuGkLzPL3ocvvwGcwRD4J+7+qOn7T36HHv5n40DGw5R+O6y9dhqmvQscR0HII9BwNq/7mGPP0Arh7lFnf/gEY+Twc/d8GzGk/dHkdfF6ByLWw/1Q1OX0FPZaa+NAVsOVbx/WlpTDnc7hjMbScB9Hvwnc/OMbM3wOD3oJb5sFtyQ2XS20sB7oAPkAksL+a+E1AD1t8KLDFysbVgavl447HnVteH5ZDly7g4wORkbC/mgMrJQW6d4eWLSE4GJ57DoqK7OtXroTevcHPz0xRUfDpp/b158/DtGn2fXTqBM88A/n5DZfTTVlgNn4GCSmQ9BRkvgdhIRAzDfLOVx6/958wdhZMGAmH1sOooTDqBcg6Zo9JeB22psP6P8M378P0MTB1EWzebY/p1wPemWPWb1tqTqRhU6G4uAFyyoKEbZB0H2Q+DWEBELMe8i45ySkbxv4VJvSFQ5NgVA8YlQpZufaYV/8Bb+6DVQ/DvqfgVi+IeQ+KrtljrhbDI/fA5AH1z6EuNgIJQBKQCYQBMUCek/i9wFhgAnAIGGWbsixuZ025Wj7ueNy55fVhIyQkQFISZGZCWBjExECekwNrwwaYMcPEf/MNvP222cfMmfaYjh1h4ULIyICDB+H++2HkSPjKVoxPnzbT4sWQlQXr1sHWrTBhQv3zKdOstLS0tDYbfP/99/zpT3/i008/5fLly3Tt2pV33nmH/v3712j7goIC/P39yd8Jfr6Vx0Q+CQN6wrIXzXxJCQQ/DNNGw4wnb4yPS4TCIvj4dfuygeMhvBusSjTzveIg7gGY/ZQ9pt/jMHwQzJtceTv++zsIexSOpcHdHatJ7JOqV0euhQFBsGxEuZxeh2kRMOPXleS0CQqvwsfjyuW0FsIDYVWsObiDlsDzUfDCYLM+vwgCFsG6UTAm1HF/6w7B9K1wIbGaPMqbW4tYJyKBAcAy23wJEAxMA2ZUEh8HFAIfl1s2EAgHVtW/OfV20+Uzt+rVLnncjagmpydd8PpQzeUxMhIGDIBly8rlFGzuMGZUcmBNnWoKy44d9mXPPw/79sEXXzj/nDZtYNEi50Vk0yZ47DEoLITmzSuPKSgAf3/Iz8/Hz8+vyrxqdQfz448/MnjwYFq0aMGnn37K119/zZIlS2jdunVtdlOlq9cg4yhER5RrpIeZTz9S+TbpRyC6wi+lmIGO8YN6w+Y98H2eOUl2HoRvs2FYZOX7LPwJ3vkI7gyC4IB65vQzZJyG6Lsq5HQXpDvprkg/6RgPENPVHn/iR8i55Bjj7wORHZ3vs7FdBTKA6HLLPGzz6U62Sa8QD+YOwVl8Y3K1fNzxuHPL68NVc5cRXe5A8fAw8+lODpRBg8w2Zd1ox4/Dli3w0EOVxxcXQ2qqKRxRUc7bkp9vutOcFZfaqtVu/v3f/53g4GDeeeed68vuvPPOhmmJzbkL5i8joI3j8oA2zvs7c36AgLY3xueU6xde+keYuMD0sTb3NF/g2pfg3r6O263YBC8uNQdQ986wfTl4tahnTpehuBQCKtyxBdwKR885yelS5fE5l+zroeqYpnYOKAYqnn8BwFEn2+Q4ic9p2KbViavl447HnVteH87ZcqpwoAQEwFEnB9ajj5rthgwxBfHnn2HSJMcuMoAjR0xBKSoCX19IS4OePZ2345VXYOLE+uVTXq3uYDZv3kz//v155JFHuP322+nTpw9r166tcpsrV65QUFDgMDWFpRvhyyOweQlkvAdLpsOUV+Hv+xzjxg03/bS7V0O3TjA6EYquNEmTRaSRuNr1YdcuWLAAVqwwYzYffACffGIKRHndu8Phw6brbPJkiI+Hr7++cX8FBTBihCk+c+c2XDtrdQdz/PhxVq5cSUJCAjNnzuTAgQM888wzeHl5ER8fX+k2ycnJvPzyyzX+jHa3gacn5FYYsMs9D4FtK92EwLaQ+4Pz+J+KYOYKSFsEI4aYZb1D4PC3sHg9RJe7Dfb3NVNIJxgYCq3vh7RdMDamxincmNMt4NkMciv8wssthEAn41CBvlXHl/2ZewnuaOUYEx5Y97Y2pHaAJ5BbYXku4KyJgbWMb0yulo87HndueX1oZ8upwoGSmwuBTv5OZ8+Gxx+Hp2xjRqGhpvtr4kR46SVzBwbg5QVdu5r/7tcPDhyAN96A1avt+7p4ER58EFq1Mnc4Lep5R1Zere5gSkpK6Nu3LwsWLKBPnz5MnDiRP/zhD6xa5Xy4MjExkfz8/OvTyZMnq/wMrxbmaY0dB8p/rpmPCq18m6hQx3iA7fvs8dd+NpNHM8cYTw8oqeIRh9JSM125WmWTq+XVHPoFwY4T9mUlJbDjOEQ5GRyMCnaMB9j+L3v8na3NyV4+pqAI9p1yvs/G5gX0A8qNQ1Jim3fWDRxVIR5gexXxjcnV8nHH484trw9e5uJffsC+pMTMOxsvuXzZXkSut9fT3i5nSkrgSrk7roICGDbMtGHzZvOIdEOq1R3MHXfcQc8KHXi/+tWv+Nvf/uZkC/D29sbb27tWjUp4FOJfhv6/goh7IOUvps9zfKxZ/0QSdGgPyVPN/LNjYOjTsGS9+QWS+hkc/AbW2Poj/XxhaF/445vQ0gc6B8LuTPiPLfDadBNz/BRs3A7DBkL71nAqFxa+a+IfGlyr5leeUxTEp0H/IIjoAClfQuE1GN/HltMH0MEPkm0Dfc9GwtB1sGQvjAiB1Cw4eBrW2P4OmjWD6QNh3h4IaWNO/NmfQ1Ar82hpmewLcP4nyM43/fGHz5jlXduAb+2+lrrlDcRjHqKJAFIwT1WNt61/AugAlL0q8SwwFFiCeZgoFTgIrLG+qTXiavm443HnlteHBNN91b8/RESYd1wKC2G87cB64gno0AGSbQdWbCy89hr06WOeQDt2zNzVxMbaC01iIgwfbt5vuXjRPNq8axds22bWlxWXy5dh/XozXzaC0b69fT/1UasCM3jwYP7nf/7HYdm3335L586d69+ScuKGwdkLMGe1GYgL7wZb37QP1GXnOP7aGBQGG+bBrJXmVjckGD5cDL262mNS50Pichg3G84XmINo/mSY9Huz3scb/uswpKTCjwVmEPDePrD3Lbi9woBinXLqBWcLYc5OMxgaHghbH7MPlmbnV8ipE2z4Pcz6HGbuMCfzh2OgV7mBwBcHm0dKJ34EF4pgSCezT59yt7hzdsK7/7TP97HdGu+Mh/sa9vmMSsUBZ4E5mIHtcGAr9oHvbBxvowcBG4BZwEwgBPgQ6GV9U2vE1fJxx+POLa8PcXD2LMyZAzk5EB5u3kkpG/jPzna8Y5k1yxT7WbPg++9NQYiNhfnz7TF5eaYwnTljHivu3dsUlwceMOszM83YDNi70cqcOGFe+qyvWr0Hc+DAAQYNGsTLL7/M6NGj2b9/P3/4wx9Ys2YN48aNq34H1Ow9GJdUzXswLmluUzdAqjW3qRtggWreg3FJNXtN0CVY9h7MgAEDSEtL4y9/+Qu9evXilVdeISUlpcbFRUREfjlq/TrNww8/zMMPP2xFW0RExI3clP8WmYiIuD4VGBERsYQKjIiIWEIFRkRELKECIyIillCBERERS6jAiIiIJVRgRETEEiowIiJiCRUYERGxhAqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbFErf+XyfVVWloKQEFhY3+yxYqaugHyi+SOx92lpm6ABQqaugENp8CWS9m1vCrNSmsS1YBOnTpFcHBwY36kiIg0sJMnT9KxY8cqYxq9wJSUlHD69GlatWpFs2bNLPucgoICgoODOXnyJH5+fpZ9TmNSTjc/d8sHlJOraKycSktLuXjxIkFBQXh4VD3K0uhdZB4eHtVWvYbk5+fnNgdQGeV083O3fEA5uYrGyMnf379GcRrkFxERS6jAiIiIJdy2wHh7e5OUlIS3t3dTN6XBKKebn7vlA8rJVdyMOTX6IL+IiPwyuO0djIiINC0VGBERsYQKjIiIWEIFRkRELOGWBWb58uV06dIFHx8fIiMj2b9/f1M3qV727NlDbGwsQUFBNGvWjA8//LCpm1QvycnJDBgwgFatWnH77bczatQo/ud//qepm1UvK1eupHfv3tdfcouKiuLTTz9t6mY1qIULF9KsWTOmT5/e1E2ps7lz59KsWTOHqUePHk3drHr5/vvveeyxx2jbti0tW7YkNDSUgwcPNnWzADcsMBs3biQhIYGkpCQyMzMJCwsjJiaGvLy8pm5anRUWFhIWFsby5cubuikNYvfu3UyZMoUvv/yS7du3c+3aNYYNG0Zhoev+C6gdO3Zk4cKFZGRkcPDgQe6//35GjhzJV1991dRNaxAHDhxg9erV9O7du6mbUm/33HMPZ86cuT598cUXTd2kOvvxxx8ZPHgwLVq04NNPP+Xrr79myZIltG7duqmbZpS6mYiIiNIpU6Zcny8uLi4NCgoqTU5ObsJWNRygNC0tramb0aDy8vJKgdLdu3c3dVMaVOvWrUvfeuutpm5GvV28eLE0JCSkdPv27aVDhw4tffbZZ5u6SXWWlJRUGhYW1tTNaDB/+tOfSocMGdLUzXDKre5grl69SkZGBtHR0deXeXh4EB0dTXp6ehO2TKqSn58PQJs2bZq4JQ2juLiY1NRUCgsLiYqKaurm1NuUKVMYMWKEw3nlyr777juCgoK46667GDduHNnZ2U3dpDrbvHkz/fv355FHHuH222+nT58+rF27tqmbdZ1bFZhz585RXFxMQECAw/KAgABycnKaqFVSlZKSEqZPn87gwYPp1atXUzenXo4cOYKvry/e3t5MmjSJtLQ0evbs2dTNqpfU1FQyMzNJTk5u6qY0iMjISNatW8fWrVtZuXIlJ06c4Ne//jUXL15s6qbVyfHjx1m5ciUhISFs27aNyZMn88wzz/Duu+82ddOAJvjXlEXKmzJlCllZWS7dD16me/fuHD58mPz8fP76178SHx/P7t27XbbInDx5kmeffZbt27fj4+PT1M1pEMOHD7/+37179yYyMpLOnTvz/vvvM2HChCZsWd2UlJTQv39/FixYAECfPn3Iyspi1apVxMfHN3Hr3OwOpl27dnh6epKbm+uwPDc3l8DAwCZqlTgzdepUPv74Y3bu3Nmo/wsHq3h5edG1a1f69etHcnIyYWFhvPHGG03drDrLyMggLy+Pvn370rx5c5o3b87u3bt58803ad68OcXFxU3dxHq77bbb6NatG8eOHWvqptTJHXfcccMPmF/96lc3TbefWxUYLy8v+vXrx44dO64vKykpYceOHW7RF+4uSktLmTp1KmlpaXz++efceeedTd0kS5SUlHDlypWmbkad/fa3v+XIkSMcPnz4+tS/f3/GjRvH4cOH8fT0bOom1tulS5f417/+xR133NHUTamTwYMH3/CI/7fffkvnzp2bqEWO3K6LLCEhgfj4ePr3709ERAQpKSkUFhYyfvz4pm5anV26dMnhF9aJEyc4fPgwbdq0oVOnTk3YsrqZMmUKGzZs4D//8z9p1arV9fExf39/WrZs2cStq5vExESGDx9Op06duHjxIhs2bGDXrl1s27atqZtWZ61atbphXOzWW2+lbdu2Ljte9sILLxAbG0vnzp05ffo0SUlJeHp6Mnbs2KZuWp0899xzDBo0iAULFjB69Gj279/PmjVrWLNmTVM3zWjqx9issHTp0tJOnTqVenl5lUZERJR++eWXTd2ketm5c2cpcMMUHx/f1E2rk8pyAUrfeeedpm5anf2f//N/Sjt37lzq5eVV2r59+9Lf/va3pZ999llTN6vBufpjynFxcaV33HFHqZeXV2mHDh1K4+LiSo8dO9bUzaqXjz76qLRXr16l3t7epT169Chds2ZNUzfpOv1z/SIiYgm3GoMREZGbhwqMiIhYQgVGREQsoQIjIiKWUIERERFLqMCIiIglVGBERMQSKjAiImIJFRgREbGECoyIiFhCBUZERCyhAiMiIpb4/4rZeKHF+7lcAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_matrix(model.transmat_)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "b8c1df4f-94fd-483a-a69b-1b7453412e38",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAGdCAYAAAAv9mXmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA490lEQVR4nO3df3DU9Z3H8WcgbIKSBPkViBCCFvkhiggmF2lNi5ypch4OndrpUC7ndYByoco5HQs3zEHt1O11rJQi2tYbQp1rJ8ABlhNLyyDgLyItlDl+FAFxlGp+HNbbAGqA7Pv++G6ymyW72U3yYZP19Zj5zvr9ft/f737e+S778rv7CWSYmSEiItLN+qR6ACIikp4UMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDiRebWfMBgM8sEHH5CTk0NGRsbVfnoREekCM+PcuXMUFBTQp0/8e5SrHjAffPABo0aNutpPKyIi3ejMmTOMHDkybs1VD5icnBwAzpyB3Nyr/ewO5aV6ACIi7jUCowi/l8dz1QOm5WOx3Nw0CxgRkc+QRL7i0Jf8IiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLiRI8NmLVroagIsrOhpAT2749de+kSPP443HijVz95MuzY0bbG74c77oCcHBg2DB54AN5668pz7dsHM2bAtdd6/+LmXXfBJ590U09AEZANlABxWgJgEzA+VH8L8FLU/i3APcBgIAM4FLX/r8C3gXFAf6AQeBgIdHL87VFPPb+ndOsH1BP0jp6wqywQCBhggQBm1v5SXY35fNi6ddjRo9j8+djAgVh9ffv1jz2GFRRg27djb7+NPfMMlp2NHTwYrikvx6qqsCNHsEOHsPvuwwoLsfPnwzVvvIHl5mJ+v1d3/Di2YQP26aexx9q6EH+pBvOBrQM7CjYfbCBYfYz618H6gv0I7BjYcrB+YIcjap4H+x7Yc2CA/SnqHIfB5oBtAzsFtgtsLNhXOhhroot66vk9pVs/6in1PQVC5wsEAh2+39OZkHj66adt9OjRlpWVZcXFxfbmm28mfGwiAVNcjFVWhtebm70A8fvbrx8xAnv66bbb5szB5s6N/RwNDd4Pae/e8LaSEmz58tjHxF06uCjFYJUR681gBWD+GPUPgs2K2lYCtrCd2ndivIDaWzaGXsiXuuEPhXrq+T2lWz/qKfU9JRMwSX9EtmHDBh599FFWrFjBwYMHmTx5MuXl5TQ0NHTpTqrFxYtw4ADMnBne1qePt75vX/vHNDV5H41F6t8fXnst9vMEQveBgwZ5jw0N8Oab3sdnd94J+flQVhb/HIm6CBwAIlqiT2g9Rkvsi6oHKI9Tn6gAkAtkdvE86snTk3tKt35APbXo6T21SDpgnnrqKebPn89DDz3ExIkT+dnPfsY111zDunXrumVAZ89Cc7P3Bh8pPx/q6to/prwcnnoKTp6EYBB27oQtW6C2tv36YBCWLIHp02HSJG/b6dPe48qVMH++9x3O7bfD3Xd75+1ST0AzENUS+UCMlqhLsj7RcXwfWNCFc0SeSz317J7SrZ+Wc6mnnt9Ti6QC5uLFixw4cICZEbcXffr0YebMmeyLcXvR1NREY2Njm6W7rV4NY8fC+PHg88HixfDQQ96dT3sqK+HIEaiuDm8LBr3HhQu9Y6dMgVWrYNw46KbsTKlGYBYwEViZ2qF0G/XU86VbP6CekpFUwJw9e5bm5mbyo24v8vPzqYtxe+H3+8nLy2tdRo0aFfc5hgyBvn2hvr7t9vp6GD68/WOGDoUXXoALF+Ddd+H4cRgwAG644craxYvhxRdh924YOTK8fcQI73HixLb1EybAe+/FHXKHhgB9gaiWqAditMTwJOvjOQd8GcgBtgL9OnGOaOrJ05N7Srd+QD216Ok9tXA+TXnZsmUEAoHW5cyZM3HrfT6YOhV27QpvCwa99dLS+M+VnQ3XXw+XL8PmzTB7dnifmRcuW7fCyy/DmDFtjy0qgoKCK6cunzgBo0d33Gc8PmAqENESwdB6rJZKo+oBdsapj6URb6qiD9iGN6WxO6gnT0/uKd36AfXUoqf31Crh6V9m1tTUZH379rWtW7e22f4P//AP9vd///cJnSPRacpZWdj69dixY9iCBd405bo6b/+8edjSpeH6mhps82ZvivIrr2AzZmBjxmAffRSuWbQIy8vD9uzBamvDy8cfh2tWrfKmKW/ahJ086c0oy87GTp2KPVZLcBZZNVgW2Hq8aYUL8KYh1oX2zwNbGlH/Olgm2JNgfwZbwZXTED8MzQzZHprVUR1ar42Y7VECdgveNMTaiOVyN8x8UU89v6d060c9pb4np9OUi4uLbfHixa3rzc3Ndv3115vf70/o+EQCxgxbs8b7PRWfz5u2XFMT3ldWhlVUhNf37MEmTPBCafBgL4Def7/t+aD9paqqbZ3fj40ciV1zDVZair36avxxWoIBY2BrwArxpgEWg9VE7CsDq4iq3wh2U6j+5tALJXJ/VYyeVoT2746xH7ypi139Q6GeekdP6daPekptT8kETIaZWSJ3Oi02bNhARUUFP//5zykuLuYnP/kJGzdu5Pjx41d8N9OexsZG8vLyCAS835RPGxmpHoCIiHuNQB4QCATI7eBNPOnpzl/72tf43//9X/7t3/6Nuro6brvtNnbs2JFQuIiIyGdH0ncwXaU7GBGR3iuZO5ge+5ddiohI76aAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLiRGbKnjkvZc/shqV6AA5kpHoAItKb6Q5GREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk702IBZCxQB2UAJsL+D+k3A+FD9LcBLUfu3APcAg4EM4FDU/r8C3wbGAf2BQuBhINDJ8bdn7VooKoLsbCgpgf1xmrp0CR5/HG680aufPBl27Ghb4/fDHXdATg4MGwYPPABvvXXlufbtgxkz4NprITcX7roLPvmkm3oiDa8T6dVTuvUD6gl6R0/YVRYIBAywAJjFWKrBfGDrwI6CzQcbCFYfo/51sL5gPwI7BrYcrB/Y4Yia58G+B/YcGGB/ijrHYbA5YNvAToHtAhsL9pU442yzWPyluhrz+bB167CjR7H587GBA7H6+vbrH3sMKyjAtm/H3n4be+YZLDsbO3gwXFNejlVVYUeOYIcOYffdhxUWYufPh2veeAPLzcX8fq/u+HFswwbs0087HnNHPffK6/QZ6ynd+lFPqe8pEDpfIBDo8P2eZANi79699nd/93c2YsQIA2zr1q1JHZ9IwBSDVUasN4MVgPlj1D8INitqWwnYwnZq34nxw25v2Ri66JcSeZF08GZdXIxVVobXm5u9APH7268fMQJ7+um22+bMwebOjf0cDQ1eb3v3hreVlGDLl3c8vnaXDnruldfpM9ZTuvWjnlLfUzIBk/RHZBcuXGDy5MmsXbu2KzdOMV0EDgAzI7b1Ca3vi3HMvqh6gPI49YkKALlAZhfPc/EiHDgAMyMG2aePt74vxiCbmryPxiL17w+vvRZnvKF720GDvMeGBnjzTe/jszvvhPx8KCuLf45EpeV1Ir16Srd+QD216Ok9tUj6PPfeey/33ntvNz39lc4CzUB+1PZ84HiMY+pi1Nd1cRzfBxZ04Ryt5zoLzc3eG3yk/Hw4HqOp8nJ46inv+5Ibb4Rdu2DLFu887QkGYckSmD4dJk3ytp0+7T2uXAlPPgm33QbPPw933w1HjsDYsV3oiTS8TqRXT+nWT8u51FPP76mF8y/5m5qaaGxsbLP0dI3ALGAisDJFY1i92guA8ePB54PFi+Ghh7w7n/ZUVnqhUV0d3hYMeo8LF3rHTpkCq1bBuHGwbp37HlzrCdepu6VbT+nWD6inZDgPGL/fT15eXusyatSouPVDgL5AfdT2emB4jGOGJ1kfzzngy0AOsBXo14lzRBsyBPr2hfqoQdbXw/AYgxw6FF54AS5cgHff9e50BgyAG264snbxYnjxRdi9G0aODG8fMcJ7nDixbf2ECfDee51uB0jT60R69ZRu/YB6atHTe2rhPGCWLVtGIBBoXc6cORO33gdMBXZFbAuG1ktjHFMaVQ+wM059LI140/p8wDa86X/dweeDqVO9j7laBIPeemkHg8zOhuuvh8uXYfNmmD07vM/MC5etW+Hll2HMmLbHFhVBQcGVU5dPnIDRo7vUUnpeJ9Krp3TrB9RTi57eU6ukpoBFATezyKrBssDW403BW4A3Za8utH8e2NKI+tfBMsGeBPsz2AqunLL3YWgWxfbQDIjq0HptxMyIErBb8Kbs1UYslxOZCZLANOWsLGz9euzYMWzBAm+acl2dt3/ePGzp0nB9TQ22ebM3RfmVV7AZM7AxY7CPPgrXLFqE5eVhe/ZgtbXh5eOPwzWrVnnTlDdtwk6e9GaUZWdjp051POaOeu6V1+kz1lO69aOeUt+T02nKbQ7GTcAY2BqwQrwpc8VgNRH7ysAqouo3gt0Uqr859EON3F8V+qFELytC+3fH2A/eNL8OXyQJTPtds8b7PRWfz5u2XFMT3ldWhlVUhNf37MEmTPBCafBgL4Def7/t+WKNt6qqbZ3fj40ciV1zDVZair36amLjTeQPRq+7Tp/BntKtH/WU2p6SCZiMUFAk7Pz585w6dQqAKVOm8NRTT/GlL32JQYMGUVhY2OHxjY2N5OXltU6HSxtJ/RR7iYxUD0BEeppGIA8IBALk5sZ/F086YPbs2cOXvvSlK7ZXVFSwfv36jgengOk9FDAiEiWZgEn692C++MUvkmQmiYjIZ1CP/csuRUSkd1PAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxIjPVA0gbGakegAOW6gF0s3S8RiI9mO5gRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuJEjw2YtUARkA2UAPs7qN8EjA/V3wK8FLV/C3APMBjIAA5F7f8r8G1gHNAfKAQeBgKdHH970rKntVBUBNnZUFIC++M0dekSPP443HijVz95MuzY0bbG74c77oCcHBg2DB54AN5668pz7dsHM2bAtddCbi7cdRd88kk39UR6Xad06wfUE/SOnrCrLBAIGGABMIuxVIP5wNaBHQWbDzYQrD5G/etgfcF+BHYMbDlYP7DDETXPg30P7DkwwP4UdY7DYHPAtoGdAtsFNhbsK3HGmczSK3uy+Et1NebzYevWYUePYvPnYwMHYvX17dc/9hhWUIBt3469/Tb2zDNYdjZ28GC4prwcq6rCjhzBDh3C7rsPKyzEzp8P17zxBpabi/n9Xt3x49iGDdinn3Yw5nS9Tp+hftRT6nsKhM4XCAQ6fL8nmXB44oknbNq0aTZgwAAbOnSozZ49244fP57MKRIKmGKwyoj1ZrACMH+M+gfBZkVtKwFb2E7tOzF+2O0tG0MX/VI3vIB6ZU8dBExxMVZZGV5vbvYCxO9vv37ECOzpp9tumzMHmzs39nM0NHi97d0b3lZSgi1f3vH4rljS9Tp9hvpRT6nvKZmASeojsr1791JZWUlNTQ07d+7k0qVL3HPPPVy4cKHLd1ItLgIHgJkR2/qE1vfFOGZfVD1AeZz6RAWAXCCzi+dJy54uwoEDMDNikH36eOv7Ygyyqcn7aCxS//7w2mtxxhu6Xx80yHtsaIA33/Q+PrvzTsjPh7Ky+OdIVLpdp3TrB9RTi57eU4ukzrMj6gPz9evXM2zYMA4cOMBdd93VLQM6CzQD+VHb84HjMY6pi1Ff18VxfB9Y0IVzRJ4r7Xo6C83N3ht8pPx8OB6jqfJyeOop7/uSG2+EXbtgyxbvPO0JBmHJEpg+HSZN8radPu09rlwJTz4Jt90Gzz8Pd98NR47A2LFd6In0uk7p1k/LudRTz++pRZe+5A+E/vdyUMv/XrajqamJxsbGNktP1wjMAiYCK1M7lG7TE3pavdoLgPHjweeDxYvhoYe8O5/2VFZ6oVFdHd4WDHqPCxd6x06ZAqtWwbhxsG6d+x5c6wnXqTulWz+gnpLR6YAJBoMsWbKE6dOnM6nlfy/b4ff7ycvLa11GjRoV97xDgL5AfdT2emB4jGOGJ1kfzzngy0AOsBXo14lzREvLnoZA375QHzXI+noYHmOQQ4fCCy/AhQvw7rvenc6AAXDDDVfWLl4ML74Iu3fDyJHh7SNGeI8TJ7atnzAB3nuv0+0A6Xed0q0fUE8tenpPLTodMJWVlRw5coTqyP+9bMeyZcsIBAKty5kzZ+LW+4CpwK6IbcHQemmMY0qj6gF2xqmPpRFvWp8P2IY3/a87pGVPPpg61fuYq0Uw6K2XdjDI7Gy4/nq4fBk2b4bZs8P7zLxw2boVXn4Zxoxpe2xRERQUXDl1+cQJGD26Sy2l3XVKt35APbXo6T21SmoKWEhlZaWNHDnSTp8+nfSxiU5TzgJbjzcFbwHelL260P55YEsj6l8HywR7EuzPYCu4csreh6FZFNtDMyCqQ+u1ETMjSsBuwZuyVxuxXO6GWSK9sqcEpilnZWHr12PHjmELFnjTlOvqvP3z5mFLl4bra2qwzZu9KcqvvILNmIGNGYN99FG4ZtEiLC8P27MHq60NLx9/HK5ZtcqbprxpE3bypDejLDsbO3WqgzGn63X6DPWjnlLfk7NpysFg0CorK62goMBOnDiRzKGtEgkYA1sDVog3Za4YrCZiXxlYRVT9RrCbQvU3h36okfurQj+U6GVFaP/uGPvBm+bX1RdQr+wpgam/a9Z4v6fi83nTlmtqwvvKyrCKivD6nj3YhAleKA0e7AXQ+++3PV+s8VZVta3z+7GRI7FrrsFKS7FXX01gvOl6nT5j/ain1PaUTMBkmJklcqcD8M///M/8+te/5je/+Q3jxo1r3Z6Xl0f//v0TOkdjYyN5eXmt0+GkB0v4ldFLZKR6ACK9XyOQhzfJKzc3/rt4UgGTkdH+n9Cqqir+8R//MbHBKWB6DwWMiERJJmCS+j2YJLJIREQ+43rsX3YpIiK9mwJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScSOqfTJbPmHT7N+zT8V/8TrdrJGlFdzAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTvTYgFkLFAHZQAmwv4P6TcD4UP0twEtR+7cA9wCDgQzgUNT+vwLfBsYB/YFC4GEg0Mnxt0c99fye1q6FoiLIzoaSEtgfp6FLl+Dxx+HGG736yZNhx462NX4/3HEH5OTAsGHwwAPw1ltXnmvfPpgxA669FnJz4a674JNPuqkn0usagXqC3tETloRnnnnGbrnlFsvJybGcnBz7m7/5G3vppZeSOYUFAgEDLABmMZZqMB/YOrCjYPPBBoLVx6h/Hawv2I/AjoEtB+sHdjii5nmw74E9BwbYn6LOcRhsDtg2sFNgu8DGgn0lzjiTWdRTD+jJ4i/V1ZjPh61bhx09is2fjw0ciNXXt1//2GNYQQG2fTv29tvYM89g2dnYwYPhmvJyrKoKO3IEO3QIu+8+rLAQO38+XPPGG1huLub3e3XHj2MbNmCfftrxmNPuGqXj6y7NegqEzhcIBDp8vyeZcNi2bZtt377dTpw4YW+99Zb967/+q/Xr18+OHDmS8DkSCZhisMqI9WawAjB/jPoHwWZFbSsBW9hO7TsxftjtLRtDF/1SN7yA1FMP6KmDN+viYqyyMrze3OwFiN/ffv2IEdjTT7fdNmcONndu7OdoaPD62rs3vK2kBFu+vOPxtbuk2zVKx9ddmvWUTMAk9RHZ/fffz3333cfYsWO56aab+MEPfsCAAQOoqanp8p1Ui4vAAWBmxLY+ofV9MY7ZF1UPUB6nPlEBIBfI7OJ51JOnJ/d08SIcOAAzIwbYp4+3vi/GAJuavI/GIvXvD6+9Fmesoc8fBg3yHhsa4M03vY/P7rwT8vOhrCz+ORKVbtcI1FOLnt5Ti05/B9Pc3Ex1dTUXLlygtLQ0Zl1TUxONjY1tlnjOAs1AftT2fKAuxjF1SdYn4izwfWBBF84ReS711LN7OnsWmpu9N/hI+flQF2OA5eXw1FNw8iQEg7BzJ2zZArW17dcHg7BkCUyfDpMmedtOn/YeV66E+fO973Buvx3uvts7b5d6Ir2uUcu51FPP76lF0gFz+PBhBgwYQFZWFt/61rfYunUrEydOjFnv9/vJy8trXUaNGtWlAV8NjcAsYCKwMrVD6TbqqfutXg1jx8L48eDzweLF8NBD3p1Peyor4cgRqK4ObwsGvceFC71jp0yBVatg3DhYt859D66l+hq5oJ4Sl3TAjBs3jkOHDvHmm2+yaNEiKioqOHbsWMz6ZcuWEQgEWpczZ87EPf8QoC9QH7W9Hhge45jhSdbHcw74MpADbAX6deIc0dSTpyf3NGQI9O0L9VEDrK+H4TEGOHQovPACXLgA774Lx4/DgAFwww1X1i5eDC++CLt3w8iR4e0jRniP0f+PNmECvPdep9sB0u8agXpq0dN7apF0wPh8Pj73uc8xdepU/H4/kydPZvXq1THrs7KyyM3NbbPEPT8wFdgVsS0YWo/1QVxpVD3Azjj1sTTiTevzAdvwpv91B/Xk6ck9+XwwdSrsihhgMOitx/kEGPC+h7n+erh8GTZvhtmzw/vMvHDZuhVefhnGjGl7bFERFBRcOXX5xAkYPbpLLaXdNQL11KKn99Qq4elfMXzpS1+yioqKhOsTnaacBbYebwreArwpe3Wh/fPAlkbUvw6WCfYk2J/BVnDllL0PQ7MotodmQFSH1msjZkaUgN2CN2WvNmK53A2zRNRTD+gpgWnKWVnY+vXYsWPYggXeNOW6Om//vHnY0qXh+poabPNmb4ryK69gM2ZgY8ZgH30Urlm0CMvLw/bswWprw8vHH4drVq3ypilv2oSdPOnNKMvOxk6d6njMaXeN0vF1l2Y9OZumvHTpUtu7d6+988479j//8z+2dOlSy8jIsN///vcJnyORgDGwNWCFeFPmisFqIvaVgVVE1W8EuylUf3Pohxq5vyr0Q4leVoT2746xH7xpfl19AamnHtBTAtN+16zxfk/F5/OmLdfUhPeVlWEVFeH1PXuwCRO8UBo82Aug999ve75YY62qalvn92MjR2LXXIOVlmKvvprYeNPuGqXj6y7NekomYDLMzBK50wH45je/ya5du6itrSUvL49bb72V7373u/zt3/5toqegsbGRvLy81ulwIldNwq/0XiQj1QOQz5pGIA8IBAIdfuWRVMB0BwWMpIwCRqTLkgmYHvt3kYmISO+mgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4kRmqgcgctWk479fb6kegAPpeJ0+o3QHIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuJEjw2YtUARkA2UAPs7qN8EjA/V3wK8FLV/C3APMBjIAA5F7f8r8G1gHNAfKAQeBgKdHH971JN6SkVPa9dCURFkZ0NJCeyP09ClS/D443DjjV795MmwY0fbGr8f7rgDcnJg2DB44AF4660rz7VvH8yYAddeC7m5cNdd8Mkn3dQT6XWNID17wrrA7/cbYI888kjCxwQCAQMsAGYxlmowH9g6sKNg88EGgtXHqH8drC/Yj8COgS0H6wd2OKLmebDvgT0HBtifos5xGGwO2DawU2C7wMaCfSXOOJNZ1JN6ctKTxV+qqzGfD1u3Djt6FJs/Hxs4EKuvb7/+scewggJs+3bs7bexZ57BsrOxgwfDNeXlWFUVduQIdugQdt99WGEhdv58uOaNN7DcXMzv9+qOH8c2bMA+/bTjMafdNUqz110gdL5AINDh+z0JJ0OU/fv3W1FRkd16663dHjDFYJUR681gBWD+GPUPgs2K2lYCtrCd2ndi/LDbWzaGLvqlbngBqSf15KSnDt6si4uxysrwenOzFyB+f/v1I0ZgTz/ddtucOdjcubGfo6HB62vv3vC2khJs+fKOx9fukm7XKM1ed8kETKc+Ijt//jxz587lueee47rrruvszVO7LgIHgJkR2/qE1vfFOGZfVD1AeZz6RAWAXCCzi+dRTx71lLju6OniRThwAGZGDLBPH299X4wBNjV5H41F6t8fXnstzlhDn6kMGuQ9NjTAm296H5/deSfk50NZWfxzJCrdrhGkZ08tOhUwlZWVzJo1i5kzo1u8UlNTE42NjW2WeM4CzUB+1PZ8oC7GMXVJ1ifiLPB9YEEXzhF5LvWknpIZR3f0dPYsNDd7b/CR8vOhLsYAy8vhqafg5EkIBmHnTtiyBWpr268PBmHJEpg+HSZN8radPu09rlwJ8+d73+Hcfjvcfbd33i71RHpdo5ZzpVtPLZIOmOrqag4ePIjf70+o3u/3k5eX17qMGjUq6UFebY3ALGAisDK1Q+k26ql3SHVPq1fD2LEwfjz4fLB4MTz0kHfn057KSjhyBKqrw9uCQe9x4ULv2ClTYNUqGDcO1q1z34Nrqb5GLrjqKamAOXPmDI888gi/+tWvyI6+j45h2bJlBAKB1uXMmTNx64cAfYH6qO31wPAYxwxPsj6ec8CXgRxgK9CvE+eIpp486im+7u5pyBDo2xfqowZYXw/DYwxw6FB44QW4cAHefReOH4cBA+CGG66sXbwYXnwRdu+GkSPD20eM8B4nTmxbP2ECvPdep9sB0u8aQXr21CKpgDlw4AANDQ3cfvvtZGZmkpmZyd69e/npT39KZmYmzc3NVxyTlZVFbm5umyUeHzAV2BWxLRhaL41xTGlUPcDOOPWxNOJN6/MB2/Cm/3UH9eRRT7G56Mnng6lTYVfEAINBb720gwFmZ8P118Ply7B5M8yeHd5n5oXL1q3w8sswZkzbY4uKoKDgyqnLJ07A6NFdaintrhGkZ0+tEp7+ZWaNjY12+PDhNsu0adPsG9/4hh0+fDihcyQ6TTkLbD3eFLwFeFP26kL754Etjah/HSwT7EmwP4Ot4Mopex+GZlFsD82AqA6t10bMjCgBuwVvyl5txHK5G2aJqCf15KSnBKYpZ2Vh69djx45hCxZ405Tr6rz98+ZhS5eG62tqsM2bvSnKr7yCzZiBjRmDffRRuGbRIiwvD9uzB6utDS8ffxyuWbXKm6a8aRN28qQ3oyw7Gzt1quMxp901SrPX3VWZptyirKys26cpG9gasEK8KXPFYDUR+8rAKqLqN4LdFKq/OfRDjdxfFfqhRC8rQvt3x9gP3jS/rr6A1JN6ctJTAtN+16zxfk/F5/OmLdfUhPeVlWEVFeH1PXuwCRO8UBo82Aug999ve75YY62qalvn92MjR2LXXIOVlmKvvprYeNPuGqXZ6y6ZgMkwM0vkTieWL37xi9x222385Cc/Sai+sbGRvLy81ulwItIFXfrT20NlpHoAEk8jkAcEAoEOv/Lo8nTnPXv2dPUUIiKShnrs30UmIiK9mwJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTmakegIh0QUaqB+CApXoADqTjdUqA7mBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4kSPDZi1QBGQDZQA+zuo3wSMD9XfArwUtX8LcA8wGMgADkXt/yvwbWAc0B8oBB4GAp0cf3vUk3pKRU/p1g/A2rVQVATZ2VBSAvvjNHXpEjz+ONx4o1c/eTLs2NG2xu+HO+6AnBwYNgweeADeeuvKc+3bBzNmwLXXQm4u3HUXfPJJN/VE+l0n7CoLBAIGWADMYizVYD6wdWBHweaDDQSrj1H/OlhfsB+BHQNbDtYP7HBEzfNg3wN7DgywP0Wd4zDYHLBtYKfAdoGNBftKnHEms6gn9ZSKnnplPxZ/qa7GfD5s3Trs6FFs/nxs4ECsvr79+scewwoKsO3bsbffxp55BsvOxg4eDNeUl2NVVdiRI9ihQ9h992GFhdj58+GaN97AcnMxv9+rO34c27AB+/TTjsecTtcpEDpfIBDo8P2eZMJhxYoVRujkLcu4ceOSOUVCAVMMVhmx3gxWAOaPUf8g2KyobSVgC9upfSfGD7u9ZWPool9K9A+GelJPPaynXtlPB2/WxcVYZWV4vbnZCxC/v/36ESOwp59uu23OHGzu3NjP0dDg9bZ3b3hbSQm2fHnH42t3SaPrlEzAJP0R2c0330xtbW3r8tprr3X67qk9F4EDwMyIbX1C6/tiHLMvqh6gPE59ogJALpDZxfOoJ496Slx39JRu/QBcvAgHDsDMiEH26eOt74sxyKYm76OxSP37Q7y3rkDoc6JBg7zHhgZ4803v47M774T8fCgri3+ORKXjdWqRdMBkZmYyfPjw1mXIkCHdNBTPWaAZyI/ang/UxTimLsn6RMfxfWBBF84ReS71pJ6SGUd39JRu/QCcPQvNzd4bfKT8fKiLMcjycnjqKTh5EoJB2LkTtmyB2tr264NBWLIEpk+HSZO8badPe48rV8L8+d53OLffDnff7Z23Sz2RftepRdIBc/LkSQoKCrjhhhuYO3cu7733Xtz6pqYmGhsb2yw9XSMwC5gIrEztULqNeuod0q2nntDP6tUwdiyMHw8+HyxeDA895N35tKeyEo4cgerq8LZg0HtcuNA7dsoUWLUKxo2Ddevc9+Caq+uUVMCUlJSwfv16duzYwbPPPss777zDF77wBc6dOxfzGL/fT15eXusyatSouM8xBOgL1EdtrweGxzhmeJL18ZwDvgzkAFuBfp04RzT15FFP8XV3T+nWD8CQIdC3L9RHDbK+HobHGOTQofDCC3DhArz7Lhw/DgMGwA03XFm7eDG8+CLs3g0jR4a3jxjhPU6c2LZ+wgTo4P+xO5SO16lFUgFz77338tWvfpVbb72V8vJyXnrpJf7v//6PjRs3xjxm2bJlBAKB1uXMmTNxn8MHTAV2RWwLhtZLYxxTGlUPsDNOfSyNeNP6fMA2vOl/3UE9edRTbC56Srd+wLsDmToVdkUMMhj01ks7GGR2Nlx/PVy+DJs3w+zZ4X1mXrhs3QovvwxjxrQ9tqgICgqunLp84gSMHt2lltLyOrVKagpYO6ZNm2ZLly5NuD7RacpZYOvxpuAtwJuyVxfaPw9saUT962CZYE+C/RlsBVdO2fswNItie2gGRHVovTZiZkQJ2C14U/ZqI5bLicx86WBRT+opFT31yn4SmKaclYWtX48dO4YtWOBNU66r8/bPm4ctXRqur6nBNm/2pii/8go2YwY2Zgz20UfhmkWLsLw8bM8erLY2vHz8cbhm1SpvmvKmTdjJk96Msuxs7NSpjsecTtfJ2TTlaOfOnbPrrrvOVq9enfAxiQSMga0BK8SbMlcMVhOxrwysIqp+I9hNofqbQz/UyP1VtJ1e3bKsCO3fHWM/eNP8uvrGpZ7UU6p66nX9JDDtd80a7/dUfD5v2nJNTXhfWRlWURFe37MHmzDBC6XBg70Aev/9tueLNd6qqrZ1fj82ciR2zTVYaSn26quJjTedrlMyAZNhZpbInQ7Ad77zHe6//35Gjx7NBx98wIoVKzh06BDHjh1j6NChCZ2jsbGRvLy81ulwIiJtJPyO1ItkpHoA3acRyAMCgQC5ufHfxZOa7vyXv/yFr3/963z44YcMHTqUz3/+89TU1CQcLiIi8tmR1B1Md9AdjIjEpTuYHi2ZO5ge+5ddiohI76aAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLiRGaqByAi0kYa/fv1rSzVA+hGjUBeYqW6gxEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEwoYERFxQgEjIiJOKGBERMQJBYyIiDihgBEREScUMCIi4oQCRkREnFDAiIiIEz02YNYCRUA2UALs76B+EzA+VH8L8FLU/i3APcBgIAM4FLX/r8C3gXFAf6AQeBgIdHL87VFP6ikVPaVbP5CmPa2FoiLIzoaSEtgfp6lLl+Dxx+HGG736yZNhx462NX4/3HEH5OTAsGHwwAPw1ltXnmvfPpgxA669FnJz4a674JNPuqkpu8oCgYABFgCzGEs1mA9sHdhRsPlgA8HqY9S/DtYX7Edgx8CWg/UDOxxR8zzY98CeAwPsT1HnOAw2B2wb2CmwXWBjwb4SZ5zJLOpJPaWip3Trp9f2ZPGX6mrM58PWrcOOHsXmz8cGDsTq69uvf+wxrKAA274de/tt7JlnsOxs7ODBcE15OVZVhR05gh06hN13H1ZYiJ0/H6554w0sNxfz+72648exDRuwTz+NPdZAwPsZBQKBDt/vSTYg/vKXv9jcuXNt0KBBlp2dbZMmTbI//OEPCR+fSMAUg1VGrDeDFYD5Y9Q/CDYralsJ2MJ2at+J8QJqb9kYeiFf6oY/FOpJPaWip3Trp9f21EHAFBdjlZXh9eZmL0D8/vbrR4zAnn667bY5c7C5c2M/R0OD19veveFtJSXY8uUdjy9ySSZgkvqI7KOPPmL69On069eP3/72txw7dowf//jHXHfddV2+k2pxETgAzIzY1ie0vi/GMfui6gHK49QnKgDkApldPI968qinxHVHT+nWD6RpTxfhwAGYGTHIPn289X0xBtnU5H00Fql/f3jttTjjDX2eN2iQ99jQAG++6X18duedkJ8PZWXxz5GspH42//7v/86oUaOoqqpq3TZmzJjuGw1wFmgG8qO25wPHYxxTF6O+rovj+D6woAvniDyXelJPyYyjO3pKt35azpV2PZ2F5mbvDT5Sfj4cj9FUeTk89ZT3fcmNN8KuXbBli3ee9gSDsGQJTJ8OkyZ5206f9h5XroQnn4TbboPnn4e774YjR2Ds2K73ltQdzLZt25g2bRpf/epXGTZsGFOmTOG5556Le0xTUxONjY1tlp6uEZgFTARWpnYo3UY99Q7p1lO69QM9o6fVq70AGD8efD5YvBgeesi782lPZaUXGtXV4W3BoPe4cKF37JQpsGoVjBsH69Z1zziTCpjTp0/z7LPPMnbsWH73u9+xaNEiHn74YX75y1/GPMbv95OXl9e6jBo1Ku5zDAH6AvVR2+uB4TGOGZ5kfTzngC8DOcBWoF8nzhFNPXnUU3zd3VO69QNp2tMQ6NsX6qMGWV8Pw2MMcuhQeOEFuHAB3n3Xu9MZMABuuOHK2sWL4cUXYfduGDkyvH3ECO9x4sS29RMmwHvvdbqdNpIKmGAwyO23384TTzzBlClTWLBgAfPnz+dnP/tZzGOWLVtGIBBoXc6cORP3OXzAVGBX5POG1ktjHFMaVQ+wM059LI14UxV9wDa8KY3dQT151FNsLnpKt34gTXvywdSp3sdcLYJBb720g0FmZ8P118Ply7B5M8yeHd5n5oXL1q3w8ssQ/W1GUREUFFw5dfnECRg9ukstRQ4icYWFhfbNb36zzbZnnnnGCgoKEj5HotOUs8DW400rXIA3DbEutH8e2NKI+tfBMsGeBPsz2AqunIb4YWhmyPbQLJHq0HptaH8gNLPkFrxpiLURy+VumPmintRTKnpKt356bU8JTFPOysLWr8eOHcMWLPCmKdfVefvnzcOWLg3X19Rgmzd7U5RfeQWbMQMbMwb76KNwzaJFWF4etmcPVlsbXj7+OFyzapU3TXnTJuzkSW9GWXY2dupU7LE6m6b89a9/3T7/+c+32bZkyRIrLS1N+ByJBIyBrQErxJsGWAxWE7GvDKwiqn4j2E2h+ptDL5TI/VWhF070siK0f3eM/eBNXezqHwr1pJ5S1VO69dMre0pg+u+aNd7vqfh83rTlmprwvrIyrKIivL5nDzZhghdKgwd7AfT++23PF2u8VVVt6/x+bORI7JprsNJS7NVX44/TWcDs37/fMjMz7Qc/+IGdPHnSfvWrX9k111xj//mf/5nwORINGC1atGhJmyWJ3zPp6UsyAZNhZtbxB2lhL774IsuWLePkyZOMGTOGRx99lPnz5yd8fGNjI3l5ea1zyEVE0l5S77I9W2Mj5OVBIBAgNzf+u3jSAdNVChgR+cz5jAZMj/3LLkVEpHdTwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBMKGBERcUIBIyIiTihgRETECQWMiIg4oYAREREnFDAiIuKEAkZERJxQwIiIiBOZV/sJW/6F5sar/cQiIqmSRm94jaFeWt7L47nqAXPu3DkARl3tJxYRSZW8VA+g+507d468vPiNZVgiMdSNgsEgH3zwATk5OWRkZDh7nsbGRkaNGsWZM2fIzc119jxXk3rq+dKtH1BPvcXV6snMOHfuHAUFBfTpE/9blqt+B9OnTx9Gjhx51Z4vNzc3bV5ALdRTz5du/YB66i2uRk8d3bm00Jf8IiLihAJGREScSNuAycrKYsWKFWRlZaV6KN1GPfV86dYPqKfeoif2dNW/5BcRkc+GtL2DERGR1FLAiIiIEwoYERFxQgEjIiJOpGXArF27lqKiIrKzsykpKWH//v2pHlKXvPLKK9x///0UFBSQkZHBCy+8kOohdYnf7+eOO+4gJyeHYcOG8cADD/DWW2+lelhd8uyzz3Lrrbe2/pJbaWkpv/3tb1M9rG71wx/+kIyMDJYsWZLqoXTaypUrycjIaLOMHz8+1cPqkvfff59vfOMbDB48mP79+3PLLbfwxz/+MdXDAtIwYDZs2MCjjz7KihUrOHjwIJMnT6a8vJyGhoZUD63TLly4wOTJk1m7dm2qh9It9u7dS2VlJTU1NezcuZNLly5xzz33cOHChVQPrdNGjhzJD3/4Qw4cOMAf//hHZsyYwezZszl69Giqh9Yt/vCHP/Dzn/+cW2+9NdVD6bKbb76Z2tra1uW1115L9ZA67aOPPmL69On069eP3/72txw7dowf//jHXHfddakemsfSTHFxsVVWVrauNzc3W0FBgfn9/hSOqvsAtnXr1lQPo1s1NDQYYHv37k31ULrVddddZ//xH/+R6mF02blz52zs2LG2c+dOKysrs0ceeSTVQ+q0FStW2OTJk1M9jG7z3e9+1z7/+c+nehgxpdUdzMWLFzlw4AAzZ85s3danTx9mzpzJvn37UjgyiScQCAAwaNCgFI+kezQ3N1NdXc2FCxcoLS1N9XC6rLKyklmzZrX5c9WbnTx5koKCAm644Qbmzp3Le++9l+ohddq2bduYNm0aX/3qVxk2bBhTpkzhueeeS/WwWqVVwJw9e5bm5mby8/PbbM/Pz6euri5Fo5J4gsEgS5YsYfr06UyaNCnVw+mSw4cPM2DAALKysvjWt77F1q1bmThxYqqH1SXV1dUcPHgQv9+f6qF0i5KSEtavX8+OHTt49tlneeedd/jCF77Q+s+I9DanT5/m2WefZezYsfzud79j0aJFPPzww/zyl79M9dCAFPxtyiKRKisrOXLkSK/+HLzFuHHjOHToEIFAgP/6r/+ioqKCvXv39tqQOXPmDI888gg7d+4kOzs71cPpFvfee2/rf996662UlJQwevRoNm7cyDe/+c0UjqxzgsEg06ZN44knngBgypQpHDlyhJ/97GdUVFSkeHRpdgczZMgQ+vbtS319fZvt9fX1DB8+PEWjklgWL17Miy++yO7du6/qP+Hgis/n43Of+xxTp07F7/czefJkVq9enephddqBAwdoaGjg9ttvJzMzk8zMTPbu3ctPf/pTMjMzaW5uTvUQu2zgwIHcdNNNnDp1KtVD6ZQRI0Zc8T8wEyZM6DEf+6VVwPh8PqZOncquXbtatwWDQXbt2pUWn4WnCzNj8eLFbN26lZdffpkxY8akekhOBINBmpqaUj2MTrv77rs5fPgwhw4dal2mTZvG3LlzOXToEH379k31ELvs/PnzvP3224wYMSLVQ+mU6dOnXzHF/8SJE4wePTpFI2or7T4ie/TRR6moqGDatGkUFxfzk5/8hAsXLvDQQw+lemiddv78+Tb/h/XOO+9w6NAhBg0aRGFhYQpH1jmVlZX8+te/5je/+Q05OTmt34/l5eXRv3//FI+uc5YtW8a9995LYWEh586d49e//jV79uzhd7/7XaqH1mk5OTlXfC927bXXMnjw4F77fdl3vvMd7r//fkaPHs0HH3zAihUr6Nu3L1//+tdTPbRO+Zd/+RfuvPNOnnjiCR588EH279/PL37xC37xi1+kemieVE9jc2HNmjVWWFhoPp/PiouLraamJtVD6pLdu3cbcMVSUVGR6qF1Snu9AFZVVZXqoXXaP/3TP9no0aPN5/PZ0KFD7e6777bf//73qR5Wt+vt05S/9rWv2YgRI8zn89n1119vX/va1+zUqVOpHlaX/Pd//7dNmjTJsrKybPz48faLX/wi1UNqpb+uX0REnEir72BERKTnUMCIiIgTChgREXFCASMiIk4oYERExAkFjIiIOKGAERERJxQwIiLihAJGREScUMCIiIgTChgREXFCASMiIk78P0tmIWkW5ZvWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_matrix(model.emissionprob_)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "a266979b-7855-4baf-816d-17ea792419b7",
   "metadata": {
    "pycharm": {
     "name": "#%%\n"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGYAAAGdCAYAAADg/QYMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAW0UlEQVR4nO2dfXBU9bnHP3ndaONuSAkJgRDkLQgIUjJJA5bYIa1anMttZ9RSCql1QJhYYWC0aJVcbWvUi61WKbfqQDqWe6OOUm0poEOBEQoBgwyB8BbMSFASFHATQgVJnvvH2YQkZsOeZBMe5fnM/CZ7znl+v3POfs5bNsuXCBERDHVEXu4NMDrGxCjFxCjFxCjFxCjFxCjFxCjFxCglurdX2NTUxMcff8w111xDREREb68+7IgI9fX1pKamEhkZvuO818V8/PHHpKWl9fZqe5zq6moGDhwYtvF6Xcw111wDQHU1eL29vfbwU1cHaWkX9ytc9LqY5suX1/v1ENNMuC/LdvNXiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiolRiloxy5bB4MEQFwfZ2bBjR/DaL76Axx6DoUOd+nHjYN26tjX19bBgAaSnw1VXwcSJsHNn8DHnzoWICHjmmTDsTFeQXsbv9wsgfj8i0nErKUFiY5EVK5B9+5DZs5GEBKS2tuP6Bx5AUlORNWuQI0eQP/4RiYtDdu26WHPHHcioUcjmzcjhw0hhIeL1IseOfXm8N95Axo1zxvz974Nvp4izH87++MP6PnVJzPPPPy/p6eni8XgkKytLSktLQ+4bipisLKSg4OJ0Y6PzJhUVdVzfvz/y/PNt5/3oR8iMGc7rs2eRqCjk739vW/OtbyG/+lXbeceOIQMGIHv3Iunpl0+M60vZK6+8wsKFCyksLGTXrl2MGzeOm2++mRMnToTlDD5/HsrKIC/v4rzISGd627aO+5w751zCWnPVVbBli/P6wgVobOy8BqCpCWbOhPvvh9Gju78v3cKtyaysLCkoKGiZbmxslNTUVCkqKgqp/6XOmI8+co7Af/2r7fz773fOpI76TJ/uXKYOHXLOrrffRq66yrkcNtfk5CC5uc74Fy4gL7+MREYiI0ZcrHn8ceR730Oampzpr8wZc/78ecrKyshrdThHRkaSl5fHtiCH87lz56irq2vTws2zz8Lw4TByJMTGwr33wl13OWdaMy+/DCIwYAB4PPCHP8D06RdrysqccYqLnZv+5caVmE8//ZTGxkaSk5PbzE9OTqampqbDPkVFRfh8vpZ2qW9h9u0LUVFQW9t2fm0tpKR03CcpCf76V2hogA8/hAMHID4ehgy5WDN0KGzeDGfOOF823LHDeZprrnn3XThxAgYNguhop334ISxa5Dwd9jY9/rj84IMP4vf7W1p1dXWn9bGxMGECbNhwcV5TkzOdk9P5uuLinDPiwgV4/XWYNu3LNd/4BvTvD6dPw/r1F2tmzoQ9e2D37ostNdW536xf72aPw4Orb2L27duXqKgoatsdzrW1taQEOZw9Hg8ej8fVRi1cCPn5kJkJWVnO7xINDc7lCWDWLEdAUZEzXVoKH30EN9zg/Pyv/3JkPvDAxTHXr3cuZRkZUFnpvOEjR14c85vfdFprYmKcszQjw9XmhwVXZ0xsbCwTJkxgQ6vDuampiQ0bNpBzqcPZBXfeCUuXwpIlzpu9e7fzC2PzFfToUTh+/GL955/Dww/DqFHwwx860rZsgYSEizV+PxQUODJmzYIbb3RkxcSEbbPDi9unhZKSEvF4PFJcXCwVFRUyZ84cSUhIkJqampD6h/J7zFep9dRTmesvld9555188sknLFmyhJqaGm644QbWrVv3pQcCo3tEiEivJmPU1dXh8/nw+78e3/avqwOfD/x+P94w7pDaDzGvdEyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUkyMUno94a8F32Vb81cCO2OUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUYmKUolbMMmAwEAdkA51klQLwGjAyUH898I92ywVYAvQHrgLygMPtan4LTASuBhK6vOXhQaWYV4CFQCGwCxgH3AwEC3f8FzAduBt4H/jPQNvbquYp4A/A/wClwDcCY37equY8cDswLyx70U3chs9s3rxZbrvtNunfv78Asnr1alf9W0J+nJSqDlsWSEGr6UaQVJCiIPV3gExtNy8b5J7A6yaQFJD/brX8MxAPyP91MN5KEF8n29e6+VEQvQjQ0NDAuHHjWLZsWbiPEcA5astwLjXNRAamg2SVsq1dPThnQ3N9FVDTrsaHc4kMNublxvUX/m699VZuvfXWntgWAD4FGoH2WU7JwIEgfWqC1Ne0Ws4larTR49/EPHfuHOfOnWuZ7omw0q8jPX7zdx1WCkQB7bJKqQWCZJWScon6lFbzQh3zcqMvrBSYALTKKqUpMB0s3DGnXT3AO63qr8UR0LqmDufpLHyBkeGlxy9lXQorBfKBTCALeAZoAAK5oswCBgCBrFLmA7nA08BUoAR4D3ghsDwCWAD8BhiOI+oRIBXnsbqZo8CpwM9GYHdg/jAg3tUehIHuPNLRQ4/LAvIcyCCQ2MDj8/ZWy3JB8tvVvwoyIlA/GmRNu+VNII+AJAcek6eAHGxXkx949G3fNl6Gx2XX0YtnzpyhsrISgPHjx/O73/2O7373uyQmJjJo0KBL9m+JXgS+BsmL1OE8eoc7etH1GbNx48YOj6r8/PyQ+od6xnxVWk+dMa7vMTfddBPSu/mmVyQqPyszTIxaTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxS1IqxsFKFWFgp7rJkHn/8ccnMzJT4+HhJSkqSadOmyYEDB1xlpFhYaWi4OmM2b95MQUEB27dv55133uGLL77g+9//Pg0NDWE7UCys1MFVyM+6devaTBcXF9OvXz/KysqYPHlyWDbIwkodupXw5/f7AUhMTAxaY2GlXaPLN/+mpiYWLFjApEmTGDNmTNA6CyvtGl0WU1BQwN69eykpKem0zsJKu0hXnhgKCgpk4MCB8sEHH7juG8pTWUngiakYpAJkDkgCSE1g+UyQxa3qt4JEgywF2Q9SCBIDUt6q5onAGG+C7AGZBnItyL9b1XwI8j7IoyDxgdfvg9RfhqcyV2KampqkoKBAUlNT5dChQ11aoYWVhoYrMfPmzROfzyebNm2S48ePt7SzZ8+GPIZlYoaGKzEdHU2ArFy5MuQxTExouHpcdtwYvYHKz8oME6MWE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUE6MUtWIsrNQFy5cvZ+zYsXi9XrxeLzk5OaxduzbsG2VhpeAqS+att96SNWvWyKFDh+TgwYPy0EMPSUxMjOzduzfkMSysNDRciemIPn36yEsvvRRy/aXEnAOJAlndbv4skP8I0icN5Pft5i0BGRt4fSTw5r3frmYyyH1KxXQ5E7OxsZHXXnuNhoYGcnKC5+S5zcS0sFIH1zf/8vJy4uPj8Xg8zJ07l9WrVzNq1Kig9W4zMQ0H12IyMjLYvXs3paWlzJs3j/z8fCoqKoLWu83EtLDSAN29Fk6ZMkXmzJkTcn2oN/972938B1zi5n9bu3k5Hdz8l7a7N2i++XcrdxmcmN/W95BwsBDIBzKBLOAZoAG4K7B8FjAAKApMzwdygaeBqUAJ8B7wQmB5BLAA+A0wHCdV9hEgFeexupmjwKnAz0Zgd2D+MCA+TPsWMm4sLl68WDZv3ixVVVWyZ88eWbx4sURERMjbb78d8hgWVhoarsT8/Oc/l/T0dImNjZWkpCSZMmWKKykilokZKhEiIr15htbV1eHz+fAD3t5ccQ9Rh/MfOPj9frze8O2R2s/KrnRMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFJMjFLUirGw0m7wxBNPEBERwYIFC8K0OQ4WVoq79KXW7NixQwYPHixjx46V+fPnh9zPwkpDo0tnzJkzZ5gxYwYvvvgiffr0CeuBch4ow7nUNBMZmN4WpM+2dvXgnA3N9VU42Zeta3w4l8hgY15uuiSmoKCAqVOnkpfX/u34MufOnaOurq5N64zOwkqDBYt+HcNKXSf8lZSUsGvXLnbu3BlSfVFREY8++qjrDbvScXXGVFdXM3/+fFatWkVcXFxIfSystIu4uSGtXr1aAImKimppgEREREhUVJRcuHDhkmNYWGlouBJTV1cn5eXlbVpmZqb89Kc/lfLy8pDGCEVMSeBNKwapAJkDkgBSE1g+E2Rxq/qtINGBN34/SCFIDEh5q5onAmO8CbIHZBrItSD/blXzIU6a+aMg8YHX74PUaxfTEbm5uWF/XBYsrLTb8b49xb2B1hGbOph3e6AFIwJ4LNCCURxoGui2mE2bNoVhM4z2qP2s7ErHxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxCjFxChFrRjLxFSIZWKCq5CfwsLCLwXgZGRkuAqvsUzM0HB9xowePZrjx4+3tC1btoT1QLFMTAfXIT/R0dGkpPRcLl5nmZgHgvT5OmZiuj5jDh8+TGpqKkOGDGHGjBkcPXq003q3YaWGgysx2dnZFBcXs27dOpYvX05VVRXf+c53qK+vD9qnqKgIn8/X0tLS0jpdh2ViBujODer06dPi9XrlpZdeClrz+eefi9/vb2nV1dUh3fyv9EzMbgXJJSQkMGLECCorK4PWeDwePB6Pq3EXAvlAJpAFPAM0AHcFls8CBgBFgen5QC7wNDAVKAHeA14ILI8AFgC/AYYD1wKPAKk4j9XNHAVOBX42ArsD84cB8a72IAx0x2p9fb306dNHnn322ZD7WCZmaLgSs2jRItm0aZNUVVXJ1q1bJS8vT/r27SsnTpwIeYxQxXxVmopL2bFjx5g+fTonT54kKSmJG2+8ke3bt5OUlBTm89iIEBHpzRXW1dXh8/nwA97eXHEPUYfzy6rf78frDd8eqfyszDAxajExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSjExSlErxsJKFWJhpbhPXzp27JjMmDFDEhMTJS4uTsaMGSM7d+4Mub+FlYaGqzPm9OnTTJo0iZiYGNauXUtFRQVPP/00ffr0CduBYmGlDq7Sl5588knS0tJYuXJly7xrr702rBtkYaUOrs6Yt956i8zMTG6//Xb69evH+PHjefHFFzvtY2GlXcOVmA8++IDly5czfPhw1q9fz7x587jvvvv485//HLSPhZV2ETc3pJiYGMnJyWkz7xe/+IV8+9vfDtrHwkq7hqt7TP/+/Rk1alSbeddddx2vv/560D4WVto1XImZNGkSBw8ebDPv0KFDpKenh3Wj7gQ+wfmFsAa4AVjHxZv3UdpegycC/ws8DDyE8+b/FRjTquYBHLlzgM+AGwNjxrWqWQK0viiPD/zcCNzUjf3pEm5Orx07dkh0dLT89re/lcOHD8uqVavk6quvlr/85S8hj2FhpaHh+hfMv/3tbzJmzBjxeDwycuRIeeGFF1z1NzGhYWGl3cTCSq8wTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSTIxSuvV/LXeF5r9kf12+j9m8H+H+C32vizl58iQAnX8f86vHyZMn8fl8YRuv18UkJiYCcPTo0bDuSHvq6upIS0ujuro6rF+SaI/f72fQoEEt+xUuel1MZKRzW/P5fD36hjXj9Xp7ZT3N+xW28cI6mhE2TIxSel2Mx+OhsLDQ9RfNr7T19Po3MY3QsEuZUkyMUkyMUkyMUnpczKlTp5gxYwZer5eEhATuvvtuzpw502mfm266iYiIiDZt7ty5X6pbtmwZgwcPJi4ujuzsbHbs6DzY5LXXXmPkyJHExcVx/fXX849/tA826Rg36ykuLv7StsfFxQWtD0pY/7VNB9xyyy0ybtw42b59u7z77rsybNgwmT59eqd9cnNzZfbs2XL8+PGW1v4fBpWUlEhsbKysWLFC9u3bJ7Nnz5aEhASpra3tcMytW7dKVFSUPPXUU1JRUSEPP/ywxMTESHl5eafb4nY9K1euFK/X22bba2pqOl1HR/SomIqKCgHaRJqsXbtWIiIi5KOPPgraLzc3V+bPn9/p2FlZWVJQUNAy3djYKKmpqVJUVNRh/R133CFTp05tMy87O1vuueeesK5n5cqV4vP5Oh0zFHr0UrZt2zYSEhLIzMxsmZeXl0dkZCSlpaWd9l21ahV9+/ZlzJgxPPjgg5w9e7Zl2fnz5ykrKyMv72IISWRkJHl5eWzb1nEIybZt29rUA9x8881B67u6HoAzZ86Qnp5OWloa06ZNY9++fZ3ua0f06IeYNTU19OvXr+0Ko6NJTEykpiZ4WMhPfvIT0tPTSU1NZc+ePfzyl7/k4MGDvPHGGwB8+umnNDY2kpzcNoQkOTmZAwc6DjapqanpsL6z7ejKejIyMlixYgVjx47F7/ezdOlSJk6cyL59+xg4cGDQdbWnS2IWL17Mk08+2WnN/v37uzI0AHPmzGl5ff3119O/f3+mTJnCkSNHGDp0aJfH7Q1ycnLIyclpmZ44cSLXXXcdf/rTn/j1r38d8jhdErNo0SJ+9rOfdVozZMgQUlJSOHGibcrYhQsXOHXqFCkpoYeFZGdnA1BZWcnQoUPp27cvUVFR1Na2DSGpra0NOm5KSoqreqBL62lPTEwM48ePp7KyMqT6Zrp0j0lKSmLkyJGdttjYWHJycvjss88oKytr6fvPf/6Tpqamljc7FHbv3g04yRwAsbGxTJgwgQ0bNrTUNDU1sWHDhjZHa2tycnLa1AO88847Qeu7up72NDY2Ul5e3rLtIdPtx4dLcMstt8j48eOltLRUtmzZIsOHD2/zuHzs2DHJyMiQ0tJSERGprKyUxx57TN577z2pqqqSN998U4YMGSKTJ09uM25JSYl4PB4pLi6WiooKmTNnjiQkJLQ8ms6cOVMWL17cUr9161aJjo6WpUuXyv79+6WwsDDkx2U363n00Udl/fr1cuTIESkrK5Mf//jHEhcXJ/v27XP1vvW4mJMnT8r06dMlPj5evF6v3HXXXVJfX9+yvKqqSgDZuHGjiIgcPXpUJk+eLImJieLxeGTYsGFy//33dxhw8Nxzz8mgQYMkNjZWsrKyZPv27S3LcnNzJT8/v039q6++KiNGjJDY2FgZPXq0rFmzJqR9cLOeBQsWtNQmJyfLD37wA9m1a1dI62mNfeyvFPusTCkmRikmRikmRikmRikmRikmRikmRikmRikmRikmRikmRin/D78Lzy3djaZJAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_matrix(np.array([model.startprob_]).reshape((7,1)))"
   ]
  }
 ],
 "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.11.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}