{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def load_data(N):\n", " \"\"\"\n", " Get MNIST images of 3 only, as in the textbook.\n", " :N:\n", " :return:\n", " \"\"\"\n", " mnist = np.load('../data/mnist.npz')\n", " data = mnist['x_train'][mnist['y_train']==3]\n", " data = data[:N]\n", " data = data.reshape([N, 28 * 28])\n", " return data" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "N = 1000\n", "X = load_data(N)\n", "\n", "# sample mean\n", "X_ = np.mean(X, axis=0)\n", "\n", "# sample covariance matrix\n", "S = 1/N * np.dot((X-X_).T, (X-X_))" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAD8CAYAAAC4nHJkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAERxJREFUeJzt3VuMXfV1x/Hfwvju8WWMPQy2AxjbQLBUB1moUlGVKiWi\nKBLkBYWHylVRnIc0aqQ+FNGHIlWVUNWkylMkR6CYKiWpBAgeolYBVaGVqgiDKBi7Nr6BPfbM+G5j\nbGwPqw+znU5gzlrD2efq//cjWTNz1uxz/rNnft7nnP/N3F0AynNDtxsAoDsIP1Aowg8UivADhSL8\nQKEIP1Aowg8UivADhSL8QKFu7OSDmRnDCYE2c3ebyffVuvKb2YNmtsfM9pnZE3XuC0BnWbNj+81s\nlqS9kh6QdETSG5Iec/ddwTFc+YE268SV/z5J+9z9gLtflvRzSQ/XuD8AHVQn/KskHZ7y9ZHqtt9h\nZlvNbIeZ7ajxWABarO1v+Ln7NknbJJ72A72kzpV/RNKaKV+vrm4D0AfqhP8NSevN7HYzmyPpW5Je\naU2zALRb00/73f2qmf2FpH+XNEvSs+7+XstaBqCtmu7qa+rBeM0PtF1HBvkA6F+EHygU4QcKRfiB\nQhF+oFCEHygU4QcKRfiBQhF+oFCEHygU4QcKRfiBQhF+oFAdXbob0zOLJ2HVqWfH3nBD/P9/dnw2\nKzSqZ8d++umnbXtscOUHikX4gUIRfqBQhB8oFOEHCkX4gUIRfqBQ9PNX6vSHz549Ozx27ty5YX3+\n/PlhfenSpU3XlyxZEh67YMGCsD5r1qywPjExEdYvXbrUsHb27Nnw2HbWo3ZJ0tWrV8P69TCGgCs/\nUCjCDxSK8AOFIvxAoQg/UCjCDxSK8AOFqtXPb2aHJJ2XNCHpqrtvbkWjmmxLWM/mrc+ZMyesz5s3\nr2Ft2bJl4bErVqwI66tXrw7ra9euDeu33357w9ott9wSHrty5cqwfuON8Z/IJ598EtbPnDnTsHbo\n0KHw2H379oX1/fv3h/WDBw82rI2NjYXHZmMIrodxAK0Y5PNH7n6iBfcDoIN42g8Uqm74XdKrZvam\nmW1tRYMAdEbdp/33u/uIma2U9Csz+193f33qN1T/KfAfA9Bjal353X2k+jgu6SVJ903zPdvcfXM3\n3wwE8HlNh9/MFprZwLXPJX1d0s5WNQxAe9V52j8k6aWqi+1GSf/i7v/WklYBaLumw+/uByT9Xgvb\nUkvWj5/NuY/68SVp+fLlDWtr1qwJj73jjjvCetRPL0kbNmwI69HjR+2WpIGBgbCejX/IxgF8/PHH\nDWvDw8PhsYODg2E9W4sg+pvI+umz8QsXLlwI6/3Qz09XH1Aowg8UivADhSL8QKEIP1Aowg8U6rpZ\nujub0pt1SWXLa0ddgVk3Ybb8dabO9NJs6mp2XrLutEWLFoX16Nxk3bPZVOmsq/D48eMNa0ePHg2P\nzX6nFy9eDOvZ9uK9gCs/UCjCDxSK8AOFIvxAoQg/UCjCDxSK8AOFKqafP6tnUzCjLZ1PnTpV677P\nnTsX1g8cOBDWo+nK2c+dTXVeuHBhWM/62qNlybMxAtl5y8YoROMIsjEGdf+e+gFXfqBQhB8oFOEH\nCkX4gUIRfqBQhB8oFOEHCnXd9PNn86cvX74c1uuMA8iWgc7m4x87diysZ33xUX93tpZANl9/1apV\nYX3x4sVhPTo3WduyOfPnz58P69Gy4dG4DUm6cuVKWO+H+foZrvxAoQg/UCjCDxSK8AOFIvxAoQg/\nUCjCDxQq7ec3s2clfUPSuLtvrG4blPQLSbdJOiTpUXc/3b5m5rK53xMTE2E96/eN+quzY7NtrrN+\n/Oz4+fPnN6xlW3Rn9Wz78HXr1oX1bJxAJFtbf3R0NKyPjIw0rJ0+Hf+5Zlt098MW3JmZXPl/KunB\nz9z2hKTX3H29pNeqrwH0kTT87v66pM8uVfOwpO3V59slPdLidgFos2Zf8w+5+7UxqaOShlrUHgAd\nUntsv7u7mTV8AWRmWyVtrfs4AFqr2Sv/mJkNS1L1cbzRN7r7Nnff7O6bm3wsAG3QbPhfkbSl+nyL\npJdb0xwAnZKG38yel/Tfku40syNm9rikpyU9YGbvS/rj6msAfSR9ze/ujzUofa3Fbakl63fN5l9n\nx0f1uved9eMvWbIkrEd96XfddVd47D333BPW77zzzrCerdsfzdn/4IMPwmOzdQ4OHz4c1sfGxhrW\nsjUWsvUfSunnB3AdIvxAoQg/UCjCDxSK8AOFIvxAoa6bpbvrdNVJ+ZbNUXdetgT13Llzw/qKFSvC\n+t133x3W77333oa1TZs2hceuX78+rA8NxdM2sp99fLzh4M906e1s6/Nsa/MLFy40rGXLrdftvu0H\nXPmBQhF+oFCEHygU4QcKRfiBQhF+oFCEHyjUddPPn8m24K4j2iJbkhYtWhTWs2mx2bTajRs3Nqxl\nU3pXrlwZ1rMxCnWWuM7GCCxcuLBWPVrSPPudtfPvpVdw5QcKRfiBQhF+oFCEHygU4QcKRfiBQhF+\noFDF9PPXnX8dzffP+oSzPuWsLz2bW/7RRx81rEXbVEvSmTNnwnr2s2XrIERboy9YsCA8Ntve+/jx\n42H95MmTDWvZWgLZ+IXrYb4/V36gUIQfKBThBwpF+IFCEX6gUIQfKBThBwqV9vOb2bOSviFp3N03\nVrc9Jenbkq51tD7p7r9sVyNbIeuvzuaWR8dnfb4XL14M69FW0pK0d+/esB718+/atSs8Nvu5M9la\nBTfffHPDWraWQLZnwJe+9KWwPjo62rB24sSJ8NjonEr5uv/R+IZeMZMr/08lPTjN7f/k7puqfz0d\nfACfl4bf3V+XFG+dAqDv1HnN/z0ze8fMnjWzZS1rEYCOaDb8P5a0VtImScck/aDRN5rZVjPbYWY7\nmnwsAG3QVPjdfczdJ9z9U0k/kXRf8L3b3H2zu29utpEAWq+p8JvZ1OVmvylpZ2uaA6BTZtLV97yk\nr0q6ycyOSPpbSV81s02SXNIhSd9pYxsBtEEafnd/bJqbn2lDW1JRX3vdfvw6/d3ZY2f9/HXn3O/Z\ns6dhbd68eeGx2VoD2Zz7Zcvi93o3bNjQsJatY5Dd9+DgYNPHL168ODx2fHw8rNdZx6BXMMIPKBTh\nBwpF+IFCEX6gUIQfKBThBwrVV0t3R11qdZfHnj17dliPlmLOuvqy+oULF8J6Nr20znnJujijba4l\nafXq1WE9mpab/VxLliwJ63V+p9nPlXWRZkt/9wOu/EChCD9QKMIPFIrwA4Ui/EChCD9QKMIPFKqv\n+vkjc+bMCesDAwNhPev3zfrqI1euXAnr2ZTfbLvoaAxCtsR03fERdcYJZGMrsmmx2XmNZFNy6y5p\n3g+48gOFIvxAoQg/UCjCDxSK8AOFIvxAoQg/UKi+6uevM2896+e/6aabwnrUX1136e5su+hLly6F\n9WgcQNa2bIvt4eHhsL5u3bqwHm2jnS3NnY1vyObUR2McsvvOxkdk27L3A678QKEIP1Aowg8UivAD\nhSL8QKEIP1Aowg8UKu3nN7M1kp6TNCTJJW1z9x+Z2aCkX0i6TdIhSY+6++n2NTXus87mX2dzx5cu\nXRrWV65c2bCW9ZVH8+0l6eTJk2H93LlzYf3y5csNa9l8/OXLl4f1W2+9NayvX78+rK9ataphLRub\nMTY2FtazrctPn27855jtGZCNrSiln/+qpL9y9y9L+n1J3zWzL0t6QtJr7r5e0mvV1wD6RBp+dz/m\n7m9Vn5+XtFvSKkkPS9pefdt2SY+0q5EAWu8LveY3s9skfUXSbyQNufuxqjSqyZcFAPrEjMf2m9ki\nSS9I+r67n5v6+tvd3cymfWFrZlslba3bUACtNaMrv5nN1mTwf+buL1Y3j5nZcFUfljQ+3bHuvs3d\nN7v75lY0GEBrpOG3yUv8M5J2u/sPp5RekbSl+nyLpJdb3zwA7TKTp/1/IOlPJb1rZm9Xtz0p6WlJ\n/2pmj0v6QNKj7Wni/6uzRHW2zHPWdTM4ONiwlnWHZVtNZ8tIZz9bVK/bxblixYqwvnjx4rAe/c6O\nHj0aHnvkyJGw/uGHH4b1kZGRhrVTp06Fx2bTsK+Hrr40/O7+X5IadbB/rbXNAdApjPADCkX4gUIR\nfqBQhB8oFOEHCkX4gUL11dLdUd9qNgUzm/45Pj7tAMXfipb2zpa3zpaojqa9Snlf+rx58xrWsn7+\numMMzp49G9YPHjzYsLZz587w2N27d4f1vXv3hvVonEA03VeKp0lL+TTtfsCVHygU4QcKRfiBQhF+\noFCEHygU4QcKRfiBQvVVP3+d+fzZds7Z3PHosbMxBtnS29ky0tl6AdE4gKyfP5u3Pjo6Gtb37dsX\n1vfs2dOwtn///vDYw4cPh/U6S55nW3RPTEyE9esBV36gUIQfKBThBwpF+IFCEX6gUIQfKBThBwpl\nnZyX3GhLrxbdd1ivu4V3NGd+YGAgPDabj5+tnZ/df7RFeHZesv7sbL5+tv59VK+7TXbWVx+t/5Ct\nu9/P8/XdPf6lV7jyA4Ui/EChCD9QKMIPFIrwA4Ui/EChCD9QqLSf38zWSHpO0pAkl7TN3X9kZk9J\n+rak49W3Punuv0zuq2c7T+uME8jWvs/q2WNn9XbKxgFk/eVRvW5fenZ8P/fV1zHTfv6ZhH9Y0rC7\nv2VmA5LelPSIpEclfeTu/zjTRhH+5h6b8Dd3POGPpSv5uPsxSceqz8+b2W5J8RYzAHreF3rNb2a3\nSfqKpN9UN33PzN4xs2fNbNo9qcxsq5ntMLMdtVoKoKVmPLbfzBZJ+rWkv3f3F81sSNIJTb4P8Hea\nfGnw58l99OzzMJ72T4+n/f2npWP7zWy2pBck/czdX6weYMzdJ9z9U0k/kXRfs40F0Hlp+G3ysvOM\npN3u/sMpt0/dmvabkuItVwH0lJm823+/pP+U9K6ka8/hnpT0mKRNmnzaf0jSd6o3B6P7KvN5GNBB\nLevqayXCD7Qf8/kBhAg/UCjCDxSK8AOFIvxAoQg/UCjCDxSK8AOFIvxAoQg/UCjCDxSK8AOFIvxA\noQg/UKh0Ac8WOyHpgylf31Td1ot6tW292i6JtjWrlW27dabf2NH5/J97cLMd7r65aw0I9GrberVd\nEm1rVrfaxtN+oFCEHyhUt8O/rcuPH+nVtvVquyTa1qyutK2rr/kBdE+3r/wAuqQr4TezB81sj5nt\nM7MnutGGRszskJm9a2Zvd3uLsWobtHEz2znltkEz+5WZvV99nHabtC617SkzG6nO3dtm9lCX2rbG\nzP7DzHaZ2Xtm9pfV7V09d0G7unLeOv6038xmSdor6QFJRyS9Iekxd9/V0YY0YGaHJG129673CZvZ\nH0r6SNJz7r6xuu0fJJ1y96er/ziXuftf90jbntIX3Lm5TW1rtLP0n6mL566VO163Qjeu/PdJ2ufu\nB9z9sqSfS3q4C+3oee7+uqRTn7n5YUnbq8+3a/KPp+MatK0nuPsxd3+r+vy8pGs7S3f13AXt6opu\nhH+VpMNTvj6i3try2yW9amZvmtnWbjdmGkNTdkYalTTUzcZMI925uZM+s7N0z5y7Zna8bjXe8Pu8\n+919k6Q/kfTd6ultT/LJ12y91F3zY0lrNbmN2zFJP+hmY6qdpV+Q9H13Pze11s1zN027unLeuhH+\nEUlrpny9urqtJ7j7SPVxXNJL6r3dh8eubZJafRzvcnt+q5d2bp5uZ2n1wLnrpR2vuxH+NyStN7Pb\nzWyOpG9JeqUL7fgcM1tYvREjM1so6evqvd2HX5G0pfp8i6SXu9iW39ErOzc32llaXT53Pbfjtbt3\n/J+khzT5jv9+SX/TjTY0aNdaSf9T/Xuv222T9LwmnwZe0eR7I49LWi7pNUnvS3pV0mAPte2fNbmb\n8zuaDNpwl9p2vyaf0r8j6e3q30PdPndBu7py3hjhBxSKN/yAQhF+oFCEHygU4QcKRfiBQhF+oFCE\nHygU4QcK9X8wOS6YYmy38gAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# mean vector, figure 12.3\n", "plt.imshow(X_.reshape([28, 28]), cmap='gray', interpolation='none')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# eigen-decompose sample covmat (SVD on X-X_ works the same)\n", "L, U = np.linalg.eig(S) # eigenvalues, eigenvectors" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/yibo/anaconda3/lib/python3.5/site-packages/ipykernel/__main__.py:3: ComplexWarning: Casting complex values to real discards the imaginary part\n", " app.launch_new_instance()\n", "/home/yibo/anaconda3/lib/python3.5/site-packages/ipykernel/__main__.py:4: ComplexWarning: Casting complex values to real discards the imaginary part\n" ] } ], "source": [ "# somehow numpy returns complex results (but with zero imaginary parts); we can safely cast to floats b/c\n", "# we know the covmat is PSD and always has real eigen-values/vectors.\n", "L = np.float64(L)\n", "U = np.float64(U)\n", "\n", "# sort the eigen-value/vector pairs by decreasing eigenvalues\n", "idx = L.argsort()[::-1] \n", "L = L[idx]\n", "U = U[:,idx]" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEJCAYAAABbtPF6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEqtJREFUeJzt3XuMXOV9xvHvY9bcvICh3rgWoDiuSFWaNk67tdqGpkSo\nBNy0JmqDcCPqSGmcC6BS5Y8g0ipWlKqoykVR1SKZQGPShJQ0IFxKUoHVipBGwIIMOJAAJaaxZeO1\nHMc2N2P71z/mOJksO+cdzzszZ9bv85FW3j3vnDm/Oetnz8z5nYsiAjMrz7ymCzCzZjj8ZoVy+M0K\n5fCbFcrhNyuUw29WKIffrFAO/5BJekrSNkm/OoRl/YuknZL2SXpa0l/UPPbAjK/Dkv6hGjtJ0s2S\nnpe0X9JmSZcOof6zJN0p6cVq2X82Y/y/Jb3SVvMPBl3T8cThH763AE8DfzqEZd0ALIuI04E/Bj4t\n6Tdne2BEjB/9An4ReBn4ejU8BvwI+H3gDOCvgdslLe21MEnrJK1LPOwfgYPAYuB9wI2z/NG8uq32\nX+61nhI5/EMWEYeBB4BfH8KytkTES0d/rL5+qYtZ/wTYBXy7ep4XI2JdRGyNiCMRcTfwQ+Cnf0gk\nfVDSk5J+Iumbkt6QU7ukBVUdfxMRByLiAeAu4Mqc57WfcfiHTNIpwGrgrcc4392S9nb4urtmvn+S\n9BLwfWAHcE8Xi1sD3Bodjv2WtBh4M/C96ufrgQ/TencxAWwHPn0sr28WbwYORcTTbdMeA2Zu+f9O\n0m5J35F0YeYyi+LwD9/fAtuAZZLGASSdIemh6nPrW2abKSLeHRELO3y9u9PCIuKjwGnA7wF3AK/W\nFSfpjbTe3m/oMD4f+AqwISK+X23hPwGsjohnI+IgcDPwW/WrIWkc2Ddj2j5ar+WojwPLgLOB9cC/\nS+rmnY3h8A+VpN8B3kvr7exPgF+rhl4C/hD4t0EsNyIOV2+bzwE+knj4lcADEfHDmQOS5gFfpvU5\n/Opq8kXAScBDR9+JAN+i9fpmzn9322OuA66refdyADh9xrQzgP1tr+vBiNgfEa9GxAbgO8DKxOuz\nyljTBZRC0snAPwMfjog9kh6j9bn/uxHxGjAtqW7+b9Laes/m2xHRzd73MdKf+f+c1o7CmcsXrS36\nYmBlVTPAWcCdEfHe1MLb36Ec3dkXEes6PPxpYEzSeRHxTDXtrVQfNTotAui8Eu3neMs/PJ8C/ici\n/qP6eTPH8Lk/Ii5t3yM/4+t1wZf0BklXSBqXdIKkd9Ha17Cp0zIk/S6tt9Bfn2X4RuBXgD+KiJfb\npj8KvFPSb1TPcbqkVar7S9aFiHiR1seUT0laIOkCWvsUvlwtZ6Gkd0k6WdKYpPcB76D1rsO64PAP\ngaQVtN7u/1Xb5M0Mdo9/0HqLvw34MfAZ4NqI2FjV9M1qR127NcAdEbG/fWK1H+BDwHJgZ1tf/X0R\n8V1af9i+IekA8CRwSaedhcfoo8AptDoPXwU+EhFHt/zzae1UnAZ2A9cAl83YQWg15It5jA5JXwI+\nExFbmq7Fjn/e8o8ISfcAFwM3SXp/w+VYAbzlNyuUt/xmhXL4zQrl8JsVaqgH+SxatCiWLl06zEWa\nFWXr1q3s3r27q2MsssIv6RLgC8AJwBcj4nVHhrVbunQpUw89lLNIM6sxuWJF14/t+W2/pBNonW99\nKXA+sFrS+b0+n5kNV85n/hXAsxHxXHUm19eAVf0py8wGLSf8Z9O6ustR26ppP0fSWklTkqamp6cz\nFmdm/TTwvf0RsT4iJiNicmJiYtCLM7Mu5YR/O3Bu28/nVNPMbA7ICf/DwHmS3iTpROAKYGN/yjKz\nQeu51RcRhyRdDfwnrVbfLW2nW5rZiMvq80fEPXR3QUgzGzE+vNesUA6/WaEcfrNCOfxmhXL4zQrl\n8JsVyjftsFpHjtPtwzyONF1C447P36yZJTn8ZoVy+M0K5fCbFcrhNyuUw29WKLf6jgPHaztukFLr\nrIRWoP/XmBXK4TcrlMNvViiH36xQDr9ZoRx+s0I5/GaFcp+/Hw4danTxtX/Bx+p/xS+9Uv/3/5VX\n6pd94ED9eN2qSa22ROksWFA/ftppvT/32FjednEuHCfgLb9ZoRx+s0I5/GaFcvjNCuXwmxXK4Tcr\nlMNvVij3+Y9KNZ3rGt6pZnjOc3ejrmk9Pl4766mJ8bHxE2vHc44D2L+/ft7Dh+vHTz65fnzhwt7G\nILnakscJzAVZL0HSVmA/cBg4FBGT/SjKzAavH3+/3hkRu/vwPGY2RP7Mb1ao3PAHcJ+kRyStne0B\nktZKmpI0NT09nbk4M+uX3PBfEBHLgUuBqyS9Y+YDImJ9RExGxOTExETm4sysX7LCHxHbq393AXcC\nK/pRlJkNXs/hl7RA0mlHvwcuBrb0qzAzG6ycvf2LgTslHX2er0bEt/pS1SDk9trrGtZ79+Y9d6q2\nVEO7rimdeu7E+Ikn1/8XWbSofvtxwgn1i6+TWq2pawnUSfXpU6u86D5/RDwHvLWPtZjZELnVZ1Yo\nh9+sUA6/WaEcfrNCOfxmhRqphkXWbZMzW1rJdlzdeO6lu1Pnj6bG685PTc2bkmoFJtbb4jM798xO\nOaX+dOFUO2134nSynLOwU+OpVuBc4C2/WaEcfrNCOfxmhXL4zQrl8JsVyuE3K5TDb1aokerzD/S2\nxrn3g85p7OY+d+I60wfp3C8/lNmvnpd7j+6a8dMXLaqdde/J9ccB5JwunL5Fd+/PPVd4y29WKIff\nrFAOv1mhHH6zQjn8ZoVy+M0K5fCbFer46Wbm9tJzjgPIvc5zYv4jY/X9bmpKTy163qGD9Q9I9fEz\nTnxPva7cO5fPn995rIRLc6d4y29WKIffrFAOv1mhHH6zQjn8ZoVy+M0K5fCbFaqAbmYlt3E7yMZv\n4rlT1zk4sa7Rnzp+IdXHz70OQs21CHburJ811eev6+MnFu0+P11s+SXdImmXpC1t086SdK+kZ6p/\nzxxsmWbWb9287f8ScMmMadcBmyLiPGBT9bOZzSHJ8EfE/cCeGZNXARuq7zcAl/W5LjMbsF53+C2O\niB3V9zuBxZ0eKGmtpClJU9PT0z0uzsz6LXtvf0QEEDXj6yNiMiImJyYmchdnZn3Sa/hfkLQEoPp3\nV/9KMrNh6DX8G4E11fdrgLv6U46ZDUuymynpNuBCYJGkbcAngRuA2yV9AHgeuHyQRQ7FIBu7qV75\nIMdT86Zk3lNgz4HO5+zv3dtLQV0vunZ8fLx+3hL6/MmXGBGrOwxd1OdazGyIfHivWaEcfrNCOfxm\nhXL4zQrl8JsVqoCGxhyQe9ps3fy5lzRPjO87UL/9qDstN7e0BQt6n7+EVl6Kt/xmhXL4zQrl8JsV\nyuE3K5TDb1Yoh9+sUA6/WaGK6XYeyfw7l7p89kDlHAeQaGgP+jbZdb32VB8/57khr5ff6O97SLzl\nNyuUw29WKIffrFAOv1mhHH6zQjn8ZoVy+M0Kddz0+XP7+I0a4O3Dc9dLbi++Tu5VxVOrrYRefY45\nnBgzy+HwmxXK4TcrlMNvViiH36xQDr9ZoRx+s0LNqT5/k738umWPcj85WVtusz01f82BAmNjmddY\nOHSw/gE5Nw3IHc+Q+n/er/9vybUv6RZJuyRtaZu2TtJ2SZurr5V9qcbMhqabP71fAi6ZZfrnI2J5\n9XVPf8sys0FLhj8i7gf2DKEWMxuinA9d10h6vPpYcGanB0laK2lK0tT09HTG4sysn3oN/43AMmA5\nsAP4bKcHRsT6iJiMiMmJiYkeF2dm/dZT+CPihYg4HBFHgJuAFf0ty8wGrafwS1rS9uN7gC2dHmtm\noynZrJR0G3AhsEjSNuCTwIWSlgMBbAU+1I9imuzj57S7k/3qxLXxU0b5OILUCf8HD3VeNwcOZD01\np5L4peXcdGB8vH58gH3+Yf2+k68gIlbPMvnmAdRiZkPkw3vNCuXwmxXK4TcrlMNvViiH36xQc+qU\n3hy5t5oepFTXaKCneKZu4Z1Ydmq97t3beSzV6lu0qH781LFEq69uAYM+pXeArcB+8ZbfrFAOv1mh\nHH6zQjn8ZoVy+M0K5fCbFcrhNyvU6Dcju5TqN6fGX3utfnz+/M5jqZZu6tTU1OnEOS3j3NOkU734\nuj4+wO7dvS871edP/lJzL0t+nPOW36xQDr9ZoRx+s0I5/GaFcvjNCuXwmxXK4Tcr1Ej1+VPnpdf1\nrFMt3VRLOOcy0qmrPA+yjw956yW3j5+6A9urr3YeS/XxFy6sH+fZRHF1v/TUwRcF8JbfrFAOv1mh\nHH6zQjn8ZoVy+M0K5fCbFcrhNytUN7foPhe4FVhM65bc6yPiC5LOAv4VWErrNt2XR8SPc4oZ5C26\nc8/3rxtP9dJT/erc8/3rxnP7+KnxnNd+zjn1887b9n/1D0i9uLoDKAq4Ln9KN2k7BHwsIs4Hfhu4\nStL5wHXApog4D9hU/Wxmc0Qy/BGxIyIerb7fDzwFnA2sAjZUD9sAXDaoIs2s/47pfbakpcDbgAeB\nxRGxoxraSetjgZnNEV2HX9I48A3g2ojY1z4WEUFrf8Bs862VNCVpajp1ILiZDU1X4Zc0n1bwvxIR\nd1STX5C0pBpfAuyabd6IWB8RkxExOTEx0Y+azawPkuGXJOBm4KmI+Fzb0EZgTfX9GuCu/pdnZoPS\nTb/i7cCVwBOSNlfTrgduAG6X9AHgeeDywZTYndzLZ6e6Ri+/3NsYpNthuWeX1rUhc69unaotdVpu\nXTtv3tbn6mfOue435J2HfRy08lKSrzAiHgDUYfii/pZjZsPiI/zMCuXwmxXK4TcrlMNvViiH36xQ\nDr9ZoUaqmZlz6e5UWzb38to58w66157Tzk49d+p05Hm7Zz2w82c2b+s8lloxub/UuuJzV8xxwFt+\ns0I5/GaFcvjNCuXwmxXK4TcrlMNvViiH36xQI9XnT6k7DmBsrP7vWO556Sed1Hms7jbUkHcMAcCC\nBfXjZ57ZeezEV/Z1HoT0OfNbM26DnZJ7kEGqV183XsD5+ine8psVyuE3K5TDb1Yoh9+sUA6/WaEc\nfrNCOfxmhSqm2Zl7x+W6lnRuHz+17HmHDtY/oO6mA6kbEqSKzzlnHupXXM6FCqCI22gPkrf8ZoVy\n+M0K5fCbFcrhNyuUw29WKIffrFAOv1mhko1QSecCtwKLgQDWR8QXJK0DPghMVw+9PiLuGVShoyzZ\np0/cjyB7ATnXpx8099pHVje/mUPAxyLiUUmnAY9Iurca+3xEfGZw5ZnZoCTDHxE7gB3V9/slPQWc\nPejCzGywjukzv6SlwNuAB6tJ10h6XNItkma9mJSktZKmJE1NT0/P9hAza0DX4Zc0DnwDuDYi9gE3\nAsuA5bTeGXx2tvkiYn1ETEbE5MTERB9KNrN+6Cr8kubTCv5XIuIOgIh4ISIOR8QR4CZgxeDKNLN+\nS4ZfkoCbgaci4nNt05e0Pew9wJb+l2dmg9LN3v63A1cCT0jaXE27HlgtaTmt9t9W4EMDqbBLue20\nutt/58p97qzX5labddDN3v4HAM0yVGRP3+x44SP8zArl8JsVyuE3K5TDb1Yoh9+sUA6/WaHcBK5k\nn3ZrNsd4y29WKIffrFAOv1mhHH6zQjn8ZoVy+M0K5fCbFUoRMbyFSdPA822TFgG7h1bAsRnV2ka1\nLnBtvepnbW+MiK6ulzfU8L9u4dJUREw2VkCNUa1tVOsC19arpmrz236zQjn8ZoVqOvzrG15+nVGt\nbVTrAtfWq0Zqa/Qzv5k1p+ktv5k1xOE3K1Qj4Zd0iaQfSHpW0nVN1NCJpK2SnpC0WdJUw7XcImmX\npC1t086SdK+kZ6p/Z71HYkO1rZO0vVp3myWtbKi2cyX9l6QnJX1P0l9W0xtddzV1NbLehv6ZX9IJ\nwNPAHwDbgIeB1RHx5FAL6UDSVmAyIho/IETSO4ADwK0R8ZZq2t8DeyLihuoP55kR8fERqW0dcKDp\n27ZXd5Na0n5beeAy4P00uO5q6rqcBtZbE1v+FcCzEfFcRBwEvgasaqCOkRcR9wN7ZkxeBWyovt9A\n6z/P0HWobSRExI6IeLT6fj9w9Lbyja67mroa0UT4zwZ+1PbzNhpcAbMI4D5Jj0ha23Qxs1gcETuq\n73cCi5ssZhbJ27YP04zbyo/Muuvldvf95h1+r3dBRCwHLgWuqt7ejqRofWYbpV5tV7dtH5ZZbiv/\nU02uu15vd99vTYR/O3Bu28/nVNNGQkRsr/7dBdzJ6N16/IWjd0iu/t3VcD0/NUq3bZ/ttvKMwLob\npdvdNxH+h4HzJL1J0onAFcDGBup4HUkLqh0xSFoAXMzo3Xp8I7Cm+n4NcFeDtfycUblte6fbytPw\nuhu5291HxNC/gJW09vj/L/CJJmroUNcy4LHq63tN1wbcRutt4Gu09o18APgFYBPwDHAfcNYI1fZl\n4AngcVpBW9JQbRfQekv/OLC5+lrZ9LqrqauR9ebDe80K5R1+ZoVy+M0K5fCbFcrhNyuUw29WKIff\nrFAOv1mh/h/UhCb60TO15AAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEJCAYAAABbtPF6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEqpJREFUeJzt3XuMXPV5xvHvg72+4IWA641xAcWhMm1JGpx26zYKpTQ0\nFGhUQG0oboQcCcWkBdS0tAoijWL1oqIohPJHSuUUGkMIFJUgKIJUQFpRUgosyICBcik1xZbBa1GD\nDQFfePvHHKNhmTm/Yc+cObP7ez7SyLvndy7vzvrZMzPvuSgiMLP8HNR0AWbWDIffLFMOv1mmHH6z\nTDn8Zply+M0y5fCbZcrhHzBJT0naIukjNW9nvqSrJb0gaZekjZJOK5n/ZyX9UNKrkp6TdFaX+VZI\nelPSd+ur/p1tLZZ0i6TXi5/j96aM/1tRy+7i8XTdNc0mDv/gfRR4BvidmrczF3gR+FXgA8CfATdJ\nWj51RklzgVuB24HFwFrgu5KO7bDebwEPVS1O0jpJ6xKzfQvYAywFPgdc1eGP5oURMVo8frpqXTlx\n+AcsIvYD9wEfq3k7r0fEuojYHBFvR8TtwP8Av9Bh9p8BfhK4IiL2R8QPgR8B57bPJOkcYCdwz9QV\nSPqCpCeLVw53SvpglfolLQJ+G/hqROyOiPto/YE6t3xJ65XDP2CSFgKrgePf53K3S9rZ5XF7D8sv\nBY4Fnuh1k7RepRxY/lDgz4E/7rDuS4EvAr8FjAFbgb/scTvdHAvsi4hn2qY9Ckzd8/+1pB2SfiTp\npIrbzIrDP3h/BWwBjpE0CiBplaT7Jd0r6QZJI1MXiojPRMRhXR6fKdtgsb7rgQ0R8V8dZnka2A78\nqaQRSafQertwcNs8fwFcHRFbpqz7g8BXgNUR8VxE7AGuBn6xx+ejm1HgtSnTXgMOafv+y8AxwJHA\neuCfJf1Uxe1mw+EfIEmfAD5L6+Xsq8DPFUMvAp+KiBOBzcAZfdzmQcB1tN47X9hpnojYC5wJ/Cbw\nEnAxcBOtP1JIWgn8OnBFh8VPBuYDDx54JQL8gNbPN7WW29vmuQS4pOTVy27g0CnTPgDsaqv7gYjY\nFRFvRcQGWm9VTu/+bFi7uU0XkAtJC4B/AL4YEa9IepTW+/77I2Jb26x7gLc7LH8n8CtdVv/vEfGe\nT/IlidZeeClwehHyjiLiMVp7+wPL/gewofj2JGA58L+tVTIKzJF0XLH+WyLis93W3baNd16hHPiw\nLyLWdZn9GWCupBUR8Wwx7XjK37YErbcr1ouI8GMAD+DrwDVt3/8N8LdT5vkQcD8w0qdt/h3wn8Bo\nD/N+DFhA66X+n9D6cHB+MXYwcETb4xvAP9F6f/8JYAfw88W8h9J65aLE9tYB6xLz3AjcACwCTqD1\nauIjxdhhwG8UNc+l1Q14HTi26d/1THn4Zf8ASFpF6+X+H7VN3kjbJ/7FB2rXAZ+Pkj30+9jmh4Dz\ngZXAS2298M8V43cWH9QdcC6wjdZ7/5OBT0fEWwAR8UZEvHTgQesl+ZsRMRkR99P6IPBmSbuBJ4FT\no0hoRX8ALCxq+h7w+xFxYM8/QutDxUlaf3wuAs6Md39AaCXUn9+RVVH02W8DLo+I97TRzOrgPf9w\nWA38EvDV4qi13226IJv9vOc3y5T3/GaZcvjNMuXwm2VqoAf5LFmyJJYvXz7ITZplZfPmzezYsaOn\nA50qhV/SqcCVwBzg7yPisrL5ly9fzsSDD1bZpJmVGF+1qud5p/2yX9IcWudbnwYcB6wuDvc0sxmg\nynv+VcBzEfF8tM7kupE+npBiZvWqEv4jaZ2NdsCWYtq7SForaULSxOTkZIXNmVk/1f5pf0Ssj4jx\niBgfGxure3Nm1qMq4d8KHN32/VHFNDObAaqE/yFghaQPS5oHnEPr5BQzmwGm3eqLiH2SLgT+hVar\n75q20y3NbMhV6vNHxB3AHX2qxcwGyIf3mmXK4TfLlMNvlimH3yxTDr9Zphx+s0z5ph1W6u1Zun84\n6L33RcnO7PzNmlmSw2+WKYffLFMOv1mmHH6zTDn8Zplyq68wW1ta1lnq951DK9D/480y5fCbZcrh\nN8uUw2+WKYffLFMOv1mmHH6zTM2aPn+Tffp9+8rH33yz2vK7d09/+dS6U+MpcxP/g8rGFy0qX3bh\nwvLxBQumv+2qZsNxAt7zm2XK4TfLlMNvlimH3yxTDr9Zphx+s0w5/GaZmjV9/qqq9MNTffxUnz41\nnlr/3r3l43UaGSkfT/Xiq6hyjEHKbOjjp1QKv6TNwC5gP7AvIsb7UZSZ1a8fe/5fi4gdfViPmQ2Q\n3/ObZapq+AO4W9LDktZ2mkHSWkkTkiYmJycrbs7M+qVq+E+IiJXAacAFkk6cOkNErI+I8YgYHxsb\nq7g5M+uXSuGPiK3Fv9uBW4BV/SjKzOo37fBLWiTpkANfA6cAm/pVmJnVq8qn/UuBWyQdWM/3IuIH\nfalqCJX1+VN99tQxBKl+9Oho+XhZLz217jlzyser2r+/vnVXuRbBbOjTVzXt8EfE88DxfazFzAbI\nrT6zTDn8Zply+M0y5fCbZcrhN8vUUJ3SO8ynUZa1zFKntR522PTXDenTYufNLXleql43PCXRh3zj\nze6/09Smf/zj6RTUo7qvWT4DeM9vlimH3yxTDr9Zphx+s0w5/GaZcvjNMuXwm2VqqJqVw9rHT42n\nli3tw0P1e3ynxquo8drbVVvtKT5tt5z3/GaZcvjNMuXwm2XK4TfLlMNvlimH3yxTDr9Zpoaqz19F\n1Z5u6loCZSr38as2vKtcuzsx/sa+eaXjb+4sX32dtw+fBafUN8p7frNMOfxmmXL4zTLl8JtlyuE3\ny5TDb5Yph98sU+6UFiodJ1Bnn74He/Z1/xueOtV/d6JPv3t3+XiVwwjmzy9fduHCatuu/YIBM1xy\nzy/pGknbJW1qm7ZY0l2Sni3+PbzeMs2s33p52f8d4NQp0y4B7omIFcA9xfdmNoMkwx8R9wKvTJl8\nBrCh+HoDcGaf6zKzmk33A7+lEbGt+PolYGm3GSWtlTQhaWJycnKamzOzfqv8aX9EBBAl4+sjYjwi\nxsfGxqpuzsz6ZLrhf1nSMoDi3+39K8nMBmG64b8NWFN8vQa4tT/lmNmgJPv8km4ATgKWSNoCfA24\nDLhJ0nnAC8DZdRY59BIN59S1AlK9+F27ysdffXX6667aCh8dLR8vO4RhZKR82YqXIsCHsZRLPjsR\nsbrL0Ml9rsXMBsiH95plyuE3y5TDb5Yph98sUw6/WabcC+lRlUt7p9ptOxOn1e7YUT5edtrtnDnl\ny6ZOm02dbVzlbOTUZb1T46k25YIF3X9nqTZh3ZeCH4bbh3vPb5Yph98sUw6/WaYcfrNMOfxmmXL4\nzTLl8Jtlatb0+Zvsq6b6zak+f2o8pey02tQpt4sWlY+njgOo8rOXnYrcy7pTvfqyn/2ww8qXLTtG\nYLaY/T+hmXXk8JtlyuE3y5TDb5Yph98sUw6/WaYcfrNMzZo+f1VVztdPSfWjUz3n1HhZP3see8oX\nrnpt78QJ/W8sOLjrWOr236lLlu/fXz5eVnrqOgfVLxs+/LznN8uUw2+WKYffLFMOv1mmHH6zTDn8\nZply+M0yNQu6lb1pso+furZ98hry+xK9+rKGeaqZnurjVzyp/uCSgxSOOGJx6bJV73dQNp76naSu\nY5C6TsJMkEyEpGskbZe0qW3aOklbJW0sHqfXW6aZ9Vsvu8PvAKd2mH5FRKwsHnf0tywzq1sy/BFx\nL/DKAGoxswGq8kb4IkmPFW8LDu82k6S1kiYkTUxOTlbYnJn103TDfxVwDLAS2AZc3m3GiFgfEeMR\nMT42NjbNzZlZv00r/BHxckTsj4i3gW8Dq/pblpnVbVrhl7Ss7duzgE3d5jWz4ZTs80u6ATgJWCJp\nC/A14CRJK4EANgPn11jjO+rs1afa2VVUvhd8lV583Seep2rbsaPr0KFHlTfblyzpfi0AqHYIQ9XL\nGMwGyf8ZEbG6w+Sra6jFzAbIh/eaZcrhN8uUw2+WKYffLFMOv1mmsjmlt6p6W4GJ24tXuRd11cJT\nPbEq/bbEsqOj5a2+1Gm5ZatPPS1795aPV719+DDwnt8sUw6/WaYcfrNMOfxmmXL4zTLl8JtlyuE3\ny9QM6Eb2R6ovW+UUz1RPOHUZ6PSlv+eVjpceB1D10typZnqV4wgS6656iELZ0zIT+vB1857fLFMO\nv1mmHH6zTDn8Zply+M0y5fCbZcrhN8tUNt3Oqu3usnPDU8cIvP56+XjJXax7UnY9gLlzE8cIVNt0\n+l7VJb38N/aV17az+1W/e1J2GEHVW3TPhuMEvOc3y5TDb5Yph98sUw6/WaYcfrNMOfxmmXL4zTLV\nyy26jwauBZbSuiX3+oi4UtJi4B+B5bRu0312RPxflWLqvAV3Suqc/LJefurS9al1p5ZPtdLLes7z\n55cvOzJS3mtPHSeQsq/kZ9u5s3zZ1POSUtbLX7SofNnZ0MdP6SVt+4CLI+I44JeBCyQdB1wC3BMR\nK4B7iu/NbIZIhj8itkXEI8XXu4CngCOBM4ANxWwbgDPrKtLM+u99vc6WtBz4OPAAsDQithVDL9F6\nW2BmM0TP4Zc0CtwMfCkiXmsfi4ig9XlAp+XWSpqQNDE5OVmpWDPrn57CL2mEVvCvj4jvF5NflrSs\nGF8GbO+0bESsj4jxiBgfGxvrR81m1gfJ8EsScDXwVER8s23oNmBN8fUa4Nb+l2dmdemlofFJ4Fzg\ncUkbi2mXApcBN0k6D3gBOLueEvsj1bpJncKZOi23zFtvVVt3qiVW1s4bGSlfNnVqa2q8iqqtvFQL\ntGz8kEPKl636cx/E29VWMADJ8EfEfYC6DJ/c33LMbFB8hJ9Zphx+s0w5/GaZcvjNMuXwm2XK4TfL\n1FCduJjqjVY55Td9G+zy8SVLpr9sqk+fuvR3FcN8amqFq373NF62/mF+XgbFe36zTDn8Zply+M0y\n5fCbZcrhN8uUw2+WKYffLFMzqttZ5znS8xLPxLwF3e/hfWii33zEEeWXv071+aveXrxMqt/dZD+8\nzmsJzITz7evmPb9Zphx+s0w5/GaZcvjNMuXwm2XK4TfLlMNvlqkZ1eevVZVmeqJRf1Bi3QfX2chP\nqXqhg9R4pQMF/N+zTt7zm2XK4TfLlMNvlimH3yxTDr9Zphx+s0w5/GaZSjZSJR0NXAssBQJYHxFX\nSloHfAGYLGa9NCLuqKvQ2lU5sb3qSe919vGrqvqz+QL5Q6uX38w+4OKIeETSIcDDku4qxq6IiG/U\nV56Z1SUZ/ojYBmwrvt4l6SngyLoLM7N6va/3/JKWAx8HHigmXSTpMUnXSDq8yzJrJU1ImpicnOw0\ni5k1oOfwSxoFbga+FBGvAVcBxwArab0yuLzTchGxPiLGI2J8bGysDyWbWT/0FH5JI7SCf31EfB8g\nIl6OiP0R8TbwbWBVfWWaWb8lwy9JwNXAUxHxzbbpy9pmOwvY1P/yzKwuvXza/0ngXOBxSRuLaZcC\nqyWtpNX+2wycX0uFM4HbYTYD9fJp/32AOgzN3J6+mfkIP7NcOfxmmXL4zTLl8JtlyuE3y5TDb5Yp\nh98sUw6/WaYcfrNMOfxmmXL4zTLl8JtlyuE3y5TDb5YpRcTgNiZNAi+0TVoC7BhYAe/PsNY2rHWB\na5uuftb2oYjo6Xp5Aw3/ezYuTUTEeGMFlBjW2oa1LnBt09VUbX7Zb5Yph98sU02Hf33D2y8zrLUN\na13g2qarkdoafc9vZs1pes9vZg1x+M0y1Uj4JZ0q6WlJz0m6pIkaupG0WdLjkjZKmmi4lmskbZe0\nqW3aYkl3SXq2+LfjPRIbqm2dpK3Fc7dR0ukN1Xa0pH+V9KSkJyT9YTG90eeupK5GnreBv+eXNAd4\nBvg0sAV4CFgdEU8OtJAuJG0GxiOi8QNCJJ0I7AaujYiPFtO+DrwSEZcVfzgPj4gvD0lt64DdTd+2\nvbib1LL228oDZwKfp8HnrqSus2ngeWtiz78KeC4ino+IPcCNwBkN1DH0IuJe4JUpk88ANhRfb6D1\nn2fgutQ2FCJiW0Q8Uny9CzhwW/lGn7uSuhrRRPiPBF5s+34LDT4BHQRwt6SHJa1tupgOlkbEtuLr\nl4ClTRbTQfK27YM05bbyQ/PcTed29/3mD/ze64SIWAmcBlxQvLwdStF6zzZMvdqebts+KB1uK/+O\nJp+76d7uvt+aCP9W4Oi2748qpg2FiNha/LsduIXhu/X4ywfukFz8u73het4xTLdt73RbeYbguRum\n2903Ef6HgBWSPixpHnAOcFsDdbyHpEXFBzFIWgScwvDdevw2YE3x9Rrg1gZreZdhuW17t9vK0/Bz\nN3S3u4+IgT+A02l94v/fwFeaqKFLXccAjxaPJ5quDbiB1svAvbQ+GzkP+AngHuBZ4G5g8RDVdh3w\nOPAYraAta6i2E2i9pH8M2Fg8Tm/6uSupq5HnzYf3mmXKH/iZZcrhN8uUw2+WKYffLFMOv1mmHH6z\nTDn8Zpn6f4OibHxW9JXlAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEJCAYAAABbtPF6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEltJREFUeJzt3XuMXPV5xvHvY4y52ca43jgWQThIpCpNi2lXVttQIEJJ\nwGoLURsUgqijojiRkrRR80cQaRUralVU5SJUtbROoTGIklIFCqUkVYJaJaQtsCBzb4AkprFl8FrU\nwSa49pq3f8zZZFh2zm+Y35w5s/t7PtJoZ85lzjuz++yZOe+5KCIws/IsabsAM2uHw29WKIffrFAO\nv1mhHH6zQjn8ZoVy+M0K5fCPmKSnJO2S9PMNL+c4STdIek7SAUk7JF086LSSDs65HZX0Fw2/htWS\n7pD0clXbB+aM/3dJh7pq+m6T9Sw2Dv/ovR14GvidhpezFPghcD5wMvBHwG2S1g8ybUQsn70BbwZe\nAf5x0OIkbZW0NTHZXwKHgbXAFcD18/zT/FhXbT87aD0lcvhHLCKOAvcBv9jwcl6OiK0RsTMiXo2I\nu4EfAL+cM23lt4G9wLdnB0j6kKQnJf1I0tckvSmnfkknVcv544g4GBH3AXcCV+Y8r/2Uwz9ikk4A\nLgfOfoPz3S1pf4/b3X3MvxZ4G/DEEKbdDNwU1b7hkq4BPgL8FjAB7Ab+pK8X1tvbgJmIeLpr2CPA\n3DX/n0naJ+k7ki7IXGZRlrZdQIH+FNgFnC9peUQcrMJ2B3AEOApcERF7umeKiN8YdIGSjgVuAbZH\nxH/nTCvpdDpfD66qHr8J+DRwTkQ8Ww27AfirQeutLAdemjPsJWBF1+NPAU/S+WrwfuCfJW2IiO9l\nLrsIXvOPkKRfBd5H5+Psj4BfqEbtA86NiPOBm6iCNaRlLgFuphOQjw1h2iuB+yLiB9XjC4HjgAdm\nP4kAX6fz+uY+/91d01wNXF3z6eUgsHLOsJOBA7MPIuL+iDgQEf8XEduB7wCb6l6j/ZTX/CMi6Xjg\n74CPRMSLkh6h873/P6vtALNWMM/HbUlfA369x9N/OyJetyVfkoAb6Gww2xQRR2rq63fa3wWu7Xq8\nGrgjIt7X67lndX96md3YFxFbe0z+NLBU0pkR8Uw17Gzqv7YEoFQd1uE1/+h8FviPiPiX6vEOur73\nS9og6X46a9yH584cERd3b3Gfc5u3hQdcD/wc8JsR8UqivuS0kn4NOJXXbuV/GHinpF+qplkp6ZLq\nn8nAIuJl4Hbgs5JOknQunW0KN1fLWSXpPZKOl7RU0hXAeXQ+dVg/IsK3hm/ARjpbz0/uGvZBOh+f\n5057GfDXQ1jm6XTWhIfofISevV1Rjf8acE0/03Y9598AN8+zrN+vXt9BOtszru+jvq3A1sQ0q4F/\nAl4G/gf4QNe4CeBBOl8D9gP/Bbyr7d/1QrqpeiOtRZKWRcTh6v57gPdExB+2XJYtcv7OPx42SPoc\nnS39h4Dfa7keK4DX/GaF8gY/s0I5/GaFcvjNCjXSDX5r1qyJ9evXj3KRZkXZuXMn+/bt62sfi6zw\nS7oIuA44BvjbiLi2bvr169cz9cADOYs0sxqTGzf2Pe3AH/slHUPneOuLgbOAyyWdNejzmdlo5Xzn\n3wg8GxHfr3ZQ+QpwyXDKMrOm5YT/VDpnf5m1qxr2GpK2SJqSNDU9PZ2xODMbpsa39kfEtoiYjIjJ\niYmJphdnZn3KCf9u4LSux2+phpnZApAT/geBMyW9VdIyOmdSuWs4ZZlZ0wZu9UXEjKSPAf9Kp9V3\nY0Qkzw9nZuMhq88fEfcA9wypFjMbIe/ea1Yoh9+sUA6/WaEcfrNCOfxmhXL4zQrlE3harVcX6fph\nCa+2XULrFudv1sySHH6zQjn8ZoVy+M0K5fCbFcrhNyuUW32VhdzSmplpZ16ApQ3+BTX53Knfdwmt\nwIX7F29mWRx+s0I5/GaFcvjNCuXwmxXK4TcrlMNvVqhF0+dvs0+f2ytPzX/oUP34gwcHnze39uXL\n68cfd1zvcStW5C07tR9Ak736xbCfgNf8ZoVy+M0K5fCbFcrhNyuUw29WKIffrFAOv1mhFk2fv+m+\nal1fN/e48yaPW0/14VOOPz5v/JKZw71HpnZCONTgyQZShTf5SxkTWa9Q0k7gAHAUmImIyWEUZWbN\nG8a/t3dGxL4hPI+ZjZC/85sVKjf8AXxT0kOStsw3gaQtkqYkTU1PT2cuzsyGJTf850bEBuBi4KOS\nzps7QURsi4jJiJicmJjIXJyZDUtW+CNid/VzL3AHsHEYRZlZ8wYOv6STJK2YvQ+8G3h8WIWZWbNy\ntvavBe6QNPs8fx8RXx9KVT3U9dqz+/yJA9tr/0sm+9V540+sn7t+fOqA/dT4VL87p19edyICyK+9\nrrZVq+rnTewgsWQR7Acw8CuIiO8DZw+xFjMbIbf6zArl8JsVyuE3K5TDb1Yoh9+sUGPVr2j09Nu5\nbaO6tlSqZZVz7m3Ir71OqmWVOiY4pxWYeu79++vH55yXPPeQ3txjpceA1/xmhXL4zQrl8JsVyuE3\nK5TDb1Yoh9+sUA6/WaHGqs+fOiy31f0AmpR7eGhdzzrVz071qxOHvh4+fmXt+LpW/MpVicOwc/ef\nqHvtbf6+x4TX/GaFcvjNCuXwmxXK4TcrlMNvViiH36xQDr9Zocaqz5+SdXruVC89Z3yqV57qtTd5\n+uxEbT+eWVY7fl/iEqypVnxdaStTh8TnnucgxyI4NXeK1/xmhXL4zQrl8JsVyuE3K5TDb1Yoh9+s\nUA6/WaEWfzOzX02epz2zj//q0vpefN1h7fsz+/Q5fXyANWtqRu7aVT9z6rz9KTnnOXCfHyTdKGmv\npMe7hq2W9A1Jz1Q/T2m2TDMbtn4+9n8ZuGjOsKuBeyPiTODe6rGZLSDJ8EfEt4AX5wy+BNhe3d8O\nXDrkusysYYNu8FsbEXuq+88Da3tNKGmLpClJU9PT0wMuzsyGLXtrf0QEEDXjt0XEZERMTkxM5C7O\nzIZk0PC/IGkdQPVz7/BKMrNRGDT8dwGbq/ubgTuHU46ZjUqymSnpVuACYI2kXcBngGuB2yRdBTwH\nXNZkkSPRYF83p08PcCjRa69rh+ceEp9qh7/5zfXjV87M3VbcJdXHTy08Nb7umgOp/TYK6PMnX2FE\nXN5j1IVDrsXMRsi795oVyuE3K5TDb1Yoh9+sUA6/WaEWTT8jdfnurNN+9/H8TTpypH58XVcqcYVt\njjuufvwpieM1lx16qX6C/TW9xpxWXT/z17XzUvMmNP33Ngpe85sVyuE3K5TDb1Yoh9+sUA6/WaEc\nfrNCOfxmhVo0ff6UNvv0qZ7viUvrj6s9cUXiuNsTasanjtlNHU/8fGL+nNOSp/r4uYfdZhyWuxj6\n+Cle85sVyuE3K5TDb1Yoh9+sUA6/WaEcfrNCOfxmhVo0ff5U3zW3z1/7/Kled+51sHPGp/r4uefu\nzjmmPvf02KlLmze4blsM+wF4zW9WKIffrFAOv1mhHH6zQjn8ZoVy+M0K5fCbFWrR9PnbPF4/W+4x\n93Xzp3rpqT597vw558dPve6EJXW1N7yPwELYDyD5CiXdKGmvpMe7hm2VtFvSjuq2qdkyzWzY+vn3\n9mXgonmGfzEiNlS3e4Zblpk1LRn+iPgW8OIIajGzEcr5YvNxSY9WXwt6XtFN0hZJU5KmpqenMxZn\nZsM0aPivB84ANgB7gM/3mjAitkXEZERMTkxMDLg4Mxu2gcIfES9ExNGIeBX4ErBxuGWZWdMGCr+k\ndV0P3ws83mtaMxtPyT6/pFuBC4A1knYBnwEukLQBCGAn8OFhFDPO59bP0vR16FP7CdTJPfd9zn4A\nufs3pMbX1Zbcv2FZ/fhFIBn+iLh8nsE3NFCLmY3QAt4tzsxyOPxmhXL4zQrl8JsVyuE3K9RYHdLb\n9Om3m3ru2kNHodl2WUrDp8c+PFP/vtV141JX4F6Se0r0uvct8bpSv9MFfQh5ZeG/AjMbiMNvViiH\n36xQDr9ZoRx+s0I5/GaFcvjNCjVWff5xVt9yrv8fmmpXL00dPpoYn9PKT9U2k3lU7ZEjvcel6j4x\nNUHOIcE5+04sEl7zmxXK4TcrlMNvViiH36xQDr9ZoRx+s0I5/GaFWlB9/nYva9z7/2TuGahzD/fP\nOXN3at5XXqkff/Ro/fhjjuk9Lrl/wv7E8frpHSgSCyib1/xmhXL4zQrl8JsVyuE3K5TDb1Yoh9+s\nUA6/WaH6uUT3acBNwFo6l+TeFhHXSVoN/AOwns5lui+LiP/NKSZ1LvTaPn9usz0x/7K6k8xnXs45\n1UtvUqpPn9tKP+GE3uOWHXqpfubU7yylrjjvA9DXmn8G+GREnAX8CvBRSWcBVwP3RsSZwL3VYzNb\nIJLhj4g9EfFwdf8A8BRwKnAJsL2abDtwaVNFmtnwvaHv/JLWA+cA9wNrI2JPNep5Ol8LzGyB6Dv8\nkpYDXwU+ERGv+bIWEUFne8B8822RNCVpanp6OqtYMxuevsIv6Vg6wb8lIm6vBr8gaV01fh2wd755\nI2JbRExGxOTExMQwajazIUiGX5KAG4CnIuILXaPuAjZX9zcDdw6/PDNrSj/9jncAVwKPSdpRDbsG\nuBa4TdJVwHPAZc2UOCS5l3uuaTvVtgGB5ctXZi267vTXueoOuYX0Ga5T40/kx71H7ttfP3NunzHj\nEt252j38vD/JdyAi7gPUY/SFwy3HzEbFe/iZFcrhNyuUw29WKIffrFAOv1mhHH6zQo3VcY2N9kZz\nz49d14xPNOqXLa8fv3rVqvplr6pvph+e6f0/PPWyku956rDa/Ylefd17k+rjp3YiSOxf0WaffyHw\nmt+sUA6/WaEcfrNCOfxmhXL4zQrl8JsVyuE3K9TiaXbmHNvdj7qedOqA/H378sYnal/WZM869/TZ\ndVJ9+pw+PriXn+A1v1mhHH6zQjn8ZoVy+M0K5fCbFcrhNyuUw29WqHIaobn7AdQdc5+aN9Urz7x8\neO38uecxyD1xf9343D59g3381OXiF4PF/wrNbF4Ov1mhHH6zQjn8ZoVy+M0K5fCbFcrhNytUslEq\n6TTgJmAtEMC2iLhO0lbgQ8B0Nek1EXFPU4U2LtUzrju2PHXcuTWizV58o9eYGJF+9pKYAT4ZEQ9L\nWgE8JOkb1bgvRsTnmivPzJqSDH9E7AH2VPcPSHoKOLXpwsysWW/oc5Ok9cA5wP3VoI9LelTSjZJO\n6THPFklTkqamp6fnm8TMWtB3+CUtB74KfCIiXgKuB84ANtD5ZPD5+eaLiG0RMRkRkxMTE0Mo2cyG\noa/wSzqWTvBviYjbASLihYg4GhGvAl8CNjZXppkNWzL8kgTcADwVEV/oGr6ua7L3Ao8Pvzwza0o/\nW/vfAVwJPCZpRzXsGuBySRvotP92Ah9upMIRWayHcLbdklqo72vb79so9LO1/z5A84xauD19M1ug\n/5bNLJvDb1Yoh9+sUA6/WaEcfrNCOfxmhSrn1N0JJfR12+D3dXx5zW9WKIffrFAOv1mhHH6zQjn8\nZoVy+M0K5fCbFUoRMbqFSdPAc12D1gD7RlbAGzOutY1rXeDaBjXM2k6PiL7OlzfS8L9u4dJUREy2\nVkCNca1tXOsC1zaotmrzx36zQjn8ZoVqO/zbWl5+nXGtbVzrAtc2qFZqa/U7v5m1p+01v5m1xOE3\nK1Qr4Zd0kaTvSnpW0tVt1NCLpJ2SHpO0Q9JUy7XcKGmvpMe7hq2W9A1Jz1Q/571GYku1bZW0u3rv\ndkja1FJtp0n6N0lPSnpC0h9Uw1t972rqauV9G/l3fknHAE8D7wJ2AQ8Cl0fEkyMtpAdJO4HJiGh9\nhxBJ5wEHgZsi4u3VsD8HXoyIa6t/nKdExKfGpLatwMG2L9teXU1qXfdl5YFLgQ/S4ntXU9dltPC+\ntbHm3wg8GxHfj4jDwFeAS1qoY+xFxLeAF+cMvgTYXt3fTuePZ+R61DYWImJPRDxc3T8AzF5WvtX3\nrqauVrQR/lOBH3Y93kWLb8A8AvimpIckbWm7mHmsjYg91f3ngbVtFjOP5GXbR2nOZeXH5r0b5HL3\nw+YNfq93bkRsAC4GPlp9vB1L0fnONk692r4u2z4q81xW/ifafO8Gvdz9sLUR/t3AaV2P31INGwsR\nsbv6uRe4g/G79PgLs1dIrn7ubbmenxiny7bPd1l5xuC9G6fL3bcR/geBMyW9VdIy4P3AXS3U8TqS\nTqo2xCDpJODdjN+lx+8CNlf3NwN3tljLa4zLZdt7XVaelt+7sbvcfUSM/AZsorPF/3vAp9uooUdd\nZwCPVLcn2q4NuJXOx8AjdLaNXAX8DHAv8AzwTWD1GNV2M/AY8CidoK1rqbZz6XykfxTYUd02tf3e\n1dTVyvvm3XvNCuUNfmaFcvjNCuXwmxXK4TcrlMNvViiH36xQDr9Zof4fI3uUH8I6KMsAAAAASUVO\nRK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEJCAYAAABbtPF6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEyhJREFUeJzt3XuMXGd9xvHvs75kk9iO43qxXSfCGNJWKS2m3VqlpDQo\nKiQRaoLaRrg0dSWKubdU/EEUWhEhqqYVF6GKRjJNGodCuChJk9JABVErGloBS2QSk0BIU0fYOPFu\nQ2Ibx07s/PrHnA2TZeZ9x3tm5szu+3yk0e6ed86c38zus+fM/M5FEYGZlWes6QLMrBkOv1mhHH6z\nQjn8ZoVy+M0K5fCbFcrhNyuUwz9kkh6QtE/SLw5hWe+UNCXpuKQbM/c9Mud2UtLftY3/k6RHJR2S\n9KCkPxlC/Wsk3Sbpx5IekfQHc8b/Q9Kxtpq/N+iaFhOHf/heCjwI/N4QlvVD4IPADbk7RsSK2Ruw\nHngK+HzbXa4FNkfEKuB3gA9K+tX5FibpGknXZO72ceBpYB3wRuC6Dv8039lW+8/Pt54SOfxDFhEn\ngbuBXx7Csm6NiH8G/u8UZ/1d4CDwn22PtScijs7+WN1ePDsu6c2S7pf0pKQvSnpBndolnVnV8ZcR\ncSQi7gZuB66s87j2Ew7/kEk6HdgGvOwU5/uCpCe63L7Q5zK3AzfFnH2/Jf29pKPAd4EDwJ3V9KuB\nt9LaIpgA9tPa4qjj54ATEfFg27RvA3PX/H8taUbS1yRdWHOZRXH4h++vgH3AZkkr2gckbZM03Wmm\niHhdRKzucntdv4qT9ELgt4BdHWp4O7AS+E3gVuB4tYZ/H7AtIh6KiKeB64Ffq1nKCuDQnGmHquXP\nei+wGdgI7AT+RdKLsZ44/EMk6RXA79PanH0S+KW2sSXV2A+aqe45VwJ3R8T/dhqMiJPVJvg5wNuA\ni4DTgG/MbokAX6L1/J6nfesFuAq4KrH1cgRYNWfaWcDhtlq+HhGHI+J4ROwCvgZcOp8nXSKHf0gk\njQP/CLw1Ih6ntQnb/r5/G60P2J7tMv8XO3wiP3v7Yh9L/SM6rPU7WErrPf8a4LY5WyJnRcSFc2do\n33qh9QHitYmtlweBpZLOa5v2MuA7iZoCUA+1Gw7/MH0A+K+I+Nfq591U7/urtf4VwGe7zRwRl7R/\nIj/ndkmneSQtrf7pLAGWSBqXtLTbMiT9Bq1N6M/Pmf4CSW+QtELSEkmvpfXP6i7gHuDVkn6luu8q\nSZdJqhXCiPgxrbcWH5B0pqQLaH2m8MlqOaslvXb2OUl6I/AqWlsd1gOHfwgkbaW1Sf/nbZN385M1\n/x8Cn4uIjmv9Gv6CVsvuqmoZT1XTZrckrp5z/+3ArRFxeM70oLWJvw/4EfAh4N0RcUdE/Detf2y3\nSDoC3A9cPPfDwnl6O3A6rc7Dp4G3RcTsmn8ZrQ8Vp4EZ4F3A5XM+ILQE+WQezZP0N8DLaW3yvwLY\nFRF/2mxVttg5/CNG0lRETDZdhy1+Dr9Zofye36xQDr9ZoRx+s0J17fkOwtq1a2PTpk3DXKRZUfbu\n3cvMzExP+1jUCr+ki4GP0dqJ5B8i4trU/Tdt2sTUN75RZ5FmljC5dWvP9533Zn+1V9rHgUuA84Ft\nks6f7+OZ2XDVec+/FXgoIh6ujuT6DHBZf8oys0GrE/6NPP8ItH3VtOeRtKM6ldTU9HTHo1XNrAED\n/7Q/InZGxGRETE5MTAx6cWbWozrh3w+c2/bzOdU0M1sA6oT/m8B5kl4kaTnwBuCO/pRlZoM271Zf\nRJyQ9E7g32i1+m5oO9zSzEZcrT5/RNxJdRJHM1tYvHuvWaEcfrNCOfxmhXL4zQrl8JsVyuE3K9RQ\nj+e3hefZRbp+GOt8bZSiLM7frJllOfxmhXL4zQrl8JsVyuE3K5TDb1Yot/oqi7WlVdeJE80te+kA\n/zpzv+8SWoH+izcrlMNvViiH36xQDr9ZoRx+s0I5/GaFcvjNCrVo+vyLuU8/duxo+g6pZvyxY/Of\nF1iea/TnmvGp8RUrkrM+u/SM5HidfRByZdf9e1oI+wks3sSYWZLDb1Yoh9+sUA6/WaEcfrNCOfxm\nhXL4zQq1aPr8gzbInvLYkUPpO8zMpMefeKL7WK7PX3M/gOyTS/Xy169Pzjq2dm160ePp/QBST22Q\nuy8sFLWegqS9wGHgJHAiIib7UZSZDV4//n+9OiIyqyYzGzV+z29WqLrhD+Arkr4laUenO0jaIWlK\n0tT09HTNxZlZv9QN/wURsQW4BHiHpFfNvUNE7IyIyYiYnJiYqLk4M+uXWuGPiP3V14PAbcDWfhRl\nZoM37/BLOlPSytnvgdcAe/pVmJkNVp1P+9cBt0mafZxPR8SX+lJVF3WOsc71deuM1+4J1z05/vh4\n97HMMfO1i6/zwtX9pWTkdmFISb2kUHifPyIeBl7Wx1rMbIjc6jMrlMNvViiH36xQDr9ZoRx+s0KN\nVMOiyVZeri30zDPdx5YtS8+bawstX706fYdM3+nZxKGtqaN9exk/diQ9nusknrOp+9jYzMH0zBl1\nf6d1ZH+nI5WszrzmNyuUw29WKIffrFAOv1mhHH6zQjn8ZoVy+M0KNVLdyNxljQe5H0AdqX0AoIez\nYy9NP6/Dh9OnqH5yX/ex3Fm/6565+6yz0uOpXRhWpWfNFrd0RfYRuhrk38NC4TW/WaEcfrNCOfxm\nhXL4zQrl8JsVyuE3K5TDb1aokerzD9Ion2o512s/fnz+j5073n7JkvR47lwFGzemx1eRuPx4bieE\nzHkO6lx9fJT/HobFa36zQjn8ZoVy+M0K5fCbFcrhNyuUw29WKIffrFCLpttZt29b90rUdeQeu85V\ntHN9/tx47pICa8aPpu+w79HuY7kntn59cngmcR6DnNNOS4+XcInu7Jpf0g2SDkra0zZtjaQvS/p+\n9fXswZZpZv3Wy2b/jcDFc6ZdBdwVEecBd1U/m9kCkg1/RHwVeHzO5MuAXdX3u4DL+1yXmQ3YfD/w\nWxcRB6rvHwXWdbujpB2SpiRNTU9Pz3NxZtZvtT/tj4gAIjG+MyImI2JyYmKi7uLMrE/mG/7HJG0A\nqL7Wu9yqmQ3dfMN/B7C9+n47cHt/yjGzYcl2KyXdDFwIrJW0D3g/cC3wOUlvAh4Brhhkkf1Qt49f\np6+bO+48d97/nFRPOne8/tmZJu3yY4nj8QGeOJIeT+1IkOnjP7w3vW46UmPRp5+enncx9PFzsk8x\nIrZ1Gbqoz7WY2RB5916zQjn8ZoVy+M0K5fCbFcrhNytUAQ2N3tRp7Qz6cs+5w0tTLa3lJzKH3OZa\ndTmZY4KfHu9+Ge2Hvpt+6NyZvXOH5aZKq3OYdC9yl5PPXY5+GLzmNyuUw29WKIffrFAOv1mhHH6z\nQjn8ZoVy+M0KtaD6/KneaK6vOspWrkyP5/r8Y8cSvfwnnjj1gk5l4Zk+/5HE4uvuH3HyZHo8dSh1\n7jDrQe8HMAoWbmLMrBaH36xQDr9ZoRx+s0I5/GaFcvjNCuXwmxVqEXQrW3LHRze5H0DdnvHYiafT\nd0g1rQfdsM7sR7Am8fjjL+l+rD/AvswluHPH+6deltxpv32JbjNbtBx+s0I5/GaFcvjNCuXwmxXK\n4TcrlMNvVqhF0K1sqdvHH+S592ufA37p8uT42OrV3eet+bpkzy+fO19Aohl/BulG/Utesjk5njsm\nP9XLz81bwvH+2b8MSTdIOihpT9u0ayTtl7S7ul062DLNrN96WS3cCFzcYfpHI2JLdbuzv2WZ2aBl\nwx8RXwUeH0ItZjZEdd4QvkvSvdXbgrO73UnSDklTkqamp6drLM7M+mm+4b8O2AxsAQ4AH+52x4jY\nGRGTETE5MTExz8WZWb/NK/wR8VhEnIyIZ4FPAFv7W5aZDdq8wi9pQ9uPrwf2dLuvmY2mbDdS0s3A\nhcBaSfuA9wMXStoCBLAXeEs/ihnkMfe5Pn6dPn/dnm792gb3uo2PZ64zX+fJZ/YRGHv0h8nxdet+\nNjmeet3q9vkzlytYELK/uYjY1mHy9QOoxcyGyLv3mhXK4TcrlMNvViiH36xQDr9ZoRbAgYe9ybXD\ncq2b3PwL4RDNQci9buMr0qffHqvTb8u0Aid+Id3qe/LJ7mO5U3fXbb8uhL8Xr/nNCuXwmxXK4Tcr\nlMNvViiH36xQDr9ZoRx+s0ItgG5kb3It46eeSo+fPJkeX7Kk+9iyZel5R7nnm+tXP/NMevzw4fT4\nuonupxXPnvY780sdO3Y0OT4+fsZ8H3qgp3IfFV7zmxXK4TcrlMNvViiH36xQDr9ZoRx+s0I5/GaF\nGqkOdO5y0HVO7Z3r49fp6+b6/Ll9DOpKPbe6x6XnxsfH0+O1Hjy3g0Sthdd76FHed6NXXvObFcrh\nNyuUw29WKIffrFAOv1mhHH6zQjn8ZoXq5RLd5wI3AetoXZJ7Z0R8TNIa4LPAJlqX6b4iIn5Up5g6\nffy6LeHcceupXn7dnnBuP4Djx9PjqXPQD7qPvzpxuD6QP2a/xoMfPZb+e0k9t9wltge4C8HI6CVt\nJ4D3RMT5wK8D75B0PnAVcFdEnAfcVf1sZgtENvwRcSAi7qm+Pww8AGwELgN2VXfbBVw+qCLNrP9O\naTtb0ibg5cDXgXURcaAaepTW2wIzWyB6Dr+kFcAtwLsj4lD7WEQErc8DOs23Q9KUpKnp6elaxZpZ\n//QUfknLaAX/UxFxazX5MUkbqvENwMFO80bEzoiYjIjJiYmJftRsZn2QDb8kAdcDD0TER9qG7gC2\nV99vB27vf3lmNii9HJj4SuBK4D5Ju6tpVwPXAp+T9CbgEeCKwZTYm7qtvtz8qfEzxtOHIuf6aatW\npBd+6Ej6f3Sqtlwr77TT0uO5jbWxmY4bfD8xM9N9LPdLWbs2OXws00VMvS6551337yV3ePooyIY/\nIu4G1GX4ov6WY2bD4j38zArl8JsVyuE3K5TDb1Yoh9+sUA6/WaFG6gTEdU7dPehTKSd7+aljaiF/\nWGtuP4DMk1uVakrnXpdc7Y9lxnM7EqQOy12/Pjlrbv+G3GHYqV7+ypXpeRfDqblzvOY3K5TDb1Yo\nh9+sUA6/WaEcfrNCOfxmhXL4zQpVQDezpW7fNrWPwVidkwFA/f0Ejh1Lj6fUPTd3plefmv/oieXJ\nWXO7EJx+eno89bLX/ZUtBl7zmxXK4TcrlMNvViiH36xQDr9ZoRx+s0I5/GaFWlDdzNTx/nUu792L\nVM956fgZyXnH1mZe5lwvPdfHzzXEU2peq/rpE/O/THZOrrScEnr1dXjNb1Yoh9+sUA6/WaEcfrNC\nOfxmhXL4zQrl8JsVKtsJlXQucBOwDghgZ0R8TNI1wJuB6equV0fEnYMqdNByPeFUvzrfy04ft87S\nzPiK9H4EdeRqX1qjTw/50wUsVLlrTCwEvewGcQJ4T0TcI2kl8C1JX67GPhoRHxpceWY2KNnwR8QB\n4ED1/WFJDwAbB12YmQ3WKb3nl7QJeDnw9WrSuyTdK+kGSWd3mWeHpClJU9PT053uYmYN6Dn8klYA\ntwDvjohDwHXAZmALrS2DD3eaLyJ2RsRkRExOTEz0oWQz64eewi9pGa3gfyoibgWIiMci4mREPAt8\nAtg6uDLNrN+y4Zck4HrggYj4SNv0DW13ez2wp//lmdmg9PJp/yuBK4H7JO2upl0NbJO0hVb7by/w\nloFU2KO6rZfcIcGL9fDQxfq86loMrbycXj7tvxtQh6EF29M3M+/hZ1Ysh9+sUA6/WaEcfrNCOfxm\nhXL4zQrlLm+lhL6uWTuv+c0K5fCbFcrhNyuUw29WKIffrFAOv1mhHH6zQikihrcwaRp4pG3SWmBm\naAWcmlGtbVTrAtc2X/2s7YUR0dP58oYa/p9auDQVEZONFZAwqrWNal3g2uarqdq82W9WKIffrFBN\nh39nw8tPGdXaRrUucG3z1Uhtjb7nN7PmNL3mN7OGOPxmhWok/JIulvQ9SQ9JuqqJGrqRtFfSfZJ2\nS5pquJYbJB2UtKdt2hpJX5b0/eprx2skNlTbNZL2V6/dbkmXNlTbuZL+XdL9kr4j6c+q6Y2+dom6\nGnndhv6eX9IS4EHgt4F9wDeBbRFx/1AL6ULSXmAyIhrfIUTSq4AjwE0R8dJq2t8Cj0fEtdU/zrMj\n4r0jUts1wJGmL9teXU1qQ/tl5YHLgT+mwdcuUdcVNPC6NbHm3wo8FBEPR8TTwGeAyxqoY+RFxFeB\nx+dMvgzYVX2/i9Yfz9B1qW0kRMSBiLin+v4wMHtZ+UZfu0RdjWgi/BuBH7T9vI8GX4AOAviKpG9J\n2tF0MR2si4gD1fePAuuaLKaD7GXbh2nOZeVH5rWbz+Xu+80f+P20CyJiC3AJ8I5q83YkRes92yj1\nanu6bPuwdLis/HOafO3me7n7fmsi/PuBc9t+PqeaNhIiYn/19SBwG6N36fHHZq+QXH092HA9zxml\ny7Z3uqw8I/DajdLl7psI/zeB8yS9SNJy4A3AHQ3U8VMknVl9EIOkM4HXMHqXHr8D2F59vx24vcFa\nnmdULtve7bLyNPzajdzl7iNi6DfgUlqf+P8P8L4mauhS12bg29XtO03XBtxMazPwGVqfjbwJ+Bng\nLuD7wFeANSNU2yeB+4B7aQVtQ0O1XUBrk/5eYHd1u7Tp1y5RVyOvm3fvNSuUP/AzK5TDb1Yoh9+s\nUA6/WaEcfrNCOfxmhXL4zQr1/4QGUyGDQ/HwAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# visualize top 4 eigenvectors, figure 12.3\n", "for i in range(4):\n", " plt.figure()\n", " plt.imshow(U[:, i].reshape([28, 28]), vmin=-1, vmax=1, cmap='seismic', interpolation='none')\n", " plt.title('$\\lambda_%d=%.2e$' %((i+1), L[i]))" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAEYCAYAAACDV/v0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+cXHV97/HXe2aTTcgvEtiEkEQTJVohapAY0FprTSXR\nh23wFjRWJW0R7EP01t7eewv1KooXK71WrtxWelFSAoqQol5SLhQjaHu1JbBggARME01isgZ2yYYE\nAknY3c/943yHTIad2RDmzMxm38/HYx5z5nvO+c7nLGHfe77nO2cUEZiZmTVSodkFmJnZyOPwMTOz\nhnP4mJlZwzl8zMys4Rw+ZmbWcA4fMzNrOIePHZMk/YWkbzS7jlokbZX0282uw6wZ2ppdgFkeIuKL\nza5hpJD0I+CbEdHSYW+txWc+ZpYrSf4j117E4WPDlqSTJX1HUo+kLZL+Y9m6z0n6Ztnr8yVtk7RL\n0mfKh7wkFSRdIunnaf0qSVPSutmSQtJySb+U9KSkT5e9/3OlbVPb6WmbUZJeLeme1OeTkr4l6fgq\nx3K9pP9e9vodknYcybEO0td7JD0q6WlJXZL+c3mfaUjyyfQz+FDZfu2SvpyO8wlJfydpbNn6pZLW\nSdqbflZLJF0B/AbwN5KekfQ3aduQdLGkTcCmsp9jW1l/P5L00bT8B5J+IukqSU9J+oWkt6b27ZK6\nJS2v8c/BhhmHjw1LkgrAPwIPATOARcCnJC0eZNtTga8BHwKmA5PSPiWfBM4BfhM4GdgN/G1FN28D\nXpve57OSXhcRvwL+Dfi9su1+H7g1Ip4HBPxl6vN1wCzgc3kea3Id8LGImADMA+4pW3cScGLqZzlw\nraTXpnVfAl4DzAdOSdt8NtWwELgB+C/A8cDbga0R8Wng/wGfiIjxEfGJsvc6BzgTOPUID/VM4GHg\nBOAm4GbgzamWD5MF3Pgj7MtanMPHhqs3Ax0RcXlEHIyIXwBfB5YNsu25wD9GxI8j4iDZL9Tymxr+\nMfDpiNgREQfIAuLciuGiz0fEcxHxEFkIvDG13wR8EECS0vvfBBARmyNiTUQciIge4CtkAZfnsQI8\nD5wqaWJE7I6IByvWfybV9M/A/wXen2q/CPjTiOiNiKeBL5a9xwXAinQ8AxHRFRE/G6Luv0x9PXeE\nx7klIv4+IvqBW8jC+vJU6/eBg2RBZMcAj8XacPVK4GRJT5W1Fcn+Cq90MrC99CIinpW0q6Kv70ka\nKGvrB6aVvX68bPlZoPQX+HeA/yVpOtlZw0CpBknTgK+SDUtNIPtjb/eRHmBFfUd6rJCdif034EuS\nHgYuiYh/S+t2R8S+sm23kf18OoDjgAeyHAKyM7diWp4F3PES694+9CaHeaJs+TmAiKhs85nPMcLh\nY8PVdrK/lOcewbY7yYbMAEjXMU6o6OuPIuInlTtKml2r44jYLen7wAfIhtZujkO3iv8i2RnW6yOi\nV9I5wN9U6Wof2S//kpMq6jvSYyUi7geWShoFfAJYRRYeAJMljSsLoFcA64EnyX65nxYRXYN0ux14\ndbW3PIL20vsdB+xNyydhI5aH3Wy4ug94WtKfSxorqShpnqQ3D7LtrcDvpAvYo8mG1VS2/u+AKyS9\nEkBSh6SlL6GWm4DzyYb3biprnwA8A+yRNIPsekk164D3SJoi6STgU0dzrJJGS/qQpEnputNesrOx\ncp9P2/0G8F7gHyJigGwo7ypJU1NfM8quK10H/KGkRcomaMyQ9Gtp3RPAq2r9gNKwYxfw4VT/H1E9\nzGwEcPjYsJSuC7yX7OL4FrK/3L9BNpmgctsNZJMKbiY7C3oG6AYOpE2+CqwGvi/paeBesovfR2o1\nMBd4PF0TKvk88CZgD9m1le/W6ONGsmtJW4Hvk13zeMnHmnwE2CppL9n1rA+VrXucbOjvV8C3gD8u\nu3bz58Bm4N607w9IZ4wRcR/wh8BV6Xj+mWw4ELKf37mSdku6usYxXkgWwLuA04B/rbGtHePkL5Oz\nkSbNmHoKmBsRW5pdT6NIegfZh0FnNrsWM5/52Igg6XckHSdpHPBl4BGyswwzawKHj40US8mGmn5F\nNkS2LHzab9Y0HnYzM7OG85mPmZk1nD/nU8WJJ54Ys2fPbnYZZmbDygMPPPBkRHQMtZ3Dp4rZs2fT\n2dnZ7DLMzIYVSduOZDsPu5mZWcM5fMzMrOEcPmZm1nAOHzMzaziHj5mZNZzDx8zMGs7hY2ZmDefw\nqbMtT+7jK9/fyM49R/rNwWZmI4/Dp8627drH1fdsZuee/c0uxcysZTl86qyg7AsyBwZ8w1Yzs2oc\nPnVWLKTwcfaYmVXl8KmzdOLDgL+qwsysKodPnXnYzcxsaA6fOvOwm5nZ0Bw+dZayh34Pu5mZVeXw\nqbMXht0cPmZmVTl86qwUPuHwMTOrqmHhI2mMpPskPSRpg6TPp/bPSeqStC493lO2z6WSNkvaKGlx\nWfsZkh5J666Wst/4ktol3ZLa10qaXbbPckmb0mN5XsdZCp/+gbzewcxs+Gvk12gfAN4ZEc9IGgX8\nWNKdad1VEfHl8o0lnQosA04DTgZ+IOk1EdEPXANcCKwF7gCWAHcCFwC7I+IUScuAK4EPSJoCXAYs\nAAJ4QNLqiNhd74MspDj3sJuZWXUNO/OJzDPp5aj0qPUbeilwc0QciIgtwGZgoaTpwMSIuDeysa0b\ngHPK9lmZlm8FFqWzosXAmojoTYGzhiyw6s5Trc3MhtbQaz6SipLWAd1kYbA2rfqkpIclrZA0ObXN\nALaX7b4jtc1Iy5Xth+0TEX3AHuCEGn1V1neRpE5JnT09PUd1jJ5qbWY2tIaGT0T0R8R8YCbZWcw8\nsiG0VwHzgZ3AXzeypor6ro2IBRGxoKOj46j68FRrM7OhNWW2W0Q8BfwQWBIRT6RQGgC+DixMm3UB\ns8p2m5nautJyZfth+0hqAyYBu2r0VXee7WZmNrRGznbrkHR8Wh4LvAv4WbqGU/I+YH1aXg0sSzPY\n5gBzgfsiYiewV9JZ6XrO+cBtZfuUZrKdC9yTrgvdBZwtaXIa1js7tdWdP+djZja0Rs52mw6slFQk\nC71VEXG7pBslzSebfLAV+BhARGyQtAp4FOgDLk4z3QA+DlwPjCWb5VaaNXcdcKOkzUAv2Ww5IqJX\n0heA+9N2l0dEbx4H6anWZmZDa1j4RMTDwOmDtH+kxj5XAFcM0t4JzBukfT9wXpW+VgArXkLJR8VT\nrc3MhuY7HNSZp1qbmQ3N4VNnnmptZjY0h0+d+cvkzMyG5vCpM892MzMbmsOnzoq+5mNmNiSHT529\nMNXa2WNmVpXDp85KU619hwMzs+ocPnV26EOmDh8zs2ocPnXmqdZmZkNz+NSZp1qbmQ3N4VNnvsOB\nmdnQHD519sJUa2ePmVlVDp86k79MzsxsSA6fOpNEQZ5qbWZWi8MnBwXJEw7MzGpw+OSgIPnL5MzM\nanD45KBQ8LCbmVktDQsfSWMk3SfpIUkbJH0+tU+RtEbSpvQ8uWyfSyVtlrRR0uKy9jMkPZLWXS1l\nl/kltUu6JbWvlTS7bJ/l6T02SVqe57FmZz4OHzOzahp55nMAeGdEvBGYDyyRdBZwCXB3RMwF7k6v\nkXQqsAw4DVgCfE1SMfV1DXAhMDc9lqT2C4DdEXEKcBVwZeprCnAZcCawELisPOTqrSh5qrWZWQ0N\nC5/IPJNejkqPAJYCK1P7SuCctLwUuDkiDkTEFmAzsFDSdGBiRNwb2djWDRX7lPq6FViUzooWA2si\nojcidgNrOBRYdSf5DgdmZrU09JqPpKKkdUA3WRisBaZFxM60yePAtLQ8A9hetvuO1DYjLVe2H7ZP\nRPQBe4ATavRVWd9Fkjoldfb09Bz1cRYKnu1mZlZLQ8MnIvojYj4wk+wsZl7F+iA7G2qKiLg2IhZE\nxIKOjo6j7qfoqdZmZjU1ZbZbRDwF/JBs6OuJNJRGeu5Om3UBs8p2m5nautJyZfth+0hqAyYBu2r0\nlQt5qrWZWU2NnO3WIen4tDwWeBfwM2A1UJp9thy4LS2vBpalGWxzyCYW3JeG6PZKOitdzzm/Yp9S\nX+cC96SzqbuAsyVNThMNzk5tuSgWfGNRM7Na2hr4XtOBlWnGWgFYFRG3S/o3YJWkC4BtwPsBImKD\npFXAo0AfcHFE9Ke+Pg5cD4wF7kwPgOuAGyVtBnrJZssREb2SvgDcn7a7PCJ68zpQ3+HAzKy2hoVP\nRDwMnD5I+y5gUZV9rgCuGKS9E5g3SPt+4Lwqfa0AVry0qo9OQfKNRc3MavAdDnJQLMjDbmZmNTh8\nclAsiH5nj5lZVQ6fHBT8IVMzs5ocPjkoyMNuZma1OHxyUCz4xqJmZrU4fHLgqdZmZrU5fHLgMx8z\ns9ocPjkoeLabmVlNDp8cFORvMjUzq8Xhk4Oiv8nUzKwmh08OCr7mY2ZWk8MnB/4+HzOz2hw+OfBs\nNzOz2hw+OZDA2WNmVp3DJwfFgofdzMxqcfjkwLPdzMxqc/jkwLPdzMxqc/jkwLPdzMxqa1j4SJol\n6YeSHpW0QdKfpPbPSeqStC493lO2z6WSNkvaKGlxWfsZkh5J666WpNTeLumW1L5W0uyyfZZL2pQe\ny/M81uyaT57vYGY2vLU18L36gD+LiAclTQAekLQmrbsqIr5cvrGkU4FlwGnAycAPJL0mIvqBa4AL\ngbXAHcAS4E7gAmB3RJwiaRlwJfABSVOAy4AFQKT3Xh0Ru/M4UAl/n4+ZWQ0NO/OJiJ0R8WBafhp4\nDJhRY5elwM0RcSAitgCbgYWSpgMTI+LeyG6gdgNwTtk+K9PyrcCidFa0GFgTEb0pcNaQBVYusq/R\ndviYmVXTlGs+aTjsdLIzF4BPSnpY0gpJk1PbDGB72W47UtuMtFzZftg+EdEH7AFOqNFXZV0XSeqU\n1NnT03PUx+fZbmZmtTU8fCSNB74DfCoi9pINob0KmA/sBP660TWVRMS1EbEgIhZ0dHQcdT+Fgr9G\n28ysloaGj6RRZMHzrYj4LkBEPBER/RExAHwdWJg27wJmle0+M7V1peXK9sP2kdQGTAJ21egrF9ls\nt7x6NzMb/ho5203AdcBjEfGVsvbpZZu9D1ifllcDy9IMtjnAXOC+iNgJ7JV0VurzfOC2sn1KM9nO\nBe5J14XuAs6WNDkN652d2nJRKOBrPmZmNTRyttuvAx8BHpG0LrX9BfBBSfPJZqFtBT4GEBEbJK0C\nHiWbKXdxmukG8HHgemAs2Sy3O1P7dcCNkjYDvWSz5YiIXklfAO5P210eEb05HScFedjNzKyWhoVP\nRPwY0CCr7qixzxXAFYO0dwLzBmnfD5xXpa8VwIojrffl8Gw3M7PafIeDHBQ8283MrCaHTw6KBeET\nHzOz6hw+OSgIn/mYmdXg8MlBwdd8zMxqcvjkoOjZbmZmNTl8cuDZbmZmtTl8clBQNuEgHEBmZoNy\n+OSgrZB9nMmTDszMBufwyUGxmIVPn8PHzGxQDp8cFOUzHzOzWhw+OSiWht18zcfMbFAOnxy8ED79\nDh8zs8E4fHJQmnDgaz5mZoNz+OSgWMh+rAMedjMzG5TDJwfF9FP1mY+Z2eAcPjl44czH4WNmNiiH\nTw58zcfMrLaGhY+kWZJ+KOlRSRsk/UlqnyJpjaRN6Xly2T6XStosaaOkxWXtZ0h6JK27Wso+WCOp\nXdItqX2tpNll+yxP77FJ0vI8j7Xwwh0OBvJ8GzOzYauRZz59wJ9FxKnAWcDFkk4FLgHujoi5wN3p\nNWndMuA0YAnwNUnF1Nc1wIXA3PRYktovAHZHxCnAVcCVqa8pwGXAmcBC4LLykKu3Q7fXyesdzMyG\nt4aFT0TsjIgH0/LTwGPADGApsDJtthI4Jy0vBW6OiAMRsQXYDCyUNB2YGBH3Rnbnzhsq9in1dSuw\nKJ0VLQbWRERvROwG1nAosOqu+MKwm9PHzGwwTbnmk4bDTgfWAtMiYmda9TgwLS3PALaX7bYjtc1I\ny5Xth+0TEX3AHuCEGn3lwrfXMTOrreHhI2k88B3gUxGxt3xdOpNp2m9sSRdJ6pTU2dPTc9T9lG4s\n6vAxMxtcQ8NH0iiy4PlWRHw3NT+RhtJIz92pvQuYVbb7zNTWlZYr2w/bR1IbMAnYVaOvw0TEtRGx\nICIWdHR0HO1h+szHzGwIjZztJuA64LGI+ErZqtVAafbZcuC2svZlaQbbHLKJBfelIbq9ks5KfZ5f\nsU+pr3OBe9LZ1F3A2ZImp4kGZ6e2XHiqtZlZbW0NfK9fBz4CPCJpXWr7C+BLwCpJFwDbgPcDRMQG\nSauAR8lmyl0cEf1pv48D1wNjgTvTA7Jwu1HSZqCXbLYcEdEr6QvA/Wm7yyOiN68DLU048IdMzcwG\n17DwiYgfA6qyelGVfa4ArhikvROYN0j7fuC8Kn2tAFYcab0vR9FnPmZmNb3k8JF0edpvHbAuIv69\n7lUNc0V/jbaZWU0vOXwi4rOSpgHzgfdJOiUiLqx/acNXW7q3m8PHzGxwRzXsFhFPkF2wz+2i/XBW\n8F2tzcxqeknhI2kW2e1u5gGvB06LiAV5FDactfn7fMzMahpyqrWkj0n6V0lPAf8OfBQYTzat+fdz\nrm9Y8oQDM7PajuTM51LgA8CTZNOixwIrIuKXeRY2nBV9V2szs5qO5EOm742ItRHx84g4D/hb4B8l\n/akkfx/QIHxXazOz2oYMj4hYX/H6TrKvJZgC/CSnuoY1n/mYmdV2tLPdDgCfkXRjnes5Jviaj5lZ\nbS9r2MwfMB2cb69jZlabr9nkoHRXa5/5mJkNzuGTA3+fj5lZbQ6fHLT53m5mZjU5fHJQ8LCbmVlN\nDp8c+MzHzKw2h08O/JUKZma1OXxyIImCHD5mZtU4fHLSVijQ77tam5kNqmHhI2mFpG5J68vaPiep\nS9K69HhP2bpLJW2WtFHS4rL2MyQ9ktZdLWVX9yW1S7olta+VNLtsn+WSNqXH8kYcb7Egn/mYmVXR\nyDOf64Elg7RfFRHz0+MOAEmnAsvIvjtoCfA1ScW0/TXAhcDc9Cj1eQGwOyJOAa4Crkx9TQEuA84k\nuyfdZZIm1//wDlcsiL5+h4+Z2WAaFj4R8S9A7xFuvhS4OSIORMQWYDOwUNJ0YGJE3BsRAdwAnFO2\nz8q0fCuwKJ0VLQbWRERvROwG1jB4CNZVsSB/mZyZWRWtcM3nk5IeTsNypTOSGcD2sm12pLYZabmy\n/bB9IqIP2AOcUKOvF5F0kaROSZ09PT0v66DaCqLPd7U2MxtUs8PnGuBVwHxgJ/DXzSwmIq6NiAUR\nsaCjo+Nl9VXwNR8zs6qaGj4R8URE9EfEAPB1smsyAF3ArLJNZ6a2rrRc2X7YPpLagEnArhp95arN\n4WNmVlVTwyddwyl5H1CaCbcaWJZmsM0hm1hwX0TsBPZKOitdzzkfuK1sn9JMtnOBe9J1obuAsyVN\nTsN6Z6e2XBUk317HzKyKo/oyuaMh6dvAO4ATJe0gm4H2DknzgQC2Ah8DiIgNklYBjwJ9wMUR0Z+6\n+jjZzLmxwJ3pAXAdcKOkzWQTG5alvnolfQG4P213eUQc6cSHo9ZW9JmPmVk1DQufiPjgIM3X1dj+\nCuCKQdo7gXmDtO8HzqvS1wpgxREXWwf+nI+ZWXXNnnBwzCrK4WNmVo3DJyc+8zEzq87hkxNf8zEz\nq87hk5OiZ7uZmVXl8MmJb69jZladwycnbYWCbyxqZlaFwycnhYK/TM7MrBqHT07aCgXfWNTMrAqH\nT07a2woc7Hf4mJkNxuGTk/ZRBfY/7/AxMxuMwycn7W1FDvT1D72hmdkI5PDJyZhRBQ74zMfMbFAO\nn5xkZz4OHzOzwTh8cpJd8/Gwm5nZYBw+OSmd+YTvcmBm9iIOn5y0t2U/Wk+3NjN7MYdPTsaMKgJ4\nurWZ2SAcPjkpnfl4urWZ2Ys1LHwkrZDULWl9WdsUSWskbUrPk8vWXSpps6SNkhaXtZ8h6ZG07mpJ\nSu3tkm5J7WslzS7bZ3l6j02SljfieF8IH5/5mJm9SCPPfK4HllS0XQLcHRFzgbvTaySdCiwDTkv7\nfE1SMe1zDXAhMDc9Sn1eAOyOiFOAq4ArU19TgMuAM4GFwGXlIZeXQ8NuPvMxM6vUsPCJiH8Beiua\nlwIr0/JK4Jyy9psj4kBEbAE2AwslTQcmRsS9kU0ju6Fin1JftwKL0lnRYmBNRPRGxG5gDS8Owbob\n156Fz76DDh8zs0rNvuYzLSJ2puXHgWlpeQawvWy7HaltRlqubD9sn4joA/YAJ9To60UkXSSpU1Jn\nT0/P0R4TAONGtwGw70Dfy+rHzOxY1OzweUE6k2nqh2Ii4tqIWBARCzo6Ol5WX+Pas/B5xuFjZvYi\nzQ6fJ9JQGum5O7V3AbPKtpuZ2rrScmX7YftIagMmAbtq9JWrUvg8e9DhY2ZWqdnhsxoozT5bDtxW\n1r4szWCbQzax4L40RLdX0lnpes75FfuU+joXuCedTd0FnC1pcppocHZqy1Xpms8zB3zNx8ysUluj\n3kjSt4F3ACdK2kE2A+1LwCpJFwDbgPcDRMQGSauAR4E+4OKIKP0W/zjZzLmxwJ3pAXAdcKOkzWQT\nG5alvnolfQG4P213eURUTnyou/HtvuZjZlZNw8InIj5YZdWiKttfAVwxSHsnMG+Q9v3AeVX6WgGs\nOOJi62DsqCKSw8fMbDDNHnY7ZklifHsbT+93+JiZVXL45GjKuNH07jvY7DLMzFqOwydHk48bze5n\nHT5mZpUcPjnymY+Z2eAcPjmaMm40ux0+ZmYv4vDJ0UkTx9D99AF/rYKZWQWHT45ee9IE+gaCzd3P\nNLsUM7OW4vDJ0WtPmgDg8DEzq+DwydG0iWMA6Hn6QJMrMTNrLQ6fHE0c00Z7W4Fuh4+Z2WEcPjmS\nxNSJ7XTv3d/sUszMWorDJ2cd49t95mNmVsHhk7OpE8Y4fMzMKjh8cuZhNzOzF3P45GzqhHb27u9j\n//P+oKmZWYnDJ2dTJ3i6tZlZJYdPzjomtgPwhIfezMxe0BLhI2mrpEckrZPUmdqmSFojaVN6nly2\n/aWSNkvaKGlxWfsZqZ/Nkq6WpNTeLumW1L5W0uxGHdvUCVn4eNKBmdkhLRE+yW9FxPyIWJBeXwLc\nHRFzgbvTaySdCiwDTgOWAF+TVEz7XANcCMxNjyWp/QJgd0ScAlwFXNmA4wEODbvt2P1so97SzKzl\ntVL4VFoKrEzLK4FzytpvjogDEbEF2AwslDQdmBgR90ZEADdU7FPq61ZgUemsKG8njBvNtInt3Hz/\n9ka8nZnZsNAq4RPADyQ9IOmi1DYtInam5ceBaWl5BlD+m3xHapuRlivbD9snIvqAPcAJlUVIukhS\np6TOnp6el39UQKEgzn/LbH7Rs8/f7WNmlrRK+LwtIuYD7wYulvT28pXpTCbyLiIiro2IBRGxoKOj\no279vmHmJAAee3xv3fo0MxvOWiJ8IqIrPXcD3wMWAk+koTTSc3favAuYVbb7zNTWlZYr2w/bR1Ib\nMAnYlcexDGb2CeMA2NH7XKPe0syspTU9fCSNkzShtAycDawHVgPL02bLgdvS8mpgWZrBNodsYsF9\naYhur6Sz0vWc8yv2KfV1LnBPOptqiJMmjaEgTzowMytpa3YBZNdyvpeu/7cBN0XEP0m6H1gl6QJg\nG/B+gIjYIGkV8CjQB1wcEaXbB3wcuB4YC9yZHgDXATdK2gz0ks2Wa5hRxQLTJ41lzWPdfOq3X0Oh\n0JC5DmZmLUsNPAEYVhYsWBCdnZ116++aH/2cK//pZ9z+ybcxb8akuvVrZtZKJD1Q9pGZqpo+7DZS\n/N6bsol3a7f0NrkSM7Pmc/g0yNSJY5h83Ch+3vNMs0sxM2s6h08DzTlxHL9w+JiZOXwaaf6syTy4\n7Sl/2NTMRjyHTwOdt2AmB/sH+M6DO4be2MzsGObwaaDXTZ/IWa+awlfv3sTje/wVC2Y2cjl8GuxL\n/+EN7DvQxzfv3dbsUszMmsbh02CzTxzHW199It/7aRcDA/6MlZmNTA6fJjjn9Bl0PfUcX/qnnzW7\nFDOzpnD4NMHvvWkG7553Etf/ZCvbdu1rdjlmZg3n8GkCSXzmvafSPqrAR1d20uup12Y2wjh8muTk\n48fyvz9yBlt37eM3/+qH/Ghj99A7mZkdIxw+TfTWV5/I6k+8jVeccBwX3fAA3/upP/9jZiODw6fJ\nXjd9It/66Jm8eup4/vzWR7j94V/hO42b2bHO4dMCjj9uNCv/8M3MnTaeT9z0Uz583VrWd+1xCJnZ\nMcvf51NFvb/P50j0DwQ3rd3G/7hrI3v39/HGmZN4x2un8s5fm8rrZ0zyl9CZWcs70u/zcfhU0Yzw\nKendd5Bb7t/OXRse56EdTxEBM44fy9tOOZHf+rWpnPHKyXRMaG9KbWZmtTh8XqZmhk+53n0H+dHG\nblY/9Cse2Labp/f3AXDypDG8YebxvGHWJF514jhmTj6OV55wHBPGjGpyxWY2kjl8BiFpCfBVoAh8\nIyK+VG3bVgmfcgf7Bli3/Ske3vEUD+3Yw8M7nmLbrmcP22bc6CJTJ46hY0I7Uye0M3XCGKZOzJYn\njhnFcaOLjBld5LjRRcaOKjI2PR83uo2ih/XM7GU60vBpa0QxrUBSEfhb4F3ADuB+Sasj4tHmVnbk\nRrcVWDhnCgvnTHmhbc9zz7Nj97Ns732WrbuepXvvAbqf3k/30wdY37WH7qe7efZg/5H1XyyUhVGR\nMel5bFlQldrHjioyqlhgVFEUCwXaCqKtqPRcoFjQ4esKolgQhYIoSBREek7LhexZEsXUrrRNsWxd\noaxN5X0UBuvv0PalbQFE1teh5ez4S21mlr8REz7AQmBzRPwCQNLNwFJg2ITPYCaNHcWksZM47eRJ\nVbd55kAf3Xv388yBPp472M+zz/ez/2A/zx7s57nn+3kuPT97sJ/9z/fz7ME+nnt+ILX3se9AHz1P\nH0jrDu3TdwzfGFXKgilb1uEhxaGVle06rP1QwFHZX5XtDuVf5TaHv0d5X4fVfCS1NFIT8rwZf0I0\n+g+XvN/td994Mp9cNDfX9xhJ4TMD2F72egdwZvkGki4CLgJ4xSte0bjKcja+vY3xHePr3m9E0D8Q\n9JUe/QMVZGoqAAAGAUlEQVTpOegbGEjPh5b7B4KBCAbi0L6l5YGAgQj6I7LXA6RtD60bCBgYOLwt\n6+fQ8kDwwvtE+X5pfQTEC/VDpFfl7US8aJvS6HSktmz50E6Rfh6VfVXuT/l2FX0N9h6Utx9NLYdt\n01jNGNJvyp9DDX7TaMAbNmJC00gKnyFFxLXAtZBd82lyOS1PSkNtxWZXYmbDzUj6kGkXMKvs9czU\nZmZmDTaSwud+YK6kOZJGA8uA1U2uycxsRBoxw24R0SfpE8BdZFOtV0TEhiaXZWY2Io2Y8AGIiDuA\nO5pdh5nZSDeSht3MzKxFOHzMzKzhHD5mZtZwDh8zM2u4EXVj0ZdCUg+w7Sh3PxF4so7l1JNrOzqu\n7ei4tqPTqrUdSV2vjIiOoTpy+ORAUueR3NW1GVzb0XFtR8e1HZ1Wra2edXnYzczMGs7hY2ZmDefw\nyce1zS6gBtd2dFzb0XFtR6dVa6tbXb7mY2ZmDeczHzMzaziHj5mZNZzDp84kLZG0UdJmSZc04f1X\nSOqWtL6sbYqkNZI2pefJZesuTbVulLQ4x7pmSfqhpEclbZD0Jy1U2xhJ90l6KNX2+Vaprez9ipJ+\nKun2VqpN0lZJj0haJ6mzxWo7XtKtkn4m6TFJb2mF2iS9Nv28So+9kj7VCrWl9/rT9P/BeknfTv9/\n1L+2eOHrhf14uQ+yr2r4OfAqYDTwEHBqg2t4O/AmYH1Z218Bl6TlS4Ar0/KpqcZ2YE6qvZhTXdOB\nN6XlCcC/p/dvhdoEjE/Lo4C1wFmtUFtZjf8JuAm4vVX+m6b32wqcWNHWKrWtBD6alkcDx7dKbWU1\nFoHHgVe2Qm3ADGALMDa9XgX8QR615fqDHWkP4C3AXWWvLwUubUIdszk8fDYC09PydGDjYPWRfdfR\nWxpU423Au1qtNuA44EHgzFapjexbd+8G3smh8GmV2rby4vBpem3ApPRLVK1WW0U9ZwM/aZXayMJn\nOzCF7Ct3bk811r02D7vVV+k/XMmO1NZs0yJiZ1p+HJiWlptSr6TZwOlkZxgtUVsa1loHdANrIqJl\nagP+J/BfgYGytlapLYAfSHpA0kUtVNscoAf4+zRc+Q1J41qktnLLgG+n5abXFhFdwJeBXwI7gT0R\n8f08anP4jDCR/XnStPn1ksYD3wE+FRF7y9c1s7aI6I+I+WRnGQslzWuF2iS9F+iOiAeqbdPk/6Zv\nSz+3dwMXS3p7+com1tZGNvx8TUScDuwjGy5qhdoAkDQa+F3gHyrXNfHf22RgKVl4nwyMk/ThPGpz\n+NRXFzCr7PXM1NZsT0iaDpCeu1N7Q+uVNIoseL4VEd9tpdpKIuIp4IfAkhap7deB35W0FbgZeKek\nb7ZIbaW/lImIbuB7wMIWqW0HsCOdwQLcShZGrVBbybuBByPiifS6FWr7bWBLRPRExPPAd4G35lGb\nw6e+7gfmSpqT/qpZBqxuck2Q1bA8LS8nu95Sal8mqV3SHGAucF8eBUgScB3wWER8pcVq65B0fFoe\nS3Yt6metUFtEXBoRMyNiNtm/p3si4sOtUJukcZImlJbJrg2sb4XaIuJxYLuk16amRcCjrVBbmQ9y\naMitVEOza/slcJak49L/s4uAx3KpLe8LaiPtAbyHbCbXz4FPN+H9v002Vvs82V9/FwAnkF2w3gT8\nAJhStv2nU60bgXfnWNfbyE7VHwbWpcd7WqS2NwA/TbWtBz6b2pteW0Wd7+DQhIOm10Y2q/Oh9NhQ\n+vfeCrWl95oPdKb/rv8HmNxCtY0DdgGTytpapbbPk/3xtR64kWwmW91r8+11zMys4TzsZmZmDefw\nMTOzhnP4mJlZwzl8zMys4Rw+ZmbWcA4fMzNrOIePmZk1nMPHbJiR9FZJlze7DrOXwx8yNTOzhvOZ\nj9kwI+kfJP1Gs+swezkcPmbDzzyy+5WZDVsOH7NhRNIYYHRE7Gl2LWYvh8PHbHg5jeyrAcyGNYeP\n2fDyejzkZscAh4/Z8OLwsWOCp1qbmVnD+czHzMwazuFjZmYN5/AxM7OGc/iYmVnDOXzMzKzhHD5m\nZtZwDh8zM2u4/w8ZMEYwswWOhQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbcAAAEaCAYAAACSFRnbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucHGWd7/HPb7rnmkkmt8kASSBAAhJuAceYFS9oNER2\nV1hf6sbVJXpY2RVcdd1dBfW17KLuEfcIHvYIRxQEXDUgoqBHxMhVXSGZQIAkXDKQQBJyGTKTTO5z\n+50/6mnSaebSc+mqnp7v+/Wq11Q/VfXUry49v36equ4yd0dERKSUlCUdgIiIyEhTchMRkZKj5CYi\nIiVHyU1EREqOkpuIiJQcJTcRESk5Sm4iIlJylNxERKTkKLklyMxONrPVZrbHzD6ddDxxMrNbzOyr\nWa/Xmtm5CcWS2LpHu7FyDg90jpjZRjN7d4whFdV7qDdJx6PklqzPAw+6+3h3vy7pYJLk7qe6+0MD\nzTcS/0Ry68h33dKrET2Hw7E5YGZ7zWx7+AdemzX9r8ysKUzfamb3mtlbe6nnITNrM7PK4cYER54j\nwz0HzWySmbmZvZRTPsPM9ptZ63BjLAZJx6PklqzjgLVJB9EXM0snHcNIKrXtGWlD3D9DOocHWNef\nu3stcDbQCHw5LPM54FvAvwMNwLHAt4H35dQ9C3gb4LnTisQ8YAswyczGZ5X/O7AZeDKRqEqNu4+J\nAfgC0Qm1B3gOWBjKHZidNd8twFfD+Ebgn4GngH3ATURvqntDPb8FJvWzzlOAh4BdRP8A3pc17QGg\nGzgI7AVOyjfmMO0s4PEw7XZgWVbcfW5TeH058EJYdh3wF1nTNob1PgUcAtLAMcBPgRZgA/DpfGLM\n2ZY+481a77sHOFY/AHqAA2Gffb6/fdzH9vy4lzpeW3cex20j8E+hvt1hW6r62OZBnT/97ec8jtug\nz++hHO+ceF53Duex745YVx/7LPtY/AfwS6AurOODebzX/wX4A3AN8Mt+5vs48Ius1+uBn2S93gTM\ny46LXs7BIZwX/wD8AvgjsCCUnQ00A98HrhmB99BGRujcG2jb6Pvcy44nn/dpr+voq/4Bz4N8Zhrt\nA3ByOFGPCa9nASeG8YGS26PhpJgO7Agn11lAFdGb+8o+1lkeTtYvAhXAu8LBOTlrnoeAvxlCzBXA\nS+FNUg58AOgk/+T2wXAylwF/GU7+o7O2eTUwE6gO86wi+odRAZwAvAic11+MOdvSb7w5/zz6rZMj\n3zD57OMjtie3jsHWGeZdEfbfZOAZ4O/6OIYbyfP86W8/D3Tc+ttn5Jfc8jrefWzjQ4RzOM99d8Sx\n6GOfZY7FTKJ/hF8BFgNd9JIQe6mjGbgUeCPRedbQx3wnEP2zLQv79SVgc9a0NqCsl7iOOH+GcF7c\nFrbpO1n77sFwfB8DLhrOe2ikz73+to3+z72NRO/pfN+nr1tHf/UPNIyVbsluoBKYa2bl7r7R3V/I\nc9n/dPft7r4F+B3wmLs/4e4HgZ8RnSy9WQDUAl939w53f4DoE+iHRyDmBUQnzLfcvdPd7wRW5lkv\n7v4Td3/F3Xvc/XaiT6zzs2a5zt03ufsB4E1AvbtfFbbjReC7wJIBYszdF/nGO5hjle8+zt6egeRT\n53Vh/7USfQKf1099+Z4//e1noN/jNpzzO7M9+RzvgeS77wY6Fj83s13A74GHibrrpgCvuntXfwGE\n62/HAXe4+yqilu5f9TZv2LY9RMfv7cB9wCtm9gbgHcDv3L1noI3O2bZ8zot5REn+KeBUM3sf0fH7\nGXA68EQvywzlPT9i514/25bPuTeY92nuOoZ8bo+J5ObuzcBngX8FdpjZMjM7Js/Ft2eNH+jldS2A\nmX0kXOTea2b3En0C2ZTz5niJ6FPU6+QuP0DMxwBbPHyUyao7L2Z2UbjDbVf4J3IaMDVrlk1Z48cB\nx2TmDfN/kejTcL77Ne94B3ms8t3Hm8hfPnVuyxrfTzgH+pDX+UM/+zkzc1/HbZjnN+R5vPOoJ599\nl8+xuNDdJ7r7ce5+aUiEO4GpeVwXXAr8xt1fDa9/FMr68jBwLlFye5ioJfqOMDycR6zZBjwvwg0u\np3A4uc0DrgY+R9StW0bUask1lPf8iJ17fW1bnudevu/T161jOOf2mEhuAO7+I3fPfKpzohMKop1Y\nkzXrUUOs/4fuXhuG9wKvADPNLHsfH0vUd5zP8v3FvBWYbmaWU3dGn9tkZscRfSr7FDDF3ScCa4Ds\nurLfQJuADeGfTWYY7+7nDxBjtoHizd0X/dWZHVu++9gHeJ1tUMdtBPW7nwc6bsM8v/M+3gPIZ9/1\nt+/780ei63QX9jWDmVUDHwLeYWbbzGwbUTfemWZ2Zh+LZZLb28L4wwyc3Ia6DRB9INlP1O33dFj3\nk+7+KFFLak0frdNBvYcGaTjHPJ//AcN6T+X5P+Z1xkRyC9/FeVf41HSQ6FNL5lPEauCvzCxlZouJ\nTuqR8BjRSfx5MysP3/f4c6KLwMON+Y9E1x8+Hep+P0d2K/a3TeOITpCWsJ6PE73h+rIC2GNmXzCz\n6lDnaWb2pgFizDZQvPluN0SfPk8I40Pdx9l15BrWcRuGPvdzmN7ncRvh83ugOPpTsH3n7ruJrgl9\n28wuNLOasI73mtk3wmwXEnVjzSVqEc0jaiX9Drioj6ofBt5JdA1wc5h3MVE3aG/dg9D/+TOQs4Cn\nPLIrrPtTYVqmu7I3eb+HhmDIxzzP/wFDPi8G8T/mdcZEciPqs/068CpR03cacEWY9hmiHb0L+Ajw\n85FYobt3hHrfG9Z7PdGF4meHG3Oo+/3Ax4BWopsL7spats9tcvd1wDeJ3izbifr4/9DPdnQDf0b0\nxtsQ4vke0d1r/e3X3H3RX7x5bXfwP4Evh66TTzO0ffxaHWb2T73EOpzjNiQD7OeBjtuInd8DxTHA\nsgXdd+7+TaLuuy8TJflNRIkhs01Lge+7+8vuvi0zAP8H+EhvXZru/jzRXY+/C6/biVpVfwj7ojd9\nnj95OCKBuftDWV2oZ9FHchvke2hQhnPMyeN/wDDPi7z+x/TGjuzCldHKzG4hutPry0nHIiKStLHS\nchMRkTFEyU1EREqOuiVFRKTkqOUmIiIlRz8kOwRTp071WbNmJR2GiMiosmrVqlfdvT6OdSm5DcGs\nWbNoampKOgwRkVHFch7zU0jqlhQRkZKj5CYiIiVHyU1EREqOkpuIiJQcJTcRESk5Sm4iIlJylNxE\nRKTkKLnF6GdPbOa/Ho3tax4iImOWkluMfvHkVm5fuSnpMERESp6SW4xSZUZnd14PkRURkWGILbmZ\n2Uwze9DM1pnZWjP7TCj/VzPbYmarw3B+1jJXmFmzmT1nZudllb/RzJ4O064zMwvllWZ2eyh/zMxm\nZS2z1MzWh2FpVvnxYd7msGxFofZBuszo7tFTGERECi3OllsX8I/uPhdYAFxmZnPDtGvdfV4YfgUQ\npi0BTgUWA9ebWSrMfwPwCWBOGBaH8ouBNnefDVwLXB3qmgxcCbwZmA9caWaTwjJXh/XPBtpCHQWR\nTpUpuYmIxCC25ObuW9398TC+B3gGmN7PIhcAy9z9kLtvAJqB+WZ2NDDB3R/16GF0twEXZi1zaxi/\nE1gYWnXnAcvdvdXd24DlwOIw7V1hXsKymbpGXLrM6OxRt6SISKElcs0tdBeeBTwWiv7ezJ4ys5uz\nWlTTgey7LzaHsulhPLf8iGXcvQvYDUzpp64pwK4wb25duTFfYmZNZtbU0tIyqO3NSJcZ3d1quYmI\nFFrsyc3MaoGfAp9193aiLsYTgHnAVuCbcceUD3e/0d0b3b2xvn5ojyNKp4wudUuKiBRcrMnNzMqJ\nEtsP3f0uAHff7u7d7t4DfJfomhjAFmBm1uIzQtmWMJ5bfsQyZpYG6oCd/dS1E5gY5s2ta8SlypTc\nRETiEOfdkgbcBDzj7tdklR+dNdtfAGvC+D3AknAH5PFEN46scPetQLuZLQh1XgTcnbVM5k7IDwAP\nhOty9wGLzGxS6PZcBNwXpj0Y5iUsm6lrxKXLyujSVwFERAouzidxnwP8NfC0ma0OZV8EPmxm8wAH\nNgJ/C+Dua83sDmAd0Z2Wl7l7d1juUuAWoBq4NwwQJc8fmFkz0Ep0tyXu3mpmXwFWhvmucvfWMP4F\nYJmZfRV4ItRREGm13EREYhFbcnP33wPWy6Rf9bPM14Cv9VLeBJzWS/lB4IN91HUzcHMv5S9yuCu0\noNKpMiU3EZEY6BdKYqQvcYuIxEPJLUapkNyiS30iIlIoSm4xKk9FvbLqmhQRKSwltxilyqLd3aUv\ncouIFJSSW4zSZZmWm74OICJSSEpuMUqHbkndVCIiUlhKbjHKtNw61S0pIlJQSm4xSqei3a2Wm4hI\nYSm5xSila24iIrFQcovRazeUqFtSRKSglNxilOmW1PfcREQKS8ktRvoqgIhIPJTcYqRuSRGReCi5\nxag8He3uTj3TTUSkoJTcYlRelkluarmJiBSSkluMMj+crJabiEhhKbnFSN2SIiLxUHKLkbolRUTi\noeQWo/J05m5JtdxERApJyS1G5eFL3B1KbiIiBaXkFiN1S4qIxEPJLUbqlhQRiYeSW4zSZbpbUkQk\nDkpuMap47ZqbuiVFRApJyS1G6pYUEYmHkluM1C0pIhIPJbcYZX5+S92SIiKFpeQWIzOjPGXqlhQR\nKTAlt5ily8rULSkiUmCxJTczm2lmD5rZOjNba2afCeWTzWy5ma0PfydlLXOFmTWb2XNmdl5W+RvN\n7Okw7Tozs1BeaWa3h/LHzGxW1jJLwzrWm9nSrPLjw7zNYdmKQu6H8pTR0aXkJiJSSHG23LqAf3T3\nucAC4DIzmwtcDtzv7nOA+8NrwrQlwKnAYuB6M0uFum4APgHMCcPiUH4x0Obus4FrgatDXZOBK4E3\nA/OBK7OS6NXAtWGZtlBHwVSWp/TzWyIiBRZbcnP3re7+eBjfAzwDTAcuAG4Ns90KXBjGLwCWufsh\nd98ANAPzzexoYIK7P+ruDtyWs0ymrjuBhaFVdx6w3N1b3b0NWA4sDtPeFebNXX9BVJWXcbBTyU1E\npJASueYWugvPAh4DGtx9a5i0DWgI49OBTVmLbQ5l08N4bvkRy7h7F7AbmNJPXVOAXWHe3LpyY77E\nzJrMrKmlpWUQW3ukynSKg53dQ15eREQGFntyM7Na4KfAZ929PXtaaIkV5X3y7n6juze6e2N9ff2Q\n66kqL+OQrrmJiBRUrMnNzMqJEtsP3f2uULw9dDUS/u4I5VuAmVmLzwhlW8J4bvkRy5hZGqgDdvZT\n105gYpg3t66CqFLLTUSk4OK8W9KAm4Bn3P2arEn3AJm7F5cCd2eVLwl3QB5PdOPIitCF2W5mC0Kd\nF+Usk6nrA8ADoTV4H7DIzCaFG0kWAfeFaQ+GeXPXXxCVarmJiBRceuBZRsw5wF8DT5vZ6lD2ReDr\nwB1mdjHwEvAhAHdfa2Z3AOuI7rS8zN0zTZ5LgVuAauDeMECUPH9gZs1AK9Hdlrh7q5l9BVgZ5rvK\n3VvD+BeAZWb2VeCJUEfBVKVT7NrfWchViIiMebElN3f/PWB9TF7YxzJfA77WS3kTcFov5QeBD/ZR\n183Azb2Uv0j09YBYVJWrW1JEpND0CyUxq0yrW1JEpNCU3GJWWZ7S99xERApMyS1mVeVlHFK3pIhI\nQSm5xawynVK3pIhIgSm5xayqvIyO7h66e4ryu+oiIiVByS1mVeXRbz8f6lLXpIhIoSi5xawyHe3y\nQ7qpRESkYJTcYpZpuR1Uy01EpGCU3GKWabnp6wAiIoWj5BYzXXMTESk8JbeYVZWr5SYiUmhKbjGr\nTIdrbvoit4hIwSi5xSzTctMXuUVECkfJLWZquYmIFJ6SW8xe+yqAkpuISMEoucWstjJ6hN6+Q0pu\nIiKFouQWs3GVUctt36GuhCMRESldSm4xG1cRtdz2KLmJiBSMklvMysqMcRUp9h5UchMRKRQltwTU\nVqXVLSkiUkBKbgkYV5lmr5KbiEjBKLklYLySm4hIQSm5JUAtNxGRwlJyS0BtZVo3lIiIFJCSWwJq\nq9RyExEpJCW3BNSqW1JEpKCU3BJQWxl9FcDdkw5FRKQkKbklYFxlmq4e12NvREQKRMktAeOrop/g\nUtekiEhhKLklIPNkgPYDnQlHIiJSmmJLbmZ2s5ntMLM1WWX/amZbzGx1GM7PmnaFmTWb2XNmdl5W\n+RvN7Okw7Tozs1BeaWa3h/LHzGxW1jJLzWx9GJZmlR8f5m0Oy1YUej8ATKwpB2C3kpuISEHE2XK7\nBVjcS/m17j4vDL8CMLO5wBLg1LDM9WaWCvPfAHwCmBOGTJ0XA23uPhu4Frg61DUZuBJ4MzAfuNLM\nJoVlrg7rnw20hToKrq46yqFKbiIihRFbcnP3R4DWPGe/AFjm7ofcfQPQDMw3s6OBCe7+qEe3Gt4G\nXJi1zK1h/E5gYWjVnQcsd/dWd28DlgOLw7R3hXkJy2bqKqi6arXcREQKqRiuuf29mT0Vui0zLarp\nwKaseTaHsulhPLf8iGXcvQvYDUzpp64pwK4wb25dr2Nml5hZk5k1tbS0DH4rs2S6JXftV3ITESmE\npJPbDcAJwDxgK/DNZMPpm7vf6O6N7t5YX18/rLrUchMRKaxEk5u7b3f3bnfvAb5LdE0MYAswM2vW\nGaFsSxjPLT9iGTNLA3XAzn7q2glMDPPm1lVQ5akyxlWk1HITESmQRJNbuIaW8RdA5k7Ke4Al4Q7I\n44luHFnh7luBdjNbEK6ZXQTcnbVM5k7IDwAPhOty9wGLzGxS6PZcBNwXpj0Y5iUsm6mr4CbWVKjl\nJiJSIOmBZxkZZvZj4FxgqpltJrqD8Vwzmwc4sBH4WwB3X2tmdwDrgC7gMnfvDlVdSnTnZTVwbxgA\nbgJ+YGbNRDeuLAl1tZrZV4CVYb6r3D1zY8sXgGVm9lXgiVBHLCZUl7P7QEdcqxMRGVNMv284eI2N\njd7U1DSsOj5846N0dvdw5yffMkJRiYgUNzNb5e6Ncawr6RtKxqyJNeXqlhQRKRAlt4TUVZezS8lN\nRKQglNwSUhdabuoWFhEZeUpuCamrLqejq4eDnXrsjYjISFNyS8jE8PuSu3THpIjIiFNyS4ieDCAi\nUjhKbgnJ/ASXfqVERGTkKbklJNNya9unbkkRkZGm5JaQ+tpKAF7deyjhSERESo+SW0Imj6vADFr2\nquUmIjLSlNwSkk6VMammQi03EZECGPCHk83sGuCpMKx1d/03HiFTayt4dY92p4jISMvnqQDNwALg\nE8ApZraNw8luJfCIEt7QTK2tVMtNRKQABkxu7n599uvwfLXTgTOATwLfMbNPuvt9hQmxdE2trWT1\npl1JhyEiUnIG/Tw3d98AbCB6OGjmgaO/JHooqAyCWm4iIoUx7BtKwtOxfzQCsYw5U8dXsL+jm/0d\nXUmHIiJSUkbkbkl3/+ZI1DPWNIyvAmB7u1pvIiIjSV8FSNDRdVFy27b7YMKRiIiUFiW3BDVkklv7\ngYQjEREpLUpuCTpqQqblpm5JEZGRpOSWoHGVacZXpdm2Wy03EZGRpOSWsKPrqtjWrmtuIiIjSckt\nYQ0TqnRDiYjICFNyS5habiIiI0/JLWFHTaiiZc8hurp7kg5FRKRkKLkl7Ki6anocWvQzXCIiI0bJ\nLWFH1UVP5H5ll7omRURGipJbwmZMqgFgc9v+hCMRESkdSm4JmxmS26ZWJTcRkZGi5Jaw6ooU08ZX\n8rKSm4jIiIktuZnZzWa2w8zWZJVNNrPlZrY+/J2UNe0KM2s2s+fM7Lys8jea2dNh2nVmZqG80sxu\nD+WPmdmsrGWWhnWsN7OlWeXHh3mbw7IVhd4PvTl2cg0v7VRyExEZKXG23G4BFueUXQ7c7+5zgPvD\na8xsLrAEODUsc72ZpcIyNwCfAOaEIVPnxUCbu88GrgWuDnVNBq4E3gzMB67MSqJXA9eGZdpCHbE7\ndnKNuiVFREZQbMnN3R8BWnOKLwBuDeO3AhdmlS9z90Phyd/NwPzw1O8J7v6ouztwW84ymbruBBaG\nVt15wHJ3b3X3NmA5sDhMe1eYN3f9sZo5uYat7Qc51NWdxOpFREpO0tfcGsKTvAG2AQ1hfDqwKWu+\nzaFsehjPLT9iGXfvAnYDU/qpawqwK8ybW9frmNklZtZkZk0tLS2D2cYBHTu5BnfY0qYfUBYRGQlJ\nJ7fXhJaYJx1HX9z9RndvdPfG+vr6Ea37uCnRHZMvqWtSRGREJJ3ctoeuRsLfHaF8CzAza74ZoWxL\nGM8tP2IZM0sDdcDOfuraCUwM8+bWFatjQ3Lb+Oq+JFYvIlJykk5u9wCZuxeXAndnlS8Jd0AeT3Tj\nyIrQhdluZgvCNbOLcpbJ1PUB4IHQGrwPWGRmk8KNJIuA+8K0B8O8ueuPVX1tJXXV5azfsTeJ1YuI\nlJz0wLOMDDP7MXAuMNXMNhPdwfh14A4zuxh4CfgQgLuvNbM7gHVAF3CZu2futriU6M7LauDeMADc\nBPzAzJqJblxZEupqNbOvACvDfFe5e+bGli8Ay8zsq8AToY7YmRknN4zn+W17kli9iEjJsagBI4PR\n2NjoTU1NI1rnl3/+NHevfoWnrlxE+OqeiEhJMbNV7t4Yx7qS7paU4OSG8ew52MX2dj0dQERkuJTc\nisSchvEAPLddXZMiIsOl5FYkTgrJTdfdRESGT8mtSEweV0H9+Eq13ERERoCSWxE5uWE865XcRESG\nTcmtiJzUMJ7nt++lq7sn6VBEREY1JbcicsaMOg50dtPcoi9zi4gMh5JbETlz5kQAnty0K+FIRERG\nNyW3IjJrSg0TqtKs3rQ76VBEREY1JbciYmacOXMiT21Wy01EZDiU3IrMmTMm8uy2PRzo0INLRUSG\nSsmtyJx93ES6e5wnNrUlHYqIyKil5FZkGmdNpszg0RdbB55ZRER6peRWZCZUlXP69DoefWFn0qGI\niIxaSm5FaMEJU3hiU5uuu4mIDJGSWxFacOIUOrudx1/WdTcRkaFQcitCb5o1mVSZ8Ud1TYqIDImS\nWxGqrUxz5ow6HlnfknQoIiKjkpJbkVp4SgNPbd7NjvaDSYciIjLqKLkVqYWnTAPggWd3JByJiMjo\no+RWpE5uGM/0idXcr+QmIjJoSm5FysxYeMo0fr/+VQ526isBIiKDoeRWxBae0sCBzm7dNSkiMkhK\nbkVswQmTqalI8Zt125MORURkVFFyK2KV6RTvmdvAr57eyqEudU2KiORLya3Ivf/sGew+0MmDz+o7\nbyIi+VJyK3LnnDiF+vGV3PX45qRDEREZNZTcilw6VcYFZx7Dg8/toHVfR9LhiIiMCkpuo8CH3jST\nzm5n2cqXkw5FRGRUKIrkZmYbzexpM1ttZk2hbLKZLTez9eHvpKz5rzCzZjN7zszOyyp/Y6in2cyu\nMzML5ZVmdnsof8zMZmUtszSsY72ZLY1vq/N3UsN43jp7Krf990t0dvckHY6ISNEriuQWvNPd57l7\nY3h9OXC/u88B7g+vMbO5wBLgVGAxcL2ZpcIyNwCfAOaEYXEovxhoc/fZwLXA1aGuycCVwJuB+cCV\n2Um0mHz8nFlsaz/Ir9dsSzoUEZGiV0zJLdcFwK1h/FbgwqzyZe5+yN03AM3AfDM7Gpjg7o+6uwO3\n5SyTqetOYGFo1Z0HLHf3VndvA5ZzOCEWlXeePI1ZU2r4ziMvEG2eiIj0pViSmwO/NbNVZnZJKGtw\n961hfBvQEManA5uylt0cyqaH8dzyI5Zx9y5gNzCln7pex8wuMbMmM2tqaYn/tvyyMuPSd85mzZZ2\nfvuMfm9SRKQ/xZLc3uru84D3ApeZ2duzJ4aWWKLNFXe/0d0b3b2xvr4+kRjef9Z0Zk2p4Zrlz9PT\no9abiEhfiiK5ufuW8HcH8DOi61/bQ1cj4W+mubIFmJm1+IxQtiWM55YfsYyZpYE6YGc/dRWldKqM\nz7x7Ds9sbec363TtTUSkL4knNzMbZ2bjM+PAImANcA+QuXtxKXB3GL8HWBLugDye6MaRFaELs93M\nFoTraRflLJOp6wPAA6E1eB+wyMwmhRtJFoWyovW+M6cze1ot3/j1c3R06c5JEZHeJJ7ciK6l/d7M\nngRWAP/P3X8NfB14j5mtB94dXuPua4E7gHXAr4HL3D3zw4uXAt8jusnkBeDeUH4TMMXMmoHPEe68\ndPdW4CvAyjBcFcqKVqrM+NL5p/Diq/u47Y8bkw5HRKQome68G7zGxkZvampKNIaPfX8Fq15q46F/\nOpcptZWJxiIikg8zW5X1da+CKoaWmwzBl/90Lgc6uvn6vc8mHYqISNFRchulZk+r5ZK3n8BPVm3m\n4ef1xAARkWxKbqPYZ949hznTarn8p0/RfrAz6XBERIqGktsoVplO8R8fPJPt7Qf5l5+v0S+XiIgE\nSm6j3LyZE/nsu0/i56tf4YeP6akBIiKg5FYSPvXO2Zx7cj1X/WIdT7zclnQ4IiKJU3IrAWVlxrUf\nmkdDXSUX39rEiy17kw5JRCRRSm4lYtK4Cm79+HwALrp5BTvaDyYckYhIcpTcSsgJ9bV8/2NvonVf\nBx+96TF27FGCE5GxScmtxJw5cyLfu6iRTa0H+MvvPMqWXQeSDklEJHZKbiXoLbOn8oOL5/PqnkN8\n6P/+kQ2v7ks6JBGRWCm5lajGWZP58SULONDZzYXf/gO/W69fMRGRsUPJrYSdNr2On196DkdNqGLp\nzSv47iMv6iGnIjImKLmVuGOn1HDXpW9h0dyj+NqvnuFjt6zUnZQiUvKU3MaAcZVpbvjo2XzlwtNY\nsWEn533rEX69ZmvSYYmIFIyS2xhhZvz1guP45d+/jemTqvm7/3qc/3HLSjbqZhMRKUFKbmPM7Gm1\n3PXJc/jS+aewYkMri659hG/8+lk9VUBESoqS2xhUkS7jE28/gQf+8R382RlHc/1DL/C2qx/k+oea\n2d/RlXR4IiLDZnpMyuA1NjZ6U1NT0mGMmDVbdnPN8ud54NkdTK2t4GNvmcVH3nwck8ZVJB2aiJQQ\nM1vl7o2xrEvJbfBKLbllrHqpjevuX8/Dz7dQVV7GB984k4+dM4sT62uTDk1ESoCSW5Er1eSW8dy2\nPXzvdy//3z9dAAAKU0lEQVRy9+pX6OjuofG4SXyocSbnn3E0tZXppMMTkVFKya3IlXpyy9ix5yB3\nPb6FO5o28WLLPmoqUiya28D5px/N20+qp6o8lXSIIjKKKLkVubGS3DLcncdfbuMnTZv59dpt7Nrf\nybiKFAtPaeC8U4/irXOmUlddnnSYIlLklNyK3FhLbtk6u3v44ws7uXfNVu5bu53WfR2kyoyzj53I\nO06q5x0nTePUYyZQVmZJhyoiRUbJrciN5eSWrau7h9WbdvHw8y08/HwLT23eDcD4qjSNx02icdZk\n3jRrMmfMqFMXpogouRU7Jbfevbr3EL9b38KKDW2s3NhK8469AFSkyjht+gROn17HqdPrOO2YOuY0\n1FKe0tcsRcYSJbcip+SWn9Z9Hax6qY2mja08/nIb615pZ19HNxAlvDccPZ5Tj5nAnGnjmT2tljkN\ntRw1oQozdWmKlCIltyKn5DY0PT3Oxp37WPNKO2u37GbNK7tZ+0o7u/Yf/umv2so0J9aP48RptZxY\nX8vMyTXMnFTNzMk1TBlXocQnMoopuRU5JbeR4+7s3NfB+u17aW7Zyws79rJ+xx6ad+xle/uhI+at\nqUgxc1JNlPAmV3NMXTUNdVUcNSEapk2o1LU9kSIWZ3LTN3IlUWbG1NpKptZW8icnTjli2v6OLja3\nHeDlnfvZ1Lafl1v3s6n1AJta9/PfL7zK/tDFmW1iTXlIdFUcNSGqd/K4CiaPq2DSuAqmjKtgUk0F\nU2orqC5PqSUoUqKU3AAzWwz8byAFfM/dv55wSALUVKQ5qWE8JzWMf900d6f9QBfb9xxk2+6DbG/P\nDIfYFsaf29bOzr0ddPXx9PHKdFmU7MZVMKGqnPFVacaHvxPC+ITqw2WZv7WVaarKU9RUpHRTjEiR\nGvPJzcxSwLeB9wCbgZVmdo+7r0s2MumPmVFXU05dTXmvyS/D3Wk/2EXbvg527uugbV8Hrfs6aN0f\n/oay9oOdvNy6nz0Hu2g/0MmeQ/k9HaE8ZVSXp6iuSFFTkc4aT1EdEmB1RYqq8hQV6TIqU2VUpMso\nD38z45XpMipSfZSnyygzI11mpLKGTFlZ2eFp2WUiY9mYT27AfKDZ3V8EMLNlwAWAklsJMDPqqsup\nqy5n1tRxeS/X0+Ps7eg6nOwOdrHnYPR3X0cXBzq6OdDRzf7O7pzxLg50drP3UBctew5xoLOb/R3d\nHOzo5lB3Dx1dPQXc2iNlklwqK+GVGa91xRoQjRpmh18bYfprZVmvQ860I5YJS4TXZC0jpefmpW/i\n2Ck1SYcxICU3mA5synq9GXhz7kxmdglwCcCxxx4bT2SSmLIyY0JVOROqypk+sXrE6nV3unqcjq4e\nOkOyO9TVQ0f34devDd2H/3b3+GtDT6ijJ7zuel0ZdPf00N1LWaaH1nHcwYHonrLwuo/pTnhBpsyz\nph1+/doMUrIq0qOjK17JLU/ufiNwI0R3SyYcjoxSZkZ5ynStTqTA9A6DLcDMrNczQpmIiIxSSm6w\nEphjZsebWQWwBLgn4ZhERGQYxny3pLt3mdmngPuIvgpws7uvTTgsEREZhjGf3ADc/VfAr5KOQ0RE\nRoa6JUVEpOQouYmISMlRchMRkZKj5CYiIiVHj7wZAjNrAV4a4uJTgVdHMJyRpNiGRrENXrHGBYpt\nqPKJ7Th3r48jGCW3mJlZU1zPMxosxTY0im3wijUuUGxDVWyxqVtSRERKjpKbiIiUHCW3+N2YdAD9\nUGxDo9gGr1jjAsU2VEUVm665iYhIyVHLTURESo6Sm4iIlBwlt5iY2WIze87Mms3s8gTWf7OZ7TCz\nNVllk81suZmtD38nZU27IsT6nJmdV+DYZprZg2a2zszWmtlniiU+M6sysxVm9mSI7d+KJbawrpSZ\nPWFmvyymuML6NprZ02a22syaiiU+M5toZnea2bNm9oyZ/UmRxHVy2FeZod3MPlsMsYV1/UN4D6wx\nsx+H90ZRxNYrd9dQ4IHoUTovACcAFcCTwNyYY3g7cDawJqvsG8DlYfxy4OowPjfEWAkcH2JPFTC2\no4Gzw/h44PkQQ+LxAQbUhvFy4DFgQTHEFtb3OeBHwC+L6ZiGdW4EpuaUJR4fcCvwN2G8AphYDHHl\nxJgCtgHHFUNswHRgA1AdXt8BfKwYYutrUMstHvOBZnd/0d07gGXABXEG4O6PAK05xRcQvdEJfy/M\nKl/m7ofcfQPQTLQNhYptq7s/Hsb3AM8QvZkSj88je8PL8jB4McRmZjOAPwW+l1WceFwDSDQ+M6sj\n+qB3E4C7d7j7rqTj6sVC4AV3f6mIYksD1WaWBmqAV4oottdRcovHdGBT1uvNoSxpDe6+NYxvAxrC\neGLxmtks4CyiFlJRxBe6/lYDO4Dl7l4ssX0L+DzQk1VWDHFlOPBbM1tlZpcUSXzHAy3A90N37vfM\nbFwRxJVrCfDjMJ54bO6+BfhfwMvAVmC3u/+mGGLri5KbAFELheifUWLMrBb4KfBZd2/PnpZkfO7e\n7e7zgBnAfDM7LenYzOzPgB3uvqqveYrgmL417Lf3ApeZ2duzJyYUX5qoe/4Gdz8L2EfUnZZ0XK8x\nswrgfcBPcqclFVu4lnYB0YeDY4BxZvbRYoitL0pu8dgCzMx6PSOUJW27mR0NEP7uCOWxx2tm5USJ\n7YfuflexxQcQuq8eBBYXQWznAO8zs41E3dzvMrP/KoK4XhM+7ePuO4CfEXVLJR3fZmBzaH0D3EmU\n7JKOK9t7gcfdfXt4XQyxvRvY4O4t7t4J3AW8pUhi65WSWzxWAnPM7PjwqWwJcE/CMUEUw9IwvhS4\nO6t8iZlVmtnxwBxgRaGCMDMjugbyjLtfU0zxmVm9mU0M49XAe4Bnk47N3a9w9xnuPovofHrA3T+a\ndFwZZjbOzMZnxoFFwJqk43P3bcAmMzs5FC0E1iUdV44Pc7hLMhND0rG9DCwws5rwfl1IdG28GGLr\nXZx3r4zlATif6C7AF4AvJbD+HxP1lXcSfXq9GJgC3A+sB34LTM6a/0sh1ueA9xY4trcSdWc8BawO\nw/nFEB9wBvBEiG0N8C+hPPHYstZ3LofvliyKuIjuDH4yDGsz53wxxAfMA5rCMf05MKkY4grrGgfs\nBOqyyooltn8j+mC3BvgB0Z2QRRFbb4N+fktEREqOuiVFRKTkKLmJiEjJUXITEZGSo+QmIiIlR8lN\nRERKjpKbiIiUHCU3EREpOUpuIqOcmf2tmbmZnZtVdlkoe0+CoYkkRslNZPQ7neiXQN4AYGY1wN8Q\n/fr9UwnGJZIYJTeR0e8Moh9PfkN4/WmiX5Tv8cM/visypii5iYx+pxA9GfkN4Uee/xL4b6LfABQZ\nk5TcREYxM5sJ7HT3F4FpwD8D/wmcBDydZGwiSVJyExndTudwEttD9Ky5W0O5rrfJmJVOOgARGZYz\nOJzc/oOoFddtZqcTJTmRMUnJTWR0O53oCea4+y+zyucSPUdNZEzS89xERKTk6JqbiIiUHCU3EREp\nOUpuIiJScpTcRESk5Ci5iYhIyVFyExGRkqPkJiIiJef/A/GAU43PUXfPAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# figure 12.4\n", "plt.figure()\n", "plt.plot(L)\n", "plt.xlabel('$i$')\n", "plt.ylabel('$\\lambda_i$')\n", "plt.title('eigenvalue spectrum')\n", "\n", "J_cum = L.sum()-np.cumsum(L) # the minimized J[m] found by PCA is the sum of the D-M smallest eigenvalues (12.18)\n", "plt.figure()\n", "plt.plot(J_cum)\n", "plt.xlabel('$M$')\n", "plt.ylabel('$J$')\n", "plt.title('sum-of-squares distortion measure for PCA with $M$ dimensions')" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEICAYAAACQ6CLfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAECZJREFUeJzt3X+MHPV9xvH3A1yQ+Fm7KK4BuwYBUYKlOsiiETqBowAG\nB2r4A2QXiE0ppmlqgRoptUBtrBSktGqoqFQcXWSCXREokol9gBE1vwpBKsJQFw5Iwo+CbMf2BdzW\nJhhsc5/+seNwmNvvrnd3dvbu+7yk0+3NZ2fm4/U9NzM7M/tVRGBm+Tms6gbMrBoOv1mmHH6zTDn8\nZply+M0y5fCbZcrhz5SkH0r6604/t8FyZkgKSUe0uyxrn3ye37pF0gzgv4G+iNhfbTfmLX+GJB1e\ndQ9WPYd/ApH0RUlPSfpfSa9I+qNi+t2SVkhaL+k3wFeLabeOmvc7krZJ+pWkPy12z08bNf+txeM5\nkrZI+rak4WKea0ct5+uS/lPSLkmbJS3v7qtgzXL4JwhJfcCDwL8BnweWAvdI+kLxlD8GbgOOBX52\n0LwXAX8JnA+cBsxpsLrfA44HTgKuA/5Z0qSi9hvgG8DvAF8Hvinpsnb+bVYOh3/i+ApwDPD9iNgb\nEU8ADwELi/q6iHg2IkYi4sOD5r0S+HFEvBIRHwDLG6xrH/C9iNgXEeuB94EvAETEUxHxcrGel4B7\ngfM68i+0jnL4J44Tgc0RMTJq2jvUts4AmxvNO+rn1HMB3jvoDbsPqP3hQdIfSnpS0q8l/R/wZ8AJ\nzfwDrLsc/onjV8A0SaP/T6cDW4vHqdM624CTR/08rY0+fgIMAtMi4njgh4DaWJ6VxOGfOJ6jtgX+\njqQ+SXOAS4H7mpj3fuDa4g3Do4B2zukfC+yMiA8lnU3tvQbrQQ7/BBERe6mF/WLgXeBO4BsR8fMm\n5n0E+CfgSeAN4D+K0kcttPLnwPck7Qb+htofFutBvsjHPkPSF4Eh4EhfjDNxectvAEi6XNKRxSm7\nvwMedPAnNoffDrgBGAbeBD4GvlltO1Y27/abZcpbfrNMdfXWSknezTArWUQ0dV1FW1t+SRdJ+oWk\nNyQta2dZZtZdLR/zF7eF/hK4ANgCPA8sjIhXE/N4y29Wsm5s+c8G3oiIt4oLTO4D5rexPDPronbC\nfxKfvgFkC5/cRPJbkpZI2ihpYxvrMrMOK/0Nv4gYAAbAu/1mvaSdLf9WPn3318l8cgeZmfW4dsL/\nPHC6pFMkfQ5YQO1WTjMbB1re7Y+I/ZL+AngUOBy4KyJe6VhnZlaqrl7e62N+s/J15SIfMxu/HH6z\nTDn8Zply+M0y5fCbZcrhN8uUw2+WKYffLFMOv1mmHH6zTDn8Zply+M0y5fCbZcrhN8uUw2+WKYff\nLFMOv1mmHH6zTDn8Zply+M0y5fCbZcrhN8uUw2+WKYffLFMOv1mmHH6zTDn8Zply+M0y5fCbZarl\nIbqtd5xxxhl1a319fcl5zz333GT9zjvvTNZHRkaS9SqtW7eubm3BggXJeffu3dvpdnpOW+GX9Daw\nG/gY2B8RszvRlJmVrxNb/q9GxLsdWI6ZdZGP+c0y1W74A3hM0guSloz1BElLJG2UtLHNdZlZB7W7\n298fEVslfR7YIOnnEfH06CdExAAwACAp2lyfmXVIW1v+iNhafB8Gfgqc3YmmzKx8LYdf0tGSjj3w\nGLgQGOpUY2ZWLkW0ticu6VRqW3uoHT78JCJuazCPd/vHcOaZZybrixcvTtavuOKKurXDDkv/fT/x\nxBOTdUnJequ/P1VbvXp1sn7TTTcl67t27epkOx0VEen/tELLx/wR8RbwB63Ob2bV8qk+s0w5/GaZ\ncvjNMuXwm2XK4TfLVMun+lpamU/1jWlwcDBZnzdvXpc6+ayJeqqvkfPOOy9Zf/bZZ7vUyaFr9lSf\nt/xmmXL4zTLl8JtlyuE3y5TDb5Yph98sUw6/Wab80d09YMOGDcl6O+f5h4eHk/WVK1cm641uCW7n\no7vPOeecZL3RuXZrj7f8Zply+M0y5fCbZcrhN8uUw2+WKYffLFMOv1mmfD9/DzjiiPTlFlOnTm15\n2fv27UvWt2/f3vKy23Xccccl60ND6WEgGn3seMratWuT9auuuipZ/+ijj1ped9l8P7+ZJTn8Zply\n+M0y5fCbZcrhN8uUw2+WKYffLFO+n78H7N+/P1nfvHlzlzrprrlz5ybrkyZNKm3dW7ZsSdZ7+Tx+\npzTc8ku6S9KwpKFR0yZL2iDp9eJ7ef9LZlaKZnb77wYuOmjaMuDxiDgdeLz42czGkYbhj4ingZ0H\nTZ4PrCoerwIu63BfZlayVo/5p0TEtuLxdmBKvSdKWgIsaXE9ZlaStt/wi4hI3bATEQPAAPjGHrNe\n0uqpvh2SpgIU39MfEWtmPafV8A8Ci4rHi4B1nWnHzLql4f38ku4F5gAnADuA7wJrgfuB6cA7wJUR\ncfCbgmMty7v9mVmwYEHd2vXXX5+ct8zP7Z88eXKyvmvXrtLWXbZm7+dveMwfEQvrlL52SB2ZWU/x\n5b1mmXL4zTLl8JtlyuE3y5TDb5Yp39JrSY0+wnrZsvQ9XaeddlrdWl9fX0s9NWvTpk11a40+0jwH\n3vKbZcrhN8uUw2+WKYffLFMOv1mmHH6zTDn8Zpnyef4eMGPGjGT9mmuuSdbPP//8Dnbzaf39/cl6\nmUO8N7qtttE1BuvXr69b27NnT0s9TSTe8ptlyuE3y5TDb5Yph98sUw6/WaYcfrNMOfxmmWr40d0d\nXVmmH909c+bMZH1wcDBZnz59eifbOSRS+lOgy/z9efjhh5P1+fPnl7bu8azZj+72lt8sUw6/WaYc\nfrNMOfxmmXL4zTLl8JtlyuE3y5Tv5+8Bjc6lN6qX6bDD0tuHkZGR0tZ9ySWXJOsXX3xxsv7II490\nsp0Jp+GWX9JdkoYlDY2atlzSVkmbiq955bZpZp3WzG7/3cBFY0z/x4iYVXzV/8gUM+tJDcMfEU8D\nO7vQi5l1UTtv+C2V9FJxWDCp3pMkLZG0UdLGNtZlZh3WavhXAKcCs4BtwA/qPTEiBiJidkTMbnFd\nZlaClsIfETsi4uOIGAF+BJzd2bbMrGwthV/S1FE/Xg4M1XuumfWmhuf5Jd0LzAFOkLQF+C4wR9Is\nIIC3gRtK7HHcGxpK/22cM2dOsn711Vcn648++mjd2ocffpict2zXXXdd3drSpUu72IkdrGH4I2Lh\nGJNXltCLmXWRL+81y5TDb5Yph98sUw6/WaYcfrNM+aO7rVTHH3983dp7773X1rIvvfTSZD3XW3r9\n0d1mluTwm2XK4TfLlMNvlimH3yxTDr9Zphx+s0z5o7utVHPnzq26BavDW36zTDn8Zply+M0y5fCb\nZcrhN8uUw2+WKYffLFM+z9+kvr6+urULL7wwOe8TTzyRrO/Zs6elnnrBtddem6zfcccdXerEDpW3\n/GaZcvjNMuXwm2XK4TfLlMNvlimH3yxTDr9ZppoZonsasBqYQm1I7oGIuEPSZOBfgRnUhum+MiL+\np7xWy9Xf35+s33LLLXVrF1xwQXLeU045JVnfvHlzsl6myZMnJ+vz5s1L1m+//fZk/aijjjrkng5o\ndP1D1cOPj3fNbPn3A9+OiC8BXwG+JelLwDLg8Yg4HXi8+NnMxomG4Y+IbRHxYvF4N/AacBIwH1hV\nPG0VcFlZTZpZ5x3SMb+kGcCXgeeAKRGxrShtp3ZYYGbjRNPX9ks6BlgD3BQRu6RPhgOLiKg3Dp+k\nJcCSdhs1s85qassvqY9a8O+JiAeKyTskTS3qU4HhseaNiIGImB0RszvRsJl1RsPwq7aJXwm8FhGj\n39odBBYVjxcB6zrfnpmVpeEQ3ZL6gWeAl4GRYvLN1I777wemA+9QO9W3s8GyenaI7k2bNiXrM2fO\nbHnZK1asSNZ3797d8rLb1eg05VlnnZWstzPE+1NPPZWsN3rd1qxZ0/K6J7Jmh+hueMwfET8D6i3s\na4fSlJn1Dl/hZ5Yph98sUw6/WaYcfrNMOfxmmXL4zTLV8Dx/R1eW6Xn+8Wz0Zdxj2bFjR7L+4IMP\n1q3deOONyXl9y25rmj3P7y2/WaYcfrNMOfxmmXL4zTLl8JtlyuE3y5TDb5Ypn+cvzJo1K1lfunRp\n3dqiRYvq1qr25ptvJusffPBBsv7MM88k6wMDA8n60NBQsm6d5/P8Zpbk8JtlyuE3y5TDb5Yph98s\nUw6/WaYcfrNM+Tx/k4488si6tcWLFyfnvfXWW5P1SZMmJetr165N1jds2FC3tm5deiyV7du3J+s2\n/vg8v5klOfxmmXL4zTLl8JtlyuE3y5TDb5Yph98sUw3P80uaBqwGpgABDETEHZKWA9cDvy6eenNE\nrG+wrHF7nt9svGj2PH8z4Z8KTI2IFyUdC7wAXAZcCbwfEf/QbFMOv1n5mg3/EU0saBuwrXi8W9Jr\nwEnttWdmVTukY35JM4AvA88Vk5ZKeknSXZLGvEZV0hJJGyVtbKtTM+uopq/tl3QM8O/AbRHxgKQp\nwLvU3gf4W2qHBn/SYBne7TcrWceO+QEk9QEPAY9GxO1j1GcAD0VEcjRLh9+sfB27sUe1YVpXAq+N\nDn7xRuABlwP+mFazcaSZd/v7gWeAl4GRYvLNwEJgFrXd/reBG4o3B1PL8pbfrGQd3e3vFIffrHy+\nn9/Mkhx+s0w5/GaZcvjNMuXwm2XK4TfLlMNvlimH3yxTDr9Zphx+s0w5/GaZcvjNMuXwm2XK4TfL\nVMMP8Oywd4F3Rv18QjGtF/Vqb73aF7i3VnWyt99v9oldvZ//MyuXNkbE7MoaSOjV3nq1L3Bvraqq\nN+/2m2XK4TfLVNXhH6h4/Sm92luv9gXurVWV9FbpMb+ZVafqLb+ZVcThN8tUJeGXdJGkX0h6Q9Ky\nKnqoR9Lbkl6WtKnq8QWLMRCHJQ2NmjZZ0gZJrxffxxwjsaLelkvaWrx2myTNq6i3aZKelPSqpFck\n3VhMr/S1S/RVyevW9WN+SYcDvwQuALYAzwMLI+LVrjZSh6S3gdkRUfkFIZLOBd4HVh8YCk3S3wM7\nI+L7xR/OSRHxVz3S23IOcdj2knqrN6z8Yip87To53H0nVLHlPxt4IyLeioi9wH3A/Ar66HkR8TSw\n86DJ84FVxeNV1H55uq5Obz0hIrZFxIvF493AgWHlK33tEn1VoorwnwRsHvXzFip8AcYQwGOSXpC0\npOpmxjBl1LBo24EpVTYzhobDtnfTQcPK98xr18pw953mN/w+qz8iZgEXA98qdm97UtSO2XrpXO0K\n4FRqYzhuA35QZTPFsPJrgJsiYtfoWpWv3Rh9VfK6VRH+rcC0UT+fXEzrCRGxtfg+DPyU2mFKL9lx\nYITk4vtwxf38VkTsiIiPI2IE+BEVvnbFsPJrgHsi4oFicuWv3Vh9VfW6VRH+54HTJZ0i6XPAAmCw\ngj4+Q9LRxRsxSDoauJDeG3p8EFhUPF4ErKuwl0/plWHb6w0rT8WvXc8Ndx8RXf8C5lF7x/9N4JYq\neqjT16nAfxVfr1TdG3Avtd3AfdTeG7kO+F3gceB14DFgcg/19i/UhnJ/iVrQplbUWz+1XfqXgE3F\n17yqX7tEX5W8br681yxTfsPPLFMOv1mmHH6zTDn8Zply+M0y5fCbZcrhN8vU/wPobzYQwSZq8AAA\nAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEICAYAAACQ6CLfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAExFJREFUeJzt3W2MXOV1B/D/H9uLX9ev9Xqx1zY2GDAGm2plVQ1CRFFT\nQFUhUoWC1IhWqM6HJGqkfCilkcKHqkJVScSHKpJTEKZKSSMlCBRIK7BaKB8SsVDHXoJrU8sGtrte\nv2Dw+7L26Ye5oMXsnLPMnZk76/P/SSvvztk7++zYf9+ZOfd5HpoZRCSfK6oegIhUQ+EXSUrhF0lK\n4RdJSuEXSUrhF0lK4RdJSuFPgORikkby0CW3ryJ5huTxJvyMb5IcIHme5JNl709ab2bVA5C22AJg\nCMBikgvM7GRx+98BeK+olfV/AP4WwB8CmNOE+5MWU/hz2ALgvwEsA3AjgF+R/F0Avw/gvwC8X/YH\nmNnPAYBkP4BVZe9PWk9P+3O4BcAuALsBbCpuexTAXwPYWNQ+QfIXJE/U+fhFW0cuLaMzfw5bADwL\nYAWAG0n+MYArATwDYAdqzwo+YWZ/1PYRStsp/Jc5klcCuAG1s/tVAP4EwB0A/hzABtSe/b1V2QCl\nMnraf/nbBOAMgAMA9gC4HcBvzOxXqL0cGDSz8YkHkPwlyVN1Pn7Z7l9AWkNn/svfLQB2W23u9gmS\nXwQwWNS24JLX+wBgZnd+3h9CciZq/55mAJhBcjaA8Uv/Y5HOoTP/5e9TATez/zSzo8WXt2CS8Dfo\nuwDOAngQwJ8Wn3+3SfctLUAt5iGSk878Ikkp/CJJKfwiSSn8Ikm1tdV3xRVX2MyZ6i6KtMr4+Dgu\nXrzIqXxvqSSSvAPAY6j1dv/JzB5xf9jMmVi2bFmZHykijqNHj8bfVGj4aT/JGQD+EcCdqE0OuY/k\nxkbvT0Taq8xr/q0A3jazA2Y2BuAnAO5uzrBEpNXKhH8lgHcnfP1ecdunkNxWrPAycPHixRI/TkSa\nqeXv9pvZdjPrN7P+K65Qc0GkU5RJ4xCAvglfr0JzloMSkTYoE/7XAFxL8mqSXQC+CuC55gxLRFqt\n4VafmY2T/CaAf0et1feEmb3ZtJGJSEuV6vOb2QsAXmjSWESkjfQOnEhSCr9IUgq/SFIKv0hSCr9I\nUgq/SFKaXC8uckpTwzuSFqf16cwvkpTCL5KUwi+SlMIvkpTCL5KUwi+SlFp9HSBqp5Wpz5gxwz02\nWkq9la2+aFm3jz76yK1Hrbwyy8ZlaBPqzC+SlMIvkpTCL5KUwi+SlMIvkpTCL5KUwi+SlPr8U1Sm\n7xvtVNTV1eXW58+f79YXLVpUt7ZkyRL32J6eHre+dOlSt75w4UK37l1HcO7cOffYY8eOufWhIX+P\nmHfeeafh+z5//rxbv3DhglufDtcJ6MwvkpTCL5KUwi+SlMIvkpTCL5KUwi+SlMIvklSaPn80tzvq\n20Z1z6xZs9z6nDlz3LrXxweANWvW1K1dc8017rFRfdWqVW59xYoVbt0b+/j4uHvs6OioW9+7d69b\nHxgYqFsbHBx0jz106JBbP3PmjFsv8++lXUqFn+RBACcBXAAwbmb9zRiUiLReM878XzSzo024HxFp\nI73mF0mqbPgNwEskXye5bbJvILmN5ADJgTJrqolIc5V92n+rmQ2RXA7gRZJ7zeyVid9gZtsBbAeA\nrq6uzp/tIJJEqTO/mQ0Vf44CeAbA1mYMSkRar+Hwk5xHcsHHnwP4MgC/fyIiHaPM0/4eAM8U67rP\nBPAvZvZvTRlVA6L3E6Ke8tjYWMP1qKe7YMECtx6tnb948WK3vnr16rq1qE8fzdePrkGIrmG48sor\n69aitQaievS4emOL/r5Pnjzp1svO9+8EDYffzA4A2NzEsYhIG6nVJ5KUwi+SlMIvkpTCL5KUwi+S\nVEdN6Y22g/aWQy7b6ou2g/ZaN2WX5o5aedGUXu/+o8f06FF/TtbIyIhbj+7f+92iNmRvb69b7+7u\ndut9fX11a8uXL3ePjdqIR44ccetRK7ET6MwvkpTCL5KUwi+SlMIvkpTCL5KUwi+SlMIvklRH9fmr\n3NZ4xowZDR8b9fFnz57t1qPrBKJrFE6dOlW3dvDgQffYs2fPuvXjx4+79aif7fXTN2/2J4V604GB\n+PoH73Ere13I5UBnfpGkFH6RpBR+kaQUfpGkFH6RpBR+kaQUfpGkOqrPH/Hmjke98miJ6eh4T9Tn\nj5bmjpaBjubce736qJ99+vTpUvX58+e7dW8+fzRnfu7cuW49usbAu0Yh2v77cliaO6Izv0hSCr9I\nUgq/SFIKv0hSCr9IUgq/SFIKv0hS06rP74n69NH68mX6/NGx0Z4A3nz8qRzvXUcQ/d7R9Q8LFy50\n6+vWrXPrW7durVvbtGmTe2y0DsKBAwfc+uDgYN1atM7BBx984NZT9PlJPkFylOTghNuWkHyR5P7i\nT3/XCRHpOFM53T0J4I5LbnsQwE4zuxbAzuJrEZlGwvCb2SsALr1O8m4AO4rPdwC4p8njEpEWa/Q1\nf4+ZDRefjwDoqfeNJLcB2AaUWydPRJqr9Lv9Vlt1s+7Km2a23cz6zay/zJtqItJcjabxMMleACj+\n9KdIiUjHaTT8zwG4v/j8fgDPNmc4ItIu4Wt+kk8DuB3AMpLvAfgegEcA/JTkAwAOAbi3lYP8WJl1\n/aN+d/R+hPezo/uO5tRHffxo/XpvTv3SpUvdY1esWOHWoz7+xo0b3fpNN91Ut9bd3e0ee+jQIbe+\ne/dut+71+YeHh+vWgHgdg+jvdDoIw29m99UpfanJYxGRNtI7cCJJKfwiSSn8Ikkp/CJJKfwiSV02\nU3qjdlsr77/sdOGo5bVmzRq3vnr16rq1q6++2j12/fr1bn3t2rVufeXKlW7da1MeOXLEPXbfvn1u\nff/+/W59ZGSkbi2aRh1t0V3ldvLNojO/SFIKv0hSCr9IUgq/SFIKv0hSCr9IUgq/SFKXTZ8/6ru2\nchWh6L7nzJnj1r1trAGgr6/Prd9www11axs2bHCPja4hiLbRjpw4caJubWhoyD322LFjbj2aCl1G\nNGVXfX4RmbYUfpGkFH6RpBR+kaQUfpGkFH6RpBR+kaQumz5/pOx1AN6c/bLbg0dLc0e83y1agjpa\nHrvs7+Y5efJkw8cC/tbkUT06Nvq91OcXkWlL4RdJSuEXSUrhF0lK4RdJSuEXSUrhF0lqWvX5y/SU\no3511Pcts0V31BOO1pCP5r1fuHChbi3q40eiOfPR1ubeegGrVq1yj120aJFbnzt3rlv3/l6iv5Mq\n+/jtusYgPPOTfILkKMnBCbc9THKI5K7i466mjEZE2mYqT/ufBHDHJLf/wMy2FB8vNHdYItJqYfjN\n7BUAx9swFhFpozJv+H2L5O7iZUHdRehIbiM5QHIgWhdNRNqn0fD/EMA6AFsADAN4tN43mtl2M+s3\ns/5WLqIpIp9PQ2k0s8NmdsHMLgL4EYCtzR2WiLRaQ+En2Tvhy68AGKz3vSLSmcI+P8mnAdwOYBnJ\n9wB8D8DtJLcAMAAHAXy9GYOJ+ptePeo3z5o1y61HfX5PNO7o5U7U5x8bG3Pr3j730eMSjT163KI9\nB7x6d3e3e2z0HlG0H4J3jcL58+fdY71rJ4DW7gPRrmsMwn/xZnbfJDc/3oKxiEgb6R04kaQUfpGk\nFH6RpBR+kaQUfpGkptWUXq+90tXV5R47f/58tx5ND43abZ7x8fFS9XPnzrl1ryUWtfKiVmD0uCxb\ntsytL1++vOFjz54969ajVqC3bHnU6ov+TqIWaNQK7ISlv3XmF0lK4RdJSuEXSUrhF0lK4RdJSuEX\nSUrhF0mqo/r8Ue/T61lHffze3l63vmDBArfuTQ+N+vDRVtRnzpxx61HPuczyaFG/OurFX3/99W79\nxhtvrFvr6elxjz1w4IBb/+CDD9z6hx9+WLcW9fmjf4vRY15mmfl20ZlfJCmFXyQphV8kKYVfJCmF\nXyQphV8kKYVfJKmO6vNHvDnSUZ8/2g66r6/PrXt93ffff989dmRkxK2Pjo669WhZce86gNmzZ7vH\nRtc/3HzzzW79tttuc+ubNm2qW4uuX4i2F3/77bfdepk+f9ml4KcDnflFklL4RZJS+EWSUvhFklL4\nRZJS+EWSUvhFkprKFt19AJ4C0IPaltzbzewxkksA/CuAtaht032vmfkN75K83mvUd124cKFbX716\ntVtfunRp3Vq0xfbQ0FCpenT/3uMS9fGvu+46t75582a3vn79erfubXW9a9cu99hXX33Vre/bt8+t\ne+soeOszAOW3Np8OpnLmHwfwHTPbCOD3AHyD5EYADwLYaWbXAthZfC0i00QYfjMbNrM3is9PAngL\nwEoAdwPYUXzbDgD3tGqQItJ8n+s1P8m1AG4B8GsAPWY2XJRGUHtZICLTxJTDT3I+gJ8B+LaZfeqi\naasteDbpomckt5EcIDlQZq05EWmuKYWf5CzUgv9jM/t5cfNhkr1FvRfApLNTzGy7mfWbWX+0eaGI\ntE+YRtbe1nwcwFtm9v0JpecA3F98fj+AZ5s/PBFplalM6f0CgK8B2EPy497MQwAeAfBTkg8AOATg\n3rKDidonXtsoWv46WuY5ekly1VVX1a0tWbLEPTbaavrEiRNuPZr66k1nXrx4sXtstDR3NFU6ms78\n2muv1a09//zz7rEvv/yyW49apN7jHrWGo2nUl0OrLwy/mb0KoN5v+qXmDkdE2kUvwkWSUvhFklL4\nRZJS+EWSUvhFklL4RZLqqKW7o22RvX531CuPlnmO+tneFt7Rsd41AgCwbt06t16Gd20EEG8fvnfv\nXrfu9fEBYOfOnXVrAwMD7rHDw8NufWxszK17vfxoSfMMfX6d+UWSUvhFklL4RZJS+EWSUvhFklL4\nRZJS+EWS6qg+f8Sbcx/N5z98+LBbf/311926t93zO++84x67YcMGtx5tHx71nL2lvaNtrqM+/p49\ne9z6m2++6dbffffdurXTp0+7x0a99Dlz5jRcjx7TaOnuy4HO/CJJKfwiSSn8Ikkp/CJJKfwiSSn8\nIkkp/CJJMZpD30xdXV0WrRPfqKgvG/V1o7o3/7u7u9s9dt68eaV+djQn37vG4fjx4+6xx44dc+vR\n9uDRnHqvV9/V1eUeW3ZtfW+HqOgaguk6X//o0aMYGxub0uB15hdJSuEXSUrhF0lK4RdJSuEXSUrh\nF0lK4RdJKpzPT7IPwFMAegAYgO1m9hjJhwH8BYAjxbc+ZGYvtGqgkagX7q0FAPh7AgDA+fPn69a8\nuf5A3DOOrrWI6t7vFv3e0X3PnTu3VN373b0+/FRM1158p5jKYh7jAL5jZm+QXADgdZIvFrUfmNk/\ntG54ItIqYfjNbBjAcPH5SZJvAVjZ6oGJSGt9ruddJNcCuAXAr4ubvkVyN8knSC6uc8w2kgMkB6Kn\noCLSPlMOP8n5AH4G4Ntm9iGAHwJYB2ALas8MHp3sODPbbmb9ZtZf9jWeiDTPlNJIchZqwf+xmf0c\nAMzssJldMLOLAH4EYGvrhikizRaGn7W3VB8H8JaZfX/C7b0Tvu0rAAabPzwRaZWpvNv/BQBfA7CH\n5K7itocA3EdyC2rtv4MAvt6SETZJ1NKKWoVVKtPSio7NsES1TG4q7/a/CmCyf0GV9fRFpDy9AyeS\nlMIvkpTCL5KUwi+SlMIvkpTCL5LUtNqiO6t2Lq8ueejML5KUwi+SlMIvkpTCL5KUwi+SlMIvkpTC\nL5JUW7foJnkEwKEJNy0DcLRtA/h8OnVsnTouQGNrVDPHtsbMfmcq39jW8H/mh5MDZtZf2QAcnTq2\nTh0XoLE1qqqx6Wm/SFIKv0hSVYd/e8U/39OpY+vUcQEaW6MqGVulr/lFpDpVn/lFpCIKv0hSlYSf\n5B0k/4fk2yQfrGIM9ZA8SHIPyV0kByoeyxMkR0kOTrhtCckXSe4v/px0j8SKxvYwyaHisdtF8q6K\nxtZH8j9I/pbkmyT/sri90sfOGVclj1vbX/OTnAFgH4A/APAegNcA3Gdmv23rQOogeRBAv5lVfkEI\nydsAnALwlJltKm77ewDHzeyR4j/OxWb2Vx0ytocBnKp62/ZiN6neidvKA7gHwJ+hwsfOGde9qOBx\nq+LMvxXA22Z2wMzGAPwEwN0VjKPjmdkrAI5fcvPdAHYUn+9A7R9P29UZW0cws2Eze6P4/CSAj7eV\nr/Sxc8ZViSrCvxLAuxO+fg8VPgCTMAAvkXyd5LaqBzOJHjMbLj4fAdBT5WAmEW7b3k6XbCvfMY9d\nI9vdN5ve8PusW81sC4A7AXyjeHrbkaz2mq2TerVT2ra9XSbZVv4TVT52jW5332xVhH8IQN+Er1cV\nt3UEMxsq/hwF8Aw6b+vxwx/vkFz8OVrxeD7RSdu2T7atPDrgseuk7e6rCP9rAK4leTXJLgBfBfBc\nBeP4DJLzijdiQHIegC+j87Yefw7A/cXn9wN4tsKxfEqnbNteb1t5VPzYddx292bW9g8Ad6H2jv//\nAvibKsZQZ1zrAPym+Hiz6rEBeBq1p4EfofbeyAMAlgLYCWA/gJcALOmgsf0zgD0AdqMWtN6KxnYr\nak/pdwPYVXzcVfVj54yrksdNl/eKJKU3/ESSUvhFklL4RZJS+EWSUvhFklL4RZJS+EWS+n/+LmPV\nChUlHQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEICAYAAACQ6CLfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAE9NJREFUeJzt3X2MlWV6BvDrAgaEYUQ+dIp8g4giZEczIcY1BrPp+tWs\nbtqYtamhjS2bdnfTTfaPWrvNmqZpTFN34x/NJmw1YrPr7iar1WzXNkraKJoakPIxQAsjMMIwMOCI\nDPI1wN0/zqsdcc59D/Oec9535rl+yYSZc593zjNn5uI959zneR6aGUQkPeOKHoCIFEPhF0mUwi+S\nKIVfJFEKv0iiFH6RRCn8IolS+McgktNJGsmuyy6fS/I0yb4a3Ma3SW4meY7k85fVZpB8meQnJLtI\n/n7e25Pam1D0AKQu2gB0A5hOssXM+rPL/w7AoayW12EAfwvgHgCTL6v9I4DzAFqzsfwryW1mtrMG\ntys1ojP/2NQG4L8B7ARwCwCQvA3AHQDezmq5mNlLZvYvAD4cfDnJZgC/C+CvzeyUmW0E8AqAR/Pe\nptSWwj823QpgK4DtAFZklz0N4C8BLM9qnyH5a5Inqnz8+gpv+0YAF8xsz6DLtiH7T0jKQw/7x6Y2\nVM62vwXgFpJfAzAJwMsA1uOyM7+Z/U4Nb3sqgJOXXXYSQEsNb0NqQOEfY0hOAnAzKmf36wH8HoB7\nAfwRKmflcQB213EIpwBcfdll0wD0D3FdKZAe9o89KwCcBrAPwA4AqwFsM7P/QuXpQIeZXRh8AMnX\nSJ6q8vHaFd7+HgATSC4ddNmXUHn9QUpEZ/6x51YA260yV/sEybsBdGS1Nlz2fB8AzOy+K70RkhNQ\n+fsZD2A8yatQea7/CcmXAPwNyT/OxvM1VF5slBLRmX/s+VzAzew/zex49uWtGCL8I/R9AGcAPA7g\nD7LPv5/V/gyV9l8vgJ8B+FO1+cqHWsxDJE0684skSuEXSZTCL5IohV8kUQ1t9Y0bN87Gjx/fyJsU\nScrFixdx6dIlDue6ucJP8l4Az6DS6/0nM3vKu/748eMxY8aMPDcpIo6+vuHP1h7xw36S41GZunkf\nKpNFHiG5fKTfT0QaK89z/lUAOs1sn5mdB/BzAA/WZlgiUm95wj8HwMFBXx/KLvsckmuzFV82X7p0\nKcfNiUgt1f3VfjNbZ2btZtY+bpyaCyJlkSeN3QDmDfp6LmqzPJSINECe8G8CsJTkIpITAXwDwKu1\nGZaI1NuIW31mdoHktwH8Oyqtvuc0c0tk9GjorL6mpiZTn1+kfvr6+jAwMDCsN/noFTiRRCn8IolS\n+EUSpfCLJErhF0mUwi+SKC3dLa6oFVzPVjE5rI5V6b73aKEzv0iiFH6RRCn8IolS+EUSpfCLJErh\nF0mUWn2Zeras8n7vPO22ixcvusdG9Wjptahl5tWjYydM8P88o5Wh8rTzUmgF6swvkiiFXyRRCr9I\nohR+kUQp/CKJUvhFEqXwiyRqzPT5i+ylR6J+dLRt+cSJE92614uPxh2NbfLkyW796quvdustLS0j\nvu3ImTNn3PrHH39ctdbf3+8ee/78ebc+Frae05lfJFEKv0iiFH6RRCn8IolS+EUSpfCLJErhF0nU\nmOnzR/Ovo3531Lf15r1H33vq1KluPeqle73y6PtPnz7dPXbOnDluff78+W69tbXVrXu3H/X5z507\n59Z7enrc+t69e6vWOjo63GO7urrcevQeg0bufj1SucJP8gCAfgAXAVwws/ZaDEpE6q8WZ/67zex4\nDb6PiDSQnvOLJCpv+A3AGyTfI7l2qCuQXEtyM8nNY+H90CJjRd6H/XeaWTfJ6wC8TvJ/zOzNwVcw\ns3UA1gFAU1NT+V8FEUlErjO/mXVn//YCeBnAqloMSkTqb8ThJ9lMsuXTzwF8FYDfPxGR0sjzsL8V\nwMtZf30CgJ+Z2b/VZFRV5FmfPnq9YWBgwK1737+pqck99qqrrnLrs2fPdusLFixw614v/oYbbnCP\nnTdvnluP3ifQ3Nzs1r37JlqXPxL12m+66aaqtSlTprjHRu8xOHjwoFuP/p7KYMT3vpntA/ClGo5F\nRBpIrT6RRCn8IolS+EUSpfCLJErhF0lUqab05pkGGbXyLly44NajVqHXsoraXdddd51bX758ea76\nokWLqtbmzp3rHhtNJz579qxb//DDD936qVOnqtaiFunMmTPd+vXXX+/WvanQH3zwgXtsZ2enWz9y\n5IhbHw2tPp35RRKl8IskSuEXSZTCL5IohV8kUQq/SKIUfpFElarPn2f57WgZ6KgeTbv1ttGOtqmO\npsVOmzbNrUfvI/AcPnzYrZ84ccKtHzhwwK1HS1x722TPmjXLPfbuu+9269FUZ+93Fv2+o7+XvNOR\ny0BnfpFEKfwiiVL4RRKl8IskSuEXSZTCL5IohV8kUaOqWem9D8Dr6Q6nHvVtvb5vNC/9/Pnzbv3Y\nsWNuPVqrYN++fVVrH330kXtsd3e3W9+zZ49bP37c36N10qRJVWt33XWXe2z0PoCoV++NLbpfTp8+\n7daj9SFGA535RRKl8IskSuEXSZTCL5IohV8kUQq/SKIUfpFEjZk+f7QWQDQ/O8+eAZ988olbP3r0\nqFv31rYH4p/Nex9BX1+fe2w0nz/62SZOnOjWV6xYUbV2zz33uMdG+xVE75/YtWtX1dru3bvdY6P3\nXiTR5yf5HMlekh2DLptB8nWSe7N//dUqRKR0hvOw/3kA91522eMANpjZUgAbsq9FZBQJw29mbwK4\n/LHjgwDWZ5+vB/BQjcclInU20uf8rWbWk31+BEBrtSuSXAtgLRA/7xaRxsmdRqu8Ulb11TIzW2dm\n7WbWrvCLlMdI03iU5GwAyP7trd2QRKQRRhr+VwGsyT5fA+CV2gxHRBolfM5P8kUAqwHMInkIwA8A\nPAXglyQfA9AF4OF6DrIR8uwZcPHiRffYkydPuvWo1x71lM+ePTviY6M+/fz58936ypUr3foDDzxQ\ntXbHHXe4x0a2bdvm1jdu3Fi11tHRUbUGxL+T6Hc+GoThN7NHqpS+UuOxiEgD6RU4kUQp/CKJUvhF\nEqXwiyRK4RdJ1Kia0ptn2m20/HU95V02fOrUqW7dW8J65syZ7rFz585168uWLXPr3pRdAFi0aFHV\nWvT73LJli1t/55133PqmTZuq1qIpu2fOnHHrY6HVpzO/SKIUfpFEKfwiiVL4RRKl8IskSuEXSZTC\nL5KoUdXnzyPqKUdTer16tEJRc3OzW4967YsXL3brS5YsqVq78cYbc33vadOmufXoPQzedGNva3EA\nePfdd936zp073XpPT0/VWrQF98DAgFvP+/dUBjrziyRK4RdJlMIvkiiFXyRRCr9IohR+kUQp/CKJ\nUp9/mLx+drT8dTSnPpozf9ttt424HvXxW1pa3HrUr4765UeOHBnxsefOnXPr0fbh3hbeUR8/Wv9h\nNPTxIzrziyRK4RdJlMIvkiiFXyRRCr9IohR+kUQp/CKJGjN9/qiPn7fP74nm80drvHv9aADo7+93\n69689eh759n+G4j75ZMmTapai/YjiN6jsH//frfe29tbtRb9PUS/s2gdg0gZ3icQnvlJPkeyl2TH\noMueJNlNcmv2cX99hykitTach/3PA7h3iMt/ZGZt2cdvajssEam3MPxm9iaAvgaMRUQaKM8Lft8h\nuT17WjC92pVIriW5meTmIvfLE5HPG2n4fwxgMYA2AD0Anq52RTNbZ2btZtYevTAmIo0zojSa2VEz\nu2hmlwD8BMCq2g5LROptROEnOXvQl18H0FHtuiJSTmGfn+SLAFYDmEXyEIAfAFhNsg2AATgA4Ju1\nGEyeXnzePv6ECf5d4fVlo9s+efKkW+/s7HTrUZ/fW78++rmiPn/0Ok3Uq1+4cGHVWrROwbx589z6\nggUL3HpXV5db90T3S9Snz/s+gEYIw29mjwxx8bN1GIuINJBegRNJlMIvkiiFXyRRCr9IohR+kUSV\nakpv1D7xWmrRuwej1kvUEvNaXtG01mhabbSE9dGjR926N/botvO2SK+99lq37k2Nvfnmm91jZ8yY\nkeu2vSXVoym7Uasv799bGejML5IohV8kUQq/SKIUfpFEKfwiiVL4RRKl8IskqlR9/ojXO/WWiAaA\nyZMnu/Wo3+1tFx0tbx1tJR05ceKEW6/nCknR9uPRFt/XXHNN1dq0adPcY6dMmeLWo9+59zuL3psR\nvQ9gLNCZXyRRCr9IohR+kUQp/CKJUvhFEqXwiyRK4RdJ1Jjp80f95mjut9ePBoBjx45VrR08eNA9\nNurT593GzLtfonnlUR8/6sUvW7bMrbe1tVWtzZkzxz32+PHjbr27u9utnzp1qmrNew8AEPf5y7DF\ndl4684skSuEXSZTCL5IohV8kUQq/SKIUfpFEKfwiiRrOFt3zALwAoBWVLbnXmdkzJGcA+AWAhahs\n0/2wmX1Uv6H6PeuoT3/LLbe49Wi7Z2+b7F27drnHRvUjR4649ajn7N0v0dr30TbYK1eudOurV692\n6979Hq2hsGPHDre+fft2t+69TyBag6Gpqcmtp9LnvwDge2a2HMDtAL5FcjmAxwFsMLOlADZkX4vI\nKBGG38x6zGxL9nk/gN0A5gB4EMD67GrrATxUr0GKSO1d0XN+kgsB3ArgXQCtZtaTlY6g8rRAREaJ\nYb+3n+RUAL8C8F0zOzn4OY+ZGckhn8CRXAtgLVDfteZE5MoMK40km1AJ/k/N7KXs4qMkZ2f12QB6\nhzrWzNaZWbuZtSv8IuURppGVU/yzAHab2Q8HlV4FsCb7fA2AV2o/PBGpl+E87P8ygEcB7CC5Nbvs\nCQBPAfglyccAdAF4uD5D/H/e1NdoS+VoauuSJUvcemtr9Zc0br/9dvfYqGXV2dnp1qOfzWvneeMG\ngMWLF+eqz5o1y617LbW3337bPfa1115z67t373brXqsvehQaLQse/T2NhlZgGH4z2wig2k/yldoO\nR0QaRU/CRRKl8IskSuEXSZTCL5IohV8kUQq/SKJG1dLd3tTWaHns/fv3u/VoSq9Xj6a9rlq1yq1H\nS1RHW3x700+bm5vdY6OluyPR8tkbN26sWtuwYYN7bEdHh1uPpkJ7U4ajLduj+2U09PEjOvOLJErh\nF0mUwi+SKIVfJFEKv0iiFH6RRCn8IokaM31+b2ltAOjq6nLrb731llsfGBioWjt9+rR77NKlS3PV\noyWuvZ+9t3fIBZY+8/7777v1aNnxTZs2uXWvVx9tbe5tsQ3EvfapU6dWrU2Y4P/pR/P1xwKd+UUS\npfCLJErhF0mUwi+SKIVfJFEKv0iiFH6RRDHqIddSU1OTRVtGj1TUl43WYZ8yZYpb98Y9d+5c99jZ\ns2e79ZaWFrcebRft9fkPHz7sHhv12qP5+tE6Ct66/dHfXtTHj3r13tr80br9eefrFzXfv6+vDwMD\nA8O6cZ35RRKl8IskSuEXSZTCL5IohV8kUQq/SKIUfpFEhX1+kvMAvACgFYABWGdmz5B8EsCfADiW\nXfUJM/uN973q2eePRH3dqO71lKN+86VLl3LVo9+Rd7y3BkIt6nnUu5dez157Wdftv5I+/3AW87gA\n4HtmtoVkC4D3SL6e1X5kZv8w0oGKSHHC8JtZD4Ce7PN+krsBzKn3wESkvq7oOT/JhQBuBfBudtF3\nSG4n+RzJ6VWOWUtyM8nN0cNbEWmcYYef5FQAvwLwXTM7CeDHABYDaEPlkcHTQx1nZuvMrN3M2qPn\n1SLSOMNKI8kmVIL/UzN7CQDM7KiZXTSzSwB+AsDfjVJESiUMPysvaz4LYLeZ/XDQ5YOnqn0dgL+l\nqoiUynBe7f8ygEcB7CC5NbvsCQCPkGxDpf13AMA36zLCGsnbbrtw4UIthzNqRFOl87S8Gjmd/EqV\ntZVXS8N5tX8jgKHuCbenLyLlplfgRBKl8IskSuEXSZTCL5IohV8kUQq/SKJG1RbdUntF9rNT6KWX\nmc78IolS+EUSpfCLJErhF0mUwi+SKIVfJFEKv0iiGrpFN8ljALoGXTQLwPGGDeDKlHVsZR0XoLGN\nVC3HtsDMrh3OFRsa/i/cOLnZzNoLG4CjrGMr67gAjW2kihqbHvaLJErhF0lU0eFfV/Dte8o6trKO\nC9DYRqqQsRX6nF9EilP0mV9ECqLwiySqkPCTvJfk/5LsJPl4EWOohuQBkjtIbiW5ueCxPEeyl2TH\noMtmkHyd5N7s3yH3SCxobE+S7M7uu60k7y9obPNI/gfJXSR3kvzz7PJC7ztnXIXcbw1/zk9yPIA9\nAH4bwCEAmwA8Yma7GjqQKkgeANBuZoW/IYTkXQBOAXjBzFZkl/09gD4zeyr7j3O6mf1FScb2JIBT\nRW/bnu0mNXvwtvIAHgLwhyjwvnPG9TAKuN+KOPOvAtBpZvvM7DyAnwN4sIBxlJ6ZvQmg77KLHwSw\nPvt8PSp/PA1XZWylYGY9ZrYl+7wfwKfbyhd63znjKkQR4Z8D4OCgrw+hwDtgCAbgDZLvkVxb9GCG\n0GpmPdnnRwC0FjmYIYTbtjfSZdvKl+a+G8l297WmF/y+6E4zawNwH4BvZQ9vS8kqz9nK1Ksd1rbt\njTLEtvKfKfK+G+l297VWRPi7Acwb9PXc7LJSMLPu7N9eAC+jfFuPH/10h+Ts396Cx/OZMm3bPtS2\n8ijBfVem7e6LCP8mAEtJLiI5EcA3ALxawDi+gGRz9kIMSDYD+CrKt/X4qwDWZJ+vAfBKgWP5nLJs\n215tW3kUfN+Vbrt7M2v4B4D7UXnF/30Af1XEGKqMazGAbdnHzqLHBuBFVB4GDqDy2shjAGYC2ABg\nL4A3AMwo0dj+GcAOANtRCdrsgsZ2JyoP6bcD2Jp93F/0feeMq5D7TW/vFUmUXvATSZTCL5IohV8k\nUQq/SKIUfpFEKfwiiVL4RRL1fyEyyauhE8bVAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEICAYAAACQ6CLfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAE/JJREFUeJzt3X2MnWWZx/Hv1XZatC1toe20tAVEypaCMuDIbkAUYtZU\n3AhmE1Z2Y9jEpWZXjWb9Q9fNRrLZ7LKbVeMfG5O6EotxdU2UhVWEYLMsaqRapNCRQl9ogZaBdtq1\n7y2d9to/zgGHOue6hnnO28z9+ySTOXOu85zn7jPz63POuZ/7vs3dEZHyTOl0A0SkMxR+kUIp/CKF\nUvhFCqXwixRK4RcplMIvUiiFfxIys3lm5mb23Bn3LzWzo2a2vwn7eNjMjpvZ4frXMyNq55jZPWZ2\nxMyeM7M/rbo/aT6Ff3LqA3YD88xs9oj7/xHYBTzRpP18wt1n1b9+b8T9/wa8AvQCfwZ81cwua9I+\npUkU/smpD3gc+DVwGYCZXQVcA/ysXmsJM5sJ/DHwd+5+2N1/CtwLfKRV+5TxUfgnpyuBjcCTwOX1\n+74I/A2wsl57jZn9wMx+0+DrB8F+/snMhszsZ2Z2ff2+S4Bhd98y4nFPUP9PSLrHtE43QFqij9rZ\ndhFwmZl9EJgB3AOs5Ywzv7v/0Tj28VngKWov7z8M/LeZ9QGzgINnPPYgMBvpKjrzTzJmNgO4lN+e\n+fuAfwb+mtpZeQqwuep+3H29ux9y9xPuvpba24kbgcPA2Wc8fA5wqOo+pbkU/snncuAo8CywCbge\neMLdH6X2dmDA3YdHbmBmPxrxqf2ZXz8a434dMGALMM3Mlo+oXUHt8wfpInrZP/lcCTzptbHavzGz\nG4CBeq2PM97vA7j7+9/IDsxsLvD7wP8Cw8CfAO8GPuXuR8zs+8Dfm9lf1NvzQWofNkoXUfgnn9cF\n3N0fHlG7EvivJuyjB/gHYAVwCngauHnEh3x/BdwF7AH2AX/p7jrzdxnTZB4iZdJ7fpFCKfwihVL4\nRQql8IsUqq2f9k+ZMsWnTp3azl2KFOXUqVOcPn3axvLYSuE3s1XAV4CpwL+7+53R46dOncr8+fOr\n7FJEAkNDQ2N+7Lhf9pvZVGpDN99PbbDIrWa2crzPJyLtVeU9/9XANnd/1t1fAb4D3NScZolIq1UJ\n/xLghRE/76rf9zpmttrMNpjZhtOnT1fYnYg0U8s/7Xf3Ne7e7+79U6aoc0GkW1RJ425g2Yifl9bv\nE5EJoEr4fwksN7O3mNl0ahM63NecZolIq427q8/dh83sE8CD1Lr67tLILZGJo1I/v7vfD9zfpLaI\nSBvpEziRQin8IoVS+EUKpfCLFErhFymUwi9SKM3eK6FsgtdOTgBrFg9bz+ql05lfpFAKv0ihFH6R\nQin8IoVS+EUKpfCLFEpdfU3Q6u6uKl1W2dRpp06dqrR9K2dnyp67ynHPnruEbkKd+UUKpfCLFErh\nFymUwi9SKIVfpFAKv0ihFH6RQqmfvy7r1436lLO+8kzWX12lzzlbEv1Nb3pTWJ81a1al+llnndWw\nlrVteHg4rB84cCCsHz58eNzPnf1Oqv7Ou4HO/CKFUvhFCqXwixRK4RcplMIvUiiFX6RQCr9IoSZN\nP3/VMfXZ9lG/8PHjxyvtO+vvXrx4cVhftGhRw9rSpUvDbS+++OKwfv7554f1hQsXhvU5c+Y0rGV9\n7UeOHAnrW7ZsCesDAwMNa5s2bQq33bt3b1jP5jno5JTmY1Up/Ga2EzgEnAKG3b2/GY0SkdZrxpn/\nBncfasLziEgb6T2/SKGqht+BH5vZY2a2erQHmNlqM9tgZhuy90ki0j5VX/a/y913m9lC4CEze9rd\nHxn5AHdfA6wB6Onp6f5PQUQKUenM7+6769/3APcAVzejUSLSeuMOv5nNNLPZr94G3gc07lsRka5S\n5WV/L3BPfaz5NOA/3P2BprSqgajvtJVj4rPnz547G/Oe9bVfd911Yb2/v3EP63nnnRduG/XDA0yb\nFv+JzJgxI6z39PSM+7mz6x9WrFgR1i+99NKGtblz54bb/vznPw/rL7zwQlifCOP9xx1+d38WuKKJ\nbRGRNlJXn0ihFH6RQin8IoVS+EUKpfCLFGrSDOnNutuqTmEdeeWVV8J6b29vWM+68latWhXWL7nk\nkoa1aOpsgEOHDoX1bFjtvn37wnrU5ZVd7n322WeH9QULFoT1qCvw+eefD7d98cUXw3o25Dc7bt1A\nZ36RQin8IoVS+EUKpfCLFErhFymUwi9SKIVfpFATqp8/6svPhodOmRL/P5ddJxBNM51dQ5D1tZ88\neTKsDw4OhvWoT/no0aPhtjt27Ajru3fvDuvZ9NvZ7yVy+eWXh/Xs+ojZs2c3rGVDmbPrPqKhyhOF\nzvwihVL4RQql8IsUSuEXKZTCL1IohV+kUAq/SKEmVD9/NH121t+cTe2dTbV87NixcT93do3BY489\nFtazaaKjaxReeumlcNuhoXiN1ejfDfk1DNGxyZYef+tb3xrWq0yPfeLEibCe/buzv7eJQGd+kUIp\n/CKFUvhFCqXwixRK4RcplMIvUiiFX6RQE6qfP5LNAZ/Vjx8/Htajufmz8fhZn/LBgwfDetafHdWz\ntmVzEUyfPj2sV5l7f+XKleG2l112WViPxutDPBfBzp07w22z6yOy4zoRpGd+M7vLzPaY2cCI+84x\ns4fMbGv9+7zWNlNEmm0sL/u/AZy5ZMzngHXuvhxYV/9ZRCaQNPzu/giw/4y7bwLW1m+vBW5ucrtE\npMXG+56/191fnVjuJaDhYnRmthpYDfk17iLSPpXT6LWRGw1Hb7j7Gnfvd/d+hV+ke4w3jS+b2WKA\n+vc9zWuSiLTDeMN/H3Bb/fZtwL3NaY6ItEv6nt/Mvg1cD8w3s13AF4A7ge+a2UeB54BbWtnIscj6\nwrMx95moPzyb8z+T9fNH1xhk+8/mn8/G42fH9dxzzw3rN9xwQ8PaBz7wgXDb7DqAPXviF5zRPAkb\nN24Mt923b19Ynwzj+dPwu/utDUrvbXJbRKSN9AmcSKEUfpFCKfwihVL4RQql8IsUakIN6Y26tLKr\nB7Ouvqy7Lhv6WmXbrLst666Lhs1WXZr8vPPOC+tXXXVVWH/Pe97TsPb2t7893Hb//jOHlLzeL37x\ni7D+k5/8pGEtG9J7+PDhsF61a7lq93Az6MwvUiiFX6RQCr9IoRR+kUIp/CKFUvhFCqXwixRqQvXz\nR7L+7GyK6SqzDM2YMSOsL1myJKwvX748rJ9//vlhPRpWm10jMHfu3LA+Z86csJ7925YuXdqwdujQ\noXDb9evXh/VsafPt27ePe9/Z1NxZP33VIeTtoDO/SKEUfpFCKfwihVL4RQql8IsUSuEXKZTCL1Ko\nSdPPn/W7Vh3XHi1VPW9evEjxhRdeGNavvfbasJ4tVX3BBRc0rGXXIGTHZdq0+E8k2z6q79ixI9x2\n27ZtYf3pp58O60NDQw1r2Xj8qvNDTAQ684sUSuEXKZTCL1IohV+kUAq/SKEUfpFCKfwihZo0/fxZ\nv2s2nr/KvPzZNQLHjx8P69kc8UePHg3rx44da1jL+rOzpaizZbCz5cOjaxDmz58fbptdH5GNuY+O\nW5Xf91i2z457N0jP/GZ2l5ntMbOBEffdYWa7zWxj/evG1jZTRJptLC/7vwGsGuX+L7t7X/3r/uY2\nS0RaLQ2/uz8CxOsmiciEU+UDv0+a2ZP1twUNL243s9VmtsHMNmTvu0WkfcYb/q8CFwF9wCDwxUYP\ndPc17t7v7v1VJskUkeYaVxrd/WV3P+Xup4GvAVc3t1ki0mrjCr+ZLR7x44eAgUaPFZHulPbzm9m3\ngeuB+Wa2C/gCcL2Z9QEO7AQ+1ozGVFnTPPs8IXvuKv2+WV/3rl27wvrDDz8c1rNx69mY+8iBAwfC\nenYdQPZW7pprrmlYu/3228Ntr7jiirC+bNmysD4w0PiclP2+3/zmN4f1qvNDdIP0r8bdbx3l7q+3\noC0i0kb6BE6kUAq/SKEUfpFCKfwihVL4RQrVVUN6q3SPVJ1KORuCGT1/NmQ36wrMloveunVrWI/a\nHg33hbzt2fZZl1fUtlWrRhsv9lvR0uMQL/8N8bTl2TDqrKtvMtCZX6RQCr9IoRR+kUIp/CKFUvhF\nCqXwixRK4RcpVFf182ei6wCqDrHs6ekJ61GfcXaNwZEjR8J6dh3AiRMnwnrUl55Nb50ZHh4O63Pm\nzAnrixcvbljLnrtq26OhzlWGj49l+6rP3w4684sUSuEXKZTCL1IohV+kUAq/SKEUfpFCKfwihZo0\n/fxRP/xY6tFS0hCP787GxL/44oth/eDBg2E9m2sgmrY8u/4hqy9YsCCs9/X1hfVo6u5sie69e/eG\n9WxK9Og6guy6jsxEWII7ozO/SKEUfpFCKfwihVL4RQql8IsUSuEXKZTCL1KosSzRvQy4G+iltiT3\nGnf/ipmdA/wncCG1Zbpvcff/a11T437+mTNnhttm/fhve9vbwnpvb2/DWjYuPVti+5lnngnr2Rzz\n0f6z+eezufEvvvjisH7dddeF9Xe84x0Na9k1Bhs2bAjrmzdvDutVrn/IlnzP+vm7Ybx+Zixn/mHg\nM+6+EvgD4ONmthL4HLDO3ZcD6+o/i8gEkYbf3Qfd/Vf124eAzcAS4CZgbf1ha4GbW9VIEWm+N/Se\n38wuBK4E1gO97j5YL71E7W2BiEwQY76238xmAd8DPu3uB0e+p3F3N7NRJy0zs9XAasjfZ4lI+4wp\njWbWQy3433L379fvftnMFtfri4E9o23r7mvcvd/d+xV+ke6RptFqp/ivA5vd/UsjSvcBt9Vv3wbc\n2/zmiUirjOVl/7XAR4BNZraxft/ngTuB75rZR4HngFta08TfirpfommaIe/yyroC+/v7G9ZmzZoV\nbpt1ly1atCisZ0OGo6GxURclwLJly8J6dlwWLlwY1qPhyuvWrQu3/eEPfxjWsyG/UXdcNqR3InTV\nVZWG391/CjQ6Eu9tbnNEpF30JlykUAq/SKEUfpFCKfwihVL4RQql8IsUakJN3R1dIZgtc50tk51N\nrx1dY7BixYpw26p96Vl/9llnndWwlvXDz507N6xPnTo1rO/YsSOsP/DAAw1rDz74YLjttm3bwnr2\nO49k14VMhqm5MzrzixRK4RcplMIvUiiFX6RQCr9IoRR+kUIp/CKFmlD9/FFf+4kTJ8Jt9+/fH9a3\nbt0a1tevX9+wNm/evHDbiy66KKy/853vDOsHDhwI63v2jDqJEgCDg4MNawADAwNhfcuWLWH90Ucf\nDetPPfVUw1r2O3EfdWa412Rj8qNl2bPnzuqTgc78IoVS+EUKpfCLFErhFymUwi9SKIVfpFAKv0ih\nJlQ/f9T3mi2pnC1z/fzzz4f1Y8eONaxlcwFk4/2zufWzcevRmPpNmzaF2z7++ONhffv27WE9W1Mg\nmoMhmocA8jH3WT9/pIR+/IzO/CKFUvhFCqXwixRK4RcplMIvUiiFX6RQCr9IoSzr7zSzZcDdQC/g\nwBp3/4qZ3QHcDrw6qfzn3f3+6Ll6eno8Wku+laL+Zsjnp4/qWX90JvsdDA8Ph/VoLoPsGoHsuTPZ\ncY3Wuc+OeSv3HdUmsqGhIU6ePDmmf9xY/mqHgc+4+6/MbDbwmJk9VK992d3/dbwNFZHOScPv7oPA\nYP32ITPbDCxpdcNEpLXe0Ht+M7sQuBJ4dU6rT5rZk2Z2l5mNOpeVma02sw1mtiG7BFdE2mfM4Tez\nWcD3gE+7+0Hgq8BFQB+1VwZfHG07d1/j7v3u3p+9RxOR9hlTGs2sh1rwv+Xu3wdw95fd/ZS7nwa+\nBlzdumaKSLOl4bfax6JfBza7+5dG3L94xMM+BMTTwIpIVxnLp/3XAh8BNpnZxvp9nwduNbM+at1/\nO4GPtaSFTZJ93pDVT5482czmvCFVhp9Onz69Ul0mr7F82v9TYLR+w7BPX0S6mz6BEymUwi9SKIVf\npFAKv0ihFH6RQin8IoWaUFN3l2qyDj+VztKZX6RQCr9IoRR+kUIp/CKFUvhFCqXwixRK4RcpVDp1\nd1N3ZrYXeG7EXfOBobY14I3p1rZ1a7tAbRuvZrbtAndfMJYHtjX8v7Nzsw3u3t+xBgS6tW3d2i5Q\n28arU23Ty36RQin8IoXqdPjXdHj/kW5tW7e2C9S28epI2zr6nl9EOqfTZ34R6RCFX6RQHQm/ma0y\ns2fMbJuZfa4TbWjEzHaa2SYz22hmGzrclrvMbI+ZDYy47xwze8jMtta/j7pGYofadoeZ7a4fu41m\ndmOH2rbMzP7HzJ4ys1+b2afq93f02AXt6shxa/t7fjObCmwB/hDYBfwSuNXdn2prQxows51Av7t3\n/IIQM3s3cBi4290vr9/3L8B+d7+z/h/nPHf/bJe07Q7gcKeXba+vJrV45LLywM3An9PBYxe06xY6\ncNw6cea/Gtjm7s+6+yvAd4CbOtCOrufujwD7z7j7JmBt/fZaan88bdegbV3B3Qfd/Vf124eAV5eV\n7+ixC9rVEZ0I/xLghRE/76KDB2AUDvzYzB4zs9Wdbswoet19sH77JaC3k40ZRbpsezudsax81xy7\n8Sx332z6wO93vcvd+4D3Ax+vv7ztSl57z9ZNfbVjWra9XUZZVv41nTx2413uvtk6Ef7dwLIRPy+t\n39cV3H13/fse4B66b+nxl19dIbn+fU+H2/Oablq2fbRl5emCY9dNy913Ivy/BJab2VvMbDrwYeC+\nDrTjd5jZzPoHMZjZTOB9dN/S4/cBt9Vv3wbc28G2vE63LNveaFl5Onzsum65e3dv+xdwI7VP/LcD\nf9uJNjRo10XAE/WvX3e6bcC3qb0MPEnts5GPAucC64CtwI+Bc7qobd8ENgFPUgva4g617V3UXtI/\nCWysf93Y6WMXtKsjx02X94oUSh/4iRRK4RcplMIvUiiFX6RQCr9IoRR+kUIp/CKF+n9mqdKA+wvU\nBQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEICAYAAACQ6CLfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFOlJREFUeJzt3XuMnOV1x/HvWdbGie8XcI1DIQgHSpBY1waqQlJICsKo\niokQaVzUuFJaJ9S5ifyRNG0FqkibViVRIlWWnOLYISk0UQK4KCl1CIEkkgMbamwDCQnGVnB8Ib6s\nd40vrPf0j3lNl83Oc9bzzsw76+f3kVa7+555Zs7O7tl33nlu5u6ISH66qk5ARKqh4hfJlIpfJFMq\nfpFMqfhFMqXiF8mUil8kUyr+04CZzTQzN7MdI46/xcxeNbP9Je//TDO7x8x2mFm/mW0ysyXD4j8w\ns6NmNlB8/HxE+1lm9oCZHS7u48/K5CPNoeI/PfQAO4GZZjZ12PF/BF4Gnil5/93Ar4A/AqYDfwd8\nw8zOH3abj7j7lOLjohHt/w04DswFbgVWmdnbS+YkJan4Tw89wP8CzwJvBzCz3wf+EPhxEWuYux92\n9zvdfbu7D7n7w8BLwKKorZlNBm4G/t7dB9z9R8BDwJ+XyUnKU/GfHhYCm4DNwKXFsbuBvwEuKWKv\nM7OHzexgnY+Howczs7nA26j9sznpn8zsN2b2YzO7ZtjxtwGD7v7CsGPPUPyTkup0V52ANEUPtbPp\n7wBvN7P3AGcCDwDrGHHmd/c/afSBzGwC8HVgnbv/rDj8KeA5ai/t3w/8l5n1uPuLwBTg0Ii7OQRM\nRSqlM/84Z2ZnAr/H/5/5e4B/Bm6ndtbtAp5v0mN1AfdSK/KPnDzu7j9x9353P+bu66hdatxYhAeA\naSPuajrQ34ycpHEq/vHvUuBVYBuwBbgGeMbdN1K7HNjq7oPDG5jZd4e9Mz/y47ujPYiZGXAPtTft\nbnb31xI5OWDF1y8A3Wa2YFj8Mt54ySAV0Mv+8W8hsNlrc7MPmtm1wNYi1sOI630Ad18y8tgYrKL2\nCuOP3f3IyYNmNgO4EngcGAT+FHgn8PHisQ6b2beBfzCzvyzyfQ+1NyOlQir+8e8NBe7uPxgWWwg8\nWPYBzOw84EPAMWB37UUAFMf+B7gLuBg4AfwMuGnEG3x/DawB9gL7gNvcXWf+ipkW8xDJk675RTKl\n4hfJlIpfJFMqfpFMtfXd/q6uLu/q0v8bkVYZGhpiaGjI4luWLH4zuwH4InAG8O/u/rnU7bu6upg+\nfXqZhxSRhL6+vjHftuHTsJmdQW2q5hJqk0eWmdkljd6fiLRXmdfgVwC/dPdt7n4cuB9Y2py0RKTV\nyhT/fGoLPJz0cnHsDcxshZn1mlmvBhSJdI6Wv/vm7qvdfbG7Lx42LFREKlam+HcC5w77/i3FMREZ\nB8oU/1PAAjN7q5lNpLaIw/rmpCUirdZwV5+7D5rZR4BHqHX1rdFMLZHxo62z+rq7u139/CKt09fX\nx+Dg4JjeXNNwO5FMqfhFMqXiF8mUil8kUyp+kUyp+EUypdV7JenEiRPJ+NDQUMseOxoOHq0NobUj\n0vTsiGRKxS+SKRW/SKZU/CKZUvGLZErFL5IpdfWNUSu7tMrOrEx1x5155pnJtmeccUYyHnWXRd1x\nqdyibsRI1D6VW/Rz5bDknM78IplS8YtkSsUvkikVv0imVPwimVLxi2RKxS+SqWz6+VvZTz9x4sRS\nj122vzvVnx3d96RJk5LxaNfXCRMmJOOpcQb79+9Ptp01a1YyHvXV//rXv64bmz17drJtJBoHMB52\np9KZXyRTKn6RTKn4RTKl4hfJlIpfJFMqfpFMqfhFMpVNP38kmvd+/PjxurHXXnst2TaaM3/s2LFk\nPOrvnjFjRt3Ym9/85mTb6667Lhn/8Ic/nIyfc845yfjhw4frxqIxCFFfevS8rl+/vm7s7rvvTrZ9\n6aWXkvEot+7uzi+tUhma2XagHzgBDLr74mYkJSKt14x/T9e6+2+acD8i0ka65hfJVNnid+B7ZvZT\nM1sx2g3MbIWZ9ZpZbw7roomMF2Vf9l/t7jvN7Gxgg5n9zN2fGH4Dd18NrAbo7u5W9Yt0iFJnfnff\nWXzeCzwAXNGMpESk9RoufjObbGZTT34NXA9sbVZiItJaZV72zwUeKOYtdwP/4e7/3ZSs6ki9ZzA4\nOJhsW3a75lQ//+TJkxtuC3D55Zcn40uXLk3GU331UW5l58xHP9uRI0fqxt70pjcl20brJBw9ejQZ\nX7ZsWd3YwYMHk22jcQD79u1LxseDhovf3bcBlzUxFxFpI3X1iWRKxS+SKRW/SKZU/CKZUvGLZMra\nOeS2u7vbp0+f3rbHGy5aPjtaajk1RTO676gbcs2aNcn49ddfn4ynprZG042j3KOpqdHPVmZqa5Rb\n1M04ZcqUurGBgYFk2xUrRh2t/rrvf//7yXg03bhV+vr6GBwcHNO64Trzi2RKxS+SKRW/SKZU/CKZ\nUvGLZErFL5IpFb9Ipjp/feExivqTo+Wxo6mrqfEQZZeYfvLJJ5Pxa6+9NhlP5b5z585k2w0bNiTj\n0dTXqVOnJuOppbsvuyw9KTT6uefMmZOM7969u24sGteRmooM8VLv0fiHTqAzv0imVPwimVLxi2RK\nxS+SKRW/SKZU/CKZUvGLZGpczedP5RrN7Y6WiY7GCfT399eNpeaNQ7wVdSRaXrvMvPW+vr5kPOrv\njvrLU1uET5s2Ldn2/vvvT8bnzZuXjM+ePbtu7OGHH062XblyZTJ+4MCBZDz6e2sVzecXkZCKXyRT\nKn6RTKn4RTKl4hfJlIpfJFMqfpFMjav5/Kl13Muuk15mTv6hQ4eSbVN93RCPMYi2g961a1fdWPS8\nROsYRP3V0b4AqXEEN9xwQ7Lteeedl4xHuaXGKOzYsSPZNjWuA8qPK+kE4ZnfzNaY2V4z2zrs2Cwz\n22Bmvyg+z2xtmiLSbGN52b8WGPkv+tPAo+6+AHi0+F5ExpGw+N39CWD/iMNLgXXF1+uAm5qcl4i0\nWKPX/HPd/eSF5m5gbr0bmtkKYAXE15ci0j6lq9Fr75TVfbfM3Ve7+2J3XxxNAhGR9mm0+PeY2TyA\n4vPe5qUkIu3QaPGvB5YXXy8HHmpOOiLSLuE1v5ndB1wDzDGzl4E7gM8B3zCzDwI7gPe1MsmTUn3x\n0SVF1I8f9dtOnDixbixauz6676ivPHqvJDWOIJrPH40xiJ63KL5s2bK6sVtuuSXZNhqj8Oqrrzbc\n/sEHH0y2PXr0aDIePW/jQfgTuHu93967m5yLiLSR3n4XyZSKXyRTKn6RTKn4RTKl4hfJ1Ljqr0h1\nr5QdOhxtuZxafjta3jrKLdo+PFoaPNVVGE0njroZb7opPW3j9ttvT8bPOuusurHJkycn20ai7rZH\nHnmkbiy1fTekp49DPGW3bNdzO+jML5IpFb9IplT8IplS8YtkSsUvkikVv0imVPwimRpX/fwpUb9s\n1Nce9XeXaRuNIViwYEEyHvW1X3XVVXVj0bTYqD960aJFyXgkNQYi+p1ES5Z/9rOfTca3bNlSNxZN\n2Y3GVpwOS9KN/59ARBqi4hfJlIpfJFMqfpFMqfhFMqXiF8mUil8kU6dNP38kNR8f4n7bVDyam33h\nhRcm42vXrk3Gzz777GQ8Nbc86s+Ofu6y885T6wlEYxA2btyYjD/22GPJ+LZt2xrKC2D69OnJ+OHD\nh5Px8UBnfpFMqfhFMqXiF8mUil8kUyp+kUyp+EUypeIXydS46udPzdmP+qOjeevRGvCDg4N1Y9F8\n/WiN92ir6RkzZiTjr7zySt1Y1J89YcKEZDxaq6DMOID+/v5kfMmSJcn4V77ylWT8wIEDdWOpLdcB\nDh48mIyX2ecB4r/HdgjP/Ga2xsz2mtnWYcfuNLOdZrap+LixtWmKSLON5WX/WuCGUY5/wd17io/v\nNDctEWm1sPjd/QlgfxtyEZE2KvOG30fNbHNxWTCz3o3MbIWZ9ZpZbyfsTyYiNY0W/yrgAqAH2AXc\nXe+G7r7a3Re7++JOeJNDRGoaKn533+PuJ9x9CPgycEVz0xKRVmuo+M1s3rBv3wtsrXdbEelMYT+/\nmd0HXAPMMbOXgTuAa8ysB3BgO/ChZiQTvSeQmv99/PjxZNuoHz+S6reN5syn1o8H+NjHPpaML1y4\nMBl/4YUX6sb6+vqSbaPn5dChQ8n4tGnTkvGbb765buy2225Lto3MmTMnGU/9zo4dO5ZsG/0tRvFo\nH4loLYN2CCvC3ZeNcvieFuQiIm2k4b0imVLxi2RKxS+SKRW/SKZU/CKZsnYOue3u7vbUkshRLqnp\nqQMDA8m2UddLasouwOzZs1v22NH00Kg7LtXVGE0tjab0lpnqDDBzZt2R3zz11FPJttF04uXLlyfj\nTz75ZN1YtPT21KlTk/Goe7fsFPJG9fX1MTg4OKahtDrzi2RKxS+SKRW/SKZU/CKZUvGLZErFL5Ip\nFb9Ipjpq6e6obzQ1DqBsv+rkyZOT8dTU1qht1CccTUeO+upToqW7o376aFnxaFnyRYsW1Y1F01on\nTZpUKp76e4naRr+TMtPPO4XO/CKZUvGLZErFL5IpFb9IplT8IplS8YtkSsUvkqmO6uePRPO7U6Jx\nANFSzql56Zdeemmy7ebNm5PxaHnsqM84NY4gmq8f9fN3daXPD6mluQHuuuuuurFonYMoHvXVp8Z2\nROM+orEVZefzdwKd+UUypeIXyZSKXyRTKn6RTKn4RTKl4hfJlIpfJFNj2aL7XOCrwFxqW3Kvdvcv\nmtks4D+B86lt0/0+dz/QulTTfadRv2xqvwCAiy66KBlfuXJl3di73vWuZNt3vOMdyfj+/fuT8WgN\n+dSc+6ifP9pi+6qrrkrG77jjjmQ8NUYhmhMfjX/Yt29fMp4aF3LkyJFk22gdhNPBWM78g8An3f0S\n4A+AlWZ2CfBp4FF3XwA8WnwvIuNEWPzuvsvdny6+7geeB+YDS4F1xc3WATe1KkkRab5TuuY3s/OB\nhcBPgLnuvqsI7aZ2WSAi48SYi9/MpgDfAj7h7m+4GPPaxduoF3BmtsLMes2st537AopI2piK38wm\nUCv8r7v7t4vDe8xsXhGfB+wdra27r3b3xe6+eDxMdhDJRVj8VqvYe4Dn3f3zw0LrgZPbpC4HHmp+\neiLSKuEW3WZ2NfBDYAtwco7lZ6hd938D+F1gB7WuvmSfVbRFdyTVbXTgQLqXcf78+cn42rVrk/Er\nr7yybizqZvza176WjEe5R1NX+/v768ai5/vyyy9Pxi+++OJkPFq2PJXbc889l2y7atWqZPzxxx9P\nxlPTtMsurR1Nha7qVe6pbNEd9vO7+4+Aenf27lNJTEQ6h0b4iWRKxS+SKRW/SKZU/CKZUvGLZErF\nL5KpcbV0d2pqa9mllqPxDqmlnqPpobfccksyPnHixGQ86jNO/exRf3b02H19faXi9957b93Yl770\npWTbaMputD14qp8/miYdbdF9OoxW1ZlfJFMqfpFMqfhFMqXiF8mUil8kUyp+kUyp+EUyFc7nb6ay\n8/lTfdLR9t3R/Ouenp5k/NZbb20oBvE211HuZfr5d+/enWy7Z8+eZPzZZ59Nxr/5zW8m4729vXVj\n0TbZ0fMWPS+pv5fUcucQj/soux5Aq5zKfH6d+UUypeIXyZSKXyRTKn6RTKn4RTKl4hfJlIpfJFPj\naj5/ao511OcbxZ9++ulkfPv27XVjL774YrLtBz7wgWT8nHPOScY3btyYjN933311Y5s2bUq2jcYB\nRGsVRFuAT5kypW5sYGAg2bbM+AZI5x6NIYh+rtOBzvwimVLxi2RKxS+SKRW/SKZU/CKZUvGLZErF\nL5KpcD6/mZ0LfBWYCziw2t2/aGZ3An8FvFLc9DPu/p3UfZWdz19GNJ8/6vcdGhpq+LEnTZqUjEe5\nRY+d6pOO1gpo9bz0VF99NJ8/el6iny31O41+35F2roNxKk5lPv9YBvkMAp9096fNbCrwUzPbUMS+\n4O7/2miiIlKdsPjdfRewq/i638yeB+a3OjERaa1Teu1jZucDC4GfFIc+amabzWyNmc2s02aFmfWa\nWW+nvlQSydGYi9/MpgDfAj7h7oeAVcAFQA+1VwZ3j9bO3Ve7+2J3X3w67G8mcroYU/Gb2QRqhf91\nd/82gLvvcfcT7j4EfBm4onVpikizhcVvtdP1PcDz7v75YcfnDbvZe4GtzU9PRFplLF19VwM/BLYA\nJ/ucPgMso/aS34HtwIeKNwfrqrKrL9LJ70dEuaXi0aVW1OVV9nkpc6nXysvETv59l9HUrj53/xEw\n2p0l+/RFpLNphJ9IplT8IplS8YtkSsUvkikVv0imVPwimRpXS3e3UicPPa4ytyof+3Tti+8UOvOL\nZErFL5IpFb9IplT8IplS8YtkSsUvkikVv0imwvn8TX0ws1eAHcMOzQF+07YETk2n5tapeYFya1Qz\nczvP3c8ayw3bWvy/9eC1RT0XV5ZAQqfm1ql5gXJrVFW56WW/SKZU/CKZqrr4V1f8+Cmdmlun5gXK\nrVGV5FbpNb+IVKfqM7+IVETFL5KpSorfzG4ws5+b2S/N7NNV5FCPmW03sy1mtsnMeivOZY2Z7TWz\nrcOOzTKzDWb2i+LzqHskVpTbnWa2s3juNpnZjRXldq6ZPWZmz5nZs2b28eJ4pc9dIq9Knre2X/Ob\n2RnAC8B1wMvAU8Ayd3+urYnUYWbbgcXuXvmAEDN7JzAAfNXdLy2O/Quw390/V/zjnOnun+qQ3O4E\nBqretr3YTWre8G3lgZuAv6DC5y6R1/uo4Hmr4sx/BfBLd9/m7seB+4GlFeTR8dz9CWD/iMNLgXXF\n1+uo/fG0XZ3cOoK773L3p4uv+4GT28pX+twl8qpEFcU/H/jVsO9fpsInYBQOfM/MfmpmK6pOZhRz\nh22LthuYW2Uyowi3bW+nEdvKd8xz18h2982mN/x+29Xu3gMsAVYWL287kteu2Tqpr3ZM27a3yyjb\nyr+uyueu0e3um62K4t8JnDvs+7cUxzqCu+8sPu8FHqDzth7fc3KH5OLz3orzeV0nbds+2rbydMBz\n10nb3VdR/E8BC8zsrWY2EXg/sL6CPH6LmU0u3ojBzCYD19N5W4+vB5YXXy8HHqowlzfolG3b620r\nT8XPXcdtd+/ubf8AbqT2jv+LwN9WkUOdvC4Anik+nq06N+A+ai8DX6P23sgHgdnAo8AvgO8Bszoo\nt3upbeW+mVqhzasot6upvaTfDGwqPm6s+rlL5FXJ86bhvSKZ0ht+IplS8YtkSsUvkikVv0imVPwi\nmVLxi2RKxS+Sqf8D6OVszKvBqBUAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAEICAYAAACQ6CLfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEFhJREFUeJzt3X2sVHV+x/H3B2VNfAzUlCLKglGTKkmvhpjG3CibVRGi\nxe1Gi1F7sUa27ZZo201K9Y81Xbdpm+rWdKvJ3WAX2q2726CATzWItbIm3Yj2ilftrg+rBcIFlTVg\nAAX59o85d73ind8MM2fmzOX3eSWTO3O+M3O+DPdzz9Oc81NEYGb5mVR1A2ZWDYffLFMOv1mmHH6z\nTDn8Zply+M0y5fCbZcrhPwpImiIpJL1z2PTTJe2VtKuEeXx42O0TSf9Y1GZJelzSLyWNSPqupGPH\neY+zJe2X9K/t9mPtc/iPDn3ANmCKpJPGTP9rYCvwUrsziIgTR2/AbwD7gH8vyvcB7wLTi14uAf54\nnLf5J+D5dnuxcjj8R4c+4H+AV4DzACRdAFwEPFfUyvRVYCewsXg8G/hRROyPiBHgP0b7GCVpMfAB\nsKHkXqxFDv/R4XxgCNgMzCmm3Q38JXBuUfsVSY9K+qDO7dEm5jcArIpPvxv+D8DvSTpe0gxgAbU/\nAKPzOxn4K+DP2vg3Wsk+t11mE1IfsJba6vh5kn4HOA54GFjJYUv+iLiy1RlJ+iK11fqbx0x+FlgK\n7AaOKea5Zkz9W8CKiNgqqdVZW8m85J/gJB0H/CafLvn7gL+ltpQ9h9r/8WslzvJG4CcR8Yti/pOo\nLeUfAk4ATgWmFD0gqQ+4FPhOiT1YCRz+iW8OsBd4C3gZmAe8FBH/TW1zYDgiDo59gaQnxtl7P3p7\nosH8fp/akn3UVGAm8N2I+Cgi3gf+GVhY1OcBs4D/kzQCfAP4qqQXW/8nWxm82j/xnQ9sLra/P5D0\nJWC4qPVx2PY+QEQsaGVGki4CZvDpXn4i4j1JvwD+UNLdwInU9glsLp4yCPxwzNt8g9ofgz9qpQcr\nj5f8E99nAh4Rz0TEe8XD8xkn/G0YAB6KiD2HTf9dajv53gXeAA4Af1r0szciRkZvwIfA/oh4t8S+\nrAXyxTzM8uQlv1mmHH6zTDn8Zply+M0y1dVDfZK8d9GswyKiqa9RtrXkl3SFpJ9JekPS8nbey8y6\nq+VDfZKOAX4OXEbttNHngesi4tXEa7zkN+uwbiz5LwTeiIi3IuJjat/iWtTG+5lZF7UT/hnAljGP\ntxbTPkPSUkmbJG1qY15mVrKO7/CLiEFq3+/2ar9ZD2lnyb8NOGPM49OLaWY2AbQT/ueBsyXNlvQF\nYDGwrpy2zKzTWl7tj4iDkv4EeJLa1VseiIhXSuvMzDqqq2f1eZvfrPO68iUfM5u4HH6zTDn8Zply\n+M0y5fCbZcrhN8uUw2+WKYffLFMOv1mmHH6zTDn8Zply+M0y5fCbZcrhN8uUw2+WKYffLFMOv1mm\nHH6zTDn8Zply+M0y5fCbZcrhN8uUw2+WKYffLFMOv1mmHH6zTDn8Zply+M0y5fCbZarlIbqtd5xz\nzjl1a5MnT06+9uKLL07W77vvvmT90KFDyXqV1q5dW7e2ePHi5Gs//vjjstvpOW2FX9LbwB7gE+Bg\nRMwtoykz67wylvxfioj3SngfM+sib/ObZard8AfwlKQXJC0d7wmSlkraJGlTm/MysxK1u9rfHxHb\nJP06sF7S/0bEs2OfEBGDwCCApGhzfmZWkraW/BGxrfi5E3gYuLCMpsys81oOv6QTJJ00eh+4HBgu\nqzEz6yxFtLYmLulMakt7qG0+/FtEfLvBa7zaP47zzjsvWV+yZEmyfs0119StTZqU/vt+2mmnJeuS\nkvVWf3+qtmrVqmT9tttuS9Z3795dZjulioj0f1qh5W3+iHgL+K1WX29m1fKhPrNMOfxmmXL4zTLl\n8JtlyuE3y1TLh/pampkP9Y1r3bp1yfrChQu71MnnHa2H+hq55JJLkvXnnnuuS50cuWYP9XnJb5Yp\nh98sUw6/WaYcfrNMOfxmmXL4zTLl8Jtlypfu7gHr169P1ts5zr9z585kfcWKFcl6o1OC27l090UX\nXZSsNzrWbu3xkt8sUw6/WaYcfrNMOfxmmXL4zTLl8JtlyuE3y5TP5+8Bxx6b/rrF9OnTW37vAwcO\nJOsjIyMtv3e7Tj755GR9eDg9DESjy46nrFmzJlm//vrrk/WPPvqo5Xl3ms/nN7Mkh98sUw6/WaYc\nfrNMOfxmmXL4zTLl8Jtlyufz94CDBw8m61u2bOlSJ901f/78ZH3KlCkdm/fWrVuT9V4+jl+Whkt+\nSQ9I2ilpeMy0qZLWS3q9+Nm5/yUz64hmVvu/D1xx2LTlwIaIOBvYUDw2swmkYfgj4llg12GTFwEr\ni/srgatL7svMOqzVbf5pEbG9uD8CTKv3RElLgaUtzsfMOqTtHX4REakTdiJiEBgEn9hj1ktaPdS3\nQ9J0gOJn+hKxZtZzWg3/OmCguD8ArC2nHTPrlobn80t6EJgHnArsAL4JrAF+DMwE3gGujYjDdwqO\n915e7c/M4sWL69ZuueWW5Gs7ed3+qVOnJuu7d+/u2Lw7rdnz+Rtu80fEdXVKXz6ijsysp/jrvWaZ\ncvjNMuXwm2XK4TfLlMNvlimf0mtJjS5hvXx5+pyus846q25t8uTJLfXUrKGhobq1Rpc0z4GX/GaZ\ncvjNMuXwm2XK4TfLlMNvlimH3yxTDr9ZpnycvwfMmjUrWb/xxhuT9UsvvbTEbj6rv78/We/kEO+N\nTqtt9B2Dxx9/vG5t3759LfV0NPGS3yxTDr9Zphx+s0w5/GaZcvjNMuXwm2XK4TfLVMNLd5c6s0wv\n3T1nzpxkfd26dcn6zJkzy2zniEjpq0B38vfnscceS9YXLVrUsXlPZM1euttLfrNMOfxmmXL4zTLl\n8JtlyuE3y5TDb5Yph98sUz6fvwc0OpbeqN5Jkyallw+HDh3q2LyvvPLKZH3BggXJ+hNPPFFmO0ed\nhkt+SQ9I2ilpeMy0OyVtkzRU3BZ2tk0zK1szq/3fB64YZ/p3IqKvuNW/ZIqZ9aSG4Y+IZ4FdXejF\nzLqonR1+yyRtLjYLptR7kqSlkjZJ2tTGvMysZK2G/37gTKAP2A7cXe+JETEYEXMjYm6L8zKzDmgp\n/BGxIyI+iYhDwPeAC8tty8w6raXwS5o+5uFXgOF6zzWz3tTwOL+kB4F5wKmStgLfBOZJ6gMCeBv4\nWgd7nPCGh9N/G+fNm5es33DDDcn6k08+Wbe2f//+5Gs77eabb65bW7ZsWRc7scM1DH9EXDfO5BUd\n6MXMushf7zXLlMNvlimH3yxTDr9Zphx+s0z50t3WUaecckrd2vvvv9/We1911VXJeq6n9PrS3WaW\n5PCbZcrhN8uUw2+WKYffLFMOv1mmHH6zTPnS3dZR8+fPr7oFq8NLfrNMOfxmmXL4zTLl8JtlyuE3\ny5TDb5Yph98sUz7O36TJkyfXrV1++eXJ1z799NPJ+r59+1rqqRfcdNNNyfq9997bpU7sSHnJb5Yp\nh98sUw6/WaYcfrNMOfxmmXL4zTLl8Jtlqpkhus8AVgHTqA3JPRgR90qaCvwImEVtmO5rI+KXnWu1\ns/r7+5P1O+64o27tsssuS7529uzZyfqWLVuS9U6aOnVqsr5w4cJk/Z577knWjz/++CPuaVSj7z9U\nPfz4RNfMkv8g8OcRcS7w28DXJZ0LLAc2RMTZwIbisZlNEA3DHxHbI+LF4v4e4DVgBrAIWFk8bSVw\ndaeaNLPyHdE2v6RZwPnAT4FpEbG9KI1Q2ywwswmi6e/2SzoRWA3cFhG7pU+HA4uIqDcOn6SlwNJ2\nGzWzcjW15Jc0mVrwfxARDxWTd0iaXtSnAzvHe21EDEbE3IiYW0bDZlaOhuFXbRG/AngtIsbu2l0H\nDBT3B4C15bdnZp3ScIhuSf3ARuBl4FAx+XZq2/0/BmYC71A71LerwXv17BDdQ0NDyfqcOXNafu/7\n778/Wd+zZ0/L792uRocpL7jggmS9nSHen3nmmWS90ee2evXqlud9NGt2iO6G2/wR8ROg3pt9+Uia\nMrPe4W/4mWXK4TfLlMNvlimH3yxTDr9Zphx+s0w1PM5f6swyPc4/kY39Gvd4duzYkaw/8sgjdWu3\n3npr8rU+Zbc1zR7n95LfLFMOv1mmHH6zTDn8Zply+M0y5fCbZcrhN8uUj/MX+vr6kvVly5bVrQ0M\nDNStVe3NN99M1vfu3Zusb9y4MVkfHBxM1oeHh5N1K5+P85tZksNvlimH3yxTDr9Zphx+s0w5/GaZ\ncvjNMuXj/E067rjj6taWLFmSfO1dd92VrE+ZMiVZX7NmTbK+fv36urW1a9NjqYyMjCTrNvH4OL+Z\nJTn8Zply+M0y5fCbZcrhN8uUw2+WKYffLFMNj/NLOgNYBUwDAhiMiHsl3QncArxbPPX2iHi8wXtN\n2OP8ZhNFs8f5mwn/dGB6RLwo6STgBeBq4Frgw4j4+2abcvjNOq/Z8B/bxBttB7YX9/dIeg2Y0V57\nZla1I9rmlzQLOB/4aTFpmaTNkh6QNO53VCUtlbRJ0qa2OjWzUjX93X5JJwL/BXw7Ih6SNA14j9p+\ngG9R2zT4gwbv4dV+sw4rbZsfQNJk4FHgyYi4Z5z6LODRiEiOZunwm3VeaSf2qDZM6wrgtbHBL3YE\njvoK4Mu0mk0gzezt7wc2Ai8Dh4rJtwPXAX3UVvvfBr5W7BxMvZeX/GYdVupqf1kcfrPO8/n8Zpbk\n8JtlyuE3y5TDb5Yph98sUw6/WaYcfrNMOfxmmXL4zTLl8JtlyuE3y5TDb5Yph98sUw6/WaYaXsCz\nZO8B74x5fGoxrRf1am+92he4t1aV2dsXm31iV8/n/9zMpU0RMbeyBhJ6tbde7QvcW6uq6s2r/WaZ\ncvjNMlV1+Acrnn9Kr/bWq32Be2tVJb1Vus1vZtWpeslvZhVx+M0yVUn4JV0h6WeS3pC0vIoe6pH0\ntqSXJQ1VPb5gMQbiTknDY6ZNlbRe0uvFz3HHSKyotzslbSs+uyFJCyvq7QxJ/ynpVUmvSLq1mF7p\nZ5foq5LPrevb/JKOAX4OXAZsBZ4HrouIV7vaSB2S3gbmRkTlXwiRdDHwIbBqdCg0SX8H7IqIvyn+\ncE6JiL/okd7u5AiHbe9Qb/WGlV9ChZ9dmcPdl6GKJf+FwBsR8VZEfAz8EFhUQR89LyKeBXYdNnkR\nsLK4v5LaL0/X1emtJ0TE9oh4sbi/BxgdVr7Szy7RVyWqCP8MYMuYx1up8AMYRwBPSXpB0tKqmxnH\ntDHDoo0A06psZhwNh23vpsOGle+Zz66V4e7L5h1+n9cfEX3AAuDrxeptT4raNlsvHau9HziT2hiO\n24G7q2ymGFZ+NXBbROweW6vysxunr0o+tyrCvw04Y8zj04tpPSEithU/dwIPU9tM6SU7RkdILn7u\nrLifX4mIHRHxSUQcAr5HhZ9dMaz8auAHEfFQMbnyz268vqr63KoI//PA2ZJmS/oCsBhYV0EfnyPp\nhGJHDJJOAC6n94YeXwcMFPcHgLUV9vIZvTJse71h5an4s+u54e4jous3YCG1Pf5vAndU0UOdvs4E\nXipur1TdG/AgtdXAA9T2jdwM/BqwAXgdeAqY2kO9/Qu1odw3Uwva9Ip666e2Sr8ZGCpuC6v+7BJ9\nVfK5+eu9ZpnyDj+zTDn8Zply+M0y5fCbZcrhN8uUw2+WKYffLFP/D+ejSPe3Y6AaAAAAAElFTkSu\nQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# figure 12.5, PCA reconstruction; the reconstruction is perfect when M=D=784\n", "x = X[0].astype(np.float64)\n", "plt.imshow(x.reshape([28, 28]), cmap='gray', interpolation='none')\n", "plt.title('original')\n", "\n", "x -= X_\n", "for M in (1, 10, 50, 250, 784):\n", " U_M = U[:, :M]\n", " z = np.dot(U_M.T, x) # coordinates in the eigen-basis\n", " x_tilde = X_ + np.dot(U_M, z) # eq (12.20), p.566\n", " plt.figure()\n", " plt.imshow(x_tilde.reshape([28, 28]), cmap='gray', interpolation='none')\n", " plt.title('$M=%d$' %M)" ] } ], "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.5.3" } }, "nbformat": 4, "nbformat_minor": 0 }