Mercurial > repos > public > sbplib

comparison +time/Rungekutta4proper.m @ 41:910a05dcdfdf

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
comparison
equal deleted inserted replaced
40:54d3ab296ba0 41:910a05dcdfdf
1 classdef Rungekutta4proper < time.Timestepper
2 properties
3 F
4 k
5 t
6 v
7 m
8 n
9 end
10
11
12 methods
13 function obj = Rungekutta4proper(F, k, t0, v0)
14 obj.F = F;
15 obj.k = k;
16 obj.t = t0;
17 obj.v = v0;
18 obj.m = length(v0);
19 obj.n = 0;
20 end
21
22 function [v,t] = getV(obj)
23 v = obj.v;
24 t = obj.t;
25 end
26
27 function obj = step(obj)
28 obj.v = time.rk4.rungekutta_4(obj.v, obj.t, obj.k, obj.F);
29 obj.t = obj.t + obj.k;
30 obj.n = obj.n + 1;
31 end
32 end
33
34
35 methods (Static)
36 function k = getTimeStep(lambda)
37 k = rk4.get_rk4_time_step(lambda);
38 end
39 end
40
41 end