diff --git a/@Ternary_model/plot_sim.m b/@Ternary_model/plot_sim.m
index d9de961455e1ab060e3dfc658c9e0651493c9246..ef2d48160b5a6566f9d319c75d90b72eb928ce6b 100644
--- a/@Ternary_model/plot_sim.m
+++ b/@Ternary_model/plot_sim.m
@@ -1,4 +1,4 @@
-function plot_sim(T, mode, nth, color)
+function plot_sim(T, mode, nth, color, start)
 %PLOT_SIM Show movie of simulation
 % mode 'mov' ... show movie
 % mode 'plot' ... show plot
@@ -31,7 +31,7 @@ elseif strcmp(mode, 'plot')
     xlabel('x [\mum]'); ylabel('volume fraction'); 
 %     N = max(max(T.sol([1, 2:nth:si(1)], :)));
     N = 1;
-    plot(T.x, T.sol([1, 1:nth:si(1)], :)/N, 'LineWidth', 1, 'Color', color);
+    plot(T.x, T.sol([1, start:nth:si(1)], :)/N, 'LineWidth', 1, 'Color', color);
 %     plot(T.x, T.sol(nth, :), 'LineWidth', 1, 'Color', color);
     plot(T.x, Ternary_model.phi_tot(T.x, T.a, T.b, T.e, T.u0, 0),...
         'LineWidth', 1, 'LineStyle', '--', 'Color', [0.97, 0.55, 0.03]);