prob_laplace.m before sending to Stefano for integral comparison.

prob_laplace.m

@@ -99,11 +99,7 @@ e = @(chi) sqrt(3/8*(chi-2));
 t = [0, 0.05, 0.1];
 T1 = Ternary_model(0, 'Gauss', [-5, b(chi, nu), 0.5, e(chi), ...
                      0.16, 0.2, 10, a], t, 2);
-% x0 = T1.x(find(T1.x>a, 1):end);
-% x0 = [x0(1:20:8000), x0(8001:2:end)];
-x0 = x0(x0<7.1);
-% x0 = x0(1:10:end);
-x0 = 5.0:0.001:6.01;
+x0 = 5.0:0.002:7.01;
 %%
 T = {};
 tic
@@ -115,9 +111,9 @@ T{i} = Ternary_model(0, 'Gauss', [-a, b(chi, nu), 0.5, e(chi), ...
 T{i}.solve_tern_frap();
 end
 toc
-save prob_laplace_X_7_7_short
+save prob_laplace_X_7_7_short_2
 %%
-ls = 0.0005:0.005:1;
+ls = 0.0005:0.005:2;
 p = nan(1, length(ls));
 parfor i = 1:length(ls)
 p(i) = int_prob(ls(i), T, x0);
@@ -141,7 +137,7 @@ N = normalization(T, x0);
 toc
 %% Write to file
 csvwrite('jump_length_7_7_3.csv', ls)
-csvwrite('prob_7_7_3.csv', p);
+csvwrite('prob_7_7_3.csv', p/N);
 %% distribution for x0 can be taken from phi_tot (steady state)
 function p = int_prob(l, T, x0)
 delta_x0 = diff(x0);