clear all
%Initial Setup

varc = 0;
G = 30; c0=1;
T=50;
IP=[-18 -13 -8 -3];
for i=1:numel(IP)
    
deltat=0.01; 
[X,Y,K,L,nmax,E1,E2,PE1,PE2,epsilon,l0] = subroutine(T, deltat, -20, IP(i), 2, c0, G, varc, 1.4);

%PLOT PARTICLE PATH
M=E1(1)
figure(1);
plot(X,Y); hold on; plot(0, 0, 'k*'); %hold off
axis equal
title('Figure 1:  Wave-packet paths as modeled by forward-Euler scheme')
xlabel('x position')
ylabel('y position')
text(-15,23,'\Delta \itt \rm= 0.01')
text(-15,21,'T = 50')
text(-15,19,'\itc_0 \rm=1')
text(-15,17,'\Gamma \rm = 30')
text(-15,15,'Wave packets released from') 
text(-15,13,'     \itx \rm = -20 and')
text(-15,11,'     \ity \rm = -18, -13, -8, -3')
PE1
PE2
epsilon
l0
end