{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as pp\n", "import sympy as sym" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsUAAAHVCAYAAAAZ9YYeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzs3Xd01fXh//HX597snUAGkIRACAl7BUGUoYAginuPumer1l2r1da2blFRcdRtrVvciGwREA17hZBAGCFAEiAhO7n38/tD68/6VQnkJu87no9zOAiEe1+Yc/TJJ59h2bYtAAAAIJA5TA8AAAAATCOKAQAAEPCIYgAAAAQ8ohgAAAABjygGAABAwCOKAQAAEPCIYgAAAAQ8ohgAAAABjygGAABAwAsy8aYdO3a0MzIyTLw1AAAAAsiyZcvKbdtOPNjHGYnijIwM5eXlmXhrAAAABBDLsra25OM4fQIAAAABjygGAABAwCOKAQAAEPCIYgAAAAQ8ohgAAAABjygGAABAwCOKAQAAEPCIYgAAAAQ8ohgAAAABjygGAABAwCOKAQAAEPCIYgAAAAQ8ohgAAAABjygGAABAwCOKAQAAEPCIYgAAAAQ8ohjtyrZt1TQ0y+22TU8BAAD4UZDpAfA/JfvrtGhTuZZu2auy6gZV1jWpqq7px++b3bYiQpzKTIxSVlKUeiRHqUdilLKSo9U1IUIOh2X6jwAAAAIMUYxWq6xt0pLN5fq6sFyLCiu0pbxGktQxKkRd4sIVEx6stPhwxYYHKzY8WNFhwdpdVa+ismotLqrQBytKfnytLnHhOm9Yus7KTVNidKipPxIAAAgwRDEOi23b+q54n15ZvEUz1+2Wy20rMsSp4d076MLhXXV0VkdlJUXJsg5+1LeqvklFe6q1cdcBfbxqpx6euVGPzy7QxL6ddMGwdB3RLaFFrwMAAHC4LNtu/3M7c3Nz7by8vHZ/34Zml0KDnO3+vv6kvsmlj1ft1CuLirW+tEqx4cE6KzdVx/VJ0cC0OAU7W3+aelFZtd74ZpveW7ZdVfXN6pkcpStGdtcZQ1KJYwAAcEgsy1pm23buQT8uUKJ41vrd+uvH6/TO1UeqS1x4u763PyivbtDLi7bozW+3a29No7KTo3XxURk6ZWAXhYe0zV806hpd+mTVTr32TbHWllRpZFZHPXB6fz5/AACgxYjinyncU61Tn16k9A4Reu/qEW0Wcv6mrtGlFxZu1rMLilTX5NK4Xsm6+KgMHdm9Q7sdtbVtW28s3ab7Pt8gh2XpzhN66ZyhaRw1BgAAB0UU/4K5+bt12at5OrF/Z009ZyBR9RtcblvvLduuKbMKtLuqQRP6JOu2iTnKTIwytmn73lrd9t5qLdlcoVE9E/XAaf3UmaPGAADgNxDFv2La/EI99MVG3T4xR9eMyTSywZvZtq35BWV64PN8bdx9QIPS4/TnSb00NCPB9DRJkttt642lW3Xf5/kKcli656Q+OmNIqulZAADAS7U0ij129wnLspyS8iSV2LZ9oqde19OuGZ2p9Tur9NDMfOWkROuYnCTTk7zGtopa/eWjtVpQUKauHSI07fzBOr5vilcdUXc4LF14ZIZG90zSre+t0i3vrtK2vbW6cVyWV+0EAAC+xZNPtLtB0gYPvl6bsCxLD58xQL07xej6N1eocE+16UnGNTa79fS8Qo1/bIHyivfqrhN6adaNozWpXyevDc30DhF64/JhOis3VVPnbNI9H6/jKXkAAOCweSSKLctKlXSCpBc88XptLTzEqed/l6uQIIeufC1PlXVNpicZ8+2WvTph6kI9PHOjjs1J0pybx+jykd0VEuT9TwAPcjr04On9deWo7nptyVb98e2VanK5Tc8CAAA+yFPl87ik2yT9apFYlnWlZVl5lmXllZWVeehtD1+XuHA9c8EQbd9Xq+vfXCFXgB1l3FfTqFvfXaWznlui2kaXXro4V89cMEQpsWGmpx0Sy7J0x/E5um1itj5etVNXvpanukaX6VkAAMDHtDqKLcs6UdIe27aX/dbH2bb9vG3bubZt5yYmJrb2bT3iiG4J+ttJfbWgoEy3vrsqIMLYtm19tLJEY6cs0PQVJbp6dKZm3TRKx+Ykm5522CzL0rVjeui+U/tpfkGZfvfS0oA++g8AAA6dJy60O0rSSZZlTZIUJinGsqx/27Z9gQdeu82dNyxde2sa9MiXBXLZth49c4CCPPBUNm+0Y1+t7vpwreZvLNPAtDg9cHo/5aTEmJ7lMecNS1dseLD++PYKnfP8N3rriuGKjQg2PQsAAPiAVkexbdt3SLpDkizLGiPpFl8J4v/6w7FZcjgsPfTFRrncth4/e6BfhbHLbevVxcV65MuNkqR7JvfW747MkNPhnRfRtcYJ/TspKixIl7/6na7+9zK9eukRPnF+NAAAMMtjt2TzddeO6SGnZen+Gfly27aeOGeQgv0gjPN3VelP76/Ryu37NSY7Uf84pa9S4yNMz2pTo3sm6oHT+uvmd1fprg/X6MHT+3vtXTQAAIB38GgU27Y9X9J8T75me7pqdKacDkv/+GyD3O4VmnruIJ89yljT0Kwn5mzSS19vUWx4sJ44Z6BOGtA5YOLw9CGp2lpRo6lzC9WtYxQPagEAAL+JI8U/c/nI7nJYlu79dL1+/5/levq8wT4VxrZta8baXbr3k/XaVVWvc49I020TchQfGWJ6Wru7cXxPbamo1YNf5KtrhwhN6tfJ9CQAAOCliOJfcOnR3eR0WLrn43U687klevKcQUrv4P2nHGwpr9HdH63Vwk3l6t0pRtMuGKzB6fGmZxnz/YNa+qtkX61ufHulOseFa2BanOlZAADAC1m23f63IcvNzbXz8vLa/X0P1Yw1pbr9/dWybem+0/pp8oDOpif9ouqGZj23oEjPLdis0CCHbj6upy4Y3tWvLhZsjYrqBp0ybZHqGt368Pcj/P6cagAA8P9ZlrXMtu3cg34cUfzbdvzwcI/l2/brnKFpumdyH4WHOE3PkiTVN7n072+2atr8Iu2tadTJAzvrzkm9lBTjWw/gaA+Few7o1GmL1Sk2TO9fM0LRYdyqDQCAQEAUe1CTy63HZxdo2vwiZSZG6anzBhm9v2+Ty61383Zo6pxN2lVVr5FZHXXLcdkawKkBv2lRYbkufHGpThnYRVPOHmh6DgAAaActjWK+vt4CwU6Hbp2Qo39fNkyVdU06+alFmjpnkypr2/epaU0utz5aWaLxUxboz9PXqHNcmN68Yrhev2wYQdwCR/XoqOuOzdIHK0r00coS03MAAIAX4UjxISqvbtCd09do5rrdigoN0vnD03XZ0d2UFN12pywUl9fore+2671lO1Re3aCclGjdOiFbx+YkBcwt1jyl2eXW2c9/o4JdB/T5DSOVlsD5xQAA+DNOn2hj63dW6ZkFRfps9U4FOR06KzdVV43K9Fhk1Te5NHPdLr357TZ9s3mvnA5Lx2Qn6dwj0nRMdpIcfvg0uvayfW+tJj2xUD1TovX2lcO5IBEAAD9GFLeT4vIaPfdVkd5fViKXbevoHh2V2zVeg9LjNSAttsUXdNU2Nmv1jkot37ZPK7bt19LNFaqqb1ZaQrjOGZquM4akKpkL6Dzmo5UluuGtlbphbJZuHN/T9BwAANBGiOJ2tquyXi8t2qJ5+Xu0aU+1JMmypOzkaA1Kj1f3jpFy27Zcti2X64fv3bb21jRq5fb9yt91QC7395+Lbh0jNTg9XqcO6qIRmR04KtxGbnp7pT5cWaJ3rjpSuRkJpucAAIA2QBQbVFnXpFXb92v5tn1avm2/Vm7bp6r65l/82OjQIPVPi9Xg9HgNSo/TwLR4JQTg0+dMOFDfpBOmfi2X29bnN4xUbDi3aQMAwN8QxV7E7bZV3disIIclh2XJ6bDktCyOAHuBFdv26Yxnl2hSv06aes5ALlwEAMDPcEs2L+JwWIoJC1ZESJDCgp0KdjoIYi8xKD1eN47L0ierdurjVTtNzwEAAIYQxQh414zpoYFpcfrbJ+u1r6bR9BwAAGAAUYyA53RYuv+0fqqqa9J9n28wPQcAABhAFAOSenWK0RWjuuvdZTu0pKjC9BwAANDOiGLgBzeMzVJ6QoTunL5G9U0u03MAAEA7IoqBH4QFO/XPU/tqc3mNps0rND0HAAC0I6IY+ImRWYk6ZWBnPbOgSJt2HzA9BwAAtBOiGPiZu07srcjQIP15+hq53e1/H28AAND+iGLgZzpGherPk3rpu+J9ejtvu+k5AACgHRDFwC84c0iqhndP0H2fb9CeA/Wm5wAAgDZGFAO/wLIs3XdqPzU0u3X/5/mm5wAAgDZGFAO/ontilK4Y2U3TV5Ro+bZ9pucAAIA2RBQDv+GaMT2UGB2qez9Zz0V3AAD4MaIY+A1RoUG6bUK2Vm7fr49WlZieAwAA2ghRDBzE6YNT1T81Vg/O2KjaxmbTcwAAQBsgioGDcDgs3X1ib+2qqtez84tMzwEAAG2AKAZaIDcjQZMHdNZzX23Wjn21pucAAAAPI4qBFvrT8TmyLOmBGdyiDQAAf0MUAy3UJS5cV47K1KerS/Vd8V7TcwAAgAcRxcAhuHp0d6XEhHGLNgAA/AxRDByCiJAg/en4HK0pqdR7y3eYngMAADyEKAYO0UkDOmtgWpwe/XKj6hpdpucAAAAPIIqBQ+RwWLrj+BztrmrQy4u3mJ4DAAA8gCgGDsOw7h00NidJz8wv0r6aRtNzAABAKxHFwGG6bWKOahqa9fS8QtNTAABAKxHFwGHKTonW6YNT9dqSrTzQAwAAH0cUA61w4/iekiVNmVVgegoAAGgFohhohc5x4bpkRIamryjRhtIq03MAAMBhIoqBVrpmTKaiQ4P00Bc8/hkAAF9FFAOtFBcRomuP6aF5G8u0pKjC9BwAAHAYiGLAAy4ekaFOsWF64It82TaPfwYAwNcQxYAHhAU7deP4nlq1fb9mrN1leg4AADhERDHgIacPTlXP5Cg9PHOjml1u03MAAMAhIIoBD3E6LN1yXLa2lNfo/eU7TM8BAACHgCgGPGh872QNSIvTE7M3qaHZZXoOAABoIaIY8CDLsnTrcdnaWVmv/yzdZnoOAABoIaIY8LCjenTQkd076Ol5happaDY9BwAAtABRDHiYZVm6ZUK2yqsb9criYtNzAABACxDFQBsY0jVe43ol6dkFRaqsbTI9BwAAHARRDLSRm4/L1oH6Zj2/sMj0FAAAcBBEMdBGenWK0eQBnfXS18UqO9Bgeg4AAPgNRDHQhm4cl6VGl1tPzys0PQUAAPwGohhoQ90To3TmkFT9Z+k2leyvMz0HAAD8CqIYaGPXj82SJE2dvcnwEgAA8GuIYqCNdY4L1wXDu+rdZdu1uaza9BwAAPALiGKgHVx7TKbCgp16jKPFAAB4JaIYaAcdo0J1yVEZ+mTVTq3fWWV6DgAA+BmiGGgnV47MVHRYkKbM2mh6CgAA+BmiGGgnsRHBumpUd83esEfLt+0zPQcAAPwEUQy0o0uO6qYOkSF69EuOFgMA4E2IYqAdRYYG6ZoxmVpUWKHFReWm5wAAgB8QxUA7u2B4V6XEhOmRmRtl27bpOQAAQEQx0O7Cgp26fmyWlm/br3kb95ieAwAARBQDRpyZm6quHSL08MwCud0cLQYAwDSiGDAg2OnQH8dlaUNplT5fW2p6DgAAAY8oBgw5aUAXZSVFacqsAjW73KbnAAAQ0IhiwBCnw9LNx/XU5rIafbCixPQcAAACWquj2LKsNMuy5lmWtcGyrHWWZd3giWFAIJjQJ0UD0uL0+KwC1Te5TM8BACBgeeJIcbOkm23b7iVpuKTfW5bV2wOvC/g9y7J0+8Rs7ays17+/2Wp6DgAAAavVUWzbdqlt28t/+OcDkjZI6tLa1wUCxYjMjhqZ1VFPzytUVX2T6TkAAAQkj55TbFlWhqRBkpb+wq9daVlWnmVZeWVlZZ58W8Dn3TYhR/tqm/TCV5tNTwEAICB5LIoty4qS9L6kP9q2XfXzX7dt+3nbtnNt285NTEz01NsCfqFfaqxO6N9JL3y9RWUHGkzPAQAg4Hgkii3LCtb3QfyGbdsfeOI1gUBzy3HZamh266m5m0xPAQAg4Hji7hOWpBclbbBte0rrJwGBqVvHSJ09NE3/+XabtlXUmp4DAEBA8cSR4qMkXSjpWMuyVv7wbZIHXhcIODeMzZLDsvTY7ALTUwAACCieuPvE17ZtW7Zt97dte+AP3z73xDgg0CTHhOmSo7rpw5Ul2lD6f07NBwAAbYQn2gFe5prRmYoODdLDMzeangIAQMAgigEvExsRrGvG9NDc/D1aurnC9BwAAAICUQx4oYtHZKhTbJj++fkGud226TkAAPg9ohjwQuEhTt1yXLZW76jUx6t2mp4DAIDfI4oBL3XqoC7q2yVGD32Rr/oml+k5AAD4NaIY8FIOh6W7TuitnZX1evHrLabnAADg14hiwIsN795B43sna9q8Qh7/DABAGyKKAS93x/E5amh280APAADaEFEMeLnuiVG6YHhXvfXtNm3cdcD0HAAA/BJRDPiAG8ZmKSo0SPd9vsH0FAAA/BJRDPiA+MgQXXdslhYUlOmrgjLTcwAA8DtEMeAjfjeiq9ITIvTPzzbIxQM9AADwKKIY8BGhQU7dPjFHG3cf0FvfbTM9BwAAv0IUAz5kUr8UDeuWoIdnbtTemkbTcwAA8BtEMeBDLMvS30/pq+r6Zj04I9/0HAAA/AZRDPiYnsnRuuzobno7b7uWbd1reg4AAH6BKAZ80PVjs9QpNkx3fbhOzS636TkAAPg8ohjwQZGhQbr7xN7aUFql15ZsNT0HAACfRxQDPmpi3xSN6pmoKbMKtKeq3vQcAAB8GlEM+CjLsnTvSX3U6HLrH5/xpDsAAFqDKAZ8WEbHSF09OlMfr9qpxYXlpucAAOCziGLAx107JlPpCRH6y0dr1djMRXcAABwOohjwcWHBTv3tpD4qKqvR818VmZ4DAIBPIooBP3BMTpJO6NdJT8zZpA2lVabnAADgc4hiwE/8/ZS+ig0P0U3vrOI0CgAADhFRDPiJhMgQ3X9aP20ordLUOZtMzwEAwKcQxYAfGd87WWcMSdW0+YVasW2f6TkAAPgMohjwM3dP7q1OseG6+Z1Vqmt0mZ4DSTUNzdpWUasNpVWqrGsyPQcA8AuCTA8A4FkxYcF66Iz+Ov+FpXpoZr7umdzH9KSA4HbbWlNSqfkby7SmpFLl1Q0qr25QRXWj6pr+9y8n8RHBSu8Qqa4JEeraIULdEyN1bHayYiOCDa0HABDFgB86qkdHXXRkV728qFjjeydrRGZH05P80r6aRn21qUzzN5bpq4IyVdQ0yrKkHolRSokNU7eOkeoQGaKO0aHqGBWq8GCnSvbXqriiVtsqarV82z59unqn3LYUEuTQhD4pOjs3TSMyO8jhsEz/8QAgoFi2bbf7m+bm5tp5eXnt/r5AIKlrdGnS1IVqbHbriz+OVHQYRyE9we22tWBTmV76eou+LiyXbX9/5HdUz0SNyU7UqKxEdYgKbfHrNTa7lb+rSh8sL9H0FSWqrGtSl7hwnZmbqjOGpCo1PqIN/zQA4P8sy1pm23buQT+OKAb817Kt+3Tms4t16qBUPXJmf1kWRx8PV32TS9NXlOjFr7eocE+1kqJDdfbQNB2bk6T+qXFyeuDIbn2TS1+u361387br6x8e2312bprumNRLseH8pQYADkdLo5jTJwA/NqRrvP5wbJamztmkgelxunB4V9OTfE7ZgQa9vqRY/166TXtrGtWnc4weO3uATujXWSFBnr1WOSzYqZMGdNZJAzprx75avbyoWC8v2qK5+Xt078l9NbFvikffDwDw/3GkGPBzbrety179Tgs3levNK4draEaC6Uk+obaxWc9/tVnPLdis+maXxuYk6/KR3TSsW0K7HnFfs6NSt7+/WutLqzSxT4ruPbmPkmLC2u39AcDXcfoEgB9V1jXplKcX6UB9sz697milxBJVv8bltvX+8h169MuN2l3VoEn9UnTLcdnqnhhlbFOTy60XFm7R47MLFBLk0J2TeunsoWmcDgMALUAUA/gfm3Yf0ClPL1JWcrTevmq4QoOcpid5nUWF5frHZxu0obRKA9PidNcJvZTrRUfWt5TX6E/vr9bSLXs1sU+KHj9noMKC+TwCwG9paRTz8A4gQGQlR+vRswZo5fb9+uvH60zP8Srb99bq8le/0/kvLFVVXZOmnjtI068d4VVBLEndOkbqzSuG664Temnm+l363YvfqrKWh4EAgCdwoR0QQCb27aQ/HNNDT80rVN8usTp/WGBfeNfkcuulr7fosdkFcliW/nR8ji4ekeHVR18dDkuXj+yu5Jgw3fTOSp313BK9culQdYoNNz0NAHwaUQwEmBvH99TanZX668frlJMSrSFdvetoaHtZvm2f/vzBGuXvOqDxvZP1t5P6qHOc74Tl5AGd1SEyRFe+vkynT1usVy89QlnJ0aZnAYDP4vQJIMA4HZaeOGeQusSF68rXlmnT7gOmJ7Wrqvom3fXhGp3+zGJV1jXpuQuH6F+/y/WpIP6vET066u2rhqvJbeuMZ5cor3iv6UkA4LOIYiAAxYYH66WLh8rhsHTuv5aqcE+16UltzrZtzVhTqnGPLtB/lm7TJSO6adZNozWhj2/f+7dP51h9cM0IJUSG6PwXlmrOht2mJwGATyKKgQDVPTFKb14xXJJ03r++0eYy/w3jXZX1uur1ZbrmjeVKjA7VR78/WndP7q2oUP84gywtIULvXX2kslOide0by7VsK0eMAeBQEcVAAOuRFKU3rxgml9vWuf/6RlvKa0xP8ii329a/v9mq8VMW6KtNZbrj+Bx99Puj1C811vQ0j+sQFapXLjlCnWLDdPmreX73uQSAtkYUAwEuKzla/7liuJpcts59/httrfCPmCrcc0BnP79Ed324VgPS4jTzj6N01ehMBTn99z97CZEheuWSI2RZli55+VvtrWk0PQkAfIb//t8BQItlp0TrjcuHqb7ZpXOf/0bb99aannTY6hpdevTLjZr0xNcq2F2tR84coNcvO0JdO0SantYuMjpG6l+/y1VpZb0uf/U71Te5TE8CAJ9AFAOQJPXqFKM3Lh+mmkaXzn5uidbsqDQ96ZDYtq0v1+3SuCkL9OTcQp3Qv5Nm3zRaZwxJDbjHIQ/pGq/Hzx6oFdv368a3V8rtbv8nlwKAryGKAfyoT+dYvXH5MEnS6c8s1utLimXiUfCHamtFjS595Ttd+foyRYUG6e0rh+uxswcqMTrU9DRjju/XSXdO6qUZa3fp/hkbTM8BAK/nH5deA/CYvl1i9dn1I3XjOyv1l4/W6dvifbr/tH5eeaeGukaXnllQpGcXFCnE6dBdJ/TSRSMyFOzH5w0fisuO7qYd++r0r4VblBofoYtGZJieBABey/v+LwfAuPjIEL100VA9s6BIj365UetKKjXtgsHKSYkxPU2SVN/k0n+WbtO0+UUqr27QyQM7685JvZQUE2Z6mlexLEt/ObG3duyr072frle/1FgNTo83PQsAvJJl4kujubm5dl5eXru/L4BD983mCl335gpV1TXp7yf31Zm55s7RbWx265287XpqbqF2VdXryO4ddPNxPZWbEZiPqm6pyromTXpioRwO6bPrRyomLNj0JABoN5ZlLbNtO/egH0cUAziYsgMNuuGtFVpcVKEBqbH647ieGpOd2G5x3ORya/ryEk2du0k79tUpt2u8bjqup0ZkdmyX9/cHy7bu01nPLdGkfp009ZyBAXfxIYDARRQD8CiX29a7edv11LxC7dhX1y5xXLD7gN5btkMfLC9ReXWD+qfG6qbxPTW6Z/sFuT95el6hHp65UQ+d0V9n5aaZngMA7YIoBtAmGpvd+mD5jv+J4xvGZWl0zyQ5Ha0P1craJn28eqfey9uuVTsqFeSwdGxOks45Ik3HZCcRw63gctu68MWlWrFtvz657mj1SIoyPQkA2hxRDKBN/TyOo0ODlJsRr2HdO2h49w7q2znmoE+Ps21bpZX1WlNSqbUllVq9o1JLNleosdmtnJRonTEkVacM6qKOUYF7azVP211Vr+OfWKjkmDBNv3aEwoKdpicBQJsiigG0i8Zmt75cv0uLCiu0dEuFNpd9/5joyBCnBneNV2J0qJyWJYdlyeGw5HRIlixt31ertSWVKq/+/lHEToelrKQoDe/eQWcMSVWfzjEcFW4jc/N369JX8nTJURm6Z3If03MAoE21NIq5JRuAVgkJcujE/p11Yv/OkqQ9B+q1dPNeLd1SobzifdpcViPbtuWybbnc+vGfk6PDNLpnkvqnxqpvl1j17hSj8BCOWraHY3OSdclRGXp5UbGO7tFRY3slm54EAMZxpBgAAlBDs0unPr1YpZV1mn3TaHXgFBUAfqqlR4p57BMABKDQIKeeOGegqhua9fdP15ueAwDGEcUAEKCykqN17Zge+nDlTs3fuMf0HAAwiigGgAB27TGZ6pEUpTunr1VNQ7PpOQBgDFEMAAEsNMip+0/rp5L9dZoyq8D0HAAwhigGgAA3NCNB5w9L18uLtmjV9v2m5wCAEUQxAEC3H5+jxOhQ/emDNWpyuU3PAYB2RxQDABQTFqy/ndRXG0qr9MLCLabnAEC7I4oBAJKkiX1TNKFPsh6fXaDi8hrTcwCgXRHFAIAf3XtyX4U4HbrzwzUy8XAnADCFKAYA/Cg5Jky3H5+jRYUV+njVTtNzAKDdEMUAgP9x7hHp6tM5Rg/MyFddo8v0HABoF0QxAOB/OB2W7pncR6WV9Xp2QZHpOQDQLjwSxZZlTbQsa6NlWYWWZf3JE68JADDniG4JOrF/Jz27oEg79tWangMAba7VUWxZllPS05KOl9Rb0rmWZfVu7esCAMy6Y1IvWZZ0/4x801MAoM154kjxEZIKbdvebNt2o6S3JJ3sgdcFABjUJS5cV4/O1GerS7V0c4XpOQDQpjwRxV0kbf/Jj3f88HMAAB931ahMdY4N098+WS+Xm1u0AfBfnohi6xd+7v/8l9OyrCsty8qzLCuvrKzMA28LAGhr4SFO3TGpl9aXVumdvO0H/w0A4KM8EcU7JKX95Mepkv7CMeUsAAAgAElEQVTPzS1t237etu1c27ZzExMTPfC2AID2cGL/ThqaEa9HZm5UZV2T6TkA0CY8EcXfScqyLKubZVkhks6R9LEHXhcA4AUs6/tbtO2tbdSTczaZngMAbaLVUWzbdrOkP0iaKWmDpHds217X2tcFAHiPvl1idXZuml5ZXKyismrTcwDA4zxyn2Lbtj+3bbunbduZtm3/0xOvCQDwLjcfl63QIIce+oJbtAHwPzzRDgDQIonRobpqdKZmrtutZVv3mZ4DAB5FFAMAWuyyo7upY1SoHpixQbbNLdoA+A+iGADQYpGhQfrjuCx9V7xPszfsMT0HADyGKAYAHJKzh6ape8dIPfhFvppdbtNzAMAjiGIAwCEJdjp028RsFe6p1nvLdpieAwAeQRQDAA7ZhD4pGpwep8dmF6iu0WV6DgC0GlEMADhklmXpjkm9tLuqQS8t2mJ6DgC0GlEMADgsQzMSNK5Xsp6dX6S9NY2m5wBAqxDFAIDDdvvEbNU0NuupuYWmpwBAqxDFAIDDlpUcrTOHpOn1b4q1fW+t6TkAcNiIYgBAq9w4vqecDktTZhWYngIAh40oBgC0SkpsmC4akaEPV5aoYPcB03MA4LAQxQCAVrt6VKaiQoL06JcbTU8BgMNCFAMAWi0+MkRXjOqumet2a9X2/abnAMAhI4oBAB5x6dHdlBAZokc4WgzABxHFAACPiAoN0rVjMrVwU7kWF5WbngMAh4QoBgB4zAXDu6pTbJgenrlRtm2bngMALUYUAwA8JizYqevHZmnFtv2as2GP6TkA0GJEMQDAo84YkqqMDhF65MuNcrs5WgzANxDFAACPCnY6dOP4nsrfdUCfrN5peg4AtAhRDADwuMn9OysnJVqPzSpQk8tteg4AHBRRDADwOIfD0i3HZau4olbv5u0wPQcADoooBgC0ibG9kjQoPU5Pzt2k+iaX6TkA8JuIYgBAm7Cs748Wl1bW642l20zPAYDfRBQDANrMUT06akRmB02bV6iahmbTcwDgVxHFAIA2dcuEbFXUNOqVxcWmpwDAryKKAQBtanB6vMbmJOm5BUWqrGsyPQcAfhFRDABoczcd11NV9c3611ebTU8BgF9EFAMA2lyfzrE6oX8nvbRoi8qrG0zPAYD/gygGALSLG8f1VH2TS8/MLzI9BQD+D6IYANAueiRF6bTBqXr9m60qrawzPQcA/gdRDABoNzeMzZJt23pybqHpKQDwP4hiAEC7SUuI0DlD0/XOd9u1raLW9BwA+BFRDABoV9cd20NBTkuPzy4wPQUAfkQUAwDaVVJMmC46MkPTV5Zo464DpucAgCSiGABgwDVjMhUVGqSHZ240PQUAJBHFAAAD4iJCdPXoTM3esFvLtu41PQcAiGIAgBmXHJWhjlGhevCLjbJt2/QcAAGOKAYAGBEREqTrx/bQt1v2akFBmek5AAIcUQwAMOacoelKSwjXQ19slNvN0WIA5hDFAABjQoIcuml8T60vrdKna0pNzwEQwIhiAIBRJw3oopyUaD365UY1udym5wAIUEQxAMAop8PSrROytbWiVu/kbTc9B0CAIooBAMYdm5OkIV3j9cTsTaprdJmeAyAAEcUAAOMsy9LtE3O050CDXllcbHoOgABEFAMAvMIR3RJ0THainplfqH01jabnAAgwRDEAwGv86fheqm5o1pNzC01PARBgiGIAgNfITonWWblpev2bYhWX15ieAyCAEMUAAK9y0/ieCnY69OAX+aanAAggRDEAwKskxYTpqlGZmrF2l/KK95qeAyBAEMUAAK9zxahuSo4J1T8+2yDb5vHPANoeUQwA8DoRIUG6+bhsrdy+X5+u5vHPANoeUQwA8EqnD05VTkq0HvwiXw3NPNADQNsiigEAXsnpsHTXCb21Y1+dXlu81fQcAH6OKAYAeK2jszpqTHainpy7iQd6AGhTRDEAwKv9edL3D/SYOneT6SkA/BhRDADwaj2To3X20HS9vmSrCvdUm54DwE8RxQAAr3fzcT0VEeLUXz9exy3aALQJohgA4PU6RoXq1gnZ+rqwXJ+v2WV6DgA/RBQDAHzCecO6qk/nGP390/WqaWg2PQeAnyGKAQA+wemwdO/JfbWrqp6L7gB4HFEMAPAZQ7rG66zcVL24cIsK9xwwPQeAHyGKAQA+5faJOYoIceruj7joDoDnEMUAAJ/SISpUt07M0eKiCn26utT0HAB+gigGAPic845IV98uMfrHZ+tVzUV3ADyAKAYA+Bynw9LfT+6r3VUNmjqHi+4AtB5RDADwSYPS43XO0DS99PUWbSitMj0HgI8jigEAPuv2iTmKiwjWze+sUpPLbXoOAB9GFAMAfFZ8ZIj+eWo/rS+t0lNzC03PAeDDiGIAgE+b0CdFpw7qoqfnFWptSaXpOQB8VKui2LKshy3Lyrcsa7VlWdMty4rz1DAAAFrqnsm9lRAZolveXaWGZpfpOQB8UGuPFM+S1Ne27f6SCiTd0fpJAAAcmriIED1wej/l7zqgJ+dwGgWAQxfUmt9s2/aXP/nhN5LOaN0cAAAOz7E5yTpzSKqeWVCk8b2TNSCNL16aYtu2KuuaVFpZr12V9ZIldY4NV6e4MMWEBZueB/yiVkXxz1wq6e1f+0XLsq6UdKUkpaene/BtAQD43l0n9tbXheW6+d1V+vS6oxUW7DQ9ye/VN7m0uKhcc/P3aHNZjUor61VaWaf6pl++G0hUaJA6x4WpU2y4clKiNXlAZ/XpHCPLstp5OfC/rIM9N96yrNmSUn7hl+60bfujHz7mTkm5kk6zW/Ag+tzcXDsvL+8w5gIA8NsWFJTpope+1VWju+uO43uZnuOXyg40aF7+Hs3asFtfbypXXZNLESFO9eoUo5TYMHWKCfv++9hwpcSGSbK1c//3sbxzf7127q9TaWW98ndVqcllKyspSqcM6qKTB3ZWanyE6T8e/IxlWcts28496Me1oGEP9kYXSbpa0ljbtmtb8nuIYgBAW7rjgzV6+7ttevuqIzU0I8H0HL/gctv6Yu0uvbJ4i/K27pNtS51jwzS2V7LG9U7W8O4JCg06tCPz+2sb9dmaUn24okTfFe+TJB3RLUFnDknVaYNT5XRw9Bit1y5RbFnWRElTJI22bbuspb+PKAYAtKXqhmadOHWhahtd+vT6o5UUHWZ6ks+qa3Tp3WXb9cLCLdq2t1YZHSJ06qBUjeudpN6dPHfaw/a9tfp41U59sHyHispq1K9LrP55al/1T+XccLROe0VxoaRQSRU//NQ3tm1ffbDfRxQDANpa/q4qnfr0YvXrEqs3rhimYCe35j8UFdUNem3JVr22pFj7aps0MC1OV4/urvG9U9r0CK5t2/psTanu/WS9yqob9LvhXXXzhGwu0MNha7fTJw4HUQwAaA8frSzRDW+t1CVHZeieyX1Mz/EJNQ3Nem5BkZ5fuFn1TW6N65WkK0dlamhGfLteDFdV36QpXxbotSXF6hAVqr+c2FuT+3figjwcspZGsSfvPgEAgFc5eWAXrdy+Xy8vKtbAtDidPLCL6Uley+229f7yHXp45kbtOdCgyQM664axPdQjKdrInpiwYP31pD46fXCq7vxwja5/c4XezduuKWcNVGJ0qJFN8G8cKQYA+LUml1vn/2up1pRUavrvRygnJcb0JK/z7Za9+vun67WmpFID0uJ094m9NKSr91yg6HLbemPpVt33+QYlRofqlUuOUGZilOlZ8BEtPVLMCVYAAL8W7HToqfMHKTosSFe/vkyVdU2mJ3mNkv11uvaNZTrruSUqr27Q42cP1PRrRnhVEEuS02Hpd0dm6K0rj1Rdo0unP7NY3xXvNT0LfoYoBgD4vaToME07f7B27KvTze+slNvd/l8l9SbNLrdeWLhZ46cs0Nz8PbpxXE/NvXmMThnURQ4vvg3awLQ4fXDNUUqICNH5LyzVp6t3mp4EP0IUAwACQm5Ggu46oZdmb9ijh7/caHqOMWt2VOqUaYv0j882aFi3BM26cbRuGJel8BDfePpfeocIvX/NCA1IjdUf/rNCzy0okolTQeF/uNAOABAwLhqRoYI91XpmfpEigp26bmyW6UntpqahWY9+WaBXFm9Rh6hQPX3eYE3ql+KTd3OIjwzR65cN083vrtL9M/K1fV+t/jq5j4K47R5agSgGAAQMy7L0j5P7qr7JpUdnFSgs2KkrRnU3PavNzdmwW3/5cK12Vtbr/GHpum1ijmLDffu+v2HBTj15ziClxoXrua82y+W2dd+p/Xwy8uEdiGIAQEBxOCw9dHp/NTS79c/PNygs2KELj8wwPatN7DlQr799sl6frS5Vz+QovX/ekV53EV1rOByW7pjUS06HpWnzi9QpNlzXB9DRf3gWUQwACDhBToceP3ugGprc+stH6xQa5NRZQ9NMz/IY27b1bt4O/eOz9apvcuvm8T111ehMhQT55+kFt07I1q6qek2ZVaCUmDC/+lyi/RDFAICAFOx06OnzB+nyV/N0+werFRrs8IuHe2wpr9GfP1ijJZsrdES3BN1/Wj+/v6evZVl68PT+Kq9u1B3T1ygxOlTH5CSZngUf459/ZQQAoAVCg5x6/sJcDeuWoJveWaWPV/nuLb7qm1x6cs4mTXj8K63dWan7Tu2nt64Y7vdB/F/BToemnT9YvTpF69o3lmvV9v2mJ8HHEMUAgIAWHuLUixcN1eD0OF3/5grdP2ODml1u07MOydz83Zrw+Fd6dFaBxvVK0pybRuu8Yelefc/hthAVGqSXLh6qjtEhuvSV71RcXmN6EnwIUQwACHiRoUH69+XDdN6wdD23YLN+99K3qqhuMD3roLZV1OryV7/Tpa/kyemw9PplR2ja+UOUFBNmepoxSdFhevWSI+S2bV308rcq94HPI7wDUQwAgL4/leK+U/vpoTP6K2/rPk1+8muv/RJ8XaNLU2YVaNxjC7SkqEJ3HJ+jL24YpZFZiaaneYXuiVF66eKh2l1Vr+vfXCFXgD/BEC1DFAMA8BNn5abp/atHyLIsnfnsEr357TbTk35U3+TSa0uKdcwj8zV1ziYd3zdFc24e49d3ljhcg9Ljde/JfbW4qELT5hWangMfwN0nAAD4mX6psfr0uqN1/VsrdMcHa/TN5grdOiFbqfERRvY0NLv0znfb9fS8Iu2qqtfQjHg9cc5ADevewcgeX3HmkFQtLizXY7MLdES3BP594TdZJp4Xnpuba+fl5bX7+wIAcChcbltT52zSMwuKJFs6b1i6/nBsD3WMCm2X9/9vDE+bX6TSynrldo3XjeN7akRmB57c1kLVDc2a/OTXqmt06fMbRiohMsT0JLQzy7KW2bade9CPI4oBAPhtO/fXaeqcTXp32Q6FBjl02dHddMWo7ooJ8/yjkm3b1vrSKk1fXqKPVu1U2YEGYriV1u2s1KlPL9bRWR314kW5/DsMMEQxAAAetrmsWlNmFejT1aWKDQ/WJUdlaFyvZPXuFNPq25/t3F+nj1bu1PQVO1Swu1rBTkvHZCfpwiO76ugeHQm5VnptSbHu/mid7jqhly4f2d30HLQjohgAgDaytqRSj3y5UfM3lkmSOkSG6KgeHTUyq6NGZiUqJfa3b4lW3+RS4Z5q5e86oPzSKq3asV95W/fJtqXcrvE6ZVAXndi/k+Ii+FK/p9i2rav/vUxzNuzRe9eM0MC0ONOT0E6IYgAA2tieqnp9XViuhZu+//bfe+KmJYQrOjRYYcEOhQU7FRr0/feSVLinWpvLa368TVhokEPZKdEam5OsUwd1UXoHMxfzBYLK2iZNmrpQDof06XUjFRvu+dNf4H2IYgAA2pFt28rfdUALN5Vp9Y5K1Te5VN/kVkPz//++2W2re8dI5aTEqFenGOV0ilZGh0g5A+zJcyYt37ZPZz27RCf076Qnzhlkeg7aQUujmFuyAQDgAZZlqVen72MX3mtwerz+cGwPPT57k04e2FnH5iSbngQvwZ2+AQBAQLl2TA/1TI7SXdPXqrqh2fQceAmiGAAABJSQIIceOL2/Sqvq9dAX+abnwEsQxQAAIOAMTo/XxSMy9Po3W/Vd8V7Tc+AFiGIAABCQbjkuW51jw3X7+6tV3+QyPQeGEcUAACAgRYYG6b7T+mlzWY2emltoeg4MI4oBAEDAGt0zUacN6qJnFxRpQ2mV6TkwiCgGAAAB7S8n9lZseLBuf3+1ml1u03NgCFEMAAACWnxkiO45qY9W76jUy4uKTc+BIUQxAAAIeJP7d9KxOUl6bHaBdlfVm54DA4hiAAAQ8CzL0j2Te6vZZevBGdy7OBARxQAAAJK6dojU5SO76YMVJVq2dZ/pOWhnRDEAAMAPfn9MDyXHhOpvn6yT222bnoN2RBQDAAD8IDI0SHcc30urd1TqvWU7TM9BOyKKAQAAfuLkgZ01pGu8HpqZr6r6JtNz0E6IYgAAgJ+wLEt/ndxHFTWNmjp7k+k5aCdEMQAAwM/0S43V2blpemVxsQr3VJueg3ZAFAMAAPyCWyZkKzzEqb9/ul62zUV3/o4oBgAA+AUdo0J1w9gsLSgo09z8PabnoI0RxQAAAL/iohEZ6pEUpb9/ul6NzW7Tc9CGiGIAAIBfEex06M4Teqm4olZvfrvN9By0IaIYAADgN4zpmagju3fQE3M26QC3aPNbRDEAAMBvsCxLd0zK0d6aRj3/1WbTc9BGiGIAAICD6J8ap8kDOutfCzdrd1W96TloA0QxAABAC9x6XLZcbluPzy4wPQVtgCgGAABogfQOEbpgeFe9/d12bdp9wPQceBhRDAAA0ELXHZulyJAgPfhFvukp8DCiGAAAoIUSIkN0zTGZmr1hj5ZurjA9Bx5EFAMAAByCS4/qppSYMN0/I5/HP/sRohgAAOAQhAU7ddNxPbVy+37NWLvL9Bx4CFEMAABwiE4fnKrs5Gg99EW+mlw8/tkfEMUAAACHyOmwdPvx2SquqNU7edtNz4EHEMUAAACH4ZjsJOV2jdfUOZtU3+QyPQetRBQDAAAcBsuydOuEbO2uatBrS4pNz0ErEcUAAACHaVj3DhrVM1HT5hfpQH2T6TloBaIYAACgFW49Llv7a5v0r4VbTE9BKxDFAAAArdAvNVaT+qXoxYWbVVHdYHoODhNRDAAA0Eo3jc9WXZNLz8wvMj0Fh4koBgAAaKUeSVE6fXCqXvtmq3burzM9B4eBKAYAAPCAG8ZlSbb05NxNpqfgMBDFAAAAHpAaH6HzhqXrnbwd2lJeY3oODhFRDAAA4CG/P6aHQpwOTZlVYHoKDhFRDAAA4CGJ0aG69OgMfbJqp9bvrDI9B4eAKAYAAPCgK0dlKiYsSFNmbTQ9BYeAKAYAAPCg2PBgXTU6U7M37NHybftMz0ELEcUAAAAedvGIDHWMCtGjX3K02FcQxQAAAB4WGRqka8b00KLCCi0uKjc9By3gkSi2LOsWy7Jsy7I6euL1AAAAfN35w9LVKTZMj8zcKNu2Tc/BQbQ6ii3LSpM0XtK21s8BAADwD2HBTl13bJaWb9uveRv3mJ6Dg/DEkeLHJN0mib8CAQAA/MSZuanq2iFCj8wskNtNKnmzVkWxZVknSSqxbXtVCz72Ssuy8izLyisrK2vN2wIAAPiEYKdDfxyXpfWlVZqxdpfpOfgNB41iy7JmW5a19he+nSzpTkl3t+SNbNt+3rbtXNu2cxMTE1u7GwAAwCecNKCLspKiNGXWRrk4Wuy1DhrFtm2Ps22778+/SdosqZukVZZlFUtKlbTcsqyUtp0MAADgO5wOSzcf11NFZTWavqLE9Bz8isM+fcK27TW2bSfZtp1h23aGpB2SBtu2zdcGAAAAfmJCnxT16xKrx2cXqLHZbXoOfgH3KQYAAGhjlvX90eId++r09nfcsMsbeSyKfzhizN2pAQAAfsHonok6IiNBU+cWqrax2fQc/AxHigEAANqBZVm6bWK2yg406JXFxabn4GeIYgAAgHaSm5GgY3OS9Oz8IlXWNpmeg58gigEAANrRrROydaChWc9+VWR6Cn6CKAYAAGhHvTrF6KQBnfXyoi3aU1Vveg5+QBQDAAC0s5vG91Szy9aTcwtNT8EPiGIAAIB21rVDpM4emqY3v92mbRW1pudARDEAAIAR14/NUpDT0pRZG01PgYhiAAAAI5JjwnTxiG76aNVO5e+qMj0n4BHFAAAAhlwzOlNRoUF6ZCZHi00jigEAAAyJjQjW1aMzNXvDHuUV7zU9J6ARxQAAAAZdclSGEqND9cCMfNm2bXpOwCKKAQAADIoICdKN43oqb+s+zVy3y/ScgEUUAwAAGHZWbqqykqL0wIx8NTa7Tc8JSEQxAACAYUFOh+6YlKPiilr9Z+lW03MCElEMAADgBY7JTtKIzA56Ys4mVdU3mZ4TcIhiAAAAL2BZlv48qZf21zVp2rwi03MCDlEMAADgJfp2idWpA7vopUVbVLK/zvScgEIUAwAAeJGbJ2TLknigRzsjigEAALxIl7hwXXp0N01fUaK1JZWm5wQMohgAAMDLXDMmUwmRIfrnZxt4oEc7IYoBAAC8TExYsG4Ym6Ulmys0b+Me03MCAlEMAADghc4blq7uHSP1j0838ECPdkAUAwAAeKFgp0N3T+6tzeU1evHrLabn+D2iGAAAwEuNyU7S+N7JenLuJpVWcou2tkQUAwAAeLG7T+wtl9vWPz/bYHqKXyOKAQAAvFhaQoSuGZOpT1eXanFRuek5fosoBgAA8HJXj85UWkK47vlonZpcXHTXFohiAAAALxcW7NTdJ/bRpj3VenVxsek5fokoBgAA8AHjeiXpmOxEPT57k/ZU1Zue43eIYgAAAB9gWZbumdxHjc1u3T8j3/Qcv0MUAwAA+IiMjpG6clR3TV9Rom+37DU9x68QxQAAAD7k2mMy1Tk2TH/5cC1PuvMgohgAAMCHRIQE6e+n9NXG3Qf01LxC03P8BlEMAADgY8b2StZpg7po2rxCrS2pND3HLxDFAAAAPujuyb0VHxmiW99bzWkUHkAUAwAA+KC4iBDdd2o/bSit0rT5nEbRWkQxAACAjxrfO1mnDOysp+YWat1OTqNoDaIYAADAh90zuY/iIkJ0y7ureQR0KxDFAAAAPiw+MkT/PLXv96dRzCsyPcdnEcUAAAA+bkKfFJ00oLOenLtJ63dWmZ7jk4hiAAAAP/DXk/ooLiJYt7y7Sg3NLtNzfA5RDAAA4AcSIkN0/2n9tb60Svd+st70HJ9DFAMAAPiJ8b2TddWo7npj6TZ9sHyH6Tk+hSgGAADwI7dOyNYR3RL05+lrlL+L84tbiigGAADwI0FOh546b5Ciw4J1zb+X60B9k+lJPoEoBgAA8DNJ0WF66txB2ra3Vre9t1q2bZue5PWIYgAAAD80rHsH3T4xWzPW7tKLX28xPcfrEcUAAAB+6oqR3TWhT7IemJGv74r3mp7j1YhiAAAAP2VZlh4+c4BS48P1+zeWq2R/nelJXosoBgAA8GMxYcF69sIhqmty6cIXlqrsQIPpSV6JKAYAAPBzOSkxevnioSqtrNfvXvpWlbXckeLniGIAAIAAkJuRoOcuHKKiPdW65JVvVdPQbHqSVyGKAQAAAsSonomaeu5Ardy+X1e9vkz1TS7Tk7wGUQwAABBAJvbtpIfOGKCvC8t1/Zsr1Oxym57kFYhiAACAAHPGkFT9dXJvfbn+/7V37zFSlWccx78/XS5FiC5SFMULeKta24qrVdRGxaiQRmpTW9qaUjExUm30j14wJMbWP1qttqmp1lolttWI1is1GMWC2ktE0XAzgCyIFaEsaIt4ARZ5+sd5V6fLDDvDwpzZPb9PcrKz57yz++yT57z7zDlnzqzjRw8tZPt2f7hHU94BmJmZmVn9fff0EWzavI1bZr3Gu5vb+dU3vsCg/n3yDis3PlJsZmZmVlBXnXMkP7nweOYsW89Ft/+Tlevfyzuk3LgpNjMzMysoSUwcfTj3XvZF3nl/K+Nv+wfPLmvLO6xcuCk2MzMzK7jTjtifx688neHNA5h0z0v87rkVRBTrOmM3xWZmZmbGIYMH8PDk0xj72WH87MmlXD19fqHuZeym2MzMzMwAGNC3id9860R+eP4x/GXhGs6++Vmmv/gvPirA3SncFJuZmZnZxyRx5dlH8tAVoxne/CmmPLKIsb9+njlL23r1JRVuis3MzMxsBycd1szDk0dz+7dHsWXbdi695yUuuXsui9/amHdoe4Ty6PhbWlpi3rx5df+9ZmZmZla7rdu2c+8Lb3Dr7OVs/LCd048YwphjhzLmMwdw6P4D8g5vpyS9HBEtXY5zU2xmZmZm1dj4YTt3/W0lMxetZcX69wE4cujAjxvkEw7el/599kJSzpF+wk2xmZmZme0xqza8z+ylbcxe2sbc19+m/aOsp2zaSwzs38Sg/k0M7NeHQf2bOHTwAG6++PO5xFltU+yPeTYzMzOzmh0+ZB8mnTGCSWeMYNPmdv6+fAOr3v6ATZvbeW/LNjZt7lja2bS5Pe9wu+Sm2MzMzMy6ZVD/Pow9YVjeYXSL7z5hZmZmZoXX7aZY0vclLZP0qqSbdkdQZmZmZmb11K3LJySdDYwHPhcRWyQN3T1hmZmZmZnVT3ePFE8Gfh4RWwAioq37IZmZmZmZ1Vd3m+KjgTMlzZX0nKSTKw2UdLmkeZLmrV+/vpu/1szMzMxs9+ny8glJzwAHltk0NT2/GTgVOBl4UNLIKHPz44i4E7gTsvsUdydoMzMzM7PdqcumOCLOrbRN0mTgkdQEvyhpOzAE8KFgMzMzM+sxunv5xGPAOQCSjgb6Ahu6G5SZmZmZWT1198M7pgHTJC0GtgITy106YWZmZmbWyLrVFEfEVuCS3RSLmZmZmVku/Il2ZmZmZlZ4borNzMzMrPDcFJuZmZlZ4bkpNjMzM7PCc1NsZmZmZoXnptjMzMzMCs9NsZmZmZkVnptiMzMzMys85fEBdJLWA2/U/Vs84OMAAAZOSURBVBfDEPwx1LVwvmrnnNXG+aqN81Ub56s2zldtnK/a5JmvwyLi010NyqUpzoukeRHRknccPYXzVTvnrDbOV22cr9o4X7VxvmrjfNWmJ+TLl0+YmZmZWeG5KTYzMzOzwitaU3xn3gH0MM5X7Zyz2jhftXG+auN81cb5qo3zVZuGz1ehrik2MzMzMyunaEeKzczMzMx24KbYzMzMzAqv1zXFki6W9Kqk7ZJaOm27VlKrpGWSzq/w/BGS5kpaLukBSX3rE3n+0t87Py2rJM2vMG6VpEVp3Lx6x9lIJF0v6a2SvI2rMO6CVHetkqbUO85GIekXkpZKWijpUUn7VRhX2BrrqlYk9Uv7amuaqw6vf5SNQ9IhkuZIWpLm/qvLjDlL0saS/fS6PGJtFF3tX8rcmmpsoaRRecTZCCQdU1I38yW9K+maTmMKXV+Spklqk7S4ZN1gSbNSLzVLUnOF505MY5ZLmli/qCuIiF61AMcCxwDPAi0l648DFgD9gBHACmDvMs9/EJiQHt8BTM77b8opj7cA11XYtgoYkneMjbAA1wM/6GLM3qneRgJ9Ux0el3fsOeXrPKApPb4RuLHCuELWWDW1AnwPuCM9ngA8kHfcOedsGDAqPR4EvFYmZ2cBT+Qda6MsXe1fwDjgSUDAqcDcvGNuhCXtn/8m+yCI0vWFri/gS8AoYHHJupuAKenxlHJzPTAYWJm+NqfHzXn+Lb3uSHFELImIZWU2jQemR8SWiHgdaAVOKR0gScA5wENp1R+Ar+zJeBtRysPXgfvzjqWXOAVojYiVEbEVmE5Wj4UTEU9HxLb07QvA8DzjaUDV1Mp4srkJsrlqTNpnCyki1kbEK+nxJmAJcHC+UfV444E/RuYFYD9Jw/IOqgGMAVZERB6fyNuwIuJ54J1Oq0vnqUq91PnArIh4JyL+A8wCLthjgVah1zXFO3Ew8GbJ96vZceLcH/hvyT/tcmOK4ExgXUQsr7A9gKclvSzp8jrG1aiuSqcYp1U4RVRN7RXRJLKjUeUUtcaqqZWPx6S5aiPZ3FV46VKSE4G5ZTafJmmBpCclHV/XwBpPV/uX56zyJlD5YJHr6/8dEBFrIXvhCgwtM6bh6qwpz1++qyQ9AxxYZtPUiHi80tPKrOt8P7pqxvRoVebum+z8KPHpEbFG0lBglqSl6ZVir7SznAG/BW4gq5MbyC47mdT5R5R5bq+qq1LV1JikqcA24L4KP6ZQNVbC89QukjQQeBi4JiLe7bT5FbJT3u+l6/4fA46qd4wNpKv9yzXWSXp/0YXAtWU2u752TcPVWY9siiPi3F142mrgkJLvhwNrOo3ZQHaaqCkdgSk3pkfrKneSmoCvAift5GesSV/bJD1Kdsq31zYs1dabpN8DT5TZVE3t9RpV1NhE4MvAmEgXlpX5GYWqsRLV1ErHmNVpf92XHU9dFoqkPmQN8X0R8Ujn7aVNckTMlHS7pCERsaGecTaKKvavQs1ZVRoLvBIR6zpvcH2VtU7SsIhYmy69aSszZjXZ9dgdhpO9Hyw3Rbp8YgYwIb1zewTZq7gXSwekf9BzgK+lVROBSkeee6tzgaURsbrcRkn7SBrU8ZjsjVOLy40tgk7X2V1E+Vy8BByl7M4mfclOwc2oR3yNRtIFwI+BCyPigwpjilxj1dTKDLK5CbK5analFxdFkK6nvhtYEhG/rDDmwI7rriWdQva/7+36Rdk4qty/ZgDfSXehOBXY2HEqvMAqnkF1fZVVOk9V6qWeAs6T1JwuPTwvrctPnu/y2xMLWWOyGtgCrAOeKtk2leyd3cuAsSXrZwIHpccjyZrlVuDPQL+8/6Y65+8e4IpO6w4CZpbkZ0FaXiU7JZ573Dnm60/AImAh2SQwrHPO0vfjyN4Vv6LIOUv71ZvA/LR03EXBNbaTWgF+SvZCAqB/mpta01w1Mu+Yc87XGWSnXBeW1NU44IqOuQy4KtXSArI3eI7OO+4c81V2/+qULwG3pRpcRMmdnIq4AAPImtx9S9a5vj7Jxf3AWqA99V+Xkb3P4a/A8vR1cBrbAtxV8txJaS5rBS7N+2/xxzybmZmZWeEV6fIJMzMzM7Oy3BSbmZmZWeG5KTYzMzOzwnNTbGZmZmaF56bYzMzMzArPTbGZmZmZFZ6bYjMzMzMrvP8BYeXHY7zoXiUAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "def f( x ):\n", " \n", " return np.cos( x ) - 0.5 * x\n", "\n", "x = np.linspace( -10, 10, 100 )\n", "\n", "pp.figure( figsize = ( 12, 8 ) )\n", "pp.plot( x, f( x ) )\n", "pp.show()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "n = 100\n", "s = np.random.normal( size = 100 )\n", "\n", "def f( mu, sigma ):\n", " \n", " retval = np.empty( mu.shape )\n", " for i in range( mu.shape[ 0 ] ):\n", " retval[ i ] = -np.prod( np.exp( -( s - mu[ i ] ) ** 2 / sigma ) )\n", " \n", " \n", " return retval" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.2230254436752806\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3YAAAHdCAYAAABR6euJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzs3Xd0nNd95//PnRn0MgCIQlSCVSQBdojqxZYUy3GRJduJyzpybK8SJ042m1+yP+9vk3NyUjZOnN3kOE52V3Ecy6mOvbIty1LUY0mWKIliA8BeUUkQbdDbzP39AYCWJZAEOIO5zzPzfp3Dg/Zono8oicIH9z7fa6y1AgAAAAD4V8B1AAAAAABAfCh2AAAAAOBzFDsAAAAA8DmKHQAAAAD4HMUOAAAAAHyOYgcAAAAAPuf5YmeM+boxpscY05Kg1/s3Y8ygMebxy3z9L40xI4m4FwAAAAAkg+eLnaRvSLo3ga/3ZUmfWugLxpgmSUUJvBcAAAAALDvPFztr7YuS+t/6OWPM2rmVtzeNMS8ZYzYu4fWekzT89s8bY4KaLX3/Jd7MAAAAAJBMIdcBrtHDkn7ZWnvCGHODpL+W9O44X/MLkh6z1nYbY+IOCAAAAADJ4rtiZ4zJl3SzpG+/pYBlzX3tAUm/v8Bf1mmtfc8VXrNK0kcl3ZnQsAAAAACQBL4rdprdPjpord3+9i9Yax+V9Og1vOYOSesknZwri7nGmJPW2nVxJQUAAACAJPD8M3ZvZ60dknTGGPNRSTKztsX5mj+01q601tZba+sljVHqAAAAAPiF54udMeafJb0q6TpjTIcx5rOSPinps8aYg5JaJd23hNd7SdK3Jd0193qX3aIJAAAAAH5grLWuMwAAAAAA4uD5FTsAAAAAwJVR7AAAAADA5zw7FbO0tNTW19e7jgEAAAAATrz55pu91tqyxVzr2WJXX1+vvXv3uo4BAAAAAE4YY84t9lq2YgIAAACAz1HsAAAAAMDnKHYAAAAA4HMUOwAAAADwOYodAAAAAPgcxQ4AAAAAfI5iBwAAAAA+R7EDAAAAAJ+j2AEAAACAz1HsAAAAAMDnKHYAAAAA4HMUOwAAAADwOYodAAAAAPgcxQ4AAAAAfC4hxc4Yc68x5pgx5qQx5osLfD3LGPOtua+/ZoypT8R9AQAAAAAJKHbGmKCkv5L0XkmbJX3cGLP5bZd9VtKAtXadpD+X9Cfx3hcAAAAAMCuUgNfYLemktfa0JBlj/kXSfZIOv+Wa+yT93tz735H0VWOMsdbaBNw/aaajMY1NRRd9vTGLf+0lXCqzhBde2usu4dolvfLSXnu5XncpmVP6n91y/cMAAADXxFqrmP3J29jct8hX+07Z6urfSl/9NRaX76rXLOp1FnGRQ/lZIQUD/v0+KRHFrlpS+1s+7pB0w+WusdbOGGMiklZI6k3A/ZPmjbP9+sTfvOY6BuBJKV2Gl/LKS7g0aIxCQaNQwCgYCCgUeOvHRqFAYMGPgwGjzGBA2ZlBZYeCyskMKCcjqOy5XzkZQeVnhxTOyVA4J0OF2RkK586+n5cZpNzDnSfnntZ475fc5kDaisashsanNTg+rYGxKUXGphUZn9bYVFRjUzOamI7OvR/V+FRU49NRTc5ENR21mo7GNB2NaSZqNR2zmp6JaSYWu/S1majVTCx2qZjFYrO1y859/Na3Vj8pcF4vO+nk2d+8Q+vK813HuGaJKHYLfYfw9n9FF3ONjDEPSXpIkurq6uJPlmD1K/L0u+9/+y7ThS3XYuRSXnYxP8W5ttddGi9kXoql/LNbrt83L/yeLem3dym/Z8vzsp74fVhq3ljMaiZmFZ17OxONvePj6Fuumf+mYnzaamomponpqCamYxqfjmpievYbkKtlCAaMinIyVF6YrcpwtlaGs1VZmK2K8OzHleFsVRflKiczuIS/c2CRzje7ToAUFhmfVtfguDoHxtU5OK6uwXF1zL3tH53S4Ni0hiamr/rnZChglJM5+0Oy3MygskJBZYRmf7iWGQwoMxRQXjCgjODs5zJCAWUEjDKCAQWDRkFjFDCzP5g0RgqY2R8PBgKzHxvNfj0w93Uzf/385+dWjOb/2qtZzI/qrvYyi/kBZqJ+JujlHy6W5me6jhCXRBS7Dkm1b/m4RlLXZa7pMMaEJIUl9b/9hay1D0t6WJKampo89/OLqqIcffbW1a5jAMCCrLWaisY0PhXV8MSMIuOz38QMjc/+RHpofPZz/WNT6hmaUHdkQgfbB9U3OvWO16ouytGasjytLcvXdSsL1FgV1oaV+coKUfgAuHVxeFItXREdPz+sYxeGdeLCiM72jmp4cuanrssMBlRVlK2qohxtqylSUW6GinIyVJSbqaLcDBXnZr5lN0NIORlB5WQGlRliaDz8KRHF7g1J640xqyV1SvqYpE+87ZrHJD0o6VVJH5H0vN+erwMArzPGKCs0+9PlotzMn/qJ25VMTEfVMzSp7si4zg9NqK1vTKcujuh076i+vbddo3PPFmcEjTZUFGjXqmLdsHqFdq8uUVlB1vL9DQFIe7GY1eHuIe053af97YM60DaozsHxS1+vKMzShooCPbCzWjXFuaoqylFVUbaqi3NUmpd1afULSAdxF7u5Z+a+IOkpSUFJX7fWthpjfl/SXmvtY5L+VtLfG2NOanal7mPx3hcAkBjZGUHVrchV3Yrcd3zNWqu2/jG1dA6ppSuiQx2D+vbeDn3z1XOSpHXl+Xr3xnLdvalCu1YV+/qhcwDeMDA6pWeOXNCLxy/qlVN96p/bVVBdlKPtdUX69M312loT1saVhQrnZjhOC3iH8erCWVNTk927d6/rGACAt5mOxtTcGdFrp/v1yqle7Tndp+moVUlepu7aWK4HdtbohtUl/KQc7/R375t9+4s/dJsDnjMwOqWnWs/rh83deuVUn6Ixq7KCLN22rlS3zP1aGc52HRNIOmPMm9bapsVcm4itmACANJIRDGhnXbF21hXr83eu1fDEtH50/KKeOXxBT7ac17ff7FBtSY4+srNWH2mqUXVRjuvIADzIWqu95wb0D3vO6cnm85qKxrRqRa4eun2N3relUg1VhZ4etAF4DcUOABCXguwMvX9rld6/tUrjU1H9W2u3vr23Q3/+7HF95fkT+sDWSv3SHWu1qbLQdVQAHjA1E9Oj+zr0jVfO6uj5YRVkhfSJG+r0kV01lDkgDhQ7AEDC5GQGdf+OGt2/o0bt/WN65JWz+qfX2/S9A12687oy/dq712vXqmLXMQE4MB2N6f++2aGvvnBSHQPjaqgq1Jce2KIPbq9SbibfkgLx4r8iAMCyqC3J1e+8f7O+8O51+vtXz+kbr5zVh//XK3r/1kp98b0bVVP8zmEtAFKPtVaPHezSnz19TO3949pWE9YffKhRd24oY3UOSCCKHQBgWRXlZurX7lqvz962Wv/7R6f18Iun9PThC/rcrav1K+9ap/ws/lcEpKqTPcP63e+16tXTfWqoKtTXP92gd11XTqEDlgH/NwUAJEVuZki/ec8Gfez6Wn35qWP6638/NftT/I9u041rVriOByCBJqaj+spzJ/Q3L51WTkZQf/ihRn18dx1HogDLKOA6AAAgvVQV5ejPf367vvPLNykYMPr43+zRHz5+WBPTUdfRACTA8QvD+uBXX9Zf//sp3be9Ws//1p36DzeuotQBy4xiBwBwoqm+RE/8+m365A11+trLZ/SBv3xZx84Pu44F4BpZa/Uvr7fpg199Wf2jU3rkM7v1Zx/dptL8LNfRgLRAsQMAOJOXFdIffmiLHvnMbg2OT+uBv/6xnm497zoWgCUanZzRr//LAX3x0WY1rSrRE//pNt2xocx1LCCtUOwAAM7dsaFMP/jCrVpbnq9f+oc39VcvnJS11nUsAIvQMzyhjz28Rz881KXffs91+uZndqu8INt1LCDtUOwAAJ6wMpytf/2lm/TBbVX68lPH9Ov/coDn7gCPO9kzogf++hWd7BnR1x5s0q++a50CPEsHOMFUTACAZ2RnBPUXP79dG1cW6k/+7agi49N6+FO7lJ0RdB0NwNu8fqZf//Gbe5URNPrWL92orTVFriMBaY0VOwCApxhj9Pk71+pPP7xVL524qM984w2NTc24jgXgLX58slf/4W9f04q8TD36+VsodYAHUOwAAJ70c9fX6n/+3DbtOd2nT//dGxqZpNwBXvDG2X597pG9Wr0iT9/5/M2qW5HrOhIAUewAAB52/44a/cXHdujNcwP69Ndf1/gUz9wBLh3qGNQv/t0bqgxn6x8+d4NK8jJdRwIwh2IHAPC0D26r0lc+tkNvtg3oN761X9EY0zIBF450D+lTf/u6ivMy9I//8QaVFXA+HeAlFDsAgOe9b2ulfvd9m/VU6wX98RNHXMcB0s75yIR+4euvKycjqH/63I2qDOe4jgTgbZiKCQDwhc/culpt/WP62stnVLciV79wU73rSEBamJiO6qG/36uxyRl991dvUW0Jz9QBXkSxAwD4xu++f7M6Bsb1e4+1qiqco7s3V7iOBKQ0a63+66PNOtQR0cOf2qUNFQWuIwG4DLZiAgB8Ixgw+srHt6uhKqzf+NYBne0ddR0JSGlfe+mMvru/U795zwb9TMNK13EAXAHFDgDgK7mZIf3vT+1SMGD0a/+8X5MzTMoElsOPjl/UHz95RD+7ZaV+7d3rXMcBcBUUOwCA71QX5ehPP7JVzZ0R/cmTx1zHAVLOxeFJ/edvHdCGigJ9+SPbZIxxHQnAVVDsAAC+9J6Glfr0zfX6+o/P6JnDF1zHAVKGtVb/7/89pJHJGX3l4zuUl8VIBsAPKHYAAN/6rz+7UQ1Vhfrt7xxU1+C46zhASvjn19v1/NEeffHejQxLAXyEYgcA8K2sUFBf/cROTc/E9P/860FZy+HlQDzO9I7qDx4/rFvXlerTN9e7jgNgCSh2AABfW12ap//vfZv06uk+fXtvh+s4gG/NRGP6z986oIyg0Zc/ulWBAM/VAX5CsQMA+N7Hr6/T7voS/dETR3RxeNJ1HMCX/s+Lp3WgfVB/dP8WVYZzXMcBsEQUOwCA7wUCRv/9gS0an4rq9x8/7DoO4Dvt/WP6y+dP6N6GlfrAtirXcQBcA4odACAlrCvP16++a51+cLBLLxztcR0H8JU/ePywjIx+9wObXUcBcI0odgCAlPHLd67RuvJ8/c73WjQ6OeM6DuALLxzt0dOHL+jX7lqn6iK2YAJ+RbEDAKSMrFBQf/LhLeqKjOsrz51wHQfwvInpqH7vB61aU5anz926xnUcAHGg2AEAUsquVSW6f0e1/u6Vs+oYGHMdB/C0h188rXN9Y/r9DzYqM8S3hYCf8V8wACDl/NbPXCcj6c+eOuY6CuBZ7f1j+qsXTup9Wyt16/pS13EAxIliBwBIOVVFOfrsrav1vQNdOtQx6DoO4En/85njMkb6nfdtch0FQAJQ7AAAKemX71yrkrxM/fcnjsha6zoO4CnHzg/rewc69eDN9ZxZB6QIih0AICUVZmfoN+5erz2n+/U8xx8AP+V/PH1M+Zkhff6Ota6jAEgQih0AIGV9fHed1pTm6Y+fPKqZaMx1HMATDrQP6unDF/TQ7WtUlJvpOg6ABKHYAQBSVkYwoP9y70ad7BnRo/s7XccBPOHPnjqmFXmZ+sVbV7uOAiCBKHYAgJT2noYKNVQV6n/9+ylFYzxrh/T2yslevXyyV7/yrnXKzwq5jgMggSh2AICUZozRF961Tmd6R/XD5m7XcQBnrLX68tPHVBnO1idvqHMdB0CCUewAACnvPQ0rta48X3/1/EnFWLVDmnrxRK/2tw3q1+9ar+yMoOs4ABKMYgcASHmBgNGvvmutjl0Y1rNHLriOAzjx8IuntLIwWx/eWeM6CoBlQLEDAKSFD2ytUl1Jrr76wknOtUPaae6I6Mcn+/SZW+uVGeLbPyAV8V82ACAthIIBff7OtTrUEdFLJ3pdxwGS6v+8eEoFWSF9fDfP1gGpimIHAEgbD+ysVmU4W199/qTrKEDStPWN6Ynmbn3ixjoVZGe4jgNgmVDsAABpIysU1EO3r9HrZ/u1r23AdRwgKf725dMKBow+cwvn1gGpjGIHAEgrP9dUq4KskB555azrKMCy6x+d0rf2tutD26tVUZjtOg6AZUSxAwCklbyskD7SVKMnmrvVMzzhOg6wrL756llNTMf00O1rXEcBsMwodgCAtPMLN9VrOmr1T6+1uY4CLJuJ6ai++eo53b2pXOsrClzHAbDMKHYAgLSzujRPd15Xpn98rU1TMzHXcYBl8WRLt/pHp/SLPFsHpAWKHQAgLT14c70uDk/qyZZu11GAZfGPe9q0pjRPN69d4ToKgCSg2AEA0tId68u0ujSPISpISUfPD2nvuQF94oY6GWNcxwGQBBQ7AEBaCgSMPnXjKu1rG1RzR8R1HCCh/um1NmWGAvrwzhrXUQAkCcUOAJC2PtJUo9zMoL7Bqh1SyOjkjB7d16n3b6lUcV6m6zgAkoRiBwBIW4XZGfrwzhr94FCXImPTruMACfGDg10amZzRJ2+scx0FQBJR7AAAae3nr6/V1ExMjx3sdB0FSIh/fK1NG1cWaGddsesoAJKIYgcASGuN1WFtqizUt9/scB0FiNvB9kE1d0b0SYamAGmHYgcASHsf3VWjQx0RHTs/7DoKEJd/fO2ccjOD+tCOatdRACQZxQ4AkPY+tKNaGUGjb+9tdx0FuGZjUzN6/FC3PrC1SgXZGa7jAEgyih0AIO2V5GXqro0V+t6BTk1HY67jANfk6dYLGpuK6oGdrNYB6YhiBwCApI/sqlHvyJReONrjOgpwTR7d36nqohxdX1/iOgoAByh2AABIuvO6MpXmZzFEBb7UMzShl09c1P07qhUIMDQFSEdxFTtjTIkx5hljzIm5t++Yq2uM2W6MedUY02qMOWSM+fl47gkAwHIIBQN6YGe1Xjjao96RSddxgCV57GCXYla6n22YQNqKd8Xui5Kes9aul/Tc3MdvNybpF6y1DZLulfQXxpiiOO8LAEDCfXRXjWZiVt/bz5l28JdH93VqW01Ya8vyXUcB4Ei8xe4+SY/Mvf+IpA+9/QJr7XFr7Ym597sk9Ugqi/O+AAAk3PqKAm2rLdJ32I4JHzl2fliHu4c44gBIc/EWuwprbbckzb0tv9LFxpjdkjIlnYrzvgAALIv7tlXp6PlhnewZcR0FWJTv7u9UMGD0gW1VrqMAcOiqxc4Y86wxpmWBX/ct5UbGmEpJfy/pF621C86SNsY8ZIzZa4zZe/HixaW8PAAACfG+rZUyRvrhoW7XUYCrisWsvn+gU3dsmB3+AyB9XbXYWWvvttY2LvDr+5IuzBW2+eK24IxoY0yhpB9K+h1r7Z4r3Otha22TtbaprIzdmgCA5KsozNb19SV6/FCX6yjAVe053afuyITuZxsmkPbi3Yr5mKQH595/UNL3336BMSZT0nclfdNa++047wcAwLL7wNZKnegZ0bHzw66jAFf0vQOdys8K6Z7NFa6jAHAs3mL3JUn3GGNOSLpn7mMZY5qMMV+bu+bnJN0u6dPGmANzv7bHeV8AAJbNvY2VChixagdPm47G9FTrBf3M5gplZwRdxwHgWCiev9ha2yfprgU+v1fS5+be/wdJ/xDPfQAASKaygizduGaFHj/Urd+8Z4OM4cBneM+rp/oUGZ/WvY0rXUcB4AHxrtgBAJCS3r+1Smd6R9XaNeQ6CrCgJ1vOKy8zqNs3MJcAAMUOAIAF3du4UsGA0eNMx4QHRWNWT7ee17s2lrMNE4Akih0AAAsqycvULetK9cPmLllrXccBfsrrZ/rVNzqln91S6ToKAI+g2AEAcBnv31qp9v5xHeqIuI4C/JQnW7qVnRHQndexDRPALIodAACX8Z7NK5URNEzHhKfEYlb/1nJed24oV25mXHPwAKQQih0AAJcRzs3QzWtL9fThC2zHhGfsaxtQz/Ck3ruFaZgAfoJiBwDAFdyzuULn+sZ0smfEdRRAkvRE83llBgN698Zy11EAeAjFDgCAK7h7U4Uk6enDFxwnASRrrf6tpVu3byhVQXaG6zgAPIRiBwDAFawMZ2trTVjPUOzgAQc7IuqKTOi9jUzDBPDTKHYAAFzFPZsqdKB9UD3DE66jIM091XpeoYC5tJIMAPModgAAXMU9DbPfRD93pMdxEqS754/06Pr6EoVz2YYJ4KdR7AAAuIrrKgpUU5zDdkw41d4/pmMXhnXXJoamAHgnih0AAFdhjNE9myv08slejU7OuI6DNPXCsdkVY6ZhAlgIxQ4AgEW4Z3OFpmZieulEr+soSFPPHenR6tI8rSnLdx0FgAdR7AAAWITr60tUmB1iOyacGJ2c0aun+litA3BZFDsAABYhY+5A6OePXtBMNOY6DtLMj0/2aioa010UOwCXQbEDAGCR7tm8UgNj09rXNug6CtLM80d7VJAVUlN9iesoADyKYgcAwCLdvqFUGUGj546wHRPJY63V80d7dPuGMmWG+NYNwML40wEAgEUqyM7QrlXF+tHxi66jII20dg2pZ3iS5+sAXBHFDgCAJbhjQ7mOnh/WhaEJ11GQJp470iNjpDuvK3MdBYCHUewAAFiCOzbMfnP9Iqt2SJLnj17QjtoircjPch0FgIdR7AAAWIJNlQUqK8hiOyaSomd4Qgc7IrprU4XrKAA8jmIHAMASGGN0+/oyvXSiV9GYdR0HKe5Hx2Z/gMA2TABXQ7EDAGCJ7riuTJHxaR3s4NgDLK+XT/aqND9LmysLXUcB4HEUOwAAlui2daUyhufssLxiMauXT/TqtvWlMsa4jgPA4yh2AAAsUXFeprbVFPGcHZbVkfND6hud0q3rSl1HAeADFDsAAK7BHRvKdLB9UAOjU66jIEW9dKJXknTbeoodgKuj2AEAcA3uuK5MMTv7DBSwHF4+0avrKgpUXpjtOgoAH6DYAQBwDbbVFCmck8F2TCyLiemoXj/bz2odgEWj2AEAcA2CAaNb15fqxeMXZS3HHiCxXj/Tr6mZmG6l2AFYJIodAADX6I4NZeoZntSR7mHXUZBiXjpxUZnBgG5YvcJ1FAA+QbEDAOAazW+Te+UUz9khsV460aum+mLlZAZdRwHgExQ7AACuUWU4R2tK8/TKqT7XUZBCeoYmdPT8sG5bX+Y6CgAfodgBABCHm9au0Gun+zQdjbmOghQxP2mVwSkAloJiBwBAHG5eW6rRqaiaOyOuoyBFvHyiVyV5mdpcWeg6CgAfodgBABCHG9eUSJJeZTsmEsBaq5dO9uqWdaUKBIzrOAB8hGIHAEAcVuRnaVNlIQNUkBDHL4zo4vCkblvHNkwAS0OxAwAgTjevXaG9Zwc0MR11HQU+t+f07MrvTWs55gDA0lDsAACI081rV2hyJqZ9bQOuo8Dn9pzuU3VRjmpLcl1HAeAzFDsAAOK0e3WJggHDc3aISyxmted0n25cw2odgKWj2AEAEKeC7AxtqQ5znh3icrxnWANj05cG8gDAUlDsAABIgFvWrdDB9kGNTM64jgKf2jP3gwFW7ABcC4odAAAJcPPaUs3ErN442+86Cnzq1dN9qinm+ToA14ZiBwBAAuxaVazMYIDn7HBNYjGr1870s1oH4JpR7AAASIDsjKB2ririPDtck6PnhzU4Nq2bKHYArhHFDgCABLl5balau4Y0ODblOgp8Zv78uhsYnALgGlHsAABIkBtWl8haae9ZzrPD0uw53afakhzVFPN8HYBrQ7EDACBBttUWKTMY0OsMUMESzD9fxzZMAPGg2AEAkCDZGUFtry3S62codli8I+eHFBmfZnAKgLhQ7AAASKDrVxerpTOiUc6zwyLtOT37g4AbKHYA4kCxAwAggXavXqGZmNX+tkHXUeATr57qU11JrqqLclxHAeBjFDsAABJoZ12RAkY8Z4dFicWsXj/TpxuZhgkgThQ7AAASqCA7Qw1VYb1+hoPKcXXHe4Y1NDGj3avZhgkgPhQ7AAAS7Pr6Eu1vG9TkTNR1FHjcG3NHY1xfX+w4CQC/o9gBAJBgu1eXaHImppbOiOso8Li9Z/tVVpCluhLOrwMQH4odAAAJNr/68hrHHuAq9p4d0O76EhljXEcB4HMUOwAAEmxFfpbWlefrDYodrqBzcFydg+NqYhsmgASg2AEAsAx2ry7R3rMDisas6yjwqL1zk1Ovr2ciJoD4UewAAFgGu+tLNDw5oyPdQ66jwKP2nh1QXmZQG1cWuI4CIAVQ7AAAWAa7V8+uwrzBeXa4jDfO9mvnqmKFgnw7BiB+/EkCAMAyqCrKUU1xjl7nOTssIDI+rWMXhtW0im2YABKDYgcAwDLZXV+iN872y1qes8NP29c2IGs5vw5A4sRV7IwxJcaYZ4wxJ+beXvZPJ2NMoTGm0xjz1XjuCQCAX+yqL1bvyJTa+sdcR4HH7D3br2DAaHtdkesoAFJEvCt2X5T0nLV2vaTn5j6+nD+Q9KM47wcAgG/sWjX78843zw04TgKveePsgBqrCpWbGXIdBUCKiLfY3Sfpkbn3H5H0oYUuMsbsklQh6ek47wcAgG+sLy9QQVZI+9oodviJyZmoDrYPcswBgISKt9hVWGu7JWnubfnbLzDGBCT9D0m/Hee9AADwlfmtdm+eG3QdBR7S0jmkyZmYmih2ABLoqsXOGPOsMaZlgV/3LfIevyLpCWtt+yLu9ZAxZq8xZu/FixcX+fIAAHjXzrpiHTs/pOGJaddR4BHzB5M3MTgFQAJddWO3tfbuy33NGHPBGFNpre02xlRK6lngspsk3WaM+RVJ+ZIyjTEj1tp3PI9nrX1Y0sOS1NTUxAgxAIDv7VpVrJiVDrZHdOv6Utdx4AFvnB3QmtI8leZnuY4CIIXEuxXzMUkPzr3/oKTvv/0Ca+0nrbV11tp6Sb8l6ZsLlToAAFLR9roiGcMAFcyy1mpf24B2rmK1DkBixVvsviTpHmPMCUn3zH0sY0yTMeZr8YYDAMDvCrMzdF1Fgd5kgAoktfWPqX90SjvrKHYAEiuuGbvW2j5Jdy3w+b2SPrfA578h6Rvx3BMAAL/ZuapYPzjYpVjMKhAwruPAofkJqTtXcX4dgMSKd8UOAABcxc66Yg1PzOhr+G+XAAAgAElEQVTkxRHXUeDYvnODys8KaX15gesoAFIMxQ4AgGXGQeWYt799QNtqwwqycgsgwSh2AAAss/oVuSrJy6TYpbmxqRkd6R7WjlqerwOQeBQ7AACWmTFGO+uKtY9il9aaOyKKxizP1wFYFhQ7AACSYNeqYp3uHVX/6JTrKHBkX9ugJLFiB2BZUOwAAEiC+efs9nPsQdra1zag1aV5Ks7LdB0FQAqi2AEAkARba8IKBQzP2aUpa632tw1qRx3bMAEsD4odAABJkJ0RVENVIcUuTXUMjKt3ZFI7OJgcwDKh2AEAkCTba4vU3Dk7QAPp5dLB5KzYAVgmFDsAAJJke12RxqaiOn5h2HUUJNn+tkHlZgZ1XQUHkwNYHhQ7AACSZPvcNMQD7YOOkyDZ9rUNzD5nGeRbLwDLgz9dAABIkvoVuSrKzdCBNopdOpmYjupw15B28nwdgGVEsQMAIEmMMdpWU8SKXZpp7oxoJmYZnAJgWVHsAABIou21RTreM6yRyRnXUZAk++YmoXLUAYDlRLEDACCJttcVyVqpuSPiOgqS5ED7oOpKclWan+U6CoAURrEDACCJttfMrtqwHTN9HGwf1LZaVusALC+KHQAASVScl6n6Fbk60M5B5emgZ2hCXZEJbasJu44CIMVR7AAASLLttQxQSRcH57bcbmfFDsAyo9gBAJBk22qLdGFoUt2RcddRsMwOtg8qGDBqqGLFDsDyotgBAJBk86s3nGeX+g52DOq6igLlZAZdRwGQ4ih2AAAk2eaqQmUGA2zHTHGxmGVwCoCkodgBAJBkWaGgNlUVaj/FLqWd7RvV0MSMtteyDRPA8qPYAQDgwI7aIjV3RDQTjbmOgmVysGO2uLNiByAZKHYAADiwvbZI49NRHb8w4joKlsnB9ohyM4NaX17gOgqANECxAwDAgUsDVNiOmbIOtA+qsTqsYMC4jgIgDVDsAABwYNWKXBXnZuggxS4lTc3EdLhriPPrACQNxQ4AAAeMMdpSU3TpOSyklqPnhzQVjWlbDcUOQHJQ7AAAcGRrdVgnekY0PhV1HQUJNr8Su42JmACShGIHAIAjW2vCisasDndHXEdBgh1oj6g0P1PVRTmuowBIExQ7AAAc2Tq3Te9QB8Uu1RzsGNS2miIZw+AUAMlBsQMAwJGV4WyVF2RR7FLM0MS0Tl0c4fw6AElFsQMAwKGtDFBJOS0dEVnLweQAkotiBwCAQ1trwjp9cVTDE9OuoyBBDs6twG6tZnAKgOSh2AEA4NDWmtlv/ps72Y6ZKpo7B1VXkqvivEzXUQCkEYodAAAOMUAl9RzqiGhLDat1AJKLYgcAgEMleZmqKc5RM8UuJfSPTqljYJxtmACSjmIHAIBj2xigkjLmt9SyYgcg2Sh2AAA4tqUmrI6BcfWPTrmOgjg1zxX0RlbsACQZxQ4AAMfmB6gcYtXO9w51RLSmNE+F2RmuowBIMxQ7AAAc21I9X+x4zs7vmjsZnALADYodAACOFWRnaE1ZHsXO53qGJ9QdmbhU1AEgmSh2AAB4wLaaIrZi+lzL/OAUih0AByh2AAB4wJbqsHqGJ3U+MuE6Cq5Rc8eQjJEaKHYAHKDYAQDgAdtqZ8sAxx74V3PnoNaW5Ss/K+Q6CoA0RLEDAMADNleGFTBSayfP2fnVoY4IB5MDcIZiBwCAB+RkBrWuPP/SAdfwlwtDE+oZnmQiJgBnKHYAAHhEY3VYLV1DrmPgGsxPNN1KsQPgCMUOAACPaKwK6+LwpC4MMUDFb5o7BhUws1tqAcAFih0AAB4xv42vmfPsfOdQZ0QbKgqUkxl0HQVAmqLYAQDgEZsrC2WM1NJFsfMTa62aOyKcXwfAKYodAAAekZcV0prSvEsHXcMfuiIT6hud4vk6AE5R7AAA8JAt1WG1dDJAxU/mt842smIHwCGKHQAAHtJYHdb5oQldHJ50HQWL1NoVUTBgtKmy0HUUAGmMYgcAgIfMr/rwnJ1/tHRGtL48X9kZDE4B4A7FDgAAD2moml31aWEypm+0dA2poYptmADcotgBAOAhBdkZWl2ap2YGqPjChblts43VbMME4BbFDgAAj2msDqu1iwEqfjA/wZTBKQBco9gBAOAxjVWF6hwcV//olOsouIqWziEZIwanAHCOYgcAgMfMH3TNeXbe19IV0erSPOVnhVxHAZDmKHYAAHjM/CAOnrPzvtbOiBoZnALAAyh2AAB4TDg3Q3UluWrlyANP6xuZVFdkgsEpADyBYgcAgAdtqQ6zYudx8wNuWLED4AVxFTtjTIkx5hljzIm5t8WXua7OGPO0MeaIMeawMaY+nvsCAJDqGqoL1d4/rsExBqh41fwh8g1MxATgAfGu2H1R0nPW2vWSnpv7eCHflPRla+0mSbsl9cR5XwAAUtr8ABWOPfCu1s4h1ZXkKpyT4ToKAMRd7O6T9Mjc+49I+tDbLzDGbJYUstY+I0nW2hFr7Vic9wUAIKU1MkDF81q6IjxfB8Az4i12Fdbabkmae1u+wDUbJA0aYx41xuw3xnzZGBOM874AAKS04rxMVRflcOSBR0XGp3Wub+zSBFMAcO2qh64YY56VtHKBL/23JdzjNkk7JLVJ+pakT0v62wXu9ZCkhySprq5ukS8PAEBqaqwupNh51PzE0kaerwPgEVddsbPW3m2tbVzg1/clXTDGVErS3NuFnp3rkLTfWnvaWjsj6XuSdl7mXg9ba5ustU1lZWXX/ncFAEAK2FId1tm+MQ1NTLuOgrdp7Zx99rGhiq2YALwh3q2Yj0l6cO79ByV9f4Fr3pBUbIyZb2rvlnQ4zvsCAJDy5leD5ksEvKOlK6LKcLZK87NcRwEASfEXuy9JuscYc0LSPXMfyxjTZIz5miRZa6OSfkvSc8aYZklG0t/EeV8AAFLefLFjO6b3tHRGeL4OgKdc9Rm7K7HW9km6a4HP75X0ubd8/IykrfHcCwCAdFOan6XKcPal89LgDaOTMzrdO6oPbKtyHQUALol3xQ4AACyjhqowRx54zJHuIVn7k7MGAcALKHYAAHjYluqwzvSOamRyxnUUzJnfGstETABeQrEDAMDDttQUylrpcBcDVLyipWtIpflZKi9gcAoA76DYAQDgYY1VDFDxmpbOiBqrC2WMcR0FAC6h2AEA4GHlhdkqL8ii2HnExHRUJ3pGLhVuAPAKih0AAB7XWM0AFa84dn5Y0ZhVYzUHkwPwFoodAAAe11gd1qmLIxqbYoCKa/NHT3CGHQCvodgBAOBxjVWFitnZMftwq6UzonBOhmqKc1xHAYCfQrEDAMDjttTMD1Ch2LnW0jnE4BQAnkSxAwDA41YWZqs0P5Pn7Bybmonp2PlhBqcA8CSKHQAAHmeMUUNVmMmYjp3oGdZUNMbB5AA8iWIHAIAPbKkO60TPiCamo66jpK3Wua2wFDsAXkSxAwDABxqrCxWNWR09P+w6Stpq6YooPyukVSW5rqMAwDtQ7AAA8IH58fpsx3SnpTOizVWFCgQYnALAeyh2AAD4QE1xjsI5GWrtYjKmC9GY1eHuIQanAPAsih0AAD5gjFFjdaFau1ixc+H0xRFNTMfUWF3oOgoALIhiBwCATzRWhXW0e1jT0ZjrKGmnZa5QMzgFgFdR7AAA8ImG6rCmojGduDDiOkraaekcUnZGQGtK81xHAYAFUewAAPCJxqrZbYAtbMdMuubOiDZVFioU5FsnAN7En04AAPhE/Yo85WUG1cpkzKSKxawOdw1pC9swAXgYxQ4AAJ8IBIwaqsJqYTJmUp3rH9PI5AwTMQF4GsUOAAAfaagu1OGuIUVj1nWUtDF/dmADEzEBeBjFDgAAH2moCmt8OqozvQxQSZaWrogygwGtLy9wHQUALotiBwCAj8yfo9bSyXbMZGntHNJ1KwuUGeLbJgDexZ9QAAD4yLqyfGWFApe2B2J5WWvV0hXhYHIAnkexAwDAR0LBgDZWFqqVASpJ0Tk4rsGxaTUwOAWAx1HsAADwmcaqQrV0RWQtA1SW2/yW10aOOgDgcRQ7AAB8prE6rOGJGbX3j7uOkvJauyIKBow2rmRwCgBvo9gBAOAz8+eptXTxnN1ya+mMaH15vrIzgq6jAMAVUewAAPCZDSvzFQoYBqgsM2utmjuH2IYJwBcodgAA+ExWKKgNFQVqYYDKsuoZnlTvyKQaq5iICcD7KHYAAPhQY3WhWjsZoLKc5ldEWbED4AcUOwAAfKixOqy+0SmdH5pwHSVltXQOyRhpUyUrdgC8j2IHAIAPzZ+rNj+OH4nX0hXRmtI85WWFXEcBgKui2AEA4EObKgsUMGKAyjJq7YywDROAb1DsAADwodzMkNaW5auVIw+WRd/IpLoiE5eOlgAAr6PYAQDgU43VYbUyGXNZzP++NlTzfB0Af6DYAQDgUw1VheqOTKh3ZNJ1lJQzf/h7Ayt2AHyCYgcAgE/Nlw5W7RKvpTOiupJchXMyXEcBgEWh2AEA4FOb5w7OZoBK4jV3RrSFwSkAfIRiBwCAT4VzMrRqRS4DVBJscGxK7f3jTMQE4CsUOwAAfKyxKsxZdgk2//vJih0AP6HYAQDgYw3VhWrrH1NkbNp1lJTRPLe1tZGJmAB8hGIHAICPzZ+z1trNdsxEaemMqLYkR0W5ma6jAMCiUewAAPCxhrkBKq1sx0wYBqcA8COKHQAAPrYiP0tV4exL564hPoNjU2rrH9OW6iLXUQBgSSh2AAD4XEN1mCMPEoTBKQD8imIHAIDPNVaFdbp3VGNTM66j+B6DUwD4FcUOAACfa6wulLXSkW6es4sXg1MA+BXFDgAAn5s/SJvz7OLH4BQAfkWxAwDA58oLslSan3VpGyGuTWRsWm39Y5eKMgD4CcUOAACfM8aosbqQASpxmp8syoodAD+i2AEAkAK2VId1omdEE9NR11F861DH3OCUKoodAP+h2AEAkAIaq8OKxqwOM0Dlms0PTinOY3AKAP+h2AEAkAK2XBqgwnbMa8XgFAB+RrEDACAFVIazVZqfqeYOit21YHAKAL+j2AEAkAJmB6iEmYx5jRicAsDvKHYAAKQIBqhcu/lCzOAUAH5FsQMAIEUwQOXaNXdGVFPM4BQA/kWxAwAgRTBA5do1dzA4BYC/UewAAEgRDFC5NgxOAZAKKHYAAKQIBqhcm/nBKVtrKHYA/CuuYmeMKTHGPGOMOTH3tvgy1/2pMabVGHPEGPMVY4yJ574AAGBhDFBZOganAEgF8a7YfVHSc9ba9ZKem/v4pxhjbpZ0i6StkholXS/pjjjvCwAAFsAAlaVjcAqAVBBvsbtP0iNz7z8i6UMLXGMlZUvKlJQlKUPShTjvCwAAFsAAlaVr6WRwCgD/i7fYVVhruyVp7m352y+w1r4q6QVJ3XO/nrLWHonzvgAAYAHzA1QOMUBlUSJj0zrXx+AUAP4XutoFxphnJa1c4Ev/bTE3MMask7RJUs3cp54xxtxurX1xgWsfkvSQJNXV1S3m5QEAwFvMD1BhxW5x5gensGIHwO+uWuystXdf7mvGmAvGmEprbbcxplJSzwKX3S9pj7V2ZO6veVLSjZLeUeystQ9LeliSmpqa7OL+FgAAwFttqQ7rpRO9mpiOKjsj6DqOp80PTqHYAfC7eLdiPibpwbn3H5T0/QWuaZN0hzEmZIzJ0OzgFLZiAgCwTBigsngMTgGQKuItdl+SdI8x5oSke+Y+ljGmyRjztblrviPplKRmSQclHbTW/iDO+wIAgMuYX33ioPKrY3AKgFRx1a2YV2Kt7ZN01wKf3yvpc3PvRyX9Ujz3AQAAi8cAlcWJjM8OTvm5plrXUQAgbvGu2AEAAI8xxmhrTZEOdQy6juJprTxfByCFUOwAAEhB22qKdPLiiEYmZ1xH8SwGpwBIJRQ7AABS0NbasKzlObsrOdQZUXURg1MApAaKHQAAKWhbTZEksR3zChicAiCVUOwAAEhBJXmZqinOYYDKZcwPTtlSQ7EDkBoodgAApKhtNUU6yIrdghicAiDVUOwAAEhR22rD6hgYV9/IpOsonnOwg2IHILVQ7AAASFFbLz1nx3bMtzvYPqhVK3IZnAIgZVDsAABIUY3VYRkjtmMu4GDH4KUBMwCQCih2AACkqPyskNaV5bNi9zYXhibUHZnQtlqKHYDUQbEDACCFba0p0qGOQVlrXUfxjIPtsyuY2yl2AFIIxQ4AgBS2rTas3pEpdQ6Ou47iGQc7BhUKGDVUFbqOAgAJQ7EDACCFbWOAyjscaB/UxsoCZWcEXUcBgISh2AEAkMI2VhYoI2gYoDInFrM61B5hcAqAlEOxAwAghWWFgtpUWahD7azYSdLp3lENT84wOAVAyqHYAQCQ4rbWhNXcGVEsxgCV+cEpOyh2AFIMxQ4AgBS3taZII5MzOt074jqKcwc7BpWfFdKasnzXUQAgoSh2AACkuPmx/gfYjqkD7YPaUh1WMGBcRwGAhKLYAQCQ4taV5asgK6T9bQOuozg1MR3Vke4hnq8DkJIodgAApLhAwGhbbZH2t6X3ZMwj3UOajlptrw27jgIACUexAwAgDeyoK9LR80Mam5pxHcWZ+cEp22uLHScBgMSj2AEAkAZ21hUrZqWDafyc3cGOiCoKs7QynO06CgAkHMUOAIA0MD9AZX97+j5nd6B9kIPJAaQsih0AAGmgOC9Tq0vz0vY5u8jYtM70jjI4BUDKotgBAJAmdtTNDlCxNv0OKp9fqeRgcgCpimIHAECa2FFXrN6RSXUMjLuOknT7zg0oYMSKHYCURbEDACBNzK9W7UvD8+zebBvQpspC5WWFXEcBgGVBsQMAIE1sXFmgnIxg2j1nF41ZHWgb1K5VHHMAIHVR7AAASBOhYEBba8Lan2YrdsfOD2t0KkqxA5DSKHYAAKSRHXXFau0a0sR01HWUpHlzrsjurKPYAUhdFDsAANLIjroizcSsWrvS56DyfecGVFaQpZriHNdRAGDZUOwAAEgjO+rmDipPo+fs3jw3oF11xTLGuI4CAMuGYgcAQBopL8hWTXFO2kzGvDg8qbb+MZ6vA5DyKHYAAKSZHXXFabNiN19gd1LsAKQ4ih0AAGlmR22RuiMT6o6k/kHl+84NKDMYUGN1oesoALCsKHYAAKSZpvrZ1as3zqb+dsw3zw1oS01YWaGg6ygAsKwodgAApJnNlYXKywzqjTP9rqMsq8mZqA51Rni+DkBaoNgBAJBmQsGAdq4q1htnU7vYtXYNaWomxvl1ANICxQ4AgDTUtKpExy4MKzI27TrKstl3bn5wSpHjJACw/Ch2AACkoetXF8taae+51F2129c2oNqSHJUXZLuOAgDLjmIHAEAa2lFbrIyg0espuh3TWnvpYHIASAcUOwAA0lBOZlCN1WHtTdHJmG39Y7owNKld9SWuowBAUlDsAABIU7vrS3SoY1AT01HXURLutdOzK5E3rqbYAUgPFDsAANLU9fUlmo5aHWgfdB0l4fac6dOKvEytK893HQUAkoJiBwBAmrp0UHkKnmf32ul+3bCmRMYY11EAICkodgAApKmi3ExdV1GQcgNU2vvH1Dk4rhtWr3AdBQCShmIHAEAau351sfadG9BMNOY6SsK8NrcCecManq8DkD4odgAApLHr60s0OhXVke5h11ES5rXTfSrKzdCG8gLXUQAgaSh2AACksd1zUyNTaTvma2f6tbu+RIEAz9cBSB8UOwAA0lhlOEfVRTnamyLFrjsyrrb+Md2whufrAKQXih0AAGlu9+oSvXG2X9Za11HiNn9+3Q2cXwcgzVDsAABIczetWaHekSkdvzDiOkrcXjvTp8LskDZVFrqOAgBJRbEDACDN3bK+VJL045O9jpPEb8/pfu1eXaIgz9cBSDMUOwAA0lx1UY7qV+T6vtj1DE3oTO8o59cBSEsUOwAAoFvWlWrP6T5N+/g8uz2cXwcgjVHsAACAbllXqtGpqA51DLqOcs1eO92n/KyQNvN8HYA0RLEDAAC6ac0KGSO9fKLPdZRr9urpPjXVFysU5NsbAOmHP/kAAICK8zLVWBXWj0/58zm7zsFxnb44qlvXlbqOAgBOUOwAAIAk6eZ1K7S/bUBjUzOuoyzZi8cvSpLu2FDmOAkAuEGxAwAAkqRb15VqOmr1+twQEj958fhFVYazta4833UUAHCCYgcAACRJ19eXKDMU8N2xBzPRmF4+2avb15fJGM6vA5CeKHYAAECSlJ0R1K66Yr180l8DVA52DGp4Yka3sw0TQBqLq9gZYz5qjGk1xsSMMU1XuO5eY8wxY8xJY8wX47knAABYPreuL9WR7iH1jUy6jrJoPzreq4ARg1MApLV4V+xaJD0g6cXLXWCMCUr6K0nvlbRZ0seNMZvjvC8AAFgGt8yVo1dO+WfV7sXjF7Wttkjh3AzXUQDAmbiKnbX2iLX22FUu2y3ppLX2tLV2StK/SLovnvsCAIDlsaU6rILskF7xybEHg2NTOtQxqNvXsw0TQHpLxjN21ZLa3/Jxx9znAACAxwQDRjevXaEfHbsoa63rOFf18slexax4vg5A2rtqsTPGPGuMaVng12JX3RYaT7Xg/ymMMQ8ZY/YaY/ZevHhxkS8PAAAS6a6NFeqKTOhI97DrKFf14vGLKsgOaVtN2HUUAHAqdLULrLV3x3mPDkm1b/m4RlLXZe71sKSHJampqcn7PyYEACAFvWtjuYyRnj1yQZurCl3HuSxrrV483qtb15UqFGTQN4D0low/Bd+QtN4Ys9oYkynpY5IeS8J9AQDANSgryNL22iI9e+SC6yhXdKJnROeHJtiGCQCK/7iD+40xHZJukvRDY8xTc5+vMsY8IUnW2hlJX5D0lKQjkv7VWtsaX2wAALCc7t5UoUMdEV0YmnAd5bJePD772AbFDgDin4r5XWttjbU2y1pbYa19z9znu6y1P/uW656w1m6w1q611v5RvKEBAMDyuntThSTpuSM9jpNc3vNHe7SuPF/VRTmuowCAc2xIBwAA77ChIl+1JTl6zqPbMftHp/TamX69p6HCdRQA8ASKHQAAeAdjjO7eVKGXT/ZqfCrqOs47PHv4gqIxq/c2VrqOAgCeQLEDAAALuntThSZnYnr5pPcOK3+ypVs1xTlq8PDUTgBIJoodAABY0O7VJSrIDunZw97ajjk0Ma0fn+zTvQ0rZcxCx+UCQPqh2AEAgAVlBAO6Y0OZnjvao1jMO8fLvnC0R1PRmN67ZaXrKADgGRQ7AABwWfdsrlDvyKQOdgy6jnLJk83nVV6QpR21xa6jAIBnUOwAAMBl3bmhXMGA0TMe2Y45NjWjfz/eo/c0rFQgwDZMAJhHsQMAAJcVzs3QzWtX6LGDXZ7Yjvni8YuamI7pvY1swwSAt6LYAQCAK7p/R7U6Bsa199yA6yh6suW8inMztHt1iesoAOApFDsAAHBF72lYqdzMoL67v8NpjsmZqJ4/0qN7NlcoFORbGAB4K/5UBAAAV5SXFdK9DSv1+KFuTUy7O6z8lZN9Gp6c4VByAFgAxQ4AAFzV/TurNTwxo+eP9jjL8Oj+ThVmh3TzuhXOMgCAV1HsAADAVd28tlQVhVl6dF+nk/v3j07pqZbzemBnjbJCQScZAMDLKHYAAOCqggGj+7ZX69+P9ahvZDLp9390X4emojF9fHdd0u8NAH5AsQMAAIty/45qzcSsHj/UndT7Wmv1T6+3aWddka5bWZDUewOAX1DsAADAomyqLNSmykI9uj+52zHfODug0xdH9TFW6wDgsih2AABg0R7YUa2D7YM6dXEkaff859fbVJAV0vu3/v/t3X2QVeV9wPHvD1aWBLsKxRfAmPCmEkkKydZm0hfEEDXOVLQ1KXZiTGumQU2aNu1MtaSZTqedNGkyznTSaULek07RSqvSqkM06rQ6ot3Wl8UYYIFExA0ECGpEYIGnf9yzmVO4lz13d+8ezuX7mTmz5zznOec+93d/PN7fvecevRumJDViYSdJkgpbumA648cFq554YUweb+++g9zb28/ShdN544SOMXlMSaoiCztJklTYmV0T+c23T2PVky/w8r6Blj/eXU9t5+ChIyz7ZS/DlKTjsbCTJElN+eii2bx28DDfWffDlj5OSonbn9zG2885jfkzTmvpY0lS1VnYSZKkpsyb1sXF55/BNx77IfsHDrfscXp+9FM27HjVb+skqQALO0mS1LQbF81m92sHubNnW0vOn1LiC9/dwNRTJ7B0wfSWPIYktRMLO0mS1LSLZk5h4bmns/K/tnDo8JFRP/9jfbtZt2UPNy+ew6ROb5oiSUOxsJMkSU2LCJYvms22Pa9z3/ofj+q5U0r83dofMOP0N/C7v+JlmJJUhIWdJEkalvfOO4vZZ0ziHx/ZTEpp1M679rkdPPPiy3xiyVw6O8aP2nklqZ1Z2EmSpGEZNy746KLZPN//Cvf1js63doeP1H5bN/uMSfzWwhmjck5JOhlY2EmSpGG7euEM5s/o4tP3rGfPawdHfL67n9rOpp0/408uPZ+O8b5NkaSinDElSdKwnTJ+HJ9//y/xyv4BPn3P+hGda//AYW57cCPzZ3Rx+YVnj9IIJenkYGEnSZJG5IKzu/jDS+byH8/2c39v/7DOkVLi1n/rZfve1/nz981j3LgY5VFKUnuzsJMkSSO2/OLZzJ/RxafuHt4lmV97dCt3PbWdTy45j3fPmdqCEUpSe7OwkyRJI5a/JPMv7lnf1F0yH+vbxWfu/wGXXXgWNy+e08JRSlL7srCTJEmj4oKzu/ijJedx77P9/PEdT7N/4PCQx2zbs4+P/fP/MmvqJL7wgQVegilJw9RR9gAkSVL7uOni2aSU+Px3N7J19z5WXvdOzuqaWLfv45t3s+LuXg4dSaz8UDendvq2RJKGy2/sJEnSqIkIPnbJXL583TvZtONVrvziozy6aRev7h/4eZ+tu17jD77dw7VfWcf+g4dZeV03M6dOKnHUklR90cw18GOpu7s79e2bC5oAAAf4SURBVPT0lD0MSZI0TM/3v8JHvtXD9r2vA7B64l/T2TGOq/etoLNjHDctnsMNvzaTiaeML3mkknRiioj/SSl1F+nrNQ+SJKkl5k3r4r5P/DqP9e3ihT37ONB7Ibv2H+La+efy8ffM4cxfqH+JpiSpeRZ2kiSpZU57wylc8bZptY1FXwXg8hLHI0ntyt/YSZIkSVLFWdhJkiRJUsVZ2EmSJElSxVnYSZIkSVLFWdhJkiRJUsVZ2EmSJElSxVnYSZIkSVLFWdhJkiRJUsVZ2EmSJElSxVnYSZIkSVLFWdhJkiRJUsVZ2EmSJElSxVnYSZIkSVLFWdhJkiRJUsVZ2EmSJElSxVnYSZIkSVLFWdhJkiRJUsVZ2EmSJElSxUVKqewx1BURPwF+VPY46pgK7Cp7ECcx418eY18eY18eY18eY18eY18eY1+eEzX2b04pnVGk4wlb2J2oIqInpdRd9jhOVsa/PMa+PMa+PMa+PMa+PMa+PMa+PO0Qey/FlCRJkqSKs7CTJEmSpIqzsGveyrIHcJIz/uUx9uUx9uUx9uUx9uUx9uUx9uWpfOz9jZ0kSZIkVZzf2EmSJElSxVnY1RER74+I5yLiSEQ0vDtORFweERsioi8ibsm1z4yIJyJiU0TcERETxmbk1RcRUyLigSx2D0TE5Dp9FkfE07llf0Rcle37ZkRsze1bMPbPorqKxD/rdzgX4zW5dnN/mArm/oKIeDybn56NiN/J7TP3m9RoDs/t78zyuC/L67fk9t2atW+IiMvGctztoEDsPxkR38/y/HsR8ebcvrrzj4opEPsPR8RPcjH+SG7f9dkctSkirh/bkVdfgdjflov7xojYm9tn3g9TRHw9InZGxPoG+yMi/j57XZ6NiHfk9lUr51NKLkctwDzgfOARoLtBn/HAZmAWMAF4Bnhrtu9fgGXZ+peAG8t+TlVZgM8Bt2TrtwCfHaL/FGAP8MZs+5vANWU/j6ouReMP/KxBu7nfwtgD5wFzs/XpQD9werZt7jcX74ZzeK7PTcCXsvVlwB3Z+luz/p3AzOw848t+TlVZCsZ+cW5ev3Ew9tl23fnHZdRi/2Hgi3WOnQJsyf5OztYnl/2cqrIUif1R/T8OfD23bd4PP/a/AbwDWN9g/xXA/UAA7wKeyNorl/N+Y1dHSun5lNKGIbpdBPSllLaklA4CtwNLIyKAS4DVWb9vAVe1brRtZym1mEGx2F0D3J9S2tfSUZ08mo3/z5n7IzZk7FNKG1NKm7L1l4CdQKH/aamOUXcOP6pP/jVZDbwny/OlwO0ppQMppa1AX3Y+FTNk7FNKD+fm9XXAOWM8xnZVJO8buQx4IKW0J6X0U+AB4PIWjbMdNRv7a4FVYzKyNpdS+k9qXwI0shT4dqpZB5weEdOoYM5b2A3fDGBbbvvFrO0Xgb0ppUNHtauYs1JK/QDZ3zOH6L+MYye+v8m+Sr8tIjpbMcg2VjT+EyOiJyLWDV4Gi7k/Uk3lfkRcRO1T3825ZnO/uEZzeN0+WV6/TC3PixyrxpqN3w3UPk0fVG/+UTFFY//b2VyyOiLe1OSxqq9w/LJLj2cCD+WazfvWafTaVC7nO8oeQFki4kHg7Dq7VqSU7ilyijpt6Tjtyhwv9k2eZxrwNmBtrvlW4MfU3vCuBP4M+KvhjbQ9jVL8z00pvRQRs4CHIqIXeKVOP3M/Z5Rz/zvA9SmlI1mzud+cInO183xrFI5fRHwQ6AYW5ZqPmX9SSpvrHa9jFIn9vwOrUkoHImI5tW+tLyl4rBprJn7LgNUppcO5NvO+ddpmrj9pC7uU0pIRnuJF4E257XOAl4Bd1L7C7cg+4R1sV+Z4sY+IHRExLaXUn7153XmcU30AuCulNJA7d3+2eiAivgH86agMuo2MRvyzywBJKW2JiEeAhcC/Yu4f12jEPiK6gHuBT2WXjAye29xvTqM5vF6fFyOiAziN2uU8RY5VY4XiFxFLqH3osSildGCwvcH84xvcYoaMfUppd27zK8Bnc8defNSxj4z6CNtXM/PGMuDmfIN531KNXpvK5byXYg7ffwNzo3YXwAnU/hGuSbVfWz5M7bdfANcDRb4BVM0aajGDoWN3zPXn2Rviwd97XQXUvQOSGhoy/hExefAyv4iYCvwq8H1zf8SKxH4CcBe13wLcedQ+c785defwo/rkX5NrgIeyPF8DLIvaXTNnAnOBJ8do3O1gyNhHxELgy8CVKaWdufa688+Yjbz6isR+Wm7zSuD5bH0tcGn2GkwGLuX/XzGj4ysy5xAR51O7UcfjuTbzvrXWAB/K7o75LuDl7MPS6uV82XdvOREX4GpqVfoBYAewNmufDtyX63cFsJHaJyYrcu2zqP1Hvg+4E+gs+zlVZaH2+5XvAZuyv1Oy9m7gq7l+bwG2A+OOOv4hoJfam9p/Ak4t+zlVaSkSf+DdWYyfyf7ekDve3G9t7D8IDABP55YF2T5zv/mYHzOHU7t89cpsfWKWx31ZXs/KHbsiO24D8L6yn0vVlgKxfzD77+9gnq/J2hvOPy6jFvvPAM9lMX4YuCB37O9n/x76gN8r+7lUbRkq9tn2XwJ/e9Rx5v3I4r6K2l2kB6i9v78BWA4sz/YH8A/Z69JL7o74Vcv5yAYtSZIkSaooL8WUJEmSpIqzsJMkSZKkirOwkyRJkqSKs7CTJEmSpIqzsJMkSZKkirOwkyRJkqSKs7CTJEmSpIqzsJMkSZKkivs/BQGr2/uIVYgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sigma = 3.0\n", "\n", "mu = np.linspace( -1.0, 1.0, 200 )\n", "f_mu = f( mu, sigma )\n", "\n", "avg = np.mean( s )\n", "print( avg )\n", "\n", "pp.figure( figsize = ( 15, 8 ) )\n", "pp.plot( mu, f_mu )\n", "pp.plot( [ avg, avg ], [ np.min( f_mu ), np.max( f_mu ) ] )\n", "pp.show()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "x**3/3 - 3*x**2/2 + 2*x" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "def f( x ):\n", " \n", " return x**3/3 - 3*x**2/2 + 2*x\n", "\n", "X = sym.symbols( 'x' )\n", "sym.integrate( ( X - 1 ) * ( X - 2 ), X )" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsYAAAHVCAYAAADywj0dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzs3Xd81fWh//H3J3snZEMSICHspRAgCogoKlbraKt11oV06O2619b+arW3t+3tbW9r663VWqvUiaNVqQsXiqKMgKwwQggJCWQnZK9zzuf3B2hRkQQ44XvG6/l45AHnnC/nvAkQ3ueTzzDWWgEAAADBLsTpAAAAAIAvoBgDAAAAohgDAAAAkijGAAAAgCSKMQAAACCJYgwAAABIohgDAAAAkijGAAAAgCSKMQAAACBJCnPqhVNTU+3IkSOdenkAAAAEifXr1zdYa9P6u86xYjxy5EgVFRU59fIAAAAIEsaYioFcx1QKAAAAQBRjAAAAQBLFGAAAAJBEMQYAAAAkUYwBAAAASRRjAAAAQBLFGAAAAJBEMQYAAAAkUYwBAAAASRRjAAAAQBLFGAAAAJBEMQYAAAAkUYwBAAAASRRjAAAAQBLFGAAAAJBEMQYAAAAkUYwBAAAwyLr73Kpv63E6Rr8oxgAAHK9Xbj/4AeCoXtxcrdP++03tqm1zOspRhTkdAAAAv1WzxekEgF9YunavhifHKD89zukoR8WIMQAAAAbNrto2FVU066szcmSMcTrOUVGMAQAAMGiWrqtUeKjRl6dnOx2lXxRjAAAADIoel1v/2FClcyZkKDUu0uk4/aIYAwAAYFC8Vlyr5s4+XTFjuNNRBoRiDAAAgEGxdN1eZSVFa05+qtNRBoRiDAAAAK+raOzQqtJGfXVGjkJCfHvR3UcoxgAAAPC6p9ZVKsRIlxX4/qK7j1CMAQAA4FV9bo+eWV+l+WPTNTQx2uk4A0YxBgAAgFet2FGn+rYeXTHTPxbdfYRiDAAAAK9auq5SGQmRmj82zekox4RiDAAAAK+pbunS2zvrdNn0HIWF+lfV9K+0AAAA8GlPr6uSx0qXF+Q4HeWYUYwBAADgFW6P1dNFlZqTn6rhKTFOxzlmFGMAAAB4xbu76rXvQJeumOl/o8USxRgAAABe8uTavUqOjdA5EzKcjnJcKMYAAAA4YdUtXXp9W60uL8hRZFio03GOC8UYAAAAJ+zJNXtlJV09y7/2Lj4cxRgAAAAnpM/t0ZPrKjV/bLpykv1v0d1HKMYAAAA4Ia8V16q+rUfXFPrvaLFEMQYAAMAJenR1ubKHRGvemHSno5wQijEAAACOW2ldm1aXNemqWcMVGmKcjnNCKMYAAAA4bo+t3quI0BC/POnu0/otxsaYh4wxdcaYrf1cN8MY4zbGfMV78QAAAOCrOnpc+vv6Kn1hcqZS4yKdjnPCBjJivETSwqNdYIwJlfQ/kpZ7IRMAAAD8wLJN+9XW49I1hSOcjuIV/RZja+1KSU39XPZvkv4uqc4boQAAAODbrLV69IMKjcuM1/QRQ5yO4xUnPMfYGJMl6VJJ9w/g2sXGmCJjTFF9ff2JvjQAAAAc8mHlAW2rbtW1p42QMf696O4j3lh893tJP7TWuvu70Fr7gLW2wFpbkJaW5oWXBgAAgBMe+6BCcZFhuuSULKejeE2YF56jQNLSQ+8UUiV9wRjjstY+74XnBgAAgI9p6ujVi5urdcXMHMVGeqNO+oYT/p1Ya3M/+rkxZomkFynFAAAAgeuZokr1uj0Bs+juI/0WY2PMk5LOlJRqjKmSdJekcEmy1vY7rxgAAACBw+2xenR1hWbmJmtMRrzTcbyq32Jsrb1yoE9mrb3+hNIAAADAp72+rVZVzV2644LxTkfxOk6+AwAAwIA9vGqPspKitWB8htNRvI5iDAAAgAEp3t+iNXua9LXTRigsNPBqZOD9jgAAADAolqwqV3R4qK6YMdzpKIOCYgwAAIB+Nbb36IVN+/WlaVlKjAl3Os6goBgDAACgX0+u3atel0c3zB7pdJRBQzEGAADAUfW5PXp0dYXmjk5VfnpgbdF2OIoxAAAAjurlLdWqbe3RjbNz+7/Yj1GMAQAAcFQPrypXbmqs5o1JczrKoKIYAwAA4HN9uLdZGysP6LrTRigkxDgdZ1BRjAEAAPC5Hl5VrvjIMH2lIMfpKIOOYgwAAIAjqm3t1stbqnVZQY7iIsOcjjPoKMYAAAA4osdWV8htra4/faTTUU4KijEAAAA+o6vXrcfX7NXZ4zI0PCXG6TgnBcUYAAAAn/H3DVVq6ujVzXMDe4u2w1GMAQAA8Aluj9WD75ZpanaiZuYmOx3npKEYAwAA4BNe31ar8sZOLT5jlIwJ7C3aDkcxBgAAwCf85d0y5SRH67yJGU5HOakoxgAAAPjY+oomra9o1qI5eQoLDa6qGFy/WwAAABzVAyvLlBgdrssKsp2OctJRjAEAACBJ2tPQode21erawhGKiQj8Az0+jWIMAAAASdKD75YpPCRE1wXJgR6fRjEGAACAGtt79Oz6Kn1pWpbS4iOdjuMIijEAAAD06OoK9bg8WhREB3p8GsUYAAAgyHX1uvXIBxVaMD5d+enxTsdxDMUYAAAgyP3r+Oc8p6M4imIMAAAQxIL1+OcjoRgDAAAEsVe2Vqu8sVNfnxdcxz8fCcUYAAAgSFlrde+K3RqVFquFEzOdjuM4ijEAAECQWrGzTturW/XNM/MVEhLco8USxRgAACAoWWv1x7dKlZUUrYtPGeZ0HJ9AMQYAAAhCq8uatGHvAX1jXp7CQ6mEEsUYAAAgKN27olSpcZG6rCDH6Sg+g2IMAAAQZDZWHtB7pQ26eW6uosJDnY7jMyjGAAAAQebeFaVKjA7X1YUjnI7iUyjGAAAAQWRnTZte31ar608fqbjIMKfj+BSKMQAAQBD509uliokI1fWnj3Q6is+hGAMAAASJisYO/XPTfl1TOEJDYiOcjuNzKMYAAABB4v53dissNESL5uQ6HcUnUYwBAACCQHVLl55dX6XLC7KVnhDldByfRDEGAAAIAveuKJW10tfPGOV0FJ9FMQYAAAhw+w506al1lbqsIEc5yTFOx/FZFGMAAIAAd++KUknSrWflO5zEt1GMAQAAAlhVc6eeKarU5QU5ykqKdjqOT6MYAwAABLB7V5TKyOiW+YwW96ffYmyMecgYU2eM2fo5j19tjNl86ON9Y8xU78cEAADAsaps6tQzRVW6YmaOhjFa3K+BjBgvkbTwKI/vkTTPWjtF0n9JesALuQAAAHCC/vhWqUJCjL51JqPFA9FvMbbWrpTUdJTH37fWNh+6uVpStpeyAQAA4DjtbezUsxuqdNXM4cpMZN/igfD2HOObJL3yeQ8aYxYbY4qMMUX19fVefmkAAAB85P/e2qWwEKNvnsm+xQPltWJsjJmvg8X4h593jbX2AWttgbW2IC0tzVsvDQAAgMOUN3ToHx/u01WzhiuDU+4GLMwbT2KMmSLpQUnnW2sbvfGcAAAAOD7/91bpwdHieYwWH4sTHjE2xgyX9A9J11prS048EgAAAI7XnoYOPb9xn64pHKF0RouPSb8jxsaYJyWdKSnVGFMl6S5J4ZJkrb1f0p2SUiT9yRgjSS5rbcFgBQYAAMDn++1rOxURGqKvz8tzOorf6bcYW2uv7OfxRZIWeS0RAAAAjsvWfS16cXO1bp2fr/R4RouPFSffAQAABIhfL9+ppJhwLWa0+LhQjAEAAALA+7sbtLKkXrecma+EqHCn4/glijEAAICfs9bq16/u1NDEKF172gin4/gtijEAAICfe21brTZWHtB3F4xWVHio03H8FsUYAADAj7k9Vr9ZvlN5abH68rRsp+P4NYoxAACAH/vHhiqV1rXrtnPHKiyUanci+OwBAAD4qe4+t37/xi5NyU7UwkmZTsfxexRjAAAAP/X4mr3ad6BLP1w4TocOWsMJoBgDAAD4obbuPt27olRz8lM1Oz/V6TgBgWIMAADghx5YWaamjl79YOFYp6MEDIoxAACAn9l/oEsPrCzTRVOHaUp2ktNxAgbFGAAAwM/8ZvlOWYnRYi+jGAMAAPiRTZUH9NyH+7RoTq6yh8Q4HSegUIwBAAD8hLVWv3hpu1LjIvTNM0c5HSfgUIwBAAD8xPLiGq0tb9L3zxmr+Khwp+MEHIoxAACAH+hxufXfr+zQmIw4XV7A0c+DgWIMAADgBx79oEIVjZ368QUTOPp5kPBZBQAA8HFNHb36w5u7NG9MmuaNSXM6TsCiGAMAAPi4e97cpY4el358wXinowQ0ijEAAIAPK61r16OrK3TlzOEakxHvdJyARjEGAADwUQe3Z9um6PBQfe+cMU7HCXgUYwAAAB/15vY6rdhZr+8uGK3UuEin4wQ8ijEAAIAP6u5z6z9fLFZ+epyuO32k03GCQpjTAQAAAPBZf36nTJVNXXpi0SyFsz3bScFnGQAAwMdUNnXqT2+X6oIpQ3V6fqrTcYIGxRgAAMDH/PylbQoxRj/+AtuznUwUYwAAAB/yTkm9lhfX6taz8jUsKdrpOEGFYgwAAOAjelxu/XRZsXJTY7Vobq7TcYIOxRgAAMBHPPReufY0dOiuL05QZFio03GCDsUYAADAB1S3dOn/3tqlcyZk6Myx6U7HCUoUYwAAAB/w85e2y+2xuvPCCU5HCVoUYwAAAIe9taNWL22u1q3z85WTHON0nKBFMQYAAHBQR49LP3m+WKPT4/T1eaOcjhPUOPkOAADAQXe/XqJ9B7r07DdOU0QYY5ZO4rMPAADgkC1VLXpo1R5dPWu4CkYmOx0n6FGMAQAAHOBye3T7PzYrNS5SP1g4zuk4EFMpAAAAHLHk/XIV72/Vn66epsTocKfjQIwYAwAAnHSVTZ367WslWjA+XedPynQ6Dg6hGAMAAJxE1lrd8fxWhRjpZxdPkjHG6Ug4hGIMAABwEv1zc7XeKanXv587VsOSop2Og8NQjAEAAE6Spo5e/eyfxZqSnajrTh/pdBx8CovvAAAATpKfvLBVLV19emzRLIWGMIXC1zBiDAAAcBK8tLlaL22u1ncXjNG4zASn4+AI+i3GxpiHjDF1xpitn/O4McbcY4wpNcZsNsZM835MAAAA/9XQ3qOfvLBVU7IT9fUz8pyOg88xkBHjJZIWHuXx8yWNPvSxWNJ9Jx4LAAAgMFhr9ZPnt6q926XfXjZVYaF8w95X9fsnY61dKanpKJdcLOkRe9BqSUnGmKHeCggAAODP/rm5Wq9srdH3zx2j0RnxTsfBUXjjLUuWpMrDblcduu8zjDGLjTFFxpii+vp6L7w0AACA76pr69adL2zVKTlJunkuUyh8nTeK8ZGWVNojXWitfcBaW2CtLUhLS/PCSwMAAPgma61+/NxWdfa69b+XTWUXCj/gjWJcJSnnsNvZkvZ74XkBAAD81vMb9+n1bbW67dyxyk+PczoOBsAbxXiZpK8d2p2iUFKLtbbaC88LAADgl6pbuvTTZds0fcQQ3Tgn1+k4GKB+D/gwxjwp6UxJqcaYKkl3SQqXJGvt/ZJelvQFSaWSOiXdMFhhAQAAfJ3bY/W9pzaqz+3Rb5lC4Vf6LcbW2iv7edxKusVriQAAAPzYn1fu1uqyJv3mK1M0MjXW6Tg4BmykBwAA4CWbKg/od6+V6IIpQ/WV6dlOx8ExohgDAAB4QUePS99Z+qHS4yP1y0smyximUPibfqdSAAAAoH8/XVasvU2devLmQiXGhDsdB8eBEWMAAIAT9NLmaj2zvkq3zM/XrLwUp+PgOFGMAQAATsC+A1360T8265ScJH377NFOx8EJoBgDAAAcp4+2ZnN7rP5wxSkKD6Va+TPmGAMAABynP7xRorV7mvTby6ZqRApbs/k73tYAAAAch7d31umet0p1eUG2vszWbAGBYgwAAHCM9h3o0vee2qhxmfH62cWTnI4DL6EYAwAAHINel0e3PL5BfW6r+66ZrqjwUKcjwUuYYwwAAHAMfvnydm2sPKD7rp6mXI58DiiMGAMAAAzQi5v3a8n75bpxdq7OnzzU6TjwMooxAADAAOyub9cPn92sU4cn6fbzxzkdB4OAYgwAANCPrl63vvXYBkWEhejeq6YpIowKFYiYYwwAAHAU1lrd9uwmldS1ackNMzUsKdrpSBgkvN0BAAA4ij+9vVsvbq7WbeeN1bwxaU7HwSCiGAMAAHyO17fV6n9f26mLpg7TN+eNcjoOBhnFGAAA4AhKatv03aUfatKwRP36K1NkjHE6EgYZxRgAAOBTmjt6tehvRYqJDNMDX+MQj2BBMQYAADiMy+3RrU9uUE1Lt/587XQNTWSxXbBgVwoAAIDD/Pyl7VpV2qj/vWyqpg0f4nQcnESMGAMAABzyxJq9WvJ+uRbNydVXpmc7HQcnGcUYAABA0ooddfrJC1s1f2waJ9sFKYoxAAAIepurDuhbj2/QhKEJ+uNV0xQWSkUKRvypAwCAoFbZ1Kkbl6xTSlyE/np9gWIjWYIVrPiTBwAAQau5o1fXPbxWfW6rpYtnKj0+yulIcBDFGAAABKXuPrdufqRIVc1deuymWcpPj3M6EhzGVAoAABB0PB6r7z+9UUUVzbr78lM0MzfZ6UjwARRjAAAQVKy1+tmL2/TylhrdccF4XTBlqNOR4CMoxgAAIKj89rUSLXm/XDfNydVNc3KdjgMfQjEGAABB4763d+uPK0p1xYwc3XHBeBljnI4EH0IxBgAAQeGRD8r1P6/u0EVTh+kXl06mFOMzKMYAACDgPbu+Sne+UKwF4zP028unKjSEUozPohgDAICA9sqWav3g2U2anZ+iP151qsI51Q6fg78ZAAAgYK3YWadvL/1Qpw4for98rUBR4aFOR4IP44APAAAQkN7aUatvPLpBYzLi9dD1MxQTQe3B0TFiDAAAAs7y4hp9/dH1GpsZr8cXzVJidLjTkeAHeOsEBBmPx6qj16XWbpdau/rU3uNSn8ujPo+Vy+1Rn9vK5fHI7bEKMUbhoUZhISEKCzUKDw1ReGiIYiJCFR8VprjIMMVHhSsijPfYAHzHS5ur9Z2lH2pydqKW3DCTUowBoxgDAcLjsapu7VZlU6dqWrpV09qtmpZu1bYe/HlDe49au1xq6+6Tx3r3tSPCQpQQFaakmAilxkUoNS5SafGRB3+Mi1RGYpSyh0QrKyma+X0ABtXzH+7T95/eqOkjhuih62coPopSjIGjGAN+ps/tUUltm3bVtqusvl27GzpUVt+hPQ3t6u7zfOLauMgwZSREKjMxStOHD1FidLgSosOVEBWuhOgwJUSFKzYyTOGhIYoI++TIcGiIkdtj1ef2HPrxXyPKHb0utXe71N5zsGi39bjU1u1Sc0ev6tt6tHVfixrae9Xe4/pM/tS4SGUNiVb2kGgNT45RbmqsRqXFKi81TkNiI07WpxFAAHqmqFI/+PtmFeam6MHrChQbSc3BseFvDODDXG6PdtW1a0tVi7bsa9HmfS3aXt2qXtfBAhxipJzkGOWlxur0USnKS4vViORYZSZGKTMxSnEO/6fQ1etWQ3uPalq7VdXcqX3NXapq7tK+A13atr9VrxXXqM/9r+HrITHhykuLU35anMYNjde4zASNHxqvpBgKM4Cje2LNXv2/57Zo7uhUPXBtgaIj+O4Ujh3FGPAhLrdHW/a16IOyRq0ua1JReZM6e92SpPjIME3MStD1p4/UpKxEjc+M1/CUGEWG+e4X/+iIUOUkxygnOUYzRiZ/5nGX26Oq5i6VNbSrrL5Du+vbtbu+Q69vr9VTRZUfX5eZEKXxQ+M1bmiCpmQlanJ2orKSojm1CoCstfrjW6X67eslmj82TfddM50pWzhuFGPAYbvr2/Xm9lq9v7tR6/Y0qeNQER6TEafLpmdr2oghmpyVqJEpsQoJsJOawkJDNDI1ViNTY3XWuH/db61VfXuPdlS3aUdNq7ZXt2l7dave3dUg16EJ0imxEZqcnagp2UmakpWoU4cnKSUu0qHfCQAnuD1Wdy3bqsdW79WXTs3Sr748hcXAOCEUY+Akc3usNuxt1hvbavX6tlqVNXRIkkalxerSaVk6LS9Vs/KSlRrEJc8Yo/T4KKXHR+mMMWkf39/jcmtHdZs2Vx3Q5qoWba5q0cqSXR8vJsxLjdW0EUNUMGKIpo8YolFpcQH3ZgLAQd19bn136Ua9Wlyjb8wbpR8uHMt3kXDCBlSMjTELJf1BUqikB621v/rU48Ml/U1S0qFrbrfWvuzlrIDfcnusVpU2aNmm/XprR52aOnoVHmpUmJei62eP1NnjM5SVFO10TJ8XGRaqqTlJmpqT9PF9nb0ubd3Xqg17m1VU3qy3dtTp2fVVkqTE6HDNGDlEhXkpmpWbognDEhRKUQb8XktnnxY9sk5FFc2688IJunFOrtORECD6LcbGmFBJ90o6R1KVpHXGmGXW2m2HXXaHpKettfcZYyZIelnSyEHIC/gNa622VbfquQ37tGzTftW19Sg+Kkxnj0vXggkZmjcmjW2EvCAmIkwzc5M1MzdZmnfw817e2Kmi8iatr2jW2j1NemN7nSQpPipMM0YmqzAvWaePStWEoQmMKAN+Zv+BLl3/8FqVN3Tq/648VRdOGeZ0JASQgYwYz5RUaq0tkyRjzFJJF0s6vBhbSQmHfp4oab83QwL+pL6tR3/fUKXnNuzTzto2hYcanTk2XV86NUvzx6WzKGSQGWOUmxqr3NRYXVaQI0mqbe3W6rJGrdnTpNVljXprx8GiPCQmXKePStXs/FTNyU/V8JQYJ6MD6MfWfS26+ZEitXe7tOTGGTp9VKrTkRBgBlKMsyRVHna7StKsT13zU0mvGWP+TVKspAVHeiJjzGJJiyVp+PDhx5oV8FnWWm3Ye0CPflCul7ZUq89tNW14kv7rkkm6cPJQ9ud1WEZClC4+JUsXn5IlSapr7dYHZY16d1eD3tvVoJe2VEuScpKjNSc/TfPGpGl2fgoj+oAPeXHzfv3HM5uUHBOhp75+miYMS+j/FwHHaCDF+EjfZ/z0uVlXSlpirf2tMeY0SY8aYyZZaz9x2oC19gFJD0hSQUGBl8/eAk6+7j63lm3ar0c+KNfWfa2KjwzTNYUjdE3hCI1Ki3M6Hj5H+mFF2VqrsoYOrSpt0Lu7GvTPTfv15Nq9CgsxmjZiiOaNOViUmXYBOMPjsfr9GyW6561SFYwYovuuma60+OBdnIzBNZBiXCUp57Db2frsVImbJC2UJGvtB8aYKEmpkuq8ERLwNU0dvXrovT16bE2FDnT2aWxGvH5+ySRdemoWJy35GWOMRqXFaVRanL522kj1ujzasLdZ75TU652d9frN8p36zfKdSouP1PyxaTprXIbmjE51/PAUIBh09Lj0vac26rVttbq8IFv/dckkn967Hf5vIF/Z10kabYzJlbRP0hWSrvrUNXslnS1piTFmvKQoSfXeDAr4gtrWbv1lZZkeX7NX3S63zp2QoRtm52pWbjLbBAWIiLAQFealqDAvRT9cOE51rd1auatBK3bW6ZWtNXq6qOrjHUXmj03X2ePTNSIl1unYQMCpbOrUzY8UqaS2TXd9cYKuP30kX2cx6Iy1/c9oMMZ8QdLvdXArtoestb8wxvxMUpG1dtmhnSj+IilOB6dZ/MBa+9rRnrOgoMAWFRWd8G8AOBkqmzp1/zu79UxRldzW6uKpw/St+aOUnx7vdDScRH1uj4rKm7ViZ53e3F6r3fUH96AenR6nBRMydM6EDJ2SncSUi2Dy8AUHf7zhJWdzBJiVJfX67lMb5XJ7dO/V0zR3dFr/vwg4CmPMemttQb/XDaQYDwaKMfxBTUu37n69RM9uqFKoMfry9Gx9c94odi+AJKmisUNvbD9YktfsaZLbY5UaF/nxlnxzR6eyC0mgoxh7lcvt0e/f2KV73y7V6PQ43X/NdOWxXgNeMNBizCQ54Ahau/t0/9u79dCqPfJ4pGsLR+gb80YpMzHK6WjwISNSYnXTnFzdNCdXLZ19erukTm9sr9PLW6r1VFGlosNDNW9Mms6dmKGzxqUrKYbdSYDPU9PSrW8/+aHWljfpqwU5+ulFExUdwRtLnFwUY+AwvS6PHl9ToXve3KXmzj5dcsow/fu5Y5WTzAgxji4xJvzjnS56XR6t2dOo14pr9dq2Gr1aXKPQEKPCvGSdOyFT507M0NBETjoEPvJOSb2+99RGdfe5dfdXp+rSU7OdjoQgxVQKQAf3IX55S43+59Ud2tvUqdn5KfrR+eM1KSvR6Wjwcx6P1ZZ9LVpeXKPXttWqtK5dkjQ1J0kLJ2Zq4aRM5aayeM9vMZXihLjcHt39RonuXbFbYzPide/V05SfztQJeB9zjIEBKqtv113LivXurgaNy4zXj74wXmeMTmX1MwbF7vp2LS+u0fKtNdpU1SJJGpsRr/MmZWrhxEyNHxrP3z1/QjE+bqV1bfr3pzdpU1WLrpiRo7u+yNQJDB6KMdCP7j63/rSiVPe/U6bIsBDdtnCsrp41QqHsKICTZN+BLr1WXKNXt9ZoXXmTPFYakRLz8UjyKTlJlGRfRzE+Zm6P1UPv7dFvXtup2IhQ/fySybpgylCnYyHAsfgOOIoVO+p057Ktqmzq0qWnZulHXxin9HgW1uHkykqK1g2zc3XD7Fw1tPfo9W21emVrjf763h79eWWZhiZG6byJmTp/UqYKRibzpg1+r7yhQ//xzCYVVTTrnAkZ+uWlkznFDj6FYoygUt/Wo588v1WvFtdoVFqsnrh5lk4flep0LECpcZG6cuZwXTlzuFo6+/TmjoMl+cm1e7Xk/XKlxkXo3EMluTAvReGhIU5HBgbM47F6fE2FfvnyDoWFGv3u8qm69NQsviMCn0MxRtB4aXO17nh+izp63brtvLG6eW6eIsIoF/A9iTHh+tK0bH1pWrY6elwfn7r3/If79MSavUqKCdeC8Rk6f1Km5oxO5Yhc+LSS2jbd8fxWrd3TpLmjU/Xrr0xhVxb4LIoxAl5zR6/uXFasf27arynZifrd5VM5sQ5+IzYyTBdOGaYLpwxTd59bK0vq9erWGi0vrtGz66sUFxmms8al6/xJmZo3Nk0xEXxZh2/o6HHpnjd36a/v7VFcVJh+9aXJ+uqMHEYVj1JkAAAgAElEQVSJ4dP4CoqA9ub2Wt3+jy060Nmrfz9njL555iiF8S1o+Kmo8FCdOzFT507MVK/Lo/d3N+jVrQe3gVu2ab+iwkM0b0yazp80VGeNT1dCVLjTkRGErLVaXlyrn/2zWPtbunV5QbZuP3+8kmM54Aa+j2KMgNTR49J//rNYTxdVaVxmvJbcMEMTh7EnMQJHRFiIzhybrjPHpuvnl3i0trxJr26tOTSaXKvwUKPTR6Vq4aRMnTMhQ6lxLHDC4Kto7NBPlxVrxc56jcuM1z1XnqqCkclOxwIGjO3aEHB21LTqW49vUHlDh74xb5S+s2A0czARNDweqw8rD2j5oW3g9jZ1KsRIBSOTdd7ETJ07IYOTHL2J7dokSU0dvbrnzV16fE2FIkJD9L1zxui600eySBQ+g32MEXSstXpqXaXuWlashOhw3XPFqTptVIrTsQDHWGu1vbrt4IEixTXaUdMmSZo4LOHjo6nHZXKgyAkJ8mLc1evWQ6v26P63d6uj16WvzsjRdxeMUUYC21/Ct7CPMYJKe49LP35ui17YuF9zR6fqd5efwt6YCHrGGE0YlqAJwxL0vXPGqLyhQ69tq9FrxbX6/ZsluvuNEg1PjtG5EzJ0zoQMTR8xhDn4GBC3x+rv66v0u9dLVNParQXjM/TDhWM1OoOFzfBvFGP4ve3Vrbrl8Q0qb+zQf5w7Rt86M18hHIQAfMbI1FgtPmOUFp8xSnVt3Xpze52WF9fokQ8q9OB7e5QUE66zxqZrwYQMnTEmTXGR/BeBT3K5PVq2ab/uXVGq3fUdOiUnSfdceapm5jKPGIGBr3rwa899WKXb/75FidHheuLmQhXmMXUCGIj0+KiPDxRp6+7Tu7sa9Ma2Wr21s07/+HCfIkJDVDgqRQvGp2v+2HTmJQe5Hpdbz66v0v3v7FZlU5fGZcbrT1dP0/mTMpmKg4BCMYZfcnusfv3qDv15ZZkK85L1x6umseoeOE7xUeH6wuSh+sLkoXK5PVpf0aw3ttfq9W21uvOFYknFGpMRp/nj0nX2uAxNG57ElIsg0dnr0pNrK/XAyt2qbe3R1Jwk3XnhRJ09Lp3vzCEgsfgOfqe1u0/fefJDrdhZr2sLR+jOL05g5TMwSMrq2/XWjjq9taNOa/c0yeWxSowO1xlj0jRvTJrOGJOq9PggXmgVoIvvKps69diaCj29rlLNnX0qzEvWrfNHa3Z+CiPE8EssvkNA2tPQoUV/W6eKxk79/JJJuqZwhNORgICWlxanvLQ4LZqbp9buPr23q0Fvbq/TOyX1+uem/ZKkCUMTNG9sms4ck6ZpI4bwRtVPeTxWq3Y36G/vV+jNHbUKMUbnTsjQTXNy2YsYQYNiDL/x7q563fL4BoWGGD160yy2YgNOsoTDplx4PFbba1r1Tkm93tlZr7+sLNN9b+9WbESoCvNSNDs/VXNGp2p0ehwjjD6usb1HL2zcr8fWVKisvkMpsRG65cx8XTVruIYlRTsdDzipKMbwC4+trtBdy4qVnxanB68rYCEQ4LCQEKOJwxI1cViivnVmvtq6+7SqtFHvldZrVWmj3txRJ0lKi4/UnPxUnT4qRYV5Kfzb9RHdfW69sb1Wz23Yp3dK6uXyWJ2Sk6S7vzpVX5g8lEORELQoxvBp1lr97vUS/d9bpTprXLruufJUtpACfFB8VLgWTsrUwkmZkqSq5k69X9qo90obtLKkXs99uE+SlJUUrVl5ySrM/agoRzOifJK4PVZr9zTp+Q/36eUt1WrrcSkjIVI3zcnVJadmafzQBKcjAo6jYcBnudwe/b/ntujpoip9tSBHv7h0EivhAT+RPSRGl8+I0eUzcuTxWO2qa9eaPY1aXdaod3bW6x8bDhblzIQoTR85RAUjhmj6iCEaPzSBOcpe1Nnr0sqSBr2xvVZv7ahTU0evYiNCtXDSUH1pWpYK81IUyu4SwMcoxvBJnb0u3frEh3prR52+fVa+vnfOGEaVAD8VEmI0NjNeYzPj9bXTRspaq9K6dq3e06S1e5q0oaJZL22uliRFh4dqak6ipo8YoinZSZqanaSMhEj+/R+DyqbOg/tSb6/Ve6UN6nV5lBAVprPGHTy85exxGYqOYKoEcCQUY/icpo5e3bhknTZXHWDnCSAAGWM0OiNeozPide2hf9/VLV1aX9H88cf975TJ7Tm4nWhafKSmZidqclaSpmQnasKwBKXHU5Y/sv9Alz7Y3agPyg6OyFc1d0mSsodE6+pZw3XOhAzNGJnMSDwwABRj+JTKpk5d99Ba7TvQpfuuma7zJmY6HQnASTA0MVoXTonWhVOGSTq4OGxbdas2Vx7Q5n0t2lzVojd31OmjrfeTYyM0LjNe4zITNG5ovMZnJmhUeqxiIgL7v7XOXpe27W/V5qoWbd3XovV7m1XR2ClJSooJV2FuihbNydXp+ewIAhyPwP4KAr+yp6FDV/1ltTp6XHps0SzNYN9MIGhFhYdq2vAhmjZ8yMf3tfe4VLyvRdurW7Wjpk3ba9r0xNoKdfd5Pr5maGKURqXFKS8tVnmpscpLi9OIlBgNTYxWRJj/jJi63B5VNneprL5dZfUd2l7Tqq37WlRa165DA+nKSIjU1OwkXXfaSBXmpWhcZjyn0QEniGIMn1Ba16ar/rJGLo/V0sWnacIwVkcD+KS4yDDNykvRrLx/7WHu9ljtberU9upW7a5rV1lDh8rq2/Xchn1q63F9fJ0xUkZ8lLKHRCtrSLSyh0QrMzFaaXERSo2LPPgRH6nYiNCTMsra6/Korq1bta3dqmnpUXVLl2pbu1XR2Kmyhg5VNHaoz/2vk2nT4iM1OStR508aqinZiZqclaj0hCA+cRAYJBRjOG5nTZuufnC1JKOliws1JiPe6UgA/ERoiFFuaqxyU2M/cb+1VvXtPSqr71BlU6eqmru070CXqpo7tb6iWS9urv54DvPhosNDlRwbofioMMVFhik+KkzxUeGKO3Q7LMQoLDRE4SFGoaFGl7R0yUh6cWWZet0eudxWLo9HfW6rPrdH7d0utXb3qbW7T23dLrV29am126Wmjt7PvHZkWIhykmOUlxqrBeMzlJcWq1FpscpLjdOQ2IjB+hQCOAzFGI4q3t+iax5co4iwED1xc6FGpcU5HQlAADDGKD0+SunxUSrM++wpmS63R00dvapv71FDe6/q23rU0N6jhrYeNXX2qr3bpfYelxrae1Xe2Km27j519Lg/Lr0fmRZxcH7vL17efthrS+EhIQoPNYqLClNCVLgSosOVEhuh3NRYxUeFKTUuUkMTo5SREKXMxChlJkQpMTqcOcGAwyjGcMzmqgO69q9rFRsRqiduLtTIT434AMBgCQsNUXpC1HFNR7DWyu2xcnmswh75oyRp6zXnKSzEKDw0hH2BAT9GMYYjNuxt1nV/XavEmHA9eXMhx8QC8BvGGIWFGoWFSgo5uKCPEzmBwMC/ZJx0W/e16Lq/rlVKXISeuLlQw5KinY4EAAAg/9m7BgFhV22brv3rGiVEh1OKAQCAT6EY46SpaOzQ1Q+uUVhoiB5fNItSDAAAfArFGCdFdUuXrvrLGvW5PXp80SwW2gEAAJ9DMcaga2jv0dUPrlFrV58euXEW+xQDAACfxOI7DKqWzj5d+9e12n+gS4/eNEuTsxOdjgQAAHBEjBhj0HT1unX9krXaXdeuB64t0IyRyU5HAgAA+FyMGGNQuD1W3176oTZWHtB9V0/TGWPSnI4EAABwVIwYw+ustfrPfxbr9W21uuvCCVo4aajTkQAAAPpFMYbXPbCyTI98UKHFZ+Tp+tm5TscBAAAYEIoxvGrZpv3671d26MIpQ3X7wnFOxwEAABgwijG8ZnVZo/7j6U2amZus/71sqkJCjNORAAAABmxAxdgYs9AYs9MYU2qMuf1zrrncGLPNGFNsjHnCuzHh63bVtmnxI0UanhKjB66drqjwUKcjAQAAHJN+d6UwxoRKulfSOZKqJK0zxiyz1m477JrRkn4kaba1ttkYkz5YgeF7Gtp7dP3D6xQZHqolN8xQUkyE05EAAACO2UBGjGdKKrXWlllreyUtlXTxp665WdK91tpmSbLW1nk3JnxVr8ujbzy6Xo0dPXrouhnKHhLjdCQAAIDjMpBinCWp8rDbVYfuO9wYSWOMMauMMauNMQuP9ETGmMXGmCJjTFF9ff3xJYbPsNbqzhe2qqiiWb/5ylROtQMAAH5tIMX4SCuo7Kduh0kaLelMSVdKetAYk/SZX2TtA9baAmttQVoaBz74u0c+qNDSdZW6dX6+vjh1mNNxAAAATshAinGVpJzDbmdL2n+Ea16w1vZZa/dI2qmDRRkBalVpg3724jYtGJ+h758zxuk4AAAAJ2wgxXidpNHGmFxjTISkKyQt+9Q1z0uaL0nGmFQdnFpR5s2g8B0VjR361uMbNCotVnd/lW3ZAABAYOi3GFtrXZJulbRc0nZJT1tri40xPzPGXHTosuWSGo0x2yStkHSbtbZxsELDOW3dfVr0tyIZI/3lawWKjwp3OhIAAIBX9LtdmyRZa1+W9PKn7rvzsJ9bSd8/9IEA5fFYfe+pjSpr6NAjN87UiJRYpyMBAAB4DSffYcDuXVGqN7bX6ScXjNfs/FSn4wAAAHgVxRgDsqq0Qb97o0SXnDJM150+0uk4AAAAXkcxRr9qWrr1naUfKj8tTr+4dLKMYbEdAAAIPAOaY4zg1ef26N+e3KDOXreWLp6m2Ej+ygAAgMBEy8FR/Wb5Tq0rb9YfrjhF+enxTscBAAAYNEylwOdaXlyjB1aW6drCEbr4lE+fAg4AABBYKMY4or2NnfqPZzZpSnai7rhwvNNxAAAABh3FGJ/R3efWNx9frxBjdO9V0xQZFup0JAAAgEHHHGN8xq9e2aHi/a3663UFykmOcToOAADAScGIMT7hrR21WvJ+uW6cnauzx2c4HQcAAOCkoRjjY3Vt3brtmc0alxmvHywc63QcAACAk4qpFJAkeTxWtz2zWe09Li1dXKiocOYVAwCA4MKIMSRJS94v1zsl9brjgvEancF+xQAAIPhQjKHt1a361Ss7tGB8uq4pHOF0HAAAAEdQjINcd59b337yQyXGhOt/vjxFxhinIwEAADiCOcZB7hcvbdeuunY9etNMpcRFOh0HAADAMYwYB7E3t9fq0dUVWjQnV3NHpzkdBwAAwFEU4yB1oLNXt/9ji8Zlxus2tmYDAABgKkWw+umyYjV39Orh62dw5DMAAIAYMQ5Kr26t0fMb9+vWs/I1KSvR6TgAAAA+gWIcZJo6enXH81s0YWiCbpmf73QcAAAAn8FUiiBz5wtb1dLVp0dvmqXwUN4XAQAAfIRmFERe2lytFzdX6ztnj9b4oQlOxwEAAPApFOMg0dDeo5+8sFWTsxL1jXmjnI4DAADgcyjGQcBaq588v1Xt3S7972VTFcYUCgAAgM+gIQWBFzdX65WtNfruOaM1NjPe6TgAAAA+iWIc4A509uqny4o1NTtRi+fmOR0HAADAZ7ErRYD75cvbdaCrT48tmsUUCgAAgKOgKQWwD3Y36umiKt08N49dKAAAAPpBMQ5Q3X1u/fi5LRqeHKPvnD3a6TgAAAA+j6kUAepPK0pV1tChR2+aqeiIUKfjAAAA+DxGjAPQrto23ffObl16apbmjk5zOg4AAIBfoBgHGI/H6kf/2KLYyDDdccF4p+MAAAD4DYpxgHly3V4VVTTrjgsmKCUu0uk4AAAAfoNiHEDqWrv1q1d26PRRKfrytCyn4wAAAPgVinEA+c8Xt6nH5dEvLp0sY4zTcQAAAPwKxThAvLerQS9trtat8/OVmxrrdBwAAAC/QzEOAL0uj+5atlUjUmK0+AyOfQYAADgeFOMA8Lf3y7W7vkN3fXGCosLZsxgAAOB4UIz9XG1rt37/RonOHpeus8ZlOB0HAADAb1GM/dx/v7xdfW6rO784wekoAAAAfo1i7MfWlDXq+Y379fV5eRqRwoI7AACAE0Ex9lMut0d3LStWVlK0vnVmvtNxAAAA/N6AirExZqExZqcxptQYc/tRrvuKMcYaYwq8FxFH8tjqCu2oadNPLhyv6AgW3AEAAJyofouxMSZU0r2Szpc0QdKVxpjPTGg1xsRL+rakNd4OiU+qb+vRb18v0dzRqTpvYqbTcQAAAALCQEaMZ0oqtdaWWWt7JS2VdPERrvsvSb+W1O3FfDiCX7+6Q919bv30oomccAcAAOAlAynGWZIqD7tddei+jxljTpWUY6190YvZcASbqw7omfVVunF2rkalxTkdBwAAIGAMpBgfaUjSfvygMSGS7pb07/0+kTGLjTFFxpii+vr6gaeEJMlaq5+/tF0psRG65SwW3AEAAHjTQIpxlaScw25nS9p/2O14SZMkvW2MKZdUKGnZkRbgWWsfsNYWWGsL0tLSjj91kFpeXKu1e5r0vXPGKCEq3Ok4AAAAAWUgxXidpNHGmFxjTISkKyQt++hBa22LtTbVWjvSWjtS0mpJF1lriwYlcZDqdXn0q1e2a3R6nK6YkdP/LwAAAMAx6bcYW2tdkm6VtFzSdklPW2uLjTE/M8ZcNNgBcdCjqytU3tip/3fBeIWFsv00AACAt4UN5CJr7cuSXv7UfXd+zrVnnngsHO5AZ6/ueXOX5o5O1ZljmIICAAAwGBh69AN/eHOX2rr79OMLxrM9GwAAwCChGPu4svp2PfpBhb46Y7jGZSY4HQcAACBgUYx93H+/skORYSH6/jljnI4CAAAQ0CjGPuyD3Y16fVutvjU/X2nxkU7HAQAACGgUYx/l8Vj94uVtykqK1k1zcp2OAwAAEPAoxj7qn5v3a+u+Vt123lhFhYc6HQcAACDgUYx9UK/Lo9++VqLxQxN00dRhTscBAAAIChRjH/TUur3a29SpHywcq5AQtmcDAAA4GSjGPqaz16U/vFmqmbnJHOYBAABwElGMfczDq8rV0N6jHy4cy2EeAAAAJxHF2Ic0d/Tq/rd3a8H4DE0fkex0HAAAgKBCMfYh97+zW+29Lt123linowAAAAQdirGPqG7p0pL3y3XpqVkamxnvdBwAAICgQzH2EX94Y5eslb63gKOfAQAAnEAx9gGlde16uqhSVxcOV05yjNNxAAAAghLF2Af87vWdig4P1S3z852OAgAAELQoxg7bUtWil7fUaNHcPKXGRTodBwAAIGhRjB129xslSowO16K5uU5HAQAACGoUYwdtrDygt3bU6ea5uYqPCnc6DgAAQFCjGDvo92+UKCkmXNedPtLpKAAAAEGPYuyQDXub9fbOet08N4/RYgAAAB9AMXbI79/YpSGMFgMAAPgMirED1lc0a2VJvRafMUpxkWFOxwEAAIAoxo74/RslSo6N0NdOG+F0FAAAABxCMT7J1lc06d1dDfr6GXmKZbQYAADAZ1CMT7K7X9+l1LgIXctoMQAAgE+hGJ9E68qb9F5pg75+xijFRDBaDAAA4EsoxifR3a+XKDUuUtcUMloMAADgayjGJ8maska9v7tR35iXp+iIUKfjAAAA4FMoxifJH1eUKjUuQlfPYrQYAADAF1GMT4JNlQf07q4G3TSH0WIAAABfRTE+Cf70dqkSosJ0TeFwp6MAAADgc1CMB1lJbZuWF9fq+tNHKj4q3Ok4AAAA+BwU40F239u7FRMRqhtm5zodBQAAAEdBMR5Eexs7tWzTfl01c7iGxEY4HQcAAABHQTEeRPev3K1QY3TzGXlORwEAAEA/KMaDpKalW88WVekrBdnKSIhyOg4AAAD6QTEeJA++Wya3tfrGGaOcjgIAAIABoBgPgqaOXj2+Zq8umjpMw1NinI4DAACAAaAYD4Ilq/aoq8+tb53JaDEAAIC/oBh7WVt3n5a8X67zJmZodEa803EAAAAwQBRjL3t8zV61drt0y/x8p6MAAADgGFCMvajX5dFD7+3R7PwUTclOcjoOAAAAjgHF2Ite2LhPdW09WsxOFAAAAH5nQMXYGLPQGLPTGFNqjLn9CI9/3xizzRiz2RjzpjFmhPej+jZrrf7ybpnGZcbrjNGpTscBAADAMeq3GBtjQiXdK+l8SRMkXWmMmfCpyz6UVGCtnSLpWUm/9nZQX/dOSb1Katt189w8GWOcjgMAAIBjNJAR45mSSq21ZdbaXklLJV18+AXW2hXW2s5DN1dLyvZuTN/3wMoyZSZE6YtThzkdBQAAAMdhIMU4S1LlYberDt33eW6S9MqRHjDGLDbGFBljiurr6wee0sdt3dei93c36obZIxURxrRtAAAAfzSQFnekeQH2iBcac42kAkm/OdLj1toHrLUF1tqCtLS0gaf0cX95t0xxkWG6ctZwp6MAAADgOA2kGFdJyjnsdrak/Z++yBizQNKPJV1kre3xTjzfV9XcqRc3V+uKGTlKiAp3Og4AAACO00CK8TpJo40xucaYCElXSFp2+AXGmFMl/VkHS3Gd92P6rodXlctIunFOrtNRAAAAcAL6LcbWWpekWyUtl7Rd0tPW2mJjzM+MMRcduuw3kuIkPWOM2WiMWfY5TxdQWrr6tHTtXl04ZaiGJUU7HQcAAAAnIGwgF1lrX5b08qfuu/Owny/wci6/8MSaverodevmM/KcjgIAAIATxBYKx6nX5dHDq/ZoTn6qJg5LdDoOAAAAThDF+Dh9dPwzo8UAAACBgWJ8HKy1+ut7ezj+GQAAIIBQjI/D6rIm7ahp0w2zR3L8MwAAQICgGB+Hh1ft0ZCYcF18ytEOAAQAAIA/oRgfo8qmTr2+vVZXzRquqPBQp+MAAADASyjGx+hv75crxBhdUzjC6SgAAADwIorxMejocempokqdPylTQxM50AMAACCQUIyPwd83VKmt26UbZnP8MwAAQKChGA+Qx2O15P1yTc1O1LThSU7HAQAAgJdRjAdo5a56ldV36IbZuWzRBgAAEIAoxgP08KpypcdH6guThzodBQAAAIOAYjwApXXteqekXtcUjlBEGJ8yAACAQETLG4C/vV+uiNAQXTlzuNNRAAAAMEgoxv1o6erT3zdU6YtThyktPtLpOAAAABgkFON+PFNUqc5et26YPdLpKAAAABhEFOOj8His/vZBuWaOTNakrESn4wAAAGAQUYyP4p1d9aps6tK1p3H8MwAAQKCjGB/FYx9UKDUuUudNzHQ6CgAAAAYZxfhzVDZ16q2ddbpyZg5btAEAAAQBGt/neHLtXhmJLdoAAACCBMX4CHpcbj21rlJnj8/QsKRop+MAAID/3969xdhVV3Ec/67OdLhfmhalTOlAQh+stRRaxhp8IGJi66VNtMapEcFgSAxETXwhPmDkSV4k0RANCgHBtBgEHUmRaNDoA9BpauWSQjIhlE5boJbSQqWXqcuH2ch4eto54nTvzdnfT3KSfWav9qzM6n/6O/syRyqBwbiN3z/7CnsOHOaa5d50J0mS1BQG4zbuf3IbA7NP5+OXzKm6FUmSJJXEYNzi+Vf2M/LSXr7y0QFmzIiq25EkSVJJDMYt7n9yG329M1izdF7VrUiSJKlEBuNJ3jo0zsObd/C5xRcw64y+qtuRJElSiQzGkzz8tx0cOHzUT7qTJElqIINxITO5/4ltLOo/m0vnnVN1O5IkSSqZwbiwadteXnj1Ta5ZPkCEN91JkiQ1jcG4cN8T2zjr1F5WXdpfdSuSJEmqgMEY2PPWIR59dhdrls7jtL6eqtuRJElSBQzGwEObd3DkaPLlwflVtyJJkqSKND4YZybrRl5m6cAsFnzwrKrbkSRJUkUaH4xHXtrLi7sP8KUrLqy6FUmSJFWo8cF4/cjLnHVKL59dPLfqViRJklShRgfjfW8fYcMzu1i15AJO7+utuh1JkiRVqNHB+LdbdnDwyL9Y6013kiRJjdfYYJyZrNu4nUX9Z7Oo30+6kyRJarrGBuOnx/axddd+hq7waLEkSZIaHIzXj7zMaTN7WLXkgqpbkSRJUg00MhgfODTO8JadfGbxXM4+dWbV7UiSJKkGGhmMH3l6JwcOH2XtoL+7WJIkSRMaGYzXbdzOgg+cyeXzZ1XdiiRJkmqio2AcESsi4oWIGI2Im9vsPyUiHij2PxURF013o9Pl+Vf2s2X7GwwNziciqm5HkiRJNTFlMI6IHuAOYCWwEFgbEQtbyq4H9mbmJcDtwG3T3eh0Wb9xO309M/j8Zf1VtyJJkqQa6eSI8SAwmpkvZuZhYD2wuqVmNXBvsf0gcHXU8HDswSNHeWjzGJ9adD6zzuiruh1JkiTVSCfBuB/YPun5WPG1tjWZOQ7sA2a3/kURcUNEbIqITbt3735vHf8fXt1/kIvPO5O1V3jTnSRpGpz/kYmHpK7Q20FNuyO/+R5qyMw7gTsBli1bdsz+k21g9hn89sYrySz9pSVJ3WjlD6ruQNI06uSI8Rgw+RDrPGDn8Woiohc4B3h9Oho8GWp4lYckSZIq1kkwHgEWRMTFEdEHDAHDLTXDwLXF9hrg8fSwrCRJkt5HpryUIjPHI+Im4DGgB7g7M5+LiFuBTZk5DNwF3BcRo0wcKR46mU1LkiRJ062Ta4zJzA3Ahpav3TJp+yDwxeltTZIkSSpPIz/5TpIkSWplMJYkSZIwGEuSJEmAwViSJEkCDMaSJEkSYDCWJEmSAIOxJEmSBBiMJUmSJMBgLEmSJAEGY0mSJAkwGEuSJEmAwViSJEkCDMaSJEkSYDCWJEmSAIjMrOaFI3YD2yp5cZgD/KOi11Z7zqSenEv9OJP6cSb15Fzqp8qZDGTmeVMVVRaMqxQRmzJzWdV96F3OpJ6cS/04k/pxJvXkXOrn/TATL6WQJEmSMBhLkiRJQHOD8Z1VN6BjOJN6ci7140zqx5nUk3Opn9rPpJHXGEuSJEmtmnrEWJIkSfovBmNJkiSJLg7GEbEiIl6IiNGIuLnN/lMi4oFi/1MRcVH5XTZPB3O5LiJ2R8SW4vH1Kvpskoi4OyJei4hnj7M/Ilc3z5QAAANvSURBVOJHxcyejojLy+6xaTqYyVURsW/SOrml7B6bJiIujIg/RcTWiHguIr7Vpsa1UrIO5+J6KVFEnBoRGyPi78VMvt+mprYZrCuDcUT0AHcAK4GFwNqIWNhSdj2wNzMvAW4Hbiu3y+bpcC4AD2TmkuLx81KbbKZ7gBUn2L8SWFA8bgB+UkJPTXcPJ54JwF8nrZNbS+ip6caB72Tmh4DlwI1tfn65VsrXyVzA9VKmQ8AnMvNSYAmwIiKWt9TUNoN1ZTAGBoHRzHwxMw8D64HVLTWrgXuL7QeBqyMiSuyxiTqZi0qWmX8BXj9ByWrgFznhSeDciJhbTnfN1MFMVLLM3JWZm4vtN4GtQH9LmWulZB3ORSUq/v2/VTydWTxaf9NDbTNYtwbjfmD7pOdjHLtQ/lOTmePAPmB2Kd01VydzAfhCcRrywYi4sJzWdAKdzk3l+lhxqvLRiPhw1c00SXHa9zLgqZZdrpUKnWAu4HopVUT0RMQW4DXgD5l53LVStwzWrcG43buO1ncrndRoenXyPf8dcFFmLgb+yLvvKFUd10r9bAYGilOVPwZ+U3E/jRERZwK/Br6dmftbd7f5I66VEkwxF9dLyTLzaGYuAeYBgxGxqKWktmulW4PxGDD5SOM8YOfxaiKiFzgHT12ebFPOJTP3ZOah4unPgKUl9abj62Q9qUSZuf+dU5WZuQGYGRFzKm6r60XETCbC1y8z86E2Ja6VCkw1F9dLdTLzDeDPHHvPRG0zWLcG4xFgQURcHBF9wBAw3FIzDFxbbK8BHk8/7eRkm3IuLdfjrWLiejFVaxj4anHH/XJgX2buqrqpJouI89+5Hi8iBpn4Wb6n2q66W/H9vgvYmpk/PE6Za6VknczF9VKuiDgvIs4ttk8DPgk831JW2wzWW3UDJ0NmjkfETcBjQA9wd2Y+FxG3Apsyc5iJhXRfRIwy8S5lqLqOm6HDuXwzIlYxcafx68B1lTXcEBGxDrgKmBMRY8D3mLhZgsz8KbAB+DQwCvwT+Fo1nTZHBzNZA3wjIsaBt4Ghuvyn0sWuBK4BnimunQT4LjAfXCsV6mQurpdyzQXuLX4T1QzgV5n5yPslg/mR0JIkSRLdeymFJEmS9D8xGEuSJEkYjCVJkiTAYCxJkiQBBmNJkiQJMBhLkiRJgMFYkiRJAuDfvB90AJBOlXcAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x = np.linspace( 0, 3, 100 )\n", "f_x = f( x )\n", "\n", "pp.figure( figsize = ( 12, 8 ) )\n", "pp.plot( x, f_x )\n", "pp.plot( [ 2, 2 ], [ np.min( f_x ), np.max( f_x ) ] )\n", "\n", "pp.savefig( 'local_min.pdf' )\n", "pp.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.9" } }, "nbformat": 4, "nbformat_minor": 2 }