% 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*(y-x0)).*exp(i*(c*y+phi0));
 x02 = -L+12/A;
 phi02=0;
 c2 = 0.75;
% u1 = A*sech(A*(y-x0)).*exp(i*(c*y+phi0));
% u=u+u1;

x0=allt(1);
x1=allt(length(allt));
y0=ally(1);
y1=ally(length(ally));

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

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