Mercurial > repos > public > sbplib

view +time/+rk/rungekutta_4.m @ 888:8732d6bd9890 feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add general Runge-Kutta class - Add a general Runge-Kutta class which time integrates the solution based on coefficients obtained from a Butcher tableau - Add butcher tableau which returns coefficents for the specified Runge-Kutta method - Remove RungKutta4proper, since obsolete
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 15 Nov 2018 17:10:01 -0800
parents 50d5a3843099
children
line wrap: on
line source

% Takes one time step of size k using the rungekutta method
% starting from v_0 and where the function F(v,t) gives the
% time derivatives.
function v = rungekutta_4(v, t , k, F)
    k1 = F(v         ,t      );
    k2 = F(v+0.5*k*k1,t+0.5*k);
    k3 = F(v+0.5*k*k2,t+0.5*k);
    k4 = F(v+    k*k3,t+    k);
    v = v + (1/6)*(k1+2*(k2+k3)+k4)*k;
end