(*^ ::[ Information = "This is a Mathematica Notebook file. It contains ASCII text, and can be transferred by email, ftp, or other text-file transfer utility. It should be read or edited using a copy of Mathematica or MathReader. If you received this as email, use your mail application or copy/paste to save everything from the line containing (*^ down to the line containing ^*) into a plain text file. On some systems you may have to give the file a name ending with ".ma" to allow Mathematica to recognize it as a Notebook. The line below identifies what version of Mathematica created this file, but it can be opened using any other version as well."; FrontEndVersion = "X Window System Mathematica Notebook Front End Version 2.2"; X11StandardFontEncoding; fontset = title, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, e8, 24, fontName, "times"; fontset = subtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, e6, 18, fontName, "times"; fontset = subsubtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, italic, e6, 14, fontName, "times"; fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, bold, a20, 18, fontName, "times"; fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, bold, a15, 14, fontName, "times"; fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, a12, 12, fontName, "times"; fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 10, fontName, "times"; fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, 12, fontName, "courier"; fontset = output, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, fontName, "courier"; fontset = message, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, 12, fontName, "courier"; fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, 12, fontName, "courier"; fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, 12, fontName, "courier"; fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, 12, fontName, "courier"; fontset = name, inactive, noPageBreakInGroup, nohscroll, preserveAspect, M7, italic, 10, fontName, "times"; fontset = header, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, 12, fontName, "times"; fontset = leftheader, 12, fontName, "times"; fontset = footer, inactive, nohscroll, noKeepOnOnePage, preserveAspect, center, M7, italic, 12, fontName, "times"; fontset = leftfooter, 12, fontName, "times"; fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "courier"; fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times"; fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, fontName, "times";paletteColors = 128; showRuler; automaticGrouping; currentKernel; ] :[font = title; inactive; preserveAspect; startGroup] Sequences :[font = text; inactive; preserveAspect] To what does the sequence a[n] = 1 + (-1)^n / n^2 converga? The answer leaps out of a graph if we plot the points ( 1/n, a[n] ). Please note that the index n increases as we move to the LEFT!!!! (If that really bothers you, plot ( - 1/n, a[n] ) instead.) :[font = input; preserveAspect] a[n_] = 1.+ (-1)^n / (n^2); spoints = Table[ {1/n,a[n]}, {n,1,128} ]; plot1 = ListPlot[ spoints ,PlotStyle->PointSize[0.02] ,PlotRange->{{-.01,1.1},{-.1,2}} ,Ticks->{{ 0.5}, {0,1,2} } ] ; :[font = text; inactive; preserveAspect] Looks like the points cluster near ( 0, 1). The limit of the sequence is the y-coordinate of the limit point, which is 1. To illustrate the epsilon-tolerances, we zoom in on the graph, and add two lines for 1 + epsilon and 1 - epsilon. In this case, epsilon is .004. :[font = input; preserveAspect] a[n_] = 1.+ (-1)^n / (n^2); spoints = Table[ {1/n,a[n]}, {n,1,512} ]; plot1 = ListPlot[ spoints ,PlotStyle->PointSize[0.02] ,PlotRange->{{-.01,.15},{.98,1.02}} ,Ticks->{{ 0.1}, Automatic} ] ; plot2 =Plot[ {1.004,.996},{x,0,.15} ]; Show[ plot1, plot2 ]; :[font = text; inactive; preserveAspect] Another more "advanced" concept is the idea of a limit superior. Some sequences don't settle down to a limit but instead careen about wildly. The following graph illustrates a sequence that has a limit superior of 2.5. That is, no matter how late in the sequence we may be, we can always find a value that is very close to 2.5, but we cannot do it for any larger value, for example, 2.6. :[font = input; preserveAspect] c[n_] = 3/2 + n / (n +50) Sin[3n]; spoints = Table[ { 1/n, c[n] }, {n,1,4096} ]; plot1 = ListPlot[ spoints ,PlotStyle->PointSize[0.01] ,AspectRatio->1/GoldenRatio ,PlotRange->{{-.0001,.02},{0,3}} ,Ticks->{{ }, {.5,1,1.5,2,2.5}} ] ; :[font = text; inactive; preserveAspect] Can you guess the technical meaning of a limit inferior? :[font = text; inactive; preserveAspect] The next sequence may be more challenging. What is the liim sup and lim inf of this sequence? :[font = input; preserveAspect] c[n_] = 2/5 + (n+50) / (n) Sin[3n]; spoints = Table[ { 1/n, c[n] }, {n,1,4096} ]; plot1 = ListPlot[ spoints ,PlotStyle->PointSize[0.01] ,AspectRatio->1/GoldenRatio ,PlotRange->{{-.0001,.02},{-1,2}} ,Ticks->{{ .01}, Automatic } ] ; :[font = text; inactive; preserveAspect] The following sequence consists of three different parts, each getting a turn at determining the sequences value. By using the same graphing trick, we can read off the set of subsequential limits. Like the lim sup and lim inf, a subsequential limit is a limiting value that we can obtain from some convergent subsequence. The lim inf happens to be the least of the set of subsequential limits. What do you think the biggest subsequential limit is called? :[font = input; preserveAspect] b1[n_] = 400 (-1)^(n) /n^2; b2[n_] = .4+ ((n^1.3)/100+2)/(5 (n^1.3)/100 +1) Sin[3 n]; b3[n_] = 1 - 20/n; spoints = Flatten[ Table[ { {1/(3n), b1[n]}, {1/(3n +1), b2[n]}, {1/(3n +2), b3[n]} }, {n,16,2048}], 1 ]; ListPlot[ spoints ,PlotStyle->PointSize[0.01] ,AspectRatio->1 ,PlotRange->{{-.001,.019},{-.4,1.1}} ,AxesOrigin->{0,0} ,Ticks->{{0.005, 0.01, 0.015}, Automatic} ] ; :[font = text; inactive; preserveAspect; endGroup] What is the lim sup and lim inf of this sequence? Describe carefully the set of subseqential limits. Use the diagram to explain how you would choose a sequence that converges to a point in the set of subseqential limits. ^*)