function P3Solve()
% Checks out Transits.
clc
T = 5000;
% Look at all possible transits in [0,T]...
Name = {'Mercury', 'Venus', 'Earth', 'Mars'};
for j=2:4
for i=1:j-1
fprintf('Transits of %s as seen from %s\n\n',Name{i},Name{j})
[StartTimes,Duration] = Transits(i,j,T);
fprintf(' Start Duration \n')
fprintf('-------------------------------------------\n')
for k=1:length(StartTimes)
fprintf('%12.6f %12.6f\n',StartTimes(k),24*60*Duration(k))
end
fprintf('\n\n')
end
end