python problem

donghaochen
ps3c-notebook-framework.ipynb

{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Problem Set 3 Notebook" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Loading word list from file...\n", " 83667 words loaded.\n", "Loading word list from file...\n", " 83667 words loaded.\n" ] } ], "source": [ "import random\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "\n", "import ps2\n", "import ps3a" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "def plot_mc_hand(hand, N, n):\n", " \n", " _, mc_scores = ps3.play_n_mc_hand(hand, N, n)\n", " # IMPLEMENT ME" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEICAYAAAB/Dx7IAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAADVtJREFUeJzt3W+MZXV9x/H3RxZEoASVaVMXdTAxttYUMZsW/8RY9AECER/YBqutNDbbpLaitTVY2mj7wECjVm0M6Rb8b9R0JRXBqAQxjU2EzgIqsFqpbnEVy1grim26Ur99cM824zCzc2Z37s58775fyWTnnHvm3u/JWd6e/e29bqoKSVIfj9jsASRJ62O4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOFWO0lOSLI7yb4kleR5yx5PkiuT/Mfw9VdJMvK53zQ8568v2bdt2De/zjl/LcnNSR5Ism89PysdiuFWV58HXg58Z4XHdgIvBs4Cfhm4EPi9dTz394C/THLcEc74I+DdwJ8c4fNIP8Vw66ga7pL/OMmXhjvRjyY5cT3PUVUHqurtVfV54H9XOOQVwFuran9VfQt4K3DJOl7iU8ABJv/DcNiq6taq+gDw9SN5Hmk5w63N8BvAecCZTO6ILwFI8oQk3z/E12+OfP5fAr64ZPuLw76xCvhz4I1Jjl/+YJLLDjXnOl5HOizbNnsAHZPeWVXfBkjyCeDpAFV1L3DaBjz/KcADS7YfAE5Jkhr5f85TVdcluRz4XeDvlj12BXDFBswpHRbvuLUZlq5L/xeT0G6kB4FTl2yfCjw4NtpL/BlwObCupRxp2gy3toxhqeTBQ3y9bORT3cXkLyYPOmvYty5VdSNwD/D7y+b800PNud7XkdbLpRJtGcNSyai77ySPBA6+xe+E4S84/2e4q34/8EdJPslkvfp1wN8s+dl9wJuq6r0jXupy4OPL5nwz8OYRMz4COAE4frKZE4GfVNWBEa8rrco7bnX1VeC/ge3Ap4fvnzg89rfAJ4AvA3cCNwz7SHIC8FjgC2NepKr+Cbj1MGd87jDXJ4EnDN9/5jCfS/p/8R9S0LEkyXOAV1XVSzd7FulwGW5JasalEklqxnBLUjOGW5KamcrbAU8//fSan5+fxlNL0kzas2fPd6tqbsyxUwn3/Pw8CwsL03hqSZpJSf5t7LEulUhSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzW+5fwJm/7IZNed19V1ywKa8rSevlHbckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzhluSmjHcktSM4ZakZgy3JDVjuCWpGcMtSc2MCneS1ya5K8mdST6c5MRpDyZJWtma4U6yHXg1sKOqngYcB1w87cEkSSsbu1SyDXhUkm3AScC3pzeSJOlQ1gx3VX0LeAtwL3Af8EBVfWb5cUl2JllIsrC4uLjxk0qSgHFLJY8GLgLOBB4HnJzk5cuPq6pdVbWjqnbMzc1t/KSSJGDcUskLgG9U1WJV/Ri4FnjWdMeSJK1mTLjvBc5JclKSAM8H9k53LEnSasascd8C7AZuA748/MyuKc8lSVrFtjEHVdUbgTdOeRZJ0gh+clKSmjHcktSM4ZakZgy3JDVjuCWpGcMtSc0YbklqxnBLUjOGW5KaMdyS1IzhlqRmDLckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzhluSmjHcktSM4ZakZgy3JDVjuCWpGcMtSc0YbklqxnBLUjOGW5KaMdyS1IzhlqRmDLckNWO4JamZUeFOclqS3Um+kmRvkmdOezBJ0sq2jTzuHcCnquolSU4ATpriTJKkQ1gz3ElOBZ4LXAJQVQeAA9MdS5K0mjFLJU8CFoH3JLk9ydVJTl5+UJKdSRaSLCwuLm74oJKkiTHh3gY8A7iqqs4GfgRctvygqtpVVTuqasfc3NwGjylJOmhMuPcD+6vqlmF7N5OQS5I2wZrhrqrvAN9M8pRh1/OBu6c6lSRpVWPfVfKHwIeGd5R8Hfid6Y0kSTqUUeGuqjuAHVOeRZI0gp+clKRmDLckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzhluSmjHcktSM4ZakZgy3JDVjuCWpGcMtSc0YbklqxnBLUjOGW5KaMdyS1IzhlqRmDLckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JamZ0uJMcl+T2JNdPcyBJ0qGt5477UmDvtAaRJI0zKtxJzgAuAK6e7jiSpLWMveN+O/B64CerHZBkZ5KFJAuLi4sbMpwk6eHWDHeSC4H7q2rPoY6rql1VtaOqdszNzW3YgJKknzbmjvvZwIuS7AM+Apyb5INTnUqStKo1w11Vb6iqM6pqHrgY+GxVvXzqk0mSVuT7uCWpmW3rObiqPgd8biqTSJJG8Y5bkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzhluSmjHcktSM4ZakZgy3JDVjuCWpGcMtSc0YbklqxnBLUjOGW5KaMdyS1IzhlqRmDLckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzhluSmjHcktSM4ZakZgy3JDVjuCWpmTXDneTxSW5OsjfJXUkuPRqDSZJWtm3EMQ8Br6uq25L8DLAnyY1VdfeUZ5MkrWDNO+6quq+qbhu+/yGwF9g+7cEkSStb1xp3knngbOCWFR7bmWQhycLi4uLGTCdJepjR4U5yCvAx4DVV9YPlj1fVrqraUVU75ubmNnJGSdISo8Kd5Hgm0f5QVV073ZEkSYcy5l0lAa4B9lbV26Y/kiTpUMbccT8b+C3g3CR3DF/nT3kuSdIq1nw7YFV9HshRmEWSNIKfnJSkZgy3JDVjuCWpGcMtSc0YbklqxnBLUjOGW5KaMdyS1IzhlqRmDLckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGbW/BdwJKmb+ctu2JTX3XfFBUfldbzjlqRmDLckNWO4JakZwy1JzRhuSWrGcEtSM4Zbkpox3JLUjOGWpGYMtyQ1Y7glqRnDLUnNGG5JasZwS1IzhluSmjHcktSM4ZakZgy3JDVjuCWpmVHhTnJekq8muSfJZdMeSpK0ujXDneQ44F3AC4GnAi9N8tRpDyZJWtmYO+5fAe6pqq9X1QHgI8BF0x1LkrSabSOO2Q58c8n2fuBXlx+UZCewc9h8MMlXD3Om04HvHubPHrZcebRf8adsyjlvsmPtnI+184Vj8Jxz5RGd8xPHHjgm3FlhXz1sR9UuYNfYF171xZKFqtpxpM/Tiec8+4618wXPeZrGLJXsBx6/ZPsM4NvTGUeStJYx4f5n4MlJzkxyAnAxcN10x5IkrWbNpZKqeijJHwCfBo4D3l1Vd01xpiNebmnIc559x9r5guc8Nal62HK1JGkL85OTktSM4ZakZrZMuI+Fj9UneXySm5PsTXJXkkuH/Y9JcmOSrw2/PnqzZ91oSY5LcnuS64ftM5PcMpzzR4e/+J4ZSU5LsjvJV4br/cxZv85JXjv8vr4zyYeTnDhr1znJu5Pcn+TOJftWvK6ZeOfQtC8lecZGzbElwn0Mfaz+IeB1VfWLwDnAq4bzvAy4qaqeDNw0bM+aS4G9S7avBP56OOf/BF65KVNNzzuAT1XVLwBnMTn3mb3OSbYDrwZ2VNXTmLyR4WJm7zq/Fzhv2b7VrusLgScPXzuBqzZqiC0Rbo6Rj9VX1X1Vddvw/Q+Z/Me8ncm5vm847H3AizdnwulIcgZwAXD1sB3gXGD3cMhMnXOSU4HnAtcAVNWBqvo+M36dmbxL7VFJtgEnAfcxY9e5qv4R+N6y3atd14uA99fEF4DTkvz8RsyxVcK90sfqt2/SLEdFknngbOAW4Oeq6j6YxB342c2bbCreDrwe+Mmw/Vjg+1X10LA9a9f7ScAi8J5heejqJCczw9e5qr4FvAW4l0mwHwD2MNvX+aDVruvUurZVwj3qY/WzIskpwMeA11TVDzZ7nmlKciFwf1XtWbp7hUNn6XpvA54BXFVVZwM/YoaWRVYyrOteBJwJPA44mclSwXKzdJ3XMrXf51sl3MfMx+qTHM8k2h+qqmuH3f9+8I9Qw6/3b9Z8U/Bs4EVJ9jFZAjuXyR34acMfqWH2rvd+YH9V3TJs72YS8lm+zi8AvlFVi1X1Y+Ba4FnM9nU+aLXrOrWubZVwHxMfqx/Wdq8B9lbV25Y8dB3wiuH7VwAfP9qzTUtVvaGqzqiqeSbX9bNV9TLgZuAlw2Gzds7fAb6Z5CnDrucDdzPD15nJEsk5SU4afp8fPOeZvc5LrHZdrwN+e3h3yTnAAweXVI5YVW2JL+B84F+AfwUu3+x5pnSOz2HyR6UvAXcMX+czWfO9Cfja8OtjNnvWKZ3/84Drh++fBNwK3AP8PfDIzZ5vg8/16cDCcK3/AXj0rF9n4C+ArwB3Ah8AHjlr1xn4MJM1/B8zuaN+5WrXlclSybuGpn2ZyTtuNmQOP/IuSc1slaUSSdJIhluSmjHcktSM4ZakZgy3JDVjuCWpGcMtSc38HxKC3d74yqf4AAAAAElFTkSuQmCC\n", "text/plain": [ "<matplotlib.figure.Figure at 0x7ff3e3d674e0>" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAEICAYAAAB/Dx7IAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAADsVJREFUeJzt3X2MZfVdx/H3h2WR52DgpqHAdkqiJLWxQCb4gCEt0AYKoTVWA0ojpmZMrAq2SqhosP5hwNimmhjjSrGtRbDykFSIFJKCDaZAZ2FpoQtNxS1daN1B5GELEaFf/7h34zDMw5ndOXPnN/N+JTd7zz1nzvlchv3sye+ec3+pKiRJ7Thg3AEkSctjcUtSYyxuSWqMxS1JjbG4JakxFrckNcbilqTGWNxa85IclOSmJDuTVJJ3zlmfJNck+a/R48+SpOO+/3i0z1+c9dqBo9cmlpnzXUnuTvJ8kp3zrJ8YrX8pyWNJzl7O/qW9LG614l7gYuD786ybAt4PvAP4SeB84DeWse9ngT9Jsmk/M/4AuA74/QXW3wA8BBwNXAnclGSwn8fUBmRxq1ejs+TfS/L10ZnoPyY5eDn7qKpXqupTVXUv8No8m/wq8Imq2lVVTwGfAC5ZxiHuAF5h+A/DPquqB6rq74En5q5L8uPAqcBVVfVyVd0MfAP4hf05pjYmi1ur4ZeAc4C3MjwjvgQgyZYkzy3y+OWO+/8J4OFZyw+PXuuqgD8Crkqyee7KJFcslnMZGZ+oqhf3I6cEwIHjDqAN4S+r6mmAJP8MnAxQVU8CR63A/g8Hnp+1/DxweJJUxy/jqaovJrkS+HXgb+esuxq4eoUz7s153H7uVxuQZ9xaDbPHpV9iWGIraQ9w5KzlI4E9XUt7lj9kOPa8rKGcjuZmZLT84jzbSouyuDU2o6GSPYs8fqXjrh5l+MHkXu8YvbYsVXUX8G3gN+fk/IPFci4j44lJjtjfnJLFrbGpqier6vBFHtfv3TbJj8z6UPOgJAfPuuTvc8BHkhyX5M3AR4HPzPrZnUku6RjrSuDyOTn/dLGcs45zwCjj5uFiDk5y0Ggf3wK2MxxHPzjJzzMc77+5838wacQxbrXiceAto+dfGv35VmAn8DfAiQyv0gC4dvQao+I8Grivy0Gq6t+SPACcuw8ZzwDunrX8MvCvwDtHyxcy/Aflv4EngQ9U1cw+HEcbXJxIQetZkp8DPlxVF407i7RSLG5Jaoxj3JLUGItbkhpjcUtSY3q5quSYY46piYmJPnYtSevStm3bnqmqTl861ktxT0xMMD093ceuJWldSvKdrts6VCJJjbG4JakxFrckNcbilqTGWNyS1BiLW5Ia06m4k/xukkeTPJLkhuXOGShJWjlLFneS44DfASar6u3AJoZfTylJGoOuQyUHAockORA4FHi6v0iSpMUseedkVT2V5M8ZfvH7y8CdVXXn3O2STAFTAFu2bNnnQBNX3L7PP7s/dl593liOK0nL1WWo5EeB9zGcbeTNwGFJLp67XVVtrarJqpocDDrdbi9J2gddhkrOBv6jqmaq6n+BW4Cf7TeWJGkhXYr7SeCnkxw6mpz1LGBHv7EkSQtZsrir6n7gJuBBhpOxHgBs7TmXJGkBnb7WtaquAq7qOYskqQPvnJSkxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTEWtyQ1xuKWpMZY3JLUGItbkhpjcUtSYyxuSWqMxS1JjbG4JakxFrckNabLZMEnJdk+6/FCkstWI5wk6Y2WnAGnqh4HTgZIsgl4Cri151ySpAUsd6jkLODfq+o7fYSRJC1tucV9IXBDH0EkSd10Lu4kBwEXAP+0wPqpJNNJpmdmZlYqnyRpjuWccZ8LPFhV/znfyqraWlWTVTU5GAxWJp0k6Q2WU9wX4TCJJI1dp+JOcijwbuCWfuNIkpay5OWAAFX1EnB0z1kkSR1456QkNcbilqTGWNyS1BiLW5IaY3FLUmMsbklqjMUtSY2xuCWpMRa3JDXG4pakxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTEWtyQ1puvUZUcluSnJY0l2JPmZvoNJkubXaeoy4C+AO6rqA0kOAg7tMZMkaRFLFneSI4EzgEsAquoV4JV+Y0mSFtJlqOREYAb4uyQPJbk2yWFzN0oylWQ6yfTMzMyKB5UkDXUp7gOBU4G/rqpTgB8AV8zdqKq2VtVkVU0OBoMVjilJ2qtLce8CdlXV/aPlmxgWuSRpDJYs7qr6PvDdJCeNXjoL+GavqSRJC+p6VclvA9ePrih5Avi1/iJJkhbTqbirajsw2XMWSVIH3jkpSY2xuCWpMRa3JDXG4pakxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTEWtyQ1xuKWpMZY3JLUGItbkhpjcUtSYyxuSWqMxS1Jjek0A06SncCLwGvAq1XlbDiSNCZd55wEeFdVPdNbEklSJw6VSFJjuhZ3AXcm2ZZkar4NkkwlmU4yPTMzs3IJJUmv07W4T6+qU4FzgQ8nOWPuBlW1taomq2pyMBisaEhJ0v/rVNxV9fToz93ArcBpfYaSJC1syeJOcliSI/Y+B94DPNJ3MEnS/LpcVfIm4NYke7f/h6q6o9dUkqQFLVncVfUE8I5VyCJJ6sDLASWpMRa3JDXG4pakxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTEWtyQ1xuKWpMZY3JLUGItbkhpjcUtSYyxuSWqMxS1JjbG4JakxnYs7yaYkDyW5rc9AkqTFLeeM+1JgR19BJEnddCruJMcD5wHX9htHkrSUrmfcnwIuB3640AZJppJMJ5memZlZkXCSpDdasriTnA/srqpti21XVVurarKqJgeDwYoFlCS9Xpcz7tOBC5LsBG4Ezkzy+V5TSZIWtGRxV9XHqur4qpoALgS+XFUX955MkjQvr+OWpMYcuJyNq+oe4J5ekkiSOvGMW5IaY3FLUmMsbklqjMUtSY2xuCWpMRa3JDXG4pakxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTEWtyQ1xuKWpMZY3JLUGItbkhrTZZb3g5M8kOThJI8m+fhqBJMkza/L1GX/A5xZVXuSbAbuTfIvVXVfz9kkSfNYsrirqoA9o8XNo0f1GUqStLBOY9xJNiXZDuwG7qqq++fZZirJdJLpmZmZlc4pSRrpVNxV9VpVnQwcD5yW5O3zbLO1qiaranIwGKx0TknSyLKuKqmq54B7gHN6SSNJWlKXq0oGSY4aPT8EOBt4rO9gkqT5dbmq5Fjgs0k2MSz6L1TVbf3GkiQtpMtVJV8HTlmFLJKkDrxzUpIaY3FLUmMsbklqjMUtSY2xuCWpMRa3JDXG4pakxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTEWtyQ1xuKWpMZY3JLUGItbkhrTZc7JE5LcnWRHkkeTXLoawSRJ8+sy5+SrwEer6sEkRwDbktxVVd/sOZskaR5LnnFX1feq6sHR8xeBHcBxfQeTJM1vWWPcSSYYThx8/zzrppJMJ5memZlZmXSSpDfoXNxJDgduBi6rqhfmrq+qrVU1WVWTg8FgJTNKkmbpVNxJNjMs7eur6pZ+I0mSFtPlqpIAnwZ2VNUn+48kSVpMlzPu04EPAmcm2T56vLfnXJKkBSx5OWBV3QtkFbJIkjrwzklJaozFLUmNsbglqTEWtyQ1xuKWpMZY3JLUGItbkhpjcUtSYyxuSWqMxS1JjbG4JakxFrckNcbilqTGWNyS1BiLW5IaY3FLUmMsbklqTJc5J69LsjvJI6sRSJK0uC5n3J8Bzuk5hySpoyWLu6q+Ajy7ClkkSR0sOVlwV0mmgCmALVu2rNRutU5NXHH7WI678+rzxnJcra71/v/Xin04WVVbq2qyqiYHg8FK7VaSNIdXlUhSYyxuSWpMl8sBbwC+CpyUZFeSD/UfS5K0kCU/nKyqi1YjiCSpG4dKJKkxFrckNcbilqTGWNyS1BiLW5IaY3FLUmMsbklqjMUtSY2xuCWpMRa3JDXG4pakxljcktQYi1uSGmNxS1JjLG5JaozFLUmNsbglqTGdijvJOUkeT/LtJFf0HUqStLAuc05uAv4KOBd4G3BRkrf1HUySNL8uZ9ynAd+uqieq6hXgRuB9/caSJC1kycmCgeOA785a3gX81NyNkkwBU6PFPUke38dMxwDP7OPP7rNcs9pHfJ2xvOcx22i/Z3/HG0Cu2a/3/JauG3Yp7szzWr3hhaqtwNauB17wYMl0VU3u735a4nte/zba+wXfc5+6DJXsAk6YtXw88HQ/cSRJS+lS3F8DfizJW5McBFwIfLHfWJKkhSw5VFJVryb5LeBLwCbguqp6tMdM+z3c0iDf8/q30d4v+J57k6o3DFdLktYw75yUpMZY3JLUmDVT3Bvxtvok1yXZneSRcWdZDUlOSHJ3kh1JHk1y6bgz9S3JwUkeSPLw6D1/fNyZVkuSTUkeSnLbuLOshiQ7k3wjyfYk070eay2McY9uq/8W8G6Glx9+Dbioqr451mA9S3IGsAf4XFW9fdx5+pbkWODYqnowyRHANuD96/n3nCTAYVW1J8lm4F7g0qq6b8zRepfkI8AkcGRVnT/uPH1LshOYrKrebzpaK2fcG/K2+qr6CvDsuHOslqr6XlU9OHr+IrCD4Z2561YN7Rktbh49xn+21LMkxwPnAdeOO8t6tFaKe77b6tf1X+iNLskEcApw/3iT9G80ZLAd2A3cVVXr/j0DnwIuB3447iCrqIA7k2wbfQVIb9ZKcXe6rV7rQ5LDgZuBy6rqhXHn6VtVvVZVJzO86/i0JOt6WCzJ+cDuqto27iyr7PSqOpXhN6l+eDQU2ou1UtzeVr9BjMZ5bwaur6pbxp1nNVXVc8A9wDljjtK304ELRmO+NwJnJvn8eCP1r6qeHv25G7iV4RBwL9ZKcXtb/QYw+qDu08COqvrkuPOshiSDJEeNnh8CnA08Nt5U/aqqj1XV8VU1wfDv8per6uIxx+pVksNGH7iT5DDgPUBvV4utieKuqleBvbfV7wC+0PNt9WtCkhuArwInJdmV5EPjztSz04EPMjwD2z56vHfcoXp2LHB3kq8zPEG5q6o2xOVxG8ybgHuTPAw8ANxeVXf0dbA1cTmgJKm7NXHGLUnqzuKWpMZY3JLUGItbkhpjcUtSYyxuSWqMxS1Jjfk/uZzNiP4CSygAAAAASUVORK5CYII=\n", "text/plain": [ "<matplotlib.figure.Figure at 0x7ff3e1bb1f28>" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEICAYAAABYoZ8gAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAFfpJREFUeJzt3X+UZ3V93/HnK8sCbTAB3NFSYFlIaCsmgnSK9uCJ2CisJhU9se0SEzEHz/YYafO7XWILBk970Bxjjg0RN3WDWgUTfzSb4ypuI5Yai+6i/EZks5IwWSIbFxECR7r47h/fu8fL8J2Z78x8d2dmP8/HOffMvZ/P597v57P3u6+5c7/3e2+qCklSO35gqTsgSTq0DH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfh50kRyb5WJL7k1SS86bVJ8k7knyrm96ZJCNu+23dNv9Vr+yIrmzdPPv5siQ3Jnkkyf1D6td19Y8n+VqSl0+r/5Ukf9OtvyXJUfN5fbXL4Nfh6gvAzwF/M6RuI/Aa4EzgBcBPA/92HtveB1yZZNUi+/h3wBbgN2aovw74KvBs4K3Ax5JMACS5ANgE/CSwDjgN+K1F9keNMPi1rHRH6b+e5PbuSPajSY6ezzaq6smq+t2q+gLw1JAmFwPvqqqpqvpr4F3AG+fxEp8BnmTwi2XBqurLVfUhYPf0uiT/CDgbuKKqnqiqjwN3AD/TNbkYeH9V3VVVDwNvZ35jUMMMfi1H/xpYD5zK4Ij8jQBJ1ib59izTz464/ecDt/WWb+vKRlXAfwauSLJ6emWSTbP1cx593F1Vj87Qz2FjeG6SZ89jHGqUwa/l6D1Vtaeq9gF/CpwFUFV/VVXHzjJ9ZMTtHwM80lt+BDhm1PP8XV+2AnuBNw2pu2q2fi6wjwf6+axZxkCvXpqRwa/lqH9e/nEGITdOjwE/1Fv+IeCxmv8dC/8Tg3Pv8zoVNaLpfaRbfnSG+gPzjyLNweDXitGd6nlslun1I27qLgYf7B5wZlc2L1W1HdgF/OK0fv7mbP2cRx9PS9I/gu/3c9gYvllV35rvONQeg18rRneq55hZpg8faJvkqN6HwkcmObp3KueDwK8mOTHJPwR+Dbi2t+79Sd44YrfeCvyHaf38r7P1s/c6P9D1cfVgMUcnObLbxteBWxl8jnB0ktcy+Lzj470xXJLkjCTHMfjr41qkERj8OlzdCzwBnAjc0M2f0tW9j8FnB3cAdwKf6srogvfZwM2jvEhV/Tnw5QX28Se6fm0D1nbzn+3VbwAmgYeBq4DXVdXe7nU/A7wTuBH4y266YoH9UGPig1ik70vyEuAtVXXRUvdFOlgMfklqjKd6JKkxBr8kNcbgl6TGHLHUHRhmzZo1tW7duqXuhiStGLfccsvfVtXEKG2XZfCvW7eOnTt3LnU3JGnFSPKXo7b1VI8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqzJzBn+TkJDcmuSfJXUl+aUibJHlPkl3ds1LP7tVdnOS+brp43AOQJM3PKNfx7wd+raq+0j0U4pYk26vq7l6bVwKnd9OLgPcCL0pyPINbxU4yeE7pLUm2dg+HliQtgTmP+Kvqwar6Sjf/KHAPg3uc910IfLAGbgaOTXICcAGwvar2dWG/ncFDtCVJS2Re39xNsg54IfClaVUnAg/0lqe6spnKh217I7ARYO3atfPp1tOs2/SpBa+7Ut1/1U8tyesu5b/1Uo1ZbViq9/ahel+P/OFukmMYPPbtl6vqO9Orh6xSs5Q/s7Bqc1VNVtXkxMRIt5uQJC3ASMGfZDWD0P9wVX1iSJMp4OTe8knAnlnKJUlLZJSregK8H7inqn5nhmZbgTd0V/e8GHikqh5k8KzT85Mc1z0Q+vyuTJK0REY5x38u8PPAHUlu7cp+k8HDoamqaxg8LPpVwC7gceAXurp9Sd4O7OjWu7Kq9o2v+5Kk+Zoz+KvqCww/V99vU8BbZqjbAmxZUO8kSWPnN3clqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY2Z80EsSbYAPw08VFU/NqT+N4DX97b3PGCie/rW/cCjwFPA/qqaHFfHJUkLM8oR/7XA+pkqq+q3q+qsqjoLuAz439Mer/iyrt7Ql6RlYM7gr6qbgFGfk3sRcN2ieiRJOqjGdo4/yd9n8JfBx3vFBXw2yS1JNo7rtSRJCzfnOf55+JfAn087zXNuVe1J8hxge5KvdX9BPEP3i2EjwNq1a8fYLUlS3ziv6tnAtNM8VbWn+/kQ8EngnJlWrqrNVTVZVZMTExNj7JYkqW8swZ/kh4GXAn/SK/vBJM86MA+cD9w5jteTJC3cKJdzXgecB6xJMgVcAawGqKprumavBT5bVX/XW/W5wCeTHHidj1TVZ8bXdUnSQswZ/FV10QhtrmVw2We/bDdw5kI7Jkk6OPzmriQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDVmzuBPsiXJQ0mGPi83yXlJHklyazdd3qtbn+TeJLuSbBpnxyVJCzPKEf+1wPo52vyfqjqrm64ESLIKuBp4JXAGcFGSMxbTWUnS4s0Z/FV1E7BvAds+B9hVVbur6kngeuDCBWxHkjRG4zrH/8+T3Jbk00me35WdCDzQazPVlQ2VZGOSnUl27t27d0zdkiRNN47g/wpwSlWdCfw34H925RnStmbaSFVtrqrJqpqcmJgYQ7ckScMsOvir6jtV9Vg3vw1YnWQNgyP8k3tNTwL2LPb1JEmLs+jgT/IPkqSbP6fb5reAHcDpSU5NciSwAdi62NeTJC3OEXM1SHIdcB6wJskUcAWwGqCqrgFeB7w5yX7gCWBDVRWwP8mlwA3AKmBLVd11UEYhSRrZnMFfVRfNUf97wO/NULcN2LawrkmSDga/uStJjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNmTP4k2xJ8lCSO2eof32S27vpi0nO7NXdn+SOJLcm2TnOjkuSFmaUI/5rgfWz1H8DeGlVvQB4O7B5Wv3LquqsqppcWBclSeM0yjN3b0qybpb6L/YWbwZOWny3JEkHy7jP8V8CfLq3XMBnk9ySZONsKybZmGRnkp179+4dc7ckSQfMecQ/qiQvYxD8L+kVn1tVe5I8B9ie5GtVddOw9atqM91posnJyRpXvyRJTzeWI/4kLwD+O3BhVX3rQHlV7el+PgR8EjhnHK8nSVq4RQd/krXAJ4Cfr6qv98p/MMmzDswD5wNDrwySJB06c57qSXIdcB6wJskUcAWwGqCqrgEuB54N/H4SgP3dFTzPBT7ZlR0BfKSqPnMQxiBJmodRruq5aI76NwFvGlK+GzjzmWtIkpaS39yVpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4JekxowU/Em2JHkoydBn5mbgPUl2Jbk9ydm9uouT3NdNF4+r45KkhRn1iP9aYP0s9a8ETu+mjcB7AZIcz+AZvS8CzgGuSHLcQjsrSVq8kYK/qm4C9s3S5ELggzVwM3BskhOAC4DtVbWvqh4GtjP7LxBJ0kE258PWR3Qi8EBveaorm6n8GZJsZPDXAmvXrh1Tt3S4WrfpU0vyuvdf9VNL8rpLNV5YujHr4BnXh7sZUlazlD+zsGpzVU1W1eTExMSYuiVJmm5cwT8FnNxbPgnYM0u5JGmJjCv4twJv6K7ueTHwSFU9CNwAnJ/kuO5D3fO7MknSEhnpHH+S64DzgDVJphhcqbMaoKquAbYBrwJ2AY8Dv9DV7UvydmBHt6krq2q2D4klSQfZSMFfVRfNUV/AW2ao2wJsmX/XJEkHg9/claTGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMaMFPxJ1ie5N8muJJuG1L87ya3d9PUk3+7VPdWr2zrOzkuS5m/ORy8mWQVcDbwCmAJ2JNlaVXcfaFNVv9Jr/++AF/Y28URVnTW+LkuSFmOUI/5zgF1VtbuqngSuBy6cpf1FwHXj6JwkafxGCf4TgQd6y1Nd2TMkOQU4Ffhcr/joJDuT3JzkNTO9SJKNXbude/fuHaFbkqSFGCX4M6SsZmi7AfhYVT3VK1tbVZPAzwK/m+RHhq1YVZurarKqJicmJkboliRpIUYJ/ing5N7yScCeGdpuYNppnqra0/3cDXyep5//lyQdYqME/w7g9CSnJjmSQbg/4+qcJP8YOA74v72y45Ic1c2vAc4F7p6+riTp0Jnzqp6q2p/kUuAGYBWwparuSnIlsLOqDvwSuAi4vqr6p4GeB7wvyfcY/JK5qn81kCTp0Jsz+AGqahuwbVrZ5dOW3zZkvS8CP76I/kmSxsxv7kpSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjRgr+JOuT3JtkV5JNQ+rfmGRvklu76U29uouT3NdNF4+z85Kk+Zvz0YtJVgFXA68ApoAdSbYOeXbuR6vq0mnrHg9cAUwCBdzSrfvwWHovSZq3UY74zwF2VdXuqnoSuB64cMTtXwBsr6p9XdhvB9YvrKuSpHEYJfhPBB7oLU91ZdP9TJLbk3wsycnzXJckG5PsTLJz7969I3RLkrQQowR/hpTVtOU/BdZV1QuA/wV8YB7rDgqrNlfVZFVNTkxMjNAtSdJCjBL8U8DJveWTgD39BlX1rar6brf4B8A/HXVdSdKhNUrw7wBOT3JqkiOBDcDWfoMkJ/QWXw3c083fAJyf5LgkxwHnd2WSpCUy51U9VbU/yaUMAnsVsKWq7kpyJbCzqrYC/z7Jq4H9wD7gjd26+5K8ncEvD4Arq2rfQRiHJGlEcwY/QFVtA7ZNK7u8N38ZcNkM624Btiyij5KkMfKbu5LUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktSYkYI/yfok9ybZlWTTkPpfTXJ3ktuT/FmSU3p1TyW5tZu2Tl9XknRozfnoxSSrgKuBVwBTwI4kW6vq7l6zrwKTVfV4kjcD7wT+TVf3RFWdNeZ+S5IWaJQj/nOAXVW1u6qeBK4HLuw3qKobq+rxbvFm4KTxdlOSNC6jBP+JwAO95amubCaXAJ/uLR+dZGeSm5O8ZqaVkmzs2u3cu3fvCN2SJC3EnKd6gAwpq6ENk58DJoGX9orXVtWeJKcBn0tyR1X9xTM2WLUZ2AwwOTk5dPuSpMUb5Yh/Cji5t3wSsGd6oyQvB94KvLqqvnugvKr2dD93A58HXriI/kqSFmmU4N8BnJ7k1CRHAhuAp12dk+SFwPsYhP5DvfLjkhzVza8BzgX6HwpLkg6xOU/1VNX+JJcCNwCrgC1VdVeSK4GdVbUV+G3gGOCPkwD8VVW9Gnge8L4k32PwS+aqaVcDSZIOsVHO8VNV24Bt08ou782/fIb1vgj8+GI6KEkaL7+5K0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0ZKfiTrE9yb5JdSTYNqT8qyUe7+i8lWderu6wrvzfJBePruiRpIeYM/iSrgKuBVwJnABclOWNas0uAh6vqR4F3A+/o1j2DwcPZnw+sB36/254kaYmMcsR/DrCrqnZX1ZPA9cCF09pcCHygm/8Y8JMZPHX9QuD6qvpuVX0D2NVtT5K0REZ52PqJwAO95SngRTO1qar9SR4Bnt2V3zxt3ROHvUiSjcDGbvGxJPeO0LeDYQ3wt0v02guSd8yr+Yob3zwc9LHN8996nJZsvx2iMfu+ZNH/1qeM2nCU4M+QshqxzSjrDgqrNgObR+jPQZVkZ1VNLnU/DpbDeXyObeU6nMe3HMc2yqmeKeDk3vJJwJ6Z2iQ5AvhhYN+I60qSDqFRgn8HcHqSU5McyeDD2q3T2mwFLu7mXwd8rqqqK9/QXfVzKnA68OXxdF2StBBznurpztlfCtwArAK2VNVdSa4EdlbVVuD9wIeS7GJwpL+hW/euJH8E3A3sB95SVU8dpLGMy5KfbjrIDufxObaV63Ae37IbWwYH5pKkVvjNXUlqjMEvSY1pLviTbEnyUJI7e2XHJ9me5L7u53FdeZK8p7vlxO1Jzl66ns9thrG9LclfJ7m1m17Vq1sxt9NIcnKSG5Pck+SuJL/UlR8u+26m8a34/Zfk6CRfTnJbN7bf6spP7W7xcl93y5cju/IZbwGz3MwytmuTfKO3387qypfH+7KqmpqAnwDOBu7slb0T2NTNbwLe0c2/Cvg0g+8jvBj40lL3fwFjexvw60PangHcBhwFnAr8BbBqqccwy9hOAM7u5p8FfL0bw+Gy72Ya34rff90+OKabXw18qdsnfwRs6MqvAd7czf8icE03vwH46FKPYQFjuxZ43ZD2y+J92dwRf1XdxODKo77+LSc+ALymV/7BGrgZODbJCYemp/M3w9hmsqJup1FVD1bVV7r5R4F7GHwL/HDZdzONbyYrZv91++CxbnF1NxXwLxjc4gWeue+G3QJm2ZllbDNZFu/L5oJ/Bs+tqgdh8B8QeE5XPux2FbP9Z1yuLu3+rNxy4FQIK3hs3Z/+L2RwdHXY7btp44PDYP8lWZXkVuAhYDuDv1C+XVX7uyb9/j/tFjDAgVvALEvTx1ZVB/bbf+n227uTHNWVLYv9ZvDPbuRbTixj7wV+BDgLeBB4V1e+IseW5Bjg48AvV9V3Zms6pGwlju+w2H9V9VRVncXg2/vnAM8b1qz7uaLHluTHgMuAfwL8M+B44D92zZfF2Az+gW8e+HOr+/lQV77ibzlRVd/s3pjfA/6A758OWHFjS7KaQSh+uKo+0RUfNvtu2PgOp/0HUFXfBj7P4Pz2sRnc4gWe3v+ZbgGzrPXGtr47dVdV9V3gD1lm+83gH+jfcuJi4E965W/oPol/MfDIgdMKK8W084evBQ5c8bOibqfRneN9P3BPVf1Or+qw2Hczje9w2H9JJpIc283/PeDlDD7DuJHBLV7gmftu2C1glp0Zxva13sFIGHx20d9vS/++XIpPlJdyAq5j8Cfz/2Pw2/cSBucP/wy4r/t5fH3/E/urGZyPvAOYXOr+L2BsH+r6fjuDN90JvfZv7cZ2L/DKpe7/HGN7CYM/iW8Hbu2mVx1G+26m8a34/Qe8APhqN4Y7gcu78tMY/LLaBfwxcFRXfnS3vKurP22px7CAsX2u2293Av+D71/5syzel96yQZIa46keSWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5Ia8/8BQVGGrB2Rr8IAAAAASUVORK5CYII=\n", "text/plain": [ "<matplotlib.figure.Figure at 0x7ff3e1b8e320>" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAEGhJREFUeJzt3X+sZHV9xvH3IysiIF2QC8Fd6mKy8Ucb+dENojSGgm0AjfCHNljbLpZ0mxQrVhtd1Naa9A8wVdHGaDeiYmsQi7SLaFSyQhpNunqX38tKoUhhBdmrFRBtReTTP+bcOi6X3XvnzO3Mfn2/ksnM+Z4z5zxz5+5zz37nztxUFZKkdj1t0gEkScvLopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mNs+jVnCT7J7kyyT1JKskpu61PkouTfL+7vDdJhtYfl2Rbkh9318ct4diV5NYkTxsa+5sknxzhcWxKckeSJ5Kcu9T7S/MserXqa8DvA99dYN0G4GzgWODFwKuAP4HBDwlgM/CPwKHAZcDmbnyxngOcM3Lyn7sZ+FPghjHsS7/ELHpNle4s/C+S3JLk4SRXJDlgKfuoqseq6pKq+hrwswU2WQ+8r6p2VtV3gPcB53brTgFWAJdU1U+q6kNAgFOXEOG9wHuSrFhK7t1V1YeragvwP332I1n0mka/C5wOHMPgjPtcgCS/muShPVx+b5H7/zUGZ8vzbu7G5tfdUr/42SC3DK1fjKuAR/j5D49fsJfHsHEJx5EWpdcZh7RMPlRV9wMk+TxwHEBV3QusHMP+DwYeHlp+GDi4m6fffd38+mctYf8F/CXw0ST/8KSVVeN4DNKieUavaTQ8r/5jBuU7To8ChwwtHwI82p3F775ufv0Pl3KAqvoicC+D1wOkibLotc/opm4e3cPl9Yvc1XYGL8TOO7Ybm1/34uHfwmEwfbSdpXsX8E7gwN0ex54ewztGOI60R07daJ/RTd0s6uw+yTMYvIgKsH/3gu5PurP2TwFvSfJFBtMsbwX+rtv2egYv4L4pyUeBP+7Gv9rt91zgr6tqzSLyXp/kVgYv/n5+aHyxj2F/BidjAZ7ePYbHquqJxdxfmucZvVp1B/DfwCrgy93t53br/p5B8d4K3AZ8oRujqh5j8KuXfwg8BPwRcHY3DnA08PUl5HgXcNiIj+ErXe6XAZu62y8fcV/6JRb/8Ii0eEm+AlxQVTsmnUVaLItekhrn1I0kNc6il6TGWfSS1Lip+PXKww8/vNasWTPpGJK0T9m2bdv3qmpmb9tNRdGvWbOG2dnZSceQpH1Kkv9czHZO3UhS4yx6SWqcRS9Jjdtr0Sf5eJJdSW4bGjssybVJ7uyuD+3Gk+RDSe7q/nDECcsZXpK0d4s5o/8kgz8CMWwjsKWq1gJbumWAM4C13WUD8JHxxJQkjWqvRV9V/wr8127DZzH4W5p012cPjX+qBv4NWJnkqHGFlSQt3ahz9EdW1QMA3fUR3fgq4L6h7XZ2Y5KkCRn3i7FZYGzBT01LsiHJbJLZubm5MceQJM0btegfnJ+S6a53deM7GXxe97zVwP0L7aCqNlXVuqpaNzOz1zd2SZJGNOo7Y69m8FdzLuquNw+NvzHJZ4CXAA/PT/EslzUbv7Ccu9+jey565cSOLUmLtdeiT3I5cApweJKdwLsZFPxnk5zH4A8gv7bb/IvAmcBdDP6o8xuWIbMkaQn2WvRV9bqnWHXaAtsWcH7fUJKk8fGdsZLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcRa9JDXOopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mNs+glqXEWvSQ1zqKXpMZZ9JLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1LheRZ/kz5NsT3JbksuTHJDkmCRbk9yZ5Iok+48rrCRp6UYu+iSrgDcB66rq14H9gHOAi4EPVNVa4AfAeeMIKkkaTd+pmxXAM5OsAA4EHgBOBa7s1l8GnN3zGJKkHkYu+qr6DvC3wL0MCv5hYBvwUFU93m22E1i10P2TbEgym2R2bm5u1BiSpL3oM3VzKHAWcAzwHOAg4IwFNq2F7l9Vm6pqXVWtm5mZGTWGJGkv+kzdvAL4dlXNVdVPgauAlwEru6kcgNXA/T0zSpJ66FP09wInJTkwSYDTgNuB64DXdNusBzb3iyhJ6qPPHP1WBi+63gDc2u1rE/B24C1J7gKeDVw6hpySpBGt2PsmT62q3g28e7fhu4ET++xXkjQ+vjNWkhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcRa9JDXOopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mNs+glqXEWvSQ1zqKXpMZZ9JLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1DiLXpIa16vok6xMcmWSbyXZkeSlSQ5Lcm2SO7vrQ8cVVpK0dH3P6D8IfKmqXgAcC+wANgJbqmotsKVbliRNyMhFn+QQ4OXApQBV9VhVPQScBVzWbXYZcHbfkJKk0fU5o38eMAd8IsmNST6W5CDgyKp6AKC7PmKhOyfZkGQ2yezc3FyPGJKkPelT9CuAE4CPVNXxwI9YwjRNVW2qqnVVtW5mZqZHDEnSnvQp+p3Azqra2i1fyaD4H0xyFEB3vatfRElSHyMXfVV9F7gvyfO7odOA24GrgfXd2Hpgc6+EkqReVvS8/58Bn06yP3A38AYGPzw+m+Q84F7gtT2PIUnqoVfRV9VNwLoFVp3WZ7+SpPHxnbGS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcRa9JDXOopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mNs+glqXEWvSQ1zqKXpMZZ9JLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcRa9JDXOopekxln0ktS43kWfZL8kNya5pls+JsnWJHcmuSLJ/v1jSpJGNY4z+guAHUPLFwMfqKq1wA+A88ZwDEnSiHoVfZLVwCuBj3XLAU4Fruw2uQw4u88xJEn99D2jvwR4G/BEt/xs4KGqerxb3gmsWuiOSTYkmU0yOzc31zOGJOmpjFz0SV4F7KqqbcPDC2xaC92/qjZV1bqqWjczMzNqDEnSXqzocd+TgVcnORM4ADiEwRn+yiQrurP61cD9/WNKkkY18hl9VV1YVaurag1wDvDVqno9cB3wmm6z9cDm3iklSSNbjt+jfzvwliR3MZizv3QZjiFJWqQ+Uzf/p6quB67vbt8NnDiO/UqS+vOdsZLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcRa9JDXOopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mNs+glqXEWvSQ1zqKXpMZZ9JLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1LiRiz7J0UmuS7IjyfYkF3TjhyW5Nsmd3fWh44srSVqqPmf0jwNvraoXAicB5yd5EbAR2FJVa4Et3bIkaUJGLvqqeqCqbuhu/xDYAawCzgIu6za7DDi7b0hJ0ujGMkefZA1wPLAVOLKqHoDBDwPgiKe4z4Yks0lm5+bmxhFDkrSA3kWf5GDgc8Cbq+qRxd6vqjZV1bqqWjczM9M3hiTpKfQq+iRPZ1Dyn66qq7rhB5Mc1a0/CtjVL6IkqY8+v3UT4FJgR1W9f2jV1cD67vZ6YPPo8SRJfa3ocd+TgT8Abk1yUzf2DuAi4LNJzgPuBV7bL6IkqY+Ri76qvgbkKVafNup+JUnj5TtjJalxFr0kNc6il6TGWfSS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcRa9JDXOopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mNs+glqXEWvSQ1zqKXpMZZ9JLUOItekhpn0UtS4yx6SWqcRS9JjbPoJalxFr0kNc6il6TGWfSS1DiLXpIaZ9FLUuMseklqnEUvSY2z6CWpcctS9ElOT3JHkruSbFyOY0iSFmfFuHeYZD/gw8BvAzuBbya5uqpuH/exJGkc1mz8wsSOfc9Fr1z2YyzHGf2JwF1VdXdVPQZ8BjhrGY4jSVqEsZ/RA6uA+4aWdwIv2X2jJBuADd3io0nuGPF4hwPfG/G+veTiRW02sXyLZL5+pj0fTH/GX+p8i+yRp/LcxWy0HEWfBcbqSQNVm4BNvQ+WzFbVur77WS7m68d8/U17RvMtv+WYutkJHD20vBq4fxmOI0lahOUo+m8Ca5Mck2R/4Bzg6mU4jiRpEcY+dVNVjyd5I/BlYD/g41W1fdzHGdJ7+meZma8f8/U37RnNt8xS9aTpc0lSQ3xnrCQ1zqKXpMbt00U/bR+1kOTjSXYluW1o7LAk1ya5s7s+dIL5jk5yXZIdSbYnuWCaMiY5IMk3ktzc5XtPN35Mkq1dviu6F/knJsl+SW5Mcs205UtyT5Jbk9yUZLYbm4rnt8uyMsmVSb7VfR++dFryJXl+93WbvzyS5M3Tkq+Pfbbohz5q4QzgRcDrkrxosqn4JHD6bmMbgS1VtRbY0i1PyuPAW6vqhcBJwPnd12xaMv4EOLWqjgWOA05PchJwMfCBLt8PgPMmlG/eBcCOoeVpy/dbVXXc0O9+T8vzC/BB4EtV9QLgWAZfx6nIV1V3dF+344DfAH4M/PO05OulqvbJC/BS4MtDyxcCF05BrjXAbUPLdwBHdbePAu6YdMahbJsZfCbR1GUEDgRuYPCu6u8BKxZ63ieQazWDf+ynAtcweIPgNOW7Bzh8t7GpeH6BQ4Bv0/0SyLTl2y3T7wBfn9Z8S73ss2f0LPxRC6smlGVPjqyqBwC66yMmnAeAJGuA44GtTFHGblrkJmAXcC3wH8BDVfV4t8mkn+dLgLcBT3TLz2a68hXwlSTbuo8Zgel5fp8HzAGf6Ka+PpbkoCnKN+wc4PLu9jTmW5J9uegX9VELerIkBwOfA95cVY9MOs+wqvpZDf7rvJrBB+S9cKHN/n9TDSR5FbCrqrYNDy+w6SS/D0+uqhMYTGmen+TlE8yyuxXACcBHqup44EdM4TRI9xrLq4F/mnSWcdmXi35f+aiFB5McBdBd75pkmCRPZ1Dyn66qq7rhqcoIUFUPAdczeC1hZZL5N/dN8nk+GXh1knsYfCrrqQzO8KclH1V1f3e9i8H88olMz/O7E9hZVVu75SsZFP+05Jt3BnBDVT3YLU9bviXbl4t+X/mohauB9d3t9QzmxSciSYBLgR1V9f6hVVORMclMkpXd7WcCr2DwYt11wGsmna+qLqyq1VW1hsH321er6vXTki/JQUmeNX+bwTzzbUzJ81tV3wXuS/L8bug04HamJN+Q1/HzaRuYvnxLN+kXCXq+YHIm8O8M5nHfOQV5LgceAH7K4OzlPAZzuFuAO7vrwyaY7zcZTCvcAtzUXc6clozAi4Ebu3y3AX/VjT8P+AZwF4P/Tj9jCp7rU4Brpilfl+Pm7rJ9/t/EtDy/XZbjgNnuOf4X4NApy3cg8H3gV4bGpibfqBc/AkGSGrcvT91IkhbBopekxln0ktQ4i16SGmfRS1LjLHpJapxFL0mN+19hMsOYH4DyOgAAAABJRU5ErkJggg==\n", "text/plain": [ "<matplotlib.figure.Figure at 0x7ff3e1b87b70>" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEICAYAAABPgw/pAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAE6NJREFUeJzt3X+QXedd3/H3BytKcBLjX2uPsGMkM2pK2sGKu5O6pGSolYDthEhtk4zTtGzAg9oplITANILQJpnyh50SEsJkCCI2FmASB2OPBKHBHmGXSacYVo5/KYor2zhGsSItxj/imsY4+faPe5as17u6Z3/p7j5+v2bunHOec+4933vu7meffe6596SqkCStfd826gIkScvDQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNC1ZiRZn+SGJA8lqSQ/MGt9klyV5NHu9qEkmbF+S5L9SZ7uplsWsO9Kck+Sb5vR9otJrl3E89iV5L4k30zyzjnW/3SSryZ5Isk1SV680H3ohclA11rzeeDfAl+dY90OYDtwAfC9wJuAfw+DPwbAHuB3gNOA3cCerr2v7wQuX3Tl33IX8B+BO2avSPJDwE5gK7AROB/44DLsUy8ABrpOiK5X/bNJ7u56ntcneclCHqOqnqmqj1bV54FvzLHJBPDhqjpcVV8BPgy8s1v3A8A64KNV9fWq+hgQ4OIFlPAh4INJ1i2k7tmq6uNVtQ/4f3OsngCurqoDVfUY8N/41nOQjstA14n0NuASYBODHvQ7AZKcl+Tx49z+Tc/H/0cMer/T7uraptfdXc/9rou7Z6zv40bgSeYJ2CHPYecSnsPZSc5YQJ16gVpST0NaoI9V1SMASf4A2AJQVQ8Dpy7D478MeGLG8hPAy7px9Nnrpte/fAGPX8B/AT6R5Left7JqpZ4DDOp8dBkeXw2zh64Taea499MMwms5PQWcMmP5FOCprlc+e930+q8tZAdV9UfAwwzG61fCXM8BFlinXpgMdI1cN+Ty1HFu7+j5UAcYvCE67YKubXrd984864XBsM8BFu4XgPcBJ896Hsd7Dj+/hOdwtKrsnWsoh1w0ct2QS6/eencK33Qor+/eWP161wv/LeA9Sf6IwfDIzwC/2m17G4M3Un8qySeAH+/a/6R73HcCH6iqjT3qvS3JPQzewPyDGe19n8N6Bp2pAC/qnsMzVfXN7jlcm+Q64AiDPx7X9nlcyR661pr7gL8FzgH+uJv/rm7drzMI2HuAe4HPdm1U1TMMTmn8EeBx4MeA7V07wCuA/7WAOn4BOH2Rz+Hmru7vA3Z186/r6vwcg7NpbgW+3N3ev8j96AUmXuBCgiQ3A++qqoOjrkVaLANdkhrhkIskNcJAl6RGGOiS1IgTetrimWeeWRs3bjyRu5SkNW///v1/XVVjw7Y7oYG+ceNGJicnT+QuJWnNS/LlPts55CJJjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY1YM1cs2rjzsyPZ70NXvnEk+5WkhbKHLkmN6BXoSX46yYEk9yb5VJKXJNmU5PYkh5Jc310nUZI0IkMDPck5wE8B41X1j4GTgMuBq4CPVNVm4DHgipUsVJJ0fH2HXNYB355kHXAyg6uRXwzc0K3fzeACvJKkERka6FX1FeCXgIcZBPkTwH7g8ap6ttvsMIOrsD9Pkh1JJpNMTk1NLU/VkqTn6TPkchqwDdgEfCfwUuDSOTad82rTVbWrqsaranxsbOj3s0uSFqnPkMvrgb+sqqmq+jvgRuD7gFO7IRiAc4FHVqhGSVIPfQL9YeCiJCcnCbAV+CJwK/CWbpsJYM/KlChJ6qPPGPrtDN78vAO4p7vPLuC9wHuS3A+cAVy9gnVKkobo9UnRqno/8P5ZzQ8Cr1n2iiRJi+InRSWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjehzkehXJrlzxu3JJO9OcnqSW5Ic6qannYiCJUlz63MJuvuqaktVbQH+CfA0cBOwE9hXVZuBfd2yJGlEFjrkshV4oKq+DGwDdnftu4Hty1mYJGlhFhrolwOf6ubPrqojAN30rLnukGRHkskkk1NTU4uvVJJ0XL0DPcl64M3A7y1kB1W1q6rGq2p8bGxsofVJknpaSA/9UuCOqjraLR9NsgGgmx5b7uIkSf0tJNDfzreGWwD2AhPd/ASwZ7mKkiQtXK9AT3Iy8AbgxhnNVwJvSHKoW3fl8pcnSeprXZ+Nqupp4IxZbY8yOOtFkrQK+ElRSWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1Ij+l6x6NQkNyT5UpKDSf5ZktOT3JLkUDc9baWLlSTNr28P/VeAz1XVPwQuAA4CO4F9VbUZ2NctS5JGZGigJzkFeB1wNUBVPVNVjwPbgN3dZruB7StVpCRpuD499POBKeA3k3whySeTvBQ4u6qOAHTTs+a6c5IdSSaTTE5NTS1b4ZKk5+oT6OuAC4Ffq6pXA/+XBQyvVNWuqhqvqvGxsbFFlilJGqZPoB8GDlfV7d3yDQwC/miSDQDd9NjKlChJ6mNooFfVV4G/SvLKrmkr8EVgLzDRtU0Ae1akQklSL+t6bvefgOuSrAceBH6UwR+DzyS5AngYeOvKlChJ6qNXoFfVncD4HKu2Lm85kqTF8pOiktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RG9LrARZKHgK8B3wCerarxJKcD1wMbgYeAt1XVYytTpiRpmIX00P9FVW2pqukrF+0E9lXVZmBftyxJGpGlDLlsA3Z387uB7UsvR5K0WH0DvYCbk+xPsqNrO7uqjgB007NWokBJUj+9xtCB11bVI0nOAm5J8qW+O+j+AOwAOO+88xZRoiSpj1499Kp6pJseA24CXgMcTbIBoJsem+e+u6pqvKrGx8bGlqdqSdLzDA30JC9N8vLpeeAHgXuBvcBEt9kEsGelipQkDddnyOVs4KYk09v/blV9LslfAJ9JcgXwMPDWlStTkjTM0ECvqgeBC+ZofxTYuhJFSZIWzk+KSlIjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiN6B3qSk5J8Ickfdsubktye5FCS65OsX7kyJUnDLKSH/i7g4Izlq4CPVNVm4DHgiuUsTJK0ML0CPcm5wBuBT3bLAS4Gbug22Q1sX4kCJUn99O2hfxT4z8A3u+UzgMer6tlu+TBwzlx3TLIjyWSSyampqSUVK0ma39BAT/Im4FhV7Z/ZPMemNdf9q2pXVY1X1fjY2Ngiy5QkDbOuxzavBd6c5DLgJcApDHrspyZZ1/XSzwUeWbkyJUnDDO2hV9XPVdW5VbURuBz4k6p6B3Ar8JZuswlgz4pVKUkaainnob8XeE+S+xmMqV+9PCVJkhajz5DL36uq24DbuvkHgdcsf0mSpMXwk6KS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEb0uUj0S5L8eZK7khxI8sGufVOS25McSnJ9kvUrX64kaT59euhfBy6uqguALcAlSS4CrgI+UlWbgceAK1auTEnSMH0uEl1V9VS3+KLuVsDFwA1d+25g+4pUKEnqpdcYepKTktwJHANuAR4AHq+qZ7tNDgPnzHPfHUkmk0xOTU0tR82SpDn0CvSq+kZVbQHOZXBh6O+Za7N57rurqsaranxsbGzxlUqSjmtBZ7lU1ePAbcBFwKlJ1nWrzgUeWd7SJEkL0ecsl7Ekp3bz3w68HjgI3Aq8pdtsAtizUkVKkoZbN3wTNgC7k5zE4A/AZ6rqD5N8Efh0kl8EvgBcvYJ1SpKGGBroVXU38Oo52h9kMJ4uSVoF/KSoJDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhphoEtSIwx0SWqEgS5JjTDQJakRfS5B94oktyY5mORAknd17acnuSXJoW562sqXK0maT58e+rPAz1TV9zC4OPRPJHkVsBPYV1WbgX3dsiRpRIYGelUdqao7uvmvMbhA9DnANmB3t9luYPtKFSlJGm5BY+hJNjK4vujtwNlVdQQGoQ+cNc99diSZTDI5NTW1tGolSfPqHehJXgb8PvDuqnqy7/2qaldVjVfV+NjY2GJqlCT10CvQk7yIQZhfV1U3ds1Hk2zo1m8Ajq1MiZKkPvqc5RLgauBgVf3yjFV7gYlufgLYs/zlSZL6Wtdjm9cC/w64J8mdXdvPA1cCn0lyBfAw8NaVKVGS1MfQQK+qzwOZZ/XW5S1HkrRYflJUkhphoEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMMdElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktSIPpeguybJsST3zmg7PcktSQ5109NWtkxJ0jB9eujXApfMatsJ7KuqzcC+blmSNEJDA72q/hT4m1nN24Dd3fxuYPsy1yVJWqDFjqGfXVVHALrpWctXkiRpMVb8TdEkO5JMJpmcmppa6d1J0gvWYgP9aJINAN302HwbVtWuqhqvqvGxsbFF7k6SNMxiA30vMNHNTwB7lqccSdJi9Tlt8VPA/wZemeRwkiuAK4E3JDkEvKFbliSN0LphG1TV2+dZtXWZa5EkLYGfFJWkRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqRFDz0PX6Gzc+dmR7PehK984kv1KWhp76JLUCANdkhphoEtSIwx0SWqEb4oOMao3JiUtv9ZPNLCHLkmNMNAlqREOueh5XojDTJ57rxbYQ5ekRiyph57kEuBXgJOAT1aVl6LTmvRC/K9klPyPaGUsuoee5CTg48ClwKuAtyd51XIVJklamKUMubwGuL+qHqyqZ4BPA9uWpyxJ0kItZcjlHOCvZiwfBv7p7I2S7AB2dItPJblvkfs7E/jrRd73RFoLda6FGmFt1LkWaoRVVmeumnfVqqpzHguu8TjPt6/v6rPRUgI9c7TV8xqqdgG7lrCfwc6SyaoaX+rjrLS1UOdaqBHWRp1roUawzuW0mmtcypDLYeAVM5bPBR5ZWjmSpMVaSqD/BbA5yaYk64HLgb3LU5YkaaEWPeRSVc8m+UngjxmctnhNVR1Ytsqeb8nDNifIWqhzLdQIa6POtVAjWOdyWrU1pup5w96SpDXIT4pKUiMMdElqxJoI9CSXJLkvyf1Jdo66HoAkr0hya5KDSQ4keVfX/oEkX0lyZ3e7bBXU+lCSe7p6Jru205PckuRQNz1thPW9csbxujPJk0nevRqOZZJrkhxLcu+MtjmPXQY+1v2c3p3kwhHX+d+TfKmr5aYkp3btG5P87Yzj+okR1jjva5zk57pjeV+SHzoRNR6nzutn1PhQkju79pEcy3lV1aq+MXjD9QHgfGA9cBfwqlVQ1wbgwm7+5cD/YfAVCB8AfnbU9c2q9SHgzFltHwJ2dvM7gatGXeeM1/urDD5IMfJjCbwOuBC4d9ixAy4D/geDz2hcBNw+4jp/EFjXzV81o86NM7cbcY1zvsbd79JdwIuBTV0GnDSqOmet/zDwX0d5LOe7rYUe+qr8ioGqOlJVd3TzXwMOMvj07FqxDdjdze8Gto+wlpm2Ag9U1ZdHXQhAVf0p8Dezmuc7dtuA36qBPwNOTbJhVHVW1c1V9Wy3+GcMPisyMvMcy/lsAz5dVV+vqr8E7meQBSvueHUmCfA24FMnopaFWguBPtdXDKyq4EyyEXg1cHvX9JPdv7nXjHIoY4YCbk6yv/sqBoCzq+oIDP44AWeNrLrnupzn/rKstmMJ8x+71fyz+mMM/nuYtinJF5L8zyTfP6qiOnO9xqv1WH4/cLSqDs1oWzXHci0Eeq+vGBiVJC8Dfh94d1U9Cfwa8N3AFuAIg3/PRu21VXUhg2/G/Ikkrxt1QXPpPqD2ZuD3uqbVeCyPZ1X+rCZ5H/AscF3XdAQ4r6peDbwH+N0kp4yovPle41V5LIG389wOx2o6lmsi0FftVwwkeRGDML+uqm4EqKqjVfWNqvom8BucoH8Tj6eqHummx4CbGNR0dHo4oJseG12Ff+9S4I6qOgqr81h25jt2q+5nNckE8CbgHdUN+nbDGI928/sZjE//g1HUd5zXeDUey3XAvwKun25bTccS1kagr8qvGOjG0q4GDlbVL89onzlm+i+Be2ff90RK8tIkL5+eZ/BG2b0MjuFEt9kEsGc0FT7Hc3o/q+1YzjDfsdsL/Eh3tstFwBPTQzOjkMEFaN4LvLmqnp7RPpbB9QxIcj6wGXhwRDXO9xrvBS5P8uIkmxjU+Ocnur5ZXg98qaoOTzespmMJrP6zXLpOxWUMziJ5AHjfqOvpavrnDP4FvBu4s7tdBvw2cE/XvhfYMOI6z2dwtsBdwIHp4wecAewDDnXT00dc58nAo8B3zGgb+bFk8AfmCPB3DHqNV8x37BgME3y8+zm9BxgfcZ33MxiHnv75/ES37b/ufhbuAu4AfniENc77GgPv647lfcClozyWXfu1wH+Yte1IjuV8Nz/6L0mNWAtDLpKkHgx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS1Ij/D6wnhV2fNN93AAAAAElFTkSuQmCC\n", "text/plain": [ "<matplotlib.figure.Figure at 0x7ff3e1ab5748>" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEICAYAAABYoZ8gAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAFfZJREFUeJzt3X+wZGV95/H3Jwi4EbKAXAkC40CWdYOWjNRdNMuuhaIISAXdtbLDZuO4ITXGaK3umqqMMSv+2N1SU2rKYEkmMgsYF92oRFKgMoVYhC1/DWSAYREhZNRxCDMK8qN0ddHv/tHnlk3TfedOd8/tO/O8X1Vdfc7zPH3O95y693PPPd19TqoKSVI7fmHWBUiSlpfBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8Gu/l+SQJJ9Ksj1JJTlzoD9J3pvk+93jfUnS178myS1Jftg9r9mLdVeSO5L8Ql/bf01y+RjbsTHJ3Ul+luS1Q/r/U5J/SPJwkk1JDu3rW53kxm4bvpHkpXu7frXD4NeB4mbg3wP/MKRvPfBK4FTgecD5wOug90cD+CzwF8CRwBXAZ7v2pXomsHbsyn/uNuD3gFsHO5K8HNgAnAWsBk4C3tk35Crgb4GnA28DPpVkbgo16QBk8GumuqP0309ye3ck+8kkT92bZVTVT6rqT6rqZuCnQ4asA95fVTuq6rvA+4HXdn1nAk8B/qSqflxVHwICvGQvSngf8M4kT9mbugdV1Yer6gbg/w7pXgdcVlV3VtVDwLvptiHJPwVOAy6uqh9V1aeBO4B/M0k9OnAZ/FoJfgM4BziR3hH5awGSrEryg0Ue/26Jy38OvaPpBbd1bQt9t9cTr11ye1//UnwGeISf/zF5gj1sw4YJtuGYJE/v+u6rqkcH+vdmG9SQiY5QpCn5UFXtBEjy18AagKr6NnDEFJZ/GPBw3/zDwGHdef7BvoX+w/di+QX8F+DSJB97UmfVvtoG6NU5ahuOm8J6dQDyiF8rQf95+R/SC7Jpegz4pb75XwIe647yB/sW+h9lL1TVdcC36b2fsC8M2wbo1TmVbVA7DH6tWN2pnscWefzmEhd1J703dhec2rUt9D2v/1M+9E433cne+yN6b6z+4sB2LLYNfzjBNjxQVd/v+k5KcvhA/zjboAYY/FqxqurbVXXYIo+PL4xNcmjfm8KHJHlqX5hfCfznJMcleSbwFuDyru9L9N4Q/o/dMt7YtX+xW+5rk2xfYr1fovem6rqB9sW24b/3bcMh3TYEOLjbhoXf0SuBi5KckuRIen9kLu+W/01gK3Bx95pX0fvj9eml1K32GPw6UNwN/Ijeee0vdNPP6vr+DPhreqG8Dbi2a6OqfkLvo56vAX4A/Dbwyq4d4ATgf+9FHX8EHDXmNlzf1f0vgI3d9Iu6Oj9P79NDNwLf6h4X9712LTAPPAS8B3h1Ve0esw4d4OKNWKTRklwPvKmq7pp1LdK0GPyS1BhP9UhSYwx+SWqMwS9JjVmR39w9+uija/Xq1bMuQ5L2G7fccsv3qmpJF+ZbkcG/evVqtmzZMusyJGm/keRbSx3rqR5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUmD0Gf5ITktyY5K4kdyZ5U9d+VJLNSe7pno8c8fp13Zh7kqwbNkaStHyWcsT/OPCWqvpV4IXAG5KcAmwAbqiqk4EbuvknSHIUvUvHvgA4nd71wof+gZAkLY89Bn9V3V9Vt3bTjwJ30bvm+QXAFd2wK+hd03zQy4HNVfVgVT0EbKZ3U21J0ozs1Td3k6wGng98FTimqu6H3h+HJM8Y8pLjgO/0ze9gxA2gk6ynu1/pqlWr9qYszcjqDdfObN3b3/OKma1b2t8t+c3dJIfRu5Xbm6vqkaW+bEjb0BsAVNXGqpqvqvm5uSVdbkKSNIYlBX+Sg+mF/ser6jNd8wNJju36jwV2DXnpDnq3rltwPLBz/HIlSZNayqd6AlwG3FVVH+jruoaf31R6HfDZIS//AnB2kiO7N3XP7tokSTOylCP+M4DfAl6SZGv3OI/eDZ1fluQe4GXdPEnmk3wUoKoeBN4NfL17vKtrkyTNyB7f3K2qmxl+rh7grCHjtwC/0ze/Cdg0boGSpOnym7uS1BiDX5IaY/BLUmMMfklqzIq85660Us3q28p+U1nT5BG/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUmD1eqyfJJuB8YFdVPbdr+yTw7G7IEcAPqmrNkNduBx4Ffgo8XlXzU6pbkjSmpVyk7XLgEuDKhYaq+rcL00neDzy8yOtfXFXfG7dASdJ0LeXWizclWT2sr7sR+28AL5luWZKkfWXSc/z/Cnigqu4Z0V/A9UluSbJ+wnVJkqZg0uvxXwhctUj/GVW1M8kzgM1JvlFVNw0b2P1hWA+watWqCcuSJI0y9hF/kqcA/xr45KgxVbWze94FXA2cvsjYjVU1X1Xzc3Nz45YlSdqDSU71vBT4RlXtGNaZ5GlJDl+YBs4Gtk2wPknSFOwx+JNcBXwZeHaSHUku6rrWMnCaJ8kzk1zXzR4D3JzkNuBrwLVV9fnplS5JGsdSPtVz4Yj21w5p2wmc103fB5w6YX2SpCnzm7uS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUmKXcenFTkl1JtvW1vSPJd5Ns7R7njXjtOUnuTnJvkg3TLFySNJ6lHPFfDpwzpP2DVbWme1w32JnkIODDwLnAKcCFSU6ZpFhJ0uT2GPxVdRPw4BjLPh24t6ruq6qfAJ8ALhhjOZKkKdrjzdYX8cYkrwG2AG+pqocG+o8DvtM3vwN4waiFJVkPrAdYtWrVBGWpBas3XDvrEqT91rhv7n4E+BVgDXA/8P4hYzKkrUYtsKo2VtV8Vc3Pzc2NWZYkaU/GCv6qeqCqflpVPwP+nN5pnUE7gBP65o8Hdo6zPknS9IwV/EmO7Zt9FbBtyLCvAycnOTHJIcBa4Jpx1idJmp49nuNPchVwJnB0kh3AxcCZSdbQO3WzHXhdN/aZwEer6ryqejzJG4EvAAcBm6rqzn2yFZKkJdtj8FfVhUOaLxsxdidwXt/8dcCTPuopSZodv7krSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYya52br6zPLm39vf84qZrVvS/scjfklqzB6DP8mmJLuSbOtr++Mk30hye5Krkxwx4rXbk9yRZGuSLdMsXJI0nqUc8V8OnDPQthl4blU9D/gm8NZFXv/iqlpTVfPjlShJmqY9Bn9V3QQ8ONB2fVU93s1+BTh+H9QmSdoHpnGO/7eBz43oK+D6JLckWb/YQpKsT7IlyZbdu3dPoSxJ0jATBX+StwGPAx8fMeSMqjoNOBd4Q5IXjVpWVW2sqvmqmp+bm5ukLEnSIsYO/iTrgPOB36yqGjamqnZ2z7uAq4HTx12fJGk6xgr+JOcAfwD8elX9cMSYpyU5fGEaOBvYNmysJGn5LOXjnFcBXwaenWRHkouAS4DDgc3dRzUv7cY+M8l13UuPAW5OchvwNeDaqvr8PtkKSdKS7fGbu1V14ZDmy0aM3Qmc103fB5w6UXWSpKnzm7uS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktQYg1+SGuPN1g8As7zRu6T9j0f8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTFLCv4km5LsSrKtr+2oJJuT3NM9Hzniteu6Mfd0N2iXJM3QUo/4LwfOGWjbANxQVScDN3TzT5DkKOBi4AXA6cDFo/5ASJKWx5KCv6puAh4caL4AuKKbvgJ45ZCXvhzYXFUPVtVDwGae/AdEkrSMJjnHf0xV3Q/QPT9jyJjjgO/0ze/o2p4kyfokW5Js2b179wRlSZIWs6/f3M2Qtho2sKo2VtV8Vc3Pzc3t47IkqV2TBP8DSY4F6J53DRmzAzihb/54YOcE65QkTWiS4L8GWPiUzjrgs0PGfAE4O8mR3Zu6Z3dtkqQZWerHOa8Cvgw8O8mOJBcB7wFeluQe4GXdPEnmk3wUoKoeBN4NfL17vKtrkyTNyJIuy1xVF47oOmvI2C3A7/TNbwI2jVWdJGnq/OauJDXG4Jekxhj8ktQYg1+SGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNWbs4E/y7CRb+x6PJHnzwJgzkzzcN+btk5csSZrEkm69OExV3Q2sAUhyEPBd4OohQ/+mqs4fdz2SpOma1qmes4C/q6pvTWl5kqR9ZFrBvxa4akTfryW5Lcnnkjxn1AKSrE+yJcmW3bt3T6ksSdKgiYM/ySHArwN/OaT7VuBZVXUq8KfAX41aTlVtrKr5qpqfm5ubtCxJ0gjTOOI/F7i1qh4Y7KiqR6rqsW76OuDgJEdPYZ2SpDFNI/gvZMRpniS/nCTd9Ond+r4/hXVKksY09qd6AJL8IvAy4HV9bb8LUFWXAq8GXp/kceBHwNqqqknWKUmazETBX1U/BJ4+0HZp3/QlwCWTrEOSNF1+c1eSGmPwS1JjDH5JaozBL0mNMfglqTEGvyQ1xuCXpMYY/JLUGINfkhpj8EtSYya6ZIOk5bF6w7UzW/f297xiZuvWvuERvyQ1xuCXpMYY/JLUGINfkhpj8EtSYwx+SWrMxMGfZHuSO5JsTbJlSH+SfCjJvUluT3LapOuUJI1vWp/jf3FVfW9E37nAyd3jBcBHumdJ0gwsx6meC4Arq+crwBFJjl2G9UqShpjGEX8B1ycp4M+qauNA/3HAd/rmd3Rt9/cPSrIeWA+watWqKZQlaRpm9a1hvzG870zjiP+MqjqN3imdNyR50UB/hrymntRQtbGq5qtqfm5ubgplSZKGmTj4q2pn97wLuBo4fWDIDuCEvvnjgZ2TrleSNJ6Jgj/J05IcvjANnA1sGxh2DfCa7tM9LwQerqr7kSTNxKTn+I8Brk6ysKz/WVWfT/K7AFV1KXAdcB5wL/BD4D9MuE5J0gQmCv6qug84dUj7pX3TBbxhkvVIkqbHb+5KUmMMfklqjMEvSY0x+CWpMQfcPXdneW9SSdPjfYb3HY/4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXmgLtkgyRN6kC/wbxH/JLUmLGDP8kJSW5McleSO5O8aciYM5M8nGRr93j7ZOVKkiY1yamex4G3VNWt3Q3Xb0myuar+z8C4v6mq8ydYjyRpisY+4q+q+6vq1m76UeAu4LhpFSZJ2jemco4/yWrg+cBXh3T/WpLbknwuyXMWWcb6JFuSbNm9e/c0ypIkDTFx8Cc5DPg08OaqemSg+1bgWVV1KvCnwF+NWk5Vbayq+aqan5ubm7QsSdIIEwV/koPphf7Hq+ozg/1V9UhVPdZNXwccnOToSdYpSZrMJJ/qCXAZcFdVfWDEmF/uxpHk9G593x93nZKkyU3yqZ4zgN8C7kiytWv7Q2AVQFVdCrwaeH2Sx4EfAWurqiZYpyRpQmMHf1XdDGQPYy4BLhl3HZKk6fObu5LUGINfkhpj8EtSYwx+SWqMwS9JjTH4JakxBr8kNcbgl6TGGPyS1BiDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXG4Jekxhj8ktSYSW+2fk6Su5Pcm2TDkP5Dk3yy6/9qktWTrE+SNLlJbrZ+EPBh4FzgFODCJKcMDLsIeKiq/gnwQeC9465PkjQdkxzxnw7cW1X3VdVPgE8AFwyMuQC4opv+FHBWkkXv0ytJ2rfGvtk6cBzwnb75HcALRo2pqseTPAw8Hfje4MKSrAfWd7OPJbl7zLqOHrb8Fcg6p29/qdU6p2t/qRP2UGsmOyfyrKUOnCT4hx251xhjeo1VG4GNE9TTW2GyparmJ13Ovmad07e/1Gqd07W/1Akrp9ZJTvXsAE7omz8e2DlqTJKnAP8YeHCCdUqSJjRJ8H8dODnJiUkOAdYC1wyMuQZY102/GvhiVQ094pckLY+xT/V05+zfCHwBOAjYVFV3JnkXsKWqrgEuAz6W5F56R/prp1H0Hkx8umiZWOf07S+1Wud07S91wgqpNR6AS1Jb/OauJDXG4Jekxhwwwb+ny0fMWpLtSe5IsjXJlq7tqCSbk9zTPR85g7o2JdmVZFtf29C60vOhbh/fnuS0Gdf5jiTf7fbp1iTn9fW9tavz7iQvX8Y6T0hyY5K7ktyZ5E1d+4rap4vUuRL36VOTfC3JbV2t7+zaT+wuBXNPd2mYQ7r2mVwqZpE6L0/y9337dE3XPrPfJ6pqv3/Qe3P574CTgEOA24BTZl3XQI3bgaMH2t4HbOimNwDvnUFdLwJOA7btqS7gPOBz9L6f8ULgqzOu8x3A7w8Ze0r3M3AocGL3s3HQMtV5LHBaN3048M2unhW1TxepcyXu0wCHddMHA1/t9tX/AtZ27ZcCr++mfw+4tJteC3xyxnVeDrx6yPiZ/T4dKEf8S7l8xErUf0mLK4BXLncBVXUTT/5uxai6LgCurJ6vAEckOXaGdY5yAfCJqvpxVf09cC+9n5F9rqrur6pbu+lHgbvofYN9Re3TReocZZb7tKrqsW724O5RwEvoXQoGnrxPl/1SMYvUOcrMfp8OlOAfdvmIxX6IZ6GA65Pckt7lKQCOqar7ofeLCDxjZtU90ai6VuJ+fmP3b/KmvlNlK6LO7hTD8+kd+a3YfTpQJ6zAfZrkoCRbgV3AZnr/cfygqh4fUs8TLhUDLFwqZtnrrKqFffrfun36wSSHDtbZWbZ9eqAE/5IvDTFDZ1TVafSuZvqGJC+adUFjWGn7+SPArwBrgPuB93ftM68zyWHAp4E3V9Ujiw0d0rZstQ6pc0Xu06r6aVWtoXeFgNOBX12knpnVOlhnkucCbwX+GfDPgaOAP5h1nQdK8C/l8hEzVVU7u+ddwNX0fngfWPjXrnveNbsKn2BUXStqP1fVA90v2s+AP+fnpx5mWmeSg+mF6cer6jNd84rbp8PqXKn7dEFV/QD4Er1z4kekdymYwXpmfqmYvjrP6U6rVVX9GPgfrIB9eqAE/1IuHzEzSZ6W5PCFaeBsYBtPvKTFOuCzs6nwSUbVdQ3wmu7TCC8EHl44fTELA+dDX0Vvn0KvzrXdpztOBE4GvrZMNYXeN9bvqqoP9HWtqH06qs4Vuk/nkhzRTf8j4KX03pO4kd6lYODJ+3TZLxUzos5v9P3BD733Ifr36Wx+n5brXeR9/aD3Dvk36Z37e9us6xmo7SR6n4i4DbhzoT565x1vAO7pno+aQW1X0fuX/v/ROwK5aFRd9P41/XC3j+8A5mdc58e6Om6n90t0bN/4t3V13g2cu4x1/kt6/67fDmztHuettH26SJ0rcZ8+D/jbrqZtwNu79pPo/fG5F/hL4NCu/and/L1d/0kzrvOL3T7dBvwFP//kz8x+n7xkgyQ15kA51SNJWiKDX5IaY/BLUmMMfklqjMEvSY0x+CWpMQa/JDXm/wPVgNxVcAbxqQAAAABJRU5ErkJggg==\n", "text/plain": [ "<matplotlib.figure.Figure at 0x7ff3e1a98780>" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "hand = ps2.get_frequency_dict('helloworld')\n", "random.seed(1)\n", "for n in (10, 100):\n", " for N in (1, 10, 100):\n", " plot_mc_hand(hand, N, n)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "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.4" } }, "nbformat": 4, "nbformat_minor": 2 }