Mercurial > repos > public > sbplib

comparison +time/Rungekutta4.m @ 75:ef5c9870f386

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Cleand time.Rungekutta4 a bit.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 23 Nov 2015 17:32:56 +0100
parents b18d3d201a71
children b5e5b195da1e c6fcee3fcf1b
comparison
equal deleted inserted replaced
74:d0df12c03aac 75:ef5c9870f386
24 obj.S = zeros(obj.m,1); 24 obj.S = zeros(obj.m,1);
25 else 25 else
26 obj.S = S; 26 obj.S = S;
27 end 27 end
28 28
29 obj.F = @(v,t)(obj.D*v + obj.S); 29 if S == 0
30 obj.F = @(v,t)(obj.D*v);
31 else
32 obj.F = @(v,t)(obj.D*v + obj.S);
33 end
30 end 34 end
31 35
32 function [v,t] = getV(obj) 36 function [v,t] = getV(obj)
33 v = obj.v; 37 v = obj.v;
34 t = obj.t; 38 t = obj.t;