Mercurial > repos > public > sbplib
comparison +time/Rungekutta4SecondOrder.m @ 345:7b5ef8b89268 feature/beams
Rungekutta bug fixes.
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Wed, 16 Nov 2016 15:39:16 -0800 |
parents | e7e73173d44d |
children | ae905a11e32c |
comparison
equal
deleted
inserted
replaced
344:61b75d6d4899 | 345:7b5ef8b89268 |
---|---|
32 default_arg('S', @(t)sparse(obj.m, 1) ); | 32 default_arg('S', @(t)sparse(obj.m, 1) ); |
33 | 33 |
34 if ~isa(D, 'function_handle') | 34 if ~isa(D, 'function_handle') |
35 D = @(t)D; | 35 D = @(t)D; |
36 end | 36 end |
37 if ~isa(D, 'function_handle') | 37 if ~isa(E, 'function_handle') |
38 E = @(t)E; | 38 E = @(t)E; |
39 end | 39 end |
40 if ~isa(D, 'function_handle') | 40 if ~isa(S, 'function_handle') |
41 S = @(t)S; | 41 S = @(t)S; |
42 end | 42 end |
43 | 43 |
44 I = speye(obj.m); | 44 I = speye(obj.m); |
45 O = sparse(obj.m,obj.m); | 45 O = sparse(obj.m,obj.m); |