WagonPackelFig.1.eps

H[n_] := Sum[1./k, {k, 1, n}]

HFast[n_] := N[EulerGamma] + PolyGamma[0., n+1]

Show[ListPlot[Table[{n, H[n]}, {n, 1, 100, 2}]],
     Plot[Log[n], {n, 1, 100}]];

WagonPackelFig.2.eps

ListPlot[Table[{n, H[n] - Log[n]}, {n, 1, 240, 3}],
     PlotRange->{0, 1}];

WagonPackelFig.3.eps

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"}}];

WagonPackelFig.4.eps

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];