Skip to content

Instantly share code, notes, and snippets.

@shivin9
Created July 12, 2017 03:19
Show Gist options
  • Save shivin9/85750faa1f6f4cbe37a61d756002a5a0 to your computer and use it in GitHub Desktop.
Save shivin9/85750faa1f6f4cbe37a61d756002a5a0 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Plots.PyPlotBackend()"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"using DifferentialEquations, Plots, PDEOperators\n",
"pyplot()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"N,M = 100,100\n",
"Δx = 1/(N-1)\n",
"Δt = 1/(M-1)\n",
"\n",
"x = collect(-10 : Δx : 10);\n",
"ϕ(x,t) = (-1/2)*sech.((x-t)/2).^2 # solution of the single forward moving wave\n",
"u0 = ϕ(x,0);\n",
"du3 = zeros(x);\n",
"temp = zeros(x);\n",
"A = LinearOperator{Float64}(1,4,Δx,length(x),:periodic,:periodic);\n",
"C = LinearOperator{Float64}(3,4,Δx,length(x),:periodic,:periodic);\n",
"\n",
"function KdV(t, u, du)\n",
" C(t,u,du3)\n",
" A(t,u,du)\n",
" @. temp = -6*u*du - du3\n",
" copy!(du,temp)\n",
"end\n",
"\n",
"single_solition = ODEProblem(KdV, u0, (0.,5.));"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"\u001b[1m\u001b[33mWARNING: \u001b[39m\u001b[22m\u001b[33mInterrupted. Larger maxiters is needed.\u001b[39m\n"
]
},
{
"data": {
"text/plain": [
"retcode: MaxIters\n",
"Interpolation: 1st order linear\n",
"t: 4-element Array{Float64,1}:\n",
" 0.0 \n",
" 0.01 \n",
" 0.02 \n",
" 0.0202497\n",
"u: 4-element Array{Array{Float64,1},1}:\n",
" [-9.07916e-5, -9.17133e-5, -9.26443e-5, -9.35847e-5, -9.45347e-5, -9.54944e-5, -9.64638e-5, -9.7443e-5, -9.84321e-5, -9.94313e-5 … -9.94313e-5, -9.84321e-5, -9.7443e-5, -9.64638e-5, -9.54944e-5, -9.45347e-5, -9.35847e-5, -9.26443e-5, -9.17133e-5, -9.07916e-5]\n",
" [-0.0123044, -0.001352, 0.0104434, -0.00980894, 0.00222443, 0.00389734, -0.00537904, 0.0029489, -0.000889321, -0.000593093 … 0.0252626, 7.18769e-5, -0.0239393, 0.0260886, -0.00715337, -0.0158951, 0.023269, -0.0115338, -0.00777372, 0.0175189] \n",
" [-0.280718, 0.44042, -0.238846, -0.156845, 0.422687, -0.342131, -0.0175155, 0.361252, -0.407939, 0.12014 … 0.269768, -0.450895, 0.252584, 0.15694, -0.432802, 0.34736, 0.0260869, -0.37536, 0.413361, -0.110848] \n",
" [-0.178055, 0.131725, 0.0232274, -0.156411, 0.154309, -0.0200335, -0.127421, 0.157255, -0.0450099, -0.104917 … 0.15682, -0.0680692, -0.0704418, 0.163932, -0.140101, 0.0104468, 0.132335, -0.176403, 0.0826366, 0.078972] "
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"soln = solve(single_solition,solver=CVODE_BDF(),dense=false,saveat=0.01,maxiters=100000)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAYAAAByNR6YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzsnXl8VOX1/z93lpBkliD7EkABEYiKoKgBV0RiFVHR1A0VK6it9tuaX6vW1kprq7bV+O23tVXBIlZrNS6odUFxoS4RQVE0IouCJIDsZJbMdmee3x9DkpnJnbn3ztyZ+zwz5/168XJMZibPvfdzzz3POec5j8QYYyAIgiAIgiAMw2L2AAiCIAiCIIoNcrAIgiAIgiAMxmbmH4/FYti+fTtcLhckSTJzKARBEARBEJphjMHr9WLIkCGwWHrGq0x1sLZv345hw4aZOQSCIAiCIIisaW1tRXV1dY+fm+pguVwuAMCmTZvQv39/M4dCCEAwGMRvfvMb3HHHHSgvLzd7OATnkF4IPZBeCs9xz0ewsT3+un2u3dzBZMHu3bsxevToLl8mFVMdrM60oMvlgtvtNnMohACUlZWhV69ecLvdZAAJVUgvhB5IL4VnnxQBKuKv3W7xHKxgMAgAaUucJDPbNHg8HlRVVaG9vZ0cLIIgCIIoEZ7+OoaL34p2/T+bL56DpebDcLGKMBAImD0EQgACgQDmzZtHeiE0QXoh9EB6KSyr9ojfglNNK1w4WErV9wSRisViQXV1NemF0ATphdAD6cVcVuyImT0E3ahphVKEBEEQBEEUlJ+vjOLetclOlWhpQiFShB0dHWYPgRCAjo4O1NfXk14ITZBeCD2QXgpLMXS+VNMKFw6W1Wo1ewiEAFitVtTW1pJeCE2QXgg9kF4IvahphVKEBnLokxHcd6IVFx7Ghd9KEARB6OSTPQzDnUC/8mKIsfDLRctlPLs52f1g8+1oDzO8s51htFtCTR++r4EQKUK/32/2EHTx69VRXPm2nPSzV1tj+NYH/Gp1NM2niFzx+/2oq6sTTi+EOZBeCD106uXY52VMfVFW/wCRNRsOsB7OVSdjnpZx/htRHPks/9dAzbaY2mi0E7tdnMI2b5jhzjXxwrz6kTGcOyLuo579GjlW+cZut6O+vl4ovRDmQXoh9GC329F60cMAgA3tJg+myNkZSJ842yVQlww128JFBKusrMzsIWhm1utRxdedfHUAiMbE7+/BI2VlZZg3b55QeiHMg/RC6KGsrAzrMKTr/3dlcAKI3EjT+Fw41GwLFw6WSCH8D3ep33T/19K99HQ33aSG4ff7UVtbK5ReCPMgvRB6SNXJde9SViJfFIl/pWpbuHCwim2GufNgiPPzfQwDHpfxaqt4DdR4pKysDA0NDUWnFyI/kF4IPaTqJERm23T+/AXfTq4QESxRaiQ6ZIZgyvX2hpUjVIwxHH2wSO+zvRTFMgKqqSH0QHoh9EA6KRzpUoSzliUXtt/6Ed9erhA1WD6fz+whaOKWlT0v9oXLo/jnxuSfSwASy7D+8FkMJnbDKBp8Ph9qamqE0QthLqQXQg+kE/N5aatYz0k1zXDhYJWXl5s9BE3sDva8+Gv3MVz5Ts8wZuI7D4SBdQdy//s7OhiCslgCNJLy8nI0Njam1cvdn0axqb10zw+RjJpeCCKRVJ282srQniZDQeSG1hos3mu11GwLFw6WzcZFt4isUBKAJCU7WABgRABryBMyLlrOd046n9hsNtTV1Snq5Z5Po7htVQyzl/PfO4UoDJn0QhCpzFnR82f/2sR3ikpUvthv9giMQc22cOFgeb1es4eQNUq5ZAnGOFRKvNzKSrYNhNfrRXV1taJefrEqbggpE0t0kkkvBJHKU9+Q8SgU1xbJCk0128KFg1VRUWH2EDSh5ExZFH5216cxhA2e+Dy/ufsLf/9pac6qKioq0NTUlFEvX+wHlm8rzfNDJKNFLwSRCd5TVHo5fqmMw56MmD2MokHNtnDhYIkQwg/KDP/+uucMJ90N6Ho0OVW1S6F+Sw+zE1KDpboq0Wazoba2VlUvb28vzfNDJKNVLwRRKqzazbBFoFr+QBSIcZyWECJF6PF4zB6CKk1p9k3SOsOZ9nJxhETNxOPxwO12Y+seT1ID18/3JV+bYpt1EtnRqRcR7AvBJz98P1Y0C4s+3ClmZP8Pn/E7bjXbwoWD5XA4zB6CKkqpQABopSbRBcPhcKC5uRnjXqrAgMe7I4RHC7ApKFF4OvUign0hzGH9AYb1BzI7ULuDBRpMnql9sXuSv26/OE4jzxkbNdvChYNltVp1f6ZDZripOYpQtDAnn6Ii5mO1WjF23Hh0HPSntvn5vfEI87FaraipqcnKvhClwdgmGWObMk/QlGx/v8ciuGaFuBO78c+IM/Z0wY1c2dHB4Ivk9gxRsy1cOFjZhPD/sT6G//0ihufSpO6MhicH67ktDN91lJ5z4fF4YJu/pOv/hz+pbCSKZSNRIjc8Hg8kSaIUIaFKarPoRJYr1HTuDQH/2MDw8Dp+01fFQjpz/vk+1qM8RCvjmyIY8oSMk17MzdEUIkXodDp1f6az7u2yt6MFaVuQLy86EzNekfHQOuXarWVtpedgOZ1O2Gov6/r/dJf9d2tiOc9MCPFxOp1obW3Nyr4QpYVSs+hOrl6R/nf3c75XHo/o3dVEacJ8IBTfii6xPGR3gGF/SNt3dzb+/myfrqH0QM22cOFgSRpCDnuDDBOejWDHwcjNmoS8bEcBop1mBEXe2MZw/Xsx/OKjnjdxKQZptOikky8Ozmy2+sjRKlUkSYLb7dalG6J02OzRbhuaEwrE71pDTlUu6LXISk7KiS90P/RDUYZ9QYYBj8vo+5i6M5Bac7dyV/ZRSDXbwoWDpaUR4GttDGv3xbuZb/YwLN7QfZIKYT+XZAgha+XZzdl9xz0KqyjMembctSaKXo+Y00fF6/VClrV505IErN4dw4gnZby8lcL4pYjX60VVVRU1GiUUGaejDmlKQoH4L1cXhz1Zu5fhmW8KfyzjVGreUlHKHq1v73591qtRHHkwkqXkvPkjDF/u7972KLX+7EfvZ+8wC9Fo1OVy6Xr/Vh3FzdeskNGUo4h8EYZXWnOPhBi5zY0Z/hVjDL9cHW+i+oQJW0i4XK4efUfCaRY5SAAWfhX/3TqVVUJEceJyudDe3q7bvhClQUinOb7qnZ6OQbZ2+IfvRXOKnOhFqWRiwnMy6t+MYmOB92/d0K7+nkRSgwmpfbHe2cGwoyP9589+LYqaZ2T0XiIf/HzK92e4iowxPLYhhj+nSQWr2RYuHCy1nCxjDD/5oPsArVLq75U/9+6OGP6xgeH7b2bn2Gz1MVz4hgzZhAnLzz7MPGYzIlidDgsAzHm78GHyp7+JQWbJBz71ReVxyAx4+Kv4haP9xEoTxhg8Ho/umg+CUOKxjT11tO5A+kleJh5cF8OJL0Txn28LY5tuzBClGfM03ysKH93AsDNhUdfv1mQ+Z5tSHMb/fpf5+mR6lP7wvRiuWhHFT5tjaFUoN1GzLVw4WD5f5tayfjm+aqOT1JBhutN9yn+6RfWKzjRRNMYw4kkZz21hOOa5wgvwvs8zj1fNv3IsjkBaGME72427gd9Q2YKGMZbXrruPr++Zmly9R70B7Jq9wBYvPWRLDZ/Ph2HDhqnaF4LQyh6FHTnUbHUqiYuy3tlRGLuUKcIjAm8lrOT8z9bM5+y4pemf10pOUqZgReckHVBeta5mW7hwsNxud8bfpzpUqYNeoUGkem+CxFV633Jony9XWT3ZWfh/u4H1AmorKe/5LAbrInVn9MOdMchZrPwss2vf8iR1qBEKYpUcbrcbjDFV+0IQWjl3Wc9IUHtY33fYHsltwv6bj6OQFkZwweuy5uisGavgjcSuw1NJvB4vpUQIlZykVbsZPlGYqPsiLG1B/u4AQ4esblu4cLCi0czpplQN1aakhc57XT1dpTew8n8t/D+R093YiTl1I2+sj9NEi4B4Xvy2Vern7MOdMdS+GMVvPtF/finVQ+ghGo2ipaVF1b4QhFY+3KW8H+0nexjeaIvpbnz9ty9juu3agoO2c+m3DMOelDX9TR4W0m72MAx5IrsFUp1lQdv8DKt2az9fD2rsU6a0EOqm5p52wxtmWLw+hgGPx3toqdkWLhwsv7/nfjMBmWHs0/E0lyfPi9YYY0n7NDHG8tZn6ksDtyhId9Mk5tT/+x3Lutt9YpRp+bYYvk7pqfbXlm5x/f3L7vOXqRdJp3P8jY4l0p3k8qAk36z08Pv9qK2tVbQvBGEUFgk49nkZM16N1+roIRBVdtrS8e6O5O/f5gfK/yFjg8pCHh4e9Ne9F806VWk7eACL1+cn8KEUiNilsEXSJW9F8YP/xp9Da/YCm/dmti08nHfFMNvb21nXUsxsVjnctzb5Yfz2DuUwoD/CYFkko/bFeNj1i30Mp/0nfzPeGgO3KNA6KXkoi27DGw4w2B+RIS2MoD3MsDPQ8z0//qD7e7cn3Dh9NPQi+dfXrKsu6kCI4bzX43+r898P34uiI2WTVT0pwpU6ZjlEcdK50TOlCAkAWLU7hpmvyVi6JWboThiJdjibzuJRBvz4/fjzJ7VAO5VT0jybjkjT+iASY2hojuJllVXwv/0kv1HeLV6GN7Zlf85tlniW5NcfF87BUnq+pvoQP/64MvP35jAmw1DqbXTrqu4LfqpOh6fPYxH8bGXPC/Hudz1/9vQ3ySfsqGdl1VUHRvP29uxqkrR+Qu8qSMZY0g3be4kMex5CzB/uim/5c8hjMl78NvloHlwXg2OxjLaEokQ9m36mziR3BsjhKjVkWUZzc7Pm3mlEcXP80rijccEbUVyRoXN7oZEA/PVgBuDfX2fvQMz/b0+dT3xOxv1fqH/nHXlyXDrRW6eWypIN8frefFnxVEeIMYb3d/b8a9+lBBo6e2tp/V5TCASSR80Yw+dZtrBnjGF/SPl3Sj5CZ7jPLDa1M0x7OYqqJXJX0brW2ZVS2ktpybBFijc5/VpjWk4po2jLg1IsAAY/kfnhN+xJGftDDMu3xbAxhy3l0s38iOIlEAigvr6+h30hCJ+BZSe+BBP2/k6Wsb/Vq62ZHZnbP471KLH4cj+DtDCi2qh60fr451r2sa4SjJb9GT9SMJSCG3p46ht9rlVEZ8AiNYL1aivDHoUUYSr+cObj4sLBSm3WlU3J0N++jMIfYT0iUjySWNR4+8fxB3+HDLx5cCnqUc9qm3ErXdozX+npSHy6l+Gi5VEcn2H5aiJrFZzbCw1okvpVSp2AVaP6fvJBFNupjIbQicvlQltbGzUaJfDohmRr6TVwr9I/p0SIMq1qP/u1nnb0pJeSf7Y1ZdV6Z1lJpv0SO3lzWwxHPitj1FP5jdr6Iwxfe1iPMo50JJaTFIJfrYovHtDaIPzmj2I4kODYKpXEKLHFn/khxoWDlRrC19tKHwBueD+GBZ/EcMlb6UX4QkoaKnE/o0LS8GG32P79dfeYOh1LLZ4zAPT/Z8/xK6U3lxxskLcvFL8B1Tj2eX3nJcYY7vo0+Xs/2BmDtDCSJNrU63rvWm033T83Mc3OmNE8tC6K3ZReFBJZlrFs2TJKERJ4KiX1ls/ITq7VFOmsjZY9d6cnTLCzXdykxAtb4vb8rFfjgzj95ShGP6Vt7z8zWHeA4bkt+o4/sQDfqFWXXDhYwWCyR7Epy1SQWj70re0sqaZnpY7VG0bSuc1Mat1VNivd2nRuZpw62zICpcUDnR3WU1ceJqJn9YzebS2MYE8wvtn2/HcpvSgiwWAQDQ0NPewLUXoUsk1Brn/rsre6nZZsusR3ctarxtgtX4Th/Dfi37WsjeGhddGuVgnBqL7a2EJhkYADaUqF0lFmTfi8UeMw6Htywul0dr3OZXVH4lYu6eCh4eTuIPDedzHYUxrO/fVL/YMb9qScVeuHY56N4Kp35Jxu4E4y7uVkUFniNSbUynVGCP18TtIIFZxOJ1paWpLsC0HkGyVr6Iswza1p1h3oft3rH9kbn2y7xIeiDBctlzG+KV6olrhNHQBc/17ycyp1p5OtPpb0HL//88Lb7he+ZZinc2I8+qn485Axhj+uNWbMXDhY3x7orjhUK3rOlc5itp+vNDcqcfJLPf/+q63ZbTXzx8/i39WiYYnwjo54mvCzffG9te48uK/Tuv0MAZlhh04H91erohmbmT66gb/ZjV6Wb2O6iyYJ84lEImhqakIkkudGegTXRGIMr2qsxTGCn61MrucBANej+a+LypVOp+hXq2N4djPDugPAw+ti+EJDOtVzMHv02d74FnODn5Dx4Jfx51JiSQzvnP9GFHtDxqWQuXCwVu4qnAHsjGBprf8pNNlsovzEprjXfYXCbu+prN7DkvL0966N4dK3ZIx/RsbVK6IYotPBfXRjrMfm24k8kEVUjkeKwVEsNcLhMBobGxEO57hGnBCax0y4d7UWSfPE4CdkTHg2kjTJvu69KD7S0FPw6hVRhKLJ+/b+8H39XerNxmhHnAsHy1VeDgD44Xv5jyodzvnO4U9+rf8Cyyy+AtGWxb44wWh3ob3epbAA4I8AW1TqwH7aHG+iZzbZ9BrrxMiCUaIwOBwONDc3w+FwmD0UwkTCJkSfv9jP8KtV8QUym7PYteLcZeY8p9bui0/Y9fLcFgbn4p5j/qsAW86lYmQZUVoHa+PGjZgyZQrGjBmDyZMno6WlRfF9jzzyCA4//HCMGjUK8+fPzyocz6LxC6N136Bc2aqzMFwEdgWAqAmzhQNh9b0g81FYnw2pNW+ZOC5lJaXom6WWIuFwGIsWLaIIVoljxq170fIofv9pDDe8H8XILFKD/9nK8MIWPuymVpQ6NvyPzq2DeMDI51VaB+u6667Dtddeiw0bNuCWW27B3Llze7xn8+bNuP322/Huu+9i06ZN2LlzJx5++GHdg7h0hYSHC+RcAcD0l/mOYmXD0i0xfLLH7FEUD6kbW3MR6iV0QTVYBGDuRsdNm7Of9Hau3CMKyx8+y7ODtWvXLqxevRpz5swBAFx44YVobW3Fpk2bkt73zDPPYNasWRg0aBAkScL111+PJ598MquBXFeA9GAnuXQE55VcbmRCHUkCpIURzH1HFq6uoFRxOBxYtmwZpQhLnH06l+sThFEoOlitra0YPHgwbLb45rqSJGH48OHYunVr0vu2bt2KESNGdP3/oYce2uM9iYRCIXg8nqR/BFFIQqEQOjo6ul53bqMSDAa7+iUpba1iOTgNXrKR4b7PY/D7/V2REZ/P19XM0uv1dr32eDyIRqNdr2OxeNGnx+MBYwyxWKzrHohGo12vZVmG1+vteu3zxVs7RyIR+P3xlvbhcLjrtdZjCoXiT5qOjo6u136/vyuFVmzHFAqFcPfdd3f9rWI4pmK8Tvk8pvYww22rxEtTEcVBQTMfd999N6qqqrr+DRs2rJB/niDwwAMP4KqrrgIQ1+OPf/xjAMBtt92G2267DQAw76bbenwuMcvw0S6G6dOnY+nSpQCAE044AW+++SYAYNy4cVi1ahUAoLq6Gl999RUAoKqqCtu3b4fX60VVVRW8Xi+2b9+OqqoqAMBXX32F6upqAMCqVaswbtw4AMCbb76JE044AQCwdOlSTJ8+HQDw2GOPYfbs2ZqP6cc//jHuvvtuAMBVV12FBx54AAAwe/ZsPPbYYwBQdMcUjUZx7733FtUxFeN1yucxael+ThD5QmIK+Y5du3Zh9OjR2LdvH2w2GxhjGDx4MN577z2MHj26631/+tOf8PXXX+PBBx8EALzyyiu466678N577yn+sVAo1DV7AeKzkGHDhgF/3gNUuI0+NoLoQfDKGKLRKCorKxEKhRCLxVBRUdE12y4vL0efJWHsDycXbvSydneTrz9MwuLaEMrKymC32+Hz+VBeXg6bzQav14uKigrYbDZ4PB44HA5YrVZ4PB44nU5IkgSv1wuXywXGGHw+H9xuN6LRKPx+P9xuN2RZRiAQgMvlgizLCAaDcDqdiEQiCIfDcDgcCIfDiEQicDgcCIVCqscUCARgsVjQq1cvdHR0wGq1olevXvD7/bDb7SgrK4Pf76djomMqqmMK2l0Y+Dh5WUSeCHiAn/RDe3s73O6ePoyigwUAp512GubOnYu5c+fimWeewT333IPVq1cnveebb77BSSedhE8++QQDBw7EeeedhxkzZuDGG2/UNDaPxxOfcZCDRRQINt+u+h61lhIXHSahabrNqCEReaIzRfiLX/wCvXr1Mns4hAns7GAYlOfm1UQJo+JgpU0RPvTQQ3jooYcwZswY3HPPPVi8eDEAYN68eXjxxRcBACNHjsRvfvMbTJ06FaNHj0b//v1x3XXXZTXOzy+0oU+BbOBfptCaMIIodmKxGNra2hCLUQ1OqbGzg+GG96NYq2F3C4JI5HADYz1pI1iFoDOCdciDe7Dvur7xARWgIeXeK2wY9qRM+fkS5NzhEl6sSx99+s+3MZyr0tdrxlAJy86mCBZB8MimdsZ9Q+lSYf+VNhz+tIw9Au23PrY38NUB9fcByD6CVUiqK7t9vOg8Gw7N496spw6W0KdcwvbLi+sB+fvjuLiU3PPS1szzif+nYY/K17cxfLKHZsa8EwwG0dDQ0FUPRPDDK1tjeH6z8ZHFaIycK1545BQreveSEBSsnddhLuMap3HxVP7XKd03mkWSsPlS9ToZJa4bq3445db4f6vKJE31OKJw20Qrzh1O7cZz4X8/j2JDu7b3rjtADhZBZMs5y6KYvdz4J+8dHxdHOvjy0RKuOrzw9ryqzLjvihzcomhopXHfmW/+dIIFpw4qMgdr+CHlhnzPbzVEcR46yWrI38oXd0/O/pJkSn0R3bze1tMID3kigpsE2vWdUKe8vByNjY0oLzfGvhDGcN273Y7VD1YYF216YUsMv/+0OO7hxhOtOP/Qwj+e/z3NilkjJHxzsfZnyb0nWPDSjJ7P1cn94+N/6xyxnks/n6D9vP9dxZ/gwsFSauyYDVYJuGRUeu/TaQdGpIT/tl3Gx8Xv2wv4+AIbbtFxcXPh8tE9z9MQgWYauVD3ahTBhI2zHmiJYkeHvu+gWCH/BAIBzJs3zzD7QhjDw191O0GLN6hHgh/5KgZpYQS3ZEjff7AzVlRbywyokHD+oRb0K/Dc4KxhFrwww4bD3MkRtGlDJNRVK1u9H9dYMHNE8nPre8MkTOoXf/8QhznW0pbln7VIEuwaH8MjyjNvE8CFg2Wx9BzGa2fpjzRZJeDJacoO02WjJHxyQc/fmXXxU9lyqQ2T+kmQstg4642zu8/VhD7p37dgUvw8f3+khMdOSz6/N4y3oI0TZ7MQ7E/Y//fGD4pj1kskY7FYUF1drWhfiPywbj/DvqCx6fN5ByNef1yrfJ9uameY+mLxOFezD+1+BkSyNE2jDFgJd3+tFRP7ApP6AW+eY4OsMJaj+wBl1vh4NyVEveYebv49d1x//c/SziV/R2V4jiYyqDLzcZp/FgDFHjV1w/QPzZbhI3dMsuLwKuUT/rOjzT8N2XrbAHB8gpDWzE7vJA2qjPeBeuoMGyyShMqEt548KO7c3TPZohjdKjbCB+1x0zfZWbDiP0Pi06tXLyxYsIB6YBWIpVtiGP+MjJNe0pf2S4wmp5LqrJ3/uoyqRyOY/LyM7zrizhzPRe3XHKHfUjx7ZrdhztbBGtdbwoqZVoztrf0zfzoh+Tl4SC8Jn8y24+ML4rXKStmVX03snqiPckv4+AIbds2x4fujzH+mvvY9K+7SWXIz0h2/Xq9/Tz3Y8ItjLBil0lvK/LMAdO1llSvWNFp22YHRVek/96cTzK/LsuZwJRIdS0lKX7yf2pDDf3X8fX843oKLD94QtxxjxeOn64tkTeyr6+1ccPenMURiDN9/s3hmvkQyHR0dqK+vN8y+EJm54GCKbp3WJe4Huf699Pdg3avJv3vhWwZPBFi9h2HwEzL6/pNf5+qIKmDhydacFlNdq2HhlhIxBpwy2IIjD9Hu4P3s6MzPwTOre46lfmTyzyb1k9C/go/pZ1WZhF8cY9Vcd927DLjwsPjx9C2XsHOODY0npj//EtR9Fy4cLKvVGAcnnYP1m2MtXZv1pmP+WHNFkTj24A/0OThaR640GWLz7bh5Qm7nP9OprbQBG7/PX+rxoa9iKHske+OcRSaXKDBWqxW1tbWG2Rcid77r6BmtWrOX4Y22eJ3VFW8n35Mt+8VdrXtjjSWrko9E7q/NTrtRlvxfozh/RPfxzB2j/dhOMXBlnl5iGlt91qQ4owMqJPz0SEtX54Ge36vuu3DhYBkVwrekuYZaRPbQSVbcd4IF14614LTBhRXD0jOtSQ5gL6v2FhK9rEgrgFRiOm62nXO0OUW3T7TAksHF819tx+g0qVmCyCe9evVCQ0MDpQg54pr/9oxWWSRgxsFI1eObGP61qXsqKK57lcwHs4x38hefak1b3nL56PjPteyOckzfeENQLTRNt2LOaAmXjZKw+FTtE+cV5xZ2kv3hed3nO5cKW0mScMYQ5edXlKn7Llw4WH6/X/nnV2u/KNOGSLCl8bC0OBaSJKHhaCseOtmKt2fmXwxsvh3vnmvF746zYGYO/asePdWqOEsaodCs9dh+2v/OAI1h3t8eZzWkoFI0nt0cwzY/Q9joKSJhGH6/H3V1dWntC5E/Ht2g/FgLKWQDU3fU+PuX8c+Go0y4JpWJJNbV1g7U9qj1ztX+7JkzWsKtaVadX3GwyPz+E9UdO5skoXcvbfbeZpHwz9NteCLNYjIeGN8bOGFA93nRGlhIdwbSOWhXHm5RtS1cOFh2u3K0Rs/gfjgu/bv1RG46eVGhr4dRjDtYeHjSIAt+OdEKa7rQWw6kFrufPljSfJPrZdEpVvx1ioWLxQKF4pnNDNX/khVn5AQf2O121NfXp7UvhHH8enXyffD4xhj++Fm0y9Hyhhn2h5hiliG1ue97Oxn8EYbbV4u9uvfqI/TbQ6dd+7PAIsVrhY7IUF/sKpNUS0guOkzMDMN4hQL+GUMltNQn3+/pnv+ptcNnpwl0VCi4AmN7A0f2kVRtCxdPxLIy5faxqTfj+xnCrIlBnI9T2jFcPUb/YZ47In+n5st64wx+ulvjkJQZida+HlrpX46uNKaZl3hSAAAgAElEQVTTLuGGGmtSfj4VXvqNGc3bOyiCxStlZWWYN29eWvtCGMeda5KdIUkCbvkohqtXxB2vQ/8to89jMt7Ypu1++duXMWzxiXtvnTFEQq90RcE6Sdd8uvPbs63T6uTmAvVeNBKHDShTOOxHTun5Q0eaJfqpjle6aOCDCkXynWVdaraFizObLsyWqs9oyoQmscgu8a2TElJhUwZKGFgppoeuBaVUoFEs+176G3e5wmbHqYLde0X3ezpvhtoBfFyLKw3ahmKbHwjIDOcuk/G/n1M0iyf8fj9qa2spRWgCqQ+WfZn7MfYgysRuhXKSSlH36TrqfG89pqcdfmBqdwH991JaGundDSTXQvxCMrl/vBXE3ittsCuEQ6udPX+Wrhg/NT6a7jworYrsXLkvRIownRdotUi4fWL3EHul6Oyuyd0/SNeFPJv0YCcuzjMLJw6QcKLGtN/wLByxkzMYiaP79vxdqmD7lHe/p1+5hKfPsOL1s7XNtvK90GDJacZF1CoXy/jPVkZb7XBGWVkZGhoaKIJlAp5I9+tfrdI/8bg3TVPRdPxtqgXeuTa8d655K0bvnmzBA1MteGGGFb+e1NMuLz2ze2zjFdonTNbRGHNEiiOReNypDpm4ccCeLD7Fikn94tHBp8/Qdq2tFglnDlV4XjFgyyU2nDKoZ+PtVFIzQC/MiD8/hIhgZcpjTkx4kA+qlPCLY7qHXGYBYvNs+Gy2La2joUe0qWhdSaeHDQa2LFArLk/cT+rPU/JveIaqRArrR1o01Rgc208qyEKDfBDJxaMnDIVqsMzjw13d98GfdDpLALA3BDz1jfq9dO1YC9h8O3443gqnXcLUQRa8M9McJ4sB+NF4K2aNUG4LdF7C3oKVCubtganaH8epZmbqoPhOHDsu59du5rqS8pJREmr6dJ/XQ13an+1KSrrpSCtGuCSsONfWtSggHWtSyo5GHWxIKkQNls/nS/u7KQPjRXr/PM2K4U4JhyWcVIZ4WE8pmvLeuVa8dY4V92VoFKZGHmrP0/bq0kpqFC8Th7klVDuA3x1nQWUWreL1/C0AGF0l6eqNko5crpkecnG+0/ET2naHG3w+H2pqajLaF0Js/jKlp604dbA5jzU9c6s7FCJcmdo1pRZkT1eIyAx1SBjEcTlMroustLYjUkLp3F6jo4lrTR+p6xo8Oa17IGq2hQsHK9Nu9wMrJcTm2zFHwcPMJMipgyw4fYhFMU9baBIvyCAdGyofotBi4y8Jkaifq3TeBYDWy+z45cTslGmRJPx7mr7P9tW43DcTnQZyTp637Fl5nvEz3eZd5GDxQnl5ORobGzPaFyL/hPN0S/zkSEvXPniptF5a+EiOFgdr++U2hH5gg0Mhkp/p45/M7o6UXDvWggodE+brM6ywFwk9XelTqTBADp/MtoPNt+OShG2A1GwLF2feZtN+9Ie54v+tdgB982w3jfLNpg6Mf9FX9TZdkaSv6nuel06n8n9qLJigELkzGqW/kOkQOuvWfndcdtL6b0J4/+GT8xPq7yzez0dxZzrNPNASxdq9lD4sJDabDXV1dbrsC6Gfv39pzuKOn9SktzG9Tegtq+XuHlwpdTmFs1JWXWt1IM7LsFpbib9r3CqGd246KvMz5boMEalFeXqWqNkWLhwsr9er+b3Th1rw5UU2tF5mV93+JlfsFgl3Hpv7KRrmjHdmP6K3vvEqNftkYJCvseF/awtz6c4/tOcYtmXI898ywYKFJ1tx2zHZje/kAoT3ZyjsqWUUFkiIMZa0PcP/fRHFjR/EMOE5GTLVaBUMr9eL6upqXfaF0M8DXxY+ajvSFS+BSIfTLqFB5YFsNBp3ZOmis1AaAH410aIY1Upk5xwbwtfYcPbw7I9LKSsiAof0guLzPjFF/GAGJyq1k4BSajkb1GwLFw5WRUWFrvePyyFUqBel1R4XmtiYbWxvCVaLVLCltWVWqUeqLlOX93KbhHljM+/B9cTpyjdCW0qvrLICqDO1Z1quBKMMQ5+I9/wJyHGL+5Pm7gfQbz+hFGKhqKioQFNTk277QugjX1UY6bYoAbRtfzZJx84VRnDxqOwN1gAN2ZgBFVLOJS+pi6w+v1CM6G66o76xRntkKrGOTc/nMqFmW7hwsHgO4Svdx6n9uArF8f0lUwo4f3KksX/zstE9v29s73iRZiJWi5RUv5YPJvaF4hLebPliP/BdAGgPA997LYoOOVlB6w5QBKtQ2Gw21NbWcm1fioF8WaTl59jS2h61ZfUAcPHI/DpYVyX00pvULz751cvDJ1vx8MlW/Gh8fu36jeMtmH2ohH7lKWnJPuYFC/TQmGHLn9fOii9oU+OdmTb8ZYoFz0437pkiRIrQ4/GYPYS0KPXCSjd7+n6eb+j+JtXqHtff0tUeI3ULHqNIN5PSG3bXiyRJeP1sG768yPjjWrGD4djnZfU3EnnB4/HA7XZzbV+KgXwE0/9TF38Iptun9RQNE02bRcKChNV69xu8Orlvee4HPn+sBfPHWvKyXVoif5lqxbNnmj/RePUs7c5NYi/EqzLsxlI3LL6gTQ13mYQba6yYfZhxOlCzLeafcQAOh8PsIaTlzKESyizJK2FunmDBS1t7FnZWOySkxrwePdU4b9nM2Mddk61JjV2NZGLf9O0rBhqc3XkhzR6T4w6R0HqpDcOeNNYh+upA8v8/szlen5Xv+kEibleam5u5ti/FQD58g3MO1hlNH2rBCGcU32bZaeNXEy04uo+E3cHOze6NSz9QOaV++uioAZvUT0KFDagpYEmQXtRsCxcRLKuV31UOkiThklHdF3jleVacNMiiWB+kFG05xsCVfsV6P/9lijVtzda0oZasutCn4+gMIXGlbRbywXObi/VK8oXVakVNTQ3X9qUYsOR5U5uJKTb0pzpKFqwWCRccZsG14yw4tr+EGQaWAyTexYfpaHppNp/lKQthNHdPtuCVs2z40wn83r9qtoULB4v3EP5Jg7pP0/ED4q+VHpFKc6NcH6VXGbRnHm90dsm/fpwFUwdlluFJA407B2qz7dRC+3zQQVnDguDxeCBJEvf2RWTafAyr9xg7YWhJSdcndkAHctvc2MjyWcbi7WhunWDBYoVNhnnl6L4SvrnYhk0XF97R0hq4L7MgbY8znhAiReh05nHHYgOYd4SEiX2tWJeQ7lGKVuWjXujR02xYsjG+sddwB/+C08qACgmeq2xwlakfU+q2Eg4b4M/SSVFLOaYW2ucDDnrflgROpxOtra3c2xeR+fs641f8pO4gMXeMBTe8HzVkYmKkjXbakXUTZ7PJ1OKi0NgtQCRFRqLE+NVsCxcRLN5385YkCcf1tyTtV6QkAMWfGaiU+wvU+6pQaHGuAODelBUk0zIs307HdWPjm8H20jArOjdNYa1RcC73okGSJLjdbu7ti8jY8mCSwgp9Sz86z4YHplrguSq3mICRdVO/mlhc9rgQDFBYGPA/GRrG8o6abeHiyERsBKjVmTLqfr7gUEnX9gjFRFWZBEeCXQ1F9E9lHzzZqmmjaQCad2nPFi5uuhLA6/WiqqpKSPsiCvkwSUpOW00fCT8ab9U8KUuHUfG2kS6UrD3OhRGunvvV5jswkU+EaDTqcrnMHoJulATwS4UZjRFC8Vxlw1N5fujzzutnW1HtAO6cEMH5IzPvYJ6KfI2+WW+5LXl5t9Fc9nY0qdM7kR9cLhfa29uFtC+iYDf4NjljiITDq/LnuGTqp6QHunuzp1pDGYYo51fNtnDhYDEBHzbXHJEskkOdUNzJfJgB5R+ustw7+IrOlIEWfHuJFXOH7MO1R+j7bDY9Zu44Nr8OrTei/POAzBCUxbsfeIQxBo/HI6R9EQWjU4R11fm1c6nd3c8eVtp21QwSz/jys63CRKuUULMtXDhYPl+WTU5MJHUj4s7TnNgw89eTLBm3lSH04fP5MGzYMPj9frOHYhihaLJD1ecxGRWLZfxzI22pkyudehHRvogAYwyLvjJWp1q2wDGSl8/iYp1XSTE8oR3OGUMteVl9XyjUbAsXDpbb7TZ7CLpJV9wmytYDIuJ2u8EYE1IvqcgHrcqIJ+MO1aMbYrjh/SiCBwt8r3xHodKX0EUx6YVHPt8HrG839jsL7WApoaVXlshRF7O55ggJwxL6cxZqRX4+ULMtXDhY0Sg9TAh1otEoWlpa0uqlvMBlarn0KOt0sHYG4v+9ekUUf/syeS73whaKYuWCml6I3MjH/dbZwT2fvDMz88C11EfeIPDKN7ORJAnr6m1oniV+XbGabeFCJcWQ8lG6Jbk4uUWE3+9HbW1tWr38cJw4Z1xmwIodmR2oe9eSg5ULanohcsPostATBkiG7nyRjlMHW7Bmtg0fna/8gJejUTQlLCq6Ybylxz6zPztafOfATBx2CScOTN+02+gt0vKFmm3hIgFdDCH8qQndxu2WeNF7w1HiPPBFIHnj3p5V4ofo2OfKCBKzxEdU6UuXvNLKcO27mWc/gkTJuYU2es4vW3zGKrSQaaG4I6fszFmtNlw00gK8GcVZ1RL+XGsBA/D0N7QFQz5wKiwK//A8LlwTVYRIEcqymMLdfUW3CP6RsFVC+Bo7Nlxsz7lnC5GMLMtobm5Oq5dbJhRWzjOGWrD87Ph1791L37VWc64A4P2dDDLtKJs1anohcuPMV4xNvfKi9M4UIZtvx6vfs8FqkWCzSHh/lhXThkh4djpFr4zkdoX2RiME2dtRzbZw4WAFAgGzh5AV/Q52pa0qi/dOIvJLIBBAfX09AoEABqWEkI+oMn7vqsQVoYnsnGND+BobLh1twRlDLXjidCuem25FVZmhfx4AsNDgVVqlRKJeCGP5YKfxujx9MB82NF1dzZSBFrx5jg2zD+PisVk0VNgkHCroblZqtoULpYjcCHDTxTZs/L4Y4UzRcblcaGtrg8vlwkfnJ5/zdLPfnx6ZvcTTrQi1SEjqS3bZaAuGOCR8d3l8TEb28sl2z0UiWS+EsUx90fiFA3dN5uJxBMlK9pzQhhCNRkUO4Y9yS+hPva4KgizLWLZsGWRZxjBn8jlPV7/xgyOMl3i6rF25TYJvrg0v1Vnx0EmURjCbRL0QxrFuf36Sedk0BM4HMUrLFxxRz7gQKcJgMGj2EAgBCAaDaGho0KWXfBTOZnoMOOzxrvvXCrSisVjJRi+EOpe9XVwOa0XKXGj+6DTbLBBECmq2hYungNMpaAKWKChOpxMtLS2KepnYrzCz3xMHFDZi+fOVMURpRp0VmfRCZMenexk+3Wv89355kXlpufUpJR5XjK80aSSli6gWTs22cOFgRSI0YyDUiUQiaGpqUtTLo6fGp6Fr0xSmG8X8sYW/Zb6mTgNZkUkvRHZc/KZx0at/ntYdOhp3iHnpwdRyA9JL4Tk2YYJ8xhA+UsVaUNMKFw5WOBw2ewiEAITDYTQ2NnbppbNFAhBfiQIA43onf8bomZGelOMfjzfm9uKkNEU4UvVC5MbGdoYNBm6NM+dwC9ous6HtMr6Kykkvhefx07tt+fJz+NJDJtS0wsWROBwO9TcRJY/D4UBzc3PX/58x1IKZw2OYktDklSdf5OcTrLj5o9yXs6fZ9pJQIVUvRG6Medr42quhDv7ETc+jwlMpaJsjNa1QBIsQhnA4jEWLFiXp5aU6G35xTPfsJ9UZSY1o6aV2QH5u/FN19Pzh4iYVECW9EIVnSCWw5ZLkuXwu7VPyDenFHJafbcXDJ4u1+lpNK1yonHLehBa01NQkui2/O86Sc/PRD87T1m8rHSvPUzYYgyqAWzV2nqcIVnZQDZZxZNvw9sLDJGy73N6jM/f9tXw+SMd98gDpxSTOGGoxpcY1F4SowaKQLKEFh8OBZcuWZdSLJEmo5CLxHef4ARb8d6YVgauTB2WVgKEaZT/r9eJaFl8otOiFSGZTO4O0MNJjI3ItWzspUcj9BXPBO9cG31wbvvz7T0kvhGaESBGGQiGzh0AIQCgUQmNjo6pe1tXHnZnpQ40P/WTzvDh5sAXlNgmhH3Q7WXceZ0VU45d9vi+LP0po1gvRzfNb4o7Vaf+JYk8wd+8o8Rs+uYCjmU8KTrsEWyxMeiF0oaYVLhysdHs/EUQi0WgUzc3NqnoZ7pTA5ttxwgDj5Z1LTVaZVcIxfYElp1ox0i0hStsM5hWteiHixBhLWpSx6+A2a7l0bk9s4Taxn4T/zrTi1bP4TA+SXgi9qGmFiylFZSU1diPUqaysRFNTk2l/f9YIKe3+hFpZM9ve9VprBIvIDrP1Ihqp/WyXbIjhDydYMf6Z7FPUqd958mAu5vSKkF4Ivaj5LlyonUKyhBZCoRAWLFhQNHqRycHKK8Wml3zxrZfhvrXRHg7/H9fmHmIVSeKkF0IvalrhIoIVi1GuhFAnFouhra2t4HrZdLENo58yvtD86ByjYURmzNKLSHztYV3aPq5/HmoWBfKwSC+EXtS0IjFm3i3g8XhQVVWF9vZ2uN1us4ZBEBmRYwxHPStj8SlWnDjQ2KBvTVMEXx7Q9t6mM6y4aCQXQWeiSHhnewynvxyvI/n+SAlPf5P8OBjfG5r1qcTq8204Ng+OG0HwgJoPw4W1pt3uCS0Eg0E0NDQUXC82i4R19XbDnSsAujaOrn8zig7KK2rGLL2IRKKaUp0rIDfnCoBQzhXphdCLmla4cLAIolTRGz9+YhM5WIRxTHuZVswRRL6gFCFBmMjJL8p4b6f2W/DBkyy4bhyfy9wJ8ZAW5rdrOZtvV38TQQiKECnCQCBg9hAIAQgEApg3b15R6UVvOe0/1lMESyvFqBczWDHTikWC7RGXDaQXQi9qWuFiFaHFwoWfR3COxWJBdXV1UelFb/z4o90Me4MMfcvFqW0xi2LUi5FoTV70L5dwymAJVx8hgTHAapGwO8Aw4PHi2sKJ9ELoRU0rlCIkCBO5/C0Z//pa3y24a45NV3E8QSgRjTHYHlF3kkI/sClumn7xm7JiYXwilCIkihkhUoQdHR1mD4EQgI6ODtTX1xeVXh7KIvVCHeC1UYx6MZInNTj2t06wKDpXAPDIKcnaFd3lJ70QelHTChcOltVa/Pl9InesVitqa2uLSi9Ou/7HUpj6IGqiGPViJFe8o76CsDJDEYnTLmHqwG79flWf/OZJ/bIemimQXgi9qGmFUoQEYTJ6V3INdQBtlyWnXhhjWN8OjO0tehyBKBRadPenEyz42dHpHyJBmWFck4y+5RJWX2DDy1tjmLks7rjtv9KG3r1Ij0TxIkSK0O/3mz0EQgD8fj/q6upKXi/bDh7+t97uudGSjfEHXeNa6mvUCekldy46LPMjotwmYfOldqy+IB69Omd49/tFc65IL4Re1LTChYNlt1MhJKGO3W5HfX096QXAqH9HcOi/Zfzxs7hDtXxbPG/4/1ZS/rAT0ktujHACh7rEcpJygfRC6EVNK5QiJAiTOWRJBAfC2X9+w/dtGPN092owWrlFaEEtRXjPZAtuOUZ/PdJ3HQw2C9CPWokQRQ6lCImiwe/3o7a2tuj08vmFubWj2+anZYVKFKtecmGLl2HFjhiiMXXN3Dwhu8fDoEpJSOeK9ELoRU0rXDQaLSsrM3sIhACUlZWhoaGh6PRS7ez5MDq+v4SPdmtznP76ZXJa0B9hcGSxOrHYKFa95MJh/45HOi8bpa4PSSotDZFeCL2oaYVShATBAanpmt8ea8GvP86+norShIQSelaskoYIIjNCpAh9Pp/ZQyAEwOfzoaampiT0MtpdWtGDfFBKetHC0i20ACITpBdCL2pa4cLBKi8vN3sIhACUl5ejsbGxKPXSclF3tt5pBy4ZJeGsanKycqGY9ZINF7yhvYVHRQn22iS9EHpR0wqlCAmCEzrTNw+dZMW14yzwhhncS7LbUJfSO0QqetKD0Xk2WEqsBosg9CJEitDr9Zo9BEIAvF4vqquri1Yve66woeNqG64dF78tlWY+f5vKxS0rBMWul3xSis4V6YXQi5pWuFhFWFFRYfYQCAGoqKhAU1NT0eqlb8rS9nKFNI27rPQefNlS7HohjIX0QuhFTStcTIdtNi78PIJzbDYbamtrS0YvZVYJT07r9rJaL7XhUg3L64k4paaXTKw/QL3S1CC9EHpR0woXDpbH4zF7CIQAeDweuN3uktLL7EMl3DDegj/XWlDtlGCRJDygIU34o/doT8JS1Es6xjZlV8tXSpBeCL2oaYULB8vhcJg9BEIAHA4HmpubS0ovZVYJf51qxf8c2R3J+uE4C76qzzxz+vs6WpJfKnp5aF0UbT6KUOVKqeiFMA41rXARC7VaS3BNMKEbq9WKmpoas4dhOpIkYZSbHqhqlIpern8vBiCG9fU2jOlNKeRsKRW9EMah5rtwEcGikCyhBY/HA0mSSC8ArPQcVaUU9NLQ3J0KXrLRmKjl/9Zy8VgoOKWgF8JYhEgROp1Os4dACIDT6URrayvpBaW3T1w2lIJe7v/C+FTwqYO5eCwUnFLQC2Esalrh4k6ihwWhBUmS4Ha7SS+EJkpNL3d9GkOHnHvq2LzW0+ZSanohckdNK1w4WNTYjdCC1+tFVVUV6YXQRCnq5fN9yd4RYwznLtO3grBE/auS1AuRG2pa4cLBcrlcZg+BEACXy4X29nbSy0GemkaLQzJR7HrxR3q6Qqnz6TY/8J+tpeoy6aPY9UIYj5pWuHCwTNwOkRAIxhg8Hg/p5SDjD8kcnv6uo+d5enRDDJvaS+P8FbteBjyuHpkKZdEOrUhPlyrFrhfCeNS0woWD5fP5zB4CIQA+nw/Dhg0jvRxE7TFw8kvJD+Ar3pZx9YooTv1PaTSdLHa9dChcxtSSkMlL9V/rcYdkOSDBKXa9EMajphUuHCylXagJIhW32w3GGOnlIDEVD+vrlBXEj2+Kf0DpwVyMlKJejl+aHLI6EM78/vNHJHtkCyZZUGkrzSLvUtQLkRtqWuHCwYpGaVsPQp1oNIqWlhbSy0HUHKzKNG2ES2WRFOlFnVHuEhGDBkgvhF7UtMKFg+X3+80eAiEAfr8ftbW1pJeDqKUI0zlY+0NAe7j460xIL+rEGLD5Ei429DAd0guhFzWtcOFgUUiW0ELnRqyklzhDKjP/viLhuXnj+8kzrVtWFv9ehaWuFy3F2jEAh7q6o1g/ruHikWAKpa4XQj9CpAhluUSKQoickGUZzc3NpJeDDKqU8Mtj0t/C5QldHB74MtmhCkaLP4JV6noJa/ChZw2PO1fLvmfFgats6FNeuinDUtcLoR81rXDhYAUCAbOHQAhAIBBAfX096SUBe4Y7eEN7+iiGpQSeo8WslzV70jvIn+2N/+7/VLbROdQJTBsaF9CMaguqykpAFBkoZr0Q+UFNK1w4WNTYjdCCy+VCW1sb6UUHz29RfhCXwqO0mPUy6fn0M+djnov/7uaPij8NbCTFrBciPwjRaJRCsoQWZFnGsmXLSC86+NrDEFFYblgKKwlJL5l5qY6K2xMhvRB6ESJFGAwGzR4CIQDBYBANDQ2kFx1EGXDNip5LiUshRUh6Ueaz2TYcuMqGI/uUgAh0QHoh9KKmFS6mME6n0+whEALgdDrR0tJi9jC4onZg5oekHANea1Pas674H66lrJebVyr357l+nAVH9y3+a58NpawXIjvUfBcuIliRSMTsIRACEIlE0NTURHpJYEa1BUdUpf/97R/HsFthkiUBOOUlGdLCCFbuKs5anVLWy5/WKl/TKO2zl5ZS1guRHWpa4cLBCodV9nMgCMR10tjYSHpJoSKLOPRDX8Xw7nfxh+2ct4uzczXppSdq3f9LGdILoRc1rXCRInQ4HGYPgRAAh8OB5uZms4fBHbkmfIr1mUt66Qk5WOkhvRB6UfNdKIJFCEM4HMaiRYtILylM6kc1NUqQXnpCDlZ6SC+EXtS0woWDRTlvQgtUI6HMA1O7W7ZfP07/LV2sZTmkl54c5iJnPB2kF0IvalqRmJYNq/KEx+NBVVUV2tvbaf8ngsiBl7fGsGo3g18G7k1T4JyJ92dZMWUgF/MtQgXGGCyL1Hs1LTrZinnvJtfXRa6xwVYKPToIogCo+TBcWNRQKGT2EAgBCIVCaGxsJL0ocM5wCxYca4Wc5YLApm+KL4xVbHq5ZWUUw/4VwTX/1bYo4XvDkh2pqjKQc5WBYtMLkX/UtMKFgxWNFucqJsJYotEompubSS8ZmKrSF6uUKDa9/HFtDG1+YPEGbc5w6j6VxZoKNopi0wuRf9S0woWDVVlZafYQCAGorKxEU1MT6SUDF420YEBF/PUHs6yZ31zklLpe7Bbgw/O6NXDSIHK+M1HqeiH0o6YVLhwsCskSWgiFQliwYAHpRYWbj47f1rUlXlNVTHrZE9QffrJbgBMGdGugaXppO9xqFJNeiMIgRIowFivOTtKEscRiMbS1tZFeVPh/R1vB5tsBAH+u5eIWN4Vi0su0l/VvQNyZItx9hQ3tV9lQaaMIViaKSS9EYVDTCq0iJIgiR1qovuz8piMtaKylCAev2BZFENVpqWPzbJAkcqoIIl8IsYqQdi8ntEC73RN6KHW9kHOlj1LXC6EfNa1w4WARBGEu938Rw44OWmbGK3RlCEI8KEVIEEWOlhQhALx1jhWnD6E5F49YF0V0b3PTWYdHEER+ECJFGAgEzB4CIQCBQADz5s0jveSJYksolbJe7jyWC9MuFKWsFyI71LSieBdu3LgRU6ZMwZgxYzB58mS0tLQofvitt97C8ccfj/Hjx6OmpgY333xzViswLBYyBoQ6FosF1dXVpJc8UWwlO8Wil39/HdMdvTprWJFdzAJQLHohCoeaVhRThNOmTcOVV16JuXPn4plnnsEf/vAHrFq1qseH16xZg6qqKowcORLBYBDTp0/HvHnzMHfuXE2DoxQhQeQfrSnCFTOtOGUwPVx444LXZSz9Vp+Htfp8G47tT04WQeQT3SnCXbt2YfXq1ZgzZw4A4MILL0Rrays2bdrU48MTJ07EyJEjAQDl5d5uFqAAACAASURBVOU45phjsGXLFt2D7Ojo0P0ZovTo6OhAfX096SVPFNvjuJT1UmzRyEJQynohskNNKz0crNbWVgwePBg2mw1AfKnv8OHDsXXr1oxf9N133+GZZ57BzJkz074nFArB4/Ek/QOAcDgMIL7ksXPZYyAQ6OqS2tHR0fXa7/d3vd/v9yMSic/OfT4fZDnejM/r9Xa99ng8XfsFeTwexGIxMMbg8XjAGEMsFusaRzQa7XotyzK8Xm/Xa5/PBwCIRCLw+/1d4+58HQqFuk52KBTqys3SMRl3TH6/H7W1tV3vKYZjKsR1ev5Mbf2tgsGgMMek5TpZrVZMmjSpq2xB5GPSi9/v4/6YeLufgsEgamtrIcty0RxTMV4nno5J1RlnKaxevZqNGTMm6WeTJ09mb775Zupbu2hvb2fHHXccu++++9K+hzHG7rjjDob4iuOkf1dccQVjjLGbbrqJ3XTTTYwxxq655hp2xx13MMYYu+iii7q+e8aMGWzhwoWMMcZOPPFE9vTTTzPGGBs/fjx77bXXGGOMDR06lH3wwQeMMcZcLhf74osv2MFUKGttbWXt7e0MAGtvb2etra2s8zR88cUXzOVyMcYY++CDD9jQoUMZY4y99tprbPz48Ywxxp5++ml24oknMsYYW7hwIZsxYwZjjLH77ruPXXTRRV3Hec0119Ax0TFxc0x4OKz6r2bWPKGOqRivk9Ixnb8soun6Jf6rHDuF62MqxutEx1Sax9R5XEqAMcaWLFnCJkyYwCZMmMDuuece5nK5WCQSYYwxFovF2MCBA9nGjRsVv8Dj8bDa2lp25513Kv4+kWAwyNrb27v+dZ7ob775hjHGWCAQYIFAgDHGWEdHBwsGg4wxxvx+f9drn8/HQqFQ1+twOMwYY8zr9XaN2ePxdL1ub29nsix3vY5GoywWi7H29nYWi8VYNBrtOjmyLHe9jkQizOPxdL32er2MMcbC4TDz+XyMMcZCoVDX62AwyPx+f9frjo4OOiaDj2n79u1sxowZXfophmMq1HXS8lBevtkv1DGpXSefz8emTZvG9u/fL/QxXfC6fgfrv1s8XB9T4nXqfG229nbu3MlmzJjB9u7dWzTHVIzXiadjWr9+fUYHS7HI/bTTTsPcuXO7itzvuecerF69ukf0y+fzoa6uDnV1dfj1r3+dOVSmQGeB2O7du9GvXz/dnydKi3A4jMceewxXXnklysrKzB6OUGgtdC+m3knFopdsitzXzLbhmL5UiKWHYtELUTj27NmD/v37py1yV3Sw1q9fj7lz52Lv3r1wu91YvHgxjjrqKADAvHnzMGvWLMyaNQu///3vsWDBAtTU1HR9tr6+Hr/85S81DY5WERJEYShFB6sYkGMM9kf0b/RMDhZB5B81H4aLTu7bt2/H4MGDzRoGIQh+vx/Tp0/H8uXL4XA4zB6OUJSig1UMenliUwxz3o7q/tx3l9swsJIcLD0Ug16IwrJjxw4MGTKE707uFI4ltFBWVoaGhgbSSx4xcb5lOMWgl4j+vs24dYKFnKssKAa9EIVFTStcOFh2e/HMmon8YbfbUV9fT3rJgodO0taqQW/HcJ4pVb3YubDq4lGqeiGyR00rXNyKnf0uCCITPp8PNTU1pJcsmFGtLaKRRcCEW0pVLxYKXmVFqeqFyB41rXDhYJWXl5s9BEIAysvL0djYSHrJgkqbtvcVUwSrGPSSja9E/lV2FINeiMKiphUuHKzOrvEEkQmbzYa6ujrSSxYMqJDwzkz1NGExOVjFoJesHCzysLKiGPRCFBY1rXDhYHW2xyeITHi9XlRXV5NesuRUDRs5H/KY/pYAvCK6XqIxhqtW6F9BSP5VdoiuF6LwqGmFCweroqLC7CEQAlBRUYGmpibSSx4J6X+ec4voetmX3TaE5GBlieh6IQqPmla4iIVSSJbQgs1m69rsmSDUEF0v2RarXzqai3mzcIiuF6LwCJEi7NxNmyAy4fF44Ha7SS+EJkTXSzb+1YQ+wCg3xbCyQXS9EIVHTStcOFjUNZfQgsPhQHNzM+mF0EQp6sXKhUUXk1LUC5EbalrhIjdntWprgkiUNlarNWnfS4LIhOh6eW6L/iWdElVgZY3oeiEKj5rvwsV8h0KyhBY8Hg8kSSK95Jn71hZHpbvoepn/rv7rQE1Gs0d0vRCFR4gUodPpNHsIhAA4nU60traSXnLgN8eq3/I/W1kc/dxLUS93TebCpAtJKeqFyA01rXBxN0rUGY/QgCRJcLvdpJcc+PUkben4rw7E01OeMEOHLGb3UZH0ciDEcCCU+3mePpQLky4kIumF4AM1rXBxN1JjN0ILXq8XVVVVpJcC8P03Zby8NYaqJTIci2V8skc8J0sUvewKMBzymJzU5PXzffrOd59ewP0ncmHOhUUUvRD8IESjUZfLZfYQCAFwuVxob28nvRSAz/cBM5d11wCdsFS8Du+i6OXoZ3ueW6WfZWLFTBt+ehQtFsoFUfRC8IOaVrhwsBgTb3ZMFB7GGDweD+nFBGQG7AuKdd5F0cvOQO7fEYjyfYwiIIpeCH5Q0woXDpbP5zN7CIQA+Hw+DBs2jPRiEv/TLNbqQhH18qP3sjvHxbTFkVmIqBfCXNS0woWD5Xa7zR4CIQButxuMMdKLSXQIliUUUS9/X6d/BadFAmoHUGF2roioF8Jc1LTChYMVjdL0i1AnGo2ipaWF9GISoj3CRdXLFW/r82R/e6wFVmqAlTOi6oUwDzWtcOFg+f1+s4dACIDf70dtbS3pxSSe28IQEqjWR1S9PL5J3zmOiXNJuEZUvRDmoaYVLhwsCskSWujciJX0Yh7//lqcp7kIennwy9yjJQL5vFwjgl4IvhAiRSjLghV3EKYgyzKam5tJLyYiUiJKBL0sXJ9713yKYBmDCHoh+EJNK1w4WIGAAeuUiaInEAigvr6e9JIjOy7Pfo93kZpci6AXIzZnJgfLGETQC8EXalrJ3tIaCDV2I7TgcrnQ1tZm9jCEZ1Bl9g91kWqpRdCLEaezOHaONB8R9ELwhRCNRikkS2hBlmUsW7aM9GIwt0zQbgasCR5BjDF4w/yGT0TQixERQYpgGYMIeiH4QogUYTAYNHsIhAAEg0E0NDSQXgzmjCHan/KXvhWFHIs7VtZFMtxLZG47X4ugF0MiWHyefuEQQS8EX6hphYsUodPpNHsIhAA4nU60tLSYPYyiY/pQfY95bwR4MKEhJgOfxe8i6MWICFYFF1ZcfETQC8EXar4LFxGsSCRi9hAIAYhEImhqaiK9GMi1Yy2QdD7lrSlv93F6OUTQixGO6a06UrxEekTQC8EXalrh4s4Mh8NmD4EQgHA4jMbGRtKLAfz86Oxv/XAMuG1VdwSragmfNSsi6CVXB+vIQ4AKG4/xQ/EQQS8EX6hphYvgssPhMHsIhAA4HA40NzebPYyi4M7jLPhkD8PNB6Mf/cuB3RpLTza2i1H0I4Jeck0RinElxEAEvRB8oea7UASLEIZwOIxFixaRXgygl1XC8nNsGOWOP+GfOsNq8oiMpxT0wun6AiEpBb0QxqKmFS4cLMp5E1qgGon8MSSH3li8IoJevJHcPCTyr4xDBL0QfKGmFYmZuMba4/GgqqoK7e3ttP8TQZjM3HdkLNmobg5mDJXw+rbk97H59nwNq2h5ZWsM5yzLbS/Csb2BdfV07gnCDNR8GC4iWKFQyOwhEAIQCoXQ2NhIeskTfz9JW5ow1bniFV71sj/EcPZrcs7OFUA9sIyEV70Q/KKmFS4crGg0d0NDFD/RaBTNzc2klzxRYZMwc3h2qcIIh096XvVy8ksyXm015nxRDZZx8KoXgl/UtEIpQoIguvCEWVZtF/r0AvZeSakqLUgLjavxmX2ohGfP5GIxOEGUHJQiJIqGUCiEBQsWkF7yiLssuwjWPg4vSTHp5dTBytflsdOKb/WnWRSTXojCIESKMBaj/eAJdWKxGNra2kgvhCaKSS+XjFR2sBz24lv9aRbFpBeiMKhphVKEBEEksdnDMPIp/WlCWkmojWxShA9MteCG93saczrnBGEeQqQIafdyQgu0231hOMxdHFGRYtLLpaO4MNVFTTHphSgMalqhu5YgCKJALPpKf/ppYl/gkF7Zr/AkCMIcKEVIEEQPskljUbpKnTNfkbFcZx+xrZfaMMwpYeZrMl5upQavBMELQqQIA4GA2UMgBCAQCGDevHmklwKw6WLxl/7zqJdsYlCdLcZSY19HHpLraIhEeNQLwTdqWuHCwbJYuBgGwTkWiwXV1dWklwIwqgjqsIpFL316xf+bmmv46HzxnWCeKBa9EIVDTSuUIiQIQhG9aUJKV6kz4xUZb+hIEf59qgXXj4/3uvp8H8PRz3av7qTzTRDmIkSKsKOjw+whEALQ0dGB+vp60guneMJ87dtSDHqJJpzSo/qIH1XkmWLQC1FY1LTChYNltVI3YkIdq9WK2tpa0kuB2H2FvhTUiS/o752VT3jUi14XicMtHosWHvVC8I2aVihFSBCEItnsS0hpq8yc9aqMZW3aTe7aC21Jkas1exje3h7D/LEWuLLc1oggCGMQIkXo9/vNHgIhAH6/H3V1daSXAuEukzB/rLgPcdH1cla11CMtOLGfhIajreRc5QHR9UIUHjWtcOFg2e006yXUsdvtqK+vJ70UkIdPFnelmuh6oexgYRFdL0ThUdMKpQgJgsiIntWElCJMT4wxWBdpT7nOGCph2dniOrgEUexQipAoGvx+P2pra0kvHMMYQ8y8OVsSvOnl8Y36zgsfZ7F04E0vBP8IkSIsKyszewiEAJSVlaGhoYH0wjEzl0V1RWnyCW96adfZxuKQXnkaCKEIb3oh+EdNK5QiJAgiI6e8JOPd73RGXyhV2IMHWqK48QPtmz3vu9KGQ3pRMTtB8IoQKUKfz2f2EAgB8Pl8qKmpIb0UmHdm6u8L9FqrdkciX/CmF4tOX4mcq8LCm14I/lHTChcOVnl5udlDIASgvLwcjY2NpJcCY5H0P+i/91o0DyPRB2960etgEYWFN70Q/KOmFS6WqNhsXAyD4BybzYa6ujqzh0EIAm96ycZRJQoHb3oh+EfNd+EiguX1es0eAiEAXq8X1dXVpBdCE7zphQtjS6SFN70Q/KOmFS7u+YqKCrOHQAhARUUFmpqaSC+EJnjTC6UI+YY3vRD8o6YVLnJzlCIktGCz2VBbW2v2MAhB4E0vehysaUPIGys0vOmF4B8hUoQej8fsIRAC4PF44Ha7SS+C8Nlec1tl8qYXPQ7Wy3X6V24SucGbXgj+UdMKFw6Ww+EwewiEADgcDjQ3N5NeBOGY58xtOMqbXvQ4WOU2imAVGt70QvCPmla4yM1ZrTRbI9SxWq2oqakxexiEIPCklwdaomj83PzeYER6eNILIQZqvgsXESwKyRJa8Hg8kCSJ9GICD57EhanQBS96YYzhxg9i+IYWp3ENL3ohxEGIFKHT6TR7CIQAOJ1OtLa2kl5M4Lpx4kWZedHLDe9T5EoEeNELIQ5qWuHCwZKoAR+hAUmS4Ha7SS+EJnjRyyPrycESAV70QoiDmla4cLCosRuhBa/Xi6qqKtILoQle9KJ3LeXys8WLFhYDvOiFEAchGo26XC6zh0AIgMvlQnt7O+mF0AQvemE6PawzhnJhlksOXvRCiIOaVri4k5leC0SUJIwxeDwe0guhCR70Eo4yyCRXIeBBL4RYqGmFCwfL5/OZPQRCAHw+H4YNG0Z6EYjFJtYf8aCXX62m+itR4EEvhFioaYULB8vtdps9BEIA3G43GGOkF4H419fmORg86GWLV180ZMEkLkxyScKDXgixUNMKF3dzNBo1ewiEAESjUbS0tJBeBMLM9Vg86EWPe8Xm23HHsVTgbhY86IUQCzWtcOFg+f1+s4dACIDf70dtbS3pxSR+cqR+c2Gmg8WDXqiaRxx40AshFmpa4cLBopAsoYXOjVhJL+bwkxouzIVmeNAL1UuLAw96IcRCiBShLJu7KSwhBrIso7m5mfRCaIIHvazdRx6WKPCgF0Is1LTChYMVCATMHgIhAIFAAPX19aQXk8jGVTCzKbbZenlnewybaFs7YTBbL4R4qGnFVqBxZIQauxFacLlcaGtrM3sYhCCYrZettNpfKMzWCyEeQjQapZAsoQVZlrFs2TLSC6EJkfTy+Om0etBsRNILwQdCpAiDwaDZQyAEIBgMoqGhgfRiElmlCA0fhXZE0svlo7kwxSWNSHoh+EBNK1ykCJ1Op9lDIATA6XSipaXF7GGULLYsvKXX2hjWH2A4onfhXS3SC6EH0guhFzXfhYtpUyQSMXsIhABEIhE0NTWRXkxihEvCuN76P7d8mznd3M3Wy44OWkEoEmbrhRAPNa1w4WCFw2Gzh0AIQDgcRmNjI+nFRK4bp99kWEzKE5qpl6DMcOsq2odQJMi+EHpR0woXKUKHw2H2EAgBcDgcaG5uNnsYhE7McrDM1EuIfCvhIPtC6EXNd6EIFiEM4XAYixYtIr0IhsWkZlhm6iVG2UHhIPtC6EVNK1w4WJTzJrRANRJiYpaRMVMvEYpgCQfZF0IvalqhFCEhDA6HA8uWLTN7GIROzEwRmqUXcrDEg+wLoRchUoShUMjsIRACEAqF0NjYSHoxkXOG6TcZVpMcLDP1Mul5alYpGmRfCL2oaYULBysajZo9BEIAotEompubSS8mMrpKwumD9XlMZkWwzNTLLtrOTjjIvhB6UdMKFynCyspKs4dACEBlZSWamprMHgahE7McLNILoQfSC6EXNd+FiwgWhWQJLYRCISxYsID0Ihhr9pqzpE4UvbRfxcU8t+QRRS8EPwiRIozFqCKUUCcWi6GtrY30YjLnH9odkrp2rLoJuXdtDLsChXeyzNLLyl3qf2/qwO5z6C4zc8dGohOyL4Re1LQiMcZM69ji8XhQVVWF9vZ2uN1us4ZBEIQOGGMY+LiM3UHgTydY8POV6g+kbZfZMMRRGo6EtFB9mf/UgRLe3xk3vWy+Pd9DIggiD6j5MFxEsGj3ckILtNs9H0iShE9n2/D+LCuiGif7Zqwk5FUvEgDqQ8ofvOqF4Bc1rXDhYBEEIRZDHBKmDLRAazIlSh5FF4EfUM0VQZQCXDhY5eXlZg+BEIDy8nI0NjaSXjhCawTLjK1jzNDLnZ+oL/HvZZUw0lWAwRC6IPtC6EVNK1w4WIEANY0h1AkEApg3bx7phSN4jmCZoZdff6ztjDx4khVPn2HFZ7MpmsULZF8IvahphYu722Lhws8jOMdisaC6upr0whE2jbVVZjhYPOvFYZdQP7I0iv5FgWe9EHyiphVaRUgQRNb4IwzOR9W3hdnwfRsOrypuh2Lx+hh+8N/kFOH6ehumvChjb0K7HFo1SBDFgRCrCDs6OsweAiEAHR0dqK+vJ71whMMu4cQB6o7Tyl2Fn8cVWi93rulZfzWmt4Q/HG8tyN8ncoPsC6EXNa1w4WBZrWSACHWsVitqa2tJL5yhxXW64p3C7+9WSL0EZIbN3uSfNRwVN69zDu92QN89l7TLK2RfCL2oaYVShARB5MQJS2V8tFvdjBRzauzXq6O4c01ygfuOy20YVBl3rt7aFsO5r0ex/0obysxoCkYQhOEIkSL0+/1mD4EQAL/fj7q6OtILZ/Da4qpQemGM4cF1PVcP9k1YwT1tqAX+q+3kXHEM2RdCL2pa4cLBstuLd2ZLGIfdbkd9fT3phTO0xsCbvinsHm+F0svKXQy7FRo62y3kTIkE2RdCL2paoRQhQRA5Mfl5Gav3qJuRMgvgmWtDlAGVWvs7CMAdH0fx20+Sncc3zrZi+lAu5q8EQeQJShESRYPf70dtbS3phTO0ztAiMWDs0zIci9XbOhhBofSS6lwBQC8uLCuhB7IvhF6ESBGWlZWZPQRCAMrKytDQ0EB64Qym0cViALb44q8veD3uZP1+TRS3r87PCsNC6CVdAqCcixbOhB7IvhB6UdMKpQgJgsiJSc9FsGav/s/de4IFP1sZj/6IusLwsCcjXU5jIqEf0GpBgih2hEgR+nwKFoogUvD5fKipqSG9cEavLB2JTucKAP61yfgC+ELoRcm5unWChZwrASH7QuhFTStcOFi0ezmhBdrtnk+ePiP3xoyXv218mtAsvZBvJSZkXwi9qGlF0cHauHEjpkyZgjFjxmDy5MloaWnJ+CWMMUybNg29e/fOapA2GxUsEOrYbDbU1dWRXjhjmLOnR2HnYOqWb738cpWyU2jj4NgJ/ZB9IfSiphVFU3Ddddfh2muvxYYNG3DLLbdg7ty5Gb/k/vvvx6hRo7IepNfrVX8TUfJ4vV5UV1eTXgRAy/6E+SbfernrU+W05v/UkIclImRfCL2oaaWHJdi1axdWr16NOXPmAAAuvPBCtLa2YtOmTYpf0NLSgqVLl+LWW2/NepAVFRVZf5YoHSoqKtDU1ER6EYBsls7EDF5vY4Zeag4B+pSb71wS+iH7QuhFTSs9HKzW1lYMHjy4K/QlSRKGDx+OrVu39vhwJBLB/Pnz8dBDD2naIDMUCsHj8ST96/weAAgGgwgG4y2RA4EAQqEQgPiO1Z2v/X4/wuFw1+vOz/p8PshyfOm31+vteu3xeBCNRv9/e/cfHVV55gH8e2cmk58zwQj+gAmh2AaBIj8qaizHg7gHrXu0thDtnrWIC9IeZbMlPd09lirYqlCQnPW0XbWwxroqB6JC3dpCu9ajVmJdfmmD2kWRkigIBphJJvN7nv0DMiHJZO68cDP33pnv5xyOQ2Zinss858lz3/ed9009TiaTEBEEAgGICJLJZCqORCKRehyPx1PdaTweTy1mi8Viqb0votFo6nEkEkmdrB2JRBAKhXhNBl9TT08P6urqoGla3lxTvrxPA53NkvUfvR0z9JpcLhcuu+yyVHxGv0/pODVrv0/5mHtGXVMkEkFdXR0SiUTeXFM+vk9Wuqbenz8kGWDnzp1SW1vb72szZ86UV155ZeBLZfny5bJ27VoREfn444+lsrJy0GvOtGLFCsGp7XD6/bnttttERGTZsmWybNkyERFZtGiRrFixQkRE5s+fL+vWrRMRkblz58r69etFROSqq66SzZs3i4jIpEmTZNu2bSIiMmbMGNmxY4eIiHg8Hmlra5PT21FIe3u7+P1+ASB+v1/a29ul95+hra1NPB6PiIjs2LFDxowZIyIi27Ztk0mTJomIyObNm+Wqq64SEZH169fL3LlzRURk3bp1Mn/+/NR1Llq0iNc0DNfk8Xjkrbfeyqtryov36ZfRfn9e/lti0Nf0/vh+ecTQa/L7/eJwOOSpp54y/H1atiOe9hqu/nXM2u9TPuaeQdc0Z84c8Xg88tBDD+XNNeXj+2S1a+q9rnQgIvKrX/1Kpk6dKlOnTpXVq1eLx+ORWCwmIiLJZFIuvPBC2b9//6BvnjVrlowdO1ZqampkzJgxomma1NTUyNGjR9P+sHA4LH6/P/Wn9x/6k08+ERGRUCgkoVBIRER6enokHA6LiEgwGEw97u7ulkgkknocjUZFRKSrqysVcyAQSD32+/0Sj8dTjxOJhCSTSfH7/ZJMJiWRSKT+ceLxeOpxLBaTQCCQetzV1SUiItFoVLq7u0VEJBKJpB6Hw2EJBoOpxz09Pbwmg6/p+PHj0tbWJpFIJG+uKV/ep4GNhojILdv7/u5ar99gTWkx9pri8bi8/fbbqdcb+T4NdQ2HupKWfp/yMfeMuia/3y9tbW0SDAbz5pry8X2y0jUdPHgwY4OVdqPR2bNnY+HChVi4cCGef/55rF69Gjt37sw4Enbw4EFMmzYNJ0+ezDxkdgZuNEqUH36+L4F/3tE3Mdi7cai2/tTQ+01jNfz3ocxrrKZWAXvn2WPD0d7rGsiuG6YSkbqz2mj0iSeewBNPPIHa2lqsXr0azc3NqecWL16Ml156yfAgifQEAoF+66/IOpZOzrwG8/4ZDjx77bnvl6WC+UIqmC+kSi9X0m7iMGHCBLS2tqb9hg0bNqT9+rhx45RGr85UUVFxVt9HhaWiogLt7e3MF4v73Q19jZT/DhdKXUCRQ0M4Yfxu7ZkMV760d6cfibtuND89aGesL6RKL1cssaOaprEwkT5N0+D1epkvFvXTKxxwasAN1X0D415333tVUZT5fXvnONDRLfCl2bj0bAxXvix5I/0Go7+9IbcjdGQs1hdSpZcrltgRjxu7UTa6urpQWVnJfLGof53qxPcvG7rJmFoFTNQ57OGx940b5RqufBlqJRnPH7Q31hdSpbzRqBk8Ho/ZIZANeDwe+P1+5otNaZqmu8O7kS3KcOWLO03VHF1m6I8gE7C+kCq9XLFEg5Xmg4xEg8gZm9WRPSVy+NYNV76kO2fxk3/kpwftjvWFVOnliiUarN4dW4ky6e7uRnV1NfPFxvQaLCOXvwxHvuztFLx4kL+A8xHrC6nSyxVLLHLnHliUDa/Xy7tLm9NtsAz8WcORL38+OniN2HM53n6ChgfrC6nS610sMYLVe8YPUSaJRAL79u1jvtiY3k4NRo5gDUe+VLoHBxjn7+S8wPpCqvRyxRINVu/hjUSZBINB1NXVMV9sTG8E68e7kwgb1LEMR754Biy1KnYC37qEnx7MB6wvpEovVzhFSLbh9Xq5y7LNZdM7vXtccIXOpw2zMRz58psBx/3MulBDkYMNVj5gfSFVtpgijMfjZodANhCPx9Ha2sp8sbF4FttcGdWvGJ0vx8OCxwfs0xVkKuYN1hdSpZcrlmiwQqGQ2SGQDYRCIdTX1zNfbCybbRo0g5a6G50vkTTNYZALsPIG6wup0ssVSzRY3NiNsuHxeNDR0cF8sbFH6/Q/cWfUCJbR+fJZmlqqt3Eq2QfrC6myxUajHJKlbMTjcWzfvp35YmMTRug3JEZOERqZL9NfHPz/+fnV3KIhX7C+kCpbTBGGw2GzQyAbCIfDaGxsZL7kOaMarFzkC88fzB+sL6RKL1c0fEZQwQAAEYZJREFUMXFntUAggMrKSvj9fn6SkKhAaOtjuq+Ru6x19ExHt6B64+C7VavFSUS5o9fDWGIEKxbTL7hEsVgMLS0tzBfKipH5cuWvBzdXLg5e5RXWF1KllyuWaLCi0ajZIZANRKNRNDU1MV8oK0bmy6c9g7/20bcssY0gGYT1hVTp5YolKkR5ebnZIZANlJeXo7W11ewwyCZU8+Vwj6AzDHy5qv/Q1FCrKMZWcAgrn7C+kCq93oUjWGQb0WgUGzZsYL4UACOWhqrmy8SWOKa8EMf7J/r/7P94L4vdUcn2WF9IlV6uWKLB4pw3ZYNrJApH0oCP3qjkS3u3wH+6Vk56vv96q4HH4wDAI1daonSSgVhfSJUt1mBxipCyUV5eju3btzNfCsC3/pj5lPpsqOTLvP8Z+ueF0jz1/cu4/1W+YX0hVbaYIoxEImaHQDYQiUTQ1NTEfLG5vd/UX/r5/MfnPoSlki89A4682fJx37RgIMrjcAoB6wup0ssVSzRYicS5361S/kskEmhtbWW+2NzU83OzOFwlXwZG9M3TI1q7jgn2dPZ/7p0sGkSyH9YXUqWXK5aoFGVlZWaHQDZQVlaGlpYWs8Mgm1DJF22Inu/Vw4MXuNdWnktUZFWsL6RKr3exxAgWh2QpG5FIBCtXrmS+UFayzZdXP03iL8cHf/03f0vi+IBvbbnOiRLuMJqXWF9IlV6uWGIEK5nkx6BJXzKZREdHB/OlQHR0C3wKe03Fk4ITEWBU6anvyTZfft+Rfo3V0h0JnPnMpSOAeV9gc5WvWF9IlV6u8CxCIsq5bM4j/OqFGv50c/b3gEvfTOAX7yWVzwdc/r8JPLw3c6G8f4YD9013wGXUSdREZHu2OIuQp5dTNnjafWF58zNBw44EjvRkdw/4h09ONUmRxKnXZ5svTp2e6eGZDjzwFSebqzzH+kKq9HLFEg0WEVE6P9uXxMXPxvGN3w8+bHmg3rH4kif1X3smvb7p36ayTBKROk4RElHOHQ0JLnxGrRF6+xYnHNBw+dY4Jo4A3qvvmwr889Ekrvp130emVaYJf7w7gRW7hp4iVJ1yJKLCYIspwlAoZHYIZAOhUAiLFy9mvuSBC0rVp9uu2JrA5VtPNWXvn+z/3OYDg+8Ts8mXpAge2sNFzcT6Qur0csUSDZbDYYkwyOIcDgd8Ph/zpYCptGVPfeTApsvXZMyXXZ8LouyvCKwvpE4vVyyxTUNxcbHZIZANFBcXY+XKlWaHQSY6c5zKHxVUutO3XO90Cu5udQBOD4qLh57iK+LCdTqN9YVU6fUulmjVe3p6zA6BbKCnpwf19fXMlzzx9OxTByZPPu/svv+yF/rWcA1skx5t61uPtf/zofOlROfMZi5wLxysL6RKL1csMYLldPJketLndDpRV1fHfMkT3/6SA/6o4MZqB/xRYMYWtUXvh7qHfu6d431jXbVbiiB3pX+d3hYNq69grhUK1hdSpZcrlmiwOEVI2SguLkZjY6PZYZCBlk425pfZwLMEd3+u/z0NOxL42T4uwKJTWF9IlS2mCIPBoNkhkA0Eg0Fcf/31zBca5E9H1HabERE2V9QP6wup0ssVSzRYRUXcZ4b0FRUVob6+nvlC/bzbKXjrqFqDtUrnaBwqPKwvpEovVyzRYLndbrNDIBtwu91YvHgx8yVPJRef3YoFf1S/ufpswHE7vzlk2v7KZFGsL6RKL1cs0WBxSJayEQwGUVdXx3zJU9rAhVRZymanhYue7b+APpFFf3XNRdzCoZCwvpAqW0wR8o6BsuF2u9HY2Mh8yWOJxS7MGZ19YxOISlYN1pm6Y4K3j+l3WK/dZInPAFGOsL6QKluMYHHOm7LBNRL5z6FpeOXvXfjotuyamzkvJ7CnU2267xdc3E5psL6QKlusweruzrChDdFp3d3dmDx5MvOlAIz3avj82y5conMG/K7PBfe8qdYwsb2idFhfSJVerliiwSopKTE7BLKBkpISNDU1MV8KxPklGnZ9w/hpOuH6dkqD9YVU6eWKJRosl4trHUify+XC9ddfz3wpIJVuDd+daEyZ2nV63RX7K0qH9YVU6eWKJRqsrq4us0MgG+jq6oLP52O+FJh40piW6NqX4/BHBT/aqT9JGLyTv2QLDesLqdLLFUs0WKWlpWaHQDZQWlqKlpYW5kuBiRs05OR2AIezPMe3zMUtGgoN6wup0ssVS9ymcUiWsuFyuVBXV2d2GJRjcYNWpXdGgOkvqh0oTYWD9YVU2WKKMBAImB0C2UAgEIDX62W+FBijRrAAIJzQf01sEW/4ChHrC6nSyxVLNFjl5eVmh0A2UF5ejtbWVuZLgampyN10nQbApbpzKeUF1hdSpZcrlmiwnE6n2SGQDTidTkyePJn5UmB+/JXclSl+wrBwsb6QKr1csUSDxSFZykYgEICmacyXAuN2aijjrB0NM9YXUmWLKcKKigqzQyAbqKioQHt7O/OlAOVqc9Cf5HC0jKyF9YVU6eWKJaqJpnHNA+nTNA1er5f5UoByNXX3oxmcHipUrC+kSi9XLNFgcWM3ykZXVxcqKyuZLwXIyAbrpblsomgw1hdSZYuNRj0ej9khkA14PB74/X7mSwEycorwphoH/qmWoxTUH+sLqdLLFUs0WMLTVykLIoJAIMB8KUC/u8GYUafzi0/991++zFEs6o/1hVTp5YolGqzu7m6zQyAb6O7uRnV1NfOlAM0ZY0yp+r9bT30c0aDN4SmPsL6QKr1cscSHn71eb9qvd3V14fDhw0gm7V0OnU4nampq4Ha7zQ7F1rxeL+8u6ZxUlZyaGhxZPPg5bgVR2FhfSNVQvUsvS5SURKL/+RXJZBKrVq3Cli1bTIrIeGVlZdi4cSPGjBljdii2lUgk8MEHH+DSSy/lZoAF7vxiYPJ5Gl4/cna/EH0VGr470YHH3++7eTv2bUuUQzIJ6wupGti7DGSJihIMBnHeeeel/r5q1Sps3boVDQ0NmD59OoqKikyM7tyFw2Hcd999eOCBB/D444/D4bDEzKztBINB1NXVoaOjQ/fOgfJbsRN47SYX7t+ZwE/2nN0I90Wl/f9e5uLC90LG+kKqgsFgxuct0WCdmcyBQABbtmxBQ0MDFixYYGJUxlq6dCmWL1+Ozs5OjBo1yuxwbIkHsRa2397gxI3bTt0x/sMlp25SvAqz7p4B92nTR7Khoj6sL6RKrxG3xFBKPB5PPT5y5AgAYPr06WaFMyx8Ph8A4MSJEyZHYl/xeBytra398oUKx9eq+8rV2itPPW6YnH0JCyzs32HdXOOA7/RZrX83hs1WoWN9IVV6uWKJBisUCqUe9y5ot/u04EC9c/p2X7BvplAohPr6+n75QoVl5QwH/mu2M7WDstupYeKIs///bf+aCw2THfjDjZYYzCcTsb6QKr1csUSDZdeN3Q4ePIjZs2ejsrIS06ZNMzucvOfxeNDR0WHbfKFzt+IrTtz+pf5l61zGniadp+HRq7mgmVhfSJ0tNhq165Cs1+vFgw8+iOeee87sUApCPB7H9u3bbZsvNDz0jo4rdQJyV36NiJPxWF9IlS2mCMPhsNkhZPTII49gyZIlqb+fPHkSI0eOBADMmjUL5eXlZoVWUMLhMBobGy2fL5RbeiNYD06L5SQOsjfWF1KllyuWWHhQUVGh+5oDAcHJqPE/e4QbGO/NXKIXL16M2tparFmzBiNGjEBzczO+/vWvo6qqyviAaEgVFRXYt2+f2WGQxeiNYDXOKMtNIGRrrC+kSq93sUSDFYtlvsP8PCz40uY4ksOwya5TA47c7sLIkqGr9IgRIzB//nw8+eSTWLZsGR577DFs2rTJ+GAoo1gshq1bt+KWW27Juw9B0Nm7brQDfzk+9IdHYrEY84V0sb6QKr3exRINVjSaeWhqZImG/be6hm0EK1Nz1auhoQE333wzJk6ciFGjRuXdNhJ2EI1G0dTUhBtvvJEFkFLWXunAv7cN3WBFo1HmC+lifSFVer2LJRqsbNYw6U3jDbdLL70U48ePx5IlS7BmzRpTYylU5eXlaG1tNTsMshiXQ8OB21x49sMk7ts1uNHiGknKBusLqdKrLZZY5K7XBVrFXXfdhXg8jvnz5wMAenp64PP5UF9fj/feew8+nw/33nuvyVHmr2g0ig0bNtgmXyh3vuDV8NWL0t+EMV8oG6wvpEovVyzRYOnNY1rFq6++irvvvjs1fFxWVoaOjg4cO3YM0WgUHR0dWLVqlclR5q9YLIaWlhbb5AvlVmKINZrMF8oG6wupssUaLKsP4X/66aeYM2cOqqqqsH37drPDKVjl5eX896chJYZYhmX1+kLWwPpCqmwxRRiJRMwOIaPRo0fjgw8+wI4dO7jLr4kikQiamposny9kjqqSwV/zJruZL5QV1hdSpZcrlmiwEomE2SGQDSQSCbS2tjJfKK2ZoxyYMbL/16774z3MF8oK6wup0ssVS0wRlpX1bQTYeyhyvu2m2ztX63BYoqe1pbKyMrS0tJgdBlnYVy90YPfnfXOFL27iMVaUHdYXUnVm75KOJRqsM4fZampqUFZWhvvvvx/33HMPfD5fqumyq1gshvXr18PtdmP06NFmh2NbkUgEq1atwr333ovi4mKzwyELkgEL3VeuXMl8oaywvpAqvSlCSzRYyWTfHafb7cbGjRvxwAMPYPny5SZGZSy3242mpqasjgWi9JLJJDo6OvrlC1EmzBfKFusLqdLLFU1k4D1f7gQCAVRWVsLv98Pr9fZ7LplMorOzEydOnLB9wjscDowePZrNFdEwW/pmAr94r69eyF3ckZuIhkemHgawyAhWOBweFJzD4cCoUaMwatQok6IiqwmHw/jhD3+Ihx9+GCUlaT4yRjRAY2Mj84WywvpCqvTWinPFNRHljeoztqVZN9PeI99EZG+WnSIkIlKVSApc/xkHwOlBIhpelp4i7O3tjh49amYYZBOhUAg/+MEPsHbtWpSWlpodDlnUlNIY/nIc+OyzOPOFssb6Qqp6e5ehxqlMHcE6cOAALrnkErN+PBEREdE5+eijjzB+/PhBXzd1BKuqqgoAcOjQIVRWVpoZCtlAIBBAdXU12tvbOaVMupgvpIL5Qqr8fj/Gjh2b6mUGMrXB6t3VvLKykglNWfN6vcwXyhrzhVQwX0jVUCe08FOERERERAZjg0VERERkMOfKlStXmhqA04nZs2fD5bLEnqdkccwXUsF8IRXMF1KVKWdM/RQhERERUT7iFCERERGRwdhgERERERmMDRYRERGRwXLaYDU0NGDcuHHQNA179+7t99z+/ftx9dVXo7a2FjNnzsS+fftyGRrZwLhx4zBhwgRMmzYN06ZNw6ZNm8wOiSyK9YRUsb5QJkP1LxlrjeTQa6+9Ju3t7VJTUyN79uzp99y1114rzc3NIiLS0tIil19+eS5DIxtIlzdE6bCekCrWF8pkqP4lU63J6QjWNddcA5/PN+jrR48exc6dO3H77bcDAObNm4f29nZ8+OGHuQyPiPIA6wkRGS1d/6JXayyxBqu9vR0XX3xxah8JTdMwduxYHDp0yOTIyGoWLFiAKVOmYNGiRTh27JjZ4ZAFsZ7Q2WJ9IRV6tcYSDRZRNl5//XW8++672L17N0aOHIk77rjD7JCIKE+wvpDRhrXBevrpp1MLBpubm4d8XXV1NQ4fPox4PA4AEBEcOnQIY8eOHc7wyOIG5k9vPhQVFeF73/se3njjDZMjJCtiPaGzwfpCqvRqzbA2WAsWLMDevXuxd+9e3HnnnUO+7oILLsCMGTPwzDPPAABeeOEF+Hw+fPGLXxzO8MjizsyfW2+9FSdPnkw9t3HjRkyfPt3E6MiqWE9IVTAYZH0hZXq1JqdH5XznO9/Byy+/jCNHjuD888+Hx+NJLQb761//ioULF6KzsxNerxfNzc2YMmVKrkIjiztw4ADmzZuHRCIBEcH48ePx6KOPYty4cWaHRhbEekIqWF9Iz1D9S6Zaw7MIiYiIiAzGRe5EREREBvt/tUUvyCFvuXsAAAAASUVORK5CYII=\" />"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"plot(x,soln(0.01))"
]
}
],
"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