Skip to content

Instantly share code, notes, and snippets.

@astropenguin
Last active February 7, 2022 09:17
Show Gist options
  • Save astropenguin/f9dd3df5854cb668e4836ff5bf811ef9 to your computer and use it in GitHub Desktop.
Save astropenguin/f9dd3df5854cb668e4836ff5bf811ef9 to your computer and use it in GitHub Desktop.
casa_atm_model.ipynb
Display the source blob
Display the rendered blob
Raw
{
"nbformat": 4,
"nbformat_minor": 0,
"metadata": {
"colab": {
"name": "casa_atm_model.ipynb",
"provenance": [],
"collapsed_sections": [],
"authorship_tag": "ABX9TyPU3JOKW1oVELB+RUw5hCa7",
"include_colab_link": true
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3"
},
"language_info": {
"name": "python"
}
},
"cells": [
{
"cell_type": "markdown",
"metadata": {
"id": "view-in-github",
"colab_type": "text"
},
"source": [
"<a href=\"https://colab.research.google.com/gist/astropenguin/f9dd3df5854cb668e4836ff5bf811ef9/casa_atm_model.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"id": "0kl40qu8OtGO"
},
"outputs": [],
"source": [
"! pip install -q casadata casatools jupyter-io"
]
},
{
"cell_type": "code",
"source": [
"# standard library\n",
"from typing import Any, Tuple\n",
"\n",
"\n",
"# dependencies\n",
"import casatools\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from jupyter_io import savefig_in_notebook"
],
"metadata": {
"id": "nvoo31LNPrhQ"
},
"execution_count": 9,
"outputs": []
},
{
"cell_type": "code",
"source": [
"def get_tau(\n",
" f_min: float,\n",
" f_max: float,\n",
" f_step: float,\n",
" pwv: float,\n",
" **atm_params: Any,\n",
") -> Tuple[np.ndarray, np.ndarray]:\n",
" \"\"\"Compute of zenith opacities at given frequencies.\n",
"\n",
" Args:\n",
" f_min: Minimum frequency (in units of Hz).\n",
" f_max: Maximum frequency (in units of Hz).\n",
" f_step: Frequency step (in units of Hz).\n",
" pwv: Precipitable water vapor (in units of mm).\n",
" atm_params: Parameters fo the ATM model.\n",
"\n",
" Returns:\n",
" freq: Array of frequencies (in units of Hz).\n",
" tau: Array of zenith opacities at the frequencies.\n",
" \n",
" \"\"\"\n",
" at = casatools.atmosphere()\n",
" qa = casatools.quanta()\n",
"\n",
" f_cent = qa.quantity((f_min + f_max) / 2, \"Hz\")\n",
" f_width = qa.quantity(f_max - f_min + f_step, \"Hz\")\n",
" f_step = qa.quantity(f_step, \"Hz\")\n",
" pwv = qa.quantity(pwv, \"mm\")\n",
"\n",
" at.initAtmProfile(**atm_params)\n",
" at.initSpectralWindow(1, f_cent, f_width, f_step)\n",
" at.setUserWH2O(pwv)\n",
"\n",
" freq = qa.convert(at.getSpectralWindow(), \"Hz\")[\"value\"]\n",
" tau = at.getDryOpacitySpec()[1] + at.getWetOpacitySpec()[1]['value']\n",
"\n",
" return freq, tau"
],
"metadata": {
"id": "_nh3Ww1COvBy"
},
"execution_count": 10,
"outputs": []
},
{
"cell_type": "code",
"source": [
"atm_params = {\n",
" 'atmType': 1,\n",
" 'humidity': 20,\n",
" 'temperature': '270 K',\n",
" 'altitude': '5000 m',\n",
" 'pressure': '560 mbar',\n",
" 'h0': '2 km',\n",
"}\n",
"\n",
"freq, tau = get_tau(220e9, 440e9, 0.5e9, 1.0, **atm_params)\n",
"plt.plot(freq * 1e-9, tau)\n",
"savefig_in_notebook()"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 298
},
"id": "NdGnPx7WP-OS",
"outputId": "3f36b136-451b-45ec-f943-f6f0f1c3f4e4"
},
"execution_count": 11,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/html": [
"<p>Download: <a download=\"figure.pdf\" href=\"data:application/pdf;base64,JVBERi0xLjQKJazcIKu6CjEgMCBvYmoKPDwgL1BhZ2VzIDIgMCBSIC9UeXBlIC9DYXRhbG9nID4+CmVuZG9iago4IDAgb2JqCjw8IC9FeHRHU3RhdGUgNCAwIFIgL0ZvbnQgMyAwIFIgL1BhdHRlcm4gNSAwIFIKL1Byb2NTZXQgWyAvUERGIC9UZXh0IC9JbWFnZUIgL0ltYWdlQyAvSW1hZ2VJIF0gL1NoYWRpbmcgNiAwIFIKL1hPYmplY3QgNyAwIFIgPj4KZW5kb2JqCjEwIDAgb2JqCjw8IC9Bbm5vdHMgWyBdIC9Db250ZW50cyA5IDAgUgovR3JvdXAgPDwgL0NTIC9EZXZpY2VSR0IgL1MgL1RyYW5zcGFyZW5jeSAvVHlwZSAvR3JvdXAgPj4KL01lZGlhQm94IFsgMCAwIDQzMiAyODggXSAvUGFyZW50IDIgMCBSIC9SZXNvdXJjZXMgOCAwIFIgL1R5cGUgL1BhZ2UgPj4KZW5kb2JqCjkgMCBvYmoKPDwgL0ZpbHRlciAvRmxhdGVEZWNvZGUgL0xlbmd0aCAxMSAwIFIgPj4Kc3RyZWFtCnichZhLrxzHDYX38yt6mWzaVawXaxkjiYDs7AjIIsjKr8SwAzgG4r+f75A10s2VZQkQ7syZahbJIs9hdb2+f9Tru6tc3/P/l+vv1z/4+/VVrzf8/+5R+Pbjozfj7w/x19z5VM7ffz4e3z4++wNLf2bFm8foV5s80dxv18cfzkcb7e6dr6x49/k8bfn0d+yIFyx+7wc29cuj1nIvs1aP/Zff7R7YenxOIL88Pn97ffbnetVyvf32UUu9qzcbdfQ9L6t3L9PG9fbrx+9slN9fb7+//vT2EVs+6tq37+HTzxYvvv/GFqvcbUwvw6a/2qKVV1tY33eZu1vNLV5+//gW1vY9WmnDt6/XW7yOolW/rdte4xzEi+8f36IVv1dpu9oq89UW/XUUzdfd6xy2n2f9/vtvbLHW7Wv1ObuN11u8joIi6ePuy/rY7DAKZt8Dv75BL/dMk73ec43hvbUd9j+07vPepbfRntbfA5+y7u2uZfXltfUa5scH5qv53dq0vp72XyAfyU+/zatpZee4d1t127bYoH4YQJ37XtMrLjx3eI98aoc57r2qzUE+W+7wYQxW622lUHfPHV4gn9jBCiVRWh171N2z3z6MwYbdsxGmv9vhPfKpHTqsQvWPBgc8O/q5w09X8lBjOdRTl+jmP99cf7v+fdn1l6uy1S8wTbVK606jOcpNZea/hYerDHdrqukv31yf/fGb//7rq2++fPP59dXPjweJturVTTXp3ma0wZq327LVhG49q4Og6klBKX71eZeyygrUqTZqYcdazgNi/eHh/e4yFsW/56gl0CUaKhxrWthV6YdUoaa1llCz1UegeLZ81PBsb+9VNLvtLn1oRdjl30h0qchkoJbmFks7SYECaiwdBLkSJeltr0CXe5odt9zCdblQRvdAsZasmS7MrYYVE7NHqSujGLW2Ay+RU9jo1mbxgI0vtfvJ2rRZAx7QOW2dcfDgDHTdNuMQcsfm9cB7jLbCdNlwWzqy6cM0Eo6MiEV6UloZK9KhDf2gR2VAB6ajfSt9QGbLiLX8GoZhtKWDTsNE0sINgy5WGyNMlNVqJJ8ShpZF50JteC6GuTxymodSl0WWOrSqfIRzMCBSEHC7LWKJ1WQmIwTeUKwlClHl4n7PcCRsTFfnCx53K7vsmqs3PwimGz1VQ7B5n7F60DnkoMbqVuus4TZdazZ7j2hqn0qj4EUmHQ4ImHLuEfrY6lxOJmBEzcP2hABO/8jB0Xv4PeGcdw0EZY08mQlRnA6SkTlKeDKJ/bmj2TnySehrSZ5BvS7aNGBCT/HVYtIXTSgapWPxMG1MikEwOr9LayWMNHQobS9Cf9+H5CYrlbYna9vmSVTJYKSBp72IEYIpaYRmPv0l26Xmaif0tSVigudqWWhO6BMWjJokIzPRxhyiLs/Fu3nCnd5/roYCIKCAEdN3fcdz5EHwxpGlQzj1PlwnhkP0DMQRjqy5CSFgi+rPPthIu1nAnUaIVulRADOC5JlbpFldcDNorB7YT39AlTD2gfnCSEEQmpHQRE8jTDPEIJF3Mtwt6tJoatRZM55WUxfR1Eb7tqb+FTyc2CMc2reQ7bkuypxTmFGugleV8YsydxR+RZQVLkLvmyMqN7LAaQXcglRWvyjzbj56wj0scgxzckyjR1GhM7CAbKjoKoy2DuoUQ50Xxc8B535kfojEr6GKkr4FDJ2JbMO7ILDIE+PdpCohBGKpg1E0LW/VKMWjyNE3JifBBsvlpEb6JuNYZE+zmLK/cSUyYUiDDrDGmYxO/gNGHKY6Jk+QDObj6MAhLWDC7BGHTQ2eQVo0Q3H3hHVAT8JB5FYUtTE5djVplJhThDu2bLi5xXEB92x+Yxhb74iF3mf6FqzTVMmEI1QDJRkw7GnOYYV/KHtQhS4hu3Le4R6tZeHegJeP2i4R9/aD7tPli9kAbyPLHNA8crtEFW7pB8P0SGlVQujR9GMzhqQMLynnzsg5ej9UoWTbDpZE8fA6qUKkgTTmatOcGqLLah8jHYTi9qEKHOwSkYDJwyEFHe9oobqC1yEFv589NPl4FJoma8+2hQ7rYQpqimFpjgP7UWieY9IISbKlCbHvEQW2ib0+YRd/9DCyULCwDXmOkGJOrqBYMViiAFCqyjQqHTJcT5ihdXN49CHzRcuSXpB7qj8jUef3mTB9luo/9GnqziEYtaWuIUGqghtezcoGXjkUaMu6d4v0Qdc9hwJgJMfOalFPTAVDDAgLzAN7sCGLJxNVEgK3SKZUqiQW99VDNwiQMWlpRoyRd6lqE3aXel+UZ+N4R+yIEgxVsCn2Jq2K+kMJLOmaIKknj5FT8JYaNzHIxnI2JFLQRAn9ora6sSAcRApKkjtFBEcVO+hKbo9plwSm2zGlEfmK4a9lKzGiMmhommEPZplpfmBXBdu1RdtVAic4Ckb6AN3cwReRbYbXmgOU5gxGxJxpbTMYpEIws9AhMSELr5oHJBE0Bw5YSgRz7QxfNPLAXSsbZ5PPlIgqMqabs9qYbUtqxNbYFieV8EqNcL2s0HAVMHlOjaDu6LcTkobe1IgZ+pPJAl0pEZHZ7ln1m6raokCdw0STLMOhqlI5IAq2LimEwJ7qoysVbJHuUVWpHKIm6D2YvBWGFLqP6hC7wZM9ULY50oFowfrBGa3wJTWCdO/BhTBXo2BFjgiG+dpM0zh1WH7r6j6iYNEy6U/MoBwfAhSHwMRDmcYMqjccI4mxIaXrcL96ZY1QPEgJg9oyxIrROUiDKkOFaJegmFlK3WlaLClfRTyUwkz03WsOLd6W463gfQZWrc4bAw0kDdpJdZR53PqFLtFOREj7liBA0na3I0EeolpXwBQRM34LyxVtqLma3c+4Kri3PJhWYixZwYut6ZZ54JhAc5zZOSDr1lqPjimquqJRW6PlzvUQG9xtMqlNU2rKm3Rn5xgreJ8xNla3EMPWFv6F6imn2497LDnDLbAG6ye8jxgyGBYfMXM0bilcZBnJcrX3mJ8azVmPvrnk0oMFmrrw6FtevUJAWlx8c+bdfOrnHBH0mfo29GJtWEt4SZZ1qwTWhzUOrHdiJUrY1UpxvnIr1ZCCr4hAT1jaoz4NTYBX0m+JT4ih+g4SDV4QulI6UUtKPOeDxj49pRN2I4E7dK9JTVI8dTPh0C2LFdxTPXVToPf2wSGylE9Yb5H4LHnY3dSo7UKN4S6AA+9UVYaF4hxfVA9RjFTVoWsXV7pcTUGkljUmKbc6DrqPqpr4dHvaIIqjqpQ//bYTpgOOfJJvsnZShbNHPotkN18Cc9fV6HCEfOo2G0agt3rkk7N0jjU8IUEjdVLjebXQcQqbESEFUe9iqK5cPG47yqehr+StocVwnRIHSe1WDryk9ZR9RKP7aaAu4ZXw6aK6Uftx4JXCpyphVm9R9DB4T+Uj7YMWjnmfuRHRD82H57h0ZauCeiofQx/Td46c5O2uKXFTw+XIaVawp8LpDVxn9k/TFEPejMgY1JUz/Efgx18fX1w/xWus8sHr8v9/UfV8Gf/ylTuPv3s3/+Prd/P8+Csv8BM9i37tmS8e/wOi5NOJCmVuZHN0cmVhbQplbmRvYmoKMTEgMCBvYmoKMjY2MAplbmRvYmoKMTYgMCBvYmoKPDwgL0ZpbHRlciAvRmxhdGVEZWNvZGUgL0xlbmd0aCAyNDcgPj4Kc3RyZWFtCnicTVG7bUQxDOvfFFzgAOtreZ4LUl32b0PJCJDCIKEvKaclFvbGSwzhB1sPvuSRVUN/Hj8x7DMsPcnk1D/muclUFL4VqpuYUBdi4f1oBLwWdC8iK8oH349lDHPO9+CjEJdgJjRgrG9JJhfVvDNkwomhjsNBm1QYd00ULK4VzTPI7VY3sjqzIGx4JRPixgBEBNkXkM1go4yxlZDFch6oCpIFWmDX6RtRi4IrlNYJdKLWxLrM4Kvn9nY3Qy/y4Ki6eH0M60uwwuileyx8rkIfzPRMO3dJI73wphMRZg8FUpmdkZU6PWJ9t0D/n2Ur+PvJz/P9CxUoXCoKZW5kc3RyZWFtCmVuZG9iagoxNyAwIG9iago8PCAvRmlsdGVyIC9GbGF0ZURlY29kZSAvTGVuZ3RoIDkwID4+CnN0cmVhbQp4nE2NQRLAIAgD77wiT1BE0P90etL/X6vUDr3ATgKJFkWC9DVqSzDuuDIVa1ApmJSXwFUwXAva7qLK/jJJTJ2G03u3A4Oy8XGD0kn79nF6AKv9egbdD9IcIlgKZW5kc3RyZWFtCmVuZG9iagoxOCAwIG9iago8PCAvRmlsdGVyIC9GbGF0ZURlY29kZSAvTGVuZ3RoIDgwID4+CnN0cmVhbQp4nEWMuw3AMAhEe6ZgBH4mZp8olbN/GyBK3HBPunu4OhIyU95hhocEngwshlPxBpmjYDW4RlKNneyjsG5fdYHmelOr9fcHKk92dnE9zcsZ9AplbmRzdHJlYW0KZW5kb2JqCjE5IDAgb2JqCjw8IC9GaWx0ZXIgL0ZsYXRlRGVjb2RlIC9MZW5ndGggMzM4ID4+CnN0cmVhbQp4nDVSOa7dQAzrfQpdIIB2zZznBal+7t+GlF8KQ7RWipqOFpVp+WUhVS2TLr/tSW2JG/L3yQqJE5JXJdqlDJFQ+TyFVL9ny7y+1pwRIEuVCpOTksclC/4Ml94uHOdjaz+PI3c9emBVjIQSAcsUE6NrWTq7w5qN/DymAT/iEXKuWLccYxVIDbpx2hXvQ/N5yBogZpiWigpdVokWfkHxoEetffdYVFgg0e0cSXCMjVCRgHaB2kgMObMWu6gv+lmUmAl07Ysi7qLAEknMnGJdOvoPPnQsqL8248uvjkr6SCtrTNp3o0lpzCKTrpdFbzdvfT24QPMuyn9ezSBBU9YoaXzQqp1jKJoZZYV3HJoMNMcch8wTPIczEpT0fSh+X0smuiiRPw4NoX9fHqOMnAZvAXPRn7aKAxfx2WGvHGCF0sWa5H1AKhN6YPr/1/h5/vwDHLaAVAplbmRzdHJlYW0KZW5kb2JqCjIwIDAgb2JqCjw8IC9GaWx0ZXIgL0ZsYXRlRGVjb2RlIC9MZW5ndGggMjQ4ID4+CnN0cmVhbQp4nC1ROZIDQQjL5xV6QnPT77HLkff/6QrKAYOGQyA6LXFQxk8Qlive8shVtOHvmRjBd8Gh38p1GxY5EBVI0hhUTahdvB69B3YcZgLzpDUsgxnrAz9jCjd6cXhMxtntdRk1BHvXa09mUDIrF3HJxAVTddjImcNPpowL7VzPDci5EdZlGKSblcaMhCNNIVJIoeomqTNBkASjq1GjjRzFfunLI51hVSNqDPtcS9vXcxPOGjQ7Fqs8OaVHV5zLycULKwf9vM3ARVQaqzwQEnC/20P9nOzkN97SubPF9Phec7K8MBVY8ea1G5BNtfg3L+L4PePr+fwDqKVbFgplbmRzdHJlYW0KZW5kb2JqCjIxIDAgb2JqCjw8IC9GaWx0ZXIgL0ZsYXRlRGVjb2RlIC9MZW5ndGggMjEwID4+CnN0cmVhbQp4nDVQyw1DMQi7ZwoWqBQCgWSeVr11/2tt0DthEf9CWMiUCHmpyc4p6Us+OkwPti6/sSILrXUl7MqaIJ4r76GZsrHR2OJgcBomXoAWN2DoaY0aNXThgqYulUKBxSXwmXx1e+i+Txl4ahlydgQRQ8lgCWq6Fk1YtDyfkE4B4v9+w+4t5KGS88qeG/kbnO3wO7Nu4SdqdiLRchUy1LM0xxgIE0UePHlFpnDis9Z31TQS1GYLTpYBrk4/jA4AYCJeWYDsrkQ5S9KOpZ9vvMf3D0AAU7QKZW5kc3RyZWFtCmVuZG9iagoxNCAwIG9iago8PCAvQmFzZUZvbnQgL0RlamFWdVNhbnMgL0NoYXJQcm9jcyAxNSAwIFIKL0VuY29kaW5nIDw8IC9EaWZmZXJlbmNlcyBbIDQ4IC96ZXJvIC9vbmUgL3R3byAvdGhyZWUgL2ZvdXIgL2ZpdmUgXSAvVHlwZSAvRW5jb2RpbmcKPj4KL0ZpcnN0Q2hhciAwIC9Gb250QkJveCBbIC0xMDIxIC00NjMgMTc5NCAxMjMzIF0gL0ZvbnREZXNjcmlwdG9yIDEzIDAgUgovRm9udE1hdHJpeCBbIDAuMDAxIDAgMCAwLjAwMSAwIDAgXSAvTGFzdENoYXIgMjU1IC9OYW1lIC9EZWphVnVTYW5zCi9TdWJ0eXBlIC9UeXBlMyAvVHlwZSAvRm9udCAvV2lkdGhzIDEyIDAgUiA+PgplbmRvYmoKMTMgMCBvYmoKPDwgL0FzY2VudCA5MjkgL0NhcEhlaWdodCAwIC9EZXNjZW50IC0yMzYgL0ZsYWdzIDMyCi9Gb250QkJveCBbIC0xMDIxIC00NjMgMTc5NCAxMjMzIF0gL0ZvbnROYW1lIC9EZWphVnVTYW5zIC9JdGFsaWNBbmdsZSAwCi9NYXhXaWR0aCAxMzQyIC9TdGVtViAwIC9UeXBlIC9Gb250RGVzY3JpcHRvciAvWEhlaWdodCAwID4+CmVuZG9iagoxMiAwIG9iagpbIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwCjYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgNjAwIDYwMCA2MDAgMzE4IDQwMSA0NjAgODM4IDYzNgo5NTAgNzgwIDI3NSAzOTAgMzkwIDUwMCA4MzggMzE4IDM2MSAzMTggMzM3IDYzNiA2MzYgNjM2IDYzNiA2MzYgNjM2IDYzNiA2MzYKNjM2IDYzNiAzMzcgMzM3IDgzOCA4MzggODM4IDUzMSAxMDAwIDY4NCA2ODYgNjk4IDc3MCA2MzIgNTc1IDc3NSA3NTIgMjk1CjI5NSA2NTYgNTU3IDg2MyA3NDggNzg3IDYwMyA3ODcgNjk1IDYzNSA2MTEgNzMyIDY4NCA5ODkgNjg1IDYxMSA2ODUgMzkwIDMzNwozOTAgODM4IDUwMCA1MDAgNjEzIDYzNSA1NTAgNjM1IDYxNSAzNTIgNjM1IDYzNCAyNzggMjc4IDU3OSAyNzggOTc0IDYzNCA2MTIKNjM1IDYzNSA0MTEgNTIxIDM5MiA2MzQgNTkyIDgxOCA1OTIgNTkyIDUyNSA2MzYgMzM3IDYzNiA4MzggNjAwIDYzNiA2MDAgMzE4CjM1MiA1MTggMTAwMCA1MDAgNTAwIDUwMCAxMzQyIDYzNSA0MDAgMTA3MCA2MDAgNjg1IDYwMCA2MDAgMzE4IDMxOCA1MTggNTE4CjU5MCA1MDAgMTAwMCA1MDAgMTAwMCA1MjEgNDAwIDEwMjMgNjAwIDUyNSA2MTEgMzE4IDQwMSA2MzYgNjM2IDYzNiA2MzYgMzM3CjUwMCA1MDAgMTAwMCA0NzEgNjEyIDgzOCAzNjEgMTAwMCA1MDAgNTAwIDgzOCA0MDEgNDAxIDUwMCA2MzYgNjM2IDMxOCA1MDAKNDAxIDQ3MSA2MTIgOTY5IDk2OSA5NjkgNTMxIDY4NCA2ODQgNjg0IDY4NCA2ODQgNjg0IDk3NCA2OTggNjMyIDYzMiA2MzIgNjMyCjI5NSAyOTUgMjk1IDI5NSA3NzUgNzQ4IDc4NyA3ODcgNzg3IDc4NyA3ODcgODM4IDc4NyA3MzIgNzMyIDczMiA3MzIgNjExIDYwNQo2MzAgNjEzIDYxMyA2MTMgNjEzIDYxMyA2MTMgOTgyIDU1MCA2MTUgNjE1IDYxNSA2MTUgMjc4IDI3OCAyNzggMjc4IDYxMiA2MzQKNjEyIDYxMiA2MTIgNjEyIDYxMiA4MzggNjEyIDYzNCA2MzQgNjM0IDYzNCA1OTIgNjM1IDU5MiBdCmVuZG9iagoxNSAwIG9iago8PCAvZml2ZSAxNiAwIFIgL2ZvdXIgMTcgMCBSIC9vbmUgMTggMCBSIC90aHJlZSAxOSAwIFIgL3R3byAyMCAwIFIKL3plcm8gMjEgMCBSID4+CmVuZG9iagozIDAgb2JqCjw8IC9GMSAxNCAwIFIgPj4KZW5kb2JqCjQgMCBvYmoKPDwgL0ExIDw8IC9DQSAwIC9UeXBlIC9FeHRHU3RhdGUgL2NhIDEgPj4KL0EyIDw8IC9DQSAxIC9UeXBlIC9FeHRHU3RhdGUgL2NhIDEgPj4gPj4KZW5kb2JqCjUgMCBvYmoKPDwgPj4KZW5kb2JqCjYgMCBvYmoKPDwgPj4KZW5kb2JqCjcgMCBvYmoKPDwgPj4KZW5kb2JqCjIgMCBvYmoKPDwgL0NvdW50IDEgL0tpZHMgWyAxMCAwIFIgXSAvVHlwZSAvUGFnZXMgPj4KZW5kb2JqCjIyIDAgb2JqCjw8IC9DcmVhdGlvbkRhdGUgKEQ6MjAyMjAyMDcwOTE1NTFaKQovQ3JlYXRvciAobWF0cGxvdGxpYiAzLjIuMiwgaHR0cDovL21hdHBsb3RsaWIub3JnKQovUHJvZHVjZXIgKG1hdHBsb3RsaWIgcGRmIGJhY2tlbmQgMy4yLjIpID4+CmVuZG9iagp4cmVmCjAgMjMKMDAwMDAwMDAwMCA2NTUzNSBmIAowMDAwMDAwMDE2IDAwMDAwIG4gCjAwMDAwMDY2NjcgMDAwMDAgbiAKMDAwMDAwNjQ3MyAwMDAwMCBuIAowMDAwMDA2NTA1IDAwMDAwIG4gCjAwMDAwMDY2MDQgMDAwMDAgbiAKMDAwMDAwNjYyNSAwMDAwMCBuIAowMDAwMDA2NjQ2IDAwMDAwIG4gCjAwMDAwMDAwNjUgMDAwMDAgbiAKMDAwMDAwMDM4NSAwMDAwMCBuIAowMDAwMDAwMjA4IDAwMDAwIG4gCjAwMDAwMDMxMjAgMDAwMDAgbiAKMDAwMDAwNTMyMSAwMDAwMCBuIAowMDAwMDA1MTIxIDAwMDAwIG4gCjAwMDAwMDQ3OTAgMDAwMDAgbiAKMDAwMDAwNjM3NCAwMDAwMCBuIAowMDAwMDAzMTQxIDAwMDAwIG4gCjAwMDAwMDM0NjEgMDAwMDAgbiAKMDAwMDAwMzYyMyAwMDAwMCBuIAowMDAwMDAzNzc1IDAwMDAwIG4gCjAwMDAwMDQxODYgMDAwMDAgbiAKMDAwMDAwNDUwNyAwMDAwMCBuIAowMDAwMDA2NzI3IDAwMDAwIG4gCnRyYWlsZXIKPDwgL0luZm8gMjIgMCBSIC9Sb290IDEgMCBSIC9TaXplIDIzID4+CnN0YXJ0eHJlZgo2ODc1CiUlRU9GCg==\" target=\"_blank\">figure.pdf</a></p>"
],
"text/plain": [
"<IPython.core.display.HTML object>"
]
},
"metadata": {},
"execution_count": 11
},
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD8CAYAAACSCdTiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deZRcZ3km8Oe999bSq7rVarXWliwjS8hgSXZHxrFJiB2DgWFsh2XMMB6H8RzBGTiHTAgE+CMwIQlMhiVwZkJiwmLAgThgBodowI4XDAdjW/IiW5K1WJZsrd0t9b5U1b33nT/uvdVVXb1XdVXXV8/vnD5da/dXV6Wn3n7vd78rqgoiIjKTVekBEBHR4mHIExEZjCFPRGQwhjwRkcEY8kREBmPIExEZbNaQF5H1IvKIiBwUkQMi8pHw9s+IyGkReTb8etviD5eIiOZDZpsnLyKrAaxW1adFpAnAPgC3AHgPgGFV/cLiD5OIiBbCme0BqnoWwNnw8pCIHAKwdrEHRkRExZu1ks97sMhGAI8BeB2APwbwhwAGAewF8FFV7Zvp+StWrNCNGzcubKRERDVq3759varavpDnzjnkRaQRwC8A/KWq3iciHQB6ASiAzyJo6fyXKZ63G8BuAOjs7Lzq5MmTCxknEVHNEpF9qtq1kOfOaXaNiMQA/AjAPap6HwCo6nlV9VTVB/B1ALumeq6q3qWqXara1d6+oA8iIiJaoLnMrhEA3wBwSFW/lHP76pyH3QrghdIPj4iIijHrjlcA1wK4HcDzIvJseNunALxXRHYgaNecAPCBRRkhEREt2Fxm1/wKgExx157SD4eIiEqJR7wSERmMIU9EZDCGPBGRwRjyRFTA9xX3PvUqMp5f6aFQkRjyRFTgX/afwcd/tB9fe/SlSg+FisSQJ6ICfngk/Es9wxUeCRWLIU9EBVrq4gCAiyPpCo+EisWQJ6IClhUcGtM3ypCvdgx5Iirg+cEO176RTIVHQsViyBNRAdcLevJs11Q/hjwRFfD8IOTHMl6FR0LFYsgTUQFvHicToqWNIU9EBaJKnqofQ56ICkQ9eap+DHkiKpBbyfus6qsaQ56ICrg5wc7+fHVjyBNRgWiefHCZIV/NGPJEVCA32F2GfFVjyBNRgbx2DXfCVjWGPBEVyK/kuaZ8NWPIE1GBvEqe7ZqqxpAnogLsyZuDIU9EBVjJm4MhT0QFfIa8MRjyRFTAZbvGGAx5IirAg6HMwZAnogIup1AagyFPRAU89uSNwZAnogKcQmkOhjwRFWAlbw6GPBEV4Dx5czDkiagAK3lzMOSJqADnyZtj1pAXkfUi8oiIHBSRAyLykfD25SLyoIgcDb+3Lv5wiagc8ufJcwplNZtLJe8C+KiqbgPwBgAfEpFtAD4B4CFV3QzgofA6ERkgb3YN15OvarOGvKqeVdWnw8tDAA4BWAvgZgB3hw+7G8AtizVIIiovz1fYlmQvU/WaV09eRDYC2AngCQAdqno2vOscgI6SjoyIKsb1FQnHyl6m6jXnkBeRRgA/AvBHqjqYe5+qKoAp3wkisltE9orI3p6enqIGS0Tl4eWEvK8M+Wo2p5AXkRiCgL9HVe8Lbz4vIqvD+1cD6J7quap6l6p2qWpXe3t7KcZMRIvM9RTxqJJnT76qzWV2jQD4BoBDqvqlnLvuB3BHePkOAD8p/fCIqBKCSt7OXqbq5czhMdcCuB3A8yLybHjbpwB8HsC9InIngJMA3rM4QySicvOUPXlTzBryqvorADLN3TeUdjhEtBS4viIRC0Ke8+SrG494JaICnu8jbrOSNwFDnogKuB578qZgyBNRAc/PmV3DkK9qDHkiKpA7T56VfHVjyBNRAU8ViRjbNSZgyBNRAddTxOxgUh3bNdWNIU9EBTxfEbMsOJZwCmWVY8gTUQHXV9i2wLKElXyVY8gTUQHP92GLBJU8166pagx5IioQrSdvs5Kvegx5Iirg+QrHCip5LjVc3RjyRFQg6snblsVKvsox5ImoQG4lz558dWPIE1EeVQ0qectiT94ADHkiyhNlui0Cx+Y8+WrHkCeiPNEyBo4tsIWVfLVjyBNRnijkoymUXLumujHkiSiPG7ZnHM6TNwJDnojy5Fbyji3wGfJVjSFPRHncvHYN58lXO4Y8EeXJq+TZk696DHkiypOdXZPtyXMKZTVjyBNRnolK3mIlbwCGPBHlcQsqeYZ8NWPIE1Ge6AhXzpM3A0OeiPK43PFqFIY8EeVxPR7xahKGPBHliU4SEiw1zHny1Y4hT0R5XK5dYxSGPBHlmZgnH0yh5Dz56saQJ6I8uT15i2eGqnoMeSLKM3lZA/bkq9usIS8i3xSRbhF5Iee2z4jIaRF5Nvx62+IOk4jKxdP8nny0I5aq01wq+W8DuGmK27+sqjvCrz2lHRYRVYqXs558LVbyGz/xr/jSA4crPYySmTXkVfUxABfLMBYiWgLy58lbNdWTH894AICvPnyswiMpnWJ68h8Wkf1hO6e1ZCMioorKPcerY9dWJT8wlgEANCacCo+kdBYa8l8DcCmAHQDOAvjidA8Ukd0isldE9vb09Czw1xFRuUxeoKyW5sn3jwYh35ys8ZBX1fOq6qmqD+DrAHbN8Ni7VLVLVbva29sXOk4iKpMo1C2Rmpsn3z+aBgA018UqPJLSWVDIi8jqnKu3AnhhuscSUXXJPRjKEoGvqJnzvPZlK3lzQn7Wv0lE5PsA3gRghYicAvBpAG8SkR0AFMAJAB9YxDESURll58nbQSUPBNMqLUglh1UWA2NRJW9Ou2bWV6Kq753i5m8swliIaAnI68nbYcj7iphdyVGVR7+BlTyPeCWiPLknDclW8jXSrukPZ9ckDPpEY8gTUZ7sKpQSzJPPvc10USUfdKLNwJAnojxT9uRrJOQHx4OQN+n1MuSJKI83aZ48gJqZRhnNIjIo4xnyRJRv8klDALMq25l42ZA35/Uy5IkoT+48+WwlXyPr10SfZQZlPEOeiPK52SNeke3Jm1TZzkSVlTwRGc7zfdiWQCS3J29O6M3EU/bkichwrq/ZcHfCKZS10pOPXqZJyzgw5Ikoj+9rtk1Tcz157nglItPlV/K1NbvGZ0+eiEzn5Vbydo3Nk2dPnohMF1TyQTTYUmOVfPhZpqzkichUnqeww2So3XZNhQdSQgx5IsqT8f3srJpaO+KVPXkiMl7GU8SdIBocu9bmyYffDXq9DHkiyuN6fs4UytqaJx/14g0q5BnyRJQv4/mIhU15p9aOeOU8eSIyXcZTxJzJPflamUIZfWfIE5GhMp6P2OQjXmukklfOriEi07meFrZramRZg6hdw3nyRGSstOdnZ9VEYZ/xaqVdw0qeiAzn+j7iYbhHUykzNVLJ+5xCSUSmy7ha85U82zVEZKzcKZSxMOxrLeQNKuQZ8kSUL+PnhnzwPV0rIR++TE6hJCJjZVzNVvDZdo1rTujNhJU8ERnPzankbSs4z2uttWvYkyciY6XdiZAHgr58rYR89DI5u4aIjOX6E+0aIGjZ1EpPXrnUMBGZLuP5cHIq+bht1U4lz1UoichkqhosUJbXrrFqZ8drLa5CKSLfFJFuEXkh57blIvKgiBwNv7cu7jCJqByihciiBcoAIObUTk9es6tQVnYcpTSXSv7bAG6adNsnADykqpsBPBReJ6IqFy1EFi01DNRWT96rxZ68qj4G4OKkm28GcHd4+W4At5R4XERUAVGYOzmVfC315H325LM6VPVsePkcgI4SjYeIKigK8/jkSt6tlZAPvnMKZQ4N5hxNu0VEZLeI7BWRvT09PcX+OiJaRNl2TcE8eXNCbyY1ueN1GudFZDUAhN+7p3ugqt6lql2q2tXe3r7AX0dE5ZCZol1TSz15Lmsw4X4Ad4SX7wDwk9IMh4gqaap2TdypjZ68qmbDvaaWNRCR7wN4HMAWETklIncC+DyAG0XkKIDfD68TUZWL2jKONWmefE2E/MRlk9o1zmwPUNX3TnPXDSUeCxFVWBTm+csaSE0cDJUb7GzXEJGRJkK+9ip5Ly/kzUl5hjwRZWWPeJ20dk0t7HjNa9cYVMoz5IkoKxPOh3cmrUJZE5W8z3YNERkuPUW7JphdY1DqTSNq0TiWsF1DRGaaOBhqUiVfA0e8RtW7bQmXNSAiM6XCME/G7OxtMUdqoifv5+yPYCVPREYaz3gAgIRTeycNiYLdZruGiEwVVfIJJ6eSty34ataiXVOJXl7Qk6/sWEqJIU9EWVEln4zlz5MHYPxKlNkdr+H+CFOmUTLkiShrqko+at2kXK8iYyqXidk1Vt71aseQJ6KsqXry0U7YsYzZIR+1o7KVvBkZz5Anogkp10fctmDlLDVcFw9iYjxjdrtGc3ryACt5IjLQeMZDIpYfC8mwdTNueCU/uV1jSMYz5IloQsr18/rxwES7xvSQL2zXmJHyDHkiykq5Xt7MGqB2evL+pHaNx5AnItOkMn7eTldgYjplyviefFTJh+0aQ14uQ56IsoJKvkbbNTkLlAFs1xCRgcanrOTDkDd9nnxYubMnT0TGSrlewY7XuqgnnzakfzGNibVrooOhKjma0mHIE1HWeMafYsdrNE/e8Eo+DPlY2K5RVvJEZJqpKvmaadfkrCefe73aMeSJKGuqSj7q0Zt+xKs36fy2nEJJRMaZqpIXESRjlvHtGuUqlERkuqkqeSBo2Zge8lElb2d78pUcTekw5IkoK+V6SEyaJw8EM2xMD/mocI9xqWEiMpGqIuX6SDrTVfJm9+Sjdo3NefJEZKLRtAdVoD7hFNyXcMzvyRce8VrJ0ZQOQ56IAAAjKRcA0DBFyNfF7RpaoCxaatiMlGfIExEAYDgM+cZEYU++Ie5gNG16yOfPruEUSiIyykgqCPGGeGEl35hwMDzulntIZRVNmcy2awzZBcGQJyIAuZX8FCGfdDA0nin3kMoq266xObuGiAw0U0++MeFgKGV2Je9NquQNyXiGPBEFRtLTh3xT0sFwyjVmZ+RUCo54NeS1Fv5rzoOInAAwBMAD4KpqVykGRUTlF/Xkp2rXNCUdqAbTLKf6EDDB5NP/mRLypajkf09VdzDgiZaWfScvIuPNfe/hRLumcHZNYyIGYKJvb6JoNk087MlnPIY8ES1RB84M4J1fexz/6+eH5/ycKMCnnF2TDG4bMniGTdSuqQ9ff8qQpZWLDXkF8ICI7BOR3aUYEBEVr3c4DQA4dHZwzs8ZSbmoj9uwwnZFrqawRWNyJR9V7vXhXzKmnLi82Obadap6WkRWAnhQRF5U1cdyHxCG/24A6OzsLPLXEdFiGUm70/bbo0re5Lny0bINLXVxAEDKNSPki6rkVfV0+L0bwI8B7JriMXepapeqdrW3txfz64hoEQ2nvCl3ugITO2OHU+bOlY9CvbkueK2mrNWz4JAXkQYRaYouA3gzgBdKNTAiKq+h8QyakjOHvMk9eVMr+WLaNR0Afiwi0c/5R1X9WUlGRURl1zeSRkt9fMr7mpPB7JqBMYMr+YwHEWQ/6Eyp5Bcc8qp6HMD2Eo6FiCqobzSDS1Y0THlfU9KBbQn6RtNlHlX5jLs+Eo6VPXG5KZU8p1ASGWgh5yedqZK3LEFrfRwXRwwO+YyHZMzOnricUyiJaMmaz0FQ0eOHUi6WN0wd8gDQ1hDHhWHDQ94JppDGbcuYM2Ex5IkM5M6zko/aMK31sWkf09oQM7pdk3J9JMKTmCcci5U8ES1dUSU/1+VX+keDHaqtM1byCVwwvV3jBP34RMxmT56Ilq75rrsS9dpbp+nJA8DyBtN78j6SOZW8KbNrGPJEBnLDSl4KVyiYUl8Y3i0ztGuWN8TRP5rJ/mzTjGc8JGJRJW+xkieipSsT9uTn2q7pDUO+rSEx7WPaGoMqv2/UzLny466fnT6ZdGxj1q5hyBMZKDPPKvT8wDhsS9DeNH3Ir2xKBo8dHC9qbEtVKuMhGU6fDCp5tmuIaIly/fm1a84MjKGjKQF7ihUoI2tb6oLH9o8VPb6lKJVTyScci5U8ES1d0Y7XubZrzg2MY9Wy5IyPWd0S3G9qyI9nvOyBUMmYjXFW8kS0VLnznF1zbmAcq5fVzfiYtoY44o6FMwNmtmuiI14BVvJEtMRF8+S9ORwUpao4O4dKXkSwtqUOp42t5HOnUNrsyRPR0pUJe/JRb34m/aMZjGU8rJ4l5AFgTUsSp/vMC3lVxbg7UcknY1zWgIiWsKhdk55D2+Z47wgATLsCZa4NbQ14uXckez5UU2Q8hSqyIV8fd7InNi+Gqlb8uAKGPNWE/tE0vvTgEYylzfgTfDZRu2YuAfNSzzAAYFN746yP3byyEQNjGfQMp4ob4BITnbu2Ph6EfGt9HEMpF+kiD4j6wgOHse3TP6/o+44hTzXhh/tO4asPHcUXHzhc6aGURTS7Zi6rUR7vGUHMFqxvnXnHKwBsXtkEADjWPVzcAJeYiyPBh1a0Cufy7IFfxS3j8Pe/OI606+Nbv365uAEWgSFPNaE3XCJ3z/NnKzyS8nCzlfzsbZWXeoaxsa0Bjj17HLxmZVDtmxby0RLKKxqDg8HawrAvdmnlaJmIJ1++WNTPKQZDnmrC8bAlcXZw3JiFp2YSLTWcmcOO14NnBnHZqqY5/dyO5gRa62N44fRAUeNbaqLVNbOVfPi9mAXZ+kbS2eLiyLmhIke4cAx5qgkvhzsXVYFTfaMVHs3iS4eVfMaduZK/MJzC6f4xbF+3bE4/V0Sws7MVz7zSX/QYl5IL2bV74nnfL4wsfN/DsbCw2LVxOc4MjGNwvDJr/jDkyXierzh5YRS7LlkOADjRa37IZ9s1s1Ty+08FFfkV61rm/LOv7GzB0e5ho07qfTGsuFtLWMlHLa23vX4VAODo+cpU8wx5Mt7pvjGkPR/Xb10JADhxYaTCI1p8bnbH68yV/BMvX4RjCV63dm6VPAB0bQw+LH9z/MLCB7jEXBxJoTnpIBbul2ipj0NkYgnmhTjWPYxkzMKNl6/CLTvWZKdnlhtDnox3vDeoqHaub0Fz0sm2bkyWbdfMMrvmwYPn8IZNbWhMOHP+2VdtaEVTwsGjh7uLGuNS0juSzu50BQDbEiyvj6NngTtePV9xtHsYm1Y0Ym1LHf7mtp24fM3cP0hLiSFPxjveE4T6pvZGbO5owtHzZs0MmUpUyc80u+ZY9zBe6hnBmy/vmNfPjtkWfueydvzboe6KH+hTKj2DqbyQB4ANbfXZYwjm68//5QAeO9KDNS2zT0tdbAx5Mt7LvSNoSjpY0RjHZR2NONI9ZNwRm5NFvfiM70/7Wh88eB4A8PuvnV/IA8A7tq9Bz1AKvzzau/BBLhGqiiPdQ7h0Zf7BYFtWNePwuYW9V+575jQA4O1XrCrJGIvBkCfjHTk/hE3tjRARbF7ZhP7RTHZqm6lylxqebpGy//fCWVyxbtmCqs3rt65EW0Mc33n8RBGjXBq6h1LoH81gS0d+yG9d1YSBsQzOD85vho3vKzKejzuvuwS37lxXyqEuCEOejOb7igNnBnFFuGNxSzgf/MAZs+Z5T5bbi3enCPmnTlzE/lMDeOeVCwuhuGPhD397Ix453INnXulb8DiXgsPhHPYtq5rzbt8avlfme0zAucFxjGf8Oa0FVA4MeTLa8d5hDKdcvD6cB76zswVx2zKizTCT3F784KSpjqqKv/7Zi2hriOM9XesX/DvuuHYjVi9L4mM/3F/VB5g992ow53/rpAPCtq9vQX3cxqNH5reDObsPiCFPtPiig3auCEO+Pu7g6k3LjZoZMpnvK071jWLb6qAyjebCR77z+Ek8daIPH79pC+riC5/W15yM4a/fdQWOdQ/jc3sOFTXmSnrw0HnsWN+SnSMfScZsXPeaFfi3g/PbwfzAwXNwLJnzUcSLjSFPRnv4xW6sbErgspUT/+HetGUlXuoZwasXzTwo6sSFEYykPdy2az0cS/DMqxPtlAcPnsdnf3oQ129dWVQVH3nj5nbced0luPvxk/jCzw/Dn8NJSpaSl3qGsf/UAN5y+dQ7SN911TqcGxzHvXtPzennPX9qAN9/8hW8u2tdwWydSpn75FiiKjOScvHo4R6886q1sHJOUP2mLe347E+BRw934/ZrNlZugIvkhTODAICuDcuxbU0z9jx/Ds++2o/9rw5gKOXiinXL8JXbdkDmepbvWXzyrVsxknLxvx85hqdOXMSfvGULuja0luznL5aobVUft/Hurqn3Tdy4rQO7LlmOP//pAbTUx/DW162a8nX5vuLBQ+fxiR/tR3tjAn/y5i2LPfw5Y8iTsb77m5MYy3gFOxc3rWjAJSsacO/eU3jf1RvyPgBM8NCh80jGLGzuaMTH3rIFt3/jyewBYP+haz3+7B3b0DCPg59m49gWPvcHr8f29S344gOH8e6/exxbOprw7q51uHXnWrQtkYo2l+8rvvDAYfz8wHl8/KYt01bdIoL/8x+vxB3ffBL/7Z6ncdWGVnRtaMX29cG+nd7hFA6dHcQjh3vwysVRbOlowt/fftWSes1SzvnCXV1dunfv3rL9PqpdB88M4p1f+zWu3rQc337/roL7f/zMKfz3f3oOn37HNrz/2ksqMMLS6R1OYSTlYkNbA372wll88HtP44O/eyk+8datAII1U05eGMUNr1256NX1aNrFfU+fxj/vO4XnXu2HYwlueO1K3LxjLXZ2tmBlUxICYCjlYlldbFHHMpW+kTQefrEb3/71CTx/egC3/dZ6/NWtr5/1g971fNzzxCu454mTONE7mj2iGAhONHLVhla888p1ePsVq7NLI5SSiOxT1a4FPbeYkBeRmwB8BYAN4B9U9fMzPZ4hT4vp6PkhfPgfn8Gbtrbjnt+8goaEjfs/fB06mgvPXer7it3f3YuHX+zGH994Gf7rGzdVbG2RhVJVPHqkB+//1lMAgN+5rB2PHenBzs4WfPfOq+e1VMFiOHJ+CP+891X8+JnT2eMS4raF5joHvcNpvHHzCuxc34KmZAxpz0f/aBpNyRhe7h3BL4/24o5rNuC6zSvQlIwhGbOyH1Cepzg3OI66mI26uA1fFa6nGBrP4JlX+/HKxVHUx2x4qmhKOOgbzeB0/xiO9wzj5MVRqAZ/zX3o916DP7hy7bw/+EZSLk5eGIXr+2itj2P1suSc1uIvRkVCXkRsAEcA3AjgFICnALxXVQ9O9xyGPC2GtOvjeO8w/uz/HsCTJ4KTM1yzqQ1ffM/2GQ/0GUm5+NgPn8Oe589hRWMcN2ztwM7OFqxuqUPCsbBtTTOaEs6iVb/jGQ/He0ZwfnAc65fXoTERw0+ePY36hIPL1zRjcCyDo+eH0TeaxvOnB7BmWR0cW9A3mkb3YApn+sdwZmA8+/M2tNXjlh1r8cHfvbSoWTOllvF87D/VjwNnBnG6bwwnL4xiWV0MvzrWizMDY5gcQR3NCQyOuRgrYlpm3LZgWUDK9bGsLoZVzUlc2t6I165uwjWXrsCVnS1Lfp9BrkqF/DUAPqOqbwmvfxIAVPVz0z1nKYS8quJ0/xjWLKvL/omWdn3E7OByz1AK7U2JvDeA7yt+cbQH29e1ZJcgjUws6ar43m9O4vqtK+d0rkwgOBJRVWFbMuUbTjU4ubAC8LOXNfufIrruKzCW9rDvZB+2rmrCobODeP26ZWiIB5Vc9C/s+j5O9Y1h9bIkrEm/L6qGxjIeDp4ZxMYVDWhKOhgcy+DrvzyOd1yxBpvaG6GYWBNlcDwDx7JQH7fh+QrX1/C7jzP942hKOojZAtdTjKRd9A6l0VIfg4igMeHA1+DIQNcLnpfxg8uuH5z82A2PHPR8RcYLLg+NZzA07uLl3hGc6hvDeMbLO9jnfVd34vqtK3H91rm3Jp44fgHfefwkHjvSg6FJJ28WARrjDpJxGwnHQsKxEHdsSM79ud99P9yW4WvIeBPbJOOFt4X3zWciyqXtDRgYy0A1WA63vTGBjuYErljXgu3rW7BzfQtEUFXBBQT/t0bSwblUow+m+rgDz1c8d6ofA6MZDI5nkMrkTGEUoKM5idGUi4yvsEVgW4KGhI3O5fVY0ZiAbQnitgVFsNhYtSsm5Iv5e24tgFdzrp8CcHURP29aX33oKO5/7gz8MPF81fzg0+CAj/qEDceyoBoEn68afgWB6oWhkXL98D+rBdsSDIxlELOC/8BDKReJ8PaI5ytSrg/HEiRjdhjKwZ+Noxkv77Dxv/jXQ3N6U0VjjDiW5L2mpWbP8+cq9rstCXbuxW0LjQkHTUkHHc1J/PalK5CMWUjGbGxoq8f2dS3YuIADUK7e1IarN7WF88vH0D0UnODh6PlhjKRcDKc8jGU8pFwPKdfPCZyJpQOia1YYtDFb4FgWHFvgWALHthALvzu2IGZZSMYsbFzRgFXNSbxycRQp18f61nrYlmA846E+buOyjiY018WMCKqpWJagKVnYm7ctwZWdrRUYkXkWvWknIrsB7AaAzs7OBf2MlU0JbOloAgSwRCCY+M8kwe9AU9LBaNqF5wf3WSKwrOA+SwDHsmCJwLaAtsYEeodScP3gQ2BZXdATHEt76GhOom8knVOhBRca4g7GMh7Srp/98HAsCw3hB0vK9bCpvRGn+kbndF5NIFjNTwTZqtWS4MNDgl8cvjZAIOHrnRhPdLvIxOtta4xjYDSDlvp4WPVNjCN63orGBPrH0hDkh4YIwjASrGutR/9oBiMpF7YleO3qZrzUM5z9MIu2eXNdUHGNpT04tsC2LDhWUFW11McwHoahYwnq4jZa6+PoH03DsSwMpTKwrSAIY/ZECNpRIEYhaUnZZr9YlqCzrR6dbfUAgOu3zn/hroWK1mgnKrViQv40gNyjKdaFt+VR1bsA3AUE7ZqF/KLbdnXitl0L+4Cg0njNyrm1oGYzud1FRIurmF3CTwHYLCKXiEgcwG0A7i/NsIiIqBQWXMmrqisiHwbwcwRTKL+pqgdKNjIiIipaUT15Vd0DYE+JxkJERCXGBcqIiAzGkCciMhhDnojIYAx5IiKDMeSJiAxW1qWGRaQHwMmy/cLKWAHA7BOIzh+3SSFuk0LcJoWibbJBVdsX8gPKGvK1QET2LnQhIVNxmxTiNinEbVKoFNuE7RoiIoMx5ImIDMaQL727Kj2AJYjbpBC3SSFuk31oXVgAAAPWSURBVEJFbxP25ImIDMZKnojIYAz5eRCR9SLyiIgcFJEDIvKR8PbPiMhpEXk2/HpbznM+KSLHROSwiLylcqNfHCKSFJEnReS5cJv8j/D2S0TkifC1/1O4HDVEJBFePxbev7GS418MM2yTb4vIyznvkx3h7SIiXw23yX4RubKyr2DxiIgtIs+IyE/D6zX7PolMsU1K+j5hyM+PC+CjqroNwBsAfEhEtoX3fVlVd4RfewAgvO82AJcDuAnA34YnQDdJCsD1qrodwA4AN4nIGwD8TwTb5DUA+gDcGT7+TgB94e1fDh9nmum2CQB8LOd98mx421sBbA6/dgP4WtlHXD4fAXAo53otv08ik7cJUML3CUN+HlT1rKo+HV4eQvAPs3aGp9wM4AeqmlLVlwEcA7Br8UdaPhoYDq/Gwi8FcD2AH4a33w3glvDyzeF1hPffINV29ulZzLBNpnMzgO+Ez/sNgBYRWb3Y4yw3EVkH4O0A/iG8Lqjh9wlQuE1msaD3CUN+gcI/H3cCeCK86cPhn1DfFJHoDMRTnex8pg+FqhT+ufksgG4ADwJ4CUC/qrrhQ3Jfd3abhPcPAGgr74gX3+RtoqrR++Qvw/fJl0UkEd5WE+8TAH8D4OMAojOht6HG3yco3CaRkr1PGPILICKNAH4E4I9UdRDBn02XIvjT/CyAL1ZweGWnqp6q7kBwnt9dALZWeEgVN3mbiMjrAHwSwbb5LQDLAfxpBYdYViLy7wB0q+q+So9lqZhhm5T0fcKQnycRiSEI+HtU9T4AUNXz4X9qH8DXMdGSmdPJzk2hqv0AHgFwDYI/JaMzj+W+7uw2Ce9fBuBCmYdaNjnb5Kaw3aeqmgLwLdTW++RaAP9eRE4A+AGCNs1XUNvvk4JtIiLfK/X7hCE/D2FP8BsADqnql3Juz+2L3QrghfDy/QBuC2cKXIJgh8mT5RpvOYhIu4i0hJfrANyIYF/FIwDeFT7sDgA/CS/fH15HeP/DatjBGtNskxej90n4ProF+e+T/xzOnngDgAFVPVuBoS8aVf2kqq5T1Y0IJiM8rKrvQw2/T6bZJv+p1O+Tos7xWoOuBXA7gOfDfisAfArAe8NpTgrgBIAPAICqHhCRewEcRDAz50Oq6pV91ItrNYC7w1lDFoB7VfWnInIQwA9E5C8APIPgwxHh9++KyDEAFxG8uU0z3TZ5WETaAQiAZwF8MHz8HgBvQ7BjfhTA+ysw5kr5U9Tu+2Q695TyfcIjXomIDMZ2DRGRwRjyREQGY8gTERmMIU9EZDCGPBGRwRjyREQGY8gTERmMIU9EZLD/D6OoDKY9+ec4AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"plt.plot(freq, 270 * (1 - np.exp(-tau)))"
],
"metadata": {
"id": "KM2WITSncyZQ",
"outputId": "1ec571e6-cf35-4d46-a17d-0a44817c01b2",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 296
}
},
"execution_count": 12,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f750344e690>]"
]
},
"metadata": {},
"execution_count": 12
},
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEFCAYAAAAPCDf9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO2dd5xcZ3nvv+/ULbO9qDdLsmTZuArLYDAuuGAIhlBi00zCvXYI5MIFLgHnQggJCaGZS0KJAQeH2JhiijHGxBh3Y8mSi2RbvZeVtu/O7uz09/5xzpmd3Z3dnV7Oeb6fjz6aPXNm9p3R6Hee+b1PUVprBEEQBGfgqvQCBEEQhPIhoi8IguAgRPQFQRAchIi+IAiCgxDRFwRBcBCeSi8AoLOzU69cubLSyxAEQagptm3b1q+17srlMVUh+itXrmTr1q2VXoYgCEJNoZQ6nOtjxN4RBEFwECL6giAIDkJEXxAEwUGI6AuCIDgIEX1BEAQHIaIvCILgIET0BUEQHERV5OkLglDdjEXi/PK54wyOR0lqTVKDAprqPLz1vCV0BPyVXmJJGRqPcteWI2xc0cam0zoqvZyCENEXBIdxqH+clnovbY2+rM7XWvMX//EMWw4NZrz/+08c5OFPXEqd113MZVYVn/nVi9y3vQe/x8X9H3ktq7sClV5S3ojoC4LDuPQrj9DW4OW5z16V1flbDw+x5dAgf3PNev7na1fhdimUUiSTmkf29PIXP9jKr184wTs2LivxyivD4HiU+7b38Lbzl3LPs8e4f3sPf33F2kovK2/E0xcEBzIUimV97h929eJxKd7/6pV43C6UUgC4XIrL1nWzsqOB+3f0lGqpFWdwPArA69Z1Ued1MRrO/r2rRkT0BUGYk5dPjLJ2QRP1vpn2jVKK85a3sbMnWIGVlYexSByAgN9NU52XYDhe4RUVhoi+IAhz8nLPKBsWNc96/xmLmjg5GmbIjIjtxljYEn0vTXUeEX1BEOxL/1iEvmCEMxY1zXrO+oXGBWHnydFyLausWJF+oxnpi70jCIJtOToYAmBlR+Os55y+wLggHOgbL8uayo0l+k1+L80S6QuCYGdOjoQBWNRaN+s5XU1+3C5Fz8hEuZZVVsYtT7/OY9o7EukLgmBTeizRb6mf9Ry3S7GgyZ86125MsXf8spErCIKNOTkaxudx0dbgnfO8Ra319AzbV/R9bhd+j1s2cgVBsDc9I2EWtdSlcvNnY2FLHSdHbSr64TiBOqOOtanOy0QsQSyRrPCq8kdEXxCEWTk5MsHC5tn9fIvFLXWcGJ5Aa12GVZWXsUicRr9Ro9Bkiv9YDUf784q+UmqZUuphpdTLSqmXlFIfMY9/Til1XCn1vPnn2rTHfFoptU8ptVspdXUpX4AgCNmTTOYmyn3BCF1N8zdTW9BcRySeZGSitjc5MzEWiRPwG/aWJfq1bPFk03snDnxca/2sUqoJ2KaUetC871at9VfST1ZKbQCuB84EFgO/V0qdrrVOFHPhgiDkTjLHSHxwPEp7Fo3ZOgK+1PmtDdk1cqsVxsJxAmak3+g3JHM8WruiP2+kr7Xu0Vo/a94OAjuBJXM85Drgbq11RGt9ENgHXFiMxQqCUBiJHEQ/nkgyGo5nJeLtjca3gaGQ/apyx6PxlNjXeQ3JjMYL8/QTSc3uk5VpXZGTp6+UWgmcB2w2D31YKbVdKXW7UqrNPLYEOJr2sGPMfZEQhIry5L5+QjUcueVCMgetsqya+TJ3ANobrEjffvZOJJakzmNE+j638XekQNG/a/Nhrv76Yzx9YKDg9eVK1qKvlAoA9wAf1VqPAt8GVgPnAj3AV3P5xUqpm5RSW5VSW/v6+nJ5qCAUjaHxKO/+3mZu/uG2Si+lLOQS6VudOLOxd9oajQuDHfvvRBNJfB5DKv1mpB+JF+ZWnxqNAPDAiycLW1weZCX6SikvhuDfqbX+OYDW+pTWOqG1TgLfZdLCOQ6kN9Zeah6bgtb6Nq31Rq31xq6urkJegyDkjRXNPr63v8IrKQ+5ePrDplWTnb1jRvo2tHei8TTRN/+OxAqL9D1uIwX2uSNDhS0uD7LJ3lHA94GdWuuvpR1flHbaW4EXzdv3AtcrpfxKqVXAWmBL8ZYsCMUjPQujlnOvsyWX7B0r0s/G3qn3uvF7XKne83YiEk+kiX5x7B0r2OgNRgpbXB5kk71zMfBeYIdS6nnz2C3ADUqpcwENHAJuBtBav6SU+gnwMkbmz4ckc0eoVtI7Jh4fmmBl5+yNxexAIifRNwS8LYtIXylFR6PPpqKfTEX4qUi/QHvHEv3RCqS4ziv6WusnMGYgT+f+OR7zBeALBaxLEMpCevOsUNT+sUkuafqWP9+aRaQP0Nbos6enH8/k6RcY6ZvfosajCeKJJB53+epkpSJXcDSjafbORMz+GTy5ePpDoRhetyLgz26Udnujz3aevtaaaCKJ3z3N3okVJ9KH8hd6iegLjib9P9x4xP6Rfi72znDIKLSar++ORVuD/eydWEKjNTM3covk6QNlH8oioi84mnRP1Rn2Tm6efjabuBbtNvT0o+bmvhXhF1P0O80q5tEJifQFoWykR/pOKNDKpThrKBTLqaVCe6OPYDhuqywoq/LWivQ9bhdulyp4I3d4IsbStgZAIn1BKCvBcAyPy7AvnBDp51KcNZxjpN9m5urbqRXDdNEHI9ovJE8/HEsQjSdZ1m6KfpkzeET0BUcTDMdZYLYOdkSkn+NGbjbpmhZWK4YhG7VisCJ6n3ua6Bdg71jBxeIW43NX7s6kIvqCoxkNx1Ktg50Q6WdbnKW1Zmg8mores8FqxTAwXv6Co1JhRfpWqiYY/n4h9s6EmfnTbQYbYu8IQhmZiCVoqvPg97gcIfrZ2jtjkTjxpM55IxfsFumb9k56pO8tLNKfMD9nnQEfLiUbuYJQViKxJH6Pm0a/xxH2TrYpm8Nm8VCuG7kAg3aK9BOl8fTBaF3x4cvWsHFl2zyPKC7ZVV0Igk2JxBP4vS7qvW5HRPrZWvq5tGCwaK03zrUuGHbAEvepol8ce6fe5+ZjV60rbIF5IJG+4GisviqNfjchKc5KYVkOLfXZ2zs+j4tGn5thG41MnMzTnxrpRwtIS7WCiwafu7DF5YmIvuBoonHD3qn3eQgVWFpfC2Tr6VsZJbmIPhh2kB1TNq3iLDA9/QLsHcvTr/OK6AtC2UlF+j43oYj9PX2dpehbGSXN9bk5wK0N3lQzMTuQStmcYe8Ux9OvBCL6gqOJxBP4PS4afM7w9LN1JfKP9L32sncyZe94XEXz9CuBiL7gWLTWqUjf7y1sc65WyN7TNyqVc41GW+vtae/MyN4pQsqmRPqCUGbiSaODot+c+lRoE61aIN3emcvqGZmI0VLvzbrDpoXd7J3MG7nuwjz9mHj6glAR0gtvCvVpa4X0jdy5ov7RcJzmHK0dmLR3st07qHYypmx6C7N3wrEESk29kJQTEX3BsViDMPxel1lw4yx7Z65MnpGJWH6iX+8jkdQEbbIpPmtxVoH2ToPXnfO3qGIhoi84lkh88qu7U+yd9IZrc7VZHp2I0VyXe+2mNVrRLhZPxjYMBXwr3N83xveeODjlIlJuRPQFxxJJy8G2RN8utsRspAv9XJH+qOnp54rVtsEum7nWfNz0qNzvcZFIauJ5FGj9430vA0YH00ohoi84lmh6pG9uqhVSaVkLZO/p52nvmJG+XVoxROKJ1Hxci0KGoy9pqy/KugpBRF9wLOmFN8Uag1ftpLdWnq3NstY6lb2TK1ZXTrvk6luRfjqp4eh5fFa6m+qKsq5CENEXHMt0ewcmo3+7kq7zsw1UCceSxBKa5rrcRb8l1XTNXvZOOpMBQu4b/5Yl9P5Xryx4bfkioi84Fisdz8jeyT96qyWm2DuziH6+1bhgP3snmkjOSK1M2Tt55OpHExqfx8Xn3nxmUdaXDyL6gmOJJsyUTY8r7T+yvdM2p9o7mc/Jt+8OgNftIuD32Eb0I7Hi2juxRHJKJlAlENEXHEsq0k+zd+we6SdLHOlbj7ONvZOYKfqWaOdj78QSSbzuyuTnW4joC44lktZXxecQ0U9ksZE7aop+Pp4+GLNybbWRW8TsnVgiiUcifUGoDFakZhRnmV/Z7W7vZJGyWWik31rvs0+kb85bSGfys5KHpx/XYu8IQqWYXpGbfsyupOv8bPZOKtLP195p8NrH048nipq9I/aOIFSQVHGW1+2c7J3k/F02R8xRifm0YQAjV98u9k4kU8pmgfaOt9ojfaXUMqXUw0qpl5VSLymlPmIeb1dKPaiU2mv+3WYeV0qpbyil9imltiulzi/1ixCEfJjSZdObf/RWS0y1dzKfMxqO0ehz5+09W/bObHsGtUSmjdzJACGfSF9Xv+gDceDjWusNwEXAh5RSG4BPAQ9prdcCD5k/A7wBWGv+uQn4dtFXLQhFIGK2uPW6lWOKs6Z02ZxFlIdD+VXjWrQ2eElqbNFpMxrPkKfvyT9PP5ZI4q1gszXIQvS11j1a62fN20FgJ7AEuA64wzztDuAt5u3rgP/UBk8DrUqpRUVfuSAUiDU1SynlGHsnm4rckYloqnFaPliPtUOnzchcop93nn4NefpKqZXAecBmYIHWuse86ySwwLy9BDia9rBj5rHpz3WTUmqrUmprX19fjssWhMKJpGVmTEZvNrd3soz0rcrafGg1vyXYodNm5pTNQuydGvD0LZRSAeAe4KNa69H0+7SxI5STgae1vk1rvVFrvbGrqyuXhwpCUUiP4grZnKslsmnDMDxRoOjbqOnanL138mzDUBOir5TyYgj+nVrrn5uHT1m2jfl3r3n8OLAs7eFLzWOCUFWkp+NNVlnaW/SnDlGZy9Mv3N6xQ65+po1cj0vhUnnaO/EaiPSVMT3g+8BOrfXX0u66F7jRvH0j8Ku04+8zs3guAkbSbCBBqBrSI32P24XbpeyfvZNekZtB8422ytHiRPo17unHE0kSST2jOMvaAwrnYQXGEkl8nsp6+tkk4l4MvBfYoZR63jx2C/BF4CdKqQ8Ah4F3mvfdD1wL7ANCwJ8XdcWCUCQisanVlsacXHtH+ulpmpk8/VA0QSyhU758PrTYxNPPNB/Xot6X38jEWCKJx1XZSH9e0ddaPwHMdmm6IsP5GvhQgesShJITTSRTXj4UPvC6Fphi72Tw9C0fvpBI3+t20VznYXC8xkU/w3xcizqPi4m8Iv0a8fQFwY5EYokp6XjGwGub2zvz9N6xfPhCPH2AziY//WORgp6j0kTjs0f6dT53XqIfrQJ7R0RfcCxGiX2aveN1Oas4K0Okb+XWFxLpA3QG/PQHazvST+/NNJ16r5tw1OYpm4JgN6YX3jjB3knMk71TDHsHoCtQ+5F+ZK5I3+smnNe4RLF3BKFiROKZ7B17i356cB/PaO+Yol+ovRPw0Vfjoh+dJ9KfyCPSj0qkLwiVY3qvdCPSt7enn27vxBOZIn3DkimGvRMMx/NKa6wWLM/eqsBNp87rZiLHTC+tde21YRAEOzG9ba7PESmbk0Ify9BmcyQUw+9xUZdB6HKhs8kPwEANZ/BYLTnqM7wX9T53zi07EkmN1kikLwiVYmb2jv09/fQe+tEMoj8citFWQLM1i86AIfr9wdq1eCzPPtMFMJ+UzZj5zarqu2wKgl2JxKfn6ds/ZTM+T6Q/XGA1rkVnwLhw1PJmbtj81lfnzVyclavoWxdZifQFoQJorY3irGkpm3aP9OMJnbK0Yhlea6G99C1SkX4Ni761UVvnyWDv5LGRG0uJvnj6glB2YgnDX51h79jc048lkzT4DBGLZdjIHSmww6ZFV5Ml+rXr6Vv2Tr1vpuj7vUamVy7TwWIS6QtC5bBsnOkpm5l8bjsRT2gaTI96Nk+/0HRNMHzwgN9T05F+yt6ZJdKH3DptWt8MMm0MlxMRfcGRZKq2NCJ9u3v6yVTkOt3T11ozFIrSUoRIHwxfv6Yj/VTKZqY8feNYLr5+MGyMj2zKc+B8sRDRFxzJpOg7y9O3Gn55XGqG6IeiCSLxJB2NhUf6YLViqOVI35ihnKk4y8royaUOYcycGRzwi+gLQtlJVVtOy96JJzVxG1s8cbMi1Ot2zfD0B8yovMPchC2UzhpvxRCOJajzuDFGikzF+raUW6RvVDsHJNIXhPJjefrpbXOtrBY7+/rxpMbtUnjdakZzuf5xQ6A7AkWK9Jt8NS76yYzpmjAZ6eeSwWPZO811xbHP8kVEX3AkVpbO9H766ffZEaPhl8Lncc2wd6xIv7OxeJH+UCiWsR6gFgjHErNWJteLvSMItUVGT9+Te0ZGrRFPGpObDHtnuugXOdI3baJaHaYyMYfoN5rCPZ5HpC/2jiBUgMwpm64p99mRWELjcavMnr4pzu1F3MgF6KvRzVzD3sks+lYGzpgp5NkwFolT53VJnr4gVILoLNk7YP9I39jIVTP2LvrHIjT5PQU3W7PoaqrtVgyReGJWT9+yaMYi2Q9/D4bjBPyV9fNBRF9wKJkGZFgXADtPz4onNB6XGenHZ3r6xbJ2IL0VQ23aO1b2TiYsiyaYQ6QfDMdorrC1AyL6gkNxrr1jRPoZN3LHI0VL14Ta779jePqZJbLRZ0X6udk7lfbzQURfcCiOzd5JzuHpj0WLVpgFxmZnvdddswVa4VgyY98dALdL0ehz5+TpB8Pxilfjgoi+4FAsP3uqp++A7J2ENrN3Mnn60aJG+lDbufpz2TtgWDy5RPqjEzGaxNMXhMqQivTTJ2e5nWLvqBkpm8mkZnA8kuqDXyyMqtza9PRD0QQN/jlE3+8hmIPoD4VitDWK6AtCRUhV5KZ7+o7I3jHsHd800R+eiJHUFNXegdpuxTAeiae8+0wE6rxZ2ztaa4ZDUVqLMJWsUET0BUcSiSdxKfC4JvuqOMHTjyXSirPik57+ZGFWke2dGhX9eCJJJJ5MFWFlosmfvb0zFokTT2raRfQFoTJE4sbUrPRmWpMVufa1d6w2DN5p2Tv9qWZrxRWlroCPwfHolIHstYBVadswy0YuGPZOtpH+0LiRz1+MATWFIqIvOJLotPm44BR7J4k7w0auFY13FKnvjkVnk5+krr1WDKGoIeZzRfq5bOQOhYzXX4yh84Uioi84kkg8MaXDJqTn6dtZ9M2Ga9M8/VOjYQAWNBff3oHay9Ufj8wv+k11HkYnsqvIHbREv8h7Jvkwr+grpW5XSvUqpV5MO/Y5pdRxpdTz5p9r0+77tFJqn1Jqt1Lq6lItXBAKIRKbGelPZu/YU/QTSWMu8GTDtUnL5dRomDqvqyhD0dOxZuXWWv+d8Yhh7zTOYe+0NfgIRuJZdREdTkX6tWHv/AC4JsPxW7XW55p/7gdQSm0ArgfONB/zLaVUZQdCCkIGLE8/HaWUMTLRpp6+JU6p4qy0i1vPSJiFzXUZB4YUwqKWOgBODE8U9XlLzbhp7zTMkb1jRe2WdTMXg6anXxP2jtb6MWAwy+e7Drhbax3RWh8E9gEXFrA+QSgJhujP/Pgbc3LtGenHzc1UYyN3qqd/ajTMQlOgi8nC5jrcLsWxoRoTfTPSn6v3vZWJY23SzsXQeBSXguYif5PKh0I8/Q8rpbab9k+beWwJcDTtnGPmMUGoKiLxREbR93nctrV3rDGQHpcLn9tFNJFEa+NCYEX6xcbjdrGopY5jQ6GiP3cpsTZy5yrOsgqtstmkPjkapqvJj9tV3G9S+ZCv6H8bWA2cC/QAX831CZRSNymltiqltvb19eW5DEHIj0g8OaUwy8Le9s5kpN9U50FrODES5ncvnaR3NMLClvqS/N4lrfU1F+lbWTlzFWe152DvnBotzUU1H/ISfa31Ka11QmudBL7LpIVzHFiWdupS81im57hNa71Ra72xq6srn2UIQt5k8vTBSNu0baSftDz9yQ3bv/qvbdz8w21EE8mU/15slrY11Jzoh6yN3LkifdPeySbSPzUaZkEti75SalHaj28FrMyee4HrlVJ+pdQqYC2wpbAlCkLxicQy2zt+j9u+nr4Z6XtcipZ6Q7B2ngym7l/e0VCS37usvZ5TwXBO82QrTTYbuVah1XBWkX6kJHsm+TBvn0+l1I+AS4FOpdQx4O+AS5VS5wIaOATcDKC1fkkp9RPgZSAOfEhrXTv/0oJjMIqzMkT6trZ3jIuZNy3ST3eYV7SXRvRXdTaiNRweCLFuYVNJfkexGTdHG87lwfs9bgJ+TyozZzbCsQQjE7GqifTnFX2t9Q0ZDn9/jvO/AHyhkEUJQqmZiCVomEX07To5y8re8bhVKkqNp7VHWNpWGtFf3RUA4GD/WM2IvtH7fv5Mm46Aj755Cs9OjliFb9Uh+lKRKziSiVgi44AMv9e+2TuxtOwdK9JP74mTaWO7GKzsbATgQP94SZ6/FIxMxLIqVFvcUs/xeTKTjps1CotbRfQFoWJMRBMZB4Ab9o49RT/d05/e+GvTqvaS/d6A30N3k58DffYT/SVt9SlRn43DA8ZFYUVHY1HWVigi+oLjSCY1kXiSeod5+pPZO4p6rzvVduKNr1jEj29+VUl/99oFAfaeCs5/YpWQtei31tMbjMxpCR4eHMfrVrWdsikItcyEmUVS73Nm9o7X7UIplaoOLUeV6PqFzew+FayZFsu5RPpaQ8/I7NH+kYEQy9oaqqIwC0T0BQeSEv0MkX6d11VTqYW5kNrINcXHsniK3WQtE+sXNhGOJTk8UBsWT7aiv7TVKGibqw7hyGCoZOmw+SCiLziOCXNARiZPv97rJhS1p+hPHxG53sykyVSvUGzOWNQMwK6T1W/xJJKaYDie1TcgS8yPDGbezE0mNQf7x1nVWR1+PojoCw4knLJ3Zop+g8/NRCxBskZsiFwIpaZBGZna15y1EIBwGfYw1nQHcCnY1TNa8t9VKMGwkXefbfaOz+Pi4CyZSYcHQ4SiidRFrxqYN09fEOzGXPZOvSmIkXgy40WhlpmYNgLw2rMW8U9vjXPtKxaW/HfXed2c1hWYUgFcrYxMZC/6LpdiVUfjrJlJ1kXujIXVI/oS6QuOwxK/TKJvCaLVZdFOTEz7huNyKd61aTmtZerxvn5hE7tOVn+kn4vog1FxfLB/LON9O3tGcSkje6laENEXHIclfnUZIvn6lOjbz9ef62JXDs5Y1MzRwYmUfVKtDIdyG2K+qquRI4OhjBO0njs6zOkLmjLuH1UKEX3BcYTnsnfMYxM2zOAJVVj0rY3jPVWer291zezIcp7t+oVNxBKafb1To/1EUvPckWEuWNE2yyMrg4i+4DgsQW+YZSMXbBrpm51FXRXKF19vbmbu7Klu0beGuHc0Zjck/szFxut66cRU62r3ySBjkTgbV4roC0JFmYgaX8Mzb+Ta2NOPJjJe6MrF4pY6mus8Ve/rD4xH8bgUzfXZ5bms6gxQ73Wz/djwlOOP7TWGQ21a1VH0NRaCiL7gOOby9K10RjsWaIWiiTn7w5capRTrFzazq8oj/cGxKB0BX9ZD4t0uxcVrOvnh04d5cl9/6vgfdvayYVEzi1tLM5EsX0T0BccxYUbxc2fv2E/0J2JGj/hKsmZBgL29Y6nZvNXIwHgka2vH4hs3nIvP7eIPu3oBo7Pm1sODXLlhQSmWWBAi+oLjmIglcLsUXvfMj791IbCl6Fc40gdY2x1gZCJG/9j806YqRb8Z6edCg8/DmYubUxbPnU8fRgPv2Li0BCssDBF9wXGMRxIE/JnFz/L0J2wo+qFoomKZOxZruo189emZLtWEEennXrtw9tJWXjw+yt5TQW5/8iDXnrWoZINpCkFEX3AcY5H4rKJvZ3snPMvgmHKytttI29zXW52+vtaa/mCUjkBu9g7AlRsWMBFLcOWtj+Fzu/j0tetLsMLCkTYMguMYC88u+nUeK9K3X/ZOKJpgcWtlRX9Bs58mv6dqI/3RcJyJWIJFeQwxv3hNJ9981/lsPjjAey9aUZVRPojoCw5kLBKn0Z9Z/FwuY8CIXYuzKm3vKKVY3W1s5lYjVl/8hXmIPsAbz17EG89eVMwlFR2xdwTHMRaJE5hj6HWDz8242DslY21Vi74xxHxRS3WlWRYTEX3BcYxF4jTNYu8ABOo8jEfsae9UOtIHo/lYXzDCgb7qE/6TKdGvjtGGpUBEX3AcY+HZ7R2ARp+HsbC9RD8aTzIRS5RlNOJ8XH3mQlobvHzyZ9srvZQZ9IyEcSnoasp9I7dWENEXHIeRvTO7+AXqPIzZLNK32gVn2zmylKzoaOTN5yxm98lg1RVp9QxP0BnwZ6zhsAv2fWWCkIFkUjMejROom8Pe8XsYt1n2Tq494kvN8vYGgpF4al3VwpHBEMvbqzPrpliI6AuOIhRLoDUE5rB3An772TuWuFaDvQOk0hmPDs4+ULwSHB4IsbKK5tmWAhF9wVFYYj6XvdPo9zAWsVf2zsiE0fagtUpE34qmZxsoXgkmoglOjoZZ2SGRviDYBsurn8veaarzMBapLtuhUKrN3lnWbqREVpPoHx405tyu6JBIXxBsQ0r058neCceSxDOMv6tVRkLVJfpNdV4WNtdV1RStQ/2G6K8U0RcE+zAUMmyOlvrZG2pZ6Zx2KtAarrJIH4yJUy8eH6n0MlLsOhlEqcmmcHZlXtFXSt2ulOpVSr2YdqxdKfWgUmqv+XebeVwppb6hlNqnlNqulDq/lIsXhFwZNkW/bY7UxSbT+rFT2ubIRIyA34OnilIRz1zSwv6+saqZUrb7ZJCVHY1VUbVcSrL5BPwAuGbasU8BD2mt1wIPmT8DvAFYa/65Cfh2cZYpCMVhaNyIeNsa5or0DdGv9arc32zvYeWnfkMoaqRGVlOUD3D2khaSGl44Wh3R/q6TwdTwdjszr+hrrR8DBqcdvg64w7x9B/CWtOP/qQ2eBlqVUtXdfUhwFMOhKErNnbpodeAM1nja5q2/3wPA8aEJBvIYDFJqLlrdYU6bOlXppTAWiXNoYJz1C5srvZSSk+93vQVa6x7z9knAmgm2BDiadt4x89gMlFI3KaW2KqW29vX15bkMQciNoZAR8bpds88/tYu94zZnvCa0pjcYoSuPHvGlJOD3sOm0dh4yRwwCFavQfeHoMFrDuctbK/L7y0nBBp82/pVy/pfSWt+mtVJk7IsAABrOSURBVN6otd7Y1dVV6DIEISsGQ1Ha57B2YHKT1/L/axWXeWGLxJL0BcN0N1eX6ANcsb6bA33jHOwfpy8YYcNnf8cf9w+UfR3PHh4C4NxlIvqzccqybcy/rUv1cWBZ2nlLzWOCUBUMh6Lz9p+x7q+2FgG5Yn2ZGQ3HGBiP0tVUfZ0jrzjDMAl++2IP+/vGmIglePbIUNnXsfXwEGu6A1W371EK8hX9e4Ebzds3Ar9KO/4+M4vnImAkzQYShIozNB6bcxMXJtMah0O1LfqWhXVkMITW0F2FnSOXtTdw8ZoOvv77vTxlRvhHBmYv2LK+ERSTSDzBloODXLy6o6jPW61kk7L5I+CPwDql1DGl1AeALwJXKqX2Aq83fwa4HzgA7AO+C/xVSVYtCDly2qd/w7cf2W9G+nOLvtftIuD3pHL6axWX6elbRUfVKPoAX377OUTjSX74x0MAHBoYn/Xcy77yCJd/5ZGi/v5nDw8zEUvwmrXOsJnnHZeotb5hlruuyHCuBj5U6KIEoZgkk5qkhn95YBc+jyurLJbWBm+qirVWsSL9g/1G5FytPeIXt9ZPmaY1W2uGaNyokA4WeYP99ztP4XO7eJVE+oJgD6Jp7RSi8WRWWSytDd5UFWutYnn6VquDJa3VOwLw4jWdqds9I2HCGWYUHxsqfp8erTW/3dHDJad3plJ17Y6IvmB7otN66HQ2ZRHp1/tsY+8cGQwR8HuqNtIHuOHC5VN+Ppoh2j9cguZsT+0f4MRImDedvbjoz12tiOgLtseyBSy6AvNnsdjB3omlXexWdzWi1Oy1CZVm3cImPvWG9Xzw0tWA0dd+OtYGr2eOGovZODYU4uYfbp1RZf0fTx6io9HHG16xMI9V1yYi+oLtmS76WUX6Dd6aj/TDscnXfVpX9TcR+8vXream154GTG7mDoei/NP9O+kfi/CffzwEgN+Tu2x96YHd/O6lUzz48ikmzEZ6RwZCPLTrFO/atBy/x979dtIR0RdsT2y6vZOFp9/e4GN4IlbT7ZXTffHVXbXRLri1wUtTnSe1mXvrg3u47bEDfOjOZ9nfN865y1oZjyZIJHOrB02alb53bTnCGZ99gBeODvOdx/bjVop3b1pR9NdRzYjoC7YnPdJ3u9S8efpgZLpobVTw1irpon/+8rYKriR7lFKs6mxk2+EhkknNqNn/aPPBQZa21fPGVxitvHKdYWxdIrYcNNqI3fr7Pdy95QjvuWgFC1uqr2itlIjoC7Ynkib67Y2+OfvuWFjVq8UuBCon4bTXfV6NiD7Aey9awUsnRvneEwcYTcugunLDglRfpFyb4U0//5HdfSxoruN/v/70whdcYzgjR0lwNOnZO9k2HbMyXXqDEc4syapKj+VduxQ11SP+7Rcs5aGdvXzxt7tId3Gu2rCQwXHjm1cwHAOyT0HtGZ4cwP7OjUt59sgwX/+zc2mZpyWHHRHRF2xPLC3i7cwybdGqXq3VSP/qWx9jIpbgpktOq7loVinF1/7sHG7+YZzH9/YD0NHo45Ur21KtGsZyjPRPjoRTtz9/3VnUeWvnIlhsRPQF25Me6Xdm2VO+q8ZFf7dZkNVS762pKN+iwefhP//iQvb3jdNU5yEcS+Bxu1ID7Wezd44OhthzKphq5AZwoG+MYCROS72XVZ2NjhZ8ENEXHED6Rm62BUp1XjdNdZ6aFX2LWh4Eo5SaMa/WsueePzrMZeu7Zzzmll/s4Kn9Azz/2StpqjOsm/u2Gz0fH/joa1nUUr1VyeVCNnIF2xPLw9MHw+LpDYbnP7HKSJ85W4rWBZVkWXsDV5+5gP/30F4+d+9LU+7b2TPK43v7SSQ1mw8MEo0nuWvzEf7t4X1cvKZDBN9EIn3B9kTyiPTBaAR2fGhi/hOrjIGxyTTTj11ZW35+Nnzpbefgce3gB08dot7n5vL13bTWe/nSA7uo97rRaD7+0xfwul30j0XYtKqdf73h/Eovu2oQ0RdsT7q9k01hlsXStnr++8RoKZZUUgbMDJfvvW9jTVTi5kpLg5evvvMcookk335kP99+ZH/qvpsvOY3V3QHu3HyExS11vHPjMl53eldqipggoi84gFjCyPs7Y1EzZy1pyfpxS9saGBiPMhFN1NRm6OC4sQ9RbYPQi0md181337eRvmCEHceHGY8k6G7y88qV7bhcinduXDb/kzgUEX3B9kTjRr76nf9jU07j8KxWxMeHQ6zpbirJ2kqBZe90NFZvV81i0dXk5/L1C+Y/UUghG7mC7bFSNn05Nupa2maI/tEa8/X7TdFvt3GkL+SPiL5geyxP3+vOzddd3t4AwOH+2cf3VSOP7+1jWXs9jTVkSQnlQ0RfsD1R09P3uXP7uHc1+Wmq87Cvb6wUyyoJhwfGeWr/AG8/f1lV988XKoeIvmB7ovEkPrcrZxG0ioP291ZnpB+NJ1n5qd9MyV75+u/34ve4uP5C2cgUMiOiL9ieaDyZs7VjsborULWRvlV49S8P7ALgv186yS+eO84HXrOKBc3OahcsZI+IvmB7Yolkzpu4Fmu7A/QFIwyNV7av/qN7+vjn3+6ccuzItJmx//bwPtZ0B/hojTVYE8qLiL5ge6Lx/EX/FWZe/47jI8VcUs7cePsW/v3RA2g92Ws4fXj47146yfZjI7xn0/K8X6vgDOTTIdieaCKJN8dNXIszq0T0Lcajk9Ow0iP9m3+4jUUtdbztgqWVWJZQQ0hxlmB7ogXYO1Y73uePDhd5VfkxNB4l4Df+2+7tHWNNd4ArzujmsT39fOGtZ6U6SwrCbIjoC7ZnaDxKcwFiuHFFGw/uPEUyqSvew2U4FKOtMc5lX3mEvmCED166mr+5Zj2ffkNFlyXUEGLvCLZnz6kx1nbn33js1Ws6GA7FeLmn8s3XhkJRth8bTvX5v/mS0yq8IqHWENEXbM3AWIT+sQjrFubfO+fi1Z0APLyrl2T60NYSs683SDiWmHJsKBRlZ48xFWvL315Ba4O0WhByQ0RfsDV7Thk59qcvyF/0u5vruGBFG199cA/XffNJcyh3aRmLxHn91x7jEz99gUTahWY4FGNnzyidAR/dTZKLL+ROQaKvlDqklNqhlHpeKbXVPNaulHpQKbXX/LutOEsVhNx56YSRdbN+UWFdMt//6pWAkcXzH08eKnBV87PLtJLu294z5SKz+1SQR/f0cebi7FtEC0I6xYj0L9Nan6u13mj+/CngIa31WuAh82dBqAhbDw2xrL2+4Kj4T85ZzK5/uIZVnY3sLIO3n75/MByaFP27Nh8hGI7xf65eV/I1CPakFPbOdcAd5u07gLeU4HcIwrxordl6eJBXrmgvyvPVed2s6Q6wt7f0bRleTpvYteXgYOr22Utb+OWHLs5pGIwgpFNoyqYG/lsppYF/11rfBizQWveY958EZMKBUBH29o7RPxblwlXFEX0w2jI8vKuXWAEFX/MxOB7lNzt6uHBVOy8cHeaT92wH4LcfeS1nLGouye8UnEOhov8arfVxpVQ38KBSalf6nVprbV4QZqCUugm4CWD58uUFLkMQZvLQzl4ALl3XXbTnXLsgQDypeenEKOcuay3a86bz/ScOMBaJ849vOYv7XjjBz587zgcvXS2CLxSFgkIVrfVx8+9e4BfAhcAppdQiAPPv3lkee5vWeqPWemNXV1chyxCEjDzwYg9nLWlmYUvxslwuW9dNc52Hr/xud0nSN/vHIty5+QhXbVjA6Qua+NhV63jiby7n3ZtWFP13Cc4kb9FXSjUqpZqs28BVwIvAvcCN5mk3Ar8qdJGCkC2xRJJ7XzjBPduO8cKxEd52fnF70bQ2+Pg/16zniX39fOKnL2TdfTMYjrHt8OCc5zy+t4+3f/spwrEE/+uKtcVYriDMoBB7ZwHwC3MwhQe4S2v9gFLqGeAnSqkPAIeBdxa+TEHIju8+foAvPbA79fOfnlf8BmTv2bScvtEw33xkP0/tH+BXH7543v71H77rOR7d08cfP305wXA8VTcQjSd5cn8/d20+woMvn2JlRwN3/PmFkpIplIy8RV9rfQA4J8PxAeCKQhYlCPny/BGjMZrXrfjMmzbQ0lD8BmRKKT521TquOGMBN3z3aa77tyf53Js3cNWGhRl78ySTmkf39AFw9a2PMRqO8/ozFhCKxtlxfIRgOE5LvZdPXHU6//OS0/B7ZLatUDqk4ZpgK/b1jXH1mQu49c/OpcFX2o/3Octaufumi/jkz7bzl//1LMva67lsXTfrFzbTGfBR73MzHonz1P6B1GNGw3EAdvaM0t3s5w1nLeSasxZy8ZpOEXuhLIjoC3nzgR88wwUr2/irS9dUeimAYZUcHgjxhrMWllzwLc5e2sqv//o13L+jh3uePc7Pth0jFE3MOO+t5y3hTWcv4mD/OO/atLxs6xOE6cgnT8ibh3b18tCu3qoR/f19YySSmjUFdNTMB6/bxXXnLuG6c5cQTyQ5FYwwMBYhHEsS8HvobJI+OUL1IKIv2IYfP3MUr1vxqtM6K7YGj9vFktZ6lrTWV2wNgjAX0mVTyIt4IlnpJUxh66FB7tx8mLecu6SoefmCYDdE9MtAOJYgGq8ukSyU8chM37pSbD4wwI23b2FJaz3/940bKr0cQahqRPTLwPrPPMC133i80ssoKsFI6XvKZ8Ovnj/O+27fwqLWeu6+6VUlSdEUBDshnn6Z2FeGzozlpJiR/s+2HWNNdyCrXjZ3bzlCa4OXV63u5Iu/3cWPthzhwlXtfOc9F9DeKFOkBGE+RPQdQDKpUQqODU3Q1ugj4M/vn/3kSJiOgA+v28VYnpH+o3v6iCeSXHGG0Xz1xeMjfOKnL9DV5OeZv339nI8dmYjxqZ/vAKDJ72E8GuemS07jE1etw+eRL62CkA3yP6VATo6EGZyj/0o5Z6pmIpHUnHbL/Xz5d7t57Zce5obbns7recKxBBf980PcYoruWFqkH4lnH/XfePsWPnDH1tTP33/iIAB9wQiv+/LDTGTIcQc4MTzBZ375YurnK87o5t4Pv4Zbrj1DBF8QckAi/QK56J8fwutW7P3CtRnvH4/Gy7yiqQyFjAvStx7ZDxjj/vLh5EgYgJ9uO8aX33EOY+HJ1zUWjuMPzF9NGkvL+BmLxKnzuPjDrskmrIcHQtzz7DH+5JzFRGIJBkNR/rh/gN/vPMUf9w+ggfdetILPX3cmZs8nQRByRES/CMQSs0fzwXBlRX+ubyG5cGJkYsrP45F42u0EHdPqoZJJPaMPzYG+8dTtl0+M8vNnjzEyEePdm5Zz5+YjAHz+vpf5u3tfmjIMfG13gA9dtoZ3blzGsvaGorweQXAqIvoFoPX81k266Gutyx6h9o9F5j1Ha80tv3iRt52/hI0rM0+ZsiJ96/yxNNFPz+QZHI/y5z94BpeCn3/w1VNe7zOHJlsLv+d7m4kmkrz3ohX83Z9s4LJ13WxY3MytD+6hu9lPW4OPriY/5yxtZWVnY06vWRCE2RHRL4CJ2PxedvqGZyiaoDHPTdR8GRibP9I/MRLmR1uO8MvnjrPzH67JeE5Pmuj3BiNTRH8kFONbj+zjzecs5q7NR3jhqNHp8qn9A7TUe9lzKsgju/u4f4cx1OTN5yzmV8+f4NJ1XXziqnUopXj9BmNj98vvmNG4VRCEIiKin8aH7nyWPzlnMdectXDGfc8fHeb0BYEpjbKGQpOCPlsUPxpOt0HiFRD9mZF+LJHErRQD41G6mvzsPRUE5r6I9aTZO79+4QSnRicvAh/98fP0BiN86YHd+D0uzl/eyo7jI7z7e5tT57Q3+rj+wmV87Mp1tDf6uOmS1cV4eYIg5IiIvsloOMZvdvTwmx09HPriG6fcNxyK8pZvPsm1r1jIt959Qep4+tSk2aL49A3PYCRO8aa1ZkcmT39gLMrf3LOdR/f0sfmWK6bUEBzoG6OzyY9bKR7e3cvStgb6ghGe3DfAGYuaqfO6+Mff7ATgtWs7OT48kfLqlYLVXQH+6U9fweB4lJ09QRa31LGqq5HTu5sy9poXBKG8iOibHBkIzXrf0UEjyn1ib/+U4yMTk5H+UCiaUfSD07JcCmFgLMInf7adD166mgtWtKW+WWw9NIhSigtWtM14TH8G0b/y1kdT6/rru56bknJ5+VcfpdHnJhJPEk/bTG2q8/CN689j7YIAP37mKG0NPt61aTkupdjXO8aa7gBul8KdJuyvXl25xmeCIGRGRN/kyOBM0X/Xd59mUUs9l6834vP0jJJgOMZHf/x86ufhUIylaZp75+bDPLq7b4oQp/vgszEcinL9bU/zlvOW8Jevm2qBfOfR/al2xj63i09fu55QNMGXf7cbj0tx+/tfSXujj+FQjOPDIQ4NhPj18ydm/I5QNMEX//QVPH90mLufOYrP4+KdG5fidil2nQyytK2B5e31rF/YTFJrlrTWc9aSFuq8Rlrmx69aN+X5Nixunvd1CYJQHYjomxxOi/R/+PRhVrQ3pCYeWf3Z0yPf/3jyEH3BSb98ODS1QvVvf2EUEu1Ns06e3NfPplXteNyzFxPd8dRhdp0M8sXf7sLjUpy7rJVgJM6OYyN8/4mDvGZNJ0/s6yeaSPL3v34ZgLOWNLP7ZJD33b5lynN5XIoVHQ18403nsbozQL3Pjdet8Hvc1PvcXH/hcv7+ujPxuV2S9y4IDsHWop9Iav7vL3fw9guWZbQ+wrEERwZDnL6giQN9k+L8mV++OMWm+P3OUwBE4km+9/gBbrhwOXdvOTLluY4PT1409qc918H+cTataueZQ4N865H93Le9h9es7cTvcTE0HmUoFGNkIsZ4JE4wHOfkaJhLTu9iZCKW8s4tXru2k++85wL2nAoSDMfZcyrIyo5GLjm9i6FQlJd7RonEkrQ3+uhu8rO8vWFeH11G9AmCs1DZ5JqXmo0bN+qtW7fOf2KOPLW/n3d9dzNNfg87/v7q1PG+YASvW3HdN5/k8ECIW65dz1f/ew+nL2hKVay21HunePYuZQjkRCyB3+Mimkjy0StO59hQiG2HhzgyGGJ1V4BIPMGxoQl8HhfffPf5/GzbMT51zXpcLmX0fH/6CPv6xojGk7Q1emlr8NFS7yXg99Dg83DGoibec9EKtIathwfRGhr9Hpa01kufeEEQpqCU2qa13pjTY+wk+smkZjwaZ3A8yoH+cb7+4B5eOGaI+A0XLmNpWwM9IxP819NH8LjUFLumM+DjVx9+DSeGJzgyEGJpWz2P7Onj9AUBfvDkIT735jM5Z2kr33lsP3tPGcO3rzlrEQC9o2G+8+gBjg2F8HlcrOho4E/PX8rqrvKO7RMEwVk4TvQf3t3LP9z3MuOROGPhOOPTmnXVe9188NLV7D4Z5LG9famMldYGL2cubuY9m1awoqORXSdHuWL9AunFLghCTZGP6Ne0p99S7+WMhc0E/B4CdR4a/R6a/B5aGrwsb2/g7KUtU4qpwrEEXrdril8Pkn0iCIJzqGnRP395G+e/e+YG7WxYKYeCIAhORRqRC4IgOAgRfUEQBAchoi8IguAgRPQFQRAchIi+IAiCgxDRFwRBcBAi+oIgCA5CRF8QBMFBVEUbBqVUH3C40usoMZ1A/7xnOQt5T2Yi78lM5D3JTCfQqLXuyuVBVSH6TkAptTXXHhl2R96Tmch7MhN5TzKT7/si9o4gCIKDENEXBEFwECL65eO2Si+gCpH3ZCbynsxE3pPM5PW+iKcvCILgICTSFwRBcBAi+oIgCA5CRL+IKKWWKaUeVkq9rJR6SSn1kQznXKqUGlFKPW/++Wwl1loulFJ1SqktSqkXzPfk7zOc41dK/VgptU8ptVkptbL8Ky0fWb4n71dK9aV9Tv5HJdZabpRSbqXUc0qp+zLc56jPicU870nOn5OanpxVhcSBj2utn1VKNQHblFIPaq1fnnbe41rrN1VgfZUgAlyutR5TSnmBJ5RSv9VaP512zgeAIa31GqXU9cC/AH9WicWWiWzeE4Afa60/XIH1VZKPADuBTDNMnfY5sZjrPYEcPycS6RcRrXWP1vpZ83YQ4x9qSWVXVVm0wZj5o9f8Mz174DrgDvP2z4ArlFIKm5Lle+I4lFJLgTcC35vlFEd9TiCr9yRnRPRLhPnV8zxgc4a7X2V+tf+tUurMsi6sAphfT58HeoEHtdbT35MlwFEArXUcGAE6yrvK8pLFewLwNqXUdqXUz5RSy8q8xErwdeCTQHKW+x33OWH+9wRy/JyI6JcApVQAuAf4qNZ6dNrdzwIrtNbnAP8K/LLc6ys3WuuE1vpcYClwoVLqrEqvqdJk8Z78GliptT4beJDJCNeWKKXeBPRqrbdVei3VQpbvSc6fExH9ImN6tPcAd2qtfz79fq31qPXVXmt9P+BVSnWWeZkVQWs9DDwMXDPtruPAMgCllAdoAQbKu7rKMNt7orUe0FpHzB+/B1xQ7rWVmYuBNyulDgF3A5crpf5r2jlO+5zM+57k8zkR0S8ipr/4fWCn1vprs5yz0PIhlVIXYvwb2PaDq5TqUkq1mrfrgSuBXdNOuxe40bz9duAP2sZVg9m8J0qpRWk/vhljf8i2aK0/rbVeqrVeCVyP8Rl4z7TTHPU5yeY9yedzItk7xeVi4L3ADtOvBbgFWA6gtf4Oxof1g0qpODABXG/nDy6wCLhDKeXGuMD9RGt9n1Lq88BWrfW9GBfKHyql9gGDGB9wO5PNe/K/lFJvxsgIGwTeX7HVVhCHf04yUujnRNowCIIgOAixdwRBEByEiL4gCIKDENEXBEFwECL6giAIDkJEXxAEoUQopW5XSvUqpV7M4txLlFLPKqXiSqm3T7vvAaXUcKama7kioi8IglA6fsDMYsTZOIKRcnlXhvu+jJEOXjAi+oIgCCVCa/0YRv58CqXUajNy36aUelwptd4895DWejsZ+uxorR8CgsVYkxRnCYIglJfbgL/UWu9VSm0CvgVcXq5fLqIvCIJQJsxmjK8GfprWFdpfzjWI6AuCIJQPFzBsdlit2AIEQRCEMmC2Wj+olHoHGE0alVLnlHMN0ntHEAShRCilfgRcCnQCp4C/A/4AfBuj8Z4XuFtr/Xml1CuBXwBtQBg4qbU+03yex4H1QACjK+8HtNa/y2tNIvqCIAjOQewdQRAEByGiLwiC4CBE9AVBEByEiL4gCIKDENEXBEFwECL6giAIDkJEXxAEwUH8f+Xn9kKjZNaiAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
"plt.plot(freq, np.exp(-tau))"
],
"metadata": {
"id": "7G-SoEQbdHHT",
"outputId": "12b32b07-f218-4ae2-bfcb-5562ee4fbf86",
"colab": {
"base_uri": "https://localhost:8080/",
"height": 296
}
},
"execution_count": 13,
"outputs": [
{
"output_type": "execute_result",
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f75035d1710>]"
]
},
"metadata": {},
"execution_count": 13
},
{
"output_type": "display_data",
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEFCAYAAADt1CyEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO29d5hkZZn3/7krV3VOMz0zPTkPMDAwgggoWUReEMMKhtVd1HWVNcuFusvP9VV/ptewLrqLuqtrQBER8RUEBUQECUNmEkzOM51zV3zeP6pOdXV3dXelrnDq/lzXXHTVOVX11OHU99zn+9zPfYsxBkVRFMW+OEo9AEVRFGVuUaFXFEWxOSr0iqIoNkeFXlEUxeao0CuKotgcV6k+uLW11SxbtqxUH68oilKRPPXUU13GmLZsXlMyoV+2bBlbtmwp1ccriqJUJCKyP9vXqHWjKIpic1ToFUVRbI4KvaIois1RoVcURbE5KvSKoig2R4VeURTF5qjQK4qi2JyS5dHnypP7evjLy13U+VzUeF3Uel3U+lw0BzwsbQnQGPCUeoiKoihlRcUJ/dP7e/nW/S9Pu/2URQ18/NI1PPRSJ1v29dIzHMLvcbKgwcfqeXW88+yluJ3C4FiEde11iEgRR68oilJ8pFSNRzZv3mxyXRkbjRmGQxGGxiIMByMMBiN0DQZ5+cQQt/x5D/2jYQDOW91KW52X3SeGeO5QPwANfveE7T941yvY1z1MJGoIR2McGxjDGMNDL3XxydeupdHv5s5nD3N8IMjlp7SztKUGgL6REHc8fZiuoSA+t5MVbTVcvH4+PrezAEdHURQlPSLylDFmczavqbiIHsDpEOp9bup97gnPX3oSuJ3CF+/ewWmLG/nxdWcBEIrEePpALx6Xgzd+51EArjx1IXc9d4Qb73ieO54+DEBjwE3fSDj5fk/s7SZmYG/XMADff3gPn3/DyXQNh/jRo/vYdWIIl0OIxOIXy5MW1vOta07jd88f4x9eswIReP5QPz97/AD7u4cxQHPAQ2PAQ4PfTa3XSa3PxaYlTbxiWTOxmOFw3yhOh1Drc1Hndekdh6IoeVOREf1MdA0Fedd/PcEXrj6F0xY3Ttm+p3OIg72jvHp1K2+4+ZFkpD+ZlW017O6MC/wnX7uW81a38sbvPJoUdYBvXXMaV522iLFwlF89fYjP/PrF5La3bl7Mnc8eJhiJ0eB3c9LCepwOoXsoRP9omP7RMMOhCNbh/8Sla/jTzk627O9NvkdzjYe/P2cZH7xgFfu7R3C7HOw6McSqebUsavQzGopyfGCMqDE0Jy4eDodeGBTFzlRNRD8TrbVefveh86bdvqKtlhVttQCsba/juUP9NAXcXLqhndXza/n877YDsHlpM7s7h9mwoJ4PXrAKgNOXNPHEvp7ke528qAEAn9vJm07vSAq9CPxiy0EA/uHVK7j+wlXUTbr7ADDGMDAa4cO/eIav3fcSAB+8YCUdTQEGx8I8urubr933Er0jYX7wl70TXnvducu585nDdA+Hks8FPE6WNAdY217HJy5dS3uDD5dD9K5AUaoc2wl9Nlh++7w6H19+80ZiMcPnf7ed15+ygCUtAQDaG3zJ/T9/9clc+o0/Jx83pWT4+NxOrti4gIO9o5y8sJ6fPn4AgGvOXJJW5AFEhIaAm49evIZHd3Xz3lcv55OvXZfc/p5zV3DNLY9NEPkPXbSaFw/384O/7MXlEL7ypo14XA56R0Ic6BnhYM8Iv3n2CPdtPc5oOMriZj9DYxGWt9bwxTeewi+ePMgfth2nMeDmHWct5djAGINjEVprvXQ0+dm0pJGB0QiLm/3TjltRlMqiqoV+SXNczA1x/8ThEJ7/7KX43U5++9wRID4fYLFmfh1fvPoUPv3rFwCo9008fN++dhPGwFfu3Zl8rtY7+yE+dXEjz3/20ikTuQ6HcOPl6/jwz5/ha28+lbXtdTQGPIyGotz6xAGWtQa4cN38Ke/39u8/xiO7ugE42DMKQN/BPi775sMArJ1fx7H+IDfe8cKM4zpreTNfe8upBDxO/rSzk6UtATYva8YYQzhq8Lh0GYaiVAJVLfSLmvwApE5TWBO8VrS+uCkw4TWNgfEo1+WcKHQiggjUpVwA6nyZHeLpsnVOX9LEwzdcOOE5v8fJ35+7fNr3aqnxTnjsdAiPfeoiXvGFPwJw2/vP5kv37ODWJw6wqNHPv1yxng0LGugcGuPhl7tY2hJgf/cI3/vzHr5y705ePj7IjmODAHzk4tU8uOMEzx3qp6PJj9fl4MzlLXz68nUc6x/jSP8YHU1+FjcF9EKgKGVCVQv9yoRXf/2Fq6Zse82aNr549SlcvWnRhOdThX46UsXdWwKxa6mduGisKeCmrc7Lk5+5mK6hIA1+N6vnxb/7/Hovl528AIAlLQHOWNqcfN1T+3uTdzYfv2QND+w8wTf/GF/D4BA41Bu/W9jdOUwwEuU3zx4hmpis9jgdnLm8mas3LeLqTYt0klhRSkhVC32D382+L70+7TaHQ3jbWUumPN+UwcrbVKEvxURoa+3UiB6grc5LW118W0fibmbyXUkq7fXx+QmXQ3jvq1ewv2eEZw70AfDL95/Nj/+6n49fupZ/+c2LyRTVb7z1VGIx2H50gAd2nODjv3yOFw73c9LCev6yq4urTluY1m5SFGXuqGqhz4VMIvpab2knMZtrZr8Ynbm8mUWNfm547dpp91mQmIhe1OTH53ZOmG9Y2VbLN6/ZBMCHL1rNn3Z28soVzVy9qSO5z2dev55P3v48P3x0X/K5e144xreuOY3LTm5nf/cIdT4Xtz5xgOWttQS8ThY2+FnZVjPjBUhRlOxQoc+SbCP6UtCSgdA3Bjw8cuOFM+6zoDEe9dd44t8nVehrUv7etKSJ5266FJmkzSLCRevmcftThwC484PncMPtz/GBnz3N9Res4tsP7Er7uQ1+NxesbePyUxZw4bp5KvqKkicq9FliTZq21k4vpplk2swlfs/4GLuGQgi52UfzEjZPIPF+qeLuniS+DdPc6axbUJ/8e+OiBv7r3a/g3C8/OEHkv/KmjaxoqyEcNRwbGOWRXd08sOMEdz57hAUNPj7z+vVcsXFhTt9BURQV+pz42XvPYkVr7bTbJ5dmKDYLE5H4/zp1If/9yD7WtNfl9D6W13/60iYAanO4U7FSWCE+79HRFGDDgnq2HR1IPn/16YsmXDiu3tRBJBrjwZ2dfOX3O7j+Z8/wo0f38bqTF/CWzR2a368oWaJCnwOvWtk64/ZcBLGQrGyr5aFPns/ipgCvXt2WFOpsOXVxIz++7kxeuaIFgFpv9gXbnA7hE5euYVlrTfK5d569lE/d8QKnLW7kf647c8rdAcQniS/ZMJ/WWg9Xf+dRntzXy5P7evnPP+/mC284hYs36ISuomSKCv0cUGrrBsZX/V6wbl5e73Pe6rbk37lOMl9/4eoJj99yRgcvHO7nTad3zHr3c9riRt54+iJes6aNjqYAn/n1C7znf7bwTxeu4qMXr9G0TUXJgNIrkg2xFgqdkqiFYxdqcojo0+FyOvji1adktK+I8PW/OS35+K7rz+Vf7nyRbz+wi8O9o3zlzRt1slZRZkGFfo64/+OvSU5m2oVyuFPxuBx86U2nsLjZz9fue4lwzPBv15ymhdsUZQZK/8u1KdaqWztRDkIP8Sj/+gtXIyJ89d6dvHJFM28/a2mph6UoZYve8yoZUy5Cb/GB81dy+pJGvvun3USisVIPR1HKFhV6JWNKnU00GRHhXa9axqHeUV48MjD7C+YQYwwDY2F6h0OE9aKjlBnl9ctVyhp/YrHYP6UpAlcqrOYvu08Mpe0oNlcYY3hyXy+/fuYwW/b1sL9nhFBkXOBbaz3ccNk6VrbVEo7GkimqilIKVOiVjBGRaYvAlYolzQHcTmFX51DRPnN35xCf+tULPLGvh1qvi1csa+KCdfNoq/XicgqDYxHu33GCG25/Pvma5hoPGzsaaK/3sWFhPa89qZ359b4ZPkVRCoftesYq1cfFX3+I7kT55W+89TQ2LcltgVgm7Dg2wJu/+9fkQrA3ndFBwDM1XuoaCrL583+c8Nza+XV0DwfpGoq3f3zlimY+fNEazl6p0b6SOdozVqlK1i+o57fPHaF3JMyn7niBez58XsHTLY0xfP/hvfyfP+ykzufmzg+ew6JEqYl0pJaKfuiT57O7cyhZnnnXiSHueeEotz5xgGu/9xhXnrqQf7liQ7KEtKIUGo3olYrn+MAY//3IPvpHw9z6xAFu+4ezOXN58+wvzIL/+es+bvrNVi5eP49/verkGUXeYuexQfZ2DXPZye1pt4+Fo3z3T7v57kO7Wdzk544PnEODX+v4KDOTS0SvWTdKxTO/3seNr1vHTVdsoM7n4hdPHizo+3cOBvnSPTs4f20bt7xzc0YiD7C2vW5akYd4JdSPXrKGH/3dmezuHObHf91XmAEryiRU6BXb4Pc4uWjdPB7ceSLZ0rAQ3PnMYUZCUW66YsOc1NY5e2ULr1nTxg8f3c9IKMLND+7ivf+zhRcP9xf8s5TqRIVesRUXrZ9Pz3CIZw/2Few9dx4fZF6dlxVzuNr5QxetomsoyOd/t52v3ruTP2w7zrW3PMbgWHjOPlOpHjISehG5TER2isguEbkxzfYlIvKgiDwjIs+LyOWFH6qizM65q+IlpB/f212w93z5xBCr589tSYszljZz+pJGfvb4AQBufN06BoMRLvvmw3z2rq0MByNz+vmKvZlV6EXECdwMvA7YAFwrIhsm7fbPwG3GmE3ANcB3Cj1QRcmEphoPK9tq2LKvtyDvZ4xh1/FBVhWhdpFV7XRFWw0nLYx35jrcN8oPH93H1d95hAGN7pUcySSiPxPYZYzZY4wJAT8Hrpq0jwGsnnENwJHCDVFRsmPz0ma27Osh34yyXz9ziOWfupvhUJTV83Pr0pUN6xNtF5sDngnZN284bSG7O4e5eZoeu4oyG5kI/SIgNY3hUOK5VD4LvENEDgF3A/+U7o1E5H0iskVEtnR2duYwXEWZnY2LGxgYi3C4bzTn9zDG8H/uewmANfNrueq0ue9Zu3Je/K7hpIX1E4T+8lMWcGpHA88UcN5BqS4KNRl7LfBDY0wHcDnwYxGZ8t7GmFuMMZuNMZvb2tqmvImiFIK1iej7peODOb/HyyeGONQ7yscvWcNvPnhuUfrUbl7axLev3cSnLl8/QeibajysX1DP9qMDed+lKNVJJkJ/GFic8rgj8Vwq1wG3ARhj/gr4gJkbqyrKHGHZLDuP5V7/5s8vxe8433RGB35PYTprzYaI8L9OXYjP7ZxwYWkKuFm/oJ7BsQjv+/FTKvZK1mQi9E8Cq0VkuYh4iE+23jVpnwPARQAisp640Ks3o5SEBr+b9nofL+cR0T+6u5vlrTUszHBxVKFxpuTrNwY8yZW+f9h2nN+/eKwkY1Iql1mF3hgTAa4H7gW2E8+u2SoinxORKxO7fRx4r4g8B9wKvNto2KGUkHn1XrqHQzm9NhYzPLmvp2xKCzf63ayZX8e2z72WRY1+fv3MYb75x5e44fbn6MnxOyrVRUZFzYwxdxOfZE197qaUv7cB5xR2aIqSOx6nI+cGIPu6hxkci7CpiPXtZ8Jqfh7wuFi/oI77th3nvm3HATg+EOSHf/cK7ZmrzIiujFVsicflmNAIJBteSJQeOKWjoZBDKgiLmwMAeF0OPnbJGh56qZMn9vaUeFRKuaNCr9gSj8tBKMeI/oVD/XhdDlbPK22D93e/ahkXr5834bklCaFvb/DxnvOWU+d18dPEalpFmQ6tR6/YErcz94j+pRNDrJpXm7RMSsVnrzxpynNW5cymgIeAx8XbX7mU/3hoN1edtpCL1s8v9hCVCkEjesWW5BPR704IfTmyoq0GgCtPjS/g+uglq1k9r5Yv3L2diDYlV6ZBhV6xJd4cI/rhYHxFbTFq2+TCqnl1PHLjhfzdOcsA8LriNe33dA5z/44TpR2cUrao0Cu2JFfrZm/XMDBejqAcWdTon5Blc+mG+bTWern9qUMlHJVSzqjQK7bE48otvXJfd1zol7XUFHpIc4bL6eDCdW08vb8wFTsV+6FCr9iSXNMrD/XGC6F1NJdmRWyuzKvz0TsSIlbAzlqKfVChV2yJ2+lgOBTl6/ftzKo2zOHeURr8buqLUMSskLTUeogZ6Bu1f836aMywu3OI7UcH6BvRlcGZoOmVii3xuOIxzL89sIu3bF6cXGg0G4d6RzJu/l1ONNd4AOgZDib/thtj4Sj/+tut3PXsEYZDUQBcDuHmt5+OMYZ7tx5n1bxa3nJGB/PqfSUebXmhQq/YEq9r/Ga1ayiYsdAf7hutKH/eoqXGC8Bn79rG99+1GZ+7OBU3i8VIKMI/3/kidzx9mDef0cFZy5up8br4wu+288W7t3OgZwTrxu0bf3iJD5y/ko9cvGZOmrlXIir0ii3xOFOFPvPb+8O9o5yzqvIqbFtR/F92dXHH04d521lLSjyiwmGM4YM/fZoHd3by9rOW8IWrT0luGwlF+cQvnwPgjg+8iuaAh2/d/zL/9sAuYgY+8dq1pRp2WaEevWJL3M7xSK5rKJjRa4aDEYZDUdor8La/tXbcrqn32yt+e+5QPw/u7ORDF67i8284ecK2K09dyLw6L3VeFxsXNbCstYav/82pXLJhPrc+cUAXkSWw1xmhKAk8rnHronMwM6E/kdivrc47J2OaS5pSfPlohWfehKOx5HfwuZ0cTmRCvX7jwilVOj0uB19+80Z6hkLJkhUiwtWbFvGHbcd5/0+e5uyVLVx37vLifokyQ4VesSWeSR59JnRWsNC7U6yqYLgyo9ixcJSP/PxZfr/1GD63g1csa+bH153F4Fg8k6jWl16uLlg7b8pzVqOWP24/zh+3H696oVfrRrElqdZNphF9JQs9wH0ffTUAN/zqeT72i2dLPJrsMMbwsdviIg8wFo7xl11ddA4GGRyLAFA3jdCno6XGM2Ge5j8e2s0nfvkc/VWQfpoOFXrFlkzOusmEE4NjQHzxUSWS2vbwjmcmt3Uub2594iB3v3CMT1++jvn18QutMfDAjuMMBuNCX+vJXOhFZMJcxZfu2cHtTx3iQ7c+U5U9d1XoFVsy0brJLOumczCIyyE0+itrsZSFzzXx53xiYKxEI8mef3/gZc5c1sx7z1uBI8WHv2/rcQbHwtR6XTmkSk7cf2VbDQ+91Mlvnj1SgBFXFir0ii1J9ayzsW5aa70Vm3vtcjomWFZPVUjtm5FQhCP9Y7xmbRsigivxHV61soWHd3Wxr2s4K9tmMtZr//Odmzl5UT1fu29n1WXjqNArtiTVnx0KRhhNrKSciZ7hUMWvKpWUKPblE0MlHEnmHOgZAWBpS3xR279fezqv37iAT1++HofAgzs78xL6Oz94Dvd99NWsmlfLP124mkO9o/xxe3WVdFahV2yJZ5KNkYlP3zcapqmmMm0bi2iK/7yns/yF/rYtB3kgUUd/aXN8RfKpixu5+W2nc/KiBt66eTEAtd7shf6zV26gKeCmo8nPmvl1AFy8fj6LGv388NG9BfoGlYGmVyq2xD2pDWBnBmUQekdCrG+vn8thzTmpOfR7ErX1y5VnD/Zxw+3PJx8vaZn6/2dpohxFLt3Crti4kCs2LpzwnNMh/O3ZS/n/79nBzmODrG2vy/p9KxGN6JWqoCsDn75/JExjoLIjeov59V52nxgq6wyT/35kPKpuCrhpSDMJvqw1Lv69w4VLi3zzGR04HcJvn6ueSVkVesWWWJFta208Va9zFuvGGEPfqH2Efl17PcOhKH0j5Zk3HosZ/vJyV/LxkmkKyS1J2Dk9w4UrR9xS6+WVK5q5+4WjBXvPckeFXrElCxriufB/f+4yALpnSbEcDEaIxgyN/sqejLVY3hoXyNkucKVi5/FBuodDOBMZTkunsdUWJxrAtDcUdm3DpRva2dM1zP7u8ra3CoUKvWJL5tX7ePFfX8s/vmYltV4XvbM0qOhPRL52ieiXJfzuEwPlJ/SxmOGzd23F43TwN5s7gPHxTsbrcvK9v93MT95zVkHHcN7qeIXSh1PuKuyMCr1iW2q9LkSExoB7VgujLyn09ojolyYiemu1bznx4pF+Ht/bw6cvX8epHY3A9NYNwCUb5he8Gczy1hoWNfp5ZJcKvaLYgqaAZ9aI3tpul4h+eYsl9OUX0d+//QQicOVpi1g1rxaAdUXOfhERNi9r4ukDvWU9YV0oVOgV29NU46F3tog+UeyqUssfTGZevZeAx1mW1s39O45z+pImmms8bF7WzJ8/eQEnL2oo+jhOX9LE8YEgR/rL766n0KjQK7anKeCmd5asjaFkhUR7CL3P5WRenbfsJmOP9Y/x4uEBLlo/Xlo4Xf58MTh9SRMAzxyojFIR+aBCr9ieTKybYatCYh5L7cuBN25aBIDDIbTWejNaP1BMHnopvgr2onXzSzwSWNNei9MhbD86UOqhzDkq9IrtaQy4GRyLzFjIyiqFG6jwptpfe8up7Pz8ZQA0+N1lV3/92YP91PtcrJlfW+qh4HU5WdlWw85jg6UeypyTkdCLyGUislNEdonIjdPs8zcisk1EtorIzwo7TEXJnaZEJk3fDKI3HIxQ43FWbOVKC4dD8CbaKDYEyk/otx7p5+RFDVNaApaKde31bD+qQo+IOIGbgdcBG4BrRWTDpH1WA58CzjHGnAR8ZA7Gqig5YWXSzJRiOTQWqXjbZjLlFtGHozF2HBssycTrdKxtr+Nw32iyXaFdySSiPxPYZYzZY4wJAT8Hrpq0z3uBm40xvQDGmOqqAaqUNfWJCdahhD2TjqFQhJocKiSWMw1+N0PBCOEyqb2+v3uEUCRW9FTKmViRWG+wv3ukxCOZWzIR+kXAwZTHhxLPpbIGWCMij4jIYyJyWbo3EpH3icgWEdnS2dmZ24gVJUssAbcya9IxNBahzmZCb6WKDpRJVH9wUt35cmCZCn1WuIDVwPnAtcD3RKRx8k7GmFuMMZuNMZvb2toK9NGKMjNWLfOZIvrhoA0j+oRlVS72zcHeuJgubiofobcuOvtsXvMmE6E/DCxOedyReC6VQ8BdxpiwMWYv8BJx4VeUkmN1J5rRurGj0PvLS+gPdI/gdTloq/OWeihJAh4X8+q87Cvz2v35konQPwmsFpHlIuIBrgHumrTPncSjeUSklbiVs6eA41SUnBm3bmaYjA3az7pp8M+ebVRMDvaOsLg5UDYZNxZLWwLs76ly68YYEwGuB+4FtgO3GWO2isjnROTKxG73At0isg14EPikMaZ7rgatKNlQ442nGw7P0DfWzhF9OXj0e7uGeXxvT7J8cjmxoMHPMZuXQcjozDbG3A3cPem5m1L+NsDHEv8Upazwupx4XA4Gp5mMNcYwHLRneiWUh3Xz08f2MxKKcuPr1pV6KFNY0ODj9y+OYYwpu7uNQqErY5WqoNbrYiiYXvCCkRjhqMmpAXU5E/DE72RGZ7iTKRYvnxhizfxaVraVfkXsZBY0+AhFYwXtYlVuqNArVUGt18VwML3gWUJoCaNd8CfKOYyUgdDvOjHEqjIUeYD2hnit+6M2tm9U6JWqoNbrmta6GQnbU+jj5RAcjIZLK/TDwQiH+0ZZPb98FkqlYrWdzNWnf3RXFz97/ECyMF45Yq97VUWZhpmsm9FQ/Afqq/CCZukIeJwlt252dw4BJJuMlBuW0B/tH836tQ/uOMHf/fBJAH751EF+8b6z8bjKL34uvxEpyhxQ65veuhlJWjf2i3v8bmfJrZuXj5e30DfXeBCBrlkayKfj9qcPMa/OyxevPoVnDvTxw0f3zsEI80eFXqkKaryuaRdMjdjUowfwe5yMldi6efnEEG6nsLS5fFbEpuJyOmj0u+kezr52//OH+ti8rIm3nbWE81a38r2H95ZNbaFUVOiVqiDgnt7CsJ7321DoAx4XI6HSese7TgyxvLUGl7N85aal1kt3lhF9z3CIgz2jbEw0OH/3q5bRORjk/u3lV9OxfI+8ohQQv8c57aTkqE0nY6E8rJtdJwbL1raxaKnxZC30H/nFswC8ckULAOevnUdTwM29W48VfHz5okKvVAU+9/RCn7Ru3Db06Ets3YyEIuzvGWHt/PqSjSETWmu9WVk3/SNh/vxSJ/94/kpOWxyP6J0O4YK183hw54kZu5mVAhV6pSrwu52EIjGiMTNlWzLrxmO/n0PAU9qI/qXjQxgD6xaUZ2qlRUuth+4sFkxtPdoPjEfzFq9Z20bfSLjsulbZ78xWlDT4EyKeLrq1e9ZNKfPodyQab69vL++IvrnGQ99IOOOJ1G1H4t/rpIUTv9crljUDsGV/T2EHmCcq9EpVYK0STSd6ltD7bZhH7y9xHv2OY4PUeJx0NPlLNoZMaKmNl07uzTCq33Z0gHl1XlprJ5ZcXtjoZ2GDj6f29xZ8jPmgQq9UBdZiqHSiNxaO4nU5cFZ4Y/B0lHoydvvRAda215V903WrG1emJZ33dg1PO8G8aUkTzx3qK9jYCoEKvVIVWKmT01k3dsy4gcTK2HCUeIHZeKXOv7zcRSzNXEWhMcaw49gg6xaUt20D2Vf63Ns1PG3J5fUL6jjYU14Nx1XolaogWclxWqG3nz8P4E98r7Fw3Hu+5c97eMcPHueBHXOf631sYIz+0TDry6gZ+HQkhX5kdnHuHQ7RNxKeQejjF7aXjpfPhKwKvVIVzGTdjIYjtlwsBeB3x3/i1qKpnz5+AEh/wSs0OxKZJ3aL6Pck2g6uaEsv9Nb33VZGmTcq9EpVMNtkrB0nYmH8AheMxCP6A4mWedNV8iwk24/FM1PWVlJEn4HQH+iJC/2S5vRCv7DBR43HyZ5EMbdyQIVeqQpm8ujHwvYVeqvsQCRqCEXGUwcHiuAf7zg6yKJGP/U+95x/Vr7UZyH0h3riVS6nyyQSERY3BzjQXT59aFXolapgpiYcwUgMr9uePwW3M57tEo7FJohYMdoL7jg2wPoyXyhl4XQIdV5XRsflcN8orbWeGctal1vDcXue3YoyiZmsm2A4hrcMa4gXApdjPKJPFbHv/mk3P3ls/5x9bjgaY0/nMGvKtNlIOur97oyFflHjzOsClrbUcDJ0A5cAABtOSURBVKBnpCjZTZlgz7NbUSbhm6F/ajASxeuyq3WTiOijMfpHJy4G+uc7X5yzzz3YM0IkZsqyR+x0NGQq9L2jLJplAdiS5gChSIzjg+XRnlCFXqkKrIg+nUcfjNg5oo8LfSRmimLXWOxNZKYsnyYzpRzJROiNMRlF9JZ/f6Qv+65Vc4E9z25FmYTb6cDlkKrz6K3J2GgsRl+aHPG5apKxpzORgjhNrnk5Uu93zbrIaTAYIRiJMb/eN+N+7ck+tNk3M5kL7Hl2K0oafG5nMs0wlWDYvtaN22FZN+kj+rmKOPd0DdMUcNMY8MzJ+88FNd7p201a9A3Hj+Fs32t+XULoB9S6UZSi4nU5CEam/pBDURtbNynplVZE70qpO7N/jlIAD/aMsKRMWwdOR00G3bh6RuLzHM01M6eMNgbceFwOjqvQK0px8bocE3LJIe652tqjn5ReWedzcdf15/LBC1YCcxfRH+4bpaPChD7gdc4a0VvVLZtmiehFhPZ6H8f6VegVpah401g34ajBmPg2O+JOSa8cGAtT73OzYWE9H75oDSJzYy3EYobDvaNlX5p4MrUeF6FobEowkErPsBXRz25Jza/3akSvKMXG43QQDE/8EVtWjt0j+khCwKxJZ4/LQUuNd04izs6hIKFojI6mSovo4wXgZqrf35uwbpoyEnqfCr2iFBuve6pHb0X4HpsK/fjKWEM4GsPjHP+eCxp8cxLRH+qN+/4VF9F743d1QzP49D3DIVyJVbSz0VLjSd4BlBp7nt2Kkob4ZOzkiD6W3GZHxlfGxghHDe4UoZ8/Rx7y4b74e3bMkmteblilqkeC0wt970iIphoPIrM3Ummq8TAwFimLRuH2PLsVJQ1el3OK/xoMR5Pb7Mi4dROP6K0IH+Yuou8cjOeOt9V5Z9mzvKixIvoZhL5nOERTILMibdaEbaZdq+YSFXqlavBUYURvRfDhWNyjT43o2xt89I2E064WzoeuoSBupyRL/1YKNVZEP4NHPzgWybgap+XjZ9qHdi6x59mtKGlIl0efFHq7rox1TIzoU+ciWmvjQtRdYCHqGgzSUuPNyN4oJ2oSvvvwDBH9UDBCrS+zbmRW5N+bQdequSajs1tELhORnSKyS0RunGG/N4mIEZHNhRuiohSGdB59KBnR29W6SUT0aTz6lpq4tdI9VNhl+l1DQVrrKmdFrIXVbnJ4hsnYobEItRlMxMK4dWNl6pSSWYVeRJzAzcDrgA3AtSKyIc1+dcCHgccLPUhFKQRpPXq7p1emFDWb7NG3WBH9UIEj+qEQrbWV5c8DSQGfadHUYDBCXaYRfYVZN2cCu4wxe4wxIeDnwFVp9vvfwJeB8kgcVZRJxNMrJ0/G2j2ijwt7NGYIRSd69JYYd81FRF+BQm/l0c9UBiG7iL6yrJtFwMGUx4cSzyURkdOBxcaY3830RiLyPhHZIiJbOjs7sx6souRDfMFUdXn01srYuHUzMY++ZQ48emMM3RUa0QfcVtZN+og+Eo0xGo5S681sMtbvduJ1OSrDupkNEXEAXwc+Ptu+xphbjDGbjTGb29ra8v1oRcmKtBF9wrpJFUA74XAIDklMxkYmevQBjwu/21lQj35gNEIoGktO9FYSDofE53GmyUKyLJ1MJ2NFhHq/e9bSx8Ugk7P7MLA45XFH4jmLOuBk4E8isg94JXCXTsgq5YbX5SQSM0RT2rvZPaKH+IRsOBaP6N2uiZkwLbWegnr0nUOVmUNv4XM7p003HQzGBTuTVbEWdT4Xg2MzV8QsBpmc3U8Cq0VkuYh4gGuAu6yNxph+Y0yrMWaZMWYZ8BhwpTFmy5yMWFFyxJpwTZ2QtfuCKYjXpI9Ep3r0AC21XroKaN1Yfn8lWjcAPreDsXD6lazWQqpMI3qIXxQqQuiNMRHgeuBeYDtwmzFmq4h8TkSunOsBKkqhsHLIU3Pp7b5gCuIRfSSNRw/QWuMpqHVT+ULvTNtAHuITsUDGk7EQvyjMtNK2WGQ0YmPM3cDdk567aZp9z89/WIpSeKyoPdWnrwahdzslUdTMpInoPWw9MlCwz+oatIS+8jx6iE+gTm/d5BLRu+kcHCrI2PLBvme3okzCEvPUUsWhSAynQ5ILi+yIyxFvuBKNpRN6L93DQYwx07w6O7qGQjgdMmtjjnLF63YyNk09eiuiz8ajr/W5kq8rJfY9uxVlEtaEayiaat1EbR3NQzyX3qqxPmUytsZDOGoYKJAYdQ0Faa7x4HBUVvkDC5/LMW1Eb1kwNdkIvdeVvBMoJfY+wxUlBcu6SZ1ss3MbQQu305FcBDTFo0946YWqm16pi6UsZsq6sS6WVqmETKhPePSxWGHumHLF3me4oqQwPhmbmnUTs3XGDcTLIFgVGdN59FC4ejedQ6GK9edhZo9+LDGJ78ui7WStz4UxMFLgCqHZokKvVA3etFk3Udt2l7JwOR3JTJIpQl9jlUEoUEQ/GKStoiP66dMrx0LZ10WyVtGWetGUvc9wRUnBmy6irwLrZmJEP3XBFED3cP4RvTEmUbmykoV+BusmHMXvdmZVftkqgFbqCVl7n+GKkoJl0YQmC72NV8XCxMnYyXcvVnOQvgIU3hoKRghGKrP8gcVMefRj4Ri+LM8VKxWz1BOy9j7DFSUFS9AnRvRR23v0bsf01o3P7cTvdtJfgHZ3lv1T6ZOxwWmsGyuiz4Zk16oZSh8XAxV6pWqwMk5Si1bFJ2Pt/TNwOSWZdTNZ6AEaA276ClBhsdJXxULcow9FYxPqIVmMhqP4ssi4gXGrLFziBuH2PsMVJYV0EX0oWg1CPz7BONmjh7h9U4ia6eOrYitZ6K3V01Mj8GAOEb07pcNXKbH3Ga4oKaQtgVAl6ZUW6coxNwbc9BdC6K2IvgLbCFr4Ehf90TQNwkfD0axSK2F8TiQc1Tx6RSkKaatXRqK2n4x1pGSJuNPcvTT6PfSN5m/ddA6FEIHmCi1/AOBPWDPpyiCMhjSiV5SyJ+nRT6peaXfrJjWIn96jL0xE3xzwVHTdICtiT5diORqOZR3RW1ZZSIVeUYqDwyHxdoKT0ivtvmAqNaJ3palB0xBw0zcazruwWddgZZc/gHF7L511EwxHs06v9GhEryjFJ94qbmLjEbt79KkFxtJG9H4PoUhs2hWhmRJfLFW5tg2kFr5LY93kMxk7TUXMYqFCr1QVXrdjUvXKKrBuUiL6dK5KYyCxaCpPn76rQpuCp+J1Tp3HsRgNR5Mefqa4dTJWUYqP1zW+ICYSjRGJGftH9CluTbrl+40FWh1b6ZUrYTxLJp3Qj+UQ0VtWmXr0ilJEPK5xj9768dk+6yZF6Z1phL4hkL/Qj4QijISiFS/06VJwAWIxw1g4hlezbhSl/PG6HMmsm1AVtBGEydZNuog+7qv352HddCfKH7TUVLZHP11Ebwl/thG90yE4HaJCryjFxJsS0Y/3i7W3dZMq7uk6PzUWIKK3Xmu9V6WSFProxKwbq1ZQtlk3kOjZqx69ohQPr8uZjNYsr97uEX2qL5/OuhmfjM1D6BN3A022jeizbzpi4XY6NKJXlGLimRDRJxpJ2NyjT820caT5qn63E4/TkVdEb9XKsSZ2K5V0PQsgv6DAo0KvKMUl1aO3fszp6r/YCecsEb2I0BBw5+XR9yeqXzZWcPkDmN2jz8XmczsdhCNq3ShK0fC6xxdMjUf09vboZZbJWICmPMsgWK9tqPCIfrxMxkShz2fi3u3SyVhFKSpel3PcuqkSjz5V3Kdrg9fo9+Rt3dR4nBVfTsIzzYKpfGw+t9OhefSKUky8Lsf4ZGy1pFc6Zo/o6/1uevNoPtI3Gqp42wbS10OC/Kwb9egVpch4Jnj0iSjN5umVqUF8Oo8eoN7vYjCPBtb9I+GKT6208KQEAxbj50puEb2mVypKEUmbR2/3rBtJzaNPv0+D381AHumVvSMhewn9pDz6pM2Xw7nicqpHryhFZYJHX43WzXQRvc/NYDCStldqJvSNhm1h3UDcail01k262jnFxN5nuKJMwutyEI0ZItFY1ayMzSTrxsqWGRzLLarvGwlXfA69hdedzqPP3bpRj15RikxqvfFguEoWTMnsWTf1CZHuz8G+icUMfXaybmaM6LUEgqKUPck86XCsehZMZfD1rIh+YDT7CdmhUISYgSa7WDfpJmOTHr2NSyCIyGUislNEdonIjWm2f0xEtonI8yJyv4gsLfxQFSV/rB9qMBKrGo8+XSGzydT7XEBuEX3fsD0WS1nEJ2MLmHXjqgChFxEncDPwOmADcK2IbJi02zPAZmPMRuB24CuFHqiiFILxWiZRgpEoHpdjWjvDLjgy+H5WTfqBHDx6q6CZXSZjUzOzLIKRGA5J33N3NjxOB7s7h/n9i0cLNcSsyeTydCawyxizxxgTAn4OXJW6gzHmQWPMSOLhY0BHYYepKIXBmngNRWIEw/ZvIwjTZ9qkUu/L3aO3VtQ22cWjT8nMsoi3nHTmFBRYTdff/5OnCzK+XMjkLF8EHEx5fCjx3HRcB9yTboOIvE9EtojIls7OzsxHqSgFwpNSndD68dqdTKybcY8+e6HvTRY0s4nQp5uMDUdznrR/+cRQIYaVFwUNZ0TkHcBm4KvpthtjbjHGbDbGbG5rayvkRytKRqRaN6EqaAwOE3vGTkfA48TpkJwieus1DX6bWDduB6HIpAVTeZwrW48MALCyrSbvseVKJiM/DCxOedyReG4CInIx8BngSmNMsDDDU5TCkhT6cIxgJPcorZKYLnc+FRGJr47NxaO3SXcpC2+aImT53P39xzvOAKDOV7rjk8lZ/iSwWkSWi4gHuAa4K3UHEdkE/CdxkT9R+GEqSmGYnHVTFdZNhr5yvc9Ffw7plb0jIWq9rmQj7ErH4xovZW1hTdznwmUnt3PphvmMhaOz7zxHzDpyY0wEuB64F9gO3GaM2SoinxORKxO7fRWoBX4pIs+KyF3TvJ2ilBTvFI/eHuI0E5kKfa71buxU0AymSa/Mc+Le73Em+86WAlcmOxlj7gbunvTcTSl/X1zgcSnKnOBJTa8M5x6lVRKZBtr1fnduWTej9hJ6b9rqlXkKvdvJaKiMI3pFsRMa0U9PfY4efe9IiEabTMTC9GWK87H5fO7SRvT2P8sVJQXrxxqMxBgL5/fjrRQymYyFeC69WjfgcTqJxMyESp7BSCyviXu/x1neHr2i2IlkUbNIjNFwlIDH/kJvRfSzBfZxjz6SXOCTKXazbtI1CM/Xow+4nYSjpmSlEFTolapivPlzlJFQlBpvFQh9IqKfTb/r/S5C0Rhj4czFKFm50kbWjTed0Odp3fgTAUWp7BsVeqWqSM2jHw1F8bszykeoaDIpgQApq2Oz8OkHg/HKlXaM6IMpXabync/xJdJ6x0o0IatCr1QVIoLH5WAsEmU4FKmOiD6h87PpfS71bvqTi6XsE9F7UoIBi7w9erdG9IpSVLwuR8KLHr+ltjOZ1LqB8XryvcOhjN87WefGJiWKIcW6iaZ69IWxbgZGI9y25SB7Ootb/0aFXqk6vC4nfQmBqvGodWNh2S+9I5lH9D2J49hca5+IPr1Hn38ePcDhvlFuuP15HtvTk98gs0SFXqk6fG4H3UNxgaqOiD7+39nkvqkmLtbWRTATehLHsdmG1o0l9JFojEjM5J1HD3BicAyIT3wXExV6peqo87k5nvjBVVN65Ww05RDRW9aNdZGwAx7n+FoLGLdw8vHorfPsWH9C6Itc4EyFXqk66nyu5A+uKqybDD16v9uJ1+VIincm9AyHcDkk2YrQDkyO6JP9YvPJo08I/VFL6Is8p6FCr1Qd9T53MlqrCusmw4heRGgKeLKejG2q8diqHeP4ZGw8Q2a8t3Du54ol7Id64434in1hVKFXqo7UH1k1RPSZCj3EJ2SzmowdDtnKn4epEX0okn9Eb61RONgzCmhEryhzTl2K0FdDRJ+pdQPxFMusJmOHQzTV2Ce1Eia2m4z/Nx7Z5+PR+9xOPC4Hxwbi1k2dRvSKMrekdvqphgVT2fQDaapxZ+3RN9toIhZSy2TEJvw33wJ4VlTvczuKXkxPhV6pOlJT2wJVUAJBkkXNZo/sGwOeLLNuwrYT+tTCd5AS0edZ0tpaVNZQgsVlKvRK1ZEa0VeFdZOFR9+csG5isdkrWEYTBc3s5tF7E+mVhcy6gXGBL3ZqJajQK1VIqj9aHR2mspuMjRkYHJu9d+zAaJiYsVcOPaRMxkYnWTfuwlg3xZ6IBRV6pQqxIvqOJn+JR1Icssl8TNa7ycCn706kYdrNuplc1KxQ1o0l9K0lKBehQq9UHZbunb2ipaTjKBZZZd3UWKtjZxd6ax+7Cb3TIbgckiaPPj+57EpcGM9d3ZbfAHPA/jNRijKJc1a18pnL1/O2s5aUeihFIRuP3io33JfBhGxPQriabObRw8S+sUmPPk/rJpB4/QVrVegVZc5xOoT3vnpFqYdRNJJZNxnsa4l2TwarY3ttat1APHqfkkefZ0T/hatP5pozF9PRFMh7fNmi1o2i2JzsFkxlbt1YJYrtGNEHPC6Gg4W1blpqvZy/dl7eY8sFFXpFsTnZWDf1PjcOydC6GQrhdzttmaJa53MxFIwfA0vwfXlaN6VEhV5RbI4ji1+5wyE013jpHAzOuu+xgTHm13vzGFn5UuN1MRSMp5gOjoXxu524s1liXGZU7sgVRcmIbIqaASxo8CVrsszEsf4x2ht8uQ6rrKn1uhhKRPKDY5Gi16YpNCr0imJzsvHoAdobfBztH511v6P9YyxssOdahFqvi6GxuHUzGAyr0CuKUt5kG9EvbPAlG2RMRzRmOD5g74h+eEJEX9kVOlXoFcXmWAF9pnrf3uBncCyS9KjT0TUUJBIzLGi0aUTvG/foB9S6URSl3MnWulmQiNKPzRDVWxH/gnp7RvTWZGwsZhgcC5ekEFkhUaFXFJvjyFLoFyai9IM9I9Pus/vEEADLWou/+KcY1HnjEfxIOKqTsYqilD+O5MrYzAR/bXsdANuODky7z/ajA3hdDpa11OQ/wDKkJiH0Q2MRBsd0MlZRlDLHWjBlmL3GPMSrLHY0+WcW+mMDrG2vw1XBueUzUZsQ9t6REGPhWHVMxorIZSKyU0R2iciNabZ7ReQXie2Pi8iyQg9UUZTcyDLpBoCTFtaz7Uh6oY9EY2w7MsC6RORvR2oTLSZL1eO10Mwq9CLiBG4GXgdsAK4VkQ2TdrsO6DXGrAK+AXy50ANVFCU3rMnYbGyWs1e0sLdrmKf290zZdv+OE/SOhLlw3fyCjbHcsGrHP7WvF6DiI/pMLlNnAruMMXsAROTnwFXAtpR9rgI+m/j7duDfRUSMMZndKyqKMme4nQ6+97ebOXVxQ8avecvmxXzz/pe59nuPs6DBh1MEEYgZONI3SkeTn4vXl6ZAVzHY2NHI8tYa/v3BXfjcDs5a3lzqIeVFJkK/CDiY8vgQcNZ0+xhjIiLSD7QAXak7icj7gPcBLFlSHbXAFaUcuGRDdtF3jdfFre99JT9+bD/DwQgxAzFjcIjw6tWtvOe8Fbb15yF+cfzWNafxvYf3cvWmhSxuruzsoqIaT8aYW4BbADZv3qzRvqKUMesX1PPFq08p9TBKxsaORr597aZSD6MgZHJJPgwsTnnckXgu7T4i4gIagO5CDFBRFEXJj0yE/klgtYgsFxEPcA1w16R97gLelfj7zcAD6s8riqKUB7NaNwnP/XrgXsAJ/JcxZquIfA7YYoy5C/gB8GMR2QX0EL8YKIqiKGVARh69MeZu4O5Jz92U8vcY8JbCDk1RFEUpBPadNlcURVEAFXpFURTbo0KvKIpic1ToFUVRbI6UKgtSRDqB/SX58OLRyqTVwYoekzToMZmKHpOpWMdkqTGmLZsXlkzoqwER2WKM2VzqcZQTekymosdkKnpMppLPMVHrRlEUxeao0CuKotgcFfq55ZZSD6AM0WMyFT0mU9FjMpWcj4l69IqiKDZHI3pFURSbo0KvKIpic1To80REFovIgyKyTUS2isiH0+xzvoj0i8iziX83pXsvuyAiPhF5QkSeSxyTf02zT1U1lM/wmLxbRDpTzpP3lGKsxUZEnCLyjIj83zTbquo8sZjlmGR9nlR2a/PyIAJ83BjztIjUAU+JyB+MMdsm7fewMeaKEoyvFASBC40xQyLiBv4iIvcYYx5L2SfZUF5EriHeUP6tpRhskcjkmAD8whhzfQnGV0o+DGwH6tNsq7bzxGKmYwJZnica0eeJMeaoMebpxN+DxP/nLCrtqEqLiTOUeOhO/Js8638V8KPE37cDF4mIFGmIRSfDY1J1iEgH8Hrg+9PsUlXnCWR0TLJGhb6AJG4rNwGPp9l8duK2/R4ROamoAysBiVvPZ4ETwB+MMZOPyYSG8oDVUN62ZHBMAN4kIs+LyO0isjjNdrvxTeAGIDbN9qo7T5j9mECW54kKfYEQkVrgV8BHjDEDkzY/Tbw+xanAt4E7iz2+YmOMiRpjTiPeY/hMETm51GMqNRkck98Cy4wxG4E/MB7J2hIRuQI4YYx5qtRjKRcyPCZZnycq9AUg4bn+CvipMeaOyduNMQPWbXuiW5dbRFqLPMySYIzpAx4ELpu0qWobyk93TIwx3caYYOLh94Ezij22InMOcKWI7AN+DlwoIj+ZtE+1nSezHpNczhMV+jxJ+IU/ALYbY74+zT7tlq8oImcSP+62PVlFpE1EGhN/+4FLgB2TdquqhvKZHBMRWZDy8Eri8z22xRjzKWNMhzFmGfE+0w8YY94xabeqOk8yOSa5nCeadZM/5wDvBF5I+K8AnwaWABhj/oP4CfqPIhIBRoFr7HyyAguAH4mIk/hF7TZjzP+t8obymRyTD4nIlcQzuXqAd5dstCWkys+TtOR7nmgJBEVRFJuj1o2iKIrNUaFXFEWxOSr0iqIoNkeFXlEUxeao0CuKohQQEfkvETkhIi9msO+rReRpEYmIyJsnbfu9iPSlK2yWLSr0iqIoheWHTF0gOB0HiKdH/izNtq8ST93OGxV6RVGUAmKM+TPx/PYkIrIyEaE/JSIPi8i6xL77jDHPk6aujTHmfmCwEGPSBVOKoihzzy3A+40xL4vIWcB3gAuL9eEq9IqiKHNIouDhq4BfplRY9hZzDCr0iqIoc4sD6EtULi3ZABRFUZQ5IlG2fK+IvAXihRBF5NRijkFr3SiKohQQEbkVOB9oBY4D/x/wAPBd4sXt3MDPjTGfE5FXAL8GmoAx4Jgx5qTE+zwMrANqiVe7vc4Yc29OY1KhVxRFsTdq3SiKotgcFXpFURSbo0KvKIpic1ToFUVRbI4KvaIois1RoVcURbE5KvSKoig25/8BE+gbQ+w28+MAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
}
}
]
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment