diff --git a/prob_laplace.m b/prob_laplace.m
index 2149d2720d4bb794ef9ccbb698703e6c2c3b8f61..549a8fa019f20e07799489a6410661636aa83f62 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -401,6 +401,7 @@ end
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 params = {-5, b(7/3, 10^-6), 0.5, e(7/3), 0, 1, 10, 7, 0, 'Constituent', 0};
 t = linspace(0, 5, 51);
+bp = 5;
 direc = 1; % 1: jump from left to right, -1: jump from right to left.
 x0 = sort(5-abs(direc)/direc*(0:0.002:4.01));
 %% Run simulations for 'delta' IC across outside