function GradeIt(t1,t2,wantSimulation)

% For grading I looked at
% GradeIt(700,1000,0) No retrograde motion CVLerr = .005 CVLtime = .397
% GradeIt(2000,3000,0) One retrograde period CVLerr = .086 CVLtime = .771
% GradeIt(0,5000,0) Many retrograde periods CVLerr = .186 CVLtime = 2.268

close all
% Obtain and display the optimal parameters...
tic
%[P1,theta1,P2,r2,theta2,P3,r3,theta3,err] = ThreeCircleCVL(t1,t2);
[P1,theta1,P2,r2,theta2,P3,r3,theta3,err] = ThreeCircle(t1,t2);
t = toc;
clc
fprintf('t1 = %4d\n',t1)
fprintf('t2 = %4d\n\n',t2)
fprintf('err = %5.3f\n',err)
fprintf('time = %5.3f\n\n',t)

fprintf(' P1: %12.6f\n',P1)
fprintf('theta1: %12.6f\n',theta1)
fprintf(' P2: %12.6f\n',P2)
fprintf(' r2: %12.6f\n',r2)
fprintf('theta2: %12.6f\n',theta2)
fprintf(' P3: %12.6f\n',P3)
fprintf(' r3: %12.6f\n',r3)
fprintf('theta3: %12.6f\n',theta3)

ErrorPlot(t1,t2,P1,theta1,P2,r2,theta2,P3,r3,theta3)
set(gcf,'position',[700 300 800 600])
if wantSimulation
pause(3)
close all
Demo(t1,t2,P1,theta1,P2,r2,theta2,P3,r3,theta3)
end