Skip to content

Instantly share code, notes, and snippets.

@josePhoenix
Created August 8, 2017 18:42
Show Gist options
  • Save josePhoenix/9a03add597d68bdeb8ecc5878226bea1 to your computer and use it in GitHub Desktop.
Save josePhoenix/9a03add597d68bdeb8ecc5878226bea1 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab inline --no-import-all\n",
"from astropy import units as u\n",
"import poppy"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEGCAYAAACXVXXgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHfZJREFUeJzt3X+0HHWZ5/H3B1BA+U0YQYhARIKoiMioODvSSo6TMEMi\nRNcwBiXCBDHIZiRzlqyCK3EYdshRcUABx4jACAibVTLuRIlIGCfiJLgXCGo0RjRBIgZM+BGEAM/+\nUd97qVzuj+6+3V1d1Z/XOX1uV3VV9XO/Xf30009XVysiMDOzctuh6ADMzGzsnMzNzCrAydzMrAKc\nzM3MKsDJ3MysApzMzcwqwMm85CTVJN1fdByNUuZaSZslrSg6njKRdKakO9L1HSU9IelVBYdlBXMy\nT9ITov/yvKSnctMfKDq+4UTEHRHxuv5pSRsk1XLTh0nqxi8T1IDjgVdGxNvHujFJO0kKSU/mHrcr\nxxxll4uI5yJit4j4TRH3L+lUST+UtFXSsjqWnynp1+nxWSxpr07E2QuczJP0hNgtInYDfgOclJv3\nL4OXl7RT56OslIOBX0XE1kZXHGXsX5d73D7SfHidUYH96BHgs8Cloy0o6Sjgi8AHgP2BbcDlbY2u\nhziZ10nSZyTdJOkGSY8DMyUdJ+mu1Cp4SNIXJL0kLd9fKZ4laa2kP0j6Qm57h0u6U9IWSZskfX3Q\nemdL+qWkxyV9StJr0n09lmLov59Jkh5I128AXgn8W6p8Pg7cmW7rr1b/NG3r+5IeTfd9naQ9c7Ft\nkPRxSfel+G6QtHPu9lMk9aVY1kp6d5q/l6SvprHYIOkiSS/axyTNBq4E/jzFdEGa/5G0vUckfVPS\nAYPG5KOS1gI/a9FjOl/SRkkPSvqbdB+H5P6X6yX9XtIDaVml286UtFzS59Jjv65/DOq4z8PS/cyS\n9Bvgu5J2kHRLimWzpDskvTa3zn6S/jWN913AobnbdhoU9w8knZ67Pd+S2SHtow+nx/VeSUeOZQwj\n4rsRcTPwUB2LzwS+GRE/iIgngAuB90l62eAFG30eGBARvgy6AA8AkwbN+wzwDHAS2YvgrsCfAm8F\ndgImAD8HzknL7wQE8C1gT+AQ4NH+7QI3A/89bWsX4M8GrbcY2B04Kt3vbWkbe5Mlsw+k5ScBD+Ti\n3ADUctOHZQ/zdv/L4cAJwEuBPwH+A1g4aBt3kVVP+6b/68x029uBzWn9HYDxwMR02xKyyutlwCuA\nu4EzhhnjM4E7ctPvBh4Gjk7j8UXg9kFjsjT9/7sOsb3+ZX4LbARuAQ4e4TH+q7Tsa4GXAzek9Q9J\nt3899xhMANYCH8rFvg34MLAj8DFgfZ371mHpfr6axmnXNI6np/vahaxaXZVb55YU38vS/vBQ/9jl\n/u/+uH8AnD7UOAN/Cfwn2f64A3AksH+67RPpcR3qsqmO/+sjwLJRlvk2cN6geU8Bbxzh8azreeBL\nOJkPOSjDJ/PbR1lvHnBzut6/M74td/tiYF66/nXgS8CBg7bRv95bc/PuyT8JgMtIyZcmkvkQcb8X\nWDloGzNy058FLk/XvwJcOsQ2DkxPzJ1z804DbhvmPgcn868BF+em9wCeAw7Kjck7RvgfBPw52QvU\n3mls7wF2HGb5a4EFuekj0n0cArwEeBY4PHf7nP5klWL/2aBYAxhXx77Vn8xfNcIy49IyL8/Fclju\n9n+kuWT+brIE+FZghxY/Z+pJ5stJRUFu3u+A/zLEsg09D3wJt1katD4/IekISd9Ob48fAy4ieyLm\nbcxd3wrslq6fR/ZEXZXaGR8atN7vctefGmJ6N5okaX9J30jthceAaxqIezzwyyE2ezCwM/C71CrY\nDFxBVqHX45XAr/snIuIx4A9kLxL91g9eKbd8RMS/R8QzEfEH4FyydyCHS5qQazNtzt1ffnv5639C\nVnH/Ojfv14NiGTw+0NhjMnB/yo5I+cfUrnmM7F0AZI/JK1Is+fjycdUtIr5L1t76EtnjdKWk3ZvZ\nVpOeIHvhy9sDeHyEddr2PKgaJ/PGDD4q5CpgNVnVtAdZD1B1bSjioYg4MyIOIKv6rpZ06GjrNRHj\nUEey/C/gaeANKe7TqTNusqTy6mHmbwX2iYi90mWPiDiqzu3+luwFAYCUZPYGHswt08hROZEuioh1\n8cKHov1HTzxEVvX3G5+7/jDZu4KDc/NeNSiWMYlUWiYfBE4E3kXWAjkszRdZ8np+UHwjHYb4JFk7\npt/+g+738xFxDPB6sjbLxwEkXaDtj+h6YogXwLG6H3hj/4Skw8ly0C9atP2e5mQ+NrsDW4An0wdW\nZ9W7oqT/Kqm/0ttMlniea0FMvyPr8fZ7GAhJ+Xm7kz3pt0gaT9YeqtdXgDMlvTN9oHaQpIkRsZ7s\nbfRCSXuk2w6T9I46t3sDcIako5R92PoPwL9HxIZ6Vpb0BklvTFXu7sDnyCrYnw+zyjfS/U1MH8Bd\n0H9DRGwj61NfLGm39CL7t8D1dcbyGdVxmF7O7mQvro+QJeK/HxTLN4FPS9pV0uvJ2lfD6QOmp2UP\nJ+vr98f1lnTZiezxf4bshYKIWJB7wRt8GfbwwTTeu5C1RXaQtIuGP0LneuA9kt4u6eVk72RvjiaO\naLIXczIfm/OAD5G9TbwKuKmBdd8KrJT0JFkvfU605ljhi8me+JslzY2Ix8kS44/SvGOBTwFvIXsh\nuhX43/VuPCJWAH8DfCGt/31eqBpnkvV5f0LWIrmZQZXhCNtdSvbk/j9kVfOryA5hq9cryBL0Y2Rt\noIOAv4qIZ4e5vyVk7YY7ySrD/0g3PZ3+fpQs2T1A9iL1NbI+ez3G57ZXj6+SvTP5LVn1OvhLVGeT\nvUv5HdmL6VdH2NZCssLgYWAR278A7ZXW30z2fz1E9nnIWMwia3f8E/DOdP1K2O4LTccBRMS9wDnA\njSm+nck+PLYW0Pbv9sx6k6Q3AD8m+wD3+TFu617g+NS7N+sIJ3PrWZJOJjtcbjfgOuCpiHhvsVGZ\nNcdtli4gaVH6IsfqomPpMXOATWRHj/wxTZuVkivzLpA+JHwCuDYiXl90PGZWPq7Mu0BE3En27VAz\ns6aU/SQ/PUPZ+Uxmp8k3FxmL9YRNEbFf0UFY/ZzMSyIirgauBlB3ntLWqqWpb5lacdxmMTOrAFfm\nZk0YfOCAVO/ZEMzaw5V5F1B2HvIfAhOVnQf8jKJjsqHlzthn1lVcmXeBiDi16BhsePUk74hwdW6F\ncjI3G6TZytsJ3YrkZG5G8wl8qO04oVsRnMytp7n/bVXhZG49pRPJ29W5FcHJ3CqviOrbCd06zcnc\nKqs/iV9yySUAnH/++UWGY9ZWTuZWGcNV4H19fR2OJOPq3DrJp8AtIZ+b5QVl2H9LmtDvjohjiw7C\n6ufK3EqnDAncrNOczK0UxpLAFyxYAMAFF1zQqnAa4naLdYKTuXWtVlXga9asacl2xsIJ3drNPfMS\nqnLPvOr7Y4kSunvmJePK3ApV9eQ9mCt0axcncytEJ5P4hRdeCMBFF13Usfs06zQnc+uYoqrw9evX\nF3K/w3F1bu3gnnkJlaVn7n1rZF2e0N0zLxlX5tZSTuBmxfDPxllLdPPPqc2fP5/58+cXHcaLdOt4\nWTm5MremlCkRPfLII0WHMCz3z61V3DMvoaJ65t5X2qcLE7p75iXjytxG5ARuVg5O5vYiVUvg8+bN\nA2DhwoUFRzI8t1tsrJzMDaheAs976qmnig6hLk7oNhbumZdQK3rmfty7V5ckdPfMS8aVeY9xEjer\nJh9nbpU3d+5c5s6dW3QYdfMLrjXDydysCzmhW6PcMy+hsfTM/XiXS4H9c/fMS8aVuZlZBTiZW+XN\nmTOHOXPmFB1GU/xOyurlo1ms8nbdddeiQxgTH39u9XDPvITcM+89BSRz98xLxm0WsxLwi7CNxsnc\nKm/27NnMnj276DDGzAndRuKeuVXevvvuW3QILeP+uQ3HPfMScs+8t3UombtnXjJus5iVjF+QbShO\n5lZ5s2bNYtasWUWH0VJO6DaYe+ZWeePHjy86hLZw/9zy3DMvIffMrV8bk7l75iXjNotZifnF2fo5\nmVvlzZw5k5kzZxYdRts4oRu4Z249YOLEiUWHYNZ27pmXkHvmNpQW98/dMy8Zt1nMKsIv1L3Nydwq\nb8aMGcyYMaPoMDrCCb13uWdulXf00UcXHYJZ27lnXkLumdtoWtA/d8+8ZNxmMasgv2j3Hidzq7zp\n06czffr0osPoOCf03uKeuVXecccdV3QIZm3nnnkJuWdujWiyf+6eecm4Mm+ApGPqWGxbRNzX9mDM\n6uSzK/YGJ/PGLAdWAiM9Mw4FDulINFaXqVOnAnDrrbcWHIlZ+ziZN2ZlRLxrpAUk3d6pYKw+J5xw\nQtEhFM7VefW5Z15C7plbsxpI6O6Zl4wPTWyCpD+T9PJ0faakz0o6uOi4zKx3OZk350vAVklvBM4D\nfglcW2xINpwpU6YwZcqUosPoCn5nVl3umTfn2YgISdOAyyPiK5LOKDooG9pJJ51UdAhdxf3zanLP\nvAmSlgNLgVnAO4CHgXsi4g0dun/3zG3MRkno7pmXjNsszXk/8DRwRkRsBA4CLi02JDPrZa7MGyRp\nR2BZRLyzwBhcmTdg0qRJACxbtqzgSLrPCNW5K/OScc+8QRHxnKTnJe0ZEVuKjsdG9/73v7/oELqW\n++fV4cq8CZK+BbwJuA14sn9+RJzboft3ZW4tNURCd2VeMq7Mm7M4XczMuoKTeRMi4muSdgVeFRFr\nio7HRlar1QC44447Co2jm7ndUn5O5k2QdBKwEHgpcKiko4GLImJqsZHZUE4//fSiQygFJ/Ryc8+8\nCZLuBt4F3BERb0rzVkfE68ewzcnAZcCOwD9HxCUjLOueubVNSujumZeMjzNvzrYhjmR5vtmNpcMd\nrwCmAEcCp0o6cgzxWc62bdvYtm1b0WGYtZWTeXPul/TXwI6SXiPpn4AVY9jeW4C1EbEuIp4BbgSm\ntSLQwTZt2kStVmPJkiUAbNy4kVqtxtKlSwFYv349tVpt4JjsdevWUavVWL58OQBr1qyhVquxYkX2\n765evZparcbKlSsB6Ovro1ar0dfXB8DKlSup1WqsXr0agBUrVlCr1VizJvuoYfny5dRqNdatWwdk\nx4LXajXWr18PwNKlS6nVamzcuBGAJUuWUKvV2LRpEwCLFy+mVquxZUv22nrTTTdRq9XYunUrANdf\nfz377bffwLHm11xzzUAPHeDLX/7ywG0AX/ziF7c7j8tll102cD50gIULF273e6KXXHIJM2bMGJhe\nsGABM2fOHJi+8MILmTVr1sD0/PnzmT179sD0vHnzmDNnzsD03LlzmTt37sD0nDlzmDdv3sD07Nmz\nmT9//sD0rFmzuPDCCwemZ86cyYIFCwamZ8yYwSWXvPAmb/r06SxcuHBgeurUqVx22WUD01OmTOGK\nK67AysfJvDkfA15H9i3QrwNbgP82hu0dCKzPTW9I8wZImi1plaRVY7ifnnTAAQdsl1BtZPkXFysP\n98ybIOl9EXHzaPMa2N57gckRcWaaPg14a0ScM8zy7plb27hnXk6uzJszv8559XoQGJ+bPijNsxbY\nunXrQNvFRuajWcrLhyY2QNIU4ETgQElfyN20B/DsGDa9EniNpEPJkvgM4K/HsD3LOfHEEwEfZz4a\nJ/JyczJvzG+BVcBU4O7c/MeBv212oxHxrKRzgO+QHZq4KCLuH0ug9oKzzz676BC6nhN5+bln3gRJ\nLyF7ISzkG6DumVur+dws5eeeeXMmA31kP1CBpKMl3VpsSDacLVu2DBy6aC/mqrwanMyb8z/Jjg3f\nDBARfcChRQZkw5s2bRrTprXlsP3ScyKvDvfMm7MtIrYMeiK4f9Glzj23I2cmLh0n8mpxMm/Odt8A\nBc5lbN8AtTY65ZRTig7BrO3cZmlO/hugNwCPAXNHXMMKs2nTpoGv/1vGVXn1+GiWEvLRLI3x+cy3\nV2ci99EsJeM2SxMkHQv8D+AQcmMYEUcVFZMN77zzzis6hK7hiry6XJk3QdIa4O+A+8id+jYift2h\n+3dlbk1pIJm7Mi8ZV+bN+X1E+Ljykug/fe7+++9fcCTFclVeba7MmyDpBOBU4HtkH4ICEBEd+ZFn\nV+aNcc+8qUTuyrxkXJk3ZxZwBPASXmizBNCRZG6NOf/884sOwaztXJk3QdKaiJhY4P27Mre6Ndle\ncWVeMj7OvDkr/Bud5bF+/fqBn6HrNe6T9w63WZrzNqBP0q/IeuYCwocmdqfTTjsN6O2euVWfk3lz\nJhcdgNXvk5/8ZNEhFMJVeW9xz7wBkn4cEceMdZkWxOGeuY2oBYncPfOScWXemNdKuneE2wXs2alg\nrD7r1q0DYMKECQVH0hmuyHuTk3ljjqhjmefaHoU15MMf/jDgnrlVm5N5Azr1dX1rrU9/+tNFh9Ax\nrsp7l3vmJeSeuQ2lxYncPfOS8XHmVnlr1qxhzZqO/+52R7kiN7dZrPLOOusswD1zqzYnc6u8iy++\nuOgQ2spVuYF75qXknrn1a2Mid8+8ZNwzt8pbvXo1q1evLjoMs7ZyZV5CrswbU9Xzmbe5veLKvGTc\nM7fKu/TSS4sOoeXcJ7fBXJmXkCvz3tahRO7KvGTcM7fK6+vro6+vr+gwzNrKlXkJuTJvTJV65h1s\nr7gyLxn3zK3yPv/5zxcdQku4T24jcWVeQq7Me08BidyVecm4Z26Vt3LlSlauXFl0GGZt5cq8hFyZ\nN6bsPfOC2iuuzEvGPXOrvMsvv7zoEJrmPrnVy5V5Cbky7x0FJnNX5iXjnrlV3ooVK1ixYkXRYTTM\nVbk1wpV5Cbkyb0wZe+ZdkMhdmZeMe+ZWeVdddVXRITSkCxK5lZAr8xIaS2Xez4979+qSZO7KvGRc\nmfeowQmjysl9+fLlABx//PEFRzK6LknkVkKuzEuoFZX5SKq2T5SlZ95lidyVecm4MrcXySeVKiT2\nRYsWFR3CqLoskVsJOZnbiKqQ2CdMmFB0CGZt52RudStrYl+2bBkAkyZNKjiSobkqt1Zwz7yE2t0z\nb0Y370fd3DPv4kTunnnJuDK3luhPSt2Y1K+77rqiQxhSFydyKyEnc2upbmzFjB8/vugQzNrOydza\npluOZV+6dCkAkydPLuT+h+Kq3FrNPfMS6saeeaM6ud91W8+8JIncPfOScTIvoSok836d2P82btwI\nwP7779/2+6qHk7m1g9ssVqhOtGK6JYlDaRK5lZCTuXWVdnyAumTJEgBOOumklmyvWU7k1k5us5RQ\nldos9RrLftoNPfMSJnK3WUrGybyEejGZ5zW6z27atAmAcePGtSOcUZUwkYOTeem4zWKl02grpqgk\nbtZJTuZWavUk9sWLFwNwyimndCSmvJJW5VZCbrOUUK+3WerVv28X1TMveSJ3m6VkXJlbZfUn082b\nNxcciVn7OZlb5e21114D1zv1TrTkVbmVkJO59ZROfEnJidyKsEPRAZgVSZKTr1WCK3MzWvfNU78w\nWFGczM0GabYV40RuRXIyNxtFPVW7E7kVzT1zswa4x27dypW5WROc0K3buDIvmKT3Sbpf0vOS/I07\nM2uKk3nxVgOnAHcWHYiZlZfbLAWLiJ+C37ab2dg4mZeEpNnA7DT5NFlFb/UbB2wqOogSmVh0ANYY\nJ/MOkLQMGOqHKD8REd+qZxsRcTVwddreKp/RrjEes8ZIWlV0DNYYJ/MOiIhJRcdgZtXmD0DNzCrA\nybxgkk6WtAE4Dvi2pO/UsdrVbQ6rijxmjfF4lYx/acjMrAJcmZuZVYCTuZlZBTiZl5RPA1AfSZMl\nrZG0VtL5RcfT7SQtkvSwJH+PoWSczMvLpwEYhaQdgSuAKcCRwKmSjiw2qq53DTC56CCscU7mJRUR\nP42INUXH0eXeAqyNiHUR8QxwIzCt4Ji6WkTcCTxadBzWOCdzq7IDgfW56Q1pnlnl+BugXawVpwEw\ns97gZN7FfBqAMXsQGJ+bPijNM6sct1msylYCr5F0qKSXAjOAWwuOyawtnMxLqsnTAPSUiHgWOAf4\nDvBT4BsRcX+xUXU3STcAPwQmStog6YyiY7L6+Ov8ZmYV4MrczKwCnMzNzCrAydzMrAKczM3MKsDJ\n3MysApzMzcwqwMnczKwCnMzNzCrAydzMrAKczM3MKsDJ3MysApzMzcwqwMnczKwCnMzNzCrAydzM\nrAKczM3MKqBlvwEqachfuZA00joNze+Gdbohhlau061xlXWdVm6vSv9L1dZp5X3cfffd34mIycNu\nsE4t/UFnSQMB918farqeZbwNb8Pb8Da6eRst3OY4WsBtFjOzCnAyNzOrACdzM7MKcDI3M6sAJ3Mz\nswpwMjczqwAnczOzCnAyNzOrACdzM7MKcDI3M6sAJ3MzswpwMjczqwAnczOzCnAyNzOrACdzM7MK\naOX5zO+PiD9GDPkbFda4ccCmooOoEI9n63gsW2uXVmyklcn8jxFxbAu319MkrfJ4to7Hs3U8lq0l\naVUrtuM2i5lZBTiZm5lVQCuT+dUt3JZ5PFvN49k6HsvWasl4yh9YmpmVn9ssZmYV4GRuZlYBdSVz\nSZMlrZG0VtL5Q9y+s6Sb0u0/knRI7rb5af4aSX/RutDLqY6x/Likn0i6V9L3JB2cu+05SX3pcmtn\nI+9OdYzn6ZJ+nxu3M3O3fUjSL9LlQ52NvDvVMZ6fy43lzyVtzt3m/TNH0iJJD0taPcztkvSFNNb3\nSjomd1vj+2ZEjHgBdgR+CUwAXgrcAxw5aJmPAlem6zOAm9L1I9PyOwOHpu3sONp9VvVS51i+E3hZ\nun52/1im6SeK/h+66VLneJ4OXD7EuvsA69LfvdP1vYv+n7p9PAct/zFgUW7a++f24/MO4Bhg9TC3\nnwj8GyDgbcCP0vym9s16KvO3AGsjYl1EPAPcCEwbtMw04Gvp+i3ACZKU5t8YEU9HxK+AtWl7vWrU\nsYyI70fE1jR5F3BQh2Msk3r2zeH8BXBbRDwaEX8AbgMmtynOsmh0PE8FbuhIZCUUEXcCj46wyDTg\n2sjcBewl6QCa3DfrSeYHAutz0xvSvCGXiYhngS3AvnWu20saHY8zyF65++0iaZWkuyS9px0Blky9\n4zk9vY29RdL4BtftJXWPSWr/HQrcnpvt/bMxw413U/tmK7/Oby0kaSZwLHB8bvbBEfGgpAnA7ZLu\ni4hfFhNhaSwBboiIpyWdRfYO8l0Fx1QFM4BbIuK53DzvnwWqpzJ/EBifmz4ozRtyGUk7AXsCj9S5\nbi+pazwkTQI+AUyNiKf750fEg+nvOuAO4E3tDLYERh3PiHgkN4b/DLy53nV7UCNjMoNBLRbvnw0b\nbryb2zfraOLvRNaAP5QXPhR53aBl5rD9B6DfSNdfx/YfgK6jtz8ArWcs30T2IdRrBs3fG9g5XR8H\n/IIRPpzqhUud43lA7vrJwF3p+j7Ar9K47p2u71P0/9Tt45mWOwJ4gPSlwzTP++fQY3oIw38A+pds\n/wHof6b5Te2bo7ZZIuJZSecA3yH7tHtRRNwv6SJgVUTcCnwFuE7SWrKG/4y07v2SvgH8BHgWmBPb\nvy3rKXWO5aXAbsDN2WfI/CYipgKvBa6S9DzZO6pLIuInhfwjXaLO8TxX0lSy/e9RsqNbiIhHJS0A\nVqbNXRQRI31YVXl1jidkz+8bI2WexPvnIJJuAGrAOEkbgE8BLwGIiCuB/0t2RMtaYCswK93W1L7p\nr/ObmVWAvwFqZlYBTuZmZhXgZG5mVgFO5mZmFeBkbmZWAU7m1pUGnYGvL38mzjFs8z2SjsxNX5S+\noGVWej400bqSpCciYrcRbt8psvMANbLNa4B/jYhbxhqfWbdxZW6lkc5Nfquk24HvSdotnfP9x5Lu\nkzQtt+wH08m17pF0naS3A1OBS1Ol/2pJ10h6b1r+BEn/L21nkaSd0/wHJH06dx9HFPLPm43CJ9qy\nbrWrpL50/VcRcXK6fgxwVPqW3E7AyRHxmKRxwF3pRxGOBD4JvD0iNknaJy1/K7nKPH3DFkm7ANcA\nJ0TEzyVdS3Yu+c+n+9wUEcdI+igwDxj4gQuzbuHK3LrVUxFxdLqcnJt/W+6rzQIulnQvsIzsNKGv\nIDsr4s0RsQmyr0ePcl8TyV4wfp6mv0b2wwL9Fqe/d5Oda8Os67gyt7J5Mnf9A8B+wJsjYpukB4Bd\n2nCf/WddfA4/Z6xLuTK3MtsTeDgl8ncC/b+XejvwPkn7AkjaJ81/HNh9iO2sAQ6RdFiaPg1Y3r6w\nzVrPydzK7F+AYyXdB3wQ+BlkZ+sE/h5YLuke4LNp+RuBv0sfdL66fyMR8UeyM9bdnLb1PHBl5/4N\ns7HzoYlmZhXgytzMrAKczM3MKsDJ3MysApzMzcwqwMnczKwCnMzNzCrAydzMrAL+P+KILUTtiwWK\nAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x113a7c2e8>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ8AAAEKCAYAAAAM4tCNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmwZVd13r/1xp7Uk1pIrQFNCDCyKTEpyDiKwmBjQlDE\nKMc2kEC1iY1dSTlVEUWCK6QqxqlyxSGGQBdRWSQ2iGADbRCWDYQoiSNoyZaQhBA0EkjdCAlarZZ6\nev2GlT/O/tjrrnPue/edO5177/pVvdpn2Pecvc85b++11157LVFVBEEQrJepYRcgCILRJBqPIAhq\nEY1HEAS1iMYjCIJaROMRBEEtovEIgqAWPWk8RORGEXlcRO5tc/4aETkqInelv/f14r5BEAyPmR5d\n548A/CGAj6+S53+r6mt7dL8gCIZMTyQPVb0NwBO9uFYQBKNBrySPTrhKRO4G8AMA/1JV76vKJCJ7\nAOxJuy8aVOFGFTHbYStcEM9kfaiqrJ2rjPTKPF1ELgLweVX96YpzWwGsqOoxEXkNgP+kqpd1cM14\n921gqz8LYCVtLwypLE1jPqVTABbT9tKQyjIK1G08BjLboqpPqeqxtH0LgFkR2TWIe48rityrTiGm\nzSz2edjnFPSWgXxzInKOiEjavjLd9/Ag7h0EQX/oic5DRD4B4BoAu0TkIIDfQSFRQ1U/AuCNAP6Z\niCwBOAngeo3lvF2xbFLKnNPu3KTB+nMYp5jcZzEIeqbz6Aeh81ibGaRW2nA6pZPyj8NGY84dX0To\nOjqh0TqPIAjGj0FO1QZ9YAplcZ0vdVIkD9aX3Sd7xEmp/7AIySMIglqE5DEGsIddcfuTgq9vKMoG\nQ0geQRDUIiSPEec0ss6jlsp8jKAEsuT2g/4QkkcQBLUIyWOM2JBSbyx2bAhlGQRbUurre3IIZZlE\nQvIIgqAWIXmMAbSs3JhSb+fBWZgTAytRf9nkUkoeXtcREkh/icZjDJh2qV1RavfHhak2qX8OQX8Z\nt+8qCIIBEZLHGMCFcKdSyoVyyy4dF1gf1puSxqI7HvSXkDyCIKhFLMkfI+h+j5IHjcb4EJeRpZNR\ne7CsywaUjeJYF0oe4Y5xfcSS/CAIBkroPMYQP+tge2geGzXDsc0pnUO1RAVkySMYDCF5BEFQi5A8\nxog5l/qXa3Ueo4Z1dOTtOJZcntB5DIZoPCYA61lr1BSlxJab9VmpyhgMjBi2BEFQi5A8xoin3T6V\njFQwziKvfxm1dR8s9yzKClMOxXz9g/4SkkcQBLUIyWMMoXk2/V1QgWo9rbP3fjKlTTNhZzm3p9RO\n1XpfrWGOPhxC8giCoBYheYwhs21SQe61vQ+QphmNsVx+2nkGWdfBKVrWL6ZoB0tIHkEQ1KInjYeI\n3Cgij4vIvW3Oi4h8UEQOiMg3ROSFvbhvUM00Wg2pFtPfafO3jObpOapgOW3ZWR/i6xsMhl5JHn8E\n4NWrnP9FAJelvz0A/kuP7hsEwZDoic5DVW8TkYtWyXItgI9rsf7/dhHZLiK7VfXRbu/tPWcHZR+e\ntlemzkPduaY8R18elpMzKjbq/bg6O+oFg3ifg1KYngfgEbN/MB2r3XhwqtEHeR5V8+tecjyl1gcG\nUCgf59057pOjfSxXJ3B6+YyU8gO1SlLvOY31DQb7f9G42RYR2YNiaBMEQYMZVONxCMAFZv/8dKyE\nqu4FsBdY3ZOY9x4VEkfGPwtKF1vQajAG5KFNUwytWFYahXHKlj3oaeRp5VH1itZPBvl/Maip2n0A\n3ppmXV4K4Ggv9B1BEAyPnkgeIvIJANcA2CUiBwH8DpLtjqp+BMAtAF4D4ACK2EP/pBf3Bdq3rIJs\nPESa0rsOCo57rZ8P6j9WXJ6m9N5ekcsP1PZycy7PpDHn9hdRfn+DeJ+9mm35pTXOK4Df6MW9giBo\nBo1TmPaKaZQ19dxfKmcfS6jPoNn2BpTHwpzKs1LKsCS0OZSnGPmuxOyzPqPmVqBb+P36xY0rGM43\nHebpQRDUYmwlD0XZQIYtNnUhtsUeR0Mj9tC03bAv28d3pV3FPHJA7CdS2q9nw3vvTOkmlPUZ3j/p\nceT6jPNCOKvz8W4XvfvFYemrxrbxsGs3vOLQin3j2Gh4WMcqP6D8APhs5pEbVw5x6POj17DB2pbS\nDSh/kP4fpqpTGEessZcfHnBYOeznEMOWIAhqMbaSh8WvgRC3D7Sf/hsHpRx79qphAbE9OyWOfted\n1+f95pHfjR9WcX8Tcn3G4d1sdPv8Jr23NLs9bImDhOQRBEEtJkLyID4c4Qpy68kegObR3lipaZ62\nOsHqMYBWacP3XlRMnkZWmPZbIcnr8342lGS7D3MauT6s3yga/3EBIL871ttPQ9tYO00LpxmSRxAE\ntZgoycP7hgByL8YHYadx7fF583v2ErxOU0y7PSzfCXPMh6RkSrP1KfO7fsd44fU3mX2Wx3tGt+kJ\nd6yp8Dthnez3Y32yAq1R/exvm2zQGJJHEAS1mCjJowpvA+G13Tw+i3IvQR2C7dmbyI9SuhHAmWmb\nY+6tKWXvuIQshbB+9J3QKx0Ipb3dKd2V0s3Iz5hSxVMppWHYYYzOLAslK78YcQnl701d2pQZldUI\nySMIglpMvOTBMSVd2bXTaM+hvBSaPQl75Kb3FieRX7jXN1jJY8Wdo5SymuTh7UZWexa8ntV1cN/r\nAhbd/ihIHVWuECz0AA+U61Plq7WpTHzjQRZc6v/Jqtz7Uxy1Slb+0zRx5ecMqqdtgdapWr/ylvX0\nAaYtnTSc3qeqn4KsUoCynFax3aR/LD/Fb83JvTNnmH3feDSpTp0Sw5YgCGoRkkcbqnoCivNscX0o\nRCD3nn7qrSnTit4wi72/VeaxV6REtcPtP+6usRYcnjwjpRy2+KHIEspKae+TpCn4cJ1++AeUV2z7\n1bA2zygSkkcQBLUIyWMN2DOcQh67el2HNQLynq82uXQBw9ODLCFLC6yDd1swjey5fMbl8bqK73d4\n37NSyilZq0diuYBCOvPKZ3pIP+HyDhqv2yB+Id8sWj18AeU4M8sYbYmDhOQRBEEtQvJYgyopgUvC\n2QvZpdLUF7BXZY/FB71g8gyj96G0RCmCEtG8Oe5N9PkM6i6YY35KLF5fwB5sEbl39vfwnvAHyQyA\n7Wmbz8lKS0DrlLL3SM/nN+xofL0mJI8gCGoRkkcN2JNscseXUNYP+JmFBeRei6bXsy5P3YVfnQQ3\n9gZafmn4HHKPwvs/nVKauf9wneVifkocvKfXGc2bbS/xbcHa1A3uzHvymfD9UJrYimzG7yUPH71u\nAWUT+ybZ+vSSaDxq4Fer2tW1XiT3/hpWkBWS/Fi9SM5G5XH0lnm0TskCueyzLgXK/kMF3eEViVWO\nfH25vAPkefR+2pZDkq3uOBsPG0CM75f7fir+JMorYpsyTd9rYtgSBEEtQvLoAirA7BqGdubfVJ5N\nIUsc3oM5e6i6YRTXEtcXADBAMD2iU/l7jikLy8OhF0MjUIp6ZkqfRl5xyyEN4fTseche0r005s3e\nF5GfKYc63O/EWKzu2qJ2a1Gs1MN35U3zfQCm0yj7IBlXQvIIgqAWIXn0APYwK8i96glzDMittPUi\n1U7C4Hj6TJR7sW7H+8dcynE9JY9ZlCUC1olKS5Zls9luJ3mcgyx5eH2BX2hnp65/nNJex4zx/k/n\nUNY5sVxWr+HDX/K9nnLHT2E8DMA6oSeSh4i8WkQeEJEDInJDxfm3i8iPROSu9PfOXtw3CILh0bXk\nISLTAD4E4FUADgLYLyL7VPWbLuvNqvrubu/XZJYAHEnb7JHYw9kZFh7zi+e4T+3/TuRezIeA7FWP\nfLG5F8vne5R2uoQVtPd/smjyVP0OyB+fNeVnOViuv21z/fVinynvxTJURaUDWiVDP8N23B0f5aX1\ndemF5HElgAOq+qCqngbwSQDX9uC6QRA0mF7oPM4D8IjZPwjg71Tke4OIXA3g2wD+hao+UpEHIrIH\nwJ4elGsosOd52h1nL2aNnfxSbb+IbhZlOwzCHnTa5OEMRSd6keem9Gx3T6A8s+Fjo1ojt3bL8nn8\nKHJP7nU9Pq7MrCkHy8VyfqvNfSy8DvU3cyhLTZxdsr5bKSWxzKfQiiD3sjx3uIPyjDuDUpj+OYBP\nqOqCiPwagJsAvLwqo6ruBbAXAESkqVEN1o31m+FFehveAchDHCs2+8bHrkylUdl6LD/bBU9eQOvq\nT6AcnoHHF5F9fXh4fCPKq2hZfypt2VhtQNkQrY5oTAXtVlNW3sM3XFXhHlh/65vFD2mC3gxbDgG4\nwOyfjzz9DwBQ1cOqyo7sYwBe1IP7BkEwRHoheewHcJmIXIyi0bgewD+2GURkt6rSPul1AO7vwX1H\niqp1DhT/KUrb6Vyg6O38kIFi/Rkmz1qrXfkbO8R4OKXsdXm9GWSRnmEauM+hE+twGOXemlDy2G2u\nQwmEkgclJg4BnkR52MdyWqrqA1T7n/XP1A/FqoJrsww0UKvyFBf0oPFQ1SUReTeAW1FIhTeq6n0i\n8n4Ad6jqPgC/JSKvQ/FengDw9m7vGwTBcBHV5o7ixknnQeZRnprd7VLqPJZRDrDNHt2uAKWBFjXQ\n/D17fHIIhTbbQqnnwpSeh2x+TkMvv6KVZfoBgPvS9tdcHmrMLwdw7hrXYfkfRh7v0kuZ94Fxfiqj\nhZILp085hj4L5ZXKnEonW5D1IPz9oy5dMXma5ku1F6hqrTWPYZ4eBEEtwjx9wFT1XBzfU2KwszFs\n3amToMSxweSl9LDZpdRn2OXt7Nkp0tlpUqDQE3gTbu8nw862kAtRzSLK06Xtpmo3mnK08xV6NrJE\nxN9zn5KDnb3yU9wsM5/5UZQ9tfOc1W9MkvFXp4TkEQRBLULyaAA0OWePbCUIShWcqfA96Syy6fWM\ny+ONumaQzb6913T22kdQnl3xi7/YM08DuMKVj1hDM5rS857WqxqQZ12OmHJU2bPwt1PuHPNSOmM5\nbc9oJSsgSxVHzT2Z8n0EqxOSRxAEtQjJo0Fw1oHj8s3IPSZ1FJQCbAwZ76Gd1pDe3Hoj8kwKj3FW\n40FTBuobeB3qVwh1DeciSz3tJI8nkSULv5iPZaC7xW8gz5xw5om6FFsGv+jQSymUPKqsZfkcrQT3\npEuDzojGo4FwevJslEMR+HUwsyj7BeE/LqdC+Y93tsnLKUtO3VIh+RLkoQ3/Kfkbiqm83maUFa6E\n15s2+VkX1oH/0BymAYWxkL0n17hwinoZ+fkcc3m9Wb/Ny4bUewRbwviFRBgUMWwJgqAWIXk0EPaK\nR5F7Z4r+7EGfk1KGcATKCkiK9TSq2o5WX6oA8KyUviyllyIbl/lwk37KdRlZuejDC1g/pd536bRL\nOczYgawgfiylLAuDZAuyRMRpZ9aX5bXBpqncfSCllNioFD6GWOxWl5A8giCoRUgeDWa1JfZiUu9t\njD0pe21rZk7phD0wTbnPT+nZqFY82v1Zs08l72dc+a5L6W6Tn4pglo+6EKvw9SEpWU7qM6zhF3Uo\nNCP3eg1Bfk7ed2uvY+JMIiF5BEFQi5A8hkjd8IhA6/Jx6hto5OSnMNnDT6M8vqeOgQZWs+Z6P3Tp\nRSnlNOom5KnUn3HXpaQwgywRcHn9d1J6jktnTTlYLmL9irYL5Xnc5AGK59pN79jN+5kEQvIIgqAW\nIXkMkW56NI7drSk1dROcgVF3fBF5ZoLm2TTyso5+fpC2fz+lv5xS2n/Q0GzLdkBS97PsKjOdum1d\nAY4l6yv25N9O6e+m9LdTeq4pB3UdNNyiLYaivUtA1s0uPvQuDdZDSByrE/48RpwzkYcpnH6k8pPD\ngSqnyzzmQ0meRv6H5TE6IX5mutGOZLk1Nw8sJy3qoovBMJtagekZ4HT6bz6S5l8fTi0XnRpzmLTd\n1IXH2Dh6gzB7jMMqGrzZBjIcFa9N+PMIgmCghOQxRtAYjAZVHAKwN9+NPOTY6c6RaeRhDyWX7Ul7\nudG5Als4DhxNYoo38abCc9t2YH5z67mTSWR4Mv2IksOPUR4qcHhFCeRh5KlZnqPQw+nXFu/bwZqE\n5BEEwUAJhWmDoJLQev7ynO/yWGMndh803KIOgCboz0deaEZ9BnULNLiyPkzPTKIM9RdUfD6c0tuQ\nfZj+tCvnvSm9/Eng6pT/mUk7uzWtcqMkM5dEhUVkqYHTsVWe2xmC8kBK+Uyquk9KYZSwvA9XoGzW\n3i6QVdBKSB5BENQiJI8esg15zN7JFCFnDqhboIGUN6oC8kI4GmhRGbQTWXqgnsHfmybou8w9+Hvq\nEqizeDaAM5JCZEMST46n6YtvJgnin6e8n5kC/u01xfZWFizxVFJM3PVV4Lo0t/oH6ffPT1385iRW\n8H6bn8j6DNaFUhDLvcvUh5KHN4rjc92ILGlRKqH6hQvlgOzdnZIfF9NRF9PJlK2996Qs8Q/JIwiC\nWoTkkeCDYGtKnYL19L3szhEuQLOSB3s6ShE0xtqMbOTEa/tYLNRHHDXXoXTCmRDOMCyb61mHOUDu\nDfmbeZQ9hVM3YZfNzyYPP7TPeCQpVv57yvP1Zxfppa8F8JwUOXQnNSsFW58o5IKrr7gTX/98cezX\nk3XYb6brXZLuw/ttRDHjAuQYNKw3yzdv6vP8lFJy8PF8p1A2gqPE9ZjZ9xHtqG+xLgn43GgCT1sY\nvjMbeJy/8zYm1hu9jwTIc3yXo+CtPSSPIAhqMTaSR9UiJhuz1MJeV1Hu7b2vUOt1m/fg79lLcHy+\nCWXT6Z0uzwaU7ROsWz8gu/a7xNzf+woVk5c9Lcf3rPcOk4fl5T0pEfF5cbZkC7LJ+VNJIXJLOvdv\n0oUvfW06cNVbgAuSG6EdNF5PHHmoSHc+C5fi5uL3HywO/UnqVt+arr9le743y3G7K6f1z8rndLmr\n55LLewy5R+cz5nNMwhMWUfY+z3fI5zprrnPKneN7ZpmqemPv7X0Z+Rvy3wDfK8+fMMfafcdWShnk\nYr7GNx5WDFTkh80Xwn9+vlw+yNMoBxfy+1PI/1i8B3+/4vKeYfJscnmZZxblIQNfuJ1y5D9EVZAi\ne/0NKIvi9lmwTmyg7PDEplQALqAcdOoFKeWg44xthUk5ADydCvin6dwHrk4bHKpc8DLg/JcW29s4\nYZrYfHbeTkOY5199JwDgDV8pDl+Xrr9tV773s47msloWTOobx+0uD7+Rk8jPib/hc+QU7jLy+1s2\nv4M5vhn5nfDaPug3r29X/U65c/Yb4/dlOwGbxzZ6fpjrnTnbOvD37KD4rfEas+aeS+iukenJsEVE\nXi0iD4jIARG5oeL8vIjcnM5/TUQu6sV9gyAYHl1LHiIyDeBDAF6FwgZnv4jsU9VvmmzvAHBEVZ8l\nItcD+D0Ab1nz2ihaVLaUdkWlF0e579ZnAcgt9qxLp5FbfO9V2/9mDq0hHoGyUnUKZU9YlDw2mbxe\nevLltHXy0hKfBZ/DlLmHDzPJXoW9z4K5N1+8906+cQswlW5CEf2fpnQr7d+pHN1xcZY4tri5WnL8\n4p/k33renS3X4/V5v41bgB1HW8vlpYkF5GfiQ3D6YeoUyj28n9a1ylCWh/D4LPI78V7RfM89Zcrh\n72k9oPnAXKyT90Y/VXEPYn3BVv3OXt9+q2Ly+Dqvh15IHlcCOKCqD6rqaQCfBHCty3MtgJvS9qcB\nvEJEatnTB0HQDHrReJyHPLMGFNLHee3yqCpDZZyJCkRkj4jcISJ39KBsQRD0icYpTFV1L4C9ADAV\nq2qDoLH0ovE4hGwxDBTrlPyqaOY5KCIzKGYu1/TToijGanaGgWM4joU57vPGNXa2pSreCNA6Jm4X\no8TqN/w9/GzLdEUeauOtvoXXZF5fLo5DbR2qZop4De/D1Aek5nWt53E+Pxo7caZm1zFgOk1fcJx/\nY0rfk94qDcBw5KHWWRXL0YM5T8r/1KHW672cdUhdxMljOZIdjbiom7AzSKyPDbht62Sn0r0ui2N/\nmvBXzbYQvp9F5HdiZ3Ls9ayeY8md89/SEso6E/9dW8/17WZbquLpeB2HjzQ4hdbZlm56514MW/YD\nuExELhaROQDXA9jn8uwD8La0/UYAX9EmOxIJgmBNupY8VHVJRN6NIszoNIAbVfU+EXk/gDtUdR+A\n/wrgv4nIARRrsa7v9PrLqNY2t+t9qozEvGa2EyOx027/aaxtJGYlGfZslAasMZA3EmOv6nudJ9E6\nMwRkexHWyUZtY17Wm7247bX9vf7W5TnnaF4qf0aavnlDqsQ3bivSn7uimDVpMUk/3sZI7JH/Czxw\nZ8vv34DW69OV4dNH82I3luvvVtSBz4cL1/xz4/NYRmsvC5TjtlQZiXnJQc02JRDrUxXIEtJKRZ4q\nIzE/i+SXRXRiJMb7WA/6axmJ+aUV3dATnYeq3oJsiMhj7zPbpwC8qRf3akfVw/LiexU8x4ffbm3L\nSZRFwtXKwX9cTpPyo61a2+Kn11iWB1H+cPxU4QLKloYsAx0Zc0XuFMoOhlkGrnG5AsBZqYBbk/XZ\na1Lj8e/SDT6c1qpcipt/MiTxa1t+cvyBO/Hdz7f+/jfRen02HsdMOVguG9CbdWCDzJWxfm3LEZPX\ndwIcpn3bHL8kbfMZ85/SDm3WWttiA3uz0VjP2hY/vKha27Lad0wG6bQ51rYEQVCLxs22DIt2qxi9\nyXIV7GGsua/3qXHYpUA50LP352H9QlBEtwZzQLEKlWXzIRUfTCnnxK3yjb0g15LwXicBLKZucGMa\nVlyQbLl/Jcn6V6Zu+zN/AFxxTTEk2br7Tlgq/Xmg9XpzaSzy1PF8bypuORT0q00XkFfefiOlrC+H\nDqyLXY/E58bnSMkByO+kjj8P5rXGhe38eawmsfZyODEoQvIIgqAW4T29QXTrw9QtTSv5MP155JCR\nq/kwvTRt98SHKQCup6MPU66iZayXw2kK97so+zD15fwegL9M296HKbF+SsOH6dqE9/QgCAZKSB5j\nRMRtKYi4LesjJI8gCAZKSB4jTsSqjVi13RKSRxAEAyXsPEaU56X0LHOMVoqUGOjgh3YeU8i9Mu0T\nvNMi6/bu/Sn95ZT+QpIYnpkcMGzZnv2dLjtlxekkIuiKmaVJ525N6R+n9Lcr7k2pgjMf1KFsNXko\nVVDfwvpS4gJy9DzrmSroDdF4DJFunNXaYQev85TL4/1jzptjbGj4z8lRx3Zkpep7UsphwUOuvJue\nzNPJNNsm1JHOoBzEig0Wr89/9lnkIdMJl5eNx2ZTdi9rb3WpdTRch0E6Ex5FYtgSBEEtQvIYIt30\naNYPCQ2qvCk9Rf9TJq/vrSny0zR+M/IQhk5afGgDmocvIk+bfsZd97qU7kbZd+xPpdR7hD+JPJxi\nuShx2FXEfuEay0Vpxy5w8z5T1kNIHKsTkkcQBLUIyaPBnIPcO1OfQV2AmtTHk2HPTIWi9frO3/E6\ndD5LgytBVrD6eCNeaplBVtj+qju30eQh3jeJ98p1GHkx32PuN9SpnEDWi7B+LJf3oq4o+9tgHupF\njiHrdIL1EZJHEAS1CMmjgbAn3Ybs5Yom4xzDb0IZzqqwV6XegObax1B2mMPFZf8rpb8AgD7B2Fvb\nuDRAa3hMGplRQiDWWRGlBkpErIN3X/AQ8jQuTe23VfyG9Trt8rD+1rKQuhxObVfF53nMHQs6IySP\nIAhqEZJHA2FPamdSvA7BelVfduesr0wg60usHcgpl4e97n6ThwvsNqAVaw7vA3AT6jGeRKu5uIVl\n4IK2/aYcLBelArvwzpvj+3vbKGx8lpSivFf7GZPHOggK1ibWtjQIKh83m/Qnq1PdObvGw4eP8KEM\nOf35GLLFpQ9zaFfeXpi2qTjlcMD7hD1pzrVrPBbMOQ4hZsw5ICs+v4/yillfzrOQQ1Ge4c7Z0AZA\nq1jN63AIddSkx905PqNJIda2BEEwUGLY0gCodKTkwZ56BeVARtYDN1D0tpQEvLGU9/+5hGxiTpGO\nQxNrEMbhAKUdPwzivY8BuCtt++lOmpw/G3low+v5EBaUGHYgD3G4mpZDGjF5WB8feImSAz/qM1D2\nNm8lIqA1gDmV0KwnyxBUE5JHEAS1CMljwNhwiWy5OXZnj0/JYRllj+h+SpRRw4Fs8MXf+0jij6E8\nHemDD500x9ijV3kb82X5vstjF7tNu3PtQiKeROsUr4Xlfgzlj5Y6E0oelKLmzL18/Bw7He7rc9zl\nXTF5fLkmmZA8giCoRUgeA4Lj/e3mmNdRVM1csMf1fjjICnIP/qg7RwMwjuWtV3BOXdIgjNLPKWT9\nBaUASkTsmU+alFKJixfX4hGMU6BeymGdKDmcQn5OLBfrQCngILIuop2Xc87GnELuHX3Qb9bJelnz\n74M6KPtPwrq0i80ySYTkEQRBLbqSPERkJ4CbUYQD+R6AN6vqkYp8ywDuSbsPq+rrurnvKGI9dvGh\ne/sM7lvTaT+74hfGTSP3mD6INanqoenwh7oJSieKbLxFn6Asg7cfeQrZzN073eHxaWQdgrfH8PGA\ntyPrIja7PNYTWDuJw9uj2PjCfraK6ZTJ49+HDWDOa05aTJfV6FbyuAHAl1X1MgBfTvtVnFTVK9Lf\nxDUcQTCOdKvzuBbANWn7JgBfBfCvurzmWMIe1Uaqp/Rg47ACuaeeRdnS0s++2Ijr68E7ybE9Mu0e\n/KwLJRLbs1PCaBcrZTNyfb1pPfUNlBhsfdnD13HmQ/3GAsoOfShN2GfWboaHz0RMOWqZYo4p3TYe\nZ6sq9XQ/RNZVeTaIyB0ovrsPqOpn211QRPYA2NNluYaGX1/CD5D/kPafwU9h8h+G/6RV4jKVhd2u\nw2DYA1/OOXOMKcvgp1gF7Q2pOAzahnLQatbbf3xLyPXkmpZvoXP4z++njYFWRTWQG4wZ5Hfihy3T\n5jjP8T1yGpzPgsOidgHTx5E1Gw8R+RJaHVKT99odVdVV1qJcqKqHROQSAF8RkXtU9btVGVV1L4C9\n6d4TtbYlCEaJNRsPVX1lu3Mi8piI7FbVR0VkN1rjLttrHErpgyLyVQAvQBHXeKyYQfb3SYUfeyoq\nlwRlUZqiMHs3TmEeR+7ZvEFVr6C5ujVp53DCDrWAcrntEMwza/J4vE8N692MksxD6C2U1KoUp6y7\nr7f1IeJWlJ+LAAAPTklEQVSHWt5rG4eaRzA50ke3CtN9AN6Wtt8G4HM+g4jsEJH5tL0LwMsQYTSC\nYOTpVufxAQCfEpF3oBhqvhkAROTFAN6lqu9E4Sz7oyKygqKx+oCqjlXjYXsx70/UL1VfRO6Z/Bjb\n5gH6EyqRPSfLtc2dX0TuTSkRcP+kS59G++XrPD6H3Ct7XQqlMvb0NmIlPad5w7Rj6I4Fl7KMFtbb\nvifv7oCw9+X1NqAsLY4rXTUeqnoYwCsqjt8B4J1p+68B/Ew39wmCoHmEeXoXWI9fQKu37tX0Gt7Q\ny/dmdeOFrBXhbB6Fsx8g62YoIXHcP4uyxEE9BCO+Vc1meH7kUiA7GaKBGl0RUAKZRfmZMg+nhB9G\ne71P3QhvfhaJKSUhQVkf4t8z9+dQNlYbV1P2ME8PgqAWIXnUgGPkTW5/CWWnM96mYRFZB8BZFW9E\n1S8T6AXk3oL3mnP71ghqyqXdzpvz9/66RCrKddrl7ceSeM7E8F6UGOzCQEoTUy4P3zP1V3b5gbcB\nGTcdSDQeNagyrAKKhoIfkQ8zYP12UgRfMsd6QSfiOhWOVYpcoBhCsDHzZecqU/4TdBosiUZCfpUq\n77ls9vncfCDuThSldYd7ayk4jyCX2Q9J/PqceZSN/za6vONCDFuCIKhFSB5rYEMesEfxU4w2KJLv\nkXyPuYDhGhGxtz/ljtt1NX4owx6d9WXv2ynMz9/79Tr2mS27cyynncYdNEvIQxtfd/9NzCLXk+lG\nt7+M8rBnFAnJIwiCWoTksQZ8QBvMNltc6gTsWNb75DhRkWdYzKA8NWs9tQNFj0iFrjeWohS13rgm\nzM+el9OeXkk7i7Jikj37KfObYUhu3kCO0hnrwudnJSQeYx3sQkMv+Y0iIXkEQVCLkDzaUPVgfEvr\np+JWkMfq3hdGU2DPTgmEPSfrdhrlqGoc7/8gpeudsqX0ReOyc1Pq47lsRjnWDKWSbs3Se03VuwcK\nKcNLTT5QuIXf2SgupgvJIwiCWoTkkWCP7BdJkWWU7Qg4bvWeuZvKErJE5MtqPwTr7AcoL5Srol1s\nFos3e7e6DqB1YaGfkWG5m9ZDe10IkJ8FPdL7QOH22Xh7m365XugHIXkEQVCLiZc8+ADaxXe1ugse\ns/oBoPkSB9mIsn0Be/0lk/Jcv6w7eT1en/qXKbQ37Wb5NqL5thF+oZ2fnVtCq0WqzeM91DdN0rJM\nfOPBl+gNvryhzyLK4nzTlKHtoFn4JpTXjHB9zbQ5TkUp3cL1WoTm9ej8lv8o20z5/D8gxftnIDc6\n650yHjRs5PzKW7uGxytT/XfY5MYjhi1BENRioiQPr6BbQdkQikY+/vgoKLA8fvXvFpR9h3gl5gJy\nz97v4YEfHs2hvBhv1qVWoe2lqKbBZ1v17fjhorh9/nYGuYf33+awCckjCIJaTJTk4f1m2OnXdj4X\nmq6cWw0/7bcR+Rn4qVXrIYuSSrvwlb3CG6zNob3XdbtgrqkGeOuBSmOvaPahL6dRfldNqXdIHkEQ\n1GIiJI/VzIS9Vn8c4ezJHMqexKq8ebH35wxHvyQPXp/3E5R9g9rpTaDQj4yTT9B2kq2djWln7j5s\nE4GQPIIgqMXYSh7TKFfOB022sy3jjLddAco9OjX5p5FnP55Gf+H1KXksIs+gtNN9CMo6gHHESmA+\nIPqMS5cwHClkbBsP+5Hx4beLhj6uUCFpV63yGJ+JDeAEFGEWBjWE4wdvA0Qx1IIPjm0bPt/gjeP7\n9B7VgPL0tVTkGSQxbAmCoBZjK3lYP5Gkyaa+/YAKSUobMyh7RvfKt2Eqjq0PUx+ywgbbZn36rdBt\nGvx+201jD5qQPIIgqEVXjYeIvElE7hORlRTcul2+V4vIAyJyQERu6OaepWujdTxMFEVPZv8mDRrB\n2WdwKv0tISvaltH+OQ4aloPlYjlZbluXKh8rk4D/rquU/oN4n91KHvcCeD2A29plEJFpAB8C8IsA\nngfgl0TkeV3eNwiCIdOVzkNV7wcAkVXbuCsBHFDVB1PeTwK4FsA3u7m3X+TWtEVDw8S/DeoEFGWP\n5U1b+Gfjx9h9u1DRS5E++PQkM8j/i0HoPM4D8IjZP5iOBUEwwqwpeYjIl5DDjVreq6qf63WBRGQP\ngD1r5fNmzNHrZOgVjb4z7YyK74naRbYbFr4c3s5jCuUobaxLU+owTAb5f7Fm46Gqr+zyHocAXGD2\nz0/H2t1vL4C9ACAia9Z9EhVma+GDDFU5N/YuFJvyHH152Gj4oOJAObRnkBnE+xzEsGU/gMtE5GIR\nmQNwPYB9A7hvEAR9pNup2utE5CCAqwB8QURuTcfPFZFbAEBVlwC8G8CtAO4H8ClVva+7Yger4acw\nZ9PfnPmr8hPRRFhOW3bWh0zqlO2wEdXmags6GbYEZRgFjutEbJR275Wca1qapi9gHRj7hHWYQ3ld\nzhMpbVodRgVVrWUSMrbm6ZPMokutibfXEzTVUxrLRVN0u8rWB5W2waWDwRHm6UEQ1CIkjzGEMxJ+\nOlZRFvWbCnUYh1PKumxE2SiM9W2KodukEJJHEAS1CMljjPCBlX0PvYjm6jjWguW2bgVYPx9Iut8e\n0IKCkDyCIKhFSB4TgPXd2oRl93Wo8r8aDJdoPMYIH5Xdm3grRlfUZLmts1+/jiPM1AfLqH5LQRAM\nmZA8xhAfMMn20KfK2UcC+vdYRjmURJghD4eQPIIgqEVIHmMAjaR8+EYfxHtUp2mB1vUsrCdT72nd\nG8cF/SEkjyAIahGSxxhAyaPKcxgwfj0w6zfn9pmGk6DBEJJHEAS1CMljDGjnum/FpeOCr5eP2RqO\ngQZDNB5jgPfN4YctJzBesD4Um8d9mNZUYtgSBEEtwg3hGOCNpsikBPb24rM3Ww9Wp64bwpA8giCo\nReg8Rhx6QgeyAnFSe1w+B7swMPQf/SMkjyAIahGSxxgw6ToPbyTmXRIE/SEkjyAIahGSx4izgnIP\nOykSB2F9vcQxbsZxTSMkjyAIahGSx4hjI6ixB560sb63LLUR8oL+EY3HiOKVhMDkNRoe1n/OHPOm\n60Hv6KpxFpE3ich9IrIiIi9eJd/3ROQeEblLRO7o5p5BEDSDbiWPewG8HsBHO8j791X1x13eL0hE\nKIL28HlMYXRDTYwCXTUeqno/AIjEKwqCSWNQOg8F8JdpodtHVXVvu4wisgfAnrS7gEK6aQq7ADRC\nekrK0V1LDSlPohHPxwW8bkSZDE0rz3Pq/nDNxkNEvgTgnIpT71XVz3V4n59T1UMi8gwAfyUi31LV\n26oypoZlb7r3HaraVpcyaKI8q9O08gDNK1MTy1P3t2s2Hqr6yroXN9c4lNLHReQzAK4EUNl4BEEw\nGvR9KlxENovIGdwG8PNo1lAkCIIadDtVe52IHARwFYAviMit6fi5InJLynY2gP8jIncD+DqAL6jq\nX3R4i7a6kSER5VmdppUHaF6ZxqY8jfYkFgRBcwkL3iAIahGNRxAEtWhM49FEU/d1lOnVIvKAiBwQ\nkRv6WJ6dIvJXIvKdlO5ok285PZ+7RGRfH8qxan1FZF5Ebk7nvyYiF/W6DOssz9tF5Efmmbyzz+W5\nUUQeF5HKiQEp+GAq7zdE5IVDLs81InLUPJ/3dXRhVW3EH4CfQmGw8lUAL14l3/cA7GpKmVCsvfou\ngEtQrMm6G8Dz+lSe/wDghrR9A4Dfa5PvWB+fyZr1BfDrAD6Stq8HcPOQy/N2AH84iG8m3e9qAC8E\ncG+b868B8EUU1vMvBfC1IZfnGgCfX+91GyN5qOr9qvrAsMth6bBMVwI4oKoPquppAJ8EcG2finQt\ngJvS9k0A/lGf7rMandTXlvPTAF4h/VvDMMjn3xFaGEA+sUqWawF8XAtuB7BdRHYPsTy1aEzjsQ5o\n6n5nMmUfNucBeMTsH0zH+sHZqvpo2v4himnwKjaIyB0icruI9LqB6aS+P8mjqksAjgI4s8flWE95\nAOANaYjwaRG5oE9l6ZRBfjOdcpWI3C0iXxSRyzv5wUD9eQza1H2AZeoZq5XH7qiqrhIU68L0jC4B\n8BURuUdVv9vrso4Qfw7gE6q6ICK/hkIqevmQy9Qk/gbFN3NMRF4D4LMALlvrRwNtPLSBpu49KNMh\nALYnOz8d63l5ROQxEdmtqo8mMffxNtfgM3pQRL4K4AUo9AK9oJP6Ms9BEZkBsA3A4R7df93lUVV7\n74+h0B0Nk55+M92iqk+Z7VtE5MMiskvXcKExUsOWhpq67wdwmYhcLCJzKBSEPZ/hSOwD8La0/TYA\nJclIRHaIyHza3gXgZQC+2cMydFJfW843AviKJs1cH1izPE6f8DoA9/epLJ2yD8Bb06zLSwEcNcPR\ngSMi51AnJSJXomgX1m7sB6WB7kAjfB2Ksd8CgMcA3JqOnwvglrR9CQpt+t0A7kMxtBhqmTRrz7+N\nonfvW5lQ6A2+DOA7AL4EYGc6/mIAH0vbPwvgnvSM7gHwjj6Uo1RfAO8H8Lq0vQHA/wBwAMWShEv6\n/J7WKs/vpu/lbgD/E8Bz+1yeTwB4FMBi+n7eAeBdAN6VzguAD6Xy3oNVZhcHVJ53m+dzO4Cf7eS6\nYZ4eBEEtRmrYEgRBc4jGIwiCWkTjEQRBLaLxCIKgFtF4BEFQi2g8giCoRTQewboRkem1cwXjTjQe\nQQkR+WxaeHgfFx+KyDER+f3ki/YqEXmJiPx1Wkz1dRE5Q0QuT9t3pUVol6Xf/oo5/lE2PskPx9+k\na3x5iFUOahBGYkEJEdmpqk+IyEYU5t9/D0Wgoreo6qeSGfi30v5+EdkK4ASA/wjgdlX945RnGsBF\nKNaSvF5VF0XkwyisGL+IYkHW1ar6EO856LoG9RnowrhgZPgtEbkubV+AYoXlMoA/TceeA+BRVd0P\n5IVVIvL/ALxXRM4H8Geq+h0ReQWAFwHYn5ZPbESxoO+lAG5T1YfSNaLhGDGi8QhaEJFrALwSwFWq\neiKtyt0A4JSqLq/2W1X9ExH5GoB/AOCWtPxdANykqu9x9/mH/Sh/MDhC5xF4tgE4khqO56KQEDwP\nANgtIi8BgKTvmEn+Qx5U1Q+iWPH7fBQL+d6Y/K/QD+uFKIYuV4vIxTze95oFPSUkj8DzFwDeJSL3\no2gkbvcZVPW0iLwFwH9OepGTKKSVNwP4VRFZROHp7N8n3cm/RuH9bQrFys7fUNXbkzL2z9LxxwG8\nahAVDHpDKEyDIKhFDFuCIKhFNB5BENQiGo8gCGoRjUcQBLWIxiMIglpE4xEEQS2i8QiCoBb/HxD/\nIcoWiOPBAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x11591e518>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"osys = poppy.OpticalSystem()\n",
"rect = poppy.NgonAperture(nsides=5)\n",
"osys.add_pupil(rect)\n",
"osys.add_detector(pixelscale=0.05 * u.arcsec / u.pixel, fov_arcsec=3)\n",
"osys.display()\n",
"plt.figure()\n",
"psf = osys.calc_psf(display=True)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAEICAYAAACDNvdHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXmcbVdR77/V3bf7zpmBkJCJBDDhgRIGEfQhTzT6wKCA\nBJAZ8xgi+hQRJwhBFBwQmQmjTEFAkaDxAREDIlMSZDAMGkKAm9wMNzd36Ntzn3p/rFV3115n7TN0\nn77dfXr9Pp/9OefsoXbttfdZtatW1W+JqlJQUFBQUFCwOhhZbQUKCgoKCgo2MoohLigoKCgoWEUU\nQ1xQUFBQULCKKIa4oKCgoKBgFVEMcUFBQUFBwSqiGOKCgoKCgoJVxJo2xCLykyLynQHKe56I3Coi\nkyJy3ADkqYgcEpFXDkK/9QYR+aSIPGWJx75KRO4QkV0icoaITA5Ip10i8oge9x2L9/C0JZznESJy\nXb/HrXUs555uBHR7ZkTk6SLyz0dWq/UDEdkqIv8kIvtF5LIlHP8+Ebl4BVQbOETkvGhrWl37JFXt\nuAA3AtPAJHAr8G5ge9x2DvBJYC+wD7gW+IW47RFAKx5ny8e7nW+lFmBTvI77D1CmAmcewWt4BLBr\ntdpwgNdxOjAFHN+w/XPAM5YoexfwiB73HYv38LTVbpMVbOtd7v97C/BOYNtq67XUewP8MfDuZJ9D\n8fr2AFcCj+/2TADPAa4alF4DvvYzQ9c82LbL7PerwPdj2/09cHSP8p8EfCH+h6/MbH8A8JW4/Wrg\nfm7bM+OxY/H3z2TsxFM6nPt9wMU96vk5YMbJvc5t63he4DjgY/HZuhF44jLuR9c+qVeP+DGquj02\n8AOBP4zrPw58CrgbcBfghcABd9zNqrrdLY/p8XwrgbsCm4G+vRgJWFb0IMo4YTkyejzP2EqfYwA4\nFbhNVfestiIbBD8f/78PAh4K/F66g4iMLPcZ7wcDfk7Pidd3H0JH/RYR+YMByh86iMj9gDcBTyH0\n3/PAG3o8/A7gNcCfZ+ROEAzYu4BjgMuAfxCRTXGXU4HvqOqCO+wHiZ14/1KuqQHPdXLPSbZ1Ou9b\nCEb4LsDTgbeJyH3iNe4UkS0D1LFnj/hn3O8/B/4ROJ7w1pV9i6IP7w34BeCbwEHgJuBFORlRl98B\nvh4b6R0EA/vP8dgrgWMy8u8V91fCm8+n4/qfILyx7Y+fP+GOuQp4JfDvBI+izfOlB48YOAN4edT9\nd936DxM8lP3AZwmdSWN7ANuiHv4t7u7AxcBHCB3QAcKb/gTwWuDmuLwWmPBtCvw+wYO4kfgmSOio\nbwVGnS6/DHytwxvnM+L3EeClhDfs2wiRk52ZY85LruPtOA8AeDWwSPUm+9oo+3VR7v54/8/u8Pb5\nW8A34r6X2bXH7S+J7X4T8Gyc50DiidPBayK8Ud+YnPdFwH9GvS8lPJufiPflk8T/Sryej0Q99hGe\ntR9xsk4A/ike92XgT7wewNmEZ30v8G3gcb2+jQN/BfyDu95XEDyUaeC05J4+B/hMbPt9wPXAQ2K7\n/TA+K7/qZP8i8NWo9w+AP3Lbzoxt/cy47dOxbZ6X6PtN4DH05xGn+1wQr+foXBt0u7dx+3Pj9d4B\n/ANwYnLOXwe+R/gPvQoYycntdK+ArfF+/ICqH5gg/Getr5oEHpTR76HAF+N92R3v0aa47fPUIwVt\nzwfwZ8B73O97A7NRp7OivveP2+4R2+EnM210ZbLuF4Dvu99C+K/9DKE/nSMY/UmCgav9jzJ6nhuf\nqYOE//KHcR5xfOa+Ftvhc8B9c/1Tt/9vsm1n1PEMt+4y4I/j94cDdxKM9YObdG/6D2b36UHIjURD\nHG/IdYQ/rwD/TTDKjwXumhz3CHo3xLvtJhPeoh6QkxF1+SKhgzuJ0DF/Bfgxgrf7aeBlDec4LT6c\nFhI5NjbmUwl/rifF38fF7VcR/iDnxO2bMjKzhjg+zE8D/jU+wG8GHpLs8yxgB5XR/Gq/7RHXXRwf\nmscSOvgtwCWxne5C6NQ/D7zCyVggvNFOAP+T8Ie9d9z+TYIHZfI/Cvx2Q5v6TvtC4L8IYecdxLfi\nXv4EJKE42g3i/yYYpKPiNZ4N3K3DQ/9Fwlv+cVGn58Rtj45tezbhxeZDDNYQ/3ts85Pjfb8GuD/h\n2fwM8Adx3xHgGbGdNhM8kWucrI8A74/38r6EjuyquG17/P00wnN5bjzXvbt1AsApBGPwMne9NwI/\nQhi6GaPdEM8T/iOjBIPzfUKnP0HodPcDW+P+jyT8X0bide8BHu3vMcFT2hqv7cnAvztdzyX8p8dY\nniHeTHjRe1RTR9jl3v5s1ONHo6w3Ub282zmvJPw3TyUY7GekcrvdK+CtwL8AJ8b2fXi8D11D04SX\n5odEuWcQnvOLOrVLcvw/kfyvcUN3wPMIL5Vboo6vysjIGeLfIRmCBP4f8Bvp/XP/oznCS90NwF+6\n52ki3rsXxna5gPA8Xuza4Nb4OUroU78LjLvn+3bCc/g54Kd6PO+DgIPJNbwE+Kj7bQ7W9wg28UV0\n7pMekdt2eJ9OG6OQGwlvL/sIf8I3AVvitpMJnch3CQ/+Z4Gz4rZHxHX73PIrDef4AfB/SDwo8obY\nx/H/Dniz+/3rxLf9zDlOo26Inwp8OdnnC1R/qKuAS7q0TZshJnh4e4ErgF/BeWMd5BwdZR3VT3vE\ndRcDn03WfZc4Vh9//xzRaFAZ4m1u+4eI3gvwu8D74/djCeM8Jzbo7TvtzwAXum3nEN6wRzLH9WuI\nf5ZgQB6Sk5d56C9wv18DvCF+fw/xrTb+PpvBGuInut8fA17vfv9f4CMNsiy6tI3Q4SwA93TbX0XV\nuT8F+Nfk+HcQjXxDe/j/7xuAze56X9rhnj4H+Jbb9mNRz+Pcuv04LySR9Qbgz/09Bk5x27dEvU6P\nv18LvC5+N2NygHofMkMXQxy37bH7Qf+G+G+AP3G/dxKiNCe7c/oo4QuBT6RyO90rguGYxUXCmv4P\nvSwEQ/Dhbu3i9v8M8QXVrbsVeLj7fQUhsvQfROOW7J8zxC8H3pes+1vgD+P31BCfSHgRHAHuSXiZ\nfWPc9khC5EXc/l+mMsRvI3G8CH3fw+L3Hye8DE0QjPRBqv96p/P+NO397PPSa43rhdCnvpvgyF0O\nnJT5Dz4iPdYvvY4JPVZVj1bVU1X1+ao6TXhSdqnqRap6T8Kb4SFCZ2e4OR5ny4ca5D+OGNIQkc+I\nyEM76HKr+z6d+b29x2u6O6Fj8vg+wdM2/LBHWR73JbxpfRX4uqrOpjuIyGjMGv6uiBwgvGBA6JCh\nv/bI6Zle2/fjOsOdqnqoYfv7gMeIyDbCi8S/qeruLudvOuc4wSNfFlT1k4Qw0JuBW0XkLSKyo8Mh\nt7jvU1TPxN2pt1V6/5eLnp7NeP//TERuiPf/+rjP8YRoz2iip/9+KvAwEdlnC/BEQsfShEe7/+9F\nqjrTILuXa1pU1TsaruuhInKViNwuIvsJRul46jh8vtiPfAT4VREZJXg87032v5/vQ4C/6KIvIrKZ\n8BK5N65aILzgeGwieFc51J5lVT1A6GSb+ob0/2XodK/uSvh/fLfb9eQgIveJ2ce3xGfoEtrbuhMm\nCS8YHjsJxsrwNkJ/9jpVnRug3MNQ1d2q+i1VbanqdwmOwOPj5rsTDKK6Q/x/9lTgd5P2PZF4n1T1\ni6o6qaqzqvpO4EvAz/dw3p6vIer2TUJ/fxPBAdmWbZkOGFhyhqr+EHgj4cb1e+zVqno+Iaz3DwQP\nbaVxM+FGepxCaMzDqvUrVFV/nPBGNQZ8WkS+LCIXJeVSTwbOJ3hVRxG8dQhvV53ao0mfdH16bafE\ndYZjoqFt266qNxEiA79MiBqkHWMTcuecI4SG+kXbdarqa1X1AYTn62zCOHC/2E0YXvE6ehwihE0N\nd1vCOXrB0wgvWo8k3P8z43ohGL4WwfsyeJ1/CPxL8oK7XVUvWqIufT/jHfBBQpTqHqp6FCE6JLWT\n1TtVCN7nUwhRjztV9eoB6PFYgrdpsn5A9R8znE7zi1jtWY4vfcdQ7xvS58j/vwyd7tWthP/HPTPH\n9XJP3koIHZ+pqjsJ+RnW1r0cfx1h+AAAEbkXwR78d/y9kzB+/Q7gEhE5pgeZObkC/A96T5JVquvY\nTf1/APX/7A+Blyftu7WDw+dld9r2HWCLiJzutt/fX4OIbBaRXxGRf4z7/yhwEeF+/FfnS2zHkg2x\niBwjIi8XkTNjxuXxBPf/i33KGReRp4jIUao6TwhFtZaqVx+4AriXiDw51gY+kdDB/+NyBcc3rRcT\n/qwvJ4zD3igiz4677CB0FHcQOv4/sWO7tMetwHEiclQXFS4D/lBEToj35aUET9fj5fFcP0kYO/2w\n2/Ye4MWEP9Df93jZlwG/JSKnxY7rlcBlqrqUe3krYQwGABF5cFzGCMZyjqU9Ix8CnhW9iW3Ay5Lt\nXwUeJyJbYsf0rCWcoxek9/9wHXq85/9AuD9bROQcQpmJ4XLgnPjcborLg0Xk3iukaz/YAexV1RkR\n+XGCh9sNnyN4hq+m95e+LETkOBF5KvB64E9VdV/c9LeEZ/NesXrhQYQx+g82iLoMeLaI3C9mAf8p\nITK0y+3zYhE5WkROIYSm/zYjp/FeqeoiIZz5WhG5W4ySPCxmF98GqIickZFp2EEYFjgkIj9CGMoC\nIMq+A/cfyuB9wGNF5Cfif+ESQmh7Km5/PWH8/jmEypg32oFR180EZ2MkGiXLgv80MCoiL4ht9xuE\nyMNnckqIyE+LyD3i91MIbf2xuPlzUf5FsY/+FULljuFtwAtE5EHxvm4XkceIyDYROVZEfjbqtklE\nnkZIcPtkt/PGCMjHgFdIqHv+SUKeyvvi/g8gvHhdRPXi+QxVvSrzotkTluMRzxHeMq8kGIv/JHQu\nz1iCrKcSDNUBwrjDihMKxPDao4HfJjy0LyaE8AZWUqOqi6r6T6r6BMIb9hfipvcQ3sZvIoQ10peX\nbHuo6rcJncQNMRSTC4dBGIe5hpBd/A1CQtsfu+23EEJtNxOSgp4bZRs+GvX9qPtjGlnB1xrO+TZC\nZ/RvhOSHg4Q/of1xJ6V7iN3wWuBJ8RpfQxhDfwdhjPBGwpvya6LsPxKRj/ciVFU/TuhQPkNIbvlU\nsstfEN6MbyPU29ZeXkTkO/GFbbl4F1VG+3WEZDqP5xESzW6N+15G+G+hqvsJY/6/SmiHWwidyETU\nsef2WAE8D/hTETlIyMrvGtmKHdd7CJGOpZatXCeBEOa/CVnZv66ql7jtbyEY+SsIxuvdwItV9coG\nnf4fwTB9lNDGp9DeJ32c8OL2H3G/d2fkdLxXhLyBbxH4F/YSXshFVQ/G/b4U/wMPzKj524Ss44ME\n7zh9EXgZ8IF4/C9ndPs6wZB8kPC8TxBybIj7PxJ4Qdz9N4CHumf/mYQhidcTon/ThDYmDnucTxiW\n2Bev/fz4gpnDA4EvisgUwfB+JbYLcVjvl4BfI/RXv0R4SbVr+CLhmXtz3P5fVC+tmwjteXtcnhv1\nsGGgxvNGPJcQjr6d0A9c6PrI3YRM9p9S1Xep6rLJiGSJBrwAEJEZQgf5OlX9o9XWpxdIYHh5n6qm\nIZ90v+8C/6epsyo4chCRvySU4jy7687rECLyLOBpqvqI1dZluRCRCwmEIj+72roUrC5E5OcIL0gT\nwM+p6meb9l0P5A9rFqq6ebV1WAmIyOMInuGnV1uXjQgROZuQsPWfhEzxZxLGlYcOMSz6fGKEYwhw\nDqGkpWCDQ1U/QYjmdcWa5pouOPIQkasIoZ4XLHF8t2D52EkIwR0ihKVfparLzl1YaxCR/00Ii/6A\n/BjruoKExJ1Hsg5fKiTwIn9HRK4XkZdktv+UiHxFRBZE5PE5GQVLRwlNFxQUFGxgSCgd+y/gUYSa\n16uBJ6nqN90+pxFeEF8EXK6qHznymg4vSmi6oKCgYGPjwcD1qnoDgIh8kJBwddgQq+qNcVuJkq0A\niiEuaIOIlDDJGoQk3zsVRGryu2DtQVWbbmFPOO+883TPnt6KPK699trrCKxkhktV9dL4/STqBCW7\nCLkJBUcIxRAXFKwDCBU11GhchJDk4VkcWvFzMS4QijiLMR4+7Nmzh2uuuaanfUVkRlVzZVAFawDF\nEBcUrCGME6i2bFaEzYTaB5uRgfjdG2IPb4itcHMhfp8luETTBN7P/QQygIL1CoXWQvfduuMm6kxh\nJ1NnEStYYRRDXFAwYIwT/lij8bcZxjkqLzWH7QSC5B3xuxnicScTKo/Y5Ke8hovUPeKFeO45KkM8\nGeXvjd+bMBrPbYbf5JvMglWEKizMdN+vO64GzpJA6XgTgRHtyYMQXNAbiiEuKBgQbKZw82DNULYI\nxmuMYPSawsRb47LNfZ+g3RCbke8WmjZfaYFqXkE7zvaxSZ9zkKiDGWLzvs3btuubbji+YKUxGI9Y\nVRdE5CLC/NCjwDtV9ToRuYQwPeflkRb0owTO7ceIyMtV9Zxln7wAKIa4oGAg8MlTI8nvUYLhG4nf\nm7rOEXds7tMMoclOPWI71yL1sHUvMnMYdfuMuvXpC4AZ9oIjDB1YaBpVvYJAAerXvdR9v5r2CRgK\nBoRiiAsKCgrWJQZniAtWF8UQFxQMAEp9zNQnVPkwcKduc84ts9Q9VqXKmrYQs/dKvR4WCrdzzVON\n6c7GxZ+rCQtU9S7mZftEMDu2eMOrhAF6xAWri2KICwoGhIXkE/oL2+4lGFHLcN5CGCO2MeflGGKT\nOUsY0z1EyJreR2fYMf1eS8ERQjHEQ4FiiAsKVhD9Gq59VMbRDPEgkrW8p73U5KpihNcYtDWorOmC\nVUYxxAUFaxTTlIzkgk4ooelhQTHEBQVrCEIg9NhG5RFvJvxRlztGbJNn+9B08XLXMcoY8dCgGOKC\nggHDE3qMUBnGWToTeuwgFGnupDLE3Zi1coa4F2atQ/E8dwIHO+g0Gs+fXksh9FgjKIZ4KFAMcUHB\nADBKM6EHBEM4TqCWzBkwIbBp7XCf6RhxyjWdqwNuURlLzzU9l8ixxKt5mklGxqkIPTa59WbkLYlr\nms4vGAUrhOIRDw2KIS4ocPAeZj9hW2OugmAkx6hCx556crT90MPHbHKLyRhz61KKy5SUo0U9gSu9\nFjPQuXPluvNRt68vxZIoyyeP9WOIl9rGBSlKstawoBjiggKHpRqGOSqvMWXI8uHc+ezR1Ta/LJAf\nC7bwsxlbb4h9aNp7xCav6Tw5zLv9W+5a1Mmya+8HxfgOCMUjHhoUQ1xQMAAoFWezebEp1/QcnY3Q\nfio6zBbB0OVC0yPxHJ3Klxbid6hC0758aRI4EM/ZhIW4TxPXdDEBawDFEA8FiiEuKBgwvLfYD6bi\nsoXO0yBayLrTNIheh6ZpEHspjVKqseCCNYbiEQ8NiiEuKCgoWJcohnhYUAzxkENE7gG8B7grwcG5\nVFX/enW1KuiEQuRR0BuKIR4WFEM8/FgAfltVvyIiO4BrReRTqvrN1VasII8thBKmpvIlGyO2TGYP\nT2/ZbYy4GPt1DtWSNT0kKIZ4yKGqu4Hd8ftBEfkWcBJQDPEKwQxnmjXdC6HHsfHTCD22UBnh5RJ6\nzFF524cIRB576Z/Qw2aaKoQeq4wyRjw0KIZ4A0FETgN+DPhSZtuFwIVHWKU1B5+B3A+M/MK+p7W3\nZgibkqTGCNSWtnSjuBTqmcwGM/qWsAV1isvNcRmPxy5GfXLduSWNpSQg6TSITSQlnbDUdi7wKIZ4\nWFAM8QaBiGwH/g74TVU9kG5X1UuBS+O+G7J/tHAvVAYHeiOr8KFjo7g0g2k1vxAMXs4Qj3dZPGHI\nJoIR7RSabrlt3jP3htTLz3Xnfrs3xAvuu+3XiyEedZ+mm+fELugXxRAPC4oh3gAQkU0EI/x+Vf37\n1dZnrcIMrg8n9wqjloQqPLxAuwFskrmY2c+XIll9sdcT8oQe3ivGyfDbFpPfTTp579pfix2L++wG\nO08asi9YIkpoemhQDPGQQ0QEeAfwLVV9zWrrs5bhw7n9YorKqCyQ967nqEg/UswSxm0nqLxFI/UY\nJ3iwKcWlZ++y/VNja/rYmK7VEU/G5RDNdcKmq70IpF6spQn1W2dcTMeAUJK1hgbFEA8/HgY8FfiG\niHw1rvt9Vb1iFXUaSswmn/3idoKRtIkfjNDDQsO5SR96McRGZZka4oM0vxgYum0vWE0Uj3hYUAzx\nkENVP0f7UGLBGoV5qgUF3VEM8bCgGOKCgjWIrYTMaV9H7Cku/exLHn6c2VNc+jriQ4RQesE6Rxkj\nHhoUQ1xQMGCY0fS1t2YMOyUnbSfUEe+kKl/KcU1bBnOnrGkrmYI81/QhwoQOe+nsgY+660lnkip1\nxKuNYoiHBcUQFxQUFKxHKMUQDwmKIS4oiDAvE+p1xL3AyC+gSq6yeYm9RzxF3gOdAI4HjqGd0MM8\nYk/okZun2PT2XitUyVqzVBSXW5yOti3FdipCD+8Rt2gn9OiHLjPNwN6QResDQatkTQ8JiiEuKKA+\nh7BnkILuIeUxqjAyVIYznY94NH4eot34GOPVFvfZxDVtRtiWdD5iW0znseQYqMqppuP5UkMs8dxG\nsZnORzzmflsYvJNvZiFu+56ydBW/bgkoY8RDgzTXo6CgIEGvhBXd0I/nN6hzLlf2oLzVlbyejYto\niHtZukBEzhOR74jI9SLyksz2CRH527j9S5EuFxE5TkT+VUQmReQNA7/EDYLiERcUEDwyT1zRT2g6\npa20ULT3GG3Sh2nyxm3aLeZN++QoH2oeoe7h5jxiY+aC9tmXTA+/pNC43nTw4WQLTZsX3cRVnZMH\nS2vjggwGNEYsIqPAG4FHAbuAq0Xk8mSGtmcDd6rqmSJyAfBq4ImE/L8/Au4bl4IloBjiggKHpRoG\nb9BGqTNrQfcs4zngNiqjfog6oYcvXxqlTurh4Q13N2Ytm32pSa+puPhzQ/VSsNS2KsZ3UBhYaPrB\nwPWqegOAiHwQOJ/6DG3nAxfH7x8B3iAioqqHgM+JyJmDUGSjohjigoIBoxOndCdMAze532N0NsRL\npbjst+suUx6uVQzMEJ8E/ND93gU8pGkfVV0Qkf3AccCeQSiw0VEMcUHBOkVK6GFh5IINAu0ra/p4\nEbnG/b40zrhWsAZQDHFBwYDhCT2gGgudoXPykxF67CCUL22lN4/YjxH34hFPEULfFpruROghUYfc\npA/FS15l9Jc1vUdVH9iw7SbgHu73ydSDM36fXSJi02ff0Ye2BR1QDHFBQURuXBd6CzNbzS1UhtMM\npc3qNEswgvszx28F7kbdEFsJU24+4l7qiD2zlhlOG8s+FM8xDtxMnvLyqKjXBPW5lVvUw9VNtdEp\nrG1z480FS8RgQtNXA2eJyOkEg3sB8ORkn8uBpwNfAB4PfFpVV6wEXESO7bRdVfeu1LlXA8UQF2x4\njFKv/TWDY4bMso2bsJVgiK2OeDNVna3PfLY/20xGnnnAW5PvOY/Yc03nxoibuKZNJ+/ZzsTzpYZ4\ngjrftc/ibkV9LChq6zvxV1s9NFQUnb6OeJaSxNU/BjNGHMd8LwI+QbjV71TV60TkEuAaVb2cMJXq\ne0XkekIg5QI7XkRuJDCzjovIY4GfTTKul4JrCY9IbsIaBc5Ypvw1hWKIC4YGPkTb73E+zOtLdqDu\nJeeQJk2ZAR51v02eX++RhprT8iRPxiHUdbXr9mPEqT5eRmrEO+njJ5fwzGOpPt0ICSyc7vdvZa6h\nHyz1fg8NBkjoEadFvSJZ91L3fQZ4QsOxpw1EibrM0wctcy2jGOKCgoKC9Yj+krXWJUREgKcAp6vq\nK0TkFOBuqvrlVVZtoCiGuGBosFTPyEK33gu29dA9KclC1/7YdIzYwsBNYW7bNk9dnvc0zdtNk7U8\ncslaxjXtST38uXJd+WyyzbzX3BjxLN3baJbq/uTGiJfi121YT9hj+Cku30R45B4JvIKQY/h3wINW\nU6lBoxjiggKWVys7B9xJZbRsjDiXNT1Ffix0H9UYszeQmwljqpYwBUvLmp6NcnOEHvsy+izGbaZD\nLmvarreJLSyV18/EEAU9YGNwTT9EVR8gIv8BoKp3ish4t4PWG4ohLigYAJQqWalT0lIn3BIXy5re\n7JY0a9pP+uDhJ31Is6Zn3GLlS52wGPfptl/BamFDGOL5SMGpACJyAkNYLl8McUHBGkM34+eTrHKh\naeOZLqHbYceGMMSvAz4K3EVEXkkonfrD1VVp8CiGuKBgwPBc0758qVMJFIRyoWMJdSBbqdcRb6Je\nvmR1vTmP2OqWffmShbutjngKOEAIP3fz4NPypeVyTRcMCBsgNK2q7xeRa4H/RXjkH6uq31pltQaO\nYogLCsgbz16TtUYJ4eQt8Xc6H3E6RpxjItgB3B04hmCIrZZ4wsnz8xH7siJfR2yGOJesZTMvGbPW\ngXiem8l74MdS1TKnY8Tp7EsH6W6YbWDP6pmLUV8mVDdC1vTPqOqVwLfduqer6t+soloDRzHEBese\nvu4W+p90YQvdCT06TfW3g2A8mwyxGZzxuG6edsO308kxj9jINIzQwwzxUqdBnKDKxjbdFuL5Un12\nUI1Ve48Y2g2x6ZVL+oLQyZhnb/vnCD36Seby12DXu/FC8cPvEQMvFZHHAS8i8Oa8nfC4FENcULCW\nMEq9c4fex0jNE7bQb47cIg31eghV8lQuq9k8vxGCsfCerYeFm0eT76mnDr1RXI5QvUiYjMWMfJ8E\n5mF6Wtv4ayF+2nEjtLd/TpZdtz+f97JH6e0Fyl68/EsTbESqzA1hiP8n8NvAV+Pvl6rqZauoz4qg\nGOKCdQ/rvL3H1qt3ZKFRM1rGqec94nn3PYXt50PZ3mh6D9XOk5Nlhn4x+W6/fQjazterR7yYkZee\nK4XpaTrnIg52nGVoN7X5PNWYtiH1iPsJT9txhhYbNLS9AcaICaM1Dwa+S5iM4tQ4D/JQBUCKIS5Y\n9zCPdamYJozf5kg4oPsY8f64v+3XFOaeI4zP5sZjDxDC0VanYQYqV77kPeSm8qWUazotX5qM59wf\nP1McjNdgIfWmcDKE9svJMCxEeRaa9u1i5CD99qpLnfN56DD8hviLwKtU9Z0isgV4NfDvwE+srlqD\nRTHEBQX1X5wMAAAgAElEQVQEQ9Atq7nTsfupz6rkxzChs8cIwVBNRxk7qSaR8JM+pMlanksa6mHj\nNFnLT4PoDfFeml9i9sbFwu+GpRrBMm3igLExPOKfUdUfAKjqNPBCEfmpVdZp4CiGuKCgoGBdYni5\npkXkPqr6beB4ETk+2dzLrJvrCsUQFxSsAJbiNS4At8VlLcGH3QvWEBRoDW2A/reAC4G/zGxTAvf0\n0KAY4oKCFUCuJrlblzlByEw5iorQY4KKb3o5FJfGM21lQlOEMPiddA/JW3Z1WvtbsAYwdGSPAap6\nYfz86dXW5UigGOKCAsIYbI4NC4I32MnwjFMn9EgnfYBq4gUj9EjlHQvcLX7uIIwRb6EaI84ReqSz\nM0F9fDg3+1Ju0odbaCcZGaMi9PATTkB+0oeDdPaax6gna6VGvXjcS4Bluw0xRGQz8Hzg4YQr/jfg\nLXF+5KFBMcQF6xpCRZSREk704rWNUSfOaMqaPkSeCnIMOJqQYGWzJ/ks5zRreiL+vj2Rc7RbjEyj\nadIHX2PcxDXdKWvavyi04rrUEB8Tl620Z017eXa9ozQnfhk5Saes6VlCG/d6z3KEKUvJvl73GHJD\nDLyH8J73+vj7ycB7gSesmkYrgGKIC1YN3pBYB9rviJfPKPYGNDWmTdhMFf5NvVhf72tc0am8cSqv\ndbNbl5JgmNx5KoOUyrGXgc3J7xwhRi/zEds5bb+W2745Xo+dI6fPJqdPei3pPMmmc2pIR6l79dCe\nge1foLpl4Yw6WTkWtH4z362N/HO4bkZdl1u3tz5wX1U92/3+VxH55qpps0JII1sFBUcMRrzhiSj6\nhR1vXmDKt9yLDv5YL69pft/0/IvJPl6HVoO8nJwW7e3Ravid6pvq2osMf84UubZI5fjrbSJR6adt\ne71fOR1MXr9I22PdGGFDq8dl/eIrIvLj9kNEHgJcs4r6rAiKR1ywqlhuxzdHM/NTLzB+Y6Of9F5W\nGppuOv4QFRFHKgvqodNJ8l7fZJQz4Y4zHRYYDNe0hYBtjPhQXJr0GXfHd+KatjZo4oq2tvPMY8vh\nmp5mcFzTS43ErAlsgDFi4Fzg8yLyg/j7FOA7IvINQFX1fqun2uBQDPGQQ0TeCTwauE1V77va+qwE\nBsGs1c+EAyn2EQaxmiZ9gDqzVi4x6WZCn2osVZb8lZv0IeWy9kg9RX9ub4htzmNL1kpxMO7rx4gN\nOUPcyZAtUCc7GQQKs1bE8Bvi81ZbgSOBYoiHH+8G3kBIeihYIfjxzaWwDSiwOy4pzBDnJpXwcxL7\n0LE3VJaV3G9m8lKOKTjCGHJDrKrfX20djgSKIR5yqOpnReS01dZjI8KHTnsJmwr18iUrHWryiM0I\np6Hp3HzE3iOeol6+1KtuPjRdsAawMULTGwLFEBdseIySn6jBjxF3Cl1vpbl8CdrHiPfSnt17AnAi\nwRDvJJT7+Dpi7xFb1rSQ94j92DLUPWKrIz5ExTW9m/Zyqomoy3bay8PS8qWZKCtX3mXw8xE3jREX\nA98nNkbW9IZAMcQFAIjIhQRKuSOGNOEGeu9XzAs0Q5VjfprpQZ55niYnTQCCyhDvz8jbTuXFbovr\ncuQgRuixJa7flcg5ATguLqkhtpcEz6xlc/I2EXr4kh5j1pqlMsRmGCVuTw3xCbQTejSRcByKsvbS\nHpYfIzCFpYbYe9f2gjCVOT7FGPVSM9++Jgd6T9zycyqbPuvqhWCIPWIRGQWu3AjsWsUQFwCgqpcC\nlwKIyIbjRSgoWJcYYkOsqosi0hKRo1R10Pl+awrFEBcsCbkQbq8wNqy0xAcqz6abTE8naZ6seck+\nLNvJw5qIcszz9F6n94hNx9mMvC3JAs0e8QihvcwrXkxkbE0+O3FN9+IR2/4LVCxffjx5PjmXhd9H\nk2tq8oj9vVuI+6bts9nJyXnEJscTc3Qi5TB5nnzFnkH/DHQKk5sOPtyflmYthaVrOf+JJWFjjBFP\nAt8QkU/hqghV9YWrp9LgUQzxkENELgMeQZhObBfwMlV9x+pqlWeFWgrEfZqhGEm2dTveGyivX9qZ\nNsnLrU8NZCdd/DX4/Xx5Uk6/JqS6i1s/kqyXZJ9O+nU693LbJncPOsnzczG3km3LgR9zXxcYfkP8\n93EZahRDPORQ1SethNzlvPErwftKvRoI3kgvnohPnjJvKPVqutUGW02teardkrVy8mzM1TxGqDxE\n/+cyL2uW9sSkaSfHy/IeqJGEQOVFduOabsqatvFYO5+d22B0nrNURCW+jtiStcxzteNzLPzTVBze\ndl8XaCb06EZR6QlYcoQpvdaDL0Z97bpy0ZR+ccTHlgfoEYvIecBfE5r17ar6qmT7BKEE8lzgDuCJ\nqnpj3PZ7wLMJTfBCVf1EXL9sDgNV/RsR2QKcoqrfWYqM9YBiiAtWDdahLwXWqU/RPQzZCQfj8UvN\nmj5IZbh6yZq22Zc8bBIIC1ebUeqWNe29QjuPP1+vWdO30/7yc3uUYYQeS82aXiD02lMMJmt6gdDm\ng8BQzLM8gKzpmBT1RuBRhDzCq0XkclX1nM7PBu5U1TNF5ALg1cATReRs4ALgHODuwJUici9VXWQA\nHAYi8hjgLwiPz+ki8qPAJar6i0uVuRZRDHHBhsciyzPmy30ZALgtLhNU8xFbxnIua9poLps84tx8\nxOblmr77afZAp4GblnlNqbzlsJcVZDA4j/jBwPWqegOAiHwQOB/whvh84OL4/SPAG0RE4voPquos\n8D0RuT7K+8KAOAwujvKuAlDVr4rIGcuUueawroZDCgqGHd7LtfHKdIF6DbFf/Fhr7lg/Blr+/EOA\n3id9OF5ErnGLL1U8Cfih+70rriO3j6oaa+lxPR67HMxnMqaHbmS8eMQFQ4F04vp+Ml4tM7lbHXFT\nFPAo6nXENmVhrs71ECEcnIZXTwLuRqjdNY/YZ03bHLxQjSF3ypq2c0I1N/MMVdjbPGLjmk693x1R\nF5tHOM2anqfyuA/F62mqLxmjtzrifqYwtKxxw4bktVDQ3k3SHlV94Apqs1K4TkSeDIyKyFnAC4HP\nr7JOA0cxxAVLwnLLl8y4pOObTfP+ptgaP62Db+rcJ2k2yj4M7Md10zluPaFHakCNiGMnzZM+mPGy\nJKnx+GnXuJVghD2z1laCEVwJistD8fwW+ga4kyq8Pkr1UmBt0232pQNRx5QYZAfthB5+7mg/3twt\nXC4EApWmlyY/lt9L+ZKf6SpNtJthHZQvAa3B+IY3Afdwv0+m/d3M9tklIsbTckePxy4Hvw78AeGx\nuAz4BPCKAcpfEyiGuKAvmFFoSrjphlHqhi81xAuEjvoQzV7OVkKHDHlmLY26jMXvTbXE26kza/kE\nK6gbYkuk8oZ4jGAsTUaOWcu3kY3tzsVjzIPcnll8fa8ZYT9G7A2xhzfEfox43F2HL/UyXbZTGS+v\nh70MpNdi94m4vkUwXndSv28mJzXEBjPEdv/sRSwHa2sbN/fXYfp4L7nJGNt9y9UR+4x3/7LUCU2J\nfis9aYYqLA4mFHA1cJaInE4wohcAT072uRx4OvAF4PHAp1VVReRy4AMi8hpCstZZwJcHohWgqlPA\nH4jIq8NPHVSu3ppCMcQFPcG8rlHqHpl5Da3kdxNGkyWtk1W3ramP8dSYI9T1EUInOEaddKKTnLHk\n09BK9rPFE32MJdvS/Uc6yPF6pHMMjybrc/Jzhji9J1B5yqlcf85Un9wykshM9bG2WOggB+rnMllj\ntLdLCn8eC03bPdfkOnqRkz5HOFnpvW6CUJfl298n060UVV0foelmGaoLInIRwdscBd6pqteJyCXA\nNap6OfAO4L0xGWsvwVgT9/sQIbFrAXhBzJgeCIeBiDwIeCchuIKI7AeeparXLve61xKKIS7oCb7O\nF+qdjL3999LZWJjUk2+Y8fSyOnkSfrs3EGloes7p20mOD5X6+t2U7Stl/DL5c24faO98/bjqHO3X\n5+XY2Kt5ivNOnsEMRb+h6VSuX5/qY4vxUedC03bMQoMcfx82Ufc6fRs3tUsK/+zY9Zq8lP+62/Nj\n99y/JFn7+evpBl8C5V9SrI1WkitWdWChaVT1CuCKZN1L3fcZ4AkNx74SeGVm/SA4DN4BPF9V/w1A\nRB4OvAu43wBkrxkUQ1xQUFCwTjEoQ7yGsWhGGEBVPyciQ5ebVwxxQV9YDgkHBA/BWJ3MK7H10Nvs\nN7NOh1nyXNOWmNTJM9pH5dH4cT4/7gj1ZK30WvZRefRGcOFnB7L9jIlqKh7jx633EcZWt1AfJ91G\nxTNtfM9QDxXnpkH0bFy4c89HHS1Z66A79z6nz2T8bePZdu6Ufcqu15K1rD09bJ3dJ8hnYFvbdGL2\nn4rntQQ2375pslanZ8jyBqadjJSDux9Pdrn/iSWjv6zp9YrPiMhbCYlaCjwRuEpEHgCgql9ZTeUG\nhWKIC1YFy2U1so6227R53bCfzp1/r8dbiQ60G2Koly+lyUiLwLcJRusYQta0JThZMpqVREH7+GbT\ny4wPTVtY15i1JgnG805C+VJquH4A3Ep7UpPJTw1xkyGaoz2TejkYZALUupvyMIMNYIjvHz9flqz/\nMcIj/8gjq87KoBjigoIBwFMvLiWtc5GQrmp1H2aIfflSSnHZzRDnKC6tjriXFxhj4ipYmxhg1vSa\nxUaYixiKIS5Y57ApFXPJRL30UVbK0o1r2sK5ueOt9reJaxoqj9S4plMv8Qyq2t0dcdlMZ67ppRji\nmbgcjIsRetyQ6HMCVR2xr2E2eTmu6b3k29xKu7pxTXcqWfMwcpNcHfFGmkhb2RBjxBsCxRAX9AU/\nPufH+czb6oa0htOP80HdYDWFILcSDF9OH6jPvjRJs4d6HHVCD5slyJNNmE6HCCHo3e74cQJ7gRlP\nIxkxo+4zsM0YThG83Xmqcdm7AKdQGeKdhDYywo2UWSstl/Fomn3JM2tNUU36sD1exySB6xrg6KiP\nN8R+XNcyne1+23jzdgK/ob9vJxLaeBt1ZjBrY6tl9oQed5DHDqpIQZMhNp2aJqEgXk/6gpGOf9vL\nVy8vBsv9TywZG2OMeEOgGOKCrhilItDwnU5KXjBL6ACbvJJxKirIToZ4ltBBpuQQEDrho8gbYutM\nzVDY+gXaJxwwxiczEibHG2IzZEYOAvWEoqPcknrEqXdtLxiboj5HURnio5LFPOKU4tIna3nCEA9P\nuOGTtSzJKU1uW4xtcxSVIbbrsZcUM1ipF2sesbWbEZ7c7uQc5eQYi5dvFzPE1mbE3+nLk913b4h9\nGVNq+Ex+et+NEsratpMhHot6NL0QCqGduv0nILzorMR4dPGIhwPFEBd0hQ+LemYtg68BHqfZCxhz\nxzcZYm+4xmk3xGmY1uvjPWIzMrY97ZC9HpvcZy40DRU71XiDjHEqQ2OyvKGwjn6equPuJCdtpxzF\n5WhsBIknMu9ocaFORmF6jLvr8kY5d11mpFImNai3EQRj2iQnvab0+MUuMpraKOcR4z473Xf/HOYM\nsbXXfNy3yRCnspr+E7bvoGef2gihaRG5lkDo8QFVvXO19VkpFENc0BV+EgUf8kwnru/Gseu9tbFk\nX8817UkpmnRJ9fH0luYNpmU8TXLsXCkT0mKyryf58Of2Uw/a+pSpaZFmnXJyFjPrxckUwoWKp7si\nGGPT015ivIxUj5T8g8w+/nrSa8nJN5geuW0+epG7FylSnXxbePn2XHR7flI5/iXO2s7LyyF33bn/\nBOR1WTY2QLIWoVzpmYQ5kq8hkHl8UlWHKh2gGOKCrpijyrL1Id90PKxTGQtxu/d4vUfsWawsszfn\nWU/G42z/3KQPFoq1MeJchvBkPNY859y4rp3DxlUtwcnLmKTy4KyzTSd9gPrY98FEJ5NzgGos1oeA\njeEq9YhHgJHEIzK90zHiXNa0nTNtIxtXH6fyVH2ylk+wgqqOOJVj48YjVKFxaG9jn6x1gOb7Zcfb\ni5x1XrlJH5ru+yxVmNjXZntD7MeIO3mxC1Hf9BnMjRGvCOf0BhgjVtXrCVzTfwQ8muAdL4rIu4C/\nVtW9q6rggFAMcUFP6IU6sBfYpPSG1BD38pq7j2p8NRc6tc60Wye6O8qxZC0/sUJuGsR9iYxJ4PsE\ng2NjupDPwE6nQbzNydlFGKeejbK6TYNo47X+5cPgPUyjaOxlGsRdTsZtUfZB6uP5TZN8TFNlYPs2\nn4vtczTN0yAa3aYla3W6Z3fEZQvtwxpm+Hp5Pr2RznnE/cCXrR1pKMNviAFE5H4Er/gXgL8D3g88\nHPg08KOrqNrAUAxxwapiuSG75RI8TLO8sbtpBjPn23fislawm3qG+HLgX5wGgUGOta5o6PgIYIOM\nEe8jcE6/RFXt/e9LIvKw1dNssCiGuKCgoGA9YgOEpoEnqGqtzF1ETlfV76nqL6+WUoNGMcQFqwIr\nwUnnI/ZEFJ1gD25ujNjClBaKbYJQL1/KlRxBnWs6lXcUoR65qXzJrm2Rahx0LyFk66/lPoTJXK2O\neCv1+ubcGHG/sy/5et0pKhKOmwkUm77dT426+PmI0zrtlNDjDtrpQieoyoX8fMRNhB776Tw8YaF6\nn6WcGyPu9RlKp0H0iVprHRshaxr4CPCAzLpzV0GXFUMxxAVdMUFVR+zHUv34p3XMneoujYjDk0T4\nMVmoxjJtTDZX/3uM08dqSnOTEkwRxgLvoH0c71gCe5Qn9LCO2SfcQJVktY8QrjXCiROBe1ARenTj\nmrakr71R5+vjtnsRDLFn1uqJa3oERpJB4lYLFlvthjjHNW3jujvj9X4z7n8mgenLCD18UhNU7eu5\npk3WD6lC2sfFNjqa6p6b/r6NTZaNW98eZXnsiPK2R1m5STWMxIV4fXeSr0e2MWu7T9Ce6Gf3vBsx\niDGgdfpPmD4DJ/cY4qxpEbkPcA5wlIh4z9e/8w4NiiEu6IgtVB2y/c7VTJo3Mk7oSNP+YSLKMQ80\nrSNOs6YPxW03UXlIY4TOODXETcxaZmChvUM+jopZy3vETVnTliA0T2WITUaO4rJT1rQQDM71iRwv\nazuO4nIExrfApmiJR8dgdFPnOuLF+aqTnp+HuWmYa1UUl2bgJep2XKZtjiNPcZnziE3WFJUhPj4u\nKaFHp6xpy4hPDfFxUVY3Zi0zxBPxHD6TX2i/592ypq3MLjWi9iyaPp3+ExDaO01mGwSGODR9b0KW\n9NHAY9z6g8CvrYpGK4hiiAs6YktctsXfTSxCZggX475piNJkmLxOzFq2bi4eM+lkGG9xqk9KwWi/\nF6mMyZw7JtWnaWo9qHtL9kICVWazeVemk/f2rY3sWtUd49vGL/YCsRWYGIfNW2FsHDZFSza2KXjC\nI2Mw4nt/oLUIrYXgGS/E2PrYbDDio1MwEhvBDKA/b5M+KbOWlQ+lLxqmu28fo+q0Nob8lJWeCtRK\nlfzUiSaj6aVpgfrLmE29uIXqJWybk9GN4nI00Sc1xF5Wt/+E6bOFASea6fCGplX1Y8DHROShqvqF\n1dZnpVEMcUFHiFvst58Ht+V+p2OWHiPJ9sN1sHF7K9kvJ8tv82UnfoHQ6aXnS+V4PSRZ56893X8k\n2Z47f66N0nW56xpNZAnB25WRYHDN8/VGOA1NH8aC23+0kpPq7K8/va60/dL7bcekOqfXlW7L3b/c\nvSA5JtXJ79d0jqbrSu+Df3bS8fdcEzfJyv0n0mseJIbVIxaRF6vqnwFPFpEnpdtV9YWroNaKoRji\ngo4wUgIfmrTOKh3nm6U5QcpCosZLDJWXZTKhCuHmZNlvz/1M1MF7bK2MjFSOX++vJfXy7NP2nXHb\n7bedy4dKrTP2Xnp6Xi9n1m235CgxZYHFRdgUFVqcCGHpEamMrUFb0IpjhwvxJPPzMD8LszPV+Wfo\nTR+7Js9eZt70nNs/lxznZdux0B7mtlBw2g4mG7fePOSUF9yPgafXQCLD7pXS7hGbrJlkSeH1MQKW\nToQe1kYDxRB7xMC34uc1q6rFEUIxxAUdYclXh+JvGytNE1OMQKMpWWuKQBRh4UJLlEk9Yj8+tz+R\nNUuYrH6e+hixkV1YJ6hUnbJl83ooFaHGDHVCj5wh9lnTnojjVkJHvBj19aFpGzP1oUpLSEoJPW4j\njDfaBBWzVGHdCWB8JixtWdMxWcuPEbda4bMTs5ZN6WhZ07dl9LEZoqxtfNhVnTyiPEuKuzWRMxqv\nO82aNvjxZkvWuo32zOlbqZ4xP75t93zebSfqs5f6i8EcIX/Bhjw8Z3maNW1JbfvJJ2tNRz2n6f6f\nMH0Gnaw1zIQeqvrx+Pk3tk5ERoDtqnpg1RRbIciQUXYWDAAiUh6KPmHTO0J9/Nt7SGYEcx3yKCHD\n2JcvpePpnlkrFwqGOid2yqxlLxTGrGWGeDf5cp0J2suXfATEG+IhTd5dUajqsqLV52wSvey47vsB\n3P9WrlXVBy7nfKsBEfkA8FzCY3c14e/x16r656uq2IBRPOKCnmDeTLe5V7uF3/wk8U1Z05ZdnOMK\nhmCs0mkZ02Qi82oOkp/fdgy4K9W0jGmyVpo1PUNVvmSGdDvVfMRmPHH65KZBnIr63EDlhZ4MnE1V\nvmRZvU0UlxZJ8GOkhpZb/HzKKcWlRRyM4vKbVDSXdyGUL6VZ02ntr91rb9R3Ud23Carypc20Z01D\n/dk5RLhf5v2mOI6qXCj1QP2QBlQecQ5WApXziNOs6abSJYOfWepIc01vkDris1X1gIg8Bfhn4CXA\ntUAxxAUbBycTOtK0bjctF/KTCewhdO4eJ1HV7fq5YC286Gc68mHKm6jmt71blOMNsXXwTZ37ZDz+\neqpM7mOBswiG2JfW+OvydcRW/3uAQH7x+aj3ucAphHKqnTTXEedC08cQehWAB1En9DhcR7wNNm8J\npUvjEzAWLdnoGIyOVqFpy5xuLbrQ9GJVvrQwC3OzoYRpZhqmD9XriG8mGH4zxOcC96QqW/OhaWvf\ntI74AKFu967AZ+M1nxuvyV5SfNa0b2Mbg7d7fivw31G3o+IxZxKen+00vzT5MWEzxDdRPYsnEJ4f\nXztuMixL2oy6Pc9W15zSmN6NejnVBHV97Lp8aHofdU7vZWNjMGttEpFNwGOBN6jq/DBG7IohLmjE\nyXHxdbtN5AV+jHgzVQe0n2DITyJ4WuaBbnYyoBoT9N7I9rjdDPFJVB5oqk+nMeKJ+N0M8d0JRBx3\nod4p9zJGPBp1GImfJ1HVk/Yamt5K3UM6OepzInDsCOw8FrbuhC07YMv2aIRd+dImf8FNLvF8SNKC\nkKi1MBeM8fQkTB+EqQNwYC9sbQU9ffThHq6dew1NWzu24vH2eSJ5Zi1DOka8I+5nDGQnuzbyDGa5\nMeKUQMMiFLc4GXcnPI/pGLF/Du0l41Cizz6qF4OTCYbYiFxyY8SpPvbMDtIYbwCP+K3AjcDXgM+K\nyKmEv/VQoRjigoKCgnWKYfeIVfV1wOvcqu+LyE+vlj4rhWKICxpxLMELaQpNp6QMxqzVomJq2u/k\n2Pinz5rOhaZnqbwwy7jFyfAecTeuafPC/Hjh8U5OJ65pT3FpTF+z8biR5LqyzFrx32W0k3NU4VlP\neGK6HDsCx54YPOLtR8PmbYHMY8QzjkA9JJEWY1tIYKEqd9oU4/StaZjZATM7YXJbDHXvhtlWOL8h\nbedx6nSaiwt1b897u4fica1MG3uKyxzXtIWKFwj3yNrGdLJQcC9c03Yuy2o2WSmzVi40bR6xPeuL\n8dh9VM+1ycoxa/nr8vqMELz+QXnEOsQUlwYRmQAeB5xG3V5dsioKrRCKIS5oxE5C6M0oF6GZRWiB\nKuQ46/Y3OTudHB8KztURW6e1mJx7ewd9cnXE1hnPxGPIyLHxyyZ+aKIsqzudjMcI1VjuTiLbVzRU\nRkc5OhbGcbUVOsyZyGy1kGmfHQQDvPNY2HEsbNsJ4xbvNotgVnycEKuWscDs4dFaAF2Isei4bjY0\n7Mg0bJ0IzFxW8rQwC5N72tvncNtEZq/0WubnYTxamdFW1U474qJJG3tj3cRbbW087drH9Nru9LJ7\nnqsjtmdwhKrMzcuye28vBk11xPZSZgxdXg8vy67Nnh3/n0j1WaTezsvGcNcRGz5GeG+9lhWg614r\nKIZ4A0BEzgP+mtBXvF1VX9Vxf0InZZ7rZuqT3ucyRH2pTEqZaAlVnsLRvJGU4tI6NEu+8ee2Y738\nnCH29ILz7hicnByNY+rVmCH2lJv+3J7ecithTBfCmO6miSqhypKnZAT0AGxZaNBnO0xsiePCZoRr\nfIxx1HNiB4xtjtRaGUPcWoSFGZiN5I7jUxWPZBTVaoXkrS3bYcueBn3GwvbN29qvZWy28vh1f2hz\n3x65dvZZ32k0xe5dro2hft9z0Qs/Jm9y02dlM7Td8xzFpWWl2zNgmeteHy/H5xf4/4TXZ8Gdf4LB\nZFEPMmu6Wx8RPdP3EHLw7gCeqKo3xm2/BzybcNkvVNVPdJIpIhcBv0nICTxBVfd0UO1kVT1vMFe5\ndlEM8ZBDREaBNwKPIkTFrhaRy1X1m43HUBkmK8+wDssbrDS7eJGqcxt38iaS9fY9xzVtHasd4yOy\nduxYZl3aCYvTJ520wMuxsiBbN0JlcNLErTG3jCTXtWlzNSmDecTeRspcWL8wD5sO1ttnDBgfD07u\n+BbYNE691x4HJnYGAwwhpjy6GTZtyRvi+em6kR4ZAw6EBok3bNNMONfYeDy3swyH71XM2B7bFPYz\njC5UbYRd+0z9fhnft2/jGtd00sbq9rV7Ae3PnW1PJ31IIyteTlp6Z89zJ65p7xnbMVDX6/C9p/k/\nYXbSzm2LvbguF4MYI+6xj3g2cKeqnikiFwCvBp4oImcDFxBmSro7cKWI3Cse0yTz34F/BK7qQb3P\ni8j/UNVvLPtC1zCKIR5+PBi43ibXFpEPAudTzXrXBuP9tcUPQ9pvP7aWrvfbvTxbny5QdaBGMynJ\nJ8kxklnnE4hbHbal56/p7UKw5hKPdTjGlrFNMGqzI0UbOBZ/t1qx5Ch6xWkIfJRw7OimcOzoYaHE\nN3HYc6wAACAASURBVJbxYHQP1y85I2wLBCMMYZtHawEWZ8JgcZxVw85z+LzOENupZSTqnfBaLxDC\n0Xa9Y5tgdKb9/tTax11z2sZNzwbUn7umZ8c/C7nnwu+X083LgIqu0uvl5UL7s+5ljToZI8n+6fO7\nLAyufKmXPuJ84OL4/SPAG0RE4voPquos8D0RuT7Ko0mmqv5HXNeLbg8HniEi3yMEyQRQVb3fEq91\nTaIY4uHHSYRpYg27gIekO4nIhcCFsDLk9AUFBYNFn6Hp40XE8zZfqqqXxu+99BGH91HVBRHZT8hX\nOwn4YnLsSfF7136nB/z8Eo5ZdyiGuACA+Ke8FGBURI2r2KKZ9n+3sS8/CYDfz38elp2styX1QnDb\nNPkkOTbdPycjvQYatvsxPVzY1eQuZI7R5PfCfJgDGGBxHEYW6uxQiwth0VYlr6ZPPH5xMSyjFiOd\nJ8SN56fbw9Dm+bbcmSw0vTjjZn2YDslbRrG1WJ3n8HkTfRaoErMW4/hwa7E6hR2LXXtDmy46eYej\nAK6N02fGL1B/7vySRj+anov0+UjvncmRhuNatMv15/TnhvrEGP65S69vUIwUi70L2rMeKS5V9fsi\n8nDgLFV9l4gYr8tQoRji4cdNBG4Fw8m0EwXV4En9jaHKxsbSELVP1rLMVU9/CBX1pc9sFeqzHfms\n6Xl3jKfNNH08mYR1oDkKxgV3jB+T83LM1tVCj616sta823fBrffXNT8D8zFyLNNx/NMb9cWQHDU3\n3d4+C8DcXEx0nob5uRDqrcf5D1QGd2IHjC0GY9tLstb8VNWY08BMOMdctM9zc/WXhsP3arqa+1hb\n9WuZn3WEITN1ikm7Nt9evv3HAIkWatGdz/b1E0p4UhV/7/10hT5r2p97LlnnnylPvgHthtg/Nzmq\nSv+sp8MzvpzK9Ldz2zIIQ2wVAgNAL32E7bNLRMYIVWB3dDm2r34nBxF5GfBA4N7Auwh/9fcBD+tX\n1lpGMcTDj6uBs0TkdMIf4QLgyZ0OsJKN2GczQ/WgeO8xLV+aThaD5zi2DrBF52kQTcaMk2WTFWyh\nnqhj589Ng2jHeH28juPUr8OP3frypfTaJNFxFJBo9xZj9VCufGl6ob19DsudDJnX05PB+B3Oj7LG\nmojMx/NTSypfMpLpuclwjtnpyLTVpM8CyGR8qciUL83Fg6ao7pO1h1JvY6i8xqbypdzzY59ThCzo\nTVT3PFe+NOuOm0pk2TmMEtXODflpEO1YT53p9ZqmSvYyjz8tXzJ97H80zWBrcAbkWffSR1wOPB34\nAvB44NOqqiJyOfABEXkNIVnrLODLhKboq99pwC8BPwZ8BUBVbxaRHUuQs6ZRDPGQI47nXAR8gtBX\nvFNVr+vl2AMEyspRqjfvboQek/G4yUSOUU1aNrNlITcRehzKyDJ+ZOuMTZ9OhB6TTg6JnC1OH/PU\n084UKkIPa48D8VwHo6zDhiYqtflQM6FHU/scJFBOWj6WtmDzPGyebSL0mAup2B0IPWrupBF6TMHM\nIZjcF853YG91Tb59DhDuTyt6zJ0IPYxK1Nr1YFTD5PoMdjJtbBEQ38bWPqbXJHXClG6EHqaL3X+T\ntYPqee5G6GEc2raQyLJMaqt999Edu65Un0FzMw7Es27oI0TkEuAaVb0ceAfw3piMtZdgWIn7fYiQ\n2LUAvEBVF+FwmVJbvyMiLwReTKDs/rqIXKGqz2lQby4afI3HbmvYb12jGOINAFW9Arii3+P2UvEH\nm3PVC9f0XuozHtnvUULn1CvX9L54rHVmd1DVcab6dJuP2DNr7SGwNI0ROswZeuea3huXkfi5ze3T\nxjW9kOeazrXPXmCiBewOQ7szh44M1/TeVnV+wx1UnMqHWc5aMNKqRyA817TNR2yyWnZNcX/vGXeb\nj3gv4R5Z25hO1s69ck1bO087WZup7kUvXNM2e1eqh7GrWd1zL1zTd9I8I9RSMSg+j1wfoaovdd9n\ngCc0HPtK4JW9yIzrU9rKTviQiLwVOFpEfg14FvD2Ho9dNyiGuKCgoGAdYoBjxGsWqvoXIvIowjv1\nvYGXquqnVlmtgaMY4oJGGCeunzmm12kQPZfyPsIg0RxLmwbRYN+nqM++1Ms0iJ7f16b9myN4PP1O\ng7gr6r2Lij+432kQ07qO7XGfqRbs2APb9xzZaRB9+/yQ4C3auGw/0yDuisdr/FyM2/qdBvFm1zbW\nprNUIeF+pkH07Ww0pf1Og7gv7rffyZphladBZPgNsYi8WlV/F/hUZt3QQFSHbmrHgmUiN99nyk6U\nyxD1GbNN2Ep9CjrfeUEV8pyiPo7q4Sd9yFFcWqh8hiq0mGKMMHeu8Q43de5QdfD7gN1UyTbbqaYL\nNM5qnD65rN6pqM8NwG1x/5OBswkDZn6SBOvgzbB7isgROkemHR/J4YxknzR3iCoMfAthgM+MxF2A\nMwhFoofZNWkf17V7bS8pe6MMu28ThGkQj6Z6YYLmlyYLBd/qdPc4jnC/rD18KDhN1jJDnIPxnfvQ\ntNfHD49MNcgweDaxpv8E5P8Xqrqskv0zRfQ1Pe57Ply7HsuXROQrqvqAZN3XC6FHwYbEXPK5VEzR\nvXPrhnRMcylYYAm1FAkmgW/H72P0Ph9xmjW7i2DgT6Ru1D0nsjfEfng47cktXGnjplAv5bJsX288\nd1OvI74tLhP0Ph9xznDOEiaSHRTuIP9S1S98Atdy4cujVgPD6kaJyPOA5wNniMjX3aYdBIrMoUIx\nxAVdMUE9NJ1LTPHJWk0d01bap0H0Cb9Q90b20260d9A+DaIZqU7JWgcTOccSPD+b0L7XZK3dhFAl\nBMN5D6rp8LwhtmQi70X60PQNwPVx/3OAH6HdEHtP3fMiH6ZJjKFpm0lJWy40TWVcfY2tvQykhvhb\ngKXSn0nwiI+lPTStrn1zyVo/jPIATojXZPe812St22h/0doR29mmm8wla9kzCJVHnN53ewZ9VAbq\nz6HPms49g4bxqFe3/4TpsxJTBw1xaPoDwD8Dfwq8xK0/qKqDznlbdRRDXNARO4BjyE+DmHpINver\nz1Q1bCUYPj/tYKc64ql4nt1UHf4EIZx8LJXRs045Nw3iLFUn6TtkibrclbohTjmgvUc5TZVxa4b4\nrlGOzZObhqbNazUP1TJxNxE6ZjPEd3HLMVSGwma+mtgcMqZtUomxicj/LJURNmgLWnGe2sPEWvMh\nc3p2phrXPBTb0zLH76AyxNY2xzg90mtJQ9PbYtvNURli38Z+9qXUEJvhm6IqUbqTurd3VypD3Gka\nRBsj3hbP48PD44SXg2Ooh9whX0c8Fc91G+3GeIvTx0+D2BSa3havKX0xWC6G1SNW1f2E96AnxUkp\n7kpo5u0isl1Vf7CqCg4YxRAXdISNgdn4nn1PO+bZ+Gljgakh3pxZPFG+vdmb57XoZPlSJTt/WlNq\n3rUne0iPmXXr/Ho/K44dZ8xNOLn+vLjf48k2qy7yHbMZ9UXq7Qn1a/K6mRHevDUYYpsByWZ28hMx\nGFqtmLS1APPx3z06F5K7AJipPPRZ6u3QpE9aHmbtm47zbk7keNm+DNonftlLk8n1552m/tz5dkkN\nn08ehPr0hbMZGbbkxoj9M2nXlRpiL8tmlkr/E37YwGqNBwlf6z6siLXIFxNSB6ybUKCMERdsHKhb\n7Ld1Mi33W5PvKdJ9F6lCt7hjOsny2/z+PrEqPU8r2d8f47mH084zd460XMQf53Ww9T5UntMpva7F\nRJYSPFyNxtVm2mn5TKwEZoQtRG3rTE6qs7/+9Lqa7m3aFqnO6XWl23q9FyTH9HN+pfN1pffBy1ok\nfx6PnKym/wTJ90FiiEPTht8E7q2qg0gPWLMohrigIzylI1QdlQ/helYjy8jNyTlE5bku0pw1beOY\nqSxLMjKGrlQf79VYONmOSbmvU33MtqVeje0/RXsWradQTMOtufHmlCrTX5dfzPsEaM2FMPP4lkDK\nAbF8aVPF3OXHiCHsvzhflS8ZHeVcq6JanKL9vDl9LDpRm8CCdmYt39a+fayNoWrPXJjb3xNPA+nD\n31upyoz8JBI+FGzXYffKX5fJ9pSbFqnIZU17fVIcSvTp9J+A9jYeBFbKuK8x/JB6NeRQohjigo6Y\nJtQF+7rdHLOWdcwHac6g3Us1lmzh5LST9oky+6h3NAuEsUwLhUM19uhZuqwztRKo3Kv0HYTO1/TJ\nhTs9xeVU1GdPImMrVfJSUx2x6W4GPWXWsmzgTU7WDFVS23gLNh1qT9bytJMGo9NMk7X8GKrVe3sW\ntFSfo6IuVv/bxA8NzYxqe+I1+Xtu+qf5BT5Zq+l+EffbSnudtj2DPjnqTurPokY5xvRl9wkqQ+wn\nnrCEtlySlT2LFrru9J8wfUqy1pJwA3CViPwTrglVe67cWhcohrigK2YZTCfiS5dGyWdNp95XCuMy\ntge3G9d0k957CR21JUZ5Qg+fAAT1rGkvb3e8Hp/NC+0Z2GZwPMXl952c/4q6H6CeNW1JZI1Z061A\nO+nD6T7Mm8uaNoPns6ZvjjoYro/HWFKcL1+CvEds2enedbmDYIByWdNNXNP7yXt5dt/9uHwnj7jp\nGfKeqTfmadZ0L2OvcwymnGo52ACG+Adxsb/BUKIQehS0IUfoUVBQMFgsl9DjdBG9pMd9n7ZOCT02\nCopHXFBQULBOMaxZ0yLyWlX9TRH5OJkgiar+4iqotWIohrhgVZEma/Xrilu8KseF3EtyzBbqdc3p\nvMtQhTz3tR0djj+WitgB8mFuT1RiJBoe96ZOcWnMWjmKSwvtpiU7Xm8fWs9RXPqZjm4BvpPIMXIR\nT3yRCydDfYw41+ZHUw9Ne739WGqaXNUEX7eL02kpLFf2zPiw+3rBkCdrvTd+/sWqanGEUAxxQU/o\nhWt6ms7ju9Ab17SfQCKHo6kza6VjxJ7VqClZa4xA7pCOEXfKmj6KOsFIyjWdM8QGT1Syg9Bmnmv6\nVCpDvJP6ZAJNY8SduKY7MWv5eaN3xHMcos41fSrtXNOGnCE2WZ5repyKWctP+tCUNW1jxLfTzDVt\n7ZIm+uXmI869OEG9bZfLNT1G+zPo9enENT0IDOsYsapeGz8/s9q6HAkUQ1zQFeN0Z9aaI3TWB2g2\nxluoGKiaCD0WnHzf2Ru2Uxkqr09u0gd7uG0WppwcS4zynjXUyRLspaFF3bBvT/TpZdIHm8N2O5Uh\nTuXYRBSpR2zX48lCssxa1Etz/CQYXi9LTPKza5k+O6hPjpBLsPIeLu66rK13xMVeUnKEHvbsWNsr\n+eiDtVG3SR/shcEy4tP7bnStnSZ9MDpSk93kpY9RZ/pq+k8Q9VgJYzzEHvGGQjHEBV3hO85R6gZ0\nNNlnlGZDbNstW9rL8YxWI27fJl1SfYx32GBy/Dma5Pgwr6dNTPdNPVB/bu+Rjya/03WpTjk5o5n1\n3miM0FBHHC2Quv3NeHqZI4nsnD5eD99WUBGW0EGOb8/0mfG81bl7kSLVKY2meD3pIsfupW9TnzVl\nbTearG+S1e0/4dcNEk1kIwXrD8UQF3SFeSz2Heqdqb39+xl5crDt5iVYBwqV4TPPpknWnNue08fk\n+nBsJzmmj/eI/HVBNb1jKmsuWUyGHZ/ziOedvCY5s9THvdMSr8OGaCE/+5KFV01/Xxvrz5Ge1+sz\n63SG5pKjnP656/IdjW+XlmuXXu+Xec7+5c0/g37/nBxrC//CleYXpM9ZDulz2vSfgLwug0DxiIcD\nxRAXdMUiVX1oL3OvNmEuyrHwc5pwA/Wx1FwnaGxPtn+qj+lr5A6T5EOLB6leLuapalO9R+zn87Xx\nS18nu58QKjV2JQtNz9I+3pySRKRy9hPCt2ZkrE2NI9sna6XesoevIfZjxD5ZyxivTI/cdR2I57QQ\ntGcg8+O6RHkHG+RYmyw4/X1I2QyxJ/TITYxgddw+5JurIzadDpC/7wtRjl1XOulDOkbcyRBrlLWc\n/8RysZ6Sy/pBU7a0oWRNF2xoLJfcwzrC5cATgwjNhrhb4pgRUDQRevis6VzizhyBaWCSdkIPP6YL\n9ReMvdTHQG8jjFvORVk2tmqTH9j4tcmzFxgzRClftxlhu37v4RnFpRFkWNa0jVcTdbPryiVr+TFn\nqAg9jDnNYzdVG6djxDlCj17JXJoMcS8e4lKyqzthUIQ3/WLIQ9MbIlvaUAxxwbqGsrovBgsEI3Zb\ntx274Ia4GHJZ04MwxN0y0g23U033uFz4F6dBwF4ECoY3NL1RsqUNxRAXFAwAVsYCzVzT5lnnXhxG\nCeVLx1AvX/LzLedC09BsiLtxTR8g0HzeQj7EOUG91MyQ45ouhnF1MMQeMQAichbwp8DZuFlIVfWM\nVVNqBVAMccGqQKgSgKDdiHSDGSFfw5mWxKRJUTn4Gtd0DmGfaGNjlKm8o6iMZ5MhtvFmP+mDn9V8\nFLgPcHcqchCrtzaiEV/z6jOI02QuC1f6drRzG3GGhdktNH008G3q7X5K1MVP+pCOo3pDbEY9jS6M\n0841nc4dbaHpqczxKWzu51ypma8j7vUZyr3MzLM+PM0hJ/QwvAt4GfBXwE8Dz6S9fH7doxjigr7g\nx+dSw9eLVyRUnp51yDbVnsmy8cKmTsZqQaH7pA+TNIdFj6Y3Qg/zZIV6CFqijNQQp4QefgxzPP7e\nSxUeNhlHx6VXQg8zwv1M+jBLNaG9JZlNx/P6+mjTZTvNY/Dm2VtNuBIMsr9vvo3NEDcReozHa2ki\n4tga9ek06YORhnSa7UioT2aB099ewCxy0IuhGyP/YtBLnsJysQEiEVtU9V9ERFT1+8DFInIt8NLV\nVmyQKIa4oKCgYB1ig3jEsyIyAvy3iFwE3ESdf2YoUAxxQU8wr2sT9VCwr+c0juNOsON9SQxU3sQC\nVXiwyaPx86GZh+5LYnw98aa2o9vleI/Tj+saM1XL7TNKnXHLrse8IrvG3NzGmpyPjJxN7vixZH3K\nrJWSmJiu3iPzsPUm18tO9RlPtvlrSYcT5pN97b6NUr/XXn+Dr9v2bF05pM+OfwYXqO6TXUO358e3\nqQ/vpyVnnWBZ+00Z9/a5UgZzAxji3yAEQ14IvAJ4JPD0VdVoBVAMcUFPyI3h+k6nFyPsj/cdMNTn\nI+42Tuy3m5Ezo2SGaIHe5RgdZlqTu5jsl8rzxBmLyf7+ZSCt7c3JablPP8Zr330IOiUO8ciFplNZ\n6Xn8BBe56/VtnatT9u3Yixz/8rCQHNvr/bIXIv+i46+jFzn+PvnnsJdn0GDtkd6Tfv4Ty8GwJ2up\n6tXx6yRhfHgoUQxxQV/w7Ebp238vWCT8o8z4paxRNkbcSaYf8zXPr2mM+FAHOZNU46SmT8pfDFUS\nUFrysxDlb47n9iQjuTFiI62YTGRNusVIRbzRG/QYsSVrTSZLqo+NZ+faN530YZJ8HbC1sV0H5HnB\nfds0wcbpbcam3BixJWt1KpcyDmo/RuyfQ9OnVyNqdcRL/U8sFYOsIxaRdwKPBm5T1fv2eey5wLsJ\naQ1XAL+hqioiFwO/RlUJ9/uqekWfsu8F/A5hHpLD9kpVH9mPnLWOYogLloTldDSdiPR7xVTymZbX\n9NKJzhKSr2xihVwCEHSeWep2KiKObXGdhX3TDGxL+tpLvf2mCCVEptMUnadB9FzdTdMg+qzvXqZB\nTF9ubon6+vKlNNPZXlQOUTFrpTgYz9c0DaJ/aepWD65UpB4pl3S/SUuLDLa22WQeaQzQ43438Abg\nPUs49s0Eg/slgiE+D/jnuO2vVHU55BwfBt4CvI3VaeIjgmKIC4YCy8keXS4zUkrtuBTcFJctVDMf\nWXa5edhp+VInQ+w9YvMgZ6nXEZuRzMEM3iCwMEBZBk9BupExKMukqp8VkdP8OhG5J/BGwoyhU8Cv\nqeq3k31OBHaq6hfj7/cAj6UyxMvFgqq+eUCy1iyGrh6roGA9o+U+LfSYLlBlzKbbfCZt7lgfzhz2\n8cVhh/axAMeLyDVuubCHU1wK/Lqqngu8CHhTZp+TqKayJn4/yf2+SES+LiLvFJFj+rg8w8dF5Pki\ncqKIHGvLEuSsaRSPuGDDY5TuXNOdQulb6cw17ccwjdAjlXcXArPWsVGWzUfciWva3qLTcXY7XxPX\n9DTVpA85rmniuY3QI60jbuKa7hTqzc1HnHJND23ccQXRx8vUHlV9YK87i8h24CeAD4scTmOcaD4i\nizcTMp01fv4l8Kw+ZViG9O+4dQoUZq2CgkGgafalXsa9rGtoIvSw2Ze6hS930PukDzZzVHq8Z8OC\niomqadKHUerMWkKI/Z3gZNnk9Z0IPfqZfckna01RjT9bItXt1NvddOk06YOffWmCYNTTEPQY7YQe\nnSZ96BbCNirRJsIUe8Hp9Rlqmn1pvWAFs7JHgH2q+qN+pYiMAtfGn5cTjO3JbpeTCaMsqOqt7ri3\nAf/YrxKqenq/x6xHFENcsCQsJ0NUqDNrQX2u3bkeZG5xnzZ+6rONrZPu1LFPUHmdTTMmmb42XV8q\nb7OTYTqZ4fQlVZaNvUjdAJIcv4Vg/Oy3N8R2Tf3WEftJIjyj1UJybvNoLUJgy+bk2HHq8y+brM2Z\n9rFrsusw/Q2+dMzauNN4vY8S+Pa1kiZDt0Qsq3HO1RH3mzXtZa7XrOk22aoHROR7IvIEVf2wBLf4\nfqr6NSA1zgdE5McJyVpPA14f15+oqrvjbr8E/Ge/eojIJuB5wE/FVVcBb1XVtEx+XaOMEQ8xROQJ\nInKdiLREpOew1JGAGYzlQt1nOvbZi/xkHO0wch1ck7zc+vT4Trr4a/D7pWPC3fTrdu5Wss2frxf9\nOp17uW3TK0tUWnfun6NBPU/raew8V6vdqWa9CSJyGfAF4N4isktEng08BXi2iHwNuA44v+Hw5wNv\nB64HvkuVqPVnIvINEfk6gSf6//Z/hbwZOJcwPv2m+H3okreKRzzc+E/gl4G3DlrwcsuXZql7Mdb5\n9RoW9CHIRdrLa/wYZhOslMd7MZ3Kl3LypuPivXvTJy1fmnX7L2ZkWMg4nQzDMoTtz7pAVb7TlDXt\n66BzWdNTyacfs15068xTzbWvsU5NZ2QYZtx6u+am8qVpumev2z3Ite8c3e+5wa4x9Yh7NVxNMo80\nBhWaVtUnNWw6r4djrwHaao9V9anL1Qt4kKre3/3+dHwxGCoUQzzEUNVvAbhki4KCgiHBevPel4hF\nEbmnqn4XQETOYAjz+oohLlg1LMf7ME+gyRPrFZPx+KVmTU9SjW0ai1enrOlJQlJTitvj/urO1y1r\nWshPg2i695o1bbRHqT6L5Gdf6idreoEwq9OgsqaNFWtQWE+JWTlsAEP8O8C/isgNhMfmVIaQ6rIY\n4nUOEbmSUPmS4g9U9WN9yLkQ6KW2cOiwXJalqYbjfZi7WwjxdmAP/WVNL5Xi8v+3d/6xllXVHf+s\nefODgWGgBowoVos/MJTwo2CLNsaAJKUNihJpa2ytlWCo1tKmDakhaRObJk2MTZP6o5lYgiQElQJC\nilQhktI/gBYaaEFAoaR1CC3gMMLMY97MvFn9Y5/lWXe/fe6vd9479563PsnNu+fec/ZZ97z37vfs\nvdf+LrOStHrEpdiWgOfcthf8aW6eVnvDFJTpe9GHqgTi24BTq5eeVNXV+O/MJCHEc46qXthSO7tI\nC/gRkb7/f68LkwqWknqPP85ez3vEPgO61CP27low2COelvXKBA7GZ9wEt3lERC5Q1e+JyKXZW28V\nEVT1lk4CWyNCiIOgBRaoly/lRR9gcB1xSRCF2tDj2OqRL1/Ke8SbKCdr5UKc94hfpbawNEOP0hf6\nVsrriEtFH0Kku6HH1/29wPeA9xfeUyCEOJgPRORDpDV9JwJ3iMjDqvorHYfVOnkG8aRzzyZ4oww9\nmuYTjye5YTUJcWmOOF9v+3rgJGpnLfOaLgmxL/hQGpr2hR+gLMT7SPO626s2ns3iMaOS0hxxSYhf\nBvY2XJ/NjGfoMcnQtd2MwGCm+Eaiz8laqvrn1dPPqeoz/j0R6Z3JRwhxj1HVW4Fbu45jGPnypUm/\nTL3ZRe6ONM4XuzfQ8IUV/PIle61U2GE7SWSsWhKURf0wtZCWjEF2VO0cQz0/3FR9yZY3jbt86RDp\nGh1g5XVaqs6Xs8PFZMPiuXja786MSprE1K5PqfpSXtJw3N9Z7spmNzqTDr/nNZ7nrYe5AW4+bgZ+\nIXvtH0jriXtDCHHQGV4Upv0itMxhL3rqXhvVXl7JSNzrdnOwye2Tt+cF0de29fO3R1x8JTcsGEy8\nEtfGpmw7P2e+MM1fR+9FXWozT/bK28mvrX0Wa8N/3pLdJjRfW38e/zsYhY/F2rSe4TSL9Pw1tfbn\nSYz72iMWkXcAPw8cl80Te1v33hBCHHRGG1941gPywmM9tnHat14ipB6b/3L3Q9MHGtqz9+xYw9ta\nWu/0oPtZamfJtbeFQQH1w5BmPDLKa9ovX7KELX8Ov12Kx8dqvXDr3fvlS769nGUGrzHUommfyeIY\nx4zDer7extQP/U/KPIluiR73iE8FLibN/Ph54ldItY97RQhxMNeYS9e0HCYNOdvwq3dr8kLcND98\nmDQ3ukw9NG1DyXnhels69FKhnb3VcZuoRdSvI95CuejDOPWID1GLal6PeC/lud2XSGK5RF3Ewsid\nyyz5q+kaWYJaPjRdcukaB7sR2Oj0fI74NuA2EXmXqt7XdTxrTQhxELC65T0HWbnkqGTBOKz3tYeU\nELaHVK0onyP289fWWx6WNX2ElRaXB6jncRdJNyAvUb6ROUxdGnGBssXlJIR4rg3z3qMfgytF5HFV\n3QtQ1TT+gqpOWk5xpgkhDoI1YBrhWSItJfrf9sNZFatxQAvWjj6vI3acYSIMoKovicjZXQa0FkT1\npSAIgjnlyJiPOWZT1QsGQEReQw87kL37QEEwC+RZwIcY3XvZzOA6YquTPMzQA/fTlwZsMvQ4UD1s\nHfEeRvfchZWGHtFDng02QI/4C8B9InIT6U/xw8BfdhtS+4QQBwHpP3xY1vSoY5sMPfKs6UXKo5JI\n0QAAEcJJREFURR+OpXbWyoU4LyLh55+HOWvlWdMHWCnEx5GGwvN1zVSxeGetkgkH1IYeo0ShZOhh\nWdMbQFBap8/JWoaqXi8iD5HqGQNcqqrf7zKmtSCEOJh71spZy5KdrBpSU8/xOMZz1lqqXjvESuHb\nWbVjbR1NMsJo01nLEr/2u89pBS/yeMxm02IY5qy1tWprmLOWr77UJOrhrDU5G2FkQlUfE5EXqNYP\ni8jPqur/dBxWq4QQB3PPAoNf7jBexSM71lyr8iVHhl87m2PDtpsoF2XwJhjLDC5D8liPdyF7vsBg\nDxhqkS+ZcZgYb6K+kbA2llnZvu9peyxOuzb+s0C9Tprq9fz6l9qyz+3P512+xjXTsBsvf9MEGzMr\nu+83HyLyAdLw9OtJifxvAh4nmX30hhDiYO6xL2/fYxv3C8qGcE20zJUrt4g8tPLQn+5/iHqJEgyK\nZj5U3NSWCf1y9ty2/TC0nW/cHvFyob38XDkWp7egzEcc7DhbKtV0zQ9Rm5sYeY941PIuT14J6ggb\no2eYsxGGpoG/AM4D7lbVs0XkfOC3Oo6pdUKIg7nHeqzT8ipp6La09hdGzxHbsK4JbD40bcJxgPIw\nMKQ51h3U/5Amdttce031iHPRL9UjtuFfMxXZX53zlepn6TNtqdo4yp3PYsvniEufyThcvW9D06U1\nyZMKaSSMJfreIwYOqeqPRWSTiGxS1XtE5G+6DqptQoiDgNV9sS8z6FCVDyf7OeImXgF+RBJImyPe\nTj236g09bDi55FttQux7utartbluM/SwrOnFhpj2VA8/R7wa8TyY/QxWzwboEe8VkR3AvcANIvI8\n6d+kV4QQB0HLTCvqizSLoscEeJjX9Lhz5OOwGgvRYO3YIIYel5DuHf8I+Cgpn/FznUa0BoQQB8GM\nYdnKTcuXprG4zJcv7Wf4cHIwH/R5eF5EFoB/VNXzSX/SX+s4pDUjhDgIWkBoLvqQzxE3fXnaOuJj\nSfPFVvRhC4OFF/LSgj5ZyxJ4fK/cCk6Y37Qv+rCHZkvNBepazaV6z77owwbomc0cfU/WUtVlETki\nIsepaqkceG8IIQ4CBqsvGT5Za1jPYytJOI922yVzkAMkUdtTaO94kgibEB9DLYJN1ZeGLV8aVn1p\nkcF1yQdYuQZ4gdrQwwuxN+Gwud5FkrAPm/tdYDBZy5im+lJQswFugPYB/ykid+HmhlX1D7oLqX1C\niIPe4HuGk7CZchlEE2Vh+NztVmrjDSj3iH3i1LZCe15wfXv2WqlHbHPF+fKlPGs6r5NsWc92rlKV\n9W3uPUsU8z1i36bVAh4mpv765IlfeZb6uEz7++4Tfe4RV9xSPXpNCHEQBMEc0udkLXPPUtXezgt7\novpS0Bum/WLyPclSr3JUQoztn1e8WWZle03mEz7T2e9baltZGWv+vHRMHpcfxm6Kp+kzlNoehr+W\nTXFPSp+FaFx6XH3pW/ZERG7uMpD1IHrEwYbH/Jab5ohHCfEig/aPh1k5R2wJU4uUlwPtr96z+Vio\nRcrmUH3W9DhzxHnWtF9HbOczc48ci9WOt6HppjniUcuulqivZ2mOOJic3GFsNYjItcDFwPOqevqE\nx54DXEfKL/w2cJWqavXeZ4BPV6HeoapXj9use37KJPHMIyHEQVDRZPc4DvuqB6xM/PLJWk09uEVS\n9vJBkjDmyVrWJkyXNe3LIJr4WtZ0k4j+hGT60ZQ1PWmC1XL2M1g9LfZ2rwO+CFw/xbFfAa4AHiAJ\n8UXAnZUd5SXAmaq6JCKvnaBNbXjeS0KIg6BlRiUuNeHFHNI/5zAh9v7PhhfhJiGe5GZDmawqUrC+\ntKVQqnqviLzZvyYibwG+BJxIul+7QlWfyPY5CdipqvdX29cDHwTuBH4P+CtVXarO8fwEIZ0pIi9T\nrQysnlNtq6runOwTzjYxRxwEc0o+N9umm1YwH+iYjynZBXxGVc8B/gT4cmGfNwC73fbu6jWAtwPv\nEZEHROSfReSd455YVRdUdaeqHquqm6vntt0rEYboEQdB63ivaWOcodztlA092u4R54Yeo3q8/tww\nvdd00C4TGnqcICIPuu1dqrqraefK3/ndwE0iP52u3TZhiJtJf9LnAe8Evikip9j8cVATQhwEDm9c\nMYnQWIEGqKslmUhaW0vUSVI5W4HXkr61dlILcb6eGGqf6XHLIJqZx8EqhldJc9DbSV8Az1G+STi6\n2s+KPth5zD7TV1+aZPh62mscrGSC6/eiqp47QdObgL2qepZ/sbKdfKjavJ00P3yy2+Vk4Nnq+W7g\nlkp4/1VEjgAnAC9MEMeGIIQ4CEj/CJaNnIvEsCQrO7YkxNYjzo0rSpaQ27PHNpqF2ItwkxD7ZUnW\nO8+NNA658+VCLAzGktcj3sxgb9xqCjch1MYh+TUedWxQZi2Xb6nqyyLyjIhcpqo3SeoWn6GqjwC5\nOL8sIueRkrU+Bvxt9da3gPOBe0Tk7aQ/4xfXKOS5JuaIg2AEbf2T5EuN1uOcq217kpjbOmcwPm2t\nIxaRG4H7gFNFZLeIXE6qdnS5iDwCPEbKgC7xKeCrwFPA06RELYBrgVNE5FHg68DvxLB0GYnrEuSI\nyIb8oxBWLjkal+0Mek37Ckl+LbD5MudsA04CfoZU582GjrdR97CtR+ytOIetI7ZepnlNL1EPTe8n\nLU96iTQ0XVrbbP7ZNkdsn8WGpv064mmGpqE2+NiIqOqq7nN2iIy94PcBeGjCoelgHYmh6SCo8J7Q\nk+LnSUu+1SbETeK+RBqzM5F7lTpZy4TYG3rYTUNTPWJl0NDDhNiStfaT1gjvobne8L5q31yIV1uo\nIeaF2yGcxfpDCHEQBMGcEkLcD0KIg6Bl2jD08BnL+fIl621b9SWP94T2FZ981nRT5naJZcLQY1ZZ\nzQhOMFuEEAfBDDKOf3MQRI+4H4QQB8EMsYOyoYdP1sq9pq1X7CsheUMPS9YqGXqUEseC+SGEuB+E\nEAdBS5jt0FGUiz4cZLjwnchkzlrjCvEoZ61h7go7qnOXij4cqLabkr2CtWdOSxwGGSHEPUZEPg+8\nn/Qd/DTwu6q6t9uoZhfLRPaZzuNyNLWhhwmxiaQt+THhalq+dAxJ+HYwWH1pS/X+aoR4icHsZ6p9\nl0jZ08OWLx3F4FIsm5v08UwyjJ5nYEevbjoia7o/xDr7fnMXcLqqngH8APhsx/HMNCZwJm62PQ65\n05VQ94pNlIe154U1F1o7Xhi8WfDx5sctZMf4G4PS5xx1PfLPYsfaYxzyzzXJ9Q3KtGXoEXRL9Ih7\njKp+123eD3y4q1jmAV/IYNKexkHqfyazbzzinnu/56bj/cPbbfqeqMXmBdnje5ne0MMPT5ceoz6T\nndPbUh5y+42Dr0fsbTmD6Qiv7v4QQrxx+ATwjaY3ReSTwCfXL5zZZFph8IJmc7peRK3oQ9MX52GS\n09VCtb/N524nzenmhh6+x5vH78XS2rYbATMe2U+aI/4JzUPwr1bHeK9puyGYdomWjzNYPXEd+0FY\nXM45InI38LrCW9eo6m3VPtcA5wKXjuP1ulEtLmeF7YyuvuTXEnv8GuJh1ZfMNSvojtVaXB4lom8c\nc9+nwuJypoke8ZyjqhcOe19EPg5cDLwvDNfng0nLCgYbl/iH7gchxD1GRC4Crgbeq6rhDzEnWAGJ\nJq/pzQwmX3l8VrN31sq9pict1BDMHpaHEMw/IcT95ouk7/G7UjlR7lfVK7sNKQiCtohkrX4QQtxj\nVPWtXcewEdlMnVAFdbLWQYYPJR5NMvTYyXhzxNYj9hnIvkc8ao7Yqi8NGyoRakMP631bIlj4HHdL\n9Ij7QwhxELSAkEw4YHjW9H6aBey46nFs1VZe9CG3uPTreaEWYUvY8haXXojNUEOqfZqEeLOLYVjW\n9H5irrIr4rr3gxDiIHD4NNZJvuTy6khe7LR6vkwS05IQWy/aPzZnr+fOWt5UA2rzBlvzm38WuyEo\nnacUk73nl0tZO94oYiuT2VxOe42DlcT16wchxEHgmPaLzSdHLVTPvWvUIQZ7qTneJMP2tTb88LPF\naOKemzqYQJaKPliMpXM1xeTXI+ev+2SwSQjxaIcYmu4PIcRB0ALL1B7S5kjlh3NHGXpodfwWBt24\nSlnT0xp6+KxpM/TYR7MwHqzayYemTYRXY+gRtEMIcT8IIQ6ClpnWdcqE8TjS3Kwla1kRCZsj9iJf\nEmIvllD3hg9QJ2vtJ7lqjeqdDptDDronhLgfhBAHwQyhwN7qEQTDiKHp/hBCHAQzSr58ySeD+WSt\n0vKlfA43X74U9IMQ4n4QQhwEa4jN947L8awcmrY54tUMTdsccT40PUnPe9LPEqw98fvoByHEQdAS\n9s9ka379kp9l6jnaJl5TPWwd8fbqYb3h1c4R++pLJu6bSKYeTdgcdf5ZbG0yhLFHV1hiXjD/hBAH\nQQuYAxXUyVV++ZKJlRfInK3uYT1hE/VtrFxHbO3ntX1tedGC299uAGx50yF3riY2V5/FD4tb+wvZ\ndghCN8TQdD8IIQ6CIJhD8jXkwfwSQhwELeCHCY8Utr3RRhPeGav003o/3sjDV1/y+/uM2nHaLLGc\nxW3n8m1A9Ia7JK59PwghDoKWsGzkZQaHpk00RxV9WKQ8H2uOWJMOTZvo+6xpb+ixyPA1wlq9bw5f\ndp4w9JgdYmi6H4QQB0HLTGvosY86oWofzfWIfSJYUz3iUc5ai6Ss6VFxLhPLnWaZtnrEInItcDHw\nvKqePuGx5wDXkf5kvw1cpaoqIt8ATq12Ox7Yq6pntRRyrwghDoIZ4iDwQuF1oVx9adg6Yu81HUOY\n/aNlQ4/rSPXLr5/i2K8AVwAPkIT4IuBOVf0N20FEvkC69wsK5DfUQRDMCfnypVIPOeg3y2M+RqGq\n95KtZBORt4jIP4nIQyLyLyLyjvw4ETkJ2Kmq96uqkoT8g9k+Avw6cOOkn2+jED3ioMSLwH+33OYJ\nVbvzwkzFazWAG5ipWMdgnuJdq1jf1EIb3yHFNw5HiciDbnuXqu4accwu4EpV/aGI/BLwZeCCbJ83\nALvd9u7qNc97gP9T1R+OGeuGI4Q4WIGqnth2myLyoKqe23a7a8U8xTtPscJ8xTvLsarqRWvVtojs\nAN4N3JQ6tEBKV5iGjxC94aGEEAdBEAQ5mygkV4nIAvBQtXk7aX74ZLfLycCzbv/NwKXAOWsa7ZwT\nU0pBEATBAKr6MvCMiFwGaZ5XRM5U1WVVPat6/JmqPge8LCLnVXPBHwNuc01dCDyhqrtXniUwQoiD\n9WLUfNSsMU/xzlOsMF/xzlOsUyMiNwL3AaeKyG4RuRz4KHC5iDwCPAZc0nD4p4CvAk8BTwN3uvd+\nkxiWHomkRLcgCIIgCLogesRBEARB0CEhxEEQBEHQISHEwbohIp8XkSdE5D9E5FYROb7rmIYhIpeJ\nyGMickREZnIJi4hcJCJPishTIvKnXcczDBG5VkSeF5FHu45lFCLyRhG5R0S+X/0NXNV1TEF/CSEO\n1pO7gNNV9QzgB8BnO45nFI+Sll7c23UgJaqlJF8CfhU4DfiIiJzWbVRDuY5kfzgPHAb+WFVPA84D\nPj3j1zaYY0KIg3VDVb+rqoerzfsZXH84c6jq46r6ZNdxDOEXgadU9b9U9SDwdZozWzunZKM4q6jq\nc6r679XzV4DHWekYFQStEEIcdMUnGFzmEEzOG4Afue2SvWCwSkTkzcDZpKIGQdA64awVtIqI3A28\nrvDWNap6W7XPNaShvxvWM7YS48QbbFwqq8ebgT+sTC6CoHVCiINWUdULh70vIh8n1T19n87AIvZR\n8c44zwJvdNsD9oLB6hCRLSQRvkFVb+k6nqC/xNB0sG6IyEXA1cAHVHWx63h6wL8BbxORnxORrSQX\no9s7jqkXVHaNfw88rqp/3XU8Qb8JIQ7Wky8CxwJ3icjDIvJ3XQc0DBH5kIjsBt4F3CEi3+k6Jk+V\n+Pb7pHJ4jwPfVNXHuo2qmQYbxVnll4HfBi6o/lYfFpFf6zqooJ+ExWUQBEEQdEj0iIMgCIKgQ0KI\ngyAIgqBDQoiDIAiCoENCiIMgCIKgQ0KIgyAIgqBDQoiDIAiCoENCiIMgCIKgQ/4fu0V33R8puUoA\nAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10dbe2160>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"poppy.display_psf(psf)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<astropy.io.fits.hdu.image.PrimaryHDU object at 0x1138226d8>]"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAEKCAYAAAASK+A+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztfXusddtV12/s1znne95HW1p6L7TGplIaEiKgjaYYC1oV\nxDYgaCRBSC4kECHR4OP6t1GJGEWMNNKoSQVJWi3RYksRUhN664WmaNtLpZZgb7297b1fv++c77z2\na/rHnGOvscaacz32Xvu5xi9Z2e+151przrHG4zfGIOccDAZDd9Hb9gAMBsN2YULAYOg4TAgYDB2H\nCQGDoeMwIWAwdBwmBAyGjsOEgMHQcZgQMBg6DhMCBkPHMdjGnxKR0RQNhjXDOUd1vmeagMHQcWxF\nEzBsBwMAfQDD8EjhsSeeM/i9GByAuXg9E+/x8xmASXictnYEhnXAhMABoQ/gCMAobMPw/hD+Qsut\nj7ggADL1sEwIAJkgiAkAXvxym4TvTwCMw3YdvmvYHkwIHAD6AE7gBcAx/KI/QnZxU0Kgh6IgAIrC\nQEMufn4tBQC/TgmBKfzinwC4Cs8vYcJgWzAhcADoq03e6ZF4fxUhwJ+7mq8d8iaEC2OZq/GYENgO\nzDFoMHQcpgkcAMbhkdXwSdhG4f0p/IWeIrvjS5OAN6kBpLQB7QeQz3ljU0COR5oDY7Vdi2MwbB60\njcpCxhNYP0bIJPwIWVSgzCTQC7/KJyCFgPYHSOcgRwl4oU9hi34TqMsTMCFgMBwojCxkMBhqwXwC\nHQNHAmLmwCohwpQ5YCrf7sOEwAFhAG//H4vn/P4wbNIpWMcxqPVJGfYrcwxK5yA7Kpk5OA7Pr8Rz\nw/ZgQuAAIIlCMcYgOwmZK7BpxiA/Z2egZgxKwpBh8zAhsOcYVWxAPirAeQNlQkCaBhDPZb4AmwBA\nWgjwfmeR/UP9ln9nUYPNw4TAnmOMvE2vFxg/TuEFwByb1wQm4fOYJqA3w+ZhIcIDQiqBiH0CqSxC\nCp/JRd8ki5CFDAuAWBah9AlYAtFmYDwBQylklIAXvNQAUtqAjgjI59IpaBd4+9gYT4CIniSiXyei\nTxPRp4jox1fdp8Fg2BxW1gSI6DUAXuOc+zgR3Qbw2wD+knPu0yW/sRuFwbBm1NUEVnYMOudeAPBC\neH5GRM8BeC2ApBAw7A5SjkHtHwCKfoCUY9CwX2g1OkBErwPwjQA+1uZ+DfXQNlmozDFoZKHDQWtC\ngIhuAXgvgJ9wzp1GPn8KwFNt/Z/BYwBPFmKi0DGKPAEWAj0UeQIp2nAdTaAObZjDg1oI8MZEoSt4\nspAJhM2jlegAEQ0B/GcAH3TO/XSN75vW2AL6AG7ACwEWBFxajB+BYohQb7HoQJ16AmXRgVg6MS9w\nLi0mS4xdhu0CFjJsCxvzCRARAfh5AM/VEQCG9hBbcH31CMSJRBLz8LlT3weKNn4qRCi1gbKNvx/T\nGuR3DJtDG9GBPwngvwP4X8huFH/POfeBkt+YJtAimCSkcwfKaMPsG1iHY1BXFZLaAFCsKsQmgQmA\ndmFkIcOi0Kg0B2K5A5oxuAxtmBmDWhuQpcXYHLDFvhmYEDAYOg6rLGQwGGrBsggPCDqBSIcIm9QT\nkLcQ3UeAnzepJ5AKEZovYPswTcBg6DhME9hznIhNawE6OqCdgto5KEuNN+k7wOHBVIiwTt+BS7EZ\nNgsTAnuMGwBuhkfNGGSaMFAUArHGI230Iow1IIkJAaYRM2NQhjAJnjBk2BzMHNhj1L14MftePm9j\nEsTCiqnndfZj2BxME9hjPESmco/hTYIxlqsxmDIHgGqykGQLxngCZWShMbwJIE0Cw2ZhPIEDAiGf\nQKRbk6fakLXZmjzVhky3JpcJRDYZ1gMjCxkMHYeRhQwGQy2YT+CAoM2BFFlImwNy4/3Iu0OMLBTL\nItRFRWIRAqBYT8DUwu3ChMAew+oJGNqACYE9xRGAW/A8AU0W0uXFNE9AlxjTZCEdKZDQEQH5XFOG\nZTqxjg7I8mLXYcysxTwM7xk2AxMCewq5oGN3d/4cKN71h8gThLQmoF9L9JEVIWGBQCg2JCFki36u\nPh+G1zGNhMduQmBzMMegwdBxmCawpzhH8W4tbfIJMnNgCn+hp+H7ExSrC7XRmlxXFZLVhqtyBy6Q\n+QbO658GQwswnsABQJYWi5GF+LlOIy6rNgw0Iwulqg1PxcYNRzVZSJYYM7QHIwsZopA9CNfZi9B6\nEm4fG6s2bNgvsMpu9f0NDBMCB4ZYPYEB6pkD0gyoQxbSnIEqcyBGFhrDsG2YENhjHKHIE5CEoVgC\nka4noEuMMVIhQiBvCvBrmUFYVU9gijxRSGYQnsN8A5uGCYE9xACeKCTJQpJsU0UWkgIgFh1oShbS\n0QGtEVSRha6QRQZuwJOFHsJMlk3BhEBLKHOmte0g44XOFYRGyLMFdWUh3ZBUEo3KKgtV9R2IVRbq\nI88UlKZGTMjwfubwi56Poe1Gpazt6OPQkY4uwoRAS2BVeBNgW3qErFRXjOWnqwM7+Asu2XosAGYo\nhgarQoR1fQJcWCQWImSegOxcvA5fAWsphiIsRLjHYF+ArjGofQIjZHTculmEQDOeQFkWYUwIaJ/A\nFTxhyKoLtQfjCXQUhHR5MR0dWDdjUEcHYoxBmwjrgxUVMRgMtWCagMFwoDDGYEchzQGODnDugK42\nrCMDq3Ql1r0HYtWG2SfAPQfMHNgNmBDYUxAynsAxvHOwTgciWVBkF7MIz+E5AiYcNgcTAi0gVrJb\n1uqXFXfawG0Ad+CFAJcX48iAzCKUHYg0R6CKNtyk70BZiHCCTBDIDkQyi1CyBi/ghcApgLNGZyUN\nWUFJXxc+hi4nO5kQWAGahZcKrcX68i2Dk/DIrcdO1MZhQi0EdI1BFgjSFJDj52MoMwe0QJACgE0C\nrmnA9Qx4PBPEIxVy36xJAMuHDGP9FyUk67GN67OvMCGwAmZIJ9rI59pOXuX/+DHW+09vPOnld6hk\nDFKIyUVeRxOI+QRiPQrl+OWWOpZVwNdHsxXlmLUQ6GKRU4sO7CEGaG4O7EIWYVNzoGt35LaxUbIQ\nEb0bwHcA+JJz7s01vm9CoCUQ8qYAb7GS46m+A7HKQk37DpQJAa1my5LjvHECkU2M9rBpIfBWeAH+\n70wIGAy7gY3yBJxzHyGi17WxL0M7qJNKzFpBm6nEuudAWSqxYTdgjsE9RVnzkTKewK40H9E8AVlU\nxJqPbBYbEwJE9BSApzb1f+sGoZiVV1VPQKbKLoMb4fEuMsdgkyzCWMhs3dWG+ZibZBGyY/BB+P5F\n6VlJg5mSQ1TXE9BZj12yVzcmBJxz7wLwLsB8AgbDLsHMgYaQnX9TZCG+q0rVmSvncHjsEs3s4hH8\n3R/h8TYyTeAmsqhArO8AawU6d0DXF2wjd0CaBDp3gFX8aRjTGFllJKmt6GhEU+1pAH9eOEw6QD7s\nyddGNlXVmkuXOia3kkpMRL8A4KMA3khEzxPRD7Wx311EL/Kc1Kbt6RgPf5UTH/uPWGtxHe+P/Wdq\nHHofcksJh6r9x/ahFydQ/h91EAtvxqom6esWG28XYGShJdFH/u6V8gnophyrZM/dCo+PoFhotGlX\nYjnusgSilGMwlUAkuQFNuhJrx+BDAPfD9x/WODcaMptSN1lhaJ+A9F0cAnPQKgsdOG6gmEMgyUI6\ndyDmFNSJNfruXJU7kBICemMhwOYPm0SSLMRC4AJeECzrDDRkMCFgWEDe9ev4BDRjMJY7EOMJVPkE\nDuHuuk+woiKGBWKVhBmxpKceihRhqPdSCUQpbcAEwO7ChMAeoo/MJyDNAY4McO6AbEMWqzRc1pC0\njk+gTkNSGXuXbcg4d4C98NIcYJ+ACY7NwIRAQ8gQoXQMxrzcQNFWliy5JmGvEYDHwvMqx6AWAjHa\ncJVPYBXHYMonIB2DUghUOQbvofm5kuzJWFERPhYZxpWOwS6FCPdOCKSq+MgJuQ71k5DF45mVJyv2\n8F1Wfp/HI4tkSJ4AF8uoqqBD8AKAhcAdZIzBWBuyKtpwEyHQtLxYXSEgt6twDJfIFi8LVsaLqF6Q\nt8MjC0XJE9B+D8kPAPKl0HhjNuN5jf9uilghGk3F3lS1o70SAnriAtmFlSePP2tTEMiFJNt6xWip\nOgzlkJ/QXIuP1eMByolDsh4Avx6q91PjA+qbA6toAvxax95jnikpNOZhbLL6kDw2Of4ybUBWLtL7\nYOGnxyP9IH1kQoo/G8LPoar/bopUWFmf0x4240/ZKyHAJyTGc5ffWYf0HKM4gXiCTFE0B7Rg0poA\n28FANXNQhtOAjGvPd0sel1xYPD7A3xG50tA2GYOxasMyd0Bu+nirFuEU2fnkxc4Cho89tuCkORBj\nDF7W+O+mkPM4ls/A49mUKdIVUpTBYEjAeAJLghl58u4aU6FTjsFlVDy2eR9FPDqgKwvp8mIpn0As\nOpDSBmJ3qqrogPYJsM2dIguxY/Ar4fvLVB1mRmeVY1D6BmTPxEOod2BkoY6gj2rasPQJxPwC63QM\nyvqCMkSYog1bWLA9mBAwGDoOYwxuELGwJUM7zdoANx/h/IFjxMOEMpogtQBpvnDqbht9Bzjqwa+l\nU1CaA+x4k9WGL5F1ILpAu81HYs5PfRybCsftIvZOCNQJEfLEa/t/ZQUfHYPWzS0k314uCs6r5wkP\nVE/2PoBXIeMJyMpCsZLjXFlImwNNqg3H0p71wl+22jCbAByLT5Uc58pC9wB8CdWmAvtMWDByHQUp\n9PSxzMSjjuDIikdtmyl1QoSbolvvvBCQxJdUKS8pBPSC47vNskKBkL/rarKQtrPl73g8sfCTJAvd\nRDkr7lUAvgreIQh4xiAXFWlCFuJx5gTBAKCe33rhAHrhNQD0lDowDyfazYF5mLHzmX/t5sBsWhQA\nWhOIkYVYGzhHXrDx+AHghcT5YTZljCykU6hTZCE9Tk0WYu1kWW1hgHgxGul70YJVnkd5vtrGTguB\nWKWcuvX8mIDDJ7puTF5Dt/mSd1veUtoJkBcCQ/jJxcKCpf4YWaUdDa6Sw01HgUz9j41LFxvlffA4\nh/ALvzcABkO/2Pt9tfD72XPAfwZkix4IQmAmns+BWRAG0wkwn3qBwCYBtx4DipqIhBSYfLyX4ThT\npCoZFQHqMwZTbchYYEq+AdeIbJrizMcsczuq6h4CRSEg51jbRVh3WgiwVNZsLiA7WWUFPNrIYuMx\naNVWVt5h6LsMI+YxZ9VZfh6DvEPxd+RC4UmbStxJHhSAyRzoD/3C7YuF78LA+LUeG3/OQsEFATCf\nA7NJphFIASAZkmVJRvwbPkagePyxc6TPJy8cvj4OxQWXakOWuubLaJPa3GCSEMJ4JGuRx8lj02NZ\nV5/EnRYCM3jJe4Hy5pKxu25b8d4xMlU91vRTJhDxmGKxaOamSxWY7yqnSN9hXPh/eaH42DjGLjWC\noXgszSKchsdxUaAtwxMoyx1IZRFeiccynsC9sKVUcT53vP8r5DVIzu2I8QSkYNUJRLIz0rLt0vk3\nch6nqh3J/cd4FuuChQhXACFfyDKVO6AF07IHz4taZhGy/XwD6b4DkpMvk5164rGNLEI+TnkXk0k5\nUgjovgNcblwKgfvi+8uAS4zFFhyPWwuBQyo3bjwBg6HjqCsELHfAYOg4dtonsMuI9R1oUk9gGRvv\nrnjkEOEN5ENqVfUEdKXhdRUVkc7Pum3I2G5+CM+dYJ4APzYB29516wlo+3sdobhdxU4LAY6tckJM\njJhT5oSTpA8OETaNEhzB8wRkZV9JxqkKW8Z8AjwmSRa6h3hp7T6AJ+B5ArKyEHMXtE+gjDG4IAz1\nMn4Ahwl75N9jSN6ABvMCGLOp5w9weJCjA9N5PIFIMwa1T4Bj8rKy0IsAnkf8+t1CnifAfI6mJcdl\nTJ7JTLIC8imah+f4FEohLX0zMWelHI+8ebCjsm0n4c4JAULm9JKxVb3gUncq5gnE8sKB7C5TJelf\nGR7lgpNkoSYdiPgOqSsLyejADXgyjL7rPQHga5AXAlITYLIQj2sx+UfAMATOR0fA4Mgv8uEw4wj0\nAj+gP8hzA5gXQAljMRYinM+CMAjPmSswmQSBcA2MwwqaXAPjcd4Dz4Kay4udIeu4dAsZB+AP1Fju\nAnhNODf8fRaMWhOIVfGJ8QTkvNFkIRZMX46fmgVGyK4RUJ8nEGMMaqHJ0YplIxYaOycEuMWWnNws\nBGKqd9ldl0+epNDy714uGcNtZHcV3fKLJ5e8y7BkT9UYlJ5zvsPIvgCELBwqhQAhowjrMWntZBEi\nPAaOb3gBMAorZxAEwmDoF3yfH/uZRsACAQivewAl3ErO+cUuhQFrB7Op5wzMpp4zMJt6gTC5BqZB\n8o6v/eurC+D4yl+fSxSFPC8IviPfRZ4divAeb1IISPq0jICU1RiUqcQ8b+S4gEwTuUI53ft2GFNK\nCJRpJ/w/UjDxmOQ66GE5U0ljZ4SAJFDU2VK5A9LmkyWu+OLLfaRMAylkNMe+h/wYpE0dyx3gz9hE\n6SNfhgvI71uPQ/6fNIH0+z1k6j3f3Vm97w+zO31MACxMAhYClD3PHYwwAYgAx9+f++9PJ3mTAvAC\notfPSEkA0J95ATEYesHQmxePRZ5PfQ2kOizPndTEtL9Dq958LRhMJpqj+H+xa5CwlHLj1SainIux\nz6V5Ivcvx5RaB6toBDsjBPggzpGdEJaCx8iaV8YSL+Q+YpqANgeqylmzcwrI1DaofU+Qj73LMenF\nzJl12rZjc4BV4PPI787Fdhzel+qtdMhN58DszN9lJxNgHFSN0Ylf9KwR9IIAWNCGB/k7f2NzIGgG\ns2nmF5hNvXbAGsBsBozDRZhM/PPxOF9K7ByZSSC7EMlzoO1hPneslTHmyPuSOGMS6jvSBpecBk1i\neig2oJw+PAvf0+NhLkIdTUA7k7n+gjSb2iqAutM8AS4dzUUymhJzYkJgmdyBMnMglTvA0F5zKQS0\nYzCmXg6QOQZfEd5jc4D9FNIc0OYT7yOXPxByB1gYsFMwly/QLxcC81n+NTsHefHPp8BkGqcNS/VW\nsvJ4crP9fRq+/xIyx2Ds+t1G0TEoHcpaK2Ckcgdk5EIKgTNk86guUrkDdcwBSWCSQqBu5MLIQh2B\nTDDSAkozBnso+lXWnUo8QT5MCBRDhFxLoGm7dkM5Dq6oSMoTn8rYayt3AGH/7ISL1RNYZ+4A4LWQ\nqnoCMndAagNAdcnxbeQO6JoCOncgVk+grDvxDRQdg8vkDqTqCbSlelflDsSiA+vOHdgZISCLaJZF\nB+QFBYqOwRgxR/IEzuEn18tIe1afQDpEuGxRkanYxsji4UBmDrwA4DlkE/MRAG8G8CSKIcJCUZEQ\nFTi+6aMCo5MsRDg8Aoaj4Bc4yhyDPVlgICUF9MFAPNdSIDCE5iJCML32zydj7x8AQojw0kcJrs59\nlOD6qlwIvAjg8wA+iSxM1wfwdchChDyHmD/BqrdMDdaHoslCMSGQChE+jzjuAng8PN4M72meQJ2i\nIprvIufyOfy8WaUYK2PrQuAEwGvhTxpQTwhoJ1yVEOCiFUAmBF4C8P8A/F8xFoKfVFII8HiaCgEt\n3aVKrIuKnMJPdF7knwyPbwpbVAjcBm49Aty4A5zcBE5uZQJgceJOwg9y7KYR0D8GhifA4NhLBQpO\ngn7wIvTCa34uMZ9mj/x8FkgBzq/63vQKvcklhrMr7xWUt32EAw+zmQXB5UPg8hy4OAUe3gdOz/JC\n4PFwHeYAfjO8/3UAvh5FnkCsqIhOkkoxBsuEgFx0Xw7/9xzyGsLXAPhqeP+NFAKy6EsdISCTseoI\ngZcBfAHNfRZAXt4bDIYOYuuawBPwmgB7vjlNltVclqDarl0lRMiefiZb8B3nyTCeJ9R47iAfHVh4\neUPITcbcgSLddjbxITQmzVyP8z6B22G/BK8VsCbA43kS3ifweNB37z4O3H7MawInt4HhnbCDheNi\nBBw/CozCfej4UX/XH90EhreCFnDkNYLewL+O3f3LeMNAUSuYXAZb4MrbAZNLYPIQGJ8D0yvgKty3\nxuf++dUYowtgdA7cuQQmp8DlmdcEzu4BdwKj68aZP98OeRPuSfi77lejyBg8BnAkSFK9XqboyMPg\ncCeTnDikKecNRweY2QlkDMaHyLTJu8jm82NYjSyUChHymNhskvsnAL+H5mhFCBDR2wH8M/jj+tfO\nuX9Y53ePwKt5r0DeHLiJfO08Ge/VQiDGANNqlPQJsEDphfcfRzaxHhfj0UJgsb6C7T068RTcwSgj\n4Og4O4PLbvEEu7oAri+B82DIsfNqrP5XnptX3AUe+yr/2d1XAnceA0ZMJbwD4GYPuPFK4ORx4ORR\n4OgucBSkBi/+o9v+uTQHFrHC47w5IB815OIHAjXwKosNTq+8AJhd+QV/fZYJA8C/vn4AXH4FuHwZ\nuPgycD7HMAzz+KY/v4Pg2aQeMHvgJ/7j4dzcF9fpFVAhwtvA0Ym/TiyouYxa7jDmnuMAZAzH6Tjj\nMVxdAFdX/n/lAmZchvGwEJBz51FkNzQgnjsgiVGLMYlHFgKaX8LmwDHyJugVvJlyH82wshAgoj6A\nnwXw7fC+kmeJ6Jedc5+u+i3Hu3keA36x8aKTXXUWhT0Dsy12s8pRV+dZZV958ZhgMxH/xbiF7K5f\noA33/F335FbmeBsdBU1gkBYCPNFmU79Nrj21d3QpJuV9Px5JUZbn5jb8Xf9OkJR3HgNGj4kv3bwD\n3H51PSEwOAJGtzMBwD6BwXH+Vik1Aw3pDwDC6rla+AQWgmB85l+XCoFHgeO7wPCLwPAUGAKjgT8s\nvkuPr4GLB9k8uQN/bfk1n6M7j/jv37ibCQF9fWLXBshfn/E1MGbn6kNgcOZZjUCeNsxzhcHzlq+f\nzh3guTwA0O8hR9lejEkoWfO5J4AxeUlzU+R4xuLcbFwIAPgWAJ91zn0OAIjoFwF8F4BKISDVffZh\nyZBXzjHY81J9MESuHp4skCkLXs6mXprT3Eta6Qgaq/3L8chNvnfMzrcT4OhmJgAGQRtg4g2PicEc\ne06okReeJ/lkApyc5yeKTApi4XMUBjXiQd1CJgBuvSYIgEe8+n/8iBcEgF/8w5tBCBxnpoDWBkho\nAForWByQcgYCfvHHtIDRLf96eBOYBI0AyAROr18UNvNTYOKPkY/35JY/B8dn+XPDz08QnKO3s+8f\nhd+PjpArqqqvzYyVmpnXAlhgMIsSCDeWs0wl13OFcYTiHDoW3+d5PDoREZrI/F38Z0jEGoyB/jyv\nNUgiEVBcR03QhhB4LXz0hvE8gD9W98+lig/xXLLchija3yw9NattNsku3mAEuCu/8PlA+8hMC9mv\nD8hX5JVsu2GYRMOhXz+DYV4ADGLmQA/A3I93PvfjnAAYzDP/AGsC/UFmI+pS5gOE/xlm6vFiBYwQ\nVPxwhx+e+Dv+6Fa26IG8ABjeDI8sAIQpUPAJlGgCg+O8WcAmAS/qySATKrHf82+m18AkaCdH58D1\n6eL4+HgXx4/ivFhEZgb587nwBQyyebMQAOHaSM1gJhQbIAiII2A4Ddd7AAym2byQcwXivULadvhs\nMe9GWGRzci6HxIKSPQNoKgTEVcYdkGtErhvemmJjjkEiegrAU5v6P4PBUA9tCIEvwDtqGU+E93Jw\nzr0LwLuAjDYsCTRs2/Bz6RAhABQI0y5ko9UxB6bjfOloIM/AYluLIUs7ywIYE06HnQD967yqyJjP\nM3PAiQqW2hzgCAF7owHkGnbIopyLbRI2Jo1zxs0YXsUe3fR3/8klMHiY3dVTd3I3ze7Gg2N/25mt\nwRyYXvvX12d5c2ByDowfej/B5DIUGzjzn4syyny8i+NHfl7w9ZqheD4n18F3FDQwJ80BcW20ObBI\ne2YlJfz3ZJpXv+VcgXhPzq3YvKNxpsX2+/m5nDIHeB7LeTETj0D+3DRFG0LgWQBvIKLXwy/+7wPw\nV+v8kDnjMoQnk4QA4fGfA8Mrbx81dQyyV5X/U26yUoz+DMiiCf3z7GLx5JlNvarJKbR1HYMcHbgM\na+LyPD+ea/H8Cv57lw+B6+DFHF96x6J3M5+KP+PFPcucdYBfbMNbfvFJx+D1LjgGQ3Tg7IvA+amP\ne136Y7wOF+HyoT8HMotOtzS/PPcqNuCvk5wLSzkGA5Hp8mEgM6l5oV9Djwf5zMZFxuccmJ0v5xiM\ntXHn/9frqAlWFgLOuSkR/RiAD8If87udc5+q81vOFLuDYtGPVBrxYB6cJGEONg0RcryXU3fFElp8\ndheZQ4fHM58D0wd+khzfAEbX7YUIz8LGTTgfhmM9FZ/duJ8VCaFeOGeLW+KpX2RXD4CTr+xtiJAP\neHwKnN4DTgNP4OH9/Dni68evbyHY5vez8z0+ASZX7YUIea4wPVfOI8YpstRmDgOyQ5o5K9f82Rzo\nj1cLEcr5K89NU7TiE3DOfQDAB5r+7j483fFYvDdDdqDrIgs9gKeivoR8haGXwyY9rHIfFwBOrnw1\nnDbJQkwbfils8ty8BB8x6T/InEbTsZ+kt86Bk4fA8BzAyRy4+SJw48W9IgstVliELPQgXJyXz7Jz\n8TKyEBifrxGwSM4C/D6Oz9ZDFpL/zfOFwWM8Qua9l9GEdZCF9HiahgeBHWAMPg9/4HKRbjJ3QDLQ\nPo8srstmAnPGo7kD4xC+QTu5Ay8gH2b5PLKSWmMAF+E2dHrm2XTp3IExcPTi1nIH8n6BdnIHuJ6A\nPj93kM0fvksvcgfC9VlX7oAez4Pw3hyZ4Nh07sAysNwBg6Hj2JmiIpZKbKnElkqcjauNVOKDqywk\nCRhWVMSKimhYUZEiDk4IGOIYwMqLGeI4iPJiVmgUVmgUVmh02UKjTce4MzhCvLNO2yHCMvV7hGKl\nI11PYCi2WMRCYorykuPMTXgBxc42r4evVvMqMSYWAjwm7oy06I40ylKdgf0rOc7mAAsBPtYhivny\nr0S8A5FOQW+j5Lj2xpctxhtYTz0B3ampaVu0GHZGCPCJ4FTMqg5EdR2D2ueF8Pk1MntQ4wayi8ct\n0WI+Ae2nqOMYnCATHHI84/AfUggMkNVW4P9H5L2FFtDLOwYXHYhE96GhcAyu0nxEEm0kNVp3IWKK\n9HSSjWf0u7J3AAAUi0lEQVQciFbDa8/EHMzzd2rZtg3wk5+Pd4C8NiDnC58fdgzqZKw6jkEZxpUL\nFWFcMsqZEgJ95OcN0DxEqLXHMfLzjb8/Ft9fFjsjBGLtl/Qm5+WCyaf2IeOrTmzyovOWgvycn8/V\n4yyMgcclxyHo6Yu7jTYLtE09U/sAiv+n9ynfZ0rqdOIX9WwKzESmJfPoZ5P8f/QH/je9PrKOQj2/\nwBdtyNTAmrQhW/QpDNoBgEXPQhYasWORTjt9DST09YY4z+wQZLVaMvjkeZS/kdcods3lmFLQx8L7\nl3NxJt6vakNWtrXlrNwpnCI7SUwbXkdD0jKcoUgblubFOhqSyjJnEL97AK9WniCz8eXdUo7rGsDR\nlQ+3HUJDUr5OTIl9gOKkf4DAqBSf8TkehX3zvJFzRt4kIM6lnje6IakcUxn4c17sfKxtNyRtAzsd\nHRjAWpP3Ya3JrTV53pdU17ltIcI1o8onEOMuSE1gmVDYXfHIzkopCOTElyFCKQTYiam1KrlA9HMJ\neafSz/XdSzrbpBDQIUIW0iwAWBNgzWiZzrscj9f2dxVtWGoC+w4TAgZDx1FXCFjugMHQceycY3Cf\nQGhOYBpjeY8uq/dcylrShtknUMYYlLaoZAvqDLtlzQHtTZe+mTLGoPYJsNPrvvj+MiBU+wRkxIVN\nlq6pqXsjBJra37zg2sodYGKOZuXJ/JtN5Q48iqIQkCQmftzl3AEm46RyBx4iI+ZsMndAO2+lN34T\nuQMMef6WzR1oMqadRR/tMwabXsgRVqMNA8ULypOLw5B8EWKCgOAX/2PIaMPsGFyaNowsOtDnghsr\nhghns8AaFNEBWYevLm1YFlKWgn6KdBLPDZQzBmM3DqAZbXgE76xsqpXwnGiTMdg0wlWFnRYCeuJW\nheOAPKnIoTgBmkrUXGgtMZaqzEaZkDNHUZDFJijEZwPxG94/j6PsDs8XNyYAZC1+Xf9e5w2wUOBF\nD/iFTjPxvOcFQa/nw4QAkid7nhirLCMuNSt5/LFd6jmhox+p6IBkmsqxDZBpNfqaNxUCffEotxRP\ngMcT08CGqCa6LYOdFgLXyKuYKSEgLyZ/X5M+llWnuJ6npLRqGrAcFyNmoujKQhxzvkY6/jxFlmXH\nmoLUQGI2u2a3DcN+BjzeKdCfAv2rTCOQvICcQFCzcR5OtKyMy7yBBWMQ+TChPFdA0Scgy2bxJo+3\nKsOQzx+bP3xO+JzLc1UWIpTjZDNFFw9tCh7zJfKCZYoiwSzGItVkoXWELvcuRNgkd6BN9JGuJ1An\nd0AKgVgWYdV/V9UTYLNJ1xQA4tqMFKbSLwC0m0rMgkD7BGQtAbnYyuoJVN0BY1mEmm9flTuQqiew\njruvvpmlcgeWhfEENghW66r8FG0dNPfgk9WOZFHWFGNQqra8KNhhJgUrULRVJeRk5dd8fFPkhV5M\nE9CMQb7T6yo+VcKxLuT1KWMMHtSkhAkBg6HzOIiiIoZq9FFs3srRlFiIMOXQTPEEUiHCKtqwVK9j\nIUL2dkva8CXaV7sN1TAhsCRkvFeq10B5ZSFWg5eZ7LIYa4wncCQ29rQDxQQivfGiX4UnoMOzsTi3\nbOElO+rorjqaJ7CMWdBHMYGoqrIQmy9t8Uv2BXspBMom6TptO249HUsllmSbqkInOpW4jtf5lcg7\nBmUCkeQISJ8ACwEZm9aCIGYvp84tQx4Xv5Z+Dy0ApIcbyMfi2Smom2qcIWubcA/Fiksx6Pb2dQp4\n8DHEIjg8pjaq98RQFiKUAmrdsNwBg6Hj2CtNoIyYE7vrtmlfjlBk5Uk1W0p17YGW3nK+I8kQWxUt\nlAuuhhohuXAgb9IcSIUIU6QiPoer5A6kzAGOifOxA0U/RE/tk+/E8nhHKI+RD5BpAjJKokOEqdwB\nLkTTD4/6e23G5zWBCSgPEa7bT7JXQoDj0HMU7dc246sxSLtWh9Mky4xRVk9AbkC1/Rmj3epW1VLA\naEefhIts63QMpnwC08ijPj98vPK3KcjvyJqCDJ4zKbJQLNlJj7stSMETq3bE491U2HKvhACQxaK3\n8b/cyVja3SzVtU+AkXIMNikd7eDtYrlPWQFHRwdYG4gVFYlpAusoKlImBDiBiM+BdAxyeTGZRXgP\n9RYDOxB5f3Ucg9InoBOIrmr+7zLYxB2+LownsIfoI51FGBMCdc2BGGOwTAjUjQ6ksghZCKSyCHdl\nkewrjCxkyKFpdEDOnliqa0wTqIoOGDYLIwsZFojd/cuEgFabY7kDc+SFQ0xQyN9IJ5xht2BCYE/B\neQMpslBZPYE6PoE6uQNNfQK6nkCKLMT5A4bNwITAkmBGWircw4hVO1q2xBh3RZLlxWRRkRRtWDsG\ndcXhth2DstJwyjGYog3HHIPL1Nfn0mJNSo5Lx2CXtBYTAg0h6ahaCDShDTdtLDlCvjVarCGpjA5w\nGi2wfHQgxSSTjkCgeXSA06m5SEeqaaiuSVgXfA6a0oZ1dKArwmAlxiARfQ8RfYqI5kT0TW0Napeh\n7WKgGHefRz6Xj3o/TRH7j5h6Lt9L/WdqHHofcktpMVX7j+1Dh+mA8v+oA53mLB/1f8gtNt4uYFVN\n4JMA3gng51oYy17AIct9b1LyjBlpyza2GCPr0iudcHz3OkYWemPGIF/ca2Qai9YGUkVFgHJNgP+f\nX6eKivAx8x0WyDQBWb2Hi3fEuhI3PV9TeM5AVcsvbQ7IcGaXwlcrCQHn3HMAQFQrEmEwGHYQ5hNY\nEswn32S7qgvxeA9px2BZGzK+++s2ZKlU4pRjMJZKLH0CTMNNOQZjbcikY7CNzL1NX599RaUQIKIP\nA3h15KOnnXPvr/tHRPQUgKcajM1QAg6tvZz4vCw6oAWBdgw2FQKx6IA0C1LRAcNuoBXGIBH9BoC/\n5Zz7rZrf75LJZTBsBcYY7AgI+RRnnU4MLFdtWGoAMTZgzByoW22YnYKSLMT5A3Z32DxW0gSI6B0A\nfga+8M19AJ9wzv3ZGr+za70CBigvOa4rC8n6gmVCQDMGGSkhoEOkKSGgqw3rykKpkuNmMqwGSyDa\nAAjpNmRl9QSWPXi28VPlxXQbslTfAclw7InnqQSiVN+BWALRHPnjTfUd0G3IYuXFHojvLwO+PqkO\nRJosxOM8iMkJMwc2ghjzDojH1134DuetL/t/8n97yC/kWIowUFzsMQGgC6XEcgf4/3kBydcx8Hcc\n8tpFaqyxY1kFet96f7IilRxz1zQQEwIrgO90seYjsQ5Eq1JQuSBp7M4m1XAelzQHpuH7Q/g7a5U5\n0FQTKCMLyYo9QLk5IJOIlmn7JcHnIUYb1tTnTVXx2UWYObCnIGQ8gWN4s6AOT0CGBvl5m+XFZHhQ\nluyq4glcIOtA1FYb8K7DfAIdhcyeqyo5XuYTAIoCgKF59imfQJ2S4zYR1gcTAgZDx1FXCFjfAYOh\n4zDH4IFBmgPaJ6B5Auv2CaR4AtIfYCrh9mFCYI/BiUM3kG9DplOJZc/EGFlIbsByqcRyi1UabpJK\nvGpUwNAMJgT2EJIxyEIgxhisKwRiIUL5qBErmlLGGIwJgViIUNcT6Fq8flswIdASNtkkVS5y7j6s\nw4OyK7H8HnMFYgQdPf6q6IDulsNddPrIIgNsWsguProzEJsTw/A9PoY2hUBVUZEumyUmBFrCJktq\n8x2UF/cYRYKPTvXljVmLsthom6nEusho3UKjXBVpLN5vEzw2QxEWItxjHKFYVCTmE9D1BCRZSPce\nYNQpNCpfy8YjkixUVmhU+gRkUZF1tQLvGown0FHEGIMDxKMDm84ijCUQWeWf9cGEgMHQcVgWoSEK\nVv21NpAqKpLiCcSKimgtoMtJOfsEEwIHAPYF8GMsRFjHHFhHyXFpCqRChNfi0bB5mBDYU1gWoaEt\nmBDYU3BUgAWA7kUoNQEtBPhRcgV0WLBJGzL5WmYTci2BAbwQkOPR2Yx638v0HzQsB0sgMhg6DtME\n9hRMxY054zg2n2owOkH+wjvkS20xmWiGemShquYjcpNj15tsVmLYHEwI7CmukTnYmH8viUKxkuPa\nKVgWHajTmrxOdEDmD8RKjjNhSJYW60In4F2C8QQOCIR8dCBVbbhOFmHTvgNlWYRl1YavYE7AdcHI\nQgZDx2GVhQwGQy2YT+CAoM2BWIgwZg7ESqYDzchCzA5M1RPQCURmDuwOTAjsOWQG4QnKyUK62rDM\nIIzxBCBeA+m+A5IxKCMEutpwiix0iXwmoWGzMCGwx5CVhTRtWBYVKassRMizBoHlacNTFIVArLIQ\nNyJhTWAYHnlMRhTaLMwnsMeYV38FQF7djj2vu586Y6n6r7r7MWwOpgnsMS6Q3Xn5jjtGee5AjCfA\nvICYOVDGE9DmgNYCyioL6dyBS5g5sC1YiPCA0EezrsS6qlCqF2GMJ5DqSqxZg3W6Ehs5aD2wEKHB\nYKgF0wQMhgOFaQIGg6EWzDF4wOA6+5IoFOMJyDAhsFxXYg4PxngCOoHIfAC7hZWEABH9FIDvhPfx\n/B8Af905d7+NgRnqgx2CzBWoQxaSVYWqHINVZCHtGJRpxFVkIVlazITDdrCST4CI/gyA/+acmxLR\nPwIA59zfrvE78wm0BK4qxFWGUiXHy1KJOTy4SqFRyRqsSiXWIUKdSmxoBxupNuyc+5B4+QyA715l\nf4ZmSNUGkI+p761aT0AWIeHfpog+sgYh71OOUW+mEWwWbfoEfhDAf2hxf4YKzJCRa/gufIx8nT8g\nq/En+xBWVRuuYw7UqTbMdQbL6gmwSXAJEwDbQKUQIKIPA3h15KOnnXPvD995Gv4av6dkP08BeGrJ\ncRoSmAI4C9sAWaFRfg4UG5JKs0AWFNE+gTLHYMwnoIuKyGrDWgjIXoTWI3C7WJknQEQ/AOCHAbzN\nOVfLpDOfwO5gHY5Bw25gIz4BIno7gJ8E8K11BYBht2Ddeg2rRgc+Cx+Vejm89Yxz7kdq/M5uGFuG\nbEe2jjZkdoG3j01FB/7wKr83GAzbhzEGDwg6i5CLimieAEcIZFUhyRgEin4BCekP4NeyoAg/aqIQ\nmx2c8mxZhLsBEwIHAC4vphmDsrIQhwfrphID5UIAqJ9KzGFCWVkoxhi0WgLbgQmBPYdmCC5bY1A3\nIZUkoVh0QPoEUkJA8hCmav9Qv+XfjWHYNEwI7DnGyJN99ALjx3n4HjcI3aQmwM9jmoDeDJuH1RM4\nILRNFgLS0QEjC+0+rAORIQoZGiyjDQN5k0AiVW04Rhu2C709bCREaNg/sCff7r4GhlUWMhg6DtME\nDhSyDVlZ8xFpDmgzoMoc0G3KY+ZArPmIdBIatg8TAgeCEYpVhVJ9B6RzsMoxWKfvQMwxqKsLVfUd\nMKGwPZgQOACMkCcMlTUk5d4DWgtoqyGp1AY4PMj1DOR4eJP/a4JgOzCfgMHQcZgmcACI1fTjOzxQ\nLBKiOwpVvdaoqwlIbYA3JN63/IHtwXgCB4RUApHMGdhG7oDMIQAsgWhTMLKQoYBUFmFPPGc0ySKc\nIS4IZBahYfMwIWAwdBzWhsxgMNTCthyDLwH4gzXs9xVh3/uAfRorsF/j3aexAusZ79fW/eJWzIF1\ngYh+yzn3TdseRx3s01iB/RrvPo0V2P54zRwwGDoOEwIGQ8dxaELgXdseQAPs01iB/RrvPo0V2PJ4\nD8onYDAYmuPQNAGDwdAQByUEiOiniOh3ieh/EtF/JKJHtj2mMhDR9xDRp4hoTkQ76c0morcT0WeI\n6LNE9He2PZ4yENG7iehLRPTJbY+lCkT0JBH9OhF9OsyBH9/WWA5KCAD4VQBvds59A4D/DeDvbnk8\nVfgkgHcC+Mi2BxIDEfUB/CyAPwfgTQD+ChG9abujKsW/AfD2bQ+iJqYA/qZz7k0A/jiAH93WuT0o\nIeCc+5BzjqnqzwB4YpvjqYJz7jnn3Ge2PY4SfAuAzzrnPuecGwP4RQDfteUxJeGc+wiAe9seRx04\n515wzn08PD8D8ByA125jLAclBBR+EMCvbHsQe47XAvi8eP08tjRRDxlE9DoA3wjgY9v4/72rJ0BE\nHwbw6shHTzvn3h++8zS8uvWeTY4thjrjNXQXRHQLwHsB/IRz7nQbY9g7IeCc+7ayz4noBwB8B4C3\nuR2If1aNd8fxBQBPitdPhPcMLYCIhvAC4D3OufdtaxwHZQ4Q0dsB/CSAv+icu9j2eA4AzwJ4AxG9\nnohGAL4PwC9veUwHASIiAD8P4Dnn3E9vcywHJQQA/AsAtwH8KhF9goj+1bYHVAYiegcRPQ/gLQD+\nCxF9cNtjkghO1h8D8EF4x9UvOec+td1RpUFEvwDgowDeSETPE9EPbXtMJfgTAL4fwJ8Oc/UTRPTn\ntzEQYwwaDB3HoWkCBoOhIUwIGAwdhwkBg6HjMCFgMHQcJgQMho7DhIDB0HGYEDAsELIGDR2DCYEO\ngYj+ExH9dshffyq895CI/gkR/Q6AtxDRNxPRbxLR7xDR/yCi20T09eH5J0KthjeE3/418f7PsRAJ\nNQg+Hvbxa1s8ZEMNGFmoQyCix5xz94joBJ4S/K3w9e6/1zn3S4Ea/Lvh9bNEdAfABYB/CuAZ59x7\nwnf6AF4H4B8DeKdzbkJE/xI+fftXAHwcwFudc7/P/7npYzXUx94lEBlWwt8goneE508CeAN8u8D3\nhvfeCOAF59yzAMBZbUT0UQBPE9ETAN7nnPs9InobgD8K4FlPg8cJgC/BF8j4iHPu98M+TADsOEwI\ndARE9KcAfBuAtzjnLojoNwAcA7hyzpU2BnbO/Xsi+hiAvwDgA0T0w/D9Sv+tcy5XvYmIvnMd4zes\nD+YT6A7uAvhKEAB/BP6OrfEZAK8hom8GgOAPGBDRHwLwOefcPwfwfgDfAODXAHw3Eb0qfPcxIvpa\neJPgrUT0en5/7UdmWAmmCXQH/xXAjxDRc/CL/Rn9BefcmIi+F8DPBL/BJbz28JcBfD8RTQB8EcA/\nCL6Fvw/gQ0TUg+9E/qPOuWeC0/F94f0vAfj2TRygYTmYY9Bg6DjMHDAYOg4TAgZDx2FCwGDoOEwI\nGAwdhwkBg6HjMCFgMHQcJgQMho7DhIDB0HH8f9uCVNNjsPYzAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10e97a5f8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"osys = poppy.OpticalSystem()\n",
"rect = poppy.RectangleAperture(width=0.25 * u.meter)\n",
"osys.add_pupil(rect)\n",
"osys.add_detector(pixelscale=0.05 * u.arcsec / u.pixel, fov_arcsec=5)\n",
"osys.calc_psf(display=True)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"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.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment