>

Data entry

>

Sunrise times are "minutes after 5AM". Sunset times are "minutes after 5PM"

> sunrise:=[190,190,190,190,190,190,190,190,190,190,190,189,189,189,189,
188,188,187,187,186,186,185,185,184,183,182,182,181,180,179,178,
178,177,176,175,174,173,172,170,169,168,167,166,165,164,162,
161,160,158,157,156,154,153,152,150,149,148,146,145,144,
143,142,140,139,137,136,134,133,131,130,128,126,125,123,122,
120,119,117,115,114,112,111,109,107,106,104,103,101,99,98,96,
95,93,91,90,88,87,85,84,82,81,79,78,76,75,73,
72,70,69,67,66,64,63,62,60,59,58,56,55,54,52,
51,50,49,48,46,45,44,43,42,41,40,39,38,37,36,
35,34,33,33,32,31,30,30,29,28,28,27,27,26,26,25,
25,24,24,24,23,23,23,22,22,22,22,22,22,22,22,
22,22,22,22,22,23,23,23,23,24,24,24,25,25,26,
26,26,27,28,28,29,29,30,30,31,32,33,33,34,35,
35,36,37,38,39,39,40,41,42,43,44,44,45,46,47,48,
49,50,51,52,53,53,54,55,56,57,58,59,60,61,62,
63,64,65,66,67,67,68,69,70,71,72,73,74,75,76,77,
78,79,79,80,81,82,83,84,85,86,87,88,89,90,91,
91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,
105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,
120,121,122,123,125,126,127,128,129,130,131,132,133,134,135,136,
138,139,140,141,142,143,144,145,147,148,149,150,151,152,153,
155,156,157,158,159,160,161,162,163,164,166,167,168,169,170,
171,172,173,174,175,175,176,177,178,179,180,180,181,182,183,
183,184,185,185,186,186,187,187,188,188,189,189,189,190,190,190]:

> sunset:=[37,38,38,39,40,41,42,43,44,45,46,47,48,49,50,
51,52,53,55,56,57,58,59,60,61,63,64,65,66,67,69,
70,71,72,73,74,76,77,78,79,80,82,83,84,85,86,
87,88,90,91,92,93,94,95,96,97,98,100,101,102,
102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,
117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,
133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,
148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,
163,164,165,166,167,168,169,170,171,172,173,174,175,176,176,
177,178,179,180,181,182,183,184,184,185,186,187,188,188,189,190,
191,191,192,193,193,194,194,195,196,196,197,197,198,198,198,
199,199,199,200,200,200,200,201,201,201,201,201,201,201,201,
201,201,200,200,200,200,200,199,199,198,198,198,197,197,196,
196,195,194,194,193,192,191,191,190,189,188,187,186,185,184,183,
182,181,180,179,178,177,176,175,173,172,171,170,168,167,166,
165,163,162,160,159,158,156,155,153,152,150,149,147,146,144,143,
141,140,138,136,135,133,132,130,128,127,125,124,122,120,119,
117,115,114,112,111,109,107,106,104,102,101,99,98,96,94,
93,91,90,88,86,85,83,82,80,79,77,76,74,73,71,
70,68,67,65,64,63,61,60,59,57,56,55,53,52,51,50,
49,47,46,45,44,43,42,41,40,39,38,37,37,36,35,
34,33,33,32,31,31,30,30,29,29,28,28,28,27,27,
27,26,26,26,26,26,26,26,26,26,26,26,26,27,27,
27,28,28,28,29,29,30,30,31,31,32,33,33,34,35,36]:

The temperatures are the average daily high and low in Indianapolis in degrees Fahrenheit.

> hightemp:=[36,36,36,36,36,36,36,36,36,36,36,36,36,36,36,
36,36,36,36,36,36,36,36,36,36,36,36,36,36,36,37,
37,37,37,37,37, 37,38,38,38,38, 38,39,39,39,39,
39,40,40,40,40, 41,41,41,41,42, 42,42,43,
43,43,44,44,44, 45,45,46,46,46, 47,47,48,48,48,
49,49,50,50,51, 51,51,52,52,53, 53,54,54,55,55, 56,
56,57,57,58,58, 59,59,60,60,61, 61,62,62,63,63,
63,64,64,65,65, 65,66,66,66,67, 67,67,67,68,68,
68,69,69,69,69, 70,70,70,71,71, 71,72,72,72,72,
73,73,73,74,74, 74,75,75,75,76, 76,77,77,77,78, 78,
78,79,79,79,80, 80,80,81,81,81, 82,82,82,82,83,
83,83,83,83,84, 84,84,84,84,84, 84,85,85,85,85,
85,85,85,85,85, 85,85,85,85,85, 85,85,85,85,85,
85,86,86,86,86, 86,86,86,86,86, 86,86,86,85,85, 85,
85,85,85,85,85, 85,85,85,85,85, 85,85,85,85,84,
84,84,84,84,84, 84,83,83,83,83, 83,83,82,82,82, 82,
81,81,81,81,80, 80,80,80,80,79, 79,79,78,78,78,
78,77,77,77,77, 76,76,76,75,75, 75,75,74,74,74,
73,73,73,72,72, 72,71,71,71,70, 70,69,69,68,68,
67,67,67,66,65, 65,64,64,63,63, 62,62,61,60,60, 59,
59,58,57,57,56, 56,55,54,54,53, 53,52,52,51,50,
50,49,49,48,48, 47,47,47,46,46, 45,45,44,44,43,
43,43,42,42,42, 41,41,40,40,40, 40,39,39,39,38,
38,38,38,38,37, 37,37,37,37,37, 37,36,36,36,36, 36]:

