Jump length distros seem to work for boundaries at 3, 5, 7.

prob_laplace.m

@@ -219,7 +219,7 @@ end
 toc
 %% Normalization factor
 % N = normalization(T, x0, 0, t_ind, direc, -params{1})
-N=sum(p)/length(p)*l_max
+N=sum(p)/length(p)*l_max;
 % N = sum(p)*0.001;
 m = sum(ls.*p/N)/length(p)*l_max;
 % should sum to one
@@ -234,13 +234,15 @@ 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
+%% Now for d = 0.5 across the entire domain for different boundary position
 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];
+t = [0, 0.05, 0.1];
 direc = 1;
-x0 = sort(5-direc*(0:0.002:3.02));
-%% Run simulations for 'delta' IC across outside
+bp = [3, 5, 7];
+for j = 1:length(bp)
+x0 = sort(bp(j)-direc*(0:0.001:3.02));
+% Run simulations for 'delta' IC across outside
 T = {};
 parfor i = 1:length(x0)
     tic
@@ -248,18 +250,20 @@ parfor i = 1:length(x0)
     T{i}.solve_tern_frap();
     toc
 end
-%% Calculate probabilities for each jump length in ls.
+% Calculate probabilities for each jump length in ls.
 l_max = 3;
-t_ind = 4;
-ls = -direc*(0.000:0.002:l_max);
+t_ind = 3;
+ls = -direc*(0.000:0.005: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);
+    p(i) = int_prob(ls(i), T, x0, direc, t_ind, bp(j), 0);
     toc
 end
 toc
+save(['prob_7_7_bp', num2str(bp(j))]);
+end
 %% 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);