nansum for probs in prob_laplace.

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=nansum(p)/length(p)*l_max;
 % N = sum(p)*0.001;
 m = sum(ls.*p/N)/length(p)*l_max;
 % should sum to one
@@ -263,6 +263,9 @@ parfor i = 1:length(ls)
 end
 toc
 save(['prob_7_7_bp', num2str(bp(j))]);
+pa = '/Users/hubatsch/Nextcloud/Langevin_vs_MeanField/data/constD/';
+csvwrite([pa, 'jump_length_7_7_bp', num2str(bp(j)), '.csv'], ls)
+csvwrite([pa, 'prob_7_7_bp', num2str(bp(j)), '.csv'], p);
 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);