diff --git a/prob_laplace.m b/prob_laplace.m
index 737bb6f3a7c7b4b27ca3fe19560a4d9161afa8f7..d38125cf2b1ead35659f7e80b2abec97c6f831f4 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -193,7 +193,7 @@ params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 0, 1, 10, 7, 0,...
           'Constituent', 0};
 t = [0, 0.05, 0.1, 1];
 direc = 1;
-x0 = sort(5-direc*(0:0.001:3.01));
+x0 = sort(5-direc*(0:0.001:4.01));
 %% Run simulations for 'delta' IC across outside
 T = {};
 parfor i = 1:length(x0)
@@ -205,19 +205,23 @@ parfor i = 1:length(x0)
 end
 % save prob_laplace_X_7_7_short_2
 %% Calculate probabilities for each jump length in ls.
-ls = -direc*(0.000:0.001:3);
+l_max = 4;
+t_ind = 4;
+ls = -direc*(0.000:0.001:l_max);
 p = nan(1, length(ls));
 tic
 parfor i = 1:length(ls)
     tic
-    p(i) = int_prob(ls(i), T, x0, direc, 3, 5, 0);
+    p(i) = int_prob(ls(i), T, x0, direc, t_ind, 5, 0);
     toc
 end
 toc
 %% Normalization factor
-N = normalization(T, x0, 0, 3, direc, 5)
+N = normalization(T, x0, 0, t_ind, direc, -params{1})
 % N = sum(p)*0.001;
-% m = sum(ls.*p/N)/3001*3;
+m = sum(ls.*p/N)/length(p)*l_max;
+% should sum to one
+sum(p/N)/length(p)*l_max
 %%
 sum((T{1}.sol(1, 1:end-1)+T{1}.sol(1, 2:end))/2.*diff(T{1}.x))
 %% Do we need to look at the left side as well?