clear all
%Initial Setup

varc = 0;
G = 350; c0=1;
T=50;
deltat=0.0009; 
[X,Y,K,L,nmax,E1,E2,PE1,PE2,l0] = subroutine(T, deltat, -20, -5, 2, c0, G, varc, 1.4);

%PLOT PARTICLE PATH
M=E1(1)
figure(4);
plot(X,Y); hold on; plot(0, 0, 'k*'); %hold off
axis equal
title('Figure 4:  Path of a wave-packet that cannot get totally sucked in')
xlabel('x position')
ylabel('y position')

PE1
PE2
epsilon
l0

prminp= (  M + sqrt(M^2-2*w*G*M/pi) )/(2*w)
prminn= (  M - sqrt(M^2-2*w*G*M/pi) )/(2*w)
nrminp= ( -M + sqrt(M^2-2*w*G*M/pi) )/(2*w)
nrminn= ( -M - sqrt(M^2-2*w*G*M/pi) )/(2*w)
figure(4)
text(-30,15,'\itM\rm = -1 < 0')
text(-30,13,'\Omega = 2.1811')
text(-30,11,'Wave-packet released from \bfx\rm = (-20, -5)')
text(-30,9,'  with \bfk\rm = (2, .45)')
text(-30,7,'\Delta\itt\rm = 0.0009')
text(-30,5,'\Gamma = 350')
text(-30,3,'T = 50')
text(-30,1,'\itc_0 = \rm1')