Skip to content

Instantly share code, notes, and snippets.

@howardjp
Created September 8, 2019 01:24
Show Gist options
  • Save howardjp/30ea1e09de2e30ee0f4e314bcaa1a5a8 to your computer and use it in GitHub Desktop.
Save howardjp/30ea1e09de2e30ee0f4e314bcaa1a5a8 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Week 2 Discussion\n",
"\n",
"Set the initial parameters"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"k1 = .01;\n",
"ec50 = .5;\n",
"n = 4;\n",
"kminus1 = 1;"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Set the totals values"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"xt = 1;\n",
"yt = 1;\n",
"zt = 1;"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Set the iniital values"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"iv = [0, 0, 0];"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Create the function to integrate"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"f = @(t,x) [k1 *(xt - x(1)) - kminus1 * x(1) * (x(3)^n / (ec50^n + x(3)^n)); \n",
" k1 *(yt - x(2)) - kminus1 * x(2) * (x(1)^n / (ec50^n + x(1)^n)); \n",
" k1 *(zt - x(3)) - kminus1 * x(3) * (x(2)^n / (ec50^n + x(2)^n))];"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Okay, let's take this thing out for a spin..."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"t0 = [0:1:1000];\n",
"[t, xa] = ode45(f, t0, iv);"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Plot me baby one more time!"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAACXBIWXMAABcSAAAXEgFnn9JSAAAAB3RJTUUH4wkIAQglFr+9agAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAwOC1TZXAtMjAxOSAwMTowODozNiUamIoAACAASURBVHic7d1/dBRVnvfxL53u/CAEAqTZCERQIBFBTNB1WDay4Dj80jkaj2edjOcwj4sRMCqgqKuMQ84wK4dFwuq4u47ACsEf+zwQ3R0clpUJyggRsRN0NQ7BMIYzgbTQhiSd0L+rnz9qpzebAImhSN+uer/+Irf7Vn/rpqlPqupW1aBoNCoAAMSbLd4FAAAgQiABABRBIAEAlEAgAQCUQCABAJRAIAEAlEAgAQCUQCABAJRAIAEAlEAgAQCUQCABAJRAIAEAlEAgAQCUQCABAJRAIAEAlEAgAQCUQCABAJRAIAEAlEAgAQCUQCABAJRAIAEAlEAgAQCUQCABAJRAIAEAlGA3dnG1tbXbtm1zuVxOp7OwsLCkpCQzM/PSXSorK99///2jR4/a7fb8/Px777331ltvjb26atWqgwcP9ux1+PBhYysHAMTXoGg0atSyKisrn332WbvdPmPGjFOnTn399dejR4/euXNnVlbWxbo88sgj+/bts9vtM2fO9Pl8tbW1kUikpKRk1apV+hvmzJlz5syZlJSUbh1ra2uNKhsAoALDAqm5uXnu3LkZGRlvvfXWuHHjRGTTpk2vvPLK3XffvX79+gt2effdd5944onrr79+69atI0aMEJE//OEPP/nJT86cObNjx45bbrlF07SpU6fedtttL7/8siFFAgCUZdg5pMrKymAwuHTpUj2NRGTlypU5OTm7d+9ub2+/YJf33ntPRJ577jk9jUTk2muvfeqpp0TknXfeEZFTp05FIpFeD/oBAEzAsECqqakRkQULFnRtnD17diQSOXLkyAW7fPTRR8nJydOnT+/aeMMNN4iI1+sVkYaGBhG56aabjCoSAKAswyY1uFyu9PR0p9PZtXH69Ok7duw4fvz47bff3rNLeXl5UlJSt8Yvv/xSRPS9In3Xyu/3r1+/vrGxUdO0nJycH//4x9dee61RZQMAFGFMIGmaFgwGx4wZ0609LS1NRJqami7Yq+tsOl1LS0t5ebmI3HfffSJy9OhRESkrK7Pb7ddcc01TU9MHH3zw5ptvlpeXz58//9Il5eXl9WtVAMCc6uvr411CL4wJpEgkIiKxs0cxNpst9mqv9u7du27dOrfbvWLFCv3AncfjEZHi4uLVq1c7HA4R2bNnz9NPP/3kk09OmTIlJyfn0gtUf/SvtLy8PAaBQWAEhEFIkL/RjbwwVtO0/nWsq6v70Y9+tHz58kAgsHHjxmXLluntpaWlO3bsKCsr09NIRBYuXLhkyZJgMPj2228bUzQAQA3GBJLD4UhKSqqrq+vWHgqFRKTnVURdbdiw4Z577qmrqyspKXnvvffuvPPO2EuTJ0++5ZZbur1/3rx5IvLVV18ZUjkAQBGGTWpwOp0tLS3dGtva2kRk2rRpF+v18MMPV1VVzZs375lnnrnqqqv68kH6mao+HgYEACQKww7Z5efnB4PBbkdpDx06JCITJky4YJft27dXVVUtWrTopZde6plGjY2Njz/+eEVFRbf2Y8eOyZ+m4eHSLH7QXMcgMALCICQIwwJJvwJp165dsZb29vb9+/dnZ2cXFBSISGtrq8vl+uyzz2JveOONN1JTU/UrYXu6+uqrq6urN27cePbs2a7t27dvF5E77rjDqMoBACow7JDd3LlzJ02aVFFRMW7cuPvuu+/06dM/+9nPfD7fmjVr9De4XK7S0tLhw4fr90VtbW09efLk8OHDH3rooZ5LmzlzZklJybJly55//vni4uIVK1ZMmjTp7NmzO3fu3Lt376xZswoLC42qHACgAsMCyWazbd68efHixWvXrl27dq3eWFpaWlRUdMH3f/755yJy7ty56urqnq+OGjVKRPT72r322mtPPPFE7KX58+evW7fOqLIBAIow8m7fIqJp2oEDB1paWvRnSfS8Mqkfzp49+8knn5w/f97hcNx4443jx4/vSy8uOwCAmITYJBocSOpIiNEHgIGREJtEnhgLAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQAoEEAFACgQQAUAKBBABQgj3eBSSwjkDE0xmK/ZiV7hiSkiQins6Q2xv0dIRkUFREhiTbs9Id40ekdu3b6gu3+sLZQ5NT7f/zN0FTWyAr3RFrcXuD/pAWezW2BE9nqMHjc3sDrb5wqt2WPTR5YtbgscNSvlPxns5QRyAS+zFWSasv3ODxeTqDbm9Qi8q1I1OzM1LGj0jtWmev9LUTkUg0mjRoUGxk/GHN3R5sPOdrag2Etai+UvmjMzLTvsP3UF9I15YLjoyIZGckXzMyLc85uO8Llx4jMzYzxW4bpK+U2xts8JzXf+l226CJWYMnZqVlpTv6vvCwFnV7g+FIVP8x1WHLzkjW/60vvKk14A9rIpKV7rhuVPrErLTvVHy370xs5N3eoKcz1OA5r4+MGDHymWl2vbu+Uk2tgcYWn158dkayPjh9X3i3kZH/PfInWnzNbYGuI58/Zkj/vpO62MjHvpONLX79pfEjUvvxH6rryNuTBmVnJOvFu73Bxha/2xtoPR+WQTIkJWniyMETs9L030vfFx4Ia9E/jU3XkW9s8R870xnbEGVnJN+cM/Q7fSeVQiBdwBfuzr3Hvq095f3q7HlPZ6jbFqpXdtugsZkpN+cM9XSERKTVH2rw+MKRaPbQZBHxhzS3N9jbMkRECq/J/MLdkZlm75Z8dtug8SNSxw5LzR6a3BGIuL2BY2fO222Drhs12G6zNbb4mtoCfSw1M83e9X+piAxJScpKd4wdlpqZZncOcWw+3Nng8fnD2sSstMxUh74ufRyNVLttSnZ6TZO3a2N2RrJe/JCUpFZfuMFzvsHjy0yzT8wa3BEM+0Nag8enZ1U/ZKU7stId2RkpYzNTwlrU3R5MsQ+q+urcdaMGZ6Un+8ORBo+v60j2Rf7ojE9P//cq6JWPzUxJtdvc3qDbG2jw+LLSHeOHp4mIpzOoj1XX7ilJtkBEu8ByexiSkjQxKy0rPXnssJTYhsYf1hZeN/JMR6jVH/J0hPr+m+22Cl+4OyZmpekjY7cNio382MyUscNSPZ3BsBa92MiPHZbS6+dmptn15WelO/RNvKupPdVum5iVlmpPGjHY/rs/tH7XkY8ZkpI0c3ymFo3qIaEnd4PHN35E6vjhaWEt6vYGYkF+QdlDk7v9EfO/Xu3ynew68vmjM8Ka5g9rDR5f/ughB79u60fxdtsgfWQm/9lgX0hrags0tvgaW/yx76T+J2Cv3/n8MUM+PdXRs33ssBT9T9JUu80f1o6d6Wzw+GIj7/YG9C2YiAyesbwf9Q+wQdFoP//zKy4vL6++vv679mpqC5T8v2NfNHfMv27kdX82+KaxGdkZKal2m76Zvpx63N6guz3YEQzrG4Xv9MddXzS2+BvP+UQkM9Vx3Z8NNnb5+v+ZjmB4SLI9e2hy7K8/o3zh7vR0BkUkOyMlOyP5O/3Z3qvYbk1mmn388LRue6WXSd+ON3jO+8PaxJGD9cExauHyp32+xpb/zuzxI1INHPmOQKTB42v49rzdNmjssNSJWWlXYuQbW/zjRqTlDEsxfORdf/Q2tvjsSYP6sc/Rq8YWf4PnfFNbQM+qKzTyqXbbxKzBhn/nGzy+L9wdrb6wntn6HlX/NokDjED6Hw0e35x/qv1RwZ+tmXuNsV9uAIivhAgkDtn9t45A5IdbP1s+K2fV7KvjXQsAWBGz7P7bM3tOTM0eQhoBQLywhyQi4vYGt33S/PunZ8S7EACwLvaQRES2HD5991Tnd53oCQAwEIEkIvJGrfu+/FHxrgIALI1AEtcfvZ7O0PzrRsa7EACwNAJJDn7devdUp7EX1gAAvisCSQ6caP3euKHxrgIArM7qgRTWoh+cOMfxOgCIO6sH0hfNnWOHpTC/DgDizvKB5O7IH5MR7yoAAJYPpDp354SR3+0O/wCAK8HqgeRq8t6cwx4SAMSf1QPpi2YO2QGAEiwdSE1tgbAWZUYDAKjA0oHU2OK7blR6vKsAAIhYPJDc3mDiPnweAEwmDo+fqK2t3bZtm8vlcjqdhYWFJSUlmZmZl+5SWVn5/vvvHz161G635+fn33vvvbfeeuvlV9Lg8U3MYoodAChhoPeQKisri4uLq6qqJk+eHAgEtmzZUlRU5PF4LtHlkUceefbZZ99///3rr78+Jydn3759Dz744AsvvHD5xZw85x83PPXylwMAuHwDGkjNzc1lZWUjR47cs2fP1q1b9+7du3Tp0tOnT2/YsOFiXd599919+/Zdf/31H3744ebNm19//fV333131KhRmzdvPnLkyGXW09ji55AdAChiQAOpsrIyGAwuXbp03LhxesvKlStzcnJ2797d3t5+wS7vvfeeiDz33HMjRozQW6699tqnnnpKRN55553LrOfYmc6JWYMvcyEAAEMMaCDV1NSIyIIFC7o2zp49OxKJXGx356OPPkpOTp4+fXrXxhtuuEFEvF7v5RQT1qJNrQH2kABAEQMaSC6XKz093el0dm3Uw+b48eMX7FJeXv6rX/2qW+OXX34pIr1Ohbg0tzcY1qJjM7kICQCUMHCz7DRNCwaDY8aM6daelpYmIk1NTRfs1XM2XUtLS3l5uYjcd999l/7EvLy8rj/W19d3/dHdHhw7LIXn8gEwq27bQPUNXCBFIhERiZ09irHZbLFXe7V3795169a53e4VK1boB+4uoVsCddN4zjd+BHO+AZhW121gQoTTQF+HpGla/zrW1dWtXbv26NGjw4cP37hx45133nmZlbT6wuNHMOcbAFQxcIHkcDiSkpLq6uq6tYdCIRFJSbnUuZwNGzZs2bIlOTm5pKTkoYceGjrUgCeOu9uDBBIAqGNA95CcTmdLS0u3xra2NhGZNm3axXo9/PDDVVVV8+bNe+aZZ6666iqjivmmI5jnZM43AKhiQGfZ5efnB4PBbqd2Dh06JCITJky4YJft27dXVVUtWrTopZdeMjCNRKTB42MPCQDUMaCBpF+BtGvXrlhLe3v7/v37s7OzCwoK9JbW1laXy/XZZ5/pP77xxhupqan6lbDGamr1D0lJMnyxAID+GdBDdnPnzp00aVJFRcW4cePuu+++06dP/+xnP/P5fGvWrIm9x+VylZaWDh8+/PDhw62trSdPnhw+fPhDDz3Uc2kzZ84sKSnpdzGezhBzvgFAHQMaSDabbfPmzYsXL167du3atWv1xtLS0qKiogu+//PPPxeRc+fOVVdX93x11KhRl1OMpzOUncFVsQCgikHRaHSAP1LTtAMHDrS0tOjPkuh5ZZIh8vLyLnEdkj+spT39gW/97FS7pZ8IBcAiLr1JVEQcnodks9nmzJkz8J/blbs9mJlmJ40AQB0W3SI3tfHgCQBQi0UDKaxFOYEEAEqxaCC524PZGcnxrgIA8D8sGkiezhAXIQGAUqwbSJxDAgClWDSQTrUFhqXGYYYhAOBiLBpIns5Q9lDOIQGAQiwaSK2+MPcNAgClWDSQ3N4As+wAQCkWDSR/WMtMY1IDACjEooHU6gszyw4AlGLdQOKQHQAoxYqB5PYGRYQLYwFAKVYMJH9IGzuMG9kBgFqsGEiezlDWEE4gAYBarBhIHcFwZiqBBABqsWIguduD8S4BANCdFQPJH9a4bxAAqMaKgdTUGuDh5QCgGitulzuDEa6KBQDVWDGQ3N5gejIXIQGAWqwYSCKS6rDoigOAsqy4XeZGdgCgIIsGEg9DAgDVWDGQRIQ7qwKAaqwYSP5whDurAoBqrBhIns5QdgY3VwUAtVgxkMJalHNIAKAaKwZSU2tgbCZ7SACgFssFUkcgkuqwsYcEAKqxXCB5OkNchAQACrJcIAEA1GS5QPKHNW71DQAKstym2e0NMOcbABRkuUACAKjJcoHk6QxlptnjXQUAoDvLBVJHIEIgAYCCLBdIAAA1WS6QmloD3OobABRkuUAKa9EUpn0DgHost2n+9jyTGgBARZYLJCY1AICarBhI3KkBABRkuU2zpzOUPZRJDQCgHMsFUliLxrsEAMAFWC6Q3N7A2GGp8a4CANCd5QKJ55cDgJosF0ju9iDnkABAQZYLJJ6HBABqYtMMAFCCtQLJH9Y4gQQAarJWILnbg9xZFQDUZPBNdGpra7dt2+ZyuZxOZ2FhYUlJSWZmZl86dnR0zJ0794UXXpg5c2bX9lWrVh08eLDn+w8fPty/Cu1J7CEBgIqM3EOqrKwsLi6uqqqaPHlyIBDYsmVLUVGRx+PpY99vv/02FAp1a6+pqfF6vcEe+lchh+wAQFmG7SE1NzeXlZWNHDnyrbfeGjdunIhs2rTplVde2bBhw/r16y/R8eTJk7/97W9feumlni9pmvbNN9/cdtttL7/8siFFur2B7IwUQxYFADCWYYFUWVkZDAaXLl2qp5GIrFy58je/+c3u3btXr149dOjQC/a69dZbz5w5c7Flnjp1KhKJ9PGgXx+xhwQAajLskF1NTY2ILFiwoGvj7NmzI5HIkSNHLtZrzZo1mzZt2rRp04wZM3q+2tDQICI33XSTUUW2+sJGLQoAYCzD9pBcLld6errT6ezaOH369B07dhw/fvz222+/YK9Ye3V1dc95Cu3t7SLi9/vXr1/f2NioaVpOTs6Pf/zja6+9tn9FtvrCWemO/vUFAFxRxgSSpmnBYHDMmDHd2tPS0kSkqampf4s9evSoiJSVldnt9muuuaapqemDDz548803y8vL58+f32v3vLy8rj/W19d7OkNDUpL6VwwAJJZu20D1GRNIkUhERGJnj2JsNlvs1X7QZ+gVFxevXr3a4XCIyJ49e55++uknn3xyypQpOTk5l+5eX1/fraUj0M9KACDhdN0GJkQ4GTntW9M0A5cmIqWlpTt27CgrK9PTSEQWLly4ZMmSYDD49ttv92+ZXBgLAGoyJpAcDkdSUlJdXV23dv26opSUfs60njx58i233NKtcd68eSLy1Vdf9WOBp9oCKdxZFQCUZNjW2el0er3ebo1tbW0iMm3aNKM+RUT0M1X9OwwY1qKcQwIANRkWSPn5+cFgsNtpm0OHDonIhAkT+rHAxsbGxx9/vKKiolv7sWPHRKTfFycxyw4A1GRYIOlXIO3atSvW0t7evn///uzs7IKCAhFpbW11uVyfffZZHxd49dVXV1dXb9y48ezZs13bt2/fLiJ33HFHP4p0e7m5KgAoyrBAmjt37qRJkyoqKl5//fVQKHTy5MlHH33U5/OtWLFCf4PL5br//vuXLFnS18pstmXLlvn9/uLi4nfffbe+vv7gwYPLly/fu3fvrFmzCgsL+1GkP6SlOjiHBAAqMuzCWJvNtnnz5sWLF69du3bt2rV6Y2lpaVFRUb+X+ZOf/OTMmTOvvfbaE088EWucP3/+unXrLrdcAIBiBkWjUQMXp2nagQMHWlpa7HZ7fn5+zyuT+uHs2bOffPLJ+fPnHQ7HjTfeOH78+L70ysvL63kd0qR1H/1HSf7ErLTLrwoAEsgFN4mqMfh5SDabbc6cOcYu0+l0Lly40JBFhbUoN1cFADVxQgUAoAQCCQCgBGsFUlNrYGwmD+gDABVZK5A4hwQAyrJWIAEAlGWhQPJ0hrhvEAAoy0KB1BGIcGdVAFCWhQIJAKAyCwWSP6yl8jAkAFCVhTbQbm8gK51bfQOAoiwUSCLCnG8AUJaFAqnVF453CQCAi7JWINmT2EMCAEVZKJBEhMfFAoCyLBRIbm+Qc0gAoCwLBZI/pMW7BADARVkokCLRKLcOAgBlWSiQ3O3B9GRuHQQAirJQIIW1aLxLAABclIUCSZhlBwAKs1Agub3BVIeF1hcAEouFNtD+kMYeEgAoy0KBJCLsIQGAsthAAwCUYKFAamrzjx2WGu8qAAAXZqFACmtRbh0EAMqyUCABAFRmoUDyhzQmNQCAsiy0gXZ7g0z7BgBlWSiQAAAqI5AAAEqwSiA1tQXGDkuJdxUAgIuySiCFI1F7EnO+AUBdVgkkAIDirBJITW3+eJcAALgUqwRSWItmptnjXQUA4KIsE0iRaGaqI95VAAAuyiqB5PYG410CAOBSrBJIYS0a7xIAAJdilUDyh7SxmVyHBADqskogeTpDPHsCAFRmlUDikB0AKM4qgRQI8+wJAFCaVbbRns5QVjrTvgFAXVYJpLAWTRrEOSQAUJdVAklEOGQHACqz0DaaQ3YAoDKrBJLbG2TaNwCozCqB5A9pqXarrCwAJCILbaM5ZAcAKrNKIIU1nhgLAEqLQyDV1tY+9thjM2fOvOuuuzZs2NDa2trHjh0dHTNnzqyuru7Hh7q9gVR7Uj86AgAGxkAHUmVlZXFxcVVV1eTJkwOBwJYtW4qKijweTx/7fvvtt6FQqB+fG9aiHLIDAJUNaCA1NzeXlZWNHDlyz549W7du3bt379KlS0+fPr1hw4ZLdzx58uTWrVvLy8v7/dH+kMYsOwBQ2YA+1buysjIYDC5dunTcuHF6y8qVK3/zm9/s3r179erVQ4cOvWCvW2+99cyZM5f50TygDwAUN6B7SDU1NSKyYMGCro2zZ8+ORCJHjhy5WK81a9Zs2rRp06ZNM2bM6N/n+sOacKcGAFDbgG6jXS5Xenq60+ns2jh9+nQROX78+MV63X777QsXLly4cGFOTk7/PtfdHhSmfQOA2gYukDRNCwaDo0aN6taelpYmIk1NTVe6AM4hAYDKBu4cUiQSEZHY2aMYm80We9VYeXl5+j9CaSPl+78wfPkAoLLYNjBRDOikBhHRNG3APqu+vl7/R4PHN3n94QH7XABQQWwbKAkSTgN3yM7hcCQlJdXV1XVr168rSklJuXIf7fYGsjOSr9zyAQCXb0AnNTidTq/X262xra1NRKZNm3blPpf7BgGA+gY0kPLz84PBYNe9SBE5dOiQiEyYMGEgKwEAqGZAA0m/AmnXrl2xlvb29v3792dnZxcUFOgtra2tLpfrs88+M/Bz3e1B5nwDgOIGNJDmzp07adKkioqK119/PRQKnTx58tFHH/X5fCtWrIi9x+Vy3X///UuWLDHwc/1hbUjyQE/fAAB8JwMaSDabbfPmzRMmTFi7du3UqVPnzp17+PDh0tLSoqKigSwDAKCgQdFodIA/UtO0AwcOtLS02O32/Pz8nlcmGSIvLy92smrLx6f/76dn9i3JvxIfBADq67pJVFYcDmTZbLY5c+YM5Ce624NM+wYAxVnifqNhLcp9gwBAcZYIJACA+qwSSGMzr+CdIAAAl88SgdTmDycN4pAdACjNEoHU6gvHuwQAQC8sEUjCITsAUJ5VAolZdgCgOEsEkqczFO8SAAC9sEQgdQQiXBgLAIqzRCCJSKrDKmsKAAmKzTQAQAmWCCR/OJJqT4p3FQCAS7FEILm93FwVAFRniUACAKiPQAIAKIFAAgAowRKB5G4PZg/lHBIAKM0SgeQPa6l2S6wpACQuNtMAACWYP5Dc3mC8SwAA9M78geQPaVyEBADqM38gNbX5uZEdAKjP/FvqsBaNdwkAgN6ZP5D8IY2n8wGA+swfSJ7OUHYGzy8HANWZP5DCWpQ9JABQn/kDCQCQECwRSGMzOWQHAKozfyC1+sIcsgMA9VkikOJdAgCgd+YPJBHhTg0AoD5LBFIKt/oGAOWZf0vd5ueQHQAkAPMHUqsvnJXuiHcVAIBemD+QRGRISlK8SwAA9MISgQQAUJ/5A6nVF85Ms8e7CgBALwgkAIASzB9IAICEQCABAJRAIAEAlGD+QPJ0BrPSuXUQAKjO/IHUEYwMSeY6JABQnfkDCQCQEAgkAIASzB9Ins4Q97IDAPWZP5A6AhHuZQcA6jN/IAEAEgKBBABQgskDqaktMHZYSryrAAD0zuSBFI5E7UmD4l0FAKB3Jg8krooFgERh8HMZamtrt23b5nK5nE5nYWFhSUlJZmbm5XRZtWrVwYMHe/Y6fPhwX+rhvkEAkCiM3EOqrKwsLi6uqqqaPHlyIBDYsmVLUVGRx+O5nC41NTVerzfYg4FlAwBUYNgeUnNzc1lZ2ciRI996661x48aJyKZNm1555ZUNGzasX7++f100Tfvmm29uu+22l19+2ag6AQBqMmwPqbKyMhgMLl26VI8WEVm5cmVOTs7u3bvb29v71+XUqVORSKTXg36XwFWxAJAoDAukmpoaEVmwYEHXxtmzZ0cikSNHjvSvS0NDg4jcdNNN/a6K+wYBQKIwLJBcLld6errT6ezaOH36dBE5fvx4/7ro+0l+v3/9+vXLli1bsmTJL37xiz/84Q9G1QwAUIcx55A0TQsGg2PGjOnWnpaWJiJNTU3963L06FERKSsrs9vt11xzTVNT0wcffPDmm2+Wl5fPnz+/16ry8vLaxs3yDx1TveYtEamvr+/PugFAYsrLy4t3Cd+NMYEUiUREJHYqKMZms8Ve7UcXfbpdcXHx6tWrHQ6HiOzZs+fpp59+8sknp0yZkpOTc+mq6uvrf7GvMRDR1s4v689aAUAi6/pXeEKEk5HTvjVNM7ZLaWnpjh07ysrK9DQSkYULFy5ZsiQYDL799tt9WX5YiyYN4k4NAJAAjAkkh8ORlJRUV1fXrT0UColISsoF7ibXly6TJ0++5ZZbur1h3rx5IvLVV18ZUjkAQBGG7SE5nU6v19utsa2tTUSmTZtmVBcR0U87XfAwYE9t/nBmmsF3owAAXAmGBVJ+fn4wGOw2ceDQoUMiMmHChH50aWxsfPzxxysqKrr1OnbsmIj08eKkVh+BBACJwbBA0i8n2rVrV6ylvb19//792dnZBQUFItLa2upyuT777LM+drn66qurq6s3btx49uzZrh+0fft2EbnjjjuMqhwAoALDAmnuG/gC+gAADq9JREFU3LmTJk2qqKh4/fXXQ6HQyZMnH330UZ/Pt2LFCv0NLpfr/vvvX7JkSR+72Gy2ZcuW+f3+4uLid999t76+/uDBg8uXL9+7d++sWbMKCwuNqhwAoIJB0WjUqGU1NzcvXrz4xIkTsZbS0tLHHntM//dvf/vb0tLS4cOHd71R96W7iMiGDRtee+21rmeM5s+fv27dusGDB1+6mLy8vPr6+gf+9fd/NSHz//z5VZe5agCQ0PRNYryr6IWRgSQimqYdOHCgpaXFbrfn5+f3vMyoH13Onj37ySefnD9/3uFw3HjjjePHj+9LJfro/3Drfy35i9F3Xp/Vj3UBANNIiEAy+IS/zWabM2eOsV2cTufChQv7Vw83VwWARGHyJ8YCABIFgQQAUILJA8nTGcxM5fETAJAATB5IHcEIF8YCQEIweSABABKFyQPJH/rONyAHAMSFyQPJ7Q3ak3j8BAAkAJMHkohkZyTHuwQAQO/MH0h2G3tIAJAAzBxITW2BeJcAAOgrMwdSOGLkbfoAAFeUmQOpIxjhBBIAJAozB5KnM5jqMPMKAoCZsL0GACjB5IHEfYMAIFGYOZA6AhHurAoAicLMgeTpDHGbBgBIFGYOJOE2DQCQOEweSNymAQAShZkDyR/SCCQASBRmDiRPZyh7KIfsACAxmDmQwlo0aRB7SACQGMwcSACABEIgAQCUYOZAikSjTGoAgERh5kBytweZ1AAAicLMgRTW2EMCgIRh5kACACQQAgkAoAQzB5LbG+RedgCQKMwcSP6QxhNjASBRsL0GACiBQAIAKMHMgeTpDGalcw4JABKDmQOpIxgZkpwU7yoAAH1i5kACACQQMwdSRyAS7xIAAH1l5kDydIbiXQIAoK/MHEgikpXuiHcJAIA+MXkgDUlhUgMAJAbTBlI4bXi8SwAAfAemDaSoeVcNAEzJzFttjtcBQAIxbSBFkocwowEAEohpA0mzJ/O4WABIIKYNJBEZOyw13iUAAPrKzIEEAEggpg2kSMqw7KHc6hsAEoZpAyma5Ei1m3btAMB82GQDAJRAIAEAlEAgAQCUYI93AX1SW1u7bds2l8vldDoLCwtLSkoyMzMv3SWcMjQ7g0kNAJAwEmAPqbKysri4uKqqavLkyYFAYMuWLUVFRR6P59K9ojZHiuUnNeTl5cW7hPhjEBgBYRAShOqb7Obm5rKyspEjR+7Zs2fr1q179+5dunTp6dOnN2zYEO/SAABGUj2QKisrg8Hg0qVLx40bp7esXLkyJydn9+7d7e3t8a0NAGAg1QOppqZGRBYsWNC1cfbs2ZFI5MiRI5foGE7J4BwSACQQ1QPJ5XKlp6c7nc6ujdOnTxeR48ePX6Jj1OZIdai+dgCAGKVn2WmaFgwGx4wZ0609LS1NRJqami7Rd8SJ99Yu/rd1Aasf1uNcrjAIjICIMAiJQOlAikQiIhI7exRjs9lir17M1zUHrlxhAADDJcBBLU3T4l0CAOCKUzqQHA5HUlJSXV1dt/ZQKCQiKSkp8SgKAHBFKB1IIuJ0Or1eb7fGtrY2EZk2bVo8KgIAXBGqB1J+fn4wGKyvr+/aeOjQIRGZMGFCnIoCABhP9UDSr0DatWtXrKW9vX3//v3Z2dkFBQXxqwsAYDDVA2nu3LmTJk2qqKh4/fXXQ6HQyZMnH330UZ/Pt2LFiniXBgAw0qBoNBrvGnrR3Ny8ePHiEydOxFpKS0sfe+yxOJYEADBcAgSSiGiaduDAgZaWFrvdnp+f3/PKJABAokuMQAIAmJ7q55AAABZBIAEAlEAgAQCUoPTNVfuntrZ227ZtLpfL6XQWFhaWlJRkZmbGuygjVVZWvv/++0ePHtWneNx777233npr1zf0ZQTMMUrnz58vLi6+7rrr1q9f3+0l0w/C2bNn//mf//mTTz759ttvb7755kWLFt18883d3tPrCib0CIhIZWXl3r17P//88xEjRkyfPv3BBx8cP3581zeY8mvQ0dExd+7cF154YebMmT1fNeSXHq8xMdukhsrKymeffdZut8+YMePUqVNff/316NGjd+7cmZWVFe/SjPHII4/s27fPbrfPnDnT5/PV1tZGIpGSkpJVq1bpb+jLCJhmlJ588slf//rXP/jBD15++eWu7aYfhN///vcPPPDAuXPnpkyZMnjwYP1r8Pd///d33XVX7D29rmBCj4CIPPzww1VVVcnJyTNmzDhz5syxY8dSU1NfffXV733ve/obzPo12L59+/PPP//qq6/+1V/9VbeXDPmlx3NMoiZy+vTpqVOn/sVf/EVjY6PeUl5enpub+9RTT8W3MKPs3r07Nzf37rvv/vbbb/WWEydOFBYW5ubmfvzxx9G+jYBpRkkfjdzc3NLS0q7tph+EcDj8/e9/f+rUqYcPH9Zbvvjii4KCgoKCgs7OTr2l1xVM6BGIRqP/9m//lpube//993u9Xr3ld7/73eTJk2fNmhWJRKIm/Ro0NjZu2bJl2rRpubm5H3zwQbdXDfmlx3dMTBVIv/zlL3Nzc7dv39618fvf//7kyZPb2triVZWBHn300dzc3Jqamq6Nv/71r3Nzc//2b/822rcRMMconT59uqCg4MEHH+wZSKYfBP03/stf/rJr4z/8wz/k5ub+x3/8h/5jryuY0CMQjUZXrlyZm5v7u9/9rmvjQw89lJub+1//9V9RM34N9D89Y3oGkiG/9PiOiakmNdTU1Mifbn8XM3v27EgkcuTIkTgVZaSPPvooOTlZf4J7zA033CAi+j3R+zICJhglTdNWrFiRkZHx85//vOerph+E999/X0Tuueeero3Lly+vr6+fP3++/mOvK5jQIxCTlJTU9cfU1FT50+NpzPc1WLNmzaZNmzZt2jRjxowLvsGQX3p8x8RUgeRyudLT051OZ9dGffN9/PjxOBVlpPLy8l/96lfdGr/88ksR0U859mUETDBKL7/88qeffrpx48Zhw4b1fNX0g3D06NH09PTRo0drmnb06NEPP/ywvb2923t6XcGEHgER0Sfy/Pu//3us5eTJkx9++GFWVlZ+fr6Y8Wtw++23L1y4cOHChTk5ORd8gyG/9PiOiXlm2WmaFgwGx4wZ0609LS1NRJqamuJRlMG6zaYTkZaWlvLychG57777+jICJhgll8v1j//4j0uXLr355pvPnz/f7VUrDMI333wzadKk/fv3r1mz5syZM3rjD37wg5///OcjRoyQPqxgoo+AiBQVFTU0NGzZsuXYsWNTpkwJhUL79u1LT0//p3/6J5vNZoWvQTeG/NLjPibmCaRIJCIiPW9zZ7PZYq+azN69e9etW+d2u1esWHHDDTfoRyouPQKJPkrt7e1PPvnklClTli9ffsE39GUFE3oQQqFQJBJpbW1dtmzZzJkzV61apWnazp079+3bd/r06V27dtlstl5XMKFHIKagoCArK+vYsWPHjh3TWwoLC/W9B9N/DXoy5Jce9zEx1SE7EdE0Ld4lDIS6urof/ehHy5cvDwQCGzduXLZsWeylvoxA4o7ST3/607a2thdffFH/H3Ix5h4EEXG73X/zN3/z2muv3XXXXUVFRW+++eZtt91WV1e3c+fO2Ht6XcGEHoF33nmntLR09OjRb7/9dn19/RdffLF69eoPP/zwr//6r2MHME3/NejJkF96HMfEPIHkcDiSkpLq6uq6tev7DSkpKfEo6orYsGHDPffcU1dXV1JS8t5779155516e19GIKFH6cCBA//5n/85f/78P/7xj9XV1dXV1R9//LGItLS0VFdX62dcTT8I+mn8pKSkbo8Ee+CBB0REH5BeVzChR0BENE37u7/7u4yMjK1bt06ZMkVEHA7HokWLHn/88T/+8Y8VFRWm/xr0ZMgvPe5jYp5DdiLidDpbWlq6Nba1tYnItGnT4lGR8fSLAefNm/fMM89cddVV3V7tywgk7ii1traKSGVlZWVlZdf2mpqaBx54ID09vba2Vsw+CDabLT093Wazdds6TJ06VUQ6Ozv1H3tdwcQdARH5/PPPvV7vvHnzhg4d2rV94cKFzz//vD7Nx9xfgwsy5Jce3zExzx6SiOTn5weDwfr6+q6Nhw4dEpEJEybEqSgjbd++vaqqatGiRS+99FLPNJK+jUDijtL06dP//n/Tp31PmTIl9m8x+yCIyJ//+Z/7fL5uB/S/+uorEUlPT9d/7HUFE3oE9NUMh8Pd2vVJLvrhXNN/DXoy5Jce3zExVSDpc+d37doVa2lvb9+/f392dnZBQUH86jLMG2+8kZqa+tRTT13sDX0ZgcQdpZycnLv+tx/+8IciMnr06Lvuuit26NLcgyAif/mXfxkOh//lX/6la+O//uu/6i/pP/a6ggk9AhMnThw5cuSBAweam5u7tuun0K6//nqxwNegJ0N+6XEekyt95e1AikQid9xxR25u7o4dO4LBYGNj46JFi3Jzc99+++14l2aAc+fO5ebmfu973/s/F/Lqq69G+zYCZhqlzs7OnndqMP0g+P3+OXPm5Obmvvjii59++unHH3/83HPP5ebm3nvvvfpdc6J9WMGEHoFoNLpz587c3NxZs2bt2rXr2LFjhw4dWrt2rd6i30zIxF+D1atXX/BODYb80uM7Jma7uWpzc/PixYtPnDgRayktLX3sscfiWJJRPvzwwwcffPBir9599936Ha/7MgKmGaXz588XFBT0vLmq6QehsbHxoYceOnnyZKzllltuefHFF/XrkHS9rmBCj4CIVFRUlJeX+3y+WMv111//0ksvxa4bNevX4Kc//enOnTsveHNVQ37pcRwTswWSiGiaduDAgZaWFv3pDD3n1JteX0bA9KNk+kHQNO3QoUNnzpyx2WzXXnvtjTfeeMH3XHoFE3oERKSjo6O6urq9vd3hcFx99dU9jymZ/mvQkyG/9HiNiQkDCQCQiEw1qQEAkLgIJACAEggkAIASCCQAgBIIJACAEggkAIASCCQAgBIIJACAEggkAIASCCQAgBIIJACAEggkAIASCCQAgBIIJACAEggkAIAS/j8mo5cUKhF+8wAAAABJRU5ErkJggg==\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Warning: MATLAB has disabled some advanced graphics rendering features by\n",
"switching to software OpenGL. For more information, click <a\n",
"href=\"matlab:opengl('problems')\">here</a>.\n"
]
}
],
"source": [
"plot(t, xa(:,3))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Try with a different initial value set."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAACXBIWXMAABcSAAAXEgFnn9JSAAAAB3RJTUUH4wkIAQgmj7bs0AAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAwOC1TZXAtMjAxOSAwMTowODozOMKitY0AACAASURBVHic7d1/VFT1/u/xj/MD0PG30JcwsqU5pKgBy6xIS/th6qlFmCtzeS99TVGKTOyW55uWeqQf60Ta0bJWgaUc69u5MHZMl4cy8BALUhvQjlJgtY5cDcgUYfg9wwz3j33PXM6AgAp7f2b7fPxxVvOZ/Zn93u+9Dy9nZs/eA9rb2wUAAFozaF0AAABCEEgAAEkQSAAAKRBIAAApEEgAACkQSAAAKRBIAAApEEgAACkQSAAAKRBIAAApEEgAACkQSAAAKRBIAAApEEgAACkQSAAAKRBIAAApEEgAACkQSAAAKRBIAAApEEgAACkQSAAAKRBIAAApEEgAACkQSAAAKfhTIDU0NMTGxhYVFWldCACg7/lTINlstgsXLrhcLq0LAQD0PZPWBfRKRUXFV199tW3bNq0LAQD0Fz8IpBkzZpw7d07rKgAA/csPAmnDhg1Op1MI8Ze//OXw4cNalwMA6Bd+EEj333+/8h9FRUUEEgDolR8E0pUJ+d2zgfWVll9PaF0IAEihvLxc6xJ6oNtAclr+I23ji/952/VaF6KliIgI+Q/B/kYT6ICgCUJERERoXULP/Om0bwCAjhFIAAApEEgAACkQSHp2jX9orqAJdEDQBD9BIAEApEAgAQCkQCABAKTgT4H0yiuvlJeX33PPPVoXAgDoe/4USAAAHSOQAABSIJAAAFIgkAAAUiCQAABS0HMgtbg8WpcAAOgtPQdSdb1T6xIAAL2l50ACAPgRAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBd0G0gCPS+sSAACXQbeBZGqt17oEAMBl0G0gAQD8C4EEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkIJJ/VWWlJTs3LnTbreHhIRMnz49MTFx+PDh3U+x2Ww5OTknTpwYOXJkTEzMsmXLbrrpph5X1Nrm6ZuKAQD9T+13SDabbdGiRbm5uRMmTGhtbc3IyIiPjz9//nw3U55++um1a9cePnx48uTJZrM5KysrLi7uyJEjPa6rut7Zd4UDAPqXqoFUVVW1cePGUaNGHThwYMeOHTk5OUlJSZWVlWlpaZeasnfv3tzc3Ntuu+2bb75JT0/fu3dvRkaGy+Vas2aNx8MbIADQD1UDyWazOZ3OpKSkMWPGKCOrV68ODw/ft2+fw+Hockp+fr4QYsWKFYMHD1ZGZsyYMWPGjOrq6tLSUnXKBgCoQNVAKi4uFkLMnTu34+DMmTPdbvfRo0e7mWg0Gjs+DAoKEkK4XK5+qBEAoA1VA8lut1sslpCQkI6DMTExQohTp051OWXGjBlCiL1793pHKioqCgoKgoODo6Ki+rNYAICq1DvLzuPxOJ3O0aNH+4wPHDhQCHH27NkuZ8XHx//0008ZGRllZWWRkZEul+vgwYMWi+Xdd981GHpI0z179hRtyPQ+LC8vv7otAAB/EhERoXUJl0e9QHK73UII77dHXkquKM92KTo6Ojg4uKysrKysTBmZPn16eHh4j2ucP3/+R/973ZVXDAD+rOO/wv0inNQ+7ftyT4377LPPkpOTw8LC9uzZU15efvLkyXXr1hUUFDz22GOXOg8CAOCP1Asks9lsNBo7nxqnnJsQGBjYeYrH43n11VeHDBmyY8eOyMhI5UUSEhKee+65M2fOZGZmdp4CAPBTqr5DCgkJqa+v9xmsq6sTQkyZMqXz8idOnKivr4+NjR06dGjH8Xnz5gkhvv/++36rFACgNlUDKSoqyul0+pxcUFhYKIQYN25c5+UtFosQoq2tzWe8qalJ/OvLJwCAPqj6N135BVJ2drZ3xOFw5OXlhYaGRkdHKyO1tbV2u/27774TQtx8882jRo3Kz8+vqqrq+DpZWVlCiIkTJ6pXOgCgn6kaSLNnzx4/fnxmZubu3btdLldFRcXKlSubm5tTUlK8y9jt9sWLF69YsUJ5+Nxzz7W1tT3++OM2m628vLyoqOiVV15JT08PDQ1NSEhQs3gAQL9S9WrfBoMhPT196dKlqampqampymBycnJ8fPylpixYsKCpqWnLli1r1671Dk6cOHHbtm3eiwkBAHRgQHt7u8qr9Hg8+fn5NTU1JpMpKiqq8y+TOmtoaCgqKnI4HGaz+cYbb/R+vteN4If/18OLl330+IS+KBkA/FtERIT8FwfQ4H5IBoNh1qxZlzVl8ODBs2fP7qd6AAAy4EQ1AIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFHQbSAM8Lq1LAABcBt0GkqnVoXUJAIDLoNtAAgD4FwIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAUCCQAgBQIJACAFAgkAIAU9B1JLm0frEgAAvaXnQKp2OLUuAQDQW3oOJACAHyGQAABSIJAAAFIgkAAAUiCQAABSIJAAAFIgkAAAUiCQAABSIJAAAFIgkAAAUiCQAABSIJAAAFIgkAAAUiCQAABSIJAAAFIgkAAAUiCQAABSIJAAAFIgkAAAUiCQAABSMKm/ypKSkp07d9rt9pCQkOnTpycmJg4fPrz7Kb/99tt777337bffXrhwYerUqQkJCVOnTlWnWgCAOtR+h2Sz2RYtWpSbmzthwoTW1taMjIz4+Pjz5893M+WHH354+OGHP/74Y7PZPHbs2K+++mrx4sV79+5VrWYAgApUDaSqqqqNGzeOGjXqwIEDO3bsyMnJSUpKqqysTEtLu9QUt9u9cuXKxsbGzMzMPXv27N69Oysry2Kx/OEPf2hqalKzeABAv1I1kGw2m9PpTEpKGjNmjDKyevXq8PDwffv2ORyOLqccOHDgzJkzK1asuP3225WRyMjIJ554orGx8euvv1apbgBA/1M1kIqLi4UQc+fO7Tg4c+ZMt9t99OjRLqccOnRICDF//vyOg6tWrSovL58zZ06/VQoAUJuqgWS32y0WS0hISMfBmJgYIcSpU6e6nHLs2DGLxRIWFubxeI4dO1ZQUHCp91IAAL+m3ll2Ho/H6XSOHj3aZ3zgwIFCiLNnz3Y569dffx0/fnxeXt6GDRvOnTunDD7wwAObNm0aOXJk92s8evRoxNbHvQ/Ly8uvvHoA8DcRERFal3B51Askt9sthPB+e+RlMBi8z/pwuVxut7u2tvapp56KjY19/vnnPR5PVlbWwYMHKysrs7OzlbmXMm3atEM7CSEA16iO/wr3i3BS+7Rvj8dzuVOqq6uffPLJjz76KC4uLj4+/pNPPrn33ntLS0uzsrL6o0IAgCbUCySz2Ww0GktLS33GXS6XECIwMLDzFKPRqPxvSkpKx/ElS5YIIY4cOdJftQIAVKfqlRpCQkJqamp8Buvq6oQQU6ZM6by8wWCwWCwGg8EnriZNmiSEaGxs7GZdA9yuqy0XAKAiVT+yi4qKcjqdPicXFBYWCiHGjRvX5ZTbbrutubnZ5xumH3/8UQhhsVi6WZexte5qywUAqEjVQFJ+gZSdne0dcTgceXl5oaGh0dHRykhtba3dbv/uu++Uh3fddVdbW9uHH37Y8XU+/fRT5SmV6gYA9D9VP7KbPXv2+PHjMzMzx4wZs3DhwsrKyvXr1zc3N2/YsMG7jN1uT05OHjFixOHDh4UQCxcu3Llz55tvvtnc3HzPPfe0trbu379/z549U6ZMiY+PV7N4AEC/UjWQDAZDenr60qVLU1NTU1NTlcHk5ORuoiUwMPDDDz9cvnz59u3bt2/frgxOmzZt69at3Z/zDQDwLwPa29tVXqXH48nPz6+pqTGZTFFRUZ1/mdTllMLCwnPnzhkMhrFjx9566609TgmPfejm/5F66OnovigZAPxbRESE/BcH0OB+SAaDYdasWZc7ZcaMGf1UDwBABnzqBQCQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkIKeA6m6vlXrEgAAvaXnQGpp82hdAgCgt/QcSAAAP0IgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkoNtAMnjcLS6P1lUAAHpLt4FkbKmrrndqXQUAoLd0G0gAAP9CIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApEAgAQCkQCABAKRAIAEApKDbQBrgcQkhWtq4JRIA+AfdBpKp1SGEqHZwSyQA8A+6DSQAgH8hkAAAUtAgkEpKSp599tnY2Ni4uLi0tLTa2trez21qaoqLi/v973/ff+UBADShdiDZbLZFixbl5uZOmDChtbU1IyMjPj7+/PnzvZy+YcOGsrKyxsbGfi0SAKA+VQOpqqpq48aNo0aNOnDgwI4dO3JycpKSkiorK9PS0nozff/+/Z9//nl/FwkA0ISqgWSz2ZxOZ1JS0pgxY5SR1atXh4eH79u3z+FwdD+3qqpq/fr1d999d/+XCQDQgKqBVFxcLISYO3dux8GZM2e63e6jR492M9Hj8aSkpAwZMmTTpk39WyIAQCOqBpLdbrdYLCEhIR0HY2JihBCnTp3qZuI777xz/PjxzZs3Dxs2rH9LBABoxKTamjwej9PpHD16tM/4wIEDhRBnz5691ES73b59+/akpKSpU6c2NTVd1krvvfdec/MF5b/Ly8svs2QA8GMRERFal3B51Askt9sthPB+e+RlMBi8z3bmcDheeOGFyMjIVatWXcFK8/LybhoZdAUTAcDfdfxXuF+Ek3qBpPB4Lu/ici+99FJdXV1mZqaSWwAAvVIvkMxms9FoLC0t9Rl3uVxCiMDAwM5T8vPzv/jii0cfffTMmTNnzpwRQrS2tgohampqioqKTCbTtGnTullj6JCAPqseANDPVH2HFBISUlNT4zNYV1cnhJgyZUrn5ZWLONhsNpvN1nG8uLh4yZIlFoulpKSkm9UFmXlTBQB+Q9VAioqKysnJKS8v7/hpZmFhoRBi3LhxnZePiYl54403Oo60tLSsX78+MjLyiSeeMBqN/V0wAEA1qgbS3Llzc3JysrOz161bp4w4HI68vLzQ0NDo6GhlpLa29qeffjKbzbfeemt4eHh4eHjHV2hqalq/fn1YWFhcXJyalQMA+puqH2rNnj17/PjxmZmZu3fvdrlcFRUVK1eubG5uTklJ8S5jt9sXL168YsWKPlljm6e9T14HANDfVA0kg8GQnp4+bty41NTUSZMmzZ49+/Dhw8nJyfHx8f20xrN1Lf30ygCAvqX2ad/XX3/9/v378/Pza2pqTCZTVFSUzy+T7r///m5+wTpo0CB+3woAuqR2IAkhDAbDrFmz1F8vAEBmnBgNAJCCngNp+EAN3v8BAK6MrgMpyKx1CQCA3tJzIAEA/AiBBACQAoEEAJACgQQAkAKBBACQAoEEAJCCngMp2GI+3+jSugoAQK/oOZAGBxobWt1aVwEA6BU9BxIAwI8QSAAAKRBIAAApEEgAACnoOZBMhgHcwhwA/IWeAyl0aEC1w6l1FQCAXtFzIBkH8A4JAPyGngMJAOBHCCQAgBQIJACAFAgkAIAUCCQAgBT0HEjDB5rqWtq0rgIA0Cs6D6TaZgIJAPyDngMJAOBHCCQAgBQIJACAFPQcSNzCHAD8iJ4DiVuYA4Af0XMgcfsJAPAjeg6kYEtAdX2r1lUAAHpFz4EUZDLwDgkA/IWeA8lkHNDmJpAAwD/oOZBuGBZ4to6P7ADAP+g5kAAAfkTngRRkMnDmNwD4BZ0HUujQAH4bCwB+QeeBBADwFzoPpCCTocHJR3YA4Ad0HkihQwLPNzq1rgIA0DOdB9LgQGOLy6N1FQCAnuk8kIIt5up63iEBgB/QeSCZDAO0LgEA0Cs6D6Rgi7m2uU3rKgAAPdN5IFkCjBea+B0SAPgBnQfSTSODTte0aF0FAKBnOg+k0CEBnNQAAH5B74E0NLDawQW/AcAP6D2QhgRwBwoA8As6D6Rgi7mh1d3Sxm9jAUB2Og8kIcQNwwOrHXyNBACy038gcV4DAPgF/QfSDcOCTtc0a10FAKAH+g+km4MH/nSeQAIA2ek/kCJDLUf+j0PrKgAAPdB/IE0NH3r8l3qtqwAA9ED/gXTLdYManG7OawAAyek/kIQQNwcPPFnVoHUVAIDuXBOBNOeWUbk/XtS6CgBAd66JQHp08nV//5lAAgCpXROBFDV68OmaFvsZTm0AAHldE4EkhFhx5+g/5lVoXQUA4JKulUB6ZvoNfz35W07ZBa0LAQB07VoJpGCL+bMlk//nJ98f/4XT7QBARtdKIAkhHpoY/Prvxt25zb56749cTAgAZGPSugBVLbs97KGJwX/48p93brO3tHluDh54y3WWIJNBCBE6JCDQ9G/x3NLmCTJdeWAHmgytV3cfpmFBprqWtqt5hat09QUYBgzwtLf3VT1X4Oo3YXCgsaHV3Vf1XAGTUbRd3fqvchOu/kgePshY23RV22A0DHB7rvxA8rQLw4CrWX8fbMJVtvHqj8OLY2ddzXR1XFuBJIQIHRLw3qMR7z0aUV3vPP5LvXIFh9M1Lb1/hbqWttrmNiHE4EDjqEFmn2d/qWtt87SHDgkYHGhs/ddfQnd7+9naViHE8IGmYUH/v+eNTvf5RpcQIshs+I/BAd5xd3t7tcPZ0uZubWvvuHzH1+8Yn8ryQWaDTz2/NjhbXB6flXpfxKd+b5Emw4DQoQHGAQN8lvcpUghxocnV0Oq+YXhgx4WVjTIZBoweFthl33yKV4oUQgRbzJYAYy+Xv1SRnTtTcbFFCNFlkZ1f3LvSS+2pLpfvfCR4Ox9oNASZDT4v0nl3d9l5ZflBAQZPu+h8eHTe2Et1/rcG17kG5+Xu7qvvvFLkDcP/rRjvi/gU2f1KDQPEyEFdHB7ddN5n99U0dbH7lANYXObu7jze2uaprnde6pi/VOd7/9fgbG3ryEFmT/tV/TXwC9dcIHmFDgmYc8sorasAADX897OHtC6hZxoEUklJyc6dO+12e0hIyPTp0xMTE4cPH979FJvNdujQoWPHjplMpqioqAULFsyYMUOdagEA6lD7pAabzbZo0aLc3NwJEya0trZmZGTEx8efP3++mynPPPPM2rVrDx06NHHixPDw8IMHDy5btuzNN99UrWb/FRERoXUJ2qMJdEDQBD+haiBVVVVt3Lhx1KhRBw4c2LFjR05OTlJSUmVlZVpa2qWm7N+//+DBgxMnTiwoKEhPT9+9e/f+/fuvu+669PT0o0ePqlk8AKBfqRpINpvN6XQmJSWNGTNGGVm9enV4ePi+ffscjq7voffll18KIV5++eWRI0cqI2PHjl2zZo0Q4rPPPlOlagCAGlQNpOLiYiHE3LlzOw7OnDnT7XZf6u3ON998ExAQEBMT03Fw8uTJQoj6eq5NBwD6oWog2e12i8USEhLScVAJm1OnTnU5ZcuWLe+//77P4Pfffy+E6PFUCACAH1HvLDuPx+N0OkePHu0zPnDgQCHE2bNnu5zV+Wy6mpqaLVu2CCEWLlzY/Rr5GlPQBCEETaADQgia4A/UCyS32y2E8H575GUwGLzP9ignJ+f111+vrq5OSUlRPri7lPLy8iutFACgAbV/h+TxXOHFM0pLS1NTU48dOzZixIjNmzc/9NBDfVsYAEBb6gWS2Ww2Go2lpaU+4y6XSwgRGBjY1aT/Jy0tLSMjIyAgIDExcfny5UOHDu3HQgEAWlD1HVJISEhNTY3PYF1dnRBiypQpl5r19NNP5+bmPvjggy+++OL111/fvyUCADSi6ll2UVFRTqfT59udwsJCIcS4ceO6nLJr167c3NyEhIRt27aRRgCgY6oGkvILpOzsbO+Iw+HIy8sLDQ2Njo5WRmpra+12+3fffac8/Pjjj4OCgpRfwgIAdEzVj+xmz549fvz4zMzMMWPGLFy4sLKycv369c3NzRs2bPAuY7fbk5OTR4wYcfjw4dra2oqKihEjRixfvrzzq8XGxiYmJqpYPgCgH6kaSAaDIT09fenSpampqampqcpgcnJyfHx8l8ufOHFCCHHx4sWioqLOz1533XX9VyoAQGUD2lW/oafH48nPz6+pqVHuJdH5l0kAgGuQBoEEAEBnat8PCQCALhFIAAApEEgAACmofS07FZSUlOzcudNut4eEhEyfPj0xMVFnN6qw2WyHDh06duyYclbIggULfK6J3psO6KNLTU1NixYtuuWWW/74xz/6PKX7Jvz222/vvffet99+e+HChalTpyYkJEydOtVnmR430K87IISw2Ww5OTknTpwYOXJkTEzMsmXLbrrppo4L6PIwaGhomD179ptvvhkbG9v52T7Z6Vr1RG8nNdhstrVr15pMpjvuuOOXX3755z//GRYWlpWVFRwcrHVpfeOZZ545ePCgyWSKjY1tbm4uKSlxu92JiYnPP/+8skBvOqCbLr3wwguff/75Aw888M4773Qc130TfvjhhyVLlly8eDEyMnLQoEHKYfDGG2/ExcV5l+lxA/26A+JfFxULCAi44447zp07V1ZWFhQU9MEHH9x+++3KAno9DHbt2vXaa6998MEH99xzj89TfbLTtexJu45UVlZOmjTpzjvvPH36tDKyZcsWq9W6Zs0abQvrK/v27bNarY888siFCxeUkZ9//nn69OlWq/XIkSPtveuAbrqkdMNqtSYnJ3cc130T2tra7rvvvkmTJh0+fFgZOXnyZHR0dHR0dGNjozLS4wb6dQfa29v/+te/Wq3WxYsX19fXKyNff/31hAkT7r77brfb3a7Tw+D06dMZGRlTpkyxWq1///vffZ7tk52ubU90FUhvv/221WrdtWtXx8H77rtvwoQJdXV1WlXVh1auXGm1WouLizsOfv7551ar9b/+67/ae9cBfXSpsrIyOjp62bJlnQNJ901Q9vjbb7/dcfBPf/qT1Wr929/+pjzscQP9ugPt7e2rV6+2Wq1ff/11x8Hly5dbrdZ//OMf7Xo8DJR/enp1DqQ+2ena9kRXJzUUFxeLf10xz2vmzJlut/vo0aMaFdWXvvnmm4CAAOWm717KjQrr6+tF7zqggy55PJ6UlJQhQ4Zs2rSp87O6b8KhQ4eEEPPnz+84uGrVqvLy8jlz5igPe9xAv+6Al9Fo7PgwKChI/OuONvo7DDZs2PDWW2+99dZbd9xxR5cL9MlO17Ynugoku91usVhCQkI6Dip/vk+dOqVRUX1py5Yt77//vs/g999/L4RQvnLsTQd00KV33nnn+PHjmzdvHjZsWOdndd+EY8eOWSyWsLAwj8dz7NixgoICh8Phs0yPG+jXHRBCKCfy7N271ztSUVFRUFAQHBwcFRUl9HgY3H///fPmzZs3b154eHiXC/TJTte2J/o5y87j8TidztGjR/uMDxw4UAhx9uxZLYrqYz5n0wkhampqtmzZIoRYuHBhbzqggy7Z7fbt27cnJSVNnTq1qanJ59lroQm//vrr+PHj8/LyNmzYcO7cOWXwgQce2LRp08iRI0UvNtDfOyCEiI+P/+mnnzIyMsrKyiIjI10u18GDBy0Wy7vvvmswGK6Fw8BHn+x0zXuin0Byu91CiM5XxjMYDN5ndSYnJ+f111+vrq5OSUmZPHmy8klF9x3w9y45HI4XXnghMjJy1apVXS7Qmw306ya4XC63211bW/vUU0/FxsY+//zzHo8nKyvr4MGDlZWV2dnZBoOhxw306w54RUdHBwcHl5WVlZWVKSPTp09X3j3o/jDorE92uuY90dVHdkIIj8ejdQlqKC0tffzxx1etWtXa2rp58+annnrK+1RvOuC/XXrppZfq6uq2bt2q/D/kUvTdBCFEdXX1k08++dFHH8XFxcXHx3/yySf33ntvaWlpVlaWd5keN9CvO/DZZ58lJyeHhYXt2bOnvLz85MmT69atKygoeOyxx7wfYOr+MOisT3a6hj3RTyCZzWaj0VhaWuozrrxvCAwM1KKofpGWljZ//vzS0tLExMQvv/zyoYceUsZ70wG/7lJ+fv4XX3wxZ86cM2fOFBUVFRUVHTlyRAhRU1NTVFSkfOOq+yYoX+MbjcaUlJSO40uWLBFCKA3pcQP9ugNCCI/H8+qrrw4ZMmTHjh2RkZFCCLPZnJCQ8Nxzz505cyYzM1P3h0FnfbLTNe+Jfj6yE0KEhITU1NT4DNbV1QkhpkyZokVFfU/5MeCDDz744osvdr6ne2864L9dqq2tFULYbDabzdZxvLi4eMmSJRaLpaSkROi9CQaDwWKxGAwGn78OkyZNEkI0NjYqD3vcQP/tgBDixIkT9fX1Dz744NChQzuOz5s377XXXlNO89H30k0BmAAAAs1JREFUYdClPtnp2vZEP++QhBBRUVFOp7O8vLzjYGFhoRBi3LhxGhXVl3bt2pWbm5uQkLBt27bOaSR61wH/7VJMTMwb/0457TsyMtL730LvTRBC3Hbbbc3NzT4f6P/4449CCIvFojzscQP9ugPKZra1tfmMKye5KB/n6v4w6KxPdrq2PdFVICnnzmdnZ3tHHA5HXl5eaGhodHS0dnX1mY8//jgoKGjNmjWXWqA3HfDfLoWHh8f9u4cfflgIERYWFhcX5/3oUt9NEELcddddbW1tH374YcfBTz/9VHlKedjjBvp1B26++eZRo0bl5+dXVVV1HFe+Qps4caK4Bg6Dzvpkp2vck/7+5a2a3G737373O6vV+uc//9npdJ4+fTohIcFqte7Zs0fr0vrAxYsXrVbr7bff/p9d+eCDD9p71wE9damxsbHzlRp034SWlpZZs2ZZrdatW7ceP378yJEjL7/8stVqXbBggXLVnPZebKBfd6C9vT0rK8tqtd59993Z2dllZWWFhYWpqanKiHIxIR0fBuvWrevySg19stO17YneLq5aVVW1dOnSn3/+2TuSnJz87LPPalhSXykoKFi2bNmlnn3kkUeUK173pgO66VJTU1N0dHTni6vqvgmnT59evnx5RUWFd2TatGlbt25Vfoek6HED/boDQojMzMwtW7Y0Nzd7RyZOnLht2zbv70b1ehi89NJLWVlZXV5ctU92uoY90VsgCSE8Hk9+fn5NTY1yd4bO59TrXm86oPsu6b4JHo+nsLDw3LlzBoNh7Nixt956a5fLdL+Bft0BIURDQ0NRUZHD4TCbzTfeeGPnz5R0fxh01ic7Xaue6DCQAAD+SFcnNQAA/BeBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJACgQQAkAKBBACQAoEEAJDC/wXHzBIS1LtARwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"iv = [1, 1, 1];\n",
"[t, xa] = ode45(f, t0, iv);\n",
"plot(t, xa(:,3))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"But then watch this."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAACXBIWXMAABcSAAAXEgFnn9JSAAAAB3RJTUUH4wkIAQgn+LHcRgAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAwOC1TZXAtMjAxOSAwMTowODozObWlhRsAACAASURBVHic7N1reFvVmTf827JkSZZsy7bsKLHdOMRxHGKCEwKkOZTQpinMwJWmZQZ46QPDA5SUcMpMaaftFJjSzvN2Ap0eYE4JLfCkM/Ql7pQmw5iSAIEkOCAHGRSQg0NsbMuyosiyDtbBsvx+WIlQZFnaR3lb/H8fuPC2lbUlb697He61VtHU1BQBAADMNtVs3wAAAAARAhIAACgEAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACgCAhIAACjCXApIwWBw7dq1R48ene0bAQAA6c2lgNTe3n727NmJiYnZvhEAAJCeerZvgJP+/v4DBw788pe/nO0bAQAAucyBgLRhwwa32z3bdwEAAPKaAwHpkUceicViRPS73/2us7Nztm8HAABkMQcC0qZNm9j/HD16FAEJAKBQzYGAJEzNn9+vDTgNI+/P9o0AAChCT0/PbN9CDgUbkGKGeTsf/d5fXT5/tm9kNi1dulT5j6Dc8CHgEyB8CERLly6d7VvIbS6lfQMAQAFDQAIAAEVAQAIAAEVAQCpkn/FBcwYfAj4BwocwRyAgAQCAIiAgAQCAIiAgAQCAIsylgPTjH/+4p6fnqquumu0bAQAA6c2lgAQAAAUMAQkAABQBAQkAABQBAQkAABQBAQkAABShkANSZCIx27cAAABcFXJAcgVis30LAADAVSEHJAAAmEMQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEQkAAAQBEKNiAVJSZm+xYAAICHgg1I6mhgtm8BAAB4KNiABAAAcwsCEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKAICEgAAKII6/0UeP378mWeesVqtNTU169evv+uuu0wmU/aXtLe3d3R0vP/++1VVVatWrbrzzjsbGxvzcrMAAJAn+e4htbe333zzzQcPHly2bFk0Gt29e/fWrVs9Hk+Wl9xzzz3f//73Ozs7L7nkEo1G88ILL2zZsuXYsWN5u2cAAMiDvAak4eHhRx99tLq6+qWXXnr66ac7Ojq2bdvmdDp37tw500tefPHFgwcPXn755W+99dauXbtefPHF3bt3T0xMfOc730kkEvm8eQAAkFVeA1J7e3ssFtu2bdvChQvZlR07djQ0NOzbt8/v92d8yaFDh4jo7rvvNhqN7MqGDRs2bNjgcrlOnDiRn9sGAIA8yGtA6urqIqJrr7029eLGjRsnJyfffvvtLC8sLi5O/VKn0xHRxAQOKQcAKBx5DUhWq9VgMNTU1KReXLVqFRGdPHky40s2bNhARC+++GLySn9//5tvvmk2m9va2uS8WQAAyKv8ZdklEolYLFZXV5d2Xa/XE9Hg4GDGV23durW3t3f37t0Oh2P58uUTExOvvPKKwWD453/+Z5UqRzR98skn//P+/ckve3p6xL0DAIC5ZOnSpbN9C/zkLyBNTk4SUXL2KInFFfbdjFauXGk2mx0Oh8PhYFfWr1/f0NCQs8R777330a88IfyOAQDmstRW+JwITvlO++abGvdf//Vf27dvX7Bgwe9///uenh673f6DH/zgzTff/Mu//MuZ8iAAAGAuyl9A0mg0xcXF01PjWG6CVqud/pJEIvGTn/ykrKzs6aefXr58OftHbr311r/+678eGBh47rnn8nDbAACQH3ntIdXU1AQCgbSLY2NjRLRixYrpP//+++8HAoG1a9eWl5enXv+zP/szIvrggw9ku1MAAMi3vAaktra2WCyWllxw5MgRIlq8ePH0nzcYDEQUj8fTro+Pj9P5yScAACgMea3T2QqkvXv3Jq/4/f5XX33VYrGsXLmSXfH5fFartbu7m4iampqqq6sPHTo0PDyc+u+88MILRHTxxRfn79YBAEBmeQ1ImzdvXrJkyXPPPbdnz56JiYn+/v777rsvHA4/+OCDyZ+xWq233HLL3Xffzb7867/+63g8ftNNN7W3t/f09Bw9evTHP/7xrl27LBbLrbfems+bBwAAWeV1t2+VSrVr16477rjjsccee+yxx9jF7du3b926daaX3HDDDePj4z/72c++//3vJy9efPHFv/zlL5ObCQEAQAEompqaynORiUTi0KFDXq9XrVa3tbVNX5k0XTAYPHr0qN/v12g0n/vc55Lje1mYr/+be++999GvLJLilgEA5ralS5cqf3OAWTgPSaVSXX311bxeYjQaN2/eLNP9AACAEhRsolpRYiIax/kUAABzRsEGJHXU7wrEZvsuAACAq4INSAAAMLcgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIgIAEAgCIUckCKxBOzfQsAAMBVIQcklz8227cAAABcFXJAAgCAOQQBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFAEBCQAAFEGd/yKPHz/+zDPPWK3Wmpqa9evX33XXXSaTKftLzpw58y//8i/vvPPO2bNnV69efeutt65evTo/dwsAAPmR7x5Se3v7zTfffPDgwWXLlkWj0d27d2/dutXj8WR5yYcffnj99df/9re/1Wg0F1100YEDB2655ZYXX3wxb/cMAAB5kNeANDw8/Oijj1ZXV7/00ktPP/10R0fHtm3bnE7nzp07Z3rJ5OTkfffdFwqFnnvuud///vd79ux54YUXDAbD3//934+Pj+fz5gEAQFZ5DUjt7e2xWGzbtm0LFy5kV3bs2NHQ0LBv3z6/35/xJS+99NLAwMDdd9995ZVXsivLly+/7bbbQqHQG2+8kaf7BgAA+eU1IHV1dRHRtddem3px48aNk5OTb7/9dsaXvPbaa0T0ta99LfXiAw880NPTc80118h2pwAAkG95DUhWq9VgMNTU1KReXLVqFRGdPHky40veffddg8GwYMGCRCLx7rvvvvnmmzP1pQAAYE7LX5ZdIpGIxWJ1dXVp1/V6PRENDg5mfNXIyMiSJUteffXVRx55xO12s4tf/vKXf/SjH1VVVWUv8e233176i5uSX/b09Ai/ewCAuWbp0qWzfQv85C8gTU5OElFy9ihJpVIlv5tmYmJicnLS5/N961vfWrt27be//e1EIvHCCy+88sorTqdz79697LUzueKKK157BkEIAD6jUlvhcyI45TvtO5FI8H2Jy+X63//7f//mN7/ZsmXL1q1b/+M//uOLX/ziiRMnXnjhBTnuEAAAZkX+ApJGoykuLj5x4kTa9YmJCSLSarXTX1JcXMz+++CDD6Zev/3224no2LFjct0rAADkXV53aqipqfF6vWkXx8bGiGjFihXTf16lUhkMBpVKlRauWltbiSgUCsl2pwAAkG95HbJra2uLxWJpyQVHjhwhosWLF2d8yeWXXx4Oh9NmmD766CMiMhgMst0pAADkW14DEluBtHfv3uQVv9//6quvWiyWlStXsis+n89qtXZ3d7Mv161bF4/Hf/3rX6f+O88//zz7VpayiiYnpL15AACQVV6H7DZv3rxkyZLnnntu4cKFN954o9PpfPjhh8Ph8COPPJL8GavVun379srKys7OTiK68cYbn3nmmccffzwcDl911VXRaHT//v2///3vV6xYsXXr1ixlFUfHZH8/AAAgnbwGJJVKtWvXrjvuuOOxxx577LHH2MXt27dnCS1arfbXv/71N7/5zaeeeuqpp55iF6+44opf/OIX2XO+AQBgbimamprKc5GJROLQoUNer1etVre1tU1fmZTxJUeOHHG73SqV6qKLLrr00ktzvqRh7XVN33jstXtWSnHLAABz29KlS5W/OcAsnIekUqmuvvpqvi/ZsGGDTPcDAABKgFEvAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQBAQkAABQhEIOSMFYfLZvAQAAuCrkgOQJTcz2LQAAAFeFHJAAAGAOQUACAABFQEACAABFQEACAABFQEACAABFQEACAABFQEACAABFQEACAABFQEACAABFQEACAABFKNiAVDwRme1bAAAAHgo2IKkmxmf7FgAAgIeCDUgAADC3ICABAIAiFHJA8oVxHhIAwJyBgAQAAIpQyAEJAADmEAQkAABQhIINSKrJ6GzfAgAA8FCwAak4FpztWwAAAB4KNiABAMDcgoAEAACKgIAEAACKgIAEAACKgIAEAACKgIAEAACKgIAEAACKgIAEAACKUMgBSa0qiiemZvsuAACAk0IOSJayElcgNtt3AQAAnMxCQDp+/Pj999+/du3aLVu27Ny50+fzcX/t+Pj4li1bvvvd73L5YXVxUXwSPSQAgLkh3wGpvb395ptvPnjw4LJly6LR6O7du7du3erxeDi+/JFHHnE4HKFQSNabBACA/MtrQBoeHn700Uerq6tfeumlp59+uqOjY9u2bU6nc+fOnVxevn///j/+8Y9y3yQAAMyKvAak9vb2WCy2bdu2hQsXsis7duxoaGjYt2+f3+/P/trh4eGHH374C1/4gvy3CQAAsyCvAamrq4uIrr322tSLGzdunJycfPvtt7O8MJFIPPjgg2VlZT/60Y/kvUUAAJgleQ1IVqvVYDDU1NSkXly1ahURnTx5MssLn3zySZvN9sQTT1RUVMh7iwAAMEvyF5ASiUQsFqutrU27rtfriWhwcHCmF1qt1qeeemrbtm2rV6+W9xYBAGD2qPNW0uTkJBElZ4+SVCpV8rvT+f3+hx56aPny5Q888ADfEoeGhr74xW9rwmfZlz09PXz/BQCAuWvp0qWzfQv85C8gMYlEgtfP/93f/d3Y2Nhzzz3H4hYvdXV1r776amOVju8LAQAKQGorfE4Ep/wFJI1GU1xcfOLEibTrExMTRKTVaqe/5NChQy+//PLXv/71gYGBgYEBIopGo0Tk9XqPHj2qVquvuOIK+W8cAADyIa89pJqaGq/Xm3ZxbGyMiFasWDH959kmDu3t7e3t7anXu7q6br/9doPBcPz4cdluFgAA8iqvAamtra2jo6Onpye183jkyBEiWrx48fSfX7Vq1T/+4z+mXolEIg8//PDy5ctvu+224uLi7MVhc1UAgDkkrwHp2muv7ejo2Lt37w9+8AN2xe/3v/rqqxaLZeXKleyKz+fr7e3VaDSXXnppQ0NDQ0ND6r8wPj7+8MMPL1iwYMuWLTmLs5RpXYFok1kv+RsBAADJ5XUd0ubNm5csWfLcc8/t2bNnYmKiv7//vvvuC4fDDz74YPJnrFbrLbfccvfdd4svDj0kAIA5JK8BSaVS7dq1a/HixY899lhra+vmzZs7Ozu3b9++devWfN4GAAAoUL7TvufPn79///5Dhw55vV61Wt3W1pa2MmnTpk1ZFgyVlpZiOREAQEHKd0AiIpVKdfXVV+e/XAAAULJCPjEWAADmEAQkAABQBAQkAABQBAQkAABQhEIOSEZtcTCaeRNxAABQmkIOSGaDxhOamO27AAAATgo5IIEw8cRUZ7/fF47noayfvzHwr0eH8lDQgZPea3d193rCchcUiSd2vPhRft7U8++O/MWz9jy0ulyB2M17TuzpcsldEBH969Ghm/ecyMPYhsM9fv3T7x04mb7dsxx+2PHxvb8/mYeNYzr7/dc//Z51ICB3QTJBQIJ03/vvU5//pfX6p9+T++/nmXeGd7z40bfae/5gPyNrQYNj0euffq/XM/6//iP99BPJ7Xjxo+ffHXlof29nv1/Wgno94duf/7Czf+xbe2VfKn778x9aB/y3P/9hnzcia0EdjrPfau95/t2RH3Z8LGtB8cTU1t+853CH/uI5u9xtr2feGX78tU/+v+6R3Z1OWQsKRie3/uY921Dgrhc+nKO7piEgwQX6vJEnDw9+9L3P93nDsjYe44mpH/7Px/vuWPGrrc1PHZG3P7G32/3VS2o+/O6aPm9E1jgRjE7+4f0zb92/+m+u+txTRwblK4iIfvpq/73r6z/87pr9H3hk7STZXSGHO/T+Q1fetHLe87YR+Qoiol+8Ofibm5Z9+N01e7pcslaprA304XfXtFqM+z/wyFdQPDH1kwN9L9zW+usblz11ZFDWN/XzNwZa5xtP/91aXzje2T8mX0HyQUCCCzzzzvBXL6lpMuu/+8WFz1plHKI5cNJr1BZf01J955oFnf1jrkBMvrKetQ4/sKFBrSr6q8vn/07OKvWZd4bXLKxorNL91eXz93a7I3F+5yNzF4xO7n3PvX1dvVFbvP4i0/Pvyvimfn3M+bUVtTq16pZV8559Z1i+gvq8kc7+sa+21rTUlpr0alnHnX73rnv7unq1qujrK2pe7pGx4WUdCKhVRdddbL6mpbpvNGIfDslXVvv77rs/v0CtKtq+rv7Zd/Ixviq5Qg5ISGoQoP099x1XzCei9YtMsjayXu7xfrW1Rq0q0qlVq+vL5SvLOhDwBCdWN5QR0bpFFbJWcy/3eG9sqyWixipdk1lvGwrKVJB10F9foW2s0hHRluXmrkEZ39Rrp0a/vKSSiDY1Vw2OReX7m3r91OjGxZUmvZqIvtpa0+E4K1NBkXjiwEfe6y42s4L2f+CRr+Py2+Oua1qqiUitKtq4uNI6KFcH3eEe7/NG2JvauLjy8GmfTAXJqpADkqEEad/8DI5FB8eiqxvKiaitzjjsj50YkatBt7fbvaW1hv3/6oayI6flCkjt77lvWjlPrSoiotUN5TZnQKbaJxJPvH5qdGNTJfuydb7R7pIrIL1o97Cqh4ha5xusA3JVc72e8KAvuqm5iojUqqKW2lL5EkPa3zvz9RXnHonL6su6nXJ9eq/3jjaZ9SycN1bpdGqVfHNje7vdN7bNY///pSWV8j3nf7Cfue5is06tIqKW2tI+b2Qu1n6FHJAKTzA66XCPy/fv7z/h2bSkijVRiWjFfMO78jS9He7xeGKK9VqIaHG1Xr4GnXUwsG5RBft/S1lJfHJKpo7L672jLbWllrIS9uW6xgr5ap/9H3g+Def15b2esEy1T4fj7DUt1SycE9GahRUy/abiiakDJ72sM0FEbXVl8kXZQ6d8yYKIqKXWYHPK8pwPjkWDscnkc95YpZNvJKDbGbzyc+Xs/43a4nqTVr7emHwQkOaMYHTyksePLftp59733DIVcehj31daqpJfrltk6jkjS/w7fNq3uqE8Wc2tv8hkd8nSFYvEE539Y2sWViSvrG6Qa3jw5R7vluU1qQXJVKV6QhOuQCxZzRm1xa3zDTLVPl2DgWQ1R0SX1ZfJNDxoHw5ZykvMBg37sqW2NBiblGlm8fVTo6lvas3CcseILM/5672j6xeZks/56oZyh3tcpqbDgZNe1pFl2haUyZ3nKQcEpDnj528MNFbq992x4of/87FMg06HP/atrr+g9pHpD/XQKd+XllQmv2ypLY1PTsnxh3r4Y59Jp072WohozcLykaAs1ZzdFWqZV5r8stVikKk7+/qp0TULK5LVHBG1LSiTabb8wElvchCS2JjnkCwB6cBH3k1LqlKvyBTRfeG43RVKrbsvayg79oksdffBj0ZTn/P6Cq2lrESOKGsbCkbiiSazPnmlZV5pj5yjKTIp5ICkVhVNTs3JZPyMfnvc9cAX6q+72BxPTMnRxu/1hOOJqdb5huSVllqDTLMgHY6zrZZPC1KrilrnG+ToJL3cc0F9SkSXLjDKEWUj8cThj30bF39allFbrNOo5EgBOHTKd9VFptQrCyt1ckRZ9ki01H4aZZvMepnmkA5+NHrV4gvelEzzVQc+8q5ZWMHmWphWi9Hhliucp44N0rlOkvRlsXyQ1DbKlZ8rlzVzVSaFHJDqTdpBX3S270Iar58ajcQTbB77mpbqwzJMTnQ4zm5suuCZtpSX9HrCkucuD45Fg9HJ9RdWqU3mUjmCX68nnNpEJaKWWkPvWekDUmf/WMu80uSIE9NqMchS+/SObmy64NNrrNLJUXenTSARkU6tkinKvt47mlZ3zzOWnB2XviD7cGh1fVnqFZmibDA66QrEUsM5K0uOfvMJV2haODfYh+XKCpFPIQekQnL447HrLjazquGyelly0rqdwbR2d32FVq0qkrw/YRsKpDVRichs0MgxwnDaG/58Y0XqFRZlJR/z7Oz3p3aPmCZzqeS1Tzwx5XCPs0zIpJZaQ69H+k+vfzSy9ML6lOSJsn3eiE6jSgvnMkXZrsHA8pTeORGpVUVGbbHkUbazf6ytzpgazomoulQzFpF+YwjWQ0q90lilYy0/ycuSFQLS3PByjzfZ0meLDCSvUh3u8dQpEKbJrJe8P3Hk9NiaheVpFxdX6yUfsgtGJ0+4QktrLnhTZoMmPjkl+WhGj3t8YaUu7aJJr5Y8ytqHQyxTOfWi2aBhiYvSltXZ75/+m5Ijyh4+7Vu/yJR2sclcKkeU7ewfm16WHFH28OmxVosx7aIcUdYViPV5IyyLPYlF2Txs3iitQg5IalWRfEvl88kViNmcgeSUb2OVzlhSLG2lEIknrAP+1Gw0pnW+UfKdvg6fHpted7fUlvZ5Jf7jYU3UtIs6taqxSid5TZexSl1aUyr5ApeMBZkNmshEQto2fiSesDkDqUkuzNIa6WfLj33iT6bmJ7GRNGmjrN0VMmqL0+puOhf8JH78Dp3yXbog/fGTI8p29o81VunS+pdE1LagrG8UAUkx6k1al3/uTetN93rv6MbFlUZtcfJKyzyDyy/l9Fhn/5hJr05rdxNRY6VO8s8wY+RrMpfKEfmmF8TKknZy0RWIeUIT04Nfk1kveY1wpG8sbbaAiIza4iazXtrg19k/1mTWpz51DBsLkrAgOtdrSf9NmQ0atapI2ijLhounX19crT91VuLfVK9nPFO3T/r5qj5vZHpBrCy5N8OVXCEHJCKao1vepjkxEkqbGm2s1EnbQ+rs96dNSzB1Fdr+USmfaVcgptOoUnP5mMYqnSsQk3bIO2PdTURNZr20Q3bWAX/bgrK02QIiaqk1SB7ObUOBjL+pxiq9tF3MXk94+ogTETWZpe/L9nkjTeb04WIispSVSPsB9nrCjdN650RUb9JKW3d7QhPB2OT0Ngrrx0gbZTP2L4locbV+ZK4l2hVyQLKUaV2BQsiye73Xl1axLq0tlXZTn/7RSMZnusms75M0ILHu0fS6W60qspSVSFtWnzfctqBs+vXqUo20UbbPG0ldApLEOhMSjhtH4ok+b6S+Qjv9W5JH2ROuUNrkPyP5Spo+b0SnVqWuFUtqrJK4i3nsE/9lDRkeifoKiTNyO/vHWmpLpz/nJMsHOGOUnXOZ34UckHRqVQH0kILRSZszkNYlb1tglHYhpHUgw/Q1EZkNJZ6QxLXP9BF8prFKJ2HT2xeOuwKxelOGuruxSidtc3h64laSSa+WsCyW0ZDc2ynVPGOJtFG2sz9D7gkRWcpLBn1RCf+yDp/2NVZlCOckQw/JOuBvrMzYdNAPjkn56R05PdZSm/mRkLYvG09M2YaCbXUZoqzZoEFAUhCjtjg/x57Kyu4KNZn1aXVQ63wpl/IFo5O2oWDGMRM56u7L6jP88RCRSa+RsCzrgL/VYpg+K0bnxmekbHdnTNxiGqt0gz7J3lTaNkippB10isQTdldoekYDEenUKrWqSMKarmswMH1oi1lYqZMwyrJte6YPF5MMJwPYnMGZnnNpe0isjTI9o4FkGMXNg0IOSGaDpiACUnB6U8tSVhKJJ6R6dzZnwKRXZxwzMRs08cSUhB/j4dO+jMNoRNRYJeV2AzPNihGR2SDl2CDris1UpUpb1kyzBcQGnaRr49uHQ42VuukZDUxjlU7CAS67K5S6s1x6QdIlUPR5w211xoxtFMmf8z5veKamg6WsRMJP7/Bp30wFzcUKsJADEhGpVUVzfdTuyOmxjH+ulrISqdp0fd5I6lanaeortFJVCp7QRMZsNGahpBl93c7gusbMf6iNlbpgdFKqB8PuCs40W0BEhhKVXbpjFPq8kYyzBURk1Eo5Njg4FplpGI2ILGVSdjHtw8GZqlRpO+jWgcz5IIyEz7kvHB8ci2Z5zoeki7LdWbtintDE3Fr6UuABqQB2D5ppOEjCSiF14/rpJOz4O9yhLHV3Y6WkzeHR8ExVqlFbLOHUTq8nXF+ROUgQ0SXzjZpMTXJh7K5g6/zM1VxjlU7CQadj/f6Z6lMiMunVUv2mWNLHTNOK9RVSTit2O4PTFwYlSficD45FLWUlM7bwTJJFPiLqG41knChlpO3L5kGBB6S5LstwkIRbT9qcwdWZUo8YCbti2evuepPWF5amoHhiyj4cap0h0YAkHeDqGgxkTNxi6solq30c5KRyGwAAIABJREFU7nFjSXHGkVUiYjWgVJMTnf3+6YuXkxqrdFLtf2Md8LMdqjJ+V9o8scOnfTP1L0nSuts2FJgpo4GI6iuknFacaRkAM+emkQo8IBlL5nZeQ5am1mKzXqr53j5vtjghYa5O9rq7vkInVUEO93hj1YxTIETUWKWXas3TTAmKjISrszv7x7JUPSRplWod9GdchMTUVUiWQHGs358xQ4xRq4qk2m/FFYi5ArEsZc0zSpZrkGvIQSfVtKLDPa5WFWVcBsBIvpRCbgUekCzlWmmzlvPMOjDjzHxjpTSbYrFOWMbUI6auQivVkPdMmcSMhFmRhz+ecaaXkSqriiXdZsxGYxqr9FI9gV2DgStn/vRYWZJ0MQfHopGJRJYhu/oKyToug2PRLMNoJN2ou8sfa6lN3449lUmvlmpz8cOnfVmec5NeLdUUZpb6gamvmGO71RR6QJLnOKy8OXJ6bKbhdUu5NG+NpUfPNGZC0uXqROKJ7HU325BNko5L9/CMM73JsiQJSK5AzGzQZOmKSTjd4vLHsow4kXTN4T5vuHV+5ox5pt6kDUalaTo43KGZMubPlVWhk2Rw1TqYq+6Wbr55poVBSVLVS33eSJZ2AxEZSoqjk0hqUIy5nmVnHfQvn5e57yLVPhQ2ZzD7H6pUfzwuf8xSVpKl7iYiS7k081WDvqilPPNcCyNVt88+nHn9VpKlrMQXjksy6GRzBmbKaGB0apUkzWHb0Ixpb0x9hTRjg6yNkr1KNWqLJWmj9I/OuCKbkarh5QvH1cVFWbpiJF3w6zmTYY/5VJKvI5RbgQckk14tx2FiOf34lb6l/2/n7c9/KDIc2odDMy6ElGizkx73eK4xE2nmYB3uUPa6m4hMerUklQLbXC4vBQWy5IPQ+S2RxP+mfOF4nzfSMu10olRSnRt7pG+srnzGaQkiMmqLB8ck2KzBOuA3aouzdMVIur6sY2Q87SCSaQVJsylJzmePiIwlakm6mLahQPaypF3zlAcFHpDmlZWczXtA6nCc/cWbA/96Q8vh077n3x0R/O/0eSNmg2amZp1RWxxPTIlvevd6wjMlbjFS9ZBy1t0kUaWQPZP4fEHSZGp0O4Pzsn56RGTUFovPdLIO+LNkEjOWcmk22rEO+KefjJWK7cMt/gO0DeXonZN0Aamzfyz7MJpUbRSbM0efjyQabw9GJ/tGcwzZ6TSqSHwundFX4AFJwsVu3P3kQP8jmxdd3WT6wabGn77aL7ghmZ9R75ylmPRqSaZ2sqzgS2qslCABiR0xkKPulmjAs7N/rC1r//J8WaLrbmeOaQki0qlVkky39Hkj2YfsSKLswaGxaM5PT5LBVbb3dvb+ZX2FVpJlpF2DgSwpdowkjTybM5BlVR/TWKVHD0lB8j+Eevi0r88b3ra2joi+cZnFF47bhgQu1D/W78+y5I3ObdwpquntcI+bdJk3DUolyVK+nJGPiOqkGIe0DQVzB4nyEk9QmoS0nHGi3qQVX80NjUVzVnPsCA+RBfV6wia9OksmMWMsKRYf/OyuUM42iiQ9JDa0lb3uJom6mLahQM6h6epSjfiRmz5vJMtqJ8ZSViLJ4GreFHhAknalNxft75254dJa9uirVUXXLTcf+Mgr7J/q7J8xo4GprxDbn7ANBXLWp3RuhZCoOBGMTvrC8ezDaERkKSsRPwsyEozlHEaTpOvscI83mfXZp6+JyFImwT7c1oEZdyBNKShHFOHC5sy20DKp3qQTP8BlGwpk2aCIkaQz0esJZx+ElKosdj5IzudckgW/PWfGF1fn+PTYFOYcyjQu9ICU3xNB4ompPV2uO65ckLzypSWVBz8aFfavDY5FslcNZoNGZIOubzTHPDmj06hEDtmxmaqcTVRJtqy2DvhzTlaRFLWP3RXMOX1NRNWlYn9NRGQd8Odsdxu1xZEJ0SNOA4Gc3SMiqq/QiiwrnpgaHItm2U2DMek14tuUR/vHsmfMMzq12Nk+Nu+bc8hBmiG7XAmKzNw6OLvAA5JaVZTPQ0EOnPRaykpS/8Y2Lanq7B8TMGLjC8ez7EPKWMpKhsQdZN7jHq/jVvuIr7uzrPxPEp/jy1aqZtl4Ikn8DNyxfn/25HKmsUoncsiuzxuJJ6ZyxglJUsxtzuBSDm2UhZU6kc8e619mXwZA549fElMQEX10Jnxx1sEGxqRXB2OiHj+7K/dUH0m0o0/OVbGM2VAyh84pLfCARPmdRnq5x7upuSr1CjvWQcANdA8HPmfSZc+IrRCdF8TOb835Y+Irhez7Wn5akOhcA3aye/a5N8ak14jc16Cz38/tTYn99KyD/pba0px1N0kRZT2hWJZdBiQsKGfK8rmCKrSRuNicGrsreHmuGTiSYg6payB3RgNJ0fBiM0Nc+rIS7kWZB4UfkCxlUp4Tk93ebveNbfPSLrYtKOvsH+P7T707GLwo1wBxY6WoGbJIPNHrCXPp9WuLVZNTouZFO7PuHp0kfgbbNhRonZ9t44nUskR2+wbHIlyqVJ1G7MnFLn+MS7uBRD/t8cSUwz3O6U2pVSK7Yt3OYJaNDVOJrFI5Zu4QUV25VuQUJpecbzo/hSnmqejzhnNmNDDix9vzqfADkrQnaWZhd4XUxUXTZy9a5pUe+8TP91879knuGlzkagb7cKjJrM/eCWPEN4cHx3JnBJEUq+U9oQku4xgkehw/GJ10+bPtAZhSkNgmEcf+JbElqyJ+Uw73eH2FlktXTHyn2ebMnQzJiHwqXIFoztQJxlJeIrLu9oRiZkPuyEdSRFkuIZakiLL5VPgBafk8Q7d0J6RlceCkd/0i0/S2+eJq/eGPfXz/tUFfjk0niQ1wiRhhcLhDXMa7iUjkKshIPOHyx7gMo5HozTWyHJGepkKnFlP72JwBjuHcUiZ2IIhj7gmdX0wjuKBBX4Rjqp6xRC0yyrr8UY51t8hZEMfIOMdPT2SQyHnoSSqRHRce/UspJuHypvAD0uqGcusA7w6KAL+zjXx9Rc3062sWVrCJDV4c7hk3DUpiw+uCO/4sI4jLT9abRE3tcO+KkeiOi20owCWfikTv7mwfDtWbOBVk1BaLTmoIW7Lu5ZNkKS8ZEdNpdoW4pEeT6OmWYHSy1xPm0r+k8zsVCS6razDAvX8pKvLlOvQklaVMVJS1DQW5TCARkU4tdsQ4nwo/ILXON/R5Ixy7/K5A7Pl3R14/Ncr3VxiMTtqGgte0VE//VkttqU6j4lXJsj+/3IsTtcU6tUpw84fL1gmMpUwrZhmpzRngWJ8SUb1J5xKavhWMTjrc41zHZ8R1XLoGAxzbwiRuvorlW3Js42uLVSLrbo6PBGsMCR5JYyl2PNooIn5TdheneR1iXTERBXX2j3FJJWWMWrWohpcz2xmAqebWdnaFH5B0alVjFaejg/716NCiHx/9yYG+v3jWvvKJt3k9LodP+7Ls2N9qMdiHeQwbeoITHAeIxSy0sjk5rYqlc3vSCJ+D7XYGOdandH5tubCC7K6QUVucc1kiI3KNWt9ohGPdTeKqVLsr2GTWc0nTIAmmFTll5zNqVZHgAS6O6Z1MXYXwWRC2DIBLmgaJe/aIzzAasfwdoY28wbFofHKKY5S1lEuzS1Z+FH5AIqKNTZU5t0t4aF/vE4c+efPey95/6MrhR9df3VR1y28/4D7I2/7emS3LM4zXMS21BhufeSzrYLZjNFMJ3kOBncvHdWy9vISIBNd0Dvc4l1xYRkyinScUa7UYudbd4nIN+rzh7IdBXFCWiIMiuexA+mlBZSWC16uyYTSO1Ryd21BK4AfYPcw1TYPERb6cB1alYiPYgp9zT4hrO5KVJXghl8sf495GkWp32vz4TASkryytOnI6W+L1z98YeP7dkdfuWcVy5NSqol9uXVKqUf39n05zLKLDcTbjeB1jKSvp4TONxD2rSqcROGRndwVz7sz4aSlqFdtiVUBBROTyRzkOLxDRvDLhsyBc9tdJEjODzWsYjcTtVNTtDF7KOfKJmf/3hCbMBg3HYTQSF/w8wQmO85dEpFIV8frzSdXnDXNZvJxk1BYLfirsw0Huz7mY/Tt49S/Zcum5Mo30mQhI17RUHz7tm+k52/+B5ycH+l7ZtjJtzuZXW5v/9a0hLinjfd5IJJ7tyOflFoPNGeB+w9YBP5fFiSQiq8oTmuCY48SYDRphNR2XTfJTiYkTXYOB7GfepDJqi+OTU8IynTyhCZNezTGck7hMDU9owmzkWnfXm4Rn2TncIY7Tb4zZIHyAK+dhg6muaCgrFlpR8epf0rnHT8hvis1fcn/OxUS+EyMhLhusJEl4GK7cPhMBSadWrVlY8cw7w9O/5XCPf2tvzwu3tU5v7TZW6e68csFPDvTl/Pf/YD+zqbkqS/XUVlfGqzXEfdS7rkI7FhEyscxrTp6I6it0wv5+HO5xdnwOx5/XqYWnww6ORTgmlzOCD6h1uEPc28JEVG0Qvrvzx2fDl3AfGywrEZxoYHeFuNenJKLpwGu4mMTlGnQ7g9l3KE5TX6ET9vg53OM6jYpX/1Lwb8o64G8y82o6zJlRu89EQCKiv7mq4RdvDKQ9aoNj0Wt32XZe37RxcWXGVz2woWFPlyvn7/JFu2fLcnOWH2gy67kftcIeUz7Jo0L+Vh0j41y2LPu0oHKBfz9cdgVNJTjXIBidtA+HeDWHjSXFwvYuO/aJn1c4F1wjeEITn/giOTd1TlKrigQfNMc9PZqpq9CGBH16vIaLichsFJ6jyHHz8k/LEpoPyeWg2FT1Jq0vIjBI2IY47ZiXZCxRB2MSnD2YB7MQkI4fP37//fevXbt2y5YtO3fu9PlyLxptb2+/9957161bd9VVVz3wwANvvvkm30I3NVc1mUsff/2T5JVeT/jqfz5+9+frblqZvtlPUmOV7qaV8548PJjlX47EEzZnIPuQLvdMPyKyDvg5rs8gEctBbM4A9wlYErFavnuY64L88wUJnAVhw2gcF2cwlnKtsBRzx8g49xQ7ErF3mW0owD1NgxEckOzDXPvljOAoy30AgGG/U2GzILahIL+mg1EjbHTLFYjx7F8KXNXOHiTu/UtCDymL9vb2m2+++eDBg8uWLYtGo7t37966davH48nyknvvvff73//+a6+9dvHFFzc0NLzyyit33nnn448/zrfo39y07Kkjgzte/Ojwad/3/vvU539pvW31/L/94sLsr7rjyvn/9tZQlj8G21CwyazncNKPlmPmt83JK6tKK+CIYl843ufNcbZFmupSjbBcA09wgle7W/AZCtw3nkgS3L/s7B/jPgUipiCOe6OlEhyQeE2BEFtGKqhKHQnGeE2BkNCOiy8c12lUHAcbmLpyrbBtG0+dDfN6zgV3xdgWIbzaKFKdbZ8HeQ1Iw8PDjz76aHV19UsvvfT00093dHRs27bN6XTu3Llzppfs37//lVdeufjii998881du3bt2bNn//79tbW1u3btevvtt3mV3lile/dvrghGJ69/+r3Ofv++Oy79uy835nzV+kWmJnPpH+xnZvqBF+1npm+oOp2lvOTESIjLfXJfr0pCV72xzgSvHpLg1fLcV/AxIv5Q+XXFSMQOobymQIiovkInLEP6bGhinpHHr4mE7pjuCU0YtcW8qjnB7e5eT5jXFAgJDX7WAT+vXxOJ2GiH4+blSWxiVcDjZx/m3fCqLkUPKZP29vZYLLZt27aFC8/1S3bs2NHQ0LBv3z6/P/PuPn/605+I6Ic//GFV1bljHS666KLvfOc7RPRf//VffG+gvkK76y9bRn/8hdfuWckxjY2Iblk179/ecs703edtI5uWVM303aTl8wyOEU6pq7yebGFxwuHmN9dCQtv4wehkn5dHih0RsWlhAX8/Z0MT1ZwziRlhLcfBsSivFDsiUhcXxSeFtLvtrhDfulvYb4rj4ToXFiRwHZJ1gOsyuySzQSMg+c3hHufVGCIinZrfpipMJJ5wuMd55ZeT0OS3EyMh7hnzjFFbLGy2L//yGpC6urqI6Nprr029uHHjxsnJyZm6O2+99VZJScmqVatSL15yySVEFAjwSKQW44ZLazv7xzLmf9uGgjq1ikuF2zKvlGMl6wlNcH+yWTuL758Q34F1OpdPLKRG4NvuJqFNb5szyLfunmcUcmK6bYjfPDmJWIckYBxSp1EJGPO0u0J8OxPCDixmbRS+VaqwFPNuZ5DXMAAJTfu2D4d0GhXfstSqIgF9WeuAn9f8JUlx0kre5DUgWa1Wg8FQU3PBjgYs2Jw8eTLjS372s5/927/9W9rFDz74gIhMJpM8t5nObNDcsKL2F28OTP/Ws9ZhLuN1RNRSa+Dy8LFkPF5PtoAjirudwXkC/ngEjZnwrbuJyGzQDPp4N735zsmT4MjHP5yToI5LJJ7whCY47oSUVFeuFbAFAPe12Els01i+yRo2Z8CoLeYfkIQEP5uT69Z8KQUJ6TQ73KG2BWV8G17CllK4AjG+3T6Rh5LkkzpvJSUSiVgsVldXl3Zdr9cT0eBg5ky2DRs2pF3xer0/+9nPiOjGG2/MXuLSpUtTv+zp6eF1w6nu/nzd1t+8t/P6ptRnLp6Y2tPleuXulVz+hSaznh1Enf2pZTU4v+Ggc+0sHrWJzRm47XIL958nosYqvYBGFt+5FsZs4F13s0DOt+4WPGby58tm3JVjJizXgFdTo7N/jPsOpEmW8pJj/by3t7cNBR7Y0MDrJax3ziafuL+KZdPwrbuN2mIBW7O7/DFeuSdEZNKrBXTFXAFOJ2OlEdAeYie5cC+L1YHhqsXuS/6fHOlbypC/HtLk5CQRJWePPr0DlSr53Zw6Ojq2bt06MDDw4IMPsoG7LHouJOiuz1mzsLyl1rCny5V6cU+Xq6XWwHGCRKdWmQ25M0qtAwG+NbiAUNHnjQhISCMivnOwfFOPGAGL2K0Dfr7p0SR0/t8+zG8VCFNv0vHt9g36eGy5lCSs28droDip3qTlu5BLWBtlnpH3MrhgdNIXifMty6RXE/8pTF6JSElGbXH/KL9HgtdJLnS+Dnx1X3vlgka+tzcr8p32nUgI3PzqxIkTN9100wMPPBCNRp944olvfetb0t5YTj/YtPAnB/pSa+QnXv/kjivnc/8XTHp1zr/ebmeQ+z6kjKWM3zGXAtbr0PkVl3z7E3xTjxgBi9hdgRjf7hEJrbvZfp18XyVgGqnPG+E7LUGCNnnyheN8B4o/LYvn8GDXYIDvE05sGSnPR8LhHue1/JYx6dUC9gpxuEMcT+FKJeDx84Ri3E9ySWLH30wVzYFtEPJ3ixqNpri4+MSJE2nXJyYmiEirzfYp79y582tf+9qJEyfuuuuuP/3pT9ddd52MNzqDTc1V9RW6x187t7R273tuIvrGZTwGvlpqcx9C4QpwPVk1yVLGb2ZeQI4TI6CmE9buFrAf0qmzYQEBScCSHU9oIp6YEhCQBJQlrN3dWKXn227gdcJFKgG5BvZh3okGJOg4O+uAn9dujUl89wNkB8UKiBPzeDYlSVDGPJ0f3pjUcd2PdRblbw6JiGpqarze9GMgxsbGiGjFihUzveqee+45ePDgV77yle9973vz5/PokUju/95y8aIfH2Un7tz8f0/su2MFr7/hJrM+51IkvpukEdG8spITLk4rnBi7KySgRiCixio9r/StobHo+MSksNqn5wy/3Z1tQ0G+s2IMixNsoIZbQby3TmAqdLwDks0ZeOQri/gWJCClyjYU5HhyeZr6Ci2vDRRYejSvLVwZSxnvDaVOjPDOmE8pK0rEtSngcI+ri4sEjENaykoO8h8bFDAGrlYVmQ2aKRXvVlT+5bUT19bWFovF0qZzjhw5QkSLFy/O+JJnn3324MGDt9566y9/+cvZjUZEVF+hfWVb2y/eHNj6m/f/838tz3LeREbVBk32pUi9nnAwyrsGZ2d3cv/5bmdQQGeC+Nd0p73hRdVC2t0ClgfZhvgtv/20LJ7JbwK2TmDMBg2vjmwwOilsuoVNMPCKE7xOUEzFdxmpfZjHCYqpzEYN31x2h3tcQN1N/DeNdfmjAsYG6dy2jTzHwHluzZdkNmjiWiHjInmW14DEViDt3bs3ecXv97/66qsWi2XlynO5aj6fz2q1dnd3sy9/+9vf6nQ6thJWCTYurjz9g7VTT3zxhhW1fF+7fpGp92y2gGRzBgTk6vCtFGxD/PbQTDKWFPOru4eCn+d8aksqs0HDaz8kXzgejE0KrFJ55jUIa6IS21mczwmhwqZAGL5Nh77RiLA3pS1W8dpoxxPicbJcqvoKrSsQ4xVlXf6ogPlLIjIbNbxG0no9YWEFWcr4TYyxsUEBSw6IqN6ki2uFRLI8y+uQ3ebNm5csWfLcc88tXLjwxhtvdDqdDz/8cDgcfuSRR5I/Y7Vat2/fXllZ2dnZ6fP5+vv7Kysrv/nNb07/19auXXvXXXfl8fbFaqktzb6/ateAkGa+sUTNa8G8JzQh8O+nnN/ReYLr7voKHa/61OEeF1bNEf9lpK5ATMD0NfGfwXYFoia9wDEW1sXknrfi8keFTSvWm7SHTuXeHDmp1xMW1sCn83socHxTkXii1xMW0LwjonlGfj2kEyMhYc853w2LeZ1+m6a+QmtXC3lu8yyvAUmlUu3ateuOO+547LHHHnvsMXZx+/btW7duzfjz77//PhGNjo4ePXp0+ndra3n3UWYXO3Q1y1IkVyAm4MnmNcDF8qmEDdlVl2q6BnnsjtE3GrmxTcjvyKgt5vWH6gnFhE1fE//ktz5vWMD0NREZS9S8jhtgPSQBBRGRSafhftwAO0FRWFl8u2KC2yh0fgE4x4DU6wnzOv02lUmvPnWW08b8jMsfE5A3SERmgyYYnYzEExzv0+WPCUgRSpY1WSLwk8+nvAYkIpo/f/7+/fsPHTrk9XrVanVbW1vayqRNmzYlJ5k2bNggcv2Q0jRW6Rzu8Zk63Q53iFceOcPmkDjOzIvpTPBN+3b5o8LqbrNBE09MBaOTHBuD7E0JKIjO7QrB9U3xPbk8Fd+JMQELAJJMejX3Nr7DPW7S8duaL4nv1ux2V/BLSzKfPZaTpYzH5lW2oYDguptvt8/hDn23VuCqU0tZSZ+Xa2uA75FLqeaVlcRL5sAcUr4DEhGpVKqrr746/+UqQUutweHOPAocT0zxPScmSadWsR0/c/6kKxAVlk9FPNO+xbS76fwAF8eA1O0MXrVY4D5SFXp1j5trRp/dFeS7rWoSO2w+51YdSQ53aEtrtlMfs+A1X2VzBgSkvZ0riOcMnN0V4nVaYyqdRsW9Lyt4Xof4D6463OPCxgaJyKRXuwJRjn8mfE8uT2UpK5kTad9zYKlUIWm1zLjnt304pC4uEjhAzHnnU2H7sDGWMi33yV6He7yxUies7iaeyW/Clt8mC+L+puzDwutT4nNWCJu+FpY3SETzjCXcB526heYNEs+6OxJPxCenBJdVX8Fjn6dTZ8OLhXaaeUXZPm+k3qQVNjZIRGYDjzXg7Og1YQXVm7RxncC+aT4hIOXVcouh25l5bWxn/5jgWrWxSs/xse4aDKxbJLChZCnnEyScAWHz5IxOXcwxU0PM9DWd61xyTQk5MRL6ytLcR43MxGzkWn2zz1lw/9JSVuLg3O3zheNLa4R3ZNksCJcftg+HWucbBLdRzAYN94w+u4v34VhJRi2PbNLDp32r64Vnr5kNGo4HSRORzRnY1Czw8VOripD2Delaag0zZX73j/I7NygV96b34FhE8JAdm0/mWKX2uMeXC+2KEZ9BJ09oQvD0NZ1bC8K19vGEJoStKWa4N71dfoHp0UxjlY57u1vA4VhpZWU8mWU6vkc1pplXVsKxIBKXzsc6zRz7zX3eSMs84Z1mk17NsSA2qyr48auv0MW1GLKDC2WZ2ba7QusahfZduC3li8QTjhHh493EZ4DLFYiJq7tLOJ6h0OcNC56+JnYGPOdp+V7PuMghO47hXPD2ToxRW8wx8kXiCftwSFjWJcMryoopiC1F4vKTkXgiMpEQsL0Tw/Y14FjWSDBWXSp8BwTuu2SJydwhothkIjEX0r4RkPKqvkIbjE1mrJU6+8cE10F1FVoum/OL7EwQn1G7Xo+ov595nPt89uGQ4KFOOj9ZxWXFpZhliQz3dPYTI2KDBMc1pI6RcQEny6XintHXPxrhexx7qnoT16YD2xJbcEHEhhy4DeRaB/xiHgnuc6XsyCXBBTXXlBbRlPKP6UNAyrdWi8HhTt96zuEe94XjfHfgTrKUcUopFtmZICJjiZpL9kQknrANBcV0Jrjvh3RiJCTmTRm1xRyPL2OzMsKyTpi6Cu0Qtzwx21BQTN3N8Xkgot6zotoNxGd1p90lPHuCiOordByDhMMdarWIWnNTb9JxDX6uEN8jly4siOvpml0DATFjg0SkHvcIOHswzxCQ8m11Q3nntPPTbEOB1vkGwZWdpZzTWJBtKChmtoBY6i2Xuntk3KRXi2l3cz86zzYUXCho64RPy+IW/NjvSExB3IfsRM7rGLXFOo2KS1mOkXEx8zpEZCkv4Xioj90lPG+Q+OypccIlcOuEJI4dl15PWKcW1b+0lGk5Lpd+8/SY4NwTRh0N8FqaPSsQkPLt0gXG6Yl2faOR9YuEn8huNpRw6bh0O4PL54mqfRqrdFxqn96z42KajUSkUxf3eTllH9mcwnO+GUuZlsvReX3eiMi6m+OQXTwx5QlNiBkIIs6z5WfHJ8TknhBRXTmnNn6vJ2zSqQXP69D5cM5lKNLuConpihGRUVvM5dMT2W4gPkezD/ujy8T9ptRRP3pIkK6ltnR6stCxfr/gpZ3EeQcXkWMmRDTPyOn8m15PWHDSLcMx+Y2lHglOj2aM2uI+DlF2JChkY6dUZoOGy6ATO/NGzNgg8XgkRM2KEdusgUv/UtwyAIZjv9nhFrVcjDgProrMuqTzvyYuUdY7PrHQJGokoLrnjxublL4UCQEp39hmDWmPoMMdaqoWlb7lC8do3iW8AAAcV0lEQVRzPtYix0yI81DGSCA2T9wfKjtoJ+egE9vwTWTdXW/Scun22YaEL21hzAZOmRp2V1Dkr4mITHpOwa/PGxb/SHDJxvaF42K6R4xOreIS/ERmeBLn57zHPb64WtQMnFFbbNKrc36Ag2NRlvsnpqySoEvkx5IHCEj5ZjZoTHq1ffjTvIbBsagnNCFyfqLepM3+WIvMhWUaq3S+cO5haLsrJLLXQkRccg3sw0HB294kzTNySgFwuEMiy+J4hkK3U8iZqmkaK3Vcgp+w49hT1Zt0XNK+e9zjS0U/Eo1V+pwDuWxskPuJixlxTB20u0IiEw1YWTnflH1YgjbKnICANAvWLzIdPv3p7o0djrPrF5kEr4Jk6it02R9r8bmwRFRv4lTNSfL3U1+Rez8kh3tcTHo0YzZoco5DskaDyDjB9sHL2Ry2DQUFb3uTZCkvyTnoxLZVFVl3c9zRJ8uewvzKyv2mRKW9nS+IU0afmD2rkrj0xuyukPj+5ZyAgDQL1i2qSN1OuNsZvHKh2LOz6k05/lZtzoDIXFjitjjRE5oQ3+4mosYqfc6yugYDItM0iKjelHuPPttQQPy8DnFb4NLrkaDuNpQU54x8dldQfN3NMaNP5FJfZmFl7pwaMbs1JnGZrGLHQopvD9VX6HI+fl2DYjM85woEpFmwur7cOuhPDt10OM5uXCx2sjHnEE2Pe/yyBrE1gkmvzjmaYR3w11doRba7iaixSselOSx+bDBn55KIPKEJScZMuAQ/VyAmMneLiBordTmDRLdTgrqbOOwe5AlNROIJ8eOQOcelSdyRS0lcdg8Sc5JLKi6Zqy5/THzDa05AQJoFLNWNLbR8/dRoPDG1WnSoqNCrs+9B0tkvakl5Us6FkK5ATHy7m4jqyrVjuZKkXYGY+HY3l8jX4x4XmaB4vix99iqVTV9LMIdUpc85tWMbCl5WL8GGm2ZDSfaZRdZuEF93c+mgi89oYHI+FWJOcklVx+FN9XrErqOYKxCQZscNK2p/Zxshot92jXy1tUaCdlZljlaqKxCVpo2fqz9x6mxYkshXb9Jmz8YWv7SFYdspZW8O20WvtWTYqUhZfkB82tu5gjh0Jjwh4SeQXlBWRY7flCsQk6TuNmrVOXNqej3j4vNciMMeffZhCXrnxGEnXE9owheJS1KW8iEgzY47rlzw8zcGfv7GwN733H+z8XPi/8HsC3fiialBX1SSlmPuySqJ2t05s4/Ep70lWcpLsleptqGAJJ/ewkrdSDBbc1iqBj7Lm8+eJN3rCYtcr3OurFwzLt3OoCRTIBzHBsXP6xCHJasnRkLi8waJw5lStqGA2aAR32adExCQZkdLbekjmxftePGjX21tFryFXarsHReWjSZ+Tp6ITHp19gadmNOjUzVW5hjdsg1JkB7N1Ffoej0zHiDkCU0MjkUlqbvNBk32N9U1EJCkK0a5OkmsVSHJs1ehyzFcbBsSu0UIk3OzBuuAv9VilKTubqzSZW+j2Ielefwaq/RZnj06N1n1megeEQLSLPr2xs9NPfHFb1xmkeRfY9vZzfS3ahsKSDJmQrkynXzhuCsQk2TQie0snqWN3+0Mik/TYLI3vR3uUH2FVqraZ/rWuqnEb3uTZDZosrRRbEMSbJ3A5FwbaxsKSDUFkv3oLwk/vWqD5uzMHZdgdNLhHpfkA8w5MdbtFLsiew5BQCoQOrXKUj5jvdDtDIpfwcc0mfVZzri0u4KWshJJ6m5LWUn2VTuuQEySySoiMhs0Wc78lipNg87Xp9nb+FLFicbKbL8p8RtPfFpQVY7+5ZSI02+nl5Ulb77HPS5V/7KxMltSgyc0YdKrJelfUq59njyhCanelPIhIBUOtilRxm/ZnNLM6xBRk7k0e+0jYTpQY1X2cUixW1sm1WVd3dk1EBCfBsnUm7RENFNZkXhCzCkkaSzlJVmibP9oRJIpEDo34DljQZ7QhKa4SKopkCZzaZayJFkAx2QfSXO4Q1K1GyjX4KqEbRTlQ0AqHE1mPUsln67PGxa/pDxZSq8nPFMbv2tQsrqbiJrMpTO1HKVKjz5fULbaxxWIiTzhIkmtKlrdUDbTm7IPS7AVTVL2KOtwj0v16VnKSyLxxExVqiRb8yUtrtZnibLWAb9UbZTsIwHWgYCEw2gttYaZHj82f/kZSbEjBKRCsrSmNOPfKpvXkWpsXadWZRlhEH86Uaoag+bkmcyVgssvTcoy01JryJKpIcny26QsTW+Re+ymaTLre8/OGGVtzoBUdTc7E2imkTQJ+5fEWl0jmd/U4FjUFYhJOIxGRDM9593OoMjtg1NlCX62oUB9hfYzkmJHCEiFxKRXe4KZD0dn8zFSFZRlJM064Jew7r5kvnGm7eysA/7V9dLUp0TUWKWbKSUknpiyDQWl6l9S1jRFCdM0iKil1pDxeaDzu9hJuPdzkznDoSqMdVCyvEFinYkZoqxtKCBtT6KldsbR6c7+MQl7SFm6fTZn8LMzXkcISIVkpjmkzn6/tFsFm/SZN2twBWKuQEzCuvuS+YbTM1Rzxz7xS1jN6dSqepM244CnfTgkVcY8s7ha33MmczV3+LRvtURTfUTUZNbPFGUln5ZoqS2dKcrahgJrFlZIVZDZoJkpyoo/EDnNTFE2npgaHItK+AG2Wowzzf4OjUXXLZLs01M+BKTC0VJb2usJT8+TPnU2LO0zPdNkVZ830moRfhD7dFn2pJFwtoCpr9C5/BlmXDr7x6QtqKW2dKapPtuQlM3hLFG2zxuRagiXWVpbmjHKBqOTEqZp0PllpBkXA/SPRqTK3GFm6suyU7gk3H7bUj7jASjWgYCEbRTlQ0AqHEZtcZNZn3rSEmMbCqyXNCAtn2eYfgo7EVkH/OsvEn7u7XRZEigcEm0ul8SWPU2/3j8akWTl/6cFlWltQ4Hp112BmLpY7CFsaRqr9PbhDL+pEyMhkSfLpWm1GGxDmR6JQb+021QbtcX1Ju30h5yVJWFXjIiW1pRmfM4Pf+yTtn+ZJcr2ecNSbUcyJyAgFZS2ujKb84LKLp6YcrjHpW3jr24oz1il9o9GJMxoICKTXm02aKa38R3ucbYBnYRlXVZfduwT//Trnf1Sjg3S+dGt6WtcOvvHpK1PichSVnJiJOMo7tj6RVI2HcyGErsrOL3p0NnvF3/oSZqMi4t94TgbXJWwoLY6Y9pfE3PsE3+jpM+5Tq3K2JR0BWKSnHAxhyAgFZSlNaWnLkzXsQ0FG6t00tbdM40Ndvb714g+2ClN24Ky6W38w6d90tanNEMbP56Ysg76xR8OksqoLW6s0jmmxYkjp8fWNUockK78XPn0N9XnjUQmpNnwLal1viEYnZw+43Lk9NhyiRYvJzVW6qY3HQ6f9tWbpNlNI6nVYnS4xzNF2THxB5ilmd6UJDan2FD+2UmxIwSkAtM632AdvOCx7nCc3dRcJW0pbNgkbY1LMDppcwYkzHxjLl1g7BpM/0M9dMp31WKJA9LqhnKbM5BW+1gHAktrDJIf1rm6oXz6Pmmvnxrd2CTxm2qrM07Ph2RdMWmrOZ1albEsuyso7XAxEV1WXza909zrCUvbbiAik15tKStJ67iwFVeSt4eu/Fx55udc0jFw5UNAKihtC8oOf+xL7bsc/Gj0K0slDkhE1FilT2vQ2ZzSnKk6raAMO7hYB/yS1whmg8Y47ZRVhzu0TLqVqklXfq48bXIinpjq9Uhz8ESqtgVlvZ5w2skaJ0akXFaVWpbtwjflCsR84bi0w8VEtP4i0/Rs7K7BgORtFCJqMpemra9yjIzXm7SSt1FaakunD9nZh0MSLuGaExCQCkpjlc6kVx846WVfRuIJmzMgYR520rrGC05hp3Nr16UvqK2urLN/LPUKm4BpMks/09tkLk2bnDh0yifHNmLT93K1DgTMBo20I06U7MteGNEPnPTKUXdfVl925PQFv6kDJ71NZr20w8VE1FJbyvYvSL1oHZB4bQOzuqHMOnBBw0vyrEumZZ4hbRIuEk/IVJaSISAVmusuNidDxeu9o20LyuSYFN3UXJkMe8x/f3hWjq5Yq8UQn5yyuz6NE/s/8KxZWCF5V4yI/vzi6tQoG09M7f/As2mJ9G9qzcKKw6cv6Mj+94ee6y42S14QEa1uKE+N6L5w3O4KbWySeHSLiFrmlaY1HQ5+NCr5MBoRqVVFTWZ9alkO97gnNCFt1iWzfJ7h5Z6zqVde7vHK8ZzXV2jrKy7IHtz/gadlnpTJ5XMCAlKhue1yyx/sZ1hT6xdvDt52uTTHW6RZs7CizxtJHizmC8c7+8euaamWo6yNTRcEP5kiHxGtX1TR4fi09rEPh4zaYjmqOUtZSZNZf/jjT4OfdSAg0/rHLy6pPHByNPnl66dGWy0GyXstRMQGUVN7fgdOylJ3E9EtqyypvbH9H3i+2lojx5v66iU1dlco+ZxH4okOx1mZmg4bmyrb33cnvzxyemzL8ho5ClIyBKRCs36Ryagt3tPlOnrabxsKfLVVlmdarSpqqzMm40SH42yrxWDSq+Uoa11jRbL2iSemDpz0ytFrIaL1i0x9o5HkaqQDH3nlaOAzf77M/OIJD/v/YHTy9V5ZOhNEdP3F5hdPnEn2xl7u8X59Ra0cBalVRRubKvd/cO5N9XrC6uIiObpiRLR+kWn/B57kANfBj0alXRKbpFOrWi2GZG/s8Mc++XotVy02Hf740yjb4TgreT6I8iEgFaCd1zXd/vyH6560/vT6JpmCBBHdtnp++3tn2P8/ceiT21bPl6mgGy6t3f+Bh8WJ598daazSydFrYdYvMu3pchFRPDH17DvD29fVy1TQxibTH94/15F95p3hGy6tlamas5SVNFbqWNPBF47v7Xb/1eVy/aa2LDf/9riL/f+/vTV0U9s8mVKW1ywsZwvsiMjhHj982nfDpbJEWSK6sW3eb4+PsP//6Wuf3HHFApkK2ri4srN/jPXG2O9LpnCuZAhIBWhTc9U7D17+1v2rb5XoONqMvnGZpaPnrG0oeOCk1+WPSXX07XRmg2ZTc9VTRwbjiaknDn3ywIYGmQoiou3r6p46MhiMTj55eFCnUcmX47RxcaVOo3r+3RFPaOKHHR/fcaVcQYKItq+r/8mBPiL6YcfH6xeZ5JuWuOHSWrsr1Nnv7/WEf/7GwB1XylV3E9EdVyz4+5dPE9FD+3rvvHKBfG/qppXz9na7+7yRw6d9tqGArM/5TSvnPfH6J0T009c+ue3y+Z+pFUhM0dTUjIdXzmlLly7t6emZ7bsocI+//slD+3qJ6D+/sfymlfPkK8jhHr9k57F4Ymr9ItNr96yU9Q/12l3dh0/7iOiVu1dKvs43VYfj7LW7uo3a4k1Lql64rVW+NxVPTK184m22Oc2b914m1TG7GT15eHDHix+ZDZq/unz+//nzxfIV5AlNLPrJUbWqSK0qev+hKyVPUEz1w46PD340+tGZ8Uc2L7p3vVydZiLq9YSX/J+3WL7lOzsul3ZWbE5UiQhIIMqBk16TXpOH1RLWgcCBj7zbPl8n3yAkE4xOPvPOcFudUfKlTtPt/8Az6It+4zKLHEmDqXzh+O5jzmtaqmWNRsze99y+cPyv5G/g93rCf7CfuWnlPAk3b53Jni6XWlUka6uLsbtCr/eO3nBpreQhdk5UiQhIAACFb05UiZhDAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARUBAAgAARZB342SpHD9+/JlnnrFarTU1NevXr7/rrrtMJtl3YgYAgHyaAz2k9vb2m2+++eDBg8uWLYtGo7t37966davH45nt+5oDli5dOtu3MPvwIeATIHwIc4TSA9Lw8PCjjz5aXV390ksvPf300x0dHdu2bXM6nTt37pztWwMAACkpPSC1t7fHYrFt27YtXLiQXdmxY0dDQ8O+ffv8fv/s3hsAAEhI6QGpq6uLiK699trUixs3bpycnHz77bdn6aYAAEB6Sg9IVqvVYDDU1NSkXly1ahURnTx5cpZuCgAApKfoLLtEIhGLxerq6tKu6/V6IhocHMz+ckxjEj4EIsKHgE+AiPAhzAWKDkiTk5NElJw9SlKpVMnvzkT5p8cDAEAqpQ/ZEVEikZjtWwAAANkpOiBpNJri4uITJ06kXZ+YmCAirVY7GzcFAACyUHRAIqKamppAIJB2cWxsjIhWrFgxG3cEAACyUHpAamtri8ViaRNCR44cIaLFixfP0k0BAID0lB6Q2AqkvXv3Jq/4/f5XX33VYrGsXLly9u4LAAAkpvSAtHnz5iVLljz33HN79uyZmJjo7++/7777wuHwgw8+ONu3BgAAUiqampqa7XvIYXh4+I477jh16lTyyvbt2++///5ZvCUAAJDcHAhIRJRIJA4dOuT1etVqdVtb2/SVSQAAMNfNjYAEAAAFT+lzSAAA8BmBgAQAAIqAgAQAAIqg6M1VhTl+/PgzzzxjtVpramrWr19/1113mUym2b4pKbW3t7/22mvvvvsuS/G44YYbNmzYkPoDXD6BwviUxsfHb7755paWlp/+9Kdp3yr4D+HMmTP/8i//8s4775w9e3b16tW33nrr6tWr034m5xuc058AEbW3t3d0dLz//vtVVVWrVq268847GxsbU3+gIB+DYDC4efPmxx9/fO3atdO/K8kvfbY+k0JLamhvb//+97+vVqvXrFkzNDR0+vTpBQsWvPDCC2azebZvTRr33nvvK6+8olar165dGw6Hjx8/Pjk5edddd337299mP8DlEyiYT+mhhx764x//+OUvf/nJJ59MvV7wH8KHH354++23j46OLl++vLS0lD0G//iP/7hly5bkz+R8g3P6EyCie+655+DBgyUlJWvWrHG73Q6HQ6fT/fu///uVV17JfqBQH4Nnn332H/7hH/793//9qquuSvuWJL/02fxMpgqI0+lsbW39/Oc/39fXx6787Gc/a25u/s53vjO7NyaVffv2NTc3f/WrXz179iy7curUqfXr1zc3Nx87dmyK2ydQMJ8S+zSam5u3b9+eer3gP4R4PP6lL32ptbW1s7OTXbHb7StXrly5cmUoFGJXcr7BOf0JTE1N/eEPf2hubr7lllsCgQC78sYbbyxbtuwLX/jC5OTkVIE+Bn19fbt3716xYkVzc/Prr7+e9l1Jfumz+5kUVED61a9+1dzc/Oyzz6Ze/NKXvrRs2bKxsbHZuisJ3Xfffc3NzV1dXakX//jHPzY3N//t3/7tFLdPoDA+JafTuXLlyjvvvHN6QCr4D4H9xn/1q1+lXvz5z3/e3Nz8P//zP+zLnG9wTn8CU1NTO3bsaG5ufuONN1IvfvOb32xubn7vvfemCvExYE3PpOkBSZJf+ux+JgWV1NDV1UXnt79L2rhx4+Tk5Ntvvz1LNyWlt956q6SkhJ3gnnTJJZcQEdsTncsnUACfUiKRePDBB8vKyn70ox9N/27BfwivvfYaEX3ta19LvfjAAw/09PRcc8017Mucb3BOfwJJxcXFqV/qdDo6fzxN4T0GjzzyyD/90z/90z/905o1azL+gCS/9Nn9TP7/9u43pIk/jgP4h5NRcVhQCjrRB/25B87SSSaYBEVirGIZYT0aP8uiIWRICaIhJBb8xMGkgqI/ZBTBptGzVBLG2EDQFGwx00BZpPlg6sx/0509+MJhd6vdg6u7+/p5PXKfuwf3fX+/+mG3706qGlJ/fz/LsqmpqRuL5M/358+fVbooJTkcjocPH4qKnz59AgDykaOcBChI6d69e0NDQ62trTt27JAepT6EwcFBlmWNRiPP84ODg16vNxKJiM5JOEBdJwAAZCPP27dvhcrExITX601JScnLywMal8Hx48ctFovFYsnMzIx7giKTrm4m9Oyy43k+Go1mZGSI6tu2bQOAr1+/qnFRChPtpgOAcDjscDgA4Pz583ISoCCl/v7++/fvX7169eDBg4uLi6KjmyGE79+/79u3r7e3t7GxcXp6mhRLSkpu3769c+dOkDFAvScAAGVlZWNjY48fPw4GgyaTaXV1taenh2XZBw8eMAyzGZaBiCKTrnom9DSkWCwGANLH3DEMIxylzLt37+7evTs1NXX9+vX9+/eTOxV/TkDvKUUikZs3b5pMpurq6rgnyBmgrkNYXV2NxWKzs7N2u72oqOjGjRs8z7tcrp6enm/fvrndboZhEg5Q1wkIzGZzSkpKMBgMBoOkUlxcTN49UL8MpBSZdNUzoeqWHQDwPK/2JfwLgUDgwoUL1dXVKysrra2tdrtdOCQnAf2m1NDQMDc353Q6yW/I79AdAgBMTU1dvHjx2bNnVqu1rKzs1atXx44dCwQCLpdLOCfhAHWdwJs3b6qqqoxGY2dn58jIyMePH+vr671eb3l5uXADk/plIKXIpKuYCT0NyWAwJCUlBQIBUZ28b9iyZYsaF/VXtLS0nD17NhAIXL58ubu7+9SpU6QuJwFdp+TxeLq6uk6cOBEKhfx+v9/v7+vrA4BwOOz3+8knrtSHQD7GT0pKEv1LsIqKCgAggSQcoK4TAACe55ubm5OTk588eWIymQDAYDDYbLaamppQKNTe3k79MpBSZNJVz4SeW3YAkJqaGg6HRcW5uTkAOHDggBpXpDzyZcDS0tK6urr09HTRUTkJ6Del2dlZAOjo6Ojo6NhYHxgYqKioYFn2w4cPQHsIDMOwLMswjOivQ05ODgAsLCyQlwkHqN8EAGB4eHh+fr60tHT79u0b6xaL5c6dO2SbD93LIC5FJl3dTOh5hwQAeXl50Wh0ZGRkY9Hn8wHAnj17VLooJT1//vz9+/c2m62trU3ajUBeAvpNKT8///9fkW3fJpNJ+BloDwEACgoKlpaWRDf0R0dHAYBlWfIy4QB1nQAZ5tramqhONrmQ27nULwMpRSZd3Uyoakhk77zb7RYqkUikt7c3LS3NbDard12Kefny5datW2tra393gpwE9JtSZmam9VenT58GAKPRaLVahVuXdIcAAIcPH15bW3v69OnG4uvXr8kh8jLhAHWdwN69e3ft2uXxeCYnJzfWyUdo2dnZsAmWgZQik65yJn/7m7f/UiwWO3nyJMdxL168iEaj4+PjNpuN47jOzk61L00BMzMzHMcVFhb+F8+jR4/W5SVAU0oLCwvSJzVQH8Ly8vLRo0c5jnM6nUNDQ319fbdu3eI47ty5c+SpOesyBqjrBNbX110uF8dxR44ccbvdwWDQ5/M1NTWRCnmYEMXLoL6+Pu6TGhSZdHUzoe3hqpOTk5cuXfry5YtQqaqqunbtmoqXpBSv11tZWfm7o2fOnCFPvJaTADUpLS4ums1m6cNVqQ9hfHz8ypUrExMTQuXQoUNOp5N8D4lIOEBdJwAA7e3tDodjaWlJqGRnZ7e1tQnfG6V1GTQ0NLhcrrgPV1Vk0lXMhLaGBAA8z3s8nnA4TP47g3RPPfXkJEB9StSHwPO8z+ebnp5mGGb37t25ublxz/nzAHWdAAD8+PHD7/dHIhGDwZCVlSW9p0T9MpBSZNLVyoTChoQQQkiPqNrUgBBCSL+wISGEENIEbEgIIYQ0ARsSQgghTcCGhBBCSBOwISGEENIEbEgIIYQ0ARsSQgghTcCGhBBCSBOwISGEENIEbEgIIYQ0ARsSQgghTcCGhBBCSBOwISGEENIEbEgIIYQ04SfQpQ158WV+iQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"iv = [1, 0.99, 1];\n",
"[t, xa] = ode45(f, t0, iv);\n",
"plot(t, xa(:,3))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"And this!"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAACXBIWXMAABcSAAAXEgFnn9JSAAAAB3RJTUUH4wkIAQgoaA7B1wAAACR0RVh0U29mdHdhcmUATUFUTEFCLCBUaGUgTWF0aFdvcmtzLCBJbmMuPFjdGAAAACJ0RVh0Q3JlYXRpb24gVGltZQAwOC1TZXAtMjAxOSAwMTowODo0MIM4q3gAACAASURBVHic7L17nFzVdee7uurUS1XVXZJKUoFKlmwJS9gKNJixGQKJmOHKxnfuh+jad2wuSWbywcTy4AfMXCdxEtvM4HzymRDbk4fvXA/wGULIHX8+QRN7wIxsA4ZrwIK0MMTCkUAK3VFJqparW/U41XXOqVNV948tilLVeezH2lvV3fv7F5S6atfjnL32Wuu31pro9Xqg0Wg0Gs3FJnKx34BGo9FoNADaIGk0Go1mTNAGSaPRaDRjgTZIGo1GoxkLtEHSaDQazVigDZJGo9FoxgJtkDQajUYzFmiDpNFoNJqxQBskjUaj0YwF2iBpNBqNZizQBkmj0Wg0Y4E2SBqNRqMZC7RB0mg0Gs1YoA2SRqPRaMYCbZA0Go1GMxZog6TRaDSasUAbJI1Go9GMBdogaTQajWYs0AZJo9FoNGOBNkgajUajGQu0QdJoNBrNWKANkkaj0WjGAm2QNBqNRjMWaIOk0Wg0mrHgIhikl19++bOf/ex11113yy233HfffdVqlf65S0tLt9xyy2//9m/Le3sajUajuSioNkgHDhy49dZbn3rqqcsvv9y27QceeGDfvn2VSoXy6V/+8pePHj3abDalvkmNRqPRqEepQTpz5sw999yzfv36J5544sEHHzx48OD+/ftPnz5933330Tz98ccf/x//43/IfpMajUajuSgoNUgHDhxwHGf//v1bt24lj9x9991btmx57LHH6vV68HPPnDnzpS996Zd+6Zfkv02NRqPRXASUGqTDhw8DwM033zz44J49ezqdzksvvRTwxG63e9ddd2Wz2f/wH/6D3Leo0Wg0mouEUoM0MzOTTqc3bNgw+ODVV18NAK+//nrAE//8z//8lVde+epXvzo1NSX3LWo0Go3mImEoW6nb7TqOs3nz5qHHU6kUAJRKJb8nzszMfOMb39i/f/8111yztLREudzOnTu536pGo9GsPI4dO3ax30II6gxSp9MBgH72qE8kEun/6yj1ev3zn//8e9/73s997nOsK47/ty+bnTt36i9Bfwn6GwD9JSyTM7o6g0TodrtMf//7v//7tVrt4YcfJnZLo9FoNCsVdQYpFotFo9HXXntt6PF2uw0AiURi9CnPPvvs9773vY985CMnT548efIkANi2DQCLi4svvPCCYRjvf//75b9xjUaj0ahAqYe0YcOGxcXFoQdrtRoAXHHFFaN/T5o4HDhw4MCBA4OPHz58+Dd+4zfS6fTLL78s7c1qNBqNRilKDdL09PTBgwePHTs2GM18/vnnAWD79u2jf3/11Vf/0R/90eAjlmV96Utfeu973/uv/tW/ikajst/wcmeVB80J+kvQ3wDoL2GZoNQg3XzzzQcPHnz00Ud/7/d+jzxSr9effvrpQqFw1VVXkUeq1erx48djsdiVV165ZcuWLVu2DL7C0tLSl770pUsvvfSWW25R+c41Go1GIxulSoG9e/dedtllDz/88COPPNJut+fm5j7zmc+0Wq277rqr/zczMzO33XbbJz/5SZVvTKPRaDQXHaUGKRKJ3H///du3b7/33nt37969d+/eQ4cO3Xnnnfv27VP5NjQajUYzhkz0ej3FS3a73WeffXZxcdEwjOnp6dHKJBR02YFGo9H0WRZbouo6JACIRCI33nij+nU1Go1GM87oalONRqPRjAXaIGk0Go1mLNAGSaPRaDRjgTZIGo1GoxkLtEHSaDQazVigDZJGo9FoxgJtkDQajUYzFmiDpNFoNJqxQBskjUaj0YwF2iBpNBqNZizQBkmj0Wg0Y4E2SBqNRqMZC7RB0mg0Gs1YoA2SRqPRaMYCbZA0Go1GMxZog6TRaDSasUAbJI1Go9GMBdogaTQajWYs0AZJo9FoNGOBNkgajUajGQu0QdJoNBrNWKANkkaj0WjGAm2QNBqNRjMWaIOk0Wg0mrFAGySNRqPRjAXaIGk0Go1mLNAGSaPRaDRjgTZIGo1GoxkLtEHSaDQazVigDZJGo9FoxgJtkDRjgeV2j1daatY6eHSh0mwrWGh20XrllKlgIQD49pGfW25XwUKzi9bMyYaChdxu71s/mXe7PQVrvXLKVPNLmXbn20d+ruZDPf6zSqlmK1gIEW2QNGPBF7574rI//LGCTWHmZOPm+1/9zN+8LnshANj30N9d9bWXyg1H9kLfPvLzff/1p//+e2/KXggA3vkHL/yT//S3Cozfnz9XuvWR1x762zOyF6q23Ku+9tKN//ll2QsBwBeeOLHvv/70ydcXZS90pNz83x78u//jL34qeyFctEHSjAWPHC5fu3Xy4LEF2Qs9c+Lcnu1rv/WTedPuSF2o2nKPzi/t273xr1+dl7oQAHznSOVDu9Z/6xXpC82cbOTTsd2F9HP/UJW91l/MnNmzfe13jlRkL/T4zyp7tq91u73ZRUv2Wt/6yfxN71733b+Xfp3/1eHyR6/YeGiursZvxkIbJE0Qv/Gtv//8Y8dlr3Joru52e5+7Ycvzb9Zkr3Xg787e/oFLCtn4kXJT6kLPvVm9duvU1cXM029I37ufOXHu3g+9a3bRkh2ieebEuY9ftemmd6/73jG5Z/xqyz1eaf3Xj1/+3JtV2QGu7x1b/MgVGz60c/0zJ85JXejQXD1pRH7vpq0zJ+tSFwKAZ06cu2V3fnchrSxojII2SBpfyg3nob8988fP/GO15Upd6Lk3q9e/M7dnx9pDc9IN0szJxk3vXnf9O3Ozi3JTVs+/Wbv+XVPXbMmWanLP3bOLltXuTm/O7MinjleWpK71/Ju1X96e274+dfSs3IUOzdWuKU5uW5fMxKOy15o5Wd+zY+0Htk4eLsnNjR2vLO3alL6mOHmk3JTtuMwuWtdunbpmy6QC44eINkgaXx599ey/eE/++nfmnnxD7nF4vuF84B2ThWzc7fakyg2Onl0q5hKFbPy9hfRr83I9pENz9fcVs9dsmTxypin1jD9Tql+zZdKITExfmpUddDpSNncXMrs2rik35Lpi5YazbV0SAHbk15Tr0tcqZOO7Nq6RLat5rdy8dutkJhEtZONSf6lSzXa7vW3rkldempF9neOiDZLGl+dnax94x+QHd657cU7uIeuV0+b05gwA5NMxqd7YkbJ5TXESALatS5aqcre5Us3aXcjk07FkLCJ1rWrLzadjALBtXfKYZGfi/N69KS3bazl2dmn7+hQA5NMxqXHI45VWJh7Np2O5lFFtydVeHik3r7w0AwDFqaRUi37kjLlrY9qITOzIpxQkxhDRBknjy+FSY3s+tW1dUrZO7Oh8c9fGNABsW5c6elbiga5itnMpAwAyiajsOGS54RA7kU/HpEbSjp1d2rlxDQBsXZucPSdx95ldtHJJI5cyilMJq92V6sseKTd3X5IGgGIuIdXKHq8sbVuXgvPXnlwre/Rsk1wSuZRRrku8p46U315ItngHF22QNL68udC64V25beuSUo1EteWWGw6xEzvyKanOxImF1vZ8CgB2FzJSP1R/7waAQjYh9Yx/9OzSro1rAKAwGZdq+Y6ebe7alCb/XcwlpAa4jp5tkm9v69qk1G+vv3cXsvFqy5W3fZt2Z3bR2pFfAwCyw4PHzi7tyKcAoDiVlJ0rxUUbJI03pZpNzsK7Nqal3jzHK61cyiCbwubJxJzkM/62tUkAKE7JNRLHK0tk6wGAXRvXSPXG+mttW5uSnYEj2xwA5NMxeTud5XZnF63iVPKthSReEicWWuRDGZGJQjYuz8Uk11shGweAzVOJUzITY7PnLOI0k/CGmjpcFLRB0nhz5Iy5u3A+r+N2e/K21NlzrWIuQf572zq5x+HjC+f37kwiCgDyPtTxSqu/dxcm46ekfShSPUPWKkzGpQaCTiy03vuWh7Rtbcp0ZDkT5bpDcvIAUMwlpMoUZxetrWuT5L9zKaNUlbVWuWHvyKeMyASQDyUzEmDabv/yK+YSyyiNpA2SxhsiSCP/vSOfkuckVcz2jvXnnQnZN+qgnZDqJB37+dL2txbalInLc1xIpippRACgOJWQe3RYtPqXRGEyLs+XLdWsbeuS5/fuqaTULEil6Qzs3RLTpaWq3Xeat61Nybay/bUK2bhsSSQi2iBpvDmx0CJhNADIp+PyBEj9mAkAGJEJy5V27m44SSNCMhMAUMgmKk2Ju8+2t87dxVxCppF427+EtxIhktaqttz+tyfVypaq9ra1b8cGK822vKDT7KJFRA0AUJyS+ktZ/es8lzLkLWS53UqzTWKDAFCcki4oRUQbJI03xystorsFgEI2Ls+ZqDTb/ZiJ1JtndtEiUSBCJhGVd/SuNNuFyfM7gtTdp9Js59Px/v8WJuPyzvjHK0t9O1HMJeR9e7PnLti7k0ZE3oeqNNv9qyKfji0sybKy86YzGBuUZ85nF61CNk78SyDXubTgKjraIGm8KVUH4jNZmcfhmt3fEaTup6bt5lKx/v8WshLXKjfsQvb8t5dLxqqWrG9vMAgJAJm4Ic/tKzecwUtC3rdXa7nbBz5UMZew2lL6GpRqdiEbJwFPANiUjUs9D/X9y1zKcDs9Sc0aZhdb/cMQvOViylhIBtogabypttz+cTgdl+hMVFvtvp1IGhG325MUnynV7OLUQHRrUu7u8/a5OxOrmLJ2BNPuTCWN/v/mUoYkO1FptnMpo3/uTsYi8izf7DmrH3ECgEI2ISkLMrvY6hsJACC9QmQsRNYauvwkKVBmF63+nQsA69OxBW2QNMsay+2WavbbOrFsfN6UmG4ZzILIO9ANLZSIRjo9KbuP5Xbdbq9/7paqnjhVG/725G1zg3t3Ph2XZ2VLVbufvwQAIzIh75LoJ5CAhKblnVHOWYNrGZEJScbvVN0mIlKC1PAGOtogaTwo1518Ota/rJOxiLzdh3Sj6f+vPIO0sNRev2YgZCftiFqq2oOfCGTuPrPnrEx8YPeZjNcsKfmqUu18YRBBnisG51V2b+/d8i6JIVfMiE5IcvtIye2gRS9OJSUJ7cp1Z/CMIu/ak4E2SBoPKs328LlbSSAIZIrfhjwkeVqDof0UZH6Bg6k+AIhOTEhaaFDzDW9FtyQFcktVe9DKymu0M3fOGnTF8mlZGsVqy00akcG1AECSla223L56AuSXUuCiDZLGg9lzrVzy7ZsnlzIkBfFL1QvyOiD3ONwaPONLjQ0OyvkAIJOISkrLzy5affUEyPxQ82+1dyKQM4QM40c6VQ8av81TiaYcnZjV7m6euuCMIsmcD+pBCPL6CJdq1uA9lTSi8jQ16GiDpPFgyJkoTiUl3jy5C/ZuiY7LcLIqLskVm120hqyspPhMteVabnfQTsizfOWGsylzQRxSkvEz7c6garm/OvpCQHQuI26fDPFbqXbBuQFkun3lhjOUGJMXb0dHGySNB3PnhrdUScWJox6SvNTOULJKnuU7Vbc3X/ihkrGIjOjWaCCoOJWQdByuttyhM34mEZVRxTwUhASpTvPiBU7zebdPwuV39GwzFbtgsy1MxiW5fdWWO3Sda1GDZnlTbjjpgSA+2fJkXNblhjNYMwEAU0lDRlqeaL4vTFbJqqSx2t1BmRNIUzqNfnvJWERiYmztBb5sPh2TYWXLDWfQSADAmbr9+s+ldDG32t0h+YmsuSS9ifcVs4MPyMv2DTbUAIBcyiCyTxlroaMNksaDcv2CnY5srzJiQadq9taRbU5SICh54RHViEwYkQk58Zlht0/SNjekfAOA4lRSUsiu0mwPBVczcSlH72rLHbKy/3Tb1Lo1ht/fi+Cd2pHgYtqdbsK44PKT1ChkKAxAkKepQUcbJI0HVas9dByW1ObLdDpJQ8VFONg6oY+k8GCpOpwYW79GSgbbtDtD6olkTFaXnaFAEAAUJqWUkc6bw8mqTCIqST1BDiUXrBWPuh38DzUkxQZpB69S1R4y57CsxvRpg6TxoFS1h7bUZCwi4+RYrg8f6CRF0kaVb/Kw3O6Qlc2nYzIqi0ePw6OnYxQstzsUCAKATDwqQ09crjsjAc/hkwTWQqN7dz4tpW3j7EhSVlJ5kOm4g/pYgrxWF+hog6TxoNJsj+rEZMSCRkt2krGIjIUstzt0FgZp+d4hOR9IE7/NX5jq66+FfhwmhdJDD47u5igMFasCQNKIyHBkyw17KOAJb0Wn0XE7vaGIcT4dkyHy9PyljMiEDLdPBtogaYapttxMIjp0Z0qKBVntYTuRNKRot0ZjJkDanmKP1SByxOFAkJwckl/CAN3KlmrDRgIAEtGILSEDV64P97koTMZlpPrKDWfISIC08qDRbJ/EHJKH27ds+qtqg6QZptpyPY7DcmJBo1llSSG7U3V7KDPRfwO4C5UbzpBzCdK2Oc+gE0iYhFuq2oXJ4Q9VzCUkpXYGm7IDQNKIyAhwjYaLQVoxgFdwVcq3Nz9SLgbnBfpSpC7oaIO0XCk3nOzvPvvk64vor0za+Qw9KKOShtyQo5E0GXju3clYRMI255FVlpQwKNWs0V8qE8evzDedzujeLelD+Z2H0Ldvu9MdPToMtk7HwrQ7bqc3FHKQ9O15ekjLqJ2dNkjLlf/nhVOm3fnCEyfQX9lTkCZjQqjnEVVSDWm5MRwIAoBCNi7Dyg5OzCNI6tTgdnuja+XTEj7UiNAA5IiJyQt6bKnRCfQkXLnubPLykNBn9HkaCfJ9ov9SpZo9pI8FmT3g0dEGably4O/O/s87rpw52UB3xgfnH/fJJKLo9arlhj2kaACpNaReYi18K9twlMn5hsY09EH/UDXL3TqyzSVjEfRUOZl2OloJkDQi6DqxUfUEyFFIe6onJBWbD04X65Mw8CMBktAGaVnidnuz56xrtkzuLqRfOWXivjgRNQw9KCMLErB3o1vZUZkTvLUp4GLandEgviTpoNvtebqY6At5nlFkBIIst+OZFStkE+hruZ2e51ro197oOBKCjLxspdn2VNlJGv2FjjZIy5KjZ5eKU4l8OlaYTMycrOO++NDcIELSkJJDGg3ig5wIg6cgLR2Pot+oC0see7eMFMhoP4j+WujGz3Obk1Hd4nZ6o2U0/feAu9Zow1OQc+1Z7nArKUIyFkE3fp650mU0gUIbpGXJK6cauwsZANi1cc3cOeTkhKdCujCJH0mz2sP9VAiZRBT3OFxtuclYZFQ9Ucjid2oY7VEkCT+JnRHBb5JWaTqj6gkZg1xHVZcEGZG0wam+fdCvPTifVvSwsjKmuowWZS8vlvFbX828etp835YsAGxfn5pdRDZIpZp3hAG9PGje9PBaCLg3qp8zIUVl5+P2odedlBv2qKIBADZ7rS68lkcGzohOyIjieqouZRgkzy47RmQCXX4y3/BQTwBALhnDTcL5XefyZsCjow3SsuRIubm7kAaAXRvXSFA6eajsZJy7/c74ZCwN4kKVpreiQUanhkrTySQ8pMPoxYmeeR2QkzAoVe1Rt4/sfbh2Yt7HQwLskF2l2c4koqPORCGbQFdqeArZCbj3lJ93XswldC87jURMu0NCKNvWpWYXW7gv7nn/FKeSMm7U0UAQACSNKO7u4ynTADnn7gA7gWtlPY0ESOgYW225nuoJArqd8KxfTsaQu0L076AhZBy8TNu7g3A+HcO9JDyPkuffgyNlKAk6Upq6B/Pyyy8/9NBDMzMzGzZsuP766++4445cLhfw96ZpPvTQQy+++OIbb7yRz+enp6c/8YlPbNu2TdX7HUeOV5aIYDqfjpGRz4jlpX6OC3rTSdNxM3HvKxDXTnjm5EFOmyI/O0FKkXbkh2Xu3Myb3pkJdMtHvj3PCwy981upZv/ydo/dYFMmjpsrrbbczEgbQAJ6FLfSbHveUOgNpUabKBJIyE5KqxVsVHtIBw4cuPXWW5966qnLL7/ctu0HHnhg3759lUrF7+/r9fq//Jf/8s/+7M9KpdKVV14ZjUb/+q//+pZbbnnxxRdVvu2xwnK7/WN4Ph1LGshd5pTlRf3sBH66xauuEySI30gjOxlqcs+1No+08wEJzW/KDdvP8KDXxpbrtp/jgrgKAFSaHvXLIKc8aLSRXR/cX8rvhpIR3pCEUoN05syZe+65Z/369U888cSDDz548ODB/fv3nz59+r777vN7yn/5L//lxIkTt9122w9/+MNvfvOb3/nOd77+9a9blvXFL35R5TsfK45XWjvyqf4tWpiMI2o6/a5pcsTDVan6HejQ0y2n6vaokB0AjChyF2RyUPDcPdH1xOW6R29QkGCQ/MpvAcCITOC6mKbT8bv8cK89PzmfjAYKfio7IzKBG4f0LNggyJipIQOlBunAgQOO4+zfv3/r1q3kkbvvvnvLli2PPfZYve5dTPPkk09Go9F/+2//bf+RD3/4w9PT03NzcwF+1crmeGVpR35N/39xxaP+eZ0IACBuqaRjnufejX8cNtueuw+6h1RueMfrQEKDS89mSCAhZGe5XU9nAiRMKqq2XM/XNCITyObcqy6NgBuHtNyu1e563lObp5D7q3oWbCwvlBqkw4cPA8DNN988+OCePXs6nc5LL73k+ZRKpVIsFjOZzOCD69evB4Bz585Je6djDfGQ+v9bnEoinogD8qK4N6qfKwYAm6cSuG2KKs22n+VDzip7decjoPfo87MT6H3zAn4pXDvhdnujc2kJ6Mmq+Yaz3udD4cYhAy4JdPysLPn5OvHM6D+NG0oN0szMTDqd3rBhw+CDV199NQC8/vrrfk/5/ve/P/iIaZqvvPLK2rVrt2/fLu+tjjOnavZgY7FCNo5qkLzDC4AtCvIroyEgB528xvkQcG1SgDOB3p02wE7gsjAyrbFPLmUgun3lhjM6iIuALkjzE0MCdhzSs7MqoZBFjkOadsfz2yMPdqPLwHlSp7LrdruO42zevHno8VQqBQClUin46adPn37xxRcPHz78/PPPNxqNP/mTP4lEQqzpzp07B//32LFj7O96HKm23CsvffuHy6dj83gHuooZdBxG1IkFnBzz6RiupMrt9vzcPhK189twWVFmJCy3a9odNRm4Ewut9xWznv/kGYnixu346kHQ61UDQna4cUjL7fhJSZMx5Em4laYz1HipvwdO3PyniAvJQ51B6nQ6ANDPHvUhdoX8awCvvPLK7/zO75D//pVf+RXiVwWzYizQEEOdFNanY6+Vm1gv7idRBQA/pRAffg2+QI7Kzm/rxNU1LDS9k1UAYEQmmg7eubvuBEixcb+9muVOgHdWD9e/DHBk0ZNVVrvrl+0zIpgdKIIDnlirEKotN5+5YK3+HvjOP3ihG1vj9aTxQrXsu9vldFGvvfbaH//4x3/5l3955513fvvb3/7oRz9qmshdrpcLpao1OLUBt6+B32RVdIJvVMQzvml3LLerJtlb9prXScCdr+p2e7jeSQDdHlxd9E4/bJ7C7NrpdntJw/uMokz2DQDJWATTIJke4y4JxRxyC/MAtw8AOjFFU1FEUGeQYrFYNBp97bXXhh5vt9sAkEiE7Bfr1q1bt27d+9///s9+9rP79+8/efLkt7/9bVnvdbwZEsIVspiyb7+qWMDWEy80237TOXGNR0BdJ2AHuJTZiQDtCWD/UuW6bUQVTfX1O6Og15AGnIdwNQgBOSTA7tkYUCCv7PgiiFIPacOGDY1GY+jBWq0GAFdcccXo3//85z9/8sknRyNv09PTADBq21YDbrc3dA4yohOII6sD+m4hb3OBNyoilustuiVk4gZiV5WAIypudMuvCSkB95eqNNsBxg+xb16l2R6KOPVBH8cV7EwgYrvdRNR7m82nMWMbfp1VCblkrGdoD+lCpqenHccZMjDPP/88AHhK5k6ePHnnnXc+8MADQ4+ThFM0iqwEXRaUG8PpkFwyVjHRDJLnkJg+iPdP0GkuGUM0sX6dVWUQ4EzgOrJ+lVUyCHYmcEvT1ISLg1tt4VanBSRl/XoX8VGuOwGOrBGd0LLvYUgF0qOPPtp/pF6vP/3004VC4aqrriKPVKvVmZmZV199FQCmp6fT6fSPfvSjoXTRE088AQDXX3+9urc+NswutobOQfkMctDJr29QJhHFjA36a9vQXTH0DIQflWbbT/qBG/UKlvMlDbTJb+WGE9AMCXd+h+12A34p3A8V4B4lDMzxvqbj3VkVsAU11ZavkB3kTKeVgVKDtHfv3ssuu+zhhx9+5JFH2u323NzcZz7zmVarddddd/X/ZmZm5rbbbvvkJz8JAJFI5F//63997ty5X/u1X3v66affeOON55577vOf//x3v/vdd7/73Xv37lX55seEUnW4fCefjlluF2tfCIjj+3Ul4V0oKDOBeKMGl6/j9sZWJvs+Vbf9MnCAOsu10mwrG/gW4EwA6i8VkCglIMYhA8obilMJxI4+5YavTAMkqEIkoTTTFYlE7r///ttvv/3ee++99957yYN33nnnvn37/J7y2c9+dnFx8b/9t//2qU99qv/g+9//fpo6pBVJqTbcWIzsF6WqvW2daIyYyIUDquIRhcsBmYlkDLNdrGJnwldfjppDqpjtweJoeVSavgcUwNYaKMvrBEjsACCTiCLWUZRqF2hi5RFQ0g4AmURUF8Z6cMkllzz++OPPPvvs4uKiYRjT09NDlUk33XTTUJLpnnvuuf3223/6058uLS3FYrFdu3YNVbyuKky7M3q4SxoRE8NUBB9RcdXYwZkJxL173vSVYhOw1gp2JgrZBGJppzJXzO30Ai4JXK1BcPMOQP2lAr493C821Mr6tVdgZcGnPxZh/ZrYssghXQQtYCQSufHGG5mesmXLli1btkh6P8uLUzV7+0ivhMJkvNJ0ANKCLx6c7C3mEicWcIYB0oxpwBryVDHb29f7HlERE2PBzgRuzCRgYh4AJGNoWRDP4eWSCLYTpIkwSqMQGiMhvgq81Vk1wN6QqgMUgxR8mlwurMao17JmqE0DwYjgZEeDq1sAb+KyX7fTQfDsRNA2h5gYC3YmAHWoQfAvhagTCw4E4RLchzSXRGtnt9Bs+3VWBdRIQEDHPAJiHEYJogAAIABJREFU01i32wuIBCRjkV5kGZQiaYO0zLDa3VFxGlZ/4lBBGta5u9IMOaIiboLBzgTgJcaCnYlcykAMeCqzE6bdCfj2MnEDS6Afmr8EvLl5wWcUxDyW2+mFqiuxPlSp5j3bkFDIxt3EJMpCUtEGaZlRqlnF3HA2G6vCrmKG5JCwjqhuN8SZyCSiKFkxCMtMIB6Hg40Ebqm8ZyqxjxGZwPLGmk4nwJlAzCGFZsWMKNqHCo4GI1aaB8yKJSSNCFZddughb1mgDdIyo9xwRjPnWM35gzeFYg4tLV+q2rlUyAG/0lQRdEIMuwc7EwSU7ZukwYMEFJNooy7UqSe6Ic4EotQluKkB4E0/Cf32CtkEVoVDteUWvEbaE4zIRC+qKPQqgjZIy4lqy/UcOoDVx37edHCLjfwIja0jEupMYO0+wc4E4lqh21x0YgKrkibYmUDsBFhu2MHOhBGZwEorlushTjOikQj1jLG+QNP2nv5OKOYSncQUykJS0QZpOUHUQaMbxKYMUsguUGvg14mZA7/BoH2wmjWQ1Fpwaaey4zBWsx+a3hNoe3egIA1x1EVoVgxRpmg6nQA7kU+j+ZcBQy7e/hukMjhlvqxUtEFaTpQbtuf1nYxFbIzLOqRUHq9xWafXSwc28solcZITnhqQQRATYzQFJSg7XWh3vtBNkJ5y3Q4IBAHeF2jaQUYCAPLpGJrbF6g1QKyVpkmMoXx7lttNGpHl0o4hAG2QlhN+1zeWzNftBAnSjOgEWlODwHY+gCeHLdWs4MxEPh3DSlbRHFFRMtihAU/E2dim0wnuAYoVtQs9fySMiBq3D9Gcz5tOdCIkMYbyoUKbIeWSsU5ctE5RAdogLSdMuxNcyi5IqWYFK0exFgqNrWMFH0pVj7KtQZIxtONwqKQKr1zMCTXYWL5sqNuHVYQ7H1asinX8N+1OMhbkTCBe56FHB6wPFXrt5VJG11gGrYO0QVpO+J3ssGpcgiffkJtHgZwPEcvtBqe+EKMcpaod/GpYzkSoAUBshxoedEIK2ZUbzpSSM4rKXEvAhCcC1uUXmlZccjquFjVocFloek/ByaUMlNxyqNYgiRQ2KdWsYDuBKGoIDmXk03EUr4UMSg+O9iSNKIo3dqpmBwc8C5Nxy0W4Hmi6nWJ1Fg9dKxmLoFzkoc4EIF5+9eFWyEMUJuM42V+zHRyabra7ulODBhmpE8xMu5M0IsHxGayqndBT6lQSZ0dYaLaDY4NJI4IyAoC4PsEfSrwdO8Ht9oLV+UkjipJTpFEtY3lINPITFP+SpskFlkEKbZEVnZhA+aVCY4NbVI1zFEQbpOWEX6kByilVZQlkteUipo4DoJE5Ya2VT8eCdx8jMoGSbgndUrGyIMFjGghYU+ZCxzSgfShTacguWA+JFbILbVVMvj3cMfAy0AZpOeEny0YRqlpuUFtixIXIkTC06SRWi7lgnxJNPRGoByFsnkogHYcVNYmhC9nhKDyJ1iDgD7CODjRdsTNxnM5V4XYCKeQwbzoBAxv7aIOkwSQ09SJCaKtvQHLF3E4vlzLU5P+DO07CW/l/8eRE6FkYkUqzHVIbFMVxxSpmSMAT8EJ2NII0lIXcbi9Yig3na2NFrSyx08FfIFbXwYrpnV0eZKKD1ttXHtogLSf8YjUoBfOlakiqHItKWF4HdS0qZ0Lc+IXqQRCx2t3gL7CI5oqpG7ETOv4Ka8LhvKmoUXql2Q7OigFeYowm3m44dcS5l5LQBmnZYLldz0Z2gORP0DSkAQDxhEGl6YTKnLAoN5y1YTcqyvYUMLy8TzIWQYlDKusESONMoNQh0cQGsdItynJIlaYTrPkGvG5SNJcf4PXNk4c2SMsGUuroV2IiHtCoWeHxZSOKcKCjOc1hyZzcTi+40QAgdYUILZUHgEIWp0ka1izdUObNcDuxKYOQQ1LrNIdffoVJhC5ZbqeXS6qabUgRb4duB6UYQCraIC0bSjUrYHco5hKCFULB8zrPrzKVEBc10JzmsAxSteUGzC8ffEuCC9E4EyjQOBNABu2IJ8YockgoUMr5UBJjoXMuAE+8E3pGQczAhVZDx6ya+EKy0QZp2RCq2BF//eDpCVhQbqniN2ql2aZxJlACniozEzRGAmUkkrpKgE74z4SVGKNyJpAuv1AjUZxKoiTGqi13BbT6Bm2QlhHletApUnyKOeXuI77NLTTbibAb1YhMiMucKO/STNwQD2Uo2xFonAksLLeTSYQnxsR7DdA4E1jQ2AmUAWPKLgly44dHnnsdFBdTKtogLRsqPn2DCOK7Ko1ODKV9C1UQP4sQGyw3Qhq39MFQaoR/KBRXjMaZwEJZHRKNkcCCxk6gDBibN53gGSsElGsvVM4HAEbrHIqLKRVtkJYNwdoqcbFTuRE+VhylcRFNs2qUPbfSbNM4EyhKjeAROwQUdb5KZ0JZJQClMyGeGKN1JgDEZy/R1AYlYxEMc+7QVA1HHVOZcoQbbZCWDbbbDcjxiB9UKXNUC0viJTtUBzqU2iAaw5ZPx1C8MTWFsZR7t3gDeLfbUynno3EmxBNjlM4EYFx+oY3sAKCQRXDFynWq/OWGv//vv7J7g+BastEGSR3PnDj33JtV7qeHBoUEr2walR0gNTUIPdAZUYQKdkpnAiVYRKM1QEmMVZrt4M6qhFwyJjgMsNxwKPducWicCZyF6JwJpOIq6oix2J1rud3Q2MZyQRskRcycbNz4f//khj9/+ejZJb5XKNXsgK1TcL645XZpetmhUG6Ep+VRCj9DOxr0EZeY0+WQ4uId0mhqg1CYXWzRGAmUxBiNM4ECpTOBc/m5tJef4Bdo2h1lXUJkow2SIr7xfOnLe9/5f+15x5/86CTfKwR3wUkYEZFzVqXZprFGKOVBpt2hcSZA2E4sLFHpBo3IhOBCpKVT6JaK4orRjIQgiIfsaBonIiXGaJ0JQeidCXEra9odmihu0hD1xlSOHJSNNkgqsNzut4/8/FffV7jt6sLBowt8LxLaUlMEypOjuEEimxflgU58LRoVhrhGgN5IqNEoAkA+HRPN/9NdEihQOhPiiTF6Z0LZLyV++VG2+l4WaIOkgmeOn8vEozvyqd2XpMt1h++mKlXtgC44uZQhIjeotpQ1+GqHtvrGIrTVNyE6MSHYYq7SdGjOwihyvnLdpjmXiEdfKUO4mbhRtRQ5E+KJMWXOBE2r7z6CV4XV7qrJwClAGyQVfO/Y4od2rQcAIzJx7dapmZN11leottzgjjuCU8zLDScTpzqici9BoO9alkmIzqRxOz3K47DgjkAzhBSQuriaTie0O9/5vxQ744cOKCLk0zE1Gbj+uxJZiNKZyKVErWy15dJLQgStLE0MfLmgDZIKjp5d+uXtOfLf126dPDTHY5AAwx74QZ9DEtSJmQ7tmIZ8Oia4lrKmBpQ7AskhCUrMKe0ECFvZMXQmxF1MSmdCPDRdbtiIw4iD0TkkDRszJ+s3vXsd+e/3bcm+etpkfQXT6QQfuARrBunbnAh6LabdUSZRNR06OxET7aRptbv0QUjBnjT0OSSRVQBg3nRoMnDiCmmr3aV0JgrZuOBaypwJyiIKQDJ+Ik8fK7RBks7Rs0uDqVS+GQSVZkhJjWDNYKdHFd0Sn1dLf5oTzzNZ7S7NKVVQMQ8KB9mRb4/mm0nHo4It5qw2VQ5JvCKb8twAAEZkQk3DRoRrz+1S3iy5pGjMs9pylQ0Yk402SNI5NFe7Zstk/3935NccPdtkfZFSNbyvjwjlujNF00NauNyBvuKkkE0IJgyUDbKbNx2aYlUQHjdg2h16tYKonbA7avrL0UdWxe0EpTNRnEoKGgmVYTSV06Rkow2SdE4stKY3Z/r/W8jGLbfLerm7Xdr8PB/0949gEL9m0eaQBKc7k6yYGjkffSComBOaw206iowEMPqyIlaW3j0V3+LpnQnRGrhmO3RcLMGICn17ZMaKziFpaDk6v/SL26YGHylk46z9GkJ3h0zcENHq0LT6BoyOPpSuGFlLZCF6IyGunrBchhySCEzOhGgrKeq0fCEbFxkOWW25lAHPhCE66oLSmRDXI5QbVBk4EI4YUzZsXC5ogySdV043pjdnBx8pZBPHK2wGybQ7wXJVwQ4uNFNcAUC8uRn9uVtwLXojIa4vp48NGpEJkXED9AuJz8eqtlzKdrGC2zdTxEnkQ9E7E+LXOf1pQNS/bNgrpm8QaIOkgNEteNfGNaz3Fc1AUrE6JHXNqpU5E2o+EQCU6zal4yJ4bqBvCZGMRQQH7TDZCZEtlX7YrmBgiikDJ0i54agpVnW7PcrY4LJAGyS5HK+0MvHo0I20dW1y7hxbIiFU8iRYokGvNVB5oBP5RJTFqiDstQCLTixpREUk5uU6wzAkZZkJQStr2p3oBK0vy70KsFwS/b/nXsvtUIUcQHjkbqXZziW1QdLQcfRsc/clmaEHt61Lst7AofcSCTJwO0n0dUiCsSDL7VLWdQr2nmAQGkwlSzWhSg5lg+wWlqhmT8D56gL+n4l887QxT7FpxbIFO30sl6Erdi5liNQ8VZq0xm9TRkg3T5+BWxZogySXV06ZuwvpoQeLuQRrEpiyKITvoFpuOEZkgvL4KRgLKtdpI2mC/SJLVUWxddI6gXLvFiwjZQrZibhipZqlLDNRbjiUQSdBV6xcd7atS1H+saCQutpy1XQJUdkDVwHaIMnl1dPm+4rZoQdzydjsYovpdcqN8CwFd0CDUmLXRyQWZLldNcLlmuWqkfMxVZyIl5GqgVLkgrNWhzY2KDjqgunaEw5NO/TVvtyrAIDb7VE6zcsCbZDk4lnDX8wlSjWb6XKnm/8W47uFKMdoDr4ZjlXgrdoOyuhWJhEVacJNbycEXQGVJZClmk0pAKNplRsA0wHfiApKBxW1HKT3LwFDtkp5XQkGexeWVk5VLGiDJJujZ5s78muGHiT7F9NJmWbXMyITfEWXrMle7oQBaYFMeSQUjc80aJuJkffDffQ2HZd+9xc9Dnd6lEeHfDomWN1Cb6eLUwmRJByTnRCh2nLp48BJI8p9nbPOWBEpwl1Jrb5BGySpVJpty+163thMtYSUA0m5z+mm3RE8UFOibOsB6qwbQaQrhNWmlWmA8DS2Us2i71omJOdrOCoV0tTFVUJKjZql6PJjMhJJIyJikNxuTxfGaqg4era5I5/yvFyKuYRp016FlEEh7s6nTHYiE+eP41etNlNwRky47NDLYUWidkwhO8FhgAxun1hiTNnMt2rLpe/wlIkLlTArOw+Vaha9kRA8o9BfEssCbZAkcuRMc/rSYUUDoZBliHJQBtm5FVwLSwxbaj4t1FacnqQREWxIQ18wKLJ9M+0Igg366Es7BXsN1CyW6JaAdFCl00zUpJR/LDIVgr7JhTiVpqNsnosCtEGSyGvzzZ0bhxNIBCa1FeVsFW4FF1O5pQilqk2/kOKTI7c3ZtqdNN0IVxDrNaCy0QB9XRqIVdKo1BxWWwzd30UMksr220y/1PijDVI43z7y85vvf/UrP5hlDcoHXCtb1ybpQzeWy5AO4cDt9igbQYLa47BIyI4pti7iuNjUpb4AQP+Xo3DI+bg/lDI7UW216bNigrLvSrOtxplg+qUE2yILNhkZN7RBCuGhvz3za//vz3bkU9957eef+e+vMz33eGVpR967EK8wySBqkH3gYio6EamkYd1SuTPYrO23C9mEiKRKzRGV1bqIDBEeHCkpFaZrW1x4qeaXqlkMBkb8Q6nJ9qlh5egFZVBuOHd/543Hbr9iz/a1lWb7nX/wwu0fuPTarZPhzwQAgCPl5q6Nw20aCEkjQr+tUwb0+UN2CvOilF3LAMCITHDrxFhH84nolJi2ORFXjCniBGIxz3JD7kDIgYXUNRpgMs8ikYByXZGRYOoSsizQHlIQn3r02KevL+7ZvhYA8unYb9+49T8+PUf53FLNThoRv5uNqdUY5fWdMDhvoWqrLT4AhoZ5k8HyiaSFKybteLQ+3AkDJo1vPh3j9lrKDXXpa7agk0AJs+B4LXqIhojeTogkxtxujz6tKILKomw1aIPky8zJxnNvVv/dL7+j/8inry8emqtRnnCPnDH93CMAyKcZpnLJLjWoNBnU2CK7D9H4Uv6xEeEfBli12FogG9EJboNUqln0tlNxDolbJM30S4kEnU7VFLUcdDu9pKHImWD6pYRcMYX+pRq0QfLlPz4997kbtgwefnMp46Z3r3vkcJnm6UfKzd2X+BqkXMqgl31TXnbcfezLdSdDfaBbv4Z/92GNpPHn5Bl1gyJ3NatSQ6TPBdN+mksZImWkiuwEozPBrWtQqUZjysCJJGVN211Jmm/QBsmPUs1+5sS5T19fHHr8f718/feOLdK8Qs1yA6RrZAuj3HCrLSoPhu/Kdrs9y1VUBUk/yE4QViOR5I12Arvjwr370A+yI3CPyaHvmNeHO9vHOj2B2xurWh5dJSVRqllqZN8rrCoWtEHy4z/9fyc/fX1x9Kras2Ptc29Wac5oxystP4kdAORSBr0ISmqkmOyP9EdvEQE6/SA7eKtYlbddLPM3xmcnyEIM5ZYCs9Ssdpd1JAefJpi+Yx4hk4hy981T5riU62wZOJG5eW63p6YwVueQVgWW233wxdMfm940+k+FbPzarVOP/6wS+iIvlxqXTAZdlMVcgnK3lXoOqphsmnKRhAFbYWw2Drx2Yt501PTkZy1WzaUM7sIR0+6wnvHVqC5F9sRyw2b2ZbnshF9XST9EImnj2S52WaANkgcPHDo9fWl2l0+ThX9+2drnZ2vBr+B2e/94zvrAO4IE4sWpJM1UJMvtup0eza7HN8GlajEXOfEF8UltENNAGuA947O2QDYiE50ez0KsBU8gELJjLUdLxiJ8/mW5YbMe8LlLmC23y7RWIZsoN3gaSo1tnwsQsLIrD22QPPjG86Xf/mfv8PvXa4rZ596sBr/C8UqrmEsE3wCFSaqmcPReObfvwnRy5B6vpzK8wKpL5P7qKk3aAbh9uPdu+qnYBG6vmuOAz2ckAMBqd9U4E6yXn+CAPiYKkwya20FO1eyVVBUL2iCNcmiu7nZ7N717nd8f7NmxtlS1ZxeDJg8dKQdpvglJI0KjyqUXjPHpR0tVO5Ng2BG4M8OsmQkQGDnIGnTiLuFSWtfpdFiNH5+Vtdpsbl/SiApYWUY7wTsMcGGpTV+UDQIhu5IqITuobU2rBm2QhvnqM//47375HQE3pBGZuP6duWAnqWKG65sLWapjEdP8Nw5Y55cDAF/rLY69O5OIcjsuauSwrLFB8sd8NU+UkdtB+IKrnmOOA+DefyvNNqsvy20nWK9zbsvHYSS4vTEtaljhELX3x6/ykDMMsvuS9KunzYA/OLHQ2u4vsSMkjMjCUvhuS39955KxqsW8fXOEMviMBF/NBF/VDmsQn/uYybqPiBgkZRpfVmcCeD+RaXdY91PuvVvZyEGOieyC7exWEtogXcB/fHruEx+4NHR7+uDOdcFCu0Nz9d2FkJAd5VXIYJB4G+YzCXW4C4k49lPufBUruZTBtyNQTgYRp9pyWZ0JbuEyhxSb28QqHNPAVpTNd7wjsOpcuMelc8hPxhxtkN6m3HAe+tszt3/g0tC/vHbrVLnhBLRaOHq2GXqh5NMxmobfparcvCVTfzkQMBIcDZC4bzaOCU980a2a5SYYvxC+czpHIIi7G5syZ6LaarP+xJlElM/4VZpsDRu5j3cqw2gqpYNq0Abpbf7gydlffV8hoJq1jxGZmL40e+SMd9TueKVVabb9VON9CpNxmiOY2+1Rxk/4uqpYbYZxPiLwDdHgcFwst2u5XSbDyX1XV5pt1ufm0zGOnU6xM8HUqYFs9BzC5UqzzXpGWb+G59sDxoaNInAcHbiNn+7UsGIpN5xHDpcHW6kGc82W7KG5uuc/lWrWro1rQvcpIzJRMcN3W/r54nyRaNPusO7dfM6EaXdYi/iMKE9/VT4joTKrzBELMm1m7Ql3I1FWZ+J8CTO7cFmllTXtDn3DRhDoecrhnfMZpJU3ewK0Qerzxf/5D5TuEeGDO9d99++900hH55eu2RI+M6k4laTpr8qxvTLBuqVyJ2BZY4PAK9+qmDxGgq8PKUe6hS/mWW44SYPtMihMxrlVdkzCTu49kVWjSOAoYSanDVbpIF/As2a5evYEN9ogAQAcPbv0rVfmf++mbfRPuXbr1NGzS543/OFS48pLM6GvQG7F0IO57HGxlaZTCGxxhAXTRAMC305XtZhv1KQR5SuVLzds1g/Flxjj+/Y4+5CyW1m+X6rT63GcUfhcsaQRUaOR4bATfFGHlRevA22QCH/w5OwnPnAp06+bSUSLU4mjZ5dG/2mmVN+2Lhn6CrmUQbNllKoW/fXN4b6YDlsogxvW2RPclOsOa+WWSCUNR5UYh53g2+ZYVwEuZwJ429mV6w6rM8HXD4nPaeaTmHPYCb6RLpUms9M8/lyEKt+XX375oYcempmZ2bBhw/XXX3/HHXfkcrngpxw4cOCHP/zhT37yE8MwpqenP/rRj95www1Y7+d4pfX4zypvfOGfsj5xz461B48uXLMlO/hgteUeOdO8dusUzSsUc4nQer1Kk0GJRI5aTPceRwMFEvJm7i7DPnuCM7budNQMwAU+Z4IrMQbsmz6nkWg4TC0HCckYVduRISrN9ge2hge3xeGbPUGOd8xl46qUbzpkh8CBAwduvfXWp5566vLLL7dt+4EHHti3b1+lElTT8+lPf/p3f/d3f/jDH77nPe/ZsmXLD37wg0984hN//Md/jPWWHnzx9Eev2Mjx035w57rR2Uilml3MJShFSpl4NLQJGKsSifWoxXGg47YTrDbs2NmlZ06cY12IozaIryyf2BXW3Yd11BBh3nQ4+jpzVLdwdO6A8+FB5kgaX2yQT37C4cjy6Ro4aoP4evta7e4K03yDYoN05syZe+65Z/369U888cSDDz548ODB/fv3nz59+r777vN7yuOPP/6DH/zgPe95z49+9KP777//kUceefzxxzdu3Hj//fe/9NJL4m/JcrsPvHiapvZolJsuWzdzsj6Uezg0V9uzfS3lKxQmE8H7IEmrMt20TH19qi1XmVCHwxWbTBocQQm322Pd5opTCfoBvn2I7VdzSmUdZAcASSPKkZavNHmUb9wVY6xr8WlqSlVbmdPsdnus/mUxl6ApSRyCu238OKPUIB04cMBxnP3792/dupU8cvfdd2/ZsuWxxx6r170l1N///vcB4Itf/OK6dee7nb7rXe/6rd/6LQD4m7/5G/G39OTri9vWJa/lihtkEtFt65KvnLqgGunY2aWdYRVIfUKVPFa7mzQi9Ocg1jtB5VxnDlds18Y1oeVcoywstdUMiSnVLA53h++Mz6Ft4U6McdTr8J1pOJwJPq+Fz+3jg0N+wsfKa/UNig3S4cOHAeDmm28efHDPnj2dTsfP3fnxj38cj8evvvrqwQd/4Rd+AQAajYb4W/qrl+c/csVG7qf/yu4NB/7u7OAjT76xSG/eQvcm1jBxIZtgyk/wNavOxKMcCQOOTg25lFGzeOoz1Ow+bpfZ5wPeWQNuh9nt46NcV9i/nLEGjhuOClxulB3y6EvmlxFKDdLMzEw6nd6wYcPgg8TYvP76655P+drXvvbNb35z6MGf/exnABAqhQjFtDsHjy786vsK3K/wkSs2PjpgkGYXrVLVplQ0AEX8odpqsx6CmAIafA1P8+k4a8KAT6LKl6wq1Zil2AQOx6U4FS6nHCI6wechOSLjz+kxHZ6cPLdwWVlhLMflxzE3T+U4EgBYeR6SOpVdt9t1HGfz5s1Dj6dSKQAolUqezxpV0y0uLn7ta18DgI997GPBK+7cuXPwf48dOzb0B6+cbhSycb4kM2F6c6baco+eXSKRpUNztWu3TtEf+tanYwuB9qPcYOsczLqVKCtlUNagCADcDnN1CwAkjUi54TBdDCqn0XBIQvLpGE/JTt3hcCZ46g3sjtvtMXfZScY41BPzprMpw3xJkOm0TEFjPqEB38HLcruhBRtDe+D4o84gdTodAOhnj/pEIpH+v4Zy8ODBP/zDPyyXy3fddRcJ3AUwaoGGeOZ4dc8OWgGCJ0Zk4uNXbXr01bO//79sA4Bv/vj0be8LGV0xSCEbP+ZVydSHNT/P2gJHWcykVLM4nAk+yg22kYOEwmScVWjHHZzhkFSV6w5reDCT4Kn2rVnM6gk+SjWb49rLpQyOcPEbP29tyam4/DhG2gOvQaKJ5w/ugcvCOKmWfXe7nKPjX3vttY9//OOf+9znbNv+6le/+qlPfUr8zTz1xrn3FbPhfxfIJ6/d/I3nS6bdOVJuHimbH59mMEihOSQOD4bpoFpuMPfdAi6RNJ+HxNku1u0qqsDlkjlxNFAg7WLVfKhKs71+jZJkVcNWY/kAIBWPTCUVDUNSNg+iXFcnHVSGOg8pFotFo9HXXntt6PF2uw0AiUTQr3jfffc98MAD8Xj8jjvu+M3f/M3JSYRiOsvtHpqr/dn//m7B15nenJnenL3hG4erLff3btrG5LCHBjoWmm2mUvZMIsqqQeLIi3LEZ/hig9waX9anECrNNk2LjT5Nh0c9UZiMn1hoMT2F9DNUk5bna1XF8Uu5HWaRCzdWu7tzA7Nck8Nxkd3oaxDT6ajpX64SpZ0aNmzYsLg4XElaq9UA4IorrvB71r/5N//mqaee+uAHP/iFL3zhkksuwXozpaqdjEU4VMWj/OX/+Z4vHvyHfDq2/7rhDFkwoeqASrO9dS3DFsl6tp03nSsvCW+7NwRHste0O2p2n3LDcbs9HrcvMsE6mr1cd0LnAnsuxPoUvs43BNZEV7XlcrzDdJzZly031DkT1VabQ7zDYZD4orh8A/rKdXWSEGUoDdlNT087jjOU2nn++ecBYPv27Z5P+Yu/+Iunnnrq13/91//0T/8U0RoBwJGyef07cyi7ZD4d+88f2Xnvh97FGhMPVQCXaravAtreAAAgAElEQVRU0c6YD0PigDiIHMkJjs2RO4fEvM3xji7lmL1UbtjKRPPK8v/KuuzwRQI4OourjOKqRKlBIhVIjz76aP+Rer3+9NNPFwqFq666ijxSrVZnZmZeffVV8r9/9Vd/lUwmSSUsLi/O1ac3MzsHuIQGOipNNhVpJhFtsqR8ld2ofJ1vOGpIuT+REZlgjXbyVbcYkQlmI8E7KJ1jx+crrgIApgsPeM05n0FSVlzVdDrKZk+AwOCPsUWpx7d3797LLrvs4Ycf3rp168c+9rHTp09/6UtfarVaX/7yl/t/MzMzc+edd65du/bQoUPVanVubm7t2rW/+Zu/Ofpq11133R133MH9ZmZKjVvem+d+OgrkIB8w3rTacpl2h3w69lq5Sf/3fI0g+TS+792UZl2I6G7ZFnL46+RZT6msxwUCx1Mst8uhG+w/l+nv+Sar5lLG3DmL6SnKuie43R6fN8bRFaJcd9S0i62YcsekXSyUGqRIJHL//ffffvvt995777333ksevPPOO/ft2+f59z/96U8B4Ny5cy+88MLov27cyN9hAQBeOdW490PvEnkFFEjUzi+XXqraUtXSVavNUW7JlzDgiA1yeEiKa4M4tznGhIHV5g/OMFXSuN0e39A8jqfw1QZxwF1stynDHEmrtlw1w1w4hn4tC1TnxC655JLHH3/82WefXVxcJLMkhiqTbrrppn6S6YYbbgitJeKm0mzvLjCf2dEh3oafQXK7zDWeTKZC2fatTNRQbbl8HQ2MKLPxczs9vrQ8uyvGmYFjHhHScIC3CR7rt8dtZcktQ78dK5sHAbwRY74Gfdog4RCJRG688Ub16w5y9OzSjnxqHHzeTNzwE3eVarYRmcDVkQ9RrvPUIWUS0XnGLZUvVc43FYKvuoVjFGm54ai5hBaWOGuDcskY0xdIOubxHR1Y1fZ8YTRgH/pVabK1O+ljRCbsDrOd4GlvyC5qKNcdDt3g+LNKJ8bOLrZ25BEE3+IEmJCKyTwcjOkOJ7F1DkEaT9FJt8fVNI95oVLV5ouZsIYHSecbvtFBylpCAKM3xtHLh8BhWpS1reLo6kvgmApRbjBPoeTD5KqBG39WqUE6Um6Oyc8Z0IyAb8wlvURY2Y4A3KEM9oKnptNRU//PXayaT8c5qn25XTEmK8sxD4LA2rMKuOZjEViPDionq/Jl4ID9Q5XrirxzxaxSgzQ+o0QCnADTZq7EZpIUcyeQ+ApjldkJvmFIuZTB2mKOb5vj6yE95s4EALD2POXuX16cSpZqDIo+7gFFHJoa7l+KNWrHMYVyWbBKDdKRcpOpA4I8ppK+I38qTeaQHZMkjzu2zjrRp9xwjMgEXxdX1roT7ht1KmmwtgFUOjeIq345l2L7UCIiF+YzCnv/cj64q7j4QtN8Fp01XVqzXDUtBxWzSg3S7GJr23gYpIDNvVS1edrMUN9CykIZIqplVoMk8qGYlE6m7WbiytItNscodyC9scdSdQm8ghoO3G4vwXUYYhW/8fUv50PlL6WSVWqQyg1nTEQNAdNX502eO5Z+A+KWw7KGMiyXZ0gMH9wfimOUFF8KhNhLpqP3mHe+AYBMnLmBQkA9eDAczoSytlXcJpZZH8vVaX78WY0Gqdxw3E5vTHJIAR5SpdnePCnxTXKHF1gbKHDHBoE941Ju8LiVwJ4QEmw0wOq4cOarGAN9VrvL50ywNs3jmxvUX0tNHDKXjDF1EWTtKjJIQPmHJ6atVXYrBVKIOiZtoPLpuN91XGm2WTcUsm1Rnp64YyYcMifub5vV+FVbLncPaVavhfvczeQZkF+Tb61NmbhfhtKTTk/RSAiO/Ggf1jMKt3/JGi6eXbSUGYlyw9YquxWCyiFaoQQc90pVi/VGItcoZeC7Zrl8gjRWRMpoWCFqbI4nGpEJJq+Fu1gVGEt3OUYUDsK0pSrL61RMoe7vTEcHZZ0a3G6PL9UH7LMoRQ5e48xqNEjKWv/SUMjGqy3X0+Gw3K7Ua46/MTZjEF+kowHTWuSb5PtQrCFcQTtBv6WKzINgdXf4o1spA1g+VNXiaeHah+nL544Ys4oaRFJ9rHFIkajDOLMaDdLsOYtpMKhs8umYZ5CtVOWZ8cyq9OXA7w37wTdZlcBUn0G8AWWTAEXsBH3CgLvFDrB3URKJbgGLnRBxxVjj2NxWlrU2SJkrRm5wrbJbISibWUJJMhYZdQLIYZ9jy8ulaLOj3Ftq0oiwVpXzpcqBMV8lUhvELB1s80sHmRxf0+5wO8qss5dE5CdMmE6Huyv2pkycafaSsrJiETkfkzdGDL82SCuEct0ZE4kdwdPh+MeqvXaNIfWwz11uyUq15apZyHL5j6ih4+RH11ITSRMJBLEaJBE7wXZ0EEtW0ccAyEWurNM8d1KWadQF31F1WbAaDdLsOas4NUYGqZBNjG6Fjtvl6yWRNGiPWtwKaWDcfQS3VPqF3A7/jcqh3eIemgcA9IkxETkf69cuYieKuQR9Rx8RQY0yc95/Bfo/VmMnVqrEDlanQSpVrW3rUhf7XbyNZz6z3LD5On0VsgnKo5aIUIc1vM69KTC1KSo3HBEjwYSyDLZIxQmTY8rd/Z0DkW9PjQ6wD0NiTNU4EpEo7pizGg3S7OJ4eUieik8FrdJEjt5M4rdyw+ZeKDrB4CEpm4oNYlr2TCLK5I1xZ+CYUNn9XQQmgb7gh2J6rsjRwYhM0Pf2Fblzx5xVZ5CqLdeIsk29k836NbF5c9gJ4A6eUB6KyQldjfxd2YFO5EZl7RgLqrLKIp3pc8kYfWKMex4EgWnOk6BGkcW/dEWc5mSMQb9TbtjcazFdfiu11TesToM0br9lJhEdnQM2d87iq7ssZOM0t5BgGQ0TQnVILCdHAOAuVmXKIZFhvnwLsSIyEiKXMvw6JY4i2K+TKQ5ZbbW5B54y3b/ibh99YkzETrCqSVdkERKsQoM0ht6uZz6mVLM5q1YjEzRjLgVvVPrtm3uyKqEwGaef2nmK90tjRTA2yFTIpSySVrXafGlLQtJgiEOK5JBYE2PK9u5qy1Uj311o8ozuXBasOoPEPRNTHsWcR7s27s7BlM8SUUgDS4s5cmoWkfPR/7HpKKtM5Gz1TUjHo/RFJ1a7KyKaZwijiUmxmQyMMjth2h2RMwrTm6y2XGWSkOiE9pBWBCJdHSWRNKKj4eNSzWKatteH8kq12qJ9iSh3OpWp8nJdKJRB7/Zxz3zjwHI73B3SkrEIvSvGPQ+CwHQUsNpd7kAFUxPuSrMtUgVPfw8KzoNg6mW3XOQnHKw6g8Q9z1gexZyHULtU5U9l09Sxi4cuKe9Ay1VUfgsApiMUSaNPLFtul9tIsDJ7zsomuVtCUCUUCSI9igj0W6rI5cfUhHthSbSNFmVizLQ7ItpdpgycoNM8zqzMTxWASJNmSRSnEkMV9Zbb5c67UF7ZggWDRpQ2B2va/MX/wHhyFCxWpafcEIpuJWORUV2lH067m4wK3aeUv1S57ojcGvRbqkphkWl3RPI6yViE8vIrN2yRKC4T5YaiXKl6Vp1BGre+QYRCNl6qvZ2SIREhvugT5SFXcFMoTiUob9RKU6ivM9PJUeVQdpHoViEbp88h1e3OZrErllIVUrOELgn6yJh4Kymm6JZIFJf+UKhyYp7KEcyKWXUGqVTjb+Yvj2Lugu5BlWY7n+HfGih3cMEOs5SrCDoTQH3ABwDT7ojsPgHj5IeYN0WD+JRbKslpiRg/+sgY9yipPpQ9T8sNhy8/SlDW54KJSrMtEsUNGNTpuZaIHnKcWXUGSVk/YyaKU8nB7MXsuda2tZzNjZIGVZBBcEulP+EKOhNAXQhCclpiOjHa/qpWu6vGnItrcHIpg9KiC65FX69q2kIeEpPVLDdsNd6YoDnPxBlE86bdGbfaFSxWoUEaO5UdAORSxmDITqREsTBJWxgrcqNuysQp61Vtl19PBSy6WxIEU6O7dbs9EctHv3MJTlYlUFp0wUuCPtxXbjjivbso7US15Yo5LmOXGCN3t65DWiGIqNfksXVtcu7c27vGiUqLr9U3gcb3LzccUZUdnSCt0myLDEovTiUptx7TEZI5ActMmnLDEVFq5NO0wmXByapMlBu2SCSNPv8vGFkFFm9MMK3I5AcLGiRK6SDJiqk5eKlnZX6qAFRWbtNTzCWOV1r9/z1eaYnMtKXxkKotoaM3fSBIfHYL9dbjCAZn6FuYV1ui6gna3acu2mOXfhyJuJ1giEOKHfCZrijBL5AyMTZvOiJR3Hw6RnlJqFRPqGd1GaSxLSjLp2OziwMGaWFpR54zh0TZ47JUFerGxpQqV5YwUNaDQ1A6SB9EEi+box9HUm6IWnRK3G5PsNGAEZ2gkZ+QDy4YMaavQxLZWxiiuE1FQy4uCqvLICmTBbOybV2qVLP7PsfR+aXpS7N8L1WcSg6mo/woNxw127d4wSDQ+Xylqi3iVgLjGGmRba6QjVNXtyg6QpGWgyJ3R9KIWC7Vh1pYQkiM0chPxJ0Jek8OZRgSjZPkdntjKMvCYrUZpHGU2AHAjnzKanfJgc5yu7LrDMj+LlifQX3uFpU5AV0ljXgwln6MtHg1G71wWXAYEtM4EqFLYpL221MmEhOfrEo/e0lkWnF/LRqDJB7FHWdWl0Ea25YbSSOybV2SXI7HKy3ueB28dXIPdinEFdL0zoTldkVS5UCdMFCpn7TcrlDBE/UuWbNEa0gL2ThlckJdy0Fhty8Tp+rf4XZ74mEA6no7m3ugBoHyOzGdFav5htVmkMY2hwQAuzamj1eWAGDmZH1Hfg3365CLNfgucjs9ZUIdq90VbKlCuSObdkfQmcgkojQZbPE53/RxSBQxMdW5WziBRJ8YEz8X5tNxyv5YytItltsV/KWMKFW+inWM5PJidRmk8dR8E3YX0s+/WQOAw6XGro38BgkAjMhEsPsiPuo7EzdMh+rcLe64UM4umzdFp5ZRFp2I954g0KwlOD0BqL0xy+2IO7KUrljVai+L2RNAWinSXef0U7sCoFmr6WiV3Uqh6XQEC+zlsT2fOjRXB4AjZ5q/+M4pkZcKvV6rllBrIqDeu6stVzy1U5xK0pR2CrbRpKfacsX300yCaiSSeFP29WtiNCXMglWxBFqD1HIFI2mUVlZ8sirlde52e4K10vSUG0L68jFndRmkUs0WLxGXxLVbp2ZO1t1u70jZvKY4KfJSRjQkO1quO2p6YQlOxWZCcCANUHcWrzQRBIoq89KU0S1Rp5k6OCY4tgoA8ukYjdNstUXDaJSQIgrBQG4uZVAe8lZqmwZYbQZpnHNIuwvpXMr4999/M5cyROXLRiS4EYDlco6jZUW8HwQ91VZb8LRB3ydGWW21+IBjeldM3EgAXWLMEusmRaAJkZ2q2YLORC4Zo6nqMzHCaLlkjOY8hOKgjy2ryyCJJ9il8rlf2vKVH8ze+YtFwdcJ3cLcbk+k8w1Q792mLZqsol+r0myr+XFRdLeUncVNW3QoO+VbFV+IEOqnkiI5QYNEL0gTb01Ld+0JjbRnomqNYzdOLFasfNATwW5dsvmdf7Z1dyH9oV3rZS9Uroum5SmjW+WGIz5ZlTqSpqgnv+kg7N1vdRZPB/+ZMp9+3nQ2ZUQXyiSiof6E1e6KfyJKOWW5juCg0zThFpxCSaCMQ4pPkxpnVuwH82Q8G9kN8i/ekxd/h6EuRc1yBXu3UGLaHWUnR/FyS8rxE2R8oshClBChgeD1QDnwGyXdQvMKpZolvp8akQnKZvOCBy/K70SwlRQhk4hSSvXG+VQtyOoySOKapWVBqEuh7HtAKValKcJF6clPGUZD0Woa0fCWgyidrk7XnZmT9dA/K9VEmxpQQmrgBF+kMBmnKcdBycABxagLq90VV3hSfi3luqKWgxeFFfvBRnG7vRU82GqQ0MIdFK0BTSTN7fbEF6Lp6KOyJ794s2oAyKdjoR8KxSAVsvF1a8JfxGoL9Z4g0GRcUJooUpamWa7ocEig668qri8HuuOU2+1ZbndspcLirCKDNM4SO1xCz48oHexptAYLS+31FLuhOCifiLIfEop/mTQioVuq6biZuOhCGzOxLkV0q9ywcXRiYaWdKGcUSoOEotSg/FrEPaToxETNCvn2UFq4jjOryCC5nd4KdnUHiU6EekgIDU9pEC+VpwSlJz9lx1hls0HFtSdAJhzS9UPKJFQED1C+PcpfCsXFpEmCLiy1xSUhNM1V3Y5QR/bxZ1Vs0ASVRZoXl9CUr2l3xMVvSSMSKkBCOdDRZLBV9uQvN2yUDHZ4ZgIj4gQAVJU0GEcHmpAdikGiycBhRUSMyETodY7iitH4WNWWKy7nG2dWk0GyFMmCLzqhITuk5EQidFY6isaXJoOtctKV5XbFN4X1a8IDnljN6WnmY6F8gbmUEWplBSer9gkN2WF9e8WpZOh1jlXFFYrKkfYXhVVkkFZ2l9xBQiPsgqPY6HG7CGFSmoQBllyFZi2301OjZccyEjQLUf6lOChnlEzcCG1viOVf0oDyS1HW261sVpFBMh1FzTfHgYArW3yuM4EubIJTiRwaM6FJcdNQzCVCa0GUqWNQJCHkhw7e6VAkIUAt0McoKw7/WrCmcRrRkN75QDe+NhQalVCpaq/gRnawqgwSStn2siCTiAZc2Vi7T3EqERoLQqlEzqdjod6t+AhXQqiHZNodIzIh/qGSsUjo7CXErGfIfCwMIwHUAn2Uyy/UmRDvzkcoTiVCP5R4/3JC6Az4FV+4sooMkjIJ8kUnn44FNFctN2x1HbgxdGI0eyXWYJHQURdYyapClmrEnBrEC0jpqbbagpNVgW6WnTIxJOCpsdWUpo0zq8ggUQ5rWRkEfFiVgjSsOH5oVlnljYqSQKLRKGJNSwkVvxG3T3whALDDPhTKL0W+ltDqbyyhQfBCltt1uz01Udy5c5aapl8Xi1VkkFZPYWw+HQ/I7qC08wGKJmlYEadM3KDJIaFsqaFxfDQPiUI6iOJMABG/BdarYrWSSsYiNGd88YVoXgprsmqoOSf/ihGaDneaV3zzs1VkkKx2d2UXOffJxKMB2R0siepUMsQgvXLaRDl259PhM2mwThuhSicsI0EDlvEL3SuxxlYVsvHgb4+078JKywdrDcp1h7IpeDBTSSM421cx2yg3VCYerrJb8cWUq8kguSt5Fv0gwbcHmW6JstDCUvCBrpdGOgGEqiewik5o2gCiXEU0ncWxChUK2YQalR2EeS3Ef0KJ4obGPBHHNAT/EFUL6dwQnYCw8KCy2ODFYhUZpEqzraY5ykWH3B5+t2vNwhmBbEQmgre5ct3ZPKkoVY6l8aVpA4h2HA5T2VkuQg9poGgPOm/iOBOhlgYxWVWYjAdnFhH9y9CiApSFSGIsuOqg3FDUlP1isboM0soWqAySNCJ+G2uliaM2LOYSwduc5XZRolvkVwteCyuUEdoGEG2bC6viIk4hyvYd+s1g2YnCZIjbV2kqzeOirFWYDIlDlqq2sigulr58bFlFBkmlDPSiE+ADIeZFw9ItOF84zUya5TVZFSiquFDmQRByKSMk3dJAUOcDcZoD3T6328tncO7B0MwiljMR6l8iBjxD0aKGFQI59q7s33KQgDgDVnuS0Juw2nKxbtTgrbnacjOJKJYzETwCQFnXMpR5EAQjMhEsfsPqIR16eC/XHax+kpm4EWzRTbuD4kyEXleVJo6oAei6mK/sU/VqMUjlhrOCp1qNUsgm/Ao8sXa6ZCwSHJ85VbPXI908hWw8QNGHKD0K1bKXaph2IiA8iDI3iLA57MqvNHFUduHOhIMmsQsVLGBJB4u5RPB1jtgZMhmLBDRrKDeclW2NYPUYJNPuqGmIOSYEHOsst4uS/8/EjfD4DNL9Y0QnAntPqMtMIE7VCj4OV5rq+jqbDo4zESqaL9cdrOkJhWxcWQfh4OscZRgSIZcyAuKQq6FwZbUYpHLDVtb6dxxIxiJ+WwNWWj70RbAqcEPXslyE8U6E0JhJ1WqLT3ElBOsaEGcb5tOxYIF+qYqTbgl9w263h+UhGZGJAEEauchxlBphphoxiptLxgLikCr7l18sVvjH67Pi1SlD+LVKs9wuVmViqE4M8UbNxIMiaYgLhfarRpTGBIfsSlW0loPBjotpdxCrW4I/FGI/yeBrGPOSMCKhalIs+Umw841V2zDOrCKDtKpCduBTYUcSsCjnrEI2HpxVxmrnAwDJWCTAICEK+sMLY+tocfzgRq6IzkQwWJ1vCMHzOxCbVQe3dCs3HKxPFBqHRIwEBMchV0Pzs9VikBAvmmVBIetd4ImlpwKqhjRoPaSDfztEgxQ6BtdyFcXxa5Y7lURTagSc8VVuc4h2IlgPWW21sa690EurXEcYaU8IjkMqU3heRFaLQUIsEV8WJIyI51GrVLMQd59gPbEyx8W0cWZP0CyEGMQPDg8ixgaDJ4tbbgdLpgEARmQioKNPtYV2LgzWQ+JWwYdcFQ6a2xfsEyOGHMaW1WKQ5s2V7+0OYkQmOj2PW8jt9rDy/wBQzCWCt1TEGzXg5Gi7OI3sQsHd5oIFFKH1KEwE7Ke4H6o4lQxwMStNtNa0wUYCS/NNCI5DIo60D4lDYkwXG3NWi0Gy2l2Ubl3LBb8dHDEFQvA7peKe5oJvVMQtNZOIBp+7lWUiEZ2J4DhkxVSnL0espAk256aN6fZBoEVHjHkG6yFXw4jRVbRHo/SpXC74tbKvNNtYvVvIKn7lQbiZCSMyMW8GOROIO0JgBS5mJjL4jI/oTCR94rcE3HLLAAkAmT2BZdGN6ES15bt31yzMPmFJI+Jn0RG7hECYgGLFz56Ai2KQXn755c9+9rPXXXfdLbfcct9991WrVconmqZ53XXXvfDCCxyLlmqrqw6pMOkt+8aNz+TTvn0nKyZm063Q2CBustcvC1JuOIgN44M7iysry+/0MH3ZgAmHxKFB1JcHKDWqLRcrrQgAhWzCz6LjGongDvpaZYfPgQMHbr311qeeeuryyy+3bfuBBx7Yt29fpVKhfO7CwkK7zTNx0u2s8DkiQ+SSMc9+J51eD3cEst+mULWQA0GBOjG7OJXEWiifjgU0Ske8igI6i593JvDO3QFuX7muaJtzOz1EZ6KQDTqjAJ7lI/hdErgqIb/yQYLK4ZAXC6UG6cyZM/fcc8/69eufeOKJBx988ODBg/v37z99+vR9990X/MS5ubkHH3zwa1/7GvfS5Ya9SoYhEfIZ79BTue4ghi4DdGK4CdhgnRhuC+RMIurnIVVbLqIzUcjG/eKQ550JpC8wOA7pdntYLQch8JJAV74F5JBKNRuxd2U+HQu4JBAPXkZ0ImDY/GoYoKPUIB04cMBxnP3792/dupU8cvfdd2/ZsuWxxx6r1+t+z7rhhhv27t37R3/0R5blW0UYiuV2V5WHlE/HPG9XXK8/ILGMayRyKSOglx3uYBEjMuGXMDiF11kVAvdu4kxgSswN6a2kCJsyvpcEcgYucL4q1gRhQiYR9XNccL+9pOG7EOiQHTqHDx8GgJtvvnnwwT179nQ6nZdeesnvWV/+8pe//vWvf/3rX7/22mu5l14NNWWDJI2IZ0VIuYFWxAeBtbHopzm/Mz5i+21CcSrpaye6vSk8KxtsznG/Pb+cIqBOCQH/egM4L8VG1JcHzVdFLMomqDHnhWw8IA65GjYxpQZpZmYmnU5v2LBh8MGrr74aAF5//XW/Z910000f/vCHP/zhD2/ZsoV76dXg7Q7hmTPHzYJsnkr4qVSxBtkRAnrZIbbf7qPGTgR4SLjOBBAPyadldbmBOfA0oGQHsf12KLjOREBXCMSibILftbcarBEAqEurdLtdx3E2b9489HgqlQKAUqmEvuLOnTvPLx1bAx/8KvrrjzmequJS1cbNIfmdHHFb5QdkQdCPGkbUV2uA60zk03G/2CCuMwEAhSwZ6pMe/SdcO5E0IgEaRVwrm0/H/Kws7vY9lTTmznnnC2qWGzpuih5ytPLs6s13nff3wOWCOoPU6XQAoJ896hOJRPr/isuxY8fIf8wuWv/kP/0t+uuPOaQcckc+1X+E5HVw0/KvnjY9/wlXDks2F89iW3TpUSHr2xUCN+CZifuqJ9CdiQChHa7KLkDLvtBs4zoTmUTU08qi1+vkUsaJBd8zypWXZrAWIj9Eue5sWzesGuUz5/09EJaJcVIt++52g5SakkCcMbyMGG3UiO71J42I3+6DuymQ3dm79wT2uTvAYOPaiQB9uQxnwtMgmXYHt11s0oj6zTzF1bkA6Y3tNQBFhkHyO6OgpzD9yqURh36NM+oMUiwWi0ajr7322tDjpK4okZDYRsF03FVVFUvIpYYnuqLXWgbokks1C7E2CAD8thj0jpMB7VtwnYkAfbkMZ8LTIMkwEn5ZENz+cgCQSRiea6FfewE/OnoKM5+OeQZyV0kWXOk2vWHDhkajMfRgrVYDgCuuuELeuuW6s6qm8xFyKWPoAI4uPUoaUb+yCXQ74TdDAd2ZWO/juKA7E+Cvxka3E+vX+NSlYZ9RAuuQHMTaIPCfHoR+7WUSUb9sH/rQPD9ftmKuijCPUoM0PT3tOM5gWBMAnn/+eQDYvn27vHUD+uGvYFJGZOj8iFuvA4GF5eidifPpmKfxm284uG1z/T6UjCOqnxq7VLORlRo+auxKE/msFqxlx21N69ejD72DcCGbCFBq4K6VT/tefiu+kR0oNkikAunRRx/tP1Kv159++ulCoXDVVVeRR6rV6szMzKuvvoq4bqlqr/i27aPkM/FWW27IzohOeB6HSQ9N3DBpJm74JQzQ926/UVLoO4Jf57dKE9nt89MaVJpt9BE7fkI7t9PDdSY2T3lLzCvNNu79nklE1bSfAH8Xc950sAY2jjNKDdLevXsvu+yyhx9++JFHHmm323Nzc5/5zGdardZdd93V/5uZmZnbbrvtk5/8JOK6bre3CnNI79m0ZunCK3u+4WxCriFNlBvO6PYtyZnwSRggOxPFXFcafXMAABsoSURBVMJzsniliTaEtE8mbni6fdWWi+tMZOJRT4U0bhkAwc/4KWs0UGkij2nIpbyTVeSEh2vR/VxM3KZfY4vSbToSidx///3bt2+/9957d+/evXfv3kOHDt1555379u2TvfRq+C2HGC1RLFUxG3z1Gb1/0LPKAJA0Ij7d+ZDdX79TPG4jO4LfcRg96+lnJGQUq/oNjUVfy09+gtvhCQDIcXY02yejL63fV/R3Z8w1sZWfQ1LtA15yySWPP/74s88+u7i4aBjG9PT0UGXSTTfdNJRk6vOVr3zlK1/5Cseip2o2YvHacqGQHT7pl2oWumHOp2OjW6qMs3AuZXi2Iq0027hy2EzcOz4jxZnwOg4T9QS6nfCMQ86bzuZJ5A9Fhsbu2rhm8EHcbqcEvynm5Qb+ZFXSHHJHIjX4ILqiAQDS8WjTy5ft9npXF7O4a40hFyEoGYlEbrzxRpUrut3eKvSQCtnhQ/HsIr7jkklESzVrsPwW5CiC1q+JnaqrkMNmElHPmInblTLBZNT4kayShFS5dyDofRK2udFsn4zx234lzJUm/igpIvIcus7LDfwwWi5lnFhojT6+SvIOK/8TwvlD9Kr4pINkElG32xs8F8u4fzzT8uWGg95frjAZHz0OW27X7fZGy9pFIFmB0bVO1WzcDBz4GL9qy5XhTHjqy2Ws5ZntQy85gLcafo9SbuCXeRSyidEBtZbbRY/i+ulcZDjoY8iq2KZNu7Pi27Z7sm1d8njl/GlLRgIWAHKp2OhOh9tZlWBEJkazIETgq2ZTQFdPAMCmTHw0C1K12tvWpTz/nptCNl6qeTgTVUuGUsOj4FeGzsVvCiXu7AmC5/glGUKDQtbDlyW9ZtCv8zFkVRgky+2gX6DLAhLNJ/89uzgcWENaIjFqkKw2/vQpz5FIpZolo3w9GYuMCu2sdhfdmfCsgpSRKjciE57NGsp15O4JAFCY9IikoY+Zh4AplBJSmJ5nlFN1G//gFfUYzW7andVQhASrxCDJcOGXBYPaKknzj/Pp2KjWoNzw6A4piGcG27Q76M4EEJH0iJUt1axiDls6GIuM7t3lhiOjJt+z4Be9NA0A1q/xEL/NN/CdZk+3j6gnJEQCPGKeuC3tzy+U9LB8q6RvEKwSg2S18UO9y4JM/O1D8ZFyc3fBY/qAIJ6NdmQYv0I2MZoqn120ZARjC5MJ78QYdibS0+2bl1OvU5xKjqgu7Uwiiu+4+KQV0UUNRmRitKGUjHgdAGyeSoyK39A7qwJAPuMRAy/VrExce0grBRk59mXB1rXJvmJHkva9kI17zaWVEjMp1ewh7bKkWsvROGS15Zp2R4YkxCOM1nBk/FKje3epir+fAkAxl/A6oyB34CaMGj8Z6gl4S/Y99GCpiu80k7Dw0FUhQ6M4nqz8bZrsLKtT1LBr05q+qOHQXO2aLfgCX89KmlLVltFlx4gMh9fnTeRGdv21huKQxPKh+9nb1qVGDRJxXHAXAmKQzAu+PdN25ezdHoMHyw05xm8qWape4PbNLlro4WLwSYxJOg+N3lOn6lJK2seQlW+Q/Lp/rgZ2bUwfPdsk/y0pkZaMRYZEQaWanYxFZByHi7lhx+V4pSVjR1ifjo00ubBkONmFbLzacodczNlFKR8qlzKGCrkkJas8VXalqhTHZVRiLqMCF7w6i5t2x2rj1y/D+Q914S9Vd3DHkYwtK98gmc4q1XwDwK6Na2YXLdPumHZHUi+sQnY4PlOq4suj31pr+EatNJ18Bn+t4lRiNDaIXlMMb2VBhr5ASUqNqeTwlLlqy5XhTHgOHpShcwGvRq7SQtPDfU9mz0lxxQAgn44PlUvPnsPvsTKerHyDVGmu0gQSABiRiWIuUarZh+Zquy9Jy6gOHu2vKumAD14D2aotV4adGG0BIGk/hZEknNvtlRvIc4P6Cw1Z2RMLra1r8T9UJhE1ohdUFpPkvwxhkWc7EimWb6TrIPrkjj7knhpaS4fsVg6Z1eHterK7kJk5WX/2H6p7tq+V8fok5jO7+PbhUZKRAE87UZdz7o5FhlyxUzVbmkFKDGZByFhVGZG0fGY0Dmlvk2CQYESsMbvYknTA3zyVGIpDyih4AoDiVKLSbF/Q9wR76lKfZGx4bGOpastwmseQVWGQ0BsgLiM+uHPdd/9+4Znj1Q9snZS0xFAkbU5aKKM4lTg1UHdCTpEyvLFt61JDBS7VlovehJSQSxmz5wYMkhwTC6QNrnvBNoc+N2hwrcEzfqXZlnQPjpYHyeg0Txg6D8kL2W1dmzxVGwpNtyWtNW6sfIO0ehSTnvzKL2z41k/mn3uzetNl6yQtMaQ1KNVsGYEgAEjHo4N2Qt5dSmTfF57xpVnZXGJuwCDJaCBNyKfjg44sSFNIA8ksDpxR5H17hezwh/r/27v3qKjKfg/gP4cBBgaRqyJqYOAWBQ2MkpSMvKe1vBxX5eo9rGOKRyNvHfVU2rJlK1tlYl6qlUtL7bI6C7DSs4oyL8gLiQLyJtSgcIRAQBphZmBmYK7nj13jOJhwzjuPM/vh+/lLnr1n7Wf/9vb57cvzPLtJw+peNsj/ts4aWqOF0RfzXHrZiX35BshsnPzvZLfFNkCO5R1FDfY780LKmRdS2E09Ehsa4NwoNGl6WMxRRESxYQrnFwa1agO7CTiC/H2cd6q+w8hoW+GBtz1Jq2/vZvS2IDZU0dppcm5S2eWJ2DCFc/RYTNMgGhly22hfcb5BRqPgY8MC6ttvzcNdq2b1HNJlxFit2jBApmmggZCQGD1TlpCMuFBGL5BE8REBzs/xa9UGRgGPjwh0vu5m1yLQHy+W/9hWt8XG7rrb5Rr/uq4nLpxJOg/y91HIZY5tqfVmha+M0bVaXHiA8zcU2D3dig1TNGluDZeuVRviIwLv/pN/ZlvOR6r2poFR5x2X+8tatXGAvECiAZKQ3Ps9Y3ARE6pwXON39VhbO02M/v+4tAg1bQZ2D9ZjwwIcY4qrWvRB/j6MbjFjw267xq9VG9nt1MgQf8c1vqpNz+hGlnodqaqWroSh7p+2iogUctnIEH/HtmrVRnY75TzvyZ/bYpL84iNunXtEVN2qHyAvkGggJKQbXUxGU4NDbJjC8WpH1WYI8vdhduXo59yfuL6je2wkq8vh+IgAx2QNXSYLu2YuPiLQue1WtelZXuPfaulUbQZGSYKIRg5ROD/dutbe7fIBWTeKDb31JK1Jw+r+km7PE+LYPnbdT+jPPjtEpNab2e2Ut+E/ITGabBEc4iMCHS2CxmiOjwhgN5VtbOit5Neq62HXpMaFBzjeV9Wqjew2JDZqjlc7tWojizlwRQlDAx0vJ2raDIksN+RouxlN4eogjrQT/113k+EdUoTS13Ev26TtiQ1TsDvPE4YG1qoN4r+Z3vZ5G/5barxDYs35wxCqNkPqKFb9y4koPiJQnAzJYrPXqo3srruTopSqtj9aBKZtt1w2KDZMUdWiJ6ImbU+IgskgJFFc+K3bPkbzsInErxWLWVbdZY5k+U7eeRhprdrArkttbGiAeJiIqKq1K3Uk6/P8j9NP1aZPHuH+WSi904BISAOnj4pHiP2axOdO5U2dicNYtd3k9IUn8UM+7Npu59s+RtOjOSRFBYlZtlXHZCYkpw0pyxp14r9r1Qam190JQwPF5ruyuZPd/SURjR0aWKc2EpHFZq9q1TO9l7XY7OItZnljZ8IwVhdDRBQXHiAOBqhVG9V68wCZpoEGQkKCeyAtZsj5Bi0Rna3rYDrpVuIw5T+au4ioqqWLadejqGA/Tfcf055WXmfbpEYF+5X+piOi8w3atJgh7DaUMExZ1aK32OxdPdaqVj3Te9mk4UGVzZ1EVPqbjt2gbCKKjwgQzz1VmyFC6cv0hXHyiCDxUeT5Bh2jjuwiR6+QskYdu8cAXoj/hKTpNg+Qb1t50OSY4ML/0WiMlvr27vTRIew2lD465GxdBxEV1mnSRzNsuxVyWVKU8nyDtknbo9abk0cEsdvWgyMHizcup652TI1luFNir5Cyxs7zDdqEoYFMx+fFhiqKr2mJ6HyDlukpkToquL6jW2O0sN4QEaWPDvn7NQ0RVTZ3Mt1WWswQcUOqNgPrnfIqAyAh4ZEdexlxoT9eaf/vX9Tpo0OYRjtpuFJjtNSqjQU1N+ePD2e3ISJKHx2S//PvP15pTx8dwvSLwzPHhFVe7+q22P5+TcP0DkkuG5QRF3q2rqOsqZPphojomeRhZ+s6NEaL6oaBXTcNIlLIZcnRg3+82l7aoGP3qk80+b7g/6q8UdWqD/LzSRrOcFsJQwM1Rkvl9a7ieu2cBFZzrHgh/hMS3APJI4K6eqz/+sUvyycPZ7ohuWzQ3ITwNV9dqVUbk6PZvun9l4lDD19sOVja8lgc20vU2DBF0nDlEwf+EaH0ZT3iZEFSxOflrR/9dH3GGIZjpYkoabiySdOz4ZurGfGhrHsVZcSHHCptyfu5bWFSJNMNzU0IL2vs3HX2tyUPDGV6jUJEz6YM21PUWNaowx0SV5o0PQN5Lrt7QyGXnVj+wMvTY/72YBTrbf3n9JgC1c2dT8azbubSYoIThgYazNYVk6OZboiItsyMPVvX8db8ONbN3N8ejKpq1de3dz85PoLphuSyQW/Njzt8sWXLzBimGyKif39kRIHqZkZcKOvXLSEB8hfTRx6+2JI9dSTTDRFR9tSRhy+2rJgcPaCGUQ6y2+19ryVBY8eOrampIaJB/3Havmu6p6sD0AeLzc46G4m6eqwKX9k929a9GXQhdj+5N7NWSnR6TEeT6M04f9uPQUggFfcmQ9Cfn7DibFv3MkNIMRtJBeeRRY8GAACp4D8hDeTPxQIASAjvCamb1acqAQDAvThPSOL7W0/XAgAA+sZ5Y63WmwdUp0kAAOniPCHhDgkAQCo4b6zVejPTORABAMBdOE9IehPukAAApIHzxhrvkAAApILzhAQAAFLBeULqttju2YwsAADwz+A8IbXqTHhkBwAgCZwnJCKS++AOCQBAAjhPSGq9CVMHAQBIAucJqctkxeSqAACSwHlC6jbbMA4JAEASOG+sWzvRqQEAQBp4TkhqvRlf5wMAkAqeE1JXjxXfLwcAkAqeE5LFZseoWAAAqeA5IbV29kQN9vd0LQAAoF94TkgaoyUkQO7pWgAAQL/wnJDQqQEAQEJ4TkgAACAhPCckDEICAJAQnhNSt9nmL+d5BwEAeIL2GgAAvALPCemmAZ0aAAAkg+eEhJkaAAAkhOeEBAAAEsJzQmrS9owcgpkaAACkgeeEZLHa8f1yAACp4DkhdZksQX6YOggAQBp4TkiYOggAQEJ4TkjdZpunqwAAAP3Fc0Jq7TQpfHneQQAAnnDbXtsHyYgIc9kBAEgFtwnJqhji6SoAAMD/AbcJyU4yfJ0PAEBCuE1IloAQJCQAAAnhNiGJ75AAAEAquG21LYrQkUMUnq4FAAD0lzQealVUVBw+fLisrCwyMjI9PT0rKyskJKSP38h85DLMGwQAIBkSuEPKz89funTpqVOnxo0b19PTc/DgwUWLFqnV6j5/GBU80Pt8jx071tNV8DwEAREgBEEivD0htbS0vP766+Hh4d9+++2hQ4cKCgpWrVrV3Ny8c+fOu//Q4h+swPfLAQCkw9ub7Pz8fJPJtGrVqpiYGLFkw4YNo0aNOnHihE6nu8sP7TJfjIoFAJAQb09I5eXlRPTEE084F2ZkZFit1gsXLtzlh3aZ3B93SAAA0uHtTXZZWZlSqYyMjHQunDRpEhFduXLlLj+0+AfjDgkAQEK8upedzWYzmUwjRoxwKQ8ICCCipqamu/w2rO6Hnau+3t1zt8d6AwHe5RKCgAgQEYIgBV6dkKxWKxE53h45yGQyx9K/cq28kF3FAADA7bz9kR0R2Wz4rBEAAP+8OiH5+vr6+PhUV1e7lJvNZiLy9/f3RKUAAIAJr05IRBQZGdnZ2elSqNVqiWjixImeqBEAADDh7QkpOTnZZDLV1NQ4FxYXFxNRXFychyoFAADu5+0JSRyBlJeX5yjR6XSnT5+OiopKSUnxXL0AAMDNvD0hzZ49e8yYMUePHv3ss8/MZnNDQ8OaNWuMRuP69es9XTUAAHCnQXa73dN16ENLS8vy5cvr6uocJdnZ2WvXrvVglQAAwO0kkJCIyGazFRYWtre3y+Xy5OTk3iOTAABA6qSRkAAAgHve/g4JAAAGCCQkAADwCkhIAADgFbx6ctX/n4qKisOHD5eVlUVGRqanp2dlZYWEhHi6Uu6Un59/5syZS5cuiV08lixZ8uijjzqv0J8I8BElg8GwdOnShISEt99+22UR90H4/fffP/zww4sXL968eTM1NTUzMzM1NdVlnT53UNIRIKL8/PyCgoLLly+HhYVNmjRpxYoVsbGxzitweRp0dXXNnj373XffnTJlSu+lbjnonooJb50a8vPzX331VblcnpaWdv369WvXrkVHR+fm5kZERHi6au7x4osvnjx5Ui6XT5kyxWg0VlRUWK3WrKysjRs3iiv0JwLcRGnTpk3Hjx+fNWvW/v37ncu5D8Kvv/66bNmyjo6OxMTEwMBA8TR45513FixY4Finzx2UdASI6IUXXjh16pSfn19aWlpbW5tKpVIoFAcOHJg8ebK4Aq+nwZEjR3bs2HHgwIHHHnvMZZFbDronY2LnSHNzc1JS0iOPPFJfXy+W5OTkCIKwefNmz1bMXU6cOCEIwsKFC2/evCmW1NXVpaenC4JQWlpq718EuImSGA1BELKzs53LuQ+CxWKZMWNGUlLS+fPnxZKqqqqUlJSUlBS9Xi+W9LmDko6A3W7/+uuvBUF47rnnOjs7xZJz586NGzdu2rRpVqvVzulpUF9ff/DgwYkTJwqCcPbsWZelbjnono0JVwlp3759giAcOXLEuXDGjBnjxo3TarWeqpUbrVmzRhCE8vJy58Ljx48LgvDyyy/b+xcBPqLU3NyckpKyYsWK3gmJ+yCIR3zfvn3Ohe+9954gCN999534Z587KOkI2O32DRs2CIJw7tw558KVK1cKgvDzzz/beTwNxEtPh94JyS0H3bMx4apTQ3l5Of05/Z1DRkaG1Wq9cOGChyrlTj/99JOfn5/4BXeHCRMmEJE4J3p/IsBBlGw22/r16wcPHrx9+/beS7kPwpkzZ4ho8eLFzoXr1q2rqamZO3eu+GefOyjpCDj4+Pg4/6lQKOjPz9Pwdxps27Zt9+7du3fvTktLu+MKbjnono0JVwmprKxMqVRGRkY6F4rN95UrVzxUKXfKycn56KOPXAp/+eUXIhJfOfYnAhxEaf/+/ZWVlbt27RoyZEjvpdwH4dKlS0qlMjo62mazXbp0qaioSKfTuazT5w5KOgJEJHbk+eabbxwlDQ0NRUVFERERycnJxONpMHPmzHnz5s2bN2/UqFF3XMEtB92zMeGnl53NZjOZTCNGjHApDwgIIKKmpiZPVMrNXHrTEVF7e3tOTg4RPfPMM/2JAAdRKisre//991etWpWammowGFyWDoQg3LhxY8yYMadPn962bVtbW5tYOGvWrO3bt4eFhVE/dlDqESCiRYsW1dbWHjx4UKVSJSYmms3mkydPKpXKDz74QCaTDYTTwIVbDrrHY8JPQrJarUTUe5o7mUzmWMqZgoKCt956q7W1df369RMmTBCfVNw9AlKPkk6n27RpU2Ji4rp16+64Qn92UNJBMJvNVqtVo9GsXr16ypQpGzdutNlsubm5J0+ebG5uzsvLk8lkfe6gpCPgkJKSEhERoVKpVCqVWJKeni7ePXB/GvTmloPu8Zhw9ciOiGw2m6ercC9UV1c/++yz69at6+np2bVr1+rVqx2L+hMB6UZp69atWq12z5494v+Qv8J3EIiotbX1+eef/+STTxYsWLBo0aIvvvhi+vTp1dXVubm5jnX63EFJR+Crr77Kzs6Ojo4+duxYTU1NVVXVli1bioqKnn76accDTO5Pg97cctA9GBN+EpKvr6+Pj091dbVLuXjf4O/v74lKMbFz587FixdXV1dnZWX98MMPTz75pFjenwhIOkqFhYXff//93LlzGxsbS0pKSkpKSktLiai9vb2kpER848p9EMTX+D4+Pi6fBFu2bBkRiQHpcwclHQEistlsb7755uDBgw8dOpSYmEhEvr6+mZmZL730UmNj49GjR7k/DXpzy0H3eEz4eWRHRJGRke3t7S6FWq2WiCZOnOiJGrmfOBhwzpw5r7zyyvDhw12W9icC0o2SRqMhovz8/Pz8fOfy8vLyZcuWKZXKiooK4j0IMplMqVTKZDKX1iEpKYmI9Hq9+GefOyjdCBDR5cuXOzs758yZExwc7Fw+b968HTt2iN18+D4N7sgtB92zMeHnDomIkpOTTSZTTU2Nc2FxcTERxcXFeahS7nTkyJFTp05lZmbu3bu3dzai/kVAulGaNGnSO7cTu30nJiY6/k28B4GIHnroIaPR6PJA/+rVq0SkVCrFP/vcQUlHQNxNi8XiUi52chEf53J/GvTmloPu2ZhwlZDEvvN5eXmOEp1Od/r06aioqJSUFM/Vy20+//xzhUKxefPmv1qhPxGQbpRGjRq14HZPPfUUEUVHRy9YsMDx6JLvIBDR1KlTLRbLxx9/7Fz45ZdfiovEP/vcQUlHID4+Pjw8vLCwsKWlxblcfIU2fvx4GgCnQW9uOegejgnrkbf3ktVqnT9/viAIn376qclkqq+vz8zMFATh2LFjnq6aG3R0dAiCMHny5H+7kwMHDtj7FwGeoqTX63vP1MB9ELq7ux9//HFBEPbs2VNZWVlaWvraa68JgrBkyRJx1hx7P3ZQ0hGw2+25ubmCIEybNi0vL0+lUhUXF7/xxhtiiTiZEMenwZYtW+44U4NbDrpnY8Lb5KotLS3Lly+vq6tzlGRnZ69du9aDVXKXoqKiFStW/NXShQsXijNe9ycC3ETJYDCkpKT0nlyV+yDU19evXLmyoaHBUfLwww/v2bNHHIck6nMHJR0BIjp69GhOTo7RaHSUjB8/fu/evY5xo7yeBlu3bs3Nzb3j5KpuOegejAlvCYmIbDZbYWFhe3u7+HWG3n3qudefCHAfJe6DYLPZiouL29raZDLZ/fff/8ADD9xxnbvvoKQjQERdXV0lJSU6nc7X1/e+++7r/UyJ+9OgN7ccdE/FhMOEBAAAUsRVpwYAAJAuJCQAAPAKSEgAAOAVkJAAAMArICEBAIBXQEICAACvgIQEAABeAQkJAAC8AhISAAB4BSQkAADwCkhIAADgFZCQAADAKyAhAQCAV0BCAgAAr4CEBAAAXuF/AYj2zGYPzALFAAAAAElFTkSuQmCC\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"iv = [0, 0.01, 0];\n",
"[t, xa] = ode45(f, t0, iv);\n",
"plot(t, xa(:,3))"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Matlab",
"language": "matlab",
"name": "matlab"
},
"language_info": {
"codemirror_mode": "octave",
"file_extension": ".m",
"help_links": [
{
"text": "MetaKernel Magics",
"url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
}
],
"mimetype": "text/x-octave",
"name": "matlab",
"version": "0.11.0"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment