diff --git a/ErrorAnalysis/ErrorAnalysis.ipynb b/ErrorAnalysis/ErrorAnalysis.ipynb
index 87588e9f32a3554d32f4067b6570c2ebe9ef4af3..d65a16d4bc9fdeb619b2b170053d2808739913d1 100644
--- a/ErrorAnalysis/ErrorAnalysis.ipynb
+++ b/ErrorAnalysis/ErrorAnalysis.ipynb
@@ -469,7 +469,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 57,
+   "execution_count": 129,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -507,6 +507,11 @@
     "        \n",
     "        self.linapprox=scipy.interpolate.interp1d(self.tvals,self.y_approx_vals,copy=True,bounds_error=True)\n",
     "        \n",
+    "    def inc_approx(self,t):\n",
+    "        if not (min(self.tvals)<=t<=max(self.tvals)):\n",
+    "            return numpy.nan\n",
+    "        return self.linapprox(t).item()\n",
+    "        \n",
     "    def invert(self,yval):\n",
     "        if not (min(self.y_approx_vals)<yval<max(self.y_approx_vals)):\n",
     "            return numpy.nan\n",
@@ -517,7 +522,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 86,
+   "execution_count": 130,
    "metadata": {},
    "outputs": [
     {
@@ -613,7 +618,21 @@
       "C:\\Users\\richa\\Anaconda3\\lib\\site-packages\\scipy\\optimize\\_minimize.py:513: RuntimeWarning: Method COBYLA does not use Hessian-vector product information (hessp).\n",
       "  'information (hessp).' % method, RuntimeWarning)\n"
      ]
-    },
+    }
+   ],
+   "source": [
+    "fdict={}\n",
+    "for rank,df in data_by_rank:\n",
+    "    df_sorted=df.sort_values(by=\"post_sparsity\",axis=0)\n",
+    "    knots=list(zip(df_sorted[\"post_sparsity\"],df_sorted[\"pre_error\"]))\n",
+    "    fdict[rank]=monotone_invert(knots)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 131,
+   "metadata": {},
+   "outputs": [
     {
      "data": {
       "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiEAAAEaCAYAAADHQod0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XlcVOX+B/DPM8MM64iggCACbojAiIKhWaa5b5hW7kViWpH+romamV3lWnltkVKLyso17VpZmkZ6vWZm5gYqKghugRub7DuzPL8/zgwOwwADDgzL9/16nRczz3nOOc8cBuY7z8o45yCEEEIIaWoicxeAEEIIIW0TBSGEEEIIMQsKQgghhBBiFhSEEEIIIcQsKAghhBBCiFlQEEIIIYQQs6AghLQojLGxjLHzjLFSxhhnjHmZu0w1YYx5aco429xlMUZLurcPQ/PaIs1dDkIIYGHuAhBiLMZYBwDfAbgB4P8AlAHIMmuhADDGIgDkcM63mrssDdVc721TYIx1AxAKYC/n/IK5y0NIW0JBCGlJHgFgB+AdzvkP5i6MjggA1wFs1UtPBWANQNHUBWqA5npvG4M1AKXO824AVgFIAUBBCCFNiJpjmgHGmE1D9pnqGi2Is+ZnnllLYSQuKOOcq8xdFiO0qHurizEmZoxJjc2v+Z0o685JCGlsFISYCGPMljG2hjF2gzFWzhi7xxiLZow56OX7nTF2hzHWizEWwxgrABCj2Repaa8OZIx9zhjLBFCsc2x7xtgGzfHlmmutZoxZ6l1jq+Y8XRhj/2GM5QJIqKP8LzDGDmnKXcEYu6W5lkwvn5gxtpwxlsQYK2GM5TDG4hhj8424RxMZY3sZY7c15U/XlLWTEcemANimeXpY8/p+13m9KQaOGarJN1QnTXuP5YyxjxhjmZrX8StjzNPAOToyxtYzxv7WlDmNMfYTY8xPs58D6AxgiOa8XFuWmvqEMMY6a8qcoTlnImNsEWOM6eXTvlc8NPetUHO/P9f/nddy3+p8z9R2b2s57zOMsVOMsTzGWBFjLJkx9rleHs4Y+4YxNp496GtykzG20MD5jH3/aX+n8xhjSxhjNwCUAxhUz3JFah7PBnBYs2uLzu8wkjH2iubxIAPl1Zbj5druEyGkdtQcYwKaf+hHAPgB2AQgGYA3gPkABjLGHuWcl+scYgPgfwAOAVgCQP+b8nYAaQBWA3DQu0Y/AF8DOA/gCQD/1KSFGCjarwCuAHgTQF0fWv+A0KSwHkAugEAALwOQA3hSJ99KzbYVQBQAKwC+AAYD+LSOa8yB8J77HEAmAB8A8wAMYIz11btH+l4DMBbASwDWQgiqMuq4Xm02Q3idqwG4QmhS+QbC6wAAMMacAJwG4AHh9cYBaA/hfgRpyvA8gA2asryrObSoposyoe/FXwA6QbhfNwFMgHAvuwNYoHeINYTf+zEASwE8CuH3kgXhd1+jerxn6nVvGWPDAXwP4A8Ab0FobuoGw+/Bvpr0zwFsATAFwMeMMSvO+Xs6+Yx9/2kthPBe2gSgBEBaPcul9YfmNb+hOddxTfpFCM1pUQBegPA70xUKod/M7lrOTQipC+ectofcALwOoAJAoF76RAAcwFydtN81aW8YOE+kZl8MAKa371XNvsV66R9p0ifopG3VpH1ej9dgayDtBc15BuqknQfwSwPvk6FrDNFcY7oRx8/W5B2hl74VQIqB/EM1+YcauMc/6eWN0KT76qR9pUkbb+DcTOfxHQC/G8jjpTl+tk7a+5q0Z3TPBeBHTbrcwHvlNb3z/gwg04j7VZ/3jMF7W8N5PwKQD8Cijnxcs43VSZMAOAUhcHBswPtP+zvNAGD/EOWK1Hk+Qv/3pLNvJ4QmKiudNGsABQC+bcjfAW200fZgo+YY05gOIBbALU31fUfGWEcAJwGUAhhu4JjPajnf55xz/eWNJ0JomtGvbXhfZ7++6DpLrsE5LwYAxpiIMWavKf8fmt2P6GTNB+DHGOtt7LkNXIMxxtpprpEA4Z/8I7UebHr69/+o5md3QLgPAJ4F8Cfn/Bf9gw38fow1EcB1zvkevXN9oHmq/61dDeALA2V10m+qqOFa9X3PGCMfgC2AMUbkTeac/6p9wjlXQKg5sobw4a9NN/b9p7WTc57/EOUy1mYA9gCe0kl7GoAM1TsiE0LqiYIQ0/CBUE2epbdlQvhn66yXP8fAP1BdfxtI8wLwN+e8TDeRc54G4UO8q5HnMYgxFswYOwzhQysPQvlvanbr9mv5J4R/wImavgyfMMYMVZcbuoY3Y+xHCN8i8/HgPrXXu0ZTSNV7nqv56aj56QThw+eiia/rBSDJQHqi5qf+7zGTc16ql6Zf1tquVd/3jDGiIbyG/Zo+HLsYYzMZYxIDea8aSEvW/Ky8fj3ef1qG3tv1KZexfoMwauYFnbRQAHfxoC8JIaSBqE+IaYgg1HqsrGF/rt5z/Q8VfTXtr+nbN6thX13XEQ4WJqU6CuGf7euan6UAxAAOQidY5ZwfZ8K8CuMhfJOdDGA+Y2wz5/zFWq4hg9CvQQHgXwCuQaiS5wD+g4cLiGu6L+JajqlpxArTe97QGo/a1HZO/X21jazRL2t9rlXTe6buE3KeyRjrB2AYgNEARgGYAeB1xtjjnHPdPjF1XqM+7z8d1d7b9SyXUTjnnDG2DcBbjDFXCPdtOIAPOOfq+p6PEFIVBSGmcR1Ae875/xrxGikAHtN06Kv8ZsuEkSX2mv0N9RSEzrLjOeeV52GM9TKUWVOLswvALsaYBYSRFXMYY+9xzg198wWED4ZOAJ7knP+ucw1rPHwtSG4N5+j2EOfMglBb08eIvPX5ME8BYKgpq7fOflNJQSO9ZzTNKoc0Gxhjr0Jo9pkJoYOnlqH3kDZNW5tRr/eficpV5bA6TrsVwheMWRACIjGoKYYQk6DmGNP4FkBvxthM/R1MGNJaV7W5MfZDmEwqXC99qc7+htJ+29Z/P7yhn1EzuqMSF+ZbuKR5WtvrrOkarxtIq69rANoxxoJ0yimF0DGzQTTfcn8AMJgxNlZ/P2NVhtMWwfhAaj+AHoyxyXrnWqKz31Qa5T2j/x7QOKf5qf8e6KV7/zRNI/+AMLJEG7Qb/f4zYbl0aWtIDP4ONYHRUQhNMs8DOMU5TzaUlxBSP1QTYhrrAIwD8A1jbAKEphkA6AHgGQjDBbc+5DW+gjDEdR1jzAfCzI6DIVQ3H9BsDXUQwodCDGPsCwgfCiEADP1Tv8IY+xPAWQgjFHpBGFaaDGEIa01OQKhd2MEY2wihX8gwCJ0Osx+i7IAQBP4bwF7G2MeatOfw8DOVvgmhyelnxtgWCB9oMgjl/hbCUGpAuBfPMcZWQegDUcQ5r+kDfi2AqQC+ZYxph+iOhzBE9lPO+aUajmuIxnrPfMUYc4Yw/PcWgI4AXoHQRPKjXt5EAP9hjH0GoR/FVAADAKzgnOdo8tTn/WeqculKgNA0GM4YKwJQCOAy5/yyTp4tAHZoHr9Sz3IRQmpi7uE5rWWDMF/GCgCXIfxDzQMQD2EkgodOvt8B3KnhHJEQqoZ71LC/PYCNEP6ZV0D4AHsbgKVevq2a89Q6VFHvmOEQ5sQohhAsbIXQOVN/OONyCHMm3Ne8zhsQRju4GHGNQAgd/QogNKH8CKFzYgqArUYcPxs1DCOFMP9FLISJq+5AqD4fjpqH6PbQO94LBoZpAnCB0OHxtuae3wOwB1WH8nYG8IvmdXFohgvXcs7OEJqwsjTlvQJhiLD+sGyD7xWd++BlxD0z9j1T4701cM5nIAwjT9OUX3tP+unl4xDmXhkPYWh3GYQmmIiHeP8Nhd6w9waWK9LAsZc198jQfmsIf9OlEJpezf4/hzbaWsPGOG+MfneEkLaOCbPJ7uScP2fusjwszcRvaQAOcc5nmLs8hLQW1CeEEELqNgVCn5HN5i4IIa0J9QkhhJAaMMaGQej3tBJCk1JjjoAjpM2hIIQQQmq2EsBjEDpdz+HUfk2ISVGfEEIIIYSYRYuoCenYsSP38vIydzEIIaRFiYuLu885dzJ3OQipSYsIQry8vBAbG2vuYhBCSIvCGNNfI4mQZoVGxxBCCCHELCgIIYQQQohZUBBCCCGEELNoEX1CCCGEtG5xcXHOFhYWXwHwB31Bbi3UAC4rlcq5QUFBmYYyUBBCCCHE7CwsLL7q1KlTbycnp1yRSERzR7QCarWaZWVl+aanp38FYKKhPBRtEtLIOAeKi81dCkKaPX8nJ6cCCkBaD5FIxJ2cnPIh1G4ZztOE5SGkzSgsBPbtA8LDge7dgWnTzF0iQpo9EQUgrY/md1pjrEHNMYSYgFoNXLgAHDokbCdOAErlg/337gH/+hewapX5ykgIIc0N1YQQ0kCZmcA33wDPPw+4ugJBQcCbbwLHjlUNQACgvBz4+GPhJyGkeWKMBc2bN89d+3zlypUuERERbvU5R05OjsjZ2blPaGiohzbt+PHjNt7e3r4eHh7+s2fP7qJWq01Z7BaNghBCGmDhQsDFRQhAvvlGCEjq8uKLgKVl45eNENIwUqmUx8TEOKSlpTW4lWDx4sWdBwwYUKib9uqrr3pGR0enpqSkXL5586bVDz/80O7hS9s6UBBCSAN4etYvv0gE5OQ0TlkIIaYhFot5aGho1po1a1wacvzx48dtsrKyJCNHjizQpqWmpkqKiopEI0aMKBaJRJg1a1b23r17HUxX6paN+oQQokOpBP77X6CiApg0qeq+GzeE/h4HDwL/+1/d5+rSBRg9GhgzBhg+HGjfvnHKTEhrM2cOuly+DBtTntPfHyWbN+N2XfmWLl2aKZfL/SIjI9N10z/77DPH9evXd9LP7+XlVXbw4MGbKpUKixcv7rJr166bMTExlTUdqampEldXV4X2uaenZ0VaWprkYV9Pa0FBCCEa9+4Bo0YBCQmAnx8wciTw++9C0HHwIHD9upCvWzdg9mxgz56qzTBWVsCQIQ8CDx8fgDFzvBJCSEM5Ojqqp0yZkr127Vpna2vrys4b4eHhOeHh4TXWZ7733ntOo0aNyuvRo4dCN53z6gN+GP1jqERBCCEab70lBCCA8NPBAVAoABsb4MknhX4gY8YAPXoIedRq4I8/hLTRo4EnngCsrc1XfkJaC2NqLBrT8uXLMwIDA32nT59+X5tWV03IqVOn7M6ePWu3ZcsW55KSEpFCoRDZ2dmpli1blqlb85Gamirt1KmTQv88bRUFIYRACDb27Kma9uSTwJIlwODBQi2Hvo0bAQlVqhLS6ri4uKhCQkJyd+3a1XHGjBnZQN01IT///PPf2scbNmzoEBsbaxsdHX0XAGxtbdVHjhyxffLJJ4t37tzZYf78+UZ0ZW8bqGMqIQD+/BMoKKia9t//AkOHGg5AAApACGnNVqxYkZ6Xl2eSL+rR0dGpr7zyipenp6e/l5dX+ZQpU/JNcd7WgGpCCAFw5Ej1NImEAg1C2pKSkpLz2sddunRRlpaWnq8tf03+8Y9/ZAPI1j5/4oknSq5du5ZggiK2OlQTQgiA1auBo0eFobQumsF5X35p3jIRQkhrRzUhhEAIPu7fFzqb7t4N9OwJyGTmLhUhhLRuFIQQAiAxEXj7bWFEzGOPARb0l0EIIY2OmmNIm5OdDahUwmOFQgg++vUD7t4FtmyhAIQQQpoKBSGkzSgsBD76COjVC/jsMyAuDujfH1i5Epg8WagNeeopc5eSEELaDvrOR1q9u3eFOT0+/xzI1wyMW7xYmKLdxQXYu5eCD0IIMQeqCSGt1qVLwvTqXbsC7733IAABhLVh3N2FmVEpACGEAABjLGjevHnu2ucrV650iYiIcKvPOXJyckTOzs59QkNDPbRpx48ft/H29vb18PDwnz17dhe1Wl3tuPfff9/pk08+6fBQL8CEUlJSJGPGjOnW2NehIIS0GgUFwPHjQq3H6NFAnz7Atm1Cvw9Dbt0ybiE6QkjbIJVKeUxMjENaWlqDWwkWL17cecCAAYW6aa+++qpndHR0akpKyuWbN29a/fDDD+30j3v99dezFixYkK2fXl9KpfJhTwEA8PLyUhw8ePCmSU5WCwpCSIvDOXDnDrB/v9Cp9JlngO7dAXt7Yf2Wf/xDmO20LjNnCvkJIc1PRATcGEOQqbaICNRZoyEWi3loaGjWmjVrXBpS5uPHj9tkZWVJRo4cWTn/cmpqqqSoqEg0YsSIYpFIhFmzZmXv3bvXofrrjXBbuXKlCwAEBwf3Cg8P7yyXy3t7eXn5Hzx40A4QAoyXXnrJ3dvb29fb29v33XffdQaAzp07y5csWeIaFBTUa/PmzQ4JCQmWgwcP7unn59c7KCio1/nz560AYNeuXfZ9+vTx6d27t++gQYO8b9++bQEAv/zyi52Pj4+vj4+Pb+/evX1zc3NFycnJ0p49e/oBwjT0o0aN6j548OCenp6e/q+88kplbdFHH33U0cvLyz84OLjX9OnTPXVrgIxBfUJIs1VSIqxse+8ekJoKxMcDFy4IW3YDvy9YWgKhoUBEhLDKLSGE6Fq6dGmmXC73i4yMTNdNr2sBO5VKhcWLF3fZtWvXzZiYmMqajtTUVImrq2tlfaynp2eF7oJ2NVEqlezSpUtXdu/ebb969Wq3MWPGXF23bp1TamqqZUJCQqJEIkFGRoZYm9/KykodFxeXDACPPvqo96ZNm1Llcnn5b7/9ZhseHu5x6tSpqyNHjiyaPn16kkgkQlRUVMfVq1d3+vLLL++sW7eu04YNG1JHjRpVnJ+fL7KxsVFnZlZd3iYxMdEmPj4+0draWt2jRw//JUuWZFhYWODDDz90PXfuXGL79u3VgwYN8vbz8yutz/1utCCEMdYFwHYAnQCoAWzinK9njO0G0EuTrT2APM5538YqB2nejh4VRqXcuyd0INUGHXfvAnl5pruOoyMwf76wuTToOw4hpC1wdHRUT5kyJXvt2rXO1tbWlZ036lrA7r333nMaNWpUXo8ePao0AHPOq+VljNVZjilTpuQCwKBBg4qXLl0qBYDffvut3SuvvJIl0awn4eLiotLmDw0NzQWA/Px80fnz5+2mTJnSXbuvoqKCAcDff/8tnTRpkntWVpakoqJC1KVLl3IAGDhwYNGSJUu6TJ06NWfGjBm53bt3r9Zp5fHHHy/o0KGDCgB69OhRduPGDcvMzEyLAQMGFGrLMXny5NyrV6/WsNqWYY1ZE6IEsJhzfo4xJgMQxxg7zDmfps3AGFsHgBbyaSUuXwaSk4WhsNpt9GggKKjmYz74APj118YrU48eQq3HCy8ANjaNdx1CSOuxfPnyjMDAQN/p06ff16bVVRNy6tQpu7Nnz9pt2bLFuaSkRKRQKER2dnaqZcuWZerWfKSmpko7depUQ0+1B6ysrDgAWFhYQKVSMUAIaBhj1aMaADKZTA0AKpUKMplMmZSUlKifZ8GCBR4LFy5MnzVrVv6BAwdkq1evdgOANWvWpE+aNCl/37599oMGDep98ODBqzY2NlUCEalUWnldsVjMFQoFMxRg1VejBSGc8zQAaZrHhYyxKwA6A0gEACaEglMBDGusMjRnp08LnSh1BQfX3kfB0DGmxLnQBKIbRBQUPHi8ZInQj6Imn30GREdXTZPJhCBErRamRU9Lq7qlppqu/IwJQUffvsI2YICwCq5YXOehhBBSycXFRRUSEpK7a9eujjNmzMgG6q4J+fnnn//WPt6wYUOH2NhY2+jo6LsAYGtrqz5y5Ijtk08+Wbxz584O8+fPz6zpPLUZMWJEweeff+40fvz4Qm1zjG5tCCDU5Li7u1ds3rzZYc6cOblqtRqnT5+2fvTRR0sLCwvFHh4eCgDYunVr5UichIQEy+Dg4NLg4ODS06dP216+fNkqODi4pK7yDB48uHj58uVdsrKyxO3bt1ft27fPoXfv3s2jOUYXY8wLQD8Ap3WSBwPI4Jxfa4oyNDdHjwLLl1dNW7as9iDE0DFN6cYNoakkP19oKtH/ee5c9WPeeUcYHpuRIczLoc/SsmFlsbISRr9oA46+fQG5HLCza9j5CCHNS1QU7kVF4Z65rr9ixYr0bdu2OZniXNHR0akvvvhi17KyMvbkk08WTJkypUEtAIsWLcq6evWqpY+Pj5+FhQV/4YUXst58880s/XzffvvtzXnz5nm+9957rkqlkk2ePDnn0UcfLV2xYsW9GTNmdHdxcano379/8a1btywB4P3333f+66+/2olEIu7t7V367LPP5t+6davOfitdu3ZVLFq0KO2RRx7p7ezsrPD29i61t7dX1XWcLpNUp9R6AcbsABwD8C7n/Eed9M8AXOecr6vhuJcAvAQAHh4eQamm/MrciJRK4MQJ4QPSQaf/c0UF8McfwoiOixeF4aE39QY/eXkJs3mKRMK3ev3t2jXgypUmfTn1IhIJNR66evYEHn8c6NQJcHWtvm3fDrzyiuHzicVCns6dATc3oZajXz8h4OjZk6ZXJ6QujLE4znl/c5fDGPHx8SkBAQH3685JmpP8/HyRvb29WqFQYPTo0T1mz559PzQ0tEqPvvj4+I4BAQFeho5v1H/jjDEJgD0AduoFIBYAngZQY28BzvkmAJsAoH///o0bKT2k3Fzg4EEhwPj1V6FWYOtWYNw4ICZGSP/vf4UmDSsr4YPUwFw1KC8XzsW5sKnVDx5zLuwzp0GDhJEl9vZA+/bCT93HX38NLFxY9ZiRI4FPP635nH37Ai++KAQZ2mBD+9PJiZpSCCGkOVu6dKnbH3/80a68vJwNGTKk4LnnnqvXkILGHB3DAHwN4ArnPEpv9wgASZzzO411/cZ27ZoQXOzfL/TTUOlVQEVEAGFhQvDg5gbMmAFMmAAMHy50kFy7tnrTSmiokF4TQ8c0pV69gJdfrnm/j48w+2i7dkJfEJlMCFxqM2CAsBFCCGl5Nm3a9FCf441ZE/IYgOcBXGKMXdCkvck5jwEwHcC3jXhtk1Mqgb/+ehB4JCfXnj8vD3jrLWDSJKHmQ39EVnCwsH6JrsGDaz+noWNMzcqqahAhkz147upa+7GjRgkbIYQQYozGHB3zJwCDg6E557Mb67qmlJcHHDokBB0xMfVrDlGrhb4QgYGG9w8bJmz10ZBjCCGEkOaKuvbpuX69ajNLQ6bht7UV5sdoV211AEIIIYRotfkgRKkETp58EHgkJTXsPB4eQEiIsA0d2vChp4QQQkhb0SaDkPz8qs0sOTVOP1MzxoQ+GtrAQy6v3u+DEEJIy8EYC5o7d27Gl19+eQcAVq5c6VJUVCSOiooyar4SsVgc1LNnz1IAcHNzq/jtt9+uA0BSUpJ06tSp3fLz8y38/f1L9uzZ87d2RtS2rs0EITduCEHHgQPAsWMNb2YZOVIIOsaPpzVICCGkNZFKpTwmJsYhLS0t3dXVtd6fEpaWlmpD06VHRES4L1iwIOOll17KnTlzpsf69es7Llu2rNokY22RyNwFaExXrgizkPr6ChNdLVoEHDlSvwCkSxfg1VeF+T/u3wd++gmYM4cCEEIIaW3EYjEPDQ3NWrNmjcn+w6vVapw8eVIWFhaWCwBz5szJ3r9/f3tTnb+la9U1IRcvAu+/X//jdJtZ+vShZhZCCGlKc/bN6XI587JJl5z0d/Yv2fzU5tt15Vu6dGmmXC73i4yMTNdNr2sBOwCoqKgQ+fv79xaLxXzJkiXpzz//fF5GRoaFTCZTaVe+9fLyqsjIyJCa6GW1eK06CElPrzsPIEweNmqUMJnY+PHCFOOEEELaHkdHR/WUKVOy165d62xtbV05t3VdC9gBwPXr1y96eXkpEhMTpSNHjuwVGBhY6uDgUG0tlZpWwm2LWnUQUlIiTCmeZ2ASWXf3B7UdTz4pTNJFCCHE/IypsWhMy5cvzwgMDPSdPn165Vo2xtSEeHl5KQDA19e3YuDAgYVnzpyxeeGFF3ILCwvFCoUCEokEKSkpUmdnZ0XTvZrmrVX3CVm+HFi16sHzRx4BVq8Gzp8XFpCLjgbGjqUAhJA9iXvw9O6noVKr8MGJD/Bz8s/mLhIhZuPi4qIKCQnJ3bVrV0dtWnh4eE5SUlKi/qYNQLKyssSlpaUMANLS0ixiY2Pt+vTpUyoSiTBw4MDCLVu2OADA5s2bO0yYMKFe66u0Zq26JgQAnn5aWN59/Pi6px0npK1RqBRYfmQ51p0UFrN+Ye8L2HVpF17o+wIm9ppo5tIRYj4rVqxI37Ztm5Ox+S9cuGA1f/58T8YYOOd47bXX0oOCgsoAYN26dXemTZvW/Z133uns5+dXsnDhQlotWINx3vybpvr3789jY2PNXQxCWpRzaefgaucKV5nh6DutMA3TfpiG47eOV0l3b+eOpPlJsJXaNkUxSSNijMVxzvubuxzGiI+PTwkICKAP51YoPj6+Y0BAgJehfa26OYaQtiqjKAOPb34cbxx5w+D+YynH0O+LftUCEADIL8tHRnFGYxeREEIoCCGkNYo6GYVSZSlO3TlVJZ1zjg9OfIDh24fXGGiM6DYCHaw7NEUxCSFtXKvvE0JIW5Ndko3o2GhYii1xNfsq8svyYW9lj/yyfITtC8NPST8ZPI6B4b0R72HJoCVgNDkOIaQJUE0IIa3MhtMbUFRRhDXD1wAA4tLicDHjIvp/2b/GAETERNg/Yz+WPraUAhBCSJOhmhBCWpH8snxsOLMBY3qMQXFFMQBg5dGVOJd2DqXK0hqP2z9jP8b1HNdUxSSEEAAUhBDSqkSfjUZeWR5O3DqBg9cPAgBO3D5R6zFrh6+lAIQQYhbUHENIK1FcUYwPT34IKwsrFFYU1prXVmILC5EFxvYYi6WPLW2iEhLSvDHGgubNm+eufb5y5UqXiIgIN2OPF4vFQT4+Pr4+Pj6+w4YN66FNT0pKkvbp08fH09PTf/z48d3KysqqtXnu3LnT/s0332xWi4b069fPp7GvQUEIIa3EZ7GfIac0B2XKslrz9e7YG062TnC2dca2SdsgYvRvgBAAkEqlPCYmxiEtLa1BrQSWlpZq7Uyqv/3223VtekREhPuCBQsyUlNTL9vb2yvXr1/fUf/YWbNm5a9Zs8bIFc9qplCYbkb48+fPJ5nsZDWg5hhCWoEyZRkif4+sM99M+UxwzrE7YTeOhB6Bk63RE0IS0qQiDkW4fXTqI5MSwnCPAAAgAElEQVTNc71o4KK0qNFR92rLIxaLeWhoaNaaNWtcNm7ceNcU11Wr1Th58qRs3759NwFgzpw52ZGRkW7Lli3L0s23YcOGDrGxsbbbt2+/9cwzz3jJZDJVfHy8bVZWluTtt9++ExYWlgsAb731lst3333XgTGG4cOH50dHR98NDg7uFRwcXHT69Gm7cePG5b388svZYWFhnnfv3pUCQFRU1K1Ro0YVHz161CYiIsKjrKxMZGVlpd66devfAQEB5bGxsVZhYWFdFQoFU6vV2LNnzw25XF5uY2PTr6Sk5PyBAwdkq1evdnN0dFQkJydby+Xykr179/4tEomwe/du+zfeeMPd0dFRKZfLS1JTUy2PHj16vfqdMIyCEEJagdcOvoZiRXGN+0d0HYEVT6xAal4qZu+bjVVDVmGo19CmKyAhLcTSpUsz5XK5X2RkZJVaCWMWsKuoqBD5+/v3FovFfMmSJenPP/98XkZGhoVMJlNJJBJt/oqMjAxpXeXIyMiQxMbGJl24cMFq8uTJPcLCwnK/++67dr/88otDXFxckkwmU2dkZIi1+fPy8sRnz55NBoCQkJCuERERGaNHjy66du2adPTo0T1v3ryZEBAQUHbmzJkkiUSCvXv3yl5//XX3Q4cO3di4caPTq6++mhEeHp5TVlbGlEpltfJcuXLF+sKFCze9vLwUQUFBPocPH7YbPHhw8cKFCz1///33JB8fn4qQkJCu9b3fFIQQ0sIlZCVgU9ymaukMDL/O+hUR/42ApYUlXO1cMWHXBAzxHIJ/PvFPM5SUkObP0dFRPWXKlOy1a9c6W1tbq7Xp4eHhOeHh4Tm1HXv9+vWLXl5eisTEROnIkSN7BQYGljo4OKj08zHG6lwvZeLEiXlisRhBQUFl2dnZEgA4fPhwu+eee+6+TCZTA8JCe9r8M2bMqCzbiRMn2l27ds1a+7yoqEicm5srysnJEU+bNq1rSkqKFWOMKxQKBgCPPvpo8Ycffuh6584d6fTp03Plcnm5fnnkcnlx9+7dFQDg5+dXcuPGDalMJlN16dKl3MfHpwIApk+fnvPVV1/Vq3qVGoMJacGKK4oxesdocFT/n7b6ydUY3WM0HnF7BGfunsG0H6bBWmKNnU/vhFgkNnA2QggALF++PGPXrl0di4uLKz8jP/vsM0dtp1PdbcyYMd20eby8vBQA4OvrWzFw4MDCM2fO2HTq1ElZWFgo1vbVSElJkTo7O9fZccPKyqryj1q7xhvnvMZ5fLSBiTZfbGzsFW3/lMzMzIsODg7qZcuWdR4yZEjhtWvXEvbv33+9oqJCBACvvPJKzr59+65bW1urx44d6/3zzz/L9M9vaWlZWR6xWAylUmmStecoCCGkheKcY+7+ubhbWL3peoL3BLw5+E0AwCNujyCrJAvxGfHYNmkbOrfr3NRFJaRFcXFxUYWEhOTu2rWrsgNpeHh4jvZDXXfTNsVkZWWJS0tLGQCkpaVZxMbG2vXp06dUJBJh4MCBhVu2bHEAgM2bN3eYMGFCXkPKNWbMmIIdO3Z0LCwsFAGAbnOMrscff7zgvffec9Y+/+uvv6wBoKCgQOzu7l4BAF988UXla0tMTJT27t27/K233socNWpU3oULF6yrn7W6gICAstu3b1smJydLAWD37t2O9X1N1BxDSAv16dlP8Z/L/6mW3t2hO3ZM3lE56mVQl0EAgMWPLqb5QEiLETU66l5dHUkb04oVK9K3bdtmdNPChQsXrObPn+/JGAPnHK+99lp6UFBQGQCsW7fuzrRp07q/8847nf38/EoWLlzYoNWCn3322YJz587Z9O3bt7dEIuEjRozI/+STT6p9C9m0adPtuXPnenh7e/uqVCo2YMCAwkGDBt1atmxZ+ty5c7tu2LCh0+DBgwu0+Xfs2OH4/fffd7CwsOBOTk6Kf//730bddzs7Ox4VFZU6ZsyYno6Ojsp+/frV3DGtBiapTmls/fv357GxseYuBiHNxl+3/8KQrUOgVFftQGZtYY2TL55EQKeAKuln7p5BkGsQNcO0MYyxOM55f3OXwxjx8fEpAQEBDfpwJuaTn58vsre3V6vVaoSGhnr07NmzbNWqVZm6eeLj4zsGBAR4GTq+0ZpjGGNdGGNHGWNXGGMJjLGFOvv+jzGWrEl/v7HKQEhrlF6UjinfT6kWgADAFxO+qBaAAEBw52AKQAghJvfxxx939PHx8e3Zs6dfQUGBOCIiol6BZGM2xygBLOacn2OMyQDEMcYOA3AB8BSAPpzzcsaYc61nIYRUupV/CxN2TcC9wuq1pa/2fxXPBzxvhlIRQtqqVatWZerXfNRHowUhnPM0AGmax4WMsSsAOgOYB2At57xcs6/BhSekLVGqlRixfQSu5Vyrtm+g+0B8NOYjM5SKEEIarklGxzDGvAD0A3AagDeAwYyx04yxY4yxR2o45iXGWCxjLDYrK8tQFkLaFAuRBT4c9SEYqg7Rc7JxwvdTvodUXOf8R4QQ0qw0ehDCGLMDsAfAa5zzAgi1Lw4ABgJYCuA7ZmDgM+d8E+e8P+e8v5MTTS1NSFFFEb4892WVOUFsJDb4YeoPcG/nXsuRbVhGBpCQYO5SEEJq0KhBCGNMAiEA2ck5/1GTfAfAj1xwBoAaQLXFfAghD9wpuIPBWwYj5loMPhn7CeY/Mh9uMjccDzuOJzyfMHfxmh+FAvjoI8DbG5g5E7hntpGehJBaNOboGAbgawBXOOdROrv2AhimyeMNQAqAhmURokPNKyc/xPm08xjw1QBcz7mOAzMOYH7wfHw85mPEvRSHQNdAM5aymTp8GAgIACIigIIC4OJFwNMT+OMPc5eMNHOMsaB58+ZVViuuXLnSJSIiws3Y469duyZ97LHHenbr1s2ve/fuftpJvJKSkqR9+vTx8fT09B8/fny3srIyw9OetkGNWRPyGIDnAQxjjF3QbOMAbAbQjTF2GcB/ALzAW8JkJYQ0ssLyQmy7sA3Dtw/HooOLAAD7k/dj8JbBEDMxTsw5gbE9xwIQ+od0squ2llbbdvMmMHkyMGoUcOVK1X0iEeBY78kcSRsjlUp5TEyMQ1paWoMGbcyaNavrkiVLMm7evJlw7ty5K25ubkoAiIiIcF+wYEFGamrqZXt7e+X69eup9l+jMUfH/Amgpmjvuca6LiEtzZ2CO3jjf2/gxys/olRZCgA4l3YOxRXF2BK/Bf069cP+GfvhKjPZquatS0kJsHYt8P77QHm1dbcEFRXA5s1AVJTh/YQAEIvFPDQ0NGvNmjUuGzdurL4eQi3i4uKsVCoVJk+eXAAA9vb2agBQq9U4efKkbN++fTcBYM6cOdmRkZFuy5YtoxEXoGnbCTE7mVSGPVf2oExZVpmWV5aHHRd34MV+L+Kj0R/BVmprxhI2c3l5Qv+PmgIQxoB584Dly5u2XKTh5szpgsuXbUx6Tn//EmzefLuubEuXLs2Uy+V+kZGR6brpn332meP69eurVT96eXmVHTx48GZiYqJVu3btVKNGjep++/ZtyyeeeKLg008/vZOVlWUhk8lUEolEm78iIyODhrJpUBBCiBllFmdi9+XdsJHYVAlCAGBsz7HYFLLJTCVrQSwsAD8/4PTp6vsGDQI2bgQCqe8MMY6jo6N6ypQp2WvXrnW2trau7JwVHh6eEx4enlPTcUqlksXGxtqdPn06sWfPnhUTJkzovnHjxo7Tpk2rtlgdY4y6IGhQEEJIEytTluGnKz/hm0vf4ND1Q1BxFbo5dENOadX/bzdyb0CpVsJCRH+mBpWVAevXA+++KzTJtG8v1IoAgKsr8MEHwsiYGpY+J82YETUWjWn58uUZgYGBvtOnT68cNFFXTYiHh0dF7969S319fSsAYOLEibmnTp2yW7hw4f3CwkKxQqGARCJBSkqK1NnZWdGUr6c5a5LJyghpq8qUZXj72NuIuxcHALiceRmBXwRi5o8zcTHjIpYMWoJL4ZeQvCAZzrbOcLByQHj/cJx88SQuvnKRApCaZGQItR9vvAEMHSrMBbJzJyCVCmnJycCsWRSAkAZxcXFRhYSE5O7atauyA2l4eHhOUlJSov528ODBmwAwZMiQ4vz8fPG9e/csAODo0aPtfH19S0UiEQYOHFi4ZcsWBwDYvHlzhwkTJlSrHWmr6D8cMSivNA+n756GmqthI7GBpdgSBuaUazRWFlYGF2LTKlWU4mLGxYc+pjHdKbiD5UeW41rONRy4egBh/cKw6NAi2FvaY9/0fZjgPQEi9uB7wP+e/x+8O3jD0sKyycrYYh04IIyG+eknYNIkIa1XLyGtc2fzlo20CitWrEjftm2b0TNlWlhYYO3atXeGDh3qDQByubxk0aJF9wFg3bp1d6ZNm9b9nXfe6ezn51eycOFCmpZCg7WE0bH9+/fnsbGx5i5Gq8U5R0peCuLS4hB3Lw5xaXE4cesESpQlZiuTdwdvJC9IrnF/8v1k+Hzq89DHNLWR3UZi++TtNLz2Yc2dC/z4I5CdTbUdtWCMxXHO+5u7HMaIj49PCQgIoA/nVig+Pr5jQECAl6F9VBPSRhSWFyIxKxGXMy/jr9t/AQCcbJ0QlxaHc2nnKvsjWIgs4O/sj8c8HsPhm4fNVt6s4iy888c7kEllkFnKIJPK0M6yXeXjzJKWt+6htYU1dkzeARc7F3MXpeU7dQoYOJACEEJaOApCWjHOOW7m3sRX577C2hNrq+2XiCTwd/bH0z5PI8gtCEGuQZC7yGFlYYWbuTfRfUN3M5RakFeWh38e/We9jvk792/0/bxvZaAis5ShnVQIXGwkNrhfbN4vWQq1AsdvHcezvs+atRwt2s6dQNeuQh+QqVPNXRpCyEOiIKSFKlWU4sr9K0jITECpshTzAufhRu6NyuYUbQ1HXpnh/k/2lvZIX5wOK4mVwf1e7b1gI7FBicI8TTI9O/TE5fDLKKwoRGF5IQrKCyofF1YU4mr2Vaz4bUWVY6wsrODZ3hOF5YXIKsnCjdwblflLFCWAGVseB3sMRvT4aPg7+5uvEC3Ix7dv43+5uTjQp8+DxPv3ged05jk8cACQSISOqFQjQkiLREFIM1emLEPy/WQkZCUgITNB+JmVgBs5NypXU5WIJHj98OvIL88HAEjFUsid5ZjqOxVBbkHo69IXT2x9AuWqB5M55ZfnI688D50khvsmiJgIE7wnICUvpdFfoyEe9h6QiCVwtHaEo3X16bZT81KxL3lftWO+n/J9jedMzUvF1B+a9tuzezt3TPWdiql+U5u0Y29Lt+jGjeqJf/1V9fnZs0BxMU1CRkgLRkFIM1GhqsDV7KuVgcblzMtIyErA9ZzrVRYzM0ShVuCpXk/hcY/HEeQWBH9nf0jFVSfk6+3UGxfSL1RJS8hMqLWD5O5ndzf8BTUyz/aeOD3XwORUJj6GNCMnTlRPe+yxpi8HIcRkKAhpYgqVAtdzrlfWbFzOuoyEzARcy7kGpVrZ4POG9QvDUK+hNe7v49IHheWF8HP2g5+TsFHTAGlRKAghpNWhIKSRqNQq3Mi9IdRo6DSjJN9PhkJt+snyEjITag1Ctj61lZoDSMs2YQLAedVmGQpCiAkxxoLmzp2b8eWXX94BgJUrV7oUFRWJo6Ki7hlz/LVr16SzZ8/2TEtLkzLGEBMTc61Xr14VSUlJ0qlTp3bLz8+38Pf3L9mzZ8/fVlZWVXqp7dy50z4hIcF6zZo16TWdv6n169fP5/z580mNeQ0KQkzkyM0jOH33dGUNR9L9pCp9MBpDl3ZdKms2Al1rXxuDAhDS7JWUAEeOAKmpgL9eLV12NtCjB+DsLDz//nshrbv5RnCR1kcqlfKYmBiHtLS0dFdX13pXTc+aNavr8uXL0yZPnlyQn58vEomEyQgjIiLcFyxYkPHSSy/lzpw502P9+vUd9VfRnTVrVj6A/Id9Ddrp4U2hsQMQgIIQk/n49Mc4cPVAo5zbTeZWpQnFz9kPvk6+aGfZrlGuR0iTSU0FfvlFGOny22/CSrjW1kBMjLB/yRIhMImPF2pB7OyAGTOAp58GRLTqRKsWEeGGjz5yNdn5Fi1KQx01GmKxmIeGhmatWbPGZePGjXfrc/q4uDgrlUqFyZMnFwCAvb29GgDUajVOnjwp27dv300AmDNnTnZkZKSbfhCyYcOGDrGxsbbbt2+/9cwzz3jJZDJVfHy8bVZWluTtt9++ExYWlgsAb731lst3333XgTGG4cOH50dHR98NDg7uFRwcXHT69Gm7cePG5b388svZYWFhnnfv3pUCQFRU1K1Ro0YVHz161CYiIsKjrKxMZGVlpd66devfAQEB5bGxsVZhYWFdFQoFU6vV2LNnzw25XF5uY2PTr6Sk5PyBAwdkq1evdnN0dFQkJydby+Xykr179/4tEomwe/du+zfeeMPd0dFRKZfLS1JTUy2PHj163dj7RkGICdzOv420wrSHPo+LrYsQZDj5VdZw+Dr5wsHawQSlJKQZUKmEicYOHBC2y5er5yktffB440ZhJdzVq4Fhw4BHHhGG5bZySmU+RCIrlJQkwcbGFyJR63/NzcXSpUsz5XK5X2RkZJVmkboWsEtMTLRq166datSoUd1v375t+cQTTxR8+umnd7KysixkMplKWzvh5eVVkZGRIdU/j76MjAxJbGxs0oULF6wmT57cIywsLPe7775r98svvzjExcUlyWQydUZGhlibPy8vT3z27NlkAAgJCekaERGRMXr06KJr165JR48e3fPmzZsJAQEBZWfOnEmSSCTYu3ev7PXXX3c/dOjQjY0bNzq9+uqrGeHh4TllZWVMqaxeCXTlyhXrCxcu3PTy8lIEBQX5HD582G7w4MHFCxcu9Pz999+TfHx8KkJCQrrW935TEGIilzIvGZ3XycapSgdR7eMONh0asYSEmFlCgjDPx4ULdefV4Lm5YDY2jVio5kepzMeffzpAKu2Eioo0ODvPgIvL8+jQYay5i9YmODo6qqdMmZK9du1aZ2tr68qhieHh4Tnh4eE5NR2nVCpZbGys3enTpxN79uxZMWHChO4bN27sOG3atGqTNTHG6py1aOLEiXlisRhBQUFl2dnZEgA4fPhwu+eee+6+TCZTA8JCe9r8M2bMqCzbiRMn2l27ds1a+7yoqEicm5srysnJEU+bNq1rSkqKFWOMKxQKBgCPPvpo8Ycffuh6584d6fTp03Plcnm1vgRyuby4e/fuCgDw8/MruXHjhlQmk6m6dOlS7uPjUwEA06dPz/nqq6+MXm8HoCDEJNzbuePPsD8R/FVwlXRHa8dqgYafsx+cbZ3NVFJCzOT2bWD4cGH123pQW1tDXHe2VqWoKB4AR0WFULuamfktCgvjKAhpQsuXL88IDAz0nT59euU0y3XVhHh4eFT07t271NfXtwIAJk6cmHvq1Cm7hQsX3i8sLBRr+2qkpKRInZ2d6xydoNtxVbvGG+e8xv592sBEmy82NvaKnZ1dlWBn7ty5HkOGDCk8fPjwjeTkZOmwYcN6AcArr7ySM3jw4OKffvrJfuzYsd7R0dEpEydOLNQ91tLSsvJcYrEYSqXSJGvPURBiAowx+Dv7Y17gvCr9NlxsXahDKCFFRUBIiHEBiLMzMG5c5VMV5xC3sb+hwsJz1dLs7PqZoSRtl4uLiyokJCR3165dHWfMmJEN1F0TMmTIkOL8/HzxvXv3LNzc3JRHjx5tFxQUVCwSiTBw4MDCLVu2OLz00ku5mzdv7jBhwgTDU1nXYcyYMQXvvvuu27x583K0zTG6tSFajz/+eMF7773n/Pbbb2cAwF9//WU9aNCg0oKCArG7u3sFAHzxxRcdtfkTExOlvXv3Lvfz88u8efOm5YULF6z1gxBDAgICym7fvm2ZnJws7dWrV8Xu3burzyxZBwpCTMRaYo1NIZvMXQxCmheVCpg1S+hYWpPAQGH47fjxQP/+QofT338XDm8Bq3ybWlHR+WppMlnto99apaioe3V1JG1MK1asSN+2bZvRTQsWFhZYu3btnaFDh3oDgFwuL1m0aNF9AFi3bt2dadOmdX/nnXc6+/n5lSxcuLBBC1k9++yzBefOnbPp27dvb4lEwkeMGJH/ySefVOtAu2nTpttz58718Pb29lWpVGzAgAGFgwYNurVs2bL0uXPndt2wYUOnwYMHF2jz79ixw/H777/vYGFhwZ2cnBT//ve/jbrvdnZ2PCoqKnXMmDE9HR0dlf369Suu72sySXVKY+vfvz+PjY01dzEIIfX1+uvABx9UT3d1FTqbjhsHuLlV2cU5h+jYMQDAGEdHfNKzJ7pbW1c/Ryt19qwcxcVVO+z26fNfODqOrPe5GGNxnPP+pipbY4qPj08JCAgw7yqTpN7y8/NF9vb2arVajdDQUI+ePXuWrVq1qsoy5/Hx8R0DAgK8DB1PNSGEkMbx9deGAxBra2D/fiAoyOBh2YoHzeUHc3Lw+Pnz6G5lhT8DW39tgEpViuLiK9XSqTmGNFcff/xxx2+//bajQqFgfn5+JREREfUKJCkIIYSYVkmJUMvx4YeG93/zTY0BCAAUqqo2cadXVEDaRvqFFBdfAlD19VtadoFU2tHwAYSY2apVqzL1az7qg4IQQojpxMQA8+cDKSmG969ZI0w0VosulpZVnm/y9oaHlZWJCti8GeoPYmfX+muASNtFQQgh5OHduwcsXAj88EPNeUJDgTfeqPNUFnozoQ5zcGgzfUIMjYyRyagphrReNO8xIaThVCrgk08AH5/aA5Bhw4BNm4AGNKtYtJGmGAAoKjI0PJdqQkjrVWcQwhgTM8YW1ffEjLEujLGjjLErjLEExthCTXokY+wuY+yCZhtX17kIIc3QuXPAwIHA//0fUFjDlAKWlsC//gX8+qvw2IAd6emIuF7zUhNtJQhRqxUoKqo+8zJ1SiWtWZ1BCOdcBeCpBpxbCWAx57w3gIEA5jPGfDX7PuKc99VsMQ04NyHEXAoLgUWLhHVcahs6P3w4cOkSsHIlIDW8VIaac6xMScEunYnM9KcNkDTTIESZr4Qit86JL41WUnIFnFedLVsicYKlZWeTXYPUjjEWNG/ePHft85UrV7pERES41XaM1v79+2U+Pj6+2s3S0jJwx44d7QEgKSlJ2qdPHx9PT0//8ePHdysrK2ueb2ozMLY55gRj7BPG2GDGWKB2q+0Aznka5/yc5nEhgCsA6K+JkJZs717A1xf4+GNArTacx9lZGAFz+DDQs2etpzuRn4+UsjJU6AQe+hOUNdeakKTZSUicmmiy8xnulNqPZl1uQlKplMfExDikpaXVu79kSEhIYVJSUmJSUlLisWPHkq2srNSTJk0qAICIiAj3BQsWZKSmpl62t7dXrl+/noY7aRgbhAwC4AdgNYB1mq2G8XfVMca8APQDcFqTtIAxdpExtpkxZnCJWMbYS4yxWMZYbFZWlqEshJCmcusW8NRTwOTJwJ07NeebNw+4ckWYJdWID8/tmhqQcp2ApqKF1IQUnitEydUS053PYKdU6g/SlMRiMQ8NDc1as2aNy8OcZ8eOHQ5DhgzJl8lkarVajZMnT8rCwsJyAWDOnDnZ+/fvb2+aErd8RkV7nPMnG3oBxpgdgD0AXuOcFzDGPgPwNgCu+bkOwBwD19wEYBMgzJja0OsTQkzg1Cng559r3u/vD3z+OfDYY0afskylwveZwvQCuoGHogXUhKhKVSi/VQ4mZbUuKlYfNdWEtEVzkpK6XC4uNunyyf62tiWbfXxu15Vv6dKlmXK53C8yMjJdN72uBex003744QfHhQsXZgBARkaGhUwmU0kkEm3+ioyMDMPtk22QUUEIY8wewCoAT2iSjgFYzTnPr+M4CYQAZCfn/EcA4Jxn6Oz/EsCBBpSbENKU+vcHOnQAsrOrpltbA6tWARERgOafrLH2Z2cjX6XC4/b2+DM/H2rOIWIMLycnV8k3+fJlvObujjEdOjzsqzCZ0uulAABewaHMUULSoX6vXR/napojpJlwdHRUT5kyJXvt2rXO1tbWlVV0dS1gp5WamipJTk62fvrppwuA6n2cAIAxRl+sNYxt99oM4DKAqZrnzwPYAqDGWYeY8NXgawBXOOdROumunPM0zdPJmvMSQpojzoGvvhKCDM6FQEM7rfrYscCnnwJduzbo1DsyMuAmlWKMoyP+zM9HhVoNxhh26zW/HsrNxQyXh6odN7mS5AfNMBXpFQ8dhJSW3oBKVVQlTSxuB2vrbg913pbKmBqLxrR8+fKMwMBA3+nTp1dOQW5sTcj27dsdxowZk2dpackBoFOnTsrCwkKxQqGARCJBSkqK1NnZ2XQ9mls4Y/uEdOecr+Kc39Rs/wJQ11/HYxCClWF6w3HfZ4xdYoxdBPAkgHoP/yWENIF794SVbV96CQgOBhIShOnYXV2B774DfvmlwQFIVkUFfs3JwUwXF1hrJicr5xy/6Ne0aMjE4ga/jMZQerW08nF5WnktOY2Tnr6tWpqdXV8wRlM5mYOLi4sqJCQkd9euXZUdSMPDw3O0HU91N0NNMTNnzqysMRGJRBg4cGDhli1bHABg8+bNHSZMmJDXdK+meTP2HV7KGHtc+4Qx9hiA0lryg3P+J+eccc776A7H5Zw/zzmXa9In6tSKEELMpawM2LFDqO3gHNi1S+jn8fvvwMaNwkgXT09g8WKh4+mUKQ2aeEzrP5mZUHKOUBeXynVhKtRqxNYw30ignV2Dr2UKnHNcnX8VheeE8pVcLan871meVgKFIrfB505L24Jbt96tlt5W+4M0FytWrEjPy8ur1yiZ5ORkaVpamnTcuHFV3sjr1q27s3Hjxk4eHh7+ubm5FgsXLqTVgjWMvcGvANiu6RsCALkAXmicIhFCmoxCARw8CCxZAly9Kjz/9Vdh9tOBA4Ft2wBv7wf5JRLA3r7m8xlpR0YGAmxtIbezw8l8oWtZjkKBOANBiJNEAi8zT9uuzFfiXvQ9SDpKIAuUoTS5FLJAGdWTE2AAACAASURBVApjC3FftQU3Tn+MQYPSIRLVr1kmO/tXJCfPM7ivU6dQUxSd1ENJSUllx5wuXbooS0tLq3fUqUWvXr0qMjMzL+qn+/r6Vly6dKn68sik7iCECfWBvTjnAYyxdgDAOS9o9JIRQhpXebnQtJKr8y1+7lzAwgJYu1YITBqhGSS5pARnCwvxYffu+C4zE2/+/TcAwOfsWYP5D/XpY/Iy1Je6VOifqMgUmvJLrpbAaYoTiq8Uo4xfhVKZA4Uis14TixUUnEVCwrPQXzUXALp0WUrDc0mbYMyMqWoACzSPCygAIaSVsLSsWssBCE0xY8YAy5Y1SgACCNO0iwDIbW0x68oVZCuVteb3tbVtlHLUh7pECEIqsiqgyFZAmaOETS8bWLpaQqkWOtKWlxvfslxaegOXLo2HWl19nhFn55no1m2taQpOSDNnbJ+Qw4yxJZr1YBy1W6OWjBDSMHfuANu3Ay+8AMTFVd2XkyPM97F0qdDccuZM9eMPHQI0tROmpuYcOzIyMNLBAVG3b0NpYPiiPkuR+TtnqkqF2gpFlqJyZIyNtw2knaRQioUgpKLC+CDk9u0PoVBUn4SxXbtH4eOzhTqkkjbD2D4h2snE5uukcdQ9QoYQ0tju3weOHgV++w04cgS4du3BPmdnIDkZOH5c2BIShHSpVFj7Zdo04D//eZB/wADg668bPOqlLsfz83GrvBzPubhgza1bjXKNxqCtCVFkKipnSbX2tobUVQq1pTCipz5BSI8eG6BWlyM9fUtlmkhkgz59DkIkonmsSNthbJ+Q5zjnJ5qgPISQuhQWAn/88SDoiI+vOe+HmtUVZDJg0CBgxgxg8GAhALG2FkbFHDsmBB/PPCPsb8ThsDvS02HDGP6Xa3g0iYQxiBlDWU3r0piJtk9IRVYFSq+WglkwWHW1gtRVCm4rvJb6BCEikQS9en0NS0t3pKa+DQDw8voXLCzamb7whDRjdQYhnHM1Y+xDAI82QXkIIfrUauD8eeDAAaGp5MwZQFW9M6NBFhZCkBEcLDzWZ2UF3L37UMNtjVWqUuH7rCwEt2uH3/MNT7Z8OCAAxSoVxl+qvqS9OalKhPutzFGiOLEYVt2tILIQQeLGgHZCN7n69AkBAMYYunZdjby8P1BQcApubi+ZvNyENHfGNjz+lzH2DKPlHAlpGsXFwL59woJw7u7CtOmRkcDJk8YHIACgVAq1HYYCEK0m+rPen52NApUKKWVl1fb529jAijH8kJWFimZWCwIABX9p+uNzIOdgDlQFKqhKVRC7PRhSXJ+aEC2VqhRFRRfg7DyFakGaAcZY0Lx589y1z1euXOkSERHhZsyx+/fvl/n4+PhqN0tLy8AdO3a0B4CkpCRpnz59fDw9Pf3Hjx/fraysrNof3c6dO+3ffPPNajOymlO/fv18GvsaxgYhEQC+A1DOGCtgjBUyxmiUDCGmlJIiTIM+dqywTsukScKU6Wn1/HBjDAgKEjqf/vqr0AG1Gdieng4HsRgp5dVnGF3bvTvGdeiAPVlZKDeis2pTyvoxC6nvpFY+5+UcFWkV4EoOkcuDf4MNCULu398LlSofnTqFmaSs5OFIpVIeExPjkJaWVq9JygAgJCSkUDuL6rFjx5KtrKzUkyZNKgCAiIgI9wULFmSkpqZetre3V65fv76j/vGzZs3KX7NmTXr1M9ePQmG6GeHPnz+fZLKT1cDYG20PYBaArpzz1YwxDwCujVcsQtoApVJYnfbAAWEK9MsPsYxS797AsGHA8OHAkCGAY/MavJapmabdzkB/k8ft7THO0REFSiV+vH8fx/Ka14zWaV/VElw45gJFgETlZTAIyc09AmvrHrCy8jR4eHr6FlhaeqJ9+6EmKm3rEXH9uttHd+6Y7HNmkbt7WlSPHvdqyyMWi3loaGjWmjVrXDZu3Hi3odfasWOHw5AhQ/JlMplarVbj5MmTsn379t0EgDlz5mRHRka6LVu2rMrwqA0bNnSIjY213b59+61nnnnGSyaTqeLj422zsrIkb7/99p2wsLBcAHjrrbdcvvvuuw6MMQwfPjw/Ojr6bnBwcK/g4OCi06dP240bNy7v5Zdfzg4LC/O8e/euFACioqJujRo1qvjo0aM2ERERHmVlZSIrKyv11q1b/w4ICCiPjY21CgsL66pQKJharcaePXtuyOXychsbm34lJSXnDxw4IFu9erWbo6OjIjk52Voul5fs3bv3b5FIhN27d9u/8cYb7o6Ojkq5XF6SmppqefTo0evG3itjg5BPAagB/H97dx7eVJn+j//9ZG3SpEu67y1QKC2l0CI7w8ggoAKKioALAgIC41yMjKgM81EGEWFGUNEPX0UHBT+iuA3OKC514aeCFtkKFFqW2kKhTbd0ydYs5/n9kbS0Tfc2TZf7dV25enpyzsn9FJreuc+zTAGwEUA1HKvj3tTWFyKEwDEx2FdfORKPL75wDJntiJgYR8IxZYrjEdazPxO8X1wMAUBVE7eSno+LA2MMMwMCIGcMHxQXd3+ALTBfcb19VIv7OpIQmX4ojJJ0cM5Re9e6quoozpyZDYnEF8OHfwGVKqXhdc1XodN9g5iY/6EhuT3I2rVri5OTk5M2bNjQoCrR1gXsAMf6MatXr9YCgFarlajVarvUucp0bGysRavVtjoESqvVSo8dO5Z96tQprzlz5gxavHix7oMPPvD5/PPP/Y8fP56tVqsFrVZbl9VXVFSIf/311xwAmDVrVtyaNWu006dP11+8eFE2ffr0+Nzc3KyUlBTz0aNHs6VSKQ4cOKB+4oknIr/66qvLr7zyStCqVau0K1euLDebzczWxNw958+fV5w6dSo3NjbWmpaWlpCenq6aNGmSYfXq1TGHDh3KTkhIsMyaNavdw+ramoSM4ZynMsZOAgDnXMcYo3FkhLTFxYuO/h2ffQb89FP7+nTUksmAm28GZs503K4ZMKDb+nJ0hbeLiiBhzGVekNs1Gkz08wMAqCUSzNBo8Gkzi9h5ir26+X+v2uG54rJYcF8bBMEIsdgbRuOlusnILBYjTp6chGHD/g1//z/UnavV7gXAERq6yM0tIO2h0WiEuXPnlm3ZsiVYoVDUdVBauXJl+cqVK1v91JCfny/NyclR3HXXXVWAY92hxhhjrd5znD17doVYLEZaWpq5rKxMCgDp6ek+DzzwQKlarRYAx0J7tccvWLCgLrbDhw/7XLx4sW6tA71eL9bpdKLy8nLxvHnz4vLy8rwYY9xqtTIAGDdunOGFF14IKygokM2fP1+XnJzscs80OTnZMHDgQCsAJCUlGS9fvixTq9X2qKiomoSEBAsAzJ8/v/zNN98Maq1t9bU1CbEyxsRwzA0CxlgQHJUR4mnnz99YWt0TEhIcfyCb01R8HTmnO3V1m3buBF56qf1xhIY6VrGdOROYOhXw8CJuHXXeYMBJvd5lPwOweUDDqYbuCQrqeUmI3jUJGfTqIIi8RLBaiwGLDLwkABgA2O0G2O0GnD49A1brjTXK7PZqnD59K5KTP4dGcws45ygqeht+fr+HQuGeOVlIx61bt06bmpqaOH/+/Lp/xLZWQvbu3es/Y8aMCrlczgEgNDTUVl1dLbZarZBKpcjLy5MFBwe3+gbn5eVVl6jUJjL1K22N1SYmtccdO3bsvEqlapDsLF26NHry5MnV6enpl3NycmRTpkwZAgArVqwonzRpkuHf//6376233jp4586debNnz26wkFNtewBALBbDZrOxphKs9mprErIDwL8BBDPGngNwD4C/dfrVSef89BMwZ45jsipPyc1teWKr225zdLjs7DndqavbNHNm25OQtDTH8TNnAqmpQA+YLbSzXr/e9G34BcHBGN4osZoVGAgpY7D2oM6pTVVCwh4Og0gqgsVSDGbwh+D4oAq73QBBMDY5Hbu3dxJ8fMYAACorf4LJdAkxMfQ22hOFhITYZ82apdu3b1/gggULyoC2V0I++ugjzaZNm+r6k4hEIowdO7b6rbfe8l++fLlu9+7dATNnzuxQx6cZM2ZUPffcc+HLli0rr70dU78aUmvixIlVW7duDX722We1AHDkyBHF+PHjTVVVVeLIyEgLALz++ut1nWPPnTsnGzp0aE1SUlJxbm6u/NSpU4rGSUhTUlJSzFevXpXn5OTIhgwZYtm/f3+7O6O1KQnhnL/LGDsO4A9wfIC5k3NOKwJ62htveDYBAYAHHwR8fByrqzb1aOpT7c6dQEhI8+cYXd/Ae63CQuDECUcy0dTQU6USmDbNkXTcdluP79vRXgLn2F3k2uFfwhg2NpG0+UokmObvj8872lemiwk1Ari1UUIkAkRyR3JotRZDVBMAW4njrdRuN0ClGobU1J9x+vStMBodb5NyeQySkw/WDcMtKnobYrEKQUH3dF9jepntgwZdb60jqTutX7++aM+ePe26tZCTkyMrLCyU3XbbbQ3+gG/btq1g3rx5Azdt2hSRlJRkXL16dYfeuO+5556qEydOKEeMGDFUKpXyqVOnVr766qsuHWh37dp1denSpdGDBw9OtNvtbMyYMdXjx4+/8uSTTxYtXbo0bseOHaGTJk2qG9r1zjvvaD788MMAiUTCg4KCrM8//3ybfu4qlYpv3749f8aMGfEajcY2cuRIQ3vb1CXlFHcbNWoUP3bsmKfD6HlKShyflgsKPBfDyJGOGTat1qYfpaVN//HtyW65BfDza5gYyWQ3tt94A2g82dYzzzhun9QeYzIBn3wCfPONo/3+/jdWq42JAWbNciQekyc7Jgzroz7UanHvedfPKyvDw7Gz8eJ5Tod0Ojx/5Qq+dv68+O9/784QW2QpteBI0JEG+yR+EkzUTQQAHDuWBmuOGrb/uw32J59EauovddUOq7UcZ8/eAYPhHEaOPAxvb8eUC3a7AUeOhCIo6F4kJPzLrfEzxo5zzke59UW6SGZmZl5KSoqHP1WR9qqsrBT5+voKgiBg4cKF0fHx8eZnnnmmQe/yzMzMwJSUlNimzm/3WGjSgwQFtTwJVXf4+OOWb13Exbneujh/3jEBl9UKWCyuicu0aUAzJfxuceUKcPVq84lVU5Wav//ddV9sLPDXvwL33+/4GWRmOhKPxMRe1am0M/7WxG01hUiE/4lpesgqAPze3x+j1Gqof/rJjZG1TVP9QcTqG8OMrdZiSERxqNHeuB1TSyrVYPjwdJhMl+oSEAAoKfkIdrueOqSSPuGll14KfO+99wKtVitLSkoyrlmzpl2JJCUhvV1CgmNdEE9xDjtrVlPxqVQtd7IcPtwxWZenHDzoSJKac+utjqnOa6uInAPvvQcEBt5IVICGI1gSEoAZM9wbdw9zqroaF0wml/1/joxEmFze4rniHpKkNdUfpDYJ4ZzDYimGUhYMg8lRzbLbGyaoYrEXVKphDfYVFb0NhWIQfH0nuilqQrrPM888U9y48tEelIT0dl984ekIWtaR+Ppim/qhFRcuuOzzk0iwNiqq1XN7dBKiciQhBkMWOLdArgwDzI4kRBBaviVuMuWiouIQ4uI2NTvKgZD+hJIQQkiXO6vXI6PatXP9U9HR8G+tegZHx9WeoKXbMRcvrgQA6JR7gNnJAACrtenVgWsVFe0BwBASsrBrAyWkl6IkhBDS5QKlUogA+EskKHPOvhgmk+FPERFtOl/UU5KQJiohErUEgmBDdfVxAEANvwzMuQwAMJt/a/ZanAsoKtoDf/+p8PJqvRpESH/Q+ychIIT0OPtLSiAA+G7ECLydkIAYuRxPx8RA2cTaMT2Zrdp1+mqxWgy9/gQEwbW/i1js2+y1KioOoaYmnxarI6QeSkIIIV1ub1ERUlUqDFep8FBoKHLGjMHDvXAOlOb6hFRU/H+uB3OAc9fEpFZR0VsQi30RGHhnV4ZIuhBjLG3ZsmV1vdKffvrpkDVr1oS39fwVK1ZEDho0KGnAgAFJixYtihKc0xP8+OOPysGDBydGR0cPq7+fUBJCCOli5wwGnNDr8WBISN0+uUgEaQdmf/Xy8IyxzfUJaTIJsUsaDNGtz2arQknJxwgOng+xWNHkMcTzZDIZP3jwoH9hYWG7uyqkp6d7Hz16VJWdnZ114cKFrFOnTnkfPHhQDQCrVq2K2blzZ35eXt7Z3Nxcr48++sin66PvnSgJIYR0qXe0WogBLKiXhHTEp8OGIesmzy7U3VQlRKQGKit/dD24RtlsElJc/AEEwYSwMLoV05OJxWK+cOHCks2bN7f7Py9jDDU1NcxsNjOTySSy2WwsPDzcmp+fL9Xr9aKpU6caRCIR7r///rIDBw74uyP+3shtHVMZY1EA9gIIhWOxu12c85frPf84gH8CCOKc0yx5hPQBAuf4P60W0zQahLS0CGAbzA4MbP0gN2sqCbEH5cBur2qwj1m8wQ3ezSYhRUVvQakcCrV6tFvi7Guyl2RHGc4alF15Te9h3saE3QlXWztu7dq1xcnJyUkbNmxosN5AawvYTZ061TBhwoTqsLCwFABYtGhRSWpqqvmHH35QhoWF1S1YFxMTYyksLGx9iFg/4c7RMTYAf+Gcn2CMqQEcZ4ylc87POROUWwBccePrE0K6iZ1ziAAcqqhAQU0N/tloddzeqqnbMTUBR132ySvSYDZcbTIJMRpzUFV1BAMGbKW5QXoBjUYjzJ07t2zLli3BCoWirvNGawvYnT17Vn7hwgWvgoKC0wAwefLkwV988YVKqVS6dACh/wc3uC0J4ZwXAih0blczxs4DiABwDsCLAJ4A8Km7Xp8Q0n025+cjXaeDSiyGWizGHT2gitEVmhodY/D5xmWf0jYOZlMJ7DV6l+ccc4OIERLyoDtC7JPaUrFwp3Xr1mlTU1MT58+fX1elb60Ssn//fr+bbrrJ4OvrKwDA1KlTKw8fPuy9bNmy8vqVj/z8fFloaKi18XX6q27pE8IYiwUwEkAGY2w2gGuc88xWzlnOGDvGGDtWUlLSDVESQjoivbwcG/Ly8GNlJb4oL4efRIKLTUzX3hu53I7xqYRR+ovLcb6yqYDZC7ZGSQjndhQV7YVGMwNyee8bHdRfhYSE2GfNmqXbt29fXTa9cuXK8uzs7HONH19++WUuAERHR1sOHz6stlqtqKmpYYcPH1YnJiaaY2JirN7e3sK3337rLQgC3n333YA77rijwnOt61ncnoQwxlQAPgbwZzhu0awH8HRr53HOd3HOR3HORwUFtWs1ZUJIN8k3m7Hg3DnUrzdfralB6rFjKDCbPRZXV3FJQsb+ArCG++TyGPgEpgEmBWyWhklIeXk6LJZr1CG1F1q/fn1RRUVFm+8WLF68WBcbG1szZMiQpMTExMSkpCTjfffdVwkAO3fuzF+xYkVsTEzMsNjY2Jq5c+dWtna9/sKtM6YyxqRwJCDvcs4/YYwlA4gDkOm8JxYJ4ARjbDTnvKiFSxFCeqBn8/LqZkSt7/GoKER6eXkgoq7l0idkwmGXYwID74RckANmL9hthQ2eKyp6GxKJBgEBM90ZJukiRqPxZO12VFSUzWQynWzp+PokEgn27duX39Rzv/vd74wXL17M6ooY+xp3jo5hAP4F4DznfDsAcM7PAAiud0wegFE0OoaQ3unV+HjobDZ8UnrjV3iKnx82xcV5MKqu06ASIjcDo107pQYG3gmZWAaYvSDwGx1TrVYdSksPIDx8OUSillcNJqS/cuftmAkAHgQwhTF2yvm4zY2vRwjpZmLGkG82w1skgowxRMnleD8xERIPTzLWVRp0TE07DnjVNHheItHA13cixGoxAm6JcCQqTsXF74HzGpqmnZAWuHN0zE8AWhyHxDmPddfrE0Lcb1N+Po7r9fgoKQlRcjlEAII6OT9IT9LgdszEn1yeDwiYBZHI8TbqPSgQZVcM4JyDMYaiorfh7T0cKtWI7gqXkF6HVtElhHTI0aoqPJefjwdDQnB3H+w8LlgF8Bru+EZkB8YfcTmm/jowYrE3ADs4t8BovITq6l8xcOCLNCcEIS2gJISQDrIIAsqsVpRZrSi1WlFms6G0dttqRbhMhlURERAzBhljLn+Mii0WLMrObrAvSCrFnqFDm33Nps7pTrXxGe12PHj+PMLlcrwSH++xeNypQX+QYWcB34azpIpECmg00+p97+04z25AYeFbYEyCkJD7uyVWQnorSkIIAVDTOKFwfq2fXNTuK7ZYUGazQW93nU2zPgZgbW4uAEAMwFsshlIshrdIBG+xGBLGcELfcEinWixGVG4uvMXiuofSeby3WIxqmw1flDc7aaPbRckdHSyfuHwZF0wmfJuSAl9J33wbaZCENHErxt9/GsTiGzOLOyohgM1WCa32/xAQMBMyWd+rEBHSlfrmuwchAAx2OzKqqhomFM1ULapbSSg6wkcsxlPR0bA7Y6l9GAUBBrsdZVbXSRP1djs2X7kC3uXRdJ2vy8vxv9ev48+RkZji33fX4brRH4Q3OzS3vtokpKTkY1itWuqQ2gsxxtKWLl2qfeONNwoA4Omnnw7R6/Xi7du3X2/L+StWrIj85ptvfAVBwO9+97uq3bt3XxWJRPjxxx+VDz/8cKzZbBZNmTKlsnZ/ff/4xz+ClEql8Oijj5a5oWntlpeXJ12xYkVU7WRs7kJJCOkVjHY7fq6qapBIyBjDUzExTR5vstvxS1UVpma2ODGvW7EW4gOAK2YzYn5pOPtmpFyO/LFjYXYmKgbnV6Nz+zeTCYtzctwderMKLRbck5WFoUolNveRYbjNsRvs8J/hj4orJ8DDCxs9K3KZ+6M2Cbl+fSek0mBoNLd2U6Skq8hkMn7w4EH/wsLCorCwMNcJcFqQnp7uffToUVV2dnYWAIwaNSrh4MGD6pkzZ1avWrUqZufOnflTpkwx/P73v4//6KOPfO69994G9/eeeOKJLpka3GazQdIF1cnY2FiruxMQgJIQ0o1Mzk//TVUiAODp2Nhmz9VaLC4Jha9YDLMgQGu1QmuxoNhiqdt2R2WjvSpsNlgFAdJ2DldljEEhFkMhFqPxCixxHp4ATMYYbtVo8ExsLBRisUdj6SwucFQcqkDhvwoR8ccI+I73hWAVUPppKa7/v+uo+K4CTMKgeO44jI3O9fWdBJms4b9ObZ8Qs/k3REaugUhEC6V2xqU1l8ILXizosrnuIx+LLBy0fVCLFQ2xWMwXLlxYsnnz5pBXXnnlWnuuzxhDTU0NM5vNjHPObDYbCw8Pt+bn50v1er1o6tSpBgC4//77yw4cOODfOAlZs2ZNuEqlsm/cuFE7evToIWlpafqffvrJp7q6Wvzaa6/lzZgxQ2+z2bBq1arIQ4cO+QDAQw89VLp+/friiIiI5AULFpR+//33Po888kjxhAkTjCtWrIguLy+XeHl5CW+++Wb+yJEjzfv27fPdsmVLmNVqFfn7+9v279+fGxUVZfv8889Vf/nLX6Jr23HkyJHs4uJiycyZM+MvXryYtWPHjoDPPvvMz2Qyia5cuSK/9dZbK1577bUCAHjxxRcDX3755dDg4GDrgAEDzDKZjO/du7fNi9NSEkI6xFg/oWih/0T9baPgsphkHbVYjHiFAlpnIlFssTTYLqqpcTmn0m7H3/PzESCRIEQmQ4hMhjSVqm47UCrFIxcuuO1nwAAESKUIkEgQKJUiUCpFQKOvLd1WCZJK8XlycoN9Xq0kLE2d0528RKJefwum5noNCncXomh3Ecy/Oeb14FaO8i/KUfhmISxFFsij5Yh7Lg6hS0JxpuAxoNG6dI1vxQA3KiEA6FZML7Z27dri5OTkpA0bNjSYxbu1BeymTp1qmDBhQnVYWFgKACxatKgkNTXV/MMPPyjDwsLq7r3GxMRY6i9o1xybzcbOnDlzfv/+/b4bN24MnzFjxoVt27YF5efny7Oyss5JpVJotdq6TwJeXl7C8ePHcwBg3Lhxg3ft2pWfnJxc891333mvXLky+pdffrlwyy236OfPn58tEomwffv2wI0bN4a+8cYbBdu2bQvdsWNH/rRp0wyVlZUipVIpFBcXN4jn3LlzyszMzHMKhUIYNGjQsMcff1wrkUjwwgsvhJ04ceKcn5+fMH78+MFJSUntWjiKkhACwLEUe6Ze36b+E6VWK0wtJBQdUW23477z5wEAEsYQLJUiWCZDiFSKRKUSQVIpXr52DTbe8M+6bsIE+Emb/31+KjcXuiamFW9MBGdCUZtENJFYNE4y/CQSiDsx/FIhFuO2gAC3n0Maqvi+Ann/k9dgX8mHjkp4wMwAhK8Ih2aGBkzMYDbnQ5/tOnN3YOAdLvtqkxC1ehRUqmFdHzjpFhqNRpg7d27Zli1bghUKRd0b3cqVK8tXrlzZbK/ws2fPyi9cuOBVUFBwGgAmT548+IsvvlAplUqXN8u2DNueO3euDgDGjx9vWLt2rQwAvvvuO58VK1aUSJ3veSEhIXUl34ULF+oAoLKyUnTy5EnV3LlzB9Y+Z7FYGAD89ttvsjvvvDOypKREarFYRFFRUTUAMHbsWP3jjz8ede+995YvWLBAN3DgQJeYJ06cWBUQEGAHgEGDBpkvX74sLy4ulowZM6a6No45c+boLly40K5yLSUhBAAgcI6048c9GsOPI0ZgqLc3/CUSiJr4JX2vuBjXLZYG+ypsthaTkBkaDUyCcCOJcCYXjRMLv2Zek/Q9PhN8IPISQTA3fJ+Nez4OMU817MNTWvqpy/ne3ilQKFz7w0ilwQBECAtb3qXxku63bt06bWpqauL8+fPr1iNorRKyf/9+v5tuusng6+srAMDUqVMrDx8+7L1s2bLy+pWP/Px8WWhoqGuv9Ea8vLw44FiTxm63MwC1E+E1WWBVq9UCANjtdqjValt2dva5xsc8+uij0atXry66//77Kz/77DP1xo0bwwFg8+bNRXfeeWflp59+6jt+/PihX3755YXGyZNMJqt7XbFYzK1WK+O8813oKQkhAIB0nQ5eIhHMXVzhaI/aSkRzpms0qLDZGiQR6lY6Cu2ZDQAAIABJREFUYO1LTOzqMEkvVfVrFa5uu4rST0rBra5vnhWHKppIQg64HNfUrRgAkMtDMWbMZXh5Nd8ZmfQOISEh9lmzZun27dsXuGDBgjKg9UpIdHS05a233gqyWq2FgiCww4cPq//0pz9pY2JirN7e3sK3337rffPNNxvefffdgD/+8Y/FzV2nJVOnTq167bXXgm6//fbq2tsx9ashgKOSExkZadm9e7f/kiVLdIIgICMjQzFu3DhTdXW1ODo62goAb7/9dl1JNSsrSz569GjT6NGjTRkZGd5nz571Gj16dOOuUC4mTZpkWLduXVRJSYnYz8/P/umnn/oPHTqUbseQ9jtQWur2BETCWIv9J1pKQABgd0KCW+MjfUPtp7PakrelxILcdbko+lcRJP4SRDwaAd+Jvsi6OwtggP80f4Q9HIbA2Q07mlqtZaio+MHl+s0lIQCgUMR2XUP6uUHbB11vrSOpO61fv75oz549bZ7oZfHixbrvv//eZ8iQIUmMMdx8882V9913XyUA7Ny5M//hhx+OM5vN7Oabb66aO3duZUdieuyxx0ouXLggT0hISJJIJPyhhx4q+etf/+oyqua9997LXbZsWczWrVvDbDYbmzNnTvm4ceNM69evv75gwYKBISEhllGjRhmuXLkiB4B//OMfwUeOHPERiUR88ODBpnvuuafyypUrrfZbiYuLsz722GOFN91009Dg4GDr4MGDTb6+vu0aFdAl5RR3GzVqFD927Jinw+jTjHY7bj51Ckerq9t0vJSxdvWfCJRK4SMW0xTWxK0MWQZcWHkBYcvDELIgBNdfv47f1v8Gu96OyD9HIubpGEjUjs9e13ZeQ8DMAHhFN30Lu6hoL7KzH2qwTy6Pwdixv/Wa/8eMseOc81GejqMtMjMz81JSUmhF9V6msrJS5OvrK1itVkyfPn3QokWLShcuXFhR/5jMzMzAlJSU2KbOp0oIAQAoxWKM9vGBUixutf9EoFQKNSUUpAexG+zIezYPBdsKwG0chrMGXP3HVRjOGOA3xQ/xr8TDO9G7wTkRqyJavGZzt2Lo/z0hN6xduzb8hx9+8KmpqWGTJ0+ueuCBBypaP+sGSkJInb66Bgjp20r/W4qLf7qImvwbw7htOhsEi4DE9xMRdG9QuxMHu92I8vIvXfa3dCuGkP5o165dBZ05n5IQQkivZL5ixqXVl1B6oOkKvmAU4BXr1aHKhU73DQShYf86iUQDX9+JHYqVENK09k3lSAghHiZYBVx54QqODj3abAICAIpBCnBbx/q8eXnFICzsETB2o29eQMAsiET0uY2QrkS/UYSQXkP3rQ6XHrsEwxlDs8cwOUPMuhhEPRkFsVfHppZXqVIQG/s0CgtfR1jYI5BI/KDRTOto2ISQZlASQgjp8QxZBlx+4jLKDzY7TQMAwP8Wf8T/bzyU8cpOv2ZZ2X8BAJGRf4K3d1Knr0cIcUW3YwghPVZNUQ1ylufg1+G/tpiAyEJlSHw/EcO/Gt4lCQjgGB3j5TUQSiVNeNdfMMbSli1bFln7/dNPPx2yZs2a8Laev3Llyoj4+Pik+Pj4pDfeeKNukaXs7GzZ8OHDE2JiYobdfvvtA8xmMw2xcqIkhBDS49gNduRtzEPGoAwUvlEINDePngiI+FMERmePRvC84C4bPmu1lkGn+5aG5PYzMpmMHzx40L+wsLDddwnef/9938zMTOW5c+eyjh8/fv7ll18OLS8vFwHAmjVrIh999FFtfn7+WV9fX9vLL7/ceIHsfouSEEJIj8HtHIX/KkRGfAbynsmDYGhh5eXRaqQdTUP8jnhIfLv2znJh4Zvg3IrQ0EVdel3Ss4nFYr5w4cKSzZs3h7T33KysLK+JEyfqpVIpfHx8hMTEROMnn3ziKwgCfv75Z/XixYt1ALBkyZKy//73v35dH33vRH1CCCE9QvlX5bi89nKLnU4BQB4jx4AtAxB8bzCYqOurFIJgw7VrO+HndzOthush2dlLogyGs11zX83J23uYMSFh99XWjlu7dm1xcnJy0oYNG4rq729tAbuRI0eaNm3aFF5dXa3V6/WiI0eO+AwdOtSs1WolarXaXrvybWxsrEWr1cq6rGG9HCUhhBCP0p/W4/Lay9B9rWvxOLGvGDF/i0HEoxEdHvXSEpMpF4JggdF4HjU1VzBo0Etd/hqk59NoNMLcuXPLtmzZEqxQKOpKca0tYHfXXXdVZWRkKG+66aYEjUZjTU1N1UskEt7U0ijNrYTbH1ESQgjxCFOuCfnP5aPorSKghbdkJmUIXxWO2P+JhTSg1TW1Oiw//1kUFb0NicQfEkkgNJpb3fZapGVtqVi407p167SpqamJ8+fPr5uIprVKCABs3bq1aOvWrUUAMGvWrLjBgwfXhIaG2qqrq8VWqxVSqRR5eXmy4OBga/e1pmejJIQQ0m2sOitKPiyB9h0tKn9qfSHRoHuCEPd8HJSDurQy7xqXtQzFxe8DAGw2R0Xm6NHBSE7+L1SqFLe+Nul5QkJC7LNmzdLt27cvcMGCBWVA65UQm82G0tJScWhoqD0jI0ORnZ2tvOuuu34TiUQYO3Zs9VtvveW/fPly3e7duwNmzpzZrvVV+jJKQgghblfxYwWu7biG0v+Wgte0Xon2GeuDgdsGwne8bzdEBxQW7oYgmBvss9v1UChoPaX+av369UV79uwJauvxFouFTZgwIQEAVCqVfc+ePbm1/UC2bdtWMG/evIGbNm2KSEpKMq5evZpWC3ZyWxLCGIsCsBdAKBwD7HZxzl9mjD0L4A7nvmIAizjn190VByHE8wynDSj5qKTV47zivDBg6wAE3dP+Rec6inM7rl//fy77Q0MXQyx2bwWG9CxGo/Fk7XZUVJTNZDKdbOn4+pRKJb98+XJWU88lJiZazpw5c74rYuxr3DlE1wbgL5zzoQDGAvgjYywRwD8558M55yMAfAbgaTfGQAjxMFuVDXajvcVjJP4SDNw+EKPPj0bw3K6b76MtOLcjMnI1pNKGUzeEh6/sthgI6a/cVgnhnBcCKHRuVzPGzgOI4Jyfq3eYN1rskkYI6emsFVbo0nUInhvcYH/VsSoUvl4I7XtaCAYBYl8x7JUNkxHlUCVCFoYgfHk4pBr3dTptiUgkQ0TEo7h69SUolSFQKAaCczuUykEeiYeQ/qRb+oQwxmIBjASQ4fz+OQALAVQCuLmZc5YDWA4A0dHR3REmIaSNBKuA8i/Lod2rrevn4Z3lDXm0HMX7inH99evQn9BDpBQheH4wwh8JhznfjHP3noM0WIrgBcEIfTAUqlRVj5iRtKzsIGpq8pCY+AGCg+dCEGyeDomQfsHtSQhjTAXgYwB/5pxXAQDnfD2A9YyxdQAeBfBM4/M457sA7AKAUaNGUbWEkB7kzO1noEtvOK/HufvPwXzJDLveDu9kb8S/Go+QB0LqZjP1Hu6N5M+S4T/NHyJpz5qs+dq1VyCTRSAw8E4AgEhEffYJ6Q5ufSdgjEnhSEDe5Zx/0sQh+wDc7c4YCCFdL+D2AJd9hlMGBNwVgJFHRmJU5ihE/DGiwXTqYi8xAm4P6HEJiMGQDZ0uHRERKyESeeaWECH9ldveDZijxvovAOc559vr7a8/5m02gGx3xUAI6RguNF981J/Rozqzusnnwh4Kg+843x5xi6Wtrl17FYzJEBa2zNOhENLvuLPmOAHAgwDOMMZOOff9FcDDjLEhcAzRzQewwo0xENLvcc5hr7bDWmKFtdQKS4kF1lJr3feNv1pKLFAmKJH2S1rdNewmO0o+KMH116+j6ucqMDmDLFwGy3ULAEAa5OjnIY+Se6qZHWKzVUGr3YPg4PmQyYJbP4H0aYyxtKVLl2rfeOONAgB4+umnQ/R6vXj79u1tmkZi5cqVEd98840fADzxxBPXly1bpgOA7Oxs2b333jugsrJSMmzYMOPHH3/8m5eXV4NM/9133/XNyspSbN68uaipa3vCyJEjE06ePOnWQoE7R8f8BKCpj0MH3fWahPQHglVwJA3NJBBN7efW9nWrqs6ohjHHCG7nuP76dWj3amGrsEExWIGB2wYi9KFQVPxQgeL3ixG6MLRH9vNoidF4EZzbodN9Dbtdj4iIP3k6JNIDyGQyfvDgQf/CwsKisLCwdvVOfv/9930zMzOV586dyzKZTKLx48cPufvuuys1Go2wZs2ayEcffVS7fPly3X333Rf98ssvBz755JMNJs65//77K+EYrNEptdPDdwV3JyAAzZhKSI9V/EExyr8sd0kyGg9zdZfMWzJRc7UGTMoQdHcQwh4Jg99kv7pbLUFzghA0p80TSvYoeXl/R3HxuxCJFFAoBkOlGuHpkEgjly6tCS8oeDGsq64XGflY4aBBLVc0xGIxX7hwYcnmzZtDXnnllWvtuX5WVpbXxIkT9VKpFFKpVEhMTDR+8sknvkuWLNH9/PPP6k8//TQXAJYsWVK2YcOG8MZJyI4dOwKOHTvmvXfv3it33313rFqttmdmZnqXlJRIn3322YLFixfrAOBvf/tbyAcffBDAGMMf/vCHyp07d14bPXr0kNGjR+szMjJUt912W8UjjzxStnjx4phr167JAGD79u1Xpk2bZvj++++Va9asiTabzSIvLy/h7bff/i0lJaXm2LFjXosXL46zWq1MEAR8/PHHl5OTk2uUSuVIo9F48rPPPlNv3LgxXKPRWHNychTJycnGAwcO/CYSibB//37fp556KlKj0diSk5ON+fn58u+///5SW39ulIQQ0g3KPi+rSyhqqxWRf45E2KLm32OrMqoci7t5iKXQgrjn4xC2JAyy4L6z8rjFokVJyQcAAEEwwWS6gIyMAUhOPgiVapiHoyOetnbt2uLk5OSkDRs2NPjla20Bu5EjR5o2bdoUXl1drdXr9aIjR474DB061KzVaiVqtdpeW52IjY21aLXaVn+htFqt9NixY9mnTp3ymjNnzqDFixfrPvjgA5/PP//c//jx49lqtVrQarV1y0lXVFSIf/311xzAsXjemjVrtNOnT9dfvHhRNn369Pjc3NyslJQU89GjR7OlUikOHDigfuKJJyK/+uqry6+88krQqlWrtCtXriw3m83MZnMtAp0/f15x6tSp3NjYWGtaWlpCenq6atKkSYbVq1fHHDp0KDshIcEya9asuPb+vCkJIaQZ7elLEbokFOHLwpu9VuXPlbj2asMPVoYsA4yXjLBX2mGrsN14VDq+Vv1c5e4mNo8BsnAZgu8N7lMJCAAUFr4JzhsuYioIJigUNDkZATQajTB37tyyLVu2BCsUCqF2f2sL2N11111VGRkZyptuuilBo9FYU1NT9RKJhHPueiuUMdbq/dHZs2dXiMVipKWlmcvKyqQAkJ6e7vPAAw+UqtVqAXAstFd7/IIFC+piO3z4sM/FixcVtd/r9XqxTqcTlZeXi+fNmxeXl5fnxRjjVquVAcC4ceMML7zwQlhBQYFs/vz5uuTk5JrG8SQnJxsGDhxoBYCkpCTj5cuXZWq12h4VFVWTkJBgAYD58+eXv/nmm+0qj1ISQvoNwSrAWtZEX4rmkoxSK7ilbX0pZBEySHwkLolE7cNw3uByTsELBSh4oaDpCzJA5NV1fSxEXiJIg6SOR+CNr7IgWYPv675qpGDi3jPCpa0EwYbr119z2R8WthRisZcHIiI90bp167SpqamJ8+fPr1torrVKCABs3bq1aOvWrUWAoxoxePDgmtDQUFt1dbW4tq9GXl6eLDg42Nr4Oo3V77ham8hwzpsdeVabmNQed+zYsfMqlarBG9jSpUujJ0+eXJ2enn45JydHNmXKlCEAsGLFivJJkyYZ/v3vf/veeuutg3fu3Jk3e/bsBkPg5HJ53bXEYjFsNhtrKsFqL0pCSJ9gvGhE1ZGqFqsVtgr3zYJZ+nEpSj++sTAmkzBI/CR1D7G32OUc1QgVItdEOo7xlTQ8XiVG2X/LcPbOs64vxgCJRuKaQDSRTNQe09Tr90dlZf9FTU3jxE+E8HAapEduCAkJsc+aNUu3b9++wAULFpQBrVdCbDYbSktLxaGhofaMjAxFdna28q677vpNJBJh7Nix1W+99Zb/8uXLdbt37w6YOXNmRUfimjFjRtVzzz0XvmzZsvLa2zH1qyG1Jk6cWLV169bgZ599VgsAR44cUYwfP95UVVUljoyMtADA66+/XrdY0rlz52RDhw6tSUpKKs7NzZWfOnVK0TgJaUpKSor56tWr8pycHNmQIUMs+/fv17S3TZSEkF5PsAm4uv0qCl8r9FgMPuN9MOTNIXXJhEghavCJpfzrcpyefrrBORJ/CUIfdPlgVUeVqsLg1wb3myqFu5jNV6DTpaO8/GvodF+7PB8QMBNeXjEeiIy0ZNCg7ddb60jqTuvXry/as2dPm28tWCwWNmHChAQAUKlU9j179uTW9gPZtm1bwbx58wZu2rQpIikpybh69erSFi/WjHvuuafqxIkTyhEjRgyVSqV86tSpla+++qpLB9pdu3ZdXbp0afTgwYMT7XY7GzNmTPX48eOvPPnkk0VLly6N27FjR+ikSZPq7ve+8847mg8//DBAIpHwoKAg6/PPP9+mn7tKpeLbt2/PnzFjRrxGo7GNHDnSteTbii4pp7jbqFGj+LFjxzwdBumhBJuAH6Q/eDQG9Wg10jLSmn2+5loNivcXN0go5GFyyCN617wavYHNVoWKikPOpCMdJtOFFo8fPvwraDTTuim67sUYO845H+XpONoiMzMzLyUlpUN/nInnVFZWinx9fQVBELBw4cLo+Ph48zPPPFNc/5jMzMzAlJSU2KbOp0oI6fVEEhGG/GsIch7O6fpr1/alaOFWhzRICllYy5035RFyRK2J6vL4iKOfR3X1r3XVjqqqXwC0bRizQjEY/v5T3RsgIX3YSy+9FPjee+8FWq1WlpSUZFyzZk27EklKQkif4DPOp03HSTSS1jtl1ntepBT1qinI+wPOOUymy9Dp0p2P72C3d2yOp7i458BY75lkjZCe5plnniluXPloD0pCSJ8gC5Uh+L7gFkd+SDQSiCT0B6c3slrLodN9V5d4mM2/dfhaIpECfn6TERHxKAICbu/CKEknCYIgMJFI1PP7CJA2EwSBwbFMS5MoCSF9gtRfisR3Ez0dBukigmBBVdXPKC93JB3V1cfQwvtYKxhUqlRoNNPg738LfH3HQySivjg90NmSkpLEoKCgSkpE+gZBEFhJSYkvgCaG+TlQEkII6VGMxgs4diwVgtDujvZ15PJo+PvfAo1mGvz8pkAmC2z9JOJRNpttaVFR0ZtFRUXD4MYV3km3EgCctdlsS5s7gJIQQkiPolAMhEgkhdCOwodYrIaf3811iYdCEU99eXqZtLS0YgCzPR0H6V6UhBBCuoXdbkRl5RHodOmIifkrJBJfl2PM5gKUlf0HjLV2u0QEH58x8Pe/Bf7+t8DHZwxEoq5ZOZQQ0n0oCSGEuFVBwQ5cvboNNTVX6vb5+IxDUNCd4JxDrz+FsrL/oLT0P9DrTwAApFLXOaK8vAbW9evw87sZUqlft7WBEOIelIQQQtrNbjfCZLoIozEHMlk4/PwmtnA0b5CAAEBR0b+g06WjrOw/zmnUGXx8xmPAgK0ICJgNxmQ4cSINfn5/qEs8FIp2L9BJCOnhKAkhhDSJcw6L5TqMxmwYjTnOh2O7pia/7riQkIXNJiE2WyUA13Vryso+g0ikhEYzDbGxGxEQcDtksuAGx0yYUArGaM0bQvoySkII6efqVzXqJxwmUw7sdn2r5xsMWais/Bkm02WYTJdgMl2C2ezYtlqbnzwxLe04vL0Tmn2eEhBC+j5KQgjpB25UNeonGrVVjSsAOj4tg15/HCdPjnd+xyCXR0GhGIjAwDlQKAbCy2sAzp9/EJzXQC6Pga/vePj73wK5PKxL2kYI6b0oCSGkD7HbTTCZLrjcPmlrVaOjEhL2Qq2+CV5esRCLvVyeV6mGQy6PglisdFsMhJDeh5IQQno5q1WHc+cWwGjM7nRVoz3k8mgolQlQKofA339qi5UNpXJIt8RECOldKAkhpIdyVDUuwmS6iKCgu5s9TiLxQUXFIXBe0+UxiETeUCqHOB8JdV8ViniqahBCOo2SEEJ6mDNn7oRef6pBVWP8+BLIZIHgnMNqLXF2AL1c95UxEXgnCiD1qxo3Eo0hkMsjaOZRQojbUBJCSDeorWoYjTnQaKZDIvFp9liT6VKDIbAAkJV1D+z2SphMlxr17WCQy6MhFqsgCKYWY2iqqqFQDIFSOZiqGoQQj6AkhJAu4hiBUugyr4bJlAOzOR+1VY0RI36Ej88YmM15ddWM2iGtJtNlGI05Ltc2GrOhVqfB13cSFIpBUCgGQqEYBC+vWIhEcuTmrseVK5sBUFWDENJ7uC0JYYxFAdgLIBSOlfR2cc5fZoz9E8AsABYAlwEs5pxXuCsOQrpa/apG/eGuJtMF2O3VrZ5/9uwc2Gw6APa6fSKREgrFICiVCRCLfVFd/UuDc0JDH8TAgf9s9pphYcsQFDQXSmU8xGLvDreNEEK6kzsrITYAf+Gcn2CMqQEcZ4ylA0gHsI5zbmOMbQWwDsCTboyDkHa7UdVonGg0rGp0hEwWiIiIFc5KhqOiIZOF1FUpioreQXa2IwlxVDWGQKGIb/GaCkVsh+MhhBBPcVsSwjkvBFDo3K5mjJ0HEME5/7reYb8AuMddMRDSmoZVjRuJhtGY06aqRkcoFIMQF/dss88HBNyGtLSTVNUghPR53dInhDEWC2AkgIxGTy0BsL+Zc5YDWA4A0dHRboyO9DelpZ/h+vX/hdGY3emqRnvI5VFQKhPg4zO2xeOk0gBIpQHdEhMhhHiS25MQxpgKwMcA/sw5r6q3fz0ct2zebeo8zvkuALsAYNSoUd3zV4L0CyUlH6C8/Eu3XFskUrp0CHV0EqWqBiGENObWJIQxJoUjAXmXc/5Jvf0PAZgJ4A+cd2Z2A0Lar7z8q05fo7aqcWOYq2PbMQJF1AVREkJI3+fO0TEMwL8AnOecb6+3fwYcHVEnc86N7np9QpqTmpqBjIy4Vo+jqgYhhLiXOyshEwA8COAMY+yUc99fAewAIAeQ7hwN8AvnfIUb4yCkAYUiFlJpEKzWEgBU1SCEEE9x5+iYnwA0NTPSQXe9JiFtNXToPkilAc7ZQqmqQQghnkAzppJ+SaOZ6ukQCCGk36NaMyGEEEI8gpIQQgghhHgEJSGEEEII8QhKQgghhBDiEZSEEEIIIcQjKAkhhBBCiEew3jBrOmOsBEB+F1wqEEBpF1ynp6N29i3Uzr6ju9sYwzkP6sbXI6RdekUS0lUYY8c456M8HYe7UTv7Fmpn39Ef2khIe9DtGEIIIYR4BCUhhBBCCPGI/paE7PJ0AN2E2tm3UDv7jv7QRkLarF/1CSGEEEJIz9HfKiGEEEII6SEoCSGEEEKIR/SZJIQxNoMxlsMYu8QYe6qJ56MZY98zxk4yxk4zxm5r4nk9Y+zx7ou6/TrTTsbYcMbYz4yxLMbYGcaYV/dG33YdbSdjTMoY2+Ns33nG2Lruj75t2tDGGMbYt872HWKMRdZ77iHG2EXn46Hujbx9OtpOxtiIev9fTzPG5nV/9G3XmX9P5/M+jLFrjLFXuy9qQjyMc97rHwDEAC4DGABABiATQGKjY3YBWOncTgSQ1+j5jwF8COBxT7fHHe0EIAFwGkCK8/sAAGJPt8kN7bwPwPvObSWAPACxnm5TB9v4IYCHnNtTALzj3NYAyHV+9Xdu+3u6TW5o52AA8c7tcACFAPw83aaubme9518GsA/Aq55uDz3o0V2PvlIJGQ3gEuc8l3NuAfA+gDsaHcMB+Di3fQFcr32CMXYnHG/kWd0Qa2d0pp3TAJzmnGcCAOe8jHNu74aYO6Iz7eQAvBljEgAKABYAVe4Pud3a0sZEAN86t7+v9/x0AOmc83LOuQ5AOoAZ3RBzR3S4nZzzC5zzi87t6wCKAfTU2T878+8JxlgagBAAX3dDrIT0GH0lCYkAcLXe9wXOffVtAPAAY6wAwEEAfwIAxpg3gCcB/N39YXZah9sJx6dKzhj7ijF2gjH2hLuD7YTOtPMjAAY4PjVfAfAC57zcrdF2TFvamAngbuf2HABqxlhAG8/tKTrTzjqMsdFwVBguuynOzupwOxljIgDbAKx1e5SE9DB9JQlhTexrPPZ4AYC3OeeRAG4D8I7zl//vAF7knOvdHGNX6Ew7JQAmArjf+XUOY+wP7gy2EzrTztEA7HCU7+MA/IUxNsCdwXZQW9r4OIDJjLGTACYDuAbA1sZze4rOtNNxAcbCALwDYDHnXHBXoJ3UmXauAnCQc34VhPQzEk8H0EUKAETV+z4S9W63OD0MZ8mac/6zs1NmIIAxAO5hjP0DgB8AgTFm5pz3xM5hnWlnAYD/j3NeCgCMsYMAUnGjPNyTdKad9wH4knNuBVDMGDsMYBQct9t6klbb6LwFcRcAMMZUAO7mnFc6qz+/b3TuIXcG2wkdbqfzex8AnwP4G+f8l26JuGM68+85DsAkxtgqACoAMsaYnnPu0rmVkL6mr1RCfgUQzxiLY4zJAMwH8J9Gx1wB8AcAYIwNBeAFoIRzPolzHss5jwXwEoDNPTQBATrRTgBfARjOGFM6+0tMBnCu2yJvn8608wqAKczBG8BYANndFnnbtdpGxligs7oDAOsA7HZufwVgGmPMnzHmD0d/n6+6Ke726nA7ncf/G8BezvmH3RhzR3S4nZzz+znn0c73oMfhaC8lIKRf6BNJCOfcBuBRON6IzwP4gHOexRjbyBib7TzsLwCWMcYyAbwHYBHnvKeWsJvUmXY6OzBuh+PN8hSAE5zzz7u/Fa3r5L/n/8LxafIsHG19i3N+utsb0Yo2tvH3AHIYYxfg6LT4nPPccgDPwtG+XwFs7KH9XjrVTgD3AvgdgEWMsVPOx4jubUHbdLKdhPRbNG1r/UilAAACsElEQVQ7IYQQQjyiT1RCCCGEENL7UBJCCCGEEI+gJIQQQgghHkFJCCGEEEI8gpIQQgghhHgEJSGEuBFj7Ijzayxj7D5Px0MIIT0JJSGEdJJz8rcmcc7HOzdj4ZjNlRBCiBMlIaTfYYx5M8Y+Z4xlMsbOMsbmMcbyGGNbGWNHnY9BzmNnMcYyGGMnGWPfMMZCnPs3MMZ2Mca+BrCXMZbkPO8UY+w0YyzeeVztmkRb4Jia+xRj7DHG2I/1J95ijB1mjA3v5h8FIYR4FCUhpD+aAeA65zyFcz4MwJfO/VWc89EAXoVjCn8A+AnAWM75SDiWZ6+/+nAagDs45/cBWAHgZc75CDjWqilo9JpPAfiRcz6Cc/4igDcBLAIAxthgAPKeOLMrIYS4EyUhpD86A2Cqs/IxqXaxNDimf6/9Os65HQngK8bYGTiWWk+qd53/cM5Nzu2fAfyVMfYkgJh6+5vzIYCZjDEpgCUA3u5UiwghpBeiJIT0O5zzC3BUMc4AeJ4x9nTtU/UPc359BcCrnPNkAI/AsVBeLUO9a+4DMBuACY6kZUorMRgBpAO4A441UvZ1uEGEENJLURJC+h3GWDgAI+f8/wC8ACDV+dS8el9/dm77Arjm3H6ohWsOAJDLOd8Bx+qpjft3VANQN9r3JoAdAH7tqQvQEUKIOzXbq5+QPiwZwD8ZYwIAK4CVAD4CIGeMZcCRnC9wHrsBwIeMsWsAfgEQ18w15wF4gDFmBVAEYGOj508DsDlX/X2bc/4i5/w4Y6wKwFtd1zRCCOk9aBVdQgAwxvIAjOKcl3bja4YDOAQggXMudNfrEkJIT0G3YwjxAMbYQgAZANZTAkII6a+oEkIIIYQQj6BKCCGEEEI8gpIQQgghhHgEJSGEEEII8QhKQgghhBDiEZSEEEIIIcQj/n880m+VWh1egQAAAABJRU5ErkJggg==\n",
@@ -630,7 +649,7 @@
     "for n,(rank,df) in enumerate(data_by_rank):\n",
     "    df_sorted=df.sort_values(by=\"post_sparsity\",axis=0)\n",
     "    knots=list(zip(df_sorted[\"post_sparsity\"],df_sorted[\"pre_error\"]))\n",
-    "    temp=monotone_invert(knots)\n",
+    "    temp=fdict[rank]\n",
     "    plt.plot(temp.tvals,temp.yvals,label=\"N={:}\".format(rank),color=colorsequence[n])\n",
     "    try:\n",
     "        plt.plot(temp.tvals,temp.y_approx_vals,label=\"N={:} increasing\".format(rank),linewidth=5,linestyle=\"-.\",color=colorsequence[n])\n",
@@ -644,6 +663,55 @@
     "plt.close()"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": 132,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "[0.86   0.8775 0.895  0.9125 0.93  ]\n"
+     ]
+    }
+   ],
+   "source": [
+    "sparsityvals=numpy.linspace(start=0.86,stop=0.93,num=5)\n",
+    "print(sparsityvals)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 138,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEaCAYAAAAcz1CnAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd0FFUbwOHf3fROSYBAiAHpSO+9KogIKkVFQPBTFIIgKiIqioJUEQsCdlAUUQEpIqIISO+d0ElICIGQQBLSNtm93x+zYMjuJpuwqXufc/Ykmbkz8+7sZN+ZuWWElBJFURTFcemKOgBFURSlaKlEoCiK4uBUIlAURXFwKhEoiqI4OJUIFEVRHJxKBIqiKA5OJQKlRBFCPCiEOCiESBVCSCFESFHHZI0QIsQU47CijsUWJWnf3g3Te5tc1HEUJ85FHYCi2EoIUR74GTgHvAikAbFFGhQghHgZiJdSLirqWPKruO7bwiCEqA4MBX6TUh4q6niKgkoESknSAvAGpkopfy3qYLJ4GTgLLMo2PQLwADIKO6B8KK77tiB4AJlZ/q4OvAOEAw6ZCNStoWJACOGZn3n22kYJUsH080aRRmEjqUmTUhqKOhYblKh9m5UQwkkI4WpredNnkpl7SQcipVQvO7wAL2Aa2qV1OhANzAfKZiu3GYgCagPrgERgs2neZEACTYGFwFXtI7q9bBngE9Py6aZtvQe4ZdvGItN6qgI/AdeBC7nE/zTwpyluPXDRtC2fbOWcgInASSAFiAf2A6E27KM+wG9ApCn+GFOslWxYNtz0nrK+Nmd5v+EWlulsKtc5y7Rb+7gBMNe0j1OAP4B7LKzDH/gYuGCK+TKwEqhvmp89JnkrFiDE9PewbOusYor5immdJ4BxgLByrASb9luSaX8vzP6Z57Dfcj1mctq3Oay3H7ALLXHcBE4BC7OVkcAS4CHgIJAKnAfG3sXxd+szfQ541fR+Mm99xnmIa7Lp92FWPsPJwAum39vmcGw9X9TfPfZ4qVtDdiCEcAM2AvWBL9AOvlpAKNBaCNFGSpmeZRFP4G+0A/9VIPsZ43doXzjvAWWzbaMJ8DXaP1ZHYJJp2sMWQvsDCAPeANxyeRtj0G5vfIyWOJoCz6N9YXbJUu5t02sR8CHgDtQDOgCf5bKNZ9BuR95KcnXQ/qFbCSEaZ9tH2b0EPAiMAGYAx9G+SPPrG7T3+R4QiHZ7Zwna+wBACBEA7Eb7Il6ElvDKoO2PZqYYhqB9YV0B3jctetPaRk334ncAldD213mgN9q+vBcYnW0RD7TPfQswHmiD9rnEon32VuXhmMnTvhVCdAN+Af4F3kK79VUdy8dgY9P0hcC3wADgIyGEu5RyZpZyth5/t4xFO5a+QEvkl/MY1y3/mt7z66Z1bTVNP4J2a+9DtCS1I9tyQ9HqUZblsO6So6gzUWl4Aa+hncU0zTa9D9pZw7NZpm02TXvdwnomm+atw/zscJRp3ivZps81Te+dZdoi07SFeXgPXhamPW1aT+ss0w4Cv+dzP1naRifTNp6wYflhprLds01fRN6vCFZmK/uyaXq9LNO+Mk17yMK6RZbfo7BwBo2FKwJglmlav6zrAlaYpjewcKy8lG29q4GrNuyvvBwzFvetlfXOBRIA51zK3Tq7fjDLNBe0M/YUoFw+jr9bn+kVwO8u4pqc5e/u2T+nLPN+QLu6cM8yzQPtSn5pfv4PiuNL1RHYxxPAPuCiEML/1gvYiXY53M3CMgtyWN9CaTrisugDJGN+1j0ry/zs5ucauYmUMhlACKETQviZ4v/XNLtFlqIJQH0hRF1b121hG0II4WvaxnG0f7QWOS5sf9n3/ybTz3tB2w9Af2CblPL37Atb+Hxs1Qc4K6Vcnm1ds01/Zj97NQKfW4g1QAjhY8O28nrM2CIB7VZoTxvKnpJS/nHrDyllBtoVlAfaF/Ct6bYef7f8IKVMuIu4bPUN4Af0zTLtMcAH88YBJZZKBPZRB+2SPTbb6yraAV8hW/l4CwdxVhcsTAtBu8+flnWilPIy2hdpNRvXY5EQoqUQ4i+0L44baPGfN80um6XoJLR/ghNCiBNCiHlCCEuX7pa2UUsIsQLtbCqB//ZTmWzbKAwR2f6+bvpZzvQzAO0L4IidtxuCVr+S3QnTz+yf41UpZWq2adljzWlbeT1mbDEf7T2sEUJECyF+FEIMEkK4WCh72sK0U6aft7efh+PvFkvHdl7istU/aHUoT2eZNhS4BPx1F+stVlQdgX3o0M7+37Yy/3q2v7P/Y2dnbb61s1BhZV5u29EW1joObUI74F8z/UxFqxheT5YTBinlVlO764fQzugeBUKFEN9IKf+XwzZ80O5zZwDvAmfQbg9ItArtuzkpsbZfnHJYxlpLHmHjuu9GTuvMPi+nFkfZY83LtqwdM7mvUMqrQogmQFegB/AA8CTwmhCivZQyax1JrtvIy/GXhdmxnce4bCKllEKIxcBbQohAtP3WDZgtpTTmdX3FlUoE9nEWKCOl/LsAtxEOtDNVst0+wxNCVEI7cw2/i3X3RavAfkhKeXs9QojalgqbrmZ+BH4UQjgDi4FnhBAzpZSWzgBB++esBHSRUm7Osg0P7v5q4LqVdVS/i3XGol21NLShbF6+UMMBS7fV6maZby/hFNAxY7rF86fphRBiFNotqEFola63WDqGbk27dVafp+PPTnHdsVguq12EdpL3FFpScqIU3RYCdWvIXpYCdYUQg7LPMLVxzu0S3hZr0Dr8jMw2fXyW+fl166wz+/HwevaCplYvt0mtPfZR0585vU9r23jNwrS8OgP4CiGaZYnTFa2yNF9MZ3u/Ah2EEA9mny+EyHo2fhPbk9kaoIYQ4tFs63o1y3x7KZBjJvsxYHLA9DP7MVA76/4z3aYZg9bi5taJk83Hnx3jyurWlYLFz9CUnDah3R4aAuySUp6yVLakUlcE9jEH6AUsEUL0RrtNBFADrV3zW9z9GcRXaM0v5wgh6qD1gOyAdum71vTKr/Vo/5jrhBCfo/1jPgxY+scKE0JsA/aitdyojdbk8RRa80prtqOdZX8vhPgUrZ6gK1pFYNxdxA5aIp4O/CaE+Mg0bTB336P3DbTbX6uFEN+ifan4oMW9FK2ZL2j7YrAQ4h20e+I3pZTWvmRnAAOBpUKIW81HH0JrvvmZlPKoleXyo6COma+EEBXQmqZeROtr8QLa7ZoV2cqeAH4SQixAu68+EGgFvCmljDeVycvxZ6+4sjqOdptypBDiJlp/jWNSymNZynwLfG/6/YU8xlX8FXWzpdLyQmtP/yZwDO2gvgEcRmuhEZyl3GYgyso6JqNdptawMr8M8CnaP5Qe7UtkCtY7lOXYjC7bMt3Q2swno31hL0KrMM3e1G4iWpvqa6b3eQ6tFUhFG7bRFK3yLRHtds4KtArDcGCRDcsPw0oTR7T28fvQOk1FoV3Kd8N689Ea2ZYPwXLnr4polZCRpn0eDSznzmamVYDfTe9LYluHssWm/ZyO1tfjZax0KMthP4TYsM9sPWas7lsL6+yH1sT5Mv91nlwONMlWTnJnh7I0tNtBL9/F8deZbE2y8xnXZAvLHjPtI0vzPdD+p1PRbgMX+XeOPV/C9CYVRVHsSggh0Zp5Di7qWO6WqXPeZeBPKeWTRR2Pvak6AkVRlNwNQKtD+KaoAykIqo5AURTFCiFEV7R6sLfRbm8VZMvAIqMSgaIoinVvA+3QGkI8I0vpvXRVR6AoiuLgSsQVgb+/vwwJCSnqMBRFUUqU/fv3X5NSBuRWrkQkgpCQEPbt21fUYSiKopQoQojsY2pZpFoNKYqiODiVCBRFURycSgSKoigOrkTUESiKUnQyMjKIiooiLS0t98JKkXB3dycoKAgXl/w9eqHUJ4L4tHjKudtj8E9FcUxRUVH4+PgQEhLCnYOuKsWBlJK4uDiioqKoVi1/zxoqtbeGpJQsOLyA3it7c+7GuaIOR1FKrLS0NMqXL6+SQDElhKB8+fJ3dcVWKhNBhiGDSdsnMf/QfJL0SYz6exTXUq8VdViKUmKpJFC83e3nU+oSQZI+iZEbR7Lq3Krb06KToxm9cTQpGSlFGJmiKErxVOoSgcFoICY5xmz68bjjTPh3AgZjTo+AVRTFkbVt2xaA8PBwfvzxR7usc/369dSuXZsaNWowY8YMi2UuXrxIly5daNKkCQ0bNmTdunW35x05coQ2bdpQv359GjRoUCCV9qUuEZRxL8P8bvMp62b+1LnNUZuZuXcmanwlRXFcmZmZVuft2LEDsF8iMBgMhIaG8scff3DixAmWLl3KiRMnzMpNnTqVgQMHcvDgQX766SdGjRp1O9bBgwezcOFCjh8/zubNm/PdMignpbLVULBvMJ90/YRnNzxLuiH9jnlLTy4lyDuIofWHFlF0ilIyhbz+e4FvI3zGQ2bTkpOTGThwIFFRURgMBiZNmsSECRN4/PHH2bRpEwA//vgjNWrUYM2aNUydOhW9Xk/58uX54YcfqFixIpMnTyY6Oprw8HD8/f158803GT58OHq9HqPRyPLly6lZsybe3t7cvHmT119/nbCwMBo3bszTTz/NihUr+PTTT2ncuDEA7dq1Y8GCBTRs2DDH97Nnzx5q1KhB9erVAXjiiSdYtWoV9erVu6OcEILExEQAEhISqFy5MgAbNmygYcOGNGrUCIDy5fP69E7blLorglsaV2jM9A7TEZhXonyw7wP+jiiVw4orSqmzfv16KleuzOHDhzl27Bg9e/YEwNfXlz179jB69GheeuklANq3b8+uXbs4ePAgTzzxBLNmzbq9nv3797Nq1Sp+/PFHFi5cyNixYzl06BD79u0jKCjojm3OmDGDDh06cOjQIcaNG8ezzz7LokWLADh9+jTp6ek0bNiQTZs20bhxY7PXrVtMly5domrVqrfXGxQUxKVLl8ze4+TJk1myZAlBQUH06tWLTz/99Pa2hBD06NGDpk2b3vF+7KnUJgKA+++5n1eav2I2XSJ5fevrHI49XARRKYqSFw0aNODvv/9mwoQJbN26FT8/PwCefPLJ2z937twJaH0eevToQYMGDZg9ezbHjx+/vZ4+ffrg4eEBQJs2bZg2bRozZ84kIiLi9nRrBgwYwNq1a8nIyOCbb75h2LBhAHTp0oVDhw6ZvW7dYrJ0G9pSC5+lS5cybNgwoqKiWLduHUOGDMFoNJKZmcm2bdv44Ycf2LZtGytXrmTjxo153IO5K9WJAGBovaE8UfsJs+nphnTG/DOGyMTIIohKURRb1apVi/3799OgQQMmTpzIe++9B9z5hXrr9xdffJHRo0dz9OhRPv/88zsqVr28vG7/PmjQIFavXo2Hhwc9evTgn3/+yTEGT09P7r//flatWsXPP//MoEGDAHK9IggKCiIy8r/vmKioqNu3fbL6+uuvGThwIKAlqbS0NK5du0ZQUBCdOnXC398fT09PevXqxYEDB/K0/2xR6hOBEIIJLSfQKaiT2bz4tHhGbRxFQnpCEUSmKIotoqOj8fT0ZPDgwbz66qu3vwiXLVt2+2ebNm0A7f56lSpVAFi8eLHVdZ4/f57q1aszZswY+vTpw5EjR+6Y7+PjQ1JS0h3Tnn32WcaMGUOLFi0oV04brSC3K4IWLVpw5swZLly4gF6v56effqJPnz5m8QQHB98+0w8LCyMtLY2AgAB69OjBkSNHSElJITMzky1btpjVL9hDqawszs5Z58ysjrMY/udwTsTdWWMfnhjOmH/G8MUDX+Dm5FZEESpK8WepIrcwHD16lPHjx6PT6XBxcWHBggX079+f9PR0WrVqhdFoZOnSpYB2r33AgAFUqVKF1q1bc+HCBYvrXLZsGUuWLMHFxYVKlSrx9ttv3zG/YcOGODs706hRI4YNG8a4ceNo1qwZvr6+DB8+3ObYnZ2dmTdvHj169MBgMPDMM89Qv359AN5++22aN29Onz59mDNnDs899xxz585FCMGiRYsQQlC2bFlefvllWrRogRCCXr168dBD9v8cSsSjKps3by7t8WCa2JRYnlr3FJeTL5vNezDkQWZ0nIFOlPqLJEXJk7CwMOrWrVvUYdzh1sOq/P39C22b0dHRdO7cmZMnT6LTFb/vCUufkxBiv5SyeW7LFr93U4ACPAOY320+Pi4+ZvP+CP+DTw9+WgRRKYpS3H333Xe0atWK999/v1gmgbtV+t5RLmqUrcHcLnNxFuZ3xb46+hW/nv61CKJSFCUvbvUHKCxDhw4lMjKSAQMGFNo2C5PDJQKAVoGtmNx2ssV5U3dNZdulbYUbkKIoShFyyEQA0LdGX0Y1GmU23SANvLL5FU7FnyqCqBRFUQqfwyYCgBcavUCfe82bcqVkpjBq4yiLg9cpiqKUNg6dCIQQTG4zmVaVWpnNu5pyldCNodzU3yyCyBRFUQqPQycCABcnFz7s8iE1ytQwm3f6+mle2fIKGcaMIohMUZTCVtyGoc7IyODpp5+mQYMG1K1bl+nTp9slpuwcPhEA+Lr68lm3z/D3MG+FsCN6B+/vel8NXa0opURJGob6l19+IT09naNHj7J//34+//xzwsPD7zqu7ByiZ7EtKntXZl63eQxfP5zUzNQ75i0/s5wq3lV4ruFzRRSdohQDk/0KYRvmw7048jDUQgiSk5PJzMwkNTUVV1dXfH19724fW1C6rwj2fQPnch5MKqv65eszu+Nsi72LPzn4Cb+fL/jx2BVFuZMjD0Pdv39/vLy8CAwMJDg4mFdfffX2OEf2VHoTwdmN8PsrsKQ/7FoINt7a6VS1ExNbTrQ4b9L2SeyLufuhLhRFsZ0jD0O9Z88enJyciI6O5sKFC8yZM4fz58/ncQ/mrsASgRCiqhBikxAiTAhxXAgx1jR9mRDikOkVLoQ4ZPeNXzsLvw4HaQRpgPUTYPWLkJme+7LAE3WeYFj9YWbTM4wZjN00lvMJ9v8gFEWxzJGHof7xxx/p2bMnLi4uVKhQgXbt2mGPcdeyK8grgkzgFSllXaA1ECqEqCelfFxK2VhK2RhYDqyw61bTEmDpE9rPrA5+D4v7wM2rNq1mXLNx3H/P/WbTE/WJjPp7FHGpcfaIVlGUXDjyMNTBwcH8888/SClJTk5m165d1KlTJ8/7MDcFVlkspbwMXDb9niSECAOqACcAhJbCBwJd7bphowF8KkHcGfN5kbvgiy7w5FIIzLmSRyd0TGs/jaspV82eZHbp5iVe/OdFvu7xNR7OOV9SKkqpYaEitzA48jDUoaGhDB8+nPvuuw8pJcOHD8+1gjo/CmUYaiFECPAvcJ+UMtE0rSPwobUhUoUQI4ARAMHBwc0iIiJs3t7KfRdI+O01hjlvsFzAxRMeWQD1H8l1XfFp8QxeN5jIJPMnmXUL7sacTnNw0jnZHJuilDRqGGqNGob6LgghvNFuAb10KwmYPAkstbaclPILKWVzKWXzgIAAm7cXGZ/ChJUnmZw5jIkZ/yNDWviSzkiBX56GTdPAaMxxfeXcyzG/23z83Mybzm28uJEP9n1gc2yKopRMahjquyCEcEFLAj9IKVdkme4MPAYss/c2q5bzZNLD9XDSCZYauvGU/g3ipPnzBwDYMhN+GQrpOQ8jEeIXwqddP8VV52o2b0nYEn4I+8EeoSuKYiM1DLV9FWSrIQF8DYRJKT/MNrs7cFJKGVUQ2x7S+h6+f6Ylfh4u7JF16aufSpgx2HLhsDXwTQ+4nvOtpyYVmvB++/ctzpu5Zyb/XLS9v4KiKEpxUpBXBO2AIUDXLM1Fe5nmPUEOt4XsoW0Nf1aFtuPeAC+iZAD99JNZb2hhufCVY/BlF4jYkeM6e1brybhm48ymSyQT/p3AsWvH7BG6oihKoSqwRCCl3CalFFLKhreai0op15nmDZNSLiyobd8S4u/FytB2dKoVQArujMwYy8eZj1kunBIHix+G/YtyXOfw+sMZUMv88jDNkEboxlCikgrkIkdRFKXAlL5aj2x83V34+unm/K99NSQ65mb2Z5R+DKnS/H4/xkxYMxbWjQeD5RFHhRC80eoN2ldpbzYvPi2eURtHkZBeNM3sFEVR8qPUJwIAZycdk3rXY2a/Brg4CdYZW9NfP5lLsrzlBfZ8AUseg5R4y+vTOfNBpw+oW868Sd2FhAuM2zwOvUFvz7egKEohKKphqCMiIujWrRsNGzakc+fOREX9d2dh8eLF1KxZk5o1a+bYSe5uOEQiuOXxFsEs+V8rynm5clyG0Dd9KvuMtSwXvvAvfNkVrp60ONvLxYt53eZR0bOi2by9MXt5Z8c7auhqRSmGiuMw1K+++ipDhw7lyJEjvP3220ycqI13Fh8fz7vvvsvu3bvZs2cP7777LtevX7/ruLJzuGGoW1Uvz6rQdvxv8V5OX4FB+jeZ4vwtjztvNi98/QJ81R36fQW1e5rNruBZgfnd5zP0j6EkZyTfMW/t+bVU8a7C6CajC+idKErharC4QYFv4+jTR82mOcIw1CdOnGDu3LmANmzFI49onV3//PNP7r///ttDWtx///2sX7/+9oB79uJQVwS3VC3nyfKRbelWpwJ6XJiQ+RzvZgzBIM1HBUSfpI1dtG2uxRFMa5WtxYedP8RZmOfUz498zsozKwviLSiKw3CEYagbNWrE8uXLAVi5ciVJSUnExcXZvPzdcshEAODj7sIXQ5vzfKfqgOBbw4MMy5hAgvS0UFrC35NhxXOQkWo2t23ltrzd5m3zxYD3dr7Hjuicm6UqimKdIwxD/cEHH7BlyxaaNGnCli1bqFKlCs7OzjYvf7ccNhEAOOkEEx+sy5wBjXB10rHV2JBH9FM4Zwy0vMDRX+DbByEx2mzWozUfZUTDEWbTM2Umr2x+hdPXT9s7fEVxCI4wDHXlypVZsWIFBw8e5P33tY6rfn5+Ni9/txw6EdzSr1kQS0e0wt/blQsykEf0U9hkaGS5cPRBbQTTqP1ms0Y3Hk3v6r3Npt/MuMmov0dxNcW2IbAVRfmPIwxDfe3aNYymcc+mT5/OM888A0CPHj3YsGED169f5/r162zYsIEePXrkbQfawOEqi61pdk85Vo1uz7OL9xF2Gf6XMZ4JcinPO1t4POXNGO3KoM+n0Ojx25OFELzb9l2upFxhb8zeOxa5knKF0I2hLOq5CC8Xr+xrVJRiz1JFbqFs1wGGod68eTMTJ05ECEHHjh357LPPAChXrhyTJk2iRYsWt5cpiEdVFsow1HerefPmsiCeymNJcnomL/98iD+PXwHgMd2/THf5CjdhpclZ2zHQfTJkGYo6IT2BIX8M4UKC+UHYvkp7Pu36Kc46lYOVkkENQ61Rw1A7EC83ZxY81YwXu9YAYIWxI0/oJ3FVlrG8wI5PzJ6I5ufmx/xu8ynnbp65t13axrTd01QfA0UpQdQw1A5IpxO88kBtPn6iMa7OOg7KmvRJn8IRYzXLC5zZoPU3iDt3e1KQTxCfdfsMdyd3s+K/nP6Fb49/W1DhK0qpp4ahti+VCHLQt3EVfn6+DQE+bsRQngH6d1hlaGu58LXTWk/kc5tuT7rP/z5mdpyJwLy519z9c1kfvr6gQlcURbGZSgS5aFy1DKtHt6NBFT/ScWVsRiizMh7HaKnzWdoNWNIPdi283fmsa3BXJrScYHHdb259k4NXDxZk+IqiKLlSicAGgX4e/Px8Gx5qEAgI5hv6MiLjZW5K89s+SAOsnwCrX4RMbeC5p+o+xeC6g82K6o16XvznRSISbX8es6Ioir2pRGAjD1cn5g1qwrju2iB1fxub8aj+PSKMFSwvcPB7+K4P3IwF4NXmr9ItuJtZsYT0BEb+PZL4NMsjnSqKohQ0lQjyQAjB2O41+WxQU9xddJyRQfTVT2GHoZ7lBS7u1J58dvkITjonpneYTgN/84G7IpMiGfPPGNIy0yysRFGUwlIch6Hu2bMnZcqUoXdv886q9qISQT481DCQX19oSyVfd27gw9CM11mceb/lwgmR2jORT6zCw9mDT7t+ShXvKmbFDsce5o1tb2CUxgKOXlEcW0kahhpg/PjxfP/993cdS05Ur6Z8uq+KH6tHt+O57/dzOPIG72QO55QM5l3nRbgIw52FM1Lg56HQ6XXKd5rA/O7zGbJuCIn6xDuK/RXxF3P3z+WV5q8U4jtRFNuE1Sn4TmV1T4aZTXPkYagBunXrxubNm/O9T22hrgjuQgVfd5aNaM0jjbVBoH40dGOw/g3ipbflBbbMgF+eprpHRT7u8jEuOhezIouOL+Knkz8VZNiKUqI48jDUhUUlgrvk7uLE3McbM75HbQB2y7r00U8lzFjV8gJhq+HrHjR3r8CUdlMsFpm+ZzpbIrcUVMiKUqI48jDUhUUlAjsQQhDapQafD2mGp6sTUbIC/fTv8qfByhAfV47CF114yKksY5qMMZttlEbG/zue43HHLSysKI7FkYehLiwqEdhRj/qV+PWFtlQp40EK7ryQ8RIfZz5quXDKNVjch2f1LvSr2c9sdmpmKqM3jib6pvmzDxTFkTjyMNSFRVUW21m9yr6sGt2OF77fz76I68zNHMBpY1U+cFmIh9DfWdiYgVg7ljdbPEdMYFu2X77zSWbXUq8RujGUxQ8uxtfVtxDfhaKYs1SRWxgceRhqgA4dOnDy5Elu3rxJUFAQX3/9td2fSaCGoS4g6ZkG3lhxjOUHtPbA9cUFvnSdQ2VhuePYzWodeNrPmdMJ58zmtarUigXdF+DiZF65rCgFTQ1DrVHDUCt55ubsxAcDGvJGrzoIAcdlNfqkv88+Yy2L5b0vbOWzixeo4GY+dPXumN1M3jlZDV2tKEVEDUOdT0KIqkKITUKIMCHEcSHE2CzzXhRCnDJNn5XTekoyIQQjOt7L1083x9vNmWv4MUj/Jj9ndrJYvlJ8OPMvhuOpczObt/rcahYeXljQIStKiaCGobavgkxtmcArUsq6QGsgVAhRTwjRBegLNJRS1gc+KMAYioWudSqyYlRbqpbzQI8Lr2WO4L2MIRgsjGBaO+UGc6KjcLIwdPX8w/NZdXZVYYSsKIoDKbBEIKW8LKU8YPo9CQgDqgAjgRlSynTTPId4onutij6sCm1Py2rlAME3hgcZnvEaidLTrGz71FTeunbN4nom75jM7su7CzhaRVGJ/xqGAAAgAElEQVQcSaHc7BJChABNgN1ALaCDEGK3EGKLEKKFlWVGCCH2CSH2xcbGFkaYBa6clytL/teKJ1ponc3+NTair34K54yBZmX7JyXz7I0Es+mZMpNxm8Zx9vrZAo9XURTHUOCJQAjhDSwHXpJSJqI1WS2LdrtoPPCzsNDVTkr5hZSyuZSyeUBAQEGHWWhcnXVMf6wB7zxcD52ACzKQR/XvsdnQyKzsi9cTePBmstn0pIwkRm0cRWxK6UiQiqIUrQJNBEIIF7Qk8IOUcoVpchSwQmr2AEag8Gp9igEhBMPbVePb4S3xcXcmES+eyRjPF5kP3VFOB0yNjaNpmvnw1JeTLxO6MZSUjJRCilpRSr/iNgx1REQEzZo1o3HjxtSvX5+FCwumwUhBthoSwNdAmJTywyyzfgO6msrUAlwByzfES7lOtQJYOaodIeU9MaJjWuZTvKJ/gXT5Xz8/V+DjK9cI0WeYLR8WH8Zr/75GptH6sLqKotypJA1DHRgYyI4dOzh06BC7d+9mxowZREfbf7SBguxZ3A4YAhwVQhwyTXsD+Ab4RghxDNADT0sHbiBfo4I3v4W2I/THA2w/G8dyY0fO6wP53HUuFcQNAMoYjcy/cpXBlSsR7+R0x/JborYwY88M3mz1psXBrBTFXj57IefxeOwhdGFXs2mOPAy1q6vr7TLp6em3h6Gwt4JsNbRNSimklA2llI1Nr3VSSr2UcrCU8j4pZVMpZcEfXcVcGU9XFg1vyZDW9wBwUNakT/oUjhir3S5TNdPAJ1dicbNwICw7tYzvTnxXaPEqSmFy9GGoIyMjadiwIVWrVmXChAkWB627W6Wvi1wJ5eKkY8oj9zGlb32cdIIYyjNQ/zarDW1ul2mUrmdGbBzCwgXUnH1z+Cvir8IMWVEKhaMPQ121alWOHDnC2bNnWbx4MVeuXMnL7rOJSgTFzJA2IXz/TEv8PFxIw40xGaOZlTHw9vzuKam8Gn/DbDmJZOLWiRy6eshsnqKUZGoY6v/K1K9fn61bt+a+0/JIJYJiqG0Nf34LbUf1AC9AMN/wCM/pX+amdAdgSGISgxKSzJZLN6Qz5p8xXEy8WMgRK0rBceRhqKOiokhNTQXg+vXrbN++ndq1a+dh79lGDUNdTFXz92LlqHa8uPQg/56O5S9jcx7Tv8tXLh8QrIvltfjrRDs7sdnrzp7J19OvM2rjKJY8uIQy7mWKKHqlNLJUkVsYHHkY6rCwMF555RWEEEgpefXVV2nQoEF+dmOO1DDUxVymwcj768L4dns4AGVIYoHLx7RxOkGKEDwTWIHjbuaD1DUMaMhHnT8iwLP0dMZTioYahlqjhqFWioyzk453Hq7P9Mca4KwT3MCHIRmv813m/XhKybwrsVTOMG8XfST2CP3X9GdrlP3vJyqKoyntw1CrK4ISZNf5OEYu2c/1FK1z2VNOfzPZeTEXXQVDAiuR5GT5AB1abyhjm47F1cnV4nxFyUlxvCJQzKkrAgfRunp5VoW2p1ZFbwB+MHRnSMZEyurd+ehqLM5Wkvp3J75j8LrBRCRGFGa4SilSEk4YHdndfj4qEZQwweU9WT6yLV3rVABgl7EeffRT8U2pwIKYq5TPNFhcLiw+jAFrBrD63OrCDFcpBdzd3YmLi1PJoJiSUhIXF4e7u3u+16FuDZVQBqNk1vqTfP7veQC8SGWOy0Kau+znrYDybPe03kGmd/XevNX6LbxcvKyWUZRbMjIyiIqKuqNNvlK8uLu7ExQUhIvLnc81t/XWkEoEJdwv+yJ5c+Ux9AYjIBns9DdvOi/hZz93PipXhkwr4w9V9anK7I6zqe9fv3ADVhSl0Kg6AgcxoHlVfnyuFf7eroBgieF+HtFPofUNP5ZEX6FqhvmopQCRSZEMXjeYRccWYZQFM5CVoiglg0oEpUDzkHL8FtqOuoG+AJySwfTRT+VASkd+uRRDbwsPtwHtaWdz9s9h1N+juJbqkCOBK4qCSgSlRlBZT1aOasvznaqjE5COK5Myn2Fc+ktMuJrGtNhreFgZwnZ79Hb6re7Hjks7CjlqRVGKA5UIShF3FycmPliXX15oS3V/rSJ4g7EFPdNnEJB4D79ciqFuut7isvFp8Tz/9/N8uO9DMgyWbycpilI6qURQCjW7pyzrxnbg2fbVEAJiKM8g/Zv8kvooiy9dZWhCotVlvz3+LUP/GEpkYqTVMoqilC4qEZRS7i5OvNW7Hj8/34Z7TI/CnGd4lKf0b/NknDPzY65SzmC5z8GxuGMMWDuAtefXFnLUiqIUBZUISrkWIeX4Y2wHhrUNAeCArEWv9OncuNmY5Zcu09o0xG12yRnJTNw6kTe3vUlKRkohRqwoSmFTicABeLo6M7lPfZY+15qq5TxIxIsXM15kZtqzzL2cyLj461aHp1h9bjUD1w7kRJz5A7cVRSkdVCJwIG3uLc/6sR1Nz0YW/GLoTB/9+7S8Xp7Fl69QxcIopgARiRE8te4pvjv+nRpmQFFKIZUIHIyXmzNTHrmPH55tRZUyHpyXlXlM/y67k7vzy6XLPGitz4Exk9n7ZhO6MZS41LhCjlpRlIKkEoGDalfDn/UvdeDJlsHoceH9zMGMTn+NV69mMCU2zmqfg62XttJ/TX92Ru8s5IgVRSkoKhE4MB93F6Y/1oDFz7Qk0M+dLcZG9EqfSdmEe1kWHUMdK30OrqVe4/m/nuej/R+RYVR9DhSlpFOJQKFTrQD+HNeRgc2DuIYfwzJeY0nK43x7KZbBVvocSCRfH/uaYX8MIyopqpAjVhTFnlQiUADwdXdhVv9GfDusBRV8Pfja8BBP6N9lwDV35sVcpayVPgdHrh1hwJoB/HHhj0KOWFEUe8k1EQghnIQQ4wojGKXodalTgQ0vdeKxplU4JqvTWz+NK0mt+PVSDK1SLY9HfzPjJq/9+xqTtk9SfQ4UpQTKNRFIKQ1A37yuWAhRVQixSQgRJoQ4LoQYa5o+WQhxSQhxyPTqlY+4lQLk5+nChwMb8+XQ5nh6+/Fa5vNMSR3J7MtJjI2/gZOVJqS/nf2Nx9c+zsn4k4UcsaIod8OmB9MIId4H/IBlwO32hVLKAzksEwgESikPCCF8gP3AI8BA4KaU8gNbg1QPpik615P1TF5znFWHogkSV/nEZR46j4tMCPAn2sXZ4jIuOhdeaf4Kg+oMQlh5MI6iKAXPrk8oE0JssjBZSim75iGgVcA8oB0qEZQ4fxy9zFu/HSMhOYWxzisY4rKaqf5l+dPb+uMuOwV1Ykq7KZR1L1ugscmMDES2R/QpimLnJ5RJKbtYeOUlCYQATYDdpkmjhRBHhBDfCCEK9ltCsYsHGwSyYVxHHmgQxJzMgYxMn8jLVyXvxsbhbqXPwZaoLfRb3Y89l/cUWFxpp05xtkcPkjZvLrBtKEppZ1MiEEL4CSE+FELsM73mCCH8bFzWG1gOvCSlTAQWAPcCjYHLwBwry424tb3Y2Fib3oxSsMp7uzH/qWbMG9SEk+6N6ZU+He+EOiyLjqGWlT4HsamxPLvhWT458Ind+xwY09OJHv8amdGXiXphJJcnT8aYoiqrFSWvbG0++g2QhHZ/fyCQCHyb20JCCBe0JPCDlHIFgJTyipTSIKU0Al8CLS0tK6X8QkrZXErZPCAgwMYwlcLQu2FlNozrRMt6NRiR8TLfpAzhm+g4BiUkWSwvkXx59EuGrx/OpZuX7BZH7MefkH769O2/b/y0jAuPPkbq4cN224aiOAJbE8G9Usp3pJTnTa93geo5LSC0WsKvgTAp5YdZpgdmKfYocCyvQStFL8DHjc+HNOPjJ5qwxrUXA9Kn8tg1Lz65EouflT4Hh2MPM2B1f/4M//Out5+8ew/x35qfi+gjIggf9BSxn85DZqhez4piC1sTQaoQov2tP4QQ7QDLA9n/px0wBOiaranoLCHEUSHEEaALoPoolFBCCPo2rsJf4zpStXYz+uinEpnYnuWXYmhupc9BUsZNXt3yKpN3TCY1M7dDyDrXaiF4dWhveabBwLXPPiP8qcGkX7iQ720oiqOwtdVQI+A7tCakANeBp6WURwowtttUq6HiT0rJigOXmLzmOG30O5nm8iXLy+pYUMYPg5UmpNX9qjOr4yxql6ud723e+OknrsychUyznHiEhwcVJ7xGmccfV01ZFYdjt1ZDQggdUFtK2QhoCDSUUjYprCSglAxCCPo1C+KvcZ3Q1+xF7/TpNIwL4tvLVwjMtPycg/MJ5xn0+yCWnlyar+ccCCEo++STVFuxAvcGDSyWkampxEx+l8jnnydTNTpQFIts6VlsBEabfk80tfxRFIsq+bnz7bAWvNyvC8/r3uavm4/wU9QV7k+23JpHb9Qzbfc0xm4ay420G/naplv1aoT8+AP+o0aBk5PFMsn/buX8w31I/OuvfG1DUUozW28NTUKrE8jeszi+4EL7j7o1VDJdupHK68uPkHx2Bx+5fMZuvxRmlitLus7y+UcFzwrM6DCDFpVa5HubqYcOcWnCBDIiLlot4/fYY1R8YyJO3t753o6ilAT27llsqcZNSilzbDlkLyoRlFxSSpbuieTj3/fxpvyCeu77GF+hPGddXS2W16FjRKMRPN/weZx1loewyI0xOZkrs2ZzY9kyq2VcqlSh8qyZeDZrlq9tKEpJYLdEYKojaCOl3G6v4PJKJYKSLzI+hdd+OUyViyuZ6LKIBeU9WebrY7V8kwpNmNlhJoHegVbL5CZp82Yuv/kWhjgrj9YUgvLPPkvAi6MRVhKTopRkdqssNtUR2DwukKJYUrWcJz8815pGvUcx2DCDvrF+fHQlFl8rfQ4OXj1IvzX9+Csi//f0fTp3pvqa1Xh362a5gJTEffklFx5/gvQzZ/K9HUUp6WztR7BBCNFPqPZ3yl3Q6QRD2oSw8KXHmVb5U84ldmX5pRiaWmn6maRP4uXNL/PezvdIy7RcJjfO5coRNO9TAt+fis7T02KZ9LAwLvTrT/x33yGtjJukKKWZrXUESYAnYADSAIFWR+BbsOFp1K2h0sdolCzeGc729cuYqpvPyrKShWX8MFo516hR5l5mdZxNzbI1871NfWQk0a9NIPXgQatlvNq2IXDaNFwqVcr3dhSluLDr6KNoHcmGAVNNX/71gfvzH57i6HQ6wfB21Xhz7Iu8EbCAunHV+ObyVSpZ6XNw9sY5nvz9CX4+9XO++hwAuFatyj1LvifgpZfA2XJFdPKOnZzv05eE33/P1zYUpSSy9YpgAWAEukop65qGjt4gpcx/O788UFcEpZvBKPl22zmu/vURzzkv5f2AMmz0snwbB6B7cHcmt52Mn5tNA+BalHr8ONGvTUB/7pzVMr69e1Np0ls4+eV/O4pSlOx9RdBKShmKdlsIKeV1QDWzUOzCSSd4tmMNBr44g+l+HzI6xoVJ1+Jxs3K//u+Lf9N/dX/2X9mf72161K9PteW/UnbwYKtlEteu5XzfR0jetSvf21GUksDWRJAhhHACJIAQIgDtCkFR7KZGBW9mvTiEzZ1+QSQ0Z2n0Fe7VW37OQUxKDM+sf4YFhxdgMFpueZQbnbs7ld56k6pffYVzhQoWy2TGxHBx2HCuTJ+BMT09X9tRlOLO1kTwCbASqGB6fvE2YFqBRaU4LGcnHc92a0Cj0CWsdH+ZLy4lMCDR8nMOjBiZf2g+/9vwP2KSY/K9Te/27ai+ehU+PXtaLRO/eDHh/fuTdvJkvrejKMWVTXUEAEKIOkA3tBZDG6WUYQUZWFaqjsAxZRqM/LhhG/ftfJk47yje8S9PkpPlcxc/V1+eb/QCfe7tk++6AykliWvXEvPeFIxJlpMPLi5UGDuGcsOHI6yMa6QoxYVdh5goaioROLawS/Ec+v51OmQs540K5Tjo7m61rKvOlQdCHqBfzX40q9gsX0NPZ0RHE/36RFL2WH/Wsmfz5gTOmIFrUJU8r19RCotKBEqpkmEwsuq3n2l9ZCK/lTXwRRlfZC5f8iG+IfSv1Z8+9/ahrHvZPG1PGo3EL1pM7Ny5Vp90pvPyouJbb+H3SF/1rAOlWFKJQCmVws5FEP/T87g6HeD1gPJctdIfICsXnQvdg7vTr1Y/WlRqgU7YWjUGaadOEz1+/B3PRs7O54EHqPTuZJzL5i3ZKEpBU4lAKbX0GQa2LJ1Jgwtzme3vwwZvL5uXDfYJpl+tfvS9ty/lPcrbtIxRryf2o4+1ZyRb+X9xCvCn8rRpeHfoYHMsilLQVCJQSr1TR/fgvPI5hFMUy328We3txQ0bK3CdhTNdgrvQv1Z/Wge2tukqIXn3HqInvk5m9GWrZcoOepIK48ej8/Cw+X0oSkFRiUBxCOn6dP78+XNCziymjjjLRi9PfvXxZo+H9Qrl7Kp4V6FfzX48UuMRAjwDcixrSEriytSpJKxabbWMa0gIlWfPwsPK4zMVpbCoRKA4lOs30ykbfwh2fQZha4hw0rHcx5tVPl7E23iV4CSc6BTUif61+tO2clucdNaXS1y/nsvvTMaYkGC5gLMz/qNG4j9iBMKGegxFKQgqESiO63oE7PkCDnxHRnoi/3h6sNzXm515uF0T6BXIozUf5dEaj1LJy/JIpBlXrnD5jTdJ3m79mU0ejRpRedZMXO+5J89vQ1HulkoEipKeBAd/gN0L4Ho4kc5OrPTxZqW3N9ecbbtK0AkdHap0oH+t/rSv0t7s8ZnSaOT6j0u5Ons20soQFMLTk4qvT6DMgAGqmalSqFQiUJRbjAY49Qfsmg8R28kA/vX04Fcfb7Z7uOfaH+GWCp4VeLTGozxW8zEqe1e+Y176+fNEj3+NtOPHrS7v3bkzgVOn4OzvfzfvRlFsphKBolgSfQh2LYBjy8GYQbSzEyu8vVnp42VTnwQAgaBdlXb0r9mfjlU74qJzAUDq9cTOn0/cF1+ClZFTncqVI3DqFHy6drXbW1IUa1QiUJScJF6GvV/Bvm8gNZ5MYJuHB7/6erPVw93qk9Ky8/fw55Eaj/BYzceo6lMVgJQDB4meMIGMyEiry/n170fF1yfilIc+EIqSVyoRKIot9ClwZJl2lXDtFAAxTk6s9PFihY83MXlo8dMmsA39avWja9Wu6FL1XJ05gxu//Gq1vEvVqlSeORPPpk3u+m0oiiVFngiEEFWB74BKaM8u+EJK+XGW+a8Cs4EAKeW1nNalEoFS4KSEsxu1eoRzGwHtAd07PNz51cebLZ4eGGy8SijnXo6+9/alX61+lNt7lsuT3sYQH2+5sE5H+RHPERAainBxsdObURRNcUgEgUCglPKAEMIH2A88IqU8YUoSXwF1gGYqESjFytUw7QrhyDLITAMg1knHb97eLPfx5pKL7VcJLSu1ZGBAD+p8vpGULf9aLedevz6VZ83E7d577zp8RbmlyBOB2YaEWAXMk1L+JYT4FZgCrAKaq0SgFEvJ12Dft7D3S7h5BdAubXe5u/OrrzebPD3ItPEqoYyrH6Mv1qbh0v2QZqWZqZsbFcaPp+xTg1QzU8UuilUiEEKEAP8C9wGdgW5SyrFCiHCsJAIhxAhgBEBwcHCziIiIAo9TUSzKTIdjK7ReyzFHb0++ptOx2seL5T7eXLTxtk7FeMnr692pEpFstYxXu3YETpuGS0XLj89UFFsVm0QghPAGtgDvA+uBTcADUsqEnBJBVuqKQCkWpISI7bBzPpxah+kR3hiBve5uLPfx5m8vTzJyOZvXGSWP7pD032bEycq/n87Pj8B3J+Obw+MzFSU3xSIRCCFcgLXAn1LKD4UQDYCNQIqpSBAQDbSUUlp96KxKBEqxE3cOdn8OB5dAxn9n99d1OlZ7e/GrjzfhrjlfJdwbLXlxjYHKVuqRAXz7PEylSZNw8vGxV+SKAynyRCC0m5yLgXgp5UtWyoSjrgiUkiz1Bhz4ThvbKOG/fgMS2G+6Stjg6YleZ/kqwTVDMvgfIz0PWP8/dK4cSOUZM/Bq2dLe0SulXHFIBO2BrcBRtKtngDeklOuylAlHJQKlNDBkwsk12m2jqDufdZyg07HWWxse+6yrq8XFG50zMup3I2WtVR0IQbnhwwl4aSw6K+tQlOyKPBHYk0oESokStQ92fgYnVoE03J4sgcNurvzq482fXp6k6e58GI5PimTEeiOtTln/n3SrVYvKs2fjXrtWQUWvlCIqEShKUbsRqTU93b8I0u58bkGiTvC7l1aXcNotyxm+lHQ6Khn+lxFPveXVGp2dKDt2NIH/G4HQ2f78ZcXxqESgKMVF+k04vFTrtRx//o5ZEjjm6sqvvt784eVJqumLPeCGJHStgXrWhyviSp0Aqk2Zzj0N2hVg8EpJphKBohQ3RiOc+VO7bRS+1Wz2TSFY5+3Fch8vTri5IYySh/dInthixNnyYKYYgatta9J0wnT8atcv2PiVEkclAkUpzi4fgd0L4egvYDC/B3Tc1YXlPt6s8/bCP1bw4moDwTk0qTAKSOvYlPqvvIN7LVV/oGhUIlCUkiDpCuz7GvZ+DSnm3/QpQrDey5OVHt403OVM7705/79KAU5d23PPmPGqQllRiUBRSpSMNDj6s9b8NDbMYpFTri5sueFLvU1ulE/KfZUe93el0ugxuNeubedgC1/sxST8Knjg6m77gH+KSgSKUjJJCec3aaOfntlgsUiSQfD3pXJUOuxGGetDFt3m3b07AaGjcK9b187BFo7MDANLJu3CkGmkyQPBNOgchIurbc+cdnQqEShKSRd7GnYvgENLITPVbPZ54cyGS+VpuN/Jeke0LLy7dcN/1Eg86pesSuUjmyLZuuzM7b89fF1p1uMe6neojLNKCDlSiUBRSouUeNj/Lez5EpIu3zFLIvinz0fs+uELOm+5Trmbua/Ou0sX/END8biv+CeEDL2BJW/tJCXRvELdy8+VZg+GUK9dZZxcVH8KS1QiUJTSJlMPJ37Tmp9ePqRNq/UgDPqJtMw0Fh/4kgs/fEXvHXqb6hC8O3XCf3QoHg0aFGzcd+HQ3xfZ/uvZHMt4l3Oj+YMh1GkbiJOTSghZqUSgKKWVlHBxl/Z8hJbPQ7UOt2ddunmJD3fMgLX/8MhOI/42JASvjh0ICA3Fo1GjAgw6fxJiU9j3ezindseQ21eVr787zXtVo3ariuhUQgBUIlAUh7b90nZm75hO9W0XeHSnEf/E3Jfxat8e/9BReDZpUvAB5tH1mGT2/h7OmX1Xbj0Gwiq/Ch60eKgaNVtURGdl1FdHoRKBoji4DEMG34d9z5cHFtDqQAqP7jRSISH35bzatsV/dCieTZsWfJB5FBd9k71rL3DuQGyuZctW8qRF72rUaFoB4aAJQSUCRVEAuJJ8hTn757Dh7Do6HZU8tsO2hODZpjUBoaF4Ns/1e6TQXYtKYs+aC1w4nOMI9gCUr+JFy97VqdbY3+GeBa0SgaIod9gbs5fpe6Zz/tppOh6TPLrDSKUbuS/n2aoV/qGjiuWDca5GJLJnzQUijsXlWjYg2IeWD1fjnvvKO0xCUIlAURQzmcZMlp1axmcHPyMlLZH2xyX9ttuYEFq0wD80FM9WLYvdF2nM+QT2rDlPZNj1XMtWrOZLy4erUbVuuWL3PuxNJQJFUayKS43jk4OfsOLMCnRGSfvjkse2G6mc+/coHs2babeMWrcudl+k0Weus3v1BaLP5J7ZAmv40fLh6gTVLlsIkRUNlQgURcnV0dijTNs9jWNxx9AZJe1OaFcIleNzX9ajaVPtllHbtsUqIUgpuXRKSwgx53OvDKlSuwwtH65O5RplCiG6wqUSgaIoNjFKI7+d/Y2P9n/E9fTrCKOkbZiWEIJyv/WOR+PG+IeG4tW+XbFLCJEn4tm9+jxXI3LvUFG1XjlaPlyNStX8CiG6wqESgaIoeZKQnsBnhz5j2allGKURYZS0OaklhKq5N87BvVFDAkJD8erQodglhPCjcexZc55rkbmPwXFPg/K0erg6AcE+hRBdwVKJQFGUfDkVf4ppu6dx4OoBAISUtD4p6bfNmOPDcW5xb9gQ/1Ej8e7UqXglBKPk/OFY9qy5QHx07qP0VWvkT8uHq+Mf5F0I0RUMlQgURck3KSXrLqxjzr45xKZqnbeElLQ6Kem/3Uhw7v25cL/vPvxHjcK7S+dilxDOHrjK3rUXuB6Tkmv5e5tWoGXvapSr7FUI0dmXSgSKoty15IxkPj/8Od+f+J5MmQloCaHlKe2WUcjV3NfhXq8e/qGj8O7atVglBKNRcmbvFfauvUBCrPkw33cQUKtFRVo8VI0yFT0LJ0A7UIlAURS7OZ9wnhm7Z7Dz8s7b04SUND+tXSFUu5L7Otzq1sV/1Eh8unVD6IrPoHBGg5GTu2LYty6cpLi0HMsKAbVbV6J5r2r4BXgUUoT5pxKBoih2JaXkn4v/MGvvLKKTo7POoPkZSf9tRqrbkhBq18Z/1Ch87u9erBKCIdNI2I7L7P8jnJvX03Msq9MJ6rQNpHmvEHzKuRdShHmnEoGiKAUiNTOVb459wzdHv0FvzPLAGClpdlZLCPfG5L4et1q1tCuEBx4oXgkhw8jxbdHsXx9OSoL5A3Gy0jkJ6rWvTLOeIXiXdSukCG1X5IlACFEV+A6oBBiBL6SUHwshpgB9TdOuAsOklNHW16QSgaIUR5FJkczeO5tNkZvunCElTc5JBmwzUuOy5WWzcqtZA/+RI/Hp0QPhVHwePZmpN3Ds30sc+DOC1KSMHMs6Oeuo37EyTXvcg5df8UkIxSERBAKBUsoDQggfYD/wCBAlpUw0lRkD1JNSvpDTulQiUJTia2vUVmbunUlEYsSdM6Sk8XntCqFWjqd6Gtd778V/5Eh8H+xZrBJCRrqBo5ujOLjhImnJOScEZxcdDToH0eSBYDx8XAspQuuKPBGYbUiIVcA8KeVfWaZNBIKllCNzWlYlAkUp3vQGPd+f+J7Pj3xOama2FjhS0uiClhBqX8p9Xa7Vq+M/8gV8e/UqVglBn5rJkU2RHPo7kvSUzI+cyNUAABDsSURBVBzLOrs50ahLEI3vD8bdy6WQIjRXrBKBECIE+Be4T0qZKIR4HxgKJABdpJRmrZKFECOAEQDBwcHNIiIishdRFKWYiUmO4cN9H/JH+B/mM6WkQbh2y6hOVO7rcg0J0RLCQw8hnJ3tH2w+padkcGhjJIc3RpKRZsixrKu7E426VaVR92DcPAr/PRSbRCCE8Aa2AO9LKVdkmzcRcJdSvpPTOtQVgaKULHtj9jJt9zTO3rDw4HkpuS9Cu0KoF5n7ulzvuYfyL7yA38O9i1VCSEvO4OBfFzmyKYrM9JwTgpunM427B9OwaxCu7oX3HopFIhBCuABrgT+llB9amH8P8LuU8r6c1qMSgaKUPLeefTDv4DxuZlgY40dK6l+UPL5dUCci5y9SAJfgYPxHPIdPzwdx8i4+vXxTk/Qc2HCRY5ujyMww5ljW3cuFJg8E06BzEC5uBX/bq8gTgdC6EC4G4qWUL2WZXlNKecb0+4tAJyll/5zWpRKBopRc11Kv8fGBj/nt7G9Wy9S9KBmy05Ua53Pp4QsINze8O3fGt1cvvDt1ROdePNrxJyekc2B9BMe3RmPIzDkhePi40KxnCPU7VMbZteASQnFIBO2BrcBRtKaiAG8A/wNqm6ZFAC9IKXOsQlKJQFFKvsOxh5m2exon4k5YLVMnUjJirx9Bp2x4IAKg8/TEu3s3LSm0bYtwLfqWOjevp7H/jwhObI/GaMj5+9XLz5VmD4ZQr11lnFzs35eiyBOBPalEoCilg8FoYOXZlXx84GNupFt/ith90c6MPhhAuSM2VCKY6Pz88H3gfnx79cKzZcsib3GUGJfK/nXhhO2MQRpz/p71LutG814h1GkbiJOT/RKCSgSKohRbCekJzDs4j59P/4xRWr+N0jbenxf2l8V9X1ie1u/k749vz5749uqFR+NGRdpzOSE2hX2/h3Nqdwy5fd36+rvTvFcItVtVQmeHhKASgaIoxd7/27v74Cju84Dj30evSELHm3kVBiEOMMQ1xDYUsMG8JaATL2fXde2xY4/rJu2knThtM01f/qmncTKJO67bP5pJJ9TTTtrQ1oYToOMdJUAGsMWbzVsiISSMeX+xThJCCOnpH7uAIEK6lTgd7D2fGeZul92952G1eti93ed39NJRvr/7++w7t6/T5Z5r+R1erh4O5TtpvRzHwMrtZIwYTqC4mEAoRJ9Jk5LWAfXymUY+LquhsuIsdPFrt9/gHKYuHsO4qUNJS+t+vFYIjDEPBFVlbfVa3t3zLhea7j7yTUZaBq+Oe5nn68eTvnUX9Zs20dbQ9Yhj7WWNHk2gJEQgFCI7GOxp6N1y8VQDH689zrG9XQ/qMGBYLlMXjyH4+BCkGwXBCoEx5oHScK2Bn3zyE352+Gc3xz7oiCDMGDGDZ0eVMO1EFk3rNlFfXo42dX3HUXvZ48cTCIUIhIrJGjWqp+F7duFkPR+tOc7xA10P+zaoII/fXVrEmMmDPX2GFQJjzAOp+otqfvDRD9h1eleXy+Zn5RMaEyJcsIiHPzlHbF2Uxl9uQ1s67wl0pz6PPUYgVEyguJjMoUO7G3q3nKuN8dGa49QevNjpco8+U8AzL03wtG0rBMaYB5aqsvnEZt75+B1ON8bRwhQI9g8SDoYpHjybrB37iUWjNO7cCa1dP6x2kwi5TzxBoCRE/sKFZAwc2M0MvDtTXcdHa6r57MhvfweSnpHGK/8ww3OraysExpgHXtP1JpZ/upz3D75/+9gHnUiXdGYVzCIcDPNU7qM0bSknVhblSkUFXd62c9uG0smbPp1AKET+VxaQHgh0MwtvTlVeZvfq45yqvHV77eR5D/P0C+M8b8sKgTHGN840nmFV1SpKq0r5vCGOFqauAdkDKCkqIRwMU9TSn/r166mLRrl64BNPny+ZmeTNmuUUhblzSMtLbIsLVeXzXzsF4cJn9bzyvRndGufACoExxnfatI09Z/cQqYqwqXbTb7e87sTEgRMJB8OUFJWQc76eWHQdsWiU5qNHPcUgOTnkz3VaXOTNmkVaduIGolFV6s410X9obrfWt0JgjPG1hmsNbKzdSKQq0uVzCO1lpmUy5+E5hINhZo6YSevxWqcolJVxrabGUwxpffuSv2ABgZIQedOnI5nJG3ugI1YIjDEpo6auhtJjpaw+tppzV87Fvd6QnCEsHruYcDBMYaCQ5qNHiUWjxMqitJyKY1i1dtL79yd/4UKnxcWTTyS9xQVYITDGpKDWtlZ2nd5FpCrC1hNb4/6CGWDy4MmEg2EWFS4iLzOPpv37nTOF9etoPd/1vf7tZQweTH7xIvqFQvSZPDlpTzNbITDGpLS65jrWHV9HpCrCoYuH4l6vT3ofFoxeQDgYZuqwqUibcuXjCmLRKPUbNtBaV+cpjsyCAufBtZIQ2RMm9GpRsEJgjDGuysuVRKoirK1ey6Wr8bW4BhiRN4KlwaUsG7uMkfkj0ZYWGnfuJFYWpX7zZtoaGz3FkVVU5D7NHCK7aIzXNDyzQmCMMXdoaWthx8kdRKoibDu5rdNWFneaOmwq4WCYBaMWkJuZS1tzMw3bthEri9Lwi1+gV696iiV74kT3aeYQWSMLvKYSFysExhjTiYtNFymrLiNyLELl5cq418vLzGNh4ULCwTBTBk9BRGhtaKShvJxYNErDjh3gscVFzuTJ7tPMi8gcOsRrKndlhcAYY+Kgqhy+dJhIZYTo8Sixa7G41x0dGM2ysctYMnYJw/KGAdBaV0f95s3EyqI07toFbZ0PW3kbEXKnTXMeXPvqV8gYMMBrOndszgqBMcZ40tzaTPln5USqIuw8tbPTQXPaS5M0ZgyfQTgYZu6ouWSnOw+ZXb9wgdjGjcSiUZoq9ngLJiODvJkznKIwfz7p+fle07FCYIwxPXG28SxrqtdQWlVKTawm7vUCWQGKxxTzbPBZJg26NRBOy+nTxNZvIFZWxtWDBz3FIllZDH7zWwx64w1v61khMMaYnlNVDpw/QKQqwvqa9TS2xH+n0I2OqCVFJTyU89DN+ddqa4mtW0esLEpzZXzfTxS89x6BRQs9xW6FwBhj7rErLVfYcmILpVWl7D6zO+71MiSDp0c+TTgYZvbI2WSm3WpFcfU3v3GKQjRKS+2JDtdPy81l3K92kJaT4yleKwTGGJNAnzd8zuqq1ZQe89YRdWCfgTc7oo4fMP7mfFXl6qHDTouLaJTrZ87c/LvAkiUUvPMjzzFaITDGmF7Qpm1UnKm42RH1amv8zxNMGjSJcDBMaEyIftn9bs7XtjanxcXaMmIbNjD87e+RP2eO59isEBhjTC9ruNbAhpoNRKoi7D+/P+71MtMymTdqHsvGLmPmiJmkp91qWKfXr4NIt5rYWSEwxpgk6klH1CVjlzgdUfsV9igGKwTGGHMfaG1rZefpnTc7ora0xf/U8ZTBUwgHwywsXEjfrL6eP9sKgTHG3GfqmuuIHo9SWlXquSPqN6d8k9cffd3T58VbCNI8bdVbAA+LSLmIHBGRQyLypjv/HRE5KiKfiMgqEemfqBiMMeZ+0i+7Hy898hIrFq/gw6Uf8uqkVxnYZ2CX611tvXrbcwj3WsLOCERkODBcVfeKSD6wBwgDI4GtqnpdRH4IoKrf7WxbdkZgjPGrlrYWtp/cTqQqwvaT2zvsiJqXmUf5C+XkZCTmOYIMT1v1QFVPA6fd9/UicgQoUNWN7RbbBTyfqBiMMeZ+d+OOoXmj5nGx6SJrq9cSqYpQ9UXVzWUWFS7yXAS8SNilofZEpBD4MnDno3h/CKy7yzrfEJEKEak4f/58YgM0xpj7wKCcQbz2pddYuXQlKxav4MUJLxLIChAOhhP6uQn/slhE+gK/BN5W1ZXt5v8d8CTwnHYRhF0aMsakqubWZrLSsro1xGXSLw25QWQCHwL/dUcReA1YDMzvqggYY0wqu9HSOpESVgjEKV/LgSOq+m67+YuA7wLPqOqVRH2+McaY+CTyjOAp4GvApyJy41nrvwX+BcgGNrmnOrtU9U8SGIcxxphOJPKuoR1ARxe1oon6TGOMMd71yl1Dxhhj7l9WCIwxJsU9EL2GROQ8UNvN1R8CLtzDcB4ElnNqsJxTQ09yHq2qg7ta6IEoBD0hIhXx3EfrJ5ZzarCcU0Nv5GyXhowxJsVZITDGmBSXCoXg35IdQBJYzqnBck4NCc/Z998RGGOM6VwqnBEYY4zphBUCY4xJcb4rBCKSLiL7RGStOz1GRHaLSKWI/I+IZCU7xntNRGpE5FMR2S8iFe68gSKyyc17k4gMSHac94qI9BeRD9whT4+IyAyf5zvB3bc3/sRE5Nt+zhlARP7cHeb2oIj8XET6+P14FpE33XwPici33XkJ38++KwTAm8CRdtM/BP5JVccBl4E3khJV4s1V1Snt7jf+a2CLm/cWd9ov/hlYr6qPAJNx9rdv81XVX7v7dgrwBHAFWIWPcxaRAuBbwJOq+iiQDryIj49nEXkU+DowDefnerGIjKMX9rOvCoGIjARKgJ+60wLMAz5wF/kPnHGTU8EynHzBR3mLSACYjdPiHFW9pqpf4NN8OzAfOKaqtfg/5wwgR0QygFycoW/9fDxPxOnGfEVVr+MM6PUsvbCffVUIgPeAvwLa3OlBwBfuPyrASaAgGYElmAIbRWSPiHzDnTfUHTf6xvjRQ5IW3b1VBJwH3ncvAf5URPLwb753ehH4ufvetzmr6ufAPwIncApAHbAHfx/PB4HZIjJIRHKBEPAwvbCffVMIRGQxcE5V97Sf3cGifrxf9ilVfRwoBv5URGYnO6AEygAeB36sql8GGvHRJZHOuNfDlwL/l+xYEs29Dr4MGAOMAPJwfr7v5JvjWVWP4Fz62gSsBw4A1ztd6R7xTSHAGQhnqYjUACtwTiHfA/q7p5YAI4FTyQkvcVT1lPt6Dufa8TTgrIgMB3BfzyUvwnvqJHBSVXe70x/gFAa/5tteMbBXVc+6037OeQFwXFXPq2oLsBKYic+PZ1VdrqqPq+ps4BJQSS/sZ98UAlX9G1UdqaqFOKfPW1X1ZaAceN5d7DWgNEkhJoSI5IlI/o33wFdxTjFX4+QLPspbVc8An4nIBHfWfOAwPs33Di9x67IQ+DvnE8B0Ecl1v+u7sZ/9fjwPcV9HAc/h7O+E72dfPlksInOA76jqYhEpwjlDGAjsA15R1eZkxncvufmtciczgP9W1bdFZBDwv8AonIPq91X1UpLCvKdEZArODQFZQDXwOs5/anyZL4B7zfgzoEhV69x5vt3HACLyFvAHOJdH9gF/hPOdgJ+P5+043222AH+hqlt6Yz/7shAYY4yJn28uDRljjOkeKwTGGJPirBAYY0yKs0JgjDEpzgqBMcakOCsExiSAiPy9iHwn2XEYEw8rBMZ0QRx2rBjfsh9uYzogIoXuWAf/CuwFlotIhdsn/q12y9WIyFsistcdE+KRDrb1dRFZJyI5vZmDMfGyQmDM3U0A/tNtbveX7lgPjwHPiMhj7Za74Db9+zFw2+UgEfkzYAkQVtWmXorbGE+sEBhzd7Wqust9/4KI7MVpa/AlYFK75Va6r3uAwnbzv4bTKO73/NQGwfiPFQJj7q4RnOFOcf6nP19VHwPKgD7tlrvxS74Vp9/TDQdxCsPIhEdqTA9YITCmawGcolAnIkPpuC9+R/YBfwysFpERiQrOmJ6yQmBMF1T1AM4v9UPAvwO/8rDuDpyziTIReSgxERrTM9Z91BhjUpydERhjTIqzQmCMMSnOCoExxqQ4KwTGGJPirBAYY0yKs0JgjDEpzgqBMcakuP8HJbf/KhKdPm8AAAAASUVORK5CYII=\n",
+      "text/plain": [
+       "<Figure size 432x288 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "plt.figure()\n",
+    "for sparsity in sparsityvals:\n",
+    "    errvals=[fdict[rank].inc_approx(sparsity)  for rank in rankvals]\n",
+    "    plt.plot(rankvals,errvals,linewidth=5,label=\"sparsity={:.2f}\".format(sparsity))\n",
+    "    \n",
+    "plt.legend()\n",
+    "plt.xlabel(\"rank\")\n",
+    "plt.ylabel(\"error\")\n",
+    "plt.title(\"error as a function of sparsity\",fontsize=\"xx-large\")\n",
+    "plt.show()\n",
+    "plt.close()\n",
+    "    "
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
diff --git a/ErrorAnalysis/error_as_rank.png b/ErrorAnalysis/error_as_rank.png
new file mode 100644
index 0000000000000000000000000000000000000000..85f1147960897bb98efdaee49731f982328f2a1a
Binary files /dev/null and b/ErrorAnalysis/error_as_rank.png differ