ClearAll;
$DefaultFont = {"Times-Roman",10};
Needs["Graphics`Graphics`"];
f[m_,b_]:= 36 (m^(5/3) b^(-3/2))/1000;
shade[m_,b_]:= ((f[1,0.1] - f[m,b])/(f[1,0.1] - f[0.1,1]))^10;
xticks = {{0.9," "},{1,1}};
yticks = Table[{0.1 i," "},{i,1,10}];
zticks = {0.5,1};
axeslabels = Graphics3D[{Text[FontForm["x",{"Times-Roman",12}],{1.16,0,0.1}],
					   Text[FontForm["y",{"Times-Roman",12}],{0,1.16,0.1}],
					   Text[FontForm["z",{"Times-Roman",12}],{0.1,0,1.25}]}];
plot = Plot3D[{f[m,b],GrayLevel[shade[m,b]]},{m,0.1,1},{b,0.1,1},
				AxesEdge->{{-1,-1},{-1,-1},{-1,-1}},
				ViewPoint->{2,2,2},
				BoxRatios->{1,1,1},
				Shading->True,
				ClipFill->None,
				PlotRange->{{0,1.16},{0,1.16},{0,1.26}},
				PlotLabel->"Figure 24.1",
				Boxed->False,
				Axes->True,
				Mesh->True,
				MeshStyle->Dashing[{0.01,0.01}],
				Ticks->{xticks,yticks,zticks}];
Show[plot,axeslabels]