{ "cells": [ { "cell_type": "code", "execution_count": 106, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from pathlib import Path\n", "import pandas as pd\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import scipy\n", "from glob import glob" ] }, { "cell_type": "markdown", "source": [ "Prior code from mva:" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 112, "outputs": [], "source": [ "def read_xy(path: str):\n", " df = pd.read_csv(path, skipinitialspace=True)\n", " df['N'] = df.index + 1\n", " df['r'] = (df.x ** 2 + df.y ** 2) ** 0.5\n", " df['cr'] = df.r.cummax()\n", " df['fd'] = np.log(df.N) / np.log(df.cr)\n", "\n", " return df\n", "\n", "\n", "def read_all(load_dir: str):\n", " paths = glob(f'{load_dir}/*.csv')\n", " return [read_xy(path) for path in paths]\n", "\n", "\n", "# We expect the format of sp data to be:\n", "# ROOT/{probability}/{RUN_ID}.csv\n", "\n", "def read_sp_dir(load_dir: str):\n", " probability = float(Path(load_dir).name)\n", " return probability, read_all(load_dir)\n", "\n", "\n", "def read_sp_full(root: str):\n", " if not Path(root).exists():\n", " raise Exception(\"Root does not exist\")\n", "\n", " raw_data = [read_sp_dir(load_dir) for load_dir in glob(f'{root}/*')]\n", " fd_data = [fd_of_dfs(dfs) for (p, dfs) in raw_data]\n", " probabilities = [p for (p, dfs) in raw_data]\n", " fd_means, fd_stds = list(zip(*fd_data))\n", " out = pd.DataFrame(zip(probabilities, fd_means, fd_stds), columns=['p', 'fd_mean', 'fd_std'])\n", " out.sort_values(by='p', inplace=True)\n", "\n", " return out" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": null, "outputs": [], "source": [ "def convergent_tail_index(series, tol):\n", " diffs = np.abs(np.ediff1d(series))\n", " for i in range(0, len(diffs)):\n", " if np.max(diffs[i:]) <= tol:\n", " return i\n", "\n", " # No convergence found\n", " return None\n", "\n", "\n", "def mean_of_tail(series, tol=0.05):\n", " tail_index = convergent_tail_index(series, tol)\n", " if tail_index is None:\n", " print(\"Ahh fuck\")\n", "\n", " return np.mean(series[tail_index:])\n", "\n", "\n", "def fd_of_dfs(dfs):\n", " fds = [mean_of_tail(df.fd, 0.001) for df in dfs]\n", " fds_clean = [f for f in fds if f < np.inf]\n", " return np.mean(fds_clean), np.std(fds_clean)" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 113, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/joshuacoles/Library/Caches/pypoetry/virtualenvs/data-analysis-B4Au_hWl-py3.10/lib/python3.10/site-packages/numpy/lib/arraysetops.py:89: RuntimeWarning: invalid value encountered in subtract\n", " return ary[1:] - ary[:-1]\n" ] } ], "source": [ "c_sticking_probability = read_sp_full(\"../data/stick-probability\")" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 121, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/joshuacoles/Library/Caches/pypoetry/virtualenvs/data-analysis-B4Au_hWl-py3.10/lib/python3.10/site-packages/numpy/lib/arraysetops.py:89: RuntimeWarning: invalid value encountered in subtract\n", " return ary[1:] - ary[:-1]\n" ] } ], "source": [ "rust_sticking_probability = read_sp_full(\"../data/rust-sticking-probability\")" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 124, "outputs": [ { "data": { "text/plain": "" }, "execution_count": 124, "metadata": {}, "output_type": "execute_result" }, { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGdCAYAAAAxCSikAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABcNElEQVR4nO3dd3gU5d7G8e/spodUSIWE3kFAFAgIiKKICkY9iPrSFLGBDSvqUazYj0flgIqIBURBQAVEEQSkSwkSegkESKEnpG+Z94/VaAQkgSS7CffnuvY67uwzs7+Zg+ztzFMM0zRNRERERDyYxd0FiIiIiJyJAouIiIh4PAUWERER8XgKLCIiIuLxFFhERETE4ymwiIiIiMdTYBERERGPp8AiIiIiHs/L3QWUF6fTSVpaGkFBQRiG4e5yREREpBRM0+TEiRPExsZisZz+Pkq1CSxpaWnExcW5uwwRERE5C/v27aNOnTqn/bzaBJagoCDAdcLBwcFurkZERERKIzs7m7i4uOLf8dOpNoHlj8dAwcHBCiwiIiJVzJm6c6jTrYiIiHg8BRYRERHxeAosIiIi4vGqTR8WERGRimSaJna7HYfD4e5SqhSr1YqXl9c5TzmiwCIiInIGRUVFpKenk5eX5+5SqqSAgABiYmLw8fE562MosIiIiPwDp9NJSkoKVquV2NhYfHx8NEFpKZmmSVFREYcOHSIlJYXGjRv/4+Rw/0SBRURE5B8UFRXhdDqJi4sjICDA3eVUOf7+/nh7e7N3716Kiorw8/M7q+Oo062IiEgpnO2dASmfa6erLyIiIh5PgUVEREQ8ngKLiIiIeLwyBZYxY8Zw8cUXExQURGRkJImJiWzbtu0f99m0aRM33ngj9erVwzAM3n777VO2Gzt2LPXq1cPPz4+OHTuyevXqspQmIiIip5CRkcF9991HgwYN8PX1JS4ujj59+rBgwQJ3l1YmZQosixcvZvjw4axcuZL58+djs9m48soryc3NPe0+eXl5NGjQgFdeeYXo6OhTtvnyyy8ZOXIkzz77LOvWraNNmzb06tWLgwcPlu1syptpwuLXYeY9UJDt3lpERETKaM+ePbRv356FCxfy+uuvs3HjRubNm0ePHj0YPny4u8srE8M0TfNsdz506BCRkZEsXryYbt26nbF9vXr1ePDBB3nwwQdLbO/YsSMXX3wx7733HkDx8LH77ruPJ554olS1ZGdnExISQlZWVvmu1vx6I8g9BHcuhti25XdcERGpEgoKCkhJSaF+/frFQ3JN0yTf5p4Zb/29raWeB+bqq6/mt99+Y9u2bQQGBpb47Pjx44SGhlZAhSc71TX8Q2l/v89pHpasrCwAwsPDz/oYRUVFrF27llGjRhVvs1gs9OzZkxUrVpxLeeXCGVYfS+4hHEd2Y1VgERERIN/moMUzP7jluzc/34sAnzP/fB89epR58+bx0ksvnRRWgEoLK+XlrDvdOp1OHnzwQbp06UKrVq3OuoDDhw/jcDiIiooqsT0qKoqMjIzT7ldYWEh2dnaJV3kzTZPvUn0ByE7bXu7HFxERqSg7d+7ENE2aNWvm7lLKxVnfYRk+fDjJycksXbq0POsptTFjxvDcc89V6HcYhsEx/zgohMKDuyr0u0REpOrw97ay+flebvvu0jiHHh8e6awCy4gRI5g9ezZLliyhTp0651RArVq1sFqtZGZmltiemZl52k66AKNGjWLkyJHF77Ozs4mLizunWk6lsEY8FIJxbHe5H1tERKomwzBK9VjGnRo3boxhGGzdutXdpZSLMj0SMk2TESNGMHPmTBYuXEj9+vXPuQAfHx/at29fYniV0+lkwYIFJCQknHY/X19fgoODS7wqglHTdY7+OakVcnwREZGKEB4eTq9evRg7duwpR/MeP3688os6B2UKLMOHD+fzzz9nypQpBAUFkZGRQUZGBvn5+cVtBg0aVKIDbVFREUlJSSQlJVFUVMSBAwdISkpi586dxW1GjhzJhx9+yCeffMKWLVu45557yM3N5bbbbiuHUzw3fpGNAQguOgi2/DO0FhER8Rxjx47F4XDQoUMHvv76a3bs2MGWLVt45513/vGmgCcq0/2scePGAXDppZeW2P7xxx8zZMgQAFJTU0sscpSWlka7du2K37/xxhu88cYbdO/enUWLFgHQv39/Dh06xDPPPENGRgZt27Zl3rx5J3XEdYfIqBiyTX+CjXw4thciq0fnJRERqf4aNGjAunXreOmll3j44YdJT08nIiKC9u3bF/+mVxXnNA+LJ6moeVg2pWXhHN+N1pY9cMtUaNq73I4tIiKe75/mEJHSKY95WLSW0BnEhQew13Td6Sk8uPMMrUVERKQiKLCcQbCfN5nWGADyMna4uRoREZHzkwJLKeQExgPgOKKhzSIiIu6gwFIKztB6AHhn7XFrHSIiIucrBZZS8I5oBECN/DRw2N1cjYiIyPlHgaUUQqLiKTS9seKArH3uLkdEROS8o8BSCnHhgew1I11vjqW4txgREZHzkAJLKfx1aLN5VIFFRESksimwlELtUH9Sfw8sBZka2iwiIlLZFFhKwc/bylFf16rUhQd3ubkaERGR0hkyZAiGYWAYBt7e3tSvX5/HHnuMgoKCcjl+vXr1ePvtt8vlWGfi2Wtje5Ci4Hg4BpbjeiQkIiJVx1VXXcXHH3+MzWZj7dq1DB48GMMwePXVV91dWpnoDkspGeENAPDP2QfVY/klERE5D/j6+hIdHU1cXByJiYn07NmT+fPnA6e+Q9K2bVtGjx4NgGmajB49mvj4eHx9fYmNjeX+++8HXAsh7927l4ceeqj4Lk5F0h2WUgqMbIB9pwVvZwHkZEJQtLtLEhERdzFNsOW557u9A+Asw0FycjLLly+nbt26pWr/9ddf85///IepU6fSsmVLMjIy2LBhAwAzZsygTZs23HnnnQwbNuys6ikLBZZSql0rmDSzJvHGITi6W4FFROR8ZsuDl2Pd891PpoFPYKmbz549mxo1amC32yksLMRisfDee++Vat/U1FSio6Pp2bMn3t7exMfH06FDBwDCw8OxWq0EBQURHV3xv4l6JFRKcWH+7DF//z9EQ5tFRKSK6NGjB0lJSaxatYrBgwdz2223ceONN5Zq3379+pGfn0+DBg0YNmwYM2fOxG53z4zvusNSSnHhASz8ffI459HdSnoiIucz7wDXnQ53fXcZBAYG0qiRa4mZiRMn0qZNGz766COGDh2KxWLB/Fu/TJvNVvzPcXFxbNu2jZ9++on58+dz77338vrrr7N48WK8vb3P/VzKQIGllKKC/dhnuO6wFGTupGx/XEREpFoxjDI9lvEUFouFJ598kpEjR3LrrbcSERFBenp68efZ2dmkpJR8iuDv70+fPn3o06cPw4cPp1mzZmzcuJELL7wQHx8fHA5H5dReKd9SDVgtBrmB8QA4jux2czUiIiJnp1+/flitVsaOHctll13GZ599xi+//MLGjRsZPHgwVqu1uO2kSZP46KOPSE5OZvfu3Xz++ef4+/sXd9qtV68eS5Ys4cCBAxw+fLhC69YdljJwhNSDAvDJ3uPuUkRERM6Kl5cXI0aM4LXXXmPHjh2kpKRw7bXXEhISwgsvvFDiDktoaCivvPIKI0eOxOFw0Lp1a7777jtq1qwJwPPPP89dd91Fw4YNKSwsPOnxUnkyzIo8eiXKzs4mJCSErKwsgoODK+Q7np2+mueSr3C9eXwP+IdVyPeIiIjnKCgoICUlhfr16+Pn5+fucqqkf7qGpf391iOhMoiuVZNMM9T1RiOFREREKo0CSxnEhwf8ZWiz+rGIiIhUFgWWMogL9yfV6RrazDHdYREREaksCixlEBcWwF4zCgD7Ya3aLCIiUlkUWMogNMCbTC/XVMy2QwosIiIilUWBpQwMwyA/uD4A1qMKLCIi55NqMqjWLcrj2imwlJFRsyEAPoVHIP+4e4sREZEK98cU9Hl5blqduRr449qdy3T+mjiujCJr1SIjJYxo4xgc3QW127u7JBERqUBWq5XQ0FAOHjwIQEBAAIZhuLmqqsE0TfLy8jh48CChoaElZtEtKwWWMooL8yfFGUO09Rgc3qnAIiJyHoiOdk1p8UdokbIJDQ0tvoZnS4GljOJrBpBiRpPAZjiy093liIhIJTAMg5iYGCIjI0usZixn5u3tfU53Vv6gwFJG8eEBLDdjADCP7EQ3BUVEzh9Wq7Vcfnyl7NTptozqhAWw+/fA4ji0w83ViIiInB8UWMrIz9vKicB6ABhHd4GGuYmIiFQ4BZaz4BVeD5tpxWrPgxPp7i5HRESk2lNgOQuxNYNJNX9fU0gdb0VERCpcmQLLmDFjuPjiiwkKCiIyMpLExES2bdt2xv2mTZtGs2bN8PPzo3Xr1sydO7fE50OGDMEwjBKvq666qmxnUonq/j5SCIDD6sciIiJS0coUWBYvXszw4cNZuXIl8+fPx2azceWVV5Kbm3vafZYvX84tt9zC0KFDWb9+PYmJiSQmJpKcnFyi3VVXXUV6enrx64svvji7M6oE8eEBpPze8ZYjmqJfRESkopVpWPO8efNKvJ80aRKRkZGsXbuWbt26nXKf//73v1x11VU8+uijALzwwgvMnz+f9957j/Hjxxe38/X1PedJZSpLXHgAq4sDix4JiYiIVLRz6sOSlZUFQHh4+GnbrFixgp49e5bY1qtXL1asWFFi26JFi4iMjKRp06bcc889HDly5B+/u7CwkOzs7BKvyhIf/ufQZqcCi4iISIU768DidDp58MEH6dKlC61atTptu4yMDKKiokpsi4qKIiMjo/j9VVddxaeffsqCBQt49dVXWbx4Mb1798bhcJz2uGPGjCEkJKT4FRcXd7anUma1aviQ7lUbAOPYHrAXVdp3i4iInI/Oeqbb4cOHk5yczNKlS8+5iJtvvrn4n1u3bs0FF1xAw4YNWbRoEZdffvkp9xk1ahQjR44sfp+dnV1pocUwDPzDapN73JdACuH4XqjVuFK+W0RE5Hx0VndYRowYwezZs/n555+pU6fOP7aNjo4mMzOzxLbMzMx/7K/SoEEDatWqxc6dp3/c4uvrS3BwcIlXZYqrGfiXjrd6LCQiIlKRyhRYTNNkxIgRzJw5k4ULF1K/fv0z7pOQkMCCBQtKbJs/fz4JCQmn3Wf//v0cOXKEmJiYspRXqf7aj0VDm0VERCpWmQLL8OHD+fzzz5kyZQpBQUFkZGSQkZFBfn5+cZtBgwYxatSo4vcPPPAA8+bN480332Tr1q2MHj2aNWvWMGLECABycnJ49NFHWblyJXv27GHBggVcd911NGrUiF69epXTaZY/11wsusMiIiJSGcoUWMaNG0dWVhaXXnopMTExxa8vv/yyuE1qairp6X9OV9+5c2emTJnCBx98QJs2bZg+fTqzZs0q7qhrtVr57bff6Nu3L02aNGHo0KG0b9+eX375BV9f33I6zfIXFx7Abufvj7U0F4uIiEiFKlOnW7MUC/0tWrTopG39+vWjX79+p2zv7+/PDz/8UJYyPILrkVAsAOaRnRhurkdERKQ601pCZ6lOmD97cd1hMXIyoKDy5oERERE53yiwnCVfLys1gsM5ZIa4NhzVYyEREZGKosByDuL+OlJI/VhEREQqjALLOYgPD2C3UyOFREREKpoCyzlwDW3+faSQ5mIRERGpMAos5yAuXHOxiIiIVAYFlnMQ//c+LKUY9i0iIiJlp8ByDuLDA0g1o3CYBhSdgJzMM+8kIiIiZabAcg7CA33w9fVjnxnp2qB+LCIiIhVCgeUcGIZBXHgA283fV6zOTHZvQSIiItWUAss5ig/3J9n5+6rV6RvcW4yIiEg1pcByjuLDA0g267nepCW5sxQREZFqS4HlHMXXDGTjH3dYDm+Dojz3FiQiIlINKbCco/jwAA4RxhEjDEyn+rGIiIhUAAWWcxQfHgDARmc91wb1YxERESl3CiznqHaoPxYDNjjquTaoH4uIiEi5U2A5Rz5eFmJCNFJIRESkIimwlIPmMcF/BpZDW8BW4N6CREREqhkFlnJwVato0gnnuBEMTjsc3OTukkRERKoVBZZycEXzKLytFjbY67k26LGQiIhIuVJgKQchAd50aVRLE8iJiIhUEAWWcnJ1qxg2Ohu43ugOi4iISLlSYCknV7SIYguujrdm5mawF7m5IhERkepDgaWchAX6EN+gGcfNQAxnkWu0kIiIiJQLBZZydPUFsSRrxlsREZFyp8BSjq5sEcXm3x8LnUhZ4+ZqREREqg8FlnJUs4YvhRGtAcjbu9bN1YiIiFQfCizlLL5lZwBCs7eDw+7makRERKoHBZZy1vnii8k2/fGliIzd6sciIiJSHhRYyllEsD/7fBsBsG3dL26uRkREpHpQYKkIMW0ByNmjfiwiIiLlQYGlAsQ07wRAdO5WcgrVj0VERORcKbBUgPBGrsDSykhh/e4MN1cjIiJS9SmwVISaDTlhDcPXsJG6abm7qxEREanyyhRYxowZw8UXX0xQUBCRkZEkJiaybdu2M+43bdo0mjVrhp+fH61bt2bu3LklPjdNk2eeeYaYmBj8/f3p2bMnO3bsKNuZeBLDICuiPQDOPSvcXIyIiEjVV6bAsnjxYoYPH87KlSuZP38+NpuNK6+8ktzc3NPus3z5cm655RaGDh3K+vXrSUxMJDExkeTk5OI2r732Gu+88w7jx49n1apVBAYG0qtXLwoKCs7+zNwsoNElAMRmJ5Ff5HBzNSIiIlWbYZqmebY7Hzp0iMjISBYvXky3bt1O2aZ///7k5uYye/bs4m2dOnWibdu2jB8/HtM0iY2N5eGHH+aRRx4BICsri6ioKCZNmsTNN99cqlqys7MJCQkhKyuL4ODgsz2lcmPuX4sx4TKOmTXYMmA9nRtHurskERERj1Pa3+9z6sOSlZUFQHh4+GnbrFixgp49e5bY1qtXL1ascD0qSUlJISMjo0SbkJAQOnbsWNzmVAoLC8nOzi7x8iRGzAUUGn6EGTns3KzhzSIiIufirAOL0+nkwQcfpEuXLrRq1eq07TIyMoiKiiqxLSoqioyMjOLP/9h2ujanMmbMGEJCQopfcXFxZ3sqFcPqzbHwCwAo3L3MzcWIiIhUbWcdWIYPH05ycjJTp04tz3pKbdSoUWRlZRW/9u3b55Y6/olP/S4ARB5bT6Fd/VhERETO1lkFlhEjRjB79mx+/vln6tSp849to6OjyczMLLEtMzOT6Ojo4s//2Ha6Nqfi6+tLcHBwiZenCWveHYAL2caGfVlurkZERKTqKlNgMU2TESNGMHPmTBYuXEj9+vXPuE9CQgILFiwosW3+/PkkJCQAUL9+faKjo0u0yc7OZtWqVcVtqiqjzkU4sBBnOUTyls3uLkdERKTKKlNgGT58OJ9//jlTpkwhKCiIjIwMMjIyyM/PL24zaNAgRo0aVfz+gQceYN68ebz55pts3bqV0aNHs2bNGkaMGAGAYRg8+OCDvPjii3z77bds3LiRQYMGERsbS2JiYvmcpbv4BnEsuBkAebvUj0VERORseZWl8bhx4wC49NJLS2z/+OOPGTJkCACpqalYLH/moM6dOzNlyhSefvppnnzySRo3bsysWbNKdNR97LHHyM3N5c477+T48eNccsklzJs3Dz8/v7M8Lc9hqZsAGzcTdngtNocTb6smFxYRESmrc5qHxZN42jwsf3Amz8IyfTCbnXUpuGMxF8aHubskERERj1Ep87DImVnquhZCbGaksn77XjdXIyIiUjUpsFS0oGiy/eOwGCbHt6sfi4iIyNlQYKkEjjjXXZYamWtwOKvFEzgREZFKpcBSCUKadgWgjbmFzWmetYSAiIhIVaDAUgksdTsD0NbYydJtB9xcjYiISNWjwFIZajaiwCcMP8PGrg3qxyIiIlJWCiyVwTAwfr/LUv/oEvYdzXNzQSIiIlWLAksl8W3bD4AbrEv5YaMeC4mIiJSFAktlaXo1hd4hxBhHSVs3193ViIiIVCkKLJXFyxdHK9ddlguPziE9K/8MO4iIiMgfFFgqUUDHIQBcaVnDz+u2uLcYERGRKkSBpTJFt+ZwUHN8DAeF6790dzUiIiJVhgJLJfNqPwCATsfncjBbj4VERERKQ4GlkoV2uBUbXjS3pLJmxSJ3lyMiIlIlKLBUtoBw9kZeBoD1tyluLkZERKRqUGBxg6CEIQB0zFnA0SytLSQiInImCixuENXmKg4aEYQauWxdNNXd5YiIiHg8BRZ3sFjZU+c6AIK3fOHmYkRERDyfAoubRHa7DYDm+es5cSzTzdWIiIh4NgUWN6nXuBU7LfWxGiY7ls5wdzkiIiIeTYHFjTKjLwXA2Pa9ewsRERHxcAosbhR+YSIAjU+sprAgz73FiIiIeDAFFjdq2q4rhwmlhpHPtlXz3F2OiIiIx1JgcSOL1cru8K4A5G6c7eZqREREPJcCi5v5tbwGgLqHl+B0ON1cjYiIiGdSYHGzJgnXUmB6E8shtievcnc5IiIiHkmBxc38AoLYXuNiAA6tmeXeYkRERDyUAosHcDbuBUCtAwvdXImIiIhnUmDxAA263AhAc+d2UvemuLkaERERz6PA4gGCI+LY5d0UgJQVmvVWRETk7xRYPER23csBCEj50c2ViIiIeB4FFg9Ru6PrsVCrgnUcOnrcvcWIiIh4GAUWDxHZqD2HLBH4G0UkL/3W3eWIiIh4FAUWT2EYHIy5DICCpOnkFzncXJCIiIjnKHNgWbJkCX369CE2NhbDMJg1a9YZ9xk7dizNmzfH39+fpk2b8umnn5b4fNKkSRiGUeLl5+dX1tKqvEZX3A7ApY4VfPLzb26uRkRExHN4lXWH3Nxc2rRpw+23384NN9xwxvbjxo1j1KhRfPjhh1x88cWsXr2aYcOGERYWRp8+fYrbBQcHs23btuL3hmGUtbQqz7duR07UaEBQzm4OLJtCZkJzooLPv+AmIiLyd2UOLL1796Z3796lbv/ZZ59x11130b9/fwAaNGjAr7/+yquvvloisBiGQXR0dFnLqV4MgxqdBsNPz9KXRbz+wzbe6NfG3VWJiIi4XYX3YSksLDzp8Y6/vz+rV6/GZrMVb8vJyaFu3brExcVx3XXXsWnTpjMeNzs7u8SrOjAu6I9pWLjYsp11639l4/4sd5ckIiLidhUeWHr16sWECRNYu3YtpmmyZs0aJkyYgM1m4/DhwwA0bdqUiRMn8s033/D555/jdDrp3Lkz+/fvP+1xx4wZQ0hISPErLi6uok+lcgTHYDTqCcCNlsW8MHszpmm6uSgRERH3qvDA8u9//5vevXvTqVMnvL29ue666xg8eLDryy2ur09ISGDQoEG0bduW7t27M2PGDCIiInj//fdPe9xRo0aRlZVV/Nq3b19Fn0rlaft/ANxoXcqaPYeZl5zh5oJERETcq8IDi7+/PxMnTiQvL489e/aQmppKvXr1CAoKIiIi4pT7eHt7065dO3bu3Hna4/r6+hIcHFziVW007Q3+YUQbR7nEspEx32/F7nC6uyoRERG3qbR5WLy9valTpw5Wq5WpU6dy7bXXFt9h+TuHw8HGjRuJiYmprPI8i5cvtO4HwP/5/ELq0TxW7j7q5qJERETcp8yjhHJyckrc+UhJSSEpKYnw8HDi4+MZNWoUBw4cKJ5rZfv27axevZqOHTty7Ngx3nrrLZKTk/nkk0+Kj/H888/TqVMnGjVqxPHjx3n99dfZu3cvd9xxRzmcYhXV9v9g9QdcZqwhmBzmbEzjksa13F2ViIiIW5T5DsuaNWto164d7dq1A2DkyJG0a9eOZ555BoD09HRSU1OL2zscDt58803atGnDFVdcQUFBAcuXL6devXrFbY4dO8awYcNo3rw5V199NdnZ2SxfvpwWLVqc4+lVYTFtILIl3mYRfa0rmJecgU2PhURE5DxlmNVkCEp2djYhISFkZWVVn/4sK8bCD0+yiUZcU/A8n9zege5NTt3vR0REpCoq7e+31hLyZK1vAosXLdlJc2Mvc35Lc3dFIiIibqHA4slqREBz12zAA6w/MS85gyK7HguJiMj5R4HF013s6nh8vdcyzIJslu087OaCREREKp8Ci6er2wUimhFAAddbf+E7PRYSEZHzkAKLpzMMuGgo4HosNH9TBoV2h5uLEhERqVwKLFVBm/6Y3oE0sRygRVEyS7brsZCIiJxfFFiqAr8QjAtcM98O9Jqv0UIiInLeUWCpKn5/LNTL8ivrN2+lwKbHQiIicv5QYKkqYi7AjOuIt+Ggj2MBi7YdcndFIiIilUaBpQoxfh/ifKvXAuZu2OfmakRERCqPAktV0uI67H7hxBpHMbfNI79Ij4VEROT8oMBSlXj5Ym0/CIB/mT/y87aDbi5IRESkciiwVDFG+8EAdLVsZMXa9W6uRkREpHIosFQ14Q3Iie2MxTCJ3j2NvCK7uysSERGpcAosVVBggmuI8w3GIn7erDlZRESk+lNgqYKM5n3I8wohxjjKnpXfurscERGRCqfAUhV5+ZLX3DXzbbO0GeQU6rGQiIhUbwosVVTNrsMA6G6sZ9m6jW6uRkREpGIpsFRRRmQz9ge1wctwkr/6U3eXIyIiUqEUWKowy+9DnNsfnc2J/EI3VyMiIlJxFFiqsJjON5NDIHHGQTYs+cbd5YiIiFQYBZYqzPAJZGd0bwB8N3xGdoGNvUdyWZd6jCXbD3GiwObmCkVERMqHl7sLkHMT1nUYTJtOm9xldBv9BRnULP6sd6toxg1o78bqREREyofusFRxdVt2YrN3K3wMBy94TyLAx0LtUH8AFmw9qCHPIiJSLSiwVANNbv8A0+LNFda1bO5fyNLHe1CvZgBFdidLth9yd3kiIiLnTIGlGvCKaYlxyUOuN98/hlFwnCtaRAEwf3OmGysTEREpHwos1UW3R6BWE8jJhPnPcGXLaAAWbj2IzeF0c3EiIiLnRoGluvDyhT7vuP553adc6EymZqAPWfk2fk056t7aREREzpECS3VSNwEuuh0A6+wHubJJCAA/6rGQiIhUcQos1U3P0RAUA0d3MYwZgKsfi2ma7q1LRETkHCiwVDd+IdD7NQDq7/iYeO/jHDiez+b0bDcXJiIicvYUWKqj5n0grhOGvYAXQucAGi0kIiJVmwJLdWQYrkdDQNeceTQw0vhxkwKLiIhUXQos1VXdBGhyFRbTwcNe09icns3+Y3nurkpEROSslDmwLFmyhD59+hAbG4thGMyaNeuM+4wdO5bmzZvj7+9P06ZN+fTTT09qM23aNJo1a4afnx+tW7dm7ty5ZS1N/u7yZwCDa6yraG3s5ic9FhIRkSqqzIElNzeXNm3aMHbs2FK1HzduHKNGjWL06NFs2rSJ5557juHDh/Pdd98Vt1m+fDm33HILQ4cOZf369SQmJpKYmEhycnJZy5O/imoJF/QH4DGvqRreLCIiVZZhnsN4V8MwmDlzJomJiadt07lzZ7p06cLrr79evO3hhx9m1apVLF26FID+/fuTm5vL7Nmzi9t06tSJtm3bMn78+FLVkp2dTUhICFlZWQQHB5/dCVVHx/Zivtsew2ljoO1J3ntqJCEB3u6uSkREBCj973eF92EpLCzEz8+vxDZ/f39Wr16NzWYDYMWKFfTs2bNEm169erFixYp/PG52dnaJl5xCWF2Mi4cC8Kj1C2au2+fmgkRERMquwgNLr169mDBhAmvXrsU0TdasWcOECROw2WwcPnwYgIyMDKKiokrsFxUVRUZGxmmPO2bMGEJCQopfcXFxFXoeVVrXR7BZA7jAksKOJV9g19pCIiJSxVR4YPn3v/9N79696dSpE97e3lx33XUMHjzY9eWWs//6UaNGkZWVVfzat093Dk6rRgR0ugeAAQVT+X5jmpsLEhERKZsKDyz+/v5MnDiRvLw89uzZQ2pqKvXq1SMoKIiIiAgAoqOjycws2SE0MzOT6Ojo0x7X19eX4ODgEi85Pe8uIyi0BtLckspvP03WVP0iIlKlVNo8LN7e3tSpUwer1crUqVO59tpri++wJCQksGDBghLt58+fT0JCQmWVV/0FhOO8+E4Ars+ezIqdh9xckIiISOmVObDk5OSQlJREUlISACkpKSQlJZGamgq4HtUMGjSouP327dv5/PPP2bFjB6tXr+bmm28mOTmZl19+ubjNAw88wLx583jzzTfZunUro0ePZs2aNYwYMeIcT0/+yr/b/RRYAmhh2cuvP0x2dzkiIiKlVubAsmbNGtq1a0e7du0AGDlyJO3ateOZZ54BID09vTi8ADgcDt58803atGnDFVdcQUFBAcuXL6devXrFbTp37syUKVP44IMPaNOmDdOnT2fWrFm0atXqHE9PSggIp/DCYQD0PPgxW9Ky3FyQiIhI6ZzTPCyeRPOwlFLeUQpeb4Gfmc9HdV5i6B26iyUiIu7jMfOwiIcJCCfrgtsB6JT6IWnH8ii0O9iUlsX0tfv5dMUebBr2LCIiHsbL3QVI5Yvq9Qj5v02kpWUPD7/7H74paIfd+eeNNrvD5PZL6ruxQhERkZJ0h+V8FBDOweauuXBed7zO115P8YTf1/QJT8WKg283aJ4WERHxLAos56m6fZ7gWFQCFsOkjWU3d/M17+Y9wTLf+8ndn8y+o3nuLlFERKSYAsv5yj+MsHvmwcPbIHEctLwB/EKJNo7xmNdU5m5Md3eFIiIixRRYzndB0dD2Vuj3MQxbiBMLV1jXkbxuubsrExERKabAIn+q2RBb074AXHF0MnsO57q5IBERERcFFinBt8ejAFxjWcnSVavcXI2IiIiLAouUFN2KA5HdsRomtTb8z93ViIiIAAoscgrBVzwOwGWFC9mza5ubqxEREVFgkVMIatyFLb5t8TEcHPvpTXeXIyIiosAip3a4nWuNoebpMzFzDrq5GhEROd8psMgptel+Hb85G+JHEUcWvO3uckRE5DynwCKnFOzvw5IY1/T9YUnjYe0k9xYkIiLnNQUWOa24TjfwteMSrKYDvnuAYzMfBafD3WWJiMh5SKs1y2ld0TKGfy1+nD2HonnYezphGz7g1y1J7Oz2XwJqhFBod1L0+6tFbDCdGtR0d8kiIlJNGaZpmu4uojxkZ2cTEhJCVlYWwcHB7i6n2rA7nCzadoidP3/GkIOv4GfY2OKM4x7bg+wxY4rbWS0GCx/uTt2agW6sVkREqprS/n7rkZD8Iy+rhZ4torh7+CPk3fotuT41aW7Zx4++o3g1ehFXtahF/VqBOJwmH/6y293liohINaXAIqUW3rQzgcN/gfrd8aGI/sc/YHzBE7zdwweAr9bs59CJQjdXKSIi1ZECi5RNSG0Y9A30fRd8gyFtHRfM6cuztRZSZHfy8bIUd1coIiLVkAKLlJ1hwIWDYPgqaNIbw2njtpwJtDV28tnKvZwosLm7QhERqWYUWOTsBcfCLV/ABTcD8GTALE4U2JmyKtXNhYmISHWjwCLnxjDg0sfBsNLBsY4Lje18tDSFQrvmaxERkfKjwCLnLrwBtL0FgMf8ZnHwRCEz1x1wc1EiIlKdKLBI+ej2KFi86GQm0d7YxvtLduNwVospfkRExAMosEj5CKsHbW8F4FHfGaQczmXOxnT31iQiItWGAouUn66PuO6ysJGLja2M+vo3Nuw77u6qRESkGlBgkfITVhfaDQBgdNA35BY5GPzxarZnnnBzYSIiUtUpsEj56voIWLxpWbSBuyOTOZ5nY8CEVaQeyXN3ZSIiUoUpsEj5Co1zTSoHPJH9Mt8GvkTT3F8ZMGElmdkFbi5ORESqKi93FyDV0BXPg+mA9ZO5wLGJz3w2sSG3AW+915/Qtn3p1KAWF9ULI8jP292ViohIFWGYplktxp6WdnlqqURZB2D5uzjXfozF7rq7sszRkuftA9lBPK1rhzAooR43XFgbwzDcXKyIiLhDaX+/FVik4uUepmDJO3j/Og6rswgHFj6z9+Q/9n+RRQ16No/k5etbExns5+5KRUSkkimwiOc5thd+fBq2fAtAvlcIDxUMY579QkL8vXn+upb0bROruy0iIueR0v5+l7nT7ZIlS+jTpw+xsa4fllmzZp1xn8mTJ9OmTRsCAgKIiYnh9ttv58iRI8WfT5o0CcMwSrz8/PRf29VOWF3o/xkM+hYimuNvz+J/fu/RN+oQWfk2HpiaxPAp67QOkYiInKTMgSU3N5c2bdowduzYUrVftmwZgwYNYujQoWzatIlp06axevVqhg0bVqJdcHAw6enpxa+9e/eWtTSpKhp0h7uXQuMrsdgL+K/xJqO6R+JlMZi7MUOrPYuIyEnKPEqod+/e9O7du9TtV6xYQb169bj//vsBqF+/PnfddRevvvpqiXaGYRAdHV3WcqSqsnrBDR/ABz0wjqVw16GXCOzzFk9/s4UJv6QwoFNdvK0adS8iIi4V/ouQkJDAvn37mDt3LqZpkpmZyfTp07n66qtLtMvJyaFu3brExcVx3XXXsWnTpoouTdzNPwxungzeAbD7Z27OnkStGj4cOJ7PnN+0DpGIiPypwgNLly5dmDx5Mv3798fHx4fo6GhCQkJKPFJq2rQpEydO5JtvvuHzzz/H6XTSuXNn9u/ff9rjFhYWkp2dXeIlVVBUS+j7LgBeK/7Li012ATB+8S6qSX9wEREpBxUeWDZv3swDDzzAM888w9q1a5k3bx579uzh7rvvLm6TkJDAoEGDaNu2Ld27d2fGjBlERETw/vvvn/a4Y8aMISQkpPgVFxdX0aciFaX1vyBhBAC9tj/LR75v0fHQdFb/uhwUWkREhHMc1mwYBjNnziQxMfG0bQYOHEhBQQHTpk0r3rZ06VK6du1KWloaMTExp9yvX79+eHl58cUXX5zy88LCQgoLC4vfZ2dnExcXp2HNVZXDDl/cDDvnl9xeI8q1oGK3R8Hb3z21iYhIhamwYc1llZeXh8VS8musVivAaW/5OxwONm7ceNowA+Dr60twcHCJl1RhVi+49Su4YyHZXZ5imbMVBaY35GTCL2/CuC6wd7m7qxQRETcpc2DJyckhKSmJpKQkAFJSUkhKSiI11TUUddSoUQwaNKi4fZ8+fZgxYwbjxo1j9+7dLFu2jPvvv58OHToQGxsLwPPPP8+PP/7I7t27WbduHQMGDGDv3r3ccccd5XCKUmVYLFCnPcFXPMbM1v/jgsIJvB/1DNSIhqO74OPeMOcRKDzh7kpFRKSSlTmwrFmzhnbt2tGuXTsARo4cSbt27XjmmWcASE9PLw4vAEOGDOGtt97ivffeo1WrVvTr14+mTZsyY8aM4jbHjh1j2LBhNG/enKuvvprs7GyWL19OixYtzvX8pIq6q1sDivDmldRmpNy8ENoNdH3w64fwv86QutK9BYqISKXS1Pzise745Fd+2nKQ69vV5j/928Kun+G7++F4KhhWuOI5V2fdU0zl73CaWC2a4l9ExNN5TB8WkbN1z6WNAJi5/gCfrdwLDXvAPcuh1b/AdLjWJfpyAOQfL7Hfh0t20/ipuczdqLlcRESqCwUW8Vjt64bxyJVNAHj2m2R+3nYQfIPgxglwzZtg9YGts+H9bpC2HoDf9h/nlXlbcZow5vst2BxOd56CiIiUEwUW8WjDezTiX+3r4DRhxOR1bE7Ldj0CuvgOuP0HCI2H43vhoyuxLXuPh6aux+F0PeXcdzSfWesPuPkMRESkPCiwiEczDIOXr29NQoOa5BY5GPrJr2RmF7g+rH0h3LUEml0LjiK85z/FU1mjaVqjgHsubQjAez/vxK67LCIiVZ4Ci3g8Hy8L4we0p2FEIOlZBQz95Ff2Hc1zfegfBv0/Z1eH5yk0vbnMmsS3Xo9zf719hAf6sPdIHt8kpbn3BERE5JwpsEiVEBLgzcdDOlAz0IfkA9lc+sYihk9ZR9K+42QV2BmwoRV9i14g068+vgWH8J/6Lz6J+ooa5PHezzuLHxOJiEjVpMAiVUZ8zQA+v6MjlzSqhcNpMue3dBLHLuPyNxeTnlVAQXgzgu5bAhfdDkDrtK/4ye9xGhxdwncbdJdFRKQq0zwsUiVtTstmwtLdfLchDZvDxGLAtLs7075umKvBrp9h9oNwbA8Ai7y60PX+iViDo91Ws4iInKy0v98KLFKlZWQV8PW6/TSMCOSqVn9be6ooj8IFL2NdNRYvnBT4ReA3YgXUiHBPsSIichJNHCfnhegQP4b3aHRyWAHwCcC394tMa/cpu53R+BUcYvcnd7L7oNYiEhGpahRYpNq7+sqreNbnUYpMKw0OLeSdt1/m8jcX8dq8reQU2t1dnoiIlIICi1R7If7evP3gQDY1uQeA570nkXsolf8t2sVTMze6uToRESkNBRY5L9Ss4Uu7m5+D2u0JNvKYWecLwDXS6OCJAneXJyIiZ6DAIucPqxckjgcvP2IOL2dUxArsTpMvV+9zd2UiInIGCixyfoloAj1HAzA07yP6WRexZeX32I/sAYfNnZWJiMg/0LBmOf84nfBpX9jzS8nthgUCIyGgJgTW/P1/I6HNza51i0REpNxpHhaRf5JzCJa8zp6t63Ae30cdyxF8OM0dFosXXPkSdLzLtVK0iIiUGwUWkVLYdzSPbq//DKaTRfe2pq53FuQdgdwjkHcYUpbAtrmuxi2vh77vgm+Qe4sWEalGNHGcSCnEhQdwWdNITCx89lsexLSBhpfBBf2YHXAdl+2/k5VNHnXdZdk0Ez7oAQe3FO+fV2Tn5blbmP2b1ioSEalICixy3huQUBeAaWv3k1/koMjuZPS3mxgxZT27j+QxIPlC9vSdDkGxcGQHfHgZbP8RgNd/2MYHS3YzYsp6Pl6WUurv3JKezUtzNpNyOLdCzklEpLpRYJHzXvfGEcSF+5OVb+PDX3Zz8wcrmLR8DwDx4QHYnSaPrPDFeecSqN8dbHkw9Rb2/PxxcTuA577bzIRfdp/x+9anHuOm91fw4S8pJI5dxvJdhyvozEREqg8FFjnvWSwGAzq67rK8NX8761KPE+TnxYeDLmLqnZ0I8LGyZu8xpm8tgAFfQ+ubwGmn3uIHGWyZR2LbWIb3aAjAi3O2/GNoWbPnKAM/Ws2JAjt+3hay8m0M+mg1U1enVsq5iohUVQosIkC/i+Lw8XL969AyNpg593XlihZRxIb681DPJgC8/P0WjhaYcP37bKh9MwCjvT/l5dBveeSKJtzfoyG+FDFuzko++ikJh7Nkf/aVu48waOJqcgrtJDSoyYonLqdvm1jsTpMnZmzkpTmbT9pHRERcNEpI5Hc/bz3IrkM5DOhUFz9va/F2m8NJn3eXsjXjBP3a1+GeSxty1X+XcKc5g0e8p7ka+YVAUR44XUOji0wrb5gD2VjnFtrFhxIT4sdLc7dQYHPStXEtPhh4Ef4+VkzT5J0FO/nPT9sBuKplNP/7vwuxWDR8WkTODxrWLFKO1u49xo3jlgPQICKQ3Ydy6dYkgk9aJ2PMfRhM5yn3e8eeyFv2foArgFzaNILxA9qXCEQA325I45FpGyiyO3n6mubc0bVBhZ6PiIinUGARKWejZmzki9/7mvh7W/nxoW7EhQfAiUzIPwa+NVxztPjUwLnkTSyLXgJgZdh1PO+8jWa1wxhzQ2t8vaynPP6UVak8OXMjPl4W5tx3CY2jNN+LiFR/modFpJw9flVTagb6APDwlU1cYQUgKAoim0FIHdejIYsVy6WPwTVvAQadjn3D3Nof89bVtfFNWwNrP4F5o+DrO+DwjuLj39IhjkubRlBkdzLyqw3YHKe+ayMicj7SHRaRMticlk3ygSxubF8Ha2n6mWyaCV8PK+7bcpLg2jB0PoTUBiAzu4Ar/7OErHwbD1zemIeuaFKO1YuIeB7dYRGpAC1ig7np4rjShRVwTef/f9PAL9T1Prg2NLwcOg2Hmo0h+wB8fgPkHQUgKtiPFxJbAfDezzv5bf/x8j8JEZEqSHdYRCqDwwa2fPD7y5/N46nw0ZVwIh3iOsHAmeDjesw0fMo65vyWTsOIQObc3/WkTroiItWF7rCIeBKrd8mwAhAa75qIzjcE9q2E6beDww7Ai9e1IiLIl12Hcnnjh21uKFhExLMosIi4U1RLuHUqWH1h+/cw5yEwTcICfXj1xtYATFyWwvrUY24uVETEvRRYRNytbmf410QwLLDuU9gwFYDLmkVxQ7vaOE14bPpvFNodbi5URMR9FFhEPEHza6HHk65/nvsIHNkFwL+vbUGtGj7sOJjDewt3urFAERH3UmAR8RSXjIS6l0BRDnw9FOxFhAX68MJ1rlFD4xbtYlNalpuLFBFxjzIHliVLltCnTx9iY2MxDINZs2adcZ/JkyfTpk0bAgICiImJ4fbbb+fIkSMl2kybNo1mzZrh5+dH69atmTt3bllLE6naLFa44X3XEOi09bDoZQB6t46hd6toHE4HX06ZiH3jTKgeg/tEREqtzIElNzeXNm3aMHbs2FK1X7ZsGYMGDWLo0KFs2rSJadOmsXr1aoYNG1bcZvny5dxyyy0MHTqU9evXk5iYSGJiIsnJyWUtT6RqC6kDfd91/fPSt2H3YnDYeLXxZub7jeL5nNF4fT0E5owEp/q0iMj545zmYTEMg5kzZ5KYmHjaNm+88Qbjxo1j165dxdveffddXn31Vfbv3w9A//79yc3NZfbs2cVtOnXqRNu2bRk/fnypatE8LFKtfPcArJ0EgZHg5QdZrjWMckw/AijEYphkxvYkfNCnePsFurdWEZFz4DHzsCQkJLBv3z7mzp2LaZpkZmYyffp0rr766uI2K1asoGfPniX269WrFytWrDjtcQsLC8nOzi7xEqk2er0MtZpA7kFXWAmMwLzsGd5qPZPhtvspNL2JSvuJTa9exnuzV7Nq9xGO5xW5u2oRkQrjVdFf0KVLFyZPnkz//v0pKCjAbrfTp0+fEo+UMjIyiIqKKrFfVFQUGRkZpz3umDFjeO655yqsbhG38gmE/pPh55eg3iXQbgCGtz/PAKk92jJ9QXP6bh5JW3MrgauHcMeyh9lrRhMZ5EuTqCA61A/nru4NTrsy9F8dyy1i0faD/LL9MG3jQxmUUK/CT09EpKwqPLBs3ryZBx54gGeeeYZevXqRnp7Oo48+yt13381HH3101scdNWoUI0eOLH6fnZ1NXFxceZQs4hkimsBNn5y0Ob5mAP930y3Y0lqT/+kNNC44wALfR5lm78bYnOtYeiKSpTsPs+dwLm/e1AbDOHndo6O5RUxfu4+fNh9kzd6jOH9/MDxj/QFaxgbTvm54RZ+diEiZVHhgGTNmDF26dOHRRx8F4IILLiAwMJCuXbvy4osvEhMTQ3R0NJmZmSX2y8zMJDo6+rTH9fX1xdfXt0JrF/Fk3rGt8L5nIXx7H167FnKL18/c7P0L22P6cOeeS5mxHuLCA05a8XnnwRMMnvgrB47nF29rFh1EgI+VdanHeWpmMrPvuwQvq2Y9EBHPUeF/I+Xl5WGxlPwaq9V1m/qP/r4JCQksWLCgRJv58+eTkJBQ0eWJVG0hdVyLJt7+IzS8DMO00zRtJgt9H+ZO63e8s2AbX6/dX9z81z1HuXHcCg4cz6dezQCe69uSpY/3YN6D3Zgw+GJCA7zZmnGCScv3uO+cREROocyBJScnh6SkJJKSkgBISUkhKSmJ1FTXKIZRo0YxaNCg4vZ9+vRhxowZjBs3jt27d7Ns2TLuv/9+OnToQGxsLAAPPPAA8+bN480332Tr1q2MHj2aNWvWMGLEiHI4RZHzQHxHV3AZOh8aXIrVtPOk9xdM9n6Zt2f8zIpdR5iXnM7/TVhFVr6NdvGhzLi3C4M716NOmGuF6PBAH0b1bgbAW/O3k/aXOzAiIu5W5mHNixYtokePHidtHzx4MJMmTWLIkCHs2bOHRYsWFX/27rvvMn78eFJSUggNDeWyyy7j1VdfpXbt2sVtpk2bxtNPP82ePXto3Lgxr732WomRRGeiYc0ivzNNWPcp5rwnMGx5ZJkBPGfeyUxbB0wTrmgRxTs3t8Pf5+QOuU6nSb/3V3A8NZmR0b9xTa2D0P1xqNPeDSciIueD0v5+n9M8LJ5EgUXkb47swjn9Dizp6wBIcjYkN6o9CV17YYm7GELjXe2KcqEwG/KPwY75FKz/Cr8jm/48jncA9JsETXpV/jmISLWnwCIi4LCRN/9l/Fa+jQVnyc+8A8BeAKbz5N0MK4vsFxDqbae9cyMYVuj7DrQb8I9fd6LAxqJth/DxstCr5ek7zYuI/EGBRUT+dDwV9q6A/b/CgTWQsRGc9j8/NyzgGwRRraH1jeQ2vJYrxm/kYFYOb/hMINGyBADbpU/j3f0R+MtQ6aO5Rfy0OZPvk9NZtvMIRQ5XAPpw0EVc0aLk/EoiIn+nwCIip2crgKz9rgnq/IJdd1v+Nl/Lmj1HeXT6b6QczuExry+51+tbAGb79+Vtr6FkF9jJLrBRYCt5hyYswJtjeTaig/2YP7IbQX7elXZaIlL1KLCIyDkzTZPkA9l8u+EA/usm8KB9IhbD5G37Dbxt/1dxu5axwVzVMpqrWkVTJyyAq/67hL1H8hjYqS4vJLZy4xmIiKdTYBGRcuV0muz7aSx1lz8FwIGur+BsN5iQAG+C/3YXZfnOw9w6YRWGAdPuSuCiepo5V0ROzWMWPxSR6sFiMah75Qjo5pq1uvbSp4g7tKRkWDm4BWbcReclA/i89tf0MZbxzvQfKLT9pb+Mww75x6Eor3JPQESqNN1hEZGyMU34ZjgkTXb1fRk829UXZvGrsGkmcPJfKQXWGvj5+rqGUNsLXBst3tDsGrhwEDToARb995PI+UiPhESk4jhs8MXNsPMn8A4EWx7FQaV5X2hyFWRs5OiOFQQe2YSvYfvn44XGQ7tB0OEO8A+r8PJFxHMosIhIxSrMgUnXQHqS633zPq5ZcaNbFzcxTZO7P1lJyrYNWCwGHZvG069zc1rVi4XD22DtJ/DbV1CY5dohPgFu+/6kEUsiUn0psIhIxcs7Cus/cz3SibnglE0OnijgoS+TWLbzSPG2i+uFMfSSBlzZIgqLowA2fwPfPeB6XHT7DxDfqbLOQETcTIFFRDxK8oEsJi5N4dsNadidrr92GkXW4J7uDenbNhbvOQ/Auk+hRSLc9MnZf5HD5pqZ9/c+MXaHk3Wpx7mgTgh+3ievnyQi7qXAIiIeKTO7gE+W7+GzlXs5UeAaPVQnzJ/HL3TQZ9m/XGHjgQ0QGnfyzvYisHoXPzJauDWTT1fs5d5LG9GhfjjkHoGPero6917xAoUtbuTeyetZsPUgrWoH89Hgi4kK9qvM0xWRM1BgERGPll1g4/OVe/nolxSO5BYB8KX/GDqaGynseB++vV8sucPhHfBJHwitCwNn8MvePG6f9Cs2h4mP1cJ/bmrDNduecD1e+t1WvzYMzxrALtO1MnxsiB8fDbmY5jH6O0LEU2geFhHxaMF+3tx7aSOWPn4Zz/VtSWyIHx8UXglA4aqJ/HduEkdyCl2N7YUw/XY4kQ77VnJsyjDu+mwNNodJRJAvRQ4nP371niusWLxwdLiHQsOXZgUbmOfzBN+3XEDzWl6kZRXQb/wKFm076MYzF5GzocAiIm7l72NlcOd6LH6sB71vGMIBI5pgcjm0bBJdXl3Is98kc2L2U5DxG/iFYlq8CdszhyGOGXRtXItfHuvB8Av9ed7rYwAWx9zGXYf+xeUFr7HQ2R5vw0HzXR8xO+AFro23kVNoZ+gna/hgyS6y8s4w3FpEPIYeCYmIR3GuHI9l3uPss9ahW+4rXGrZwMc+rwOwsdv7fLdsPU86xuPEoKjfFPxa9Mb8/AaMXQtJcjbgxqLncGDF18vChMEX0dXxK3x7H+QdxgyoybiIf/PatkgArBaDi+qG0bN5FD2aRRIW4E2h3UmR3Umh3UmIvzfRIerzIlKR1IdFRKqmwhPwVgsozGZnp1eIXD2GYGcWH9t78Zx9MADvBn9Kn6J54BsM7QfD8nfBy48fu07j3h9ysFgMJgy6iG5NIlzHzNoPU2+F9A2YhpXVTR7m3+ld2H4wlwAKqG+kU9fIxAc7BibG75PgFeBDnfrNuO6yS2jZsK67rohItabAIiJV17wnYeXY4rf5NVvyRNhbfLfpCHXCAvjqjvZEz7oJUlf8uc9Vr0Knu9lzOBerxSAuPKDkMW358O39sPEr1/voC7DnHMYrJ61UJeUaNSC8PoEX3QId79FSAiLlRIFFRKquoynwTjvAdK1XdNcSqNWYIzmFBPh44e9jhZyD8MGlkH0A6neDgd+cOUSYJqx4D+Y/A6bzz+0BNTHDG4FP4O8jpg0wDPJzjlN0OIUQ+5ESh3HUvQTrDe9DSJ1//j6nA/b8ApmbXXeCfALP4mKIVG8KLCJStX19B2ycBonjoO2tp25zZBckfw0X3Q6BtUp/7IyNrld4Q6jVGALC/7H5vozDzFi4jKwtC3jYMpVAoxCnbwiWvv+FlteXbGyacGAdJE931ZaT6dp+4WDo+07paxQ5TyiwiEjVZi+EExkQ5jl9R9buPcaLn37Hs7b/0Nay27Wx6TXgGwS5hyD3oKvm3EN/7uQXAgVZYFjg7mUQ1cI9xYt4KAUWEZEKsP9YHnd+vJLeRz/hXus3WI2T/wq1W/3ZHtqVhV5d+SanOY/lvsEVrMTR4DKsg2ae8rg2hxNvq/rFyPlHgUVEpIKcKLBx3xfrydq+nJ7WtWSbgRw2QzhCMIfMEHabMeTz53DoeCOTn3wewcdwsKLzB3S64iYMw6DA5uD75HS+/HUfq1OOMiihHk9f0xwvBRc5jyiwiIhUILvDyRs/bufrdfsB8LIYWH9/RdTwpXFUEE2iatAkKojjeTayv32cW+zfsM1Zh9Gx79MwOoRvktKK11P6Q9fGtXjvlgsJCfB2x2mJVDoFFhERD1KQfQTznXb427N40jaUKY7LAdfCj/3axxEb6scz32wi3+agQa1APhpyMfVraVSRVH8KLCIinmbleJj3ODleYbzccApXX9SEzg1rYsGEnAy2px3m6a83cCQnnxBfKw/1u5yuLeu5u2qRCqXAIiLiaRw2+F8nOLITGl/pGl10aDsc2QH2gpOaZ5kBjA5+gTadLue6trUJC/RxQ9EiFUuBRUTEE22d41om4O8MK3j5gcWKaVgoKirE15nPcTOQ/kX/JsVSjytaRjHyiiY0jKhR+XWLVBAFFhERT2Sa8MsbcGQ3RDSBWk0hoimE1gWr15/tinKxT7oOr7RfOWqEcX3Bv9lrRhMZ5Mus4V2Ite1zzdi7byX4BIFfMA6fYPbmemGN70DdHredeSZeEQ+gwCIiUtXlH4NJfSBzI0U16jDEeJ7kQ3ZGh8zmettcDKf9tLuaGBj1LoE2t0CLvq7HT8DavUd5dNpv1AkP4NUbWxMT4l9ZZyNySgosIiLVQc5BmHgVHN2FPaQuJ7KOEUY2AM7GV2Hp9ghZBTZe++ZXjh89TJQliystq+lk2fLnMXxqYPZ9j0+z2/HC7M3Yna6/9sMCvHmjXxsubx7ljjMTARRY3F2OiEj5Ob7PFVqyXXO+7DRr85xtILXbX8ODPZsw4KNV7DyYQ60aPnx6e0e+STrAnF9WcZ1lGbcHraRmQSoAb9j68Z4jkd6tYth3LI/kA67gc3uX+jzRuxk+XhYKbA52ZOawNSObo7lF2J0mNocTu8PE38dKv4vqEBnkd9pSRcpKgUVEpDo5vBMWPg91u7Ag8BqGTd6A04QgPy9OFNiJCfHj8zs6FnfIHb94F698vxULTp73+4IBzAFgZ/TVNBw6kSLDh1e/38bEZSkANKgVCAbsOZyL04SWxh7ijUwyzTAOEsZBM5QivGlQK5BpdydQs4Zv6erOO+paX+kvayjtOZyLYUDtUH/N6isKLCIi1dmnK/bwzDebAIgPD2DyHR2JCw8o0ebLX1MZNWMjThPuDFjME+ZHWEw7xHWE/pOhRgQ/bc7kkekbOJ5nA8CKg6f9pnMb35z0nUcJYb2jAalBbbnpxv4E1rsIrP8wI2/aevj8X5B3GHq/Bh3vYl5yOvdMXodpgrfVID48gPq1anBBnRDu6FqfAB+v0x9PqqUKCyxLlizh9ddfZ+3ataSnpzNz5kwSExNP237IkCF88sknJ21v0aIFmza5/mUbPXo0zz33XInPmzZtytatW0tdlwKLiJxvJvyym3Wpx3i2T0uigk/9mGbJ9kP8vO0gd3dvSNThlfDVINfq0f5hkDACOt5FZqE3i7cdIs7nBO1/fRif/ctdO8e0gfzjcCIdHEUnHdv0DsBo0AMufwYim5X8cNfP8OUAKMop3pR16Yt0W9yUrHwbVouBw1ny5+eSRrWYMPgi/LytpzyX3EI7gb4KNNVNhQWW77//nmXLltG+fXtuuOGGMwaWrKws8vPzi9/b7XbatGnDfffdx+jRowFXYJk+fTo//fRTcTsvLy9q1apV6roUWERESuHwDvhyIBz6vVOufxh0vg+iL4BvhkNOJvjUgL7vQqsbXG1M0zVi6WgKaRsXsXnl91zIFsKN38OIYYUOw+DSJ1zHS/4aZtwFThvU7+Y69or3AHjWNpikmJv46u4EDp0oZM/hPLZnnuDNH7eRW+TgsmaRjB/QHh+vPx8VHTpRyGPTN7Bo+yHuuKQ+j1/VTI+SqpFKeSRkGMYZA8vfzZo1ixtuuIGUlBTq1q0LuALLrFmzSEpKOttSFFhERErL6XCFisWvumbd/auI5tD/M6jV+LS7L991mNsmrqKBcy+v1ZxD65ylrg/8w6H5tbDuM8CEFolwwwdg9WHNxIe4aN/HABzp+gI1L78fcg9D5ibI3MS+jEzuXxfFels8vVvF8O4t7fCyWli8/RAPf7WBwzmFxd/fuWFN3r2lXen70ZQDh9PkeF5RpX7n+cJjA0ufPn0oLCzkxx9/LN42evRoXn/9dUJCQvDz8yMhIYExY8YQHx9/2uMUFhZSWPjnH+Ds7Gzi4uIUWERESsth/zO4HN0FF9wM174FPmdedHFecgb3Tl6L04Q+Nbbxkv9kgk/8JfxcPAx6vwoWK+tSj9Fv/HJGWqYy3Otb1+eBEZB76KTj7jDrMNPemYKm12MJr8eEpa5OwU2jgri5Qxyv/7CNvCIH3YIzebPZNiLqNHLd3alA2zNPcP8X69mWeYKHejbhvssaYRhGhX7n+cQjA0taWhrx8fFMmTKFm266qXj7999/T05ODk2bNiU9PZ3nnnuOAwcOkJycTFBQ0CmPdap+L4ACi4hIWTnscCINQk//H4mnsnBrJs98s4n9x/Kx4uCxWisY5PUT/hfdCl0eBMPgRIGNq9/5hX1H87muTQxvR3yHsfSt349gQFg9iGoJhgHbfwTHn/8huslZl7XOJgQ0TKDPtdfjG1STzOVTOLr0I5o7d/xZyPXvQ5ubS1d0QTZsnAbZB1wjmPKPuv43JA56vQQB4cVNTdPk81WpvDh7M4V2Z/H2a1rH8Hq/C9RBuJx4ZGAZM2YMb775Jmlpafj4nH4Rr+PHj1O3bl3eeusthg4deso2usMiIuJ+BTYHHy1N4X8/7yS3yAGAr5cFL4uB1WLgNCGn0E6dMH/mPtCVYF8v2LscvHwhohn4/mVdpIIs2PIdh5d/RvjBVViMv/08GRYwXcHBjpXNzngusKTgsPphvXOhK/j8hWmaTFq+h2N5Nu7v0QCv5K9g/rOQe/DUJ1P7Ihj0DfjW4GhuEY9N/42ftmQC0L1JBF0b1+LVeVuxOUxaxgbzwaCLqB2qmYLPVWkDS6XFQ9M0mThxIgMHDvzHsAIQGhpKkyZN2Llz52nb+Pr64uurZ4kiIu7k521leI9G9Gtfh9d/2Mb0dfsptDsp/EsbH6uF/97clmC/34dA1+tymoOFQLsB1Go3gPQDewk5tIaAjLWwbxWkb3B14q3VFC4ciKV1f76cn86x9SPozm8UThmA7z2Lwc/1g2eaJi/O2cJHS1NoZezmxnVfUDffNTKV8IbQqKerg3BAOHj7u9ZlOrCGrE/6Mzb6JWb8dojDOUX4WC083rsZt3Wuh8VicEGdUO7+fC2b0rIZ9O5cXr+1Exc2rF1xF1iKVdodlkWLFtGjRw82btxIq1at/rFtTk4O8fHxjB49mvvvv79UtajTrYiI+2Xl28gptONwmNicThxOk9AA73OfHdeW7xrBFFrX9fgIsDmc3DfhJ/6ddje1jSMUNO6D362fYQIvz93C4qVLuNNrDjdYfsFimBRZAvC57HHodC94/fkfzr/tP87SRfMYsvMBAihgjqMD99nup15EEO/c3I5WtUNKlLL/0BFWTHiEGwpmspn6ZPabQ89Wsed2fuexCnsklJOTU3zno127drz11lv06NGD8PBw4uPjGTVqFAcOHODTTz8tsd/AgQPZsWMHK1euPOmYjzzyCH369KFu3bqkpaXx7LPPkpSUxObNm4mIiChVXQosIiLnn6x8G0++O5H/5I7Cx3Bgu+Ilvk6tQczmj+hu/a243QzHJbxiu4XhfS9hcOd6ABzPK+LVeduY+msqpgldLBv52Od1fLCzr34/Im4Zj9/f+6mkrnIN/z7yZx+aR+z30OXGEVzfruTq2D9tzuTl77cQ6u/N+wMvIiJITwVOpcIeCa1Zs4YePXoUvx85ciQAgwcPZtKkSaSnp5Oamlpin6ysLL7++mv++9//nvKY+/fv55ZbbuHIkSNERERwySWXsHLlylKHFREROT+F+Hvz6O0DeGvsVp4wJ+I9/yluBrCCEwuWFn2h8/3s3xbMwfnbGf3dJmrW8CG/yMGY77dyNNc1Id51bWP5V/sOWAqbw9dDiEuZBh9vh8gWrkdINRvC/l9h5TjAhBrROOt2xrJpBg9YpnPZlwkcz7NxW5f67Duax3PfbS7u/wJw0/sr+PyOjurzcg40Nb+IiFR5K3Ye5sgnA7jWuoIc04/Uuv+ixfWPukYh4erT8sw3m/hs5d4S+zWJqsGLia3pUP/P0UGs+xS+vR84zc9j2wGuEUVWH8x32mLkZPK07TY+d1xBr5ZRLN5+iAKbEy+LwaCEevywKYMDx/OJDfHjs7+s9yQuWktIRETOK7PWprB43gw6dr2Cm7tdcNLnDqfJiCnr+D45A39vKw/2bMztl9TH+1Sz5h7d7eroe2SX63V0l2uEUvcnoHHPP9ut/hDmPkKud03an3iDAlyPfTrWD+fFxFY0jgoi7Xg+Az5axe5DudQM9OHToR1oGRty8neepxRYRERE/qbI7mT+5kzaxYcSWx6PZ+xF8F57OJ5KUrORjD7Sk0EJdbm+Xe0Sk8sdySlk0MTVbErLJsjPi//e3JbLmkXBsb2uZQsKc+DCgRCfQJ7Nwc9bD9GpQfhpZ9bddzSPV+ZtpXXtEG7vUr/EUgZ/lZFVQIi/N/4+p16fyRMosIiIiFSGpCkw6x7XMOkHNriGZ59CdoGNoZN+5dc9x4jiKG9Gz6dL9lwMp624zfGQ5vw3pyeTcy+iZkgQHw666KRRStsyTjBo4ioys12DxxtGBPLCda3o3OjP9feSD2Tx7sId/LApk7o1Xat51wkruZq3p1BgERERqQxOB/wvAQ5vg+6PQ48n/9yetd+1ZlJhNhSewJaXxfrVS2iTORNfwxVUMiM6YwmrS/D2Gfj+PoPNITOED+zX8JWlN6/d3IFeLaMBWLv3GLdP+pWsfBsNIgLJzrcXr7PUt00sN1xYm0+W7+HnbSWXPagd6s8XwzoRX7OUocXpBEvlLDCpwCIiIlJZNn8DXw1yrXTd4NLf+73sLrHUwN8lW1vwQt6NrDKbAxDKCQb7LGKY3wJqFLlm4z1g1uRt+400umIYTWJCuffzdeTbHFwYH8rEIRdjGAZv/riNz1buxTQhmFzaWXbS3rKDK4NTaWjfyXJHc27PvZfIkECmDOtE/Vr/sFZU7mH44SnXGlO1GkPjK6FJL6jTAawVM9esAouIiEhlMU344FJITyq53eoDNaLANxh8g1wv/zBoczO2epfywS8p/HfBDmwOJze1j+PhXk2IDLDCb19i/vwSRnYaANucdXjXfj0rnC1o1aQR4wZc+OdaRjmHOLD8C7JXT6G5fcspy5vpcy0PZd9KZJAvU4Z1omFEIGlZBezIPEHK4VyaRtWgc+4CmDfKtb7S30/PL5Ts2t0JueY5CK9fjhdOgcXd5YiIyPnm8E5I+hyCYlzzttRs5FpU0fLPHV7Ts/IptDmp9/c7H7Z8WP0hhT+/jq89u3izGRKPUac9RLaEfSth189gOv7cL6w+xHWAOhe7ZgWe8zAA7/rfw5vHuhLoY8UE8n5f+6mOcZCXvT6im3Wja//IlnDVGNdq2tt/wL79R7wKj+PEIPveZEIjS06Qd64UWERERKqD/ONkfP8aAbu/Jzhn96nbxF4IrftBy0QI/tsyAb+8CQuexzSsPF1jNJMPNQQgxFrIo4Hfc5NtFj5mEYWmN1/VuJXLb3+R2JrBFNgcvPHDNj5etou27KBz4AF63fbMSZ2Az5UCi4iISHVTkAVp6+HAWsjcBLWaQKt/Qa1Gp9/HNGHm3fDbVEzfYH7rNZ3o3M1Ern4V40Q6AEcjOzH44M1sLIgkPNCHB3s2ZtKyPew+nAvATRfV4elrW/y5gGU5UmARERERF3shfNLHtfK1xdu18jW4FpO88kVo3od9x/K5Z/Jakg/8+fgpKtiXV264gB7NIiustNL+flfOmCURERFxHy9f6D8ZQuNdYcWnBlz+LAxfDS36gmEQFx7A9Ls7c2vHeKwWgxsurM2PD3av0LBSFrrDIiIicr44ngpbZkOrGyAo+rTNCu0OfL0qZ3bcClutWURERKqo0HhIuPeMzSorrJSFHgmJiIiIx1NgEREREY+nwCIiIiIeT4FFREREPJ4Ci4iIiHg8BRYRERHxeAosIiIi4vEUWERERMTjKbCIiIiIx1NgEREREY+nwCIiIiIeT4FFREREPJ4Ci4iIiHi8arNas2magGuZahEREaka/vjd/uN3/HSqTWA5ceIEAHFxcW6uRERERMrqxIkThISEnPZzwzxTpKkinE4naWlpBAUFYRjGWR8nOzubuLg49u3bR3BwcDlWKH+na115dK0rj6515dG1rjwVea1N0+TEiRPExsZisZy+p0q1ucNisVioU6dOuR0vODhY/wJUEl3ryqNrXXl0rSuPrnXlqahr/U93Vv6gTrciIiLi8RRYRERExOMpsPyNr68vzz77LL6+vu4updrTta48utaVR9e68uhaVx5PuNbVptOtiIiIVF+6wyIiIiIeT4FFREREPJ4Ci4iIiHg8BRYRERHxeOdlYBk7diz16tXDz8+Pjh07snr16n9sP23aNJo1a4afnx+tW7dm7ty5lVRp1VeWa/3hhx/StWtXwsLCCAsLo2fPnmf8/0b+VNY/13+YOnUqhmGQmJhYsQVWI2W91sePH2f48OHExMTg6+tLkyZN9PdIKZX1Wr/99ts0bdoUf39/4uLieOihhygoKKikaquuJUuW0KdPH2JjYzEMg1mzZp1xn0WLFnHhhRfi6+tLo0aNmDRpUsUWaZ5npk6davr4+JgTJ040N23aZA4bNswMDQ01MzMzT9l+2bJlptVqNV977TVz8+bN5tNPP216e3ubGzdurOTKq56yXutbb73VHDt2rLl+/Xpzy5Yt5pAhQ8yQkBBz//79lVx51VPWa/2HlJQUs3bt2mbXrl3N6667rnKKreLKeq0LCwvNiy66yLz66qvNpUuXmikpKeaiRYvMpKSkSq686inrtZ48ebLp6+trTp482UxJSTF/+OEHMyYmxnzooYcqufKqZ+7cueZTTz1lzpgxwwTMmTNn/mP73bt3mwEBAebIkSPNzZs3m++++65ptVrNefPmVViN511g6dChgzl8+PDi9w6Hw4yNjTXHjBlzyvY33XSTec0115TY1rFjR/Ouu+6q0Dqrg7Je67+z2+1mUFCQ+cknn1RUidXG2Vxru91udu7c2ZwwYYI5ePBgBZZSKuu1HjdunNmgQQOzqKioskqsNsp6rYcPH25edtllJbaNHDnS7NKlS4XWWd2UJrA89thjZsuWLUts69+/v9mrV68Kq+u8eiRUVFTE2rVr6dmzZ/E2i8VCz549WbFixSn3WbFiRYn2AL169Tpte3E5m2v9d3l5edhsNsLDwyuqzGrhbK/1888/T2RkJEOHDq2MMquFs7nW3377LQkJCQwfPpyoqChatWrFyy+/jMPhqKyyq6SzudadO3dm7dq1xY+Ndu/ezdy5c7n66qsrpebziTt+G6vN4oelcfjwYRwOB1FRUSW2R0VFsXXr1lPuk5GRccr2GRkZFVZndXA21/rvHn/8cWJjY0/6l0JKOptrvXTpUj766COSkpIqocLq42yu9e7du1m4cCH/93//x9y5c9m5cyf33nsvNpuNZ599tjLKrpLO5lrfeuutHD58mEsuuQTTNLHb7dx99908+eSTlVHyeeV0v43Z2dnk5+fj7+9f7t95Xt1hkarjlVdeYerUqcycORM/Pz93l1OtnDhxgoEDB/Lhhx9Sq1Ytd5dT7TmdTiIjI/nggw9o3749/fv356mnnmL8+PHuLq3aWbRoES+//DL/+9//WLduHTNmzGDOnDm88MIL7i5NysF5dYelVq1aWK1WMjMzS2zPzMwkOjr6lPtER0eXqb24nM21/sMbb7zBK6+8wk8//cQFF1xQkWVWC2W91rt27WLPnj306dOneJvT6QTAy8uLbdu20bBhw4otuoo6mz/XMTExeHt7Y7Vai7c1b96cjIwMioqK8PHxqdCaq6qzudb//ve/GThwIHfccQcArVu3Jjc3lzvvvJOnnnoKi0X/jV5eTvfbGBwcXCF3V+A8u8Pi4+ND+/btWbBgQfE2p9PJggULSEhIOOU+CQkJJdoDzJ8//7TtxeVsrjXAa6+9xgsvvMC8efO46KKLKqPUKq+s17pZs2Zs3LiRpKSk4lffvn3p0aMHSUlJxMXFVWb5VcrZ/Lnu0qULO3fuLA6FANu3bycmJkZh5R+czbXOy8s7KZT8ERRNLZtXrtzy21hh3Xk91NSpU01fX19z0qRJ5ubNm80777zTDA0NNTMyMkzTNM2BAweaTzzxRHH7ZcuWmV5eXuYbb7xhbtmyxXz22Wc1rLmUynqtX3nlFdPHx8ecPn26mZ6eXvw6ceKEu06hyijrtf47jRIqvbJe69TUVDMoKMgcMWKEuW3bNnP27NlmZGSk+eKLL7rrFKqMsl7rZ5991gwKCjK/+OILc/fu3eaPP/5oNmzY0LzpppvcdQpVxokTJ8z169eb69evNwHzrbfeMtevX2/u3bvXNE3TfOKJJ8yBAwcWt/9jWPOjjz5qbtmyxRw7dqyGNVeEd99914yPjzd9fHzMDh06mCtXriz+rHv37ubgwYNLtP/qq6/MJk2amD4+PmbLli3NOXPmVHLFVVdZrnXdunVN4KTXs88+W/mFV0Fl/XP9VwosZVPWa718+XKzY8eOpq+vr9mgQQPzpZdeMu12eyVXXTWV5VrbbDZz9OjRZsOGDU0/Pz8zLi7OvPfee81jx45VfuFVzM8//3zKv3//uL6DBw82u3fvftI+bdu2NX18fMwGDRqYH3/8cYXWaJim7pOJiIiIZzuv+rCIiIhI1aTAIiIiIh5PgUVEREQ8ngKLiIiIeDwFFhEREfF4CiwiIiLi8RRYRERExOMpsIiIiIjHU2ARERERj6fAIiIiIh5PgUVEREQ8ngKLiIiIeLz/B/NP71wHtg+CAAAAAElFTkSuQmCC" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(c_sticking_probability.p, c_sticking_probability.fd_mean)\n", "plt.plot(rust_sticking_probability.p, rust_sticking_probability.fd_mean)\n", "plt.legend([\"C\", \"Rust\"])" ], "metadata": { "collapsed": false } }, { "cell_type": "code", "execution_count": 125, "outputs": [ { "data": { "text/plain": "" }, "execution_count": 125, "metadata": {}, "output_type": "execute_result" }, { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAA9hAAAPYQGoP6dpAABUiklEQVR4nO3de3xTVbo38N/eadILbVMKlLbSUkCxFBQVqAKiwAEpKsh7HGGGEXDEG6BY9eiAo7QwM4KDOuCIqKMcBEVAVERAEGWQo4AFhBlqKYpQirYFAXujtEmz1/tHmrS5tUmae37fzyfndO+snb2yZejDWs96liSEECAiIiLyE9nfHSAiIqLwxmCEiIiI/IrBCBEREfkVgxEiIiLyKwYjRERE5FcMRoiIiMivGIwQERGRXzEYISIiIr+K8HcHnKEoCsrKyhAXFwdJkvzdHSIiInKCEAI1NTVITU2FLDse/wiKYKSsrAxpaWn+7gYRERG54fTp0+jWrZvD94MiGImLiwNg/DLx8fF+7g0RERE5o7q6Gmlpaebf444ERTBimpqJj49nMEJERBRk2kqxYAIrERER+RWDESIiIvIrBiNERETkV0GRM0JERORtBoMBer3e390IKiqVChEREe0uu8FghIiIwl5tbS1++uknCCH83ZWgExMTg5SUFGg0Grc/g8EIERGFNYPBgJ9++gkxMTHo0qULi2s6SQgBnU6HX375BSdPnsQVV1zRamGz1jAYISKisKbX6yGEQJcuXRAdHe3v7gSV6OhoqNVqnDp1CjqdDlFRUW59DhNYiYiI0HYtDLLP3dEQi8/wQD+IiIjCXp2uERlztiBjzhbU6Rr93Z2gwmCEiIiI/IrBCBERkQcYlOaVOAUnL1gcU+sYjBAREbXTtsJyjHrpS/PxPf+7Hzc+vxPbCsu9fu+Kigo88sgj6NmzJyIjI5GWloZx48bhiy++8Pq9PYXBCBERUTtsKyzHjHe+xZnqBovzFVX1mPHOt14NSEpKSjBgwADs3LkTixcvxpEjR7Bt2zaMGDECs2bN8tp9PS1sl/bW6RqRNW87AKBowRjEaML2URARkZsMisD8T4pgb0JGAJAAzP+kCKOzkqGSPb9aZ+bMmZAkCQUFBejQoYP5fN++fXHvvfd6/H7eEr4jI7qLKImajJKoyYDuor97Q0REQajg5AWUV9U7fF8AKK+qR8HJCx6/94ULF7Bt2zbMmjXLIhAxSUhI8Pg9vSV8gxEiIqJ2OlvjOBBxp50rjh8/DiEEMjMzPf7ZvsZghIiIyE1Jcc5VHHW2nStCaR8dBiNERERuyu6RiBRtFBxlg0gAUrRRyO6R6PF7X3HFFZAkCcXFxR7/bF9jMEJEROQmlSwhb1wWANgEJKbjvHFZXkleTUxMxJgxY7Bs2TJcvGib+1hZWenxe3oLgxEiIqJ2yOmXguV3X4ek+EiL88naKCy/+zrk9Evx2r2XLVsGg8GA7OxsfPDBB/jhhx9w9OhRvPzyyxg8eLDX7utpXM9KRETUTjn9UjD08s64Kv8zAMDKPwzCsCu6eGVEpKWePXvi22+/xV//+lc88cQTKC8vR5cuXTBgwAAsX77cq/f2JAYjREREHtAy8Mjukej1QMQkJSUFr7zyCl555RWf3M8bGIwQERF5QIwmAiWLbvN3N4ISc0aIiIjIrxiMEBERkV8xGCEiIiK/YjBCREREfsVghIiIiPyKwQgRERH5FYMRIiIiT9BdBPK1xpfOtjw7OcZghIiIiPyKwQgREZEnKIbmn0/tsTz2gnvuuQeSJEGSJKjVavTo0QNPPfUU6uvrPfL5GRkZWLJkiUc+qy0MRoiIiNqraBOwLLv5+N3fAEv6Gc97UU5ODsrLy3HixAn8/e9/x+uvv468vDyv3tMbGIwQERG1R9EmYP1UoKbc8nx1ufG8FwOSyMhIJCcnIy0tDRMmTMCoUaOwY8cOAPZHNq655hrk5+cDAIQQyM/PR3p6OiIjI5GamorZs2cDAIYPH45Tp07hscceM4++eBODESIiIncpBmDbHwEIO282nds2x+tTNgBQWFiIPXv2QKPRONX+gw8+MI+m/PDDD9i4cSOuuuoqAMCHH36Ibt26YcGCBSgvL0d5eXkbn9Y+4btRXos/GHLpXiBzNCCr/NghIiIKOqf2ANVlrTQQQPXPxnY9hnn89ps3b0ZsbCwaGxvR0NAAWZad3r23tLQUycnJGDVqFNRqNdLT05GdbZxqSkxMhEqlQlxcHJKTkz3eb2vhOTJStAlRbwwxH0atn+STuT0iIgoxtWc8285FI0aMwOHDh/HNN99g2rRp+MMf/oA777zTqWvvuusuXLp0CT179sT999+Pjz76CI2NjV7pZ1tcCkYWLlyIQYMGIS4uDklJSZgwYQKOHTvW6jXfffcd7rzzTmRkZECSJJ9l5jrUNLcn1VoOOQkfzO0REVGIie3q2XYu6tChAy6//HL0798fK1aswDfffIO33noLACDLMoSwnD7S6/Xmn9PS0nDs2DG8+uqriI6OxsyZM3HTTTdZtPEVl4KRL7/8ErNmzcK+ffuwY8cO6PV63HLLLbh40XFxl7q6OvTs2ROLFi3yyVBPq5rm9gQErFNxJAjj7J6P5vaIiCgEdB8CxKcCNr9VTCQg/jJjOy+TZRlPP/00nnnmGVy6dAldunSxyPWorq7GyZMnLa6Jjo7GuHHj8PLLL2PXrl3Yu3cvjhw5AgDQaDQwGHzz+9ClYGTbtm2455570LdvX/Tv3x8rV65EaWkpDh486PCaQYMGYfHixfjtb3+LyMjIdne4XZrm9hz/kWkxt0dERNQWWQXkPN90YPvPXABAziKf5STeddddUKlUWLZsGUaOHInVq1fj//7v/3DkyBFMmzYNKlVzP1auXIm33noLhYWFOHHiBN555x1ER0eje/fuAIyrcXbv3o2ff/4Z586d82q/25UzUlVVBcCY6BIMlJoK59qtvJ2lfImIyDlZ44GJq4A4q9H/+FTj+azxPutKREQEHn74Yfztb3/DnDlzcPPNN+P222/HbbfdhgkTJqBXr17mtgkJCfjnP/+JoUOH4uqrr8bnn3+OTz75BJ06dQIALFiwACUlJejVqxe6dOni1X5LwnpCyUmKomD8+PGorKzEV1995dQ1GRkZyM3NRW5ubqvtGhoa0NDQYD6urq5GWloaqqqqEB8f7053AQDffb0FfXdMbrudIR19n90HaDq4fS8iIgoO9fX1OHnyJHr06IGoqKh2fFA1sCjN+PPvNwC9RobFKs3Wnl91dTW0Wm2bv7/dHhmZNWsWCgsLsXbtWnc/wqGFCxdCq9WaX2lpaR753OMxV6FMJEJxEH4pAigTnXAcqR65HxERhZGWgUf3IWERiHiKW8HIww8/jM2bN+Nf//oXunXr5uk+Ye7cuaiqqjK/Tp8+7ZHPTYrvgPn6qQBgE5CYjufrpyAJVajT+Wd5ExERBSlNByC/yvjiyLpLXApGhBB4+OGH8dFHH2Hnzp3o0aOHVzoVGRmJ+Ph4i5cnZPdIxH/ibsJMfS4qYJnnUoFOmKnPxX+UnsiWiz1yPyIiImqbSxVYZ82ahTVr1uDjjz9GXFwcKiqMCaFarRbR0dEAgKlTp+Kyyy7DwoULAQA6nQ5FRUXmn3/++WccPnwYsbGxuPzyyz35XdqkkiXkjcvCjHfqsaNhIAbJxUhCJc4iAfuVTCiQsVz9d6gkt9JoiIiIyA0ujYwsX74cVVVVGD58OFJSUsyvdevWmduUlpZarGsuKyvDtddei2uvvRbl5eV44YUXcO211+K+++7z3LdwQU6/FCy/+zp0jovCPiULm5Qh2KdkIUkbgyW/yUSOar9f+kVERBSuXBoZcWbhza5duyyOMzIynLrOl3L6pWBoejQKF+fgLBKgnfQ6hmWlo6GuGtjs794REZE/BNrvqmDhiecWthvlqWQJg1VHAQB1GQlQyd7dHpmIiAKTqRCYTqczpxyQ8+rq6gAAarXa7c8I22CEiIgIMBYKi4mJwS+//AK1Wg1ZDs89ZF0lhEBdXR3Onj2LhIQEi+qurmIwQkREYU2SJKSkpODkyZM4deqUv7sTdBISEtq99xyDESIiCnsajQZXXHEFdDqdv7sSVNRqdbtGREwYjBAREcG46227ysGT2zgxRkRERH7FYISIiIj8isEIERER+RWDESIiIvIrBiNERETkVwxGiIiIyK8YjBAREZFfMRghIiIiv2IwQkRERH7FYISIiIj8isFIS4rB/KNcutd4rLsI5GuNL91FP3aOiIgoNHFvGpOiTYja+pT5MGr9JCA+FRi1wI+dIiIiCn3hG4xoOiCjfg0A4NiJncCHf4AEYdmmuhz48H4/dI6IiCh8hP00jQwF6h1PAxCQbN4VthcQERGRR4V9MJItF0OuKWulBQMSIiIibwr7YCQJlf7uAhERUVgL+2DkLBL83QUiIqKwFrbBSIwmAiWLbsO7Cx7DGXSC4mA2RhFAhUiAgbM1REREXhG2wYhJwakqzNNNAQCbgMR0nKe/BwVKHx/3jIiIKDyEfTBytqYe25VszNDnogKJFu9VoBNm6HOxXcnGWSSgTtfop14SERGFrvCtM9IkKS4KALBdycaOhoHIlouRhEqcRQIKlEwoTfEaE12JiIi8I+yDkeweiUjRRqGiqh4KZOxTsizelwAk4xyy5WI0+KeLREREIS3sp2lUsoS8ccYAxLromek4T70aKokZrERERN4Q9sEIAOT0S8Hyu69DUpzG4nyyNgpLfpOJHNV+P/WMiIgo9IX9NI1JTr8UDE2PRuHiHJxFArSTXsewrHQ01FUDm/3dOyIiotDFYKQFlSxhsOooAKAuIwEq2Xa3GiIiIvIsTtMQERGRXzEYISIiIr9iMEJERER+xWCEiIiI/IrBiIvqdI3ImLMFGXO2sDw8ERGRB3A1jat0F1ESNRkAUKcrBTRaP3eIiIgouHFkhIiIiPyKwQgRERH5FYMRIiIi8isGI0RERORXDEbaohjMP8qley2OiYiIqP24mqY1RZsQtfUp82HU+klQYlP82CEiIqLQw2DEAdX3nwKbHoIEYXFeqq3wU4+IiIhCE6dpHFDvzAcgYL1vr0VwwikbIiKidmMw4oBcW952m58KfNATIiKi0MZgpB2ki2f93QUiIqKgx2CkJU0HZNSvwW91zzjVXHRIAsD9aoiIiNqDwYgdBUomlLhUOxkjRqIpbUTplu3DXhEREYUmBiN2KJBxoM8fIYSAIqzeE4AAcMjQC5BVfukfERFRKGEw4sAjh7thhj4XFUi0OF+BTpipfxQz9Y/CYIpUmnbyLYmaDOgu+qG3REREwYt1Rhw4U92A7cjGjoaByJaLkYRKnEWCcQqnKYY7WFqF4f0S/NtRIiKiIMdgpIUYTQRKFt2Gjw//jEfXHgZgnLLZp2TZbf9Lrd6HvSMiIgpNnKaxIykuyql2XWLVXu4JERFR6GMwYkd2j0SkaKMcrKUBJChIwTkMSNf6tF9EREShiMGIHSpZQt4449SMbTl44//NU6+GSnYUrhAREZGzGIw4kNMvBcvvvg5JcRqL813jNViuXoIc1X4/9YyIiCi0MIG1FTn9UjA0PRqFi3NwFgnQTnodA9K1iHuJgQgREZGnMBhpg0qWMFh1FABQl5Hg384QERGFIAYjLjIoAnsNfYwjJSWVGJYV7+8uERERBTUGIy7YUXwOi7b9gO6GO41F0Nasx59i++OpW3rhDn93joiIKEgxGHHSNsMgfPrhKnygXoVUzQXz+bKGRCz4cCoiVYOY1EpEROQGrqZxgkFI2GW4Gq+qlyAZFyzeS8YFvKpegl2Gq2EQTUt9dReBfK3xxb1qiIiIWsVgxAkFypWYrf4IAGBdWsR0PFu9EQXKlT7uGRERUfBjMOIEA2SkShdsAhETWQJSpfMw8HESERG5jL89nZAkVXq0HRERETVzKRhZuHAhBg0ahLi4OCQlJWHChAk4duxYm9e9//77yMzMRFRUFK666ips3brV7Q77Qy+pzKPtiIiIqJlLwciXX36JWbNmYd++fdixYwf0ej1uueUWXLzoOElzz549+N3vfofp06fj0KFDmDBhAiZMmIDCwsJ2d95XVBJwSaihCPvvKwK4JDRQNU3j1Okaze+1/JmIiIhsSUIIB79i2/bLL78gKSkJX375JW666Sa7bSZNmoSLFy9i8+bN5nM33HADrrnmGrz22mtO3ae6uhparRZVVVWIj/dtkbG62irEvJAOABCQICAsIjgFgAQJEoyPse5/SgHAfE3d/5QiJpa7+xIRUfhx9vd3u3JGqqqqAACJiYkO2+zduxejRo2yODdmzBjs3bvX4TUNDQ2orq62eAUC3fjXgNgUi3NS/GXG80REROQWt4MRRVGQm5uLoUOHol+/fg7bVVRUoGvXrhbnunbtioqKCofXLFy4EFqt1vxKS0tzt5seZeg9FvUP7DEf109cByn3CAy9x/qxV0RERMHN7WBk1qxZKCwsxNq1az3ZHwDA3LlzUVVVZX6dPn3a4/dwm6wy/6ikD7Y4JiIiIte5VQ7+4YcfxubNm7F7925069at1bbJyck4c+aMxbkzZ84gOTnZ4TWRkZGIjIx0p2sBq2XuCfNIiIiImrk0MiKEwMMPP4yPPvoIO3fuRI8ePdq8ZvDgwfjiiy8szu3YsQODBw92radEREQUklwaGZk1axbWrFmDjz/+GHFxcea8D61Wi+joaADA1KlTcdlll2HhwoUAgEcffRQ333wzXnzxRdx2221Yu3YtDhw4gDfeeMPDX4WIiIiCkUsjI8uXL0dVVRWGDx+OlJQU82vdunXmNqWlpSgvLzcfDxkyBGvWrMEbb7yB/v37Y8OGDdi4cWOrSa/BxqAI7DX0wceGwSgoqYTBUUESIiIisuHSyIgzJUl27dplc+6uu+7CXXfd5cqtAoemAzLq1wAAijQdbHbh3VZYjryPC3FG/6zxxDuF6BqnwXzDIOSo9vu6t0REREHHrQRWMtpRfA65G4ohQcENcjGSUImzSMD+mkzMQC6WYwnsl4IjIiIiEwYjbjIICc9tP4Fb5ALkqVchVbpgfq9MJGKBfgrm66fgMwdTNnW6RmTN2w4AKFowBjEa/qcgIqLwxF173VSgZOKai19huXoJknHB4r1kXMCr6qW4Wj6Bg6VVfuohERFRcGAw4qaz0CJPvQoAIEuW75mO89Sr8UtNvY97RkREFFw4N9CGGE0EShbdZj6u0xn//+Uos5iasSZLQCrO44pLRwBc6eVeEhERBS+OjLipj1zqVLsrYy+23YiIiCiMMRhxk/XUjCNSbNe2GxEREYUxBiPtoMSmQMB+VGIqyaJ0y/Zhj4iIiIIPg5F20I/MhwTYBCQCEiTTKe7qS0RE1CoGI+1g6D0WmLgKItZyB2IR53hHYiIiIrLEYKS9ssaj/oE95sP6ietw8Q+78Z0hHQeUK/B9wQ4YGhv92EEiIqLAxqW9ntBiKuY/P5ai+/qZ6Ks6bzyxezrO7H4KpwbMgUX2iO4iSqImAwDqdKWARuu7/hIREQUQjox40CFDLww88CS6iPMW57uI8xh44EkcMvTyU8+IiIgCF4MRDzEIIFU2BiGOKrKmyOdhaHvjYyIiorDCYMRDipV0dJUqHdYfkSUgWapEsZLu244REREFOAYjHnJJimx3uzpdIzLmbEHGnC2o0zHplYiIwgODERfFaCLs/hwtGpy63tl2RERE4YKraVyl6QDkV9mczpRLcUYkoAvsT9UoAjiLBGTKpWA4QkRE1IwjIx6ikoAypRMAY+DRkum4XOkEVWt72jQt9y2JmgzouMEeERGFBwYjHnSt6kccGLgYv0idLM6flTrhwMDFuFb1o596RkREFLgYjHhYvxGTEJO731yB9fBNb6HLM9+j34hJ/u4aERFRQGLOiBeoIiLQV1UKAKjLHg1VBB8zERGRIxwZ8RFDYyP3qyEiIrKD/2T3gUPb30bq3vlW+9V0wnFDL+aREBFR2GMw4mWF/1qHgQeeNB60WEnTRZxHF/k8Dhl64Ur/dI2IiCggcJrGiwwC6H5wEYA29qvhlA0REYUxBiNeVKykoyvOt7lfzY/f/st+A91FIF9rfLHuCBERhSgGI17k7H419ZUVTn8m968hIqJQw2DEi5zdhyYqIdnLPSEiIgpcDEa8KFMuxRl0sikPb6IIoFx0xIXO18HgqBEREVGIYzDiRSoJODVgDgDH+9Xk66fh3jVHcePzO7Gj+JxFm5bTMJySISKiUMVgxMv6jZiEfw952Wa/mgp0wgx9LrYr2cbjqnrkbijGNsMgf3STiIjIb1hnxAeuHTMNNdffjiMvjsC7hlEoQQoKlEwoLWJBAWMZkvn6KRgtH/BbX4mIiHyNwYiPqCIiUIsYrFX+CwAgQ8ENchGSUImzSDAHJ+XojAIlE/0dfZDuIkqiJgMA6nSlgEbrmy9ARETkJQxGfOgsEgAAY+QC5KlXIVW6YH6vTCRivn4qtivZ5nZEREThgDkjPpSESoyRC7BcvQTJuGDxXjIuYLl6CcbIBUhCpX86SERE5AccGfGAmFgtkF9l/BlAXW2V3XbZ8lH0lMsB2C8PrwhgvnoluqASzlUoISIiCn4MRnxIJQFdWxn1kCUgmaMiREQUZjhNQ0RERH7FYISIiIj8isGIF8RoIuz+bCJgfxtfAQkNIgIfG25AQUklS8QTEVFYYM6InwhIkCAsjgGB2fqHjVVZ3ylE1zgN5hsGIUe13+Hn1OkakTVvOwCgaMEYu8EPERFRIOPIiB/oxr8GEWu5U2+5SMQM3WxUIRbj5T24QS7CuZp6zNDntl4ivqkIWknUZEB30cs9JyIi8jz+M9oPDL3HwnD5LYh5qQcUAcyKyAcu/Yp56tU2hdAW6Kdgvn4KPuOUDRERhSiOjHiDpoOx7kh+lfFne2QVAOAbpQ+US5VYpl5qtxDaq+qluFo+gYOl9muX2FOna0TGnC3ImLOFu/0SEVHAYzDiZ2ehRZ56FQD7hdAAIE+9Gr/U1Pu4Z0RERL7BYMTPLkcZUqULNoGIiSwBqdJ5XHHpiG87RkRE5CMMRvysj1zqVLsrY5mcSkREoYnBiJ85GhGxJsV29W5HiIiI/ITBSABQYlMcF0JrWkSjdMv2YY+IiIh8h8FIANCPzG8qeWYZkAhIkEynmlbfuKuutgrI1wL5Woe7ChMREfkDg5EAYOg9Fpi4yqYQmohLdnAFERFR6GAwEiiyxqP+gT3mw/qJ63DxD7vxnSEdB5Qr8H3BDhgaG2FQBPYa+uBjw2DH+9ewKisREQURVmANJC2mYv7zYym6r5+JvqrzxhO7p6Ni91P4G+7Bh/pnjeec3L+GiIgokHFkJAAdMvTCwANPoos4b3E+SZzHC+JFjJELzOfO1uja3r+GiIgogDEYCTAGAaTKxiCktYqsMhQAaNr3V2C+for9KRsiIqIAx2AkwBQr6egqVbZZkTVbLjafE5BRjs4u7V9DREQUKBiMBJhLUqRT7ZJQaXPul1q9h3tDRETkfQxGAky0aHCq3eXyT7hBLjJP1wBAl1i10/fhzr5ERBQouJrGR2I0ERY/OwoAMuVSnBEJ6AL7UzVCAJIEzI7YiNnYiDKRiAX6Kfi30gsD0rXOd6hp+S8A1OlKAY0WdbpGZM3bDgAoWjDGos9ERETewt82vqLpAOS3yOnQ2c/vUElAmaETusiVUIRlEquwk5+ajAt4Vb0U/1Z6QSVP8GiXGZwQEZEvcJomAF2r+hEHBi7GL1InyzckNJeHbyJLgNR0DRRD+27MYmlEROQHDEYCVL8RkxCTux/fGdJxTLkMABxspdccoMg/FThoQUREFLg47h7AVBER6Ksqdbq9dPGsF3tDRETkHRwZCSFHa2JY+IyIiIIOR0b8RdMBGfVrAABFmg5tNldiUyDXVsBUc9XiPQFUIRbPf/4jSvftwFO39MIdnuijnRU3REREnubyyMju3bsxbtw4pKamQpIkbNy4sc1rli1bhj59+iA6OhpXXnklVq1a5U5fw5p+ZD4AQFhljoimFTcdpVq8p3kO7zc8iE8/XMW9aoiIKGi4HIxcvHgR/fv3x7Jly5xqv3z5csydOxf5+fn47rvvMH/+fMyaNQuffPKJy50NJTGaCJQsug0li25zasmsofdYYOIqiNjkVtsZl/ouwS7D1TAIRymvREREgcPlaZqxY8di7NixTrdfvXo1HnzwQUyaNAkA0LNnT+zfvx/PP/88xo0b5+rtw1vWeNR3uxEnF9+EVPlXJKDW7lJfRQCz1RtRoFyJ/v7pKRERkdO8nsDa0NCAqKgoi3PR0dEoKCiAXm9/L5WGhgZUV1dbvKiJrMKviENHyTYQMTdp2kzPwPxkIiIKAl7/bTVmzBi8+eabOHjwIIQQOHDgAN58803o9XqcO3fO7jULFy6EVqs1v9LS0rzdzaCSJFV6tB0REZE/eT0YefbZZzF27FjccMMNUKvVuOOOOzBt2jTjzWX7t587dy6qqqrMr9OnT3u7m0Gll1Tm0XZERET+5PVgJDo6GitWrEBdXR1KSkpQWlqKjIwMxMXFoUuXLnaviYyMRHx8vMWLmqkk4JJQw1FJEUUANSIam5UbUFBSCYMiYFAE9hr64GPDYPM5t+guAvla40t3kbv/EhFRu/mszoharUa3bt0AAGvXrsXtt9/ucGSE2hYt6SEgQYGwiCgVAJCAN/S34pRIwdk16/F4ZD8okFGpf9bY6J1CdI3TYL5hEHJU+9vXEdYiISKidnI5GKmtrcXx48fNxydPnsThw4eRmJiI9PR0zJ07Fz///LO5lsj333+PgoICXH/99fj111/x0ksvobCwEG+//bbnvkWY0o1/Deqd+UBtuflcpYgFADyh/sB8rkxJxAL93aiU45GESpxFAvbXZGIGcrEcS3CTC/es0zUipsXPRERE7eVyMHLgwAGMGDHCfPz4448DAKZNm4aVK1eivLwcpaXN+6kYDAa8+OKLOHbsGNRqNUaMGIE9e/YgIyOj/b0Pc4beY2G4/BbEvNQDigD+qZqI+w3rbdol4wKWq1+2WH1TJhKxQD8F8/VT8BlLyBMRkR+5HIwMHz4cQjj+5bVy5UqL4z59+uDQoUMudywcxcRqgfwq488A6mqr2r5IVgEAvlEyMQ6fG0/ZqT1i/Z/MWBxtKWboc3GwtArD+yW0s/fN6nSNyJq3HQBQtGCMU0XdiIgofDFpI0QYICNVumATiJjYK44GAHnq1filpt6znWnKIymJmmxMeCUiImoFg5EQ4U5NEVNxtCsuHfF8h1rgihsiImoNg5EQ0Z6aIlfGcvSCiIj8h8FIiFBJwKWoJIe1R1ojxXb1fIeIiIicxGAkhMi3/BmSJBlrjbTgKD4xJbUq3bK92S3mkBARUasYjIQQQ++xkCauAmJTbN6zDkgEmpNa5Z8KAMXg9f4RERHZw2Ak1GSNR/0De8yH9RPXQbrrbQjrAEVq/k8ftX4SxJJ+UH3/qa96aVNWnoiIwhcLQASwlvU5XKrV0VR7BACU9MFArBb1aTch5qUeAIyjIkIoaLnaV1SXQbPpoXb2mIiIyHUcGQkXTQHKJaGGELb/4WUAQgjUCQ0OKJfj+4Id0On0ntlcj4iIqBUcGQky1qMlrtTtMAjjBnuOyBIQAx0GSseB3dNR/mUiVhqmYruSDbxTiBTtccwZnYE72vUNjKz3uInReOBDiYgoKHFkJJBpOhjLw+dXGX9upx9Fqkvtu+IClquXYIxcAACoqKpH7oZibDMMandfiIiITBiMhJGzIsGl9qaS8YvUb2KwXAjTouH5+ikwCAd159vDiaRWVnMlIgo9DEbCiAoKykSiS4XRZAnoKNXiPc1z+CpyNm6RC1COzihQMr3XUSIiCisMRsJItnwML+v/HwC4Vak1ucW0zVkkeLZz9tgbKWEBNSKikMNgJIyoJIHhqv9gpj4XFUh0+fqWO/0m4VcP9w4W0y6cgiEiCh9cTRNmclT70XDH67hz22B0rzuCrriAeerV6Igac7DRGlkCUnEeXeXzaPB+d9tUV1uFmBfSjT//TyliYrV+7hEREbmKIyPBxgMrbEZndsZns2/Ao6oPMFJ1CF9d8UcArk3dHBJXeL32CEdKiIjCA0dGwpRKljBYdRQAUDdhHf79dTxS985HV5x36voXGidh3zuF6BqnwXzDIOSo9nuzu0REFMI4MkIAgGvHTENM7n58Z0iDTkQ43OlXEcA5EY+uuIAb5CKcq6nHDH0ua48QEZHbODJCZqqICPRVnQYACEjGnX1bhCWKACQAnaVqLNW8CgAoE4lYoJ+C+fop+CxQysXrLgLPNRV4e7rMIwXjiIjIezgyEuTc3kyvDbrxr0HEJlucs5ffmowLeFW9FNMjtmLnZ5/A0Oj/3A7mmhARBReOjAQ7U0Krhxl6j4Xh8lsQ81IPHFAuRw/pjN0VN6bj+yK2AYXbcKYwDyXXzYFi6IOzSIC2pBID0rnChYiIHGMwQo417fQbLXToJNc4dUkXcR5dDj6JGYZc8wZ7THIlIqLWcJqG2tRHLnW6rWmk5K/qFbhD/opJrkRE1CaOjFCbnCmGZt2+MyyTXOfrp2COfjpeOlmJm/vGe6GXjtXpGpE1bzsAoGjBGI/m1hARUftxZIScosSmQNhNYW2bcU+bpbheLsa97xbixud3YkfxOQ/3sBX29rOx3vfGiR2DiYjIOxiMkFP0I/MhAW4FJC33tJGhoKKqHrkbijltQ0REABiMhCZNB2TUr0FG/RpzjY32LgE29B4LTFxls9zXWbIEpErncY9qGyQoAID5+ikwCPdGW3yhTteIjDlbkDFni0tLhN29jogoXDEYIedljUf9A3sAACVKVwjh2n42ADBP/Q6+ipyNMfI+dJfP4keRArl0L6AYvNBhx1iLhIgocDCTj1zTtNw3Qz6DggGL0f3gIqf3szEx5pC8DMk0KLJ+EhCfCtXwPA931vO4SzARkedxZITc1m/EpKb9bNJxQLkcOk1Hp3JK7K7OqS6HZtNDnu9kAOC0DRFR6zgyEi6sK7XqPFO11bifjbEOSUPO65A2PQQByWJPG3skm4DEsr1cuhdKt2yP9NFtTatwAKBOVwpoOApCROQNHBkhj2lvkmvL+CRq/SREvTHE9h6KwF5DH3xsGIyCkkoY/Lw5H0c9iIjajyMj5FlZ41Hf7UbEvNSj3R8l1VZYHG8rLEfex4U4o3/WeCIQSs3bGz2xOcddg4mIWsNgJATFaCJQsug2/3WgKckVgFNTNo60vO7br7Zi1q5YGKwG887W6DADuViOJbjJvd4CMI5wxLT8WdOODyMiIpdwmoa8zjqpVbgRm9y47yH8X+RsjJELrD7b+H/n66f4fcqGiIjcw2CEvEo3/jW7OSTuBCTGJcFL7AQkMsrRGQdLPZOUS0REvsVpGvIqQ++xMFx+C2Je6oEDyhV4oXESOqIGz6pXIxUXXPosWTIWWfuregWi9DqcQSIKlEwoTTH1L7V6b3yF9rOTV8LN+4iImvFvwHDVVDIeAIq8nWDZlEOiFxHYp2QBALY3DEK2XIyuuIB56tXoiBqndge2vyPwVGxXstElVu21r0BERN7DYIR8JlsuRtc4Dc7W6KBANgcm9XoNlquXQBEOCqK1wjR1s84wHNnSQ4DSzQs9JyIib2LOCPmMShJ4ekxPAJY1RbYr2Zihz0UFEl3+TFkyvn4XsQsx7/8WYkk/qL7/1EM99h3WKyGicMZgJEyZlv+WLLrNp/kKozM7Y/nd1yEpznLt7HYlGzc2vIzf6Z7GryLW5Q34TER1GTSbHsQxw2UBUxiNiIhax2CEHGoZpHgyYMnpl4LPHxmE99R/xsOqD83nFcjYq/TDHP19xmM3YggZxpU6cXI9HtPPwj3vFGLUP/Zjm2GQh3rvJU1JriVRkwHdRfM55GuNL9M5IqIQxGCEzGJitcb9a/KrvL4brUqWMFh1FI9FfICucRqPTdsAxmmbVOk8suViAE2F0fS5gRWQ2As+iIjCFIMRcsy0uV5+lddKmreWR3Jjw8v4re4ZPKqbifMizuWRkiRUAmBhNCKiQMdghFzijakbR3kkphU3Hys34mn9dOM5F2KJy+WfcINchAg04nq5GIPk7/Hjgc8BxeCRfvuas0muTIYlomDDpb0UEHL6pWBoejQKF+dghzIAKwy3QgLMu9NsV7IxU5+LeepVbRZLEwKQJGB2xEbMxkYYhASV1PRJuwEcfhaqm581t5dL90Lplu2V70VERG3jyAgFDFMeyTz1O1j6m0ybkZJ/x92E/bdvh0EA7+mHQwjbkRJ7ZeZl6436qsug+WSG+TBq/SREvTHEU1/DY1qOatTpGv2eZ8IRFyLyFo6MUEAandkZI3t3QuHiHJxFArSTXsewrHQ01FVDtQXoKF/EDAcjJZJV4TTrY3uk2gqbcwZFYK+hj/H+JZUYlhXfnq/kNyw9T0SBjn8rUcAyjZQAQF1GAlQtyrPmqPaj4Y7Xcee2wehedwRDpELMVm90KvAAAOtmEoR5ekc+9TW2KVcjb9NRnNE3Tee8U4gU7XHMGZ2BO9r4bOsgZkC6F1cm6S4Cz6Uaf366zPlEY3evIyLyAgYjFLSaR0/yoZbaP21gCmSi3v8drhaJuEY/FdvRnEtSUVWP3A3FiFQPQo5qv93P2FZYjryPCy2CmK5xGsw3OL6mPep0jYhp+bMGthvzATYb9dlgcEJEfsScEQpqptGTgfIPHv1c0543Y+QC8zlT5sl8/RQYhO0QzI7ic5jxzrc4U6OzOB+QdU6IiAIIgxFyjRO1R7xVubUtSmwKbCdg3GOaEfqregXukL/CDXIRZCiQoKC7fBafGgbh+4IdMDQaR2QMQsJz209Yp8oCMAYxAgJz9NOx72Rglqe3SZYFE1aJyHc4TUMhQz8yH5GbHoKABKlFWGDKBXF07IgsAZ1RjaWaVwEAF0QsACBRqjU22F2AM7s74cvGbBSLdJxp0Dn6KAAyKhGPe991PvfEm6ynd4iI/IkjI9Q+mg7IqF+DjPo1fs8zMPQeC0xcBRGbbHFesfPH3N4S4LZ0RC06otbiXBdxHreqCvCD6Ob055hyTzhtQ0RkxJERahfT7r8BI2s86rvdiJOLb8brhvE4iwQcUHpjoPw9klCJs0hAR9TgWfXqNounWbM3miJLxloni9Rvolofg2+ULABAtlxsvl+BkmkREAkYJ5Pm66dgtHygHV/Wx6wTY5kIS0QewmCEPM+UV9IkBs0FunySQyKrcByXYZPSXMhsX1OQYLK9YRCul4vwqvplaFELuR2pJrJkHDV5T/Oc7VQOgDKRiPn6qdiuNK/MEQDK0RkFSib6u39rv6urrULMC+nGn/+nlDVMiMgtnKahkGTaJM8RBTL2Kv0wR3+f8dhDOaX2pnLsrcwxOYsEz9yYiCiIMRihkJQtFyMF51tdWxOHWmxXsjFDn4sKJHrkvpJkO51jGnXJU6+GDMXivR+Uy1BQEpgrbDxWfl53EcjXGl8ufA5X8xCFDwYjFJJUkkCeehUAe9VWja+Jqt0AjJvw3djwMn6rewaP6mbivIjz2EiJiSwBqdJ53KPaZhGQvGL4b9zzTiFG/WN/WCW0MtAgopYYjFDIylHtxxI7G+4la6Ow5DeZGCUfNJ9TIGOfkoWPlRvxtH668ZydTfjcWYXT0jz1O/gqcrbNlE2oFEazV6/EKRw9IQprzDYj77NKaDWdy6hfAwAoMq24sD7ngZ1pW9twL/IT41ROBTpZFCszTd3kWW3Cp4tMwCVdIxKsckJcZcohmaHPNSe1GlfYCMzXT8FngThl0w52S9a7gyt1iEIWR0Yo5JlKxt+h2ovsFhvutTaV85mSjWENL6PAcCUAoH7iOkTOOQFV7n/wnSENOhFht9qqMxzlkAjIKEdnHCytcnAluYOjJ0SBj8EIhbXWpnJe+k0WslXHAABK+mBAVkEVEYG+qtPQSI0wjmVYhjECcCpIMeWQZMvFxmMouEEuwnh5D3BqL6AYPPDtwoBVki0DD6LgxGka8gt3i6VZ73vjiV84rU3lYLPj63TjX4N6Zz6k2nLzOSk6EUIIoP5Xp+6dIxdgkFSMyRE7kSI1TQntB8Sx+VANz2vP1wo5HpvuIaKA4/LIyO7duzFu3DikpqZCkiRs3LixzWveffdd9O/fHzExMUhJScG9996L8+fPu9NfCic+LDXvaCqnNYbeY1H/wB7zcf3EdcCTx3Fp5rdO3/eeiM/whHoDkq2qwYrqMmg2PYhjhsvwsWFw4C7/dZenlg27cy97ybJuJtASkWe4HIxcvHgR/fv3x7Jly5xq//XXX2Pq1KmYPn06vvvuO7z//vsoKCjA/fff73JniQKOrDL/aJrKaXnOehrHfN4qrrCpTdLUppNcg52Ga/HamvUYvXQvlur/n8vBiUER2GvoE9BBjdurcAL8XkTkHJenacaOHYuxY8c63X7v3r3IyMjA7NmzAQA9evTAgw8+iOeff97VWxMFJZtdhOHersFljYmYL6bi7/q7gHcc7/5rCj7OIgE/7D6F9w+dwRn9s8Y3W7ku4FjvhePONfb2z/EkrvAh8givJ7AOHjwYp0+fxtatWyGEwJkzZ7BhwwbceuutDq9paGhAdXW1xYsoGOnGv2azi7C72+CYlgQ/E7EKN8hFOFtVZ7P777bCcoz6x378Tv8sHtU/gld2n8aZGp3F53DXYFt2R0usp27sTOVwlIXIM7wejAwdOhTvvvsuJk2aBI1Gg+TkZGi12laneRYuXAitVmt+paWlebubFABMSa0li24LmQ3XrPNK2kOWjK/7IrZhreYv+L+m4mnz9VNgEBJ2FJ/DjHe+tQk+rJnGaEzXERH5m9eDkaKiIjz66KOYN28eDh48iG3btqGkpAQPPfSQw2vmzp2Lqqoq8+v06dPe7iYFILvBiamAWn5V8AyJt8gh8aRkXMCr6iW4Wj6BF/V3In/rcadrn0hQ0F0+ix9FCuRSLiW2x1OjHnW1VeYRlbpa1pAhssfr//xcuHAhhg4diieffBIAcPXVV6NDhw4YNmwY/vKXvyAlJcXmmsjISERGRnq7a0Q+p8SmQK6tgHPVSFonS8aS9X9Vr8Cf9Xfjyvr/oACZUKz+jSFDQbZcjCRU4iwS0BE1eFa9GqmmpcTrJ0GJtf3fIVny1NLiOl0jsuZtBwAULRgTMqOARO3h9f8V1NXVISLC8jYqlfFfiqK9G30QBRn9yHxEbnrINqlVOJfUas0myVUkYoH+blQiHkmoRHepwrKGCezvryPVVrR5r5aJsdqSSgzLine9wyHOqUDD10m2REHA5WCktrYWx48fNx+fPHkShw8fRmJiItLT0zF37lz8/PPPWLXKWGZ73LhxuP/++7F8+XKMGTMG5eXlyM3NRXZ2NlJTUz33TYiaeKMwmqcYeo8FJq6C2PqURbG0M1Ii3tWPQKnoinnq1eiIGjhR6sSGMcn1ZYvAxm7wYfXZEgSEAC5Bg4KT5zGsb3zorMoJAhwtoXDn8p/4AwcOYMSIEebjxx9/HAAwbdo0rFy5EuXl5SgtbV6Gd88996CmpgavvPIKnnjiCSQkJGDkyJFc2kueYW/DvUCXNR713W5EzEs9ABiLpXVIHYQhL96Gy1GG8iHzkbj3CSgQLid1yVLbNUwcjcBIEhADHXavfRGPa26HAhmVpuBjt23elmlVTqR6EHJU+13saYjiqAeRW1wORoYPH97q9MrKlSttzj3yyCN45JFHXL0VUeiyKpamAjBYdRQAUDd0HaS0LhBbnwJajJ44y53pnpbmqd/BfcpWzNdPxXZkO2xn3GnYuCpntHygfTclS6xfQmGGY4EU1Jza48be6In1uUArAW41emIM/y3zTLzJVNNkhj4X2xVjQGKdCFugGJNly9EZBUom+vukZ+HB2WTZutoqxLyQbvz5f0oRE8uRGApODEaIAlWL0RPd+NdtNuXz6q2bVurkqVdjR8NAjJYPIE+9qnkFDiyTZdVSI+TSvVBSBzR/RuleIHO0T/ob9JyoNutM4OFu7glzVsjf+CeOwpL1iEpd63XC/M7QeywMl99iHilp0HSEuuFXt5JcAedW78gSkIrz+If6ZdwqF9iMydgky66fBCE1Z7lErZ8ExKdCdfOzzZ/JACUw2QRDEZwmIp/yetEzomBlvSrH71oWT8tZBEmSoFg1cWYSx9V98m5TFUBqqv5q0R17wYyw6lF1GTSfzDAfRq2fBCzpB9X3n7rWCXKOL3dDJvIgBiNEQcjQeyykiasAO8XKrGMN6+MKdMJM/Wws0N/d7n7YLhF2QnU5NJscV2Amz6rTNSJjzhZkzNmCOl2j5yrC2tmrx5n7E9nDYIRCj59LxvtsRCVrvMW+N/UT10G6620IqwCl5UZ99RPX4T+/2Y1vO9yElYYclIlEl0dKXGUboDh3Q1Odk48Ng1FQUglDY/MvMlMJe5s23v4y1DonAxQiawEw9kwUAIKxXglgs0QYsVrUp91kUcNE6ZZtPlbSByMnVouh3WNRuDgHPymdkaK6AAW+/ZeJzaqgFnvjyKV7sc3QD3mbjpqLrI1ZswJ9NKsR19Qmav0kXIpOxl90d2Ndi0JsXeM0mG9g3RMLTiTHeordVUBBUHuFCbz+xydO5IhphCXYWAcodqhkyVzXpKFppY47NU3a64ByBaK+3oCsI4vN56LWT8LVIhHXNNU5GSMXYLl6SXNhkyaRdRVYiBeQpLoTp0QKziIB+2syMQO5WI4luMnn3ybEWdU+8dRePabrLYIBNLSZQMsAIrTwvx4RnKxX4uA6ez+72safrFfqOMvd/XRa+tLQH4/tedxmLsdU52Smfjbmqd8BYD+BVgjgCfUH5nNlIhHz9VMwRz8dL52sxM19uX+Op1gHHwHJulgcwFVBQSLw/mYkCmAxsVrzaEkMEDrz4i1HU1zYWdg6IGnr2HwPAVQgEZMjdhpvb92dpjonf1GvRGep2uH9rT/bGMQsxQx9Lu59Nx4p2uN4alR3JLfY4O+abvH4Nzf8a5UnRz2cYjWVUwf49v7kd0xgJfI0PyfQtpd+ZD4AQFgNV7QMTb664Q3M1D+KciRatFHs/JVivXuEKcf0vcYRSJEuOKyVIktoNRBxdA0ALFK/icFyIc5U1eGJD45iqeFOfGEYgNfWrEf287stjm9atAM7jp5p/oym5FjyDo+t5nH3/s6s7nEiEZerhDyLIyMU8tydgglXjnYWFnEpkGqMx9fdmIPabjfhzo1D0L3uiLlE/AGlNwbK35uPO6IGz6pXIxXNlVsr0AkL9L9HluydZEpZAjqiFu9pnsMFEQsASJRqm7+fkKCSmiOkCw2xwGaYp4qi1k+CiE+Fanhe2zezSrxVujney4d8z+cjPOQ2BiNE7RGsSa5tsbOzcMtVOQCQ0y8FQ9OjUbg4HzuUAdhkGAIJwD4ly+KjtjcMstjTpiNqME+92qK0vLd0RK3NOdlq+sleG1Fd1nYtlKJNiNr6lPkwav0kKHbqvpBnWJfDB+DeKh17q3tszvFXo69xmobIDwI9qRWA7aqclhVgm5hW5cxTv4Olv8lEUpztPz0VyNinZGGTMgRa1GKZeimS4f1ABDDmlNgUZrNzbH1OBiCEQIMw/reRS/cCjc17BkTsWQKsn2qzV5BUW+GhnpM7Wk6XtGfqxFOf0577h9sUUID+LUhEwWZ0ZmeM7N0JhYtzcBYJ+GHw37Bs92kAxnwTGQry1KsAOCgl7yJPrOZpjSwBkTD+Ioiy2ndHs+cl65XGAIz1U8z9aprCMRVmMyXMDkgPvDobRP7GYITI20J1KseOlvVL6m7qjn7pXZD3cSHO1OiQLRd7bGpmS2M2xqoKIIRnAhunKIpF9OHotqYASf6pANuU/sbvz8Js4cl6qbGmA+ujOMCnQERe05xXkgO15Lnh5tXKLfhEGYJ56lUWybHe5OoozI4DRZhdZLBZIH22RsfCbCHKOq/Fm4FGqAU1wd17olARrOXondBytMQZ50UcOqLG7oiHsT5JJxQomeiqjcH+Ub/HHVsGAbAp0Op37xYbzNNTLRN4jSuOirHLcDWSCnag/00T/N1V8iU3y+NbBx/BUGbfFQxGiAKUdZKru4lsgZQs21pBNVOg8Wf93XhVsxQCksUeNqbjn5ROWHH31RiWlY6Gumpgi6mF1BSQtLzG9E4ze4XZAM/lnygCqEIsJCgYK39jXNrcYnrKYmnx7rdwZvdTODVgDrgoOLB5LRiwM5XjqJ1H9hgK0Cq1XE1DRD7juKCaBAnAT0onTJo8Hcpdqyx2GwYAEZcMSQKyVd8jOyMBKquhE93412yukaITIaISLDshW/619yti8Sti3f1Kln1symHpKBnrnLyqXooUq2kk66XFXcR5DDzwJDY3Zre6QzH5WVMwUBI12enKy86syrHbxo17Ofrstlbl+HvlkAmDESJ/CPIqre4yF1RrK9DoewfqH9hjfr9+4jrU37/H+uNsPtv6Gjx5HPUzD1mcu/To9xbHBb/5Brep38Jvdc/gUd1MnBdx5iqxbbGuLmtPW0uLTTHVdaof8Zh+Ft5bswJ1SwaZ349aPwlY0g9S8WZ8Z0jHAeUKfF+wwyJgMTGt3DEHNc5+ESI/4zQNUZAIpOmWdnGioBoAp3YftmGvNkorn6OkD0ZOrBZDu8eicHE+KkRHPCc/gMXKi1CsVupYT++Yfs//Xf/f+EPEZ0hAbZuBh8NuS0AqzmOWaiMei9hgkwAjqssQtXkm+pq+yu7pxumd654yT+98+9U2PHlAi/IWK3dStMcxZ3QG7nCuGxTkgjmvJIj/RiMKHWFXst6dQMOLWibZyre9gpkfGmxW6iiQoYJiPq5AJ8zXT0EVYvG49KFH+vFQxCeQYGc0BU2jMC3OdxHnkXTwj+ZzN+57AB+IRMyXp2K7YgxRKqrqkbuhGJHq5qXE1nVPhmXFsxZKqPBUXokfMBghChZO1isJmREUPxmd2RmRkx/CnRsHW+y7863ojeuk5n13Tsf2x//c3hNZm8Z77N4dpAaH79mb3rGeJjLuWrwEM/S52K5km7NTntbfi0tCjR93n8IH35ahu+FO4/dYsx6PR/aDAhmVrtZCsdqXB5mj7QY6FFis9+sJFPybiihQhVGxtEDTct+ds0iAdtLruKZbPI6+uNB8bFrNE7Pl53bfz91qsvYCFEUYdy2u1sfgGyULCmRcgBaPNT6MMV99gA/Uq5CqaR7xKVMSsUB/NyrleHOgtb8ms9VaKKrvP4V6Z775OGr9JFyKTsZfdHdjndU00VOjuiOZoy7UBgYjREHM3VGQmFitOdCJAfyylXugs6gmm5EAABbH1qt5rJciN59ve2mxJ8vat9y1uEwkYr7eOG0zRi7AcvUSm/bG0ZSXLfpQJjpibeNI7DFk4aZTe4A+t1hcY9xE0PK7RtZVYCFeQKWca54mKq+qx2MfHAPQ+qgLp4mIwQhRMOPoSUCxrY1iS5Esc0+8ub+OadpmSaMxydZePoq96Z4U/IrH1R8YD95/GyI+FarheS1aWC/Obh6VyVOvxo6GgVCaFmtaF32zHnXZVljOkvnEYIQo5NgLUBiweJ1u/GtQ78y32MlXir8MDTc/g8hPZgAwrhzSJ1+H0r//F6IlHXrKzu3y6+7oiSnQeFzdeoJtW6uARHUZNJsebH6/lful4jxyIzZgj9IPHVFjU/StTCRigX4K5uun4Kmic3j8w2KboO1cTT1WyrcgRTqPK5vyUbyJIzP+x2CEiMgDDL3HwnD5LRZLlqMyR8NQVwt8YmyjpA+GCkBflesrHdyd3vHEyIsM29U8rZkdsRGzsdFuHZZkXMCr6qVYYcjBpi3nIKE3RIuSV2PkAuSpVzUHMOsn2RmZccAqqVZJHWBxbC+o4chMYGAwQkR2c088VY4+rNirc9KKtsrjVyEWM/WzkYCLxhEGq2quns43aY2793FU5O2+iG24D9tQFtmUQIt4jJIOYnrEpzZPwzgy85DFOZuVO9IuRH36R/P7UesnQUiyxTGsgpodxeeQu8F2ZIabGfoegxEiIj/Rj8xH5KaHbHJNTAXV5ujvw16lHwBge8Mg/Jd8EPfJW/BP9e+g1lXaDVACiTMBjL0EWpt8FACKEKiHBtGSDjuKz+G57SfMoxlj1qzAzZolth8uFMvj6nJoNj2EY4bLUCS6Y8HW43bzeiQoyJaLuJmhDzEYISKnEmGtV+BQ+5nL4299yiLXxFRQzbQqRQIgIONO1Ve4XvU9Km67A7kbivFZwyBky0V4Vf0ytKi1u9NxW/y927G9BFpH7aKhwz7DlcjdUAwJCm6Qi9EVFzBPvdq8L1BLtt9LQACIk+vxRMNMKHZK6ttMEznYzJA1VTyLwQgR2efMSh1NB2TUrwEAFDXtsWN9TG2wKo//1Q1vGMu6NzTnPyRrozBndAZyNhtzGEZndsbyu69D3seF2FvTD3P092G5egkUuLfhmL8DElemgEpEMh5WfYjJETuRIrk+KiQDSJWak2wLlEzzyh9Hy5+7iPPocuBJHFJ64VrVjzYjMyy9334MRojIY+yWtbcTsJCVFrkl192Yg8+GCBQuzrEpsIbNzZc0F2Yztisa8hKyjiwGWoywtKX+rvcQBR2UT+dYjMz4MhfFVb+N+BKAc6MprTEl2RpX99yNKsRikfpNh8ufFQGkyOextXGg3TwTe6X3yXkMRojIbT7dUyeMghrrgmvWBdbsthu6DvWD/595hEU35HFo9vy9adSjZe2T5vwUpftQIFaL+rSbEPNSD/xb6YnPDdcZRx1a5KI4Ck4UAVQgERKArrhgd5rIU6Mu1n3wVLBkL2fFHlkCklGJ95XhdvNMTN9zvn4KRssHPNO5MMJghIiCFlf8WGkxwtI4JBeabtfa5KOIuGRINeV2r+svn0DJ+A347213mffl6S5V4LGIDTY5GaYk2/n6qQBgnCayk7fhCd4cqXE2Z8UkDvUO3xMAytEZf2+8EwNLKjGgW2zzfRzs32Ovpokz+SihlrPCYISIKFRZ5aPUT1wHpVu2+die0ZmdMbJ3JxQuzkeF6IgFmlz8UN/NZhdj6yTbGfpcY+Jny9U98ZdB11T0rUTpiu7SGQhYBizOBBruBiLeqMVyFglttnnF8N8Ys2YF+mhWI67pnGn/nr82TEZJi40KT8VchbwWNU3s1T2xzkcJxZwVBiNE5FU+ncohW9a1T5zQcvpHvvUV5G5oxI6GgRjUoqx7y8RPAPhMycbnDQPxnvrPyFYdsyn6liGfQcGAxeh+cBG64nzzzSTbqRwFAASwS7kaI1X/cfurK7Asvd8eximpTihQMgHYlrk3nc+Wi5vrpVh9MeP+PS9B0jSfK9MnYoHBOLrU4KDuSUVVPR7dUIwS1f9DFTpgxYZiyE2ricxl9qsynctZsS4M1y3bcVsfYjBCRH5nb9mwzRSMvQutV/zo2i577+7mguGq5cqdfTVZ5vOycb2xmWnFT/bmYwDsF33rN2ISDMPuxHcvDsclKRIRw+fgqiQVpG1zLaaSEJeCHypllIhkAM4HI4oALiAef9bfjTNIxAGlNwbK32OIVIjZ6o1ufPvmzwWA9xqH43Z5H7pLFTareS4I45RMolRrPmdv/x7rKSFjRdol+B/9g/iXg7onAsbgZ6/oiyRU4hE7q4mMibhTLXJWnCkMp3RIbu6fqUptG8X6vIH/SySioODMCAtroXiH9cod7aTXcU23eBx9cWyrK37sUUVEmMvh12WPhipWi7r0m22mkq58qQcuGOJQJhKR7CA5tiVTwPAn/b3mqSMA2KdkoQCZ+E3Ebqc+x55KGAONJ0ybB8I2qOiIWjjD0UqdZ9Xvory+E76BMeBrOepib48fR0HNDH0u/t54J7D7FN4/dKbNwnDSxeb9kaKaSu9LOc8DWeOd+j6ewmCEiAKTt3YktvO5HC1pm/UKHwBOrfhxioOppGz5GJ6OeAILG1+wSY61zgexzmFpSYGM+fqpeE2zxGZKyNFqHyGAS1DjtcZxyI2w3Wiwrc0FXSFLxmDmPc1zdkdY7CXYOgpq8tSrcWPDUii7T5unclwpDCeqy4D1UyFNXOXTgIT/qyOi0ObMkmB3Ax/rz9ZddLeXZIdKEhgydjJmfthok0BbjkSs1Y/Af6kOoWbIHEzbFWORw9JMYLrqUwz73VNQpOsgffpHy52VYRuQCEiQJIFIocfv1V8CXlolZI+jERZngh3TrsnZcjG0qLWsJOskc+n9T55EdOZtPpuyYTBCRD7ny5EIJtAGt9GZnRE5+SHcuXGwebnxWSSgtMNVmKe8gP7yCdTdmINl3eqMq1BqdOZrk+M1yL/0PHJU+1GX8RxUsXegLm2Y5c7Koh7Cuuhb0/JnlQRjsq0PC8DZCzpcHXV5QPUJRsj/tpt/4gxZAqIvVcBQ8jVUPX2zVSCDESLyPW9NwXgR81H8pzlnJd+cozIgXYu4l/bbaZPjsA0A2ymhFkXfAOeWP7vLV5VtR6r+DaD9MdSPJ35EbwYjRER+FEYVX4OBvZwVd9rY5cbyZ3d5MyDx9GefFQno7bmPaxWDESIiO5xavWPTxnb6iVVig5sSmwK5tgJwYtLDpl5K0yUrDDn4XBnYvCoGrm/w1+a9PRiImGqqqDKGeuYDneDOBo9ERERhQT8yH4AxqbUlm9AkOhGI6mhx6qzUCTP0ufhL41TsU7JwOO5m7L99OwBAJ1TmYMUZ1k2tj9tTpbYlU59eVk9Hdq8u7n2oGxiMEBEROWDoPRaYuAoiNtnifMvj+onrgCeP49LMby3Odcjdj3tU27FU/Q+svLsfvvrjSIzu0xUAoJEMkCTJpj6sgDNjMLDpj7MUAZwT8XhUNxMv6n+DCiRavF+BTpipz8XwCfe2b7m2izhNQ0TkKV5MzGUCrR+1scePudqsVe6JCq3XYtGNfw3qnflAy+qz0YnG4Yr6X5vPtdjjx979nWWvMNwywwSLAmunY/vj2buuQk6/FJc+u70YjBAR+RITY4OTF5JcDb3HwnD5LZZLjTNHo662yuacaY8fe/d3Nq+lAp2wwKownAIZ+5Q+mK76FP/9uwcwLCvdpyMiJgxGiIh8yKm6J24GLM7Ub2G12QBjHeTYGWFpq/CYfmQ+Ijc9ZCzW1iIgaZlQWz9xHf5j6IdDm44CrdVi8UMgAjAYISLyLiembpzaFDAIMNDxD3Ney9anrIq3pUCqMR4r6YORE6vF0O6xbddi8QP+aSEi8jeb3YedXP5rM4LSYNPEJkDwVF4Lp5sCSxt5LSZu12LxMgYjREREocCHxds8jUt7iYjCnPXoiSenW5jHQs7gf3UiolDRjt2HLaeJuPsw+RaDESKiUBZgmxIGQ70UjtT4Hp8yEVGAcWr5rwvtXGYvgHE3YdWNYMjXwQCDD//jUyciosDSjukm64CJK36CA4MRIiJqk/UojN3dhz2Vs+LkNb4MNLj7sncxGCEiIpd5bYrIg+z10akpmQDLswkHDEaIiCjoBEMwRM5jMEJERIGPFV9DGoMRIiIKHw5WCnlixQ/zStzHYISIiAKeM9MynLoJXgxGiIiIvMHO1JKnRk9CrTaKy3vT7N69G+PGjUNqaiokScLGjRtbbX/PPfdAkiSbV9++fd3tMxERkW+ZpnLyq4w/Wx/bYRqpKVl0W3PA4OR19n721PfIqF9jDJICKPfG5WDk4sWL6N+/P5YtW+ZU+6VLl6K8vNz8On36NBITE3HXXXe53FkiIqKA5USgYc3ZwCPUNxx0ubdjx47F2LFjnW6v1Wqh1WrNxxs3bsSvv/6KP/zhD67emoiIKDzZS7JtI+k2GPYBMvF56PTWW29h1KhR6N69u8M2DQ0NaGhoMB9XV1f7omtERES+5ecCa4ESsPg0GCkrK8Onn36KNWvWtNpu4cKFmD9/vo96RUREFEC8GKAE6oojl3NG2uPtt99GQkICJkyY0Gq7uXPnoqqqyvw6ffq0bzpIREQUrNzIWQkUPhsZEUJgxYoVmDJlCjQaTattIyMjERkZ6aOeERERkT/5bGTkyy+/xPHjxzF9+nRf3ZKIiIiCgMsjI7W1tTh+/Lj5+OTJkzh8+DASExORnp6OuXPn4ueff8aqVassrnvrrbdw/fXXo1+/fu3vNREREYUMl4ORAwcOYMSIEebjxx9/HAAwbdo0rFy5EuXl5SgtLbW4pqqqCh988AGWLl3azu4SERFRqJGEEMLfnWhLdXU1tFotqqqqEB8f7+/uEBERkROc/f3t09U0RERERNYYjBAREZFfMRghIiIiv2IwQkRERH7FYISIiIj8isEIERER+RWDESIiIvIrBiNERETkVwxGiIiIyK8YjBAREZFfubw3jT+YKtZXV1f7uSdERETkLNPv7bZ2ngmKYKSmpgYAkJaW5ueeEBERkatqamqg1Wodvh8UG+UpioKysjLExcVBkiS3P6e6uhppaWk4ffo0N9zzMj5r3+Gz9h0+a9/hs/Ydbz5rIQRqamqQmpoKWXacGRIUIyOyLKNbt24e+7z4+Hj+4fYRPmvf4bP2HT5r3+Gz9h1vPevWRkRMmMBKREREfsVghIiIiPwqrIKRyMhI5OXlITIy0t9dCXl81r7DZ+07fNa+w2ftO4HwrIMigZWIiIhCV1iNjBAREVHgYTBCREREfsVghIiIiPyKwQgRERH5VcgFI8uWLUNGRgaioqJw/fXXo6CgoNX277//PjIzMxEVFYWrrroKW7du9VFPg58rz/qf//wnhg0bho4dO6Jjx44YNWpUm/9tqJmrf65N1q5dC0mSMGHCBO92MIS4+qwrKysxa9YspKSkIDIyEr179+bfI05y9VkvWbIEV155JaKjo5GWlobHHnsM9fX1Pupt8Nq9ezfGjRuH1NRUSJKEjRs3tnnNrl27cN111yEyMhKXX345Vq5c6d1OihCydu1aodFoxIoVK8R3330n7r//fpGQkCDOnDljt/3XX38tVCqV+Nvf/iaKiorEM888I9RqtThy5IiPex58XH3WkydPFsuWLROHDh0SR48eFffcc4/QarXip59+8nHPg4+rz9rk5MmT4rLLLhPDhg0Td9xxh286G+RcfdYNDQ1i4MCB4tZbbxVfffWVOHnypNi1a5c4fPiwj3sefFx91u+++66IjIwU7777rjh58qTYvn27SElJEY899piPex58tm7dKv70pz+JDz/8UAAQH330UavtT5w4IWJiYsTjjz8uioqKxD/+8Q+hUqnEtm3bvNbHkApGsrOzxaxZs8zHBoNBpKamioULF9ptP3HiRHHbbbdZnLv++uvFgw8+6NV+hgJXn7W1xsZGERcXJ95++21vdTFkuPOsGxsbxZAhQ8Sbb74ppk2bxmDESa4+6+XLl4uePXsKnU7nqy6GDFef9axZs8TIkSMtzj3++ONi6NChXu1nqHEmGHnqqadE3759Lc5NmjRJjBkzxmv9CplpGp1Oh4MHD2LUqFHmc7IsY9SoUdi7d6/da/bu3WvRHgDGjBnjsD0ZufOsrdXV1UGv1yMxMdFb3QwJ7j7rBQsWICkpCdOnT/dFN0OCO89606ZNGDx4MGbNmoWuXbuiX79+eO6552AwGHzV7aDkzrMeMmQIDh48aJ7KOXHiBLZu3Ypbb73VJ30OJ/743RgUG+U549y5czAYDOjatavF+a5du6K4uNjuNRUVFXbbV1RUeK2focCdZ23tj3/8I1JTU23+wJMld571V199hbfeeguHDx/2QQ9DhzvP+sSJE9i5cyd+//vfY+vWrTh+/DhmzpwJvV6PvLw8X3Q7KLnzrCdPnoxz587hxhtvhBACjY2NeOihh/D000/7osthxdHvxurqaly6dAnR0dEev2fIjIxQ8Fi0aBHWrl2Ljz76CFFRUf7uTkipqanBlClT8M9//hOdO3f2d3dCnqIoSEpKwhtvvIEBAwZg0qRJ+NOf/oTXXnvN310LObt27cJzzz2HV199Fd9++y0+/PBDbNmyBX/+85/93TXygJAZGencuTNUKhXOnDljcf7MmTNITk62e01ycrJL7cnInWdt8sILL2DRokX4/PPPcfXVV3uzmyHB1Wf9448/oqSkBOPGjTOfUxQFABAREYFjx46hV69e3u10kHLnz3VKSgrUajVUKpX5XJ8+fVBRUQGdTgeNRuPVPgcrd571s88+iylTpuC+++4DAFx11VW4ePEiHnjgAfzpT3+CLPPf1p7i6HdjfHy8V0ZFgBAaGdFoNBgwYAC++OIL8zlFUfDFF19g8ODBdq8ZPHiwRXsA2LFjh8P2ZOTOswaAv/3tb/jzn/+Mbdu2YeDAgb7oatBz9VlnZmbiyJEjOHz4sPk1fvx4jBgxAocPH0ZaWpovux9U3PlzPXToUBw/ftwc8AHA999/j5SUFAYirXDnWdfV1dkEHKYgUHCLNY/yy+9Gr6XG+sHatWtFZGSkWLlypSgqKhIPPPCASEhIEBUVFUIIIaZMmSLmzJljbv/111+LiIgI8cILL4ijR4+KvLw8Lu11kqvPetGiRUKj0YgNGzaI8vJy86umpsZfXyFouPqsrXE1jfNcfdalpaUiLi5OPPzww+LYsWNi8+bNIikpSfzlL3/x11cIGq4+67y8PBEXFyfee+89ceLECfHZZ5+JXr16iYkTJ/rrKwSNmpoacejQIXHo0CEBQLz00kvi0KFD4tSpU0IIIebMmSOmTJlibm9a2vvkk0+Ko0ePimXLlnFpr6v+8Y9/iPT0dKHRaER2drbYt2+f+b2bb75ZTJs2zaL9+vXrRe/evYVGoxF9+/YVW7Zs8XGPg5crz7p79+4CgM0rLy/P9x0PQq7+uW6JwYhrXH3We/bsEddff72IjIwUPXv2FH/9619FY2Ojj3sdnFx51nq9XuTn54tevXqJqKgokZaWJmbOnCl+/fVX33c8yPzrX/+y+/ev6flOmzZN3HzzzTbXXHPNNUKj0YiePXuK//3f//VqHyUhOL5FRERE/hMyOSNEREQUnBiMEBERkV8xGCEiIiK/YjBCREREfsVghIiIiPyKwQgRERH5FYMRIiIi8isGI0RERORXDEaIiIjIrxiMEBERkV8xGCEiIiK/YjBCREREfvX/Af6jSdLJxXMoAAAAAElFTkSuQmCC" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.errorbar(c_sticking_probability.p, c_sticking_probability.fd_mean, c_sticking_probability.fd_std, fmt=\"o\")\n", "plt.errorbar(rust_sticking_probability.p, rust_sticking_probability.fd_mean, rust_sticking_probability.fd_std, fmt=\"o\")\n", "plt.legend([\"C\", \"Rust\"])" ], "metadata": { "collapsed": false } }, { "cell_type": "markdown", "source": [ "Hmmm even with a different RNG I still get a dip around 0.6, I think theres something hiding there" ], "metadata": { "collapsed": false } } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.6" } }, "nbformat": 4, "nbformat_minor": 0 }