  function P3()
% Checks out BestEccentric
clc
disp(' k       xc     phi(degrees)      maxErr(Days-Hours-Minutes)')
disp('------------------------------------------------------------------')
for k=1:20
    [xc,phi,maxErr] = BestEccentric(k);
    phi = 180*phi/pi;
    Days = floor(maxErr);
    maxErr = (maxErr-Days)*24;
    Hours = floor(maxErr);
    maxErr = (maxErr-Hours)*60;
    Minutes = floor(maxErr);
    fprintf('%3d  %10.6f   %10.6f       %5d %2d %2d\n',k,xc,phi,Days,Hours,Minutes)
end