Integration to obtain jump size works. Need normalization.

prob_laplace.m

@@ -93,14 +93,16 @@ end
 %% Solve integrals
 a = 5;
 nu = 10^-6;
-chi = 7/3;
+chi = 7.7/3;
 b = @(chi, nu) nu^(1/3)*sqrt(chi/(chi-2));
 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, 3);
 x0 = T1.x(find(T1.x>a, 1):end);
-x0 = x0(1:10:end);
+x0 = [x0(1:20:8000), x0(8001:2:end)];
+x0 = x0(x0<7.1);
+% x0 = x0(1:10:end);
 % x0 = 0.0:0.001:0.1;
 %%
 T = {};
@@ -113,6 +115,7 @@ 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
 %%
 ls = 0.0005:0.005:1;
 p = nan(1, length(ls));
@@ -122,7 +125,7 @@ end
 %% Do we need to look at the left side as well?
 figure; hold on;
 plot(ls, p);
-% plot(ls, q);
+% plot(ls, q/50000);
 %%
 int_prob(0.01, T, x0+a)
 %%