diff --git a/examples/odeintw_demo_complex.py b/examples/odeintw_demo_complex.py
index 7e17736..3829c55 100644
--- a/examples/odeintw_demo_complex.py
+++ b/examples/odeintw_demo_complex.py
@@ -33,7 +33,7 @@ def bjac_cols(y, t, c):
 
 
 def bjac_rows(y, t, c):
-    return np.row_stack((np.r_[0, np.diag(c, 1)], np.diag(c)))
+    return np.vstack((np.r_[0, np.diag(c, 1)], np.diag(c)))
 
 
 c = np.array([[-20+1j, 5-1j,      0,       0],
diff --git a/odeintw/__init__.py b/odeintw/__init__.py
index bfd47cc..22035a9 100644
--- a/odeintw/__init__.py
+++ b/odeintw/__init__.py
@@ -6,4 +6,4 @@
 from ._odeintw import odeintw
 
 
-__version__ = "1.0.1.dev0"
+__version__ = "1.0.1.dev1"
diff --git a/odeintw/tests/test_odeintw.py b/odeintw/tests/test_odeintw.py
index c89070b..0e4c8f3 100644
--- a/odeintw/tests/test_odeintw.py
+++ b/odeintw/tests/test_odeintw.py
@@ -108,7 +108,7 @@ def system3_bjac_cols(y, t, c):
 
 
 def system3_bjac_rows(y, t, c):
-    return np.row_stack((np.r_[0, np.diag(c, 1)], np.diag(c)))
+    return np.vstack((np.r_[0, np.diag(c, 1)], np.diag(c)))
 
 
 def test_system3():