Normalization works.

prob_laplace.m

@@ -124,7 +124,7 @@ p(i) = int_prob(ls(i), T, x0);
 end
 %% Do we need to look at the left side as well?
 figure; hold on;
-plot(ls, N*p);
+plot(ls, p/N);
 % plot(lx, xt);
 % plot(ls, q/50000);
 %%
@@ -183,15 +183,15 @@ end
 
 function p = normalization(T, x0)
 delta_x0 = diff(x0);
-p = 0;
-parfor i = 1:length(delta_x0)
+logi = T{1}.x<5;
+delta_x = diff(T{1}.x(logi));
+p_i = zeros(1, length(delta_x0));
+for i = 1:length(delta_x0)
+    sol = T{i}.sol(3, logi);
+    sol = (sol(1:end-1)+sol(2:end))/2;
     x = (x0(i)+x0(i+1))/2;
     p_x0i = T{1}.phi_tot(x, T{1}.a, T{1}.b, T{1}.e, T{1}.u0);
-    for j = 1:length(delta_x0)
-        xj = (x0(j)+x0(j+1))/2;
-        p_j = @(j) interp1(T{i}.x, T{i}.sol(2, :), xj);
-        pj = (p_j(i)+p_j(i+1))/2;
-        p = p + delta_x0(i)*delta_x0(j)*p_x0i*pj;
-    end
+    p_i(i) = delta_x0(i)*sum(sol.*delta_x)*p_x0i;
 end
+p = sum(p_i);
 end
\ No newline at end of file