Skip to content

Instantly share code, notes, and snippets.

@shivin9
Last active June 28, 2017 06:27
Show Gist options
  • Save shivin9/cd94aac140f9cf13183d88db334cc163 to your computer and use it in GitHub Desktop.
Save shivin9/cd94aac140f9cf13183d88db334cc163 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 25,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"using BenchmarkTools, DifferentialEquations, Plots\n",
"pyplot();"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"heateqsystem (generic function with 1 method)"
]
},
"execution_count": 32,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"function heateqsystem(t, u, du)\n",
" du[1] = -3.5*u[1] + 4*u[2] - 0.5*u[3] - 3*1/length(u)\n",
" @fastmath @inbounds for i in 2 : length(u)-1\n",
" du[i] = u[i+1] - 2*u[i] + u[i-1]\n",
" end\n",
" du[end] = -3.5*u[end] + 4*u[end-1] - 0.5*u[end-2] + 3*(-1)/length(u)\n",
" scale!(du, length(u)^2)\n",
"end"
]
},
{
"cell_type": "code",
"execution_count": 33,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"test (generic function with 1 method)"
]
},
"execution_count": 33,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"function test(N)\n",
" h = 1 / (N-1)\n",
" x = collect(0 : h : 1)\n",
" u0 = -(x - 0.5).^2 + 1/12\n",
" tspan = (0., 1.)\n",
" prob1 = ODEProblem(heateqsystem, u0, tspan)\n",
" sol1 = solve(prob1, dense=false, saveat=[0.03, 1.0])\n",
" return x, sol1, prob1\n",
"end"
]
},
{
"cell_type": "code",
"execution_count": 34,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"x, sol1, prob1 = test(100);"
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAYAAAByNR6YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xl8FPX9x/H3zGxiQi5QDoEgggiBiKCAEkRBUdGq4EG0tlqpIlqqFtLWqrTVVlG0ErXU1goWr1olItarIt6i4RJECLdcichN9s7uzsz390dKflUDTJZs5vudfT8fDx+PVEv4bl7Z7DfzmZ3RhBACRERERNRs9Mb+5YYNGzBkyBD07NkTgwYNQlVVVaN/eMuWLRg+fDgKCgrQv3//lC6UiIiISBWNbrBuuukmjB8/HuvXr8dvfvMbjB07ttE/nJ+fj/vuuw8vvPBCKtdIREREpBTtuyPCXbt2oUePHti3bx98Ph+EEOjYsSMWLFiAHj16NPpJPvzwQ0ycOBFffPFFk/5y27axfft25OXlQdO05B8FERERUQsQQiAYDKJTp07Q9UaPUwEAfN/9F9XV1ejYsSN8vvr/pGkajjvuOGzbtu2gGyynYrEYYrFYw//++uuv0adPnyP6nEREREQtrbq6GoWFhQf97wffeqXAAw88gIKCgoZ/DmyuSktL4ff7MWHCBEyYMAF+vx/XXnst7rjjDvj9fowePRpTpkyB3+/HOeecgz//+c/w+/0YOHAgnn76afj9fvTq1Qtz5syB3+9Hx44dMX/+fPj9fuTm5mLhwoXw+/0AgNWrV6O6uhpA/Rdn9erVAAC/34+FCxciNzcXfr8f8+fPR8eOHeH3+zFnzhz06tULfr8fTz/9NAYOHAi/348///nPOOecc+D3+zFlyhSMHj0afr8fd9xxB6699lrPP6adO3eid+/e+NWvfuWZx+TFTo09pnvuuQcTJkzA8OHDXXtMNbtr8ZsZr6Pz2D/hgjl7kfdYNfDYHuCxPdAe3YXc+7/ExXP34vQpb2LQL/+OF77cj9Ky+zD21tuxv7a2WTrt2leLVl2K8Pf5y/HU8v3A2RNw9ZxtOONfO4E/rmxYDx7ZiROf3oOL/rkJmWePx4LNtfjPO+597z355JOYMGECevbsqdz3nhefT019TG+88QYmTJjgqcfkxU4He0wff/wxACAvL+/Qmx4hhHjmmWdEv379RL9+/cTUqVNFXl6eSCQSQgghbNsWHTp0EBs2bBAH88EHH4h+/fod9L8fUFdXJ/x+f8M/1dXVAoDYsmWLEEKIaDQqotGoEEKISCQi6urqhBBChMPhho9DoZCIxWINH8fjcSGEEMFgsGHNgUCg4WO/3y9M02z42LIsYdu28Pv9wrZtYVmW8Pv9QgghTNNs+DiRSIhAINDwcTAYFEIIEY/HRSgUEkIIEYvFGj6uq6sT4XC44eNIJOL5xxSNRsWtt94qamtrPfOYvNipscdUW1srJk2aJPbs2dNij2nbnoCYs8kS4z+Kiz4v1QntybjAk3FRMCsmznszIe5cGBPPrg6LVXttEQjL0Wlv1BKfbDfFn5eFxIQFpjh1TlwYM2ICT8ZFq3/ExdC5deLupab4qCYu9vtb7nsvEAiISZMmid27dyv3vefF51NTH1MwGBSTJk0SO3fu9Mxj8mKngz2mb775RgBoWPvBfO8cLAAYPnw4xo4di7Fjx+Lll1/G1KlTsXTp0oNu0pI9BysQCKCgoAB+vx/5+flN+rNEJDfTFli0S+CdrwXm1wgs2i1gC6BnATCso4aS9joGt9fQqzWgK3QOZsQUWLZHYOEugc92CnywXaA2DuRnACM6azivs4aRhTq656vzmIjIOad7l0Y3WOvWrcPYsWOxd+9e5OfnY9asWejbty8AYNy4cRg1ahRGjRqFSCSCnj17IhaLwe/3o3379rj22mvxwAMPNGmRO3bsQIcOHZJ8qOSWaDSKW2+9FdOnT0d2drbby6EmSFW7hC3w4XaBlzcLzN1iY3cd0DoTOLezhvMLdZzXWcPxed7aeFi2wNI9Au/U1P9TuUvAEkC/o4Ex3XSM6a6jqHXzPmY+99Tm5X7xeBxbt26FZVluL+WIGIaBrl27IjMz83v/befOnTj22GMPu8H63knuANCrVy9UVlY2+gdmzpzZ8HGrVq1QU1PT1HV/z6HOwid56bqOwsJC9lNQc7YzbYF3vxao2GTj1a0C+2JAtzxgbE8dlx2v4bR2GgzdW5uq/2XoGk5vr+H09sDvTgUC8fqN1pwtNh780sbvPrdR3Aa4opuOH3bX0bvNkX8t+NxTm1f71dTU4Ec/+hEikYjbS2kWrVq1wr/+9S907tz5W//eabdGj2C1FI4IidS1rlZg1nobz26w8U2kfvQ3ppuOMd109D8GvPQKgDqzfkT68iYbr20T8MeB09truL6njqtO0FCQya8ReYNt27j55puxY8cO/PGPf0RWVpbbSzoidXV1+N3vfoeOHTviiSee+NamyunepdEjWC0tEolwg6WgSCSC6667Ds888wxatWrl9nKoCZJtF4wLzN5Uv7H6dKdAm6OAH5+g46e9dJzCTdX3ZPk0jOqqYVRXHTFL4I1tAv9YZ+Nnn1qYWAlc0U3DT3vqGN5Ja9J5aHzuqc2L/fbs2YNly5ZhypQpnrmzyy233ILJkydj7969aNeuXcO/d3qETooNlmEYbi+BkmAYBkpKSthPQU1ttzkgML3KxlPrbAQTwPmFGl48x8DorhqyfNxUOXGUoeGKbhqu6Kbj67DAcxts/GOdjec3WihqDfyiWMe1J+rIyTj815PPPbV5sV9tbS0AHPK6UKo58Fj279//rQ2W024cERJRo4QQ+HiHwKMr68dbrTOBm3rruLm3juNyualqDkIIfLJD4M9VNuZuESjIBMYX6fh5Hx1d+DUmhaxduxbXXHMNnn/+eRQVFbm9nGZxsMfkdO8ixRl24XDY7SVQEsLhMEaOHMl+CjpUO8sWePErG6fONTH8DQvr/QJPDDVQ/SMf7h9kcHPVjDRNw1kddbx8rg9fXeXD9T11PLHGRrcXTfzwPRNf7m38918+99TGfu7ZsmULhg8fjoKCgqRHmU67SbHBysjIcHsJlISMjAyUlpayn4Iaa2fa9WOr4pdNXP2+hQ7ZGt650MCqMT7cWKSjFUeBKXV8noaHBxuo+ZEPj5XoWLxboN8rJi57x8SyPd/eaPG5pzb2c09+fj7uu+8+vPDCC0l/DqfdpNhgNXadCZJfZmYmxo0bx34K+t92Cbv+xOuiChM/+dBCr9YaFl9q4O0LfTivUOeJ6y0sN0PDz4sNrLvSh6eHGajaLzBgromL3zaxaJcNgM891bFf6j388MMYP358w/+ura1F27ZtAQBDhw5FTk5O0p/baTcpTnIPh8M8B0tB4XAY5557Lt59990j+mallhcOhzHi3PMwfsZ7uHelD1tCwOXHa6gYYeCUttxQySBD13BdTw0/7qFh9iaB+5ZbGPxvCxd2sXF33zgmjhnB556i0uVn56ZA/V0OmlvrTBz2Tgnjxo1Dz5498dBDD6F169aYNWsWRo8ejaOPPvqI/36nI0IpNljcxaspMzMTZWVl7KegT/ZkYteN/8YNlT6M7qrhtZEG+h7NjZWMfLqGH/XQcFV3DS9vFvjtUgslb2Vg2M3PYWciE93dXiA1WTr87NxTJ3DibBN2Ct5GZ2jAjmt8aJt18J9ZrVu3xpgxY/CPf/wDkyZNwt/+9je89NJLzfL3K3UEi3NoNR04j4DUsWKvwG8WW5hXAwxu3wbPnq5j6LFSnClAh2HoGq46QcPl3TQ8ucbGH5Z1RfErwC9OsnBHPx2tj+IGWRXp8LOzbZaGDVf6UnYE61CbqwNuu+02jBo1Cr1790a7du1wyimnNMvf73TPIsUGKxQKcUSooFAohNNPPx2LFi1Cbm6u28uhQ9gVFbhzsYVZ6wVOLACeH5rAlDGD0H/RIgBsp5IMXcN1XaOYfs05uP7xjzC9yocZa23cN1DH+CLd07cl8op0+dnp9g3Pi4qK0L17d4wfPx4PPfRQs33eUCjk6P8nxa+uql9SP11lZWWhvLyc/SRm2QKPV1noNdvE3K0C04foWDXGh6tOzMQjbKesrKwsPDb1Xtx3WgY2XuXDpV01TPjUxmn/NrFwp+328ugw+LOz5dx4440wTRNjxowBUH8V9sLCQpSWlmL16tUoLCzEnXfe2aTP6bSbFEewfD4plkFN5PP5MHLkSLeXQQdRudPGzz+18MVeYFyRhvsHGf9/WF1nO5X973Ovow94apgPNxbZ+PlnFkpes3B9TxtTTzPQLptHs2TEn50t54MPPsCECRMaxnqtWrVCTU3NEX1Op3sWKY5gBYNBt5dASQgGgygsLGQ/yeypE7j+IxNDXrNgaBoWjjbw5JnfPiGU7dTWWL/BHXQsHu3DX8/QMXerQM/ZJv622oLt3s066CD4/Eu97du3o6ioCMuWLcPEiROb9XM77SbFoaPs7Gy3l0BJyM7ORkVFBftJQgiBik0Ct3xmwRTAE0N1jOvV+Dk5bKe2g/UzdA0/62NgTDcddy6xMOFTG//6SmDmmQZ6tubRLFnw+Zd6nTp1wtq1a1PyuZ12k+IIFkeEavL5fCgpKWE/CXwTEbh8voWr3rdwVkcNq8f4cFNv46AnPLOd2g7Xr122hpln+fDBRQa2R+qvCP+nFRbMVLxnnpqMzz+1KTUiDAQCbi+BkhAIBJCfn89+LhJC4On1NvpUmPhsl0DFCAMvn+vDsa0OfbSC7dTmtN/wTjq+vMKHCX103LHERsm/Lazcx02W27z4/NP1+u1EIpFweSXN58BjOfDYDnDaTYrts5evZOtlOTk5qKysZD+XfBMRuP4jC2/XCFzbQ8MjJQaOcXBtGIDtVNeUfq18GqYNNlDaTcP1H1sYMNfE3afquKMfL+ngFi8+/zp16oTMzEzMmDEDN954o/LXt0wkEpgxYwYyMzPRqVOnb/03p900Idw7AzIQCKCgoAB+v5/XwSJqgle32Bj3sYUMHZh5loGLjpPiYDRJLmYJ/HGZjakrbAxpr+G5sw0cn8dNFjWPhQsXoqysDPF4Cq4u6oLMzEyUl5dj8ODB3/r3TvcuUmywqqurUVhY6NYyKEncILe8cEJg0kILM9YKjO6qYeZZhqMrGn8X26ntSPst2GHjmg8s7I8Bfx1q4Mc9uEFvSV5+/oVCIWzfvh22rfb12HRdR6dOnRq9EGxNTQ26dOmixgZr//79aN26tVvLoCTZto3t27ejU6dO35tRU/NbstvGj9+38HUEeLTEwLheGjQtuaMPbKe25ujnjwvc8qmF5zcKXH2Chr+eYfB2Oy2Ezz+11dbWok2bNofdYElRNtkXCXKXpmnIz89nvxSzhcBDKywM+beFgkwNyy/z4cYi/Yi+7myntuboV5Cp4bmzfXjhbANvVde/07CSV4FvEXz+qc1pNyk2WLzYmpqCwSAKCgrYL4X2xwRGv2PhN4tt/OpkHZ+Nbp7rGbGd2pqz39U9dKy43IcuORrOet3CoystuDjYSAt8/qnNaTcpRoS1tbUoKChwaxmUJCEEgsEg8vLy+JtYCizdbaP0XQv+BPDc8OY9kZ3t1JaKfglb4M7FNqattHFFNw1PnWWgIJPfG6nA55/a/H4/WrdurcaIkL8tqUkIgUAgwH7NTAiBv622cMZrFtpl148Em/tdgmyntlT0y9A1PDzYwCvnGphfIzBwrokVe/n9kQp8/qnNaTcpNlihUMjtJVASQqEQunTpwn7NKJwQuOaD+lucjC/S8cklBrqm4G30bKe2VPa7rJuOZZf7kJcBDP63iX+s43lZzY3PP7U57SbFiNCLb1UlaqqtQYHR75jYGACeOsvAVSdI8fsPpak6U+C2yvpLgtzSR0d5iY4MXpiUyPHeRYqf4JZlub0ESoJlWaiqqmK/ZvDxNzYGvmoikAAWjvalfHPFdmpriX5ZPg1PnunDE0N1PLHGxsi3LOyt40irOfD5pzan3aTYYIXDYbeXQEkIh8MoKSlhvyP0xGoLI9600LeNhiWX+nDS0ak/SsB2amvJfjf1NvDeRQZW7hcY9KrJexk2Az7/1Oa0G0eERC6JWwK/qLTxxBqbIxiS3pagwKX/HWE/N9zAZd2k+P2cqMUpNSI0TdPtJVASTNNEZWUl+yVhf0xg5H8sPLXOxowzDUw/w2jRzRXbqc2Nfsfnafh0lA8XdtFw+bsW7l/O62Uli88/tTntJsUGKxqNur0ESkI0GkVpaSn7NdHmgMCQ1+pHLe9fZGBcUcs/DdlObW71y8nQMHuEgXtO1TF5qY3xn1hI2NxkNRWff2pz2o0jQqIWtGS3jYvnWcjLAP5zgQ8nFnAkSGp6Zr2NcR9bGNFZQ8UIA3m8KCmlCY4IKeVM08S8efPYz6HXttoY/oaF7nkaKke5u7liO7XJ0O+6njrevtBA5U6BM1838XWYR7KckqEfJU+pEWFdXZ3bS6Ak1NXVoaysjP0ceLzKwmXzLVzYRcP7Fxlol+3ub/tspzZZ+o3orOPTUT7si9VflJTvMHRGln6UHKfdOCIkSiEhBO5cYuPBFTbK+ur40+k6dN57jDxme1jg4nkmvgoA/z7fwPBOUvzuTpQSSo0IE4mE20ugJCQSCVRUVLDfQVi2wPhPLDy4wsYjg3VMG2xIs7liO7XJ1q9TjoaPLvbh9PYaLnjbwmtbeXudQ5GtHzWN025SbLDi8bjbS6AkxONxlJeXs18jYpbAVe9bmLVe4JlhBib2Ndxe0rewndpk7JeXqeH1kQYuOU7D5fMtPLuem6yDkbEfOee0G0eERM0slBC4bL6FT3YIzB5hYFRXKX6PIWoRli3ws0/r72FYPljHJMl+uSA6UkqNCLmLV1M8HsfMmTPZ73/srRMY8aaFRbsE5l0o7+aK7dQmcz9D1/D3oQbu6KejbKGN3y7hBUm/S+Z+dHhOu0nx059zaDXxPIJv+yYicNbrJjYHBT682IdhHaV4ejWK7dQmez9N0/DAaQYeOk3HlC9s3PKZDZubrAay96NDc9qNI0KiZlATEjjnTRNRC3j3Bz70ai3HyexEbpu5tv6K7zf00vD3M+V5owdRspQaEcZiMbeXQEmIxWIoLy9P+37bQgLD3jARs4GPLlZjc8V2alOp37giHU8PM/CP9QLXf2TB4q11lOpH3+e0mxQbLMuy3F4CJcGyLFRWVqZ1vy1BgWGvmxCo31x1z5d/cwWwnepU6/eTnjqeG27guY0C131kwUzzTZZq/ejbnHbjiJAoSV8FBM5+w0SmAXxwkQ9dctXYXBG5pWKTjavftzCmm4bnzjaQofM5Q+rhiJBSLhaL4Z577knLfutr68eC2b76I1eqba7SuZ0XqNqvtLuO2SMMzNkscPX7FuJWeh7JUrUf1VNqRGjbvCCdimzbRk1NTdr1+yogcPabJvIzgA8v9qFzjlqbKyB923mFyv0u76ZjznkGXttav8lKx3Ghyv3I+Z6FI0KiJtgaFDjrDRNZRv2Rq2Nbqbe5IpLB61ttXD7fwpXdNTw73IDBcSEpQqkRIe8orqZ0uyP81+H6SzH4NOD9i9TeXKVbO6/xQr9Luur41zkGXtwkcOMnVlpdJ8sL/dKZ026+FK+DyBN2RgRGvGkiYQOfXKLmWJBINmO663jWBq79wEKWYePxM3RovE4WeQRHhESHsaeu/t2Ce2PAxxf70KOALwBEzemptTbGfWJh0kk6pg3mJovkptSIMBqNur0ESkI0GsW4ceM83a82JnD+WyZ2RoH3fuCdzVU6tPMyr/W7oUjHX4boeGSVjclLvX/it9f6pRun3aQYEeq6FPs8aiJd11FYWOjZfhFT4OJ5FraEgA8v8qF3G29srgDvt/M6L/b7ebGBOgv41SIbrTOB2/sZbi8pZbzYL5047cYRIVEjErbApe9Y+PAbgfd+YGBwB/4gJGoJv11iYcoXNp46y8D1vfi8I/koNSKMRCJuL4GSEIlEUFpa6rl+tqi/Z9r8rwVeOdebmyuvtksXXu5370AdNxXpuPETC69u8ea40Mv90oHTblK8chiGdw8Fe5lhGCgpKfFUPyEEfrnQxj83Cjw73MDILlI8RZqdF9ulEy/30zQNj5+h4/LjNfzwfQsfbvfeJsvL/dKB024cERL9j/uXW5i81MZfhuj4eTF/+BG5JWYJXPS2hcW7BT662IdT2nrnHEhSm1IjwnA47PYSKAnhcBgjR470TL8n19S/g+nuU72/ufJau3STDv2OMjTMPc9ArwINF7xtYoPfOxciTYd+Xua0mxQbrIyMDLeXQEnIyMhAaWmpJ/q9ttXGzz618PM+Ou4+VYqnRUp5qV06Spd+eZka3rrAQJtMYOR/TOyKemOTlS79vMppN44IKe0t2W1j2OsWLuyiYfYI3hONSDZbggKD/22ia66GDy420MrH5yi5hyNCSrlwOIySkhKl+20K1F/rqt8xGp4/O302V15ol87Srd/xeRreHOnDqv0CP3rfgmWrfSQr3fp5jVIjwszMTLeXQEnIzMxEWVmZsv321gn84G0T+RnAa+cbyE6j34pVb5fu0rHfgHb1R5hf3yYwsdKGi8OXI5aO/bzEaTeOCCkt1ZkC571lYa1foHKUd26BQ+R1f19j4eYFNqadrqPsZG+/GYXkpNSIMBQKub0ESkIoFEJxcbFy/WwhcN1HFpbuEXjtfCMtN1eqtqN66dzvpt4G7uin45eLbFRsUvMaWenczwucdpPiXoRZWVluL4GSkJWVhfLycuX63bnYRsUmgZfPNVDiwau0O6FqO6qX7v2mDNKxNSRw7YcWOucAQxR7Hqd7P9U57cYRIaWVp9fb+OlHFsoH65jUl+MFIlXFrP+O+WsFllzqQ9e89DsSTe5QakQYDAbdXgIlIRgMorCwUJl+C3bYGP+JhRuLNEw8SYpvfdeo1o6+jf3qL0T6ynkGcjOAS94xEYyrc9I7+6nNaTcpXmWys7PdXgIlITs7GxUVFUr02xwQuGy+hTM6aPjLEAOalt6/7arUjr6P/eq1zdLw+vk+bAkC13yozuUb2E9tTrtxREieF4gLDHnNRJ0FLBrtwzFZ6b25IvKat7bZuOQdC7/qq+PB0zn6p9RSakQYCATcXgIlIRAIID8/X+p+li1w9fsWqkPA6+dzc3WACu3o4Njv235wnI6HT9fx0Jc2nl4v/zsL2U9tTrtJ8S7CnJwct5dAScjJyUFlZaXU/W5fbOPtGoG3Rhro3YabqwNUaEcHx37fN/EkHav3C4z/xEKPfGDosVIcP2gU+6nNaTeOCMmzZq2zcf3HFh4r0XHbSRwbEHld3BI4/z8WqvYLLOU7CylFOCKklAsEAtA0Tcp+i3fZuHmBhRt6abi1WIpvc6nI3I4Oj/0al2loePlcAzk+4PJ3TURNOU96Zz+1Oe0mxStPbm6u20ugJOTm5qK6ulq6fjsjApe/a+GUthoeP4PvGGyMrO3IGfY7uLZZGuae58Oa/cBNn1hS3rOQ/dTmtJsUGyy+AKpJ0zTk5+dL1S9hC5S+Z8GygTnnGjjKkGdtMpGxHTnHfod2SlsNM88y8NxGgelV8p30zn5qc9pNig0WL7ampmAwiIKCAqn6/XKhjcqd9bfB6ZzDH14HI2M7co79Du9HPXSU9dVRttDGR9/ItcliP7UpdaHRvLw8t5dAScjLy4Pf75em39PrbUyvsvHnITrOkPgdRDKQrR01Dfs58+BpOoZ11FD6roVtIXlGheynNqfdpHgVknFGTocnhEAgEJCi39Ld/39S+829pfi2lppM7ajp2M8Zn67hpREGWvmAy+db0pz0zn5qc9pNileiUCjk9hIoCaFQCF26dHG93546gcvnW+h3NG+D45Qs7Sg57Odc2ywNr5znQ9V+gZ9/arm9HADspzqn3aS40CivgaWm/Px8138Ds2yBH79vIWrVn9Se5ePmygkZ2lHy2K9pTm2r4e9DDVz3kYWhx9q4vpe7xxbYT21O9yxSHMGyLDl+q6CmsSwLVVVVrvab8oWN+V8L/OscA4W53Fw5JUM7Sh77Nd1PeuoYX6Tj559a+GKvy78Ysp/SnHaTYoMVDofdXgIlIRwOo6SkxLV+82ts3PO5jT8M0HFuZym+lZXhdjs6MuyXnMdKdPRuDYx514Q/7t4mi/3U5rQbb5VDSqoJCZwy18SAthreusCAzvOuiMiBTQGBU+eaOKeThjnn8pxNajqlbpVjmqbbS6AkmKaJysrKFu+XsAWufM9CtgE8fzY3V8lwqx01D/ZLXvd8Dc8ONzB3i8AjK925Phb7qc1pNyk2WNFo1O0lUBKi0ShKS0tbvN9vFtlYukdg9ggDbbO4uUqGW+2oebDfkRnVVcftJ+u4fbGNBTtafpPFfmpz2o0jQlLKnM02xrxr4c8lOm49yXB7OUSkKNMWGPGmhY0BgeWX+9A+m7+skTMcEVLKmaaJefPmtVi/zQGB6z+yUNpNwy3FUnzrKqul21HzYr8j59M1vDjCgCmAn3xowW7BYw3spzalRoR1dXVuL4GSUFdXh7Kyshbpl7AFfvi+hWOygBln8cTUI9WS7aj5sV/z6NhKw7PDDMyrEZj2ZcuNCtlPbU67cURISvjNIgvlK218OsrAae2l+L2AiDzijsUWpn1pY8EoA6fz5wsdhlIjwkQi4fYSKAmJRAIVFRUp7zev2sZDX9p44DSdm6tm0lLtKDXYr3ndO1DHgHYafviehdpY6o85sJ/anHaT4tUqHo+7vQRKQjweR3l5eUr7fRMRuPZDCxcUaijrK8W3qye0RDtKHfZrXhm6hn+dbWB/HBj/iZXy29iwn9qcduOIkKRlC4Hz37JQtV9gxRV8lw8RpVbFJhtXvmfhyTMN3FjEX+iocUqNCLmLV1M8HsfMmTNT1m/qFzbe3y7w/NkGN1fNLNXtKLXYLzVKu+u4qUjHbZ9ZWLUvdcce2E9tTrtJscHiHFpNqTyPoHKnjd9/bmPyKTpG8D6DzY7ngKiN/VLnkRIdPfKBq94zETVTs8liP7U57cYRIUknEBfo/4qJY7M1fHyJAZ/Oo1dE1HJW7xcYMNfEuF46pp/BCxrTtyk1IozwRlgVAAAgAElEQVTFYm4vgZIQi8VQXl7e7P1u+8zCnrr6+wxyc5UaqWpHLYP9UqtPGw1/Ol3HX1bb+E91818fi/3U5rSbFBssy7LcXgIlwbIsVFZWNmu/ik02ntkgMH2Ige753FylSiraUcthv9T7eR8dFxRq+OlHFnZFm3fQw35qc9qNI0KSRk1I4ORXTIzopGH2CF6tnYjctSMi0HeOiZL2Gv59Pn8mUT2OCCnlYrEY7rnnnmbpZwuB6z6y0MoH/H0of5ClWnO2o5bHfi3j2FYanjrLwOvbBJ5c23yjQvZTm1IjQttuuXtAUfOxbRs1NTXN0q98Zf0lGZ4ZZuDoLG6uUq0521HLY7+WM6qrjvFFOiZV2lhX2zwDH/ZTm9NuHBGS677YK3DaqyZuK9bx8GC+Y4eI5BJOCJw610RehobPRhnINPhLYDpTakTIO4qrqTnuCF9nCvz4fRN9WgNTBknx7ZgWmqMduYf9WlZOhoZ/nm1gxV6BPyw78qNO7Kc2p934ikau+t3nNjYGgOfP9uEo/lZIRJIa2E7H3afqmLrCxqJdHO3R4XFESK5ZsMPGWa9bePA0Hb/ux9EgEcnNtAWGvGbBHxdYfrkPrXz8pTAdKTUijEajbi+BkhCNRjFu3Lik+oUSAtd9aGFIBw1lfaX4NkwrR9KO3Md+7vDpGp4dbmBbCLhrSfJHsdhPbU67SfHKputSLIOaSNd1FBYWJtXv9kU2dkSBp4cZMHi19hZ3JO3IfeznnqLWGu4fpOOxVTY+2J7cJov91Oa0G0eE1OLeqbEx8j8WHj9Dx4Q+HA0SkVpsIXDOmxa2BAVWXuFDXiZ/SUwnSo0II5GI20ugJEQiEZSWljapX21M4PqPLZzbWcPNvaX49ktLybQjebCfu3RNw6yzDOypA8oWNv12N+ynNqfdpHiFMwwexVCRYRgoKSlpUr9fVFoIxoF/nGVA59XaXZNMO5IH+7mvW76G8sEGZq4TeGtb00aF7Kc2p904IqQW8+oWG5fNt/D0MAPX9ZRib09ElDQhBH7wtoUV+wRWXeHjXSjShFIjwnA47PYSKAnhcBgjR4501G9fncDNCyyM6qrhJyfyh5DbmtKO5MN+ctA0DTPPMhAxgUlNGBWyn9qcdpNig5WRkeH2EigJGRkZKC0tddRv0kILMQv42xm8kbMMmtKO5MN+8uico+GRwQae3eB8VMh+anPajSNCSrm3ttm4aJ6FWcMMjOVokIg8RgiBC/5jYXWtwKoxPhTwXYWexhEhpVw4HEZJSckh+/njAjctsDCyUMN1HA1Kw0k7khf7yUXTNDx5poHaeP01/g6H/dSm1IgwMzPT7SVQEjIzM1FWVnbIfrcvslEbB548k6NBmThpR/JiP/l0zdPw0Gk6nlxr472vD73JYj+1Oe3GESGlzHtf2zj3LQt/O0PHzbygKBF53IELkG4LCXx5hQ+5Gfyl0ouUGhGGQiG3l0BJCIVCKC4ubrRfKCFw4ycWhnfUMJ4XFJXOodqR/NhPTrqmYeaZBnZEgMmHuFch+6nNaTcpXvmysrLcXgIlISsrC+Xl5Y32m7zExo4IMJMXFJXSodqR/NhPXj0K6u9VOL3KxoIdjW+y2E9tTrtxREjN7tMdNs583UL5YB0T+3I0SETpxbIFznzdwt6YwIrLfcjy8ZdML1FqRBgMBt1eAiUhGAyisLDwW/1iVv1ocFA7DbcWS/HtRY1orB2pg/3kZuganjrLwOYgMOWL7x/FYj+1Oe0mxStgdna220ugJGRnZ6OiouJb/R5aYWODH5hxpgFD529tsmqsHamD/eTXu42Gu/rrmPqFjVX7vj0oYj+1Oe3GESE1m7W1Av3mmPjlyTruH8TRIBGlt5gl0P8VE20yNSwYxfNRvUKpEWEgEHB7CZSEQCCA/Px8BAIB2EJg/CcWjssFfneKFN9WdAj/247Uw35qOMrQ8ORQA5W7BJ5Y8/+jQvZTm9NuUrwS5uTkuL0ESkJOTg4qKyuRk5ODf6wT+GSHwN+HGsjmCZ3S+992pB72U8eZHXWML9Jxx2IbX4frB0bspzan3TgipCO2IyLQu8LEpcdrmDXM5/ZyiIikUhur/xk5uL2GuefzZ6TqOCKklAsEAtA0DRM+rkOGDjx8Os+7UsWBdnzuqYn91NL6KA3Thxh4davA3M02+ylOqRFhbm6u20ugJOTm5mLW5zsxt9qHR0sMHJPF0aAqcnNzUV1dzeeeothPPVd003DJcRpu+cyClZnDfgpz2k2KDRZvAqymsAn8fk0BRhZquPoENlSJpmnIz8/nc09R7KceTdPw+BkGAglg8lLBfgpz2k2KDRYvtqam3y6qQ3VtDA/2i/IHhWKCwSAKCgr43FMU+6mpS66GewfoeGKNjYKTz2E/RSl1odG8vDy3l0BNtGqfwF/XZeC3A3w4uSMPc6smLy8Pfr+fzz1FsZ+6binWcfLRGvo//Cla5fBnp4qcPu+k2GC5+EZGSoIQAhM+tXBCPjC24z72U5AQAoFAgO0UxX7q8ukaHh+i44v9+reujUXqcPq8k2KDFQqF3F4CNcEzG+qvefWnU2LocXwX9lNQKBRCly5spyr2U1vfnDDw8UxMXmpjR4SbZNU4fd7xOljUJHvrBIoqTJzfWcM/z+H1XIiIkrGvTqAXf5YqSanrYFmW5fYSyKE7l1hI2MC0wQYsy0JVVRX7KYjt1MZ+arMsC998tRoPDtLwwlcC733NUaFKnD7vpNhghcNht5dADlTutDFjrcCUgTqObaUhHA6jpKSE/RTEdmpjP7Ud6HdFxwjOPFbDhE8txCyOClXh9HnHESE5YtoCA+eayNA1LBxtwNB5WQYioiNVtU+g/ysm7j5Vx29P5d0wVKDUiNA0TbeXQIfxlyobX+4D/jZUb9hcmaaJyspK9lMQ26mN/dT2v/2Kj9ZQ1lfHlC9sbArwKJYKnD7vpNhgRaNRt5dAh7AjIvD7z23c1FvHwHb//y0TjUZRWlrKfgpiO7Wxn9q+2+/3p+pomwVMWshz6lTg9HnHESEd1tgPTbyxTWD9lT4czfsNEhE1u9lf2bjqfQtvXWDgwi5SHPugg+CIkJrFZzttPLNB4P5Bxvc2V6ZpYt68eeynILZTG/uprbF+pd01nN1Rwy8+4wnvslNqRFhXV+f2EqgRli1wy6cWBrTVcEOv7x+5qqurQ1lZGfspiO3Uxn5qa6yfpmmYPsTApiDwyEpetkFmTp93HBHSQT2x2sLPPrVROcrA4A5S7MWJiDytrNLC39faWFfqQ2EuT8mQkVIjwkQi4fYS6Dv21glMXmrjpz21g26uEokEKioq2E9BbKc29lPbofrdPUBHXgbwq0U84V1WTp93jb5ybtiwAUOGDEHPnj0xaNAgVFVVNfqH33//fZx22mno06cPiouLcfvtt8O2m35oMx6PN/nPUGpNXmLDEsADgw5+XZZ4PI7y8nL2UxDbqY391HaofgWZGh46zcBLmwQ+2M5RoYycPu8aHRGec845+MlPfoKxY8fi5ZdfxoMPPoglS5Z87w8vX74cBQUF6N69O+rq6nDuuedi3LhxGDt2rKO/nCNCOS3bU39R0UdLdNx2Ei98R0TUkmwhcObrFgJxgWWX+5DBCztLJekR4a5du7B06VJcc801AIArrrgC1dXV2Lhx4/f+8CmnnILu3bsDALKystC/f39s2bKlyYvlb2HysEX9ie3FbYAJfQ49QY7H45g5cyb7KYjt1MZ+ajtcP13T8JchBqr2A39dzaNYsnH6vPveK2h1dTU6duwIn6/+7t6apuG4447Dtm3bDvmJduzYgZdffhkXX3zxQf8/sVgMgUDgW/8AQDAYBFB/Zv6Bs/Oj0ShisRgAIBKJNHwcDocbHlw4HG6YhYZCoYa3TgaDwYaPA4FAw40ZA4EAbNuGEAKBQABCCNi23bAOy7IaPjZNs2FdpmkiFAoBqJ+9HrgPUTweb/g4FoshEok0fHzgQmSqPaZn1iZQuUtg+hADViJ+yMeUSCTw4osvNnweWR+TFzsd6WMKhUKoqKhAbW2tZx6TFzsd7DFFIhFUVFRg//79nnlMXux0sMcUjUZRUVGBvXv3HvQx9T8GGF+k4fdLLeyOCukfkxc7HewxHfg8hyW+Y+nSpaJnz57f+neDBg0S77333nf/rw38fr8YOHCgmDZt2kH/P0IIcffddwsA3/vn2muvFUIIMWnSJDFp0iQhhBA33HCDuPvuu4UQQowZM6bhc59//vlixowZQgghBg8eLGbPni2EEKJPnz7i7bffFkII0blzZ/HZZ58JIYTIy8sTq1atEv8dhYrq6mrh9/sFAOH3+0V1dbU48GVYtWqVyMvLE0II8dlnn4nOnTsLIYR4++23RZ8+fYQQQsyePVsMHjxYCCHEjBkzxPnnny+EEGLatGlizJgxDY/zhhtuUO4xBeO2aP1krejwu/c885i82ImPiY+Jjyk9HtOKr2oEHtkpbvrY9Mxj8kKn1atXNzyeQ4EQQjzzzDOiX79+ol+/fmLq1KkiLy9PJBIJIYQQtm2LDh06iA0bNjT6CQKBgCgpKRH33nvvIf8iIYSoq6sTfr+/4Z8DX+itW7cKIYSIRqMiGo0KIYSIRCKirq5OCCFEOBxu+DgUColYLNbwcTweF0IIEQwGG9YcCAQaPvb7/cI0zYaPLcsStm0Lv98vbNsWlmU1fJFM02z4OJFIiEAg0PBxMBgUQggRj8dFKBQSQggRi8UaPq6rqxPhcLjh40gkotxjmrzYFEfNjImqHc4eU11dnZg6dWrD55fxMXmxU3M8Jr/fL6ZNmyb27t3rmcfkxU4He0zBYFBMmzZN7NmzxzOPyYudDvaYQqGQmDZtmti1a9dhH9PUJSGhz4iLz3cmpH5MXux0sMe0Y8cO5xus7xo2bJiYNWuWEEKIiooKMWDAgEb/cDAYFEOGDBF/+MMfDvmXHMyBXe0333yT1J+n5rM5YIujnoqLyYtNx38mHA6LMWPGNHyTkzrYTm3sp7am9ItbtiiaHRfDX08I27ZbYHV0ON98842jDVaj7yJct24dxo4di7179yI/Px+zZs1C3759AQDjxo3DqFGjMGrUKEyZMgX33HMPiouLG/5saWkpJk+e7Gg8yXcRyuPKd00s2Fl/v8HcDL5jhYhIFv+ptvGDty3MOdfA5d2kuHxlWnO6d5HiSu67du1Cu3bt3FpG2vv4GxvD3rDwzDADP+np/Mkbi8XwwAMP4M4778RRRx2VwhVSc2M7tbGf2pLp94O3TaytFVg9xocsH38JdtPu3bvRvn17Na7knszFSal5WLbALyotnNZOwzUnNu1Ja9s2ampq2E9BbKc29lNbMv3KBxuoDgGPrmJztzntJsURLI4I3TNzrY0bP7Hw2SgDJbzfIBGRtCZVWpi5zsb6K33o2IpHsdyi1L0IeUd4dwTiApOXWvhxDy2pzVVjd4QnNbCd2thPbcn2+/2pOo7SgclLeJ9CNzntJsUGi9xx33IboQQw9RD3GyQiIjm0OUrDvQN1PL1eYOlujgplxxFhmtoUEOhdYeKu/jruHsANFhGRCkxb4JRXTLQ5SsNHFxvQNI4KW5pSI8IDl7inlnPnEgtts4BfnZz8t0A0GsW4cePYT0Fspzb2U9uR9PPpGh4+3cAnOwT+vdW14yNpzWk3KTZYui7FMtJG5U4bszcJTBlkIOcIrnml6zoKCwvZT0Fspzb2U9uR9hvZRcfIQg23L7IQt7jJamlOu3FEmGaEEBjymoWoKfD5ZT4YOg8vExGpZuU+gf6vmHhksI7bTuJpHi1JqRHhgTtiU+q9vFlg4S6BaYONI95cRSIRlJaWsp+C2E5t7Ke25ujX92gN1/fU8IdlNvbHeBSrJTntJsUGyzC4+24JMUvgN4stXNRFw4jOR57eMAyUlJSwn4LYTm3sp7bm6nfvQAMxC5iynO8obElOu3FEmEamfWnhN4ttrLzCh95tOBokIlLdvcss3LfcxppSH7rn8+d6S1BqRBgOh91eguftrRO4b7mNG4v0ZttchcNhjBw5kv0UxHZqYz+1NWe/sr462mbVvzOcWobTblJssDIyMtxegufdu9yGJYA/DGi+5BkZGSgtLWU/BbGd2thPbc3ZLydDw5RBBmZvEqjcyVFhS3DajSPCNLDBL9CnwsQfB+q4sz/P2SAi8hJbCAyYayLL0PDZKF58NNU4IqQGdy2x0LEVMPGk5s0dDodRUlLCfgpiO7Wxn9qau5+uaZh2uoGFuwRe2cJ3FKaaUiPCzMxMt5fgWYt22Xh5s8C9Aw1k+5r3t5rMzEyUlZWxn4LYTm3sp7ZU9Duns44LCjXcudhCwuYmK5WcduOI0MOEEBj+hoX9cYHlvKgoEZGnrdhbf5/Cv56h4+Y+PB0kVZQaEYZCIbeX4ElvVQt8vENg6qAjv6hoY0KhEIqLi9lPQWynNvZTW6r69TtGwzU9NNyzzEYowaNYqeK0mxQbrKysLLeX4DmWLXDHYgvDOmq4sEtqjlxlZWWhvLyc/RTEdmpjP7Wlst8fBxrYHwMeWcl3FKaK024cEXrU0+tt/PQjC4tGGzitvRT7aCIiagG/XGjhybU2vrrKh/bZPDWkuSk1IgwGg24vwVOipsDvlloY001L6eYqGAyisLCQ/RTEdmpjP7Wlut9d/XUYGnDvMh7FSgWn3aTYYGVnZ7u9BE/5S5WNbyLA/YNSe5JjdnY2Kioq2E9BbKc29lNbqvsdk6Xhjn46nlhj46sAz8Vqbk67cUToMfvqBE54ycTVJ+j461C+i4SIKB1FTIGes00M7aDhxRE+t5fjKUqNCAOBgNtL8IypK2wkbOD3p6Y+bSAQQH5+PvspiO3Uxn5qa4l+rXwa/nCqgZc2CSzdzVFhc3LaTYoNVk5OjttL8ITqkMCfq2z88mQdx7ZK/YmNOTk5qKysZD8FsZ3a2E9tLdXvup4aercG7ljMDVZzctqNI0IPufFjE69uFfjqKh/yM/nOESKidPfqFhuXzbfw7g8MjOgsxTEV5XFEmGbW1wrMWi9wV3+9xTZXgUAAmqaxn4LYTm3sp7aW7De6q4bT2mm4a4kNF4+neIpSI8Lc3Fy3l6C8339ef0Pnn/VuuaS5ubmorq5mPwWxndrYT20t2U/TNNw/SMfi3QL/3soNVnNw2k2KDZamcZx1JJbvEXhpk8DdpxrIauYbOh+KpmnIz89nPwWxndrYT20t3W9EZx0jOmmYvMSCxRtBHzGn3aTYYPFieUfmt0stnJgPjO3Zsj9sg8EgCgoK2E9BbKc29lObG/3uH6RjdS3wwlfcYB0ppS40mpeX5/YSlLVgh423qgXuHWjAl4IbOh9KXl4e/H4/+ymI7dTGfmpzo99p7XVc2lXD3Z9biFvcZB0Jp92k2GDxxLvkCCFw52Ib/Y8BSru3/KhACIFAIMB+CmI7tbGf2tzqd99AA1uCwIy1vGzDkXDaTYoNVigUcnsJSnq7RmDBToEpAw3oLpyLEQqF0KVLF/ZTENupjf3U5la/4qM1XNNDw73LbYQT3Jwny2k3KTZYvAZW09lC4K4lFoZ20HBhF3dOdM3Pz4cQgv0UxHZqYz+1udnvDwMM7IsB06t4FCtZTrtJscGyLMvtJSjn5U0CX+ytP3HRrXcSWZaFqqoq9lMQ26mN/dTmZr9u+RrGF+l4cIWN2hiPYiXDaTcpNljhcNjtJSjFtAV+97mFC7toOLOjewnD4TBKSkrYT0Fspzb2U5vb/SafoiNmAQ9/yaNYyXDaTYpbbPMwd9P8c6PAej/wwtmGq+vgzWbVxXZqYz+1ud2vYysNtxTreKzKxsS+Otpm8XpqTaHUiNA0TbeXoIyELfDHZRYu7aphQDt3nxSmaaKyspL9FMR2amM/tcnQ7/Z+9S//D63gUaymctpNig1WNBp1ewnKeHq9wKZg/YmKbotGoygtLWU/BbGd2thPbTL0a5ul4RfFOv5SZWNHhOdiNYXTbppw8UIqTu9ITfVilsCJL5ko6aDhpRFSTHeJiEhR+2MCx//LxE976Xi0xP1f2lXhdO8ixREsHuZ25ql1NmrCwN2nyvFEME0T8+bNYz8FsZ3a2E9tsvRrc5SGsr46nlhj4+swj2I5pdSIsK6uzu0lSC9qCkxZbuNHPTT0aSPHCYl1dXUoKytjPwWxndrYT20y9ZvYV0crH3D/FzwXyymn3TgiVMSjKy38apGNNaU+nFggxwaLiIjU98AXFu7+3MaGK33omsfXl8NRakSYSCTcXoLUwgmBB1bY+MmJmlSbq0QigYqKCvZTENupjf3UJlu/W4t1tM4E7lvOC9c64bSbFBuseDzu9hKk9vhqG/vqgN+dIse5VwfE43GUl5ezn4LYTm3spzbZ+uVmaPhNPx2z1gt8FeC5WIfjtBtHhJILxgW6vWhiTDcdT5wp1waLiIi8IWIKnPCiifMLNTwznO9SPxSlRoSy7OJlNL3KRjBRf2sD2cTjccycOZP9FMR2amM/tcnYr5VPw139dTy/UWB9LY9iHYrTblK8assyh5ZNMC4wbaWNcUU6uuTKc+7VAbKdR0DOsZ3a2E9tsva7sUjHsdk8F+twnHbjiFBiD3xh4Z7PbWy8yiflBouIiLxl+ioLExfaWMt3rB+UUiPCWCzm9hKkE4wLPPyljRt6yXn0CqjvVl5ezn4KYju1sZ/aZO53Y5GODjyKdUhOu0mxwbIshvyux1fXn3t1Z38pEjXKsixUVlayn4LYTm3spzaZ+2X5NNzRT8c/Nwps9PNcrMY47cYRoYRCifr7Q5V21/G3oXznIBERtZyoKdD9RRMXdNEwaxjfUfhdHBEq7PEqGwHJj14B9d3uuece9lMQ26mN/dQme79sn4Y7+ut4bgOvi9UYpUaEts17IB0QSgg8vNLG9T11HCfpuVcH2LaNmpoa9lMQ26mN/dSmQr/xRTraZQFTeC7W9zjtxhGhZP60wsLkpbwnFBERuevAPXDXX+lD93y+Hh2g1IhQhjuKyyCcEPjTlzZ+2lNXYnMl0x3hqWnYTm3spzZV+t3UW0dbHsX6HqfdpNhgUb2/rbGxPwbcJfm5V0RE5H3Zvvp7FD67QWAzz8VqMo4IJREx6+85OLqrhifP5Ls2iIjIfZH/vqPwkuM0zDiLr02AYiPCaDTq9hJcN2Otjb11wJ391LksQzQaxbhx49hPQWynNvZTm0r9Wvk0/PpkHc9sENgW4lEswPmeRYoNlq5LsQzXxCyBh1bYuPZEDd0UOpFQ13UUFhamfT8VsZ3a2E9tqvW7qbeO/AzgTyvkfddjS3LajSNCCfx9jYWfLai/91PP1upssIiIKD1MWW7h3uU2tvzQh2NbpffrlFIjwkgk4vYSXJOwBaZ+YeOqEzTlNleRSASlpaVp3U9VbKc29lObiv1uKdaRZQDTVvIoltNuUmywDEOd846a2z83CmwJAXf1V+9rYBgGSkpK0rqfqthObeynNhX7FWRquLVYx99W29hTl97nYjntxhGhiyxboHeFieI2Guaez3dnEBGRvPbU1d8nd+JJOu4bpM7msLkpNSIMh8NuL8EVFZsFNgSAyadIkaHJwuEwRo4cmbb9VMZ2amM/tanar22Whp/10TG9ykZtLH2PYjntJsUre0ZGhttLaHG2EJiy3MIFhRoGtpMiQ5NlZGSgtLQ0Lfupju3Uxn5qU7nfL/vqiNnAX6rS91wsp904InTJq1tsXDbfwoJLDJxxrJobLCIiSj+3fmrhha9sbL3ah9wMtd6c1Rw4IpSYEAL3LbcxvKOm9OYqHA6jpKQk7fp5Adupjf3Upnq/X/fTEUwAT6xJz6NYSo0IMzMz3V5Ci5pXI/D5HoHfKnru1QGZmZkoKytLu35ewHZqYz+1qd7vuFwN152o4eEvbUTN9DsXy2k3jghdcOZrJhICqBxlQNPS7/AqERGp7auAQM/ZJh4r0XFLcXq9o1CpEWEoFHJ7CS1mwQ4bC3YK3NVfV35zFQqFUFxcnFb9vILt1MZ+avNCvxPyNVzVvf4oVsJOr6NYTrtJscHKyspyewkt5oEvbBS3AS4+Tu3NFVDfrby8PK36eQXbqY391OaVfnf0M7A1BLz4VXptsJx244iwBa3YK9D/FRPPDTdwzYlS7G2JiIiSdvHbJjYHBVaO8UFXfCrjlFIjwmAw6PYSWsTUFRaOzwV+eII3vgmDwSAKCwvTpp+XsJ3a2E9tXup3Z38dq2uB17emz1Esp92k2GBlZ2e7vYSU2+gXmL1J4Nf9dPh0b2ywsrOzUVFRkRb9vIbt1MZ+avNSvzOO1XHmsRoeWGHDxYFYi3LajSPCFnLTJxZe3Wpjyw99yPZ5Y4NFRET01jYbF82z8MFFBoZ3kuK4TUopNSIMBAJuLyGltocFnl5vY+JJuqc2V4FAAPn5+Z7v50Vspzb2U5vX+l3YRUO/o+vfxJUOnHaTYoOVk5Pj9hJS6pFVNrIMYEIfKb7czSYnJweVlZWe7+dFbKc29lOb1/ppmoY7+ht452uBz3d7f0zotBtHhCm2PyZw3L9M3NJHxwOnpdfF2IiIKD2YtkBRhYlTjtFQca7P7eWkFEeEkni8yoZpAxNPkuJL3awCgQA0TfN0P69iO7Wxn9q82M+na7j9ZANzNgusr/X2USylRoS5ubluLyElIqbAY1U2ru+lo0Mr75x7dUBubi6qq6s928/L2E5t7Kc2r/a7rqeGY1sBD31pub2UlHLaTYoNluq3jDmYp9ba2B8DftVXii9zs9M0Dfn5+Z7t52Vspzb2U5tX+x1laJh0ko5nNwh8HfbuUSyn3aR45ffCxda+y7QFylfauKq7hm753noSHRAMBlFQUODJfl7HdmpjP7V5ud9NvXW08gGPrfLuOwqVutBoXl6e20todi9vFtgSAn7dz7sntufl5cHv93uyn9exndrYT21e7pefqeHm3jqeWGPDH/fmUSyn3aTYYHnt6q9CCHCCtHwAACAASURBVDy0wsJ5nTX0P8abR6+A+scZCAQ81y8dsJ3a2E9tXu93W7GOmAU8ucabR7GcdpNigxUKhdxeQrN6b7vA8r3A7f2k+PKmTCgUQpcuXTzXLx2wndrYT21e79cpR8O1J2p4dJWNmOW9TaTTblLsALx2DayHVtg45RhgRCfvHr0C6rsJITzXLx2wndrYT23p0O9XJxvYHgFe2Oi9DZbTblJssCzLO2/pXL5HYP7XArf3Mzz3DpHvsiwLVVVVnuqXLthObeyntnToV9Raw+iuGv70pQXbY6NQp92k2GCFw2G3l9BsHv7SwvG5wJhu3t5cAfXdSkpKPNUvXbCd2thPbenS7/Z+OtbUAm9u89YGy2k33iqnGW0JCvR4ycSjJTpuKfbuuweJiIicGPqaCQ3AJ6O8c/scpW6VY5qm20toFo+stNE6E/hpTym+rClnmiYqKys90y+dsJ3a2E9t6dTv9n46FuwU+Gynd95R6LSbFDuBaDTq9hKO2N46gZnrbNxSrCMnw/vjQaC+W2lpqSf6pRu2Uxv7qS2d+l18nIai1sCfVnhng+W0mxTH7LxwsbW/rrYhBPDzPlLsWVtEXl4eampq3F4GJYHt1MZ+akunfrqm4dcnGxj3sYV1tQK9Wqt/AEKpC42qfpg0agpM/+9Nndtlq//N45Rpmpg3b57y/dIR26mN/dSWbv1+3KP+JtAPe+Qm0EqNCOvq6txewhF5boPAnjpgkkdv6nwwdXV1KCsrU75fOmI7tbGf2tKt31GGhtuKdTy3UWBnRP13FDrtxncRHiFbCPSpMFHcRsOc86SYuBIREUllf0ygywsmfnmyjj8MUPtd9kq9izCRSLi9hKS9tU1gnR/45clSfClbVCKRQEVFhdL90hXbqY391JaO/docpeH6Xjr+utpG1FT7KJbTblLsCuLxuNtLSNq0lTYGt9cwpIMUX8oWFY/HUV5ernS/dMV2amM/taVrv4kn6dgXA57doPY7Cp1244jwCCzbIzBgromKEQbGdE+/DRYREVFTjHnXxMp9AmtKfdAVvZ2cUiNCVXfx07600C0PuOx4Nb9JjlQ8HsfMmTOV7ZfO2E5t7Ke2dO73y7461vvVvn2O025SbLBUnENXhwRe2iQw8SQdhp6eG6x0PI/AK9hObeyntnTuV9JBR0l7DdNWqjsmdNqNI8Ik/XqRhZlrbVT/yIfcNLlyOxER0ZGas9nGmHctLL3UhwHt1Hv9VGpEGIvF3F5CkwTiAk+usXFTbz2tN1exWAzl5eXK9SO2Ux37qS3d+13aVUP3PGDaSjUvPOq0mxQbLMtS64v81DobERO4tViKL59rLMtCZWWlcv2I7VTHfmpL936GrmHiSTpmbxLYFlLvXCyn3TgibCLTFujxkomzjtXw7Nm8sCgREVFThRL1Fx69oZeOhwerdeFRjghTZM5mga0hoOxktb4hUiEWi+Gee+5Rqh/VYzu1sZ/a2A/IzdBwU28dM9baCMTVOoql1IjQttV5N8Ejq2yc00lD/2PS99yrA2zbRk1NjVL9qB7bqY391MZ+9W4t1hEx60+7UYnTbhwRNkHlThtDXrPw+vkGLu4qxd6UiIhIWdd8YOKznQIbrvQpc8kjpUaEqtxR/NFVNk7MB35wnBrfBKmWbneE9xK2Uxv7qY39/t+kkwxsDgKvbVVnTOi0mxQbLBVsCwnM2Szwi5N0ZS/vT0REJJMB7TQM7aDh0VVqjQmd4IjQodsXWZjBC4sSERE1qwMXHv38Mh9ObSv/66tSI8JoNOr2Eg4plBCYsdbGjUXpfWHR74pGoxg3bpz0/ej72E5t7Kc29vu20V01dM0FHlXkwqNOu0mxwdJ1KZZxUM9usBFMALek+YVFv0vXdRQWFkrfj76P7dTGfmpjv2/z6RpuK9bx4iaBbyLyn4vltBtHhIdhC4Gi2Sb6H6Nh9rm8sCgREVFz88cFCl8wMfEkHfcOlPs6k0qNCCORiNtLOKj/VAtsCACT+krxpZJKJBJBaWmp1P2ocWynNvZTG/t9X0Gmhp/21PHEGhtRU+6jWE67SbFrMAx5d6uPrrJxWjsNg9vz3KvvMgwDJSUlUvejxrGd2thPbezXuNuKdeytA17YKPcGy2k3jggPYeU+gZPnmHjhbANX95BiL0pERORZo98x8VVAYOUVPmiSXhJJqRFhOBx2ewmNemyVhc45wJjuckZ2WzgcxsiRI6XtRwfHdmpjP7Wx38FNPElH1X7gve3yHsVy2k2KDVZGRobbS/ie3VGB5zcK3NJHR4Yil+9vaRkZGSgtLZWyHx0a26mN/dTGfgc3vKOGk48GHlkp74VHnXbjiPAgpiy3MGV5/YVFj8niBouIiKglzFpn4/qPLay/0ocTC+R7/eWI8AgkbIG/rrZxzYkaN1eHEA6HUVJSIl0/Ojy2Uxv7qY39Du3qEzS0zQL+UiXnUSylRoSZmZluL+Fb5m4W2B4Bbi3mOzwOJTMzE2VlZdL1o8NjO7Wxn9rY79CyfBrGF+mYtd5GMC7fuVhOu3FE2Iihr5nI0IEPLuaFRYmIiFpaTUjg+BdNPFqi4xbJDnYoNSIMhUJuL6HBsj0Cn+4UuI23xTmsUCiE4uJiqfqRM2ynNvZTG/sdXmGuhsuP1/CXKhu2e8eBGuW0mxS7iKysLLeX0GB6lYXjcoFLuvLcq8PJyspCeXm5VP3IGbZTG/upjf2cue0kHev8wPwauTZYTrtxRPg/dkcFuvzLxB8H6Li9n1yHJImIiNKJEAID5pro2ErDmxfIc8qOUiPCYDDo9hIAADPW2tABjOslxZdFesFgEIWFhdL0I+fYTm3spzb2c0bTNNxabOCtaoENfnmOYjntJsVOIjs72+0l1F+aYY2NH/fQcDQvzeBIdnY2KioqpOhHTcN2amM/tbGfcwcu2fD4anku2eC0G0eE/zX7KxtXvW9hxeU+nHwMN1hEREQymLzEwvQqG1//yIe8TPdfn5UaEQYCAbeXgOlVdv0l+rm5ciwQCCA/P1+KftQ0bKc29lMb+zXNz3rriJjAsxvkOIrltJsUG6ycnBxX//5lewQW7BS4lZdmaJKcnBxUVla63o+aju3Uxn5qY7+m+b/27j04yvre4/jneTbEBNiEKrcQEkKAkIRi8AKC6DngjXgOHFoVb7WKEqTKmc6QM9pzZupAW3uGYwfOnKmtRsC7RQG13tEqFmyJrVQtcs8FSAIBRGFv2duzz/f8EbMVyeXZZHef5xc/rxlmFrKb/T15i/78/Z59no5LNvzaIZdssNqNW4QA7txqYMtRQcNNGcjgjZ2JiIgc5U/HTFz+WgybK12YU2DvYgi3CC06GRKsbxDcW6ZzcpUgr9cLTdO4zK0gtlMb+6mN/RI3c4SGC85rP53HbkptEQ4ePNi29163rz3WolJH/CiUMnjwYDQ3N9vaj3qH7dTGfmpjv8Rpmoal5e2XbDjotXeb0Go3R8wqNM2elaOYKXhkr4mbizUM5aUZEqZpGnJycmzrR73HdmpjP7WxX+/cMl7DkHOAR/bau4pltZsjJlh2XWztzWbBYT/w7zy5vVd8Ph9yc3N5sTwFsZ3a2E9t7Nc7AzM03FmiY91+E0HDvlUspS406na7bXnf3+wxMW2YhouHOeLHoBy32w2Px2NbP+o9tlMb+6mN/XrvnjIdp8LAC432TbCsdnPEzMKODzIeOC14u0WwtNwRPwIliQi8Xq8t/ahv2E5t7Kc29uu98bkaKgs0PLzbtO3nZ/V9HTG78Pv9aX/PR/aaGJoF3FjMPfDe8vv9KCgosKUf9Q3bqY391MZ+fbO0XMffTgr++rk9Eyyr3Rxxe+p0XwMrEBU8ccDEPWU6sjI4weqtnJwc/h+YothObeynNvbrm8rRGsa6gd/sNnHJ8PSvE1mdszhiBSsWi6X1/Z6rF/iiwI/KHHH4yorFYti9e3fa+1HfsZ3a2E9t7Nc3Ll3DPWU6XmgUnAimf6JqtZsjZhiBQCBt7yUi+M2eGOYWahjj5upVXwQCAcyYMSOt/Sg52E5t7Kc29uu7uybq0DVg3f70X7LBardOJ1h1dXW49NJLUVJSgqlTp2L37t2dvri2thZTpkzBlClTMGnSJCxZsgThcDjhwaZzi/DPxwU7vwRPbk+CjpuV2nmbI+odtlMb+6mN/fruvCwNt4zT8OheEzEzvatYfdoiXLJkCe6++24cOHAAP/nJT7Bw4cJOX1xRUYGPPvoIn376KT777DOcOHECv/3tbxMerGEYCb+mt36zx0RJLnBVPlev+sowDNTW1qa1HyUH26mN/dTGfsmxtNyFJj/welN6J1hWu501wTpx4gR27NiB2267DQBw/fXXo7m5GfX19We9eODAgRgwYAAAIBKJIBgM9urKtMFgMOHX9EZrm2BTo+Dech06r6DbZ8FgEAsWLEhbP0oetlMb+6mN/ZLjomEaLhmu4eE96d0mtNrtrAlWc3Mz8vLykJHR/gFDTdNQWFiIpqamTr/BoUOHUFFRgaFDhyI3Nxf33ntvl28WDofh9XrP+AUg/l6hUAihUCh+AB3bjW1tbfHHgUAAkUgk/jgajQJo/9hkx6zS5/PFH3u93vgJaQ//PYhMF3D7eC1+DRLTNOPjiMVi8ceGYcSv1moYRvxjmdFoNL7/GolE4o/D4TDa2trijzsCpPqYvF4vTNM847oq6Tomt9uNuro6ZGZm9ptj6o+dOjumzMxMtLS0QNf1fnNM/bFTV8eUlZWFlpYWaJrWb46pP3bq6piys7PR0tICEek3x2RXp6XlOt49Ivj7sfQd06BBg2CJfMOOHTukpKTkjD+bOnWqvPfee9986hl8Pp/Mnz9f1q9f3+Vzli9fLgDO+vWDH/xARESWLVsmy5YtExGRRYsWyfLly0VE5IYbbpBVq1aJiMg111wja9asERGR6dOny4YNG0REpLy8XDZv3iwiIvn5+bJ9+3YREXG73bJr1y6JxkzB/zTKrW+dFo/HIwDE4/FIc3OzdPwYdu3aJW63W0REtm/fLvn5+SIisnnzZikvLxcRkQ0bNsj06dNFRGTNmjVyzTXXiIjIqlWr5IYbbogf56JFi1J+TCIiAKS5udmWY4pGozJnzhx54IEH+s0x9cdOnR3TQw89JJs3b5arr7663xxTf+zU1TGtX78+/p795Zj6Y6eujmnbtm2yefPmfnVMdnUKGaYMfOSUTPrZW2k7pp07d8aPpzsQEXnqqaekoqJCKioqZOXKleJ2uyUajYqIiGmaMmLECKmrq+v2G4mIrF+/XubOndvl10OhkHg8nvivjh90Q0ODiIgEg0EJBoMiItLW1iahUEhERAKBQPyx3++XcDgcfxyJRESkfYLXMWav1xt/7PF4xDAMeflgTPBYRHYcN8Q0TfF4PGKapsRisfgPyTCM+ONoNCperzf+2OfziYhIJBIRv98vIiLhcDj+OBQKSSAQiD9ua2tL+TF1PI7FYrYck8/nk9LSUjl58mS/Oab+2KmzYzp58qSUl5fLsWPH+s0x9cdOXR3Tl19+KeXl5dLa2tpvjqk/durqmE6dOiXl5eVy5MiRfnNMdnaq/nNYhjwZlkDUTMsxdcxdeppgaSJnX+1s1qxZWLhwIRYuXIhNmzZh5cqV2LFjx1mrX/X19RgzZgwGDBiASCSCH/7whxg/fjx++ctfWlo983q9yM3NhcfjSfmnKea8acATBT6c74hrqxIREVESNHgF418w8Pg/uXDnxNRfIcDq3KXTkdTU1KCmpgYlJSVYuXIlnnjiifjXqqqq8OqrrwIAtmzZggsuuAAVFRW44IILMGLECDzwwAMJD7Zj/zNVGryCd44I7uGFRZMqGo1i48aNKe9Hycd2amM/tbFfco3L0TBndPslG9LBardOV7DSpWMWePToUeTl5aXsfe7/Swxr95s4cmsGsnlrnKQJBAK46qqr8O6771o/6Y8cge3Uxn5qY7/ke+WQie/9IYa/fT8DFw5N7X/nW1tbMWrUqB5XsBwxwUrlFmHIEIz+nYHbJ+hYPcOVkvcgIiIi+ximYOzzBipHa1jzT6k9FahPW4Tp1vFxyVTYdFDwRZj3HUyFSCSCtWvXprQfpQbbqY391MZ+yZeha1hcquN3DQJPJLXrRla7OWLWkcp96Ef2mrhylIaSIdwaTDaeR6AutlMb+6mN/VKjaqKOcAx4+kBqz8VS6hysVG0R7vxCUPGSgU1XuXD9WEfMJYmIiChFbnjXwJ5Tgt03ZPTqzjJWKLVF2JsbRFvxyF4TeQOBfxvD1atUCIfDWL16dcr6UeqwndrYT23slzr3lOnYexrYdix1a0dWuzligtVxCfpk8kUEz9abqJqoY4DOCVYqxGIx1NbWpqQfpRbbqY391MZ+qXPFKA0lucAjKbw/odVu/XaL8NE9MSzdbuLQzRkoGMwJFhER0bfB/34Ww0/+aqL5lgyMGJj8//5/q7cIRQSP7DUxr1Dj5CqFwuEwVqxYwWVuBbGd2thPbeyXWndM0OHSgMdTdLK7UluEppncH0LtCcHOL3lphlQzTRMtLS1J70epx3ZqYz+1sV9qnZul4aZiDTV7TcTM5G/SWe3WL7cIF/7RwLZjgvqbMqCn6FMERERE5EwfHjcx49UY3qp0obIguYstSm0RhkKhpH2vU2HBC42CxaU6J1cpFgqFUF1dndR+lB5spzb2Uxv7pd4lwzWcfy5Qk4L7E1rt5ogJVjI9W2fCMIE7S/rdoREREZEFmqbh7lIdrzUJjgbs2ajrV1uEIoLJLxooHaJh01WpvRcREREROZcnIhj1nIH/qtDx0wuTdy9ipbYIg8FgUr7P9uOC3aeAJaWOOKx+LxgMoqqqKmn9KH3YTm3spzb2S4/cTA03F2tYsz+5J7tb7eaImYiuJ2cYj+0zUewGrsznuVfpoOs6Ro8enbR+lD5spzb2Uxv7pc/dZTqa/MA7R5I3wbLard9sEZ4Kty8FLr9Qx39OSd5SIBEREalJRHDBSwaK3Bp+f01yTh1Saouwra2tz9/jGZ7cnnZtbW1YsGBBUvpRerGd2thPbeyXPpqmYUmZjtebBEeSdLK71W6OmI24XH1bcRIR1Ow18b0iLSWXxafOuVwuzJgxo8/9KP3YTm3spzb2S69bx+s4xwU8vj85l2yw2q1fbBH++ZiJy16L4Q//4sJV+Y6YMxIREZFDVG0z8IcjgsabMuDS+7YQo9QWYSAQ6NPra/aaGJfTfhdtSp9AIIA5c+b0uR+lH9upjf3Uxn7pt+Srk93fbun7mpLVbo6YYA0YMKDXr/0yJNhwULB4Iq/cnm4DBgzAggUL+tSP7MF2amM/tbFf+l08VMOU84CafX3fJrTaTfktwv/bFcN9fzHRcmsGhmdzgkVERERne3RPDEu3m2i6JQP5g3o/X/hWbBGKCB7ba+J7YzROrmwQCAQwY8YMLnMriO3Uxn5qYz973DpeR7YLWNfHk92V2iLMzMzs1etqTwj2nG6/kBilX2ZmJqqrq3vdj+zDdmpjP7Wxnz1yMjXcPE7Duj5e2d1qN6W3CO/aauCPrYL6mzJ4/hURERF168PjJma8GsPb17pwzejeLc4otUXo9/sTfo03InihUbCIJ7fbxu/3Y9KkSb3qR/ZiO7Wxn9rYzz6XDNfw3e8Aa/pwsrvVbo6YYGVlZSX8mvUNJkIxYCGv3G6brKwsrF69ulf9yF5spzb2Uxv72UfTNFSV6njlsOBEsHcbeFa7KbtFOPVlAyMHAq/NSc69hYiIiKj/+yLUfu/i/56q4z/OT/xq+kptEfp8voSe/+kXgh0nBVUTHTH8by2fz4fRo0cn3I/sx3ZqYz+1sZ+9zsvScN1YDWv2mejNGpPVbo6YoWRnZyf0/LX7TOQNBP61kOde2Sk7OxsbN25MuB/Zj+3Uxn5qYz/7LZ6oY78H+PPxxCdYVrspt0UYNAR5zxm4t1zHf0/ljTKJiIgoMaYIJrxg4PKRGp6cldipRkptEXq9XsvP3XRQ4IkAi7g9aDuv14ucnJyE+pEzsJ3a2E9t7Gc/XdOwaKKODY2C0+HE1pmsdnPELGXQoEGWn7t2n4krRmkYl8PtQbsNGjQItbW1CfUjZ2A7tbGf2tjPGRaW6IiY7VclSITVbkptER44LZi40cDvZrtwy3hHzA2JiIhIUfPfMdASEPzt+9ZvvN0vtwjX7jfxnXOA7xdx9coJvF4vNE3jMreC2E5t7Kc29nOOqok6Pj4JfHzS+lqTUluEgwcP7vE5kZjgqToTPxyvIyuDEywnGDx4MJqbmy31I2dhO7Wxn9rYzzmuLdAwamD76UdWWe3miAmWZuFWN683CU4EgapSRwyZ0N4tJyfHUj9yFrZTG/upjf2cI0PXcGeJjufqTbQZ1laxrHZzxGzFykW71u43MW2Yhsnn8h9Ip/D5fMjNzeXF8hTEdmpjP7Wxn7PcNVGHNwpsarQ2wVLqQqNut7vbr7f4BW+3CC/N4DButxsej6fHfuQ8bKc29lMb+zlLcY6G2XkaHj9gbZvQajdHzFh6+iDjU3UmztGBm8dx9cpJRARer7dXtxoge7Gd2thPbeznPItKdWxtFdR7em5itZsjJlh+v7/Lr5kieHy/iQXFGnIyOcFyEr/fj4KCgm77kTOxndrYT23s5zzXFWnIzQSesLCKZbWbIyZY3V1HYluroNHHK7c7UU5ODkTE8m2OyDnYTm3spzb2c57sDA23jtPx5AEThtn9CpXVbo6YtcRisS6/9vh+E+NzgMtHcvXKaWKxGHbv3t1tP3ImtlMb+6mN/ZzprokajrYB77R0P8Gy2s0RE6xAINDpn3sigk0HBXdN1PlxVgcKBAKYMWNGl/3IudhObeynNvZzpouGajj/XGDd/u63Ca12S+wW0inS1XLb8w0mwiZwxwRHzAPpG3izUnWxndrYT23s50zaVzeA/o8PTXweFAzL7nxhR6ktQsMwOv3zdfsF147WMGoQV6+cyDAM1NbWdtmPnIvt1MZ+amM/5/rBeB26BjxT1/UqltVujphgBYPBs/7ssy8FH33Oa185WTAYxIIFCzrtR87GdmpjP7Wxn3Odl6Xhe0Ua1u03u7wcg9Vumth4IY7u7ki9rDaG3zWYaLk1AwN0rmARERFR6r3TYmLOWzF8ON+FS4afvcjT3dzl6xyxPPTN5bZwTPDMVzd25uTKuQzDwNtvv81lbgWxndrYT23s52xXjtJQMKjrk92V2iIMhUJn/P61w4Ivwu33ByLnCoVCqK6uPqsfOR/bqY391MZ+zubSNdw5UcfzDYJA9OxNPqvdHLlFeO1bBk5HgNr5jviQIxEREX2LHPQKil8w8OQ/u3BHyZmLPUptEUaj0fhj3thZHdFoFBs3bjyjH6mB7dTGfmpjP+cbm6PhylEaHu9km9BqN0fMYiKRSPzxM/UmslzAjcU898rpIpEIVq9efUY/UgPbqY391MZ+arhzoo5txwSN3jM3+qx2c9QWoYigdKOBacM0PDOb24NERERkjzZDMPJZA9WTday4yBX/c6W2CDtmg385ITjgARaWOGJY1INIJIK1a9fy/8IUxHZqYz+1sZ8aBmZouKlYw1MHTJhfW4uy2s0RM5mO/cwnDwgKBgGzR3F7UAU8j0BdbKc29lMb+6ljYYmOQ35gW+s/JlhWuzlmi3DAQDfynjPw7+U6Hpzq6vnFRERERCkkIijZYGDmCA1Pzmo/dUmpLcJwOIxXDgs8EZz1cUhyrnA4jNWrVyMcDts9FEoQ26mN/dTGfurQNA13TNCx6aDA/9U1sax2c8RsJhaL4akDJmaO0DAhl9uDqojFYqitrUUsFrN7KJQgtlMb+6mN/dRye4mONgN48WD7BMtqN0dsEe49ehqT3hiIRy9zYXGpI+Z8RERERACAK98wYArw/twMtbYInzsQQabOa1+pJhwOY8WKFVzmVhDbqY391MZ+6llYouOPrYKDXlFri/D5Rg3XFWnIzeQESyWmaaKlpQWm2fkNMcm52E5t7Kc29lPPdUUaBg9ovxi61W6O2CLE/53EO9d9B1ePdsR8j4iIiOgMi7YaeL9V8PG1bfjOkCFqbBGOzBZcwWtfKYd3hFcX26mN/dTGfmq6o0THQR+wtVmhLcIbxwpcOidYRERE5EyXjdRQ7AaeP2htvuKILcIdh0/josJcu4ZBRERE1KOffxzDQ385jcC9Q3vcIrT1jsodc7vcyOfwermCpZpgMIj77rsPv/rVr5CdnW33cCgBbKc29lMb+6nr+yMEy30+AP+Yw3TF1hWsxsZGjBs3zq63JyIiIuqVhoYGFBcXd/l1W1ewzj33XABAU1NT+6cJSSlerxcFBQVobm7udpmUnIft1MZ+amM/tXk8HhQWFsbnMF2xdYKl6+3n2Ofm5vIfMoXl5OSwn6LYTm3spzb2U1vHHKbLr6dpHERERETfGpxgERERESWZa8WKFStsHYDLhVmzZiEjw9bdSuol9lMX26mN/dTGfmqz0s/WTxESERER9UfcIiQiIiJKMk6wiIiIiJKMEywiIiKiJEvLBKuurg6XXnopSkpKMHXqVOzevbvT561btw4TJkzAuHHjsHjxYkSj0XQMj3pgpd+WLVswbdo0lJeXY9KkSbj//vthmqYNo6Wvs/p3D2i/7cMVV1yBIUOGpHGE1B2r/T777DPMmjULZWVlKCsrw0svvZTmkVJnrPQzTRPV1dUoLy/H+eefj9mzZ6O+vt6G0dLX/fjHP0ZRURE0TcOnn37a5fO6nbdIGsyePVueeOIJERHZuHGjXHzxxWc9p7GxUfLy8qS1tVVM05R58+bJww8/nI7hUQ+s9Pv444+loaFBRESCwaDMnDkz/hqyj5V2HVatWiVVVVWSm5ubptFRT6z0CwQCMnbsWPnggw9ERMQwDDlx4kQ6h0ldsNLv5ZdflmnTpkkkEhERkV/84heyYMGCdA6TOrF161Zpbm6WMWPGyCeffNLpHKSJswAAA2hJREFUc3qat6R8gnX8+HFxu90SjUZFRMQ0TRkxYoTU1dWd8byHHnpIlixZEv/9G2+8ITNnzkz18KgHVvt909KlS2X58uVpGCF1JZF2u3btkssvv1zq6+s5wXIIq/3WrFkjt9xyix1DpG5Y7ff73/9eKioqxOv1immact9998myZcvsGDJ1orsJVk/zlpRvETY3NyMvLy9+rQhN01BYWIimpqYzntfU1IQxY8bEf19UVHTWcyj9rPb7umPHjmHTpk2YO3duuoZJnbDaLhqNYvHixaipqYHL5bJjqNQJq/327NmDc845B3PnzsWUKVNw++234/PPP7djyPQ1VvvNmzcPs2bNwsiRI5GXl4f33nsPP//5z+0YMiWop3kLT3KnpPJ6vZg3bx7uv/9+XHzxxXYPhyz42c9+huuuuw5lZWV2D4V6wTAMvPvuu6ipqcEnn3yC/Px83HPPPXYPiyzasWMHdu3ahSNHjuDo0aO48sor8aMf/cjuYVESpHyCVVBQgNbWVhiGAaD9RNqmpiYUFhae8bzCwkIcPnw4/vtDhw6d9RxKP6v9AMDn86GyshLz589HdXV1uodK32C13datW/HrX/8aRUVFuOyyy+D1elFUVMRVEJsl8u/O2bNnIz8/H5qm4bbbbsOHH35ox5Dpa6z2e/rpp+MfLtF1HXfccQfef/99O4ZMCepp3pLyCdbw4cNx4YUX4tlnnwUAvPjiixg9ejTGjx9/xvOuv/56vPrqqzh27BhEBI8++ihuvvnmVA+PemC1n9/vR2VlJSorK/HTn/7UjqHSN1ht98EHH+Dw4cM4dOgQ/vSnPyEnJweHDh3CsGHD7Bg2fcVqvxtvvBEfffQRvF4vAODNN99ERUVF2sdLZ7Lar7i4GFu2bEEkEgEAvP766/jud7+b9vFS4nqct6Tu1LB/2Ldvn0yfPl0mTJggF110kezcuVNERBYtWiSvvPJK/HmPPfaYFBcXS3Fxsdx1113xT1WQvaz0e/DBByUjI0MqKirivx588EE7h01i/e9eh4MHD/Ikdwex2u/pp5+WSZMmyeTJk6WyslKamprsGjJ9jZV+oVBIqqqqpLS0VCZPnixXX311/BPZZJ+7775b8vPzxeVyyfDhw2XcuHEikti8hfciJCIiIkoynuRORERElGT/D5uenTxE5JpCAAAAAElFTkSuQmCC\" />"
]
},
"execution_count": 37,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"plot(x, sol1(0.09))"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAYAAAByNR6YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xl4FeXd//HPzJyEBHISUBGB4IqIIAZBlAgKuEDdQJb4tJWiLbi0tVqo1adYC7bVx7rEttTHtqiIy4MFWcSl4AbIElFkBxcUkURkEeHs28x8f3/kx6lggDmQZOaefF7XxXUdjcF7ePfY+8z35D6aiAiIiIiIqN7obi+AiIiIyG+4wSIiIiKqZwG3F3Awtm1j27ZtCAaD0DTN7eUQERERQUQQiUTQrl076PrB71N5doO1bds2dOjQwe1lEBEREX1HdXU1SktLD/p1z26wgsEgAODTTz9F69atc/pe0xacP8fEyUENP+6k49ZlFgI68HgfA5eVcirakJLJJO69915MmDABBQUFbi+HcsB2amI3NbGbunbt2oWOHTtm9ykH49kN1r6xYDAYRHFxcU7f++RHNj7NWJjeL4BzjtMw4DTBTxZZGLFU8PMuOh48X0fzAMeODSE/Px/NmjVDcXEx/6OhGLZTE7upid3UlUwmAeCwb1/SvHpMQzgcRklJCUKhUE4brKQpOH26iQvaaPjXJf/ZP4oI/nejjTuW2zg5CDw/IIAex3GTRURERM453Z949g7WPolEIqcN1uMf2vgqDvyhp7Hf39c0DT/vauDidjquW2Di/Dkm/nCujl+frcPQudGqL4lEAr/4xS8wadIkFBYWur0cygHbqYnd1OTXbul0Gl988QUsy3J7KUfFMAycdNJJyM/P/87XEomEo9/D8xusQ71D/0CRtOD+1TZ+3ElDp5Z1b5rObKXh3SEBTPjAxvj3bbxWLXimv4GTg9xk1Qdd11FaWppTN/IGtlMTu6nJj91qamrwwx/+EPF43O2l1IvmzZtj2rRpaN++/X5/32kzX40I7/3Awv+ssfHptQGUFh1+w/TOVzZ+tNDCnhTwWB8DIztqPBKCiIgoR7Zt45ZbbsH27dvx+9//Xvn3lSWTSdxzzz1o27Yt/v73v++3qfLNiDAejzvaYH2dFDyyzsbPu+iONlcAcFFbHWuHa/jFUgujFlp4ZauGx/sYOKaAm6wjFY/Hcf3112Pq1Klo3ry528uhHLCdmthNTX7r9vXXX2PlypW477770L17d7eXUy9uvfVW3H333di9e/d+pxk4vUPn+XuThmEc/h8C8D+rbQDAb7rndkkl+RqeGRDACxcbeL1GcPYsE299aee8TqplGAbKy8sddyPvYDs1sZua/NZt7969AHDIc6FUs+9a9uzZs9/fd9rM8xusZs2aHfafqY4KHtto41fddBx3hHef/us0HeuGB3BGiYZLX7Pwq3ctJE1PTk89rVmzZhg3bpyjbuQtbKcmdlOT37rZdu2NCb9sGIH/XMu+a9vHaTPPb7Bisdhh/5k/rrIRzAPGdTu6yykt0vDGFQYqe+v42wYb571kYt033GTlIhaLYdCgQY66kbewnZrYTU3s5o4tW7agf//+KCkpOeJRptNmnt9g5eXlHfLrn4cFT31s464yHcH8o3/vlK5pGNvNwPvXBCACnDvbxKPrLNje/FkAz8nLy0NFRcVhu5H3sJ2a2E1N7OaO4uJi/PGPf8T//d//HfHv4bSZ5zdYdZ1B8W33rbZwTAHw0zPr91LOPlbD+9cE8PMuOsa9a2PgaxZqotxkHU5+fj7GjBlz2G7kPWynJnZTE7s1rIcffhg33XRT9q/37t2L4447DgDQt29ftGjR4oh/b6fNPL/BOtStuM/Cgqc/Efx3mY4WefX/k38FAQ2V5QbeuMLARyFBt5kmpn/GN8AfSiwWQ3l5OW97K4jt1MRuamoK3TaHBSu/rv9fm8OHv9kxZswYzJkzJ/vm+ylTpmDIkCE45phjjvq6nDbz/DENh9op/nGVhdYFwC31fPfqQJe217F2mIZbllj4r7ctvLLVxqQ+BkrqYSTpN/n5+Rg3bhxflSmI7dTEbmrye7evk7UfW2c3wODH0IDtIwOH/KG2li1bYsSIEXjqqacwduxYPP744/jXv/5VL/9+p808v8E62KxzU0jwzCbBo711FDbCBzcfU6DhX5cYuPpTwc+XWnhnu4ln+xu4sK3nbwI2qn3vKyD1sJ2a2E1Nfu92XIGGTdcGsDdd/793y3w4OjHgtttuw+DBg3HmmWeidevWOOecc+rl3+/0PVie32BFo9E6Dxr9w0oLJxQCN3VuvA2Opmn40ekaLjxBw6iFFvq9YuGuMsG9PXXkG7ybBdT2Ov/887F8+XIUFRW5vRzKAdupid3U1BS6nVrs7v8vdu7cGaeeeipuuukmPPjgg/X2+0ajUUf/nOdvv9R13P7HewXPfyYY311HQSPcvTrQyUENC640cH8vHQ+vtdH7JRMf7uEb4IHaXpWVlcp/TEJTxHZqYjc1sVvjuPHGG2GaJkaMGAGg9hT20tJSVFRUYOPGjSgtLcVvfvObnH5Pp808fwcrEPjuEn+/0kK75sCYRrx7dSBD1/Df3Q1c1l7HyIUmesw28fD5On7WRW/Sn2cYCAQwaNAgt5dBR4Dt1MRuamK3xrFgwQL87Gc/y471mjdvjpqamqP6Peval9TF83ewIpHIfn+9cY9g2meCu7vraOaBsVzP1ho+GBrA6DN03LrMxpXzLWyPN927WZFIBKWlpd/pRt7HdmpiNzWxW8Patm0bOnfujJUrV+KXv/xlvf7eTpt5foNVWFi431//fqWFDkXAT87wztKbBzT8rY+B175nYOXXtcc5zNnSNI9zKCwsxIwZM77TjbyP7dTEbmpit4bVrl07fPTRR1i2bBmCwWC9/t5Om3lnl3IQ374Vt/4bwfTNgru7G558U/nlHWo/z7BPGw1D37Bw4zsmopmmdTcrEAigvLzc8S1U8g62UxO7qYnd1OWbEWE4HM4+vnelhZOKgB+f4b3N1T6tCzXMvszAExcamPaZoPssE+/uaDp3s8LhMIqLi/frRmpgOzWxm5r81k3Xa7cTmUzG5ZXUn33Xsu/a9nHazPNb533H2a/ZLXjxc8GTFxnI0727wQJqj3MY3VlDv7YaRi600PdlC789R/Dbc3QEPL72o9WiRQtUVVUd1ccQkDvYTk3spia/dWvXrh3y8/MxefJk3Hjjjcp/xmImk8HkyZORn5+Pdu3a7fc1p800EW9+inE4HEZJSQlCoRCKi4sx7A0Ta3YLPro24PkN1reZtuC+VTb+sMrGua01PNvfwOkl6qyfiIjIiXfffRfjxo1DOt0Ap4u6ID8/H5WVlejdu/d+f//A/cnBeH6DVV1djV0F7dFjtomn+xm4vpPnp5p1eneHjZELLWyPA38uNzD6DM2Xxzk4/R8eeQ/bqYnd1OTXbtFoFNu2bYNtq/3WGF3X0a5duzoPga2pqUGHDh3U32Dt2bMH179XhA/3CDZWBJQesUUzgrFVFp74WDD4JA1PXGigdaG611MX27axbds2tGvX7jtza/I2tlMTu6mJ3dS1d+9etGrV6rAbLM9XXbUbmPuF4Hc9DKU3VwBQlKdh8kUBzL7MwLIdtcc5vLZV7V3+gTRNQ3FxsS/vzvkd26mJ3dTEbupy2szzG6z7PkjjjBLgB6f553+E15xce5xDj+M0XDnfws+XWoibnryRmLNIJIKSkhIenqcgtlMTu6mJ3dTlm4NG39puYEIPA4bid68OdEJzDa8OMvBYHx1PfWyj52wTK79Wf5MVDAYRCoXq/WA3anhspyZ2UxO7qctpM89vsM4oAa491V+bq300TcPPuhhYOTSA5gHg/Dkm/me1BctWd6MlIgiHw/DoW/voENhOTeymJnZTl9Nmnt9g3X5G2nd3rw50ZisNVYMDuONsHXe/b6P/Kxa2RNR80kWjUXTo0AHRaNTtpVCO2E5N7KYmdlOX02be/ynCvXvRsqTE7eU0msVf2fjRQgvfpIDH+hgY2dGfxzkQERGpyOkRG56/gyWKn6WRqwvb6lgzPIBrTtIwaqGF779t4ZukJ/fAdbIsCxs2bIBlWW4vhXLEdmpiNzWxm7qcNvP8BisWi7m9hEZXkq/hmQEB/OtiA298KTh7lok3v1RjoxmLxVBeXt4ku6mO7dTEbmpiN3U5beb5EaHfTrnNVU1UcMMiC29tE4w9S8f9vXQUBDgyJCIicoNvRoSmabq9BFeVFml4/QoDlb11PLbRRq85Jtbu9uSeGEBtr6qqqibfTUVspyZ2UxO7qctpM89vsBKJhNtLcJ2uaRjbzcCKoQFoGtBrjolH1lqwPXjzMZFIoKKigt0UxHZqYjc1sZu6nDbjiFAxSVPw2xU2HllnY0BbDVP7G+hQxJEhERFRY+CI0KcKAhoe7m3gzSsMfBIWnD3TxL8+884b4E3TxPz589lNQWynJnZTE7upyzcjwmQy6fYSPOmS9jrWDgtgYKmG779t4UcLTITS7t+MTCaTGDduHLspiO3UxG5qYjd1OW3GEaHiRATPfyr4+VILLZsBz/Y3cFFbz++biYiIlOSbEWEmk3F7CZ6maRpGnl57OOlJRRr6v2LhN+9ZSFvu7JszmQxmzJjBbgpiOzWxm5rYTV1Om3l+g5VOp91eghJODmpYcKWB+3vpeHitjd4vmfhwT+NvstLpNCorK9lNQWynJnZTE7upy2kzjgh9aOXXgusWmNgSAR4+X8fPuuj8PEMiIqJ64JsRIXf3uetxnIYPhgYw5gwdty6zccU8C1/FG2cfnU6n8cQTT7CbgthOTeymJnZTl9Nmnt9gcT59ZJoHNEzqY+C17xlYtVvQ7UUTc7Y0/HEOfF+ButhOTeymJnZTl9NmHBE2AbsSghsXW3jpC8HoMzT8udxAUR5HhkRERLnyzYgwlUq5vQTltS7UMPsyA09caOCFzwTdZ5mo2tEwd7NSqRQqKyvZTUFspyZ2UxO7qctpM89vsCzLcnsJvqBpGkZ31rF6WADHFWjo+7KFCR9YyNj1ewPTsixUVVWxm4LYTk3spiZ2U5fTZhwRNkGmLbhvlY0/rLLR8zgNzw0wcHoJR4ZERESHwxEhHVRA1zChp4ElVxv4JlU7Mpz8kY362GunUilMnDiR3RTEdmpiNzWxm7p8MyK0be98kLHf9G6jY9WwAK7rqOGmxRauecPCrsTRbbJs20ZNTQ27KYjt1MRuamI3dTltxhEhAQBe2mJjzGILhgY8dZGBK070/N6biIio0flmRMhPGm8cQ07WsW54AD2O03DlfAs/W2Ihbua+9+YnxKuL7dTEbmpiN3U5beb5DRY1nhOaa3h1kIHH+uh4+hMbPWaZ+GCXJ29wEhEReRpHhFSnj/bWfp7h2t3AvT113FWmw9D5k4ZERNS0+WZEmEgk3F5Ck9S5pYaqwQHcWabjtyts9HvFwufhw+/FE4kExowZw24KYjs1sZua2E1dTpt5foOl655fom/lGxru62Vg0VUGamKCslkmpn5y6OMcdF1HaWkpuymI7dTEbmpiN3U5bcYRITkSSgtuW2bhmU2CilM0/L2vgWMKODIkIqKmxTcjwng87vYSCEBJvoap/QP418UG3twm6DbTxJtffvcskHg8joqKCnZTENupid3UxG7qctrM8xsswzDcXgJ9y7Wn6Vg7LIAzW2q47DUL46osJL91nINhGCgvL2c3BbGdmthNTeymLqfNOCKkI2KL4C/rbfzmfRunFwPPDwjg7GM5MiQiIn/zzYgwFou5vQSqg65pGNvNwPvXBKBpQK85Jh5ZayESjWHQoEHspqBYjO1UxG5qYjd1OW0WaOB1HLW8vDy3l0CH0O0YDe8NCeC3K2zcsdzGK1/k4fIR17ObgvLy8lBRUcF2imE3NbGbupw244iQ6s3bX9oYtchCLAM83tfA90/z/A1SIiKinHBESI3u4vY63r08DeOjt/CDty2MXGBib8qT+3eqQywWQ3l5OZ9zimE3NbGbupw28/wGKz8/3+0lUA7aFOXjf8+JYMqFwMtf1B5Ouuir7x7nQN6Tn5+PcePG8TmnGHZTE7upy2kzjgipwXwREYxaaGHxdsGdZTp+31NHvsGfNCQiInX5ZkQYjUbdXgLlIBqNomvXrohGozgpqOHtKw38Ty8dlets9H7JxMY9ntzPE/ZvR+pgNzWxm7qcNvP8BqugoMDtJVAOCgoKUFlZme1m6Bru6m7g3SEBJCyg52wTf9tgHfLzDMkdB7YjNbCbmthNXU6bcURIjSZuCu5abuNvG218r1TDU/0MtG3OkSEREanjqEaEmzZtwgUXXIBOnTqhV69e2LBhQ53fbNs27rjjDpx11lno3LkzRo8ejXQ6nf361q1bcfXVV+OMM85Aly5dMGnSpJwvJBKJ5Pw95J5IJILS0tI6uzUPaJjUx8Br3zOwereg24smZn/ON8B7xaHakXexm5rYTV1Om9W5wbr55ptx00034ZNPPsFdd92FG264oc5vfvLJJ7Fy5UqsXLkSH374IXRdx1/+8hcAgIhg6NChGDVqFD7++GNs3LgR1157bc4XUlhYmPP3kHsKCwsxY8aMQ3a7vIOOdSMCuPAEDcPetDB6kYlI2pM3UpsUJ+3Ie9hNTeymLqfNvjMi3LlzJzp27IhvvvkGgUAAIoK2bdtiyZIl6Nix437ffOutt6Jdu3YYP348AGDWrFmYOHEi1q5dizfffBO/+93vsGzZsiO6AI4I/U9EMOUTwW3LLLQpBJ4bYKC8jeffFkhERE3YEY8Iq6ur0bZtWwQCtZ+io2kaTjzxRGzduvU739yzZ0/MnTsX4XAYmUwG06dPx5YtWwAAGzduROvWrfH9738f55xzDoYOHYrNmzcfdCGpVArhcHi/X0Dthg8AkskkkskkACCRSCCVSgEA4vF49nEsFsuOKGOxGDKZDIDad/ybpgmg9tbevsfhcBiWZWUf27YNEUE4HIaIwLbt7Dosy8o+Nk0ze4vQNM3sTxRkMpnsAWTpdDr7OJVKIR6PZx8nEgnfXtPOnTsRDAYRDocPe02apuG/ShNYMRg4vlBD35dN3PN+BhlbPHVNfuxU1zWFw2EEg0Hs2rXLN9fkx04HXtNXX32Vfc755Zr82OnAa9q9ezeKi4uxbds231yTHzvVdU3btm2DI3KAFStWSKdOnfb7e7169ZK33nrrwH9UbNuWCRMmSPfu3aW8vFzuueceadWqlYiIPPLII1JUVCTr168XEZHHH39cevbs+Z3fY58JEyYIgO/8GjlypIiIjB07VsaOHSsiIqNHj5YJEyaIiMiIESPkkUceERGRgQMHyuTJk0VEpHfv3jJ9+nQREenSpYvMmzdPRETat28vy5YtExGRYDCYXR8Aqa6ullAoJAAkFApJdXW17PsjWr9+vQSDQRERWbZsmbRv315ERObNmyddunQREZHp06dL7969RURk8uTJMnDgwOyfxYgRI7LXOXr0aN9e0+233y6jRo0S0zRzuqaMZUvrH/9Z9H8m5bzZGWnTrY9nrsmPneq6JtM0ZeDAgfLQQw/55pr82OnAa7rssstk4sSJYpqmb67Jj50OvKYXXnhB1q9f76tr8mOnuq6pbdu22TUfCkREpk6dKmVlZVJWViYPPPCABINByWQyIlK7iWrTpo1s2rTpkL+RiMi0adOkb9++IiIyY8YMufDCC7Nfi0ajommapNPpOr83mUxKKBTK/tr3B719+3YREUkkEpJIJEREJB6PSzKZFBGRWCyWfRyNRiWVSmUf7/t3RSKR7PWEw+Hs41AoJKZpZh9bliW2bUsoFBLbtsWyrOwfoGma2ceZTEbC4XD2cSQSERGRdDot0WhURERSqVT2cTKZlFgsln0cj8d5TQe5piXb0tLxhbQ0fzIlj69PZ9eu8jX5sROvidfEa+I1NdVrqqmpcb7BOlC/fv1kypQpIlK7UTrYnadEIiHffPONiIjs2rVLysrKZO7cudmFnXLKKVJTUyMitTvNfbtOJ/btaqurqx1/D7nv269GjlQkbcuN72QE/0zL4PkZ2Rm363GFdDD10Y4aH7upid3Ute8G0OHa1XkO1scff4wbbrghOyOeMmUKunXrBgAYM2YMBg8ejMGDB2PHjh3o378/dF2Hbdu4/fbbccstt2R/n9dffx133nknRAQlJSV47LHHsr/P4ex7E9mePXvQsmVLR99D7rNtG9u2bUO7du2g60f3hvWXttgYs9iCoQFPXWTgihP5BviGVJ/tqPGwm5rYTV179+5Fq1atDvsmd88fNLp3716UlJS4vRxySKT2DerBYBCadvSHiG6PC0a/Y+G1asFPz9TxcG8dzQM8nLQh1Hc7ahzspiZ2U1coFELLli3V/yxCHsKmlkgkgpKSknrrdkJzDa8MMvC/fXQ8/YmNHrNMrNjFw0kbQn23o8bBbmpiN3Ud1UGjXhIMBt1eAuUgGAwiFArVazdN0/DTLgZWDgugRR5Q/pKF+1dZsGxP3nxVVkO0o4bHbmpiN3U5beb5DZZHJ5h0EPKtM03qW+eWGqoGB3BnmY57PrDR7xULn4f5v4/60pDtqOGwm5rYTV1Om3l+g7XvQDFSQzQaRYcOHRqsW76h4b5eBhZdZeDLmKBslompn9j8j1Q9aOh21DDYTU3spi6nzTz/Jnd+VA4dTDhd+zE7UzcJRpyi4e99DRxbwDeLEhFRwznij8rxmn3H6pMaLMvChg0bGqVbcb6Gp/sHMP0SA29tE5w908QbNXwD/JFqzHZUf9hNTeymLqfNPL/B2vfZQqSGWCyG8vLyRu1WcaqOdcMD6NJKw8B/WxhbZSFpevLGrKe50Y6OHrupid3U5bQZR4TkG7YIJm2wcdd7NjoWA88PCKDsWI4MiYio/vhmRLjv06tJDaZpoqqqypVuuqbh9rMMrLgmAEMDzptj4uG1FmxvvobwHDfb0ZFjNzWxm7qcNvP8BiuRSLi9BMpBIpFARUWFq93OOkbDe9cE8IuuOu5cbuPSVy1UR7nJOhwvtKPcsZua2E1dTptxREi+9vaXNq5fZCGaAR7va+D7p3n+NQUREXkYR4TkCtM0MX/+fM90u7i9jrXDA/heBw0/eNvCyAUm9qY8+ZrCdV5rR86wm5rYTV2+GREmk0m3l0A5SCaTGDdunKe6tWqmYdrFATw/wMDLX9QeTrpwG49zOJAX29HhsZua2E1dTptxREhNyhcRwfWLLLzzleDXZ+v4/bk6mhn8SUMiInLGNyPCTCbj9hIoB5lMBjNmzPBst5OCGt66wsAD5+l4dL2N3i+Z2LjHk68xGp3X21Hd2E1N7KYup808v8FKp9NuL4FykE6nUVlZ6eluhq7hzjIDy4cEkLKAnrNNTFpvNfnPM1ShHX0Xu6mJ3dTltBlHhNSkJUzBXe/ZmLTBxqBSDVP6GWjbnCNDIiKqm29GhNzdqyWdTuOJJ55QplthQMNfLzDw7+8ZWLNb0O1FE7M+b5pvgFetHdViNzWxm7qcNvP8BovzabWo+r6C73XQsW5EABeeoGH4mxZGLzIRSXvy5m6DUbVdU8duamI3dTltxhEh0beICKZ8IrhtmYU2hcBzAwyUt/H86xAiImokvhkRplIpt5dAOUilUqisrFS2m6Zp+MkZOtYMD+D4Qg19X7bwuxUWMrYnX4fUK9XbNVXspiZ2U5fTZp7fYFmW5fYSKAeWZaGqqkr5bqcVa1h8tYEJPXTcv9pGn7kWPtnr702WX9o1NeymJnZTl9NmHBESHcbynTZGLrCwLQ482tvAjZ01aBp/0pCIqCniiJBckUqlMHHiRF91O/94HauGBTCyo46bl1gY8rqFnQlPvi45Kn5s1xSwm5rYTV2+GRHadtP8kXlV2baNmpoa33UrytPwjwsNvDTQwLs7Bd1mmnjlC39do1/b+R27qYnd1OW0GUeERDnaEReMfsfCq9WCW87U8fD5OlrkcWRIRNQU+GZEyE8aV0tT+IT4Ns01vDzIwON9dEz9xEaP2SZW7FL/VWhTaOdH7KYmdlOX02ae32AReZGmabili4FVwwII5mkof8nCfassmE3gOAciIjo8jgiJjlLGFtz7gY3/WWOj9/Eanu1v4NRijgyJiPzINyPCRCLh9hIoB4lEAmPGjGlS3fJ0DX/sZWDRVQa2xQTdZ5mY+okNj752Oaim2M4P2E1N7KYup808v8HSdc8vkb5F13WUlpY2yW59T6g9AX7YyRpuWGSh4i0Lu5PqbLKacjuVsZua2E1dTptxREjUAF7cbOPmJRaaGcDT/QwMLOV/RImI/MA3I8J4PO72EigH8XgcFRUVTb7biFN1rB0ewFmtNAz6t4VfVllImJ58LZPFdmpiNzWxm7qcNvP8BsswDLeXQDkwDAPl5eXsBqB9Cw3zLjfw53Idf//QRq85Jtbs9u4mi+3UxG5qYjd1OW3GESFRI1j/jeC6BSY+2gvc10vHuG46dH6eIRGRcnwzIozFYm4vgXIQi8UwaNAgdjvAWcdoeO+aAG47S8edy21c8qqFrVFvvbZhOzWxm5rYTV1Om3l+g5WXl+f2EigHeXl5qKioYLc6NDM0PHS+gTevMPBpWHD2TBMvfOadE+DZTk3spiZ2U5fTZhwRErlgT0rws6UWXvhM8MPTNDzWx0DLZhwZEhF5HUeE5IpYLIby8nJ2O4xWzTRMuziA5wcYeGVr7d2shdvcvZvFdmpiNzWxm7p8MyLMz893ewmUg/z8fIwbN47dHPphx9rjHE4t1nDxqxbuWm4hZblzU5nt1MRuamI3dTltxhEhkQdYtuCRdTZ+u8JG11bA8wMC6NKKI0MiIq/xzYgwGo26vQTKQTQaRdeuXdktR4au4c4yA8uHBJC2gJ6zTUxab8FuxNc/bKcmdlMTu6nLaTPPb7AKCgrcXgLloKCgAJWVlex2hM45TsOKoQHc1FnHbVU2rphnYVuscTZZbKcmdlMTu6nLaTOOCIk86vUaGzcsspCygMkXGhh2iudfDxER+Z5vRoSRSMTtJVAOIpEISktL2a0eDCzVsW54AP3bahj+poWfLDIRSTfc6yG2UxO7qYnd1OW0mec3WIWFhW40xXr6AAAgAElEQVQvgXJQWFiIGTNmsFs9ObZAw4uXGpjSz8CMzwXdZ5lYtqNhjnNgOzWxm5rYTV1Om3FESKSIzWHBjxZaeHen4O7uOu7poSNP508aEhE1Jt+MCMPhsNtLoByEw2EUFxezWwM4tVjDoqsM3NtDx/2rbfSZa+GTvfX3+ojt1MRuamI3dTlt5vkNVosWLdxeAuWgRYsWqKqqYrcGEtA1/LaHgWWDDexNCc6ZbeIfH1qojxvRbKcmdlMTu6nLaTOOCIkUFcsI7lhu4+8f2rjqRA1PXmTg+EKODImIGhJHhOSKcDgMTdPYrRG0yNPweF8DcwcaWL5T0G2miVe+OPI3wLOdmthNTeymLt+MCIuKitxeAuWgqKgI1dXV7NaIrj6p9jiH81pruPp1C7csthDL5H5jmu3UxG5qYjd1OW3m+Q2WpnHkoRJN01BcXMxujaxNcw1zBxp4vI+OZzbZOGeWifd35XY3i+3UxG5qYjd1OW3m+Q0WD2FTSyQSQUlJCbu5QNM03NLFwKphAZTka7jgJQt/XGnBtJ3dzWI7NbGbmthNXb45aDQYDLq9BMpBMBhEKBRiNxed0VLDsiEG/ru7jgkrbfR7xcLm8OE3WWynJnZTE7upy2kzz2+wPPpDjnQQIoJwOMxuLsvTNfzhXAPvXGXgq7igbJaJpz+xD9mF7dTEbmpiN3U5beb5DVY0GnV7CZSDaDSKDh06sJtH9DlBx+phAYw4RcOPF1moeMvC7mTd/3FgOzWxm5rYTV1Om/EcLKIm4sXNNm5eYqGZATzdz8DAUs+/viIi8hzfnINlWZbbS6AcWJaFDRs2sJsHjThVx9rhAZzVSsOgf1u4fZmFhPmf11dspyZ2UxO7qctpM89vsGKxmNtLoBzEYjGUl5ezm0e1b6Fh3uUG/lyu4x8f2Th3tonVu2s3WWynJnZTE7upy2kzjgiJmqgN3wiuW2Bi417gvnN1jOumw9B5Jg8R0aH4ZkRomqbbS6AcmKaJqqoqdlNA12M0LL8mgF+epeOu92xc/KqJ2Ys+YDvF8DmnJnZTl9Nmnt9gJRIJt5dAOUgkEqioqGA3RTQzNDx4voG3rjTwWUgwfO3JmPphyu1lUQ74nFMTu6nLaTOOCIkIALAnJfj5UgvTPhP84DQN/9vHQMtmHBkSEX0bR4TkCtM0MX/+fHZTUNCwcH3mLTxzEfBateDsmSYWbsvt8wyp8fE5pyZ2U5dvRoTJZNLtJVAOkskkxo0bx24K2tduaPsU1g4P4NRiDRe/auHO5RZSlidvdBP4nFMVu6nLaTOOCImoTpYtqFxn4+4VNrq0BJ4fEEDXYzgyJKKmzTcjwkwm4/YSKAeZTAYzZsxgNwUd2M7QNfy6zMB71wSQsYGec0z8db0F25uvyZosPufUxG7qctrM8xusdDrt9hIoB+l0GpWVleymoIO1636shhVDA7i5s47bq2xc/m8L22LcZHkFn3NqYjd1OW3GESEROfZ6jY0bFllIWcDkCw0MO8Xzr9GIiOqVb0aE3N2rJZ1O44knnmA3BTlpN7BUx7rhAfRvq2H4mxZ+sshEJO3J12hNBp9zamI3dTltVucGa9OmTbjgggvQqVMn9OrVCxs2bKjzm23bxh133IGzzjoLnTt3xujRo/f7F//pT39Cly5d0L17d/Tu3RvvvfdezhfC+bRa+L4CdTltd2yBhhcvNTCln4EZnwvKZplYup3HObiFzzk1sZu6nDarc0R48cUXY9SoUbjhhhvw4osv4k9/+hPef//973zz5MmTMW3aNMybNw95eXm46aab0KlTJ/z617/G6tWrMWTIEGzYsAFFRUV47rnn8Ne//tXxJosjQiLv2xwW/GihhXd3CsZ31/G7Hjry+HmGRORjRzwi3LlzJ1asWIGRI0cCAIYPH47q6mp8+umn3/nmNWvW4NJLL0V+fj40TcPll1+OZ599FgCgaRoymUz2U6f37t2L0tLSnC8kleLHdqgklUqhsrKS3RR0JO1OLdaw6CoDv++p44HVNvrMtfDJXo4MGxOfc2piN3U5bfadDVZ1dTXatm2LQCAAoHajdOKJJ2Lr1q3f+eaePXti7ty5CIfDyGQymD59OrZs2QIAKCsrw9ixY3HKKaegtLQUjz76KCZNmnTIBYfD4f1+Achu0JLJZPZwr0Qikb3AeDyefRyLxbIjylgslr2NF41GsyevRiKR7ONwOAzLsrKPbduGiCAcDkNEYNt2dh2WZWUfm6aJSCSSfRyNRgFgvw1lOp3OPk6lUojH49nH+z7HyI/XFIvFsGTJEliW5Ztr8mOnuq7JsiwsXrw4+887vSbYFu4+x8DrlyQRSgnOmW3izytjsCzL9WvyY6cDrykcDmPp0qWwLMs31+THTgdeUzKZRFVVFUKhkG+uyY+d6rqmUCgER+QAK1askE6dOu3393r16iVvvfXWgf+o2LYtEyZMkO7du0t5ebncc8890qpVKxER2bx5s/Tu3Vu+/PJLERGZNGmS9OnT5zu/xz4TJkwQAN/59aMf/UhERMaOHStjx44VEZHRo0fLhAkTRERkxIgR8sgjj4iIyMCBA2Xy5MkiItK7d2+ZPn26iIh06dJF5s2bJyIi7du3l2XLlomISDAYlPXr18v/H5NKdXW1hEIhASChUEiqq6tl3x/R+vXrJRgMiojIsmXLpH379iIiMm/ePOnSpYuIiEyfPl169+4tIiKTJ0+WgQMHiojII488IiNGjMhe5+jRo3lNvCZfXtN7azbILYtNwT/TcunsiGzaof41+bETr4nXxGs68mtq27Ztds2HAhGRqVOnSllZmZSVlckDDzwgwWBQMpmMiNRuotq0aSObNm065G8kIjJt2jTp27eviIg89NBDcuONN2a/Fo1GBYCkUqk6vzeZTEooFMr+2vcH/cUXX4iISCKRkEQiISIi8XhcksmkiIjEYrHs42g0mv39o9GopNNpERGJRCLZ6wmHw9nHoVBITNPMPrYsS2zbllAoJLZti2VZ2T9A0zSzjzOZjITD4ezjSCQiIiLpdFqi0aiIiKRSqezjZDIpsVgs+zgej/v2mvbu3St333139u/54Zr82Kmua0omkzJ+/Pjsv+torulfH0ak9TNpaf1MWl7YGGGnBrym3bt3y29/+1tJJpO+uSY/djrwmiKRiEyYMEG+/vpr31yTHzvVdU2bN292vsE6UL9+/WTKlCkiIjJjxgzp2bNnnd+cSCTkm2++ERGRXbt2SVlZmcydO1dERGbOnCldunTJ/kG88MIL37kzdij7drXbt293/D3kvng8LqNHj87+D5rUUd/ttsdsuWpeRvDPtNz8jinRtF0vvy/tj885NbGburZv3+5og1XnTxF+/PHHuOGGG7B7924UFxdjypQp6NatGwBgzJgxGDx4MAYPHowdO3agf//+0HUdtm3j9ttvxy233LJv9Ijx48dj9uzZaNasGVq0aIFJkyahZ8+ejkaX/ClCIvWJCP75kY1x79po3xx4boCB8473/PF7REQH5XR/4vmT3Hfs2IHjjz/e7eWQQ8lkEuPHj8f999+PgoICt5dDOWjIdp/sFYxcaGHl14IJPXT8pruOAI9zqBd8zqmJ3dS1c+dOtGnTRv2T3IlIfZ1aalg62MD47jomrrRx0csWPgt78rUdEVG98PwdLI4Iifxl2Q4bIxdY2JUE/lpu4IZOGjSNd7OISA2++SzCfedXkBoSiQTGjBnDbgpqrHYXtNGxZlgAFado+Mk7Fka8aeHrpCdf5ymBzzk1sZu6nDbz/AZL1z2/RPoWXddRWlrKbgpqzHbBfA1P9Qtg5qUGFn4lOHumifnV/DzDI8HnnJrYTV1Om3FESESu2hYT/HiRhde/FNzWVccD5+koDHBkSETe5JsR4b7j7kkN8XgcFRUV7KYgt9q1a6Hh35cb+Gu5jn98ZOPc2SZW7/bk6z5P4nNOTeymLqfNPL/BMgzD7SVQDgzDQHl5ObspyM12uqbhF2cZ+OCaAPJ04Lw5Jh5cY8GyudE6HD7n1MRu6nLajCNCIvKUlCW4Z4WNh9fauKithmf6GzixiCNDIvIG34wI9306NqkhFoth0KBB7KYgr7RrZmh48HwDb19pYHO49g3w//cp3wB/MF7pRrlhN3U5beb5DVZeXp7bS6Ac5OXloaKigt0U5LV2/dvpWDs8gCs7aLhugYUfvm1iT8qTN9xd5bVu5Ay7qctpM44Iicjzpn1q46dLLQTzgGf6GxjQzvOvDYnIpzgiJFfEYjGUl5ezm4K83O4HHWvvZnUs1nDJqxZ+vdxCyvLka8NG5+VudHDspi7fjAjz8/PdXgLlID8/H+PGjWM3BXm93YlFGt660sCD5+v4y3ob580xsf4bbrK83o3qxm7qctqMI0IiUs6a3YLrFpj4NAz86Twdv+iqQ+fnGRJRI/DNiDAajbq9BMpBNBpF165d2U1BKrUrO1bD+9cEcMuZOn5ZZeN7/7bwZcyTrxUbnErd6D/YTV1Om3l+g1VQUOD2EigHBQUFqKysZDcFqdauMKDhz+UG5l9uYP0eQbeZJl7c3PSOc1CtG9ViN3U5bcYRIREpb3dScNNiC7O2CG7opOEv5QaK8zkyJKL655sRYSQScXsJlINIJILS0lJ2U5DK7Y4t0PDipQam9DPw4ueCspkmlmxvGnezVO7WlLGbupw28/wGq7Cw0O0lUA4KCwsxY8YMdlOQ6u00TcMNnXSsGRZA+xYa+r1i4bfvW8j4/PMMVe/WVLGbupw244iQiHzHsgUPrLEx8QMb3Y/V8NwAA2e05MiQiI6eb0aE4XDY7SVQDsLhMIqLi9lNQX5qZ+ga7j7HwLIhBsIZwTmzTDy+0YJHX08eFT91a0rYTV1Om3l+g9WiRQu3l0A5aNGiBaqqqthNQX5s16u1jpVDA7i+k46fLbVx9XwLO+L+2mT5sVtTwG7qctqMI0IiahJe3WrjJ+9YEAGevMjA1Sd5/vUlEXkQR4TkinA4DE3T2E1Bfm935Yk61g0PoPfxGga/buHmxRZiGU++vsyJ37v5FbupyzcjwqKiIreXQDkoKipCdXU1uymoKbQ7vlDDSwMN/KOvgec+tXHOLBPv7VT7OIem0M2P2E1dTpt5foOl8fPFlKJpGoqLi9lNQU2lnaZpuOlMHauGBtCymYYL5lr4/UoLpqLHOTSVbn7Dbupy2szzGywewqaWSCSCkpISdlNQU2vXqaWGpYMNjO+u496VNi562cJnYfU2WU2tm1+wm7p8c9BoMBh0ewmUg2AwiFAoxG4Kaort8nQNvz/XwOKrDexICLrPMvHUx7ZSxzk0xW5+wG7qctrM8xsslf5DR7W9wuEwuymoKbe7oI2O1cMCqDhFw+h3LAx/08LXSTX+HJpyN5Wxm7qcNvP8Bisajbq9BMpBNBpFhw4d2E1BTb1dMF/DU/0CmHmpgUVfCbq9aGJ+tfffAN/Uu6mK3dTltBnPwSIiOsC2mODHiyy8/qXg1i46HjxfR2GAb0YmIh+dg2VZlttLoBxYloUNGzawm4LY7j/atdDw78sN/KVcx+SPbfScbWLV1558LcpuimI3dTlt5vkNViwWc3sJlINYLIby8nJ2UxDb7U/XNNx2loEPrgmgmQGc/5KJP622YHnsOAd2UxO7qctpM44IiYgOI2UJJnxg48E1Ni48QcMz/Q2cFOTIkKgp8s2I0DRNt5dAOTBNE1VVVeymILY7uGaGhgfOM7DgKgNbooKzZ5p4/lNvHOfAbmpiN3U5beb5DVYikXB7CZSDRCKBiooKdlMQ2x1ev7Y61gwL4OqTNIxcYOGHCyzsSbm7yWI3NbGbupw244iQiOgIvPCZjZ8usVCUB0ztZ+Di9p5/vUpE9YAjQnKFaZqYP38+uymI7XLz/dN0rB0ewOnFGi55zcId71pIWY3/epXd1MRu6vLNiDCZTLq9BMpBMpnEuHHj2E1BbJe7DkUa3rzSwEPn6/jrBhvnzTGx/pvG3WSxm5rYTV1Om3FESERUD9bsFly3wMSnYeCBXjpuO0uHrvEnDYn8xjcjwkwm4/YSKAeZTAYzZsxgNwWx3dEpO1bDimsC+OmZOsa+a2PQvy18GWv416/spiZ2U5fTZp7fYKXTabeXQDlIp9OorKxkNwWx3dErCGh4tNzA65cb2LhH0G2miRc3N+znGbKbmthNXU6bcURIRNQAvkkKbl5i4cXPBdefruGvFxgozufIkEh1vhkRcnevlnQ6jSeeeILdFMR29euYAg3TLzHwdD8Ds7YIymaaWLK9/u9msZua2E1dTpt5foPF+bRa+L4CdbFd/dM0Ddd3qj2ctH0LDf1esXD3+xbS9XicA7upid3U5bQZR4RERI3AsgV/WmNjwgc2yo7V8NwAA51bcmRIpBrfjAhTqZTbS6AcpFIpVFZWspuC2K5hGbqG8ecYqBpiIJoR9Jhl4vGN1lF/niG7qYnd1OW0mec3WJZlub0EyoFlWaiqqmI3BbFd4zi3tY6VwwL48Rk6frbUxlXzLeyIH/kmi93UxG7qctqMI0IiIpe8utXGT96xIAI8cZGBwSd5/jUvUZPHESG5IpVKYeLEieymILZrfFeeqGPd8ADK22gY8rqFmxabiGZye83LbmpiN3X5ZkRo2w17SB/VL9u2UVNTw24KYjt3HF+oYc5lBv55oYHnP619b9Z7O503YDc1sZu6nDbjiJCIyCM2hQQjF1j44GvB73roGN9dR0DnTxoSeYlvRoT8pHG18BPi1cV27ju9RMOSwQbuPkfHvSttXPiyhc/Ch34NzG5qYjd1OW3m+Q0WEVFTkqdruLengSVXG9iZEHSfZeKpj+2jPs6BiBoXR4RERB4VSQt+WWXhqU8EQ0/W8M8LDRxXwJEhkZt8MyJMJBJuL4FykEgkMGbMGHZTENt5TzBfw5P9Aph1qYF3vhJ0e9HEvOr932DLbmpiN3U5beb5DZaue36J9C26rqO0tJTdFMR23jX0FB3rRgRQdqyGy+dZ+MVSC3GzdvjAbmpiN3U5bcYRIRGRIkQEj2208evlNk4OAs8PCKDHcRwZEjUm34wI4/G420ugHMTjcVRUVLCbgtjO+zRNw61dDXwwNIACA+j9kok/vJ/EiIpr2U0xfL6py2kzz2+wDMNwewmUA8MwUF5ezm4KYjt1dGmlYfmQAMZ10zFhtY5VFz+ImgS7qYTPN3U5bcYRIRGRwhZ9ZWPUQgt7U8BjfQxc11GDpnFsSNRQfDMijMVibi+BchCLxTBo0CB2UxDbqenc4gROe+5aXNHewo8WWvjB2xb2pDz5upm+hc83dTlt5vkNVl5enttLoBzk5eWhoqKC3RTEdmrKy8vDD4deian9dEy72MD8GsHZM028/SU/487L+HxTl9NmHBESEflIdVRw/UILC74S/Kqbjvt66WhmcGRIVF84IiRXxGIxlJeXs5uC2E5NB3brUKThzSsNPHy+jkkbbPSaY2LdN558Hd2k8fmmLt+MCPPz891eAuUgPz8f48aNYzcFsZ2a6uqmaxp+dbaB964JwBbg3NkmHl1nwfbmwKJJ4vNNXU6bcURIRORjSVPwm/dt/Hm9jUvba3i6n4H2LTgyJDpSvhkRRqNRt5dAOYhGo+jatSu7KYjt1HS4bgUBDY+WG3jjCgMb9wi6zTQxYzPfAO82Pt/U5bSZ5zdYBQUFbi+BclBQUIDKykp2UxDbqclpt0vb61g3PIBL2mm49i0L1y80EU57coDRJPD5pi6nzTgiJCJqQkQEz24S3LrMwrHNgGcHGOh7gudfaxN5xlGNCG+77TacfPLJ0DQNq1evPuS/6Mknn8Tpp5+O0047DTfeeCMymYyjrzkViURy/h5yTyQSQWlpKbspiO3UlGs3TdMwqpOONcMCaN9CQ79XLNz9voW05cnX2r7F55u6nDarc4M1YsQILFmyBCeddNIhv/nzzz/HPffcg8WLF+PTTz/Fjh078M9//vOwX8tFYWFhzt9D7iksLMSMGTPYTUFsp6Yj7XZKsYZFVxn4Q08dD66xccFcCx/t5SarsfD5pi6nzercYF100UUoLS097De/+OKLGDx4ME444QRomoZbbrkF06ZNO+zXchEIBHL+HnJPIBBAeXk5uymI7dR0NN0MXcP4cwy8OySAaEbQY5aJ/91owaPvHPEVPt/U5bTZUQ3et27dut9drpNPPhlbt2497NfqkkqlEA6H9/sFADt37gQAJJNJJJNJAEAikUAqlQIAxOPx7ONYLIZ0Op19vG8kGY1GYZomgNpbe/seh8NhWJaVfWzbNkQE4XAYIgLbtrPrsCwr+9g0zewtQtM0sz9RkMlksgeQpdPp7ONUKoV4PJ59nEgkfHtNO3fuRDAYRDgc9s01+bFTXdcUDocRDAaxa9cu31yTHzsdeE1fffVV9jl3pNfU4zhg4aUJ3NBJw8+X2hj0agrb48JODXhNu3fvRnFxMbZt2+aba/Jjp7quadu2bXBEDuGkk06SVatWHfTrt956q9x///3Zv96wYYN06NDhsF+ry4QJEwTAd36NHDlSRETGjh0rY8eOFRGR0aNHy4QJE0REZMSIEfLII4+IiMjAgQNl8uTJIiLSu3dvmT59uoiIdOnSRebNmyciIu3bt5dly5aJiEgwGJT169eLiAgAqa6ullAoJAAkFApJdXW17PsjWr9+vQSDQRERWbZsmbRv315ERObNmyddunQREZHp06dL7969RURk8uTJMnDgQBEReeSRR2TEiBHZ6xw9erRvr+n222+XUaNGiWmavrkmP3aq65pM05SBAwfKQw895Jtr8mOnA6/psssuk4kTJ4ppmvVyTU9/sEPwcLUc90xa/vrO5+zUQNf0wgsvyPr16311TX7sVNc1tW3bNrvmQ4GIyNSpU6WsrEzKysrkqaeeyn7xcBusBx98UG6++ebsX7/66qvSp0+fw36tLslkUkKhUPbXvj/o7du3i4hIIpGQRCIhIiLxeFySyaSIiMRisezjaDQqqVQq+zidTouISCQSkUwmIyIi4XA4+zgUColpmtnHlmWJbdsSCoXEtm2xLCv7B2iaZvZxJpORcDicfRyJREREJJ1OSzQaFRGRVCqVfZxMJiUWi2Ufx+NxXhOvidfEa/LsNX26IySD52cE/0zLqDfjEknbyl+THzvxmty5ppqaGkcbrEMe03DyySdjzpw56N69e51f37x5M/r27YuVK1eiTZs2GDJkCAYOHIhbb731kF9zYt+PQVZXVzt6Pxh5A4/XUBfbqamhuokInvhY8MsqC+2aA88NMHD+8TzOob7w+aaumpoadOjQ4ciOabj55ptRWlqKmpoaDBo0CB07dsx+bcyYMZg7dy4A4NRTT8W9996LPn36oGPHjmjdujVuvvnmw34tF0VFRTl/D7mnqKgI1dXV7KYgtlNTQ3XTNA03dtaxelgAxzTT0GeuhXs/sGDafAN8feDzTV1Om3n+oNG9e/eipKTE7eWQQyKCSCSCYDAITePnnamE7dTUGN0ytuC+VTb+uMrGua01PNffQMcS/m/kaPD5pq5QKISWLVuq/1mEPIRNLZFIBCUlJeymILZTU2N0y9M1TOxpYPHVBr5OCrrPMvHkRzaPczgKfL6p66gOGvWSYDDo9hIoB8FgEKFQiN0UxHZqasxu5W10rBoawPdP0zBmsYVhb1j4OslN1pHg801dTpt5foPFV0hqkW+daUJqYTs1NXa3YL6GJy4KYNalBhZvF3R70cS/q+1G+Xf7CZ9v6nLazPMbrH0HipEaotEoOnTowG4KYjs1udVt6Ck61o0IoOxYDVfMs3DrUgtxk5sFp/h8U5fTZp5/kzt/hJWIyLtEBI9ttPHr5TZOCQLPDQigx3F80zb5l9P9iefvYO07Vp/UYFkWNmzYwG4KYjs1ud1N0zTc2tXAyqEBFBhA75dMPLDagsXjHA7J7W505Jw28/wGa99nC5EaYrEYysvL2U1BbKcmr3Q7s5WGd4cE8KtuOsa/b2PAqxa2RLjJOhivdKPcOW3GESEREdWrd76yMWqhhT0p4LE+Bq7rqPGsJ/IN34wI9316NanBNE1UVVWxm4LYTk1e7HZRWx1rhgcw+CQNP1po4QdvW9iT8uRredd4sRs547SZ5zdYiUTC7SVQDhKJBCoqKthNQWynJq92K8nX8OyAAKZdbGB+jeDsmSbe/pLHOezj1W50eE6bcURIREQNqjoquGGRhbe3CX7VTcd9vXQ0MzgyJDVxREiuME0T8+fPZzcFsZ2aVOjWoUjDG1cYePh8HZM22Og1x8S6bzz52r7RqNCN6uabEWEymXR7CZSDZDKJcePGsZuC2E5NqnTTNQ2/OtvA+9cEIAKcO9vEo+ss2N4cojQ4VbrRdzltxhEhERE1qqQpGP++jUfX27iknYan+xkoLeLIkNTgmxFhJpNxewmUg0wmgxkzZrCbgthOTSp2KwhoqCw38MYVBj7cK+g208T0z5rWG+BV7Ea1nDbz/AYrnU67vQTKQTqdRmVlJbspiO3UpHK3S9vrWDc8gMvaa/ivty2MWmAilPbkUKXeqdytqXPajCNCIiJylYjg2U2CW5dZOKYZ8Gx/Axe29fzrf2qifDMi5O5eLel0Gk888QS7KYjt1OSHbpqmYVQnHWuGBVDaQkO/VyyMf99C2vLk6/964YduTZXTZp7fYHE+rRa+r0BdbKcmP3U7pVjDoqsM3NdLx0NrbJTPNfHRXn9usvzUralx2owjQiIi8pwPdgmuW2BiaxR4uLeOn56p8/MMyRN8MyJMpVJuL4FykEqlUFlZyW4KYjs1+bVbz9YaVg4L4Mdn6Pj5UhtXzrewPe7J+wFHxK/dmgKnzTy/wbIsy+0lUA4sy0JVVRW7KYjt1OTnbs0DGh7rY+DVQQZWfl17nMNLW/xxnIOfu/md02YcERIRkeftSghuXGzhpS8EY87Q8Gi5gaI8jgyp8XFESK5IpSrqReAAAB0/SURBVFKYOHEiuymI7dTUVLq1LtQw+zIDky80MO0zQfdZJt7doe7drKbSzY98MyK0bXWfQE2RbduoqalhNwWxnZqaUjdN0zCms47VwwI4rkBD35ctTPzAgml7chBzSE2pm984bcYRIRERKce0BfetsvGHVTbOba3huf4GOpZwZEgNzzcjQn7SuFr4CfHqYjs1NdVuAV3DhJ4Glg428HWydmT4xEc2PHrP4Duaajc/cNrM8xssIiKigzn/+NqR4Q9O03DjYgtD37CwK6HGJov8jSNCIiLyhTlbbIx5x0JAB6b0M3B5B95DoPrnmxFhIpFwewmUg0QigTFjxrCbgthOTez2H9ecrGPdiADOOVbDFfMs3LrUQtz05D0EdlOY02ae32DpuueXSN+i6zpKS0vZTUFspyZ221/b5hpe+56BSRfoePJjGz1nm1j5tfc2WeymLqfNOCIkIiJf+nBP7ecZrvsG+MO5On59tg5D508a0tHxzYgwHo+7vQTKQTweR0VFBbspiO3UxG4Hd2YrDe8OCeDXZTrGv29jwKsWtkS8cU+B3dTltJnnN1iGYbi9BMqBYRgoLy9nNwWxnZrY7dDyDQ339zKw6CoDW6OCs2eaeHaT+8c5sJu6nDbjiJCIiJqEUFrwi6UWnv1UcO2pGh7vY+CYAo4MKTe+GRHGYjG3l0A5iMViGDRoELspiO3UxG7OleRreGZAAP+62MAbXwrOnmXirS/d+agadlOX02ae32Dl5eW5vQTKQV5eHioqKthNQWynJnbL3bWn6Vg7LIDOJRoufc3Cr961kGzk4xzYTV1Om3FESERETZItgr+st/Hf79k4oyXw/IAAuh3DkSEdGkeE5IpYLIby8nJ2UxDbqYndjpyuaRjbzcCKoQEAwLmzTVSutWA3wn0HdlOXb0aE+fn5bi+BcpCfn49x48axm4LYTk3sdvS6HaPhvSEB3NpVx6+W27jsNQs10YbdZLGbupw244iQiIjo/3vrSxvXL7IQywD/6Gvg2tM8fx+CGplvRoTRaNTtJVAOotEounbtym4KYjs1sVv9uqR97RvgB5Zq+K+3LYxaYCKUrv/7EOymLqfNPL/BKigocHsJlIOCggJUVlaym4LYTk3sVv+OKdDwwsUGnulv4KUvBGUzTSz+qn6Pc2A3dTltxhEhERHRQWyJCEYttLBku+CuMh339tSRb/AnDZsy34wII5GI20ugHEQiEZSWlrKbgthOTezWsE4OalhwpYH7e+l4eK2N3i+Z+HDP0d+XYDd1OW3m+Q1WYWGh20ugHBQWFmLGjBnspiC2UxO7NTxD1/Df3Q0svyaAhAX0mG3ibxuso/o8Q3ZTl9NmHBESERE5FDcFdy638dhGG98r1TCln4ETmnNk2JT4ZkQYDofdXgLlIBwOo7i4mN0UxHZqYrfG1Tyg4W99DLz2PQOrdgu6zTQxZ0vub4BnN3U5beb5DVaLFi3cXgLloEWLFqiqqmI3BbGdmtjNHZd30LFueAB92mgY+oaFMe+YiGacD4TYTV1Om3FESEREdIREBE99LLi9ysIJzYHn+hvo3cbz9y7oKHBESK4Ih8PQNI3dFMR2amI3d2mahtGddaweFsBxBRr6vmxh4gcWMvah712wm7p8MyIsKipyewmUg6KiIlRXV7ObgthOTezmDR1LNCy52sA95+j44yobF75sYVPo4JssdlOX02ae32BpGn86QyWapqG4uJjdFMR2amI37wjoGib0NLDkagNfJwXdZ5mY/JFd53EO7KYup808v8HiIWxqiUQiKCkpYTcFsZ2a2M17erepHRn+8DQNNy22cM0bFnYl9t9ksZu6fHPQaDAYdHsJlINgMIhQKMRuCmI7NbGbNxXlaZh8UQBzLjOwbEftcQ6vbf3PcQ7spi6nzTy/wfLoDznSQYgIwuEwuymI7dTEbt425OTa4xx6HKfhyvkWfr7UQtwUdlOY02ae32BFo1G3l0A5iEaj6NChA7spiO3UxG7ed0JzDa8OMvBYHx1TPrbRY5aJxVtj7KYop80CDbyOo8YzsNRSXFzMV2SKYjs1sZsaNE3Dz7oYuLidjusWmLj0zQLcv8pEiyLP3+egAzjdl3i+rGVZbi+BcmBZFjZs2MBuCmI7NbGbWjq31FA1OIA7ugF3v2+h3ysmtkS4QVaJ0+ea5zdYsVjM7SVQDmKxGMrLy9lNQWynJnZTT76h4TdnJlDwt6tQHbFx9kwTz26q+zgH8h6nzzV+VA4REZFLQmnBbcssPLNJUHGKhr/3NXBMAc/G8jLffFSOaZpuL4FyYJomqqqq2E1BbKcmdlPTvm4tdAtT+wfwr4sNvLmt9jiHN/9fe3cfHFV973H8c/YsT0KCeJHCJEBMwEoAE3mSGLEE5cGWlArEq9XaXgEpt1xmyK3YmcoAithCydTBtgJiK+qt1/AkBQq9Co2kRBANQlCvEB4SFAiDspvd7NM553v/yGXLQ0LOCuzZ3+bzmsnMxizJ9/A2zs/z2z3nC6vlb0COsfu7lvALrEAg4PQIFINAIICioiJ2UxDbqYnd1HRptwezXNg/0Y1+N2oYvcVEcYWJoJGQG0ytnt3fNW4REhERJQhLBC9UWfjFHgu3dgbeKHDj9n/hlmEi4RYhOcIwDGzbto3dFMR2amI3NTXXzaVpmD1Qx94H3NA0YOgGA0v3m7AS81xIq5Q0W4TBYNDpESgGwWAQxcXF7KYgtlMTu6mppW4Db9KwZ4Ib/9HfhZ/vtjB6i4laHxdZicDu7xq3CImIiBLY9i8sPFZmwh8BXrpbx79mJfy5kaSWNFuEkUjE6REoBpFIBKWlpeymILZTE7upKZZuo9Ia72c4Jl3DQ9tN/GiHAU84Ic+NtAp2f9cSfoEVDoedHoFiEA6HUVJSwm4KYjs1sZuaYu3WpZ2GN0fpeG2kjo3HBbevNfDeSV7OwQl2m3GLkIiISCHH6wU/+ruJ8lOCp3JcWDDYhbY632kYL1e1RThr1ixkZGRA0zTs27fvij9o1apV6Nu3L7KysjBt2rToqbPt27dj2LBhyM7ORv/+/TFnzhxYVuyrbf5fmVrC4TBefvlldlMQ26mJ3dR0Nd16p2jY8T0di4a6sPSAheFvG/j064Q8V5KU7DZrcoE1efJklJeXo3fv3lf8w0ePHsXcuXOxc+dOHD58GKdPn8aKFSsAAF26dMGbb76JTz75BB9++CF27dqF1atXx3gYfA2Wavh6EHWxnZrYTU1X2013afhFro73J7gRMIFB6w387qDJ+xnGgd1mV9wizMjIwIYNG5Cbm9vk15csWYLq6mq89NJLAIAtW7Zg0aJFKC8vv+y5M2fORNeuXTF//nxbg3GLkIiIqGUNhuCp3RZe/MTCuHQNr3xHR48buGV4vcTlXYQ1NTUXneXKyMhATU3NZc87deoU1qxZg/Hjxzf7vUKhELxe70UfAHDu3DkAjdedOH/tiUAggFAoBABoaGiIPvb7/dFTd36/P7rK9Pl80QuD1dfXRx97vV6Yphl9bFmNdzP3er0QEViWFZ3DNM3oY8MwUF9fH33s8/kANK5qz99lOxwORx+HQiE0NDREH5+/zH4yHpPH48Gvf/3r6D9LhmNKxk5NHVMoFMKvfvWr6M9KhmNKxk6XHtNXX32FxYsXIxQKJc0xJWOnS4/J5/OhpKQEZ8+evepjaqdZeG6gH1vG6ag8KxhQGsH6oxY7XadjOnv2LGyRK+jdu7dUVlY2+/WZM2fKokWLop8fPHhQevbsedFzPB6PDBkyRJYuXXqlHyXz5s0TAJd9PPzwwyIiMnv2bJk9e7aIiEyZMkXmzZsnIiKTJ0+Ofu8xY8bIypUrRURk+PDh8tZbb4mISHZ2tmzdulVERNLS0mTXrl0iIpKSkiJVVVUiIgJAamtrxePxCADxeDxSW1sr5/+KqqqqJCUlRUREdu3aJWlpaSIisnXrVsnOzhYRkbfeekuGDx8uIiIrV66UMWPGiIjI0qVLZfLkydHjnDJlStIe08yZM6Vv377i9/uT5piSsVNTx+T3+yUtLU2ef/75pDmmZOx06THde++9MnjwYPH7/UlzTMnY6dJjeu2112Ty5MnSr1+/a3pMm9/bI+1n/0WwIixj3zgmt+UMZqdrfEzdu3ePznwlEBF59dVXJScnR3JycuSVV16JfrGlBdbixYtl+vTp0c83b94s+fn50c+9Xq/k5eXJs88+e8UhRESCwaB4PJ7ox/m/6FOnTomISCAQkEAgICIiDQ0NEgwGRUTE7/dHH/t8PgmFQtHH4XBYRETq6+slEolEZzr/2OPxiGEY0cemaYplWeLxeMSyLDFNM/oXaBhG9HEkEhGv1xt9XF9fLyIi4XBYfD6fiIiEQqHo42AwKH6/P/q4oaGBx8Rj4jHxmHhMPKbrckwej1de/tSUjq+EJfO/QlJxylT+mBKp04kTJ+wvsJrT0gKrurpaevToISdPnhTLsqSwsFCWLVsWHeauu+6SBQsWXHGA5pxf1dbV1X2jP0/OCAaDMm/evOi/rKQOtlMTu6kpHt0OnbPkzg0R0VeGZd5eQ8Kmdd1+VmtSV1dna4HV5Guwpk+fjvT0dJw4cQJjx45Fnz59ol+bOnUqNm7cCADIzMzEggULkJ+fjz59+uDmm2/G9OnTAQAvvPAC9uzZg3Xr1iE3Nxe5ubl47rnn7O1bXuCbXNqBnGNZFk6cOMFuCmI7NbGbmuLRrU9nDeWFOube4cLCSgsj/mLikIfvMrxadpvxQqNERERJbnedhUd3mPiyAfhtno6p39agaXyn4TeRNPci5B3i1dLSHeIpcbGdmthNTfHudmc3FyonuvFIHw1P7DTxg/8xcSaQkOdXEp7dZgm/wCIiIqKr16mNhhUj3NgwWseu04IBaw1sruHW8vXCLUIiIqJW5lSDYMp7JrbUCmb0c+E3w124wc0tQzuSZovw/AXCSA2BQABTp05lNwWxnZrYTU1Od+t+g4ZNY3X8Pt+FP31uYdA6A3vP8GyWHXabJfwCy+VK+BHpAi6XC+np6eymILZTE7upKRG6aZqGGdk6PproRqc2GvLeNvFcpQnTSsiNrYRhtxm3CImIiFq5sCl45iMLz39sIa+bhtdG6rgllVuGTUmaLcLz9xMiNTQ0NKCoqIjdFMR2amI3NSVat7a6hoVDdZSN1/GFX5CzzsCrnzfeK5AuZrdZwi+wdF13egSKga7ryMvLYzcFsZ2a2E1Nidrt7u4u7JvkxgMZGn5SZuLBd02cDXKRdSG7zbhFSERERJcpPWJhermJDjrwp+/oGJ2e8Odk4iJptgj9fr/TI1AM/H4/xo4dy24KYjs1sZuaVOhWlOnCgUluZHfRMOavJmZXmAgaCXlOJq7sNkv4BVabNm2cHoFi0KZNGxQVFbGbgthOTeymJlW6pXXUsO1+HSXDXfjDpxaGbjDw8dnWvciy24xbhERERNSiA18JHtlh4H/PAYuGujB7oAuuVng/Q24RkiP8fj/y8vLYTUFspyZ2U5OK3QbepOGDH7gxa4ALP99t4b7NJmp9CXmO5rpKmi3Ctm3bOj0CxaBt27YoLi5mNwWxnZrYTU2qdmuna1hyp453v6vjkFdw+1oDb1a3rivA223GLUIiIiKK2dchwYxyE/99RPBIHw0v3qXjxnbJv2WYNFuEPp/P6REoBj6fD/3792c3BbGdmthNTcnQrUs7DX8epeP1Ah1/Od54cdKyk8l/Nstus4RfYLVv397pESgG7du3R0lJCbspiO3UxG5qSpZumqbhkT4u7J/kRkYnDQWbTPxij4mwmZCbY9eE3WbcIiQiIqKrZlqC3+y3MPdDC/27AG8UNF5DK9kkzRZhfX290yNQDOrr65Gens5uCmI7NbGbmpKxm+7S8FSujvcnuBEygcHrDSyrMpPufoZ2myX8AqtDhw5Oj0Ax6NChA0pLS9lNQWynJnZTUzJ3G9RVw4cPuDH12y7MqrBw/1YTX/qTZ5Fltxm3CImIiOi62Fpr4d/KTEQsYMUIHRNvSfjzOi1Kmi1Cr9fr9AgUA6/Xi9TUVHZTENupid3U1Fq6jevpwoHJbozormHSOyYeLzNQH07I8zq22W2W8Ausjh07Oj0CxaBjx46oqKhgNwWxnZrYTU2tqVvX9hrWjdax6h4dpUcFuesM7Dqt7uUc7DbjFiERERHFRbVX8KMdJnafEfwy14W5g1xo41LrnYbcIiRHeL1eaJrGbgpiOzWxm5paa7esVA3vFeqYP8iFRfss5G808fm5hDzP06yk2SLs1KmT0yNQDDp16oTa2lp2UxDbqYnd1NSau7ldGuYO0vGP7+v4OiS4Y72BFZ9aylzOwW6zhF9gaZpapw5bO03TkJqaym4KYjs1sZua2A24s5sLlRPdeLSPC9PLTUz4m4m6QOIvsuw2S/gFVjJdhK01qK+vR+fOndlNQWynJnZTE7s16tRGw/IROt4eo6OiTjBwrYHNNYn9AvikudBoSkqK0yNQDFJSUuDxeNhNQWynJnZTE7td7Pu9XTgwyY0hXTWM32ZiRrkJfyQxz2bZbZbwCyxV9mSpkYjA6/Wym4LYTk3spiZ2u1z3GzRsGqvj9/kuvPq5hUHrDew9k3hns+w2S/gFls/nc3oEioHP50PPnj3ZTUFspyZ2UxO7NU3TNMzI1lE50Y2UNhry3jbxXKUJ00qchajdZrwOFhERESWcsCl45iMLz39sIa+bhtdG6rgl1fk3BSTNdbBM03R6BIqBaZo4ePAguymI7dTEbmpit5a11TUsHKqjbLyOL/yCnHUGXv3c+cs52G2W8Assv9/v9AgUA7/fj7y8PHZTENupid3UxG723d3dhY8nuTExQ8NPykw8+K6Js0HnFll2m3GLkIiIiJSw5oiFJ8pNdNCBP31Hx+j0+J8nSpotQsMwnB6BYmAYBioqKthNQWynJnZTE7t9M5MzGy/nkN1Fw5i/mphdYSJoxPc8kd1mCb/ACgQCTo9AMQgEAigqKmI3BbGdmthNTez2zaV11LDtfh2/zXPhD59aGLLBwMdn47fIstuMW4RERESkpKqvBI/sMPDZOeC5oS4UD3TBdZ1vP8QtQnKEYRjYtm0buymI7dTEbmpit2tjwE0a9vzAjVkDXJiz28J9m03U+q7veaOk2SIMBoNOj0AxCAaDKC4uZjcFsZ2a2E1N7HbttNM1LLlTxzvf1XHI23g/wzerr98V4O024xYhERERJYWvQ4J//4eJN6sFP8zS8Lt8HTe2u7ZbhkmzRRiJRJwegWIQiURQWlrKbgpiOzWxm5rY7fro0k7Dn0e58UaBjk01gtvXGvj7l9f2bJbdZgm/wAqHw06PQDEIh8MoKSlhNwWxnZrYTU3sdn39sI8L+ye5kZmqYdRmE0/tNhEyr82Gnd1m3CIkIiKipGRagqUHLDy910L/LsAbBY3X0LoaSbNFyNW9WsLhMF5++WV2UxDbqYnd1MRu8aG7NMzJ0bF7ghthExi83sCyKvOq7mdot1nCL7C4P60Wvq5AXWynJnZTE7vF1x1dNex9wI1pt7kwq8LC/VtNnGz4Zossu824RUhEREStxrZaCz8pMxGxgBUjdEy8JbZzTUmzRRgKhZwegWIQCoVQUlLCbgpiOzWxm5rYzTlje7pwYLIb9/TQMOkdE4+XGagP2z/XZLdZwi+wTNN0egSKgWmaqKioYDcFsZ2a2E1N7Oasru01rL1Pxyv36Cg9KshdZ2DXaXuXc7DbjFuERERE1GpVewWP7jCx54zgl7kuzB3kQhtX8+805BYhOSIUCmH+/PnspiC2UxO7qYndEkdWqoadhTrmD3Jh0T4L+RtNfH6u+XNPSbNFaFnX735CdO1ZloUTJ06wm4LYTk3spiZ2Syxul4a5g3Ts+r6Or0OCO9YbWP5p05dzsNuMW4RERERE/88XEfzn+xZWfGZhfC8Nq+7R0a3DP7cMk2aLkHcaVwvvEK8utlMTu6mJ3RJXpzYalo/QsXGMjt11goFrDWw6/s+zVnabJfwCi4iIiCjeCnu7cGCSG0O7aij8m4kZ5Sb8EfubftwiJCIiImqGiGD5pxaK37fQsxOwfIgfBVldWlyfuOM4Y0zOr/vq6uocnoRiEQgE8OSTT2LJkiXo0KGD0+NQDNhOTeymJnZTyw/TgaH3CabtNHHvGg8AtHg/w4Q9g3XkyBFkZWU5PQYRERHRZaqrq5GZmdns1xP2DNZNN90EAKipqUHnzp0dnobs8nq96NmzJ2pra7m1qxi2UxO7qYnd1OXxeNCrV6/oOqU5CbvAcrkaX3/fuXNn/sunoNTUVHZTFNupid3UxG7qOr9OafbrcZqDiIiIqNXgAouIiIjoGtPnz58/3+khmqPrOkaOHAm3O2F3MqkJ7KYutlMTu6mJ3dRlp13CvouQiIiISFXcIiQiIiK6xrjAIiIiIrrGuMAiIiIiusYcX2AdOnQId911F2699VYMHToUBw8ebPJ5q1atQt++fZGVlYVp06YhEonEeVK6kJ1u27dvx7Bhw5CdnY3+/ftjzpw5sCyrie9G8WT3dw5ovBXEqFGjcOONN8ZxQmqK3W4HDhzAyJEj0a9fP/Tr1w/r1q2L86R0ITvdLMtCcXExsrOzcfvtt6OgoACHDx92YFo6b9asWcjIyICmadi3b1+zz7vi2kQcVlBQIH/84x9FRKS0tFSGDBly2XOOHDkiPXr0kJMnT4plWVJYWCgvvvhinCelC9np9tFHH0l1dbWIiAQCAcnPz4/+GXKOnXbnLV26VKZOnSqdO3eO03TUHDvd/H6/3HLLLbJz504RETEMQ+rq6uI5Jl3CTrf169fLsGHDJBwOi4jIs88+K0VFRfEcky5RVlYmtbW10rt3b6msrGzyOS2tTRxdYJ0+fVpSUlIkEomIiIhlWfKtb31LDh06dNHzFi9eLNOnT49+vnnzZsnPz4/rrPRPdrtd6mc/+5nMmzcvDhNSc2JpV1VVJSNGjJDDhw9zgeUwu91WrlwpDz/8sBMjUhPsdtuwYYPk5OSI1+sVy7LkySeflNmzZzsxMl3iSgusltYmjm4R1tbWokePHtHrSGiahl69eqGmpuai59XU1KB3797RzzMyMi57DsWP3W4XOnXqFNasWYPx48fHa0xqgt12kUgE06ZNw/Lly6HruhOj0gXsdvvkk0/Qrl07jB8/Hrm5uXjsscdw5swZJ0Ym2O9WWFiIkSNHonv37ujRowfeffddPPPMM06MTDFoaW3i+GuwKPl5vV4UFhZizpw5GDJkiNPjkA0LFizAxIkT0a9fP6dHoRgYhoF33nkHy5cvR2VlJdLS0jBjxgynx6IW7N27F1VVVfjiiy/w5Zdf4t5778VPf/pTp8eiq+ToAqtnz544efIkDMMA0PiC2pqaGvTq1eui5/Xq1QvHjx+Pfn7s2LHLnkPxY7cbANTX12PcuHGYMGECiouL4z0qXcJuu7KyMixbtgwZGRm4++674fV6kZGRwbMhDonlv5UFBQVIS0uDpml49NFH8f777zsxMsF+t9WrV0ffTOJyufDjH/8YO3bscGJkikFLaxNHF1jdunXDoEGD8PrrrwMA1q5di/T0dPTp0+ei502aNAkbN27EqVOnICJ46aWX8NBDDzkxMsF+N5/Ph3HjxmHcuHF4+umnnRiVLmG33c6dO3H8+HEcO3YM5eXlSE1NxbFjx3DzzTc7MXarZ7fbgw8+iA8++ABerxcAsGXLFuTk5MR9Xmpkt1tmZia2b9+OcDgMANi0aRMGDBgQ93kpNi2uTa7fS8Ps+eyzz2T48OHSt29fGTx4sOzfv19ERKZMmSJvv/129HkrVqyQzMxMyczMlMcffzz6bgtyhp1uCxcuFLfbLTk5OdGPhQsXOjk2if3fufOOHj3KF7knALvdVq9eLf3795eBAwfKuHHjpKamxqmRSex1CwaDMnXqVLnttttk4MCBMnr06Og7sMkZTzzxhKSlpYmu69KtWzfJysoSkdjWJrwXIREREdE1xhe5ExEREV1j/wcgg3QhGmyvRgAAAABJRU5ErkJggg==\" />"
]
},
"execution_count": 46,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"first_order_coeffs = [-1.83333, 3.0, -1.5, 0.333333]\n",
"derivs_start = []\n",
"derivs_end = []\n",
"\n",
"for i in 0:1/99:1\n",
" push!(derivs_start, sum(first_order_coeffs.*sol1(i)[1:4])*100)\n",
" push!(derivs_end, sum(first_order_coeffs.*sol1(i)[end-3:end])*100)\n",
"end\n",
"# derivs\n",
"# plot(x,derivs_start)\n",
"plot(x, derivs_end)"
]
}
],
"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
}
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment