From 6db52d799a9382fe3065cd92aa7e75ef7b44dedf Mon Sep 17 00:00:00 2001 From: Ksero Date: Fri, 28 Oct 2016 15:29:18 +0200 Subject: [PATCH] Fix multivariate Runge-Kutta ODE solver typo The multiplication by the time-step was missing in this method (but not the others) --- src/Numerics/OdeSolvers/RungeKutta.cs | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/Numerics/OdeSolvers/RungeKutta.cs b/src/Numerics/OdeSolvers/RungeKutta.cs index bbf82490..77a65d71 100644 --- a/src/Numerics/OdeSolvers/RungeKutta.cs +++ b/src/Numerics/OdeSolvers/RungeKutta.cs @@ -117,7 +117,7 @@ namespace MathNet.Numerics.OdeSolvers { k1 = f(t, y0); k2 = f(t, y0 + k1 + dt); - y[i] = y0 + 0.5 * (k1 + k2); + y[i] = y0 + dt * 0.5 * (k1 + k2); t += dt; y0 = y[i]; } @@ -153,4 +153,4 @@ namespace MathNet.Numerics.OdeSolvers return y; } } -} \ No newline at end of file +}