> lowtemp:=[20,20,20,20,20, 20,20,20,20,20, 20,20,19,19,19,
19,19,19,19,19, 19,20,20,20,20, 20,20,20,20,20, 20,
20,20,20,20,21, 21,21,21,21,21, 21,21,22,22,22,
22,22,23,23,23, 23,23,24,24,24, 24,25,25,
25,26,26,26,26, 27,27,27,28,28, 28,29,29,29,30,
30,30,31,31,32, 32,32,33,33,34, 34,34,35,35,36, 36,
36,37,37,38,38, 38,39,39,40,40, 40,41,41,41,42,
42,43,43,43,44, 44,44,44,45,45, 45,46,46,46,47,
47,47,47,48,48, 48,49,49,49,49, 50,50,50,51,51,
51,52,52,52,53, 53,53,54,54,54, 55,55,56,56,56, 57,
57,57,58,58,58, 59,59,59,60,60, 60,61,61,61,61,
62,62,62,62,62, 63,63,63,63,63, 63,64,64,64,64,
64,64,64,64,64, 64,64,64,64,65, 65,65,65,65,65,
65,65,65,65,65, 65,65,65,65,65, 65,65,65,64,64, 64,
64,64,64,64,64, 64,64,64,64,64, 63,63,63,63,63,
63,63,62,62,62, 62,62,61,61,61, 61,60,60,60,60, 59,
59,59,59,58,58, 58,58,57,57,57, 56,56,56,55,55,
55,55,54,54,54, 53,53,53,52,52, 52,51,51,50,50,
50,49,49,49,48, 48,48,47,47,47, 46,46,45,45,45,
44,44,44,43,43, 42,42,42,41,41, 41,40,40,39,39, 39,
38,38,38,37,37, 36,36,36,35,35, 34,34,34,33,33,
33,32,32,31,31, 31,30,30,30,29, 29,29,28,28,27,
27,27,26,26,26, 25,25,25,25,24, 24,24,23,23,23,
23,23,22,22,22, 22,22,21,21,21, 21,21,21,21,20, 20]:

Setting up data for a plot

The 'plots' package loads several special plotting commands.
The display and textplot commands are part of this package.

> with(plots):

The next four commands set up the data so that Maple can plot them.

> rise:=[ seq([n,sunrise[n]],n=1..366)]:

> sset:=[ seq([n,sunset[n]],n=1..366)]:

> hight:=[ seq([n,hightemp[n]],n=1..365)]:

> lowt:=[ seq([n,lowtemp[n]],n=1..365)]:

Sunrise plots

The time that the sun rises has a period of one year. The graph looks like this;

> plot(rise,style=point);

To display several plots at once, assign variable names to the plots.
riseplot is the plot of the times that the sun rises in Terre Haute.
Since the times of sunrise are peiodic, we will try to match a sine curve to the data.
sinerise is the plot of a sine function.

curvelabels are some gory, unimportant Maple details for labelling the curves.
The display command displays the plots on a single set of axes.

> riseplot:=plot(rise,style=point):

> sinerise:=plot( 0+ 1*sin((x-0)/1),x=0..366,color=blue):

> curvelabels:=textplot([[280,10,`sine curve`],[60,160,`sun rise`]],color=gold,align=RIGHT):

> display({riseplot,sinerise,curvelabels});

> sinerise:=plot( 105+ 85*sin((x+95)*(2*Pi/365.25)),x=0..366,color=blue):

> display({riseplot,sinerise});

Other plots

Here is a plot of the sunset times.

> plot(sset,style=point);

A plot of the daily high temperatures in Indianapolis.

> plot(hight,style=point);

A plot of the daily low temperatures in Indianapolis.

> plot(lowt,style=point);

>