Skip to content

Instantly share code, notes, and snippets.

@shivin9
Last active July 28, 2017 08:48
Show Gist options
  • Save shivin9/ed786b6f97965824ab89e5af5c99b31e to your computer and use it in GitHub Desktop.
Save shivin9/ed786b6f97965824ab89e5af5c99b31e to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"using DifferentialEquations, Plots, BenchmarkTools"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"pyplot()\n",
"using Revise"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[1m\u001b[36mINFO: \u001b[39m\u001b[22m\u001b[36mRecompiling stale cache file /home/shivin/.julia/lib/v0.6/PDEOperators.ji for module PDEOperators.\n",
"\u001b[39m"
]
}
],
"source": [
"using PDEOperators"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"advection_left_dirichlet (generic function with 1 method)"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"function advection_left_dirichlet(t, u, du)\n",
" h_inv = length(u) - 1\n",
"# du[1] = 1.5*u[1] - 2*u[2] + 0.5*u[3]\n",
"# du[2] = +0.5*u[1] + 0*u[2] - 0.5*u[3]\n",
" du[1] = -1.5*u[1]\n",
" du[2] = +2.0*u[1] - 1.5*u[2]\n",
" @inbounds for i in 3 : length(u)\n",
" du[i] = -0.5*u[i-2] + 2.0*u[i-1] - 1.5*u[i]\n",
" end\n",
" scale!(du, h_inv)\n",
"end"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"N,M = 100,5000\n",
"Δx = 1/(N-1)\n",
"Δt = 1/(M-1)\n",
"a = 1\n",
"\n",
"x = collect(0 : Δx : 1)\n",
"ϕ(x,t,a) = sin.((π*x-a*t)) # solution of the single forward moving wave, a is the speed of wave\n",
"u0 = ϕ(x,0,a);\n",
"\n",
"A = UpwindOperator{Float64}(1,2,Δx,length(x),true.|BitVector(length(x)),:Dirichlet0,:nothing);"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [],
"source": [
"single_solition1 = ODEProblem(advection_left_dirichlet, u0, (0.,5.));\n",
"soln1 = solve(single_solition1,SSPRK22(),dense=false,saveat=0.001,dt=Δt);\n",
"\n",
"single_solition2 = ODEProblem(A, u0, (0.,5.));\n",
"soln2 = solve(single_solition2,SSPRK22(),dense=false,saveat=0.001,dt=Δt);"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAYAAAByNR6YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzsnXd4FFXbh+/Z2fTe6L0X6R3pvYMKdhBFUREQoiBSRQUEaSIIioqgfiooYkFAqrTQBSkCAgLBENLLJtns7sz5/shLXn2lbGCTmd3MfV25roWd3fmd/e2ZefaU55GEEAIDAwMDAwMDAwOXYdJagIGBgYGBgYGBp2EEWAYGBgYGBgYGLsastYDrqKpKXFwcQUFBSJKktRwDAwMDAwMDg38hhCAzM5MyZcpgMt18nEo3AVZcXBzly5fXWoaBgYGBgcEtkWQJv9KBBFYMIaBCMAHlg/EvE4hPmC+S/N8briPHjj3ThmJVUKwO1FwHjhwHql1FMoEkm5BkCZNsQpJNyH5mvIK98Q72Qfb95+05NyWH7KsWcuIsZMdZyI7LJPNCGkqOo6ibb/AfYmNjKVeu3E2f102AFRQUBMC5c+eIiorSWI1BQbBarUyfPp1p06bh6+urtRyDAqAn7zJtgj3xgpgEwa/JgqPJgnRb3nPlA6BmqETlIInKQVApMO9x+UAI9OKOR71VIUjKgUsWwYVMuJgp+DNT8GcmnEwVZNrzjqsYCA0jJBpHSrQtLdEwXEI2aTvSrifvhOLAFnuO3AsnsP15Cnvcn6iWdAAkX3+8SlfCXLI85rASyKFRyGFRmMOiMAWGIt1iBOCW51QV1OxM1Mw0HMnXcCTFYU+KQ0m6iiM5DjUjDQBTYAje5avjVb463hWq412hFiafgn1emTYLJ5JPcyL5NCdTznAp8woONS+wKeEXSeXgClQIKkuUXyRRfuFE+UUQ6RdBgJf/Td8z379JN/bP6rCSZssgzZpOXNY1/sq6yhXLVf6yXCUu6xp21Y6ERIWgstQMq0at//yVDSyNSTJW/xQmiYmJVKtWLT9uuRm6CbCuXyCDgoIIDg7WWI1BQfD29sbHx4fg4GDNL/QGBUNL7xyq4GCiYPNfgs1XBPsSBA4BZQOgWaTEuGoSTSMlmkRJRPoWXjATGgLVgM7/8/+qEJzPgEOJgkNJgsNJgjlnBdNOQbgPdCoj0bWsiS5lJaoEF32wpaV3QgjsV/8k9/QRcs8dI/f8CURuDpKPP8FV6+LVaQDeZavgVbYqcnjJwlv2ERoGVLjhU0pmGrZLp/P+Lp/Ftv8nbL9kYZO98Kl6D761m+JbqwnmUhX/pS9XsXEk/jcOXT3K0YQTnE+9iEBQOrAkDUrUpW/1blQJq0TV0EoEegfckfTb+RdMMCUoceO2qQpXMq9yKukMJ5POcDLpNJsv/4JAEO4bSvMyjWlRpgnNSjUkyCfwjvQZ3Byr1Qrc/oedpJc0DRkZGYSEhJCenm4EWAYGHopNEWz5S7D6gsp3lwRpNgjxvh6sSHQpa6Ja8J2PSBU2dlWwP+G/QeGBRIEioEYIPFjFxINVTNwTpl/9d4NQVWyXz5Dz2x5yju1BSb6K5OWDd5W6+FRvgE+1BniXr44ky1pLvSFCVXEkxGI98yu5pw9jPXcM7Dbk0Ch8azdFrdecX31z2fXXAfbHHSbHYaVUQAkalaxHw5L30LDEPZQKvHHAowey7dmcSjrL4fhj7Is7zIW0S8iSibqRtWhZtgkdKtxL2aDSWsv0CJyNV3QXYMXHx1OyZMl/PZ+ZmcnVq1dRVVUDda5DlmUqVqyIt7e31lJcRk5ODqNGjeLdd9/Fz89PazkGBaAovLOrgq3/Caq+vZgXVNX8T0DSq7xE0ygJs8bTbXdKuk2wPU7w3aW8tqXboFbof4Ktyibqhhdeu4rCOyEEtktnyDm8jZzf9qKkJ2EKDMWvXiv86t+LT/X6SGb3vJYJWy4Z535lx8mf2JZ2iuO+NhSTRFUpmLblW9GhXh8qhZQvtGC5sP1LyEpkX9wR9scd5tDVo1iVXOpE1qRrpfZ0rHgvYb6hLj9nceHatWuUKlXqtgGWbqYIr/O/K/JVVWXWrFl8++23GilyPf7+/nzxxReULVtWaykuwWQyUa5cuVvupjDQJ4Xp3fkMwfLTKivOqiTkQPVgGFnXs0Z5QrwlBlSSGFDJxPtt8ka2Vl9QeeeEyutHVJpESjxb28QjVSUCvVzb3sL0TrGkk31oG1n7NuKIv4QcEolfg3vxa9AG78p1kEz6HKVyBiEEZ1LO8dP5rWy9uBOLPYt6FWvzbFB1GsdlEnBsP+qB75F3HCSjaScCW/VEDo10uY7Cvm6WCIiiX/Xu9KveHasjlz1X9rP54k4WH/6IxYc/pGnpRnSr3IH25VvhJXsVigZPxVnPdDeC9b8R4YwZM1i3bh2jRo2iUaNGeHm59xfBarUyZcoUSpcuzbJly4ygxMCjsKuC7y8Jlv2usuUvQag3DKlu4smaJhqEe0ZQ5Qy5imBjbF6A+VOsINALHqtm4tnaJhpG6PMzEEKQe/ZXsmI2knM8BhD41WtFQIvu+NRs5NZBFeRNof10fivrz2/mQtolovwi6F6lIz2rdKZccJn844SqkHvuONlHtpNz5BeEw4ZfvdYEtOmLT7X6bv8dTrNm8MvlPWy+uJPjiacI9w2lT7Vu9Kvegyj/CK3luQVuO0V49epVSpUqlf9/nTp1YvTo0QwZMkRjha5j06ZNTJo0iQ0bNnjEjsns7GyeeOIJVq5cib//zXfNGOgPV3mXmitYclJlySmV+BxoVSJv5GZQFQl/s3vfkO6WyxbBh6dVPjqjEpcNLUpIjKtvYkDFu9uJ6CrvhMNO9pEdZG7/BsfVi5hLVSCgZQ/8m3ZCDnT/aaSErES+OfMjP5z7GasjlzblW9C7alealmqAfJugUbVmkX1wK5Y9P+KIv4y5RHkC2/TBv0U3TD53N62nh+vmxfTLrDu7gY0XtmFTbLQt34r7avSiQYm6bh9IFibx8fGULl3a/aYI5b8tkIyPjwegUaNGWskpFK7nzUhNTfWIAEuWZVq1avUP7wzcg7v17q8swYLjKu+fVnGoMLSGiRF1TNQrxLVH7kaFQInXm8pMbWzix8uCd06oDNyiUCMExteXeby6hI9c8M/rbr1Tc7LIivkJyy/foaQn4VunOaH3P+8RozQAZ1PO89Xv69h+aQ9+Zl/6Ve/B/TV7U8Lf+ek+k28AgW37EdCmL7bzx7Hs/pG0dR+QsfEzAtsNILBdP0z+t96qfzP0cN2sFFKBMc2e5ZmGg9l0YTvfnl3Pi1smUTO8Gk/Ue5DWZZt7xHfB1Tjrme5GsP4eEZ4+fZrHH3+czz77jFq1amms0HV4arsMig9n0wRzflNY9YfA3wwv1DExuq6Jkv7GxdgZ9ieozD6msu6ioLQ/RNczMbyWiSDvwv/81KxMMretwbL7R4TDhn+TTgR1vB+v0pUK/dxFwe9JZ/not885ePUopQJKMKhWP3pV7Yz/LXJSFQRHaiKWbWuw7NuIZDIT2KYPgR3uRw5y/9E+IQQHr/7KZye/5ljCSaqGVmLwPQ/SvkIrI7fW33B2ilB3n1hWVpbWEu6Iixcv0qFDB0JCQmjYsKHWcoqUrKwsunfv7rbeFWcK6l2sRfDULw5qf+3gp1jBzGYmLj9iZkYz2QiuCkCLEibWdjVzapCZ7uUkXj2oUuUrB++cUMhVnPvNW1DvVGs2GRs/5+obT2DZ9T0B9/am9NSVhD8a7RHB1bnUP5n4ywye2zSOhKwkpt37Mp/3W8bAWn1dFlwBmMOiCH1gBKWnriSwTR8su3/g6utDSFu7DMWS5vT76PG6KUkSzcs0ZlHXmSzqMpNQ3xBe2z2HoT+OYvOfO1BURWuJusBZz3Q3Reiui9iDg4N58803SU9PZ9KkSVrLKVK8vLwYNGiQ23pXnHHWu2SrYNZRlcWnVIK9YGFLE8Nrm+5oasvgv9QKlfi4vZnpTQTTjyhE71NZcFzl9SYyj1W79RotZ70Ttlwse34kc8tq1NxsAu/tQ1CXhzxixAXgYnosK377gh2X91A2qDSTW4+lU8W2t11fdbfIQWGE9H2KoM6DsOz8nswd35C1/2eCujxIYPsBmLxvnfxV79fNBiXrMr/k65xMOsOnx1fz5t4FfPn7Op5vNJSmpYvXIML/4qxnxhRhAZk7dy5nz57lgw8+ACAtLY1q1apx9uxZwsPD2bFjB2PGjOHo0aM3fQ89tsvA4EZk2QULT6jMOaaiAuPqmxh7T9FMZRVHfk8VTD6ksPai4J4wmNlMpk8F6Y7WwQhVJfvQNtLXr0DNTCWgRXeCuj2KOcz9131C3m64j3/7P344t4kovwieqP8w3St3xKzRbkfFkk7m5i+x7P4BU2AIIT0H49+8q9vvvrzOicTTLD2yghNJp2leujHPNx5KldCKWsvSBGenCHU3gpWVlXXbTO4XMvKSFbqaUG9uW/Li6aefpkaNGsyZM4fQ0FBWrFhB//79CQ8Pd70gNyErK4suXbqwZcsWAgLurGyEgTbczDshBF//KYjep5CQAyPqmJjY0ESUnxFYFSa1wyS+6Wpmf4LKhAMq/X5W6FpW4t3WMjVD//nZ36rf2S6fIe2bpdguncavYTtCeg/FHFUGT8ChOlh3dgMrjn8BAp5vNJQBNXrhrXEuJzkwhND7niWwbT/S139C6pcLydzxLaEDhuNbq8m/jne36+Y9UbVY3O0tdsbG8P7RVQz7aQw9q3RiWIPHiPArXvc/t50ivF2G8ySroPpqB2ohjLvJEsQ/br5l3bPQ0FAGDhzIxx9/zNixY1m6dClfffWV68W4Ed7e3kRHR3tUdvriwo28+z1VMGqvwtY4Qf+KEgtaylTWoNZecaZFCRPbekusvyx4MUah3jcOouuZmNzIlJ+w9EbeKZlppK9fQfb+n/EqXYmokXPwqVZfq2a4nINXj7L48IdcSr9Cn2rdeLrBY4T6hmgt6x+YI0sT8cSr2Do+QNp3y0laNgm/hm0JHfDsPxKWuuN1U5Ik2ldoTeuyzfjuj42sOvEV2y/t4emGjzOges9Cn5bVC856prsA63Zzm5G+En88aC60ESxnisqOHj2afv36Ubt2baKiojwujURBub6WwMD9+Lt3FrvgjSMq84+rVAyC9d1lelXQ3T6YYoMkSfSpKNGlrMSc31RmHVX5/JzK/JYyAytL//BOqAqWXT+QsWEVmEyEPjCCgFa9dFsXsKAkZSfzzqHl7IyNoX5UHT7oOZ8a4VW0lnVLvCvUIGrkHHIObydt3QfEz3qG4J6DCWzbH0mW3fq66SV7MbBWX7pX7sgHxz7l3UMfsunCNqKbP0+tiOpayyt0nF2DpbsAy2Kx3HaKUIvK9X+nVq1aVKlSheHDhzNnzhxNtegBi8VCixYt2L9/P4GBRuV2d+K6d9PWHiT6sBcpVpjW2MTL9U34FvMEoXrB1ywxtbHM4Gomxu5TeHCrQreyEvMbW3mwcwv2fv811u+WYY89S0CrXgT3fgI54NbXUHdBFSo/ntvMsl8/wUf2Yeq9L9GpYlu3yc0kSRL+TTvhW6c56T+tJP275WQf2EzooFHYoyq4/XUzyCeQl5o/T48qnZh/YCnPbRzHgBq9eLrBYwR663/a806xWCxOHae7n6e+vrfeeaEXnnnmGRwOBwMHDgTysvKWK1eOQYMGcerUKcqVK8err76qscqiwdfXl/nz57uNdwb/JVvyodTEjTy004sG4RK/DzIzubFsBFc6pHKwxLpuZtZ3l/k9TdB2vYnJw4aR+d44hDWbqNHzCHtwlMcEV1cy4hi7dQrzDrxH+wqtWdVnMZ0rtXOb4OrvmPwDCRv4AiXGvgOymcR3orFtXMmCt2d7xHWzbmRN3u8xjxGNn2Tjha0M/uEFdsXu01pWoeGsZ7obwTKbdSfphmzfvp0RI0bkDxX6+/tz5coVjVVpg9lspnv37lrLMCgg6y6qPLcbcpVSrGwvM7j6ne1WMyhaelUwcazFWS6sXEh4RizfVRpIr8GPUirCR2tpLsGhKqz+fR0rjn9JpF848zu/TpNSDbSW5RK8K9SgxNiFWHZ+R/r6T6gfGoVSuwrmynW0lnbXmE0yD9buT4cKrZl/8H0m75xF98odGdX0aYK83XOE7mY4G6fobgQrMzNTawm3JC4ujlq1anHkyBHGjBmjtRxdkJmZSbly5XTvnUEeSVbBo9sc3LdZoXGYgu+sFtxXOssIrtwAoThI/3EFWYujKR0o8fgpP+ZEPUa9700sOamg6iPrzh3zV+ZVRm1+leXHPuO+Gr1Y0XuRxwRX15FMMkEd7ifw+Tn8evoPEhe9RNr3HyLshbCwWANKBEQxq/0kJrQcze4r+3ly/WgOXr152iJ3xNl7ne6Gi/z87q6AZmFTpkwZTp8+rbUMXeHn58eaNWt0750BbL6iMmSHQq4Kn3WUebAiHFqxxPDODbAn/kXKqrew/3WB4F5D8Gt3H3PbHqFWQ4nJR0yM3Kuy7pJgVQeZ0m6WVV8IwcYL23jn0AeE+YbybtdZ3BPl2TkCAytUI/S5WQRmxWLZ9DnWkwcIf+xlvCvU0FraXSNJEj2rdqZxqfq8FbOIl7dNY0D1njzXeCh+ZvefEnX2eqm7AMtdpggN/ovZbKZVq1ZayzC4BXZVMPmgypzfVLqWlViZfxM2Gd7pHCEE2ft/Jm3tUuTgcEqMmY93hZoA+d4tuRfuqyQxeLtC/W8crGzvPjtAM3IzmXvgPX65vJceVTrxYtNnXFraRq+YzWZa3XsvAAH1WpLy+VwSFo4lpO9TBLa/D8nkHv7dipIBUczrPJ3vzm5g6a+fcDD+KK+1Ga/7HaC3w22nCDMyMrSWYFBAMjIyCA4ONrzTKeczBG2+V5h/XGVOcxMbe/53hMPwTt+oWZmkfDKD1C8X4NeoHSXGLckPrv7Xuy5lTfz2gJnmURK9NylExzhf11ArjsT/xlPrX+RI/G+81mY8r7Z6sVgEV/BP/7xKV6LEmAUEth9A+nfLSf5wGoolXWuJLsEkmbivZm8+7vUOAV7+jNg0jm/PrEcnRWTuCGevl7oLsNwho63BPwkICCAmJsbwTof83zmVRmsdJOcK9vSTGddAxvS3tVaGd/ol989TXHv7eaxnjxI+dBLhj0Rj8vnv1MSNvIvyk/ixu8zCViaWnFJp9Z2Ds2n6u5GpQuWT418SvXUq5YLL8HHvd+hY8V6tZRUp/+ufZPYitP8zRAx/Hduls1x7ewS5537TWKXrKBdchiXdZtO3encWHvqAabtmk2lzLt2B3nD2eqm7AEv2kMR4xQlZlqlbt67hnY7IVQTDdzl4bLtCv4oSR+4z07zEv7u74Z3+EEJg2fU9iYvHI4eVoOT4pfg3bPuv427mnSRJvHiPzL7+ZrIc0PhbB6vPq0Ul/7ak52bwyvY3+OS3Lxla72HmdZpOCf/I27/Qw7iZf351mlNy3BLMkWVIXDKBjI2fIVRFI5WuxVv24sWmw3mj7QQOx//G0z+N5feks1rLKjDOXi91F2C561TFtm3baN68OXXq1KFu3bqMHz8eVdXPRa0wycjIQJIkt/XO04i1CNr9oLDqD8FH7WQ+62gm+CbFmQ3v9IWaayXlszmkffMegff2IWrknJsWZ76dd40iJQ7fZ6ZfRYmHtimM26/gKIwaYwXg9+Q/eOanaM6k/MGcjlMZWv/hYlNe5X+5lX9yaCRRL7xFULdHyNj0fyR9MBU1y3N2aber0IoPey0gzDeEkZtf5ZszP7rVlKHbThG6a0bbsLAwvvzyS06dOsXhw4fZu3cvq1at0lpWkRAYGEhsbKzbeudJbI9TafKtg/gcwZ6+Zp6qeesubninH+wJV0hYOAbr8RjCh0wg9P7nkOSbL6Z1xrtAL4nPO8osaGliwXGVbj8pJOYU/Y1MCMF3Zzcw6ucJhPuFsrznApqXaVzkOvTE7fyTTDIhPQcT+dyb2C+f5dr80djj/ixilYVH6cCSvNt1FvfV6MmiQ8uZve9dbIpda1lO4ez1UncBlt5z8cydO5fhw4fn/zstLY3IyEgqVqxIlSp5OyN8fX1p2LAhFy9e1Ehl0SJJEsHBwbr3zpMRQjD/N4WuPynUD5c4NMBMk6jb+2F4pw9yjseQMH80KA5KRL+Df+MOt32Ns95JksSYejJbe8ucTBM0+dbBwcSiG123KXbe2reI+QeX0adaNxZ1nUXJgBuPyhUnnPXPt2ZjSrz0LiYfXxIWjiX76K4iUlj4eMlejGzyNK+2epGtF3fy4uaJJGUnay3rtjh7vdRdgKX3ZJVPP/0069atIy0tDYAVK1bQv39/wsPD84+Jj4/n66+/pk+fPlrJLFIyMzMJCQnRvXeeSrZD8Oh2hZf2q7xUL2+XYJSfcxcAwzttEUKQsflLkj+anncjjX4Hr1IVnXptQb1rX9rE4QFmygRItPleYcWZwg+yUq1pRG+dwraLu5jUeixjmj2Lt+xcoVxPpyD+mSNKEfXiAnzrtiDlkxmkr//EY9ZlAfSo0ol3u80iMSeZ4Rtf4mTSGa0l3RK3TTQaFBR022McSVdRc1y/+8DkF4g5svQtjwkNDWXgwIF8/PHHjB07lqVLl/LVV1/lP5+RkUHfvn0ZP348TZs2dblGPRIUFER6erpT3hm4lvhsQb+fFU6mClZ3lhlUpWC/mQzvtEM4bKR+tYjsg1sI6v4Ywd0fK1DuozvxrlygxC99ZEbuUXhqp8KZdMHMZqZ/7Cx1FedTLzLxlxnYFBsLu86gbmRNl5/DnSmofyYfX8KHTMBSvhrpP6zA/td5wodMwOTrGTuAa0VU54Me85m66y1e3DyRsc2eo3e1rlrLuiHOeqa7AOt2C90USzrxM4aBKIRfXyYTpV//Ajkw5JaHjR49mn79+lG7dm2ioqJo1KgRkBfV9ujRg/79+xMdHe16fTpFCEFGRgaBgYHGVFMRcjxF0GeTA4cKu/qaaRxZ8M/e8E4bFEsayR+9ji32D8IHv4J/k44Ffo879c5HlvigrUztUJWX96v8kS74tKOMvwsLfO+5coA39syjbFBpZrafaUwJ3oA78U+SJII6DcKrTBWSP5lJ4jsvETH8jZtuhHA3wv1CWdD5DRYdWs6c/Yu5nPEXzzYagknS12SbswvydRdgWSwWQkNDb/q8HBhCqUkfFdoI1u2CK4BatWpRpUoVhg8fzpw5c4A83T169KBHjx5MnjzZ5dr0jMVioXz58qSnpxMcHKy1nGLBhliVh7YqVAmCH7ubKRd4ZzdHw7uix371IknLpyHsNqJGzsGnUu07ep+78U6SJKLry1QNlnh0u0L7HxS+7373JXaEEHz5+7e8/+sq7i3XnEmtx+LvZZRhuhF3459vrSaUGDOfpPenkLBwDJHPTMe7XLVCUlq0eMlevNRiBBVDyrP48EfEZyUwsfUYfGRvraXlY7E4F3/oLsBy5ot2u2m8ouCZZ55h5MiRDBw4EIB33nmHAwcOkJWVxdq1awEYNGgQkyZN0lJmkRAcHOxWW2zdnSUnFUbHqPQqL/FFJ5lArzu/KRreFS3W04dJXjEDc0RJIka9jTmsxB2/lyu861/JxK6+En03OWixzsGP3c3Uj7iz75OiKiw8+D7fn9vE43UHMqzBY7obedATd+ufV6mKlBi7kKTl00hc9DLhT7yKX90WLlSoLQNr9aVkQBRv7JlH9JYpzGg/iVBfffwIdDYg1t23X1HcY+He9u3bGTFiBF5eeQs2J02ahN1u5+jRo/l/xSG4gjzPTp486TbeuSuqEETHKIzcqzK6rol1Xe8uuALDu6Ik+9A2kj6Yik/VukSNnndXwRW4zrvGkRIHBpiJ9IV7f3Dw85WCL7/IdeQyddds1p/fzPgWI3mm4WAjuLoNrvBPDg4nauTb+NRoRPKH07Hs+t6FCrWnbfmWLOzyJlcyr/LCz69wJfOq1pIA5+MU3fWArKwsrSXckri4OGrVqsWRI0cYM2aM1nJ0QVZWFq1atdK9d+6MTREM3q6w8ITK4tYmFrSSkU13v2bG8K5oyNz+DSmfzcG/aScihk3D5Hv39fZc6V3ZAImdfc20KyXRZ5PClwXI/J6Rm0n01qkcvPorM9pP0u3CZL3hKv9MPr5EPDWZwHb9SfvmPdK+W+5Ro9J1ImuytPscTJLEiE3jOZF4WmtJTnvmllOEWlKmTBlOn9beYD1hFAsuXLLsgoFbFLbGCb66g52Ct8LwrnARqkr6Dx9h2f4NQV0eIrj3UJdtJnC1d4FeEuu6yQz7ReHRbQpJVsHIurfOsh5vSWDc9umk56azsMub1DF2CjqNK/2TTDKh9z2LHF6S9G+XoWZlEvbQi0geUgKrTFAplnSbzaRfZhK9dQqvt51Ay7JNNNPjtlOEDodDawkGBcThcBATE2N4VwgkWwWdf1LYfU2woYdrgyswvCtMhMNO6udvY9mxltD7nyekz5Mu3alZGN55mSQ+6SATXc/EqL0qUw8pNx0NOZ96kRd+fgWbYmNJt9lGcFVACsO/oPYDCH98PNmHtpD8yQyE3eay99aaYJ8g5nZ6jSalGjDxlxlsu6hdwlVnPdNdgJWTk6O1BIMCkpOTw6BBgwzvXEysRdD2BwcXMgQ7epvpXNb13dXwrnBQbVaSPnyN7KO7CR/yKoHt+rv8HIXlnUmSmNtSZk5zE2/8qvL8bhXlf2oYnkw6w+jNEwnzDeW97nMoH1zWpRqKA4Xl3/VpaOvpQyR9MAXVmu3S99cSH7MPb7SbQKeKbXl9zzy+/2OjJjqc9Ux3U4RGwkP3IygoiCtXrmgtw6M4kybo+pMDkwS7+5qpEVo4OaoM71yPas0m6YOp2K+cI/LZ1/Gt0ahQzlPY3o1rIBPpK/HMLoXkXMHnHWW8ZYlfrx3n1R1vUi2sCm91mEygt2ckuixqCtM/v7otiHpuJknLp5K45BUin30DOfDm6Y/cCbPJzMTWLxLkE8C8A0vJyLXwWN0HijSPn9PJYQtZR4ExpircD4fDwaZNmwzvXMTJFEH7Hx0EecHefoUXXIHhnatRszNJWjoRe9wFIp+fWWjBFRSNd0/WNLG2q8z3lwQPbFHYFXuY8dsdE414AAAgAElEQVRfp05kTd7uNM0Iru6CwvbPp+o9RI18GyU1kcRFL6OkJRXKebTAJJkY3eQZhtZ7mOXHPmXZryuLdGG/204RWq1WrSUYFBCr1Up0dLThnQs4lizosN5BST/Y0SevblxhYnjnOhRLGolLJuBI/IuoF2bjU7lOoZ6vqLzrV9HE991kDl/dx+SdM2lUsgGzOkzGz+xbqOf1dIrCP+9yVYl6cR7CbiNh8TgcqYmFdq6iRpIknqz/CKOaPM2Xv3/LgoPvoxZGhZcb4KxnuguwAgMDtZZwR8TExNCwYUMaNmxI3bp1efbZZ8nNzdVaVpEQGBjIyZMn3dY7vXA4UdBxvYMKgbCtt9npgs13g+Gda1DSk0lcPB4lPZmokXPwLl+90M9ZlN6ZHLuo4TWPFKU5By0v41CNgs13S1H55xVVlqiRc0BVSVw8DkfKtUI9X1EzsFZfxrV4ge//2MiCA8uKJMhy1jPdBVh2u11rCXdEgwYNOHjwIEePHuX48eMkJCTw3nvvaS2rSLDb7axZs8ZtvdMD+66pdP7JQY0Qia29zET4Fs16AsO7u8eRmkjiu+MQOdlEjXobrzKVi+S8ReXdhvNbmbF3Ad2rdGBxl2gOJcv03KiQafOcXEtaUJR9zxxRKi/IQiLx3XE4kuML/ZxFSZ9q3RjXciQ/nPuZeQeWFnqQ5axnuguwbDZ9byudO3cuw4cPz/93WloakZGRWK3W/KzuNpuNnJycYlM812azMX/+fN17p1d2x6t026BQL0zi554yoT5F970xvLs7HGmJJC4ej1AcecFVyfJFdu6i8G7The3M3vcufat145WWo2hfxoufe8ocSxZ026CQlmsEWXdKUfc9c3hJokbOQZLNeUFWUlyRnLeo6F21CxNajWb9uc3M3b+kUIMsZz3T3S7CgIDbL5qMy4zHYnd95ulArwDKBJW65TFPP/00NWrUYM6cOYSGhrJixQr69+9PeHg4Fy9epH///pw/f57evXszYsQIl2vUIwEBAcTExGgtwy3Ze02lxwaFZlESP3S/+9I3BcXw7s5R0pNJWjIBVIWokXMwR9z62uFqCtu7zX/+wlv7FtGrahfGNn8uv/RNq5ImtvaGbhsUemxU2NxTJsi7ePyYdCVa9D1zWBRRI+eQuOQVEt8dT+TI2XhFeU6KjR5VOiEhMSvmHRShMr7FSGST65OtOhOngA4DrNtFhmnWDB774flCiU5lycTa+1fesqBkaGgoAwcO5OOPP2bs2LEsXbqUr776CoBKlSpx7NgxLBYLjz/+OGvXruXhhx92uU69YbPZWLVqFUOGDMHbWz8Vz/XOoUSVnhsUmkRK/NhdJqCIgyswvLtTlMxUEpdMQNhyiRr1dpEHV1C43m27uIuZMQvpVrkDL7cY8a+6gk2jTGzuKdFpvYPemxQ29NDm++vOaNX35NDI/wRZE0hcPJ4So97GHFmmyM5f2HSv0hGTJDEz5h0AXmk5yuV1Md12BOt2c5uhvsF83ndpoY1gOVOte/To0fTr14/atWsTFRVFo0b/3IodGBjIww8/zOeff14sAqzrawkeeeQR4ybtJNenWOqEaRdcgeHdnXB9t6BqzcoLriJLa6KjsLz75fJe3tw7ny6V2jG+xcib3pyaREls7CnTbYNCv58Vfuwu42c2gixn0bLvySERRI2cTeK740hcMuE/xcejilRDYdK1cgckycSMvQvwNnkR3fx5ly7ZcXYNlu4CLGeG3m43jVfY1KpViypVqjB8+HDmzJkDwLlz56hYsSJeXl7YbDa+/fZb6tevr6nOoiIgIIBNmzZpLcNtOJUq6PKTgypBEht6aDu9YnhXMJSsDJLeexU1K4OokXM0nV4pDO92xe5j+u65dKzYhgktR992eqVVSRPru0OPDQoPbFH4tquMj2wEWc6gdd+Tg8OJHPEWiYteIum9CUSNehs5OFwzPa6mS6V22BQbs/e9i6/ZlxGNXVeqytkpQt0tcneX1AbPPPMMDoeDgQMHArBt2zYaNWpEgwYNaNSoESVLlmTKlCkaqywacnNzmT9/vtt4pyVn0wSd1zso4w+binhB+40wvHMeNdtC0tJXUTJSiBoxq0gXtN8IV3t3IO4Ir+1+m3blW/JqqzFOr11pV9rE991ltsUJHtqqYFeNhe/OoIe+Zw6LIuqF2ag2K4lLJ6JkeVbh915Vu/Bi0+GsPv0dnxz/0mXv66xnuguwFEXRWoJTbN++nREjRuTvHBw+fDgnTpzg2LFjnDx5kkWLFuHrWzwS8SmKQkxMjNt4pxV/Zgg6/+QgzAc2F2EqhltheOccqs1K0vKpOFKuETliFl6lK2ktyaXeHU/4nck7Z9G8dCMm3xuNuYALg7uUNbG2i8xPsYLHtys4jCDrtuil75kjSxM14i3UzFSSlk1CzXH98hstub9mb4Y3HMwnx7/ki1NrXfKeznqmuwDL399fawm3JC4ujlq1anHkyBHGjBmjtRxd4O/vz5o1a3TvnZbEZwu6bnDgK8PW3mZKFEESUWcwvLs9QnGQ8skM7H9dIHL4G3iXqaK1JMB13p1L/ZMJO96gdkQNXmszDrPpzlaO9KpgYnVnmW/+FDy/WynS0iXuiJ76nlfJ8kQ+PwtH8tW8AtG5nlX8/bG6AxlcdxDLfl3JurMb7vr9nPVMdwGW3qcqypQpw+nTp9m7d69RmPo/5Obm8tprr+neO61Itwl6bnSQ7cgbuSrtr4/gCgzvbodQVVL+bx7WM78SMWwqPpVqay0pH1d4dyUjjnHbXqNMUClmdpiEj9nnrjQNqGTio3YyH54RTD5UNGVL3BW99T3vslWIenYG9riLJH/8BsLhWcmHhzV4jIE1+7Lg4DI2Xdh+V+/ltlOEqmp0SndDVVWuXLlieHcDrA5B/58VLmbCzz3NVArST3AFhne3QghB2tql5BzZQfjg8fjWbKy1pH9wt94lZCfx0ra8gs1vd5xGgJdrRlKeqGFibgsTM4+qLDxuTD3fDD32Pe+KNYl4ehq5546T8sV8hI603S2SJDGyyTB6Ve3C7H2L2PfX4Tt+L2c9012A5efnl/9YlvPWAXhaIdrrWzxNJt19/HeEn58fH3744T+8MwCHKnhkm8KBBMGP3WXuCddXcAWGd7ciY+NnZO3+gdAHR+PfsJ3Wcv7F3XiXZs3g5a3TEEIwr9N0Qn1DXKrtpfoyrzQwMXafymd/eM5N2pXote/5Vm9A+ODx5BzZQfp3yz1qqleSJF5qPoIWZZowbddsTiWduaP3cdYz3aVpsFqtBAfn5aKqWLEi/v7+TJ06lRdeeIFy5crlB13uit1uZ/ny5Xh7e1OmjGckd7NarUycOJGZM2cWm4X9t0MIwXO7FX64LPium8y9pfQZTBve3RjLzu/I3PQ5wX2eJLBVT63l3JA79S7HYWXCjtdJz81kcbdZlAgonPxHs5qZSMwRPPmLQrhP3hotg/+i577n37AtquUF0r5ejCkojOAuD2otyWWYTTLT2owjeusUJux4gyXdZlM+uGDpVpwd9NFdgPV3vL29+eKLL5g+fTqTJk3SWo7L8Pb2Zv78+YVeRd1AOyYeVPnojGBVB5nexo3Frcg+tpu0b5cR2PEBgjp7zo0FwKEqTN89l4vpsbzTZUaBbywFQZIk3m8rk5yrMHCLwpbe0Lqk0RfchcA2fVAyU8n48WPkoFACWnTTWpLL8DX78FaHKYz8eQIvb3uN97rPJsLP9TnAJKGT8b+MjAxCQkJIT0/PH8G6jqqqJCcnk5qaqqv56jvBZDJRpkwZI7jyYBafVBi1V2V+SxNj67n3iGtxI/fCSRLfm4BfvdaED34FyUOm8SFvVHX+gaWsP7+Zme0n07JskyI5b45D0GODwvFUQUw/MzVD9TdVbnBjhBCkrV5E1v5NRAybhl/dFlpLcikJWYmM2PQKIb7BvNNlBoHeziUQvVW88g+ETkhPTxeAiI+P11qKQQHJzs4Ww4YNE9nZ2VpL0ZzvLirCtNwmxu51aC3FKQzv/ost/rL469WB4tqil4Vqz9Vazm0pqHefHl8j2n3WT/z4x8+FrOzfpFhVUXu1TVT+wiauZatFfn494i59T1UcIvGj6eLKy/1E7sXTWstxORdSL4leqx8Ro3+eKHIdNqdeEx8fLwCRnp5+y+N09/PMUxZ+FydMJhPlypUr9t4dTFR5eKvCgIoSc1u6x2dheJeHkplK0gdTMAWHETlsGpJZ/3UZC+LdpgvbWX7sU4bWe5je1boWgbp/EuYj8VMPM9kO6LtJIduhi4kTTXGXvieZZCIefwWvclVJWj4NR3K81pJcSuXQCsxqP5lTSWeYvW+RU4v6nfXMLaYIDQz0zp8Zgpbf59UX3NbbKHrrTqi5VhIXj0NJT6bE2IWYw0poLcmlHLp6lPHbX6d7lY6MbzHSpUVvC6wlUaX9jwrdy0ms6Swjm4x+4i4oljQSFoxFMntR4sX5mPw9a5nLtku7mb77bZ6o9xBP1X/0lsc6G6/oLnTOzs7WWoJBAcnOzmbQoEHF1rvUXEGvTQ6CvOD7bu4VXBV374SikLJqFo6EK0QOf8OtgitnvDufepEpO9+iSakGvNT8eU2DK4CmUSa+6iTz3SXBy/vdez3t3eJufU8ODCVy+OsomSkkr3jT4xKRdqrYhuENB7Py+FdsvLDtlsc665nuAix3T8NQHJFlmVatWhVL73IVwX2bFRJy4KfuZqJ0UgLHWYqzdwBp697H+vtBIoZOwrtcVa3lFIjbeZeSk8arO96kTFAppre98xI4rqZPRRPvtjax8ITKohPFNxGpO/Y9r5LliRw2jdwLJ0ld7dx0mjvxaJ0H6F21C2/vX8LRayduepyznhlThAYGd4gQgsE7FL7+U7C1l35zXRncGMvuH0j7egmhD44msHUvreW4lFzFxpgtk4m3JPB+z7mU8I/UWtK/GLdfYd5vKuu6yfSraPQddyL70DZSPptDcK8hBHe79XSau+FQHYzfPp2zKRd4r/tsKgSX+9cxbjtFmJXlWZW8iwNZWVl079692Hn31jGVz88JPmnvvsFVcfXOeuYIaWuXEtj+PrcNrm7mnRCC2THvci71T2Z2mKTL4ApgdnMT91eWeHSbwm/JuvidX6S4c9/zb9qJ4J6DyfhpFdmH766un94wm8xMb/sKEX5hTNj+BmnWjH8d46xnursreHl5aS3BoIB4eXkxaNCgYuXduosqEw+qTG1s4uGquutGTlMcvbNfiyV5xQx8azYhpN/TWsu5Y27m3aoTq9l6aScTW42hdkR1jdTdHpMksbK9TPUQ6Pezg4Sc4hVkuXvfC+r2KP7NupDyxXxsl+6s5IxeCfIO5K0OU8h25DB550xsyj/XmznrmTFFaGBQQI4lC+793kGP8hKrO8uYNF44bOA8SlYGiQvGgNlMiRcXYPJzLrGgu3B9J9Sw+o8xpJ57ZKGPtQiarXNQPVhiS28ZH9noT+6CcNhIXPwKjpRrlIxehByqz9HSO+Vk0hle3DyRbpU7MO5vO3CNKUKDIiMrK4tWrVoVC+8ScgT9fnZQIwRWtnf/4Ko4eSccdlJWvImak0Xk09PdPrj6X+9+TzrLrJh36FqpPYPvGaSxOucpHyixrqvMwSTB87sVj1s4fTM8oe9JZm8inpqCZDKR9NF0hC1Xa0kupW5kTV5u8QLrz29h7dn1+f/vtlOE3t76T/Bn8E+8vb2Jjo72eO9yFcH9mxVyFfium5kAL/cOrqD4eCeEIO2b98j98xQRT03BHFlaa0l3zd+9S8pOZtLOmVQPq8y4ltrmuroTWpY08WFbmRVnBQuOF4/0DZ7S9+TgcCKGTcMRf5mUrxZ6XIDco0onHqo9gCWHP+LQ1aOA83GKMUVoYOAEQgie2qnwxXnBL31kWpTQ3W8Tg1twfcdg2CPRHlW0FsCm2Hlx80QSs5P5oOd8wv1CtZZ0x7x6QGHObyo/dJPpZRRJdyuyf91JysqZBPd5kuAuD2ktx6UoqsKEHW/ye/JZlvWYS7AIcM8pQovForUEgwJisVioW7euR3v3zgmVT84KPmzrWcFVcfAu9/wJ0tYuI7Bdf48Krq57N2dv3o7BN9q/6tbBFcCMZiZ6l5d4eJvC6TRd/PYvNDyt7/k3akdQt0fIWP8JOSf2aS3HpcgmmaltXiLUJ5iJO94kIS3Rqdfp7k7h6+urtQSDAuLr68v8+fM91rvtcSov71d5ub6Jx6vrrsvcFZ7unSMtkeQVb+JduQ4h/Z/RWo5L8fX15bEZT7H58i+81GKErncMOotJkvi8o0y5ALhvs4MMm+cGWZ7Y94J7DMb3nlakfDobe/wlreW4lCDvQGZ2mERyTipLT6906jW6u1uYzfrINmzgPGazme7du3ukd7EWwUNbFTqUlpjVTHfd5a7xZO+E3Ubyx28gmb2IGDoJSfasNp5MPsOW7D08ULMPPap00lqOywjylvi2q5m4LBj6i4Kqj1UsLscT+55kMhH++DjksBIkf/Q6ao77LuC/ERWCyzHl3pc4kvibU8fr7o6RmZmptQSDApKZmUm5cuU8zjurI29Ru78ZvuwsY/bAwrSe6p0QgtQ1i7FfvUjEsKnIQe49dfa/JGQlMmXnW1jOpTG4+kCt5bicmqESn3aU+fai4K2jnrno3VP7nsnHj4hhU1EyU0n5fC5C9Sz/WpZtwoymrzp1rO4CLD8/P60lGBQQPz8/1qxZ41HeCSEYsUfhRKpgbVczkb6eF1yBZ3oHkLX7B7IP/EzYg6PxLu/+U2d/J9eRy+Sds/AxezOj0ySCAoK0llQo9KtoYmpjE5MPqWyM9aybNHhu3wPwiipL+OPjsZ6IIXPraq3luJxaUc5dU4xdhAYGN2DpKYURe1RWtpcZUkN3v0MMbkHu+eMkLplAYJu+hN7/nNZyXIoQgrf2LWL7pd0s7jabGuFVtJZUqKhC0G+Twp5rgkP3maka7Jk/dDyV9J9Wkbn5CyKffRPfWk20luMy3DbRaEbGv+v+GOibjIwMgoODPca7vddUXoxRGVnH5PHBlad5p6Qnk/zJzP8sanffMjg348fzm9l4YRvRzZ+nlDnSo7y7ESZJ4rOOMpG+eYves+y6GA9wCZ7W925EcI/H8KnZhJRVb+FIjtdajstw1jPd3T0CAtw7u3JxJCAggJiYGI/w7lq2YOAWhRZREvNb6a57uBxP8k4oCsmr3gJJImLoRI9b1H425TyLDn5Av2rd6VGlk0d5dytCffIWvV/IgOG7PCfTe3HwTzLJRAx+BckvgOQVb3hMpndnPdPdHUSWZa0lGBQQWZapW7eu23unqIJHtyuoAlZ3kfHywEXt/4uneAeQvv4TbH+eJOKJichBYVrLcSmZuRam7HyLyqEVGNk0b2TOk7y7HfeES3zYTub/zgve/90z1mMVF/9MAUFEPDkF+7VYUte86xEBsrOe6S7A8uThUk8lIyMDSZLc3rvXjqjsuCr4spNMaX/PD67Ac7zLOR6DZdsaQvo8hU/Ve7SW41JUoTIjZiFZ9mymt30FHzmvTIeneOcsD1c1MaKOiRdjVA4nuv9Nujj5512uKmEPjib74Bay92/SWs5d47ZThIGBgVpLMCgggYGBxMbGurV3G2JV3vxV5Y0mJjqU0V23KDQ8wTtHUhwpn8/Ft15rAjs+oLUcl/N/J9cS89dBJrUeS+nAkvn/7wneFZT5LU3UD5cYtNVBaq57B1nFzb+AZl0IaNWT1G/ewxZ3QWs5d4WznunuTuJuRUoN8jwLDg52W+8uWwSPb1foVV5iQkPddYlCxd29E7ZcklfMwBQYTPgj0W7bjptxJP43PvrtcwbXHUSrsk3/8Zy7e3cn+MgSqzvLpObCk7+493qs4uhf6H3P4RVVjpQVM1Ct2VrLuWOc9Ux3dxNPS7pWHMjMzCQkJMQtvbMpgge3KgR6waoOMqZidLED9/YOIG3tUuwJsUQMnYzJ37NGApJzUnh9z1walbyHJ+s/8q/n3d27O6VysMSqDjLfXRLMO+6+67GKo3+Stw/hQyeiZKSS+tU7bhsgO+uZ7gKsoCDPTJrnyQQFBZGenu6W3o0/oHIkSbC6s0yEhyYTvRXu7F3WwS1k7dtI2AMv4F2uqtZyXIqiKryxZz4mycSUe19CNv17Ua07e3e39K1oYnx9ExMOqOyOd88gq7j651WiHGEPjyHn11/I2vuT1nLuCGc9012A5a4RbXFGCEFGRobbeffNnyrvnFCZ18JEixK66wpFgrt6Z78WS9qaxfg360JAy+5ay3E5n55cw7GEk0y59yXCfG9c5sddvXMVM5qZuLekxENbFRJz3O8zKM7++TdqR0CbvqStXYYt9g+t5RQYZz3T3V3FYrFoLcGggFgsFsqXL+9W3l3MFAzbqTCwssTIurrrBkWGO3on7DZSVs5CDo0kdOALWstxOUevnWDl8a944p6HaFSy3k2Pc0fvXInZJPFFJxm7Ck+4YVHo4u5f6IBn8CpTieRPZrpdUWhnPdPdncUok+N+BAcHI4RwG+/squCRbQqh3rC8rVysFpn+L+7mHUDad8uxJ8QS/sRETD6eVcctzZrOG3vmUb9EHQbfM+iWx7qjd66mTEDeeqwNsYKFbrYeq7j7J5m9iRg6CTUrg9SvFrrVSJ6znukuwFIURWsJBgVEURROnjzpNt5NO6xyMDEv31WoT/ENrsD9vMs5tpus3T8QOmA43mU9qw6fKlRmxbyDXXUwpXX0Dddd/R13866w6FHexMv1TUw4qHIo0X2CLMM/MEeUyluPdXQX2fs2ai3HaZz1THcBVlaWew0VGuR51qpVK7fwbstfKm8dVXmzqYmWJXX39S9y3Mk7R3I8KV8uxK9BGwLu7aO1HJez5vT37Is7zKTWY4n0j7jt8e7kXWEzo6mJhhESD29TyLC5x0iI4V8e/g3bEtCqJ2lrl2GPv6y1HKdw1jPd3WGK63CpO3O9YKnevbuWLRi8XaFzWYnxDXT31dcEd/FOKA5SPp2NydefsIfGeNy07qmkM7z/6yoeqXM/Lco0duo17uJdUeAt563HSsiB53e7R34sw7//EnLfs8jhJUlZNQtht2kt57a47RShw+HQWoJBAXE4HMTExOjaO1WIvIWwwKfFMN/VzXAH7wAyflqF7fJZwp941ePyXWXZs3l9zzxqhlfl6QaPOf06d/GuqKgaLPF+m7x6hSv/0H+AZfj3X0zevoQ/MQF7whXSvv9Qazm3xVnP8gOsP/74g9atW1OjRg2aNWvGyZMn/3XwihUraNiwYf5fZGQk999/PwAXL15EluV/PH/+/PkCC8/JySnwawy0JScnh0GDBunau/nHVTZdEaxqL1OqmNQZdAZ38M569iiZ29YQ3OsJfCrV1lqOy1l48H3SrRlMbfMyZpPZ6de5g3dFzSPVTDxZQ+KFPQqn0/QdZBn+/RPvMlUI7f8MWbu+J+dEjNZybomznkniP2OpnTp1YsiQIQwdOpSvv/6a2bNnc/DgwVu++J577mH69Ok88MADXLx4kYYNG5KWlnZHgjMyMggJCSE9Pd0YMjVwKYcSVVp9pzC2nok5LTy7cr2noWZlcm3O85hLlCPy+ZlIJt0Nut8VWy7u5I0985jUeizdKnfQWo5HkGUXNPnWga8M+weY8ZGNH1TughCC5I+mY7twkpLjlyKHRmot6YY4G6+YABISEjh06BCPP/44AA888ACxsbGcO3fupi/cv38/CQkJ9OvXz6XCjeFS98PhcLBp0yZdepdlFzy2XaF+BLzZ1LNuzq5Az94JIUhd/Q6q3UrYYy95XHAVb0lgwYFldK7Yjq6V2hf49Xr2TksCvCS+7GzmVBpMPqTfXYWGf/9GkiTCHolG8vIh+dPZCFWfOywLNEUYGxtL6dKlMZvzhqclSaJChQpcvnzzFf0fffQRgwcPxsvLK///srKyaNasGY0bN+b111+/5VbG3NxcMjIy/vEH5I+AWa1WrFYrkDccl5ubC0B2dnb+46ysLGw2W/5ju90O5CUBu/4BZGZm5j/OyMjI15SRkYGqqv/Ipquqar4ORVHyHzscjvzaQw6HIz/JmN1uz99NYLPZ8h/n5uaSnZ2d//j6cKKntiktLY3o6GhSU1N116aX9qnEWmB5CxveslSsfbpRm1JSUhgzZgxWq1V3bUrd9SM5x3YTdP8L2LwDPMqnbGsOb+yeR6B3AC80eDK/HQVpk9VqZcyYMfnn1bpNeupP9cMEU+rZmfebyuZYhy7bZLVaefHFF/M1FEefbtQmOSCY4EeisV04Qea2b3TZJqd3fgohxKFDh0SNGjXE32nWrJnYunWruBEWi0UEBQWJkydP5v+f1WoV165dE0IIkZycLLp06SJmz559w9cLIcS0adME8K+/wYMHCyGEGDt2rBg7dqwQQohhw4aJadOmCSGEGDhwoJg3b54QQohu3bqJ5cuXCyGEaNmypVi9erUQQog6deqIjRs3CiGEKFu2rNi7d68QQoigoCBx4sQJ8Z9pUREbGyvS09MFINLT00VsbKz4z0ciTpw4IYKCgoQQQuzdu1eULVtWCCHExo0bRZ06dYQQQqxevVq0bNlSCCHE8uXLRbdu3YQQQsybN08MHDgwv53Dhg0z2qRBm8r3fk7wgU08/ekBj2mTJ/p0ozZ9vniB+GNUN5H8f/M8pk1/92n4stGi3af9xNFrJzymTbrzSZJExx/souQnWSKwRDnPaJMn+nSTNs3s0VjERvcWO75apbs27du3L78NtwIhhLh27ZoICgoSdrtdCCGEqqqiZMmS4o8//rjhi1asWJHfyJvxf//3f6JPnz43fd5qtYr09PT8v+sfdGxsrBBCiJycHJGTkyOEECI7O1tYrVYhhBBZWVn5jy0Wi8jNzc1/bLPZhBBCZGZm5rclIyMj/3F6erpwOBz5jxVFEaqqivT0dKGqqlAUJf8Dczgc+Y/tdrvIyMjIf5yZmY0VhcYAACAASURBVCmEEMJmswmLxSKEECI3Nzf/sdVqFVlZWfmPs7OzPbpNGRkZYvXq1SI9PV03bbqSYRORK3NF3412kZtr+HSzNqWlpYnPP/88//300CbVYRdX544Sf01/QijWbI/z6WTiadHh8wFi6cEVd9Umm80mVq1alf/ZuNt3ryh8upypitBPbOK+DdlCVVVdtclms4mVK1fmv2dx9ulGbUpPSRbxc0aIuBlPi4yUJF216dq1a84HWEII0b59e7FixQohhBBr1qwRTZo0uemL2rRpkx/xXefatWv5wqxWqxg4cKCYMmXKLU/+d65HtXFxcU6/xkAfWCwW0bJly/wvu9aoqip6bbCLEp/axLVsVWs5ukZv3gkhRNr6T0Ts2J7C+ufvWktxOVm2LPHIuuHi2Q0vCbtiv6v30qN3euSrc4rgA5tYeUbRWso/MPy7Pba4P0XsS31E6tdLtJbyD+Li4pwKsPJ3EZ45c4ahQ4eSnJxMcHAwK1asoF69ejz99NP069cvfzH7mTNnaNq0KXFxcQQFBeVPNa5du5apU6ciyzIOh4NOnToxd+5cfHx8nJqqNHYRGriK904pvLBHZX13mV4VPGthtKeTe+EEie+OJ7jn4wR3e1RrOS5nzr7FbLu0iw97LaRcUGmt5RQbhmx3sO6S4Nj9ZioHG7sK3YnMX9aR/u0yIp99E9/aTbWWAzgfr+QHWFpzXXBiYiKRkfrcmmlwY2w2G6tWrWLIkCF4e3trquX3VEHjbx08VdPEknuNlAy3Q0/eqdZsrr09Ajk4nKhRbyPdphafu7H3ykFe/eVNXm7xAn2rdbvr99OTd3onwyZo8I2DsgESO/rImE3aB1mGf84hVJWk9ydjv3qRkq8sQw7QfgAmKSmJqKgo59I06InrK/gN3Ae73c6aNWs0986mCB7b7qBSELzdQndfbV2iF+8A0r9bjpqZRvhjL3tccJVmzeDt/YtpVaYpfap2dcl76sk7vRPsLfFZR5mYBMHsY/pI3WD45xySyUT4oy+Bw07a6kW6KIPkrGe6G8EypggN7pRphxVm/qqyr7+ZJlHa/0I1cJ6ck/tJXj6N0AdHE9i6l9ZyXIoQgmm7ZnM04QQrer9LhF+Y1pKKLRMPKsz9TeXAADMNI4xrhDuRfXQXKZ/MIOyRaAJa3P0I8N1QoESjeuJ6bgoD9yE3N5f58+dr6t3BRJUZv6pMbmQygqsCoAfvFEs6qV8uxLd2MwJa9dRMR2Gx5eIv/BIbQ3Tz510aXOnBO3djWmMTtUNhyA4HuYq2YwuGfwXDv2Fb/Jt1IW3tMhwp1zTV4qxnuguwbpWc1ECfKIpCTEyMZt7lOARP7FBoGCExsZHuvtK6RmvvhBCkfb0YodgJe3gMkocV4U7ITmLhwQ/oWqk9HSrc69L31to7d8RHlljVwczpNJh+RNupQsO/ghN6//OY/PxJ/XIBQtXOP2c9M6YIDdyel/YpLDmlcniAmbrhnnWD9nSyD28n5dPZhA95Ff/GBS8Xo2dUofLytte4/P/snXd4FNX6x78zs6mb3YRmQJq0AAklAYIEVCxIRIooBLuCIKjXK8hPr14LVdF7Va7tqhRFsAHBQhUUEaUECRB6b5JQQoBke515f3/E5KJAMpvs7pQ9n+fJ8yxky3fy2ZN9854zZyxFmNP/XZhiEpSOxPiDaQUiXt4qYf1AAVnJ7I8yLeHevxXnPnoRSUP/hoTrBiqSgU0RMsKGx+PBpEmTFHH3y2kJ/9kl4ZVuPCuuaoCS7sSycyhd9F/EZfTWXXEFAIsPfo+tZ3bguaynQlJcKelO6/yjM4/MBhweXivC4VOmx8D81YzYdl1h7DUAliWz4T93SpEMmp0ilBRs+zFqhiRJKCoqCrs7m5cw4hcRvZI5PN1BdW9lTaCUOyLChfn/ARcVjTpDnwzra4eDIttpfFjwKe5MuR2ZjdJD8hpKudMDBp7D3N4CCh3AP/OV+fkxfzUncdBI8Ka6uPDlW4pcEFquMzZFyNAsY9aJ+OKwhB1DDGjFNg/UFPa871G24B3UGz0VcamZSscJKhJJeOrHF3HBdQEf938HcYZYpSMxrsA7u0WMy5Ow+nYBtzRmf6RpCc+RXSh5/x9IHDQKppuGhPW1NTtFWHHla4Z2cLvdGD9+fFjdrSyUMHO/hDev5VlxVQuUcOcvPQvLd7MQ3yNbd8UVAHx9YBl2lezFcz2eCmlxpYQ7vfH3NB43NuIw4hcRVm94ew3MX+2IadURCTcMhmX5p/CdORHW15brTHUFFoNRHRYv4dF1Ivo05jCmPXsLawkiQun8t8HHxSPpjtFKxwk6RdZTmLX9M9zVtj86J6cpHYdRDTzHYU5vAaVe4Nnf2FSd1kjsPxyGusm48OWbIBWejcmmCBmaY/Q6P746Qtg9xIDmJta90hKOTStROv9tNjXIUBUV1y9lU4Xaw/v7AZx9+2mYb38I5lvvCctranaK0OVyKR2BESAulwujRo0Ki7vVJyXM2k/4d3eeFVdBIJzu/KUlKPtuJuK799VdcQWEb2qwgnC60zuPtS+fKhz1qwh7mM4qZP6CQ3TztjDdPBTWlV+EbapQrjPVFVg8r7pIjGrgeR5NmjQJuTublzDqVxE3NmJTg8EiXO6ICKUL3gEXE4ekwWxqMBiEy10kwHMcZt8g4KwbeH5zeKYKmb/gYc6+H4Z6ySj9anpYziqU64xNETI0w982iPj0oIRdQwxoyRa2awrHbz+g9KvpqPfoZMSlXat0nKAiSiLGrn4RF1ylbGpQ41ScVbh2gIDejVjhoyU8R/eg5L1nkHjHozDdeFdIX0uzU4ROp1PpCIwAcTqdyMnJCam7tackfLBXwuuZPCuugkg43PnLSlD27QzEZ/bRXXEFAN8eXIFdJfvCNjVYQTjcRRp/T+NxXTKHkb+KcPpD23tg/oJLTMs0JFx/B6zL58JfEtoNSOU6U12BJQiC0hEYASIIArKyskLmzuEjjPxVxPUNOfwtTXVvWU0TandEhLKF74GLjkHSnWNC8hpKcsp2BrO2f4Y7U8J/1mCo3UUiPMfhk94CTjqAF0O8ASnzF3zM/YeDN9fBhRBfq1CuMzZFyFA94/JEzNgnYecQA9oksu6Vlqi41mC9URMR1yFL6ThBhYjwf2smoMh6Gp8OeA/xUXFKR2IEibd2inj2NwnrBgro1ZD9Uacl3Id24Nx/n0PS0CeRcN2AkLyGZqcIHQ6H0hEYAeJwOJCdnR0Sd3nFEt7dLeHVbjwrrkJAKN2JdgvKvvkIcek36K64AoAVR3/C1jM78cy1TyhSXIXSXaQzrgOPa68qnyp0h2iqkPkLDbFtOsPYsz8sSz+G/0JxSF5DrjPVFVhRUVFKR2AESFRUFHJycoLuziuWnzXYrQGHsexagyEhVO4AwPLdDBCJSBryeNCfW2nOOc/jg62f4LaWN6P71V0UyRBKd5GOwHOYfb2AozZg2vbQTDUxf6EjcdAj4OMSULrgHYRikk6uMzZFyFAtU7aJmLpNwtY7DehUj3WvtIR73xacm/ES6tz3fzB2v1XpOEGFiPDSr69hz7n9mDfgvzDHmJSOxAgRk7aKeLVAwra7DOhYl/0O0hL/+x30DIzd+wT1udkUISNsOBwOZGVlBdXd3lLCKwUSnuvMs+IqhITCneR2onThu4hp2wXxmcH9xaYGfinciPVFv2FstzGKFlehcMf4M/9M59EmEXh0nQhRCm4vgvkLLbHtuyGu602wLJ4J0V4W1OfW7BRhdHS00hEYARIdHY3x48cHzZ1E5VODLUzASxmqe4vqimC7AwDr8k8hOSyoM+wpcJy+imOrx4a382fguibX4sZmPRXNEgp3jD8TI5RvQLr5LOH9vcGdKmT+Qk/S4DEAESzfzgjq88p1prpPLzYfrT2CvZbgw70S8s4SZl0vINagrw9otRFsd55je2FfvxTm2x+GoV7DoDynmvjvtk/gE/14OnOM4sUjW8MTHnom83gilceL+RKO24LXxWL+Qo9gSkLi4DFwbv0Zrr35QXteuc5UV2DZ7XalIzACxG63Iy0tLSjuCu2E5/MljGnH4wa2k3LICaY78ntRuuBtRDVNQcINdwQhnbrIP70dK4+uweNdhqN+fD2l4wTVHaNqpmXyqBMDPLZeDNqiaeYvPMRn3oKYtl1QlvseJE9wrvso15nqPsFiY9llJrRGbGwspk+fXmt3RIQnNogwRQH/ulZ1b01dEix3AGBbnQv/2ZOoe884cLy+Nk90+z2YvvlDpCd3QP9W6li0H0x3jKoxR3P4sJeAVUWEL48Ep8Bi/sIDx3Gok/N3SHYLrCvmBeU55TpT3aeYwWBQOgIjQAwGA7Kzs2vtbuFRwrIThA96CUiMZlOD4SBY7nxni2D9cT5Mt+Qg6uoWQUqnHubumo9zzvN4pvsTik8NVhAsdwx5DGjO4+6WHMZuFFHiqn2RxfyFD0P9RjD3ewD2XxfDe+JA7Z9PpjPVFVg2m03pCIwAsdlsaNKkSa3clXoIY/NE3HUNh8HXqO5tqVuC4Y6IUJb7HoSk+jDfem8Q06mDI6XHsWDfd3iwwzA0NTdWOk4lwXDHCIx3sgRIAJ79Taz1czF/4SWh912IatwCpfPfAYn+Wj2XXGeq+ySLi2OXm9AacXFxyM3NrZW7f26W4PQD7/bU19SS2gmGO2f+angO7UCdYX8HFx0TxHTKI0oi3vztv2hqbox7U+9UOs6fCIY7RmAkx3P4d3cBcw8R1p6q3VmFzF944QQBde4eB9/p47D9/E2tnkuuM9UVWKxdqj0MBgOysrJq7C6vWMKM/eWXw2lsVMf0S6RQW3ei3QLL4lmI73oTYtsqs6N5KFl8aCX2nj+IZ7r/DVGCus72qq07Rs14pC2HXskcxqwX4RFrPlXI/IWf6KZtkNB7MGyrvoD//JkaP49mpwitVqvSERgBYrVaYTaba+TOJxFGrxPRrT6HJ1JV93bUPbVxBwCWxbNAJCFx8JggJ1Oes85zmLX9MwxqnY2OV7VXOs4l1NYdo2bwHIcZ1wk4agX+taPmXSzmTxnM/R4EbzSh7Ov/1viMULnOVPeJZjQalY7ACBCj0Yi8vLwaufvPLgl7y4AZ1wsQeNa9Cje1cec+tAPO/NVIGjgSgikpBOmU5Z38mYg1xGJ0xkNKR7kstXHHqB1pdTk825nHqwUSDpbV7EOa+VMGPiYOSXc9AffefLh2bqjRc8h1proCSxDYGhytIQgC0tLSAnZ3zEqYtFXC2DQeXeqz4koJauqOfF6ULXwX0S3TEH9tdojSKce6wk1YX/Qbnur2KEzRCUrHuSw1dccIDi9l8GhsBB7fULO9sZg/5Yjr1BOxHXqg7JsPIbmdAT9erjPVFVisXao9rFYrOI4LyB0R4cmNIurFAlO6qe5tGDHUxB0AWFcvgP9CcfnlcHh9+XP6nHgnfyayru6m+OVwqqKm7hjBId7A4YNeAtacInxxOPACi/lTlqS7Hge57DXaG0uzU4QJCer8a5FxZRISElBYWBiQu6+PEVYUEt7vKSAhinWvlKIm7nzFhbCtXgjTzUMR1bB5CNMpw5yd82H12jA2c7Rq9ry6HDVxxwgutzUt3xtr/CYRF9yBFVnMn7IY6ibDfNuDsK9bAm/hoYAeK9eZ6gosNf9CY1wejuNgNptlu7N4CU/liRjcnMMdbM8rRQnUHRGh7OsPICTV0+WeV4dLj+HrA0vxcMe70SghWek4VRKoO0Zo+E+WAK8EPLc5sL2xmD/lSeg9GFENm6F04bsgSb4/uc5U9+nGNl3THjabDYmJibLdTdgiweple16pgUDduQp+hedgAZKGPKG7Pa8kkjB984doam6MYe3Ufy3FQN0xQkOjeA6vZfKYfYCwsVj+WYXMn/JwggFJw56Cr/AQHBuWy36cZjcaNZlMSkdgBIjJZILFYpHlbvt5wvt7JUzswqNpAvvLTWkCcSe5HSj7bgbiOvVCXGr3MKQLL8uPrMaecwcwPvNx1e15dTkCcccILaPb8eVbzawX4ZfkTRUyf+ogpkUqjFn9YFn2KUTLeVmPketMdQVWsK5UzggfRASr1VqtO4kIT6wX0S4RGNdRdW+9iESuOwCwrvgM5HYi8c7HwpAsvJS5LZhRMBe3tbwZnZPTlI4ji0DcMUKLwHP48DoeOy8AH+yV18Vi/tRD4sBHwEVFoWzJbFn3l+tMdZ9ydrtd6QiMALHb7WjatGm17j49SMg7S/jgOgFRbM8rVSDXnbfoCOzrlsB82/0w1GkQpnTh48OCTwEAj2cMVzRHIMh1xwgP3RrwGNOex0tbJJx2Vv8BzPypBz7ehMSBI+Ha+jPch3ZUe3+5zlRXYJnNZqUjMALEbDaDiKp0d95N+MdvIh5ozaF3I9W97SIWOe5IklCW+x4MyU2R0Ftd1+MLBjuK92Dl0TUYnf4QkmITlY4jGznuGOHl1W48YgTgmU3VL5hm/tRFfGYfRLdIRdmi/1Z7MWi5zlT3SSeKtb9KOSO8iKKIPXv2VOnuhXwJPgl441q2sF1NyHHn/G0VvL/vR52cv4MT9HXdNJ/ow/T8D5Favy0GtL5V6TgBIccdI7zUjeXwxrUCvjxC+Lmai0Ezf+qC43kkDX0S/rNFsP/ybZX3letMdQWWw+FQOgIjQBwOB7Kysq7obvNZCbP2S3ilG4+G8WxqUE1U5060W2BZ+gniM/sgplWHMKcLPbn7l+KE9STGZz4GnlPdr8Mqqc4dQxkealN+MegnNojwVnExaOZPfUQ3bomE6wfCuvIL+MtKrng/uc5U9xuFtUu1R8UFSy/nTpQIj28Q0bke8Di7mLPqqModAFiWfQIiQuKgUWFOFnpKnOcxb/cC3JlyO9rUbal0nICpzh1DGXiufIf3Q5bya61eCeZPnZj7PQQuJg6W72Ze+T5anSL0+6ue+2SoD7/fj7y8vMu6m7FfwrZzwIe9BBjYwnbVUZU7z/H9cG5ahcT+D+vyYs4fbJuDWEMMRnTS5oapVbljKEunehyeSuMxpUDCCfvlu1jMnzrh44xIvGMUXNvXwX1g22XvI9eZ6gosl8uldARGgLhcLuTk5FzirsRFeDFfwqi2HHokq+6txsCV3ZEkouzr/yKqcSsYe96uULrQUVC8C2t+X4fHMoar9mLO1XEldwx1MKkrj8RoYPwVFrwzf+olvuvNiG7VsXzBu997yfflOlPdpx7bdE17mEwmFBUVXeLuhXwRHAe81p0tbFcrV3Ln+O0H+AoPIWno38Dx+vLnl/x4J38m0uq3Rd8WNyodp8ZcyR1DHZijObx5rYCvjxFWn7x0qpD5Uy8cx6HO0L/Bf/40bD9/c8n3NbvRKGuXag+/349Vq1b9yV1+iYSPDxCmduVRP5ZNDaqVy7mTHDZYl80pX9jeIlXBdKHh24Mr8Lu1COM0uLD9Yi7njqEu7m3F4fqGHJ7aKML3lx3emT91E9XoGiTcMBi2H7+Cv/TPC941O0XodruVjsAIELfbjfHjx1e6k4jw5AYJHesCY9qr7i3GuIi/ugMAy4q5INGPxEEjFUwWGs67SjFn51cY1DobKRpc2H4xl3PHUBccx+G9ngIOWID3dv+5i8X8qR/zbfeDi4mH5S87vMt1prpPv4QEba6HiGQSEhKwZ8+eSnefHiRsLiG835MtbFc7f3XnLTwEx8blMPd7CIKpjsLpgs+Mgrkw8AJGdr5f6Si15q/uGOqkcz0Oj7fnMWnbn3d4Z/7UDx9rROLAR+Aq+OVPO7zLdaa6Asvn8ykdgREgPp8Pubm58Pl8KPMQnt8s4r5WHK5nO7arnovdkSSh7OsPYGjYHAnXDVQ6WtDZdXYfVh37GaPTH4I5RvvrXi52x1A3U7ryiOaB5zf/b8E786cN4rvdgujm7WD55kPQHxuMynWmuk9Ar/fSFfsMdeP1ejF9+nR4vV5M3CrBJbId27XCxe6c+avhPb4PdYY8AU7Qlz9REvH2lhloV68Nbm/VR+k4QeFidwx1UzeWw2vdBcw7RNhYXD5VyPxpg/Id3v8G35nf4diwDID8OkV1BZbRaFQ6AiNAjEYj8vLycMwTj//ulfByBo+rjWxqUAtUuIvjAcvSTxDX5UbEtO6kdKygs/zIjzhcegzjuo3W9ML2i6lwx35naoNHUjh0rc/hyQ0iRImYPw0R3bQNjD1ug+X7zyDay2Q7U91vGlbNaw+v14tZs2bjyQ0+tDID4zqo7m3FuAJerxezZ89G6Yp5IK8bSTrcsd3msWP2js9xW8ub0b5+itJxgkaFO/Y7UxsIPIf3e/IoOA/MPiAxfxrD3P9hAIB1+VztdrDYfLT28Pl8eH/TGfxyhsO7WQKiBda90go+nw/rF+fCtWEZTLfeAyGpvtKRgs7c3fPhFX0Ynf6g0lGCClvDoz16JPMYnsLhhXwJZ+zMn5YQEpKQ2O9BODathOv3A7Ieo7oCi7VLtQcXE4/SW57DoOYcspuq7i3FqIL4+Hi80b8bhDpXwXTjXUrHCTrHLYX45sAKPNRhGOrF1VU6TlAxGo1YtWoV+52pMV7PFOCTgDf2xTJ/GsPYawCiGjaH94fPZd1fdZ+GHo9H6QiMAHltmxen7CJey7j8JSEY6sW2fT08+7fC2H8EuKhopeMEFSLC+1tno6GxAYa2G6R0nKDj8Xgwffp09jtTYyTHc5jQhceH+0Q8+/anzJ+G4AQBSUOegPfEIVn3V12BJYrsQ1pLnLAT3trNofXx5bjGyNxpCfJ7YV/6MQ64BQhtuyodJ+jkndyC/NPb8USXRxAtRCkdJ+iIooi8vDz2O1ODPJXGo0UCMM/dGX4/86clYlp3gjHn77Luq7oCKz4+XukIjAB4brOIxGgO+ZPuYu40hm3tt5DKzuKGSR/obprCJ/rw/raP0bVhZ/Rq0l3pOCEhPj4eubm5bNxpkGiBw3+yDDhbrwNWl8QqHYcRIEndbpJ1P9UVWKxdqh3Wn5Ew/whhSoaEt6ZNZu40hGg5D9sPXyEuqz9e/WiO7tx9fWAZztiL8WTXkeA4fZ504fF4MGnSJN25ixT6JPvQyn0E4zf54RGp+gcwVIPcMae6AkuSLr3qOEN9SEQYmyeiW30O97WQUFRUxNxpCMvST8BFxSD2lrt15+6Cqwxzdy3AoDb90DKpudJxQoYksXGnZYgkpB/4HL/bgXd2M4daQu6Y44hIFaWz1WpFYmIiLBYLzGaz0nEY1fDJAQkjfxWxcZCArGTV1emMKvAc34eSt59G0rCnkNDzdqXjBJ1/b3of6wo34YtBH+rikjgMfTMuT8THByQcGmZAw3h9dlv1htx6RXWfjOzK4urH6iX8M1/E/a05ZCXz7KrwGoIkCWXffIioJq1g7JGtO3eHS49hxZHVGN7pHt0XV3pzF2lU+Hsu1YcYHnghny121wpyx5zqCiyG+nmlQILdV76fC0NbOLeuge/EQSTd+Rg4Xl/+iAgfbPsETUyNcEeb25SOw2DIok4M8Eo3HnMOEvJL2FShnmBThIyAOGwhpC7yY0IGj5e66OsDWu9IHjeKp41C9DXtUG/ES0rHCTp5J7fg+bVTMa33i7o9c5ChT0SJ0OVbPxKiOKwfKOj2xAy9oNkpQpfLpXQERhX8Y7OIhnHA/3X631vH5XJh1KhRzJ3Ksa3JhWi3IHHgyMr/04s7v+THB9vmICO5I3o2zlQ6TljQi7tI5WJ/As/hPz0EbCwm5B5VRc+DUQVyx5zqCiyeV10kxh+sPSXh2+OE17sLiDP87y8snufRpEkT5k7F+MtKYF+zCAm9B8NQv1Hl/+vF3dLDP6DQehJPdHkkYv7614u7SOWv/m5uzGNgMw7PbRbh9rMiS83IHXNsipAhC1EiZH7nR7TAIW8Qa2FrjQufvwH3/q1o+NLH4GP1tamo3evAfUseQ1bjbvhn1lil4zAYNeZgGSFtkR9Tu/F4Pp0twVArmp0idDqdSkdgXIZ5hwgF54H/9OAvKa6cTidycnKYO5Xi/f0AnFt+gvn2hy4prvTg7vM9i+DxezCq8wNKRwkrenAXyVzOX0oSh7+l8Zi2XUKxUxW9D8ZlkDvmVFdgCQKr2tWG3Ud4cYuIu1tyl93zShAEZGVlMXcqhIhQ9t0MRF3dAsYe2Zd8X+vuTtuLsWj/EtybehcaxNdTOk5Y0bq7SOdK/iZk8DDwwISt7IxCtSJ3zLEpQka1TNgi4t87JRzIMaC5iU0Nagnntl9wYd5rqP/Ea4hNyVA6TtCZtO7f2FWyD58P+hBxBnZNN4Y+eGe3iPGbJGy/y4COddnvXLWh2SlCh8OhdATGRRTaCW/ulDC+I3/F4srhcCA7O5u5Uxnk9cCy9GPEduhxxeJKy+52l+zHzyc2YFTnByKyuNKyO0bV/p5I5dHaDIzfJEIlPRDGRcgdc6orsKKiopSOwLiIF/JFmKKB5ztf+a0SFRWFnJwc5k5l2H75FqLlPBIHjbrifbTqjojwUcGnaF2nBfq2uFHpOIqgVXeMcqryF8VzeONaAatPElYUsgJLbcgdc2yKkHFFNp+VcO1iETOvF/BoO9XV4owqEO1lODP1ERiv7Yukux5TOk7QWVe4CS/9+hrevHkyMhulKx2HwQg6RIQ+K0ScdBB2DTUgimdThWqBTREyagURYfwmCR3rAo+kVD2wHQ4HsrKymDsVYV35BcBxMGXfV+X9tOjOL/kxo2AuMhulR3RxpUV3jP9RnT+O4zC9h4CDFuCjvWzBu5rQ7BRhdHS00hEYAL49TthQTHjrWgFCNX85RUdHY/z48cydSvAVF8KxcTnMfe+BYKy6G6xFd8sO/4gi22mMSX9Y6SiKokV3jP8hx1/nehxGpHCYvE2CxauKySYG5NcpbIqQcQlesXyzu9ZmDt/3MygdhxEg52ZPhu/kETR8YTa4KH19+Dp9Tty35DF0b9QFL/Qcp3QcBiPknHQQUhb68VQaj9e6sy051IBmpwjtRHXADAAAIABJREFUdrvSESKej/ZJOGoD3rhW3mC22+1IS0tj7lSA58guuHfnIXHACFnFldbczd/3HRxeJ0Z1vl/pKIqjNXeMPyPXX2Mjh2c68fjPbgkn7Kroh0Q8csec6gqs2NjIO91aTZR5CFO2SXgkhUMHmfuvxMbGYvr06cydwpAkoWzxLEQ1S0FcRm9Zj9GSu3PO81iw9zsMbTcQVxkbKB1HcbTkjnEpgfh7thOPpGjgxXwxDMkY1SF3zKmuwDIY2JSUkry2XYJLBKZ0k9+KNhgMyM7OZu4UxlXwC3wnDiLpjkfBybwYqZbczdk1HzGGaNyXNkTpKKpAS+4YlxKIv4QoDlO7Cfj8MGFrCetiKY3cMae6AstmsykdIWI5biO8s0fCs514NIqXf0qwzWZDkyZNmDsFIZ8XlmVzENshCzGtOsp+nFbcHbecwIojq/FQh7thik5QOo4q0Io7xuUJ1N+IFA6pScAzv7HNR5VGrjPVFVhxcXFKR4hYXswXUScaeKZTYG+LuLg45ObmMncKYl+3GKLlHBIHPhLQ47TibkbBPDQ0XoU72tymdBTVoBV3jMsTqD/DH5uPrj1NWH6CFVhKIteZ6gos1u5WhvwSCV8eIUztJiAhKrAN7QwGA7Kyspg7hRAdVlh/mA9jz/6ISm4a0GO14G5H8R5sPJmPUZ3vR5TAdi2vQAvuGFemJv76NeVwy9Ucnt0swi+xIkspNDtFaLValY4QcRARntkkoUOd8jZ0oFitVpjNZuZOIWyrFwAkwZwd+Jl1andHRPho+1y0rdsKNzW/Tuk4qkLt7hhVUxN/HMfhzR4CDpQBs/ezzUeVQq4z1RVYRqNR6QgRx5LfCb+eIbwhY1PRy2E0GpGXl8fcKYD/QjHsvy6B6eahEExJAT9e7e7WFW7C3nMH8FjGcPCc6n5dKYra3TGqpqb+0utxeKgNh4nbJNjY5qOKINeZ6n5jCQLbSC2c+CXC8/kibrmaQ3aTml3rShAEpKWlMXcKYP3+M/DxCUi48a4aPV7N7vySiJnb56F7owx0adhJ6TiqQ83uGNVTG39TuwmweoE3d7EulhLIdVZZYB06dAg9e/ZESkoKMjMzsWfPnkvuvHbtWsTFxSE9Pb3yy+VyVX7/448/Rps2bdCqVSs8+uij8Pl8AQdn7e7wMucgYX8Z8K/uAjiuZgWW1WoFx3HMXZjxnjoK55afYM6+H3xMzRY6q9ndiiOrUWQ7jdHpDykdRZWo2R2jemrjr2kCh6fSeLy1U8IZJ+tihZuApwjHjBmD0aNH4+DBg3juuecwfPjwyz6gbdu22L59e+VXxWr6Y8eO4eWXX8a6detw+PBhFBcXY+bMmQEHT0hgp2CHC4ePMHGriHtbcejaoOZXak9ISEBhYSFzF2asS+fAUP9qGLNqfmadWt25/G7M2fkV+lzTG23qtlQ6jipRqzuGPGrr7/l0HtECMHkb62KFG7nOeAA4e/YstmzZggceeAAAMGTIEBQWFuLw4cOyX3DRokUYNGgQGjZsCI7j8Nhjj+Grr74KOHhNuyiMwHl7t4RzbuDVADYVvRwcx8FsNjN3YcR9aAfc+/Jh7j8cnFDzs8jU6m7R/iWweW0Y2fk+paOoFrW6Y8ijtv7qxHB4MZ3HrP0SDpSxLlY4keuMB4DCwkI0atSo8tRDjuPQrFkznDhx4pIHHDlyBF26dEFmZiY++OCDyv8/ceIEmjdvXvnva6655rKPr8Dj8cBqtf7pCwBKSkoAAG63G263GwDgcrng8XgAAE6ns/K2w+GA1+utvF0xJWm32+H3+wGUbwhWcdtqtUIUxcrbkiSBiGC1WkFEkCSpMocoipW3/X5/5cZifr+/8jpEPp8PDocDAOD1eitvezweOJ3OytsV06hqOqbTdh/+tV3EE6k8msaLtTqmkpISJCYm4uzZs8xTGI6JiFC2ZDaEJm0Q1/m6Wh1TcXExEhMTYbPZVOOpzG3Fl3u+Qf8Wt6JRQrJmPV18TEDw33s2mw2JiYm4cOGCbo5Jj56udEwV/kpLS2t8TH9L5dE4Hng2z6OKY9Kjp8sdU1lZGWRBRLRlyxZKSUmhi8nMzKSffvrpT/9nsViorKyMiIgKCwupY8eOtGDBAiIievLJJ2natGmV992zZw81bdqUrsTEiRMJwCVfDzzwABERPf300/T0008TEdHIkSNp4sSJREQ0dOhQeuutt4iIqG/fvjRr1iwiIurRowctXLiQiIhSU1Np5cqVRETUuHFj2rhxIxERmUwm2r17N1H5NrhUWFhIFouFAJDFYqHCwkL640dCu3fvJpPJREREGzdupMaNGxMR0cqVKyk1NZWIiBYuXEg9evQgIqJZs2ZR3759iYjorbfeoqFDh1Ye58iRI1V3TEMXFRL37jkqcUm1PqZx48aRxWKhRx55hHkKwzE5Cn6hwrHZ9MitvYJyTK+++ipJkqQaT+9tmUW3fHYX3TpQ254uPqZgePrrMUmSRJmZmZW/g/VwTHr0dKVjkiSJrr76atqwYUOtjunl5fsJM7204Yyo+DHp0dPljikvL6/yGKoCRETFxcVkMpnI5/MREZEkSZScnEyHDh2q8sHTpk2jJ598koiI/v3vf9OYMWMqv7d8+XLq1avXFR/rdrvJYrFUflX8oE+dOkVERC6Xi1wuFxEROZ1OcrvdRETkcDgqb9vtdvJ4PJW3vV4vERHZbLbKY7FarZW3LRYL+f3+ytuiKJIkSWSxWEiSJBJFsfIH5vf7K2/7fD6yWq2Vt202GxEReb1estvtRETk8Xgqb7vdbnI4HJW3nU6nqo7psEWiqNlempDnCMoxORwOKiwsJLvdzjyF+JisZaV0euoIKv7wxaAck81mo8OHD5MoiqrwdLj4KN385V00Z/tXmvYUjveeKIp04MCBP7nU+jHp0dOVjkkURdq3b1/lfWp6TB6vlzoudFOvxT7yer3MUxiO6dy5c/ILLCKi3r1705w5c4iIKDc3l7p27XrJnU+dOkWiKFa+YM+ePenjjz8mIqIjR45Qo0aN6PTp0yRJEg0cOJDee++9Kl/8Yiqq2sLCQtmPYdSMu1f76OrPveTwSUF5vov/ImGEFtu6pVQ47jbynDwSlOdTm7up69+iO79+mJw+l9JRVI/a3DECI5j+Vp4QCTO99N0xMQjJGNVR0RCqzl3lWYQzZszAjBkzkJKSgtdffx1z5swBAIwaNQpLliwBAHz99dfo2LEjOnfujB49euDWW2/FiBEjAAAtW7bE5MmT0atXL7Ru3RoNGjTAmDFj5M1TXoTZbA74MQz55JdIWHCUMKWrgHhDcBbHms1mEBFzF2IkjxvWVV8gvtstiL46OGfWqcndwQtH8ePxXzC8472IM8QqHUf1qMkdI3CC6a9vEw59GnN4Pp9dQiccyHXGEanjstxWq7VywWadOnWUjqNLiAi3LBdR7CLsGGKAoQa7tl8OURSxf/9+tGvXjm16GEKsP86HdeXnaPjCbBjqNQzKc6rJ3bNrJuG0vRifDngfBp69j6pDTe4YgRNsf9vOEbp+68fM6wU82k51e4jritLSUtStWxcWi6XKYkt1FirOCmAEnx+KCD+fJrzWXQhacQWUO8vKymLuQojksMH2Uy4Set4etOIKUI+7bWd2YvPpAoxOf4gVVzJRiztGzQi2vy71OdzbisOkrSKcflX0TXSLXGeqK7BYuzs0SFR+SZxeyRwGNgvuvjkVFyxl7kKH9aeFgCTC1PfeoD6vGtwRET4qmIvUeim4vmkPxXJoDTW4Y9ScUPib2k3AWRfw3m62+WgoketMdQVWxV4TjOCy4Ahh+3ng9e580Dcm9Pv9yMvLY+5ChFh2DvZ1i5Fw450QTMGdPleDu7UnNuLAhcMYk/Ew2zQzANTgjlFzQuGvlZnDmPY8Xt8hodTDulihQq4z1RVYF1/bkBEcvCLhpS0iBjTjcF3D4Ct3uVzIyclh7kKE9YcvwUfFwHTT0KA/t9Lu/JIfs3Z8hh5Xd0V6cgdFMmgVpd0xakeo/L2cwcMnAa9vZ12sUCHXmeoKLJPJpHQE3TH7gIRjNmBaZmjWtphMJhQVFTF3IcBXchKOTSth6nM3+Dhj0J9faXfLDv+IU7Yz7ILONUBpd4zaESp/yfEcnu7I4909Ek46WBcrFMh1proCi7W7g4vDR5iyTcIDrTl0rBua6Re/349Vq1YxdyHAumIeBFNdJFw3MCTPr6Q7p8+Fubvmo2+LG9GqzjVhf32tw8adtgmlv2c78TAagMnbxKA/N0PDU4QV1w1iBIe3d0so9QBTanlB56pwu90YP348cxdkvEWH4Sr4Babb7gcXHROS11DSXe7+JbB57XikE7ugc01g407bhNKfOZrDixk8PjlA7ELQIUCuM9Xtg1XdvhIM+Zx3E1rO92NEWx5vZ7FT37XGuRkvw3/uFJKfnwlOZ/sclbktuHfxGAxo3Rd/6/qI0nEYDN3h9hPa5vrRvQGH3D4GpePoCrn1iuo6WBVXsmbUnte2S5AAvJgeWs0+nw+5ubnMXRDxHNkF9758mPs/HNLiSil3c3cvAMdxuD8t+Av3IwU27rRNqP3FGjhM6Spg0TFCfglb8B5M5DpTXYHl9XqVjqALTtgJ7++V8ExHHg3iQnvqu9frxfTp05m7IEFEsCybg6gmrRDX6bqQvpYS7k7aTmPJoVW4L3UIkmJZt7qmsHGnbcLh74HWHNLqAM9vlqCSySpdINeZ6gosozH4Z0pFIpO3ijBHAeM7hl6x0WhEXl4ecxck3Pvy4T22F4n9R4DjQ+tPCXezd3yBOjGJGNouNAv3IwU27rRNOPwJPIdpmQLWnCKsPskKrGAh15nqCiz211jtOVBG+PQQ4aUMHqbo0G/c6PV6MXv2bOYuCJAkwbp8LqJbdkBMu64hf71wu9t//hDW/L4OIzrdi1hDaBbuRwps3GmbcPkb2IxD1lUcXshnXaxgodkOFltPUHte3iKiiREY0z48etlakODh2rEevpNHkNh/eFh2NQ+nOyLCjIJ5aG5uguyWN4f89fQOG3faJlz+OI7DtEweW84Rvj3OCqxgINcZO4tQZ1RcUf2TGwSMaKu6+plRBSSKKH59NAz1r0b9MVOVjhN0Np/ahmd/noxpvV9ArybXKh2HwYgoslf4Uegg7BpigMCzS1LVBs2eRejxeJSOoGlezBfRLgl4sE34BpDH48H06dOZu1rizF8Nf8lJmG8P367m4XInkYQZ2+ehY4NU9GzcPaSvFSmwcadtwu3v1Uwe+8qAzw+roqeiaeQ6U12BJYps59ma8utpCSuLCFO7CjCE8S8UURSRl5fH3NUC8nthXfk54tKvR3TTNmF73XC5+/HYLzhcegyPsQs6Bw027rRNuP11a8BjSAsOE7eK8IisyKoNcp2xKUKdQES4fqkIl0jIH2wAzz7ENIXtl+9g+W4mkp+fgajkpkrHCSoe0YsHlzyBlHqt8MoN/1Q6DoMRsewrJXT42o93sng8maavzYvDCZsijDC+LyRsKCZM6yaEvbjyeDyYNGkSc1dDJI8Lth/nIz6zT9iLq3C4++7g9zjnOo/R6Q+G7DUiETbutI0S/trX4fBQGw6vFEhw+FTRW9Ekmp0ilCS242ygSER4IV/EDQ059G0S/s6VJEkoKipi7mqI/dfFkFwOmG+7P+yvHWp3No8dn+1eiP6tbkUzc5OQvEakwsadtlHK38QuAi54gHf3sPdNTZHrjE0R6oAFRyTcs0bE+oECejVUXc3MqALJacPpKcNhzLwFSUOeUDpO0PmoYC6+O7gCXwz6CPXi6igdh8FgAHhqo4jPDkk4eo8BdWLYcpJA0ewUIbsyfGD4JcLLW0T0b8opVlyF8qrwese2ZhEg+WHqe68irx9Kd8WOEny9fymGtb+DFVchgI07baOkvxfTeXgl4I2drItVE+Q6U12BxQiMuQcJh6zAK5lswaLWEK0XYP/1OyT0vhOCSX8FyCc7v4QxOh73tB+sdBQGg3ERyfEcxnXg8c5uCcVOVUxi6RI2RahhPCIhZaEfPa7isOAWg9JxGAFS9vUHcGxZg0YvzwEfb1I6TlA5UnocI1eMw9jM0bgz5Xal4zAYjL9Q6iG0mO/H8BQeb2exP9ADQbNThC6XS+kImmHWfglFDmByV2UHh8vlwqhRo5i7APBfKIZ94wqYbh6qaHEVKncfFcxFY1MjDGzdN6jPy/gfbNxpG6X91Ynh8ExHHh/ulVBoV0WfRTPIdaa6AovnVRdJlTj9hFcLJDzYmkO7JGUXKfI8jyZNmjB3AWBd+QX4uAQk3KDs9Fko3G09swObT2/D6PQHYeBZZzVUsHGnbdTgb2wHHuZo4JUCthYrEOQ6Y1OEGuWNHSJeyJdwcJgBLczsLBAt4SsuRPHrY5B05xgk3HCH0nGCikQSxqx8BgbegA/6/ovt2s5gqJy3dop4frOE/cMMaMU+S2Sh2SlCp9OpdATVY/USXt8hYVQ7XhXFldPpRE5ODnMnE+v38yAk1YexZz+lowTd3U/H1+HghSN4PGM4K65CDBt32kYt/p5I5dEgDpi8lV1ySS5ynamuwBIEttiuOt7eLcHhB17KUIc+QRCQlZXF3MnAW3gIru3rYL7tAXCGaKXjBNWdV/Rh9o7PcV2Ta9HpqtQgpGNUBRt32kYt/uIMHF7K4PH5YcLeUlVMaKkeuc7YFKHGuOAuP/NjZFse09mZH5rj3IyX4T9/GsnPzQCnsw/GhfsW48OCTzGn/7u4JlFf11NkMPSM948z0rs34LCwD1s3WR2anSJ0OBxKR1A1b+yUIBLwfLp61DkcDmRnZzN31eA5shvuffkw3/6QaoqrYLmzex34bHcubm/VhxVXYYKNO22jJn/RAoeJXQTkHiMUnFNFz0XVyHWmnk/pP4iKilI6gmopdhLe3SNhbAceV8WpZ31LVFQUcnJymLsqICJYls9BVJNWiOt0ndJxKgmWuy/3fgOP6MGIjvcEKRmjOti40zZq8/dgGw4picAEtharWuQ6Y1OEGmJcnohPD0o4xq4fpTlcezfj/MwJqDd6KuJSM5WOE1RKnOdx/5LHMKzdHRiV/oDScRgMRg2Zf0TCvWtEbBwkICtZdf0X1cCmCHVGkZ3w0T4J/9eRV11x5XA4kJWVxdxdAZIkWJd9iuiWHRDbvpvScf5EMNzN2fklYg2xuDftriAmY1QHG3faRo3+hrXk0KEOMGEr2xerKjQ7RRgdrfyZVWpk2nYJRkP5xnBqIzo6GuPHj2furoBr+6/wnTqKxAEjVLd1QW3dHbecwPdH1+ChDsNgjIoPcjpGVbBxp23U6I/nOEzuKmD1ScKvp1mRdSXkOmNThBrgdxuhzUI/pnbl8Vy6OhZHM+RBoh/Fr42G4aomqD96itJxgs4/176C45ZCzBvwPqIEdawlYTAYNYeI0PVbP0xRHNYOEFT3R6Ea0OwUod1uVzqC6nilQERSNPBkmup0ASh3lpaWxtxdBsdvP8B/7hTM/YcrHeWy1MbdjrN7sPFkPkZ1vp8VVwrAxp22Uas/juMwpauAX88Q1pxSRf9Fdch1prpP7NjYWKUjqIojVsKcg4TnO/MwRqnzL4nY2FhMnz6dufsL5PXAuuoLxHW5EdGNWyod57LU1B0R4aOCuWhbtxVuaq6esyIjCTbutI2a/fVvxqF7Aw4vb5GgkkkuVSHXmeoKLIOBbXJ2MVO2ibgqDng8VXWqKjEYDMjOzmbu/oJ9/VJItjIk9ntI6ShXpKbu1hVuwt5zBzAm42HwnHrfm3qGjTtto2Z/5V0sHnlnCSuLWIH1V+Q6U91vRpvNpnQE1XCgjPD5YcIL6TziDOrsXgHlzpo0acLcXYTkcsC2egGMPbJhaHC10nGuSE3c+SURM7fPQ/dGGejasHMI0zGqgo07baN2f32bcOiVzGEC62Jdglxnqiuw4uLilI6gGiZvE3F1PDCqreo0/Ym4uDjk5uYydxdhW/sNJJ8H5r73KR2lSmribuXRn1BoO4XR6ertzEUCbNxpG7X7q+hibTlHWHqCFVgXI9eZ6j651dguVYLdFwjzjxBeyuARq+LuFVDuLCsri7n7A9FeBvvab5Bw/SAISfWVjlMlgbrziF58ums+bm5+PdrUVee6skiBjTttowV/NzfmcWMjDhO2iJBYF6sSzU4RWq1WpSOogsnbRDRPAEakqE7RJVitVpjNZubuD2w/LgA4DqZbhikdpVoCdffdwe9xwVWKkZ3U3ZmLBNi40zZa8TelK48dF4Bvj7MCqwK5zlT36W00GpWOoDg7zhMWHSO83EVAtKDu7hVQ7iwvL4+5A+C/UAz7+mUw3TwUglH9+7kF4s7hc+KLPbm4vdWtaGJW77qySIGNO22jFX/XN+Jxa2MOE7eKECVWZAHy6xTVFViCwDbSnLhVRCsz8FAb9RdXQLmztLQ05g6AdeUX4OOMSOh9p9JRZBGIuwX7voPL78HDHe8OQzJGdbBxp2205G9KVx57SoHcY6zAAuTXKaorsNTeLg01284RFv9OmJAhwMBro8CyWq3gOC7i3fnOnIAzfzXMfe8FH6POhat/Ra67MrcFC/ctxl0pt6NBfL0wpWNUBRt32kZL/nok8+jXlMOUbayLBWh4ijAhIUHpCIoyaauINmbgvtbaKK6AcmeFhYUR7876/TwISfVh7NlP6Siykevu8z254Dke96UNCVMyRnWwcadttOZvUhce+8qAhUdZgSXXmeoKrEi+7tGWEglLTxAmdNFO9wood2Y2myPanffEAbh2rIe534PgDOq5eGt1yHFX7CjBdwe/x93tByMxRv3ryiIFNu60jdb8db+KR/+mHKYUsC6WXGeqK7DUuulaOJi0VULbRODeVtoYcBXYbDYkJiZGtDvLsk9haNgM8d1uVjpKQMhx9+mu+TBGGZHTbmAYkzGqg407baNFf5O68thfBsyP8C6WZjcaNZlMSkdQhM1nJSwvLO9eCRrqXgHlziwWS8S6cx/cDs/BAiTe/jA4Xv0LVi+mOne/W4qw8ugaPNghB/FR8WFOx6iKSB93WkeL/ro14DGwWflaLH8Ed7HkOlNdgRWpW/JP2iahXRJwd0ttFVdAuTOr1RqR7ogIlmVzENWsLWI79lQ6TsBU527Ozi/RIK4eBrW5LczJGNURyeNOD2jV38QuAg5agK+OaCt3MJHrTHUFlt1uVzpC2NlULOH7wvIzB7XWvQLKnTVt2jQi3bl3bYTvxAEkDhihmbUUF1OVu6Nlv2PtiY14qOMwRAtRCqRjVEUkjzs9oFV/XRtwGNQ8srtYcp2prsAymyNvEe2kbRLaJwHDNNi9AsqdEVHEuSNJhGXFXMSkZCA2JV3pODWiKndzdy1AsrEBslvcpEAyRnVE6rjTC1r2N6mLgMNW4IvDkVlgyXWmugJLFEWlI4SVvGIJq4oIEzW49qoCURSxZ8+eiHPn3LYW/jMnkNh/uNJRasyV3JV3rzbgwQ45iGLdK1USqeNOL2jZX0Z9DoObc5haEJldLLnOVFdgORwOpSOElUlbJaTVAXI02r0Cyp1lZWVFlDsS/bB+/zliO2QhunlbpePUmCu5m7trARoar8JtLbV1VmQkEYnjTk9o3d+krgKOWIHPDkVegSXXmeoKLC22S2vKxmIJP5ws717xGly/U0HFBUsjyZ1j0yqIF84g8faHlI5SKy7n7kjpcaw9sQEPdRgGAy/vqvGM8BOJ405PaN1f53oc7rqGw6vbI6+LpdkpQr/fr3SEsDF5q4QOdYAhLbRbXAHlzvLy8iLGHfm8sP3wFeIyeiPq6hZKx6kVl3M3d3d59yq7JVt7pWYibdzpDT34m9ilvIv1eYStxZLrTHUFlsvlUjpCWMj7o3s1QePdK6DcWU5OTsS4s29YBtF2AeZ+Dyodpdb81d2R0uP45cRG1r3SAJE27vSGHvx1+qOL9UqErcWS60x1BZaWNl2rDZO3la+90nr3Cih3VlRUFBHuJI8LttULEN/9VkQ1aKx0nFrzV3dzdy9Ao4Rk1r3SAJE07vSIXvxN+KOLFUlnFGp2o1Ett0vlsumPMwcnZGi/ewWUO1u1alVEuLP/8h0klxPm7PuVjhIULnZX0b16kHWvNEEkjTs9ohd/netxuPOayDqjULNThG63W+kIIWfyNgmpScBQDZ85eDFutxvjx4/XvTvJaYNtzSIk9OwHQ52rlI4TFC52N3fX/PLuVYsblY7FkEGkjDu9oid/FV2sLyOkiyXXmeoKrISEBKUjhJTfzkpYWaSPtVcVJCQkYM+ePbp3Z/v5a0D0w3TrPUpHCRoV7s76z+OXwjzWvdIQkTLu9Iqe/KXXi6x9seQ6U12B5fP5lI4QUib/sWv7UB2svarA5/MhNzdX1+5EWynsv3yLhBvugGCuq3ScoFHh7ovdi5Ac34B1rzREJIw7PaM3fxP+2N09Eq5RKNeZ6gosr9erdISQsfnsH9cc1PCu7ZfD6/Vi+vTpunZnW7MI4AWYbs5ROkpQ8Xq9mD7rbaw5sR53pw5m3SsNEQnjTs/ozV9GfQ53NOcwNQKuUSjXmeoKLKPRqHSEkDF5m4R2SUCOjrpXQLmzvLw83boTbaVwrF+GhBsGgzdq+4yfv2I0GnH3qw/DFJ2A/q1uVToOIwD0Pu70jh79Tegi4JAVmK/zLpZcZ6orsPRSzf+V/BIJKwrLzxzUU/cKKHc2e/Zs3bqzrVkECAJMve9UOkrQKbaWYMmB7zG4dT/EGmKUjsMIAL2PO72jR39d6nMYFAFrsTTbwdLLfPRfmbJNQttEYJhOzhy8GL2tJbgY0Vb2R/fqDt11rwDgmwPLIPpF3N78FqWjMAJEz+MuEtCrv4ldBBy06LuLpdk1WHpql1aw7Rxh2QnCyzpbe1WB0WjEqlWrdOnOtia3fO1V77uUjhJ0nD4nlh3/EcM6DkbDOslKx2EEiJ7HXSSgV39d6nMY0Kz8GoWiTrsrxAqsAAAgAElEQVRYmp0i9Hg8SkcIOlO2iWhjBu7WYfcKKHc2ffp03bmr7F711mf3aunhH+D2uVG2rlh37iIBvY67SEHP/iZ04bG/DMg9ps8CS64z1RVYoigqHSGobD9PWPw74aUMAQYddq+Acmd5eXm6c2f7eZFuu1de0YeF+xbj5mbXY/uGbbpzFwnoddxFCnr2l9mAR7+m5WcUSqS/IkuuM45IHUdvtVqRmJgIi8UCs9msdJygMXS1HwXnCAeGGXRbYOkR0VaGM1MfRkLvO5HYf7jScYLO8sM/4o3f/ou5A95H88QmSsdhMBg6I69YQs8lInJvETC0pep6ObVCbr2iuqPWU7t09wXC18cIL+q4ewWUO5s0aZKu3Nl+XgRwAhJu1F/3SpREfLX3G1zftAcaxjbQnbtIQY/jLpLQu7+sZB63Ni4/o1BvXSzNThFKkqR0hKDxSoGI5gnAg230W1wB5c6Kiop04060l8GxfikSbhgEwaifbmoF64t+Q6HtFO5LvUt37iIJ5k7bRIK/CV147LwALPldXwWWXGdsijBE7C0ldFjkx0fXCRjdXnV1LKMKypZ8DMeGZWg4Ya7uCiwiwmOrnkWcIRZv93lF6TgMBkPn3LzMjzIvYeudBnA6uf6uZqcI9XBlcQB4tUBEEyMwPEUfb6iq0NNV4UWHtbx7db0+u1e7S/Zj//lDuKd9+aapenIXaTB32iZS/E3syqPgPLDshCp6OUFBrjPVFVh64EAZYf5Rwj/TeUQL+i+w9IRj3VKACAk36m/XdgDIPbAEzcyN0f3qDKWjMBiMCKB3Ix43NOQwZZsElUyYhQ02RRgCHl7rx0+nCEfuNiCGFViaQfK4cWbKQ4jrciPqDHlC6ThB54z9LO5dMgZju43G4JR+SsdhMBgRwk8nJfRZIWLFbQL6NdV+X0ezU4Qul0vpCLXiiJXwxWHCc534iCmuXC4XRo0apXl3jk0rIbnsMOnwzEEA+Pbgchij4pHd8qbK/9OLu0iEudM2keTv5qs59EzmMFknXSy5zlRXYPG86iIFxLQCEQ1igVHttH0cgcDzPJo0aaJpdyT6YV/7DeK73AhDvYZKxwk6Tp8Lyw7/iIGt+yLOEFv5/3pwF6kwd9omkvxxHIcJGTx+O0tYfVL7BZZcZ2yKMIgctxHaLPDjX915jO8kKB2HEQCO/NUo/eJNJP/jQ0Rd3ULpOEHnmwPL8f7W2Zh/x0xcZWygdBwGgxFhEBGuXSwiVgB+HWhQOk6t0OwUodPpVDpCjfnXDglJMcCYCNuWwel0IicnR7PuSJJg+ykXsanddVlcSSRh0YGl6N2s5yXFldbdRTLMnbaJNH8cx+HlDB7rzhB+Oa3tvb/kOlNdJSAI2uz8nHQQPjkg4f868jBGRcbaqwoEQUBWVpZm3bn3bob/zO8w9RmmdJSQsOnkVpy0ncbQdoMu+Z7W3UUyzJ22iUR/A5pxSK8HTN2m7QJLrjM2RRgkxm4U8dlhCb/fY4ApOrIKLK1z9p3xAICrxk5XOEloeHr1y3D73fjwtjeUjsJgMCKcr49JGLpaxMZBArKSVdfjkYVmpwgdDofSEQLmjJMwc7+EcR34iCyuHA4HsrOzNenOc3Q3vMf2wtTnbqWjhIQjpcexrXgnhrYbeNnva9ldpMPcaZtI9XfnNRzS6gBTC7TbxZLrTHUFVlRUlNIRAuatXRKieeCpNNX9OMNCVFQUcnJyNOnOtnohDA2bI7Z9ptJRQsKiA0vRIL4eejfrednva9ldpMPcaZtI9cdzHF5MF/B9IWFLiTaLLLnO2BRhLTnnJlzzlR/jOvB4JTNy5tL1gO/UMRT/+3HUuf8ZGDP7KB0n6JS6yzDs21EY0ele3Jc2ROk4DAaDAQAQJULqIj/aJ3H4rq/2zigMeIrw0KFD6NmzJ1JSUpCZmYk9e/Zccuc1a9age/fuSE1NRVpaGv7xj39UXlX6+PHjEAQB6enplV9HjhwJOLjW2qX/2VV+/OM6Rmb3Cih3lpWVpTl3tp+/hlDnKsR3uVHpKCFh6aEfwHEcBrTue8X7aNUdg7nTOpHsT+A5vJAuYPHvhJ3nVdHjCYiApwjHjBmD0aNH4+DBg3juuecwfPjwS+5cp04dzJ8/H3v37sXWrVuxceNGzJs3r/L7JpMJ27dvr/xq1apVwMGjo6MDfoxSlHoI7+2R8EQqj/qxkbf2qoLo6GiMHz9eU+5E6wU4t61Fwg13gBO09xdUdfglEUsOr0Sfa3rDHGO64v206I5RDnOnbSLd332tObQwAa8UiEpHCRi5zngAOHv2LLZs2YIHHngAADBkyBAUFhbi8OHDf7pzRkYGWrZsCQCIjY1Feno6jh8/HsTY2lqD9e5uCX4J+L8I7l4B2lxLYN+wHJxggPHabKWjhISNRZtR4jyPO1Nur/J+WnTHKIe50zaR7i+K5/DPdAGLjhH2lWqriyXXGQ8AhYWFaNSoEQyG8r/kOY5Ds2bNcOLEiSs+8MyZM1i0aBEGDBhQ+X8OhwOZmZno0qULpkyZAlG8cmXq8XhgtVr/9AUA586dAwC43W643W4A5df98Xg8AMo3+Kq47XA44PV6K2/7fD4AgN1uh9/vBwDYbLbK21artTKT1WqFJJVfF8lqtYKIIElSZQ5RFCtv+/1+2Gy2ytt2ux1WL+Gd3RKGtxKRHM/B6/VWtg09Hk/lRmQej6fyukVqPyYA8Pl8lcch95jOnTuHtLQ0lJSUaOKY7GWlcGxYhphut8D9RxNXb56+PbgCqfXaoklcoyqP6ezZs2jfvj3sdrvqj0mPnmpzTHa7He3bt0dpaalujkmPnq50THa7He3atUNZWZlujilQTw+2IjSOl/DqdlFTx2SxWCALIqItW7ZQSkoKXUxmZib99NNPdDksFgt169aN3nrrrcr/c7vdVFxcTERE58+fpz59+tC//vWvyz6eiGjixIkE4JKv+++/n4iInn76aXr66aeJiGjkyJE0ceJEIiIaOnRo5ev27duXZs2aRUREPXr0oIULFxIRUWpqKq1cuZKIiBo3bkwbN24kIiKTyUS7d++mPxb2U2FhIVksFgJAFouFCgsL6Y8fCe3evZtMJhMREW3cuJEaN25MREQrV66k1NRUer3AT4aZLsq4eSAREc2aNYv69u1LRERvvfUWDR06tPI4R44cqYljIiJauHAh9ejRI6BjGjt2LK1cuZJGjBihiWOanNOXCsdm04zXp+jS0+DhQ+iGzwfRY2/+vdpjGjJkCD366KPk8/lUfUx69FTbY/L5fNSuXTv66quvdHNMevR0pWPy+XxUv359+vXXX3VzTDXydNPjxM/y0ootBzRzTOvXr688hqoAEVFxcTGZTCby+XxERCRJEiUnJ9OhQ4cueYDVaqWsrCyaOnVqlU/85Zdf0oABA674fbfbTRaLpfKr4s1z5swZIiJyuVzkcrmIiMjpdJLb7SYiIofDUXnbbreTx+OpvO31eomIyGazVR6L1WqtvG2xWMjv91feFkWRJEkii8VCkiSRKIqVPzC/31952+fzkdVqrbxdXGqjBvO89OhaL9ntdiIi8ng8lbfdbjc5HI7K206nU/XHZLPZiIjI69X3Mbndbjr12mgqmTlBN8f0V09vbvwvDcp9kKwOm26OSY+e2DGxY2LHJNKZ8xZK/sxLI9d6NXNMFy5ckF9gERH17t2b5syZQ0REubm51LVr10vubLPZqGfPnjR58uRLvldcXFwZzO1209ChQ+nll1+u8sUvpqKqLSoqkv0YpXh7l5+EWV46apGUjqIKrFYrNW7cuHJAqBnXwe1UODabXAe2KR0lJDi8Tuq34B6aWTBP1v215I7xZ5g7bcP8/Y83dvgparaXfrdp4zO1qKhIVoFVuTp7xowZmDFjBlJSUvD6669jzpw5AIBRo0ZhyZIlAIB33nkHmzdvxjfffFO5FcOrr74KAFi/fj0yMjLQuXNndOnSBQ0bNsSLL74ob57yIuLi4gJ+TDjxiIR/75DwQGsOLcyRe+bgxcTFxSE3N1f17gDA/su3MDS6BjFt0pWOEhJ+PLYWLr8bg9rcJuv+WnLH+DPMnbZh/v7HY+15mKKAN3ZoY+NRuc7YRqMBMmOfiMfXS9iXY0DbJFZgaQl/ySmcmTYSde4eC2MPeQWIliAijFg+Fo1NDfFq7xeUjsNgMBiyeWWbiFe2Szh+jwEN49X92arZaxFWnEWgRnwS4fXtEoa15FhxdRFWqxVms1nV7gDAvm4x+Hgz4rvcpHSUkLCzZC+OWX6vdmuGi9GKO8alMHfahvn7M0+m8Yjhyy89p3bkOlNdgWU0GpWOcEW+PEw4bgdeSGeXxLkYo9GIvLw8VbuTXA44fvsBxl63g4uOUTpOSPju4Ao0MzdG14adZT9GC+4Yl4e50zbM359JiuHw9zQeH+6VcM6tiom1KyLXmeoKLEFQZ/EiSoRp20Xc0ZxDp3qse3UxgiAgLS1Nte4AwPHbKpDfh4ReA6q/swY577qAX07k4Y42/cBx8t+fWnDHuDzMnbZh/i5lXEceBOCd3eruYsl1proCS63t0kXHCActwIsZqvuRKY7VagXHcap1R5II+6+LEZ9xA4TEekrHCQlLD/+AKN6A7JaBTX+q3R3jyjB32ob5u5T6sRweb8/j3d0Syjzq7WJpdoowISFB6QiXIBHhlQIR2U04ZDZQ3Y9McRISElBYWKhKdwDg3rMZ4oViJNwwWOkoIcEviVh26Afc2uJGmKIDc6B2d4wrw9xpG+bv8vxfJ/7/27vz8CiqtG3gd1V31k53s4UQskFIghAgkU0CKKAIqCBuKLiM6LCojDhmFEGdV3ADRVHHZUBF/fwYBNlBEHxVRJagMiRAgoQAIQv7lnR676467x+ZZABDUkm6u041z++65jKhq6ueyg3DwzmnT8ElAx/u53cUS2lm3HULjZneCJR1JQz5F4AXafSqToIgwGQycZkdAFi3rUNo0jUITUxTuxS/2FH+K844zuGOtFsa/V7esyNXRtlpG+VXt9hIARM6i3hnnwyrh89RLKWZcdcx1Dx/iBeMMbyaK+OGdgIGtuPux8WFqqoqmM1m7rIDAM/pcrgKd8MwMDjXXgHA6qJvkd6mM1Jadmz0e3nOjtSPstM2yu/KpmWIqHQDC37ncxRLaWbcdQxGo1HtEi7xv8cYdp1lNHpVD6PRiMrKSu6yAwDb9vUQDWZEZt6gdil+UW45jn+f3IPRqY0fvQL4zo7Uj7LTNsrvyhKjBPwpVcBbe2U4vfyNYinNjLuugZN9T2u9liujT7SAoXE0jHsl7KKno/NEdjmrt2boNxxCSKja5fjF2kObYA4zYnDSgCa9n9fsSMMoO22j/Oo3PVOH007g84P8jWIpzYy7BstqtapdQq1tJ2X8fJLhhUyR5snrYbVakZCQwFV2AODYvRnMZYehv/KNN7XE5XXh28M/4JbkmxCma1oDyWt2pGGUnbZRfvVLNQu4N1nAG3tkeGS+mlClmXHXYPH0mJzXcmV0awmMSqLmqj4mkwmMMa6yY4zBuu0bhHftC33rdmqX4xebS7fD4q7CqJThTT4Hj9kRZSg7baP8GvZ8pg4l1upNvnmiNDPuGixJktQuAQDw7zMMG8sZns/UQaTRq3pJkoSCggJusgMA99Hf4Tl2GFEDR6ldit+sKdqI3u0yEW9q3+Rz8JgdUYay0zbKr2HdWwm4PUnA7DwJEkejWEoz467BstlsapcAAHgtT0KKCbg3mZqrhthsNmRlZXGTHQDYtq2Drk0swjr3VLsUvyg6fwT7zxZidFrzHlrNY3ZEGcpO2yg/ZV7IFFFYCaw8yk+DpTQzvZ/raDQehksLzjOsOsrw6fU66ERqsBrC2wNLpaoK2PO2wTxyPASRu39D+MSaoo2IjmiN/nF9m3Ue3rIjylF22kb5KdO3rYihcTJey5VwT0eBi/XQmp0i9Hq9apeA2XskJBiAh1LVD1ILvF4vcnJyuMgOAGw7NwKiAEPfYWqX4hc2jx3fH92CkSnDoBeb9xwz3rIjylF22kb5KfdCpog954ENZXyMYinNjLsGy+FwqHr9wxaGrw4zPNtDRKiOGiwlHA4HxowZo3p2QPVzB207NiDy2sEQDcG5v8x3RzbDLblxW8rNzT4XT9mRxqHstI3yU25QrID+MQJey5W52NZCaWbcNVhqb7r2xh4JbcKBCddw96PhltFoRHl5uerZAYBz/6+QLpxGVJDu3M4Yw5qijegf3xfRkc1/cDVP2ZHGoey0jfJTThAEvJApIuc0w08n1G+wNLvRqJrDpeVWhi8OMmR3ExGhp9ErpbxeLzZt2sTFULd16zqEJHYO2ucO7j2zH8WVpU3euf1yPGVHGoey0zbKr3FuSRBwbevq7ZPUptkpQqfTqdq1394nw6AHHu/K3Y+Fa06nE9nZ2apmBwCeM8fgKtyNqAG3qVqHP605uBFxxlj0atfDJ+fjJTvSeJSdtlF+jSMIAl64VocfjjP8clrdJktpZtx1ElFRUapc94yDYcHvMqZ2E2EKpdGrxoiKikJBQYFq2dWw7dgAMdKIyGsHqVqHv5x3VGBL2Q6MTr0FouCbP7q8ZEcaj7LTNsqv8e7sIOCaFsDreeo2WEoz467B8ng8qlz3vXwZogBMTefuR8I9j8eDZcuWqZYdADC3C7ZfNiHyumEQQsNUq8Ofvjn0HXSCiFuTb/LZOXnIjjQNZadtlF/jiYKA6Rk6rC1h2HdevbVYSjPjrptwu90Bv2alm+GD/TIe6yKidTiNXjWW2+3GvHnzVMmuhj3vZzC7NWifO+iVJaw9tBFDO9wAY5jv/sXLQ3akaSg7baP8mub+FAFJUcDsPPV2wFeaGXcNlsFgCPg1P9ovw+EFsrtz9+PQBIPBgJycHFWyq2Hbvh5hnXsiJDpOtRr8aUf5rzhjP4c70ny7voyH7EjTUHbaRvk1TYgoYFqGiKVHGA5VqjOKpTQz7jqKQHfzdi/DO/tkPJImor2BRq+awu1249NPP1XtX2LusiK4Sw4gakBwbs0AACsPrke3NtcgrVWyT8+rdnak6Sg7baP8mu6RNBHR4cCbe9UZxdLsCFag56MXHpBx3gVMy+DuR6EZaq8lsG1fD12LNghPv06V6/vb0cpS5J7ahzvSfD/9qXZ2pOkoO22j/JouQi/gb91FfHGQodwa+FEspZkJjIdtUQFYLBaYzWZUVlYG7HmEbokhZakXg2IF/P8h3D2WkSgg26048dIDMA69F6bhD6hdjl+8+9sC/FS6HV/fsRChuhC1yyGEENVVuRkSv/JifJqId7Ka98iwxlLar3A3bONyuQJ2rUWHGMpswPSMwIYTbFwuF+bNmxfQ7GrYfvseTPLC0G9EwK8dCHaPHZuObMaolOF+aa7UzI40D2WnbZRf8xhDBUztJmLB7zLOOAI7TqQ0M+4aLEkKzJyqJDPMyZNwZwcB6a1o7VVzSJKEnJycgGVXgzEG2/b1iOjRHzpz8x8bw6NNxT/BJbkwKnW4X86vVnak+Sg7baP8mm9qughRqN5mKZCUZnbVThF+fVjGfT9K+PUOHfpEc9dnEgWcRXtw9sPn0GbKGwhPzVC7HJ9jjGH8+ieRaIrHKzdMV7scQgjhzjM7JXxaKKNknB7mAG0STlOE9WCM4fU8CUPjBGqufMDlcmHmzJkBH+q2bVsHfUwiwlJ889gY3uSe2oejlWW40w+L22uolR1pPspO2yg/38juLsLhrd5uKVA0O0Uoy/7/IX1bxrDnPPBCJne3r0myLKO8vDwg2dWQKs/BsS8HUQNugyAE5xTv6oMb0MGcgGtjuvvtGmpkR3yDstM2ys832hsEPJIm4p19MuzewEzIKc3sqpsiZIxh4DoJjAHbb9cF7V/Owc6ycRGqflyG2FmLIUYE30Z9p+1nMXb1RDzZawLu7By8D68mhJDmOmJhSP3ai/eyRPwl3f8fWtPsFKG/nyy+9STDjlMMz2eK1Fz5SKCfCs+8Hli3r0dk76FB2VwBwKrC9QjXh2NY8hC/XifQ2RHfoey0jfLznWSTgHGdBMzdK8Mj+3/MSGlm3DVY/vZ6nozurYDbEqm50ipH3lbIVRcQdcPtapfiF3aPA2sPbcLIlGEwhESqXQ4hhHBveoYOpVZg8SEuJuUAXGVThLvPMvRa5cXiITqMS7nqesugcWreUxDDIxH9xGy1S/GL5QfW4aPdn+Gr0R8jxhCtdjmEEKIJo7/zorCCoeAePXSi/wZRNDtF6HA4/Hbu2XkSOpmAMck0euVLDocDEyZM8Gt2NVxHD8BTWoioG0b7/VpqkGQJywvXYXDiwIA0V4HMjvgWZadtlJ/vzcgQUVgJrC7x77iR0sy4a7BE0T8lHahgWFHM8FyGDno/drZXI1EUER8f77fsLmbduga61rEI79rH79dSw9bynThhPYX7ugSmgQxkdsS3KDtto/x8r1+MiBvbC3g9T4I/J+eUZnbVTBE+ssWL78oZjozVI0xHDZYWSZXncOLlh2Ee9SiMg+9Suxy/eGLTNISIIXjv5tfULoUQQjTn+2Mybt4gYdMtOgyL90/zqtkpQrvd7vNzllQxLCpi+Ft3kZorP7Db7RgzZoxfsruYdccGCDo9DH2H+fU6ask/cwAFZwtxb4BGr4DAZUd8j7LTNsrPP25qL6BPtIDX8/y3v5jSzLhrsHQ63+9h8dZeGaZQYFIX7m43KOh0OmRlZfkluxrM64ZtxwZE9h0KMTLKb9dR09e/r0aiKQ5Zcb0Dds1AZEf8g7LTNsrPPwRBwPOZIracYNhxyj9NltLMgn6K8LSDIekrL6ZniHipF/1G1irbrh9wYdFcxEz/GCHtEtUux+eOV53EA+sex9N9JuP21BFql0MIIZolM4buy73oaBTwzQi9z8+v2SlCm83m0/O9my9DJwBPpnN3q0HDZrNh+PDhPs+uBmMM1i2rEda5Z1A2VwCwrHAtjKEGDO/o341FL+fv7Ij/UHbaRvn5jygImJ6pw/oyhr3nfD+GpDQz7rqOkJAQn52r0s3wYYGMx7qIaBVOa6/8JSQkBGPGjPFpdhdzlxyAp6woaLdmqHJZ8e3hHzA69RaE6cMCem1/Z0f8h7LTNsrPv8Z2EpAUBczZI/n83EozC+opwtl5Emb+W0bxWD3aG6jB0qpzX86Bu7QQ7Z5fCCEIP9L8r4Ll+GLvEiy941O0imihdjmEEBIUPiyQMDVHRuEYPVLMvusBrvopQoeX4d18GePTRGqu/MxmsyErK8svQ93eC2fgyNuKqOtHB2Vz5fK6sKLwG9zccZAqzZU/syP+RdlpG+Xnf492FtEmHJi717eL3TU7RRgaGuqT83xWKOOsE5iWwd0tBp3Q0FBkZ2f7LLuLWTcvhxgeCUO/4T4/Nw/WHtqECmcl7k+/R5Xr+zM74l+UnbZRfv4XoRfwdDcRXxyUcdzmu8k6pZkF5RShR2ZIXepFVoyAr270/ScISGBIVRU4+fLDMN50L0wjHlC7HJ9zel0Yt2YSrmvfC9OzpqpdDiGEBJ1KN0PiYi8mdREx9zrf7CSg2SlCq9Xa7HMsOcxQYgVmZNC2DIFgtVqRnp7uk+wuOe+WVYAoIur62316Xl6sKfoWla4q/KnbvarV4K/siP9RdtpG+QWGOVTAlHQR/9wv47zTN+NJSjPjrsEKDw9v1vtlxjAnT8JtCQJ6tKa1V4EQHh6OefPmNTu7i8l2K6zb1sEw4DaIBqPPzssLh9eJxQUrMSL5RrQ3tlOtDn9kRwKDstM2yi9w/tpNhMSAD/b7Zi2W0sy4a7D0+uZN6a0tYdhfAczI5O7WgpZer8fw4cObnd3FrNu/AfN6YBx0p8/OyZNVhRtg9djwULcxqtbhj+xIYFB22kb5BU7bCAETOov4R74Mm6f5o1hKM+OuC6mqqmryexljmJ0n44Z2Aga04+7WglZVVRXi4+Obld3FZLcT1p9WwXDdcOjMrX1yTp7YPXYs+X0lbk2+CbFRMarW4uvsSOBQdtpG+QXWMz1EVLqBTw40fxRLaWbcdSERERFNfu/m4wy/nmE0ehVgERERWLZsWbOyu5gtZyNkRxWMN6rzyTp/W1m4AXaPAw+qPHoF+D47EjiUnbZRfoGVZBRwf4qAt/bJcEnNG8VSmhl3nUhzhktn75FxbWtgeDytvQokvV6PrKwsnwx1M68H1s3LEdlzCPSt1Vub5C82jx1Lf1+NkSnDEGOIVrscn2ZHAouy0zbKL/Cey9DhmA1YVNS8BkuzU4QWi6VJ79t1Rsb3xximZ+ggCNRgBZLFYoHJZGpydhez7/oRUsVZGIeq98k6f1px4Bs4vU48oNK+V5fzZXYksCg7baP8Aq9rSwF3JAl4Y48ESW56k6U0M+4aLIPB0KT3zc6TkWoC7u5IzVWgGQwG5OTkNDm7GkyWUPXD1wjv0R8h7ZJ8VB0/rG4blh5YjVGpwxEdycfaMl9lRwKPstM2yk8dMzJFFFmAlUeb3mApzYy7Bkuna/zeVQcqGFYdZZiWoYNOpAYr0HQ6HdLT05uU3cUce7bDe+YYTEPv81FlfFmyfxXckgcPpN+tdim1fJUdCTzKTtsoP3X0bSvixvYCZudJaOo+60oz467Baspw6Rt7JMRGAg+lUnOlBovFAkEQmjXUzSQvLN9+ibDOPRGa2NmH1fGh3HIcS35fhbFd7kDriFZql1PLF9kRdVB22kb5qWdGpojcc8B35U1rsDQ7RRgVFdWo48usDIuKGP7WXUSYjhosNURFRaGsrKzR2V3Mtn09vGeOwXz7BB9WxgfGGN7b9TFaR7TiZu1VDV9kR9RB2Wkb5aeem9oL6N1GwJw9TduyQWlm3DVYjV2g/vY+GaZQYFIX7m7lqiEIAkwmU5M/XCDbq2DZuAiR1w1DaFyyj6tT39bynfj1RC6m9p6AcH2Y2uVcornZEfVQdtpG+alHEATMyBTx0wmGnaca32QpzQ8BB7kAABM6SURBVIy7rqQxm66ddTJ8ckDGk+kiokLoN6laqqqqYDabm7xhnmXTYjDJC/OtD/u4MvU5vE58sGshstr3Rv+4vmqX8wfNzY6oh7LTNspPXXd0EHBNi+rtnRpLsxuNGo3Knzv3j/zqH8yT6dzdxlXFaDSisrKyUdnV8Jwuh3XrWhiH3gediZ+1Sb6yKH8ZLjgr8GTvCVz+S7U52RF1UXbaRvmpSxQEPJehw9oShvzzjVuLpTQz7joTpav6q9wM7xfImHSNiNbh/P3FdTVhjMFisTTpExmVaxdCZ24dlM8cLLMcw5LfV+P+9LsRZ4xVu5w6NSc7oi7KTtsoP/Xd30lAgqH6g3KNoTQz7hosq9Wq6LgFB2TYvMDfunN3C1cdq9WKhIQExdnVcBbtgTM/B+aRj0AI5WttUnPVLGyPjmyN+7vepXY5V9TU7Ij6KDtto/zUF6oT8EwPEV8dZii2KG90lWbGXXdiMpkaPMYlMczbJ+OhFAHxUTR6pTaTyQTGmKLsajBZQuWqBQhNugYRPQf7rziV/FyWg99O5GFqr4kI42xh+8Wakh3hA2WnbZQfHyZcI6JlGPDWPuVrsZRmxl2DJUkND9V9WcRw0g5My6AN2nggSRIKCgoUZVfD/uv38Bw/AvOdk7lcm9QcNo8dH/x7IfrH9UH/+D5ql1OvpmRH+EDZaRvlx4dIvYCn0kV8VijjlF3ZKJbSzLhrsGw2W72vSzLDm3sk3NVRQOcWwfUXs1bZbDZkZWU1mF0N2WlH5YYvEHHtIIR16OLn6gKLMYY3d34Am8eOqb0nql1OgxqbHeEHZadtlB8/pqSL0IvAewXKRrGUZsZdg9XQ0NvyYoZDFmAGjV5xo+aBpUqHTStWzQdzOmAe9aifKwu81Qc34KfS7Xiu35OIjYpRu5wGNTY7wg/KTtsoP360DBPwWBcRHxbIqHQ3PIql2SlCr9d7xdcYY5izR8LNcQJ6RdPoFS+8Xi9ycnLqza6GPfdn2H/5Di3ufgL6Vvw3II1x4FwRPtj9Ge7uPBKDEvurXY4ijcmO8IWy0zbKjy9PdxPhlID5+xsexVKaGXcNlsPhuOJrm8oZ8s5VP0eI8MPhcGDMmDH1ZgcA3vOncGHpe4jIvAGRfW8OUHWBUeWy4qWtbyKlZUc8fu14tctRTGl2hD+UnbZRfnxpbxAwPk3EO/kyHN76R7GUZiYwTjbhsFgsMJvNqKysvOLw26B1XjglYOdoXdAtjA52TJJw5sNpkC6cQcyzH0KMDJ7N9RhjePHn2cg7nY9Pb3lHE1ODhBBCLnWokqHzMi8+6C/i8a5XXoakpF8BOBzButLQ245TMn4+yTAjU6TmijNerxebNm2qd9i06n+XwF38O1o9NC2omisAWHZgLbaV/4Lns57SXHOlJDvCJ8pO2yg//qSYBYzpKGDuXhle+cpjT5qdInQ6nXX++pw8GV1aALcnUXPFG6fTiezs7Ctm5zpSAMumf8E4bBzCkrsFuDr/KjhbiPm5/w/3dbkDA+KvU7ucRmsoO8Ivyk7bKD8+PZehQ3EV8PWRKzdYSjPTxBRh/nmG7iu8+GKQDg+ncdcTknrIditOzX0CuhZtEP2XuRB0wfPpz+KKUjz9w4toH9UO/7j5dehFvdolEUIIaaZbvvXimJ1hz136OmfMNDtF6PF4/vBrc/ZISDAA96fQ6BWPPB4Pli1b9ofsmCzjwtJ3ITtsaPXgtKBsrlqFt8Trg17UbHN1pewI/yg7baP8+DUjU8S+88D60rrHn5Rmxl2D5Xa7L/m+2MKw5DDDMz1EhIjUYPHI7XZj3rx5l2THZBkXvn4Pjr070Grc09C3bqdihb51cXM176ZX0CJcu/vY1JUd0QbKTtsoP35d305AVlsBc/bUvWWD0sy4nyKcsl3C10dklIzTI1JPDZYWMFlGxbL3Ydu5ES3v/xsMfYaqXZLPBFNzRQghpG7rSmTc/p2En0fqcH3spWNRmp0ivLgzPGVn+KxQxtR0kZorjrndbnz66adwu91gjKFixUfVzdXYp6m54tzF2RFtoey0jfLj222JAtJbArPrGMVSmhl3DdbFc5vvFcjQi8CUrtyVSS5Ss5bA7XajYuU/Ydv+DVre9xQM1w1TuzSf+f3swaBrrgBaB6JllJ22UX58EwUB0zN0+LaMYc+5Syf6lGbG7RRhpZshcbEXk7qImHtd8CyODlaMMVSuWgDrz6vR4t6piOp/q9ol+YQkS1i8fyU+3/sV0lolY87g/wma5ooQQsiVeWWG1KVe9IsR8NWN//0gk2anCF0uF4Dq5wE5pernAxG+OSrOYfsLj1Y3V/f8JWiaq9O2M8j+4X+wcM+/MLbrnfhg2Jyga65cLhfmzZtX++eOaAdlp22UH//0ooBneoj4+gjDYct/x6KUZsZd9yJJEhxehnfyZTycJqC9gdZe8cx5MBcX3nkKMdYTMNz7FKIGjlS7JJ/YXLIdj2x4CsesJ/DO0FcwKfMhzW7FUB9JkpCTkwNJktQuhTQSZadtlJ82PNpZROsw4K29/12LpTQzLqcIF5cbMGWHjANj9Eg1U4PFI9nthOWbz2H9eQ3CUjPR8v5s6Fu2VbusZjtpPY3P9i7GpuLNGJTYH8/0fQKmsOB6tA8hhBDlXs+V8HKujKNj9WgXKWh3itDmcOHNvTLu6ShQc8Upd2khTr/1F1hzvoX5zsdgfPQlvPreR5oe6i63HMcbO9/H/Wsfw87juzCt35OYNXBa0DdXLpcLM2fO1HR2VyvKTtsoP+14oquIUBF4N796FEtpZtzNeawpYSiuAlYMpYXtPGGMwXUwF9af18C5/1eExKcg5pkPERKTAIfDgfLycshy3Zuy8ay4ohSLCpbhx5JtaBFmxuRrH8btqcMRoQ9Xu7SAkGVZs9ld7Sg7baP8tKNFmIDHu4j4aL+M6Rmi4sy4myJM//Is4tuYsfEW7nq/q5LsdsL+2w+wbl0D78lShLTviKgbRiOyz1AIOm1mdNJ6GtvKf8G28l+Qe2ofYiKjMS79LtzaaSjCdKFql0cIIYQzJ+wMHZd48VJPEVOSbY2bIiwqKkL//v2RlpaGPn36oKCgoM43LFy4EKmpqejUqRMmTpx4yX4Q9b2mVMEFYHoGdzOXVxXJWgF73s+48PX7ODHzQVQs/xD6tvFoM+UNtH32Ixj6jbikueL9qfBeWcLB84fx+d6v8OcNf8V9aybin7lfIEQMwYysp/Cv2/+JO9NuvSqbK96zI1dG2Wkb5actsZECxqeKeDdfxgWrssxq/5acPHkyJk2ahPHjx2P58uUYP348fvvtt0sOLi4uxt///nfs3r0bMTExGD16ND7++GNMmTKl3tcao1drhkGxtPYqUJjHDe+5k/CcLIHr0F64Du2F92QJAEAfHQfDdcMRNXCUJp4lyBiD1WNDSWUZii4Uo+j8ERy6UIziihK4ZQ+iQgzoF9cLD6Tfg+va94QhJFLtkgkhhGjEsxkiPimUsaRYWY8iMMbY6dOnkZKSgvPnz0Ov14MxhtjYWGzbtg0pKSm1B8+dOxeHDx/G/PnzAQAbNmzA66+/jm3bttX7mhI1U4SL917AuO4tGnnbhDEGSB4wjwfM6wHzusFcDsh2K2SHFbKtqvq/9ipIF07De/YEpHMnIVWerT2Hvk17hKX0QFhqBsI6dYeuRZuA34ckS3DLHrglNzySB27JA7fshs3jgM1tg9Vjh81jh81tg8VVhTOOczhjP4ez9ur/OqXqxYd6UY8O5gSktOyI1JbJSG2ZjK5t0hCiCwn4PRFCCAkO4370YkdxJUontGlwilAPAGVlZYiNjYVeXz2gJQgCEhMTUVpaekmDVVpaiqSkpNrvO3TogNLS0gZfq4vL5bpkJX5lZSX6vHUjlv/6MFbuEoCalWEClH9d01Sq8XVT6m0UVseXrPpLxcvoBAiiAIg6CC11YK11EHVtAFEHJuog6vUAKwAryYdQJlY3bYxBEGu+BgRRAJNlMFQ/SkBmDIzJOH/+Alq2bAlBEAChuuFj//nhyEyuLVpiMhgYGGOQmASZMchMgsRkyP95TYkIXTiMoVFoHd4SJr0RvaIz0CLUDKMYiZQ2yYiNiIHk9iIqKgoejwdutxsOmwOV7kp4PB4YDAa4XC5IkoTIyEi4XC7IsoyIiIjaIfvw8HA4HA6IooiwsDDY7XbodDqEhYXBZrMhJCQEoaGhsNlsCA0NRUhICKxWK8LDw6HX61FVVYWIiAjo9XpYLBYYDAbodDpYLBZERUVBEARUVVXBaDRWj75ZrTCZTJAkCTabDSaTCV6vFw6HA0ajEV6vF06n85J7MhgMcLvdTb6n8+fP44UXXsC8efMgy3JQ3FMw5lTXPUmShL/+9a+YO3cuTCZTUNxTMOZ0pXvyeDx46qmn8Pbbb8NoNAbFPQVjTpff0yNtq7BkV/Xyp4aWsKu2Snn27NmYNWvWH379tzqOJYT41+LFi9UugTTRkiVL1C6BNMPSpUvVLoE00blz52A2m6/4uh4AEhIScOLECXi93topwtLSUiQmJl5ycGJiIg4fPlz7/dGjR2uPqe+1usyYMQPZ2dm131dUVCApKQmlpaX1Fkz4Y7FYkJCQgLKysnqHSwl/KDvtouy0jfLTrsrKSiQmJqJVq1b1HqcHgLZt26Jnz55YtGgRxo8fjxUrViA+Pv6S6UEAuPvuuzFw4EDMnDkTMTExmD9/PsaOHdvga3UJCwtDWFjYH37dbDbTbzaNMplMlJ1GUXbaRdlpG+WnXaJY/44Hta8uWLAACxYsQFpaGubMmYPPP/8cADBhwgSsXbsWAJCcnIxZs2ZhwIABSElJQXR0NCZPntzga4QQQgghVxPuNhptaFU+4Q9lp12UnXZRdtpG+WmX0ux0M2fOnBm4suqn0+kwePDg2k8zEu2g7LSLstMuyk7bKD/tUpIdNyNYhBBCCCHBgp5JQwghhBDiY9RgEUIIIYT4GDVYhBBCCCE+FvAGq6ioCP3790daWhr69OmDgoKCOo9buHAhUlNT0alTJ0ycOBEejyfAlZLLKcnuxx9/RN++fdG1a1ekp6dj2rRpkGVZhWrJxZT+uQOqH/9w4403okULeiYoD5Rmt2/fPgwePBhdunRBly5dsHLlygBXSi6nJDtZlpGdnY2uXbuiR48eGDJkCA4dOqRCteRiU6dORYcOHSAIAvLy8q54XL29CguwIUOGsM8//5wxxtiyZctY7969/3DMkSNHWGxsLDtx4gSTZZmNGjWKffDBBwGulFxOSXa7d+9mhw8fZowx5nA42IABA2rfQ9SjJLsab7/9NpswYQIzm80Bqo7UR0l2NpuNdezYkW3dupUxxpjX62WnT58OZJmkDkqyW7VqFevbty9zu92MMcZeeeUVNmbMmECWSeqwZcsWVlZWxpKSklhubm6dxzTUqwS0wTp16hQzGo3M4/EwxhiTZZnFxMSwoqKiS45788032eTJk2u/X79+PRswYEAgSyWXUZrd5aZMmcJeeumlAFRIrqQx2eXn57Prr7+eHTp0iBosDijN7pNPPmHjxo1To0RyBUqzW716NcvIyGAWi4XJssyeffZZ9vTTT6tRMqlDfQ1WQ71KQKcIy8rKEBsbW7tvhCAISExMRGlp6SXHlZaWIikpqfb7Dh06/OEYElhKs7vYyZMnsXz5cowcOTJQZZI6KM3O4/Fg4sSJWLBgAXQ6nRqlkssozW7//v0ICwvDyJEjkZmZiT/96U84c+aMGiWT/1Ca3ahRozB48GC0a9cOsbGx+OGHH/Dyyy+rUTJppIZ6FVrkTvzCYrFg1KhRmDZtGnr37q12OUSBWbNm4a677kKXLl3ULoU0ktfrxffff48FCxYgNzcXcXFxePzxx9Uuiyiwa9cu5Ofn49ixYzh+/DhuuukmPPbYY2qXRXwgoA1WQkICTpw4Aa/XC6B6MW1paSkSExMvOS4xMRElJSW13x89evQPx5DAUpodAFRVVWHEiBEYPXo0srOzA10quYzS7LZs2YL3338fHTp0wMCBA2GxWNChQwcaCVFRY/4/c8iQIYiLi4MgCHjwwQexc+dONUom/6E0uy+//LL2QyWiKOLhhx/G5s2b1SiZNFJDvUpAG6y2bduiZ8+eWLRoEQBgxYoViI+PR0pKyiXH3X333Vi7di1OnjwJxhjmz5+PsWPHBrJUchml2VmtVowYMQIjRozAiy++qEap5DJKs9u6dStKSkpw9OhRbNu2DSaTCUePHkV0dLQaZRMoz+7ee+/Fb7/9BovFAgDYsGEDMjIyAl4v+S+l2SUnJ+PHH3+E2+0GAHzzzTfo1q1bwOsljddgr+K/pWF1O3DgAOvXrx9LTU1lvXr1Ynv37mWMMfbnP/+ZrVmzpva4jz/+mCUnJ7Pk5GT26KOP1n7CgqhHSXavvvoq0+v1LCMjo/Z/r776qpplE6b8z12N4uJiWuTOCaXZffnllyw9PZ11796djRgxgpWWlqpVMvkPJdk5nU42YcIEds0117Du3buzm2++ufaT2EQ9kyZNYnFxcUyn07G2bduyTp06McYa16vQswgJIYQQQnyMFrkTQgghhPjY/wGhgF/GLW+I0gAAAABJRU5ErkJggg==\" />"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"plot(x, [ϕ(x,0.0,a),soln1(0.2),soln2(0.3)])"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 0.6.0",
"language": "julia",
"name": "julia-0.6"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "0.6.0"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment