Skip to content

Instantly share code, notes, and snippets.

@ischurov
Created November 27, 2020 15:32
Show Gist options
  • Save ischurov/17c8a1bb55021327ddebaaf7d8146370 to your computer and use it in GitHub Desktop.
Save ischurov/17c8a1bb55021327ddebaaf7d8146370 to your computer and use it in GitHub Desktop.
numerical integration of ode in Python.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "import matplotlib.pyplot as plt\nimport numpy as np\n%matplotlib inline",
"execution_count": 1,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "from scipy.integrate import odeint",
"execution_count": 2,
"outputs": []
},
{
"metadata": {},
"cell_type": "markdown",
"source": "$$\\dot x = y,\\quad \\dot y = x^2-1$$"
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "def f(X, t):\n x, y = X\n return np.array([y, x ** 2 - 1])",
"execution_count": 3,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "T = np.linspace(0, 3.5, 200)\nfor x0 in np.linspace(-3, 0, 11):\n trajectory = odeint(f, np.array([x0, 0]), T)\n plt.plot(trajectory[:, 0], trajectory[:, 1])\n plt.xlim(-5, 5)\n plt.ylim(-5, 5)",
"execution_count": 4,
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": "/Users/user/miniconda3/envs/featurevis-py3.6/lib/python3.6/site-packages/scipy/integrate/odepack.py:247: ODEintWarning: Excess work done on this call (perhaps wrong Dfun type). Run with full_output = 1 to get quantitative information.\n warnings.warn(warning_msg, ODEintWarning)\n"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOydd3gU5dqH79maZHdTN72HhN4JvQWk2I9djxVR8LOLvRyP6Dl27KKi2EAROYKKFaVKDR1CC0kgvddt2Wyb749BMAbIBhKSwNzXtdfCzDvvvJtkf/vs8z5FEEURGRkZGZnOi6K9FyAjIyMjc3rIQi4jIyPTyZGFXEZGRqaTIwu5jIyMTCdHFnIZGRmZTo6qPW5qNBrFhISE9ri1jIxMJ6PW5qSgxkaAr5q4YL+mA2xVUJsP+nDwjwLA5XGRU5uDWqkmMSARAeEMr7pt2LZtW6UoiqF/P94uQp6QkMDWrVvb49YyMjKdiNWZ5dz++VYuTwjis1uH4KNWNh5QuA0+PR/iL4EbF4NCiSiK3L3iblSlKr6++Gu6BHZpn8W3AYIg5B3vuOxakZGR6ZBsyK7kjvnb6BZh4MObU5uKuKUcvr4RDBFw1SegkM4vOLCAtUVrmTFoxlkl4iejXSxyGRkZmZOx+XA1t32+lUSjjvm3DcXfR914gMsBi26B+hq47TfwCwYgszqT17e+zpiYMVzf/fp2WHn7IAu5jIxMh2J7fg23frqZqEAf5t82lGCdpumg356C/A1wxVyI7AtAvauex/54DH+tP/8Z+R8E4ezwi3uDLOQyMjIdhozCOm75ZDOhBi0Lpg0j1KBtOmjHl7D5Qxh+D/S9+ujh17a+Rk5dDnMmzCHYJ/gMrrr9kX3kMjIyHYJ9xSZu+iSdAF81C6YNI9zfp+mgom3w4wxIHAMTnj16eEX+Cr7O/JopvaYwInrEGVx1x0AWchkZmXZnX7GJmz5Ox1et5Ktpw4gK9G06qK4IvrpeCjO86jNQSg6FMmsZz2x4hh7BPbhvwH1nduEdBFnIZWRk2pUd+TVc9+FGtCoFC6YNI/Z4seIOKyz8JzgscP1C0IUA4Pa4eXLdkzjcDl4Z8wpqpbrptecAspDLyMi0GxtzqrhxbjpBOg2L/m84iUZd00EeD3x7B5TslsIMw3sdPfXJnk/YXLqZJ4Y8QUJAwplbeAdD3uyUkZFpF1ZllvN/87cRF+zHl7cPJex4PnGAVf+F/T/ApOeh6+Sjh7eVbePdne9yQcIFXJZ82RladcdEFnIZGZkzzi8ZJdy3cAddww0nDjEE2LUQ1r4GA2+G4XcfPVxjr+HRPx4l1hDLv4f/+5wKNTwespDLyMicURZtKeDxJbsZEBfEJ1MGE+B7Ar92fjosvRcSRsOFr8ERsfaIHp5c9yS19lpmXzQbvUZ/BlffMZGFXEZG5owgiiJvLM/i7RVZjE4x8sGNg9BpTyBBVTnS5mZADFwzD1THLPbP9n7GuqJ1/Gvov+ge3P0Mrb5jIwu5jIxMm+NweXhiSQaLtxdy9aAYXriiD2rlCWItLBXwxZXSv2/45mj6PcCO8h28vf1tJsVP4ppu15yBlXcOZCGXkZFpU0x2J3d9sZ112ZXMmNCV+85LPrFP22GFBdeAuRSm/Aghx4pe1dpreWTNI0TqIpk5YuY57xf/K7KQy8jItBnFtfVM/WwL2eUWZl3dj6sGxZx4sNsF/5sCJTvh2i8hJvXoqT/94tX2auZfOB+DxtD2i+9EyEIuIyPTJqQfquLuBduxOz18eutgRqc06YdwDFGEn2ZA1m9w0evQ/cJGpz/a/RFri9by5NAn6RXS6wSTnLvIQi4jI9OqiKLI5xty+e9P+4kL8WPh9EEkhzVjQa95BbbPg9EPweDbGp3aULyB2Ttnc1HSRVzX7bo2XHnnpdWEXBAEJbAVKBJF8eLWmldGRqbzYHe6efLbDJZsL2JCj3Bev7Zf01rifyd9Dqx+AfpdD+OfbnSqxFLCY388RpfALvx7mBwvfiJa0yK/H9gP+LfinDIyMp2Egmobd325nYyiOmZM6Mq945NRKJoR3p0L4JdHofvFcOk7R2PFARxuBw+teQinx8kbaW/gpz5ODRYZoJWEXBCEGOAi4HngwdaYU0ZGpvPw/c4i/vXtHgDm3pzKhJ7hzV+0byl8fzckjoUrPz5azfBPXtnyChmVGbyR9sY5XUfFG1rLIn8TeBQ4oSNMEITpwHSAuLi4VrqtjIxMe2K2O3nm+70s2VHEoPgg3ry2//GrF/6dnJWw+DaIToXrFoC6cZ2VH3J+OFpffEL8hDZa/dnDaQu5IAgXA+WiKG4TBCHtRONEUfwQ+BAgNTVVPN37ysjItC87C2q576sdFNbYeGBCCveMS0Z1oiSfv5K/CRbeAMZucMMi0DZOsd9TuYeZG2YyOGIw9w+8v41Wf3bRGhb5SOBSQRAuBHwAf0EQvhBF8cZWmFtGRqaDUe9w8+aKg8xde5gIfx8W3TGc1AQvW6vlp8MXV4F/FNy0BHyDGp2usFVw/8r7MfoaeW3sa6gUcmCdN5z2T0kUxSeAJwCOWOQPyyIuI3N2siG7kie+zSCvysZ1g2N54sIeJy569Xfy06XUe30Y3PKj9PwXHG4HM1bPwOw0M/+C+QT5BJ1gIpm/I3/cycjINEutzcELP+9n0dZCEkL8WDBtKCO6GL2f4K8iPuUn8I9sdFoURZ5Pf55dFbuYNXYW3YK7tfIrOLtpVSEXRXE1sLo155SRkWk/nG4P8zfm8daKLCwNLu5M68L956Xgo1Z6P0kzIg6wMHMhS7KWML3vdCYnTD7OJDInQ7bIZdoVp9uD3elGBJSCgFIhIAigUSrk5I92RBRFlu8v58Wf93Oo0sqoZCP/urgH3SNamCZyeC18dZ3UMPkEIr6ldAsvb36ZtJg07u5/93EmkWkOWchlWh23R6S4tp5DlVYOVVgorKmn0tIgPcwOqm0O6h1u7E43Ls/xA5gEAXQaFXqtCp1WicFHTahBS6hBS9iR51C9ljB/H6IDfTHqNbLwtwKiKLLxUBVvLc8i/XA1XUJ1fDIllXHdwlr+8838BRbdAsFJcNO3xxXxIksRD65+kHj/eF4c/SIKQW4jfCrIQi5zWng8IocqLWzPr2VnQS0782vJrrDgcHmOjvFVKzEaNBj1WuJD/BgYH4ivWoWvRoGPSomPWokgSB8AblFEFKHB6cbS4MbS4MTa4Kau3klBtY1teTVUWx1N1qHTKIkN9iM+xI/4EB1xwX4kGXWkhBtkkfcCURT5I6uSd1ZksTWvhlCDlmcv7cX1Q+NOXDf8ZOz6Gr67E6L6N6kp/ic2p437Vt6HW3Tz9vi35U4/p4Es5DItpsrSwJqDFazKrOCPgxXU1TsBMPio6B8byKiUBJKMOpJC9SQada0upE63h0pLAxXmBspMDRTW2MirspFfbSO73MKqzIpGHyRBfmq6hhuOPPRH/x10oj6R5xB2p5ufM0r4bEMuuwvriArw4bl/9OKa1NiW+cH/SvqH8MsjkDhGSvbRNs0TdHlcPPLHI+TU5jD7vNnE+8ef5is5t5GFXMYrqq0Olu4s4vtdxewsqEUUwajXMrFnOEMTgxkQF0SSUdd8bY1WQK1UEBngS2SA73HPezwipSY7hyqsHCwzk1VuJrPUzHc7ijA3uI6OM+q19Ig00D3CQPcIf7pHGkgO06NVnaKAdSIKqm18mZ7Poq0FVFsdJBl1vHhFH64cGINGdYruDVGElf+FtbOg20Vw1SdNMjalYSIvb36ZPwr/4OlhTzMyeuRpvhoZWchlTojL7WHlgXK+2VbIqsxynG6RnpH+zJjQlXHdwugV5X9GhLulKBQCUYG+RAX6MirlWIicKEoCn1lqJqvMQmaZJPCfb8w7asErFQJdQnVHhb3HkecIf59O756ptDTwS0YJP+wqYXNuNUqFwMQe4dw0PJ4RXUJO7/W5GqS6KRn/kzreX/RGk9opfzJ/33wWZi5kSq8pcru2VkIWcpkm2J1u/re1gI/WHia/2oZRr2XKiASuHBTT8qiFDoQgCEct+bRux5JRXG4PuVU2DpSaOFBi5kCpiW15NSzdVXx0TICvmu4RBnpE+ksWfKQ/XcP1+Gk67ltIFEWyyy38kVXJqgPlbMipxCNCSpieByd25erUmBN+q2kRtmop5T5/A5z3DIya0aiK4V9ZkbeCWVtnMTF+IjMGzTj9e8sAIIjimS97kpqaKm7duvWM31fm5FgaXHy67jCfbcilyuqgf2wg/zc2iQk9wr2roXGWUVfv5GCZmQMlJvaXSs+ZpWasDjcgaVVCiI5u4Qa6H3HRJBr1xAX74as58+4Zl9tDZpmZ3YV1bM+rYV12JSV1dgCSQnVc2DuSS/pF0S2iFdukVeVIPTZr8+HyD6D3lSccmlGRwdRlU+ka1JWPJ3+Mj6qp20Xm5AiCsE0UxdS/H++45oTMGcPtEVm0tYDXfsuk0uJgfPcw7hiTxJDE4E7vTjgdAnzVDE4IZvBf6oh4PCKFNfXs/4v1nllqZtm+Uv5qE4X7a4kP1hEX4kdcsB/h/lrCDD5HwyeDdZpT+nB0e0SqrQ4qzA3kV1s5VGnlcIWVnAoL+0pM2J2SiyjQT83ILkZGpxgZlWIkJqgNankfXAaLp4FCCTcvhfjhJxxaZCninpX3EOIbwlvj35JFvJWRhfwcZ312Jf/5cR8HSs2kxgcx95bB9I8NbO9ldVgUCkES5xA/JveKOHrc5nCRVWYht8pKfpWNvGob+VU2/jhYQbm5oek8AgT6aY7EyavQa6UwTIUgoBBAIQi4PCL1Tine3uZwU2tzUG118PfQ+zCDlkSjjhuGxtM3JoD+sYHEBfu13YewxwNrX4NVz0NEb6lRctCJo07qGuq4e/ndOD1OPj3vU4y+LUjtl/EKWcjPUUx2J8/9sI9vthUSG+zLezcM5ILeEee0BX46+GlU9IsNpN9xPgTtTjcV5gYqLA2Um6TnCpOdKqsDa4MLS4Mbm8OFtcGFRwSPKOL2iKiUCnxUCoJ1GqIClATp1ITqtRgNWox6LXHBfiQYdei1Z/BtbK+D7+6CAz9Cn2vgkrdAc2Jr3+a0cc+Ke8g35/PBhA9ICkw6c2s9h5CF/BxkXVYlj36zizJzA/eMS+ae8cmnHjMs0yw+ailZyauGCx2Zgi2weCrUFcHkF2HYnSfc1ARwup08uOZBdlfuZtbYWQyJHHIGF3tuIQv5OYTD5eH5n/bx+cY8uoTqWHznCNmNItM8Hg+sf1NypfhHwdRlEDv45JeIHp5a/xTri9bzzPBnmBg/8Qwt9txEFvJzhHKTnTu/3M62vBqmjkzk0fO7yVa4TPPU5sP398DhNdDzMsmV4nvyD39RFHkx/UV+OfwLDwx8gKu6XnWGFnvuIgv5OcDOglrumL8VU72Ld68fwMV9o9p7SVIWoMclJZK4HdJDUEoREEo1qHxPmFAicwbweGDbJ/D7M9L/L3kLBt5yUlfKn7y/6/2jCT9Te09t44XKgCzkZz1LdxXz8KJdhAdoWXLXCHpEnqGEHrcLKg9Kj5rDUH0Y6grAWiElkNiqwGU/+Rwag9QKzDfgyHMwBMRAQCwExh579pU7ybQqVTnww/2QuxaS0uCSt08alfJXvtz/Je/vep/Lki/jwUEPypvnZwhZyM9ivtqcz5PfZjAkIZgPbhzUtkWibNXSGz93HRTvgNI94Ko/dt7PKImBfwxE9ANdiCTUKq30UKhA9EhWuscFDhvYa6G+5tijNAMO/tr0A8DPCKHdIbTbkeeuEN5HuoeM9zissPZ12PA2KLUtssIBvjrwFS9tfonz4s7jmeHPyCJ+BpGF/Cxl7tpD/Pen/YzrFsr7Nw5qG394+QHYuwSyfoPinYAIGj1E9oPUqdJzWA8ISgCfVvomIIqSVV9bIFn4tXlQmQUVmZDxDTTUHRsbECutIao/RA6QnnVyDHMTRBH2fgu/PQ2mQuh7LUx8DgwRzV97hK8PfM0L6S+QFpvGq2NelZsmn2Hkn/ZZyNsrsnj994Nc2CeCN68dcOrV7I6HuQx2zIc9i6F8HyBA3DBIe0L6Gh49UPJxtxWCILUM04dBzKDG50QRLGVQvl+y3kt2Sh8wB348NsY/BmJSIXYoxA2FiL5tu96OjChCzkpY8Zz0swrvDVd+BPEjWjTNosxF/Df9v6TFpPH62NdRn6s/z3ZEFvKzjM/WH+b13w9yxcBoXrmyb+vUSBFFKNwKm+fA3u/A44TYYXDBq9DzH2AIP/17tAaCIFmRhgjoMu7YcXsdlOw+Iuw7oGAz7PtOOqfyhehBEDtEEvfYIcdtgnBWIYpwaBX88RrkrYOAOLjsfckSV7Tsm9s3B7/hP5v+w5iYMbyW9pos4u2ELORnET/tLuHZH/cxqWc4r17VD2VrlJjN2yjFD+euBa0/DL5dehiTT3/uM4VPACSOlh5/UlcEhZulxsAF6ZJf2HOkVrmx6xFRPyLsISmgOAuKhrmdkgtl/dtQlgH6COnDeNAt0j5FC1l8cDHPbnyW0dGjeSPtDTRKuVFHeyFXPzxL2JhTxS2fbKZvTABf3D709H3iJbth+UzIWQG6MBj9IAy48bjdXs4KHDYo3g75mySLvSBd2mwF8Ak8YrEPgZghkgWv7URtySqzYcc82LlA2l8wdoMR90Lfa05JwAE+3/s5s7bOYlT0KN4c9yZa5anNI9My5OqHZzHZ5Ramz9tKXIgfc29JPT0Rt5tg1QuSG8UnUNr0GjztpPU0zgo0fpAwSnqAFEddlS0JeuFmSdyzfpPOCQrJn/ynxR47BALjvY7uOCPU5MH+H2Df99L6BSV0PR8GTYHkCaf8DUMURd7Z8Q4fZXzE5ITJvDjqRdmd0gGQLfJOjs3h4rLZ66m0OPjh3lFEB55Go4D9P8DPj4C5FAbfBuOfbjaL75yivgYKt0niXpAORdvAYZHO+QZBRB9p8zSij/Qwdj1zG6kOK+RvhEOrpUdphnQ8og/0ugL6X9+iKJTj4fa4eSH9BRYdXMRVXa/iX0P/hbKFPnWZ00O2yM9CRFHkySUZZJVbmDd1yKmLuMMGy56AbZ9Jb/xrv2waESIjiXXKBOkB4HFLkTsFm6F0t+SO2jL3WJy7UgshXSA4SXoOSYagRKleiT781NwzbpcUcll9SPrGUJohbeBWHJDi8JUa6ZvChGeh56XSvVsBp9vJU+ue4pfcX7it923cP/B+OU68AyELeSdmweZ8vttZzIwJXRmdEnpqk5Ttg29ulYRg5P0w7l+gkjetvEKhPGZ9/4nbdUxgS3dLMe6VB6UmDB5n4+s1eknQfQJAowO1r/RQaqSSBS4HuBuk5/pqsJRLGbH85Vu0LhSiBkD3i6Uw0Ljhre4Gq2uo46HVD5Fems6Dgx7k1t63tur8MqePLOSdlH3FJp5duo8xXUO5d/wpRpBk/grfTJVE5MYlkHxe6y7yXESpgrDu0qPv1ceOu11SAlPNYSkW31IqCbO5FBpM4KyXNiKd9VL9GaVG+kBVHsl8DU6SLG19GATGQXAXycrXhbapbz7flM/dK+6myFLEC6Ne4JIul7TZvWROHVnIOyEut4dHF+/C31fFm9f2b3kne1GE9A9g2ZOST/efC8E/sm0WKyOhVEFwovToJGwt3coDqx9AQGDupLkMDB/Y3kuSOQGykHdCPlx7iD1FJt67YSDBLa2fIoqw7CnYNFv6On7Fh5JF3obUu+rJrcul0FJIkbmIivoKzA7z0YfNZUNAQBAElIIShaBAr9ETqA0kUBtIgDaAcL9wovXRROmjCPMLQyGcBXHdHZhvs77luU3PEWuIZfb42cT6x7b3kmROgizknYzscgtvLs/igt4RXNinhVa0KMIvj0mhhUPvhMkvtHqiiyiK5Jpy2VK6he3l29lftZ9cUy4e0XN0jK/KF4PGgL/G/+gzgFt0I4oiTo+TYksx+6r2UddQR4O7cc9LtUJNrCGWlKAUUgJTSAlKoXtwdyJ1kfIG3Glid9l5afNLLM5azLDIYbyW9trR349Mx0UW8k6ExyPy2OLd+KqVPPuPXi27WBThl0dh84dSMsjE/7Sab9XtcbO9fDvL85azIn8FZbYyAEJ9Q+ll7MWkhEmkBKYQa4gl2hDdYmGod9VTai2l2FJMkaWIQkshh+sOs6dyD8tylx0dZ/Q10i+039FHz5Cecrf2FlBgKuDBNQ9yoPoA0/pM467+d8nFrzoJ8m+pE/H9riK25dXw6lV9CTO0UKCWz2x1ES+2FPNt9rcsyVpCua0crVLLyKiR3NHvDoZEDCHOENcqFrKvypfEgEQSA5r6l61OK9m12eyr2sfuit3sqtjFivwVAKgUKvoa+zI0cijDIofRJ7QPaoWcvHI8luct5+n1T6MQFMw+bzZjYsa095JkWoCcENRJsDvdjJ+1mhC9lu/vHtmyDc6tn8KPD0ilZS96/bRFPLM6k7kZc/kt7zdEUWRk9EguS76M0dGj8VO3fwZoVX0VGZUZbC/fzuaSzeyr2oeIiK/Kl0HhgxgVPYoxMWOINch+X5PDxMubX2ZpzlJ6h/TmtbTXiNJ3gA5SMsflRAlBspB3EmavyubVZZl8NW0Yw7u0oGFC1nJYcI0UWnjdV6fVPu1w3WHe2PYGqwpWoVPruLbbtVzb7doO/8ava6hjS+kWNpVsIr0knVxTLgBJAUmMjRnLmJgx9A/rf865EdYXreffG/5NVX0V0/pOY3qf6XK6fQenzYRcEIRYYB4QAXiAD0VRfOtk18hC3jIqLQ2kvbqaYUkhzL2lye/wxFTlwJyxEJwAt/5yygWvzA4z7+18j4UHFqJVaZnSawrX97i+026C5Zvy+aPwD9YUrmFr2VZcHhf+Gn9GRo9kbMxYRkWPIkAb0N7LbDOq6qt4a/tbfJv9LV0CuvD8qOfpZfRuz8X06zK0XZLQpqS08SpljkdbCnkkECmK4nZBEAzANuAyURT3negaWchbxjPf7+GL9Hx+mzGGLqFepnU77fDxBKgrhDvWSr0tT4H1Ret5ZsMzVNRXcEXKFdzd/26MvmdPlx2Lw8LGko2sKVjD2qK1VNurUQpKBoYPZGzMWMbFjiPOP669l9kquDwuvs78mtk7ZlPvquemXjdxd/+7va5caN24kfzpd2BIG0vMO++08Wpljkeb1VoRRbEEKDnyb7MgCPuBaOCEQi7jPRXmBhZuKeCqgTHeizjA709LaeL/XHhKIt7gbuCVza+w6OAiEgMSmZ82n76hfVs8T0dHr9EzMX4iE+Mn4hE9ZFRmsKZgDasKVjFr6yxmbZ1FYkAiabFpjIsdR19j305XKEoURdYWreWNbW+QXZvN8MjhPD7kcZICva/DYt+/n8J770ObkEDk88+34WplToVWdQoKgpAADADSW3Pec5nPNhzG4fZwx9gWFD868LMUoTLsbuh2QYvvWWwpZsbqGeyr2seUXlO4Z8A950S9aYWgOBq6eN/A+yg0F7KmcA2rC1Yzf+98Pt3zKUHaIEbHjCYtNo0RUSPQqds2mep0EEWRdUXreH/X+2RUZhCtj+bNtDcZHze+RdFEjoIC8qdNR6HXE/vhHJT+ndOldjbTapudgiDogTXA86IoLjnO+enAdIC4uLhBeXl5rXLfsxmz3cmIl1YyOsXIezd4WY2wvhZmD5VqcExb2eICWDvLd3LvyntxeVw8P+p5xseNP4WVn32YHWbWF69ndcFq1hauxeQwoVaoGRIxhLTYNNJi04jQnV6Z2NbC7rLza+6vLDywkL1Ve4nSRTG973Qu7XJpizczXZWV5F5/A566OuIXfIm2S5c2WrWMN7Rp1IogCGrgR2CZKIqvNzde9pF7x5w1Obz4ywGW3jOSvjFe1gX/4QHY/rkk4lEDWnS/VfmreOSPR4jQRTD7vNnE+8efwqrPflweFzvKd7CmYA2rC1eTZ5KMkm5B3RgeNZzU8FQGhA84o5vBoiiyv3o/Px36ie+yv8PkMJEUkMRNPW/iH13+cUrRKG6Lhbybb8ZxOJf4Tz/Bt3//Nli5TEtoy81OAfgcqBZF8QFvrpGFvHmcbg8jX1pJ13ADX9w+1LuLctfDZxfC8Htgcsv8mD8d+okn1z1Jz+CezJ4wm2Cftm1ALHo8uGtq8NhseOrrERsaEJRKBI0GQatF6e+Pwt+/U6TcH647zOqC1awpXMPuit04PU4UgoJuQd1IjUhlYNhAugd3J1of3aqvx+a0sbtyN+uL1vN73u8UWYpQCSrGx43nuu7XkRqeesr38zQ0UDD9DmzbthH73mz0Y+QEoY5AWwr5KGAtkIEUfgjwpCiKP5/oGlnIm+fXPaX83xfbmHtzKhN6etGl3uOGD0ZJHWvu2tSiQlgr8lbw0JqHGBA2gNnnzW7VpB7R7aYhOxv7nr3Y9+6lITsbZ0kJztJScDpPeq2gVqMMCUEdEYEmIUF6dEnCt3dvVBERHVLk7S47GZUZbCndwtayrewq34XD4wDAoDHQPbg7XYO6EmuIJUoXRZReeujV+hO+HrvLTo29hgJzAYfqDnGo7hB7K/eyr2ofLtGFSqFieORwJsZPJC02jSCfoNN6DaLTSeGMGViWryDq1VcIuEQuXdtRaMuolXVAx3tHdXK+3pJPuL+WtG5eNozY9ZXUrebqz1ok4huKN/DwHw/Ty9iLd897t1VE3G0yYVnzB5Y1a7CuXYu7rg4AhU6HtmtXfPv1w//8yajCwlHodCj8fBG0WvB4EBsa8DQ48JjqcFVW4qqoxFlSgnXjRuq+++7oPZShRnx798G3X190w4bh07s3gqr9E3p8VD4MjhjM4IjBgBT9c7D6IPur95NZncmB6gMsyVpCvau+0XUKQYFOrZMEHQEPHjweD2anuclYP5Uf3YK7MaX3FAaFD6J/aH/0mtZpBi263RQ/9hiW5SsIf+opWcQ7Ce3/ly/ThOLaetYcrODuccmolF5UJ3TWw8rnpe7uPS/z+j6H6g7x0OqHSApI4v0J759WBIYoitTv2EHt14sw/forYkMDyuBg9Glp6EaOwKd3HzQJ8QinUW3RY7XSkJ1NfcYe7Bm7qc/Yg2XVKip4C4XBgG7YUHQjR6IfNx51eNgp36c10Sq19AntQ5/QY12ERFGkpuEtPUIAACAASURBVKHmaBGwUmspdQ11WJ1WLE6pB6iAgFKhRK/WE+QTRJA2iCh9FIkBiYT7hbfJtxHR46Hkyacw/fwLYY88TPBNN7b6PWTaBlnIOyD/21qIR4RrUr2M/970PpiL4cqPvK6jUtdQx70r7kWj1PDu+HdPeWNOFEWsf/xBxVtvY9+3D4VOR8DllxF42WX49O17WsL9dxQ6Hb79+uHbrx9wAwCumhpsGzdi2bAB64YNmH9fDjOfxbd/fwyTJmGYNBFNTEyrraE1EASBYJ9ggn2C6W3s3d7LAaTfY+nMZ6n7/nuM991LyG23tfeSZFqAXGulg+HxiIx+ZRVJoTrm3+bFJqfDCm/0hphUuOF/Xt1DFEXuXXkvG4o38MnkT+gfdmrRCLbtOyh/5RXqd+5EHRNDyPRpBFx0EQpd+8RWi6KIIycH8++/Y/r9dxr27QfAN3UQgZdfjmHy+Sj1HTfuu70QRZGyF16kZv58Qu64g9AH5MbKHZU285HLtC47Cmopqq3n4cldvbzgC6kx7+iHvL7HVwe+Yk3hGh4f8vgpibjbYqXijTeoWbAAVVgYEc8+S+AVlyOo27fgkiAIaJOT0SYnY7zzThwFBZh+/oW6776j5Kl/Ufrf5zFMnEDglVfhN2SwLFZIIl7+0svUzJ9P8JQpsoh3UmQh72As21uKWikwvrsXkSpuF2x4V2rKGzfMq/mza7J5betrjIkZw/Xdr2/x+mzbtlH0yCO4SkoJuvFGwh64v90s8ObQxMZivGM6IdOnYd+1i9rvvsP008+Ylv6ANiWFoBtuIOCSizvs+tsa0eOh9LnnqF34tfS7fOxRWcQ7KXLjww6EKIr8uqeUEV2MBPh6Yd3u/Rbq8mGkV+H7eEQPMzfORKfW8dyI51r0phVFkep588i7ZQqCWk38l18Q8dSTnUIEBUHAt39/ImfOJGXtH0S+8AKCWk3pzJlkpY2j7MUXcZxjmcaiy0XJE09Qu/BrQqZNI/ypJ2UR78TIQt6B2F9iJr/axgW9vUz13vIRBHeBrud7NXxx1mJ2Vezi4cEPE+LrfU1z0e2m7D//oeyFF9GPHUviN9/gN7BzdlRX+PgQeMXlJCz+hvivFqAfM4bqLxeQc/4FFN53P/UZGe29xDZHdDgoeuhh6r5fSugD9xP20IOyiHdyZCHvQPy6pwSFgHcJQBWZUJAOg27xqoFyrb2WN7a9wZCIIVyS5H1ssOh0UvTgQ9Qs+IqQ228j5p23URpOra55R0IQBPwGDCD6tVkkr1xByPTpWDdtIvfqa8ibciuW9etpj0CAtsZjs1Fw772Yly0j/InHMf7f/7X3kmRaAVnIOxDL95eTGh+MUe9FpcEd80Ghgn7/9GrujzI+wuq08sSQJ7y2vkSXi6JHH8W8bBlhjz1G2MMPt2o4YUdBHRZG2IwHSF65grBHHsFx6BAFt91O7pVXSTHxbnd7L7FVcFVUkHfTzVjXriPiuWcJvuWW9l6STCshb3Z2EGqsDvaVmHhoohfRKi4H7FoouVT0zSe+lFhK+OrAV1za5VKSg5K9Wo8oipT+97+Yf/mVsEcfJeTWKV5d5y0OuwtztZ16kwOXw4PL6UH0iKh9lGh8Vfjq1RhCfFCpz1ztb6VeT8htUwm66UZMS5dSNfdjih6YgSYxkZDp0wm4+KJ2j8w5VRqysymYfgeu2lpi3puNIS2tvZck04rIQt5BSD9cBcCIZC9814dWg7UCBniXeTdn9xwEBO7uf7fX66n+7PMjG2G3EzL1Vq+vOx6iKFJZaCF/bxVlh02U55mx1jY0f6EA+kAtgeF+hCf4E57oT2SXQHz0bSumCo2GwKuuIuDyyzH//juVH8yh5IknqHz3XUJuv42AK65Aoe089dmtm9IpvPdeBB8t8fPm4dvbu7ZuMp0HWcg7CBtyqvDTKL0rV7t/KWj9oUvztcIrbBUszVnKFSlXeF0v27ZlC+Wvvoph0iRCZ8zw6prjUVdRz/4NxWRuKsVSIwl3YLgf0V0DCY7S4R/ii1+ABpVGiUqjQEDA0eDCaXdjMzkwVdZTV1FPVZGF7b/lI3pEECA8wZ/43iEk9Q8lJLp1aowcD0GpxP/88zFMnoxlzRqq3v+A0mefo/K99wmeOpWga69B4dd6BcZaG1EUqf78c8pfnYUmIYHYOXPQxES397Jk2gBZyDsIG3KqGJIYjLq52ioeN2T+DCmTQNW8VbjgwAJcHhe39PTOH+qqrqbooYfRxMZKYXqn4BOvKrKw5adccnaUIwBxvUIYckkicb1C8GgV7Cs2sbPYRGFpJRVZDVRbHbg90saiWqkg0E9NiE7Dw5O7YfCRrG+nw01FvpmizBry9lSx+cfDbP7hMMZYPd2GRtB1SAR+/i1rouEtgiBgSEtDP3YstvR0Kj+YQ/nLL1M1Zw7Bt9xM0A03dLiuOR6rleJ//QvzL79imDiByBdfRKlvuw89mfZFFvIOQLnZTna5hasHeVETJH8j2KqgR/ORJ/Wuer7O/JoJ8ROI9feubkvZf5/HVVND4pwPWpzObrc62fRdDnvXFaPWKhk0OZ7eY6Mpd7n4YVcJK+YfIqOojj+DQfRaFaEGLcE6DSqFtAFrdbgoqq2n2urgsQu6H51brVESlRxIVHIggy9KxGZykL2tjMxNpaz/JpuN3+aQMjic/hPiMMa0jWAJgoBu2DB0w4Zh27GDqg/mUPHW21R9/AlB119P8JRbUAW3bR13b7AfOEDRww/jOHSY0IceJOT22+XwwrMcWcg7ADvyawEYnOiFCBz4CZRaSJ7Q7NDlecsxO8z8s7t3kS3mlasw/fwzxvvuxadHD6+u+ZO8vVWs+GwfdquLfuNjGXRBPBsKarht0Q7SD1cjCDAgNpD7z0uhX0wgvaMDCDWcup/Zz19D33Gx9B0XS3WxlT1ri9i/oYTMTaXE9ghiyCVJRCQFnPL8zd5/wAD85nyAff9+Kud8SNVHH1E9bx5B115D8NSpqMO9CCFtZUSXi6q5H1MxezbKgADiPp6LbvjwM74OmTOPXDSrA/D6b5nMXp3D3mcn49NclMZ7I0AfCjd/3+y8U5dNpcxaxo+X/9isReZxODh00cUofLQkLl6MoPHOTSF6RDb/eJitP+cSEq3jvCk9yfe4ePaHvewurCM60JcbhsVx1cAYwvx9vJrzVLFbnexdW8SuFQXUm50k9jMy7LIuBEe2ffZpw6FDVM35kLoff0RQKAi4/HJCpt2OJtbLCpaniT3zIKX//jf1u3ZhOP98Ip75N6qg02swIdPxkItmdWD2FJtIDtU3L+LWSijfC72fbnbOAnMBW0q3cN+A+7z6Wl0zfz7OggJiP57rtYi73R5WztvPwfQyuo+IZPAVXXhleSZfbMonMsCHV67sy+UDo5v3+7cSPjo1g85PoE9aDLtXFrD9t3wWPpdO33GxDLk0EY1P2/25a5OSiHr5JYz33kPV3LnULV5C7eLF+F90Icbp09Emexf22VLctbVUvDubmq++QqnXEzVrFv4XXSi7Us4xZCHvAOwpqmNUirH5gbnrpOfE5vsn/pb7GwAXJV3U7Fi3xUrVhx+hGzMa/ciRza8D8Lg9LPtwD4d3VTL0H0kEDQrh6o82cbDczG2jEnlwYld02vb589L4qEi9MJFeY6JJX3qYXasKyN5ezpjrupLU38uOS6d675gYImfOxHjnnVR/8ik1ixZh+uFHdCNGEHDF5RgmTGiV0EVXTQ3V8+ZRM/8LPDYbQdddS+h996EM9LJJt8xZhSzk7Uy5yU65uYHeUV74c3PXgloHUQOaHboyfyW9QnoRpY9qdmzt1wtx19URes893iwZURRZ/WUmh3dVMvraFNxd9Fz1wUYEQeCzW4cwtmvzYlnndLHHUk+B3UGFw0XDkagVg0pBuEZNlFZNN50PAepT/xP11WtIu74b3YdFsHpBJr98kEH34RGMvrZrm1rnAOrwcMKfeJyQO6ZT8+UCapcsofihh1H4+2OYNBFDWhq6ESNaFL4oulzYtm2nbsmSo12YDJMnY7z7Lny6eln2WOasRBbydmZvsQmA3tFeCHlBOsQOAeXJE2LKrGXsrtzNfQPua3ZK0eWi+vN56EYMx7dvX6/WvPP3AvZvKCH1wgScSXpumpuOUa9l/m1DiA85sT96p8nG0vJafq+qI8vmRUIQcKExgE/6JHo19kREJAVw9ROpbP0pl22/5FKcVcuk23sTntD2IYOq4GBC770H4113YktPp/bb7zD/uoy6b6R9CN++ffHp0wefnj1QR0SgCg+XXFuiiMdqxVlSSkNWFvY9GVg3bMRdW3u0C1PQ9dfLAi4DyELe7uRUSD0au4Y3EzLnaoDy/TCieXFeVyS5YMbFjmt2rGXNGlzl5UTMfKb5xQKlh+rY9F0OSQNCCRpi5Oo5mwgzaFl0x/DjbmaKosivlXW8k1/OdpMNtSAwIlDP1RHB9NZpCbVb0dttCDYrCqWSBo0PFr0/VWof9lvthGha509UqVQw9NIkYnsG8/sne/l21nbG3diNbsMiW2X+5hCUSnQjRqAbMQLR6cS2bTuW1aux7dhOzZdfIjocJ71eFRWJbsxoDOPPQz96VKcoHyxz5pCFvJ3JrbLi76Mi0K+ZDcbyfeBxQWS/ZufcXLoZo6+RLoFdmh1bu+h/qMLC0I9p3u/ucrpZ/uk+dEFahlyTzFUfp6NWKpg3dehxRTzbZufRzEI21FpI8tXyfEo0aQ4zZds2krt7Bxm5h3E5jm+Za3z9CI1PJK53X4r6DiQqpVurFOyKSg7kmicHs+zDPSz/bD9VRVaGX94FQXHmNgcFtVpqFD1MauUnOp048vJwlpbhKi9HdDlBEFD4+qGOjEATH4/K6MUeisw5iyzk7Uxele2k7oijlOySnpsRclEU2VK6hcHhzbcyc5aWYlm7lpA7piOomv9T2PZrHnUV9Vz6QH+eW3aAvCorC6YNIy6kqZ93UWk1j2UW4qMQeDElmhGlh9j2wVy+y8oEQSAyuSt9J5xPSEwsfv6BaHU6PG43Dns95spKakoKKc0+yKbFX7Pxm6/wDw2n19jx9J90EX4Bp7eh56vXcMn9/Vn3dRY7fs+n3uJg3I3dUZyh6Jq/I6jVR1vUycicCrKQtzP51Tb6eOMfL9kN2gAISjj5fOZ8KuorSI1oEmraBPOKFeDxEHDJpc2OtdTY2bEsn5TB4RwUnXy/s5gZE7oyLKlxkS9RFJmVW8pruWUMD9TxUqgfuz5+nR/3ZRAYHknazdPoMToNP3/vknXsFguHtm9m39pVbFy8kC3fL6bvhPMZdtU/8dWfel10pVLBmH92xS9Aw+YfDuOwu5l8e692E3MZmdNBFvJ2xOn2UFRTz8V9vfDTVmVBaFdoxsreW7kXgH6hzbtgLKvXoImPR5vU/Gbi9mVS0aoBF8Vz+eebSQnTc2daU9fNy4dLeTOvjGsjgrmjroBlT72JQqHgvNvuos/4SSi9sPz/io9eT88x4+k5ZjzVxYVs/u4bdiz7kX3rVjPm+in0HjfxlGOmBUFg8EVSfPm6/2Wxct4Bzrulxxl1s8jItAay+dGOlNTacXlE4oO9cK1U50JQ84K7v3o/GoWGpMCkk47zWK3YNm1C70VdapvJwb51xXQfHsHPhyoprKnnmUt6oVE1/vP5uqSaN/PKuCEymFvyM/jlzZcwxsZx8yvv0H/ShScUcbfbRoOjkoaGCjwe1wnXERwVw/l3PcBNL7+NMSaO3+a8zdLXXsBusTT7Gk5Gv/NiGXppEpnppWxYkn1ac8nItAeyRd6OlJntAEQENJO67moAUyEEeyHkVfvpGtQVteLkIYrWLVsQnU70aWObnfPAphLcLg89x8Xw+Lx0hiQEM/JvddOzrHYePVjA6CA9U8pz+O2T9+mSOpSL7n8UtaZxAozLZaGs7Ecqq1ZhMu3G4Sj/y1kFPj6R+Pv3IzHhHvT6bk3WExqXwDX/foFtP33H2q/mseDph7nyiZkEhHnZ6/Q4DLogHpvJwc7lBRhjDXQbeupzycicaWQhb0eqLFLIWbCumYiV2nwQPV5Z5Fm1WaTFpjU7zr57NygU+PZrfvN0//oSIpMD2FJtpszUwEtX9G3kzvCIIg8eKMBPoeBZvciyV94itmcfLn7gcVR/6ajjdtvJy5tDfsEnuN0WfHyiCQkehZ9fEiqV5O9ucFRgsx3GZNp50nUJCgWpl1xBRHJXvn/1v3z19CNcO/MlgiJPrd62IAiMvDqZqiILq784QGisgeAoOcRPpnMgC3k7Um2VhDxE35yQ50nPzWx0Wp1Wqu3VxBqaL9RUvzsDbUpKs5mFVUVWasts9J/QjRc3FxAd6MuYv2VuLi2vZYvJyhsp0Wx+4xl8/f255MEnGom41ZrN7oz/w2Y7TGjo+cTHTcPfv99p1wSJ6dGba599mUXPPsE3zz/NP597FX2wF12WjoNSqWDS7b1Y+NxmVs7fzxWPDEIh+8tlOgGyj7wdqbFJQh7UXAy5tVJ6bqY/Z4G5AIA4Q9xJx4miiD0jA58+vZtdY26GdO/AZH825FRy+YBolH8RN7co8lpuKV39fEjetYGqwnwmTb8XX8OxrEmTaTdbtl6J02liQP959O0zm4CA/q1W2MkYG8+VTz6HzVTHj2+9guc0miXrArSMvi6FssMm9qwpbJX1yci0NbKQtyNVFgc6jdK7qocAfie3NP8U8hjDyRtUuMrLcdfV4dOzZ7NrLNhXTWicgfSSOjwiTOzZuM72+hoLWbYG7okOZvO3XxPXpz9JAwcfPW+3F7Nz11TU6kCGDP6O4GDvinK1lPCkZCZOu4eiA3vZ/P03pzVXSmo40d0C2fpzLg77iTdfZWQ6CrJrpR2ptjYQ3JxbBcBWCQo1+Jw89rrCVgFAmN/JLXdnUREgVeo7GR6PSHm+mR7DI/nxUBVBfuomMe//K6vGX6UgOTuDArOJYZdfc/ScKIrs3fcwHo+T1EGf4ePTuICXPbsW65ZSHLl1uE0OEARURh+0XQLRDY5AE9WyTj89R48je8tG0r9dRM8x4/A3nvzncCIEQWDYZV1Y/PI29vxRxMBJ8ac0j4zMmUK2yNsRS4MbvdaLjvDWSskab8YVUdsgdRoK1J4889FZVAyAOvrkG4M1pVZcDW7CEgzsLKilf2xgI5+xRxT5vdLE+cYActauJCgyipiefY6er6j8jdradFKSH8fP79hGrcfuonL+PirnZtCQVYMmIQBDWiz60dGogn2xbimj/O0dVP/vIB5Hy9wkaTffjsftZsvSJS267u9EJAYQmRzAvrXFtEfzFRmZliBb5O2Iw+1pEot9XOx14Nt8Wnq1vZoAbQAqxcl/rX9a5OrIkyci1ZbZAPAP9yO73ML5vRqH5O232ql1uRmm11K4fy8DL7y0kd87L+9D/PwSiYy8+ugx0emmYm4GzmIL/ucnYBgVjfC3n4Gn3oVpdQGWPwpxllkJndYHhZe1zf2NYfQYNZY9q39n9PW3oPHx9eq649FrVBTLP9tP2WFTm7aNk5E5XWSLvB1xuNxovUkJdzWAqvlmBKYGEwGa5gXHVVWFQq9vNmLFXCXFuZuUIh6RJjVhdpgkoU+oLsPjdhHX61gZ3Pr6fEymnURFXoPiLx8sdb/m4iy0EHJDD/zTYnG4nWzcuJFvvvmGJUuWkJWVhcJXReAFiYTc2BNnsYWaxVnNvqa/0nPMebgaGsjbvaNF1/2d+N5GEKDwQPVpzSMj09a0ipALgnC+IAiZgiBkC4LweGvMeS7gcHlpkbsbpIbLzWB329F6IfiivR6Fb/OWqqWmAZVGQUWDE4DooMbX5NU3oBJAUyptshrjE46eq6lJl44Zzzv2MiwOLJtK0A2OwLeXkZqaGt577z2WLVtGYWEhubm5VFVVHR3v2ysE//Fx1O+upCHf1Ox6/yS6e0/UWh8K9mZ4fc3x8NGrCYnSUZLj/b1lZNqD03atCIKgBGYDE4FCYIsgCEtFUdx3unOf7TjcHgK9EXKXwyuL3OFxoFE0v3nqsdUj+DUv5E67C42vCvORyI1Av8b+/OIGJ1FaDdbCcpRqNfqgY1E1FmsmSqVfI9+4PbMG3CK64ZJL54cffsButzN16lTi4qSQyb/7o/WjYzCtLqR+VwXaOO8aQShVKkJi46gqzPNq/MkICPOjpsR62vPIyLQlrWGRDwGyRVE8JIqiA1gI/KMV5j3rcbg8aLxxrbgbQNm8QDvdTjRejPPY7Si88B07HR5UGiXWBknI/f7Wds3u8eCnVOCot6H10zXyjzsdNajVIQjCsdfnqrCBQkAdocNut3Po0CGGDRt2VMRBihjZvXs3ixcvxu12o9Aq0UTpcJY2FdOy8l/I2HMfbnfTmua6wGDqzWYAfqusY/reXMyuk2+crsuq5O4vt1NlOTafn78Gm/nkTR9kZNqb1hDyaKDgL/8vPHKsEYIgTBcEYasgCFsrKipa4badn9ggP2KCvNiMC4htNqsTIEIX0WwyEIAq1Ii2y8mLaoEkYsGRuhMm7jg9IipBwNfgjzG2cYieKLoaiTiAoFGiiTUgKAQ8Hg8AanXTqJ2cnBxycnJQKo/E1ysERHfTyBGzKYOKit9QHOdbiEKhOJoYtN9iZ2l5Lepmon4Olpn5KaMExd/GyR3pZTo6rRG1cry/8ibvOlEUPwQ+BEhNTZXjuYCPpwxufhDANZ97Nez5Uc97NS5y5kyvxo26OgWAH3dL4Yr2v1m0BpUSk9XOqOtubnKtWh2E01mNKIpHhdB/fBz+46UPGj8/PyIjI9myZQu9e/cm8Ej397y8PPbu3UuvXr0AEF0enKU2fHs3TYZqaChDowk5rtBa62qP1jyvdrnwUQj4NPPtp65e2gvw9z324dJgc6HxaSZhS0amnWkNIS8E/lrcIwYoboV5ZToIBh9J2P4Uuj8xqlVUOFyNxPpPdLpkXC4zdnsRvr7HTzy6+OKLmTdvHu+++y6JiYk0NDSQn59PUFAQEyZMAMC6rQzR7sKvb+P6LqIoUlu3FX//pg2jRY+H6sJ8UoaOAOCg1U4Xv+b3GA5VWokK8GlUgqCu3EZA6KmHMMrInAlaw7WyBUgRBCFREAQNcB2wtBXmlekgRB0ps1tcW9/oeIKflnqPh6IGZ5NrAgKlDkVV1X+ccN7o6GimT59Ov379MJvNuN1uxowZwx133IHBYMBRZKHu58NoEgPQpjSOo6+t3YLdXogxZHyTectzD2G3Woju3gunR2S7yUY/w8lDLQF2FdTSJ+ZY+KbD7qKy0IIx9tQ7EcnInAlO2yIXRdElCMI9wDJACXwiiuLe016ZTIchKlCySAuqbY2O99VLx3ebbcT4NPZT63Xd8PNLpqT4f0RH/fOEfuaQkBAuueSSRsdEj4h1axm1S3NQ+KkIvq5bo+tFUeTw4bdQqQIJD7+4yZz71q5CqVKRNHAwf9SYqXO5mWw8eXz9wTIz+dU2bht1LMqm8EANHrdIXM/gk14rI9PetEpmpyiKPwM/t8ZcMh0PnVZFQogfe4oax1P30vviq1CwptrMhaGNLWZBEIiLncKBzH9RWbWS0L/Ek/8d285yEEF0e3CW26jfW4W7yo4m3p/gf3ZHFdDYLVJc/DU1tZvo1u0/KJWN3R71ZhN7Vv1O8uDh+Br8+Wz3IYLVStKCT25Vf7+zCIUAF/Q5lr26f0MJvgY1kcmn1+xZRqatkVP0ZbyiT0wgW3Mbb176KBWMDzHwa2UdL3aNaRLtERl5FfkFn5F54GkChgxAozm+ZVv70yE85iPuGaWANjGAgMkJ+PY2NumfWV29gcyDMwkOGkV01LVN5tq4+Csc9nqGXXkdu8w2fq8y8VhiBFrFib2IlgYXX2zKZ3z3cMIMkhupttxGXkYlAyfHo/Qm1l9Gph2RhVzGK4YnhfDDrmKyyy2khB+zbi8NC+SnijqWV5mY9Df3hUKhpnevN9iy9Up27Z7GgP6fo1I1rWgYft9APHYXgkJAGaBtUnvlT8orlrF37wz8/BLo3fstpFy0YxTu28OOX3+k38QLCYqJ4/ptBwnVqJgabTzpa/t8Qy519U7uHnesmXT694dQqhX0GXfyCpEyMh0B2dSQ8Yq0blLUyIoD5Y2OX2gMJEqrZk7B8XMDDIae9On9FmZzBtu2X4vNlttkjNKgQR3qhyrE97gi7nSayMycSUbGXeh13Rg4YAFqdWN3h6mynB/efInAsAjG3DCFN3PL2G2u5z/J0QSoT2yvFNbYeGdlFpN6hjMgLgiA4qwasreV039CHLqA5qNdZGTaG1nIZbwiKtCX/rGBLNle2CiNXq0Q+L/YUNbXWlhedfyaJKGhk+jXdy52ewnpmy8k59DrOByVzd6zwVFJbu57bNw0nsKiL4iNvZVBgxY2cdFYa2tY8uJMXA4Hlz36NKstDmbllnJVeBD/CDuxf9vtEXl8cQYCAs9cKsWtO+wuVny+H3+jDwMmNZ9cJSPTEZBdKzJec01qLE9+m8GOgloGHrFeAaZEG5lXXMXTWYUMD+yGTtk0gSYkZAxDh/xEdvZL5ObOJi/vI4KChhIUOAxfv3jJwhY9OBxVWG051NakU2fagSi6CA4eTXKXRzAYejWZt7aslCUvPoO5upIrHnuGA7ogpu8+RB+9L690iz1pVuZ7q7JZl13Ji1f0ITrQF1EUWfNVJqYqO5c/OBCNj/z2kOkcCO1RND81NVXcunXrGb+vzOlhtjsZ+dJKhncJYc5NqY3Orasxc/XOHK4MD+KdHnEnFVCrNYfi4q+prFqNzZZznBECBkNPgoPHEBlxOTpdl+OMgcM7tvLzO7MA+MejT5MZGsu0vbnEaDUsGZCMUXNiIf51Twl3frmdf/SL4o1rpf6hO5fns/6bbIZcksjgixJPeK2MTHshCMI2URRT/35cNjlkvMbgo2bKyETeXpHFvmITPaOOVSMcFWTgoYQIZuWW0k3nw73x4SecR6f7//buoKtSaAAAC/BJREFUPTjL8k7j+PeXkzlBEkwghSQEQhAEoqyBloOjI5SCUpjZ3altZ9VWW3e6dasdXFdhZsfZnW137VTdnXbbZVqrrVbRokPbrQq0uLq0ChE5GCScD0ESCSBJSEJOv/0joUHJgZCQJ7e5Pn/xHnjea+4JF3fu53mfu4DCwuUUFi6nubmG+vrDNDVXY8QQH59BUlIesbGJXf79htpa/vfpJ3h3w1qy8vL5/LIVrG6O45+272dKahJPF43vtsTf2n+Cbz23lem56Xz3L4swM/aUVLJx9V7GT8+ieFH+JY2PSFRU5NIrd80Zxy/+dJCHf13Kqr/9zEdm3t/OH8W+ugb+dX/bjaf+LjerxxtOxcUN63TJpDONDfVsffV/2LzmV5ytr2Pm0r/mqqVf4Nv7K3ilqpoFVw7nR1ePJSWu63ujbNxbxdd/XkJuRhI/vWMGSQmxHNhexfondvKpgjTmf/XqCy55FBnsVOTSK2nJ8TywcBIPvbiDF0rK+cKMjtvsxJrxn5PH0uTOv+x7n311DXynMKfHm1V1x92p3L+XnW/8gdLXfk9jfR3jphcz69bbWJeQxt1b9nGmpZV/njCar+VkXXAt+/leLa3g7599h3FXpvCLr80kIyWBPZsrWf/kTjJzU1n8zWuIT9ANsiQ8KnLptVuLc1mz9SgP/6aU6/IzKMjquDY8PsZYOSWfRw5U8PihSjadPsMjE3OZnXHh9eOdqT11klPHjnL80EEq9+3m8M4d1J6oIiY2jqtmzWXaglsoSc/mi4cqKTtTzqz0FL5TmMPk1K5vbOXu/HDDXr6/bjdFOek8+ZUZZKQksH3DEd54fg+jJ6Rz8zemkZCkfw4SJp3slEtScbqBRf/xOiNSEnjxG3NIS77wvuKvnazmH8rKOdLQyI0Zw7hn7Ehmp6d2O2v+zaPfZfdbGwFISc9g9MTJFBR/moQp17KmpolVFSc50tDIxORElo0bxZKs9G6Xbz6sa+TB1Tt4pbSCpdeO5t//qoh4jNdX7Wbn/73PuGsyWXDXFOI0E5cAdHWyU0Uul+zN/Se4/aebuDY3nSfvnEFyJycYz7S08LPyKn505DgnmpoZc0U8N2elMSs9laJhyXzqinhizyviY3vKqKurpzZzFIfjEympruO1kzWUnWnAgOszUrljTCaLMtO6/Q8B2nb8WfbCVk7UNvLAwqv4+vXjqTnZwNqflFJ5oJrrFo5l5pLxxGhNXAKhIpfL4tfb3ue+597hL/IyeOKrMxieeOHMHKCupZVXqk6zuuIUGz+soaG17ecuwYzMhDiSYmJoxTnd3MLp5hbObQh0RYzx6bQUbhgxnKUj0y+4y2JnqmrP8m8v7+JXb5dTkJXC47dOZ+qY4ezeVMnrz5bhDjfdPpkJ143st3EQGQgqcrlsfrfjGPc+9w55I5L579uKmTCy+/Xws62tbK+pZ9eZeg7WN1LV2ExDaysGpMXFkh4fx/ikK5iUmsjE5ESSLvJkaUNTC7986zCPrd9NQ1MLd80dz73zCmmtb+aNVXvYt+WDP1+ZMjxTm0VIeFTkcln9ad8J7vnlFhqaWlh+y2S+NCNvwJYszja3sGrzEf5rwz4qqhu4vjCTh5dMoSArlT0llbz2TBktTa3MWJzP9AVjtZQiwdIXguSymlVwJb/91lyWPb+NFS+9y0tbjvLgokkU51++TRnKT9Xx7KbDrNp8hKraRmbkZ/Dordcwu6DjboeJqfFk5Q3jxi9fRfqonncJEgmRZuTSr9ydF94u55FXyqiqPcvcCZncPmssN00aSVwfric/p6r2LGtLK3n53WNs3Nt24615k0fxldn5zC7ofCPmzvYUFQmRllZkQNU1NvPUHw/x1B8PUlHdwIiUBOZPHsmcCZlMz80gd0RSj+Xq7lRWn2XnsdNsOnCKTQdOsPXIh7Q6jL0ymc8XjeaLM3PJydBMW4YGFblEormlld/v+oDf7TjGH977gJqzzQCkJ8eTk5FE9vAkRqTEExtjxJhR39TCyTONnKht5GDVmT+/Pz7WmDYmjbmFWSyams2k7GGaZcuQozVyiURcbAyfm5LN56Zk09TSSllFDdvKP+Tdo9UcO11P+ak6dhxtpNWhtdVJjI9lREoCGSkJTM9Lp3BkKhNHDaMoJ50kfWlHpFMqchkw8bExTB2TxtQx3e9oLyK9ox2CREQCpyIXEQmcilxEJHAqchGRwKnIRUQCpyIXEQmcilxEJHAqchGRwKnIRUQCpyIXEQmcilxEJHAqchGRwKnIRUQC16ciN7PvmdkuM9tuZi+ZWXp/BRMRkYvT1xn5OmCquxcBu4GH+h5JRER6o09F7u5r3b25/eGbQE7fI4mISG/05xr5ncDL/Xg8ERG5CD3uEGRm64HsTl5a4e5r2t+zAmgGnunmOHcDdwPk5eVdUlgREblQj0Xu7vO7e93M7gAWA/O8m52c3X0lsBLaNl/uZU4REelCn/bsNLOFwD8CN7h7Xf9EEhGR3ujrGvkPgGHAOjPbamY/7odMIiLSC32akbv7hP4KIiIil0bf7BQRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHAqchFRAKnIhcRCZyKXEQkcCpyEZHA9UuRm9n9ZuZmltkfxxMRkYvX5yI3s1zgs8DhvscREZHe6o8Z+WPAA4D3w7FERKSX+lTkZrYEOOru2y7ivXebWYmZlRw/frwvHysiIueJ6+kNZrYeyO7kpRXAcmDBxXyQu68EVgIUFxdr9i4i0k96LHJ3n9/Z82Y2DRgHbDMzgBxgi5nNdPeKfk0pIiJd6rHIu+LuO4CR5x6b2UGg2N2r+iGXiIhcJF1HLiISuEuekX+cu+f317FEROTiaUYuIhI4FbmISOBU5CIigVORi4gETkUuIhI4FbmISOBU5CIigVORi4gETkUuIhI4FbmISOBU5CIigVORi4gETkUuIhI4FbmISOBU5CIigVORi4gEztwHfh9kMzsOHBrwD/6oTEDb0rXRWHTQWHTQWHQYLGMx1t2zPv5kJEU+GJhZibsXR51jMNBYdNBYdNBYdBjsY6GlFRGRwKnIRUQCN5SLfGXUAQYRjUUHjUUHjUWHQT0WQ3aNXETkk2Ioz8hFRD4RVOQiIoFTkQNmdr+ZuZllRp0lKmb2PTPbZWbbzewlM0uPOtNAM7OFZlZmZnvN7MGo80TFzHLNbIOZvWdmpWZ2b9SZomZmsWb2jpn9NuosnRnyRW5mucBngcNRZ4nYOmCquxcBu4GHIs4zoMwsFvghsAi4GviSmV0dbarINAPL3H0y8Bngm0N4LM65F3gv6hBdGfJFDjwGPAAM6bO+7r7W3ZvbH74J5ESZJwIzgb3uvt/dG4HngKURZ4qEux9z9y3tf66hrcDGRJsqOmaWA9wC/CTqLF0Z0kVuZkuAo+6+Leosg8ydwMtRhxhgY4Aj5z0uZwiX1zlmlg9MB96KNkmkHqdtstcadZCuxEUd4HIzs/VAdicvrQCWAwsGNlF0uhsLd1/T/p4VtP1q/cxAZhsErJPnhvRvaWaWCqwG7nP36qjzRMHMFgMfuPvbZnZj1Hm68okvcnef39nzZjYNGAdsMzNoW0rYYmYz3b1iACMOmK7G4hwzuwNYDMzzofcFg3Ig97zHOcD7EWWJnJnF01biz7j7i1HnidAcYImZ3QwkAsPN7Gl3/5uIc32EvhDUzswOAsXuPhjucDbgzGwh8Chwg7sfjzrPQDOzONpO8s4DjgKbgS+7e2mkwSJgbTObp4CT7n5f1HkGi/YZ+f3uvjjqLB83pNfI5SN+AAwD1pnZVjP7cdSBBlL7id57gFdpO7n3/FAs8XZzgNuAm9p/Fra2z0hlkNKMXEQkcJqRi4gETkUuIhI4FbmISOBU5CIigVORi4gETkUuIhI4FbmISOD+H/kpfWureXEfAAAAAElFTkSuQmCC\n",
"text/plain": "<Figure size 432x288 with 1 Axes>"
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "from scipy.integrate import ode",
"execution_count": 6,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "def f(t, X):\n x, y = X\n return np.array([y, x ** 2 - 1])",
"execution_count": 44,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "solver = ode(f)",
"execution_count": 61,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "T = np.linspace(0, 5, retstep=True)\ndef get_trajectory(solver, initial_value, T, max_norm):\n solver.set_initial_value(initial_value)\n traj = [initial_value]\n for t in T:\n if t == 0:\n continue\n x = solver.integrate(t)\n if np.linalg.norm(x) > max_norm:\n break\n traj.append(x)\n return np.array(traj)",
"execution_count": 87,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "for x0 in np.linspace(-3, 0):\n traj = get_trajectory(solver, np.array([x0, 0]), np.linspace(0, 20, 300), 10)\n plt.plot(traj[:, 0], traj[:, 1])\n traj = get_trajectory(solver, np.array([x0, 0]), np.linspace(0, -20, 300), 10)\n plt.plot(traj[:, 0], traj[:, 1])",
"execution_count": 95,
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd5gcxYH38W91mpw2B2lXYbXKEkJCSCIjcrbNOeF0L7bPr304h/PZxoE723dO57P92jiCfcYBk0EgIYIIAqGMctrV5jw5d6j3j1lkgsCcUQDUn+dpTaunZzrs7q9qqqp7hJQSl8vlcr35Kcd7B1wul8t1bLiB73K5XCcIN/BdLpfrBOEGvsvlcp0g3MB3uVyuE4R2vHfgldTU1MhJkyYd791wuVyuN4yNGzeOSilrD/fc6zrwJ02axIYNG473brhcLtcbhhCi6+Wec5t0XC6X6wThBr7L5XKdINzAd7lcrhOEG/gul8t1gnAD3+VyuU4QbuC7XC7XCcINfJfL5TpBuIHvcrlcx1kyX+Zrd+8gW7KO6nZe1xdeuVwu15vd2v2jfPrPWxnLlThrei3L2iIIBIZqHPFtuTV8l8vlOg5Kls03V+zi3b9ch9+jcsdHT2NSQ4533/duvrfhe0dlm24N3+VyuY6x/cMZPv6HLewcSHPNqS18+dJZPNR7P9949Bt4VA+nNZ92VLbrBr7L5XIdI1JK/ufpLv7tvl0EPBq/fN8iTmsP8x/P3MBt+25jQd0CPur9V3KrHeT7JUKII7p9N/BdLpfrGBjNlvj8X57l4d3DnNVey3f+YR5Zu5933/cR9if388GZH2L27uXc+kwPp0+pxio76B71iO6DG/gul8t1lD2ye5jP/WUr6aLF1y6fxfuXTeKejnv4t6f/DZ/m44cLf0r/XQbfG+thV9DmvPMbj3jYgxv4LpfLddQUzUrH7G+f6mJGQ4jff3AJE6tVvvrgDeze0sekqUv44sTrWPXrPv6iZRg1HD534XQunN1wVPbHDXyXy+U6Cnb0p/jkH7ewbzjLtadP5nMXTqc328l1N/2IGVvPxa8tpmpAZb3WzW+MEl6vxm+vWcQZ0w773SVHhBv4LpfLdQQ5juRXT3TynZV7iPp1fnftYk5vq+GOPXfy4K1bWGi2syKoc3lR46Ah+aYoMqcmwI0fPJXmqO+o7psb+C6Xy3WEDKaKfObWLTy5f4wLZ9fzrbfOw2tYfOX+G1AfbWBqc5774g7Xlrz8j7fMJmzetaCZr71tLh7tyLfZv5gb+C6Xy3UE3L9tgC/esQ2z7PCz09o4//yp7M8f4Ls3/5z5XXPY07qJoeFm3p+bzbc8RVICvvvWuVx9Sssx28fXHPhCiOnAn563aApwvZTyv563ztnAXUDn+KLbpZTfeK3bdrlcruMtV7L4+j07uG19L/+gBbhKyyPSX+U3K+bx7FbBIr2eu+t3c0HvPCw5gX/VCtQHDO64djGzmyLHdF9fc+BLKfcAJwEIIVSgD7jjMKs+LqW87LVuz+VyuV4vtvQkue7WB5D2Gr4u38GcyDZ6Z/6Ge0eqmLi5kdZJW/lTxsN1yRnc1LCN7ekMZ9Qv44fvWUjErx/z/T3STTrLgQNSypf91nSXy+V6o7MdyY8f3sNPN93MsrLDNYkrUGf9ietzXkI7Tuei5Mk8M/UZkiPVfMwT49+rn6IYeZrW6lZ+9a5PoqrH5zZmRzrw3wn84WWeWyqE2Ar0A5+VUu443EpCiA8DHwZoaTl2bVsul8v1anSMZPnkbasYLdzBR4Yu4qxYmbXzf8tNgy1cnJjAWbUZvle/jTP6p9Na38/1yuMogf0srDqdn170veMW9gBCSnlk3kgIg0qYz5ZSDr3ouTDgSCmzQohLgB9KKaf9rfdctGiR3LBhwxHZP5fL5XotTNvhxsf28ZMNv2FhucA1o5cQnb6KGx2HzozDR8ZOx5r2KDcmPHwiM4t767ewwbsVRcvwsXmf4p8WvO+I3xvncIQQG6WUiw733JGs4V8MbHpx2ANIKdPPm18hhPh/QogaKeXoEdy+y+VyHRXb+1J86vZVZOw7eP/w+Zzv97F74U18rXcWs8olvmG3cc/UhzgwGuSLShPfrltHIrgevwhz44W/5eSGecf7EIAjG/jv4mWac4QQDcCQlFIKIRZTuQ//2BHctsvlch1xRdPm+w/u5ubtv2W+k+ZjAx+gpW0tv/MO8mhfM+9NNLOwuZtvyy0sGGzhspoxvshaZHAP0wILuOnyHxHxHNuROK/kiAS+EMIPnA/80/OWfQRASvkz4Grg/wohLKAAvFMeqbYkl8vlOgrWdYzx2bsepKjewbvGzuVitYGhRb/jM0NtGAW4Ib6URPsqvhD38M+ZGWxs2sL31b0oepx3Tv0wXzztYyji9fUdU0esDf9ocNvwXS7XsZYpmnxrxU5u3fd75ipjvL337cxs2cFdVVv580g95+WDvEtp5taGdexKePmEnMQPQtvp8W9Gx8cPzv0+Z7cuOW77f6za8F0ul+sN7aFdQ/zrvauxvHdwdfJMLrPPobjgVr6crGcoKfn0yElMm7Sdr1tbmds/iffWxfkXuQ4zuJNGfQa3XPlTagM1x/swXpYb+C6X64Q3li3xtbu3s7L3T8zQBrl6/7Wc1NDLoxN+xy/7JzNbJvhM8i10zbqXfx0N86HkXHa1rucG2YHiHeKCpnfxn8u/gKoc/fvhvBZu4LtcrhOWlJK7tvTztfsfRgnfzmWpZVxeOB/vnHv4Tklj+4ifaxItnBeeyK9aVzPcX8MXxQR+3LiBbm0TmlD4t6X/xeXty4/3obwqbuC7XK4TUn+ywL/esYW1o7czzdfLW/b9H06pSrNl1q/4Ud90mujnG0NXoUx7hM9kPZzbO4PTGwb4EuuxQpup1Sbzuyt+xoRQ0/E+lFfNDXyXy3VCcRzJ79d18e2HHsEfvY8rE0u5pHAhkZkP8kslw5qBai5NhblKb2PVzHt5cDjMJzKzWTdxHf+u7Efx9XBu41V8d/n16Oqxvx/Oa+EGvsvlOmEcGMnyudvWsSN3Kydpgqv2foi5dYPsnvtjru+bj49RPt97BRPanuEGa5CG/nquVxv5bsMz9HjWoak215/6bd4249LjfSh/FzfwXS7Xm16maPKjh/Zx87Y7aQxs5AMDV3C2aESbfze/psSq3imck9G42r6I3fPv5gvD1bwnOYNIUwefcTZiR9ZRozfzm0v/H1Oik4/34fzd3MB3uVxvWo4juX1zH99avQbTfydnFqdzef/HmNq8h02tP+HGrqXo2j4+0X8ac1t3c6PxOANDQb6emcuKietZLQ6g+rs5s/5CvnfeDfi0o/sVhEebG/gul+tNaUtPkuvv3sTu4u20G3EuOXANp3p1igt/x4+LIR7pms1ZxRzvyF/O2Lx7+ORYDWfEG3l3ROWbdc8w4n0aQ5F8Zcm3eNv0N8dXebiB73K53lSGM0X+8/7d3LlvJTVVj3Hl6PlckJ9P06SNPNH4BL/sPB2vZxsfHz6NxfUD/Ln5PtaMerkuPpvspA18uhyH8FaaPG386pL/ZmJ44vE+pCPGDXyXy/WmULYcbl57kB+ueQpZdQeL1QYu3vkJTormGVl8I/+RaGVdz2SW2YO8Z+Dt2LPu4V/SYapH/fyHOZvftqxlndiLEhrlqsnv5frTP4WuvLFG4fwtbuC7XK43vEf2DPONe7fQJ1cwMdbBeR1v5zSnidj0x1ldtY2bDy7GH1zPx4au4syqAqvn/JFbRqp4e6KFeU29fMneQMr7DF4lwPfP+RlntSw73od0VLiB73K53rAOjua44d6dPNr7GLH6lZzbt4zzuz7JrNoROtt/xH8PzmX7YA2LZT/v67oW34z7+WZJkB4McUN6MdsnP8Lny3FEZB/twUX86tLvE/PGjvdhHTVu4LtcrjecbMnixw/v59dPb8Sou4e5Pp0LdlzHQtWPZ/b93Bbo4w+di6n1b+af+9/CGSGHp+b9nl8NNXFOpoZLIvDDpifYy04Uf54Pzvo0H1/0gWPyjVTHkxv4LpfrDcNxJHdu6eNbD2wm5VlFS203Zxy8imXlyUxs2seuyXfz8+7F9BRGOaOU4+rhD+Gfvor/tBz6BmN8dnQhxbbH+VQpjxlYT0Rt4MYLf8nculnH+9COCTfwXS7X656Ukod3D/OdVTs4UHyQ2qpNXNKznLPSb2VGVZLU/F9ws6VzZ8ciWo09XNf1Nk6tSfHUSb/jF4OtLCvo/F8xgd9PWcMa2Yka6eHM+kv47vKv4tf9x/vwjhk38F0u1+va2gOjfGflLrYlHyNa+yhnphdy1vbPMi8gEHNWsKZqB/9zYDk57xYuzHi5rPhePLPv4z9KBt3DIa4bnkvdpO183tpG2rMOQxV8dem3eUv7G/P2CK+FG/gul+t1aXN3gu+u2sPT/U8TaljFydoEztr2SeapAWomb2XflHv4ffdpbDpwMnPV3VzZ8S7mN/bwxIyb+VX/ZE4uO3yrtIiH2x/nO6UkatUWWnzT+cUlP6Q52Hy8D++4cAPf5XK9ruzsT/P9B/fwcMcWgo0raY/AWXvfx0lmA1Nrhxic9kt+na/h3u1XEok8ybsSZ3O+cxrOgjv4di5Cx3CID43OYNaEg3ybdXQpz6KGE7xz2j/yhSUfR1NO3Ng7Ul9ifhDIADZgvfj7FEWl6/uHwCVAHviAlHLTkdi2y+V6czgwkuUHD+7lvp07CTasprmhjzO63srJmRnMiZjkZ93K/cEufr/vYlKBjSxU9/G2A9cybeIuHqq/id/0T2O2k+HbmeXsnvEw12VtZGQdUa2WH513Myc3LDjeh3jcHcmi7hwp5ejLPHcxMG18OhX46fijy+U6wfUm8vxw9T5u27IXX92j1E7ayim9F7J45L3M9at42p9ie+sqbjlwLttGIkyIrOIdnW/nLCNMauHvuSFVy4GRAO9LTGJp/Qg/jqxhi9WBGh1gefPlfOusL59QHbOv5Fh9trkS+K2UUgJPCyGiQohGKeXAMdq+y+V6nRnOFPnJw/u5Zf0+tNhT1LY9zfSBRZy2/yvM1Dw0NnfRO+OP3DU2mfu3XUmg6iEuSC3iqv5P0Dh5PauqN3BTfzszZYpvJi4jMWMFH0/pFH1P4Vf9fOvM/+b8Secc78N8XTlSgS+BVUIICdwopfz5i55vBnqe9//e8WUvCXwhxIeBDwO0tLQcod1zuVyvFx0jWX7xeCe3bdmLiDxBbMpGpg+ezJLNX2SG4qO1Jkl88i3craf4467LyQafZr6xk8v2fJQF4SJ9p/yK68da6R2J8N74VM6sSfDbKQ/wSGEYNdbJ/Nhp/OiCb1LlrTreh/q6c6QC/zQpZb8Qog54UAixW0r52POeP9zla/JwbzReWPwcYNGiRYddx+VyvbFIKdnYleDGxzp4aO9+PNVPEp60iZmDizh16xdpx8fUaJ7UpD+zrmY7f9h7KXudMVrDq3nrwbdwhmwiMGMVf/KOcEfvVBbYWb6RvILs9BV8MuMho2/C8Ao+u+irXDPrbW/6K2b/Xkck8KWU/eOPw0KIO4DFwPMDvxd4/j1GJwD9R2LbLpfr9ct2JKt2DPLzxzvY0n+QUP3jVE3ZwszBJZx68CtMlV7aIxaZifeytWktd3aew1N9VxGtfpiL+k/nvJ7rmDphL5tbfsIvehZh5kf58Gg7y2qS3NJ2PyuLSdSqHUwOzOanF36XCaEJx/uQX9dec+ALIQKAIqXMjM9fAHzjRavdDfyzEOKPVDprU277vcv15lUo2/xlYw+/fKKT7nQ3Vc1PUDXlWWYNLePU/V9jkvQwMyzINz7Mnkmrubf3FFZv+BBG1YOcqmU4b9d1zA8XyZxyEz/IxXiyZwrLSinenXormfaVfDLrJaltRfPk+dCc6/jogmtRFfV4H/br3pGo4dcDd4x/hNKAW6SUDwghPgIgpfwZsILKkMz9VIZl/uMR2K7L5XqdGcuW+O1TXfzu6S6SVjd1E58kVv0sc4bPZNHg12mxPcwKq5TrnqKz7R5WDc1gxYaPYIafoC12P8s7r2axqCEy82Eeiu7mt12LCHu384mhKzilbog/tt/HvYUCWvUjNHhb+cl5v2ZG9YzjfdhvGKIycOb1adGiRXLDhg3HezdcLtffsL0vxS3PdHPbxh4sYz/NrRvIWbuYN3oWCwfPp9EymBvRcaKbGZ5xB2uSddzZcQkZ/2bqgjs4resyzsnPo7Wxg71TbufXfYvYL/s5O9PA2/OnkGy/j++mIsRDa1H0FG9vv4bPn/pJPKrneB/6644QYuOLr4V6zol7yZnL5XpNUgWTu7f08acNPWwfGMYb20zV9PUo+RIndV/AjOH30YzKzIgBoV2MzriddUWV2569hiHlINHGmzlz4FQu6P4Cs8MFMgtv4iZb4e6OU5jo2cPnev+Bec37+VPz3dyTtVCrV1BrNPODc/+bBfXuRVR/DzfwXS7Xqyal5OmOOH/e0MOKbQOYaj91zRupmr6emlQzpx94B/WDrbR6VdojGuXoFsba72e7necve67mgJkhXPcnTh2dwbnbPsc83YO3/TGeqn+Gmw6cSd63jcszAa5MvZ2x+XfyqXgdo+xGjaT4h7b38PlTP4FX8x7v0/CG5Qa+y+X6m4bSRf6ysZdbN/RwcCxNqHo3jTM2kDD30ZpYzNKDXyaYCDEtqDGpGnI1zzAw4362F2DF/ivYnlEJNjzA/HQtZ+74Z+YTpnHCHvZOuZvf953Mps5ZLKCLt3W8hylT1vPHwL2sSArUqvupNSbwX8t/xEl1Jx3v0/CG5wa+y+U6rJJls2bPCH/e0MMje0ZwlCSTJj9LY9NarEKZuQMXM3PgQwSKGrNiBvVVJumGx+huX8mzmQArdlzNroxBsH4VMxQ4fe81zDUbmFY7xFDb77gxV82KPedQH9jEh4Yv5iyvl86Ft/CJkSaS7EGNpLi67T18wa3VHzFu4LtcrkPyZYs1e0a4f/sgD+8eJmtmqardw9S5Oxkqb0fLtHDp6Pup6plMVArm1XoJ+XKkJq6kc+oqNidquf/ZD7AvaxNqfJTJ3jTLuq7i5PxUZkQL5Kbcwj3eEf7QcQ7lwHqWl7NcOXwtVe2P8Es1yZqEilZzH3WeCfzgXLdWf6S5ge9yneBSeZOHdg/xwPZB1uwdoWSXiFYfoHHaTsacLRhFH+19y7li5Br0lIcWv8q0Gh3dSpFseYCOlgfZMNrCA5s+RmchT7jxQZpCIyzruZJFqbnMCjo4sx5kY8M6/mf/BXTgZ5ZnM2/pfAcn1cTZsPAmvjQ0iXxsC5ovzbva38enT7nOrdUfBW7gu1wnoJFMiQd3DvHAjkHW7h/Fcmxqa3ton72LIXs9ZbNEXWIJlyW/hNYfJaoIZtZ4qNUcSoEOUu1rSMSeYv3QLB5Y/1l6ykmijffQFImzpP8SFowtYJZPxTtlIwen3sOtfQt5fNdy6iNruab7Ys6Xk7Dn3sN3Sh42JAR63Z00eCfy/XN/wvza+cf79LxpuYHvcp0AsiWL9QfjPH1gjLUHxtjen0JKi+aGQU46qYNB+2lSpTh2up2rcx8l3D0BUYSpUYMpTRpaIUe29nF629cwKrtZO3AGj+65gWG7n2jTbUww0yzpvYzZo3OY5deJNnUyMP3P3Jat4a5t/4ASeZizLI0r9l7H5NatrK69iZsGpmBXP4yhFfnArA/ysZM/4o6rP8rcwD+BSCmxHEnJciiZduXRcrAdB0UIFCFQFYEQvGBeHZ/36ioeTXFvTPUGUCjbbOxKsPbAKE91jPFsbwrbkRieNK0Te5h90l6GzW2krTz+sUlcYL6TxuF2yqOCao/C7Ho/0XwZUx8mM+0JknWP0p31sqbzcp7smUPZs4/q5t8zsZhnaddltI3OYmZApa6pj9Fpt3O7Lblt15Wk/BuZ7dnIxfvfz6JwkaFTfs1X4hPZk7bQ6u+hLTyL/zz7Btpj7cf7lJ0Q3Ctt34AcRzKSLTGaLTGWLRPPlRnLlRnLlojnyoxmy8RzJRJ5k3zZGg94h5Jl47zGH7ciwKer+AwNv6GOz6v4x6egRyPs04n4dMLe8UefTtinHVoW9esEPZpbcBwhUkp6EwV2DqTZ3pdiXUeczT0JTFuiKQ5tLSNEazpIym305TtQHIWZ5UUsyp9NaKAJKy2IqDC93k8dEiVfJl+/g/TMx0nqG9g2OpvHBq5iy1AEb9VmovXr8GdUlg5cQVu8nel+jfrIIGNT72SDPsat+6+gV+mlxb+Xc7su53RnAqH2B7ndM8TtYzUo1Y+gqwqfPPnjvGfWu9174Bxhr3SlrRv4r1NSVkL94GieztEsHaM5Do7m6BzN0TWWp2Q5L3mNqkhiIZtIoEQwUMDrKaFrFkIxURQLhFmZlDISCykqj0IIBILKXayV8RtXj8+PL1fQEdIDjo7t6Di2hm3rWJaGaWmYpka5bFAoecjkNTIFXrFw0VVBzG9QFahMsYBB1Yv+XxMwqAoaVAc8xPw6mqocjVP9hmLZDgdGcuzoT7GjP83O/jQ7B9KkCiYAilJm6sQE9bVD2MZBOrPPkjOzeB0/Z3AR7cmFiO4wVsEhaijMaPRTazuIbJlSuIvczI2komtJlXM8PbScR3rPYSCfJ9awHiW8jliqjjOH38LEsRZmBDTqAqOMtd3NjmAHtx24nF1Fk5qqJzm173QuSC9lcuN+NrY+wM975zMWfgrV18PC2qV888yv0RRsOs5n883JDfzXOSklHaM5NnYl2NydYHtfms7RHNmSdWgdXYXmGpPaqjSBQBzFSGCLNGWZJmcnSJcTpEpJLGm9wpYACX4RwC+CBAiiKwYgkUgcHBzp4Dz3P+kgpaSkFCjIPAUnf/hvNjgMj+ohqIfwa0G8ahCPEkQXfnSCKE4Q2w5glf2USl7yBS+ZnIdETiOVf2lBBiAERH061UEP1QGD6ucKgoBBdeCvBcWhAsNvYGhvzAKibDn0JQv0xPP0JPL0xCvzXfEce4eylMcLe48mmNpUorp6AOHtJm7tpTfXgS1tFEdlvlzCnPJiahOtlPoVbFNSFdCY2einyrRRUmXK/mHyM7eQqnmCvNPDnsRMNo5eylO9Eylr3TRMfIaCspXJ8fmcNnY5tckoM4M6Nb4E8an30lG1lTs7LubpRDWh+lXMjbdx0eDFzA0XiE+7g5tzVawtWBjVa/BrIb6y9F+4dMql7qe7o8gN/NeZXMlia2+STV0JNnUn2dSdIJmv1NBCXkH7xByRcBzNM0pZGSJp9tOf76ZgFQ69h1fx0qhMoE40Uy3riDo1BMwIPjOIXvSjFHUoq0gTHBPsssQ2Hazy4QP11RAKaIaKoglUXSB0EJpE6BLhcXAMC0srUdaKlLUCeTVLTqTJihQpESfOCCNykIyVPvz7IwgbYSKeGCEthleNYBBBkUGkFcQy/RSLAXJ5L6mch3hWkMq/fAEX8mrE/AYRX6UZKezTiY43Nz23LOLTCXp0fIZKwKPi1zX8HpWAoeHVj0x/hZSSgmmTyJskcmWSeZNEvkwyX64sy5dJ5Mr0p4r0xvMMpIs8/89SVwWNMZvqWJJoNIHuHaHAAF2ZvSRKicqxqmFOUc6kLT+f0Gg9xT6wTYkmYEq9j+aQTrhoQ6aMZWQozHqWdMNTZJ3tDOZq2Rh/C0/2zGIoZxGq3kW0/hkKpQEWjJ7D/OGziZYM5kQMolqCxNSVdNWtY8XBs3l0aAZazUraygbnd7+NRWoYX9tqVoU6+Z/eOdjVq1GMMS5svZQvLfkCMW/sNZ9P1ytzA/84k1KyezDDA9sHeWj3EDv704eaO6bUK0xqHMUT6iYl97E/tfNQsAsErcYUponZNFtTiBbq8GRCWHGNfNxEHqbNxGMoREM6Ib+O11DQVQVdAU0RaIpApdJTryBREC+ssctD/4zvN6AIHMAGbPncVOn8tRyJaUuKlkOx7JAvWuQLNoWi/Yrnw+PX0AMKqg/w2TieMqZRoGjkyGpJUmqcMYYYYYBBu5+MefgCwqN6iHmriBgxgloUrxJBJ4xwAjhWANP0USp5KZa85AuVpqZ0wSZVMLFeRWeGEOAf76NQFVHpvFbFoU7syqSgKmDZkrLtYNoOZeuvkzm+/JUEPRoRv6AmUqQ6UiQYKKAbKYpikLjZQ1+261CwAwSUIDOVk2izZtNYnIwnGSHTa2GbDgiY1OinJeohYtooY0Wk42BGBym27yFbtZmMs52cabA1cRFPDSxj55CB5uthYusOMup6/NkQy8YuY+LgXOpQmF3lwasOkGi7n+7YFu7vOos1AwsQVY/Soo5xRvdVLCm30jxhBztaV/CrnsUc1PehR7bS4JvAN06/nqVNS//m+XYdGW7gHweOI9nck2TVjkEe2DFI11geIWB+i0JTYzeO5yADxd10pg8gkahCZY53AfPtJTRkpqDE/RRHbYrZv9ZgNU2hodZLbcQgbKh4BeiORDUdRNGCvIl8hRq8RCJVE+mtTEIFRKU5pxL8kucHvnA0MDVESYWygrB0BK+imUQB4dHAUMFQcHQVR1MwBZSBkiMpmA75kkO2aJHOm+QyFs5hQlhRBN6QjhFUUAIO0mtieYsUjRw5PU1ajRNXhhmRAwzbA4yVxrCcw9f6VaES8USIeqL4tSAe1Y8h/OiKHw0fGj4U/AjHi+MYWJaCaalYtoqQGlJq4FQeHSmwbQdbViZFgK6JSuGqClTVBqWEFCVQiiBKaFoJRS2DUsIWedLmKMnSGKOFkRcE+nMiRpTpntlMYTp1xRYCqSqcUYPMUBnHqpwrj64wsd5HY9ggqoA6VkQWLBxhUp7WRWHiDlKeZyhZfWTLfvZkLmDr6Cms7w1jkqZxwg70yEaSpT7a0gtYlryUwEA1rV6F6WEDxzhIon0FPaGdPHBwOWsGTkZEH6fe18mpXRezLDeX6dVxRtr/zJ+zNaxOhtBqVqOpNh+a9yE+OO9ad6jlMeYG/jFiO5KnO8Z4YPsgK3cMMpwpoauChVMFDY37GLY3sH1sCxJJUAux2HM608sLiCWaKPepZONlAHSPysRGPzVBnYiu4LUdtJyJkyxVqtjjpNfEqc5hhzM4oTS2L4XpSWLpcUwlgUMeWxYqk5PDtvPA39+kAyCEjqJ4UGfKxOIAACAASURBVIQHVQmgiSAqQVTpR3UCKLYf1fKjmD7UUgClGETJBVAyfkTKiygefiSw0BWEXwefhmMoWKpCWUBRQsFyyJVsMgWbVKZMLmPi2C9fOHgCKmoAhNfB8ZWxjBJFPUdey5BVUqRFnDRJUsRJO0myVpZsOUvRLr6mc/NqeVUvIT1EvaeRBrWZWhqJmbUEijGMXAAyOuUE5BMm9vM6530hnQkNfhqCOmEp0bNl5FgRJEgczOYhzNYu8rFdJO31OE6eRKme3ZnL2TQ8i619GjYmNbWdVDU8y7C1mVC+mjMylzNxYA6egkJ7WKdFFxSDe0nOvJ9eYx8ruy7mkb6TELEnqA7uYGHveSyLL2FOyMaafB+rg33ccvBUSlWrUX19nFy7mG+cfj2t4dZjcj5dL+QG/lGWLVncuqGHm9YepGssj09XWTzNobp+D73ldeyK7wBgemAm58jLqR2cSqrDolyo1ESDYZ0pzUHq/RqBkgXD+UM1dSkk1OcwG0cwY0OUfL0U1W4K5kEs+4VNHVJC0YmRcyaTsxspOxFKjp+iFaBoeylaXgqWQcEysBwViUBKgSMrr3UQlUcJmgK66mAoNppqoys2mmKhKya6YuLTSvi0HB41g1dJ4VGSGCKOThyPMoqmHL5gUdUAuhpDV6LoVKE5MTQzilaOoObDKNkgajqISAQgJw77VffCp6EEdaS3UjiYivLXTw6WQ/a5Tw5Zk0LWxHyFJiZFE3j9Oh6/huZVUXQJ4/0S6E5l0hyk5mBjYWNjY+HgoCgCoQiEUml+wxZgKUgbFEdFkzqq1FFtFVlSsAtQzjuUchbF7AvD/Dkev0a02ktNWCdqqAQVgce0UQsWTqII4zV7GXCwpgxSrN9P3reLjLUN287iSMFg6WT2Zs5lfX8Lu4ckiBITm7uI1O5i0NyCXbaZnz6DhYlzUYaCVOuCObVewsUSudrNpGY8Qo/sY3XPxTzSOw8l+iSh6Cbm9Z/G6cPnMdtQCbSs49nmNdzceSad2kGM2FOEjShfXvpFLpp0kdspexy5gX+U9MTz3Lz2IH9a30OmZLGg1c+s9r3syq5iX3IvAPPCCzjTvoTq/kmM7i1iWw6BiM7MSWFqDQVv1sQZzldq7gLkhDzl1i7ysd0UlA7y5Q5sO3domwVnImPWSQwVp5EoVhEvBBjNG4zkFIYzNkXz5WrwNppeIuAz8RkWugaKIhHCQQgJQqJQmRdCwbFVTFtgWiqmpVC2FExToWQqOLYHeOWx02GvSsyvEvVJIl6LsKdE2FMgqGcJ60lCxghBdYCA0oMm+yr78CKaFsHQajGUanRZg27F0Mox1EIENRc+VDDItIDDhCcKKH4dxa8jPSpyvHnJUgSmEJRlpX29aEqKZZuS6VA0HUplG7PsYJVszFJl/nD9Ja9IgKopaLqCoil4vCq+gE7AP94prCl4FTAA3XbQHYlSsnHSZZx0+QXHIGodzIYhylX9lAI9FNQOssVdSGniSEHcXkJHdik7R1vZ0qeRKdkItUhbaxf+2E56i5sxbYvphQUsyV+Av6sB3YT2Kg8tHgXMBOlJj5Oc9AgHMgarey/l6cGpaLG1RGLrmTFwCqcNXcAs1Utdw1662u7g1uGZrElG8NQ9AFqOd7S/nY8v/DhhI/y/O0+uI84N/CNISsmGrgS/fqKTlTsGUYRg+RwPVU3reWzgHtLlNHMi81huXUW0t4WRfXkcWxKKGsyeGqVBAXrSyKKNVIDJWUotneTDu8k4WymWegFwRIwR8zSGiu30ZRvpTgU5MApjOfPQvhiaRXW0QDSUIxDIoRsp0FKYIkHRSVFychTtPHkzR9kpHdHz4FE9+LQAPjWAofjRhQ/tecMuHcuPafoplXxk8x7SOQ+prAfTemmTjqEp1AY1qgOCar9Fta9EzJchaiQIG0OE1AGC6kEUpw8pzZe8XtMiePQ6dLUanWp0O4ZmRlCLYdRiCCUXQs0GEGkfMm/j5K3DfnI4RFT6IBRDAU1BaApoAqE+71ERlUJAykrntqx08o6/HFm2K1PJxilV5l+uNU14NUQEnFgaK5LAjiQx/UMUlE5y5n5KpYFD69qimlFzMX2F+ewancDmPoVE3gIcJtTHaWzopWzspju3Hcd2mFU8hcWF5fh7GrALkiafyoxqL8FcmWKwi/SsNSTCT7B9bAoP9b2VbaNBArVr8Uc2Mn1gIUsGL6Jd+JhUM8Jo25950NL4S/cSrKoHUf1dzKyaw1eXfpnZNbNf9e+O6+hyA/8I2dqT5Ov37GBTd5KIT+fCk01K/kd5rG81Dg7n11zM6YnLGd1kU8pZhGMeZk2NVEK+uxLyeATO/GEyTc+QsB6jVB4CQCp1DJjncSA1l23D1WzrLx8ab+332ExsyFAVHUPzDlIQfQwVD5IsxV+wf4qj0kQrE8RkqmQNfieMzwngsfzolgfN8qCWdYStIqRAOoAz3mziVP4voFKrVCWoEqlIUCRScZCqg9QsTK2EqZYoqXkKSp68yJAX2crQSzHMsDNAQeQOO2Y/qIeIeqoJqlV4RARVRpBWGLMcpFAIksr4GU16yRymOT3oUakL6dQGBTUBk2pfgZg3TdSIEzYGCap9+OjCskaQsvzSN0DBMKrQtSiaGkEjXJlkGNUJodoBFNOLMA1EyYNielAsL8L0IEwDLFGZzErzjRCikvKC8bS3kYqFVGzwSPA4SMNEevLYRg5bz2Fr2fHHHJaIU3IGKZb6KZdHXrCnQmgIfQZD5ZPpybbRmaxm74hOx2jx0Aivxuo8LRP6EN599BS3ki6nUByVRfZZzEufhre3FisvqfaqzGj0U10wkeU8udYtpKc+TpJdbBg6ldW9l3MwWyLSsBbVv5kZg4tZOngxbfhoixRJT76LdaFO/rjvIga9OzFiTxMywnzulM9wZduVKOKNeb3Dm9VRDXwhxETgt0ADlTrMz6WUP3zROmcDdwGd44tul1J+42+99+sl8MeyJb6zcg9/2tBDdcDDW5YV2FP8C5tHNuHX/FxdfQ3Te5bStyWL40jaZ1czs9qD2JtAFizwKMj5I2Sa1xM3H6FUHkIIDwlxJc+OLWLbYIzt/SXKtoMiJO0TCzTU9WEZ+xkuHWAw31cZSQOEiDBHLGSyPZNYsQ6jEETJGlgZQSnz0rZqBfCoAr+v0pTgMxQ0TaCMDzNUBIjxRwVAVFqXHEfiOBLbqQzDtB1J2a40fRRLDoWihem8fBewogr0gILmB7w2jr9MyZMjp6dIaWOMiSEG6aXP6SJP9iWvj3lixDy1BNVqPKIKxYlilyMUiyHS2QDxpI/hjIP9oqYWIaAm6KE2qFETgGq/RcxXImxkCekpgnqcoB4noA7jEcPYVhLTSuA4hysgXpkQKkJUPrFIaSHlKw9HfT5NC6FqVVjKZJLmZEaLTQznqxjIBOlLqXQnLEazf92n+oigtSlOMDxAWT3IYHEfQ4UBkDCZ6SyxllMfn4rVa2CVnb+GvGlDtkwxtp/szPWkQk8ylDd4cvBiHu89hZQzQnXzk9jaNmaNLOXUwYuZYvuYHrbJTXiQPU1P8cf957OtbOKtW4lQC7xz+jv52IKPEvFE/tfnzHX0He3AbwQapZSbhBAhYCNwlZRy5/PWORv4rJTysv/Nex/vwLdsh1ue6ea7K/eQL9u8Y0mYfOhOVnevpN5Xz7tCHyS2u42BPWk0Q2Hegjoma2DtjoMjUefZJKc8ypj1EKVSP2CQ1S5l08hZPHLAT3e8iKrA9AklGuv7sDz76MptY6w4CsBEbTILndNpKE4imKnGGdXJx/865NBnKNREDaIBnZCh4lcFhiPRyzaULCjZhzr5XonEQQq70vEolVc39PI5igCPCh4VqVdG11iKwERStP86wiaVs0hkypiH2R/Dr+KJKIigjeUvUPBkSBtjjCqDDCjd9DgdJMz4S14X88So9tYR1mvwiio0GUOaYcqlELl8kGTGy3BaEM+VOdyvuaqIQxdihX0qIQ8EPZKQxyHosfGqFrpSxlBLGEoZQ7PxjHdegz0e8JXrkhWhgdCwbJ2yo1GyNUqWRsk2MB0vedMgUdBI5GEs7xy6B9KLrweoDxs0V5eIRpIY/lEcfYCh0j4Opg/gyErxOkVvZ76zhInJGah9YQqJSkHTVOOlrcpDNG8isiamL0F+9kaS1WvIO73siM9n7dAVrO+tQgvuo37iBormAU4aPpeThs5hoqMzK6xQaniMnikruavnFB4enoJRvwLh7WV+7QK+suRLTK+a/up/P1zH3DFt0hFC3AX8WEr54POWnc0bLPDXH4xz/V072DWQZmlbhJPnbOcvB36D5Vj8n4n/xIRnF9G7M0kgorNwfi11uTLmgRRCV2BxlnjzfYykHgCgpF/MptFzeaQjzP7hPIqAhW0WtQ076So9SXemC4AmbQJLxXJaM7NQB8JkBis1PEVAY62PhqhBzFDxWQ5quoR83lWmUthYoRROTQo7lsTypXGM8SYENYstMlikscnhSBOJieOYSGkiefG4daVSe0Ubr8XqqMKHKvyo0oeCD8XxoTheNDuIaoZQikHUgh+R8yOyXpSUH7UUfGnhoSuIgI7jUbF0lbICeVuSKTukcibxtEk2+6J2elEZkuiJKIiQhekvkPckSRqjjKoD9NNFr9lFspx8yc/Rp/mo8dUS81QT1KrwiRgqIaQVwDT9mGU/5bKPQtFHrqCRLlikCibpovWSTw+vhVdXqAl6qArqRINF/N48hjeN0DIoepqCHGKk2EN3pusFw0Or9GpO1pbSVppLJNmINaCTGx++G/CotE8IUO9R8CYrvw+WN0Fh1k4ydRvIOFuIFyNsSVzFw10nMZAtE6vbir/maex8nlNHLmH6wGImCJWZER2r+hkG2+/kodHJ3NV1Jmb0UfTIFqq9tXxm0ae4bMpl7uibN4BjFvhCiEnAY8AcKWX6ecvPBm4DeoF+KuG/42Xe48PAhwFaWloWdnV1HbH9ezVM2+Hf79vFTWsP0hTx8s4zS6we+RmdqU7OaVjOZan3c+CxJKqmsOzMZmr7M5i9WURIw1naz0jkThKptahqgJT6f7hn/yIe2pMC4KRJkuYJe+m31rIvuQeB4LTAuZySWY7RFyPdV6mJGrrC1JYgzUGdQNGCscKhmrrtKWK3DGPW9VMK9VBSBykxQNkcQlKp6UkJRdtLzopSdBooODXk7WryZoSSE8BydCypYdoqpqNhOSq2VCrNOsKpTDgI4aApDh7VxKMW8Sh5PGoeXclhKBl8agKfMoxPGcWnFXlxFgihYWg16KIaXVZXRtmUYmj5GGo6hpoII0YDKLb+wtd5VETIwPGqlHWVIpA1HdIFi7GMSTJRwnrRaCTNUAjEPOhhIGBS9uUpjn9SGFOHGKafAauXkeLIC25R8XyKUAjoAUJ6iKARJKAF8WkBPKoPgTZ+AzkNRehoQkcV6vi9hxzk+ONzo54sWaBo58enHDkzR6acYaw4dqim/vztNvobaddn02q1U11owpeOYsc10kOlyhW0QCRqMKnBT51XxZ81K78XEuxYkvyMHWSi68jYz1KwPOxIXsC6odPY3OcHY4SW1s3kjKfxZ0OcnXgbDb3ttOgKM8M6VvhZRmfcyRPZIHd0XE7c8yzemjVoquAf53yAa+dci1/3/z1/Tq7j4JgEvhAiCKwB/l1KefuLngsDjpQyK4S4BPihlHLa33rPY13DH82W+OjvN/FMZ5z3LmugEP4TK7vuZ0JgAh+NfYHEQwaZeJEZi2qZFzIobRlGCejIcwfpVX9BNrcbXa+n1/knbt85hQ1dacI+jfNOTjKqrmTzyHoAFgZO5fTipfgP1pPoLYKAKS0hWqMGYcuB4TyYDlI42G3DFCfso+DvIC8OUCx3A+BIQcKcStycTbw0kZFiLSO5EAMZg/6UpGC++OcqQSkhlDJCMTF0G1110DUbXbNRFIHjVMbl2874uHxHYFs6+bKCYxtIxwCpw2GafDQFIr7KMMwqn03MVyLmzRIxkoT0EYJaP2HtIEG15yVj9A2tBkNtwCPrMcx69EINaqYKLR5FDIUQxRcOARWGghLx4Pg1LF2lKAQ5W5IpWiRyFslUiXzqpc04iirwRwx8YR0tCMJvYxklTKNAXs+SV9PktQw5kSYtkpWLsuwMeTOP6ZiU7TJlp4xpm5Sd8qHgVoWKIpRDkyY0AkaAoB4koAUIK1GCMkJQRog61UTManzFMFrei8xqlNMO6dHiC+5zFIgaNNb5qfNrxBQw0mWc8YB3VBOrvY/ChN1kPJvIlfdiOSoH88t5ZvhsnjwYomiXqG/cQ6h2E0Ol3UxOzeGsxFvwD1YzyacyPahihrYTn3U3W0omf9l/NV32GIGGB3DUBBe0XsCnF32a5mDz3/GX5DqejnrgCyF04F5gpZTy+69i/YPAIinl6CutdywDf1tvin/63QbGcmW+cHkN9wx+k85UJx+e+lEmbDqFnh0Jqhr9nLmoDrFpCCdvYSzTGJz4O0biK/F4J7Gn+M/8YXOM/SM5mqIGZ5zUz/7yvexN7KbeU8/btH+kqnMKox2VP9wJE4PMqPcRTpexh/OVe1ROSFOecoBsZBtpayOWna6MtbbmM1g6mZ7MZDoSMfaMQK70XEDYBPxpaquyhEJpDG8CqaawyFBy0uTtNFkrheW8dEjj38Or+gloIbxqEEOE0AggHD/SDmGbQQoFP+msn3jaS7kU4Plj9hUBdSGN+pCkNlCixp+l1jdCldFHVN9PQOxBFc8fQirwGPV41GY8srFSIORq0NI1aGNR5IgHii8sQIReKRBEQMf2qJRVQQlBwZFkyjaZvEU6XSafMSnmzFccoqnqCoZXRfeolRvHqQJFVVAUMT4vkFLi2BLbkji2Mz7vUC5YlHKHv2UEVJqpglEPsbBBzKcR1RV8toOaLWOPFJDPfYrxgt0Wp9zYQca3mVR5E45ToGwHOFC4mC0jC1nXHSJTsgmFB2iZtJ1h52nsks2S9EXMGTwDLa3THtaZbCiUgjsYm30P+5w4dxy4mq1pD6GmFThGF9Nj0/mXxf/CoobD5oXrDeBod9oK4GYgLqX85Mus0wAMSSmlEGIx8BegVf6NjR+rwL99Uy9fvH0bNUEPH7vE5sc7vooiFL7e9h/03CopFSyWLZ9IXX8GsyuD3uIjd8Y6uuM/Q0obGfoUP3lmJhu7Usxo9HHy7L1sSt1JX7aPqcE2rpYfxNkSIzNaJFbjZe6kMDUlC7s3UwmbKTly0zcxpj1EoXQQgJzTRmfhQnaMTmNDr0E8V2ln93rytDQliEaGcIweEmYXo6UhnPERIoqjUm01MIHJxGQtEaeKoBXFawYwyn4US0XY6vgQQwVpURma+ZznZoVEaCB0OT4c08ZRLWzNxNSLFLQceSVNRkmSIsGoGGSYfgp6pnJ18KG3E8Q81YT1OnyiGtWpwipHyOfDJFNhhuJ+LPuvY/MVAQ1hgwlRaAwVqA8kqfMNUGV0ElF3gd37gp+dovjweSfiUZrwOA0YpTq0fDVauhp1LIpMKjiZ8ktDXROoQQMlZCD8Gs5zHc4CTMBEUHYkpjV+cVbZwbIlluPgSHBsOT45iPHwV5Xxm6upCpoCHkPFp/1/9t48Pq6ybv9/n2X2mcxkm+xNmjRp0nTfC10oUGjLVkClspRFgYKIqIioqAXUR1QERZFdBVmKIC1bi0DbdC9Nt3RJmqbZl8k2mX2fc35/nDZpaBF95Pt7ngdzvV7zes05556zJtd9n899fa6PiEES0IsCuhPeR3IsiRBOoHijJL2xYYUDBJuMkJsgltVKOOUYQbkOf+QQiqKFoYJMpDF4HntdJexqgWhCwZHipqiogYBUTXe4lbxwCeeGrsDWko8xAeMzjWQlkoTtRxioXEeT2sXbTZezoycXc/Z7CNYaMkyZ3Dnl61xaculIQZL/4/h/TfhzgS3AQYZUet8HRgGoqvqEIAh3ALcBCSAMfEtV1e2ftu//14SfVFR+9m4tz25tYlZxKvOnH+GpQ49R4ijhbtuD1Lzei9mu58LFhSTeb0GQBcRFfprFRwgG67GkLOK9jht5flc/KUaZL80Ps9X9FK3+VqY6pnNR+DoCe3WEfXFyC21MKbCib/RoiT9ZSSITDjFg3YQvtA8ALxdS3XcuH7U5aOjVRrnpKRHGFHYgmOvpjtXRF9F0+1JSx3imUZaYSFokB2PQBl49Ua86jNx0ApgkgRSLjMUoo5c1MtKddM8UhKHY+yl/C4oAigJxVVMrJZIQVxTCMYVQNEkwnCCuatuHQQC9RTM6U6wxYqagFkeXXXTSQrN6DL80MEyjn27MJE2fg0l0IiUyiUfS8frtdPVZGQgMnxjIshkYlSaSmxIjx+rDae4mXd+KXa5DjTcNEuPg9evSMBpy0UtZ6FUn+kQGcjRdm0cI2BB9VlQfKMG49lz+2YlaERA0awUEAVVR/ilFFAJIKXpEh07zQbIPELN0EjG0EqaJYKSB+AlFkiDIoJ9MY/BsDvcVs7fdSItb+7vISQ9SVNiAR/yI9uBxDAkTC6OXM6ZrGolemUyDyHinCVsgQtB5gIFx62lPdrOu+TKqOiswZGxCl7odvSRz0/ibuL7y+pE4/ecEI4lXZ4Cqqty35hAv7mrlujm5RB2vsK75XRYVXMClvV/h8EYXeWV25lWkEa5qR863ELlgL8fbH8JoyKZduI9Hq2S6vBEun25FTXuLD9veo9BWyC26e+jZCNFQguKxDiZkmjRNfiyJUBnDU/E+PYG1KEqMmDiRg94r2diUS60rjCyqjC8ZwJ5eT3/yEC2BBgDyGc2s5LnkBIrR9acQ6VUGE6VSzTJOhx6HScYqCRiTClI4AZHkMALT3DJjKFIYRYqiyBHQJ7XkqkHHTEWTZiZkhJgOQdE+YsKImDSeqIx1CgQ0wzODREIWiQsCEVUlENdi6v3eGL5wcpheXzaI6B0C2GOEzV4GDD24dC00qnW0JhuHdQYZpkycxgKsYg5S0kkklIHHa6ejz0h/YLi6yGkzUJimJ9eeJMcaJNPkJt3YQZqhBZ3STDTaOcymYvB8ZDsGgxOdLh3dyYQsxY6k2k4YwhkREwbEhJaAJSg6UAUEVQRV0JLWJAVVjIOc1BxJxThJOURS5ycp+UmIARJ4iSa7iUQ6iEa7UE8pViNJFgymsfTFJtHoG8MxdwZHXBKNfSEArAaBccVuLI56uhP7aQs0ISoS89TFTPTMRWmyQFylPMtEkV5ECoQJ5O9hYOw6XPE+3m9bxoetkxAdOzE7q0ioIZaNWcYdU+7AaXb+6/9AI/hfixHCPwMe+/AYD79fz83zi+gxPsXGto18reJOsrdPo71ugEnz8yhXkkQO9WOckkrP+Ofp6n4dR9oi1jTfygs7OynPsbJwRgNrmp8hmozy1dyvkblnIj2NfvLH2JlZYCVZ0wuKijAljLvkHXp96wARj7iCNxtms7khjKLCuIIY+aOO0BTZjCvUiR4D88QLKA9Ox9iZQbBHIweHWWaU00SGUcISVxA9kUEHzYTOTyKzl0RmH4mUPhJ6D3F5gIToJqb2k1C8/GNPgU+DiCzakAQrMjZkxY6cTEWO2pHCNm3E7LUjDdiRI/bhkkyzjGLREdOJhFXwRhXcgTgud2SYNl9nlDCliyiOKGHLAL3GDlqkeuoSB/EmhmSXelFPvm0UGYZ8LEIOQsJJMJCG22OntU+lLzDcSsJqkMlzGMmx68lOSZBliZBu8mM3uLHrerFI7QhKH/G4h3jcTSLh/ZcSqT4NkmRGJzswGLNRxEK8iUJ6Qlm0+1Np8Zhp7EvS2Bsc9M5Pt+oZWxDCkdpGVK6n1luNP+ZHJ+iYJ19IpXsONKYQDymkW2TGZ5tx+KMkk0ECFTtx571HTyTEB+1X8mHzRLDux563gYjaz9y8udw19a4RPf3nFCOE/zG8Wt3GPa/VsGxyLumFb/Nq/at8d8q9mNaPpafVz8JlxaQe7ifuCmJanEKj5UF8vn2kZt/Jr7bNYvvxfq4/O5MuwzPs6NrB2ZlzWea7mabNXvQmmfnzcrAedZMciCJPF+gZ/RK93r8jCEb6xJv5W90Udjb5cZhFZk9oxy1VUTtwAEEVWGhcylTPQmJHjcTCSXSyQGmBlXyLDrMviuqLoaKStHqIF3UQzWwmZDpGWGk+QegQV2R6w1n4EoUEErn44pl4o6l4YxYiCR2hmEwoLhGKCySSgpZZq4KiauEgWRLQSyo6CXSSiklWsOgTWPRRTHIUsxwkRe/DKvdjkbswix049P1YdKEhPxlBRi9lYSALQzwHXTgL2ZOB3JuG1JuGqOgHn4dg06NYdUR1Iv4kuMMJXO4IHs9QpqkoCdgyDcjpSaIpPgZMLtp0xzmmHKY12ELyFHJOM6YxylZIuj4fo5AF8QyioVQ8fjsuj0LHQAhf5HTPfLtJd6J0olYi0WGCFKOKSZfEJCcwynGMchSdlEREQRJOmM+hkFAlTe6qyCRUmVBMIhDTE4jp8EVE3KEEXZ4Ind4w/o8duyDNRLFTJtXRi8HShU89xpGB/bgjWmgnx5DLfJZQ5J5AtElHxJ/AqBeZUGQjJ6kg9EeIm9z4p27Fbf07fWGRjV3X8EFTOYrpKBkFHxJQ26hMr+Rb077FzJyZn9n/0gj+92GE8E/Bxroevvp8NWeVpHP2tAP8/sBj3DjuRsqrz6eppo/FK8qxbOskGYhhuFKgznc38bgXo/Mhvvu2GZc3wrcvSuXt7p/SHmjnnrL7UNbl4XGFGD8ziwoJYrVuJKeB8LnVtHr/gKoqePV38Od949jXFiDDJjNvcgt1kTW0B9oYYxrLheGrMDfm4u+OIutEJpbayRNBdAUhoaKYwkQrGgll1eAT9hJL9KKq4Ik56YjMpS1YSqc/nTaviQ6vOhTJESOIshebNYjVHEGv1wpxSCeLcQhatigktYlaRFAl1JMfRYeSNBCP64nG9ESjegIhE7GYfaoXOAAAIABJREFUGTVhOyHT1GDSCWTZwGmJkWEJ4DT3kmFoJ1WuI1V/HIN0UiUkYNTlYRIKMcTy0fty0PVlI7VnDPPLF606VLuBqFHGr0J/KEFXXxjPCQ940DoCu9OEIROS9hBeSy9d+mYa1TqaA030R/qHPf90Yzp5tjyyTLnYZCcG0iHhIBFLIRiy4g/qBssOuk+UI/y0ilX/CBa9hMOsx2HWkeswkZmiYDQPIOh7SYg9+JIdNHjqafY1D/4mx5zDbOMCSkOTMbky6G8Ik4grGIwS40bbyJVF5K4gaiJJtKQR35gqBtRt9IbT2eS6ng8bC8DYRHbRBjzKMQpsBdw55U4uKLpgxPfmPwAjhH8C+9s8fPmpnZQ4LVx3fi8/+ejHXDz6Yi5sv57DmzuZ94UxZB91E+sKYr3BTk3XjYiiAZ/5Ub7zRh9mg8wdFyV58sj9yKLMqlG/pPG1KIIocOHSQoQtHajxJNK5Udqsv8UfOITBdj5vNt3Aq3vdZNv1zJvSzIHAX+kKdjHBOplLgjfg3ysTjyQpKrRRnmbA2BNCDcZRHBEikw7gte3CHzuAShJvvIhjgSUccRdT222mJ5AEVAwmN9kZA9hS3IiGbsJqJ554N5GkFgNGESjuKWBMVw4ZPjPGaAwlOQDq6TFtDToE0YwqWogaDbhtCVoyXTRkN5KQTynKIVmw6TIwC+mISirJmJ1wyEG/147HawfFONg21y5RmJog3+Yh29JOtqGOdN0+dOJJLx0Bk6EAk1CCKToag7cAuSsHOkwIJ8M+IkjpJlS7gYhewptU6fHH6ekO4nNHhzoCUesIbFkGBEeMsM2L2+DCpWuhLdpCe6AdV9A17M0ANBfQNGMa6cZ00kxpQ9WxRDM6wYyIEQHpRCKWVixSUZIgaglrKgkiSpCkGiGuhgjEffSEe+gN9dIT6iEQH/INEhDIteZSljKWCnUKWYEidH0p9DaECJ+wSLZnGCkrsJIjCYjtfs1l1RYlNHk//SnrCceb6AqVsaHrWjY1OtCZOskr3kRvsganycmtk27l8tLL0YnDk9tG8PnFCOGjFSk5/+EqdLLAvVcI/GD7N5mRPYOb4/ey+60WpiwqoCwQI1LnxvblLA4HbyWR8CM5/8xNL7Qyxmll2fwWflfzS8Y4xnCX4X4Oru0lNdvMojnZRDa2ImWaiCzeS5PrV8hyKu38kEc2G+kPRPnCbB0d8vPU9B1gin06F3lX4NkrkoglqRyXRplRRDnuRRVVEpNb8eZvwR2tQlXj9CfOYr/7fHZ3ZHOsV7PCzc7sIzurA4zN9MSO4o97kBMi5x+YQXbPAKifZocsYDKkYzCmIupMiKIOQRBQlSRKMkYiHiQc7iOWCJzJ9PKMSJiLqRndy8GCwycmgiFF5yBNn4+JPJLRTDyedNq7HURjJu0sBChM01OcFmO0o59863FyDNXolLrB/erkVCzGcszJMRgDRei78xHabSinhnwsMnKWhaRNT0gW8cQUerxR3N1hfL3hYVp4k02HPdNESqYRya6SsIQImjx4dL24xR76Y324I27cETcD0QGCsSCBeGDQwO6fhVVnxaa34TQ7cZqdZBoyySQXZzQfsz8VdUCPpzNCf3tgsBiKxa6noMTOqBQ9tlCcZLNPs1Y2CCQndeDJqcId2UhSidEcWcL7rRewq0XCbO1mVMkWOmLV2A12vjr+qywvX45RNn7KWY7g84YRwgd++s4Rnt7SxAs3T+C+6uvIMGXw44yH2fJ8I2UznExLNRD6yIXt0lzqdXcTCNaTNupPXP+XAVLNem6/2M8Du77PgrwFXNl3O0c2uSiakM7sLBPh6m70FTb6pq+mq/tV7GmL+Uvt9aw50Etlnonpk/bzZssLWHQW7kj5HuFNdsL+GJUT0imVBZRGLxgEYnPr6ba8TCjSSFzNojZ4LR82lXGoK4okJRlb1EVK+lE6Y3sZiPaTPeBkya48VOVjNsmijoycWbhsKkHhs/XB/zgyFBvpgSTu7v1EEmcuNo5goqMwjQ/Kdwx2BBlGJ5mGYgxKPpFALq7eTNp6h0ah2SkGKnOgNG2AIttxsg27UWKHB/3wZdmOzVKJhbGYgsXo+0ZBp5mEKzxUDEUUkDNNyFlmklYdIUnEF9cmi729Ybw9YQKe6Gnz2CabDrPdgMWux2jRYbToMFhk0CsoUhJBp4CsDla+kgRJUy/FRUiKSEkZNSoSDSQIB+KEfFH8/REC7uhpHU9ajoWsPCvZVhlrJIHSHiDepb11iTY9VHrx5eykP/EBkWgHipBGXXAFb9WPpb4nTprDTUHxVhrD27HpbKyoXMG1Fddi1Vs/oyc8gv9r+I8n/KMuP0t/u4UvTssnJf9tVh9dzfPnvMjuR92k5Vg4f3YWvneasJ6TQ0vOQ/T1byK3+HFuXi0RjiW5/yojP9r1DSZkTOA25YfseauViefkUR6KE60fwLTARrPzITzej7A7v85/bZnBnhYPNy9MY3/019QN1HFp9hXMaLiEjoM+nPkW5pY6SB7oAb1IfG49LvNLhCKNJOUJVHV/hTcO6gjFkpTkhCgoqqEhtAlvzMNZjVMpqxsel87InExTxumFRZIGhYPj89hpnfVv38OPY2yilgk9x0g7nkBUzhwXzk04EF0NeP3NZ9wezC5mzeQtxEVtEjPDlEm+uRyLUkzQn09LZyqt/do2QYCyLCvjs1XK0vsoTqklhd0EQ0cH5Y06XToptolYpHLMkRIM/YXgMhDvDpIcOKXjk0V0mSZkpxkpzUjcIBERBQJxlUA4TsAbI+SJEvLFiATjRAJxYv+gTOInQW+SMVp1mG06bOkmbOlGUqw67LKAWVGhP0Ks1Uei90TugCygH2VDGTOAP62a/sgHBEPHEASJpP48trqW8tZhA+5QnNG5PjILNlPn24ZJNnHtuGtZMW7FiGXxCP6zCV9VVZY/tZOj3X6euCmLWzes4KqxVzH76OUc3enii9+YTPTFIxhG2/Et2EBT828YVfwg33prFA09AR6+Jouf7PsaGaYM7s94lO1/aaZsppOpFpnwvl4syzKp404ikTZ0GT/nu2/b6A9G+dZFZl5pvZ9QIsR9BT/DtUYmEVM4e24OGa0+kv0RpJlJOkc9hce3C1U3ns09t/DaAYlwPMn8CX6i1vUcHthDSsTGFRvSGRqKijhKz6NNHhi8zoQEW2aVc1RXfto9cIQ83Fb7Iiv9f8MgnK5O+WfRoabzhvkcXiq5lMb0ok9sNzZRx1mHj2DwnN4RpClWdB3HCAU6TtsmyPnsnevngLkGAJ2oY2zqOLL15QjRYnp6c6lpjQ4qbFLNOqaOsjMhJ05pWhv5xn2EQzUEgw2cvFcGQzY223hspnGY4mMw+kdBn5lkT5h4b5jkQGTYCF/Qi0ipRiS7AcmmR7LrkSwniqtLIookkFBVkmhJV6qqavkQgoqEgISKqICUVCCcIBmIkxiIkByIknCHUfxD9haiRYd+lA2pQEck8xh+/V76+jcQjrQCAtaUabRELue9hgI2HtUkqXPKI8jpG9jbvxmLzsLV5VezYtwKHEbHv/g0R/B5xX804a/Z18Fdq/fz02WVvOv+Ae3+dp6b/DLrH6ljyqJRjI3GCR9xY/+ak+qjl5KRsYg3W1by9JYmfnN1Kb+vv4O4EueR0ifZ8XQHOaV2zqlIJbCpHeuiHBoc38fvP0hK3jPc+GIAk17i5sUBnjzyU9KMaXzP/nNq13iwZxg5f2YWsa0diKl6wud/RKvvCQRBpoMf8MiWVPqDMRaMj5B0rGN/3y7K/WXM3jI0Ms3MnkFj6tDzCmUneX7slcOud27LLl5rvmdwWQGeLL2Un+fcTlQw/Vv38iS+MvAX7q57kdRYaNj69fJMvjXl+7jNqaf9Zmy8jgV7DiNGT0/bL+qN09934LT1gpxP7UUJqpN7SCgJREGkPLWcUvtkTImx9PXlcaA1TGOfFgLRyyKT8u1MG2VjXJaXMSn1KLGD+P2HCIWaOMnsOl06Nts4rNZyrKaxGBOF6PxZqB5IDkRIDERJ+qIkfbEz2zL8KxBBshuQU41IaUZ0TjOCUyVsacAX34/XuwePZw+qGkMU9aSmnkVIuoDNrWW8ecBNpzdChtXAORODuHXvsrtnO1adlasrrh4Z0Y/gjPiPJfxANMHCX20i125kxQU9rNrxYx6c8xOif80h7I/zhRXleJ4/Qsr5o2jK/AkDnl2kF6/l8j/UcuXUfMx5b/DGsTd4bu7z7HnMgzXVwJJz8wm83Yh5hpOOssfo6V1HzujfcvNqE9GEwvevhAd238349PF8JXwvRz7ooaDcwVk5FiL7etBPtNI17g/0uT9Eb13MC7VX8/daDxW5BsaP3857ba+RpWazeN2QTj1j9Hk0GTWNvSIrvHbWfNxCBgC6RIwd275MPpoP3TNjlnJf3nf/jbv+38e4xGH+XPMDCvxDbx5rdPP4+swfEZf1w9qmqm6u3FmFHBs+JWxVjNjajhMIdfFxyAVT2X9uKwf6a4grcWRRZmLGRCZnzMSqVNLVk8GeVi+HO7wkFBVBgPLsFGYUpTK1wER5Ri8WoQ6//zCBwFECwWOnlEIUMBkLMFtGYzIVYjYVYjTmYdDnoFMzkWIW1LCCGkmgxhXtc2rhdFFA0IkIeglRLyKYJZKGIDGxl3CkmXComUCwHr//MOHwSctvEau1nLTUOST1Z7O1NYe3a3rZ3+ZBFGBuaQZTx7rY53uDvT17cBgcXDfuOpaXLx8pFj6CT8R/LOGv3t3Kd18/yMu3TOfe3V+iMKWQu40/ZcvqY1xw0zhsW9tREyryCg81R1ZSUnIv9743ntouP4/flMltG67n2nHXMrnmIhr39fKlr08k9PwRDCUOvPPW0dL2BHmF3+c7746loSfAQ1en89O9d2jH0f2M3WtaqTwrh3GxBLEGD4aFRppSH8QfqEWf8UN+sL4AlzfCVfOS7A7+jo5ABzfuuBB1QFOoZGbPpDFVI5WkOcHTM74weG1373uSu30vAfCFWT9jq/HsT7+foWrGNbVh7vvXilgkRAFfjkR3rpUq8zn/0m9/1PUwt9e/Obj8R8NSvjfzOyAOD/dc0LuB4iOnT/oW9Sv09+w9bb0g5cJXy9kZ3EmtuxYAu8HOnJw5zMiag00Zz9FOqG4eYG/rAKGYFoPPc5iYVpjKtMJUJhfYGJXSRyx6nGCwgWDwGOFQC6Fw8xksGLR6uLKcopUnFE0Iok7zu1G16ldJJUIi4SeR8BGL9Z9WcN1oLMBmG4fNOg5byiQG4mPZWO/n/SPdVLcMoKpQnm3jsilOrOmHWdP4EvUD9WSZs7ih8gauKL1ixO9mBJ+K/1jCv+aZnXQMhFm1XODrG77OYwsfo+0pPRa7gUVzc/CuPU7q9WPYP/BlRNFAj+FpvrH6IA9cNo717h/QFezi6coXee+xOqYvLaSkP0y02YflNgt7jnyBnJwrebpmOWv3d/DL5YU8fvQbiKLIQ3mPs+PPbZRMyWSGXUd4fy+WZWkcUW4jHvcQsf2ae95MopdFbrxwgGfrfk62IZvzXtdIUEDCVzEF0Pxvnlxw+eA1vbj1m5yX3EuDLYe5U18543VnKD1cumcb+pBG7Hk9QXz9tae00GFNOHD6AmT43FjDAUQljqAkSUgCIYOZgZQM3LYMfCaBmNo9bP9GUyZSTjkuw8nJxjiHJ+WwxTr3U5/Jlr3LKfUPjd7HzHqbgNE2rM1CbxVj9w98/KfkecHXeYa/B9FGybe/yO7IHrZ1bBtMtqpIq2Be/jzOzpmHFB/FnhYv1S1u9rQM0O3TQmUGWWR8np2J+dpnXI6d4gwzqjJAJNJJJNJFNOoiFusjFu8/Qeh+ksnwiTq2cU2XL0iIogFZl4Is2dAbMjHonRgM2ZjNRZhMhfQEVKqbB9h+vJ+tDb20ubX7V5GTwoWVWcwvt7B3YB0v175MT7iHMY4x3FB5A0tHL0UnjejoR/DP4T+S8Ht8EWb/14fcsXAM3cbn2Na5jdUz3mTNL/dz7opy0vf1oCYU1KvaOXzkm1ROeIbLn1FIt+q54YIe7t+5ip/O+RmBl50kYkmuuKoUz0t1pCwtot7ybUKhJqTs17nm2YPcfk4xzfJj7Onew+OTn2X3E31kFtg4d2IawU3tWC5wcsx2D6FQI7qs5/jKi30UpJpZfPYR/lT3OHNMMxn7ukaqdvsY2nO1CTj3GHg1bxkAlx57j6c6f0aHOY1pM9447XqvaH8b5/EEaQkTiWPbUFEwJzOZe2Q3sqKNbpvKklSWeyg0/2tSTUWFfX4LbU1WxtYOxeAHrKlUFxcTF7QJRYs1n978XJKCAvoE62dOpVkq/sT9bt57NWV+bfJ2QLVSOX8tijhcbXRN3Rps3af/Nqetl0Dg9GpoenMuCx68g13uXWzp2MKB3gMoqkKaMY25eXM5p+AcZmfPxhuS2N/mYW+Lh5p2D4c6vURO+M/rZZGSTCulTiujMywUZZjJc5jJsRvJtBkw6j7ZPjipaD4+nZ4wLf0hjvcGqO3yc7jTS5dXK11oNcjMLk5nflkG55Y7CaptvFz3Mu80vkMkGWF2zmxWjFvB3Ly5IyUFR/Av4z+S8J/b2sQDbx/h7TtncNPGi7mo+CLObfsyhzZ3cP33puN+bD/2JUU0pj6I338YJfMNVjy3m99fPYWnGm/DorPwg5RfsunFoyy9pRLjB60IOhHlS83UHr2H0rKfc8trWfgjCX70JZW7t9zFt6fdjfWdSnx9Ya68rhz/S3WYpjtpL3uEvr4PyRj1B77yiqpN7C7t45d7HuSi9IvJfOEgAM7smRw/EcLZPTeXPZLmeeKqWgBA7vwNKMJwsrm16g1Mih756A4EYP7RXqwRH41jkiye2o0kQlyVuKz0EfbmTRr2W1FNMjZaT27YhRQViCs6BFkhbNTTZB6FS8457b5OddWw+ug92NBGpwkFtu124mzSiLonxUn1aG3ELstWAmMqiQtJErY4z0z94hmf1dmxbby+4/uDy8tG/4qdo2YMazNKaWfpltP/FgoMWXhrtqAqp4eDJl14E9OvPpdtndvY3L6ZrR1b8cV8yKLM9KzpnFNwDvPz51NgKyCRVDjeG6S2y8eRLh/13X6OdQfo9IZPq5xl0klYDDJmvYR0QqkTTSgEo4nTPHpEAUZnWE68RTiYUZRKRU4KCTXKBy0f8OrRV9nfux+jZOSi4ou4puIaSlM/tRjcCEbwifiPJPxlv99GLKHw9UsD3LP5Hp5d9CyHH42RNTqFueUOfO+1kH53KTsOnkNB/gqe2n8R6w65eOlrBVyz7ip+OPuHyG+VEAkmWHZ5Me6X6ki7oZy9/cswGLI4EnuYH795hMeuruT3x1ZilIz8JPN3VL1Qz3nXlePY2YmaVBGudXG47i6Kir/H19eOpcMT5sdf0vPg7m8yM2sGZc91ApCeMYnmTO21/YX55xEUbIhKgs4t5xESdRTP+2DY9a2seoOcoJ5A6w7y/ClMatzH8fEJLh7fA8DkKatxpWQDcE3tWmw9Q885JaEnpa8f70A9/1CCIsoY0ksIpaYRkobi0YqcpGr6eI4aNLdFW8THto+uw6lqI/3NtWlkHtAyPKvHjKfHonUO1qKz6TJFiesEnj3rsjMesrNqwaDH5o+ct/BUxTXDtktqgts2v0uC0/1tbrnuZl7+/q1n3O+tT7yM0W5if89+qtqrqGqvosnbBECxvZj5+fOZnz+fyc7Jw2wIIvEk7QMhOj0Rurxh+gIxPKEYgWiSSDw5WOjcIIuYT/jmZFj15NhNjEo3U5huxiBrnbSqqtT01fBmw5usa1qHP+6nMKWQL5Z9kWVjlo0obkbwmeA/jvDb3CHm/WIj9y4p53DyUY70H+GP41fz1m9quPDm8Vg3tyEaZeKX1lJ39AdMmvwGC3/bxaJxWeSXbOCFwy/w1qL1rL3/MLMuHU1Rf5hYqx/jyiT7DlzH+MrfcN1LNiwGmS+f18GDOx/kyYVPUfdYAnummUWznPjWNeO4roj9nuXo9ekciT3KqrdqefTL43j82K2YZBPzX9ZCLXpdCv1jygDYNXcU+6SpoCi4tizELxsoPfvvg9d2w6G/YewXsdVWkxlOY0b9bnzn+5mV4ecbzm+zuuJS8pJtXLy1GgGBvKAOb+uOQXsEc8LG7Po6jPHIP30/o5KO3WWV+PT+wXU6Yxru0UPhmpYyG+tyzgNgbvsuXjuuSUMP9FnQf6AR2bbyyXgNfkDEXzEVgN5RBl4fveS0Y558qwFYt+C33Mjwt5PpuhZ+ERvFq5vXDlsvqSI/vP9HfPinN9i/7tnT9vuF+35N4QTtXrf6WtncvplN7ZvY061JP606K7NyZnFW7lnMyZlDvi3/3wqrKKrCkf4jfNDyAe81v0d7oB2jZOS8wvO4YswVzMieMRK2GcFnin9E+KenZ34OcLBDkzDOHZPOHzd9xCUll9BS40bWi+Rlm+jvDuG4rIT67l9gMhWxpzMTf6SNiydm8/ND65mdO5uegxohlkzIIPjkAayzcujqfRZJMtOfmEV990c8eFkl7zY9TbG9mMyeYvb5D7PoxiICfzuGodRBn+Utot0uCop/w1efPMbcMRm0K+/QFexilfBdmtEmXU+SfbQgqpE94NqyEGAY2d9a9QZ57iT+7r0sPtBIR24DFcu7+H7abVw2YTlLe/7Oyqo1pCaMJI7tAcDuNTO3+fjgPkZd4cIy4czujyEMmDlzfH90rI1Dm52k9Wnjb4VGqoIJwpIWSpnWXUJh1hoA3Hk6shdUAXB3zVPcnfEiAOE3VExRAa8ljW1oHXlqdA4rW9egmBM8dYoKKXtBFS8cu4tFnftYUnUnLqDuGzWcc2IytzpeyLkCrLvjWoqiDn7z9O8ASAoKq1atYmH5HL69+m0iwQiP33wDalIzLXvtJ98C4MofPELRxFKuHXct1467lkAswK4uLe6/vXM7H7Z+qJ2HJZupzqlMzJxIZXolpamlWHSWM94jgLgSp9nbzMG+g1S7qtnRtYO+cB+SIDErZxa3TLyFRYWLRqwPRvA/gs8l4Xd6tBCC3ZIglAhRYCvA1xfG7jSjujUi142y4WuoISfni6yv7cdqkHGmu+kKdvG1yV/j+F97cBba0PeECCZUDBMd9LSsJyPjfN4+2I8kCkwvEfjFuj18bfLXOLa9G7NdT4Ys0h+IY56RRYPrDRyOWWxsTMcX6eabF4xi5eY7WFK0hObHNbIvHLOMQ2h1Wv9YfBUAz3x0H8AgaQKsrFpDducAfm8jSw8cJ3yhlwtSg2TP2wiiyMqqNehUCWNdNQlg8YFGRFRc42NULNc0+s9aL2bh1LvhlBFlSbSBkr52Mr1e5KiK32akP9XKR6lTCYvDJYA357/Mgz1PAFo1xL6qJE6X1gHsKrMRdWsdizU+k5Uda1BR+dWCW/gVtzCt6wDvXH6n1vZtlaUH3PhNKWxhBzYgOnYWK6vWEEwXeGG8Fu65rvRRHMUe6rZqy+W/mUjbvZ2849rHyiaNMJccDHCr4yD3r1pFT2Mnjz//FAAb63awcdUOVq1axbde0u71M3fdj7drNwCv//SbANz14lokWcKqt3Je4XmcV3geqqrS4mthZ9dOdrt2U+2q5t2mdwfvQ4Ypg2xzNnajHYNoQFEVAvEAfeE+2gPtJBQtjp9mTGNm9kzm5c9jQf6CkZDNCP7H8bkk/C5vBJNOIpjUiC7HkoN7IIot1UDypMOiLUIyGcJkzKPTEyY/1URXSFOMlKWWsdnVQ+XZucTa/AhGmZijncRxL5mZi9jwfg9zitPZ278VgPNzLuT9w82Mn59HuKYXwSCRLOgmtK+RUQU38dO1nZRn2+hK7CGSjHCxfQm7OQLAIZ1G9sZZQ6nxF4e30GrJGFz+yo43yAsa8J0g+8iFXqamBsleUMX0UDXTd7eT3dlH0NvMzKYAGb5uQku8TLMHqWCo45gSPMDKU0IgtoQOfedxYsGewXVOoASBc5ytdKcPJ/yGsXlkV2j7soe9HBUuBcCXEJn12n4AXKm57OUjbIA6RiNxRVR4at4VZOdU8eCeX3PzxWtJKlD/Kiw94KOqciIc3XXC2XE8K6vW8MQCTZ3kkRyMWbCehqrFAOh+nssl9/WzrEgme6N2zCc9o6nf/jovn3Ulq1atYu3jq9nXo8lQV61axapVqwD46qM/BuDZb/4QT6dWR/jRay5j8e0PULlg6uB1CoJAkb2IInsRy8uXo6oq3aFuavtrOe49TquvlZ5QD56Ih5gSQxIkzLKZstQyzh11LmMcYxifMZ6ilKKRcM0I/lfhMyF8QRAWA78BJOAZVVV//rHtwontS4EQcIOqqqdn03xG6PSEyXUY6Q5per5sSzatAx5yiu0kvVEEnUhU6AXAYMylwxMh12HCFdRkfmlCJoloF9Y0A8kOP7LDQDSqacdNxlG0uts5qySDJm8TNr0NWygdJdFEXlkqsQ0tGEocDPi2ACCZz2VvazV3X1DGBy1/I9uSTfVDf9JOVBwi+UeN5wDwo6bVAMyc/vrgNl1MwNe6nXOPdBCTYcrJkT0wfXc7WW1dBAMdg6P6iuXaRHDOvA2oosSX6t4lrVvr6Gy1ezg5UZseTmFG/X7EM0zc9qV4iBY5SZ5w21SBcnUmY45qYZu3zpox2JG4qhaQsryTjrAO1sLSAfhgwkRiDbuwi3q8YyeysmoNa86exQ+nfYv/it1K446lVCzvpPaVXBYcrsGVmsfeUWCrrcZfMZ2VVWvYPy2HndZZBDCx87aPmP0HTbUk/SQdVnlxLZzMjTvXsC5cxMZoCft765icWc5lt1/FOV1uHnnyt8Bw0gf4yiMPEotGeWyFZkux/vEfEQndy7QlZ84jEASBbEs22ZZsFrLwjG1GMIL/C/i3y98IgiABvweWAOOALwuCMO5jzZYApSc+twB/+HeP+4/Q6T1J4C4AMmQn0WBCI3BPBMlhIBbVthkNOXSE/T3uAAAgAElEQVR5tQ7CFXRhkk2IQQMA1lQjSU8UyWEgcoLwY2QRiiW1DiXYTY4lh8CJMJG2/yhyqoFwpBNZTqE7oO1rbHYKTd4mKtMrURUtMWjaeV8+7dxvb3182PItVdroHsAYjzDpC530q7bBMA5AKNDB9JYQIirlJ8g+e0EVqiixsmoNad0xbLV7sdVWkxU0s/TAcZYeOE5xoJryqzqoWN552qd84REurDky2LbAa8JU9xG22mocMZlLtu/m1qo3Bo/1/azbyTPFBzub8w/WkB2woigxbLVavH7Ztl1YlCAhvYWSWe8ADLbPHuig2K2pY062n7xnKEFr2ZEwyfv6BpcHerVtf5y9bHDd4kNDE9H2nDTuufPbg8st+xuG3Ve9wcC3V7+N5URd101/GjZGGcEIPpf4LOqdzQQaVFVtVDVjkleAj2vuLgOeVzXsBByCIJwu8v4soCis7Ps5S9XNSG9t4LZ3VcLNHmKBt1ASfWzqqOag2Mq7fQF+zXdRQgI/iT/MDOUg2a9uYcUWmeb9B4gF1iLpIrw/sJsGtYtX+nQ8J9yOp7WRP+geoVToZPpf9rJod4Lj1duJBd5CTUb4QD1Al+rmqYF8/iZcS7R+A0/rfkWu4GfZC41MPTw0KXqw8+inXo6IgLd1GwJDhSymnf3a4Pe0Bm0fTo9GgAJwYdnvAa2zALDWVgMKC490Mq2hBs9ZQSqWdzJ/fjdzprxI9oKq0z5P5q8YJP+uKVEmNh9i6QEtRp88vpNsj4KAMNjpPFd+1eA5lX5JI/Gpx4cM0UolTSJ63Zb3AQgahyYtc6/QOt/ytqGiJydxcv8Akjwkl+x74uLB766Fk89478xpQxm8f1zzlzO2WfmHhwe/b/jzmjO2GcEIPi/4LAg/D2g7Zbn9xLp/tQ0AgiDcIghCtSAI1b29vf/yySgInMV+Jiu1pPQEmHtEIRFPosSPoST6cSXddNBPT8LAHmEmrmQaF0s7KVMbyW4PMfVYglgkiBI/jqAEaBV66VEGaE062MFcErEoS6TdFNJBcXOE8qYY0ZAbJX4MMRmjSeqhX/XTkMxlL5PRRQdYJO0lXdfPtKNJitqHqjRFlfg/uJJT7gkgCKdow6UhI7Jk8nR5ZZfVCWidxcnfAxji2mR2frbmctmmZNJqzz/jMZ8vHho55406vQyi6O8/bd1RRduXfIa/Kjn8ydm9Zvl01ZCTfzzBaU16/+H2k7ALn6yo+TiCnk8o4DKCEXxO8FkQ/plmpT4eFP5n2mgrVfUpVVWnq6o6PTMz818+GVEUsGePpsLsY/qExejiKvl5GnmIUpC00VmEDHHGZWlFQWL2dNDbqDD5mFC5kHSvyphpJQDE417sGamELEkqnDMJYSCvRJNQFuk8FI6ZRmHITOnMMdrBdXFkWSZkSjA2YyJuIYsp48cDkK14MecVUBwdcjmsyB/zqdejAilZk1HUIQ38C9u/M/g9XDYNgIhuyPq4Zq+W0Xpy4tNfoUly108qYcDqJPC3TGpfySWa9OKqWnDGT/3uS9nQZ6f2lVyUtWkEDTbenaTdF7ujjM6CzGHHyB5wMVbUJqBrX8kFwJU2avCcavWanLJr4ukl9xpe1dr3nUgUA+hBI/T9c878Ihhb/tfB71/c9rcztlEUBe+Jmr2XzrzgjG3W/eGlwe8X3XHNGduMYASfF3wWk7btQMEpy/lA53+jzWeHlHwYaEIu1QhE9HjRm0z4+3uxOwtobm4m16CNkjsjMSbZ88HXgS5nLEoggNmokae/rw+73Y7X6yXPoI2w28UUKiQDeNvQ5eQQOXwYW7pGfgF3PykpKXi9XnINOvrjCULWXMwA/Q3oC0YRPVoPNq3/q3v/T1Axbdipb0pbzDnu9YPL6+dOYsnWA9i6NW/7d6qzuGi6FuPeOruEuTuPY5aMbBiXy9IDx6l9JZeK5Z24qhaQvaCKJxYs49bNf8NfMZ3RPj072A7YkKQ8zntjO7XKJxdE0duyeHfSUOjFkVZBW5YFPxDMT/BCiaabf2D7r7gl/haKCkdXa+RdVzCFxjRtxHyyw+kqM7A2VVPbnEysWn3MyUQgLur4aLRlWPuBfImdeq1j/lVuG6waSsYqLNfM5dY2bWFLTEsA+1X+kOGaoig88MADg8tTl5512vX99oavEg9r4aSiKRcjSp/skTOCEXwe8FkQ/m6gVBCE0UAHsBy4+mNt3gTuEAThFWAW4FVV9XTD888K9jxo3oouRyOfRFcXtvRM/H19pI+pJBqNkiFoYYSOaFxr721Hztc6CF0giKTT4e/vJSUlg87OTqYatQ6iLRqnwp4P/Y3oChaSdLux6fQIgkjnsTpycnJobm5mzrnaiPLvcSvLMiug9k1s512P6/4HWPLQY6x76VFOfcmZRhd7yGH5hO/hqlrPy0e/zpfHPkaLNBqV/aSkV7B+Eiw9cJz+idIgoZcXtQDjST16iHcnlXBhTRO1r+Qy6osuXFULWJn+XZ6cfwUAX9n/OroUjUzzEg42in7i8U8OY9hsRSTycgkLWhjKKyn89ex5DAjpAPxi689YkXwPgDc351LaqYXU1k8qBnykO8ppztE6jKfmX4IiSJS7jrLp6C2A9iYwEWjOGccRpxbyOUn2HVMl3rJdMngu17587eD3/tsOkw58YdsatsaKhtqUagqaxt11PP/OkJPofd/7wbDr+mjtRra8NBS7L55+KZd/55ZPvA8AqqIQa2kh2tBAvLWNeLcLxetFicYQRBHRakXOSEdXMApDaSmGslJEvf4f7nMEI/j/G/824auqmhAE4Q7gPTRZ5nOqqh4WBGHlie1PAO+iSTIb0GSZN/67x/2HSMmDqBddmkY28S4XtoxM/P29FNu18I4uFMQgCnREYmDPh4496Gdq0wqxxiYcWTn0NDVQeOF49u7dizMcwCQKfNjv44KSc2HfX7Bd+T16H4bE1m0UVI7n6PbNzL7lLg4fPozT3UOOQccbPQMsG38lbPwJthUP4ZIk0g7uGjzVbMmBK+lhRtUu9pwIj3SSwUJXDWgCEp5ccLmmtvFaeW/SFC782z7EK12DpD8jczfTgAKvyHuAJZnJgr/uBOC+K37NE/0PATB2xlq8Zk0KKqgKl3r8ZNYmkOOnR9yS9hg7y7PZYpw/uE5MJnhuy70sVrXkpcMBE7Vva52q05zOu5O0fafaimnNT8MPHJ9m4X3rIqRkHNfWcwHY2Gsn+0MLClqYCaKkOifTmi6TFFSenj9kBy2oSbo2nzu43PjlrTjs1hMa/CLtceOnfuE8lKTCAw8OjeoBfvj9+5D02p/5R29+wJYXHx22/eqf/YGckgI+DiUaJbxvH6GPdhPat5dIzUGU4NBchmg2IzkcCEYjJJMkAwGSAwOgaAMJQa/HOHECltlzsC6Yj7GyEkH8LCKoIxjBfx+fiQ5fVdV30Uj91HVPnPJdBb72WRzrn8KJiUhJCiLo9cQ7OkjJyKSrvo7UVK38nsvlosJiYrsnAAWzYc+fMNqDSA4H/r+/R+msOex646/Mv+l2BEHg+JHDXJBVzFu9Hn5aeQXy7qcxxGsxlJfjfecdxt5yE+8/9Rh2nYTBYOBgTQ2XTZrNM+29tIxdRuHmXyDvfxzHFVcw8MorFFx+HW31mwke+gBOjGrP9uxnm2MyUxe8jqtqAc1bzqNonpbi/8SCZaxkDc7mVt6dVMLS1+HA/CCuqgXMmPQSTyyYwY1712CwT8fa7R+Mt5+/Vo8/qY3Q/xS4mjnWIdIKqgZe1F/AfmcF9bYiAkYr+cEuStwtLPDu4cFdQx0TQEtcR82OTGo7NZL3WDPZOUmbkxAQEcfOxSuG8KKwev4CvEIqqQH3YPimqj8F5/tWsoB3J08ENUhKymg68tLxA3+cf+GwMoyXBN/l6eqHBpej93Zx1o566GweXPfdzFa+Me7iYTp7AIdg4a4fa3Mda3/9HA27hsf5x579JS6+c8XgsqqqxBobCWzeQnDrVkLV1ajRKIgihvKx2C+7FGNlJYayMvRFRUi24R7+AGo8Tqy9nejRo4QP1BDavZu+3/+evt/9DjkrC9sFF5CydAmmyZNHErJG8D+Cz2WmLU4tDUBo3YFp8mQCVVUUrrqPmg/WE+91YbfbOXToEJfNW8T9xztpmryI0bIRoW4NtsUX4l2zlrKVt7Lzb6vpqNlLUVERhw4dYtmEqazt8bDZOo5zU/Kh5lVSLlpK78O/ZlS6E9lgoHrta0ybNoPt27ezbNoMXhRF7u0ReemsOxG2/IrML7yK7733mNlUMyhbyuhw0ZeXzYQDzWxboEkMx8x8m4aPLmbt4Vu4rFKzC3hiwTJuEdeQFyzjXWB6TYjazV1sSFyLTVaomL6GAUsqX2x+i/QWrROpMjiIRzUXS1NbEYbaHZhOqUU7jz3MYw+9KaBLgOOUMrW15A5+j0l6PqqYgy+znfoTc+nOtIn4s7SwhSol+eDsIo4LZaAoPLv5e1zEThIKbFmfS4YPZOvQnIAzdwrHUyL4pARPzh1SBJ1E47bzMCe0+YWq7Ou5auxNsKN+WJum6aU89PAaHmB4Dt89d91NJBTm4asu5uMoO+saLvmGlv+gRKOEPtpNYONGAlVVxDu0TGt9SQmOq76EZc4czDNmIFn/Od8bQafDMHo0htGjSVmszVUkBgYIbt6M7/338axezcALL6ArHIVj2TLsl1+OLjv7U/Y6ghF8dvhcumWiqvD7WWDJYMB8A64f/5iCV1/hT4/8jDEzZiOVjmfHjh1c/fU7mbe/mbuLsvn2R9+B1h2E5v+FlutvJPeXv2Tt5vXIBgPlV1zDm2++yYqbvsKSFg8z7RZe6F8NVT8n8YXXafzKfRhKS3FdeQlbX3meS+6+jzertpKSkoK6+DJ+fLyLZ8ucXLR6EYgS3px76PzxTzFc/kXeaNTISh49iwGj5p45qHzxudi/7yr2pRWzZMIfBy9vWmgPM3a3kd3mIhhoZ259HylhL/KlbkrNEf5kXMq9s7S6thf0b6D40FCcvjSaQX/7TuIxzz99O3VGB8bcSXQaht4OYhaF7dPGUidUAlDS3ci2Oi1SN6BINLyZhTUCIb2VTRVZABhMTvqKNOVO/YRUNqQt4OP44MB1jPe0AtCGkxkL/npam6cCHezds/u09bdfewsfPv8KnUc3nLZt/rV3M+OSc0j09hLYvJnApk0Etm1HDYUQjEYss2djPWcB1nnz0OV9TDGcTIC/C3wdEOyF0P/H3nmHSXKV5/53KnZ17slpZ3NO2lVOrIQCihhJ2AgQWViIcEkS2IRrwICvjW2MrzHxEi1kFCwhQAhJCImVdrWsVlptDrMzOzl1zhXP/aNHu9ogAcYCbOZ9nnqme6q6uqq66j3nfOd73y8LThmcKsiZlFLNBCMCVgoirY1RZmIOPK+Wr18qUXrwIQr33kt161ZQFKIbNpC6/jVEzj9/NuQzi/8S/NHZIwPw2N/Bzz+L99ZNHLziNTS98Y08E9M5uOUJrvnrf+Dr/+8bXHnllXxWS5J2PH4R3Y+4803I199N39s/0yDwqy5l4/e+xXUf/yzf//FPaGtro/LyK/h0/zjfW97By7/3cgg3kUu+m4lPfJL2v/1b7v35jwDBujfexA/uu4/zL7iQz0U72V+p82BnkcV3XAcdq5kqvILM175B/k/fwKYDmwDwF59JVTuW9DXPZeSJiwmArueZqQFcnP05i3YWmJO2yU/vJOW0cvbeRuzevDrDgkhjIvRd7bdy97Jje7sXlB9j/tgE+qSBEpwYXnBDErfbYaCtg43GscS8YKqfz+z5AheKhpfNVKDy5MZ2Fs9Mw29bfBaT4YaGIta2irHmRirmjrWdbEqeedKf6+kt19BVzwLwc3Udrz3vn07Y5u1P/geqfSIpLiPM6N5fnPB/PdTBTV/+IkF/H+XHHqP86GPUd+wAQOvoIHrhBUQ3bCBy1lkooRC4dZjeCxO7YGovTO+D7CHID4P0T3rcLwqhQmpuY8TZuRa610PPGRCK4wwNkb/rbvJ3342fyaDPmUPqda8jed21qPHZAuWz+M/jj5PwM4fg/66HS/6a4e/sob5/P/oX/oG7P/u/uer9f8FDT20nEomgXno1HzowwgNreznlq6dA92mkK5cy/U9foOPL/8od//5NwskUy657PT/+8f288tpr+V81HU9KfhHZh3H3W5CX/R2HP/cQ7tQkyqc/yQ/+9R9ZctZ5OL2L2blzJxtedQ23VDU8KflZ9AAt97wVueACxnctonDfj5m65FU8NdWoehWacxrTMxGEH55zGqN6Yz7ijo3v4WXBDr6z4GI+NOfjJ5zujZvuQXMFHaNpKsXDgMrZfRlSlYZAqn+5x9krMjTpJxJXP51o+Kj4qDJAxadVnJi9U7RVHh9OMf+po73WqWQ3T89tIqARC4q2rWa8uWEFUehUuH3JK1/wJ3p7/lv89bNHRy4vX/IV9nQuO2abNfVnOWfLwAmf7fITlA787KT7ffkbb2Gh7lDZuJHyo4/hTU+DEFhr1jR68RdeiLl4EWJ6H4xshdFtMLYdpvYcJXYtBC1LoHkRNM2HZG8j3TfaClYTmDHQw6BojV6+b4NdhloOypON0UB2ANL7YXI3ZGasHYQC7atg4YWw8OXIztMp/fwxsrd9j9q2bYhwmOQ119D0pjdi9Pae9PxmMYsXwx8n4QN89UKQPoWuDzF2663M+c63+e7X/5nOxUtpOnMDjzzyCNe/9W1cO5RnSTjEPfWfIh78KME136L/1i8jNA3nwx/k/n/9Ry5628388vAopVKJtTe8mRv2DHPL3HZu2XgTjDxF/WVf4vB7P4U5bx7Tr3s1G++8jVOvvo4BF0ZGRjj7T6/n5qkacU3lh+JJuh68FdmylHT1KtL/799wTj2Th72GV4webiE7dx4AtajCt09tkKbqu4w+fjEA/2fFa/mn1nec9LRf13cv8UY4mjlZj9zk9iPKN4MWejNl5o/tR3+RHHxoJI3mokkGOheTjqj48qjyWTcTqN0rmTYbE8JuSLJt3Ty2G+tedJ9vLNzO320/Mp/Ph1rfw3dWvPqE7d669Z4jRdifj+T+nfjBiard5euv5IyURmXjRqrPPAOehxKNEjn3XKIXXED0zFPQqn0wtAWGtzRI3ml45GOloGs9dJ0CHWugYzWk5lEJYMR2mbBdph2XtONR9H1Knk/VD3ClxJcNVaGuCAwhiGsqcU2l1dDoMHW6TIO5loHplGHsaRjcDIc3No4h8MCMw6KLYflV1Lx55G6/k8L994PnEbvkEprf9lastWtPON9ZzOKF8MdL+Ju/CD/9CMGbH+XAtTcSPfdcBs8/kyfvvp1rPvLX/MdDPyMej2NdeQ23Hhjli0s6ue4HrwK7RGn1PzDy7vfR+sEP8Mj4ANOH+7nsw5/ku7ffzrp163hg0Vruncrz3fkWF99zHfgO5RWfZfiWvyJy3nnsP201O372AGdf/0aeGZsmm82y/Iqr+XRdp+z73BEbZv0DN4NmUIi+gYmv3UcA3L/4qLK0uux0fNH4fYYWxri/p1FRSvE9xh5vvB6KtHLGaXedcOrH45LMI8zbn0d1jw2JmFIjFKgEMkAKgS8CXBHgiRPtDqq6zuCSJFub11IVv34Bj4e2v5nVhaO99CsXf4FtXSf637xQyCZ1cC+ed6K9w+IFZ3OqV6CyaRN+pjGSMZctI3r+eURPXYmVLCLGtsLQ5kYvWwZHe9hzzoA5Z2J3nUq/2cn+qs2Bap2Bqk1/zWaw5pD3ThwNCSCqKoRVBU0INCEIAF9K6kFA0fPxjnukFKA7ZLAkHGJVzGJV1OJUM6BrfDPs/wkceKAxN6CHYclluD2Xkds0TO77dxIUi4TPOIPmt7+dyHnnzmb3zOJX4o+X8MtT8M/rYN55pIsXNsI0//qv3Hn3dzCtMGtveDv33nsvV1x5FX8l4ozaDps7ckS+exW87FaG7xqn8uSTJL/2ZW7/h08zZ+UaEme8jCeeeILzXn4Rn4+0s7tc40c9PqvueBXEO8nFbmLis58jesUVPN2Vom/rk6y46DIm9AiDg4OsOH8DX4l3srdc52OxIu/Y+hGUyZ24XZczvtmksvmXjK05ne2iEc8WQqW47GivebrH4O6FVxx5v+HwJr4/+JcATIXinHbGXTjC/M9fs98Sp9Wf5rZtHyXhHU33uanjL/nB0stO2PbS9CMs2H1i6ChZF3gDW0/qx7FhzjpSO7fiDDasrNXmZiLnnE1kzUIinT56aScMboJ8Yz16GHpOh96zKXSdzo7ECrbbKnvKNfZU6vRV6/gzj4AAekIGCyyTeZZBT6ixdJg6zZpCXAYYnovvebiui+8fbRCEEOi6jqZp+JpOSdWY9gJG6w79NZv+qs2+Sp0Dz/u+blPnjESE8xNhLqnspvXgfbDnB1DNQKSVYMk15IZbyd55P97kJOaK5bTc9A5il1w8O8E7ixfEHy/hAzzxz/DQxwmu+y4DtzSsh/2/+ig//MLfsuENb2PndJ6pqSle9uYbuXbPMG/vaeVTOz4Be+7FueI2Bm76KHpnB+V3vp2HvvVVVrzsIirtc9i1axcvf9U1fLiuU/B8Hm6eouvO10DbctLuq5j+4tcJnXIKgxefz7aH7mfO6nVoy9eyY+cuOucvYMcpZ3N3tsIiU+G20o+Y+8vPI1WLonoZ0w/1446M8sy6cxkPGtJ/CZSXH/0NJQH3n30aw8bROG+8WuQnW/+chYwf+cy3Fr6Cj3Z/mED819oGbKhu5NaD3+G0/LGpkl+yXsVn1r0bT9dP+Mx5hU2s2D55xNTteSdD00DfkfTR56Mt0suZo4fwBxuZO0o0Snj9OsLLu4l0eJjyIGJ4C9QaDSThZug9m3L3mexpPoVNoflsr7jsKNUYs4+a1XWbOiujFiuiFotCOp2+Q7JeoVYokM/nKRaLlEolSqUSlUqFer3Ob/qsWJZFJBIhmUySTCZJpVLEmpvJx5IcCBR+WayyJV9m0mmE1haFTS5JWVxfepolffcg9v8EAhfZdSaFyimkf7Idd2gYY9FCWm66ifgVVyBm7SBmcRz+uAnfd+HL54FTpXzK5xm++T20vPe9bMyPM7p/D1d/9DN8+3u3c8opp/DEsvXcNp7h3xfGOP/OV0LgUln7OYY+8HHC69YxcsVFPHnPHZx57Wvoq/kMDQ2x4c9ey83pOgqCe6L9LL7/nWBEKHa9j7G/+ypaczOlG9/Moz+8CyuRYN4lV/PM/oP4vk/Lhov4ppZgsO5wo5nhQ/1fI973E6SRJC8vIf1QH97EJAeXrOegddQdMpSYw3RX+zGnObbM4L72K44/e1pL0/zl7q/wOvuh3+46Pg+blJX8Xe9bebL31GPKJT4fS919bHh6N0r9REJKuBoMPEvg105YF9FaeNmhg4h8g8DVRAJrzXLC85KEW2uE5AHE1E6YEZPRtBCn50wOt65jS2I1j4p2tpdqDcuMGSy0TNbGwywzNbrdOs3lPPX0NJlMhkwmQy6XO4bMFUUhFosRi8WIx+NEIhEsy8KyLEzTRNd1FEWgCgWkRNIQbknA9Txs26ZarVKtVimXy+TzefL5PLXa0fM1DIOOjg66e3rwO7rps+Jsrjo8kSvjSEmTrvLqqM9bMg8zb8/tiEwf0myiKC4ks3Ec+9BhjPnzabn5HQ3i1/5nSmpm8Zvjj5vwAQY2wrcbYZqRext52M3f/ia3fe6TLFx/BuaqU9m0aRMXXXU1n5JRDtVs7u/xWXbHn0C0g0L7+xn7+KeIXXkluxZ2s/uxn3H+m/6crQPDZDIZTrniav7GMxmtu3ytrcorfvrnUJmituqjjPzDXfiFAvL117NlapjMyBDzzjyPWksXff39JNvaGVl/Dj9wBCU/4GZtjPcPfJN4/4NIdMriHPJ9BuVte/Gk5KE1S5EcJTPNSJBbuPiEU/Y1yY613WyJnvHbX78XwcX5R1iwK4fin7ynGfIVrOE+vOd64MdhiZNi0d6jv7Exfy7Wwg6sDpVwdBrD3o2ozRQ+0Sz8zlOYaDuFnak1PBJexuOORf/zrJfnhgzWRi0WKgGdtRKJ7DSlqUkmJycplY46jmqaRktLC83NzaRSKaIhEy0IUDyHoFqmks9Ryeeol4rUyiXq5RJOrYZTr+HZL2z1rOo6umESisawYnGseJxwIkm8pY1QsgnfMLFRyJdKjI+PMz4+fiQ01NzcTOf8BUx1zWW7ZvGzfIWKH9Cuq7xXHuC6wTtJHPopUkJJnkt6q4PdP4wxbx4t73rnbI9/FsAs4Tdw99thz724197HoRveRXjdOsavupTHv/9dzrj2egZdOHToEJf82fX8r7xH3Zc83DpNxx1/1kjVtK9i+gv/QuKG1/Ok5jHwzFOsufxPGBMGAwMDrD77HL7XNp8n8hU+1Kry/ic/iBjZgjfvaia2xSk99HP0xQsZv+witm16DNUwmXP+RQxX6mQyWSLtHUysO4v7HIWSH3A1Y7wv8xDLD92LUsviik6KxaWUBjxq+w4TBAGbVq2jqB4fAxdYXWuZSvz6PT5PhSAcoPiy8cYF7Te8LRJ1EON9+PUXFnQtKIVY2r/7SEBH7+kk1NNEqF3DihUJKYdQnYkj5+G1LCHTupb9yRU8EVnKA+ocDtT9I5Zz3abOCstgvnRprxSJpycojY+RTqcJZjxtVFWltbWV1pYWYmELkwClXqWeSVOYmqAwNUFxeprAPzZjSdV1IskUVix+hLwNy0I3Q+ihEJphNs5DSoSqEgQBvuvgOQ6ubVMvl6iVitRKxUbjkTu2wQtFYzT39NLU04vZ2oFjhJjM5hgaGsJxHIQQdM6bR3HhMp6xkmws1XGk5CIlz4cyP2b1ge8jagVK9dWkt+vYgxMYixbS+u53E7v00tkY/x8xZgkfoDQB/3I6tCwmF3krE5/6DPHrrmVnTyu7H32YC992M1sHRkin02x43Q3cOFKkWdd4ILSL2L03ItZTJmQAACAASURBVJdczuShFeRuu53wBRdwcN0Knn3kp8w75TSMlevY9vQzLFi8mO3rzuV7k3nWRzS+kb+Xji1fAN2i1PwmJm57HG9yEnH5ZeyJmwzs3Yluhek+90LGbZ90Oo2eTJFffSpbwikO1F3iuNxib+OKyZ/TPboR4VTwvDBleznVTJTaYAFndIoA2DN3OUNJ50Uvg6pamMkegnAc29CoaxJHOTYjR5FgBCohDwzbhWoBuzh20hDM8RCozMlJlg0fQptRoarJOGZ3C2argRGzCYXSmOIwqjJTL1fRsJsWMdW0ggPxpTwZXsRPtV4O+kfz/btMnaWGwhzPoaWUIzI1Tm18lHK5fGSbWCxGW1sbcSuEKX0oF6lNT5IfHyU/OU7wvEnWUCRKsqOTRHsnsXiCiG4SDiQh18Os1lGKRYJsDj+fxy+XCEplgnKZoFpFui7SdY8YpR2BpiF0HSUaQY3GUGIx1Hgcra0N0dJEPRqhHjIpCchXS2RGR8gMD2JXGxlIhmXRvmAJ0bkLcKwoo9NpnisCFOvoIrNsNZtDcfbWXFJBjU+VH+WVh76HURiilJvH9O4YzmgGc+lSWt/3XqIXXDCb1fNHiFnCfw57fwR3vAEWX8rU5FlkvvJVmm6+mSeqGQ4/+wyXvudWHt66DcdxWHf9Ddw4MM2ySIg7nUeJP/QXyLYV5PTXMfn5L2EumE/+Da/lsXu/T7KjiwVXXsejjz9BJBrFuuBSvlIXpB2Pv4iXedezn0Ed2oTffgbp4eXkfvwLpOviX3QhB1vi9O/ZgWaGaFt/JlUrxsjEJIGUqIuWMTB3MVvRSXsBZuDwBmcf1xW2sGRqK+HsAQQSz1aouQuwa03YRQ17uo4zlkE6jdBPKRShr3Mu6ZiKK6q/4iK9MFRpkrANujMZutPDR6vnqCp6SwI9FUaLaRixAMOqYqhT6Ookmnn0HnMSvWQTCxmKzmNXeAGbQvN4VOmkPFPRSwXmmTrzlYAOu0Iyn8WcGKE2OXFMr72luZlENIolJKJWwclMURgdpjg9eeS7FFUj1dlFqr2TRCRKFIWo7WHlC4iJCdyRUdzxcWTtxIZMiURQW5pRk8mj5B2LIswQiq6CpiA0FaEIpO9D0Mj6lK5HUK3il8sEpTJ+Po83PY2XTsPzGhwUBb27G33eXNy5veSTMbKew8T4CNODAyAlesiibdlKtI4esrbH2ERj9CPmL2Zg/lI2Sp2a73FzeQvvGrqN1NQuitNdTO9O4E4VsNato/X97yNyxksb1pvFHxZmCf/52Pp1+PEHkeveyPiWGIX/uIfm//0xHty5lczIMJe+7y+47+FHME2T7quv5UMjeZK6yl2Rfhbe/05QDSpLPsLI33yt0Xu65f089JN78FyHFVdcw6F8mYmJCXqWLGXnqtP5fqZMj6HyFedx1m/7PKI8iZc6lezkcrIPbEZWa3jnn8tgRzP9Q4dw63ViHd1El69huu6QyxeQgDdnPtm5C9lnxdnrNIgv7pW5yu7jktp+lpf20VboxyoOIWSADMCzFby6hSub8Pw4nmMQuCq+KwgcCOwAKUH6EgKJDAKEqswsIBSB0ASKKVB1iaJ5qJqLqttoooRGDjXkopkBYob9pVCxox2UIp1MRns5bPWwx+ziKbWdp/ROqmrDDVMAPYZKNwFtbp1UuUAkPQljIwTO0Rh5LBYjFYthaQqqXcfLpymPDlN6HrGrmkaqq4dUqomEYRFzPMKFEuboON7QEN7U1DG3gBKPo/d0Y3R1obck0BImekxBszx0w0bVqihuvpEeWcuBXQKn0hBqeSeWlDwBQmkoca3UkUWGmvBEM54bxSkKnIyNM1nEHhzG7uuDGaM4JZFAWbmCwrw5TOsKo2PDFKYaRJ+avxhz7kJybkA6m8XTDfKrT+WpVAeH3YCrS0/z8bHbmTO+lfxIK+k9Cbx8lch559H2gfcTWrHiP/nQzOK/E2YJ/3j87K9h498jz/8ww98/TOXxJ2j63N/yw5/ei12tcs6b38HDm7fgui5rrn4VnyhL0o7HV9sdLn3wJsgP4Zz2MYb/5SGcw4cxXn0dOyxB//ZtNPfOo+P8i9m2czcA7edfyG1mir1Vm8W65O/rv+CMnV9FFEbwEqvITa8k99hu/EwW2ZQif/YZDGmS0cOHkEC0cw7WvIXUNJOpbI4gCHA0HWfOfErtXYxFEhyUCvmZ5G4jcFhaH+Usd4yFXoZuN0OHk6GpPk20nsXwKuhOGdUpIeSJ4qrjIREEehjXTGAbMWp6jIoeI282kTZSjGophrQkfWqSPWobk2YTvjg6f9CiKnQqklbfIVGrYJXymOkp1KmJIyEfgGg0QiwcxlIEqmPjFTJUx0epZo6qe1VNI9neSSqeJK5oROs24UweY3AYb2TkCGkCqE1NGHN7MTqbMZosjIRED9sYeh7VGW/445QnOWk1znBTI70z3AxWisCI4ukRHD2CrYXxVINAqASKSoCKBFTpofsOWuCiBw6aU8aw86j1QqPRqExDcQyOr2OcmEOQWoLtdlDPh6hPOtT2D2IfbFgxiFAIb+0q0nO6GLWrTAwfBimJ9MzDnLeYqUqNar1OrrOXg0tW8QwG5xZ38Znx77Fk5Alyh1vJ7I3iV2ziV19N63vfi9Fz0nLSs/gfglnCPx5Swg/eBdtvI7jkcwx+/iHsvj6iH/lLHtz8CIXJcc567VvYOT7F2NgYp114EV+PdfBkocIH201u2fYRxKFH8OdfyXRfL7m77kNraaH2+tew+dmtlLNplm64hHwsRd+hfkKWhXbmeTwYbmZHpU6PJvics4mX7f4aavYQUg1T0c4j329Q3rIL6Ti4vT3kli9mylAZmxjDc2wUwyQybxFaUyuOqpOrVLBtBwnUdQOnuQ2npY1yLEnODFNUNHIIsoHEPf5nlhIrsNGl2/DPkQGCAEVKHEXHVnQcYeAJ9aSpl5aApCJICknc94h4DiG7RrhaQS/mUbLTWNXKMaRuhUJELIuQqqD6LlTLuPkM1fExfPtoWEUPhUg2tZCwIsRRiJSrWJNp9MODyMLR9FSh6xi9czC6WzFawxgJiRmuYGjTqLVhKIxwDKGrBiR6CBJzqEa7KUQ6yITamDCaGTZaGNKSjChxMp4k6/pkXY+y71MP/vPPSFRVSOkaKV2lWVVYTIllXpp5ziTd5WFaiocIZ/Yj0gePjh70CF5yDdVqF9UJheqBCeyDhwDwe7rJrVrOiOIzNjaMlBBbtAy/pYOpQolCJM7QylN5KpxgZWE3/2fsNlYObyJzsIXsvhCgkHrd62h+x01oM7UhZvE/C7OEfzL4Ltx+PRx6BO/cTzDytc3Unn2W2I1v45dumYHtT7Hy5a+g3NzJrl27WLJiBU+vPpN/m8hxTszky/kf0Pbk50HRqM19KxN376a+Zy/mOWczsH4Vzz7xKKquM/ecCymGIvQfHkQ3DKzTz+YXqU6eLNUxgJtFP6+fepA5h36MsIv4RifF2lrKh30qu/qRtRq+rlFds4pMezNZ3yWdncZzGkSvxRJYnT2osQSBEcKWUKnb1J+XOigBW9Op6SaBbiBDIQLDRBomUtdBCCSiQY1CIJAonofiueC6CNdFODZqvYrpuRiei3rcfaNpGpZpYmoquiLQAh/sOl65iJPLYGfTiOBoDFvVdOKpFLFQmKhQidguVq5AaGwCdXTsGGmW1tqC0d2O0RbFSAqMcA3TyKD7w4jKcZUyI634qXmU4nOZDnczEu6m3+xkj9HGHpKMOB7TzokeQpqAZl1rkLOmklQV4gqEkZgywAgCjMBHD3w02bhaQjYWpEQqCr5Q8BUFTyi4ikJNKFQUlbIUlCRkXJ8x22HquO+3FMGysM65MsNZlf0sy+2iPf0s+sTOIyMCN7KScnk+5cMulWcPIut17GSCzLrVjKgB6elJsMKEl64hFwgKCAaWrOHp1m4WF/bx2bHvsmpgM9P72igcVFGiEVre8Q5SN9yAYv7+lNmz+K/HLOG/EOwy3PUWOPggwZobmHgqQeE/7iG8YQOD55zK1vt/QOfSFXRuuJRfPP4ELS0tiA2X8n9zdUq+z4eSDu/c9Tm0vgeRrSvJ+Vcw/e17CRwHXnEJ/a0JDjy7DYC5Z51PPdXGwf4BhBCElq3icO9CNnoK065Pm/C41d3O5WM/pXl4I8KrIaVOVaymkm+lfLCA3T8MUhIA9a4OyvPmUIxFKONTrFYoF4/2fqVQkLqOHk+iR+Oo4QjCDCFUtbFOKMgZHxghBMfeBhIhBIoEZNAI/QQBBD7C8wjsGn6tilut4JZKCM+GIDiWpHWdSDROJGQRESqW6xGq1DBzBYzJKYxc/pjt1VQSvb0ZozWCkdIwInUMPYvBKKozyTFHF2nDTs4jH5/HeKSHfquH3UYnT6vt7PONEzxwLEWhO6TTqam0KJJU4BF3bSy7jlkpoZVLiEoJu1Y7IpgKjs/A+S9AKBQiEokQTiTxU03YsSRlK8KkZjIkFfps95jGYJEB1/hDXFB8hqWTTxIZ3YoIXAJpUpFrKY0lKD1zmKBSodLazPiKxRyulbAdB33OfLzmDjKuT9+8pezoWciSwi4+O/RNFg/sYHJ3G5Uh0Lu7afvgB4hdfvlsRs//EMwS/osh8OHnn23E9LtPI6f+KZP/+EWMOXOo3fRWfnb37YRiMVb+yZ/xy937KBaLLFq3ni0LVnFnukiXofFVbTenPvEpRHEEd86VZPrayP/kMWS9jtjwMoYW9rJ319N4tk33KadBRy+j2RylUgnNMGH1Ova39bCx3ggfxHG5wevjyuI2Vkxuxpra1ThUT6POQur1dmpZjfpwAWdkgufY2heCeiKG3dWBHY/hhkxsTcEWAjvwqbsOnu81FrfhCfOroGoamq6jqRqaqqFrGqaqoQsFI5Donodmu5h1G71YQk9n0EsVtOMaACUWRW9JojdF0OMqesRHt2oYagZdjqIelz0URDuoJOaRjfUyGumhz+xml9HOVrWDA75+xI8GwBCCnpBOp6rQik/KtYnVK4TLRYxCDjeXpVQs4p3kfEOhENFolHA4jKnr6JqKrihoikCREiF98Hyk5zYavCBABgEy8JF+gJQBQgiEojZET4rSWISK0HXQNAKh4gOelNiuR6VaJV8oUKkcawgXDoeJdHRit3WSjqfo10LscSRTbuO4W7B5o3eQy4tPs3TkUYzcIQIfKu4qimMpSjtG8GyH6YW9jHa3M5HPIOMptHlLyfgBfd0L2DlvGavyz/DpoW/Qtb+PyZ1t2JkAa+1a2j78YcLrX9ztdBZ/+HjJCF8I8TngasABDgFvkVKeoLwRQhwGSvDcfX/ygzkevxPCfw6774V73wlmjMqyjzH62S8jHQf9ve/m0ac3kRsbZf5pZ2IsWc3Tz+7AMAy6LriYbylR9lTqXBZX+ZvJO+h89ltgF/HazyM3tZTc/ZvwCwWU9esYW72M/Yf7qOSzaGaI9lNOx0+2MDQ5hW3bKOEI7qJlTLZ2sUc16bNnPFZkkevq+zijdohFhX20pHehVhvq08AH12vCkR24dhynauDkA7yijV+q4hdKSPfkxC4BX2mEc577z3OvhKRBeC9yyYRhoMYjqFELNayjRnS0qIJm+mh6FU0to4ssmpJG1Y+9zwKriXqsi0Kki6lwN0OhTvYbHezUWnlKaSGDccz2LbpKt6bSJgKaPId4vUK4VEDPZQgy01Sfl48PjVFLLBYjGo0QNk1MVUFDIlwXWa/hV0o4pTx2sSGOqpfLyF8xia2oKoqqzfw9ukgpCTwP3/MIfB/fc5EvMkJQVJVoUwuR5hbMZBNKNI40Q9golGp10pkMjtPQUwhFwZozl0r3XMbiTWz3FYYdD6Rkgz/GW0qbOXv8URLTOwk8QbG0hMLhGNX9Y5Qtk9HVSxn0bWxVR1uwjLyis69jLjsXrOCs3GY+OfB14rummN7djFcJiF1+Ge233HJi1a9Z/LfBS0n4lwKPSCk9IcTfAkgpP3yS7Q4Dp0kp07/J/n+nhA8NC91/fx0Ux3BP/ygjX3+C+q5dhM4+i/Hzz2Trow+BhJWX/wnDts/g4CDtXV0Uz7qAbxZs8p7PhWHJJwoPsmTnNxGlcYKmVeSrZ5L56Xa88XFEJELt3LMYa04wcLgPu1ohFEvQtGodTjhGplSmOGMBIJqaqS1YxkiyhcNCY8Dx8QGkZJGf5XJ3gDX2KHPrY7RXRkmURzCLI4jnZYJICYEn8B0NXyYJhIVEJwg0pNQIAhUZNHqmAgkigCMR/QBFDVCEg8BGoYoiayhKHdUIUE4i5pVCxYu0UQu3UrLayIZamDKaGdGbOKS3sUdr5lmliZJiHfO5qKrQpSm0CUnKc4jbVcKVEno+i5qewi0dqyhWFIVEPE4kbGFp2gyZO1Cr4Bbz1LNpytn0SW0QzHCESDJFOB4nZIYwNR1TKJiImVFLgOo4qLaDWq2h2Daa6zUaC8chcB2k44DnN66bqoAycw0VBTQVJRzBi4TxrBCeaeAZOp6q4pg6NSQV16FSrVDOZyllMscofcPJFPHuXvSWNjzDouz5TKUzR0YoQWsHxfmLGUq08KwHlUCyyJ3mf1We5NKxB0hO78IpG+SzyyjstakVK4z2djLYlqLkedCzgFIkyZ7OueyZv5TLpx7kI33fInjWJrM/DopO841vo/nGG1HC4f/UozSL3x9+JyEdIcQ1wKullK8/ybrD/HcgfGjUK737bXDoEWTPWeSci5n+1l0E1Srm9X/GLgMObttCoqOT+ZdczTP7DlAul2nrnUt2zWnc6yiM2C7LTcFn7M2ctfsbKOn9SLOJinEuxcEQpc07CEolSKUon3MGo2GD4eHDOLUqEkj2zifUOx/HCDOdLxwx3fIVlaC7l1p7J9lYihHNZCKAtHe0N6lIn14nwxJZoCco0eGXaPOLtHhFUm6BiFfDCGbSB30bNXDQfAdF+kihEszE94OZ176i4WhhamqYqmZRUS3KSoiiapFV44ypMUaVGEMiwogao6RGkOJYWX9UEbRoCikkicAj4dlEalVClRJaPovIphG1Y0M6qqoSi0aJWCFMVUWTPsKu41fKuLkMlfQkdrl0zGeEUIikmojGE0SsMGHNwEIQclyMah2jUETLZiGbw89mkS/iiSN0HSUWQ4lFUcIRhKGj6HpDbKVrCENHKI3ePb7fGB34DV0DgSSo1ggqZfxypaHQLZfhJM+aEomgtrXhdndQa22hHLYoSp98uUB2YhzXbmTu6FaY5MIl6C3t1BSdiUwWx3EIFAV30TKGO+eyTQ2R8yWra4O8r/AoF478BKs0Tnk6SW50DuUDGaaaEwwu7GXKcwg6e6kmW9nZMZe+ufN57di9vHv/v1N8Wqc4aKG1tdJ2663Er7pqNr7/3wi/K8L/IfB9KeW/nWTdAJCjEUX4ipTyqy+ynz8H/hygt7f31MEZ3/PfKYIAtv8bPPwJqOXwlt7A1DMhCj/4EVpbG84Nr2XLnmfIjY/SuWwF8VXr6RudIJ/Pk0il8E8/lwe0KLsqdVo1hb8IdnH5+IOkDj2IcCsEVjsVcRbFQ4LSlh3Ieh1pmthrVpLv7mBaeoyPDR/JxIl2dhPq6EZEkzhCoVCrUS4fjf96ioKXasZPtWDHk5TDUaq6QU3RqAiFMoJSAAU/4FdH7V8YGhBWBWEhCAmwkESlJCJ9LM/FdG2Meg2tXkOrltFKBUQ+hx6cWEgkHLaIhCxCeqN3rngusl7FLxexcxmq6Wnkcf42ZjhCNJkiYoUbVggSQrZLqFzByObRJqfwp6fBdU/4PrWpCa21Fa0pgRYLoUY0VEtBMyVqyEfTPRStjkoFRZZRvCK4lUadW7cGbvU3q2ur6A0vft0CI4w0ogR684wAzsKzdbyqwC15eEUHZyKDc3gQWT8q7BKJBO6CeZR7OsmFDNKVIunx0cYcAoL4gsUYnXMoS4V0LkeAoNDZw/SCZTxrRsl6PheVdvDe3M84dfhB3KxHbnwB+b02GU3l8OK5jCngt3VTae5ke/dCRjvbedfId3nNsz9h6pk4dkbFWruG9o99DGv16l///Gfxe8NvRfhCiIeBjpOs+qiU8gcz23wUOA24Vp5kh0KILinlmBCiDXgIeI+U8sSq08fh99LDfz5qucaE7tavg5Wi1nsjE3dso75nD8b69UyetZ7d+3dTTE+R6Oik88zzGS3XGRsbI2RZhE87i83Jdh4v2bhSslgPeJ+znYvHHybe/zDCtwlCXVSClVSnI1T60th9DVl9YIWor1lJvq2ZoiLIVUrk09NH4syqFSHc2YMeTyCsMIFm4EiouS6VShXPP5GcJOApjZ57oAgCVUWqGmgaUtVQhED4PgQ+0vMaEwRBgBJIDN9DfZEYtxCCkGliGgaGpqIpAlVK1MBHeA5+tYJXLmHns9j53HNJoEeg6jrRRJJIOEpYN7AkWK5PqFrFyBbQp6cRk9PHiKugIUzS29rQWpIz1g46WgT0kIumV9BEAU2mEfU0vIixG2YCrAQylCAwk7hmHFcP46gmthrCUUO4aghPC+EpOj6CQEqeu8qqUFAFGL6N4dfRvTq6V8PwaoSCOoZdQqulEZV0ozCPf/zoQiDjPXjaHGynGacSws4F2GMF6gf6jzQEMpmgunwphZYUk77D5MQoge+jWhHii5fjRJNMF0u4vk+hs5fhBct4Wg8Tcoq8K/sQN4zdRzw9SH6khdyhJJmKy+H53YyEDJzWDortvWzrXUI9pfPx/i9x+lPbmdqZwq8LktddR+sH3o/W1PTC13EWv3e8pD18IcSbgHcAF0kpf6VRixDiE0BZSvn3v2rb3zvhP4eJnXD/rTC0GdlxCnnvYtJ3/gxvchJt/nxKF29gb3qCif6DmJEoc8+9kKJhcfBQPwCJrm5KS1exI5LiyXIdT8JK3eW99ac5Z+JRmse2IKqNEn2+0UXFX041Habal8YeGD5i0hWYJvWF86h2tFGKWJQDj4pdp1wq4DyvZygBFAUtHEWLRNCsCGrIQhgmQjdQVA0pBPK5YfpMHn4jB39mD7Jht9AwiAkarpC+R+C6BE6dwLbx6jW8agWvWgHfO+kEr6ppWJEo4XAESzcICZWQH2DYDkalilEooqezqNPpRk778yBCIbTWZvTmBHrSQotr6OEAzayja0U0kUF1JhDOiVWz0CxkrB0v0kEt3EIp1EzeaCJjJJnWU0zoScbUJKNKlDFhkfEk+ZlatS9V3poCJHWVpKrSKeosCAosDIosdKaYUx2mrTxMrDiIlhtAPM9OWhpJbG0ZtUoLtWlBfTCDPTAEQYBvGpRWryDT1sx4tUQxl0EqCtbcRYi2LjKVGjUEY3MW0N+7mIOobMht44Pp+1k/9hjlEYPsQDfZtMNAbyeDiQh2SxfZjl6emrecVmOKT+3/Ik1b0mQPRlHCYVrf/wFS118/a8X8B4qXctL2MuAfgQ1SPq/C9bHbRABFSlmaef0Q8Ckp5QO/av9/MIQPDQLceSc8+HEoTyBbV1GULyP7yAHqe/agJhI4V15GnyY5tONpFEWhY+VajJ75ZOoOY+MNgVC8q4fi0lU8G06wpWwTAAbwSmWSqyu7OCX7NK3jW1AqjcsZCAtHWYDttFMvmdhTDvbwNN70sdMhrmngtbdhNzdhx6O4IQNXUXCFwCXA8QNc38NxHTzPJZASGfgEQUAQBEg/IJgJvSiqijqTjaKqKoqioqgKqqKi6zq6oqErAh2BFkg0P0DzfHTXQ6/W0Msl1FwRNZdDdU/SEKgqWlMKLRVHS1ho0UavXLN8dLOGphTQmUZxpxDyuCCUUJGxdpxIJ5VIOwWrjXSolTGjhWG9mX41yUElxWFpknF9TjYmeY54U5pGUlOIC4gKSRSJ4TfEVbrnovteQ2zlew0hmeeD782MfgJUQBWN0U0gBD4CXygEigKahmIYeJreEL2pGlWhUhEKVaXxN+MFjNgORe/YowwpguWqzXnBJGfU+lla3E97dg/G1B7EzMggEDGqcgWVbIpKXx67v1ERrN6cIrd8CWOGwmRm6gj509pJulwlFwozsWQ1z6baCdWnee/Uj7h++B784Trp/m4y4z6HetoYbE7gtHYz0TmPXy5cwbn2U3x4+zeoPSmoTpmYSxbR8YlPEl6//jd7jmbxkuOlJPw+wAQyM/96Ukr5DiFEF/B1KeUVQogFwD0z6zXge1LKz/w6+/+DIvzn4NZgxx3w5Jdgei8y3EYtdTXZp8uUHn28YZF70csZ6Wji8OgghalJhKLQvmIN5pwFZG2HsfGGGZbV3II9fzFTTW0cUE12Vx18QEVyuZjiFdV9LK8O0F08RDx7ALU8ceQwfBnDEV14fhLXsfDqBm45wCs4uNkSXq54UhfIF8Nzd8JvOj2nRMIoYQslHEKNhtAiJmpERTVBM31U3UZVqmhqCU3kUf004iRF0mUogRtppx5ppxBuJ2O2MmG2MGy0MKA2s19Nso8YmZOE0jUBrYZOq6bQJCAuA6K+S9hzGgIru4ZeraBVSvilIrVq9ZgKVL8Kz9Wr1TRtphFsTExLKY9Uy3rudRAEOI5zJLXyhaCqKtFoFD2RxIknqUXjVKwIlVCEKVWn3/UZrB+dj0iIgEuCMc63+1lf2sfcqW0Y6X0AeEGCiruMylSEyv4pvHSWuqGTWb2ciXiYyVyaQNUwehdix5soOS5DHb0cXLiCCXzeNPkA7xm7C2toisyhDiZHBX2dLQy1NOG093Codym75i/kDVP3cMPm+0g/E8OrKiReeTVtt96K1tr6a1/LWby0mBVevRSQEvp/3iD+gw+CauJ0XUm2P0XhgV8QVCqIWAzv3DOZbG/h8PgI+clxhFBoX7YSvWceVRQmM9kjxGPEE3jzF5Nu6aBPtzjo+JT8o8Q4nwoXeqOcZg+xtHqYluoE0dokofI4SiV9Qly8IZDVCUQMX0YICOMHJkFggNAbSSNSnPC3kWoYIPBRhIcQPgIXaOu3OwAAIABJREFUITwELorqoMpyY3JTcRAvUmsjsFJ4oSbqVhNVs4mSmSRrNDFlNlI1B9UkfUqSvcTJcGIdXE1Aq67Roio0CUlC+kRdm4hTx6xV0ctF1FIBWchTrVROWndWURTC4TCWZTXmGXQNXVXRBIggQPge0nUak8SuC56L/1zYyq7j1uv4jk3gB0dGRXJmee41QqDqOpquo+o6qtbI6lF1HaHpaFYEPRxGDYVBNwhUDQ9w/IC661K3HUrlMvXnT9oKQaS5Gbe9m2KymUw4xqhq0OcGFGbUxMtkkdfae3hZYTsLJn+JmTuElGDXWyiWFlM6UMMZmaKua2RWLmU8EWEynyGIxAgtWEZBKkxYUfoWrWRfPMVl0xv50Pgd9AweJH2ghfERk76uFobbWqh1zmXX/JWkO+N89OCXWP34XjL7YyihEK0fuIXUa2fDPH8ImCX8lxrTB2DLl2H798CrESQXU+FUSgNQenI7QaEI4TD+2Wcy1dXKYHqS7NgIAKph0rRoGUZ7F45mkC4Uj+ThIwRaSxtuexflZDMZK8KYojPgBsc0BAB64LFc5lke5FjoZ+lxc8SDGhG/StSrYnkVLLeC6VUw3TKq7zYaiOesE2RwzHup6ASqga+aeKqOp5h4qoGrGLiqSU2LUNaiFLUIeTVCTg2TUcKkhcWEEmFIiXNYRHDFyStvRVWFFk2hSREkZEAscAm7DpZdw6xV0Mol1EKeoJDDPUlPWVGUhk2BZREyDIwZYZUS+AjXIbDrBNUyTqlAvVikVsy/aGlCIRT0UKOalW6ajcpWuoFuGA2lsaIiaPgMiUAiBA11raIc6e37vo/nz4ivggBfSgIkvu9j16qNcoml4gsqnCPJFJG2dsymVpRIDE8zqAeSUrVGvlA4UgpRUVW03vkUOnsYiabYJxVGnca6hV6G19V384r8VhaMPIpiF7ErUYrlZZT6fezBScqmzsTyxQwZgprroHTMwW/tIi0Ffb2L2dO1gOXFHfzl2L+zpu8ppvc2Mz4eYv+cNkbb2yl3L+Sphatos6b5q2f+BW1TncpkiNCyxXR8+m+wVq18kYdlFi81Zgn/d4VqFnbeBft+CIefAOkjoz1U9LMoDeqUtuzEz2QbWTHLl1GeP4dcxGKqlGd6ZHhGpq+QmLcAq7MHGYrgCIVy3T7mgZeAkmxCpppwo3Hq4Si1kEVZMymqGlkEOR8qQUD1t3B6/HUQURSiqiCiCKJILCSWDAgFjVRNy3Mw6zX0ehWlUkYpF5GFAtI9ebgjFAo1SNw0MGZ64WrgN0RPdg2vUsItFqjls9RLxZPmtquahhVPYkWjWFaYkG5gqhomAt3z0VwX3XHRbQe1bqNXqohKGVlq5Mz75fJJUzt/G4hQCCUSQUQjyEgYLxLFb0rix2O4lkVFgYrnUq5VKOZzlDLpY5S/RjhCcu4CjNZ2PDNC2XGZzmaPiLH8RIpa7wKmm9vZr1scdnz0wOXV/7+9O4+Purr3P/46s0/2PZkkkz0kQAi7gKC4Q21rf15rS2tXa/W63aqtVWu1aq+97bXbVWtbRVu3WqkLKCCyySJ7yEYWsu/7nkwy+5zfHxMUKbggMIGc5+MxznyTL9/5fEd45+R8z/ccRznfHNpDfts2jKOduGwGhkdyGKr24egYoCc2ko6sNNrtI3jMIWhTsxgSemoSUqjInEbSaC0PtrzA7NoD9FRE09IbRKU1nu4EC32pORRk5LLMvpMbdv6LgcIgvE4tkSu+RuxdP0YbGnpKP0Pl01GBHwijfVD9jn+Vrbqt4HUiTTHYzYuwdQRhbxzAXlnz4c0/SYmMTc1hMCqMXpeDvt5u7EfdXarR6wlNtGKMiUeEhOHR6nD7JE63mzGHk7Gx4w+QkoBPbwCzGWk0gdE/U6bPYPRPC3DUfv7XH3YMaaUPjceN8HjA7Rq/09SJcLvxOZ1oXU40nzCmRa/XYzaZ/AGu06ETAi0SjdfjXyrQ5cA3NorHNoxjcADHyNBxQ1xnMBIUHu4f7WM0YdLqMQmBwePD6HKjH7NjGBlBOzAE/f34BodOHNoajX8Fq+BgNCEhaILNaM1GNEFGNGYDWpMWjV6D0PgQGi9CIxFa0OgEQjv+EBLp81+8xevxX0zW6AAtCB1SCrxOHz6XD5/Ti8/l9S864/Lhc3rwjtj9K2F1d/uXSzyGiIzEk5qCKzEeR3gYQ0LSbxumt7Ptg99UhFZHqDUNQ1wCXnMoNreb/gH/0FN3VAxDGbnUh8dwyAsen4/FYzX8wLaPxZ3bCB2owz4QxFBvFsMVI4x5PHRmWGmODGXY64WkNMaCw6mNSqAsK48EZxM/a3meedUFdFdE0TAcSmVKAgPxSTRk5VObauGuxr8z/70SBmqC0UaGkfDAL9SkbAGgAj/QnCNQswkq3/Y/u/xdNtIcj0M7DbstmrE2F/bDjR+OvtFo8CUn4Uy2YI+KxGbUM+L1MGgbYri399/mfRFaLabIaAxhEehCwtCYgxB6A1KjGR+C6Y9m/+JW0t/dcMxcOfKD/xw5KGg1GjRC+B9IhOSoqYH9s2gKrxfpdvln0XTa8YyN4Rmz4RoZwec5cUvZFBKCOSgEk8mMUa/HpNVhkgKDx+ufjG10DP3wCNq+AejvhxP8UNMEBaGNjkYXGYY2LAhdiBGtWYvWBFqDF63Bg1bnQqe1oxU2NAwjHINgHwTPZ7uwfcyHDhqdf1oFnwe8H3+R9t8Yw5HmCLwiEo83BI8nDI8nFI/LgHvEh6ujH2djC97+oxZAN5txpaf6h+aGBjGIj77+XsbGZ0o1RcUQkj4FT3AYvSM27A4HLq0OW1o2nYkplOmC6Pf6mDNWy52DW7mw7V30I/0Md8cy1BbPaF0fXREhtGSm0O3z4ItOwJuQTG1oNIcypxPtbeee5hdYUHWAjkNR1LgjqEpNZMCSQsnU+Rgjxnhk3+PodjlwDBgIXjCPhEf+G0Nq6sl/zspnogJ/IvF6oLsCWvdDy/hjoAEAKfS4g6djdybjGjHg6nfj7BzE1dL+kRE3MjgId4oVT3gYrpBgXCYDTp0WB+D0ebC7nIyNjeJyOHC7nMdtMZ8SQvj7uA1GjCYzRoMRg06PXqPBIDToJei8XnRuDwa3B53dgd42im5gEE3/ABwzW+QHNBq0kZHoIsPRhYegDTWiC9ajNYPO5EWrc6LT29GJYbQMonH1ffxNVfrg8aUGI/zPpvCjXkeAMQwMQWAI9u9rOOahM4FW/9GAF5p/XxhGSv86C14neFz+Z6/L/9o96u/ysw+MP/d/9HmsD4bb/CtjfaT2IDxGKy5PHM6xEFwjBpy9Tpyt3XjGR3tJrRZPTjZDKcn0GbV09vcwOuxfGtMYm4DZmo7TYKZvxIbX52MoJoHu7OkUmsKwedx8eegAt/RvZnr7djxDPga70xg8DP0eL80ZVlr1GlxhkcikdJpCoziUlUeo7OJnTX9n9qEi2stiqNaHU52RSk/KFAqm5HO5632+v+V1+kuDQeqIufVWom+4wT+DqHJaqcCf6Gzd0HoAWvZBywHoKPEHxDipMeIxpuHyxuO0h+Ia1uAecOEZtuMdsuHpHzjhvDBSo0GaTf5FT0wmfEYDPpMRn16P12BA6MZHVUj8F2vHb7gS4O+ukBKtzz8xmcbhRONwgt2OZmwM4XR98hBOIT5YAFwbbEIbGoQ2xIg2SOvPUYMXrd6FTudvfWvFEFpvP8I5cILjaSAoBoJj/EsQBsd8dPvI48jXzZGgMxz/WBOR0waDzTDYBANNHz4PNPofR/5eaHR4grOxO63Y+43YW0Zx1DTjG/PPx+SIj2U4M43+0GC6RocYGxlGCg2hWbkQY6HXNorD66UvwUpbRg4l+mDMriG+27+d7/RsJKH7EMOtYfQ3JjDUY6clOZ6maP+wUV9iGs0RsRRlzSDR28zD9X/BWtRES1UM5RGRNGZk0JKVT2NqAvfUPkPm1gZGWs0YM1Kw/Oa3aoqG00wF/tlGShjphL5a6K/zP/fV+R/99f+2LqpEg08XhVdE4ZFheDxBeN0mvG49Pq9AesHnkUi3D5/bh3T7kG4vPrcX6fH6W6pC+PtahWD8ptujXguEQYfGoEPoBRqdQKMDoZNotD6E1odGeNAYvGi1LrTCjhZ/14nGO4TwfszC3xr9UWEd5Q/po0M7OAaCY/2PIwGu+ZhxoOcyn8//22DnIegshY5S/7PNv0CM9IHTl4LdnYZ9wIy9oR9XS5v/B4A1if6cTNqFl57uTqQQmKwZaOKT6Lc7sUvoSE6nJSWLCq2J7NF67ul5m8ta1+Nq99HfmsJQnYP2qDDqUy0Mmcz4EtOoj7ZwMGM6sx2lPFDzDMFFA9Q1xXEoKZ6ulHRKp84nJqKPe7f/Fec+DR6HlqjrVhB7191qJs7TRAX+ucTnhaEWf2tvtNf/GOv1dwWM9vq7Bo68/rhujlNFo/d3h+jHH8bQD7tMTOH+LpOjXx/ZPhLqxtDjrpmrfAYjXUf9ECj2jxAbXy/BrbViG83C1gqj5U1IuwOH2cRAXg5doUF0DvTilaBNSEJnsTLo8mDT6GhNz6Hcko7DbeOGng3c0P4mQZ1d9DclMFCro81kpC7VwmBQMN6kdGrirJSk53Dp8PvcffjveAo9VPTHU5ZupSs9l+Kp+XxtcB3LNuxgqDYIfVwkCY/+hpALLgjwh3fuUYE/WXk9/ouSHqd/gewPnh0f/ZrPP8oEMT6vuxDHbI+/1pvH+7mDPnytVX2yE46U0HMYGnZCw3Zo2gX2AXxesDvTsA1asNX7b8hyazX0pSbTm5JI++gQbo8HfVIaMi6JgTE7bVHxNE2ZQaXBzBV9u/lJ9xpy2gsYbAynrzaSVnTUpVnoNwfjSU7ncEI6lWmZfLV3Hf9Ztoq+g0ZKZQI1menUTp1LX0oEv9j/BMHbh3GN6An/wuXEPfiwWlD9FFKBryiTmc8HXWXQuBMadkDTbnAO47absbnzGGk1MlrehNfnozc3k7b4aH/L32jGlJmLTWukT6unPmMq5XFWEm313NH1Fle2bmC0VkNvTRytUk9tmoX+0DCcljQqkzKoT0nmh22vcO2Bd2k4FEVpZAKNU3I4mL+I+RRz/TuvMVRuRhsSRPzPHyLsqi+rIZyngAp8RVE+5PX4R4lVrvUPFR5qxuPUMTw6neEGPfaadux6HV3Tp9Bs1jFiH4PoeERiGoMeH/VxydRnTmUQFzd1rObGltdwVbvpqY6jzaunLi2RnohIxhLTKUrNYSzGyAM1TzF9z2EOt8RTlJFKU+5MarMzuKfyaSxbO3D0GQheOBfLr3+LPuF4s7Ern5YKfEVRjk9Kf79/5dtQ8Rb01eCyaRkeyma4XoejrZ/+sGA6p6TT5nPj1OrQJqXhCA6n3RxCdVYeTcFBXN++mltb/4Wv2kFPTRxtbj1V6Yn0RcUyZM1ib9pUUg1tPFT6JzR7HRS7LFTkTKE8bwFJ4Z3816YXGCnyT98df//PCf/qtaq1f5JU4CuK8un0VEHlW/7w7yzFMahjeDCb4Vqw9w/TZYmlPTWJ7rERZHQCMjGVDq2Byszp1EVG8o2OddzZ8iqamhF6quNo8JmoSklgOCaBzrRc9qdP5Uuj27j5wCu0FodyMMpKQ04e5Xl5/LBtFdPeqWKsx0jw/Hwsj/2fau2fBBX4iqJ8dgONUP4mFL6I7KtjbDCSgY40Rsp7GTLqaJmWTbPXiScoDE1qFt06ExWpOdTEW7i6awM/aX0FY/UgnZVx1GiDqbEmYIu2UJ+dT601kdvaX+LS99+ntDmJoowMDs9YiDbRwz0bn8FxUAN6Awn3/Yzwr61Qrf3PQAW+oignT0po3gOFL0D5atwjTgY70xms0mCzOWhJT6Y53MyYzoCwZtJnDKY8OYvq5DS+0LOFnza/jLmsn47DMVQFh1FvtWCzpHAoaybeGMEDFU8S/X4P+93JVEzLo3juAlZ0rWPh2oP+1v68Gf7WvsUS6E/irKACX1GUU8Mx5J8RtuhFZGsRIx2hDLQmMdI4QltsBM0pFgYF+Cyp2EKjKLOkUpWazWW927i/7m9oS2y01kVTFR1Bo9XKcHIm+7NmMFVXy727n6G9KIgDlnTKZp+PLS2EBzc+ibsA1dr/DFTgK4py6nUegsIXofRVnF02BpoTGKzV0m3Q0ZSWTKdegyfGgismkfJ4K5Xp2fxH51rurHkZR4mkuTWKKkssrSmptKdPoygjmxu7XmXZ1u0UdiVRNCOPwrmLuaZvA+evK2CsW7X2Pw0V+IqinD5uBxxeCwXP4a3dzXBrFP2NMfQNOmjKSKbZrMcdFY8j3kqxJZ0Gq5XrW1/lhqo3GSwxUd8XSVl6In2JqZRNnYsrVsN/F/4B8b6dPRGZlM5ZxGh6CPdt+Suug4DWQMLPVGv/RFTgK4pyZrQdhN1PIstXM9Jqprc+iYEeB/XpSTSHmHDFWBiNS6bQmkVPXAR3NL/A1RUb6SqJ4LAnmop0K/0pWeyfOocLPAe4aesrlNXHUTA1j5K5i7h6cBPnrS9irNtIyKLZWB57HF1MTKDPekJRga8oypk10Ah7nkIWvshos5fe+mT6ut3UpyTQFB6CO8bCYJyVg2m5uMMF9zU8w5KyAzSXxVJqjqYxLZWWKTOpTU/mp/XPkrGtjvc1mZTOWcBoZhg/2bIS10ENIshE4q9+TegVywN9xhPGaQt8IcRDwA+BI5N4/0xKuf44+y0H/g/QAiullL/+NMdXga8oZ7mxfih4Drn3r4w1DNJXZ6GnU1KXHEdzZDiuuER641MoSJ9GmHGQh+r+QsbBempq4yhN8M+4WZK3gOjwfu7f/WcaS0PYlzGd4vOWcE3fu8xcdwjngIHwKy8m/pHH0IYEB/qMA+50B75NSvnbj9lHC1QDlwOtwAHgG1LKik86vgp8RTlHuB1waBXsfgJ7VQO91bF0deioT4qlOToSZ1wSHZY0DmRMZ4qvlkcqnsKw10bFYDxlGVa603Mpzp3BNSMbWLbpPXaOZXJw/vk4M0K4a/1zOMs06KJCSfrDkwSdtyDQZxtQgQ78RcBDUspl49v3AUgp/+eTjq8CX1HOMT4f1GyE3U/gKNlLX1U0He1G6hJiaI6LxhWXTHNyJoXpOVw1tJnbCl+muyCEMkMsNVOyaMyeRa81jAeLnsS+x8OutGkcWHwJ1zWvJvedGtyjOqKvu4bYex5EGM6ihW9OodMd+N8DhoEC4MdSyoFj9vkqsFxKecP49reBBVLK205wzBuBGwFSUlLmNjU1nXR9iqJMYC0H4L3/xlm8k76aeNqb9dQlxdIUE4UrMY1K6xQaUxK5q+lZLtm3h7rDsZRYLLRm53IgbyELKOJ761axZziNXedfhDZJy+1rnsdRq8eQEkvSEysx5UwJ9FmecZ8r8IUQm4HjTWhxP7AX6MW/9PUvAYuU8vpj/vy1wLJjAv88KeXtn1S4auEryiTQsBO2/hLHoYP0VMbT3qWnMs1CV3Q0o9YsDqRPwxw6xqMVjxOxt4/K3jiKc7JoyJlFY1Yyvyh+At4fZVvWbArPX8Jt5S8Qs7UHr0dH3O03E3XjrYhJtEraGRmlI4RIA9ZKKfOO+brq0lEU5eNJCTWbYOsjjJZW0VUWR7PdxOH0RAai4+hNzWFP5nQucu7jrsK/019gpNCURG1ODoX5i5mlK+P69a+yaySTHUsuJTW8m6+vfpuxNiOmvEySH38afWJioM/yjDidXToWKWXH+Os78XfVrDhmHx3+i7aXAm34L9p+U0pZ/knHV4GvKJOMzweVa5BbH2WkqIXOiljqdMFUpyRii0uiNmM6VWkp3Nr2Elfu20ZldQIHU1JozJ1NbW4G95c+hW/XGO/lzqN2fj537fsbuj0OpFaP5aGHCL/6q4E+w9PudAb+i8As/F06jcBNUsoOIUQi/uGXV47vdyXwR/zDMp+TUj76aY6vAl9RJimvBw6tQm75FQMH+2ivjqE6NJz6lERGLekUZs/AFwW/PPw4sbs72T9i5fCULIpnXUCuqZYfvPMq2x1T2L34YhZqDrH0rb3Yew2EXr4Ey//88ZwevqluvFIU5ezkcUHh83g3/5b+g6M0N0RRER9NW1IyA6k57MmawVxZxj37VtJeEERBbAp102ZRPT2Heyufxr17lPemLaJvppX/2vw8zlIN2pgwrE8+g3nmzECf3WmhAl9RlLObawwOrMS98Q/0Fnip74qmPDme3iQrzRl5lGZkclPHP1m+ZRuF3UmUTsni0MwlZIY08b0Nq9jqmcq+Cy7mqz0bSd/QgMehI+6WG4i65c5z7oKuCnxFUc4NjmHY9Uec7zxFV0kI1bZIKtKSGUywUpY7G6J9PFrwe1y7PewPS6Umbxa103O58/Cz2Pc52DRrKfoMA99/8zUcLXrM+Rkk/+nv6GJjA31mp4wKfEVRzi399fDuzxnduYm24jgO68OpSk9hIDWHfTn5fHlsC9/csoZDTRaKc3Moy1+ENbSNb739GpuMeRy4aCm3F75M8O4RhElP0m/+l5DLvhDoszolVOArinJuqt2CXHcP/XvaaayO4VBiPO0pqdTmzKE7OZxHS39P8I5B9pvSKZs1l8q8adx38Cl6SzWsv2A500wtLHtjG64hPZFfXU78g7856+/QVYGvKMq5y+uG/c/gXvdrOvdrqRqOoTzTSq81i33T5rHYV8At2/9BVWUse6dO4+C8pcwTpSxfs4m1KYtomJ/Hj9evxHNYoE+JJuUvL2LISA/0WZ00FfiKopz7bD2w5WFs76yipSSG0qBoGjIzaJoyi/o0Cw/WPEXiey3s0mVSPHch7blJPPjeE5S2JfDusi9xXfNakja34ZNaEn9+H+Ff/3agz+ikqMBXFGXyaC/C99bd9G+ppKYxluKURLpTszmQt4BcYx1371zJ4fIo9uXMoGD+hVw9uJGMjZW8OecKohM9XPPmOpzdeoIvnk/SY0+hDQkJ9Bl9JirwFUWZXKSE0lW43niQjvedlLksVGSl0ZadT1l2Fne1/I1pWyvZRjaF85fgSTNxx/qVbBYzKLtwAf+183lksRdNdCgpTz6DedasQJ/Rp6YCX1GUyck5Ajt/x8hrT9NQHElhdAKtGdkU5y0kKryfB3f+iapDEezLzadkwfncXPcy9gIX71y4jOX2vWS/U4vHoSP21huJvuWOs2INXRX4iqJMbn11+NbeR+9buylvT6AkI5W2rDwKp83gxo5Xmf1OETtMWRScdxGJsf38v7fX8HbyUpgaynVvr8beYsA0fzopTz6LNjw80GfzsVTgK4qiAFRvxPWPO2nZ4aCIRKqnTKEibwEi1sMjm/9AVVUk+2bOoXZuPncVrKSsNY6DFy3hh8Wr0Ox3ISKCSP3zc5hnTtwuHhX4iqIoRzhHkJsfYeTNF6gqt1CQlExrdh5F02fyn+3/IHNDNTticjmw8GIWaYrJ3VjCmjnLOF8cZsa6cv+0DHfdQtQPbpuQXTwq8BVFUY7VvA/vqlvp2NJDyYiFstxsKmcsJCx6mLs3/JXCtmT2zV9I3zQLt295lg26udhmJ/Gd1a/jbNcRtDgf6+PPoQmeWDNvqsBXFEU5HrcDdv4W26t/oro4hgOJyTROnUl53nTurXoa87Y+tqfmU7TgAr7WtRbHwTG2XL6cmwpWYThoRxMTTNozL2DMnR7oM/mACnxFUZSP01mGd9UtdLzTROFoIuW5OZTMWsJUUw0/WP8qu2wZ7Dn/IswJHq5e/QarZn2R+aKO2WtL8Xq1WO69g4hv3RToswBU4CuKonwyrwf2PoXtpf/lcFEk+5PTqMubS2NOCg8XPs7oHh9b8hdQNn8eP9n5F3a6c+mal8n1r6/C3a0l5LJ5JP1uJRqjMaCnoQJfURTl0+qrw/uv22hfU0GBPZnyadM4OGcxl7t28cX1m9lkzGPn0iu4cnAbwbu7WPOFr3Db7pcxlDjRJoWR9uw/MaQFbi4eFfiKoiifhZRQ+Dy25x6iojCMfakZVOUvYCA9gl/u/CNVVZFsWXgJ2lQ9K1a/xktz/x8XOA4xdUMVPqEl6dFHCPtyYNbPVYGvKIpyMobb8b5+J22r9rLPlUpF3gwOzlnEd/reJGdjORsS51G6eCE/3r2Srd4ZOPJi+Poba3H16wi7ZhmJDz2G0OvPaMmncxHzV4Gc8c0IYFBK+W93JAghGoERwAt4TlTMsVTgK4oScFJCxWpsT/+UQwVB7EvPpmzWEgwJLu7Z8hd29mayfekyLh/bhX7fANuuuJRb3nsF32GJLjOetGf/iT4h4YyVe0Za+EKI3wFDUspHjvO9RmCelLL3sxxTBb6iKBPGWD/eNT+l9cXN7PamUZE/h5KZs/lx/XP4tg+xcfZSyA7h6rfXsGreVXylZzux27uRBj3W3/+RkIsuOyNlnvbAF/7bzZqBS6SUNcf5fiMq8BVFORfUbsb21I8o2a9nT9ZUiuYsJSOkiW+ueZ2NwfkUL13CbfueZ6N+NrGpXi5evQvXsJ7w61dg+ckDp33R9DMR+BcCvz/hmwjRAAwAEvirlPLpjznWjcCNACkpKXObmpo+d32KoiinlGMI7+t30PryNnaQSdns86iZPoUHd/0f1bVhbF+wmEu0FbiKx6i8cA7f2/g6zgYdujk5ZL385mmdkuFzBb4QYjNwvA6o+6WUa8b3+TNQK6X83QmOkSilbBdCxAGbgNullDs+qXDVwlcUZUIr/Re2p++msCyK/dNns3fmEvqCg/nhP56nPjoE1zwrl2/Zytp5y/hm07ukODoI+ca30S9/FLSn52LuaW3hCyF0QBswV0rZ+in2fwiwSSl/+0n7qsBXFGXCG2hCvnEj9oaDrB27lD3R+WyduoBHN75BZ2sNRUvm892ad9hons9FcTUs8RXiTJiL8bpXIDT+lJfzcYF/KjqTLgMOnyjshRDBQojQI6+BK4CyU/C+iqIogReZivj+eoKuuJdrwzfzsHYVT5l6WPrYb7gJFfLhAAAG80lEQVT2jntYvnUPr4YtJDemk/pDZt7wXIamowTXEwuQTbvPaKmnIvBXAK8c/QUhRKIQYv34ZjzwvhCiBNgPrJNSbjgF76soijIxaLRw4U8QN2wkLEjL4v03E7rv94RfdCHL//YSK9p7qa/WU79kClFNTTxn/xIjLon82xfxvP+Ef+jnGaBuvFIURTmVnDZ49z4ofAEss+CalcjwNJofe4xtOzZTev4sbhp7l032WVwaUUGuqMeR9UVM1z4Nxs+/YLq601ZRFOVMq3wb3rodPE5Y/j8w57sMb9vG7l89QnWUkQsTWylzpZAQ5eMS7V6Izkaz4mWInfK53vZ09+EriqIox5r6Zbh5D1jPg7d/BP+8jrDz8rn0xVe4QBNBXUkEyyMjCJpxHa+3zMDZ04Tvr0uhYs1pK0kFvqIoyukSZoFvvQnLfgW1m+DPi9Dbypn14ktc8KVr6HmtiIQX/8W8bzzAa71L6bLpYNV34PD6Tz72SVCBryiKcjppNLDoVvjhVjBHwUv/gdj8AHG334z1mWfw9PbivPteLr3ganaZvsG2rnQ276jB43af+lJO+REVRVGUf5cwA258D867CfY+Bc9cQsiUSNJXv4l5+nSGHv4lCx1axMJb6W1rOy1346rAVxRFOVP0Zrjyf+G612C0B56+GH3DG6Q89yzRN/8nI2vWkPyvNVy14vtodbpT/vYq8BVFUc607Mvh5t2QeQlsuBfx6grirv86Kc+uxDswyMAzK0/L26phmYqiKIEiJRQ8B+/eD4Yg+PrLeIKyEAYD2vDwkzqkGpapKIoyEQkB838AN+2AhHyISkcXG3vSYf9JTn0nkaIoivLZxE6B76w+7W+jWviKoiiThAp8RVGUSUIFvqIoyiShAl9RFGWSUIGvKIoySajAVxRFmSRU4CuKokwSKvAVRVEmiQk9tYIQogdoOsk/HgP0nsJyTqWJWttErQtUbSdjotYFqraT9WlqS5VSxh7vGxM68D8PIUTBieaTCLSJWttErQtUbSdjotYFqraT9XlrU106iqIok4QKfEVRlEniXA78pwNdwMeYqLVN1LpA1XYyJmpdoGo7WZ+rtnO2D19RFEX5qHO5ha8oiqIcRQW+oijKJHFOB74Q4pdCiFIhRLEQYqMQIjHQNQEIIR4TQhwer+1NIUREoGs6QghxrRCiXAjhE0IEfGiaEGK5EKJKCFErhLg30PUcTQjxnBCiWwhRFuhajiaEsAoh3hNCVI7/v/xRoGs6QghhEkLsF0KUjNf2cKBrOpoQQiuEKBJCrA10LUcTQjQKIQ6NZ9lJr/t6Tgc+8JiUMl9KOQtYCzwY6ILGbQLypJT5QDVwX4DrOVoZ8B/AjkAXIoTQAn8CvgBMA74hhJgW2Ko+4u/A8kAXcRwe4MdSyqnAQuDWCfS5OYFLpJQzgVnAciHEwgDXdLQfAZWBLuIELpZSzlLj8E9ASjl81GYwMCGuUEspN0opPeObe4HkQNZzNCllpZSyKtB1jDsPqJVS1kspXcA/ga8EuKYPSCl3AP2BruNYUsoOKWXh+OsR/AGWFNiq/KSfbXxTP/6YEP8uhRDJwBeBlYGu5XQ5pwMfQAjxqBCiBbiOidPCP9r1wDuBLmKCSgJajtpuZYIE19lCCJEGzAb2BbaSD413mxQD3cAmKeVEqe2PwE8BX6ALOQ4JbBRCHBRC3HiyBznrA18IsVkIUXacx1cApJT3SymtwMvAbROlrvF97sf/6/fLZ6quT1vbBCGO87UJ0Ro8GwghQoDXgTuO+W03oKSU3vFu1mTgPCFEXqBrEkJ8CeiWUh4MdC0nsFhKOQd/9+atQogLT+YgulNb05knpbzsU+76D2Ad8IvTWM4HPqkuIcR3gS8Bl8ozfDPEZ/jMAq0VsB61nQy0B6iWs4oQQo8/7F+WUr4R6HqOR0o5KITYhv86SKAvfC8GrhJCXAmYgDAhxEtSym8FuC4ApJTt48/dQog38Xd3fubrbGd9C//jCCGyj9q8CjgcqFqOJoRYDtwDXCWlHAt0PRPYASBbCJEuhDAAK4C3AlzThCeEEMCzQKWU8veBrudoQojYI6PShBBm4DImwL9LKeV9UspkKWUa/r9nWydK2AshgoUQoUdeA1dwkj8gz+nAB3493lVRiv9DmijD054EQoFN48Os/hLogo4QQlwthGgFFgHrhBDvBqqW8QvbtwHv4r/wuEpKWR6oeo4lhHgF2APkCCFahRA/CHRN4xYD3wYuGf/7VTzecp0ILMB74/8mD+Dvw59QQyAnoHjgfSFECbAfWCel3HAyB1JTKyiKokwS53oLX1EURRmnAl9RFGWSUIGvKIoySajAVxRFmSRU4CuKokwSKvAVRVEmCRX4iqIok8T/B+/zuSfdvmnVAAAAAElFTkSuQmCC\n",
"text/plain": "<Figure size 432x288 with 1 Axes>"
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.6.10",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "numerical integration of ode in Python.ipynb",
"public": false
}
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment