diff --git a/prob_laplace.m b/prob_laplace.m
index 4ae50c08dcd660e35559328b87043354ee9f2376..2e94c11115ddde092877096e31708b75938fb264 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -189,12 +189,10 @@ end
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % t = linspace(0, 10, 9000); % Entropy
 t = linspace(0, 10, 21); % Jump size distro
-tic
 params_mov_bound = {-10, b(7/3, 10^-6), 0.5, e(7/3),...
                     0, 1, 20, 7, 0.5, 'Client', 0};
 T_mov = Ternary_model(0, 'phi_tot', params_mov_bound, t, 0);
 T_mov.solve_tern_frap();
-toc
 norm_fac = 1/sum(diff(T_mov.x).*...
                  (T_mov.sol(1, 1:end-1)+T_mov.sol(1, 2:end))/2);
 s = '~/Nextcloud/Langevin_vs_MeanField/Data_Figs_FokkPla/';
@@ -225,6 +223,7 @@ end
 % csvwrite([s, 'Mov_Bou_Entr.csv'], [T_mov.t; s_dot])
 %% MOVING BOUNDARY JUMP LENGTH DISTRIBUTION
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+load entropy_mov_bound.mat
 direc = 1; % 1: jump from left to right, -1: jump from right to left.
 x0 = sort(10-direc*(0:0.02:4.01));
 t_snap = repmat((0.25:0.5:9.75)', 1, length(x0));
@@ -235,7 +234,7 @@ F = {};
 t = linspace(0, 10, 21);
 v = T_mov.v;
 params = params_mov_bound;
-params{end-1} = 0;
+params{end-2} = 0;
 for i = 1:si_t(2)
     tic
     parfor j = 1:si_t(1)