{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Functional form\n", "\n", "This notebook demonstrates the use of the `FunctionalForm` model to describe an exponential increase/decay from one layer to the next. The code in `functional.py` can be used as a starting point to change the exponential profile to a different analytic form." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# necessary imports\n", "%matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "from refnx.reflect import Slab, SLD, Structure\n", "from functional import FunctionalForm" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAaoElEQVR4nO3dfZAcd33n8fd3n59XWj1ZD5bWsg3Eh20QC7HjhyImpsxDcaTiq0vqSEgCqO6Ku5C6cDl8oVKV5KhL6qog5C7F4UoukAQCOWzOiZMAxoAPH4dsSX6WbCzLerSkXe3TzOw8z/zuj+5d7a53tT2rme6ens+raunph5357o/1R7/9dfevzTmHiIjEV1vUBYiIyOUpqEVEYk5BLSIScwpqEZGYU1CLiMRcRyPedPPmzW50dLQRby0ikkiHDh266JzbstK+hgT16OgoBw8ebMRbi4gkkpmdXG2fhj5ERGJOQS0iEnMKahGRmFNQi4jEnIJaRCTmFNQiIjGnoBYRibmGXEct0kxOT2X50SsXuZgpUq44KtXqwr43XjXE+27aHmF1Eju5GTjxQ5h8BaolqFa8L4Cufrj9N+v+kQpqaWkPHDrDf3zgWcrVpfOym3nLD9y8Q0Etl5x/Hv7qgzA3scJOg4GtCmqRenptJsd9Dz7H2OhGPvPzN7JzQy+d7W20t1nUpUkcOQcPfBSsHT78MOx4K3T0QFv7pX/ZG0RBLS3rqwdOUapW+a/33szVI31RlyNxd+KHMHEUfv6LcM0doX60TiZKy/rWC+e5/brNCmkJ5shD0DUAN/zz0D9aQS0taSZb5Nh4hp++ZiTqUqRZnD4Au8agszf0j1ZQS0s6fGoagLftUVBLAIU0XHgBrr4lko9XUEtLOnouDcCNu4YjrkSawsRL4Kqw/aZIPl5BLS3p1GSWzQPdDHTrfLoEMHXcW45cG8nHK6ilJZ2YnGN0k04iSkCTrwAGG0cj+XgFtbSkk5NZ9mzqj7oMaRZTr8Dw1dDZE8nHK6il5eRLFc6n8uxRj1qCmjkFG/dE9vEKamk546kCANuHo+kdSRNKn4PB6KYSUFBLy5nI5AHYMtgdcSXSFJyD9AUY3BZZCYGC2sw2mNk3zOxFMztqZrc2ujCRRplIez1qBbUEkp+BSgEGroqshKDXJn0e+JZz7l4z6wI0uCdNS0EtNUlf8JaDMQ5qMxsG7gR+FcA5VwSKjS1LpHEm0gXaDDb1K6glgMx5bzkQ76GPa4AJ4C/M7Ckz+zMz03VN0rQmMkVG+rs1nakEkxn3ljEP6g5gH/AF59xbgTngU8sPMrP9ZnbQzA5OTKw0qbZIPEykCxr2kOCyk96yf3NkJQQJ6jPAGefcAX/9G3jBvYRz7n7n3JhzbmzLli31rFGkrqazRUb6O6MuQ5pFbhow6IluXpg1g9o5dx44bWZv9De9CzjS0KpEGmgmW2RDb1fUZUizyE55Id3WHlkJQa/6+HfAV/wrPo4Dv9a4kkQaazZXZrhPPWoJKDcNfdFOhxsoqJ1zTwNjjS1FpPGcc8zmigz3KqgloNwU9G6MtATdmSgtJVusUKo4NiioJajsFPRG26NWUEtLmcmVANigoQ8JKgZDHwpqaSmzWS+oNfQhgeWmNfQhEqaZnHdT7bCu+pAgqhUopKBnQ6RlKKilpcz3qDX0IYEUvGdr0jMUaRkKamkps/4Y9ZCGPiSIQspbdiuoRUKTKZQBGOzRQ20lgLwf1OpRi4RnPqj7uxTUEoB61CLhy+TL9HW1a+Y8CUY9apHwzRXL9HerNy0BLfSoo5uQCRTU0mLS+TKDCmoJKj/rLdWjFgnPXKHMgE4kSlDzQa0xapHwZAplnUiU4AopaO+Czp5Iy1BQS0tJ59WjlhoU0pH3pkFBLS1mrlhmQGPUElQhA90DUVehoJbWkskrqKUGxQx0DUZdhYJaWkumoMvzpAaFtHrUImEqlL2HBuj2cQmsmIEuBbVIaDJ57/ZxDX1IYBqjFgnXXKECoKEPCU49apFwpQveFKfqUUtghQx062SiSGg09CE1cc7vUfdHXYmCWlrHXNEPap1MlCCKc4DT0IdImNLqUUstihlvqZOJIuGZP5mooJZACn5Q64YXkfBk5k8mauhDgij6D7ZVj1okPPMnE/s62yOuRJrCQo9aQS0SmkyhwkB3B216DJcEoTFqkfBlCiWNT0twGqMWCZ83IZOGPSQgjVGLhG9+6EMkEI1Ri4Qvky/pig8JrqigFgndXKGi5yVKcIUMdPZDW/QxGX0FIiHJFPR0F6lBMR4PDQAFtbSQuaIebCs1KMRjilNQUEuLcM6RyesxXFKDYjweGgAKamkRhXKVctVp6EOCK8TjwbYAgX5rzewEkAYqQNk5N9bIokTqba6gmfOkRsU0DO2MugogYFD7ftY5d7FhlYg0UMYPag19SGAxeboLaOhDWkRmoUetOxMloEK66YLaAd8xs0Nmtn+lA8xsv5kdNLODExMT9atQpA70YFupWRMG9e3OuX3Ae4CPm9mdyw9wzt3vnBtzzo1t2bKlrkWKXKmMHmwrtSgXoFJorqB2zp31l+PAN4F3NLIokXrL6OkuUov5eT66h6Ktw7dmUJtZv5kNzr8G3g083+jCROppTicTpRaFlLeMSY86yG/tNuCbZjZ//Fedc99qaFUidTb/dBfdmSiBFOanOG2SoHbOHQduDqEWkYZZuDxPkzJJEDELal2eJy1hrlCmt7Oddj2GS4KYD+qY3JmooJaWoAmZpCbqUYuEL53XFKdSg5idTFRQS0uY0/MSpRbqUYuEb07PS5RaFDOAQVd/1JUACmppEWk93UVqUUh7N7tYPE4+K6ilJXhDHwpqCShG83yAglpahIJaalJIKahFwpYulBlUUEtQ6lGLhKtYrlIsV9WjluAU1CLhSuW9KU6HezsjrkSahoJaJFyzOQW11EhBLRKulB/UQ70a+pCA5i/PiwkFtSSeetRSk2rFu+GleyDqShYoqCXxFNRSk/yst+zdGG0diyioJfFS/kMDhnoU1BJAbtpbKqhFwnNpjFpBLQEoqEXCl8qV6O5oo6dTs+dJAApqkfDN5krqTUtwC0E9Em0diyioJfFmsiWdSJTgslPeUj1qkfBMzRXZ1N8VdRnSLOZ71D3D0daxiIJaEu/iXIHNA91RlyHNIjfthXR7fG6QUlBL4k1mimwaUI9aAspOxmp8GhTUknDFcpXZXIlN/epRS0CZCzCwLeoqllBQS6JNZ4sA6lFLcJlxGNgadRVLKKgl0S5mCgBsVlBLUJnzMHhV1FUsoaCWRJvMzPeoNfQhAZTy3lwf6lGLhOdCKg/AFgW1BDE37i01Ri0SnjPTOcxgx4beqEuRZpC+4C0V1CLhOTOd46qhHro69KsuAcyc9JYbdkdbxzL67ZVEOzOdZddG9aYloKnj3nLjaKRlLKeglkQ7M51jp4Y9JKipV2FwB3TG63dGQS2JlS9VODebY/dIX9SlSLOYfhVGrom6itdRUEtiHTmXourghh3xmVxHYqxahfEjsPn6qCt5HQW1JNbzZ71n3924S0EtAUwe866h3jkWdSWvo6CWxHrq1Awj/V3sGO6JuhRpBmee9Ja73h5tHSsIHNRm1m5mT5nZw40sSKQe8qUKjxy5wF1v2oqZRV2ONIMjD8HQTtj8hqgreZ1aetSfAI42qhCRevrqgVNkCmU++JadUZcizeDcs3Dsu3DTv4S2+A00BJoZ28x2Ae8DPgP8+4ZWJA3z7JkZDp+cJlMoU6w4cC7qktbFAVXnKFcd1aq3rFQvrU/NFfnuUa83fdt1m6ItdvYMlAt+W/vtvdDul1uvx7GsfayrQjkHxbmlX6UsVCvr+IFrsbzuNbYt2b7SttW+/zLHARRScPTvoX8z3PpvA1UetqCPMPhj4LeBwdUOMLP9wH6A3bvjdVdPq3POcd+Dz/G1J08v2d7MIwJtZrS3GR1t3nLx657Odj78M6P81rvfGHzYwzkvnOrp0Jfg2/+pvu8ZqhB+QRb+/7GA2xZtX2nbet6zoxtG74C7fx/6I/6HfRVrBrWZvR8Yd84dMrN3rnacc+5+4H6AsbGx5uyqJdT/fvosX3vyNB+74xo+esdeNvZ10dluGrtd7OsfghcbcPpl9A7Y9yv+ykoBM7++fN8K65fbV9P7sHRfVz909nnLrgHo6vPW29pr+EGlkYL0qG8DPmBm7wV6gCEz+2vn3IcaW5rUy5d+dJI3bBvgvvf8FG1tLRDOEz/x/pwNqjgHL/0jvOEe2PMz9aujrRNu/kXoi9djnaT5rBnUzrn7gPsA/B71JxXSzePUZJZnTs/w6fe1SEhfOAJfuHV93/uu34Vt/6y+9YjUQXwesysNcfDkFAC3X7854kpC8upj3vIX/tx7knRQvSMKaYmtmoLaOfcD4AcNqUQa4vCpaQa6O7h+66rngeOvXPBOyuWm1z727GEY3g033tv4ukRCoh51wj13Zpabdg3T3szDHmeehCf/DIav9s7QX05bO7ztV0MpSyQsCuoEc85x/OJc89/0cf45b/nRR2EwXk/eEAmDgjrBZrIl0vkyezbFfJrP6ZNw6ser73/pn6B/q0JaWpaCOsFOTM4BMLqpP+JK1vCPn4SXv3P5Y970/nBqEYkhBXWCnZz07rQb3RzzHvX4i14Q3/37qx8zvCu8ekRiRkGdYGdncgDs2hjjoC7lYPY07Ptl2HRt1NWIxJKCOsHGU3mGejro6YzBrcClHBx7FKqlpdtT5wCnkBa5DAV1go2nC2wdismk+U9/Ff7hMhMvXnVzeLWINBkFdYKNpwtsHVzjuuOwTB2Hjl742PdePylR96DGoEUuQ0GdYOPpPPt2b4y6DM/MSdiwG7bdEHUlIk0nfo8ykLpwzjGeilGPeuaUF9QiUjP1qBMqlS9TKFfZOhjyGHX6ArzwTe/JIYtNHo/l051FmoGCOqFmskUARvq7wv3gA/8DHv/syvt27gu3FpGEUFAnVCpXBmCotzPkDz7rTZ70rx9fut3aoGco3FpEEkJBnVCzOe965aGekP8vTr0GQzugd0O4nyuSYDqZmFCpvBfUw30h96jT52Fwe7ifKZJwCuqESi30qBXUIs1OQx8JNd+jbtgY9fnn4bm/BbfogfOuCsU0DF7VmM8UaVEK6oSazZVobzP6uxo0z8eP/hs8+zXvbsPFuodh19sb85kiLUpBnVCpXJmhng5s+e3a9ZK54F0X/bFHG/P+IrJAY9QJlcqXGntp3twEDGxt3PuLyAIFdUKlcqXGnkjMjEP/lsa9v4gsUFAn1GyuxFBvg0a2qhXIXlSPWiQkCuqESuXLDDdq6CM76V3h0a+gFgmDTiYmVF2HPk4/Ad/7A6j6Ey2VvGcxMqChD5EwqEedUHU9mfjiw3Bi0dwdnX1w3d2w+9b6vL+IXJZ61AlUKFfIl6r1G/qYm4SBbfBr/1Cf9xORmqhHnUALM+fVa0Km7EXo21yf9xKRmimoE6jut4/PXYT+TfV5LxGpmYI6geo+IZN61CKRUlAn0MJc1PUco+5XUItERScTEyiV98aoh9d7w8v0SXjtsPe6WvFmxFOPWiQyCuoEuuKhj4c+Did+uHTbyDVXWJWIrJeCOoGu+GRi+hxc93Pw7v/srbd3wcjeOlUnIrVSUCfQbK5EV0cbPZ3rnIs6Owl73wlbf6qudYnI+uhkYgJ5c1GvszddKUNuWmPSIjGyZlCbWY+ZPWFmz5jZC2b2e2EUJuvn3T6+zj+WctPesk/XTYvERZD/mgvAXc65jJl1Ao+b2T85537c4Npkna5oQqbspLfUDS4isbFmj9p5Mv5qp//lLvMtErFUrrT+eT6yF72letQisRHo72MzawcOAdcBf+qcO9DQquSKpPJldm/qr+2bDv8lvPyI9yxEUFCLxEigk4nOuYpz7i3ALuAdZvbm5ceY2X4zO2hmBycmJupcptTCG/qocYz6//4JHH8M8ikYvQNGrm1McSJSs5r+a3bOzZjZ94F7gOeX7bsfuB9gbGxMQyMRcc4xu56hj9wU3PgL8P7PNaYwEVm3IFd9bDGzDf7rXuBu4MUG1yXrlCtVKFddbTe7VKve1R69GxtXmIisW5Ae9Xbgy/44dRvwt865hxtblqzXpbmoawjqQsp7BmLvSIOqEpErsWZQO+eeBd4aQi1SB5duH69hVCs35S37FNQicaQ7ExNmforTmsaos/5NLupRi8SSgjph1jVz3sLdiApqkTjSpEwJE3jmvEIGjn0XKiU484S3TScTRWJJQZ0wgR9se+gv4DufvrTe0QODVzWwMhFZLwV1wgR+DFdm3Jtn+t/8CKzN6013D4ZQoYjUSkGdMKlcib6udjrb1zj9UEhBzwbYfH0odYnI+ulkYsKk8gFnzsvPQs9w4wsSkSumoE6YVK4c7BpqBbVI01BQJ8xs0Lmo87PQM9T4gkTkiimoE2YmV2JDX9faB6pHLdI0FNQJM5stsqFPY9QiSaKgTpjpbIkNQW4fz6cU1CJNQpfnJUi+VCFXqrCxf5Whj3PPwMRLkHoNKgUFtUiTUFAnSOpyEzKlz8MX77y03tYB224MqTIRuRIK6gSZ8YN6xTHq+Wchvu+zcOO/gK4BaNPIl0gzUFAnyEzWD+reFYY+8ilvuek6XZYn0mTUpUqQmWwRWKVHXUh7S4W0SNNRUCfIfI96xTHqgt+j7lZQizQbBXWCzOQu06POK6hFmpWCOkFmsiU62oyB7hVOPRRmvaWGPkSajoI6Qbzbxzsxs9fvzKegvRs6usMvTESuiII6QWazpdUfaltIqTct0qQU1AkyOVdg42oTMuVTGp8WaVIK6gQZTxfYOrRsaKOYhUNfhhceVI9apEkpqBNkIlVg62DP0o2P/RH8/W/Axmtg7NejKUxErojuTEyIXLFCulBmy+CyHvXcBAzthN94ClY6ySgisacedUKMp/MAbF0e1PPzTiukRZqWgjohxtMFALYOLRv6KKShezCCikSkXhTUCTGe8oN6eY+6kNbVHiJNTkGdEGdnsgDsGO5duqOQUo9apMkpqBPixGSWjX2dDC+f50NDHyJNT0GdECcn59i9qf/1OwppXT8t0uQU1Alx4mKW0U19SzdWSlDKaoxapMkpqBMgX6rw2myOPct71PMPC1BQizQ1BXUCPHd2FufgzTuWBXJ20ltqjFqkqcUqqCvlMq5ajbqMpnP45DQA+/ZsvLTx+A/gyx/wXm+6NvyiRKRuYnULeeEPdtBnBarOqNBGlbaFZdXaqNDOxY6tTI68ja13/jp733xL1CXHwoFXp9g90sfmgW6Ym4Tv/A488zcwshc+9ADsVjuJNLM1g9rMrgb+EtgGOOB+59znG1HMM3s/CuUiuApUK5irLnpdgWqZ/swJ9l14gK5vfJ0nv/1z7Lr3v7B9zxsaUU5TGE/leewnE+y/7Wr48RfgB38IxTm48z/AHZ+Ezp6130REYs2cc5c/wGw7sN05d9jMBoFDwAedc0dW+56xsTF38ODB+la6yOzUOC8+8BluPvMVDDi89YNsueMj7L3h7bR1xOqPhMZxjlxqii9+/UG6Tj/O/qEf0zF3Hvb+LNzzh7D1TVFXKCI1MLNDzrmxFfetFdQrvNlDwH93zj2y2jGNDup5F069zNlvfpobpx6h0yoUXCdp66dg3ZTpxBkY85MRuUX/u/hFMEEPN789DYfhFr7TlmxjyTG26BNW38aS9+uhSB/ebeNV2mi79p1wy8fhundpAiaRJlS3oDazUeD/AG92zqWW7dsP7AfYvXv3206ePLnugmt18cJZXv1/D1K9cBSXT9FRydHuyjjncMCqseUH2mr7X7fdlr5Yun/xJxkYC/E6/1leS7d5S7t0rPOPZ1Ece/uXf/98fBuVtk6yPVex6/qbeePb79ZNLSJNri5BbWYDwGPAZ5xzD17u2LB61CIiSXG5oA50eZ6ZdQIPAF9ZK6RFRKS+1gxqMzPgz4GjzrnPNr4kERFZLEiP+jbgl4G7zOxp/+u9Da5LRER8a17L5px7nMucjxMRkcaK1S3kIiLyegpqEZGYU1CLiMScglpEJOZqvoU80JuaTQBh3Zq4GbgY0mc1E7XLytQuK1O7rCzMdtnjnNuy0o6GBHWYzOzganfztDK1y8rULitTu6wsLu2ioQ8RkZhTUIuIxFwSgvr+qAuIKbXLytQuK1O7rCwW7dL0Y9QiIkmXhB61iEiiKahFRGKuaYPazO4xs5fM7JiZfSrqesJmZv/TzMbN7PlF20bM7BEze9lfbvS3m5n9id9Wz5rZvugqbxwzu9rMvm9mR8zsBTP7hL+91dulx8yeMLNn/Hb5PX/7NWZ2wP/5v25mXf72bn/9mL9/NNIfoMHMrN3MnjKzh/312LVLUwa1mbUDfwq8B7gB+CUzuyHaqkL3JeCeZds+BTzqnLseeNRfB6+drve/9gNfCKnGsJWB33LO3QDcAnzc/71o9XYpAHc5524G3gLcY2a3AH8EfM45dx0wDXzEP/4jwLS//XP+cUn2CeDoovX4tYtzrum+gFuBby9avw+4L+q6ImiHUeD5Resv4T0xHmA78JL/+ovAL610XJK/gIeAu9UuS9qkDzgM/DTeHXcd/vaF/6aAbwO3+q87/OMs6tob1B678P7xvgt4GG9K59i1S1P2qIGdwOlF62f8ba1um3PunP/6PLDNf91y7eX/WfpW4ABql/k/758GxoFHgFeAGedc2T9k8c++0C7+/llgU6gFh+ePgd8Gqv76JmLYLs0a1LIG5/2z35LXXvoPYn4A+E3nXGrxvlZtF+dcxTn3Frwe5DuAN0VbUfTM7P3AuHPuUNS1rKVZg/oscPWi9V3+tlZ3wcy2A/jLcX97y7TXKg9ibvl2meecmwG+j/cn/QYzm3/K0+KffaFd/P3DwGS4lYbiNuADZnYC+Bre8MfniWG7NGtQPwlc75+d7QJ+Efi7iGuKg78DPuy//jDeGO389l/xr3K4BZhdNBSQGJd5EHOrt8sWM9vgv+7FG7c/ihfY9/qHLW+X+fa6F/ie/5dIojjn7nPO7XLOjeJlyPecc/+KOLZL1IP5V3AS4L3AT/DG2n4n6noi+Pn/BjgHlPDG0T6CN172KPAy8F1gxD/W8K6SeQV4DhiLuv4GtcnteMMazwJP+1/vVbtwE/CU3y7PA7/rb98LPAEcA/4X0O1v7/HXj/n790b9M4TQRu8EHo5ru+gWchGRmGvWoQ8RkZahoBYRiTkFtYhIzCmoRURiTkEtIhJzCmoRkZhTUIuIxNz/B8Kpm7j5pBYiAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAaoElEQVR4nO3dfZAcd33n8fd3n59XWj1ZD5bWsg3Eh20QC7HjhyImpsxDcaTiq0vqSEgCqO6Ku5C6cDl8oVKV5KhL6qog5C7F4UoukAQCOWzOiZMAxoAPH4dsSX6WbCzLerSkXe3TzOw8z/zuj+5d7a53tT2rme6ens+raunph5357o/1R7/9dfevzTmHiIjEV1vUBYiIyOUpqEVEYk5BLSIScwpqEZGYU1CLiMRcRyPedPPmzW50dLQRby0ikkiHDh266JzbstK+hgT16OgoBw8ebMRbi4gkkpmdXG2fhj5ERGJOQS0iEnMKahGRmFNQi4jEnIJaRCTmFNQiIjGnoBYRibmGXEct0kxOT2X50SsXuZgpUq44KtXqwr43XjXE+27aHmF1Eju5GTjxQ5h8BaolqFa8L4Cufrj9N+v+kQpqaWkPHDrDf3zgWcrVpfOym3nLD9y8Q0Etl5x/Hv7qgzA3scJOg4GtCmqRenptJsd9Dz7H2OhGPvPzN7JzQy+d7W20t1nUpUkcOQcPfBSsHT78MOx4K3T0QFv7pX/ZG0RBLS3rqwdOUapW+a/33szVI31RlyNxd+KHMHEUfv6LcM0doX60TiZKy/rWC+e5/brNCmkJ5shD0DUAN/zz0D9aQS0taSZb5Nh4hp++ZiTqUqRZnD4Au8agszf0j1ZQS0s6fGoagLftUVBLAIU0XHgBrr4lko9XUEtLOnouDcCNu4YjrkSawsRL4Kqw/aZIPl5BLS3p1GSWzQPdDHTrfLoEMHXcW45cG8nHK6ilJZ2YnGN0k04iSkCTrwAGG0cj+XgFtbSkk5NZ9mzqj7oMaRZTr8Dw1dDZE8nHK6il5eRLFc6n8uxRj1qCmjkFG/dE9vEKamk546kCANuHo+kdSRNKn4PB6KYSUFBLy5nI5AHYMtgdcSXSFJyD9AUY3BZZCYGC2sw2mNk3zOxFMztqZrc2ujCRRplIez1qBbUEkp+BSgEGroqshKDXJn0e+JZz7l4z6wI0uCdNS0EtNUlf8JaDMQ5qMxsG7gR+FcA5VwSKjS1LpHEm0gXaDDb1K6glgMx5bzkQ76GPa4AJ4C/M7Ckz+zMz03VN0rQmMkVG+rs1nakEkxn3ljEP6g5gH/AF59xbgTngU8sPMrP9ZnbQzA5OTKw0qbZIPEykCxr2kOCyk96yf3NkJQQJ6jPAGefcAX/9G3jBvYRz7n7n3JhzbmzLli31rFGkrqazRUb6O6MuQ5pFbhow6IluXpg1g9o5dx44bWZv9De9CzjS0KpEGmgmW2RDb1fUZUizyE55Id3WHlkJQa/6+HfAV/wrPo4Dv9a4kkQaazZXZrhPPWoJKDcNfdFOhxsoqJ1zTwNjjS1FpPGcc8zmigz3KqgloNwU9G6MtATdmSgtJVusUKo4NiioJajsFPRG26NWUEtLmcmVANigoQ8JKgZDHwpqaSmzWS+oNfQhgeWmNfQhEqaZnHdT7bCu+pAgqhUopKBnQ6RlKKilpcz3qDX0IYEUvGdr0jMUaRkKamkps/4Y9ZCGPiSIQspbdiuoRUKTKZQBGOzRQ20lgLwf1OpRi4RnPqj7uxTUEoB61CLhy+TL9HW1a+Y8CUY9apHwzRXL9HerNy0BLfSoo5uQCRTU0mLS+TKDCmoJKj/rLdWjFgnPXKHMgE4kSlDzQa0xapHwZAplnUiU4AopaO+Czp5Iy1BQS0tJ59WjlhoU0pH3pkFBLS1mrlhmQGPUElQhA90DUVehoJbWkskrqKUGxQx0DUZdhYJaWkumoMvzpAaFtHrUImEqlL2HBuj2cQmsmIEuBbVIaDJ57/ZxDX1IYBqjFgnXXKECoKEPCU49apFwpQveFKfqUUtghQx062SiSGg09CE1cc7vUfdHXYmCWlrHXNEPap1MlCCKc4DT0IdImNLqUUstihlvqZOJIuGZP5mooJZACn5Q64YXkfBk5k8mauhDgij6D7ZVj1okPPMnE/s62yOuRJrCQo9aQS0SmkyhwkB3B216DJcEoTFqkfBlCiWNT0twGqMWCZ83IZOGPSQgjVGLhG9+6EMkEI1Ri4Qvky/pig8JrqigFgndXKGi5yVKcIUMdPZDW/QxGX0FIiHJFPR0F6lBMR4PDQAFtbSQuaIebCs1KMRjilNQUEuLcM6RyesxXFKDYjweGgAKamkRhXKVctVp6EOCK8TjwbYAgX5rzewEkAYqQNk5N9bIokTqba6gmfOkRsU0DO2MugogYFD7ftY5d7FhlYg0UMYPag19SGAxeboLaOhDWkRmoUetOxMloEK66YLaAd8xs0Nmtn+lA8xsv5kdNLODExMT9atQpA70YFupWRMG9e3OuX3Ae4CPm9mdyw9wzt3vnBtzzo1t2bKlrkWKXKmMHmwrtSgXoFJorqB2zp31l+PAN4F3NLIokXrL6OkuUov5eT66h6Ktw7dmUJtZv5kNzr8G3g083+jCROppTicTpRaFlLeMSY86yG/tNuCbZjZ//Fedc99qaFUidTb/dBfdmSiBFOanOG2SoHbOHQduDqEWkYZZuDxPkzJJEDELal2eJy1hrlCmt7Oddj2GS4KYD+qY3JmooJaWoAmZpCbqUYuEL53XFKdSg5idTFRQS0uY0/MSpRbqUYuEb07PS5RaFDOAQVd/1JUACmppEWk93UVqUUh7N7tYPE4+K6ilJXhDHwpqCShG83yAglpahIJaalJIKahFwpYulBlUUEtQ6lGLhKtYrlIsV9WjluAU1CLhSuW9KU6HezsjrkSahoJaJFyzOQW11EhBLRKulB/UQ70a+pCA5i/PiwkFtSSeetRSk2rFu+GleyDqShYoqCXxFNRSk/yst+zdGG0diyioJfFS/kMDhnoU1BJAbtpbKqhFwnNpjFpBLQEoqEXCl8qV6O5oo6dTs+dJAApqkfDN5krqTUtwC0E9Em0diyioJfFmsiWdSJTgslPeUj1qkfBMzRXZ1N8VdRnSLOZ71D3D0daxiIJaEu/iXIHNA91RlyHNIjfthXR7fG6QUlBL4k1mimwaUI9aAspOxmp8GhTUknDFcpXZXIlN/epRS0CZCzCwLeoqllBQS6JNZ4sA6lFLcJlxGNgadRVLKKgl0S5mCgBsVlBLUJnzMHhV1FUsoaCWRJvMzPeoNfQhAZTy3lwf6lGLhOdCKg/AFgW1BDE37i01Ri0SnjPTOcxgx4beqEuRZpC+4C0V1CLhOTOd46qhHro69KsuAcyc9JYbdkdbxzL67ZVEOzOdZddG9aYloKnj3nLjaKRlLKeglkQ7M51jp4Y9JKipV2FwB3TG63dGQS2JlS9VODebY/dIX9SlSLOYfhVGrom6itdRUEtiHTmXourghh3xmVxHYqxahfEjsPn6qCt5HQW1JNbzZ71n3924S0EtAUwe866h3jkWdSWvo6CWxHrq1Awj/V3sGO6JuhRpBmee9Ja73h5tHSsIHNRm1m5mT5nZw40sSKQe8qUKjxy5wF1v2oqZRV2ONIMjD8HQTtj8hqgreZ1aetSfAI42qhCRevrqgVNkCmU++JadUZcizeDcs3Dsu3DTv4S2+A00BJoZ28x2Ae8DPgP8+4ZWJA3z7JkZDp+cJlMoU6w4cC7qktbFAVXnKFcd1aq3rFQvrU/NFfnuUa83fdt1m6ItdvYMlAt+W/vtvdDul1uvx7GsfayrQjkHxbmlX6UsVCvr+IFrsbzuNbYt2b7SttW+/zLHARRScPTvoX8z3PpvA1UetqCPMPhj4LeBwdUOMLP9wH6A3bvjdVdPq3POcd+Dz/G1J08v2d7MIwJtZrS3GR1t3nLx657Odj78M6P81rvfGHzYwzkvnOrp0Jfg2/+pvu8ZqhB+QRb+/7GA2xZtX2nbet6zoxtG74C7fx/6I/6HfRVrBrWZvR8Yd84dMrN3rnacc+5+4H6AsbGx5uyqJdT/fvosX3vyNB+74xo+esdeNvZ10dluGrtd7OsfghcbcPpl9A7Y9yv+ykoBM7++fN8K65fbV9P7sHRfVz909nnLrgHo6vPW29pr+EGlkYL0qG8DPmBm7wV6gCEz+2vn3IcaW5rUy5d+dJI3bBvgvvf8FG1tLRDOEz/x/pwNqjgHL/0jvOEe2PMz9aujrRNu/kXoi9djnaT5rBnUzrn7gPsA/B71JxXSzePUZJZnTs/w6fe1SEhfOAJfuHV93/uu34Vt/6y+9YjUQXwesysNcfDkFAC3X7854kpC8upj3vIX/tx7knRQvSMKaYmtmoLaOfcD4AcNqUQa4vCpaQa6O7h+66rngeOvXPBOyuWm1z727GEY3g033tv4ukRCoh51wj13Zpabdg3T3szDHmeehCf/DIav9s7QX05bO7ztV0MpSyQsCuoEc85x/OJc89/0cf45b/nRR2EwXk/eEAmDgjrBZrIl0vkyezbFfJrP6ZNw6ser73/pn6B/q0JaWpaCOsFOTM4BMLqpP+JK1vCPn4SXv3P5Y970/nBqEYkhBXWCnZz07rQb3RzzHvX4i14Q3/37qx8zvCu8ekRiRkGdYGdncgDs2hjjoC7lYPY07Ptl2HRt1NWIxJKCOsHGU3mGejro6YzBrcClHBx7FKqlpdtT5wCnkBa5DAV1go2nC2wdismk+U9/Ff7hMhMvXnVzeLWINBkFdYKNpwtsHVzjuuOwTB2Hjl742PdePylR96DGoEUuQ0GdYOPpPPt2b4y6DM/MSdiwG7bdEHUlIk0nfo8ykLpwzjGeilGPeuaUF9QiUjP1qBMqlS9TKFfZOhjyGHX6ArzwTe/JIYtNHo/l051FmoGCOqFmskUARvq7wv3gA/8DHv/syvt27gu3FpGEUFAnVCpXBmCotzPkDz7rTZ70rx9fut3aoGco3FpEEkJBnVCzOe965aGekP8vTr0GQzugd0O4nyuSYDqZmFCpvBfUw30h96jT52Fwe7ifKZJwCuqESi30qBXUIs1OQx8JNd+jbtgY9fnn4bm/BbfogfOuCsU0DF7VmM8UaVEK6oSazZVobzP6uxo0z8eP/hs8+zXvbsPFuodh19sb85kiLUpBnVCpXJmhng5s+e3a9ZK54F0X/bFHG/P+IrJAY9QJlcqXGntp3twEDGxt3PuLyAIFdUKlcqXGnkjMjEP/lsa9v4gsUFAn1GyuxFBvg0a2qhXIXlSPWiQkCuqESuXLDDdq6CM76V3h0a+gFgmDTiYmVF2HPk4/Ad/7A6j6Ey2VvGcxMqChD5EwqEedUHU9mfjiw3Bi0dwdnX1w3d2w+9b6vL+IXJZ61AlUKFfIl6r1G/qYm4SBbfBr/1Cf9xORmqhHnUALM+fVa0Km7EXo21yf9xKRmimoE6jut4/PXYT+TfV5LxGpmYI6geo+IZN61CKRUlAn0MJc1PUco+5XUItERScTEyiV98aoh9d7w8v0SXjtsPe6WvFmxFOPWiQyCuoEuuKhj4c+Did+uHTbyDVXWJWIrJeCOoGu+GRi+hxc93Pw7v/srbd3wcjeOlUnIrVSUCfQbK5EV0cbPZ3rnIs6Owl73wlbf6qudYnI+uhkYgJ5c1GvszddKUNuWmPSIjGyZlCbWY+ZPWFmz5jZC2b2e2EUJuvn3T6+zj+WctPesk/XTYvERZD/mgvAXc65jJl1Ao+b2T85537c4Npkna5oQqbspLfUDS4isbFmj9p5Mv5qp//lLvMtErFUrrT+eT6yF72letQisRHo72MzawcOAdcBf+qcO9DQquSKpPJldm/qr+2bDv8lvPyI9yxEUFCLxEigk4nOuYpz7i3ALuAdZvbm5ceY2X4zO2hmBycmJupcptTCG/qocYz6//4JHH8M8ikYvQNGrm1McSJSs5r+a3bOzZjZ94F7gOeX7bsfuB9gbGxMQyMRcc4xu56hj9wU3PgL8P7PNaYwEVm3IFd9bDGzDf7rXuBu4MUG1yXrlCtVKFddbTe7VKve1R69GxtXmIisW5Ae9Xbgy/44dRvwt865hxtblqzXpbmoawjqQsp7BmLvSIOqEpErsWZQO+eeBd4aQi1SB5duH69hVCs35S37FNQicaQ7ExNmforTmsaos/5NLupRi8SSgjph1jVz3sLdiApqkTjSpEwJE3jmvEIGjn0XKiU484S3TScTRWJJQZ0wgR9se+gv4DufvrTe0QODVzWwMhFZLwV1wgR+DFdm3Jtn+t/8CKzN6013D4ZQoYjUSkGdMKlcib6udjrb1zj9UEhBzwbYfH0odYnI+ulkYsKk8gFnzsvPQs9w4wsSkSumoE6YVK4c7BpqBbVI01BQJ8xs0Lmo87PQM9T4gkTkiimoE2YmV2JDX9faB6pHLdI0FNQJM5stsqFPY9QiSaKgTpjpbIkNQW4fz6cU1CJNQpfnJUi+VCFXqrCxf5Whj3PPwMRLkHoNKgUFtUiTUFAnSOpyEzKlz8MX77y03tYB224MqTIRuRIK6gSZ8YN6xTHq+Wchvu+zcOO/gK4BaNPIl0gzUFAnyEzWD+reFYY+8ilvuek6XZYn0mTUpUqQmWwRWKVHXUh7S4W0SNNRUCfIfI96xTHqgt+j7lZQizQbBXWCzOQu06POK6hFmpWCOkFmsiU62oyB7hVOPRRmvaWGPkSajoI6Qbzbxzsxs9fvzKegvRs6usMvTESuiII6QWazpdUfaltIqTct0qQU1AkyOVdg42oTMuVTGp8WaVIK6gQZTxfYOrRsaKOYhUNfhhceVI9apEkpqBNkIlVg62DP0o2P/RH8/W/Axmtg7NejKUxErojuTEyIXLFCulBmy+CyHvXcBAzthN94ClY6ySgisacedUKMp/MAbF0e1PPzTiukRZqWgjohxtMFALYOLRv6KKShezCCikSkXhTUCTGe8oN6eY+6kNbVHiJNTkGdEGdnsgDsGO5duqOQUo9apMkpqBPixGSWjX2dDC+f50NDHyJNT0GdECcn59i9qf/1OwppXT8t0uQU1Alx4mKW0U19SzdWSlDKaoxapMkpqBMgX6rw2myOPct71PMPC1BQizQ1BXUCPHd2FufgzTuWBXJ20ltqjFqkqcUqqCvlMq5ajbqMpnP45DQA+/ZsvLTx+A/gyx/wXm+6NvyiRKRuYnULeeEPdtBnBarOqNBGlbaFZdXaqNDOxY6tTI68ja13/jp733xL1CXHwoFXp9g90sfmgW6Ym4Tv/A488zcwshc+9ADsVjuJNLM1g9rMrgb+EtgGOOB+59znG1HMM3s/CuUiuApUK5irLnpdgWqZ/swJ9l14gK5vfJ0nv/1z7Lr3v7B9zxsaUU5TGE/leewnE+y/7Wr48RfgB38IxTm48z/AHZ+Ezp6130REYs2cc5c/wGw7sN05d9jMBoFDwAedc0dW+56xsTF38ODB+la6yOzUOC8+8BluPvMVDDi89YNsueMj7L3h7bR1xOqPhMZxjlxqii9+/UG6Tj/O/qEf0zF3Hvb+LNzzh7D1TVFXKCI1MLNDzrmxFfetFdQrvNlDwH93zj2y2jGNDup5F069zNlvfpobpx6h0yoUXCdp66dg3ZTpxBkY85MRuUX/u/hFMEEPN789DYfhFr7TlmxjyTG26BNW38aS9+uhSB/ebeNV2mi79p1wy8fhundpAiaRJlS3oDazUeD/AG92zqWW7dsP7AfYvXv3206ePLnugmt18cJZXv1/D1K9cBSXT9FRydHuyjjncMCqseUH2mr7X7fdlr5Yun/xJxkYC/E6/1leS7d5S7t0rPOPZ1Ece/uXf/98fBuVtk6yPVex6/qbeePb79ZNLSJNri5BbWYDwGPAZ5xzD17u2LB61CIiSXG5oA50eZ6ZdQIPAF9ZK6RFRKS+1gxqMzPgz4GjzrnPNr4kERFZLEiP+jbgl4G7zOxp/+u9Da5LRER8a17L5px7nMucjxMRkcaK1S3kIiLyegpqEZGYU1CLiMScglpEJOZqvoU80JuaTQBh3Zq4GbgY0mc1E7XLytQuK1O7rCzMdtnjnNuy0o6GBHWYzOzganfztDK1y8rULitTu6wsLu2ioQ8RkZhTUIuIxFwSgvr+qAuIKbXLytQuK1O7rCwW7dL0Y9QiIkmXhB61iEiiKahFRGKuaYPazO4xs5fM7JiZfSrqesJmZv/TzMbN7PlF20bM7BEze9lfbvS3m5n9id9Wz5rZvugqbxwzu9rMvm9mR8zsBTP7hL+91dulx8yeMLNn/Hb5PX/7NWZ2wP/5v25mXf72bn/9mL9/NNIfoMHMrN3MnjKzh/312LVLUwa1mbUDfwq8B7gB+CUzuyHaqkL3JeCeZds+BTzqnLseeNRfB6+drve/9gNfCKnGsJWB33LO3QDcAnzc/71o9XYpAHc5524G3gLcY2a3AH8EfM45dx0wDXzEP/4jwLS//XP+cUn2CeDoovX4tYtzrum+gFuBby9avw+4L+q6ImiHUeD5Resv4T0xHmA78JL/+ovAL610XJK/gIeAu9UuS9qkDzgM/DTeHXcd/vaF/6aAbwO3+q87/OMs6tob1B678P7xvgt4GG9K59i1S1P2qIGdwOlF62f8ba1um3PunP/6PLDNf91y7eX/WfpW4ABql/k/758GxoFHgFeAGedc2T9k8c++0C7+/llgU6gFh+ePgd8Gqv76JmLYLs0a1LIG5/2z35LXXvoPYn4A+E3nXGrxvlZtF+dcxTn3Frwe5DuAN0VbUfTM7P3AuHPuUNS1rKVZg/oscPWi9V3+tlZ3wcy2A/jLcX97y7TXKg9ibvl2meecmwG+j/cn/QYzm3/K0+KffaFd/P3DwGS4lYbiNuADZnYC+Bre8MfniWG7NGtQPwlc75+d7QJ+Efi7iGuKg78DPuy//jDeGO389l/xr3K4BZhdNBSQGJd5EHOrt8sWM9vgv+7FG7c/ihfY9/qHLW+X+fa6F/ie/5dIojjn7nPO7XLOjeJlyPecc/+KOLZL1IP5V3AS4L3AT/DG2n4n6noi+Pn/BjgHlPDG0T6CN172KPAy8F1gxD/W8K6SeQV4DhiLuv4GtcnteMMazwJP+1/vVbtwE/CU3y7PA7/rb98LPAEcA/4X0O1v7/HXj/n790b9M4TQRu8EHo5ru+gWchGRmGvWoQ8RkZahoBYRiTkFtYhIzCmoRURiTkEtIhJzCmoRkZhTUIuIxNz/B8Kpm7j5pBYiAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# some materials\n", "si = SLD(2.07, 'Si')\n", "poly1 = SLD(2.0, 'poly1')\n", "poly2 = SLD(4.0, 'poly2')\n", "d2o = SLD(6.36, 'd2o')\n", "\n", "# increase from one to the other\n", "poly1_l = poly1(100., 3.)\n", "poly2_l = poly2(100., 3.)\n", "\n", "ff = FunctionalForm(200., 30., 5., poly1_l, poly2_l)\n", "\n", "struct = si | poly1_l| poly2_l | d2o(0, 3)\n", "struct2 = si | poly1_l | ff | poly2_l | d2o(0, 3)\n", "\n", "struct2.contract = 1.\n", "\n", "plt.plot(*struct.sld_profile())\n", "plt.plot(*struct2.sld_profile());" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAhfElEQVR4nO3deZRkZ3nf8e9bS9fS60xPz2ik2ZCRBJIQWiYCjCCgsEks8YmVww7BJLJPgAMH+ziAsxgCOUlsjOAkh0QJkDjG7HAw0rGwWAx2iCTPoJEstCEJzWi0zUxPL9VV1XXr3vvmj3urp6ene/rWTNVdun6fc1pVXbe66ul7Ss88/bzve19jrUVERNIrl3QAIiJyekrUIiIpp0QtIpJyStQiIimnRC0iknKFfrzoli1b7J49e/rx0iIiG9L+/fuPWWunVjvWl0S9Z88e9u3b14+XFhHZkIwxB9c6ptaHiEjKKVGLiKScErWISMopUYuIpJwStYhIyilRi4iknBK1iEjK9WUetUiWtFyPHz94lMen6ziuj+dbOpf/veicMV5/2faEI5RUsRYO3QFP3wOt+eB76wMWilW45kM9f0slahlox+sOb735Dh56tnbKMWPgTS88V4laTrAWvvPbcO/XVj8+co4StUiv3fSDh3ns2AKff/uVvOzCKSrFPDkDxpikQ5M0+sV3giR9zYfhRb8D1UkwueBf9T5+ZpSoZWAttFy+se8w//jy87juBaqaJYJ9X4SJXXDtv4FcfEN8GkyUgXXg0CzNtscbX3hu0qFIFjgNOPgzuPSGWJM0KFHLANt/cAZj4IpdE0mHIlnw9D1gPdh5dexvrUQtA+vAEzNcuHWUsXIx6VAkC57cH9yee2Xsb61ELQPr4HSD524bSToMyYrpXwaDh6PbYn9rJWoZSL5vOTzbZMdEJelQJCtmDwUDiQlQopaBdGyhheP67NikRC0RzT6hRC0SpydmmgDs2FRNOBLJBGthTolaJFZPznYStSpqiaB+FNxFGFeiFonNkflFALaOlROORDJh4dngdvScRN5eiVoG0kzDIZ8zjJW1OFciaEwHt9XNibx9pERtjJkwxnzTGPOgMeYBY8xL+h2YSD/NNNpsqhZ1TQ+JZilRTyby9lHLic8Ct1lrbzDGDAEagZFMm6k7TFSHkg5DsqJxPLhNa6I2xowDLwf+GYC11gGc/oYl0l8zDYfNStQSVSdRVzYl8vZRWh/PAY4CXzLG3G2M+Z/GmOE+xyXSV7ONNhNVLR2XiBrTUBqHfDKfmSiJugBcCXzeWnsFUAc+svJJxpgbjTH7jDH7jh492uMwRXrreN1hkypqiaoxndhAIkRL1IeBw9baO8Pvv0mQuE9irb3ZWrvXWrt3amqqlzGK9JS1ltlGm03DStQSUdoTtbX2GeAJY8xF4UP/CLi/r1GJ9FHD8XA8X60Pia41D+XxxN4+6qyPDwBfDmd8PAa8p38hifRXbdEFYFRzqCWqVg3GdyT29pE+qdbaA8De/oYiEo+FVpCoR0pK1BJRawFKo4m9vVYmysCph4l6eEiJWiJq1WBIiVokNkuJWhW1ROH74KiiFomVWh/SlXYdsFBKbjcgJWoZOHWnU1HnE45EMqG1ENyqohaJz0LLA1RRS0StWnCrHrVIfNSjlq50ErUqapH41FsuxkB1SK0PicBRohaJ3ULLZXiooGtRSzRLFbUGE0ViU2+5GkiU6Jx6cDukRC0Sm3rLU39aoms3gtticvulKFHLwKk7rmZ8SHTtYMd6isntWK9ELQOn4XhUimp9SESdinoouf1SlKhl4LTaHmUlaomq3YRcIbHdXUCJWgZQs62KWrrQbibanwYlahlAi22fiuZQS1ROPdH+NChRywBqtj3KRX30JaJ2U4laJG6L6lFLN9oNtT5E4qZELV1RRS0SL9fzaXtWg4kSnQYTReK16PoAStQSXbuhilokTk0nuBa1BhMlMrU+ROK12O4kalXUEpFaHyLxUqKWrmnWh0i8FtvqUUuX1PoQiVczrKi1MlEisRbcJhTKiYahRC0D5UTrQx99icB3wfpQVKIWiU1TPWrphrsY3KqiFolPp6JWj1oicVvBrRK1SHxa4WBiSYlaoujs7lIoJRqGErUMlJYbVNSlgj76EoEqapH4tcIl5ENK1BKFetQi8XO8MFHn9dGXCFRRi8TPcZWopQtLFbV61CKxabk+Q/kcuZxJOhTJArU+ROLnuL760xLdUutDFbVIbJSopStuOD1P1/oQiY8Ttj5EIklJRV2I8iRjzONADfAA11q7t59BifSL46mili6kpEcdKVGHXmmtPda3SERi0HI9LXaR6FJSUesTKwNFPWrpSkoq6qifWAv8lTFmvzHmxtWeYIy50Rizzxiz7+jRo72LUKSHWkrU0o2MLXi5xlp7JXAd8D5jzMtXPsFae7O1dq+1du/U1FRPgxTpFQ0mSlfcRciXwCQ77z7SJ9Za+2R4ewT4DnB1P4MS6RdV1NKV9mLi/WmIkKiNMcPGmNHOfeA1wH39DkykHxzXp1TQJU4lIq+VikQdZdbHNuA7Jij9C8CfW2tv62tUIn3ieL5mfUh0bitofSRs3URtrX0MeGEMsYj0nWZ9SFfcFhSGko5C0/NksGgwUbripaOi1idWBkrL9VRRS3Suo4paJG7BYKI+9hKRKmqR+OlaH9IV10nFrA99YmVg+L6l7VklaonOa0FerQ+R2Cztl6hELVGpohaJV0v7JUq3VFGLxKuzsW2pqJWJEpGbjpWJStQyMDqtj5IqaonKc1RRi8SpU1GrRy2RqaIWiZcStXTNczSPWiROLdcDNJgoXdC1PkTidWIwUR97icD3wW+rohaJk6PpedINzwluVVGLxKelBS/SDS/cL1EVtUh8NJgoXXE7FbUStUhsOisTtRWXRLJUUav1IRKbpcFEVdQShRsmalXUIvFR60O60hlMVEUtEh9H86ilG6qoReKny5xKV5YqaiVqkdi02upRSxeWKmq1PkRi43g+OQMFtT4kCs2jFomf42q/ROmCq5WJIrFrub4GEiU6VdQi8Qt2INdiF4lIKxNF4tdq+xpIlOi0MlEkfo6nRC1d0Dxqkfg5rqfBRIlOKxNF4qdZH9IVVdQi8dOsD+mKViaKxM9xfW3DJdG5LcgVIJf8Zyb5CERi4niqqKULKdmBHJSoZYCoRy1dSckO5KBELQMkSNRa8CIReS1V1CJx02CidMV1sldRG2Pyxpi7jTG39DMgkX5paTBRuuEuZrKi/iDwQL8CEek3x/VUUUt0npOKOdQAhShPMsbsAF4PfAr4cF8jOkue53H/HbfROHQ3XrOGsT7Geic9x2KCO2bpP8Ft5+7ScbP0nWX54dyK563+c53vzbLjxix7/2XvaU68yPIfwKx4reXHCvkc50+NMD42Adsvg3MuO+m4nExLyKUrbisVqxIhYqIGbgJ+Hxhd6wnGmBuBGwF27dp11oGdiWajzqOfvZ4XtA6ccsy3QQLLGRtzVDE69wq4/tOw46qkI0kda61mfUh3slRRG2PeAByx1u43xrxiredZa28GbgbYu3dvItnw3i9/jBe1DnDX8/4Vz33luxjfvJVcPo/J5U/f47EWa214d7VbS+cXsn74OH7nR8P/sOLnThzv/LS1Nniv8BWWfmzpNe2J12TZa/r+qcetZa7p8Iffu5+jR57h669uMfZ3n4MvXQf/5Ga45DcinLHB4foW32pjW+mC24JiJekogGgV9UuBNxljrgfKwJgx5s+ste/ob2jdcZoLPP/Jb/DzkZdz9Vs+1t0PG7PUYshS42AL8PE3b+baT/+E/+NexPt++wb46tvgm78VtGcuflPSIaaG44b7JWowUaLyWlCZSDoKIMJgorX2o9baHdbaPcBbgB+lLUkDPHjHbYxRx1z57qRDidX5UyNcvWcz37vnKRiehHd8E867KkjWD92WdHip0UnUqqglMjc9rY8N86mtP/gDHFvg+S9+bdKhxO6aC7bw0LM15hfbUBoNkvU5l8LX3gE/+SMIWyeDzPHCRK0FLxKVuwiFctJRAF0mamvtX1tr39CvYM7G5uP38EjxQsrVNcc7N6wrdk1gLRw4NBs8UB6Hd34Hnnc9/PiT8OiPEo0vDZYqag0mSlS61kdvWd9ne/txauMXJh1KIl64cwKAew/Pnniwsgle88ngfu3p2GNKm5YbTNHU9DyJzF1U66OXnn7iMcZokNt2cdKhJGKsXGTraInHpxsnHyiNBbet+fiDSpmWKmrplnrUvTVz8O8BqJ53acKRJGfP5DAHp+snP1gK20CtWvwBpYxaH9I1VdS91TzyKACbdlyUcCTJ2T1ZPbWizuWhOAyLqqiXpudp1odE4fvgt9Wj7iV35glcm2PrubuTDiUxe7YMc7TWot5yTz5QHlPrg+WzPjbER176zUvPfomwQRJ1ofYkR80khWIx6VASs2NTsILqqdnmyQdKo2p9AK12WFFrep5EsbSxbQan56VVpfkMs8WtSYeRqK2jwQfqSK118oGSKmpQRS1dWkrU6bgo04b41G5yj1CvbE86jERtHQv+RDtSWzz5gCpqQIOJ0iVPFXVPWd9n0j+OU9mWdCiJ2joaJur5FRV1eUyDiWgetXSpU1FrMLE36guzlEwbRga79TFSKlAp5ldpfYxC/SgcuhOc+uo/PAA686iVqCUSN/zLVIOJvTF75EkACqODnaiNMWwdK52aqLdcBM3j8MXXwDfek0xwKXDi6nkaTJQIXCe4VaLujYXpYHl0aWKwWx8QtD+OzK/oUf/6B+CD98LwFMwdTiawFFBFLV1RRd1bjdlnABjZfE7CkSRvcrjETMM5+UFjYNNueO6rBnr2R6vtYQwUclm64rgkxlOPuqecuSMAjG85L+FIkjdRLTLbaK9+sDTYg4otN9gv0WhPSYnC1YKXnrILnUStino8TNTWrrITWmeF4oBemzpI1OpPS0RK1L2Vax5n3lYpDqVjvmOSJipDOJ5Ps+2derA8DlhwBnNOdcv11J+W6LQysbcKziw1M3ibBaxmUzVYQr9q+6NzydMBbX+0XF/7JUp0Sz1qrUzsiSFnjnpeiRqCHjWskajLg31t6pbra79EiU4VdW+V3Xka+bGkw0iF8Urwr/9s0zn14KBX1G31qKUL6lH3VsWbp1UcTzqMVOhU1HOrVtThORrYitpT60Oi0zzq3hrxa7hDStSwrPXRPE2P+tbfg89fA//jWnhyf4zRJaszPU8kEi/8q1Q96h7wfUbtAm55U9KRpMJE2Po4ZdELwObz4ar3wPbLYHxHkKR/9TcxR5gcTc+TrriLwWKXlMy7LyQdwNlwG7MUjIXyRNKhpEJlKE+pkFu99ZEvwBtvCu5bC5+YhMW5WONLkuP6usSpROe2UjOQCBmvqOuzRwEwFVXUHaddndhhTNCzHqBErXnU0hW3lZpNAyDjibpROw5AYViJumOiMrT6rI+VyuMDNbCoWR/SFVXUvdOcDxJ1cWRzwpGkx3i1yMx6FTUMYEWtBS/SBa+VmoFEyHiibi3MAFAeUUXdMV4pMr/arI+VBi5Rq/UhXVBF3TvtxiwAlTFV1B1j5SK1RXf9Jw5cotZgonRBPerecetBRT0yPplwJOmhivpU1locTc+TbriLqqh7xTbn8K1hdGwi6VBSY6xSoNZy8fxVLnW6XHl8YJaTO552d5EueU5qViVCxudR05qjRpWxUjHpSFJjrByci9pim4nqaf50q0xAux7MpwYYPRfefxcUK/0PMmbahku65i6man1GphN1rjXHghlmPCWrh9JgrBIk6vmme/pEffnbwWuD78KRB+GhW2H+KZj8tZgijU+rrY1tpUuuKuqeKTg1GrnhpMNIlfEwUc+t16ceOxde+bHg/sPfDxJ1c7a/wSWk5QYbKZR0mVOJyl1MVaLO9Ce36NZYzI8kHUaqjJWDf3vnFyMMKHZ0VnY2j/chouQttT40j1qicltK1L1SchdwCkrUy51ofXSTqMPpjY2Nmagd9ailW14rNTuQQ8YTddVfoF3UpgHLLSXqbirqapioN3pFrel5ElXWFrwYY8rGmLuMMfcYY35hjPl4HIFFMWzreENK1MtF7lEvVx4HzIatqFvhZr+qqCWylC14iTKY2AKutdYuGGOKwN8aY/7SWntHn2M7Leu5jNDAlrRpwHLDQ3lyJpj1EVkuHyTr5kz/AktQp6LWykSJxPfBb6eq9bFuorbWWmAh/LYYfq2zmqL/FuvzVODEFlMCgDGGsUqxu9YHBO2Ph2+D+hGoboHr/nNwDesNQK0P6YrbDG5TtKYgUolhjMkbYw4AR4DbrbV39jWqCBbmjgGQrypRrzRWjriMfLlLfxOGhuHwPtj3BTj2cH+CS0Bnep4qaomk3UnU1WTjWCbSJ9da61lrLwd2AFcbYy5d+RxjzI3GmH3GmH1Hjx7tcZinaoSXOC1UJ/r+XlkzXil216MGuPZfw/vuhN/4fPB9Y7r3gSWk6QSJujqkiloiaDeC26xV1B3W2lngx8DrVjl2s7V2r7V279TUVI/CW9timKiHdInTU4xVCsxHuYLeaqrhkvINlKgXw8HEslYmShTtcAfyLCVqY8yUMWYivF8BXg082Oe41rW4ECTqsjYNOMUZtT46hrcEt41jvQsoYc0wUVdUUUsUSxV1elofUUaLtgP/2xiTJ0jsX7fW3tLfsNbnNoIZClVdi/oUY+UzGEzs2ICLX5pOMJhYUUUtUSz1qNMzjzrKrI97gStiiKUrfrhpwPD4lmQDSaHx6hn0qDsKQ1Aa21Ctj2bbY6iQI5/TxbskghRW1JkdBvcbs8G1qCdUUa80Vi6w2PaXZjt0rToJ9aPgeye+MqzpuKqmJbp2+qbnZXaibG5xhhpVxguZ/RX6prOMvLboUho5gwQ1PAX3fSv4AjB5eNvX4YJX9TDK+DTbnhK1RJfC6XmZzXJ5Z46aGUGzqE/V2Txgrtlmy8gZrK567afgsb8O7vse/OQ/wjP3ZjhR+5qaJ9GlcHpeZhN1sT1PIz+adBipNH4mV9BbbufVwVfHzz6X6Z510/E0NU+iczvT89JTUWe2R11qz9PM64JMqxmrdK5JfYZzqVeqbgl61hnVbLuamifRpbCizmyirno1nKIq6tV0Wh9nXFGvNDwJ9ezOq2466lFLFzo96ixd5jStqn4Nd0gd6tWc0aVOT6e6JdMLYJptXxW1RNduQKECKdqLNZuJ2lpGbR1flzhd1RltHnA6w1ugnt0e9aJmfUg32s1UtT0go4OJi415ysbDVDWHejWlQo6hfK67a1KfTnUSFp6F733oxGPFKrziI1BO/ziBWh/SFSXq3qjNHqMM5HTlvFUF16Qu9K6i3nNNMKf6wVuD73032LZr90vg+W/szXv0UcPRYKJ0oVWDoXTtxZrJRF2fOcYUUBxWRb2WsTO51OlaLrou+OqYfxr+5HmwcKQ3r99ni+pRSzecBSila6JCJnvUjflgYKs0qkS9lrO6gt56OlfYy8CUvbbn43g+VbU+JKrWApTSVVFnMlEv1oIru1XGdEGmtQTbcfWoR71SvhhcZS8DFfVCeA5Gy5n841GS4CykrvWRyUTdDq9FPTKhRL2WsXKhfxU1wMjWYH/FlFtoBYl6JJxbLrKuVk2tj17ww2slj25Sol7LeKWPrQ8ILty0kP7WRy2sqEdKqqglIg0m9oZtzuLaHJXhiaRDSa3OTuTWWkw/Ju6PbA1mgfy3l538eH4I3vhZOOeUbTUT0amo1fqQSKxN5WBiJj+9ucUZFswwE7lM/kEQi7FykbZn+zfj4Yp3BPNNrT3xmPXgl38Fv/pJahJ1LZyiqIpaInFbwfTTlA0mZvLTO9Q6zlxugomkA0mxzoWZ5prt/iTqX7s2+FrOWvjUdph/qvfvd4ZO9Kgz+VGXuDkLwe1QuirqTJakVWeaWkG7j5/OeK+XkUdhDIyeA7Vn4nvPddQ060O60aoFtymrqDOZqEfcGRaHJpMOI9V6fgW9qEa3Q+3peN/zNJZ61CXN+pAIOok6ZYOJmUzU4/4sbkUzPk6n5xdmivzGKUvUiy75nKFczORHXeLWaX1oMPHsNBt1Rk0ThrcmHUqqjZVP9KhjNbodjn8b/uiCkx/P5eENnzl5KXoMaottRkqF/sx8kY0nnPpLJV2t1cwl6uPPHuY8ID+mRH06m6pDAMzUY07UV74rnA3in/z4PV8J9mGMO1G3XM34kOiaYaJO2ZU5M/cJnj8WJOrypu1Jh5JqE9UixbzhSK0V7xtPXQRv+JNTH3/iLpg9FG8swFyjvTSwKrKuzt6glXQl6sw17prHgv/Zh6f2JBtIyhljmBopcaS2mHQogYldiSTqY3WHyZGh2N9XMqpxPFi0NTScdCQnyVyido8fBGBi+/kJR5J+U2NljsZdUa9lYhfMHATPBd9b8eWv//NnaHqhxZaRUt9eXzaY5vGgmk7ZmEbmWh929glqtsKmzZr1sZ6toyWeON5IOozApt3g1ODfrzKtcmgE3ncnjO/o+dtOLzhMDquilogaM6nrT0MGE/XQwpMcK2xjNGX/4qXR1tES+w/OJB1G4LK3nFieu1z9KNx1Mzx1d88TdcNxabY9JlVRS1SdijplMpeox1pPUStrIDGKraNljtcdHNdnqJBwl2t4El724VMfX5wPEvX0Iz1/y+kFB0AVtUTXmA4GxFMmUz1qr+2w0ztMY/yC9Z8snDtRBuDpuWbCkZxGeQxGtsGxPiTqepioNZgoUVgbXKdm7LykIzlFpirqpx67j53Gg60XJx1KJuzZEoxcPz7dYPdkukaxT7LlQjjwZ8Fc67Xs/S14/R939bLTC8FA6mZV1BLF4mywMnF8Z9KRnCJbifqh/ewEzrnwqqRDyYTdm6sAHJyuA1PJBnM6r/pDeOgv1z7+y+/Dw7d1nagPhQOpOzZVzyI4GRhzh4PbPgxqn61MJWoO/YyGLbHzwsuTjiQTpkZLVIp5Hj+Wkpkfa9mxN/haS2kUfvDvoDnT1dLex4/VGSkV2KLWh0SxlKhVUZ8x6/ucN/0zHqpczhVFjeJHYYxh92SVx44tJB3K2elsQvDVdwQ97fVc+pvwght4fLrBni1VXedDollK1OpRn7H777iNS+wzPHvh7yQdSqZcet44P3rwSP+25IrDzhfBc14eVNStudM/d+4wTD8aJuo6l543Hk+Mkn3P3gfliWBwO2VSlah9z8fz2niug9du4/tt/LbDU4ceZfT23+MYEzz/tf886TAz5cpdm/jm/sMcnG4sDS5mTmkU3v29aM/9m0/DDz9B84m7sccf54rn7YGZx9f/uaHRYAqhDK6n7oZzL0/dqkRIWaJufWIbFeOw8hI640CNCode+yUuGU7XdWLT7qrdQU/3/z02nd1E3Y3d1wBQ+cIr+GkJ2B9+rSdXhA/sD1ZQyuBpN+HZ++HX3590JKtaN1EbY3YCfwpsAyxws7X2s/0I5sD5/wJ8H5srQq4A+QLkCpRHN3P+i97IJVvO6cfbbmgXbhvh/KlhvrX/MG+9elfS4fTfzqvhLX/Ot3/2C+781XE+8aZLKBXX2TOy3YBbfxf+701wwWt7F0u+AM/5h5DX1ftS78FbwW/D+a9MOpJVRamoXeB3rbU/N8aMAvuNMbdba+/vdTAvefd/6PVLDjxjDG+7ehefvPUBvnLXId68dye5XPr+tOsZY/jF2DX8wa/yvOaSbZT+wRXRfu7er8O+LwZfvfSCfwov/pdLsYV3Tr7f82MsO7bK85YfMwYKleBqcYVSKv/s77vGcfjpHwcXDtvzsqSjWZWx1nb3A8Z8F/gv1trb13rO3r177b59+842NumRtufzzi/cyR2PHWe0XGC8UmSokCOX4f8prbX4Fnxr8XyLteD5Fte3HAuvmPfd97+U8yYq0V6wOQvHH+ttkD//U9j/pd6+Zj+ZXHCBrGIVhqrBX7XdvcAZvGmYf6w99f5Sblp5nzUeP8P7rVqw0cXbvw7PfdUZ/A69YYzZb61ddZ5qV4naGLMH+ClwqbV2fsWxG4EbAXbt2nXVwYMHzzhg6b2253PLvU/x84Oz1FsuLc9f9oHPFkswgyVnDHkDuVznviGXg52bq7x5787kL8bke/CrnwYXo1qZJGDtBBL5GGf/mtYHdzFYkec0wKlDux7crtyl53S6LPhOxGXW/ovhlPsr/ho45T4RnrPK/dIovOAG2HbJGfwOvdOTRG2MGQF+AnzKWvvt0z1XFbWISHdOl6gjXZTJGFMEvgV8eb0kLSIivbVuojbBKokvAA9Ya1fZDE9ERPopSkX9UuCdwLXGmAPh1/V9jktERELrDutaa/+WMxvOFRGRHsjUxgEiIoNIiVpEJOWUqEVEUk6JWkQk5bpeQh7pRY05CsS1NHELcCym98oSnZfV6bysTudldXGel93W2lX3zOtLoo6TMWbfWqt5BpnOy+p0Xlan87K6tJwXtT5ERFJOiVpEJOU2QqK+OekAUkrnZXU6L6vTeVldKs5L5nvUIiIb3UaoqEVENjQlahGRlMtsojbGvM4Y85Ax5hFjzEeSjiduxpgvGmOOGGPuW/bYZmPM7caYX4a3m8LHjTHmc+G5utcYc2VykfePMWanMebHxpj7jTG/MMZ8MHx80M9L2RhzlzHmnvC8fDx8/DnGmDvD3/9rxpih8PFS+P0j4fE9if4CfWaMyRtj7jbG3BJ+n7rzkslEbYzJA/8VuA64GHirMebiZKOK3f8CXrfisY8AP7TWXgD8MPwegvN0Qfh1I/D5mGKMW2cj5ouBFwPvCz8Xg35eWsC11toXApcDrzPGvBj4T8BnrLXPBWaA94bPfy8wEz7+mfB5G9kHgQeWfZ++82KtzdwX8BLg+8u+/yjw0aTjSuA87AHuW/b9Q8D28P524KHw/n8H3rra8zbyF/Bd4NU6Lyedkyrwc+BFBCvuCuHjS/9PAd8HXhLeL4TPM0nH3qfzsYPgH+9rgVsILumcuvOSyYoaOA94Ytn3h8PHBt02a+3T4f1ngG3h/YE7X+GfpVcAd6Lz0vnz/gBwBLgdeBSYtda64VOW/+5L5yU8PgdMxhpwfG4Cfh/o7OQ7SQrPS1YTtazDBv/sD+Tcy3Aj5m8BH7LWzi8/NqjnxVrrWWsvJ6ggrwael2xEyTPGvAE4Yq3dn3Qs68lqon4S2Lns+x3hY4PuWWPMdoDw9kj4+MCcrzU2Yh7489JhrZ0FfkzwJ/2EMaazy9Py333pvITHx4HpeCONxUuBNxljHge+StD++CwpPC9ZTdR/B1wQjs4OAW8B/iLhmNLgL4B3h/ffTdCj7Tz+rnCWw4uBuWWtgA3jNBsxD/p5mTLGTIT3KwR9+wcIEvYN4dNWnpfO+boB+FH4l8iGYq39qLV2h7V2D0EO+ZG19u2k8bwk3cw/i0GA64GHCXptf5B0PAn8/l8BngbaBH209xL0y34I/BL4AbA5fK4hmCXzKPD3wN6k4+/TObmGoK1xL3Ag/Lpe54XLgLvD83If8G/Dx88H7gIeAb4BlMLHy+H3j4THz0/6d4jhHL0CuCWt50VLyEVEUi6rrQ8RkYGhRC0iknJK1CIiKadELSKSckrUIiIpp0QtIpJyStQiIin3/wEArd1e675O/QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAhfElEQVR4nO3deZRkZ3nf8e9bS9fS60xPz2ik2ZCRBJIQWiYCjCCgsEks8YmVww7BJLJPgAMH+ziAsxgCOUlsjOAkh0QJkDjG7HAw0rGwWAx2iCTPoJEstCEJzWi0zUxPL9VV1XXr3vvmj3urp6ene/rWTNVdun6fc1pVXbe66ul7Ss88/bzve19jrUVERNIrl3QAIiJyekrUIiIpp0QtIpJyStQiIimnRC0iknKFfrzoli1b7J49e/rx0iIiG9L+/fuPWWunVjvWl0S9Z88e9u3b14+XFhHZkIwxB9c6ptaHiEjKKVGLiKScErWISMopUYuIpJwStYhIyilRi4iknBK1iEjK9WUetUiWtFyPHz94lMen6ziuj+dbOpf/veicMV5/2faEI5RUsRYO3QFP3wOt+eB76wMWilW45kM9f0slahlox+sOb735Dh56tnbKMWPgTS88V4laTrAWvvPbcO/XVj8+co4StUiv3fSDh3ns2AKff/uVvOzCKSrFPDkDxpikQ5M0+sV3giR9zYfhRb8D1UkwueBf9T5+ZpSoZWAttFy+se8w//jy87juBaqaJYJ9X4SJXXDtv4FcfEN8GkyUgXXg0CzNtscbX3hu0qFIFjgNOPgzuPSGWJM0KFHLANt/cAZj4IpdE0mHIlnw9D1gPdh5dexvrUQtA+vAEzNcuHWUsXIx6VAkC57cH9yee2Xsb61ELQPr4HSD524bSToMyYrpXwaDh6PbYn9rJWoZSL5vOTzbZMdEJelQJCtmDwUDiQlQopaBdGyhheP67NikRC0RzT6hRC0SpydmmgDs2FRNOBLJBGthTolaJFZPznYStSpqiaB+FNxFGFeiFonNkflFALaOlROORDJh4dngdvScRN5eiVoG0kzDIZ8zjJW1OFciaEwHt9XNibx9pERtjJkwxnzTGPOgMeYBY8xL+h2YSD/NNNpsqhZ1TQ+JZilRTyby9lHLic8Ct1lrbzDGDAEagZFMm6k7TFSHkg5DsqJxPLhNa6I2xowDLwf+GYC11gGc/oYl0l8zDYfNStQSVSdRVzYl8vZRWh/PAY4CXzLG3G2M+Z/GmOE+xyXSV7ONNhNVLR2XiBrTUBqHfDKfmSiJugBcCXzeWnsFUAc+svJJxpgbjTH7jDH7jh492uMwRXrreN1hkypqiaoxndhAIkRL1IeBw9baO8Pvv0mQuE9irb3ZWrvXWrt3amqqlzGK9JS1ltlGm03DStQSUdoTtbX2GeAJY8xF4UP/CLi/r1GJ9FHD8XA8X60Pia41D+XxxN4+6qyPDwBfDmd8PAa8p38hifRXbdEFYFRzqCWqVg3GdyT29pE+qdbaA8De/oYiEo+FVpCoR0pK1BJRawFKo4m9vVYmysCph4l6eEiJWiJq1WBIiVokNkuJWhW1ROH74KiiFomVWh/SlXYdsFBKbjcgJWoZOHWnU1HnE45EMqG1ENyqohaJz0LLA1RRS0StWnCrHrVIfNSjlq50ErUqapH41FsuxkB1SK0PicBRohaJ3ULLZXiooGtRSzRLFbUGE0ViU2+5GkiU6Jx6cDukRC0Sm3rLU39aoms3gtticvulKFHLwKk7rmZ8SHTtYMd6isntWK9ELQOn4XhUimp9SESdinoouf1SlKhl4LTaHmUlaomq3YRcIbHdXUCJWgZQs62KWrrQbibanwYlahlAi22fiuZQS1ROPdH+NChRywBqtj3KRX30JaJ2U4laJG6L6lFLN9oNtT5E4qZELV1RRS0SL9fzaXtWg4kSnQYTReK16PoAStQSXbuhilokTk0nuBa1BhMlMrU+ROK12O4kalXUEpFaHyLxUqKWrmnWh0i8FtvqUUuX1PoQiVczrKi1MlEisRbcJhTKiYahRC0D5UTrQx99icB3wfpQVKIWiU1TPWrphrsY3KqiFolPp6JWj1oicVvBrRK1SHxa4WBiSYlaoujs7lIoJRqGErUMlJYbVNSlgj76EoEqapH4tcIl5ENK1BKFetQi8XO8MFHn9dGXCFRRi8TPcZWopQtLFbV61CKxabk+Q/kcuZxJOhTJArU+ROLnuL760xLdUutDFbVIbJSopStuOD1P1/oQiY8Ttj5EIklJRV2I8iRjzONADfAA11q7t59BifSL46mili6kpEcdKVGHXmmtPda3SERi0HI9LXaR6FJSUesTKwNFPWrpSkoq6qifWAv8lTFmvzHmxtWeYIy50Rizzxiz7+jRo72LUKSHWkrU0o2MLXi5xlp7JXAd8D5jzMtXPsFae7O1dq+1du/U1FRPgxTpFQ0mSlfcRciXwCQ77z7SJ9Za+2R4ewT4DnB1P4MS6RdV1NKV9mLi/WmIkKiNMcPGmNHOfeA1wH39DkykHxzXp1TQJU4lIq+VikQdZdbHNuA7Jij9C8CfW2tv62tUIn3ieL5mfUh0bitofSRs3URtrX0MeGEMsYj0nWZ9SFfcFhSGko5C0/NksGgwUbripaOi1idWBkrL9VRRS3Suo4paJG7BYKI+9hKRKmqR+OlaH9IV10nFrA99YmVg+L6l7VklaonOa0FerQ+R2Cztl6hELVGpohaJV0v7JUq3VFGLxKuzsW2pqJWJEpGbjpWJStQyMDqtj5IqaonKc1RRi8SpU1GrRy2RqaIWiZcStXTNczSPWiROLdcDNJgoXdC1PkTidWIwUR97icD3wW+rohaJk6PpedINzwluVVGLxKelBS/SDS/cL1EVtUh8NJgoXXE7FbUStUhsOisTtRWXRLJUUav1IRKbpcFEVdQShRsmalXUIvFR60O60hlMVEUtEh9H86ilG6qoReKny5xKV5YqaiVqkdi02upRSxeWKmq1PkRi43g+OQMFtT4kCs2jFomf42q/ROmCq5WJIrFrub4GEiU6VdQi8Qt2INdiF4lIKxNF4tdq+xpIlOi0MlEkfo6nRC1d0Dxqkfg5rqfBRIlOKxNF4qdZH9IVVdQi8dOsD+mKViaKxM9xfW3DJdG5LcgVIJf8Zyb5CERi4niqqKULKdmBHJSoZYCoRy1dSckO5KBELQMkSNRa8CIReS1V1CJx02CidMV1sldRG2Pyxpi7jTG39DMgkX5paTBRuuEuZrKi/iDwQL8CEek3x/VUUUt0npOKOdQAhShPMsbsAF4PfAr4cF8jOkue53H/HbfROHQ3XrOGsT7Geic9x2KCO2bpP8Ft5+7ScbP0nWX54dyK563+c53vzbLjxix7/2XvaU68yPIfwKx4reXHCvkc50+NMD42Adsvg3MuO+m4nExLyKUrbisVqxIhYqIGbgJ+Hxhd6wnGmBuBGwF27dp11oGdiWajzqOfvZ4XtA6ccsy3QQLLGRtzVDE69wq4/tOw46qkI0kda61mfUh3slRRG2PeAByx1u43xrxiredZa28GbgbYu3dvItnw3i9/jBe1DnDX8/4Vz33luxjfvJVcPo/J5U/f47EWa214d7VbS+cXsn74OH7nR8P/sOLnThzv/LS1Nniv8BWWfmzpNe2J12TZa/r+qcetZa7p8Iffu5+jR57h669uMfZ3n4MvXQf/5Ga45DcinLHB4foW32pjW+mC24JiJekogGgV9UuBNxljrgfKwJgx5s+ste/ob2jdcZoLPP/Jb/DzkZdz9Vs+1t0PG7PUYshS42AL8PE3b+baT/+E/+NexPt++wb46tvgm78VtGcuflPSIaaG44b7JWowUaLyWlCZSDoKIMJgorX2o9baHdbaPcBbgB+lLUkDPHjHbYxRx1z57qRDidX5UyNcvWcz37vnKRiehHd8E867KkjWD92WdHip0UnUqqglMjc9rY8N86mtP/gDHFvg+S9+bdKhxO6aC7bw0LM15hfbUBoNkvU5l8LX3gE/+SMIWyeDzPHCRK0FLxKVuwiFctJRAF0mamvtX1tr39CvYM7G5uP38EjxQsrVNcc7N6wrdk1gLRw4NBs8UB6Hd34Hnnc9/PiT8OiPEo0vDZYqag0mSlS61kdvWd9ne/txauMXJh1KIl64cwKAew/Pnniwsgle88ngfu3p2GNKm5YbTNHU9DyJzF1U66OXnn7iMcZokNt2cdKhJGKsXGTraInHpxsnHyiNBbet+fiDSpmWKmrplnrUvTVz8O8BqJ53acKRJGfP5DAHp+snP1gK20CtWvwBpYxaH9I1VdS91TzyKACbdlyUcCTJ2T1ZPbWizuWhOAyLqqiXpudp1odE4fvgt9Wj7iV35glcm2PrubuTDiUxe7YMc7TWot5yTz5QHlPrg+WzPjbER176zUvPfomwQRJ1ofYkR80khWIx6VASs2NTsILqqdnmyQdKo2p9AK12WFFrep5EsbSxbQan56VVpfkMs8WtSYeRqK2jwQfqSK118oGSKmpQRS1dWkrU6bgo04b41G5yj1CvbE86jERtHQv+RDtSWzz5gCpqQIOJ0iVPFXVPWd9n0j+OU9mWdCiJ2joaJur5FRV1eUyDiWgetXSpU1FrMLE36guzlEwbRga79TFSKlAp5ldpfYxC/SgcuhOc+uo/PAA686iVqCUSN/zLVIOJvTF75EkACqODnaiNMWwdK52aqLdcBM3j8MXXwDfek0xwKXDi6nkaTJQIXCe4VaLujYXpYHl0aWKwWx8QtD+OzK/oUf/6B+CD98LwFMwdTiawFFBFLV1RRd1bjdlnABjZfE7CkSRvcrjETMM5+UFjYNNueO6rBnr2R6vtYQwUclm64rgkxlOPuqecuSMAjG85L+FIkjdRLTLbaK9+sDTYg4otN9gv0WhPSYnC1YKXnrILnUStino8TNTWrrITWmeF4oBemzpI1OpPS0RK1L2Vax5n3lYpDqVjvmOSJipDOJ5Ps+2derA8DlhwBnNOdcv11J+W6LQysbcKziw1M3ibBaxmUzVYQr9q+6NzydMBbX+0XF/7JUp0Sz1qrUzsiSFnjnpeiRqCHjWskajLg31t6pbra79EiU4VdW+V3Xka+bGkw0iF8Urwr/9s0zn14KBX1G31qKUL6lH3VsWbp1UcTzqMVOhU1HOrVtThORrYitpT60Oi0zzq3hrxa7hDStSwrPXRPE2P+tbfg89fA//jWnhyf4zRJaszPU8kEi/8q1Q96h7wfUbtAm55U9KRpMJE2Po4ZdELwObz4ar3wPbLYHxHkKR/9TcxR5gcTc+TrriLwWKXlMy7LyQdwNlwG7MUjIXyRNKhpEJlKE+pkFu99ZEvwBtvCu5bC5+YhMW5WONLkuP6usSpROe2UjOQCBmvqOuzRwEwFVXUHaddndhhTNCzHqBErXnU0hW3lZpNAyDjibpROw5AYViJumOiMrT6rI+VyuMDNbCoWR/SFVXUvdOcDxJ1cWRzwpGkx3i1yMx6FTUMYEWtBS/SBa+VmoFEyHiibi3MAFAeUUXdMV4pMr/arI+VBi5Rq/UhXVBF3TvtxiwAlTFV1B1j5SK1RXf9Jw5cotZgonRBPerecetBRT0yPplwJOmhivpU1locTc+TbriLqqh7xTbn8K1hdGwi6VBSY6xSoNZy8fxVLnW6XHl8YJaTO552d5EueU5qViVCxudR05qjRpWxUjHpSFJjrByci9pim4nqaf50q0xAux7MpwYYPRfefxcUK/0PMmbahku65i6man1GphN1rjXHghlmPCWrh9JgrBIk6vmme/pEffnbwWuD78KRB+GhW2H+KZj8tZgijU+rrY1tpUuuKuqeKTg1GrnhpMNIlfEwUc+t16ceOxde+bHg/sPfDxJ1c7a/wSWk5QYbKZR0mVOJyl1MVaLO9Ce36NZYzI8kHUaqjJWDf3vnFyMMKHZ0VnY2j/chouQttT40j1qicltK1L1SchdwCkrUy51ofXSTqMPpjY2Nmagd9ailW14rNTuQQ8YTddVfoF3UpgHLLSXqbirqapioN3pFrel5ElXWFrwYY8rGmLuMMfcYY35hjPl4HIFFMWzreENK1MtF7lEvVx4HzIatqFvhZr+qqCWylC14iTKY2AKutdYuGGOKwN8aY/7SWntHn2M7Leu5jNDAlrRpwHLDQ3lyJpj1EVkuHyTr5kz/AktQp6LWykSJxPfBb6eq9bFuorbWWmAh/LYYfq2zmqL/FuvzVODEFlMCgDGGsUqxu9YHBO2Ph2+D+hGoboHr/nNwDesNQK0P6YrbDG5TtKYgUolhjMkbYw4AR4DbrbV39jWqCBbmjgGQrypRrzRWjriMfLlLfxOGhuHwPtj3BTj2cH+CS0Bnep4qaomk3UnU1WTjWCbSJ9da61lrLwd2AFcbYy5d+RxjzI3GmH3GmH1Hjx7tcZinaoSXOC1UJ/r+XlkzXil216MGuPZfw/vuhN/4fPB9Y7r3gSWk6QSJujqkiloiaDeC26xV1B3W2lngx8DrVjl2s7V2r7V279TUVI/CW9timKiHdInTU4xVCsxHuYLeaqrhkvINlKgXw8HEslYmShTtcAfyLCVqY8yUMWYivF8BXg082Oe41rW4ECTqsjYNOMUZtT46hrcEt41jvQsoYc0wUVdUUUsUSxV1elofUUaLtgP/2xiTJ0jsX7fW3tLfsNbnNoIZClVdi/oUY+UzGEzs2ICLX5pOMJhYUUUtUSz1qNMzjzrKrI97gStiiKUrfrhpwPD4lmQDSaHx6hn0qDsKQ1Aa21Ctj2bbY6iQI5/TxbskghRW1JkdBvcbs8G1qCdUUa80Vi6w2PaXZjt0rToJ9aPgeye+MqzpuKqmJbp2+qbnZXaibG5xhhpVxguZ/RX6prOMvLboUho5gwQ1PAX3fSv4AjB5eNvX4YJX9TDK+DTbnhK1RJfC6XmZzXJ5Z46aGUGzqE/V2Txgrtlmy8gZrK567afgsb8O7vse/OQ/wjP3ZjhR+5qaJ9GlcHpeZhN1sT1PIz+adBipNH4mV9BbbufVwVfHzz6X6Z510/E0NU+iczvT89JTUWe2R11qz9PM64JMqxmrdK5JfYZzqVeqbgl61hnVbLuamifRpbCizmyirno1nKIq6tV0Wh9nXFGvNDwJ9ezOq2466lFLFzo96ixd5jStqn4Nd0gd6tWc0aVOT6e6JdMLYJptXxW1RNduQKECKdqLNZuJ2lpGbR1flzhd1RltHnA6w1ugnt0e9aJmfUg32s1UtT0go4OJi415ysbDVDWHejWlQo6hfK67a1KfTnUSFp6F733oxGPFKrziI1BO/ziBWh/SFSXq3qjNHqMM5HTlvFUF16Qu9K6i3nNNMKf6wVuD73032LZr90vg+W/szXv0UcPRYKJ0oVWDoXTtxZrJRF2fOcYUUBxWRb2WsTO51OlaLrou+OqYfxr+5HmwcKQ3r99ni+pRSzecBSila6JCJnvUjflgYKs0qkS9lrO6gt56OlfYy8CUvbbn43g+VbU+JKrWApTSVVFnMlEv1oIru1XGdEGmtQTbcfWoR71SvhhcZS8DFfVCeA5Gy5n841GS4CykrvWRyUTdDq9FPTKhRL2WsXKhfxU1wMjWYH/FlFtoBYl6JJxbLrKuVk2tj17ww2slj25Sol7LeKWPrQ8ILty0kP7WRy2sqEdKqqglIg0m9oZtzuLaHJXhiaRDSa3OTuTWWkw/Ju6PbA1mgfy3l538eH4I3vhZOOeUbTUT0amo1fqQSKxN5WBiJj+9ucUZFswwE7lM/kEQi7FykbZn+zfj4Yp3BPNNrT3xmPXgl38Fv/pJahJ1LZyiqIpaInFbwfTTlA0mZvLTO9Q6zlxugomkA0mxzoWZ5prt/iTqX7s2+FrOWvjUdph/qvfvd4ZO9Kgz+VGXuDkLwe1QuirqTJakVWeaWkG7j5/OeK+XkUdhDIyeA7Vn4nvPddQ060O60aoFtymrqDOZqEfcGRaHJpMOI9V6fgW9qEa3Q+3peN/zNJZ61CXN+pAIOok6ZYOJmUzU4/4sbkUzPk6n5xdmivzGKUvUiy75nKFczORHXeLWaX1oMPHsNBt1Rk0ThrcmHUqqjZVP9KhjNbodjn8b/uiCkx/P5eENnzl5KXoMaottRkqF/sx8kY0nnPpLJV2t1cwl6uPPHuY8ID+mRH06m6pDAMzUY07UV74rnA3in/z4PV8J9mGMO1G3XM34kOiaYaJO2ZU5M/cJnj8WJOrypu1Jh5JqE9UixbzhSK0V7xtPXQRv+JNTH3/iLpg9FG8swFyjvTSwKrKuzt6glXQl6sw17prHgv/Zh6f2JBtIyhljmBopcaS2mHQogYldiSTqY3WHyZGh2N9XMqpxPFi0NTScdCQnyVyido8fBGBi+/kJR5J+U2NljsZdUa9lYhfMHATPBd9b8eWv//NnaHqhxZaRUt9eXzaY5vGgmk7ZmEbmWh929glqtsKmzZr1sZ6toyWeON5IOozApt3g1ODfrzKtcmgE3ncnjO/o+dtOLzhMDquilogaM6nrT0MGE/XQwpMcK2xjNGX/4qXR1tES+w/OJB1G4LK3nFieu1z9KNx1Mzx1d88TdcNxabY9JlVRS1SdijplMpeox1pPUStrIDGKraNljtcdHNdnqJBwl2t4El724VMfX5wPEvX0Iz1/y+kFB0AVtUTXmA4GxFMmUz1qr+2w0ztMY/yC9Z8snDtRBuDpuWbCkZxGeQxGtsGxPiTqepioNZgoUVgbXKdm7LykIzlFpirqpx67j53Gg60XJx1KJuzZEoxcPz7dYPdkukaxT7LlQjjwZ8Fc67Xs/S14/R939bLTC8FA6mZV1BLF4mywMnF8Z9KRnCJbifqh/ewEzrnwqqRDyYTdm6sAHJyuA1PJBnM6r/pDeOgv1z7+y+/Dw7d1nagPhQOpOzZVzyI4GRhzh4PbPgxqn61MJWoO/YyGLbHzwsuTjiQTpkZLVIp5Hj+Wkpkfa9mxN/haS2kUfvDvoDnT1dLex4/VGSkV2KLWh0SxlKhVUZ8x6/ucN/0zHqpczhVFjeJHYYxh92SVx44tJB3K2elsQvDVdwQ97fVc+pvwght4fLrBni1VXedDollK1OpRn7H777iNS+wzPHvh7yQdSqZcet44P3rwSP+25IrDzhfBc14eVNStudM/d+4wTD8aJuo6l543Hk+Mkn3P3gfliWBwO2VSlah9z8fz2niug9du4/tt/LbDU4ceZfT23+MYEzz/tf886TAz5cpdm/jm/sMcnG4sDS5mTmkU3v29aM/9m0/DDz9B84m7sccf54rn7YGZx9f/uaHRYAqhDK6n7oZzL0/dqkRIWaJufWIbFeOw8hI640CNCode+yUuGU7XdWLT7qrdQU/3/z02nd1E3Y3d1wBQ+cIr+GkJ2B9+rSdXhA/sD1ZQyuBpN+HZ++HX3590JKtaN1EbY3YCfwpsAyxws7X2s/0I5sD5/wJ8H5srQq4A+QLkCpRHN3P+i97IJVvO6cfbbmgXbhvh/KlhvrX/MG+9elfS4fTfzqvhLX/Ot3/2C+781XE+8aZLKBXX2TOy3YBbfxf+701wwWt7F0u+AM/5h5DX1ftS78FbwW/D+a9MOpJVRamoXeB3rbU/N8aMAvuNMbdba+/vdTAvefd/6PVLDjxjDG+7ehefvPUBvnLXId68dye5XPr+tOsZY/jF2DX8wa/yvOaSbZT+wRXRfu7er8O+LwZfvfSCfwov/pdLsYV3Tr7f82MsO7bK85YfMwYKleBqcYVSKv/s77vGcfjpHwcXDtvzsqSjWZWx1nb3A8Z8F/gv1trb13rO3r177b59+842NumRtufzzi/cyR2PHWe0XGC8UmSokCOX4f8prbX4Fnxr8XyLteD5Fte3HAuvmPfd97+U8yYq0V6wOQvHH+ttkD//U9j/pd6+Zj+ZXHCBrGIVhqrBX7XdvcAZvGmYf6w99f5Sblp5nzUeP8P7rVqw0cXbvw7PfdUZ/A69YYzZb61ddZ5qV4naGLMH+ClwqbV2fsWxG4EbAXbt2nXVwYMHzzhg6b2253PLvU/x84Oz1FsuLc9f9oHPFkswgyVnDHkDuVznviGXg52bq7x5787kL8bke/CrnwYXo1qZJGDtBBL5GGf/mtYHdzFYkec0wKlDux7crtyl53S6LPhOxGXW/ovhlPsr/ho45T4RnrPK/dIovOAG2HbJGfwOvdOTRG2MGQF+AnzKWvvt0z1XFbWISHdOl6gjXZTJGFMEvgV8eb0kLSIivbVuojbBKokvAA9Ya1fZDE9ERPopSkX9UuCdwLXGmAPh1/V9jktERELrDutaa/+WMxvOFRGRHsjUxgEiIoNIiVpEJOWUqEVEUk6JWkQk5bpeQh7pRY05CsS1NHELcCym98oSnZfV6bysTudldXGel93W2lX3zOtLoo6TMWbfWqt5BpnOy+p0Xlan87K6tJwXtT5ERFJOiVpEJOU2QqK+OekAUkrnZXU6L6vTeVldKs5L5nvUIiIb3UaoqEVENjQlahGRlMtsojbGvM4Y85Ax5hFjzEeSjiduxpgvGmOOGGPuW/bYZmPM7caYX4a3m8LHjTHmc+G5utcYc2VykfePMWanMebHxpj7jTG/MMZ8MHx80M9L2RhzlzHmnvC8fDx8/DnGmDvD3/9rxpih8PFS+P0j4fE9if4CfWaMyRtj7jbG3BJ+n7rzkslEbYzJA/8VuA64GHirMebiZKOK3f8CXrfisY8AP7TWXgD8MPwegvN0Qfh1I/D5mGKMW2cj5ouBFwPvCz8Xg35eWsC11toXApcDrzPGvBj4T8BnrLXPBWaA94bPfy8wEz7+mfB5G9kHgQeWfZ++82KtzdwX8BLg+8u+/yjw0aTjSuA87AHuW/b9Q8D28P524KHw/n8H3rra8zbyF/Bd4NU6Lyedkyrwc+BFBCvuCuHjS/9PAd8HXhLeL4TPM0nH3qfzsYPgH+9rgVsILumcuvOSyYoaOA94Ytn3h8PHBt02a+3T4f1ngG3h/YE7X+GfpVcAd6Lz0vnz/gBwBLgdeBSYtda64VOW/+5L5yU8PgdMxhpwfG4Cfh/o7OQ7SQrPS1YTtazDBv/sD+Tcy3Aj5m8BH7LWzi8/NqjnxVrrWWsvJ6ggrwael2xEyTPGvAE4Yq3dn3Qs68lqon4S2Lns+x3hY4PuWWPMdoDw9kj4+MCcrzU2Yh7489JhrZ0FfkzwJ/2EMaazy9Py333pvITHx4HpeCONxUuBNxljHge+StD++CwpPC9ZTdR/B1wQjs4OAW8B/iLhmNLgL4B3h/ffTdCj7Tz+rnCWw4uBuWWtgA3jNBsxD/p5mTLGTIT3KwR9+wcIEvYN4dNWnpfO+boB+FH4l8iGYq39qLV2h7V2D0EO+ZG19u2k8bwk3cw/i0GA64GHCXptf5B0PAn8/l8BngbaBH209xL0y34I/BL4AbA5fK4hmCXzKPD3wN6k4+/TObmGoK1xL3Ag/Lpe54XLgLvD83If8G/Dx88H7gIeAb4BlMLHy+H3j4THz0/6d4jhHL0CuCWt50VLyEVEUi6rrQ8RkYGhRC0iknJK1CIiKadELSKSckrUIiIpp0QtIpJyStQiIin3/wEArd1e675O/QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "si = SLD(2.07, 'Si')\n", "poly1 = SLD(4.0, 'poly1')\n", "poly2 = SLD(2.0, 'poly2')\n", "d2o = SLD(6.36, 'd2o')\n", "\n", "# decay from one to the other\n", "poly1_l = poly1(100., 3.)\n", "poly2_l = poly2(100., 3.)\n", "\n", "ff = FunctionalForm(200., 30., 3., poly1_l, poly2_l)\n", "\n", "struct = si | poly1_l| poly2_l | d2o(0, 3)\n", "struct2 = si | poly1_l | ff | poly2_l | d2o(0, 3)\n", "\n", "struct2.contract = 1.\n", "\n", "plt.plot(*struct.sld_profile())\n", "plt.plot(*struct2.sld_profile());" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "si = SLD(2.07, 'Si')\n", "poly1 = SLD(3.0, 'poly1')\n", "d2o = SLD(6.36, 'd2o')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAa3ElEQVR4nO3de5Cd9X3f8fd379qbdiUtuoKEAJmbEZK3RA6YTiA4GF9oGo+HdOLQ2hNNXdratT0NaWY8cTzJtJ0mxWk9dmTXjp2QOKlbSuoO1CSQOthGIIEE4iKQQAIJXVaXvZ2ze66//vE8u1qt9vKcc57nPM959vOa0TznPM/Zsx8Oq49++3tu5pxDRESSqynuACIisjAVtYhIwqmoRUQSTkUtIpJwKmoRkYRrieJNV61a5TZt2hTFW4uIpNLevXvPOOcG5toWSVFv2rSJPXv2RPHWIiKpZGZH59umqQ8RkYRTUYuIJJyKWkQk4VTUIiIJp6IWEUk4FbWISMKpqEVEEi6S46hlaSuUyjzz5lkOnx4nky+RK5a9DTMuqfueNb18+Ka1MSWUWDg342dg1uOp7YEez/H1cz5mnq+f63EZ8lnIj0M+4/8Zn/F8HEqFOd53lrYuuO1zC38OVVBRS6gOD43zqT95jqNns3NuN/OWH71pnYo6rYZeh2/9IuRGmS7C1LKLn3ZfpqKWZJvIl7j/288yWSjxjV/bzuCmFXS3t9De0oSZLf4GEr2f/hc4/GS03+P8USgX4fYvAnbhX+cFHzP3evCfV/p4oe8312sM2rq9EXFb16zH/vOWtio/kNqpqCU0f7XnHY6dn+AvfmMH779qZdxxZLbCBDz5e9C5Enoj/G2mcwXs+Azc8hvRfY8lRkUtoflf+45z3dpelXRUnn4Ijj1X/ddPjkBxAj72Vbj6F0OLJdFTUUsoToxM8MLbw/zm3dfGHSWdChPw5Fe80XDnqurf56o7YdMHwssldaGillA8f3QYgFuv1mg6Eif2e/O+H3kIrr0n7jRSZypqCcX+Y8O0tTRx7ZreuKM0hqd+H579ZvDXl/LecsM/iCaPJJqKWkKx751hbljXS1uLzqEK5JVHYVk/XHVH8K9ZeTV0z3ldeUk5FbWE4vVTY9zzXh0XHUgxB2cPwa2fhTu/FHcaaQAqaqnZyESB4WyBTSs7444Sr/NH4MhPFn9d5rQ333zZ9ZFHknRQUUvN3vbPQrxiRVfMSWL22IPw+mPBXmtNsH57tHkkNVTUUrOj5zIAbFzqI+qzb8A1H4R7/tPir23rhi4dISPBqKilZkenR9RLuKjLJe/U6Ws/DP0b404jKaNd9FKzY+ezrOxqo6t9Cf+7P3ocygVYsTnuJJJCgf5mmVkf8C3gRrzLYX3KOfezCHNJAzk9mmN1b0fcMaL16ANwcIH551LRW/ZfWZ88sqQEHQJ9FXjcOfdxM2sDlvDvuDLb0HiOgZ72uGNE6+Bj0LMWrtgx/2vaexbeLlKlRYvazJYDtwP/FMA5lwfy0caSRjI0lmPL6p64Y0Qnn4HsWXj/A/CBL8SdRpagIHPUVwJDwHfM7AUz+5aZXXIclpntNLM9ZrZnaGgo9KCSTOWyY2gs5SPqkWPecvkV8eaQJStIUbcA24GvO+e2ARngwdkvcs7tcs4NOucGBwZ0mutSMTxRoFh2XJbmoh5+x1v2XR5vDlmygsxRHwOOOed2+89/wBxFLUvT0FgOoPFH1Lkx+NN/DBPn5t4GsHxDfTOJ+BYtaufcSTN7x8ze45w7CNwJvBJ9NGkE00Xd3eBFffpVOPasd63m7ssu3b58A/Sur38uEYIf9fGvgIf9Iz7eBP5ZdJGkkZzNeEW9stGLevS4t7z738OaG+PNIjJLoKJ2zu0DBqONIo1oZKIAQF9na8xJajT6rrfsXRdvDpE56MxEqclw1ivq5ctSUNStnd41okUSZgmf8ythGJko0NXWTGtzg/ybf/4InDl06foT+73RtFndI4ksRkUtNRnOFhprNP3wJ+DMwbm3bbm7vllEAlJRS01GJgr0NkpROwfDb8N7PwG37Lx0+8CW+mcSCUBFLTUZnSg0zo7E/DgUJ7yjOi7XTWKlcTTIxKIk1fBEvnGmPsZPe8vu1fHmEKmQilpqMjJRoG9ZW9wxgpkq6i5d4kAai4paajIyUWB5o0x9jJ/ylhpRS4PRHLVUbbJQYrJQTubUx9//IZw6cPG6c295SxW1NBgVtVRtdCKhJ7uUy/DkV6CjDzpXXLzt6rugUzeVlcaiopaqDSe1qCeHwZXhH/4m7PjncacRqZnmqKVqib3OR+aMt9TIWVJCRS1VG/Gv89HbkbCizp71lrOnPUQalIpaqpbJe3fe7u5I2AzaVFF3rYo3h0hIVNRStUyuBEB3e9KKWlMfki4qaqla1h9Rd7Y1x5xklqkR9TJNfUg6JGwoJI1kakTd2Rbjj9HBx+EHn4Jy8cK6cgFau6CtM75cIiFSUUvVsvkiHa1NNDfFeA3n43uhkIVbP3vx+jXvjSePSARU1FK1TL5IV5yjafCmOTpXwF1fjjeHSIQ0Ry1Vy+ZKdLbHPD89cU5z0ZJ6KmqpWnJG1Dq6Q9JNRS1Vy+ZL8R/xkT2nopbUU1FL1cZzRbriPoZ6ao5aJMW0M1Gqls2VuKynvb7f1DkoTEw9UVHLkqCilqpl8jGMqH/4b2Dvdy5epzu2SMqpqKVq2Xyp/jsTT74EK6+GbZ/0nje1wNZfrW8GkTpTUUvVMrli/Q/PmzgPa7fCbZ+r7/cViZF2JkpViqUyuWK5/iPqifOak5YlR0UtVckWpq7zUccRdbns3b1lWX/9vqdIAqiopSpZ/4JMdd2ZmBv1brGlopYlRkUtVcnEcYnTifPeUkUtS4yKWqqSyXlFXdc5ahW1LFE66kOqMn0t6iiP+hh+Gw4+5p3kAnD2kLfURZhkiVFRS1Wm7u4S6Yj67/8A9v7JxetaOqB/Y3TfUySBVNRSlUy+DjsTs2dh5TXw6R9dWNfSoTu3yJKjopaqZKfmqKOc+pgc9eajddy0LHHamShVmRpRR3q/xNwodCyP7v1FGkSgv2VmdgQYA0pA0Tk3GGUoSb6pEXWkh+dNjsCKzdG9v0iDqGQ49AvOuTORJZGGksmXaGtporU5wl/KJjWiFgFNfUiVsvkiXVGOpp3zpj7ae6P7HiINImhRO+BHZrbXzHbO9QIz22lme8xsz9DQUHgJJZHGc8Vo56eLk1DKQ4eKWiRoUd/mnNsOfAh4wMxun/0C59wu59ygc25wYEAXck+7bK4U/REfoKkPEQLOUTvnjvvL02b2CHAL8OMog0myZfIRjKhHjl04+3D0hLdsV1GLLPo3zcy6gCbn3Jj/+IPA70aeTBItm49gRP1nvwJDr128rmdNuN9DpAEFGRKtBh4xs6nX/7lz7vFIU0niZXJFVnSFfIbg2Em47mOw4zPe89ZlsPbmcL+HSANatKidc28CW+uQRRpINl+iO8zTx52D3Bis2gIbfz689xVJAR2eJ1XJ5ovhnuxSyIIrQXtPeO8pkhIqaqlKJlcK94JM00d56HA8kdlU1FKxUtkxUSiFO6LO+UWtE1xELqGilopFci3q3Ji3VFGLXEJFLRXL5iO4u8vkiLfU1IfIJVTUUrFI7pc4PfWhnYkis6mopWLTI+pQ56g19SEyH93hRSo2PaKu9aiPoYOw5zvgynDqgLdOUx8il1BRS8VCG1E//z3Y/XXo6POer38ftGnqQ2Q2FbVULOMf9VHzmYm5UeheA188GEIqkfTSHLVULJubOuqj1qIe085DkQBU1FKx8emjPmqc+lBRiwSiopaKTZ3wUvP1qCdHVdQiAaiopWKZfInWZqOtpcYfn9yYjvIQCUBFLRXLhnW/xNyYjpsWCUBFLRXL5Evh3IFcc9QigaiopWLZfLH2Iz7KZe/wPI2oRRalopaKZXIhjKgLGcBpRC0SgE54kYplq70DeTEPT/0eTA5DYcJbp6IWWZSKWiqWyZVY19dW+Reeegl+8pB3ynhLOyy/AtbdHHI6kfRRUUvFMvkiXdVcizo37i3vexg23RZuKJEU0xy1VCyTK1U39ZHPeMu2rnADiaSciloqls0Xq9uZmPdH1LpCnkhFVNRSkXLZkc2Xqjs8b7qoNaIWqYSKWioyUfCunFfViDqnohaphopaKjJ1LerqRtSaoxaphopaKjJ1Leqq56hbO6EpxHstiiwBKmqpSKaWS5zmx6GtO+REIumnopaKZKZG1NUcR53PaNpDpAoqaqlITXcgz41Du0bUIpXSmYlSkanbcAW+sW25DMVJcP7V8jT1IVIxFbVUZOo2XIFH1H/+CTj0xIXnW+6OIJVIuqmopSLj/hx1d9CdiadfgXXb4YZfBjO46s4I04mkk4paKnJhjjrgzsR8Bta/D2791xGmEkk37UyUimRyRdpbmmhpDvijU5iAts5oQ4mknIpaKjKeK1awI7EEpRy06pA8kVqoqKUimVwx+I7E6VPGNaIWqUXgojazZjN7wcx+GGUgSbbxXCl4URey3rJVRS1Si0pG1J8FXo0qiDSGTK5IdyU7EkFnI4rUKFBRm9kG4MPAt6KNI0mXqeTGthpRi4Qi6Ij6IeDfAuX5XmBmO81sj5ntGRoaCiObJFBFOxOn7jSuohapyaJFbWYfAU475/Yu9Drn3C7n3KBzbnBgYCC0gJIs3s7ESqc+VNQitQgyor4V+JiZHQG+D9xhZn8WaSpJrIx2JorU3aJF7Zz7LefcBufcJuA+4Enn3K9FnkwSxzlHJl/B1EfeL2rtTBSpiY6jlsCy+RLOVXBBpoI/9aERtUhNKrrWh3Pu74C/iySJJF4myJXzzh6G3X/s3c3l1AFvneaoRWqiizJJYFN3d1nwOOp9D8Ozfwy9G7yCvvYj0L68TglF0klFLYFNXzlvoeOoc+NeMX/+5TqlEkk/zVFLYGOT/t1dOhYo6kJGUx0iIVNRS2CjkwUAli9rnf9F+ax2HoqETEUtgY1OeEXd27FAUReyGlGLhExFLYGN+lMfvQuOqDO6/rRIyFTUEtjUiHrBE140ohYJnYpaAhudLNDT3kJzk83/Is1Ri4RORS2BjU4UF572AP+oD019iIRJRS2BjU4W6Fno0DzQiFokAipqCWx0ohBgRJ3ViFokZCpqCWx0srjwoXnlslfUGlGLhEpFLYF5I+oFpj6Kk95SR32IhErX+pDARicL84+ox05euFqejqMWCZWKWgIplMqMTRbp65yjqH/2Nfi//+7C89519QsmsgSoqCWQ89k8ACu72i7deOYN74p5v/JN6L8SBrbUOZ1IuqmoJZDzGe+sxP65ijqfgc5+2PJLdU4lsjRoZ6IEcjaTA2DFXEVdyEJbd50TiSwdKmoJZGpEPWdR58d1SJ5IhFTUEsi5hUbUeZ3kIhIlFbUEcjbj7Uzs75xnjlpFLRIZFbUEcj6Tp7ejhdbmOX5kChlNfYhESEUtgZzN5Oee9gCNqEUipqKWQE6P5RjoaZ97o+aoRSKlopZATo5Msnb5sks3lMu6BrVIxFTUsqhy2flF3XHpxuKEt9QctUhkVNSyqHPZPPlSee6izme8pUbUIpFRUcuiTgx7ly9dM3vqo5iHE/u9xypqkcjoWh+yqBMj3vTGur4ZI2rn4L++D4bf9p4v3xBDMpGlQUUti3p32CvqNTOnPvLjXklv+yR84Auw4sqY0omkn6Y+ZFFvncnQ3d7CQPeMw/MmR73lhkGVtEjEVNSyqMNDGa4a6MLMLqycHPGWHcvjCSWyhKioZVGHh8bZPDDrMqYqapG6UVHLgs6M5zgxMsn1a3sv3jBV1O0qapGoqahlQS8eGwZg6+V9F2/I+XPUGlGLRE5FLQv66aGztDU3ceP6eUbUKmqRyKmoZV7OOf72tdP83OYVdLbNOpJzcthbdvRe8nUiEq5Fi9rMOszsWTPbb2Yvm9mX6xFM4veTQ2d560yGj25dd+nGyRFoWQYt81xRT0RCE+SElxxwh3Nu3MxagafN7DHn3DMRZ5MYHR+e4EuPHmD98g4+dk07HNsLJ1+Eky/BiX1w4kXoXh13TJElYdGids45YNx/2ur/cVGGCpNzjtee/zHD+/43LWPHsGKOJoo0uRI4N/0fYv4j89fZ9JYL65n12ultTH0N3qnV877Of+xmbndzv3b6bWZvn5UXh3MXb7dLci32njPzAa5Mvljme5ZnXfMITf85f+Ebt/fCmpvg/Q/AdR9FRKIX6BRyM2sG9gJXA19zzu2e4zU7gZ0AV1xxRZgZqzaZHee1r32CmzM/oeyMM9ZPzjooWTNlmnGY305TFWoXff2F5/52m67aOb5mxrY533Pm13rPHeCsacbz+d/D7JI6vfi9zX/drP8ONzuHzZPHZnwvoKejlf6BPppWboDedbD8clhzI/RtnPEeIlIPgYraOVcCbjazPuARM7vROXdg1mt2AbsABgcHEzHiPrDr02wf/ynPbH6AG+79PJf1rYo7kohIxSo66sM5Nww8BdwdSZoQvfbc3zA4/Di719/Pjvt/nx6VtIg0qCBHfQz4I2nMbBlwF/BaxLlqlv1/f8R5etj6T3437igiIjUJMqJeCzxlZi8CzwFPOOd+GG2s2gyfOcl7x57mtdUfobNbJ2SISGMLctTHi8C2OmQJzeGfPsL7rMSKn/vVuKOIiNQsnWcmvvEjztDHNVtvizuJiEjNUlfUzjnWjb3I2z3baGpujjuOiEjNUlfU777zFms5Q2n9YNxRRERCkbqiPnHgxwCseI+mPUQkHVJX1Pmju8m5Vq64YUfcUUREQpG6ou4/t48jbVfT2tax+ItFRBpAqoq6XCqxMX+Y8303xh1FRCQ0qSrqk2+/QaflsNU3xB1FRCQ0qSrqocMvALB8400xJxERCU+qinri+EsArNvSUCdSiogsKFVF3Xr2ICcYoHf5irijiIiEJlVF3Z85xKmOzXHHEBEJVWqKulTIs6F4jGz/lrijiIiEKjVF/e6bL9NmRZpXXxd3FBGRUKWmqM+8tR+Afh3xISIpk5qizr37CmVnrL9ma9xRRERClZqibj3/OieaLqOruzfuKCIioUpNUa/IvMnpjivjjiEiErpUFHWxkGd96RiTfdfEHUVEJHSpKOqTR16hzUo06YgPEUmhVBT1mTe9Iz76dMSHiKRQKoo6d0JHfIhIeqWiqFvPHuRE02V064gPEUmhVBR1f+YthnTEh4ikVMMX9eREhvWld8j1vyfuKCIikWj4oj768s9osxKtG2+JO4qISCQavqhH33gGgLU33BpzEhGRaDR8UbecfIGTrGTN+k1xRxERiUTDF/Xq0QO803k9ZhZ3FBGRSDR0UQ+dOsY6d5LS2u1xRxERiUxDF/Xbzz0GwMobfiHmJCIi0Wnooi4ffJxz9LJ56+1xRxERiUzDFnV2fJjrR5/m8IrbaW5ujjuOiEhkGrao9/+fXXTZJF23/HrcUUREItWQRT185hRXvfp1DrZex3W33BV3HBGRSLXEHaASrlzm8Mu7KT76OTa7EUbv+S7W1JD/1oiIBJaooj70le20lScxHEYZc/4SRxNlelyGqy3HuFvGyz//ENu2aSeiiKTfokVtZpcD3wNWAw7Y5Zz7ahRhhrs2Y64I1jTjj00/LjR30bT6Orbcfh/bBtZGEUFEJHGCjKiLwBecc8+bWQ+w18yecM69EnaYwc//IOy3FBFpeItO8DrnTjjnnvcfjwGvAuujDiYiIp6K9sSZ2SZgG7B7jm07zWyPme0ZGhoKKZ6IiAQuajPrBv4H8Dnn3Ojs7c65Xc65Qefc4MDAQJgZRUSWtEBFbWateCX9sHPuf0YbSUREZlq0qM27fuh/A151zv1h9JFERGSmICPqW4FPAneY2T7/zz0R5xIREd+ih+c5554GdFV+EZGY6PxrEZGEM+dc+G9qNgQcDf2NF7cKOBPD961VI+ZuxMyg3PXUiJkhvtwbnXNzHjIXSVHHxcz2OOcG485RqUbM3YiZQbnrqREzQzJza+pDRCThVNQiIgmXtqLeFXeAKjVi7kbMDMpdT42YGRKYO1Vz1CIiaZS2EbWISOqoqEVEEi41RW1md5vZQTM7ZGYPxp1nPmZ2xMxe8k/F3+OvW2FmT5jZG/6yPwE5v21mp83swIx1c+Y0zx/5n/2LZrY9Ybl/x8yOz3UJBDP7LT/3QTP7pZgyX25mT5nZK2b2spl91l+f2M97gcxJ/6w7zOxZM9vv5/6yv/5KM9vt5/tLM2vz17f7zw/52zfFkRvnXMP/AZqBw8BmoA3YD1wfd655sh4BVs1a9x+BB/3HDwL/IQE5bwe2AwcWywncAzyGd6mBHcDuhOX+HeCLc7z2ev9npR240v8Zao4h81pgu/+4B3jdz5bYz3uBzEn/rA3o9h+34l1bfwfwV8B9/vpvAJ/xH/8L4Bv+4/uAv6x3ZudcakbUtwCHnHNvOufywPeBe2POVIl7ge/6j78L/KP4oniccz8Gzs1aPV/Oe4HvOc8zQJ+ZxXJTy3lyz+de4PvOuZxz7i3gEN7PUl25+e+ilNjPe4HM80nKZ+2cc+P+01b/jwPuAKbuBTj7s576f/AD4E7/iqJ1lZaiXg+8M+P5MZJ7uzAH/MjM9prZTn/daufcCf/xSbwbCSfRfDkb4fP/l/40wbdnTC0lLvesuyg1xOc9x52fEv1Zm1mzme0DTgNP4I3uh51zxTmyTef2t48AK+samPQUdSO5zTm3HfgQ8ICZ3T5zo/N+x0r8MZONktP3deAq4GbgBPAHsaaZx0J3UUrq5z1H5sR/1s65knPuZmAD3qj+2ngTLS4tRX0cuHzG8w3+usRxzh33l6eBR/B+UE5N/erqL0/Hl3BB8+VM9OfvnDvl/+UsA9/kwq/cick9z12UEv15z5W5ET7rKc65YeAp4P1400dTl32emW06t799OXC2vknTU9TPAdf4e27b8Cb9/zrmTJcwsy4z65l6DHwQOICX9X7/ZfcDj8aTcFHz5fxr4Nf9oxF2ACMzfmWP3az521/G+8zBy32fv2f/SuAa4NkY8s13F6XEft7zZW6Az3rAzPr8x8uAu/Dm158CPu6/bPZnPfX/4OPAk/5vN/UVxx7MKP7g7Ql/HW++6bfjzjNPxs14e773Ay9P5cSb8/pb4A3gb4AVCcj6F3i/uhbw5uw+PV9OvD3pX/M/+5eAwYTl/lM/14t4f/HWznj9b/u5DwIfiinzbXjTGi8C+/w/9yT5814gc9I/65uAF/x8B4Av+es34/3DcQj470C7v77Df37I3745jtw6hVxEJOHSMvUhIpJaKmoRkYRTUYuIJJyKWkQk4VTUIiIJp6IWEUk4FbWISML9f8VpJdVTh9sSAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAa3ElEQVR4nO3de5Cd9X3f8fd379qbdiUtuoKEAJmbEZK3RA6YTiA4GF9oGo+HdOLQ2hNNXdratT0NaWY8cTzJtJ0mxWk9dmTXjp2QOKlbSuoO1CSQOthGIIEE4iKQQAIJXVaXvZ2ze66//vE8u1qt9vKcc57nPM959vOa0TznPM/Zsx8Oq49++3tu5pxDRESSqynuACIisjAVtYhIwqmoRUQSTkUtIpJwKmoRkYRrieJNV61a5TZt2hTFW4uIpNLevXvPOOcG5toWSVFv2rSJPXv2RPHWIiKpZGZH59umqQ8RkYRTUYuIJJyKWkQk4VTUIiIJp6IWEUk4FbWISMKpqEVEEi6S46hlaSuUyjzz5lkOnx4nky+RK5a9DTMuqfueNb18+Ka1MSWUWDg342dg1uOp7YEez/H1cz5mnq+f63EZ8lnIj0M+4/8Zn/F8HEqFOd53lrYuuO1zC38OVVBRS6gOD43zqT95jqNns3NuN/OWH71pnYo6rYZeh2/9IuRGmS7C1LKLn3ZfpqKWZJvIl7j/288yWSjxjV/bzuCmFXS3t9De0oSZLf4GEr2f/hc4/GS03+P8USgX4fYvAnbhX+cFHzP3evCfV/p4oe8312sM2rq9EXFb16zH/vOWtio/kNqpqCU0f7XnHY6dn+AvfmMH779qZdxxZLbCBDz5e9C5Enoj/G2mcwXs+Azc8hvRfY8lRkUtoflf+45z3dpelXRUnn4Ijj1X/ddPjkBxAj72Vbj6F0OLJdFTUUsoToxM8MLbw/zm3dfGHSWdChPw5Fe80XDnqurf56o7YdMHwssldaGillA8f3QYgFuv1mg6Eif2e/O+H3kIrr0n7jRSZypqCcX+Y8O0tTRx7ZreuKM0hqd+H579ZvDXl/LecsM/iCaPJJqKWkKx751hbljXS1uLzqEK5JVHYVk/XHVH8K9ZeTV0z3ldeUk5FbWE4vVTY9zzXh0XHUgxB2cPwa2fhTu/FHcaaQAqaqnZyESB4WyBTSs7444Sr/NH4MhPFn9d5rQ333zZ9ZFHknRQUUvN3vbPQrxiRVfMSWL22IPw+mPBXmtNsH57tHkkNVTUUrOj5zIAbFzqI+qzb8A1H4R7/tPir23rhi4dISPBqKilZkenR9RLuKjLJe/U6Ws/DP0b404jKaNd9FKzY+ezrOxqo6t9Cf+7P3ocygVYsTnuJJJCgf5mmVkf8C3gRrzLYX3KOfezCHNJAzk9mmN1b0fcMaL16ANwcIH551LRW/ZfWZ88sqQEHQJ9FXjcOfdxM2sDlvDvuDLb0HiOgZ72uGNE6+Bj0LMWrtgx/2vaexbeLlKlRYvazJYDtwP/FMA5lwfy0caSRjI0lmPL6p64Y0Qnn4HsWXj/A/CBL8SdRpagIHPUVwJDwHfM7AUz+5aZXXIclpntNLM9ZrZnaGgo9KCSTOWyY2gs5SPqkWPecvkV8eaQJStIUbcA24GvO+e2ARngwdkvcs7tcs4NOucGBwZ0mutSMTxRoFh2XJbmoh5+x1v2XR5vDlmygsxRHwOOOed2+89/wBxFLUvT0FgOoPFH1Lkx+NN/DBPn5t4GsHxDfTOJ+BYtaufcSTN7x8ze45w7CNwJvBJ9NGkE00Xd3eBFffpVOPasd63m7ssu3b58A/Sur38uEYIf9fGvgIf9Iz7eBP5ZdJGkkZzNeEW9stGLevS4t7z738OaG+PNIjJLoKJ2zu0DBqONIo1oZKIAQF9na8xJajT6rrfsXRdvDpE56MxEqclw1ivq5ctSUNStnd41okUSZgmf8ythGJko0NXWTGtzg/ybf/4InDl06foT+73RtFndI4ksRkUtNRnOFhprNP3wJ+DMwbm3bbm7vllEAlJRS01GJgr0NkpROwfDb8N7PwG37Lx0+8CW+mcSCUBFLTUZnSg0zo7E/DgUJ7yjOi7XTWKlcTTIxKIk1fBEvnGmPsZPe8vu1fHmEKmQilpqMjJRoG9ZW9wxgpkq6i5d4kAai4paajIyUWB5o0x9jJ/ylhpRS4PRHLVUbbJQYrJQTubUx9//IZw6cPG6c295SxW1NBgVtVRtdCKhJ7uUy/DkV6CjDzpXXLzt6rugUzeVlcaiopaqDSe1qCeHwZXhH/4m7PjncacRqZnmqKVqib3OR+aMt9TIWVJCRS1VG/Gv89HbkbCizp71lrOnPUQalIpaqpbJe3fe7u5I2AzaVFF3rYo3h0hIVNRStUyuBEB3e9KKWlMfki4qaqla1h9Rd7Y1x5xklqkR9TJNfUg6JGwoJI1kakTd2Rbjj9HBx+EHn4Jy8cK6cgFau6CtM75cIiFSUUvVsvkiHa1NNDfFeA3n43uhkIVbP3vx+jXvjSePSARU1FK1TL5IV5yjafCmOTpXwF1fjjeHSIQ0Ry1Vy+ZKdLbHPD89cU5z0ZJ6KmqpWnJG1Dq6Q9JNRS1Vy+ZL8R/xkT2nopbUU1FL1cZzRbriPoZ6ao5aJMW0M1Gqls2VuKynvb7f1DkoTEw9UVHLkqCilqpl8jGMqH/4b2Dvdy5epzu2SMqpqKVq2Xyp/jsTT74EK6+GbZ/0nje1wNZfrW8GkTpTUUvVMrli/Q/PmzgPa7fCbZ+r7/cViZF2JkpViqUyuWK5/iPqifOak5YlR0UtVckWpq7zUccRdbns3b1lWX/9vqdIAqiopSpZ/4JMdd2ZmBv1brGlopYlRkUtVcnEcYnTifPeUkUtS4yKWqqSyXlFXdc5ahW1LFE66kOqMn0t6iiP+hh+Gw4+5p3kAnD2kLfURZhkiVFRS1Wm7u4S6Yj67/8A9v7JxetaOqB/Y3TfUySBVNRSlUy+DjsTs2dh5TXw6R9dWNfSoTu3yJKjopaqZKfmqKOc+pgc9eajddy0LHHamShVmRpRR3q/xNwodCyP7v1FGkSgv2VmdgQYA0pA0Tk3GGUoSb6pEXWkh+dNjsCKzdG9v0iDqGQ49AvOuTORJZGGksmXaGtporU5wl/KJjWiFgFNfUiVsvkiXVGOpp3zpj7ae6P7HiINImhRO+BHZrbXzHbO9QIz22lme8xsz9DQUHgJJZHGc8Vo56eLk1DKQ4eKWiRoUd/mnNsOfAh4wMxun/0C59wu59ygc25wYEAXck+7bK4U/REfoKkPEQLOUTvnjvvL02b2CHAL8OMog0myZfIRjKhHjl04+3D0hLdsV1GLLPo3zcy6gCbn3Jj/+IPA70aeTBItm49gRP1nvwJDr128rmdNuN9DpAEFGRKtBh4xs6nX/7lz7vFIU0niZXJFVnSFfIbg2Em47mOw4zPe89ZlsPbmcL+HSANatKidc28CW+uQRRpINl+iO8zTx52D3Bis2gIbfz689xVJAR2eJ1XJ5ovhnuxSyIIrQXtPeO8pkhIqaqlKJlcK94JM00d56HA8kdlU1FKxUtkxUSiFO6LO+UWtE1xELqGilopFci3q3Ji3VFGLXEJFLRXL5iO4u8vkiLfU1IfIJVTUUrFI7pc4PfWhnYkis6mopWLTI+pQ56g19SEyH93hRSo2PaKu9aiPoYOw5zvgynDqgLdOUx8il1BRS8VCG1E//z3Y/XXo6POer38ftGnqQ2Q2FbVULOMf9VHzmYm5UeheA188GEIqkfTSHLVULJubOuqj1qIe085DkQBU1FKx8emjPmqc+lBRiwSiopaKTZ3wUvP1qCdHVdQiAaiopWKZfInWZqOtpcYfn9yYjvIQCUBFLRXLhnW/xNyYjpsWCUBFLRXL5Evh3IFcc9QigaiopWLZfLH2Iz7KZe/wPI2oRRalopaKZXIhjKgLGcBpRC0SgE54kYplq70DeTEPT/0eTA5DYcJbp6IWWZSKWiqWyZVY19dW+Reeegl+8pB3ynhLOyy/AtbdHHI6kfRRUUvFMvkiXdVcizo37i3vexg23RZuKJEU0xy1VCyTK1U39ZHPeMu2rnADiaSciloqls0Xq9uZmPdH1LpCnkhFVNRSkXLZkc2Xqjs8b7qoNaIWqYSKWioyUfCunFfViDqnohaphopaKjJ1LerqRtSaoxaphopaKjJ1Leqq56hbO6EpxHstiiwBKmqpSKaWS5zmx6GtO+REIumnopaKZKZG1NUcR53PaNpDpAoqaqlITXcgz41Du0bUIpXSmYlSkanbcAW+sW25DMVJcP7V8jT1IVIxFbVUZOo2XIFH1H/+CTj0xIXnW+6OIJVIuqmopSLj/hx1d9CdiadfgXXb4YZfBjO46s4I04mkk4paKnJhjjrgzsR8Bta/D2791xGmEkk37UyUimRyRdpbmmhpDvijU5iAts5oQ4mknIpaKjKeK1awI7EEpRy06pA8kVqoqKUimVwx+I7E6VPGNaIWqUXgojazZjN7wcx+GGUgSbbxXCl4URey3rJVRS1Si0pG1J8FXo0qiDSGTK5IdyU7EkFnI4rUKFBRm9kG4MPAt6KNI0mXqeTGthpRi4Qi6Ij6IeDfAuX5XmBmO81sj5ntGRoaCiObJFBFOxOn7jSuohapyaJFbWYfAU475/Yu9Drn3C7n3KBzbnBgYCC0gJIs3s7ESqc+VNQitQgyor4V+JiZHQG+D9xhZn8WaSpJrIx2JorU3aJF7Zz7LefcBufcJuA+4Enn3K9FnkwSxzlHJl/B1EfeL2rtTBSpiY6jlsCy+RLOVXBBpoI/9aERtUhNKrrWh3Pu74C/iySJJF4myJXzzh6G3X/s3c3l1AFvneaoRWqiizJJYFN3d1nwOOp9D8Ozfwy9G7yCvvYj0L68TglF0klFLYFNXzlvoeOoc+NeMX/+5TqlEkk/zVFLYGOT/t1dOhYo6kJGUx0iIVNRS2CjkwUAli9rnf9F+ax2HoqETEUtgY1OeEXd27FAUReyGlGLhExFLYGN+lMfvQuOqDO6/rRIyFTUEtjUiHrBE140ohYJnYpaAhudLNDT3kJzk83/Is1Ri4RORS2BjU4UF572AP+oD019iIRJRS2BjU4W6Fno0DzQiFokAipqCWx0ohBgRJ3ViFokZCpqCWx0srjwoXnlslfUGlGLhEpFLYF5I+oFpj6Kk95SR32IhErX+pDARicL84+ox05euFqejqMWCZWKWgIplMqMTRbp65yjqH/2Nfi//+7C89519QsmsgSoqCWQ89k8ACu72i7deOYN74p5v/JN6L8SBrbUOZ1IuqmoJZDzGe+sxP65ijqfgc5+2PJLdU4lsjRoZ6IEcjaTA2DFXEVdyEJbd50TiSwdKmoJZGpEPWdR58d1SJ5IhFTUEsi5hUbUeZ3kIhIlFbUEcjbj7Uzs75xnjlpFLRIZFbUEcj6Tp7ejhdbmOX5kChlNfYhESEUtgZzN5Oee9gCNqEUipqKWQE6P5RjoaZ97o+aoRSKlopZATo5Msnb5sks3lMu6BrVIxFTUsqhy2flF3XHpxuKEt9QctUhkVNSyqHPZPPlSee6izme8pUbUIpFRUcuiTgx7ly9dM3vqo5iHE/u9xypqkcjoWh+yqBMj3vTGur4ZI2rn4L++D4bf9p4v3xBDMpGlQUUti3p32CvqNTOnPvLjXklv+yR84Auw4sqY0omkn6Y+ZFFvncnQ3d7CQPeMw/MmR73lhkGVtEjEVNSyqMNDGa4a6MLMLqycHPGWHcvjCSWyhKioZVGHh8bZPDDrMqYqapG6UVHLgs6M5zgxMsn1a3sv3jBV1O0qapGoqahlQS8eGwZg6+V9F2/I+XPUGlGLRE5FLQv66aGztDU3ceP6eUbUKmqRyKmoZV7OOf72tdP83OYVdLbNOpJzcthbdvRe8nUiEq5Fi9rMOszsWTPbb2Yvm9mX6xFM4veTQ2d560yGj25dd+nGyRFoWQYt81xRT0RCE+SElxxwh3Nu3MxagafN7DHn3DMRZ5MYHR+e4EuPHmD98g4+dk07HNsLJ1+Eky/BiX1w4kXoXh13TJElYdGids45YNx/2ur/cVGGCpNzjtee/zHD+/43LWPHsGKOJoo0uRI4N/0fYv4j89fZ9JYL65n12ultTH0N3qnV877Of+xmbndzv3b6bWZvn5UXh3MXb7dLci32njPzAa5Mvljme5ZnXfMITf85f+Ebt/fCmpvg/Q/AdR9FRKIX6BRyM2sG9gJXA19zzu2e4zU7gZ0AV1xxRZgZqzaZHee1r32CmzM/oeyMM9ZPzjooWTNlmnGY305TFWoXff2F5/52m67aOb5mxrY533Pm13rPHeCsacbz+d/D7JI6vfi9zX/drP8ONzuHzZPHZnwvoKejlf6BPppWboDedbD8clhzI/RtnPEeIlIPgYraOVcCbjazPuARM7vROXdg1mt2AbsABgcHEzHiPrDr02wf/ynPbH6AG+79PJf1rYo7kohIxSo66sM5Nww8BdwdSZoQvfbc3zA4/Di719/Pjvt/nx6VtIg0qCBHfQz4I2nMbBlwF/BaxLlqlv1/f8R5etj6T3437igiIjUJMqJeCzxlZi8CzwFPOOd+GG2s2gyfOcl7x57mtdUfobNbJ2SISGMLctTHi8C2OmQJzeGfPsL7rMSKn/vVuKOIiNQsnWcmvvEjztDHNVtvizuJiEjNUlfUzjnWjb3I2z3baGpujjuOiEjNUlfU777zFms5Q2n9YNxRRERCkbqiPnHgxwCseI+mPUQkHVJX1Pmju8m5Vq64YUfcUUREQpG6ou4/t48jbVfT2tax+ItFRBpAqoq6XCqxMX+Y8303xh1FRCQ0qSrqk2+/QaflsNU3xB1FRCQ0qSrqocMvALB8400xJxERCU+qinri+EsArNvSUCdSiogsKFVF3Xr2ICcYoHf5irijiIiEJlVF3Z85xKmOzXHHEBEJVWqKulTIs6F4jGz/lrijiIiEKjVF/e6bL9NmRZpXXxd3FBGRUKWmqM+8tR+Afh3xISIpk5qizr37CmVnrL9ma9xRRERClZqibj3/OieaLqOruzfuKCIioUpNUa/IvMnpjivjjiEiErpUFHWxkGd96RiTfdfEHUVEJHSpKOqTR16hzUo06YgPEUmhVBT1mTe9Iz76dMSHiKRQKoo6d0JHfIhIeqWiqFvPHuRE02V064gPEUmhVBR1f+YthnTEh4ikVMMX9eREhvWld8j1vyfuKCIikWj4oj768s9osxKtG2+JO4qISCQavqhH33gGgLU33BpzEhGRaDR8UbecfIGTrGTN+k1xRxERiUTDF/Xq0QO803k9ZhZ3FBGRSDR0UQ+dOsY6d5LS2u1xRxERiUxDF/Xbzz0GwMobfiHmJCIi0Wnooi4ffJxz9LJ56+1xRxERiUzDFnV2fJjrR5/m8IrbaW5ujjuOiEhkGrao9/+fXXTZJF23/HrcUUREItWQRT185hRXvfp1DrZex3W33BV3HBGRSLXEHaASrlzm8Mu7KT76OTa7EUbv+S7W1JD/1oiIBJaooj70le20lScxHEYZc/4SRxNlelyGqy3HuFvGyz//ENu2aSeiiKTfokVtZpcD3wNWAw7Y5Zz7ahRhhrs2Y64I1jTjj00/LjR30bT6Orbcfh/bBtZGEUFEJHGCjKiLwBecc8+bWQ+w18yecM69EnaYwc//IOy3FBFpeItO8DrnTjjnnvcfjwGvAuujDiYiIp6K9sSZ2SZgG7B7jm07zWyPme0ZGhoKKZ6IiAQuajPrBv4H8Dnn3Ojs7c65Xc65Qefc4MDAQJgZRUSWtEBFbWateCX9sHPuf0YbSUREZlq0qM27fuh/A151zv1h9JFERGSmICPqW4FPAneY2T7/zz0R5xIREd+ih+c5554GdFV+EZGY6PxrEZGEM+dc+G9qNgQcDf2NF7cKOBPD961VI+ZuxMyg3PXUiJkhvtwbnXNzHjIXSVHHxcz2OOcG485RqUbM3YiZQbnrqREzQzJza+pDRCThVNQiIgmXtqLeFXeAKjVi7kbMDMpdT42YGRKYO1Vz1CIiaZS2EbWISOqoqEVEEi41RW1md5vZQTM7ZGYPxp1nPmZ2xMxe8k/F3+OvW2FmT5jZG/6yPwE5v21mp83swIx1c+Y0zx/5n/2LZrY9Ybl/x8yOz3UJBDP7LT/3QTP7pZgyX25mT5nZK2b2spl91l+f2M97gcxJ/6w7zOxZM9vv5/6yv/5KM9vt5/tLM2vz17f7zw/52zfFkRvnXMP/AZqBw8BmoA3YD1wfd655sh4BVs1a9x+BB/3HDwL/IQE5bwe2AwcWywncAzyGd6mBHcDuhOX+HeCLc7z2ev9npR240v8Zao4h81pgu/+4B3jdz5bYz3uBzEn/rA3o9h+34l1bfwfwV8B9/vpvAJ/xH/8L4Bv+4/uAv6x3ZudcakbUtwCHnHNvOufywPeBe2POVIl7ge/6j78L/KP4oniccz8Gzs1aPV/Oe4HvOc8zQJ+ZxXJTy3lyz+de4PvOuZxz7i3gEN7PUl25+e+ilNjPe4HM80nKZ+2cc+P+01b/jwPuAKbuBTj7s576f/AD4E7/iqJ1lZaiXg+8M+P5MZJ7uzAH/MjM9prZTn/daufcCf/xSbwbCSfRfDkb4fP/l/40wbdnTC0lLvesuyg1xOc9x52fEv1Zm1mzme0DTgNP4I3uh51zxTmyTef2t48AK+samPQUdSO5zTm3HfgQ8ICZ3T5zo/N+x0r8MZONktP3deAq4GbgBPAHsaaZx0J3UUrq5z1H5sR/1s65knPuZmAD3qj+2ngTLS4tRX0cuHzG8w3+usRxzh33l6eBR/B+UE5N/erqL0/Hl3BB8+VM9OfvnDvl/+UsA9/kwq/cick9z12UEv15z5W5ET7rKc65YeAp4P1400dTl32emW06t799OXC2vknTU9TPAdf4e27b8Cb9/zrmTJcwsy4z65l6DHwQOICX9X7/ZfcDj8aTcFHz5fxr4Nf9oxF2ACMzfmWP3az521/G+8zBy32fv2f/SuAa4NkY8s13F6XEft7zZW6Az3rAzPr8x8uAu/Dm158CPu6/bPZnPfX/4OPAk/5vN/UVxx7MKP7g7Ql/HW++6bfjzjNPxs14e773Ay9P5cSb8/pb4A3gb4AVCcj6F3i/uhbw5uw+PV9OvD3pX/M/+5eAwYTl/lM/14t4f/HWznj9b/u5DwIfiinzbXjTGi8C+/w/9yT5814gc9I/65uAF/x8B4Av+es34/3DcQj470C7v77Df37I3745jtw6hVxEJOHSMvUhIpJaKmoRkYRTUYuIJJyKWkQk4VTUIiIJp6IWEUk4FbWISML9f8VpJdVTh9sSAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# decay to backing solvent\n", "poly1_l = poly1(100., 3.)\n", "d2o_l = d2o(0, 3)\n", "\n", "ff = FunctionalForm(200., 30., 5., poly1_l, d2o_l)\n", "\n", "struct = si | poly1_l | d2o(0, 3)\n", "struct2 = si | poly1_l | ff | d2o_l\n", "\n", "struct2.contract = 1.\n", "\n", "plt.plot(*struct.sld_profile())\n", "plt.plot(*struct2.sld_profile());" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAdXElEQVR4nO3deZCcd33n8fe3e+7RHNLMaDQ6R7J12NiWbQbb2C4WbAg2GG/CQi1HEsiS1W4t4UyKhUrtbrEVqvYIIQ4ssMqGABsg4UqCHS4Hm4WkwEYyMtiWL8mXpJHmnp7u6WO6+7d/PD3SzHikeQZN93P051Wlerqf55nu7zzS89Gvf/17fo855xARkfBKBF2AiIhcmIJaRCTkFNQiIiGnoBYRCTkFtYhIyDVU40V7e3vd4OBgNV5aRCSWDh8+POac61tuW1WCenBwkEOHDlXjpUVEYsnMnjvfNnV9iIiEnIJaRCTkFNQiIiGnoBYRCTkFtYhIyCmoRURCTkEtIhJyVRlHLRIlT5ye4cFnJ5ieLVAqQ6lcPrtt76ZOXn/VQIDVSejMnIZnfgzTz0O5BK7sLQGa2uHm96/5Wyqopa59/PtP8Mn7nn7RejNveef+zQpqOefJ78FX3wHF7DIbDdZtVFCLrKWfHBvnk/c9zb+6disfeM1u+jqaaUwkSCQs6NIkjLJT8I3fhb498IY/g769kGiERPLc/+xVoqCWuvWZ/3eMjR3NfOw3rqClMRl0ORJ2D30B8im481MwcFVN31pfJkpdyuSL/PTYOP/y6s0KafHnie/AwP6ahzQoqKVO/fT4OIVSmVfu3Rh0KRIFuWl44UG49NWBvL2CWurSwyemSRi8dMf6oEuRKDj9S3Al2P7yQN5eQS116cnTMwz2tqvbQ/w585i37H9JIG+voJa69OSZGfb2dwRdhkTFyKPQ0g0dwQzVVFBL3cnNlXh2PMNuBbX4NfokbLy86sPwzkdBLXXn1FSWsoPBnragS5GomHoe1u8I7O19BbWZdZvZ183scTM7ambB9KiLrIHh6RwAA12tAVcikVCag5lT0LU1sBL8XvByF/Bd59ybzKwJUFNEIuvklHf575ZuBbX4kDrlzefRtS2wElYMajPrAl4BvBPAOVcACtUtS6R6hqe8FnV/V3PAlUgkTJ/wlt3BBbWfro+dwCjwl2b2czP7P2bWXuW6RKrm1FSWvo5mmhs0NE98mA/qAFvUfoK6AbgW+Ixz7hogA3x46U5mdsDMDpnZodHR0TUuU2TtDKdyDHS1BF2GREX6jLdcF9xVrH6C+gRwwjn3QOX51/GCexHn3EHn3JBzbqivr28taxRZU+PpPH3r1O0hPmVGIdkMzZ2BlbBiUDvnTgMvmNneyqpbgceqWpVIFY2nC/Ssawq6DImKzKjXmg5oDDX4H/XxHuBLlREfx4HfqV5JItXjnGM8k2dDu1rU4lN6BNqD7SXwFdTOuSPAUHVLEam+VK7IXMnRqxa1+JUZgc4tgZagKxOlrkxkvJGl6voQ39KjgbeoFdRSV8bTeQB61PUhfjgHs2MKapFaGkt7LeoN7WpRiw/5GSgXoW1DoGUoqKWuzHd9KKjFl+yEt2wN9gYTCmqpK6ncHADdbY0BVyKRkJ30lq1qUYvUzHR2joaE0ao7u4gfs2pRi9TcdHaOrtZGLMCLFyRC5lvU6qMWqZ35oBbx5WzXh1rUIjWTys7RqaAWvxTUIrWnoJZVmZ2Apg5IBvtvRkEtdUVdH7IquWlo7Q66CgW11BcvqP3ORSZ1LzcNLV1BV6GglvrhnCOVK6pFLf4pqEVqK1MoUSo7BbX4p6AWqa3prHdVYmeLglp8UlCL1Nb0rBfUalGLbwpqkdqan+dDQS2+lMuQTymoRWrpbNeHglr8KMwATkEtUkvzQa0WtfiSm/aWCmqR2kmpRS2rMR/UzZ3B1oGCWupIKjuHGXQ064IX8UEtapHam87O0dnSSCKhKU7FBwW1SO1NZ+fo1OXj4peCWqT2dPm4rIqCWqT2NHOerIq+TBSpvVSlj1rEl9x0ZS7q4LvLFNRSN1I5BbWsQi4cVyWCglrqSCpbpKtNQS0+5aZCE9S+2vRm9iwwA5SAonNuqJpFiay1QrFMdq5EZ0vwH2MlIkIyIRP4DOqKVznnxqpWiUgVzU/IpKsSxbfcFHRuCboKQF0fUidSmotaVis7Hfjdx+f5DWoHfN/MDpvZgeV2MLMDZnbIzA6Njo6uXYUia0ATMsmq5aagpTvoKgD/QX2zc+5a4Hbg3Wb2iqU7OOcOOueGnHNDfX19a1qkyMVK5YoAujJR/CkVvbmoQ3AHcvAZ1M65k5XlCPC3wHXVLEpkranrQ1bl7FWJ3YGWMW/FoDazdjPrmH8M/BrwSLULE1lL6vqQVclNecuQtKj9fA7sB/7WzOb3/7Jz7rtVrUpkjWnUh6xKdspbhqRFvWJQO+eOA/trUItI1aSyRZqSCZobNNBJfMhNesuIjfoQiTRvitNGKp8MRS5svkUdkq4PBbXUhVROc1HLKmQrLeqQdH0oqKUuaOY8WZWQfZmooJa6kKp0fYj4kp2ChlZoaA66EkBBLXVCd3eRVclNheaLRFBQS53wuj7URy0+ZadC0+0BCmqpA865ypeJalGLT9mp0HyRCApqqQOZQom5kqNbQS1+5abUohappYl0AYAN7U0BVyKRoRa1SG1NzCqoZZXUohaprYlMHlBQi0/FAhTSGvUhUkvj6vqQ1Zit3HGwvTfYOhZQUEvsTarrQ1YjU7lDVXt4boCioJbYG88UaEwa65o1jlp8UFCL1N5kpsCG9ibNnCf+ZOa7PhTUIjUzli7Q0x6OORskAtIj3lJBLVI7p6dzbOpqCboMiYrMKCSbobkj6ErOUlBL7J1J5ejvVFCLT+kRWLcRQtRVpqCWWMsXS4xnCmxSUItfqZPQuSXoKhZRUEusjaS8i102damPWnyafgG6tgZdxSIKaom1M6kcgLo+xJ9yGaZPKqhFaunUtBfUA12tAVcikZAZgfKcglqklo6PpjGDHT1tQZciUTD1grdUUIvUzrHRDFu6W2lpTAZdikTB2BPesmd3sHUsoaCWWDs+mmZX37qgy5CoGDnqjaHesDPoShZRUEtslcqO46MZdvW2B12KRMXo49C3BxLh+gSmoJbYeuxUiuxciWu2dwddikRBuQynjkD/FUFX8iIKaomtB54ZB+D6nT0BVyKRMPKoNxf1zlcEXcmL+A5qM0ua2c/N7J5qFiSyVr7/6Bl29bZrng/x5/Fve8tdrwy0jOWspkX9PuBotQoRWUs/PT7Og89O8LbrtwddikTB7AT87M/h0ldD5+agq3kRXzOpm9lW4PXAx4APVrUiqZrnx2f52bMTTM4WKJYdzoHDAeDc4n2dW7zeLdo2v+78+8yvdOf5meV+7kX7LHjBc/usXNdsocjdD59isKeNf/2ybcgqOAdnHoEzj0J2Eoq5xdvOPVlm/dK/FLfMNudv29nnq9m2TD1+3t+V4dgPIDcNr/pDwsjvLS/+FPgQcN55/8zsAHAAYPt2tWLC5nP/9Ax/9A+PUXYr77tW5icfs7PPbdHzxfss3vlC+5zvdefXNSSNf7G3j4/cfhkdLY0X+2vUj5OH4e/eDaO1/OA8/5dp539+oW1nn/vZdoH37N0Nd34Stlx7Mb9M1awY1GZ2BzDinDtsZq88337OuYPAQYChoaEaxoGs5Ohwio99+yi37NvIh2/fx8bOFhoTXq/XwpkcVxuIi/YJ0ZSQ8is4cRi+cAe09cAb7oIdN3mPG5deer/MP5iF6/2EqP6trJqfFvVNwJ1m9jqgBeg0s79yzv1mdUuTtfLFnzxHS0OCP37zfrrbdINXWSKfhq+907ujybvuhY7+oCuSJVb8MtE59xHn3Fbn3CDwFuA+hXR0lMqOex87wyv3bVRIy/Lu/xhMPw9vPKiQDimNo465p0ZmGEvnuXXfxqBLkTA6cQh++hl42e/C9huCrkbOw++XiQA4534I/LAqlUhVHB1OAXDllq6AK5HQKebhW++BjgG49b8EXY1cwKqCWqLn6PAMTQ0Jdmq+C1nqOx+CkcfgrX8DLZ1BVyMXoK6PmDs6nGJP/zoakvqrlgWO3QeHPw83fxD23hZ0NbICnb0x9/zELDt7Nc2nLDF+zFve8B+CrUN8UVDHWLnsGJ7OsVlzXchSee+7C5rPew2bhIiCOsbGMwUKxTKbu3W/QFkin4ZEAzTo7uxRoKCOseHpLAADalHLUoW015rWVYKRoKCOsVNT3oQ6alHLi+RnoEndHlGhoI4xtajlvPIz0KwvmaNCQR1jE5kCCYP1unRclprv+pBIUFDH2HimwPq2JhIJ9UPKEvkZaFKLOioU1DE2mSmwoV2taVlGPq2ujwhRUMfYeKbAegW1LCc/o66PCFFQx9hkpkCPglqWMzuuUR8RoqCOsQm1qGWph/8aPnszlPKwfkfQ1YhPmj0vpsplx+SsWtSyxI8/DoVZuP1/wEt/J+hqxCcFdUylcnOUnYbmyRL5Gbj0Vrj+3wVdiayCuj5iKpUtAtDZqrtwywL5NDRr7umoUVDHVDrvBfW65mTAlUholMtQ0GiPKFJQx9S5oFaLWioKaW+pC10iR0EdU5lKULerRS3z5oNaLerIUVDH1LkWtb4vlor8jLdUUEeOgjqmzgZ1i4JaKs4Gtb5MjBoFdUyd6/pQUEvF2aBWH3XUKKhjar5F3d6koJYKdX1Els7imErnirQ1JUlqilPJTcPU8zDymPdcQR05CuqYyhSK6vYQz+fvgNO/8B5bElrXB1uPrJrO5JhK50sa8SGe6Rdg92vhZe+Czs3Q0hV0RbJKOpNjKp2bU1ALOAe5FGy6Eva8Nuhq5Fe04peJZtZiZg+a2cNm9qiZfbQWhcnFyeRLuthFoJABV1IrOuL8jPrIA7c45/YDVwO3mdkNVa1KLlo6X1SLWrwvEgFaNHY6ylY8k51zDqhce0pj5Y+rZlFy8RTUAiwIarWoo8zXOGozS5rZEWAEuNc598Ay+xwws0Nmdmh0dHSNy5TVyuQ16kNQUMeEr6B2zpWcc1cDW4HrzOyKZfY56Jwbcs4N9fX1rXGZslpqUQsA+ZS3VFBH2qquTHTOTQH3A7dVpRpZE3OlMvliWS1qOdeiblZQR5mfUR99ZtZdedwKvAZ4vMp1yUXIaOY8KZfhhQfh5GHvuVrUkebnTB4AvmBmSbxg/6pz7p7qliUXQ1OcCsfvh796o/e4oVVBHXF+Rn38ArimBrXIGklr5jxJnfKWb/4CDFwFDbrJcZTpTI6hjOailuykt7zkFo2hjgFNcxpD6XwJ0I1t61puypuASTPlxYKCOobSOXV91L3spDdLnmma2zhQUMeQRn3I2aCWWFBQx5BGfQjZSWjbEHQVskYU1DGkUR+iFnW86EyOoUy+SHNDgsak/h+uK+kR+PHHoZiHiWdg4+VBVyRrREEdQ5rno049+T144LPQ1gONrTB4c9AVyRrR2RxDac2cV5+yE97yfQ9rWF7M6LNxDGXUoq5P2UlINELTuqArkTWmoI4hdX3UqdkJjZ2OKQV1DHldH7oqse5kJzQkL6YU1DHk3dhWLeq6k53SkLyYUlDHUDpfpEMTMtWf2QloVYs6jnQ2x1A6V6S9SX+1dWF2AmbHK4/HYLNmJI4jnc0xUyo7snPq+qgLpSLcdTXkp8+tW7cxsHKkenQ2x0ym4F0+rq6POpCd9EL6mt+EXa/yRntcckvQVUkV6GyOGU1xWkfmuzx2vQqufFOwtUhV6cvEmMloQqb6MX8lYltPsHVI1SmoY2Z+5rwOBXX8zbeoNXY69hTUMaMpTuvIrFrU9UJBHTPnuj50ZWLszbeoNXY69tTsipn5G9t2NDcGXIlURbkEqZPe46nnoaEFmtqCrUmqTkEdM+ncHKAWdWz9wwfh8OfPPe/eHlgpUjsK6pjJFLwWtfqoY2rsKejdAze9z3ve/5Jg65Ga0NkcM+l8kYaE0dygrx9iKTMGffu8i1ykbuhsjpl0rsi6lgZMcxLH0+wYtPcGXYXUWKha1A/f/zXMlUkkk5glSSQSJJNJrKmdTZfup6NL326vJJPXhEyxVS55Q/LaFNT1ZsUz2sy2AV8E+gEHHHTO3VWNYnb/8N20WX7ZbUWX4IG+X+elBz5LQ1NzNd4+FnR3lxjLTgJOLeo65OeMLgK/75x7yMw6gMNmdq9z7rG1LubUG79JqThHuVSmXC5RLntLl50m//h3uX7sm/z0c63c8O8/vdZvHRvpvNf1ITGUGfOWusCl7qx4RjvnhoHhyuMZMzsKbAHWPKgv3X+B29u/5m08cNfbGRr+Cmde+AD923av9dvHwkyuSM+6pqDLkLUy9QJ850NQzEGuMp2pWtR1Z1VfJprZIHAN8MAy2w6Y2SEzOzQ6OrpG5S227c7/RIOVOX7vwaq8fhzM5OboaNHFLrFx/IfwxLe9qxAt4U1jOrA/6Kqkxnx/RjazdcA3gPc751JLtzvnDgIHAYaGhtyaVbjA5p37eLTxCjae/MdqvHwsqI86ZtJnvOW/+R40tgZbiwTGV4vazBrxQvpLzrlvVrekC0tvfQWXlI4zeuZEkGWEVipXpFN91PGRGYWmDoV0nVsxqM0bkPsXwFHn3J9Uv6QL67zs1QCcOqJW9VL5YolCsay7u8RJekS31xJfLeqbgN8CbjGzI5U/r6tyXee188obmHNJZp97KKgSQmsmN38bLvVRx4aCWvA36uOfgNBc5tbS2s7xhu20jz8SdCmhcy6o1aKOtLnsucfpM7DxsuBqkVCI5Bk92rGP3VP/HHQZoTNTmTlPLeoI+9H/hPv+aPG6S14VTC0SGpEM6tKGPWyY+g6pqVE6u/uCLic01KKOgZMPQcdmuP5AZYXBFW8MtCQJXiTP6NaBvXAchp/+JZ1DtwRdTmica1FH8q9VAFKnoP9yuPkDQVciIRLJ2fN6d3hz8E6feDzgSsIlVWlRd6rrI7pmTkPHQNBVSMhEMqgHBvdRdAnmRp4MupRQUddHxJWKkBlRUMuLRPKMbmhq4WSyn6bpZ4IuJVTSlaDWlYkRUS7D+NPgvLvykBkDV4ZOBbUsFtkzeqJlO+uzzwVdRqjM5OZoa0rSkIzkB6X68+BB+O5/fPF63QdRlohsUGc7d3JJ5gilUolkUjdyBa/rQ90eETL6OLR0wRsWTO/e2Aa7NBxPFovsWW09u2k7nWf45DMMbL806HJCYSavmfMiJXUSunfAS34j6Eok5CL7Gbl9kxfOY89r5Me8VFYt6kiZPgldW4OuQiIgskHdt30fAOnTTwVcSXhMZAr0tOumAZGROgGdW4KuQiIgss2vni2XMOeSlMY08mPeRKbASzZ3Bl2GLHX0bvjxx8EtnKbdeXds6VJQy8oiG9SJhkZOJzfSNKORHwDOOSZmC2zQbbjC55dfh7GnYceNi9fv2wZ7Xx9MTRIpkQ1qgMnmLXRldQMBgEzBm4t6Q5uCOnQmn4Vt18Hbvxp0JRJRke2jBsh3bKe/NEypXJU7f0XKZKYAwAb1UYfP1HOwfkfQVUiERbpFbRt20T2S4eSZYbYMbA66nECNK6iDVyp6XxAuVJiF7CSsHwykJImHSAd1a/+l8DiMPvd43Qe1WtQhcPd74ciXlt+24ZLa1iKxEumg7tm+F4CZ4SeB+p7udCydBxTUgTp5GDZfC9f928XrG1pgz2uDqUliIdpBvXUfZWeURzWW+vR0DoD+zpaAK6lTpSKMH4OXvxuuflvQ1UjMRDqoE81tDCc30jJ9LOhSAndqOkdPexMtjZr3pKpKRTjxIJTmFq9Pj0B5Dnr3BFOXxFqkgxpgrGWQnuyzQZcRuFNTWTZ3twZdRvw9/GX41nvOv33TlbWrRepG5IN6tnMXuzMPUSrOkWyo3wmJhqezDPa0B11G/J04BC3d8JYvv3hbcwcMXFXzkiT+Ih/UiU1X0XL6Kzz31C/YcdlLgy4nMMNTOW68pDfoMqKvWIBC+vzbh4/AwH4YvKlmJYlEPqj79rwMjsDIkw/WbVBPZArM5ItsXa+uj4viHHzm5d5dVy7kxvfWph6RisgH9bY9V5N1TZROPBR0KYE5OpwCYN8mTch0USaOeyG9/62w+Zrl97EEXP7rNS1LJPJBnWxo5Onmy+mbOBR0KYE5G9QDHQFXEkLpUXji2969CFcy/LC3vPG90H95desSWYXIBzVAZvONXPnspxk7c4re/vq7QvGxUyn6OprpXdccdCnhc//H4PBf+t+/cwv07atePSK/glgEde81r4dnP83xH32Z3jf/QdDl1FSp7PjRU2Ncv3ND0KVU1zM/8uZvXq0nvwt7boc7PuFv/5YuSER6rjKJoRWD2sw+B9wBjDjnrqh+Sat3yZU38szf76D36P+lVHw/yYZY/P/jy+HnJhlL5/m1l/QHXUr1PPNj+MIbfvWff/VHoXNg7eoRqTE/ifZ54FPAF6tbyq/OEgnGr30PQ4f+gJ/8+e9x7Tv/mObWdUGXVXUTmQL/7TtH6Wpt5NbLFgR1sQDTL9S2mJGj8IP/CnPZtX/t3DS0roff+jtIrPLKy2STrhaUyDPnVp7L2cwGgXv8tqiHhobcoUO1/XLPlcv87FO/zXUTdzPnkozbBkqWpESSshlgq3/NNZzm+sIv5ZZ9uNLPlsoOAzZ1tSy++3j6jDe1Zq1t2AXbbqjOa++9HS6/szqvLRICZnbYOTe03LY16yMwswPAAYDt27ev1cv6f/9Egpf93hf55T/fzczR+2jInMbKRRKUMFfGi7sLh/WyW22F7cv+yIV/aKXXMZ9v1NyQoL+r9cV3Ht90lXdBRkMNx1VbAi69Fdpi3lcuEoDYtKhFRKLsQi1qfb0tIhJyCmoRkZBbMajN7CvAT4C9ZnbCzN5V/bJERGTeil8mOufeWotCRERkeer6EBEJOQW1iEjIKahFREJOQS0iEnK+LnhZ9YuajQLPrfkLL68XGKvRe0WJjsvydFyWp+OyvFoelx3Oub7lNlQlqGvJzA6d72qeeqbjsjwdl+XpuCwvLMdFXR8iIiGnoBYRCbk4BPXBoAsIKR2X5em4LE/HZXmhOC6R76MWEYm7OLSoRURiTUEtIhJykQ1qM7vNzJ4ws6fN7MNB11NrZvY5Mxsxs0cWrNtgZvea2VOV5frKejOzP6scq1+Y2bXBVV49ZrbNzO43s8fM7FEze19lfb0flxYze9DMHq4cl49W1u80swcqv//fmFlTZX1z5fnTle2Dgf4CVWZmSTP7uZndU3keuuMSyaA2syTwv4DbgcuBt5rZ5cFWVXOfB25bsu7DwA+cc7uBH1Seg3ecdlf+HAA+U6Maa60I/L5z7nLgBuDdlX8X9X5c8sAtzrn9wNXAbWZ2A/DfgU845y4FJoH5KYzfBUxW1n+isl+cvQ84uuB5+I6Lcy5yf4CXA99b8PwjwEeCriuA4zAIPLLg+RPAQOXxAPBE5fH/Bt663H5x/gP8PfAaHZdFx6QNeAi4Hu+Ku4bK+rPnFPA94OWVxw2V/Szo2qt0PLbi/ed9C3AP3i1NQ3dcItmiBrYALyx4fqKyrt71O+eGK49PA/2Vx3V3vCofS68BHkDHZf7j/RFgBLgXOAZMOeeKlV0W/u5nj0tl+zTQU9OCa+dPgQ8B5crzHkJ4XKIa1LIC5/23X5djL81sHfAN4P3OudTCbfV6XJxzJefc1XgtyOuAfcFWFDwzuwMYcc4dDrqWlUQ1qE8C2xY831pZV+/OmNkAQGU5UllfN8fLzBrxQvpLzrlvVlbX/XGZ55ybAu7H+0jfbWbzd3la+LufPS6V7V3AeG0rrYmbgDvN7Fngr/G6P+4ihMclqkH9M2B35dvZJuAtwLcCrikMvgW8o/L4HXh9tPPrf7syyuEGYHpBV0BsmJkBfwEcdc79yYJN9X5c+sysu/K4Fa/f/iheYL+pstvS4zJ/vN4E3Ff5JBIrzrmPOOe2OucG8TLkPufc2wnjcQm6M/8ivgR4HfAkXl/bHwZdTwC//1eAYWAOrx/tXXj9ZT8AngL+EdhQ2dfwRskcA34JDAVdf5WOyc143Rq/AI5U/rxOx4WrgJ9XjssjwH+urN8FPAg8DXwNaK6sb6k8f7qyfVfQv0MNjtErgXvCelx0CbmISMhFtetDRKRuKKhFREJOQS0iEnIKahGRkFNQi4iEnIJaRCTkFNQiIiH3/wHmO0sLMsHzMwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAdXElEQVR4nO3deZCcd33n8fe3e+7RHNLMaDQ6R7J12NiWbQbb2C4WbAg2GG/CQi1HEsiS1W4t4UyKhUrtbrEVqvYIIQ4ssMqGABsg4UqCHS4Hm4WkwEYyMtiWL8mXpJHmnp7u6WO6+7d/PD3SzHikeQZN93P051Wlerqf55nu7zzS89Gvf/17fo855xARkfBKBF2AiIhcmIJaRCTkFNQiIiGnoBYRCTkFtYhIyDVU40V7e3vd4OBgNV5aRCSWDh8+POac61tuW1WCenBwkEOHDlXjpUVEYsnMnjvfNnV9iIiEnIJaRCTkFNQiIiGnoBYRCTkFtYhIyCmoRURCTkEtIhJyVRlHLRIlT5ye4cFnJ5ieLVAqQ6lcPrtt76ZOXn/VQIDVSejMnIZnfgzTz0O5BK7sLQGa2uHm96/5Wyqopa59/PtP8Mn7nn7RejNveef+zQpqOefJ78FX3wHF7DIbDdZtVFCLrKWfHBvnk/c9zb+6disfeM1u+jqaaUwkSCQs6NIkjLJT8I3fhb498IY/g769kGiERPLc/+xVoqCWuvWZ/3eMjR3NfOw3rqClMRl0ORJ2D30B8im481MwcFVN31pfJkpdyuSL/PTYOP/y6s0KafHnie/AwP6ahzQoqKVO/fT4OIVSmVfu3Rh0KRIFuWl44UG49NWBvL2CWurSwyemSRi8dMf6oEuRKDj9S3Al2P7yQN5eQS116cnTMwz2tqvbQ/w585i37H9JIG+voJa69OSZGfb2dwRdhkTFyKPQ0g0dwQzVVFBL3cnNlXh2PMNuBbX4NfokbLy86sPwzkdBLXXn1FSWsoPBnragS5GomHoe1u8I7O19BbWZdZvZ183scTM7ambB9KiLrIHh6RwAA12tAVcikVCag5lT0LU1sBL8XvByF/Bd59ybzKwJUFNEIuvklHf575ZuBbX4kDrlzefRtS2wElYMajPrAl4BvBPAOVcACtUtS6R6hqe8FnV/V3PAlUgkTJ/wlt3BBbWfro+dwCjwl2b2czP7P2bWXuW6RKrm1FSWvo5mmhs0NE98mA/qAFvUfoK6AbgW+Ixz7hogA3x46U5mdsDMDpnZodHR0TUuU2TtDKdyDHS1BF2GREX6jLdcF9xVrH6C+gRwwjn3QOX51/GCexHn3EHn3JBzbqivr28taxRZU+PpPH3r1O0hPmVGIdkMzZ2BlbBiUDvnTgMvmNneyqpbgceqWpVIFY2nC/Ssawq6DImKzKjXmg5oDDX4H/XxHuBLlREfx4HfqV5JItXjnGM8k2dDu1rU4lN6BNqD7SXwFdTOuSPAUHVLEam+VK7IXMnRqxa1+JUZgc4tgZagKxOlrkxkvJGl6voQ39KjgbeoFdRSV8bTeQB61PUhfjgHs2MKapFaGkt7LeoN7WpRiw/5GSgXoW1DoGUoqKWuzHd9KKjFl+yEt2wN9gYTCmqpK6ncHADdbY0BVyKRkJ30lq1qUYvUzHR2joaE0ao7u4gfs2pRi9TcdHaOrtZGLMCLFyRC5lvU6qMWqZ35oBbx5WzXh1rUIjWTys7RqaAWvxTUIrWnoJZVmZ2Apg5IBvtvRkEtdUVdH7IquWlo7Q66CgW11BcvqP3ORSZ1LzcNLV1BV6GglvrhnCOVK6pFLf4pqEVqK1MoUSo7BbX4p6AWqa3prHdVYmeLglp8UlCL1Nb0rBfUalGLbwpqkdqan+dDQS2+lMuQTymoRWrpbNeHglr8KMwATkEtUkvzQa0WtfiSm/aWCmqR2kmpRS2rMR/UzZ3B1oGCWupIKjuHGXQ064IX8UEtapHam87O0dnSSCKhKU7FBwW1SO1NZ+fo1OXj4peCWqT2dPm4rIqCWqT2NHOerIq+TBSpvVSlj1rEl9x0ZS7q4LvLFNRSN1I5BbWsQi4cVyWCglrqSCpbpKtNQS0+5aZCE9S+2vRm9iwwA5SAonNuqJpFiay1QrFMdq5EZ0vwH2MlIkIyIRP4DOqKVznnxqpWiUgVzU/IpKsSxbfcFHRuCboKQF0fUidSmotaVis7Hfjdx+f5DWoHfN/MDpvZgeV2MLMDZnbIzA6Njo6uXYUia0ATMsmq5aagpTvoKgD/QX2zc+5a4Hbg3Wb2iqU7OOcOOueGnHNDfX19a1qkyMVK5YoAujJR/CkVvbmoQ3AHcvAZ1M65k5XlCPC3wHXVLEpkranrQ1bl7FWJ3YGWMW/FoDazdjPrmH8M/BrwSLULE1lL6vqQVclNecuQtKj9fA7sB/7WzOb3/7Jz7rtVrUpkjWnUh6xKdspbhqRFvWJQO+eOA/trUItI1aSyRZqSCZobNNBJfMhNesuIjfoQiTRvitNGKp8MRS5svkUdkq4PBbXUhVROc1HLKmQrLeqQdH0oqKUuaOY8WZWQfZmooJa6kKp0fYj4kp2ChlZoaA66EkBBLXVCd3eRVclNheaLRFBQS53wuj7URy0+ZadC0+0BCmqpA865ypeJalGLT9mp0HyRCApqqQOZQom5kqNbQS1+5abUohappYl0AYAN7U0BVyKRoRa1SG1NzCqoZZXUohaprYlMHlBQi0/FAhTSGvUhUkvj6vqQ1Zit3HGwvTfYOhZQUEvsTarrQ1YjU7lDVXt4boCioJbYG88UaEwa65o1jlp8UFCL1N5kpsCG9ibNnCf+ZOa7PhTUIjUzli7Q0x6OORskAtIj3lJBLVI7p6dzbOpqCboMiYrMKCSbobkj6ErOUlBL7J1J5ejvVFCLT+kRWLcRQtRVpqCWWMsXS4xnCmxSUItfqZPQuSXoKhZRUEusjaS8i102damPWnyafgG6tgZdxSIKaom1M6kcgLo+xJ9yGaZPKqhFaunUtBfUA12tAVcikZAZgfKcglqklo6PpjGDHT1tQZciUTD1grdUUIvUzrHRDFu6W2lpTAZdikTB2BPesmd3sHUsoaCWWDs+mmZX37qgy5CoGDnqjaHesDPoShZRUEtslcqO46MZdvW2B12KRMXo49C3BxLh+gSmoJbYeuxUiuxciWu2dwddikRBuQynjkD/FUFX8iIKaomtB54ZB+D6nT0BVyKRMPKoNxf1zlcEXcmL+A5qM0ua2c/N7J5qFiSyVr7/6Bl29bZrng/x5/Fve8tdrwy0jOWspkX9PuBotQoRWUs/PT7Og89O8LbrtwddikTB7AT87M/h0ldD5+agq3kRXzOpm9lW4PXAx4APVrUiqZrnx2f52bMTTM4WKJYdzoHDAeDc4n2dW7zeLdo2v+78+8yvdOf5meV+7kX7LHjBc/usXNdsocjdD59isKeNf/2ybcgqOAdnHoEzj0J2Eoq5xdvOPVlm/dK/FLfMNudv29nnq9m2TD1+3t+V4dgPIDcNr/pDwsjvLS/+FPgQcN55/8zsAHAAYPt2tWLC5nP/9Ax/9A+PUXYr77tW5icfs7PPbdHzxfss3vlC+5zvdefXNSSNf7G3j4/cfhkdLY0X+2vUj5OH4e/eDaO1/OA8/5dp539+oW1nn/vZdoH37N0Nd34Stlx7Mb9M1awY1GZ2BzDinDtsZq88337OuYPAQYChoaEaxoGs5Ohwio99+yi37NvIh2/fx8bOFhoTXq/XwpkcVxuIi/YJ0ZSQ8is4cRi+cAe09cAb7oIdN3mPG5deer/MP5iF6/2EqP6trJqfFvVNwJ1m9jqgBeg0s79yzv1mdUuTtfLFnzxHS0OCP37zfrrbdINXWSKfhq+907ujybvuhY7+oCuSJVb8MtE59xHn3Fbn3CDwFuA+hXR0lMqOex87wyv3bVRIy/Lu/xhMPw9vPKiQDimNo465p0ZmGEvnuXXfxqBLkTA6cQh++hl42e/C9huCrkbOw++XiQA4534I/LAqlUhVHB1OAXDllq6AK5HQKebhW++BjgG49b8EXY1cwKqCWqLn6PAMTQ0Jdmq+C1nqOx+CkcfgrX8DLZ1BVyMXoK6PmDs6nGJP/zoakvqrlgWO3QeHPw83fxD23hZ0NbICnb0x9/zELDt7Nc2nLDF+zFve8B+CrUN8UVDHWLnsGJ7OsVlzXchSee+7C5rPew2bhIiCOsbGMwUKxTKbu3W/QFkin4ZEAzTo7uxRoKCOseHpLAADalHLUoW015rWVYKRoKCOsVNT3oQ6alHLi+RnoEndHlGhoI4xtajlvPIz0KwvmaNCQR1jE5kCCYP1unRclprv+pBIUFDH2HimwPq2JhIJ9UPKEvkZaFKLOioU1DE2mSmwoV2taVlGPq2ujwhRUMfYeKbAegW1LCc/o66PCFFQx9hkpkCPglqWMzuuUR8RoqCOsQm1qGWph/8aPnszlPKwfkfQ1YhPmj0vpsplx+SsWtSyxI8/DoVZuP1/wEt/J+hqxCcFdUylcnOUnYbmyRL5Gbj0Vrj+3wVdiayCuj5iKpUtAtDZqrtwywL5NDRr7umoUVDHVDrvBfW65mTAlUholMtQ0GiPKFJQx9S5oFaLWioKaW+pC10iR0EdU5lKULerRS3z5oNaLerIUVDH1LkWtb4vlor8jLdUUEeOgjqmzgZ1i4JaKs4Gtb5MjBoFdUyd6/pQUEvF2aBWH3XUKKhjar5F3d6koJYKdX1Els7imErnirQ1JUlqilPJTcPU8zDymPdcQR05CuqYyhSK6vYQz+fvgNO/8B5bElrXB1uPrJrO5JhK50sa8SGe6Rdg92vhZe+Czs3Q0hV0RbJKOpNjKp2bU1ALOAe5FGy6Eva8Nuhq5Fe04peJZtZiZg+a2cNm9qiZfbQWhcnFyeRLuthFoJABV1IrOuL8jPrIA7c45/YDVwO3mdkNVa1KLlo6X1SLWrwvEgFaNHY6ylY8k51zDqhce0pj5Y+rZlFy8RTUAiwIarWoo8zXOGozS5rZEWAEuNc598Ay+xwws0Nmdmh0dHSNy5TVyuQ16kNQUMeEr6B2zpWcc1cDW4HrzOyKZfY56Jwbcs4N9fX1rXGZslpqUQsA+ZS3VFBH2qquTHTOTQH3A7dVpRpZE3OlMvliWS1qOdeiblZQR5mfUR99ZtZdedwKvAZ4vMp1yUXIaOY8KZfhhQfh5GHvuVrUkebnTB4AvmBmSbxg/6pz7p7qliUXQ1OcCsfvh796o/e4oVVBHXF+Rn38ArimBrXIGklr5jxJnfKWb/4CDFwFDbrJcZTpTI6hjOailuykt7zkFo2hjgFNcxpD6XwJ0I1t61puypuASTPlxYKCOobSOXV91L3spDdLnmma2zhQUMeQRn3I2aCWWFBQx5BGfQjZSWjbEHQVskYU1DGkUR+iFnW86EyOoUy+SHNDgsak/h+uK+kR+PHHoZiHiWdg4+VBVyRrREEdQ5rno049+T144LPQ1gONrTB4c9AVyRrR2RxDac2cV5+yE97yfQ9rWF7M6LNxDGXUoq5P2UlINELTuqArkTWmoI4hdX3UqdkJjZ2OKQV1DHldH7oqse5kJzQkL6YU1DHk3dhWLeq6k53SkLyYUlDHUDpfpEMTMtWf2QloVYs6jnQ2x1A6V6S9SX+1dWF2AmbHK4/HYLNmJI4jnc0xUyo7snPq+qgLpSLcdTXkp8+tW7cxsHKkenQ2x0ym4F0+rq6POpCd9EL6mt+EXa/yRntcckvQVUkV6GyOGU1xWkfmuzx2vQqufFOwtUhV6cvEmMloQqb6MX8lYltPsHVI1SmoY2Z+5rwOBXX8zbeoNXY69hTUMaMpTuvIrFrU9UJBHTPnuj50ZWLszbeoNXY69tTsipn5G9t2NDcGXIlURbkEqZPe46nnoaEFmtqCrUmqTkEdM+ncHKAWdWz9wwfh8OfPPe/eHlgpUjsK6pjJFLwWtfqoY2rsKejdAze9z3ve/5Jg65Ga0NkcM+l8kYaE0dygrx9iKTMGffu8i1ykbuhsjpl0rsi6lgZMcxLH0+wYtPcGXYXUWKha1A/f/zXMlUkkk5glSSQSJJNJrKmdTZfup6NL326vJJPXhEyxVS55Q/LaFNT1ZsUz2sy2AV8E+gEHHHTO3VWNYnb/8N20WX7ZbUWX4IG+X+elBz5LQ1NzNd4+FnR3lxjLTgJOLeo65OeMLgK/75x7yMw6gMNmdq9z7rG1LubUG79JqThHuVSmXC5RLntLl50m//h3uX7sm/z0c63c8O8/vdZvHRvpvNf1ITGUGfOWusCl7qx4RjvnhoHhyuMZMzsKbAHWPKgv3X+B29u/5m08cNfbGRr+Cmde+AD923av9dvHwkyuSM+6pqDLkLUy9QJ850NQzEGuMp2pWtR1Z1VfJprZIHAN8MAy2w6Y2SEzOzQ6OrpG5S227c7/RIOVOX7vwaq8fhzM5OboaNHFLrFx/IfwxLe9qxAt4U1jOrA/6Kqkxnx/RjazdcA3gPc751JLtzvnDgIHAYaGhtyaVbjA5p37eLTxCjae/MdqvHwsqI86ZtJnvOW/+R40tgZbiwTGV4vazBrxQvpLzrlvVrekC0tvfQWXlI4zeuZEkGWEVipXpFN91PGRGYWmDoV0nVsxqM0bkPsXwFHn3J9Uv6QL67zs1QCcOqJW9VL5YolCsay7u8RJekS31xJfLeqbgN8CbjGzI5U/r6tyXee188obmHNJZp97KKgSQmsmN38bLvVRx4aCWvA36uOfgNBc5tbS2s7xhu20jz8SdCmhcy6o1aKOtLnsucfpM7DxsuBqkVCI5Bk92rGP3VP/HHQZoTNTmTlPLeoI+9H/hPv+aPG6S14VTC0SGpEM6tKGPWyY+g6pqVE6u/uCLic01KKOgZMPQcdmuP5AZYXBFW8MtCQJXiTP6NaBvXAchp/+JZ1DtwRdTmica1FH8q9VAFKnoP9yuPkDQVciIRLJ2fN6d3hz8E6feDzgSsIlVWlRd6rrI7pmTkPHQNBVSMhEMqgHBvdRdAnmRp4MupRQUddHxJWKkBlRUMuLRPKMbmhq4WSyn6bpZ4IuJVTSlaDWlYkRUS7D+NPgvLvykBkDV4ZOBbUsFtkzeqJlO+uzzwVdRqjM5OZoa0rSkIzkB6X68+BB+O5/fPF63QdRlohsUGc7d3JJ5gilUolkUjdyBa/rQ90eETL6OLR0wRsWTO/e2Aa7NBxPFovsWW09u2k7nWf45DMMbL806HJCYSavmfMiJXUSunfAS34j6Eok5CL7Gbl9kxfOY89r5Me8VFYt6kiZPgldW4OuQiIgskHdt30fAOnTTwVcSXhMZAr0tOumAZGROgGdW4KuQiIgss2vni2XMOeSlMY08mPeRKbASzZ3Bl2GLHX0bvjxx8EtnKbdeXds6VJQy8oiG9SJhkZOJzfSNKORHwDOOSZmC2zQbbjC55dfh7GnYceNi9fv2wZ7Xx9MTRIpkQ1qgMnmLXRldQMBgEzBm4t6Q5uCOnQmn4Vt18Hbvxp0JRJRke2jBsh3bKe/NEypXJU7f0XKZKYAwAb1UYfP1HOwfkfQVUiERbpFbRt20T2S4eSZYbYMbA66nECNK6iDVyp6XxAuVJiF7CSsHwykJImHSAd1a/+l8DiMPvd43Qe1WtQhcPd74ciXlt+24ZLa1iKxEumg7tm+F4CZ4SeB+p7udCydBxTUgTp5GDZfC9f928XrG1pgz2uDqUliIdpBvXUfZWeURzWW+vR0DoD+zpaAK6lTpSKMH4OXvxuuflvQ1UjMRDqoE81tDCc30jJ9LOhSAndqOkdPexMtjZr3pKpKRTjxIJTmFq9Pj0B5Dnr3BFOXxFqkgxpgrGWQnuyzQZcRuFNTWTZ3twZdRvw9/GX41nvOv33TlbWrRepG5IN6tnMXuzMPUSrOkWyo3wmJhqezDPa0B11G/J04BC3d8JYvv3hbcwcMXFXzkiT+Ih/UiU1X0XL6Kzz31C/YcdlLgy4nMMNTOW68pDfoMqKvWIBC+vzbh4/AwH4YvKlmJYlEPqj79rwMjsDIkw/WbVBPZArM5ItsXa+uj4viHHzm5d5dVy7kxvfWph6RisgH9bY9V5N1TZROPBR0KYE5OpwCYN8mTch0USaOeyG9/62w+Zrl97EEXP7rNS1LJPJBnWxo5Onmy+mbOBR0KYE5G9QDHQFXEkLpUXji2969CFcy/LC3vPG90H95desSWYXIBzVAZvONXPnspxk7c4re/vq7QvGxUyn6OprpXdccdCnhc//H4PBf+t+/cwv07atePSK/glgEde81r4dnP83xH32Z3jf/QdDl1FSp7PjRU2Ncv3ND0KVU1zM/8uZvXq0nvwt7boc7PuFv/5YuSER6rjKJoRWD2sw+B9wBjDjnrqh+Sat3yZU38szf76D36P+lVHw/yYZY/P/jy+HnJhlL5/m1l/QHXUr1PPNj+MIbfvWff/VHoXNg7eoRqTE/ifZ54FPAF6tbyq/OEgnGr30PQ4f+gJ/8+e9x7Tv/mObWdUGXVXUTmQL/7TtH6Wpt5NbLFgR1sQDTL9S2mJGj8IP/CnPZtX/t3DS0roff+jtIrPLKy2STrhaUyDPnVp7L2cwGgXv8tqiHhobcoUO1/XLPlcv87FO/zXUTdzPnkozbBkqWpESSshlgq3/NNZzm+sIv5ZZ9uNLPlsoOAzZ1tSy++3j6jDe1Zq1t2AXbbqjOa++9HS6/szqvLRICZnbYOTe03LY16yMwswPAAYDt27ev1cv6f/9Egpf93hf55T/fzczR+2jInMbKRRKUMFfGi7sLh/WyW22F7cv+yIV/aKXXMZ9v1NyQoL+r9cV3Ht90lXdBRkMNx1VbAi69Fdpi3lcuEoDYtKhFRKLsQi1qfb0tIhJyCmoRkZBbMajN7CvAT4C9ZnbCzN5V/bJERGTeil8mOufeWotCRERkeer6EBEJOQW1iEjIKahFREJOQS0iEnK+LnhZ9YuajQLPrfkLL68XGKvRe0WJjsvydFyWp+OyvFoelx3Oub7lNlQlqGvJzA6d72qeeqbjsjwdl+XpuCwvLMdFXR8iIiGnoBYRCbk4BPXBoAsIKR2X5em4LE/HZXmhOC6R76MWEYm7OLSoRURiTUEtIhJykQ1qM7vNzJ4ws6fN7MNB11NrZvY5Mxsxs0cWrNtgZvea2VOV5frKejOzP6scq1+Y2bXBVV49ZrbNzO43s8fM7FEze19lfb0flxYze9DMHq4cl49W1u80swcqv//fmFlTZX1z5fnTle2Dgf4CVWZmSTP7uZndU3keuuMSyaA2syTwv4DbgcuBt5rZ5cFWVXOfB25bsu7DwA+cc7uBH1Seg3ecdlf+HAA+U6Maa60I/L5z7nLgBuDdlX8X9X5c8sAtzrn9wNXAbWZ2A/DfgU845y4FJoH5KYzfBUxW1n+isl+cvQ84uuB5+I6Lcy5yf4CXA99b8PwjwEeCriuA4zAIPLLg+RPAQOXxAPBE5fH/Bt663H5x/gP8PfAaHZdFx6QNeAi4Hu+Ku4bK+rPnFPA94OWVxw2V/Szo2qt0PLbi/ed9C3AP3i1NQ3dcItmiBrYALyx4fqKyrt71O+eGK49PA/2Vx3V3vCofS68BHkDHZf7j/RFgBLgXOAZMOeeKlV0W/u5nj0tl+zTQU9OCa+dPgQ8B5crzHkJ4XKIa1LIC5/23X5djL81sHfAN4P3OudTCbfV6XJxzJefc1XgtyOuAfcFWFDwzuwMYcc4dDrqWlUQ1qE8C2xY831pZV+/OmNkAQGU5UllfN8fLzBrxQvpLzrlvVlbX/XGZ55ybAu7H+0jfbWbzd3la+LufPS6V7V3AeG0rrYmbgDvN7Fngr/G6P+4ihMclqkH9M2B35dvZJuAtwLcCrikMvgW8o/L4HXh9tPPrf7syyuEGYHpBV0BsmJkBfwEcdc79yYJN9X5c+sysu/K4Fa/f/iheYL+pstvS4zJ/vN4E3Ff5JBIrzrmPOOe2OucG8TLkPufc2wnjcQm6M/8ivgR4HfAkXl/bHwZdTwC//1eAYWAOrx/tXXj9ZT8AngL+EdhQ2dfwRskcA34JDAVdf5WOyc143Rq/AI5U/rxOx4WrgJ9XjssjwH+urN8FPAg8DXwNaK6sb6k8f7qyfVfQv0MNjtErgXvCelx0CbmISMhFtetDRKRuKKhFREJOQS0iEnIKahGRkFNQi4iEnIJaRCTkFNQiIiH3/wHmO0sLMsHzMwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# use of the `reverse` keyword reverses the slab representation of\n", "# this specific component\n", "si = SLD(2.07, 'Si')\n", "poly1 = SLD(1.0, 'poly1')\n", "poly2 = SLD(4.0, 'poly2')\n", "d2o = SLD(6.36, 'd2o')\n", "\n", "# decay from one to the other\n", "poly1_l = poly1(100., 3.)\n", "poly2_l = poly2(100., 3.)\n", "\n", "ff = FunctionalForm(200., 30., 3., poly1_l, poly2_l, reverse=True)\n", "\n", "struct = si | poly1_l| poly2_l | d2o(0, 3)\n", "struct2 = si | poly1_l | ff | poly2_l | d2o(0, 3)\n", "\n", "struct2.contract = 1.\n", "\n", "plt.plot(*struct.sld_profile())\n", "plt.plot(*struct2.sld_profile());" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAiM0lEQVR4nO3deZCjd33n8fdPrau71df09Iw9l8cHDjY2PmiDjU0AZwlgHBN2vbCUOXaXLYeCZO0Am8K7lWwOkqrsbjhCEYIrkATCsgQcF8YsZh2OMgSwmfHF2OPxOeM5PFff3VJLeqTf/vE80nT3qGekaUm/R3o+r6quR3okS995/Mx3vv07jbUWEREJr5jrAERE5NSUqEVEQk6JWkQk5JSoRURCTolaRCTk4q340PXr19vt27e34qNFRLrSzp07j1trx2q91pJEvX37dnbs2NGKjxYR6UrGmH2rvaamDxGRkFOiFhEJOSVqEZGQU6IWEQk5JWoRkZBTohYRCTklahGRkGvJOGqRTpL3SvzwqWPsnVig4JUplS2V5X9/5axB3vbKsx1HKKFiLbz4c3jpMcjP+s9tGbCQ6IPrbm/6VypRS6RNLhR4950/Z8+RuZNeMwZuumyTErWcYC3c/Vvw+Ndrv545S4lapNk+/c9P8/zxeT5/y5W87sIxehM9xAwYY1yHJmH0xN1+kr7uI/CaD0LfKJiY/696C+8ZJWqJrPm8xzd2HODtl2/mrZeqapY67PgSDG+D638fYu3r4lNnokTWoy9OkyuW+I3LNrkORTpBIQv7fgqX3NzWJA1K1BJhO/dNYQxcsW3YdSjSCV56DGwJtr667V+tRC2R9ej+KS7cMMBgOuE6FOkEB3f6x01Xtv2rlaglsvZNZLlgY8Z1GNIpJp7xOw8HNrb9q5WoJZLKZcuB6RxbhntdhyKdYvpFvyPRASVqiaTj83kKXpktI0rUUqfp/UrUIu20fyoHwJaRPseRSEewFmaUqEXa6uB0JVGropY6LBwDbxGGlKhF2ubo7CIAGwbTjiORjjB/xD8OnOXk6+tK1MaYYWPMN40xTxljdhtjrml1YCKtNJUt0BMzDKY1OVfqkJ3wj33rnHx9vXfpZ4D7rLU3G2OSgBr2pKNNZYuM9CW0pofUp5qoR518/WkTtTFmCPhV4N8DWGsLQKG1YYm01tRCgeG+pOswpFNkJ/2jo0RdT9PHucAx4G+NMY8YY/7GGNO/8k3GmFuNMTuMMTuOHTvW9EBFmmkqW2CdErXUq5Koe0ecfH09iToOXAl83lp7BbAAfHzlm6y1d1prx62142NjY00OU6S5prNFhvs0dVzqlJ2A1BD0uLln6knUB4AD1toHg+ffxE/cIh1rcqHAiCpqqVd2wllHItSRqK21h4H9xphfCU79GvBkS6MSaSFrLdPZIiP9StRSJ8eJut5RH78DfDUY8fE88B9aF5JIa2ULJQqlspo+pH75WUgPOfv6uhK1tfZRYLy1oYi0x9yiB8CAxlBLvfJzMLTF2ddrZqJEznzeT9SZlBK11Ck/D6kBZ1+vRC2RsxAk6v6kErXUKT8HSSVqkbapJmpV1FKPchkKqqhF2kpNH9KQ4gJgIeVuNyAlaomchUKlou5xHIl0hPy8f1RFLdI+8/kSoIpa6pSf849qoxZpH7VRS0MqiVoVtUj7LOQ9jIG+pJo+pA4FJWqRtpvPe/Qn41qLWupTrajVmSjSNgt5Tx2JUr/Cgn9MKlGLtM1CvqT2aalfMesfE+42tlKilshZKHga8SH1K/o71pNwt2O9ErVETrZQojehpg+pU6WiTp60sVXbKFFL5OSLJdJK1FKvYg5icWe7u4AStURQrqiKWhpQzDltnwYlaomgxWKZXo2hlnoVFpy2T4MStURQrlgindCtL3Uq5pSoRdptUW3U0ohiVk0fIu2mRC0NUUUt0l5eqUyxZNWZKPVTZ6JIey16ZQAlaqlfMauKWqSdcgV/LWp1Jkrd1PQh0l6LxUqiVkUtdVLTh0h7KVFLwzTqQ6S9Fotqo5YGqelDpL1yQUWtmYlSF2vBy0E87TQMJWqJlBNNH7r1pQ5lD2wZEkrUIm2TUxu1NMJb9I+qqEXap1JRq41a6uLl/aMStUj75IPOxJQStdSjsrtLPOU0DCVqiZS851fUqbhufalDSCrqujaOM8bsBeaAEuBZa8dbGZRIq+SDKeRJJWqpR0jaqBvZ4fON1trjLYtEpA0KpSBR9yhRSx1CUlHrbpVIKXhK1NKAakXdGW3UFvh/xpidxphbWxmQSCvlvTLJnhixmHEdinSCDmv6uM5ae9AYswG43xjzlLX2gaVvCBL4rQDbtm1rcpgizVHwymqflvpVmz46oKK21h4MjkeBu4FX13jPndbacWvt+NjYWHOjFGkSJWppiBcMzwv7Wh/GmH5jzEDlMfDrwK5WBybSCoWg6UOkLiGpqOtp+tgI3G2Mqbz/f1tr72tpVCItUiipopYGdEobtbX2eeCyNsQi0nJ5r6TJLlK/kFTUumMlUtRGLQ0JSUWtO1YiJa9ELY3QhBeR9lNnojTEW4SeFBi34+51x0qkqKKWhhQXnbdPgxK1REzBK5OKa4lTqVMpr0Qt0m6FUlmjPqR+Xt5v+nBMd6xEikZ9SEO8PMSTrqNQopZoUWeiNKSkilqk7fJeSRW11M8rqKIWaTe/M1G3vdRJFbVI+2mtD2mIV9CoD5F2KpctxZJVopb6lfLQo6YPkbap7peoRC31UkUt0l557ZcojVJFLdJelY1tUwnNTJQ6eZqZKNJWlaaPlCpqqVepoIpapJ0qFbXaqKVuqqhF2kuJWhpWKmgctUg75b0SoM5EaYDW+hBprxOdibrtpQ7lMpSLqqhF2qmg4XnSiFLBP6qiFmmfvCa8SCNKwX6JqqhF2kedidIQr1JRK1GLtE1lZqK24pK6VCtqNX2ItE21M1EVtdTDCxK1KmqR9lHThzSk0pmoilqkfQoaRy2NUEUt0n5a5lQaUq2olahF2iZfVBu1NKBaUavpQ6RtCqUyMQNxNX1IPTSOWqT9Cp72S5QGeJqZKNJ2ea+sjkSpXydW1MaYHmPMI8aYe1sZkEir+DuQa7KL1KlDZybeBuxuVSAirZYvltWRKPUL0czEeD1vMsZsAd4G/CnwkZZGtEZFz+OXP72P7IFdlPILGMoYWwYs2FX+I7PsQK1n1bPmdO9YcX7Z+1f5L87wM/uScbav7yOTrOt/Y+NGL4BX/GZrPtuBQkmJWhoQonHU9f4N/zTwe8DAam8wxtwK3Aqwbdu2NQd2JmamJjn4uRu40lPh3zT3nwOxNjcXnPNaePvnmv6xBa+kzkSpX4hmJp42URtjbgSOWmt3GmPesNr7rLV3AncCjI+Pr1a7ttSTX7mdq4p7eOTyP+T81/5rBoZGMLE4mJj/s0Z2SUlu7dLzSx4veWH5+dN/zrLvquP9+yYW+Ng3HmMqW+T+j7y+uZW1LcPPPgvH9jTvM+sx+xI88g/+/69E3xo+yMCV74ONF1fPaNSHNKTDKuprgZuMMTcAaWDQGPMP1tr3tDa0xkwdPcirJu7lkbG3c9U7frcl32FWeezKxX19/OE7ruSdX/gZ9+2e5OZXbWnuF/zqf2nu59UjPw9f/HV48ltr/5zZg/Cur5w4pVEf0ogQzUw8baK21t4B3AEQVNQfC1uSBtj9wF281pQYff2trkNpq6u2j7BlpJd7Hz/U/ETtQioDH/rp2j/n27fD41+H73y0eup9k4fpiRn4zl21/5ueFFz3u5AZW/v3S+fz8hCLQ8z9P+4t6oVqv/S+H3CMdZx3yTWuQ2krYwyvv3CMex49RLlsicXCUOuHwOW3wNP3wRN3V09ds1gkZgw8scptn52Aoc1wzYfbFKSEWkh2IIcGE7W19kfAj1oSyRqdNf8k+zKvZMxEL1FdunmIrz74Ivsms5y7vt91OOGw9Sr46FPLTr3r0w9wzmgfX3jveO3/5lOXwv6HlKjFF5IdyKFLKurJ40fZZI+wf8O/dR2KE5dsHgJg18EZJepT8DsTTzGCZetVsPte+OyrTv1B/RvgPXdBci2dnRJ6pXxnVtRhdfiZnawD+rZd4ToUJy7cOEDMwNNH5lyHEmqn7Uy8+kP+aBNbXv09C8fghQfg8C9h22uaH6SEh1dQRd1M2UP+r7jrzrnUcSRuJOMxzh7q5eBUznUooZb3yqQSp0jUW8Zhy9+c+kOmX4RPXwpHdilRdztvURV1MxUnXqBoe9i45TzXoTizeaSXA0rUp1TwSmsfnje0FVKDfidlZfjWUrE4XHoz9I6s7XvEvVIhFGOooUsSdWJuP0dj69mcSLgOxZktI738/LkJ12GEWlOmkBsD218He74De39c+z3FLFx729q+R9zz8qGYlQhdkqgzuYNMpTaz2XUgDm0Z6ePw7EHNvluFtbZ51+ZdX4H8bO3XPvcaOPb02r9D3FNF3Vxj3hGeG7jOdRhObRnupWzhyOwiW9dpNMJKXtlStk3a2DbWs3rTxvoLYeKZtX+HuOflIdHrOgqgCxK1V8gzyjR7Mptch+LU2ID/L/+x+bwSdQ0FL9gv8VSdic0wegHsugt+/Bcnv9Y3Cle+f/kSjBJepTz0DruOAuiCRD159CAbgJ7Bs1yH4tRoxm9Lm5iv0cEl1UTd8rU+tl8HO/8Ovv/HtV/fejVseHlrY5Dm8NT00TQzx/azAUgNn+06FKfWZ/wbamI+7ziScCqUgkTd6h1eLr0ZLrqJkxY/P/QIfOnNMPm8EnWn8BYhnnYdBdAFiXrh+EEAMuu7YEGiNVjX71fUx5Woa6pW1O3oaK01SWL9hf5xam/rv1+aI0RrfXT88ID81CEAhjZEO1GnEz0MpOIcV9NHTXmvBOBuh5feEX/8tRJ15/AW1fTRLHbuMAAjY1EenOcbzSSZWFCiriXfzoq6FmNg5Bx4+Mvw1Ir9oV/3UbjqA27iktWpjbp5YrkJpskwnAzHBXVpfSalNupVtLXpYzVvuAP2/N/l5/bcB09/T4k6jFRRN08iP8mMGWLYdSAhMNKfZP9k1nUYoVQdnudyh5eXv83/Wepr74apfW7ikdWVy1Auqo26WVKFaRZ6hlyHEQqD6QRzi57rMELpxKiPkN3yg5th5oDrKGSlUnj2S4QuSNR93jS5xLDrMEJhsDfObK7oOoxQyheDirrVw/MaNbQF8jOwuMqUdHGjurGthuc1RaY8w6HkK1yHEQqD6QRzeY9S2fp7A0pVaCvqoWC00j2/DakB/7GJ+Wtjb7jIXVxRV03UWpRp7axlyM5SSmtJSYChXn/1wLnFIsN94bjBwiIUnYm1bLkKRl8GB3acODd7yB/K9+Y/dRdX1JVUUTdNKTdDghLl3lHXoYTCYJCoZ3OeEvUKzsdRr2bkHPidHcvP/eUVfrIWdyoVtToT125h6ggApk+JGmAw7f+7O7uoduqVKuOoQ5eoaxnYBHMvuY4i2rxF/6jOxLWbmz4OQGJAiRqWVtRK1CudWD0vZJ2JtQyerYraNS+YOKZEvXa5mWMApDJK1OB3JgLMKFGfpLMq6rNh7jBYe/r3SmuErKLu6Dbq/NwkAOmh9Y4jCYehvqCiVtPHSfLFEsZAvBNGwwxu8juzPnulPwIE4JJ/A2/8r27jipKS2qibprjg7xHYPzzmOJJwqLZR5zTpZaW85++XaDph0f6Xvw0uezdsugLOvszv2HryW66jihYvXBNeOrqiLi9MATA4rIoaoD8ZJ2ZUUdfiJ+oOaJ8GGN4G7/jrE8+/8zHY9U138URRyBJ1R1fU5dwUCzbFYKbfdSihEIsZBtIJtVHXkPdKndE+XUtmA+SmTiQPab2QzUzs0DvX17M4zZzJaBbeEppGXlveK7d+v8RWyWzwjwvH3MYRJdU26nDMR+jQO9fXU5hhPjbgOoxQGepNMKuFmU6S98qt3y+xVTIb/eP8EbdxREnIKuqObqNOFmfJ9Qy6DiNUBtX0UVO+2EFt1CtVKupnf+A3gYA/hG+j1rhpmZC1UXd0ou71ZplPRnsLrpUyqTgT81qTeqW8V+rcpo+hbf4wvR9+4sS5WAI+vg+S6p9pCY2jbp6+8hzFpNaiXmognWBOoz5OUhme15EyY/DhX0DOnzfAM/fDA/8D5o/CunPdxtatSsHMxJC0UZ82URtj0sADQCp4/zettf+91YHVY8DOUU5p5bylBtJx5vJqo14p75UZDqbYd6T1F5x4nJv2jwvHlahbxVv0J7uEZNx9PRV1HrjeWjtvjEkAPzHGfNda+/MWx3ZKpXyWNEVsryrqpQbScebzHuWyJabRMFUFrxy+JU7PVH+wZEL2uNs4upmXD01HItQx6sP65oOnieDH+SIEs8GCTLG+dY4jCZdMKo61kC2WXIcSKh09jnql/mAm7oISdct4+dBsGgB1Ds8zxvQYYx4FjgL3W2sfrPGeW40xO4wxO44da/14z/mpowDE+5WolxpIn9g8QE7o6FEfK/UFM3E1rrp1Oq2iBrDWlqy1lwNbgFcbYy6p8Z47rbXj1trxsbHWr72RnfHX+UgOKFEvlQnW+5jXWOplOnrCy0rJPkj0+xV1uez/SHOV8qHpSIQGJ7xYa6eBHwJvaUk0DcjP+b/29Q5qnY+lBqqbByhRL9VVTR/gj63++efgj0f8n/vucB1RdwlZRV3PqI8xoGitnTbG9AJvAv685ZGdRnHeH6rUryVOlxlIBRW1Rn4sk++mzkSAGz8J+3/hP37sa3DwYbfxdJuQtVHXM+rjbODvjTE9+BX4P1pr721tWKfnBUucZkY2OI4kXNRGfTJrLYVOWj2vHudf7/8AHH3S/5Hm8RY7q6K21j4OXNGGWBpishMUbA+DgxpHvZTaqE9WKHXQ7i5nom8UshOuo+gupQIkel1HUdWxd24sN8WMGSDWqQvttEiljXpOibqqo7bhOhN9o/4aIGUNyWyayoSXkOjYOzdRmGY+pgWZVupPBolabdRV+WIHbWx7JvrWgS3D4ozrSLqHVwjNOh/QwYk6VZxmQSvnnaQnZsik4mr6WCLv+ZVmqlt/++qrzFScdBtHN/EWlaiboc+bYTEx7DqMUMqk4upMXKLa9NEt46hXqszOPbYbpvdDSf/v18zLhypRd+zqef2lWYrJYddhhFJlvQ/xFbq9jTpzln/8+nv840U3wbu+4i6eblDKh6qNujMTtbUM2jnKaY34qCWTjqszcYkTnYld2ka98RXwzq/4bdQP3QmTL7iOqPN12oSXMMrNT9NrShgtyFSTNrhdLh8sUNW1FbUxcPFN/uN9/wJ7f+I2nm4QsgkvHXnnzkz6e8fFM6OOIwmnAbVRL1OpqLtqZuJq0sMn1quWM1MuQ7kYqqaPjrxz5yf9lfMSA61f/KkTDaQ16mOprm/6WKp3GApzUNL//zPm5fyjJrysTXbaT9S9Q0rUtfijPvQXtaIyPC8aFXWwkUZ+1m0cnaxYSdR9buNYoiPv3MrKef3DStS1DKQT5IolvJKWvwTIFfxE3ZeMQEWdHvaPld3KpXHFYHNoVdRr4837iXpo3VmOIwmn6nofGqIHwGLQmZju1pmJS/UO+0fNUjxzxWAHciXqtbELE5StYWBES5zWovU+lssFibo3ShX14rTLKDpbtaIOT9NHRw7PM4tTzJl+hno6MvyWq6xJrUTtyxX8JqDeKFXU937EfxxLwA3/EzZd7jCoDlNtow7POOqOrKhTi8eZjmmyy2oqa1Kr6cOXK5ZIxmP0RGFX9tEL4PJbYP2F/hogBx6CvT92HVVnUUXdHJnCMWbjGkO9mky16UNjqQFyBS8a1TRATwJ+86/8x+Uy/Mmo2qsbVdTwvKYY8CZYTGtnl9UMqDNxmVyxFJ1EvVQsBqlBJepGaXje2pVLZUbLU5T6N7oOJbQqbdTa4NaXK5ajMTSvlvSQEnWjNDxv7aYmj5A0HmZQQ/NWU22jVqIG/HHUkRiaV4sSdeO8yvA8VdRnbPLwfgBSw5scRxJe6YTfcaY2al+u6EVjaF4tStSNU0W9dgtHngegf8O5jiMJL2OM1qReIleIaBs1KFGfiUobdYiWOe24RF08/hwAQ5tf5jiScNN6HyfkiuUIV9TDsKh1PxpSzEK8118+NiQ6LlGXJ19gwaZZP6amj1MZSCeUqAOLUR31Aaqoz0QxF6pmD+jARJ2c3c+R+FnEunWj0ibRmtQnRLrpo3fEX/b0j9bBJ86C3fe6jij8QpioO27Cy/r8Pib61exxOgPpOIdnF12HEQrZQoQ7E6+4BWzJ3/D2x38BLz0KF93oOqpwy89BMuM6imU6qixdnJ9iq32J7LqLXIcSeto38YTFKLdRD26CN3wcfu33g8kvaq8+rcI8pAZcR7FMRyXq53c9CEDvtisdRxJ+GvXhK5bKFEpl+qLa9LFUekgbCtQjPw8pVdRnbGqPv2nnea+81nEk4ZdJJZhbLGKtdR2KU5VJP5Vp9ZGW1nTyuhTm1fSxFqMHvs+zPecxtGGr61BCbyAdp1iy1f0Co6ryW0UmmK0ZaekhNX3UIz+npo8ztW/3Ti4s7GZy65tch9IRtHmAr/Lnz6RUUWuBpjqFsDMxVHevVyzgeUU8r0ipWKRYLFL2iswc2w/fvo0sac6/4XbXYXaEpSvojQ2EZ9v7dqtU1Gr6wG/6OKpEfUrWhrIz8bR3rzFmK/BlYCNggTuttZ9pRTDFT2ym1xROOr8ByNkkT7/uM1y2QRNd6pFJaWEmOLEmtypqNPmlHl4eyl7oOhPruXs94KPW2oeNMQPATmPM/dbaJ5sdzGPnf9Af89mThFgPsVgc05Mg1jvI9qtu4LKNW5r9lV1rQJsHAEvbqJWoSQ36v9aXy/5a1XKywrx/THZYRW2tfQl4KXg8Z4zZDWwGmp6or37fnzT7IyOrUkHORXyI3pxGfZyQHgJbhm+8H8b/I5z/RtcRhU9+zj+GrKJu6J9VY8x24ArgwZZEI00zGIxyiHpnYrWNOqVRH1x0I1z8dth9D9z9QdfRhFMlUYesM7HuRG2MyQB3Abdba08a42OMudUYs8MYs+PYsWPNjFHOwFCvn5imsye3+UfJ/KJHT8yQTuhXfdadB+/8Mrz2P8PitOtowqnS9BGyzsS67l5jTAI/SX/VWvtPtd5jrb3TWjturR0fGxtrZoxyBgZ74yR6DBML0U7Uc4tFMqk4JkRLVjqXGvB3MSlF+7etmrKT/rF3xG0cK5w2URv/Dv8isNta+8nWhyTNYIxhtD/FxHzedShOzeU9jfhYqfJrfWHObRxhlAsSdd86t3GsUE9FfS3wXuB6Y8yjwc8NLY5LmmA0k+T4fLQr6plssdoMJIFKR1l+3m0cYZSd8I+94UrU9Yz6+Amg3xs70GhGFfXxhQKjmaTrMMKlWlErUZ8kO+kPD072u45kGfWwdLH1qqiZmM+zPhPdmZk1VTrKVFGfLDfpV9Mh69NQou5i6zMpjs/nI72C3sR8gdF+VdTLqI16ddmp0LVPgxJ1VxvtT5L3ypFdlzpb8MgVS4yqol5OFfXqKhV1yChRd7FNw/6+b4emo7kl10TQ7KOKeoVqZ6Iq6pNkJ6AvXEPzQIm6q20Z8RP1gams40jcqIwhV2fiCpV1LNSZuJy1MHsIBje7juQkStRdbHM1UeccR+JGZcTLOlXUy1Uqam0isNzitP+P11D4NiZRou5iY5kUqXgsshX1i5P+n3vLSJ/jSEImnoJ158PPPgsTz7mOJjxmDvjHofCt0qlE3cWMMWwe6WX/ZDQr6r3HF8ik4qxX08fJ3husBHHXfwIv2kM4q6qJWhW1tNn5YxmePhLNTqO9E1m2r+/TOh+1jGyHmz4Lhx6Gn7ZkH5DOU03UaqOWNrt08xDPH1+I5AYCeycWOGc0XDPMQuXit8PLb4Qff/JEkoqyI7sgPQyZja4jOYkSdZe7dPMQAE8eilbHUbbgsX8yy/lj4VpXOHTe/Gf+ZgL3/4HrSNw79Ahsujx0sxJBibrrXbZ1mJiBf3n2uOtQ2urRF6cpW7hi67DrUMJt5By49jbYdRfs/YnraNwp5uDIk7DpCteR1KRE3eXW9Se55vxR7nnsEOVydKaS3/fEYZLxGK8+N3yzzELn2tv9DrRv3x7dzW+f+g6Ui3BeOLcnU6KOgHeOb2XvRJa/+tGzlCKQrJ84NMM3dhzgrZecRb/Woj69ZB+8469h6gX4+5vgxQf9DXCjIjsJD/wvGN4G21/nOpqaTCsW7BkfH7c7duxo+ufKmbHW8qGvPsx3dx2mP9nDcF+SVCJGLIRtcbVYa7EAFsrB47K1WEvwYylbsFhKZTgerJj3rd++ls3BNHqpw5774O7f8id+xHuhbxQS6SVvCO6X6n2z8v4Jckk1p5zp81Z85ornS8/l5/x2+lv+ES74V7hijNlprR2v+ZoSdTSUy5b7njjMQy9MMrtYJO+VV/yFCCeLxWDAQMwYDBAz/hhxY8DgH2PB41gMtq7r413jW7UY05lYnIU934XDj/uVZilYz3zVxFc5t0oSP+Pnp3rP6V4/zfOV51IDcOnNsPEVuKRELSIScqdK1GqjFhEJOSVqEZGQU6IWEQk5JWoRkZBTohYRCTklahGRkFOiFhEJOSVqEZGQa8mEF2PMMWBf0z+4tvVAtJaGq4+uS226LrXputTWzutyjrV2rNYLLUnU7WSM2bHabJ4o03WpTdelNl2X2sJyXdT0ISISckrUIiIh1w2J+k7XAYSUrkttui616brUForr0vFt1CIi3a4bKmoRka6mRC0iEnIdm6iNMW8xxuwxxjxrjPm463jazRjzJWPMUWPMriXn1hlj7jfGPBMcR4Lzxhjzl8G1etwYc6W7yFvHGLPVGPNDY8yTxpgnjDG3Beejfl3SxpiHjDGPBdflj4Lz5xpjHgz+/F83xiSD86ng+bPB69ud/gFazBjTY4x5xBhzb/A8dNelIxO1MaYH+BzwVuBi4N3GmIvdRtV2fwe8ZcW5jwPft9a+DPh+8Bz86/Sy4OdW4PNtirHdPOCj1tqLgauBDwf3RdSvSx643lp7GXA58BZjzNXAnwOfstZeAEwBHwje/wFgKjj/qeB93ew2YPeS5+G7LtbajvsBrgG+t+T5HcAdruNycB22A7uWPN8DnB08PhvYEzz+AvDuWu/r5h/gW8CbdF2WXZM+4GHgNfgz7uLB+erfKeB7wDXB43jwPuM69hZdjy34/3hfD9yLv7Fi6K5LR1bUwGZg/5LnB4JzUbfRWvtS8PgwsDF4HLnrFfxaegXwILoulV/vHwWOAvcDzwHT1loveMvSP3v1ugSvzwCjbQ24fT4N/B5QDp6PEsLr0qmJWk7D+v/sR3LspTEmA9wF3G6tnV36WlSvi7W2ZK29HL+CfDXwcrcRuWeMuRE4aq3d6TqW0+nURH0Q2Lrk+ZbgXNQdMcacDRAcjwbnI3O9jDEJ/CT9VWvtPwWnI39dKqy108AP8X+lHzbGxIOXlv7Zq9cleH0ImGhvpG1xLXCTMWYv8H/wmz8+QwivS6cm6l8ALwt6Z5PAvwPucRxTGNwDvD94/H78NtrK+fcFoxyuBmaWNAV0DWOMAb4I7LbWfnLJS1G/LmPGmOHgcS9+u/1u/IR9c/C2ldelcr1uBn4Q/CbSVay1d1hrt1hrt+PnkB9Ya28hjNfFdWP+GjoBbgCexm9r+2+u43Hw5/8a8BJQxG9H+wB+e9n3gWeAfwbWBe81+KNkngN+CYy7jr9F1+Q6/GaNx4FHg58bdF14JfBIcF12AX8QnD8PeAh4FvgGkArOp4Pnzwavn+f6z9CGa/QG4N6wXhdNIRcRCblObfoQEYkMJWoRkZBTohYRCTklahGRkFOiFhEJOSVqEZGQU6IWEQm5/w/CAMX6P0cmUAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAiM0lEQVR4nO3deZCjd33n8fdPrau71df09Iw9l8cHDjY2PmiDjU0AZwlgHBN2vbCUOXaXLYeCZO0Am8K7lWwOkqrsbjhCEYIrkATCsgQcF8YsZh2OMgSwmfHF2OPxOeM5PFff3VJLeqTf/vE80nT3qGekaUm/R3o+r6quR3okS995/Mx3vv07jbUWEREJr5jrAERE5NSUqEVEQk6JWkQk5JSoRURCTolaRCTk4q340PXr19vt27e34qNFRLrSzp07j1trx2q91pJEvX37dnbs2NGKjxYR6UrGmH2rvaamDxGRkFOiFhEJOSVqEZGQU6IWEQk5JWoRkZBTohYRCTklahGRkGvJOGqRTpL3SvzwqWPsnVig4JUplS2V5X9/5axB3vbKsx1HKKFiLbz4c3jpMcjP+s9tGbCQ6IPrbm/6VypRS6RNLhR4950/Z8+RuZNeMwZuumyTErWcYC3c/Vvw+Ndrv545S4lapNk+/c9P8/zxeT5/y5W87sIxehM9xAwYY1yHJmH0xN1+kr7uI/CaD0LfKJiY/696C+8ZJWqJrPm8xzd2HODtl2/mrZeqapY67PgSDG+D638fYu3r4lNnokTWoy9OkyuW+I3LNrkORTpBIQv7fgqX3NzWJA1K1BJhO/dNYQxcsW3YdSjSCV56DGwJtr667V+tRC2R9ej+KS7cMMBgOuE6FOkEB3f6x01Xtv2rlaglsvZNZLlgY8Z1GNIpJp7xOw8HNrb9q5WoJZLKZcuB6RxbhntdhyKdYvpFvyPRASVqiaTj83kKXpktI0rUUqfp/UrUIu20fyoHwJaRPseRSEewFmaUqEXa6uB0JVGropY6LBwDbxGGlKhF2ubo7CIAGwbTjiORjjB/xD8OnOXk6+tK1MaYYWPMN40xTxljdhtjrml1YCKtNJUt0BMzDKY1OVfqkJ3wj33rnHx9vXfpZ4D7rLU3G2OSgBr2pKNNZYuM9CW0pofUp5qoR518/WkTtTFmCPhV4N8DWGsLQKG1YYm01tRCgeG+pOswpFNkJ/2jo0RdT9PHucAx4G+NMY8YY/7GGNO/8k3GmFuNMTuMMTuOHTvW9EBFmmkqW2CdErXUq5Koe0ecfH09iToOXAl83lp7BbAAfHzlm6y1d1prx62142NjY00OU6S5prNFhvs0dVzqlJ2A1BD0uLln6knUB4AD1toHg+ffxE/cIh1rcqHAiCpqqVd2wllHItSRqK21h4H9xphfCU79GvBkS6MSaSFrLdPZIiP9StRSJ8eJut5RH78DfDUY8fE88B9aF5JIa2ULJQqlspo+pH75WUgPOfv6uhK1tfZRYLy1oYi0x9yiB8CAxlBLvfJzMLTF2ddrZqJEznzeT9SZlBK11Ck/D6kBZ1+vRC2RsxAk6v6kErXUKT8HSSVqkbapJmpV1FKPchkKqqhF2kpNH9KQ4gJgIeVuNyAlaomchUKlou5xHIl0hPy8f1RFLdI+8/kSoIpa6pSf849qoxZpH7VRS0MqiVoVtUj7LOQ9jIG+pJo+pA4FJWqRtpvPe/Qn41qLWupTrajVmSjSNgt5Tx2JUr/Cgn9MKlGLtM1CvqT2aalfMesfE+42tlKilshZKHga8SH1K/o71pNwt2O9ErVETrZQojehpg+pU6WiTp60sVXbKFFL5OSLJdJK1FKvYg5icWe7u4AStURQrqiKWhpQzDltnwYlaomgxWKZXo2hlnoVFpy2T4MStURQrlgindCtL3Uq5pSoRdptUW3U0ohiVk0fIu2mRC0NUUUt0l5eqUyxZNWZKPVTZ6JIey16ZQAlaqlfMauKWqSdcgV/LWp1Jkrd1PQh0l6LxUqiVkUtdVLTh0h7KVFLwzTqQ6S9Fotqo5YGqelDpL1yQUWtmYlSF2vBy0E87TQMJWqJlBNNH7r1pQ5lD2wZEkrUIm2TUxu1NMJb9I+qqEXap1JRq41a6uLl/aMStUj75IPOxJQStdSjsrtLPOU0DCVqiZS851fUqbhufalDSCrqujaOM8bsBeaAEuBZa8dbGZRIq+SDKeRJJWqpR0jaqBvZ4fON1trjLYtEpA0KpSBR9yhRSx1CUlHrbpVIKXhK1NKAakXdGW3UFvh/xpidxphbWxmQSCvlvTLJnhixmHEdinSCDmv6uM5ae9AYswG43xjzlLX2gaVvCBL4rQDbtm1rcpgizVHwymqflvpVmz46oKK21h4MjkeBu4FX13jPndbacWvt+NjYWHOjFGkSJWppiBcMzwv7Wh/GmH5jzEDlMfDrwK5WBybSCoWg6UOkLiGpqOtp+tgI3G2Mqbz/f1tr72tpVCItUiipopYGdEobtbX2eeCyNsQi0nJ5r6TJLlK/kFTUumMlUtRGLQ0JSUWtO1YiJa9ELY3QhBeR9lNnojTEW4SeFBi34+51x0qkqKKWhhQXnbdPgxK1REzBK5OKa4lTqVMpr0Qt0m6FUlmjPqR+Xt5v+nBMd6xEikZ9SEO8PMSTrqNQopZoUWeiNKSkilqk7fJeSRW11M8rqKIWaTe/M1G3vdRJFbVI+2mtD2mIV9CoD5F2KpctxZJVopb6lfLQo6YPkbap7peoRC31UkUt0l557ZcojVJFLdJelY1tUwnNTJQ6eZqZKNJWlaaPlCpqqVepoIpapJ0qFbXaqKVuqqhF2kuJWhpWKmgctUg75b0SoM5EaYDW+hBprxOdibrtpQ7lMpSLqqhF2qmg4XnSiFLBP6qiFmmfvCa8SCNKwX6JqqhF2kedidIQr1JRK1GLtE1lZqK24pK6VCtqNX2ItE21M1EVtdTDCxK1KmqR9lHThzSk0pmoilqkfQoaRy2NUEUt0n5a5lQaUq2olahF2iZfVBu1NKBaUavpQ6RtCqUyMQNxNX1IPTSOWqT9Cp72S5QGeJqZKNJ2ea+sjkSpXydW1MaYHmPMI8aYe1sZkEir+DuQa7KL1KlDZybeBuxuVSAirZYvltWRKPUL0czEeD1vMsZsAd4G/CnwkZZGtEZFz+OXP72P7IFdlPILGMoYWwYs2FX+I7PsQK1n1bPmdO9YcX7Z+1f5L87wM/uScbav7yOTrOt/Y+NGL4BX/GZrPtuBQkmJWhoQonHU9f4N/zTwe8DAam8wxtwK3Aqwbdu2NQd2JmamJjn4uRu40lPh3zT3nwOxNjcXnPNaePvnmv6xBa+kzkSpX4hmJp42URtjbgSOWmt3GmPesNr7rLV3AncCjI+Pr1a7ttSTX7mdq4p7eOTyP+T81/5rBoZGMLE4mJj/s0Z2SUlu7dLzSx4veWH5+dN/zrLvquP9+yYW+Ng3HmMqW+T+j7y+uZW1LcPPPgvH9jTvM+sx+xI88g/+/69E3xo+yMCV74ONF1fPaNSHNKTDKuprgZuMMTcAaWDQGPMP1tr3tDa0xkwdPcirJu7lkbG3c9U7frcl32FWeezKxX19/OE7ruSdX/gZ9+2e5OZXbWnuF/zqf2nu59UjPw9f/HV48ltr/5zZg/Cur5w4pVEf0ogQzUw8baK21t4B3AEQVNQfC1uSBtj9wF281pQYff2trkNpq6u2j7BlpJd7Hz/U/ETtQioDH/rp2j/n27fD41+H73y0eup9k4fpiRn4zl21/5ueFFz3u5AZW/v3S+fz8hCLQ8z9P+4t6oVqv/S+H3CMdZx3yTWuQ2krYwyvv3CMex49RLlsicXCUOuHwOW3wNP3wRN3V09ds1gkZgw8scptn52Aoc1wzYfbFKSEWkh2IIcGE7W19kfAj1oSyRqdNf8k+zKvZMxEL1FdunmIrz74Ivsms5y7vt91OOGw9Sr46FPLTr3r0w9wzmgfX3jveO3/5lOXwv6HlKjFF5IdyKFLKurJ40fZZI+wf8O/dR2KE5dsHgJg18EZJepT8DsTTzGCZetVsPte+OyrTv1B/RvgPXdBci2dnRJ6pXxnVtRhdfiZnawD+rZd4ToUJy7cOEDMwNNH5lyHEmqn7Uy8+kP+aBNbXv09C8fghQfg8C9h22uaH6SEh1dQRd1M2UP+r7jrzrnUcSRuJOMxzh7q5eBUznUooZb3yqQSp0jUW8Zhy9+c+kOmX4RPXwpHdilRdztvURV1MxUnXqBoe9i45TzXoTizeaSXA0rUp1TwSmsfnje0FVKDfidlZfjWUrE4XHoz9I6s7XvEvVIhFGOooUsSdWJuP0dj69mcSLgOxZktI738/LkJ12GEWlOmkBsD218He74De39c+z3FLFx729q+R9zz8qGYlQhdkqgzuYNMpTaz2XUgDm0Z6ePw7EHNvluFtbZ51+ZdX4H8bO3XPvcaOPb02r9D3FNF3Vxj3hGeG7jOdRhObRnupWzhyOwiW9dpNMJKXtlStk3a2DbWs3rTxvoLYeKZtX+HuOflIdHrOgqgCxK1V8gzyjR7Mptch+LU2ID/L/+x+bwSdQ0FL9gv8VSdic0wegHsugt+/Bcnv9Y3Cle+f/kSjBJepTz0DruOAuiCRD159CAbgJ7Bs1yH4tRoxm9Lm5iv0cEl1UTd8rU+tl8HO/8Ovv/HtV/fejVseHlrY5Dm8NT00TQzx/azAUgNn+06FKfWZ/wbamI+7ziScCqUgkTd6h1eLr0ZLrqJkxY/P/QIfOnNMPm8EnWn8BYhnnYdBdAFiXrh+EEAMuu7YEGiNVjX71fUx5Woa6pW1O3oaK01SWL9hf5xam/rv1+aI0RrfXT88ID81CEAhjZEO1GnEz0MpOIcV9NHTXmvBOBuh5feEX/8tRJ15/AW1fTRLHbuMAAjY1EenOcbzSSZWFCiriXfzoq6FmNg5Bx4+Mvw1Ir9oV/3UbjqA27iktWpjbp5YrkJpskwnAzHBXVpfSalNupVtLXpYzVvuAP2/N/l5/bcB09/T4k6jFRRN08iP8mMGWLYdSAhMNKfZP9k1nUYoVQdnudyh5eXv83/Wepr74apfW7ikdWVy1Auqo26WVKFaRZ6hlyHEQqD6QRzi57rMELpxKiPkN3yg5th5oDrKGSlUnj2S4QuSNR93jS5xLDrMEJhsDfObK7oOoxQyheDirrVw/MaNbQF8jOwuMqUdHGjurGthuc1RaY8w6HkK1yHEQqD6QRzeY9S2fp7A0pVaCvqoWC00j2/DakB/7GJ+Wtjb7jIXVxRV03UWpRp7axlyM5SSmtJSYChXn/1wLnFIsN94bjBwiIUnYm1bLkKRl8GB3acODd7yB/K9+Y/dRdX1JVUUTdNKTdDghLl3lHXoYTCYJCoZ3OeEvUKzsdRr2bkHPidHcvP/eUVfrIWdyoVtToT125h6ggApk+JGmAw7f+7O7uoduqVKuOoQ5eoaxnYBHMvuY4i2rxF/6jOxLWbmz4OQGJAiRqWVtRK1CudWD0vZJ2JtQyerYraNS+YOKZEvXa5mWMApDJK1OB3JgLMKFGfpLMq6rNh7jBYe/r3SmuErKLu6Dbq/NwkAOmh9Y4jCYehvqCiVtPHSfLFEsZAvBNGwwxu8juzPnulPwIE4JJ/A2/8r27jipKS2qibprjg7xHYPzzmOJJwqLZR5zTpZaW85++XaDph0f6Xvw0uezdsugLOvszv2HryW66jihYvXBNeOrqiLi9MATA4rIoaoD8ZJ2ZUUdfiJ+oOaJ8GGN4G7/jrE8+/8zHY9U138URRyBJ1R1fU5dwUCzbFYKbfdSihEIsZBtIJtVHXkPdKndE+XUtmA+SmTiQPab2QzUzs0DvX17M4zZzJaBbeEppGXlveK7d+v8RWyWzwjwvH3MYRJdU26nDMR+jQO9fXU5hhPjbgOoxQGepNMKuFmU6S98qt3y+xVTIb/eP8EbdxREnIKuqObqNOFmfJ9Qy6DiNUBtX0UVO+2EFt1CtVKupnf+A3gYA/hG+j1rhpmZC1UXd0ou71ZplPRnsLrpUyqTgT81qTeqW8V+rcpo+hbf4wvR9+4sS5WAI+vg+S6p9pCY2jbp6+8hzFpNaiXmognWBOoz5OUhme15EyY/DhX0DOnzfAM/fDA/8D5o/CunPdxtatSsHMxJC0UZ82URtj0sADQCp4/zettf+91YHVY8DOUU5p5bylBtJx5vJqo14p75UZDqbYd6T1F5x4nJv2jwvHlahbxVv0J7uEZNx9PRV1HrjeWjtvjEkAPzHGfNda+/MWx3ZKpXyWNEVsryrqpQbScebzHuWyJabRMFUFrxy+JU7PVH+wZEL2uNs4upmXD01HItQx6sP65oOnieDH+SIEs8GCTLG+dY4jCZdMKo61kC2WXIcSKh09jnql/mAm7oISdct4+dBsGgB1Ds8zxvQYYx4FjgL3W2sfrPGeW40xO4wxO44da/14z/mpowDE+5WolxpIn9g8QE7o6FEfK/UFM3E1rrp1Oq2iBrDWlqy1lwNbgFcbYy6p8Z47rbXj1trxsbHWr72RnfHX+UgOKFEvlQnW+5jXWOplOnrCy0rJPkj0+xV1uez/SHOV8qHpSIQGJ7xYa6eBHwJvaUk0DcjP+b/29Q5qnY+lBqqbByhRL9VVTR/gj63++efgj0f8n/vucB1RdwlZRV3PqI8xoGitnTbG9AJvAv685ZGdRnHeH6rUryVOlxlIBRW1Rn4sk++mzkSAGz8J+3/hP37sa3DwYbfxdJuQtVHXM+rjbODvjTE9+BX4P1pr721tWKfnBUucZkY2OI4kXNRGfTJrLYVOWj2vHudf7/8AHH3S/5Hm8RY7q6K21j4OXNGGWBpishMUbA+DgxpHvZTaqE9WKHXQ7i5nom8UshOuo+gupQIkel1HUdWxd24sN8WMGSDWqQvttEiljXpOibqqo7bhOhN9o/4aIGUNyWyayoSXkOjYOzdRmGY+pgWZVupPBolabdRV+WIHbWx7JvrWgS3D4ozrSLqHVwjNOh/QwYk6VZxmQSvnnaQnZsik4mr6WCLv+ZVmqlt/++qrzFScdBtHN/EWlaiboc+bYTEx7DqMUMqk4upMXKLa9NEt46hXqszOPbYbpvdDSf/v18zLhypRd+zqef2lWYrJYddhhFJlvQ/xFbq9jTpzln/8+nv840U3wbu+4i6eblDKh6qNujMTtbUM2jnKaY34qCWTjqszcYkTnYld2ka98RXwzq/4bdQP3QmTL7iOqPN12oSXMMrNT9NrShgtyFSTNrhdLh8sUNW1FbUxcPFN/uN9/wJ7f+I2nm4QsgkvHXnnzkz6e8fFM6OOIwmnAbVRL1OpqLtqZuJq0sMn1quWM1MuQ7kYqqaPjrxz5yf9lfMSA61f/KkTDaQ16mOprm/6WKp3GApzUNL//zPm5fyjJrysTXbaT9S9Q0rUtfijPvQXtaIyPC8aFXWwkUZ+1m0cnaxYSdR9buNYoiPv3MrKef3DStS1DKQT5IolvJKWvwTIFfxE3ZeMQEWdHvaPld3KpXHFYHNoVdRr4837iXpo3VmOIwmn6nofGqIHwGLQmZju1pmJS/UO+0fNUjxzxWAHciXqtbELE5StYWBES5zWovU+lssFibo3ShX14rTLKDpbtaIOT9NHRw7PM4tTzJl+hno6MvyWq6xJrUTtyxX8JqDeKFXU937EfxxLwA3/EzZd7jCoDlNtow7POOqOrKhTi8eZjmmyy2oqa1Kr6cOXK5ZIxmP0RGFX9tEL4PJbYP2F/hogBx6CvT92HVVnUUXdHJnCMWbjGkO9mky16UNjqQFyBS8a1TRATwJ+86/8x+Uy/Mmo2qsbVdTwvKYY8CZYTGtnl9UMqDNxmVyxFJ1EvVQsBqlBJepGaXje2pVLZUbLU5T6N7oOJbQqbdTa4NaXK5ajMTSvlvSQEnWjNDxv7aYmj5A0HmZQQ/NWU22jVqIG/HHUkRiaV4sSdeO8yvA8VdRnbPLwfgBSw5scRxJe6YTfcaY2al+u6EVjaF4tStSNU0W9dgtHngegf8O5jiMJL2OM1qReIleIaBs1KFGfiUobdYiWOe24RF08/hwAQ5tf5jiScNN6HyfkiuUIV9TDsKh1PxpSzEK8118+NiQ6LlGXJ19gwaZZP6amj1MZSCeUqAOLUR31Aaqoz0QxF6pmD+jARJ2c3c+R+FnEunWj0ibRmtQnRLrpo3fEX/b0j9bBJ86C3fe6jij8QpioO27Cy/r8Pib61exxOgPpOIdnF12HEQrZQoQ7E6+4BWzJ3/D2x38BLz0KF93oOqpwy89BMuM6imU6qixdnJ9iq32J7LqLXIcSeto38YTFKLdRD26CN3wcfu33g8kvaq8+rcI8pAZcR7FMRyXq53c9CEDvtisdRxJ+GvXhK5bKFEpl+qLa9LFUekgbCtQjPw8pVdRnbGqPv2nnea+81nEk4ZdJJZhbLGKtdR2KU5VJP5Vp9ZGW1nTyuhTm1fSxFqMHvs+zPecxtGGr61BCbyAdp1iy1f0Co6ryW0UmmK0ZaekhNX3UIz+npo8ztW/3Ti4s7GZy65tch9IRtHmAr/Lnz6RUUWuBpjqFsDMxVHevVyzgeUU8r0ipWKRYLFL2iswc2w/fvo0sac6/4XbXYXaEpSvojQ2EZ9v7dqtU1Gr6wG/6OKpEfUrWhrIz8bR3rzFmK/BlYCNggTuttZ9pRTDFT2ym1xROOr8ByNkkT7/uM1y2QRNd6pFJaWEmOLEmtypqNPmlHl4eyl7oOhPruXs94KPW2oeNMQPATmPM/dbaJ5sdzGPnf9Af89mThFgPsVgc05Mg1jvI9qtu4LKNW5r9lV1rQJsHAEvbqJWoSQ36v9aXy/5a1XKywrx/THZYRW2tfQl4KXg8Z4zZDWwGmp6or37fnzT7IyOrUkHORXyI3pxGfZyQHgJbhm+8H8b/I5z/RtcRhU9+zj+GrKJu6J9VY8x24ArgwZZEI00zGIxyiHpnYrWNOqVRH1x0I1z8dth9D9z9QdfRhFMlUYesM7HuRG2MyQB3Abdba08a42OMudUYs8MYs+PYsWPNjFHOwFCvn5imsye3+UfJ/KJHT8yQTuhXfdadB+/8Mrz2P8PitOtowqnS9BGyzsS67l5jTAI/SX/VWvtPtd5jrb3TWjturR0fGxtrZoxyBgZ74yR6DBML0U7Uc4tFMqk4JkRLVjqXGvB3MSlF+7etmrKT/rF3xG0cK5w2URv/Dv8isNta+8nWhyTNYIxhtD/FxHzedShOzeU9jfhYqfJrfWHObRxhlAsSdd86t3GsUE9FfS3wXuB6Y8yjwc8NLY5LmmA0k+T4fLQr6plssdoMJIFKR1l+3m0cYZSd8I+94UrU9Yz6+Amg3xs70GhGFfXxhQKjmaTrMMKlWlErUZ8kO+kPD072u45kGfWwdLH1qqiZmM+zPhPdmZk1VTrKVFGfLDfpV9Mh69NQou5i6zMpjs/nI72C3sR8gdF+VdTLqI16ddmp0LVPgxJ1VxvtT5L3ypFdlzpb8MgVS4yqol5OFfXqKhV1yChRd7FNw/6+b4emo7kl10TQ7KOKeoVqZ6Iq6pNkJ6AvXEPzQIm6q20Z8RP1gams40jcqIwhV2fiCpV1LNSZuJy1MHsIBje7juQkStRdbHM1UeccR+JGZcTLOlXUy1Uqam0isNzitP+P11D4NiZRou5iY5kUqXgsshX1i5P+n3vLSJ/jSEImnoJ158PPPgsTz7mOJjxmDvjHofCt0qlE3cWMMWwe6WX/ZDQr6r3HF8ik4qxX08fJ3husBHHXfwIv2kM4q6qJWhW1tNn5YxmePhLNTqO9E1m2r+/TOh+1jGyHmz4Lhx6Gn7ZkH5DOU03UaqOWNrt08xDPH1+I5AYCeycWOGc0XDPMQuXit8PLb4Qff/JEkoqyI7sgPQyZja4jOYkSdZe7dPMQAE8eilbHUbbgsX8yy/lj4VpXOHTe/Gf+ZgL3/4HrSNw79Ahsujx0sxJBibrrXbZ1mJiBf3n2uOtQ2urRF6cpW7hi67DrUMJt5By49jbYdRfs/YnraNwp5uDIk7DpCteR1KRE3eXW9Se55vxR7nnsEOVydKaS3/fEYZLxGK8+N3yzzELn2tv9DrRv3x7dzW+f+g6Ui3BeOLcnU6KOgHeOb2XvRJa/+tGzlCKQrJ84NMM3dhzgrZecRb/Woj69ZB+8469h6gX4+5vgxQf9DXCjIjsJD/wvGN4G21/nOpqaTCsW7BkfH7c7duxo+ufKmbHW8qGvPsx3dx2mP9nDcF+SVCJGLIRtcbVYa7EAFsrB47K1WEvwYylbsFhKZTgerJj3rd++ls3BNHqpw5774O7f8id+xHuhbxQS6SVvCO6X6n2z8v4Jckk1p5zp81Z85ornS8/l5/x2+lv+ES74V7hijNlprR2v+ZoSdTSUy5b7njjMQy9MMrtYJO+VV/yFCCeLxWDAQMwYDBAz/hhxY8DgH2PB41gMtq7r413jW7UY05lYnIU934XDj/uVZilYz3zVxFc5t0oSP+Pnp3rP6V4/zfOV51IDcOnNsPEVuKRELSIScqdK1GqjFhEJOSVqEZGQU6IWEQk5JWoRkZBTohYRCTklahGRkFOiFhEJOSVqEZGQa8mEF2PMMWBf0z+4tvVAtJaGq4+uS226LrXputTWzutyjrV2rNYLLUnU7WSM2bHabJ4o03WpTdelNl2X2sJyXdT0ISISckrUIiIh1w2J+k7XAYSUrkttui616brUForr0vFt1CIi3a4bKmoRka6mRC0iEnIdm6iNMW8xxuwxxjxrjPm463jazRjzJWPMUWPMriXn1hlj7jfGPBMcR4Lzxhjzl8G1etwYc6W7yFvHGLPVGPNDY8yTxpgnjDG3Beejfl3SxpiHjDGPBdflj4Lz5xpjHgz+/F83xiSD86ng+bPB69ud/gFazBjTY4x5xBhzb/A8dNelIxO1MaYH+BzwVuBi4N3GmIvdRtV2fwe8ZcW5jwPft9a+DPh+8Bz86/Sy4OdW4PNtirHdPOCj1tqLgauBDwf3RdSvSx643lp7GXA58BZjzNXAnwOfstZeAEwBHwje/wFgKjj/qeB93ew2YPeS5+G7LtbajvsBrgG+t+T5HcAdruNycB22A7uWPN8DnB08PhvYEzz+AvDuWu/r5h/gW8CbdF2WXZM+4GHgNfgz7uLB+erfKeB7wDXB43jwPuM69hZdjy34/3hfD9yLv7Fi6K5LR1bUwGZg/5LnB4JzUbfRWvtS8PgwsDF4HLnrFfxaegXwILoulV/vHwWOAvcDzwHT1loveMvSP3v1ugSvzwCjbQ24fT4N/B5QDp6PEsLr0qmJWk7D+v/sR3LspTEmA9wF3G6tnV36WlSvi7W2ZK29HL+CfDXwcrcRuWeMuRE4aq3d6TqW0+nURH0Q2Lrk+ZbgXNQdMcacDRAcjwbnI3O9jDEJ/CT9VWvtPwWnI39dKqy108AP8X+lHzbGxIOXlv7Zq9cleH0ImGhvpG1xLXCTMWYv8H/wmz8+QwivS6cm6l8ALwt6Z5PAvwPucRxTGNwDvD94/H78NtrK+fcFoxyuBmaWNAV0DWOMAb4I7LbWfnLJS1G/LmPGmOHgcS9+u/1u/IR9c/C2ldelcr1uBn4Q/CbSVay1d1hrt1hrt+PnkB9Ya28hjNfFdWP+GjoBbgCexm9r+2+u43Hw5/8a8BJQxG9H+wB+e9n3gWeAfwbWBe81+KNkngN+CYy7jr9F1+Q6/GaNx4FHg58bdF14JfBIcF12AX8QnD8PeAh4FvgGkArOp4Pnzwavn+f6z9CGa/QG4N6wXhdNIRcRCblObfoQEYkMJWoRkZBTohYRCTklahGRkFOiFhEJOSVqEZGQU6IWEQm5/w/CAMX6P0cmUAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# use of the `reverse` keyword reverses the slab representation of\n", "# this specific component\n", "si = SLD(2.07, 'Si')\n", "poly1 = SLD(4.0, 'poly1')\n", "poly2 = SLD(1.5, 'poly2')\n", "d2o = SLD(6.36, 'd2o')\n", "\n", "# decay from one to the other\n", "poly1_l = poly1(100., 3.)\n", "poly2_l = poly2(100., 3.)\n", "\n", "ff = FunctionalForm(200., 30., 3., poly1_l, poly2_l, reverse=True)\n", "\n", "struct = si | poly1_l| poly2_l | d2o(0, 3)\n", "struct2 = si | poly1_l | ff | poly2_l | d2o(0, 3)\n", "\n", "struct2.contract = 1.\n", "\n", "plt.plot(*struct.sld_profile())\n", "plt.plot(*struct2.sld_profile());" ] } ], "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.9.5" } }, "nbformat": 4, "nbformat_minor": 2 }