ClearAll;
$DefaultFont = {"Times-Roman",10};
Needs["Statistics`DataManipulation`"];
Needs["Graphics`Graphics`"];
bods = Table[16.5 + k,{k,1,10}];
heads = {2.77,2.97,3.24,3.52,3.76,3.99,4.33,4.67,4.78,5.04};
datapoints = Map[Thread,Transpose[{bods,heads}]];
logdatapoints = Map[Thread,Transpose[{Log[bods],Log[heads]}]];
G[t_]:= Fit[Drop[logdatapoints,-2],{1,t},t];
c0 = Coefficient[G[t],t,0]; c1 = Coefficient[G[t],t,1];
W[t_]:= Exp[c0] t^c1;
plot = Show[{ListPlot[logdatapoints,PlotRange->{{2.85,3.35},{1,1.65}},
	AxesOrigin->{2.8,0.95},PlotLabel->
	"Log of blue whale head length versus log of body length"],
	Plot[G[t],{t,2.8,3.3},PlotStyle->Dashing[{0.001,0.01}]]}];
Show[GraphicsArray[{{plot}}],PlotLabel->"Figure 22.4"]