Skip to content

Instantly share code, notes, and snippets.

@ryonsd
Created December 17, 2019 02:52
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save ryonsd/ced10ac16190897221c036b0c54b1f28 to your computer and use it in GitHub Desktop.
Save ryonsd/ced10ac16190897221c036b0c54b1f28 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"from scipy import stats\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"ガンベル分布は連続確率分布の一つで,極値分布のタイプI型に相当"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"母数:$\\mu \\in \\mathbb{R}, \\eta>0$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"平均:$E(x)=\\mu+\\gamma\\eta\\ (\\gamma=0.577\\cdots$, オイラー定数) \n",
"分散:$V(x)=\\frac{\\pi^ 2 \\eta^ 2}{6}$ \n",
"最頻値:$\\mu$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"確率密度関数"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$$f(x) = \\frac{1}{\\eta} \\left\\{ -\\left( \\frac{x-\\mu}{\\eta} \\right) \\right\\} \\exp \\left[ -\\exp \\left\\{ -\\left( \\frac{x-\\mu}{\\eta} \\right) \\right\\} \\right]$$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"累積分布関数"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$$F(x)=\\exp\\left[-\\exp \\left\\{-\\left(\\frac{x-\\mu}{\\eta} \\right) \\right\\} \\right]$$"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOydd3iUVfb4PyedNEpIAgklAUKHoCAWiqBSRNBVLKDYFXdXf991ZXfFtaDY1rbqrl2xstgbogKKoqJID72F0EJIJ6SROuf3xzuBmUlCCjOp9/M888zMPffe99yZ5D1z77n3HFFVDAaDwWCoLV6NrYDBYDAYmhfGcBgMBoOhThjDYTAYDIY6YQyHwWAwGOqEMRwGg8FgqBPGcBgMBoOhThjDYWhURORBEZlfz7Y3iMiKk8i/FZHrq6orIvki0qM+120qiMgoEdnZ2Hq4ExHZKiJjGlsPw8kxhsNQZ0Rkn4gcs99800TkbREJbmy9XFHVC1X1nWpkwaqaBGDX/5H6Xsfh88gTkRwR+U1E/igiHv3/UtVfVLWPix4X1KcvEYkREbV/pxWPje7TtsprVvrcVXWAqi735HUNp44xHIb6MkVVg4HTgWHAfa4VxKK1/I1NUdUQoDvwL+BuYF7jqlQv2tmNarCqxje2MoamSWv5pzZ4CFU9BHwLDAQQkeUi8qiI/AoUAj1EJEpEFopItogkisitLt0EiMiH9l/s60Xk+A1LRGaLyB67bJuIXOrSVkTkBRE5KiI7ROR8B8FyEbmlKr3tv657ichM4BrgH/Zf2V+JyN9F5FOX+v8Rkedr8XkcVdWFwFXA9SJS8bn4i8jTInLAPkt7RUTa2GVjRCRZRGaJSLqIHBaRGx2uPck+9jwROSQif3NsZ3/9HtAN+Mo+jn+IyNci8v9cxrGpis/wpLguJzrMTnzs75eLyMMi8qtdx6Ui0tGh/kj7LCxHRA7alw0rfe72usdnTfbP7DkRSbE/nhMR/9p8ZgbPYgyH4ZQQka7AJGCDQ/G1wEwgBNgPfAAkA1HA5cBjInKeQ/1LgI+BDsAC4AsR8bXL9gCjgLbAQ8B8Eens0PZMe52OwBzgMxHpUFv9VfU14H/Ak/Zf2VOA+cBEEWlnH6MPMA14tw79rraPeZS96F9Ab2AI0AuIBh5waNLJPsZo4GbgRRFpb5fNA26zz2gGAj9Ucb1rgQPYZ4Kq+iTwDjCjoo7dIEcDX9d2HHXgauBGIALwAyqMW3esHxb/BcKxxp9Qzefuyr3AWfY28cBwnGe2J/vMDB7EGA5DfflCRHKAFcBPwGMOsrdVdauqlmH9c48A7lbVIlVNAN4ArnOov05VP1HVUuDfQADWDQNV/VhVU1TVpqofAruxbiAVpAPPqWqpXb4TuOhUBqaqh4GfgSvsRROBTFVdV8euUoAOIiJYhvSvqpqtqnlYn9c0h7qlwFz7OL4B8oE+DrL+IhKqqkdUdX0tr78Q6C0icfb31wIfqmrJSdpk2mcGORUzm1rylqruUtVjwEdYN3uwDMr3qvq+fWxZ9r+B2nAN1meSrqoZWD8crnWQn+wzM3gQYzgM9eUPqtpOVbur6p/tN4wKDjq8jgIqbpYV7Mf6lVipvqraODE7QUSuE5GEipsZ1i/ujg5tD6lzpM79FW1PEcdf6zOA9+rRRzSQjfVLOxBY5zCOxfbyCrLshraCQqBiw8FUrFndfhH5SUTOrs3FVbUI+BCYYfc1Ta/FODrav9d2qvp0ba5jJ7Ua3btizQjrQxTW91mB63d7ss/M4EGM4TB4AscbecWv7hCHsm7AIYf3XSte2G9wXYAU+zLH68AdQJiqtgO2AOLQNtr+i96x75RT0LeCL4DBdh/FZKxllVojImdgGY4VQCZwDBjgcFNua99cULNyqmtU9RKsZaAvsH7R13Yc72D9cj8fKFTVlXUZh50CLMNXQac6tD0I9KxGVlNo7hSszQYV1Oe7NXgAYzgMHkVVDwK/AY+LSICIDMZaj3Y8uzFURC6z+xLuBIqB34EgrJtLBoDd+TnQ5RIRwP+JiK+IXAH0A76po5ppgNOZDvuv9U+wfC6rVfVAbToSkVARmYzl15mvqpvts6jXgWdFJMJeL1pEJtSiPz8RuUZE2tqX8nIBWx3GsdJe/xnqN2sCSABGi0g3EWkL3FOHtv8DLhCRK0XER0TCRKRiGauSvi68D9wnIuF2Z/sDOP/dGBoJYzgMDcF0IAbr1+LnwBxV/d5B/iXWLqQjWGvYl9nXrbdh3fBWYt1kBgG/uvS9CojD+lX/KHC5qmbVUb95WD6EHBH5wqH8Hfs1a3PD/UpE8rB+Yd+L5atx3OVzN5AI/C4iucD31H49/lpgn73dH7FmEFXxONaN1tU/8a59HPW66arqd1hLXpuAdcCiOrQ9gLXMNgtr2S4By9EN1X/uFTwCrLVfdzOw3l5maGTEJHIyGKpGRLoBO4BOqprb2PrUFxG5DpipqiMbWxdDy8DMOAyGKrD7Wu4CPmjmRiMQ+DPwWmPrYmg5GMNhMLggIkFYvoRxWGdDmiV2H0oG1jLfgkZWx9CCMEtVBoPBYKgTZsZhMBgMhjrh09gKNAQdO3bUmJiYxlbDYDAYmhXr1q3LVNVw1/JWYThiYmJYu3ZtY6thMBgMzQoR2V9VuVmqMhgMBkOdMIbDYDAYDHXCGA6DwWAw1IlW4eMwGAyep7S0lOTkZIqKihpbFUMdCQgIoEuXLvj6+tZcGWM4DAaDm0hOTiYkJISYmBicAxYbmjKqSlZWFsnJycTGxtaqjVmqMhgMbqGoqIiwsDBjNJoZIkJYWFidZorGcBgMBrdhjEbzpK7fm1mqMtRIRl4xP+/KICO/mIFRbTmnZxheXgLHciDxe8g5ABH9oed54OPX2OoaDAYPYwyH4aR8sPoAj3y9nfziExk647u05eWhKUSt+CcUZJyoHBYHl7wI3c5sBE0NBkNDYZaqDNXy9q97mf3ZZuK7tmXR/xtJwgPjeOaKePpm/0jkkpkUtukENy6Gu/fDVfPBVgrvXgy7v6+5c4PB0GwxhsNQJb8lZvLgV9sY3z+Sd24czsDotrQL9GNq1zz+xX/Y6tWH8UfuJrXdadCmHfSbArcsg4694aPrIH17Yw/BYDglFi9eTJ8+fejVqxf/+te/qq0XExPDoEGDGDJkCMOGDWsQ3Q4ePMjYsWPp378/AwYM4Pnnn6+yXm3HUGdUtcU/hg4dqobaU1BcqiOfWKbnPvmDFhaXnRCUlaq+Okb1iVjdk5Sk/e7/Vqe+9KuWldtO1DmaovpUnOp/TlctLmh45Q2NxrZt2xpbBbdRVlamPXr00D179mhxcbEOHjxYt27dWmXd7t27a0ZGRoPql5KSouvWrVNV1dzcXI2Li6ukX13GoFr19wes1SruqWbGYajEaz8ncTD7GE9eHk8bP+8Tgo3vQ8p6uPBJesTG8vAlA1m7/wgLVjnEQQvtDJe9DlmJ8POTDa+8odUzZswYduzYAUBWVhYDBw6scx+rV6+mV69e9OjRAz8/P6ZNm8aXX355yro98cQTXHvttUyYMIHu3bvz9NNP16ufzp07c/rppwMQEhJCv379OHToUIOMAYxz3ODC0WOlzFuxlwkDIhke2+GEoLwUfn4Kok6DgVMBuOz0aL5IOMSTi3cycWBnwkP8rbo9zoUhM+C3/8LgaRDRtxFGYmhMHvpqK9tS3Jtxt39UKHOmDKixXmJiIr179wZg06ZNDBo0yEk+atQo8vLyKrV7+umnueCCCwA4dOgQXbt2PS7r0qULq1atqvJ6IsL48eMREW677TZmzpxZrW6bN2/G29ubb775hqSkJKZOncrf/va3OuvnyL59+9iwYQNnnum8KaUuY6grxnAYnHj7133kFZXxl/N7Ows2fQg5+2HSU2Df8y0iPHTxAMY9+zMv/LCbhy5x+GU3/mHYvhCWzYXpJmupoWHYv38/0dHReHlZiymbNm1i8ODBTnV++eUXt15zxYoVREdHk56ezrhx4+jbty+jR4+usu7mzZv57LPP8Pb2xtvbmw4dOlSqUxf98vPzmTp1Ks899xyhoaH1HkNd8ajhEJGJwPOAN/CGqv7LRT4aeA4YDExT1U/s5WOBZx2q9rXLvxCRt4FzgaN22Q2qmuDJcbQWSstt/G/Vfsb0Cad/lMMfoSqsetU6qxE33qlNj/BgrjqjKwtWH+DmkT3oFhZoCQI7wIj/gx8egQOrzBbdVkZtZgaeYOPGjU6GYt26dVx11VVOdWrziz46OpqDBw8elyUnJxMdHV3lNSvKIyIiuPTSS1m9enWVhqO0tJSMjAx69uwJVD0bqq1+Ff1NnTqVa665hssuu6xKvWo7hjpTlePDHQ8sY7EH6AH4ARuB/i51YrCMxrvA5dX00wHIBgLt79+urm51D+Mcrx3fbErR7ncv0u+2pjoLDq5VnROquuq1KtulHj2mcfd+o7M/3eQsKM5XfaKH6ntTPaSxoSnRFJzjDz30kM6YMUNVVXft2qWhoaG6f//+OvdTWlqqsbGxmpSUdNyxvGXLlkr18vPzNTc39/jrs88+W7/99ltVVT3vvPM0OTn5eN3Nmzfreeedd/z9nDlz9PXXX6+zbqqqNptNr732Wv3LX/5yymOooKk4x4cDiaqapKolwAfAJS5Ga5+qbgJsJ+nncuBbVS30nKoGgAWrDxDVNoCxfSOcBeveBN8gGHxVle0iQwO4fGgXPl2XTHquQ7wbvyA464+Q+B2kbvGg5gaDxcaNG7HZbMTHxzN37lz69+/PO++8U+d+fHx8eOGFF5gwYQL9+vXjyiuvZMCAE7OoSZMmkZKSQlpaGiNHjiQ+Pp7hw4dz0UUXMXHiRGw2G4mJiU5LUZs3b3aaDSUkJBAfH1+vcf7666+89957/PDDDwwZMoQhQ4bwzTffOOlW0xhOiaqsiTseWDf8NxzeXwu8UE3dt6l+xvEDMNml7k5gE9Zyln817WYCa4G13bp1q9bKGizSc4s0dvYifXrJDmdBaZHqY11VP/vjSdvvzcjX2NmL9IlvtzsLCrNVH41S/fRWN2tsaGo0hRlHr169js8AGpPNmzfrX//615PW6dmzpxYWFjaQRjXTVGYcp4yIdAYGAUsciu/B8nmcgbWMdXdVbVX1NVUdpqrDwsMr5Vo3uPDtlsPYFKbERzkLEpdB8dHjO6mqI6ZjEOf3i+SjtQcpKXOYQLZpD0Ouhq2fQ0GWBzQ3GCzy8vIQEUJCQhpbFQYOHMi///3vauU5OTn4+fnRpk2bBtTKfXjScBwCujq872IvqwtXAp+ramlFgaoethvDYuAtrCUxwyny1cYUekcG0zvS5Z9uy6fQpoO1xbYGrj6zG5n5JSzdluosGHojlJfARrO7yuA5QkJC2LVrV2OrUSvatWvHtm3bGluNeuNJw7EGiBORWBHxA6YBC+vYx3TgfccC+ywEseIA/wEwi+enSHpeEWv2HWHyYJfZRlkx7FpshRPxrjkz2Oi4cKLbtWHBqgPOgsj+0PVMWPe2tUPLYDA0azxmOFS1DLgDa5lpO/CRqm4VkbkicjGAiJwhIsnAFcCrIrK1or2IxGDNWH5y6fp/IrIZ2Ax0BB7x1BhaC8t3WBFuL+gX6SzY/yuU5EOfSbXqx9tLuPrMbvy2J4ukjHxn4bCbrNPk+9y7h95gMDQ8HvVxqOo3qtpbVXuq6qP2sgdUdaH99RpV7aKqQaoapqoDHNruU9VoVbW59Hmeqg5S1YGqOkNVXe5QhrqybEcandsG0K+zyzLVriXgEwCxVR9mqoorhnXBx0t4f7XLrKP/JRDQDta+5QaNDQZDY9KkneMGz1NcVs6K3ZmM7RvhnAVM1TIcsaPBL7DW/UWEBDCufySfrEumuKz8hMC3jeUk3/4VFGS6cQQGg6GhMYajlbNu3xEKSso5r4/L2Y3sJDiyt9JJ8dpw5bCuHCks5eddLgbitBlWzo6tn5+CxgaDobExhqOV8+ueTHy8hLN6hjkLkpZbzz3Pq3OfI+M60j7Ql4UbU5wFkQMgvJ+1U8tgMDRbjOFo5fyamEV813YE+7uELdv7E4R2gQ496tynr7cXkwZ15rttqRQ4pJwFYNDlcGAl5BysurHBYGjyGMPRisktKmVTcg7nuM42bDbY+4t1dsPR71EHLhkSTVGpje+3pzkLBl1uPZtZh8HQbDGGoxWzZm82NoWzXQ1H+lY4ll2n3VSuDOvenqi2AXyZ4LJc1T4GupwBmz+pd98GQ0Nw0003ERERUWMiKI+lZz0JjZ061hiOVszvSVn4+Xhxerf2zoJ9K6znmFH17tvLS5gyJIqfd2VwpKDEWTjoCkjbDOk76t2/weBpbrjhBhYvXnzSOuXl5dx+++18++23bNu2jffff79BToT7+PjwzDPPsG3bNn7//XdefPHFStf1pG7GcLRi1uw7QnyXtgT4ejsLDqyEdt2g7anF7r84Pooym/LNlsPOggGXgnjBFjPrMLgfd6SOBRg9enSViZYcqWt6VpM61tCsOVZSztaUo9wyysX5rQoHfofYmmNT1UT/zqHEdgxi8ZZUrjmz+wlBcAR0H2Gd6TjvvlO+jqEJ8u1sSN3s3j47DYILa15ucUfq2NpS1/SsJnWsoVmzMTmH0nJlWHeXZaoj+yA/zS0Z+0SECQM68cYvSRwtLKVtoEO8q34Xw7d/h4xdEN67+k4MhjrQGKlj64JJHWto1qzbfwSAoa6G48Dv1nO3s91ynQkDInnlpz38sDONS0/rckLQ9yLLcOz4CsJnueVahiZELWYGnsBdqWNrS13Ss7ak1LHGcLRS1u0/QlxEMO0C/ZwFB1eBf6h1UM8NxHdpR6fQABZvSXU2HG2jIXoobF8Eo4zhMLiHhIQEioqsLJS7d+/myy+/5JFHnOOgunPGccYZZ7B792727t1LdHQ0H3zwAQsWWOkDzj//fN59993jN+udO3fSr9+J/6vqMgDWRj9V5eabb6Zfv37cddddddbtVDHO8VaIqrIpOYf4ru0qC1M2QNRp4OWePw0vL2H8gEh+2pXBsZJyZ2G/KZCy3hwGNLgNd6WOBZg+fTpnn302O3fupEuXLsybN++4rKb0rCZ1bAt4DB06tKpMia2WlJxC7X73In3nt73OgtIi1YfCVJc+4Nbrrdidod3vXqTfbj7sLMjYrTonVHXly269nqFxMKljT2BSxxpaHJuSjwIwKLqtsyBtqxWEMOo0t15veGwH2gX6snSrS2bAjr2sJbEdi9x6PUPrxKSObTiM4WiFbE4+io+X0K+zyy6MlA3Ws5sNh6+3F+f3jeT77WmUlduchf0mWwmjCrPdek1D68Okjm04jOFohWw6dJTekSGVD/6lbLDyi7fr5vZrjusfQW5R2fHdXMfpcyGoDRK/d/s1DQaDZzCGo5WhqmxOzmFwl7aVhSkJ1myjnoENT8bIuHB8vYUfdqQ7CzqfBkERVm5zg8HQLPCo4RCRiSKyU0QSRWR2FfLRIrJeRMpE5HIXWbmIJNgfCx3KY0Vklb3PD0XEz7VfQ/UkHznGkcJSBrkajtJjkL7N7ctUFQT7+3BmbBjLXA2Hl5eVLCrxeygv9ci1DQaDe/GY4RARb+BF4EKgPzBdRPq7VDsA3ABUtbn4mKoOsT8udih/AnhWVXsBR4Cb3a58C2bzIcsxPjjaZStu2lbQco8ZDoDz+kaQmJ7P/qwCZ0HvCVB01DpDYjAYmjyenHEMBxJVNUlVS4APgEscK6jqPlXdBNiq6sAVsZJinwdURMd7B/iD+1Ru+WxKPoqvt9C7U7Cz4LhjfIjHrn1+Pys9baXlqp5jwcvXynFuMBiaPJ40HNGA48muZHtZbQkQkbUi8ruIVBiHMCBHVSvSylXbp4jMtLdfm5GRUVfdWyybD+XQt1Mo/j5VOMaDwiHUPSEJqqJ7WBA9w4MqGw7/EIgZaQyHwdBMaMrO8e6qOgy4GnhORHrWpbGqvqaqw1R1WHh4uGc0bGaoKpuSj1b2b8CJE+MecIw7cn6/SFYlZZPvmlK290TI3AnZSR69vsFgOHU8aTgOAV0d3nexl9UKVT1kf04ClgOnAVlAOxGpiLFVpz5bOweyC8krKqt88K+kEDJ2QGfPLVNVMLZPBCXlNlbsznQW9B5vPe9a6nEdDAbDqeFJw7EGiLPvgvIDpgELa2gDgIi0FxF/++uOwAhgm/0I/I9AxQ6s6wH3ZCZpBWw/bEXc7O968C9jh3WWolPlSJ3uZlhMe0ICfPhhh0su8g49oGNv2G2WqwxNg9qmjo2JiWHQoEEMGTKEYcOGNYhuRUVFDB8+nPj4eAYMGMCcOXOqrNfsUsfa/RB3AEuA7cBHqrpVROaKyMUAInKGiCQDVwCvishWe/N+wFoR2YhlKP6lqhXHLO8G7hKRRCyfx4nIY4aTsjM1DxHoHekSkiHd/tFGuikA2knw9fbi3N7h/LAjA5tNnYW9J1hpa4srh5Q2GBqa2qSOreDHH38kISGBtWvXelgrC39/f3744Qc2btxIQkICixcv5vfff3eq02xTx6rqN6raW1V7quqj9rIHVHWh/fUaVe2iqkGqGqaqA+zlv6nqIFWNtz/Pc+gzSVWHq2ovVb1CVYs9OYaWxI7UXLp3CKSNn4tjPG0r+LSB9jENosf5/SLIzC8+vjX4OL0nQnkJJC1vED0MLZOGTB1bV9yVOlZECA62dkaWlpZSWlqKuPgnTepYg1vYmZpHn05VBIBL2woRfcHLu7LMA5zbOwIvgWU70p1Du3c9EwLaWqfI+01pEF0MnuGJ1U+wI3uHW/vs26Evdw+/u8Z6DZk6Fqyb+Pjx4xERbrvtNmbOnFltXXemji0vL2fo0KEkJiZy++23m9SxBvdTVFrOvqwCJsdHVRamb4O4CQ2mS4cgP07r1p7lO9O5a5xD2lhvX+h5Puz+zsp97uEdXoaWR2Okjl2xYgXR0dGkp6czbtw4+vbty+jRo6us687Usd7e3iQkJJCTk8Oll17Kli1b6j27qivGcLQSdqflY1Po6zrjyM+AggyIdD3U71nG9gnn6aW7yMgrJjzE/4Qgbjxs/QxSN0Hn+iW5MTQ+tZkZeIKGTh0LHM/wFxERwaWXXsrq1aurNBzuTh1bQbt27Rg7diyLFy92MhwmdazhlNmRmgtUYTjS7fsRGsAx7siYPhE8vXQXP+3K4PKhDille9n/MXYvNYbDUGcaOnVsQUEBNpuNkJAQCgoKWLp0KQ888ADg2dSxGRkZ+Pr60q5dO44dO8Z3333H3Xc7G2uTOtZwyuxMzSPA14vuYUHOgjT7LouIhjUcA6JCiQjx58edLqfIg8Mh6nRrucpgqCMNnTo2LS2NkSNHEh8fz/Dhw7nooouYOHGix1PHHj58mLFjxzJ48GDOOOMMxo0bx+TJk51082TqWDPjaCXsSM0jLiIEby8Xv0HaVivUSHDDnq4XEcb0CWfxllTKym34eDv8hokbDz8/aSV3CnTvrhZDy2bTpk2sX7/eLVkA33///WplFfm9wTJWrmzbto2pU6c6ZfibPn0606dPP/7+VHwSgwcPZsOGDTXqNmnSJCZNmlSva5wMM+NoJeyobkdV+laIaFj/RgVj+1jJnTYczHEWxI2zDiTu+aFR9DI0T0zq2IbDGI5WQFZ+MZn5xZX9G7ZySN/R4P6NCkbEdcTHS/jRNehh1GkQGGb5OQyGWmJSxzYcxnC0AnamWjs0Ks04juyDsmONNuMIDfBlaPf2/LjTJXqxl7flJE/8Hmy1irhvMBgaEGM4WgE7qjMcaRU7qhrHcACM7RvB9sO5pB4tchbEjYfCrBN5QgwGQ5PBGI5WwM7UPDoE+REe7O8sSN8GCIT3bRS9wPJzACx33V3V8zwQL7NcZTA0QYzhaAUkZuTTKyK4UiwbMnZCu27gF1R1wwagd2QwUW0DKm/LDewAXc4whsNgaIIYw9HCUVUS0y3DUYnMXRDep+GVckBEGNM3ghW7Mykpc/FnxI2DlPWQn151Y4PB0CgYw9HCycwv4eixUnqFuxgOWzlk7rZyYDQyY/tEUFBSztp92c6COHtyp8RlDa+UwWCoFmM4WjiJ6fkAlWccOfuhvLjRZxwA5/QMw8/bq/JyVafBENzJLFcZDE0MYzhaOIkZ1RiODPt+946NbziC/H04s0eHyttyRSDuAtizDMrLqm5sMBgaHGM4Wjh70vMJ9POmc9sAZ0HmTus5vPGXqsAKepiYns/B7EJnQdx4KDoKyWsaRzFDq6S2qVnBc+lZa0N5eTmnnXba8ThVDaWbRw2HiEwUkZ0ikigis6uQjxaR9SJSJiKXO5QPEZGVIrJVRDaJyFUOsrdFZK+IJNgfQzw5huZOYno+PcOr2lG1C4IioE37xlHMhbF9rFhZlbbl9hgDXj5mucrQoNQmNSt4Nj1rbXj++eedIu42lG4eMxwi4g28CFwI9Aemi4jrSbMDwA2Aa6zfQuA6eyrZicBzIuKQKo6/q+oQ+yPBIwNoIVS/o2pnk/BvVBDbMYjuYYEsd12uCmgLXc8y0XINtcYdqWNrk5oV6p6e1V2pY8HKr/H1119zyy23VClvrqljhwOJqpoEICIfAJcAx02equ6zy5z2YarqLofXKSKSDoQDLtHwDCcjr6iU1NyiyoZD1ZpxDLq86oaNgIgwtk8EH6w5QFFpOQG+Dmls48bB93MgNwVCq8hgaGhypD72GMXb3Zs61r9fXzr985811nNX6tiaUrNC3dOzujN17J133smTTz5ZZd366FYXPGk4ooGDDu+TgcqffA2IyHDAD9jjUPyoiDwALANmq2rxqSjaUknKKACgp+tW3Px0KD7apGYcAGP6hPP2b/tYtTebc3s7hHmPG28Zjt3fwdDrG09BQ5PHnaljPZGa1V2pYxctWkRERARDhw5l+fLlp6RTfWjS+ThEpDPwHnC9qlbMSu4BUrGMyWvA3cDcKtrOBGYCdOvWrUH0bWpUuxW3wjHeBM5wOHJWjzACfL34cUe6s+GI6AehXSw/hzEczYLazAw8gSdSx1aXmhXqlp7Vnaljf/31VxYuXMg333xDUVERubm5zJgxg/nz59dLt7riScNxCOjq8L6LvaxWiEgo8DVwr6oe90qp6mH7y2IReQv4W9ggxtEAACAASURBVFXtVfU1LMPCsGHDtG6qtwwSM/Lx8RK6hwU6CzIqdlS5d8ahquw6sosN6RvYn7ufY2XHaOPThq4hXYmPiKdfh354SfVutQBfb87uEWZ3kDuEehexlqs2fwxlJeDj51a9DS0Hd6WOrU1qVjh5elZPpo59/PHHefzxxwFYvnw5Tz/9tJPRqEm3U8WThmMNECcisVgGYxpwdW0aiogf8Dnwrqp+4iLrrKqHxfJU/QHY4l61Ww6J6fnEdAzC19vlZp25C/xCIKSzW65TUl7Cp7s/ZcH2BezL3QdAG582BPkGUVBawLGyYwBEBEZwWdxlzOg3g7b+bavsa2zfCH78cit7MwuI7egQQytuPKx7Cw7+DrGj3aK3oeWxceNGAgICiI+PZ/DgwcdTx95///116ufw4cNcf/31lJeXY7PZuPLKK522vE6aNIk33niDqKio4+lZy8vLuemmmxgwYECtU8dedNFFpz5oF2rSzS2oqscewCRgF5Z/4l572VzgYvvrM7B8HwVAFrDVXj4DKAUSHB5D7LIfgM1YBmM+EFyTHkOHDtXWyNinftTb3l1bWfD2FNXXxrrlGqsPr9ZJn07SgW8P1Ku/vlo/3vmxpuSlqM1mU1VVm82mh/MP68LEhfqn7/6kg94epOcsOEff2vyWFpUVVepvf2aBdr97kc77JclZUJSnOrej6pJ73aK3wf1s27atsVXQXr16aW5ubmOroZs3b9a//vWvJ63Ts2dPLSwsbCCNaqaq7w9Yq1XcUz3q41DVb4BvXMoecHi9BmsJy7XdfLtRqKrP89ysZoukpMzG/uxCJg2qYlaRuQt6jD2l/lWVVza+wssbX6ZrSFdevuBlRkaPrFRPROgU1IkpPacwpecUdmbv5Ln1z/HMumf4cs+XPDn6SeLaxx2v3y0skJ7hQfy4M52bRsae6Mg/GLqPsBzk4x+pdB2DwaSObTjMyfEWyoHsQsptSs8Il5DpxfmQdxg6xlXdsBaU2cr454p/8tLGl5jcYzIfT/m4SqNRFX069OHlC17mpfNfIrsom+lfT+fjXR871RnbJ4JVSdkUlriEGYkbDxk74Mj+eutuaLmY1LENhzEcLZQke4yq2I4uO6qyEq3nsF716rfcVs49v9zDoqRF3DHkDh4d+SiBvoE1N3RhVJdRfHrxpwyNHMrclXN5as1T2Owb58b2jaCk3MZviVnOjY5HyzWHAQ2GxsQYjhbK3kzrDEdsmMuM4xQMh6ry+OrHWbxvMXcNvYvb4m+r8jRtbenYpiMvnf8S0/tO591t7zL7l9mUlJcwLKY9QX7elaPlhvWE9jHmFLnB0Mg06XMchvqzN7OAsCA/2gb6Oguy9gACHWKrbHcyPtr5ER/u/JAbB97IjQNvdIue3l7e3DP8HjoFdeLZdc9SWFrIs2OeZUSvjizfmYGqnjBOItasY/17UFoEvgEn79zQ4Dh9X4Zmg+UHrz1mxtFCSXLdzlpBViK07Qq+dXPKbc/azhNrnmBk9EjuPP1ON2lpISLcNPAm7j/rfn5K/om7f7mbc/uEcSjn2PFDjMeJGw9lx2D/CrfqYDh1AgICyMrKqvNNyNC4qCpZWVkEBNT+h5iZcbRQ9mYWHI8460RWorXkUwfyS/KZ9dMs2ge057GRj530EN+pcGWfKykqK+KptU+hXX2A0fy4M524SIddMjEjwSfAWq7qVfmUr6Hx6NKlC8nJyWRkZNRc2dCkCAgIoEuXShtcq8UYjhZIXlEpGXnFlR3jqpbhiJ9Wp/4eW/UYKfkpvDnhTdoHeDYM+3UDrqOwrJAXE14kKkb5cUdHZo52MHS+bawDgLuXwoVPeFQXQ93w9fUlNrbuS6CG5odZqmqBHHeMuy5VFWRAcW6dHOO/pfzGV0lfccugWzg98nR3qlkttw2+jcviLiOvzWLWZy8lr6jUuULceMhOsvtrDAZDQ2MMRwukwnD0CK9uR1XtlqqKy4t59PdH6R7anVsH3+pOFU+KiHDfWffRr91QfDt9ynsbfnCuULFEZZI7GQyNgjEcLZCkjAJEoFsHl/MVddyK+/qm1zmQd4B7z7wXf29/N2t5cny9fHl1/HNIaRjzdj9MRqHDunmHWCuyrzEcBkOjYAxHC2RvZgFd2rdxToYEluHw9rN2VdXAwdyDzNsyj4t6XMTZUWd7SNOT075NO4YG3klJ+TFm/TSLUpvDklXceNi3AkoKGkU3g6E1YwxHCyQpM7+yYxwgMxE69AAv78oyF/6b8F98vXyZNXSWBzSsPRf1OY1jhy9lQ/oGnl/3/AlB3DgoL4G9PzeecgZDK8UYjhaGqrI3o4Ae1Z3hqMUy1Y7sHXy791tm9JtBeGAVW3obkHP7hFOWexoDgifyzrZ3+G6//dR4t7PBL9icIjcYGgFjOFoYGXnFFJSUV95RZSu3diLVwnA8v/55Qv1CuWHgDZ5Rsg5EhAQwKLotJemTGdRxEPf/ej/Jecng4w+x51qGwxw4MxgaFGM4WhhJ1W3FzTkAttIaDcea1DWsOLSCWwbdQqhfqKfUrBNj+4STcCCfB4Y/jiDc88s9lNnKrOWqowdOZDQ0GAwNgjEcLYzqt+LazzzUYDhe2fgK4W3Cmd53uifUqxdj+kZgU9id4st9Z91HQkYCb2x+wzIcYHZXGQwNjDEcLYy9mQX4+XgR1dYlFlUttuJuytjE6tTVXD/gegJ8mk4Awfgu7Wgf6MuPO9O5qMdFTIqdxCsbX2FTSTZEDDCGw2BoYIzhaGEkZeQTGxaEl5dLhNKs3eDfFoI6Vtv2zS1vEuoXyuW9L/ewlnXD20s4t3c4P+3MwGZT7j3rXiICI5j9y2wKeo6BAyuhKLex1TQYWg0eNRwiMlFEdopIoojMrkI+WkTWi0iZiFzuIrteRHbbH9c7lA8Vkc32Pv8jJoazEyeNituxlxWavKp2OUksO7CM6X2nE+RbRftGZmzfCLIKSth86CihfqE8NvIxkvOSeYpssJVB0vLGVtFgaDV4zHCIiDfwInAh0B+YLiL9XaodAG4AFri07QDMAc4EhgNzRKQiut7LwK1AnP0x0UNDaHaUlds4kFVIrKt/Aywfx0mWqd7a+hYB3gFc3e9qD2pYf0bHhSPC8eROwzoN48aBN/Jp6q/8FhpmlqsMhgbEkzOO4UCiqiapagnwAXCJYwVV3aeqmwCbS9sJwHeqmq2qR4DvgIki0hkIVdXf1Qr6/y7wBw+OoVmRfOQYZTatfIaj9BgcPVit4cgozGBR0iIujbuUDgEdGkDTutM+yI/Turbjx50nQo/8ecifiW0by5yO7chPNNtyDYaGwpOGIxo46PA+2V52Km2j7a9r7FNEZorIWhFZ21ryA1S7oyo7yXquJrjhx7s+ptxWzox+Mzyp3ikztk8Em5JzSM8rAsDf25+HRzxMupbxjF8xHE5oZA0NhtZBi3WOq+prqjpMVYeFhzfu6eeGYk+GlS2vUriRzN3WcxUzjtLyUj7e9TEjo0fSLbSbp1U8JcYNiEQVvtuWdrwsPjye63tfySehIaxMmNeI2hkMrQdPGo5DgGM0vS72slNpe8j+uj59tnj2ZhbQto0v7SvlGbdvxe1Qecbx/YHvyTyW2aTObVRHn8gQYsICWbI1zan8z2f8jRj1YU7azxSUmqCHBoOn8aThWAPEiUisiPgB04CFtWy7BBgvIu3tTvHxwBJVPQzkishZ9t1U1wFfekL55she+46qShvNsvZASBT4Vw58uGD7ArqFdGNE9IgG0rL+iAgTBnRi5Z5Mch2SOwX4BPBw90tIFRv/XjGnETU0GFoHtTIcIhIhIpeKyO0icpOIDBc5eeJpVS0D7sAyAtuBj1R1q4jMFZGL7f2eISLJwBXAqyKy1d42G3gYy/isAebaywD+DLwBJAJ7gG/rOOYWy97Mgsr+Dag2z/i2rG0kZCQwre80j+URdzfjB3SitFz5cUe6U/mQ027mutw8PjqwhN8P/95I2hkMrYOT5hwXkbHAbKADsAFIBwKwdjL1FJFPgGdUtcrTV6r6DfCNS9kDDq/X4Lz05FjvTeDNKsrXAgNPpndrpLCkjMNHi6qPitv/kkrFH+38iADvAC7pVVnWVDmtazsiQvxZsjWVS4Y47Ito15U7/LuzXHN58LcH+ezizwj0Day+I4PBUG9q+pk5CbhVVc9Q1Zmqep+q/k1VLwbisYzJOI9raaiRE3nGXZajCrPhWHYlx3hhaSHf7v2WCTETmkwww9rg5SWM6x/J8p0ZFJWWO8kC+k3hwdQUDuUf4oWEFxpJQ4Oh5VPTctPfVfVANbIyVf1CVT/1jGqGurC3uqi41cSoWrJvCYVlhVwWd1lDqOdWJgzoRGFJOb/sznQW9J3CsKJirmo3iPnb5rMxY2PjKGgwtHBOulRVgYi8B9yhqkft72OAeap6vudUM9SFvRmW4YjpWLs8458nfk5MaAynRZzWEOrVmvK8PIq2bKF4dyJlGemUZWSi5dbMwisgAJ/wcPqHhzM0P4dlGzoyrn/kicbhfSCsF3fm5LE8MII5v87hoykf4eft10ijMRhaJrUyHMAKYJWI3IV14O7vQOPmFDU4sTezgM5tAwj0c/lKsxLBywfadz9elJSTxIb0DcwaOqvyDqwGRlUp3r6d3KVLyV/2A8W7d58Q+vri07Ej4mttL7YVFlKelQWqPALYvheSPu5D8KhRBI8eRZvTT0f6TiZ45Qs8cM173L7ibt7Y/AZ/HvLnxhmcwdBCqZXhUNWKHU8/ApnAaaqa6lHNDHUi6WQ7qtrHgPeJsx2fJ36Oj/gwpeeUhlPQBVtJCblff0P2u+9SvH07eHkROGwY4X/5PwIGDyagb1+827dHvJxXU7WsjLK0NH5dspLvv/yZK33SKX7rLbJefx2fyEhCRw+lXQGMzs/loh4X8frm17mg+wX0bt+7kUZqMLQ8artUdS1wP9a5icHANyJyo6qaReQmgKqSlJHPlPioysLMRKeDf6W2UhbuWci5Xc8lrE1YA2ppoeXlHP3iSzKef56y9HT843rR6cE5hEyYgE/79jW2Fx8ffKOjGT7jD9yWFIT3GV154PwYCn7+maNfLSL786Vkl0UQtPNZ7pz1ICt9f2POr3OYP2k+3l7eDTBCg6HlU9vN+1OBkar6vqreA/wReNtjWhnqRFZBCblFZfQId9lRZbNB9h7oGHe8aGXKSrKLsrmkZ8NvwS3csIG9l03l8L334tOpE13feIPYhQtpP21arYyGI4F+PpzbO5xvt6RCYBChkybR9eWXiPtpOeEX9qE4JYcjf76L5z4KQtdvYv72+R4alcHQ+qiV4VDVP6hqusP71Vghzw1NgGqDG+YmQ1mRk2N8UdIi2vq3ZWT0yAbTz1ZcTNpTT7H/mhmU5+US/e9niPnwA4JHjjglH8uU+CjS84pZsy/7eJlPWBgd/zqbXpNT6XTrZIKzjjFngY2gWU+x9/fv3DEcg6HVc1LDISL32XNjVEJVS0TkPBGZ7BnVDLUlyR7csGd1wQ3tM46C0gJ+PPAjE2Mm4uvtEs/KQ5Ts38++K64ke96btLv8cnosXEjopEluccqf3y+CNr7eLNqU4izoMhxpH037zgfouXQJgbNup2u6jWM3/h+H58yhPCfnlK9tMLRmappxbAa+EpFlIvKUiPxDRB4QkfdEZDMwBVjleTUNJyMpowA/by+i2588z/iyA8soKi9ico+GsfV5y5ez9/IrKEtLo+trr9J57kN4B1eOl1VfAv18OK9fBN9uTqWs3CGli5cXDLgUEpfhZSuk+613cOjNf/LNMOHIx5+w58JJ5Hz6GWrydxgM9aImw3G5qo7Aije1FfAGcoH5wHBV/auqto5kF02YpMwCuocF4l0pz3gi+IVAsHXW4as9X9EluAvx4fEe1yn73XdJ/tOf8e3ahZhPPyV49GiPXGfK4M5kFZSwMinLWTDwMrCVwo6vAbg0/hq2zzibB28JQrpFc/jeezl4882UHj7sEb0MhpZMTYZjqIhEAddgRbZ9FSvr3hqgzckaGhqOpIz8qrfiZu4+nmc8vTCdVYdXMbnnZI+e3VBV0p/5N2mPPU7IBRcQs2ABfl1qm7+r7ozpE0Gwvw+LNroYgKjTrW3IW6zABiLCnLPnsDdC+c8tkUQ+OIfChI0kTbmYnM+/MLMPg6EO1GQ4XgGWAX2BtQ6PdfZnQyNTVm7jQHZh5RhVYI+Kay1Tfbv3WxTlotiLPKaL2mykPjCHrNdfp91VVxH93LN4BQR47HoAAb7ejOsfyeKtqZSUOSxXicCAyyDpJyiwQpN0C+3GHafdwfKUn1h9Zjt6fPkF/n37cPiee0i+4/9RduSIR3U1GFoKNcWq+o+q9gPeVNUeDo9YVe3RQDoaTkLykWOUlmvlGUdJoT3PuOUYX5S0iEEdBxHTNsYjeqgqqXPnkvPxx4TddhudHpyDeDfMuYnJgztz9FgpKxJdVk0HXgZaDttPpIG5pt81DAwbyOOrH6cgPJju77xDxD/+QcHPP7P3sqkcSzDpZw2Gmqjtdtw/eVoRQ/1IyrTvqKouz3jHXuw+spsd2Tu4qIdnZhuqStqjj5HzwYeE3Xor4Xf+pUFDmYyKCyc0oIrlqsiBluHc8tnxIh8vHx4a8RC5xbk8ueZJxNubsJtupPuCBYi3N/tmXEvWm2+ZpSuD4SQ0j+w9hmpJsgc37OG6VJVVkWc8jq+TvsZbvJkYM9EjOmS++BJH5s+nww03EH7XXxs8/pWfjxcTB3Zi6bY051DrIjBwKuxbAXknIuT0bt+bWwbfwqKkRfyc/DMAbQYNJPazTwkZO4b0J58k+fY7KD96tEHHYTA0F4zhaOYkZRbQLtCX9kEuEWAzra242qEHi/ct5qzOZ3kkxEjOp5+S+cILtP3DH4i4+x+NFjRx8uAo8ovLWL7TOTMgAy8DFLZ+7lR866Bb6dm2Jw///jD5JdaszTs0lOj//IfIf95D/i+/sPeyqRTt2NFAIzAYmg/GcDRzkjLyq8n6txtCo9mat49D+YeYEDPB7dfO/2UFhx+YQ9CIEXR+eG6jRto9p2cY4SH+fLr+kLMgvA90Ggwb33cq9vP246ERD5FWkMZz6587Xi4idLjuOmL+Nx8tK2Pf9KvJXby4IYZgMDQbPGo4RGSiiOwUkUQRmV2F3F9EPrTLV9nzfCAi14hIgsPDJiJD7LLl9j4rZBGeHENTJymjoHKMKji+o2rx3sX4ePlwXrfz3Hrd4qS9HLrrLvzj4oh+/vnjoc8bCx9vL/4wJIofd6STlV/sLBxyNRzeCGnbnIrjw+OZ0X8GH+78kHVp65xkbQYPJvaTjwno04dDd/6V9OefR202DAaDBw2HiHgDLwIXAv2B6SLS36XazcARVe0FPAs8AaCq/1PVIao6BLgW2KuqjttdrqmQO8bQam3kF5eRnldcOeufKmQmoh16smT/Es6JOoe2/m3ddt3yvDySb78d8fGh64sv4B1cxYynEZg6tAtlNmXhRpcQJIOusHKSbFxQqc0dQ+4gOjiaOb/NoaisyEnmEx5Ot3ffoe3lU8l6+RXL75Gf78khGAzNAk/OOIYDiaqapKolwAeAa0jWS4B37K8/Ac6Xyusd0+1tDS5UZP2rtKOqIAOKj7IxOJTUglS3OsXVZiPlb3+n5OBBop9/Dt9ozx3uqyt9O4UyICqUT9cnOwuCOkLceNj0EZSXOYkCfQN58JwH2Z+7n1c2vlKpTy8/Pzo//DCR991H/s8/s2/aNEqSD1WqZzC0JjxpOKKBgw7vk+1lVdZR1TLgKODqwb0KeN+l7C37MtX9VRgaAERkpoisFZG1GRktMypKxVbcSktV9hhVS0oz8fXyZUzXMW67Zta8eeT/9BOR98wmaPhwt/XrLqae3oUth3LZmZrnLIifDvlpkPRjpTZndT6Ly+Iu4+2tb7Mta1sluYjQYcY1dJv3BmXpGeybNo1jmzZ5aggGQ5OnSTvHReRMoFBVtzgUX6Oqg4BR9se1VbVV1ddUdZiqDgsPD28AbRuepIwCRKB7mEue8czd2IClR7YyInoEIX4hbrle4YYNZDz3PCETJ9L+6qvd0qe7uWRIFD5eUnnW0XsCtGkPCZWXqwBmDZtFh4AOzPltDqW20irrBJ11FjHvL8ArIID9111P7tKl7lbfYGgWeNJwHAK6OrzvYi+rso6I+ABtAcdoddNwmW2o6iH7cx6wAGtJrFWSlFlAl/Zt8PdxOaGdtZuEwGDSi7LctkxVfvQoh2bNwrdz50bfQXUywoL9Gds3gs83HHKOmOvjDwMvt4IeHqscVj3UL5R7z7qXHdk7eHvL29X279+zJzEffmA5zf9yJ1nz3jSHBQ2tDk8ajjVAnIjEiogflhFY6FJnIXC9/fXlwA9q/y8UES/gShz8GyLiIyId7a99gcnAFlop1lbcKnZUZSaypH0E/t7+blmmUlVS7r2XsvQMov/9DN4h7pnBeIqpp3chI6+YXxIznQVDpkN5caUzHRWc3+18xncfz8sbXybpaFK1/fuEhdHtnbcJmTCB9KeeIvXBh9CysmrrGwwtDY8ZDrvP4g6skOzbgY9UdauIzBWRi+3V5gFhIpII3AU4btkdDRxUVcf/YH9giYhsAhKwZiyve2oMTRlVZW9mQZVRccuzdrPUD0ZFjyLI99R3PB1ZsID875cRcdddtBk8+JT78zTn9Y2gfaAvn6xzWa6KOh069ql0psORe868hzY+bZjz6xxsWv32W6+AAKL//Qxht95KzocfcvCPfzI7rgytBo/6OFT1G1Xtrao9VfVRe9kDqrrQ/rpIVa9Q1V6qOtzRSKjqclU9y6W/AlUdqqqDVXWAqv5FVctphaTlFlNYUl758F95KesLU8ikzC2H/oqT9pL+5FMEjRpFhxuur7lBE8DPx4uL46P4blsaOYUlJwQi1qzj4KoT2RFd6NimI3cPv5uEjAQ+2HHyzXzi5UXErLvo9PBcClauZP/V15j8HoZWQZN2jhuqpyJdbKUdVUf2syTQnwDxYXSXU0uepOXlHL7nHiQggM6PPIJ4NZ8/l6vO6EZJma3yrCP+autMx9q3qm07pccURkSN4Ln1z5GSn1JtvQraX3EFXV97ldKUFPZdNY2i7dtPVX2DoUnTfO4EBif2ZNqDG7osVZVl7uS7oEBGhw8h0Dewqqa1JuvNNzm2cSOd7rsP38jmdUC/f1Qop3drx4JVB5yd1yGR0HcyJPwPSo9V2VZEeODsBwCYu3JurZzfwSNG0P1//wMvL/ZfM4P8X1a4ZRwGQ1PEGI5myt6MAtr4ehMZ4pwoaV3yL2R7ezOh55RT6r9o1y4y//NfQsaNI3Sy55I/eZIZZ3UnKbOAlXtc0sqecTMU5VTrJAeICo7iztPv5NeUX/lk9ye1ul5An97EfPgBvt26cfCPf+TIxx+fivoGQ5PFGI5mSmJGPj0jgvByyTO+JGM9bWzKqB4X1rtvLS0lZfZsvIKDrYRMTXTrbU1MGtSZdoG+zF+131kQM8rK07H2zZO2n9Z3Gmd3Ppun1jzFvqP7anVN38hIus+fT9A555B6/wOkP/uc2a5raHEYw9FMSUzLIy7CeVtsma2M74tTGSNBtPGpf0r4rHlvUrxtO50efBCfMPeHYm8oAny9uWJoF5ZuTSM91yEOlQgMuwmS18Dh6k+Ae4kXj4x8BD9vP2b/Mrvag4GueAcH0fWlF2l3xeVkvfoqKX//B7aSkpobGgzNBGM4miF5RaWkHC2iV4SzY3zN4TUcEWV8SK969128dy+ZL71EyIQJhE4Yf6qqNjpXn9mdMpvy4ZqDzoIh08GnDaydd9L2EYERPHj2g2zN2srLCS/X+rri60unuXMJv/NOchct4uDNt5jEUIYWgzEczZA99uCGcS6GY2niF7Sx2RgZNaJe/aoqqXMeRPz9ibz3n6esZ1MgtmMQI3t15P3VByi3OSwZtWlvZQfc9DEU5Z60jwu6X8Afev2BeVvmsT5tfa2vLSJ0/ONtRD31FIUJCey7+hoTINHQIjCGoxmyO80K4BcXeWKpqsxWxvfJPzGm8BgBkQPq1e/Rzz6ncPVqIv72N3wjmtcuqpMx46xupBwt4scdLhH4h90EpQWw6cMa+5g9fDZRQVHc88s95JXk1VjfkbZTJtPtjTcoy7AHSNzcaoMdGFoIxnA0QxLT8/Hz8aJr+xN+jNWpq8kpK2BCQSGE961zn2VZWaQ9+SRthg6l3RWXu1PdRuf8fpFEhvrz7u8uTvLo06FzPKyZZ+UwOQlBvkE8Pupx0grTeHzV43XWIejM4VaARD8/9l93HXk/Vo7SazA0F4zhaIbsTrfSxfp4n/j6lu5bSiDejLD5Q3DdZwtpj/8LLSyk89yHmtVBv9rg6+3FtWd15+ddGexKc5gtiMCZf4KM7bDnhxr7GRIxhJmDZ/JV0ld8mfhlnfWoCJDo36MHybffQfaCqiP1GgxNnZZ1h2gl7E7Pc3KMl9pKWXZgGeeqPwHhfa0bYh0oWLmS3EWLCJs5E/+ePd2tbpPg6jO7E+DrxZsr9joLBk6F4E6w8sVa9TNz8EyGRQ7j0VWPsidnT5318AkPp/t77xJ87rmkzX2YtCefMilpDc0OYziaGYUlZSQfOea0FXdN6hpyinOYcDQbwvvUqT8tKSH1kUfx7dqVsJm3ulvdJkOHID8uO70Ln204RKZjTnIfPxh+C+xZBuk1hwrx8fLhidFP0ManDbOWz6KwtLDOungFBtLlhf/S/urpZL/5JofumoWtuLjmhgZDE8EYjmZGUkYBqhAXeWLGsXTfUgJ92jAyJ73O/o3s9+ZTsmcPkf+8By9/f3er26S4aUQsJWU25rv6OobeZG3N/f2lWvUTERjB4yMfJ+loEo+vrru/A0C8vYm8/34i/v538hYv5sCNN1F25Ei9+jIYGhpjOJoZielWcMOKrbiltlK+P/A9YzoMxF+p04yjgTPnAwAAIABJREFUNC2dzBdfJHjMGELGjvWEuk2KXhHBnNc3gvdW7qeo1CGoclCYda5j44eQX7s0w+dEn8Otg2/li8Qv6uXvAGu7btjNNxH93LMUbdnC/mnTKTlwoF59GQwNiTEczYzd6Xn4eAndw6zghmsOr+Fo8VEm+EdaFeow40h/6im0rIzIf97jCVWbJLeMiiWroISPXaPmnnU7lJfAqtof8vtT/J9Oyd9RQejEiXR7+y3Kjx5l37TpHEtIqHdfBkNDYAxHM2N3Wj4xHYPw87G+uiX7lxDkG8SIYyXgFwKhUbXqp2D1asshfsvN+HXr5kmVmxRn9wjjtG7teGX5HkodU8t27AX9L4HVr1eZWrYqHP0ddy2/i4LSgnrrFXj66XR/fwFeQUHsv/4Gcr/7rt59GQyexhiOZkZier7TMtWyA8sY03UM/pm7Ibx3rXZUaWkpaQ8/gm9UFGG3tlyHeFWICHeM7cWhnGMsTHDJtTFqFhTnwpraJ5WMCIzgidFPsC93H//85Z8nzRpYE/6xsVY+8759OfR/fyH7nXfq3ZfB4EmM4WhGFJWWsy+r4PhW3NWHV1vLVN0nQMaOWi9THVmwgOLduy2HeJv6B0NsrpzXN4K+nUJ4aXkiNscwJJ0HQ9x4+P1lKKn97OGszmcxa+gsfjj4A69uevWUdPPp0MHKZ37BBaQ9/i9SH3nU5DM3NDk8ajhEZKKI7BSRRBGZXYXcX0Q+tMtXiUiMvTxGRI6JSIL98YpDm6Eistne5j/SXGN+14PE9HxsCn06WVtxl+yzlqnOadcH8tMgon+NfZRlZJDx3xcIGjmS4PPP97TKTRIR4faxvdiTUcCSranOwlF/g8IsWFe3X/vX9r+WyT0m81LCS/xwoObDhCfDKyCA6OeepcP113Nk/nwOzryN8pzaLZ8ZDA2BxwyHiHgDLwIXAv2B6SLieme7GTiiqr2AZ4EnHGR7VHWI/fH/27vv8Kiq9IHj3zM1vfcEUiAgHQQBKYKgAtIUFHFVEF0RK6u/ddeydl13LaCirIIVRQURFBVEmiKgQJAaaYGASUhCes/U8/vjDiRgEAKZTEjO53nmmTv3nnvnvQyTd+49bVqt9f8D7gCSXY/h7jqHpmZPtjYYX4fogBO3qS5vdTnm/P1agajOZzzGsZdfwWmxEPnYoxfsPBsN4eou0SSF+fLa6gMnX3W07qPN17HhVbCefR8NIQRPXvoknUI78chPj5xXZTm4mus+8jDRzz9HxZYtpE+4AcuBuudJV5TG5s4rjt5AmpTykJTSCnwGjD2lzFjg+E+7RcDQP7uCEEJEAwFSyl+kNjvOPOCahg+9adqTXYaXUUdCqC+bsjdRai1lWMIwyHENmhf554mj8tdfKfnqK0KnTMGcmNgIETddep1g+hXJ7M0p49td2SdvvPxR7Qpuyzv1OqaXwYtXL38VL4MX09dOp9T656Puno2g8eOJ//BDnJWVHL5hImWrV5/3MRXlfLkzccQCtSdByHStq7OMlNIOlADHZw5KFEJsE0L8KIQYWKt87XaUdR0TACHEVCFEihAiJS/v7NrmN3V7sktpH+mPXidYcXgFfkY/+sX0g9zd2rAZvmGn3Vfa7eQ88yyGqCjCpt3ZiFE3XaO7xtA+0p+Zq/Zjr93CKr4ftBkK62eCpX4j4Ub5RjFj8AyyyrL4x4//wO48//oJn4t7kLjoc0yuMa7yZs9WswoqHtVUK8ezgdZSyh7Ag8AnQoiA+hxASjlHStlLStkrPDzcLUE2Jikle3NK6RAdgMVhYdWRVQxtPRST3qQljjPcpipasADL3r1EPvxPdD4+jRR106bTCR64sh2H8ir48tQWVkMeg6pC+OWtunf+Ez0je/JY38fYcHQDz296vkH+yBujooj/+CMCxowm//VZZP3tAZwV5978V1HOhzsTRxbQqtbrONe6OssIIQxAIFAgpbRIKQsApJRbgYNAO1f5uDMcs1nKLbVQVGmjQ3QAP2b8SLmtnJFJI8Fhg7x98CdzcNgLC8l77XV8Lu2L/7BhjRh10zesUyRdYgN5ddV+LPZavclje0L7kbBxFlQW1vu417W7jts638ai/Yt4P/X9BolV5+VFzH//qw1TsnKlNjGU6mmueIA7E8cWIFkIkSiEMAETgaWnlFkKTHYtXweskVJKIUS4q3IdIUQSWiX4ISllNlAqhOjrqguZBJzbeA8XmOMV4xdF+bMsfRlh3mH0juoN+fu1Hs+RXU6777FXXsFZWUnUv/7VoivE6yKE4KFh7cksquKjn08Zw2rIY2Atg3Uvn9Oxp188nWEJw5i5dSYrDq9ogGhrhilp9fbb2LKzSR9/nar3UBqd2xKHq87iXmAFsAdYKKVMFUI8I4QY4yr2LhAqhEhDuyV1vMnuZcBOIcR2tErzaVLK4z/77gbeAdLQrkSWu+scmpI9OVriiAmFdZnrGJ4wHL1OD7mpWoHT3Kqq2r6dki8WEzJpUrMdMv18XdYunMvahTNrTRrFldaaDZGdoPtNsHkOFNS/lZRO6Hh+wPN0D+/Ooz89Wq9pZ8/Eb+AAEhd/gal1azLvuZdc1/AxitIY3FrHIaVcJqVsJ6VsI6V83rXuCSnlUtdytZTyeillWyllbynlIdf6L6SUnVxNcS+WUn5d65gpUsrOrmPeK1tILeGe7DJig7zZnPsDNqdNu00FkLML9CYIbfuHfaTDQc6zz2GIiCDs7rsbOeILy2NXd6Cs2sbrq9NO3jDkX6A3wqqnzum4Zr2Z14e8TrRfNPeuuZd9hfvOP1gXU1wc8Z/MJ2jiDRS++x6/3zoF27FjZ95RUc5TU60cV06xJ1urGF+Wvoz4gHg6hbrqNHJ3az3G9cY/7FP8+edUp6YS8Y9/oPfzbeSILyzto/yZ0KsVH/1ymMP5tSqd/aOg/3TYsxSO/HxOxw72CmbOlXPwNngzbdU0MkozzrzTWdKZzUQ/9RQxL/6XqtRU0seNp2LT5gY7vqLURSWOC0C1zcGhvHJaR1jZkrOFkYkja+oqclPr7L9hLyzk2MxX8enTh4CRVzdyxBemB69qh0mv49lvfjt5Q7/7wD8avvsnOB1173wGMX4xzLlyDjanjakrp5JX2bBNxAPHjCFx4QL0/v78PmUKeW++qW5dKW6jEscF4ECuNtRIhTEFieTqJFciKM/TOqrVUb9xbMYMnBUVRP3rMVUhfpYi/L2YfkUyq/ceY9VvuTUbTL5w1XOQvQNS3jvn47cJasPsobMpqC7gzlV3UmIpaYCoa5iTk0n4/HMCRo0kf9YbHJl8K7ajR8+8o6LUk0ocF4DjLar2lv1I59DOxAfEaxtyd2nPp1xxVG3fTsmiL7QK8eTkxgz1gjelfyLJEX489XUqVdZaVxedx0PiZbDm2bOe7KkuXcO78trlr3Gk5Ah3fH9HgycPvZ8vsS++SMyL/8Wydy+Hxl5D6fIW0X5EaUQqcVwAdh8twdc3n0Ol+2uuNkD7BQwQVdMUVzocZD/zjKoQP0dGvY5nxnYms6iK2T/UqigXAq5+RRu/auUT5/Uel8ZcyquXv0pacRpTV05t8OQBrltXXy7BlJRI1gMPcvTRx1SHQaXBqMRxAdiRUUxEzC70Qs+IxBE1GzJTICQJfEJOrCpasADLb3uIfPifqkL8HF3aJpRrusfw1o8H2Z9ba8iR8HbQ717Y8Qkc+uG83mNg3EBevfxV9hftZ9rKaQ0yrtWpTK1akfDxx4TeNY2SJUtIHzeeqp07G/x9lJZHJY4mrtrm4LfsIspNmxgYO5Aw71rjUWX9qvVwdrEXFpL36mv49O2L/4gRdRxNOVuPj+qIv5eRhxbtxFF79NxB/4SQNrD0/nrN2VGXy+IuY+bgmewt2sud399JcXXDD50ujEYipk8nft6HOK1WDk+8UesQarE0+HspLYdKHE3cnuxSnN77qXYWc01yrYGAS49C2dGTEseJHuKPqx7i5yvUz8xTYzqxI6OY99an12wwesPYN6D4CKx+9rzfZ3CrwcwcPJP9RfuZsmJKg7e2Os7nkktIWvoVgeOupWDuO6SPV1cfyrlTiaOJ25FRjDFwC4GmYC6Lu6xmQ5arF7IrcVRu20bJF4sJvXWy6iHeQEZ3jebKjpG8/P0+DuWV12yI7we9p8Kmt865b0dtg1sNZvYVs8kqz2Lyd5PJKnfP8Gt6f39innuOVnPn4CyvcF19zMBptZ55Z0WpRSWOJm5LRiZG/72MaTsKo65WJ7+sFNAZIKprzZDpkZGE3XWX54JtZoQQPHdNZ7yMev62YDu22kOvD30SglrDkjuh+vzrJ/pE9+Gdq96hxFLCpOWTznsiqD/jN3AgSV8vdV19zCV93Diqduxw2/spzY9KHE3c1oLVIBxc2/bakzdkbdWa4Rq9KPxwHpY9e4h85BF0vqpCvCFFBnjxn3Fd2JlZwsyV+2s2mP1g3FwoyYBlDzXIe3UN78r7w9/HKZ3csvwWtuRsaZDj1qWuq4/sp5/GUdrwlfRK86MSRxNWXGGl1LCecFNbkoNr9cdwOiFrG8T2xJqZSd6sWfgNGYL/sKs8F2wzNqJLNBN6xfG/Hw/yy6GCmg2t+8Bl/4Cdn8GuRQ3yXu2C2zH/6vmEe4czdeVUvjn0TYMc93T8Bg4k6ZuvCZl0C8ULFnLw6pGUfP2NmihK+VMqcTRhi/esR++Vy/DWp8yOm78frGXImIvJeepphE6nKsTd7MnRnYgP8eFvn20nv7xWi6TLHoK43vDNA+c0gm5dYvximDdiHj0ievDIT48wZ+cct/4h1/v5EfnIIyQu+hxjdDRHH3qI32+7DUt6+pl3VloklTiasK8PLUI6vJjcrY7bVEDpnkoq1q8n/MEHMUZHeyDClsPXbOCNv1xMUaWV+z/dVtNEV2+A694FnR4W3HLeTXSPCzQH8tYVbzEyaSSzts3i4Z8eptpe3SDHPh2vjh1J+OxTIp94nOpdu0kfM5Zjr72mOg4qf6ASRxNVWF3Iwaqf8bH2JdLvlFlzs7Zid/qT+8YHeHXrSvCNEz0TZAvTOTaQZ6/pzMaDBcxYWWt49KDWMP5dOPabduXRQFcHJr2JFwa8wP097md5+nImLZ9Ednl2gxz7dIReT8hf/kKb5cvwHzaMgv+9xcHhIyhevATpdJ75AEqLoBJHE7XkwBIkdnoE1dGRL3MLuakxOMrKiH7mWYRe3/gBtlATerXixt6teHPtQZbvqvVHvO1QuPxR2LlAa6bbQIQQ3NH1DmYNmUVGWQYTv53o1krz4wzh4cS+/BLxn36CISaa7Ecf5fB111O5xf3vrTR9KnE0QU7p5NM9C7FXJHF50ikj31aXULb1AKWpZYRNnYpX+3aeCbIFe3J0J3rGB/PAwu3szKzV23vg3+GiUbDiUdjfMFPFHjeo1SDmj5xPgCmAO76/g3d2vYNTuv8KwKdHDxI+/ZSYl17CXlTEkVsmkXnf/ViPHDnzzkqzpRJHE7Qucx25VUexFfWld2LISdvsu1aTvTkAc5tWhN01zUMRtmxeRj1v39KTMD8zt3+YwtHiKm2DTgfj5miDTi66DXJ2N+j7JgUm8enIT7ky/kpe+/U17lp1FwVVBWfe8TwJnY7A0aNos+xbwqffT/mGDRwcOYrsJ57ElpPj9vdXmh63Jg4hxHAhxD4hRJoQ4uE6tpuFEAtc2zcJIRJc668UQmwVQuxyPQ+ptc8PrmNudz0i3HkOnvBB6gd4EUYQF5MYVtMvQ0pJziv/w2nTEfPflxDGP876pzSOMD8z7916CdVWB1Pe31IzV7nJF278DMz+8MkEKG642f4A/Ex+vHjZizxx6ROk5KRw/dfX80v2Lw36Hqej8/Ym7K67aPPdcoInTKB4yRIOXjWM3BdewJ6f3ygxKE2D2xKHEEIPvAmMADoCNwohOp5S7HagSErZFpgJ/Ne1Ph8YLaXsAkwGPjplv5tc85F3l1I2q0mWd+fvZmvuVpzFA+iTGH5SE9vSb5dRtj2TsP6BeHXu5sEoFYB2kf68fUtP0vMruO2DLVRaXTPuBcTATZ+DpRw+ugYqGvaPqhCC69tdzycjP8HX6Msd39/Bfzb/hyp7VYO+z+kYIyKIeuJx2ixfTsDoURR+PJ+0K6/i2IyZOIobfqBGpelx5xVHbyBNSnlISmkFPgPGnlJmLPCha3kRMFQIIaSU26SUx6cuSwW8hRBmN8baZHyY+iG+Bj8KcnvQp9ZtKlvuMXKeeQbvUBuh41VHv6aiX9swXr+xO9szipn28a9Y7K7Jn6K6wF8WQEkmfDwOqht+zo32Ie1ZOHohN3W4ifl75jPh6wnszm/Y22N/xhQXS8zzz5P0zdf4DxlCwdy5pA29gmMvv4w9zz2DNSpNgzsTRyxQ+zo907WuzjJSSjtQAoSeUmY88KuUsvY40O+7blM9Lk7T600IMVUIkSKESMm7QP4TZ5Vn8f2R7+kcMAycZvq31YZQlw4HR//+d6Slmug+hYjkIWc4ktKYhneO5oVxXVi3P4+7Pv6VapsrecRfChPmafPCf3QtVDX8r3FvgzcP936YuVfNpdpRzc3LbublLS9Taats8Pc6HXNiIrGvvEzil1/iN2gQBe+9T9rQK8h+8imsv//eaHEojadJV44LITqh3b66s9bqm1y3sAa6HrfUta+Uco6UspeUsld4eLj7g20A81LnoUOHo6g/sUHeJLnqN/Jn/4/KLVuIGtcRc4gRWl/q4UiVU91wSWv+fW0X1uw9xp0fba1JHu2Gackje6d226qqyC3v3ze6L1+M+YJxyeP48LcPGfvVWH7I+MEt73U6Xu3bETvjFdosX0bgtddSsngxB4ePIOv//k713r2NGoviXu5MHFlAq1qv41zr6iwjhDAAgUCB63UcsASYJKU8MZaDlDLL9VwGfIJ2S+yCl1uRy6L9ixiZNIqUQ04GtA1DCEHFL5vInz2bwLFjCQpNg/j+YPTydLhKHf7SpzX/GdeFdQfyuPX9zZRU2bQNF42EGz7Srjw+GKXNpeIGAaYAnrj0CeaNmIef0Y/71tzHA2sfIKeicVs+meLjiX76KdqsXkXIlFspX7uW9Guu5cgtkyj9/nuk3d6o8SgNz52JYwuQLIRIFEKYgInA0lPKLEWr/Aa4DlgjpZRCiCDgW+BhKeWG44WFEAYhRJhr2QiMAhrvpq4bHW+XPyDsRsqq7QxsF4a9oICjDz2EKSGBqL/dro1R1UbdpmrKJvZuzcwJ3dl6pIgJb/1c01S3/QittVXRYXjnSji2x20x9IjowcJRC5l+8XTWZ61nzJdjeGvHW41WeX6cMSKCyIceou3aNUQ89HdsWVlk3T+dtCuvIn/OXOxF7rn6UtzPbYnDVWdxL7AC2AMslFKmCiGeEUKMcRV7FwgVQqQBDwLHm+zeC7QFnjil2a0ZWCGE2AlsR7timeuuc2gsORU5fHHgC65JvobdR/ToBPRLCOboP/6Jo6SE2FdnostcpxVuO9SzwSpndE2PWD6Y0pujxVWMm72RvTmuocrbDoUpy8Bpg/eGweH1bovBqDfy1y5/ZfHYxfSP6c+b299k1JJRfJX2VaN0HKxNHxhI6O2302bl98S9+QamhHjyZswgbfDlHH30MSq3bVOj8V5gREv4wHr16iVTUlI8HcZpPfvzsyxOW8yya5dx69wDBPkYmVW+kYK57xD19NME3zABPh6vjb56/zZQo+BeEPZkl3Lr+5uptDiYeUN3rugYqW0oOgLzr9OuPkbNhB43uz2WrblbeXnLy+wu2M1FIRfxYM8H6Rvd12MjKlsOHKBw/nxKln6NrKzE1KYNQePGETh2DIawMI/EpPyREGKrlLLXqeubdOV4S3Ck9AiL0xYzPnk8Nksg+3LLuLFsDwVz3yFo4g1a0rCUQfo67V65ShoXjA7RASy5uz+tQ33467wUXlqxVxtVNzgeblsBrfvCV/fA0vvA5t6Rb3tG9mT+yPm8eNmLlFpKmbpyKrd+dysbj270yK99c3Iy0U89RfK6dUQ/9yx6f3+OvfQSBwZfTsa991K2di3SZmv0uJSzo644POy+1fexJXcL31z7DV+mlLLg4+95/ee38O7Smfj330OYTJC6BD6/FW5dBgn9PR2yUk/VNgdPfpXKgpQM+rUJ5fUbexDmZwanA9Y+Dz+9AlFdtdZXIYluj8fisLD4wGLe3fUuuZW5dAvvxl3d7qJfTD+PzuliOXiQ4i8WU/LVVzgKCtAHBeE/bBgBI6/Gp1cvhE79zm1sp7viUInDgzZkbWDaqmk82PNBpnSewu0vfMmdC/9NSLA/CQsX1FyyL5ys3Q//v33a/A/KBWlhSgaPf7mbAG8jL1zbpebW1b7vYMlUkMDwF6D7XxrlytLqsPJl2pfM3TWXnIocOoV2YlLHSVyZcOXJ89s3MmmzUf7TT5R+86125VFVhSEykoARIwgYeTVenTurScsaiUocTSxx2Jw2xi8dj8PpYMnYJWQfzmP/hIlEymraL1qAOSlJK1hdCi8nw8WT4OqXPBu0ct5+O1rKgwu3szenjHEXx/LkqE4E+hi1+o4ld8HvGyF5GIx+DQIaZ3Ium8PGlwe/ZF7qPA6XHibCJ4IbL7qR69tdT6A5sFFiOB1nZSVla9ZS+u23lK9fDzYbhpho/Idegf/Qofj06okwqB9T7qISRxNLHO/vfp8ZW2fwxpA3GBDcky3X3Yz37wcJnP02CYP71RTc/il8OQ1uXwmtmkWXlRbPanfyxpoDvPnDQcL8TDw9pjPDOkUipITNb8Oqp8Fggqueh+43aaPuNgKndLI+az0f/fYRv2T/gpfei5FJIxmfPJ7OYZ7/le8oKaFs1SrKVq2mYuNGpMWCPjAQv8GD8btiKL6X9kPv53vmAylnTSWOJpQ4DhYfZMLXE+gf25+Zvf9NxtQ7Kd+2nQWj7ubZl+45ufCHY7Rfo9N3qIrxZmZXZgkPLdrB3pwyBrQN48nRHUmO9Ndaz311D/z+M8RcDCNehFaXNGps+4v2M3/PfJanL6fKXkVycDLjk8czKmmUx69CQLsSKd+wgfJVqyn74QecJSVgNOLTvTu+AwbgO6A/Xh06qHqR86QSRxNJHDanjVuW3UJWeRaLr/qEyumPUbH1V/7d8y9c87dJXNczrqZwwUGYdTEM+Rdc9pDnglbcxu5w8vEvR5ixcj8VVge39I3nviFtCfUxwq6FsPJJKM+Bbjdq/w8C48580AZUbi1n+eHlLN6/mN0FuzHpTAxqNYgRiSMYGDsQL4PnRzGQdjuVW3+lYv1PlK/fgGWP1rlSHxKCb79++A7oj++l/TBGNrsZGNxOJY4mkjje2vEWb25/kxk9nqbt8wuo3rWb78dO4z2v9mx6dChexlrTwK54TJuG9IFU8I/yXNCK2xWUW3hl5X4+3fw73kY9t/ZL4I6BSQQbLFqrq5/f1ApePBkGPqgN3d7I9hXuY0naEr5L/46C6gJ8jb4MaTWE4YnDuTT6Uoz6pjE/jD0vj4qNGynfsIGKDRtxFGiTXRnjW+PTqxc+vS7B55JeGGNjPX77ralTiaMJJI5N2Zu4c+WdXBMwkFvmpmPLyMD7uRcY+gtM7pfA46NqTVdiKYOZnSFpkNZMU2kR0o6V8drqNL7ZeRRfk4FJl8Zza78EIpzHtASy7WMQeug5GfreBSFJjR6j3WknJTeF5enLWXlkJWXWMnyNvvSL6cfgVoMZGDuQYK/gRo+rLtLpxLJ3LxWbNlOZkkJVSgqOEm2Ie0NkpCuR9MS7WzfMyclqcrRTqMTh4cSRWZbJxG8n0uOYL/d+XgmVVcTNns1Lub58/MsRfnhoMHHBPjU7bHgNVj4Bf10NcX/43JRmbn9uGa+tPsCyXdkYdIJRXWO4rX8iXXyLYd1LsONTrR9I+6u1BJIwwCN1YDaHjY1HN7I2Yy3rMteRV5WHQNAtvBuDWg1iYOxAkoOT0YmmUdcgnU4saWknkkjFli048rSJtoTZjFfHjnh37YJX1654d+2KMS6uRV+VqMThwcRRaavk5uU3k7jxCH9dZsMYEUncm29SFhPPgP+uYUy3GF66vtaMfrYqeLUrRHaESV95LG7F844UVPDhxiMsTMmg3GKnW6sgru8Zx5gkQcDueZDyHlQWQHgH6HETdL0B/DxzL98pnewp2MMPmT/wY8aP7CnU6hqCzEFcEnUJvaN60zuqN4mBiU3mj7GUEltmJlU7d1K9cxdVu3ZRnZqKtGjT/+iDgvDq2AFz+4swt2+H10UXYU5K0jrmtgAqcXgocVTaKnlg2V30+DiFQbuc+PTtS+zMGRiCg3lk8S4+T8lg5YODTppbnB9fgrXPwZTlEN/v9AdXWoyyahufp2SyMCWDvTllmA06hnWKYnSnYAZV/4Bp58eQuQV0Bki+CjqP1569AjwWc25FLptzNrMpexObczaTXZENQJh3GD0je9I1rCvdIrrRIaQDJn3T+UMsbTYsBw5QtXMXVbt2Ytm7D8uBA0ira155gwFzUhLm9u3xuqg95nbtMCUmYYyJbnatuFTi8EDiqLJX8cKcSQyZt5uIUkH4tGmE3X03wmDgt6OljJr1E5P7JfDk6E41O5UehVk9oe0V2hwOilKLlJLdWaV8vjWDpTuOUlxpw9uoZ3D7cK6Lr2BA+QrMqYu0llh6EyQNhg6jtVtavp4bPFBKSWZ5JpuzN7M5ZzPbj23naIU2L4lRZ6RDSAe6hnelU1gnLgq+iITABAy6ptOxT9rtWI8coXrvXiz79lO9T3u259TMdSK8vDDFx2NKSsScmIQpMVFbTkhA53th9i9RiaORE0deVhprHplC18352CKDaTtjFj49ewJaB7BrZ28gu6Satf83WOs5DCAlzL8eDv8E92yC4IRGjVm5sNgcTjanF/Ld7hxWpOZwrMyCQSe4uFUA4yKyGez8hcis7xHFv4PQQXQ3LZEkDYZWfcDo7dH48yrz2Jm3kx35O9iZt5MFMK5DAAAN7ElEQVTU/FSqHdpgjyadiTZBbWgf0p72we1pH9KedsHtmkQfktrsRUVY09KwpKdjPZSONT0dS3o6tsxMcNYMX68PD8MUG4cxLg5jXCymVq0wHn8dFdlke7+rxNFIicNZWcmet1+h+sNPMVklFeOH0PuRl9D51FR8v7B8D2//eIi3b+nJsE61mtluehuW/wNGvAR9pjZKvErz4HRKtmUUs2pPLhvS8tmVVYKU4GvSMT62iKuN2+hQ9SsBBdsRTjvozdpIBK36aI0vYnuBn2enWLY5bRwuOczewr3sL9rPvsJ97CvaR2F14YkyIV4hJAQkkBCYQEJAAvEB8SQEJtDKv5VHx9c6ldNqxXbkCBZXMrFm/I4tMwtbZia27OyTkgp6PcboaFcSicIQGYkxKhJDpPYwRkaiDwnxyG0wlTjcnDjsRUUULlhA9vtzMZVUktrei47/eoGOlww/qdyirZn8/fMd3NSnNc9f26VmQ9oqmD8Bkq+EiZ822jATSvNUUmnj50MFbEjLZ1N6AQeOlWuJhCrGBB/mKu99dLbuJLTiADrpmh89KB5iL4bIThDRCSI6aOs8/H8xvyqfvYV7OVB0gMOlhzlccpjDpYdPSih6oSfKN4oYvxiifaOJ9YutefaLJso3qskkFmmzYcvJwZaZiTUzsyahZGZiy83FnpcHDsfJOxmNGMPDTyQTQ2gohrBQ9CGhGEJDtGfXa52vT4M1PlCJww2JQzqdVG3dSvFXSyn+einCYmVHgiD9uku48+aZhHiFnFR+0dZM/rFoB/3ahPH+lEsw6l1fyAOrYMHNENoWblsOZv8Gj1Vp2cqqbezMLGF7RjHbfi9md1YJOaXVeGGhs0inr/kQ/c2Hae88SIgt+8R+0uiLiLhISyIhSdrt0+BEbfh3b8/21SixlHCk9MiJZJJVnsXR8qMcrThKXmUekpq/bTqhI8w7jHDvcMK9wwnzCTvx+sR6n3BCvUI93pFROhzY8wuwH8vFlpODPfcY9txcbLk1y/bCQpylpXXuL7y8MISEoA8NxRASQuTjj2OKiz2nWFTiaKDE4ayooGLTZirWr6ds7Vrs2dlYTTp+6iDZclkkt415ikGtBp20T7XNwYvf7eO9DekMaBvG27f0xNdsAIcdNrwKa/+t/cq7ebHHbxcoLUdxpZX9ueXsyyllX24Z+3LKOHCsHFtlKe1EJu11GbQXGXQ2ZJEsMgmSxSftbzMF4giMRx8UiyEwFhEQBf4x2igHAa5nryCP9S/JqcghqyKL7PJsssqzyKnIIb8qn7yqPPKr8imqLjopuRwXYAogyBxEoDmQQHPgH5dNrmevQAJMAfgaffEz+jV6yzBptWIvKsKen4+jsBB7QQGOguPPBdgLC3EUFBA3+02MUec28oRHEocQYjjwGqAH3pFS/ueU7WZgHtATKABukFIedm17BLgdcAD3SylXnM0x63KuiUM6HFjT06n+7TeqU3+jKnU3Vdt3gN2O3aRnd7zgxw5OcnvGc1OP2xjdZjRmvfnE/lVWB19tz2LWmjSyiqu4tV8Cj1x9EWYcsPdr+PFFyNsLncZpw2h7sOmkohxXUmUjo7CSIwWVHCms4PeCSjKKKiktKcJUmkGY7SitRS7xIpfW4hgRopgoUUiQqPjDsRzCgMUYhM0chMMcgtM7GLxD0PmFYvALxegXisk3EL1XAJj8wewHJj/tqtvsDzp9HRE2DJvTRmFVIfnV+eRX1iSUwupCii3FlFhKKLGUnFgut5X/6fGMOiN+Rj98jD74Gf20hGLyw9fgi6/J98Q2L70XXgavk57NBjPeBm/MevMftnkZvDzWgbLRE4cQQg/sB64EMoEtwI1Syt9qlbkb6CqlnCaEmAhcK6W8QQjREfgU6A3EAKuAdq7d/vSYdTnXxJF60/Xotu4GwGbUkRVlYEesje2JgqNJAQxKuoLhCcPpFdGbMouTogormcVV/Ha0lO0ZxWxMO4bJWkKfKMl9vXzpJA7D0W1wcA1Ul0BIG7jiKa25ZBPpEKUoZ1JhsZNbWk1uqYXc0mpySqspqrBSXl6KLMtFX56NuSoXX0se3vZiAmUZwaKcYFFGMNpyEOUYheOM72URZqp1Plh13th1Zhw6M069GYfejNSbcRq8cOrMSIPXidfovZAGMzq9EfRGdAYjQq89dAYTOr1BezYYT5QR+uPLBoReKyN0enQ6PUIn0OkMOJCUOSopt5dT5qikxFZOub2CKkc1lY5KKh3VlNsrqbBVUm6vpNJWSbmtggp7BRU27VFlrzqnf3OTzoRZb8aoN2LUaQ+T3vTHZf3Jr016E/d0v4cIn3PrFHq6xOHONmC9gTQp5SFXAJ8BY4Haf+THAk+5lhcBbwitVmcs8JmU0gKkCyHSXMfjLI7ZYL6IO0xOtI7sCNAFQLTDRjuL5M5qO20PHsN4cD6CjymWTnRIgpCEIOmOZJKQ+Oqq0Hk5oRgt9QEExGlt6juN0/pqqEpw5QLjazaQFO5HUrjfGcs6nZKyajvFVVbKLXaKLQ6yLHbKq21YKoqxVxRiryrFXlmKtJSBtRy9rQKDvQKDvRyToxKjvQKTswqDzYpRWjA6LRhlJWasmLHhI6wnlr2wYjqLhHSuQs5cpE4OKXCiw46gSuioFoIqnQ6LTlAtBNWuddWu1xZRa70OLKIKqxDYhMAqwCYENlzPQmATUOF6rlmvPY/0v4yILkMb8F/BvYkjFsio9ToT6HO6MlJKuxCiBAh1rf/llH2P1+6c6ZgACCGmAlMBWrdufU4n0K1jPwYUHCBAGhGVegQCKXTgo+OYrwChQwgdRoMek1GP2WjA22Qk2NeMyWAAr0DwCdUefhEQ2Rl8Q88pFkW5EOl0gkAfY01fpZOcW4XtcXaHE6vDic0hsTucVDklZU6J3WbDbq3GYbfisNuw26xIuxW73YbTbsVpt+OwW3E6bEiHDRw2cNq0lkxOG8JpRzhtOKVT61vldIB0IJ1Se5ZOqPWQzuPLDoR0attdZQVOhJQI6QApXa+dcLxuRUpAIpB4S4kXIJCuzdK13bXsKodTnlglXOuRWm3NiddIbWIwIDG07Xn9O9elafY6aQBSyjnAHNBuVZ3LMcZOfq1BY1IUpeEY9DoM+rqu2L0BVV/oTu68T5IFtKr1Os61rs4yQggDEIhWSX66fc/mmIqiKIobuTNxbAGShRCJQggTMBFYekqZpcBk1/J1wBqp1dYvBSYKIcxCiEQgGdh8lsdUFEVR3Mhtt6pcdRb3AivQms6+J6VMFUI8A6RIKZcC7wIfuSq/C9ESAa5yC9Eqve3APVJq3VvrOqa7zkFRFEX5I9UBUFEURanT6ZrjqragiqIoSr2oxKEoiqLUi0ociqIoSr2oxKEoiqLUS4uoHBdC5AFHznH3MCC/AcO5EKhzbhla4jlDyzzvcz3neCnlH4bsbhGJ43wIIVLqalXQnKlzbhla4jlDyzzvhj5ndatKURRFqReVOBRFUZR6UYnjzOZ4OgAPUOfcMrTEc4aWed4Nes6qjkNRFEWpF3XFoSiKotSLShyKoihKvajEcRaEEE8JIbKEENtdj6s9HZO7CCGGCyH2CSHShBAPezqexiCEOCyE2OX6bJvlaJhCiPeEEMeEELtrrQsRQqwUQhxwPQd7MsaGdppzbtbfZSFEKyHEWiHEb0KIVCHEdNf6Bv2sVeI4ezOllN1dj2WeDsYdhBB64E1gBNARuFEI0dGzUTWay12fbXNt3/8BMPyUdQ8Dq6WUycBq1+vm5AP+eM7QvL/LduD/pJQdgb7APa7vcIN+1ipxKLX1BtKklIeklFbgM2Csh2NSGoCUch3anDe1jQU+dC1/CFzTqEG52WnOuVmTUmZLKX91LZcBe9Amd2/Qz1oljrN3rxBip+vyt1ld0tcSC2TUep3pWtfcSeB7IcRWIcRUTwfTiCKllNmu5Rwg0pPBNKKW8F1GCJEA9AA20cCftUocLkKIVUKI3XU8xgL/A9oA3YFs4BWPBqs0tAFSyovRbtHdI4S4zNMBNTbXlM0toW1+i/guCyH8gC+Av0kpS2tva4jP2m1Tx15opJRXnE05IcRc4Bs3h+MpWUCrWq/jXOuaNSllluv5mBBiCdotu3WejapR5AohoqWU2UKIaOCYpwNyNyll7vHl5vpdFkIY0ZLGfCnlYtfqBv2s1RXHWXD9Qx93LbD7dGUvcFuAZCFEohDChDYH/FIPx+RWQghfIYT/8WXgKprv53uqpcBk1/Jk4CsPxtIomvt3WQghgHeBPVLKGbU2NehnrXqOnwUhxEdol7YSOAzcWet+YbPiap74KqAH3pNSPu/hkNxKCJEELHG9NACfNMdzFkJ8CgxGG147F3gS+BJYCLRGm3ZggpSy2VQmn+acB9OMv8tCiAHAT8AuwOla/ShaPUeDfdYqcSiKoij1om5VKYqiKPWiEoeiKIpSLypxKIqiKPWiEoeiKIpSLypxKIqiKPWiEoeiKIpSLypxKIqiKPWiEoeieIAQ4hLXQHtert7rqUKIzp6OS1HOhuoAqCgeIoR4DvACvIFMKeULHg5JUc6KShyK4iGu8cC2ANVAPymlw8MhKcpZUbeqFMVzQgE/wB/tykNRLgjqikNRPEQIsRRtlsVEIFpKea+HQ1KUs6Lm41AUDxBCTAJsUspPXHO9bxRCDJFSrvF0bIpyJuqKQ1EURakXVcehKIqi1ItKHIqiKEq9qMShKIqi1ItKHIqiKEq9qMShKIqi1ItKHIqiKEq9qMShKIqi1Mv/Ay+Dsq4Do4ttAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# 確率密度関数\n",
"x = np.arange(-5,20,0.1)\n",
"pdf1 = stats.gumbel_r.pdf(x, loc=0.5, scale=2.0)\n",
"pdf2 = stats.gumbel_r.pdf(x, loc=1.0, scale=2.0)\n",
"pdf3 = stats.gumbel_r.pdf(x, loc=1.5, scale=3.0)\n",
"pdf4 = stats.gumbel_r.pdf(x, loc=3.0, scale=4.0)\n",
" \n",
"plt.plot(x, pdf1, label=(r'$\\mu=0.5, \\eta=2.0$'))\n",
"plt.plot(x, pdf2, label=(r'$\\mu=1.0, \\eta=2.0$'))\n",
"plt.plot(x, pdf3, label=(r'$\\mu=1.5, \\eta=3.0$'))\n",
"plt.plot(x, pdf4, label=(r'$\\mu=3.0, \\eta=4.0$'))\n",
"plt.xlabel(\"x\")\n",
"plt.ylabel(\"f(x)\")\n",
"plt.title(\"Probability Density Function\")\n",
"plt.legend()\n",
"# plt.savefig(\"../img/gumbel_pdf.png\")"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3hUxfrA8e9sSe+FhCRA6IQWSsAGwpUiIvaGiv6uotix93q5ioCooCKKggULIoKgooCKClykQ0InQCC992y2zu+PXSQJgbTd1Pk8zz5bzjkz72Y3592Zc86MkFKiKIqitF2apg5AURRFaVoqESiKorRxKhEoiqK0cSoRKIqitHEqESiKorRxKhEoiqK0cSoRKE4lhHhFCPFFA7bfJ4QY6cSQmqRuIcStQoi1FZ5LIUQ3Z5TtKK9ECNHFWeU1J0KID4QQLzZ1HG2JSgSthBDiFiHEdscOIl0I8bMQYlhTx3UuQohPhRCvVnxNStlHSvmHk+uJduyISxy3TCHEj0KIMXWtu0JZunOtJ6X8Uko51gnhI4T4QwhxV5XyfaSUx5xRfpW6koQQhgp/qxIhRISz66lQ37+FEBsrvialvFdK+V9X1amcSSWCVkAI8RgwB5gOhAEdgfeBq5oyrmYoQErpA8QC64AVQoh/O7uSmpJEC3CFI9GcuqU1dUCKi0kp1a0F3wB/oAS44RzrfAq8WuH5SCClwvMk4EkgHigFFmJPKD8DxcCvQGB121bYfrTj8SvAFxWWfQtkAIXAX0Afx+tTADNgcsT/Q8WygAjAAARVKGsgkAPoHc/vBA4A+cAaoNNZ3n80IAFdldefADIBTTXvYyiwHShyrPOW4/WTjrJKHLcLgH8Dm4C3gVzgVcdrGyvUJYGpwDHHe3ijQr1V/2b/xAu8BliBckd971Uor1uF78DnQDZwAnihQtn/BjYCsx1/p+PAZef4rvzzN6jyem0+96WOOIqBfUBchXU7AMsdMeYC7wExjvdldby3grN8X+8GEoE8YBUQUeXvei9wBCgA5gGiqf8vW9pNtQhavgsAD2BFA8u5DhgD9ACuwJ4EngNCsbccp9az3J+B7kA7YCfwJYCUcoHj8Sxp/9V5RcWNpP1X6GZHXKfcAiyTUpqFEFc54rvWEeMG4Os6xrbcEVfPapbNBeZKKf2Arth3cgAXO+4DHHFvdjw/D/tOPgz7zrs61wBxwCDsrbU7awpQSvk89vf2oKO+B6tZ7V3syaALMAK4HbijwvLzgENACDALWCiEEDXVXQ9XAkuAAOw77PcAhBBa4EfsSSoaiASWSCkPYN+Jb3a8t4CqBQohLgFeB24E2jvKWFJltQnAEKC/Y71Lnf3GWjuVCFq+YCBHSmlpYDnvSikzpZSp2Hc8W6SUu6SU5diTzMD6FCqlXCSlLJZSGrH/aowVQvjXcvOvgJsBHDuuiY7XwL4DeV1KecDx3qcDA4QQneoQ3qkuj6BqlpmBbkKIEClliZTy75rKklK+K6W0SCkNZ1lnppQyT0p5EntX3s11iLVajp3sROBZx985CXgTuK3CaieklB9JKa3AZ9h3qGHnKPZ7IUSB4/Z9HcLZKKVc7ahnMfYuOLC3riKAJ6WUpVLKcinlxrOWUtmtwCIp5U7Hd+hZ4AIhRHSFdWZIKQscf9f1wIA6xKygEkFrkAuEOKFfOrPCY0M1z33qWqAQQiuEmCGEOCqEKMLelQD2X6a18R32f/r22H+J27AnKYBOwNxTOyzs3QYC+6/N2jq1bl41yyZjbx0dFEJsE0JMqKGs5FrUV3GdE9h3jg0VAugd5VUsu+LfIePUAyllmePhuT7Pq6WUAY7b1XWIJaPC4zLAw/G97IA9GdXnx0oEFd6blLIE+3e+2vfnqLfO39W2TiWClm8zYATO9Q9bCnhVeB7egPoqleX4RRp6lnVvwd4FMhp710X0qc0c9+cc+lZKmQ+sBW5ylLVESnlqm2Tgngo7rAAppaeU8n91eC/XAFnYu02q1n1ESnkz9q6jmcAyIYT3OWKuzTC+HSo87sjpFklNn8+5ys7B3nqp2BLqCKTWIp66qMvnXlUy0PEsP1Zq+rulUeG9OT6DYJz//to0lQhaOCllIfASME8IcbUQwksIoRdCXCaEmOVYbTcwXggRJIQIBx5pQJWHsf/Su1wIocd+YNL9LOv6Yk9Sudh3ItOrLM/E3q99Ll9h7/O+ntPdQgAfAM8KIfoACCH8hRA31OYNCCHChBAPAi9j71KxVbPOJCFEqGNZgeNlG/aDnbZaxF2dJ4UQgUKIDsDDwDeO13cDFwshOjq6zZ6tst1Z/06ObpilwGtCCF9H19hjQL2v5TiLunzuVW0F0oEZQghvIYSHEOIix7JMIEoI4XaWbb8G7hBCDBBCuGP/Dm1xdIEpTqISQSsgpXwT+z//C9h3VMnAg8Cp/t3FwB7sXTNrOb0Dqk9dhcD9wMfYf5WVAilnWf1z7M36VGA/ULWffSHQu4a+6FXYDzZnSCn3VIhjBfZf6ksc3U57gctqCL9ACFEKJADjsZ9ptegs644D9gkhSrAfOJ4opTQ4ulZeAzY54j6/hjorWgnswL7j/wn7+0dKuQ77ZxLvWP5jle3mAtcLIfKFEO9UU+5D2D+HY9jPEPoKONv7qpc6fu5Vt7ViPwGhG/azrlKwt/IAfsd+hlGGECKnmm1/BV7E3k2Yjv3A/cSGvBflTOJ0S1tRFEVpi1SLQFEUpY1TiUBRFKWNU4lAURSljVOJQFEUpY1rcYNjhYSEyOjo6KYOQ1EUpUXZsWNHjpSy2ms/WlwiiI6OZvv27U0dhqIoSosihDhxtmWqa0hRFKWNU4lAURSljVOJQFEUpY1TiUBRFKWNU4lAURSljXNZIhBCLBJCZAkh9p5luRBCvCOESBRCxAshBrkqFkVRFOXsXNki+BT7CI5ncxn2USW7Y5+/dr4LY1EURVHOwmXXEUgp/6oynVxVVwGfOyYa+VsIESCEaC+lTHdVTIqiOJfVJjFbbVhtEotNYj11s1qxmMuRFitWqwmr1YrNasZqsWCznrqZsVmtSMdjabM/xmYBaQWbFWwWhM2KlDYkEmw2pLQhkPbXpARpf4yUcOresRybBE4tq7AeFdc9vVwg/1kuTs2ZU2GEZllxHh3H66LS3DqywlQ7px6f3soe96nH9iWntpcV1rUgkdjfs1XasDmeRQy4kv5xY53z4VXQlBeURVJ56r4Ux2tnJAIhxBTsrQY6duzYKMEpSkslpaTUZKXQYKbIYKbUaKHMZKXMaMFUXoylNB9bWSHW8mJspmKksQxpKkWYy9BaytBaDOisBtxsBvQ2A3pbOTqbCZ00o5On7i3opRk9ZtyEBTcsuGHGw3GvE2fM9dNs2IByISjTCMqFwCQERsfNVOX+3K+BWQgsQmAB+32Fx2YBVgQWAZZT90Jgwb6s0nYIrBqwCnHO2CcdLm91iaDWpJQLgAUAcXFxagIFpU0pNVrIKTGSU2Iku9hETomR3BIThUWFWIoy0ZRlozPk4G7Mxcuch68lj0CK8KMMf1FKEKVEizL8KMVNWGtVpwUdRuGBUeOBSeOBReeOVbhh07ph03ghNXpsWjfMGjdMWjekVo9N447UuiG1bqB1Q2r0oNUjNFrQ6NBodKDVIjQ6hPb0c41jHaHRodFqQehAqwON/bHQagFBOWYMtnKKrQbKbAZKrAZKreWUWssotRoot5kotxkxWI2U28oxWI0YbOWUW8sxOG7lViMGa3mDPg93jRtuWnf0Gj16rR6d0KHTOG5Ch06j/eexh+b0Mr1Gf3q9f9bVodfo0Gn1aIV9O43QotNoHfc6NEKDVmjRCi2DwwY3KPazacpEkErlOVyjUPOQKm2MlJLMIiMp+WWkFhjst3wDqflllOWnoylMJtiSRaTIJkLkEilyGCxyiBQ5+Imyasssd/PB6BaIxS0Am3s40iMAm4c/hZ4BaL0C0HkHovMOxN3LH62HD+i9wM379L2bNzqtHh3g7aL3bbKayCvPI688j2zHfZ7Bfp9bnkt+ST755fkUm4spMhZRbCrGIi3nLNNd646XzgsvvReeOk+89F746gIJc7xWaZnOfu+h88Bd646b1q3S/dle02v0iBp+tdeVlBJpMGAtLMRaVIS1sBCb495aWIS1KN/+vKiYgOs6wgXdnFo/NG0iWAU8KIRYApwHFKrjA0prZbHaOJZTytGsEhKzSjiaXcLR7FKOZRfRzpxGN5FKV5FGV5HGMF06XUQaPrLUfjqHYzZfs94Xi28Uwr8n+uCx4B8B3u3Apx14h9rvvULw0Hvg0UTvU0pJobGQzLLM07fSyvdZZVmUmEuq3d5N40aQZxBBHkEEugfSwbcDvm6++Ln72e/d/PBz8/vnNT+9H37ufnjrvdFpmkcHh7RYsObnY8nLw5qbiyU3D2ue4z4/z/E875+dvrWoCMzmsxeo1aL19UXj54fPiBEuidllfzkhxNfASCBECJGCfaJwPYCU8gNgNfZ5YxOBMuAOV8WiKI3JYrVxMKOYhNRC9qYWsjetiIPpRZgsFrqLVPqJ41zomcx92iQ6647hpjH8s63NJwxNaE8IGQHB3SGwE/h3gIAO6D387f9ATcxkNZFakkpKcQopJSkkFyeTUmy/Ty1JxWAxVFpfIzSEeIYQ7hVO14CunN/+fII9gwnyCPrnFuwRTKBHIN56b6f/4nYWaTZjycnBkpmJOSsLS2YWlqxMzJmZ9sc5OVhzc7EWFFRfgFaLLigIbXAw2sAA3MPD0fr5ofX3Q+vvj8bPD62fP9oAf7R+fmgcjzXerv+btLg5i+Pi4qQafVRpTgwmK7tO5rM1KY/tSfnsPJlPmcmKG2bO8zjBeN8k4jQH6VS2FzdzkX0jvTe07w/tYyG8H4TGQEg38PBv2jdTQZm5jONFxzlWcIyjBUc5WniUYwXHSClJwSZPHwz20HoQ5Rtlv/lEEeETQZhXGGHeYYR5hRHiGdJsfq2fi81gwJyaijk1FVNKCuaU1H+em7MysebkVjqDCAC9Hn1oKLqwMHShoehCgtEGBaMLDrLfBwWiDQ5GFxSExs8PoWm6a3iFEDuklHHVLWv+n46iNDNSSo5ml7D+YDZ/HM5i6/E8zFaJEDAitIxZHQ8wxLyd0NwtaCzlUAyE9IB+10DHCyByMAR3tR8MbQaklGSWZXIw7yAHcg+wP28/h/MOk1aa9s86Oo2OaL9oegb15LLOl9HJrxNRvlF08O1AsEdws/0VX5W1sBDT8eMYjydhSkrCnJyMKTUFc2oa1pycSusKNzf0kZHoIyJw79UTfVgYunZh6MLa2R+HhaENDGzSnbuzqESgKLVgs0m2JeXxU0I6vx/MIiXf3v3Ro503z8UaGGXbRGT2X2hzj0AREBgNg26HziOg4/ngHdKk8VdUZCoiPjue3Vm72ZuzlwN5B8grzwNAIOjs35nY0Fiu7X4tXQO60iWgCx18O6DXNIeOqZpJmw1zairGI4mYjh/HlHQc4/HjmI4nYc3NPb2iVou+fXv0UVH4jByBW2Qk+qgo9JFR6CMj0YWGtIqdfG2oRKAoZyGlJCG1kBW7UlmdkE5mkREPvYZh3UJ5Ok4wsvx3fBNXwf4k0OghehgMmQzdxth/8TeTX8mpJansyNzB7qzd7MraxdGCo0gkGqGhW0A3Lo66mJigGHoH96ZHYA+89F5NHXKtWUtKMB4+jPHQIcoPHcJ46DDGw4exlZb+s442KAi3zp3x+ddI3Dt3xq1zZ9yiO+PWIQqhbxnJzdVUIlCUKgoNZlbuTuXrrckcSC/CTadhZI9QruoTwGjb/3CPnwMb/gahhS4jYPgTEDMBPAObOnQACo2FbM3Yyt9pf7M5fTPJxfbrNn30PsSGxnJp9KUMaDeA/iH9W9ZOv7AQw969lCckYEjYi/HgQcypp8841/j54dGzJ/7XXIN7zx54dO+OW+fOaP2bz3GX5kolAkVxOJpdwscbjrN8ZwpGi40+EX789+q+XN1F4rvrI1j7ORiL7GfzjPkvxN4MPtVOAduopJTsz93P+uT1bErdxP68/dikDS+dF0PDh3JrzK0MCR9CV/+uaJvJcYma2EwmjAcOYIhPwJAQT3l8AqakpH+Wu3XujGdsfwJuuAH3Xj3x6NkTXXh4izlW0dyoRKC0eduS8vjwz2P8eiATN52GawdGMun8TvTVnoRNr8La5fazRfpcDUPush/wbeIdjtlqZmvGVtYnr2d98nqyyrLQCA39Q/pzT/97uCDiAvqG9G0x/fq2sjIMe/ZQtm07Zdu3Y9izB2k0AqALDcUjtj/+11yDZ/9+ePTpg9bPr4kjbl1UIlDarD3JBcxee4gNR3II9NIzdVR3bju/E6HlSbD+Idi/Etx8YOg9cP69ENC041xZbVa2Z25n9fHVrDuxjmJTMZ46Ty6MuJBLOl7CxZEXE+AR0KQx1pbNYKBs+3bKtm6lbOs2DPv2gcUCGg0evXoROPEmPAcNxnNALPqwsKYOt9VTiUBpc45kFvPGmkOs3Z9JoJee58fHMOn8TniWpcKvj0D8EvtwCxc/BRc8AJ5Nu3M9kHuAVUdXsSZpDdmGbLx0XozqOIqx0WM5v/35eOia6jri2pNSYjx0iNJNmyjZuBHD9h1Isxn0ejz79iX4jjvwGhKH58CBaH19mzrcNkclAqXNKDFamPvrYT7ZlISnXstjY3pw57DO+Ggs8L+3YMOb9qGJz78fhj3apKd8lphKWH18NcsOL+NA3gH0Gj3DI4czvst4Lo66GE+dZ5PFVlu2sjJKNm6k5Pf1lGzaiDXbfp6+e/fuBN56K97DhuE1eBAaz+b/Xlo7lQiUVk9Kyao9abz20wGyS4zcFNeBJy/tSbCPOxz5FX5+EvKOQcyVcOl0COhQc6Eusj93P0sOLuGXpF8wWAz0COzBs0Of5fIul+Pv3vzPfrHk5FC8fj0lv/1O6f/+hzSZ0Pj743PRRXhfdBHewy5SXT3NkEoESquWVVTO09/Fs/5QNv2j/FlwexwDOgSAoQCWT7V3AwV3h0nLoduoJonRJm38mfwnn+//nO2Z2/HUeTK+83iu634dfUP6NvszYSzZ2RT9/AtFP/+MYfdukBJ9ZCQBE2/C95JReA0epM7Xb+ZUIlBarR/j03jh+70YTFZevqI3t18QjVYj4Oh6WPkAFGfAiKft1wHo3Bo9vnJLOSsTV7L4wGJOFJ0g3DucJ+Ke4Nru1+Lr1rz7ya2FhRSvW0fhTz9RtmUr2Gy49+xJyIMP4DtqFO49ezb7BKacphKB0uqUGC08tzyBVXvSiO0QwFs3xtI11AcsJljzImz5wD72z13r7OP+NDKj1ciyw8tYmLCQbEM2fYL7MOviWYzuNLpZn+4pzWaK//iDwhXfU7JhA5jN6Dt1JOTee/AbPx73bs4fJ19pHCoRKK1KYlYJ936xg2PZJTw2pgf3j+yKTquBojRYejukbIPz7oXRr4C+cQ9SGq1Gvjv8HQsTFpJlyCIuLI6ZF88kLiyuWf96Nh47RsGy7yhcuRJrbi66du0ImjQJv8svx6NP72Ydu1I7KhEorcYve9N54tt43HUavrjrPC7s6jjrJ2kjfPtvMBvghs/sF4Y1Ipu0sTJxJe/tfo+ssiwGtRvE68NfZ2j7oY0aR13YDAaKfv6FgmXLMOzcCTodPiNHEHD99fgMG4bQqV1Ha6I+TaXFk1Ly9rrDvPN7IrEdAph/6yAiAhy/9rd+BD8/DUFd4N8/QWjPRo1te8Z2Zm2bxYG8A/QP6c9rw17jvPDzmu2vaHNGBvlffkXB0qVYCwtxi46m3ZNP4H/VVehCms8IqopzqUSgtGhmq41nlyewbEcKN8ZF8d+r++Ku09qHhPhtGmx8C3pcBtcuAI/GG5YguSiZt3a8xa8nfyXcO5yZw2dyWefLmmUCkFJi2LWbvMWfU7x2HUiJ76hLCJx0G15DhzTLmBXnUolAabHKTBbu/3InfxzK5pHR3Xl4VHf7TstqhlVTYc9XMPjfMP5N0DbOV91oNbIwYSEfJ3yMTqPjoYEPcXvv25vl1b/SaqV43TpyFy6iPCEBja8vQbffTuCtt+IWFdnU4SmNSCUCpUXKKzVxxydbSUgt5PVr+3HzUMc4QKYyWHobJP4KI5+DEU812gBx2zK2MW3zNJKKkhjfeTxPxD1BqFfTj05alTSbKVz1A7kffYQpKQm36GjCXnqRgKuuQuPt3dThKU1AJQKlxSkoMzHp4y0czS7hw9viGNPbcaWq2QBf32Q/OHzFXHtroBGUmEp4Y/sbLD+ynCifKD4c/SEXRl7YKHXXhc1opODbZeQuWoglLR333jFEzpmD75jRCG3LGJ5acQ2VCJQWpdBg5raFW0nMLuHj2+O4uIfjF7e5HJbcAsc3wDUfQOzERolna/pWXtj0ApllmdzZ907ujb232Y0DJM1mCpavIOf997FkZuI5eDDt//MfvIcNU/3/CqASgdKCFJWbuX3RVg5lFPPhbYNPJwGL0X6NwNHf4cr3GiUJlFvKmbtzLl8c+IJov2gWX7aY/qH9XV5vXUirlaLVq8l+9z3MJ0/iOWAAETNn4n3+eU0dmtLMqESgtAgGk5U7PtnGvtRCPpg0mH/1amdfYLXAsjvhyBqYMAcG3ebyWI4WHOWJP58gsSCRW3rdwiODH2l2rYCSv/4i643ZGI8cwT0mhqgP5uMzYoRqASjVUolAafasNslDX+9i18l85t0yiNGnjglICT8/BQd/hMtmQdwdLo/l+8Tvmb5lOp46z2Z5LMB45AiZM2dRunEj+k4diZzzNr5jxyI0mqYOTWnGVCJQmjUpJS+v2suvBzKZdlUfLuvX/vTCjW/D9oVw0SNw3j0ujcNgMfDq36+y6ugqhoYPZcbwGc3qjCBLXh7Z77xDwdJv0fj40O6Zpwm65RaEW+MPpqe0PCoRKM3axxuO88XfJ7lnRBduvyD69IK938Fv/4F+N8Col10aQ2pJKo+sf4RDeYe4P/Z+pvSf0mwmgZdWK/lLlpA9Zy62sjICb76ZkAcfQBcY2NShKS2ISgRKs/Xn4Wxe//kA4/uF8/SlvU4vSN8D3z9gn0T+qnngwm6PbRnbePyPx7HYLMwbNY/hUcNdVlddGfbuI+OVVyjfuxevC84n/IUXcO/atanDUloglQiUZul4TikPfbWTHmG+zL4hFo3GcZCzJBuW3ApewXDj56Bzd1kMXx/8mplbZ9LRryPv/Osdov2jXVZXXViLi8me+w75X32FNiiIiNmz8bt8vDoQrNSbSgRKs1NmsjDl8+1oNYKPbo/Dy83xNbVa7KOIlmbDnWvAp51L6rdJG7O3z2bx/sWMiBrBjOEz8HHzcUlddVX8xx9kvPQyluxsAm++mdBHHkbr13hjKCmtk0oESrPz4vf7SMwu4YvJ59EhyOv0gj9nwImNcPUHEDHAJXUbrUae2/Aca0+s5ZZet/DUkKeaxfEAa2EhmdNfp3DlSty7dydq3nt49uvX1GEprYRKBEqz8u32ZL7bmcLDo7pzUbcKwx4n/gZ/zYaBk2DAzS6pu9BYyNTfp7IzaydPxD3B7b1vbxbdLf+0AnJzCb7vXkLuuw+NOhtIcSKXJgIhxDhgLqAFPpZSzqiyvCPwGRDgWOcZKeVqV8akNF9HMot5ceVeLugSzNRR3U8vKM6A5VMgtBdc9oZL6s4uy2bKuimcLDrJ7BGzuTT6UpfUUxfWklIyX3uNwhUr7K2A99/Hs2+fpg5LaYVclgiEEFpgHjAGSAG2CSFWSSn3V1jtBWCplHK+EKI3sBqIdlVMSvNlsth45JvdeLnpmDtxgH2SebBfNLbyQTCV2CeWcfM6d0H1kFaSxt1r7ybbkM380fObxcxhhoQEUh9/AnNKCsH33kPI/ferVoDiMq5sEQwFEqWUxwCEEEuAq4CKiUACp450+QNpLoxHacbe+e0I+9KK+PC2wbTzqzB2/45PIHGd/crhdr3OXkA9JRUmcfe6uyk1l/LR2I+IDY11eh11IW02chcuJHvuO+hCQ+n0+Wd4xcU1aUxK6+fKRBAJJFd4ngJUHe3qFWCtEOIhwBsYXV1BQogpwBSAjh07Oj1QpWntOJHP+38kcsPgKC7tE356Qe5RWPM8dPkXDLnb6fUeLTjK5DWTkUgWXbqIXkHOTzR1Yc7MIu2Zpynb/De+l15K+2n/Qevv36QxKW1DUw9AcjPwqZQyChgPLBZCnBGTlHKBlDJOShkXGtp8LutXGq7cbOXJb/fQ3t+Tl67ofXqBzQbf3w9avUsuGksqTOKutXchhOCTcZ80eRIo/XsLx6+9FsPuPbR/9b9EznlbJQGl0biyRZAKdKjwPMrxWkWTgXEAUsrNQggPIATIcmFcSjPyzm9HOJZTypd3nYevh/70gh2LIPlvuHo++Dt32sTkomQmr52MTdpYdOkiuvh3cWr5dSGlJG/hQrLeehu3zp3p9Pln6upgpdG5skWwDeguhOgshHADJgKrqqxzEhgFIISIATyAbBfGpDQj+9OK+PCvY9wwOKryqaKFqbDuFXuXUKxzTxVNK0lj8trJmKwmPhr7EV0Dmm6nay0pIXXqw2TNfhPfsWOJ/uYblQSUJuGyFoGU0iKEeBBYg/3U0EVSyn1CiGnAdinlKuBx4CMhxKPYDxz/W0opXRWT0nxYrDae/i6eQC89z18ec3qBlPDT4yCtcMUcp843nGPI4a61d1FiLmHh2IX0COzhtLLrypiYSMqDD2FKTraPFPp//9csrllQ2iaXXkfguCZgdZXXXqrweD9wkStjUJqnL/4+QUJqIe/ePJAArwqnRR78CQ7/DGNfhcBop9VXYirh/l/vJ8eQw8djPyYmOKbmjVyk5M8/SX3scYSnJ50+/QSvIUOaLBZFgaY/WKy0QTklRt5cd5jh3UOY0L/C/AKmMvjlWWjXG867z2n1mawmHln/CEfyj/DWyLeabEpJKSV5n31G8n33o+/Ukc7LvlVJQGkW1BATSqOb+fNBys1WXr6iT+XukI1vQ+FJ+4VjWud8NW3SxnMbn2NLxhamD5vOsMhhTim3rqTZTMZ/X6Vg6VJ8x4wmYuZMNF7OvzhOUepDJQKlUe08mc+3O1K4Z0QXurWrMKJn3kPJRgYAACAASURBVDHYNBf6Xg/RzttZz94+mzVJa3gi7gmu6HqF08qtC2tBASkPP0LZli0ET5lC6CMPq6kjlWZFJQKl0Ugp+c8P+wnzc2fqJd0rL1z7Imh0MPa/Tqtv6aGlLN6/mEkxk/i/Pv/ntHLrwpyWxsm7p2A+eZL2M14n4OqrmyQORTkXlQiURrM6IYM9yQW8cX1/vN0rfPVObrFPQP+v58Evwil1bU7bzPQt0xkeOZwn4p5wSpl1VX74MMl3T8FWWkqHhR/jPbTpxzBSlOqo9qnSKEwWG7PWHKRXuC/XDoo6vUBKWPci+ITBBQ84pa7jhcd5/M/H6ezfmVkXz2qS+QTKtm/nxKTbwGaj05dfqCSgNGsqESiN4qstJziRW8bTl/U6PbIo2E8XTd4CI58FN+8G11NoLOTB3x5Er9Hz3qj3mmRmseJff+XknZPRBQcTveRrPHr2bPQYFKUuVCJQXK643Mw7vydyQZdgRvaoMFaU1QK/vgIhPWDgbQ2ux2qz8vSGp0krTWPuv+YS6ePcoSlqo2DZMlKmPox7TC86ffUl+sjGj0FR6kolAsXlPvzzGHmlJp4d36vy6aK7PofcIzD6FaecLjp/z3w2pW7i2aHPMqCda6ayPJe8L78k/YUX8b7wQjp98gm6wMBGj0FR6kMdLFZcKrOonI83HuPK2Aj6RwWcXmAsgfWvQ8cLoOf4BtfzR/IffBj/IVd3u5obetzQ4PLqKnfhIrLeeAOfUaOIfPstNYmM0qKoRKC41PvrE7FYJU+MrdJPvu1jKM2Cm75o8HhCJ4pO8OyGZ+kd3Jvnz3u+UcfskVKSM38+Oe+8i+9l44icNQuh19e8oaI0IyoRKC6TWVTO19uSuW5QFB2DK1xFayqF/70LXS+BjlXnKqqbcks5j/7xKDqNjrdHvo2HzqPmjZxESkn223PIXbAA/6uuov301xDaxj9DSVEaSiUCxWU+/PMYVpvkgX91q7xg+ydQlgMjnmlwHbO2zeJI/hHmj55PhI9zrkGoDSkl2W+9Re5HHxNw442Ev/KyulpYabFUIlBcIqu4nC+3nOCagZFVWgNl9qEkOo9ocGtgTdIavj38LXf0vaPRxxDKeW+ePQlMvInwl19WQ0grLZpKBIpLfPTXMcxW25mtgZ2f2Y8NjPi0QeWnFKfwyv9eoX9Ifx4a+FCDyqqrnA8XkDNvHv7XXUv4Sy+pJKC0eKotqzhdTomRL/4+yVUDIukcUuEiMXM5bJwD0cMhuv7TUJhtZp766ykEgpkXz0SvabyDs7mffEr222/jd+UVtJ82TXUHKa2CahEoTvfRhmOUW6zVtAY+h5IMuO6jBpU/b9c8EnISeHPEm0T5RtW8gZPkffklWTNn4jtuHBHTp6sDw0qroX7OKE5VWGZm8eYTTOgfUXmYaasZNs2xXzcQPbze5e/M3MmivYu4rvt1jI0e64SIa6dgxfdk/vdVfEaPIvKNWQid+g2ltB4qEShO9eXWE5SZrNw3osok7PtWQFEqDHus3tcNlJpLeW7jc0T6RPLkkCedEG3tFK9fT/oLL+B94QVEvvWWuk5AaXXUzxrFaUwWG59uSmJYtxB6R/idXiCl/bqBkJ7QbXS9y5+1bRbppel8Ou5TvPUNH6CuNsp27CD1kUfxiIkh8p131RXDSqukWgSK06zak0ZWsZG7hneuvOD4X5ARbx9mup4HV9efXM/yI8u5o88dDGw30AnR1qz80GH7/MLt29NhwYdofRon+ShKY1OJQHEKKSUfbzhGjzAfRlQcYRRg83vgHQr9b6pX2Xnlebyy+RV6BvbkgQHOmbOgJqaUVJLvuguNhwcdF36MLiioUepVlKagEoHiFBsTcziYUcxdw7tUPq8++xAcWQtDp4C+fsM/zNw6kyJTEdOHT0evdX3/vCU/n+TJk7EZjXT4+CM1lLTS6qljBIpTfLThOKG+7lw1oMowD5vfA50nxE2uV7l/Jv/J6uOruT/2fnoE9nBCpOdmMxpJuf8BzOnpdPz0Ezx6uL5ORWlqqkWgNNihjGL+OpzN/13QCXddhXPrS7Jgzzcw4GbwDq5zuSWmEqb9PY1uAd24q99dToy4etJmI/3ZZzHs2kXErJl4DRrk8joVpTlQLQKlwRZuPIaHXsOt53WqvGDbQrCa4Pz69eu/veNtssuyeXvk243SJZQ9Zy5Fq38m9PHH8Bs3zuX1KUpzoVoESoMUlJlYuTuNawZGEehd4dRKqxl2fALdx0BIt7MXcBbbM7az9PBSJvWeRP/Q/k6MuHr5335L7oIFBNxwA8F3ub71oSjNiUoESoMs25GC0WLjtvOrtAYO/gglmTDk7jqXWW4p55XNrxDpE8mDAx50UqRnV7JpExmv/Afviy4i/KUX1SBySpujuoaUerPZJIv/PkFcp8DKF5ABbP0YAjpBt1F1Lnf+nvmcKDrBgjEL8NJ71bxBAxiPHSP14Udw79qVyLlz1FXDSpukWgRKvW1IzOFEbhm3XVClNZB1AE5shCGTQVO3gdkO5R3is32fcU23a7gg4gInRnsma2EhKffdj3Bzo8MH89H6+NS8kaK0Qi5NBEKIcUKIQ0KIRCFEtdNRCSFuFELsF0LsE0J85cp4FOdavPkEIT5ujOsbXnnBtoWgdYcBk+pUnk3aeG3La/i5+fF43ONOjPRM0mIh9bHHMaWlEfXuO+gjGm92M0VpblzWNSSE0ALzgDFACrBNCLFKSrm/wjrdgWeBi6SU+UKIdq6KR3GulPwyfj+YyX0ju1Y+ZdRYDHuWQN9r63zK6A9Hf2BX1i6mXTgNf3d/J0dcWdbsNyndtInw/07Da/Bgl9alKM2dK1sEQ4FEKeUxKaUJWAJcVWWdu4F5Usp8AClllgvjUZzo660nAbh5aMfKC+KXgqkYhtTtzJtCYyFv7XiL2NBYrupW9WviXAUrvifv008JvPVWAm+4waV1KUpL4MpEEAkkV3ie4nitoh5ADyHEJiHE30KIak/eFkJMEUJsF0Jsz87OdlG4Sm0ZLVa+2ZbMJb3CiAqscDBXSnu3UPtYiKzbr+z3dr1HgbGAF85/AY1w3dfSsHs3GS+9hNf55xP2zNMuq0dRWpKmPlisA7oDI4GbgY+EEAFVV5JSLpBSxkkp40JDQ6suVhrZmn2Z5JSYmHR+ldZA8hbI2mcfTqIOp2Duz93P0sNLmdhzIr2Cejk52tMs2dmkPDQVXXg4kW+reQUU5RRXJoJUoEOF51GO1ypKAVZJKc1SyuPAYeyJQWnGlm5LJirQk4u7V0nKOxeDmy/0u77WZdmkjdf+fo1A90AeGOi6kUWl2Uzqo49hLS4m6r130QUGuqwuRWlpXJkItgHdhRCdhRBuwERgVZV1vsfeGkAIEYK9q+iYC2NSGig5r4yNiTncMLgDGk2FX/3GYvssZH2vBbfaj9u/4sgK4nPieTzucfzc/GreoJ6y3p5D2fbttJ/2Hzx69nRZPYrSErksEUgpLcCDwBrgALBUSrlPCDFNCHGlY7U1QK4QYj+wHnhSSpnrqpiUhvt2RwpCwPVxVSaN37sczKUw8LZal1VoLGTOzjkMajeICV0mODnS04rWrCVv0SICb7kZ/yuvrHkDRWljXHplsZRyNbC6ymsvVXgsgcccN6WZs9oky7YnM6xbCJEBnpUX7loMob0gKq7W5X0Y/yFFpiKeO+85lw3rYDx2nPTnnsOjf3/aPVPtpSyK0uY19cFipQXZlJhDWmE5Nw3pUHlB1kFI2WZvDdRyh55UmMTXB77mmm7X0DPINV01ttJSUqY+hHBzI2ruHDXfsKKchRprSKm1b7YnE+ClZ0zvsMoLdi0GjQ5iJ9a6rDe3v4m7zp0HB7pmUDkpJekvvYzp2HE6fvwR+vbtXVKPorQGqkWg1Ep+qYl1+zK5ekBk5SuJLSb7lcQ9LwPvkFqVtTltM3+k/MHd/e4mxLN229Q53i+/ouinnwidOhXvCy90SR2K0lqoRKDUyopdqZistjO7hQ7/AmU5MPD2WpVjtVl5Y/sbRPpEMql33cYiqi1DQgKZM2fiM3IkwVPqPgy2orQ1KhEoNZJSsnR7Mv2j/IlpX+UUz12LwTei1sNNL09czpH8Izw2+DHcte5Oj9VaXEzqY4+jCw0hYuYMhEZ9xRWlJuq/RKnR/vQiDmYUc8PgKqeMlmRB4q/2YwO1GG662FTMe7veY1C7QYzpNMbpcUopyXj5ZcxpaUTOfhOtv2sHrlOU1kIdLFZqtGJnKnqtYEL/KkM1JywDaav1QeKP4j8ivzyf90e/75LTRQu+/dY+5/Cjj+I1aKDTy1eU1kq1CJRzstokK/ekMbJnu8pzEgPEfwPtB0Bozad/ppWk8cWBL7ii6xX0Ce7j9DjLDx8m87XpeF94IcF3qzmHFaUuat0iEELEAcOBCMAA7AXWnRpCWmmd/nc0h+xiI9cMrDJwbPYhSN8Nl75eq3Lm7Z6HQPDQwIecHqPNYCD1scfQ+PoSMWumOi6gKHVU43+MEOIOIcRO7BPIeAKHgCxgGPCrEOIzIUTHc5WhtFwrdqbi66Hjkl5V5gyK/waEtlYDzB3KO8QPR3/g1phbCfcOr3H9usqcPh3T0WNEzpqJLsQ1p6MqSmtWmxaBF/YZxAzVLRRCDMA+YuhJZwamNL0yk4Vf9mVwZWwEHvoKB4NtNoj/Frr+C3xqnlTunV3v4OPmw+R+k50eY+FPP1Hw7TKCp0xR1wsoSj3V2CKQUs47RxJwk1LullL+5vzQlKa2dl8mZSbrmd1CJzdD4Unof1ONZWzL2MZfKX9xV7+7nD79pCk5mYyXXsZz4EBCpzq/y0lR2opad6YKIf4QQkRXeD4U+1DTSiu1YlcqkQGeDIkOqrwg/hvQe0Ovy8+5vZSSOTvm0M6rHbf0usWpsUmLhbSnngaNhsjZbyB06gQ4Ramvuvz3vA78IoR4B/uUk5cBd7gkKqXJZRcb2XAkm/tGdq0874C5HPZ9DzFX1DjvwG8nfyM+J55pF07DQ+fh1PhyFizAsGsXEW+8gT6y6gyoiqLURa0TgZRyjRDiXmAdkAMMlFJmuCwypUmt2pOGTXJmt9CRNWAshNhzdwtZbBbm7pxLV/+uXNH1CqfGZoiPJ2fe+/hNmID/Fa6bx0BR2oq6dA29CLwLXAy8AvwhhDh334DSYn2/K5V+kf50a+dbeUH8UvAJh84jzrn9isQVJBUl8fCgh9FpnNdtYystJfXJJ9GFtSP8pRedVq6itGV1OeE6GBgqpdwspfwQuBR4xDVhKU0pMauYhNRCrq7aGijLg8Nr7KeMnmNICYPFwPzd8xnYbiAjO4x0amyZM2ZgPplMxIwZaP1cN7WlorQltU4EUspHKp49JKU8IaV0/oAxSpNbsSsVrUZwZWyVISX2rQCbucazhb488CXZhmweHfyoU4eSKP71V/uponfdhffQoU4rV1HautpcUPaREKLfWZZ5CyHuFELc6vzQlKYgpeT7XWkM6xZCqG+V0UHjl0JoDIRX+3UA7PMQL0xYyMgOIxnYznnj/Zizskh/4UU8evcm9CHXTGajKG1VbTpv5wEvOpLBXiAb8MB+EZkfsAj40mURKo1q58kCUgsMPD62R+UFhSmQ/Ddc8uI5p6P8dN+nlJpLnTqUhJSS9Oeex1ZeTsTsNxBqyklFcaoaE4GUcjdwoxDCB4gD2mMfa+iAlPKQi+NTGtmP8Wm46TRnTke5b4X9vu+1Z902rzyPLw98ybjocfQI7HHW9eoq/4svKd24kfCXX8K9Sxenlasoil2NiUAI0VFKeVJKWQL84fqQlKZis0lWJ6Qzskcovh76ygv3LrePNBp09h3xooRFGK1G7htwn9NiMiYmkjV7Nj4jRhAwsfZzIiuKUnu1OVj8/akHQojvXBiL0sS2JeWRWWRkQtWDxHnHIG3nOVsDWWVZLDm0hAldJtDZv7NT4pFmM2lPPY3G25v2r73qkjkMFEWp3TGCiv99ql3eiv0Qn4anXsvomCoDyZ3qFupzzVm3/TjhY6w2K/fG3uu0eHIWLKB8/34i35mrRhVVFBeqTYtAnuWx0opYrDZ+Tsjgkph2eLlV+X2wdwVEDYWA6kcbTy9JZ9nhZVzd/Wo6+Haodp26MuzbR878D/C74gr8xo51SpmKolSvNi2CWCFEEfaWgafjMY7nUkqpruppBf4+lkduqYkr+revvCD7MGQmwLgZZ932w/gPAbin/z1OicVmMpH+zLPoAgMJf+F5p5SpKMrZ1easoZpnJVdavB/j0/B20zKyZ9VuoeWAgN5XV7tdclEyKxNXckPPG5w26UzOu+9hPHKEDh9+oCagV5RGoOb0UzBZbPy8N4OxfcIrT0Ajpf1soU4XgV/7arf9IP4DtBotd/e72ymxGHbvJnfhQvyvvw6fEecez0hRFOdQiUBhU2IOhQYzE6p2C2Xth5xD0Lf6g8THCo/x47EfmdhzIqFeoQ2Ow2YwkPbMs+jCwwh75pkGl6coSu2o2TwUfohPw89Dx/DuVXbme78DoYGYq6rd7oPdH+CudefOfnc6JY6st9/GlJREx08/Qevj45QyFUWpmWoRtHHlZivr9mVyaZ9w3HQVvg6nuoU6Xww+Z/7aP5x/mJ+TfmZSzCSCPILOWF5XpVu2kv/5YgJvvRXv889vcHmKotSeSxOBEGKcEOKQECJRCHHWtr4Q4johhBRCxLkyHuVMfx7Optho4YqqF5Gl74b849D3umq3e3/3+/joffi/Pv/X4BisJaWkP/cc+k4daff4Yw0uT1GUunFZIhBCaLEPWHcZ0Bu4WQjRu5r1fIGHgS2uikU5ux/j0wnyduPCrsGVF+z9DjQ66HXmDGD7cvfx28nfuL337U6ZkD7rjTcwp6UR8frraLy8Glyeoih148oWwVAgUUp5TEppApYA1XU2/xeYCZS7MBalGgaTld8OZDKubzg6bZVuoX3fQ9dLwOvMbp95u+bh7+7PpN6TGhxDyYaNFHzzDUF33IHXoEENLk9RlLpzZSKIBJIrPE9xvPYPIcQgoIOU8qdzFSSEmCKE2C6E2J6dne38SNuo3w9mUWaynnm2UMo2KEyGPmeOLbQ7azcbUjfw7z7/xtfN94zldWEtKiL9hRdw69qV0IenNqgsRVHqr8kOFgshNMBbwOM1rSulXCCljJNSxoWGNvw0RcXuhz1phPq6c17nqt1Cy0HrDr3Gn7HN+7vfJ8gjiFt63dLg+jNfm44lJ4eIGTPQuLvXvIGiKC7hykSQClQceCbK8dopvkBf4A8hRBJwPrBKHTBuHCVGC+sPZXF5v/ZoNRXGFbRZ7YPMdR8DHpX7/3dk7mBz+mbu7HsnXvqG9eUX//YbhStXEnLPFDz79W1QWYqiNIwrE8E2oLsQorMQwg2YCKw6tVBKWSilDJFSRkspo4G/gSullNtdGJPi8Ov+TIwW25ndQic3Q0lGtSONvr/7fUI8Q7ix540NqtuSn0/6Sy/jHhNDyL3OG61UUZT6cVkikFJagAeBNcABYKmUcp8QYpoQ4kpX1avUzo/xabT392BQx8DKC/YuB50n9BhX6eWt6VvZmrGVyX0n46nzbFDdGdOmYS0qImLG62raSUVpBlx6ZbGUcjWwusprL51l3ZGujEU5rdBg5s/D2fz7wmg0FbuFrBbYvxJ6jgP301f2SimZt3se7TzbcX2P6xtUd9Hq1RT//AuhjzyCR8+eDSpLcR2z2UxKSgrl5epkvpbGw8ODqKgo9Hp9zSs7qCEm2qA1+zIwWyUT+le5iCzpLyjLOeNsoS0ZW9iZtZPnznsOD51Hveu1ZGeT8Z9pePTvT/Bdk+tdjuJ6KSkp+Pr6Eh0drWaGa0GklOTm5pKSkkLnzrWfKVANMdEG/RifTscgL/pHVbkYbO9ycPOxHyh2kFIyb9c8wrzCuK579VcZ14aUkvSXXsZWXm7vEtKp3yDNWXl5OcHBwSoJtDBCCIKDg+vcklOJoI3JKzWxKTGHy/u3r/xPbjHBgR+g53jQnz4G8L+0/7E7ezdT+k/BTVv//vzC71dSsn49oY8+gnsXNeNpS6CSQMtUn89NJYI25pe9GVht8syzhY79AeUFlSaoP3VsIMI7gmu6nX2+4pqY09PJfO01POMGE3T77fUuR1EU11CJoI35MT6NLiHe9G5fZYbRfcvB3d8+rITDhtQNJOQkMKX/FPTa2h94qkhKSfrzzyNtNiKmT0do1FdOUZob9V/ZhmQXG/n7WC4TqnYLmcvh4E8QMwF09it8T7UGIn0iubJb/c/2LViyhNL/bSbsqSdx69ixoW9BURQXUImgDfl5bzo2CROqDjl99DcwFlU6W+iP5D/Yn7ufe2PvRa+pX2vAdPIkmbPewPvCCwm46aaGhK4oDfLLL7/Qs2dPunXrxowZM866XnR0NP369WPAgAHExTXOIAfJycn861//onfv3vTp04e5c+dWu15t30N9qFM32pAf96TTI8yHHmFVBovbuxw8g6CLfY5gm7Qxb/c8Ovp2ZEKXM4ehrg1ptZL27HMInY72019TBx6VJmO1WnnggQdYt24dUVFRDBkyhCuvvJLevc8YFR+A9evXExIS0mjx6XQ63nzzTQYNGkRxcTGDBw9mzJgxleKr63uocwxOKUVp9jIKy9l2Io9HR/eovMBUBod+hn7Xg+M4wO8nf+dQ/iGmD5uOTlO/r0jeZ59j2LGDiJkz0IeHNzR8pQn954d97E8rcmqZvSP8ePmKPjWuN3LkSD744AN69epFbm4uI0aMYO/evXWqa+vWrXTr1o0ujrPVJk6cyMqVKxu8E505cyZ79+4lKyuLgwcP8tBDD/HEE0/UuZz27dvTvr395A1fX19iYmJITU2tFJ+r3sMpKhG0ET8lpCMlZ54tdGQtmEv/OVvoVGsg2i+ayzpfVq+6jEeOkD1nDj6jR+F3pRpNRKm/xMREevSw/3iJj4+nX79+lZYPHz6c4uLiM7abPXs2o0ePBiA1NZUOHU6PfxkVFcWWLdXPgyWEYOzYsQghuOeee5gyZcpZY0tISECr1bJ69WqOHTvGddddd0YiqE18FSUlJbFr1y7OO++8Sq/X5T3Uh0oEbcSP8Wn0bu9Hl9Aqk8LvWw7eodBpGADrTqwjsSCRmcNn1qs1IM1m0p55Fo23N+1feUV1CbUCtfnl7gonTpwgMjISjeNMs/j4ePr3719pnQ0bNji1zo0bNxIZGUlWVhZjxoyhV69eXHzxxdWum5CQwPLly9FqtWi1WoKCzpzEqS7xlZSUcN111zFnzhz8/Pxq3sCJVCJoA1Lyy9h1soCnxlUZ28dYAofXwsBbQavDarMyf/d8uvp35dLoS+tVV86CBZTv20fknDnoGrGfVWl99uzZU2nHv2PHDm6qctJBbX5xR0ZGkpx8eo6slJQUIiMjz9jm1LoA7dq145prrmHr1q3VJgKz2Ux2djZdu3YFqm+t1Da+U+Vdd9113HrrrVx77ZkTQtXlPdSHSgRtwE/x6QBM6FflbKHDv4DF8M/ZQmuS1nC08ChvjHgDrUZb53oM+/aRM/8D/CZMwG9c/RKJopyye/fuf4ZKOHLkCCtXruTVV1+ttE5tfnEPGTKEI0eOcPz4cSIjI1myZAlfffXVGeuVlpZis9nw9fWltLSUtWvX8tJL9jEyR40axeeff/7PzvfQoUPExMRUijU2NvaMMmsTn5SSyZMnExMTw2OPPdag91Bf6vTRNuDH+HRio/zpGFxlMpm9y8G3PXS8wN4a2DOfbgHdGNtpbJ3rsJlMpD/zDLrAQMJfeN5JkStt2Z49e7DZbMTGxjJt2jR69+7NZ599VudydDod7733HpdeeikxMTHceOON9Olzurtr/PjxpKWlkZmZybBhw4iNjWXo0KFcfvnljBs3DpvNRmJiYqWun4SEhEqtlbMlgtrYtGkTixcv5vfff2fAgAEMGDCA1atXV4qtpvfQYFLKFnUbPHiwVGrveHaJ7PT0j3LBn0crLzAUSDktRMrVT0sppVyZuFL2/bSvXJu0tl71ZM6eLff37CWL//yzoSErzcD+/fubOgTZrVs3WVRU1NRhyISEBPnoo4+ec52uXbvKsrKyRoqoZtV9fsB2eZb9qmoRtHI/Jdi7hS6verbQwdVgNUHfazFbzby/+31igmIY1XFUneso27mL3IWLCLjhBnzOcmBNUeqiuLgYIQS+vr41r+xiffv25a233jrr8oKCAtzc3PD0bNiETU1JJYJW7oc9aQzuFEhEQJUv6b7l4N8Boobw3ZHvSC1JZeqgqWhE3b4S1pJS0p55Bn14OO2eftqJkSttma+vL4cPH27qMGolICCA/fv3N3UYDaISQSuWmFXCwYziM68dKMuDo79Dn6spsxj4MP5DBocN5qKIi+pcR+aM1zEnJxMxayZaH28nRa4oSmNSiaAVW7UnDSFgfL8qieDAKrBZoM+1fH3wa3IMOUwdOLXO5/wXrVtH4bLvCL77brwaaVwWRVGcTyWCVkpKyfe7UrmoawhhflWml4xfCsHdKQrpyqK9ixgeOZxBYYPqVL45K4uMF1/Co3dvQh98wImRK4rS2FQiaKV2JRdwMq+MqwZUuXagIBlObIL+N/Lpvs8oMhUxddDUOpUtbTbSn3vePu3k7DcQbvWfuUxRlKanEkEr9f2uVNx1Gsb1rTLg295lAOT0GM0XB75gXPQ4egX1qlPZ+V9+RenGjYQ9/ZSadlJRWgGVCFohs9XGj/HpjO4dhq9HlbkE4pdC1FAWJq/DZDXxwIC6desYjxwha/ZsvEdcTMDEiU6MWlGUpqISQSu08UgOeaUmrh5QZSySjL2QtZ+0mHF8c+gbru52NdH+0bUu12YykfrkU2i8vYl4Tc0xoCithRprqBVasSuVAC89I3qEVl6QsBQ0Ot4xpaARGu6NvbdO5WbPmYvx4EGi3p+nBpRTlFZEtQhamRKjhbX7M7i8X3vcdBU+XpsNSYa8TAAAFa9JREFUEpaxr8uF/HTyV27rfRvh3rWfMKZkwwbyFi0i4Kab8L3kkpo3UJRm5P/bu/PwqMqzj+PfO5NJggmQBAhLQiAKyL6YGEAoiq9YBKvS4lZpEbSgL1gpRUFRsWjFpVKpIMomKC516YKACK+IC7IFDGsEEQNkIYSEEMg6mXnePzJgEgIJMEsyc3+uK1dmznnmeB8Pk3vOMuc3evRooqKi6Nq163nHuTMO8lzqQlSlNgIfs3r3EYptDm7rVeWw0MH1mPx0Xm4AkSGR3Nf1vlov03b0KBmTpxDcvj3NH5vi4oqVcr97772XVatWnXfM6TjITz/9lD179vDee+955BvDp6Mq9+zZw8aNG5kzZ85Z/11316aHhnzMf5IziIloQHxsROUZOz/gy0YRbDmVyuO9HycsKKz6BVRh7HYyHp2Mo7CQ6LeWEBASUvOLlG/5dAoc2enaZbboBjfV/KnWFVGVAAMGDCA1NfW8Yy40DlKjKlWdlJVfzDc/ZPPgdVcQEFDhRK6tiLLd/2VmTAvahjVneIfhtV5mzrx5FG7cSMu/Pktwu3ZuqFqpc3NFVGVtXWgcpEZVqjrp421pOAwMj29deUbKJ/wryM5PjiJeiZ+ANcBa/QKqKExKIvvV2TQaOpTG1aQmKT9Ri0/u7uCNqMoLoVGVtSQig4FZgAVYYIx5vsr8icD9QBmQDYw2xhx0Z02+yhjDh0lpJMZFEte08s3fTm1bzJzISK6KuorrW9fuRG/Z8eOkT3oEa0wMLf6i2cPK81wVVVlbFxIHqVGVtSQiFmAOMAhIA7aIyDJjTMUzHN8BCcaYQhF5EHgRuPPspamabEk9zk/HChg3sMrhm9yfWHRiF7nhjZmdMKlWf9CNw0HmY49TlpND2/fewxJWu/MJSrmSq6Iqa+t8cZAaVXnxEoH9xpgDxphS4H3g1ooDjDFfGGMKnU83AjFurMen/XPLYcKCAxnSrfIloUeS5vN2o4bcFH0t3Zqd/YmlOjlvvMGpdetoPnkyDbq6MA5PqQvgqqhKgLvvvpu+ffuyd+9eYmJiWLhw4Zl5NcVBalTlJfwAwyk/HHT6+e+A2ecZPxt44hzzxgBJQFJsbOylZbj5oPyiUtPxiU/NlI+3V55hLzMT3+hs4hd3M2kn02q1rJNffW32dOxk0h55xDgcDjdUq+oDjar8mUZVeoiIjAASgJeqm2+MmWeMSTDGJDRr1qy6IX5txY5Mimx2bk+ofJJ407Z5rA4O4L5WA4kOq/l4YmlaOhmTJhHcoQMt//IXPS+gvEajKj3LnSeL04GKf5linNMqEZEbgKnAtcaYEjfW47P+mXSY9lFh9GodfmaazWHj+T2LiC5zMOoX02tchqOkhPSHH8Y4HMT8YxYB9fgftar/NKrSs9y5R7AFaC8icSISBNwFLKs4QER6AW8AtxhjjrqxFp/1Q9ZJvjuUxx0JrSt9gn9/x0L2m2IebZJISEjjGpeT9eyzFO/eTasXnieoTRt3lqyUqmPc1giMMWXAeOAzIAX4wBizW0Smi8gtzmEvAWHAhyKSLCLLzrE4dQ4fJB0mMEAYdtXPh36OFR3jtZ3z6VdYxMA+f65xGcc/+IC8Dz+iyQNj9T5CSvkht36PwBizElhZZdpTFR5f2IW+qpJim50Pt6YxqHNzmoYFn5n+961/p9heyuTAaKTF+W+yVbhlC0eeeZbQ/v1p9tBD7i5ZKVUH1YmTxeriLNueQV6hjd/3bXtm2sbMjSz7cRkjT5wgLvHB876+9PBh0h76I0ExMUTPfBmxWNxcsVKqLtJGUE8ZY3hrQyodmofR5/Ly65uLy4qZvmE6sRLEAyVW6HzbOV9vP3WKww8+iDGG1nNfw+Lhr7QrpeoObQT11HeH89iVns/v+rY9c5J47va5HD55mGmZ6YRc9XuwVn+nUGO3k/HnSZSmHiRm1isEtW3rwcqVUnWNNoJ66u0NB2kYHMivnbkDKTkpLNm9hF9f1pbE4hJIGF3t64wxZP31OU59+SUtnphKaJ8+nixbKVUHaSOoh7JPlrBiRya/iY8hNDiQMkcZT294mvDgxkxM3QNXDoHw1tW+NmfBAo6/+y6Ro0cToeHzSim0EdRLS75NxeZwMPKatuXPdy9hT84eHmval8YFx6B39VnEJz75hOyXZ9JoyBCiJtV8WalSvqK2UZVt27alW7du9OzZk4SEBI/UVlxcTGJiIj169KBLly5Mmzat2nEaVanOKCgp460NqfyycwvimoaSkpPC7OTZDIq9gRt3r4LoeGjb/+zXbdhAxuNTuax3b1o+PwMJ0E2v/EdtoipP++KLL0hOTiYpKcnNVZULDg5m7dq1bN++neTkZFatWsXGjRsrjdGoSlXJ+1sOk19cxthrL6e4rJgpX08hIjiCp5okIrmL4Panoco9goq2bydt3HiC4+KIefUfBAQFead4VS+9sPkFvs/93qXL7BjZkcmJk2sc58moygvlqqhKESHMeat3m82GzWY76z5fGlWpzrDZHSz8+gCJcZH0io3g+c3Pc+DEAd644XXCVz4JEXHQ6VeVXlOcksKhP4zB0rQprefP18tEVb3iyahKKP+jfOONNyIijB07ljFjxpxzrCujKu12O/Hx8ezfv59x48ZpVKU6t+U7Msg4Ucwzt3Vlffp63kl5hxGdRnCNzUD6Vhg6EwJ+/lJYyY8/cmj0fQSEhtLmzUVYm0d5sXpVX9Xmk7s7eCOq8ptvviE6OpqjR48yaNAgOnbsyIABA6od68qoSovFQnJyMnl5eQwbNoxdu3bVeD7DlbQR1BNldgevrt3Plc0b0j3Wwp0rnqRdeDsevuphWDocQqOg52/PjC89dIhD944Ci6W8Cbgw1k4pT/B0VCVwJoEsKiqKYcOGsXnz5mobgaujKk8LDw9n4MCBrFq1qlIjqLdRlcq1/pOcwYHsAl67pyePrZ9Cfmk+c2+YS8ihTZD6NQx+Aazlt44uOXCAQ/eOwthsxL61RL8wpuolT0dVFhQU4HA4aNiwIQUFBaxevZqnniq/NZo7oyqzs7OxWq2Eh4dTVFTEmjVrmDy58l5YfY6qVC5iszuY9fk+ukY3Yr/tYzZlbmJq76lcGdEB1j4LjaIh/l4Air//noMjfodxOIhdsoQQ5/FVpeobT0dVZmVl0b9/f3r06EFiYiJDhw5l8ODBbo+qzMzMZODAgXTv3p2rr76aQYMGcfPNN1eqrd5GVbrrJz4+/iKC2+q3dzYeNG0mLzezvv236bq4q5m2flr5jL2fGTOtkTFbFhpjjCncvt18n9jb7Lv2OlN84ID3Clb1nkZV/kyjKpXXFZXaeXXtD3RpU8A7B2bQKbITj/V+DBx2WPsMhLeBniMo+PZbDo0ajaVRI9osXUpwXJy3S1fqomlUpWdpI6jj5n75I0dOZVMQMY8waxivXv8qwZZg+O5tOLIDrn+SvP9+wqExY7G2akWbpW8TFKMnhlX9plGVnqUni+uww7mFvP5VCq2ufJ+CshMsvmkxzUObQ9Fx+Hw6pnUfstdmkPP664T260f0rFewOL+YopRStaWNoA6bvnwH1pZLOWV+4uVfvEyXJs6TQ1/MwHEyj8zU1uSvfZ3w24fT4qmnEKvVuwUrpeolbQR11Bd7j/B13mysjb/nyb7TuKGN85rjI7so+fxN0rdeQUnWtzSbOJEmf7j/rK+kK6VUbWkjqINyC4qZuHYq1sbbeajnwwzvMLx8Rlkp+S+OInNNEyTUQusF8wnr18+7xSql6j1tBHVMmaOMO//1R8ou28ywtqMY0+N+AByFhRydcBfHvzpFSPs2xLzxJtZWrbxcrVLKF2gjqEOKy4q5d/kEjjg2EN/wHqZfOxGAwm3fkTHpT9gysojoF0PU3E/0DqJKKZfRRlBH5BTlMHb1ePbm7Say5Hbmj3gUR1ERx+bMIWfRm1gvcxD7qwBCn/kYtAkopVxIG0EdsDd3Lw+t/SOZp7IJyB7JW6MeoGjNGrJefIGyjEzCu4cR1ekQljErIURvI62Uci39QpkXGWP4aN9H3LPyHrJPFVB8aCzz4wfApPGkT5iAJawhsWN60bLzPiy3v1aePqaUuiC1jYIE98ZB1sRut9OrV68z9xnyZG26R+AlucW5PLfpOT5L/YwI6Qpb+zHnZDIN35tFcaNGNH/iCSKa7EY2zYEBj0LX33i7ZKXqpdNRkGFhYdhsNvr3789NN91Enz59Ko07HQe5Zs0aYmJiuPrqq7nllltclgJWk1mzZtGpUyfy8/PPmufu2rQReJgxhhU/reCFzS9wynaKPtnX02vVUa5Lf42ABg2IvP9+IkeNInDjDNj0BiSOgese83bZyo8dee45SlJcG1UZ3KkjLR5/vMZxroiqrE0UJFx4HKSroiqhPF9gxYoVTJ06tdr7GmlUpQ9JPprMzK0z2ZG5jeFpMSR8GUqb9NXYg0Noev/9RI4eRWBoMCz/E+x4H/qOhxufPSuDWCl/4aqoypqiIOHC4yBdGVU5YcIEXnzxxWrHXkxtF0obgQckH01m4Y4FZG5ax/X7gpiYEkxw/kEyQpuSfvcYBk4YjaVxY8jeBwtGwtEUGDgVBjyiTUB5XW0+ubuDK6Mq3REF6aqoyuXLlxMVFUV8fDzr1q27pJouljYCNym0FbJm73K2fbaUht/t5859QpN8B2WBZWyO6sim3gMYN+m39GgTCbYi+PIl+GZmecrYiI+h3f94exWU8ip3RFWeKwoSLiwO0pVRlevXr2fZsmWsXLmS4uJi8vPzGTFiBEuXLr2o2i6GNgIXOnEsnW3rPiRj8zqCdu6nwyE7HezgsFrI7diTOQ07szbySm7u046XhnSmsaUEkhbB1zPhxGHo9KvyyMnGehtppVwVVVmbKEg4fxykO6MqZ8yYwYwZMwBYt24df/vb3yo1gZpqcwW3NgIRGQzMAizAAmPM81XmBwNvAfFADnCnMSbVnTVdKmMM9uPHKUo/TFbKNo7u2UbhD3sJPphF02OltACigBOtGnHkhnh2NUtkcUEEJ00gfS9vwoe/bEuXku2wdjHs+BBKT0LLnnDbXIj7hZfXTqm6Y/v27YSEhNCjRw+6d+9+JqryySefvKDlZGZmMnLkSOx2Ow6HgzvuuKPSJZpDhgxhwYIFtGrV6kwcpN1uZ/To0XTp0qXWUZVDhw699JWuoqbaXEXKE8xcT0QswD5gEJAGbAHuNsbsqTDmf4HuxpgHROQuYJgx5s5qF+iUkJBgkpKSLrouYwymtBRTUoIpKcHh/G0vKqKsuJCiEzkU5WZTnJtNaV4utrw8bMdzceTkEJidR4OcAgJtjjPLswtkRQZwvHkjclu3ISuqG/sD23HwRCGhpoArgvO4sVUp/RofIyL/eziyC+wlYL2sfA8g4T5onajnAlSdkpKSUukTrze0b9+ebdu2eT2lbNeuXSxatOi8KWXt2rVj586ddSalrLrtJyJbjTEJ1Y135x5BIrDfGHPAWcT7wK1AxSifW4GnnY8/AmaLiBg3dKePHr6Z9p//SFBZ7V9jBWzBUBICeaFworFQFG0woYaQBg4aXmanVQM7Xe1lBJHHZewj2PYp2ICKd4HIAHLDoUU36D0GrrgeYq8Ba4hrV1IpH6FRlZ7lzkYQDRyu8DwNqHrN1pkxxpgyETkBNAGOVRwkImOAMQCxsbEXVYzEtWFHzzQcFsERKDgsggkUHJYATKBgLAKBAWC1EBAciCXEijXISnBAIA2NlWbGSjQWEMFqEQItFgItAVgtFkoDLQQGWzFhjSAsHIIbQXBY+e/G0dAoBkKb6qd+pWpJoyo9q16cLDbGzAPmQfmhoYtZxm8mzHFpTUop5Svcea+hdKB1hecxzmnVjhGRQKAx5SeNlVJKeYg7G8EWoL2IxIlIEHAXsKzKmGXASOfj4cBad5wfUEpdOH0r1k8Xs93c1giMMWXAeOAzIAX4wBizW0Smi8gtzmELgSYish+YCExxVz1KqdoLCQkhJydHm0E9Y4whJyeHkJALuxDFbZePusulXj6qlKqZzWYjLS3tzBe6VP0REhJCTEwMVqu10nRvXT6qlKqnrFYrcXFx3i5DeYgG0yillJ/TRqCUUn5OG4FSSvm5eneyWESygYMX+fKmVPnWsh/wx3UG/1xvXWf/cLHr3MYY06y6GfWuEVwKEUk611lzX+WP6wz+ud66zv7BHeush4aUUsrPaSNQSik/52+NYJ63C/ACf1xn8M/11nX2Dy5fZ786R6CUUups/rZHoJRSqgptBEop5ef8rhGIyNMiki4iyc6fId6uyV1EZLCI7BWR/SLiF3d2FZFUEdnp3LY+eXdCEVkkIkdFZFeFaZEiskZEfnD+jvBmje5wjvX22feziLQWkS9EZI+I7BaRh53TXb6t/a4ROP3dGNPT+bPS28W4g4hYgDnATUBn4G4R6ezdqjxmoHPb+ur15YuBwVWmTQE+N8a0Bz7HN2/pvpiz1xt89/1cBvzZGNMZ6AOMc76HXb6t/bUR+INEYL8x5oAxphR4H7jVyzUpFzDGfAXkVpl8K7DE+XgJcJtHi/KAc6y3zzLGZBpjtjkfn6Q81yUaN2xrf20E40Vkh3NX0+d2oZ2igcMVnqc5p/k6A6wWka0iMsbbxXhQc2NMpvPxEaC5N4vxMJ9/P4tIW6AXsAk3bGufbAQi8n8isquan1uBucAVQE8gE3jZq8UqV+tvjLmK8kNi40RkgLcL8jRn3Ku/XBfu8+9nEQkDPgYmGGPyK85z1bb2yWAaY8wNtRknIvOB5W4ux1vSgdYVnsc4p/k0Y0y68/dREfk35YfIvvJuVR6RJSItjTGZItISOOrtgjzBGJN1+rEvvp9FxEp5E3jHGPMv52SXb2uf3CM4H+f/uNOGAbvONbae2wK0F5E4EQkC7gKWebkmtxKRUBFpePoxcCO+u32rWgaMdD4eCfzXi7V4jC+/n0VEKM91TzHGzKwwy+Xb2u++WSwib1O+G2mAVGBsheNtPsV5Kd0rgAVYZIz5q5dLcisRuRz4t/NpIPCuL66ziLwHXEf57YizgGnAf4APgFjKb9N+hzHGp06snmO9r8NH388i0h/4GtgJOJyTH6f8PIFLt7XfNQKllFKV+d2hIaWUUpVpI1BKKT+njUAppfycNgKllPJz2giUUsrPaSNQSik/p41AKaX8nDYCpS6RiFztvOlZiPPbzbtFpKu361KqtvQLZUq5gIg8C4QADYA0Y8wML5ekVK1pI1DKBZz3c9oCFAPXGGPsXi5JqVrTQ0NKuUYTIAxoSPmegVL1hu4RKOUCIrKM8hS4OKClMWa8l0tSqtZ8Mo9AKU8Skd8DNmPMu86s6G9F5HpjzFpv16ZUbegegVJK+Tk9R6CUUn5OG4FSSvk5bQRKKeXntBEopZSf00aglFJ+ThuBUkr5OW0ESinl5/4fErB2rylu3IgAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# 累積分布関数\n",
"x = np.arange(-5,20,0.1)\n",
"cdf1 = stats.gumbel_r.cdf(x, loc=0.5, scale=2.0)\n",
"cdf2 = stats.gumbel_r.cdf(x, loc=1.0, scale=2.0)\n",
"cdf3 = stats.gumbel_r.cdf(x, loc=1.5, scale=3.0)\n",
"cdf4 = stats.gumbel_r.cdf(x, loc=3.0, scale=4.0)\n",
" \n",
"plt.plot(x, cdf1, label=(r'$\\mu=0.5, \\eta=2.0$'))\n",
"plt.plot(x, cdf2, label=(r'$\\mu=1.0, \\eta=2.0$'))\n",
"plt.plot(x, cdf3, label=(r'$\\mu=1.5, \\eta=3.0$'))\n",
"plt.plot(x, cdf4, label=(r'$\\mu=3.0, \\eta=4.0$'))\n",
"plt.xlabel(\"x\")\n",
"plt.ylabel(\"F(x)\")\n",
"plt.title(\"Cumulative Distribution Function\")\n",
"plt.legend()\n",
"# plt.savefig(\"../img/gumbel_cdf.png\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.6.9"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment