Mercurial > repos > public > sbplib
annotate +time/Rungekutta4.m @ 31:d1f9dd55a2b0
Fixed bug in step() for CdiffNonlin.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Mon, 28 Sep 2015 21:56:02 +0200 |
| parents | b18d3d201a71 |
| children | ef5c9870f386 |
| rev | line source |
|---|---|
| 0 | 1 classdef Rungekutta4 < time.Timestepper |
| 2 properties | |
| 3 D | |
| 4 S | |
| 5 F | |
| 6 k | |
| 7 t | |
| 8 v | |
| 9 m | |
| 10 n | |
| 11 end | |
| 12 | |
| 13 | |
| 14 methods | |
| 15 function obj = Rungekutta4(D, S, k, t0, v0) | |
| 16 obj.D = D; | |
| 17 obj.k = k; | |
| 18 obj.t = t0; | |
| 19 obj.v = v0; | |
| 20 obj.m = length(v0); | |
|
13
b18d3d201a71
Fixed initialization of step counter in timesteppers.
Jonatan Werpers <jonatan@werpers.com>
parents:
0
diff
changeset
|
21 obj.n = 0; |
| 0 | 22 |
| 23 if S == 0 | |
| 24 obj.S = zeros(obj.m,1); | |
| 25 else | |
| 26 obj.S = S; | |
| 27 end | |
| 28 | |
| 29 obj.F = @(v,t)(obj.D*v + obj.S); | |
| 30 end | |
| 31 | |
| 32 function [v,t] = getV(obj) | |
| 33 v = obj.v; | |
| 34 t = obj.t; | |
| 35 end | |
| 36 | |
| 37 function obj = step(obj) | |
| 38 obj.v = time.rk4.rungekutta_4(obj.v, obj.t, obj.k, obj.F); | |
| 39 obj.t = obj.t + obj.k; | |
| 40 obj.n = obj.n + 1; | |
| 41 end | |
| 42 end | |
| 43 | |
| 44 | |
| 45 methods (Static) | |
| 46 function k = getTimeStep(lambda) | |
| 47 k = rk4.get_rk4_time_step(lambda); | |
| 48 end | |
| 49 end | |
| 50 | |
| 51 end |