clear all
%Initial Setup

varc = 0;
G = 350; c0=1;
T=50;
%PLOT CONVERGENCE TEST FOR Rmin

dts = [0.1 .05 .025 .015 .01 .0075 .005 .003 .0015];
for i=1:numel(dts);
  dts(i)
  [X,Y,K,L,nmax,E1,E2,PE1,PE2,l0] = subroutine(T, dts(i), -20, -5, 2, 1, G, varc, 1.4);
  R=sqrt(X.^2+Y.^2);
  Rmin(i)=min(R);
  ind=find(R==min(R));
  TC(i)=ind*dts(i)
  M=E1(1)
  w=E2(1)
end
figure(6)
plot(dts,Rmin, 'r+')
title('Figure 6:  Convergence test for minimum radius for path shown in Figure 4')
xlabel('timestep')
ylabel('minimum radius')
figure(7)
plot(dts,TC, 'g+')
title('Figure 7:  Convergence test for time to reach minimum radius for path shown in Figure 4')
xlabel('timestep')
ylabel('time to reach minimum radius')