diff --git a/examples/second_example.cpp b/examples/second_example.cpp index 364f9e2..3edb79b 100644 --- a/examples/second_example.cpp +++ b/examples/second_example.cpp @@ -17,7 +17,9 @@ int main(int argc, char *argv[]) { vector alpha(n * m, 1); vector field(n * m, 0); - field[0] = 2000; + field[n * 19] = 2000; + field[n * 19 + 19] = 2000; + // for (int i = 1; i<20; i++) { // for (int j = 0; j<20; j++ ) { // field[i] = 0; @@ -40,12 +42,13 @@ int main(int argc, char *argv[]) { input_param.setDomainSize(n, m); BoundaryCondition bc(n, m); boundary_condition bc_constant; - bc_constant.type = BC_TYPE_CONSTANT; + bc_constant.type = BC_TYPE_CLOSED; bc_constant.value = 0; bc.setSide(BC_SIDE_LEFT, bc_constant); - bc.setSide(BC_SIDE_TOP, bc_constant); bc.setSide(BC_SIDE_RIGHT, bc_constant); bc.setSide(BC_SIDE_BOTTOM, bc_constant); + // bc_constant.value = 2000; + bc.setSide(BC_SIDE_TOP, bc_constant); input_param.setBoundaryCondition(bc); // int iterations = 1000; diff --git a/proto/FTCS.ipynb b/proto/FTCS.ipynb index dd4ca8a..f3efad0 100644 --- a/proto/FTCS.ipynb +++ b/proto/FTCS.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 12, + "execution_count": 109, "metadata": {}, "outputs": [], "source": [ @@ -17,7 +17,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 124, "metadata": {}, "outputs": [], "source": [ @@ -36,7 +36,7 @@ "delta_x = domain_size['x'] / grid_size['x']\n", "delta_y = domain_size['y'] / grid_size['y'] \n", "\n", - "alpha_x = np.ones((grid_size['x'],grid_size['y']))\n", + "alpha_x = np.random.random((grid_size['x'],grid_size['y']))\n", "alpha_y = np.ones((grid_size['x'],grid_size['y']))\n", "\n", "# time dimension\n", @@ -51,26 +51,25 @@ }, { "cell_type": "code", - "execution_count": 29, + "execution_count": 125, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQUAAAD4CAYAAADl7fPiAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAANR0lEQVR4nO3df6hf9X3H8edrSbSQ6dQ5U39krXRBiGVkJcSVuaGztVFkaUfpImPLNkFXJqwwGG6DWrp/OoaTDcWSdkE7WrVsyxpo/BGygRVaa5T4q9WZSYq5TZO16Uxdu7roe3/ck3I/N9+b3H1/3O/33jwfcPme8/l8vud8Dl945Zzz/ea8U1VI0nE/Ne4JSJoshoKkhqEgqWEoSGoYCpIay8c9gV7OyJn1NlaOexrSkvU//Ddv1I/Tq28iQ+FtrOSKXDPuaUhL1hO1e84+Lx8kNQYKhSQbk7yUZF+S23r0n5nkwa7/iSTvHGR/kkav71BIsgy4G7gOWAvcmGTtrGE3Ad+vql8A7gT+qt/9SVoYg5wpbAD2VdUrVfUG8ACwadaYTcB93fI/Atck6XlzQ9JkGCQULgZenbF+oGvrOaaqjgGvAT/ba2NJbk6yJ8me/+XHA0xL0iAm5kZjVW2tqvVVtX4FZ457OtJpa5BQmAJWz1i/pGvrOSbJcuBngO8NsE9JIzZIKDwJrElyaZIzgM3AjlljdgBbuuUPA/9a/l9taaL1/eOlqjqW5FbgEWAZsK2qXkjySWBPVe0A/h74hyT7gCNMB4ekCZZJ/If77JxX/qJRGp0najdH60jPbwIn5kajpMlgKEhqGAqSGoaCpIahIKlhKEhqGAqSGoaCpIahIKlhKEhqGAqSGoaCpIahIKlhKEhqGAqSGoaCpIahIKlhKEhqDFIhanWSf0vyjSQvJPnjHmOuSvJakr3d38cHm66kURuk6vQx4E+q6ukkZwFPJdlVVd+YNe4rVXXDAPuRtID6PlOoqoNV9XS3/APgm5xYIUrSIjOUewpdNelfAp7o0f3eJM8keSjJ5SfZhmXjpAkwyOUDAEl+Gvgn4GNVdXRW99PAO6rq9STXA/8CrOm1naraCmyF6Ue8DzovSf0Z6EwhyQqmA+HzVfXPs/ur6mhVvd4t7wRWJDl/kH1KGq1Bvn0I0xWgvllVfzPHmLcfLz2fZEO3P2tJShNskMuHXwF+B3guyd6u7c+Bnweoqk8zXT/yo0mOAT8CNltLUppsg9SSfBzoWXZqxpi7gLv63YekhecvGiU1DAVJDUNBUsNQkNQwFCQ1DAVJDUNBUsNQkNQwFCQ1DAVJDUNBUsNQkNQwFCQ1DAVJDUNBUmPgZzQKHvn23nmN+8BF60Y6D2kYPFOQ1DAUJDUGDoUk+5M815WF29OjP0n+Lsm+JM8mec+g+5Q0OsO6p3B1VX13jr7rmK71sAa4Arine5U0gRbi8mET8Lma9jXgnCQXLsB+JfVhGKFQwKNJnkpyc4/+i4FXZ6wfoEfNScvGSZNhGJcPV1bVVJILgF1JXqyqx/6/G7FsnDQZBj5TqKqp7vUwsB3YMGvIFLB6xvolXZukCTRoLcmVSc46vgxcCzw/a9gO4He7byF+GXitqg4Osl9JozPo5cMqYHtXLnI58IWqejjJH8JPSsftBK4H9gE/BH5/wH1KGqFMYmnHs3NeXZFrxj0Nacl6onZztI70LPvoLxolNQwFSQ1DQVLDUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNQwFSY2+QyHJZV2puON/R5N8bNaYq5K8NmPMxweesaSR6vvBrVX1ErAOIMkyph/bvr3H0K9U1Q397kfSwhrW5cM1wH9U1beGtD1JYzKsUNgM3D9H33uTPJPkoSSXz7UBy8ZJk2HgR7wnOQP4NnB5VR2a1Xc28FZVvZ7keuBvq2rNqbbpI96l0Rr1I96vA56eHQgAVXW0ql7vlncCK5KcP4R9ShqRYYTCjcxx6ZDk7enKRyXZ0O3ve0PYp6QRGahsXFc/8v3ALTPaZpaM+zDw0STHgB8Bm2sSS1JJ+gnLxkmnIcvGSZo3Q0FSw1CQ1DAUJDUMBUkNQ0FSw1CQ1DAUJDUMBUkNQ0FSw1CQ1DAUJDUMBUkNQ0FSw1CQ1DAUJDUMBUkNQ0FSY16hkGRbksNJnp/Rdl6SXUle7l7PneO9W7oxLyfZMqyJSxqN+Z4p3AtsnNV2G7C7q+Owu1tvJDkPuB24AtgA3D5XeEiaDPMKhap6DDgyq3kTcF+3fB/wwR5v/QCwq6qOVNX3gV2cGC6SJsgg9xRWVdXBbvk7wKoeYy4GXp2xfqBrkzShhnKjsavlMNCz4q0lKU2GQULhUJILAbrXwz3GTAGrZ6xf0rWdoKq2VtX6qlq/gjMHmJakQQwSCjuA498mbAG+1GPMI8C1Sc7tbjBe27VJmlDz/UryfuCrwGVJDiS5CfgU8P4kLwPv69ZJsj7JZwGq6gjwl8CT3d8nuzZJE8qycdJpyLJxkubNUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNQwFSQ1DQVLDUJDUMBQkNU4ZCnPUkfzrJC8meTbJ9iTnzPHe/UmeS7I3yZ4hzlvSiMznTOFeTiz1tgt4d1X9IvDvwJ+d5P1XV9W6qlrf3xQlLaRThkKvOpJV9WhVHetWv8Z0kRdJS8Aw7in8AfDQHH0FPJrkqSQ3n2wjlo2TJsPyQd6c5C+AY8Dn5xhyZVVNJbkA2JXkxe7M4wRVtRXYCtN1HwaZl6T+9X2mkOT3gBuA3645KspU1VT3ehjYDmzod3+SFkZfoZBkI/CnwG9U1Q/nGLMyyVnHl5muI/l8r7GSJsd8vpLsVUfyLuAspi8J9ib5dDf2oiQ7u7euAh5P8gzwdeDLVfXwSI5C0tBYS1I6DVlLUtK8GQqSGoaCpIahIKlhKEhqGAqSGoaCpIahIKlhKEhqGAqSGoaCpIahIKlhKEhqGAqSGoaCpIahIKlhKEhqGAqSGv2WjftEkqnu+Yx7k1w/x3s3Jnkpyb4ktw1z4pJGo9+ycQB3duXg1lXVztmdSZYBdwPXAWuBG5OsHWSykkavr7Jx87QB2FdVr1TVG8ADwKY+tiNpAQ1yT+HWrur0tiTn9ui/GHh1xvqBrq0ny8ZJk6HfULgHeBewDjgI3DHoRKpqa1Wtr6r1Kzhz0M1J6lNfoVBVh6rqzap6C/gMvcvBTQGrZ6xf0rVJmmD9lo27cMbqh+hdDu5JYE2SS5OcAWwGdvSzP0kL55RVp7uycVcB5yc5ANwOXJVkHdOl5vcDt3RjLwI+W1XXV9WxJLcCjwDLgG1V9cIoDkLS8Fg2TjoNWTZO0rwZCpIahoKkhqEgqWEoSGoYCpIahoKkhqEgqWEoSGoYCpIahoKkhqEgqWEoSGoYCpIahoKkhqEgqWEoSGoYCpIa83lG4zbgBuBwVb27a3sQuKwbcg7wX1W1rsd79wM/AN4EjlXV+qHMWtLInDIUmC4bdxfwueMNVfVbx5eT3AG8dpL3X11V3+13gpIW1ilDoaoeS/LOXn1JAnwE+PUhz0vSmAx6T+FXgUNV9fIc/QU8muSpJDefbEOWjZMmw3wuH07mRuD+k/RfWVVTSS4AdiV5sStYe4Kq2gpshelHvA84L0l96vtMIcly4DeBB+caU1VT3ethYDu9y8tJmiCDXD68D3ixqg706kyyMslZx5eBa+ldXk7SBDllKHRl474KXJbkQJKbuq7NzLp0SHJRkp3d6irg8STPAF8HvlxVDw9v6pJGwbJx0mnIsnGS5s1QkNQwFCQ1DAVJDUNBUsNQkNQwFCQ1DAVJDUNBUsNQkNQwFCQ1DAVJDUNBUsNQkNQwFCQ1DAVJjYl8yEqS/wS+Nav5fGAp1o9YqscFS/fYlsJxvaOqfq5Xx0SGQi9J9izFClNL9bhg6R7bUj2u47x8kNQwFCQ1FlMobB33BEZkqR4XLN1jW6rHBSyiewqSFsZiOlOQtAAMBUmNRREKSTYmeSnJviS3jXs+w5Jkf5LnkuxNsmfc8xlEkm1JDid5fkbbeUl2JXm5ez13nHPsxxzH9YkkU93ntjfJ9eOc47BNfCgkWQbcDVwHrAVuTLJ2vLMaqqurat0S+N77XmDjrLbbgN1VtQbY3a0vNvdy4nEB3Nl9buuqameP/kVr4kOB6UrV+6rqlap6A3gA2DTmOWmWqnoMODKreRNwX7d8H/DBhZzTMMxxXEvaYgiFi4FXZ6wf6NqWggIeTfJUkpvHPZkRWFVVB7vl7zBddHipuDXJs93lxaK7LDqZxRAKS9mVVfUepi+N/ijJr417QqNS0999L5Xvv+8B3gWsAw4Cd4x1NkO2GEJhClg9Y/2Srm3Rq6qp7vUwsJ3pS6Wl5FCSCwG618Njns9QVNWhqnqzqt4CPsMS+9wWQyg8CaxJcmmSM4DNwI4xz2lgSVYmOev4MnAt8PzJ37Xo7AC2dMtbgC+NcS5DczzoOh9iiX1uy8c9gVOpqmNJbgUeAZYB26rqhTFPaxhWAduTwPTn8IWqeni8U+pfkvuBq4DzkxwAbgc+BXwxyU1M/1f4j4xvhv2Z47iuSrKO6cuh/cAt45rfKPgzZ0mNxXD5IGkBGQqSGoaCpIahIKlhKEhqGAqSGoaCpMb/AZlpw/ulk6QYAAAAAElFTkSuQmCC", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa4AAAGdCAYAAABKG5eZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAkE0lEQVR4nO3df3RU9Z3/8dfFhAlyklGUJDMQQuQglB+HQkACKj+kBENFWKmg7oGwtrZsqRVTTiG2HnH/aLCtLgdBWbv8kLWLnG4IZDdsJRzzQ0tgQRLXWsS4RpKVpBw4JQO4DAE+3z/8ZuqYmcGBmSSf+Hycc8/x3vv+fPKeTyKv3Jk7E8cYYwQAgCV6dXUDAABEg+ACAFiF4AIAWIXgAgBYheACAFiF4AIAWIXgAgBYheACAFgloasbiJUrV67oxIkTSk5OluM4Xd0OACAKxhidPXtWXq9XvXpFvqbqMcF14sQJZWRkdHUbAIDr0NTUpIEDB0as6THBlZycLEm6S7OVoMQu7gYAEI1LatPb2hP4tzySHhNc7U8PJihRCQ7BBQBW+f+fmvtVXurh5gwAgFUILgCAVeIWXC+99JKysrKUlJSk7OxsvfXWWxHrq6qqlJ2draSkJN12223auHFjvFoDAFgsLsG1Y8cOLV++XD/72c9UW1uru+++W3l5eWpsbAxZ39DQoNmzZ+vuu+9WbW2tnnrqKf34xz9WcXFxPNoDAFjMiccfkpw4caLGjRunl19+OXDsG9/4hubNm6eioqIO9StXrlRpaamOHj0aOLZ06VK9++67qqmp+Upf0+fzye12a5rmcnMGAFjmkmlTpXartbVVKSkpEWtjfsV18eJFvfPOO8rNzQ06npubq/3794ccU1NT06F+1qxZOnz4sNra2kKO8fv98vl8QRsAoOeLeXCdOnVKly9fVlpaWtDxtLQ0tbS0hBzT0tISsv7SpUs6depUyDFFRUVyu92BjTcfA8DXQ9xuzvjyvfjGmIj354eqD3W8XWFhoVpbWwNbU1PTdXYMALBBzN+AfOutt+qGG27ocHV18uTJDldV7dLT00PWJyQk6JZbbgk5xuVyyeVyxaZpAIA1Yn7F1bt3b2VnZ6u8vDzoeHl5uSZPnhxyzKRJkzrU7927V+PHj1diIjdaAAD+Ki5PFRYUFOif//mftXnzZh09elRPPvmkGhsbtXTpUkmfP823ePHiQP3SpUt1/PhxFRQU6OjRo9q8ebM2bdqkFStWxKM9AIDF4vJZhQsXLtTp06f1D//wD2pubtaoUaO0Z88eZWZmSpKam5uD3tOVlZWlPXv26Mknn9SGDRvk9Xq1bt06zZ8/Px7tAQAsFpf3cXUF3scFAPbq0vdxAQAQTwQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqMQ+uoqIiTZgwQcnJyUpNTdW8efN07NixiGMqKyvlOE6H7YMPPoh1ewAAy8U8uKqqqrRs2TIdOHBA5eXlunTpknJzc3X+/Pmrjj127Jiam5sD29ChQ2PdHgDAcgmxnvD3v/990P6WLVuUmpqqd955R1OmTIk4NjU1VTfddFOsWwIA9CBxf42rtbVVktSvX7+r1o4dO1Yej0czZsxQRUVFxFq/3y+fzxe0AQB6vrgGlzFGBQUFuuuuuzRq1KiwdR6PR6+88oqKi4u1c+dODRs2TDNmzFB1dXXYMUVFRXK73YEtIyMjHg8BANDNOMYYE6/Jly1bprKyMr399tsaOHBgVGPnzJkjx3FUWloa8rzf75ff7w/s+3w+ZWRkaJrmKsFJvK6+AQCd65JpU6V2q7W1VSkpKRFr43bF9fjjj6u0tFQVFRVRh5Yk5eTkqL6+Pux5l8ullJSUoA0A0PPF/OYMY4wef/xxlZSUqLKyUllZWdc0T21trTweT4y7AwDYLubBtWzZMv3rv/6rdu/ereTkZLW0tEiS3G63+vTpI0kqLCzUp59+qm3btkmS1q5dq8GDB2vkyJG6ePGiXnvtNRUXF6u4uDjW7QEALBfz4Hr55ZclSdOmTQs6vmXLFi1ZskSS1NzcrMbGxsC5ixcvasWKFfr000/Vp08fjRw5UmVlZZo9e3as2wMAWC6uN2d0Jp/PJ7fbzc0ZAGChbnFzBgAA8UBwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsEvPgWr16tRzHCdrS09MjjqmqqlJ2draSkpJ02223aePGjbFuCwDQQyTEY9KRI0dq3759gf0bbrghbG1DQ4Nmz56txx57TK+99pr+8Ic/6Ic//KH69++v+fPnx6M9AIDF4hJcCQkJV73Kardx40YNGjRIa9eulSR94xvf0OHDh/XrX/+a4AIAdBCX17jq6+vl9XqVlZWlhx56SB9//HHY2pqaGuXm5gYdmzVrlg4fPqy2traw4/x+v3w+X9AGAOj5Yh5cEydO1LZt2/TGG2/oN7/5jVpaWjR58mSdPn06ZH1LS4vS0tKCjqWlpenSpUs6depU2K9TVFQkt9sd2DIyMmL6OAAA3VPMgysvL0/z58/X6NGj9a1vfUtlZWWSpFdffTXsGMdxgvaNMSGPf1FhYaFaW1sDW1NTUwy6BwB0d3F5jeuL+vbtq9GjR6u+vj7k+fT0dLW0tAQdO3nypBISEnTLLbeEndflcsnlcsW0VwBA9xf393H5/X4dPXpUHo8n5PlJkyapvLw86NjevXs1fvx4JSYmxrs9AIBlYh5cK1asUFVVlRoaGnTw4EF95zvfkc/nU35+vqTPn+JbvHhxoH7p0qU6fvy4CgoKdPToUW3evFmbNm3SihUrYt0aAKAHiPlThf/7v/+rhx9+WKdOnVL//v2Vk5OjAwcOKDMzU5LU3NysxsbGQH1WVpb27NmjJ598Uhs2bJDX69W6deu4FR4AEJJj2u+EsJzP55Pb7dY0zVWCw1OMAGCTS6ZNldqt1tZWpaSkRKzlswoBAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAViG4AABWIbgAAFYhuAAAVol5cA0ePFiO43TYli1bFrK+srIyZP0HH3wQ69YAAD1AQqwnPHTokC5fvhzY/+Mf/6iZM2fqwQcfjDju2LFjSklJCez3798/1q0BAHqAmAfXlwNnzZo1GjJkiKZOnRpxXGpqqm666aZYtwMA6GHi+hrXxYsX9dprr+nRRx+V4zgRa8eOHSuPx6MZM2aooqIinm0BACwW8yuuL9q1a5fOnDmjJUuWhK3xeDx65ZVXlJ2dLb/fr3/5l3/RjBkzVFlZqSlTpoQd5/f75ff7A/s+ny+WrQMAuinHGGPiNfmsWbPUu3dv/fu//3tU4+bMmSPHcVRaWhq2ZvXq1Xr22Wc7HJ+muUpwEqPuFQDQdS6ZNlVqt1pbW4Pudwglbk8VHj9+XPv27dP3vve9qMfm5OSovr4+Yk1hYaFaW1sDW1NT07W2CgCwSNyeKtyyZYtSU1P17W9/O+qxtbW18ng8EWtcLpdcLte1tgcAsFRcguvKlSvasmWL8vPzlZAQ/CUKCwv16aefatu2bZKktWvXavDgwRo5cmTgZo7i4mIVFxfHozUAgOXiElz79u1TY2OjHn300Q7nmpub1djYGNi/ePGiVqxYoU8//VR9+vTRyJEjVVZWptmzZ8ejNQCA5eJ6c0Zn8vl8crvd3JwBABbqFjdnAAAQDwQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCoEFwDAKgQXAMAqBBcAwCpRB1d1dbXmzJkjr9crx3G0a9euoPPGGK1evVper1d9+vTRtGnT9P7771913uLiYo0YMUIul0sjRoxQSUlJtK0BAL4Gog6u8+fPa8yYMVq/fn3I87/85S/1wgsvaP369Tp06JDS09M1c+ZMnT17NuycNTU1WrhwoRYtWqR3331XixYt0oIFC3Tw4MFo2wMA9HCOMcZc82DHUUlJiebNmyfp86str9er5cuXa+XKlZIkv9+vtLQ0Pffcc/rBD34Qcp6FCxfK5/PpP//zPwPH7r33Xt18883avn37V+rF5/PJ7XZrmuYqwUm81ocEAOgCl0ybKrVbra2tSklJiVgb09e4Ghoa1NLSotzc3MAxl8ulqVOnav/+/WHH1dTUBI2RpFmzZkUc4/f75fP5gjYAQM8X0+BqaWmRJKWlpQUdT0tLC5wLNy7aMUVFRXK73YEtIyPjOjoHANgiLncVOo4TtG+M6XDsescUFhaqtbU1sDU1NV17wwAAayTEcrL09HRJn19BeTyewPGTJ092uKL68rgvX11dbYzL5ZLL5brOjgEAtonpFVdWVpbS09NVXl4eOHbx4kVVVVVp8uTJYcdNmjQpaIwk7d27N+IYAMDXU9RXXOfOndNHH30U2G9oaFBdXZ369eunQYMGafny5frFL36hoUOHaujQofrFL36hG2+8UY888khgzOLFizVgwAAVFRVJkp544glNmTJFzz33nObOnavdu3dr3759evvtt2PwEAEAPUnUwXX48GFNnz49sF9QUCBJys/P19atW/XTn/5U//d//6cf/vCH+stf/qKJEydq7969Sk5ODoxpbGxUr15/vdibPHmyXn/9df385z/X008/rSFDhmjHjh2aOHHi9Tw2AEAPdF3v4+pOeB8XANiry97HBQBAvBFcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrEFwAAKsQXAAAqxBcAACrRB1c1dXVmjNnjrxerxzH0a5duwLn2tratHLlSo0ePVp9+/aV1+vV4sWLdeLEiYhzbt26VY7jdNguXLgQ9QMCAPRsUQfX+fPnNWbMGK1fv77Duc8++0xHjhzR008/rSNHjmjnzp368MMPdf/991913pSUFDU3NwdtSUlJ0bYHAOjhEqIdkJeXp7y8vJDn3G63ysvLg469+OKLuuOOO9TY2KhBgwaFnddxHKWnp0fbDgDgaybur3G1trbKcRzddNNNEevOnTunzMxMDRw4UPfdd59qa2sj1vv9fvl8vqANANDzxTW4Lly4oFWrVumRRx5RSkpK2Lrhw4dr69atKi0t1fbt25WUlKQ777xT9fX1YccUFRXJ7XYHtoyMjHg8BABAN+MYY8w1D3YclZSUaN68eR3OtbW16cEHH1RjY6MqKysjBteXXblyRePGjdOUKVO0bt26kDV+v19+vz+w7/P5lJGRoWmaqwQnMerHAgDoOpdMmyq1W62trVfNi6hf4/oq2tratGDBAjU0NOjNN9+MKrQkqVevXpowYULEKy6XyyWXy3W9rQIALBPzpwrbQ6u+vl779u3TLbfcEvUcxhjV1dXJ4/HEuj0AgOWivuI6d+6cPvroo8B+Q0OD6urq1K9fP3m9Xn3nO9/RkSNH9B//8R+6fPmyWlpaJEn9+vVT7969JUmLFy/WgAEDVFRUJEl69tlnlZOTo6FDh8rn82ndunWqq6vThg0bYvEYAQA9SNTBdfjwYU2fPj2wX1BQIEnKz8/X6tWrVVpaKkn65je/GTSuoqJC06ZNkyQ1NjaqV6+/XuydOXNG3//+99XS0iK3262xY8equrpad9xxR7TtAQB6uOu6OaM78fl8crvd3JwBABaK5uYMPqsQAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYJWog6u6ulpz5syR1+uV4zjatWtX0PklS5bIcZygLScn56rzFhcXa8SIEXK5XBoxYoRKSkqibQ0A8DUQdXCdP39eY8aM0fr168PW3HvvvWpubg5se/bsiThnTU2NFi5cqEWLFundd9/VokWLtGDBAh08eDDa9gAAPVxCtAPy8vKUl5cXscblcik9Pf0rz7l27VrNnDlThYWFkqTCwkJVVVVp7dq12r59e7QtAgB6sLi8xlVZWanU1FTdfvvteuyxx3Ty5MmI9TU1NcrNzQ06NmvWLO3fvz/sGL/fL5/PF7QBAHq+mAdXXl6efvvb3+rNN9/U888/r0OHDumee+6R3+8PO6alpUVpaWlBx9LS0tTS0hJ2TFFRkdxud2DLyMiI2WMAAHRfUT9VeDULFy4M/PeoUaM0fvx4ZWZmqqysTA888EDYcY7jBO0bYzoc+6LCwkIVFBQE9n0+H+EFAF8DMQ+uL/N4PMrMzFR9fX3YmvT09A5XVydPnuxwFfZFLpdLLpcrZn0CAOwQ9/dxnT59Wk1NTfJ4PGFrJk2apPLy8qBje/fu1eTJk+PdHgDAMlFfcZ07d04fffRRYL+hoUF1dXXq16+f+vXrp9WrV2v+/PnyeDz65JNP9NRTT+nWW2/V3/zN3wTGLF68WAMGDFBRUZEk6YknntCUKVP03HPPae7cudq9e7f27dunt99+OwYPEQDQk0QdXIcPH9b06dMD++2vM+Xn5+vll1/We++9p23btunMmTPyeDyaPn26duzYoeTk5MCYxsZG9er114u9yZMn6/XXX9fPf/5zPf300xoyZIh27NihiRMnXs9jAwD0QI4xxnR1E7Hg8/nkdrs1TXOV4CR2dTsAgChcMm2q1G61trYqJSUlYi2fVQgAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwCsEFALAKwQUAsArBBQCwStTBVV1drTlz5sjr9cpxHO3atSvovOM4Ibdf/epXYefcunVryDEXLlyI+gEBAHq2qIPr/PnzGjNmjNavXx/yfHNzc9C2efNmOY6j+fPnR5w3JSWlw9ikpKRo2wMA9HAJ0Q7Iy8tTXl5e2PPp6elB+7t379b06dN12223RZzXcZwOYwEA+LK4vsb15z//WWVlZfrud7971dpz584pMzNTAwcO1H333afa2tqI9X6/Xz6fL2gDAPR8cQ2uV199VcnJyXrggQci1g0fPlxbt25VaWmptm/frqSkJN15552qr68PO6aoqEhutzuwZWRkxLp9AEA35BhjzDUPdhyVlJRo3rx5Ic8PHz5cM2fO1IsvvhjVvFeuXNG4ceM0ZcoUrVu3LmSN3++X3+8P7Pt8PmVkZGia5irBSYzq6wEAutYl06ZK7VZra6tSUlIi1kb9GtdX9dZbb+nYsWPasWNH1GN79eqlCRMmRLzicrlccrlc19MiAMBCcXuqcNOmTcrOztaYMWOiHmuMUV1dnTweTxw6AwDYLOorrnPnzumjjz4K7Dc0NKiurk79+vXToEGDJH3+tN3vfvc7Pf/88yHnWLx4sQYMGKCioiJJ0rPPPqucnBwNHTpUPp9P69atU11dnTZs2HAtjwkA0INFHVyHDx/W9OnTA/sFBQWSpPz8fG3dulWS9Prrr8sYo4cffjjkHI2NjerV668Xe2fOnNH3v/99tbS0yO12a+zYsaqurtYdd9wRbXsAgB7uum7O6E58Pp/cbjc3ZwCAhaK5OYPPKgQAWIXgAgBYheACAFiF4AIAWIXgAgBYheACAFiF4AIAWIXgAgBYheACAFiF4AIAWIXgAgBYJW5/j6urlHz4nlKSry+PZ3m/GZtmAOBr4I0Tddc9h+/sFd18+1er5YoLAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYBWCCwBgFYILAGAVggsAYJUe8xeQjTGSJN+5K9c91yXTdt1zAMDXhe/s9f+72/5vd/u/5ZH0mOA6e/asJClz3CcxmO3jGMwBAF8PN98eu7nOnj0rt9sdscYxXyXeLHDlyhWdOHFCycnJchwnZI3P51NGRoaampqUkpLSyR1eO/rufLb2Tt+di75jxxijs2fPyuv1qlevyK9i9Zgrrl69emngwIFfqTYlJaXbfLOiQd+dz9be6btz0XdsXO1Kqx03ZwAArEJwAQCs8rUKLpfLpWeeeUYul6urW4kKfXc+W3un785F312jx9ycAQD4evhaXXEBAOxHcAEArEJwAQCsQnABAKzS44LrpZdeUlZWlpKSkpSdna233norYn1VVZWys7OVlJSk2267TRs3buykTj9XVFSkCRMmKDk5WampqZo3b56OHTsWcUxlZaUcx+mwffDBB53UtbR69eoOXz89PT3imK5e63aDBw8OuX7Lli0LWd9V611dXa05c+bI6/XKcRzt2rUr6LwxRqtXr5bX61WfPn00bdo0vf/++1edt7i4WCNGjJDL5dKIESNUUlLSaX23tbVp5cqVGj16tPr27Suv16vFixfrxIkTEefcunVryO/BhQsXOqVvSVqyZEmHr5+Tk3PVebtyvSWFXDfHcfSrX/0q7Jydsd7Xo0cF144dO7R8+XL97Gc/U21tre6++27l5eWpsbExZH1DQ4Nmz56tu+++W7W1tXrqqaf04x//WMXFxZ3Wc1VVlZYtW6YDBw6ovLxcly5dUm5urs6fP3/VsceOHVNzc3NgGzp0aCd0/FcjR44M+vrvvfde2NrusNbtDh06FNR3eXm5JOnBBx+MOK6z1/v8+fMaM2aM1q9fH/L8L3/5S73wwgtav369Dh06pPT0dM2cOTPwuZ2h1NTUaOHChVq0aJHeffddLVq0SAsWLNDBgwc7pe/PPvtMR44c0dNPP60jR45o586d+vDDD3X//fdfdd6UlJSg9W9ublZSUlKn9N3u3nvvDfr6e/bsiThnV6+3pA5rtnnzZjmOo/nz50ecN97rfV1MD3LHHXeYpUuXBh0bPny4WbVqVcj6n/70p2b48OFBx37wgx+YnJycuPV4NSdPnjSSTFVVVdiaiooKI8n85S9/6bzGvuSZZ54xY8aM+cr13XGt2z3xxBNmyJAh5sqVKyHPd4f1lmRKSkoC+1euXDHp6elmzZo1gWMXLlwwbrfbbNy4Mew8CxYsMPfee2/QsVmzZpmHHnoo5j0b07HvUP7rv/7LSDLHjx8PW7Nlyxbjdrtj21wEofrOz883c+fOjWqe7rjec+fONffcc0/Ems5e72j1mCuuixcv6p133lFubm7Q8dzcXO3fvz/kmJqamg71s2bN0uHDh9XW1jV/2qS1tVWS1K9fv6vWjh07Vh6PRzNmzFBFRUW8W+ugvr5eXq9XWVlZeuihh/Txx+E/Vb87rrX0+c/Na6+9pkcffTTshzO36+r1/qKGhga1tLQEranL5dLUqVPD/rxL4b8PkcbEW2trqxzH0U033RSx7ty5c8rMzNTAgQN13333qba2tnMa/ILKykqlpqbq9ttv12OPPaaTJ09GrO9u6/3nP/9ZZWVl+u53v3vV2u6w3uH0mOA6deqULl++rLS0tKDjaWlpamlpCTmmpaUlZP2lS5d06tSpuPUajjFGBQUFuuuuuzRq1KiwdR6PR6+88oqKi4u1c+dODRs2TDNmzFB1dXWn9Tpx4kRt27ZNb7zxhn7zm9+opaVFkydP1unTp0PWd7e1brdr1y6dOXNGS5YsCVvTHdb7y9p/pqP5eW8fF+2YeLpw4YJWrVqlRx55JOKHvQ4fPlxbt25VaWmptm/frqSkJN15552qr6/vtF7z8vL029/+Vm+++aaef/55HTp0SPfcc4/8fn/YMd1tvV999VUlJyfrgQceiFjXHdY7kh7z6fDtvvxbszEm4m/SoepDHe8MP/rRj/Tf//3fevvttyPWDRs2TMOGDQvsT5o0SU1NTfr1r3+tKVOmxLtNSZ//T9xu9OjRmjRpkoYMGaJXX31VBQUFIcd0p7Vut2nTJuXl5cnr9Yat6Q7rHU60P+/XOiYe2tra9NBDD+nKlSt66aWXItbm5OQE3Qhx5513aty4cXrxxRe1bt26eLcqSVq4cGHgv0eNGqXx48crMzNTZWVlEYOgu6y3JG3evFl/+7d/e9XXqrrDekfSY664br31Vt1www0dfpM5efJkh9942qWnp4esT0hI0C233BK3XkN5/PHHVVpaqoqKiq/851m+KCcnp0t/G+rbt69Gjx4dtofutNbtjh8/rn379ul73/te1GO7er3b7+CM5ue9fVy0Y+Khra1NCxYsUENDg8rLy6P+0xq9evXShAkTuvR74PF4lJmZGbGH7rLekvTWW2/p2LFj1/Tz3h3W+4t6THD17t1b2dnZgTvE2pWXl2vy5Mkhx0yaNKlD/d69ezV+/HglJibGrdcvMsboRz/6kXbu3Kk333xTWVlZ1zRPbW2tPB5PjLv76vx+v44ePRq2h+6w1l+2ZcsWpaam6tvf/nbUY7t6vbOyspSenh60phcvXlRVVVXYn3cp/Pch0phYaw+t+vp67du375p+cTHGqK6urku/B6dPn1ZTU1PEHrrDerfbtGmTsrOzNWbMmKjHdof1DtJVd4XEw+uvv24SExPNpk2bzJ/+9CezfPly07dvX/PJJ58YY4xZtWqVWbRoUaD+448/NjfeeKN58sknzZ/+9CezadMmk5iYaP7t3/6t03r++7//e+N2u01lZaVpbm4ObJ999lmg5st9/+M//qMpKSkxH374ofnjH/9oVq1aZSSZ4uLiTuv7Jz/5iamsrDQff/yxOXDggLnvvvtMcnJyt17rL7p8+bIZNGiQWblyZYdz3WW9z549a2pra01tba2RZF544QVTW1sbuPtuzZo1xu12m507d5r33nvPPPzww8bj8RifzxeYY9GiRUF31f7hD38wN9xwg1mzZo05evSoWbNmjUlISDAHDhzolL7b2trM/fffbwYOHGjq6uqCfub9fn/YvlevXm1+//vfm//5n/8xtbW15u/+7u9MQkKCOXjwYKf0ffbsWfOTn/zE7N+/3zQ0NJiKigozadIkM2DAgG693u1aW1vNjTfeaF5++eWQc3TFel+PHhVcxhizYcMGk5mZaXr37m3GjRsXdFt5fn6+mTp1alB9ZWWlGTt2rOndu7cZPHhw2G9svEgKuW3ZsiVs388995wZMmSISUpKMjfffLO56667TFlZWaf2vXDhQuPxeExiYqLxer3mgQceMO+//37Yno3p+rX+ojfeeMNIMseOHetwrrusd/tt+F/e8vPzjTGf3xL/zDPPmPT0dONyucyUKVPMe++9FzTH1KlTA/Xtfve735lhw4aZxMREM3z48JgHcKS+Gxoawv7MV1RUhO17+fLlZtCgQaZ3796mf//+Jjc31+zfv7/T+v7ss89Mbm6u6d+/v0lMTDSDBg0y+fn5prGxMWiO7rbe7f7pn/7J9OnTx5w5cybkHF2x3teDP2sCALBKj3mNCwDw9UBwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKxCcAEArEJwAQCsQnABAKzy/wBVnTUJwLaO3QAAAABJRU5ErkJggg==", "text/plain": [ - "
" + "
" ] }, - "metadata": { - "needs_background": "light" - }, + "metadata": {}, "output_type": "display_data" } ], "source": [ "# initialize environment with start values\n", "C_t = np.zeros((grid_size['x'],grid_size['y']))\n", - "C_t[5,4] = 2000\n", + "C_t[19,0] = 2000\n", + "C_t[19,19] = 2000\n", "\n", "\n", "plt.imshow(C_t, vmin=0, vmax=2000)\n", @@ -79,7 +78,7 @@ }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 114, "metadata": {}, "outputs": [], "source": [ @@ -229,7 +228,7 @@ }, { "cell_type": "code", - "execution_count": 32, + "execution_count": 115, "metadata": {}, "outputs": [], "source": [ @@ -328,19 +327,17 @@ }, { "cell_type": "code", - "execution_count": 33, + "execution_count": 126, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUEAAAD8CAYAAADpLRYuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAAaO0lEQVR4nO3dfaxddZ3v8ffHysOIcAGrHZ4UxltN1IxVGmDiw4WrQG2M4MRx2pkIKHeqI9zozUzmopOI0ZB474hmHAzcOjTAjRYZEW1mKlC5ZjomghQGgfIwFITYWtuUGsAnpO3n/rF+G7anZ5+zztr7dD+szytZOWv/9trrgYZvfr/fd631lW0iItrqRcM+gYiIYUoQjIhWSxCMiFZLEIyIVksQjIhWSxCMiFZLEIyIoZJ0gqTvSXpA0mZJHyvtR0vaIOmR8veo0i5JX5K0RdK9kt7cta/zy/aPSDq/1vFzn2BEDJOkY4BjbN8t6XDgLuBc4AJgt+3PSboEOMr2/5S0HPjvwHLgVODvbZ8q6WhgE7AUcNnPybZ/PtPx0xOMiKGyvd323WX9GeBB4DjgHODastm1VIGR0n6dK7cDR5ZAejawwfbuEvg2AMtmO/6LB3kxg3KwDvGhHDbs04iYWL/hl/zWz6qffZx9xmF+cvfeWtvede+zm4HfdDWttr166naSTgTeBNwBLLK9vXz1M2BRWT8O+EnXz7aWtl7tMxrJIHgoh3Gq3jHs04iYWHf4tr73sWv3Xu645fha2x50zKO/sb10pm0kvRS4Efi47aelF2K0bUual7m7DIcjoiGz1/tqLbORdBBVAPyq7W+W5h1lmNuZN9xZ2rcBJ3T9/PjS1qt9Rn0FQUnLJD1csjSXTPP9IZK+Xr6/o3R1I2ICGNiHay0zUdXluxp40PYXur5aB3QyvOcD3+5qP69kiU8DnirD5luAsyQdVTLJZ5W2GTUeDktaAHwZOJNq7H2npHW2H+ja7ELg57b/s6QVwP8C/rTpMSNitOxj9l5eDW8BPgDcJ+me0vZJ4HPADZIuBJ4A3l++W0+VGd4C/Ar4IIDt3ZI+C9xZtvuM7d2zHbyfOcFTgC22HwOQdD1V1qY7CJ4DfLqsfwO4QpKc+3Iixp4xz9UY6s66H/v7QK8kzX7JgRI/LuqxrzXAmrkcv5/hcJ1MzPPb2N4DPAW8bLqdSVolaZOkTc/xbB+nFREHgoG9uNYyykYmO1zS5asBjtDRo/1fLSIAZp3vGwf9BME6mZjONlslvRj4T8CTfRwzIkaEgb0TMLPVz3D4TmCxpJMkHQysoMradOvO7rwP+H+ZD4yYHPtqLqOscU/Q9h5JF1OloBcAa2xvlvQZYJPtdVRp7/8raQuwmypQRsQE8BjM99XR15yg7fVU6erutk91rf8G+JN+jhERo8mG58Y/Bo5OYiQixo3Y2/POlvGRIBgRjRjYl55gRLRZeoIR0VrVzdIJghHRUgae8/i/iCpBMCIaMWLvBLyNL0EwIhrb5wyHI6KlMicYES0n9mZOMCLaqnqzdIJgRLSULX7rBcM+jb4lCEZEY/smYE5w/PuyETEUVWLkRbWW2UhaI2mnpPu72r4u6Z6yPN6pPyLpREm/7vruqq7fnCzpvlLc7UvqrtvZQ3qCEdHQQBMj1wBXANd1Gmw/X5RN0uVU5Tk6HrW9ZJr9XAn8BVXx9vXAMuA7Mx04PcGIaKSTGKmzzLoveyPVO0f3U3pz7wfWzrSPUpv4CNu3l5c3XwecO9uxEwQjorG9Vq2lT28Ddth+pKvtJEn/LulfJb2ttB1HVfCtY7rib/vJcDgiGjHiOdcOIQslber6vLoUV6tjJb/bC9wOvNL2k5JOBr4l6fV1T2SqBMGIaKSTGKlpl+2lcz1GKdD2x8DJzx/Xfhaqury275L0KPAaqsJux3f9fLrib/tpPByWdIKk70l6QNJmSR+bZpvTJT3VlcX51HT7iojxY+oNhfscDr8TeMj288NcSS+XtKCs/wGwGHjM9nbgaUmnlXnE84Bvz3aAfnqCe4C/sn23pMOBuyRtsP3AlO3+zfa7+zhORIyoQT0xImktcDrVsHkrcKntq6mKs01NiLwd+Iyk56iK2X3Ediep8lGqTPPvUWWFZ8wMQ3/V5rZTjc2x/YykB6kmIacGwYiYQDYDu0XG9soe7RdM03YjcGOP7TcBb5jLsQcyJyjpROBNVPfmTPVHkn4E/BT4a9ube+xjFbAK4FBeMojTioh5VCVG8tgckl5KFZU/bvvpKV/fDbzK9i8kLQe+RTV+30/JFK0GOEJHT0D5lojJNwkvVe3rCiQdRBUAv2r7m1O/t/207V+U9fXAQZIW9nPMiBgNRuxzvWWUNe4JluzL1cCDtr/QY5vfp7rJ0ZJOoQq6TzY9ZkSMlknoCfYzHH4L8AHgvs6DzcAngVcC2L4KeB/wl5L2AL8GVpTHWSJizFV1h1scBG1/H2Z+j47tK6geio6IiaO8Xj8i2qsquZnscES0lK12D4cjIlJoKSJaq3qfYOYEI6K1UnIzIlqsukUmPcGIaKk8OxwRrZfi6xHRWtWrtDIcjogWy5xgRLRW9RaZDIcjoqWqx+YSBCOitSajJzj+VxARQ7MP1VpmI2mNpJ2S7u9q+7SkbV3VKpd3ffcJSVskPSzp7K72ZaVti6RL6lxDeoIR0ciAs8PXUL1277op7V+0/fnuBkmvo6pC93rgWOC7kl5Tvv4ycCawFbhT0rppKmD+jgTBiGhsUMNh2xtLwbY6zgGuL0XYfyxpC3BK+W6L7ccAJF1ftp0xCGY4HBGNzLHGyEJJm7qWVTUPc7Gke8tw+ajSdhzwk65ttpa2Xu0zSk8wIhoxsKd+T3CX7aVzPMSVwGfLoT4LXA58aI77mNUgSm4+DjwD7AX2TL3QUpDp74HlwK+AC2zf3e9xI2L45jM7bHtHZ13SV4B/Lh+3ASd0bXp8aWOG9p4G1RM8w/auHt+9i6rW8GLgVKrofuqAjhsRwzLP5TQlHWN7e/n4XqCTOV4HfE3SF6gSI4uBH1LVPFos6SSq4LcC+LPZjnMghsPnANeVKnO3SzpyysVFxBga5EtVJa0FTqeaO9wKXAqcLmlJOdTjwIcBbG+WdANVwmMPcJHtvWU/FwO3AAuANbY3z3bsQQRBA7dKMvB/bK+e8n2vycrfCYJlonQVwKG8ZACnFRHzbVA9Qdsrp2m+eobtLwMum6Z9PbB+LsceRBB8q+1tkl4BbJD0kO2Nc91JCZ6rAY7Q0alNHDHi8lLVwva28nenpJuo7tfpDoIzTWJGxJgyYs++8b/Lrq8rkHSYpMM768BZvDB52bEOOE+V04CnMh8YMRkG9djcMPXbE1wE3FTdBcOLga/ZvlnSRwBsX0U1Pl8ObKG6ReaDfR4zIkaBMxymPJ7yxmnar+paN3BRP8eJiNGTOcGIaL0EwYhoLSP2TkBiJEEwIhob9aRHHQmCEdGIkxiJiLZzgmBEtNf8vkDhQEkQjIjG0hOMiNayYe++BMGIaLFkhyOitUyGwxHRakmMRETLeQLe/JkgGBGNTcJwePwf/IuIoaiywy+qtcym1BXeKen+rra/k/RQqTt8k6QjS/uJkn4t6Z6yXNX1m5Ml3Sdpi6QvlWqXM0oQjIjG7HpLDdcAy6a0bQDeYPsPgf8APtH13aO2l5TlI13tVwJ/wQsVLqfucz8JghHRmK1ay+z78UZg95S2W23vKR9vpyrN0ZOkY4AjbN9e3mN6HXDubMdOEIyIRky9AFiC4EJJm7qWVXM83IeA73R9PknSv0v6V0lvK23HUVWz7OhUtpxREiMR0dgcksO7bC9tcgxJf0tVX/irpWk78ErbT0o6GfiWpNc32Tf00ROU9Nquicl7JD0t6eNTtjld0lNd23yq6fEiYsQYvE+1lqYkXQC8G/jzMsTF9rO2nyzrdwGPAq+hqmLZPWSuVdmycU/Q9sPAknKiC8rBbppm03+z/e6mx4mI0TWft8hIWgb8DfBfbP+qq/3lwG7beyX9AVUC5DHbu0tn7DTgDuA84B9mO86ghsPvoMrWPDGg/UXEGBjUzdKS1gKnU80dbgUupcoGHwJsKHe63F4ywW8HPiPpOWAf8BHbnaTKR6kyzb9HNYfYPY84rUEFwRXA2h7f/ZGkHwE/Bf7a9ubpNioTpasADuUlAzqtiJgvg3x22PbKaZqv7rHtjcCNPb7bBLxhLsfuOzss6WDgPcA/TfP13cCrbL+Rqlv6rV77sb3a9lLbSw/ikH5PKyLmmwGr3jLCBnGLzLuAu23vmPqF7adt/6KsrwcOkrRwAMeMiBEwwJulh2YQw+GV9BgKS/p9YIdtSzqFKug+OYBjRsTQ9Zf5HRV9BUFJhwFnAh/uavsIgO2rgPcBfylpD/BrYEUnzR0RE2AC/m/uKwja/iXwsiltV3WtXwFc0c8xImJEeTLeIpMnRiKiubb3BCOi7dITjIg22zfsE+hfgmBENNO5T3DMJQhGRGOTcK9HgmBENJcgGBGtluFwRLSZ0hOMiNayoO2PzUVEy6UnGBGtliAYEa2WIBgRrTUhN0un7nBENCbXW2bdj7RG0k5J93e1HS1pg6RHyt+jSrskfUnSFkn3Snpz12/OL9s/Iun8OteQIBgRzbnmMrtrgGVT2i4BbrO9GLitfIbqbfaLy7IKuBKqoElVoOlU4BTg0k7gnEmCYEQ0NqieoO2NwO4pzecA15b1a4Fzu9qvc+V24EhJxwBnAxts77b9c2AD+wfW/WROMCKaqz8nuFDSpq7Pq22vnuU3i2xvL+s/AxaV9eOAn3Rtt7W09WqfUYJgRDRTf6gLsMv20saHquoUzUsuutZweC6TltP8ds4TlRExJgY3JzidHWWYS/m7s7RvA07o2u740tarfUZ15wSvof6k5fOaTlRGxHjQvnpLQ+uATsfpfODbXe3nlSzxacBTZdh8C3CWpKNKnDmrtM2oVhCc46Rlt0YTlRExJgbUE5S0FvgB8FpJWyVdCHwOOFPSI8A7y2eA9cBjwBbgK8BHAWzvBj4L3FmWz5S2GfUzJ9hr0rJbo4nKiBh9dTO/ddhe2eOrd0yzrYGLeuxnDbBmLsceSGJkEJOWklZR3fPDobxkEKcVEfOt5U+M9Jq07FZ7otL2attLbS89iEP6OK2IOGDmNzFyQPQTBHtNWnZrNFEZEeNhUDdLD1PdW2RqT1pKWirpH6H5RGVEjAHPe3b4gKg1JzjHSctNwH/r+jznicqIGBMj3surI0+MRERzCYIR0WajPt9XR94iExGtlp5gRDQ3AT3BBMGIaMajn/mtI0EwIppLTzAi2kpMRmIkQTAimksQjIjWGoNH4upIEIyI5pIYiYg2S08wItotQTAiWmsM3hVYRx6bi4jGBvE+QUmvlXRP1/K0pI9L+rSkbV3ty7t+8wlJWyQ9LOnsfq4hPcGIaG4APUHbDwNLACQtoHr7/E3AB4Ev2v589/aSXgesAF4PHAt8V9JrbO9tcvz0BCOisXl4qeo7gEdtPzHDNucA19t+1vaPqarOndL0GhIEI6KZuvVFqt7iQkmbupZVPfa6Aljb9fliSfdKWtNVs3ygVSwTBCOiEc1hAXZ1CqmVZfV++5MOBt4D/FNpuhJ4NdVQeTtw+XxcR4JgRDQ32Gpz7wLutr0DwPYO23tt76Mqst4Z8tauYlnHrEGwdEN3Srq/q+3vJD1Uuqk3STqyx28fl3RfyexsanqSETGaBlxtbiVdQ+FOSd/ivUAnBq0DVkg6RNJJwGLgh02voU5P8Bpg2ZS2DcAbbP8h8B/AJ2b4/Rm2l9he2uwUI2JkDagnKOkw4Ezgm13N/7t0ou4FzgD+B4DtzcANwAPAzcBFTTPDUOMWGdsbJZ04pe3Wro+3A+9regIRMaYG+FJV278EXjal7QMzbH8ZcNkgjj2IOcEPAd/p8Z2BWyXdNUM2CABJqzqZo+d4dgCnFRHzbrBzgkPR183Skv4W2AN8tccmb7W9TdIrgA2SHrK9cboNS7ZoNcAROnrE/7NFBEzGCxQa9wQlXQC8G/hz29P+p7C9rfzdSXUHeOMbGiNiBE1AT7BREJS0DPgb4D22f9Vjm8MkHd5ZB87ihexOREyAAWeHh6LOLTJrgR8Ar5W0VdKFwBXA4VRD3HskXVW2PVbS+vLTRcD3Jf2IKn39L7ZvnperiIgDz1QvVa2zjLA62eGV0zRf3WPbnwLLy/pjwBv7OruIGFkptBQRkSAYEW2m6XOiYyVBMCKaGYPMbx0JghHRWOYEI6LVBvXY3DAlCEZEc+kJRkRrjcGN0HUkCEZEcwmCEdFWuVk6IlpP+8Y/CiYIRkQzuU8wItpuEm6RSbW5iGhucDVG9ivKJuloSRskPVL+HlXaJelLkraUYm9v7ucSEgQjorEBv09walG2S4DbbC8GbiufoSrNubgsq6jqEzeWIBgRzRiw6y3NnANcW9avBc7tar/OlduBI6eU55yTBMGIaEz76i3Awk4htbJMLbw2XVG2Rba3l/WfUb2oGeA44Cddv91a2hpJYiQiGpnjfYK7Zqk9vl9Rtu4vbVuan7sS0xOMiGbqDoVrDId7FGXb0Rnmlr87y+bbgBO6fn58aWukTo2RNZJ2Srq/q+3TkraVTM49kpb3+O0ySQ+XLM4l020TEeNrEImRGYqyrQPOL5udD3y7rK8DzitZ4tOAp7qGzXNWZzh8DVVhpeumtH/R9ud7/UjSAuDLwJlUY/Y7Ja2z/UDDc42IUTOYAeoi4CZJUMWkr9m+WdKdwA2luNsTwPvL9uupahltAX4FfLCfg9cptLRR0okN9n0KsKUUXELS9VRZnQTBiAkxiFm6XkXZbD8JvGOadgMX9X/kSj9zgheXGxXXdG5inGJOGRxJqzqZo+d4to/TiogDwsBe11tGWNMgeCXwamAJsB24vN8Tsb3a9lLbSw/ikH53FxEHwCQUX290i4ztHZ11SV8B/nmazQaawYmIETQB1eYa9QSn3J39XqpMzlR3AoslnSTpYGAFVVYnIiZEK3qCktYCp1Pd8b0VuBQ4XdISqlmBx4EPl22PBf7R9nLbeyRdDNwCLADW2N48HxcREUPQlldp2V45TfPVPbb9KVXquvN5PVU6OyImjACNeNKjjjw2FxGNaQLmBBMEI6KZtgyHIyKm19drskZGgmBENDbqmd86EgQjorn0BCOitZzscES03fjHwATBiGgut8hERLslCEZEaxmYgOLrCYIR0YhwhsMR0XL7xr8rmGpzEdFMZzhcZ5mBpBMkfU/SA5I2S/pYae9Z0E3SJ0oBt4clnd3PZaQnGBGNDWg4vAf4K9t3l6pzd0naUL7br6CbpNdRvZ/09cCxwHclvcb23iYHT08wIpobQN1h29tt313WnwEeZIZ6RFQF2663/aztH1NVnTul6SUkCEZEQ4Mrvt5RKlu+CbijNE1X0G1ORdxmkyAYEc3Mrdrcwk41ybKsmro7SS8FbgQ+bvtp5qGg23QyJxgRjc1hTnCX7aU99yMdRBUAv2r7mzBjQbeBFnGbtSdYuqE7Jd3f1fb1rozN45Lu6fHbxyXdV7bb1PQkI2JEDWA4LElUJTsetP2FrvZeBd3WASskHSLpJGAx8MOml1CnJ3gNcAVwXafB9p92nejlwFMz/P4M27uanmBEjCgD+waSHX4L8AHgvq4O1SeBldMVdLO9WdINwANUmeWLmmaGoV6hpY1lsnI/JYK/H/ivTU8gIsbVYN4sbfv7VHWbpupZpM32ZcBlfR+c/hMjbwN22H6kx/cGbpV013QTod0krepMmj7Hs32eVkQcEAPODg9Dv4mRlcDaGb5/q+1tkl4BbJD0kO2N021oezWwGuAIHT3a/9UiomSHW/zYnKQXA38MfL3XNra3lb87gZvo44bGiBg1Bu+rt4ywfobD7wQesr11ui8lHVYegUHSYcBZvJDdiYhJMAHD4Tq3yKwFfgC8VtJWSReWr1YwZSgs6VhJncnMRcD3Jf2IKn39L7ZvHtypR8RQdbLDdZYRVic7vLJH+wXTtP0UWF7WHwPe2Of5RcQoG/FeXh15YiQimksQjIjWsmFv43uUR0aCYEQ0l55gRLRagmBEtNfoZ37rSBCMiGYMHvEboetIEIyI5ibgsbkEwYhoxp6IkpsJghHRXBIjEdFmTk8wItpr9F+OUEeCYEQ0M7jX6w9VgmBENGLAE/DYXOoOR0QzHtxLVSUtk/SwpC2SLjkAZ/+89AQjojEPYDgsaQHwZeBMYCtwp6R1th/oe+c1pCcYEc0Npid4CrDF9mO2fwtcD5wz7+dejGRP8Bl+vuu7/sYTU5oXApNYv3hSrwsm99om4bpe1e8OnuHnt3zX31hYc/NDJW3q+ry6FFcDOA74Sdd3W4FT+z2/ukYyCNp++dQ2SZtsLx3G+cynSb0umNxrm9Trmivby4Z9DoOQ4XBEDNs24ISuz8eXtgMiQTAihu1OYLGkkyQdTFXEbd2BOvhIDod7WD37JmNpUq8LJvfaJvW6hsL2HkkXA7cAC4A1tjcfqOPLE/DYS0REUxkOR0SrJQhGRKuNRRAc5iM180nS45Luk3TPlHuoxo6kNZJ2Srq/q+1oSRskPVL+HjXMc2yix3V9WtK28u92j6TlwzzH6M/IB8GuR2reBbwOWCnpdcM9q4E6w/aSCbjv7Bpg6n1jlwC32V4M3FY+j5tr2P+6AL5Y/t2W2F5/gM8pBmjkgyBDfqQm6rG9Edg9pfkc4Nqyfi1w7oE8p0HocV0xQcYhCE73SM1xQzqXQTNwq6S7JK0a9snMg0W2t5f1nwGLhnkyA3axpHvLcHnshvnxgnEIgpPsrbbfTDXUv0jS24d9QvPF1b1Yk3I/1pXAq4ElwHbg8qGeTfRlHILgUB+pmU+2t5W/O4GbqIb+k2SHpGMAyt+dQz6fgbC9w/ZeV0V3v8Lk/bu1yjgEwaE+UjNfJB0m6fDOOnAWcP/Mvxo764Dzy/r5wLeHeC4D0wnsxXuZvH+3Vhn5x+aG/UjNPFoE3CQJqn+Hr9m+ebin1JyktcDpwEJJW4FLgc8BN0i6EHgCeP/wzrCZHtd1uqQlVMP7x4EPD+v8on95bC4iWm0chsMREfMmQTAiWi1BMCJaLUEwIlotQTAiWi1BMCJaLUEwIlrt/wPp94uV7Lx+gQAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgYAAAGiCAYAAACcWg7FAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABC7UlEQVR4nO3df1hUZd4/8PeoMCgLk4gwM4lEXkolxCKYgJWYBVJqpq2afRHLxXwyW0Iu09qKfdonrF3TK00zF8UfFD77KGqri8IqEIHmL1p/xeJGgsZIecmM2ArI3N8/XM46MhwYZwZwzvu1130tc87n3HPP0ZwP96+jEkIIEBEREQHo1d0NICIiop6DiQERERFJmBgQERGRhIkBERERSZgYEBERkYSJAREREUmYGBAREZGEiQERERFJmBgQERGRhIkBERERSZgYEBEROVFGRgZGjhwJLy8v+Pn5YfLkyaioqLCIEUIgPT0der0effv2RWxsLE6dOmUR09jYiAULFsDX1xeenp6YNGkSzp8/bxFz+fJlJCYmQqPRQKPRIDExEfX19Ta1l4kBERGRExUVFWH+/Pk4ePAg8vPzcf36dcTFxeHq1atSzAcffIAPP/wQq1atwuHDh6HVavHEE0/gypUrUkxKSgpyc3ORk5ODkpISNDQ0YMKECWhpaZFiZs6cifLycuTl5SEvLw/l5eVITEy0rcGCiIiIukxdXZ0AIIqKioQQQpjNZqHVasXSpUulmGvXrgmNRiM++eQTIYQQ9fX1ws3NTeTk5EgxFy5cEL169RJ5eXlCCCFOnz4tAIiDBw9KMWVlZQKA+Pbbbzvdvj725UE9h9lsxg8//AAvLy+oVKrubg4REdlACIErV65Ar9ejVy/ndWZfu3YNTU1NDqlLCNHm+0atVkOtVsteZzQaAQA+Pj4AgKqqKhgMBsTFxVnUM2bMGJSWluKll17C0aNH0dzcbBGj1+sREhKC0tJSxMfHo6ysDBqNBqNGjZJioqKioNFoUFpaiuDg4E59LpdJDH744QcEBAR0dzOIiMgONTU1GDRokFPqvnbtGoICfwFDXUvHwZ3wi1/8Ag0NDRbH3nnnHaSnp7d7jRACqampePjhhxESEgIAMBgMAAB/f3+LWH9/f5w7d06KcXd3R//+/dvEtF5vMBjg5+fX5j39/PykmM5wmcTAy8sLAPAwnkIflVs3t4aIiGxxXTSjBLulf8udoampCYa6FlQdDYS3l329EqYrZgRFnENNTQ28vb2l4x31Frzyyiv4+9//jpKSkjbnbu19sNYjcatbY6zFd6aem7lMYtD6ofuo3JgYEBHdiYT1LzZH8/bqZXdiINXl7W2RGMhZsGABdu3aheLiYoteEa1WC+DGb/w6nU46XldXJ/UiaLVaNDU14fLlyxa9BnV1dYiJiZFiLl682OZ9f/zxxza9EXK4KoGIiBSlRZgdUjpLCIFXXnkF27dvx/79+xEUFGRxPigoCFqtFvn5+dKxpqYmFBUVSV/6ERERcHNzs4ipra3FyZMnpZjo6GgYjUZ8/fXXUsyhQ4dgNBqlmM5wmR4DIiKizjBDwAxhdx2dNX/+fHz22WfYuXMnvLy8pPF+jUaDvn37QqVSISUlBe+99x6GDh2KoUOH4r333kO/fv0wc+ZMKXbOnDlYuHAhBgwYAB8fH6SlpSE0NBSPP/44AOD+++/H+PHjkZycjLVr1wIA5s6diwkTJnR64iHgxB6D1atXIygoCB4eHoiIiMCXX34pG19UVISIiAh4eHjg3nvvxSeffOKsphERkYKZHfS/zlqzZg2MRiNiY2Oh0+mksnXrVilm0aJFSElJwcsvv4zIyEhcuHAB+/bts5hzsXz5ckyePBnTpk3D6NGj0a9fP3zxxRfo3bu3FJOdnY3Q0FDExcUhLi4ODz74IDZv3mzT/VEJIexLm6zYunUrEhMTsXr1aowePRpr167Fn/70J5w+fRqDBw9uE19VVYWQkBAkJyfjpZdewldffYWXX34Zn3/+OaZOndqp9zSZTNBoNIhVTeYcAyKiO8x10YxCsQNGo7HTY/a2av2e+KFikEMmH+qDzzu1vd3FKYnBqFGjMGLECKxZs0Y6dv/992Py5MnIyMhoE//6669j165dOHPmjHRs3rx5+Oabb1BWVtap92RiQER05+rKxKDm27sdkhgE3HfBJRMDhw8lNDU14ejRoxabMABAXFwcSktLrV5TVlbWJj4+Ph5HjhxBc3Oz1WsaGxthMpksChERUUda5xjYW1yVwxODn376CS0tLVY3amhvgwWDwWA1/vr16/jpp5+sXpORkSE9JEKj0XBzIyIiIgdw2uRDWzdqsBZv7XirJUuWwGg0SqWmpsbOFhMRkRKYIdBiZ3HlHgOHL1f09fVF79692/QO3LxRw620Wq3V+D59+mDAgAFWr+nMftRERES36urlincah/cYuLu7IyIiwmITBgDIz89vd4OF6OjoNvH79u1DZGQk3Nw4kZCIiKirOGUoITU1FX/605+wfv16nDlzBq+99hqqq6sxb948ADeGAWbNmiXFz5s3D+fOnUNqairOnDmD9evXIzMzE2lpac5oHhERKViLEA4prsopOx9Onz4dly5dwn//93+jtrYWISEh2LNnDwIDAwHc2Maxurpaig8KCsKePXvw2muv4eOPP4Zer8dHH33U6T0MiIiIOsv872JvHa7KKfsYdAfuY0BEdOfqyn0Mvj3jDy879zG4csWM++6/6JL7GPBZCUREpCitKwvsrcNVMTEgIiJFaRE3ir11uComBkREpCicYyDPaRscERER0Z2HPQZERKQoZqjQgvZ34u1sHa6KiQERESmKWdwo9tbhqjiUQERERBL2GBARkaK0OGAowd7rezImBkREpChMDORxKIGIiIgk7DEgIiJFMQsVzMLOVQl2Xt+TMTEgIiJF4VCCPA4lEBERkYQ9BkREpCgt6IUWO38vbnFQW3oiJgZERKQowgFzDATnGBAREbkGzjGQxzkGREREJGGPARERKUqL6IUWYeccAxd+VgITAyIiUhQzVDDb2WFuhutmBhxKICIiIgl7DIiISFE4+VAeEwMiIlIUx8wx4FACERERKQB7DIiISFFuTD608yFKHEogIiJyDWYHbInMVQlERESkCEwMiIhIUVonH9pbbFFcXIyJEydCr9dDpVJhx44dFudVKpXV8oc//EGKiY2NbXN+xowZFvVcvnwZiYmJ0Gg00Gg0SExMRH19vU1tZWJARESKYkYvhxRbXL16FWFhYVi1apXV87W1tRZl/fr1UKlUmDp1qkVccnKyRdzatWstzs+cORPl5eXIy8tDXl4eysvLkZiYaFNbOceAiIgUpUWo0GLn0xFtvT4hIQEJCQntntdqtRavd+7cibFjx+Lee++1ON6vX782sa3OnDmDvLw8HDx4EKNGjQIArFu3DtHR0aioqEBwcHCn2soeAyIiottkMpksSmNjo911Xrx4Ebt378acOXPanMvOzoavry+GDx+OtLQ0XLlyRTpXVlYGjUYjJQUAEBUVBY1Gg9LS0k6/P3sMiIhIUVocsCqh5d+rEgICAiyOv/POO0hPT7er7o0bN8LLywtTpkyxOP78888jKCgIWq0WJ0+exJIlS/DNN98gPz8fAGAwGODn59emPj8/PxgMhk6/PxMDIiJSFLPoBbOdOx+a/73zYU1NDby9vaXjarXarnoBYP369Xj++efh4eFhcTw5OVn6OSQkBEOHDkVkZCSOHTuGESNGALgxifFWQgirx9vDoQQiIqLb5O3tbVHsTQy+/PJLVFRU4Ne//nWHsSNGjICbmxsqKysB3JincPHixTZxP/74I/z9/TvdBiYGRESkKK1DCfYWZ8jMzERERATCwsI6jD116hSam5uh0+kAANHR0TAajfj666+lmEOHDsFoNCImJqbTbeBQAhERKYoZtq8qsFaHLRoaGnD27FnpdVVVFcrLy+Hj44PBgwcDuDGR8c9//jOWLVvW5vp//vOfyM7OxpNPPglfX1+cPn0aCxcuRHh4OEaPHg0AuP/++zF+/HgkJydLyxjnzp2LCRMmdHpFAuCEHoOMjAyMHDkSXl5e8PPzw+TJk1FRUSF7TWFhodWNHb799ltHN4+IiKjLHTlyBOHh4QgPDwcApKamIjw8HG+//bYUk5OTAyEEnnvuuTbXu7u7429/+xvi4+MRHByMV199FXFxcSgoKEDv3r2luOzsbISGhiIuLg5xcXF48MEHsXnzZpvaqhLCsc+OHD9+PGbMmIGRI0fi+vXrePPNN3HixAmcPn0anp6eVq8pLCzE2LFjUVFRYTGJY+DAgRYfWI7JZIJGo0GsajL6qNwc8lmIiKhrXBfNKBQ7YDQaLb4HHKn1e2LNsZHo+wv7Osz/1XAd/zXisFPb210cPpSQl5dn8XrDhg3w8/PD0aNH8eijj8pe6+fnh7vuusvRTSIiIpLczpbG1upwVU7/ZEajEQDg4+PTYWx4eDh0Oh3GjRuHAwcOyMY2Nja22ViCiIiI7OPUxEAIgdTUVDz88MMICQlpN06n0+HTTz/Ftm3bsH37dgQHB2PcuHEoLi5u95qMjAzpIREajabNJhNERETWmKFySHFVDp9jcLP58+dj9+7dKCkpwaBBg2y6duLEiVCpVNi1a5fV842NjRZbT5pMJgQEBHCOARHRHagr5xgsPxLjkDkGr0WWco6BLRYsWIBdu3ahuLjY5qQAuLG/85YtW9o9r1arHbLDFBERKYtjtkR23TkGDk8MhBBYsGABcnNzUVhYiKCgoNuq5/jx49KmDURERNQ1HJ4YzJ8/H5999hl27twJLy8v6cENGo0Gffv2BQAsWbIEFy5cwKZNmwAAK1aswD333IPhw4ejqakJW7ZswbZt27Bt2zZHN4+IiBTOLFQw27vBkZ3X92QOTwzWrFkDAIiNjbU4vmHDBsyePRsAUFtbi+rqaulcU1MT0tLScOHCBfTt2xfDhw/H7t278eSTTzq6eUREpHBmBwwlmDmU0HmdmcuYlZVl8XrRokVYtGiRo5tCRERENuKzEoiISFEc89hl9hgQERG5hBao0GLnPgT2Xt+TuW7KQ0RERDZjjwERESkKhxLkMTEgIiJFaYH9QwEtjmlKj+S6KQ8RERHZjD0GRESkKBxKkMfEgIiIFKVF9EKLnV/s9l7fkzExICIiRREOeGyy4HJFIiIiUgL2GBARkaJwKEEeEwMiIlIUPl1RnuumPERERGQz9hgQEZGitDjgscv2Xt+TMTEgIiJF4VCCPNdNeYiIiMhm7DEgIiJFMaMXzHb+Xmzv9T0ZEwMiIlKUFqFCi51DAfZe35O5bspDRERENmOPARERKQonH8pjYkBERIoiHPB0RcGdD4mIiFxDC1RosfMhSPZe35O5bspDRERENmOPARERKYpZ2D9HwCwc1JgeiIkBEREpitkBcwzsvb4nc91PRkRERDZjYkBERIpihsohxRbFxcWYOHEi9Ho9VCoVduzYYXF+9uzZUKlUFiUqKsoiprGxEQsWLICvry88PT0xadIknD9/3iLm8uXLSExMhEajgUajQWJiIurr621qKxMDIiJSlNadD+0ttrh69SrCwsKwatWqdmPGjx+P2tpaqezZs8fifEpKCnJzc5GTk4OSkhI0NDRgwoQJaGlpkWJmzpyJ8vJy5OXlIS8vD+Xl5UhMTLSprZxjQERE5GQJCQlISEiQjVGr1dBqtVbPGY1GZGZmYvPmzXj88ccBAFu2bEFAQAAKCgoQHx+PM2fOIC8vDwcPHsSoUaMAAOvWrUN0dDQqKioQHBzcqbayx4CIiBSldfKhvQUATCaTRWlsbLztdhUWFsLPzw/Dhg1DcnIy6urqpHNHjx5Fc3Mz4uLipGN6vR4hISEoLS0FAJSVlUGj0UhJAQBERUVBo9FIMZ3BxICIiBTFDJW0LfJtl3/PMQgICJDG8zUaDTIyMm6rTQkJCcjOzsb+/fuxbNkyHD58GI899piUaBgMBri7u6N///4W1/n7+8NgMEgxfn5+ber28/OTYjqDQwlERES3qaamBt7e3tJrtVp9W/VMnz5d+jkkJASRkZEIDAzE7t27MWXKlHavE0JApfrPfIebf24vpiPsMSAiIkURDliRIP7dY+Dt7W1RbjcxuJVOp0NgYCAqKysBAFqtFk1NTbh8+bJFXF1dHfz9/aWYixcvtqnrxx9/lGI6g4kBEREpit3DCA54OmNHLl26hJqaGuh0OgBAREQE3NzckJ+fL8XU1tbi5MmTiImJAQBER0fDaDTi66+/lmIOHToEo9EoxXQGhxKIiEhRumPnw4aGBpw9e1Z6XVVVhfLycvj4+MDHxwfp6emYOnUqdDodvv/+e7zxxhvw9fXFM888AwDQaDSYM2cOFi5ciAEDBsDHxwdpaWkIDQ2VVincf//9GD9+PJKTk7F27VoAwNy5czFhwoROr0gAmBgQERE53ZEjRzB27FjpdWpqKgAgKSkJa9aswYkTJ7Bp0ybU19dDp9Nh7Nix2Lp1K7y8vKRrli9fjj59+mDatGn417/+hXHjxiErKwu9e/eWYrKzs/Hqq69KqxcmTZoku3eCNSohhEMfBZGeno7f/e53FsdunjVpTVFREVJTU3Hq1Cno9XosWrQI8+bNs+l9TSYTNBoNYlWT0UfldlttJyKi7nFdNKNQ7IDRaLSYzOdIrd8TT+97EW6e7nbV1Xy1CTvj1ju1vd3FKT0Gw4cPR0FBgfT65mzmVlVVVXjyySeRnJyMLVu24KuvvsLLL7+MgQMHYurUqc5oHhERKdjtbGlsrQ5X5ZTEoE+fPu3u3nSrTz75BIMHD8aKFSsA3BgjOXLkCP74xz/KJgaNjY0WG0mYTCa72kxEREROWpVQWVkJvV6PoKAgzJgxA9999127sWVlZRY7OQFAfHw8jhw5gubm5navy8jIsNhUIiAgwGHtJyIi13UnrEroTg5PDEaNGoVNmzZh7969WLduHQwGA2JiYnDp0iWr8QaDoc36Sn9/f1y/fh0//fRTu++zZMkSGI1GqdTU1Dj0cxARkWtiYiDP4UMJNz8kIjQ0FNHR0RgyZAg2btwozcK81a07MrXOh5TbqUmtVjtsIwkiIiK6wenLFT09PREaGirt3nQrrVbbZsVCXV0d+vTpgwEDBji7eUREpDCO+I3flXsMnL7zYWNjI86cOSPt3nSr6Ohoi52cAGDfvn2IjIyEmxuXHRIRkWNxKEGewxODtLQ0FBUVoaqqCocOHcKzzz4Lk8mEpKQkADfmBsyaNUuKnzdvHs6dO4fU1FScOXMG69evR2ZmJtLS0hzdNCIiIuqAw4cSzp8/j+eeew4//fQTBg4ciKioKBw8eBCBgYEAbuztXF1dLcUHBQVhz549eO211/Dxxx9Dr9fjo48+4h4GRETkFAL270Pg0J0BexiHJwY5OTmy57OystocGzNmDI4dO+bophAREbXBOQby+KwEIiJSFCYG8vjYZSIiIpKwx4CIiBSFPQbymBgQEZGiMDGQx6EEIiIikrDHgIiIFEUIFYSdv/Hbe31PxsSAiIgUxQyV3fsY2Ht9T8ahBCIiIpKwx4CIiBSFkw/lMTEgIiJF4RwDeRxKICIiIgl7DIiISFE4lCCPiQERESkKhxLkMTEgIiJFEQ7oMXDlxIBzDIiIiEjCHgMiIlIUAUAI++twVUwMiIhIUcxQQcWdD9vFoQQiIiKSsMeAiIgUhasS5DExICIiRTELFVTcx6BdHEogIiIiCXsMiIhIUYRwwKoEF16WwMSAiIgUhXMM5HEogYiIiCTsMSAiIkVhj4E8JgZERKQoXJUgj0MJRESkKK2TD+0ttiguLsbEiROh1+uhUqmwY8cO6VxzczNef/11hIaGwtPTE3q9HrNmzcIPP/xgUUdsbCxUKpVFmTFjhkXM5cuXkZiYCI1GA41Gg8TERNTX19vUViYGRERETnb16lWEhYVh1apVbc79/PPPOHbsGN566y0cO3YM27dvxz/+8Q9MmjSpTWxycjJqa2ulsnbtWovzM2fORHl5OfLy8pCXl4fy8nIkJiba1FYOJRARkaLc+I3f3jkGN/7fZDJZHFer1VCr1W3iExISkJCQYLUujUaD/Px8i2MrV67EQw89hOrqagwePFg63q9fP2i1Wqv1nDlzBnl5eTh48CBGjRoFAFi3bh2io6NRUVGB4ODgTn029hgQEZGitE4+tLcAQEBAgNRtr9FokJGR4ZA2Go1GqFQq3HXXXRbHs7Oz4evri+HDhyMtLQ1XrlyRzpWVlUGj0UhJAQBERUVBo9GgtLS00+/NHgMiIqLbVFNTA29vb+m1td4CW127dg2LFy/GzJkzLep+/vnnERQUBK1Wi5MnT2LJkiX45ptvpN4Gg8EAPz+/NvX5+fnBYDB0+v2ZGBARkaKIfxd76wAAb29viy9vezU3N2PGjBkwm81YvXq1xbnk5GTp55CQEAwdOhSRkZE4duwYRowYAQBQqdoOkQghrB5vD4cSiIhIURw5lOBIzc3NmDZtGqqqqpCfn99hwjFixAi4ubmhsrISAKDVanHx4sU2cT/++CP8/f073Q4mBkRERN2sNSmorKxEQUEBBgwY0OE1p06dQnNzM3Q6HQAgOjoaRqMRX3/9tRRz6NAhGI1GxMTEdLotHEogIiJlceRYQic1NDTg7Nmz0uuqqiqUl5fDx8cHer0ezz77LI4dO4a//OUvaGlpkeYE+Pj4wN3dHf/85z+RnZ2NJ598Er6+vjh9+jQWLlyI8PBwjB49GgBw//33Y/z48UhOTpaWMc6dOxcTJkzo9IoEwAk9Bvfcc0+bDRhUKhXmz59vNb6wsNBq/LfffuvophEREQGOGEawcSjhyJEjCA8PR3h4OAAgNTUV4eHhePvtt3H+/Hns2rUL58+fxy9/+UvodDqptK4mcHd3x9/+9jfEx8cjODgYr776KuLi4lBQUIDevXtL75OdnY3Q0FDExcUhLi4ODz74IDZv3mxTWx3eY3D48GG0tLRIr0+ePIknnngCv/rVr2Svq6iosBhPGThwoKObRkRE1C2PXY6NjYWQuUjuHHBjWWRRUVGH7+Pj44MtW7bY1rhbODwxuPULfenSpRgyZAjGjBkje52fn1+b9ZpERETUtZw6+bCpqQlbtmzBiy++2OFSifDwcOh0OowbNw4HDhzosO7GxkaYTCaLQkRE1JGeuiqhp3BqYrBjxw7U19dj9uzZ7cbodDp8+umn2LZtG7Zv347g4GCMGzcOxcXFsnVnZGRY7DYVEBDg4NYTEZFLap0jYG9xUSrR0cCGHeLj4+Hu7o4vvvjCpusmTpwIlUqFXbt2tRvT2NiIxsZG6bXJZEJAQABiVZPRR+V2220mIqKud100o1DsgNFodOiGQTczmUzQaDS4J/Mt9OrnYVdd5p+v4fs57zq1vd3FacsVz507h4KCAmzfvt3ma6OiojqcPNHegyqIiIjkdMfkwzuJ0xKDDRs2wM/PD0899ZTN1x4/flzasIGIiMihumEfgzuJUxIDs9mMDRs2ICkpCX36WL7FkiVLcOHCBWzatAkAsGLFCtxzzz0YPny4NFlx27Zt2LZtmzOaRkRERDKckhgUFBSguroaL774YptztbW1qK6ull43NTUhLS0NFy5cQN++fTF8+HDs3r0bTz75pDOaRkRECueIVQWuvCrBKYlBXFxcu5s1ZGVlWbxetGgRFi1a5IxmEBERWefCQwH24kOUiIiISMKHKBERkaJwKEEeEwMiIlIWrkqQxcSAiIgURvXvYm8drolzDIiIiEjCHgMiIlIWDiXIYmJARETKwsRAFocSiIiISMIeAyIiUhZHPDaZyxWJiIhcA5+uKI9DCURERCRhjwERESkLJx/KYmJARETKwjkGsjiUQERERBL2GBARkaKoxI1ibx2uiokBEREpC+cYyGJiQEREysI5BrI4x4CIiIgk7DEgIiJl4VCCLCYGRESkLEwMZHEogYiIiCTsMSAiImVhj4EsJgZERKQsXJUgi0MJREREJGGPARERKQp3PpTHxICIiJSFcwxkcSiBiIjIyYqLizFx4kTo9XqoVCrs2LHD4rwQAunp6dDr9ejbty9iY2Nx6tQpi5jGxkYsWLAAvr6+8PT0xKRJk3D+/HmLmMuXLyMxMREajQYajQaJiYmor6+3qa1MDIiIiJzs6tWrCAsLw6pVq6ye/+CDD/Dhhx9i1apVOHz4MLRaLZ544glcuXJFiklJSUFubi5ycnJQUlKChoYGTJgwAS0tLVLMzJkzUV5ejry8POTl5aG8vByJiYk2tZVDCUREpCgqOGCOwb//32QyWRxXq9VQq9Vt4hMSEpCQkGC1LiEEVqxYgTfffBNTpkwBAGzcuBH+/v747LPP8NJLL8FoNCIzMxObN2/G448/DgDYsmULAgICUFBQgPj4eJw5cwZ5eXk4ePAgRo0aBQBYt24doqOjUVFRgeDg4E59NvYYEBGRsrQuV7S3AAgICJC67TUaDTIyMmxuTlVVFQwGA+Li4qRjarUaY8aMQWlpKQDg6NGjaG5utojR6/UICQmRYsrKyqDRaKSkAACioqKg0WikmM5gjwEREdFtqqmpgbe3t/TaWm9BRwwGAwDA39/f4ri/vz/OnTsnxbi7u6N///5tYlqvNxgM8PPza1O/n5+fFNMZTAyIiEhZHLgqwdvb2yIxsIdKZblpkhCizbE2zbglxlp8Z+q5GYcSiIhIWYSDioNotVoAaPNbfV1dndSLoNVq0dTUhMuXL8vGXLx4sU39P/74Y5veCDlMDIiIiLpRUFAQtFot8vPzpWNNTU0oKipCTEwMACAiIgJubm4WMbW1tTh58qQUEx0dDaPRiK+//lqKOXToEIxGoxTTGRxKICIiRemOnQ8bGhpw9uxZ6XVVVRXKy8vh4+ODwYMHIyUlBe+99x6GDh2KoUOH4r333kO/fv0wc+ZMAIBGo8GcOXOwcOFCDBgwAD4+PkhLS0NoaKi0SuH+++/H+PHjkZycjLVr1wIA5s6diwkTJnR6RQJwGz0GjtikwZpt27bhgQcegFqtxgMPPIDc3Fxbm0ZERNSxbhhKOHLkCMLDwxEeHg4ASE1NRXh4ON5++20AwKJFi5CSkoKXX34ZkZGRuHDhAvbt2wcvLy+pjuXLl2Py5MmYNm0aRo8ejX79+uGLL75A7969pZjs7GyEhoYiLi4OcXFxePDBB7F582ab2qoSQtj08f7617/iq6++wogRIzB16lTk5uZi8uTJ0vn3338f//M//4OsrCwMGzYMv//971FcXIyKigqLD3izsrIyPPLII3j33XfxzDPPIDc3F2+//TZKSkosll3IMZlM0Gg0iFVNRh+Vmy0fiYiIutl10YxCsQNGo9Fhk/lu1fo9cc/v/we9PDzsqst87Rq+/+2bTm1vd7E5MbC4WKWySAyEENDr9UhJScHrr78O4MYWjv7+/nj//ffx0ksvWa1n+vTpMJlM+Otf/yodGz9+PPr374/PP/+8U21hYkBEdOfq0sTgXQclBm+5ZmLg0MmHndmkwZqysjKLawAgPj5e9prGxkaYTCaLQkRE1JHWOQb2Flfl0MRAbpMGuc0VDAaDzddkZGRY7DYVEBBgR8uJiIgIcNJyxdvZpMHWa5YsWQKj0SiVmpqa228wEREphwO3RHZFDl2uePMmDTqdTjp+8wYM7V0nt7GDNe09qIKIiEiWA3c+dEUO7THozCYN1kRHR1tcAwD79u2zaUMGIiKizuAcA3k29xjYu0kDAMyaNQt333239BSq3/zmN3j00Ufx/vvv4+mnn8bOnTtRUFCAkpISB3xEIiIi6iybE4MjR45g7Nix0uvU1FQAQFJSErKysrBo0SL861//wssvv4zLly9j1KhRbTZpqK6uRq9e/+msiImJQU5ODn7729/irbfewpAhQ7B169ZO72FARETUaRxKkGXXPgY9CfcxICK6c3XlPgb3vvUeetu5j0HLtWv47t03uI8BERERuTbXe4iSqteNYg9hdkxbiIgczd5/33qsXl3XPc+hBFmulxgQERHJYWIgy1VTTyIiIroN7DEgIiJFccQ+BK68jwF7DIiIiEjCxICIiIgkHEogIiJl4eRDWUwMiIhIUTjHQB4TAyIiUh4X/mK3F+cYEBERkYQ9BkREpCycYyCLiQERESkK5xjI41ACERERSdhjQEREysKhBFlMDIiISFE4lCCPQwlEREQkYY8BEREpC4cSZDExICIiZWFiIItDCURERCRhj4E1KuZLBECYu7sFdwb+90J3GE4+lMfEgIiIlIVDCbKYGBARkbIwMZDFPkAiIiKSsMeAiIgUhXMM5DExICIiZeFQgiwOJRARETnRPffcA5VK1abMnz8fADB79uw256KioizqaGxsxIIFC+Dr6wtPT09MmjQJ58+fd0p7mRgQEZGitA4l2Fs66/Dhw6itrZVKfn4+AOBXv/qVFDN+/HiLmD179ljUkZKSgtzcXOTk5KCkpAQNDQ2YMGECWlpaHHJPbsahBCIiUpYuHkoYOHCgxeulS5diyJAhGDNmjHRMrVZDq9Vavd5oNCIzMxObN2/G448/DgDYsmULAgICUFBQgPj4eNvbL4M9BkRERLfJZDJZlMbGRtn4pqYmbNmyBS+++CJUKpV0vLCwEH5+fhg2bBiSk5NRV1cnnTt69Ciam5sRFxcnHdPr9QgJCUFpaanDPxMTAyIiUhbhoAIgICAAGo1GKhkZGbJvvWPHDtTX12P27NnSsYSEBGRnZ2P//v1YtmwZDh8+jMcee0xKMgwGA9zd3dG/f3+Luvz9/WEwGOy5E1ZxKIGIiBRF9e9ibx0AUFNTA29vb+m4Wq2WvS4zMxMJCQnQ6/XSsenTp0s/h4SEIDIyEoGBgdi9ezemTJnSbl1CCIteB0dhYkBERHSbvL29LRIDOefOnUNBQQG2b98uG6fT6RAYGIjKykoAgFarRVNTEy5fvmzRa1BXV4eYmJjbb3w7OJRARETK4sChBFts2LABfn5+eOqpp2TjLl26hJqaGuh0OgBAREQE3NzcpNUMAFBbW4uTJ086JTFgjwERESlKd+x8aDabsWHDBiQlJaFPn/989TY0NCA9PR1Tp06FTqfD999/jzfeeAO+vr545plnAAAajQZz5szBwoULMWDAAPj4+CAtLQ2hoaHSKgVHsrnHoLi4GBMnToRer4dKpcKOHTukc83NzXj99dcRGhoKT09P6PV6zJo1Cz/88INsnVlZWVY3f7h27ZrNH4iIiEhWN/QYFBQUoLq6Gi+++KLF8d69e+PEiRN4+umnMWzYMCQlJWHYsGEoKyuDl5eXFLd8+XJMnjwZ06ZNw+jRo9GvXz988cUX6N27923cAHk29xhcvXoVYWFheOGFFzB16lSLcz///DOOHTuGt956C2FhYbh8+TJSUlIwadIkHDlyRLZeb29vVFRUWBzz8PCwtXlEREQ9TlxcHIRom0307dsXe/fu7fB6Dw8PrFy5EitXrnRG8yzYnBgkJCQgISHB6jmNRmMxBgIAK1euxEMPPYTq6moMHjy43XpVKlW7mzsQERE5lAs/68BeTp9jYDQaoVKpcNddd8nGNTQ0IDAwEC0tLfjlL3+Jd999F+Hh4e3GNzY2WmwkYTKZHNVkohtUnJtL5Ir4dEV5Tv2X79q1a1i8eDFmzpwpu5zjvvvuQ1ZWFnbt2oXPP/8cHh4eGD16tLRUw5qMjAyLTSUCAgKc8RGIiIgUxWmJQXNzM2bMmAGz2YzVq1fLxkZFReH//b//h7CwMDzyyCP43//9XwwbNkx2LGXJkiUwGo1SqampcfRHICIiV9RNyxXvFE4ZSmhubsa0adNQVVWF/fv3d3rzh1a9evXCyJEjZXsM1Gp1hztMERER3YpDCfIc3mPQmhRUVlaioKAAAwYMsLkOIQTKy8ulzR2IiIioa9jcY9DQ0ICzZ89Kr6uqqlBeXg4fHx/o9Xo8++yzOHbsGP7yl7+gpaVFesCDj48P3N3dAQCzZs3C3XffLT1s4ne/+x2ioqIwdOhQmEwmfPTRRygvL8fHH3/siM9IRET0H1382OU7jc2JwZEjRzB27FjpdWpqKgAgKSkJ6enp2LVrFwDgl7/8pcV1Bw4cQGxsLACguroavXr9p7Oivr4ec+fOhcFggEajQXh4OIqLi/HQQw/Z2jwiIiJZHEqQZ3NiEBsba3WThlZy51oVFhZavF6+fDmWL19ua1OIiIjIwfisBCIiUhYOJchiYkBERMrCxEAWEwMiIlIUzjGQxz1fiYiISMIeAyIiUhYOJchiYkBERIqiEgKqTqyg66gOV8WhBCIiIpKwx4CIiJSFQwmymBgQEZGicFWCPA4lEBERkYQ9BkREpCwcSpDFxICIiBSFQwnyOJRAREREEvYYEBGRsnAoQRYTAyIiUhQOJchjYkBERMrCHgNZnGNAREREEvYYEBGR4rjyUIC9mBgQEZGyCHGj2FuHi+JQAhEREUnYY0BERIrCVQnymBgQEZGycFWCLA4lEBERkYQ9BkREpCgq841ibx2uiokBEREpC4cSZHEogYiIiCRMDIiISFFaVyXYWzorPT0dKpXKomi1Wum8EALp6enQ6/Xo27cvYmNjcerUKYs6GhsbsWDBAvj6+sLT0xOTJk3C+fPnHXVLLDAxICIiZWnd4MjeYoPhw4ejtrZWKidOnJDOffDBB/jwww+xatUqHD58GFqtFk888QSuXLkixaSkpCA3Nxc5OTkoKSlBQ0MDJkyYgJaWFofdllacY0BERIrSHfsY9OnTx6KXoJUQAitWrMCbb76JKVOmAAA2btwIf39/fPbZZ3jppZdgNBqRmZmJzZs34/HHHwcAbNmyBQEBASgoKEB8fLx9H+YW7DEgIiK6TSaTyaI0NjZajausrIRer0dQUBBmzJiB7777DgBQVVUFg8GAuLg4KVatVmPMmDEoLS0FABw9ehTNzc0WMXq9HiEhIVKMIzExICIiZREOKgACAgKg0WikkpGR0ebtRo0ahU2bNmHv3r1Yt24dDAYDYmJicOnSJRgMBgCAv7+/xTX+/v7SOYPBAHd3d/Tv37/dGEfiUAIRESmKI4cSampq4O3tLR1Xq9VtYhMSEqSfQ0NDER0djSFDhmDjxo2Iioq6UZ9KZXGNEKLNsVt1JuZ2sMeAiIjoNnl7e1sUa4nBrTw9PREaGorKykpp3sGtv/nX1dVJvQharRZNTU24fPlyuzGOxMSAiIiUpRtWJdyssbERZ86cgU6nQ1BQELRaLfLz86XzTU1NKCoqQkxMDAAgIiICbm5uFjG1tbU4efKkFONIHEogIiJF6epVCWlpaZg4cSIGDx6Muro6/P73v4fJZEJSUhJUKhVSUlLw3nvvYejQoRg6dCjee+899OvXDzNnzgQAaDQazJkzBwsXLsSAAQPg4+ODtLQ0hIaGSqsUHMnmHoPi4mJMnDgRer0eKpUKO3bssDg/e/bsNhs5tI6hyNm2bRseeOABqNVqPPDAA8jNzbW1aURERD3O+fPn8dxzzyE4OBhTpkyBu7s7Dh48iMDAQADAokWLkJKSgpdffhmRkZG4cOEC9u3bBy8vL6mO5cuXY/LkyZg2bRpGjx6Nfv364YsvvkDv3r0d3l6bewyuXr2KsLAwvPDCC5g6darVmPHjx2PDhg3Sa3d3d9k6y8rKMH36dLz77rt45plnkJubi2nTpqGkpASjRo2ytYlERETt6+JnJeTk5MieV6lUSE9PR3p6ersxHh4eWLlyJVauXNn5N75NNicGCQkJFjMsrVGr1VY3cmjPihUr8MQTT2DJkiUAgCVLlqCoqAgrVqzA559/bmsTiYiI2tUdGxzdSZwy+bCwsBB+fn4YNmwYkpOTUVdXJxtfVlZmsXEDAMTHx8tu3NDY2NhmYwkiIiKyj8MTg4SEBGRnZ2P//v1YtmwZDh8+jMcee6zd3aCAG8s05DZ3sCYjI8NiU4mAgACHfQYiInJhZuGY4qIcviph+vTp0s8hISGIjIxEYGAgdu/eLe0DbY2tmzssWbIEqamp0muTycTkgIiIOtbFcwzuNE5frqjT6RAYGIjKysp2Y7RarezmDtao1epObSRBRER0MxUcMMfAIS3pmZy+wdGlS5dQU1MDnU7Xbkx0dLTFxg0AsG/fPqds3EBERETts7nHoKGhAWfPnpVeV1VVoby8HD4+PvDx8UF6ejqmTp0KnU6H77//Hm+88QZ8fX3xzDPPSNfMmjULd999t/Swid/85jd49NFH8f777+Ppp5/Gzp07UVBQgJKSEgd8RCIiopvYuXOhVIeLsjkxOHLkCMaOHSu9bh3nT0pKwpo1a3DixAls2rQJ9fX10Ol0GDt2LLZu3WqxUUN1dTV69fpPZ0VMTAxycnLw29/+Fm+99RaGDBmCrVu3cg8DIiJyOC5XlGdzYhAbGwshkynt3bu3wzoKCwvbHHv22Wfx7LPP2tocIiIiciA+K4GIiJSFqxJkMTEgIiJFUQkBlZ1zBOy9vifjY5eJiIhIwh4DIiJSFvO/i711uCgmBkREpCgcSpDHoQQiIiKSsMeAiIiUhasSZDExICIiZeHOh7KYGBARkaJw50N5nGNAREREEvYYEBGRsnAoQRYTAyIiUhSV+Uaxtw5XxaEEIiIikrDHgIiIlIVDCbKYGBARkbJwHwNZHEogIiIiCXsMiIhIUfisBHlMDIiISFk4x0AWhxKIiIhIwh4DIiJSFgHA3n0IXLfDgIkBEREpC+cYyGNiQEREyiLggDkGDmlJj8Q5BkRERCRhjwERESkLVyXIYmJARETKYgagckAdLopDCURERCRhYkBERIrSuirB3tJZGRkZGDlyJLy8vODn54fJkyejoqLCImb27NlQqVQWJSoqyiKmsbERCxYsgK+vLzw9PTFp0iScP3/eIffkZkwMiIhIWVrnGNhbOqmoqAjz58/HwYMHkZ+fj+vXryMuLg5Xr161iBs/fjxqa2ulsmfPHovzKSkpyM3NRU5ODkpKStDQ0IAJEyagpaXFIbelFecYEBEROVFeXp7F6w0bNsDPzw9Hjx7Fo48+Kh1Xq9XQarVW6zAajcjMzMTmzZvx+OOPAwC2bNmCgIAAFBQUID4+3mHtZY8BEREpiwN7DEwmk0VpbGzs8O2NRiMAwMfHx+J4YWEh/Pz8MGzYMCQnJ6Ourk46d/ToUTQ3NyMuLk46ptfrERISgtLSUkfcFQkTAyIiUhYHJgYBAQHQaDRSycjI6OCtBVJTU/Hwww8jJCREOp6QkIDs7Gzs378fy5Ytw+HDh/HYY49JiYbBYIC7uzv69+9vUZ+/vz8MBoNDbw+HEoiIiG5TTU0NvL29pddqtVo2/pVXXsHf//53lJSUWByfPn269HNISAgiIyMRGBiI3bt3Y8qUKe3WJ4SASmXv2ktL7DEgIiJlMTuoAPD29rYoconBggULsGvXLhw4cACDBg2SbaJOp0NgYCAqKysBAFqtFk1NTbh8+bJFXF1dHfz9/W36+B1hYkBERIrS1csVhRB45ZVXsH37duzfvx9BQUEdXnPp0iXU1NRAp9MBACIiIuDm5ob8/Hwppra2FidPnkRMTIztN0EGhxKIiEhZunhL5Pnz5+Ozzz7Dzp074eXlJc0J0Gg06Nu3LxoaGpCeno6pU6dCp9Ph+++/xxtvvAFfX18888wzUuycOXOwcOFCDBgwAD4+PkhLS0NoaKi0SsFRmBgQERE50Zo1awAAsbGxFsc3bNiA2bNno3fv3jhx4gQ2bdqE+vp66HQ6jB07Flu3boWXl5cUv3z5cvTp0wfTpk3Dv/71L4wbNw5ZWVno3bu3Q9tr81BCcXExJk6cCL1eD5VKhR07dlicv3Xnptbyhz/8od06s7KyrF5z7do1mz8QERGRLLNwTOkkIYTVMnv2bABA3759sXfvXtTV1aGpqQnnzp1DVlYWAgICLOrx8PDAypUrcenSJfz888/44osv2sQ4gs09BlevXkVYWBheeOEFTJ06tc352tpai9d//etfMWfOHKuxN/P29m6zRaSHh4etzSMiIpLHpyvKsjkxSEhIQEJCQrvnb921aefOnRg7dizuvfde2XpVKlW7Oz5Z09jYaLGRhMlk6vS1REREZJ1TVyVcvHgRu3fvxpw5czqMbWhoQGBgIAYNGoQJEybg+PHjsvEZGRkWm0o4ozuFiIhckSM2N3LdHgOnJgYbN26El5eX7OYMAHDfffchKysLu3btwueffw4PDw+MHj1aWr9pzZIlS2A0GqVSU1Pj6OYTEZEr6uKHKN1pnLoqYf369Xj++ec7nCsQFRVl8XjJ0aNHY8SIEVi5ciU++ugjq9eo1eoOd5giIiIi2zgtMfjyyy9RUVGBrVu32nxtr169MHLkSNkeAyIiottidsBQgA2rEu40ThtKyMzMREREBMLCwmy+VgiB8vJyaccnIiIihxFmxxQXZXOPQUNDA86ePSu9rqqqQnl5OXx8fDB48GAAN1YI/PnPf8ayZcus1jFr1izcfffd0lOofve73yEqKgpDhw6FyWTCRx99hPLycnz88ce385mIiIjoNtmcGBw5cgRjx46VXqempgIAkpKSkJWVBQDIycmBEALPPfec1Tqqq6vRq9d/Oivq6+sxd+5cGAwGaDQahIeHo7i4GA899JCtzSMiIpLHfQxkqYRwjU9nMpmg0WgQ22sK+qjcurs5RERkg+uiGYXm7TAajRaPMXak1u+Jx++ehz697Ju8ft3ciIILnzi1vd2Fz0ogIiJlYY+BLD52mYiIiCTsMSAiImURcECPgUNa0iMxMSAiImXhUIIsDiUQERGRhD0GRESkLGYzADs3KDJzgyMiIiLXwKEEWRxKICIiIgl7DIiISFnYYyCLiQERESkLn64oi0MJREREJGGPARERKYoQZgg7H5ts7/U9GRMDIiJSFiHsHwrgHAMiIiIXIRwwx8CFEwPOMSAiIiIJewyIiEhZzGZAZeccAc4xICIichEcSpDFoQQiIiKSsMeAiIgURZjNEHYOJXC5IhERkavgUIIsDiUQERGRhD0GRESkLGYBqNhj0B4mBkREpCxCALB3uaLrJgYcSiAiIiIJewyIiEhRhFlA2DmUIFy4x4CJARERKYsww/6hBNddrsihBCIiUhRhFg4ptlq9ejWCgoLg4eGBiIgIfPnll074dPZjYkBERORkW7duRUpKCt58800cP34cjzzyCBISElBdXd3dTWvDZYYSWsd7rovmbm4JERHZqvXf7q4Yu78uGu0eCriOG+01mUwWx9VqNdRqdZv4Dz/8EHPmzMGvf/1rAMCKFSuwd+9erFmzBhkZGXa1xdFcJjG4cuUKAKBEfGH3hlZERNQ9rly5Ao1G45S63d3dodVqUWLY45D6fvGLXyAgIMDi2DvvvIP09HSLY01NTTh69CgWL15scTwuLg6lpaUOaYsjuUxioNfrUVNTAy8vL6hUKqsxJpMJAQEBqKmpgbe3dxe38Pax3V3vTm0729212G7HEULgypUr0Ov1TnsPDw8PVFVVoampySH1CSHafN9Y6y346aef0NLSAn9/f4vj/v7+MBgMDmmLI7lMYtCrVy8MGjSoU7He3t495j8GW7DdXe9ObTvb3bXYbsdwVk/BzTw8PODh4eH097Hm1iTCWmLRE3DyIRERkRP5+vqid+/ebXoH6urq2vQi9ARMDIiIiJzI3d0dERERyM/Ptzien5+PmJiYbmpV+1xmKKEz1Go13nnnHatjQD0Z29317tS2s91di+2mzkpNTUViYiIiIyMRHR2NTz/9FNXV1Zg3b153N60NlXDlfR2JiIh6iNWrV+ODDz5AbW0tQkJCsHz5cjz66KPd3aw2mBgQERGRhHMMiIiISMLEgIiIiCRMDIiIiEjCxICIiIgkLpcY2PpYy6KiIkRERMDDwwP33nsvPvnkky5q6Q0ZGRkYOXIkvLy84Ofnh8mTJ6OiokL2msLCQqhUqjbl22+/7aJWA+np6W3eX6vVyl7T3fe61T333GP1/s2fP99qfHfd7+LiYkycOBF6vR4qlQo7duywOC+EQHp6OvR6Pfr27YvY2FicOnWqw3q3bduGBx54AGq1Gg888AByc3O7rN3Nzc14/fXXERoaCk9PT+j1esyaNQs//PCDbJ1ZWVlW/wyuXbvWJe0GgNmzZ7d5/6ioqA7r7c77DcDqfVOpVPjDH/7Qbp1dcb+p53KpxMDWx1pWVVXhySefxCOPPILjx4/jjTfewKuvvopt27Z1WZuLioowf/58HDx4EPn5+bh+/Tri4uJw9erVDq+tqKhAbW2tVIYOHdoFLf6P4cOHW7z/iRMn2o3tCfe61eHDhy3a3brpyK9+9SvZ67r6fl+9ehVhYWFYtWqV1fMffPABPvzwQ6xatQqHDx+GVqvFE088IT1QzJqysjJMnz4diYmJ+Oabb5CYmIhp06bh0KFDXdLun3/+GceOHcNbb72FY8eOYfv27fjHP/6BSZMmdVivt7e3xf2vra116Na2Hd1vABg/frzF++/ZI/8wnu6+3wDa3LP169dDpVJh6tSpsvU6+35TDyZcyEMPPSTmzZtncey+++4Tixcvthq/aNEicd9991kce+mll0RUVJTT2tiRuro6AUAUFRW1G3PgwAEBQFy+fLnrGnaLd955R4SFhXU6vife61a/+c1vxJAhQ4TZbLZ6vifcbwAiNzdXem02m4VWqxVLly6Vjl27dk1oNBrxySeftFvPtGnTxPjx4y2OxcfHixkzZji8zUK0bbc1X3/9tQAgzp07127Mhg0bhEajcWzjZFhrd1JSknj66adtqqcn3u+nn35aPPbYY7IxXX2/qWdxmR6D1sdaxsXFWRyXe6xlWVlZm/j4+HgcOXIEzc3NTmurHKPRCADw8fHpMDY8PBw6nQ7jxo3DgQMHnN20NiorK6HX6xEUFIQZM2bgu+++aze2J95r4Mbfmy1btuDFF1/s8GEm3X2/b1ZVVQWDwWBxT9VqNcaMGSP7GNf2/hy689GvRqMRKpUKd911l2xcQ0MDAgMDMWjQIEyYMAHHjx/vmgbepLCwEH5+fhg2bBiSk5NRV1cnG9/T7vfFixexe/duzJkzp8PYnnC/qXu4TGJwO4+1NBgMVuOvX7+On376yWltbY8QAqmpqXj44YcREhLSbpxOp8Onn36Kbdu2Yfv27QgODsa4ceNQXFzcZW0dNWoUNm3ahL1792LdunUwGAyIiYnBpUuXrMb3tHvdaseOHaivr8fs2bPbjekJ9/tWrX+nbX2Ma3t/Dt316Ndr165h8eLFmDlzpuxT/u677z5kZWVh165d+Pzzz+Hh4YHRo0ejsrKyy9qakJCA7Oxs7N+/H8uWLcPhw4fx2GOPobGxsd1retr93rhxI7y8vDBlyhTZuJ5wv6n7uNyzEmx9rKW1eGvHu8Irr7yCv//97ygpKZGNCw4ORnBwsPQ6OjoaNTU1+OMf/9hl22smJCRIP4eGhiI6OhpDhgzBxo0bkZqaavWannSvW2VmZiIhIUH2GfA94X6353Ye49pTHv3a3NyMGTNmwGw2Y/Xq1bKxUVFRFhP9Ro8ejREjRmDlypX46KOPnN1UAMD06dOln0NCQhAZGYnAwEDs3r1b9ou2p9xvAFi/fj2ef/75DucK9IT7Td3HZXoMbuexllqt1mp8nz59MGDAAKe11ZoFCxZg165dOHDgAAYNGmTz9VFRUd2azXt6eiI0NLTdNvSke93q3LlzKCgowK9//Wubr+3u+926AsTWx7i29+fQ1Y9+bW5uxrRp01BVVYX8/HzZ3gJrevXqhZEjR3brn4FOp0NgYKBsG3rK/QaAL7/8EhUVFbf1970n3G/qOi6TGNzOYy2jo6PbxO/btw+RkZFwc3NzWltvJoTAK6+8gu3bt2P//v0ICgq6rXqOHz8OnU7n4NZ1XmNjI86cOdNuG3rCvb7Vhg0b4Ofnh6eeesrma7v7fgcFBUGr1Vrc06amJhQVFck+xrW9P4eufPRra1JQWVmJgoKC20oMhRAoLy/v1j+DS5cuoaamRrYNPeF+t8rMzERERATCwsJsvrYn3G/qQt0169EZcnJyhJubm8jMzBSnT58WKSkpwtPTU3z//fdCCCEWL14sEhMTpfjvvvtO9OvXT7z22mvi9OnTIjMzU7i5uYn/+7//67I2/9d//ZfQaDSisLBQ1NbWSuXnn3+WYm5t9/Lly0Vubq74xz/+IU6ePCkWL14sAIht27Z1WbsXLlwoCgsLxXfffScOHjwoJkyYILy8vHr0vb5ZS0uLGDx4sHj99dfbnOsp9/vKlSvi+PHj4vjx4wKA+PDDD8Xx48el2ftLly4VGo1GbN++XZw4cUI899xzQqfTCZPJJNWRmJhosSrnq6++Er179xZLly4VZ86cEUuXLhV9+vQRBw8e7JJ2Nzc3i0mTJolBgwaJ8vJyi7/zjY2N7bY7PT1d5OXliX/+85/i+PHj4oUXXhB9+vQRhw4d6pJ2X7lyRSxcuFCUlpaKqqoqceDAAREdHS3uvvvuHn2/WxmNRtGvXz+xZs0aq3V0x/2mnsulEgMhhPj4449FYGCgcHd3FyNGjLBY9peUlCTGjBljEV9YWCjCw8OFu7u7uOeee9r9D8dZAFgtGzZsaLfd77//vhgyZIjw8PAQ/fv3Fw8//LDYvXt3l7Z7+vTpQqfTCTc3N6HX68WUKVPEqVOn2m2zEN1/r2+2d+9eAUBUVFS0OddT7nfrMslbS1JSkhDixpLFd955R2i1WqFWq8Wjjz4qTpw4YVHHmDFjpPhWf/7zn0VwcLBwc3MT9913n8MTHLl2V1VVtft3/sCBA+22OyUlRQwePFi4u7uLgQMHiri4OFFaWtpl7f75559FXFycGDhwoHBzcxODBw8WSUlJorq62qKOnna/W61du1b07dtX1NfXW62jO+439Vx87DIRERFJXGaOAREREdmPiQERERFJmBgQERGRhIkBERERSZgYEBERkYSJAREREUmYGBAREZGEiQERERFJmBgQERGRhIkBERERSZgYEBERkeT/AygWtIhHw0BIAAAAAElFTkSuQmCC", "text/plain": [ - "
" + "
" ] }, - "metadata": { - "needs_background": "light" - }, + "metadata": {}, "output_type": "display_data" } ], @@ -369,13 +366,13 @@ }, { "cell_type": "code", - "execution_count": 37, + "execution_count": 127, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "324063b36c0147808a548646753cb24a", + "model_id": "77835379be7343d0aff52d275a86be97", "version_major": 2, "version_minor": 0 }, @@ -392,7 +389,7 @@ "" ] }, - "execution_count": 37, + "execution_count": 127, "metadata": {}, "output_type": "execute_result" }