From 63ea12bc3dc413a34a146b35307aa3bf32bd5641 Mon Sep 17 00:00:00 2001 From: Krishna Vedala <7001608+kvedala@users.noreply.github.com> Date: Wed, 10 Jun 2020 12:50:09 -0400 Subject: [PATCH] update independent variable in the main loop --- numerical_methods/ode_forward_euler.c | 1 + numerical_methods/ode_midpoint_euler.c | 1 + 2 files changed, 2 insertions(+) diff --git a/numerical_methods/ode_forward_euler.c b/numerical_methods/ode_forward_euler.c index 9d98255ea..5df0067b6 100644 --- a/numerical_methods/ode_forward_euler.c +++ b/numerical_methods/ode_forward_euler.c @@ -116,6 +116,7 @@ int main(int argc, char *argv[]) { fprintf(fp, "%.4g,%.4g,%.4g\n", x, y[0], y[1]); // write to file forward_euler(dx, &x, y, dy); // perform integration + x += dx; // update step } while (x <= X_MAX); // till upper limit of independent variable /* end of integration */ diff --git a/numerical_methods/ode_midpoint_euler.c b/numerical_methods/ode_midpoint_euler.c index 6c4087515..5be3c0900 100644 --- a/numerical_methods/ode_midpoint_euler.c +++ b/numerical_methods/ode_midpoint_euler.c @@ -124,6 +124,7 @@ int main(int argc, char *argv[]) { fprintf(fp, "%.4g,%.4g,%.4g\n", x, y[0], y[1]); // write to file midpoint_euler(dx, &x, y, dy); // perform integration + x += dx; // update step } while (x <= X_MAX); // till upper limit of independent variable /* end of integration */