(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 9.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 157, 7]
NotebookDataLength[ 54532, 1268]
NotebookOptionsPosition[ 51825, 1171]
NotebookOutlinePosition[ 52205, 1188]
CellTagsIndexPosition[ 52162, 1185]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["Graphing and Curve Fitting", "Title",
CellChangeTimes->{{3.5927605207222347`*^9, 3.592760525364723*^9}}],
Cell[CellGroupData[{
Cell[TextData[{
"Introduction\n",
StyleBox["\tThis notebook will help you get started with some basic tasks in \
", "Text"],
StyleBox["Mathematica", "Text",
FontSlant->"Italic"],
StyleBox[": importing a text file containing data, doing simple \
manipulations of that list, graphing the data, and performing simple curve \
fitting.\n\tThe first cell below is how you import the tab-delimited text \
file (which is what you get if you save a spreadsheet file as .txt) that \
contains your data. Note that I\[CloseCurlyQuote]ve chosen to give the \
imported data the name ", "Text"],
StyleBox["rawData", "Text",
FontColor->RGBColor[0, 0, 1]],
StyleBox[", and ", "Text"],
StyleBox["Mathematica", "Text",
FontSlant->"Italic"],
StyleBox[" will refer to it by that name from now on.", "Text"]
}], "Section",
CellChangeTimes->{{3.5927605432845783`*^9, 3.5927607574195223`*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"rawData", "=",
RowBox[{"Import", "[",
RowBox[{"\"\\"", ",", "\"\
\""}],
"]"}]}]], "Input",
CellChangeTimes->{{3.592743233708932*^9, 3.5927432585574083`*^9}, {
3.592743353798915*^9, 3.5927433589778223`*^9}, {3.5927608494740143`*^9,
3.592760863423894*^9}, {3.5927614603715897`*^9, 3.592761476199666*^9}, {
3.72657366237144*^9, 3.726573663369505*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "2.688094259`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "6.651050631`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "10.97636448`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"3", ",", "12.93577531`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"4", ",", "14.60505538`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"5", ",", "15.27245987`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"6", ",", "14.0891266`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"7", ",", "12.7060682`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"8", ",", "10.44592111`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"9", ",", "6.65659726`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"10", ",", "12"}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{3.59276153727271*^9, 3.726573575311041*^9,
3.72657366560567*^9}]
}, Open ]],
Cell[TextData[{
"The first argument to the Import command is the full path to the data file \
you wish to work on. The one I sent with this notebook is called \
\[OpenCurlyDoubleQuote]TestData.txt\[CloseCurlyDoubleQuote]. You can get the \
path to the file\[CloseCurlyQuote]s location on your machine by using the \
\[OpenCurlyDoubleQuote]Insert...File Path\[CloseCurlyDoubleQuote] from the \
menu bar and selecting the file in the dialog box that comes up. Note that, \
unless you put a semicolon at the end of a ",
StyleBox["Mathematica",
FontSlant->"Italic"],
" command, you will get some kind of output, even if it\[CloseCurlyQuote]s \
only en echo of the input line. In this case, it\[CloseCurlyQuote]s a \
listing of rawData. "
}], "Text",
CellChangeTimes->{{3.5927609330207663`*^9, 3.592761074318637*^9}, {
3.592761488422447*^9, 3.592761522788701*^9}, {3.592761566991819*^9,
3.5927616327844276`*^9}},
TextAlignment->Left],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"graph1", "=",
RowBox[{"ListPlot", "[",
RowBox[{"rawData", ",", " ",
RowBox[{"PlotMarkers", "\[Rule]",
RowBox[{"{",
RowBox[{"Automatic", ",", " ", "Medium"}], "}"}]}], ",", " ",
RowBox[{"Frame", "\[Rule]", "True"}], ",",
RowBox[{"FrameLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"Independent", " ", "Variable"}], ",", " ",
RowBox[{"Dependent", " ", "Variable"}]}], "}"}]}], ",",
RowBox[{"PlotLabel", "\[Rule]",
RowBox[{"Style", "[",
RowBox[{"\"\\"", ",", "Blue", ",", "Large", ",",
RowBox[{"FontFamily", "\[Rule]", "Times"}]}], "]"}]}]}],
"]"}]}]], "Input",
CellChangeTimes->{{3.592743371284417*^9, 3.592743381465045*^9}, {
3.592761782500079*^9, 3.592761816137885*^9}, {3.592761937782362*^9,
3.5927619483007107`*^9}, {3.5927620980639563`*^9, 3.592762103879722*^9}, {
3.592762159603216*^9, 3.592762176197979*^9}, {3.592762243690099*^9,
3.592762245694525*^9}, {3.5927629925094347`*^9, 3.592763029585919*^9}, {
3.592763275136916*^9, 3.59276330951936*^9}, {3.5927633434379663`*^9,
3.592763346349267*^9}, {3.592763378812046*^9, 3.592763414418799*^9}, {
3.592763453344495*^9, 3.5927635186704197`*^9}}],
Cell[BoxData[
GraphicsBox[{{}, {
{RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[6],
AbsoluteThickness[1.6], GeometricTransformationBox[InsetBox[
StyleBox["\<\"\[FilledCircle]\"\>",
StripOnInput->False,
FontSize->Medium], {0., 0.}], {{{0., 2.688094259}}, {{1.,
6.651050631}}, {{2., 10.97636448}}, {{3., 12.93577531}}, {{4.,
14.60505538}}, {{5., 15.27245987}}, {{6., 14.0891266}}, {{7.,
12.7060682}}, {{8., 10.44592111}}, {{9., 6.65659726}}, {{10.,
12.}}}]}, {}}, {}, {}, {{}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0, 0},
DisplayFunction->Identity,
Frame->{{True, True}, {True, True}},
FrameLabel->{{
FormBox[
RowBox[{"Dependent", " ", "Variable"}], TraditionalForm], None}, {
FormBox[
RowBox[{"Independent", " ", "Variable"}], TraditionalForm], None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotLabel->FormBox[
StyleBox["\"A Typical Plot\"",
RGBColor[0, 0, 1], Large, FontFamily -> Times, StripOnInput -> False],
TraditionalForm],
PlotRange->{{0, 10.}, {0, 15.27245987}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.59274338282514*^9, 3.5927618215932837`*^9, 3.592761950664105*^9,
3.592762106235001*^9, 3.592762177926197*^9, 3.592763033055614*^9,
3.592763311275777*^9, 3.5927633881546507`*^9, {3.592763466382751*^9,
3.592763519825111*^9}, 3.726573672991252*^9}]
}, Open ]],
Cell[TextData[{
"Note that you use ListPlot to make a scatter plot of data; Plot is used to \
graph a mathematical function. I\[CloseCurlyQuote]ve added a few bells and \
whistles just to jazz it up; feel free to consult the help pages for more. \
To make a bare-bones plot, you really only need ListPlot[rawData]. By the \
way, you can make that right arrow symbol with the Typesetting menu on the \
Writing Assistant palette, or you can just type - followed by > and let ",
StyleBox["Mathematica",
FontSlant->"Italic"],
" put them together."
}], "Text",
CellChangeTimes->{{3.592761326237598*^9, 3.592761426480568*^9}, {
3.592763559651469*^9, 3.59276361235321*^9}, {3.5927636548888063`*^9,
3.592763708325561*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"fixedData", "=",
RowBox[{"Take", "[",
RowBox[{"rawData", ",", "10"}], "]"}]}]], "Input",
CellChangeTimes->{{3.5927433971051197`*^9, 3.592743423494978*^9}, {
3.5927434596695623`*^9, 3.592743459812367*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "2.688094259`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "6.651050631`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"2", ",", "10.97636448`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"3", ",", "12.93577531`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"4", ",", "14.60505538`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"5", ",", "15.27245987`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"6", ",", "14.0891266`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"7", ",", "12.7060682`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"8", ",", "10.44592111`"}], "}"}], ",",
RowBox[{"{",
RowBox[{"9", ",", "6.65659726`"}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{3.592743424862162*^9, 3.592743461711575*^9,
3.726573677748814*^9}]
}, Open ]],
Cell[TextData[{
"Since that last point is obviously an outlier, we used the Take function to \
take the first 10 data point pairs (thus leaving out that last suspicious \
point) and made a new list called ",
StyleBox["fixedData",
FontColor->RGBColor[0, 0, 1]],
". Here\[CloseCurlyQuote]s the graph."
}], "Text",
CellChangeTimes->{{3.592763767643784*^9, 3.592763858447545*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"graph2", "=",
RowBox[{"ListPlot", "[", "fixedData", "]"}]}]], "Input",
CellChangeTimes->{{3.592743428716177*^9, 3.592743437558111*^9}}],
Cell[BoxData[
GraphicsBox[{{}, {{},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
AbsoluteThickness[1.6],
PointBox[{{0., 2.688094259}, {1., 6.651050631}, {2., 10.97636448}, {3.,
12.93577531}, {4., 14.60505538}, {5., 15.27245987}, {6., 14.0891266}, {
7., 12.7060682}, {8., 10.44592111}, {9.,
6.65659726}}]}, {}}, {}, {}, {{}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0, 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotRange->{{0, 9.}, {0, 15.27245987}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.592743470765345*^9, 3.726573682024439*^9}]
}, Open ]],
Cell[TextData[{
"Since the graph looks vaguely like a parabola, I\[CloseCurlyQuote]m going \
to choose that functional form to fit the data. First I define a function \
named ",
StyleBox["fitFun",
FontColor->RGBColor[0, 0, 1]],
" as a general parabola with parameters A, B, and C."
}], "Text",
CellChangeTimes->{{3.592763871263115*^9, 3.5927639477319307`*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"fitFun", "=",
RowBox[{
RowBox[{"A", " ",
RowBox[{"x", "^", "2"}]}], " ", "+", " ",
RowBox[{"B", " ", "x"}], " ", "+", " ", "C"}]}]], "Input",
CellChangeTimes->{{3.592743544739332*^9, 3.592743593575551*^9}, {
3.592743663988764*^9, 3.592743702403207*^9}, 3.5927640045328817`*^9}],
Cell[BoxData[
RowBox[{"C", "+",
RowBox[{"B", " ", "x"}], "+",
RowBox[{"A", " ",
SuperscriptBox["x", "2"]}]}]], "Output",
CellChangeTimes->{3.7265736899750423`*^9}]
}, Open ]],
Cell[TextData[{
"Now I\[CloseCurlyQuote]m going to let ",
StyleBox["Mathematica",
FontSlant->"Italic"],
" perform the curve fitting and put the best fit parameters in a list called \
soln (so I can refer to the three paramater values in a shorthand way later). \
Note that FindFit wants four inputs: the data list to fit, the function to \
fit with, the list of adjustable parameters, and the name(s) of the \
independent variable(s)."
}], "Text",
CellChangeTimes->{{3.592764025673869*^9, 3.59276417429897*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"soln", "=",
RowBox[{"FindFit", "[",
RowBox[{"fixedData", ",", "fitFun", ",",
RowBox[{"{",
RowBox[{"A", ",", "B", ",", "C"}], "}"}], ",", "x"}], "]"}]}]], "Input",
CellChangeTimes->{{3.592743544739332*^9, 3.592743593575551*^9}, {
3.592743663988764*^9, 3.592743702403207*^9}, 3.5927640045328817`*^9}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"A", "\[Rule]",
RowBox[{"-", "0.5075951015681828`"}]}], ",",
RowBox[{"B", "\[Rule]", "5.023244382731825`"}], ",",
RowBox[{"C", "\[Rule]", "2.564511982399986`"}]}], "}"}]], "Output",
CellChangeTimes->{3.59274359476974*^9, 3.5927437056267033`*^9,
3.726573694920907*^9}]
}, Open ]],
Cell[TextData[{
"So now my best fit parameters are known, but I\[CloseCurlyQuote]d like ",
StyleBox["Mathematica",
FontSlant->"Italic"],
" to construct the finished function for me as a complete function of x. \
The /. command means \[OpenCurlyDoubleQuote]evaluated at\
\[CloseCurlyDoubleQuote], so the following command tells ",
StyleBox["Mathematica",
FontSlant->"Italic"],
" to take the general function ",
StyleBox["fitFun",
FontColor->RGBColor[0, 0, 1]],
" and put in the parameter values found in ",
StyleBox["soln",
FontColor->RGBColor[0, 0, 1]],
"."
}], "Text",
CellChangeTimes->{{3.59276419421035*^9, 3.5927642977182093`*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"y", "=",
RowBox[{"fitFun", "/.", "soln"}]}]], "Input",
CellChangeTimes->{{3.5927436273723383`*^9, 3.592743629141824*^9}, {
3.592743713707353*^9, 3.5927437193388844`*^9}}],
Cell[BoxData[
RowBox[{"2.564511982399986`", "\[VeryThinSpace]", "+",
RowBox[{"5.023244382731825`", " ", "x"}], "-",
RowBox[{"0.5075951015681828`", " ",
SuperscriptBox["x", "2"]}]}]], "Output",
CellChangeTimes->{3.592743721380932*^9, 3.726573699468677*^9}]
}, Open ]],
Cell["\<\
Now graph it using Plot, since it\[CloseCurlyQuote]s a function. Note that I \
give it a list in curly brackets containing the name of the independent \
variable and the starting and ending values for the plot.\
\>", "Text",
CellChangeTimes->{{3.5927643325816193`*^9, 3.592764374298957*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"graph3", "=",
RowBox[{"Plot", "[",
RowBox[{"y", ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "9"}], "}"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.592743734573678*^9, 3.592743761441757*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwV13k81NsbB/CpUEnMzPc7Y48RyfJFuYXCfORWbqvKUq4oIVSWSHu6SJay
RJRUosVWhChEVLYo6ldIiy0hoZKmlH6nv+b1fn23c55znuc5I3DxWe82mcVi
TZvEYv35HfXTGBRdDTc/vyPSXiVUDK5Ryk6K592FoQMdL4NTxSBfoMMeSAkU
qm+Sse4qE8No9XzP5pQwYYz4qMzpr2LotC7iPUxJFPLPvXpX5y6OVp7oQ3nK
VWHXbe0BXysJbBDJtZSmFAn3VJTYFE6bhjjPzFUlKdXCY6LMm+zDMxASevqh
XmyLUEZ7xGXnVg72qkxXn/Zfp7C9ZO8Tw8s0VihKuZsqvxdOD7l/2bZIFsum
CZbrGH4UHlA//L9vvQpY/bVfKs73k5Dbl7Bbv0IZUs4NDz1DR4UO4ptz8+1V
oSDHEfZEfBO2X/2WaVokwP8eqd90+feH8KLharXgQTXM3le5u9zzlzD4cv/+
Em11xMQk8mXUWLi3T9U7b4kGCiN/PZMJn4R/V1+4LL90DmZP1jIJkpqCM3lh
VMhKTfSNlSaPBIrB0y7E6aj+XHzbTl1u+i2ONc27GCGthRuslJP//TcV0eui
Phl2aaGUf8L/r/5pMBmqi8uK0wY3E5PsfCQxSyLTqs1MBxW3HaP2t89AoPBc
VGSjDor/Nrs7w30m3liyfW/a6cJOtOaJd7806tO7JBLqdaHwOnwbPSyN6Z/r
lwgadPGmsurbnVFpOBqWHbzeqAu3qIWq4iwZzI+501PVpIuAWbP8zvFlcCVb
Jb7zhS7i//7IrbWUQcnes1uHunTRFBtlJ7ggg2OOlg+X/tDFFbpnyFOLjayN
vtc+aDBgxVon7tdlQ0zZWd5Gk4HDjLtmEfpszL3Q7VA6l4HM5MQTGQvYePTi
+qowHQb7h5drv7dgI0iw/MPMeQxW12e7uTqwcVPpJN2zmMHXIP9XTlFsCLLc
ZXOsGViPvw3xjmYjY+rMQ0PrGGQFrtI5EsfGd8E/efobGDjv0jhwPomN0RPB
x2/YMqhzaJF9lc6Gk+mPkfMODM4vWLxhUxkbItavXEtXBks/TK7f8JENE4H1
m6xABskLPPboj7Bh/MXLpXAvg5GgRtUZX9hIK9d8XbaPXKfP7K0SsVFgcbyo
/gCDITNGY54YBy2i16K2IwzOxNoflVbkgKMzUFodxmBgQc7COisOnC4eXtx2
moHwKKfr8koOsi+0xFcmMkioDzx5dA0HbknPujKSyHVnix4jGw7oPTddA86S
6+Ev4q46c+DTKSaaSGFg1s76GBLIgcoRyc1d6Qxij9qmm1/mYF1NrGpXHoO/
mIMjM65xcHyat/uVmwza2lLN2jI52KwrvL49n4Ga4WCLfy4HY4p3zfoLGBS+
C5bKLOGgOVnTq6OIQeuKvD1UM3m/bPHIhTIGKrwZVn2/ONCTvFfoWMPgQaXB
6VssLhouysuzahl4eNt1BU/hYobc9iNpxHk1lw4pT+diT97rlT11DJYcML65
nuZi5KjJ538bGLi/dVO4q83Ff6nZOxWbGdzIvDcYZ8+F4sv9cu1tDFa+XSCb
4sDF2+6ltz1fMuijsy2ubuaiY8t7nW/EakdPJ5Zs4+J0YVa09CsGiXZe6Pbh
4n7V/CX6bxgETaET/grnYtkBq63GXQyUTCLKzaO4qCreYV5IXOI90WcVzcV+
D9dhvW4Go219ZpsTuDg7+0WDag+ZT97d98dSuUhjle4Ze0f22+bti1tuc/Hc
NjLIs5/Bx1Ov3DpLuZhaI5b8lDiydl3sh3Iy32yb8EUDJD6Gpu9+P+Did2NA
h8QHBsaSnJi5zVx07S86EjfIQFBc0rW/n4uHO4X+3sMMygf1Z4YOcpHD4Tx6
QOyodsUoepiLJ1q3x+RHyPxOxpxI+8pFbUFY3T1iSVfXhY8mUQjW09jM+szg
M1s6UkmBAr9iu+nyUTIe3qHbn5UoyErOtg8jVlMY6K1VoXAoZe/OB8TrZldb
7tGgsDdpMMzsK1mvBUd+NRpQ0OwVPzh7jIHVokGdy4YU1uyM8nIg7jB3cDiw
kMIrkyybWGIZq4XFc0wpqGcXyH0n9nYY8g1aTsG/96tz+TcGEs6OF21XULh/
ban6R+IL2+obdVZT2Jen160gYtC486p263oK7DMN1v7EukFOPQZOFH5VGH2Q
/U7iF9LAnbqVQkDZulBzYsfwRRavt1FYyxKXcyWOiuNfiPCkYORkpJ9DPJD+
2L4rgMK93Y+UDH4wCM4wDbu9l8KWrNzINcQK17MKow9Q2KbFfNlB/E9RGGfx
UQrJv74UpBN3lIwKOSEUGtabSZUT76tw8X5/jNw/rOLcQpxRa/4oPopC5UTe
qMQ4AzTmfPeMpnDcz9lImbi1WWEu4igU6FwLmE88tX0sdDCRjC9w8ttNxBfe
uhZUnaWQIMGX2km8sOdp55kUCuveVBkeJnb/eMN8aRqFiW19u1OIJz4p7VK8
QmGQszUyizhxLPLcp2sUMNU/pZhYb1xUV5NFYcxyVuZ94oe/3UXnr1Pg1rrk
PSZ2FHs+JyCPfC8e+a3Eo9MsbVcUkPW/Wnyjk/jEzJshqkUUqlm1V/qJ1bkq
+WO3yfqk+SWNEJfxT3Y0lFI4GVUYMka8QXFcOr2cwtSyWK9x4gEVT7P9lRR6
DSVW/SYOVm/ZsfYBBe0v9NzJP0m8tZYma9RQePa9bGIKcT5TUDteR6Fj+cQT
MeIV8wXfmhso2L5uPffHnQtjNDKeUNAtWeXy5/7dD1bNK31K4aL0ivcs4snr
Jc0ePyfzkzEx+UW+d+ptjVVnK4XwUvmNImK1XcdsRttJ/KhBy8/E+T+WbJn6
lsJuubyfA8SW4aydCl0U+mvdgruI/8cr38u8o+CsMfP5n/i4ph8MQR+FVubq
aOOfeBmYxGz4QGF1l0FnJXFo+Viy+xCJ15LryYXE9KrCq/s/UTBfqyi4Snyl
zS//xCiFoqkHAxL/rO92/fKL3ygY7GxIOEZcPTpYl/+DQvYhmWB/4l4Zj85W
Fo3US173VxIHntf4+GEKDYWiYNmFxBI63aIJCRpJgSeMVYg1lzmxNWbSiBkP
6PxI9uPtZ4pKxmwaUvc2eD4jttraprmSohFloF5eTOxxyEboJ0/jOC/58SFi
0XTuylAlGlMjhZGbiSOSntglqZDrWS+4ZsRZ+Su872oQD7dFikg+fejD+enz
aURH/MXfSnwocCJDaQENGaftJsbEUmJlhfrGNBLeRGlJE+uqGDXYmtM4153p
W0Ty2duWGb+0gsbNXybJIpL/E10DErdW0zCZ/0O9mjjaN4Nba02jgptz6BRx
btRs7SE7Gk0tPac0iUcq5Tct2kZj1C/LwZLUl6NrW1xXu9MIkhCtmkrMfp3g
u8WTxozlRlL1pH7NE8mEh/nQUGuI/LSSeLeeRPHTgzTse4/ZW5B69zX5M70j
gUaLQ0Ps509k/SN9q3ySaJw2XlNzmbjwwJBPQDIN6XOvn9sSBzgM1B9OpWEw
pH6wkNTbz/JdR2NzaBg32vZuI/V5+Ezz4K0HNBy8IYwi9fxx+NrkkhoaNUnJ
Y9rEN/Y1Lq+opzG8dHJYHan/3hvr0mqbaPgZ/NCdQjwoW7mx/RUNI/MT5zxI
/+hPzHvA+kqj2dRkfJz0m9owvd3iIhpDN3cYRhBnBOaoSI7TqM7OtuQRe9hl
HKAm8XCZsRnUIv2ql5dqMGcmD5pey52Wkv7WnRCTslKDB9Pqogd/k374+pR3
QKItDxcl75xze072z7p3Gts38iCXOUy3/o/kE9uxxehfHiYZ5574h3hZ9D+L
2rbwQGvW7dV6RvpDhDpLaScPDis8TF83kXw/+jLqUjAP9lc2BNGkn8/xXpae
k8uD7PU3HyoryfmPuWtzOJ+HeZsMF9LELz8YSqy5xcOdr+wgt3sMYjwFnsMl
PLBZ26XEK8j+dP2pN7+ah5ebjvCMyXmj/t/8kuJXPLwROHTuIucRn39mPa2a
zkf2lGpnuSxynojr2dIkxcdvk5lSJpkkH9uyhl/L8GFyIjJiYwbp7x5GM7/z
+PAK1RTFXyX5FG5tpa/GR8b9gfSfaeT8Vx1cfm4RH7mHExecPsfgkeX7rIAd
fDS/DbbojiL1Ouq6SbA3H1rxJ+k3kQx8n/nXxvjx4eHy+HZLBAP+NlZv1l4+
Um0+B9QcZ+ASLC/oDOFj5Pj90jMhDL7fW5m0OoUPvnZUyOSDJH/N80LmNPJx
wDF3/LMHqS82PnxWEx+DBVFd7dsZVO3Qy2x7ykdKaXTlfXcGTmdzHp9o5eNG
+zSPWHI+TRjNkP/STcbv1GinvIXEN/tSbsUP8rxPUqm4Pek/svHtG+fKIqdC
P439N+n/n/YYRoXIgk/7dUsrMrjm9XW7XZgsavfPSyuTJ/21OyBFECGLH89c
tnnKkfx87i9+O1oWA25yHyt4DPxK/Fp6zspi6Pdj5a1sBhtDdh00z5OFr5TX
qL84Aw2e2/1Pr2QRdtw+sH1IF/eMbTZsXCiHVnF2481yXSiV0zWOIjkEHM+b
iLTSJf+Hppu4XJeHdCpyDtTpYG3o24kLfgrwvajY5MLoYM7D/qpXcxXR+3id
jVukNprOSn0LfqUIry8bN/i91EJqDidoy0UlmDYPJgdN14LikRXiQltlSPpF
j/y3YC5KT7q9tVGYhd3abj//NtOE/ryzWs/aZ6FBLkv52cI5WBPP2r0+VgVp
tktKl83TQJ28xnExbVXgQIeC9xR1WNRts9lyTBVbKOudh1vVkF51Y9PIG1Xs
Y7TLHdMFkG0fmKYEAQL8y4/8tFBF51MY1SYLIDqSa1D0WBkpc1fyN6cI4KE9
lGBTpYyNR2y/fjovwBllNZ+RW8po0vIqULwkwL52z2Oq55VRfjRe3+eaAMa3
9pm47FBGsl6vJr9QAIOO+IodU5WxITJKzrWRmL5nmWSmBOmO099Ej8n37lxz
7zFQQv2C1BcnmwQIvxb6W09dCUs6CxOKnwkg5zArvkxSCQbGb9gzXgoA5WaL
Oy2KkOo1mJ7fS55X+OuusZ8iahcv7lvWJ0CHma9ws6siQuKW1rT3C5CRccko
yF4RP0wdjol/FCDVtiei1EwRffEhkzZ9EcB3Zv8tnqQiLvef7Pg4St7n/eS9
xi8FOAvPVASPCWAtlx1jOKKAFwM5h3O+k/uXLLFY/kIBp1DsaDFOxnP/14rV
dQpYk1i5+MVPAVqjbtRblylAcvCRwo4JAaxSbO6ty1XAQ4sX33//FiB2dFjH
Ok0B/we/gERu
"]]}}, {}, {}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0, 2.564512905036692},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotRange->{{0, 9}, {2.564512905036692, 14.992224216261036`}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{{3.592743747866268*^9, 3.592743762731365*^9},
3.726573706379355*^9}]
}, Open ]],
Cell["\<\
Not I\[CloseCurlyQuote]d like to see how well I did, so I use the Show \
command to superimpose the last two plots.\
\>", "Text",
CellChangeTimes->{{3.592764391306087*^9, 3.592764430128737*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Show", "[",
RowBox[{"graph2", ",", "graph3"}], "]"}]], "Input",
CellChangeTimes->{{3.5927437673194027`*^9, 3.5927437773936777`*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {{},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[
0.012833333333333334`], AbsoluteThickness[1.6],
PointBox[{{0., 2.688094259}, {1., 6.651050631}, {2., 10.97636448}, {3.,
12.93577531}, {4., 14.60505538}, {5., 15.27245987}, {6., 14.0891266}, {
7., 12.7060682}, {8., 10.44592111}, {9.,
6.65659726}}]}, {}}, {}, {}, {{}, {}}}, {{{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwV13k81NsbB/CpUEnMzPc7Y48RyfJFuYXCfORWbqvKUq4oIVSWSHu6SJay
RJRUosVWhChEVLYo6ldIiy0hoZKmlH6nv+b1fn23c55znuc5I3DxWe82mcVi
TZvEYv35HfXTGBRdDTc/vyPSXiVUDK5Ryk6K592FoQMdL4NTxSBfoMMeSAkU
qm+Sse4qE8No9XzP5pQwYYz4qMzpr2LotC7iPUxJFPLPvXpX5y6OVp7oQ3nK
VWHXbe0BXysJbBDJtZSmFAn3VJTYFE6bhjjPzFUlKdXCY6LMm+zDMxASevqh
XmyLUEZ7xGXnVg72qkxXn/Zfp7C9ZO8Tw8s0VihKuZsqvxdOD7l/2bZIFsum
CZbrGH4UHlA//L9vvQpY/bVfKs73k5Dbl7Bbv0IZUs4NDz1DR4UO4ptz8+1V
oSDHEfZEfBO2X/2WaVokwP8eqd90+feH8KLharXgQTXM3le5u9zzlzD4cv/+
Em11xMQk8mXUWLi3T9U7b4kGCiN/PZMJn4R/V1+4LL90DmZP1jIJkpqCM3lh
VMhKTfSNlSaPBIrB0y7E6aj+XHzbTl1u+i2ONc27GCGthRuslJP//TcV0eui
Phl2aaGUf8L/r/5pMBmqi8uK0wY3E5PsfCQxSyLTqs1MBxW3HaP2t89AoPBc
VGSjDor/Nrs7w30m3liyfW/a6cJOtOaJd7806tO7JBLqdaHwOnwbPSyN6Z/r
lwgadPGmsurbnVFpOBqWHbzeqAu3qIWq4iwZzI+501PVpIuAWbP8zvFlcCVb
Jb7zhS7i//7IrbWUQcnes1uHunTRFBtlJ7ggg2OOlg+X/tDFFbpnyFOLjayN
vtc+aDBgxVon7tdlQ0zZWd5Gk4HDjLtmEfpszL3Q7VA6l4HM5MQTGQvYePTi
+qowHQb7h5drv7dgI0iw/MPMeQxW12e7uTqwcVPpJN2zmMHXIP9XTlFsCLLc
ZXOsGViPvw3xjmYjY+rMQ0PrGGQFrtI5EsfGd8E/efobGDjv0jhwPomN0RPB
x2/YMqhzaJF9lc6Gk+mPkfMODM4vWLxhUxkbItavXEtXBks/TK7f8JENE4H1
m6xABskLPPboj7Bh/MXLpXAvg5GgRtUZX9hIK9d8XbaPXKfP7K0SsVFgcbyo
/gCDITNGY54YBy2i16K2IwzOxNoflVbkgKMzUFodxmBgQc7COisOnC4eXtx2
moHwKKfr8koOsi+0xFcmMkioDzx5dA0HbknPujKSyHVnix4jGw7oPTddA86S
6+Ev4q46c+DTKSaaSGFg1s76GBLIgcoRyc1d6Qxij9qmm1/mYF1NrGpXHoO/
mIMjM65xcHyat/uVmwza2lLN2jI52KwrvL49n4Ga4WCLfy4HY4p3zfoLGBS+
C5bKLOGgOVnTq6OIQeuKvD1UM3m/bPHIhTIGKrwZVn2/ONCTvFfoWMPgQaXB
6VssLhouysuzahl4eNt1BU/hYobc9iNpxHk1lw4pT+diT97rlT11DJYcML65
nuZi5KjJ538bGLi/dVO4q83Ff6nZOxWbGdzIvDcYZ8+F4sv9cu1tDFa+XSCb
4sDF2+6ltz1fMuijsy2ubuaiY8t7nW/EakdPJ5Zs4+J0YVa09CsGiXZe6Pbh
4n7V/CX6bxgETaET/grnYtkBq63GXQyUTCLKzaO4qCreYV5IXOI90WcVzcV+
D9dhvW4Go219ZpsTuDg7+0WDag+ZT97d98dSuUhjle4Ze0f22+bti1tuc/Hc
NjLIs5/Bx1Ov3DpLuZhaI5b8lDiydl3sh3Iy32yb8EUDJD6Gpu9+P+Did2NA
h8QHBsaSnJi5zVx07S86EjfIQFBc0rW/n4uHO4X+3sMMygf1Z4YOcpHD4Tx6
QOyodsUoepiLJ1q3x+RHyPxOxpxI+8pFbUFY3T1iSVfXhY8mUQjW09jM+szg
M1s6UkmBAr9iu+nyUTIe3qHbn5UoyErOtg8jVlMY6K1VoXAoZe/OB8TrZldb
7tGgsDdpMMzsK1mvBUd+NRpQ0OwVPzh7jIHVokGdy4YU1uyM8nIg7jB3cDiw
kMIrkyybWGIZq4XFc0wpqGcXyH0n9nYY8g1aTsG/96tz+TcGEs6OF21XULh/
ban6R+IL2+obdVZT2Jen160gYtC486p263oK7DMN1v7EukFOPQZOFH5VGH2Q
/U7iF9LAnbqVQkDZulBzYsfwRRavt1FYyxKXcyWOiuNfiPCkYORkpJ9DPJD+
2L4rgMK93Y+UDH4wCM4wDbu9l8KWrNzINcQK17MKow9Q2KbFfNlB/E9RGGfx
UQrJv74UpBN3lIwKOSEUGtabSZUT76tw8X5/jNw/rOLcQpxRa/4oPopC5UTe
qMQ4AzTmfPeMpnDcz9lImbi1WWEu4igU6FwLmE88tX0sdDCRjC9w8ttNxBfe
uhZUnaWQIMGX2km8sOdp55kUCuveVBkeJnb/eMN8aRqFiW19u1OIJz4p7VK8
QmGQszUyizhxLPLcp2sUMNU/pZhYb1xUV5NFYcxyVuZ94oe/3UXnr1Pg1rrk
PSZ2FHs+JyCPfC8e+a3Eo9MsbVcUkPW/Wnyjk/jEzJshqkUUqlm1V/qJ1bkq
+WO3yfqk+SWNEJfxT3Y0lFI4GVUYMka8QXFcOr2cwtSyWK9x4gEVT7P9lRR6
DSVW/SYOVm/ZsfYBBe0v9NzJP0m8tZYma9RQePa9bGIKcT5TUDteR6Fj+cQT
MeIV8wXfmhso2L5uPffHnQtjNDKeUNAtWeXy5/7dD1bNK31K4aL0ivcs4snr
Jc0ePyfzkzEx+UW+d+ptjVVnK4XwUvmNImK1XcdsRttJ/KhBy8/E+T+WbJn6
lsJuubyfA8SW4aydCl0U+mvdgruI/8cr38u8o+CsMfP5n/i4ph8MQR+FVubq
aOOfeBmYxGz4QGF1l0FnJXFo+Viy+xCJ15LryYXE9KrCq/s/UTBfqyi4Snyl
zS//xCiFoqkHAxL/rO92/fKL3ygY7GxIOEZcPTpYl/+DQvYhmWB/4l4Zj85W
Fo3US173VxIHntf4+GEKDYWiYNmFxBI63aIJCRpJgSeMVYg1lzmxNWbSiBkP
6PxI9uPtZ4pKxmwaUvc2eD4jttraprmSohFloF5eTOxxyEboJ0/jOC/58SFi
0XTuylAlGlMjhZGbiSOSntglqZDrWS+4ZsRZ+Su872oQD7dFikg+fejD+enz
aURH/MXfSnwocCJDaQENGaftJsbEUmJlhfrGNBLeRGlJE+uqGDXYmtM4153p
W0Ty2duWGb+0gsbNXybJIpL/E10DErdW0zCZ/0O9mjjaN4Nba02jgptz6BRx
btRs7SE7Gk0tPac0iUcq5Tct2kZj1C/LwZLUl6NrW1xXu9MIkhCtmkrMfp3g
u8WTxozlRlL1pH7NE8mEh/nQUGuI/LSSeLeeRPHTgzTse4/ZW5B69zX5M70j
gUaLQ0Ps509k/SN9q3ySaJw2XlNzmbjwwJBPQDIN6XOvn9sSBzgM1B9OpWEw
pH6wkNTbz/JdR2NzaBg32vZuI/V5+Ezz4K0HNBy8IYwi9fxx+NrkkhoaNUnJ
Y9rEN/Y1Lq+opzG8dHJYHan/3hvr0mqbaPgZ/NCdQjwoW7mx/RUNI/MT5zxI
/+hPzHvA+kqj2dRkfJz0m9owvd3iIhpDN3cYRhBnBOaoSI7TqM7OtuQRe9hl
HKAm8XCZsRnUIv2ql5dqMGcmD5pey52Wkv7WnRCTslKDB9Pqogd/k374+pR3
QKItDxcl75xze072z7p3Gts38iCXOUy3/o/kE9uxxehfHiYZ5574h3hZ9D+L
2rbwQGvW7dV6RvpDhDpLaScPDis8TF83kXw/+jLqUjAP9lc2BNGkn8/xXpae
k8uD7PU3HyoryfmPuWtzOJ+HeZsMF9LELz8YSqy5xcOdr+wgt3sMYjwFnsMl
PLBZ26XEK8j+dP2pN7+ah5ebjvCMyXmj/t/8kuJXPLwROHTuIucRn39mPa2a
zkf2lGpnuSxynojr2dIkxcdvk5lSJpkkH9uyhl/L8GFyIjJiYwbp7x5GM7/z
+PAK1RTFXyX5FG5tpa/GR8b9gfSfaeT8Vx1cfm4RH7mHExecPsfgkeX7rIAd
fDS/DbbojiL1Ouq6SbA3H1rxJ+k3kQx8n/nXxvjx4eHy+HZLBAP+NlZv1l4+
Um0+B9QcZ+ASLC/oDOFj5Pj90jMhDL7fW5m0OoUPvnZUyOSDJH/N80LmNPJx
wDF3/LMHqS82PnxWEx+DBVFd7dsZVO3Qy2x7ykdKaXTlfXcGTmdzHp9o5eNG
+zSPWHI+TRjNkP/STcbv1GinvIXEN/tSbsUP8rxPUqm4Pek/svHtG+fKIqdC
P439N+n/n/YYRoXIgk/7dUsrMrjm9XW7XZgsavfPSyuTJ/21OyBFECGLH89c
tnnKkfx87i9+O1oWA25yHyt4DPxK/Fp6zspi6Pdj5a1sBhtDdh00z5OFr5TX
qL84Aw2e2/1Pr2QRdtw+sH1IF/eMbTZsXCiHVnF2481yXSiV0zWOIjkEHM+b
iLTSJf+Hppu4XJeHdCpyDtTpYG3o24kLfgrwvajY5MLoYM7D/qpXcxXR+3id
jVukNprOSn0LfqUIry8bN/i91EJqDidoy0UlmDYPJgdN14LikRXiQltlSPpF
j/y3YC5KT7q9tVGYhd3abj//NtOE/ryzWs/aZ6FBLkv52cI5WBPP2r0+VgVp
tktKl83TQJ28xnExbVXgQIeC9xR1WNRts9lyTBVbKOudh1vVkF51Y9PIG1Xs
Y7TLHdMFkG0fmKYEAQL8y4/8tFBF51MY1SYLIDqSa1D0WBkpc1fyN6cI4KE9
lGBTpYyNR2y/fjovwBllNZ+RW8po0vIqULwkwL52z2Oq55VRfjRe3+eaAMa3
9pm47FBGsl6vJr9QAIOO+IodU5WxITJKzrWRmL5nmWSmBOmO099Ej8n37lxz
7zFQQv2C1BcnmwQIvxb6W09dCUs6CxOKnwkg5zArvkxSCQbGb9gzXgoA5WaL
Oy2KkOo1mJ7fS55X+OuusZ8iahcv7lvWJ0CHma9ws6siQuKW1rT3C5CRccko
yF4RP0wdjol/FCDVtiei1EwRffEhkzZ9EcB3Zv8tnqQiLvef7Pg4St7n/eS9
xi8FOAvPVASPCWAtlx1jOKKAFwM5h3O+k/uXLLFY/kIBp1DsaDFOxnP/14rV
dQpYk1i5+MVPAVqjbtRblylAcvCRwo4JAaxSbO6ty1XAQ4sX33//FiB2dFjH
Ok0B/we/gERu
"]]}}, {}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0, 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotRange->{{0, 9.}, {0, 15.27245987}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.592743778726776*^9, 3.726573711541669*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell["\<\
If you want to cut to the chase, follow the steps below. FindFit is good for \
learning, but NonlinearModelFit is the way to go...
\
\>", "Subsubsection",
CellChangeTimes->{{3.726574442498281*^9, 3.726574486615632*^9}}],
Cell["\<\
But if you want more quantitative information about how good your fit is, you \
need to learn a slightly more complicated (but much more powerful) fit \
routine: NonlinearModelFit
\
\>", "Text",
CellChangeTimes->{{3.7265734706338453`*^9, 3.726573512519113*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"soln2", "=",
RowBox[{"NonlinearModelFit", "[",
RowBox[{"fixedData", ",", "fitFun", ",",
RowBox[{"{",
RowBox[{"A", ",", "B", ",", "C"}], "}"}], ",", "x"}], "]"}]}]], "Input",
CellChangeTimes->{{3.72657352270905*^9, 3.726573566429303*^9}}],
Cell[BoxData[
TagBox[
RowBox[{"FittedModel", "[",
TagBox[
PanelBox[
TagBox[
RowBox[{"2.564511982399986`", "\[VeryThinSpace]", "+",
RowBox[{"5.023244382731825`", " ", "x"}], "-",
RowBox[{"0.5075951015681828`", " ",
SuperscriptBox["x", "2"]}]}],
Short[#, 2]& ],
FrameMargins->5],
Editable -> False], "]"}],
InterpretTemplate[
FittedModel[{
"Nonlinear", {$CellContext`A -> -0.5075951015681828, $CellContext`B ->
5.023244382731825, C -> 2.564511982399986}, {{$CellContext`x},
C + $CellContext`B $CellContext`x + $CellContext`A $CellContext`x^2}}, {
1}, {{0, 2.688094259}, {1, 6.651050631}, {2, 10.97636448}, {
3, 12.93577531}, {4, 14.60505538}, {5, 15.27245987}, {6, 14.0891266}, {
7, 12.7060682}, {8, 10.44592111}, {9, 6.65659726}},
Function[Null,
Internal`LocalizedBlock[{$CellContext`A, $CellContext`B,
C, $CellContext`x}, #], {HoldAll}]]& ],
Editable->False,
SelectWithContents->True,
Selectable->True]], "Output",
CellChangeTimes->{3.7265737168488417`*^9}]
}, Open ]],
Cell["\<\
The information returned by NonlinearModelFit is different from that returned \
by FindFit, so you have to do the graph in a slightly different way.\
\>", "Text",
CellChangeTimes->{{3.72657376877359*^9, 3.726573798627615*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"graph4", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"soln2", "[", "x", "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "9"}], "}"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.7265738017114563`*^9, 3.726573824714617*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwV13k81NsbB/CpUEnMzPc7Y48RyfJFuYXCfORWbqvKUq4oIVSWSHu6SJay
RJRUosVWhChEVLYo6ldIiy0hoZKmlH6nv+b1fn23c55znuc5I3DxWe82mcVi
TZvEYv35HfXTGBRdDTc/vyPSXiVUDK5Ryk6K592FoQMdL4NTxSBfoMMeSAkU
qm+Sse4qE8No9XzP5pQwYYz4qMzpr2LotC7iPUxJFPLPvXpX5y6OVp7oQ3nK
VWHXbe0BXysJbBDJtZSmFAn3VJTYFE6bhjjPzFUlKdXCY6LMm+zDMxASevqh
XmyLUEZ7xGXnVg72qkxXn/Zfp7C9ZO8Tw8s0VihKuZsqvxdOD7l/2bZIFsum
CZbrGH4UHlA//L9vvQpY/bVfKs73k5Dbl7Bbv0IZUs4NDz1DR4UO4ptz8+1V
oSDHEfZEfBO2X/2WaVokwP8eqd90+feH8KLharXgQTXM3le5u9zzlzD4cv/+
Em11xMQk8mXUWLi3T9U7b4kGCiN/PZMJn4R/V1+4LL90DmZP1jIJkpqCM3lh
VMhKTfSNlSaPBIrB0y7E6aj+XHzbTl1u+i2ONc27GCGthRuslJP//TcV0eui
Phl2aaGUf8L/r/5pMBmqi8uK0wY3E5PsfCQxSyLTqs1MBxW3HaP2t89AoPBc
VGSjDor/Nrs7w30m3liyfW/a6cJOtOaJd7806tO7JBLqdaHwOnwbPSyN6Z/r
lwgadPGmsurbnVFpOBqWHbzeqAu3qIWq4iwZzI+501PVpIuAWbP8zvFlcCVb
Jb7zhS7i//7IrbWUQcnes1uHunTRFBtlJ7ggg2OOlg+X/tDFFbpnyFOLjayN
vtc+aDBgxVon7tdlQ0zZWd5Gk4HDjLtmEfpszL3Q7VA6l4HM5MQTGQvYePTi
+qowHQb7h5drv7dgI0iw/MPMeQxW12e7uTqwcVPpJN2zmMHXIP9XTlFsCLLc
ZXOsGViPvw3xjmYjY+rMQ0PrGGQFrtI5EsfGd8E/efobGDjv0jhwPomN0RPB
x2/YMqhzaJF9lc6Gk+mPkfMODM4vWLxhUxkbItavXEtXBks/TK7f8JENE4H1
m6xABskLPPboj7Bh/MXLpXAvg5GgRtUZX9hIK9d8XbaPXKfP7K0SsVFgcbyo
/gCDITNGY54YBy2i16K2IwzOxNoflVbkgKMzUFodxmBgQc7COisOnC4eXtx2
moHwKKfr8koOsi+0xFcmMkioDzx5dA0HbknPujKSyHVnix4jGw7oPTddA86S
6+Ev4q46c+DTKSaaSGFg1s76GBLIgcoRyc1d6Qxij9qmm1/mYF1NrGpXHoO/
mIMjM65xcHyat/uVmwza2lLN2jI52KwrvL49n4Ga4WCLfy4HY4p3zfoLGBS+
C5bKLOGgOVnTq6OIQeuKvD1UM3m/bPHIhTIGKrwZVn2/ONCTvFfoWMPgQaXB
6VssLhouysuzahl4eNt1BU/hYobc9iNpxHk1lw4pT+diT97rlT11DJYcML65
nuZi5KjJ538bGLi/dVO4q83Ff6nZOxWbGdzIvDcYZ8+F4sv9cu1tDFa+XSCb
4sDF2+6ltz1fMuijsy2ubuaiY8t7nW/EakdPJ5Zs4+J0YVa09CsGiXZe6Pbh
4n7V/CX6bxgETaET/grnYtkBq63GXQyUTCLKzaO4qCreYV5IXOI90WcVzcV+
D9dhvW4Go219ZpsTuDg7+0WDag+ZT97d98dSuUhjle4Ze0f22+bti1tuc/Hc
NjLIs5/Bx1Ov3DpLuZhaI5b8lDiydl3sh3Iy32yb8EUDJD6Gpu9+P+Did2NA
h8QHBsaSnJi5zVx07S86EjfIQFBc0rW/n4uHO4X+3sMMygf1Z4YOcpHD4Tx6
QOyodsUoepiLJ1q3x+RHyPxOxpxI+8pFbUFY3T1iSVfXhY8mUQjW09jM+szg
M1s6UkmBAr9iu+nyUTIe3qHbn5UoyErOtg8jVlMY6K1VoXAoZe/OB8TrZldb
7tGgsDdpMMzsK1mvBUd+NRpQ0OwVPzh7jIHVokGdy4YU1uyM8nIg7jB3cDiw
kMIrkyybWGIZq4XFc0wpqGcXyH0n9nYY8g1aTsG/96tz+TcGEs6OF21XULh/
ban6R+IL2+obdVZT2Jen160gYtC486p263oK7DMN1v7EukFOPQZOFH5VGH2Q
/U7iF9LAnbqVQkDZulBzYsfwRRavt1FYyxKXcyWOiuNfiPCkYORkpJ9DPJD+
2L4rgMK93Y+UDH4wCM4wDbu9l8KWrNzINcQK17MKow9Q2KbFfNlB/E9RGGfx
UQrJv74UpBN3lIwKOSEUGtabSZUT76tw8X5/jNw/rOLcQpxRa/4oPopC5UTe
qMQ4AzTmfPeMpnDcz9lImbi1WWEu4igU6FwLmE88tX0sdDCRjC9w8ttNxBfe
uhZUnaWQIMGX2km8sOdp55kUCuveVBkeJnb/eMN8aRqFiW19u1OIJz4p7VK8
QmGQszUyizhxLPLcp2sUMNU/pZhYb1xUV5NFYcxyVuZ94oe/3UXnr1Pg1rrk
PSZ2FHs+JyCPfC8e+a3Eo9MsbVcUkPW/Wnyjk/jEzJshqkUUqlm1V/qJ1bkq
+WO3yfqk+SWNEJfxT3Y0lFI4GVUYMka8QXFcOr2cwtSyWK9x4gEVT7P9lRR6
DSVW/SYOVm/ZsfYBBe0v9NzJP0m8tZYma9RQePa9bGIKcT5TUDteR6Fj+cQT
MeIV8wXfmhso2L5uPffHnQtjNDKeUNAtWeXy5/7dD1bNK31K4aL0ivcs4snr
Jc0ePyfzkzEx+UW+d+ptjVVnK4XwUvmNImK1XcdsRttJ/KhBy8/E+T+WbJn6
lsJuubyfA8SW4aydCl0U+mvdgruI/8cr38u8o+CsMfP5n/i4ph8MQR+FVubq
aOOfeBmYxGz4QGF1l0FnJXFo+Viy+xCJ15LryYXE9KrCq/s/UTBfqyi4Snyl
zS//xCiFoqkHAxL/rO92/fKL3ygY7GxIOEZcPTpYl/+DQvYhmWB/4l4Zj85W
Fo3US173VxIHntf4+GEKDYWiYNmFxBI63aIJCRpJgSeMVYg1lzmxNWbSiBkP
6PxI9uPtZ4pKxmwaUvc2eD4jttraprmSohFloF5eTOxxyEboJ0/jOC/58SFi
0XTuylAlGlMjhZGbiSOSntglqZDrWS+4ZsRZ+Su872oQD7dFikg+fejD+enz
aURH/MXfSnwocCJDaQENGaftJsbEUmJlhfrGNBLeRGlJE+uqGDXYmtM4153p
W0Ty2duWGb+0gsbNXybJIpL/E10DErdW0zCZ/0O9mjjaN4Nba02jgptz6BRx
btRs7SE7Gk0tPac0iUcq5Tct2kZj1C/LwZLUl6NrW1xXu9MIkhCtmkrMfp3g
u8WTxozlRlL1pH7NE8mEh/nQUGuI/LSSeLeeRPHTgzTse4/ZW5B69zX5M70j
gUaLQ0Ps509k/SN9q3ySaJw2XlNzmbjwwJBPQDIN6XOvn9sSBzgM1B9OpWEw
pH6wkNTbz/JdR2NzaBg32vZuI/V5+Ezz4K0HNBy8IYwi9fxx+NrkkhoaNUnJ
Y9rEN/Y1Lq+opzG8dHJYHan/3hvr0mqbaPgZ/NCdQjwoW7mx/RUNI/MT5zxI
/+hPzHvA+kqj2dRkfJz0m9owvd3iIhpDN3cYRhBnBOaoSI7TqM7OtuQRe9hl
HKAm8XCZsRnUIv2ql5dqMGcmD5pey52Wkv7WnRCTslKDB9Pqogd/k374+pR3
QKItDxcl75xze072z7p3Gts38iCXOUy3/o/kE9uxxehfHiYZ5574h3hZ9D+L
2rbwQGvW7dV6RvpDhDpLaScPDis8TF83kXw/+jLqUjAP9lc2BNGkn8/xXpae
k8uD7PU3HyoryfmPuWtzOJ+HeZsMF9LELz8YSqy5xcOdr+wgt3sMYjwFnsMl
PLBZ26XEK8j+dP2pN7+ah5ebjvCMyXmj/t/8kuJXPLwROHTuIucRn39mPa2a
zkf2lGpnuSxynojr2dIkxcdvk5lSJpkkH9uyhl/L8GFyIjJiYwbp7x5GM7/z
+PAK1RTFXyX5FG5tpa/GR8b9gfSfaeT8Vx1cfm4RH7mHExecPsfgkeX7rIAd
fDS/DbbojiL1Ouq6SbA3H1rxJ+k3kQx8n/nXxvjx4eHy+HZLBAP+NlZv1l4+
Um0+B9QcZ+ASLC/oDOFj5Pj90jMhDL7fW5m0OoUPvnZUyOSDJH/N80LmNPJx
wDF3/LMHqS82PnxWEx+DBVFd7dsZVO3Qy2x7ykdKaXTlfXcGTmdzHp9o5eNG
+zSPWHI+TRjNkP/STcbv1GinvIXEN/tSbsUP8rxPUqm4Pek/svHtG+fKIqdC
P439N+n/n/YYRoXIgk/7dUsrMrjm9XW7XZgsavfPSyuTJ/21OyBFECGLH89c
tnnKkfx87i9+O1oWA25yHyt4DPxK/Fp6zspi6Pdj5a1sBhtDdh00z5OFr5TX
qL84Aw2e2/1Pr2QRdtw+sH1IF/eMbTZsXCiHVnF2481yXSiV0zWOIjkEHM+b
iLTSJf+Hppu4XJeHdCpyDtTpYG3o24kLfgrwvajY5MLoYM7D/qpXcxXR+3id
jVukNprOSn0LfqUIry8bN/i91EJqDidoy0UlmDYPJgdN14LikRXiQltlSPpF
j/y3YC5KT7q9tVGYhd3abj//NtOE/ryzWs/aZ6FBLkv52cI5WBPP2r0+VgVp
tktKl83TQJ28xnExbVXgQIeC9xR1WNRts9lyTBVbKOudh1vVkF51Y9PIG1Xs
Y7TLHdMFkG0fmKYEAQL8y4/8tFBF51MY1SYLIDqSa1D0WBkpc1fyN6cI4KE9
lGBTpYyNR2y/fjovwBllNZ+RW8po0vIqULwkwL52z2Oq55VRfjRe3+eaAMa3
9pm47FBGsl6vJr9QAIOO+IodU5WxITJKzrWRmL5nmWSmBOmO099Ej8n37lxz
7zFQQv2C1BcnmwQIvxb6W09dCUs6CxOKnwkg5zArvkxSCQbGb9gzXgoA5WaL
Oy2KkOo1mJ7fS55X+OuusZ8iahcv7lvWJ0CHma9ws6siQuKW1rT3C5CRccko
yF4RP0wdjol/FCDVtiei1EwRffEhkzZ9EcB3Zv8tnqQiLvef7Pg4St7n/eS9
xi8FOAvPVASPCWAtlx1jOKKAFwM5h3O+k/uXLLFY/kIBp1DsaDFOxnP/14rV
dQpYk1i5+MVPAVqjbtRblylAcvCRwo4JAaxSbO6ty1XAQ4sX33//FiB2dFjH
Ok0B/we/gERu
"]]}}, {}, {}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0, 2.564512905036692},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{
"DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" ->
AbsolutePointSize[6], "ScalingFunctions" -> None,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotRange->{{0, 9}, {2.564512905036692, 14.992224216261036`}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.726573825723515*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Show", "[",
RowBox[{"graph2", ",", "graph4"}], "]"}]], "Input",
CellChangeTimes->{{3.7265738308310823`*^9, 3.7265738461618137`*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {{},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[
0.012833333333333334`], AbsoluteThickness[1.6],
PointBox[{{0., 2.688094259}, {1., 6.651050631}, {2., 10.97636448}, {3.,
12.93577531}, {4., 14.60505538}, {5., 15.27245987}, {6., 14.0891266}, {
7., 12.7060682}, {8., 10.44592111}, {9.,
6.65659726}}]}, {}}, {}, {}, {{}, {}}}, {{{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
1.], LineBox[CompressedData["
1:eJwV13k81NsbB/CpUEnMzPc7Y48RyfJFuYXCfORWbqvKUq4oIVSWSHu6SJay
RJRUosVWhChEVLYo6ldIiy0hoZKmlH6nv+b1fn23c55znuc5I3DxWe82mcVi
TZvEYv35HfXTGBRdDTc/vyPSXiVUDK5Ryk6K592FoQMdL4NTxSBfoMMeSAkU
qm+Sse4qE8No9XzP5pQwYYz4qMzpr2LotC7iPUxJFPLPvXpX5y6OVp7oQ3nK
VWHXbe0BXysJbBDJtZSmFAn3VJTYFE6bhjjPzFUlKdXCY6LMm+zDMxASevqh
XmyLUEZ7xGXnVg72qkxXn/Zfp7C9ZO8Tw8s0VihKuZsqvxdOD7l/2bZIFsum
CZbrGH4UHlA//L9vvQpY/bVfKs73k5Dbl7Bbv0IZUs4NDz1DR4UO4ptz8+1V
oSDHEfZEfBO2X/2WaVokwP8eqd90+feH8KLharXgQTXM3le5u9zzlzD4cv/+
Em11xMQk8mXUWLi3T9U7b4kGCiN/PZMJn4R/V1+4LL90DmZP1jIJkpqCM3lh
VMhKTfSNlSaPBIrB0y7E6aj+XHzbTl1u+i2ONc27GCGthRuslJP//TcV0eui
Phl2aaGUf8L/r/5pMBmqi8uK0wY3E5PsfCQxSyLTqs1MBxW3HaP2t89AoPBc
VGSjDor/Nrs7w30m3liyfW/a6cJOtOaJd7806tO7JBLqdaHwOnwbPSyN6Z/r
lwgadPGmsurbnVFpOBqWHbzeqAu3qIWq4iwZzI+501PVpIuAWbP8zvFlcCVb
Jb7zhS7i//7IrbWUQcnes1uHunTRFBtlJ7ggg2OOlg+X/tDFFbpnyFOLjayN
vtc+aDBgxVon7tdlQ0zZWd5Gk4HDjLtmEfpszL3Q7VA6l4HM5MQTGQvYePTi
+qowHQb7h5drv7dgI0iw/MPMeQxW12e7uTqwcVPpJN2zmMHXIP9XTlFsCLLc
ZXOsGViPvw3xjmYjY+rMQ0PrGGQFrtI5EsfGd8E/efobGDjv0jhwPomN0RPB
x2/YMqhzaJF9lc6Gk+mPkfMODM4vWLxhUxkbItavXEtXBks/TK7f8JENE4H1
m6xABskLPPboj7Bh/MXLpXAvg5GgRtUZX9hIK9d8XbaPXKfP7K0SsVFgcbyo
/gCDITNGY54YBy2i16K2IwzOxNoflVbkgKMzUFodxmBgQc7COisOnC4eXtx2
moHwKKfr8koOsi+0xFcmMkioDzx5dA0HbknPujKSyHVnix4jGw7oPTddA86S
6+Ev4q46c+DTKSaaSGFg1s76GBLIgcoRyc1d6Qxij9qmm1/mYF1NrGpXHoO/
mIMjM65xcHyat/uVmwza2lLN2jI52KwrvL49n4Ga4WCLfy4HY4p3zfoLGBS+
C5bKLOGgOVnTq6OIQeuKvD1UM3m/bPHIhTIGKrwZVn2/ONCTvFfoWMPgQaXB
6VssLhouysuzahl4eNt1BU/hYobc9iNpxHk1lw4pT+diT97rlT11DJYcML65
nuZi5KjJ538bGLi/dVO4q83Ff6nZOxWbGdzIvDcYZ8+F4sv9cu1tDFa+XSCb
4sDF2+6ltz1fMuijsy2ubuaiY8t7nW/EakdPJ5Zs4+J0YVa09CsGiXZe6Pbh
4n7V/CX6bxgETaET/grnYtkBq63GXQyUTCLKzaO4qCreYV5IXOI90WcVzcV+
D9dhvW4Go219ZpsTuDg7+0WDag+ZT97d98dSuUhjle4Ze0f22+bti1tuc/Hc
NjLIs5/Bx1Ov3DpLuZhaI5b8lDiydl3sh3Iy32yb8EUDJD6Gpu9+P+Did2NA
h8QHBsaSnJi5zVx07S86EjfIQFBc0rW/n4uHO4X+3sMMygf1Z4YOcpHD4Tx6
QOyodsUoepiLJ1q3x+RHyPxOxpxI+8pFbUFY3T1iSVfXhY8mUQjW09jM+szg
M1s6UkmBAr9iu+nyUTIe3qHbn5UoyErOtg8jVlMY6K1VoXAoZe/OB8TrZldb
7tGgsDdpMMzsK1mvBUd+NRpQ0OwVPzh7jIHVokGdy4YU1uyM8nIg7jB3cDiw
kMIrkyybWGIZq4XFc0wpqGcXyH0n9nYY8g1aTsG/96tz+TcGEs6OF21XULh/
ban6R+IL2+obdVZT2Jen160gYtC486p263oK7DMN1v7EukFOPQZOFH5VGH2Q
/U7iF9LAnbqVQkDZulBzYsfwRRavt1FYyxKXcyWOiuNfiPCkYORkpJ9DPJD+
2L4rgMK93Y+UDH4wCM4wDbu9l8KWrNzINcQK17MKow9Q2KbFfNlB/E9RGGfx
UQrJv74UpBN3lIwKOSEUGtabSZUT76tw8X5/jNw/rOLcQpxRa/4oPopC5UTe
qMQ4AzTmfPeMpnDcz9lImbi1WWEu4igU6FwLmE88tX0sdDCRjC9w8ttNxBfe
uhZUnaWQIMGX2km8sOdp55kUCuveVBkeJnb/eMN8aRqFiW19u1OIJz4p7VK8
QmGQszUyizhxLPLcp2sUMNU/pZhYb1xUV5NFYcxyVuZ94oe/3UXnr1Pg1rrk
PSZ2FHs+JyCPfC8e+a3Eo9MsbVcUkPW/Wnyjk/jEzJshqkUUqlm1V/qJ1bkq
+WO3yfqk+SWNEJfxT3Y0lFI4GVUYMka8QXFcOr2cwtSyWK9x4gEVT7P9lRR6
DSVW/SYOVm/ZsfYBBe0v9NzJP0m8tZYma9RQePa9bGIKcT5TUDteR6Fj+cQT
MeIV8wXfmhso2L5uPffHnQtjNDKeUNAtWeXy5/7dD1bNK31K4aL0ivcs4snr
Jc0ePyfzkzEx+UW+d+ptjVVnK4XwUvmNImK1XcdsRttJ/KhBy8/E+T+WbJn6
lsJuubyfA8SW4aydCl0U+mvdgruI/8cr38u8o+CsMfP5n/i4ph8MQR+FVubq
aOOfeBmYxGz4QGF1l0FnJXFo+Viy+xCJ15LryYXE9KrCq/s/UTBfqyi4Snyl
zS//xCiFoqkHAxL/rO92/fKL3ygY7GxIOEZcPTpYl/+DQvYhmWB/4l4Zj85W
Fo3US173VxIHntf4+GEKDYWiYNmFxBI63aIJCRpJgSeMVYg1lzmxNWbSiBkP
6PxI9uPtZ4pKxmwaUvc2eD4jttraprmSohFloF5eTOxxyEboJ0/jOC/58SFi
0XTuylAlGlMjhZGbiSOSntglqZDrWS+4ZsRZ+Su872oQD7dFikg+fejD+enz
aURH/MXfSnwocCJDaQENGaftJsbEUmJlhfrGNBLeRGlJE+uqGDXYmtM4153p
W0Ty2duWGb+0gsbNXybJIpL/E10DErdW0zCZ/0O9mjjaN4Nba02jgptz6BRx
btRs7SE7Gk0tPac0iUcq5Tct2kZj1C/LwZLUl6NrW1xXu9MIkhCtmkrMfp3g
u8WTxozlRlL1pH7NE8mEh/nQUGuI/LSSeLeeRPHTgzTse4/ZW5B69zX5M70j
gUaLQ0Ps509k/SN9q3ySaJw2XlNzmbjwwJBPQDIN6XOvn9sSBzgM1B9OpWEw
pH6wkNTbz/JdR2NzaBg32vZuI/V5+Ezz4K0HNBy8IYwi9fxx+NrkkhoaNUnJ
Y9rEN/Y1Lq+opzG8dHJYHan/3hvr0mqbaPgZ/NCdQjwoW7mx/RUNI/MT5zxI
/+hPzHvA+kqj2dRkfJz0m9owvd3iIhpDN3cYRhBnBOaoSI7TqM7OtuQRe9hl
HKAm8XCZsRnUIv2ql5dqMGcmD5pey52Wkv7WnRCTslKDB9Pqogd/k374+pR3
QKItDxcl75xze072z7p3Gts38iCXOUy3/o/kE9uxxehfHiYZ5574h3hZ9D+L
2rbwQGvW7dV6RvpDhDpLaScPDis8TF83kXw/+jLqUjAP9lc2BNGkn8/xXpae
k8uD7PU3HyoryfmPuWtzOJ+HeZsMF9LELz8YSqy5xcOdr+wgt3sMYjwFnsMl
PLBZ26XEK8j+dP2pN7+ah5ebjvCMyXmj/t/8kuJXPLwROHTuIucRn39mPa2a
zkf2lGpnuSxynojr2dIkxcdvk5lSJpkkH9uyhl/L8GFyIjJiYwbp7x5GM7/z
+PAK1RTFXyX5FG5tpa/GR8b9gfSfaeT8Vx1cfm4RH7mHExecPsfgkeX7rIAd
fDS/DbbojiL1Ouq6SbA3H1rxJ+k3kQx8n/nXxvjx4eHy+HZLBAP+NlZv1l4+
Um0+B9QcZ+ASLC/oDOFj5Pj90jMhDL7fW5m0OoUPvnZUyOSDJH/N80LmNPJx
wDF3/LMHqS82PnxWEx+DBVFd7dsZVO3Qy2x7ykdKaXTlfXcGTmdzHp9o5eNG
+zSPWHI+TRjNkP/STcbv1GinvIXEN/tSbsUP8rxPUqm4Pek/svHtG+fKIqdC
P439N+n/n/YYRoXIgk/7dUsrMrjm9XW7XZgsavfPSyuTJ/21OyBFECGLH89c
tnnKkfx87i9+O1oWA25yHyt4DPxK/Fp6zspi6Pdj5a1sBhtDdh00z5OFr5TX
qL84Aw2e2/1Pr2QRdtw+sH1IF/eMbTZsXCiHVnF2481yXSiV0zWOIjkEHM+b
iLTSJf+Hppu4XJeHdCpyDtTpYG3o24kLfgrwvajY5MLoYM7D/qpXcxXR+3id
jVukNprOSn0LfqUIry8bN/i91EJqDidoy0UlmDYPJgdN14LikRXiQltlSPpF
j/y3YC5KT7q9tVGYhd3abj//NtOE/ryzWs/aZ6FBLkv52cI5WBPP2r0+VgVp
tktKl83TQJ28xnExbVXgQIeC9xR1WNRts9lyTBVbKOudh1vVkF51Y9PIG1Xs
Y7TLHdMFkG0fmKYEAQL8y4/8tFBF51MY1SYLIDqSa1D0WBkpc1fyN6cI4KE9
lGBTpYyNR2y/fjovwBllNZ+RW8po0vIqULwkwL52z2Oq55VRfjRe3+eaAMa3
9pm47FBGsl6vJr9QAIOO+IodU5WxITJKzrWRmL5nmWSmBOmO099Ej8n37lxz
7zFQQv2C1BcnmwQIvxb6W09dCUs6CxOKnwkg5zArvkxSCQbGb9gzXgoA5WaL
Oy2KkOo1mJ7fS55X+OuusZ8iahcv7lvWJ0CHma9ws6siQuKW1rT3C5CRccko
yF4RP0wdjol/FCDVtiei1EwRffEhkzZ9EcB3Zv8tnqQiLvef7Pg4St7n/eS9
xi8FOAvPVASPCWAtlx1jOKKAFwM5h3O+k/uXLLFY/kIBp1DsaDFOxnP/14rV
dQpYk1i5+MVPAVqjbtRblylAcvCRwo4JAaxSbO6ty1XAQ4sX33//FiB2dFjH
Ok0B/we/gERu
"]]}}, {}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0, 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotRange->{{0, 9.}, {0, 15.27245987}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.726573848835876*^9}]
}, Open ]],
Cell["\<\
Looks like the same good fit we got with FindFit, but now we can quantify \
that. (Check out the help page for NonlinearModelFit and go down a few pages \
into the Details and Options tab to find Properties.)\
\>", "Text",
CellChangeTimes->{{3.726573864161867*^9, 3.726573930110243*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"soln2", "[", "\"\\"", "]"}]], "Input",
CellChangeTimes->{{3.726573937930924*^9, 3.7265739494774837`*^9}}],
Cell[BoxData[
StyleBox[
TagBox[GridBox[{
{"\<\"\"\>", "\<\"Estimate\"\>", "\<\"Standard Error\"\>", "\<\"t\
\[Hyphen]Statistic\"\>", "\<\"P\[Hyphen]Value\"\>"},
{"A",
RowBox[{"-", "0.5075951015681828`"}], "0.013042951290943865`",
RowBox[{"-", "38.91719674830206`"}], "1.925545130712004`*^-9"},
{"B", "5.023244382731825`", "0.12193590016984932`", "41.19577889476971`",
"1.2949242343675163`*^-9"},
{"C", "2.564511982399986`", "0.2356410490185396`", "10.883129204700735`",
"0.000012214375143753844`"}
},
AutoDelete->False,
GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
GridBoxDividers->{
"ColumnsIndexed" -> {2 -> GrayLevel[0.7]},
"RowsIndexed" -> {2 -> GrayLevel[0.7]}},
GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
GridBoxSpacings->{
"ColumnsIndexed" -> {2 -> 1}, "RowsIndexed" -> {2 -> 0.75}}],
"Grid"], "DialogStyle",
StripOnInput->False]], "Output",
CellChangeTimes->{3.72657395032309*^9}]
}, Open ]],
Cell["\<\
Note that now we get an idea of the errors (at least the statistical errors) \
in our fit parameters.\
\>", "Text",
CellChangeTimes->{{3.726573957525022*^9, 3.726573983412018*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"soln2", "[", "\"\\"", "]"}]], "Input",
CellChangeTimes->{{3.7265739984176598`*^9, 3.726574008043118*^9}}],
Cell[BoxData["0.9995160449817018`"], "Output",
CellChangeTimes->{3.726574008820772*^9}]
}, Open ]],
Cell[TextData[{
"The value of the square of the correlation coefficient (that\
\[CloseCurlyQuote]s what ",
Cell[BoxData[
FormBox[
SuperscriptBox["r", "2"], TraditionalForm]],
FormatType->"TraditionalForm"],
"is called) gives another quantitative idea of the \
\[OpenCurlyDoubleQuote]goodness of fit\[CloseCurlyDoubleQuote]. The closer \
this value gets to 1, the better the fit. In the physical sciences, we \
typically want this value to be bigger than 0.9 for any firm conclusion."
}], "Text",
CellChangeTimes->{{3.7265740334029083`*^9, 3.726574180836123*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"soln2", "[", "\"\\"", "]"}], ",",
RowBox[{"Filling", "\[Rule]", "Axis"}]}], "]"}]], "Input",
CellChangeTimes->{{3.7265741879480677`*^9, 3.726574211940154*^9}}],
Cell[BoxData[
GraphicsBox[{{}, GraphicsComplexBox[CompressedData["
1:eJxTTMoPSmViYGCQA2IQDQEf7B/si+kRW7zfHirg0NDoV9hTdns/hMvhcMCe
2fJ+8E2ovIBDA19P36TlB6DyIg4MDBuObl+/ESov4ZDQcHPvKdZLUL6Mw4Eb
5a+/3LwKVa/g0CBYH+wgdBjKV3JgYNReusPmOFS9igPDg1pNAf/4/TD3MaAC
B1QuBxpfAI0vgsaXQOPLoPEV0PhKaHwVNP7gDT8AyhFeaQ==
"], {{{}, {}, {},
{RGBColor[0.368417, 0.506779, 0.709798], Opacity[0.3], LineBox[{12, 2}],
LineBox[{14, 4}], LineBox[{17, 7}], LineBox[{18, 8}],
LineBox[{20, 10}]},
{RGBColor[0.368417, 0.506779, 0.709798], Opacity[0.3], LineBox[{11, 1}],
LineBox[{13, 3}], LineBox[{15, 5}], LineBox[{16, 6}],
LineBox[{19, 9}]}}, {{},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[
0.012833333333333334`], AbsoluteThickness[1.6],
PointBox[{21, 22, 23, 24, 25, 26, 27, 28, 29,
30}]}, {}}}], {}, {}, {{}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{0., 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImagePadding->All,
Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
(Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
Part[#, 1]],
(Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
Part[#, 2]]}& )}},
PlotRange->{{0., 10.}, {-0.4291106325636278, 0.3957441384090963}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.726574212951767*^9}]
}, Open ]],
Cell["\<\
The fit residuals are just the amounts by which the curve fit misses the data \
value at each point. You want theses residuals to be small compared to the \
y-values in the graph, and you want them to be ideally scattered about zero \
(some above, some below). If these values are large, and especially if they \
show a smooth trend (smoothly going up or down, or making some smooth curve), \
this means that your fit function is probably not the correct one. Choose \
another one and see what the residuals look like.\
\>", "Text",
CellChangeTimes->{{3.7265742211071987`*^9, 3.726574382208156*^9}}]
}, Open ]]
}, Open ]]
}, Open ]]
},
WindowToolbars->"EditBar",
WindowSize->{1323, 867},
WindowMargins->{{4, Automatic}, {Automatic, 4}},
FrontEndVersion->"11.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (July 28, \
2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 111, 1, 92, "Title"],
Cell[CellGroupData[{
Cell[715, 27, 886, 19, 140, "Section"],
Cell[CellGroupData[{
Cell[1626, 50, 432, 8, 32, "Input"],
Cell[2061, 60, 898, 26, 32, "Output"]
}, Open ]],
Cell[2974, 89, 945, 17, 68, "Text"],
Cell[CellGroupData[{
Cell[3944, 110, 1251, 25, 54, "Input"],
Cell[5198, 137, 2257, 52, 277, "Output"]
}, Open ]],
Cell[7470, 192, 728, 13, 68, "Text"],
Cell[CellGroupData[{
Cell[8223, 209, 238, 5, 32, "Input"],
Cell[8464, 216, 839, 24, 32, "Output"]
}, Open ]],
Cell[9318, 243, 383, 8, 49, "Text"],
Cell[CellGroupData[{
Cell[9726, 255, 163, 3, 32, "Input"],
Cell[9892, 260, 1580, 37, 245, "Output"]
}, Open ]],
Cell[11487, 300, 367, 8, 30, "Text"],
Cell[CellGroupData[{
Cell[11879, 312, 320, 7, 32, "Input"],
Cell[12202, 321, 175, 5, 32, "Output"]
}, Open ]],
Cell[12392, 329, 519, 10, 49, "Text"],
Cell[CellGroupData[{
Cell[12936, 343, 344, 7, 32, "Input"],
Cell[13283, 352, 335, 8, 32, "Output"]
}, Open ]],
Cell[13633, 363, 654, 17, 49, "Text"],
Cell[CellGroupData[{
Cell[14312, 384, 200, 4, 32, "Input"],
Cell[14515, 390, 267, 5, 32, "Output"]
}, Open ]],
Cell[14797, 398, 303, 5, 30, "Text"],
Cell[CellGroupData[{
Cell[15125, 407, 236, 6, 32, "Input"],
Cell[15364, 415, 6867, 126, 239, "Output"]
}, Open ]],
Cell[22246, 544, 205, 4, 30, "Text"],
Cell[CellGroupData[{
Cell[22476, 552, 160, 3, 32, "Input"],
Cell[22639, 557, 7059, 128, 245, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[29735, 690, 232, 6, 83, "Subsubsection"],
Cell[29970, 698, 274, 7, 68, "Text"],
Cell[CellGroupData[{
Cell[30269, 709, 280, 6, 32, "Input"],
Cell[30552, 717, 1081, 27, 44, "Output"]
}, Open ]],
Cell[31648, 747, 237, 4, 30, "Text"],
Cell[CellGroupData[{
Cell[31910, 755, 272, 7, 32, "Input"],
Cell[32185, 764, 6817, 125, 239, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[39039, 894, 160, 3, 32, "Input"],
Cell[39202, 899, 7037, 128, 245, "Output"]
}, Open ]],
Cell[46254, 1030, 299, 5, 49, "Text"],
Cell[CellGroupData[{
Cell[46578, 1039, 147, 2, 32, "Input"],
Cell[46728, 1043, 1038, 23, 81, "Output"]
}, Open ]],
Cell[47781, 1069, 191, 4, 30, "Text"],
Cell[CellGroupData[{
Cell[47997, 1077, 141, 2, 32, "Input"],
Cell[48141, 1081, 88, 1, 32, "Output"]
}, Open ]],
Cell[48244, 1085, 575, 12, 50, "Text"],
Cell[CellGroupData[{
Cell[48844, 1101, 243, 5, 32, "Input"],
Cell[49090, 1108, 2067, 46, 229, "Output"]
}, Open ]],
Cell[51172, 1157, 613, 9, 68, "Text"]
}, Open ]]
}, Open ]]
}, Open ]]
}
]
*)