diff --git a/Chapter 11.ipynb b/Chapter 11.ipynb new file mode 100644 index 0000000..7fc89b5 --- /dev/null +++ b/Chapter 11.ipynb @@ -0,0 +1,256 @@ +{ + "cells": [ + { + "cell_type": "raw", + "id": "9e3ff2ec-8487-4937-8659-4cc62805df44", + "metadata": {}, + "source": [ + "Diffusion avec mécanisme lacunaire" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "id": "ef9bd780-0b14-49a1-9687-ccbeb1504d99", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGzCAYAAABzfl4TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAdj0lEQVR4nO3df5BVdf348ddlkYV09wooArIQGWgKbCZKhFYmadRolqESTkRaU7P+inGmmCnBUpemNMvMX5U2FZE5kj9mCEkFplFGwJjBahQKcwMVbWQXtula7Pn88fm6n+8m6LnLe7n3wuMxc2e8l3v3vuY1d9znnHv23kKWZVkAACTQr9IDAAAHDmEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIyyAA8quXbtiwYIF8ZGPfCSGDBkShUIh7r777kqPBQcNYQEcUF555ZX4xje+EX/+85+jubm50uPAQad/pQcAeF1nZ2cceuih+/QzRowYES+88EIMHz481q1bFyeffHKi6YA8HLGAGrdy5cooFApxzz33xDXXXBNHH310NDQ0xKc+9alob2+PUqkUV155ZQwbNiwOO+ywmDt3bpRKpTf8nJ///Odx0kknxaBBg2LIkCFx4YUXRltbW4/7bNq0Kc4777wYPnx4DBw4MEaNGhUXXnhhtLe3R0TEc889t9e3HgqFQixcuLD7+sKFC6NQKMSf/vSn+PSnPx2DBw+OU089tax59qS+vj6GDx+ec3tAao5YwAGitbU1Bg0aFF/96ldj8+bNcfPNN8chhxwS/fr1i1dffTUWLlwYa9asibvvvjvGjh0bV199dfdjr7vuuvj6178e559/flxyySXx8ssvx8033xzvf//74w9/+EMcfvjh8dprr8VZZ50VpVIpLrvsshg+fHhs3bo1HnroodixY0cUi8VezT1z5swYN25cXH/99ZFlWe55gCqVATXtscceyyIimzBhQvbaa6913z5r1qysUChkM2bM6HH/qVOnZmPGjOm+/txzz2V1dXXZdddd1+N+GzduzPr37999+x/+8IcsIrJf//rXe51ly5YtWURkd9111xv+LSKyBQsWdF9fsGBBFhHZrFmzetwv7zx5rF27dq/zAH3DWyFwgPjMZz4ThxxySPf1KVOmRJZl8bnPfa7H/aZMmRJtbW3xn//8JyIi7rvvvujq6orzzz8/Xnnlle7L8OHDY9y4cfHYY49FRHQfkVi+fHn885//TDb3F7/4xR7X884DVCdvhcABYvTo0T2uvx4CTU1Nb7i9q6sr2tvbY+jQobFp06bIsizGjRu3x5/7eqyMHTs25s2bFzfeeGP84he/iNNOOy3OOeecuOiii3r9NsjrP/f/l3ceoDoJCzhA1NXVlXV79v/OZ+jq6opCoRDLli3b430PO+yw7v++4YYb4rOf/Wzcf//98fDDD8fll18era2tsWbNmhg1alQUCoU9Ptfu3bv3OvegQYN6XC9nHqD6CAs4yB1zzDGRZVmMHTs2xo8f/5b3nzhxYkycODG+9rWvxeOPPx7Tpk2L2267La699toYPHhwRETs2LGjx2P+9re/9dk8QHVxjgUc5D75yU9GXV1dXHPNNd1HMV6XZVn84x//iIiIjo6O7vMyXjdx4sTo169f95+vNjY2xhFHHBGrV6/ucb8f/vCHyecBqpMjFnCQO+aYY+Laa6+N+fPnx3PPPRfnnntuNDQ0xJYtW2Lp0qXxhS98Ia666qp49NFH49JLL42ZM2fG+PHj4z//+U/87Gc/i7q6ujjvvPO6f94ll1wSixYtiksuuSQmT54cq1evjmeffTb5PG/mBz/4QezYsSO2bdsWEREPPvhg/P3vf4+IiMsuu2yfzgkB3pywAOKrX/1qjB8/Pr773e/GNddcExH/e9LnmWeeGeecc05ERDQ3N8dZZ50VDz74YGzdujXe9ra3RXNzcyxbtize+973dv+sq6++Ol5++eW4995745577okZM2bEsmXLYtiwYUnneTPf+c53erz9ct9998V9990XEbHPJ5sCb66Q/fexRgCAXnKOBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACCZ/f45Fl1dXbFt27ZoaGjY6/cKAADVJcuy2LlzZ4wcOTL69dv7cYn9Hhbbtm17w7ctAgC1oa2tLUaNGrXXf9/vYdHQ0BAR/ztYY2Pj/n56AKAXOjo6oqmpqfv3+N7s97B4/e2PxsZGYQEANeatTmNw8iYAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACSz3z/HYl/19lPAsyztHLXArsrTm33ZVXkOxn3ZVX52lV8178oRCwAgGWEBACRTdlhs3bo1Lrroohg6dGgMGjQoJk6cGOvWreuL2QCAGlPWORavvvpqTJs2LU4//fRYtmxZHHnkkbFp06YYPHhwX80HANSQssLiW9/6VjQ1NcVdd93VfdvYsWPf9DGlUilKpVL39Y6OjjJHBABqRVlvhTzwwAMxefLkmDlzZgwbNixOPPHEuPPOO9/0Ma2trVEsFrsvTU1N+zQwAFC9ClmW/49PBg4cGBER8+bNi5kzZ8batWvjiiuuiNtuuy3mzJmzx8fs6YhFU1NTtLe39+pr06v5T2yqjV2Vx5+b5ue1lZ9d5WdX+VViVx0dHVEsFt/y93dZYTFgwICYPHlyPP744923XX755bF27dp44oknkg6214G98HKzq/IIi/y8tvKzq/zsKr9qDouy3goZMWJEHH/88T1ue9e73hXPP/9876YEAA4oZYXFtGnT4plnnulx27PPPhtjxoxJOhQAUJvKCosvf/nLsWbNmrj++utj8+bNsXjx4rjjjjuipaWlr+YDAGpIWWFx8sknx9KlS+OXv/xlTJgwIb75zW/GTTfdFLNnz+6r+QCAGlLWyZsp7OvJmwDA/tcnJ28CALwZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAk07/SA5SrUOjd47Is7Ry1wK7K05t92VV5DsZ92VV+dpVfNe/KEQsAIBlhAQAkU1ZYLFy4MAqFQo/Lcccd11ezAQA1puxzLE444YT43e9+938/oH/NnaYBAPSRsqugf//+MXz48Nz3L5VKUSqVuq93dHSU+5QAQI0o+xyLTZs2xciRI+Md73hHzJ49O55//vk3vX9ra2sUi8XuS1NTU6+HBQCqWyHL8v/xybJly2LXrl1x7LHHxgsvvBDXXHNNbN26NZ5++uloaGjY42P2dMSiqakp2tvbo7GxsfyBq/hPbKqNXZXHn5vm57WVn13lZ1f5VWJXHR0dUSwW3/L3d1lh8d927NgRY8aMiRtvvDEuvvjipIPtjRdefnZVHmGRn9dWfnaVn13lV81hsU9/bnr44YfH+PHjY/PmzfvyYwCAA8Q+hcWuXbviL3/5S4wYMSLVPABADSsrLK666qpYtWpVPPfcc/H444/HJz7xiairq4tZs2b11XwAQA0p689N//73v8esWbPiH//4Rxx55JFx6qmnxpo1a+LII4/sq/kAgBpSVlgsWbKkr+bI7WA8Sae37Ko89pWfXeVnV/nZVX7VvCvfFQIAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGT6V3qAchUKvXtclqWdoxbYVXl6sy+7Ks/BuC+7ys+u8qvmXTliAQAkIywAgGT2KSwWLVoUhUIhrrzyykTjAAC1rNdhsXbt2rj99ttj0qRJKecBAGpYr8Ji165dMXv27Ljzzjtj8ODBb3rfUqkUHR0dPS4AwIGpV2HR0tISH/vYx2L69Olved/W1tYoFovdl6ampt48JQBQA8oOiyVLlsRTTz0Vra2tue4/f/78aG9v7760tbWVPSQAUBvK+hyLtra2uOKKK2LFihUxcODAXI+pr6+P+vr6Xg0HANSWssJi/fr1sX379njPe97Tfdvu3btj9erV8YMf/CBKpVLU1dUlHxIAqA1lhcUZZ5wRGzdu7HHb3Llz47jjjouvfOUrogIADnJlhUVDQ0NMmDChx22HHnpoDB069A23AwAHH5+8CQAks89fQrZy5coEY+R3MH7ZTG/ZVXnsKz+7ys+u8rOr/Kp5V45YAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgmf6VHqBchULvHpdlaeeoBXZFX/Hays+u8rOr/Kp5V45YAADJCAsAIJmywuLWW2+NSZMmRWNjYzQ2NsbUqVNj2bJlfTUbAFBjygqLUaNGxaJFi2L9+vWxbt26+NCHPhQf//jH449//GNfzQcA1JBClu3bqRxDhgyJb3/723HxxRfv8d9LpVKUSqXu6x0dHdHU1BTt7e3R2NhY9vNV8wkr1cau6CteW/nZVX52lV8ldtXR0RHFYvEtf3/3+hyL3bt3x5IlS6KzszOmTp261/u1trZGsVjsvjQ1NfX2KQGAKlf2EYuNGzfG1KlT41//+lccdthhsXjx4vjoRz+61/s7YlE5dkVf8drKz67ys6v8qvmIRdmfY3HsscfGhg0bor29Pe69996YM2dOrFq1Ko4//vg93r++vj7q6+vLfRoAoAbt8zkW06dPj2OOOSZuv/32XPfPWzx7o2jzsyv6itdWfnaVn13lV81HLPb5cyy6urp6vNUBABy8ynorZP78+TFjxowYPXp07Ny5MxYvXhwrV66M5cuX99V8AEANKSsstm/fHp/5zGfihRdeiGKxGJMmTYrly5fHhz/84b6aDwCoIWWFxY9//OO+miO3g/G9tN6yK/qK11Z+dpWfXeVXzbvyXSEAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEimf6UHKFeh0LvHZVnaOWqBXZWnN/s6WHflxZWfVeVnV/lV864csQAAkhEWAEAyZYVFa2trnHzyydHQ0BDDhg2Lc889N5555pm+mg0AqDFlhcWqVauipaUl1qxZEytWrIh///vfceaZZ0ZnZ2dfzQcA1JBClvX+VI6XX345hg0bFqtWrYr3v//9e7xPqVSKUqnUfb2joyOampqivb09Ghsbyx+4ik9YqTZ2VR4nb5bBiys3q8rPrvKrxK46OjqiWCy+5e/vfTrHor29PSIihgwZstf7tLa2RrFY7L40NTXty1MCAFWs10csurq64pxzzokdO3bE73//+73ezxGLyrGr8jhiUQYvrtysKj+7yq+aj1j0+nMsWlpa4umnn37TqIiIqK+vj/r6+t4+DQBQQ3oVFpdeemk89NBDsXr16hg1alTqmQCAGlVWWGRZFpdddlksXbo0Vq5cGWPHju2ruQCAGlRWWLS0tMTixYvj/vvvj4aGhnjxxRcjIqJYLMagQYP6ZEAAoHaUdfJmYS9ni9x1113x2c9+NtfPyHvyx95nKPshEeHknnIcjLuKcPJmWby4crOq/OwqvwPm5M19+MiLZKpghJphV+WxrzJYVm5WlZ9d5VfNu/JdIQBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASKZ/pQcoV6HQu8dlWdo5aoFdlac3+7Kr8hyM+7Kr/Owqv2relSMWAEAywgIASKbssFi9enWcffbZMXLkyCgUCvGb3/ymD8YCAGpR2WHR2dkZzc3Nccstt/TFPABADSv75M0ZM2bEjBkzct+/VCpFqVTqvt7R0VHuUwIANaLPz7FobW2NYrHYfWlqaurrpwQAKqTPw2L+/PnR3t7efWlra+vrpwQAKqTPP8eivr4+6uvr+/ppAIAq4M9NAYBkhAUAkEzZb4Xs2rUrNm/e3H19y5YtsWHDhhgyZEiMHj066XAAQG0pOyzWrVsXp59+evf1efPmRUTEnDlz4u677042GABQe8oOiw9+8IORVfAbXw7GL5vpLbsqj33lZ1f52VV+dpVfNe/KORYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEimf6UHKFeh0LvHZVnaOWqBXZWnN/uyq/IcjPuyq/zsKr9q3pUjFgBAMsICAEimV2Fxyy23xNvf/vYYOHBgTJkyJZ588snUcwEANajssPjVr34V8+bNiwULFsRTTz0Vzc3NcdZZZ8X27dv7Yj4AoIaUHRY33nhjfP7zn4+5c+fG8ccfH7fddlu87W1vi5/85Cd7vH+pVIqOjo4eFwDgwFRWWLz22muxfv36mD59+v/9gH79Yvr06fHEE0/s8TGtra1RLBa7L01NTfs2MQBQtcoKi1deeSV2794dRx11VI/bjzrqqHjxxRf3+Jj58+dHe3t796Wtra330wIAVa3PP8eivr4+6uvr+/ppAIAqUNYRiyOOOCLq6uripZde6nH7Sy+9FMOHD086GABQe8oKiwEDBsRJJ50UjzzySPdtXV1d8cgjj8TUqVOTDwcA1Jay3wqZN29ezJkzJyZPnhynnHJK3HTTTdHZ2Rlz587ti/kAgBpSdlhccMEF8fLLL8fVV18dL774Yrz73e+O3/72t284oRMAOPgUsmz/fn1LR0dHFIvFaG9vj8bGxv351ABAL+X9/e27QgCAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAk0+ffbvrfXv88ro6Ojv391ABAL73+e/utPldzv4fFzp07IyKiqalpfz81ALCPdu7cGcVica//vt8/0rurqyu2bdsWDQ0NUSgUkv3cjo6OaGpqira2Nh8V/hbsKj+7Ko995WdX+dlVfn25qyzLYufOnTFy5Mjo12/vZ1Ls9yMW/fr1i1GjRvXZz29sbPTCy8mu8rOr8thXfnaVn13l11e7erMjFa9z8iYAkIywAACSOWDCor6+PhYsWBD19fWVHqXq2VV+dlUe+8rPrvKzq/yqYVf7/eRNAODAdcAcsQAAKk9YAADJCAsAIBlhAQAkIywAgGQOmLC45ZZb4u1vf3sMHDgwpkyZEk8++WSlR6o6q1evjrPPPjtGjhwZhUIhfvOb31R6pKrV2toaJ598cjQ0NMSwYcPi3HPPjWeeeabSY1WlW2+9NSZNmtT9SX9Tp06NZcuWVXqsmrBo0aIoFApx5ZVXVnqUqrRw4cIoFAo9Lscdd1ylx6paW7dujYsuuiiGDh0agwYNiokTJ8a6dev2+xwHRFj86le/innz5sWCBQviqaeeiubm5jjrrLNi+/btlR6tqnR2dkZzc3PccsstlR6l6q1atSpaWlpizZo1sWLFivj3v/8dZ555ZnR2dlZ6tKozatSoWLRoUaxfvz7WrVsXH/rQh+LjH/94/PGPf6z0aFVt7dq1cfvtt8ekSZMqPUpVO+GEE+KFF17ovvz+97+v9EhV6dVXX41p06bFIYccEsuWLYs//elPccMNN8TgwYP3/zDZAeCUU07JWlpauq/v3r07GzlyZNba2lrBqapbRGRLly6t9Bg1Y/v27VlEZKtWrar0KDVh8ODB2Y9+9KNKj1G1du7cmY0bNy5bsWJF9oEPfCC74oorKj1SVVqwYEHW3Nxc6TFqwle+8pXs1FNPrfQYWZZlWc0fsXjttddi/fr1MX369O7b+vXrF9OnT48nnniigpNxIGlvb4+IiCFDhlR4kuq2e/fuWLJkSXR2dsbUqVMrPU7VamlpiY997GM9/r/Fnm3atClGjhwZ73jHO2L27Nnx/PPPV3qkqvTAAw/E5MmTY+bMmTFs2LA48cQT484776zILDUfFq+88krs3r07jjrqqB63H3XUUfHiiy9WaCoOJF1dXXHllVfGtGnTYsKECZUepypt3LgxDjvssKivr48vfvGLsXTp0jj++OMrPVZVWrJkSTz11FPR2tpa6VGq3pQpU+Luu++O3/72t3HrrbfGli1b4rTTToudO3dWerSq89e//jVuvfXWGDduXCxfvjy+9KUvxeWXXx4//elP9/ss+/1r06HWtLS0xNNPP+293Tdx7LHHxoYNG6K9vT3uvffemDNnTqxatUpc/Je2tra44oorYsWKFTFw4MBKj1P1ZsyY0f3fkyZNiilTpsSYMWPinnvuiYsvvriCk1Wfrq6umDx5clx//fUREXHiiSfG008/HbfddlvMmTNnv85S80csjjjiiKirq4uXXnqpx+0vvfRSDB8+vEJTcaC49NJL46GHHorHHnssRo0aVelxqtaAAQPine98Z5x00knR2toazc3N8b3vfa/SY1Wd9evXx/bt2+M973lP9O/fP/r37x+rVq2K73//+9G/f//YvXt3pUesaocffniMHz8+Nm/eXOlRqs6IESPeEPLvete7KvLWUc2HxYABA+Kkk06KRx55pPu2rq6ueOSRR7zHS69lWRaXXnppLF26NB599NEYO3ZspUeqKV1dXVEqlSo9RtU544wzYuPGjbFhw4buy+TJk2P27NmxYcOGqKurq/SIVW3Xrl3xl7/8JUaMGFHpUarOtGnT3vAn8c8++2yMGTNmv89yQLwVMm/evJgzZ05Mnjw5TjnllLjpppuis7Mz5s6dW+nRqsquXbt6lP6WLVtiw4YNMWTIkBg9enQFJ6s+LS0tsXjx4rj//vujoaGh+3ydYrEYgwYNqvB01WX+/PkxY8aMGD16dOzcuTMWL14cK1eujOXLl1d6tKrT0NDwhvN0Dj300Bg6dKjzd/bgqquuirPPPjvGjBkT27ZtiwULFkRdXV3MmjWr0qNVnS9/+cvxvve9L66//vo4//zz48knn4w77rgj7rjjjv0/TKX/LCWVm2++ORs9enQ2YMCA7JRTTsnWrFlT6ZGqzmOPPZZFxBsuc+bMqfRoVWdPe4qI7K677qr0aFXnc5/7XDZmzJhswIAB2ZFHHpmdccYZ2cMPP1zpsWqGPzfduwsuuCAbMWJENmDAgOzoo4/OLrjggmzz5s2VHqtqPfjgg9mECROy+vr67LjjjsvuuOOOisxRyLIs2/85AwAciGr+HAsAoHoICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAk8z869ike6n3a7gAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGzCAYAAABzfl4TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAe1UlEQVR4nO3df5DU9X348ddyyEH0bgWUX+EgaEHDzyqn9IJJmkh1qIOxTdBQbC5Em0l7RA11JmGmEZwaj0wn1tRSfyXVjiklxpb4Y0oIMQLNKCNgmIGkVUgxUhTQVO7gOjmTu8/3j3699grofo73sbvweMzsjLvsj9e8Zsd7zmc/t1fIsiwLAIAEBpR7AADg1CEsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFUJW2bNkSixcvjilTpsSZZ54Z48aNi2uvvTZeeumlo+77/PPPx5/8yZ/EzJkz44wzzohCoVCGieH0ICyAqvTVr341/vEf/zEuv/zy+PrXvx6f/exnY9OmTXHxxRfHzp07e933n//5n+Mb3/hGFAqFOO+888o0MZweCv4IGXCydXR0xJlnnnlCz/Hss89GY2NjDBo0qOe2Xbt2xbRp0+ITn/hEfOtb3+q5/cCBA1FfXx9DhgyJxYsXx8qVK8P/+qB/OGIBVWLDhg1RKBTi0Ucfjdtvvz3e+973Rl1dXXziE5+Itra26OzsjFtuuSVGjBgRZ511VixatCg6OzuPep5vfetbMXPmzBgyZEgMGzYsPvnJT8bevXt73WfXrl3x8Y9/PEaNGhWDBw+OsWPHxic/+cloa2uLiIiXX345CoVCPPzww0c9f6FQiOXLl/dcX758eRQKhfjpT38af/AHfxBDhw6Nyy67LNc8x/KBD3ygV1REREycODGmTJkS//qv/9rr9pEjR8aQIUPe9TmBEzew3AMA+bS2tsaQIUPiS1/6UuzevTvuueeeOOOMM2LAgAHx5ptvxvLly2Pz5s3x8MMPx4QJE+K2227reexXvvKV+PKXvxzXXntt3HjjjfH666/HPffcEx/60Ifixz/+cZx99tnx1ltvxZVXXhmdnZ3x+c9/PkaNGhX79u2Lp556Kg4dOhTFYrFPc8+fPz8mTpwYd955Z8/RglLmySPLsjhw4EBMmTKlTzMCCWRAVXjmmWeyiMimTp2avfXWWz23L1iwICsUCtncuXN73b+pqSkbP358z/WXX345q6mpyb7yla/0ut+OHTuygQMH9tz+4x//OIuI7Dvf+c5xZ9mzZ08WEdlDDz101L9FRLZs2bKe68uWLcsiIluwYEGv+5U6Tx6PPPJIFhHZN7/5zePep6WlJfO/Pug/PgqBKvOpT30qzjjjjJ7rs2bNiizL4jOf+Uyv+82aNSv27t0bv/71ryMi4p/+6Z+iu7s7rr322njjjTd6LqNGjYqJEyfGM888ExHRc0Ri3bp18V//9V/J5v7c5z7X63qp85Tq3/7t36KlpSWampqiubk52dxAPj4KgSozbty4XtffDoGGhoajbu/u7o62trYYPnx47Nq1K7Isi4kTJx7zed+OlQkTJsSSJUvirrvuir//+7+PD37wg3H11VfH9ddf3+ePQd5+3v+t1HlKsX///rjqqquiWCzGY489FjU1NX2eEzgxwgKqzPF+aB7v9uz/n8/Q3d0dhUIh1q5de8z7nnXWWT3//bWvfS0+/elPx+OPPx7f//7346abborW1tbYvHlzjB079rjfA9HV1XXcuf/vyZN55nknbW1tMXfu3Dh06FD8y7/8S4wZM6akxwH9Q1jAaeL888+PLMtiwoQJMWnSpHe9/7Rp02LatGnxZ3/2Z/Hss8/G7Nmz47777os77rgjhg4dGhERhw4d6vWYn//85/02z7H88pe/jHnz5sVLL70UP/jBD2Ly5Ml9eh4gHedYwGni93//96OmpiZuv/32o77DIcuy+MUvfhEREe3t7T3nZbxt2rRpMWDAgJ5fX62vr49zzjknNm3a1Ot+f/M3f5N8nuPp6uqK6667Lp577rn4zne+E01NTSW/NtB/HLGA08T5558fd9xxRyxdujRefvnluOaaa6Kuri727NkTa9asic9+9rNx6623xg9/+MNYvHhxzJ8/PyZNmhS//vWv45FHHomampr4+Mc/3vN8N954Y6xYsSJuvPHGaGxsjE2bNh3z67RPdJ7j+dM//dN44oknYt68efGf//mfvb4QKyLi+uuv7/nvn//85/HII49ERMTWrVsjIuKOO+6IiIjx48fHH/7hH5Y8N/DOhAWcRr70pS/FpEmT4i//8i/j9ttvj4j/PunziiuuiKuvvjoiImbMmBFXXnllPPnkk7Fv3754z3veEzNmzIi1a9fGb/3Wb/U812233Ravv/56PPbYY/Hoo4/G3LlzY+3atTFixIik8xzP9u3bIyLiySefjCeffPKof//fYbFnz5748pe/3Ovf377+4Q9/WFhAQr7SGwBIxjkWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGRO+vdYdHd3x6uvvhp1dXXH/XsDAEBlybIsDh8+HGPGjIkBA45/XOKkh8Wrr7561F9hBACqw969e2Ps2LHH/feTHhZ1dXUR8d+D1dfXn+yXBwD6oL29PRoaGnp+jh/PSQ+Ltz/+qK+vFxYAUGXe7TQGJ28CAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAyJ/17LE5UX78FPMvSzlEN7CqfvuzLrvI5HfdlV6Wzq9JV8q4csQAAkhEWAEAyucNi3759cf3118fw4cNjyJAhMW3atNi6dWt/zAYAVJlc51i8+eabMXv27PjIRz4Sa9eujXPPPTd27doVQ4cO7a/5AIAqkissvvrVr0ZDQ0M89NBDPbdNmDDhHR/T2dkZnZ2dPdfb29tzjggAVItcH4U88cQT0djYGPPnz48RI0bERRddFA8++OA7Pqa1tTWKxWLPpaGh4YQGBgAqVyHLSv/lk8GDB0dExJIlS2L+/PmxZcuWuPnmm+O+++6L5ubmYz7mWEcsGhoaoq2trU9/Nr2Sf8Wm0thVPn7dtHTeW6Wzq9LZVenKsav29vYoFovv+vM7V1gMGjQoGhsb49lnn+257aabbootW7bEc889l3Sw4w7sjVcyu8pHWJTOe6t0dlU6uypdJYdFro9CRo8eHZMnT+512/vf//545ZVX+jYlAHBKyRUWs2fPjhdffLHXbS+99FKMHz8+6VAAQHXKFRZf+MIXYvPmzXHnnXfG7t27Y9WqVfHAAw9ES0tLf80HAFSRXGFxySWXxJo1a+If/uEfYurUqfHnf/7ncffdd8fChQv7az4AoIrkOnkzhRM9eRMAOPn65eRNAIB3IiwAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZAaWe4C8CoW+PS7L0s5RDewqn77sy67yOR33ZVels6vSVfKuHLEAAJIRFgBAMrnCYvny5VEoFHpdLrzwwv6aDQCoMrnPsZgyZUr84Ac/+J8nGFh1p2kAAP0kdxUMHDgwRo0aVfL9Ozs7o7Ozs+d6e3t73pcEAKpE7nMsdu3aFWPGjInzzjsvFi5cGK+88so73r+1tTWKxWLPpaGhoc/DAgCVrZBlpf/yydq1a+PIkSNxwQUXxGuvvRa333577Nu3L3bu3Bl1dXXHfMyxjlg0NDREW1tb1NfX5x+4gn/FptLYVT5+3bR03luls6vS2VXpyrGr9vb2KBaL7/rzO1dY/F+HDh2K8ePHx1133RU33HBD0sGOxxuvdHaVj7AonfdW6eyqdHZVukoOixP6ddOzzz47Jk2aFLt37z6RpwEAThEnFBZHjhyJn/3sZzF69OhU8wAAVSxXWNx6662xcePGePnll+PZZ5+N3/u934uamppYsGBBf80HAFSRXL9u+h//8R+xYMGC+MUvfhHnnntuXHbZZbF58+Y499xz+2s+AKCK5AqL1atX99ccJTsdT9LpK7vKx75KZ1els6vS2VXpKnlX/lYIAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSGVjuAfIqFPr2uCxLO0c1sKt8+rIvu8rndNyXXZXOrkpXybtyxAIASEZYAADJnFBYrFixIgqFQtxyyy2JxgEAqlmfw2LLli1x//33x/Tp01POAwBUsT6FxZEjR2LhwoXx4IMPxtChQ9/xvp2dndHe3t7rAgCcmvoUFi0tLXHVVVfFnDlz3vW+ra2tUSwWey4NDQ19eUkAoArkDovVq1fHCy+8EK2trSXdf+nSpdHW1tZz2bt3b+4hAYDqkOt7LPbu3Rs333xzrF+/PgYPHlzSY2pra6O2trZPwwEA1SVXWGzbti0OHjwYF198cc9tXV1dsWnTpvjrv/7r6OzsjJqamuRDAgDVIVdYXH755bFjx45ety1atCguvPDC+OIXvygqAOA0lyss6urqYurUqb1uO/PMM2P48OFH3Q4AnH588yYAkMwJ/xGyDRs2JBijdKfjH5vpK7vKx75KZ1els6vS2VXpKnlXjlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACCZgeUeIK9CoW+Py7K0c1QDu8qnL/s6XXflzUV/8LYqXSXvyhELACAZYQEAJJMrLO69996YPn161NfXR319fTQ1NcXatWv7azYAoMrkCouxY8fGihUrYtu2bbF169b46Ec/Gh/72MfiJz/5SX/NBwBUkUKWndipHMOGDYu/+Iu/iBtuuOGY/97Z2RmdnZ0919vb26OhoSHa2tqivr4+9+tV8gkrlcau8nHyZg7eXPQDb6vSlWNX7e3tUSwW3/Xnd5/Psejq6orVq1dHR0dHNDU1Hfd+ra2tUSwWey4NDQ19fUkAoMLlPmKxY8eOaGpqil/+8pdx1llnxapVq+J3f/d3j3t/RyzKx67yccQiB28u+oG3Vekq+YhF7u+xuOCCC2L79u3R1tYWjz32WDQ3N8fGjRtj8uTJx7x/bW1t1NbW5n0ZAKAKnfA5FnPmzInzzz8/7r///pLuX2rxHI+iLZ1d5eOIRQ7eXPQDb6vSVfIRixP+Hovu7u5eH3UAAKevXB+FLF26NObOnRvjxo2Lw4cPx6pVq2LDhg2xbt26/poPAKgiucLi4MGD8alPfSpee+21KBaLMX369Fi3bl38zu/8Tn/NBwBUkVxh8c1vfrO/5ijZ6fhZWl/ZVT72lYNl0Q+8rUpXybvyt0IAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJDMwHIPkFeh0LfHZVnaOaqBXeXTl33ZVT6n477sqnR2VbpK3pUjFgBAMsICAEgmV1i0trbGJZdcEnV1dTFixIi45ppr4sUXX+yv2QCAKpMrLDZu3BgtLS2xefPmWL9+ffzqV7+KK664Ijo6OvprPgCgihSyrO+ncrz++usxYsSI2LhxY3zoQx865n06Ozujs7Oz53p7e3s0NDREW1tb1NfX5x+4gk9YqTR2lY+TN0vnvVU6uyqdXZWuHLtqb2+PYrH4rj+/T+gci7a2toiIGDZs2HHv09raGsVisefS0NBwIi8JAFSwPh+x6O7ujquvvjoOHToUP/rRj457P0csyseu8nHEonTeW6Wzq9LZVekq+YhFn7/HoqWlJXbu3PmOURERUVtbG7W1tX19GQCgivQpLBYvXhxPPfVUbNq0KcaOHZt6JgCgSuUKiyzL4vOf/3ysWbMmNmzYEBMmTOivuQCAKpQrLFpaWmLVqlXx+OOPR11dXezfvz8iIorFYgwZMqRfBgQAqkeukzcLxzlb5KGHHopPf/rTJT1HqSd/HH+G3A+JCCf35HE67irCyZt5eG+Vzq5KZ1elO2VO3jyBr7xIpgJGqBp2lY99lc6uSmdXpbOr0lXyrvytEAAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJDOw3APkVSj07XFZlnaOamBX+fRlX3aVz+m4L7sqnV2VrpJ35YgFAJCMsAAAkskdFps2bYp58+bFmDFjolAoxHe/+91+GAsAqEa5w6KjoyNmzJgRK1eu7I95AIAqlvvkzblz58bcuXNLvn9nZ2d0dnb2XG9vb8/7kgBAlej3cyxaW1ujWCz2XBoaGvr7JQGAMun3sFi6dGm0tbX1XPbu3dvfLwkAlEm/f49FbW1t1NbW9vfLAAAVwK+bAgDJCAsAIJncH4UcOXIkdu/e3XN9z549sX379hg2bFiMGzcu6XAAQHXJHRZbt26Nj3zkIz3XlyxZEhERzc3N8fDDDycbDACoPrnD4rd/+7cjK+NffDkd/9hMX9lVPvZVOrsqnV2Vzq5KV8m7co4FAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSGVjuAfIqFPr2uCxLO0c1sKt8+rIvu8rndNyXXZXOrkpXybtyxAIASEZYAADJ9CksVq5cGe973/ti8ODBMWvWrHj++edTzwUAVKHcYfHtb387lixZEsuWLYsXXnghZsyYEVdeeWUcPHiwP+YDAKpI7rC466674o/+6I9i0aJFMXny5LjvvvviPe95T/zt3/7tMe/f2dkZ7e3tvS4AwKkpV1i89dZbsW3btpgzZ87/PMGAATFnzpx47rnnjvmY1tbWKBaLPZeGhoYTmxgAqFi5wuKNN96Irq6uGDlyZK/bR44cGfv37z/mY5YuXRptbW09l7179/Z9WgCgovX791jU1tZGbW1tf78MAFABch2xOOecc6KmpiYOHDjQ6/YDBw7EqFGjkg4GAFSfXGExaNCgmDlzZjz99NM9t3V3d8fTTz8dTU1NyYcDAKpL7o9ClixZEs3NzdHY2BiXXnpp3H333dHR0RGLFi3qj/kAgCqSOyyuu+66eP311+O2226L/fv3x2/+5m/G9773vaNO6AQATj+FLDu5f76lvb09isVitLW1RX19/cl8aQCgj0r9+e1vhQAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIpt//uun/9fb3cbW3t5/slwYA+ujtn9vv9r2aJz0sDh8+HBERDQ0NJ/ulAYATdPjw4SgWi8f995P+ld7d3d3x6quvRl1dXRQKhWTP297eHg0NDbF3715fFf4u7Kp0dpWPfZXOrkpnV6Xrz11lWRaHDx+OMWPGxIABxz+T4qQfsRgwYECMHTu2356/vr7eG69EdlU6u8rHvkpnV6Wzq9L1167e6UjF25y8CQAkIywAgGROmbCora2NZcuWRW1tbblHqXh2VTq7yse+SmdXpbOr0lXCrk76yZsAwKnrlDliAQCUn7AAAJIRFgBAMsICAEhGWAAAyZwyYbFy5cp43/veF4MHD45Zs2bF888/X+6RKs6mTZti3rx5MWbMmCgUCvHd73633CNVrNbW1rjkkkuirq4uRowYEddcc028+OKL5R6rIt17770xffr0nm/6a2pqirVr15Z7rKqwYsWKKBQKccstt5R7lIq0fPnyKBQKvS4XXnhhuceqWPv27Yvrr78+hg8fHkOGDIlp06bF1q1bT/ocp0RYfPvb344lS5bEsmXL4oUXXogZM2bElVdeGQcPHiz3aBWlo6MjZsyYEStXriz3KBVv48aN0dLSEps3b47169fHr371q7jiiiuio6Oj3KNVnLFjx8aKFSti27ZtsXXr1vjoRz8aH/vYx+InP/lJuUeraFu2bIn7778/pk+fXu5RKtqUKVPitdde67n86Ec/KvdIFenNN9+M2bNnxxlnnBFr166Nn/70p/G1r30thg4devKHyU4Bl156adbS0tJzvaurKxszZkzW2tpaxqkqW0Rka9asKfcYVePgwYNZRGQbN24s9yhVYejQodk3vvGNco9RsQ4fPpxNnDgxW79+ffbhD384u/nmm8s9UkVatmxZNmPGjHKPURW++MUvZpdddlm5x8iyLMuq/ojFW2+9Fdu2bYs5c+b03DZgwICYM2dOPPfcc2WcjFNJW1tbREQMGzaszJNUtq6urli9enV0dHREU1NTucepWC0tLXHVVVf1+v8Wx7Zr164YM2ZMnHfeebFw4cJ45ZVXyj1SRXriiSeisbEx5s+fHyNGjIiLLrooHnzwwbLMUvVh8cYbb0RXV1eMHDmy1+0jR46M/fv3l2kqTiXd3d1xyy23xOzZs2Pq1KnlHqci7dixI84666yora2Nz33uc7FmzZqYPHlyuceqSKtXr44XXnghWltbyz1KxZs1a1Y8/PDD8b3vfS/uvffe2LNnT3zwgx+Mw4cPl3u0ivPv//7vce+998bEiRNj3bp18cd//Mdx0003xd/93d+d9FlO+p9Nh2rT0tISO3fu9NnuO7jgggti+/bt0dbWFo899lg0NzfHxo0bxcX/sXfv3rj55ptj/fr1MXjw4HKPU/Hmzp3b89/Tp0+PWbNmxfjx4+PRRx+NG264oYyTVZ7u7u5obGyMO++8MyIiLrrooti5c2fcd9990dzcfFJnqfojFuecc07U1NTEgQMHet1+4MCBGDVqVJmm4lSxePHieOqpp+KZZ56JsWPHlnucijVo0KD4jd/4jZg5c2a0trbGjBkz4utf/3q5x6o427Zti4MHD8bFF18cAwcOjIEDB8bGjRvjr/7qr2LgwIHR1dVV7hEr2tlnnx2TJk2K3bt3l3uUijN69OijQv79739/WT46qvqwGDRoUMycOTOefvrpntu6u7vj6aef9hkvfZZlWSxevDjWrFkTP/zhD2PChAnlHqmqdHd3R2dnZ7nHqDiXX3557NixI7Zv395zaWxsjIULF8b27dujpqam3CNWtCNHjsTPfvazGD16dLlHqTizZ88+6lfiX3rppRg/fvxJn+WU+ChkyZIl0dzcHI2NjXHppZfG3XffHR0dHbFo0aJyj1ZRjhw50qv09+zZE9u3b49hw4bFuHHjyjhZ5WlpaYlVq1bF448/HnV1dT3n6xSLxRgyZEiZp6ssS5cujblz58a4cePi8OHDsWrVqtiwYUOsW7eu3KNVnLq6uqPO0znzzDNj+PDhzt85hltvvTXmzZsX48ePj1dffTWWLVsWNTU1sWDBgnKPVnG+8IUvxAc+8IG4884749prr43nn38+HnjggXjggQdO/jDl/rWUVO65555s3Lhx2aBBg7JLL70027x5c7lHqjjPPPNMFhFHXZqbm8s9WsU51p4iInvooYfKPVrF+cxnPpONHz8+GzRoUHbuuedml19+efb973+/3GNVDb9uenzXXXddNnr06GzQoEHZe9/73uy6667Ldu/eXe6xKtaTTz6ZTZ06Nautrc0uvPDC7IEHHijLHIUsy7KTnzMAwKmo6s+xAAAqh7AAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDL/D+INf3btKVOIAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGzCAYAAABzfl4TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAepElEQVR4nO3df5DU9X348ddyyEH0bgUE4cKBxIJG+VEjSikmaSLRoVZjm6JSnFwItpPOKRjGmYSZKDhVjk4nltYYFdtqpy0l1inROEOQGoHJKKNAmEIyo5BCvYAKOnoH18lp7j7fP/L12gs//OzxPnYXHo+ZnXH3dm9f85odec5nP7dbyLIsCwCABAaUewAA4PQhLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBXDauP/++6NQKMSkSZOO+tlzzz0XCxYsiEmTJkVNTU1ccMEFp35AOAMIC+C08Itf/CKWL18eZ5999jF/vnr16li9enUUi8VoaGg4xdPBmaPgS8iAcujo6DhuBPTFLbfcEocOHYqurq54++23Y9euXb1+fuDAgRgxYkScddZZ8Qd/8Aexa9eu2LdvX7LnB37NEQuoIhs3boxCoRBPPvlk3HvvvfHxj3886urq4o//+I+jra0tOjs7484774yRI0fGOeecE/Pnz4/Ozs6jfs8///M/x+WXXx5DhgyJYcOGxS233BKtra297rN79+740pe+FKNGjYrBgwfHmDFj4pZbbom2traIiNi3b18UCoV44oknjvr9hUIhli1b1nN92bJlUSgU4mc/+1n8yZ/8SQwdOjSuuuqqkuY5kc2bN8dTTz0VK1euPO59Ghoa4qyzzsr9O4G+GVjuAYDStbS0xJAhQ+Kb3/xm7NmzJx588ME466yzYsCAAfHuu+/GsmXLYsuWLfHEE0/E+PHj45577ul57P333x9333133HTTTXHbbbfFoUOH4sEHH4zPfOYz8ZOf/CTOPffceP/99+Paa6+Nzs7OuOOOO2LUqFGxf//+ePbZZ+O9996LYrHYp7nnzJkTEyZMiOXLl8eHB0vzzHMiXV1dcccdd8Rtt90WkydP7tNcQEIZUDVeeOGFLCKySZMmZe+//37P7XPnzs0KhUI2e/bsXvefMWNGNm7cuJ7r+/bty2pqarL777+/1/127tyZDRw4sOf2n/zkJ1lEZP/2b/923Fn27t2bRUT2+OOPH/WziMiWLl3ac33p0qVZRGRz587tdb+885zId77znaxYLGYHDx7MsizLPvvZz2aXXnrpCR9z3XXX9doLkI63QqAKffnLX+51WH/69OmRZVl89atf7XW/6dOnR2tra/zqV7+KiIh///d/j+7u7rjpppvi7bff7rmMGjUqJkyYEC+88EJERM8RifXr18f//M//JJv7a1/7Wq/reec5nnfeeSfuueeeuPvuu2PEiBHJ5gT6zlshUIXGjh3b6/qHIdDY2HjU7d3d3dHW1hbDhw+P3bt3R5ZlMWHChGP+3g9jZfz48bF48eJ44IEH4l/+5V/i05/+dNxwww1x66239vltkA9/7/+Vd57j+da3vhXDhg2LO+64o88zAWkJC6hCNTU1Jd2e/f/zGbq7u6NQKMS6deuOed9zzjmn57+//e1vx1e+8pV4+umn47nnnouFCxdGS0tLbNmyJcaMGROFQuGYz9XV1XXcuYcMGdLreinz/Kbdu3fHqlWrYuXKlXHgwIGe23/5y1/GBx98EPv27Yv6+voYNmzYcX8HkJ6wgDPIhRdeGFmWxfjx42PixIkfef/JkyfH5MmT41vf+la8+OKLMXPmzHjkkUfivvvui6FDh0ZExHvvvdfrMf/93//db/P8X/v374/u7u5YuHBhLFy48Kifjx8/PhYtWnTCvxQB0nOOBZxB/uiP/ihqamri3nvv7TmK8aEsy+Kdd96JiIj29vae8zI+NHny5BgwYEDPn6/W19fHeeedF5s3b+51v+9+97vJ5zmWSZMmxdq1a4+6XHrppTF27NhYu3ZtLFiwIPcsQBqOWMAZ5MILL4z77rsvlixZEvv27Ysbb7wx6urqYu/evbF27dr4sz/7s7jrrrviRz/6Udx+++0xZ86cmDhxYvzqV7+Kf/qnf4qampr40pe+1PP7brvttlixYkXcdtttMW3atNi8eXO89tpryec5lvPOOy9uvPHGo27/8AjFb/7sP//zP+OZZ56JiIg9e/ZEW1tb3HfffRERMXXq1Lj++utzzw0cn7CAM8w3v/nNmDhxYvz1X/913HvvvRHx65M+r7nmmrjhhhsi4tf/0F577bXxgx/8IPbv3x8f+9jHYurUqbFu3br4nd/5nZ7fdc8998ShQ4fiqaeeiieffDJmz54d69ati5EjRyadJ4Xt27fH3Xff3eu2D683NTUJC0jER3oDAMk4xwIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQzCn/HIvu7u44cOBA1NXVHfe7BgCAypJlWRw+fDgaGhpiwIDjH5c45WFx4MCBo76BEQCoDq2trTFmzJjj/vyUh0VdXV1E/Hqw+vr6U/30AEAftLe3R2NjY8+/48dzysPiw7c/6uvrhQUAVJmPOo3ByZsAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJDMKf8ci5PV108Bz7K0c1QDuypNX/ZlV6U5E/dlV/nZVX6VvCtHLACAZIQFAJBMyWGxf//+uPXWW2P48OExZMiQmDx5cmzdurU/ZgMAqkxJ51i8++67MXPmzPjc5z4X69atixEjRsTu3btj6NCh/TUfAFBFSgqLv/zLv4zGxsZ4/PHHe24bP378CR/T2dkZnZ2dPdfb29tLHBEAqBYlvRXyzDPPxLRp02LOnDkxcuTIuOyyy+Kxxx474WNaWlqiWCz2XBobG09qYACgchWyLP8fnwwePDgiIhYvXhxz5syJV155JRYtWhSPPPJINDU1HfMxxzpi0djYGG1tbX362vRK/hObSmNXpfHnpvl5beVnV/nZVX7l2FV7e3sUi8WP/Pe7pLAYNGhQTJs2LV588cWe2xYuXBivvPJKvPTSS0kHO+7AXni52VVphEV+Xlv52VV+dpVfJYdFSW+FjB49Oi655JJet33yk5+M119/vW9TAgCnlZLCYubMmfHqq6/2uu21116LcePGJR0KAKhOJYXF17/+9diyZUssX7489uzZE6tXr45Vq1ZFc3Nzf80HAFSRksLiiiuuiLVr18a//uu/xqRJk+Iv/uIvYuXKlTFv3rz+mg8AqCIlnbyZwsmevAkAnHr9cvImAMCJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgmYHlHqBUhULfHpdlaeeoBnZVmr7sy65Kcybuy67ys6v8KnlXjlgAAMkICwAgmZLCYtmyZVEoFHpdLr744v6aDQCoMiWfY3HppZfGf/zHf/zvLxhYdadpAAD9pOQqGDhwYIwaNSr3/Ts7O6Ozs7Pnent7e6lPCQBUiZLPsdi9e3c0NDTEJz7xiZg3b168/vrrJ7x/S0tLFIvFnktjY2OfhwUAKlshy/L/8cm6deviyJEjcdFFF8Ubb7wR9957b+zfvz927doVdXV1x3zMsY5YNDY2RltbW9TX15c+cAX/iU2lsavS+HPT/Ly28rOr/Owqv3Lsqr29PYrF4kf++11SWPym9957L8aNGxcPPPBALFiwIOlgx+OFl59dlUZY5Oe1lZ9d5WdX+VVyWJzUn5uee+65MXHixNizZ8/J/BoA4DRxUmFx5MiR+PnPfx6jR49ONQ8AUMVKCou77rorNm3aFPv27YsXX3wx/vAP/zBqampi7ty5/TUfAFBFSvpz01/84hcxd+7ceOedd2LEiBFx1VVXxZYtW2LEiBH9NR8AUEVKCos1a9b01xy5nYkn6fSVXZXGvvKzq/zsKj+7yq+Sd+W7QgCAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkMzAcg9QqkKhb4/LsrRzVAO7Kk1f9mVXpTkT92VX+dlVfpW8K0csAIBkhAUAkMxJhcWKFSuiUCjEnXfemWgcAKCa9TksXnnllXj00UdjypQpKecBAKpYn8LiyJEjMW/evHjsscdi6NChJ7xvZ2dntLe397oAAKenPoVFc3NzXHfddTFr1qyPvG9LS0sUi8WeS2NjY1+eEgCoAiWHxZo1a2L79u3R0tKS6/5LliyJtra2nktra2vJQwIA1aGkz7FobW2NRYsWxYYNG2Lw4MG5HlNbWxu1tbV9Gg4AqC4lhcW2bdvi4MGD8alPfarntq6urti8eXN85zvfic7OzqipqUk+JABQHUoKi6uvvjp27tzZ67b58+fHxRdfHN/4xjdEBQCc4UoKi7q6upg0aVKv284+++wYPnz4UbcDAGcen7wJACRz0l9CtnHjxgRj5HcmftlMX9lVaewrP7vKz67ys6v8KnlXjlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACCZgeUeoFSFQt8el2Vp56gGdlWavuzLrkpzJu7LrvKzq/wqeVeOWAAAyQgLACCZksLi4YcfjilTpkR9fX3U19fHjBkzYt26df01GwBQZUoKizFjxsSKFSti27ZtsXXr1vj85z8fX/ziF+OnP/1pf80HAFSRQpad3Kkcw4YNi7/6q7+KBQsWHPPnnZ2d0dnZ2XO9vb09Ghsbo62tLerr60t+vko+YaXS2FVpnLyZn9dWfnaVn13lV45dtbe3R7FY/Mh/v/t8jkVXV1esWbMmOjo6YsaMGce9X0tLSxSLxZ5LY2NjX58SAKhwJR+x2LlzZ8yYMSN++ctfxjnnnBOrV6+O3//93z/u/R2xKB+7Ko0jFvl5beVnV/nZVX6VfMSi5M+xuOiii2LHjh3R1tYWTz31VDQ1NcWmTZvikksuOeb9a2tro7a2ttSnAQCq0EmfYzFr1qy48MIL49FHH811/7zFczyKNj+7Ko0jFvl5beVnV/nZVX6VfMTipD/Horu7u9dbHQDAmaukt0KWLFkSs2fPjrFjx8bhw4dj9erVsXHjxli/fn1/zQcAVJGSwuLgwYPx5S9/Od54440oFosxZcqUWL9+fXzhC1/or/kAgCpSUlj8/d//fX/NkduZ+F5aX9lVaewrP7vKz67ys6v8KnlXvisEAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJDCz3AKUqFPr2uCxLO0c1sKvS9GVfdlWaM3FfdpWfXeVXybtyxAIASEZYAADJlBQWLS0tccUVV0RdXV2MHDkybrzxxnj11Vf7azYAoMqUFBabNm2K5ubm2LJlS2zYsCE++OCDuOaaa6Kjo6O/5gMAqkghy/p+KsehQ4di5MiRsWnTpvjMZz5zzPt0dnZGZ2dnz/X29vZobGyMtra2qK+vL33gCj5hpdLYVWmcvJmf11Z+dpWfXeVXjl21t7dHsVj8yH+/T+oci7a2toiIGDZs2HHv09LSEsVisefS2Nh4Mk8JAFSwPh+x6O7ujhtuuCHee++9+PGPf3zc+zliUT52VRpHLPLz2srPrvKzq/wq+YhFnz/Horm5OXbt2nXCqIiIqK2tjdra2r4+DQBQRfoUFrfffns8++yzsXnz5hgzZkzqmQCAKlVSWGRZFnfccUesXbs2Nm7cGOPHj++vuQCAKlRSWDQ3N8fq1avj6aefjrq6unjzzTcjIqJYLMaQIUP6ZUAAoHqUdPJm4Thnizz++OPxla98JdfvyHvyx/FnKPkhEeHknlKcibuKcPJmKby28rOr/Owqv9Pm5M2T+MiLZCpghKphV6Wxr/zsKj+7ys+u8qvkXfmuEAAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJDOw3AOUqlDo2+OyLO0cVcGyStKXdZ2hq/LSKoFd5WdX+VXyrhyxAACSERYAQDIlh8XmzZvj+uuvj4aGhigUCvH973+/H8YCAKpRyWHR0dERU6dOjYceeqg/5gEAqljJJ2/Onj07Zs+enfv+nZ2d0dnZ2XO9vb291KcEAKpEv59j0dLSEsVisefS2NjY308JAJRJv4fFkiVLoq2trefS2tra308JAJRJv3+ORW1tbdTW1vb30wAAFcCfmwIAyQgLACCZkt8KOXLkSOzZs6fn+t69e2PHjh0xbNiwGDt2bNLhAIDqUnJYbN26NT73uc/1XF+8eHFERDQ1NcUTTzyRbDAAoPqUHBa/93u/F1kZv/HlTPyymT6zrJJYV352lZ9d5WdX+VXyrpxjAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZAaWe4BSFQp9e1yWpZ2jGtgV/cVrKz+7ys+u8qvkXTliAQAkIywAgGT6FBYPPfRQXHDBBTF48OCYPn16vPzyy6nnAgCqUMlh8b3vfS8WL14cS5cuje3bt8fUqVPj2muvjYMHD/bHfABAFSk5LB544IH40z/905g/f35ccskl8cgjj8THPvax+Id/+Idj3r+zszPa29t7XQCA01NJYfH+++/Htm3bYtasWf/7CwYMiFmzZsVLL710zMe0tLREsVjsuTQ2Np7cxABAxSopLN5+++3o6uqK888/v9ft559/frz55pvHfMySJUuira2t59La2tr3aQGAitbvn2NRW1sbtbW1/f00AEAFKOmIxXnnnRc1NTXx1ltv9br9rbfeilGjRiUdDACoPiWFxaBBg+Lyyy+P559/vue27u7ueP7552PGjBnJhwMAqkvJb4UsXrw4mpqaYtq0aXHllVfGypUro6OjI+bPn98f8wEAVaTksLj55pvj0KFDcc8998Sbb74Zv/3bvx0//OEPjzqhEwA48xSy7NR+fUt7e3sUi8Voa2uL+vr6U/nUAEAf5f3323eFAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEim37/d9Dd9+Hlc7e3tp/qpAYA++vDf7Y/6XM1THhaHDx+OiIjGxsZT/dQAwEk6fPhwFIvF4/78lH+kd3d3dxw4cCDq6uqiUCgk+73t7e3R2NgYra2tPir8I9hVfnZVGvvKz67ys6v8+nNXWZbF4cOHo6GhIQYMOP6ZFKf8iMWAAQNizJgx/fb76+vrvfBysqv87Ko09pWfXeVnV/n1165OdKTiQ07eBACSERYAQDKnTVjU1tbG0qVLo7a2ttyjVDy7ys+uSmNf+dlVfnaVXyXs6pSfvAkAnL5OmyMWAED5CQsAIBlhAQAkIywAgGSEBQCQzGkTFg899FBccMEFMXjw4Jg+fXq8/PLL5R6p4mzevDmuv/76aGhoiEKhEN///vfLPVLFamlpiSuuuCLq6upi5MiRceONN8arr75a7rEq0sMPPxxTpkzp+aS/GTNmxLp168o9VlVYsWJFFAqFuPPOO8s9SkVatmxZFAqFXpeLL7643GNVrP3798ett94aw4cPjyFDhsTkyZNj69atp3yO0yIsvve978XixYtj6dKlsX379pg6dWpce+21cfDgwXKPVlE6Ojpi6tSp8dBDD5V7lIq3adOmaG5uji1btsSGDRvigw8+iGuuuSY6OjrKPVrFGTNmTKxYsSK2bdsWW7dujc9//vPxxS9+MX7605+We7SK9sorr8Sjjz4aU6ZMKfcoFe3SSy+NN954o+fy4x//uNwjVaR33303Zs6cGWeddVasW7cufvazn8W3v/3tGDp06KkfJjsNXHnllVlzc3PP9a6urqyhoSFraWkp41SVLSKytWvXlnuMqnHw4MEsIrJNmzaVe5SqMHTo0Ozv/u7vyj1GxTp8+HA2YcKEbMOGDdlnP/vZbNGiReUeqSItXbo0mzp1arnHqArf+MY3squuuqrcY2RZlmVVf8Ti/fffj23btsWsWbN6bhswYEDMmjUrXnrppTJOxumkra0tIiKGDRtW5kkqW1dXV6xZsyY6OjpixowZ5R6nYjU3N8d1113X6/9bHNvu3bujoaEhPvGJT8S8efPi9ddfL/dIFemZZ56JadOmxZw5c2LkyJFx2WWXxWOPPVaWWao+LN5+++3o6uqK888/v9ft559/frz55ptlmorTSXd3d9x5550xc+bMmDRpUrnHqUg7d+6Mc845J2pra+NrX/tarF27Ni655JJyj1WR1qxZE9u3b4+WlpZyj1Lxpk+fHk888UT88Ic/jIcffjj27t0bn/70p+Pw4cPlHq3i/Nd//Vc8/PDDMWHChFi/fn38+Z//eSxcuDD+8R//8ZTPcsq/Nh2qTXNzc+zatct7uydw0UUXxY4dO6KtrS2eeuqpaGpqik2bNomL39Da2hqLFi2KDRs2xODBg8s9TsWbPXt2z39PmTIlpk+fHuPGjYsnn3wyFixYUMbJKk93d3dMmzYtli9fHhERl112WezatSseeeSRaGpqOqWzVP0Ri/POOy9qamrirbfe6nX7W2+9FaNGjSrTVJwubr/99nj22WfjhRdeiDFjxpR7nIo1aNCg+K3f+q24/PLLo6WlJaZOnRp/8zd/U+6xKs62bdvi4MGD8alPfSoGDhwYAwcOjE2bNsXf/u3fxsCBA6Orq6vcI1a0c889NyZOnBh79uwp9ygVZ/To0UeF/Cc/+cmyvHVU9WExaNCguPzyy+P555/vua27uzuef/557/HSZ1mWxe233x5r166NH/3oRzF+/Phyj1RVuru7o7Ozs9xjVJyrr746du7cGTt27Oi5TJs2LebNmxc7duyImpqaco9Y0Y4cORI///nPY/To0eUepeLMnDnzqD+Jf+2112LcuHGnfJbT4q2QxYsXR1NTU0ybNi2uvPLKWLlyZXR0dMT8+fPLPVpFOXLkSK/S37t3b+zYsSOGDRsWY8eOLeNklae5uTlWr14dTz/9dNTV1fWcr1MsFmPIkCFlnq6yLFmyJGbPnh1jx46Nw4cPx+rVq2Pjxo2xfv36co9Wcerq6o46T+fss8+O4cOHO3/nGO666664/vrrY9y4cXHgwIFYunRp1NTUxNy5c8s9WsX5+te/Hr/7u78by5cvj5tuuilefvnlWLVqVaxaterUD1PuP0tJ5cEHH8zGjh2bDRo0KLvyyiuzLVu2lHukivPCCy9kEXHUpampqdyjVZxj7Skisscff7zco1Wcr371q9m4ceOyQYMGZSNGjMiuvvrq7Lnnniv3WFXDn5se380335yNHj06GzRoUPbxj388u/nmm7M9e/aUe6yK9YMf/CCbNGlSVltbm1188cXZqlWryjJHIcuy7NTnDABwOqr6cywAgMohLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQzP8DMsh9FZibk64AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGzCAYAAABzfl4TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAfOklEQVR4nO3de4yV9Z348c9hkAFkZgTlNmVA6oLWclkEZSlub7IS1sXadlEpXSnq7nYzXigxaUlWwaw6bDatVsuKlxY3u2HRuqVaE4rUCqRRokDJgt0ouLBOQUEbnYFpdqzM8/ujP2czFfQ5w3c458DrlZzEcziXTz6ZOO8855lzClmWZQEAkECvUg8AAJw8hAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAq3rZt2+Lyyy+PQYMGRf/+/WPcuHFx7733drnP008/Hdddd12MGzcuqqqq4uyzzy7NsHCS613qAQCOx9NPPx2zZ8+OSZMmxa233hoDBgyIV199NX796193ud+qVavi0UcfjQsuuCDq6+tLNC2c/Aq+hAwohba2tjj99NOP6zlaW1tj7Nix8alPfSoef/zx6NXr2Adh9+/fH4MHD47TTjst/uIv/iJ27twZe/fuPa7XBz7IWyFQQTZs2BCFQiEee+yxuP322+NjH/tY1NTUxF/+5V9GS0tLtLe3x8KFC2PIkCExYMCAWLBgQbS3t3/gef7t3/4tJk+eHP369YtBgwbF1VdfHc3NzV3us2vXrvjyl78cw4YNi759+8aIESPi6quvjpaWloiI2Lt3bxQKhXjkkUc+8PyFQiGWLl3aeX3p0qVRKBTiV7/6VXzlK1+JgQMHxsUXX1zUPEezatWqOHDgQNx5553Rq1evaGtri46OjqPet76+Pk477bSPfE7g+HgrBCpQU1NT9OvXL771rW/F7t2747777ovTTjstevXqFW+//XYsXbo0Nm/eHI888kiMHj06brvtts7H3nnnnXHrrbfGlVdeGddff328+eabcd9998WnP/3p+OUvfxlnnHFGvPvuuzFz5sxob2+PG2+8MYYNGxb79u2Lp556Kt55552oq6vr1txz5syJMWPGxF133RXvHyzNM8+x/OxnP4va2trYt29fXHHFFfHKK6/E6aefHn/1V38Vd999d/Tt27dbcwLHIQMqxrPPPptFRDZu3Ljs3Xff7bx97ty5WaFQyGbNmtXl/tOmTctGjRrVeX3v3r1ZVVVVduedd3a5344dO7LevXt33v7LX/4yi4jshz/84TFn2bNnTxYR2cqVKz/wbxGRLVmypPP6kiVLsojI5s6d2+V+eec5lgkTJmT9+/fP+vfvn914443Zf/zHf2Q33nhjFhHZ1VdffczHXXbZZV32AqTjrRCoQNdcc02Xw/pTp06NLMvi2muv7XK/qVOnRnNzc7z33nsREfGjH/0oOjo64sorr4y33nqr8zJs2LAYM2ZMPPvssxERnUck1q1bF7/97W+Tzf31r3+9y/W88xzL4cOH47e//W1cc801ce+998aXvvSluPfee+Nv//ZvY/Xq1bFr165kswP5CAuoQCNHjuxy/f0QaGho+MDtHR0dnedF7Nq1K7IsizFjxsTgwYO7XP7rv/4rDh48GBERo0ePjkWLFsXDDz8cZ511VsycOTOWL1/e+TzdNXr06C7X885zLP369YuIiLlz53a5/Stf+UpERDz//PPHNS9QPOdYQAWqqqoq6vbs/5/P0NHREYVCIdauXXvU+w4YMKDzv7/97W/H1772tXjiiSfi6aefjptuuimamppi8+bNMWLEiCgUCkd9rSNHjhxz7vdD4H3FzHM09fX18dJLL8XQoUO73D5kyJCIiHj77bc/9PFAesICTiHnnHNOZFkWo0ePjrFjx37k/cePHx/jx4+Pv//7v4/nnnsupk+fHitWrIg77rgjBg4cGBER77zzTpfH/M///E+PzfOHJk+eHOvXr499+/bFueee23n7/v37IyJi8ODBRT8ncHy8FQKnkC996UtRVVUVt99+e+dRjPdlWRa/+c1vIuL3nw/x/nkZ7xs/fnz06tWr889Xa2tr46yzzopNmzZ1ud8///M/J5/nWK688sqIiPj+97/f5faHH344evfuHZ/97GdzzwKk4YgFnELOOeecuOOOO2Lx4sWxd+/euOKKK6Kmpib27NkTa9asib/5m7+JW265JX7+85/HDTfcEHPmzImxY8fGe++9F//6r/8aVVVV8eUvf7nz+a6//vpYtmxZXH/99TFlypTYtGlTvPLKK8nnOZZJkybFtddeGz/4wQ/ivffei8985jOxYcOG+OEPfxiLFy/u8gmb//mf/xlPPvlkRETs3r07Wlpa4o477oiIiIkTJ8bs2bOLXSdwFMICTjHf+ta3YuzYsXH33XfH7bffHhG/P+nz0ksvjcsvvzwifv+LdubMmfGTn/wk9u3bF/3794+JEyfG2rVr40/+5E86n+u2226LN998Mx5//PF47LHHYtasWbF27drOcxxSzfNhVqxYESNHjoyVK1fGmjVrYtSoUXH33XfHwoULu9xv27Ztceutt3a57f3r8+fPFxaQiI/0BgCScY4FAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIJkT/jkWHR0dsX///qipqTnmdw0AAOUly7I4dOhQ1NfXR69exz4uccLDYv/+/R/4BkYAoDI0NzfHiBEjjvnvJzwsampqIuL3g9XW1p7olwcAuqG1tTUaGho6f48fywkPi/ff/qitrRUWAFBhPuo0BidvAgDJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMif8cyyOV3c/BTzL0s5RCeyqON3Zl10V51Tcl13lZ1f5lfOuHLEAAJIRFgBAMkWHxb59++KrX/1qnHnmmdGvX78YP358bNmypSdmAwAqTFHnWLz99tsxffr0+NznPhdr166NwYMHx65du2LgwIE9NR8AUEGKCot//Md/jIaGhli5cmXnbaNHj/7Qx7S3t0d7e3vn9dbW1iJHBAAqRVFvhTz55JMxZcqUmDNnTgwZMiQmTZoUDz300Ic+pqmpKerq6jovDQ0NxzUwAFC+ClmW/49P+vbtGxERixYtijlz5sSLL74YN998c6xYsSLmz59/1Mcc7YhFQ0NDtLS0dOtr08v5T2zKjV0Vx5+b5udnKz+7ys+u8ivFrlpbW6Ouru4jf38XFRZ9+vSJKVOmxHPPPdd520033RQvvvhiPP/880kHO+bAfvBys6viCIv8/GzlZ1f52VV+5RwWRb0VMnz48Dj//PO73PaJT3wiXnvtte5NCQCcVIoKi+nTp8fLL7/c5bZXXnklRo0alXQoAKAyFRUW3/jGN2Lz5s1x1113xe7du2PVqlXx4IMPRmNjY0/NBwBUkKLC4sILL4w1a9bEv//7v8e4cePiH/7hH+Kee+6JefPm9dR8AEAFKerkzRSO9+RNAODE65GTNwEAPoywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJLpXeoBilUodO9xWZZ2jkpgV8Xpzr7sqjin4r7sKj+7yq+cd+WIBQCQjLAAAJIpKiyWLl0ahUKhy+W8887rqdkAgApT9DkWn/zkJ+NnP/vZ/z1B74o7TQMA6CFFV0Hv3r1j2LBhue/f3t4e7e3tnddbW1uLfUkAoEIUfY7Frl27or6+Pj7+8Y/HvHnz4rXXXvvQ+zc1NUVdXV3npaGhodvDAgDlrZBl+f/4ZO3atXH48OE499xz4/XXX4/bb7899u3bFzt37oyampqjPuZoRywaGhqipaUlamtrix+4jP/EptzYVXH8uWl+frbys6v87Cq/UuyqtbU16urqPvL3d1Fh8YfeeeedGDVqVHznO9+J6667Lulgx+IHLz+7Ko6wyM/PVn52lZ9d5VfOYXFcf256xhlnxNixY2P37t3H8zQAwEniuMLi8OHD8eqrr8bw4cNTzQMAVLCiwuKWW26JjRs3xt69e+O5556LL37xi1FVVRVz587tqfkAgApS1J+b/vrXv465c+fGb37zmxg8eHBcfPHFsXnz5hg8eHBPzQcAVJCiwmL16tU9NUdup+JJOt1lV8Wxr/zsKj+7ys+u8ivnXfmuEAAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJNO71AMUq1Do3uOyLO0clcCuitOdfdlVcU7FfdlVfnaVXznvyhELACAZYQEAJHNcYbFs2bIoFAqxcOHCROMAAJWs22Hx4osvxgMPPBATJkxIOQ8AUMG6FRaHDx+OefPmxUMPPRQDBw780Pu2t7dHa2trlwsAcHLqVlg0NjbGZZddFjNmzPjI+zY1NUVdXV3npaGhoTsvCQBUgKLDYvXq1bFt27ZoamrKdf/FixdHS0tL56W5ubnoIQGAylDU51g0NzfHzTffHOvXr4++ffvmekx1dXVUV1d3azgAoLIUFRZbt26NgwcPxgUXXNB525EjR2LTpk3xve99L9rb26Oqqir5kABAZSgqLC655JLYsWNHl9sWLFgQ5513Xnzzm98UFQBwiisqLGpqamLcuHFdbjv99NPjzDPP/MDtAMCpxydvAgDJHPeXkG3YsCHBGPmdil820112VRz7ys+u8rOr/Owqv3LelSMWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIpnepByhWodC9x2VZ2jkqgV0Vpzv7sqvinIr7sqv87Cq/ct6VIxYAQDLCAgBIpqiwuP/++2PChAlRW1sbtbW1MW3atFi7dm1PzQYAVJiiwmLEiBGxbNmy2Lp1a2zZsiU+//nPxxe+8IV46aWXemo+AKCCFLLs+E7lGDRoUPzTP/1TXHfddUf99/b29mhvb++83traGg0NDdHS0hK1tbVFv145n7BSbuyqOE7ezM/PVn52lZ9d5VeKXbW2tkZdXd1H/v7u9jkWR44cidWrV0dbW1tMmzbtmPdramqKurq6zktDQ0N3XxIAKHNFH7HYsWNHTJs2Lf73f/83BgwYEKtWrYo///M/P+b9HbEoHbsqjiMW+fnZys+u8rOr/Mr5iEXRn2Nx7rnnxvbt26OlpSUef/zxmD9/fmzcuDHOP//8o96/uro6qquri30ZAKACHfc5FjNmzIhzzjknHnjggVz3z1s8x6Jo87Or4jhikZ+frfzsKj+7yq+cj1gc9+dYdHR0dHmrAwA4dRX1VsjixYtj1qxZMXLkyDh06FCsWrUqNmzYEOvWreup+QCAClJUWBw8eDCuueaaeP3116Ouri4mTJgQ69atiz/7sz/rqfkAgApSVFh8//vf76k5cjsV30vrLrsqjn3lZ1f52VV+dpVfOe/Kd4UAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACCZ3qUeoFiFQvcel2Vp56gEdlWc7uzLropzKu7LrvKzq/zKeVeOWAAAyQgLACCZosKiqakpLrzwwqipqYkhQ4bEFVdcES+//HJPzQYAVJiiwmLjxo3R2NgYmzdvjvXr18fvfve7uPTSS6Otra2n5gMAKkghy7p/Ksebb74ZQ4YMiY0bN8anP/3po96nvb092tvbO6+3trZGQ0NDtLS0RG1tbfEDl/EJK+XGrorj5M38/GzlZ1f52VV+pdhVa2tr1NXVfeTv7+M6x6KlpSUiIgYNGnTM+zQ1NUVdXV3npaGh4XheEgAoY90+YtHR0RGXX355vPPOO/GLX/zimPdzxKJ07Ko4jljk52crP7vKz67yK+cjFt3+HIvGxsbYuXPnh0ZFRER1dXVUV1d392UAgArSrbC44YYb4qmnnopNmzbFiBEjUs8EAFSoosIiy7K48cYbY82aNbFhw4YYPXp0T80FAFSgosKisbExVq1aFU888UTU1NTEG2+8ERERdXV10a9fvx4ZEACoHEWdvFk4xtkiK1eujK997Wu5niPvyR/HnqHoh0SEk3uKcSruKsLJm8Xws5WfXeVnV/mdNCdvHsdHXiRTBiNUDLsqjn3lZ1f52VV+dpVfOe/Kd4UAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACCZ3qUeoGiFQvcel2Vp56gAVlWc7uzLropzKu7LrvKzq/zKeVeOWAAAyQgLACCZosNi06ZNMXv27Kivr49CoRA//vGPe2AsAKASFR0WbW1tMXHixFi+fHlPzAMAVLCiT96cNWtWzJo1K/f929vbo729vfN6a2trsS8JAFSIHj/HoqmpKerq6jovDQ0NPf2SAECJ9HhYLF68OFpaWjovzc3NPf2SAECJ9PjnWFRXV0d1dXVPvwwAUAb8uSkAkIywAACSKfqtkMOHD8fu3bs7r+/Zsye2b98egwYNipEjRyYdDgCoLEWHxZYtW+Jzn/tc5/VFixZFRMT8+fPjkUceSTYYAFB5ig6Lz372s5GV8htfTsVvm+kmqyqOfeVnV/nZVX52lV8578o5FgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkhAUAkIywAACSERYAQDLCAgBIRlgAAMkICwAgGWEBACQjLACAZIQFAJCMsAAAkhEWAEAywgIASKZ3qQcoVqHQvcdlWdo5KoFdFac7+7Kr4pyK+7Kr/Owqv3LelSMWAEAywgIASKZbYbF8+fI4++yzo2/fvjF16tR44YUXUs8FAFSgosPi0UcfjUWLFsWSJUti27ZtMXHixJg5c2YcPHiwJ+YDACpI0WHxne98J/76r/86FixYEOeff36sWLEi+vfvHz/4wQ+Oev/29vZobW3tcgEATk5FhcW7774bW7dujRkzZvzfE/TqFTNmzIjnn3/+qI9pamqKurq6zktDQ8PxTQwAlK2iwuKtt96KI0eOxNChQ7vcPnTo0HjjjTeO+pjFixdHS0tL56W5ubn70wIAZa3HP8eiuro6qqure/plAIAyUNQRi7POOiuqqqriwIEDXW4/cOBADBs2LOlgAEDlKSos+vTpE5MnT45nnnmm87aOjo545plnYtq0acmHAwAqS9FvhSxatCjmz58fU6ZMiYsuuijuueeeaGtriwULFvTEfABABSk6LK666qp4880347bbbos33ngj/viP/zh++tOffuCETgDg1FPIshP79S2tra1RV1cXLS0tUVtbeyJfGgDopry/v31XCACQjLAAAJIRFgBAMsICAEhGWAAAyQgLACAZYQEAJCMsAIBkevzbTf/Q+5/H1draeqJfGgDopvd/b3/U52qe8LA4dOhQREQ0NDSc6JcGAI7ToUOHoq6u7pj/fsI/0rujoyP2798fNTU1USgUkj1va2trNDQ0RHNzs48K/wh2lZ9dFce+8rOr/Owqv57cVZZlcejQoaivr49evY59JsUJP2LRq1evGDFiRI89f21trR+8nOwqP7sqjn3lZ1f52VV+PbWrDztS8T4nbwIAyQgLACCZkyYsqqurY8mSJVFdXV3qUcqeXeVnV8Wxr/zsKj+7yq8cdnXCT94EAE5eJ80RCwCg9IQFAJCMsAAAkhEWAEAywgIASOakCYvly5fH2WefHX379o2pU6fGCy+8UOqRys6mTZti9uzZUV9fH4VCIX784x+XeqSy1dTUFBdeeGHU1NTEkCFD4oorroiXX3651GOVpfvvvz8mTJjQ+Ul/06ZNi7Vr15Z6rIqwbNmyKBQKsXDhwlKPUpaWLl0ahUKhy+W8884r9Vhla9++ffHVr341zjzzzOjXr1+MHz8+tmzZcsLnOCnC4tFHH41FixbFkiVLYtu2bTFx4sSYOXNmHDx4sNSjlZW2traYOHFiLF++vNSjlL2NGzdGY2NjbN68OdavXx+/+93v4tJLL422trZSj1Z2RowYEcuWLYutW7fGli1b4vOf/3x84QtfiJdeeqnUo5W1F198MR544IGYMGFCqUcpa5/85Cfj9ddf77z84he/KPVIZentt9+O6dOnx2mnnRZr166NX/3qV/Htb387Bg4ceOKHyU4CF110UdbY2Nh5/ciRI1l9fX3W1NRUwqnKW0Rka9asKfUYFePgwYNZRGQbN24s9SgVYeDAgdnDDz9c6jHK1qFDh7IxY8Zk69evzz7zmc9kN998c6lHKktLlizJJk6cWOoxKsI3v/nN7OKLLy71GFmWZVnFH7F49913Y+vWrTFjxozO23r16hUzZsyI559/voSTcTJpaWmJiIhBgwaVeJLyduTIkVi9enW0tbXFtGnTSj1O2WpsbIzLLrusy/+3OLpdu3ZFfX19fPzjH4958+bFa6+9VuqRytKTTz4ZU6ZMiTlz5sSQIUNi0qRJ8dBDD5VklooPi7feeiuOHDkSQ4cO7XL70KFD44033ijRVJxMOjo6YuHChTF9+vQYN25cqccpSzt27IgBAwZEdXV1fP3rX481a9bE+eefX+qxytLq1atj27Zt0dTUVOpRyt7UqVPjkUceiZ/+9Kdx//33x549e+JP//RP49ChQ6Uerez893//d9x///0xZsyYWLduXfzd3/1d3HTTTfEv//IvJ3yWE/616VBpGhsbY+fOnd7b/RDnnntubN++PVpaWuLxxx+P+fPnx8aNG8XFH2hubo6bb7451q9fH3379i31OGVv1qxZnf89YcKEmDp1aowaNSoee+yxuO6660o4Wfnp6OiIKVOmxF133RUREZMmTYqdO3fGihUrYv78+Sd0loo/YnHWWWdFVVVVHDhwoMvtBw4ciGHDhpVoKk4WN9xwQzz11FPx7LPPxogRI0o9Ttnq06dP/NEf/VFMnjw5mpqaYuLEifHd73631GOVna1bt8bBgwfjggsuiN69e0fv3r1j48aNce+990bv3r3jyJEjpR6xrJ1xxhkxduzY2L17d6lHKTvDhw//QMh/4hOfKMlbRxUfFn369InJkyfHM88803lbR0dHPPPMM97jpduyLIsbbrgh1qxZEz//+c9j9OjRpR6ponR0dER7e3upxyg7l1xySezYsSO2b9/eeZkyZUrMmzcvtm/fHlVVVaUesawdPnw4Xn311Rg+fHipRyk706dP/8CfxL/yyisxatSoEz7LSfFWyKJFi2L+/PkxZcqUuOiii+Kee+6Jtra2WLBgQalHKyuHDx/uUvp79uyJ7du3x6BBg2LkyJElnKz8NDY2xqpVq+KJJ56ImpqazvN16urqol+/fiWerrwsXrw4Zs2aFSNHjoxDhw7FqlWrYsOGDbFu3bpSj1Z2ampqPnCezumnnx5nnnmm83eO4pZbbonZs2fHqFGjYv/+/bFkyZKoqqqKuXPnlnq0svONb3wjPvWpT8Vdd90VV155Zbzwwgvx4IMPxoMPPnjihyn1n6Wkct9992UjR47M+vTpk1100UXZ5s2bSz1S2Xn22WeziPjAZf78+aUerewcbU8Rka1cubLUo5Wda6+9Nhs1alTWp0+fbPDgwdkll1ySPf3006Ueq2L4c9Nju+qqq7Lhw4dnffr0yT72sY9lV111VbZ79+5Sj1W2fvKTn2Tjxo3Lqqurs/POOy978MEHSzJHIcuy7MTnDABwMqr4cywAgPIhLACAZIQFAJCMsAAAkhEWAEAywgIASEZYAADJCAsAIBlhAQAkIywAgGSEBQCQzP8D056buVHeRycAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "\n", + "\n", + "n = 7\n", + "timesteps = 8\n", + "\n", + "\n", + "lattice = np.zeros((n, n))\n", + "center = (n//2, n//2)\n", + "radius = 8\n", + "for i in range(n):\n", + " for j in range(n):\n", + " if np.sqrt((i-center[0])**2 + (j-center[1])**2) <= radius:\n", + " lattice[i,j] = 1\n", + "\n", + "\n", + "x, y = center\n", + "lattice[x, y] = -1\n", + "\n", + "\n", + "colors = np.where(lattice == 1, 'blue', 'white')\n", + "colors[x, y] = 'red'\n", + "\n", + "vacancy_pos = None\n", + "\n", + "for t in range(timesteps):\n", + "\n", + " dx, dy = np.random.choice([-1, 0, 1], 2)\n", + " \n", + " # Check if the new position is valid\n", + " if 0 <= x+dx < n and 0 <= y+dy < n:\n", + " if lattice[x+dx, y+dy] == 1:\n", + "\n", + " lattice[x, y] = 1\n", + " x, y = x+dx, y+dy\n", + " if vacancy_pos is None:\n", + " neighbors = []\n", + " for i in [-1, 0, 1]:\n", + " for j in [-1, 0, 1]:\n", + " if (i, j) != (0, 0) and 0 <= x+i < n and 0 <= y+j < n and lattice[x+i, y+j] == 1:\n", + " neighbors.append((i, j))\n", + " if neighbors:\n", + " dx_vac, dy_vac = neighbors[np.random.randint(len(neighbors))]\n", + " vacancy_pos = (x+dx_vac, y+dy_vac)\n", + " lattice[vacancy_pos[0], vacancy_pos[1]] = 0\n", + " else:\n", + " x_vac, y_vac = vacancy_pos\n", + " lattice[x_vac, y_vac] = 1\n", + " vacancy_pos = None\n", + " elif np.random.random() < 0.8:\n", + " \n", + " if vacancy_pos is None:\n", + " neighbors = []\n", + " for i in [-1, 0, 1]:\n", + " for j in [-1, 0, 1]:\n", + " if (i, j) != (0, 0) and 0 <= x+i < n and 0 <= y+j < n and lattice[x+i, y+j] == 1:\n", + " neighbors.append((i, j))\n", + " if neighbors:\n", + " dx_vac, dy_vac = neighbors[np.random.randint(len(neighbors))]\n", + " vacancy_pos = (x+dx_vac, y+dy_vac)\n", + " lattice[vacancy_pos[0], vacancy_pos[1]] = 0\n", + " x, y = x+dx, y+dy\n", + " \n", + "\n", + " colors = np.where(lattice == 1, 'blue', 'white')\n", + " colors[x, y] = 'red'\n", + " if t % 2 == 0:\n", + "\n", + " plt.scatter(*np.indices((n, n)), marker='s', edgecolors='none',\n", + " c=colors.flatten(), s=100)\n", + " plt.title(f\" mesure {t*10+1}\")\n", + " plt.show()\n" + ] + }, + { + "cell_type": "raw", + "id": "c5441288-b5d6-4591-af32-57bb8efcd699", + "metadata": {}, + "source": [ + "Random walk" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "id": "2bde4492-1d76-4f74-8a27-9b120500f612", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "La distance entre le point de départ et le point d arrivée est 47.34.\n", + "La vitesse de la particule est 593.89 m/s.\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAHHCAYAAAC/R1LgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAABqrUlEQVR4nO3dd1xV9f8H8NcFmSKICKKJ2xRx5grcWzP3Ss2VZimSWzPLVc7MkQstc5TmzG0qmnvkFvc3CwQHbkBRhnB+f7x/XLiACngv547X8/G4D+Dcc+993ytyX/czNYqiKCAiIiKyAFZqF0BERESUUxh8iIiIyGIw+BAREZHFYPAhIiIii8HgQ0RERBaDwYeIiIgsBoMPERERWQwGHyIiIrIYDD5ERERkMRh8iIxE7969UaxYMVUe+8CBA9BoNDhw4IBe71ej0WDChAl6vU9TVr9+fdSvX1/tMnJcsWLF0Lt3b+3Py5cvh0ajwenTp9UriiwWgw+RgT148ACDBw9G2bJl4eDgAA8PD9SoUQOjR4/Gs2fP1C7vre3cudPows2ECROg0Wi0FxsbGxQrVgxffPEFIiMj1S7PKMyYMQMajQbnzp3TOa4oClxdXaHRaBASEqJzXWxsLOzs7NCtW7ecLJVIr3KpXQCROXv8+DGqVauG6OhofPLJJyhbtiwePXqE4OBgLFq0CAMGDICTkxMA4KeffkJSUpLKFWfdzp07sWDBggzDz4sXL5Arl3p/ZhYtWgQnJyfExMRg3759mDdvHs6ePYsjR46oVpOxqF27NgDgyJEjqFKlivb45cuXERkZiVy5cuHo0aMoXry49rpTp04hPj5ee1siU8TgQ2RAS5cuRVhYGI4ePQo/Pz+d66Kjo2Fra6v92cbGJqfLMzh7e3tVH79jx47Inz8/AOCzzz7DRx99hLVr1+LkyZOoUaOGqrWprVq1arC3t8eRI0cQEBCgPX706FG4ubmhWrVqOHLkCD7++GPtdcmBkcGHTBm7uogM6N9//4W1tTXef//9dNc5OzvrBIO0Y3xCQ0Oh0Wgwc+ZMLFiwACVKlICjoyOaNm2K8PBwKIqCb7/9FoULF4aDgwPatGmDx48f6zzGq8bYpB1zkZHDhw+jU6dOKFKkCOzs7ODl5YWhQ4fixYsXOjUvWLBA+1jJl9c9/rlz59CiRQs4OzvDyckJjRo1wokTJ3TOSR4DcvToUQwbNgzu7u7InTs32rVrhwcPHry27tepU6cOAPl3Sfb48WOMGDECFSpUgJOTE5ydndGiRQtcuHBB57bJ46DWrVuHyZMno3DhwrC3t0ejRo1w48aNdI+1ZMkSlCxZEg4ODqhRowYOHz6cYU33799H3759UaBAAdjb26NSpUpYsWKFzjn6+F1Iy9bWFtWrV8fRo0d1jh89ehS+vr6oVatWhtflzZsX5cuXBwDMnDkTfn5+cHNzg4ODA6pWrYoNGza89nFf5cmTJ6hRowYKFy6M69evZ+s+iDKDLT5EBlS0aFEkJibi119/Ra9evbJ1H6tWrUJ8fDwCAgLw+PFjzJgxA507d0bDhg1x4MABjB49Gjdu3MC8efMwYsQI/PLLL3qpff369Xj+/DkGDBgANzc3nDx5EvPmzcOtW7ewfv16ANKKcufOHQQFBeHXX399431evnwZderUgbOzM0aNGgUbGxssXrwY9evXx8GDB1GzZk2d8wMCAuDq6orx48cjNDQUc+bMwaBBg7B27dpsPafQ0FAAgKurq/bYf//9h82bN6NTp04oXrw47t27h8WLF6NevXq4cuUKChUqpHMf06ZNg5WVFUaMGIGoqCjMmDED3bt3x99//609Z+nSpfjss8/g5+eHIUOG4L///kPr1q2RL18+eHl5ac978eIF6tevjxs3bmDQoEEoXrw41q9fj969eyMyMhKDBw/WeWx9/y7Url0bhw8fRmhoqDZ0Hz16FP369UONGjUwfvx4REZGIm/evFAUBceOHYOvry+srOQz89y5c9G6dWt0794d8fHxWLNmDTp16oTt27ejZcuWmf53efjwIZo0aYLHjx/j4MGDKFmyZKZvS5RlChEZTEREhOLu7q4AUMqWLat8/vnnyurVq5XIyMh05/bq1UspWrSo9ueQkBAFgOLu7q5z/pgxYxQASqVKlZSEhATt8a5duyq2trZKbGys9hgAZfz48ekeq2jRokqvXr20P+/fv18BoOzfv1977Pnz5+luN3XqVEWj0Sg3b97UHvP391de9ack7eO3bdtWsbW1Vf7991/tsTt37ih58uRR6tatqz22bNkyBYDSuHFjJSkpSXt86NChirW1dYavX2rjx49XACjXr19XHjx4oISGhiq//PKL4uDgoLi7uysxMTHac2NjY5XExESd24eEhCh2dnbKpEmTtMeSXyNvb28lLi5Oe3zu3LkKAOXixYuKoihKfHy84uHhoVSuXFnnvCVLligAlHr16mmPzZkzRwGg/Pbbb9pj8fHxiq+vr+Lk5KRER0dr63nb34WM7NixQwGg/Prrr4qiKMrdu3cVAMrBgweVp0+fKtbW1sqOHTsURVGUS5cuKQCUyZMna2+f9nckPj5eKV++vNKwYUOd42l/35L/fU+dOqXcvXtX8fHxUUqUKKGEhoa+tl4ifWBXF5EBFShQABcuXMDnn3+OJ0+eIDAwEN26dYOHhwe+/fZbKIryxvvo1KkTXFxctD8nt4p8/PHHOgOHa9asifj4eNy+fVsvtTs4OGi/j4mJwcOHD+Hn5wdFUdLNBMqMxMRE7NmzB23btkWJEiW0xwsWLIhu3brhyJEjiI6O1rlN//79dbrO6tSpg8TERNy8eTNTj1mmTBm4u7ujWLFi+OSTT1CqVCn8+eefcHR01J5jZ2enbcFITEzEo0eP4OTkhDJlyuDs2bPp7rNPnz46Y7OSu8/+++8/AMDp06dx//59fP755zrn9e7dW+ffEZCB4Z6enujatav2mI2NDb744gs8e/YMBw8e1Dlf378Lfn5+sLKy0o7dOXr0KGxsbFC9enU4OTmhYsWK2u6u5K+px/ek/h158uQJoqKiUKdOnQxft4zcunUL9erVQ0JCAg4dOoSiRYtm6nZEb4PBh8jAChYsiEWLFuHu3bu4fv06fvzxR7i7u2PcuHFYunTpG29fpEgRnZ+T3/hSd5mkPv7kyRO91B0WFobevXsjX758cHJygru7O+rVqwcAiIqKyvL9PXjwAM+fP0eZMmXSXeft7Y2kpCSEh4frHE/73JO7qDL7HDdu3IigoCCsXr0a77//Pu7fv6/zZg0ASUlJmD17NkqXLg07Ozvkz58f7u7uCA4OzvB5vqmm5FBWunRpnfNsbGx0Al/yuaVLl9YGr2Te3t469/Wqx37b34W8efPCx8dHJ9xUqVJF+xr5+fnpXGdra6szKHz79u14//33YW9vj3z58sHd3R2LFi3K9O9Hjx49cP/+fRw8eBDvvPNOpm5D9LYYfIhyiEajwbvvvouAgAAcOnQIVlZWWLVq1RtvZ21tnaXjmWlFSkxMfOP1TZo0wY4dOzB69Ghs3rwZQUFBWL58OQDk2LT7t3mOAFC3bl00btwYXbt2RVBQEBwcHNC9e3ed+qdMmYJhw4ahbt26+O2337B7924EBQXBx8cnw+f5tjW9DUP8LtSuXVs7hT3t7EM/Pz+cPHkSCQkJOHLkCKpWraodkH/48GG0bt0a9vb2WLhwIXbu3ImgoCB069Yt069F+/btERkZiblz52bqfCJ94OBmIhWUKFECrq6uuHv3rkEfx9XVNd2CffHx8W983IsXL+J///sfVqxYgZ49e2qPBwUFpTs3dVfU67i7u8PR0THDGTvXrl2DlZVVupYLfXJycsL48ePRp08frFu3Dh999BEAYMOGDWjQoEG61rfIyEjtVPisSO6u+eeff9CwYUPt8YSEBISEhKBSpUo65wYHByMpKUmn1efatWs692VItWvXxqJFi7B3716cO3cOI0eO1F7n5+eHFy9eYMeOHfjvv//QoUMH7XUbN26Evb09du/eDTs7O+3xZcuWZfqxAwICUKpUKYwbNw4uLi748ssv9fOkiF6DLT5EBvT3338jJiYm3fGTJ0/i0aNHGXb76FPJkiVx6NAhnWNLlix5Y4tPcgtC6k/uiqJk+Mk8d+7cAPDGFZGtra3RtGlTbNmyRTu7CgDu3buH1atXo3bt2nB2dn7tfbyt7t27o3Dhwpg+fbpOXWlbKNavX5/tsVLVqlWDu7s7AgMDER8frz2+fPnydK/RBx98gIiICJ1Zai9fvsS8efPg5OSk7Vo0pOQxO7NmzUJCQoJOi0+xYsVQsGBBzJgxQ+dcQF43jUaj87sUGhqKzZs3Z+nxv/nmG4wYMQJjxozBokWL3uKZEGUOW3yIDOjXX3/FqlWr0K5dO1StWhW2tra4evUqfvnlF9jb2+Orr74y6OP369cPn3/+OTp06IAmTZrgwoUL2L179xtbMsqWLYuSJUtixIgRuH37NpydnbFx48YMx4xUrVoVAPDFF1+gWbNmsLa21rampPXdd98hKCgItWvXxsCBA5ErVy4sXrwYcXFx2jdXQ7KxscHgwYMxcuRI7Nq1C82bN8eHH36ISZMmoU+fPvDz88PFixexatWqdONxsvIY3333HT777DM0bNgQXbp0QUhICJYtW5buPvv374/Fixejd+/eOHPmDIoVK4YNGzbg6NGjmDNnDvLkyaOPp/1aRYoUgZeXF44fP45ixYqlm77v5+eHjRs3QqPRoFatWtrjLVu2xKxZs9C8eXN069YN9+/fx4IFC1CqVCkEBwdnqYbvv/8eUVFR8Pf3R548eXQWTSTSN7b4EBnQZ599hk8//RQXLlzAuHHjMHDgQKxZswZNmzbFsWPH0q1bo2+ffvopRo8ejUOHDmH48OEICQlBUFCQtpXmVWxsbLBt2zZUrlwZU6dOxcSJE1G6dGmsXLky3bnt27dHQEAAdu3ahR49eujMUErLx8cHhw8fRvny5bX3W7RoUezfv9/gr0Wy/v37w8XFBdOmTQMAfPXVVxg+fDh2796NwYMH4+zZs9ixY8dbdbv1798fCxcuxJ07dzBy5EgcPnwYW7duTXefDg4OOHDgALp3744VK1Zg+PDhePz4MZYtW5ZuDR9DSm7JSbu6OABt2Clbtizc3Ny0xxs2bIilS5ciIiICQ4YMwe+//47p06ejXbt22aohMDAQH330Efr06YMtW7Zk6z6IMkOj5MSIPCIiIiIjwBYfIiIishgMPkRERGQxGHyIiIjIYjD4EBERkcVg8CEiIiKLweBDREREFoMLGKaRlJSEO3fuIE+ePJleip+IiIjUpSgKnj59ikKFCqXb+Dc1Bp807ty5Y9D9goiIiMhwwsPDUbhw4Vdez+CTRvIS8eHh4QbfN4iIiIj0Izo6Gl5eXm/c6oXBJ43k7i1nZ2cGHyIiIhPzpmEqHNxMREREFoPBh4iIiCwGgw8RERFZDAYfIiIishgmG3ymTZsGjUaDIUOGaI/FxsbC398fbm5ucHJyQocOHXDv3j31iiQiIiKjYpLB59SpU1i8eDEqVqyoc3zo0KHYtm0b1q9fj4MHD+LOnTto3769SlUSERGRsTG54PPs2TN0794dP/30E1xdXbXHo6KisHTpUsyaNQsNGzZE1apVsWzZMhw7dgwnTpxQsWIiIiIyFiYXfPz9/dGyZUs0btxY5/iZM2eQkJCgc7xs2bIoUqQIjh8//sr7i4uLQ3R0tM6FiIiIzJNJLWC4Zs0anD17FqdOnUp3XUREBGxtbZE3b16d4wUKFEBERMQr73Pq1KmYOHGivkslIiIiI2QyLT7h4eEYPHgwVq1aBXt7e73d75gxYxAVFaW9hIeH6+2+9So2Fvj1V6BDB6B+ffn6669ynIiIiDLFZFp8zpw5g/v37+O9997THktMTMShQ4cwf/587N69G/Hx8YiMjNRp9bl37x48PT1feb92dnaws7MzZOlvb+tWoHdv4MkTwMoKSEqSr3/8AQweDKxYAbRqpXaVRERERs9kWnwaNWqEixcv4vz589pLtWrV0L17d+33NjY22Ldvn/Y2169fR1hYGHx9fVWs/C1t3Qq0bQtERsrPSUm6XyMjgTZt5DwiIiJ6LZNp8cmTJw/Kly+vcyx37txwc3PTHu/bty+GDRuGfPnywdnZGQEBAfD19cX777+vRslvLzZWWnoAQFEyPkdRAI1GzrtzB9BjNyAREZG5MZkWn8yYPXs2PvzwQ3To0AF169aFp6cn/vjjD7XLyr7166V761WhJ5miyHkbNuRMXURERCZKoyhvele1LNHR0XBxcUFUVBScnZ3VLaZDB2Dz5pRurdexspIusY0bDV0VERGR0cns+7dZtfiYnUePMhd6ADnv8WPD1kNERGTiGHyMmZubtORkhpUVkC+fYeshIiIycQw+xqxt26y1+LRrZ9ByiIiITB2DjzHr1AlwdZVZW6+j0ch5HTvmTF1EREQmisHHmNnby+KEAJRXhZ/k4ytWcCo7ERHRGzD4GLtWrRC5fDOiNHkBAErymJ/kr3nzAlu2cOVmIiKiTDCZBQwt1YsXQPOFrRHhfgdnxmyA26FNMnsrXz4Z09OxI1t6iIiIMonBx4glJcmCzMHBwKFD9nCr9jEw+GO1yyIiIjJZDD5GbPx4YN06WZOwWjW1qyEiIjJ9HONjpFauBL77Dpg2DWjfXu1qiIiIzAODjxE6fBjo1w/45BNg1Ci1qyEiIjIfDD5G5sYNGbNcqxawaNGbl/AhIiKizGPwMSJPngAtW8pOFRs3Ara2aldERERkXji42UjEx8tm7A8fAidOcNstIiIiQ2DwMQKKAgwcCBw5AuzdC5QurXZFRERE5onBxwjMnAksXQosXw7Urat2NUREROaLY3xUtmkTMHo08NVXQK9ealdDRERk3hh8VHTmDNC9u+w68e23aldDRERk/hh8VHLrluwrWqGCbKxuxX8JIiIig+PbrQqePQM+/BCwsZGN1R0c1K6IiIjIMnBwcw5LTAS6dgX++w84ehTw9FS7IiIiIsvB4JPDRo4Edu4Etm+Xbi4iIiLKOQw+OSgwEJg9G5g3D2jRQu1qiIiILA/H+OSQv/4CBg1KuRAREVHOY/DJIQcPygrNX3+tdiVERESWi8Enh/j7y+ytH35QuxIiIiLLxeCTQzw8gCFDgPnzgbt31a6GiIjIMjH45KARIwA7O2DyZLUrISIiskwMPjkob16Zzr5kCRAaqnY1RERElofBJ4d98QXg6gpMmqR2JURERJaHwSeHOTkBY8bI/lzXr6tdDRERkWVh8FHB558DhQoB48erXQkREZFlYfBRgb098M03wNq1wIULaldDRERkORh8VNKnD1CypAQgIiIiyhkMPiqxsQEmTAC2bQNOnFC7GiIiIsvA4KOirl2BcuW4jQUREVFOYfBRkbU18O23wL59wP79aldDRERk/hh8VNauHVC1KjB2rGxiSkRERIbD4KMyjQb47jvg+HFg5061qyEiIjJvDD5GoFkzoHZtmeGVlKR2NUREROaLwccIaDSycem5c8Aff6hdDRERkfli8DESdesCTZsC48YBiYlqV0NERGSeGHyMyHffAVevAqtWqV0JERGReWLwMSLVqwNt28rChvHxaldDRERkfhh8jMy33wKhocAvv6hdCRERkflh8DEy5cvLis7ffgu8eKF2NUREROaFwccITZwI3LsH9OoF7NgBxMSoXREREZF5YPAxQqVKAbNnAydPAh9+CLi6Ag0bAtOny5R3rvVDRESUPRpF4UYJqUVHR8PFxQVRUVFwdnZWtRZFAf75B9i9G9izR/bziokBPDyAJk1k4cOmTYECBVQtk4iISHWZff9m8EnDmIJPWnFxwLFjEoJ275bWHwCoVCklBNWuDdjZqVsnERFRTmPwySZjDj5p3b8PBAVJENqzB4iIABwdgXr1JAg1awaUKSMrQxMREZkzBp9sMqXgk5qiAMHBKa1Bhw/LWkBeXimtQY0by3ghIiIic8Pgk02mGnzSev4cOHgwZXzQ1auAlZUskpjcGlSjBpArl9qVEhERvT0Gn2wyl+CTVliYdIvt3g3s3Qs8eQK4uACNGklrULNmQLFialdJRESUPZl9/zaZ6eyLFi1CxYoV4ezsDGdnZ/j6+uLPP//UXh8bGwt/f3+4ubnByckJHTp0wL1791Ss2LgUKQL07QusWwc8eACcOAEMGybrBfn7A8WLA+++CwQEANu2Ac+eqV0xERGR/plMi8+2bdtgbW2N0qVLQ1EUrFixAt9//z3OnTsHHx8fDBgwADt27MDy5cvh4uKCQYMGwcrKCkePHs3S45hri8/rREbKVPndu+USGgrY2AC1aqW0BlWuLF1lRERExsgiurry5cuH77//Hh07doS7uztWr16Njh07AgCuXbsGb29vHD9+HO+//36m79MSg09qigLcuJEyNuivv2TtIHf3lLWDmjQBChZUu1IiIqIUmX3/NsmhrYmJiVi/fj1iYmLg6+uLM2fOICEhAY0bN9aeU7ZsWRQpUuSNwScuLg5xcXHan6Ojow1au7HTaIDSpeUyaJDMDDt+PKU1aPVqOa9ixZTWoNq1AXt7desmIiLKDJPqvLh48SKcnJxgZ2eHzz//HJs2bUK5cuUQEREBW1tb5M2bV+f8AgUKICIi4rX3OXXqVLi4uGgvXl5eBnwGpsfWVtYFmjIFOHNG1g5atQqoUgX47Tdp/cmXD2jRQrbZuHJFWo2IiIiMkUkFnzJlyuD8+fP4+++/MWDAAPTq1QtXrlx5q/scM2YMoqKitJfw8HA9VWue3N2Bbt2A5cuBO3eACxeASZOAly+BMWMAHx/dgdSPHqldMRERUQqT6uqytbVFqVKlAABVq1bFqVOnMHfuXHTp0gXx8fGIjIzUafW5d+8ePD09X3ufdnZ2sOMeD9mi0UiXV8WKwIgRsnbQoUMyNmjDBuCXX2Ql6bt3AQscLkVEREbIpFp80kpKSkJcXByqVq0KGxsb7Nu3T3vd9evXERYWBl9fXxUrtCyOjoCvr3SPPXggm6nOmcPQQ0RExsNkWnzGjBmDFi1aoEiRInj69ClWr16NAwcOYPfu3XBxcUHfvn0xbNgw5MuXD87OzggICICvr2+WZnRR9sXHA4sWAd9+C7x4AYwcKZc8edSujIiIKIXJBJ/79++jZ8+euHv3LlxcXFCxYkXs3r0bTZo0AQDMnj0bVlZW6NChA+Li4tCsWTMsXLhQ5arNn6JIt9aYMUBICPDJJ8DEiUChQmpXRkRElJ5Jr+NjCJa+jk9WHD4srTp//w188AEwfTpQvrzaVRERkSUyuy0ryHhcuwa0bQvUrQskJAD79gE7dlhm6ImLA7ZuBX78Ue1KiIgoM0ymq4vUd+8eMGEC8NNPQOHCsp7PRx9Z3lYW8fGy4evatcCWLUDympf9+skAbyIiMl4MPvRGMTHADz8A338P5MolXVr+/pa1WnN8vLRsrVsHbN4s+5t5e8tGr7dvyzgnhh4iIuPH4EOv9PIlsGwZMG4c8Pix7Nz+1VeyUrMlSEiQzVvXrgU2bQKePEnZwb5zZ1msUaOREFikiNrVEhFRZjD4UDqKImN2Ro+WLSi6dQO++w4oXlztygzv5UvgwAFp2fnjD1l5ulQpYOBACTsVKkjYSS0sjMGHiMhUMPiQjtOnZabWgQNA/frAihVAtWpqV2VYiYmy4vTatcDGjcDDh0CJEsCnn0rYqVw5fdhJLSwMqFMnx8olIqK3wOBDAGQNnrFjgd9/B8qVA7Zvlynqr3vDN2WJicCRI9Kys2GDbL5arJisQ9S5M/Dee5l/7mzxISIyHQw+Fu7xY2DyZGD+fMDNTWZs9e4tg5jNTVIScOyYtOxs2ABERABeXkCPHhJ2qlfPetCLjpaBzl5eBimZiIj0zAzf3igzYmMl7EyeLONavv5aZijlzq12ZfqVlAScOCEtO+vXy47y77wDdO0qYadGjbebjh8eLl/Z4kNEZBoYfCxMUpJ0Z40dC9y6BfTvD4wfDxQooHZl+qMospp0cti5dQsoWBDo1EnCjq+v/tYeCguTrww+RESmgcHHgvz1lwxcPntWVl7etQsoW1btqvRDUWRg9rp1cgkLAzw9gY4dJezUqmWYhRbDwgBrawlWRERk/Bh8LMDly8CoUcDOnUDNmjKDyRxmISkKcO6cjNlZtw4IDQXc3VPCTp06EkoMKTxcus7McUwUEZE54p9rM3bnjiw+uGyZrMGzbp2EAlOeqaUowIULKS07//4L5M8PdOggYadu3ZwNIZzRRURkWhh8zNDTp8CMGbLNhKMjMHs28PnngK2t2pVlj6IAFy+mhJ1//pHVo9u3BxYtAho0UK/FhcGHiMi0MPiYmXXrgEGDgAcPZAG+778HXFzUrip7Ll9OCTvXrgF580rYmTcPaNgQsLFRu0IJPr6+aldBRESZxeBjZvbulT2lAFmTZ+tWWZDQx0f3a/786tb5KlevpoSdK1cktLVtK61XjRsbV6tVYqLMGGOLDxGR6dAoiqKoXYQxiY6OhouLC6KiouDs7Kx2OdkSHy/dQVeuSKvJ5cvy/f/+J2v2AICHR/ow5OOjTiD63/9Sws7Fi0CePBJ2OncGmjQB7OxyvqbMSF4TqGNHoFUroHx52bHdwUHtyoiILE9m378ZfNIwh+DzKvHxwI0bKUEoORSlDkTu7hkHInd3/dZy40ZK2LlwAXByAlq3lrDTrBlgb6/fxzOE2FigZ0/g5Eng5k05ptEAJUvKa1a+fMrXd9813gBHRGQOGHyyyZyDz6skJOi2ECV//d//5DpAWoJ8fNKHIg+PzD/Of//JgoLr1slaQrlzS0tJ585A8+am3VLy9GnK63bpUsrXO3fkemtrCT+pA5GPD1C6NKfCExHpA4NPNlli8HmVhIT0LURXrgDXr+sGouQglLaFSKORtXWSw87p0xJuPvxQws4HH8isM3P25ElKy1rqQPTggVxvawuUKaPbOuTjI8sPGHoNIiIic8Lgk00MPm+WkCDr56QeP3T5sm4gcnOTbTCuXJFuq5YtJey0bGl++4Flx4MH6cPQpUuy4SkgAdHbO30gKlLEtNdhIiIyFAafbGLwyb6XL6WFKDkIhYfLGjsffigDlun1FAW4ezd9ILp8GXj2TM7Jk0da1VJ3l5UvL1tmMBARkSVj8MkmBh8yNooi6wWlDURXrwIvXsg5efNKAJowAWjUSM1qiYjUkdn3bw6rJDJyGg1QtKhcPvgg5XhiIhASIiFo1CjgyBHg8WP16iQiMgUG2K+aiHKCtTVQqpQsB/DPP8CsWUCnTmpXRURk3Bh8iEzYlCnSvTVtGjB0qNrVEBEZPwYfIhM1cyYwdiwwaRIwerTa1RARmQYGHyITNHcuMHIk8PXXwDffqF0NEZHpYPAhMjELFwJDhsiA5kmT1K6GiMi0MPgQmZCffwb8/SX4TJvGtXuIiLKKwYfIRKxYAfTvDwwcKDO4GHqIiLKOwYfIBKxeDfTpA/TrB8ybx9BDRJRdDD5ERm79eqBnT6BXLyAwELDi/1oiomzjn1AiI7Z5M9CtG9Cli4zvYeghIno7/DNKZKR27JAd7du1k/E91tZqV0REZPoYfIiM0O7dQPv2QMuWwKpVQC7uqkdEpBcMPkRG5q+/gLZtgSZNgLVrARsbtSsiIjIfDD5ERuTQIaBVK6BePWDDBsDWVu2KiIjMC4MPkZE4dgz44APg/feBTZsAe3u1KyIiMj8MPkRG4ORJoEULoGpVYOtWwMFB7YqIiMwTgw+Riu7flwUJmzUDypcHtm8HcudWuyoiIvPFuSJEOezpU1mfZ9UqYO9eWYW5bVtZpydPHrWrIyIybww+RDkgPh7YtUu2nti6FXjxAqhTB1iwAOjYEXBzU7tCIiLLwOBDZCBJScCRI9Kys3498OQJULEiMH480LUrUKSI2hUSEVkeBh8iPVIUIDhYWnZ+/x0IDweKFgU+/1y2nihfXu0KiYgsG4MPkR6EhkrYWb0auHxZuq66dJGw4+vLPbaIiIwFgw9RNj14IF1Yq1bJGjyOjjJIecYMWXWZKy4TERkfBh+iLHj2DNiyRVp2du+WGVnNmkn4adOGU9GJiIwdgw/RGyQkSMhZvVpCz/PnQK1asv5Op05A/vxqV0hERJnF4EOUgaQk6b5KnpH16JEMTP7mG+Cjj4BixdSukIiIsoPBhyiVixdTZmTdvAl4eQH9+gHduwMVKqhdHRERvS0GH7J4N29K0Fm9WoJPvnxA584yI6tWLc7IIiIyJww+ZJEePgQ2bJCurCNHZFPQNm2AKVOApk0BW1u1KyQiIkNg8CGLERMj20WsXi3bRyiKhJxff5Vp6E5OaldIRESGZjLBZ+rUqfjjjz9w7do1ODg4wM/PD9OnT0eZMmW058TGxmL48OFYs2YN4uLi0KxZMyxcuBAFChRQsXJSk6IAQUHAypWyMWhMDODnB8yZIzOyPDzefB8PHgADBsj3zs5vvuTJk/I91/IhIjIuJhN8Dh48CH9/f1SvXh0vX77EV199haZNm+LKlSvI/f+LpwwdOhQ7duzA+vXr4eLigkGDBqF9+/Y4evSoytWTGoKDgcGDgQMHgHLlgK++kj2yihfP2v1YW8smo3v2AHFxWbutvT3QqhWwbl3WbkdERIahURRFUbuI7Hjw4AE8PDxw8OBB1K1bF1FRUXB3d8fq1avRsWNHAMC1a9fg7e2N48eP4/3338/U/UZHR8PFxQVRUVFwdnY25FMgA3n4UKadL1kClC4NzJ4NNG8uiw2+jZgY4K+/gO3bgR07gNu3ZcHCWrVkp/XatWX15pMngUmTpKWoSBFg1iygQwf9PDciIspYZt+/TabFJ62oqCgAQL58+QAAZ86cQUJCAho3bqw9p2zZsihSpMhrg09cXBziUn2Mj46ONmDVZEgJCcDChcCECdLF9cMPgL+//rqbcueW1ptWreT+L1yQALRjBzBunBxLbeFCoG9fDpQmIjImJjlRNykpCUOGDEGtWrVQ/v+3u46IiICtrS3y5s2rc26BAgUQERHxyvuaOnUqXFxctBcvLy9Dlk4Gsns3ULEiMHSobA76zz/AkCGGG2Oj0QCVKwNjx8rg6A8+SH/OxImyK/sffwBPnxqmDiIiyhqTDD7+/v64dOkS1qxZ89b3NWbMGERFRWkv4eHheqiQcso//0gLTPPmQIECwNmzQGAg4O5u+McOCZEWnTJlgDNngLlzgRcvZDzQgQNAjx7AiRPSzeXmJhuXzpkD3Lhh+NqIiChjJhd8Bg0ahO3bt2P//v0oXLiw9rinpyfi4+MRGRmpc/69e/fg6en5yvuzs7ODs7OzzoWMX1QUMHIk4OMjiw6uXw/s3y+tMIZ28ybQvz/w7rsy3uf774H//gO++EIGM9vYAPXqyfErV4B//5VuN2trYPRoGXdUpgwwbJiMGYqPN3zNREQkTCb4KIqCQYMGYdOmTfjrr79QPM3UnKpVq8LGxgb79u3THrt+/TrCwsLg6+ub0+WSgSQmAkuXSuhYuFDG1ly9CnTs+PaDl9/k1i1g4EAJLps3A1OnSuAZOlQWQHyVEiWAgABZO+jRI7ltvXrA2rVAo0ayyWmnTsDy5cC9e4Z9DkREls5kZnUNHDgQq1evxpYtW3TW7nFxcYHD/7/rDBgwADt37sTy5cvh7OyMgIAAAMCxY8cy/Tic1WW8jhyR6elnz8reWdOmAaka/Qzmzh0JOUuWyBo9I0fKoOm3XfBQUYDz51NmiZ08KcerVwdatpRLlSrcMoOIKDMy+/5tMsFH84qP88uWLUPv3r0BpCxg+Pvvv+ssYPi6rq60GHyMT3i4dBH9/jtQrZqMpfHzM/zj3r0LTJ8uY4YcHYERI6TlJk8ewzze/fvAn39KCNq9G4iOBgoWlIHTH34ING7M1aWJiF7F7IJPTmHwMR7Pn8s4menTZRXkadOAnj0N3wJy7x4wY4Z0pdnby1icwYOlhpySkCAtXDt2SIvQ9esyLb5ePQlBLVsCJUvmXD1ERMaOwSebGHzUpygyWHnkSCAiQsbQfPWV4YPHgwcStObPlwHKQ4fKlPg0KySo4saNlDWDDh6UAdFly0oA+vBDWUSR22MQkSVj8MkmBh91nTsnrSuHD8tu6TNnAqVKGfYxHz2Sx5k3T1qTBg+WVh5XV8M+bnY9fQrs3ZsShCIiABcX2XD1ww9l+vz/7+JCRGQxGHyyicFHHffvA19/Dfz8M+DtLevdNGli2Md8/Fi2k5g7V1qZvvgCGD5c1twxFUlJEhaTu8ROnZLtM/btYwsQEVkWBp9sYvDJWfHx0rU0caK0tkyaJKsdG/JNOzJS9u+aMwd4+RIYNEgGLufEooeGdvgw0LChdNF9/73a1RAR5Ryz36uLTN/OnTKO5sYNCTsTJ8qaNoZy6xaweLF0acXHy5o8o0YBHh6Ge8ycVqeOBJ6hQwFfX6B9e7UrIiIyLgw+lOOuXZMxNH/+CTRoAGzYAFSoYJjHSkqS1ZEXLgS2bpWFBvv1k+nxWVjlwKQMHgwcOwb07g2ULy+LPRIRkchW8Nm3bx/27duH+/fvIykpSee6X375RS+FkfmJjJSurHnzAC8v2byzbVvDrLj85ImshLxokeznVb68PO7HHxtuHR5jodHI6tbVq8tA5xMnONiZiChZloPPxIkTMWnSJFSrVg0FCxZ85cKCRMmSt5kYO1Y28Zw0Sbpi7O31/1inT0vrzpo1Mn6nY0fgl19kurcl/armySPBskYN4LPPZAd5S3r+RESvkuXgExgYiOXLl6NHjx6GqIfMzMGD0vVy4YIsPjh1KlCokH4f48UL2fdq4UKZ1VSkiMwQ69tXdmy3VOXKAT/9BHTrJitdDxyodkVEROrLcvCJj4+HX07sF0Am7eZNWYBw/XppdThxAqhZU7+P8c8/sp3EsmXStdW8uYzj+eAD2QmdgK5dgePHZZZX1ar6/zcgIjI1WV78v1+/fli9erUhaiEzEBMjO6aXLStbLqxcKW+8+nrDfflSdjdv1kwG7S5fLi07N27IYOlWrRh60po5U/Y469hRVqcmIrJkWW7xiY2NxZIlS7B3715UrFgRNmkWXJk1a5beiiPTcuqUTJ++dUvWxpk6VX+bakZEyOKGixfL/desCaxYAXTqJDO16NVsbYF164D33pNd7f/8k+GQiCxXloNPcHAwKleuDAC4dOmSznUc6GzZnj1LWXgwMBA4f162UWjSRFoccmXxt01RgEOHZGbWxo1y3927AwMGyJs4ZV7hwrK7fdOmsl7SpElqV0REpA6u3JwGV25+e//+C+zZAwQFydYJ0dGy0WfDhilBqESJV98+OlpmIS1cCFy5Il1aAwfK4Ghj3T/LVEydKhu+bt8uG5wSEZmLHNmy4tatWwCAwoULZ/cujA6Dj369fCldYMlB6MQJmd5esqQEoKZNZRHDvHmB4GAJO7/9BsTGyho/AwZIYGJjon4kJcnreuQIcOYMULy42hUREemHwYJPUlISvvvuO/zwww949uwZACBPnjwYPnw4xo4dCyurLI+XNioMPoYVFQUcOJAShP75J/05X30lLTzvvJPj5VmEJ0+k6zFvXuDoUcOsp0RElNMMtlfX2LFjsXTpUkybNg21atUCABw5cgQTJkxAbGwsJk+enP2qyey5uABt2gCVKgHOzsC0aenPmTcPuHQppVusdGm2+OiTq6uMmfL1lR3plyxRuyIiopyT5RafQoUKITAwEK1bt9Y5vmXLFgwcOBC3b9/Wa4E5jS0+hpOUBOzeLd1ZO3bI6sK9e8sGpd7e0gV25oy0BAUFSWvEy5dA0aIp3WKNGgH58qn9TMzDsmXAJ5/IytZ9+qhdDRHR2zFYV5e9vT2Cg4PxbpqdD69fv47KlSvjxYsX2avYSDD46N+zZxJ2AgOBkBDAx0cW1Ova9fV7SD17Jis/J3eLXb0qLT9Vq6a0Bvn5yXRtyp7+/WUg+fHjwP9P1iQiMkkGCz41a9ZEzZo18eOPP+ocDwgIwKlTp3DixInsVWwkGHz0b+5cCTqpOTpKl4urq7TgpP76qu9jYmSg9F9/AXv3Ag8fSnCqVy8lCHl7s1ssK2JjZR+zf/+VPb0GDZINZImITI3Bgs/BgwfRsmVLFClSBL6+vgCA48ePIzw8HDt37kSdOnXernKVMfjoX1IScO2aDKp9/Fi+vur71McSEjK+P2dnGZh79276c955RwJQ27Yyloje7P59YMYM2dcrJkYWhRw6VLYaISIyFQadzn7nzh0sWLAA165dAwB4e3tj4MCBKKTv3SdVwOBjHBRF3oRfF4xe9f2TJ7Ka8+3bEpAoc54+lS1A5s6VFiBfXwlA7dplffFJIqKcliPr+JgjBh/Tl5Qkg6I59id7EhNl8Pns2bL0QJEiQEAA0K8fgyQRGS+9Bp/g4GCUL18eVlZWCA4Ofu25FStWzHq1RoTBhyjFuXPSArR6tQTJTz6RKfClSqldGRGRLr0GHysrK0RERMDDwwNWVlbQaDTI6GYajQaJiYlvV7nKGHyI0ouIkJl5ixYBjx4BrVpJN1i9ehxMTkTGQa/B5+bNmyhSpAg0Gg1u3rz52nOLFi2a9WqNCIMP0avFxgKrVkk32OXLMgV+yBDgo48AOzu1qyMiS5bZ9+9M7S9RtGhR7c7rN2/exDvvvIOiRYvqXN555503hiIiMm329kDfvsDFi7K2UqFCsghl0aKy4/uDB2pXSET0elneWKtBgwZ4/PhxuuNRUVFo0KCBXooiIuOm0QCNG8sg6KtXZebXtGmyBlC/frLlCBGRMcpy8FEURdv6k9qjR4+Q+3XL8BKRWSpbVsb+3LoFTJgA7NoFVKgg6ynt3Cmz7IiIjEWmV+do3749ABnA3Lt3b9il6tBPTExEcHAw/Pz89F8hEZmEfPmAL78Ehg8HNmyQcUAtWwJlygCDBwM9e75+ixIiopyQ6RYfFxcXuLi4QFEU5MmTR/uzi4sLPD090b9/f/z222+GrJWITICNjezD9vffwJEj0vqTvBXGmDHSMkREpJYsL2A4ceJEjBgxwmy7tTiri0j/QkOBefOAn38Gnj9P2RajenW1KyMic8GVm7OJwYfIcJ4+BZYtk0UR//sP8POTANS2LbfFIKK3o9fg895772Hfvn1wdXVFlSpVMhzcnOzs2bPZq9hIMPgQGV5iIrB9u4wDOnhQpsMnb4vh4qJ2dURkijL7/p2pz1ht2rTRDmZu27atXgokIstlbQ20aSOXs2eBOXNk/M+ECSnbYpQsqXaVRGSO2NWVBlt8iNRx965sixEYKNtitG4t3WB163JbDCJ6M72u3JxaeHg4bqWalnHy5EkMGTIES5YsyV6lREQAChYEvv0WCAsDliwB/vkHqF8fqFoVWLkSiI9Xu0IiMgdZDj7dunXD/v37AQARERFo3LgxTp48ibFjx2LSpEl6L5CILIuDQ8rqz7t3A56eQK9eMg7ou++4LQYRvZ0sB59Lly6hRo0aAIB169ahQoUKOHbsGFatWoXly5fruz4islAaDdC0qaz+fOWKjAeaMgUoUgT49FPZJJWIKKuyHHwSEhK0A5337t2L1q1bAwDKli2Lu3fv6rc6IiIA3t4y9ic8HBg3TsJQ+fISjP78k9tiEFHmZTn4+Pj4IDAwEIcPH0ZQUBCaN28OALhz5w7c3Nz0XiARUTI3N5n9FRIC/PYb8Pgx8MEHgI+PBKPnz9WukIiMXZaDz/Tp07F48WLUr18fXbt2RaVKlQAAW7du1XaBEREZkq0t0L07cOoUcPgwUK4c4O8v22J89RVw+7baFRKRscrWdPbExERER0fD1dVVeyw0NBSOjo7w8PDQa4E5jdPZiUxTSEjKthgvXgCdO8t0+GrV1K6MiHKCwaazA4C1tTVevnyJI0eO4MiRI3jw4AGKFStm8qGHiExX8eLArFmyCerMmcDx47IXWO3awMaNslo0EVGWg09MTAw++eQTFCxYEHXr1kXdunVRqFAh9O3bF8/ZwU5EKnN2BgYPlnWA/vhDVonu2BEoVQpYs0bt6ohIbVkOPsOGDcPBgwexbds2REZGIjIyElu2bMHBgwcxfPhwQ9RIRJRl1tZAo0ayAaqrq+wQf/iw2lURkdqyPMYnf/782LBhA+rXr69zfP/+/ejcuTMemPjqYhzjQ2T6QkKAH38Eli6V8T5duwJDhgDvvad2ZURkKHrdpDS158+fo0CBAumOe3h4sKuLiFSjKDKuZ9YsYNMmIG9e2fHd3x8oVEjt6ojIWGS5q8vX1xfjx49HbGys9tiLFy8wceJE+Pr66rU4IqI3efkSWLsW8PUFatUCLl4EFiyQxQ4nT2boISJdWW7xmTNnDpo1a4bChQtr1/C5cOEC7O3tsXv3br0XSESUkchImbo+b55sbNqwIbBtmyxoaJWt+apEZAmytY7P8+fPsXr1aly9ehUA4O3tje7du8PBwUHvBeY0jvEhMm7//ZcyficuDujWTdbr+f/PYURkoQwyxufEiRPYtm0b4uPj0bBhQ/Tr1++tCyUiehNFAY4eBWbPBjZvlllaQ4YAAwcCBQuqXR0RmZJMB58NGzagS5cucHBwgI2NDWbNmoXp06djxIgRhqyPiCxYQgKwYYMEnlOngLJlgUWLgI8/Bhwd1a6OiExRpnvCp06dik8//RRRUVF48uQJvvvuO0yZMsWQtRGRhXryBJgxAyhRQrqy8uaVHdkvXwb692foIaLsy3TwuX79OkaMGAFra2sAwPDhw/H06VPcv3/fYMWldejQIbRq1QqFChWCRqPB5s2bda5XFAXjxo1DwYIF4eDggMaNG+Off/7JsfqI6O3cuQN88YVsNvrNN0DTpkBwMLBnD9CiBQctE9Hby/SfkefPn+sMFrK1tYW9vT2ePXtmkMIyEhMTg0qVKmHBggUZXj9jxgz8+OOPCAwMxN9//43cuXOjWbNmOlPvicj4JCUBgYGAtzewejUwYoTM1Fq6FKhQQe3qiMicZGlw888//wwnJyftzy9fvsTy5cuRP39+7bEvvvhCf9Wl0aJFC7Ro0SLD6xRFwZw5c/D111+jTZs2AICVK1eiQIEC2Lx5Mz766COD1UVE2ZfcfXXsGPDpp8D06TJ4mYjIEDIdfIoUKYKffvpJ55inpyd+/fVX7c8ajcagwed1QkJCEBERgcaNG2uPubi4oGbNmjh+/DiDD5GRiY0FpkwBpk2TsTwHDwJ166pdFRGZu0wHn9DQUAOW8fYiIiIAIN12GgUKFNBel5G4uDjExcVpf46OjjZMgUSkdeAA8NlnsqfWV18BY8YAdnZqV0VElsDihwpOnToVLi4u2ouXl5faJRGZrcePgb59gQYNAA8P4Px5YMIEhh4iyjlmE3w8PT0BAPfu3dM5fu/ePe11GRkzZgyioqK0l/DwcIPWSWSJFAX4/XcZvLxxI7B4sXRtlSundmVEZGnMJvgUL14cnp6e2Ldvn/ZYdHQ0/v7779dunmpnZwdnZ2edCxHpT2io7J/VrRtQrx5w9aoMZubUdCJSQ6bH+Ny5cweFVN7m+NmzZ7hx44b255CQEJw/fx758uVDkSJFMGTIEHz33XcoXbo0ihcvjm+++QaFChVC27Zt1SuayEK9fAnMnQuMGwe4uQFbtwKtWqldFRFZukx/5vLx8cHq1asNWcsbnT59GlWqVEGVKlUAAMOGDUOVKlUwbtw4AMCoUaMQEBCA/v37o3r16nj27Bl27doFe3t7NcsmsjhnzgA1agAjR8oU9cuXGXqIyDhkenf2hQsXYvTo0WjevDkWL16MfPnyGbo2VXB3dqLse/YMGD8emDNHFh786SegenW1qyIiS5DZ9+9Mt/gMHDgQwcHBePToEcqVK4dt27bppVAiMg87dwLly8smolOnyqaiDD1EZGyytHJz8eLF8ddff2H+/Plo3749vL29kSuX7l2cPXtWrwUSkXoURfbMWroUKFRI9tAqUkT3q52drLa8fj3QpAnw11+yICERkTHKUvABgJs3b+KPP/6Aq6sr2rRpky74EJF5UBRZXHDaNJmFBQDh4RJswsKAp0/T3+b5c+Drr9MHpCJFZBsKjSZnnwMRUVpZSi0//fQThg8fjsaNG+Py5ctwd3c3VF1EpLKJEyX0/PADMGyY7nXXrgEffQRcuAC4uwODB0voCQuTcHTiBHDrFpCQkHIbR0fdMFSyJDBoEJAnT84+LyKybJkOPs2bN8fJkycxf/589OzZ05A1EZHKJk+W4DN1qm7oiYuTbq3JkyW87N0LNGqU8X0kJQH37kkQSg5Ee/fKWCAAcHAA2raVRQ2JiHJKpoNPYmIigoODUbhwYUPWQ0QqmzlTuqsmTgS+/DLl+JEj0uX1zz/AqFFyjoPDq+/HygooWBDw9ARevABWrgT+/FO2qhg0CBgwAMif3/DPh4gotUwHn6CgIEPWQURG4McfZe2dsWNlUDMAREZKAFq8GKhZEzh7Vqaqv0lCggx4njVL1vXx8QF+/llWcObSWkSkFo5MJiIAMg198GAJPt9+K8c2bAACAoCYGGD+fODzzwFr69ffT2SkrN/z448yzqdJE2DXLqBpUw5uJiL1MfgQEZYuBQYOlOAzfboEFn9/YNs2oE0bCT1v6uUOCZEtKpYulbFA3bvL+KDMtA4REeUUBh8iC7dypWwrMWCAjO/58UcZv+PsDPzxB9Cu3etvf+KEdGdt3AjkzSvhyd9fxvcQERkbBh8iC/b770CfPkDfvhJ+/PyA06clBE2ZAri4ZHy7xERgyxaZ6n7sGFC6tLQK9ewJ5M6ds8+BiCgrGHyILNSGDUCPHkCHDtJSU706ULaszN7y88v4NtHRwIoVshfXf/8BdesCmzfLBqRWmd4Ah4hIPQw+RBZoyxaga1fAzQ04eRLYulWmr48cCdjapj//7FkgMBBYvRqIjQU6dQLWrgWqVcv52omI3gaDD5GF2blTgsvLl8D9+0CDBsCePcC77+qeFxMj4SYwUDYcfecdCUZ9+755oDMRkbFi8CGyIEFBQPv2ssaOq6uM0endW3ea+eXLsmbPypXStdW8ubQQffABwK35iMjU8c8YkYXYvx9o3VqmmnfrBsyeLasoA9J9tXGjBJ7Dh+X4wIEy4Ll4cXXrJiLSJwYfIgvRq5dsHxEYCDRrJsf++QdYsgRYtgx49Aho2FC6t9q2zXisDxGRqWPwIbIQe/fKOB1bW2ndCQyUY66uMqW9f3+gTBm1qyQiMiwGHyILYW8PTJsm+2VFRMiU9ZUrgY4dX7/ZKBHlnMhIWRS0QAGZRJA/P7d60TcGHyIzlpgoO6IvXiyzuXLnlrV7PvsMqFhR7eqIKJmiAGvWAEOGyGzLZHZ2EoBed/Hw4DpaWcHgQ2SG7t6VPbN++gkICwPee0+6trp2BZyc1K6OiFL791+ZTLBnjyw1MWGCLCdx65buJSxMVkq/fRuIj0+5vY2NdGO/Khh5eUkL0ps2GLYUDD5EZiIpCfjrLwk4W7bIH8Nu3WRHdS40SGR84uNlSYlJkySYbN8OtGyZcn316hnfLikJePhQNxSFh6d8f/q0fI2NTbmNtTVQqFDGoSj5+4IFLWPJCgt4ikTm7eFDYPly6c66cQPw8ZGp6h9/LFtREJHxOXZMupyvXgWGDpVWnszuc2dlJd1bHh7SmpsRRQEeP04fipIvwcFy/Plz3fv19Hx1MCpcWMKTqc/4ZPAhMkGKIntqLV4MrF8vxzp1kmnptWpxMCSRsXryBBgzRv7v1qghrTOVK+v/cTQa2ZLGzQ2oVCnjcxQFiIrKOBjduiWzPsPDgadPdW+XPPA6+dK8OfDhh/p/DoaiURRFUbsIYxIdHQ0XFxdERUXB2dlZ7XKIdERGAr/+Kt1ZV64ApUrJp8bevWX2BxEZJ0WRNbKGDJFWlqlTpRvaFMbdREdnHIxu3ZJW5v/+Aw4devXmxjlXZ+bev9niQ2TkFEU+FQYGAr//LttNtG0L/Pij7LPF2RxExu2//2Tw8u7dQIcOwNy5MhjZVDg7A+XKySWtly+BunWB7t2B8+cBF5ccLy/L+CeTyEg9eyazsqpVkybxvXuBsWNlZsf69UCjRgw9RMYsIUHWzvLxkbE8W7cCGzaYVuh5k1y5gFWrZOX3gQPVriZz+GeTyMgEBwP+/jKI8LPP5Ov27fKpcexYmXlBRMbt+HEZePz11/L/+fJloFUrtasyjOLFpUV69Wrgt9/UrubN2NVFZARevJBWnMBA+YNZsCAweDDQrx9QtKja1RFRZkVGyuDlwECZjn7qFFClitpVGV63bsCuXdLq4+sLlCypdkWvxhYfIhVdvw4MGyZN3716yeKCGzcCN28C337L0ENkKpIHL5ctK10/P/4oH2IsIfQkmz8fcHeX8T4JCWpX82oMPkQ5LD4eWLdOdkIvW1b2y+rXT3ZK37MHaN9eFh8kItMxcCDw0UcyHu/KFSAgwDRmbOmTs7OEvtOnZVFGY8XgQ5RDQkKAr76SBcG6dJHZEKtWyZTQGTNkajoRmaZy5eQDy5UrMp7HUr3/PjBxIjB5MnDwoNrVZIzBh8iAXr6U7SNatJA+74ULJfRcuiTrXnTrJrumE5FpCwgALlwAihSRBf06dZI9tSzRl18CderI6vFPnqhdTXoMPkQGcPu2fOopVkzW3Hn8WDYNvXNH+v59fNSukIj0zdsb2LdPZjYdPixd2bNmyQcgS2JtLQutPnsG9O8v45+MCYMPkZ4oiixQ1q6dDEr+/nvZcPDMGeDvv4E+fQBHR7WrJCJD0mhkcO+1azJhYcQIoGpV2ZvLkhQpIuuQbdggW+kYEwYfIj3ZuFGauDdvlt2TmzUDSpSQvW6uXzfuWQ5EpF9588osp1OnZFPPWrVkEsOjR2pXlnM6dgQ++US6Af/3P7WrScG9utLgXl2UXfHxwI4dskLrtWsSdq5dk31uAFnhtGRJaf5OfSlTBnB1Vbd2IjKcxERgyRJZ3ydXLmD6dGkBtoSV1589k4UcnZ2l1cuQO7tn9v2bwScNBh/SJ0UBIiIkACVfkgPRzZsp53l4pA9EZctKc7GlTYklMlf37gEjR8r4Fz8/YNEioGJFtasyvDNnZFHDoUMl9BkKg082MfhQTnn+XJp/U4eh5O9fvJBz7O2B0qXTB6J335XFDonI9Bw8CAwYIP//v/hCJkLkyaN2VYb1/ffAqFGy52CjRoZ5DAafbGLwIbUlJcm4oLSB6No14O7dlPMKF864lahQIRlgSUTGKz4emD1bQo+rKzBnjoyJMdf/u0lJQNOmMhTgxg3AwUH/j8Hgk00MPmTMoqIkDKUNRP/8kzJ42slJxg2lDUSlSnHNICJjc/Om7Mu3ZYtMiJg/33wXM92/X1asv35dWq31LbPv39yklMiEuLjIkvg1augef/kSCA3VDUPXrsmmgcmzSKysZBfljEJR/vzm+0mTyJgVLSozQbdtk9lP5cvLIOjRo83vg0pyd97z5+rWwRafNNjiQ+bm4cP0LUTXrgH//SfNzwBQoQIQHKxunUSW7vlz4LvvgJkzJRAtWCDdQ+bi8mUJdkePyuBufWOLDxEBkNac/PllHZHU4uKAf/+VT5cXLqhTGxGlcHQEpkwBevQA/P2l66tTJxkL9M47alf39pIXcFW7xccCVhEgoozY2cnGisnrbBCRcUi99cWhQ9Id/d13sg3G48dqV5d9xhJ82OJDZMEUBTh7Fhg+XO1KiCi15K0vWrYEvv4amDQpZQJDoULSZZT6Uq4ckDu3ujW/CYMPEakuNBSIjGSLD5GxSt76YtYsWffn0qWUy5Yt0g2WPFK3RAnZADl1ICpTRlp3jUHyFPbkdcrUwuBDZMHOnpWvDD5Exs3WNiXMpBYTI2vjpA5EK1cCt2/L9dbWMnU8bQtRyZI5vyp8rlzyPNjiQ0SqOXsWKFgQ8PRUuxIiyo7cuYFq1eSSWmSkzKJKHYh+/DFleQt7exlLlDYQeXkZdmkLR0cGHyJS0blzbO0hMkd588pMztSzORUFuH9fNwxdugRs2iSTHABZaydtGCpfXvYT1AcGHyJSjaLI5oGffaZ2JUSUEzQaoEABuaTeLyt5m5zUYejUKekyi4uTc9zd04chHx9ZVDUrGHyISDV378qnP7b4EFk2KytZMLFoUZlFluzlS1nrK3UgCgqShRWTFz/18kofiLy9X70Xl4MDgw8RqYQDm4nodXLlkllhZcoAHTqkHI+NldXgUwei9etlB3ZAWpZKlUofiEqXlhYfzuoiIlWcPQvkyyef2IiIYl/GYv3l9dh8fTMePX8EN0c3tC3TFp18OsE+V8rGYfb2QKVKcknt6VPgyhXdQPTTT0BEhFxvYyNrESW3FqnFLFduXrBgAYoVKwZ7e3vUrFkTJ0+eVLskIqNz9qy09nBzUiLaen0rCv1QCD0398Tma5tx8OZBbL62GT0390ShHwph2/Vtb7yPPHmAmjWBFi2A99+XlqLUg6KTF2BMG5hymtkFn7Vr12LYsGEYP348zp49i0qVKqFZs2a4f/++2qURGRXO6CIiQEJP2zVtERkbCQBIUpJ0vkbGRqLNmjbYen1rutsqiiysuHQp0Lu3rA/0zjvARx8Be/bINPtly4ADB2RQde3aMkZITWa3O3vNmjVRvXp1zJ8/HwCQlJQELy8vBAQE4Msvv3zj7bk7O1mChw9llsaaNUCXLmpXQ0RqiX0Zi0I/FEJkbCQUvDoOaKBBXvu8CB9yB9cv2+PwYdk77MgR4N49GSBdqRJQp45catdOWR/s+XOgbl2ZTHH6tP6mxqdlkbuzx8fH48yZMxgzZoz2mJWVFRo3bozjx4+rWBmRcTl3Tr6yxYfIsq2/vB5PYp+88TwFCp7EPkH+uhsQe/Jj2NkBNWoAfftK0PHzAzLKGooi51y9Chw9arjQkxVmFXwePnyIxMREFChQQOd4gQIFcO3atQxvExcXh7jkhQogiZHI3J09K/3xJUuqXQkRqWnz9c2w0lhpu7VeS7FCyQ83IfCHj1G9eub2AJs6VVqW168HKld+63L1wuzG+GTV1KlT4eLior14cYoLWYCzZ4EqVaR5mogsT0QE8PvvwIG/H2Uu9ACAJgnuXo9Ru3bmQs+WLcDYscD48UDHjm9Xrz6Z1Z+9/Pnzw9raGvfu3dM5fu/ePXi+YjOiMWPGICoqSnsJDw/PiVKJVLNqFbBjh8y+ICLL8OgRsHEj0LSpzOQsWBDo1g14fMsNSMpcFLDSWCGfY75MnXvxIvDxx0D79sC4cW9Tuf6ZVfCxtbVF1apVsW/fPu2xpKQk7Nu3D76+vhnexs7ODs7OzjoXInMUGSl/6D7+GGjbFvj6a7UrIiJDiYoCNm8G6tWToJM/v7S6BAWlnNOlCzCoSVvAKnMtPklKEtqVbffG8x4+BNq0AUqUAFasML6WZbMa4wMAw4YNQ69evVCtWjXUqFEDc+bMQUxMDPr06aN2aUSqOXQI6NFD/hiuXg107ap2RUSkTzExwN69wPTpwKvm8vToIR9+Ug9Ejn3ZCat+GJzpWV0dy72+zyohAejUSTY93b8fcHLK7jMyHLMLPl26dMGDBw8wbtw4REREoHLlyti1a1e6Ac9EliA+XvrXp0+XmRcrV8p+PERk2mJjgV27gO++k82GM9KvH9CrF147ENk+lz1WtF2BNmvaQANNhuFHA1nldEXbFTorOGdkyBCZ4v7XX8b7t8bs1vF5W1zHh8zF9etA9+7AhQvApEnAqFGAtbXaVRFRdiQkyGDhSZNk/ExGBg4EPv0UqFAh6//Xt17fit6be+NJ7BPtLK/kr672rljRdgValWn12vsIDAQGDACWLJE6cppFruNDRLJuxpIlwNChsg/X8eOyeioRmY6XL4F16yToXL+e8TlDh0rQKFXq7beeaV2mNe4Mv4MNVzZg07VNePz8MfI55kO7su3QsVzHN7b0HDgABAQAgwapE3qygi0+abDFh0zd1q0ysLBwYVkllb28RKZBUWQczqpVGV8/erSEi3feydm63iQkRLrTKleW7rdcKjWpZPb928jGWhPR22rUSPrZ79+XJeQXLJCxPkRk3OLjdUPPN98ADx5IIFIUYNo04ws9T58CrVsDefMCa9eqF3qygsGHyMzkzg3Mni0bBzZvLp8Qvb1lNldSJtcpI6KcZ2eXEnIURbq58udXu6pXS0qSFqqbN2X8kZub2hVlDoMPkZkqWhRYvhwIDgbKl5eBzu+9B/z5p/xRJSJ6G+PHS9f66tWAj4/a1WQegw+RmStfXj6NHTki+3N98AHQoAFw4oTalRGRqVq7VqbST50KfPih2tVkDYMPkYWoVUsWMty+HXj8GPD1Bdq1A65cUbsyIjIlZ84AffrIYoijRqldTdYx+BBZEI0GaNkSOHcO+PVX4Px5WfPjk0+AsDC1qyMiYxcRIVve+PgAP//89tPo1cDgQ2SBrK1lz65r14A5c6QV6N13geHDZTNDIqK04uJk09HERNkHzMFB7Yqyh+v4pMF1fMgSPX0qM8G+/142FBw1SqbE586tdmVElNMURbrDQ0NljZ6QEPn+1CmZLHHwIFCzptpVppfZ928GnzQYfMiSPXgATJ4MLFoEuLoC48bJfj+2tmpXRkT69PSpbqhJ+/3TpynnOjkBxYvLpX9/6S43Rgw+2cTgQyR//MaPl3FAxYvL7I0uXaQ1iIiM34sXKSEmo2Dz+HHKufb2QLFickkOOKm/z5fPNMbyMPhkE4MPUYpLl4CvvgK2bZPl6KdOBZo1M40/gkTmLD5eJiRkFGpCQoB791LOzZVL1vV6VbDx8DCPDzUMPtnE4EOU3tGjwJdfylpA9erJ0vnvv692VUTmKzERuHUr42ATGgrcvp2yEruVlezNlxxm0gacd97J+m7tpoi7sxOR3iSvAbRjBzBmTMoaQJMny3YYRJQ1SUnSKvOqcTZhYbJDezJPz5QwU6eObrDx8uI4vKxgi08abPEher3ERFmi/ptvgPBwoHdvYMIE+eNLRBk7fx5YuhS4cUPCzc2bQGxsyvVubhm31hQvLt1Upjp1PCexqyubGHyIMicuDli8WAY+R0cDgwZJa5CpbFRIlBNOngS+/VbWyvLykv3y0gabYsVkOxl6Oww+2cTgQ5Q1T58Cs2YBM2dyDSCiZIcPy4eCPXuAMmWAsWOBrl1loDEZRmbfv81gHDcRqSlPHpn6/u+/0u01cSJQsiSwcCGQkKB2dUQ5R1GAffuA+vWBunWBu3dlM8/Ll4EePRh6jAWDDxHphYcHMHcu8L//AU2bSteXtzewZk3K7BMic6QowM6dgJ8f0Lgx8OyZbOlw/jzQubNlzKgyJQw+RKRXxYoBK1cCFy5I8OnaFahWDdi9W94giMxFUhKwaZP8frdsKV29f/4pWzu0aWMea+OYI/6zEJFBVKggCx8ePgw4OgLNmwMNGwJ//612ZURvJzFRWjIrVZJNO11cpIvryBH5PecCn8aNwYeIDKp2bQk/W7cCDx/Kwoft28vO8ESmJCEBWLECKFdOWjILF5aw89dfEuoZeEwDgw8RGZxGA7RqJWMeVqwAzp4FfHxkA9Rbt9Sujuj14uKAJUtkdlbv3tKFe/KkdGvVqqV2dZRVDD5ElGOsrYGePYHr12UK/JYtQKlSwMiRupsmEhmDFy+A+fPld/Tzz2Usz/nzMnC5enW1q6PsYvAhohxnZwcMHixT4L/8EggMBEqUAKZMAWJi1K6OLN2zZ8APP8jv5ODBMj398mVg3ToZ10OmjcGHiFTj7CzbXfz7L9Crl3zv7c3xP6SO6GgJ38WKSSBv2VJaJ3/9lXvSmRMGHyJSXfIaQNeuSRhK/oRNlBMeP5ZFOIsWlQU4O3eWPbV+/lm6uci8MPgQkdEoUQLYvx8oUEDCT3Cw2hWRObt/X1p2ihYFvv9eBi6HhMiq40WLql0dGQqDDxEZFXd3mR7s5QU0aACcO6d2RWRu7twBhg2TLq0FCwB/fyA0FJg9GyhUSO3qyNAYfIjI6Li5yYJwJUvK+iinT6tdEZmDmzcl5JQoAfzyi8wmDA0Fpk2T7layDAw+RGSUXF2BoCCgbFnZ/4grPlN2/fuvrBlVqpRsGjpunISgiRMlZJNlYfAhIqPl4iJ7fJUvDzRpAhw9qnZFZEquXpVd0d99F9i+HZg6VVp4vvpKfrfIMjH4EJFRc3YGdu0C3nsPaNYMOHRI7YrI2AUHy8wsHx8ZLD9njgxaHjECcHJSuzpSG4MPERk9Jydgxw7Z56t5cxn8TJTW6dNA27ayyOCpU7Iw5r//AgEBgIOD2tWRsWDwISKTkDu37PZet64sLLdnj9oVkbE4dgxo0UK2kbh6FVi+HPjf/4D+/WWVcKLUGHyIyGQ4OMg+SY0aAa1bAzt3ql0RqUVRpBurYUPZKDQ8HPj9d+DKFVkF3MZG7QrJWDH4EJFJsbcHNm6ULq927aQViCyDogBhYTIzq04dCT2RkfL7EBwMfPSRbIRL9Dq51C6AiCir7Oxkw8hu3YD27eX7du3Uror07c4dGbeT+vLggVxXs6bM1PrgA0CjUbdOMi0MPkRkkmxtpWujRw+gUyf5vlMntaui7Lp3L33IiYiQ6zw8gGrVgAED5GvVqlxhmbJPoyiKonYRxiQ6OhouLi6IioqCs7Oz2uUQ0Wsoinzq79AByJtX3jz56d/4PXwInDmjG3Ju3ZLr3Nwk3CRfqlYFChfmvyu9WWbfv9niQ0Qm6fRpYNQoGeBarx4wcybfHI3RkyfA2bO6ISc0VK5zcZFw0717StApWpT/jmRYDD5EZFJCQoCxY6Vrq1w5GdzcsiXfLI1BdHT6kPPvv3Kdk5O03nTokBJySpbkvxvlPAYfIjIJjx4BkycD8+cD+fMDP/0E9O4N5OJfMVU8ewacP68bcq5fl+scHWWl7VatUrqr3n0XsOI8YjIC/JNBREbtxQvgxx9ln6WkJGD8eGDIEFnQkHLGixfAhQu6IefqVfn3sLcHKleWvdTGjJGgU7Ysp5WT8WLwISKjlJgI/PYb8PXXMrvn88+Bb76RGT5kOHFxsiZOcsA5cwa4dEn+PWxsZDuIOnWAoUMl5JQrx8UCybQw+BCRUVEU2Y5i1Ch5A+7USbq4SpdWuzLzk5AgoSZ1S87Fi3I8Vy6gQgXZBiJ5Gnn58twCgkwfgw8RGY1z5yTw7N0L1K4NHD8uG5PS23v5UrqnUoecCxekhcfKSnYyr1YN6NtXvlasKN1YROaGwYeIVHfzpnRp/fabjA/ZskUGxnLGT/YkJspA49Qh5/x5Gauj0chrnHoaeaVKHDNFloPBh4hU8+QJMGWKDF7Olw9YvBj45BPzmamlKNLSEh8vl4SElO9fd+xtzg0NlSnlMTFSQ+nSEm46dpSvVaoAefKo+rIQqcpM/rwQkSmJjQUWLJCxOwkJsi7PsGGy1ktOi4gAfv5ZWkMMEVD0wdpatuiwtZWBxMnfp74kH/f0lJlvySEnb1791EBkLhh8iCjHJCUBq1dL0Ll9G/jsM2DcOKBAAfVqunRJalAU2f/J1TXjQGFrCzg7Zz6A6Ou4jQ2nhhPpE4MPEeWIvXtl4PK5c7Kj+p49QJkyalcFNG4sqz937y5dQOvXA97ealdFRIbCdTSJyKAuXACaN5cF7hwcgKNHgY0bjSP0JGvZEjh1SsYW1agBbNqkdkVEZCgmE3wmT54MPz8/ODo6Iu8rOq3DwsLQsmVLODo6wsPDAyNHjsTLly9ztlAiAgCEh8uWElWqyP5af/wBHDkC+PmpXVnGSpcGTpyQkNa+vXTHJSaqXRUR6ZvJBJ/4+Hh06tQJAwYMyPD6xMREtGzZEvHx8Th27BhWrFiB5cuXY9y4cTlcKZFli4wERo+WIPHnn8C0adJ9VKiQzHAyZk5OwLp1wPTpUveHHwKPH6tdFRHpk0ZRFEXtIrJi+fLlGDJkCCIjI3WO//nnn/jwww9x584dFPj/kZKBgYEYPXo0Hjx4AFtb20zdf3R0NFxcXBAVFQVnZ2d9l09kVh49Ao4dk+6ro0elRSc1a2vdVhM3N2lN6dIFqFfPuKetBwUBH30ks6I2bZIF/YjIeGX2/duI/+xkzfHjx1GhQgVt6AGAZs2aYcCAAbh8+TKqVKmiYnVEpicuTtaDSQ41R48CDx5k7rbVq8ulSBHAy0suNjbA5s3A2rWys7qHh6wt07mzrNJsbDOXmjSRhf/atwd8fYGlSyUIEZFpM5vgExERoRN6AGh/joiIeOXt4uLiEBcXp/05OjraMAUSGRFFAW7c0A01V69m7T4cHIBateRSu7YMCn5TI+n778su66dPSwBatw5YuBAoWFD25OrSRc6xMpJO+OLF5bX57DOga1cZAD19unG3VBHR66n63/fLL7/E9OnTX3vO1atXUbZsWYPVMHXqVEycONFg90+ktsREoEMH2QYiKypXTgk2tWpJq40+tpDQaFJahGbMkAHFa9fKOKAff5THSQ5B1aurv22FoyOwcqXUMmyYTMdfuxZwd1e3LiLKHlXH+Dx48ACPHj167TklSpTQGZ/zqjE+48aNw9atW3H+/HntsZCQEJQoUQJnz559ZVdXRi0+Xl5eHONDZiM2VlpnUvPykjDj5yetNRUqqN+KkZQkY4TWrgU2bADu35cWl86dJQRVrqx+CDp4UOqxs5NxP1WrqlsPEaUwiTE+7u7ucNfTxyZfX19MnjwZ9+/fh4eHBwAgKCgIzs7OKFeu3CtvZ2dnBzs7O73UQGSM7O2la8vYWVkBdevKZe5cCRlr18p2EtOnyyyxLl0keJQvr04IqlcPOHNGWtBq1QICA2XKPhGZDiPpSX+zsLAwnD9/HmFhYUhMTMT58+dx/vx5PHv2DADQtGlTlCtXDj169MCFCxewe/dufP311/D392ewITIxuXIBjRoBS5YAd+8Cu3ZJy9T8+TK7yscHmDgRuHYt52srXFhCWY8eQJ8+gL+//vbkIiLDM5np7L1798aKFSvSHd+/fz/q168PALh58yYGDBiAAwcOIHfu3OjVqxemTZuGXFlow+d0diLjFR8vW12sXStjlp4+lW66xo2BOnWkFeb/G3xzxJIlwKBBMrB7/XoZpE1E6sjs+7fJBJ+cwuBDZBpiY6Ul6I8/gEOHgJs35XiZMtI6VLu2hKESJQzbLXb8uEzLVxTZisPX13CPRUSvxuCTTQw+RKYpPFwGRydfLl6UMOLpmRKCateWrjJ9D+SOiJCZaH//LTPTPvtM/YHYRJaGwSebGHyIzMOTJ9Iac/iwBKGTJ6WrLE8eaZVJbhWqWVOmrL+t+Hhg+HAZh/TJJ8CCBTKwnIhyBoNPNjH4EJmn2FhZODG5RejoUdlXLFcumZae3CpUqxaQP3/2H2fFCmnxqVhRur68vPT2FIjoNRh8sonBh8gyJCUBly+ntAgdPgzcuiXXeXvrjhMqVixrXVdnzshWFy9eyKDnevUM8hSIKBUGn2xi8CGyXDdvprQIHT4swQiQneVTjxOqUOHNe4s9eCB7ex08CPzwA/DFFxz3Q2RIDD7ZxOBDRMkeP07Zdf7IEdmrKyFB9iRLXvW6Th3ZziLt6tgA8PIl8OWXEny6d5fp7/oYT0RE6TH4ZBODDxG9yosXEn6SW4SOHQOio2Xn+WrVUlqEatUC8uVLud2aNUDfvsC778r0++LF1XsOROaKwSebGHyIKLMSE4FLl3THCd25I9f5+OiOE4qKAtq1kwHVv/8ONG2qaulEZofBJ5sYfIgouxQFCA1NCUFHjgBXr8p1hQtLGNq7V86bMgUYNYrjfoj0xSQ2KSUiMicajXRjFS8ue3kBwMOHuuOENBppKfryS9mHbM4cVUsmsjgMPkREBpQ/P9CmjVwA4PlzWUzx8GGgZEl1ayOyRAw+REQ5yNERqF9fLkSU86zULoCIiIgopzD4EBERkcVg8CEiIiKLweBDREREFoPBh4iIiCwGgw8RERFZDAYfIiIishgMPkRERGQxGHyIiIjIYjD4EBERkcVg8CEiIiKLweBDREREFoPBh4iIiCwGgw8RERFZjFxqF2BsFEUBAERHR6tcCREREWVW8vt28vv4qzD4pPH06VMAgJeXl8qVEBERUVY9ffoULi4ur7xeo7wpGlmYpKQk3LlzB3ny5IFGo1G7HKMQHR0NLy8vhIeHw9nZWe1yTAJfs6zh65U1fL2yhq9X1pjq66UoCp4+fYpChQrByurVI3nY4pOGlZUVChcurHYZRsnZ2dmk/hMYA75mWcPXK2v4emUNX6+sMcXX63UtPck4uJmIiIgsBoMPERERWQwGH3ojOzs7jB8/HnZ2dmqXYjL4mmUNX6+s4euVNXy9ssbcXy8ObiYiIiKLwRYfIiIishgMPkRERGQxGHyIiIjIYjD4EBERkcVg8KFXCg0NRd++fVG8eHE4ODigZMmSGD9+POLj43XO0Wg06S4nTpxQsXL1ZOY1A4Dg4GDUqVMH9vb28PLywowZM1SqWH2TJ0+Gn58fHB0dkTdv3gzPyeh3bM2aNTlbqJHIzOsVFhaGli1bwtHRER4eHhg5ciRevnyZs4UaqWLFiqX7XZo2bZraZRmVBQsWoFixYrC3t0fNmjVx8uRJtUvSK67cTK907do1JCUlYfHixShVqhQuXbqETz/9FDExMZg5c6bOuXv37oWPj4/2Zzc3t5wu1yhk5jWLjo5G06ZN0bhxYwQGBuLixYv45JNPkDdvXvTv31/lZ5Dz4uPj0alTJ/j6+mLp0qWvPG/ZsmVo3ry59udXvembuze9XomJiWjZsiU8PT1x7Ngx3L17Fz179oSNjQ2mTJmiQsXGZ9KkSfj000+1P+fJk0fFaozL2rVrMWzYMAQGBqJmzZqYM2cOmjVrhuvXr8PDw0Pt8vRDIcqCGTNmKMWLF9f+HBISogBQzp07p15RRi7ta7Zw4ULF1dVViYuL0x4bPXq0UqZMGTXKMxrLli1TXFxcMrwOgLJp06YcrcfYver12rlzp2JlZaVERERojy1atEhxdnbW+Z2zVEWLFlVmz56tdhlGq0aNGoq/v7/258TERKVQoULK1KlTVaxKv9jVRVkSFRWFfPnypTveunVreHh4oHbt2ti6dasKlRmvtK/Z8ePHUbduXdja2mqPJX+ievLkiRolmgR/f3/kz58fNWrUwC+//AKFS5Bl6Pjx46hQoQIKFCigPdasWTNER0fj8uXLKlZmPKZNmwY3NzdUqVIF33//PbsB/198fDzOnDmDxo0ba49ZWVmhcePGOH78uIqV6Re7uijTbty4gXnz5ul0czk5OeGHH35ArVq1YGVlhY0bN6Jt27bYvHkzWrdurWK1xiGj1ywiIgLFixfXOS/5TSoiIgKurq45WqMpmDRpEho2bAhHR0fs2bMHAwcOxLNnz/DFF1+oXZrRiYiI0Ak9gO7vl6X74osv8N577yFfvnw4duwYxowZg7t372LWrFlql6a6hw8fIjExMcPfn2vXrqlUlf6xxccCffnllxkOFk19SftLfvv2bTRv3hydOnXS6RvPnz8/hg0bhpo1a6J69eqYNm0aPv74Y3z//fc5/bQMSp+vmSXIzuv1Ot988w1q1aqFKlWqYPTo0Rg1apRZ/Y7p+/WyNFl5/YYNG4b69eujYsWK+Pzzz/HDDz9g3rx5iIuLU/lZUE5hi48FGj58OHr37v3ac0qUKKH9/s6dO2jQoAH8/PywZMmSN95/zZo1ERQU9LZlGhV9vmaenp64d++ezrHknz09PfVTsMqy+nplVc2aNfHtt98iLi7OLPYT0ufr5enpmW4Wjrn9fqX1Nq9fzZo18fLlS4SGhqJMmTIGqM505M+fH9bW1hn+fTKn3x0GHwvk7u4Od3f3TJ17+/ZtNGjQAFWrVsWyZctgZfXmRsLz58+jYMGCb1umUdHna+br64uxY8ciISEBNjY2AICgoCCUKVPGbLq5svJ6Zcf58+fh6upqFqEH0O/r5evri8mTJ+P+/fvaWThBQUFwdnZGuXLl9PIYxuZtXr/z58/DysrKfGYsvQVbW1tUrVoV+/btQ9u2bQEASUlJ2LdvHwYNGqRucXrE4EOvdPv2bdSvXx9FixbFzJkz8eDBA+11yel/xYoVsLW1RZUqVQAAf/zxB3755Rf8/PPPqtSstsy8Zt26dcPEiRPRt29fjB49GpcuXcLcuXMxe/ZstcpWVVhYGB4/foywsDAkJibi/PnzAIBSpUrByckJ27Ztw7179/D+++/D3t4eQUFBmDJlCkaMGKFu4Sp50+vVtGlTlCtXDj169MCMGTMQERGBr7/+Gv7+/mYTFLPr+PHj+Pvvv9GgQQPkyZMHx48fx9ChQ/Hxxx+bzYeOtzVs2DD06tUL1apVQ40aNTBnzhzExMSgT58+apemP2pPKyPjtWzZMgVAhpdky5cvV7y9vRVHR0fF2dlZqVGjhrJ+/XoVq1ZXZl4zRVGUCxcuKLVr11bs7OyUd955R5k2bZpKFauvV69eGb5e+/fvVxRFUf7880+lcuXKipOTk5I7d26lUqVKSmBgoJKYmKhu4Sp50+ulKIoSGhqqtGjRQnFwcFDy58+vDB8+XElISFCvaCNx5swZpWbNmoqLi4tib2+veHt7K1OmTFFiY2PVLs2ozJs3TylSpIhia2ur1KhRQzlx4oTaJemVRlE4J5SIiIgsA2d1ERERkcVg8CEiIiKLweBDREREFoPBh4iIiCwGgw8RERFZDAYfIiIishgMPkRERGQxGHyIyKItX74cefPmfeN5Go0GmzdvNng9RGRYDD5ElCMSExPh5+eH9u3b6xyPioqCl5cXxo4d+8rb1q9fX7vLtr29PcqVK4eFCxfqpa4uXbrgf//7n/bnCRMmoHLlyunOu3v3Llq0aKGXxyQi9TD4EFGOsLa2xvLly7Fr1y6sWrVKezwgIAD58uXD+PHjX3v7Tz/9FHfv3sWVK1fQuXNn+Pv74/fff3/ruhwcHDK1QaWnp6fF73VFZA4YfIgox7z77ruYNm0aAgICcPfuXWzZsgVr1qzBypUrYWtr+9rbOjo6wtPTEyVKlMCECRNQunRpbN26FYBs3NmmTRs4OTnB2dkZnTt3xr1797S3vXDhgnZjSmdnZ1StWhWnT58GoNvVtXz5ckycOBEXLlzQtjAtX74cQPqurosXL6Jhw4ZwcHCAm5sb+vfvj2fPnmmv7927N9q2bYuZM2eiYMGCcHNzg7+/PxISEvTwShJRdnF3diLKUQEBAdi0aRN69OiBixcvYty4cahUqVKW78fBwQHx8fFISkrShp6DBw/i5cuX8Pf3R5cuXXDgwAEAQPfu3VGlShUsWrQI1tbWOH/+PGxsbNLdZ5cuXXDp0iXs2rULe/fuBQC4uLikOy8mJgbNmjWDr68vTp06hfv376Nfv34YNGiQNigBwP79+1GwYEHs378fN27cQJcuXVC5cmV8+umnWX6+RKQfDD5ElKM0Gg0WLVoEb29vVKhQAV9++WWWbp+YmIjff/8dwcHB6N+/P/bt24eLFy8iJCQEXl5eAICVK1fCx8cHp06dQvXq1REWFoaRI0eibNmyAIDSpUtneN8ODg5wcnJCrly54Onp+coaVq9ejdjYWKxcuRK5c+cGAMyfPx+tWrXC9OnTUaBAAQCAq6sr5s+fD2tra5QtWxYtW7bEvn37GHyIVMSuLiLKcb/88gscHR0REhKCW7duZeo2CxcuhJOTExwcHPDpp59i6NChGDBgAK5evQovLy9t6AGAcuXKIW/evLh69SoAYNiwYejXrx8aN26MadOm4d9//32r+q9evYpKlSppQw8A1KpVC0lJSbh+/br2mI+PD6ytrbU/FyxYEPfv33+rxyait8PgQ0Q56tixY5g9eza2b9+OGjVqoG/fvlAU5Y236969O86fP4+QkBDExMRg1qxZsLLK3J+wCRMm4PLly2jZsiX++usvlCtXDps2bXrbp/JGabvTNBoNkpKSDP64RPRqDD5ElGOeP3+O3r17Y8CAAWjQoAGWLl2KkydPIjAw8I23dXFxQalSpfDOO+/oBB5vb2+Eh4cjPDxce+zKlSuIjIxEuXLltMfeffddDB06FHv27EH79u2xbNmyDB/H1tYWiYmJr63F29sbFy5cQExMjPbY0aNHYWVlhTJlyrzxuRCRehh8iCjHjBkzBoqiYNq0aQCAYsWKYebMmRg1ahRCQ0OzdZ+NGzdGhQoV0L17d5w9exYnT55Ez549Ua9ePVSrVg0vXrzAoEGDcODAAdy8eRNHjx7FqVOn4O3tneH9FStWDCEhITh//jwePnyIuLi4dOd0794d9vb26NWrFy5duoT9+/cjICAAPXr00I7vISLjxOBDRDni4MGDWLBgAZYtWwZHR0ft8c8++wx+fn6Z7vJKS6PRYMuWLXB1dUXdunXRuHFjlChRAmvXrgUg6wc9evQIPXv2xLvvvovOnTujRYsWmDhxYob316FDBzRv3hwNGjSAu7t7hmsFOTo6Yvfu3Xj8+DGqV6+Ojh07olGjRpg/f36W6yeinKVRsvOXhoiIiMgEscWHiIiILAaDDxEREVkMBh8iIiKyGAw+REREZDEYfIiIiMhiMPgQERGRxWDwISIiIovB4ENEREQWg8GHiIiILAaDDxEREVkMBh8iIiKyGAw+REREZDH+DyPgqf4H8NeWAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "import random\n", + "import matplotlib.pyplot as plt\n", + "import math\n", + "\n", + "f=0.1 #facteur de corrélation\n", + "\n", + "k_B = 1.38e-23 \n", + "T = 300 \n", + "M_w = 18e-3 \n", + "N = 6.022e23 \n", + "\n", + "\n", + "v = math.sqrt(8 * k_B * T / (math.pi * M_w / N))\n", + "\n", + "dt = 0.1 \n", + "n = 50 \n", + "\n", + "step_size = v * dt * f / math.sqrt(2)\n", + "\n", + "x = 0\n", + "y = 0\n", + "\n", + "path_x = [x]\n", + "path_y = [y]\n", + "\n", + "for i in range(n):\n", + " dx = random.uniform(-1, 1) * step_size\n", + " dy = random.uniform(-1, 1) * step_size\n", + " x += dx\n", + " y += dy\n", + " path_x.append(x)\n", + " path_y.append(y)\n", + "\n", + "\n", + "distance = ((x - 0) ** 2 + (y - 0) ** 2) ** 0.5\n", + "\n", + "\n", + "print(f\"La distance entre le point de départ et le point d arrivée est {distance:.2f}.\")\n", + "print(f\"La vitesse de la particule est {v:.2f} m/s.\")\n", + "\n", + "\n", + "plt.plot(path_x, path_y, 'b-', linewidth=1.0)\n", + "plt.plot(0, 0, 'go', markersize=8.0)\n", + "plt.plot(x, y, 'ro', markersize=8.0)\n", + "plt.xlabel('X Position')\n", + "plt.ylabel('Y Position')\n", + "plt.title(' Simulation Random Walk')\n", + "plt.show()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "d32f5365-a3ba-467e-9c8c-67e6d225cb20", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "e6000f8c-8713-45d2-81d2-7da227778ca2", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.8.10" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}