diff --git a/prob_laplace.m b/prob_laplace.m
index 0ae9bb89a328e71493ceaeff0b0724305616be36..6ae42f3f5d9c4b03b3dc9b64392007df9cd4c936 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -304,22 +304,22 @@ end
 %% MOVING BOUNDARY JUMP LENGTH DISTRIBUTION
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % load entropy_mov_bound_in_to_out.mat
-direc = -1; % 1: jump from left to right, -1: jump from right to left.
-x0 = sort(10-direc*(0:0.02:4.01));
+direc = 1; % 1: jump from left to right, -1: jump from right to left.
+x0 = sort(10-direc*(0:0.005:4.01));
 % t_snap = repmat((0.25:0.5:9.75)', 1, length(x0)); % midpoint
 t_snap = repmat((0.:0.5:9.5)', 1, length(x0)); % starting point
 si_t = size(t_snap);
-x0 = sort(10-T_mov.v*t_snap-direc*repmat(0:0.02:4.01, [si_t(1), 1]), 2);
+x0 = sort(10-T_mov.v*t_snap-direc*repmat(0:0.005:4.01, [si_t(1), 1]), 2);
 %% Run simulations for 'delta' IC across outside
 F = {};
 t = linspace(0, 1, 21);
-params = params_mov_bound;
 for j = 1:si_t(1)
     tic
     parfor i = 1:si_t(2)
+        params = params_mov_bound;
+        params{8} = x0(j, i);
+        params{1} = -10+params{9}*t_snap(j);
         F{i, j} = Ternary_model(0, 'Gauss', params, t, 0.1);
-        F{i, j}.x0 = x0(j, i);
-        F{i, j}.a = -10+F{i, j}.v*t_snap(j);
         F{i, j}.solve_tern_frap();
     end
     toc
@@ -336,7 +336,7 @@ for j = 1:n_T
     tic
     F1 = F(:, j);
     bp = bound(j);
-    for i = 1:length(ls)
+    parfor i = 1:length(ls)
         p(i, j) = int_prob(ls(i), F1, x0(j, :), direc, j, bp, 10, T_mov);
     end
     toc