{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Exact diagonalization of atom in a magnetic field (Breit-Rabi)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import scipy.constants as cnst\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define constants\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "hbar = cnst.hbar\n", "h = cnst.h\n", "bohr_mag = cnst.value('Bohr magneton')\n", "a_hfs = 3.4e9\n", "gJ = 2.0023\n", "gI = 1e-3" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define electron spin operators for $|J=1/2, m_J = [+1/2, -1/2] >$" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "Jz = hbar/2 * np.array([[1,0],[0,-1]])\n", "Jplus = hbar * np.array([[0,np.sqrt((1/2- (-1/2))*(3/2+ (-1/2)))],[0,0]])\n", "Jminus = hbar * np.array([[0,0],[np.sqrt((1/2+ (+1/2))*(3/2- (+1/2))),0]])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define nuclear spin operators $|I=3/2, m_I = [3/2, 1/2, -1/2, -3/2]>$" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "Iz = 3 * hbar / 2 *np.array([[1,0,0,0],[0,1/3,0,0],[0,0,-1/3,0],[0,0,0,-1]])\n", "Iplus = hbar * np.array([[0,np.sqrt((3/2-1/2)*(5/2+1/2)),0,0],[0,0,np.sqrt((3/2-(-1/2))*(5/2+(-1/2))),0],[0,0,0,np.sqrt((3/2-(-3/2))*(5/2+(-3/2)))],[0,0,0,0]])\n", "Iminus = hbar * np.array([[0,0,0,0],[np.sqrt((3/2+3/2)*(5/2-3/2)),0,0,0],[0,np.sqrt((3/2+1/2)*(5/2-1/2)),0,0],[0,0,np.sqrt((3/2+(-1/2))*(5/2-(-1/2))),0]])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define Hamiltonian in basis $|I, m_I>|J, m_J>$" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "def Hamiltonian(B):\n", " H = (a_hfs * h / (2*hbar**2) * (np.kron(Iplus, Jminus)+np.kron(Iminus, Jplus)) + a_hfs * h / hbar**2 * np.kron(Iz, Jz) +\n", " bohr_mag / hbar * B * (gI*np.kron(Iz, np.identity(2))+gJ*np.kron(np.identity(4), Jz)))\n", " return H" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot eigenvalues as a function of magnetic field" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "Bs = np.linspace(0, 2, num=500)\n", "\n", "ws = []\n", "for B in Bs:\n", " w, v = np.linalg.eig(Hamiltonian(B))\n", " ws.append(w)\n", "ws = np.array(ws)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Energy/GHz')" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAERCAYAAABl3+CQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9eXTb6Xnf+3lBggQJcAO47xR3kdRKrZREeTwe2xPHrp3FM0knthN7OoknTm4cu0tOm9TtucltcnOb3jTJnSZp4rQnSZ2kbZKmmXHsEbVrRtJoo7jv+wZuAIj1994/XhCkNCQIUZRIUe/nHBwSP/zww0vMCF887/M830dIKdFoNBqNZj1M270AjUaj0exstFBoNBqNJipaKDQajUYTFS0UGo1Go4mKFgqNRqPRREULhUaj0WiismuFQgjxh0KISSHEvS241gEhxBUhRKsQ4o4Q4vOrHvsDIcTt8PG/EELYHvf1NBqNZichdmsfhRDiDOACvi2lrH/Ma1UBUkrZJYTIB24AtVLKOSFEqpRyIXzebwKTUspfe9z1azQazU5h10YUUsrzgHP1MSFEuRDi74UQN4QQF4QQNTFeq1NK2RX+fRSYBLLC95dFQgBJwO5UXo1G89yya4ViHd4CflZKeRj4ReB3HvUCQoijQALQs+rYfwbGgRrg/92apWo0Gs3OIH67F/C0COcOTgLfUV/+AUgMP/Y54FtrPG1ESvnxVdfIA/4E+IKU0lg+LqX8khAiDiUSnwf+8xP5IzQajWYb2LU5CgAhRCnwt1LKeiFEKtAhpczb5LVSgXPAr0opv7POOc3AN6SUn9rcijUajWbn8dxsPYVzCX1CiB8BlVMQQuyP5blCiATgv6MS499ZdVwIISqWfwd+EGjf8sVrNBrNNrJrIwohxJ8CZ4FMYAL4ZeD7wO8CeYAZ+DMp5VpbTg9f6x+jtpNaVx3+InAHuACkAgK4Dfz0coJbo9FodgO7Vig0Go1GszU8N1tPGo1Go9kcu7LqKTMzU5aWlm73MjQajeaZ4caNG9NSyqy1HtuVQlFaWsr169e3exkajUbzzCCEGFjvMb31pNFoNJqoaKHQaDQaTVS0UGg0Go0mKlooNBqNRhMVLRQajUajiYoWCo1Go9FEZduEQghRJIR4VwjRFp4c93NrnCOEEP9BCNEdniB3aDvWqtFoNM8z2xlRBIGvSylrgePAV4UQex8655NAZfj2OsqnSaPRaDSrMDwBPLcnWWwZeiLX37aGOynlGDAW/n1RCNEGFAD3V532GZRjqwSuCiHShRB54edqNBrNc4mUkuDUEt52J0ttTvwD82BAXFoCtlOFiDix8UUegR3RmR2eG3EQuPbQQwXAaokcDh/TQqHRaJ4rZNDA1z+Pt83JUruT0IwXAHOulZTmIiy1dhIKUxCmrRUJ2AFCEZ4895fAz69hz73WX7ym3a0Q4nXU9hTFxcVbukaNRqPZDkIuP96OWbztTryds0hfCOIFlvJ0Uk4XYKmxE59ueeLr2FahEEKYUSLxX6WUf7XGKcNA0ar7hcDoWteSUr6FmolNY2Oj9k7XaDTPHFJKAmNuJQztTvxDiyDBlJJA8v4sLDV2EivSMSXEPdV1bZtQhCfC/QHQJqX8zXVO+2vgTSHEnwHHgHmdn9BoNLsJGQjh7ZlX4tDmJDTvA8BcaCP1o8VYah2Y862oj8z1rhHA88EHBIZHSP/cZ7d8jdsZUTQBrwF3hRC3wsf+BVAMIKX8PeDvgJeBbsADfGkb1qnRaDRbSmjex1I4avB1zyEDBiLBRGJlBqkvFmOpsROXkhD1GsHZWdwXLuA6dw7XxUsYCwuYUlNJ+/QPIuK39qN9O6ueLrJ2DmL1ORL46tNZkUaj0TwZpCEJjLgi4hAYcQEQl5FIcmMOSbUOEsvSEOb1OxaklPg6O3G9ew5XSwtLt2+DYRDncJDy4ovYmpuxNp3ccpGAHZDM1mg0mt2I4Qvh65qNiIPhCoCAhJJUUj9RSlKtnfjs5KhbSsbSEu6rV3G1tOBqOU9wTO28W+rqyHzjDWwfOYulrg5herItcVooNBqNZosIOr3h3oYZfL3zEJIISxyWqgwstQ4sVRnEWc1RrxEYHcXV0sLiuXN4rl5D+nyI5GSsJ0+Q8tWfwXr6DOac7Kf0Fym0UGg0Gs0mkSGJf2gh0tsQnPAAEJ+VhO1kvqpSKk1FxEXZUgqFWLp9G9e5FlznzuHr7ATAXFRE+o/+KLbmZpKPHsGUED1nIaXE6w6QZIt+3mbQQqHRaDSPgLEUxNupKpS8nbMYniCYBIllqVgb92CptWPOTIp6jdD8PK6LF3Gda8F94QKhuTmIiyP58GGyv/ENbB85S0JZWdRtKYBgIMRo5xz992YYuDuNlPDavz2x4fMeFS0UGo1GEwUpJcHpJSUM7U58/couw2SNx1JjV7eqDEyW9T9OpZT4e3pUhdK5FjwffAChEHHp6diaz6hE9KlTxKWmbrge16yPgXvT9N+dYbjdSdBvEG82UViTQUlDJtKQu9PCQ6PRaHYSq+0yvO1Ogg/bZdTYSSiKbpdh+Hx43ntfiUNLC4HhYQASa2pwfPnL2M42k7RvHyIuevOcYUgm+xfovzvNwL0ZpodUxVSK3ULNiTxKGzIpqEon/gk24Wmh0Gg0GqLbZdhOhe0yMqLbZQQmJnG1nMPVch735cvIpSWExYL1+HElDs1nMOflbbgWnyfA4H0nA3dnGGidwesKIEyC3D2pnPhsOSUNDux50ZvwthItFBqN5rlESklg3IO3fQZv2+bsMqRh4L17N1Kl5LvfBkB8fh7pn/1HKhF97BgmS3SBkVIyO+ah/940A3dnGOuZRxoSi9VMcb2d0vpMivbasWxQMfWk0EKh0WieG2TAwNs7p7aU1rLLqLFjzrdF3VIKuVy4L15SW0oXLhCamQGTiaSDB8n6hV/AdraZxMrKmBLRI51zDNyZpv/eDIvh7S1HoY1DLxVT0pBJTlkqpljcYKWE6U6Y6YGal2N/Q2JEC4VGo9nVhBbCdhltq+wyzKvsMqrtxKVGLyn19fWpprdzLXiuX4dgEFNqKrbTp7GdVYno+IyMDdfimvXSf3eGgXurEtEJJgpr7Bz+RAkl9Q5sG2xvRQh4YeAidL6tbnMDkJgG3+yFuF1i4aHRaDRPAmlIAqMultoesstID9tl1NhJ3JMe3S7D78dz40akSsk/MABAYmUFji9+AdvZsyQdOLChXYZhSCb6Fhi4q6KGmeFwItphofZkPiUNDpWINseYiF4Yha53lDD0noOAB+ItUNYMTV+Dyo9vuUiAFgqNRrMLMHwhfN2zShw6nBiLYbuM4rBdRo2d+JzodhnB6WlcLedxtbTgvnQJw+1GJCSQfOwYGT/xGrbmsyQUFmy4Fq87wNB9J/33phm858TrVonovPI0TnyunNL6TDLyoq9l5Q8LwchN6HobOv8exu+q42lFsP9VqPoElJ0Gc/S+jcdFC4VGo3kmidhltDvx9cwpu4zEOCzVsdllSMPAe79NVSmda8F7V30Ix2dnk/oDP6C2lI4fx5ScHHUdUkqcY24G7s7Qf3ea8d4FlYi2mSmpd1DS4KCo9hES0Utz0PN9FTl0vQOeGRAmKDoGH/1lJQ7ZtfCUKp5AC4VGo3lGkIbEP7gQmRMdscvITMJ2Ih9L7cZ2GYbbjfvKlUi+ITg1BUKQtG8fWT/3NWxnz5JYU7NxItqvEtH9d1WV0qJTJaIzi2wc+ngxpQ2ZZJc+YiJ6OdcweAVkCJIyoOJFtZ1U8VFItsf+Zm0xWig0Gs2ORdllhHsbOpwrdhmlqVh/YA+WmgzMWdG/8fuHhiI+Sp733kMGAphsNqynTmE724zt9GniHY4N17Lo9DIQtsoYbp8lGFiViP5kCSX1mdgyEmP7wyKJ6HfUltKcyoGQXadyDVWfgILGJ5Jv2Aw7YxUajUYTJjDliUx7i9hlJMdjqbZjqbVjqczAlBTFLiM87c11rgVXSwv+nh4AEsrKyPjxH8d29izJhw8hzNG3ggxDMtE7H/ZRmmEmnBRPzbRQeyqf0noH+ZtKRL8Dve+unYhOL9r4OtvAds/M/kPgU8CklLJ+jcfPAv8T6Asf+isp5bee3go1Gs2TRoYMfH0LkTnRweklAMy5yaScKcJSu7FdxlrT3jCbsR5pJOPzyoE1oaRkw7V43QEG78/Qf2eGwfsz+NxBTCZBXkUaJz9XQUmDg4zczSSi34bxO+p4JBH9cSg9DQnRI6KdwHZHFH8E/Dbw7SjnXJBSfurpLEej0TwNQu4A3g5neEspbJcRJ0gsT8fWlL+hXcaHpr3dugVSEpeZScrHwtPeTjYRZ7NGXYeUEueoO+KjNN4zj5SQlGKmtCGT0gbVEZ0YJYJ5AO88dH8vnIj+LnimH0pEfxyy9z7VRPRWsK1CIaU8L4Qo3c41aDSaJ4+UkuCEJ9Lb4B9cCNtlmEnet8ouI3H9bZx1p73V15P5Mz+D7exZLHV7N5z2FvSHGO6YVVVK96ZxOVV3dmaRjcOfLKWkwUFOSWrUCGbVHwbTXSrP0PWOSkQbwR2ViN4KtjuiiIUTQojbwCjwi1LK1rVOEkK8DrwOUFxc/BSXp9Fo1uIBu4x2J6G5sF1GgY2UF4pJqt3YLmO9aW+2ppPY3vwq1tOnMWdvPO1t0emNNL0Nt88SChjEJ8ZRVJPBkZfLKKl3YE2PMREd9EF/uCO6622Y7VfHs+vg5M8qcSg8smMS0VvBTv9LbgIlUkqXEOJl4H8AlWudKKV8C3gLoLGxUT69JWo0mmVCCz687WpOtK9r9kG7jBeKsdRkEJe6/gfyhtPezjaTfGTjaW9GyGC8byHS2+AcdQOQmpVE3alwR3RlBnFRurMfYGHsoY5o90oielkcdmgieivY0UIhpVxY9fvfCSF+RwiRKaWc3s51aTQaRVS7jMM5JNVubJex5rS3+HiSDx0i+5vfxHa2OaZpb15XgIFW5aM02DqDzxNORFemcfKHKihtcJC+QXd2BMOA0Zvh3oa/fygR/cozlYjeCna0UAghcoEJKaUUQhwFTMDMNi9Lo3muMfwhfF1zka5oY9G/Ypfx8VKSaqPbZaw77S0jA1tzs+qIbmracNqblJKZETcDYWvu8d6VRHTZ/kxK6jeRiO75vipf7Xpn1ySit4LtLo/9U+AskCmEGAZ+GTADSCl/D/hh4KeFEEFgCXhFSqm3lTSap0xw1hvpiPb1zkFwlV1GeBRonG397SA17e29yJZSYGQEgMTaWhxf+TIpZ89iaWjYcNpbwB9ipH02MiPaNavyHlnFKRx+uZTS+kyyS6KX0kZYTkR3reqINoJgSYfKj+2aRPRWIHbj525jY6O8fv36di9Do3lmkYbEP7QYTkTPEBgP22U4LMpHqcZOYll0u4x1p72dOIHt7Fk17S03d8O1LMwshXMNM4x0qkS0OTGOolo7JQ0OlYhOe8RE9HK+YTbcopVdB1Uv7cpEdKwIIW5IKRvXeuz5ezc0Gs2aGN6wXUbbarsMSCxNI+0HytRQnyh2GetNezPn55P+2c+qRPTRoxtOezNCBuO9KzOilxPRaVlJ1J3Op7Qhk/yK9EdPRHe9Az3vPpSIfnPXJ6K3Ai0UGs1zTGB6CW/bTNguYwEMuWKXEd5SimaXEVpYwH3pktpSunCBkNO5Mu3t679AytmzJFRUbJhAXnL5GWx1MnB3msH7zlWJ6HSafjiP0oZM0nNiTBwbIRi5sVK++oA19+eVj9JzlIjeCrRQaDTPETJk4OtfiPQ2LNtlxOckk3KmAEuNnYTi9ZvNIono5WlvN2+qRHRaGtbTp1VH9KmmDae9SSmZHnYxEJ72NtG3KhF9IIvSemXNnRBrInppdqUjuvsfwtbccSoR/eKvqKjhKVtz7ya0UGg0u5yQOxDeUprB2zmL9D5kl1FtJ96+/naQ4fWuJKJbWlYS0dXVOH7qp7CdbSZp//4NE9F+b5Dh9lkGwltK7nk/ANklKTS+XEpJQybZxY+QiJ68H44a3oGhayANSLKHE9EvqUR00sbjSTUbo4VCo9llROwywg6sq+0ykuozVW9DRUZUu4zA2FgkanBfvYr0ehFJSVhPnMDxla+oRHRe3oZrmZvwRHINo11zGCFJgiWOor0qCV1S7yB5g3nVEfwe6DuvtpO6vgvzQ+p47j44/XUVNRQcAlOMbq6amNFCodHsAmTAwNc7FxGHR7XLkMHgSkd0S8uDHdE//MPYmptJPnoEU2L06qJgIMRo11xkS2l+Sm1tZeRZ2f9CESX1DnIr0oiLUi31ALP9K30N/Rcg6AWzFco/Ame+oSKH1I0FS/N4aKHQaJ5RQgt+vB3h3obuWaQ/bJdRkU7KC0UkVduJi1I2GpydxR3uiHZdvIgxP7/SEf2Nb2D7yNmYOqIjA33uzTDc7iToN4gzmyisyWD/R5U4pGbGONM5FIDBq+HehndgukMdt5dD408qYSg5CfExlsNqtgQtFBrNM4KUksCIK9IRHRgO22WkJZJ8KEcN9dmThlhnkI6UEl9HRyRqWLp9GwyDOIeDlI98ZKUjOiUl6jpW+ygN3FsZ6JPisFB7Io/iegeF1RnEJ8S4BeSaVAnozrdV+apvHkxmKG2Cw19UHdGO8pjfJ83Wo4VCo9nBGP4Qvu5VdhkLYbuMohRSP16CpcaBOcogHcPjUdbc51pwnT9PcHwcAEtdHZlvvIHtbDOW+voNrbmXFv0MLvsorS5f3dRAHwPGbq00vY1+AEhIyYO6z6ioYc9ZSIwuWJqnhxYKjWaHsWyX4W134u1ZZZdRFbbLqI5ulxGZEd3SomZE+/2YrFasJ09i+9k3Y7LmllIyPeSKJKIn+lVCPCk1gbIDWZTUOR7RR2lBjf9czje4JwEBhY3wkV9SXdG5+3T56g5FC4VGs81E7DLanXjbVuwy4hwWbMfysNQ6SCxNRcSv/a1fBgJ4btwMD/Rpwd/bC0BCaSkZr76qOqIPH0ZsYM3tXwoy1O5UW0qtM3jmVfSSU5rK0U+pmQ1ZG4wkXVnUej5KaVD+UbWdVPEiWDMf7c3SbAtaKDSabSBil9Eetstwr7LLeLkMS62d+MykdbdygtPTuM5fwNXSgvvSJQyXC2E2k3zkCBmvfD6mGdFSynD5qtpSGutW5auJyfEU7bVTUu+geO8jlK8GvDBwMRw1rB7osxdOvKnEofDoc+mj9Kyj/4tpNE8JZZehTPZ8favsMqoyVCK6MgNTsnnN50rDwNvaGtlS8t67B0B8djapn/yE6og+cQKTNfqM6GAgxEjncvnqNAvTXgDs+VYOvBguX92ThinW8tX5kZUKpb4WCHggPgnKzoQH+rwE6Xri5LOOFgqN5gkRscsI5xuCU6vsMk4XYKm1k1CUiohbO2oIuVy4L11WW0rnzxOangYhSNq3j6yf+xq25mYSa2tjL1+9O81w+yzBgEF8gonCGjsHXyqhpN5BSpTO7AcXFYSR6ysd0RNKsEgvhgM/Hh7ocwrMMZbDap4JtFBoNFtIVLuME/lYata3y5BS4u/rj+QaPDduQCCAKTUV26kmFTWcPk28Pfp8hFDIYKJ3PrKlFBkDmmmh9lQ+pfUO8qvSiV+njPZDeJyryle/p3yVRBwUn4CPfUt1RGdV60T0Lma7Bxf9IfApYFJKWb/G4wL4LeBlwAN8UUp58+muUqNZHyklwUlPZBSofyBsl2EL22XU2EmsTMeUuPY/NcPvx/Pe+xFxCAwOApBYWYHjCz+BrbmZpIMHEfHR/6l6FvwM3p9h4K4qX/UvBTHFCfIr06k9mUdJ/SOMAZVSRQrLUcPw+8pHKTlTOa9WvgTlL0BS+iO/X5pnk+2OKP4I+G3g2+s8/kmgMnw7Bvxu+KdGs23IoIGvd56lthm87U5C4Slr5nyrssuosWMuWN8uIzAxERaG87ivXEF6PIjERJKPH8P+xS+Q0tyMuaAg+hoMydTQYiRqmAwLVHJaAuWHsiitz6SwJiN291WfS+UYut5RPkoLyviPvANhq4yPQ/5B2KDfQrM72VahkFKeF0KURjnlM8C3w+NPrwoh0oUQeVLKsaeyQI0mTGjRH2l683U9ZJdxtoikmvXtMmQoxNLtO5Fcg69NDfSJz88j7TOfVltKx45hSoq+r+9bCjJ036lmRLc6WQo33+WWpXLsB8soqc8ks8gWW9QA4OxdqVDqvwghPySkQPlZOPvPlQtrysYT6DS7n+2OKDaiABhadX84fOxDQiGEeB14HaC4WFdZaB4PKSWBUTfetpm17TJq7FjK17fLCM3N4bp4SZWvXrhAaG4O4uJIOniArK//gkpEV1ZG/VCXUjI75qH/3jSD92YY657HMFT5anGdcl4trrOTFKX57gGCftXPsNwRPdOljjsq4ejrakup+ATEx3g9zXPDTheKtf4VrTnkW0r5FvAWqJnZT3JRmt1JVLuMl0rUKNA865of7lJKfJ1dkVzD0gcfKB+l9HSsZ9RAH9upU8SlpUVdQ8AfYqRjNmKytzijylcdhTYOvFRMab2DnLLU2MtXFyfC20lvQ8858C9CXIKqTDryZdURbd/zqG+V5jljpwvFMLB6mG0hMLpNa9HsQoJzYbuMNifennkIGjHbZRhLS7ivXcN17pzyURpVgW5ibS2O17+iEtH79m040GduwsNA6wyD92YY6ZwjFDSIT4yjqCaDw59Q5au2jBjLV42Q8k5ajhrGbqnjKfnQ8EMqaihrhkRbzO+RRrPTheKvgTeFEH+GSmLP6/yE5nGQhsQ/vKiEoc1JYFyVjiq7jFwstXYSS9PWtcvwDw7iajmP6/x5PNeuIf1+RHIy1hMnsL3xBrbmZsw5OVHXEPSHGOmaUwZ7q2Y2pOckU99cQEmdg/zKdOLMMUYNHqcaA9r93VVjQE1QeARe+JeqtyGnXpevajbNdpfH/ilwFsgUQgwDvwyYAaSUvwf8Hao0thtVHvul7Vmp5lnG8Abxds0qceiYxXAHwAQJJWG7jBo78Vlr22UYfj+e99/Hff48rpbz+Pv7gWUfpVewnj6jBvps4KM0P7WkhKF1hpGOcNOb2URBeGZDcZ2DtKwYm9QMA8bvqOqkrndUA5w0INkBFR9TSejyFyA5er+FRhMr21319OoGj0vgq09pOZpdRHB6SU17a3fi65uHkEQkxWOpziBpA7uMwNiY8lE6v6p8NSGB5GPHyPjxH8d25vSGPkrLk94G7zkZaJ1hbkIZ/aVlJ7H3VD4lj9r05p1Xsxq6vqsiB9eEOp5/CM58U20p5R/QY0A1T4SdvvWk0cSEDBn4BxYio0AjdhnZydhOFZBUYyeheG27DBkIsHTrFq5w1BAZA5qfr8pXz5xR5avJyVHXsDC9EjUMd8xGJr0VVKXTcLaA4joH6dnRr7GyKAmTbSt9DUNXH3RfrXwJKj4Ktuh24RrNVqCFQvPMEnIH8HXOKnHomEV6g8ouY08atuN5akvJsfZ2TnBqCteFiypquHQJY3FRjQE9fFiNAW0+Q0J5edTy1VDAYLRnJdcwG7YHT820UHsyn+I6OwXVGZhjnfTmc0Hf+VVNb8PqeE4DnPyaEofCI9p9VfPU0f/HaZ4Zlu0yvO1qTvQDdhl1DpJq17fLkKEQ3rt3I1GDt7UVgPisLFI+/pKKGk6eJM4WvRpo2WBvsHWGofZZgr4QpnhBQVUGdacLKKl3kJa9vj34Q38QzHSHheEdGLj8UNPbP1UzG1LzN/N2aTRbhhYKzY5m2S5jubch5FR9BeZ8KykfKSKp1rGuXUZwdhb3xUsqalhuejOZSDpwgKyf/3lszWdIrKmJHjUEDcZ65iPisGywl2K3UHMsl5J6h4oaEmOMGgJLqgt6WRyWZzZk1cCxf6KihqLjuulNs6PQQqHZcYQW/Xg7VNSwbJdBvAlLRTopzYVqS2kNuwwpJb62tkjUsHT7tmp6y8jA1nwG65kz2JqaiEuPbmbnml2OGpwMtTsJeEMRg72aE8pgL+b50KDEYLlCqe88BL1qZsOeZjWzoeJjkBE9Oa7RbCdaKDTbTsQuIxw1BIYWAYhLSyD5YLYaBbonDdMae/2hxUXcl6/gOt+C+/wFglNTAFgaGsj86Z/G1nwGS11d1Ka3UMhgvGeewVbVDT0zoqIGW0YilUdyKKlzKIM9S4z/XIK+sFVGWBymVXIc+x44/EVVvlpyCswxNtFpNNuMFgrNtrDaLsPb7iS0bJdRmELqx0qw1K5tlyGlxN/dHYkaPDdvQjAYmdlgPXMG26lTxGdGn8XsnvNFuqGH2pz4vSFMJkFeRRonPldOSZ0De/7adh1rMj+iyla7vgu958DvgrhEZZXR+FNKHBzlm3y3NJrtRQuF5qkRnPNFhMHbPafsMhLisFSlY6lxKLuMlA/vzRseD+6r13Cdb3nQKqO6GseXvoSt+QxJBw5EndmwnGsYuj/DQKuTmbDJnzUtgYrD2RTXOyiqscduyx0KwNB7KxVKkyo5TloR7Pt82CrjNCREH02q0TwLaKHQPDEesMtodxIYC9tl2MN2GTV2EsvWtsvwDwysWGW8996KVcbJsFXGmTOYc6NbYC9ML6ntpFYnIx2zBHwrUcPxf7SHkvpMHAWPEDUsTiiLjK53VPObbx5M8eFJb/9GiYOe9KbZhWih0Gwpyi5jZUtpxS4jlbRPlmGpXdsuw1hawvPee6q34cJ5AgNq0lvCnj1k/NiPqajh8OGoVhkBf4jRzjkGW9WUt+Vu6BS7hapjuRTvtVNY/QjDfEIBNd2t+x9U1DB+Rx235cLeTyth2HMWLKmP+jZpNM8UWig0j01wZikyCvRDdhk1dixVH7bLUPOh+3BfuIDr/AU877+vogaLBeuxY9hf+wnV9FZUtM6rrsxrGLyvhGE07Ly63A1df6aA4jp77CNAIZxr+Ad1621RUYOIg+Lj8NF/pfoacvfpqEHzXKGFQvPIyJAM22XMPGSXkYStKWyXUfJhuwzD7Va23OfP475wkcCIGreZUF5OxquvYj1zmuTGRkyJa0+KAzXlbbjdyWCrk8HWGVzhMaQZucp5tXivnfzKdOJj7YYO+mDw6oo4TN5Xx1MLoO4zqnR1T7OyztBonlO0UGhiwvAE8HbOqsjhIbsM65Ug570AACAASURBVPE8ktawy5BS4uvqUlHDhYt4btyAQABTcjLJJ07g+MqXsZ46TULh+vOhpSGZHnapCqXWGcZ7F5CGxGyJo6jGTuPLdor22kldx6pjTWYHwpbc31NRQ8ANJjOUnFS5hooXIbtWRw0aTRgtFJo1kVISnFpSo0BX22VYlV2GpcaOperDdhmhxUXcV65ExCE4Pg5AYlUV9p94DdvpMyQfOoiIkmtYWvQzeN/J4P0Zhu47WVoMAJBZZOPgS8WU1NnJ2ZNGXKxT3gJeGLiohKHruysjQNOLYf8rqnS19LQe5qPRrIMWCk0EGTTw9c3jbXvILiNP2WVYauwkFKY8YJchpcTX3o7r/AXcFy7guXVL9TXYbFhPnsT25lexnjoVtUIpFDKY6Ftg6L7aTpocXAQJFquZor12SursFO11kJwao62FlODsDVty/4OyzAgurfQ1HPkpFTU4KnTUoNHEwHYPLvoE8FtAHPD7Uspfe+jxLwK/DoyED/22lPL3n+oidznLdhneNiferjmkP/SgXUa1nfj0B3MGobk53JcvqwqlixcITU0DkLi3FsdP/iS2M6dJ2r8fYV573oOUkvnJJYbanAy1ORnumCXgDSEE5JSlcfRTZRTXOcgqTsG0hofTmvjd0HchnGv47oqHkqMCDn9B5RpKTkJCjDbfGo0mQkxCIYT4l8AfSSmHVh17XUr51mZfWAgRB/xH4GOo2djvCyH+Wkp5/6FT/1xK+eZmX0fzIFJKAmPuSG+Df1h9e49LTSD5YJbqbShPf8AuQxoG3tZWXBcu4D5/gaU7d8AwMKWlYWs6ifX0GaxNJzFnrz8bwesOMNw+GxGHxRkVraQ4LFQeyaGoxk5hTQYW69rissYfAlMdK+M/l51XzclQdgZOvKmiBnvZY71fGo0m9ojiZ4FXhRBflVK+Gz72BrBpoQCOAt1Syl6A8FzszwAPC4XmMTH8IXw94d6GtrBdBmAuSiH1xbXtMgKTk7gvX8Z98RLuy5cJOZ0gBJb6ejLfeAPr6VMk7du3rodSKGgw0TfPUNssg/edTA0sICUkWOIoqM7g4MeKKdprJ22dEaRr4nEqe4ye76uGt+V5DVm1cPR1JQwlJyF+/aopjUbz6MQqFCOoD/HvCCH+Qkr568Djbu4WAEOr7g8Dx9Y474eEEGeATuD/WB3VaNYnOO+LCMMDdhmV6Vhq7Viq7Q/YZRheL+4bN5QwXLoUmfIW53BgbWrCduY01qYm4u1rz2GWUjI34VERw30nI51zBHzL20mpNL5cSlGtneyy1NiT0JGGt+8pcRj9AJW8SIOyZjjzi0oc0tfvtdBoNI9PzDkKKeWgEKIZ+F0hxHeAR6hHXJO1hEY+dP9vgD+VUvqEEG8Afwy8sObFhHgdeB2guLj4MZf27BGxywiLwwN2GUdzsdQ+aJchpcTb2RkRBs/160ifD2E2k3T4MNm/+HWsTU0kVlcjTGt/sHtdAYbanRFxWO5pSM0Md0LX2imoTidxndnUH/4jwknonu+rW995Za4n4qCwEc7+cyh/AfIP6ilvGs1TJNZ/bdcBpJRe4EtCiK8Chx/ztYeB1V8FC4HR1SdIKWdW3f1PwP+13sXC+ZK3ABobGx8WnF2J4Qvi65oL9zY4MVwB5cBakkraJ0ux1DoesMsIOp24L13GfUmJw7Ild0J5ORmvfB5rU5NqeFtnNnTAF2KsZ47h9lmG22eZGlL5jYSkeAprMjj8STtFtRmkZT1CwnhpTgnCsjjMDajj6SWw70eVMJSehqToMyQ0Gs2TIyahkFJ+5aH7/xGViH4c3gcqhRBlqK2tV4AfW32CECJPSjkWvvtpoO0xX/OZJzizpGZEtzvx9YbtMixhu4zaB+0ypN+P5733cV+8iPvSJbz3VfonLi2N5JMnsJ06hfXkScx5eWu+VihoMNG/wHD7LCMds4z3zmOEJKY4QU5ZKkd+oIzivXayS1IwxbydFITRmyvCMHwdZEiN/yw7owb5lL+gZjfo0lWNZkcQVSiEEHf58HZQBCnlvs2+sJQyKIR4E3gbVR77h1LKViHEt4DrUsq/Br4mhPg0EAScwBc3+3rPKit2GU687TMEJ8N2GVlJ2JrySapxROwypJT4e3txX76ixOH995EeD8THk3RgP1k//3NYm5qw7N27ZhJ6uQt6uH2W4Q4no93zBH0hEJBVlML+F4oorMkgryI99tGfy9tJfS0qAd3XAt55QEDBITj9C0oYCo9AXIxbVBqN5qkipFx/l0YIsTyfUQD/C3h59eNSyoEnt7TN09jYKK9fv77dy9g0EbuM9rBdxlLYLqMsDUuNXdllZKoUUWBsDPeVq7ivXsFz9RrByUkAzCXF2JqasJ46RfLRo8TZPtx1vJyAVsIwy0jnLD53EFDeSYXVGRTUZFBQ9Qhlq6DsuPtalD1GXwvMh+sPUguUKJS/oFxXk9dOjGs0mqePEOKGlLJxrceiRhSrhUAI4dupwvCss2KX4WSpfUbZZRhhu4xaO5ZaB5bKdEyWeNXsdu2CEoYrV/H39wMQZ7djPX6M5OPHsZ44sabrqpSS2XEPo11z6tY5i3telcra7ImU7c+isDqDwuoMrOmPUGLqnYf+SyviMBXeIbSkq+E9TT8Hez6iJrzp7SSN5plDl45sExG7jPCc6NDMKruM5iIstcouQ/q8eG7cZPq3r+C+clXlGaREJCeTfKSR9M9/HuvJEyRWVn6oOkkakplRV1gU5hjtnov4JiWnJZBfma6EoSaD1MxH6GcIeGH4vZWIYeSmyjPEJyk77v2vKMfV3H1ginGLSqPR7Fg2ylEcWnU3SQhxkFVlrVLKm09qYbuRkMuPt30Wb/uMssvwrbLLOF2IpcaOKVngvXOHhb/9OzxXrrJ06xYyEACzmaT9+8h886tYT5wgqaHhQxYZRshgasgVEYWx7jl8HrWVlGK3UFznIL8ynfzK9EdrdAsFVA9D/0UlDINXIehVZasFh1WeoawZio7qZjeNZheyUUTxf6/6fRz4zVX3Jev0NGgUEbuMcJWSP1xOakpNIHm/ssswF1rw3b+H5/13mfmd91m6fRvp84EQJNbWkPHaa1hPHCf58OEPla163QEm+hYY751nvHeeib4FAr4QAGnZSZQfzCK/Mp28yvRHs+EO+mDkhtpOGrioZkMH1LQ4svfC4S+piKGkSU9302ieAzbKUXzkaS1ktyADIbw983jbZvC2OwnNP2iXkViajH+snaXr/4v5v7yO984dFTGYTFhqash45RWSjx4h+fBh4tJXegekIXGOuiOiMN47z+y4+vAWAhyFNqqP50YiBmvaI3yz93tUB/TAJSUOw+9DSDXPkVMPB1+D0iYoPgm2rC17rzQazbPBRltPqUCOlLIrfP9HWOnIfltKOfGE1/dMEJr3qQqlNie+njlkwEAkmEiszMB6PAHp6cfXdg3nW9dZam2FYBDi4rDU1ZHxE6+RfCQsDCkpgIpEXLM+pj6YYnJwgamBRSb6FyLbSInWeHL3pFF1NJfc8jSyS1JIsDxCusk9raKEoWtqG2nkBhgBECbIbYAjXw4LwwldmaTRaDbcevoN4DIQnvTCrwL/GyUWJ1HGgM8d0pAERlwshaOGwGjYLiM9kcTyBKRvEF/3+8x9+yaBoXBpqNlMUl0djp/8SZKPHCHp4EHibNaIKAz0LDI1OMXkwAJTg4uRpLMwCex5yZQfyiZ3Txq5e1IfbQa0EVLjPYfeU7fh91RfA4ApHvL2w/GfVnMaio/rkZ8ajeZDbCQUR4B/sur+opTyZwGEEBef2Kp2IOvZZcSlhoizjuHvv4brnQuqwQ2Iy8ok+cBBMl59laQDB7DU7cXrEzhHXYyPuZn5n0M4R9w4R134vSqvsCwKJQ2ZZBenkFWcQmahLfb5z1KqnoXRWzB2S3U9j9xQfkkA1iwoOgaHvqB+5h8A8+Nadmk0mt3ORkIRLx/syHtt1e+73nwn6PSqUaDtaksJAxBB5NIg/t5r+PveU/OW4+OxVFeT9tnPYao7QKCgGjc2xqaXmJ9aYuF7S8z96XuRKAHU9pEj30b1sVzs+VYyix5RFAwD5gdh7PaKMIzegiWnelzEQc5eVapaeFRVJGWU6j4GjUbzyGwkFIYQIldKOQ4gpbwHIIQoQH1s7io8c7OE2mfwfDBCYDQAAZUQDnqm8U514XIO4vYvQF4xVB4gePRlfEkZeI1E3AsBPJM+giMGMBi5pi0jkbSsJMr2ZWLPt2HPt2LPt5KcmhDb9pGUMD8MU+0w2Rb+eR+mOpVIgdpCyqqFmpch74ByV82p09GCRqPZEjYSil8H/kYI8XXgg/CxQ6jcxa8/yYU9bYKuJf7D//O7JMsECkMO0vwOfP44JgMSt5EGaY3qtowb4vwmrOkmbOkmckpSsO7PxJqeSFp2MmmZSaRmWjaOEII+cE2Gb+MwO6DGeM6Ff84OqHnPy9hyIKsGDr2mfubuC4uC5Um8LRqNRrNheex/EUJMA/8WqEP1TrQC/0pK+b+fwvqeGjLBRIUvlbH4Be7GDyDNAyQkB8kLuagMzVIsJ0gxLZBoWiIx2UxiihWzzYpItKlv7gnJagynkQhjEsZQ0QCoiiKfS+UK/C4133lpFlwTYYO8h0iwqW0iR0V4nOceyK5VwqCrkDQazVNmI1PAV4F3HpoLseN5XFPApaUluru76ezooKurE6/Pj0lAiT2BqrQAVZZZHKFJ8Myo7Z/AkupFCHhUhCAEIFZ+muIgMQUSrEoEEqxqvoItB6zZYMtWv9uyIL1UiYHOJWg0mqfIpk0BgRLU+FMz8D1Uaex7Mpq67AKSkpJoaGigoaGBUCjE8PAwnZ2ddHZ28nbvFG+TjsNRTlVVFVVVVRQXFxO3zuxojUajedaJGlFEThIiBXgR+ARwFDVA6O/ZoU13T9Jm3Ol00tXVRWdnJ/39/YRCISwWCxUVFVRVVVFRUUHyOhPiNBqNZqcSLaKISSjWuOBe4JPAS1LKjz/m+racpzWPwufz0dvbG4k23G43QgiKiooi0UZWVlbszXEajUazTTy2UAgh/hL4A+DvpZRbVhYrhPgE8FuoCXe/L6X8tYceTwS+jZrPPQN8XkrZv9F1t2NwkWEYjI6ORkRjfHwcgPT09IholJaWEh+vnd01Gs3OYyuE4kXgS8Bx4DvAH0kp2x9zUXFAJ/AxYBg1Q/tVKeX9Vef8DLBPSvmGEOIV4LNSys9vdO2dMOFufn4+skXV29tLMBjEbDZTXq5yG5WVlaSEvZ00Go1mu9myrSchRBrwKvBLwBDwn4D/IqUMRH3i2tc6AfzK8taVEOKfA0gpf3XVOW+Hz7kihIhHWZ1nbZRM3wlCsZpAIEBfX18k2lhYWAAgPz+f6upqqqqqyM3N1VtUGo1m23icqqfVF3EA/xhl4/EB8F+BU8AXgLObWFcBSmyWGQaOrXeOlDIohJgHHMD0Gut7HXgdoLi4eBPLeXKYzebI9pOUkomJiYhovPvuu7z77rukpKREzikrKyMhIWG7l63RaDRAjEIhhPgroAb4E+AHpZRj4Yf+XAix2a/ua319fjhSiOUcdVDKt4C3QEUUm1zTE0cIQW5uLrm5uZw5cwaXy0V3dzcdHR3cvXuXGzduEB8fT1lZWWSLKj1919tqaTSaHUysEcVvSym/v9YD64UqMTAMFK26XwiMrnPOcHjrKQ1wbvL1diQ2m40DBw5w4MABgsEgAwMDkWijq0u5u+fk5ESijYKCAkwPzcbWaDSaJ0msyezPrXF4HrgrpZzc1AurD/5O4KPACCqZ/WNSytZV53wVaFiVzP6clPJHN7r2TstRbAYpJdPT0xHRGBwcREpJcnIylZWVVFVVUV5ejsWiPZ40Gs3jsxU5ip8CTgDvhu+fBa4CVUKIb0kp/+RRFxXOObwJvI0qj/1DKWWrEOJbwHUp5V+jSnL/RAjRjYokXnnU13lWEUKQlZVFVlYWTU1NeDweenp66OzspKOjg9u3b2MymSgpKYkkxO127QOl0Wi2nlgjir8BvrzchS2EyAF+F/gycF5KWf9EV/mI7IaIIhoP24pMTU0BkJmZGdmiKioq0rYiGo0mZraij+KulLJh1X2B2naqF0J8IKU8uHXLfXx2u1A8jNPpjIhGf38/hmFoWxGNRvNIbMXW0wUhxN+imu0Afgg4L4SwAnNbsEbNY2C32zl+/DjHjx/H5/NFtqi6urq4d++ethXRaDSPRawRhQA+h+qbEMBF4C93qovs8xZRrIe2FdFoNLHyWFtPYauNt6WULz6JxT0JtFCszVq2IgkJCQ/Yithstu1epkaj2QYea+tJShkSQniEEGlSyjXGsWmeFdLS0mhsbKSxsRG/309/f38k2mhrawOgoKAgEm1oWxGNRgOxbz39N5Qh4HcB9/JxKeXXntzSNs+TiigMQ7LgDTDt8uN0q9tSIMiS38DjD+INhPAHw+a6QiDUD8xxJqwJcVgT4yO3tCQzWSmJZNoSSIzf3uokKSXj4+MR0RgZGQHQtiIazXPEVlQ9fWGt41LKP37MtT0RNisUPaO9eP3JjLtCDC/4GJ71MuRcYnjOw8SCj1m3n6Cx9WmZVEs8WSmJZKUkUpiRTLE9fHOonw5rwlP9Zu9yuSJbVD09Pfj9/gdsRaqqqkhLS3tq69FoNE+eLXGPFUIkAcVSyo6tXNyTYDNCEQwZ/MZrn8IcWrHHMADDZAKTGZGYjMmSjDnJSqLVSrIthZSMDFLsDtIcDtKzMnFkZpGe6cAU7l+QUiIl+EMGHn8Ity+I2x/E7Qsy5wkwtehjatHHtMvHlMvHxIKP4VklSqtJscRTmW2jKieFypwUqnLU79kpiU9cQFbbinR0dDA3p4rccnNzI6KRn5+vbUU0mmecrYgofhD4DSBBSlkmhDgAfEtK+emtXerWsCmhMIL863//S/iDs/iDc7j9c3j9S5gkxIcESSEzGaSSIpOwBOIw+Qz8LhdGKPTAdYTJRGpWNuk5eeFbLum5+aTn5pGRV0BcDBVG3kCI4VkPg04PAzMeeqfcdE4s0jmxyKxnxdE9LcnM3rxUGgrTqC9Io6EgjRJ7MibTkxGP9WxFrFZrxFZkz5492lZEo3kG2QqhuAG8AJxbbq57uAlvJ7FVOYp53zy98710z3XTNdtFu7Oddmc7S8ElACwmC3uTK6iOL6U0Pp88w068K8j8xDhzE+PMTYzic0dSOpji4rHnF+AoKiFz1S0tOwcRwzdyKSXTLj9dYdHomHBxf3SetrFF/CGVG0lJjKeuIJWGgjT2F6VzqDiD/PSkx34v1sLj8dDd3U1nZyfd3d14vV5MJhOlpaWRaEPbimg0zwZbIRTXpJTHVndhCyHuSCn3bfFat4QnWR4bMkL0L/Rzf+Z+5NbubMcT9ABgt9jZl7WP/Vn72Z+1nwpLCb7pOZxjI8wMDTA9PMjM0ADzkxORa8YnJpJVXErOngpy9lSSs6cCR0FRZAtrI/xBg86JRe6NzHN3ZJ57I/O0jS9GEut5aRYOlWRwqDiDwyUZ7M1LJSF+a7eKQqEQQ0NDkWhjelqNDNG2IhrNs8FWCMUfAN8D/hmqK/trgFlK+cZWLnSreNp9FCEjRM98D7enbnN78ja3p27Tv9APQLyIp8pexf6s/RzKPsThnMNkJWfh9y4xMzzI9NAA04MDTPb1MNHXQ8CropX4hESySlbEI6+iGnt+QUyRByjxaB9f4MbALDcH57g5MMvInLp2YryJ/YXpHCxJp7HEztFSO2nJ5i19T6LZilRXV1NeXq5tRTSaHcRWCEUyavzpS6jO7LeBfyOl9G7lQreKndBwN+ed4870HW5N3uLO1B3uTN+JbFmVpJbQmNPI4ZzDNOY0kmfLA0AaBrPjo0z0dDHR181E74PiYbHayKusJr+qlryqGvIqqkhIiv3Ddmx+iZsDc9wcnOXGwCyto/MEQhIhoCY3lWNldo7vsXO0zIHdunWlsF6vl97e3ohweDwehBAUFxdHoo3MzEzds6HRbCNbNjP7WWEnCMXDBI0gHc4Ork9c5/r4dW5M3mDRvwhAga0gIhqNOY0UphRGPjSlYeAcHWGsq53RrnZGO9qYGRkCKRHCRGZxCflVNeRX1ZJfVUtaTuxNct5AiNtDc1zrc3Ktb4YbA7N4A2q7qirHxrEyB8f22DlW5iArJXFL3ofVtiIdHR1MTKgtuIyMjIholJSUaFsRjeYpsxURRRXwi0Apq7q5pZQvbNEat5SdKBQPEzJCdM11KdGYuMGNiRvM+mYByE7OpjGnkWN5xziWd4wCW8EDz/V53Ix1dTDa2cZoZztjXR34l1SOxJZhp6C2nsLaeor21mMvKIpZOPxBgzvDSjiu9irh8PhVVdeeLCvH9zg4scfByXIHDtvWCMf8/Hwk0ujr69O2IhrNNrEVQnEb+D3gBhCpB5VS3tjkguzAn6OEpx/4USnl7BrnhYC74buDsZbjPgtC8TCGNOid61URRzjqmPHOAFBoK+RY3jGO5x3naN5R7JYHK4mkYTAzPMhIRxvDbfcYbruHy6mem5SSSmFtPYW1dRTubSCzuASTKbaEciBkcG9kXkUcvTNc759l0RcEoDYvlaZyB00VmRwts2NNfPwIwO/309fXFxGOxcVwxKVtRTSaJ86WlMdKKQ9v4YL+HeCUUv6aEOKfARlSyn+6xnkuKeUjf518FoXiYaSU9Mz1cG38GlfHrnJ9/DqugAuAqoyqiHAczjmM1Wz90HPnJycYvn+X4bZWhtvuRqqsEq1WCqr3KvHYW09OWUXM1VXBkMHdkXku98xwqXua6wOz+IMG8SbBweJ0TpZn0lSRyYGi9MeuqlrPViQ1NZXKykqqq6spKyvDbN7aJLxG87yyFULxK8Ak8N+BSNuwlNK5yQV1AGellGNCiDxUf0b1Guc9t0LxMEEjyP2Z+1wbU8Jxa/IWfsNPvIinPrM+sk21P2s/CXEfTkQvTE8x0naPobZ7DLe1Mjs6DEBCUjKFe+sprttPccN+MguLY66s8gZCXO+f5VLPNJe7p7kzMo+UkJwQx9EyO6cqMjlZnklNbspjNwEuLi5GejZW24rs2bMnskWlbUU0ms2zFULRt8ZhKaXcs8kFzUkp01fdn5VSZqxxXhC4BQSBX5NS/o8o13wdeB2guLj48MDAwGaW9szgDXq5NXWLa2PXuDZ2jdaZVgxpYImzcDjnMCfyT3Ay/yQV6RVrbtW452YZun+XodY7DN67zdz4GABJqWkU1+2juH4/xfX7Hyk5Pu8JcKV3hss901zqnqZnSjUb2q0JnCh3cLoik9NVWRQ8ZgPgsq1IR0cHnZ2d2lZEo9kCtqXqSQjxD0DuGg/9EvDHMQpFvpRyVAixB/g+8FEpZc9Gr70bI4qNWPAvcH38OtfGrnFl7Ap980rbs5KyOJF/gqb8Jo7nH/9QfiPy/OlJBu/dYejebQbv3cY1q4LFlMwsiuv3U1K/n6K6fdjsjpjXND7v5VL3dDjimGF8QVVTl2dZOV2ZxZmqTI7vcZCcsPn8hpSSqampyBbV0NDQh2xFysvLSUzcmuS7RrNb2bRQCCG+KaX8d+Hff0RK+Z1Vj/2fUsp/sckFxbT19NBz/gj4WynlX2x0/edRKB5mzDXGlbErXB69zNWxq8z71CiRWnstJ/NPcjL/JAeyD6y5TSWlZHZshMG7txlsvc1Q6128LpVYthcUUVyvIo6ivfuwxFiRJKWke9LF+a5pzndOca1vBm/AwBwnaCyxc6Yqi9OVmezNS32sbSptK6LRbI7HEYqbUspDD/++1v1HXNCvAzOrktl2KeU3HzonA/BIKX1CiEzgCvAZKeX9ja6vheJBQkaINmcbl0cvc3n0MrcnbxOUQZLik2jMaaSpoIkT+ScoSy1bc5tJGgaTA30M3rvN0L3bDLe1EvB5EcJEbnklJfsPUtJwgLzKmphMD0HlN24MzHK+c4rzXdO0jS0AkGlL4FRFJqcrlXBkp27eYHA9W5GsrKyIaBQWFmpbEY2GxxOK1d5Okd/Xuv+IC3IA/w0oBgaBH5FSOoUQjcAbUsovCyFOAv8fyu3bBPx7KeUfxHJ9LRTRcQfcvDf2HpdHL3Nl7AoDCyqfk2vN5WT+SU7kn+B47nHSLelrPj8UDDDW3cnAnVsM3P2A8a5OpDRISEqiqG4fJQ0HKNl3iIy8/JjzG5OLXi6Go42L3dNMu/wA1OSmRKKNI6V2LObNf6jPzMxE5mysthVZ3qKqqKggKenJGChqNDudHRdRPGm0UDwaw4vDXBm7wpXRK1wdvcpiYBGBoM5Rx8kCtU21L2sfZtPapahet4uhe3cYuPsB/Xc+YH5iHIDUrOyIaBQ37CfJlhLTegxD0ja+wPnOaS50TXG9fxZ/yCAx3sSxPQ7OVGZypiqLymzbpnsqvF4vPT09dHZ20tXVpW1FNM89jyMUIdToUwEkAZ7lhwCLlHJHFrFrodg8QSNI60wrl0fUNtXd6buEZAib2cbxvOOcKjhFU0ETuda16hQUc+NjSjRuf8BQ6x18HjcIQe6eCkr2HaJ030HyqqqJi4/tfx+PP8i1XictnVNc6JqKVFPlplporsribHUWTZWZpFo297+jYRiMjIxEtqi0rYjmeUR7PWk2zYJ/gffG3uPiyEUujlxkwqM+RCvSKyKicSj70JpJcQAjFGK8p5P+2x8wcPcWY13tSMPAnGihqK6Bkn0HKdl3EHt+Yczf3kfmlrjQOcX5rikudE2z6A0SZxIcLs6guTqL5qos6vJTNx0NzM3NRbaoent7CYVCEVuR6upqKioqtK2IZtehhUKzJSx3i18cucjF0YvcnLhJwAiQFJ/EsdxjEeEoTClc9xo+j5vB1jsqv3HnZqR/w+bIpHTfIcoOHKK44QAWa2wfxIGQwQeDc7R0TnKuY4rWUZUUz0pJpLlKicbpykzSkzfnhruerUhhYWEk2sjJydFbVJpn6nO/ygAAIABJREFUHi0UmieCJ+DhvfGVaGPEpWw2SlNLI6LRmNOIJX79yqX5yXEG7tyi/85NBu/exudxI0wm8iprKNt/iNIDh8kpK4+5W3xy0cv5zmnOdUxyoWua+aUAJgEHitI5W53N2eos6vPTNlWCK6VkbGwsIhqjo6OAshVZFg1tK6J5VtFCoXniSCkZWBjg0uglLoxc4Pr4dXwhH4lxiTTmNnIq/xSnCk5Rklqy7rdvIxRirKuD/ts36Lt1k4m+bpCSpJRUSvcfonT/IUr2HcSa/qHezDUJGZJbQ3O0dEzS0jkVsRhxWBM4E442zlRlbXr2xuLiYmSLqqenh0Ag8ICtSFVVFampqZu6tkbztNFCoXnqeINebkzciEQbyxP/CmwFnCpQonE09yjJ5vUHL3kW5hm48wH9t27Qf+cDPPPKqiO7rJyyA4cp3X/okXo3Zlw+LnSpaON81zROtx8hYF9heiQpvr8wnbhNRBvBYJD+/v5ItPGwrUh1dTV5eXnaVkSzY9FCodl2hheHuTRyiYujF7k2do2l4BJmk5lDOYci0UZ5evm60YY0DCb7e+m/fZO+WzcY7WxDGgYJScmUNByg9ICKOFIzs2Naj2FI7o7Mc65jiv+/vTMPjiu77vN3sRIgQewNYiFWgiBIgsRKDheA5Iw8M5rYViaWU5IdR3JsT2RbsrI45aim4kysKJJjV9nyFkmZUpWdOLa8RLbsSJal4QKQQ84ABECCCwCC2NgAiB3EDvRy88d9/djEAI1uDgAC4PmqWOx+y32nHx/79L3nnN+51D5E84MJvBoSYiOpLkzlrDXbeJqGTYFkRXwzjfz8fJEVETYV4iiETcWiZ5GmoSZ7ttEx0QFAWmyaPds4nn6cuKiV6y4WZmfobblB143rdDc3MjU6DBiJkbzScnKPVpBVfJiIqOCWlcZnFrncMWI5jmFGpo1I8uHM3dZsw0HZ3gQiwkOfEfjLity7d4+FhQXCw8OfkBVJTAxuOU0Q1gtxFMKm5uHMQ670XeFK/xWu9l9l2jVNhIqgLK2M6sxqTmeeXlEFF8wv+LE+pxXbuI7z7i08LhcRUdHsPXiY3NIKco9WBF0p7vVq7gxMcql9mIttQzT2TuDxauJ2RFBTaJaozhSl4ogLXV7E4/HQ29trzzZGR02DKZEVEZ414iiELYPL6+Lm8E0u912mzllH23gbYORFqjOrqc6s5nj68YCxDdfCPM47t+zZxviAycaKd6SRe7SC3NIKsg+VEBWz8hj+PJpzcaXDxDYutg0zNGVmGyWZ8ZwrSuXcAQdHnjK2MTo6ajuNnp4evF4vMTEx7Nu3T2RFhA1FHIWwZfHNNur66rjaf5VZ9yyRYZFUplVSnWUcR6BMKoCJwYd032ik+0Yjvbdu4JqfIyw8gqzig+SWVpJfVhl0b3GtzWzjYtswF1qHaOwdx6tNzw1fQLymMJXEp8ikWklWJCcnx55tJCcnS82GsC6IoxC2BS6Pi8ahRuqcdVzuu8z9R6Y1yd64vZzOPE11ZjVVe6oC1m143C762+7S1WyWqUZ6uwHTdyO/rJLc0kqyDx8hakdwv+InZhepvTfCxdYhLrYPMzazSJiCsuxEzhWZ2MbTVImvJCuSlJRkO43s7GyRFRHWDHEUwrakb7qPy87L1PXV8d7Ae8x75okOj+bYnmP2bCNQlTiYFrHdN67T1dRAT4uZbYRHRJBZfNhyHBVBy4t4vJqbzgkutJnYxk2n6QHiiIvmbFEq54ocT61JNTExYTuNrq4uW1bEt0RVWFjIzp07Vx9IEFZAHIWw7VnwLNDwsIG6vjpqnbU8mHoAQF58nh0Qr0irWFGTCsxso6/1jpltNDUw6uwFrNiGtUS191AJkdHBBbGHpxa41D7MhbYhatuHmZp3ExGmqMxN5FyRgxcPONj3FAq4i4uLdHZ22o5jenoaEFkR4cMhjkJ47uiZ7KHOWWdXiS96F4mJiOGF9Bfs2UYgBVyAyeEhupob6Gq+Tk9LM+6FBcIjI9l7sIS80gryyipJTM8Myh63x0tj7wQX2oa40DpE60OjGZWZEMO5A2a2caIg9LawXq+Xhw8fiqyI8KHZdI5CKfWTwFtAMXBMa73st7pS6lXgq0A48LbW+ivBjC+OQvBn1jVL/cN6e7YxMGOECAsTC+1MqqOOoyv22wBwu1z03b1NV3M9XU3XGet3ApCQlk5eWSV5pRVkHSohMiq4Irr+iTkTEG8b4krHCLOLHqIiwnghP9lkUhU5yE0JfSlpOVmRyMhIW1aksLBQZEWEZdmMjqIY07nu68CvLucolFLhQDvwI4ATqAc+Ka1QhQ+D1prOR532bKNxsBG3dhMXGceJjBNUZ5llqpSYlIDjTAw+pKu5ge7m6/Teuol7cYGIyCj2HiqxHEclCXvSg7Jpwe2hvmvczDbahui0+m3kpey0YxvH8kLv7udyuejp6aG9vZ22tjYePTIxk/T0dHu2IbIigo9N5yjsiyt1kZUdxQngLa31K9b7LwBorb+82rjiKIRgmV6c5trANbtuY2huCIDipGJ7iaokpYTwsJW/pN2LizjvtFiZVA2MD5jln8T0TLNEVVpB1sGSoKvEe0Zn7NnG1fujLLi9xESGc2pfMmeLHJw74CAzIbTaCq01Q0ND9hKV0+lEa82uXbvsVrAiK/J8s1UdxceBV7XWP2+9/xnguNb6s6uNK45CeBq01rSPt1PXV0eds47m4Wa82kt8dDynMk5RnVXNqYxTJO4ILLcx/rCfrqbrdDc38OB2C27XIhHR0WQfOkJeaSV5ZRXEOwLHR3zMLXq41jnKhbYhzrcO4RyfA6AoLY6zVmyjIieRyBClRWZmZmxZkY6ODpEVEZ6No1BK/RBY7n/Dm1rrv7WOucjKjuIngVeWOIpjWuvPrXC9N4A3ALKzsyt6enrW5HMIzy+PFh5xtf8qdX2mbmNsfgyFoiS1hNOZp6nJrKE4uZgwtfKXtGthngd3WuhqMrMNXz/xpIwse4kqs/gQEUEEm7XW3B+e5kKrmW283zWG26uJi46gen+K6bexPxXH7tCkRVaTFSkqKiIrK0uWqLY5W3VGIUtPwqbBq73cGb1jF/u1jLSg0STvSOZU5ilqsmo4mXEyoJCh1prxgX66mxvobGqwNakio3eQXXLUWqaqZHdqcAq4U/MurnSMctGKbQxOPhYyPFfk4GyRg9K9oUuLjIyM2AFxf1kR3xJVQUGByIpsQ7aqo4jABLNfAvowweyf0lrfXm1ccRTCejM2P2akRZx1XOm/wuTipC1kWJNZQ3VWNfnx+QFrGVzz8/TevmnXbUwOm+rr5Kxs8spM3UZG0cGg+m1orbk7MGWn3y6VFjl3wMGZwlTiY0NLk/XJirS1tXHv3j3m5uY+ICuSkhI48C9sDTado1BKvQ78PpAKTADNWutXlFIZmDTY16zjXgN+F5Me+02t9ZeCGV8chbCRuL1ubg7fpNZZS21fLffG7wGmSVN1ZjU1WTWrSov4FHC7mhvoamrAefc2Xo+bqJhYco+UkVdeRV5pRdDd/SZmF02xX6vp7jc+6yI8TFGRnci5A6bYb39aaMV+Xq8Xp9NpL1ENDZnAv8iKbA82naNYb8RRCM+SgekBOyB+beAa8555doTv4Hj6cWqyaqjOrCZ9V+DU2cW5WXpu3aCrsZ6upgamx8cASMsvNLON8kr25BcG1UvctIQd53zrEOdbh7k7MAk8LvZ78YCDkwUpIaffLicrEh0dTUFBgciKbEHEUQjCM2LBs0D9w3oz23DW0jdtJM/3JeyjJquGmqwajqYeJSJs5V/hWmuGe7rotJzGwL02tPYSszuevNIK8suryDlSxo6du4KyaeDRHBdahznfaor95lweoiPCOFmQzIsHTPptVmJwEuz251xYoKur6wOyInv37rVnGw6HQ2RFNjHiKARhE6C1putRF7XO2ieK/XZH7bbTb09nnl41/XZuapLu5ut0NjXQfaOR+ekpVFgYmUUH7dhG8t7A0us+5l0e3usa40KrSb/tHZsFYH/aLrNEZaXfhtLZz19WpK2tjYEBUwkfHx9vO43c3FyRFdlkiKMQhE3I1OIUV/uv2o7DP/22JtPMNg4kHQj4he/1ehi4105XUz2dTQ0Md3cCj2XT88qqyD58JCghQ5N+O2M7jfpuk367e0cENftTeanYwZn9DpJC7LUxOTlpZ1F1dnaKrMgmRRyFIGxyfOm3tc5a6px13Bq9BYAjxmEqxLOqOZF+ImBnP4CpsRG6mkxAvOdmM66FeSNkaBX75ZdXkZAWXLHf5LyLy/dGON86xMW2IUamF1EKyvYm2EtUB9ND67Xhcrno7u62l6hEVmTzII5CELYYI3MjXO67TK2z1u4j7t/ZryarhpzdOQHH8AkZdjbV09VUb0uL+Ir98suryDxwkPCI1ZeAvF5NS98jzreamg1fr409u3fY6ren9qWwMzr4jKelsiIPHhhp+F27dtlOIz8/n6ggpU+ED4c4CkHYwri8LpqHmu2AeOcjs7yUszvHTr+tTKskMjzwF/74QB9dTVax350WPG43UTEx5JSUmSrxskp2JSYFZdPQ1LzdDrbu3gjTC26iwsM4np/Ei1b6bU5yaBlPK8mK5OXl2Y4jISEhpDGF4BFHIQjbCOeU067ZqB+oZ9G7SGxErFG/zTTLVI7YwNXdi/Nz9N66SVdjPZ3NDUyPjgDgyC0gv9zENvbsKyQsgBiiPZbbS0P3mEm/9VO/zU/dyYtWg6bK3CSiIoJfTvLJirS1tdHe3s7YmEkPdjgcttMQWZG1RRyFIGxTfL02fI7j4YzRkvKp39Zk1XA4+XBA9VutNSO93Sb9trmB/rZWk34bt5vc0gryyyrJOVpOzK6V5Un86RmdsWo2hnivc4xFj5dd0RFUF6Zw7oCDs0WpOOJC06MaGRmxl6h6e3s/ICuyb98+duwIbUzhScRRCMJzgNaaexP37IC4T/02MTrxCT2q+Oj4gOPMTU/RfaPRBMWbrzM/NYlSYWQUHbAD4inZuUEFsWcW3FzpGLHVb316VEey4u12sCWZ8YSFoEc1NzfH/fv3aW9vt2VFwsLCyM7OFlmRD4E4CkF4Dnm08Ih3+9+l1lnL5b7LTCxMEK7COZp61C7225ewb9X024cd9+z026Gu+wDsSk4hv9TENbJLjhK1Y3WRQK01dwYm7fTbpgcTaA0pu6I5W2QqxKsLU4jbEXx9RSBZkaKiIltWJDw8tKrz5xFxFILwnOPxemgZabFrNlrHWgFI35luB8SPpR8jJiLwF/702KgtYth9swnX/BzhERFkHSwxdRvlVSTuyQjKprGZRS61G1mRS21DTM67iQhTVOUm2em3Bak7Q0q/HR8ft2s2/GVF9u3bZy9RiazI8oijEAThCQZnBh+n3w5cZc49R3R4NFV7quzZRuauzIBjeNwu+lrv0NloZhvjVh/xxPRMq0K8KuheG26Pl8beCZN+2zpE2+AUANlJsbbTOB5iO9iFhQU6OzvtJarp6WmUUmRlZYmsyDKIoxAEYUUWPYs0DDZQ56yj1llL71QvAAXxBXZAvNRRSmRY4C/8iYcDdDY10NXcwIPbN02vjR0x5JQctdNv45KCix04x2e5YKXfXukY8WsHm2Kn3+6JDz547fV6GRgYsJeoRFbkg4ijEAQhaHome+yajYbBBtxeN3GRcZzIOEFNVg2nM0+THJMccAzTa+OGqdtobGBqdBiA1Jw88suryCurIr1wf1Dpt/MuD1fvj9qZVH0Tph1scfpuXrTUb0v3JobUoMlfVuT+/fu43e4nZEX2799PXFxwWV7bBXEUgiA8FTOuGa71X6O2z2RSDc8No1AcTjlsxzZWawertWb0QY+ZbTQ10Nd2B+31smNXHLlHy8kvqyS3tIKYuNX1nrTW3Buatp3G9Z5xPF5NYmzk4wZN+1NJiA2+mttfVqStrY3JSSPDnpGRYTuNPXv2bPuaDXEUgiB8aLTWtI612jUbLcOmHWxKTIrpIZ5Vw4n0E+yKCix3Pj89TffNx+m3c5OPUCqM9MIi8surQkq/fTTrovaeWaK62D7M2MwiYQoqchJ58UBayA2a/GVF2tracDpN3OV5kBXZdI5CKfWTwFtAMXBsuVao1nHdwBTgAdwrfYiliKMQhPXH1w621lnLlf4rTC1OEaEiqEiroDqrmjNZZ8iNzw04hvZ6eXj/Hp1NDXQ2vm+n38Ylp5JXVkF++bGg1W89Xs0N54Sdfnu7/4MNmk7kpxATFXxAfGZm5oklqu0sK7IZHUUx4AW+zgo9s63juoFKrfVIKOOLoxCEjcXtdXNj+IYd2+iY6AAgOy7bzqIKRo/Kl37b2VhPz80mXAvzRERGsfdQCXnlVeSXVRHvSAvKpsHJedtpXO4YYXbRr0FTsZltZCasXv9hf0a3m97eXjsgvt1kRTado7AvrtRFxFEIwrajf7rfdhrvP3yfBc+CrUd1JusM1VnVpMQEzoByu1w4794yelSN9UwMmkyl5Kxss0RVVkVGUTFhQRTTLbg9vN81Zsc2ekZNg6YDe+I4d8DBSwcclGWHFhD3lxXp6elBa01sbKwtK1JQULClZEW2sqPoAsYBDXxda/2NAGO9AbwBkJ2dXdHT07P2BguCEDJz7jneH3ifS85L1DprGZwdBOBQ8iF7tnEw+eCqAfHxgT6rHWw9zru38Xo8RO/cSe5R0w4292g5sbsDy5P4xuocmeH83ScbNCXERnL2KQPiK8mK5OTk2LON5OTAmWLPmmfiKJRSPwSW65Dyptb6b61jLhLYUWRorfuVUg7gB8DntNa1q11bZhSCsDnRWtM+3m7PNm4M30CjSd6RbMc1TmScYGdk4OrphdlZelqa6LxuhAxnH02AUiYgXmYC4qk5eUEFsSfnXdS1j/BO6yAX20xAPDxMUZGdyIvFpmaj0BF8QNzj8TwhKzI8bFKDk5OTbaexGWVFtuyMYsmxbwHTWuvfXu1YcRSCsDUYnx/nct9l6px1XO6/bALiYSYgfibrTFANmrTXy2BnB51NZolqsNPER2w9qvIqcg4fJTKIZSD/gPg7d4e4M/A4IP5SsakQP5GfHFKF+Pj4uO00uru7PyArUlhYSGxs4M6FG8GWdBRKqZ1AmNZ6ynr9A+A3tNb/sNq44igEYevh9rpNg6a+Wmof1HL/kcmAytmdYy9RVTgqVg+Ij4/R1dxAV6OfHpXVDjbf6uwX7wiuHezAozkutA5z3qoQn3N5rArxZDv9NpQKcX9Zkfb2dmZmZmxZEZ+IYWpq6jORFdl0jkIp9Trw+0AqMAE0a61fUUplAG9rrV9TSuUD37ZOiQD+j9b6S8GML45CELY+yzVo2hm5k5MZJ+0GTcEExJdrB5uclW23g83YX0x4xOotXOddHq51Pq4Qd46bCvGD6buNrEixg6NZCUEHxFeSFUlISLCXqHJycjZMVmTTOYr1RhyFIGwvZl2zvDfwnpltOGsZmjVy4oeTD5vZxt4aipMCV4gDVkDc1GyYgLib6NidpkK8vIrc0oqgA+J2hfjdIa73mgrxpJ1RfpLpqcTHBP8lPzk5aTuNzs5OW1akoKDAXqJaT1kRcRSCIGwbtNa0jbdR66zlkvPSExXiNVk11GTW8ELGC0EFxHtbmu3Yhh0Q37ef/PJjIQXEJ2YXudT+uEJ8YtZFRJiiMjfREjFMC0ky3eVy0dXVZTuO5WRF0tPT13SJShyFIAjbFl+F+CXnJd7te5cp1xSRYZFUplXasY3s3dkBx9BeL4Nd9+3024f37wGwKynZWqI6FlJAvKl3nHcsyfTWh09Kpr94wMHx/CSiI4ILiGutGRwctJ2GT1YkLi7OrtlYC1kRcRSCIDwXuLwuExC30m87H3UCkLs7l5qsGs5knaEsrWxVyfSZiXGjfNtkKsQX56yA+MES8qz024S04ALifRNzdp8Nn2R6bFQ4p/elmEyqIgeO3cEHxKenp+no6KC9vZ2Ojg4WFxeJiIiwZUXKy8ufKvVWHIUgCM8lD6Ye2E6j/mE9Lq+LXZG77ArxYCTTHzdoep/OxgbGB/oASMrIIr/iGPlllWQUHQwqID636OFq5wjv3DWOo//RPACHM3fbWVRHQugh7i8r0tbWhtaaz3/+80+1JCWOQhCE555Z1yzXBq6ZdrDOOobmhmzJdN8SVXFS8apfsuMP++lqrOd+Yz3OO7fsgHjOkTKr10Zl0AHx1odTdhZVU+84Xr8e4i8dcHA6hB7iWmtmZ2efutWrOApBEAQ/fJLpl5yXqHPW0TJiAuKOGAfVWSb19kT6CWIjAxfCLc7N0tPSTGdjA11N9cxMjJuAeMF+8spNbMORmx/UL/zleohHhiuO5SVxrsjBS8Vp5KWsX79vcRSCIAgBGJ0btXuIv9v/LtOuaSLDIp/oIb43bm/AMbTXy1B3p9VD3AqIa82uxCTTCra8ipySUqJ2rK5Y6/Z4ud4zbs827g1NA5CXstMOiFflJhEVsXZKteIoBEEQgsTlddE02GSn33ZPdgOQF59ny4oE00N89tGEkUy//j7dN5tYnJslPCKCrIMltvptwp70oGx6MDbL+dYh3mkd4tr9URY9XnZFR1BdmMK5AyYgnhoX/aE+tzgKQRCEp6R3sveJHuIur4u4yDhOZp60e4gn7UgKOIbH7TYBcatmY7zfpLgmZWTZfTYyDwQXEJ9ddHOlY5TzrYOcbx1icHIBgKNZ8bx4II1fPldARHjoMw1xFIIgCGuAfw/xWmctI3MjKBQlqSXUZNZwZu8ZihKLVo1JTDwcsJ2G804LHrebqJhYu0I8r7SC2PjVO+dprbndP2kvUU3Nu3jn3599qs8mjkIQBGGN8Wovd8fumtnGg1pujd4CwBHroDrTSKYfTz++ekB8fo6elma6GuvpampgenwMlGJPQaG9ROXIKwgqID7v8oSkbOuPOApBEIR1ZmRu5ImA+IxrhqiwKKrSq6jJNAHxrLisgGNorRnq7rS7+g3cbwet2ZmYRF5pJfnllSYgHrP2suTiKARBEDYQl8dF41CjnX7rC4gXxBfYWVSljlIiwgLHJOyAeFMD3c3XnwyIW5lUiXsy1sRmcRSCIAjPkJ7JnicC4m6vm7ioOE5lnLID4ok7EgOO4XG76W+7Q2dTA52N9Yz1PQAgMT3TLFGV+wLiTydLLo5CEARhkzC9OM21gWv2bGN0fhSF4kjqETv9dn/i/tUD4oMPbRHDB7dv4nG7iY1P4F//jz8mbDtoPSmlfgv4MWARuA/8rNZ6YpnjXgW+CoRjGhp9JZjxxVEIgrAV8Govd0fvcsl5iVpnLbdHbwOQFptmL1EdTz9OTETgIj3X/Dw9t24wOTxI+Ud//Kls2YyO4mXgvNbarZT6TQCt9a8tOSYcaAd+BHAC9cAntdZ3VhtfHIUgCFuR4dnhJwLis+5ZosOjn6gQz9yVuS7X3nSO4gkDTFvUj2utf3rJ9hPAW1rrV6z3XwDQWn95tTHFUQiCsNVZ9CxyffC6HdvoneoFYF/CPttpHE09umpAPFg2u6P4O+BbWuv/vWT7x4FXtdY/b73/GeC41vqzK4zzBvAGQHZ2dkVPT8/6Gi4IgrCBdD/qtuMa1wev49Zudkft5lSmFRDPOE3CjtWL9FYikKNYG1e0/EV/CCzX2eNNrfXfWse8CbiBP11uiGW2rejVtNbfAL4BZkYRssGCIAibmNz4XHLjc/nUoU8xtTjF1f6rRjK9r47vdX2PMBVGmaOMt19+e81mGT7WzVForT8SaL9S6lPAjwIv6eWnNU7AX64xC+hfOwsFQRC2JnFRcbyc+zIv576MV3u5PXKb2r5ahmeH19xJwDo6ikBY2Uy/BpzRWs+ucFg9UKiUygP6gE8AP7VBJgqCIGwJwlQYJakllKSWrN811m3kwPwBEAf8QCnVrJT6GoBSKkMp9V0ArbUb+CzwfeAu8Bda69vPyF5BEITnlmcyo9Ba71thez/wmt/77wLf3Si7BEEQhA/yrGYUgiAIwhZBHIUgCIIQEHEUgiAIQkDEUQiCIAgBEUchCIIgBEQchSAIghCQZ671tB4opYaBpxV7SgFG1tCctULsCg2xKzTErtDYjnblaK1Tl9uxLR3Fh0Ep1bCSMNazROwKDbErNMSu0Hje7JKlJ0EQBCEg4igEQRCEgIij+CDfeNYGrIDYFRpiV2iIXaHxXNklMQpBEAQhIDKjEARBEAIijkIQBEEIyHPjKJRSryql2pRSHUqp/7jM/mil1Les/e8ppXL99n3B2t6mlHplg+36d0qpO0qpm0qpd5RSOX77PFY/j2al1Hc22K5PK6WG/a7/8377PqWUumf9+dQG2/U7fja1K6Um/Pat5/36plJqSCl1a4X9Sin1e5bdN5VS5X771vN+rWbXT1v23FRKvauUOuq3r1sp1WLdr4YNtuusUuqR37/Xr/vtC/gMrLNd/8HPplvWM5Vk7VvP+7VXKXVBKXVXKXVbKfX5ZY5Zv2dMa73t/wDhwH0gH4gCbgAHlxzzS8DXrNefAL5lvT5oHR8N5FnjhG+gXeeAWOv1L/rsst5PP8P79WngD5Y5NwnotP5OtF4nbpRdS47/HPDN9b5f1tg1QDlwa4X9rwHfw/SCfwF4b73vV5B2nfRdD/iozy7rfTeQ8ozu11ng7z/sM7DWdi059seA8xt0v9KBcut1HNC+zP/JdXvGnpcZxTGgQ2vdqbVeBP4c+NiSYz4G/LH1+q+Al5RSytr+51rrBa11F9BhjbchdmmtL+jH7WKvYXqHrzfB3K+VeAX4gdZ6TGs9DvwAePUZ2fVJ4M/W6NoB0VrXAmMBDvkY8CfacA1IUEqls773a1W7tNbvWteFjXu+grlfK/Fhns21tmsjn68BrXWj9XoK0/Uzc8lh6/aMPS+OIhN44PfeyQdvsn2MNm1YHwHJQZ67nnb583OYXww+diilGpRS15RS/3Qgui1dAAAHY0lEQVSNbArFrp+wprh/pZTaG+K562kX1hJdHnDeb/N63a9gWMn29bxfobL0+dLAPyqlriul3ngG9pxQSt1QSn1PKXXI2rYp7pdSKhbzZfvXfps35H4psyxeBry3ZNe6PWPPpBXqM0Ats21pXvBKxwRz7tMS9NhKqX8BVAJn/DZna637lVL5wHmlVIvW+v4G2fV3wJ9prReUUp/BzMZeDPLc9bTLxyeAv9Jae/y2rdf9CoZn8XwFjVLqHMZRnPbbfMq6Xw5Mf/tW6xf3RtCI0R6aVkq9BvwNUMgmuV+YZacrWmv/2ce63y+l1C6Mc/o3WuvJpbuXOWVNnrHnZUbhBPb6vc8C+lc6RikVAcRjpqDBnLuedqGU+gjwJvDjWusF33Zteoyjte4ELmJ+ZWyIXVrrUT9b/idQEey562mXH59gybLAOt6vYFjJ9vW8X0GhlDoCvA18TGs96tvud7+GgG+zdkuuq6K1ntRaT1uvvwtEKqVS2AT3yyLQ87Uu90spFYlxEn+qtf6/yxyyfs/YegReNtsfzMypE7MU4QuAHVpyzC/zZDD7L6zXh3gymN3J2gWzg7GrDBO8K1yyPRGItl6nAPdYo6BekHal+71+HbimHwfOuiz7Eq3XSRtll3VcESawqDbifvldI5eVg7P/hCcDje+v9/0K0q5sTNzt5JLtO4E4v9fvAq9uoF17fP9+mC/cXuveBfUMrJdd1n7fj8idG3W/rM/+J8DvBjhm3Z6xNbu5m/0PJiOgHfOl+6a17Tcwv9IBdgB/af2neR/I9zv3Teu8NuCjG2zXD4FBoNn68x1r+0mgxfqP0gL83Abb9WXgtnX9C8ABv3P/lXUfO4Cf3Ui7rPdvAV9Zct56368/AwYAF+YX3M8BnwE+Y+1XwB9adrcAlRt0v1az621g3O/5arC251v36ob17/zmBtv1Wb/n6xp+jmy5Z2Cj7LKO+TQmwcX/vPW+X6cxy0U3/f6tXtuoZ0wkPARBEISAPC8xCkEQBOEpEUchCIIgBEQchSAIghAQcRSCIAhCQMRRCIIgCAERRyFsCZRSWin1v/zeRyijXvv3G2jDWaXUSb/3n1FK/csQzvdXr21WSuUqpSqVUr8XxHWX/ZyWYmmK3/uvWyqjzcqoDs/5Xe/jwdpqjfURpdTfhHKOsD15XiQ8hK3PDHBYKRWjtZ4DfgTo22AbzgLTmGIqtNZfC/H8Oa116ZJt3cBaSlIfByq01h5LE+jvl7mmIISEzCiErcT3MNWnsES5Uyl1zOqn0GT9XWRtj1VK/YUlXvgtZXqNVFr7ppVSX7KE564ppdKs7alKqb9WStVbf05ZX7qfAf6t9eu8Win1llLqV61z9imlfmiN1aiUKgjmA/nPFpRSO5Xph1BvfY4PqKIqpZKVUv9o7f86fjo+SqlioF0/qW+19PxCpdT3LeG6WqXUfmv7J5Tpr3BDKXVhmfNeUEpdta57RSlVGMznE7YH4iiErcSfA59QSu0AjvCkemYrUKO1LgN+Hfhv1vZfAsa11keAL/JYkwqM1MI1rfVRoBb4BWv7V4Hf0VpXAT8BvK217ga+Zm0v1VrXLbHtT4E/tMY6ianuXUqM3zLQt5fZ/yamv0EVpg/Jbymldi455j8Dl63P+R2MBIePjwL/sMy4/nwD+CWtdQXwBeAP/MZ9ybL/9WXOuwuctq77ReC/rnIdYRshS0/ClkFrfdP6Zf9J4LtLdscDf2z90tVApLX9NOaLH631LaXUTb9zFgHf2v91zHIWwEeAg0rZP9Z3K6XiVrLL2peptf62dZ35FQ5dbunJn5eBH/fNUjCyMtlLjqkB/pl1nf+nlBr32/cK8LMB7EzAaAD9td9n830HXAH+RCn1l8BygnMJ1v6gZkrC9kIchbDV+A7w25h4QbLf9i8CF7TWr1vO5KK1fTmJZR8u/VjDxsPj/w9hwAkrFmLj9+W6lEDXCAUF/ITWum3JddOWHPcB3R1l+iMkaEvBNMD4Iys4q1/AxDd+FLhhKcr68yXg+1rrP1JK7WP1mYuwjZClJ2Gr8U3gN7TWLUu2x/M4uP1pv+2XgX8OoJQ6CJQEcY1/xIjSYZ3n+2KdwrShfAJt+gI4ldUMSZn+67FBXGcp3wc+pyyPpJRaTga9Fvhpa/9HMWqgYJaqPhBbWGLnODCglHrdOj9MPe6Rna9NV7T/hBEJXNrYZqX7KzwHiKMQthRaa6fW+qvL7PrvwJeVUlcwfZV9/BGQai05/RpGffPRKpf5FaDSCoDfwQSxwTRret0XzF5yzs8Av2Jd512MTHaofBGzZHZTKXXLer+U/wLUKKUaMUtVvdb2YOITYCT0P6OU8qmc/qi1/XeUUi0Y1dEfaq1vLTnvNzExkyuhfCBheyDqscK2RikVDkRqreet9fV3gP3a9FveNliO47jW2vWsbRG2HxKjELY7scAFZbqDKeAXt5uTANBalz9rG4Tti8woBEEQhIBIjEIQBEEIiDgKQRAEISDiKARBEISAiKMQBEEQAiKOQhAEQQjI/wfgwiqZoFLkIQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure('Breit-Rabi_Formula')\n", "plt.plot(Bs, ws[:, 0])\n", "plt.plot(Bs, ws[:, 1])\n", "plt.plot(Bs, ws[:, 2])\n", "plt.plot(Bs, ws[:, 3])\n", "plt.plot(Bs, ws[:, 4])\n", "plt.plot(Bs, ws[:, 5])\n", "plt.plot(Bs, ws[:, 6])\n", "plt.plot(Bs, ws[:, 7])\n", "plt.xlabel('Magnetic Field/Tesla')\n", "plt.ylabel('Energy/GHz')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.4" } }, "nbformat": 4, "nbformat_minor": 2 }