Maple Tool 2

This tool is designed for use with Activity 2 in Section 5.4.

>	with(stats[statplots]):with(plots):

We plot the sunrise and sunset data first for reference.

>

Months:=[1,2,3,4,5,6,7,8,9,10,11,12]; Rise:=[446.3,421.8,379.5,332.4, 296.2,283.2, 296.7, 323.2, 351.3, 380.9,414.0, 442.2]; RisePlot:=scatterplot(Months,Rise,color=navy, thickness=2, view=[0..13,0..1200], xtickmarks=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12], ytickmarks=[0,200,400,600,800,1000, 1200]):

>	Set:=[1031.9,1065.7, 1097.6, 1127.7, 1157.2, 1177.6, 1174.0, 1143.4, 1097.4, 1049.9, 1015.5, 1008.7]; SetPlot:=scatterplot(Months,Set,color=red, thickness=2, view=[0..13,0..1200], xtickmarks=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12], ytickmarks=[0,200,400,600,800,1000, 1200]):

>	display(RisePlot, SetPlot);

We calculate the monthly averages of time of daylight and plot that data.

>	Daylight:=Set-Rise;

>	DaylightPlot:=scatterplot(Months,Daylight,color=navy, thickness=2, view=[0..13,400..1000], xtickmarks=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12], ytickmarks=[400,600,800,1000]):%;

Now enter values for , , , and  (standing in for ), and compare the graph of the model function to the daylight data.

>	A:=(894.4+566.5)/2; B:=(894.4-566.5)/2; C:=Pi/6; delta:=3.1; f:=t->A+B*sin(C*(t-delta)); fPlot:=plot(f(t),t=0..12): display(fPlot,DaylightPlot);

>

>

>