Mercurial > repos > public > sbplib

comparison +time/Rungekutta4.m @ 991:a99f00896b8e feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make Rungekutta4SecondOrder native second order
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 09 Jan 2019 10:17:00 +0100
parents 0585a2ee7ee7
children
comparison
equal deleted inserted replaced
990:1066bb31bc95 991:a99f00896b8e
24 v = obj.v; 24 v = obj.v;
25 t = obj.t; 25 t = obj.t;
26 end 26 end
27 27
28 function obj = step(obj) 28 function obj = step(obj)
29 t = obj.t;
29 v = obj.v; 30 v = obj.v;
30 dt = obj.dt; 31 dt = obj.dt;
31 32
32 k1 = obj.F(t, v); 33 k1 = obj.F(t, v);
33 k2 = obj.F(t + 0.5*dt, v + 0.5*dt*k1); 34 k2 = obj.F(t + 0.5*dt, v + 0.5*dt*k1);