From 52c60121978710a2b4f806d87e764445ee0760dd Mon Sep 17 00:00:00 2001 From: Arthur DANJOU Date: Mon, 24 Mar 2025 16:07:17 +0100 Subject: [PATCH] add: TP1/2 --- M1/Numerical Methods/TP1.ipynb | 75 ++++++++++++++-------------------- M1/Numerical Methods/TP2.ipynb | 28 +++++++------ 2 files changed, 47 insertions(+), 56 deletions(-) diff --git a/M1/Numerical Methods/TP1.ipynb b/M1/Numerical Methods/TP1.ipynb index 714ccca..3594122 100644 --- a/M1/Numerical Methods/TP1.ipynb +++ b/M1/Numerical Methods/TP1.ipynb @@ -170,8 +170,8 @@ { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:26.900678Z", - "start_time": "2025-03-24T13:04:26.896193Z" + "end_time": "2025-03-24T13:14:56.636016Z", + "start_time": "2025-03-24T13:14:56.632263Z" } }, "cell_type": "code", @@ -185,85 +185,85 @@ ], "id": "61588e94b56ea779", "outputs": [], - "execution_count": 22 + "execution_count": 41 }, { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:27.128889Z", - "start_time": "2025-03-24T13:04:27.123995Z" + "end_time": "2025-03-24T13:14:57.036853Z", + "start_time": "2025-03-24T13:14:57.033785Z" } }, "cell_type": "code", "source": [ - "def euler_explicit_sir(N, T, u0, f, args):\n", + "def euler_explicit_sir(N, T, y0, f, args):\n", " t = np.linspace(0, T, N + 1)\n", " y = np.zeros((N + 1, len(y0)))\n", " y[0] = y0\n", " h = T / N\n", " for i in range(N):\n", - " y[i + 1] = y[i] + h * F(y[i], t[i], *args)\n", + " y[i + 1] = y[i] + h * f(y[i], t[i], *args)\n", " return t, y" ], "id": "a32f0facb958a1a0", "outputs": [], - "execution_count": 23 + "execution_count": 42 }, { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:27.315627Z", - "start_time": "2025-03-24T13:04:27.312627Z" + "end_time": "2025-03-24T13:14:57.487407Z", + "start_time": "2025-03-24T13:14:57.484369Z" } }, "cell_type": "code", "source": [ - "def heun_sir(N, T, u0, f, args):\n", + "def heun_sir(N, T, y0, f, args):\n", " t = np.linspace(0, T, N + 1)\n", " y = np.zeros((N + 1, len(y0)))\n", " y[0] = y0\n", " h = T / N\n", " for i in range(N):\n", - " k1 = h * F(y[i], t[i], *args)\n", - " k2 = h * F(y[i] + k1, t[i + 1], *args)\n", + " k1 = h * f(y[i], t[i], *args)\n", + " k2 = h * f(y[i] + k1, t[i + 1], *args)\n", " y[i + 1] = y[i] + 0.5 * (k1 + k2)\n", " return t, y" ], "id": "ca061c36b09f7200", "outputs": [], - "execution_count": 24 + "execution_count": 43 }, { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:27.489907Z", - "start_time": "2025-03-24T13:04:27.485432Z" + "end_time": "2025-03-24T13:14:57.869526Z", + "start_time": "2025-03-24T13:14:57.865654Z" } }, "cell_type": "code", "source": [ - "def runge_kutta_order_4_sir(N, T, u0, f, args):\n", + "def runge_kutta_order_4_sir(N, T, y0, f, args):\n", " t = np.linspace(0, T, N + 1)\n", " y = np.zeros((N + 1, len(y0)))\n", " y[0] = y0\n", " h = T / N\n", " for i in range(N):\n", - " k1 = h * F(y[i], t[i], *args)\n", - " k2 = h * F(y[i] + 0.5 * k1, t[i] + 0.5 * h, *args)\n", - " k3 = h * F(y[i] + 0.5 * k2, t[i] + 0.5 * h, *args)\n", - " k4 = h * F(y[i] + k3, t[i + 1], *args)\n", + " k1 = h * f(y[i], t[i], *args)\n", + " k2 = h * f(y[i] + 0.5 * k1, t[i] + 0.5 * h, *args)\n", + " k3 = h * f(y[i] + 0.5 * k2, t[i] + 0.5 * h, *args)\n", + " k4 = h * f(y[i] + k3, t[i + 1], *args)\n", " y[i + 1] = y[i] + (k1 + 2 * k2 + 2 * k3 + k4) / 6\n", " return t, y" ], "id": "afc8571b14a6b464", "outputs": [], - "execution_count": 25 + "execution_count": 44 }, { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:27.725395Z", - "start_time": "2025-03-24T13:04:27.640769Z" + "end_time": "2025-03-24T13:14:58.376791Z", + "start_time": "2025-03-24T13:14:58.220423Z" } }, "cell_type": "code", @@ -298,13 +298,13 @@ "output_type": "display_data" } ], - "execution_count": 26 + "execution_count": 45 }, { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:28.172102Z", - "start_time": "2025-03-24T13:04:27.807260Z" + "end_time": "2025-03-24T13:22:36.716443Z", + "start_time": "2025-03-24T13:22:36.310338Z" } }, "cell_type": "code", @@ -356,7 +356,7 @@ "output_type": "display_data" } ], - "execution_count": 27 + "execution_count": 48 }, { "metadata": {}, @@ -389,8 +389,8 @@ { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:04:30.150912Z", - "start_time": "2025-03-24T13:04:29.984979Z" + "end_time": "2025-03-24T13:46:24.490267Z", + "start_time": "2025-03-24T13:46:24.315085Z" } }, "cell_type": "code", @@ -444,20 +444,7 @@ "output_type": "display_data" } ], - "execution_count": 29 - }, - { - "metadata": { - "ExecuteTime": { - "end_time": "2025-03-24T12:48:05.902735Z", - "start_time": "2025-03-24T12:48:05.890888Z" - } - }, - "cell_type": "code", - "source": "", - "id": "aeb67c5158f41a0b", - "outputs": [], - "execution_count": null + "execution_count": 53 } ], "metadata": { diff --git a/M1/Numerical Methods/TP2.ipynb b/M1/Numerical Methods/TP2.ipynb index 8bba08b..14c1961 100644 --- a/M1/Numerical Methods/TP2.ipynb +++ b/M1/Numerical Methods/TP2.ipynb @@ -17,13 +17,12 @@ { "metadata": { "ExecuteTime": { - "end_time": "2025-03-24T13:10:01.794956Z", - "start_time": "2025-03-24T13:09:58.777694Z" + "end_time": "2025-03-24T14:45:40.568860Z", + "start_time": "2025-03-24T14:45:37.895550Z" } }, "cell_type": "code", "source": [ - "\n", "import numpy as np\n", "from sklearn.neural_network import MLPClassifier\n", "from sklearn.datasets import make_classification\n", @@ -42,7 +41,7 @@ " y_pred = model.predict(X_test)\n", " accuracies.append(accuracy_score(y_test, y_pred))\n", "\n", - "print(f'Accuracy: {np.mean(accuracies) * 100:.0f}%')\n", + "print(f'Mean Accuracy: {np.mean(accuracies) * 100:.0f}%')\n", "print(f\"Max accuracy: {np.max(accuracies) * 100:.0f}%\")\n", "print(f\"Min accuracy: {np.min(accuracies) * 100:.0f}%\")" ], @@ -52,21 +51,26 @@ "name": "stdout", "output_type": "stream", "text": [ - "Accuracy: 93%\n", - "Max accuracy: 100%\n", - "Min accuracy: 80%\n" + "Mean Accuracy: 92%\n", + "Max accuracy: 96%\n", + "Min accuracy: 86%\n" ] } ], - "execution_count": 21 + "execution_count": 28 }, { - "metadata": {}, + "metadata": { + "ExecuteTime": { + "end_time": "2025-03-24T14:37:53.507776Z", + "start_time": "2025-03-24T14:37:53.505376Z" + } + }, "cell_type": "code", - "outputs": [], - "execution_count": null, "source": "", - "id": "96b6d46883ed5570" + "id": "96b6d46883ed5570", + "outputs": [], + "execution_count": null } ], "metadata": {