
% Ploting lines on top of trajectories of two solitons generated by NLS
%
% This is the IC:
 L=100;
 A = 0.25;
 x01 = -L+5/A;
 phi01=1;
 c1 = 1.1;
% u = A*sech(A*(x-x0)).*exp(i*(c*x+phi0));
 x02 = -L+12/A;
 phi02=0;
 c2 = 0.75;
% u1 = A*sech(A*(x-x0)).*exp(i*(c*x+phi0));
% u=u+u1;

t0=allt(1);
t1=allt(length(allt));
x0=allx(1);
x1=allx(length(allx));

M = max(max(abs(allu)));

figure(5)
hold on
plot3([t0 t1],[x01 x01+c1*x1],[1 1]*M,'k--','LineWidth',2);
plot3([t0 t1],[x02 x02+c2*x1],[1 1]*M,'k--','LineWidth',2);
%set(gca,'FontSize',16)

