(*^
::[ 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 = "Macintosh Mathematica Notebook Front End Version 2.2";
MacintoshStandardFontEncoding;
fontset = title, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, e8, 24, "Times";
fontset = subtitle, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, e6, 18, "Times";
fontset = subsubtitle, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeTitle, center, M7, italic, e6, 14, "Times";
fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, bold, a20, 18, "Times";
fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, bold, a15, 14, "Times";
fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, a12, 12, "Times";
fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 10, "Times";
fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, L-5, 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, R65535, L-5, 12, "Courier";
fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier";
fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, B65535, L-5, 12, "Courier";
fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, 12, "Courier";
fontset = name, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, 10, "Geneva";
fontset = header, inactive, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = leftheader, inactive, L2, 12, "Times";
fontset = footer, inactive, noKeepOnOnePage, preserveAspect, center, M7, 12, "Times";
fontset = leftfooter, inactive, L2, 12, "Times";
fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 10, "Times";
fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, 12, "Times";
paletteColors = 128; automaticGrouping; currentKernel;
]
:[font = input; preserveAspect]
ClearAll;
$DefaultFont = {"Times-Roman",10};
Needs["Graphics`Graphics`"];
njags = 100;
mean[c_]:= Gamma[1+1/c];
mu2 = mean[2];mu5 = mean[5];
l2 = Graphics[{Dashing[{0.001,0.01}],Line[{{mu2,0},{mu2,f2[mu2]}}]}];
l5 = Graphics[{Dashing[{0.001,0.01}],Line[{{mu5,0},{mu5,f5[mu5]}}]}];
poly[a_,b_,f_]:= Join[Table[{a+i(b-a)/njags,f[a+i(b-a)/njags]},
{i,0,njags}],{{b,0},{a,0}}];
shade[a_,b_,c_,f_]:= Graphics[{GrayLevel[c],Polygon[poly[a,b,f]]}];
f[m_,z_]:= m z^(m-1) E^(-(z)^m);
f2[z_]:= f[2,z];f5[z_]:= f[5,z];
M2 = (Log[2])^(1/2);M5 = (Log[2])^(1/5);
T2[z_]:= (z - mu2)f2[z];
ushading = Show[l2,Axes->True,
PlotRange->{{0,2.75},{0,0.9}},
Ticks->{{0.5,{mu2,FontForm[m,{"Symbol",9}]},1.5},{0,0.25,0.5,0.75}},
PlotLabel->"Probability density",
AspectRatio->0.92/GoldenRatio,
AxesLabel->{"x/s","s f"}];
lshading = Show[shade[0,mu2,0.65,T2],shade[mu2,2.75,0.65,T2],
Graphics[Text[FontForm["+",{"Times-Roman",24}],{1.35,0.1}]],
Graphics[Text[FontForm["-",{"Times-Roman",30}],{0.4,-0.15}]],
Axes->True,
PlotRange->{{0,2.75},{-0.35,0.25}},
Ticks->{{1.0,1.5,2.0,2.5},
{-0.2,0,0.2}},
PlotLabel->"Turning density",
AxesLabel->{"x/s","T"}];
plotu:= Plot[f2[z],{z,0.001,3}];
plotl:= Plot[T2[z],{z,0.001,3}];
upper = Show[ushading,plotu];
lower = Show[lshading,plotl];
Show[GraphicsArray[{{upper},{lower}}],
PlotLabel->" Figure 26.2\n
Mean of Weibull distribution (c = 2)"]
^*)