(*^
::[ 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["Statistics`DataManipulation`"];
Needs["Graphics`Graphics`"];
wheat = {{0,0}, {5,6}, {10,11}, {15,15}, {25,18},
{30,16}, {35,14}, {40,10}, {45,6}, {51,0}};
maize = {{12,0}, {15,9}, {20,20}, {25,29},
{30,35}, {35,38}, {45,30}, {50,7}, {51,0}};
w0[t_]:= t(t+17)(51-t)(5025 - 127 t + t^2);
m0[t_]:= t(t-12)(51-t)(1884 - 71 t + t^2);
approxw[t_] = Fit[wheat,{w0[t]},t];
approxm[t_] = Fit[maize,{m0[t]},t];
s = Abs[FindMinimum[-approxm[t],{t,37,37.1}][[1]]];
wheatscaled = {{0,0}, {5,6/s}, {10,11/s}, {15,15/s}, {25,18/s},
{30,16/s}, {35,14/s}, {40,10/s}, {45,6/s}, {51,0}};
maizescaled = {{12,0}, {15,9/s}, {20,20/s}, {25,29/s},
{30,35/s}, {35,38/s}, {45,30/s}, {50,7/s}, {51,0}};
lpw = ListPlot[wheatscaled];
lpm = ListPlot[maizescaled];
appw[t_] = approxw[t]/s;
appm[t_] = approxm[t]/s;
Print[Factor[appw[t]]];
Print[Factor[appm[t]]];
Print[s];
plot1 = Plot[appw[t],{t,0,51}];
plot2 = Plot[appm[t],{t,12,51}];
upper = Show[plot1,lpw, AxesLabel->{"T"," "},GridLines->Automatic,
PlotLabel -> "(a) Net photosynthesis versus temperature in wheat",
PlotRange->{{0,52},{0,1.05}}];
lower = Show[plot2,lpm,AxesLabel->{"T"," "},GridLines->Automatic,
PlotLabel -> "(b) Net photosynthesis versus temperature in maize",
PlotRange->{{0,52},{0,1.05}}];
FindMinimum[-appw[t],{t,20,30}]
FindMinimum[-appm[t],{t,37,37.1}]
Show[GraphicsArray[{{upper},{lower}}],PlotLabel->"Figure 8.1"]
^*)