bugfix: EE4 was not returning solution

chichilalescu · Dec 8, 2014 · 3619410 · 3619410
1 parent e827dc0
commit 3619410
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/pyNT/ode.py b/pyNT/ode.py
@@ -76,6 +76,7 @@ def EE4(self, h, nsteps, X0):
             X4 += self.rhs(X4)*h/4
             X4 += self.rhs(X4)*h/4
             X[t] = 32*X4/3. - 13.5*X3 + 4*X2 - X1/6
+        return X
     def EE4_alt(self, h, nsteps, X0):
         X = np.zeros((nsteps+1,) + X0.shape,
                      X0.dtype)