{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import numpy.fft as fft\n", "import matplotlib.pyplot as plt\n", "from lmfit import Model, Parameters\n", "from lmfit.models import ConstantModel\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [], "source": [ "def cavity(x, x_0, I_0, gamma, bg):\n", " return I_0/(1 + (4*(((x-x_0)/gamma)**2))) + bg\n", "\n", "def fit_cavity(freq, trans, trans_err):\n", " mod = Model(cavity)\n", " \n", " #set parameter and make init guess\n", " p = Parameters()\n", " p.add('I_0', 1)\n", " p.add('x_0', 138)\n", " p.add('gamma', 6)\n", " p.add('bg', 0.1)\n", " \n", " return mod.fit(trans, x=freq ,params=p, weights = 1/(trans_err),method = 'least_squares')" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [], "source": [ "filename = 'cavity_laser_scan.dat'\n", "data = np.genfromtxt(filename)\n", "f_data = data[:,0]\n", "trans = data[:,1]/np.max(data[:,1])\n", "trans_err = data[:,2]/np.max(data[:,1])\n", "\n", "fit_result = fit_cavity(f_data,trans,trans_err)" ] }, { "cell_type": "code", "execution_count": 24, "metadata": { "scrolled": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[Model]]\n", " Model(cavity)\n", "[[Fit Statistics]]\n", " # fitting method = least_squares\n", " # function evals = 6\n", " # data points = 42\n", " # variables = 4\n", " chi-square = 60.3012410\n", " reduced chi-square = 1.58687476\n", " Akaike info crit = 23.1906888\n", " Bayesian info crit = 30.1413673\n", "[[Variables]]\n", " I_0: 1.02519988 +/- 0.01546722 (1.51%) (init = 1)\n", " x_0: 137.764780 +/- 0.03880287 (0.03%) (init = 138)\n", " gamma: 3.92481417 +/- 0.05064578 (1.29%) (init = 6)\n", " bg: 3.7073e-04 +/- 3.8460e-04 (103.74%) (init = 0.1)\n", "[[Correlations]] (unreported correlations are < 0.100)\n", " C(I_0, gamma) = -0.678\n", " C(gamma, bg) = -0.664\n", " C(I_0, x_0) = -0.237\n", " C(x_0, gamma) = 0.171\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAt8AAALJCAYAAABshU9CAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8GearUAAAgAElEQVR4nOzde5xdZX33/c8vk9NADkNCwEwOJgKmoEhGo6hRoYoN9oCRWoMFxVbF+rS192ONQttHvbUKFq2td+0teGjxgGIpjai0KYoolYiEBgmHhoRTyCRAMiHMABMymbmeP/beYWeYSeaw91r78Hm/Xnnt2ddea69fwmbyzTXX+l2RUkKSJElS9U3IuwBJkiSpWRi+JUmSpIwYviVJkqSMGL4lSZKkjBi+JUmSpIwYviVJkqSMGL4lqcFExKKISBExcQTHvisi/mvQ2LcjYmUV63ttRGwa47lTIuJ/ImJOpeuSpCwYviUpZxHxYETsi4ijB41vKIboRRnW8hLgFOB71bpGSummlNKSsms+GBFnjPDcZ4CvARdWqz5JqibDtyTVhgeAt5eeRMTJwBE51PE+4FuptndguxI4PyKm5F2IJI2W4VuSasM3gHeWPT8f+HrpSUTMjIivR8TOiHgoIv4qIiYUX2uJiM9GxK6IuB/4rfI3Lp771YjYERGdEfHXEdEyTB1vAn466Pz3RsQ9EdETEXdHxEuL4xdGxH1l428pjk+JiD0R8eKy95gTEb0RcUxEnB4R24rj3wAWAt+PiCcj4sMR8cOI+NNBNdxRev+U0jbgceCVI/3DlaRaYfiWpNrwC2BGRJxYDMbnAN8se/3/ADOBFwCnUQjqf1B87b3AbwMdwDLgrYPe+5+B/cDxxWN+A3jP4AIi4khgMbCpbOz3gI8XrzcDOAvoKr58H/DaYl3/G/hmRMwtLg25hrKZfOBtwE9TSo+VXzOl9A5gK/A7KaVpKaW/Aa4Aziur4RRgHvDDslPvobA8RpLqiuFbkmpHafb7jRTCZWdxvBTGL0op9aSUHgQ+B7yj+PrbgL9LKT2cUtoNXFx6w4g4FvhN4H+llJ4qht/PF99vsLbiY0/Z2HuAv0kp3ZoKtqSUHgJIKf1LSml7SmkgpXQVsBl4RfG8Kwdd4/eLYyNxLfDCiDih+PwdwFUppX1lx/SU1StJdeOwd8JLkjLzDeBnFGafv142fjQwCXiobOwhCrPBAO3Aw4NeK3l+8dwdEVEamzDo+JI9xcfpwN7i1wsozHA/R0S8E/ggsKg4NK1YK8BPgCMi4lTgUWAp8G9Dvc9gKaW9EXEVcF5E/G8KM+iDZ/Onl9UrSXXD8C1JNSKl9FBEPEBhpvrdZS/tAvooBOm7i2MLeXZmfAeFkEzZayUPA88AR6eU9h/m+k9FxH3AC4GdZecfN/jYiHg+8GXgDcC6lFJ/RNwORPG9+iPiuxSC86PAD1JKPYPfp3TpIcauoPCPkf8Cnk4prRv0+okUZv8lqa647ESSasu7gdenlJ4qG+sHvgt8KiKmF4PvB3l2Tfh3gQ9ExPyIOIqyNnwppR3AfwKfi4gZETEhIo6LiNOGuf51FNaUl3wF+FBEvCwKji9e/0gKoXknQET8AfDiQe91JbAKOJdDLzl5lMJa9gOKYXuAQsD+RvlrETEPmEVhnbwk1RXDtyTVkJTSfSml9UO89KfAU8D9FGaDr6TQ7xoKM9BrgV8B/03hZsdy7wQmU5g1fxy4Gpg7TAmXA+dGcY1KSulfgE8Vr9cDrAFmpZTuphCM11EIzycDPx/0e7mlWHM78O+H+G1fDPxVsUPKh8rGv158328OOv73gSuKN3ZKUl2J2m7lKknKWkRcCXw3pbQm5zreCVyQUnpN2dgUCv/IeN3gzimSVA8M35KkmhMRRwA3AP+YUvr64Y6XpHrhshNJUk2JiBUU1pI/ysjbE0pSXXDmW5IkScqIM9+SJElSRpqqz/fRRx+dFi1alHcZkiRJamC33XbbrpTSnKFea6rwvWjRItavH6qDlyRJklQZEfHQcK+57ESSJEnKiOFbkiRJyojhW5IkScqI4VuSJEnKiOFbkurcqsvWseqydXmXIUkaAcO3JEmSlJFcw3dEnBkRmyJiS0RcOMTrUyLiquLrt0TEouL4oojojYjbi7++lHXtkiRJ0mjl1uc7IlqALwJvBLYBt0bEtSmlu8sOezfweErp+Ig4B/gMsKr42n0ppaWZFi1JkiSNQ54z368AtqSU7k8p7QO+A7x50DFvBq4ofn018IaIiAxrlCRJkiomz/A9D3i47Pm24tiQx6SU9gNPALOLry2OiA0R8dOIeO1wF4mICyJifUSs37lzZ+WqlyRJkkapXm+43AEsTCl1AB8EroyIGUMdmFK6PKW0LKW0bM6cOZkWKUmSJJXLM3x3AgvKns8vjg15TERMBGYCXSmlZ1JKXQAppduA+4AXVr1iSZIkaRzyDN+3AidExOKImAycA1w76JhrgfOLX78VuCGllCJiTvGGTSLiBcAJwP0Z1S1JkiSNSW7dTlJK+yPiT4C1QAvwtZTSXRHxCWB9Sula4KvANyJiC7CbQkAHeB3wiYjoAwaAP0op7c7+dyFJkiSNXG7hGyCldB1w3aCxj5Z9vRf4vSHO+1fgX6teoCRJklRB9XrDpSRJklR3DN+SJElSRgzfkiRJUkYM35JUx9Zs6GTD1j3c8sBull9yA2s2DO7YKkmqJYZvSapTazZ0ctE1G9nXPwBA555eLrpmowFckmqY4VuS6tSlazfR29d/0FhvXz+Xrt2UU0WSpMMxfEtSndq+p3dU45Kk/Bm+JalOtbe1jmpckpQ/w7ck1anVK5bQOqnloLHWSS2sXrEkp4okSYeT6w6XkqSxW9kxD4APX30H+/oHmNfWyuoVSw6Mj8Sqy9YBcNX7XlWVGiVJBzN8S1IdW9kxj2//citggJakeuCyE0mSJCkjhm9JalJu0CNJ2TN8S1ITcoMeScqH4VuSmpAb9EhSPgzfktSE3KBHkvJhtxNJakLtba10DhG03aCnMazZ0MmlazexfU8v7aNsQTmecyUdnuFbkprQ6hVLuOiajQctPXGDnsZQWs9f+m9bWs8PHDZEj+fc0vlZh37/saB6k2v4jogzgb8HWoCvpJQuGfT6FODrwMuALmBVSunB4msXAe8G+oEPpJTWZli6KiivGRrP9S+2ZlaJDXpUXWP9f+9Q6/kPd/54zs0j9NfjPxbq8dx6q3e851ZbbuE7IlqALwJvBLYBt0bEtSmlu8sOezfweErp+Ig4B/gMsCoiTgLOAV4EtAM/iogXppQO/o6Rs3r8wGV9bl4zNJ5b+3+xQfPtvpj179cNemrXeP7fG896/vGcm0for7d/LNTjufVW73jPzUKklPK5cMSrgI+nlFYUn18EkFK6uOyYtcVj1kXEROARYA5wYfmx5ccd6prTp09PL3vZy6rx23mOJ2f/Gl0vOJPUMunAWPT3Mfv+/2Ba1/94btHDHRfQP2Xmc8ZbnnmCBRsuP+Q1Pbf6547nmuP5PJXOf6z9NcQRbUzc181RW382ovPq2d3buwE4qX1GJueN91wd3pOzf43HF76O/skzaBnF57jevlcAPHjqhyDiuS+kxKJbPluVc8dzzXr8M663vwfq8dxK+elPf3pbSmnZUK/l2e1kHvBw2fNtxbEhj0kp7QeeAGaP8FwAIuKCiFgfEev7+voqVPrhPb7wdQeFDoDUMonHF77Oc8v0Tx76L/zhxj0323PHc83xfJ5KwX3CkUcREfRPmUnXC87kydm/dthzpVpR+hz3T5kJo/wcj+f/vaO2/ozoP/jvu+jv46itP6vquS37ukc1Xolzx3PNevt+nNe59VbveM/NQsPfcJlSuhy4HGDZsmXpxhtvzOS6iy/84ZDjA1Nmcrgamunc5ZfcMGTHhXlHHXHYa3pu9c8dzzXH83lafskNpEHXTS2TmPLyt3Hjha8/5LmVkNdyl7Fet9mW59SL8XyOx/P/Hgy15GspKztWjqjusZ47+Ef9ULiJ9+JVp7Hy879flXPHc816+36c17n1Vu94z62UGOonMkUjmvmOiOURcWTx6/Mi4m8j4vnjrKsTWFD2fH5xbMhjistOZlK48XIk5+ZquHZdI2nj1Uznrl6xhNZJLQeNjbTjgudW/9zxXHM8n6dK9KBeddm6A6FUysN4Psfj+X8PCutaf37h63ngkt/i5xe+flTrXMd67sqOeVx89snMa2slgHltrVx89skjOn+s547nmvX2/Tivc+ut3vGem4WRznz/X+CUiDgF+HPgKxS6kJw2jmvfCpwQEYspBOdzgMH/TL0WOB9YB7wVuCGllCLiWuDKiPhbCjdcngD8chy1VNx42ng107mlb5BjuSnPc6t/7niuOZ7P03h7UK/Z0MmGrXvY1z/A8ktuyPQu9zxmoZ3xrk3j+RyP5/+9PK3smDfmGsd67njOg/r5fpzXufVW73jPzcKIbriMiP9OKb00Ij4KdKaUvloaG9fFI34T+DsKrQa/llL6VER8AlifUro2IqYC3wA6gN3AOSml+4vn/iXwh8B+4H+llP79cNdbtmxZWr9+/XhKHpV66jqS97lqTONpazXkj5JHMKM1nnNLxhOg8zq3mWT9Z7xmQ+eYWjJW4rMoqT5FxLA3XI40fP8U+A/gD4DXAY8Bv0opnVzJQqst6/AtaezGGtyHXevX1srPR7BefKxBq8TwXX1Z/hmPN0A7KSE1p0OF75EuO1lFYUnIu1NKj0TEQuDSShUoSYON9UfJ41lnWwpa+/oHgNrrDavsjaePNIxvGYakxjSiGy5TSo+klP42pXRT8fnWlNLXq1uaJI3eeG70PFTQUv0r3QtwywO7WX7JDazZcPj79Ctx868klTtk+I6InojoHuJXT0QcvommJGVsPHe5G7Qa13A/1ThcAB/PP+YkaSiHDN8ppekppRlD/JqeUqqNTuWSVKbUemxyS+Hb22hajxm0GtdYf6pR6y3LJNWfUW2yExHHAFNLz1NKWytekSSN08qOeXz7l4VvT6O5KW887RFV28b6U43SP9rGcxOuJJUbUfiOiLOAz1Hoqf0Y8HzgHuBF1StNksZuLJ0wDFqNa7w9t/0MSKqUEd1wCXwSeCVwb0ppMfAG4BdVq0qScrKyYx4dC9s4dfGsUe8MqNrl8hFJtWKk4bsvpdQFTIiICSmlnwBD9i6UpGY1lm4aGp2x/hmP514ASaqkka753hMR04CfAd+KiMeAp6pXliTVF3uEV994/4zHei+AJFXSSGe+3wz0Av8vhZ0u7wN+p1pFSVK9GW+PcGfND68Sfdivet+rDN6ScjWime+UUvks9xVVqkWSasJYwpk7a1affdglNYIRzXwP2mxnb0T0u8mOJD3LnTWrzz7skhrBSLeXP7DZDtAK/C7wj1WtTJLqiDtrVp8dSyQ1gpGu+T4gFawBVlShHkmqS+6sWX12LJHUCEa6yc7ZZU8nUGgzuLcqFUlSnXJnzeqzY4mkejfSVoPlnU32Aw9S6IAiSRond9aUpOYx0m4nf1DtQiSpmTmjK0nN4ZDhOyL+D5CGez2l9IGKVyRJkiQ1qMPdcLkeuA2YCrwU2Fz8tRSYPNaLRsSsiLg+IjYXH48a5rjzi8dsjojzy8ZvjIhNEXF78dcxY61FkiRJysohZ75TSlcARMT7gdeklPYXn38JuGkc170Q+HFK6ZKIuLD4/CPlB0TELOBjFG7uTMBtEXFtSunx4iHnppTWj6MGSZIkKVMjbTV4FDCj7Pm04thYvZlnd8q8Alg5xDErgOtTSruLgft64MxxXFOSJEnK1Ui7nVwCbIiInwABvA74+Diue2xKaUfx60eAY4c4Zh7wcNnzbcWxkn+KiH7gX4G/TikNuTY9Ii4ALgBYuHDhOEqWJEmSxmek3U7+KSL+HTi1OPSRlNIjhzonIn4EPG+Il/5y0HuniBj2ps5hnJtS6oyI6RTC9zuArw9T++XA5QDLli0b7XUkaVTsVFJ9/hlLqmeH63byayml/4mIlxaHSjPR7RHRnlL67+HOTSmdcYj3fTQi5qaUdkTEXOCxIQ7rBE4vez4fuLH43p3Fx56IuBJ4BcOEb0mSJKlWHG7m+4MUlmx8bojXEvD6MV73WuB8CstZzge+N8Qxa4FPl3VC+Q3gooiYCLSllHZFxCTgt4EfjbEOSZIkKTOH63ZyQfHx1yt83UuA70bEu4GHgLcBRMQy4I9SSu9JKe2OiE8CtxbP+URx7EhgbTF4t1AI3l+ucH2SJElSxY1ozXdE/B7wH8VlHn9Foef3J1NKG8Zy0ZRSF/CGIcbXA+8pe/414GuDjnkKeNlYritJkiTlaaStBv+/YvB+DXAG8FXgS9UrS5IkSWo8Iw3f/cXH3wIuTyn9kHHscClJqn+rLlvHqsvW5V2GJNWVkYbvzoi4DFgFXBcRU0ZxriRJkiRGHqDfRqH7yIqU0h5gFrC6alVJkiRJDWhE4Tul9DSFXtyvKQ7tBzZXqyhJkiSpEY2028nHgGXAEuCfgEnAN4Hl1StNkpqLOzdKUuMb6bKTtwBnAU8BpJS2A9OrVZQkSZLUiEYavvellBKFXS0pbnQjSZIkaRQOG74jIoAfFLudtEXEe3FXSUmSJGnUDrvmO6WUijtcfhDoprDu+6MppeurXZwkqTGV+oO7zl1SsxnRDZfAfwN7Ukq2F5QkSZLGaKTh+1Tg3Ih4iOJNlwAppZdUpSpJkiSpAY00fK+oahWSJElSExhR+E4pPVTtQiRJkqRGN9JWg5IkSZLGyfAtSZIkZcTwLUkatTUbOtmwdQ+3PLCb5ZfcwJoNnXmXJEl1IZfwHRGzIuL6iNhcfDxqmOP+IyL2RMQPBo0vjohbImJLRFwVEZOzqVyStGZDJxdds5F9/QMAdO7p5aJrNhrAJWkE8pr5vhD4cUrpBODHxedDuRR4xxDjnwE+n1I6HngceHdVqpQkPcelazfR29d/0FhvXz+Xrt2UU0WSVD/yCt9vBq4ofn0FsHKog1JKPwZ6yseK292/Hrj6cOdLkipv+57eUY1Lkp6VV/g+NqW0o/j1I8Cxozh3NoXdNvcXn28D5lWyOEnS8NrbWkc1Lkl6VtXCd0T8KCLuHOLXm8uPSyklIFWxjgsiYn1ErN+5c2e1LiNJTWP1iiW0Tmo5aKx1UgurVyzJqSJJqh8j3eFy1FJKZwz3WkQ8GhFzU0o7ImIu8Ngo3roLaIuIicXZ7/nAsHf5pJQuBy4HWLZsWdVCviQ1i5UdhR82fvjqO9jXP8C8tlZWr1hyYFySNLy8lp1cC5xf/Pp84HsjPbE4U/4T4K1jOV+SNH4rO+bRsbCNUxfP4ucXvn5Uwds2hZKaWV7h+xLgjRGxGTij+JyIWBYRXykdFBE3Af8CvCEitkXEiuJLHwE+GBFbKKwB/2qm1UuSxsQ2hZKaXdWWnRxKSqkLeMMQ4+uB95Q9f+0w598PvKJqBUqSquJQbQpdtiKpGbjDpSQpM7YplNTsDN+SpMzYplBSszN8S5IyY5tCSc0ulzXfkqTmZJtCSc3O8C1JytTKjnl8+5dbAbjqfa/KuRpJypbLTiRJkqSMGL4lSZKkjBi+JUmSpIwYviVJkqSMGL4lSZKkjBi+JUmSpIwYviVJkqSM2OdbkjQm9uiWpNFz5luSJEnKiOFbkiRJyojhW5IkScqI4VuSJEnKiDdcSpIy582akppVLjPfETErIq6PiM3Fx6OGOe4/ImJPRPxg0Pg/R8QDEXF78dfSbCqXJEmSxi6vZScXAj9OKZ0A/Lj4fCiXAu8Y5rXVKaWlxV+3V6NISZIkqZLyCt9vBq4ofn0FsHKog1JKPwZ6sipKkiRJqqa8wvexKaUdxa8fAY4dw3t8KiLuiIjPR8SU4Q6KiAsiYn1ErN+5c+eYipUkSZIqoWo3XEbEj4DnDfHSX5Y/SSmliEijfPuLKIT2ycDlwEeATwx1YErp8uIxRMTOiHholNeqZ0cDu/IuQg3Dz5Mqzc+UKsnPkyptPJ+p5w/3QtXCd0rpjOFei4hHI2JuSmlHRMwFHhvle5dmzZ+JiH8CPjTC8+aM5jr1LiLWp5SW5V2HGoOfJ1WanylVkp8nVVq1PlN5LTu5Fji/+PX5wPdGc3IxsBMRQWG9+J0VrU6SJEmqgrzC9yXAGyNiM3BG8TkRsSwivlI6KCJuAv4FeENEbIuIFcWXvhURG4GNFH4k8NeZVi9JkiSNQS6b7KSUuoA3DDG+HnhP2fPXDnP+66tXXUO5PO8C1FD8PKnS/Eypkvw8qdKq8pmKlEZ7r6MkSZKkschr2YkkSZLUdAzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYM35IkSVJGDN+SJElSRgzfkiRJUkYm5l1Alo4++ui0aNGivMuQJElSA7vtttt2pZTmDPVaU4XvRYsWsX79+rzLkCRJUgOLiIeGe81lJ5IkSVJGDN+SJElSRgzfkiRJUkaaas23JEmSRq+vr49t27axd+/evEupKVOnTmX+/PlMmjRpxOcYviVJknRI27ZtY/r06SxatIiIyLucmpBSoquri23btrF48eIRn+eyE0mSJB3S3r17mT17tsG7TEQwe/bsUf80wPAtSZKkwzJ4P9dY/kwM35IkSVJGDN+SJEmqaQ8++CAvfvGLx/0+N954IzfffPOQrz3zzDOcccYZLF26lKuuuor3vOc93H333QB8+tOfHve1S7zhUpIkSU3hxhtvZNq0abz61a9+zmsbNmwA4Pbbbwdg1apVB1779Kc/zV/8xV9UpAZnviVJklTz9u/fz7nnnsuJJ57IW9/6Vp5++mkAbrvtNk477TRe9rKXsWLFCnbs2AHAF77wBU466SRe8pKXcM455/Dggw/ypS99ic9//vMsXbqUm2666cB7P/bYY5x33nnceuutLF26lPvuu4/TTz+d9evXc+GFF9Lb28vSpUs599xzx/37iJTSuN+kXixbtiytX78+7zIkSZLqyj333MOJJ55YePLvF8IjGyt7geedDG+6ZNiXH3zwQRYvXsx//dd/sXz5cv7wD/+Qk046iT/7sz/jtNNO43vf+x5z5szhqquuYu3atXzta1+jvb2dBx54gClTprBnzx7a2tr4+Mc/zrRp0/jQhz70nGvceOONfPazn+UHP/gBAKeffjqf/exnWbZsGdOmTePJJ58csraD/myKIuK2lNKyoY532YkkSZJq3oIFC1i+fDkA5513Hl/4whc488wzufPOO3njG98IQH9/P3PnzgXgJS95Ceeeey4rV65k5cqVudU9mOFbkiRJI3eIGepqGtzWLyJIKfGiF72IdevWPef4H/7wh/zsZz/j+9//Pp/61KfYuLHCs/Vj5JpvSZIk1bytW7ceCNlXXnklr3nNa1iyZAk7d+48MN7X18ddd93FwMAADz/8ML/+67/OZz7zGZ544gmefPJJpk+fTk9Pz6ivPWnSJPr6+iry+zB8S5IkqeYtWbKEL37xi5x44ok8/vjjvP/972fy5MlcffXVfOQjH+GUU05h6dKl3HzzzfT393Peeedx8skn09HRwQc+8AHa2tr4nd/5Hf7t3/7tOTdcHs4FF1xwYBnLeHnDpSRJkg5pqJsKVTDaGy6d+ZYkSZIyYviWJEmSMmL4lqRm1tcL3dvzrkKSmobhW5Ka2A/+7v088XevhCa6/0eS8mT4lqRmlRKv2HszMweegKe78q5GkppCTYbviPhaRDwWEXcO83pExBciYktE3BERL826Rkmqd2tuuo23PPUXLN77TZb/3a2s2dCZd0mSGsiqy9ax6rLnbn7T7GoyfAP/DJx5iNffBJxQ/HUB8H8zqEmSGsaaDZ1ctHYHncwhMYHOngEuumajAVxS3fj4xz/OZz/72WFfX7NmDXfffXeGFY1MTYbvlNLPgN2HOOTNwNdTwS+AtoiYm011klT/Ll27id7+g/8K6O3r59K1m3KqSFIjWbOhkw1b93DLA7tZfskNufzD3vBdWfOAh8uebyuOPUdEXBAR6yNi/c6dOzMpTpJq3fY9vaMal6SRWrOhk4uu2ci+/gEAOvf0Vuwna5/61Kd44QtfyGte8xo2bSpMFnz5y1/m5S9/Oaeccgq/+7u/y9NPP83NN9/Mtddey+rVq1m6dCn33XffkMfloV7D94illC5PKS1LKS2bM2dO3uVIUk1oP3KY8bbWbAuR1HAuXbuJ3r7+g8Yq8ZO12267je985zvcfvvtXHfdddx6660AnH322dx666386le/4sQTT+SrX/0qr371qznrrLO49NJLuf322znuuOOGPC4P9Rq+O4EFZc/nF8ckSSOwet6dtPLMQWOtk1pYvWJJThVJahTV+snaTTfdxFve8haOOOIIZsyYwVlnnQXAnXfeyWtf+1pOPvlkvvWtb3HXXXcNef5Ij6u2eg3f1wLvLHY9eSXwREppR95FSVJdSImVT3yTixfcyuSWCUBiXuzi4re8mJUdQ67gk6QRG+4naNX6ydq73vUu/uEf/oGNGzfysY99jL17947ruGqryfAdEd8G1gFLImJbRLw7Iv4oIv6oeMh1wP3AFuDLwP+TU6mSVH923QtPbGXlK46nY2EbHzvmJn4+5QOsXDI178okNYDVK5bQOqnloLFK/GTtda97HWvWrKG3t5eenh6+//3vA9DT08PcuXPp6+vjW9/61oHjp6bFK9EAACAASURBVE+fTk9Pz4Hnwx2XtYm5XfkQUkpvP8zrCfjjjMqRpMay+frC4wlvhP/uZGfLsYXnex6CI4/Ory5JDaH0E7QPX30H+/oHmNfWyuoVS8b9k7WXvvSlrFq1ilNOOYVjjjmGl7/85QB88pOf5NRTT2XOnDmceuqpBwL3Oeecw3vf+16+8IUvcPXVVw97XNYiNdGWwsuWLUvr16/PuwxJytcVZ8GTj8If31J4/sid8KXl8Hv/DC96S66lSapN99xzDyeeeOKoziltsHPV+15VjZJqxlB/NhFxW0pp2VDH1+TMtySpSp55Eraug1dc8OxYW/H+9T1b86lJUkNq9NA9VjW55luSVCUP/Az698EJv/Hs2NSZMLXN8C1JGTB8S1Iz2XI9TJ4GCwfNSLUtNHxLOqRmWqo8UmP5MzF8S1KzSKlws+Xi02Di5INfM3xLOoSpU6fS1dVlAC+TUqKrq4upU0fXKco135LULHZugicehtd+8LmvtT0f7ruhENAjsq9NUk2bP38+27ZtY+fOnXmXUlOmTp3K/PnzR3WO4VuSmsWWYovB49/43NfaFkLf0/B0l+0GJT3HpEmTWLx4cd5lNASXnUhSs9j8nzDnxGe7m5RrW1h43PNQtjVJUpMxfEtSM3imBx5aByecMfTrB8K3674lqZoM35LUDB74GQz0HdxisJy9viUpE4ZvSWoGm/+z0GJwwSuHft1e35KUCcO3JDW6lGDzj+AFpz+3xWA52w1KUtUZviWp0e38H+jeBicM0eWknOFbkqrO8C1JjW7zfxYeh2oxWK7t+YXw7SYaklQ1hm9JanSbr4djToKZ8w59XHmvb0lSVRi+JamR7e2Grb84/JITsNe3JGXA8C1JjeyBnxZaDB5uyQnY61uSMmD4lqRGtvl6mDwdFg7TYrCcvb4lqeoM35LUqFKCLT+C406HlkmHP95e35JUdYZvSWpUj90N3Z0jW3JSYrtBSaoqw7ckNarN1xcejz9j5OcYviWpqgzfktSotvwIjn3x4VsMlrPXtyRVleFbkhrR3m7Yum50s95gr29JqjLDtyQ1ovtvhIH9I+vvXc5e35JUVYZvSWpEW66HKTNgwamjO89e35JUVYZvSWo0KcHmH8ELTh9Zi8Fy9vqWpKoyfEtSo3n0LujZPvolJ2Cvb0mqMsO3JDWaLaUWg2MI32C7QUmqIsO3JDWazT+CY0+GGXPHdr7hW5KqxvAtSY1k7xOFFoMnjLLFYDl7fUtS1Ri+JamR3H8jpP6xLzkBe31LUhUZviWpkWy+HqbMhAWvGPt72OtbkqrG8C1JjSKlwpbyx50++haD5ez1LUlVY/iWpEbx6J3Qs2N8S07AXt+SVEWGb0lqFJtLLQbHcbMl2OtbkqrI8C1JjWLz9fC8cbQYLGe7QUmqCsO3JDWC3j3w8C3jX3JSYviWpKowfEtSIyi1GDzhNyrzfvb6lqSqMHxLUiPYfH1hrfb8l1fm/ez1LUlVYfiWpDq36ks3s/uO6+AFvw4tEyvzpvb6lqSqMHxLUp17/v77mTWwu3JLTsBe35JUJYZvSapzS/bdXfjiBadV7k3t9S1JVWH4lqQ6N7t/J/uZCNPbK/em9vqWpKowfEtSnZs10MXjLbNgQoW/pdtuUJIqzvAtSXXuqP4udk84uvJvbPiWpIozfEtSnZvdv5PdLbMr/8b2+pakijN8S1I9S4lZA13sbqnSzLe9viWpogzfklTP9j7B1LS3estOwF7fklRBhm9Jqmc9OwCqtOzEXt+SVGmGb0mqZ93bAaq07MRe35JUaYZvSapjazbuYvnev+eKzrksv+QG1mzorNyb2+tbkipuYt4FSJLGZs2GTi669Qh6mQZA555eLrpmIwArO+ZV5iKHaDe46rJ1AFz1vldV5lqS1ASc+ZakOnXp2k30Dhz8bby3r59L126q3EXs9S1JFWX4lqQ6tX1P76jGx8Re35JUUYZvSapT7W2toxofE3t9S1JFGb4lqU6tXrGEVp45aKx1UgurVyyp3EXs9S1JFWX4lqQ6tfLko7l44pdpm1BYZjKvrZWLzz65cjdbgr2+JanC7HYiSfWqZwcrJ97MjpmncuMRK6rTdcRe35JUUc58S1K9quYGOyX2+pakijJ8S1K9KoXvCVXYWr6c7QYlqWIM35JUr3p2ANDVMqe61zF8S1LFGL4lqV51b4dJR9IbR1T3Ovb6lqSKMXxLUr3q3g4z5kJEda9jr29JqhjDtyTVq+7tMKO9+tcZotf3mg2dbNi6h1se2M3yS25gzYbO6tchSQ3A8C1J9apnB0zPMnwX1n2v2dDJRddsZF//AACde3q56JqNBnBJGgHDtyTVo4GBQvieMbf61xrU6/vStZvo7es/6JDevn4uXbup+rVIUp0zfEtSPXpqJwzshxkV3M1yOIN6fW/f0zvkYcONS5Ke5Q6XklSPego9vpk+tzo7Ww5W1m6wva2VziGCdntba/XrkKQ658y3JNWj7kKP70xuuISDwvfqFUtondRy0Mutk1pYvWJJNrVIUh0zfEtSPeou3tyYWfh+ttf3yo55XHz2yUxuKfwVMq+tlYvPPpmVHRksgZGkOueyE0mqRz07IFrgyCrvbllS3uv7yKNZ2TGPb/+yMBOeybIXSWoQznxLUj3q3g7T58KElsMfWwlD9PqWJI2e4VuS6lFpd8usDOr1LUkaG8O3JNWjnh2Fme+sDOr1LUkaG8O3JNWj7u3Z9PguGdTrW5I0NoZvSao3e7th35PZLjuBg9oNSpLGxvAtSfWmp9jje3pGbQZLDN+SNG41G74j4syI2BQRWyLiwiFeXxgRP4mIDRFxR0T8Zh51SlLmsu7xXVLW61uSNDY1Gb4jogX4IvAm4CTg7RFx0qDD/gr4bkqpAzgH+Mdsq5SknBzY3TKHZSelXt+SpDGpyfANvALYklK6P6W0D/gO8OZBxyRgRvHrmcD2DOuTpPz0FL/d5bHsBOz1LUnjUKvhex7wcNnzbcWxch8HzouIbcB1wJ8O9UYRcUFErI+I9Tt37qxGrZKUre7t0DoLJk3N9rr2+pakcavn7eXfDvxzSulzEfEq4BsR8eKU0kD5QSmly4HLAZYtW+ZCRUn1r3tH9uu94Tm9vt1WXpJGr1ZnvjuBBWXP5xfHyr0b+C5ASmkdMBU4OpPqJClP3Z35hG97fUvSuNVq+L4VOCEiFkfEZAo3VF476JitwBsAIuJECuHbdSWSGl/Wu1uWs92gJI1LTYbvlNJ+4E+AtcA9FLqa3BURn4iIs4qH/Tnw3oj4FfBt4F0p2f9KUoPbvw+e2pnPzDcYviVpnGp2zXdK6ToKN1KWj3207Ou7geVZ1yVJuSptsJNb+H4+3HdDodd3RD41SFIdq8mZb0nSMPLa3bLEXt+SNC6Gb0mqJ93FHt95LjsBe31L0hgZviWpnhwI3znecAmu+5akMTJ8S1I96dkBE1sLLf/yMKjXtyRpdAzfklRPSj2+87rZ0V7fkjQuhm9Jqid57W5ZznaDkjRmhm9Jqic92/PbYKfE8C1JY2b4lqR6MTBQIzPfzy+Eb/c1k6RRM3xLUr14ugsG+mogfNvrW5LGyvAtSfWiJ+ce3yX2+pakMTN8S1K9KPX4zmt3yxJ7fUvSmBm+Jale5L3BTom9viVpzAzfklQvurdDtMC0Y/Otw17fkjRmhm9Jqhc9OwrBe0JL3pXYblCSxsjwLUn1ont7/ktOSgzfkjQmhm9Jqhfd2/PvdFJir29JGhPDtyTVi54d+Xc6KbHXtySNieFbkurBMz3wTHcNzXzb61uSxsLwLUn1oHtH4bHmwrfrviVpNAzfklQPSrtbTq+VGy7t9S1JY2H4lqR60F0jW8uX2OtbksbE8C1J9aDWwjfYblCSxsDwLUn1oGdHYaZ5UmvelTzL8C1Jo2b4lqR60L0dZszLu4qDzVwAex6217ckjYLhW5LqQS3tblnStgD6noLex/OuRJLqhuFbkupBz47a6XRSMrPY8eSJh/OtQ5LqiOFbkmpdfx88+VjtLTs50G7Q8C1JI2X4lqRa1/MIkGpv2cmBme9t+dYhSXXE8C1Jte5Am8Eam/k+YjZMbHXZiSSNguFbkmpdre1uWRJRWHpiu0FJGjHDtyTVuu4dhcda2mCnZOZ8Z74laRQM35JU67o7YeJUaD0q70qea+YC13xL0igYviWp1pXaDEbkXclztS2Ap3ZCX2/elUhSXTB8S1Kt695Rm0tOAGYuLDw6+y1JI2L4lqRa191Zw+F7fuHRdd+SNCKGb0mqZSkV+nzXWqeTEjfakaRRMXxLUi17ejf0P1N7Pb5LprdDTHDmW5JGyPAtSbWsu7PwWGu7W5a0TCwEcGe+JWlEDN+SVMt6ij2+p9fomm8oLD3xhktJGhHDtyTVsgMz3zUcvmcugCfc5VKSRsLwLUm1rHtHYU31tGPzrmR4bQugezsM9OddiSTVPMO3JNWynu1w5DGFtdW1auZ8GNhf6MoiSTokw7ck1bLu7bW95ATKNtrxpktJOhzDtyTVslre3bLEXt+SNGKGb0mqZT31MPNd2uXSmy4l6XAM35JUq/Y9BXufqN3dLUsmHwmts2w3KEkjYPiWpFrVXezxXesz31BYeuKyE0k6LMO3JNWqeujxXTJzgTdcStIIGL4lqVbVw+6WJTOLu1ymlHclklTTDN+SVKu6txceZ9T4mm8oLDvZ9yT0Pp53JZJU0wzfklSrurfD1JmFGxpr3cxiu0GXnkjSIRm+JalW9eyojyUn8Gy7QW+6lKRDMnxLUq2qh90tS9pKu1zablCSDsXwLUm1qnt7faz3BjhiNkxsddmJJB2G4VuSalF/Hzz5aP0sO4ko9vp2l0tJOhTDtyTVoicfBVL9LDuBwrpvl51I0iEZviWpFtXT7pYlbrQjSYdl+JakWtRT7PE9vU7WfENh2clTO6GvN+9KJKlmGb4lqRYd2GBnXr51jMaBXt8uPZGk4Ri+JakWdW+HlilwxKy8Kxk5N9qRpMMyfEtSLerZAdOfV+giUi/aiuHbjXYkaViGb0mqRd3b62vJCRTaIsYEZ74l6RAM35JUi+ppg52SlomFAO6ab0kaluFbkmpNSvW1tXy5tgUuO5GkQ5iYdwGSpIJVl60D4Kp3LIH+Z+pnd8tyMxfAw7/IuwpJqlnOfEtSrTnQZrDOlp1AYZfL7u0w0J93JZJUkwzfklRr6rHHd0nbAhjYDz2P5F2JJNUkw7ck1Zp63N2yZObCwqMdTyRpSIZvSao13TuAKPT5rjf2+pakQzJ8S1Kt6e6EacdAy6S8Kxm9mfMLj858S9KQDN+SVGt6dtTnkhOAyUdC6yzDtyQNw/AtSbWmHne3LGevb0kaluFbkmpNPe5uWW7mAme+JWkYhm9JqiGT017Yu6d+l51AMXxvK+zUKUk6iOFbkmrIrP6uwhf1vuxk35PQ+3jelUhSzTF8S1INmdW/q/BFvS87AZeeSNIQDN+SVENmDRRnvqe351vIeBxoN7gt3zokqQbVbPiOiDMjYlNEbImIC4c55m0RcXdE3BURV2ZdoyRVypoNnWzYuofPPXIKy/f+PWseqNlvz4fXVtzl0o4nkvQcE/MuYCgR0QJ8EXgjsA24NSKuTSndXXbMCcBFwPKU0uMRcUw+1UrS+KzZ0MlF12xkX/8AEHQyh4u+vwUmtbKyow7Xfh8xGya2uuxEkoZQq1MrrwC2pJTuTyntA74DvHnQMe8FvphSehwgpfRYxjVKUkVcunYTvX39B4319vVz6dpNOVU0ThGFpSd7tuZdiSTVnFoN3/OA8imTbcWxci8EXhgRP4+IX0TEmUO9UURcEBHrI2L9zp07q1SuJI3d9j29oxqvC20LXPMtSUOo1fA9EhOBE4DTgbcDX46ItsEHpZQuTyktSyktmzNnTsYlStLhtbe1jmq8LrjRjiQNqVbDdyewoOz5/OJYuW3AtSmlvpTSA8C9FMK4JNWV1SuW0Dqp5aCx1kktrF6xJKeKKqBtATy1E/rqePZekqqgVsP3rcAJEbE4IiYD5wDXDjpmDYVZbyLiaArLUO7PskhJqoSVHfO4+OyTmdqSCAaY17qfi88+uT5vtiw50Ot78LyJJDW3mgzfKaX9wJ8Aa4F7gO+mlO6KiE9ExFnFw9YCXRFxN/ATYHVKqSufiiVpfFZ2zOPs5+3igann8fNzJtd38Iay8O1Nl5JUriZbDQKklK4Drhs09tGyrxPwweIvSap7s/uLN4XPqOMNdkraiuF7mF7fqy5bB8BV73tVVhVJUk2oyZlvSWpGDbG7Zcn0uRAT7HgiSYMYviWpRhzV30Ufkwqb1NS7lkmFf0TY8USSDmL4lqQaMbt/J4+3zIIJDfKtuW2BW8xL0iAN8h1ekurfrIEudk9ogFnvkpkLvOFSkgYxfEtSjZjVv4vdLQ20GdjM+dC9HQb6865EkmqG4VuSakFKzOrvYndLA818ty2Agf3Q80jelUhSzTB8S1It2LuHKTzTYMtOFhYevelSkg4wfEtSLejeDkBXoy07AdsNSlIZw7ck1YI9hRsTu1qOzrmQCjqw0Y43XUpSieFbkmrBrs0AbJ+4IOdCKmjykdA6y2UnklTG8C1JtWDXvXDkHL72/t/Iu5LKGqLX95oNnWzYuodbHtjN8ktuYM2GzpyKk6TsGb4lqRZ0bYHZJ+RdReXNXHDQmu81Gzq56JqN7OsfAKBzTy8XXbPRAC6paRi+JakW7LoXjm7U8P0wpATApWs30dt3cN/v3r5+Ll27KY/qJClzhm9JytvTu+HprsYM320LYN+T0Ps4ANv39A552HDjktRoDN+SlLfizZYc/cJ866iGQe0G29tahzxsuHFJajSGb0nK2657C4+NOPM9s9i9pdjxZPWKJbROajnokNZJLaxesSTryiQpFxPzLkCSml7XZmiZDG3Pz7uSymsr7nJZ7HiysmMeAB+++g729Q8wr62V1SuWHBiXpEZn+JakvO3aDLOOgwkthz+23hwxGya2HtTre2XHPL79y8LGO1e971V5VSZJuXDZiSTlbde9cPTxeVdRHRGFdd9utCNJgOFbkvLV3wePP9iYN1uWDLHRjiQ1K8O3JOVp9wMwsL+xw3ep17ckyfAtSbnqKrYZbMTdLUtmLoCndkKfvbwlyfAtSXk60GawQdd8Q2HZCcATbiEvSYZvScrTrs0w7ViYOjPvSqrnQK/vrfnWIUk1wPAtSXnatbmx13vDs7tcetOlJBm+JSk3KRXbDDbwem+AGe0QEw5sMS9JzczwLUl5eboL9u5p7JstAVomwfR2O55IEu5wKUn5OXCzZYMvO4Hn9Pp2Z0tJzcqZb0nKSzN0Oilxl0tJAgzfkpSfXZth4tRnu4E0spkLoLsTBvrzrkSScmX4lqS87NoMs4+HCS15V1J9bQsKO3n2PJJ3JZKUq6qG74j404g4qprXkKS61VUM383gQK9vl55Iam7Vnvk+Frg1Ir4bEWdGRFT5epJUH/Y/A48/2Bw3W0JZ+LbdoKTmVtXwnVL6K+AE4KvAu4DNEfHpiDiumteVpJq3+35IA43f47uktMX8Hne5lNTcqr7mO6WUgEeKv/YDRwFXR8TfVPvaklSzdm0uPDZL+J58JLTOctmJpKZX1T7fEfFnwDuBXcBXgNUppb6ImABsBj5czetLUs0qtRls9A12ys2c77ITSU2v2pvszALOTik9VD6YUhqIiN+u8rUlqXZ1bSns+jhlWt6VZKdtIXTdl3cVkpSraq/5/tjg4F322j3VvLYk1bRd9zbPkpOSmQsKy05SyrsSScqNfb4lKWspFdZ8N134ng/7noTex/OuRJJyY/iWpKw9+Rg80908bQZL2mw3KEmGb0nKWulmy6ab+XajHUkyfEtS1rqKbQabqdMJFG64BNhj+JbUvAzfkpS1XZth0hEwY17elWTriNkwsdWZb0lNzfAtSVnbdS/MPg4mNNm34Ihir2/Dt6Tm1WTf+SWpBuza3Hw3W5a0LXDZiaSmZviWpCz19cKerc0bvp35ltTkDN+SlKXd9wMJZh+fdyX5mLkQntpZ+EeIJDUhw7ckZelAm8Emnfk+0Ou7M986JCknhm9JytKuUpvB4/KtIy8Hen1vzbcOScqJ4VuSsrRrcyGATj4y70ryMXN+4dFdLiU1KcO3JGVp173Nt7NluRntEBPseCKpaRm+JSkrKUHXlubb2bJcyySY3m7HE0lNy/AtSVnp2QH7nmzumW8oLD1x5ltSkzJ8S1JWDnQ6afLw3bbAmW9JTcvwLUlZKXU6adY2gyUzF0B3Jwz0512JJGXO8C1JWdm1GSZPg+lz864kX20LYGA/9DySdyWSlDnDtyRlpWtzYWfLiLwrydeBXt+2G5TUfAzfkpSVXZtdcgJl4dt135Kaj+FbkrKw76lC2Gz2my3h2Y129rjLpaTmY/iWpCx03Vd4NHzDlGnQepQz35KakuFbkrJwoM2gy06AwkZDj/1P3lVIUuYM35KUha4tQMCsF+RdSW1o74BH7rDdoKSmY/iWpCzsuhfaFsKk1rwrqQ3tHYXdPru25F2JJGXK8C1JWdh1r+u9y7V3FB63b8i3DknKmOFbkqptYKBww6XrvZ919Akw6UjDt6SmY/iWpGrr7oS+p535LjehBea+xPAtqekYviWp2ro2Fx5nG74P0t4BO+6A/v15VyJJmTF8S1K17SqGb5edHKy9A/b3wq5NeVciSZkxfEtSte26F6bMgGnH5F1JbfGmS0lNyPAtSdW2a3NhvXdE3pXUllnHweTpsP32vCuRpMwYviWp2nZtdsnJUCZMgPalznxLaiqGb0mqpmd6oGc7zD4+70pqU/tSeGQj9PflXYkkZcLwLUnVVNrB0ZnvobV3QP8z8Ng9eVciSZkwfEtSNR3odGKbwSHNXVp4dOmJpCZh+Jakatq1GWICzHpB3pXUplkvgCkzDd+SmobhW5Kqade9cNQimDgl70pqU4Q3XUpqKoZvSaqmri3ubHk47R3w6F2w/5m8K5GkqqvZ8B0RZ0bEpojYEhEXHuK4342IFBHLsqxPkg5roL8Qvl3vfWjtHTDQB4/dnXclklR1NRm+I6IF+CLwJuAk4O0RcdIQx00H/gy4JdsKJWkEnngY9u81fB+OO11KaiI1Gb6BVwBbUkr3p5T2Ad8B3jzEcZ8EPgPszbI4SRqRXbYZHJG2hdA6y/AtqSnUavieBzxc9nxbceyAiHgpsCCl9MNDvVFEXBAR6yNi/c6dOytfqSQNZ9e9hUfD96FFFGa/Dd+SmkCthu9DiogJwN8Cf364Y1NKl6eUlqWUls2ZM6f6xUlqaqsuW8eqy9YVnnRthqltcMTsfIuqB+1LCxvt9PXmXYkkVVWthu9OYEHZ8/nFsZLpwIuBGyPiQeCVwLXedCmppuzaXJj1jsi7ktrX3gED+wtdTySpgdVq+L4VOCEiFkfEZOAc4NrSiymlJ1JKR6eUFqWUFgG/AM5KKa3Pp1xJGsKue73ZcqS86VJSk6jJ8J1S2g/8CbAWuAf4bkrproj4RESclW91kjQCe5+AJx81fI/UjHlw5BzDt6SGNzHvAoaTUroOuG7Q2EeHOfb0LGqSpBGz08noHLjp8va8K5GkqqrJmW9JqntdmwuP7m45cu0dsPMe2Pd03pVIUtUYviWpGnbdCxMmwqzFeVdSP9o7IA3AIxvzrkSSqsbwLUnVsOteOGoRtEzKu5L6MXdp4dF135IaWM2u+ZakerNmQycbtu5hX/8Ay1vexOr2BazMu6h6MmMuTHue4VtSQ3PmW5IqYM2GTi66ZiP7+gcA6Oxv46LOV7FmQ+dhztRB3OlSUoMzfEtSBVy6dhO9ff0HjfUOtHDp2k05VVSn2jsKS3ae6cm7EkmqCsO3JFXA9j1Db4s+3LiG0d4BJNhxR96VSFJVGL4lqQLa21pHNa5htBdvutxhv29JjcnwLUkVsHrFElontRw01jqphdUrluRUUZ2adgzMmO+6b0kNy/AtSRWwsmMeF599MpNbJhAk5k3s5uKzT2Zlx7y8S6s/7UsN35IaluFbkipkZcc8Oha2seGIP+bnL/+5wXus2jugawvsfSLvSiSp4gzfklRBRw708P+3d/dxVpd1/sdf19wy3DkgoMyICd6QqCBJZuJ9N2hZTlaL1na7m7VlW2071dRutbsVbmPbzW71s6222twsi51cs6jQNMk7DAE1MRAVZsQbcBDlwNxdvz++BxiYMwMMM+d7zpnX8/E4j3PO93u+zMe65pz3XOe6qe1ph0knpF1K8do97ntlunVI0jAwfEvSEKrr2pA8mHR8uoUUs6lzk3uHnkgqQYZvSRpCdV0bkwf2fA/emMOh9mjDt6SSZPiWpCFU17WRLiqg9kVpl1Lc3OlSUokyfEvSEKrr2simijoor0i7lOJWNxeefRQyz6ZdiSQNKcO3JA2h08c9w1HHzU67jOJXt2vct5vtSCothm9JGirdnbDlESdbDoWpc5J7h55IKjGGb0kaKptWQU8XHHFy2pUUv5oJMHGG4VtSyTF8S9JQWXszEGDGeSkXUiKmnuqwE0klx/AtSUNl3dJkuMSYSWlXUhrq5sLWx+GFZ9KuRJKGjOFbkobCjq2w4W447hVpV1I6nHQpqQQZviVpKDxyK8RuONbwPWScdCmpBBm+JWkorFsKVeNg2ulpV1I6Ro2Hw4+HJ+z5llQ6DN+SdKhiTCZbTj8HyivTrqa0uNOlpBJj+JakQ7V5bTIx0PHeQ69uLjzXCtue3H1o4TV3sPCaO1IsSpIGz/AtSYdq7dLk3vA99HZNunToiaQSYfiWpEO1bilMPBYmHJN2JaXnyFOA4NATSSXD8C1Jh6JrJzx6u73ew6V6LEyeafiWVDIM35J0KB6/Azq3u8TgcNo16TLGtCuRpENm+JakQ7F2KZRVwjFnpV1J6aqbC88/CdueSLsSSTpkhm9JXu2negAAIABJREFUOhRrl8LRZyTDIzQ8eu102bKilRWPt3PX+i3Mv+pmWla0plubJB0kw7ckDdZzT8BTDzjee7gdcTKEclqWr6dp8Wo6unsAaG3P0LR4tQFcUlExfEvSYK27Obl3vPfwqhoNU06kec0kMp3de53KdHbTvGRNSoVJ0sEzfEvSYK1bCmOmJD2zGl51p9LWOSbnqbb2TJ6LkaTBM3xL0mD0dMO6W5IhJ2W+lQ67qadSxzM5T9XV1uS5GEkaPD8xJGkwnrgPMlsccpIvdS+hseLH1JTvvdxgTWU5jQtmplSUJB08w7ckDcbam4EAx56fdiUjwxEn0VB1N4tevJ6q8uSjq762hkWXnkLD3PqUi5OkA1eRdgGSVJTWLYWpc2DMpLQrGRkqR8GUWTTEpfzo6HkA/Pi9L0+5KEk6ePZ8S9LB2rEVNtztEoP5Vjc3Ge7jTpeSipjhW5IO1vrbIHY73jvf6uZC5lkmdz+ZdiWSNGiGb0k6WGuXQtU4mHZ62pWMLNmdLo/tfDjlQiRp8AzfknQwYkzC9/RzoLwy7WpGlimzoLyKGZ1/TrsSSRo0w7ckHYzNa2Hr43DcBWlXMvJUVMERJxm+JRU1w7ckHYy1S5N7x3uno24up4T1/Pg9L0u7EkkaFMO3JB2MdUth4gyYOD3tSkamurmwcys8uz7tSiRpUAzfknSgunbCo7fDca9Mu5KRKzvpkrYV6dYhSYNk+JakA/X4HdC53SEnaZr8YqgeD+tuTrsSSRoUw7ckHai1S6GsEo45K+1KRq7ySpj1enjw59CxPe1qJOmgGb4l6UCtuxmOPgOqx6Zdycg253LoeB7W3JR2JZJ00AzfknQgtm2CJ+93S/lCcPSZcNg0WPmjtCuRpINm+JakA7FrjLHjvdNXVgazFyb/n2xzq3lJxcXwLUkHYu1SGDMFjjg57UoEMOcyiD2w+vq0K5Gkg2L4lqT96elOelmPvSDpdVX6Jh0P9afByuvSrkSSDoqfIpK0P0/cB5ktjvcuNHMuhydXw6b7065Ekg6Y4VuS9mdtdrz3jPPTrUN7O+lSKKuAVfZ+Syoehm9J2p91S2HqqTB2ctqVqLcxh8PxC2DV9cnQIEkqAoZvSRrIjq2w4W6HnBSqOQvh+U3wyO/SrkSSDojhW5IGsv42iN0uMVioTrgQRh0Gq36cdiWSdEAM35K0j4XX3MHCa+5InqxdClXjYNrp6Ral3Cqqk7Hff/o/2Lkt7Wokab8M35LUnxiT8d7Tz4HyyrSrUX/mXA6d25MALkkFzvAtSf3ZvA7aH4fjLki7Eg1k2ukwYbprfksqCoZvSeqlZUUrKx5v5671W5j/9Qdo6TrT8d6FLoRkx8v1t8HW1rSrkaQBGb4lKatlRStNi1fT0d0DQGumgqbuK2h5rCrlyrRfsxcCEVb/JO1KJGlAhm9JympesoZM597rRWdiFc1L1qRUkQ7YxOkw7Yxk6EmMaVcjSf0yfEtSVlt75qCOq8DMuQyefgieWJl2JZLUL8O3JGXV1dYc1HEVmJMaoLzKiZeSCprhW5KyGhfMpKayfK9jNZXlNC6YmVJFOig1E2DmRXD/T6G7M+1qJCknw7ckZTXMrWfRpacwqjwS6KG+ppNFl55Cw9z6tEvTgZp9GbzwNKy7Oe1KJCknw7ck9dIwt54PTvoj60f9JcuumGHwLjbHvRJGHw4rf5R2JZKUk+FbkvYxZ+e9tJdNgCNOTrsUHayKKjj5TfDQTZBpT7saSerD8C1JvXVs5yzuo/bkBVDmW2RRmrMQunfCgz9PuxJJ6sNPFknq7d7/gh3tMO9daVeiwap7CUw6AVb9OO1KJKkPw7ck7dKZgWVfhWPOhqPPSLsaDVYIyY6Xjy2DZx9LuxpJ2ovhW5J2+eMP4Pkn4dyPp12JDtXshcn9Krebl1RYDN+SBNC1E27/Chx9JhxzVtrV6FDVTku+wVj5I7ebl1RQDN+SBLDih7CtDc79WDJsQcVvzmWwZR1sXL770MJr7mDhNXekWJSkka5gw3cI4cIQwpoQwtoQwidynP+7EMKDIYRVIYSlIYQXpVGnpBLQ1QG3fxmOeinMOC/tajRUTnw9VNTAqmS7+ZYVrax4vJ271m9h/lU307KiNeUCJY1EBRm+QwjlwNeBi4BZwOUhhFn7vGwFMC/GOBv4KfDF/FYpqWSsug62bkjGetvrXTpGjYcXvxbu/xktyx+jafFqOrp7AGhtz9C0eLUBXFLeFWT4Bk4H1sYYH4kxdgDXAZf0fkGM8ZYY4/bs0zuBo/Jco6RS0N0Jt10NdXOT3RFVWuZcDplnaf7lajKd3XudynR207xkTUqFSRqpCjV81wMbej3fmD3Wn78CfpnrRAjhihDC8hDC8qeffnoIS5RUElZfD+2P2etdqmacB2Om0PZC7tNt7Zl8ViNJBRu+D1gI4S+BeUBzrvMxxm/FGOfFGOdNnjw5v8VJKmw93Umv95GnwAkXpl2NhkN5Bcz+C+rC5pyn62pr8lyQpJGuUMN3KzCt1/Ojssf2EkJ4JfAp4PUxxp15qk1Sqbh/cbIaxjmucFLSZi+ksfw6asp79jpcU1lO44KZKRUlaaQq1PB9D3B8CGF6CKEKuAy4ofcLQghzgWtIgvdTKdQoqZj1dMNtzTBlFrz44rSr0XA68hQa6ray6PBfUVWefOzV19aw6NJTaJg70IhGSRp6FWkXkEuMsSuEcCWwBCgHvhtjfCCE8M/A8hjjDSTDTMYC14ekx+rxGOPrUytaUnF58OfwzBp403ehrFD7ITQkQoA5C2n4zae5pe58NlXU8+P3vjztqiSNUCGOoJ2/5s2bF5cvX77/F0oqbT098P/mQ08XvP9OKCtPuyINt+eegC/PgrP/Hi74VNrVSCpxIYR7Y4zzcp2zu0fSyLPmF/DUg3BOo8F7pBg/NVn5ZNV1yR9fkpQSw7ekkSVGuPVfYeKxcNKlaVejfJp9GbQ/Do/elnYlkkYww7ekkeXhX8Gm1XD2R5Nl6DRynPg6OGwa/PIT0NWRdjWSRijDt6SRI0a49YtQ+yKY/RdpV6N8qxoNr7kanv4T/OFraVcjaYQyfEsaOdYuhbY/Znu9K9OuRmmYeSHMuiT5I2zzurSrkTQCGb4ljQwxwq1XJcMO5lyedjVK04X/ChXVcONHknYhSXlk+JY0MjzyO9h4D5z1YaioSrsapWn8VHjlZ2D9rbDyurSrkTTCGL4ljQy3NcO4Opj7trQrUSE47d1w1Omw5JPwwua0q5E0gjjVX1LJaVnRysd+uoqO7h7qa2toPA0aHlsGF30xGW4glZXB674K15wNv/4HeMM3065I0ghhz7ekktKyopWmxavp6E42Umltz9B083O0VFwEL3l7ytWpoBwxC878W1j5P/DIrWlXI2mEMHxLKinNS9aQ6eze61gmVtLcfRlU1qRUlQrWuR+DCdOTyZedO9KuRtIIYPiWVFLa2jO5j2dcWlA5VNbA674CW9bB769OuxpJI4DhW1JJqavN3bvd33GJGeclW8/f/hV46qG0q5FU4gzfkkpK44KZ1FSW73WsprKMxgUzU6pIRWHB56F6HPzfh6CnJ+1qJJUww7ekktIwt55Fl55C/bhyAj3Uj+pg0aWzaZhbn3ZpKmRjJsGrPwcb7oQ/fj/taiSVMJcalFRyGuZMpeH+K6H1Xvjwahh1WNolqRic+hZY+SP4zWdg5kUw7si0K5JUguz5llR6fvcFWHczXPCPBm8duBDg4q9A1w741SfSrkZSiTJ8Syot9y9OdrOc+zZ46V+nXY2KzaTj4Jy/hwf+Fx7+NZCsHX/Cp37JMZ/4BfOvupmWFa0pFympmBm+JZWOJ1ZCy/th2hnw2i8lPZnSwZr/YZj8YvjFR2m555G+mzYtXm0AlzRohm9JpeH5p+BHb4HRh8PC/3YbeQ1eRVUy/GTr4zTfuKLvpk2d3TQvWZNScZKKneFbUvHr6oAfvw22b4bLroWxU9KuSMXuRS+H095J287cf8T1t5mTJO2P4VtScYsRbvposkRcw9eh7tS0K1KpeOVnqSt7NucpN22SNFiGb0nF7e7/hD/+AM7+ezj5jWlXo1JSM4HGl4+lhp17H64sd9MmSYNm+JZUvB75XbIk3MzXwPmfSrsalaCGixtYdNQd1IdnCEB9bQ2LLj3FTZskDZqb7EgqTlsegevfCZNOgDdcA2X2JWgYhEDDW6+k4RtnwPRz4fIfuYqOpEPip5Wk4rNzW7KyCcDl/wOjxqdbj0rbhBfB+Z+Eh38Jt3whmWcgSYNkz7ek4tLTA4uvgGcehrcthokz0q5II8EZH0ja3G1fhJ5OeMVn7AGXNCiGb0nF5ZbPw5qb4KIvwozz0q5GI0VZGVz8VSirgNu/DN2d8OrPGcAlHTTDt6Ticf/P4PdXw0veDqdfkXY1GmnKyuC1/5YE8Dv+A3q64cJFBnBJB8XwLak4tN0HLR9Ito5/jVvHKyUhJN+6lFXCnV9PhqBc1OyEX0kHzPAtqSC1rGjlYz9dRUd3D/WHVdHY8z0adm8dX5V2eRrJQoAFn4fyClj21WQIysVfMYBLOiCGb0kFp2VFK02LV9PR3QNA69YOmngjnPkBGtw6XoUgBHjlPyU94L+/OhmC8vqvQVl52pVJKnD+mS6p4DQvWUOms3uvYxmqab5rR0oVSTmEABf8A5z7Cbjvh9Dy/iSES9IA7PmWVHDa2jMHdVxKTQhwflMyCfOWzyVjwN/wrWRIiiTl4LuDpIJTd1g1rVt39j1eW5NCNdIBOLcxCdy//Sz0dMEbvwPllWlXJakAGb4lFZanHqIxXEsTryND9e7DNZXlNC6YmWJh0n6c9ZFkDPivP5UMP3nTf9Gy+mmal6yhrT1DXW0NjQtm0jC3Pu1KJaXI8C2pcNz/M/j5B2moGgOveBvN92JoUXE588pkCMqvPk7LNZ+madO5ZDqzE4fbMzQtXg1gW5ZGMMO3pPR1d8JvPg13fiNZx/vN36Nh/FQaXpV2YdIgnPE+KK+g+WcVZOjZ61Sms5vmJWsM39II5monktK1bRN8/3VJ8H7Z38A7b4TxU9OuSjo0L/1r2piU85QTh6WRzZ5vSel5dBlc/07oeCGZoHbKm9KuSBoydbWjac0RtJ04LI1s9nxLyr8Y4Q//nvR4jzoM3rPU4K2S07hgJjWVe2+6U1PeQ+Orjk+pIkmFwJ5vSfm1cxv8/APw4M/hxNfBJd+AUePTrkoacrvGde9e7aRiG418n4aVwIu+BpMM4dJIZPiWNGxaVrTuvczameNoWPV+2LwOXvUvcOYHk01KpBLVMLd+z+TKGOG+bljySfjmmXDux+DMD0FFVbpFSsorh51IGhYtK1ppWrya1vYMkewyazdtoOXZGfCOG2D+3xq8NbKEAHPfClfeAy++GG7+HHzrPNh4b9qVScojw7ekYdG8ZA2Zzu69jmWoprn8r+CYs1KqSioAY6fAm/8LLr8OdrTDt18Bv2qCnc+nXZmkPDB8SxoW/S2n1vZcZ54rkQrUzIvg/XfCS/86WWrzGy+HP/827aokDTPDt6ShFSM8dBN15e05T7vMmtTLqPHw2qvh3UugsgaufSMsvgJe2Jx2ZZKGiRMuJQ2dR34HS/8ZWu+lcdwbaNr2JjLde8Z111SW07hgZnr1SYXq6DPgfb+H338Jfv9vsPa3tMxspvlPtbS170gmLC+Y6c6YUgkwfEs6dBvuTkL3o7+H8UfB6/+dhjlvgVVP7r3aieFB6l9FNZz/SZjVQMu1X6fpjkCGHUB2wvLi1QD+DklFzvAtafA2rU5WbHj4VzBmMlz4rzDvXUmIYJ9l1iQdmCNm0dxx6e7gvUums5vmJWv8nZKKnOFb0oD6rNW9YCYN0zJwy+fhgcXJDpWv+DSc/l6oHpt2uVJJaGvf0c/x7fDgDfDi10JZec7XSCpshm9J/dq1VveuJQNb2zM0Xb8cyr9Fw6gVcPbfJxvl1NSmXKlUWupqa2jNsWJQXVk7/OQDUHt08gfvS96W/AEsqWi42omkfuVcq7unguayd8OHVsIr/tHgLQ2DxgUzqancu2e7prKcxjeeCwt/CIdNg19/Cv5tFvzy47DlkZQqlXSw7PmWlNv2LclX3PTdhbJt5ygYOzn/NUkjxK5x3bknLE+DE18HbffBnd+Ee74Dd10DM18DZ/xNsomVu8dKBSvEGNOuIW/mzZsXly9fnnYZUuHq2gkPL4FVP4aHlzB/+9W00jdk19fWsOwTF6RQoKQ+tm2Ce74Ny78L2zfDkafAGe+Hk9+4e/KzpPwKIdwbY5yX65w939JIF2OyVOCq6+D+xcl212OPgJe9l8bqmTQt3brX0BPX6pYKzLgj4YJ/gLM/Cquvhzu/SctP/5vmn3TT1jOBunHlNF54Mg2nTUu7UkkYvqURo8+qJWdNoKFrCay8Dp5dDxU1cOLFMOcymH4elFfQAHBYjtVOXOpMKjyVNfCSt9PCBTQtXkmmJzncuq0nmSi98loazpkHx5wDFVXp1iqNYA47kUaAfVctAahhJ4sqvk3D8eUw5/JkDGn1uBSrlDQU5l91c86VUurDMyyr/luoHg/HvypZrvC4VyVb3EsaUg47kUaq7k5ovZfmnz9BpnPvnq4M1TSP/jAN73h1SsVJGg5tOYI3QFucBG/5CTx0Izx0E9z/MyirhBnnJkF85muSISyShpXhWyolMcJTf4L1t8Ijv4NHb4eO52nbcW3Ol7c915nf+iQNu37XCK+tgRMugBMWwMVfgY33JEH8TzfCjR9Jbke9lJaxl9H8yDTanutyqJk0DAzfUhHJudvkjAiPZMP2+lvh+SeTF0+cAbP/AmacR93Pq2nd2tHn36urrclr/ZKGX+OCmX2Hme07UbqsHI4+I7m96l/g6YfgoRtpuWcdTWsnkqELyG6s9dM/QvtjNJxzOpQbG6RD5W+RVCRy7jb5k7uT3SYr/gCjJ8GM87K3c5Md8LIad+YY8+2qJVJJGniN8BxCgCknwpQTaf7DzWTYu9c8011G85KHabjjzVB/WhLYp70Mjnqp48WlQXDCpVTItm+BJ1bCE/cxf0k9rZ1j+rykvqaLZe89FqbMgrL+N63N2WvuV8mSepn+iV+QKxUEIuvPvgU23AVPPgCxB0IZTDkJpp2+O5C3rC+n+dcP+z6jEc8Jl1IB6TcEv/AMPHFfsmvdE/clobv98d3XtXX2M247UwFHnrzfn9swt94PQUkD6n+8+Gh47ZeSJzueg9bl8PhdsOHOZFOu5d+hpetMmrquIEMyubu1PUPTz1ZCT49rjEu92PMt5VEydGQVmc6e3cdqQheLxv6Ehs4b97xw4gyYOgemnpq9n8P8r92Xe/kwd5uUNERyLktaWc6iS0/p/4/3nm548gHmf2s9rZm+fXr1YTPLjv1BdmjLScn9ESfBmEnD9Z8hpc6eb2mIHdAQjp3bYPNaeGZtcr/5zzSveCWZ7sP2elkmVtC8s4GGC89IwvaRp0BNbZ+feUCTqCTpEBz0eHFIJm9OnU1bZkPO021xIlSMSlZV+eMP9pwYMzkZLjdlFkw5kZb2GTTfvZO2rTsdsqKSZs+3RqzBjoHO2TNUDovmPEVDzcps4P4zPL+p11UBao9m+qYvEAl9/s0ArL/qtcNWsyQNt34399n17VyM8PxT8NSDyZKoTz2Qvf8TLZlTaep6Dxmqd19XU9bFopPaaDh5Akw4JrmNPSKZINqL74sqRAP1fBu+NSId8FernTvguVbYugHaN8DWjcy/5ThaO0b3+TfreZplEz4Lhx8Pk46Hw4/d83jCdKgctf8PJ0kqUoMasgLQ08P8q35D63NdfU7V8zTLRn1oz4GKmj1BfOJ0Wp4/kaYVh5PZ8yMP7GdKw8xhJypZg+3xaF6yZq8PCIBMZzfNLXfS8HBLEra3boQXnt7nykBbxw9z/pttTIaPPzrgz3XoiKRSNaghKwBlZbTlCN6QfV+98l549lF4dn32/lHYsh7W30bzti+wb3dGprOb5sW307BhGYyvg3FT976vmbC799xec6XB8K3UDdXwj9b2DE2LV0HnCzRMJ9ls5vknk685dz9Onre1N0GO4R9tO6uTzSYOOwqOnA2HTUse77qNr6fu6tv73z1uPwb94SRJRWCwqyoNuCvnpOOS275ipK3pppz/XlvnGFhzU44OFJLx5+OOpCWeR9NTF5DpSaJQsjrLfdD+OA2nHZNMCC2v7Ldmg7sGy/CtvRzKm8lgrs0doFcD2aDa3ZmsdZ3ZAts3J4+3b4bMFpp/ewyZzlF7/XuZzh6aF/+Bht5fUwKUVyVjBcdOgdqjqdu0g9adfcNyXW0NXHnPgDUfau+1S/5J0t4G9b4awsBLIzauha6OZP7Nc0/Atra97pvvO2N38N4l05V0jjTcemFyYFRtMjF0zKTsbTKMmUzL5mk0rZhApjvpxNnd+QMH3HlkcB+5HPNdgoZ0IuEBjp0b8NrZR8DO52BHO+zYutdt/v+NpzXTt2ehvuxZlo1rSq7rx/Qd1/YzeTGy/i2ZJGjvCty9vmY81P/WXdf7xilJQ2coOnDgwN/LB9xQ6E1PJXsvvPB09vbMnufbNzN/x1doZXKfa+vDMyybtCgJ7TUTsrfaPfejamnZNImmu6r2GadexqJLZw/754+fXfnjhMusfIfvNH5BDuWNqN/JgOMrWPbuumTpvJ3PQ8eu++ezx7Yx/w9zc/Yk14fNLKv+YL8/c/qOHxLpuytjILL+3N/D6MOTN63RE7OPs/ejJzL/S3cc0uRF34QkqfgN9r180BPge7qZ/slf9R/cz1wCmWch0569PZvcOl9Ifu6Or/Yf3Cf8E1SPh1Hjk/vq8VA9Lvt8HC1P19G0+ggy3Xs+N2sqYNGrj6Th1KlQNQYqx0B534ENh9rBZuA/OIbvrHyG72HrRd51bYzJkIzO7dC1Azoz0Jlh/rcfo3VbT59/s76mk2XnPgQdLyTXdGxP3gg6tmePvcD09R/rpye5h/Wj/rL/gitHM33bf/bfC33hAzDqsJy3+d9cQ+vWnX3rPYAAfai915KkkWtYOqsG+uzq6oAd7Uz/3D39B/ezliadWjuey3ZuPZd9nNzmZ76cO7j3WRVmFFSNTcJ49n7+o39Fa9e4vtfWdLJswSaorMmG95rsbTRUjqblz500/fpJMl17qh7SPLOf64s19Bu+s/IZvvv9xRzdnTTyrh3QtTN73/vxTubf9ypaO8f0vbbsWZbVfiYJ2l0ZiH1Ddv89ydkAXTEq+YWqGpO9H538lVw1mvkPvTn3zx3dw7I3RageC1Xjkr/Cq8dmf6HHQnnFoHsQHP4hSUpLQX3LvL9OpxiZ3nTTAENlnk460zpeSL6Z7nh+r+fTH3zPoDrY+u2pL9vMssO/kM0VNcl9xSioHJUsCVk5ivmrLsq9NG9NF8teuzl5fXkVVFT3ejwKKqpoebiDpt9u3if0l7Go4WQaTpvW//9OFEbnnEsNpqAtxy8WQNv2AL9s7HUkZBtttuFVVNPWeUnua3tqYdbrk9C8q7HvumUbet3iMlpf6Htt3WE18PEtyU5k/Wjsp7E2vu5UmDVwYx3sJMRDXf3DyYuSpMEa7GfIoXx2DXrS/v4mmM5714CX1/UT+usOq4EPre/1rfj27LfpyeO27/cN7ABtPRPhuFfu+fZ9Vyfijueg62noytDWcWnuazNlcONHBqy3ecdXyewT+jOdPTRfv5SGG/8uyU3llVBenQ3tVcl9eRXNG99Ppmv8Ptd2J5NpCyAzGL6HSb+/IOOr4UPr9oTtsoo+u3X1+wtSOxou/vKAP7exo58AfeGJAwZvOLQ3k0O9thB+GSRJOlBFFdwHuvbCE5N5VUzMeV1d7QCZ5JL/GPBn9h/4R8GVa7Lf+Hck990dvUYC7KTte32/3QdoYxKc9eHkdd2d0L3rvmP3sbYcw2ug/47RfCvY8B1CuBD4KlAOfDvGeNU+56uBHwCnAZuBhTHGR/NdZ3/6beQXnZQsVzSYaw/glyvNnmRDtCRJ+5dGcB/stcMT+GfBuCMHvHbA0P+KTw987eP9Xbv//TjyoSDHfIcQyoGHgVcBG4F7gMtjjA/2es37gdkxxveFEC4D3hBjXDjQvzsSVjuRJEkaSsW2eluhj/ku1PD9cuCzMcYF2edNADHGRb1esyT7mjtCCBXAJmByHOA/aKSs8y1JklTsirkTsxgnXNYDG3o93wi8rL/XxBi7QghbgcOBZ3q/KIRwBXAFwNFHHz1c9UqSJGkIlepQ2L5r0pWYGOO3YozzYozzJk/uu1SOJEmSlC+FGr5bgd6LOB6VPZbzNdlhJ4eRTLyUJEmSClKhhu97gONDCNNDCFXAZcAN+7zmBuAd2cdvAm4eaLy3JEmSlLaCHPOdHcN9JbCEZKnB78YYHwgh/DOwPMZ4A/Ad4L9DCGuBLSQBXZIkSSpYBRm+AWKMNwE37XPs070e7wDenO+6JEmSpMEq1GEnkiRJUskxfEuSJEl5YviWJEmS8sTwLUmSJOWJ4VuSJEnKE8O3JEmSlCeGb0mSJClPDN+SJElSnhi+JUmSpDwxfEuSJEl5YviWJEmS8sTwLUmSJOVJiDGmXUPehBCeBh5Lu448mgQ8k3YRKhm2Jw0125SGku1JQ+1Q2tSLYoyTc50YUeF7pAkhLI8xzku7DpUG25OGmm1KQ8n2pKE2XG3KYSeSJElSnhi+JUmSpDwxfJe2b6VdgEqK7UlDzTaloWR70lAbljblmG9JkiQpT+z5liRJkvLE8C1JkiTlieG7iIUQvhtCeCqEcH+vY80hhIdCCKtCCP8bQqjtda4phLA2hLAmhLAgnapVqPppT/+SbUv3hRB+HUKoyx4PIYSvZdvTqhDCS9KrXIUqV5vqde6jIYQYQpiUfW6b0oD6eY/6bAihNfsedV8I4TW9zvmZp3719/4UQvhgNkc9EEL4Yq/jQ9aeDN/F7XuJMVpDAAAD0klEQVTAhfsc+w1wcoxxNvAw0AQQQpgFXAaclL3mGyGE8vyVqiLwPfq2p+YY4+wY46nAjcCns8cvAo7P3q4AvpmvIlVUvkffNkUIYRrwauDxXodtU9qf75GjPQFfjjGemr3dBH7m6YB8j33aUwjhfOASYE6M8STg6uzxIW1Phu8iFmO8Ddiyz7Ffxxi7sk/vBI7KPr4EuC7GuDPGuB5YC5yet2JV8PppT8/1ejoG2DVD+xLgBzFxJ1AbQpian0pVLHK1qawvAx9jT3sC25T2Y4D2lIufeRpQP+3pb4CrYow7s695Knt8SNuT4bu0vRv4ZfZxPbCh17mN2WPSgEIInw8hbADeyp6eb9uTBiWEcAnQGmNcuc8p25QG68rsUKXvhhAmZI/ZnjQYJwBnhxDuCiHcGkJ4afb4kLYnw3eJCiF8CugCrk27FhW3GOOnYozTSNrSlWnXo+IVQhgNfJI9f8RJh+qbwLHAqcATwJfSLUdFrgKYCJwBNAI/CSGEof4hhu8SFEJ4J3Ax8Na4ZyH3VmBar5cdlT0mHahrgTdmH9ueNBjHAtOBlSGER0nazR9DCEdim9IgxBifjDF2xxh7gP9kz1AA25MGYyOwODv87W6gB5jEELcnw3eJCSFcSDKW8vUxxu29Tt0AXBZCqA4hTCeZ1HR3GjWqeIQQju/19BLgoezjG4C3Z1eoOAPYGmN8Iu8FqqjEGFfHGKfEGI+JMR5D8kH3khjjJmxTGoR95gW8Adi1coWfeRqMFuB8gBDCCUAV8AxD3J4qhqBQpSSE8CPgPGBSCGEj8BmS1U2qgd9kvym5M8b4vhjjAyGEnwAPkgxH+UCMsTudylWI+mlPrwkhzCT56/8x4H3Zl98EvIZk0sl24F15L1gFL1ebijF+p5+X26Y0oH7eo84LIZxKMnn3UeC9AH7maX/6aU/fBb6bXX6wA3hHdgTBkLYnt5eXJEmS8sRhJ5IkSVKeGL4lSZKkPDF8S5IkSXli+JYkSZLyxPAtSZIk5YnhW5IkScoTw7ckSZKUJ4ZvSVIfIYSXhhBWhRBGhRDGhBAeCCGcnHZdklTs3GRHkpRTCOFzwCigBtgYY1yUckmSVPQM35KknEIIVcA9wA7gTLfnlqRD57ATSVJ/DgfGAuNIesAlSYfInm9JUk4hhBuA64DpwNQY45UplyRJRa8i7QIkSYUnhPB2oDPG+D8hhHLgDyGEC2KMN6ddmyQVM3u+JUmSpDxxzLckSZKUJ4ZvSZIkKU8M35IkSVKeGL4lSZKkPDF8S5IkSXli+JYkSZLyxPAtSZIk5cn/B1Zl7uvvmZuBAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "print(fit_result.fit_report())\n", "plt.rcParams[\"figure.figsize\"] = (12,10)\n", "fit_result.plot()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# export simulation result" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [], "source": [ "f_list = np.arange(115,165,0.1)\n", "\n", "result_I0 = fit_result.params['I_0'].value\n", "result_x0 = fit_result.params['x_0'].value\n", "result_gamma = fit_result.params['gamma'].value\n", "result_bg = fit_result.params['bg'].value\n", "\n", "fit_result = cavity(f_list, result_x0, result_I0, result_gamma, result_bg)" ] }, { "cell_type": "code", "execution_count": 29, "metadata": {}, "outputs": [], "source": [ "out_text = np.c_[(f_list, fit_result)]\n", "file_name = 'cavity_laser_fit.dat'\n", "np.savetxt(file_name, out_text, delimiter='\\t', fmt='%f', header=\"freq(MHz)\\ttransmission\")" ] }, { "cell_type": "code", "execution_count": 30, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 30, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsIAAAI/CAYAAAB9Hr8eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8GearUAAAgAElEQVR4nOzdeZScZ3Xv+99TY3dXV8+T1C1Zs215NrIwNjYm2ImBBBNOmHISSOBCBkhuhpVckpORnHNz74FkXcghIWRYkBGIDwEHTMxkMDY2tmzZkmVbssZWd6vHqp6qu6q6qp77R/Xbaskaeqiq9616v5+1tKyuLnVteXlV//z0fvY21loBAAAAfhNwuwAAAADADQRhAAAA+BJBGAAAAL5EEAYAAIAvEYQBAADgSwRhAAAA+FLIrRfu6OiwW7ZscevlAQAA4BNPP/30uLW28/zHXQvCW7Zs0b59+9x6eQAAAPiEMebUhR6nNQIAAAC+RBAGAACALxGEAQAA4EsEYQAAAPgSQRgAAAC+RBAGAACALxGEAQAA4EsEYQAAAPgSQRgAAAC+RBAGAACALxGEAQAA4EsEYQAAAPgSQRgAAAC+RBAGAACALxGEAQAA4EsEYQAAAPgSQRgAAAC+RBAGAACALxGEAQAA4EsEYQAAAPgSQRgAAAC+RBAGAACALxGEAWCN0gt5/dk3Duttf/mYnjyRcLscAMAqEYQBYI3+7tET+ovvHNXh4Rn9/pefV75g3S4JALAKBGEAWIN8wepfftiv23e062Nvv0GHR2b0b/tOu10WAGAVCMIAsAbfeWlUg5Pz+tlbr9Abr+3Rru5GffnZQbfLAgCsAkEYANbgP58fVmtDWHdf3S1jjF63q1PPnJrUfDbvdmkAgBUiCAPAGjzTn9SeLW0KBYtvo7dt71A2X9DTp5IuVwYAWCmCMACs0sRsRifGU3rVFa1Lj+3d2qZQwOixY+MuVgYAWA2CMACs0v7+SUk6JwjHoiHduKlFjx+bcKssAMAqEYQBYJWe7k8qHDS6rrf5nMdv2NSil4anGaMGAFXiskHYGPP3xphRY8zzF/m8McZ80hhz1BhzwBhzc+nLBADv2N+f1O6NzaoLB895/MruuNILBZ1OzLlUGQBgNVZyIvxZSfde4vNvlLRz8dcHJf3V+ssCAO96eWRWV/fEX/H4lYuPvTQ8U+mSAABrcNkgbK19RNKldofeJ+kfbNETklqMMRtKVSAAeEkyldVEKqvtnY2v+NzO7uJjR0YIwgBQDUrRI9wrafk6pYHFxwCg5hwbm5Ukbe+KveJzDZGQNrc16DBBGACqQkUvyxljPmiM2WeM2Tc2NlbJlwaAknCC8I7OV7ZGSMX2iCO0RgBAVShFEB6UtGnZx32Lj72CtfYz1to91to9nZ2dJXhpAKiso6OzioQC6m2tv+Dnr+yO6/h4SpkcG+YAwOtKEYQfkPSexekRt0qastaeKcHXBQDPOTaW0raOmIIBc8HPb+uMKV+wGkzOV7gyAMBqhS73BGPMv0q6S1KHMWZA0h9KCkuStfbTkh6U9CZJRyXNSfr5chULAG47Ojqr6/qaL/r5zW0NkqT+xJy2XeBCHQDAOy4bhK21777M562kD5WsIgDwqGyuoIHknO67ceNFn+MEYWYJA4D3sVkOAFbozNS8ClbatBh2L6QzHlU0FFA/QRgAPI8gDAArNLDY99t3kYtykmSM0aa2Bp1O0CMMAF5HEAaAFXIuwG1qvfiJsFRsj+BEGAC8jyAMACs0kJxTwEg9zXWXfN7mtgadTsypeIUCAOBVBGEAWKGB5Lw2NNcrHLz0W+emtgbNZHKanFuoUGUAgLUgCAPACg0k59XbcvH+YIczOeIU7REA4GkEYQBYoYHk3CUvyjmc5wxNcmEOALyMIAwAK5DNFTQ8nV5REN6w2EN8Zipd7rIAAOtAEAaAFRieSqtgpb7LTIyQpOb6sOrCAQ1PcSIMAF5GEAaAFRhcbHPYuIIeYWOMNjTXcyIMAB5HEAaAFRiZLobanuboip7f01SnYYIwAHgaQRgAVsAJwt1Nl54h7NjQXMeJMAB4HEEYAFZgeDqtWCSoeF14Rc/vaa7TyHRa+QJLNQDAqwjCALACI9NpdV9mo9xyG1rqlStYTcxmylgVAGA9CMIAsALDU2n1rLAtQpI2NDFCDQC8jiAMACswMp1ZVRDuYZYwAHgeQRgALqNQsBqZTqtrNSfCi0GYWcIA4F0EYQC4jMRcVrmCVU/TykanSVJbLKJw0OjMNCfCAOBVBGEAuAxnHnDPKi7LGWPU2RjV2AyX5QDAqwjCAHAZq50h7OhsqiMIA4CHEYQB4DKGp1d/IiyJE2EA8DiCMABchhNmOxpX3iMsSV1NBGEA8DKCMABcxthMZvHy2+reMrviUU2kslrIF8pUGQBgPQjCAHAZYzMZda7yNFiSOuPFPzPOdjkA8CSCMABcxvhsRh3xyKr/XFe82FNMewQAeBNBGAAuY2x2fSfCo9MEYQDwIoIwAFyCtbbYGhFffRDuWvwzY7RGAIAnEYQB4BJmMzmlFwprCsLOlAlOhAHAmwjCAHAJax2dJkmRUECtDWGNzrBmGQC8iCAMAJcwPpuVpDWdCDt/bpTLcgDgSQRhALgE50R4rUG4ozGqCXqEAcCTCMIAcAlji20Na5kaIUntjVElUtlSlgQAKBGCMABcwthsRsGAUWvD6ucIS1J7LKIJgjAAeBJBGAAuwVmvHAiYNf359lhEM+mcMrl8iSsDAKwXQRgALmFiNrumiRGOtsbiSXIytVCqkgAAJUIQBoBLmEhl1R5bW1uEJLXHiiF6nAtzAOA5BGEAuIREKqu29QThxRNhLswBgPcQhAHgEtYdhBf/7ESKE2EA8BqCMABcRCaX12wmp47G9bdGTMxyIgwAXkMQBoCLcNoZ2mJrvyzXVB9SKGAYoQYAHkQQBoCLcE5x19MaYYxRWyyiBCfCAOA5BGEAuAjnRLh9Ha0RxT8fpUcYADyIIAwAF3G2NWKdQZjtcgDgSQRhALgIJ7yuZ46wVDxRZnwaAHgPQRgALiKRyigYMGqqC6/r67TFIkyNAAAPIggDwEUkUlm1NkQUCJh1fZ2OxqhmMzmlF/IlqgwAUAoEYQC4iInZ9a1Xdjg9xrRHAIC3EIQB4CLWu1XO0U4QBgBPIggDwEUkUlm1rXN0mnR2/Nr4LCPUAMBLCMIAcBETqdK0RjhrljkRBgBvIQgDwAUs5Auaml9YCrHr4ZwqMzkCALyFIAwAF5B0lmmUoDUiHg0pEgywVAMAPIYgDAAXUKplGpJkjFmcJUyPMAB4CUEYAC6gVOuVHWyXAwDvIQgDwAWU8kRYWtwuRxAGAE8hCAPABSQW2xhKdSLc0RjVRIrWCADwEoIwAFxAIpWVMVJLQ+lOhBNMjQAATyEIA8AFTKSyam2IKBgwJfl67Y0RpbJ5pRfyJfl6AID1IwgDwAWUar2yw+k1ZrscAHgHQRgALmCixEG4dbHFYnJuoWRfEwCwPgRhALiARInWKzucUM0INQDwDoIwAFxAqVsjnEt3yTmCMAB4BUEYAM6TL1gl58pzIpzkRBgAPIMgDADnmZzLytrSzRCWpOb6sIyRkvQIA4BnEIQB4DxOWG0tYRAOBoya6sKapDUCADyDIAwA53HCamuJlmk42mIRJTgRBgDPIAgDwHmcyQ6lbI2QpJYGToQBwEsIwgBwHmfWb0tDuKRft7UhwtQIAPAQgjAAnCdZptaI1oaIkilaIwDAKwjCAHCe5NyCIsGAGiLBkn7d1oYwJ8IA4CEEYQA4TzKVVWssLGNMSb9uayyiuWxe6YV8Sb8uAGBtCMIAcJ7kXLbkbRHS2VaLSSZHAIAnEIQB4DyTcwslvygnFVsjJNYsA4BXEIQB4DyJMp0ItzSwZhkAvIQgDADnmZzLlnSrnMOZS8yaZQDwBoIwACxjrdXk3MJSG0Mp0RoBAN5CEAaAZWYyOeUKltYIAPABgjAALOOE1JYyBOFIKKDGaIjWCADwCIIwACzjhNRytEZIxbXNk7RGAIAnEIQBYJml9cpluCwnFWcJJwjCAOAJBGEAWMY5rS1Hj7BUDNi0RgCANxCEAWCZRKq8rRGttEYAgGcQhAFgmcm5rAJGaqorVxCOKMHUCADwBIIwACyTnMuqpSGiQMCU5eu3NkQ0k85pIV8oy9cHAKwcQRgAlknOLailTG0RktQaK37tSfqEAcB1BGEAWCaZypbtopx0dj4xfcIA4D6CMAAskyzTemVHm7NdjhNhAHAdQRgAlpmcK/eJcDFkc2EOANy3oiBsjLnXGHPYGHPUGPORC3x+szHmYWPMfmPMAWPMm0pfKgCUX3IuW7ZlGtLZRR20RgCA+y4bhI0xQUmfkvRGSbslvdsYs/u8p/2epC9aa2+S9C5Jf1nqQgGg3OazeaUXCmW9LEdrBAB4x0pOhPdKOmqtPW6tzUr6vKT7znuOldS0+PtmSUOlKxEAKiNZ5q1yklQfCSoaCiy9FgDAPaEVPKdX0ullHw9IevV5z/kjSd8wxvyKpJiku0tSHQBUUCWCsCS1xSJK0iMMAK4r1WW5d0v6rLW2T9KbJP2jMeYVX9sY80FjzD5jzL6xsbESvTQAlIYz27ecUyOk4gg1WiMAwH0rCcKDkjYt+7hv8bHl3i/pi5JkrX1cUp2kjvO/kLX2M9baPdbaPZ2dnWurGADKxJnkUM7LclIxaNMaAQDuW0kQfkrSTmPMVmNMRMXLcA+c95x+SW+QJGPM1SoGYY58AVQVZ5JDOS/LScWgTRAGAPddNghba3OSPizpIUkvqjgd4pAx5qPGmLcsPu03JX3AGPOcpH+V9HPWWluuogGgHJJLrREVOBGmRxgAXLeSy3Ky1j4o6cHzHvuDZb9/QdLtpS0NACorOZdVPBpSOFjeXUOtDRFNzS+oULAKBExZXwsAcHFslgOARclUVi2x8rZFSMUgXLDSdJoLcwDgJoIwACxKzi2UvS1CklpjrFkGAC8gCAPAosm5bEWCcAvb5QDAEwjCALCoeCJc/tYIZ83yJJMjAMBVBGEAWJRMZZdOa8vJOXWmNQIA3EUQBgBJC/mCZjK5yrRGLPYIT9IaAQCuIggDgJatV67A1Ih4NKRQwLBUAwBcRhAGAGkplFbiRNgYo5YGtssBgNsIwgAgLW16q0QQLr5OWMkUrREA4CaCMADo7CizlgpMjZCk1lhECU6EAcBVBGEA0NlRZq2xyp0IMz4NANxFEAYAael0tq1irRERFmoAgMsIwgCg4tSIaCig+kiwIq/XGosomcrKWluR1wMAvBJBGABUvCxXqYtyUrE1Ilewms3kKvaaAIBzEYQBQMXLcpW6KCdpaYMdSzUAwD0EYQBQcY5wZU+EI0uvCwBwB0EYAFQMpG0VmhghSW2LG+wSKYIwALiFIAwAKrYoVGK9soPWCABwH0EYgO8VClaTtEYAgO8QhAH43nR6QQV79pS2EprrwzLm7GpnAEDlEYQB+J6z2KKtgq0RwYBRc32YpRoA4CKCMADfcy6sVfJEWHK2y3EiDABuIQgD8L3JxTBayR5hSWppCHNZDgBcRBAG4HvOiXBbhYNwW0OE8WkA4CKCMADfc05lWyrYIywVWzEmaY0AANcQhAH4XnIuq1DAKB4NVfR1Wxu4LAcAbiIIA/C95FxWLQ0RGWMq+rqtsYjmF/JKL+Qr+roAgCKCMADfS6YW1NpQ2bYIiaUaAOA2gjAA30vMZdUaq+xFOUlL4TuZoj0CANxAEAbge8X1yi6cCC+Gby7MAYA7CMIAfC85t1DxGcLS2daIBEEYAFxBEAbga9ZaJVMut0YwOQIAXEEQBuBrs5mccgXrSmuEs9J5kqUaAOAKgjAAX3OWabjRGhEJBdQYDdEaAQAuIQgD8DVnxbEbQViSWhrCS2EcAFBZBGEAvubM8G2t8HplR2tDhDnCAOASgjAAX1sKwi6dCLfGIlyWAwCXEIQB+JqzzMK1INwQVpLLcgDgCoIwAF+bnMsqYKSmelojAMBvCMIAfC0xl1VzfVjBgHHl9VsawppJ57SQL7jy+gDgZwRhAL7m1lY5R9vSmmX6hAGg0gjCAHzNra1yjqWlGrRHAEDFEYQB+FrxRNid/mCJNcsA4CaCMABfm5zLLp3KusFpy+DCHABUHkEYgK8lUtmlPl03OG0ZjFADgMojCAPwrflsXplcQS20RgCALxGEAfiW047Q5mJrRH04qEgowGU5AHABQRiAbyUW2xHc7BE2xqitIbJUCwCgcgjCAHzLmd3r5tQIqbhUg9YIAKg8gjAA30o4rREuXpaTipMjaI0AgMojCAPwLSd8utkaIRWDeIIgDAAVRxAG4FvJVLEdwc2pEc7rs2IZACqPIAzAt5JzWcXrQgoH3X0rdFojCgXrah0A4DcEYQC+lZzLLm12c1NLQ1gFK82kc26XAgC+QhAG4FuJVHZps5ubnMt69AkDQGURhAH41uTcguuj0yQtnUonCcIAUFEEYQC+5aXWCEmMUAOACiMIA/CtZMobQXipNSLF5AgAqCSCMABfyuYKSmXznmiNcOYYcyIMAJVFEAbgS07o9MJluaa6kIIBo0SKIAwAlUQQBuBLzoQGL7RGGGPU2hDhshwAVBhBGIAvOVvlvNAaIUntsYgmZgnCAFBJBGEAvpT0UGuEJLXGwpwIA0CFEYQB+FLSQ60RktQei2qCHmEAqCiCMABfSsw6J8LeaI1ojYWVJAgDQEURhAH4UmIuq8ZoSNFQ0O1SJEltsagm5xeUL1i3SwEA3yAIA/ClRCq7tMjCC9oawrKWWcIAUEkEYQC+5Lkg3BiVJGYJA0AFEYQB+NLEbFbtXgrCDc6aZYIwAFQKQRiALyXnsp4ZnSZp6XSaIAwAlUMQBuA71lpNpLx1ItzeuBiE6REGgIohCAPwnVQ2r2yu4Kke4ZbFDXcJtssBQMUQhAH4ztkZwt4JwtFQUPFoiKUaAFBBBGEAvuO0H3ipNUIqBnPWLANA5RCEAfhOIpWRJE+1RkjFergsBwCVQxAG4DsTs86JcNTlSs5FEAaAyiIIA/Adp/2gNRZ2uZJzEYQBoLIIwgB8ZyKVVSQYUGM05HYp53CCsLXW7VIAwBcIwgB8JzFbXK9sjHG7lHO0xSLK5Aqay+bdLgUAfIEgDMB3EilvbZVzsGYZACqLIAzAdxJz3toq52DNMgBUFkEYgO8kUlnPjU6Tzi74YM0yAFQGQRiA7zg9wl7jnFKzZhkAKoMgDMBXMrm8ZjI5TwZh50SY7XIAUBkEYQC+Mjm3IMl7W+UkqakupFDAaIIeYQCoCIIwAF85u1XOe0HYGKPWWERJgjAAVARBGICvOBMZvDg+TSoGdE6EAaAyCMIAfMWZyODFE2FJam3gRBgAKmVFQdgYc68x5rAx5qgx5iMXec47jDEvGGMOGWP+pbRlAkBpJGYzkrzZIyxJbY0R5ggDQIWELvcEY0xQ0qck3SNpQNJTxpgHrLUvLHvOTkm/I+l2a23SGNNVroIBYD0SqayMkVoavBmEaY0AgMpZyYnwXklHrbXHrbVZSZ+XdN95z/mApE9Za5OSZK0dLW2ZAFAaE6msWurDCgaM26VcUHssqqn5BS3kC26XAgA1byVBuFfS6WUfDyw+ttwuSbuMMY8ZY54wxtxbqgIBoJSSc95cpuFob2TNMgBUSqkuy4Uk7ZR0l6R3S/obY0zL+U8yxnzQGLPPGLNvbGysRC8NACs3MZtVeyzqdhkX1bEYhMcXe5kBAOWzkiA8KGnTso/7Fh9bbkDSA9baBWvtCUlHVAzG57DWfsZau8dau6ezs3OtNQPAmiVSWbXGwm6XcVHtjcWQPsGaZQAou5UE4ack7TTGbDXGRCS9S9ID5z3nyyqeBssY06Fiq8TxEtYJACUxkcqqo9G7J8LOWLeJFCfCAFBulw3C1tqcpA9LekjSi5K+aK09ZIz5qDHmLYtPe0jShDHmBUkPS/ota+1EuYoGgLXI5QtKzmWXTl29iBNhAKicy45PkyRr7YOSHjzvsT9Y9nsr6TcWfwGAJyXmsrL2bB+uFzXVhRQJBjROEAaAsmOzHADfcE5ZvdwaYYxRe2NEE1yWA4CyIwgD8A0nCHt1vbKjvTHC1AgAqACCMADfcMJlR9y7J8JScakG2+UAoPwIwgB8YykIe3iOsKTF1giCMACUG0EYgG+Mz2YVDho11a/onrBrOhqjGp/NqHgPGQBQLgRhAL4xMZtReywqY4zbpVxSR2NEmVxBqWze7VIAoKYRhAH4xvhsRh1xb1+Uk7S0AprJEQBQXgRhAL4xkcouhUwva1+cc8wsYQAoL4IwAN+YmPX2emWHUyMj1ACgvAjCAHzBWqux2Yynt8o5nBNhJkcAQHkRhAH4wmwmp2yusBQyvawt5gRhToQBoJwIwgB8YbwK1is7oqGgmupCLNUAgDIjCAPwBed0tb0KgrB0dpYwAKB8CMIAfGFpq1wVtEZIbJcDgEogCAPwhWpqjZCKs4QnUpwIA0A5EYQB+IJzuupcRPM6ToQBoPwIwgB8YXw2o5aGsMLB6njba2+MKjGXVS5fcLsUAKhZ1fEdAQDWaSKVqZq2CKnYy2ytlJxbcLsUAKhZBGEAvjA+k1V7lbRFSGd7mekTBoDyIQgD8IXxKjsRbo+xXQ4Ayo0gDMAXxmeqY72yw5l3zCxhACgfgjCAmpfNFTSdzlXNMg1J6owXax2bIQgDQLkQhAHUvESqumYIS1JTXUiRUIAgDABlRBAGUPPGl9YrV09rhDFGXfEoQRgAyoggDKDmjVXZemVHZzyqUYIwAJQNQRhAzXNOVbvidS5XsjqdjZwIA0A5EYQB1DwnTDoX0KpFV1N06TQbAFB6BGEANW9sJqN4XUh14aDbpaxKZ2OdEqmssjnWLANAORCEAdS80Zm0uqrsNFg6e4LNdjkAKA+CMICaNzqdqbr+YElL4Z0+YQAoD4IwgJo3Npupuv5g6eyJ8Og0QRgAyoEgDKCmWWsXT4SrNwhzYQ4AyoMgDKCmzWZyml/Iq6up+oKwswmP1ggAKA+CMICaNlqlo9MkKRIKqLUhrNGZtNulAEBNIggDqGnVukzD0cmaZQAoG4IwgJo2uhSEq+9EWCoGeIIwAJQHQRhATRudLrYVVGNrhFSse5QgDABlQRAGUNPGZjOKhAJqrg+7XcqaOK0R1lq3SwGAmkMQBlDTxqYz6myMyhjjdilr0hWPKpMraCaTc7sUAKg5BGEANW10pjqXaThYqgEA5UMQBlDTxmaqc5mGo5NZwgBQNgRhADVtdCZdlcs0HE7tzBIGgNIjCAOoWdlcQcm5BXU2VucMYUlLtXMiDAClRxAGULPGZxdnCFfxiXBTfUiRUEBjswRhACg1gjCAmlXtyzQkyRijzsaoxrgsBwAlRxAGULOqfZmGozMe5UQYAMqAIAygZjnhsStevT3CUvFEe2Say3IAUGoEYQA1a3Q6I2Ok9saI26WsS3dTnUZojQCAkiMIA6hZozMZtTVEFA5W91tdT3OdpuYXlF7Iu10KANSU6v7uAACXMFblW+Uc3U3F1o7hKdojAKCUCMIAatbIdHopRFazHicI0ycMACVFEAZQs85MpbWhuQaCcHPxVJsLcwBQWgRhADUpmytoIpWpiRNh5+9whtYIACgpgjCAmjQ6k5a1qokT4XhdWLFIkB5hACgxgjCAmuS0EXTXQBCWin8PWiMAoLQIwgBqktNGUAsnwlLxwhyX5QCgtAjCAGqS00bQUwM9wlLx7zFCawQAlBRBGEBNGp5Kqy4cUHN92O1SSqK7uU6jMxkVCtbtUgCgZhCEAdSk4em0eprqZIxxu5SS6GmqU65gNZ5i1TIAlApBGEBNGp5Kq6dG+oOlsyPURqYIwgBQKgRhADXJORGuFU6o58IcAJQOQRhAzSkUrEam0+pprne7lJJhzTIAlB5BGEDNScxltZC36mmKul1KyXQ0RhQwYnIEAJQQQRhAzVkanVZDJ8KhYECd8SgnwgBQQgRhADXnbBCunR5haXGWMEEYAEqGIAyg5jinprV0WU4qTo4YpjUCAEqGIAyg5gxPpRUMGHXGa6dHWCqecNMaAQClQxAGUHOGp9PqbIwqGKiNZRqO7qY6zaRzmsvm3C4FAGoCQRhAzam1ZRqOpRFqtEcAQEkQhAHUnFpbpuHY0EwQBoBSIggDqDm1eiK8saU4Dm6IIAwAJUEQBlBTZtILms3kajIIO3+nocl5lysBgNpAEAZQU0ZqdHSaJNWFg+pojGowSRAGgFIgCAOoKQOLIbG3tXa2yi3X21qvoSmCMACUAkEYQE0ZXGwb6G2p0SDcUrf0dwQArA9BGEBNGZqcVzBg1FVjyzQcG5vrNTQ5L2ut26UAQNUjCAOoKYPJefU01SkUrM23t40t9UovFJRIZd0uBQCqXm1+pwDgW4OT8zXbHyyd7X0emmSEGgCsF0EYQE0ZTM6rr0b7g6Wzvc/0CQPA+hGEAdSMXL6g4en00uKJWrS0VIMgDADrRhAGUDOGp9Mq2NodnSZJrQ1h1YUDBGEAKAGCMICa4SyaqNXRaZJkjFFvSz2tEQBQAgRhADVjaYZwDZ8IS8X2CE6EAWD9CMIAasZQjS/TcBRPhJkaAQDrRRAGUDMGJ+fVHouoLhx0u5Sy2thSr/HZjNILebdLAYCqRhAGUDMGkrU9Q9jhnHifmeJUGADWgyAMoGYMTs7XfFuExAg1ACgVgjCAmmCt1ZBPgjBLNQCgNFYUhI0x9xpjDhtjjhpjPnKJ5/0XY4w1xuwpXYkAcHmJVFbphUJNL9Nw9DTXyRhOhAFgvS4bhI0xQUmfkvRGSbslvdsYs/sCz4tL+j8l/bDURQLA5fhldJokRUIBdcWjS3OTAQBrs5IT4b2Sjlprj1trs5I+L+m+CzzvTyT9v5K4vQGg4vwyOs3R21KvAYIwAKzLSoJwr6TTyz4eWHxsiTHmZkmbrLVfK2FtALBiTijs88GJsCRtbmtQf2LO7TIAoKqt+7KcMSYg6c8l/eYKnvtBY8w+Y8y+sbGx9b40ACwZnJxXQySo5vqw26VUxOa2Bp2Zmlc2V3C7FACoWisJwoOSNi37uG/xMUdc0rWSvmuMOSnpVgW4QO0AACAASURBVEkPXOjCnLX2M9baPdbaPZ2dnWuvGgDOczoxr02tDTLGuF1KRWxqa1DBcmEOANZjJUH4KUk7jTFbjTERSe+S9IDzSWvtlLW2w1q7xVq7RdITkt5ird1XlooB4AL6Eyltbm9wu4yK2dxW/LvSHgEAa3fZIGytzUn6sKSHJL0o6YvW2kPGmI8aY95S7gIB4HKstepPzC2FQz9wQv/pJEEYANYqtJInWWsflPTgeY/9wUWee9f6ywKAlRubySi9UNAVPjoR7o7XKRIMcCIMAOvAZjkAVe/UYhj004lwIGDU11av0wRhAFgzgjCAqtc/4b8gLDFCDQDWiyAMoOqdSszJGKmv1YdBeIIgDABrRRAGUPVOJ+a0sblekZC/3tI2tzVoOp3T1NyC26UAQFXy13cNADXp1ETKd20RUnGWsMQINQBYK4IwgKrnt9FpDmYJA8D6EIQBVLVUJqfx2ayvlmk4OBEGgPUhCAOoak4I9NMMYUdjNKT2WIQgDABrRBAGUNX6fThDeLlNbQ3MEgaANSIIA6hqzviwK9piLlfiDmYJA8DaEYQBVLX+xJya6kJqbgi7XYorNrc1aHByXrl8we1SAKDqEIQBVLVTiTld0e7P02CpGITzBavByXm3SwGAqkMQBlDV+idSvpwY4djaWfyfgOPjKZcrAYDqQxAGULXyBauB5LxvL8pJ0taOYhA+MUYQBoDVIggDqFoDyTnlClZbfHwi3B6LqKkupBOcCAPAqhGEAVSt44unoNs6G12uxD3GGG3tbNTx8Vm3SwGAqkMQBlC1jo0Vw9+2Dv9elpOKf39aIwBg9QjCAKrW8fGUmuvDaotF3C7FVds6YhqaSms+m3e7FACoKgRhAFXr+NistnXGZIxxuxRXOZMj6BMGgNUhCAOoWsfHUtrW4d/+YMfS5AiCMACsCkEYQFWaSS9odCajbZ3+7g+WlgdhLswBwGoQhAFUJef0cztBWA2RkDY01y1N0QAArAxBGEBVYnTaubZ2xNguBwCrRBAGUJWOj80qYKQrfLxMY7mtHTEdH5uVtdbtUgCgahCEAVSl4+Mp9bU2KBoKul2KJ2zrbNR0OqdEKut2KQBQNQjCAKrS8bEUF+WW2cbkCABYNYIwgKpTKFidGGd02nLO5Aj6hAFg5QjCAKrO8HRa8wt5ToSX6WutVzholtZOAwAujyAMoOqcnRhBEHaEggFt7Yjp6AhBGABWiiAMoOocX1wcsZ3RaefY2R3XkdEZt8sAgKpBEAZQdV4emVU8GlJXPOp2KZ6yqyuu04l5zWVzbpcCAFWBIAyg6hwemdGunriMMW6X4im7uosn5EdHaY8AgJUgCAOoKtZaHRmZ0a7uuNuleM7OxX8nR+gTBoAVIQgDqCpjMxlNzi3oym76g8+3pb1BkWBAL4/QJwwAK0EQBlBVDi+GvF09nAifLxQMaFtnTEcIwgCwIgRhAFXl8HAx5F1Ja8QF7eqO0xoBACtEEAZQVY6MzKijMaL2RiZGXMiu7kYNTs5rNsPkCAC4HIIwgKpyeGSWi3KX4FyYo08YAC6PIAygahQKVi8zMeKSdi0FYdojAOByCMIAqsbg5LzmsnldyUW5i9rc1qBoKMCFOQBYAYIwgKrhXJTjRPjiggGjHV2NOsJSDQC4LIIwgKqxNDqNGcKXtKs7To8wAKwAQRhA1TgyMqPelnrF68Jul+Jpu7rjOjOV1uRc1u1SAMDTCMIAqsaRkVlOg1fgmo1NkqQXzky7XAkAeBtBGEBVyOYKOjo6w0a5FdjtBOEhgjAAXApBGEBVODIyo4W81bUbm90uxfM6GqPqborqEEEYAC6JIAygKjinm9f2EoRX4pqNzZwIA8BlEIQBVIXnh6bUGA3pirYGt0upCtdsbNLRsVmlF/JulwIAnkUQBlAVDg1Na/eGJgUCxu1SqsLuDU3KF+zS7GUAwCsRhAF4Xr5g9cLQ9NIlMFzeNYu91EyOAICLIwgD8LwT47OaX8jTH7wKm9rqFa8L6dDQlNulAIBnEYQBeN6hpYtynAivlDFGuzc0MTkCAC6BIAzA854fnFIkFND2TpZprMbujU166cyM8gXrdikA4EkEYQCed2hoWlf3xBUO8pa1GtdsbNb8Ql4nxlNulwIAnsR3FQCeZq3V84NT2s0ijVVzVi3TJwwAF0YQBuBpA8l5Tadz9AevwY6uRkVDAR0YIAgDwIUQhAF4mnOaeQ0nwqsWDgZ0bW+znj096XYpAOBJBGEAnrb/9KQiwYCu3hB3u5SqdOOmFj0/OKWFfMHtUgDAcwjCADxtf/+krt7YpGgo6HYpVenGTS3K5ApsmAOACyAIA/CsXL6ggwNTumlTi9ulVK0bF//d7ac9AgBegSAMwLMOj8xofiGvmzYThNeqr7Ve7bGInu0nCAPA+QjCADxr/2J4u3lzq8uVVC9jjG7c1KJnTyfdLgUAPIcgDMCz9vdPqj0WUV9rvdulVLUbN7Xo2FhKU/MLbpcCAJ5CEAbgWftPJ3XT5hYZY9wuparduNhacmCA9ggAWI4gDMCTpuYWdHwspZtoi1i36/uKQfg5LswBwDkIwgA86dnF00smRqxfc31Y2zpjLNYAgPMQhAF40v7+pIyRrutjo1wp3LSpVc/0T8pa63YpAOAZBGEAnrS/f1K7uuKK14XdLqUm3LKlVYlUVsfGUm6XAgCeQRAG4Dn5gtUzp5K6+Qr6g0tl79Y2SdKTJxIuVwIA3kEQBuA5LwxNayaT063b2twupWZs7YipozGqp04ShAHAQRAG4Dk/PDEh6ewpJtbPGKO9W1s5EQaAZQjCADznhycS2tzWoA3NLNIopb1b2jQ4Oa+B5JzbpQCAJxCEAXhKoWD11MmEXs1pcMnt3douSbRHAMAigjAATzkyOqPJuQW9elu726XUnCt74orXhWiPAIBFBGEAnvLD48WQxolw6QUDRrdsaSMIA8AigjAAT/nhiQltbK5TXyv9weWwd2ubjo2lND6bcbsUAHAdQRiAZ1hr9eSJhF69rV3GGLfLqUnOJA7n5B0A/IwgDMAziieVWdoiyui63mbFoyE9enTc7VIAwHUEYQCe8dhiOHvNdi7KlUs4GNCt29v16NExt0sBANcRhAF4xiNHxnRFe4OuaI+5XUpNu2Nnh04n5nVqIuV2KQDgKoIwAE/I5gp6/PiE7tjZ4XYpNe+1O4r/jh95mfYIAP5GEAbgCc/0JzWXzeuOnZ1ul1LztnbE1NtSr0dfpj0CgL8RhAF4wiNHxhQMGN1Gf3DZGWN0x84O/eDYhHL5gtvlAIBrCMIAPOH7L4/r5s0titeF3S7FF167s0Mz6ZwODE65XQoAuIYgDMB1E7MZPT80pTtpi6iY27d3yBjpUfqEAfgYQRiA6x49Oi5rpTt2EYQrpTUW0XW9zXrkCH3CAPyLIAzAdY8cGVdLQ1jX9Ta7XYqvvP7KLj3Tn1QilXW7FABwBUEYgKvyBauHD4/qzp2dCgZYq1xJd1/drYKVHn5p1O1SAMAVKwrCxph7jTGHjTFHjTEfucDnf8MY84Ix5oAx5tvGmCtKXyqAWuScSN6zu9vtUnzn2t4mdTdF9a0XR9wuBQBccdkgbIwJSvqUpDdK2i3p3caY3ec9bb+kPdba6yXdL+l/lrpQALXpmy+MKBw0uutK+oMrzRijN1zdrUeOjCmTy7tdDgBU3EpOhPdKOmqtPW6tzUr6vKT7lj/BWvuwtXZu8cMnJPWVtkwAtchaq2++MKJbt7UzNs0l91zdrVQ2ryeOJ9wuBQAqbiVBuFfS6WUfDyw+djHvl/T19RQFwB+OjaV0YjylH6UtwjWv2d6u+nBQ33qB9ggA/lPSy3LGmJ+RtEfSxy7y+Q8aY/YZY/aNjTGyB/C7by6Gr7sJwq6pCwd1x84OffvFEVlr3S4HACpqJUF4UNKmZR/3LT52DmPM3ZL+m6S3WGszF/pC1trPWGv3WGv3dHbSDwj43TdfGNZ1vc3a0Fzvdim+dvfubg1NpXVoaNrtUgCgolYShJ+StNMYs9UYE5H0LkkPLH+CMeYmSX+tYghmDg+AyxqdSWv/6UndfTWnwW675+puBQNGXzt4xu1SAKCiLhuErbU5SR+W9JCkFyV90Vp7yBjzUWPMWxaf9jFJjZL+zRjzrDHmgYt8OQCQJH394LCsld50XY/bpfheayyi23d06KsHhmiPAOAroZU8yVr7oKQHz3vsD5b9/u4S1wWgxv3Hc0O6qieund1xt0uBpB+/foN++/4DOjg4pev7WtwuBwAqgs1yACpucHJe+04l9RM3bHS7FCz6sd09CgeNvnqA9ggA/kEQBlBxXzswJKl4CglvaG4I646dnfragTO0RwDwDYIwgIr7j+fO6Ia+Zl3RHnO7FCzz5us2aHByXvtPT7pdCgBUBEEYQEWdGE/p4OAUbREedM813YoEA/qP54bcLgUAKoIgDKCinJD1ZtoiPKepLqwfuapLDzw7pIV8we1yAKDsCMIAKqZQsLr/6QHduq2NJRoe9fY9fZpIZfXwS4yEB1D7CMIAKubJkwn1J+b0zls2Xf7JcMWduzrV0RjV/U8PuF0KAJQdQRhAxXxx32nFoyHdew1tEV4VDgb0kzdt1HdeGtXEbMbtcgCgrAjCACpiJr2gBw+e0U/cuFH1kaDb5eASfupVm5QrWH35WS7NAahtBGEAFfHVA2eUXijoHXtoi/C6K3viur6vmfYIADWPIAygIr7w1Gnt6m7UDX3NbpeCFXj7nk168cy0nmWmMIAaRhAGUHaHhqb07OlJvWPPJhlj3C4HK/DWGzcqFgnqHx4/6XYpAFA2BGEAZfe5H5xUfTiot7+KtohqEa8L62039+mrB85waQ5AzSIIAyirZCqrrzw7pLfe1KvmhrDb5WAVfvY1VyibK+gL+067XQoAlAVBGEBZfWHfaWVyBb33tivcLgWrtKs7rlu3temfn+hXvmDdLgcASo4gDKBs8gWrf3z8lG7d1qareprcLgdr8J7XbNHg5Ly+/eKI26UAQMkRhAGUzbdeHNHg5Lx+7rYtbpeCNbpnd7c2Ntfpbx894XYpAFByBGEAZWGt1ae/d0x9rfW6++put8vBGoWDAb3vtVv15ImEnulPul0OAJQUQRhAWTxxPKH9/ZP6hddtVyjIW001e/fezWquD+vT3z3mdikAUFJ8dwJQFn/1vWPqaIzo7a/qc7sUrFMsGtJ7XnOFvvniiI6OzrpdDgCUDEEYQMk9PzilR46M6X2v3aq6cNDtclAC771tiyLBgP7mkeNulwIAJUMQBlByf/XdY4pHQ/qZWxmZVis6GqN65y2b9KX9AxpIzrldDgCUBEEYQEm9eGZaXzt4Ru+57Qo11bFAo5b80l3bZYzRX3z7qNulAEBJEIQBlNSffeOI4nUhffCO7W6XghLb0Fyvn967Wfc/M6CT4ym3ywGAdSMIAyiZ/f1JfevFEf3CndtYp1yjfvn12xUOGn3i2y+7XQoArBtBGEDJfPwbh9Uei+jnb9/qdikok654nd77mi368rODenlkxu1yAGBdCMIASuKxo+N67OiEfvn1OxSLhtwuB2X0C6/brsZISH/69ZfcLgUA1oUgDGDdcvmC/uSrL6ivtV7/9dWb3S4HZdYWi+jDP7JD33lpVN9/ecztcgBgzQjCANbt80+d1kvDM/rdN13N3GCf+Lnbt2hzW4P++1dfVC5fcLscAFgTgjCAdZmaX9Cff/OI9m5t0xuv7XG7HFRINBTU77zxKh0emdEX9p12uxwAWBOCMIB1+cS3XlZyLqs/+PHdMsa4XQ4q6N5re7R3a5s+/tBhJVJZt8sBgFUjCANYs+cHp/TZH5zQu27ZrGt7m90uBxVmjNFH77tGM+mc/vTBF90uBwBWjSAMYE1y+YJ+50sH1RaL6iP3XuV2OXDJVT1N+sCd2/RvTw/oieMTbpcDAKtCEAawJp/9wUkdHJzSH71lN8szfO5Xf2SnNrXV67/9+0Flcnm3ywGAFSMIA1i104k5/dk3juhHrurSm6/b4HY5cFl9JKg/ue9aHRtL6ZNsnANQRQjCAFYlX7D6zS8+p1DA6E/eei0X5CBJuuvKLr1jT5/+6rvH9Ex/0u1yAGBFCMIAVuVvvn9cT55M6I/vu0a9LfVulwMP+f0f360NzfX6zS8+p/ksLRIAvI8gDGDFXhia1p9947DedF2PfvKmXrfLgcfE68L62Nuv14nxlP5vpkgAqAIEYQArksrk9Cv/+oxaGiL6H2+9jpYIXNBt2zv0gTu26h+fOKWvHzzjdjkAcEkEYQCXZa3V7/77QZ0YT+kT77pRrbGI2yXBw37rx67SjZta9Nv3H9CpiZTb5QDARRGEAVzWvzzZr688O6TfuGeXbtve4XY58LhIKKD/9dM3yRjpQ//yjNIL9AsD8CaCMIBLeqY/qT9+4AXduatTv3zXDrfLQZXoa23Qn73jRj0/OK3f/dJBWWvdLgkAXoEgDOCiBifn9cF/eFo9zXX6xDtvVCBAXzBW7p7d3fr1u3fpS/sH9bffP+F2OQDwCgRhABc0l83pA5/bp8xCXn/33j30BWNNfuVHduhN1/XoT7/+oh5+adTtcgDgHARhAK+Qyxf0q/+6Xy8NT+uTP32TdnbH3S4JVSoQMPr422/QVT1N+tC/PKMDA5NulwQASwjCAM7hTIj41ouj+uP7rtXrr+xyuyRUuYZISJ/9+VvU2hDR+z77FJMkAHgGQRjAOT720GF9cd+AfvUNO/Wzt17hdjmoEV1NdfqH9+9VrmD1nr9/UiPTabdLAgCCMICzPvntl/WX3z2md+/drF+/e6fb5aDGbO9s1N//3C0an8no3Z95QqOEYQAuIwgDkCR96uGj+vNvHtHbbu7Vf3/rtWyOQ1ncvLlVn3vfXg1Pp/Wuv3lCozOEYQDuIQgDPmet1Se//bI+9tBhvfXGjfrYT92gIGPSUEZ7trQVw/BUungyTBgG4BKCMOBjhYLVHz1waOkk+ONvJwSjMm7Z0qbP/vxenZlK651//YT6J+bcLgmADxGEAZ/K5gr6tS88q889fkofuGOrPv5TNygU5C0BlbN3a5v+8f17lZzL6m1/9ZgODky5XRIAn+G7HuBD0+kFvf9zT+mB54b0kTdepd9909VsjYMrXnVFm+7/xdsUDQX1zs88ru8dGXO7JAA+QhAGfOb42Kx+8lOP6fFjE/qf/+V6/eLrtnMxDq7a0dWof//l27SlPab3f/Yp/eMTp2StdbssAD5AEAZ85LuHR3Xfpx5Tcm5B//R/vFrvuGWT2yUBkopzhr/wC7fqjp0d+v0vP6/fuv+A0gt5t8sCUOMIwoAPFApWn3r4qN732afU19qgr3zodt26rd3tsoBzxOvC+rv33qJffcNO3f/0gN7+6cc1kOQSHYDyIQgDNW50Oq33/P2T+thDh/XG6zbof//Sa7SprcHtsoALCgSMfuOeXfqb9+zRyfGUfvwvHtV/Pj/sdlkAahRBGKhh33lpRPd+4vvadyqh/+dt1+l/vfsmNURCbpcFXNY9u7v1lQ/frr7Wev3iPz2t/+v+A0plcm6XBaDGEISBGjSdXtDvfOmg3vfZfeqKR/XVX3mt3rV3M5fiUFW2dTbqS790u375ru364tOn9eZPfl9Pn0q6XRaAGkIQBmrMt14Y0T1//j194al+feCOrfryh27Xjq6422UBaxIJBfTb916lz3/gVi3krX7q0z/QH37lec2kF9wuDUAN4GekQI0Ynkrrfzz4ov7juSFd1RPXZ352j27Y1OJ2WUBJvHpbux769Tv18YcO63OPn9RDh0b0x/ddox+7psft0gBUMePWrMY9e/bYffv2ufLaQC1JL+T1d4+e0KcePqpcwepDd+3QL921XZEQP/BBbXr29KQ+8r8P6KXhGb1uV6d+781Xa2c3P/UAcHHGmKettXte8ThBGKhO1lp9/flh/enXX9TpxLx+dHe3fu/Nu7W5nYkQqH0L+YI+94OT+sS3X9ZcNq//+urN+rW7d6ktFnG7NAAeRBAGaoS1Vt9/eVwfe+iwDg5OaWdXo/7wJ67Ra3d2uF0aUHGJVFb/37eO6J9/2K+GcFDvv2Or3vfarWqqC7tdGgAPIQgDVc5aq8ePTeiT33lZTxxPqLelXr9290795E29CgVpg4C/vTwyo49/47AeOjSi5vqwPnjnNv3cbVsUi3IVBgBBGKha+YLVQ4eG9envHdOBgSl1xqP68Ot36F17NykaCrpdHuApBwem9OffPKyHD4+puT6sn7l1s9572xZ1xevcLg2AiwjCQJWZz+b15WcH9ZlHjuvEeEpb2hv0wTu3620396ouTAAGLmV/f1J//b3jeuiFYYUDAf3kTb36wJ3btKOr0e3SALiAIAxUiWNjs/rnJ/p1/9OnNZ3O6breZv3i67br3mt7FAywEANYjRPjKf3t94/r/qcHlMkVdOeuTv303s16w9VdCtNSBPgGQRjwsEwur2+9MKp//uEp/eDYhMJBox+7pkc/c+sVevXWNjbCAes0MZvRPz3Rr88/1a8zU2l1xaN6x55Neuctm7SpjUkrQK0jCAMeUyhYPd2f1JeeGdTXDgxpOp1Tb0u9fvrVm/WOPZvUGY+6XSJQc3L5gr57eEz/+mS/Hj48Kivp1q3tuu/GjXrjdRvUXM+0CaAWEYQBD7DW6vnBaf3noTP6yrNDGkjOqz4c1L3X9uitN/XqtTs6aH8AKmRocl7/tm9AX3l2UMfHU4oEA3r9VZ2678Ze3XVlpxoiTJwAagVBGHBJLl/QUyeTeujQsL75wogGJ+cVMNLtOzr0tpt79aO7exjxBLjIWquDg1P68v4h/ceBIY3NZBQNBXTHzg7ds7tbb7i6Wx2N/IQGqGYEYaCCkqmsHj06ru8dGdO3XxxRcm5BkVBAd+7s0I9e06O7r+5mAxbgQfmC1Q+PT+gbL4ws/Y+rMdKrNrfqnt3det2VnbqyO07fPlBlCMJAGWVzBe3vT+r7L4/rkZfHdHBwStZKTXUhvf6qLv3YNT163a5OTn6BKmKt1QtnpvXNxVB8aGhaktQZj+q1OzqKv3Z2qLuJGcWA1xGEgRKay+b07OlJPXUiqX2nEnrmVFKpbF7BgNFNm1p0x85O3bGrQ9f3NrP1DagRQ5PzevTouB59eVyPHR3XRCorSdrZ1ahbt7Vrz5ZW7d3apg3N9S5XCuB8BGFgHUan09p/elL7Tib05MmkDg1OKVewMka6sjuuvVvbdPuODr1me7ua6rh1DtS6QsHqxeFpPXZ0XI8endDTJxNKZfOSpN6Wet2ypVV7trTpps0t2tUdZ2Yx4DKCMLBCYzMZPT84pQMDUzo4OKmDg1Mamc5IkiKhgG7sa9EtW4vf5G7e3Mq4JQDK5Qt6aXhGT51MaN/JpJ48mdDYzNn3jd0bmnR9X7Ou623W9X0t2tHVyIQYoIIIwsB5srmCjo/P6vDwjI6MzOjw8IwODU3rzFRakmSMtK0jpuv7WnRdb7Nu2NSsa3ubFQ2x3hjApVlrdToxr+cGJnVgYFIHBqb0/ODU0qlxfTioqzfEdWVPk67sbtSunriu6mniEi1QJgRh+FZ6Ia/TiTkdG5vVkZFZHR6Z0ZHhGZ0YTylXKP73HwwYbeuI6eplpzbX9DarkcttAEqkULA6Pp7SwcFiMD40NK0jIzOanFtYek5HY1RX9jRqV3dcV/XEta2zUVs7YmqPRZhUAawDQRg1bSFf0EByXifHUzo+ntLJ8ZROLP4ampqX85+5MdLmtgbt7IovfbO5sieurR0xTnoBVJy1VmMzGR1e/KnU4eGZ4v+sj8wovVBYel68LqStHbFX/LqiPUZ7FrACFwvCHHehKmRyeZ2ZTGtwcl4DyTkNJOc1mJwv/nNyXmem5lVY9v908bqQtnXEtGdLq7Z29C1909jR1ci2KACeYYxRV1OduprqdMfOzqXHCwWr08k5HR9P6cRYSicniv9jv+9kUg88NyR73vtdX2uD+lrr1ddar96W+nM+bq4Pc5oMXASJAK6y1mo6ndPIdHrxV0Yj02mNLv5+eDqtM1PzGp3JnPPGHzDShubiG/6rt7apt7Vem9satK0zpi3tMbXxY0QAVSwQMLqivXji+/orz/1ceiGv/sScToyndGoitXQo0D8xpx8cHV/qQ3Y0RkPa0FynnuY6dTfVqbspqp4m5/fFxzsao1zegy8RhFFyTrhNpLJKpDKamM1qIpVVIpXVxGxWozNpjU5nNDJTDL/Lf/znaKoLLb1Jv3ZH59LJhnPK0dNcxzgiAL5UFw5qV3dcu7rjr/ictVZT8wsaSJ796dlAcl5Dk/Mamcno5ZFxjc1mlC+c2xYZMMVFIT1NdeqM16mjMaL2xojaY9Fz/9kYUVtDhPnoqBkEYVxSLl/QdDqnqfmFV/6aOzfgTiwG30Qqq4X8hXvPY5Fg8ceA8ahu6GtRd1NU3Ys/FuyOR9XTXKeueJ3qI/TrAsBqGWPU0hBRS0NE1/Y2X/A5+YLVxGxm6aduzk/khqfSGp4utqAdGJhUIpVdulB87mtILfVhtTdG1R4rhuPm+ohaGsJqqQ+rpSF89uOGsFoWf18X5n0d3kMQrmHWWs0v5DWbyWk2nVMqk9dMZkGpTF6pTE4zmZyml0Lt2YA7Ob+w9PhsJnfJ12iMhtS2+EbY21Kn63qb1BaLqqMxorZY8VdHY3Tp97wRAoC7goGzfcnX6cJhWSr2KU+nFzQ+m9XEbEYTqeI/x2ezmnB+2jeb/f/bu/sQy+o6juPv7znn3juzs+uutSvZuqT4EGbYYiYSBPpHZRINUeGSkD2AGdlfUWCCSRJYFlJYguWyGJJYVCyykUJQf5kuoZaLxlKau2mu2j7eufeeh29/nN+9c2acWWd2Hu7OPZ8XXM7v/M7v3PuD78ydzzn33DM8/8oxjkxlHG7PHZz7Wkk0CMYbQ2g+Y7zB+lbChrHy/3efEgAACNxJREFUsb7VYH1ob2gloV2OWd9KdPmGLDsF4dNElhdMpTlTvbxcpjntXk6nVy77fVNh/UQ3KwNuP+T2Mo51Mk50s0HIPdHNOMl70sBYI2LjeGPw2LppjIvP3jCjr/rYtK5889o43tCdFkRERlQUTZ9dvuCs9W853t1p93IOT6Ucbvc40k5DO+XwVFgP7cPtlBdfb3O8m3G0U550WchNrCaa8YxwvGEsYaKZsK4ZM96MwzJhotJeF9rrKuMmmslgvC6zqzcF4TlkeUE36z9yev12Wq6fvH96v25aDIJtu5fTSWe2272MqV5OJy3o5W++TvZkIoOJVnnEPDE4ai6/EDERjpzXh/4Z42b0x5wxpo+rRERk6cyMifC3Zuum8UXt2w/RxzoZx7spxzpZaJcne/ph+Xil/1g341gn5dWjXdppRrs7feJoMRqxMd6ImWiV4biVxIw1IsbCcrDeiGklYdmYtS2JaQ32mWO/sN5KIppxRKQz26eNWgXhff85yh2P7KOXTwfVQWjNikGwnf0lglPRSiKaSVQefTbKo9LxRsR4M+bMdc3ySLRRHpmONeLKuMqy3w7bB+NCv+6KICIio6AaomFsSc9VFD446TTVy2mn5aWB5SeqGVNpzonu9MmodprT7ma0Q7ublpmgk+a8djyjm5UnrDppeRKrf8JrKeLIaMYRjdhoJjHN2GiG3NCIp5et/noc0QghuplEg/Fzjk2mxzciI4kjkthIIiOJytdM4qhcj+foC/s0wrYkspEO7gsKwmZ2DfAjIAZ+7u53ztreAh4A3g+8Dlzn7i8s71SXzgyyomC8EbNpvEGrUf6wtJL+0VrZbib9dkQrHAHO3x+F/unnaMaRQqqIiMgQRFE1VK+MovDypFpa0MnyEJKLGaG5H6Y7aU4nK+imOb28IM2cXl5+qpzmTjcrSPMirJfLXlge72aVbT5jW3/8ya7LXi6RUYbjyIgjozEI19Mhe0ZfCNmNEK7jSt89n71sxee7GG/5U2JmMfAT4MPAAeBJM9vt7vsqw74E/M/dLzCzHcD3gOtWYsJLcfHZZ/Crmz447GmIiIjIGhZFxlhUflK7keH+Z79+KC9DdjEjbKe5k+VOWpSfdqd5QZY7WdFfzuorwvgQsKv7pEVBPmuf/vNW9ymfo9x+PMvI8vJ5sqJY0HXgq20hh0tXAPvd/Z8AZvYQMAlUg/AkcHto/xq4x8zMh/X/m0VERERqoBrKZfEW8lXJrcBLlfUDoW/OMe6eAUeAty/HBEVEREREVsKq3jPEzG40s71mtvfQoUOr+dIiIiIiIjMsJAgfBLZV1s8JfXOOMbME2Ej5pbkZ3P0+d7/c3S/fsmXLqc1YRERERGQZLCQIPwlcaGbnmVkT2AHsnjVmN3BDaH8a+KOuDxYRERGR09lbflnO3TMzuxn4A+Xt03a6+7Nm9h1gr7vvBu4HfmFm+4E3KMOyiIiIiMhpa0E32XP3PcCeWX23Vdod4DPLOzURERERkZWjf7AtIiIiIrWkICwiIiIitaQgLCIiIiK1pCAsIiIiIrWkICwiIiIitaQgLCIiIiK1pCAsIiIiIrWkICwiIiIitaQgLCIiIiK1pCAsIiIiIrWkICwiIiIitaQgLCIiIiK1pCAsIiIiIrWkICwiIiIitaQgLCIiIiK1pCAsIiIiIrWkICwiIiIitWTuPpwXNjsEvDiUF6+nzcBrw56ErCjVuB5U53pQnetBdV4973L3LbM7hxaEZXWZ2V53v3zY85CVoxrXg+pcD6pzPajOw6dLI0RERESklhSERURERKSWFITr475hT0BWnGpcD6pzPajO9aA6D5muERYRERGRWtIZYRERERGpJQXhEWBmO83sVTP7e6XvLjN7zsyeMbPfmtmmyrZbzGy/mT1vZh8dzqxlseap8x2hxk+Z2aNm9s7Qb2b241DnZ8zssuHNXBZjrjpXtn3dzNzMNod11XmNmuf3+XYzOxh+n58ys2sr2/S+vQbN9/tsZl8Lf6OfNbPvV/pV51WmIDwadgHXzOp7DHivu18K/AO4BcDM3gPsAC4J+/zUzOLVm6oswS7eXOe73P1Sd98OPALcFvo/BlwYHjcC967WJGXJdvHmOmNm24CPAP+udKvOa9cu5qgzcLe7bw+PPaD37TVuF7PqbGZXA5PA+9z9EuAHoV91HgIF4RHg7n8G3pjV96i7Z2H1ceCc0J4EHnL3rrv/C9gPXLFqk5VTNk+dj1ZWJ4D+Rf+TwANeehzYZGZnr85MZSnmqnNwN/BNpmsMqvOadZI6z0Xv22vUPHX+CnCnu3fDmFdDv+o8BArC9fBF4PehvRV4qbLtQOiTNcrMvmtmLwHXM31GWHUeIWY2CRx096dnbVKdR8/N4TKXnWZ2ZuhTnUfLRcCHzOwvZvYnM/tA6Fedh0BBeMSZ2a1ABjw47LnIynD3W919G2WNbx72fGR5mdk64FtMH+TI6LoXOB/YDrwM/HC405EVkgBvA64EvgE8bGY23CnVl4LwCDOzzwMfB6736fvkHQS2VYadE/pk7XsQ+FRoq86j43zgPOBpM3uBspZ/NbN3oDqPFHf/r7vn7l4AP2P6Y3HVebQcAH4TLml6AiiAzajOQ6EgPKLM7BrK6wk/4e7tyqbdwA4za5nZeZRfsnliGHOUpTOzCyurk8Bzob0b+Fy4q8CVwBF3f3nVJyhL5u5/c/ez3P1cdz+X8o/oZe7+CqrzSJl1ffcngf6dBvS+PVp+B1wNYGYXAU3gNVTnoUiGPQFZOjP7JXAVsNnMDgDfprxLRAt4LHzi8ri73+Tuz5rZw8A+yksmvuru+XBmLosxT52vNbN3U55ReBG4KQzfA1xL+WWLNvCFVZ+wnJK56uzu988zXHVeo+b5fb7KzLZTfiHyBeDLAHrfXrvmqfNOYGe4pVoPuCF8aqs6D4H+s5yIiIiI1JIujRARERGRWlIQFhEREZFaUhAWERERkVpSEBYRERGRWlIQFhEREZFaUhAWERERkVpSEBYRERGRWlIQFhEREZFa+j8xfg68U1RKlwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(f_list, fit_result)" ] }, { "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.6.9" } }, "nbformat": 4, "nbformat_minor": 2 }