--- a/+time/CdiffTimeDep.m	Wed Apr 10 22:22:46 2019 +0200
+++ b/+time/CdiffTimeDep.m	Wed Apr 10 22:40:55 2019 +0200
@@ -3,7 +3,7 @@
         D
         E
         S
-        k
+        dt
         t
         v
         v_prev
@@ -15,8 +15,8 @@
         % Solves u_tt = Du + E(t)u_t + S(t)
         % D, E, S can either all be constants or all be function handles,
         % They can also be omitted by setting them equal to the empty matrix.
-        % CdiffTimeDep(D, E, S, k, t0, n0, v, v_prev)
-        function obj = CdiffTimeDep(D, E, S, k, t0, n0, v, v_prev)
+        % CdiffTimeDep(D, E, S, dt, t0, n0, v, v_prev)
+        function obj = CdiffTimeDep(D, E, S, dt, t0, n0, v, v_prev)
             m = length(v);
             default_arg('E', @(t)sparse(m,m));
             default_arg('S', @(t)sparse(m,1));
@@ -25,7 +25,7 @@
             obj.E = E;
             obj.S = S;
 
-            obj.k = k;
+            obj.dt = dt;
             obj.t = t0;
             obj.n = n0;
             obj.v = v;
@@ -38,13 +38,13 @@
         end
 
         function [vt,t] = getVt(obj)
-            vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u))
+            vt = (obj.v-obj.v_prev)/obj.dt; % Could be improved using u_tt = f(u))
             t = obj.t;
         end
 
         function obj = step(obj)
-            [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E(obj.t), obj.S(obj.t));
-            obj.t = obj.t + obj.k;
+            [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.dt, obj.D, obj.E(obj.t), obj.S(obj.t));
+            obj.t = obj.t + obj.dt;
             obj.n = obj.n + 1;
         end
     end