{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2018-11-04T07:56:17.196018Z", "start_time": "2018-11-04T07:56:17.125812Z" } }, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "ExecuteTime": { "end_time": "2018-11-04T08:07:31.861611Z", "start_time": "2018-11-04T08:07:31.859084Z" } }, "outputs": [], "source": [ "b200 = np.loadtxt('b.txt')" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "ExecuteTime": { "end_time": "2018-11-04T08:07:37.311599Z", "start_time": "2018-11-04T08:07:36.926098Z" } }, "outputs": [ { "data": { "text/plain": [ "Text(0,0.5,'Reflectance (%)')" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjEAAACBCAYAAAAxOuuoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl0VOX9P/D3nSWzZJJA9hUCSUhQZJElCQQRfyoEUbEsIqBUe/z97HKkB0//aKse28r3tKen2G/rObWlVltFUUBBggYRRMlCgJAo+xJCFhLIPpNZM8v9/ZFmkkBCMpNZMpn36xzOIXfu8nme3Nz53Oc+93mEtLQ0EUREREQBRuLvAIiIiIjcwSSGiIiIAhKTGCIiIgpITGKIiIgoIDGJISIiooDEJIaIiIgCEpMYIiIiCkhMYoiIiCggMYkhIiKigMQkhoiIiAISkxgiIiIKSExiiIiIKCAxiSEiIqKAxCSGiIiIAhKTGCIiIgpITGKIiIgoIDGJISIiooDEJIaIiIgCEpMYIiIiCkhMYoiIiCggMYkhIiKigMQkhoiIiAISkxgiIiIKSExiiIiIKCDJ/B0AEbmuqqrKb8cWRREAIAhCv58Hcqd1hrP9cPTs504xDvezW9dJT08fUWxE5F1siSEiIqKAxCSGiIiIAhKTGCIiIgpITGKIiIgoIDGJISIiooDEJIaIiIgCEpMYIiIiCkhMYoiIiCggMYkhIiKigMQkhoiIiAISkxgiIiIKSExiiIiIKCAxiSEiIqKAxCSGiIiIAhKTGCIiIgpITGKIiIgoIDGJISIiooDEJIaIiIgCEpMYIiIiCkhMYoiIiCggMYkhIiKigMQkhoiIiAKSkJaWJvo7CCIiIiJXsSWGiIiIAhKTGCIiIgpITGKIiIgoIDGJISIiooDEJIaIiIgCEpMYIiIiCkhMYoiIiCggMYkhIiKigMQkhoiIiAISkxgiIiIKSExiiIiIKCAxiSEiIqKAJBvOSosWLUJubi6SkpJQXV2NrVu33nF9tVoNq9XqkQCJiIhobJPL5TAajS5vN6wkRqvVorCwEKmpqZg8efId11Wr1Vi/fr3LgRAREVHw2r59u8uJzLCSmMrKSgBAZGTkkOv2tMBs37496FtjwsLC0NnZ6e8w/Gqs1YEQooRo7QJEh0vbuV4PAtSxiQAE2I1NgCCBVBXt/NRh0UG0dkIamgAI3U+FRYcFDmMLJMrxEGRq57p2QyMEqRIS5fg+23dAtBoh1STCZjaiS9fmUnncIVUoEBGfArPJ3B1Dlw5i161l6ILD2NyvDMbmBpfr23UCVNEJECQSiKINDsNNCCHhkISEOdcY8PfgLEM8IEj7l0ExDoI8tHd7ww0I0hCoxyf21oGlA6LVAKkmybmeaDfDYWqFRBUFQars3V5/HYI8FBLFuN5l5lbAbu0+fk9MVj1EixYSdSwEibx7n14uk8Vkh2i1uFTjvro2SBRqOCyu3+H7giBTQKNWoVPX4e9QPCZEKaDLLLq0jVwux/r1693KGYaVxLhDqVRCKpV6a/cBwWq1QqlUDr3iGDam6kAix33/Lws3L+tw8asalzZ1tR6mPDIHKantsHQpIHZaAEECQRPbu4JZCtEiQAiPcX7RwGaGaOiCoIoC+nxRiTojBJkKUPfZ3giIVjmEiFioVQYU7RoPY8N1l8rkCokMWPR0OKyycADh/y2DbJAyWJxlkMutaGqYiLN7yr0WGwCkL5mFielaWLqUgMMGsdMEQREJ9En8RH0XAOGW30NPGWJ7y2A3Q9T3lCG8d3udEYJMCaj71IFRgGiVQYjos0+rAaLRCiE0GnAmoyJEbSeEkAigT8IBgx2iowtCWJ/tu5QQTSIETTQgVXQv82KZJCFq2OQJKH3nrEvJpi+uDWETJmLRahO+3R0K3bVqrx7LZYIE2RvvgcUkQ+VH3j2/fUUdFYrslePw9bbrgAt5jFwud/uYXktiOjs72RIzxloh3DGW6mDmqoVouwnExDSgepyI1jr9sLd1pR7GJ6mRmNyKfW9cgb3fHWT9AGu3DHMZAFwbYFkTkmbPwOzlKnz+hgFweKfFI3fd3bh4og4NZ84NUA+Dl0GQKbD0+TgoIsxoqffO9SQiRoqk2AvY979G2E19f6eD1ePIfg9hYW0D1EGzC/usGmD5jRHFNNIy3fd/k5CcNw3nPy8eZP+38/a1QZAAC+/X4ruTk3FvficK/myF2GX22vFclblkPmxWEaHyBmhiu9BY5VpL1mgjCMD81TKU7qpFp861lq+RJDF8O4loGCJT4zEhTY/jH5WjtFCGvHXpkKo0Hj+ORKHE/EdV+PbfZ29JYLznevl3aLrUgLtyvXNXnJI7GxFxGlwsc708os2Coo+uY/6qRMg1ER6PTRKiwsINGSj+tOOWBIZcUfrhGWRONyAsMc7foTjNzE9F/dk2nCv4BjXft2HGQq/ds7tMHSFBZlYtSj84j/I9NzD30Rgoxsf4O6wRmblqIVo74tBc49tHd0xiiIYgkQK5S7pw9L3LsJuN6KhtQvWlMGSvvdfjx5q3dh7qGiZB22z3+L7v5NTn1zF5TiyisjI9ul9lZAzm3G9D0Y56wGFzax+dbXZUXUlB9trpHo0NAOasnouG67Fob+zy+L6DibWzA6W765G7tAuC4O9ogLjp05B0VzROF3efc6cPXEP8lPGIn3G3nyMDIJVj4bpJKNt7ExZdG8x6EefPpGDBU1MAjILKc0PMXVORMsmAiv21Pj/2sJIYiUQCmUwGqVQKQRCc/ycKBnPX5aH22ji0N/TeqX//2TFEjDcjec40jx0nac5MREaZ8f3+7zy2z+Fy2IFjB5RY8HgopEr10BsMU86j4Th9XAl948j625zdX4rQUCsmLpjjociA+Jn3ICbBiIpPyjy2z2B283IL2lojMOOJXL/GIVOHI3eZDMW7252djUURKP1Chpx8mVda9Fwx4/Ec6ExJuHmt9/Ho5a+OAQAyHszxV1hukyo1yH1UieI9Bp+1Hvc1rCRm2bJlePPNN7Fy5UpkZGTgzTffxKZNm7wdG5HfJc6eiahYC04fuqXPgN2Ko+9X4d6cG5ArRn73FBKqwtyHRBz9+AZg909fspYLF1BXFYo5azyTKEy6Nxp2YwuuHPJAkiDacfSjGsyar4NCNfL6lskFzLtfh6JdrRBtgd0XYTSpONCECRlmRGdl+S2GeU+k4so5Ddqv9u87pL1WjctnNMjxQovecEVNycTEDD2O7zh1yyciij+4iLvv7YAqLLAaCGbeL8HVSgNaL130y/GH9ZCwoKAABQUF3o6FaFSRh4/HvIdFHN7eBNhuf9ygb2rFmWPRWPDMbBzZVg6XuuPfYu4SOU4fakRnfdMIIh65yj1lWPrjyYidqEBTjftf7prEZMxcEoX9b172WGzGG404tb8NOU9E45sPWjCi+l6RjPNFTdBeM3ksPgLsZiOK9xqx4PFQ7L8WCrvZ4NPjJ6QroZLWomSfbsDPz+4/hod+fDeSs8JQf8G3LxxIpEDuMhGl+8ywm2/vf2XpaEHZLh0W/GAcvnpP6/bjV19KuvcuRKd04MA/L/ktBvaJIRrE3AcFnCu1QFc7+OvUV0+2QhoSgrQHst0+zuRFc6GIiEbVcf8mMAAg2rr7/uQ8Ng7y8HFDbzAQiQx5Tybi5BEZujo927xce84CISwNmQ+7/8giNW8uwhOTcPnk6HlTZSxpvXARNWdtmPmAb4dWUEbGIGdVOkr2mgZ/1Vu04+i/z+Deh5RQRUb5NL4ZD0bhWvl1NJ87P+g6jVe7YEQa7nl09D9WCgkfj7lL5Cj5QjaS+4kRYxJDNIDU3LuhCJXj0pELQ6wpovjDK5iebYImPtHl46hjEzBjfhdK942eN2P0HXZcvpKG3LXu9feZlj8TnVoF6kq9M/ZF6UeXMXW2AWHJE1zeVhUVi1kLLTi6owZ+vfKOcd/t/x7RiRIk3JPuoyMKmP/UFJw5EQaT9s4tiGaDA2e+T8WCdZnOQRa9LX7mdCROm4izRUO3TB3feR6Tp+oxPi3NB5G5S0DuurtxviIUnfW+78zbF5MYoluoYuIx634ZjhUOr++Fua0Jpw47kJPv4oEEKfKemoiKoyEwtfi/Faav84VlUCjtmHTfPJe2C4+SInXiNZR94L3BuywdrTj5lRR5q2Ncfplj/qoEfHdMBePNRu8ER91EoPRLNbKXh0IePn7o9Uco48FsiCJw+dCxYa1/9ZuTsNkEZD7s/RYPuSYCOUslKNrdDnEYebPN2IljBTYs+ME4CNLR+RU9YXokZIIRFw8Mr769aXTWEJG/CFLkPTUJlcUKmJoHGkBsYDVlp9HZYsDUxcO/88yYGwpDmx7Xik64E6l3iXYUfViNmfONUGiGNxCVEKLEfT+8C6X7TLCZvfuWQv2JU2i5fA1Tc4b/JtXkmSqYW67h6tfHvRgZ9dDV1uBCpQbzn7ob3nx1WB0uwdS7qlHywQUMv3VNROkH5zH1Xj1Co8KGXn0Ecp68C5fOhELrwojBN8+cQd13dbjnPpUXI3OPKnIcZi4SUPRuhQ+mAxkakxiiPjJyY2Exiqj+1vUvuvLDIjKzNRg3aehm4IiUFGRlq1C246w7YfqEqfkGTuypwvzHQjGcL6HZK+fhRsM4tNb5pjPnqcJWpOcmIzIjY8h1w5InYPqSyThe0Ak+RvKdC4WlCJFbMXFWgncOIJFh4dOZOPG5DpaOVpc2tejaUPpxFeY/KvNajpU0JQRy6xWc2+96i8X3XzUhZVo84qaPgrFteghSLHw6C9+dTILFNDr+jpjEEP2XJnYcsmZ1oeQ/7r0SbDPocOwLO/JWjocgVwy6niBTIO/JOJw4GgubdXRcCAZTf8EEW8gETM2/c0fauOnTEJ9sQPlu37Uq2W1AWaEUC54IgyTkDnesUjny1sTj5OEQWC3+v3MMKqID3/7rOGYssEAR7vlWhXsey4XeFIPGK+69ZXazygCtMR7Tn8jzcGSAKjoWcx9PRcnuVkB0ffBK0QGU7AdylskhU4cPvYEP3PVIDkxGOWpKfT+W1WCYxBABgCwE9z2dgZNHI2Hrcj+xuPHdadxsUGHG0kmDrjPrB/PQfEONG5Wn3T6OLx37pA6ZM40InzBxwM+lMmDeojYU7Wrz+dw0TefO4/o1DeasnjvoOtMfz0ZHmxL1Jyp8GBn1sBhFfFcWhUU/nNE7oaQHRKanY3KWHmU7RvZ7Ld9Xh0kZBkROmeKhyAAIEix4Kg3fnxgPs8H9xLn9ahWqzmmQvXam52JzU3hCFKZM0+PYh9/7O5R+mMQQAZi5IhutTSo0fj/ycU1O7ihCQkITIpNuf8U0MiUcSRM7cfLjUdgPZhBWXTuOHxCRtzoWgvT2/jHzVmXi0nEtOqr9M0twxSdliBpXj5iU22OLiJFiQmIVjn906+Bi5Es1Jy7BbJZhar5nOtIKEiBnqQ2lBRbYjCMb78Vu0qNknxkLVmggVXimtShzSQ5sVimuHhl5/6sz+0oRqriBhHTfvrLelyAAucscKP6oGla91m9xDIRJDAW9qLQJmJCmx8mdJz2yP9EhomS/iIUbsiBV9XYalCpCkLtchm/fPQNHV2ANstZw6js0XWnGXbn9E4UJubMRkTAOF0/4b9Rb0WZB0Yc1yP1BXL8h5aUKFRaun4ziXS2c3HEUKP3wDDJn6BGenDzifc1ckoT6M81oOnPOA5EBzefOo/ZMJ6Z7YJLI0AgJsqbWo/SDc/BI/yvRjqPvX8HcZeOhjPTPJJGzVi1AS0skmqsHHkTQn4Y9d9KaNWvwpz/9CVu3bsWGDRsgk42eGUGJ3CWRAjkP6nH0g2uwe/CNGu31dly9qEH22lnOZfOeykFtbQq0zaN/JM6BnCq4hkmzohGd2d3sroyMwewRTu7oKZ1tdlRVTew3pPzs1XPR0BCP9hv+mcaB+rN2dqD0kybkPmwc0SSRcdOnIWlaAk4XeTZx/v6Ly4hPD0P8jKnu70QqR97aFBz7tBEWXZvHYjN1OnD+bCrmr/X9JJExd01F0iQTKr5o8Olxh2tYSczSpUuRkZGB3/72t3j11VeRmJiIFStWeDs2Iq+buzYXNVc1aK9r9/i+T+8rQ8Q4M1LmzUTM9LswLtKM04WeuXP0B4cdKD2oxvwVGkiVauQ8HoUzx1UjntzRU87uPwZVqBWpeXO7J3eMN6LiE75OPZrcvNiIliYVZqxwb4RrWWj35I5Fu9sh2jybOIsiUFKoQE6+wu1JImesyIHWPLHf5I6e0jMGTsZDvhvNt2dyx5I9+lHbejys5pS8vDzs2rULWm33s7CCggL86Ec/wu7duyEOMnqPJmECrDYbHF1aiBYdpJpEZ6cu0d4Fh/EmJMrxEOQa5zb2zgYI0hBI1NHOZQ5zO0SrHtKwFOcy0WaCw9QCiSoagkzVZ/s6CHINJMrewZUcxhaI9i5Iw3pHUxWtejjM7ZCo4yBIQ/670A67vgGCIhySkN4T2G64CQCQhsb17nOYZVKqVbChakyVydXfk1JihkMVOyrLFJ2Zjqh4G8o+Gv54MC6xW/HtB9fwf56UQBoah6/ea/fb5I6e0nrhAmqnLcTSn6RDV1eDy4eGGtHYh0Q7inbU4OGNEXDYuvDNLj0ndxyFKr9qxyObYtHRNA2Oy1fgMLVCoo6FIO15o0+EvbMeQkgYJIreqS/sxibMfjwdV86p0XG10iuxdU8SuQDz18/A8Q9KIMiUkKh6pydwmFoh2kyQhvU+EhNtRjhMrYhIvwcT0/XY/7/emkeoe5LI/Bfi0HxaAZsQ6xxxWLRb4DA2QaKKgiDrHTvJ3nl9gDK0QbQZByyDRBULQdb7e5g6qxPV31vQesl/cyMNZcgkRqVSITIyErW1vUML19bWIjQ0FOPHj0db28BNZjlPxEEUHIBFAdEsgxAeDwj/PZzdAlFvhaCKBUJ6Xx0TdRZApoCg7h1TQDRJga4QCBF9xhmwGSAaHBBC4wBZaO+6WgMQEt69355lRgA2C4TwPtt36SCaBAiaeKDnD0e0QdSZISgjAUVk7/b67mxf0PTZ3qUyRY3BMrnye2octWVyyCUo39eBMJUCwOCvRI+I0YTz30ohqPWAth1hYd4dWMsXqg6ehfqJabjwzXW3yuPVOjDqcabQBFkI4Gi1j9r6Hq1x+Ur5FzZMe0BE1sxYiEYRgiYOkP73Rkd0QNQZAcV4CMreL1/R4IC+uRY1R2xerb/ab89g3MpZuG91JCALhaDuvTESjRLAZuh/nbJ2QjSKQChw6oAdarkEkLsW37DLY7egcu91ZOeHdcfQM22C3QRRb4Ogju13bFFnAmRqCOr43mUmCWC9tQz6238PcEBfa8CFg3VeP1/l8uENqDmQIZMYpbK7R7TJ1NuUZDQa+302kG+2VcJq7XvXOdB75YP1ch7o2dtAw5gPtL0WQN0AywcaCGmw7Qd6y2KgYeHvXKawsDB0dvb0nB8bZepv6DJ118FAd+ujtUyed+0sEBYmgdk8ViYcNOPE9lK3tuz/N+EdNYPPrzcq+KIORjvz+Qs4Ud+3Hgb7G7/mu6D6OP5+8X//1w6gfoA1hnvtHpqr58P1K93/gIFeLR8ohnYAAz3y9VwZRspud30cnR5DJjE9F16VSgWDoXskTrVa3e+zvnoyqvXr17sdFBEREQUXuVx+S+PH0IZMYkwmE9ra2pCSkoKWlhYAQEpKCgwGA9rbb+8MaTQasX37dpcDISIiouAkl8udT3lcMayOvUVFRcjPz8fVq1dht9uxfPlylJaWDtqp151AiIiIKDi52/AhpKWlDTkaj0QiwapVq5CdnQ2JRILy8nLs2LEDNg+/4kZEREQ0XMNKYoiIiIhGG047QERERAGJSQwREREFJCYxREREFJA8Ootj3w7AgiDg1KlTQdcBeOPGjZg7d26/Mr/xxhuoqanxY1Tet2jRIuTm5iIpKQnV1dXYunWr8zOFQoF169Zh+vTpsNlsKC4uxp49e/wYrXfcqQ42b96MSZMm9RvU6ZVXXhmTg57JZDKsXbsWmZmZCAsLg06nw5EjR3D48GEAwXGdGKoOguk68dRTT+Gee+6BSqWCxWJBeXk5PvnkE9jt9qA4F3rcqR6C6XwAul+nfuWVV6DRaLB582YA7l8XPJrE9J0o0m634yc/+QlWrFiBXbt2efIwo94333yDnTt3+jsMn9JqtSgsLERqaiomT57c77M1a9ZAo9Hg17/+NVQqFTZt2gStVouvv/7aT9F6x53qAAA+/fRT55fYWCaRSKDVavGXv/wFLS0tSEpKwosvvgitVovy8vKguE4MVQdA8Fwnjhw5gt27d6OrqwsajQbPP/88lixZgs8//zwozoUed6oHIHjOBwB49NFH0dHRAY2md04+d88Fjz5OysvLwxdffAGtVgu9Xo+CggLk5uZCGMm86xQQKisrUVlZCZ1O12+5XC7HvHnzsHfvXhiNRrS2tuLgwYNYsGCBnyL1nsHqINh0dXVh3759aG5uhiiKqK+vx+nTp5GWlgYgOK4TQ9VBMGlsbERXV5fzZ1EUERMTAyA4zoUed6qHYJKSkoJp06ahsLCw33J3zwWPtcS4O1HkWJSTk4OcnBxotVqUlJTg0KFDgw4MONbFxcVBKpWivr53/pHa2lokJCRAEISgqpdly5bhkUceQWtrKw4dOoSysjJ/h+QTEokE6enp+PLLL4P2OtG3DnoE03ViyZIlyM/Ph1KphF6vxyeffBKU58JA9dAjGM4HiUSCDRs2YMeOHf2Wj+Rc8FgS4+5EkWPN4cOHsXv3bhgMBqSmpuL555+HKIo4dOiQv0PzC6VSia6uLjgcDucyk8kEqVSKkJAQWCwWP0bnO59++qnzTiwrKwvPP/88LBYLKisr/R2a1z355JMwmUw4duwYwsO7Z0MPtutE3zoAgu86ceDAARw4cADx8fHIzs6GTqcLyu+MgeoBCJ7z4cEHH8T169dx6dIlTJkyxbl8JOeCxx4n9Z0ossedJoocq+rq6qDX6yGKIqqrq3HgwAHMnj3b32H5jdlsRkhICCSS3lNNpVLBbrf3a1od66qrq2E2m+FwOHDu3DkcPXo0KM6LlStXIi0tDX/9619ht9uD8jpxax0AwXuduHHjBurr67Fx48agPBd69K0HIDjOh+joaNx///3YvXv3bZ+N5FzwWEuMqxNFBou+LRDB6ObNm7Db7UhOTnY2FaakpKCxsXHMNZW6IhjOi9WrVyMrKwtvvPEGDAYDgOC7TgxUBwMJhvOhhyAIiI2NDbpz4VY99TCQsXg+ZGRkQKPR4NVXXwXQ/QafUqnEH/7wB/z97393+1zwaMfenokiIyIioNFohpwociyaPXu2s/lrwoQJWLp0KSoqKvwclfdJJBLIZDJIpVIIguD8v9VqxfHjx/HYY49BpVIhKioKDz30EIqLi/0dsscNVgcqlQrTpk2DXC6HIAjIzMzEwoULx/R5sWbNGueXt16v7/dZsFwn7lQHwXKdUCgUyM3Ndd5hJyYm4pFHHsG5c+cABM+5MFQ9BMP5cPLkSbz88svYsmULtmzZgvfeew8WiwVbtmxBTU2N2+eCR+dO4kSRwEsvvYSkpCRIJBJ0dHSgpKQEBw8eHHN/lLdavnw5li9f3m/ZpUuXsHXr1n7jxNjtdhQVFY3JcWIGq4N//OMf+OlPf4r4+HgAQGtrKw4fPoySkhJ/hOl1kZGR+J//+R9YrdZ+4+JcuXIFb775ZlBcJ4aqg2C5ToSEhODHP/4xJkyYAKlUis7OTlRUVGDfvn2wWq1BcS4AQ9dDsJwPfU2ZMgUvvPDCgOPEuHIucAJIIiIiCkicdoCIiIgCEpMYIiIiCkhMYoiIiCggMYkhIiKigMQkhoiIiAISkxgiIiIKSExiiCjgbd68GRs2bPB3GETkY0xiiMaIjRs3YtOmTf4OY9iioqLw1ltvIS0tzd+hEFGAYhJDRF4llUr9HQIRjVEemwCSiEa3uXPn4oEHHkB8fDzsdjuqq6uxc+dONDU1Aeh+JHPz5k1s376933a/+93vUFZWhoKCAgDAnDlzsGTJEsTHx0On06GiogKfffaZc1byzZs3o7m5GR0dHcjLy4NEIsEvfvGLIeOLiorCli1bsG3bNmRnZyMrKwtarRYFBQUoKytzrhcZGYn169cjIyMDer0eX3755W37kkgkWLZsGXJychAREYHm5mZ8/fXXOHr0KABg8uTJeOmll7Bt2zZUVlYC6B4GfdOmTfjb3/6GM2fOuFHDRORrTGKIgoRMJsPnn3+OxsZGqFQqLF++HD/72c/wm9/8Bna7Hd9++y02bNiAXbt2wWKxAAAyMzMRFRXlnLAzNzcXq1atwscff4yqqiqMGzcOa9euhUajwbvvvus81uzZs3H8+HH8+c9/hiAILsW5YsUK7NmzBzt37sSCBQvwzDPPoLq62plsvfDCC3A4HNi6dStsNhtWrlyJlJQU5+cA8PTTTyMlJQXbt29HU1MTUlNTsX79etjtdpSUlODq1asoKCjA008/jZqaGnR1deG5557D4cOHmcAQBRA+TiIKEqWlpTh9+jRaWlpQV1eHf//734iNjUVqaioAoKKiAlarFXPmzHFuk5eXh7Nnz6K9vR1A9ySXe/bsQVlZGVpaWnDlyhXs2LEDOTk5UKvVzu20Wi0+/PBDNDY2oqGhwaU4jxw5gvLycjQ3N+Ozzz6D1WpFZmYmACArKwsTJkzAv/71L1y7dg319fV4++23+z2yioqKQnZ2NrZt24bz58+jtbUV5eXl+Oqrr7B48WLneoWFhaitrcVzzz2HZ599Fh0dHfj0009drlci8h+2xBAFieTkZCxfvhzJycnQaDTO5ZGRkaiqqoLdbkdpaSny8vJQXFyM0NBQzJw5E//85z8BABqNBlFRUVi9ejVWrlzp3L6npSUmJgY1NTUAgNraWrdn4K2vr3f+3+FwoLOzE2FhYQCAhIQEdHZ29mt10ev1uHnzpvPniRMnQiKR4Je//GW//UqlUjgcDufPoijinXfewWuvvQaJRILXX3+93+dENPoxiSEKAnK5HC+++CKqqqrwn//8BzqdDgDw6quvQibrvQwcPXoUDz74IJKSkpCZmQmDwYDTp08D6E08JNZ0AAACl0lEQVRWPv74Y1y8ePG2Y/S01gBwPo5yh81m6/ezKIqQSCTOGIZKjnri/OMf/+jsp9N3X32lpKRAoVAA6E7mWlpa3I6biHyPSQxREEhISEB4eDj27t2LGzduAOju3NqTHPRobm7GxYsXkZeXhylTpqC0tNTZOtHZ2Ym2tjbExcWhqKjI52UAgIaGBoSHhyM2NtbZGhMaGoq4uLh+rUBAd1LSk4ANJDw8HBs3bkRhYSFCQkLw7LPP4vXXX4fBYPB+QYjII9gnhmgMUSgUSE5O7vcvLi4Ora2tsFqtWLx4MaKjo5GZmYk1a9YM+Pjk6NGjyMvLQ0JCwm3Jyt69e7F48WLk5+cjMTERcXFxmDFjBtatW+eT8l24cAF1dXV49tlnkZqaiuTkZDz33HOw2+3OdZqbm1FcXIwNGzYgOzsbMTExSEpKwvz58/Hwww871/vhD3+IpqYmFBQUYM+ePdDpdHjmmWd8Ug4i8gy2xBCNIZMnT8bLL7/cb9mNGzfw2muv4Z133sGKFSswf/58NDY2YufOnfj5z39+2z4qKythMplQV1eH1tbWfp+VlZXBbDZjyZIlyM/Ph91uR0tLCyoqKrxarr7eeustrF+/Hi+99BL0ej0OHjwIuVzeb533338fDz30EPLz8xEdHQ2z2YyGhgYcOXIEAPDwww9j4sSJeP311yGKIux2O95++2386le/wuLFi/H111/7rDxE5D4hLS3Nvd53RDQmqdVq/P73v8c777zj0+SEiMhVbIkhIgDdA8SFh4dj2bJl0Gq1zkHgiIhGKyYxRAQASE9Pd462++6777r9ijQRka/wcRIREREFJL6dRERERAGJSQwREREFJCYxREREFJCYxBAREVFAYhJDREREAen/A5/3VAlSYwtnAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVgAAACkCAYAAADSSzaWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xt4VNW9+P/37Jk9yWQyCbkQSCAQQgIBFEQuCYiI1gsoaqgoFiz4Oy2VQ/1pS3u0PbXoQamVI+BR6rHHU7W1ICKoKD0IEiRyCUECVYiQEG5JSIDcMzOZydz294+YMUNuE5jJTJL1eh6eh+zrZyXwyZ611/os1fDhwxUEQRAEn5MCHYAgCEJvJRKsIAiCn4gEKwiC4CciwQqCIPiJSLCCIAh+IhKsIAiCn4gEKwiC4Cea7rzZLbfcwpQpUxg0aBBnz55lzZo17R4bEhLC/PnzGTt2LA6Hg/379/Pxxx97vV8QBCHQ1NHR0c91182ioqIoKiqipqaGqKgocnJy2j12/vz5RERE8PLLL5Obm8ucOXMAOHfunFf7BUEQAq1bn2D/+c9/AhAdHd3hcbIsM3nyZP7zP/+ThoYGGhoa+Pzzz7nlllv44osvOt3fnrCwMOx2u0/bJAiCIMsyDQ0NrbZ3KcEOGDCA2NhYZFnGZDJRUlJCY2Ojz4JseR+1Wk1paal7W3FxMfHx8ahUqk73K0rr2b9hYWEsWLDA57EKgiAArF+/vlWS7TTBxsTEcMsttzB58mQiIiI89rlcLoqKiti7dy95eXk+CzQ0NBSbzYbL5XJvs1gsqNVqtFptp/vbSvrNT67r16/3+inWYDAgZazAfPRNHJXfXmOreh6DwYDRaAx0GAHRl9sOvbT9ai397nwFU87LOGrPtHtYV9suyzILFixoM690mGB/+MMfMn36dAoKCti6dStnz56lpqYGu92OXq9n0KBBjBgxgszMTGbNmsVf//pXSkpKvA6sPVarFa1WiyRJ7iSq0+lwOp3YbLZO93fEbrd7nWDtdjuS3YYDbZ/sWujK96q36ctth97Zfk3EcBwuCUvVaXA52j3Ol23vMMGGhoby7LPPUldX12qf0Wjk5MmTnDx5kk8++YQbb7yR+Ph4nyTYS5cu4XQ6GTx4MMXFxQAkJiZSXl6Ooiid7vclxW5BJYf59JqCIHQ/TcwIHDVnOkyuvtbhONgNGza0mVzbcuTIEQ4dOtTxzSQJjUaDWq1GpVK5/34lu93OoUOHuO+++9DpdMTExHDHHXewf/9+r/b7kuKwoNLofH5dQRC6lyZ6BI7qgu6959WcJEkSAwYMQKVSuZ8mvXH33Xcze/Zs99fr1q2jsLCQNWvW8Pjjj1NUVMRnn30GwKZNm5g/fz5/+MMfcDqd7Nu3z2OEQGf7faXpCVYkWEHo2VRoolNpPLOzW+/a5QSblJTEz372M2RZRpIk7HY7f/nLXzh16lSn527bto1t27a1uW/dunUeXzc2NvL222+3e63O9vuK4mhAEk+wgtCjqSMSUWlCcVQXdet9u5xgH3roITZs2MDx48cBuO2221iwYAHPPfecr2MLCordAuIJVhB6NE3MCJz1pSiO1mNV/anTWgQ///nPiYyMdH8dFhbm8bR66tQpwsPD/RNdEFAcVtEHKwg9XFP/a2H337ezA/Lz83n66afZtm0bBw4c4ODBg/zmN7/hm2++Qa1Wc+ONN3LgwIHuiDUgFEcDkniCFYQeTY5No+H4hm6/b6cJds+ePRw7dowFCxYwadIk3n33XUpKSkhLS0OlUrFp0yb3FNjeSLGLUQSC0JNJ+gFIoVHYK090+7296oOtqqri1VdfZdq0afz6179mx44dbNmyxd+xBQXF0QAiwQpCjyXHpOE0laM0ejfk1Je6VA923759rFq1ijFjxrBs2TJiY2P9FVfQUOxWJDHRQBB6LE3sKOyVJwNz784OSEhIYOHChQwcOJCSkhL++te/8vrrr5Oens6yZcvIysoiKyurO2INCDHRQBB6Nk1sGpb89wNy706fYBcuXMiZM2d48cUX+frrr5k/fz4Aubm5vPjiiyQnJ/PUU0/5PdBAUewN3000UAU6FEEQukgKi0Oti8FeFZgn2E4T7MCBA9m9ezeXLl1i165dDBgwwL3PaDTy5ptvsnNn986O6E6KwwKAShMa4EgEQegqTWwaTtNFFGtNYO7f2QHnzp0jMzOT/fv3M2bMmDaLufTqUQTuBKtz/10QhJ5BjkkL2NMrePEE+9e//hWVSsXcuXMxGAxs2ND9Y8kCSbF/l2DFWFhB6HE0saNwBOgFF3jxBFtTU8Obb77ZHbEEJ5cdxeVApREjCQShJ5F0sajDYnEE6xOsLMtdulhXj+8pmipqiT5YQehJNP1H4TRfwmWpClgMHSbY559/nltvvRWdruOPx8OGDeOxxx7j9ttv92lwwUJxNIihWoLQw8ixY7BXBHappw67CF5//XUefPBB5syZw6lTpzh37hy1tbXuJWPi4+NJTU3FYDDw+eef8/nnn3dX3N1KrGogCD2P3H8MDcf+HtAYOkywxcXFrF69miFDhpCens71119PTEwMWq0Wo9FISUkJu3bt4tChQ35ZXTZYiMkGgtCzqA2DkEL7YQ/wYqVe1SIoLi52r33VF4mCL4LQs2j6j8FRdx7FFtiVca9qyZirJUkSc+fOJT09HZVKxZEjR9i4cSMOR+tFyF555RWPr2VZpry8nBdeeAGARYsWMWnSJI9z165dy/nz530et+IQy8YIQk/S1P+aH+gwujfBzpw5k9TUVFasWIHT6WTp0qVkZmayefPmVsf+4he/8Pj6mWee4fDhwx7bsrOz+eCDD/waMzQnWNEHKwg9gkpCjh1F47nA10jpUjWtazVt2jS2b99OXV0dJpOJbdu2MWXKFFSqjuf5JyUlER8fT05OTjdF6kmxmVHJ+oDcWxCErtH0Swa1jL2qe1eQbTOW7rqRTqcjOjraoy+3uLgYvV5PVFQU1dXV7Z47depU8vPzWy0hnpGRQUZGBnV1dRw4cICsrCwURfF57C6bEY229y6LIwi9iab/mKbFDZ2Bf/HebQk2NLRpoL7F8v18/oaGBo99bZFlmUmTJvHOO+94bN+9ezdbtmzBbDaTlJTE4sWLURSl09KJBoMBu93uddwGgwGV2oE6NBKDweD1eb1FX2xzs77cdui57VcPHItSU3BN8Xfl3I4mWHU5wQ4aNIi4uDiOHz+O3W5Ho9HgdDo7fXK0Wq1A05Os2WwGmhZQbLmvLRMnTsRms3Hs2DGP7S2Lzpw9e5YdO3aQnp7eaYI1Go1eJ1iDwYDRaETWVaDXhGE0BvaNZHdrbn9f1JfbDj24/eoQoiKTMR3biOMq4+9q232SYMPCwliyZAkpKSkALF++nMrKSh5++GGsVmubL6paslgsVFdXk5iYSGVlJQCJiYmYzWZqatovJXbTTTeRk5ODy+Xq8Pqd7b8WLpsRlTacppqwvu+CEATBN+TY0SgOK47a04EOBejCS64HHngAl8vF7373O2w2m3v7kSNHGD16tFfX2LdvH7NmzSIyMpLw8HBmz55NTk5Ou0+/AwYMIDk5uc1VaydMmODuWhgyZAgzZ87k6NGj3janSxSbCZVKEiMJBCHIyQPGYq84Dor/Hri6wusn2NGjR/Pf//3frZ42L1++THR0tFfX+OyzzwgPD2f58uVIkkReXh4ff/wxgHulhJblEG+66SaKioq4fPlyq2vNmDGDBQsWIEkStbW1ZGdns2vXLm+b0yWu7wYrq7QGFLvZL/cQBOHayXFjsRRuDXQYbl4nWL1e734p1VJISIjXb+5dLhebNm1i06ZNrfa1VWf2ww8/bPdaq1ev9uqePuGworgcSFoDLvPF7ruvIAhek/QDUOvjsF8+1vnB3cTrLoLi4mLGjBnTavuUKVM4e/asT4MKRorN9F0/rCAIwUiOG4ujrjhgy8O0xesn2H/84x8sWbKEfv36IUkS6enpJCQkMHbs2FbTWnsjl82IKkQkWEEIVnLc2KB6eoUuPMGeOHGC119/neTkZFwuF3fddRcRERG89tprnD4dHG/s/ElpNCJpe+a4QEHo9SQZOXYU9svfBDoSD10aB1tQUEBBQeCnnwWCYhddBIIQrDQxIwEFR3VhoEPx4PUT7NChQ0lKSmq1PSkpiSFDhvgypqDkajSiEk+wghCUtAPGNa1e4GpdmS+QvE6w8+bNIyYmptX26OhoHn74YZ8GFYwUm0l0EQhCkJIHjsd+yT/j4K+F1wk2ISGhzVqrxcXFxMfH+zSoYOQSowgEIShJhgTU+gHYLvbgBKsoSpuLH4aFhXVabrA3UGxGJJFgBSHoaAeMx1FzBqWxrvODu5nXCfb06dNtrhp7++23c+bMGZ8GFYwUm+iDFYRgJA8cH5RPr9CFUQSffPIJy5Yt49///d8pKChAURTS0tKIi4tjzZo1/owxKHzfRSAKvghCsFBpw9FEp9Jw7N1Ah9ImrxNscXExL730EnfddZd7RldxcTHvvPMOZWVlfgswWCg2o7vgi6hHIAjBQR4wDpe1Bmed79fi84UujYMtLy9vVfi6rxAFXwQh+GgH3Ig9SLsH4CoKbqvV6qYq/1e82OqopmuvIAq+CEJwkTTIcddjOrwu0JG0y+sEGxsby49//GNSUlLaHDWwdOlSnwYWjBSbSdQjEIQgIfcfg4KCvfLbQIfSLq8T7COPPIJOp+Ott96itrbWnzEFre9XNhAEIdC08ZOwX/xn0M3easnrBJuUlMSqVav6xAut9oiCL4IQJFRq5PgbMf/zrUBH0iGvx8HW1NQgSV4f3is1FXwRCVYQAk0Tm4ZK0gZd9awreZ0xP/roIzIzM90rwfZFTQVfRBeBIASaNn5SU3J12jo/OIC87iKYO3cuERERrFq1irq6OhwOz36PZ599ttNrSJLE3LlzSU9PR6VSceTIETZu3NjqWgCLFi1i0qRJHvvWrl3rrofQlWv5imIzoY4Y7LfrC4LgDRXa+Ak0HG+9zFSw8TrB5ubmXvPNZs6cSWpqKitWrMDpdLJ06VIyMzPbXfI7OzubDz74wCfX8gWXzYRGPMEKQkBpolNRyXpsQVg960pdWjLmWk2bNo3NmzdTV9dUlGHbtm385Cc/YcuWLV4vnOiPa3lLFHwRhMDTJkxuWprbYQ10KJ3q8kSDq6XT6YiOjqa4uNi9rbi4GL1eT1RUFNXV1a3OycjIICMjg7q6Og4cOEBWVpa7qldXr9XMYDBgt9u9jttg+P6llkrtQAqJ8NjW2/Wltl6pL7cdgrX9KjSDM3AWfejX+LpybVmW293ndYJVq9XMnDmTSZMmERMTg1qt9tjf2USD0NBQACwWi3tb8zLgzfta2r17N1u2bMFsNpOUlMTixYtRFIWsrKwuX6slo9HodYI1GAwYjUb312p1BRGyHqPRRF8o+HJl+/uSvtx2CN72a/qPwaAOpf7sfnA2+uUeXW17RwnW61EEs2fP5uabb+bLL79EURS2bt3K/v37MZvNvP/++52eb7U2Pc63rCnbPCKheV9LJSUlmEwmFEXh7Nmz7NixgwkTJlzVtXzl+4Iver/dQxCE9oUMmtJUmtBPydXXvE6wEydOZP369ezevRuXy0VeXh7r169n+/btJCcnd3q+xWKhurqaxMRE97bExETMZrNXdQxcLpfPrnW1vi/4IvphBaHbSTJywiRsFw4EOhKveZ1gIyMjKS0tBaCxsdH99PjNN99w/fXXe3WNffv2MWvWLCIjIwkPD2f27Nnk5OS0+VJqwoQJ7o/7Q4YMYebMmRw9evSqruUzDiuK04YUGum/ewiC0CZ5wDhQXNgvBffkgpa87oOtra3FYDBQU1NDVVUVqamplJSUkJCQ4PF02ZHPPvuM8PBwli9fjiRJ5OXl8fHHHwMwf/58ADZsaBrbNmPGDBYsWIAkSdTW1pKdnc2uXbu8upY/uSw1SKFRfr+PIAieQgZPwVb2FSjOQIfiNa8TbEFBAePGjaO4uJgDBw4wb948Jk6cSEJCgtdjZF0uF5s2bWLTpk2t9jUn1marV6++6mv5k8tajaRrvbquIAj+o9LokAfcgPHgy4EOpUu8TrDr1693/33fvn00NDSQkpJCbm4ue/fu9UtwwchlqULSRQc6DEHoU7SDMnA11uOoPBnoULrE6wQbFRXl8QLpyJEjHDlypM19vZnLUoPaENzLlEv6gchx1yP3H4PaMAiVWguSjKuxDpepHGd9KbbLX+OsOUNfGG4m9HwhQ6ZjK9lLT/v36nWCXblyJU8//XSr8WF6vZ6VK1f2iYLb0PQEK/cfE+gw2iSFxxOWNhc5YSLO+lIcFfnYyg6hOKzgsqMKjUYdHo8mZiShI+5FsTdgu3CQxnNf4DSWBjp8QWiTZEhAE52CKe/1QIfSZdc8k0ur1fq1wEqwaeqDDbIuApWEbvQ8QpPvwn7pn9R/8btOE6ZKo0OOG4t2yHQibl2Jo/oU1lP/wN4D5ncLfUtI4nTsFd/iaqgIdChd1mmCvfvuu91/v+OOOzwG8kuSxPDhw/tUEW6XpRoptB+o1EHxNlMlhxM+6eeowxMw7v8DjupCr85THBZsZbnYynKRwuIIHXY74RN/jtN8CWvhJ9jKDvk5ckHwgkpNSOI0GvKDv3JWWzpNsBkZGe6/jx8/3mOcqcPhoLKyko8++sg/0QUhl6WpzoGkiw74b1QpNBrDTf+OYqunLns5SmPdVV3H1XCZhvwNWE59QmjynYTd8C+EjrgPzm8H4z4fRy0I3pMHjAW1jK38cKBDuSqdJtjly5cD8Mtf/pI///nP7jn/fZViM3432SAqoAlWJYcRPuXfcJrKMH31Gri8L2DTHsVmwnLyQ6yndxCacjehoxcRMeQOLCc2Y798zAdRBxl1CCqpqaaGorh6RHWmviZkyAxsFw4GfWHt9njdB7t27Vp/xtGjuCw1gR0LK8mET/4lOBsxHf6TT5JrS4rdjOXEB2gu7cMZfyvhk3+Bo/YcloKPcFQc9+m9/Eklh6GOGILaMAi1IQEprD+SLgZJF4VKE+ZOrs0UlxPFbsZlrUVqrEJXV4qj9hzO6iJc1vYrtAn+IelikQfeQH328kCHctW8TrCLFy+mpKSEzz77zGP7XXfdRWJiIv/7v//r8+CCVaDHwupv+ClSaCT1e5/3b9ELu6mp6+D0dnSp92JIX4aj9izWU59iv/Q1wTVkRoU6IhFNzAg00alo+iWjDh+I4rThNJY1/ak9h608D8Vag8tuRrFbvl+RVGoq4iPJ4Ui6KEKjElGHJxCSOA0pNApnQyX2S183/ak47vNfakJrIcNuw1FdhLPufKBDuWpeJ9jU1FS2b9/eavvx48e57bbbfBpUsAvkSALt4Klo4ydQl/17FFv3lJNTrDU0HPsbllOfokudTfjEx3FZKrEUbf/u41sgKhupUEcMRhM7Gjl2FJqYNFSyrml4WnUhlsKtOGrP4DKW05VfBM2vLbXVBkzfDUmUdDFoYtOQB9yAfsJjoND0grB4L46aIt83TQBJJmTIDBqOvRvoSK6J1wlWp9PR2Nj6P5Ldbu9zCyG6LNWoDQndfl9JF0PY2IU05L+Hy1Te7fdvSrTvYjn5ESHDfkDYqLmEXTcfW8l+Gkv24aw949f7S+HxyLGj0cSOQo4dhUobjrOuGHvltzSez8ZRVYDi8P07ApelClvJfmwl+5uWix5wAyGJN2GY9jucdcVYz+zAdiE3KEaV9BbaQZNBcWEr/yrQoVwTrxNsRUUFo0aNoqLC88XOqFGjqKys9HlgwcxlqUaOu66b76pCP/5nOKqLaDyX1c339qTYTVgLt2I9tQ154A2EDL2ViJuX47JWYy87jL3iOPbqwmt7aaQOQdMvCU2/5KaP/DEjkUIicNSX4qg8gfnrt3FUnkSxm3zXMG8oTuwX87BfzEMVGkXosB8Qdv0j6NIewFLwMbbS/aB4V/xIaF/osDtoPL/n+y6cHsrrBJudnc2cOXOQZZkTJ06gKAqjR49m9uzZbN261Z8xBp2msbDd20UQMnQG6shE6nb/tlvv2yHFib08D3t5HiptONqBE5DjJ6AfegsqdQhOYynO+hKc9aW4rDW4rLUodjO4nCiKE5Vai0oTikobjhQajaSLQR0ej9owCEnfH1wOnHXFOKpPYS49gKOqEMVWH+hWuynWGiwnNmMp/JTQYbcTdt2P0I24l4bj74kJG9dA3S8Zdb9hmA69GuhQrlmXEqzBYOD+++/ngQceAJrGwX7++efs2bPHX/EFJZf1u8kGkqZbfsOq5DB0o+ZiObH5qse6+ptiM9FYnE1jcTaoJNSRSWj6DUMdMRh5wDik0H5IoVGoNJ5L+iiOxu/e3NfgslThNJZiK8vFWV+K03ihZ3zsdjZiLfoH1nNZ6FLuIXzS4ziqCjAf+zsuU9+ZhOMrupR7sF3I7RUjN7o0VXbbtm3s2LGDhISm/sfy8nJstp45Pu1auCxVAN02FlY3cg4uay2N577w+718QnHhrD3TTp+squkXk0pqehPfmz5OO6xYTm6h8fwedGN+ROStK7EWfoLl1Kc9/qNud5H0ccgJE3v00KyWvF7RoJlarUZRFEpLS/tkcoWmpzXFaeuWkQRSeAIhw27/7m1qMA2LulpKU2J1Nvau5NqCy1KF+fA6TLmvoB0yncgZL6Du1/mySgKEDr8bR8W3PXpoVkteJ1hZllm0aBFr1qzh6aefJiqqqar//PnzmTVrlt8CDFbd1Q8bdt187OVHcFT1rDqYAtgvf03d7t9gr/iWiJt/j27kD5tqWAhtUoVEEDLkZixF2wIdis943UVw3333MXjwYNasWcPjjz/u3p6fn88999zT5hjZK0mSxNy5c0lPT0elUnHkyBE2btzYqhqXRqPh4YcfZuTIkRgMBurr69mzZw+7d+92H7No0SImTZrkce7atWs5f757fvO5LP4fC6uJTkXufx11WU/59T6CHzkbaTj2N2wXjxA+fjHygLGYDv8p4HUsglHosDtxGstwVOQHOhSf8TrBjh8/nrfffpvTp097bC8vLyc2Ntara8ycOZPU1FRWrFiB0+lk6dKlZGZmsnnzZo/jJEmirq6OV199lcrKSgYNGsQTTzxBXV0deXl57uOys7P54IMPvG2CT3XHZANd2g9pLNmLq+GyX+8j+J+j4jh1X/w7+vE/JWLGCzT88y+iYlkLKllPSPIdNPzzrUCH4lNedxFERES0uWqBWq1Gkry7zLRp09i+fTt1dXWYTCa2bdvGlClTUKlUHsfZbDY+/fRTKioq3P29x44dY/jw4d6G63dNT7D+q0egiUlDE5OGtaBvDYHrzRS7GdOh/8Jy4gP0Ny4h7Pofiy6D74Sm3N00oaOX/dLx+gn24sWLDB8+nOpqz6ETN9xwAxcuXOj0fJ1OR3R0NMXFxe5txcXF6PV6oqKiWl23JUmSSElJYefOnR7bMzIyyMjIoK6ujgMHDpCVldXpst0GgwG73ft55AaDoc3tKpcRtWFgu/uvlXrMg7jKc9BrGsFP9/CGv9rXE/it7ZW5OI+UE3LdYrQxKTiP/y801vrnXteg2372cjia4XfhzH8HgyG8e+7Zia60XZbldvd5nWB37tzJww8/TEhICAAjRoxg+vTp3Hrrrbz55pudnh8a2jT+0WKxuLc1lz5s3teeefPmYbFYOHjwoHvb7t272bJlC2azmaSkJBYvXoyiKGRldTzLyWg0ep1gDQZDqyVymmkun8YwYh5GU4PPx2pqYkdhiEym/qvXcVm6p95AWzpqf2/n97Yb81FVPIN+whI0E3+D6dB/eV0svTt0589eN2Y2KuMF6s8GR+3hrra9owTrdRfB4cOH2bJlCzNnzkSr1fLII48wYcIE3n33Xb755ptOz29eCUGn07m3NdcwaLlKwpUeeOABhg8fzmuvvYbT+X0iKykpwWQyoSgKZ8+eZceOHUyYMMHb5lwzp7EUlaRBHT7Q59fWpd5LY/Fe93hboXdS7CZMB1fTeP4LDDf9lpCkHwQ6pG7XNN34diwnNnd+cA/UpYkGOTk55OTkoNfrkSSpS1neYrFQXV1NYmKiu3ZBYmIiZrO53RVpH3zwQdLS0li7di1ms7nD67tc3TumUnFYcDZUojYMbppx5CPqyCQ0/cdgzvo3n11TCGYKlhObcdQVEz7+Z6gjEpvGPPeEGWw+EDZqLo7as72zoDtXMdEAwGw2X9XHh3379jFr1iwiIyMJDw9n9uzZ5OTktNlv+tBDD7mTq8nUuqDHhAkT3F0LQ4YMYebMmRw92r3zv531pagjEn16zdDUe7CVfYXLLEYO9CX2skPU712BPGAchqlPo9IGR1+kP6n7JaMdfFOPL0nYkQ6fYJ944gmvL/Tqq50XZvjss88IDw9n+fLlSJJEXl4eH3/8MdA0YQFgw4YNREdHc9ttt2G323nhhRfc5xcVFbFu3ToAZsyYwYIFC5AkidraWrKzs9m1a5fX8fqC01iCOmKwz64n6ePQJkymPvtZn11T6Dmc9cXUZz9L+OQniJj+H5hy1/j001FwUaG//sc0Fn/Za2ZttaXDBFtb69s3my6Xi02bNrFp06ZW+zZs+H7VyOrqapYsWdLhtVavXu3T2K6Gs74UbcJkn13v+2mC53x2TaFnUWz1GA/8Ef24f8Fw83LMh/+E/XLn7zh6Gu3gKUiGBCyHevdSVB0m2AMHDnDmzJlu79/sKZz1paj1A0Adcs1V/ZunCRpz1/goOqHHcjkwH/0fQo33EJ7+Sxry36PxzM7Oz+shVHIYYWMexlrwMUpj8JSf9IcO+2CXLVvmftP//PPPo9fruyWonsJpKkNxOVEbBl3ztUKH3d7rpgkK18Za9A9MX60jbNSDhI19tNdMSgi7bgEuax3WXvRLoz0dJtiGhgb3NNiYmJhWM676PJcDl/nitb/oUmsJGXY71tOd13MQ+hb7xTzq9z7f9PJryr+hknv2Q44cNxbt4KmYj77ZJ0ZKdNhF8PXXX7Ns2TJ3X+zTTz/dbnfBs8/2zRczjvoSNBGDuZbCjSGJN6M4bU3rOgnCFZz1xdR/+RyGyU8SMf05jLlre2Qhb5VGh/6Gn2At/ARnfXGkyruGAAAUF0lEQVTnJ/QCHSbYv//973z77bfExcVx7733cvTo0TYXPuzLnPWlyLFp13AFFaHDZ9F4ekef+I0uXB2lsY76/S+iv+FfiJj+LObDr2O//HWgw+qSsLGP4rKZsBR+EuhQuk2HCVZRFHf1qrS0NLZv3+4x1VVoSrChw+646vPl+AmoQgxYz/eQ1QqEwHHZMR/5M6EpdxOe/gssJz/EeurTQEfllZCkH6CNv5G67OV96kHC65lca9d+P5xCr9d3OrOqr3DWlyCFRqLSRlzVgny61HtpPLf72lZgFfoUa9H/4awrQT/x52gikzD9882g/vej7jeMsOsWYD76PwFZbj6QvJ7JpVKpmD17NqtXr2bVqlXul1+ZmZncfPPNfgsw2LkaKlAcjVc14UDT/zrUEYOxnv7MD5EJvZm94hj12cuRwgcSOf05n4xk8QdVSAThk/5/Gs9/ge3Cwc5P6GW8TrB33XUXU6ZMYdOmTa2KrkyZMsUvwfUMCo7qU8j9r+vymboR99F4PjtoV4oVgpur4TL1e/8DR80ZIqb/B9rBUwMdkidNKIaMf8NlrqAh/71ARxMQXifYjIwM1q9fT25ursdIgrKyMuLi4vwSXE9hu5iHNv7GLp2jiU5FE52KtegffopK6BOcNsxH/4eG439Hf8NP0N/4GGg6Lv/ZLSQNhsm/AMB06JU+u6qu1wk2OjqaixcvttrudDrRarU+DaqnsZcfQW0YhBQe7/U5oan3YivNESUJBZ9oPL+Huuzfo44YQuSMF9BEpwYuGEkmfOLjSGH9MR58GcXRd1+Me51gq6qqGDy4dT/jqFGj2ky8fYnLWo2j5gzaeO/q0WqiUpDjxmLpIW+AhZ7BZSyj/svnsF88iuGm3xF23QJQd+/Dj0oThmHKU6gNCRj3/6HPd395PYogOzubhx56yL0aQFxcHGPGjOH+++8P2MKDwcRWfhht/ESspzpbclhF2HULaCzO7nNvVIVu4LLTcHw9trKv0I//KZEDx9Nw7F3sl/w/ZlbSxRKe/gtwOanf+zyKrW+uhtGS1wl2z549hIeHs2TJEmRZ5vHHH8fhcLB9+3ZycnL8GWOPYLuYR9joh1CFRqFY2y4gDqAdlIHaMAhjL68iJASWo7qQuj3PoEudTfikJ7BX5NOQ/57ffqlrE9IJu+FfcFSewHTkjaAeNtadurSiwbZt29ixYwfx8fFIkkRZWRk227VMEu09XMYynKZytANvpPFcO+uCqbXoRj+EpXBrr68iJAQBpw3LyQ9pLN5L2HXzibztj9hKD2Ap+MhnBd1VoVGEjZ6HNmFSU9Wvs91bkznYdSnBAtjtdo+VYaFpAcTCwuBZsC1QbOV5aBMmtptgdan3gqL0iSpCQvBwNVRgOvRfqPslo0v7IZE/+E/sF4/SeDYLe8VxoOOVmNui0kYQOvwuQofPxFFzhvovn8NZX+L74Hs4rxNsSEgILpfLY0XWxMRE5syZQ1paGkuXLvVLgD1JY8leQpPvRJs4DVuJ5wqZ2kEZhKbOxnhwNbi8XzZcEHzFWXsG08GXUUcOJSTpB4RPfhLF0YCt/Aj2i0dx1BSh2NufoakKiUATPYKQxJuQB9yA01SO6at12C9171JNPUmnCTYyMpLFixeTnJyMy+UiKyuLTz/9lB/96EdkZGRw7NgxXn75Za9uJkkSc+fOJT09HZVKxZEjR9i4cSMOR+sxcp0d25VrdReXsQzz1++gH/f/4TRewFl7FgBN9Aj04xfT8PU7OCqOByw+QQBw1p2n4eu3sORvRB54A/LAGwmf+HNUsg5nQ2VTCU7Fht5iBE0oklaPFBaHWh+Hq7EeW3kexn0rcdQUBbopQa/TBJuZmUloaCibNm3ixhtv5I477iAlJYXq6mpWrFjBpUuXvL7ZzJkzSU1NZcWKFTidTpYuXUpmZiabN7desrezY7tyre5kK9mLJnIohsm/wFL0D9SGQWgTJmM9vZ3G4uyAxiYILSmOBmylB7CVHsCMCil8AJrIJKSw/kj6KBSXhGKpwll3HpdlH46aU2Ixzi7qNMGOHDmSv/zlL5w+fZqjR4/y4osvcuLECbZt62w4UmvTpk1j8+bN1NU1jY3btm0bP/nJT9iyZUurlWU7O7Yr12pJlmWv45VluUvHN7MXbkYbFkH40Kk4jRdwFG7GUfzlVV0rkK62/b1Bn2x7YxXK5SqcgNpgoOGKlaPVgLoPfE+6+rPv6FivuggqKioAqKurw263u0sYdoVOpyM6OtrjBVlxcTF6vZ6oqCiqq6u9PtZisXh9rWbN34QFCxZ0OfarYweaVwRN/u6PIAi9lSzLHu+owIsEq1KpPGoPKIpyVUOzQkOb5ke3rCfb0NDgsc/bY5ufUL25Vsv969evb/UNEARBuFayLLtzUEtejSL42c9+5n55JMsyjz76aKtE9eqrr3Z4Dau1aeCxTqdz15JtXlCxeZ+3x3blWi219Q0QBEG4Vu09uHWaYA8e9KzhmJt7detGWSwWqqurSUxMpLKyEmga5mU2m6mpqenSsYqieH0tQRCEQOk0wf7tb3/z2c327dvHrFmzOHPmDE6nk9mzZ5OTk9PmS6nOju3KtQRBEAJBNXz48G7LSC3HrkqSRF5ennvs6vz58wHYsGFDp8d6s18QBCHQujXBCoIg9CVe14MVBEEQuqbLxV56E1mW+f3vf094eDjLli0DeuYU3a5YtGgRkyZN8oh37dq1nD9/Huj97Qe4/vrruffee4mLi6OxsZFdu3bx+eef9/q2v/LKKx5fy7JMeXk5L7zwAtD7f/aRkZHMmzePESNGAHD69Gnee+89amtr/db2Pt1F8MMf/pCkpCQGDx7sTrB3330348ePZ926de4puGfOnHFPwe1sf7BbtGgRDQ0N7RZJ7+3tHz16NAsXLuTtt9/m1KlTaLVaoqOjKSsr6/Vtv9IzzzzD4cOH+eyzplWNe3v7lyxZgkql4q233sLlcvHjH/+YsLAw1q1b57e299kugsTERK677jr3P65m06ZNY/v27dTV1WEymdi2bRtTpkxBpVJ5tb+n6+3tv++++/i///s/CgoKcLlcWK1WysrKgN7f9paSkpKIj4/3KJbf29sfGxtLXl4ejY2N2O12Dh06xKBBTcud+6vtfbKLQJIkHnnkETZu3Oix3R9TdINRRkYGGRkZ1NXVceDAAbKyslAUpde3X6vVMmTIEA4fPsyzzz6LXq/nzJkzvP/++1it1l7d9itNnTqV/Px8dy2P3v6zB8jKymLChAkcO3YMl8tFeno6x44d82vb+2SCvf3227lw4QKFhYXu/hjwzxTdYLN79262bNmC2WwmKSmJxYsXoygKWVlZvb79YWFhSJJERkYG69ato76+ngcffJDHHnuMP//5z0DvbXtLsiwzadIk3nnnHfe23v6zh6Y+16lTp7J69WoASktLeeWVV/za9j7XRRAbG8uMGTPYsmVLq30tp+A2a2+Kblv7e4KSkhJMJhOKonD27Fl27NjBhAlNq+H29vY3x/jFF19QVVWF3W5n69atJCUluett9Na2tzRx4kRsNhvHjh1zb+vtP3uVSsWTTz7J2bNnefLJJ3nyySc5efIkjz/+uF/b3ucSbGpqKuHh4SxfvpyXXnqJxx57jNDQUF566SXi4+PdU3CbtZyC23IKb1v7e6KWhXw6a19Pb7/VaqWqqqrd2X69ue0t3XTTTeTk5PSpn31YWBgxMTHs3r0bu92O3W7niy++IDk5GUmS/Nb2PtdFcPjwYfLz891fJycns2jRIlauXInZbO71U3QnTJhAfn4+VquVIUOGMHPmTPbs2ePe39vb/+WXX3Lbbbdx4sQJjEYj9913H+fOnaOurq7Xtx1gwIABJCcntzkFvje332w2c/nyZWbMmOGuZX3rrbdSXV3t1//3fXqYFjQt2LhkyZI2x8H2xim6v/rVrxg0aBCSJFFbW8uBAwf4/PPP3f9Qenv7VSoVmZmZTJ06FZVKxenTp9m4cSM1NTW9vu3w/dDENWvWtNrX29sfHx/P3LlzGTp0KCqVipKSErZs2UJJSYnf2t7nE6wgCIK/9Lk+WEEQhO4iEqwgCIKfiAQrCILgJyLBCoIg+IlIsIIgCH4iEqwgCIKfiAQr9HgjRozgjTfeoF+/foEOxW327NmsWLHiqs4NCQnhj3/8I0OHDvVxVG3f66WXXnJXlRJ8SyRYoUNTp07lT3/6EyEhIR7bn3nmmXa3L1q0qDtDDKh+/frxxhtveBQNulZ33nkn58+fdxdB96fGxkaysrKYO3eu3+/VF4kEK3To5MmTqNVqjwSi1+tJSEigrq6uze0nT54MRKi9gkaj4ZZbbmHv3r3dds+cnBxGjBhBQkJCt92zr+hztQiErqmuruby5cukpaW5qy+lpaVRVlbGqVOnWm2XJImCggIAYmJimDt3LsOGDSMsLIyKigp27txJbm4u0FTE+IEHHuCpp57Cbre773nnnXdy22238dvf/hZFUejfvz9z5swhLS0NRVE4f/48mzdvdhfKbktn50yZMoVHHnmENWvWMG/ePAYOHEhZWRkbNmzwqPuZlpbGgw8+SFxcHJcuXeL999/nV7/6FW+99RaHDh3ij3/8I4B7qnVVVRW/+93v3OePGzeOOXPmEBUVxblz53j33XeprKxsN+4xY8YgyzLffvute1tMTAwrV67kzTffJD09nbS0NOrq6ti2bZv7ewnwxhtvsHHjRlJSUrj++usxm81s3ryZb7/9lvnz5zN27FiMRiMffvghR48edZ9nNBo5ffo0kydP5uOPP243NqHrxBOs0KmTJ08ycuRI99cjR46koKCAgoKCVtvLy8upra0Fmvr3Tp48yWuvvcaKFSvYu3cvCxcudD/1Hj58GLVazbhx4zzul56eTm5uLoqiYDAY+PWvf43RaOTll19m1apVXLp0iWXLlhEeHt5mvN6e01yXYNOmTaxcuZKGhgYWL16MJDX9t+jXrx9Lly7l7Nmz/OEPf+CDDz5o9VF65cqVQFNye+qpp3jxxRfd+yIjI5k+fTpvvfUWq1atQqfTsXDhwg6/1yNGjKCkpMSj0lWzzMxMcnNzef7558nLy2PhwoXExcV5HDNr1iyOHz/OCy+8wLFjx3j00Uf56U9/yokTJ1i5ciXHjx/n0UcfRa/Xe5x37tw5j5+l4BsiwQqdKigoICEhAYPBAHyfYAsLC4mPj/fY3rJ7oKysjOzsbC5cuEBlZSV79uzh+PHjTJ48GWgqH/j111+TkZHhPicxMZFBgwZx8OBBAG655Raqqqp47733KCsrcz9FWiwW93Wu5O05kiSxadMmioqKuHTpEp9++in9+/enf//+7usYjUbWr19PeXk5BQUFbN261eNeRqMRaCrAXF9fj8lkcu/TaDS8/fbbFBcXc+HCBXbu3ElKSgoaTfsfHGNiYty/oK60Z88e8vLyqKio4JNPPsFut7dKiocPH+bgwYNUVFTw6aefotVqqaioICcnx70tJCSE5ORkj/NqamqIjY1tNy7h6oguAqFTzUlz5MiRFBUVERsbS2FhIVarlQsXLri3DxgwgA8//NB9nizL3HPPPYwdO5bIyEjUajWyLLu7EAByc3P513/9VyIiIqivrycjI4Pz589TXl4OwNChQxk6dGirFVG1Wm2rp7dm3p7jcrkoLS11f92c2AwGA5cuXSI+Pp5z5855lKQ7c+aM19+35vWbWl5fkiQMBkO7dURlWfaonN9Sy1hdLhdGo9H9y62tY0wmE06nkwsXLri3NTQ0YLfbW51nt9uRZdnrtgneEQlW6JTZbKa0tJS0tDQ0Gg0lJSXuSu4FBQXu7U6n0yN5PvDAA4wbN47Nmzdz8eJFbDYbc+fO9agMn5+fj8lkYvLkyezevZuJEyd6LESpUqk4efJkq/XTgHYTkbfnKIrikTxblmz0hStL2XlzfZPJ1Orje0fXu/JaTqez1XltbbtysT69Xu/xy0DwDZFgBa8UFBQwfvx41Gq1RzdAYWEh8+bNQ61Wc/78eY8lNFJTUzl06BB5eXlA03/quLg498dqaEoShw4dIiMjg4sXL6LX6/nqq6/c+8+fP8+UKVOora31eBHWkas5py3l5eVMmjQJlUrlTo5XfrRuTl6+Wlm1uLiYGTNm+ORaXTFo0KBuGRbW14g+WMErJ0+eJDY2lhtuuMHjKfXUqVNERUW12g5w6dIlxo0b514iesGCBW1OBjh48CCDBw/mvvvucz/RNtuzZw+SJLFkyRJSUlKIiYlh+PDh3H///a2S3bWc05bs7GwiIiKYP38+AwcOZMSIEdx///0ex5hMJqxWK6NHjyYiIsK9VtPVys/Pp3///kRFRV3TdboqJSWF48ePd+s9+wLxBCt45dSpUzgcDmRZpqioyL3darVSXFzMsGHDOHHihMc5H3zwAY888gi//OUvsVqt7N27lyNHjrhfIjW7cOECxcXFDBkyhO3bt3vsMxqNrFq1iszMTPf6afX19RQVFbmXnL7S1ZzTltraWl5//XUefPBBMjIyuHTpEh9++CFPPPGE+8lYURTee+897r33Xm6//XZqa2s9hml11cWLFykoKCA9Pd2jq8SfRowYQUhICIcPH+6W+/UlYkUDQeiClJQUfv3rX7NixYoOx+Fe6z1++tOf8vvf//6auji89fjjj1NYWMjOnTv9fq++Rh0dHf1coIMQhGA1ffp01Go1kiQxfPhwHn74YS5fvuzXp8vq6mqsViuNjY1+f/EUEhJCVFQUu3btanPsrXBtxBOsIHQgMzOTyZMnYzAYqK+v58SJE3z00UeYzeZAhyb0ACLBCoIg+IkYRSAIguAnIsEKgiD4iUiwgiAIfiISrCAIgp+IBCsIguAn/w/wWpvD18dfeQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import qxplots as qxp\n", "import makeqx as mkq\n", "import matplotlib.pyplot as plt\n", "#import jtplot submodule from jupyterthemes\n", "from jupyterthemes import jtplot\n", "# currently installed theme will be used to\n", "# set plot style if no arguments provided\n", "jtplot.style()\n", "tfi200 = mkq.digitize_qx(b200[1:], dlevels=2)\n", "qxp.plot_qx(tfi200)\n", "Rs = mkq.tmm_eval_wsweep(qx=tfi200, cthick=b200[0], inc_ang=0, lam_low=400, lam_high=800, lam_pts=100, n_subs=1.52)\n", "lams = np.linspace(400, 800, 100, endpoint=True)\n", "plt.figure(figsize=(5,2))\n", "#/plt.plot(lams, 100*Rs0)\n", "plt.plot(lams, Rs)\n", "plt.ylim(0,1)\n", "plt.axhspan(1, 4)\n", "plt.xlabel('Wavelength (nm)')\n", "plt.ylabel('Reflectance (%)')\n" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "ExecuteTime": { "end_time": "2018-11-04T08:07:43.891632Z", "start_time": "2018-11-04T08:07:43.887885Z" } }, "outputs": [ { "data": { "text/plain": [ "(265.5202634234206, 0.3274814296523358)" ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" } ], "source": [ "b200[0], np.mean(Rs)" ] } ], "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.7.0" }, "latex_envs": { "LaTeX_envs_menu_present": true, "autoclose": false, "autocomplete": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false } }, "nbformat": 4, "nbformat_minor": 2 }