Mercurial > repos > public > sbplib

comparison +time/Rungekutta4SecondOrder.m @ 13:b18d3d201a71

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed initialization of step counter in timesteppers.
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 22 Sep 2015 08:41:28 +0200
parents 48b6fb693025
children e7e73173d44d b5e5b195da1e
comparison
equal deleted inserted replaced
12:c71095442c57 13:b18d3d201a71
18 function obj = Rungekutta4SecondOrder(D, E, S, k, t0, v0, v0t) 18 function obj = Rungekutta4SecondOrder(D, E, S, k, t0, v0, v0t)
19 obj.D = D; 19 obj.D = D;
20 obj.E = E; 20 obj.E = E;
21 obj.S = S; 21 obj.S = S;
22 obj.m = length(v0); 22 obj.m = length(v0);
23 obj.n = 0;
23 24
24 I = speye(obj.m); 25 I = speye(obj.m);
25 O = sparse(obj.m,obj.m); 26 O = sparse(obj.m,obj.m);
26 obj.M = [O, I; D, E*I]; % Multiply with I to allow 0 as input. 27 obj.M = [O, I; D, E*I]; % Multiply with I to allow 0 as input.
27 28