{"nbformat": 3, "nbformat_minor": 0, "metadata": {"name": "", "signature": "sha256:ea5488a343dfb8e7e303aaba95236c270acad097309de7a371845457a049a5c4"}, "worksheets": [{"cells": [{"cell_type": "markdown", "source": ["# Interpolation with Radial Basis Functions"], "metadata": {}}, {"input": ["import numpy as np\n", "import numpy.linalg as la\n", "import matplotlib.pyplot as pt"], "cell_type": "code", "prompt_number": 10, "collapsed": false, "language": "python", "metadata": {}, "outputs": []}, {"input": ["plot_x = np.linspace(-3, 3, 200)"], "cell_type": "code", "prompt_number": 15, "collapsed": false, "language": "python", "metadata": {}, "outputs": []}, {"input": ["np.random.seed(20)\n", "centers = np.random.randn(10)*0.05 + np.linspace(-1.5, 1.5, 10)\n", "centers = np.sort(centers)\n", "centers"], "cell_type": "code", "prompt_number": 147, "collapsed": false, "language": "python", "metadata": {}, "outputs": [{"text": ["array([-1.45580534, -1.15687342, -0.81545651, -0.6171631 , -0.2209083 ,\n", " 0.19465148, 0.54697347, 0.78440928, 1.19182151, 1.52032072])"], "prompt_number": 147, "output_type": "pyout", "metadata": {}}]}, {"input": ["radius = 0.3"], "cell_type": "code", "prompt_number": 176, "collapsed": false, "language": "python", "metadata": {}, "outputs": []}, {"input": ["\n", "def radial_basis_function(x, i):\n", " return np.exp(-(x-centers[i])**2/radius**2)\n", "\n", "pt.plot(plot_x, radial_basis_function(plot_x, 3))"], "cell_type": "code", "prompt_number": 177, "collapsed": false, "language": "python", "metadata": {}, "outputs": [{"text": ["[]"], "prompt_number": 177, "output_type": "pyout", "metadata": {}}, {"text": [""], "png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEACAYAAACuzv3DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGQBJREFUeJzt3XuQXGWZx/Hvk5lkCXKJJJB7MhEjJEDwgiHugg4CGkDE\nS21p0N1VlxLdwlVLd/FSpbGsdb1tYSnKRgstV0QoUbdiKWAKbW9AALlmkwkTk0kmGTIkgQCBBHJ5\n9o+3h0x6ZrpPd59z3r78PlVT9OWdc54myW+eec95zzF3R0REmte42AWIiEh9FOQiIk1OQS4i0uQU\n5CIiTU5BLiLS5BTkIiJNrmKQm9n3zWzQzB4pM+abZtZrZg+Z2avSLVFERMpJ0pH/AFg61ptmdjHw\ncnefD3wQuC6l2kREJIGKQe7ufwSeLDPkrcAPi2NXA5PMbGo65YmISCVpzJHPBPqHPd8KzEphuyIi\nkkBaBzut5LnW/YuI5KQzhW1sA2YPez6r+NoRzEzhLiJSA3cvbZaPkEZHvhL4RwAzWwLsdvfBMYpp\n2a/Pf/7z0Wtohs+3caMzZYqzZo1z223O9OnOc8+1xmdr1C99vub+SqJiR25mPwHeAEwxs37g88D4\nYjCvcPdfm9nFZrYBeBZ4f6I9S1u6+WZ497vhtNPC17x58Oc/wwUXxK5MpHlVDHJ3X5ZgzFXplCOt\nrlCAD33o8PM3vhHuuENBLlIPrexMSXd3d+wSMpXG59u/H+68E17/+sOvnX8+/Pa3dW+6Lvqza26t\n/vmSsKRzMHXvyMzz2pc0prvvDt34gw8efm3fPjjxROjvh0mT4tUm0qjMDM/hYKdIIoUCnHfeka8d\ndRQsWQK//32UkkRagoJcclMowGi/Bb/xjfGnV0SamYJccnHwYJgfP/fcke+dc06YdhGR2ijIJReb\nNsEJJ4SvUgsXQk8P6BCKSG0U5JKLnh5YsGD09yZPDnPlAwP51iTSKhTkkot16+DUU8d+f8GCMEZE\nqqcgl1z09JQP8oULYe3a/OoRaSUKcslFuakVUEcuUg8FuWTOXVMrIllSkEvmduwI/z3xxLHHKMhF\naqcgl8wNTatYmUXGM2aE5fq7duVXl0irUJBL5ipNq0AIeXXlIrVRkEvmKp2xMmTBgjBWRKqjIJfM\nPfoonHJK5XHz5oUVoCJSHQW5ZK6vD7q6Ko/r6oLNmzMuRqQFKcglU+6wZQvMmVN57Ny5CnKRWijI\nJVNPPgnjxiW7aYSCXKQ2CnLJ1ObNIaCTmDkTBgfDLeFEJDkFuWRq8+Zk0yoA48fD1KmwbVu2NYm0\nGgW5ZGrLluQdOYSxfX2ZlSPSkhTkkqlqplZA8+QitVCQS6aqDXKdgihSPQW5ZKqaOXJQRy5SCwW5\nZKqWOXIFuUh1FOSSmb174amnYNq05N+jIBepnoJcMrNlC8yaFRYEJTVnDvT3w6FD2dUl0moU5JKZ\naufHASZOhOOPDwuDRCQZBblkpr+/+iCHsMJzYCD9ekRalYJcMjMwEEK5WjNmKMhFqqEgl8wMDIRQ\nrpaCXKQ6CnLJTD1BruutiCSnIJfM1BrkmiMXqY6CXDKjqRWRfCjIJRMHD8Ljj1e3GGiIglykOhWD\n3MyWmlmPmfWa2dWjvD/FzG4zswfNbI2ZvS+TSqWpPP44nHBCuMZ4tRTkItUpG+Rm1gFcCywFFgLL\nzGxBybCrgAfc/ZVAN/BfZtaZQa3SRGqdVgE48UTYvRuefz7dmkRaVaWOfDGwwd373H0/cBNwWcmY\nx4Djio+PA3a5+4F0y5RmU0+QjxsX7hS0fXu6NYm0qkpBPhPoH/Z8a/G14b4HnGZmA8BDwEfTK0+a\nVT1BDjpzRaQalaZAPME2PgM86O7dZnYysMrMznT3Z0oHLl++/MXH3d3ddHd3V1GqNJN6g1zz5NKu\nCoUChUKhqu+pFOTbgNnDns8mdOXD/S3wHwDu/lcz2wScAtxXurHhQS6tbWAAzjqr9u9XkEu7Km1y\nv/CFL1T8nkpTK/cB882sy8wmAO8CVpaM6QEuADCzqYQQ35i4amlJaXTkWt0pkkzZjtzdD5jZVcDt\nQAdwvbuvM7Mri++vAL4E/MDMHiL8YPh3d38i47qlwdV6wawhM2ZAT0969Yi0soqnCbr7rcCtJa+t\nGPZ4J3Bp+qVJM1NHLpIfreyU1O3fD088Ec4Hr9X06Tr9UCQpBbmkbnAwhHhHR+3bmDpVdwkSSUpB\nLqkbHAxBXI/Jk8ONm/fvT6cmkVamIJfUpRHkHR0hzHfsSKcmkVamIJfUpRHkoOkVkaQU5JI6BblI\nvhTkkjoFuUi+FOSSOgW5SL4U5JK6tIJ82jSdSy6ShIJcUqeOXCRfCnJJnYJcJF8KcknVgQPhNm1T\nptS/LQW5SDIKcknVjh3hpsv1LM8foiAXSUZBLqlKa1oFQlf/5JOhyxeRsSnIJVVpBnlnJ7z0pbBz\nZzrbE2lVCnJJVZpBDppeEUlCQS6pSjvIdS65SGUKckmVOnKR/CnIJVUKcpH8KcglVQpykfwpyCVV\nCnKR/CnIJVUKcpH8KcglNQcPwq5d4cbLaVGQi1SmIJfU7NwJxx8P48ent00FuUhlCnJJTdrTKhC6\n+127QrcvIqNTkEtqBgfDAp40jR8PkyZpmb5IOQpySU0WHTloekWkEgW5pEZBLhKHglxSoyAXiUNB\nLqlRkIvEoSCX1CjIReJQkEtqsgryadMU5CLlKMglNVl25LomucjYFOSSikOHwo2XTzop/W1rakWk\nPAW5pOKJJ+DYY2HChPS3rSAXKa9ikJvZUjPrMbNeM7t6jDHdZvaAma0xs0LqVUrDy2paBUKXv3Nn\n6PpFZKTOcm+aWQdwLXABsA2418xWuvu6YWMmAd8G3uzuW81sSpYFS2PKMsjHj4fjjkv/yooiraJS\nR74Y2ODufe6+H7gJuKxkzOXAz9x9K4C766oYbSjLIAdNr4iUUynIZwL9w55vLb423HzgBDP7nZnd\nZ2b/kGaB0hwU5CLxlJ1aATzBNsYDrwbOB44G7jKzu929t97ipHkoyEXiqRTk24DZw57PJnTlw/UD\nO919L7DXzP4AnAmMCPLly5e/+Li7u5vu7u7qK5aGNDgIJ5+c3fYV5NIuCoUChUKhqu8x97GbbjPr\nBNYTuu0B4B5gWcnBzlMJB0TfDPwNsBp4l7uvLdmWl9uXNLdLLoEPfQguvTSb7X/pS/D00/DlL2ez\nfZFGZWa4u5UbU7Yjd/cDZnYVcDvQAVzv7uvM7Mri+yvcvcfMbgMeBg4B3ysNcWl9eUyt9GqyTmRU\nlaZWcPdbgVtLXltR8vzrwNfTLU2aSR5B/vjj2W1fpJlpZafUzT2EbJZBftJJmiMXGYuCXOq2ezcc\ndVT4yooOdoqMTUEudcu6G4fDUys6Xi4ykoJc6pb1/Dgc7vh37852PyLNSEEudcsjyEHTKyJjUZBL\n3RTkInEpyKVuCnKRuBTkUrfBwWzuDFRKQS4yOgW51C3PjlyLgkRGUpBL3fIKci0KEhmdglzqpjly\nkbgU5FIXdwW5SGwKcqnLnj1gBscck/2+FOQio1OQS13y6sbhcJBrmb7IkRTkUpc8g/yYY0L3v2dP\nPvsTaRYKcqlLnkEOml4RGY2CXOqiIBeJT0Eudck7yE86SYuCREopyKUu6shF4lOQS13yuKnEcApy\nkZEU5FIXdeQi8SnIpS4KcpH4FORSFwW5SHwKcqnZ3r3wwgtw3HH57VNBLjKSglxqNtSNm+W3TwW5\nyEgKcqlZ3tMqAMcfD88/H34bEJFAQS41ixHkZloUJFJKQS41ixHkoOkVkVIKcqmZglykMSjIpWYK\ncpHGoCCXminIRRqDglxqFivIp02Dxx7Lf78ijUpBLjXbvj2Eat6mT1eQiwynIJeaPfZYCNW8KchF\njqQgl5rs2QMHDuS7PH+IglzkSApyqclQN57n8vwhQ0Hunv++RRpRxSA3s6Vm1mNmvWZ2dZlxrzWz\nA2b2jnRLlEYUa1oF4CUvgQkTYPfuOPsXaTRlg9zMOoBrgaXAQmCZmS0YY9xXgNuACD2a5C1mkIOm\nV0SGq9SRLwY2uHufu+8HbgIuG2XcR4BbgB0p1ycNSkEu0jgqBflMoH/Y863F115kZjMJ4X5d8SXN\nXLYBBblI46gU5ElC+RvAp9zdCdMqmlppAwpykcbRWeH9bcDsYc9nE7ry4V4D3GTh9IUpwEVmtt/d\nV5ZubPny5S8+7u7upru7u/qKpSEoyEWyUSgUKBQKVX2PeZlzuMysE1gPnA8MAPcAy9x93RjjfwD8\n0t1/Psp7Xm5f0lxOPx1uvBEWLYqz/xtugF/9Cn7ykzj7F8mLmeHuZWc6ynbk7n7AzK4Cbgc6gOvd\nfZ2ZXVl8f0Vq1UpTGRiAGTPi7V8duchhZTvyVHekjrxl7NsXbrm2b1+cBUEAa9fC298O69fH2b9I\nXpJ05FrZKVXbvj3/my6XUkcucpiCXKoW+0AnwKRJ8MIL8OyzcesQaQQKcqlaIwS5mbpykSEKcqla\nIwQ5wMyZ4aCrSLtTkEvVYp+xMmTmTNi2LXYVIvEpyKVqW7fCrFmxqwg1bC1dnibShhTkUrVGCXJ1\n5CKBglyq1ihBro5cJFCQS1XcFeQijUZBLlXZvRs6O+HYY2NXoqkVkSEKcqnK1q0hQBvB9OkwOAgH\nD8auRCQuBblUZdu2xphWgXDfzsmTQ5iLtDMFuVSlUebHh8ycqXlyEQW5VKXRglwHPEUU5FKlRgty\nHfAUUZBLlRotyNWRiyjIpUqNdNYKhCBXRy7tTkEuVWm0jlwHO0UU5FKFPXvg+efhhBNiV3LYrFnQ\n3x+7CpG4FOSS2FA3HvMWb6Vmzw5TK4cOxa5EJB4FuSS2eTPMnRu7iiNNnBhu+7Z9e+xKROJRkEti\nfX3Q1RW7ipG6ukJtIu1KQS6JNWqQz50bflsQaVcKckls8+bGDHJ15NLuFOSSWF9f482RgzpyEQW5\nJNaoUyvqyKXdKcglkX37YNcumDEjdiUjqSOXdqcgl0T6+8M55B0dsSsZaSjI3WNXIhKHglwSadT5\ncYBjjoGjj4YdO2JXIhKHglwSadT58SFz52qeXNqXglwSadRTD4d0dWmeXNqXglwSUUcu0rgU5JJI\nI8+RA8ybB5s2xa5CJA4FuSSycSO87GWxqxjby18Ovb2xqxCJQ0EuFT3zDOze3Vh3Bio1f76CXNqX\nglwq2rABTj4ZxjXw35aurnAp2337Ylcikr9E/zTNbKmZ9ZhZr5ldPcr77zGzh8zsYTP7s5ktSr9U\niaW3N3S8jayzE+bM0Ty5tKeKQW5mHcC1wFJgIbDMzBaUDNsIvN7dFwFfBL6bdqEST28vvOIVsauo\nTPPk0q6SdOSLgQ3u3ufu+4GbgMuGD3D3u9z9qeLT1UAD3Z5X6vXoo43fkYPmyaV9JQnymcDw29tu\nLb42ln8Gfl1PUdJYmmFqBRTk0r46E4xJfCkiMzsP+ADwd6O9v3z58hcfd3d3093dnXTTElEzBfkv\nfhG7CpH6FAoFCoVCVd9jXuGScWa2BFju7kuLzz8NHHL3r5SMWwT8HFjq7htG2Y5X2pc0nt27w53q\nn34azGJXU97GjXDeeVqqL63FzHD3sv/6kkyt3AfMN7MuM5sAvAtYWbKjOYQQf+9oIS7Na6gbb/QQ\nh3DWyuCgTkGU9lMxyN39AHAVcDuwFrjZ3deZ2ZVmdmVx2OeAlwLXmdkDZnZPZhVLrprlQCeEUxDn\nzoW//jV2JSL5SjJHjrvfCtxa8tqKYY+vAK5ItzRpBM0yPz5kwQJYtw5OOy12JSL5aeC1etII1qxp\nrlA8/fRQs0g7UZBLWY88AmecEbuK5E4/PdQs0k4U5DKmvXthyxY45ZTYlSR3xhnqyKX9KMhlTGvX\nhvnx8eNjV5Lc/Pnhh8/evbErEcmPglzG1GzTKgATJoRrrqxbF7sSkfwoyGVMa9Y0X5CDDnhK+1GQ\ny5iasSMHzZNL+1GQy5geeSR0t81GHbm0GwW5jGrXLnj22bDsvdmcfjo8/HDsKkTyoyCXUT3wACxa\n1BzXWCk1bx489xw89ljsSkTyoSCXUa1eDWefHbuK2pjB4sVwj674I21CQS6jWr0aliyJXUXtzj47\nfAaRdqAglxHc4e67m7cjB3Xk0l4U5DJCX19YzTmrie+8ungx3HsvHDoUuxKR7CnIZYShbrwZD3QO\nOfFEmDIF1q+PXYlI9hTkMkKzz48PWbxY8+TSHhTkMkKzz48PWbIE7rordhUi2VOQyxGefjqsinzt\na2NXUr/zzoM77ohdhUj2FORyhEIhdLJHHx27kvqdcQbs2QMbN8auRCRbCnI5wm9+A296U+wq0mEG\nF1wAq1bFrkQkWwpyOcJvfgNvfnPsKtLzpjcpyKX1mbvnsyMzz2tfUptNm+B1r4OBARjXIj/iBwbC\nRbR27ICOjtjViFTPzHD3sicDt8g/V0nDqlVw4YWtE+IAM2aEr3vvjV2JSHZa6J+s1OsXv4BLLold\nRfre9ja45ZbYVYhkR1MrAsDgIJxyCmzbBi95Sexq0rV2bZgr37KltX7bkPagqRVJ7Oab4dJLWy/E\nARYuhMmT4Y9/jF2JSDYU5ALAj38M73lP7Cqyc/nlcOONsasQyYamVoRHH4Vzzw3TKp2dsavJxubN\n8JrXhOmVVljsJO1DUyuSyDXXwBVXtG6IA8ydC+ecA9//fuxKRNKnjrzNbd8e5pB7euCkk2JXk627\n7gpTLL29rf1DS1qLOnKp6Jpr4L3vbf0Qh7DYac4c+OlPY1ciki515G1syxZ49avh/vtDwLWDVavg\nwx+GRx6BiRNjVyNSmTpyGZN7CLSPfax9QhzCytWzzoIvfjF2JSLpUUfepm64Ab76VbjvPpgwIXY1\n+dq+HRYtgltvDWeyiDQydeQyqjvvhI9/HH74w/YLcYBp0+C668LS/S1bYlcjUr+KQW5mS82sx8x6\nzezqMcZ8s/j+Q2b2qvTLlLTcfz+84x3wox/Bq9r4T+qd7ww/zC66CPr7Y1cjUp+yQW5mHcC1wFJg\nIbDMzBaUjLkYeLm7zwc+CFyXUa0NrVAoxC6hLHe4/vpwrfHvfAeWLq3u+xv989Xi4x+H978fzjyz\n0NK3hGvFP7vhWv3zJVGpI18MbHD3PnffD9wEXFYy5q3ADwHcfTUwycympl5pg2vUv0wHD8Ltt4dT\n7771Lfj970NHXq1G/Xz1MINPfhIuuaTABz4Q/r/cc0/4oddKWvHPbrhW/3xJVFoWMRMY/ovnVqD0\n/uqjjZkFDNZdnVTl4MGwzH7z5rDo5c47wwG96dPD2SmXX66r/43m5JPhu9+Fb387/D/q7AxXSzz7\nbHjZy2DePJg6NQS/SCOqFORJe5PSv+Kjft9b3jJsQMmIZnte+lpfH/z2t/nXtHdvuMHwnj3w1FMw\nZQp0dYUAWrIkdJynnjqydjnSxInh/9UnPhFuQlEowMqV4a5JmzbBM8/AcceF67QcfXS4SuTEieGu\nQ2bVfeVt/Xr4y1/y329eWv3zJVH29EMzWwIsd/elxeefBg65+1eGjflvoODuNxWf9wBvcPfBkm21\n2C+sIiL5qHT6YaWO/D5gvpl1AQPAu4BlJWNWAlcBNxWDf3dpiCcpREREalM2yN39gJldBdwOdADX\nu/s6M7uy+P4Kd/+1mV1sZhuAZ4H3Z161iIi8KLeVnSIiko1cz2Ewsy8WFw09aGZ3mNnsPPefJTP7\nmpmtK36+n5vZ8bFrSpOZ/b2Z/Z+ZHTSzV8euJy1JFrw1KzP7vpkNmtkjsWvJgpnNNrPfFf9erjGz\nf41dU1rM7CgzW13MyrVm9p9lx+fZkZvZse7+TPHxR4Az3f2K3ArIkJldCNzh7ofM7MsA7v6pyGWl\nxsxOBQ4BK4BPuPv9kUuqW3HB23rgAmAbcC+wzN3XRS0sJWZ2LrAH+B93PyN2PWkzs2nANHd/0MyO\nAf4CvK2F/vyOdvfnzKwT+BPwSXf/02hjc+3Ih0K86BhgZ577z5K7r3L3Q8Wnqwnn0rcMd+9x90dj\n15GyJAvempa7/xF4MnYdWXH37e7+YPHxHmAdMCNuVelx9+eKDycQjlE+MdbY3JeHmNl/mNkW4J+A\nL+e9/5x8APh17CKkotEWs82MVIvUoXhm3asITVRLMLNxZvYgYXHl79x97VhjU7/hlZmtAqaN8tZn\n3P2X7v5Z4LNm9ingGproLJdKn6045rPAC+7edPdsT/L5WoyO9LeA4rTKLcBHi515Syj+hv/K4vG2\n282s290Lo41NPcjd/cKEQ2+kybrWSp/NzN4HXAycn0tBKaviz65VbAOGH3CfTejKpUmY2XjgZ8AN\n7v6/sevJgrs/ZWa/As4CCqONyfuslfnDnl4GPJDn/rNkZkuBfwMuc/d9sevJWKss7npxwZuZTSAs\neFsZuSZJyMwMuB5Y6+7fiF1PmsxsiplNKj6eCFxImbzM+6yVW4BTgIPAX4EPu/vjuRWQITPrJRyU\nGDogcZe7/0vEklJlZm8HvglMAZ4CHnD3i+JWVT8zuwj4BocXvJU9zauZmNlPgDcAk4HHgc+5+w/i\nVpUeMzsH+APwMIenyT7t7rfFqyodZnYG4aqy44pfP3L3r405XguCRESamy5qKiLS5BTkIiJNTkEu\nItLkFOQiIk1OQS4i0uQU5CIiTU5BLiLS5BTkIiJN7v8B8bTeRRUyWhAAAAAASUVORK5CYII=\n", "output_type": "display_data", "metadata": {}}]}, {"input": ["def f(x): return x**3 - 3*x\n", "\n", "pt.plot(plot_x, f(plot_x))"], "cell_type": "code", "prompt_number": 178, "collapsed": false, "language": "python", "metadata": {}, "outputs": [{"text": ["[]"], "prompt_number": 178, "output_type": "pyout", "metadata": {}}, {"text": [""], "png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEACAYAAACnJV25AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHUBJREFUeJzt3XmYVOWVx/HvEVSi4jLiABIScCQqyqBGsWPMUAEXXKIS\n425coxgS1Gg0qIm0OjoqYghkXAMqLhiU6IiKYZGKKyiryK4BQVkUZBEFm4Yzf7zVdttpmuqu5Vbd\n+n2e5z5dy616Tz3iqVPvfRdzd0REJD62izoAERHJLiV2EZGYUWIXEYkZJXYRkZhRYhcRiRkldhGR\nmMkosZtZWzObYGazzOw9M7si9fi/mdlYM5tvZmPMbPfshCsiIttimYxjN7NWQCt3n25muwBTgFOB\ni4CV7n6Xmf0O2MPd+2YlYhERqVdGFbu7L3f36anb64E5QBvgZODR1GmPEpK9iIjkQdb62M2sHXAI\nMAlo6e4rUk+tAFpmqx0REalfVhJ7qhtmJHClu39e8zkPfT1at0BEJE+aZvoGZrY9Iak/5u7PpR5e\nYWat3H25mbUGPqnjdUr2IiKN4O5W3/OZjooxYAgw290H1njqeeCC1O0LgOdqvzYVXGyPfv36RR6D\nPp8+Xyl+vjh/Nvf06uFMK/YfAucB75rZtNRj1wN3ACPM7BJgEXBGhu2IiEiaMkrs7v46W6/6j87k\nvUVEpHE08zRHEolE1CHklD5fcYvz54vzZ0tXRhOUMmrYzKNqW0SkWJkZnsuLpyIiUniU2EVEYkaJ\nXUQkZpTYRURiRoldRKRIvPdeeucpsYuIFInBg9M7T8MdRUSKwOrVsM8+sGaNhjuKiMTC0KFw4onp\nnauKXUSkwG3eDB06wPDhUFamil1EpOi99BK0aAFHHJHe+UrsIiIFbtAg6NMn/fPVFSMiUsDmzIEf\n/xg+/BB23FFrxYiIFL3Bg+HSS0NST5cqdhGRArVqFey7L8yeDa1bh8dUsYuIFLH77oOePauTerpU\nsYuIFKCNG6FdOxg/Hg48sPpxVewiIkXq8cfh0EO/mdTTlXFiN7OhZrbCzGbWeKzczD4ys2mpo0em\n7YiIlIotW+Cee+Caaxr3+mxU7A8DtRO3A/e4+yGp4+UstCMiUhJGjw6jYLp1a9zrM07s7v4asLqO\np+rtAxIRkbrdfTf89rdgjcyiuexj72NmM8xsiJntnsN2RERiY8oUeP99OOOMxr9H0+yF8w33Abek\nbt8KDAAuqX1SeXn517cTiQSJRCJH4YiIFIcBA+DKK2H77cP9ZDJJMpls0HtkZbijmbUDRrl7p3Sf\n03BHEZFvWrwYDj4YFi6E3Xar+5zIhjuaWc3h9D2BmVs7V0REgoED4aKLtp7U05VxxW5mw4GuQAtg\nBdAPSAAHE0bHLAR6ufuKWq9TxS4ikrJqVVhzfcYMaNt26+elU7Fr5qmISAG46SZYtgweeqj+85TY\nRUSKwLp1YT/TiRPDol/10ZICIiJF4N574dhjt53U06WKXUQkQl9+Gar1cePgoIO2fb4qdhGRAjdk\nCJSVpZfU06WKXUQkIhUVoftl5Eg4/PD0XqOKXUSkgA0dCh07pp/U06WKXUQkAl99FcatP/MMdOmS\n/utUsYuIFKghQ6BTp4Yl9XSpYhcRybONG0Pf+rPPNrwbRhW7iEgBeughOOSQ7PetV1HFLiKSRxs2\nhGp91Kiwp2lDqWIXESkwDzwQKvXGJPV0qWIXEcmTzz8P1fq4ceHCaWOoYhcRKSD33APHHdf4pJ4u\nVewiInnw6adwwAHwzjvQvn3j30fL9oqIFIjf/AYqK2Hw4MzeR4ldRKQAfPhhuFg6eza0bJnZe6mP\nXUSkAJSXQ+/emSf1dDXNTzMiIqVp+nQYPRrmzctfmxlX7GY21MxWmNnMGo/9m5mNNbP5ZjbGzHbP\ntB0RkWLjDtdcE/Yz3W23/LWbja6Yh4EetR7rC4x19+8B41P3RURKyosvhg2qL7ssv+1mnNjd/TVg\nda2HTwYeTd1+FDg103ZERIrJpk1w7bXQvz80zXOnd64unrZ09xWp2yuAPF0yEBEpDA89BG3awAkn\n5L/tnH+PuLubWZ3jGsvLy7++nUgkSCQSuQ5HRCTn1qyBm2+GMWPA6h2YuG3JZJJkMtmg12RlHLuZ\ntQNGuXun1P25QMLdl5tZa2CCu+9f6zUaxy4isXTVVfDFF6Fqz7Z0xrHnqmJ/HrgAuDP197kctSMi\nUlDeew+eeCJMRopKxhW7mQ0HugItCP3pNwH/B4wAvgMsAs5w9zW1XqeKXURixR26d4ef/hR+/evc\ntKElBURE8ujpp+HWW2Hq1NyNhFFiFxHJky++gI4dYdgw6No1d+1orRgRkTy5/XY48sjcJvV0qWIX\nEcnQrFmQSMCMGbD33rltSxW7iEiObdkCl18exq3nOqmnS4ldRCQDQ4eG5QN69Yo6kmrqihERaaQV\nK8L+pWPHQufO+WlTo2JERHLo3HPDejB33ZW/NqOceSoiEmsvvAATJ8K770Ydyb9SYhcRaaA1a8IF\n08cfh513jjqaf6WuGBGRBrr4YmjWDO69N/9tqytGRCTLRo+GCRMKswumihK7iEia1q4NwxqHDoXm\nzaOOZuvUFSMikqbzz4eddoL7748uBnXFiIhkyYgRYRTMtGlRR7JtqthFRLbho4/g+98PQxwPPzza\nWLRWjIhIhrZsgQsvhD59ok/q6VJiFxGpx8CBsGED9O0bdSTpU1eMiMhWvP02nHQSTJoE7dtHHU2g\nrhgRkUZaswbOOgseeKBwknq6clqxm9kiYB2wGdjk7l1qPKeKXUQKkjucdhq0bQt/+lPU0XxTIQx3\ndCDh7p/luB0RkawZPBgWL4bhw6OOpHHyMY693m8WEZFC8uabcNtt8NZbsOOOUUfTOLnuY3dgnJlN\nNrNLc9yWiEhGli6FM84ISwbss0/U0TReriv2H7r7MjPbCxhrZnPd/bWqJ8vLy78+MZFIkEgkchyO\niEjdKirgZz8Ly/GeeGLU0VRLJpMkk8kGvSZvwx3NrB+w3t0HpO7r4qmIFIxf/hKWL4eRI2G7Ah4v\nGOlwRzPbycyap27vDBwLzMxVeyIijXX//ZBMwqOPFnZST1cuu2JaAs+aWVU7T7j7mBy2JyLSYGPG\nQHk5vP467Lpr1NFkh2aeikjJmjULfvzj0P3yox9FHU16NPNURGQrVqwIywXcc0/xJPV0KbGLSMlZ\nvx5OPjlsnHHeeVFHk33qihGRklJRAT/5SVgu4KGHwIpsCmU6XTFK7CJSMrZsgXPPhY0b4emnoWkR\n7iFXCGvFiIgUBHe46qowu/Tvfy/OpJ6uGH80EZHAHa6/Hl57DSZMgGbNoo4ot5TYRSTW3OEPf4DR\no+GVV2D33aOOKPeU2EUk1m65BZ57LlTqe+4ZdTT5ocQuIrHkDjffDH/9a1guYK+9oo4of5TYRSR2\n3OHqq0OVnkxCy5ZRR5RfSuwiEiubN0OvXjB7dkjse+wRdUT5p8Qu/+Krr+Djj+Gzz2Dduurj88+h\nsjKcY1Z9fOtbYfGkXXeF5s3D3xYtwk/fJk2i/SxSWjZuhJ//HFavDot77bJL1BFFQ4m9RK1fHyqa\n994LCyEtXAhLloTjs89g771Dcq6ZrJs3h+23Dz9zax4bNoSkX/UFsHYtfPpp2OW9ZcvwXnvvDd/9\nLuy7L3ToEI7vfjfeY4klv1atglNOgW9/G154If5DGuujmaclwD0k79dfh0mTwvHhh7DffnDQQXDg\ngWEbsLZt4TvfCck4G5V2RUXYuGDp0vAL4MMPYcECeP/98Hf58tDuf/5n9dGpU4ih2KZ5S7Q++ACO\nPx5++lO4/fZ4rKm+NVpSoIStXAkvvghjx8L48bDTTvBf/wVlZXDEESGhR10tf/UVzJsH775bfcyY\nEb4QjjgiHGVl0KVLafaTSnpefRXOPDOsqd6rV9TR5J4Se4lZtgxGjAhjdqdOhe7dQxXTvXtxbcy7\ndGn1L4uJE2HKlNCVU1YWlldNJOA//kNVfalzhz//Gf77v+Gxx+DYY6OOKD+U2EvAhg0hkQ8bFpLg\nKaeEn6PHHBMuasZBZWW4HjBxYvWUcAgJvupQoi8tGzaEPUqnTg3//oupcMmUEnuMLVwI994LDz8M\n3/8+XHABnHpq6HKJO/fQp5pMVh8QvsyOOy78LZUZhqVo3rzQ9bL//jBkCOy8c9QR5ZcSewy9+Sbc\neSe88QZcdFGoWkqpWqmLe7gYO2ZMWLXvH/8I/9Mfd1w4ysqiv54gmXOHRx6B664L3S+XXVaav9Ii\nTexm1gMYCDQB/uLud9Z6Xok9Te7hAuhtt8GiReEf9gUXlEZ13hgVFeGLryrRL1wY9rU8/ng48cTQ\nXy/F5bPP4Ne/DhfX//rXcPG/VEWW2M2sCTAPOBr4GHgHONvd59Q4R4k9DW+9Bb/7HXzyCdxwA5x9\ndhhLLulbsSKMDnrpJXj55fAL56STwnHoofEeGhcHzz8ffpmedhrccYcKmigT+w+Afu7eI3W/L4C7\n31HjHCX2esybFxL522+H1enOP1+zOLOhsjJ0Z40aFSaxrFkTqviTToKjjy7dmYqFaOXKsDHGxIkw\ndGgYrivpJfZc1SptgCU17n+Ueky2Yf166NsXfvjDMH57/vzQl66knh1Nm4YE0b8/zJkTRtl06gT/\n+7+hi6ZHjzCEbtGiqCMtXZs3w/33Q8eOYVmKGTOU1BsqV5eU0irFy8vLv76dSCRIJBI5CqfwucPI\nkWFFuq5dYeZMaN066qjib9994corw7FuXeiyeeEFuPXWkFSqumx0ATY/Jk0KfenNmsG4cWE2cqlL\nJpMkq4Z+pSlXXTFlQHmNrpjrgS01L6CqK6ba8uWhD3H+fLjvPlUnhWDLFnjnnZDkR40Ka+gcdxyc\ncEL4W0pre+fD/Plw443hmtLtt4eFvEpxxEs6ouyKmQx0MLN2ZrYDcCbwfI7aKlru8OST0Llz+Nk5\ndaqSeqHYbruwpMGtt8L06WG5g27d4G9/CwuYlZWFax+TJ4cvAWmcDz+Eyy8PXY+HHRYS/PnnK6ln\nKpfDHY+nerjjEHf/n1rPl3TFvmZNGIc7e3YYm3vYYVFHJOmqqAgLqr30UjhWrQpDKU84IUxrL4U9\nNTM1b14Y4fL883DppXDttZpUli5NUCpQb70F55wT+m779y/t5UXj4J//DBslv/RSuBh78MEhwXfr\nBocfruGpVdzDZtKDB4d5Bn36hEMLvDWMEnuBcQ+zRv/4R3jwwbCui8TLhg1hiYNx48KksoULw8Jl\n3bqFxdg6dSq9cfMrV8Lw4eH6UZMm4eLouedqaGljKbEXkM8/hwsvDOuSP/10WPtc4m/lyrBo2Suv\nhES/enWYBZtIhH7lgw6K51DWDRvCr5hhw8IX3Yknhi6Xrl3Vf54pJfYCsWBBWKDryCPDGOkdd4w6\nIonKkiUhySeToUtu2bLQXXPkkeEoKyvePvqVK0Myf/bZ8CV22GFw3nlhxuiuu0YdXXwosReA8ePD\nMgC33BI2AVC1IjWtWhVmVr75ZjgmTw4TpTp3Dn31nTuHo02bwvu389lnYdz5hAmh6+mDD8IvkZ49\n4Sc/0cXQXFFij9gjj4R1XkaMCD9BRbalshLmzg2zLadPr/67ZUtI8Pvv/819Y9u3hx12yG1M7vDR\nR9X7486YERL68uVhyehEIizH0KWLLhTngxJ7RNzDNl2PPRa2pzvggKgjkmLmHpLojBlhnHfNfWOX\nLAmbjldtGN66dfi7557Vm5BXHVVfAFWVvxls3Bhm3NbcjHzFCli8OLz34sUhqe+2W/X+uJ06hTH+\nHTvG8/pAoVNij8DmzWF8+syZYcZiy5ZRRyRxtmlTSMRLl37zWL26OlFXHZs2hS8JCH/dw1Dbmsm/\neXP493+v3ti8bdtwlNpmFoVMiT3PKirCMK61a8MFJP3PICLZFuWSAiXnyy/DuPTKylCpK6mLSFSU\n2LNg/fowpXyvvcIYdQ1nFJEoKbFn6IsvwtIAHTqEUTBa2lVEoqY+9gxs2BCSetu2YYeXUpsqLiL5\npz72HNq4Mcwmbd0ahgxRUheRwqGKvREqK+H008NkjCefVPeLiORPOhW7UlIDuYfdjr74Iuyuo6Qu\nIoVGaamBfv/7MANw/PjcT+UWEWkMJfYGGDQInnkm7J7TvHnU0YiI1E2JPU2jRoWtvN58UxsZi0hh\nU2JPw7RpcPHFoU+9XbuooxERqV9OBumZWbmZfWRm01JHj1y0kw8ffwwnnxy29TriiKijERHZtlxV\n7A7c4+735Oj98+LLL0NS/9Wv4Gc/izoaEZH05HJaTYHt99Iw7mGPxo4dw2YZIiLFIpeJvY+ZzTCz\nIWZWdLs4DhwIc+bAAw8U3pZkIiL1aXRXjJmNBVrV8dSNwH3ALan7twIDgEtqn1heXv717UQiQSKR\naGw4WfXKK3DnnWEvyp12ijoaESllyWSSZDLZoNfkfEkBM2sHjHL3TrUeL8glBRYvDns3PvEEdO8e\ndTQiIt8U2SJgZta6xt2ewMxctJNtmzbBmWfC1VcrqYtI8cpJxW5mw4CDCaNjFgK93H1FrXMKrmK/\n9lqYPTtMRtJqjSJSiLTnaQO88AL07g1Tp4Zd30VECpFWd0zTkiVwySUwcqSSuogUv5LvcNi8Gc45\nB37zGzjqqKijERHJXMkn9jvvDBtmXHdd1JGIiGRHSfexT54MJ5wAU6aEfUtFRAqd9jytx5dfwnnn\nhTXWldRFJE5KtmL/1a9g7Vp4/PHIQhARaTCNitmKMWPC8MYZM6KOREQk+0quK2bdurBq44MPwu5F\ntzSZiMi2lVxXzOWXQ2Ul/OUveW9aRCRj6oqpZfx4ePFFeO+9qCMREcmdkumKWb8efvGLsL76brtF\nHY2ISO6UTFfMVVfB6tXw6KN5a1JEJOvUFZMyeTI89RTMmhV1JCIiuRf7rpjKSujVKywdsOeeUUcj\nIpJ7sU/sf/4z7LornH9+1JGIiORHrPvYlyyBQw6BN96A/fbLaVMiInlR8mvFXHEF9OmjpC4ipSW2\nF09ffhlmzoThw6OOREQkv2JZsVdUwJVXwsCB0KxZ1NGIiORXoxO7mZ1uZrPMbLOZHVrruevNbIGZ\nzTWzYzMPs2EGDoQOHeCkk/LdsohI9DLpipkJ9AQeqPmgmXUEzgQ6Am2AcWb2PXffkkFbaVu6FO66\nCyZOzEdrIiKFp9EVu7vPdff5dTx1CjDc3Te5+yLgfaBLY9tpqOuug8sug333zVeLIiKFJRcXT/cG\natbLHxEq95x78034xz9gzpx8tCYiUpjqTexmNhZoVcdTN7j7qAa0U+eA9fLy8q9vJxIJEolEA96y\nVgMOV18Nt98Ou+zS6LcRESkoyWSSZDLZoNdkPEHJzCYA17j71NT9vgDufkfq/stAP3efVOt1WZ2g\n9NRTcPfd8PbbsF0sx/qIiOR3glLNRp4HzjKzHcysPdABeDtL7dRp40bo2xcGDFBSFxHJZLhjTzNb\nApQBL5rZaAB3nw2MAGYDo4HeuV474E9/goMPhq5dc9mKiEhxKPq1Yj79FA44IFw4/d73shCYiEgB\nS6crpugTe58+YAaDBmUhKBGRAhf7xP7Pf8Lhh8PcubDXXlkKTESkgMV+dcd+/ULFrqQuIlKtaFd3\nnDkTxoyBBQuijkREpLAUbcV+441w/fVhdyQREalWlBX7G2/AjBkwYkTUkYiIFJ6iq9jd4YYb4Oab\ntda6iEhdii6xT5gAy5fDz38edSQiIoWpqBK7O5SXwx/+AE2aRB2NiEhhKqrEnkyGav2ss6KORESk\ncBVVYq+q1psW5SVfEZH8KJrEnkyGbe/OPjvqSEREClvRJHZV6yIi6SmKxJ5MwscfwznnRB2JiEjh\nK4rEfvPNYaapqnURkW0r+MT+6quwZAmcd17UkYiIFIeCT+yq1kVEGqagE/ukSfDBB6rWRUQaIpM9\nT083s1lmttnMDq3xeDsz22Bm01LHvY1t46674JprYPvtG/sOIiKlJ5MOjplAT+CBOp57390PyeC9\nmTcPXnsNhg3L5F1EREpPoxO7u8+FsE1TLgwYAL17w8475+TtRURiK1eXJNub2TRgLfB7d3+9IS9e\ntgyeeQbmz89NcCIicVZvYjezsUCrOp66wd1HbeVlS4G27r461ff+nJkd6O6fpxvUoEFw7rnQokW6\nrxARkSr1JnZ3P6ahb+juFUBF6vZUM/sA6ABMrX1ueXn517cTiQSJRIJ16+DBB2Hy5Ia2LCISP8lk\nkmQy2aDXmLtn1KiZTQB+6+5TUvdbAKvdfbOZ7QO8Chzk7mtqvc7ravvuu2HKFBg+PKOwRERiycxw\n93ovbjY6sZtZT2AQ0ILQlz7N3Y83s9OAm4FNwBbgJnd/sY7X/0tir6iAffaBUaPgkIzG1IiIxFNO\nE3um6krsDz8cKvUxYyIJSUSk4KWT2Atmov6WLdC/PwweHHUkIiLFrWCWFBgzBnbYAbp1izoSEZHi\nVjCJfdAguPJKyNF8JxGRklEQfezz58NRR8HixdCsWSThiIgUhXT62AuiYh88GC69VEldRCQbIq/Y\n166F9u1h5kxo0yaSUEREikZRVOyPPALHHqukLiKSLZFW7JWVzn77wWOPwQ9+EEkYIiJFpeAr9tGj\nYY89oKwsyihEROIl0sQ+aBBccYWGOIqIZFOkXTGtWjmLFsGOO0YSgohI0Sn4rphevZTURUSyLdKK\nfdkyp1Vd23iIiEidim51RxERqV/Bd8WIiEj2KbGLiMSMEruISMwosYuIxIwSu4hIzDQ6sZtZfzOb\nY2YzzOxvZrZbjeeuN7MFZjbXzI7NTqgiIpKOTCr2McCB7t4ZmA9cD2BmHYEzgY5AD+BeMyu5XwbJ\nZDLqEHJKn6+4xfnzxfmzpavRCdfdx7r7ltTdScC3U7dPAYa7+yZ3XwS8D3TJKMoiFPd/XPp8xS3O\nny/Ony1d2aqkLwZeSt3eG/ioxnMfAVptXUQkT5rW96SZjQXqmvR/g7uPSp1zI1Dh7k/W81aaYioi\nkicZLSlgZhcClwLd3X1j6rG+AO5+R+r+y0A/d59U67VK9iIijZCztWLMrAcwAOjq7itrPN4ReJLQ\nr94GGAfsq4VhRETyo96umG0YDOwAjLWwU8Zb7t7b3Web2QhgNlAJ9FZSFxHJn8hWdxQRkdyIdHy5\nmd2amuA03czGm1nbKOPJtvomccWBmZ1uZrPMbLOZHRp1PNlgZj1SE+sWmNnvoo4n28xsqJmtMLOZ\nUceSbWbW1swmpP5NvmdmV0QdUzaZWTMzm5TKl7PN7H+2em6UFbuZNXf3z1O3+wCd3f0XkQWUZWZ2\nDDDe3beY2R0A7t434rCyxsz2B7YADwDXuPvUiEPKiJk1AeYBRwMfA+8AZ7v7nEgDyyIz+xGwHhjm\n7p2ijiebzKwV0Mrdp5vZLsAU4NSY/ffbyd2/NLOmwOvAb9399drnRVqxVyX1lF2AlVs7txjVM4kr\nFtx9rrvPjzqOLOoCvO/ui9x9E/AUYcJdbLj7a8DqqOPIBXdf7u7TU7fXA3MI82piw92/TN3cAWgC\nfFbXeZFP9Tez28xsMXABcEfU8eRQzUlcUpjaAEtq3NfkuiJlZu2AQwgFVWyY2XZmNh1YAUxw99l1\nnZfJqJh0A6l3kpO73wjcmBr//kfgolzHlE1ZnMRVkNL5fDGikQQxkOqGeQa4MlW5x0aqB+Dg1PW6\nv5tZwt2Ttc/LeWJ392PSPPVJirCi3dbnS03iOgHonpeAsqwB//3i4GOg5gX8tnxzeQwpcGa2PTAS\neNzdn4s6nlxx97Vm9iJwGJCs/XzUo2I61Lh7CjAtqlhyITWJ61rglKqZuTFW70y4IjEZ6GBm7cxs\nB8Iqpc9HHJOkycKEmiHAbHcfGHU82WZmLcxs99TtbwHHsJWcGfWomGeA/YDNwAfAL939k8gCyjIz\nW0C4yFF1geMtd+8dYUhZZWY9gUFAC2AtMM3dj482qsyY2fHAQMKFqSHuvtUhZcXIzIYDXYE9gU+A\nm9z94Wijyg4zOwp4FXiX6m6169395eiiyh4z6wQ8SijItwMec/f+dZ6rCUoiIvES+agYERHJLiV2\nEZGYUWIXEYkZJXYRkZhRYhcRiRkldhGRmFFiFxGJGSV2EZGY+X8LKTugrrJxuAAAAABJRU5ErkJg\ngg==\n", "output_type": "display_data", "metadata": {}}]}, {"cell_type": "markdown", "source": ["Let's build a Vandermonde matrix at the centers:"], "metadata": {}}, {"input": ["\n", "nodes = centers\n", "\n", "V = np.array([\n", " radial_basis_function(nodes, i)\n", " for i in range(len(centers))\n", " ]).T"], "cell_type": "code", "prompt_number": 179, "collapsed": false, "language": "python", "metadata": {}, "outputs": []}, {"cell_type": "markdown", "source": ["Find the coefficients:"], "metadata": {}}, {"input": ["coeffs = la.solve(V, f(nodes))"], "cell_type": "code", "prompt_number": 180, "collapsed": false, "language": "python", "metadata": {}, "outputs": []}, {"cell_type": "markdown", "source": ["Find the interpolant:"], "metadata": {}}, {"input": ["\n", "interpolant = 0\n", "for i in range(len(centers)):\n", " interpolant += coeffs[i] * radial_basis_function(plot_x, i)\n", "\n", "pt.figure(figsize=(8,8))\n", "pt.ylim([-5,5])\n", "pt.plot(plot_x, interpolant, label=\"Interpolant\")\n", "pt.plot(plot_x, f(plot_x), label=\"$f$\")\n", "pt.plot(centers, f(centers), \"o\")\n", "pt.legend(loc=\"best\")"], "cell_type": "code", "prompt_number": 185, "collapsed": false, "language": "python", "metadata": {}, "outputs": [{"text": [""], "prompt_number": 185, "output_type": "pyout", "metadata": {}}, {"text": [""], "png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHaCAYAAADYNZpeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8jdfjB/DPSSKJJCLLDBFib1XUjr1Lqy2t1my/1SFG\nqZJaX6UUtarVZfzqi7ZKW9SuGKWE2hKJWJGQLRJk3vP746JCEje549zxeb9eeZEnN/d+XK1PznnO\ncx4hpQQRERGZlp3qAERERLaIBUxERKQAC5iIiEgBFjAREZECLGAiIiIFHIz9AkIILrMmIiKbI6UU\nhX3dJCNgKSU/dPiYNm2a8gyW8MH3ie8T3yvreJ+6/dANWyO2Kv9zGeNDF5yCJiIiJS6lXEI1z2qq\nYyjDAiYiIpPL1eTiWuo1+Hv4q46iDAvYjAQGBqqOYBH4PumG75Pu+F7pxpDvU0xaDHxcfODs4Gyw\n57Q0Qte56mK/gBDS2K9BRESWZe/lvZgWMg37h+1XHcUohBCQT1mEZfRV0AURotBcZIb4gxQRGYqt\nn/8FFBYwwH/QLQl/YCIiQ7qUcgkBngGqYyjFc8BERGRyl25xBMwCJiIik+MUNAuYiIgUiEqOYgGr\nDkBERLYlNSMV93LuoaxrWdVRlGIBW7DAwEB8//33qmMQERXJ5VuXUc2zms0v7mQB58Pf3x979ux5\n6uNUF6AQwiD/AYeEhKBy5coGSERE9HRcAa3FAs6HrsWmb/lpNBq9vp+IyBJxAZYWC7gAUkqsWrUK\nbdq0wYQJE+Dl5YVq1aph+/btAIDg4GAcOHAA77//PkqVKoWgoCAAQHh4OLp06QJvb2/Url0bP//8\n88PnHDp0KN555x307NkTbm5u2Lt3L4YOHYqRI0eia9eucHd3R2BgIK5du/bwew4dOoRmzZrBw8MD\nzZs3x+HDh/PNGxUVhY4dO8LHxwdlypTB66+/jtTU1Idf9/f3x4IFC9CoUSN4eHhg4MCByMzMxJ07\nd9CjRw/ExsaiVKlScHd3x82bN43xlhIRAeACrAdYwAV4MLo9evQoateujaSkJHz44YcYMWIEAGDW\nrFlo27Ytli1bhrS0NCxZsgR37txBly5d8PrrryMhIQHr16/Hu+++i7CwsIfPu27dOkyZMgXp6elo\n06YNAGDt2rWYOnUqEhMT0bhxYwwaNAgAkJycjF69emHMmDFITk7GuHHj0KtXL6SkpOSbOTg4GDdu\n3EBYWBiio6Mxffr0PH+en3/+GTt27MDly5dx+vRprFq1Cq6urti+fTsqVqyItLQ03L59G+XLlzfG\nW0pEBIDXAD9g1gUshP4f+qpSpQpGjBgBIQQGDx6MGzduID4+/uHXH93Na8uWLahatSqGDBkCOzs7\nNG7cGC+++GKeUXC/fv3QsmVLAICTkxMAoHfv3mjTpg0cHR0xa9YsHD58GNevX8fWrVtRq1YtDBo0\nCHZ2dhg4cCBq166N33///YmcAQEB6NSpE0qUKAEfHx+MHTsW+/bty/OYoKAglC9fHp6enujTpw9O\nnjz5xJ+BiMjYOAWtpXQryqcxh154dDTo4uICAEhPT0fZstrl84+eB7569SqOHDkCT0/Ph8dycnIw\nePDgh4+tVKlSnud//Jirqyu8vLwQGxuLGzduwM/PL8/jq1SpgtjY2CdyxsXFYfTo0Th48CDS0tKg\n0Wjg5eVV4J+lZMmS+T4PEZEx5WhyEJ0abdO3IXzArEfA5u7xRVh+fn5o3749UlJSHn6kpaVh2bJl\nBT6HlBLR0dEPP09PT0dycjJ8fX1RsWJFXL16Nc/jr169Cl9f3yeeZ/LkybC3t8fZs2eRmpqKH374\nQedFXrZ+KQARmc611Gso51bOpm9D+AALuAC6TMuWK1cOUVFRDz/v3bs3IiIisGbNGmRnZyM7Oxuh\noaEIDw8v9Dn/+OMP/PXXX8jKysKUKVPQsmVL+Pr6okePHoiIiMC6deuQk5ODH3/8EeHh4ejdu/cT\nz5Geng5XV1e4u7sjJiYG8+bN0/nPWq5cOSQlJeH27ds6fw8RUXFcTL6I6l7VVccwCyzgAjy4FOnx\n0eGjn48ePRobNmyAl5cXxowZAzc3N+zcuRPr16+Hr68vKlSogEmTJiErKyvPcz7+fK+99hpmzJgB\nb29vnDhxAmvWrAEAeHt7Y8uWLViwYAF8fHwwf/58bNmy5YmpZQCYNm0a/vnnH5QuXRp9+vRB//79\nCx3ZPpqldu3aePXVV1GtWjV4eXlxFTQRGU1kUiSqe7KAAUAYewGOEELm9xr3b1Zs1Ne2BMOGDUOl\nSpUwc+ZM1VEKxb8vIjKEsdvHwtfdF+NbjVcdxaju/5tZ6Pk9joAVY6kRkS25mMIp6AdYwIoZajtJ\nIiJLEJkUyQK+j1PQpBP+fRGRvnI1uXCd7YqUiSkoWaKk6jhGxSloIiIyG9G3o1HGtYzVl6+uWMBE\nRGQSkUmRqOFVQ3UMs8ECJiIik+A1wHmxgImIyCRYwHmxgImIyCQikzkF/SgWMBERmQRHwHmxgImI\nyOhyNbm4fOsyArwCVEcxG2Z9O0JzsnDhQiQmJsLPzw9vv/226jhERBbl+u3r8C7pDZcSLqqjmA0W\nsA5SU1Px008/YfHixXB1dVUdh4jI4nD6+UksYB0cOXIEjRs3RvPmzVVHISKySBeTL3IB1mPMuoDF\nDP33SJbT9Ns+8ciRI1i8eDEqVqyITZs24YUXXtA7ExGRrYlM5h7QjzPrAta3PA2hRYsWKFmyJMaM\nGYN69eqpjkNEZJEuJl9Ey0otVccwK1wFrYOwsDDUrVtXdQwiIosVmRyJGt6cgn4UC/gp4uLi4OPj\nw1sGEhEVk0ZqcCnlEgI8eQnSo1jAT3HkyBG0bt1adQwiIosVczsGns6ecHXkVSSPMsg5YCGEPYBj\nAK5LKfsY4jlVO378OL799lt4eXlhwIABquMQEVksTj/nz1CLsEYDOA+glIGeTzl7e3tUqlQJLi4u\naNSokeo4REQW62LyRVT35Arox+ldwEKISgB6ApgFYJzeicxE48aN0bhxY9UxiIgsHjfhyJ8hzgEv\nBDABgMYAz0VERFaGU9D502sELIToDSBeSnlCCBFY0OOmT5/+8PeBgYEIDCzwoUREZGVsYQQcEhKC\nkJCQIn2PkLL4m10IIWYDeANADgBnAO4AfpFSDn7kMTK/1xBCQJ/XJtPi3xcRFYdGauA22w3xE+Lh\n5uimOo7J3P83s9DrV/WagpZSTpZSVpZSVgUwEMCfj5YvERHZtti0WJR2Lm1T5asrQ18HzCESERE9\nZAvTz8VlsL2gpZT7AOwryvdwdykiIuvGuyAVTNnNGHg+kYjI+kUkRbCAC8CtKImIyGguJF1ALZ9a\nqmOYJRYwEREZzYXEC6jlzQLODwuYiIiMIjs3G1duXeEirAKwgImIyCgu37oMX3dfODk4qY5illjA\nRERkFJx+LhwLmIiIjOJCEgu4MCxgIiIyiguJXAFdGBYwEREZBUfAhWMBExGRUfAa4MKxgImIyOBu\nZdzCvex7qOBWQXUUs8UCJiIig7uQeAE1vWtyz/9CsICJiMjgOP38dCxgIiIyOF4D/HQsYCIiMjiu\ngH46FjARERkcp6CfjgVMREQGlavJRVRyFO8D/BQsYCIiMqhrqdfg4+IDV0dX1VHMGguYiIgMitPP\numEBExGRQXEFtG5YwEREZFBcAa0bFjARERkUp6B1wwImIiKD4hS0bljARERkMOlZ6UjJSEHl0pVV\nRzF7LGAiIjKYiKQIVPeqDjvBenkavkNERGQwnH7WHQuYiIgMhiugdccCJiIig+EKaN2xgImIyGA4\nBa07FjARERmElBKRyZGo6V1TdRSLwAImIiKDiL4djdJOpVHaubTqKBaBBUxERAZxPuE86papqzqG\nxWABExGRQbCAi4YFTEREBnE+4Tzq+NRRHcNisICJiMggOAIuGhYwERHpTUrJAi4iFjAREentRvoN\nODk4wdvFW3UUi8ECJiIivXH0W3QsYCIi0tv5hPOo68MCLgoWMBER6Y0j4KJjARMRkd5YwEXHAiYi\nIr1IKXEu4RwLuIhYwEREpJeEuwkAgLKuZRUnsSwsYCIi0suD6WchhOooFoUFTEREeuEWlMXDAiYi\nIr1wAVbxsICJiEgvLODiYQETEZFeWMDFwwImIqJiS7qbhLvZd+Fbyld1FIvDAiYiomILSwzjCuhi\nYgETEVGxcfq5+FjARERUbCzg4mMBExFRsbGAi48FTERExcYCLj4WMBERFUtqRipuZdyCX2k/1VEs\nEguYiIiKJTwxHLV8asFOsEqKg+8aEREVC29BqB8WMBERFcvZ+LOoV6ae6hgWiwVMRETFcib+DBqW\na6g6hsViARMRUbGciTuDBmUbqI5hsVjARERUZPF34pGRk4FK7pVUR7FYLGAiIiqyM3Ha6WfuAV18\nLGAiIiqyM/GcftYXC5iIiIrsTNwZNCjHAtYHC5iIiIrsdPxpjoD1xAImIqIiydXk4nzCedQvW191\nFIvGAiYioiK5lHIJZVzKoLRzadVRLBoLmIiIiuRMPM//GgILmIiIiuR03Gk0LMsdsPTFAiYioiLh\nCNgwWMBERFQk3ILSMFjARESks7vZdxF9Oxo1vWuqjmLxWMBERKSzc/HnUMu7FkrYl1AdxeKxgImI\nSGc8/2s4LGAiItIZz/8aDguYiIh0xi0oDYcFTEREOntwG0LSHwuYiIh0EpcehxxNDiqWqqg6ilVg\nARMRkU4eLMASQqiOYhVYwEREpJPTcTz/a0h6F7AQorIQYq8Q4pwQ4qwQIsgQwYiIyLycief5X0My\nxAg4G8BYKWU9AM8BeE8IUccAz0tERGaElyAZlt4FLKW8KaU8ef/36QDCAPAMPRGRFcnR5CAsMQz1\ny9ZXHcVqOBjyyYQQ/gCaADhiyOclytHkIP5OPG5l3IKDnQOcHZxRwa0Ct8MjMpHwxHD4lvJFKadS\nqqNYDYMVsBDCDcAGAKPvj4Qfmj59+sPfBwYGIjAw0FAvS1bqyq0r+DX8VxyKPoRjsccQfTsaPi4+\n8HD2QI4mB/ey7yHhbgLKuZZDkwpN0Lpya3Su1hlNyjfhCk0iIzhx4wSaVGiiOobZCgkJQUhISJG+\nR0gp9X5hIUQJAFsAbJNSLnrsa9IQr0HWZ//Wrdi5ZAkcMjOR4+SEwPdG4lrlFHwZ+iWupl7F8zWf\nR4eqHdCsYjNU96oOezv7PN+fnZuN6NvROBZ7DH9d+wtbI7dCIzUYWH8g3m76Nqp4VFH0JyOyPuN2\njENZ17L4qM1HqqNYBCEEpJSFjgb0LmChHW6sBpAkpRybz9dZwPSE/Vu3Ysfo0ZgVFfXw2Ktedgh7\noQnaNfwENR06o4S9A0qVAry9tR9eXtpf3d2B/Aa5UkqcjjuN1adWY/Wp1WhfpT2mtJuCtH9i8xR9\n16AgtOvVy4R/WiLL12F1B3zU+iN0q95NdRSLYKoCbgNgP4DTAB482SQp5fb7X2cB0xM+7tYNn+zc\n+cTxnj7d4NtvO+ztAY0GuH0bSE4GkpL+/fXePcDT899CfvxXb2/AxSMdR7K+x69/TkfXbdlYGX/n\n4WsEBwSg2+LFLGEiHUkp4fWZFy68fwFlXcuqjmMRdClgvc8BSykPght6UBGkpQGXz9/N92vN62Vg\n+reFf39WFpCSkreUH/312jUgKckNycmjUfPYZqxM25Pn+2dFRWHK0qUsYCIdXbl1BS4lXFi+BmbQ\nVdBET/PHH8CbEyNQpcQ/+X4919n5qc/h6AiUK6f9eJrpgTnAviePXzyVgbg43Z6DyNaduHkCTcpz\nAZahceRKJiElMGMGMHTmTtx7tS0Cg4YhOCAgz2MmBwSgy6hRBn3dHCenfI8nlYxFvfoaLFyoneom\nooKduMECNgaOgMnoMjOB4cOBv++uAV4Yj98G/Ix2Vdphf42tmLJ0KewzMpDr7Izuo0YZfFq4a1AQ\ngqOi8iz2mlDVDym97dGwQT+sX/N/CAnxwA8/aBd3EdGTTtw8geFNhquOYXUMchlSoS/ARVg2LScH\nGDgQiHT/Bon1/oudb+xAvbL1TJph/9at2PVI0XcZNQrPde+C8TvH44/IbWgavgVnQmph+3bAz8+k\n0Ygsgu/nvjg47CCqelZVHcVimGQVtA4hWMA2SqMBRowAjuZ8i/RnZmHPkN2o7lVddaw8Vp5YiY/2\nfIQXc3/E7u8CsX8/UKGC6lRE5iP+TjxqfVELyR8mc5ObItClgHkOmIxm6lTgUOoGJDeaZpblCwDD\nmgzDuv7rsNF+AJoOWY9OnYCEBNWpiMzHiRsn0Lh8Y5avEbCAySh++w34dtdeJD/3Lra9/odZlu8D\nHat2xO43duOAywfwe/Fr9O8PZGerTkVkHrgC2nhYwGRwERHA8A8uIaffq/jx5fVoXL6x6khP1aBc\nA+wbug8RZeYiscYCjB+vOhGReWABGw8LmAwqKwvo/1oanIY+j/92moKOVTuqjqSz6l7VsW/oPmQ0\n+AprLy3EmjWqExGpx5swGA8XYZFBfTxF4vvbL6N3Zy980/trizxvdC31Glp/2wEZP3TG4NJXUMqO\ne0iTbUrLTEP5BeWR+lEqHOx41WpRmGQrSqIHQkOBxX8th/+Ll/BFj/9ZZPkCgF9pPyys/DF+iX4L\nC07nPjwefP9aYpYw2YpTcadQr0w9lq+RcAqaDCIrCxgQdAai41T88up6ODnkvwOVpTi5Yj3WpeTm\nOTYrKgq7li5VlIjI9LgDlnGxgMkg5szPQEL7gVjaez5qetdUHUdvDpmZ+R63z8gwcRIidU7c5Plf\nY2IBk96uXgU+PTwDbWrVweBGg1XHMYiC9pDW5WYRRNaCK6CNiwVMehs8KRT2z67AqpeXWex538d1\nDQp64mYRr3o5oPlb1vEDBtHTZOVm4ULiBTQo10B1FKvFM+ukl83bMvF3uWH4ps9ClHOznnv7PVho\n9eBmEUkZztgt3BBz6zv00LzCRSlk9c7Fn0NVz6pwKeGiOorV4mVIVGy5uYDvgE9RqdVfCB272WpG\nvwUZ8VYu9pTrjX5tamFR90Wq4xAZ1YoTK/Dn5T+x5kVeEF8cvAyJjGrBd9eQVGsBDg0NtfryBYA5\ns+1Ru8labPJsjiblm2BI4yGqIxEZTWhMKJpVbKY6hlXjOWAqlvR0YNrhsRhWNwjVvGzjFmVlygDT\nPvREhX2/Yvyu8Tgac1R1JCKjCY0NRTNfFrAxsYCpWN5fuBMOvqewZMCHqqOY1LvvAumX6mG4z7fo\n/1N/3Ey/qToSkcFl5GTgfMJ5i9jH3ZKxgKnIklNysSbhA8ztNB/ODrZ1WY6DA7B0KfDjjH4YXH8E\nBmwYgBxNjupYRAZ16uYp1PKpxQVYRsYCpiIbtmQVvF088U6HvqqjKNGhA9C8OWD/1xSUsCuBmftm\nqo5EZFChsTz/awpchEVFEpNwB5vTp+Ln/ptsYuFVQebPB5o0scf2g2vQ949n0N6/vUXd+YmoMEdj\njqKtX1vVMaweR8BUJK9/sRBVRFv0f6656ihK+fkBY8YAcz4uj9X9VuONTW8g/k686lhEBhEaG4rm\nvrb9/7gpsIDpqfZv3YqPu3XDxJatkf5/MzC2cifVkczC+PHAiROAuNwFQxoNweBNg6GRGtWxiPSS\nmpGK6NRo1CtbT3UUq8cCpkLt37oVO0aPxic7d2Lu34cQeiUHcYvnYv/WraqjKVeyJPD550BQEDCl\nzX+RlpWG+Yfmq45FpJfjN46jUflG3O3NBFjAVKidS5Zg1v374D7A2/L9q29foFIl4OuvHLCu/zos\nOLwAoTGhqmMRFVtoTCiaV+T0symwgKlQvC1f4YQAFi8GZs0CnDL8sKT7Eryx6Q3cy76nOhpRsXAD\nDtNhAVOheFu+p6tTBxg8GJg8GRhQfwCaVGiCSXsmqY5FVCxHY45yAZaJsICpUF2DgjDU0zPPsckB\nAegyapSiROZp2jRg2zbg6FFgWc9l2HB+A/Ze3qs6FlGRxKXHIT0rHQGeAU9/MOmNZ9mpUHXbBOKX\nTrlwjG6Nis4OyHV2RvdRox7ero+03N2B2bOBUaOAw4e98G2fbzH0t6E4PfI0SjuXVh2PSCehsaF4\ntuKzNn2NvynxdoRUqJfnL8HB6yG4sWij6ihmT6MBWrUC3n4bGDYMGLllJDJzM7Gy70rV0Yh0MnXv\nVGikBp90/ER1FIuny+0IOQVNBbqXmY1NcQvw364fqY5iEezstPtET54MpKYC87vOx/6r+/Fr+K+q\noxHphFtQmhYLmAr04eoNcMuqhrd6ckGGrpo1A3r1AmbMANwc3bC632q8u/VdpNxLUR2NqFBSSu09\ngLkC2mRYwJQvjUbi+3ML8X6zMaqjWJzZs4EffgDOnQPa+LXBi3VexAc7P1Adi6hQl29dhrODMyqW\nqqg6is1gAVO+Fm/8GzmOSZg2sLfqKBanbFntCHj4cCAnB/i006fYc3kPdkXtUh2NqEAc/ZoeC5jy\nNWfvIrxSJQglHOxVR7FII0cCpUpp75pUyqkUlvdajv9s+Q/Ss9JVRyPKF8//mh4LmJ6weV80Et13\nY/HQYaqjWCw7O2DFCmDBAuDsWaBHjR5o69cWH//5sepoRPniBhymxwKmJ0z4cTnalH4d3m7uqqNY\nND8/YO5cYNAgICMDWNhtIX489yMORx9WHY0oj+zcbJy4eQLPVnxWdRSbwgKmPM6cz0JEqe/x+Wsj\nVUexCsOGAbVra29d6O3ijcXdF2PE7yOQmZP/HttEKpyOO40qpavAw9lDdRSbwgKmPMZ89RuquNRB\nU786qqNYBSGAb74B/vgD2LQJeLnuy6juVR3zDs1THY3ooUPRh9C6cmvVMWwOC5geiokB9t9bjuBu\nHP0aUunSwLp12h2yIiMFlvZYioV/L8TF5IuqoxEBAP6K/gutKrdSHcPmsIDpoamLL6CE71kMbvaC\n6ihWp0UL7S0L+/YFPEQVTGw9Ee/98R64TSuZg0PRh1jACrCACQCQkgL8L+JrDG00HI72jqrjWKW3\n3gI6dtQuygpqNhaxabH4+fzPqmORjYtOjUZGTgaqe1VXHcXmsIAJALD0yyyg4RqMCxyhOopVW7QI\nSEsDZs4ogeW9lmPsjrFIzUhVHYts2IPRL++AZHosYMK9e8DnW7agfrk6/CnYyEqUADZsANauBa4f\nbo2e1Xtiyt4pqmORDeP0szosYMLq1YDTcysxqvVw1VFsQpky2hXR778PvFp2Dn469xOOxR5THYts\n1KHrXAGtCgvYxuXkAJ9+EYu7PgfxUt2XVMexGY0bA8uXA0Ne9sbEpnPw7tZ3oZEa1bHIxtzJuoPz\nCefRtGJT1VFsEgvYxq1dC9g/8wMGNOgPV0dX1XFsSv/+wHvvAf/3wWBA2mPVyVWqI5GNCY0NRcNy\nDeHs4Kw6ik1iAduw3Fzgk1kSufVXYVhj7vuswsSJQJPGdij55xcI/jMYtzJuqY5ENuRQ9CG0qsTz\nv6qwgG3Yzz8DzlWPw8E5i4swFBFCOxVtF9cUZZL7YHrIdNWRyIYcij6E1n48/6sKC9hGaTTAJ58A\nAf3+h0ENBvESBIUcHYFffgHSf5+FFaFrcTb+rOpIZAM0UoND0YfQslJL1VFsloPqAKTGr78Czi45\nOJy+DnMa7FMdx+Z5eQE7NpbBMyOnYVD0a+h9tQJKZGYix8kJXYOC0K5XL9URycpcSLwAD2cPVChV\nQXUUm8UCtkFSAjNnAv0++BNb7/mhlk8t1ZEIQI0awPR2lXH08wuYdfvMw+PBUVEAwBImg+L0s3qc\ngrZBW7ZoSzjKdQ0GNRikOg49IuXQMvx4OyvPsVlRUdi1dKmiRGSt/or+iwuwFGMB25gHo98Jk+9i\nc8RmDKw/UHUkeoRDZv73CbbPyDBxErJ23AFLPU5B25gdO4A7dwCHulvRLLMZyrmVUx2JHpHj5JTv\n8VxnXqdJhpN4NxE30m+gftn6qqPYNI6AbYhGAwQHA9OnAxvCfsKAegNUR6LHdA0KQnBAQJ5jE/wC\n0GXUKEWJyBodjj6MFr4tYG9nrzqKTeMI2IZs2KC97rT78+l4c+FOLO+1XHUkesyDhVZTli6FuHcX\nu2OP47bPG5jbgwuwyHD2X92PNn5tVMewecLYNwQXQkjedFy97GygXj3gyy+BpPI/YtWpVdg2aJvq\nWPQUv4ZtxmsrJmJWpVMYO7qE6jhkJZp92wyfd/0cbau0VR3FagkhIKUsdIMFTkHbiBUrAD8/oHNn\n4MdzP+KVuq+ojkQ66Fu7N5pUr4gpv32DiAjVacgapGakIjwxHM19m6uOYvNYwDYgNVV73nfuXCAt\nMw27L+1Gv9r9VMciHQgh8FW/z2EX+F8MGZkCDW+YRHo6eO0gmvs2h5ND/gv+yHRYwDbgk0+Anj2B\npk2BLRFb0MavDTxLeqqORTpqWK4hBjTui8uVP8Hq1arTkKULuRKCwCqBqmMQWMBWLyICWLkSmDVL\n+/nG8I3oX6e/2lBUZJ90nInM2qsx/tNIJCaqTkOWLORqCAL9A1XHILCArZqUwLhx2lvelS8P3Mu+\nh51RO/F8redVR6MiKudWDhPbjofnyx9iwgTVachS8fyveWEBW7ENG4BLl4CgIO3nuy7tQpPyTVDG\ntYzaYFQsY54bgxyfk/j9zF6EhqpOQ5bowLUDaOHbgud/zQQL2EqlpACjRwPffQc82FxpU/gmvFD7\nBbXBqNicHZzxWde5cH1xHEaPyQWv7qOiCrnC6WdzwgK2UhMmAC+8ALS6v9VrjiYHmy9s5upnC/dy\n3ZdRuZwrYsqsxk8/qU5DloYFbF5YwFZo+3Zg507g00//Pbb/6n74e/ijikcVdcFIb0IILOq+EHda\nfIzxk9PAezSQrm5l3MKFpAtoVrGZ6ih0HwvYysTFAcOHAz/8ALi7/3v81/BfOf1sJZr5NkPP2p1R\nsuscLOduoqSjg9cO8vyvmWEBWxGNBhgyBBgxAmjf/t/jUkr8fuF3rn62IrM7zUZ85eWY9cVVpKer\nTkOWgNPP5ocFbEXmzNHuejV1at7jZ+PPAgBvPWZFKrlXwphWo+DS9yMsXqw6DVkCFrD5YQFbiU2b\ngK++An5LjUKvAAAgAElEQVT5BSjx2J79WyK2oE/NPhCi0H3BycJMaDUBmeUOYN76w0hJUZ2GzBnP\n/5onFrAVOHkS+M9/tCVcseKTX98csRl9avUxfTAyKldHV3zWbTYc+4zFZ/O4STQV7MBVXv9rjljA\nFi4sTLvP81dfAc8+++TX4+/E43zCebSv0v7JL5LFe73h66jgm4ule9cjLk51GjJXIVdC0MG/g+oY\n9BgWsAW7cEF7e8G5c4GXXsr/Mdsit6Fztc78yddK2Qk7fNH7c9h1mYSZn95THYfMVMjVELT35w/h\n5oYFbKFCQ4GOHYGZM4E33ij4cZsjNqNPTU4/W7O2VdqiXfVnsSJsEa5dU52GzE3S3SRcTL7I/Z/N\nEAvYAv3887/TzsOHF/y4rNws7L60Gz1q9DBdOFJiUe+5QMsFmDyb89CU157Le9DWry0c7R1VR6HH\nsIAtyJ07wHvvAePHa3e6ev4pl/Ueij6EGt41UNa1rGkCkjLVvapjSOPB2JA4DVevqk5D5mRX1C50\nqdZFdQzKBwvYQuzeDTRpAty+rV313KTJ079nW+Q29KjO0a+tmN1tCuzrb8T4z86qjkJmQkqJnZd2\noksAC9gc6V3AQojuQohwIUSkEGKiIULRv0JDge7dgXfeAebN024x6emp2/duj9rOArYhniU9Edwu\nGL/dncBRMAEAIpMjkaPJQR2fOqqjUD70KmAhhD2ALwB0B1AXwKtCCP5N6+nWLeD//g947jng5ZeB\n3r2Bc+eAvn11f46Y2zG4fvs6F17YmPHt34GbXxTeWbBddRQyAw+mn7kJj3nSdwTcHMBFKeUVKWU2\ngPUAilATBAAZGcDff2tHuD16AH5+wIYNwKRJQFQU8P77gGMR109sv7gdXap1gb2dvXFCk1lytHfE\n4l7zsBPjEXU5R3UcUmzXpV3oGtBVdQwqgIOe3+8LIPqRz68DaKHnc1o0KYGsLODuXe2iqQe/3rkD\npKUBN24AMTH/fkRGAlevArVqAW3bam+k8NNPQKlS+uXYdnEbLz+yUa83ex4fb12IYUtWYP/C/6iO\nQ4rkaHIQciUE3/T5RnUUKoC+BSx1eVCFCtpievhNMu+v5nDMkM9dogTg6gq4uGh/ffB7Nzfte+Hr\nCzRsqB3tBgQANWsWfYRbmOzcbOy5vAfLei4z3JOSxRBCYMXAz9FlZS+cjRyI+jXcn/5NZHX+vv43\n/D38jX4VhJRAfDxw7Zr2IzoaSEnR3hjmwcfdu9qBSVYWkJ2d9/caTd5/Pwv6fWHHLJW+BRwDoPIj\nn1eGdhScx2uvTQcACAG0ahWIVq0C8eCUxKOnJlQfM9Tz2CleWx4aG4oqpaugnFs5tUFImU51nkF9\n524Y/M1c/DNvluo4pIAxroKQEjh/HggJAY4e1V6RERGhHVz4+Wk/KlcGvLyAKlWA0qW1Hy4ugJOT\ndnDi6Jj3V/v7Z8ke/Te1oN8Xdky1Q4dCcPhwyMPPP//86d8jpB4/PgghHABcANAJQCyAowBelVKG\nPfIYqc9rUNHN3DcTqZmpmN91vuoopNCZKzFotLwhDrzxD1rXq6I6DpnYM18/g0XdF6FdlXZ6PY9G\nA+zfrz01tmkTULIk0KGDdpFokyZAnTraWT7KSwgBKWWhPxroNQKWUuYIId4HsAOAPYDvHy1fUmP3\n5d2Y1GaS6hikWAN/X7Qq8T6G/jAZkXP+pzoOmdDN9Ju4fOsyWlZqWeznSE8HVqwAvvgCcHYGXn1V\nW8Q1ahgwqI3Te7JUSrlNSllLSlldSvmpIUJR8d3JuoPjscfR1q+t6ihkBta8MwFRuSH47dhR1VHI\nhHZc3IFOVTuhhH2Jpz/4MVlZ2tKtUQM4cEBbwqdOaa/KYPkaFnfCsjIHrh3AMxWegasj54QI8K/o\nhu5OMzFy0zjwVJDt2HaxeOd/Dx8GGjUCtmwBtm3T7jvfpo15nGO1RixgK7Pn0h50rtZZdQwyIyuD\nhiAhNR3LD/yiOgqZQI4mB7su7UL36t11/p7sbGDCBODFF7V3WNu+HWjc2IghCQAL2OrsubwHnap2\nUh2DzEi5svZ4ufQCfLRrIjJzMlXHISM7GnMUvqV84evuq9Pjb97U3to0PBw4c6bge4uT4bGArUji\n3UREpURx+0l6wtKxnXDvWl3M3PmF6ihkZFsjtqJnjZ46PfbUKaBZM6BzZ+C33wAfHyOHozxYwFZk\n7+W9aOvXtlgLL8i6+fgAQyrOw4Ijc5B4N1F1HDKizRGbddoF79AhoGtXYMECYNo09fsX2CK+5VaE\n089UmNnjakOeGYAJW/6rOgoZydVbV3Ej/Qaeq/RcoY/780/tzV1WrwZeecVE4egJLGArsvvSbnSq\nxgKm/JUpA/yn5jSsP7sOFxIvqI5DRrAlYgt61uhZ6E1YQkOBAQO0K5y7675Oi4yABWwlrt66irSs\nNNQvW191FDJj0z8sA7u/P8Q7mz5UHYWM4GnTz+HhQJ8+wPffA4GBpstF+WMBW4k9l/egY9WOsBP8\nK6WCeXkBwZ1H4cjlM9h7ea/qOGRAaZlpOBR9qMDbDyYmAj17AnPmAM8/b+JwlC/+a20ldl/azfO/\npJNxQc5w+Wsu3t44DrmaXNVxyEB2X9qN5yo9B3enJ+9+lZ2tPdf78svA0KGmz0b5YwFbASkl/rz8\nJzfgIJ04OwPzhr2Em9Eu+L9TP6iOQwbye8TvBU4/f/CB9u999mwTh6JCsYCtwLmEc3B1dIW/h7/q\nKGQhBg8WqHh2Acb/8THuZN1RHYf0lKPJweYLm9G3dt8nvvbLL8DWrcDatf/e+o/MAwvYCuy+tBud\nq3L0S7qzswO+n/EcMi60w+x9vG2lpTtw9QD8PfzhV9ovz/ErV4B33gHWrwc8PNRko4KxgK1AyJUQ\ndKjaQXUMsjCtWwOd5KdY+NcSxNyOUR2H9LAxbCNerPNinmM5OcBrrwETJ2p3uyLzwwK2cBqpwcFr\nB3n7QSqWZbOrQB7/D0b//rHqKFRMUkr8euFXvFD7hTzH584F3NyAsWMVBaOnYgFbuPDEcLg7ueu8\n8TrRoypXBj5qPQlbwrbjn9gTquNQMRyLPQbXEq6oU6bOw2NnzgCLFmmv9+UWk+aLfzUW7sDVA2hb\nhaNfKr7JH7jD68w0DP7fB7xnsAV6fPo5O1t7qdGcOdofsMh8sYAt3IFrBzj9THopUQJYO/5NXIiJ\nw48nt6iOQ0UgpcTG8I15pp/nz9duOzp8uMJgpBNh7J94hRCSP1Ubj99CP+wevBs1vWuqjkIWrkfQ\nNvwT9yZGpNSDY1YWcpyc0DUoCO169VIdjQpwOu40nl/3PC6PvgwhBC5dApo3B44dA/z9VaezbUII\nSClFYY9xMFUYMryrt64iMzcTNbxqqI5CViCoTS5Wv5mC2Wm7Hh4LjooCAJawmfrp3E94pd4r9/+x\nB959F5gwgeVrKTgFbcEeTD8LUegPWUQ6+ev7pVifdi/PsVlRUdi1dKmiRFQYKeXDAga0dzeKiQHG\njVMcjHTGEbAFO3CV53/JcBwyM/M9bp+RYeIkpItTcaeQK3PRtEJT3Lmj3W5y3TrtOX2yDBwBW7AD\n17gCmgwnx8kp3+O5zs4mTkK6+OncT3ilrnb6+dNPgXbtgDZtVKeiomABW6jEu4mISYtBo3KNVEch\nK9E1KAjBAQF5jk2oEoAuo0YpSkQFeXT6+dIlYPly7cYbZFk4BW2hDl47iJaVWsLejrurk2E8WGg1\nZelS2N27h90x/+C27yuY24MLsMzNsdhjEEKgcfnG6N9fu9tVpUqqU1FR8TIkCzVuxzh4l/RGcLtg\n1VHISu2JCkGvb4ci2OM8pnzkojoOPWLM9jHwdPZEJ4dpGDQIuHBBe7tBMh+6XIbEKWgLxfO/ZGyd\nAgLRte5zmHNwLg4fVp2GHsjR5GD92fV4rcEgTJgAzJzJ8rVULGALlJ6VjrCEMDT3ba46Clm5L/vN\nh/1zy/DSW5eQkqI6DQHa24/6e/jjzL7quHsXGDRIdSIqLhawBTocfRhNKjSBswN/7CXjquReCR+1\nHwfnvuPw5psAzyap978z/8PAuoMwaZJ24ZU9l4FYLBawBeL+z2RKH7T8ACh7FqfSd+Crr1SnsW3p\nWenYfGEzsk4OQKVKQLduqhORPljAFogFTKbk5OCExT0WQdMtCFNnZOHkSdWJbNev4b+iRcVWWPRJ\nWXz2GcBN8CwbC9jCZOVmITQmFK0qt1IdhWxI75q9Uad8dXSZshgDBgDp6aoT2aaVJ1fC4/IwtG8P\nNG2qOg3piwVsYU7dPIUArwCUdi6tOgrZmEXdFmHX3blo0i4W77+vOo3tuXLrCk7dOI1dy57HrFmq\n05AhsIAtzN/X/0YL3xaqY5ANquFdAyOfHYnsTmNx9Cjwww+qE9mWVSdXoXLqqxg0wAnVqqlOQ4bA\nArYwf8f8jecqPac6Btmo4LbBOBl/DKMWb8e4cdoNIMj4NFKD74+vQtQvwzBpkuo0ZCgsYAvz93UW\nMKlTskRJLOu5DPPD3sOU/97DgAEAb5ZkfHsv70VmqgeGdW+CihVVpyFDYQFbkIQ7CUi6m4TaPrVV\nRyEb1r16dzSt0BRxNWehRg3gww9VJ7J+yw5/j/QDwzFxouokZEgsYAtyJOYImvk2g53gXxuptbDb\nQnx9fDnGzw3Dxo1ASIjqRNYr/k48/ojchjcavcHRr5Xhv+QW5O/rf+M5X04/k3q+7r6Y2n4qJu5/\nB8uXS4wYAdy5ozqVdVp6YBXk+RcwfaKn6ihkYCxgC8Lzv2RO3mv2HtKz0pFU6f/Qpg24OMjA9m/d\niuBu3XD47anoevQ8Iv/ZqjoSGRjvB2whcjW5CI0NRYtKvASJzIO9nT2W916OXmt74cDsHujQvCz6\n9wfat1edzPLt37oVO0aPxqyoqPtHjiB49GgA/963mSwfR8AWIjwxHGVcysDHxUd1FKKHnq34LAY3\nHIxpf4/G8uXgVLSB7Fyy5JHy1ZoVFYVdS5cqSkTGwAK2EJx+JnM1o8MMhMaEAjU3o2VLIDhYdSLL\n55CZme9xe17zZVVYwBbiSMwRFjCZJZcSLvi2z7d49493MfOzVPz8M3D4sOpUli3HySnf47nOvAWp\nNWEBWwiOgMmcdajaAT2q98Cc4xMxbx7w/vtAbq7qVJar3sA3McAj7z/PkwMC0GXUKEWJyBi4CMsC\npGWmISolCg3LNVQdhahAn3X5DPW/rI8BL4TAzS0Q334LjBypOpVl+iY8DrFdWmLKbTfYZ2Qg19kZ\n3UeN4gIsK8MCtgChsaFoXL4xHO0dVUchKpCHswe+7PUl/rPlLfxv0Sn06e6Cl14CfLhusEhiYjXY\nn7UYv0xaiX5N2qiOQ0bEKWgLcOT6EW7AQRbh+VrP49mKz2Jt3GS8+iowebLqRJbn7YW/oYybJ/o2\nbq06ChkZC9gC8A5IZEmW9VyGDec3oOOIP7F5MxAaqjqR5YiNldie/ilm9fgIQgjVccjIWMBmTkqp\nvQcwN+AgC+FV0gvfPf8dRu0ZhqmzU/Hee4BGozqVZRj1+V64l7mNYS37qY5CJsACNnNXbl2BvbBH\nZffKqqMQ6ax79e7oWb0n/vYYjRIlgBUrVCcyf7GxwO/JczC140TecMVG8G/ZzD24/pfTUWRp5nWd\nh7+iD+KFyZvw8cdAaqrqROZt3PxjcK4chnfbDlIdhUyEBWzmjlw/gha+nH4my+Pm6IbV/VZjQfi7\n6NgnHrNmqU5kvmJjgY1JM/BR2w95tYMNYQGbuWM3juHZis+qjkFULK39WmNoo6FIaDUC36+QuHhR\ndSLzNHb+UTj5n8QHgW+pjkImxAI2Y7maXJy8eRLPVHhGdRSiYpvRYQZSc+LQcswSfPih6jTmJzYW\n2JgyDR+3D4azA7eatCUsYDMWkRSBsq5l4VmSN+Imy+Vo74j1L63HUadZ+PvacYSEqE5kXoLmHYaL\nXxjGth+uOgqZGAvYjB2/cRxNKzRVHYNIb9U8q+GLnl9A8+IABE24zX2i77tyReK39EmY1uFjnvu1\nQSxgM3Ys9hgLmKzGK/VeQd+GnXGj6UisXClVxzELb87bAg/fBAS1G6o6CinAAjZjx28cR9OKLGCy\nHou6LYRHzbP44H8rcPu26jRqnT2fgxCHD/FF38/gYMdt+W2RkNK4P4kKIaSxX8Ma5Wpy4THXA1fH\nXIVXSS/VcYgMJiwhDK1HtUC7o/XR2M8ROU5O6BoUZHN3+mn69te45fsjLk7Zw+v8rZAQAlLKQv9i\n+WOXmXqwAIvlS9Ym4eglvHbABV/EHgYua48FR0UBgM2U8P6jt3Cy9HTsG7SF5WvDOAVtprgAi6zV\nziVL8EVsXJ5js6KisGvpUkWJTG/wqilo6dUXbQL4/7gt4wjYTB2PZQGTdXLIzMz3uH1GhomTqLFi\n2wlEu/+EQyPPq45CinEEbKa4AIusVY6TU77Hc52tfxOKXI0GY3e9h8G+s1DRw1t1HFKMBWyGNFLD\nHbDIanUNCkJwQECeYwM97VDlte6KEpnOuyu/RHa2wNcjuekGcQraLEUkRcDHxYcLsMgqPVhoNWXp\nUthnZOB2rjP+jGuKfXHz0flWP/h7+KsNaCRhcZfw3cUZ+KbzQTiW4NiHWMBm6VjsMU4/k1Vr16tX\nnhXPblOB3yLLo9uabvhr+F/wcfFRmM7wNFKDPt+8idqJH2JE31qq45CZ4I9hZogLsMjWBAcDOX+N\nQh3ZH73W9sKdrDuqIxnU7N3LcDX2LjZOGKc6CpkRFrAZ4iVIZGucnIBVq4BDn8xCVde6eGXDK8jO\nzVYdyyBO3jyJmfv/izc9/4daNe1VxyEzwgI2Mw8WYHEKmmxNs2bAW28K3Fn/DSCB/2z5Dyx9F707\nWXfQe9VAeB1ZjIVTA57+DWRTWMBmhguwyJZNnQrciCmBtnE/ITwxHEHbgiy2hKWUGLbxbSSfbomf\nprwGG7jKioqIBWxmjsfy+l+yXU5OwIYNwMK5rphabTuOxh7FmO1jLLKE5x2ajz2nwvGG55do21Z1\nGjJHLGAzw/O/ZOv8/YGVK4G33iiNFR124PD1w3j/j/ehkRrV0XS2LXIbPv1zEcrs+RUL55VUHYfM\nFAvYzLCAiYCePYExY4CXentgfc/dOJtwFoM3DbaIhVnHYo9h0IYhkD/9hN9+qAQXF9WJyFyxgM2I\nRmpw4sYJ7oBFBGD8eKB/f+ClPu5Y33s7UjNT0Xd9X6RlpqmOVqALiRfQc00fOGz7FsuDW6MWL/ml\nQrCAzUhUchQ8S3rC24V7xBIBwMyZQMeOQMe2JfFp442o7F4ZbVe2xfXb11VHe0JUchQ6re4Ku72z\n8VG/vhg4UHUiMnfcCcuMnIo7hUblGqmOQWQ2hADmzwdq1wY6ti+BpUuXo3qDBWjxXQv8+NKP0JxJ\nxc4lS+CQmYkcJyd0DQpSck/h8wnn0WlVV9gd/BgjWwzDOO63QTpgAZuRUzdPoXH5xqpjEJmdN98E\n6tcH3n5bwM1tPILG1cOr03ujz04HfBmT9PBxwVFRAP7db/rOHSA8HAgL0/568yaQnAxkZwMlSwIe\nHkBAAFCrFtCyJVCunG559m/d+rD4E2UGfqlyAXevLMHHfd7A+PEG/+OTlRLGXt4vhJCWeAmBCs+v\nex5DGg1B/7r9VUchMku5ucCaNcC8eUDpS+3w170DTzzmpSrdkF57O8LCgIQEoEYNoE4d7Sja1xfw\n9AQcHYF797RlfPGitqAPHwbKlwcCA4H27bUfFSo8mWH/1q3YMXo0Zt0vewB41ak82kz/Du99ZPrR\nN5knIQSklKKwx3AEbEZOxZ1Co/KcgiYqiL09MGSI9mNiCzvg6JOPcbfLwND3taXr76/9Hl3k5gKn\nTgH79gHr1wPvvQd4eQGtWmnLOyAAcHAA/pixBN89Ur4AsC7zJqaELAVYwFQELGAzkXIvBcn3klHN\ns5rqKEQWoYSHU77H3f2z0bt30Z/P3h545hntx9ixgEYDnDsHHDkCXLgAHDsGSAlobibk//0ZGUV/\nUbJpeq2CFkLME0KECSFOCSE2CiFKGyqYrTkddxoNyjaAneDCdCJddA0KQnBA3v2V3/P1wS/+5/D6\nxtcRmRSp1/Pb2QENGmjPP8+bByxeeR1ubwxGose5fB+fy70mqYj0/dd+J4B6UspGACIATNI/km3i\nCmiiomnXqxe6LV6MKd26YXr79pjSrRsGfL0K57+KRm2f2mj5fUv0WdcHOy7uQK4mt9ivcy7+HN76\n/S00/Koh/Er74eM5a54o/skBAegyapS+fySyMQZbhCWEeAFAfynl648d5yIsHQz/bTia+zbHyGdH\nqo5CZBXuZt/F2jNr8dWxr3D99nX0q9UPXQO6olXlVqhQKp/VVfdl5WbhdNxp7IrahY3hG3H99nW8\n8+w7GPnsSJR1LQtAuxBr19KlsM/IQK6zM7qMGqXk8icyX7oswjJkAW8GsE5Kufax4yxgHTT9pim+\n6PEFWlZuqToKkdWJSo7CpvBNCLkSgkPRh+Bo74iqnlVRsVRFONo7QkAg8W4ibqbfxMXki6juVR3t\nqrRD/zr90bZKWzjYcbkMFY1BClgIsQtA+Xy+NFlKufn+Y4IBPCOlfOL6GSGEnDZt2sPPAwMDERgY\n+PT0NiRHkwP3T90RPyEebo5uquMQWTWN1OBG2g1cuXUFN9JvIDs3GxqpgbeLN8q5lkNN75pwdXRV\nHZMsTEhICEJCQh5+PmPGDOOPgIUQQwG8BaCTlPKJZYAcAT/dufhz6PdjP0SO0m/RCBERmQejXwcs\nhOgOYAKA9vmVL+mGC7CIiGyPvquglwJwA7BLCHFCCPGlATLZnFM3WcBERLZGrxGwlLKGoYLYslNx\np/Bus3dVxyAiIhPirg9mgFPQRES2hwWsWPydeGTkZMCvtJ/qKEREZEIsYMUenP8VotDFckREZGVY\nwIpx+pmIyDaxgBXjLQiJiGwTC1gxXoJERGSbWMAKZeZkIjI5EvXK1lMdhYiITIwFrFBYYhiqeVaD\nswPvI0pEZGtYwApx+pmIyHaxgBU6HXcaDcs1VB2DiIgUYAErdDbhLBqUbaA6BhERKcACVuhc/Dku\nwCIislEsYEVS7qUgNTOVW1ASEdkoFrAi5xLOoW6ZurAT/CsgIrJF/NdfkXPx51C/TH3VMYiISBEW\nsCJn48+iflkWMBGRrWIBK3IugQuwiIhsGQtYEY6AiYhsGwtYgfg78cjWZKOCWwXVUYiISBEWsALn\n4s+hftn6EEKojkJERIqwgBU4l8AV0EREto4FrMDZ+LNcgEVEZONYwApwARYREbGATUxKqb0EqQxH\nwEREtowFbGKxabFwtHdEGdcyqqMQEZFCLGAT4+iXiIgAFrDJ8fwvEREBLGCTe3ANMBER2TYWsImd\nTTjLKWgiImIBm5JGanA+4TyvASYiIhawKV1LvYbSTqXh4eyhOgoRESnGAjYhLsAiIqIHWMAmdC6e\nlyAREZEWC9iEziZwBExERFosYBPiJUhERPQAC9hEcjW5CE8MR50ydVRHISIiM8ACNpErt66grGtZ\nuDm6qY5CRERmgAVsIuGJ4ajtU1t1DCIiMhMsYBNhARMR0aNYwCYSlhiGOj48/0tERFosYBPhCJiI\niB7FAjYBKaV2BMwV0EREdB8L2AQS7yZCSokyLmVURyEiIjPBAjaBsMQw1PapDSGE6ihERGQmWMAm\nEJ4YzgVYRESUBwvYBLgAi4iIHscCNgEuwCIiosexgE2AI2AiInocC9jI7mbfxc30m/D38FcdhYiI\nzAgL2MgikiJQ3as6HOwcVEchIiIzwgI2srCEME4/ExHRE1jARsZLkIiIKD8sYCMLT+ICLCIiehIL\n2Mg4BU1ERPlhARtRriYXkcmRqOVdS3UUIiIyMyxgI7py6wrKupaFq6Or6ihERGRmWMBGxAVYRERU\nEBawET24CxIREdHjWMBGxBEwEREVhAVsRNwDmoiICsICNhIpJaegiYioQCxgI0m8mwgpJcq6llUd\nhYiIzBAL2EgejH6FEKqjEBGRGWIBGwkXYBERUWFYwEZyIfECanrXVB2DiIjMFAvYSCKSI1DLh1tQ\nEhFR/ljARhKRFIEaXjVUxyAiIjPFAjaCHE0Ort66igCvANVRiIjITLGAjeDKrSuoUKoCnB2cVUch\nIiIzxQI2Ak4/ExHR07CAjSAyKZIroImIqFAsYCOISIpgARMRUaFYwEYQmRzJKWgiIioUC9gIOAIm\nIqKnYQEbWEZOBm6m30QVjyqqoxARkRljARtYVHIU/D384WDnoDoKERGZMRawgXH6mYiIdMECNjBe\nA0xERLpgARtYZDKvASYioqdjARsYp6CJiEgXehewEOIDIYRGCOFliECWLiIpAjW8OQVNRESF06uA\nhRCVAXQBcNUwcSzb7czbSMtKQ8VSFVVHISIiM6fvCPhzAB8aIog1iEzS7oBlJzizT0REhSt2Uwgh\n+gK4LqU8bcA8Fo3Tz0REpKtCd4sQQuwCUD6fLwUDmASg66MPL+h5pk+f/vD3gYGBCAwMLEpGixGZ\nHImaXlyARURka0JCQhASElKk7xFSyiK/kBCiPoA9AO7eP1QJQAyA5lLK+MceK4vzGpbo9Y2vo3O1\nzhjaeKjqKEREpJAQAlLKAgemQDGnoKWUZ6WU5aSUVaWUVQFcB/DM4+Vra3gJEhER6cpQq4VsY4hb\nCCklN+EgIiKdGeSOAVLKaoZ4HkuWeDcRAOBd0ltxEiIisgS8XsZAHox+hSh0yp+IiAgAC9hgeP6X\niIiKggVsILwLEhERFQUL2EC4AIuIiIqCBWwgHAETEVFRsIANQEqJi8kXuQ0lERHpjAVsADfSb8DN\n0Q3uTu6qoxARkYVgARtAVHIUAjwDVMcgIiILwgI2gIvJFxHgxQImIiLdsYANICqFI2AiIioaFrAB\nsICJiKioWMAGEJUcxSloIiIqEhawAUSlRKG6V3XVMYiIyIKwgPV0K+MWsnKzUMaljOooRERkQVjA\neieP7V8AAAUASURBVHpwCRLvgkREREXBAtYTL0EiIqLiYAHriSugiYioOFjAeuIuWEREVBwsYD1x\nBTQRERUHC1hPUSm8BpiIiIqOBayHjJwMJNxJQGX3yqqjEBGRhWEB6+FyymX4lfaDvZ296ihERGRh\nWMB64PQzEREVFwtYDxeTL3IFNBERFQsLWA9RyVwBTURExcMC1gM34SAiouJiAeuB54CJiKi4WMDF\nlKvJxdVbV1HVo6rqKEREZIFYwMV0/fZ1+Lj4oGSJkqqjEBGRBWIBFxO3oCQiIn2wgIuJlyAREZE+\nWMDFFJXMBVhERFR8LOBi4iVIRESkDxZwMfESJCIi0gcLuBiklNopaI6AiYiomFjAxZB4NxEl7EvA\ns6Sn6ihERGShWMDFYKzzvyEhIQZ/TmvE90k3fJ90x/dKN3yfDIsFXAzGWgHN/7h1w/dJN3yfdMf3\nSjd8nwyLBVwMUSlRqOZRTXUMIiKyYCzgYrh86zKqebKAiYio+ISU0rgvIIRxX4CIiMgMSSlFYV83\negETERHRkzgFTUREpAALmIiISAEWMBERkQImKWAhxEwhxCkhxEkhxB4hRGVTvK6lEULME0KE3X+v\nNgohSqvOZK6EEC8LIc4JIXKFEM+ozmNuhBDdhRDhQohIIcRE1XnMlRBihRAiTghxRnUWcyaEqCyE\n2Hv//7mzQogg1ZnMkRDCWQhx5H7XnRdCfFro402xCEsIUUpKmXb/96MANJJSvmn0F7YwQoguAPZI\nKTVCiDkAIKX8SHEssySEqA1AA+BrAB9IKf9RHMlsCCHsAVwA0BlADIBQAK9KKcOUBjNDQoi2ANIB\n/J+UsoHqPOZKCFEeQHkp5UkhhBuA4wD68b+pJwkhXKSUd4UQDgAOAhgvpTyY32NNMgJ+UL73uQFI\nNMXrWhop5S4ppeb+p0cAVFKZx5xJKcOllBGqc5ip5gAuSimvSCmzAawH0FdxJrMkpTwAIEV1DnMn\npbwppTx5//fpAMIAVFSbyjxJKe/e/60jAHsAyQU91mTngIUQs4QQ1wAMATDHVK9rwYYD+EN1CLJI\nvgCiH/n8+v1jRHoTQvgDaALtIIEeI4SwE0KcBBAHYK+U8nxBj3Uw4IvuAlA+ny9NllJullIGAwgW\nQnwEYCGAYYZ6bUvytPfp/mOCAWRJKdeaNJyZ0eW9onzx4n4yivvTzxsAjL4/EqbH3J/FbHx/Dc8O\nIUSglDIkv8carICllF10fOha2PDI7mnvkxBiKICeADqZJJAZK8J/U5RXDIBHFzpWhnYUTFRsQogS\nAH4BsEZK+avqPOZOSpkqhNgK4FkAIfk9xlSroGs88mlfACdM8bqWRgjRHcAEAH2llBmq81iQQrd7\ns0HHANQQQvgLIRwBDADwu+JMZMGEEALA9wDOSykXqc5jroQQPkIIj/u/L/n/7dwhTkNREIXhf9gC\n3QOKHTSpQlSyDNgBhg3UNPWEEEQNC0BjkCRYwgbqugByEC+kqnV0HuH/9E3uuJObmTvABQfy7lhT\n0E/AGfAFfAJXSTa/fvEfU1UfDI37n6b9a5LrxpJGq6ougRUwAbbAW5J5b1XjUVVzYMkwBHKX5OB3\niP+qqtbADDgFNsBtkvveqsanqqbAC/DOrsVxk+S5r6rxqapz4IHhcXsCPCZZ7D3vLmhJko7PTViS\nJDUwgCVJamAAS5LUwACWJKmBASxJUgMDWJKkBgawJEkNvgEAZ1SORNsO7wAAAABJRU5ErkJggg==\n", "output_type": "display_data", "metadata": {}}]}, {"cell_type": "markdown", "source": ["* Play around with the radius of the RBFs\n", "* Play with node placement"], "metadata": {}}], "metadata": {}}]}