{ "cells": [ { "cell_type": "code", "execution_count": 109, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "\n", "def f(x):\n", " return np.sin(np.pi * x)\n", "\n", "x = np.linspace(0,1,100)\n", "x = x**2\n", "\n", "h = np.diff(x) # we define the vector h as the difference of successive x[i] using the built in diff function\n", "\n", "f = f(x)" ] }, { "cell_type": "code", "execution_count": 110, "metadata": {}, "outputs": [], "source": [ "A = np.zeros((len(f),len(f))) #we create the matrix A by initializing to 0" ] }, { "cell_type": "code", "execution_count": 111, "metadata": {}, "outputs": [], "source": [ "for i in range(1,len(f)-1):\n", " A[i,i-1] = 2/(h[i-1]*(h[i-1] + h[i])) #we populate the diagonal (minus the extreme ends), the one to the left and right with a,b,c respectively\n", " A[i,i] = -2/(h[i-1] * h[i])\n", " A[i,i+1] = 2/(h[i]*(h[i] + h[i-1]))\n", "\n", "A[0,0] = 1 # we denote the far corners of our matrix to be 1\n", "A[len(f)-1,len(f)-1] = 1\n" ] }, { "cell_type": "code", "execution_count": 112, "metadata": {}, "outputs": [], "source": [ "u = np.linalg.inv(A).dot(f) # we solve our system" ] }, { "cell_type": "code", "execution_count": 113, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 113, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAv70lEQVR4nO3deVxU5fs+8IsdcmMQYRAQyEAhdwXFJRUFxA3XUlNREHC30oL0U7YnpV8zd0c03Jc0QVHZUlEDm1QQFYRRVEYcEFBxAdme3x8WvwyQbWaeGeZ+v173HzOc4VxP6lydWc7RAsBACCFEY2nzDkAIIYQvKgJCCNFwVASEEKLhqAgIIUTDUREQQoiG0+UdoCFyc3Nx584d3jEIIUSt2NjYwMzMrMr9alkEd+7cgbOzM+8YhBCiVsRicbX300tDhBCi4agICCFEw1EREEKIhqMiIIQQDUdFQAghGk4uReDp6Ym0tDRkZGQgKCio2m3WrFmDjIwMJCcno3v37vV6LCGEEMVijRltbW0mkUiYnZ0d09PTY0lJSczR0fGVbby8vNjx48cZANa7d2+WmJhY58dWN2KxuFGZaWhoaDRxanrubPQRgYuLCyQSCTIzM1FaWop9+/bB29v7lW28vb2xY8cOAMCFCxdgbGwMoVBYp8fK01suPeHmN01hv58QQhRFz9AA3kEfQGAhlPvvbnQRWFpaIisrq/K2VCqFpaVlnbapy2P/4e/vD7FYDLFYDFNT0wZl7djfFV4LAmHazqpBjyeEEF76jPfGO1Pfg6CtChaBlpZWlfsYY3Xapi6P/YdIJIKzszOcnZ2Rl5fXoKxnwvagvLQMQ/x9GvR4QgjhQVdfH4N9p0IivoRbF5Pk/vsbXQRSqRTW1taVt62srJCdnV2nberyWHl6kl+AhINH0HPkMJhYtVXYfgghRJ5cxo5EK7M2iN28XSG/v9FFIBaLYW9vD1tbW+jp6WHSpEmIiIh4ZZuIiAhMnz4dANC7d288fvwYMpmsTo+Vt1O/7AarqKD3CgghakFHVxduftOQefkKMi78pZB9NLoIysvLMX/+fERFRSE1NRUHDhzA9evXERgYiMDAQADA8ePHcevWLUgkEohEIsydO/e1j1WkwtwHuHD4KJy9R8DYvOpZ+AghRJX08h4OgYUQMZu2KXQ/3D/SVN9p7MdHjYXmLORSPBu3bAn3tdDQ0NDUNNq6OmzpiV/Zoj2hcvl9Cvv4qDp6JMuBODwSLmNHomWbhn0CiRBCFK3nCE+0trJEtIKPBjSyCADg9607oK2jg8Ezp/KOQgghVWjr6GDILB9Ir99Aavx5xe5Lob9dhRXcu4+Lx07CdeIYtGhtwjsOIYS8otuwIWhj2w4xCvqk0L9pbBEAQJwoDDp6uhjoM4V3FEIIqaSlrY0h/jNwP+Mmrp2KV/j+NLoI8u5KcflEDPq+Nw7NBMa84xBCCACgi/tgCNvbIWbz9hq/ZCtPGl0EABC75RfoGRrgnWmTeEchhBBoaWlhaMAMyG5m4krMKaXsU+OLIDfzDq5E/47+UybgjVYtecchhGi4twcPQFuHtxAn+gWsokIp+9T4IgCAmC2/wLBZMwyY+h7vKIQQDece6IsHt+8i6WSc0vZJRQBAlnETV2JOYcCUiTBs0Zx3HEKIhnJ8px+snDogVhSGivJype2XiuBvMZu3w6hlC/SfMpF3FEKIhvKY7Yu8LCkuHY9S6n6pCP6WfSMDV0/FY+C0STBo9gbvOIQQDdOxfx+06+yEONEOVJQp72gAoCJ4Rcym7XijVUv0mzSBdxRCiIbxmO338ouuR08ofd9UBP8ivZ6G6/HnMchnMvSNjHjHIYRoCAdXZ9h07YS40B0oLytT+v6pCP4jZvN2NBMYo+9743hHIYRoCI/Zfi9Phnkkksv+qQj+4+6Va7hxPhGDZkyBnqEB7ziEkCauvXMP2PXoit9Dd6K8tJRLBiqCakRv3IYWrU3gOnEs7yiEkCbOY7YvHv99wSxeqAiqcTs5BemJYgye+T50DeiogBCiGG/27Ia3XHri1LZdKCsp4ZaDiqAGMZu2oWUbU/QZP5p3FEJIE+UeOBOFeflIPBTONQcVQQ1uXUyCRHwJbr7ToKuvzzsOIaSJsenaCQ6uLji9fTdKi19wzUJF8Boxm7ahlXkb9B43incUQkgT4zHbD08LHiLh4G+8o1ARvI7kz4u4dTEJbn7ToKOnxzsOIaSJaNfZCR3798HpX3ajpKiYdxwqgtrEbN4GY6E5XMaO5B2FENJEeMx5eTRwft9h3lEAUBHUKj1BjNtJKS+PCnR1ecchhKg5605OcBzQF2d27EVJURHvOACoCOokemMoTNpawHnMCN5RCCFqzmO2L549eozzew/xjlKJiqAObvxxAXeSr2LILB86KiCENJiVU0c4DeyHM2F78eL5c95xKlER1FH0plCYWFqg12gv3lEIIWrKY7Yvnj8uxLm9B3lHeQUVQR2lnUvE3ZTrGOI/A9q6OrzjEELUjKWjA94ePABnduzFi2eqczQAUBHUS/TGULS2aoteo4bzjkIIUTMes33xvLAQ5/ao1tEAQEVQL6ln/8Ddq9cxNMCHjgoIIXXWtoM9OrkNRPzO/Sh++ox3nCqoCOopeuM2tLayRM+Rw3hHIYSoCffZvigqfIKzuw/wjlKtRhWBQCBAdHQ00tPTER0dDWNj42q38/T0RFpaGjIyMhAUFFR5/w8//IDU1FQkJyfj8OHDaNWqVWPiKEVq/HlkXUuFe+BMOioghNTKwqE9ugwdhPhd+1H85CnvONVqVBEEBwcjLi4ODg4OiIuLQ3BwcNUdaGtj/fr18PLygpOTEyZPngxHR0cAQExMDDp16oSuXbsiPT0dn376aWPiKA0dFRBC6spjtp9KHw0AjSwCb29vhIWFAQDCwsIwZsyYKtu4uLhAIpEgMzMTpaWl2LdvH7y9vQG8LILy8nIAQGJiIqysrBoTR2munzmHrOtpGOo/A9o6dFRACKle2w726OI+GPG79qOo8AnvODVqVBGYm5tDJpMBAGQyGczMzKpsY2lpiaysrMrbUqkUlpaWVbbz9fXFiRMnatyXv78/xGIxxGIxTE1NGxNbLqI3hMK0nRV6jqKjAkJI9f55byB+137eUV6r1q/JxsTEQCgUVrl/2bJlddqBlpZWlfsYY6/cXrp0KcrKyrB79+4af49IJIJIJAIAiMXiOu1bka6fOYesa6kYGjADF4+dREVZOe9IhBAV0raDPboMHYSo9SKVfW/gH7UWgbu7e40/y8nJgVAohEwmg1AoRG5ubpVtpFIprK2tK29bWVkhOzu78vb06dMxcuRIDBkypL7ZuYvaEIpZ61ei58hhEB+J5B2HEKJCPOb44XlhIeJV+L2BfzTqpaGIiAj4+PgAAHx8fBAeXvVya2KxGPb29rC1tYWenh4mTZqEiIgIAC8/TRQUFITRo0ejSEXOwlcf9AkiQkh1LDs6oPOQv783oOJHA/9gDR0TExMWGxvL0tPTWWxsLBMIBAwAs7CwYJGRkZXbeXl5sRs3bjCJRMKWLl1aeX9GRga7e/cuu3z5Mrt8+TLbuHFjnfYrFosbnFne4/hOP7YqJYG5jBnJPQsNDY1qzMw1K9jX56OYYfNm3LP8e17z3Mk/nBwXw2UW7Q1lS0/8yrR1dbhnoaGh4TuWjg5sVUoCGxo4k3uW/05Nz530zWI5iN4QitZWlnAeTecgIkTTec6Z9fIMo2rw3sA/qAjkIPXsH7hz5RqGBsyk6xUQosGsnDpUnmFUFc8pVBMqAjmJ3rgVJpZ0FTNCNJnH30cDqvwt4upQEchJ2rlE3Em+iqEBM6Cjp8c7DiFEyaycOuLtQf1xOmyPyl1voDZUBHJ0cr0IAgsheo8bxTsKIUTJhs33x7NHj1XyegO1oSKQo/SEP5F5KRlD/H2gq6/POw4hRElsunaC44C+OLV9l9odDQBUBHJ3cr0IxuZm6DNhNO8ohBAlGTbPH0/yC3B+7yHeURqEikDOJH9ehER8CW5+06FrYMA7DiFEwd7s2Q0Ori74fdtOlKjhGRIAKgKFiFovQiuzNuj77ljeUQghCuY5zx+FD/KQcOA33lEajIpAAW5dTEJ6ohhuftOgb2TIOw4hREHecumJt5x7IG5rGEqLX/CO02BUBAoStU6EFq1N0G/SeN5RCCEKMmx+AB7JcpD4awTvKI1CRaAgt5NTkHouAYNnToXBG2/wjkMIkbMOfXvDrnsXxG4JQ1lJCe84jUJFoEBR60RoJjDGgKnv8o5CCJGzYfMDUHDvPv787SjvKI1GRaBAWddScfVUPAb5TIFRyxa84xBC5MTxnX5o19kJMZu3o7ysjHecRqMiULCT60QwatkCA30m845CCJGTYfP8kXdXir+OHucdRS6oCBTsfroESVFxGPD+u2gmMOYdhxDSSJ2HDoKVUwdEb9rWZK5VTkWgBFHrRdA3NMTgmVN5RyGENIKWtjaGzQ9Azq3buBQZxTuO3FARKEFu5h1cioxG/8kT0MK0Ne84hJAG6jHcA8L2dji5XgRWUcE7jtxQEShJ9MZQaOvqYGjADN5RCCENoK2rA4+5friXmo6UmFO848gVFYGS5Evv4c8jx9BngjcEFkLecQgh9eQydhRMra1wYu1mMMZ4x5ErKgIlit20HWAM7rN9eUchhNSDrr4+3ANn4nZSClLP/sE7jtxRESjRo5xc/LH/N/Qa7QVTG2vecQghddT3vXEwNjfD8Z838Y6iEFQEShYXGoayklIMm+fPOwohpA70jYzg5jcN6Ql/4qb4Eu84CkFFoGRP8x/i7K796O7lDguHt3jHIYTUYsDUd9GitQmO/7yZdxSFoSLg4HTYHjwvLITXgkDeUQghr2HUsgUGz3gfV0/FI+vqdd5xFIaKgIOiwic4tW033h7UH7ZdO/OOQwipwSCfKTBo3gwn123hHUWhqAg4ObfnAJ7kF8Br0WzeUQgh1WjeWoABU99F8slY3E+/yTuOQlERcFJSVIyYzdvxlnMPdOjbm3ccQsh/uAfMhK6+Pk6uF/GOonBUBBwlHjyCfGk2vBbNhpaWFu84hJC/mVi1RZ+JY3Dh8FHk3ZXyjqNwVAQclZeVIWrDVlg7dURn98G84xBC/uY5dxYqysoRs2kb7yhK0agiEAgEiI6ORnp6OqKjo2FsbFztdp6enkhLS0NGRgaCgoKq/Hzx4sVgjKF1a807IdulyCjcz7gJr/kB0NbR4R2HEI1n4dAePUZ44uzuAyh8kMc7jlI0qgiCg4MRFxcHBwcHxMXFITg4uOoOtLWxfv16eHl5wcnJCZMnT4ajo2Plz62srODu7o47d+40JoraYhUVOLF2M8zsbODsPZx3HEI03vCFc1D89ClObd/FO4rSNKoIvL29ERYWBgAICwvDmDFjqmzj4uICiUSCzMxMlJaWYt++ffD29q78+erVq/HJJ580uZM41ce1U2dxOykFHnNnQdfAgHccQjSWXY+ucBrYD7+H7kRR4RPecZSmUUVgbm4OmUwGAJDJZDAzM6uyjaWlJbKysipvS6VSWFpaAgBGjRqFe/fu4cqVK7Xuy9/fH2KxGGKxGKampo2JrZIif9oAY3Mz9J88gXcUQjTWiEVz8Dj3Ac7tOcg7ilLp1rZBTEwMhMKqp01etmxZnXZQ3adhGGMwMjLCsmXL4OHhUaffIxKJIBK9/BiXWCyu02PUya2LSUg9+weGzJqOxEPhKH7ylHckQjSK4zv9YNejKw5+FYLS4he84yhVrUXg7u5e489ycnIgFAohk8kgFAqRm5tbZRupVApr6/9/pk0rKytkZ2ejffv2sLOzQ3JycuX9ly5dgouLC3JychqyFrUX+dNGfHQwDG6+03B8zUbecQjRGFra2hjxwRw8uH0Xf/52lHccpWvUS0MRERHw8fEBAPj4+CA8PLzKNmKxGPb29rC1tYWenh4mTZqEiIgIXL16Febm5rCzs4OdnR2kUil69OihsSUAvLzQ/aXIKLwz9T20NGvDOw4hGqPHcA9Y2LfHiXVbmswF6eujUUWwYsUKuLu7Iz09He7u7lixYgUAwMLCApGRkQCA8vJyzJ8/H1FRUUhNTcWBAwdw/XrTPXlTY51ctwVa2lrwnOPHOwohGkFHTw/D5gcg63oarkT/zjsON0zdRiwWc8+gyBn9ySL2Y9I5ZmZnwz0LDU1Tn3emT2KrUhKYfe9e3LMoemp67qRvFquguC2/4MXzIoz4YA7vKIQ0aUYtW8A9YCZSzyUg48JfvONwQ0Wggp49eoxT23ahk9tA2HXvwjsOIU3WEL/pMGzRHJGr1/OOwhUVgYqK37UPj3MeYOTi+byjENIkCSyE6P/+RFw8eqLJn2a6NlQEKqq0+AWiNohg27UzOg8ZyDsOIU3OsPkBAICTa5v2RWfqgopAhYnDj0MmuYURH8yFti6dkI4QebHs6IAeIz1xdtd+PMqp+v0nTUNFoMIqystxbPUGtLFtB9cJY3jHIaTJGPHhXBQVPkFc6E7eUVQCFYGKS40/D8mfF+Exxw+GzZvxjkOI2nNwdUGHvr0Rs3k7ncrlb1QEauDoqrVobiLAYN9pvKMQota0tLUx8qN5yJfewx/7D/OOozKoCNSA9PoNXDx2EgOnTYKxedUzvBJC6qbnyGGw7OiA42s2oby0lHcclUFFoCZO/LwZ0AK8Fs7mHYUQtaRvZAivhYG4e/U6kqPieMdRKVQEauLhfRnid+5Hr9FesHLqwDsOIWpnkM8UGJubIeKHNRp9IazqUBGokbitYXiSX4BRSxbyjkKIWmlp1gaDZk5FcvTvyLxc+4WwNA0VgRp58ew5ojZsxVvOPfD24AG84xCiNrwWBEBHVwfHNPxUEjWhIlAzFw5FQHYzE6M+mg8d3VqvK0SIxrN0dECv0cNxdtcBFEizecdRSVQEaqaivBxHV61FG9t26EfXNyakVqOXLMTzR48RK/qFdxSVRUWghtLOJuDG+US4z56JN1q15B2HEJXVye0dvOXSEyfXi1D89BnvOCqLikBNRaxcC8NmzeA5dxbvKISoJB1dXYz8aD5kNzNx4VAE7zgqjYpATckkt5D4azhc3x0L8zdtecchROX0nTQebWyscXTlz6go17zrENcHFYEaO7lehJLnRRj98SLeUQhRKW+0agmP2b64cT4RaecSecdReVQEauzZw0eI3rQNHfv3geOAvrzjEKIyhs0PgEGzNxD+48+8o6gFKgI1d37vr8jNvIPRHy+kj5MSAsDC4S24ThyD8/sOIedmJu84aoGKQM2Vl5Uh/Mc1MLOzQb8p9HFSQsZ++hGKCp8gemMo7yhqg4qgCUg7m4DUs3/AY7YfmrcW8I5DCDfdPIegfa/uOP7zJhQVPuEdR21QETQR4T+sgb6hIYYvnMM7CiFc6BsZYuTi+ZBev4ELh4/yjqNWqAiaiAe37yJ+1370HjcK1p2ceMchROkG+06DwEKIIyv+D6yigncctUJF0ITEbN6Gwrx8jA3+EFpaWrzjEKI0JpYWGDzzfVw6Hk1nF20AKoIm5MWz54hcvQE2XTuh12gv3nEIUZpRixegorwCx/5vHe8oaomKoIm5ePQEbienYMSH82DYojnvOIQonH3vXujiPhhxW8PwOOcB7zhqiYqgiWGM4fC3K9FMYEznISJNno6uLsZ8+hHypfdwJmwv7zhqi4qgCbqXmo6EA7+h/+QJsHBozzsOIQrzzvRJELa3w2/fr0ZZSQnvOGqrUUUgEAgQHR2N9PR0REdHw9jYuNrtPD09kZaWhoyMDAQFBb3ys/nz5yMtLQ1Xr15FSEhIY+KQfzmxdguKCp9g7NLFvKMQohDGQnO4B/ri6u9nkBp/nncctccaOiEhISwoKIgBYEFBQWzFihVVttHW1mYSiYTZ2dkxPT09lpSUxBwdHRkANmjQIBYTE8P09fUZANamTZs67VcsFjc4syZN73Gj2KqUBNZz5DDuWWho5D0+q79n3/95igkshNyzqMu85rmz4b80LS2NCYUv/xCEQiFLS0ursk2fPn3YyZMnK28HBwez4OBgBoDt37+fDRkyRJ6LofnXaGlpsQW7trAvTkcyo5YtuOehoZHXdOzfh61KSWBuftO5Z1Gnqem5s1EvDZmbm0MmkwEAZDIZzMzMqmxjaWmJrKysyttSqRSWlpYAAAcHBwwYMACJiYk4ffo0evXqVeO+/P39IRaLIRaLYWpq2pjYGoMxhsPfrEQz41bwWhDIOw4hcqFrYICxSxcjN/MOzoTt4R2nSaj1dJUxMTEQCoVV7l+2bFmddlDdF5sYYy93rqsLgUCAPn36wNnZGQcOHMCbb75Z7e8RiUQQiUQAALFYXKd9E+BeWjrO7f0V/adMxJ+/HYX0+g3ekQhpFDffqTC1tsKmWQtQXlbGO06TUGsRuLu71/iznJwcCIVCyGQyCIVC5ObmVtlGKpXC2tq68raVlRWys7Mrf3b48GEAL5/cKyoqYGpqiry8vHovhNQsar0I3TyHYPxnn+Dn9/3p6/dEbbW2soSb3zRcPh6NjAt/8Y7TZDTqpaGIiAj4+PgAAHx8fBAeHl5lG7FYDHt7e9ja2kJPTw+TJk1CRMTL64ceOXIEbm5uAAB7e3vo6+tTCShA8dNniPjxZ7Tr5IS+743jHYeQBhu79COUl5YhYuVa3lGanAa/8WBiYsJiY2NZeno6i42NZQKBgAFgFhYWLDIysnI7Ly8vduPGDSaRSNjSpUsr79fT02M7d+5kKSkp7OLFi2zw4MGNesOD5vUTsGk1+zYhlrU0q9uns2hoVGm6eLixVSkJbMDU97hnUddRyKeGVHAxNK+Z1laWbIX4NJu+6lvuWWho6jNGLVuw5aeOsQ/2bWPaOjrc86jrKORTQ0S95EvvIWbLdnT1cKNrHBO1MuLDuWhm3AoHv1iBivJy3nGaHCoCDXN6+27IJLcw7n9LoG9kxDsOIbV6s2c3uE4Yg/id+3EvLZ13nCaJikDDlJeV4eCXITBpa4FhCwJ4xyHktXT19TFxeTDypfcQvXEr7zhNFhWBBrqddAXn9x3CgCkTYf22I+84hNRoyKzpMLOzwaGvf0RJUTHvOE0WFYGGOr5mI57kFeDdLz+Ftq4O7ziEVGH+pi3cZk3HxWMnceOPC7zjNGlUBBqq+OkzHP5uJdp2sMcgn/d5xyHkFVpaWpi4PBgvnj1H+A9reMdp8qgINNjV3+ORHP07POb4wszOhnccQir1mTgGdj26IuLHn/Hs4SPecZo8KgINd/i7lSgpKsZ7Xy2Dljb9dSD8CdoKMfKjeUhP+BN/RRznHUcj0L98Dfc0/yGOhKyGbbfO6D95Au84hODdL5cCAA4s/55zEs1BRUBw6VgUrsefh9fC2TCxass7DtFgfSaOgUMfZxxdtQ4P78t4x9EYVAQEAPDrVyGoKC/He18urfbU4YQomqCtEKMWz0d6wp9IPHiEdxyNQkVAAACPcx4g4oc1eMulJ/pOGs87DtFA9JIQP1QEpNKfR44h9ewfGPHBXLS2tuIdh2iQPhO86SUhjqgIyCsOfLECFWVlmPzN/+hTREQpBBZCjFqyAOmJYnpJiBP6l05eUZj7AL+tWA27Hl3xztT3eMchTZyWlhbe/fJTAMCB5d9xTqO5qAhIFRePnsDV38/Aa2EghG9Vfw1pQuSh3+QJcHB1efmSUDa9JMQLFQGp1sEvQ1D89BmmfL8cOrq1XtqakHozb2+HkR/Nw7XT5+glIc6oCEi1nhY8xMEvvodlRwd4zJ3FOw5pYnT09PD+ii9Q/PQZvSSkAqgISI2unT6HC4ci4OY7FXbdu/COQ5oQr/kBsOzogP2ff4enBQ95x9F4VATktcJ/WIOH92WY8v0XMGzejHcc0gS079UdA2dMwR8HfkNq/HnecQioCEgtXjx/jt3BX6CVeRuM/+wT3nGImjNs0RyTv/sc+XelOLryZ95xyN+oCEit7iRfRfSmbegx3AM9Rw7jHYeosfHLlqBlG1PsDv6CrjimQqgISJ3EicJw8+JljFu2BK2tLHnHIWqo58hh6DHCE9GbtiHrWirvOORfqAhInbCKCuwJ/hKsogLvh3xJHykl9WJmZ4Pxn30MifgS4kRhvOOQ/6AiIHX2SJaD/cu/g02XtzF80RzecYia0DUwwLSV36C0+AV2B38BVlHBOxL5DyoCUi8psadxft8hDJoxBY7v9OMdh6iBMcEfoK3DW9iz9CsU5j7gHYdUg4qA1FvEjz/jXmo6Jn/7GVqZt+Edh6iw7l7ucJ0wBnFbd+DG+UTecUgNqAhIvZWVlGDHx/+Djp4upv3wNbR1dXhHIirI1MYaE5YHIfNSMk6u28I7DnkNKgLSIHl3svDrlyGw69EVIxbN5R2HqBhdfX1M//EblJeWYdcnn6OivJx3JPIaVASkwS6fiMG5PQcxaMYUdB46iHccokLGfPohLB0dsHfZ13iUk8s7DqlFo4pAIBAgOjoa6enpiI6OhrGxcbXbeXp6Ii0tDRkZGQgKCqq8v2vXrkhISMDly5chFovh7OzcmDiEg4gff8ad5KuY9PX/YGpjzTsOUQF9JnhXvi9Ap5BQH6yhExISwoKCghgAFhQUxFasWFFlG21tbSaRSJidnR3T09NjSUlJzNHRkQFgUVFRbNiwYQwA8/LyYqdOnarTfsVicYMz08h/jM3N2JdnjrMlh3cxfSND7nlo+E27Lm+zkEvxzH/jaqalrc09D82rU9NzZ6OOCLy9vREW9vLLIWFhYRgzZkyVbVxcXCCRSJCZmYnS0lLs27cP3t7eAADGGFq2bAkAaNWqFbKzsxsTh3DyKCcXu4KWw/xNW7z39f94xyGctGhtghn/9z0eyXKwK2g5fV9AjTTq66Hm5uaQyV5eVUgmk8HMzKzKNpaWlsjKyqq8LZVK0bt3bwDABx98gKioKKxcuRLa2tro27dvjfvy9/dHQEAAAMDU1LQxsYkCZCSKEbl6A0YtWYB7qen4PXQH70hEiXR0dTF91bcwbNEcomkfoqiwkHckUg+1FkFMTAyEQmGV+5ctW1anHWhpaVW5jzEGAJgzZw4+/PBDHD58GBMnTkRoaCjc3d2r/T0ikQgikQgAIBaL67Rvolynw/bA0qkDvBYGIjs9A2lnE3hHIkoy+uOFeLNnN+z8+DPcT7/JOw5pgAa/3pSWlsaEQiEDwIRCIUtLS6uyTZ8+fdjJkycrbwcHB7Pg4GAGgD169OiVbR8/ftyo17lo+I+eoQH76EAY++Z8NGtj2457HhrFj7P3cLYqJYGNWryAexaa149C3iOIiIiAj48PAMDHxwfh4eFVthGLxbC3t4etrS309PQwadIkREREAACys7MxcOBAAICbmxsyMjIaE4eogNLiF9i+KAjlZWXwW/sjjP5+D4g0TW/26o4Jy4ORnihG5E8beMchjdDgdjExMWGxsbEsPT2dxcbGMoFAwAAwCwsLFhkZWbmdl5cXu3HjBpNIJGzp0qWV9/fr14/99ddfLCkpiSUmJrIePXo0qtVoVGdsu3VhIRfPsDmh65iOri73PDTyH1Mba/b1uSj2SfheZtSyBfc8NLXPa547+YeT42JoVGh6jPRkq1IS2LtfLuWehUa+80arliz42AH25ZnjzMSqLfc8NHWbmp476aTyRGEuHYuCma0N3ANn4sHtOzi1fTfvSEQOdPT0MGPNChgLzbDJbwEKpPSxb3VHRUAUKmq9CG1s22HkR/Px8H4Okk7G8o5EGundLz5F+57dsfPjz3A7OYV3HCIHVAREoRhj2Lv0K7QwNcHkbz/Dk/wC3BRf4h2LNJD7bF/0Gu2FE2s3U6k3IXTSOaJwZSUl2L4wCHl3pZi5JgRC+/a8I5EGcH13LIbN84c4PBKxW37hHYfIERUBUYqiwicQzfkIJc+LELBxNQQWVb+kSFRXt2FDMW7ZElw7dRYHvviedxwiZ1QERGkeyXIgmvsh9I0MEbhlDZq3FvCOROqgQ78+mPLdcmReSsaOjz9DRRldW6CpoSIgSnU//Sa2zl2MlmZtELh5DYxatuAdibyGbbcumLH6e9zPuIltCz5G2YsXvCMRBaAiIEp3OzkFv3wQBDM7G8zasAr6Rka8I5FqWDi0h9/6H/E4JxeiOR+i+Okz3pGIglAREC7SE8TY+fHnsH7bEb5rf4CeoQHvSORfzNvbIXDLzygpKsbmgEV4WvCQdySiQFQEhJurv5/Bvs++QXvnHvBbu5LKQEUI7dtjTug6VJSXY9OsBXh4X8Y7ElEwKgLC1aVjUdi77Cu0d6EyUAUWDu0xZ+taVJSVY4PvPDy4fZd3JKIEVASEu8oycO4Ov3UroW9kyDuSRmrbwR5ztq5DWWkpNvjORd6drNofRJoEKgKiEirLoFd3BGz6CYYtmvOOpFEsOzpg9ta1KCkuxoYZc5F3V8o7ElEiKgKiMi5FRmPnx5/BurMT5m5bT98zUJI3e3XHnNB1ePH8OTb4zkO+9B7vSETJqAiISrkScwqh8xbDtJ015odthqAtfQNZkbq4D0bg5p/w+EEe1vvMoTOJaigqAqJy0hPE2BywEM0ErbBgxxa07WDPO1KT1G/yBExb+Q2yrqVh3fTZeCTL4R2JcEJFQFTSneSrWD9jLhirwLywjeg4wJV3pCbFa+FsjFu6GNdPn8Um/4UoKizkHYlwREVAVJYs4ybWTJmFvDtS+K39Ea7vjuUdSe3p6utj8nefY6i/DxIOHkHYR8votBGEioCotsIHeVg/Yw5SzyZgwmefYPQni6Cto8M7lloyNjfDvLCN6DXq5fUEfv0qBBXldAI5QkVA1EBJURG2LwpC/M79GDht0sszl5rQJ4rq482e3fDB/u0ws7XBtoWf0PUEyCuoCIhaYBUVCP/hJ+z59EvYdOmEDw/8gnadnXjHUgv9Jo3HbNFaFBU+wZopfrh26izvSETFUBEQtXLx2EmsnRaA8tIyzPtlI/pNnsA7ksoyaPYGJn/7OcYtW4Ibf1zAmil+yM28wzsWUUFUBETt3EtLx0+TZiI9UYxxSxfDd+2PaCYw5h1Lpdh27YzFv+5AjxEeiNqwFdsWfEynkSY1oiIgaun540KEzluC375fBQdXZyw5tBMOri68Y3GnraMDz7mzMC9sIwAtrJ8xF9EbQ8EY4x2NqDAqAqLWzu35FT9N9sOzR48RuGUNxn/2CQybN+MdiwtTG2vMC9sIjzl+uHgsCqsmTMPtpCu8YxE1oMs7ACGNJcu4iZ8m+8FrfgDemfYe3h7UH4e/XYWrv5/hHU0pdA0MMNTfB4Nnvo+S4mLs/PgzJJ2M5R2LqBmmbiMWi7lnoFHNsXLqyD46EMZWpSQwn9XfM0FbIfdMipyO/fuwT48fZKtSEtiU75ezFq1NuGeiUd2p6bmTjghIkyK9noafpvhi4PTJ8JjtB8cBrji7az/itu5oUm+Wtra2wogP5qCrhxtyM+9go998SP68yDsWUVNUBKTJqSgrx6ltu3ApMgpeC2bDzW86nMeMRPTGUFw4fBTlpaW8IzaYsdAc7oEz4TxmBMpLy3D85004/csetV4T4U8LLw8N1IpYLIazszPvGERNWDo6YPSShXjLpSce5eTi9PbdSDwUjtJi9TnHTgvT1hgyazpcJ44BACQcPII4URie5BfwDUbUyuueOxv8epNAIGDR0dEsPT2dRUdHM2Nj42q3Cw0NZTk5OSwlJaVBj//v0HsENA0ZB1dnNmfberYqJYF9eeY4G+Lvw5qbCLjnet1Yd3Jik7/9nIVcPMN+uHyWTVgexIyF5txz0ajn1PTc2agjgpCQEBQUFCAkJARBQUEQCAQIDg6ust2AAQPw9OlT7NixA507d6734/+LjghIY9h264KhAT5wHNAX5aVluHoqHom/HkFG4l8q8Xl7XQMDdPVwQ//JE9CusxOKnz3DX+HHEb/rAPKz6BKSpOEUckSQlpbGhMKXn8oQCoUsLS2txm1tbGyqHBHU5/H/HjoioJHHmNnZsFFLFrCv4k+wVSkJbOmJQ2zU4gXMrnsXpqWtrdQs+kZGrKvnEDZt5TfsuwtxbFVKAvskfC/rN2k8M2j2Bvf/VjRNYxRyRPDw4UMIBP//LJAFBQUwMTGpdlsbGxscO3bslSOC+jze398fAQEBAABTU1PY2dk1NDYhr9DR00PnIQPRa7QX7Hv3gq6+Pp7kF+D6mfO4+ddl3E5Kkfv/ievo6sLS0QG23bugfa/u6ODaG3qGBniSX4CU2NNIiorDTfElue6TkJqOCGr91FBMTAyEwqrXjV22bJl8ktWRSCSCSCQC8HIxhMhLeWkpkk7GIulkLAyavYGO/V3Rye0ddBk6CL3HjQIAPMkvwN0r13Bfcgv5WfeQnyVFftY9FOblv/ac/roGBhBYmMPU2gqtrS3R2toSlh0d0K6TE/QMDQAAeVlSJB4Kx5WYU8i8fAWsokIp6ybkH7UWgbu7e40/y8nJgVAohEwmg1AoRG5ubr123tjHEyJvL549R3JUHJKj4qClrQ3zN21h260zbLt1hk2XTug4wBU6uq/+syl98QIlz4vw4nkRyktLoWdkCAMjI+gbGUFH79Vti589Q+6tO/jj4G+4ffkKMi9fwZO8fGUukZAqGvU9goiICPj4+CAkJAQ+Pj4IDw9X6uMJUSRWUQGZ5BZkkltI/PXl301tHR0YW5jD1NoSra2s0Ly1APpGhjB44w3oGxlCV08PJUXFeFFUhJLnRSgpKsYjWQ7ys+4hL0uKpwUPOa+KkOo1+I0HExMTFhsby9LT01lsbCwTCF5+FM/CwoJFRkZWbrdnzx6WnZ3NSkpKWFZWFvP19X3t42sberOYhoaGpv6jkDeLeaGPjxJCSP3V9NxJp6EmhBANR0VACCEajoqAEEI0HBUBIYRoOCoCQgjRcFQEhBCi4agICCFEw6nl9whyc3Nx586dBj3W1NQUeXl5ck6k2mjNmoHWrBkas2YbGxuYmZlV+zPu33ZT5mjit5JpzZoxtGbNGEWsmV4aIoQQDUdFQAghGk7jimDLli28IygdrVkz0Jo1gyLWrJZvFhNCCJEfjTsiIIQQ8ioqAkII0XBNtgg8PT2RlpaGjIwMBAUFVbvNmjVrkJGRgeTkZHTv3l3JCeWvtjVPmTIFycnJSE5Oxvnz59GlSxcOKeWrLn/OANCrVy+UlZVh/PjxSkwnf3VZ78CBA3H58mVcvXoVp0+fVm5ABahtzS1btkRERASSkpJw9epVzJgxQ/kh5Sw0NBQ5OTlISUmpcRt5P39x/1ysvEdbW5tJJBJmZ2fH9PT0WFJSEnN0dHxlGy8vL3b8+HEGgPXu3ZslJiZyz63oNbu6ujJjY2MGgA0bNkwj1vzPdnFxcSwyMpKNHz+ee25FrrdVq1bs2rVrzNramgFgbdq04Z5b0Wv+9NNP2YoVKxgAZmpqyvLz85menh737I2ZAQMGsO7du7OUlJRqfy7v568meUTg4uICiUSCzMxMlJaWYt++ffD29n5lG29vb+zYsQMAcOHCBRgbG0MoFPKIKxd1WXNCQgIePXoEAEhMTISVlRWHpPJTlzUDwIIFC3Do0CHk5uZySCk/dVnvlClTcPjwYWRlZQEAHjx4wCOq3NRlzYwxtGjRAgDQvHlzFBQUoKysjEdcuTl79iwKCgpq/Lm8n7+aZBFYWlpW/kMAAKlUCktLy3pvo07qux4/Pz+cOHFCGdEUpi5rbtu2LcaOHYtNmzYpO57c1WW9Dg4OEAgEOHXqFP766y9MmzZN2THlqi5rXrduHRwdHZGdnY2UlBQsWrQIjDFlR1UqeT9/6cojlKrR0tKqct9//2LUZRt1Up/1DBo0CH5+fujfv7+iYylUXdb8008/ISgoCBUVFcqKpTB1Wa+uri569uyJIUOGwMjICAkJCUhMTERGRoayYspVXdbs6emJpKQkuLm5oX379oiJiUHXrl3x5MkTZcVUOnk/fzXJIpBKpbC2tq68bWVlhezs7Hpvo07qup7OnTtj69at8PLyeu2hpzqoy5p79eqFffv2AXh5sq7hw4ejrKwM4eHhSs0qD3X9e52Xl4fnz5/j+fPniI+PR9euXdW2COqy5pkzZ2LFihUAgJs3byIzMxMdO3aEWCxWalZlUsTzF/c3RuQ9Ojo67ObNm8zW1rbyDSYnJ6dXthk+fPgrb7ZcuHCBe25Fr9na2pplZGQwV1dX7nmVteZ/z/bt29X6zeK6rLdjx44sNjaW6ejoMCMjI5aSksLefvtt7tkVueYNGzaw5cuXMwDMzMyMSaVS1rp1a+7ZGzs2NjY1vlmsgOcv/gtWxHh5ebEbN24wiUTCli5dygCwwMBAFhgYWLnNunXrmEQiYVeuXGE9e/bknlnRaxaJRKygoIBdvnyZXb58uUmcubEuf87/jLoXQV3Xu2TJEnbt2jWWkpLCFi1axD2zotdsYWHBoqKi2JUrV1hKSgp7//33uWdu7OzZs4dlZ2ezkpISlpWVxXx9fRX6/EWnmCCEEA3XJD81RAghpO6oCAghRMNRERBCiIajIiCEEA1HRUAIIRqOioAQQjQcFQEhhGi4/weUVmSuAUf7ZAAAAABJRU5ErkJggg==", "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2022-04-01T21:45:38.729824\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.4.3, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt #we plot the result\n", "\n", "plt.plot(x,u)" ] }, { "cell_type": "code", "execution_count": 115, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAasAAAEgCAYAAAAKZlx2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACA+klEQVR4nO29eZgc1XU2/k4v0/sy+6oNMRISGCSwBBhwjEEIgY0AL59IbBFIQLbRz7FNHIRjx4njj2DyKRg7QIiCPwsbo2AWMw4MSODdZhkhoQVpkEbrjGZfe+/p7fdH9a3unt5qufdWj756n6cfmOmqmquqW/ec8573nFsFIAUdOnTo0KGjgmHQegA6dOjQoUNHOejGSocOHTp0VDx0Y6VDhw4dOioeurHSoUOHDh0VD91Y6dChQ4eOiodurHTo0KFDR8VDN1Y6dOjQoaPioRsrHTp06NBR8TBpPQAd/2+hoaEBV1xxBVpbWxEOh3Hw4EHs3r0bqZRem65Dh47iqILewUIHB3zsYx/Dli1bUFtbi71792JkZARWqxVLlizB4sWL8dxzz2Hr1q3w+/1aD1WHDh0VCN1Y6eCChx56CD/84Q/R19eX953RaMQnPvEJGI1GvPDCCxqMTocOHZUO3Vjp0KFDh46Khy6w0MEVTz31FNxut/jzggUL8Prrr2s4Ih06dMwF6MZKB1f84Q9/wNtvv41169bhr//6r7Fz5058//vf13pYOnTomANI6R/9w/NzxRVXpGZmZlIDAwOppqYmzcdD89PV1aX5GPSP/jkbP3pkpYMrPve5z+FHP/oRNm7ciB//+Md45ZVXcOGFF2o9LGqor6/Xegg6dJy10Nxi6p//dz4vvvhiqqGhQfx51apVqT179mg+Llqf7u5uzcegf/TP2fjR1YA6NIfZbEYsFtN6GFTQ3d2NVatWaT0MHTrOOug0oA4u+Pu//3vU1NQU/C4Wi+Hqq6/GjTfeyHlUOnTomCvQ2y3p4IIDBw7gl7/8JSKRCPbs2YPR0VFYrVZ0dHRgxYoVeP311/HAAw9oPUwdOnRUKHRjpYMLPv3pT+PKK6/E17/+dYyMjKClpQU+nw8//elPcffddyMSiWg9RB06dFQwdGOlgwsuueQSzJ8/H3/xF3+Bq6++Ouc7m82mGysdOnSUhG6sdHDBf/zHf+DVV1/FOeecg927d4u/r6qqQiqVwuLFizUcnQ4dOuYCNJck6p//dz6PPfaY5mNg+dGl6/pH/7D56GpAHVzxpS99icl1165di56eHhw9ehT33Xdf3vdLly7Fn/70J0QiEdx7772Szq2pqcHOnTtx5MgR7Ny5E16vl8nYdejQIQ2aW0yaH7PVkmpd2qH5OGZ/rvnr21Mrrr9W83Fkfy5ae03qwjVXaz6O7E/jogWpz/7TN1Imi0XyOQaDIdXb25tatGhRymw2p957773UsmXLco5paGhIffjDH05997vfTd17772Szv3e976Xuu+++1IAUt95+v+mvrf9ybJj4RlZrVh7Teq6L/6V5s8s+1PX3paqtlk1H0f2p6alOfW5h76Tsrldmo+FfEwWS6qmpVnzccylz1kXWa3ZdCfufe4pfPof7oOpulrr4QAAGhbOxw1/8wV8/l//GReuubr8CRxQVVWFW+7/Gv78X76NhoXztR6OiD/beBsuvfWTuPzT6yWfs3r1avT29uLEiROIxWLYsWMH1q/PPX90dBS7d+/OKz4ude769euxfft2nLv6EjguPA+uiy9AlaEyXpnlf3YlPv9/vou1X/przL/wfK2HAwCwupz4+otP497nf1IxYwKAj2y4FSvXrcGHb7pB66GIuP6eu/DNnS/irscfhs3t0no4AIAll6/C136+HV/+6TZ0XFZ5he2V8eZRxHlXXIbg1DQu/8zNuPwzN2s9HADAhddmDNT5H7tKw5Fk0L78PLjqamG2WPDZf7wfVVVVWg8JRpNJNOZX3/k5yc5GW1tbzqaO/f39aGtrU31uU1MThkdGcNPffhnxmRmEkUTHpR+W+s9hivOuvAzJRALRUAgfu/3PtR4OAGDhRRfAbLXA4fHgrx/dWhFzqqqqCiuuvxYAsPqWT2g8mgyWXL4KU8MjOO/Ky7By3RqthwMAuPjGtaif346mxYtwaQXdK4Kzylg5arxoW7YEv33qGUwPj6Jt2VKthwQAuOi6j+PE3v3Y+8pOdFxWGYvdso9+BMlkEq/84D9wziUrMO9Dy7UeEpZecRnsHjd+/aOfwtPYgEtv/aSk8wotiqlUisq5bectQduyJfjF974PS5Wh4JjuuusudHd3o7u7m1sj2/bl5+H4nn3Y3dmF8668rCIivoUrL0QiHsdrjz8Jh9eDmrYWrYeEBRd9CLWtLTi2ey9al5yL9uXnaT0kWF1OtCw5F2899xKCk1NoO2+J1kMCAMy7YDl639mDD/70NuZfeIHWw8mD9jOcIs5dfQkA4OjbuzF07DiaO87ReERAXXsb2pYtwf5dv8aRt3bD09iApnMWaj0sLLvqIzh94H3se+0NAKiIMa1cdy2CU9Po+uETGDvdj8WrLpZ0Xn9/P+bNmyf+3N7ejoGBAdXnDg8PY+mHVwIApo6dQiNMWHL56rxrbNu2DatWrcKqVaswNjYm6e+qgcFkROvSc9F/qAd97x+GxW5H3TxpkSRLLFp5EQY+OIpT+w8CAFo6tC9HWHnDGsQiUTy95duIRaIVEV0tWnEhDAYDju/eizM9R9C2XHun2uKwo3HRApw+eAgn9x1AXXsrXHW1Wg8rB2eVseq49MMI+wPof78HQ8dOoGnRQs2piAUXCdz9B396G8ff3QsAmvP5ztoazP/Qchz+3Z8wcWYQ8ZkZNC5aoOmYAOLZvYtEPI6Rk6dQP79d0nnd3d3o6OjAwoULYTabsWHDBnR2dqo+t7OzE3924zrEolHctGYNTh7qgd3j1jzH0HTOIpgtFvQf+gADPUcBQHPv3GAyYv4Fy3Fiz34M954AADSfq72zeN6Vl6Hnj29hengUh373Ryz76Ee0HhLO+fAKxGMxnDpwCP2HP0BLx2IYTdqWvLYvPw8GgwF9Bw7h9L73AQALLqqs6OqsMlYtHYvRf6gHyUQCQ0ePo9pmRW1bq6ZjqpvXjmQyifG+M5gYGEQykUBtq7b0CFnYjr/7HpKJBEZP9WlurAwmI2pbWzB68jQAYOxUv2RjlUgksHnzZrz22ms4fPgwnn32WRw6dAibNm3Cpk2bAAj5p76+Pnzta1/DN7/5TfT19cHlchU9FwAefPBBNC9ehDqrHWuuXYMXnv4ZAKBWY3prXprK6j8kOGWJWBytSzs0HVPbeUtRbbPixHv7EQ2FMN4/oHlkZTAaUdPSjKHe4wCAkROn4G1qhMFk1HRc51y8An0HDyMejeLM4SMwmc1oWrxI0zHNT6cBTh88jP7DHyAei2FBBYlkgLOsg0VNazN6/vAWAGDomDBBm89dhPH+M5qNqa69Db6RUcRnZgAAU8MjqNHYWJFcArkvIydOab6w1DQ3w2g2YayvHwAw1tcPi90OV30d/GPjZc/v6upCV1dXzu+eeOIJ8f+Hh4dz6L5y5wLAxMQETk+M4fi77+Gu+/9JNAi1ba04c/iI5H8bbbSffx4igSDGTvUhlUph6NhxtJ2nrbFauOJDAICTe/cDAIaOHtM8svI0NsBoMmHizCAAYOLMIAxGI7xNjeLvtEDT4kV495fCfOs//AEAoH3ZUgx8cFSzMc2/YDlGT/Uh7PMBAIZ7T6C5AmjcbJw1kZXRbIansQGTg0MAgOFjhIrQ9obXz2vDWF/GWE4ODKG2XVtjVdvagngsBt+oYARGTp5C3bw2TamI+gWCIRk71Zfz34YFhQ0MD1idjhzPfOKMkMuqa9c2P9S+bCn6D38gCkFGT55G/TxpUSgr1M9vR2jaB9+okLMb7D2OxoULNJ1TJAKeHBjM+a+WzqLV5YTN5cTkgLBOjZ/uRyQY1Dxv1XzuOTnGcmJgEDUtzRqOKB9njbGqaWkCANFjigZDmBgYRPO52obXdfPaMNGfSfZPnBnUnAasbW3G1OAwUskkACGyMppMmibpCeU3dlqIrEZPC8ZKy0WYUDND6RxMJBBEcGoade3aUctVBgNal3bgTNojBwQHyNvSpGl+1tvUgKnhEfHnoaPHYDSb0KAhvUyMVXZklf17LUAMwETaqU6lUhjqPY5mjWlAT1MDpoaGxZ8nB4dQ06obKyYg3hKJrAAhutKSC662WeFuqBepLUDw7gg9oRVq2lowMZChQUaOnwIANC5aqNGIBKMUDYXgH58AAEwNDiMei4kRlxYgRmnsdKYOa7z/jKbGylnjhdlqEY06IMx5s8UCZ23hzS15wNPUiOksYzWYjka1pJdrWluQTCbFRXhqaBjJZFLTiIH8bRJZAcDU0AjcDXxKHgrB5nbBYrfnOBuTg0OwOhywud2ajWs2zhpjVZv2AiazuOih3hNoXLRAsxqU2jRdNJ5FA04MCLy5p7lRkzEBAg2YfZ+IqEFLkUX9/PacBTiZSGCif0CyyIIF3A0NAIDp4VHxdxP9A+Jz1QLuRmFRI3QbAEwOCouxlp6wYKwy92n0hOAANWj4/GrbWuAbGUUiHgcAJOJx+EZGtY2syDo1mHn/fCNj8DQ1aDUkeJqEtWhqKMtYpY0pYawqAWeNsappbUEiHsf0SOaFGTvdB7PFAo9GXgvxwHOMFaEiNKICTRYL3A31GD+ToSajoRCmR0Y1zQ/NNlaAQAlqaaw8jQ2IBIOIhkLi78b7B1DT2qyZA0QMaK6x0jYXYzSZ4K6vy4msEvE4gpNTcGpYq1PT2pInpJgcGNI0Z1XT0oxYNIrA+KT4u+mRUVjsdlgcdk3G5E07ztnPTzRWFUQFnkXGqhnTw6NIJhLi7wilpNULQ3JAuQILbY0V8ZQmB3JfYt/oGFz12i0s3uamHAoXEPJWmkZWjfXwjeQW+Y73n4EpLebRZEwNdQCQMy5xYWnWxgsm0d5UVmQFAL6xcbjq67QYEgDhHZuYNc8nBrTNGdek88XZIA62VnPKm543s3NWACpKZHFWGavZi13GWGnD5de2tiDsD4hyUECQricTCc1a0ZAXdeLMrHs1NgFXnTYLi8VuR7XNCv9orkR9vO8MLHa7Zs/P01CfE6kDEMUyWuWtCEvgG8/cq0ggiLDPr5kX7G3K98wB4f1za2SsDEYjPE0NObkhQDBWnqYGGIza1FrVtOSvUz5irJq0SQ14mxoFijTr/QtOTmEmHNEjKxaoLRDyB9LGSqu2IYVqhJLxBKaGR7SLrIhCapbH6R+f0PA+CX83ewEGIN47Z60243I3NuTQbQBE+lQr+bq7sQH+8Qkk44mc308ODmnmBYs5j1nGKjA+oZ2j0ZRbY0UweWYQRpNJsyimprU5z4ASh0grkYW3uRH+sXFRHUwwNTSsR1a0YTAZc2qsCAITAi+s2SJcVwv/xETe7ycHhlDTps0kIDVW/lmLsH9sHM66Gk3kz4Qq8o/l3ivx+WmkcvMUoAGnhoaRiMdRO0+byMpdX5dnQIH089PsPhEhSq6x8o2NaxatE2dwNt09QShTDZgNU3U13PV1omydgDxPrQyop6kxR1xBMDkwqEdWtOFpFML62R7LTDiCaCikWc7KWVuTk0gl0LLWylVXi8DEZF5Xcv/4BIwmE+we/lLVjLGaFVlpmHO0e9wwVVfn0YDJeAJTQyOo06iNl7uxHr4CHT2CU9Nw1Hr5DwiAp7kRkWAQkUAw5/f+sQlY7DZY7PyFA7NrrAi0FDiR3N5soz4TjiDs88PTqFFk1dSYFxUDwPTomGY0biGcFcaKRE6F2vL4xyc088yJYZiN6WGhrkKLKMZR40VwYirv9wENDYM7TQMS4ySOKX3vtIgYiJR4ulAUM65dFOOuz4/2ACAwMQVnjTZj8s6SrROINK4Gc8qV3qpl9vMjIgKvBpJsZ40XAAo6sNMjo3BrKLDIFlcQBCem4EiPuRJwVhgrR/olDUxO5X0XGJ/U5GUhUcrsBRgQxmk0mWB1ObmPy1njRWAy/2XxpcephSflqqtDIh5HaGo65/dhnx+JWFwTwyBKxAsYBq1e4iqDAa762oI0YGByElanA0azmfu4PE0NedECkHE+3BqoTJ21XkRDIcSj0ZzfJ2IxhHw+TRzYUuuUb3RMkxIbm9uNapu1oLMRmJiE2WLRTFI/G2eFsXKm6Y9ggUVYK+GAI/0yFPKiglNTwjFeD88hCX+z1otgQaOuXWTlqq9DYDyfmgSEF0aL5yfmYUbyF+Hg1LQmxspZWwOD0VjQWAXTht6pARVYjEbypwUzWsjXHV4vAgUYBEA7Z6PUOjU9MqpJYbArLYAp5lQD0Cxin42zw1iR8LoQvTUxqYlnTjy3QpMgOCksLJq8MDU1Be8TyYNoUWvlqq/NUwIS+McntImsxE4R+eMKTE6Kc44nxBqrQsYqTZk6vF6eQ0pHe3V5uT0gI5jRxFjVekWncDY0czZKrFPTI2Nw1ddxLzYn96GQASUMjFa50Nk4K4yVo6YGsUgUM+Fw3nf+8Qk4arzcJwGJUArlrEhkw3thMZrNsDodBWnAiD+A+MyMJlFMqW1AtHI23PV1CE5NIxGL5X0XnJiCqbqau3CgUPcKAtEL5ryw2D1uGE0mMTLPRnBqGslEQhtmw+spyCAA2jkbpdYp38gojCYT97lO1qCCjv64tmrc2TgrjJWztnAeBhDoLYPBAEcNX8rNWSqyIjQg9zF5hb9f5CUWKFMtcla1ebJ1Aq1oQEdNYboUyHp+nA0DUYsVzKMRB4gzZUOobMIWZCOVTGr2/Jw1NQXHBGhJA9YUXadE5SvnOSWuCQWiUK3mVDGcFcbKUeMtmLQEMpOA9wvjEiOrEjQg58iKcM/FuHyhiwXf+1RVVQVXbW1Bow5oRwM6PB6Epn0FvyP3j7d3TsQvhSjToJhf8HIcEWD3pI3VVGHD4B+b0IYGrPEUd2AntctZFVunyP0j95MXyBpUyLCT+6dFHrQQzgpj5aypKSjHBrQrDHbWCiF/NBjK+24mHEYsGuUeWZXip4F0ZMU5Z2X3uGE0m0rSgNU2K6ptNr7j8rrLR1a8nY26WoFam9W9AgBC0z4kEwnui7DD607//cLGyjc+zn1OmSwWWOz2opFVYHISJrOZuxrXUeMtmBYAIDpGvOscHbVeRIJBcSfzbMQiUURD4YqRr58dxkpKeK1BZFWoewVBcGqa/2KX9pCKRqFj45oYdaBwbi/797zH5fB6ECyyAGuZH5ot7ydIpVKaCAfs3tKRVUADNa6zjFMW1CgyLuVUE2PPWyFcKrcHaJczLoSzwliV8lgyXDDfG16sewVBcHJKg8iqDA04wV+MQjZ3C/v8Bb8nNCrvHnN2jwehqcI0IFlweDsbdo+7KDUJCHOK9wLsSNNWxe6VX4PO6+S9KmpANZJkO4rUOALaRValcnsAmVO6saICs9UCi91W1DvQSuXmqiuehwE0iqxqvEjE44j4ixiGdMslnt6dze0CAISKGCstnA2TxYJqm7XoYhcNhRCfmeEusLC5XQj5ShgrDVouOWo8iMdiOXt+ZcM3NgGT2cx1x9lyudmgBpLscutULBJFLBqFnfPOvI4S4jRAoEx16TollJuYgDaFwc7amqLRHpCOrHiH/LVeBKemCxbfApl7yNdYCXmDcJFFWCoNuHbtWvT09ODo0aO47777Ch7zyCOP4OjRo9i3bx9WrlwJAFiyZAn27t0rfqanp/Hlr30VAPDxq65Cf3+/+N26desy49LA4ywXWQUmJjUYk6coNQlkis15drEgkVWxRZhEEjyjUCnrVGjKJ9KqvCCFBqwU6bpJ6wGoRTnRAJB+iXnnPGo8RYsSAW0KE4WQf6ro92Qh5OkF29ORVdgfKPi9qLwr8cIYDAY8+uijWLNmDfr7+9Hd3Y3Ozk4cPnxYPGbdunXo6OhAR0cHLr30Ujz++OO47LLLcOTIEdFwGQwGnDlzBr/50x/xF3++HqmZGB5++GFs3bo1729q4WxIoQG5j8nrKRqBAhkqjuecInR3KYFF9nF8xuQFUHqdCvl8lUcDVlB/wLkfWZURDQD85c9mqwWm6uqieRhAWFhsbhfXTeCcJXJ7QCZvxPOFsbqIsSp8rxKxGML+QEnKdPXq1ejt7cWJEycQi8WwY8cOrF+/PueY9evX46mnngIAvP322/B6vWhuzt3+4JprrsGxY8fgCwuUVmomvyCYIDg5xZUeqTIYYHeXiawmp2D3erjmHB3e4hJ/IOOE2Dgq7xw1npJ0d0blxs+wS1mnglPTsHv5vXtmq0B3l6YBp4T+gBp0zp+NOW+syokGAKESmycNSLzIYnkYQJiYBoNBzNnwQKlCVyAzXp5jsrtdiASDBeXYBJFAAFaXo+j3bW1t6OvrE3/u7+9HW1ub7GM2bNiAZ555RvQkU9EZbN68Gfv27cOTTz4Jb5bBDExOiQb0rrvuQnd3N7q7u1Ffz6YZKVnsy0VWBoOBq7Nh97hLRlaRQNpYufkZK2dNDULTvqJ0N8BfOCBlnQpN+7jmrJxlIlBA250PZuMsMFZE+TNV9JjQtA82F98FGCiucAOyWy7x9O5qSnp2oiKJo7GyuV2IFKEACcL+QMnnV2irldkLVbljzGYzbrrpJvz85z8XF/sfb/svLF68GCtWrMDg4GAOHZitvNu2bRtWrVqFVatWYWwsv7sEDZAxlRRYaFAY7PCWzlmRiNnq5BhZeT0lGQRAWC94RlaS1qmpaa7rgZQxZYqV+e9zNxtz3ljZ3W4kk0lEZ238lo2w3w+L3QajiU+KTlxYSnnBU3yb2RpMRsELLvESEy+Y58S0uVwlI1BAeH6laKT+/n7MmzdP/Lm9vR0DAwOyjlm3bh327NmDkZERccE43XsMyWQSqVQK27Ztw+rVq8XjA2kat5LmFHFEeOYY7F5P0YJgAIj4hfeSp7NIhESlEPb5+TIIHjeSiUTJdYp3zkpkgErMKUKl8oyMi2HOGyvimZcK+Qlvzqti3SZGVuW9YF7ydUIvBEtMzFQyiZDPx/UltrldRfNVBBFf6ciqu7sbHR0dWLhwIcxmMzZs2IDOzs6cYzo7O7Fx40YAwKWXXorp6WkMDWV2lr7tttvwzDPPABAW4EgwiIYsSu+WW27BwYMHxZ9JrRUv9ZYkB4izsbLY7TCZzQgWqbECgPjMDGKRKN+clbc03Q0IlDdPA2pzuRAus06FpqZhqq7m1q2FPBPipBZCZu3kd6+KYc6rAW0uZ0lqBMhQETaXs+wkpgF7mdohgH8zW2KoSxlQ4Xs/58jKicnBoZLHhAMBNLvOKfp9IpHA5s2b8dprr8FoNOJHP/oRDh06hE2bNgEAnnjiCbzyyiu44YYb0Nvbi1AohDvuuCMzBpsNa9asEY93pAuCH3roIaxYsQKpVAonT54UvwcydJzD6ynaKoom5ERWvHIxRAxQikYChOfHs7VROdUkkI7WOedmyzll2YXBhTqz0wahZsO+4saKrGE8UwPFMPeNldtVVPZMQB4GL0+qXFcGAKI3yiuyIv/2UhMTSHucnCOrgQ+OljwmUiZnBQBdXV3o6urK+d0TTzyR8/PmzZsLnhsOh3OEEXavG8HpaWzceEfB48mYAMDqKC78oAlirMJSIitOKkVCl5bKWQHk+fEzVjaXs6xhCPtKU8u0YXU5S64HQMZYObyegtvM0wb595e6V5EsR19rzHka0O52l50EvHlXu0fIo5USDsSjUURDIW6RVWZiljHsnBVJwvMrL7CwOh0FRRIs4ChT6AoAkaCQe7A4ORkr4gCVeH6JWAyRYFBsgcR8TGLH9fLROq/FzmSxwGyxlF0Twr4AzBYLTBYLl3HZ3e6y716Qc8slq8uJZDKJmVDxKG4mHEEiFq8IGnDOGyspHgtv3lVKHg3g23JJSh4NECIrXi+LwWiE1emQRE0ajEZU2/lw+fYyCjcgE1nZOBkrm0dwypKJ4hJ/AIgEgiVl/jQhRlYlBBZAuvSAkxpQdMrKOUCc6S05kRWvPKjN5UQkUH6dCvl8FUEDznljZS/TLw3IzVnxgJQxAYQe4bTYlekUQRCa5iewsKYX+vJqQL40rtBxvfTzI5EVr1yM3eOWNKfCPj83w1CuYSxB2B/gNqekUFtAxmnjuSaUNVZEJs6J2bA6naJasxR407jFMOeNlU3CJOCfsyo/JiCdeObtcUrh8nktLBKoLSCjVuJlGIRorxy1nDZWHHNW5UQDgBBZ2XjNKQm5WSBN43IzCtLmlOgAcTIMNpcEY8WZBrS5y+f2gPTz0yMrdcjw06Un5kw4jEQ8ztGLkriw+IP8DKjLhXgshlgkWvK4sM8Hk9nMRT6boWykRlbsn5/FbofBaCxbqBwNhZBMJvlGVpKMVQBWTrlZq8Mh3Idy1KQ/wM2AWss0RiYQu7VweH6m6mqYrZayBjQ+M4NoKMyt5ZLV6USkRN0XQbk6R16Y08ZK6mIHpF8YbhGDxMjK7+eWX7BJ4MwBvl0s7J7SfQEJwuLCwn5M4pwqUXtCEA0EYXHw6Zkm1QEK8zQMTofExS4g9stkDVH1WjaySs8pD785JY3G5VcYLEU1CaTLWThvXVIIc9pY2SWKBgDSsodfUXC5PAyQToZzpAElGStiGDi8MGIT2wqKrKwSVZOAYNC4RetSI6u0cpIHrE4HosHC+1hlQ2y5xMExy3Txl+YA8XDKiJNcLloH0qIrTmpOqTkrnjRuKcxpYyWlYSyBEMXwowGlRnvcFhYJ9WgA58hKQvE0wLcZKnEepCws0WCIn7MhUbRTro8iTVgkRlbi8+MwLqtUNaBoQHnMc+nrVGiaY2QlOWel04CqIVXhRo7h8bJY7HYYzSZJ0V4kEIDRZOKSH7K7XEW3TMiGSLlxCPszHmc5MUNaYMHBMEitRwPSzgYHgYXFYYfRZJKcszKaTTBb2dcPWR2Okq16CIjh4PP8XOnaoOLbuwBAMp5ANBTi4pSJeTQJ718kEORSu1dVVZV+fhIiK5+fa01aMcxxYyUtmSocw8c7ELdpl5hfAPio3KR0+gD47mlFRB8z4UjJ4xLxOKKhMFfPXNIiXGbrEloQu1dIYhD4GQY5OSuAD40rNQ8D8HNgpezCQBANhrg4QNV2GwxGo7Q55eP3/EphThsrOZMgzKlxpU0itQXwLSy1upwVRwNKFaIA/Go9CC0riQYMBLksLDaJuT0ga07xUE467GK9WSlkaFxOc0rCswP4lWlkRB8Snl+Qj2iHiHCksVKV0XJpThsrqXU65BgeEYwsAxrgV1gq1eOMhkKCzJ/LSyxN9AHwyzlKVZMBwvPjMaZMtCchiuFYk2aTIX0GwCU/K8cB4s22lMujAfwiK1I3JY3G5b8payHMcWOV3mW2TJ0HwG9PKzk0kljsypiyMVVXS+qXRhANhbh4d1Kq+gkigSAnA+pAfGYG8ZmZssdGA0EuCzBZvCSJGcj+URxoQItUNSDHonxZNCDHyEpKHg0QjJXRbGIu85eTmxUjK91YKYfN6ZRE1wBZSXrGnhRZ5CMSXmJelI0czw4QXhguVISE5p4EYb+fj2FwSqNLASGKMVssMJrNbMfkIsaqciibapsNBoNBkgHlWZQvZedpAl47DMiiu0kbL8ZzXY7qlXe7s2KY08bK4rDLWuwA9i+xxS4s8qV2BM2MiY8BldpqiSAaDIn/DpYQaMDyQhSAXzLc5pLuAEV5LSwkspKRH2I9p8TcnoQxAaTBLh8aV86aUGnRHolUWb9/NhkKRT1nRQFWp9DuRQp4URHkJZYyLjHxzHixkyPxB4QFiFd+QYoQBeCXX7C6nJK6VwBAmFN/wIwXLIMG5GSspDhlAL9aHSmbsRKQOVVlYLsMyoqsAnwcIEITS5Ou8xPIlMKcNlYWh13Gy8KnsNTisCOZSJSVYwPpvWLiceY5K0WRFQ8aUIYXHAnw6aMotV8akB3FsDZWDsRjMWl5tHSvPtZzypI20FINO48+mBZHuq+jRLpb3JOM8VyXM8+Jk8t6THI6tSTS5SU6DagCVodDUm4I4FexbnVISzoT8BAOyJE+A3xowGqbLV08LV0NaDSbUG2zMh2XHIVixgtmH8VIdcoAPpSbksiKV7Qg1YCSsTOPjF3SiqeBzJyyMB6TzelELBKVJPoAKqOLxdw2VhLVSAC/wkSLwy6ZmgT4vMRWuZEVBzWg1B5uBLyKXa0uBZEV82S4Q/ICDPBZWOTmrHjsaWWRaUDF/BDjuW6x2yU71WIelHVkJbHVEgHPrYOKYU4bK4vdLn1hEZOEjF8YGWMCBHqE9QJMepNJyXkApMEu4wVYRj0aAEQ4bekgJxkeCQgLEHMD6nAgGpARrXOYU6KxkmHYeTiKgAzRB0eBjFSnOsLJgEqtkSPguYFmMcxtYyVDYCHUObCXz8qJ9oD0/kOscx4yaocAPjSgnOJpgI981mA0ynOAeEVWMkQfAJ82UJnaL+nMBg+jDkDy+8dDeWcwGWG2WiSvU5nIijENKIPuBnQaUBXk1HkQhP3sQ1m5NKDgcVZWtBcNhWA0mZg2rrTJpQGJYSgikFm7di16enpw9OhR3HfffQWPeeSRR3D06FHs27cPK1euFH9/4sQJ7N+/H2/t7gaQqT2pqanBzp07ceTIEezcuRNerzfneryUW0JkJSdaZy/JFim3kETD7hPoboPRyG5MMmochePYPz+LXV4EOhMKC+expgGdTsl5NIBfOUspzFljRR6mnCiGx55WUjsZE4Q5bBMi14CSe8qSNyeLqVTpeqaPYv7zMxgMePTRR7Fu3TosX74ct912G5YtW5ZzzLp169DR0YGOjg7cfffdePzxx3O+v/rqq3H9Jz8BIBPFbdmyBW+88QaWLFmCN954A1u2bMk5JxGLIRaNMu8WYXFK68FHEPZzoHEdDsyEI0jGy3ePATKtxVgKB0TRh8R7xSNnZZW5TqVSKaE/IIdyFqkKRYBfOUspzFljJTfBC/DpLydfYMF+Z1erwy6PmhRfYnaT0yaXBvQRNWf+vVq9ejV6e3tx4sQJxGIx7NixA+vXr885Zv369XjqqacAAG+//Ta8Xi+am5tzxyTWngTEc7Zv3w4A2L59O26++ea8v81jSwe5+YVIgMOcckpXuAFZOWOGpSNiQb7UyIqDGlB0quU4i4EQrIyjGKvTIctY8SpnKQW2jfIYgkxMeWIGDpSbQ6bAIr3YVVVVIZVKsRmTXV4ejYciiUS4Uu9VqZxVW1sb+vr6xJ/7+/tx6aWXlj2mra0NQ0NDSKVS2LlzJyaRwD4kxL/V1NSEoaEhAMDQ0BAaGxvz/nYkEMQFF12ELd0ChVhfXy/p3yMHUjc5JAj7A8znlNTtQbLHBLDtWSiXBuQRWWUYIHk0PHMHyCW9LRWQoQFZzqlymPORlewohnHOSnadlT8Ag8HA/IWJKKABWY6J0BCpZFLS8UQgUsgzr6qqyvvd7Beq1DFXXHEFLrnkEvzDP38HALDs3A5JYwKEKKZ/eBCrVq3CqlWrMDY2JvlcKTBZLDCZzTKjGPZzyup0yGQ12LeBsjoFIZHU2qFEPI5YNMo2spIZ7QHstwkxms0wWy2yIytA2AdLK8x9YyW3WJKhx2K2WmAwGmUbUICt/NmilAZkSEXIpZGA4v0B+/v7MW/ePPHn9vZ2DAwMSD5mcHAQADCTNpzLlywFAAwPD4tUYXNzM0ZGRvL+diQQZBotWJ0KGAQO3fzly+nZqzktdnnzHGBPb8mV0wNpGpAlBS9jZwiCTLcP7fJWc9ZYKaEBWSta5GzlQMCj8ahcYxXl0IZGycIS9hUuoO7u7kZHRwcWLlwIs9mMDRs2oLOzM+eYzs5ObNy4EQBw6aWXYnp6GkNDQ7Db7XCmF3V3bQ0A4ND+/eI5t99+OwDg9ttvx0svvZT3t1l3ixD7AiqIYpjmh2RHVqROjq3AQs6YAPbCgYzoQ8b7F2IbWcltEgDwEV2Vw5zNWSmhAVl3ZlCiUOTRzFZoSyXDgJKJyXRMCoxVkcgqkUhg8+bNeO2112A0GvGjH/0Ihw4dwqZNmwAATzzxBF555RXccMMN6O3tRSgUwh133AFAyEu9+OKLAIBTxgR6MYNXfvlLAMCDDz6IZ599Fn/1V3+F06dP4zOf+Uze344EAozVZGTn4sqLrOTlrMhGowwjKwVzirUDq4gGDLAdU2aXYHmOPqBtZDVnjZUYXsukAQ0GA6ptNsyEw+zGJFNmDLBdWKrtNmWRFcuX2CG9oJsgUqIwsaurC11dXTm/e+KJJ3J+3rx5c955J06cwIoVKwAAN371S7jqLz4ryrEnJiZw7bXXlh4TcxpQgerVx8FYOR3yRANB9t3g5dYTAuzzQ3JFH4DgVLN0FEnEHZETWYW0j6zmLg3ocCARi0vuygCwFw7IraAHsnfhZPMSk+JpOWPiUZhocUjvl0YQZtz0Vwk1yVq6TgyOnNwsKaBmSgM67LK6aiQTCQ6UqXwHKBpkq7yzOuyIx2KSRR9Aek4xFceQhr9KIivdWMmGEn6aVNuzuuEWBcaKNWWjhJpMpVLsKVMlOSvGdXJya+SATLTOKgolHrYcw0AMG6uIr9pmhdFkkmVAAfb9ARUJLAJBpjVNFpndRwDhXTWZzcy2thcFFjLaLfHqBl8Kc9ZYKQn5WScJReWWnPwQ42aomaJEuQtLkHmxpGzD4GNb7Kokj0acDVbeuRLVK7mvrGTGFpl9AQlY7/Yst/YLACKMa5qUMAis99nKRFYyHCCdBlQOJSE/a0k26QMmZ8GLR6NIxOIMqUn5kRU5njUNKDuyCgRgtlpgNJvZjElm8TSQ3QaKrbGSswjPhCNIJhIMFzv5eTSguJqTFpQ4QFHGYgarkjGJTjWbe2Vzu5BMJkW6Xwp4dYMvhblrrBSG1wC7UFaJ6ANIJ1RZUZN2+QlegK1y0mg2w2Q2KzKgADvvTm7xNJDVDJURvWV1OIRN8uJxWedFQ2F21KSCEg1yPKvIqipNxcpnW4Kw2G3MGuwqcoACrNMVwg7rcjpRxKPCHNSyme2cNVZKFhbWW0ZbnQ7BY5GpNGSpSFJS50GOZx7tyaQmWT8/JdEekZSz8oKtLqfsCAYgDhBjalJuZOX3MxN9WNKUp+znx7gzgyI5PaHcWEXrCpS4gPb9AeessZK71TfAXtFisdtlhdYELGs9lMjpyfGsFjsloo/s4yvKWDGPrOyyergRsFxYLAojK5Y5KyUNYwH2+0cJOSv50R7AMF2hII8GaN95fc4aK3UCC3YepyIvOMiurkJJHo0cz9IoAAqoSTHnyOpeKYms0mpOVlGo0siKoQNEulDIbZcVYbgdjmJqkkM5i1IakN2aIH+eA9rvaTV3jZVTfuJyJhxGMpmsKM8cILw5W8NQSZX9agwowDIytikQ7bAt6rY6HLK6VxBUYmQVDYVgMAo759Ifk9J5ztgwqKABma5TSmhAxuUs5TAnjVVVVZXsdi+AUD80EwozpQHljgkQvDuWE1NJHi0aZJfzUErZRBi+xNU2q9CEWDZlE0IymWRGA8rtwUfAsr+cmjwowOb5iQZULt0dYKsQrlbiAAXY04CKIyvdWMlDtcJkKjmHVRGgEjk9wJgGdCjLo0VCIaGLvIm+Skq1F8zg+SlVTaZSKaaGXUntEJB2gBiqAZUoFCMMaVzFJRohdpFVpnuMvOdH3tdKowFZ116Ww5w0VkrrPAC2oaxSjyXCkAa0KuWnRe+OwcKioLln9vFsPHNl0R4gGFEro27iSrZSAUhkzGieK4z2ZkhkzEB5p6R7DMBWJq40Nytubc9QjavTgJwg9uBT5HGymwRydwkmYGpAVUR7ABvhgGqBBUtjpcCwh/0BdpGVggQ9wHZOKY32WOZilHSPAdiKrpRK/IVz2NLwSsekGyuZEBc7pbUCDOkRpYaBZeJZiRfMMj9E2tvIrbNKxIXGxSxzHoqpZQY5KyWbeWaPyVRdDaOJ/sYKcrecIWCpvKtE0Y6S7UEIWBkGo8kEU3W1Ium6rgZUACX90giYRjFKBRZEqsrCu7PbZe3oSsDU47TbEYtGxa045IDVC6NmYQkH2Eiy1UR7LBdhq9OhuPYLYJRzTHc3l7MLA8B2a3ulDIJwDptu/mrmVCQYFIydhb5TLQVz0lhlkuGVE8qaLBYYzSbFlA3Ajt6SG8EAWTkrVmNScJ8AdsrJTFcNBS8xIxqQRAtKvWCAnbFSSiMBQDUDYyWMSeGcYrQlh1LRByD0LGRl1AGluVnSWUOb6GpOGiul0llyDpNogcIkYGUYFHl2zA2osoWFVR9FtV4wCxpQzWKXeX4svHOH7I7rAGOnTKGQCACzzQ5VGYYQm8hK7doJsCvKL4c5aayUFiUC7CIrNWNiSQMqr1avzMhKoAHZ3Cfh+gqen5+NpJeUaMyEVXjBjKIYpUpcgKVTJn9MALvO6+rWBFZ0twpjpfE2IXPUWKnjXc0W+vVDqkJ+ht3gFVerM8xZqfKCmTkbwjUV1aQFg0JRMeU5pSbaY0oDKijIB4BkPCHkh1gYUIWqSQCIsIpi1KwJjKM9pSmU7Gvwxpw0VlanA7Go/KJEgF0omzGgleNxGk0mmC2WikvQK+mYT8BKIEOMupxtEwgye1rRpQKV1qMBGaUl7XuVyc0qjGIYOhuqIitGbIuS7jFAugVbhZVosO6sUQ5z0lgpVd0B7OqHMu1eKscwqJmYRCXFyjOvxMhK+ZjSLzFlT7gS1YBkk0nF718oxGQ7DjUCiyijHQaU7sIACDSgyWymvtGoOgeICCz0nJVkKG1rBLATDlhVRFasclZqFjtyHiuZuHI1YOFuH2vXrkVPTw+OHj2K++67r+C5jzzyCI4ePYp9+/Zh5cqVAID29nb86le/wq2f+QyaGxrw5S9/WTz+29/+Nvr7+7F3717s3bsX69atK3jdcLrRLO3ISqz9UiPaofz8LKKxki9dB9h11lA3pyrQAWLUBkqdkEinAWVD2CVY+QIMMIhinMojK1bd4JU29yRgyZurUW7NbtdjMBjw6KOPYt26dVi+fDluu+02LFu2LOeYdevWoaOjAx0dHbj77rvx+OOPAwDi8TjuvfdevPGb36C35wPcc889Oec+/PDDWLlyJVauXImurq6CYyILdyVFVqwWlsxWHEqfX5iNQEYNDchKIaxQiAJkNdhlxbYoKWdhvM9WOcxJYyXkPNTRgLRvuFWFmgwAk27wakJ+ch6TDucKOlETFOr2sXr1avT29uLEiROIxWLYsWMH1q9fn3Pe+vXr8dRTTwEA3n77bXi9XjQ3N2NoaAh79+6FxWFHyB/A4cOH0dbWJmtMJDK20TZWNhtikSiSCfnF04lYDPFYjHoUo6YvJ8BmTlVVVamrsyICGcpb26sTErFjW5QW5BNKU4+sZMDiVBFZMUo8q1GTAWx25lVLA0YC9BvsKu1ETVAoMm5ra0NfX5/4c39/f57BKXeMxW6HMZXCypUr8fbbb4u/37x5M/bt24cnn3wSXq+34JhuuuFGAMD3tm5FfX29on9XIaipRwOIYaA7p9R0jwHYCGSqbeldGFTmsVmsCWqoSXINmlBjQEmDXT1nJQNKe5MB2VXYlA1DOuRXoiYD2HicanooAmwWFnFMamncLCqpqqoq77jZz6HcMTaXA5de8mF85Stfgd/vBwA8/vjjWLx4MVasWIHBwUFs3bq14Jie+r8/AgD8y78+hLGxMZn/ouKgYqxo56xU1A4BZEx0BRYiBa80j81I5aa0ewzANrKqtDklFXPSWKktKgXY0IBKxwSwNVZqohhm0Z7Sl7hAYWJ/fz/mzZsn/tze3o6BgYGc80odYzKZ0NTaiuNHe/Hiiy+Kx4yMjCCZTCKVSmHbtm1YvXp14TFVoGcOsHE21NKALLbDUVPPBLBTuVnsyrrHAOy2LrGqnVMadl6fk8ZKaVEiwDC8VsGZA2nDQP0lVl6tTs5jp5qkJ5Dp7u5GR0cHFi5cCLPZjA0bNqCzszPnvM7OTmzcuBEAcOmll2J6ehpDQ0MAgCeffBJxpLD7rbdzzmlubhb//5ZbbsHBgwcLjonklah75hXoAFlVStdZ5GbFXqFK1wRWhkGNnJ6ZAVW5Tmm4p9WcM1YGoxHVNqviUDaR7sxMv85KuRoJYNNlWXXOikFhotIdeQkiBSLjRCKBzZs347XXXsPhw4fx7LPP4tChQ9i0aRM2bdoEAHjllVdw/Phx9Pb2Ytu2bfjSl74EALjiiivwuc9/HjAaccfnN+ZI1B966CHs378f+/btw9VXX42vfvWrRcfFKjJWSm0B6cJS2g4QKciPxZSNKUR/Oxwxj6ZCtAMwotwUqwHZUZNq5hSrpr9SQH+zG8bI5DyUGwYmiWeFW3EQsMkv2BGfmVHU6YOMyepwoKqqSnEuLn9MaqO9wl5wV1dXnrT8iSeeyPl58+bNedf74x//CIfbhf/95uv4/r9txe+e2iF+RyIxSeMKMXh+djsmzgwqPj8SDMHT3ERxROpYDSA3Mo5FolTGJDplitkW+pGV0WyGyWxW7JSxUt5ZHHZMDCifU9FQCDUtzeUPZIA5F1mp2SWYIMLCMDgdivMwAJtu4mo6fQCZhYWorWhATSfq7PNovsRqI1AgXT/EYGFRqi4F2MwpNdQWkBUZ22g+P+U1jgAbMYNaupuV8k51HlTPWUmHWoUbwEjlpiKZCqQpGxYTU03ITwwDxXHR6KqRfR0qY1JZj0bOZRFZqaGWWakBaUVWtJDZnqdyCnDVOmUAm+enptUZoKsBZYFGZMUq8azK42TQC0y18kdsA0VxYVFbqJyONGiKUdRSk+Rc+g6Q8uJpgOE8VzEmst0JC8OgvKsG/d6OaiX+QJqGp+zAqinIB9h1+5CCOWesLCqlswC7/JDSZCqQ1QuMasRAibKhahjSxdMKOlEDQCqZTFNuDCgbVdE63cQz6aigVrlFXcygMrJioca1OBxIxOOIR5XlwJLxBGIRuluXqKUBAfoCJ7UF+WRMZiv9LZakYO4ZKyr5Bboei5hMVSmwAOh7nGrpUnIdWiCqSTWCDdo0Lo05FaFc7EqLmgToPz9VThlR3lE2DGruE0BfjUuLBqQZxahpYkug5W7Bc85YWSmF11Q5c1E6q8JjEWs9KEYMqmu/6Hd+Vls7BNCv1VErpwfSY2IQgdKJjOkueDQWu2rKNK4apgWg3w2eDrVMN7KiZUABbXYLnnPGisYNp11FTyvaA2jTgOoMg5h4pnyv1Bor2jSu2q4aAH0HSIys1KgBGSwsqhP0jKJ1KnOKZhRDITKOUB6TSE2qoXEZbbEkBXPOWIlRjGpFCz3Khka0x4bLV6kmY9D0V61CEUhvQ15pNGAoBLPFAuS3IFQ3pgqicQ0mIf+lbk6xUQNSoQEZsC1qc+u0nVdAJQ3IqNuHFMw9Y+VwIBoKK9o2gSAaDIqdMGiAbnhNz5OqttsqNOdBwwumO6ZkIqGqSJX8mwwGOq9UxjNXs9jR7YJA6EQ1zy8ZTyAWpStmsDgcqu4TwC4/pJoGZKBQpLEmaKEInHPGiopnTlnlRqWrBmWPJaP8UX6vSM87qsWSKhWKAP0+ijTyaOT8KgMdlRSdxY6us0FD4UbOp88gVFZkZbHbFe9FJo6JcjkLlRSKTgNKh9VhV8W5AvSpCLUNY4HsSUDHMNBY7Mj5lUYDUhfIONTVDgHATIhyZFWBxiqjJlP//lXTpOFVCokA+jVNail4gH45C61oL/taPDHnjJXF4VAlxwbYvcRqJgHtXmA0RAMAA4+TghfMQrqu2jMPkciKrrFS1XS0Ap0ycj7tXIxaw0B7o1ErDQaICJwoGVG1BfnZ52rRxWLuGSsnhciKMu8qCixUvDDJRAIz4Qj1MVWc8o4S5Wax2wtuqKjlmAD6OSs1vQFnKHf7oBqtU5bTqx5TKER1a3uLysbWAIOco8qCfGFMOg0oGWp2CSbI3HDKlJuKhQWgG8XQKAAE6NIjard3yR4TAGpUEi1qEqAYWVHIeRAHiPY8V09v0auTo9GVAWAwp5wORFSyGrT3tKJRkJ+IxwWBjAZb2885Y0XLiyLXogGy2KWSSVXXoZmLoZYMp0i5kYVAtcyYQX5P/X0SHBWaNKBaowDQfX6ZaL1yclZq+wISiIaBYhSjep5TFl3RYBAANu3qpGDOGSu1vckA+jVNagslCWjKZ2l5wTS5fCsFObZwPt1iVyo0YNqLNhjpvFI0ch4AXQeIJg1I69nR6B4DMFBOUqSWaa4JVNYpjXYLnnvGSmXXZyAjE68kLwqgTANSqIkB6E5MC4WC7uzzqRnRCqUBK80LpmWsaHaQodEpIvt8qmwLhftErkUDNOhuQHBg9XZLZWA0m2GqrlYdWZEEIzUvyqk+jwZU5sJCM2dFQ+EG0JfPqm3OCmRq0mjSgFQcIIrdPqyO9Jb2CneeJqDqANHKzTKglmnl1mmKrqg5QHpRcGlk8jDqJgHZhbPiuGDKL3EymVSl/AHo7qpMM48G0DFWRpNJcIAovcS01IBq9x0ioE0D0pjnM6Ewta1LqEdWFOZ6lcFAiVqmX85SaeuUHMwpY0UrWiDXoBnF0IqsaEYxamTPBNFgUFjQLZWzsNDMOVKdU+EwtciKtBVTC5rdPmjOc4COYaAhxwYozylKQiJh7za6awItB0hvt1QGNHbfJKCqvKNQQQ9QFjNQDPmF69FYWCovZ0WjBQ0BzciKas6qwoRENNud0aMB6XWLoDmnaK4JQvF0ZTn6cjCnjFWmkzGdUJaax0KRBqRVmEgtmUqxJo3WS1zIgK5duxY9PT04evQo7rvvvoLnPfLIIzh69Cj27duHlStXpsck/LseeuBf8s6tqanBzp07ceTIEezcuRNer7fsuCotZ0WbWqbpANEYF22FKQ1JPS0GgVyj8pwNul1tpGJOGasMZVNZYgaaAguA0gtDmbKh4nFS2KIAECifZDIpGhqDwYBHH30U69atw/Lly3Hbbbdh2bJlOeesW7cOHR0d6OjowN13343HH38cAGBLOyz/+M1v5Z27ZcsWvPHGG1iyZAneeOMNbNmypeS4KjZnRanbBzUaiWLOUaQBVVKmRCBDZ0z0GKBIMEjFqSY5Qmo0oNNBrYOMVHAxVn/30jNY9/9tUn0dGvtGEdDyWEiCno4XRagI9ZPTSqHdC0B/YUnE44hHlW/FQZC9M+/q1avR29uLEydOIBaLYceOHVi/fn3O8evXr8dTTz0FAHj77bfh9XrR3NyMCy8WIqxTx47nnbt+/Xps374dALB9+3bcfPPNJccUDQVRRaHOilZXBoCuA0SjxhGgP6eioZCqrgwEUUq7PWdqvyqHcqOd7weEecoTVQCKPuWuri7U19er/iMtS84F4nEMHj+p6jqOGi9qWpsxeKQXiZg6+WxtWyuq7VYMHT1e8rj6+nqMjY0V/d5gNKL1vA5MDQ4jMDGpakw2twt189ow3HsCMZULetPiRTBWVWGgt/S/rxyqbVY0nrMQY6f6EQkEVF3L29wEu9eNgZ6jZY8td99blpyLSCCIyYFB1NTUwO1249SpUwCA2tpaOJ1OnD59Wjz+3HPPxdDQEALpf8OSJUvQ398Ph8cNS61XvOfZ565YsQLvvfeeeI3ZP2ePtaGhAZGqFNwWG/bt3SvxjhSGwWRC69JzMTkwhODklKprie/MB71lJefl7nnrknMRDgQwOTCkakxmiwVN5y7CeN8ZhH1+VdeqaW2Gw+1Gf88RVdcBcueUGthcLtTNb8PwsZOIRSIljy13z+vntcFYXY3hYydUjcloNqNlyWJMnhlEcGpa1bUsDjtqGxsxdOq06q49ADA2NoZ169ZJOjbF+vN3Lz2T2v7eW6qv82cbb0ttPfBmyuKwq77Wp771d6l//M3LZY/r7u4u+X1tW0tq64E3U6tuvlH1mM678rLU1gNvpuZfeL7qa32j67nUTyjc88ZFC1JbD7yZWnH9taqvteG730x9c+eLko4td9/v69yR+vz/+W4KQOrTn/50atu2beJ3n/vc51I/+MEPco7/n//5n9QVV1wh/vz666+nLr744tQXv/vt1NYDb6ZqWpvzzp2cnMy5xsTEROk59c2vp4aCftX3qX5+e2rrgTdTl3zietXXWnnDdamtB95MNS5aoPqeP/D2G6mbvv5l1WPyNjelth54M7X6lk+qvtbnHvpO6rF9b6q+DiCsUxu3/m/V17nkk+tSWw+8maqb1676nt/2wD+kvtH1nOoxNS1elNp64M3URWuvoXKvyo2bxYcLDRgJBKG8HWcGJLymIskOBKnQbSzCayrjstthTKm+DNVmmrSEKEBukre/vx/z5s0Tv2tvb8fAwEDO8cWO8af/fWRc2ecODw+jubkZANDc3IyRkZGSY6L1b6OtUATUK+9I7VDF0YB2O+j0SadHuYn1hCpbQAH0ZOI08/1agYuxigaDSFSpXzlpdA0miIRCMFstqpV3NPNoNNurWJ0OKi8xbZk4rQU9u36ou7sbHR0dWLhwIcxmMzZs2IDOzs6c4zs7O7Fx40YAwKWXXorp6WkMDQ1heFygYJobGvLO7ezsxO233w4AuP322/HSSy+VHFMkFEKqCqrnFC2Jv3ANOnOKVu1Q9jVo5ayoGiuKog+1zXUB4flRWQ8oOtVawcTjj0QCQYxMTKi+jtXhoCIaAHJfmFK8+X/+53+WvA4LL1it8o6IPva+uVv9mChW0ctRk5W779FQCI7WFgBAIpHA5s2b8dprr8FoNOJHP/oRDh06hE2bBFHPE088gVdeeQU33HADent7EQqFcMcddwAQcnLJeBxdL7+Scy4APPjgg3j22WfxV3/1Vzh9+jQ+85nPlB6TxDlVDrRUk4C8KKbUPac5z5OJBGKRKJViZYvdjjOnTpc/UAKioSActV7V17E4HJKFROXmeSQYgqm6GkazGYlYTPmYKM4poPy4WYCPsQoGEZpRrwCzUGhiS5AdMZRaWLZt21Z6TJQVisI1VXrB6fP/+NvfqR4TqaKnRdn4RsclHVvuvs/2gru6utDV1ZVzzBNPPJHz8+bNm/PH5HAg7A9g6cqr8r6bmJjAtddeK2m8QGZr+3JzqhxotRUD5HVmKHXPaTIIAKhtE2Jx2HH0/cMURkRai9FQTUp3ysrP80yxshphhBitU1o/y42bBbjlrKjkPBx0OHOAXi6GZnhNqwCXZh6NXIcWb05zTFSoSUpV/UC2A6RuwavEnBX1OUXLAXLYVTdGJpihJF2nPc/JNdWA9vPTApxyVnQWOxq7BBNQ4/IpeiyJWAzxmRnVNKAY7VG7V3QWFqHfHUVjxZmaLIcIJeGAuEklpd6AgHpqmbZnTmtPK5qiHWpziqaQSNyAUaUDS7Grhlagaqw+/elP4+DBg0gkErjkkkvE30cCQRiMRhw+cgQ9PT247rrrCp5frr0NE49llie1Y8cO7N27F3v37sWJEyewt0jNzIkTJ7B//3585d57c65HY1xKJ+a3v/1t9Pf3Y8ezzwIAPrT8/ILHSWlNlA1a7VWKvcQPPfQQDh8+jH379uGFF16Ax+MpeD6553v37sVffn4jzBYLDCa1YgZlc6rQPYzOiowLtXeSApLLoaF6ze720d7ejl/96lc4dOgQDh48iC9/+ct5x//Zn/0ZpqamxHfgW9/6ljAmUTRA0QGyy5vn2c+/u7s7rVC05Tw/pfccSIuuKM2pbEdxyZIl4v3cu3cvpqen8Td/8zc55xS77yIDRCGySsTiiM/MlDzuySefxPDwMA4cOCD+TmrbMbnrilxQNVYHDx7Erbfeit/9LjdXUuNyAQCu/NhHcf311+Oxxx4r2JamXHsbWhX0QHEuf8OGDVi5ciVWrlyJ559/Hi+88ELRa1x99dX4vz95CjPhCJIJGuL8NG+uYmI+/PDD2HTPlwAAf/zNb/O+l9KaaDZoeJxCkthU0DDs2rULF1xwAS666CIcOXIE999/f9HrXH311Vi5ciUe+bd/AwDZC95sKPGCi91DYlwsdnvR9k6SxuRwUJ1ThN6Kx+O49957sXz5clx22WW45557Cj773//+9+I78M///M/imAB6TllEIQ1Inv+qVavyFIpq7nn2ddTOqdk9+I4cOSLez0suuQShUAgvvvhi3nmF7rsYWalNV0jM9//4xz/G9ddfn/M7KW3HlKwrckHVWPX09ODIkfxK8gvTHr6huhonT55Eb28vVq9enXdcufY2NHYJJpDisXz2s5/FM888U/I6tHrwieMKqt+Fk0zsQly+lNZE+WNST+Nm8jD592rXrl1IpBfmt956C+3t7WWvR4/ekk8DFruH2Q5QsfZOrMZUCsTZGBoaEpmCQCCAw4cPo62tTfKYAHo04EworFpgMXtMau45QHlOFTHq11xzDY4dO5bTXUXSmNTmHO3S1qnf//73mJil3JbSdkzJuiIXXHJWtW6B1iELXn9/f8GXpKmpCUNDQiuXoaEhNDY25nxPqygRgLi1fTHv7qqrrsLw8DB6e3sLfp9KpbBz507cdOstsJroiSqV0CPZ2Lx5Mx78138FAFgK1Pu0tbWhr69P/LnYs8gZE4VkuNQE75133pmn6iMg93z37t248vLLc66rFFYFNGCxe0gMscVuU3SfCWjS3UDh57dgwQKsXLkSb7/9dt7xl19+Od577z288sorWL58OQAGakAFOavs53/XXXfl5WHU3HOAXrFyKWdjw4YNRR3gQvddrL1UGVmpmVPl1mVA/b2XAtnGateuXThw4EDe56abbip6TlW6f1S28k5uYa/JYknTSOpfll27duFPv/89AODr999f8N9w2223lYyqrrjiClxyySV48513UOfx4qqr8qXPSlCOHil1/x9//HEsXrwY//KvDwEAvnn/N/LOL9QpudyzoLGnzo9/8hMAwL98938XnTff+MY3EI/H8fTTTxe8Brnn69atww3XrQVAo9hVfmRc7B5mK7eU3OfsMdE0VrNzjg6HA88//zy+8pWvwO/Pldjv2bMHCxYswIoVK/DDH/4Qv/jFL4QxScx5yBqTzDmV/fzvueceXHrFRwBkhChq7jlAU3RV+PmZzWbcdNNN+PnPf573XbH7TksNKEdOrwRq770UyA4J1qxZI/uPDJ8ZQO2KpSLvXagdDpBpbzM0NJTX3sbqpKdmIf+G7+35HZ7Y/hReeSSX2zYajbj11ltzRCKzMTgoNLusqjZhcmwcq1evxu/TBlANosEQ6tpay469FAiXv/rii/O+k9KaqNCY1NKAf/O1r2LzU0/gS3dvwpE338n7fuPGjfjEJz6Ba665pug1yD0fHR3Fn373OziuvoyKJFvuS1zsHsYiUSAlXFPJfRbHRGl7EILsbh8mkwnPP/88nn766YJ5k2zj1dXVhcceewx1dXX0qUkF0Xr283/xxRdxwUUXCddKGxk191y4DqXSgyLlEOvWrcOePXsKtu0qdt8nJicBqG/BpsQpIyi1LhOovfdSwIUG3NX1KgDA6XVj4cKF6OjowDvv5C9YpdrbEGqMFg0BFBcOXHvttejp6cGZM2cKnme32+F0OoX/dzpR763BwYMH6YxJRRRD+HmLw45UIoFDB9/PO0ZKa6LZEFtTqVBJlcp5rF27Fvfddx9uuukmhItsT55zz+12XLZqdc51lYC025LrAJW6h6lkEha7vWh7JylgaRiefPJJHD58GA8//HDBY5uamsT/X7VqFQwGA8bHx6mKm4D03l/pPZakYPbzv+6669A3KCyG5F6puefCdTICGaUwWy0wmgozQKXYmmL3nVZRvhoaUErbMSXrilxQNVY333wz+vr6cPnll+Pll1/Gq68KRupgeluF7zzwAF599VXcc889SKapwW3btokRzIMPPog1a9bgyJEjWLNmDR588EHx2iSyoitmKBwxFOKVW1pa8PLLLwMQJtYf/vAHvPfeezhn6RIMnzmD1157jcqYlKqkAEECvn//fmy84w4YUYWvfvWreWPPbk10+PBhPPvss2J7oWKgoZIqpSb793//d7hcLuzatQt79+4VVVzF7vk777yD3/3q1+kxKX+JrQoVbsXuIWnvZLHb8corr+D48ePo7e3Ftm3b8KUvfUny9WnTgMQpu+KKK7Bx40Z8/OMfFyXS69atw6ZNm8Sxk/KT9957Dz/4wQ+wYcMGYUyUDeiMTMMw+/m//PLLOHTkAwDAzZ8UqGQ19xzILz1QgmK5WZvNhjVr1uSoi6Xcd3ItGgInKc/vZz/7Gd58800sXboUfX19uPPOO4uuy2rXFSVg3trdaDKlth54M3XNXbcrvsY5H16Z2nrgzdS5qy+hNq6/feGnqb/8/oOqrvGPv3k59alv/R21Ma29567U1gNvpqqqqhRfg9a2AuSz+pZPCttotDSrvoa3uYnKmBw13tTWA2+mrrjt04qvQXMrDvLpmxwXty5R+vmHNzpTn/3H+6mN6VPf/Hrqn377iqprbNr2g9T/95P/pDYmcRuN9jbF17j8s7ekth54M+Wqr6MyJmdtjeo5VTcvPac+uY7avbqvc0fqcw99R9U1/vmPr6Vuuf9r1MakxYcLDZiIxxGLRlW1NqKtRgLo1A/NrqlQCxq7cCpRuJUCjW7wtKXPNNoI0R4TACQTSSqUDa0WQgC9eU773QPUzSnancQjFOYUi+7mVGhAim3FtAIXYwWk+wOqCGVFgQXVl1idys1gEjh3mtQkHcPAxoCqeX60e5PFZ2aQiMVV3yeaYwKAVDKpqiamqqqKugNEozMDizwaua5SWBwOJJNJzBTJc8pFPBpFIq52TtFfp9SunSYLyaPpxkoSosGQqsiKhcBCbbcIMiYWhkHtC0PbMyfXVQqr3Y5YJEqtKwMgLAhqCjhpb5sACMZKnejDCoDRnFKVc6SdR6PglNntVFpSZUP1nGLBAKmMrKwMDKgW4BpZqQqvRek63UmgZmKKk4ABPaIqCqUdWYXo0IA0I1BAfc9CmltxECSTSVWdGVh45jQ6M7CiAdVEobQl/oB6ypTFnFK7a4XolFF8flqAn7EKBlVVYVscDiQTCcyEI/TGpNKAlmprpBRUaEDaRaVlun1IAW0aCVC/TQgLw0Ck60rBYkwzFCi3asqGgVYelDa1FQ2GVOWLWTkbVOaUTgNKQ1Ql78pkYqoNrxm03adCj1CnbOjkrCrt+bHYNiGZTKqMiumPSe3WJdU2GwwGA/VoAVDXRoiJAxQKU6GWaa8JqiIrnQaUh4jKG06ziS2B3MLE2SAvWqXRgNV2G/V6NKCyoj1APWWTeYnp5T1SyaRYbKxoTEwWO3UqN/LeRgL066xUOUAs5lRI+RY9AJs5lb21vaoxUabheYNvzkrlYkebc1VrGMi/p5LEDBkvmN6YEnGhJ9zZSANGQ2Gk0gXqNECK3ZXmrVgoFNXOqVId85UilUqpz8UwogFVsS1pp5rmnMre2l7RmCpw48UPf/jD2LdvHywWC+x2Ow4ePIjzzy+8/x4BvXbhZaA2lLU66e0STJD9EvvHJ8ocXWBMZBJQVSiqowHJPaY9MVULZBwOTPTT7RWmmgZkYEBTabWj1W5HxB+QPyZxl+DKiaxY7TIbCVJIDVSaA8SIQQCEf29walr+mBitCWqwe/dudHZ24rvf/S5sNht++tOf4v3389vDZYObsYoEhFoPo9mMRCwm+3yaW9oTqK31EAUWFCcBkXdXkhcsXE9l6QEjL1hVfojBwkI8avXPr3JqmjI0IGU1Z0Cl6IpBoSsVB4i6UVfXBooFA0QD3/nOd9Dd3Y1IJFJw1+rZ4JqzApSHshYHOxpQ7SSgbhhULMKsvGAa+SHaLwupk3tN4ZbbFocD9d4aHD58GPv27cMLL7wAj8cDQNjvKRQKiT30pO46m1RrrBg8P7WdGVgl6NWKrmiXaADpRgEVNs9p0YC0a9LUora2Fk6nEy6XC9Z0fWEpcM1ZAcrVPywEFmpfYqvDIURCcXqFroA67070olgYK4X3iUVXBjImo8mEX/3614q23LY47JgcG8MFF1yAiy66CEeOHMH9998vnnfs2DFxm/EvfvGLksaUSqYA0DAM9BYWtZ0ZWM0pNaIrIoxisSaYzGZVYgZ2kZXS5+fATDhCtSCfBv7zP/8T3/rWt/D000/je9/7XtnjOXawIN6B8oiBfmRFIdpjoLBRkx9iIX0G1HWDJ3UrLJRbAPDsC88DkL/ltsVux1D/GSTSL/Fbb72F9vZ2VWNSHVkxEH0A6grgxe70LGhApXOKQW4PUF9AbXU4qN8ncj01QrBKk61//vOfRzwexzPPPIMHH3wQq1atwtVXX13yHO6RlVJPyuKk75mr9VhoN4wlUNOaqhIpG1ZFieR600FByCB3y23BC84sLHfeeSe6urrEnxctWoQ9e/bgN7/5Da688sqi47jrrrvQ3d2N7u5ueNM0oprIisXCoobGZTanVORBWSncCFWmXM3JRvRBrq0ELKI9tfjJT36CT33qUwAEB++yyy7Dr3/965LncBRYEJWb/MnJoigRyCrAVbyw0DeggEoakEEPRYDWYqdsTLt27RI3lszGw09vF65f5vkV23I7O7/wjW98A/F4HE8//TQAYVfa+fPnY2JiAhdffDF+8Ytf4Pzzz8/bBh4Q9mTbtm0bAODdvXsAkwoHiIHoA1BH41ocdurdYwB1DVpZOUCZAmoVDBADuhRQnkJhoVDUAvyk6+lJoMSTEjlzikWJQFZhooooJkJZXAEIL7Grvk7RuSxfYq1EA2vWrCn4+yWXr8L5f/EptM6bh8EjvbK33HZfegGiwRA2btyIT3ziE7jmmmvEY2ZmZjAxIZQz7NmzB8eOHcOSJUvw7rvvlhxrKpFUZayqGTRnBdQ5GyzyjQAlB4hRtK60ZyGT3CxZp1Q4GyzWKd7gTwMq8FhE6SzlyCqVSqmKYgR+urIiK6uT7rYJ4pjSCsVCUUo5sEvQC9e74aZPAJC35fYvf/lLWOw2zGtpxX333YebbroJ4ax7Vl9fD4NBeD0WLVqEjo4OHD9+vOyYUqmUUHpQYQuLWtEOC2oyEgjCYDSi2lZeCZY/pnQerYIoN7JlEPXaPbK1vQqn+myIrPgLLJzyJ4HS7celQC09wqKFiaqcFSvPPO1sKGnyyUr0Qa530cUXy95y+9ipkwCAT95wA1wuF3bt2pUjUf/oRz+K/fv347333sNzzz2HL3zhC5icnJQ2rlBY5ZxiYBjUzHNGG/dlCuAVsC2sSjRU7DDAioIHVEbGOg0oD0Q6qcQ7yHjmlTUJBDVghS0srBY7wuUrKCFglqBP/zsf+rd/wzsv/jLnu8HBQdx4443iz11dXTniCU9TAwDgn771D3jrufxo7IUXXsALL7ygbFxqImO7HWPBvvIHyoSqxc5JX+EGZKncnA74x8ZlncuyRCP7+nLAyikD1Of3zgZjxS2yApS3V8m0EGKQH1JLAzKJ9oIwmkwwWeQ32GWpJgOUSXpZbFIJqOvMwMozJ9esNOWWIF2vLOmzqtQAw4J8QJnoimV380qkcXmDr7FSWFeR2TyschYWwrWzogEBZZQpq9ovKi9xBXnB4o6uzKjlypE+AyrpbsY0oBLKm0jLaVPeaspZWHY3V7p2VlVVnTVqQK7GSmkbIZaRVTQYglXBdt+saIjsayrZhpwVP62mwS45h7boI5lIIBaJKlJJsdzqOxpSZhjEhYWRGtBoNsFUXS37XFa5WTXlLKzk9IlYDPFYTJlTxjJaVxgZE6M+13cJBjSIrJRJ1xl6wQrD68zEZGFAlXfWYCf6UO4Fk2gvlUrRHpbire1ZesFK5xSrTh+AOuGAlUFBPqByTjEy6oAQramZU5WUW2elmtQC3HNWSr2oRCyOeDRKf0xKqUmGkZVaeovpmJREMQxpCPXORuVQy6zo0uxrKr1XTBiEgDqFMKsFOBJU1u6M1fY8gHKBBUvRB2/wpQEVRlYsmtgSVCI1qWZLAJYtoADlNCAzY1WhhkFdgp7hnJI5LpPFAqPJxDQ3q9SBrbg5xZgGrLTcLG/wpwEV1lSw4lyjoZCibcjZeuYqaUAWyi2VOSuWzkblKbeU1VmxjvYA+XOKRD0sxpSIx4Wco1IakGG0rvTdAxjlsQNBYT9Ak7xqI5a5Wd7gTwMqCfkZ7BJMoDRiYNVVQ7imsjEZTUICncVLTDaFVOpsVBwNyChBDwhzotpmle8AsaSWFcr8xUJXRs9PaTkL88hKoWgnHosp2lxWypgABc+PYW6WN/jTgA4Hqgzy/iyreiZAeTNbtjSSsvoTlmMi161EGlCxAWWUoI8qVLkRL5hVb0A1Y2K12CkXXdmZ0KWAOjEDi+JpQLkYhWW0zhucIyvCm8urQWFVOwSo8FjE2i/64yKLlXIvqvKimEqjAVkn6AH5wgHS3SUcCFAfk1KBjIWhaIBcV1m7JRuzaE/NPGcXgaqNrHRjJQtKK9ZZtXsBsrcEUEYDspgEShvsislUVouwmpwjq5dYjQFlNSa/YGysTqes82zp4yMMjJXixc7OjpoUrqswsmLUAxNQRwOyc8pUsi16zkoelO7LwppGAuQnni0OO+IzM0jE4yyGpWi3YFY7uhJEAkFYXfIWYIA1jRtStFEeyzwaWdjlLsLEuEX8LGq/lAlkWNOASjf1ZLomKBbIsJvnaiIrVrlZ3uCeswIU8K48aEAFhoFlVbiSzuusQ/6w3w+bTGNlMLLZNoEgGgzBYDDI7gbPJbKSzSAIDlB8Zob6mJLxBGJR+d0+LE72Agu5oiuDyQizxcKMQVDa7YOpAVWYB2WZm+UN7u2WAHkvcVVVFVPPXOk2BaybQyqjAdOUDSsa0B+QTW3xEH1k/x2pYJmgF3NWMg271elk7gDJjqwYdmoBlDEIrBojEyiNQnk41bJLDxjmZnmj4mlA1r2tlDaNtTjY1X4BCmlAxgtLWIFyS5T4++nnYQA19BabTh+AisjKxdhYqXCA2OWH5M8pm4vsG8VmToV9ynKOTJ1qMbdeOSkU3qh4gQXLeiYgmwaUK+llNzGBdGFihdGAEb8fNrdL1jnkhWehcANUeJwudjRuJrKSb9hZ3SdAWWRlSdc4sujrCAg7KZgtFhjNZuljIkIiRs8vnHY25FLeFZlbP0s6rgMaRVZyFmGWuwQDEEQSsbhCGomhsVKysDDPWQVhMptl7bNFXnh2hkFFzpFRtBeLRJGIx+VHVgzHBJCaJpnRAuPFLqKgWwvrORX2+4W/I9Mxq7bbmK0Jwtb2YT2y4gUl3gHLTsYESugRHgILJYXKLBWK4kssw+MUvWBWNKCCYtdqmw0Go5EZjQQoU05aXU6mnQaUlB6wXuzEbULkOLCiapIRDaggsqq22WAwGJg+v6gCMQrL3CxvcDVWgncQkjkx2UZWgLIuy8xfYgVbX7A2oOTacl5iMb/ALEGfzi/IoNyIEQkzvldK6gnDDGTrBJFAQDY1aWHY6gxQVj9E/g2sKNOIT4lTRpoEsGVbZKs5GdY48gZXYwXI35KDNT8NKKfc2L7EIdmNK1kn6DPFrpXoBUunbGyMRR9A2jDI9IJtTifTaC/sD4iFx1LBnAZUUM7Cup5QyZwSnWqGUYwg85ffLkunARVCbi831kWJgPwuy1UGA/PEpRJJNusEvbKXmG0UYzUJifn7v/VN7Ny5E16vt+Bxa9euRU9PD44ePYq//uIXAAgG5dvf/jb6+/uxd+9e7N27F+vWrRPP2bJlC44ePYqenh5cd911ssYldGaQmaB3so+MlZQesKYmhb8j3wFiNaeioZDQtNktI7Li0INPqUBGl64rhNzGlRbGAgtybTk0IOltWGnGyuZ0Mo8WAHn1Q1anA/FYjMnGmQDwd3/7t0jF4/iP//ovvPHGG9iyZUveMQaDAY8++ijWrVuH5cuX4+Nr1gDIdIp4+OGHsXLlSqxcuRJdXV0AgGXLlmHDhg04//zzcf311+Oxxx6DQUYD5ohfnhdssdthMBgYG6sAjGYTzFbpAhlWm3mKY1IgurK52M4pIB2FyonWeVDLMjvUG81mmMzms2JLe0ALYyXzhrPcNoFArsdCxs+SBlSknGRMAxKBhTwakF1fRwBYv349AlPTsLmc2L59O26++ea8Y1avXo3e3l6cOHECsVgMv3vzjwBK5zzWr1+PHTt2YGZmBidPnkRvby9Wr14teVzRYFAW5WZlXDsEZCJjOdEVaxopqqicha1TBsjv1iLmQdP5LhaQv06dPU1sAa1yVjIXO0FeTn+PGHFMQXnKLR6djJXUf1mdDrYKNz8RWMjzOFlSk01NTQhO+2BzuzA0NITGxsa8Y9ra2tDX1yf+PJUeD7lXmzdvxr59+/Dkk0+KNOLsc/r7+9HW1lZwDHfddRe6u7vR3d2N+vr69LXl5WZZ5/YAZQIZ5jSggg71VpeTqaMIpCMrGdJ1lk2ICeQ2CjibmtgCmuSs5EVWrBVugDAx5VGTHOT0inJWTtF7ZgGRy5ejvHM6VTdm3bVrFw4cOJD3uemmmwAIC3ypBbiqqirnZ6NF6PkWCQTx+OOPY/HixVixYgUGBwexdevWgucAKFoYu23bNqxatQqrVq3C2NhY+tqlxzQbZP4xpZFkCmSMZjNM1dVMWQ2yqacsB9bhYNLsNxsRnzwaUIysGL9/StapsyWykrdHMgXIz1mxVd0BglTV6nDAYDQimUiUPZ6XnB6QbqyqqqqYt4AChOcn5yW2OO2qvc016RxTIQwPDyM5MwOry4nm5maMjIzkHdPf34958+aJP3vr65FKpTATCmMk6xlu27YN//M//1PwnPb2dgwMDEgecyQQgqm6GkazWRIrQCIrlpSpmHOU+P7xopHkpgZYMwiAQAM21i+QfDxxTFgrJ4lCWEotpShEYUyZ8oImOSu5+SHWLwt5mFIVSbyUP4B0Lt/iSCfomXP58qJQm5MtZdPZ2YnWhkbYnE7cfvvteOmll/KO6e7uRkdHBxYuXAiz2Yyl5y/HTCiMVCqF5uZm8bhbbrkFBw8eFK+7YcMGVFdXY+HChejo6MA777wjeVxyDYNNjKwY5qxEmbi0iI/XluhyHVghN8t+nstVvUYCQaSSSWZjksu2kPGzzKPxBPfIKhoIwmgSFEmxSHk1D5fIKj3xbW4nwj5f2eNZ9ysUri1vYlo5cOZAmnKTk6CnQAOWwoMPPoh/6/oF/PZqLLPV4TOf+QwAoKWlBf/1X/+FG2+8EYlEAps3b8Zrr70Go9GI/z60ByGLMPUfeughrFixAqlUCidPnsSmTZsAAIcOHcKzzz6LQ4cOIR6P45577kFSxkIUyTIMwcmpsseLW3Gw7A0o21ix3R6EQO5uwUJkxdqB9ctTvbocogCJFaJZG8WGpsuvUzYO1CRPaEADZiIGScbK6YB/dJzpmOS2V+EisJDZTZxHVwbh+gFZ9SesKZuJiQn899M/w8f+8s9x7bXXir8fHBzEjTfeKP7c1dUlytLveORB1LYLYomNGzcWvfYDDzyABx54QNG4MpGVRC+Y4caLBGGZ3T640YByIysuNGAAFrtNMuVmc7m4UPCAjGjdTebU2RFZ8RdYkEVYMm/uYE5DkDBZssdpZ99VIxlPIBaJSqYBeXRlANLCARmRlcVhZ0ptAcKYTGaz5PohHtLniMwoxupyIplMYibMbqM8YnSkPj8eQiJyfTkqN6uDbfE7kLUmSHRgrU4Hc7pNrkLYSmhAxveKFzSQrsvPxbCmIeR2ZrA67YjHYkzl9IC8/B4vGjDsk74BY7VN8ExZK7fkPj/WnSIA+XVypB6N1VYcgNCbU04XC9Y7HhBEZURWFrtdaELMfE7J6w/IJbKSKfO3uZyIhsJIxsuLxuYCNBFYAHIUSewFFtk5KymwcBgTIK8NFPEA2VMR0iXZ4pg4qDkBGQsL4x58QIbOk2wYGLfKIpDTzJaX9FlO018r48bIBGQDRskOrIs9NRmR6ZTZXC7mDAJP8KcBZfCuBqMR1TYrlwJAYUzSFxbWExOQ1wZKrNPhoAa0OB0F65DyxkRoJA5jAuRExg7m90l+ZMW2+wiBoOaUJ7BgXVQakbH1hdg9hvXzE/e0ku4AsZ5TobT4S2qxss3lZC764AlNOlgA0mTimbb7fBKXkj1ztwshDnJQOV2WeSl/Iv4ADAaDJHoyE+2xr4kR/p5U75ztvlGA/F2xWbelIpAjZuAlsIimI6sqCb0X+c0pIkaRGlmxN1ZitCfZWLnOGiUgUOE0IC/OPJVMyqqrsLvdXGoX5PQCszqdzJt7Atkqt/KGnYfCDZAXWZksFpjMZqb7RgFAIhZDLBqVZUC50YBS80MOB2bCEUmF8qrGJAoHbGWPtXKbU9Kp5WqbNZ2bZbsmJGIxzIQjsLvdko7nIafnCc0iKykeJy81EpDOxUgN+d0ubsZK6mZrVifbLdEJMh5n+XvFLWclY0w2DvVMBLJyMRzaigHy6uSEXWbZ52blrAk2DjWOgLycFestS7IR9vllRVZ6zkoFErEY4jMz0mgkTpEVII/Lt7mcXIyVXBqQy2JHxChSjFUF5qx40UjC35DeIJmbsyFnTBzadwFZrcUkzHULpxKNmXAYiXhc0jwnx/B4fiGfT2bOSjdWqiCVNycTkweXL2dLAG6RVSAk2Qu2OvnQSMTjlGLYMx4n23HFZ2YQi0YlRcbEAWJNAwLylJM2Dp3EAZkycV6qVxmiKxunOQWko1AJhoFHE1uCsM8Pu0RjZeXkVPOCNsZKYv0QVxpQYpdlk8UCs8XCRWARmvbBbLXAZClf7MrPM5cu8xcLXUPsCl0JpOYceUZWoWlpXrCpulrobs7BgIYDAaEZqtlc9liLk32rMyArjy2BBuQ/p2REVlycRWk0IKlx1CMrlYgGpLW650oDBgKS6BHi1fDwWEj/LymeFOuNFwlEGlBSZMW+0JUgIrHBrpVjzio0NS0pGc5zTHKUr3a3W1IPOrUQmzZLNAy85pTU/oAig8BjTZBorIgzqQssVELqlgA8GsYSlNsTicAmGiv2LzGpq7B7pC14XBY7GcWudo9b/DewhtTIipdCERAWFqnPDuBDbWVk/hKf39Q06yGJLIUkA+pxIzjNfkyAkmidj8BCigMkdlzXIyt1kLpbMM/Nw8jWF+WKXckk4EUDAtKMlWAY2I9JzA9JXVg4LHaAUMQpy9ng4HGGpn2we9xl5xQvOTYZEyBnTnFwymSMyebhE+0Bacqt4mhAH6xOYe89SWPSIyt1kLpbsJVTnQcgTEyD0YjqMrUemcWOR85jOudvFoPRZILV4eDiBQNpwy7hJXZ4PAjzWlgkJsPtXg8SsTgXLzg07YPBaCzrmIkLCwcGgcyRcoah2mYVcrMcnl88GkUsEpUUMdjdbm5zKiI1snI6kYjFMROOMB+T1ChUbGLr0yMrVZCsBuSwlxWB1FyMnScNKNHjJN/z8jil5ofsPL1gqQbU6+FGTUp+fl6PcDwPyk3qmNx85xSJQsuhEueUjcNmkAQkL1bOMdNzVpQgOWflsHOhAIGsYteyk4CfwIL8jXIep2iseFFuEre250VNAtJLD3hSkxmBjLTnJ2WTRrUgY3J4PCWPs3F2gEI+n/g3S8Hh9SDIzVj5xT2tSoFnPZMYWZWZU3rOihIigSDMVkvZSWBx8qnzAKRvwMiTBowEgkjE49I9c25iBn/ZNkJVVVVCD0Vu+QVBkm2qri55nMPr4WbUpUYxjhovAHARDoT9ASSTyYqL1qVEVlUGA6wuJ9fICigvRrG5XdwimIwDW9pZtHIqnuYJjXJW0rZs56omE7ng8pFVJBDkkkcD0uofiQsLNzGDhMjKkk4Cc6MmJe6Ca/d6xFwga5C/U9ZYeT2CEeGw71AqmUTEH5BhrPgZ9nJjsjqdMBgMHI2VtPyQ3cPPAQpL7Lxuc7kQDYUl7XI8V6CZwAIoX7Hu4DgJJC92bhc3AwpIe4kJpcNTJSXl2QEcxySx5ZLD40Fwih+1BZQ3Vs4aL4JTUxxGJEDKnKrEyIr3mKT2B+RJTYak5qw45tF4QTMaECi/TYjd6+EWLUjND/FqtUQgpRdYJmfFK4opv9ss72hPshfs5VM7BADh6fScKrsIexCc5DMmgBiG0jkrTYyV1NwsNyGRtD2t7F6+cnpA4jp1FlGAgEbGitzEUrwr75xHyOdHMpmEw1sm8cx5jxhJHmdajs2jQzYgFK9a7DYYTMVrPciYucmMJXjBRI7Nq6g0EY8jEgxKyFnxoyYBITcmJYqJpSXlPBCa9qHaZi3ZWow3NSllTyuD0Sh0+uDkAMVnZhCLlK9ztLmcZ1W+CtDIWBHKw17CMFhdLhgMBm6eeSqZRHjaJya7i4F7ZCWRHuG1AANZHmeJ6EpcWBhTpjU1Ndi5cyee/+//BgDUNTcVPG7t2rV4e88eAMBHL/+I+PsdO3Zg79692Lt3L06cOIG9e/cCABYsWIBQKCR+9/jjjysanyQa1+vlGlmFJVLLvBxFIIsyLeHAOrycaUAJ0TphPXitU4BUtoXv8+MBbYxV+sUsZRjEiclxEgSnpstGVnbOxkpKexWetSdApmN5KY+TF2WzZcsWvPHGG7j6qo8CANZ/6lN5xxgMBjz66KO4e/OXAABXrLoUy5YtAwBs2LABK1euxMqVK/H888/jhRdeEM87duyY+N0Xv/hFReMLT/vLU25eN/+clbcMjcR5TklRTlZizop3tAdIa2brrPUiMDHJaUR8oFFklTZWJQwDibp4RgzBqemS0R5AtrTn+xLb3K6S7VXsHn5V/WRM5O8Wg41TZLV+/Xps375dTCZfsOKivGNWr16N3t5eTKWFPb//1a+wfv36vOM++9nP4plnnqE6vnKUm6m6GlaHg3/Oyu0uuY08dwdIirFKO228nEUpe1oRIREv0Q4g0JPlHFhnbQ0Ck7qxUo1ELIZIMFjSWIlqMq6R1VTJMRlMRljsdu40IFCaihDUSBzvU/olcJaKjD0eQeLPWI7d1NSEoaEhzIQjiMdiMBdol9XW1oa+vj440gvh4Ok+tLW15Rxz1VVXYXh4GL29veLvFi1ahD179uA3v/kNrrzyyqJjuOuuu9Dd3Y3u7m7U19fnfFeOBtTKKQNKzynexkqqAxT2B7iVjQDloxie3UcIyo3J6nTAVF2NwPjZZaxKV+UyRHByGg6vt+j34kvM0WMJTk6jffl5Rb/n7dkBmcjE5nYV5cXtHjf6D3/AbUyBdKcFZ6236DE0O0Xs2rULzc3Neb//+7//+5yfg5NTiHjq844jjWTJnIqHI3lbTNx22205UdXg4CDmz5+PiYkJXHzxxfjFL36B888/H/4CxZ/btm3Dtm3bAADd3d0535UzVsQ54tG9gkDMD3mL5zXsHjf6Dh7mNyYJ3T4EA8rPKADle06SdAVPZyPk86Fp8aKi3ztqagBk3tOzBdoZq6kpOGpKRFbEY+E5CaamS7ahIWPiOQmkcfkebrJ1AAhOTAHIvBSFQLOge82aNUW/Gx4eRnNzM4aGhhD1BzAdzldE9vf3Y968eaKxaq6rw8DAgPi90WjErbfeiksuuUT83czMDCYmJgAAe/bswbFjx7BkyRK8++67ssZeTpItGiueCXqJlFulRVa8oz2gfBsvuwYMULndgl21aWOl56zoIDTlKx1ZedxIJhJc5ZfBqSmYrRZU26wFv3eSSTA+wW1M5V5ik0UYL8+XOBoKIRaNlqQBeeXROjs7cfvttwMAvHYHhifG847p7u5GR0cHWubNQzQYwobP/i90dnaK31977bXo6enBmTNnxN/V19fDkM7pLFq0CB0dHTh+/Ljs8YWnfTCaTbDYC3drISIjnotduc7rmTnFb0zRUAiJWLxkf0DeCkWgfOd1nl38CQgNWCzn6EgzHoEJfusUD2hmrIJTUyUVScSL4rEjqDgmolIsYkSddbUA+EZWYhFgkZdYLL7lTI8EJ6fEl6IQeHnBDz74INasWYMjR46g0VuLyaDg3LS0tODll18GACQSCWzevBk33rIeXqcTzz77LA4dOiReY8OGDXnCio9+9KPYv38/3nvvPTz33HP4whe+gEkFCetyzoYmNGCZMRGvnbdhCPl8pUVXmkRW5WhAfl38Ccj6U+xeEaeaMCBnCyo2Z+UowaezAlnw7V43JgeH8r53ahBel11YOMt5CQITU3CWowE5jGliYgLXXnstAGD9330Fq2/5BAAh53TjjTeKx3V1daHtk9fCWVeL7z/wQM417rjjjrzrvvDCCzkydqXIbrlUaE7xbGIrjqlM53VRNMDbWJUTo2hgrKSMiSeFCwD+MYE9cNXXFVyLnGdpzkrTyMrmchbtvM6z1ZI4pnKRVW0NkokEX0lvmfYqROHGk0YCBEVgsTq5qqoqTRaWwOSkqIQqBJ6dBgikOBu8mtgSlOu8rpUDVMowiB1tOEcx/vEJOGq8RUtHeHbxJ/CNCsbK3ZAvJgKEdSoSCCI+M8NzWMyhobEqzZtrwU+Twsxii7CztgbBqWmkkkluY0omEgI9Uls4itEsspqcKqoGtDjsXDuui2OaKC2pt2uwsIilB0XmOe8mtkD5zuu8N14kCExMiuzFbIhd/DkKiQAhP20wGIpSbjy7+BP4xsYAAO6GuoLfO+tqzjpxBVABxqqYYRCaQ3JeWMRi5SILS602k8A3Og53feGJqUWdByAYq2LPjqgE+UfGU8LfL2JEeXbHJihHuTm8fJvYEgSnpos/P85tjQh8o2PFowVCl/KOYtKUG8lXzwbPLv4EYmRVX/xe6caKIoJlkoSuulqRm+WFTDNbb+Ex1dZoUmjnHx0v+hKTsXIXWExMweooTLkRj883OsZ1TOTZOGvzFxZx4z6NDKirrkhkXOPl7pQBacPQWHhOudKOkZ+j6hUQDIPD64HRbM77jsx//xjvOSXcA1cxY1Xj4R4Zx6NRhH3+4pHVWdi9AqjQyMrmdsNUXS16ELxQrpmtVpNgenS0qLHyNDUg5PNx645NECjRxYKM1cfZ2Sg1Jrubb2NkgkQ8Dv/4BNxNDQW/d3DeHoRgemQUnsbCY3I31CMwMYlELMZ1TP60c1OIRSBzanqEr7EiBttVYEyOGi9M1dWYHh7lOiZAeLcKjQkQxnW2da8ANDRWpAagkMdCPAbekRVQmh7RigYUIqvCE9PdUA8f5xcYKE25kcXGzzuySj+bQmOya1BkTjA9XNwwaOGZA6XH5Gms5x4VAxl6y1VgrosOEOdx+UtEVuT+TY9oYKxGx4qmBpw1NWedEhDQMrKamEIiHi8YMZCH4ONMQwDFPU6jyQSb26WJsZoeHYOpuhq2AopAT2ODJi8LMVaF5OvuhjrEYzHuUUwkEEQ8FhMr+LORKTuY4jomoPicMlkssNjtmkVWZouloMjC3dCAaU2MlfA3PYXWhIZ6xNL0F09EgyHEItGKM1b+8YmCRt3mdsFoNuk5K5pIpVLwj40X5M1d9cLE0CKymhoagacAZePQsIUJiVA8Be6VVl5woERk5aqv0+TZAYITVKgNlJYLSzFjpUVfOQJyHzxNjXnfuRu1idbJPC7kwGo1zwHBMDgL5Bw9zcK982lBA46OFRRYaFELyguaGStA4J8LelH12oT8ADA9PAJPQ0NeKxMt+21NjxaWqlZVVcFVX6cJZ04ilIKRVX0d93wjQTH5sze9KE8Pj/AeEqZHRuGsrckTDpC+cjy7V4hjSt+H2Y5ZlcEAV12tZg5QMpEoGDG46us0MaCA4DQXi6ySySR84/znun90HNU2K6xOR87vRdWkLrCgC0GRlO9xuurrMBOOIBrks017NqaGR2A0m/IWPC1ppEwRYO69ctbWwGgyabKwhH0+zIQjBSMGV0M9/Bq8wICQCy1krNxNDYiGQlx7uBEQz3t2ZEzobk0cIHFMheeUFhFoKpkUxCgFIgZPYwN3wQ6Bf2KiqLEKjE9wLegmyNRa5d4rMvf9usCCLorRI+6GOs0WO/ISe2fRIy5xYeGfR/MXKQIkFKoWCwsATA4OoaY1f+sOLSOrqeFReJvzqS1vU6MmESiQRbnNmuu1ba0AgIkzA3nnsIZvbAzJZDJvTMSgakW5+cYKi4ncDfWazXP/+ERB5Z1W+WIgS4wyy4iSsg0tonXW0DayGhmD3eOGyWLJ+b2rTrvFLkOP5C54tW0tSCaTmBwc5j6mmXAEYX8gz4vyNKapLY3okYmBQdS2teT8zmgSolLeSsDsMXkaG/LqvzyNDZobq9ksQl17K+IzM5rQW8l4AoGJyTxjRaJ3zSi30fG8yKraZoPV6dDMgAbGJ+HwevJSA54mLY1V4ciK5JB1Y0UZvtG0xznrhrvq+RcEE0yljZW3Ofclrp/fjqmhYe61JwSFqvs1j6zODKK2NddYEU9PK8pm8swgAMDb0pTze3ejdp550ciqvRUTZwa57iyQjemR0byclTinNDIM06OjeaKrTJG5RjTg2DgMRmNe/Z6WDpAoqZ/NttTXITg1jUQ8rsWwmEJzGhBA3uT0NDVq5kUFJ6cQn5mBtyl3saub14bx02eKnMUeheoqPA31SCaTXPfXysbEwCAcNd6cvZoIXaLVwkI6m2cb0aqqKngaGkRHhDfCPn/B/B4xVlrBV6DWSus55R8dz2scSyJSrdaEiQFhTtVksQimtOxfK2MV9vkRi0bzotD6+e0Y79NunWIJzdWAQG5k5ajxwuZyYux0v1bDwtTwiChLJaif146xPu3GVEjmLyZ4E/wTvEAmisnOW2nVaolgosCYBCWeCT6NIiuA5Gdzn19dWyvG+7VbWKaGRwrQgPWazinf6DgMBkNOLz6x7lKj5zeRfkZ17W3i77QshSDwFWgW0LBgPkZPndZoRGyhuRoQQE4rmvr57QCgqbGaHs6lR6wuJ5y1NdqOaSS/rkJLagsQIisAqMmKYurnzwMATA5oEzH4RseQiMVzIiuSf9QqsgKERc09a07ZPW5NI6vpkVGhZVBWzljrOSWq3LJYBOKkaUUtk3le194q/o6sD1reK/9Ybn7PVF0Nb0sTxk71aTYmltDUWIV9foR8PjQsmC/+rm6e4L1o6XFOD4/kqAHryZg0jKx8o2MwWy05jX89jQ2aJcKBzEucLbJoWrwI/vEJ7t0rCJKJBKaGh3PGJC4sGlE2gBAVZEcxdWkloJbznEQq2cyGp0HbOTU1JAiYsp+ft7lJEBlx7l5BEItEMT0yWjiy0tAB8o2OiQ0UAGHtNBgMGNXQqWYJTY0VAAwePYaWc88Rf66f145kMonxfv5yXgJCj1RVVQEA6ualoz0NjdXwsRMAgJYl54q/q2lp1jRaCIxPIhaJ5kQxzYsXYaj3uGZjAgQqsKYlQwPWp5+fVtEekN+Lr7ZdO9k6QaaLhTCuqqoq1M1rK7ijMS8MHz+FZCKB1qUd4u+aFy/C8ImTmo0JACb6B8RnBlQGDTg1PIKa1mZRpdiwQGA1Rk/qNCATDB09juZsYzW/HVOD2qnuAKHlkqm6Gq60x0kWOy0Tl/2HegAA884/D4BgqGxuFwaPHNNsTIAQXWXnh5oWLxINq1aYHBzKSYa3LVuCqaFhzaI9AJgcGobZYhEVnZnISkNjRQqD0yxCbXsbrE4HzvQc1WxM8WgUIydOoTXLKWtZci4GP+jVbEwAMNZ3JocGrGtvRdjn16RxAUH/oQ9gsdvRuFBgpoixGjut04BMMHj0GGxuF7zNgvpOayEDAPS9fxgAsPCiCwAIBtQ3OoaZcESzMQWnpjHeP4D285cBAFrPEzzPMx8c0WxMQG6tlbepEVanA0NaG6szg3A31IvtjVqXdmi6AAOZObXgwvMBCJFVyOdDxB/QbEzj/QOIz8ygfdlSAIJRB4ABjefUwJFetCwVjJWzrgbu+joMHNHWWE30n4GnqVGcU+3nn4f+wx9oOqa+g4cAAPM/tByAkC/2j09o0qWFBzQ3VkNHhcigpWMxAMEwaClkAIQoJhoKYfGHVwIQuGCtxwQI4yKRVdvSDiSTSfH+aYXsWqumdISsdWQ13j8Ag8GAhoXzYbJY0LhoAc70aLsAnzl8BPGZGSy46EMABGOlZVQFAPGZGZzcdxCLV18MQJhTiVgcg0e1pXEHPjiK2tYW2NxuMcIa1NhYkTlV194Ko9mM1qUdogOiFUZPnkbYH8C8C9LGasG8s1ZcAVSAsRpM5zeaO86Bo8YLR40X4xobhmQ8gZPvHcA5aWNVCdEeIHjnde1twkt83hKMnerTNNoDgKFjJ+Co8aJuXjuaFy8SfscxZ1VTU4OdO3fiyJEj2LlzJ7xeL47t3gsAWHL5KrR0LIbRZMLN16zB8PAwDhw4UPZ8gi1btuDo0aPo6enBddddp2qciVgM/Yc+EKP1urZWTGhsrADg2Dvvou28JbC5XWhbtgTDx09oSsEDwECa8mtdshgtHZVhrES2ZcWFaF1yLkxmM/oOamusUqkU+t4/LEZWDQvmYfQspQCBCjBWEX8Ak4NDaOlYjAuuvgoAcPTt3RqPCji2ey9al5yL+vnt8DQ1VEShXf/7JG+1FK1LOzDwgbbUFgAc/v2fAADLP/oRUQkYmvZx+/tbtmzBG2+8gSVLluCNN97Ali1bMDU0jKFjJ3DeFZehLU2X/uKpn+L666+XdD4ALFu2DBs2bMD555+P66+/Ho899hgMBnWvy8l9B9B+/nlwN9SjfsE8zaNiAOjt3gODwYBzLllREXQpAAyko+DWpR1oWXIupodHNc03AsDIiVPwjY6h49JLMO8CgYrX2lgBwOkDh9C6pAOOGi88jQ0YPakbK6YY6j2OtvOW4KK112DsdL/mlA0AHE9753f+8F+RTCTw/m9+r/GIIHLkH7rmY6hrb62IhWWifwBDx07ggo9/FPM/tJy7EnD9+vXYvn07AGD79u24+eabAQAf/PEtnHPJCixccSHCPj9e+0UnJgo0IS52/vr167Fjxw7MzMzg5MmT6O3txerVq1WN9dS+gzBbLLj27r+EwWDA+7/9g6rr0cCp/e8jFoli5fXXwtPYUBHvnn98Av7xCbSe14HWpedi4Ki2URVBb/cenLtaMFb+8QlNVZMEfQcPwWg24dq7/hIAcPrA+9oOiCEqwli9/5s/oPncc7D0I5di385faT0cAMDpg4cRi0TRdM5C/P7pZzVX3QFCXdrY6X585H/dilg0ip4/vKn1kAAAh3/7R5y7+hK0dCzG7s4urn+7qakJQ0PCojE0NITGdHPfnj+8BbPFgotvvA6n9hd/gYud39bWhr6+jJfa39+Ptra2gte466670N3dje7ubtQX2N6C4OS+gwCAyz97C6aGR3DmsPaGIRGL4cR7+3HR9dcCAE7tP6jxiAT0vX8YF3z8o2g7b4noOGqN3rd3w91Qj5Xr1lREVAUIkRUAfPTz/wtTwyMiBX42oiKM1ZvPvoj3fy1ELvtee0Pj0QhIxGI4tnsvJgeH8Oq/b9N6OCL+8MxzePd/XsVDN/95RdCAAHDgV78FAPz2qWewu/MV6tfftWsXDhw4kPe56aabip5z/N33EPL5MHryNH7+j/8i+2+SGrtsFGs4u23bNqxatQqrVq3C2FjxglrfyCiOvPkODAYD9r6yS/aYWOGDP7wFg8GAX/6fH+J0CcPOE10/eAJmiwVTwyP4w89+rvVwAACHfvtH9L1/GO+9+jp+8b3vaz0cAEJhMKHid7/0ClLJpMYjYotUJXyqbbbUkstXaz6O7I/d404562o0H8dc+NS2t6aqqqq4/92enp5Uc3NzCkCqubk51dPTI37nbWpMma0W8ecFCxakDhw4IOn8LVu2pLZs2SIe9+qrr6Yuu+yysuPp7u4ue4zJYil7DM+PwWhM1S+Yp/k4Zn/mX3h+qumchZqPYy58GhbOTxnNZs3HwfJTEZEVAMyEwzjy5jtaDyMHoWkfAmfhjpssMNE/oMlWF52dnbj99tsBALfffjteeukl8bup4RHEIlFF53d2dmLDhg2orq7GwoUL0dHRgXfeoTM/49HSY+KNZCJRkZLn0/vfx/Dxk1oPY05g9ORpzVWcPKC5xdQ/+kfpp7a2NvX666+njhw5knr99ddTNTVCJNzS0pJ6+eWXxeN+9rOfpQYGBlIzMzOpvr6+1J133lnyfACpb3zjG6ne3t5UT09P6vrrr5c0HimRlf7RP/pH/qcq/T86dOiggO7ubqxatUrrYejQcdahYmhAHTp06NChoxh0Y6VDhw4dOioeurHSoUOHDh0VD91Y6dChQ4eOiodurHTo0KFDR8VDN1Y6dOjQoaPioRsrHTp06NBR8dCNlQ4dOnToqHjoxkqHDooo1chWhw4dyqF3sNChQ4cOHRUPPbLSoUOHDh0VD91Y6dChQ4eOiodurHTo0KFDR8VDN1Y6dOjQoaPioRsrHTp06NBR8fj/ActsS0iQqP0/AAAAAElFTkSuQmCC", "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", " \n", " 2022-04-01T21:46:20.461457\n", " image/svg+xml\n", " \n", " \n", " Matplotlib v3.4.3, https://matplotlib.org/\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "\n" ], "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 115, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from sympy import *\n", "import math\n", "\n", "init_printing(use_unicode=False, wrap_line=False)\n", "x = Symbol('x')\n", "\n", "plot(integrate(integrate(sin(pi*x),x),x))\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "interpreter": { "hash": "a7b5a1112396ad393f5588006c146c6954be07cba002f50e558ad6b83bb9990c" }, "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.4" } }, "nbformat": 4, "nbformat_minor": 2 }