diff +time/+expint/Magnus_mp.m @ 509:566304a5733e feature/quantumTriangles

changed the magnus midpoint rule
author Ylva Rydin <ylva.rydin@telia.com>
date Fri, 02 Jun 2017 17:43:59 +0200
parents 49bc2929dfb1
children 4ef2d2a493f1
line wrap: on
line diff
--- a/+time/+expint/Magnus_mp.m	Fri Jun 02 11:33:04 2017 +0200
+++ b/+time/+expint/Magnus_mp.m	Fri Jun 02 17:43:59 2017 +0200
@@ -5,7 +5,7 @@
 
 if isa(D,'function_handle')
    % v = krylov(k*D(t +k/2*t),v);
-   v = expm(k*D(t +k/2*t))*v;
+   v = expm(k*D(t +k/2))*v;
 else
    %v = krylov(k*D,v);
    v = expm(k*D)*v;