(*^ ::[ 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 = "NeXT Mathematica Notebook Front End Version 2.2"; NeXTStandardFontEncoding; fontset = title, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e8, 24, "Times"; ; fontset = subtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e6, 18, "Times"; ; fontset = subsubtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, italic, L1, e6, 14, "Times"; ; fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, bold, L1, a20, 18, "Times"; ; fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, bold, L1, a15, 14, "Courier"; ; fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, L1, a12, 12, "Times"; ; fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 10, "Times"; ; fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, L1, 12, "Courier"; ; fontset = output, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier"; ; fontset = message, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, L1, 12, "Courier"; ; fontset = name, inactive, noPageBreakInGroup, nohscroll, preserveAspect, M7, italic, B65535, L1, 10, "Times"; ; fontset = header, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, L1, 12, "Times"; ; fontset = leftheader, 12; fontset = footer, inactive, nohscroll, noKeepOnOnePage, preserveAspect, center, M7, italic, L1, 12, "Times"; ; fontset = leftfooter, 12; fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12, "Courier"; ; fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; paletteColors = 128; automaticGrouping; currentKernel; ] :[font = section; inactive; Cclosed; preserveAspect; startGroup] WagonPackelFig.1.eps :[font = input; preserveAspect] H[n_] := Sum[1./k, {k, 1, n}] :[font = input; preserveAspect] HFast[n_] := N[EulerGamma] + PolyGamma[0., n+1] :[font = input; preserveAspect; endGroup] Show[ListPlot[Table[{n, H[n]}, {n, 1, 100, 2}]], Plot[Log[n], {n, 1, 100}]]; :[font = section; inactive; Cclosed; preserveAspect; startGroup] WagonPackelFig.2.eps :[font = input; preserveAspect; endGroup] ListPlot[Table[{n, H[n] - Log[n]}, {n, 1, 240, 3}], PlotRange->{0, 1}]; :[font = section; inactive; Cclosed; preserveAspect; startGroup] WagonPackelFig.3.eps :[font = input; preserveAspect; endGroup] a[i_] := 1/i; rects = Table[Rectangle[{i-1, 0}, {i, a[i]}], {i, 10}]; lines = Table[Line[ {{i-1, 0}, {i-1, a[i]}, {i, a[i]}, {i,0}}], {i, 10}]; Plot[1/x, {x,.5,10.5}, PlotPoints->50, PlotStyle->AbsoluteThickness[1.5], Prolog->{{GrayLevel[0.5], rects}, {AbsoluteThickness[.5], lines}}, Axes->Automatic, PlotRange->{{0,10.5},{0,1}}, DefaultFont->{"Helvetica", 6}, AxesOrigin->{0, 0}, Ticks->{Range[0,10], {ToExpression@#,#}& /@ {"1/2","1/3","1/4","1"}}]; ;[s] 1:0,0;464,-1; 1:1,7,6,Courier,1,10,0,0,0; :[font = section; inactive; Cclosed; preserveAspect; startGroup] WagonPackelFig.4.eps :[font = input; preserveAspect; fontSize = 10; fontName = "Courier"; endGroup] a[i_] := 1/i; nPieces = (200); step = (20); rects = Table[Rectangle[{i-1, 0}, {i, a[i]}], {i, 2, 10}]; lines = Table[Line[{{i-1, 0}, {i-1, a[i]}, {i, a[i]}, {i, 0}}], {i, 2,10}]; {gray, black, white} = GrayLevel /@ {.5,0,1}; pieces = Table[Polygon[Join[Table[{j/step, a[j/step + i]}, {j, 0, step}], {{1, a[i+1]}, {1, a[i]}}]], {i, nPieces}]; image = Plot[1/x, {x, 1.01, 10.5}, PlotPoints->(50), MaxBend->1, PlotStyle->{black, AbsoluteThickness[1.5]}, Prolog->{Text[FontForm["y = ", {"Times-Italic", 12}], {2.3, 0.7}, {-1,0}], Text[FontForm["1", {"Times", 11}], {3.55, 0.745}, {-1,0}], Text[FontForm["/", {"Times", 12}], {3.73, 0.7}, {-1,0}], Text[FontForm["x", {"Times-Italic", 12}], {3.88, 0.65}, {-1,0}], {white, Rectangle[{0,0}, {1,1}]}, {AbsoluteThickness[1], Line[{{11,0},{0,0}, {0,1}, {1,1}}]}, {pieces}, {white, Polygon[Join[Table[{i,a[i]}, {i,1,10,.1}], {{10,0}, {1,0}}]]}, {AbsoluteThickness[.5], lines}, {gray, rects, black, AbsoluteThickness[.5], lines} }, Axes->Automatic, PlotRange->{{0,11},{0,1.07}}, AxesOrigin->{0,0}, DefaultFont->{"Helvetica", 6}, Ticks->{Range[0,10], {ToExpression@#,#}& /@ {"1/2","1/3","1/4","1"}}, DisplayFunction->Identity]; Show[image, DisplayFunction->$DisplayFunction]; ;[s] 13:0,0;27,1;32,2;42,3;46,4;415,5;419,6;779,7;784,8;904,9;909,10;1034,11;1038,12;1357,-1; 13:1,7,6,Courier,1,10,0,0,0;1,7,6,Courier,1,10,65535,0,0;1,7,6,Courier,1,10,0,0,0;1,7,6,Courier,1,10,65535,0,0;1,7,6,Courier,1,10,0,0,0;1,7,6,Courier,1,10,65535,0,0;1,7,6,Courier,1,10,0,0,0;1,7,6,Courier,1,10,65535,0,0;1,7,6,Courier,1,10,0,0,0;1,7,6,Courier,1,10,65535,0,0;1,7,6,Courier,1,10,0,0,0;1,7,6,Courier,1,10,65535,0,0;1,7,6,Courier,1,10,0,0,0; ^*)