Mercurial > repos > public > sbplib
diff +time/Rungekutta4proper.m @ 41:910a05dcdfdf
Added Rk4 as it should be.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Thu, 05 Nov 2015 16:06:06 -0800 |
parents | +time/Rungekutta4.m@b18d3d201a71 |
children | e1a05acc1b5d b5e5b195da1e |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/+time/Rungekutta4proper.m Thu Nov 05 16:06:06 2015 -0800 @@ -0,0 +1,41 @@ +classdef Rungekutta4proper < time.Timestepper + properties + F + k + t + v + m + n + end + + + methods + function obj = Rungekutta4proper(F, k, t0, v0) + obj.F = F; + obj.k = k; + obj.t = t0; + obj.v = v0; + obj.m = length(v0); + obj.n = 0; + end + + function [v,t] = getV(obj) + v = obj.v; + t = obj.t; + end + + function obj = step(obj) + obj.v = time.rk4.rungekutta_4(obj.v, obj.t, obj.k, obj.F); + obj.t = obj.t + obj.k; + obj.n = obj.n + 1; + end + end + + + methods (Static) + function k = getTimeStep(lambda) + k = rk4.get_rk4_time_step(lambda); + end + end + +end \ No newline at end of file