  function ShowDailyAngle(t1,t2)
% Graphically displays the daily angle from Day t1 to Day t2..
close all
set(gcf,'position',[100 100 1200 800])
% Compute and plot the daily angles
psi = DailyAngle(t1,t2);
plot(((t1+1):t2)',psi)
step = min(200,floor((t2-t1)/25));
set(gca,'xTick',t1:step:t2,'yTick',-.01:.001:.015)
grid on
shg