{ "cells": [ { "cell_type": "code", "execution_count": 38, "id": "c4d9dbc2", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, \"résolution approchée de -u''=f\")" ] }, "execution_count": 38, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEICAYAAACzliQjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA9qUlEQVR4nO3dd3wUdf7H8dcnjSSUBEhoSSAJvbfQq4oKngqoKEVRUQHFep7K3XlnP/VO/WEFsYAVREVBBJWiIJ2AobfQQw09hdT9/v6Y5W6NCdlAktnsfp6Pxz7I7szsvHcTPjv7nZnPiDEGpZRSvsHP7gBKKaXKjxZ9pZTyIVr0lVLKh2jRV0opH6JFXymlfIgWfaWU8iFa9H2AiDQWkQ0iEneJzxMrIkZEAi5y+REi8tOlZKjInO9dIzfmu1lE5otIcCmt92kR+bQ0nqs8iMjzInJcRI7YncUbadH3ciISBrwH3GSM2VOO6/3DB4Qx5jNjzFXllaEiEpH2wF3AIGNMlt15SkpE9opIrPPnqSJyRwmXjwEeBVoYY+qUfkJ1UVtsyrOJSIAxJg/AGHMG6GtvIu/h+t6WBWPMb8DVZfX8FUAD4IQx5pjdQbyVbul7CecW1hMisgHIEJEAEekqIstF5LSIrBeRvi7z3yEiu0UkTUT2iMgI5+N+IvKkiOwTkWMi8rHz20JR6+znct91GGGJ89/TIpIuIt2c61zqMn93EVkjImec/3Z3mfaLiDwnIsucGX8SkYgiclQXkTkikioip5w/Rxd4rhdFZLVzXbNEpIZz2vlvJKNF5JCIHBaRRwu8pq9E5FMROQvcISL1RGS2iJwUkWQRucdlfn8R+ZuI7HLmXuvcej2vn4jsdOZ8W0TEZdlRIrLVOe1HEWngMq2Zc8jnpIhsF5GbC3svnPPGichi5/rnAxEFphf5d1Ec59b78y73+4pIirvLF/Pc/YD5QD3n38zU0nheVYAxRm9ecAP2AklADBACRAEngGuwPtyvdN6PBCoDZ4GmzmXrAi2dP48CkoF4oAowE/jEOS0WMECAyzr7uWR4Gvi0sHmdj90BLHX+XAM4BdyG9Y1zmPN+Tef0X4BdQBPn6/kFeKmI114TuBEIBaoCXwLfukz/BTgItHK+9q8LyTnNOa01kHr+dTlfUy4wyPk+hgCLgXeAYKCdc/4rnPM/BmwEmgICtHV5TQaYA4QD9Z3L9XdOG+R835s7348ngeXOaZWBA8CdzmkdgOPnf2eFvB8rgNeASkBvIM3l9Rb5d+Hm39lU4HmX+32BFDeW6wmcvsCtZ0meT2+XUCvsDqC3UvpFWgV4lMv9J3AWa5fHfgRudxaR085CGVJgnoXAfS73mzqLXgClW/RvA1YXWPcK4A7nz78AT7pMuw/4wc33oh1wyuX+L7h8YAAtgBzA3yVnM5fp/wY+cHlNS1ymxQD5QFWXx14Epjp/3g4MLCKXOV/cnPdnAOOdP88D7nKZ5gdkYg133AL8WuC53gWeKmQd9YE8oLLLY5+7/F6K/Ltw872dykUU/RL8HZfq8+ntjzcd3vEuB1x+bgAMcX6FPy0ip7G2tuoaYzKwCslY4LCIfC8izZzL1QP2uTzPPqyCX7uUsxZcz/l1Rbncdz16IxPrm8cfiEioiLzrHJI6izW0FC4i/i6zub43+4BAfj/sUXB6vSKm1QNOGmPSisgdg/UNpShFvaYGwOsuv6uTWN8UopzTuhT4XY4ACtvRWQ/rAy+jQL7zivy7KPhEYh1tle68zbvAa1IViBZ97+LaMvUA1hZduMutsjHmJQBjzI/GmCux/rNvwzrCB+AQVmE47/yW49FC1peBNaRynmsRKq59a8H1nF/XwWKWK8yjWN9IuhhjqmENaYBVNM9zHVevj/Xt5fgFph9yue/6Wg4BNUSkahG5DwANS/oCnMuNKfD7CjHGLHdOW1xgWhVjzL2FPM9hoLqIVC6Qz3U9Rf5duDLW0VZVnLcBzocv9Dsvkoj0cvkAKezWy53nUZdOi773+hS4TkSudu5cDHbudIsWkdoicr2zMGQD6VhDFmCNbT/i3BlYBfgX8IUp/IiVJGCoiASKSAJwk8u0VMCBtW+gMHOBJiIyXKydzrdgDbvMuYjXWhU4h7XTuAbwVCHz3CoiLUQkFHgW+MoYk+8y/R/ObwwtscbOvyhsRcaYA8By4EXne9oG6xDLz5yzvA88J9a5ESIibUSkphuvYRLwV+f6EZEwERninDYH6726zfleB4pIJxFpXki+fUAi8IyIBIlIT+A6l1mK/LtwIyNYv/NrRKSGiNQBHnZnIWPMry4fIIXdfnVz/eoSadH3Us7iNBD4G1YBPoC1k9HPeXsUa6v1JNAHa8wc4EPgE6whkj1AFvBAEav5B9ZW7SngGayx4/PrzwReAJY5hxG6Fsh3ArjWmeME8DhwrTHGdevbXROwdrAeB1YCPxQyzydY49FHsHbAPlhg+mKsHakLgVeMMRc6iWwY1r6AQ8A3WGPr853TXsMaq/8Ja2f5B85sF2SM+QZ4GZjuHKLaBAxwTksDrgKGOtd5xDlvpSKebjjQBet3+xTwsct6LvR34Y5PgPVY+3N+oogPR+W5xBi9iIrybiLyC9aOzPcLmRaL9eEWWMS3GaW8im7pK6WUD9Gir5RSPkSHd5RSyofolr5SSvkQj2y4FhERYWJjY+2OoZRSFcbatWuPG2Mii5vPI4t+bGwsiYmJdsdQSqkKQ0QKnuFeKB3eUUopH6JFXymlfIgWfaWU8iFa9JVSyodo0VdKKR+iRV8ppXyIFn2llPIhHnmcvlK2OnMQTu+HtMOQkw45meT6hzA3sB8Bfn5c7VhCQGYqBIdBaA0IrQlV60D1WLuTK1UsLfrKNzny4dgWOPQbHNsK2Wdh4NvWtNn3w65Fv5s9hSgeyrKurjg75FXamO2/f76ojnCPc5nv/wJ+AVCvHdRrDzUbgZ8/SnkCLfrKNxgD4rx64q+vwtIJVqEHCAiBOq3A4QA/P+j9GLsa3clX23P5fmcGaY4gOjeJ4ZOeTcjLN0xY+iprko9QM+Acg5oEM7hZCA1q1fjfuk4kw4FVsGqidT+oCnQZC1f8449ZlCpnWvSV98rLgeQFsH0u7JxvbYmHRUH1OGh9E8R0hegE676fH7n5DuatP8TUZYZ1+/2pUqkSN3Vpxh3dY4mN+N8lZy9rVovkY2lMWbaXd9cdZMKWfLrE+XFnj8Nc2aIO/iO/tb5JHN9hfZNISfzf0M+5U/BON2jSH1oMhNie4B9oy9ujfJNHtlZOSEgw2ntHXbQzKbDkFdj8DWSdhkph0OhyuOxJiGj0h9lPZuQwbfV+PlmxjyNns2hQM5Tbu8UyJCGaqsEXLshnMnP5InE/Hy3fx8HT54gKD+H27g24JaE+YaGFLHsmBX56Enb8CLmZEFIdWgyCno9A9YLXiVfKfSKy1hiTUOx8WvSVV8g6AxnHoWZDSDsCb3Wytqbb3AzxfQvdmt56+CxTl+3l26SDZOc56Nkogjt7xNK3aS38/Uo2/JKX72DB1qN8uGwvq/ecJCTQnxs7RnFH9zga1aryxwVyMmHXQtgyC7Z9D2OXWtlP7YOQcGsnsVIloEVf+YYTu2DF27B+OkR1gDvmWI/nZUPAH68bnu8wLNh6lCnL9rBy90mCA/0Y3D6aO3vE0qR21VKJtOngGaYu38vspEPk5Dvo1TiCUT3i6NMkEr/CPkxyMiDIOXz0xW2w62foeDt0vRfCokslk/J+WvSVd0vdDkv+A5u+Br9AaHUjdL7HKvyFOHMuly8TD/DRir0cOHmOemHB3NYtlmGdYwgPDSqTiMfTs5m2aj+frNzHsbRs4iMqc3v3WG7sGE2VSkXsTjv0Gyx/yxqaEoGWN0CfxyGicZlkVN5Di77ybivehkUvQKdR0O0BqFq70Nl2pabz0fK9fLU2hcycfDrFVufOHnFc1aI2Af7lc25iTp6DeZsO8+Gyvaw/cJqqlQK4uVMMt3eLpX7N0MIXOr0fVk6CtVOg633/O/JHqSJo0VfeJT0VFr8EMV2scfqcTMg9B5Vr/mFWh8OwZGcqU5btZfGOVIL8/bi2bV1G9YijVZS9Y+Xr9p9iyrK9zNt4mHxjuKJZbUb1iKVbw5pIYYdxpqdaw1TB1aydvzvnw2V/s04KU8qFu0VfD9lUni0vG1a8Bb/+n3W0SxXnFn1QqHVzkZGdx8x1KUxZvpfdqRlEVq3EI/2aMLxLfSKr/nF83w4d6lenQ/3qHLmmOZ+u3Mfnq/ezYOtRmtWpyh3dYxnUPorgQJcTuaq4XP3u2FZI/BA2fgmX/R0SRoG//hdWJaNb+spz7VsOsx+wTnZq+ie48plCx7YPnMzk4xV7mb7mAGlZebSJDuPOHrH8qXU9ggI8u71UVm4+s5MO8eGyPWw7kkb10ECGda7Pbd0aUDcs5I8LHN0CP4yHPYuhVgu45j/Wsf7K55Xqlr6I9AdeB/yB940xLxWY3gyYAnQA/m6MecX5eAzwMVAHcACTjTGvl+SFKB+WeQIceXDrTGh0xe8mGWNYufskU5btYcHWo4gI/VvVYVSPWDrUr174UIkHCg705+ZOMQxJiP7v65m0eBfvLtld+Oup3QJGzoJtc+DHv1sfiFr0VQkUu6UvIv7ADuBKIAVYAwwzxmxxmacW0AAYBJxyKfp1gbrGmHUiUhVYCwxyXbYwuqXvwzZ+ZR1v33Ws1a4gP+d3h16e3zKesnwvWw+fLX7LuAJy+5tL7jkICLaO8tn4FQSGQrNr7AuubFVqO3JFpBvwtDHmauf9vwIYY14sZN6ngfTzRb+Q6bOAt4wx8y+0Ti36PijzJHz/KGyeCbG94Pbvftef5siZLD5ZuZdpqw9wMiOHprWrcmePQsbAvUhh+yhGdKnPiC4Nfr+PwhiY+ifYtwxaD4E/vaond/mg0hzeiQIOuNxPAbpcRKBYoD2wqojpo4HRAPXr1y/p06uKbOcCmDUOMo/D5U9Cj0f+W/BLfLSLF6lcKYDbusUyoksDfk0+zpRle5iwYCfv/Lzr90cjiVhDPr++Botftpq93fAe1O9q90tQHsidol/Y/6wS7f0VkSrA18DDxpizhc1jjJkMTAZrS78kz68qsJN74PObIaIJjJgBddtax7UnHfzdce23d4+98HHtXszPT+jTJJI+TSJ/d97BzHUHf3/eQd8noOHlMPNumHIN3L/Gau2glAt3in4KEONyPxo45O4KRCQQq+B/ZoyZWbJ4ymudbz1QIw6GfwGxvTieLUxbuPO/Z7DGRVTmmetbXvgMVh/TMLIKzw5sxaNXNeXLxANMXb6X+z5bR1R4CLd1a8DQTm0JH/Or1c/nfMHPz9VOnuq/3BnTD8DakXsFcBBrR+5wY8zmQuZ9GpcxfbG+f38EnDTGPOxuKB3T93LJC+CbsTD4XWh0BZsPnWHKsr3MXn+InDwHvZtEcmf32KJ71aj/yncYFm49ypRle1mx+8QfewmlrIWv7nQO95R4VFZVIKV6Rq6IXANMwDpk80NjzAsiMhbAGDNJROoAiUA1rEMz04EWQBvgV2Cj83GAvxlj5l5ofVr0vZQxsPQ1WPgcJrIZSzu8ypsb/Ap0pYylUa3SaXzmawrrGnp/y2y6rH4QOZMC/V+ETnfrBVy8lLZhUJ4lN8s60WrjDJJrXc3oM3ey+4yj+P7zqsQKXh+gZQ0HE0MnU//4Emg3Aq6dAAFl02RO2UfbMCiPcnTlNGpvnMEExy1M2H89XeLCefy6OK5sUbvEvevVhdWoHMS4yxoxunc8P2w6wpRle+izfzSPVarFfUmfcSKsJTUvG2d3TGUT3dJXZSr58AmemZvMrztT6RCwl4Zte3FHj1ha1tPjyMvT+gOnmbJsD8c3LWBFfjP6Nq3DU9c2p35EIRd4URWSbukr26VvXUiVGWPI5K88emVvhne5kppVPKPxma9pGxPOhKHtOXa2OZ+u2s/cZYnkvXMX54a/TUijXnbHU+XIs7tRqQrLsXEmlb64hbP5lXhqSDceuKKxFnwPUKtaMH++sgkvXt8Mk5dDwGeDMVvn2B1LlSMt+qr0rZyEfD2K3xzxrLn8c9q0aGF3IlVApw4dWNT9Ezbn18fMGAlJn9sdSZUTLfqqdG34En54gp/yE/i8yesM79PG7kSqCHddlcCkBv/Hivzm8O29sEnPnfQFOqavStXR6Kv42G8UP1W5lm+GdPL6/jgVmZ+f8NKwbgx6/Z+MyvuC6+r1prrdoVSZ0y19demMgWWvk5t+gvtnbObD3Kt5+7bO2jqhAggPDeKN27rwfNYQHv52N47sDFg71fqdKq+kRV9dGocD5jwM8//JohlvsmbvKV66sbXVAkBVCG2iw3nq+hYs3pHK4i8mwHcPwfx/aOH3Uroppi6ewwGz74ekz9jVbAxjkhK4rWsDBraLsjuZKqHhneuzdu8pRiUZlra6lajlb1pF/6rntW2Dl9Giry6OMTD3UUj6jFOd/szA1V1oG1OVJ69tbncydRFEhBcGt2bzobNcm3w9v7YNosqKt6yJWvi9ig7vqItz7hTsWkRut4cYtqMvgf5+vDOiA5UCvPMqVr4gJMifibd2INcBtx68kfxOo2H9dEg7Ync0VYq06KuSMwZCa2Du+YUnTg9m+7F0JgxtT1S4d1yj1pfFR1bh3ze1ISnlDM/mjoR7l0G1unbHUqVIi74qmaUT4Nv7ID+PaRvTmPnbIR68vDF9mkTanUyVkmta1+WunnF8tHI/s3c7rA/5Bc/A2o/sjqZKgRZ95b61U2HBU5CXxcaDZ3l69mZ6NY7gwSsa251MlbLxA5qR0KA647/ewM7Dp+DIBuuons3f2B1NXSIt+so923+AOY9Ao36c7v8W905LIqJKEK8Pba+tkb1QoL8fbw3vQGiQP/dO30jGoCkQ0wW+vse68pmqsLToq+KlJMKXd0CdNjhumsqfv97C0bNZvD2iAzUq68U4vFWdsGDeGNqe3anpjP9uF2b4dIhsBjNuh8Mb7I6nLpIWfVW8nHSIaAQjvmTiiqMs2naMJ//Ugvb19aR9b9e9UQSPXtWU79Yf4uPfzsCIGVA5Ek7utjuaukh6nL4qWn4u+AdCfF8YvYRlu0/y6k/bua5tPUZ2a2B3OlVO7u3TkHX7TvH891toHd2NDuNWQYCzTbYxegx/BePWlr6I9BeR7SKSLCLjC5neTERWiEi2iPylwLQPReSYiGwqrdCqHORkwof9YfV7ABxJy+HBab8RH1mFl25orY3UfIifn/Daze2oXS2Y+z9bx8ls5+9+09fw2RDIy7E3oCqRYou+iPgDbwMDgBbAMBEp2CD9JPAg8EohTzEV6H9pMVW5MgZm3QcH10K1KHLzHdz/+TrO5eYz6dYOVNZGaj4nLDSQSbd25HhGDg9N/418h7GKffJ8+O5B7dNTgbizpd8ZSDbG7DbG5ADTgYGuMxhjjhlj1gC5BRc2xizB+lBQFcXil61D8/o9Dc2u4aV520jcd4qXbmxDo1raSM1XtYoK45nrW/LrzuO8sXAntBsGff8G66fBstftjqfc5E7RjwIOuNxPcT5WqkRktIgkikhiampqaT+9ctemmfDLi9B2GPR4iLkbD/PB0j3c0T2W69vWszudstnQTjHc2CGaNxbt5Jftx6DP49BiECx4Gnb8aHc85QZ3in5hg7el/l3OGDPZGJNgjEmIjNSzO22Tfgzqd4frXmfX8Qwe+3I97euH87drtJGashqzPT+oFU1rV+XhL5JIOX0OBr0DdVpbw4HK47lT9FOAGJf70cChsomjbNd1LNwxh0yHP/d+upZKgf68PbwDQQF6dK+yWI3ZOpKfbxj32Tqy/YJh1I9w2d/sjqbc4M7/5DVAYxGJE5EgYCgwu2xjqXKVlwOf3wI7fgLAiB9//2YTO4+l8/rQdtTTRmqqgLiIyvxnSBvWp5zh+TlbISjUmnBwLXwzFvLz7A2oilRs0TfG5AH3Az8CW4EZxpjNIjJWRMYCiEgdEUkB/gw8KSIpIlLNOW0asAJo6nz8rrJ6Meoi/fQk7PgBctIA+GzVfr757SAPX9GEXo11qE0Vrn+rutzTK45PVu5jVtJB68FjW60duwuftjWbKppbx94ZY+YCcws8Nsnl5yNYwz6FLTvsUgKqMrZhBqx+F7qOg1Y3siHlNM9+t4U+TSJ54PJGdqdTHu7x/s1IOnCa8V9vpHndajRpfyscXAfL37R69TS/zu6IqgAdqPVlRzbB7AehQQ+48hlOZeRw76friKxaiQm3tMNPG6mpYpxvzFa5UgBjP11LenYe9H8R6nWwWnCf2GV3RFWAFn1ftukrCAmHm6bgkAAemZHEsTSrkVp1baSm3FS7WjBvDmvP3uMZPPH1Box/ENz8MfgFwMqJdsdTBWjR92VXPAWjf4GqtXn752R+2Z7KP69tQbuYcLuTqQqmW8OaPHZ1M77fcJipy/dCeAzc9RMMeNnuaKoALfq+aMMMSN1hNcqqWoelO4/z2oIdDGxXj1u7aiM1dXHG9omnX/PavPD9VtbuOwURjcHP37rG7vZ5dsdTTlr0fc3BtdZY62JrC+zwmXM8OP03GkVW4UVtpKYugYjw6s1tqRcewv2fr+NEerY1Yf4/rR78RzbaG1ABWvR9S9YZ+GoUVK0D1/yHnDyHdXJNbj4Tb+1IaJA2UlOXJiwkkHdGdOBERg4PTU+yGrNd/S9r39FXo6zurcpWWvR9hTHw3cNw+gDc+AGE1uDFeVtZt/80L9/Uhka1qtidUHmJVlFhPDewJUuTj/P6gh1QOQJumAzHd8KPf7U7ns/Tou8rNs+0bpf/Hep34bv1h5iybC939ojl2jbaSE2Vrls61WdIx2jeWJTMz9uOWRfi6fkwrJ0KW+fYnM63adH3FU3/BAP+DT0eIflYOuO/3kCH+uH8dYA2UlNl47lBrWhetxoPf5HEgZOZcNnfoddfILaH3dF8mhZ9b5eXDdlpEBgMXcaQkev4XyO1EdpITZWd4EB/Jo7ogMNhGPf5OrKNH1zxDwipbvV70v48ttD/8d7u5xdgUk/IOosxhr99s5Hk1HTeGNqeumHaSE2VrdiIyrxyc1s2pJzh2e+2WA9mp8MHV8KSf9sbzkdp0fdme5fBsjcg/jIIrsanK/cxK+kQf+7XhJ6NI+xOp3zE1S3rMKZPvLORXwpUqgK1msOSVyAl0e54PkeLvrfKOmu1uK0RB1c9T9KB0zw7ZwuXNY1k3GXaSE2Vr8euakqXuBr8deZGth9Js87UrVoXvhmjh3GWMy363uqH8XA2BQZP5lReEOM+W0etqsH8nzZSUzYI8PfjzeHtqRocyL2friWNUBj0NpxIti61qMqNFn1vlHsOTu2FXo/iiErg4S+SSE3LZuKtHQgP1UZqyh61qgbz1rD27DuZaTVmi+sDXcbCniW6tV+O9BRMbxQYArd/B8bBm4uSWbwjlecHtaJNdLjdyZSP6xJfk8evbsqL87bx4bK93NXvaUCso8tUudAtfW9iDCydYDW48vNnya7TTFi4g8HtoxjRpb7d6ZQCYHTveK5qUZsX524l8eA5q+Bnp8Fvn9odzSdo0fcmm76GBU/BppkcOn2Oh6b/RuNaVXhhcCttpKY8hojwnyFtiaoewrjP13E8PRsSP4RZ4/Rs3XKgRd9bpKfC3McgqiM5He/hvs/WkZtvtJGa8khhIYFMHNGR05m5PDjtN/I73wt1WsP3f4Zzp+yO59XcKvoi0l9EtotIsoiML2R6MxFZISLZIvKXkiyrSsm8xyAnHQa+zQvztpN04DT/vqkNDSO1kZryTC3qVeO5Qa1YvusEry3aDde/BRnH4acn7Y7m1Yot+iLiD7wNDABaAMNEpEWB2U4CDwKvXMSy6lJtnQObv4E+jzP7cBgfrdjHXT3juKZ1XbuTKXVBNyfEcEtCDG//vIuFZ+pCjwetsf1di+yO5rXc2dLvDCQbY3YbY3KA6cBA1xmMMceMMWuA3JIuq0pB/W7Q61F2NrqL8V9vIKFBdcYPaGZ3KqXc8szAlrSsV41HvkjiQJsHof1tUD3O7lhey52iHwUccLmf4nzMHZeyrHKHwwGVa5LR82/cO30joUH+vDW8A4H+urtGVQxWY7aOANz7xWayrnndOpNclQl3KkNhh30YN5/f7WVFZLSIJIpIYmpqqptP7+N2zof3r8CcOcj4mRvZ7WykVidMj3lWFUv9mqG8dnM7Nh08yzPfbbEOTJg2DPavsjua13Gn6KcAMS73o4FDbj6/28saYyYbYxKMMQmRkZFuPr0Py8m0jnTIyeDTjRl8t/4Qj17VlO6NtJGaqpj6tajNvX0bMm31fmZtPmldU3f2/ZCbZXc0r+JO0V8DNBaROBEJAoYCs918/ktZVl3I4pfh9H52dH6WZ+clc0WzWtzbp6HdqZS6JI9e2YRu8TV5Ys5uDvR4AY7vgGUT7I7lVYot+saYPOB+4EdgKzDDGLNZRMaKyFgAEakjIinAn4EnRSRFRKoVtWxZvRifcXQzrHiL7NbDuWNhILWrBfPazdpITVV8Af5+vDGsPdWCAxm5JIzc5oPh19fgxC67o3kNMcbd4fnyk5CQYBITtc92kb6+B7NrIfeFv8vCfXl8fW93WkeH2Z1KqVKzes9Jhr23khsbB/DykbuQRv1gyBS7Y3k0EVlrjEkobj49xKMiuv5NPmv8OvN25/D09S214Cuv0zmuBuP7N2PG9lzmtnwVrn3N7kheQ4t+RZJ5EnIy+GX3Wf6x2o8bOkQxrHNM8cspVQHd3SuO/i3r8ODKKqw5aqxr6moL5kumRb8imfc4ue/04i/T19C0dlVeGNRaG6kpryUi/HtIG2Kqh/DQp6vInXw5zP+n3bEqPC36FcWun2Hjl3yV3YVsRwATb+1ISJC/3amUKlPVggOZeGtHTmbDgvRYzJr34eBau2NVaFr0K4K8HJj3OCeConj61FX8Z0gb4iIq251KqXLRvG41XhjUmsdOXEd6YE2Y8wg48u2OVWFp0a8IVr8Lx3fwWPowRvZqSv9W2khN+ZYbO0ZzXedmjM8YDofXw1o9kudiadH3dMaQvv0XFpv2pNW/gsf7ayM15Zueuq4F++pcySpakrXuC+tKcarE9OoaHi49J5/rT9xPfkAaM7SRmvJhwYH+TLw1gZFvPEzVrEhm5DkIDtT9WiWlFcSDmWNb+dcXi9h7IpOXhvegdjVtpKZ8W0yNUP4xtA8bDmfy4jer4fR+uyNVOLql76kc+Zz85A7uPJNG9FWz6Nawpt2JlPIIlzerzbi+8Vy7fAipR2oT+cAC0EOX3aZb+h5q3/x3qJm2jUW172Bsn0Z2x1HKo/z5qmYsq3EDkScTSVn6md1xKhQt+h7o7ImjhK94iXV+rRh6x8PaSE2pAvz9hEGjxrNN4gha9E8y0s7YHanC0KLvgXZ/9SSVTSaVB75KWOUgu+Mo5ZEiqoWSd9XL1DIn2DHzObvjVBha9D3Muew89h45zuJq19G0bVe74yjl0Vp1u5pfgy8jfV8SuXl6wpY7tOh7mC/XpfDwubupOli7CirlDsd1b3Jb5sN8t+Gw3VEqBC36HiTvQCI/L15Ih/rhdIrTo3WUckfvFtE0rV2Nb35eiTm2ze44Hk8P2fQUjnwyvhrH05kn2XbtL9o9Uyk3iQhjejcgYdYozs6IIWzcIj2E8wJ0S99DmKTPCTuzjY9CRnJli3p2x1GqQrmuXQyfBd1E2PF1sGWW3XE8mhZ9T5CdTs78Z1nnaESTy0fqIZpKlVCgvx+1et/FVkcM2T/8A/Ky7Y7ksdwq+iLSX0S2i0iyiIwvZLqIyBvO6RtEpIPLtIdEZJOIbBaRh0sxu/dY/gaVzh3jrcBRDO4YbXcapSqkW7rEMcFvJJXS9sPq9+yO47GKLfoi4g+8DQwAWgDDRKRFgdkGAI2dt9HAROeyrYB7gM5AW+BaEWlcaum9xJFMmJHXh069+lMpQBtIKXUxqlQKoHG3gSzMb8/p1BS743gsd7b0OwPJxpjdxpgcYDowsMA8A4GPjWUlEC4idYHmwEpjTKYxJg9YDAwuxfxe4V9nB/Cs/ziGd6lvdxSlKrTbu8cyzjzGvx0j7I7isdwp+lHAAZf7Kc7H3JlnE9BbRGqKSChwDaBX8j7v2DZS18xkzoaDDO9Sn7CQQLsTKVWhRVatxA0d6/PV2hRObl8KJ/fYHcnjuFP0C9urWPDqBYXOY4zZCrwMzAd+ANYDeYWuRGS0iCSKSGJqaqobsbzAgqeo8sMDhPudY1SPOLvTKOUV7ukVT3B+OlW+uAkWPmt3HI/jTtFP4fdb59HAIXfnMcZ8YIzpYIzpDZwEdha2EmPMZGNMgjEmITIy0t38Fde+FbDjBybmXssV7ZpQJ0x75StVGuIiKtOzVTxTHQNg80w4lGR3JI/iTtFfAzQWkTgRCQKGArMLzDMbGOk8iqcrcMYYcxhARGo5/60P3ABMK7X0FZUxsOAp0gMjmJxzNaN7x9udSCmvMqZ3Q97MuoasgDDd2i+g2KLv3AF7P/AjsBWYYYzZLCJjRWSsc7a5wG4gGXgPuM/lKb4WkS3Ad8A4Y8yp0nwBFdL2eXBgFRPybqBn8xga165qdyKlvErbmHBaxkfznhkIuxbCniV2R/IYbrVhMMbMxSrsro9NcvnZAOOKWLbXpQT0TobDNToz9VBPpvdpaHcYpbzSmD4NGTvlckZGLifsbMERad+lvXdskNd4AENmhdC2QTAJsTXsjqOUV+rbJJK4OjW52fEKP7TpW+jRJr5I2zCUp9wsSPyQeUn7SDl1jjE6lq9UmRERxvSJZ/uxTH7edsQaVnVoz30t+uUp8QOY8wjLfp5Lw8jK9Gte2+5ESnm1a9vUo15YMCt++gqmDYX10+2OZDst+uUlOw2WvMKpOj2YfjyWMb0bamM1pcpYoL8fd/WK573DcWREtIFfXrS+cfswLfrlZdUkOHeS/3MMpVbVSgxsr+2TlSoPQzvFEBYSxOTAkXDmACR+aHckW2nRLw/nTsPyNzlbvx8f76/JqJ5x2lhNqXJSuVIAI7s14I299ciM6gFLX4OcDLtj2UaLfnnISIWajXiHm6laKUAbqylVzm7vHkugvx8fBw+HwFA4tdfuSLbRol8eIhqzf/B3TN5ZheFd61MtWBurKVWeIqpUYkjHaF7bVpNjdy6H2i3tjmQbLfplbecCyDjO+8v24O8n2lhNKZvc0yuePIeDqStSIPccHFxrdyRbaNEvSxnHYcZIsr8fz4zEAwxuH0XtatpYTSk7xEZUZkCrunyych+53z4An94IWWftjlXutOiXpWWvQ945plUaQlauQxurKWWz0b3jScvKY07IQDh3Cla/a3ekcqdFv6ykHYXV75HX8kYmrBeubFGbRrW0sZpSdmobE063+Jq8vCEUR5MBsPxNyDpjd6xypUW/rCybAPk5zKp2K6czcxnbR7fylfIEY/s25MjZLBbVucsq+Csn2h2pXGnRLwvGQOYJHG2H8do6BwkNqtOxgTZWU8oT9G4cQbM6VXk5KQjT7DpISbT+z/oILfplQQRumMx3DZ7g4OlzjNH2yUp5DBFhbJ+G7DyWzi8tn4MRX1r/Z32EFv3SlnEcju/EGMOkJftoVKsKVzSrZXcqpZSLP7WpS1R4CBOXHbEKfnqqz4zta9Evbcteh3e6sXLDNrYePsvo3vHaWE0pDxPo78fdveJYvfck67fthNfb+MzYvhb90pRxAtZ8AC0H8ebqM9SuVomB7bSxmlKe6JZOMYSHBvL26jMQ39cq+tlpdscqc1r0S9PKdyA3g+1NRrN81wlG9dDGakp5qtCgAEZ2bcD8rUc50GocZJ22Ntq8nBb90nLuNKyeDM2v542NAdpYTakKYGT3WIL8/Xh7RzVoeAWseAtyMu2OVabcKvoi0l9EtotIsoiML2S6iMgbzukbRKSDy7RHRGSziGwSkWki4p19CFISIT+XQ23vZ97Gw4zo2oCq2lhNKY8WUaUSNyfEMHPdQU4lPGQdiLFnid2xylSxRV9E/IG3gQFAC2CYiLQoMNsAoLHzNhqY6Fw2CngQSDDGtAL8gaGllt6TNO4Hj27jnW2hBPj5cWePWLsTKaXccHevOPIcDt7dWxseSoKm/e2OVKbc2dLvDCQbY3YbY3KA6cDAAvMMBD42lpVAuIjUdU4LAEJEJAAIBQ6VUnbPcfYwGMPx/BC+TEzRxmpKVSANalZmQOu6fLZyH2khUdaDXnyRFXeKfhRwwOV+ivOxYucxxhwEXgH2A4eBM8aYnwpbiYiMFpFEEUlMTU11N7/9cjLh3V4w/x98vHwvOfkORmvLBaUqlLG9G5KWnce01fvh5xdhYg/Iz7U7Vplwp+gXdpB5wXOWC51HRKpjfQuIA+oBlUXk1sJWYoyZbIxJMMYkREZGuhHLQ6z7CDJSORd/NR+t2MeVzWvTMLKK3amUUiXQOjqMHo1q8sHSPeTWaQun9sDGL+2OVSbcKfopQIzL/Wj+OERT1Dz9gD3GmFRjTC4wE+h+8XE9TF42LHsDGvRg2pFozpzL1ZYLSlVQY3o35OjZbL5JbwV1WsOvr4Ij3+5Ypc6dor8GaCwicSIShLUjdnaBeWYDI51H8XTFGsY5jDWs01VEQkVEgCuAraWY314bZkDaIfJ6PMIHS/fQKbY6HRtUtzuVUuoi9GocQfO61Zj86x4cvf4CJ5Jhy7d2xyp1xRZ9Y0wecD/wI1bBnmGM2SwiY0VkrHO2ucBuIBl4D7jPuewq4CtgHbDRub7Jpf0ibJP0GdRpzZz05hw8fY6xupWvVIVlNWKLJ/lYOotMZ6jZCFZ530VWxHhgS9GEhASTmJhod4zi5WRizh5kwKeHyXcYfny4t/bZUaoCy8t30Oc/v1AvPJgvr6sEYTFQpWI0TBSRtcaYhOLm0zNyL4Yx1lhfUCiLT4Sx7UiaNlZTygsE+PtxT6841uw9xdq8uApT8EtCi/7F2Lcc3kqAY9t4d/Fu6lQLZmC7gkexKqUqopudjdgmLd4Nx5Phg6vhUJLdsUqNFv2LsWwCZJ1lQ0Y4K3af4K6ecQQF6FuplDcIDQpgZLdY5m85yu5zoXBsi9Uy3UtopSqpo5th50/QZSyTlh+kanAAQzvHFL+cUqrCuL1bA4ID/Zi0KhUS7rSO4jm5x+5YpUKLfkktex0CK7O/4TDmbTrCrdpYTSmvU9PZiO2b3w6S2nIU+AVYHTi9gBb9kjh9ADZ+BR3v4N01pwjUxmpKea27e8aT7zC8v/4ctLkFfvvUuqxiBadFvySq1YObPuREm3v4cm0KN3aMolZVbaymlDeqXzOUa1rX5fOV+0lPuA+ufBaCKtsd65Jp0S8JP39oOYgpG3PIzXdwdy9trKaUNxvbx2rE9mlyJegyBoJC7Y50ybTou2vlRFj8bzKycvl4xV6uaqGN1ZTydq2iwujZKIIPl+4hOzcP1k612q9UYFr03ZGbBUtegZREpiemcDYrT1suKOUjxvSJ51haNrOSDlsFf+GzkJ9nd6yLpkXfHRtnQOZx8rrcxwe/7qZzXA3a19fGakr5gp6NImhZrxrvLtmFo+s4OHMAthbsOVlxaNEvjjGw4m2o3ZrZZxpy6EwWY/UiKUr5DBFhTJ+G7ErNYEF+e6gRbx2+6YF9y9yhRb84uxZC6jZMt/t4d8kemtSuQt8m3tePQylVtGta1SG6egjv/roXut4HB9fCgdV2x7ooWvSLUzkS2g5jcVBvth9NY0zvhtpYTSkfYzVii2ftvlOsqz4A4vuCqZgXWNGiX5y6bWHwJCYtPUDdsGCua1vP7kRKKRsMSYimemgg7yw/AiNnQYOKeRFALfoXsuFLOLmbpAOnWbn7pDZWU8qHhQYFcHv3WBZsPcrOo2mQeRJ2/Wx3rBLTClaU9GMw6z5Y/ibvLt5FteAAhnaub3cqpZSNRnaLJTjQj8lLdsNPT8IXt0HWGbtjlYgW/aKseR/yczjQ9E5+2HyE27o1oEqlALtTKaVsVKNyELckxPBt0kGOt7gDctJg3Sd2xyoRLfqFyT1nFf0mA5i4SQj09+P27rF2p1JKeYC7e8XjMPBeclVo0BNWTapQJ2u5VfRFpL+IbBeRZBEZX8h0EZE3nNM3iEgH5+NNRSTJ5XZWRB4u5ddQ+jbMgMwTnGp3D1+tTeHGDtHaWE0pBUBMjVD+1Loun63aT0bH0dbJWtu+szuW24ot+iLiD7wNDABaAMNEpEWB2QYAjZ230cBEAGPMdmNMO2NMO6AjkAl8U2rpy0r6UYjuxPv765Gb7+CeXnF2J1JKeZDRveNJz87jkxPNoUZDOLjO7khuc2dLvzOQbIzZbYzJAaYDAwvMMxD42FhWAuEiUrfAPFcAu4wx+y45dVnr8zjpI77nk5X76d+yDvHaWE0p5aJVVBi9GkfwwfL9ZN/9M1z1nN2R3OZO0Y8CDrjcT3E+VtJ5hgLTilqJiIwWkUQRSUxNtfFCBSd2ATA98SBns/IY3VtbLiil/mhM74akpmXz7Wbn0TvnTtkbyE3uFP3CTj8t2HTigvOISBBwPfBlUSsxxkw2xiQYYxIiIyPdiFUGTuyCNzuSt+oDPli6hy7aWE0pVYQejWrSKqoa7y7ZjWPdp/BqMzh72O5YxXKn6KcArlf+jgYOlXCeAcA6Y8zRiwlZbta8D37+/JjXnsNnsrR9slKqSCLCmN4N2Z2awdLcppCXDYkf2h2rWO4U/TVAYxGJc26xDwUK9hWdDYx0HsXTFThjjHH9yBvGBYZ2PEJ2Ovz2KabFQF5fnUbT2lXp29SmbxxKqQphQKs6xNQI4f/W5mCaXAVrp1jF34MVW/SNMXnA/cCPwFZghjFms4iMFZGxztnmAruBZOA94L7zy4tIKHAlMLOUs5euDdMh+yzr6tzCjqPpjOkTj4g2VlNKFe18I7bf9p9mR4PhkJEKm7+1O9YFifHAntAJCQkmMTGxfFf6/pWQn8PN5kVSTp1j8eOXEeiv564ppS7sXE4+PV5eRIfoaryfPg5Ca8BdP5V7DhFZa4xJKG4+rWrn3fYNm7u/yuq9pxjVM04LvlLKLSFB/tzeLZYF24+zv8+rMOQjuyNdkFY2sK6AU6kKbyYJYSGBDNPGakqpEhjZrQEhgf68sT0cqhU8RcmzaNE/tRcm9uDgxiX8uOUIt3VtQGVtrKaUKoHqlYO4pVMMs5IOkpqcCJ8Mtjr1eiAt+mveh9RtfLw5RxurKaUu2l0943AY+Gr9cdi1CNZOtTtSoXy76OdkwrpPyGr8J6ZszOWmjtFEVq1kdyqlVAUUUyOUa9vU5a31kBt3uXXMfn6u3bH+wLeL/sYZkHWamQHXkOtwMLqXtlxQSl280b3jycjJ56cq10PaYdha8JQm+/lu0TcGVk0mv1ZrXtxSnQGt6hAbUdnuVEqpCqxlPasR29NbonCEx8KqyXZH+gPfLvq9/sz8OveQlpXPmN7ackEpdenG9mlIakYuqxqMhubXgcNhd6Tf8d2i7+dHTvMbeHpbNF3ja9A2JtzuREopL9C9YU1aR4Xx910tye86Dvw8q8x6VpryknEcfn2NH9Zs4chZbaymlCo9IsKYPvHsPp7Bwo37IGkaZJ21O9Z/+WbR/+0TWPgMs5Ym0axOVfo00cZqSqnS079lHerXCGX+z4vg27Gw4Qu7I/2X7xV9hwMSp3AqsjMLT1TXxmpKqVJnNWKL48sjtcio0QoSp1j7ET2A7xX9XYvg9D6m5lxBVHgI17apZ3cipZQXGpIQQ83KlZghV8GxzXBgld2RAF8s+okfkBsSwTtHm3OXNlZTSpWR4EB/bu8ey78PtiI/qKrHXGDFtyqeIx+Mgx8r9Sc0JIRbOsUUv4xSSl2k27o2gMDKLK/cz7ocqyPf7kg+VvT9/Enu9wEPHB3AyG7aWE0pVbbON2Ibc3QQh276Dvz87Y7kQ0U/PxdOH+D9X3cT5O+vjdWUUuXi7l5xZBPEh8v2QtYZ20/W8p2iv20O5vU2JP/2C0MSoomooo3VlFJlL7p6KNe1qcum1YswrzaH3YtszeM7RT/xQ84G1iYpP467e2pjNaVU+RnduyHrcqI5R5B1+KaN3Cr6ItJfRLaLSLKIjC9kuojIG87pG0Skg8u0cBH5SkS2ichWEelWmi/ALcd3wp4lTM3uw9WtorSxmlKqXLWoV42uTeoxI68PZvtcOHPQtizFFn0R8QfeBgYALYBhItKiwGwDgMbO22hgosu014EfjDHNgLbA1lLIXTKJH5IvAXyS3YcxfXQrXylV/sb2ief9rL7WSVrrPrYthztb+p2BZGPMbmNMDjAdGFhgnoHAx8ayEggXkboiUg3oDXwAYIzJMcacLr34bsjPw2z8kkXShSYN42kTHV6uq1dKKYBu8TWpEdWYVf4dMOs+su3wTXeKfhRwwOV+ivMxd+aJB1KBKSLym4i8LyKFjq2IyGgRSRSRxNTUVLdfQLH8A/iu25c8e+4mxmhjNaWUTUSEMb0b8kzmTSztMsm2wzfdKfqFNaYp2ESiqHkCgA7ARGNMeyAD+MM+AQBjzGRjTIIxJiEysvQaoDkchjdWn6VKncb0bhxRas+rlFIl1b9VHTJrNOeVpACMTb143Cn6KYDrqavRwCE350kBUowx55tOfIX1IVA+Tuzi9KT+SOo2xmpjNaWUzfz9hHt6xXPm4DaOf3w7nC1YSsueO0V/DdBYROJEJAgYChS88ONsYKTzKJ6uwBljzGFjzBHggIg0dc53BbCltMIXa93HhB1bTUi1mlzTum65rVYppYpyU8dowkKCiNwzC5I+L/f1F1v0jTF5wP3Aj1hH3swwxmwWkbEiMtY521xgN5AMvAfc5/IUDwCficgGoB3wr9KLfwH5ueSu+5SF+e0Z3LujNlZTSnmE4EB/+vXoxvL8FuSs+ajcz9B1q/mMMWYuVmF3fWySy88GGFfEsklAwsVHvEg7fiDw3HG+C7ibl7WxmlLKg9zWrQHPL+5H97Q3YO8SiO9bbuv22s3fjBUfctjUIL7bQEKDtLGaUspzhIcGUb3jDZwyVchcWb4tl7226C/Ibc0kx2BGdtfDNJVSnueOPs2Ymj+AxLQa5bperyz6R89m8dj+bjg6jqKmNlZTSnmgqPAQ9re+n7EH+3M6M6fc1ut9Rd+Rz6pZEwl0ZHJPL225oJTyXKN7x5OVk8vPP3xdbtfQ9bqin7FtPtfvepo/x+6jfs1Qu+MopVSRmtetxhNRGxm8YSzZe5aXyzq9rugf/fldTpiqdOl/m91RlFKqWO2uvJV0E0zKwnfLZX1eVfSzTx+mfupiVla9ilYNSq+Vg1JKlZXOTWNYGtKXqIPzyM88Xebr86qiv+2HyQSQT60+99gdRSml3CIiVO12F8HksHV+2V9gxWuKvsNhOJW8ik0BLUhI6Gp3HKWUclvXnv1I9ovlyJZlZd6IzWvOWsrMzeeH5i9xWWwIrbSxmlKqAvH392PjldNZdzSfnnkOggPLru2y2NXe80ISEhJMYmKi3TGUUqrCEJG1xphiW954zfCOUkqp4mnRV0opH6JFXymlfIgWfaWU8iFa9JVSyodo0VdKKR+iRV8ppXyIFn2llPIhHnlyloikAvsKPBwBHLchjrs8PR94fkZPzweasTR4ej7w/IyF5WtgjCm206RHFv3CiEiiO2eb2cXT84HnZ/T0fKAZS4On5wPPz3gp+XR4RymlfIgWfaWU8iEVqehPtjtAMTw9H3h+Rk/PB5qxNHh6PvD8jBedr8KM6SullLp0FWlLXyml1CXSoq+UUj7Eo4q+iPQXke0ikiwi4wuZ3kxEVohItoj8xUMzjhCRDc7bchFp62H5BjqzJYlIooj0LM987mR0ma+TiOSLyE3lmc+57uLex74icsb5PiaJyD89KZ9LxiQR2Swii8sznzsZReQxl/dvk/N3XcOD8oWJyHcist75Ht5ZXtlKkLG6iHzj/D+9WkRaFfukxhiPuAH+wC4gHggC1gMtCsxTC+gEvAD8xUMzdgeqO38eAKzysHxV+N++nDbANk97D13mWwTMBW7ytIxAX2BOef8NliBfOLAFqO+8X8vTMhaY/zpgkSflA/4GvOz8ORI4CQR5WMb/AE85f24GLCzueT1pS78zkGyM2W2MyQGmAwNdZzDGHDPGrAFy7QiIexmXG2NOOe+uBKI9LF+6cf6FAJWB8t6TX2xGpweAr4Fj5RnOyd2MdnEn33BgpjFmP1j/dzwwo6thwLRySWZxJ58BqoqIYG0snQTyPCxjC2AhgDFmGxArIrUv9KSeVPSjgAMu91Ocj3mSkma8C5hXpol+z618IjJYRLYB3wOjyinbecVmFJEoYDAwqRxzuXL399zN+dV/noi0LJ9ogHv5mgDVReQXEVkrIiPLLZ3F7f8rIhIK9Mf6kC8v7uR7C2gOHAI2Ag8ZYxzlEw9wL+N64AYAEekMNKCYDU1PKvpSyGOedjyp2xlF5DKsov9EmSYqsNpCHvtDPmPMN8aYZsAg4LmyDlWAOxknAE8YY/LLPk6h3Mm4DqvXSVvgTeDbsg7lwp18AUBH4E/A1cA/RKRJWQdzUZL/z9cBy4wxJ8swT0Hu5LsaSALqAe2At0SkWtnG+h13Mr6E9eGehPXt+DeK+TYSUCrRSkcKEONyPxrrE9aTuJVRRNoA7wMDjDEnyikblPA9NMYsEZGGIhJhjCmv5lLuZEwAplvfqokArhGRPGPMt+WS0I2MxpizLj/PFZF3yvF9dOc9TAGOG2MygAwRWQK0BXaUQ77z63f3b3Eo5Tu0A+7luxN4yTkcmiwie7DGzVeXT0S3/w7vBHAOQ+1x3opWXjsl3NhpEQDsBuL4306LlkXM+zT27MgtNiNQH0gGuntovkb8b0duB+Dg+fuekrHA/FMp/x257ryPdVzex87A/vJ6H93M1xxrrDcACAU2Aa086T10zheGNVZe2QN/xxOBp50/13b+X4nwsIzhOHcuA/cAHxf3vB6zpW+MyROR+4EfsfZaf2iM2SwiY53TJ4lIHSARqAY4RORhrL3ZZ4t63vLOCPwTqAm849xSzTPl1K3PzXw3AiNFJBc4B9xinH8xHpTRVm5mvAm4V0TysN7HoeX1PrqTzxizVUR+ADYADuB9Y8ym8sjnbkbnrIOBn4z1jaTcuJnvOWCqiGzEGmp5wpTfN2J3MzYHPhaRfKyjte4q7nm1DYNSSvkQT9qRq5RSqoxp0VdKKR+iRV8ppXyIFn2llPIhWvSVUsqHaNFXSikfokVfKaV8yP8DhPs8f95Y7wIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from math import *\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "#Nb de points intérieurs\n", "N=7\n", "h=1/(N+1)\n", "d=np.ones(N-1)\n", "A=(2*np.eye(N,N)-np.diag(d,1)-np.diag(d,-1))/h\n", "\n", "#Exemple de résolution\n", "f=h*np.ones(N)\n", "U=np.linalg.solve(A,f)\n", "\n", "#Tracé\n", "x=np.linspace(h,h*N,N)\n", "plt.plot(x,U)\n", "\n", "M=100000\n", "xref=np.linspace(h,h*N,M)\n", "y=(xref-xref**2)/2\n", "plt.plot(xref,y,'--')\n", "plt.title(\"résolution approchée de -u''=f\")" ] }, { "cell_type": "code", "execution_count": 37, "id": "b41cf0c1", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, \"comparaison entre u et f quand u-eps*u''=f\")" ] }, "execution_count": 37, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEICAYAAABPgw/pAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAr+UlEQVR4nO3deZgc1Xnv8e87+66RRjPaN4QACSyxiNXYYDsxmxNsx3bAOw7GJHaSm5vr4CROvMaxr+Mk1wGMsYMdnMR4IwRsFi8xYMwmyQgEEgjtowVpNJp97+n3/lHVoqdnRtMjTXd19/w+z9PPdPeprnqruvqtM6dO1TF3R0RE8l9R1AGIiMjUUEIXESkQSugiIgVCCV1EpEAooYuIFAgldBGRAqGEXuDM7AEz+0DUcUwXZvZaM3vZzLrN7K1RxzMeM7vUzPZGHUcuMrNKM7vPzDrM7AdRxzMZSugFzt2vcPd/izqOdOR6kjGzb5vZ5yeY7LPAze5e4+73ZCGsgmZmS83s05OcflfS611mtnSSi30HMAdocPd3TvKzkVJCz2NmVhJ1DNmWB+u8BHgh6iDynZldYGZ/DZSEr19vZn+VpcUvAba6eyxLy5s67q5H+AAWAXcDLUArQU0LggPfJ4HdwCHgTmBGWLYUcOA6oBloA24EzgWeA9oT8wmn/yDwa+BfgA7gReBNSeXXAVuALmAH8JGkskuBvcBNwCvAd4CZwI/DmNvC5wuTPvMwcH34/GTgkXC5h4HvJU13EbAuLFsHXJQyj8+FcXcBPwVmH2M7vgXYGK7748DqpLJdwP8Jt00H8D2gAqgG+oA40B0+5gOfBn4I/DvQCVwPzAD+FTgA7AM+DxSPE8u3gc+nbsNjxH4a8DPgCPAS8K7w/RuAIWAwjO2+MT67PYy/L5ymfIxpzgJ+E27H7wF3JeIL943HUqZ34OTw+VXAM+F2aAY+nTTd0nDaDwB7wu/3r5PKK8Nt0QZsBj4+3nZImlfJWPvRZLZb0ndwW1jeRbAPLgnLDPgngt9VR7hfnBGWXR1+5kXg74GapH3ot5Lm/2ng35Ni35Wyvy2dRA74TPgdD4Xf4R9EnZcm84g8gFx5AMXAs+HOVU2QZC4Oyz4EbANOAmoIkv53knYgD3fYCuDNQD9wD9AELAh31kvC6T8IxIA/A0qB3w935Flh+VXA8nBHvwToBc4Oyy4NP/sloDz8kTYAvwdUAbXAD4B7ktbr6A8R+C7w1wQHqOT1m0XwQ38fQY3o2vB1Q9I8tgOnhMt8GPjiONvx7HB9zw+36QfCH1V5WL4LeJogWc8iOHjdmLR+e1Pm9+nwx/XWMO7KcNt+PfyemsL5fWSceL5Nmgk9nF8zwUG1JFyXw8DpY81rnHnsIinZpJSVEVQKEt/9O8J1SzehXwq8JtwOq4GDwFtT9sNvhNtoDTAArAzLvwj8Ktzmi4Dnj7EdEvNKK6Gnud26gNcT7Lf/L7GewGXABqCeYJ9fCcwLy36XVxP654HqsbYxSQl9gu/mEwSVjDEfk51fLj7U5PKq8wiSzMfdvcfd+939sbDsPcA/uvsOd+8G/hK4JuXf/8+Fn/kp0AN8190Pufs+gh/SWUnTHgL+2d2H3P17BDWaqwDc/Sfuvt0DjxDUhl+X9Nk48Cl3H3D3PndvdfcfuXuvu3cBf0dwIBjLEMG/k/NT1u8q4GV3/467x9z9uwQ/ot9J+uy33H2ru/cB3wfOHGcZHwa+7u5PufuwB+33A8AFSdN81d33u/sR4L5jzCvhCXe/x93jQB1wBfC/wu/pEMFB+JoJ5pGOtxDU7r4VboffAD8iSLxT4QKCRJ747n9I8N9QWtz9YXff5O5xd3+O4ACd+l1/JtwvniWooKwJ338X8HfufsTdm4GvnvDavCqd7fYTd3/U3QcIKhUXmtkign2ylqCGb+6+xd0PmNkFBAetPyT4L+anwJ+eSJDu/kV3rx/vcSLzzhVK6K9aBOz2sdvN5hPUrBJ2E9RE5iS9dzDped8Yr2uSXu/zsCqQNL/5AGZ2hZk9aWZHzKwduBKYnTRti7v3J16YWZWZfd3MdptZJ/AoUG9mxWOsx18Q1IKeNrMXzOxD46xfIqYFSa9fSXrem7I+yZYAf25m7YkHwbadfxzzSmhOmX8pcCBp/l8nqKmfqCXA+SmxvweYOwXzhmAbjPXdp8XMzjezX5pZi5l1EDTtzU6ZbLxtO5+R2zHt5Y4RxwthL55uM3sd6W23o8sOK0VHCCoW/wPcDNwCHDSz282szt2fdPfPE/xHSngw+MLxxjxdKKG/qhlYPM5Jt/0EO23CYoId7eAY06ZjgZlZyvz2m1k5Qc3mH4A5Ya3hfoIknJB6e8w/B04Fznf3OoJ/a0n5TPBB91fc/cPuPh/4CHCrmZ08xvolYtp3HOvWTFATTK79VIW1/omMd+vP5PebCWr8s5PmX+fup4/z2R6C5qiEYyXnZuCRlNhr3P0PJ4gvXQcY+7sfM1YzS431P4F7gUXuPoOgmW/U93yMZS8aZ7mpesK/Y243dz893C417v4rJt5uJC/bzGoImn72h/P7qrufA5xO0Kz38aRl7XL3T48RX7rf6VFm9ldJB6JRj3TmkeuU0F/1NMFO/0UzqzazCjN7bVj2XeDPzGxZuDN+geCE4vGeBW8C/sTMSs3snQTthvcTtLGWE5zgjJnZFQRt8sdSS/AfQLuZzQI+Nd6EZvZOM1sYvmwjSFDD4bJPMbN3m1mJmf0+sIrgBOtkfQO4MaxNWrgtrzKz2jQ+exBoMLMZ403g7gcI/v3+ipnVmVmRmS03s/GamTYCV5rZrDBB/q9jLP/HBNvhfeF3U2pm55rZyqT4TkpjPcbzBEFF4E/C7fx2gqa+hGeB083sTDOrIGjLTVYLHHH3fjM7D3j3JJb9feAvzWxmuA/88XgTunsLwcH8vWZWHP4nt/wY855ou0HwHVxsZmUEJ9ifcvfmcLrzzayUIFH3E+yTx7KRoMmz1MzWkmaTmLt/IelANOqRzjxynRJ6yN2HCdqMTyboJbCX4IQlwB0EPUoeBXYS7HTj/iDS8BSwguDE0d8B7wjbwruAPyH48bUR/GDvnWBe/0xwEuww8CTw4DGmPRd4KqyN3Av8qbvvdPdWgnbQPyfo3fMXwFvc/fBkV8zd1xO0o98crsM2gpN96Xz2RYKD547wX/f540z6foKD3+ZwGT8E5o0z7XcIEuUuggPB946x/C6CA+g1BLXHV3j1BDQEPWtWhbHdk846pcx/EHg7wfZoI9i/7k4q30rQj/3nwMvAYymz+CPgs2bWBfwtwX6Srs8QNLPsJNgO35lg+g8T1JRbCWrOj483YRrbDYL/Lj5F0NRyDkGTDATnRL5BsD12h8v7hwli+xuCA0xbuF7/OcH004aNbM6TTDOzDxL0Frg46lgkemb2bYLeJp+MOpZMmQ7rmCtUQxcRKRBK6CIiBUJNLiIiBUI1dBGRAhHZjY5mz57tS5cujWrxIiJ5acOGDYfdvXGsssgS+tKlS1m/fn1UixcRyUtmNu5VvmpyEREpEEroIiIFQgldRKRAKKGLiBQIJXQRkQIxYUI3szvM7JCZPT9OuZnZV81sm5k9Z2ZnT32YIiIykXRq6N8GLj9G+RUEdw5cQTDu4tdOPCwREZmsCfuhu/ujZrb0GJNcDdwZjsLypJnVm9m88L7VmfHSg7Bvw8j3ikrg0puC55v/G15J+YeitBJe97+D55t+CC0vjSyvrIcLPxo83/ifcGTnyPKaJjjvw8Hz9d+Czv0jy2csgHM+GDx/6uvQk3Ln2VknwZnXBs8f/xfo7xxZ3nQanPF7wfNH/wFiA8TdGYzFGYzH6Zt1Bu1LLmMwFqdh/VcYHh5mOO7Bw53W+tUcnHMJ8dggp229LRhjEPBw3NhXZp7D3pnnUxTr5TW7vhW8nzReQ3P9BRyoP4vyoQ7W7B15N1IHds+6mIN1Z1A1eJgz9v1wRBnA9tlv4HDNqdT0vcLpB+9JKg+m2Nr4ZlorT6K+bw8rD91Pqs1NV9FeuYiGnm2cevjno8o3zX0r3eVzaep+keVHHh5RZhjPznsXfWWzmNv5HEvbRt/p9ZmF72GopJb57RtY1L4u+cMAbFj0QYaLK1jc9iTzOjaO+vy6pTfgVsyS1keZ0/nCiGEl3IpZv/QGDOOkll/Q0L11xGdjxRVsXHIdACcffJCZPSP3rcHSWjYtfi8AKw7cR13f3hGjVvSVNbBlUXAn51P3/Yjq/pHjqvRUzOWlBW8H4PTmu6gYahtR3lW5iG3zg9EEz9h9J2WxkWM5tNcsZ9fcywB4zc5vUjw8GGyaMIgjtaexZ86bMIM1274GHh+x+VrqV7O/8XVYfIjVO75BkUGRGcVFRpEZ7U3n0jX/tZTF+1my5XaKwveLzSgphviyN1B60sVUD3dS8vRto7Y9p1wGC9dC10FY983R5St/B+athvZm+M2do8vP+L3g99W6HZ69a3T5mmugYTkc3Awv/Nfo8nM+ADMWwv6N8OJPRpefdwPUNELz0/Dyz0aXX/jRIL/segx2PAJ182HtdaOnmwJTcWHRAkYObbU3fG9UQjezGwhq8SxefKwBUyaw7eejv9iS8lcT+ksPjP7iqhpeTegv/NfoL2bmklcT+nPfhx0Pjyyfc8arCX3jf8DelIuiFp3/akJf/y1oeXFk+clvgjOvJTYchye+RnHXqwcEBzZUv54vPDqXzr4h7u76J2o9GDSmLHw8MHwJN8UqAdhafjMlxEfM/rHhy/hsbAblDLKl/OukemT4EF+J1TKLTtaXj/5R/E+sm68NV7HIDvJI2b+OKr9ve4zvDJex0nbzwTHKv7e9hLvjxjn20pjl39pex0PxYS4pepYPlN4xqvzmHU08Fu/nyqIned8Y5V/avphnfAXvLP4V15aMLv/ktlN42RdyXfEveWfJv48q/7OXV7Of2Xy0+Be8reQHo8qv33oendTwFyU/5S3F940qv3br64lRwmdK7ufK4pEHnEFKeNeLwfgaXym9l98uGnkb8yPUcs3mCwG4rfS/OKto5L7T7I1c+/xaAO4s/SFnFY2sjLzoi/n8c6sBuLvsLs607SPKN/gKPr8xGEviobL/YIWNHGjq0fhqPvubYHyKx8rvZD6tI8rvj5/Hp4eCAYU2ln+LOnpHlP9g+BI+FQsGBdpa/o1R+963hy/js7FZ4b53O6lueX4/X4mVh/veLaPKv/zEQb423Bvue18eNQbTvS8PsW1pLaewi9954stH3z862ayTgoTeuQ8e/TKjzFsdJPQjO8cuX3JhkNAPvzR2+Yo3Bwn9lU1jl5/x9iCh79swdvnZ7w8S+p4ngvKFazOW0NO6OVdYQ/+xu58xRtlPgL/3cMBhM/sF8BfuviF12mRr1671Qr9StG9wmOf2tvPs3naebe5gyyud7GntJRYfuc0bqstorC2nsbacuopSaitKqCkvoTZ8XlVWTHlpEWXFxZSXFFGW9CgvKaKsuGhErceMo7WjImNkWREUW/DaUn44o14n/bJGlyV/zsYtS/1s6rSTlbq/jrX7jrVHj7Wfp74z9rzG+Fwa97M7kXmNjmvi2MecVwQxxN2JDQePoXg8+DscJxZ3YsNxhoadWNL7/bE4PQMxegZidA/E6B0cpjt83dY7REvXAIe7B2jtHiD5Z2MG82dUcsqcGs5cNJMzF9dz5sJ6ZlSVjhFVYTGzDe6+dqyyqaih72XkWIULCccKnI4OdPTxwKZXeHhrC0/uaGUwFtRmFs2q5PR5M7j89Lksm13NstnVLJhZyeyackqL1dkoXaMOHmkfH07sQCLRGo47bb2D7GvrY1drDzsPB48tBzp5eGsL7sG+sGZhPW88rYk3nz6H0+bWRR121k1FQr8X+JiZ3QWcD3RktP08B8WG4zz4wit8b10zj207jDuc1FjNe89fwsUrGlizsJ6GmvKJZyQiYyouMmbXlDO7ppw1i+pHlHX1D7FpbwdP7zrCL19q4Z9+vpV//NlWzlhQx7vWLuLtZy+kpjyy21Zl1YRNLmb2XeBSYDbBILmfAkoB3P22cATzmwl6wvQC14XjSh5TITS5DMedH6xv5taHt7PnSC8LZ1by9rMX8razFrBsdnXU4YlMSy1dA9y/6QDfW9fM5gOd1FeV8qHXLuO61y6ltiL/m2SO1eQS2QAX+Z7Qn9rRyqfv28yWA52sWVTPH126nN9eOYeiIv1rL5IrfrOnjVt/uY2fbzlEU205n3zLKn5n9bwTPpcTJSX0KTQQG+ZLD7zEHb/eyYL6Sj551UouP2NuXu8gIoVuY3M7f3PP82za18EbT2viH965hlnVZVGHdVyU0KfIKx39XH/nOp7f18kHLlzCJ65YSWVZcdRhiUgahuPOvz2+iy8+8CKzqsv42nvP5qzFM6MOa9KOldDVvSJNL77Sydtu/TU7W3r4xvvX8pmrz1AyF8kjxUXGhy5ext1/dBFlJUW8+xtP8csXD0Ud1pRSQk/DS690cc3tTxJ35/s3Xshvr5oTdUgicpzOWDCDH/3hRSxvqub6O9fz4POF0ylPCX0Cuw738J5vPkV5SRE/+MhFnD5/RtQhicgJaqwt564bLmTNwhn8yXc38vi2wxN/KA8ooR9DV/8QH/q3dcTd+Y/rz2dxQ1XUIYnIFKkpL+GOD57L0tlVfPjO9Ww71BV1SCdMCX0c7s7Hf/Acu1t7ufU9Z3NyU23UIYnIFKuvKuPfPnQelWXFfOQ7G+geiEUd0glRQh/Ht369iwdfeIW/vOI0LjipIepwRCRD5s2o5KvXnsXOwz387T1jDvuQN5TQx7DrcA//96EXeeNpTfzBxcuiDkdEMuyi5bP52BtXcPcz+/jFloMTfyBHKaGncHdu+tFzlBYX8YW3vUYXDIlMEx97w8mcNreWv/qvTXT1D0UdznFRQk/xk00HeGrnEf7yipXMnVERdTgikiVlJUV86fdWc6hrgFt+uX3iD+QgJfQk/UPDfPGBF1k5r47fP3fRxB8QkYKyZlE9bz9rIXc8tpPmI70TfyDHKKEn+c4Tu9nb1scnr1pJsW6yJTItffyyUykqgi8/9NLEE+cYJfRQ/9AwX390BxefPJvXnjw76nBEJCJzZ1TwgYuW8uPn9rPzcE/U4UyKEnro++ubOdw9wEffcHLUoYhIxK6/+CRKi4v42sPbog5lUpTQCUYc+vojOzhnyUwuOGlW1OGISMQaa8u55txF3P2bfexv74s6nLQpoQM/33KQfe193PD6k9RNUUQAuP51JzHszn8+tSfqUNKmhA78+5N7mD+jgjed1hR1KCKSIxbNquKNpzZx17rmo4O957ppn9B3Hu7hsW2Hufa8xZQUT/vNISJJ3nvBEg53D/DTza9EHUpapn0Gu+vpPZQUGb9/nvqdi8hIrz+lkUWzKvnu0/nR7DKtE3o87tz77H4uOaWRplpdFSoiIxUXGW87ayGPb2/lUGd/1OFMaFon9HW7jnCgo5/fPXN+1KGISI763TXzcYf7nsv9kY2mdUL/72f3U1larCHlRGRcJzfVcMaCOu7duC/qUCY0bRP60HCc+zcd4M2nz6GqrCTqcEQkh129ZgHP7u1gV45fOTptE/q6XUdo7x3iytfMizoUEclxV64O8sTPNuf2vdKnbUL/xZZDlJUU8boVum+LiBzbgvpKTptby89zfPCLaZnQ3Z1fbDnIRcsb1NwiImn5rZVzWL+7jY7e3B38Ylom9B2He9jV2qsrQ0UkbW9c2cRw3Hl466GoQxnXtEzoiTED37hSvVtEJD1nLqynobqMX2xRQs8pv3r5MCuaalhQXxl1KCKSJ4qKjEtOaeTX2w4Tj3vU4Yxp2iX0wVic9bvaNIiFiEzahcsbaO0ZZOuhrqhDGdO0S+jP7m2nb2iYC05qiDoUEckzFy4P8sYT21sjjmRsaSV0M7vczF4ys21m9okxymeY2X1m9qyZvWBm1019qFPj8W2tmKGBLERk0hbOrGLxrCoez9eEbmbFwC3AFcAq4FozW5Uy2UeBze6+BrgU+IqZlU1xrFPiiR2HWTWvjvqqnAxPRHLcRcsbeGpHK8M52I6eTg39PGCbu+9w90HgLuDqlGkcqLVguJ8a4AgQm9JIp0D/0DC/2d3ORcvV3CIix+fC5Q109sfYvL8z6lBGSSehLwCak17vDd9LdjOwEtgPbAL+1N1HDfFhZjeY2XozW9/S0nKcIR+/Tfs6GByOc+5SNbeIyPE5f1lQIVy/+0jEkYyWTkIfa5DN1P81LgM2AvOBM4Gbzaxu1Ifcb3f3te6+trGxcZKhnrhn9rQBcNbimVlftogUhrkzKpg3o4Jn9rRHHcoo6ST0vUDycD4LCWriya4D7vbANmAncNrUhDh1Nja3s2hWJY215VGHIiJ57KzF9TzT3BZ1GKOkk9DXASvMbFl4ovMa4N6UafYAbwIwsznAqcCOqQx0Kjyzp52zFql2LiIn5qxFM2k+0kdL10DUoYwwYUJ39xjwMeAhYAvwfXd/wcxuNLMbw8k+B1xkZpuAXwA3ufvhTAV9PA509HGgo5+zFtdHHYqI5LlEHtnY3B5pHKnSutWgu98P3J/y3m1Jz/cDb57a0KbWxrC9S+3nInKizlgwg5Ii45k9bTk14tm0uVL0meZ2ykqKWDVv1LlaEZFJqSgtZuW8upw7MTptEvrz+zpYObeWspJps8oikkGrF87g+f0duOfOBUbTIru5O5sPdLJq/oyoQxGRArFqfh1d/TH2tvVFHcpR0yKhH+jop713iFXzaqMORUQKxMqw+Xbzgdy5YnRaJPTEJbqr5qv9XESmxmlzazEjp24BMC0S+pYDnZjBqXOV0EVkalSVlbBsdjVbVEPPrs0HOlnaUE1NuQaEFpGps2penZpcsm3zgU51VxSRKbdqfh172/ro6BuKOhRgGiT0rv4hdrf2slInREVkiiVOjOZKs0vBJ/SXD3UDaj8Xkam3MswrWw/mxhijBZ/Qt4UJfUVTTcSRiEihmVNXTk15ydE8E7WCT+jbD3VTVlzEwpmVUYciIgXGzFjeWM32FiX0rNh2qJtls6spKS74VRWRCCxvqlENPVu2t3RzsppbRCRDTm6q4WDnAJ390fd0KeiE3j80zJ4jvSxXQheRDDm5Mcgv23Ogll7QCX1Xaw9xRzV0EcmYRH7Z3tITcSQFntAT7VrLG6sjjkRECtXiWVWUFltOtKMXfEI3g+WNqqGLSGaUFBextKFaCT3Ttrf0sKC+korS4qhDEZECdnJTDTtyoOtiQSf03a09LJut5hYRyawlDdU0t/UyHI929KICT+i9LJ5VFXUYIlLgljRUMTTs7G+PdvSigk3oHb1DdPQNsaRBCV1EMmtJWHHcc6Q30jgKNqEnNuziWWpyEZHMWtyghJ5Ru48EfUJVQxeRTJs3o5LSYmN3qxJ6RiQ2rNrQRSTTiouMRTOr2HMk2ouLCjah72ntZXZNOdUadk5EsmBxQ5Vq6Jmy+0gPi2fplrkikh2LZ1Wxp7UX9+i6LhZsQt/T2suSBp0QFZHsWDyriq6BGG290d11sSAT+kBsmAOd/Wo/F5GsSVQgd7dG145ekAl9b1sf7urhIiLZsyQHui4WZELf1xZcrbVwphK6iGRHYpjLfRFeLVqQCT1x+e38+oqIIxGR6aKqrIT6qtJIL/9PK6Gb2eVm9pKZbTOzT4wzzaVmttHMXjCzR6Y2zMnZ395HkcGcOiV0Ecme+TMq2d/eH9nyJ+ykbWbFwC3AbwN7gXVmdq+7b06aph64Fbjc3feYWVOG4k3L/o5+5tRVUKqBoUUki+bXV7K3Lbfb0M8Dtrn7DncfBO4Crk6Z5t3A3e6+B8DdD01tmJOzv72PeTNUOxeR7JpfX5HzTS4LgOak13vD95KdAsw0s4fNbIOZvX+sGZnZDWa23szWt7S0HF/Eadjf3sf8el1UJCLZNb++ks7+GF390fRFTyeh2xjvpV4KVQKcA1wFXAb8jZmdMupD7re7+1p3X9vY2DjpYNPh7uzv6GeBErqIZFmiInmgI5p29HQS+l5gUdLrhcD+MaZ50N173P0w8CiwZmpCnJzWnkEGY3E1uYhI1s0P805UzS7pJPR1wAozW2ZmZcA1wL0p0/w38DozKzGzKuB8YMvUhpqeV7ssqoYuItmVyDtR9XSZsJeLu8fM7GPAQ0AxcIe7v2BmN4blt7n7FjN7EHgOiAPfdPfnMxn4eJTQRSQqTbXlFBdZZDX0tO4t6+73A/envHdbyusvA1+eutCOT+LIqIQuItlWUlzEnNpy9nfkbpNLXtnf3kdFaREzq0qjDkVEpqH59ZU53YaeV/Z3BF0WzcbqnCMiklnz6ysju59LwSX0Vzr6matL/kUkIvNmVHCwcyCSgS4KLqEf7BzQPVxEJDKNteUMxuJ09GX/4qKCSujuTkvXAE215VGHIiLTVFNYoTzUNZD1ZRdUQu/oG2JwOH50g4qIZNucsEJ5qFMJ/YQcDDegaugiEpVEhfJgZ/YvLiqohH6oK9iAakMXkagkKpRqcjlBh1RDF5GIVZeXUFNecrSCmU0FldAPhhuwqU4JXUSi01Rbrjb0E3Woc4Da8hKqytK6o4GISEY01ZWrhn6iDnX106jauYhErKm24mgnjWwqrITeOcCcWp0QFZFozQlr6Nm+WrSwEnrXgNrPRSRyTbUV9A/F6RqIZXW5BZPQ3Z2Dnf3qsigikUtULA9luS96wST0zv4YA7G4uiyKSOSawqbfbPd0KZiEnjgSNiqhi0jEEjX0g1nu6VIwCb2lK3FRkZpcRCRaTRHdz6VgEvrhnkEAZteURRyJiEx3NeUllJUUcSTMS9lSMAn9SHdwJGyoUZOLiETLzJhdXUarEvrxae0ZpMigvlJjiYpI9BpqymntVpPLcTncPcis6jKKijSWqIhEb5Zq6MfvSM8ADdVqbhGR3NBQU0ZrtxL6cWntHqRBJ0RFJEfMrimntSe7g0UXTkLvCZpcRERywazqMvqH4vQODmdtmYWT0LsHmK0eLiKSIxrCCmY2m10KIqEPxuJ09seObkARkaglKpitPdnr6VIQCT3ReV990EUkVyTO6amGPkmJI6Da0EUkVyTykWrok5Q4AuqyfxHJFYlu1IdVQ5+cxBFQTS4ikisqy4qpLivO6v1c0kroZna5mb1kZtvM7BPHmO5cMxs2s3dMXYgTS9TQ1eQiIrlkVk1ZVi//nzChm1kxcAtwBbAKuNbMVo0z3ZeAh6Y6yIm09gxSWmzUVZRke9EiIuNqqC7P6uX/6dTQzwO2ufsOdx8E7gKuHmO6PwZ+BByawvjS0todXPZvpvu4iEjumF1TlnNt6AuA5qTXe8P3jjKzBcDbgNumLrT0HenRZf8iknsaqss5kmO9XMaq9qbenOCfgZvc/ZjXuJrZDWa23szWt7S0pBnixHTZv4jkolk1ZRzpGcza/VzSSeh7gUVJrxcC+1OmWQvcZWa7gHcAt5rZW1Nn5O63u/tad1/b2Nh4fBGPob13iPoqJXQRyS31laUMDTs9WbqfSzpnEdcBK8xsGbAPuAZ4d/IE7r4s8dzMvg382N3vmbowj62td5CZVRrYQkRyy8ywotneO0hNeeY7bUxYQ3f3GPAxgt4rW4Dvu/sLZnajmd2Y6QAnEo87HX1DGqlIRHLOjLCi2d47lJXlpXXIcPf7gftT3hvzBKi7f/DEw0pfZ/8Q7qjJRURyzqs19Owk9Ly/UjSxoerV5CIiOSaRl9p6s9N1Me8TemJDzVQNXURyTCKht/ephp6WxIaaoRq6iOSY+sqwySVLV4vmf0JXDV1EclRZSRHVZcWqoafraBu6ermISA6qrypTG3q62nqHMIM6JXQRyUH1VaV0qJdLetp7B6mrKKW4SDfmEpHcU19Vqhp6uoLL/lU7F5HcVF9Vpjb0dLX1DuqiIhHJWTOrSnVhUbp02b+I5LL6yjLaeweJxzN/x8W8T+i6MZeI5LL6qlLiDl0DsYwvK+8Tum6dKyK5rD7pjouZltcJPTYcp6s/ppOiIpKzZmbxjot5ndA7+nRRkYjktmzeoCuvE3pbeMSbqeHnRCRHJZpcOrLQdTGvE3qiTWqGaugikqMSLQhtWbhBV14n9KNNLjopKiI5KlHhzMbFRXmd0Dv7w1vnqoYuIjmqpDi442JXv7otHlNnX7CBaisyP/iqiMjxqqsspVM19GNLbCAldBHJZXUVpUdbFDIpvxN6/xAVpUWUlxRHHYqIyLjqKkuOtihkUl4n9K7+GHUVaj8XkdxWV1FK14Bq6MfU2T+kgS1EJOfVVqiGPqHOvpjaz0Uk59VVqg19Qp39Q2pyEZGcV1cR9HJxz+wtdPM7ofepyUVEcl9dZQlxh57B4YwuJ78Ten+MOjW5iEiOqw1bEjLdFz1vE7q706WToiKSBxJNw5m+WjRvE3r/UJyhYVcbuojkvLrKoCUh0ydG8zahJzaMermISK6rU5PLsSU2jJpcRCTXJfKUaujjSGwYnRQVkVyXaEnI9MVFaSV0M7vczF4ys21m9okxyt9jZs+Fj8fNbM3UhzpSZ3hyQTV0Ecl1iYTeFXUN3cyKgVuAK4BVwLVmtiplsp3AJe6+GvgccPtUB5rqaJOLToqKSI4rLymmorToaEU0U9KpoZ8HbHP3He4+CNwFXJ08gbs/7u5t4csngYVTG+ZoR2voanIRkTyQuFo0k9JJ6AuA5qTXe8P3xvMHwANjFZjZDWa23szWt7S0pB/lGHRSVETySW1FSU6cFLUx3hvzhgRm9gaChH7TWOXufru7r3X3tY2NjelHOYbO/iHKiosoL8nb87oiMo0EoxZF3+SyF1iU9HohsD91IjNbDXwTuNrdW6cmvPF19ceoqyzBbKzjjYhIbqmrKI3+pCiwDlhhZsvMrAy4Brg3eQIzWwzcDbzP3bdOfZijdfbpTosikj+CW+hmtoY+4RlFd4+Z2ceAh4Bi4A53f8HMbgzLbwP+FmgAbg1rzDF3X5u5sIOTorVqPxeRPFFXUZLxk6JpdRFx9/uB+1Peuy3p+fXA9VMb2rEFNXT1cBGR/FAbDhTt7hlrKs7bM4pd/UPUlCuhi0h+qK0oYWjYGYjFM7aMvE3oPQPDSugikjcSV4t2D2SuHT1vE3r3QIwaNbmISJ6oLgsTegZPjOZlQo/HnZ7BmGroIpI3alRDH1vv0DDuKKGLSN5I5Csl9BQ94QapVkIXkTyRSOg9SugjJcbl02hFIpIvqlVDH1viCKcmFxHJF+rlMo5uNbmISJ45WkNXL5eRulVDF5E8U1VajJna0EdJHOGU0EUkXxQVGdVlJXQpoY/UMxgmdJ0UFZE8UlNeohp6qi7V0EUkD1WXF+ukaKqegRglRabRikQkr9RUlNI9MJyx+edlRuweiFFdrtGKRCS/1JaX0J3BUYvyNqGruUVE8k11eTE9qqGP1N2vhC4i+aemvFRt6Kl6BnXrXBHJPzU6KTpad39MV4mKSN6pqSiheyCGu2dk/vmZ0Adi1Cqhi0ieqS4vYTieuWHo8jahV5cXRx2GiMikJCqiXRm6n0teJvRgPNHSqMMQEZmU6gzfEz3vEno87mG3RdXQRSS/ZHrUorxL6L1DQR9O9XIRkXyjhJ4icadF9XIRkXxzdKBotaEHdC90EclXR9vQB5XQASV0Eclf6uWSQuOJiki+Ui+XFF1qQxeRPFVVFgxDp5OiocbaMq58zVxm15RHHYqIyKSYGb+7Zj7LG2syM/9M3VNgImvXrvX169dHsmwRkXxlZhvcfe1YZXlXQxcRkbGlldDN7HIze8nMtpnZJ8YoNzP7alj+nJmdPfWhiojIsUyY0M2sGLgFuAJYBVxrZqtSJrsCWBE+bgC+NsVxiojIBNKpoZ8HbHP3He4+CNwFXJ0yzdXAnR54Eqg3s3lTHKuIiBxDOgl9AdCc9Hpv+N5kp8HMbjCz9Wa2vqWlZbKxiojIMaST0G2M91K7xqQzDe5+u7uvdfe1jY2N6cQnIiJpSieh7wUWJb1eCOw/jmlERCSD0kno64AVZrbMzMqAa4B7U6a5F3h/2NvlAqDD3Q9McawiInIME14/7+4xM/sY8BBQDNzh7i+Y2Y1h+W3A/cCVwDagF7huovlu2LDhsJntPs64ZwOHj/Oz+UrrPD1onaeHE1nnJeMVRHal6Ikws/XjXSlVqLTO04PWeXrI1DrrSlERkQKhhC4iUiDyNaHfHnUAEdA6Tw9a5+khI+ucl23oIiIyWr7W0EVEJIUSuohIgcjphD4db9ubxjq/J1zX58zscTNbE0WcU2midU6a7lwzGzazd2QzvkxIZ53N7FIz22hmL5jZI9mOcaqlsW/PMLP7zOzZcJ0nvJ4ll5nZHWZ2yMyeH6d86vOXu+fkg+Aipu3ASUAZ8CywKmWaK4EHCO4lcwHwVNRxZ2GdLwJmhs+vmA7rnDTd/xBcxPaOqOPOwvdcD2wGFoevm6KOOwvr/FfAl8LnjcARoCzq2E9gnV8PnA08P075lOevXK6hT8fb9k64zu7+uLu3hS+fJLhvTj5L53sG+GPgR8ChbAaXIems87uBu919D4C75/t6p7PODtSamQE1BAk9M6MpZ4G7P0qwDuOZ8vyVywl9ym7bm0cmuz5/QHCEz2cTrrOZLQDeBtyWxbgyKZ3v+RRgppk9bGYbzOz9WYsuM9JZ55uBlQQ39tsE/Km7x7MTXiSmPH9NeC+XCE3ZbXvzSNrrY2ZvIEjoF2c0osxLZ53/GbjJ3YeDylveS2edS4BzgDcBlcATZvaku2/NdHAZks46XwZsBN4ILAd+Zma/cvfODMcWlSnPX7mc0KfjbXvTWh8zWw18E7jC3VuzFFumpLPOa4G7wmQ+G7jSzGLufk9WIpx66e7bh929B+gxs0eBNUC+JvR01vk64IseNDBvM7OdwGnA09kJMeumPH/lcpPLdLxt74TrbGaLgbuB9+VxbS3ZhOvs7svcfam7LwV+CPxRHidzSG/f/m/gdWZWYmZVwPnAlizHOZXSWec9BP+RYGZzgFOBHVmNMrumPH/lbA3dM3Tb3lyW5jr/LdAA3BrWWGOex3eqS3OdC0o66+zuW8zsQeA5IA58093H7P6WD9L8nj8HfNvMNhE0R9zk7nl7W10z+y5wKTDbzPYCnwJKIXP5S5f+i4gUiFxuchERkUlQQhcRKRBK6CIiBUIJXUSkQCihi4gUCCV0EZECoYQuIlIg/j8w5ZA+kFmqhwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from math import *\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "#Nb de points intérieurs\n", "N=1000\n", "h=1/(N+1)\n", "d=np.ones(N-1)\n", "A=(2*np.eye(N,N)-np.diag(d,1)-np.diag(d,-1))/h\n", "B=h*(2*np.eye(N,N)/3+np.diag(d,1)/6+np.diag(d,-1)/6)\n", "\n", "#Exemple de résolution\n", "eps=0.001\n", "f=h*np.ones(N)\n", "U=np.linalg.solve(B+eps*A,f)\n", "\n", "#Tracé\n", "x=np.linspace(h,h*N,N)\n", "plt.plot(x,U)\n", "\n", "M=100000\n", "xref=np.linspace(h,h*N,M)\n", "y=np.ones(M)\n", "plt.plot(xref,y,'--')\n", "plt.title(\"comparaison entre u et f quand u-eps*u''=f\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.10.12" } }, "nbformat": 4, "nbformat_minor": 5 }