(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 48099, 1366]*) (*NotebookOutlinePosition[ 50268, 1430]*) (* CellTagsIndexPosition[ 49835, 1416]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Mathematica", FontSlant->"Italic"], " Tutorial" }], "Title"], Cell[CellGroupData[{ Cell["Algebraic Manipulations", "Section"], Cell[BoxData[ \(\(Expression1 = EnzS[t]\ \((1 + \(k1b + k2f\)\/\(\(k1f\)\(\ \)\))\) \[Equal] S[t];\)\)], "Input"], Cell[TextData[{ "In ", StyleBox["Mathematica", FontSlant->"Italic"], ", every object (from equations, to numerical solutions and graphics) is \ represented as an expression. All expressions have a 'Head', that tells you \ what type of expression the object is. In ", StyleBox["Mathematica", FontSlant->"Italic"], ", one can find out information of commands by typing ", StyleBox["?Command", FontWeight->"Bold"], ". So, to find out what head does, type and evaluate:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?Head\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"Head[expr] gives the head of expr.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"Head", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384087771-5332558"] }, Open ]], Cell[TextData[{ "Now find out what head ", StyleBox["Expression1", FontWeight->"Bold"], " has." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Head[Expression1]\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \(Equal\)], "Output"] }, Open ]], Cell[TextData[{ "In ", StyleBox["Mathematica", FontSlant->"Italic"], ", functions can be written in the ", StyleBox["Prefix", FontWeight->"Bold"], " form as well as the ", StyleBox["Postfix", FontWeight->"Bold"], " form. As an example for showing the two (equivalent) ways to input \ expressions, evaluate and compare the following:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(\(Head[Expression1]\)\(\[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(Expression1 // Head\)}], "Input"], Cell[BoxData[ \(Equal\)], "Output"], Cell[BoxData[ \(Equal\)], "Output"] }, Open ]], Cell[TextData[{ "What does it mean that objects are expressions? To look at how ", StyleBox["Expression1", FontWeight->"Bold"], " is represented, use the command ", StyleBox["FullForm", FontWeight->"Bold"], ":" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Expression1 // FullForm\)], "Input"], Cell[BoxData[ TagBox[ StyleBox[\(Equal[ Times[Plus[1, Times[Power[k1f, \(-1\)], Plus[k1b, k2f]]], EnzS[t]], S[t]]\), ShowSpecialCharacters->False, ShowStringCharacters->True, NumberMarks->True], FullForm]], "Output"] }, Open ]], Cell[TextData[{ "As another (visual) representation of the same expression, use ", StyleBox["TreeForm", FontWeight->"Bold"] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Expression1 // TreeForm\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Equal\"\>", "[", RowBox[{ InterpretationBox[GridBox[{ {"\<\"|\"\>"}, { RowBox[{"\<\"Times\"\>", "[", RowBox[{GridBox[{ {"\<\"|\"\>"}, { RowBox[{"\<\"Plus\"\>", "[", RowBox[{"1", ",", GridBox[{ {"\<\"|\"\>"}, { RowBox[{"\<\"Times\"\>", "[", RowBox[{GridBox[{ {"\<\"|\"\>"}, {\("Power"[k1f, \(-1\)]\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ",", GridBox[{ {"\<\"|\"\>"}, {\("Plus"[k1b, k2f]\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}]}], "]"}]} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}]}], "]"}]} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ",", GridBox[{ {"\<\"|\"\>"}, {\("EnzS"[t]\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}]}], "]"}]} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ColumnForm[ {"|", "Times"[ ColumnForm[ {"|", "Plus"[ 1, ColumnForm[ {"|", "Times"[ ColumnForm[ {"|", "Power"[ k1f, -1]}], ColumnForm[ {"|", "Plus"[ k1b, k2f]}]]}]]}], ColumnForm[ {"|", "EnzS"[ t]}]]}], Editable->False], ",", InterpretationBox[GridBox[{ {"\<\"|\"\>"}, {\("S"[t]\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ColumnForm[ {"|", "S"[ t]}], Editable->False]}], "]"}], TreeForm[ Equal[ Times[ Plus[ 1, Times[ Power[ k1f, -1], Plus[ k1b, k2f]]], EnzS[ t]], S[ t]]], Editable->False]], "Output"] }, Open ]], Cell[TextData[{ "How do we extract elements of expressions? Use commands ", StyleBox["Last", FontWeight->"Bold"], " and ", StyleBox["First", FontWeight->"Bold"], ". Find out what they do below, and find out what are the ", StyleBox["Last", FontWeight->"Bold"], " and ", StyleBox["Head", FontWeight->"Bold"], " of ", StyleBox["Expression1", FontWeight->"Bold"], "." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(Expression1 // First\[IndentingNewLine]\), "\[IndentingNewLine]", \(Expression1 // Last\)}], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \(\((1 + \(k1b + k2f\)\/k1f)\)\ EnzS[t]\)], "Output"], Cell[BoxData[ \(S[t]\)], "Output"] }, Open ]], Cell[TextData[{ "The command ", StyleBox["Together", FontWeight->"Bold"], " is useful in cancelling factors in a common denominator. Find out what it \ does, and factor the left-hand-side of ", StyleBox["Expression1", FontWeight->"Bold"], ", using combination of ", StyleBox["First", FontWeight->"Bold"], " and ", StyleBox["Together", FontWeight->"Bold"], ". " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(?Together\)\[IndentingNewLine]\), "\[IndentingNewLine]", \(Together[Expression1 // First]\)}], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ RowBox[{"\<\"Together[expr] puts terms in a sum over a common \ denominator, and cancels factors in the result.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"Together", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384087850-7840154"], Cell[BoxData[ \(\(\((k1b + k1f + k2f)\)\ EnzS[t]\)\/k1f\)], "Output"] }, Open ]], Cell[TextData[{ "For comparison, we see that applying ", StyleBox["Together", FontWeight->"Bold"], " on ", StyleBox["Expression 1", FontWeight->"Bold"], " fails to get the desired result." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Expression1 // Together\)], "Input"], Cell[BoxData[ \(\((1 + \(k1b + k2f\)\/k1f)\)\ EnzS[t] \[Equal] S[t]\)], "Output"] }, Open ]], Cell[TextData[{ "The problem arises from the fact that the term ", Cell[BoxData[ \(\((1 + \(k1b + k2f\)\/k1f)\)\)]], "is a subexpression, i.e., not at the first level of ", StyleBox["Expression1", FontWeight->"Bold"], ". To obtain the desired result, we should map ", StyleBox["Together", FontWeight->"Bold"], " to all subexpressions. Find out what ", StyleBox["MapAll", FontWeight->"Bold"], " does, and apply it to ", StyleBox["Expression1", FontWeight->"Bold"], ". Complete the following:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?MapAll\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"MapAll[f, expr] or f //@ expr applies f to every \ subexpression in expr.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"MapAll", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384087862-1699713"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(Together\)\(//@\)\(\ \ \)\(Expression1\)\(\ \ \ \)\)\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \(\(\((k1b + k1f + k2f)\)\ EnzS[t]\)\/k1f \[Equal] S[t]\)], "Output"] }, Open ]], Cell["\<\ The command Thread is useful in writing down a system of equations, \ equating a list of left-hand-side to a list of right-hand-side. Simply \ equating the lists do not work:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(LHS = {a, b, c};\)\), "\[IndentingNewLine]", \(\(\(RHS = {1, 2, 3};\)\(\[IndentingNewLine]\) \)\), "\[IndentingNewLine]", \(LHS \[Equal] RHS\)}], "Input"], Cell[BoxData[ \({a, b, c} \[Equal] {1, 2, 3}\)], "Output"] }, Open ]], Cell[TextData[{ "To equate the individual terms, one needs to apply ", StyleBox["Thread", FontWeight->"Bold"], "; see what happens:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Thread[{a, b, c} == {1, 2, 3}]\)], "Input"], Cell[BoxData[ \({a \[Equal] 1, b \[Equal] 2, c \[Equal] 3}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Perturbation Analysis", "Subsection"], Cell["\<\ Suppose one has algebraic equation obtained from perturbative \ analysis, e.g.,\ \>", "Text"], Cell[BoxData[ \(\(PertEqn = \((La\ + \ Lb\ \[Epsilon])\)*\((Lc + Ld\ \[Epsilon])\) + Lf\ \[Epsilon]\^2 \[Equal] \ Ra\ + \ \((Rb + Ra)\)*\[Epsilon]\ + \ Rc\ \[Epsilon]\^2;\)\)], "Input"], Cell["\<\ One would like to equate equal powers of left and right-hand sides.\ \ \>", "Text"], Cell["Firstly, obtain the LHS of PertEqn below:", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(PertEqnLHS\ = \ \ \ \ PertEqn // First\[IndentingNewLine]\), "\[IndentingNewLine]", \(PertEqnRHS\ = \ \ \ \ PertEqn // Last\)}], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \(Lf\ \[Epsilon]\^2 + \((La + Lb\ \[Epsilon])\)\ \((Lc + Ld\ \[Epsilon])\)\)], "Output"], Cell[BoxData[ \(Ra + \((Ra + Rb)\)\ \[Epsilon] + Rc\ \[Epsilon]\^2\)], "Output"] }, Open ]], Cell[TextData[{ "The terms of each order in ", "\[Epsilon] ", "can be obtained using the command ", StyleBox["CoefficientList", FontWeight->"Bold"] }], "Text"], Cell[BoxData[ \(\(?CoefficientList\)\)], "Input"], Cell[TextData[{ "Using ", StyleBox["CoefficientList", FontWeight->"Bold"], ", obtain a list of equations matching LHS and RHS, for each order in ", StyleBox["\[Epsilon]", FontWeight->"Bold"], ". Hint: use ", StyleBox["Thread", FontWeight->"Bold"], " to equate each term in the LHS and RHS lists. Complete the following:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(LHSCoefficientList = CoefficientList[\ PertEqnLHS, \ \[Epsilon]]\), "\[IndentingNewLine]", \(RHSCoefficientList = CoefficientList[\ PertEqnRHS, \ \[Epsilon]]\ \[IndentingNewLine]\), \ "\[IndentingNewLine]", \(OrderEqn\ = \ \ Thread[\ \ LHSCoefficientList\ \[Equal] \ \ \ RHSCoefficientList\ ]\)}], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \({La\ Lc, Lb\ Lc + La\ Ld, Lb\ Ld + Lf}\)], "Output"], Cell[BoxData[ \({Ra, Ra + Rb, Rc}\)], "Output"], Cell[BoxData[ \({La\ Lc \[Equal] Ra, Lb\ Lc + La\ Ld \[Equal] Ra + Rb, Lb\ Ld + Lf \[Equal] Rc}\)], "Output"] }, Open ]], Cell[TextData[{ "\nNow, let's extract the RHS unknowns, based on the following hint using \ the command ", StyleBox["Cases", FontWeight->"Bold"] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?Cases\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"Cases[{e1, e2, ... }, pattern] gives a list of the ei that \ match the pattern. Cases[{e1, ... }, pattern -> rhs] gives a list of the \ values of rhs corresponding to the ei that match the pattern. Cases[expr, \ pattern, levspec] gives a list of all parts of expr on levels specified by \ levspec which match the pattern. Cases[expr, pattern -> rhs, levspec] gives \ the values of rhs which match the pattern. Cases[expr, pattern, levspec, n] \ gives the first n parts in expr which match the pattern.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"Cases", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384088001-7593916"] }, Open ]], Cell[TextData[{ "As an example on using ", StyleBox["Cases", FontWeight->"Bold"], " to extract all symbols from a expression, evaluate below:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Cases[{a + b, \ c*d/f + f, \ 2.0*5 + g}, \ _Symbol, \ Infinity]\)], "Input"], Cell[BoxData[ \({a, b, c, d, f, f, g}\)], "Output"] }, Open ]], Cell[TextData[{ "Try a call to Cases, apply it to ", StyleBox["RHSCoefficientList", FontWeight->"Bold"] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(RHS\[UnderBracket]Unknowns\)\(=\)\(\ \ \ \)\(Cases[ RHSCoefficientList\ , \ _Symbol, \ Infinity]\)\(\ \)\)\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \({Ra, Ra, Rb, Rc}\)], "Output"] }, Open ]], Cell["Remove redundant elements, by using the command Union:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?Union\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"Union[list1, list2, ... ] gives a sorted list of all the \ distinct elements that appear in any of the listi. Union[list] gives a sorted \ version of a list, in which all duplicated elements have been dropped.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"Union", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384088126-3864133"] }, Open ]], Cell[TextData[{ "Apply ", StyleBox["Union ", FontWeight->"Bold"], "to ", StyleBox["RHS\[UnderBracket]Unknowns", FontWeight->"Bold"], StyleBox[", and over-write the result:", FontVariations->{"CompatibilityType"->0}] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(RHS\[UnderBracket]Unknowns = \ \ \(\(RHS\[UnderBracket]Unknowns\)\(//\)\ \(Union\)\(\ \ \ \)\)\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \({Ra, Rb, Rc}\)], "Output"] }, Open ]], Cell[TextData[{ "Let's now solve ", StyleBox["OrderEqn", FontWeight->"Bold"], ", with respect to ", StyleBox["RHS\[UnderBracket]Unknowns", FontWeight->"Bold"], ", using the command ", StyleBox["Solve", FontWeight->"Bold"], ":" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?Solve\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"Solve[eqns, vars] attempts to solve an equation or set of \ equations for the variables vars. Solve[eqns, vars, elims] attempts to solve \ the equations for vars, eliminating the variables elims.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"Solve", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384088145-4635598"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(RHS\[UnderBracket]Soln = \ Solve[\ OrderEqn, \ \ RHS\[UnderBracket]Unknowns\ ];\)\), "\ \[IndentingNewLine]", \(\(\(%\)\(//\)\(ColumnForm\)\(\ \ \ \ \ \ \ \)\)\)}], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ InterpretationBox[GridBox[{ {\({Ra \[Rule] La\ Lc, Rb \[Rule] \(-La\)\ Lc + Lb\ Lc + La\ Ld, Rc \[Rule] Lb\ Ld + Lf}\)} }, GridBaseline->{Baseline, {1, 1}}, ColumnAlignments->{Left}], ColumnForm[ {{Ra -> Times[ La, Lc], Rb -> Plus[ Times[ -1, La, Lc], Times[ Lb, Lc], Times[ La, Ld]], Rc -> Plus[ Times[ Lb, Ld], Lf]}}], Editable->False]], "Output"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Numerical Integration", "Section"], Cell[TextData[{ "To numerically integrate ordinary differential equations, the function ", StyleBox["NDSolve", FontWeight->"Bold"], " is used." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?NDSolve\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"NDSolve[eqns, y, {x, xmin, xmax}] finds a numerical solution \ to the ordinary differential equations eqns for the function y with the \ independent variable x in the range xmin to xmax. NDSolve[eqns, y, {x, xmin, \ xmax}, {t, tmin, tmax}] finds a numerical solution to the partial \ differential equations eqns. NDSolve[eqns, {y1, y2, ... }, {x, xmin, xmax}] \ finds numerical solutions for the functions yi.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"NDSolve", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384088179-1846413"] }, Open ]], Cell[TextData[{ "Lets now solve the ODE system\n\n", Cell[BoxData[ FormBox[ StyleBox[\(dy\/dt = \(-y\) + 1\), FontColor->RGBColor[0, 0, 1]], TraditionalForm]]], "\n\nsubject to the initial condition ", StyleBox["y[0]=2", FontColor->RGBColor[0, 0, 1]], ". \n\nComplete the following" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?D\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"D[f, x] gives the partial derivative of f with respect to x. \ D[f, {x, n}] gives the nth partial derivative of f with respect to x. D[f, \ x1, x2, ... ] gives a mixed derivative.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"D", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384088181-7830048"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ODE = D[\ \ y[t], \ t]\ \[Equal] \ \(\(-\ y[t]\)\(\ \)\(+\)\(\ \)\(1\)\(\ \ \ \)\)\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["y", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(1 - y[t]\)}]], "Output"] }, Open ]], Cell["\<\ Also, we need the initial condition. Complete below:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(IC\ = \ \ \(\(y[0]\)\(\[Equal]\)\(2\)\(\ \)\)\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ \(y[0] \[Equal] 2\)], "Output"] }, Open ]], Cell[TextData[{ "\nNow, combine ", StyleBox["ODE", FontWeight->"Bold"], " and ", StyleBox["IC", FontWeight->"Bold"], " into a single List, and call it ", StyleBox["ODE\[UnderBracket]IC", FontWeight->"Bold"] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(ODE\[UnderBracket]IC = \ \ {ODE, \ IC}\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["y", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Equal]", \(1 - y[t]\)}], ",", \(y[0] \[Equal] 2\)}], "}"}]], "Output"] }, Open ]], Cell[TextData[{ "Now, using ", StyleBox["NDSolve", FontWeight->"Bold"], ", compute the solution from t=0 to t=5." }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(NumericSoln = \ \ NDSolve[ODE\[UnderBracket]IC, \ y, \ {t, 0, 5}\ ]\)], "Input", Background->RGBColor[0, 1, 1]], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"y", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 5.`}}, "<>"]\), False, Editable->False]}], "}"}], "}"}]], "Output"] }, Open ]], Cell["Lets now plot the solution.", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(?Plot\)\)], "Input"], Cell[BoxData[ RowBox[{"\<\"Plot[f, {x, xmin, xmax}] generates a plot of f as a function \ of x from xmin to xmax. Plot[{f1, f2, ... }, {x, xmin, xmax}] plots several \ functions fi.\"\>", " ", ButtonBox[ StyleBox["More\[Ellipsis]", "SR"], ButtonData:>"Plot", Active->True, ButtonStyle->"RefGuideLink"]}]], "Print", CellTags->"Info3384088259-2927757"] }, Open ]], Cell["Complete the following:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[\ y[t] /. NumericSoln, \ \ \ {t, 0, 5}];\)\)], "Input", Background->RGBColor[0, 1, 1]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.190476 -0.581875 0.592597 [ [.21429 -0.00178 -3 -9 ] [.21429 -0.00178 3 0 ] [.40476 -0.00178 -3 -9 ] [.40476 -0.00178 3 0 ] [.59524 -0.00178 -3 -9 ] [.59524 -0.00178 3 0 ] [.78571 -0.00178 -3 -9 ] [.78571 -0.00178 3 0 ] [.97619 -0.00178 -3 -9 ] [.97619 -0.00178 3 0 ] [.01131 .12924 -18 -4.5 ] [.01131 .12924 0 4.5 ] [.01131 .24776 -18 -4.5 ] [.01131 .24776 0 4.5 ] [.01131 .36628 -18 -4.5 ] [.01131 .36628 0 4.5 ] [.01131 .4848 -18 -4.5 ] [.01131 .4848 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .01072 m .21429 .01697 L s [(1)] .21429 -0.00178 0 1 Mshowa .40476 .01072 m .40476 .01697 L s [(2)] .40476 -0.00178 0 1 Mshowa .59524 .01072 m .59524 .01697 L s [(3)] .59524 -0.00178 0 1 Mshowa .78571 .01072 m .78571 .01697 L s [(4)] .78571 -0.00178 0 1 Mshowa .97619 .01072 m .97619 .01697 L s [(5)] .97619 -0.00178 0 1 Mshowa .125 Mabswid .0619 .01072 m .0619 .01447 L s .1 .01072 m .1 .01447 L s .1381 .01072 m .1381 .01447 L s .17619 .01072 m .17619 .01447 L s .25238 .01072 m .25238 .01447 L s .29048 .01072 m .29048 .01447 L s .32857 .01072 m .32857 .01447 L s .36667 .01072 m .36667 .01447 L s .44286 .01072 m .44286 .01447 L s .48095 .01072 m .48095 .01447 L s .51905 .01072 m .51905 .01447 L s .55714 .01072 m .55714 .01447 L s .63333 .01072 m .63333 .01447 L s .67143 .01072 m .67143 .01447 L s .70952 .01072 m .70952 .01447 L s .74762 .01072 m .74762 .01447 L s .82381 .01072 m .82381 .01447 L s .8619 .01072 m .8619 .01447 L s .9 .01072 m .9 .01447 L s .9381 .01072 m .9381 .01447 L s .25 Mabswid 0 .01072 m 1 .01072 L s .02381 .12924 m .03006 .12924 L s [(1.2)] .01131 .12924 1 0 Mshowa .02381 .24776 m .03006 .24776 L s [(1.4)] .01131 .24776 1 0 Mshowa .02381 .36628 m .03006 .36628 L s [(1.6)] .01131 .36628 1 0 Mshowa .02381 .4848 m .03006 .4848 L s [(1.8)] .01131 .4848 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(2)] .01131 .60332 1 0 Mshowa .125 Mabswid .02381 .04035 m .02756 .04035 L s .02381 .06998 m .02756 .06998 L s .02381 .09961 m .02756 .09961 L s .02381 .15887 m .02756 .15887 L s .02381 .1885 m .02756 .1885 L s .02381 .21813 m .02756 .21813 L s .02381 .27739 m .02756 .27739 L s .02381 .30702 m .02756 .30702 L s .02381 .33665 m .02756 .33665 L s .02381 .39591 m .02756 .39591 L s .02381 .42554 m .02756 .42554 L s .02381 .45517 m .02756 .45517 L s .02381 .51443 m .02756 .51443 L s .02381 .54406 m .02756 .54406 L s .02381 .57369 m .02756 .57369 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .60332 m .06244 .49453 L .10458 .39851 L .14415 .32577 L .18221 .26871 L .22272 .21929 L .26171 .18068 L .30316 .14744 L .34309 .12159 L .3815 .10134 L .42237 .08384 L .46172 .07019 L .49955 .05948 L .53984 .05018 L .57861 .04292 L .61984 .03665 L .65954 .03177 L .69774 .02795 L .73838 .02464 L .77751 .02205 L .81909 .01983 L .85916 .0181 L .89771 .01675 L .93871 .01558 L .97619 .01472 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{35, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg00fMWI000Q0=WIf@030000003IfMT0fMWI0:H0fMWI0`00 003S0=WIf@008@3IfMT00`000000fMWI0=WIf@2T0=WIf@800000iP3IfMT00240fMWI00<000000=WI f@3IfMT0X@3IfMT300000>P0fMWI000Q0=WIf@030000003IfMT0fMWI09h0fMWI0`00003[0=WIf@00 8@3IfMT00`000000fMWI0=WIf@2L0=WIf@800000kP3IfMT00240fMWI00<000000=WIf@3IfMT0V@3I fMT300000?00fMWI000Q0=WIf@030000003IfMT0fMWI09L0fMWI0P00003c0=WIf@008@3IfMT00`00 0000fMWI0=WIf@2E0=WIf@800000m@3IfMT00240fMWI00<000000=WIf@3IfMT0TP3IfMT300000?L0 fMWI000Q0=WIf@800000T@3IfMT200000?X0fMWI000Q0=WIf@030000003IfMT0fMWI08h0fMWI0P00 003l0=WIf@008@3IfMT00`000000fMWI0=WIf@2<0=WIf@800000oP3IfMT00240fMWI00<000000=WI f@3IfMT0RP3IfMT200000?l0fMWI0@3IfMT00240fMWI00<000000=WIf@3IfMT0R03IfMT200000?l0 fMWI0`3IfMT00240fMWI00<000000=WIf@3IfMT0QP3IfMT200000?l0fMWI1@3IfMT00240fMWI00<0 00000=WIf@3IfMT0Q03IfMT200000?l0fMWI1`3IfMT00240fMWI00<000000=WIf@3IfMT0PP3IfMT2 00000?l0fMWI2@3IfMT00080fMWI1`0000050=WIf@800000103IfMT6000000L0fMWI00<000000=WI f@3IfMT0P03IfMT200000?l0fMWI2`3IfMT000D0fMWI00<000000=WIf@3IfMT01P3IfMT2000000D0 fMWI00D000000=WIf@3IfMT0fMWI000000070=WIf@030000003IfMT0fMWI07l0fMWI00<000000=WI f@3IfMT0o`3IfMT;0=WIf@001@3IfMT00`000000fMWI0=WIf@0>0=WIf@040000003IfMT0fMWI0000 00L0fMWI00<000000=WIf@3IfMT0O@3IfMT200000?l0fMWI3P3IfMT000D0fMWI00<000000=WIf@3I fMT03`3IfMT00`000000fMWI0=WIf@070=WIf@030000003IfMT0fMWI07`0fMWI00<000000=WIf@3I fMT0o`3IfMT>0=WIf@001@3IfMT00`000000fMWI0=WIf@0@0=WIf@030000003IfMT0fMWI00H0fMWI 0`00001j0=WIf@800000o`3IfMTA0=WIf@001@3IfMT00`000000fMWI0=WIf@0A0=WIf@030000003I fMT0fMWI00D0fMWI00<000000=WIf@3IfMT0N@3IfMT00`000000fMWI0=WIf@3o0=WIfA40fMWI0005 0=WIf@030000003IfMT0fMWI0140fMWI00<000000=WIf@3IfMT01@3IfMT00`000000fMWI0=WIf@1g 0=WIf@800000o`3IfMTD0=WIf@001@3IfMT00`000000fMWI0=WIf@0A0=WIf@030000003IfMT0fMWI 00D0fMWI00<000000=WIf@3IfMT0M@3IfMT200000?l0fMWI5P3IfMT00080fMWI0P0000000`3IfMT0 00000=WIf@0=0=WIf@040000003IfMT0fMWI0=WIf@8000001`3IfMT00`000000fMWI0=WIf@1d0=WI f@030000003IfMT0fMWI0?l0fMWI5P3IfMT000@0fMWI0P00000?0=WIf@@00000203IfMT00`000000 fMWI0=WIf@1b0=WIf@800000o`3IfMTI0=WIf@008@3IfMT00`000000fMWI0=WIf@1a0=WIf@030000 003IfMT0fMWI0?l0fMWI6@3IfMT00240fMWI00<000000=WIf@3IfMT0K`3IfMT200000?l0fMWI703I fMT00240fMWI00<000000=WIf@3IfMT0KP3IfMT00`000000fMWI0=WIf@3o0=WIfA`0fMWI000Q0=WI f@030000003IfMT0fMWI06d0fMWI00<000000=WIf@3IfMT0o`3IfMTM0=WIf@008@3IfMT00`000000 fMWI0=WIf@1[0=WIf@800000o`3IfMTP0=WIf@008@3IfMT00`000000fMWI0=WIf@1Z0=WIf@030000 003IfMT0fMWI0?l0fMWI803IfMT00240fMWI0P00001Z0=WIf@030000003IfMT0fMWI0?l0fMWI8@3I fMT00240fMWI00<000000=WIf@3IfMT0I`3IfMT200000?l0fMWI903IfMT00240fMWI00<000000=WI f@3IfMT0IP3IfMT00`000000fMWI0=WIf@3o0=WIfB@0fMWI000Q0=WIf@030000003IfMT0fMWI06D0 fMWI00<000000=WIf@3IfMT0o`3IfMTU0=WIf@008@3IfMT00`000000fMWI0=WIf@1S0=WIf@800000 o`3IfMTX0=WIf@008@3IfMT00`000000fMWI0=WIf@1R0=WIf@030000003IfMT0fMWI0?l0fMWI:03I fMT00240fMWI00<000000=WIf@3IfMT0H@3IfMT00`000000fMWI0=WIf@3o0=WIfBT0fMWI000Q0=WI f@030000003IfMT0fMWI05l0fMWI0P00003o0=WIfB`0fMWI000Q0=WIf@030000003IfMT0fMWI05h0 fMWI00<000000=WIf@3IfMT0o`3IfMT/0=WIf@008@3IfMT00`000000fMWI0=WIf@1M0=WIf@030000 003IfMT0fMWI0?l0fMWI;@3IfMT00240fMWI00<000000=WIf@3IfMT0G03IfMT00`000000fMWI0=WI f@3o0=WIfBh0fMWI000Q0=WIf@030000003IfMT0fMWI05/0fMWI00<000000=WIf@3IfMT0o`3IfMT_ 0=WIf@008@3IfMT2000005/0fMWI00<000000=WIf@3IfMT0o`3IfMT`0=WIf@008@3IfMT00`000000 fMWI0=WIf@1I0=WIf@030000003IfMT0fMWI0?l0fMWI<@3IfMT00240fMWI00<000000=WIf@3IfMT0 F03IfMT00`000000fMWI0=WIf@3o0=WIfC80fMWI000Q0=WIf@030000003IfMT0fMWI05L0fMWI00<0 00000=WIf@3IfMT0o`3IfMTc0=WIf@008@3IfMT00`000000fMWI0=WIf@1F0=WIf@030000003IfMT0 fMWI0?l0fMWI=03IfMT00240fMWI00<000000=WIf@3IfMT0E@3IfMT00`000000fMWI0=WIf@3o0=WI fCD0fMWI000Q0=WIf@030000003IfMT0fMWI05D0fMWI00<000000=WIf@3IfMT0o`3IfMTe0=WIf@00 8@3IfMT00`000000fMWI0=WIf@1D0=WIf@030000003IfMT0fMWI0?l0fMWI=P3IfMT00240fMWI00<0 00000=WIf@3IfMT0D`3IfMT00`000000fMWI0=WIf@3o0=WIfCL0fMWI000Q0=WIf@030000003IfMT0 fMWI0580fMWI00<000000=WIf@3IfMT0o`3IfMTh0=WIf@008@3IfMT00`000000fMWI0=WIf@1A0=WI f@030000003IfMT0fMWI0?l0fMWI>@3IfMT00240fMWI00<000000=WIf@3IfMT0D03IfMT00`000000 fMWI0=WIf@3o0=WIfCX0fMWI000Q0=WIf@800000D03IfMT00`000000fMWI0=WIf@3o0=WIfC/0fMWI 000Q0=WIf@030000003IfMT0fMWI04h0fMWI00<000000=WIf@3IfMT0o`3IfMTl0=WIf@008@3IfMT0 0`000000fMWI0=WIf@1=0=WIf@030000003IfMT0fMWI0?l0fMWI?@3IfMT00240fMWI00<000000=WI f@3IfMT0C03IfMT00`000000fMWI0=WIf@3o0=WIfCh0fMWI000Q0=WIf@030000003IfMT0fMWI04/0 fMWI00<000000=WIf@3IfMT0o`3IfMTo0=WIf@008@3IfMT00`000000fMWI0=WIf@1;0=WIf@030000 003IfMT0fMWI0?l0fMWI?`3IfMT00240fMWI00<000000=WIf@3IfMT0BP3IfMT00`000000fMWI0=WI f@3o0=WIfD00fMWI000Q0=WIf@030000003IfMT0fMWI04T0fMWI00<000000=WIf@3IfMT0o`3IfMU1 0=WIf@000P3IfMT7000000D0fMWI0P0000060=WIf@@000001`3IfMT00`000000fMWI0=WIf@180=WI f@030000003IfMT0fMWI0?l0fMWI@P3IfMT000D0fMWI00<000000=WIf@3IfMT01P3IfMT2000000P0 fMWI00<000000=WIf@3IfMT01P3IfMT00`000000fMWI0=WIf@170=WIf@030000003IfMT0fMWI0?l0 fMWI@`3IfMT000D0fMWI00<000000=WIf@3IfMT0403IfMT00`000000fMWI0=WIf@060=WIf@030000 003IfMT0fMWI04H0fMWI00<000000=WIf@3IfMT0o`3IfMU40=WIf@001@3IfMT00`000000fMWI0=WI f@0<0=WIf@H000001`3IfMT00`000000fMWI0=WIf@150=WIf@030000003IfMT0fMWI0?l0fMWIA@3I fMT000D0fMWI00<000000=WIf@3IfMT0303IfMT01@000000fMWI0=WIf@3IfMT0000000P0fMWI0`00 00150=WIf@030000003IfMT0fMWI0?l0fMWIA@3IfMT000D0fMWI00<000000=WIf@3IfMT03@3IfMT0 10000000fMWI0=WIf@0000080=WIf@030000003IfMT0fMWI04@0fMWI00<000000=WIf@3IfMT0o`3I fMU60=WIf@001@3IfMT00`000000fMWI0=WIf@0>0=WIf@030000003IfMT0000000P0fMWI00<00000 0=WIf@3IfMT0@`3IfMT00`000000fMWI0=WIf@3o0=WIfDL0fMWI00050=WIf@030000003IfMT0fMWI 00h0fMWI00<000000=WIf@000000203IfMT00`000000fMWI0=WIf@120=WIf@030000003IfMT0fMWI 0?l0fMWIB03IfMT00080fMWI0P0000000`3IfMT000000=WIf@0@0=WIf@800000203IfMT00`000000 fMWI0=WIf@110=WIf@030000003IfMT0fMWI0?l0fMWIB@3IfMT000@0fMWI0P00000A0=WIf@800000 203IfMT00`000000fMWI0=WIf@100=WIf@030000003IfMT0fMWI0?l0fMWIBP3IfMT00240fMWI00<0 00000=WIf@3IfMT0@03IfMT00`000000fMWI0=WIf@3o0=WIfDX0fMWI000Q0=WIf@030000003IfMT0 fMWI03l0fMWI00<000000=WIf@3IfMT0o`3IfMU;0=WIf@008@3IfMT00`000000fMWI0=WIf@0n0=WI f@030000003IfMT0fMWI0?l0fMWIC03IfMT00240fMWI00<000000=WIf@3IfMT0?@3IfMT00`000000 fMWI0=WIf@3o0=WIfDd0fMWI000Q0=WIf@030000003IfMT0fMWI03`0fMWI00<000000=WIf@3IfMT0 o`3IfMU>0=WIf@008@3IfMT00`000000fMWI0=WIf@0l0=WIf@030000003IfMT0fMWI0?l0fMWICP3I fMT00240fMWI0P00000l0=WIf@030000003IfMT0fMWI0?l0fMWIC`3IfMT00240fMWI00<000000=WI f@3IfMT0>P3IfMT00`000000fMWI0=WIf@3o0=WIfE00fMWI000Q0=WIf@030000003IfMT0fMWI03X0 fMWI00<000000=WIf@3IfMT0o`3IfMU@0=WIf@008@3IfMT00`000000fMWI0=WIf@0i0=WIf@030000 003IfMT0fMWI0?l0fMWID@3IfMT00240fMWI00<000000=WIf@3IfMT0>03IfMT00`000000fMWI0=WI f@3o0=WIfE80fMWI000Q0=WIf@030000003IfMT0fMWI03P0fMWI00<000000=WIf@3IfMT0o`3IfMUB 0=WIf@008@3IfMT00`000000fMWI0=WIf@0g0=WIf@030000003IfMT0fMWI0?l0fMWID`3IfMT00240 fMWI00<000000=WIf@3IfMT0=P3IfMT00`000000fMWI0=WIf@3o0=WIfE@0fMWI000Q0=WIf@030000 003IfMT0fMWI03H0fMWI00<000000=WIf@3IfMT0o`3IfMUD0=WIf@008@3IfMT00`000000fMWI0=WI f@0e0=WIf@030000003IfMT0fMWI0?l0fMWIE@3IfMT00240fMWI00<000000=WIf@3IfMT0=@3IfMT0 0`000000fMWI0=WIf@3o0=WIfED0fMWI000Q0=WIf@030000003IfMT0fMWI03@0fMWI00<000000=WI f@3IfMT0o`3IfMUF0=WIf@008@3IfMT2000003@0fMWI00<000000=WIf@3IfMT0o`3IfMUG0=WIf@00 8@3IfMT00`000000fMWI0=WIf@0c0=WIf@030000003IfMT0fMWI0?l0fMWIE`3IfMT00240fMWI00<0 00000=WIf@3IfMT00=WIf@005@3IfMT01@000000fMWI0=WI f@3IfMT0000000L0fMWI0P00003o0=WIfHh0fMWI000F0=WIf@040000003IfMT0fMWI000000L0fMWI 0P00003o0=WIfHh0fMWI000G0=WIf@030000003IfMT0fMWI00L0fMWI00<000000=WIf@3IfMT0o`3I fMV=0=WIf@00603IfMT00`000000fMWI0=WIf@060=WIf@<00000o`3IfMV=0=WIf@006@3IfMT00`00 0000fMWI0=WIf@050=WIf@030000003IfMT0fMWI0?l0fMWIS@3IfMT001T0fMWI00<000000=WIf@3I fMT01@3IfMT00`000000fMWI0=WIf@3o0=WIfHd0fMWI000I0=WIf@030000003IfMT0fMWI00D0fMWI 00<000000=WIf@3IfMT0o`3IfMV=0=WIf@00503IfMT010000000fMWI0=WIf@3IfMT2000000L0fMWI 00<000000=WIf@3IfMT0o`3IfMV=0=WIf@005@3IfMT4000000P0fMWI00<000000=WIf@3IfMT0o`3I fMV=0=WIf@008@3IfMT00`000000fMWI0=WIf@3o0=WIfHd0fMWI0000\ \>"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.432098, 0.922776, \ 0.0194004, 0.00623578}}] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"5.2 for X", ScreenRectangle->{{0, 1400}, {0, 1050}}, WindowSize->{1139, 942}, WindowMargins->{{17, Automatic}, {Automatic, 0}}, ShowSelection->True, ShowCellLabel->False, Magnification->1.5, StyleDefinitions -> "Report.nb" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{ "Info3384087771-5332558"->{ Cell[2654, 88, 285, 8, 55, "Print", CellTags->"Info3384087771-5332558"]}, "Info3384087850-7840154"->{ Cell[9096, 314, 353, 9, 78, "Print", CellTags->"Info3384087850-7840154"]}, "Info3384087862-1699713"->{ Cell[10561, 372, 326, 9, 55, "Print", CellTags->"Info3384087862-1699713"]}, "Info3384088001-7593916"->{ Cell[14506, 526, 753, 14, 193, "Print", CellTags->"Info3384088001-7593916"]}, "Info3384088126-3864133"->{ Cell[16180, 584, 468, 11, 101, "Print", CellTags->"Info3384088126-3864133"]}, "Info3384088145-4635598"->{ Cell[17496, 638, 447, 10, 101, "Print", CellTags->"Info3384088145-4635598"]}, "Info3384088179-1846413"->{ Cell[19042, 693, 663, 13, 170, "Print", CellTags->"Info3384088179-1846413"]}, "Info3384088181-7830048"->{ Cell[20119, 726, 427, 10, 101, "Print", CellTags->"Info3384088181-7830048"]}, "Info3384088259-2927757"->{ Cell[22563, 830, 405, 10, 78, "Print", CellTags->"Info3384088259-2927757"]} } *) (*CellTagsIndex CellTagsIndex->{ {"Info3384087771-5332558", 48788, 1385}, {"Info3384087850-7840154", 48902, 1388}, {"Info3384087862-1699713", 49017, 1391}, {"Info3384088001-7593916", 49133, 1394}, {"Info3384088126-3864133", 49251, 1397}, {"Info3384088145-4635598", 49369, 1400}, {"Info3384088179-1846413", 49487, 1403}, {"Info3384088181-7830048", 49605, 1406}, {"Info3384088259-2927757", 49723, 1409} } *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1776, 53, 94, 4, 125, "Title"], Cell[CellGroupData[{ Cell[1895, 61, 42, 0, 102, "Section"], Cell[1940, 63, 137, 3, 76, "Input"], Cell[2080, 68, 504, 13, 87, "Text"], Cell[CellGroupData[{ Cell[2609, 85, 42, 1, 46, "Input"], Cell[2654, 88, 285, 8, 55, "Print", CellTags->"Info3384087771-5332558"] }, Open ]], Cell[2954, 99, 117, 5, 41, "Text"], Cell[CellGroupData[{ Cell[3096, 108, 83, 2, 46, "Input"], Cell[3182, 112, 39, 1, 69, "Output"] }, Open ]], Cell[3236, 116, 371, 12, 64, "Text"], Cell[CellGroupData[{ Cell[3632, 132, 141, 3, 92, "Input"], Cell[3776, 137, 39, 1, 69, "Output"], Cell[3818, 140, 39, 1, 69, "Output"] }, Open ]], Cell[3872, 144, 240, 8, 41, "Text"], Cell[CellGroupData[{ Cell[4137, 156, 56, 1, 46, "Input"], Cell[4196, 159, 277, 8, 69, "Output"] }, Open ]], Cell[4488, 170, 143, 4, 41, "Text"], Cell[CellGroupData[{ Cell[4656, 178, 56, 1, 46, "Input"], Cell[4715, 181, 3035, 73, 179, "Output"] }, Open ]], Cell[7765, 257, 418, 17, 64, "Text"], Cell[CellGroupData[{ Cell[8208, 278, 160, 3, 92, "Input"], Cell[8371, 283, 71, 1, 91, "Output"], Cell[8445, 286, 38, 1, 69, "Output"] }, Open ]], Cell[8498, 290, 406, 15, 64, "Text"], Cell[CellGroupData[{ Cell[8929, 309, 164, 3, 92, "Input"], Cell[9096, 314, 353, 9, 78, "Print", CellTags->"Info3384087850-7840154"], Cell[9452, 325, 73, 1, 91, "Output"] }, Open ]], Cell[9540, 329, 218, 8, 41, "Text"], Cell[CellGroupData[{ Cell[9783, 341, 56, 1, 46, "Input"], Cell[9842, 344, 85, 1, 91, "Output"] }, Open ]], Cell[9942, 348, 547, 17, 95, "Text"], Cell[CellGroupData[{ Cell[10514, 369, 44, 1, 46, "Input"], Cell[10561, 372, 326, 9, 55, "Print", CellTags->"Info3384087862-1699713"] }, Open ]], Cell[CellGroupData[{ Cell[10924, 386, 122, 2, 46, "Input"], Cell[11049, 390, 87, 1, 91, "Output"] }, Open ]], Cell[11151, 394, 198, 4, 64, "Text"], Cell[CellGroupData[{ Cell[11374, 402, 191, 4, 115, "Input"], Cell[11568, 408, 62, 1, 69, "Output"] }, Open ]], Cell[11645, 412, 154, 5, 41, "Text"], Cell[CellGroupData[{ Cell[11824, 421, 63, 1, 46, "Input"], Cell[11890, 424, 76, 1, 69, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12003, 430, 43, 0, 62, "Subsection"], Cell[12049, 432, 103, 3, 41, "Text"], Cell[12155, 437, 229, 4, 51, "Input"], Cell[12387, 443, 93, 3, 41, "Text"], Cell[12483, 448, 57, 0, 41, "Text"], Cell[CellGroupData[{ Cell[12565, 452, 207, 4, 92, "Input"], Cell[12775, 458, 119, 2, 73, "Output"], Cell[12897, 462, 84, 1, 73, "Output"] }, Open ]], Cell[12996, 466, 171, 6, 41, "Text"], Cell[13170, 474, 53, 1, 46, "Input"], Cell[13226, 477, 358, 11, 64, "Text"], Cell[CellGroupData[{ Cell[13609, 492, 395, 9, 115, "Input"], Cell[14007, 503, 72, 1, 69, "Output"], Cell[14082, 506, 51, 1, 69, "Output"], Cell[14136, 509, 120, 2, 69, "Output"] }, Open ]], Cell[14271, 514, 164, 5, 64, "Text"], Cell[CellGroupData[{ Cell[14460, 523, 43, 1, 46, "Input"], Cell[14506, 526, 753, 14, 193, "Print", CellTags->"Info3384088001-7593916"] }, Open ]], Cell[15274, 543, 164, 5, 41, "Text"], Cell[CellGroupData[{ Cell[15463, 552, 103, 2, 46, "Input"], Cell[15569, 556, 55, 1, 69, "Output"] }, Open ]], Cell[15639, 560, 123, 4, 41, "Text"], Cell[CellGroupData[{ Cell[15787, 568, 184, 3, 46, "Input"], Cell[15974, 573, 50, 1, 69, "Output"] }, Open ]], Cell[16039, 577, 70, 0, 41, "Text"], Cell[CellGroupData[{ Cell[16134, 581, 43, 1, 46, "Input"], Cell[16180, 584, 468, 11, 101, "Print", CellTags->"Info3384088126-3864133"] }, Open ]], Cell[16663, 598, 249, 9, 41, "Text"], Cell[CellGroupData[{ Cell[16937, 611, 160, 3, 46, "Input"], Cell[17100, 616, 46, 1, 69, "Output"] }, Open ]], Cell[17161, 620, 264, 11, 41, "Text"], Cell[CellGroupData[{ Cell[17450, 635, 43, 1, 46, "Input"], Cell[17496, 638, 447, 10, 101, "Print", CellTags->"Info3384088145-4635598"] }, Open ]], Cell[CellGroupData[{ Cell[17980, 653, 243, 5, 69, "Input"], Cell[18226, 660, 474, 12, 69, "Output"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[18761, 679, 40, 0, 102, "Section"], Cell[18804, 681, 165, 5, 41, "Text"], Cell[CellGroupData[{ Cell[18994, 690, 45, 1, 46, "Input"], Cell[19042, 693, 663, 13, 170, "Print", CellTags->"Info3384088179-1846413"] }, Open ]], Cell[19720, 709, 332, 10, 187, "Text"], Cell[CellGroupData[{ Cell[20077, 723, 39, 1, 46, "Input"], Cell[20119, 726, 427, 10, 101, "Print", CellTags->"Info3384088181-7830048"] }, Open ]], Cell[CellGroupData[{ Cell[20583, 741, 180, 5, 46, "Input"], Cell[20766, 748, 180, 5, 69, "Output"] }, Open ]], Cell[20961, 756, 77, 3, 64, "Text"], Cell[CellGroupData[{ Cell[21063, 763, 112, 2, 46, "Input"], Cell[21178, 767, 49, 1, 69, "Output"] }, Open ]], Cell[21242, 771, 242, 10, 64, "Text"], Cell[CellGroupData[{ Cell[21509, 785, 104, 2, 46, "Input"], Cell[21616, 789, 275, 8, 69, "Output"] }, Open ]], Cell[21906, 800, 135, 5, 41, "Text"], Cell[CellGroupData[{ Cell[22066, 809, 142, 3, 46, "Input"], Cell[22211, 814, 224, 6, 69, "Output"] }, Open ]], Cell[22450, 823, 43, 0, 41, "Text"], Cell[CellGroupData[{ Cell[22518, 827, 42, 1, 46, "Input"], Cell[22563, 830, 405, 10, 78, "Print", CellTags->"Info3384088259-2927757"] }, Open ]], Cell[22983, 843, 39, 0, 41, "Text"], Cell[CellGroupData[{ Cell[23047, 847, 115, 2, 46, "Input"], Cell[23165, 851, 24894, 510, 291, 3453, 241, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)