Mercurial > repos > public > sbplib
diff +time/Cdiff.m @ 978:1a30dbe99c7c
Refactor CdiffImplicit to take input arguments in the right order
author | Jonatan Werpers <jonatan@werpers.com> |
---|---|
date | Mon, 07 Jan 2019 16:15:49 +0100 |
parents | c9009d5a3101 |
children | 7a5e770974ed |
line wrap: on
line diff
--- a/+time/Cdiff.m Mon Jan 07 15:56:35 2019 +0100 +++ b/+time/Cdiff.m Mon Jan 07 16:15:49 2019 +0100 @@ -5,8 +5,7 @@ G k t - v - v_prev + v, v_prev n end @@ -32,9 +31,9 @@ obj.G = G; % Rewrite as AA*v_(n+1) + BB*v_n + CC*v_(n-1) = G(t_n) - obj.AA = A/k^2 + B/2/k; + obj.AA = A/k^2 + B/(2*k); obj.BB = -2*A/k^2 + C; - obj.CC = A/k^2 - B/2/k; + obj.CC = A/k^2 - B/(2*k); obj.k = k; obj.v_prev = v0; @@ -70,4 +69,4 @@ obj.n = obj.n + 1; end end -end \ No newline at end of file +end