Moving boundary works for v=10^-15.

@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(3, 5.7, 500), linspace(5.705, 6.1, 6000),...
-       linspace(6.105, 8, 2250), linspace(8.05, 12, 750)];
+T.x = [linspace(3, 5.7, 500), linspace(5.705, 6.05, 16000),...
+       linspace(6.06, 7.1, 1500), linspace(7.15, 10, 750)];
 end

@Ternary_model/plot_sim.m

@@ -10,8 +10,8 @@ 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([-1, 1]);
-        xlim([0, 10]);
+        cla; xlim([-T.a-1, -T.a+3]); ylim([0, 1.2]);
+        xlim([3, 10]);
         ax = gca;
         ax.FontSize = 12;
         xlabel('position'); ylabel('probability [not normalized]');
@@ -20,7 +20,7 @@ if strcmp(mode, 'mov')
         plot(T.x, T.sol(i, :), 'LineWidth', 2);
         text(abs(-T.a+2), 0.7, num2str(T.t(i)));
         shg
-        pause();
+%         pause();
 %         print([num2str(i),'.png'],'-dpng')
     end
 elseif strcmp(mode, 'plot')