changeset 41:910a05dcdfdf

Added Rk4 as it should be.
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 05 Nov 2015 16:06:06 -0800
parents 54d3ab296ba0
children f121bf58c1b9
files +time/Rungekutta4proper.m
diffstat 1 files changed, 41 insertions(+), 0 deletions(-) [+]
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