diff --git a/prob_laplace.m b/prob_laplace.m
index 636960acaec166e6c251a1f26d4500ed704b12dc..1d23534442660391c2aef6a4eed68e5b8ef29483 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -234,6 +234,32 @@ cd /Users/hubatsch/Nextcloud/Langevin_vs_MeanField/Data_Figs_FokkPla/jump_length
 csvwrite('jump_length_7_7_lb01.csv', ls)
 csvwrite('prob_7_7_lb01.csv', p);
 
+%% Now for d = 0.5 across the entire domain
+params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 'const', 0.5, 10, 7, 0,...
+          'Constituent', 0};
+t = [0, 0.05, 0.1, 1];
+direc = 1;
+x0 = sort(5-direc*(0:0.002:3.02));
+%% Run simulations for 'delta' IC across outside
+T = {};
+parfor i = 1:length(x0)
+    tic
+    T{i} = Ternary_model(0, 'Gauss', [params(1:7), {x0(i)}, params(9:end)], t, 0.2);
+    T{i}.solve_tern_frap();
+    toc
+end
+%% Calculate probabilities for each jump length in ls.
+l_max = 3;
+t_ind = 4;
+ls = -direc*(0.000:0.002:l_max);
+p = nan(1, length(ls));
+tic
+parfor i = 1:length(ls)
+    tic
+    p(i) = int_prob(ls(i), T, x0, direc, t_ind, 5, 0);
+    toc
+end
+toc
 %% Mean cross jump length chi 7/3
 run_jump_lengths(-5, 1, 7/3, 1, 4, '73_1_01', 10^-6, 3, 0.5);
 run_jump_lengths(-5, 1, 7/3, 1, 4, '73_1_10', 10^-6, 4, 0.5);