Moving boundary closer to working, still oscillating around boundary.

@Ternary_model/create_mesh.m

@@ -29,6 +29,6 @@ function create_mesh(T)
 % %                                    ceil((T.x(ind+20)-T.x(ind-20))*200)),...
 % %            T.x(ind+21:end)];
 % % end
-T.x = [linspace(0, 5.5, 1000), linspace(5.505, 6.2, 5000),...
-       linspace(6.2005, 8, 3000), linspace(8.005, 12, 2000)];
+T.x = [linspace(3, 5.7, 500), linspace(5.705, 6.1, 6000),...
+       linspace(6.105, 8, 2250), linspace(8.05, 12, 750)];
 end

@Ternary_model/plot_sim.m

@@ -10,7 +10,7 @@ end
 if strcmp(mode, 'mov')
     figure(20); hold on;
     for i = 1:nth:length(T.t)
-        cla; xlim([-T.a-1, -T.a+3]); ylim([0, 1]);
+        cla; xlim([-T.a-1, -T.a+3]); ylim([-1, 1]);
         xlim([0, 10]);
         ax = gca;
         ax.FontSize = 12;

prob_laplace.m

@@ -42,7 +42,7 @@ p_in = @(D_p, D_m, ga, x0, x, t) 1./(sqrt(pi*t)*(sqrt(D_m)+...
 x_left = linspace(-4, 0, 1000);
 x_right = linspace(0, 4, 1000);
 %% Moving boundary
-t = linspace(0, 1, 100);
+t = linspace(0, 1, 300);
 tic
 T_mov = Ternary_model(0, 'Gauss', [-6, b(7/3, 10^-15), 0.5, e(7/3),...
                    0, 1, 300, 7, 0.2], t, 0.02);