diff --git a/scripts/run_experiments.jl b/scripts/run_experiments.jl
index e78258cadd0974d27ec19a9c408415913027da13..1b31e0b4284bde0971b3c9421049d0decf8f564a 100644
--- a/scripts/run_experiments.jl
+++ b/scripts/run_experiments.jl
@@ -210,7 +210,9 @@ function experiment_convergence_test(ctx)
 
     #noise_σ = 0.1
 
-    #f = clamp.(fo .+ noise_σ * randn(Float64, size(fo)), 0., 1.)
+    #seed = 42
+    #rng = MersenneTwister(seed)
+    #f = clamp.(fo .+ noise_σ * randn(rng, Float64, size(fo)), 0., 1.)
     #saveimg(joinpath(ctx.outdir, "input_noisy.png"), f)
 
     #prob = DenoiseProblem(f; λ, β)
@@ -386,7 +388,9 @@ function experiment_convergence_denoising(ctx)
 
     noise_σ = 0.1
 
-    f = clamp.(fo .+ noise_σ * randn(Float64, size(fo)), 0., 1.)
+    seed = 42
+    rng = MersenneTwister(seed)
+    f = clamp.(fo .+ noise_σ * randn(rng, Float64, size(fo)), 0., 1.)
     saveimg(joinpath(ctx.outdir, "input_noisy.png"), f)
 
     prob = DenoiseProblem(f; λ, β)
@@ -411,7 +415,10 @@ function experiment_convergence_inpainting(ctx)
     M = (2,2)
     overlap = (5,5)
 
-    f_mask = rand(Float64, size(fo)) .> 0.5
+
+    seed = 42
+    rng = MersenneTwister(seed)
+    f_mask = rand(rng, Float64, size(fo)) .> 0.5
 
     prob = InpaintProblem(fo, f_mask; λ, β)
     saveimg(joinpath(ctx.outdir, "input_mask.png"), f_mask)