(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 21104, 638] NotebookOptionsPosition[ 19101, 564] NotebookOutlinePosition[ 19529, 581] CellTagsIndexPosition[ 19486, 578] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Fitting an exponential curve", "Subtitle"], Cell["\<\ The following data were taken March 26, 2008 in a PS210 lab experiment \ demonstrating \"Newton's Law of Cooling.\" These data were collected using \ the glass thermometer.\ \>", "SectionFirst"], Cell[CellGroupData[{ Cell[TextData[{ "The {x,y} coordinates represent {time(sec), temperature (", Cell[BoxData[ FormBox[ RowBox[{ SuperscriptBox["\[InvisiblePrefixScriptBase]", "o"], "C"}], TraditionalForm]]], ")} ." }], "Subsection"], Cell[BoxData[ RowBox[{ RowBox[{"data0", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "84"}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", "54"}], "}"}], ",", RowBox[{"{", RowBox[{"60", ",", "47"}], "}"}], ",", RowBox[{"{", RowBox[{"90", ",", "42"}], "}"}], ",", RowBox[{"{", RowBox[{"120", ",", "39"}], "}"}], ",", RowBox[{"{", RowBox[{"150", ",", "36"}], "}"}], ",", RowBox[{"{", RowBox[{"180", ",", "34"}], "}"}], ",", RowBox[{"{", RowBox[{"210", ",", "32"}], "}"}], ",", RowBox[{"{", RowBox[{"240", ",", "30"}], "}"}]}], "}"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"data", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"30", ",", "54"}], "}"}], ",", RowBox[{"{", RowBox[{"60", ",", "47"}], "}"}], ",", RowBox[{"{", RowBox[{"90", ",", "42"}], "}"}], ",", RowBox[{"{", RowBox[{"120", ",", "39"}], "}"}], ",", RowBox[{"{", RowBox[{"150", ",", "36"}], "}"}], ",", RowBox[{"{", RowBox[{"180", ",", "34"}], "}"}], ",", RowBox[{"{", RowBox[{"210", ",", "32"}], "}"}], ",", RowBox[{"{", RowBox[{"240", ",", "30"}], "}"}]}], "}"}]}], ";"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"data", "=", "data0"}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"data", ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"PointSize", "[", "0.020", "]"}]}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "100"}], "}"}]}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointSize[0.02], PointBox[{{0., 84.}, {30., 54.}, {60., 47.}, {90., 42.}, {120., 39.}, { 150., 36.}, {180., 34.}, {210., 32.}, {240., 30.}}]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{0, 100}, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.441884903709178*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f1", "=", RowBox[{"FindFit", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"a", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "tau"}]]}], " ", "+", "b"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "60"}], "}"}], ",", "tau", ",", RowBox[{"{", RowBox[{"b", ",", "25"}], "}"}]}], "}"}], ",", "t"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"a", "\[Rule]", "49.61115309662235`"}], ",", RowBox[{"tau", "\[Rule]", "45.011995376002616`"}], ",", RowBox[{"$Failed", "\[Rule]", "32.919531547473404`"}]}], "}"}]], "Output", CellChangeTimes->{3.441884903848166*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "FindFit assumes that each data point has the same error bar which, in this \ case, it does (~1 ", Cell[BoxData[ FormBox[ SuperscriptBox["C", "o"], TraditionalForm]]], ")." }], "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], "/", "tau"}]]}], " ", "+", "b"}], "/.", "f1"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "250"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "100"}], "}"}]}], ",", RowBox[{"Epilog", "\[Rule]", RowBox[{"Prepend", "[", RowBox[{ RowBox[{"Point", "/@", "data"}], ",", RowBox[{"PointSize", "[", "0.015", "]"}]}], "]"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.4418849107083607`*^9, 3.441884913466683*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwd0X8w1Hkcx/Ftb9MPFRUnkkJNolLUVGT3tT/stmspbJLBxrZrc/Jtr6Y2 +nV+xE4ov6Zfm3YJlzo5yklxYu8STY6rK4mJ7kpXulAi5D7bH+95z2Oef73n HUBF2stpNBqdjJ9xMx8p+OKxgTeQ6IcKHcu9WBV6nvvBwu8h0Z70fKDNYTXo rRV9pmpITjuv8LzjwvpTz6AV1cZDkhZly7lyn/VCN+KbeHockv3js+RF7axB 3eDgAKMWEvEZqWa8Dt/o+gXLPO0QOJJc9mOWLebq/ntv+MsPgSGMAIkJDYv3 eZsmZgIBtgXsT45PsVKWWRY0OBv+FpsbF0ykY435YX/1LSa2LKdZvu5ahA3B z36q3auDX27llMMl/mDyR9s6hV3wtcssPZLmDXZLZF+cjxw+FQ/HRrJqwYtt YXyqZkN0wLMtgZEGfvd8WneVKYSUzF6vWQJh6sZYQb8LNlGvTs7q64WY7lVn c9MDgvPWY8WP2+BH/b5EKrwG/nMTwcymcmzRzP5iPr0GfOe8n6v41+HfUNZR VGUD72x1CoO6hsB51arRRAd4T/IqP7vuNiRRKVnT3MLB+yH71qVid2y9aLjk u2w2eFYuY5+0PQgaiKE2e1HglllG16fbIHhl0vKUszxwQ9Y0uoWkY3uySf6R OQZwp7i5rvYSI+SROH/FngPg1N4psNuwBKFm9Jy5lnpwEnYu3Vx8F2EyR2bl E9IFHzWD8hMI11V0nOoxgGMlejCpeh2kT6bSHVt5YL9NHK7lGLAj58We6tgq sBtm1OlaCxDhFtd+ZJ4GbN2xy/HWDojoFOTpWleBfdSK89F+OiKThpwujr4H O6LP7JeHxZDZh16pFsWALcw7K6w/BVnT8An34hqw3SwchhL9sHOXWNbo1AS2 ndr1xNZcyBkJkzbM7wd7pkfJwMrpkF99OzptQR8wlmBRQe+FQmTW94FbCLzL 5l543QzFB6UXU+0GdFu1uafvRFS2Orngj3jgUafS2eMulK7RWf+c6wSaJ0T5 SZ1QPskcO77PEqinot/t7cWu/SopRzUNqFYv3JEnR7T9cbuwQ21AuSGRnmqN 6HuH8iWfbwBXTZvke47hO4XLwihtIFCkao37Mgcx5rOG28V3Af3kuBvBKsSU ljQrFg0B2h9sMnj3sFvkFNkSxQLOpOaOxOdi90BH0Onb4UDONQYz7zpiT/G5 8WvXA5kH20Mvs0B5vEyrud0LZEgTFl+IIP48IB1RGu1UcU4JytMsmL7W2MNb uB2Xidc78ktfGW2rdL5PrLnZontJHFr56+hMUBuXv+xO+pt4+5ur0qOgvPay J4u7iSVBaqsJ4pSiYs0O4sBx12wNsXaG8LfnRhf0BtwiNjzOYBod0B+mnQyK +W2s1Sqj/dO2ba8kvnl+mUUXsV99yuEk4ge0Zv9wo3cFrqaIexS7079280VO VaBYpqvLTcg/MnzDarq3EYc2egw/I/YZmRpH7mepVnSuDTX6YtemAuLk7KOq r51/IwTEpeEN/3YQi3KiBNbE9QbF0hCjPTkOPsSPnafKnhILe+xQSPzmZInO 2IWpo/mOxBMfxc++dten/Lr/AX0WKKU= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, Epilog->{ PointSize[0.015], PointBox[{0, 84}], PointBox[{30, 54}], PointBox[{60, 47}], PointBox[{90, 42}], PointBox[{120, 39}], PointBox[{150, 36}], PointBox[{180, 34}], PointBox[{210, 32}], PointBox[{240, 30}]}, PlotRange->{{0, 250}, {0, 100}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{{3.441884903967348*^9, 3.441884914810042*^9}}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Using the fitted parameters from FindFit, calculate the ", Cell[BoxData[ FormBox[ SuperscriptBox["\[Chi]", "2"], TraditionalForm]]], "." }], "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"chisqmin", "=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], RowBox[{"Length", "[", "data", "]"}]], FractionBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"a", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}], "/", "tau"}]]}], " ", "+", "b"}], "/.", "f1"}], ")"}]}], ")"}], "2"], SuperscriptBox["1.0", "2"]]}]}]], "Input"], Cell[BoxData["48.625622900913214`"], "Output", CellChangeTimes->{3.441884904109703*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"chisq", "[", "\[Tau]_", "]"}], ":=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], RowBox[{"Length", "[", "data", "]"}]], FractionBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"a", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}], "/", "\[Tau]"}]]}], " ", "+", "b"}], "/.", "f1"}], ")"}]}], ")"}], "2"], SuperscriptBox["1.0", "2"]]}]}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Plot the ", Cell[BoxData[ FormBox[ SuperscriptBox["\[Chi]", "2"], TraditionalForm]]], " as a function \[Tau] to determine the uncertainty in the time constant (\ \[Delta]\[Tau])." }], "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"chisq", "[", "\[Tau]", "]"}], ",", RowBox[{"{", RowBox[{"\[Tau]", ",", RowBox[{ RowBox[{"(", RowBox[{"tau", "/.", "f1"}], ")"}], "-", "2."}], ",", RowBox[{ RowBox[{"(", RowBox[{"tau", "/.", "f1"}], ")"}], "+", "2."}]}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV13c8VX0cB3AaKkWTys7ISIRQSn7h/u659eBeI7kNCqE0jJDITBpSyYrK SCKENiqlMkpReITEI0WlISkNPZ/zl9f7D+f+zvf3Xcd25+Z57iIiImNERERn sn8JJ+LgdD1FU8JT3jYhsucDXGJJSrbCGp4LC2bCLwL5FimEN6/3eVKSGuFE ar7c6TEEu9Rl+Q/CdXmFxQWEp+S49OyrJMKJnmKuZS5HeApUVj5EhXBi4v1U 3H4TnsyC99FBHoRzOLFZKr+c8GbMvjXBvohwkvaK6vkeIzzR+UViRaKEk5Hr 2GlnRJhB8WyfngzCyc9IsStsI0zP643ySU2Ec9V3pYhlM2GamW6/CQ6Ec+fw LNOhTMLUrL07r3EC4dSF2JR23SdM+ayKMQn7CafZxP8z7SVMUeNAMA/nf2V2 //aDFsKcmyZnU+VCOO/kQ+5HKBImVcXzrTghnCFz53UOSwgTnzV5wuojhDO6 eEFZtg1hYoZP9Ab7EDpx1EWw8BZh9nFHKr0SCJ059U5Qw2/CBLT3+E9aQqj8 SMy/xScIszPOcUlzJaEaxcZ/xecTxkvx3c3K+YTqF828HnKIMK6PgtV/xxBq Gu2Yu4ASZmNF7fiAt3BB/SYdN/iiXXnHINxsIasXCaf2vdshT+iK+Trxxnfg QP1LN/bBNeMDLXHejfoWraF4Hpl8mbtBkzAb8tfHM0WEmidIvo2fQpj1WVu+ fdxIKDfzofrgM8IIS9pq/FLh+tORz/6DMx9U1unAI/4vS7/Axyv++N0ilBGo JPhNg30aLKwjCOWJhI9+tYb1D46KXyN0lfOy5q+PCON0bZV0x0lCreRKooYQ /7UPf08dkIaZ2M6mJvjGznVNnbC/y9IrPXDex76YdPix5KfdY+C4GVX/qhJq vddb+I3AdgNhsxoItWlT1/uGeDv2fPzEMyPUNjn91fA1WELp84t38L3c7h9F sMjakgft8EBpz88cwqz5mr6iN4NQO4uat6Mn4TZ6y+EV/Gnos5gvfEGTkZAi 1IGxEpHWgalqaCB+z/HHHyXD84RxODRYaN9H6FrVicrGZ+Cw7MpRXdhmpurS RNjf5XT0BzhHU8M0Gt7wp7hShFAngcMizmZY74Gp6FlChXmFxEGBMPbtB5uF zwjdINzg4p8MG/etukzgyP1TygZhnbJ2zWI4v+Cm6D+wWnL+mvHwz9/TjyI+ 9jN3jduD821MO3Uvdwdh7D6l5cQjns4draptlnBeZ+yEY4RudrbtN/sMz59h VIH82XwgODFmNayg7S82Br6UtfIx3s9Omj/Ab4H/fjnltA4Wu5GQKCTUNePY at+HhLHtG1MQcpFQt+76wpx0uCDid7s9oR5ujO9kHmyqHX4qH45vvKOTBBv+ 7bixAi5fN0WAfLLVeWUa8oRQzxnbz6cEw0r3phypgStPtKvlw+Nnv3xwg9Ct ci8tCfJT0DBucWwhzGw57roWrk1+smgx7PepM+YcXGW89WE0/GhMUP1y+Nqp AkEkodv2qBcKt8Onv1h1/ibUu9lX2v8JvKNzsEaf0O0iv1xTZGAv/aKgZFg7 urh8C+x6Ypsz3n97VPJq0VF47e4+t1FCd+hVhMXrwuaPPzBWhO48KtaXdwye Wz7WYx6hvjStrNOWMPyGpKKr5+GdK0sp7of/uLh9+VM4+W1eYQlc07y3divc b5ASivrh310xJ5VLqN/hRwGy9XDpkdtZpwj1b/il58TAyYoXnqJ/+I9kalT+ gBO6JxsUELpbmaukkQsfK+haMhn2S5D8IQYfcu2uXkNogLTWhyTUG3+fDTl6 l9DAdcLcJjPYU9bJOBeOFjm77CPsHtWugnoJLDiflI185m8efrLqEfxncL8v +hd/w+8Jn6IIDco45Dq9DLaPeqiCeg6qXSTcsxW2NXvqnwV/aRF0z4X5Yjwr 5OMeC2VSHASvLj1qhv/f01umYIP+zbeIFV7SJjRY0kXq2ht4ZXJvw3PYWExC AfXDNyvpqjGCYwW/Bobg5ePVOssI3avV/+IIGz9DjXcm22C7+IYhF9jg8JeU ITjEsGY9G0+9EU+ZErg+4voC1ANf52uWzhVCQ3bNPVmnDWu6Mxq2cPJXAzs2 PhrKqr864DuPG18i3/nz+yduSSQ0dGq45Gf0Y77qyf8Kz8HGToXB82CVPS+c v8POBv+MvQEre7Y4xMJFbw5Kv4WV3P67Jgm3VKpn7IMVfb58nwiPpj7UkoYV Dk5QNyB0n7XVWFPUJ1+u2UJ0DhwwP/NhB2txr/D78BkRwveDZXmJaZnwwOVQ 1yx4bu+odQWhYVJxsgMmrC1MbuH3wky3lAU0wnPyQ3zjYHfiJOLJWrZKiPwJ i5v7/ZAIPDt5qvM4+MrXpJnJrOVcvMfCHfWLT6Pf8aUvXvPWJDR83Pnn8x+y 5kxfdRzWDvcp3gBL9e0alIHthdOWsvcjldhs9QIOMSiqOsLaagWD/Ak/J2Fl pcJa8uIT9JfwR2/et7D5M6tNvsYb/lp5yIWtl1mXTqpMIDRC9pRGfz/reMka 9KcIC79qv3DWQXE5iE/ENiv3P3NYe0/LCoET1MfFFLP2Si4wgctFsqdxWfso F2JeRfzXtjK1k3VkyTHMv0jxy10qAazPWK5E/kXqx+0rlGBd1ZaP+RIp9JA3 Osf6q/+1cDiSlFcuZ8+/cIZgOZwvI+Q9Z73rsrsT3Pj1x3M2/6VuO77sgkfq kzew9S4t9TcH+Rw1L9fwTSrrwPwLy2BeeNOuRax71jZgXkf5CH1Hatj4O4lP RP+KSl08PcqZdevt5egPUX1vrZOOsvc3rGuOfh49rXJASY2974T373fBS04d zquAZZZdXJYNx1rVVLxn8yVPb98fuFjdg0ay+bTjB4P7j24VHd+AfseXX3HX XZHQ/epXzHvY/qLwW0hNYZu4bu8uNl/faXb5wYEeYcOBbD53/wo+DFfL3Jp0 ns33T+eEeP+YLRHGi8bC6lHi7H4Rc7wiIIm9f42ONS6YXzEV36/+mQVrrfS8 in3nwMwdi+uw7/AX6nRZnYEr1+m7sb9vkDzHEvvAgfdJPnVs/hpan3yC+oqV bizRE8LG01uG8b6x2xnd0W2wSfnFYmVCD8oYLkxh42OentCNfecgZ5f339mw ZamJK/r9wV0XC9wzYPo8p3AGXK28QJ/tH6uX1qmj3x3yn6rxuIntZwnc576E Hq7vVxVl4+MVPJqGeXt4RNXVIxv2/uf1n2FCj6g6Zz1h+8dOraPBjXBwy7w0 U3j3JC0nzIs49SpFQ7YfRdhF3sG8OBpxWsaLjVe63Jk65PfRglanBswP/tn8 fD7r1pmpRmz/yLYMtMTz43UOzRlbC+enOUzfDLcHSp3+AN9skediXztuKJj2 TB9uC6upQz0e33DhZyj2B35nuFqlBBwj8lpDHP7v0G1x5N/xluLroZi3/He3 M3iyhJ4ImrZRE/sn/1em5qSXhCaUPb0Yjv1CoChYclYL7lFPXIB6Eah8IxXY 705ODtvXUgdrZEc24Twn1+vaaqP/CPQVrqThfk/+OfqjFfUr4MQ/1lAiNIlY U91Sdr6uv8HBvpPkeV63DfNE4HfOoQ31lXRsdO5+9FNB0LC2pBfcXfihHfuQ IOpGxB/sS8lREgkHZsCp1xtT0R9SHjzu6sQ+JqjO+PCoidC01av3xulhP1gg t9ePdeB4MSlJWH+J5qQUOOv2sXTEy3apc7NkFfxTL6cA8bflPlcfmUVoet7s p4/lYHfeNR72tzNir1UnY17aZgWsv419LONeyONY9F87tfpp0tiXMj4aOU5D v7LT8VxjiP0sc+6n7hT0Szvjqcm2r2GfzcMX8P52TKT4HGtCs5R4StXV8Lbh onrso9n7pHaLIR52VxU6arFvnl9aJB8di32Nv1XcMx52S7PRsYKFDuo1qIfz 8bHhrXh/e7fV9V2tcO/m1wtOw0HrLjdi/uQen53/DM+zzxzf04H+faEvYonS S3go6lbDDkIvJtvb3kJ+OJxTzZ4fTGjx0M8D3/A+jmrfoqo/E1oiP8719ThY r9do6gDMkTR7thA27Xy/RwxOnPe9CPPK0eGzs1MPoaWGPA8PJTjmsNckzJfL u1Noqzvcn/mlKIfQq8NG48uwP6y9zrcsxnnLvvtEhYti/99VuOYH9qFy5Zci q1FPwnCHSjvcV7k1Eyq9GD42qf858rc8RzGgAPNJeCkx5Dnmd4VtvVdrA/xp 6dl3mC+3ijT5eqjfdbuu3vvuRmilW7dcD/aX9SEhd8di/j14xr9GQ/F9c93s rrIloU+HhMFqyFfXjYLRtNOEvuixFguTIoznEe6PDOxrr/Pz6mZdJ8yOjqnC iT8J/XhzkUKVFmF2u2/QrcL3z0iGQn4h5lmoblHbtiLCHfvEuLoknjD7FUZk 21UJV6LWpFBXnTBx0Slvbc0Jd056Tm0/5nkyf1jDI4dwla8vCyzH92Om7JES zDvuwvyl6+sGCFOwretMNSHcJVvnidtjPtwwXTS7zoJwzX0MmG99hKnqtRHY NBGuFb/wdhP6e+NS2yG9dMJ1DLR98MiAMJ180VOdRoS7aWt34vIVhHl/Z9My bTHC9TZqTevH9+3Pv5mFSVGEG8D7aXCpifAmdSSUKNYSboRh9JhSZ3xfX3g0 5+UV2OiZ3pcsOL9VeuMDeInSpkVv4IJe6ZDFsEnFnUvecPFfhSO6sNm3kII9 8E0Da5Hd8Kotw+cT4Np0LZMdsDPv3alq+N3OfwctYZeUuW3i8IdQjdFieNMb 3lxr+OOhkNwg2C0yNxnnkxnMWfBEC/YqdzvZDf9sy9sZAftrdx79RXiyk6ny mEB4916JJytgCbv0pjVwQK2pRCQ8dZOMeQy8xz3tyER4ZojckCi878zag1Kw bOnK/Ptw2IcDNWth+TuNztZwhMn1CemwYr1rsQYc9a9UjDKs0ndsfgO8X83y gTusNqz1TwEc4+c/7gKsMbZasgU+cDfb4gOsNX3Ldhv44NRnUbqwtuJET6X/ AYKFR0w= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{43., 49.}, PlotRange->{{43.011995376002616`, 47.011995376002616`}, { 48.625622939166036`, 50.26889127627496}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{3.441884904473227*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Find the uncertainty in the time constant (\[Delta]\[Tau]).\ \>", "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f2", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"chisq", "[", "\[Tau]", "]"}], "\[Equal]", RowBox[{"chisqmin", "+", "1"}]}], ",", RowBox[{"{", RowBox[{"\[Tau]", ",", RowBox[{ RowBox[{"(", RowBox[{"tau", "/.", "f1"}], ")"}], "+", "1"}]}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"\[Tau]", "\[Rule]", "46.596907373262994`"}], "}"}]], "Output", CellChangeTimes->{3.441884904616152*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[Delta]\[Tau]", "=", RowBox[{"Abs", "[", RowBox[{ RowBox[{"(", RowBox[{"tau", "/.", "f1"}], ")"}], " ", "-", RowBox[{"(", " ", RowBox[{"\[Tau]", "/.", "f2"}], ")"}]}], "]"}]}]], "Input"], Cell[BoxData["1.5849119972603773`"], "Output", CellChangeTimes->{3.441884904705284*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Print", "[", RowBox[{"\"\<\[Tau] = \>\"", ",", RowBox[{"(", RowBox[{"tau", "/.", "f1"}], ")"}], ",", "\"\< \[PlusMinus] \>\"", ",", "\[Delta]\[Tau]"}], "]"}], "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{ "\"\<\!\(\*SuperscriptBox[\(\[Chi]\), \(2\)]\)/dof = \>\"", ",", "chisqmin", ",", "\"\< / \>\"", ",", RowBox[{ RowBox[{"Length", "[", "data", "]"}], "-", "3"}]}], "]"}]}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"\[Tau] = \"\>", "\[InvisibleSpace]", "45.011995376002616`", "\[InvisibleSpace]", "\<\" \[PlusMinus] \"\>", "\[InvisibleSpace]", "1.5849119972603773`"}], SequenceForm[ "\[Tau] = ", 45.011995376002616`, " \[PlusMinus] ", 1.5849119972603773`], Editable->False]], "Print", CellChangeTimes->{3.441884904806321*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"\\!\\(\\*SuperscriptBox[\\(\[Chi]\\), \\(2\\)]\\)/dof = \"\>", "\[InvisibleSpace]", "48.625622900913214`", "\[InvisibleSpace]", "\<\" / \"\>", "\[InvisibleSpace]", "6"}], SequenceForm[ "\!\(\*SuperscriptBox[\(\[Chi]\), \(2\)]\)/dof = ", 48.625622900913214`, " / ", 6], Editable->False]], "Print", CellChangeTimes->{3.4418849048116417`*^9}] }, Open ]] }, Open ]] }, Open ]], Cell[TextData[{ "So, the final result is \[Tau] = (49.0 \[PlusMinus] 1.3) sec. with a ", Cell[BoxData[ FormBox[ SuperscriptBox["\[Chi]", "2"], TraditionalForm]]], "/dof = 23/6.\nThe ", Cell[BoxData[ FormBox[ SuperscriptBox["\[Chi]", "2"], TraditionalForm]]], "/dof is a little high and probably reflects a situation where the actual \ \"setting error,\" or uncertainty in the temperature (\[Delta]T), is probably \ underestimated, and is more likely somewhere between 1 ", Cell[BoxData[ FormBox[ SuperscriptBox["C", "o"], TraditionalForm]]], " and 2 ", Cell[BoxData[ FormBox[ SuperscriptBox["C", "o"], TraditionalForm]]], "." }], "Subsection"], Cell[TextData[{ "Question: What happens if you increase the setting error in the ", Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ SuperscriptBox["\[Chi]", "2"], "'"}], "s"}], TraditionalForm]]], " from 1", Cell[BoxData[ FormBox[ SuperscriptBox["C", "o"], TraditionalForm]]], " to 2", Cell[BoxData[ FormBox[ SuperscriptBox["C", "o"], TraditionalForm]]], "? More specifically, how does the ", Cell[BoxData[ FormBox[ SuperscriptBox["\[Chi]", "2"], TraditionalForm]]], "/dof change? And, how does the uncertainty in the time constant (\[Delta]\ \[Tau]) change? Does this make sense?" }], "Subsection"] }, Open ]] }, WindowSize->{785, 678}, WindowMargins->{{4, Automatic}, {Automatic, 4}}, Magnification->1.25, FrontEndVersion->"6.0 for Mac OS X PowerPC (32-bit) (May 21, 2008)", StyleDefinitions->FrontEnd`FileName[{"Creative"}, "PastelColor.nb", CharacterEncoding -> "UTF-8"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 48, 0, 38, "Subtitle"], Cell[641, 25, 205, 4, 42, "SectionFirst"], Cell[CellGroupData[{ Cell[871, 33, 229, 8, 45, "Subsection"], Cell[1103, 43, 689, 22, 66, "Input"], Cell[1795, 67, 624, 20, 66, "Input"], Cell[2422, 89, 74, 2, 46, "Input"], Cell[CellGroupData[{ Cell[2521, 95, 368, 10, 66, "Input"], Cell[2892, 107, 403, 10, 323, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3332, 122, 458, 14, 54, "Input"], Cell[3793, 138, 276, 6, 46, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[4118, 150, 214, 7, 64, "Subsection"], Cell[CellGroupData[{ Cell[4357, 161, 678, 20, 78, "Input"], Cell[5038, 183, 2272, 49, 322, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[7359, 238, 179, 6, 47, "Subsection"], Cell[CellGroupData[{ Cell[7563, 248, 763, 24, 84, "Input"], Cell[8329, 274, 88, 1, 46, "Output"] }, Open ]], Cell[8432, 278, 798, 25, 84, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[9267, 308, 222, 7, 66, "Subsection"], Cell[CellGroupData[{ Cell[9514, 319, 357, 12, 46, "Input"], Cell[9874, 333, 5519, 97, 298, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[15442, 436, 89, 2, 42, "Subsection"], Cell[CellGroupData[{ Cell[15556, 442, 362, 12, 46, "Input"], Cell[15921, 456, 144, 3, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16102, 464, 239, 7, 46, "Input"], Cell[16344, 473, 88, 1, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16469, 479, 453, 11, 73, "Input"], Cell[CellGroupData[{ Cell[16947, 494, 373, 8, 25, "Print"], Cell[17323, 504, 407, 9, 30, "Print"] }, Open ]] }, Open ]] }, Open ]], Cell[17769, 518, 675, 20, 110, "Subsection"], Cell[18447, 540, 638, 21, 78, "Subsection"] }, Open ]] } ] *) (* End of internal cache information *)