diff --git a/proto/FTCS.ipynb b/proto/FTCS.ipynb index f47d159..b3febbf 100644 --- a/proto/FTCS.ipynb +++ b/proto/FTCS.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -17,7 +17,7 @@ }, { "cell_type": "code", - "execution_count": 113, + "execution_count": 2, "metadata": {}, "outputs": [], "source": [ @@ -42,23 +42,23 @@ "max_stable_time_step_x = delta_x**2 / (2 * np.max(alpha_x))\n", "max_stable_time_step_y = delta_y**2 / (2 * np.max(alpha_y))\n", "max_stable_time_step = max_stable_time_step_x if max_stable_time_step_x <= max_stable_time_step_y else max_stable_time_step_y\n", - "max_stable_time_step = 0.25\n", + "max_stable_time_step = 0.1\n", "\n", "C_t = np.zeros((grid_size['x'],grid_size['y']))" ] }, { "cell_type": "code", - "execution_count": 114, + "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 114, + "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, @@ -83,7 +83,7 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 4, "metadata": {}, "outputs": [], "source": [ @@ -96,7 +96,7 @@ }, { "cell_type": "code", - "execution_count": 51, + "execution_count": 5, "metadata": {}, "outputs": [], "source": [ @@ -152,12 +152,12 @@ }, { "cell_type": "code", - "execution_count": 115, + "execution_count": 38, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAT4AAAD8CAYAAADub8g7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAAsTAAALEwEAmpwYAAAYa0lEQVR4nO3dfawd1X3u8e8TY0BxoOC4cXgriVoLiaLitJZpFNqaEIixuCGpuKndKoWEyiEKUpBa9dJWgoj8k6ol6PaSG18HLENFHO5N4gQ1DuDSSA5SQjCWSWxeioOIsHFsGRNe8gL4nOf+MXPSne29zxnPnnPO3meeTzQ6M2vWzFobKz+tl5lZsk1ERJu8abYrEBEx0xL4IqJ1EvgionUS+CKidRL4IqJ1EvgionUS+CJiVkk6S9K3JT0uabekT5XpCyVtlfR0+ffUPtdfVeZ5WtJVlcrMc3wRMZsknQacZnuHpJOAR4EPAlcDh21/VtINwKm2/0fXtQuB7cAywOW1f2D7xcnKTIsvImaV7f22d5T7rwBPAGcAVwB3ltnupAiG3d4PbLV9uAx2W4GVU5V5XAP1btzxOsEnsmC2qxExZ/2Sn/G6X9Mg93j/RQv8wuGxSnkf/cFru4FfdiStt72+O5+kdwDvAh4GFtveX576CbC4x63PAJ7rON5bpk1qKAPfiSzgAl0829WImLMe9oMD3+PQ4TEevv/MSnnnn/ajX9peNlkeSW8Bvgpcb/tl6b/ism1LamxcLl3diKjJjHm80jYVSfMpgt7dtr9WJh8ox/8mxgEP9rh0H3BWx/GZZdqkBgp8klZKekrSnnLwsfv8CZLuKc8/XDZjI2IOMDCOK22TUdG0uwN4wvbnOk7dC0zM0l4FfKPH5fcDl0o6tZz1vbRMm1TtwCdpHvB54DLgXGCNpHO7sl0DvGj7d4BbgX+sW15EDJ/xiv+bwnuAjwDvlbSz3FYBnwUukfQ08L7yGEnLJN0OYPsw8BngkXK7uUyb1CBjfMuBPbafKSvzZYpZmMc78lwBfLrc/wpwmyQ5z9BEjDxj3qjQjZ3yPvZDQL+JlqMG+21vB/6q43gDsOFYyhykq1tlNuVXeWwfAV4C3trrZpLWStouafsbvDZAtSJiJhgYw5W2YTM0s7rl1PZ6gJO1cPj+S0XEUaYavxtWgwS+KrMpE3n2SjoO+A3ghQHKjIghYWBsREetBunqPgIskfROSccDqylmYTp1zspcCfxHxvci5o7xituwqd3is31E0nUUU8fzgA22d0u6Gdhu+16KKep/lbQHOEwRHCNiDvCQjt9VMdAYn+0twJautBs79n8J/PdByoiI4WTDG6MZ94ZnciMiRo0Y6/sUynBL4IuIWgyMp8UXEW2TFl9EtErxAHMCX0S0iIE3PJofeErgi4hajBgb0S/bJfBFRG3jTlc3IlokY3wR0UJiLGN8EdEmxReYE/giokVs8brnzXY1akngi4jaxjPGFxFtUkxuNNPVlbQBuBw4aPu8Mu0e4JwyyynAT20v7XHts8ArwBhwZKplLCGBLyJqa3RyYyNwG3DXRILtP/tVSdItFEtX9HOR7UNVC0vgi4hampzcsL2t3/Kz5fKTHwbe20hhZEHxiBjAmFVpG9AfAQdsP93nvIEHJD0qaW2VG6bFFxG1GPGGK4eQRZK2dxyvLxcYq2INsGmS8xfa3ifpbcBWSU/a3jbZDRP4IqKWY5zcOFRl0qFbuUjZnwJ/0Lce9r7y70FJmynW/J408NXu6ko6S9K3JT0uabekT/XIs0LSSx2ro9/Y614RMXpMtW7ugF3d9wFP2t7b66SkBZJOmtgHLgV2TXXTQVp8R4C/tr2jLPhRSVttP96V7zu2Lx+gnIgYUk1NbkjaBKyg6BLvBW6yfQfFAmWbuvKeDtxuexWwGNhczH9wHPAl2/dNVd4gq6ztB/aX+69IegI4A+gOfBExB9k09jiL7TV90q/ukfY8sKrcfwY4/1jLa2SMr5yGfhfwcI/T75b0GPA88De2d/e5x1pgLcCJvLmJakXENComN1r6ypqktwBfBa63/XLX6R3A2bZflbQK+DqwpNd9yhme9QAna+GILmES0S6j+iHSgWotaT5F0Lvb9te6z9t+2far5f4WYL6kRYOUGRHDwYhxV9uGTe0WX/k09R3AE7Y/1yfP2ykePLSk5RSB9oW6ZUbEcBnVFt8gXd33AB8BfihpZ5n298BvAdheB1wJfELSEeAXwGrb6cZGzAHFurotC3y2H4LJv0lj+zaKF48jYs5RPj0fEe1SLC/Z0lndiGgnW+3r6kZEZLGhiGiV4nt8GeOLiFbJ8pIR0TLF4yxp8UVEi7T6Xd2IaK8sKB4RrVJ8lipd3YhomYzxRUSrFF9nSVc3IlqkeGUtgS8iWmV0W3yjWeuIGArjqNI2FUkbJB2UtKsj7dOS9nWs0riqz7UrJT0laY+kG6rUO4EvImqZmNVtaHnJjcDKHum32l5ablu6T0qaB3weuAw4F1gj6dypCkvgi4jaxv2mSttUbG8DDteownJgj+1nbL8OfBm4YqqLEvgiopZjXHNjkaTtHdvaisVcJ+kHZVf41B7nzwCe6zjeW6ZNKpMbEVGLgSPVJzcO2V52jEV8AfhMWdRngFuAjx3jPXpqYnnJZ4FXgDHgSPePKxcl+p8UCwD/HLja9o5By42I2Teds7q2D0zsS/oi8G89su0Dzuo4PrNMm1RTLb6LbB/qc+4yirV0lwAXUETxCxoqNyJmyzQvHSnpNNv7y8MPAbt6ZHsEWCLpnRQBbzXw51Pdeya6ulcAd5Wrq31P0ildPygiRlCTHyKVtAlYQTEWuBe4CVghaWlZ1LPAx8u8pwO3215l+4ik64D7gXnABtu7pyqvicBn4AFJBv6P7fVd5/sNPv5a4CsHO9cCnMibG6hWREy3plp8ttf0SL6jT97nKYbOJo63AEc96jKZJgLfhbb3SXobsFXSk+XU9DEpA+Z6gJO1MGvvRgy5Vn+I1Pa+8u9BSZspnqvpDHy1Bh8jYrgZcWR8NJ+IG6jWkhZIOmliH7iUowcg7wX+UoU/BF7K+F7E3NDUK2szbdAW32Jgc/HECscBX7J9n6RrAWyvo+h7rwL2UDzO8tEBy4yIYeCWdnVtPwOc3yN9Xce+gU8OUk5EDJ9Wj/FFRHsl8EVEqxgxNqKTGwl8EVHbME5cVJHAFxG1uK2TGxHRbk7gi4h2md6PFEynBL6IqC0tvohoFRvGxhP4IqJlMqsbEa1i0tWNiNbJ5EZEtJBH9MuZCXwRUVu6uhHRKsWsbjPv6kraAFwOHLR9Xpn2T8B/A14HfgR81PZPe1z7LJOs9NjLaL5hHBFDwa62VbARWNmVthU4z/bvAf8J/N0k119ke2nVtXsT+CKiNluVtqnv423A4a60B2wfKQ+/R7FsRSMS+CKiFlMt6JWBb5Gk7R3b2mMs7mPAt/pWpVjp8dGq980YX0TUdgyTuoeqdkO7SfoH4Ahwd58sx7zSY+0Wn6RzJO3s2F6WdH1XnhWSXurIc2Pd8iJiyBg8rkpbXZKuppj0+ItyGYujq9Gx0iMwsdLjpGq3+Gw/BSwtKzePYsnIzT2yfsf25XXLiYjhNZ2Ps0haCfwt8Ce2f94nzwLgTbZf6Vjp8eap7t3UGN/FwI9s/7ih+0XECGhqVlfSJuC7wDmS9kq6BrgNOImi+7pT0roy7+mStpSXLgYekvQY8H3gm7bvm6q8psb4VgOb+px7d1mp54G/sb27V6ZyUHItwIm8uaFqRcR0afJdXdtreiTf0Sfv8xRL1vZd6XEqA7f4JB0PfAD4fz1O7wDOtn0+8L+Ar/e7j+31tpfZXjafEwatVkRMNwNWtW3INNHVvQzYYftA9wnbL9t+tdzfAsyXtKiBMiNiCDT4APOMaqKru4Y+3VxJbwcO2Lak5RSB9oUGyoyIWTfYjO1sGijwlbMolwAf70i7FsD2OuBK4BOSjgC/AFb3m5KOiBE0ov9vHijw2f4Z8NautHUd+7dRzMxExFzjfJ0lItqojS2+iGi7tPgiom3GZ7sC9STwRUQ9E8/xjaAEvoiobVSf0Ujgi4j6EvgionXS1Y2ItlFafBHRKha08ZW1iGi5tPgionUS+CKidRL4IqJVRvgB5qyrGxG1ydW2Ke8jbZB0UNKujrSFkrZKerr8e2qfa68q8zwt6aoq9U7gi4j6XHGb2kZgZVfaDcCDtpcAD5bHv0bSQuAm4AKKZSVv6hcgOyXwRURtTbX4ygXAD3clXwHcWe7fCXywx6XvB7baPmz7RWArRwfQo2SMLyLqqz7Gt0jS9o7j9bbXT3HNYtv7y/2fUCwl2e0M4LmO471l2qQS+CKinurdWIBDtpfVLqpYt6exOeRKXd2ZHniMiBHR3BhfLwcknQZQ/j3YI88+4KyO4zPLtElVHePbyAwOPEbEaNB4ta2me4GJxtJVwDd65LkfuFTSqWVsubRMm1SlwDfTA48RMSIaavFJ2gR8FzhH0l5J1wCfBS6R9DTwvvIYScsk3Q5g+zDwGeCRcru5TJvUIGN80zbwGBHDr+qMbRW21/Q5dXGPvNuBv+o43gBsOJbyGpncaGLgUdJaYC3Aiby5iWpFxHRr4ZsbjQ482l5ve5ntZfM5YYBqRcSMmd7JjWkzSOCbtoHHiBgNTT3APNOqPs4yowOPETECPO2zutOm0hjfTA88RsSIGMLWXBV5cyMi6kvgi4i2GcbxuyrydZaIaJ20+CKivhFt8SXwRUQ9Hs4Z2yoS+CKivrT4IqJNxOhObiTwRUR9CXwR0SpD+jpaFQl8EVFfJjciom3S4ouI9kngi4hWGdJv7VWRV9YiorYmvscn6RxJOzu2lyVd35VnhaSXOvLcOEi90+KLiPoaaPHZfgpYCiBpHsVX2jf3yPod25cPXmICX0QMYBpeWbsY+JHtHzd+5w7p6kZEPVXX2yhahYskbe/Y1va562pgU59z75b0mKRvSfrdQaqeFl9E1KJyq+iQ7WWT3k86HvgA8Hc9Tu8Azrb9qqRVwNeBJdWL/3Vp8UVEfc2usnYZsMP2gaOKsV+2/Wq5vwWYL2lR3WpPGfgkbZB0UNKujrR/kvSkpB9I2izplD7XPivph+UszPa6lYyI4dTwKmtr6NPNlfR2SSr3l1PErhfq1rtKi28jsLIrbStwnu3fA/6T3k3TCRfZXjpVMzciRlBDLT5JC4BLgK91pF0r6dry8Epgl6THgH8BVtuuPac85Rif7W2S3tGV9kDH4ffKSkVEmzT4IVLbPwPe2pW2rmP/NuC2ZkprZozvY8C3+pwz8ICkRyeZxQFA0tqJGZ83eK2BakXEtGt2jG/GDDSrK+kfgCPA3X2yXGh7n6S3AVslPWl7W6+MttcD6wFO1sIh/E8VEd1G9SMFtVt8kq4GLgf+ol9f2/a+8u9Biiexl9ctLyKG0Ii2+GoFPkkrgb8FPmD7533yLJB00sQ+cCmwq1feiBhNDc/qzpgqj7NsAr4LnCNpr6RrKAYZT6Lovu6UtK7Me7qkLeWli4GHylmY7wPftH3ftPyKiJh5pvgQaZVtyFSZ1V3TI/mOPnmfB1aV+88A5w9Uu4gYWllsKCLaKYEvItpG9Z8hnlUJfBFRz5DO2FaRwBcRtWWMLyJaZxo+RDojEvgior60+CKiVYb04eQqEvgior4EvohokzzAHBGtpPHRjHwJfBFRT57ji4g2GtXHWbLKWkTU19yaG5MuTKbCv0jaUy5y9vuDVDstvoioreHJjYtsH+pz7jKKdXSXABcAXyj/1pIWX0TUY8Cutg3uCuAuF74HnCLptLo3S+CLiNo0Xm0DFk0sJlZu3YuPTbUw2RnAcx3He8u0WtLVjYhajvE5vkNTrK1deWGyJqTFFxH1VO3mVujqVliYbB9wVsfxmWVaLVXW3Ngg6aCkXR1pn5a0r5yB2SlpVZ9rV0p6qpyJuaFuJSNiODWx2FDFhcnuBf6ynN39Q+Al2/vr1rtKV3cjxeJCd3Wl32r7n/tdJGke8HngEor++COS7rX9eM26RsSwaWZWdzGwWRIUMelLtu+TdC2A7XXAFor1fPYAPwc+OkiBVRYb2ibpHTXuvRzYUy46hKQvU8zMJPBFzBFNPM7Sb2GyMuBN7Bv45OClFQYZ47uufJBwg6RTe5w/plkYSWsnZnze4LUBqhURM8LAmKttQ6Zu4PsC8NvAUmA/cMugFbG93vYy28vmc8Kgt4uIGTCqC4rXepzF9oGJfUlfBP6tR7ZGZ2EiYgiN6CprtVp8XU9Mf4ijZ2AAHgGWSHqnpOOB1RQzMxExR8zZFp+kTcAKiiev9wI3ASskLaXo5T8LfLzMezpwu+1Vto9Iug64H5gHbLC9ezp+RETMgrn8WSrba3ok39En7/MUU84Tx1sopqEjYo4RoCGcuKgir6xFRG0a0TG+BL6IqGcud3UjInpr7JNTMy6BLyJqG8YZ2yoS+CKivrT4IqJVnFndiGij0Yx7CXwRUV8eZ4mI9kngi4hWMTCiC4on8EVELcLp6kZEC42PZpMvq6xFRD0TXd0q2yQknSXp25Iel7Rb0qd65Fkh6aWOBc5uHKTqafFFRG0NdXWPAH9te0e52tqjkrb2WJjsO7Yvb6LABL6IqK+BwFcuE7m/3H9F0hMU6/NM28Jk6epGRE3NLSg+oVzR8V3Awz1Ov1vSY5K+Jel3B6l5WnwRUc/EKmvVLJK0veN4ve31nRkkvQX4KnC97Ze7rt8BnG37VUmrgK8DS2rVmwS+iBjAMYzxHbK9rO99pPkUQe9u21/rPt8ZCG1vkfS/JS2yfehY6wzV1tzYAFwOHLR9Xpl2D3BOmeUU4Ke2l/a49lngFWAMODLZD4+IEdTAGJ8kUSxn8YTtz/XJ83bggG1LWk4xTPdC3TKrtPg2ArcBd00k2P6zjgrdArw0yfUX1Y3KETHEDIw3Mqv7HuAjwA8l7SzT/h74LQDb64ArgU9IOgL8Alht14+6VRYb2lYOOB6ljNQfBt5btwIRMaqa+QKz7Yco1i6aLM9tFA2wRgw6q/tHFM3Pp/ucN/CApEclrZ3sRpLWStouafsbvDZgtSJiRjQ8qztTBp3cWANsmuT8hbb3SXobsFXSk7a39cpYzvCsBzhZC4fvv1RE/DoDYy17ZU3SccCfAvf0y2N7X/n3ILAZWF63vIgYNgaPV9uGzCBd3fcBT9re2+ukpAXl6ydIWgBcCuwaoLyIGDYj2tWdMvBJ2gR8FzhH0l5J15SnVtPVzZV0uqQt5eFi4CFJjwHfB75p+77mqh4Rs2piVrfKNmSqzOqu6ZN+dY+054FV5f4zwPkD1i8ihtkQtuaqyJsbEVFfAl9EtIoNY2OzXYtaEvgior60+CKidRL4IqJdhnPGtooEvoiox+AhfDi5igS+iKhvRF9ZS+CLiHrskV1eMoEvIurL5EZEtI3T4ouIdhnODxBUkcAXEfU09+n5GZfAFxG1GPCIvrKWBcUjoh439yFSSSslPSVpj6Qbepw/QdI95fmH+60DVFUCX0TU5nFX2iYjaR7weeAy4FxgjaRzu7JdA7xo+3eAW4F/HKTeCXwRUV8zLb7lwB7bz9h+HfgycEVXniuAO8v9rwAXl6s81jKUY3yv8OKhf/dXftyVvAiYi+vzztXfBXP3t82F33X2oDd4hRfv/3d/ZVHF7CdK2t5xvL5cYAzgDOC5jnN7gQu6rv9VHttHJL0EvJWa/w5DGfhs/2Z3mqTttpfNRn2m01z9XTB3f9tc/V3HyvbK2a5DXenqRsRs2wec1XF8ZpnWM0+5wuNvAC/ULTCBLyJm2yPAEknvlHQ8xUJm93bluRe4qty/EvgPu/7T00PZ1e1j/dRZRtJc/V0wd3/bXP1ds6Ics7sOuB+YB2ywvVvSzcB22/cCdwD/KmkPcJgiONamAYJmRMRISlc3IlongS8iWmckAt9Ur7OMKknPSvqhpJ1dzziNHEkbJB2UtKsjbaGkrZKeLv+eOpt1rKPP7/q0pH3lv9tOSatms45x7IY+8FV8nWWUXWR76Rx4Lmwj0P1c1w3Ag7aXAA+Wx6NmI0f/LoBby3+3pba3zHCdYkBDH/io9jpLzDLb2yhm2zp1vmZ0J/DBmaxTE/r8rhhxoxD4er3OcsYs1aVpBh6Q9KiktbNdmWmw2Pb+cv8nwOLZrEzDrpP0g7IrPHJd+LYbhcA3l11o+/cpuvGflPTHs12h6VI+bDpXnp36AvDbwFJgP3DLrNYmjtkoBL4qr7OMJNv7yr8Hgc0U3fq55ICk0wDKvwdnuT6NsH3A9piLRWW/yNz7d5vzRiHwVXmdZeRIWiDppIl94FJg1+RXjZzO14yuAr4xi3VpzEQwL32IuffvNucN/Str/V5nmeVqNWExsLn8pNhxwJds3ze7VapP0iZgBbBI0l7gJuCzwP+VdA3wY+DDs1fDevr8rhWSllJ03Z8FPj5b9Yt68spaRLTOKHR1IyIalcAXEa2TwBcRrZPAFxGtk8AXEa2TwBcRrZPAFxGt8/8BNnHxTGa4MTYAAAAASUVORK5CYII=", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAT4AAAD8CAYAAADub8g7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAAsTAAALEwEAmpwYAAAcgklEQVR4nO3df6xc5X3n8ffHPwDhQIF44/CrJGotJIqK01qmUeiuCYEYiw1Jxab2VikkVA5RkILUqktbCSLyT6ouQdslG3IDlqEiDrtJnKDGAVwayUFKCMYyiflVHESEDbFlTPiRX3Dv/e4f59xoMp6593nOnLlz5p7Pyzq6M+c8c84zc+d+/ZznOc/5KiIwM2uTRaOugJnZfHPgM7PWceAzs9Zx4DOz1nHgM7PWceAzs9Zx4DOzkZJ0pqTvSHpC0uOSPlWuP0XSDknPlD9P7vP6K8syz0i6MumYvo7PzEZJ0qnAqRGxW9IJwKPAB4GrgCMR8VlJ1wMnR8T/6HrtKcAuYDUQ5Wv/OCJenu2YbvGZ2UhFxIsRsbt8/BrwJHA6cDlwZ1nsTopg2O39wI6IOFIGux3AurmOuaSGetfuGB0bx7Fs1NUYLY26Agtcy090fsXPeSN+PdC37P0XLouXjkwllX30h79+HPhVx6qJiJjoLifpHcC7gIeBFRHxYrnpp8CKHrs+HXi+4/n+ct2sGhn4jmMZ5+ui+nesIUUT1d9w1iJHvmGK6RFHvpjOKFt/XR+OBwfex+EjUzx8/xlJZZee+uNfRcTq2cpIegvwNeC6iHhVHX+vERGSavsgfKprZhUFUzGdtMxF0lKKoHd3RHy9XH2w7P+b6Qc81OOlB4AzO56fUa6b1UCBT9I6SU9L2ld2PnZvP1bSPeX2h8tmrJktAAFME0nLbFQ07e4AnoyIz3VsuheYGaW9Evhmj5ffD1wi6eRy1PeSct2sKgc+SYuBzwOXAucAGyWd01XsauDliPh94BbgH6sez8yaZzrx3xzeA3wEeK+kPeWyHvgscLGkZ4D3lc+RtFrS7QARcQT4DPBIudxUrpvVIH18a4B9EfFsWZmvUIzCPNFR5nLg0+XjrwK3SlL4GhqzsRcEb+b0VfbbT8RD9B/OO6qzPyJ2AX/V8XwzsDnnmIOc6qaMpvymTERMAq8Ab+21M0mbJO2StOtNfj1AtcxsPgQwRSQtTdOYUd1yaHsC4ESd0rxPysyOMlf/XVMNEvhSRlNmyuyXtAT4HeClAY5pZg0RwNSY9loNcqr7CLBS0jslHQNsoBiF6dQ5KnMF8O/u3zNbOKYTl6ap3OKLiElJ11IMHS8GNkfE45JuAnZFxL0UQ9T/ImkfcIQiOJrZAhAN7b9LMVAfX0RsB7Z3rbuh4/GvgP82yDHmtGhxclEtTi9LxsyJ5P3mzBxZlN4Y17BmpDRh9kjGDIuckwlNJ7ZDck5QUvcJxFRO2bRpYUUdMsoOKALeHM+415zBDTMbN2JqTCeVO/CZWSVBVoO8URz4zKwyt/jMrFWKC5gd+MysRQJ4M8bzBk8OfGZWSSCmxvTOdg58ZlbZdPhU18xaxH18ZtZCYsp9fGbWJsUdmB346pU4DStnGpoWZ/ySli6tf79ZU+aG9L5yDGsqXI6cKWM5ZVOngeVMLZucTC6rRcOZWhapNwat4cLjCPFGZHynG6S5gc/MGm/afXxm1ibF4EY9ZxuSNgOXAYci4txy3T3A2WWRk4CfRcSqHq99DngNmAIm50pjCQ58ZlZZrYMbW4BbgbtmVkTEn//mSNLNFKkr+rkwIg6nHsyBz8wqqXNwIyJ29ks/W6af/DDw3loOhhOKm9kApkJJy4D+FDgYEc/02R7AA5IelbQpZYdu8ZlZJYF4M5JDyHJJuzqeT5QJxlJsBLbOsv2CiDgg6W3ADklPRcTO2XbowGdmlWQObhxOGXToViYp+zPgj/vWI+JA+fOQpG0UOb9nDXyVT3UlnSnpO5KekPS4pE/1KLNW0isd2dFv6LUvMxs/Qdpp7oCnuu8DnoqI/b02Slom6YSZx8AlwN65djpIi28S+OuI2F0e+FFJOyLiia5y342IywY4jpk1VF2DG5K2AmspTon3AzdGxB0UCcq2dpU9Dbg9ItYDK4BtZd6ZJcCXI+K+uY43SJa1F4EXy8evSXoSOB3oDnxmtgBFUNvlLBGxsc/6q3qsewFYXz5+Fjgv93i19PGVw9DvAh7usfndkh4DXgD+JiIe77OPTcAmgOM4Pv3gOZnAcqahLcn4aI5J268yMqeRc/ycKWs509CGNWVtSNPQlDG9LDUjWt40tPTPK95ILpqeEQ6Iqfm7UKMY3GjplDVJbwG+BlwXEa92bd4NnBURr0taD3wDWNlrP+UIzwTAiTplTFOYmLXLuN6IdKBaS1pKEfTujoivd2+PiFcj4vXy8XZgqaTlgxzTzJohENORtjRN5RZfeTX1HcCTEfG5PmXeTnHhYUhaQxFoX6p6TDNrlnFt8Q1yqvse4CPAjyTtKdf9PfC7ABFxG3AF8AlJk8AvgQ2Rk+7ezBqryKvbssAXEQ/B7PekiYhbKSYem9mCI9963szapUgv2dJRXTNrpwi171TXzMzJhsysVYr78bmPz8xaxekl66e0D3RoWdYSp6EBKHUq3JL0ukZG2bzsbQ2YspZjCNPQADSZluVMGZ9BznVamk4vHakZ4eZZcTlLA74jFTQ38JlZo7V6rq6ZtZcTiptZqxS3pfKprpm1jPv4zKxViruz+FTXzFqkmLLmwGdmrTK+Lb7xrLWZNcI0SlrmImmzpEOS9nas+7SkAx1ZGtf3ee06SU9L2ifp+pR6O/CZWSUzo7o1pZfcAqzrsf6WiFhVLtu7N0paDHweuBQ4B9go6Zy5Djb+p7o5swtyZnlkJQZK228szfi4M8rGkoy6ZryvGNLMDeXcizZnNsZUxmyIxPeW8wnkvK+s2RgZv4fkhEcZE2Jm3U19WdZ2lknLcq0B9pXZ1pD0FeBy5sj26BafmVWSmXNjuaRdHcumxMNcK+mH5anwyT22nw483/F8f7luVuPf4jOzkQhgMr3FdzgiVmce4gvAZ8pDfQa4GfhY5j56qiO95HPAa8AUMNn95sqkRP+LIgHwL4CrImL3oMc1s9Eb5qhuRByceSzpS8C/9ih2ADiz4/kZ5bpZ1dXiuzAiDvfZdilFLt2VwPkUUfz8mo5rZqMy5NSRkk6NiBfLpx8C9vYo9giwUtI7KQLeBuC/z7Xv+TjVvRy4q8yu9n1JJ3W9ITMbQ3XeiFTSVmAtRV/gfuBGYK2kVeWhngM+XpY9Dbg9ItZHxKSka4H7gcXA5oh4fK7j1RH4AnhAUgBfjIiJru39Oh9/K/CVnZ2bAI7j+BqqZWbDVleLLyI29lh9R5+yL1B0nc083w4cdanLbOoIfBdExAFJbwN2SHoqInbm7qQMmBMAJ+oU5941a7hW34g0Ig6UPw9J2kZxXU1n4KvU+WhmzRaIyenxvCJuoFpLWibphJnHwCUc3QF5L/CXKvwJ8Ir798wWhrqmrM23QVt8K4BtZV6CJcCXI+I+SdcARMRtFOfe64F9FJezfHTAY5pZE0RLT3XLaSLn9Vh/W8fjAD6ZtWNlTL3JmVq2KCM/wJKcKWOJ+82ZhrY0JzHRsKaspe82R15Snoz3NpkxvS2xXGRMQ1PG9Domc5JkZXwX0mswsFb38ZlZeznwmVmrBGJqTAc3HPjMrLImDlykcOAzs0qirYMbZtZu4cBnZu0y3JsUDJMDn5lV5hafmbVKBExNO/CZWct4VNfMWiXwqe7IKCcD1eKMiy2zyiZmWcuYWpZVNqOuWdPbhpRljawsa0OahJVYB01nTHOcTM+clvO9be492jy4YWYtlPN/WJM48JlZZT7VNbNWKUZ165mrK2kzcBlwKCLOLdf9E/BfgTeAHwMfjYif9Xjtc8yS6bGX8ZxhbGaNEJG2JNgCrOtatwM4NyL+EPgP4O9mef2FEbEqNXevA5+ZVRahpGXu/cRO4EjXugciYrJ8+n2KtBW1cOAzs0qCtKBXBr7lknZ1LJsyD/cx4Nt9q1Jkenw0db/u4zOzyjIGdQ+nnoZ2k/QPwCRwd58i2ZkeK7f4JJ0taU/H8qqk67rKrJX0SkeZG6oez8waJiCmlbRUJekqikGPv4g+eQA6Mz0CM5keZ1W5xRcRTwOrysotpkgZua1H0e9GxGVVj2NmzTXMy1kkrQP+FvgvEfGLPmWWAYsi4rWOTI83zbXvuvr4LgJ+HBE/qWl/ZjYG6hrVlbQV+B5wtqT9kq4GbgVOoDh93SPptrLsaZK2ly9dATwk6THgB8C3IuK+uY5XVx/fBmBrn23vLiv1AvA3EfF4r0Jlp+QmgOM4vqZqDSBnutZQMsINaXpbal0BFg/pf/OcUx8NaXpb4uebM7Us6zszjO/XPKtzrm5EbOyx+o4+ZV+gSFnbN9PjXAZu8Uk6BvgA8P96bN4NnBUR5wH/G/hGv/1ExERErI6I1Ut17KDVMrNhC4ocpClLw9RxqnspsDsiDnZviIhXI+L18vF2YKmk5TUc08waoMYLmOdVHae6G+lzmivp7cDBiAhJaygC7Us1HNPMRm6wEdtRGijwlaMoFwMf71h3DUBE3AZcAXxC0iTwS2BDvyFpMxtDY/rXPFDgi4ifA2/tWndbx+NbKUZmzGyhCd+dxczaqI0tPjNrO7f4zKxtpkddgWoc+Mysmpnr+MaQA5+ZVTau12iMf+DLmc4z4ulHkZNZK+c/0py6ZkxDy6lvDi1K/2vJ+sPKyDanqcRztJzfWcZ3UTnTF3Modb81/W4d+MysdXyqa2Ztk3MPiSZx4DOzakJ5d9ppEAc+M6vOLT4zax0HPjNrHQc+M2uVMb6A2Xl1zawyRdoy536kzZIOSdrbse4USTskPVP+PLnPa68syzwj6cqUejvwmVl1kbjMbQuwrmvd9cCDEbESeLB8/lsknQLcCJxPkVbyxn4BspMDn5lVVleLr0wAfqRr9eXAneXjO4EP9njp+4EdEXEkIl4GdnB0AD2K+/jMZpEzbW88e7sGlN7Ht1zSro7nExExMcdrVkTEi+Xjn1Kkkux2OvB8x/P95bpZOfCZWTXpp7EAhyNideVDFXl7ahtDTjrVne+ORzMbE/X18fVyUNKpAOXPQz3KHADO7Hh+RrluVql9fFuYx45HMxsPmk5bKroXmGksXQl8s0eZ+4FLJJ1cxpZLynWzSgp8893xaGZjoqYWn6StwPeAsyXtl3Q18FngYknPAO8rnyNptaTbASLiCPAZ4JFyualcN6tB+viG1vFoZs2XOmKbIiI29tl0UY+yu4C/6ni+Gdicc7xaBjfq6HiUtAnYBHAcx9dRLTMbthbO3Ki14zEiJiJidUSsXqpjB6iWmc2b4Q5uDM0ggW9oHY9mNh7quoB5vqVezjKvHY9mNgZi6KO6Q5PUxzffHY9mNiYa2JpL4Zkb/Qwhb54y9pl19Jy6ZtwqPCcbWpaMFkDOaVLO55tah6x9NkGkfrh1DcfWs5v55sBnZpU1sf8uhe/OYmat4xafmVU3pi0+Bz4zqyaaOWKbwoHPzKpzi8/M2kSM7+CGA5+ZVefAZ2at0tDpaCkc+MysOg9umFnbuMU3KtMZn3zO9KOcslOJ/+1Np//3qOmMa8tzPoOMb+qwZmtl/bFkvDdNZZRNfXM5v7PU70HmfhvNgc/MWqWh99pL4SlrZlZZHffjk3S2pD0dy6uSrusqs1bSKx1lbhik3m7xmVl1NbT4IuJpYBWApMUUd2nf1qPodyPissGP6MBnZgMYwpS1i4AfR8RPat9zB5/qmlk1qfk2ilbhckm7OpZNffa6AdjaZ9u7JT0m6duS/mCQqrvFZ2aVqFwSHY6I1bPuTzoG+ADwdz027wbOiojXJa0HvgGsTD/8b3OLz8yqqzfL2qXA7og4eNRhIl6NiNfLx9uBpZKWV632nIFP0mZJhyTt7Vj3T5KekvRDSdskndTntc9J+lE5CrOraiXNrJlqzrK2kT6nuZLeLknl4zUUseulqvVOafFtAdZ1rdsBnBsRfwj8B72bpjMujIhVczVzzWwM1dTik7QMuBj4ese6ayRdUz69Atgr6THgn4ENEdUvsZ+zjy8idkp6R9e6Bzqefr+slJm1SY03Io2InwNv7Vp3W8fjW4Fb6zlaPYMbHwPu6bMtgAckBfDFiJjot5NylGcTwHEcn3zwrKCfUXYY049yplQxOaQpTTnT2xYPpws4J3PZ0D6z1LI534OcLHo5U9ZyfmfzrcFVm81AgU/SPwCTwN19ilwQEQckvQ3YIempiNjZq2AZFCcATlx0yph+nGbtMq43Kaj8X7qkq4DLgL/od64dEQfKn4corsReU/V4ZtZA9Y7qzptKgU/SOuBvgQ9ExC/6lFkm6YSZx8AlwN5eZc1sPNU8qjtvUi5n2Qp8Dzhb0n5JV1N0Mp5Acfq6R9JtZdnTJG0vX7oCeKgchfkB8K2IuG8o78LM5l9Q3Ig0ZWmYlFHdjT1W39Gn7AvA+vLxs8B5A9XOzBrLyYbMrJ0c+MysbXIuTWoSBz4zq6ahI7YpHPjMrDL38ZlZ6wzhRqTzopmBLyASp+koZ+rP1FR62ZzsWpNp+w2l370s4z5neenQFqVfupk1bS9Hxm6z+pAypqzpzcm0com/2+L4Od+vnGx3DW5WNbhqs2lm4DOz5mvoxckpHPjMrDoHPjNrE1/AbGatpCbfMmsWDnxmVo2v4zOzNhrXy1mcZc3Mqqsv58asiclU+GdJ+8okZ380SLXd4jOzymoe3LgwIg732XYpRR7dlcD5wBfKn5W4xWdm1QTFxfMpy+AuB+6KwveBkySdWnVn49/iy/lQM2YixGTalf0ASpyRkTMbI+dqfU0vTi+bMXuEjLJZs1JyfmdZM3Pqn21D4gwPADK+M1llM2Ycpc54qmtQIqOPb3nXKexEV/KxuRKTnQ483/F8f7nuxbwaF8Y/8JnZSGRex3d4jtzayYnJ6uBTXTOrJvU0N6GFn5CY7ABwZsfzM8p1laTk3Ngs6ZCkvR3rPi3pQDkCs0fS+j6vXSfp6XIk5vqqlTSzZqoj2VBiYrJ7gb8sR3f/BHglIiqd5kLaqe4WiuRCd3WtvyUi/me/F0laDHweuJjifPwRSfdGxBMV62pmTVNPX+EKYFvZ/7wE+HJE3CfpGoCIuA3YTpHPZx/wC+CjgxwwJdnQTknvqLDvNcC+MukQkr5CMTLjwGe2QNRxOUu/xGRlwJt5HMAnBz9aYZA+vmvLCwk3Szq5x/Z+ozA9SdokaZekXW/y6wGqZWbzIoCpSFsapmrg+wLwe8AqiuHkmwetSERMRMTqiFi9lGMH3Z2ZzYNxTShe6XKWiDg481jSl4B/7VGs1lEYM2ugJt8dehaVWnxdV0x/iKNHYAAeAVZKeqekY4ANFCMzZrZALNgWn6StwFqKK6/3AzcCayWtojjLfw74eFn2NOD2iFgfEZOSrgXuBxYDmyPi8WG8CTMbgYV8W6qI2Nhj9R19yr5AMeQ883w7xTD08GRMacqahrYoY7pW6j6zpqFlTNXKSXKTMw0t4zPISo6UISvhUVZiosTPLOM7k/P9Glbiq/kkQA0cuEjhKWtmVlnWvOsGceAzs2oW8qmumVlvtd1yat458JlZZU0csU3hwGdm1bnFZ2atEh7VNbM2Gs+458BnZtX5chYzax8HPjNrlQCaOalkTs0NfJH2iUZOZq1FGdmq3kguihIzW0XONKXJ0WdO06IGpGTJmZKYM7UrNRvZsLKh5UyFy5q2N3+RSIRPdc2shRo6j3guDfgv3czG0sypbsoyC0lnSvqOpCckPS7pUz3KrJX0SkeCsxsGqbpbfGZWWU2nupPAX0fE7jLb2qOSdvRITPbdiLisjgM68JlZdTUEvjJN5Ivl49ckPUmRn2doicl8qmtmFdWXUHxGmdHxXcDDPTa/W9Jjkr4t6Q8GqblbfGZWzUyWtTTLJe3qeD4REROdBSS9BfgacF1EvNr1+t3AWRHxuqT1wDeAlZXqjQOfmQ0go4/vcESs7rsfaSlF0Ls7Ir7evb0zEEbEdkn/R9LyiDicW2dIy7mxGbgMOBQR55br7gHOLoucBPwsIlb1eO1zwGvAFDA52xs3szFUQx+figtR7wCejIjP9SnzduBgRISkNRTddC9VPWZKi28LcCtw18yKiPjzjgrdDLwyy+svrBqVzazBgvQLwWf3HuAjwI8k7SnX/T3wuwARcRtwBfAJSZPAL4ENEdWjbkqyoZ1lh+NRykj9YeC9VStgZuOqnjswR8RDzJGvKiJupWiA1WLQPr4/pWh+PtNnewAPSArgi92dmZ0kbQI2ARzH8ckfaNY0sAw5Wc6S65AzXWxx+pS1rK9eRua0sZPR+khuLAwpG1rONLSs7/h8TyFr6ZS1jcDWWbZfEBEHJL0N2CHpqYjY2atgGRQnAE7UKeP5aZq1SQA584gbpPJ1fJKWAH8G3NOvTEQcKH8eArYBa6oez8yaJoqbIqQsDTPIBczvA56KiP29NkpaVk4/QdIy4BJg7wDHM7OmqfkC5vkyZ+CTtBX4HnC2pP2Sri43baDrNFfSaZK2l09XAA9Jegz4AfCtiLivvqqb2UjNjOqmLA2TMqq7sc/6q3qsewFYXz5+FjhvwPqZWZM1sDWXwjM3zKw6Bz4za5WIvMt9GsSBz8yqc4vPzFrHgc/M2qWZI7Ypxj/wTWdktsq4kDKm6r9HqzKmizXi66Qh3ae2gRe09hPD+sPO+Qya2qqKvL+pJhn/wGdmozOmU9Yc+MysmoixTS/pwGdm1TX1NHwODnxmVlm4xWdm7dLMGxCkcOAzs2rqu/X8vHPgM7NKguHdAX3YnFDczKqJ+m5EKmmdpKcl7ZN0fY/tx0q6p9z+cL88QKkc+MysspiOpGU2khYDnwcuBc4BNko6p6vY1cDLEfH7wC3APw5Sbwc+M6uunhbfGmBfRDwbEW8AXwEu7ypzOXBn+firwEVllsdKGtnH9xovH/63+OpPulYvBwbLzzvqftje3SGDv6/mWqjvbSG8r7MG3cFrvHz/v8VXlycWP07Sro7nEx1ZF08Hnu/Yth84v+v1vykTEZOSXgHeSsXfQyMDX0T8p+51knZFxOpR1GeYFur7goX73hbq+8oVEetGXYeqfKprZqN2ADiz4/kZ5bqeZcoMj78DvFT1gA58ZjZqjwArJb1T0jEUiczu7SpzL3Bl+fgK4N8jOSv80Rp5qtvHxNxFxtJCfV+wcN/bQn1fI1H22V0L3A8sBjZHxOOSbgJ2RcS9wB3Av0jaBxyhCI6VaYCgaWY2lnyqa2at48BnZq0zFoFvruks40rSc5J+JGlP1zVOY0fSZkmHJO3tWHeKpB2Snil/njzKOlbR5319WtKB8ve2R9L6UdbR8jU+8CVOZxlnF0bEqgVwXdgWoPu6ruuBByNiJfBg+XzcbOHo9wVwS/l7WxUR2+e5Tjagxgc+0qaz2IhFxE6K0bZOndOM7gQ+OJ91qkOf92VjbhwCX6/pLKePqC51C+ABSY9K2jTqygzBioh4sXz8U2DFKCtTs2sl/bA8FR67U/i2G4fAt5BdEBF/RHEa/0lJ/3nUFRqW8mLThXLt1BeA3wNWAS8CN4+0NpZtHAJfynSWsRQRB8qfh4BtFKf1C8lBSacClD8Pjbg+tYiIgxExFUVS2S+x8H5vC944BL6U6SxjR9IySSfMPAYuAfbO/qqx0znN6ErgmyOsS21mgnnpQyy839uC1/gpa/2ms4y4WnVYAWwrbym2BPhyRNw32ipVJ2krsBZYLmk/cCPwWeD/Sroa+Anw4dHVsJo+72utpFUUp+7PAR8fVf2sGk9ZM7PWGYdTXTOzWjnwmVnrOPCZWes48JlZ6zjwmVnrOPCZWes48JlZ6/x/bLJC8XAYli8AAAAASUVORK5CYII=", "text/plain": [ "
" ] @@ -171,7 +171,7 @@ "source": [ "records = []\n", "\n", - "for i in range(1000):\n", + "for i in range(10):\n", " C_t = simulate(C_t)\n", " records.append(C_t)\n", "\n",