--- a/+time/+expint/Magnus_mp.m~	Thu Oct 05 18:04:23 2017 +0200
+++ b/+time/+expint/Magnus_mp.m~	Fri Oct 06 13:51:21 2017 +0200
@@ -4,14 +4,18 @@
 function v = Magnus_mp(v,D, t , k,matrixexp,tol)
 
 if isa(D,'function_handle')
- switch matrixexp
-     case 'expm'
-        v = expm(k*D(t +k/2))*v;
-     case 'Arnol'
-    v = time.expint.expm_Arnoldi(-D(t +k/2),v,k,toler,100);
+    Omega = D(t +k/2);
 else
-   %v = krylov(k*D,v);
-   % v = expm(k*D)*v;
+    Omega = D;
 end
 
+switch matrixexp
+    case 'expm'
+        v = expm(k*Omega)*v;
+    case 'Arnoldi'
+        [v, i] = time.expint.expm_Arnoldi(-Omega,v,k,tol,1000);
+    otherwise
+        error('No such matrix exponential evaluation')
+        
+end
 end
\ No newline at end of file