,
,
, and
(standing in for
), and compare the graph of the model function to the daylight data.
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:=???; B:=???; C:=???; delta:=???; f:=t->A+B*sin(C*(t-delta)); fPlot:=plot(f(t),t=0..12): display(fPlot,DaylightPlot); |
| > |
| > |