diff --git a/tutorials/adaptivepd.py b/tutorials/adaptivepd.py
index 0fceb45..fa29d79 100644
--- a/tutorials/adaptivepd.py
+++ b/tutorials/adaptivepd.py
@@ -21,6 +21,7 @@
 import pyproximal
 
 plt.close("all")
+np.random.seed(10)
 
 
 def callback(x, f, g, K, cost, xtrue, err):
diff --git a/tutorials/constrained.py b/tutorials/constrained.py
index 4359974..3c1af50 100644
--- a/tutorials/constrained.py
+++ b/tutorials/constrained.py
@@ -27,6 +27,7 @@
 import pyproximal
 
 plt.close("all")
+np.random.seed(10)
 
 ###############################################################################
 # Let's start defining the terms of the quadratic functional
diff --git a/tutorials/denoising.py b/tutorials/denoising.py
index 46c4c95..b845625 100644
--- a/tutorials/denoising.py
+++ b/tutorials/denoising.py
@@ -29,6 +29,7 @@
 import pyproximal
 
 plt.close("all")
+np.random.seed(10)
 
 ###############################################################################
 # Let's start by loading a sample image and adding some noise
diff --git a/tutorials/lowrankcompletion.py b/tutorials/lowrankcompletion.py
index 841a877..a49e914 100644
--- a/tutorials/lowrankcompletion.py
+++ b/tutorials/lowrankcompletion.py
@@ -44,8 +44,8 @@
 
 import pyproximal
 
-np.random.seed(0)
 plt.close("all")
+np.random.seed(0)
 
 ###############################################################################
 # Let's start by loading a sample image
diff --git a/tutorials/nonlinearconstrained.py b/tutorials/nonlinearconstrained.py
index ec7814e..afee690 100644
--- a/tutorials/nonlinearconstrained.py
+++ b/tutorials/nonlinearconstrained.py
@@ -51,6 +51,7 @@
 import pyproximal
 
 plt.close("all")
+np.random.seed(10)
 
 ###############################################################################
 # Let's start by defining the class for the nonlinear functional