(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 68795, 1816] NotebookOptionsPosition[ 65208, 1697] NotebookOutlinePosition[ 66037, 1727] CellTagsIndexPosition[ 65867, 1720] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ StyleBox[ RowBox[{ "Using", " ", "Mathematica", " ", "to", " ", "Solve", " ", "Differential", " ", "Equations"}], "Title"]], "Input", CellChangeTimes->{{3.4466800280615077`*^9, 3.446680033876565*^9}}], Cell[BoxData[ StyleBox[ RowBox[{ StyleBox[ RowBox[{ RowBox[{"NOTE", ":", " ", RowBox[{ "I", " ", "have", " ", "a", " ", "lot", " ", "of", " ", "difficulty", " ", "with", " ", "Mathematica", " ", "solving", " ", "differential", " ", RowBox[{"equations", ".", " ", "Frequently"}], " ", "I", " ", "find", " ", "that", " ", "if", " ", "I", " ", "make", " ", "a", " ", "typo", " ", "in", " ", "entering", " ", "the", " ", "equation", " ", "that", " ", "it", " ", "is", " ", "NOT", " ", "enough", " ", "to", " ", "fix", " ", "the", " ", "typo"}]}], ",", " ", RowBox[{ "but", " ", "rather", " ", "I", " ", "need", " ", "to", " ", "save", " ", "the", " ", "file"}], ",", " ", RowBox[{"close", " ", "Mathematica"}], ",", " ", RowBox[{"and", " ", "then", " ", RowBox[{"reload", "!"}]}]}], FontFamily->"Arial", FontColor->RGBColor[0.5019607843137255, 0., 0.]], "\[IndentingNewLine]"}], "Subsection"]], "Input", CellChangeTimes->{{3.446716340828904*^9, 3.4467163898482275`*^9}, { 3.4467167470987997`*^9, 3.4467167495146604`*^9}, {3.4467203899738207`*^9, 3.4467203916690245`*^9}, {3.4467204476000147`*^9, 3.4467204504953375`*^9}}], Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{ "In", " ", "these", " ", "notes", " ", "we", " ", "solve", " ", "various", " ", "types", " ", "of", " ", "differential", " ", RowBox[{"equations", ":", " ", RowBox[{ "with", " ", "and", " ", "without", " ", "initial", " ", "conditions"}]}]}], ",", " ", RowBox[{"systems", " ", "of", " ", "differential", " ", "equations"}], ",", " ", RowBox[{ RowBox[{"...", ".", " ", "We"}], " ", "also", " ", "solve", " ", "difference", " ", RowBox[{"equations", "."}]}]}], "Subsection"]], "Input", CellChangeTimes->{{3.446720392953122*^9, 3.446720432265059*^9}}], Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{"(*", " ", RowBox[{"NOTE", ":", " ", RowBox[{ "Anything", " ", "between", " ", "parentheses", " ", "star", " ", "and", " ", "star", " ", "parentheses", " ", "is", " ", "a", " ", "comment"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ "To", " ", "compile", " ", "a", " ", "line", " ", "type", " ", "shift"}], "-", "enter"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ "YOU", " ", "DO", " ", "NOT", " ", "NEED", " ", "TO", " ", "CLEAR", " ", "A", " ", "VARIABLE", " ", "BEFORE", " ", "USING", " ", "IT"}], ",", " ", RowBox[{"BUT", " ", RowBox[{"IT", "'"}], "S", " ", "OFTEN", " ", "A", " ", "GOOD", " ", "IDEA", " ", "SO", " ", "MATHEMATICA", " ", RowBox[{"DOESN", "'"}], "T", " ", "USE", " ", "AN", " ", "OLD", " ", "VALUE"}]}], " ", "*)"}]}], FontFamily->"Courier New", FontColor->RGBColor[0., 0., 1.]]], "Input", CellChangeTimes->{{3.446720491325732*^9, 3.4467204934427843`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Print", "[", RowBox[{ "\"\\"", RowBox[{ "Hyperlink", "[", "\"\\"", "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{ "\"\\"", ",", " ", RowBox[{ "Hyperlink", "[", "\"\\"", "]"}]}], "]"}], ";"}]}], "Input", CellChangeTimes->{{3.4467167506630244`*^9, 3.4467167900901885`*^9}, { 3.446716820443714*^9, 3.4467169365905776`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\<\"If you want to see my template for Mathematica, you can go to \ the following address (you can just click it below and your browser will open \ it!): \"\>", " ", TagBox[ ButtonBox[ PaneSelectorBox[{ False->"\<\"http://www.williams.edu/go/math/sjmiller/public_html/math/\ LaTexMathematica/MathematicaIntroVer6.nb\"\>", True-> StyleBox["\<\"http://www.williams.edu/go/math/sjmiller/public_html/math/\ LaTexMathematica/MathematicaIntroVer6.nb\"\>", "HyperlinkActive"]}, Dynamic[ CurrentValue["MouseOver"]], BaseStyle->{"Hyperlink"}, BaselinePosition->Baseline, FrameMargins->0, ImageSize->Automatic], BaseStyle->"Hyperlink", ButtonData->{ URL[ "http://www.williams.edu/go/math/sjmiller/public_html/math/\ LaTexMathematica/MathematicaIntroVer6.nb"], None}, ButtonNote-> "http://www.williams.edu/go/math/sjmiller/public_html/math/\ LaTexMathematica/MathematicaIntroVer6.nb"], Annotation[#, "http://www.williams.edu/go/math/sjmiller/public_html/math/\ LaTexMathematica/MathematicaIntroVer6.nb", "Hyperlink"]& ]}]], "Print", CellChangeTimes->{{3.446716859694132*^9, 3.4467168772271895`*^9}, 3.4467169380904818`*^9, 3.446720505193985*^9, 3.4467597423832054`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"In case that link doesn't let you download the code but rather \ opens the code in the browswer window, just go to: \"\>", "\[InvisibleSpace]", TagBox[ ButtonBox[ PaneSelectorBox[{ False->"\<\"http://www.williams.edu/go/math/sjmiller/public_html/math/\ handouts/latex.htm\"\>", True-> StyleBox["\<\"http://www.williams.edu/go/math/sjmiller/public_html/math/\ handouts/latex.htm\"\>", "HyperlinkActive"]}, Dynamic[ CurrentValue["MouseOver"]], BaseStyle->{"Hyperlink"}, BaselinePosition->Baseline, FrameMargins->0, ImageSize->Automatic], BaseStyle->"Hyperlink", ButtonData->{ URL[ "http://www.williams.edu/go/math/sjmiller/public_html/math/handouts/\ latex.htm"], None}, ButtonNote-> "http://www.williams.edu/go/math/sjmiller/public_html/math/handouts/\ latex.htm"], Annotation[#, "http://www.williams.edu/go/math/sjmiller/public_html/math/handouts/\ latex.htm", "Hyperlink"]& ]}], SequenceForm[ "In case that link doesn't let you download the code but rather opens the \ code in the browswer window, just go to: ", Hyperlink[ "http://www.williams.edu/go/math/sjmiller/public_html/math/handouts/latex.\ htm"]], Editable->False]], "Print", CellChangeTimes->{{3.446716859694132*^9, 3.4467168772271895`*^9}, 3.4467169380904818`*^9, 3.446720505193985*^9, 3.446759742388088*^9}] }, Open ]] }, Open ]], Cell[BoxData[""], "Input", CellChangeTimes->{{3.4467204849853177`*^9, 3.446720485370058*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "THESE", " ", "LINES", " ", "SOLVE", " ", "A", " ", "LINEAR", " ", "DIFFERENTIAL", " ", RowBox[{"EQUATION", ".", " ", "THE"}], " ", "FIRST", " ", RowBox[{"DOESN", "'"}], "T", " ", "HAVE", " ", "ANY", " ", "INITIAL", " ", "CONDITIONS", " ", "AND", " ", "THUS", " ", "GENERATES", " ", "ALL", " ", "SOLUTIONS"}], ";", " ", RowBox[{"THE", " ", "SECOND", " ", "HAS", " ", "AN", " ", "INITIAL", " ", RowBox[{"CONDITION", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", "y", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "x", "]"}], " ", "+", " ", RowBox[{"y", "[", "x", "]"}]}], " ", "==", " ", "x"}], ",", " ", RowBox[{"y", "[", "x", "]"}], ",", " ", "x"}], "]"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "x", "]"}], " ", "+", " ", RowBox[{"y", "[", "x", "]"}]}], " ", "==", " ", "x"}], ",", " ", RowBox[{ RowBox[{"y", "[", "0", "]"}], " ", "\[Equal]", " ", "1"}]}], "}"}], ",", " ", RowBox[{"y", "[", "x", "]"}], ",", " ", "x"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.4466800423652797`*^9, 3.4466800651743665`*^9}, { 3.4466801130853624`*^9, 3.4466801139212465`*^9}, {3.4466801899525127`*^9, 3.446680192761904*^9}, {3.446680323877535*^9, 3.446680325520008*^9}, { 3.4467155272589054`*^9, 3.446715583014126*^9}, {3.4467170583659863`*^9, 3.446717059470408*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], "\[Rule]", RowBox[{ RowBox[{"-", "1"}], "+", "x", "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", "x"}]], " ", RowBox[{"C", "[", "1", "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446680063701804*^9, 3.446680066156726*^9}, 3.4466801158381157`*^9, {3.4466801695378036`*^9, 3.446680193163245*^9}, 3.446680326717197*^9, 3.446680411024301*^9, 3.446715585795198*^9, 3.4467162323544016`*^9, 3.446720507207528*^9, 3.446759750418824*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], "\[Rule]", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", "x"}]], " ", RowBox[{"(", RowBox[{"2", "-", SuperscriptBox["\[ExponentialE]", "x"], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", "x"], " ", "x"}]}], ")"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446680063701804*^9, 3.446680066156726*^9}, 3.4466801158381157`*^9, {3.4466801695378036`*^9, 3.446680193163245*^9}, 3.446680326717197*^9, 3.446680411024301*^9, 3.446715585795198*^9, 3.4467162323544016`*^9, 3.446720507207528*^9, 3.446759750436401*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"IT", "'"}], "S", " ", "NOT", " ", "ALWAYS", " ", "EASY", " ", "TO", " ", "PULL", " ", "FUNCTIONS", " ", "OFF", " ", "OF", " ", "MATHEMATICA", " ", "WHEN", " ", "IT", " ", "SOLVES", " ", RowBox[{"SOMETHING", ".", " ", "SOMETIMES"}], " ", "THE", " ", "EASIEST", " ", "WAY", " ", "IS", " ", "TO", " ", "JUST", " ", "CUT", " ", "AND", " ", "PASTE"}], ",", " ", RowBox[{"WHICH", " ", "WE", " ", "DO", " ", RowBox[{"BELOW", ".", " ", "TO"}], " ", "DEFINE", " ", "A", " ", "FUNCTION"}], ",", " ", RowBox[{ "WE", " ", "HAVE", " ", "AN", " ", "UNDERSCORE", " ", "AFTER", " ", "EACH", " ", "VARIABLE", " ", "AND", " ", "A", " ", "COLON", " ", "BEFORE", " ", "THE", " ", "EQUAL", " ", RowBox[{"SIGN", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "IF", " ", "YOU", " ", "HAVE", " ", "A", " ", "SEMICOLON", " ", "AFTER", " ", "A", " ", "LINE", " ", "IT", " ", "SUPPRESSES", " ", "THE", " ", "OUTPUT"}], ";", " ", "HOWEVER"}], ",", " ", RowBox[{ "IT", " ", "ALWAYS", " ", "SUPPRESSES", " ", "OUTPUT", " ", "WHEN", " ", "DEFINING", " ", "A", " ", RowBox[{"FUNCTION", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"solution", "[", "x_", "]"}], " ", ":=", " ", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", "x"}]], " ", RowBox[{"(", RowBox[{"2", "-", SuperscriptBox["\[ExponentialE]", "x"], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", "x"], " ", "x"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{{3.446715593253705*^9, 3.4467157322399263`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "SOMETIMES", " ", "WHEN", " ", "IT", " ", "SOLVES", " ", "AN", " ", "EQUATION"}], ",", " ", RowBox[{"THE", " ", "RESULT", " ", "IS", " ", "NOT", " ", RowBox[{"SIMPLIFIED", ".", " ", "YOU"}], " ", "CAN", " ", "FIX", " ", "THIS", " ", "BY", " ", "USING", " ", "EXPAND", " ", "OR", " ", RowBox[{"SIMPLIFY", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Expand", "[", RowBox[{"solution", "[", "x", "]"}], "]"}], "\[IndentingNewLine]", RowBox[{"Simplify", "[", RowBox[{"solution", "[", "x", "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.4467157619470096`*^9, 3.4467157726806974`*^9}, { 3.4467162500788527`*^9, 3.446716251381504*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"2", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "x"}]]}], "+", "x"}]], "Output", CellChangeTimes->{3.446715773365224*^9, 3.4467162523667927`*^9, 3.4467205108137426`*^9, 3.4467597613839426`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"2", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "x"}]]}], "+", "x"}]], "Output", CellChangeTimes->{3.446715773365224*^9, 3.4467162523667927`*^9, 3.4467205108137426`*^9, 3.44675976139859*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "The", " ", "commands", " ", "below", " ", "allow", " ", "you", " ", "to", " ", "plot", " ", "the", " ", "solution"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"solution", "[", "x", "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "2"}], ",", "5"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.446759777512793*^9, 3.446759800206653*^9}}], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVl3c81d8fx+3NxVUphNJSFBUavI9IkmRFSaUQDSNlfCUjaWhIZCsjM8nO yieSLTObzM+9rnmvS8j4fX5/vR/Px+P9OI/zPu/X633Okb3uZGzLxsLCoszK wvL/OB7W73DObBk79D66YiebNFJ988PaSmkJK/tb5B1jKokelW7S62H7iyE+ 4cM5c5tQZr6ZzepHJhbLjKuppZKRmmaXYpATA8tPPxEfXCCMohgtW1ZsGdio SFZVW5owWkms57ppycAaN7Avbo0VRhgXNnDiNAMz0yy1b/YXRjoNSa8XtjGw qo4/E++NhJGJ+b3pyx107Mar7bH1/STk6CCUuU+djvX1/86avCqEPkZqK9Tx zWKqEndPCoXxI93BzAZr1lmstz/dX96PH3WdVe5anpnBrAelq87f4Ud29r3n r9XNYCWGfoe6T/CjJ7EKRoo+M1i5Z3r521k+9IO79VQNdRrrnkhlqzXgQxp9 Ww6vlExhopGcRlwbedHBgAySzfUJzCe+s3q+kQsJ7WKv7DGcwMbPipSUfONC 1JqLHkYwgbXo5y28/MyF3gtwD2lITmDN+GNFo1dciCf0Ws7mDhpWu3KQefUs F+pL2Hi+6TQN4+AY9Qtu4kRPvvtGHVUax6rtlj9u6eNA3f9Md4iwUDBV0c7X PmLsaAtTb7vgFI7l3ZEnVXKyI8tJJMvbg2MaUbdkxP6yocG+fVKseQR/hp6O bjaEf+PYQLfDsX5RjsmxD2yI6ZPH0dI0hsV5O0bLKbIhIQ6x0TcfRrEeLZuw eWNWdEKgLUEYhrGK2BuyKx7rsHHzmY2GCsPYhYgHz32ur8PEjh/PgySGsQ2r dcZi+usQCrnOQotDGFUmYviu9DrgLiEaAtlDWGVzEO1f1Rq86Dbp4d4+hHFZ OelpbV6D9uR2kXWuQSyJt6Tpct0K/HqUwdob24tpniysq7JcBl/6lQKBB72Y 0g2e8AK9ZVC2ErmtcaEX6xP/NJWntgzv1N3a40V7MfvPJlv6xZbBchFS7J/2 YF8LRd71NC4BzaFNf8G5G4vWUVX3P7EEXBeXw0W0O7GQP0dizqssAuzXVdSl tWLBY+9krugsQJneEbfHla3Y5ItzXm4qC6BxQ77s+/tWTKLOvCZ2J8GxAgZH TFqxw0qyQpu5CeZvcZAvbcHSgs9fKq+ZB3XKxc8Cr5uxtxwi1i0G8+Cm0Xkh wboBW3FTtc+1ZYKbaxF1i3gDNinpJxtpTnBGtEdIQz0WF7uv5OVpJrhusYp4 fLge+8Lz1zFWgQn3F6idN7jrsMjap0zrhTlwyVw22/upGrMwyLf3DpwDh61b z+cxKrC68KOu2aUMKEzhKTmaUYHdikjREM9mAIfSnEy5bQWW7uSpGZjEgGit msnGrnIMvU/sCXrNgFp7l8eUsu8YR/R6psw1Bpy+8cru8MZv2EeXn4EmPMR6 QQNybKfysGQ3MbH2K3Q4l5A7K1SciwV/FT0fZkoHPO9ZqYRCLmY0d0nYWo8O Yj3KpofJOdhU4M8BCRU6/NEV4m4p+4KREq8d5ROig4HL1crYmWQsfqcTox2b hSwx7pU+7QBs8/gOlX/ysyCvFVDxztsXW96ddUVTdhb0RrhEY53dsTTVIo6Q TbPQfH/tC8sVI6x4bOzzBY5Z8HLnbrac8gXek+W3u0tnoJ2V36A+KBpIpJwb 1ntmwNf6iY6fcxbcsF/tYSVPQ3ZsgYIebzYUCkbFmLBNw3AnLkZOyIb1Nb2c YvoUaOmfGvnYlgPnRjek1DRPAddhHt8alXyY5NxnOhw0BS+4AotIq0Uw+7Hv B6voFNS362xsqiyH6Z+NRVE7JqH1lMV0I6kCtveJ1X3fNAk9xQ5VDRYV0G/q O73EOwnUuFC3utkKePZM9lTS9ARwOAx3/JSshEd1ASx/CifgOKd3eOn9Kjhf 2DXbbjgBmYfzxdPl6uG5UPma92sa5KfWzKY61UOnAb8j4xENSiX6alKK6yH2 8+sId3ca1LGy/5dk2AD23ZXzBVY0oDQadsd5NUJ+6cAXnkM0kLkxGRne3gQu DeyKM/3jcLo++NWJE60Qc1p7xE9jHJIvhQsds20FATvpkeLD48AxGRN08Fkr RDj/UGNTGIfvAqnBcr9agW3v1qoSiXFQMSh7x2XRBuESYkHdS1SQa6HF1t1t hxHP0E1FhVR4dG1264/QdngjZNzalkWFQfr8h5Kv7XDXPcFuJZUK0aKsCRmr 7bA78Du7eyQVRE03Jb9+/hukEo/sF3lAhfUOrUzj+A4wv1Rf14qoYGl3ev+Z yg6ICEEZ1UeoUPzXIEuL0gGDLw/Ta5Wp4CZukXNIoRMMLtk9W5CjwuRF54KN RZ3Q5HVSZoaXCj19MWU9zV2QZ/4tybyDAvnD87+us/QC9nbFuN2DAtpVpaWG O3uha8BrkOFCgfY0/3QN/V5w/CyxXcaBAkxnkSdbInrB8Dv+MfkaBZTX96m3 7u+D+8oHUq+foUCmhM0nzav9UKd3TD5IlgIa6/KR+wP6wSWCzPJOkgKNw/Qn Up/64c4FWba0TRSYTPO5vrTQD5ckTW7+E6TAXrXozdmvB2A89UqP7j8cUk1b n8qU/QGpiSxeny4c1NQiXYVG/4Dj87Ey63YcRgWzlNY4BiGC85qkWTMO9+hu 0pPyg3CypOeLdQ0Obws5/lW7DYLefynmskU4WMhT6WXPB4G1FHNwycdhW0w9 NT9mEA4wFyaas3HI8XnbkVgxCIHHeFgL0nFo0ZHJ8REaAqmryVnSsThEFLKn uckOQZdXgHxvJA5W8pQPDoeG4CP388bEMBxmBTNfXbIYgi82XU6Wb3AQ/n3s pmryEORSwvRvPcbBwNpcZkZ9GH71PP5Md8TB/3thfMK5YSg8F5iWdAeHQqkt 282uDUO7uGOOzS1if119O749Hgbwtmtlv4HDvP71vS/qh6FxZTXu52Uc5NN/ ZGj0D8O2iuKLLZdwuMK9Q5ExPQwN3T3TlIs41JRTDlwUHQHvGZ2ww2Y4xKg4 qOy6OAL2r94pG58j6gv59bXn1gjcc5TclXcWBy76/iOvvUbAxtiqRkYfB+dP jGPzH0bg+PnApc2ncUjiMf2Wmj0CdkqZD5NP4dBjm69h+WMEHJ/yhavr4KAl 46H5Ax8BkWp+lxdaOHg87K5wWxwB/DzP7hMncMjoOaotzzcK/yVNXmLXxGHj u1WdYIVR2K0IO5M0cKDx+utfvT8KoZ5+ibZHcZC2G20UfTIKc1slUpyO4GBa efJcVfgosBaHLvur4VDmzWOsUDIKqfJWnvUqOMz13mwbbBiFLh6z6LXDOOw+ Um8aOjAK9hP35o8RHDz3yvwf6xiwar/V+H0QhyrDma5M8hgcrD23TYng5c+G Ftd3jEGGw5fd4co47OfP6d2gOgaaDzfpcxNsY0++XKs7BhyW9Cd+SjhE/rw/ 4GUxBn6P09o4Cf61rePqgTtjsPHJl32hB3Bg91UdGnk4BkfY3V7uI1itP+J6 eNAY/PjiMvVrPw4OR5dH9OKJ9c9t1XtAcEL4Jdu1nDHQzmXGKBHcwSzFsyvH ILvo/iBdEQd+4632th1jIHhGhlRKMPriMy5OHYPWC17b3xB8X2DoVsPSGFT2 Cos7Epx288SkDz8Omq32lPMED1QlOhyUIvatIfryFMFkOc4ZnIhXdrxn0SKi rt8N5yiEA/YiDnQJfjhQTT9rjMMf1ywtc4Jzju25x2pD6G7MjtuZYEpEIDPP FYe8DuuAtwRLLky42j/FIVBTu7CMYEOTs38lCJ9s2pARxSQ4ICvTo4nwlaeq 5q5DRH3GX52zVr/h0MWljvsTLFOmTN3bgsNdRUX6AMHTlUxpi1FCp3pGmieJ 8yutLzB/9hcHEU5KUwHBga0eQQV8FPAQPhR2kDj/C91Hq0elKLCyy+9VCcFM vExFU5sCt32MmZNE/yqmfB2dzClwRMfLLIzo9xvmieTYWxSQ3Eb/pXcIh71s 1RuW31DgBc9Bajuhl0We52d3f6RAifvPB+mEvqpIZwLMvlLg2XcF/heqOFyX amLm9FNASVaF3YHQ5wG54H1DsxSw219yxZHQ75q8iQ2Jgwrq63+iPY/hEKXW 2XZbngrpQ6txueo42EMkf5Q6FdYO/Ls0SOhfReeSVo0hFb4XxA1sIvrQYvIn R86dCrefb5JNIfzz3iKeZhxIhSYhgfY1wl93rllv84ulwpLKzKFrJ3HgcaIE 9/+gwvYdDPykLqGPwBnnMOFx2H3m98XXBjgIBWenVm4fh7hGobqjhjj0hd8b ZKiMA7vNg3y6EQ7uSX/Pnbs8DqzTnaPe53HI/L6uyJM+Dg3anQsbiHnjVV1+ Q+XbOHDWVgdvuIqD3i//9zbN41DRJuMmfQ0HvJdbqHxhHGxrHdiMbXGQ+kua +k+LBitnaA+ViHlHW22RSzWjgYiCV+EjZ2LecYRadtykwe+gCP4BF8LfouIN ym9owCl5aiTHHYeXCjKfJvpocFhHNqvdj/Cn9YGbl90m4HJoybO5cBy0kchR 5ecTIHxb7+vvKBxeSTL4uGMmYK44/mQ5Ma9lfudmZJVPwP73CbfzEnHQ0VGd ZRechD2Gzlp8WcT9sFvDPT1xEpL6uqyc6ojz4JDW9SmYhNQibHW0EYcdQ+vi prVEPntVojWhy8LI8uLV6Um4bT4//4C4jwb4Tq4ZHp2CxbbJreIUwhdTZwL+ Eu8myY7WqSVOCtyv3WfWMDIFWz7IJK/yUuBbkuCu+AUif4JTjE+IAoZXmmr0 JKfhWIf0yvGNFHBrMhZ4bzcNAs4Ljny7KFCRfTFEe20aaDrb7oMeBRRKcptr NszCriVjjZkwCvBIW3sHb50FM2XH/NRoCow+Et1nsWsW3nonqN6Mo0CU3t2n E2qz8Flq7Ox6GpHfo6ghZDkLRfcOmPl8o8DwYtonk4RZWG3L8akfo0CYStyT AUU6mHNV0v87RoW7UecOpajS4chUkIegJhX019aGnBAdWIcH3qboUIG96rI6 qzEdXBlR92aMqOB8XoK53ZUOedKCfo32VNC7H3btZgkd9qSZZGLhxHsm++Vx 5ikGiJf/rLdepkLyrTsnPY0YoOu877U+6zgYbNc3YLnEgCsHFp8e5xmH96H8 VgKODJDYHCGsvHEcNNwD/eVCGTD8zGPNRXkcvI4/qzMdZMCtqNUI3TvjsFjl fyHfYw5atf983jE2DrO9nvfd0pjAnpLqUEWjgUFgZIVFDhNCL1vV1TFokKFW KAwlTPipIRrSvkyDm+/mM7gbmRAgaWm5yDcBwwbOYxGzTJgdSff03TsBv8tt zUpU5yHteqH+G4cJKE4xVGOpmoeXaq7NkfMT4H9v58rz4QXivg3ozxebgpfF GSnOwksw5BHzzPbhDOQe0pXOdPgHA/9C1NiU6PBQMh73SliF4l1cNjUpDMiJ 5lY3YmdByRuXq/0uMmGIJSOicT8rmqx7ne0gsQBDkvtjRHXZ0HZ11uWvLX/B ubWUk92DHX1Ym9GLL1uCZo09m0pDONABPSdLych/MMnmXzyQxomUsn9bT7uu QrgV5UdGJxca5Mi+k/tkHciu8q5adG7E46awklbGgrZu93gZsI0XmXK8zPRs Y0XLTz+eeQN8aOCn68OnxP85WLsshsOcH7nX8Rhd6GdHX6OGTeSeCKAP3l3C q0sc6EvKWI5AvCD6s5EhPMfJhVQM5tzj2oRQbs+eBB95bnR8pzvFgEUY/Xwi fyX2DA+qYhnFs18KI5E7NgnVFrxIfLdZJE1OBJE+vio48IgPtaZ//cCdK4Ko kWmHpML50a1uktNFDVFU2dJh458igGpNHei+v0WR63+5mowmQRS8t7j65XUy Cpazlz09JoRax++HstmSkRDvmXkZmhASTd1/zcOOjJCmye/FaSH0Vi5p2foO GQl05xelLgqhUKk3CsfcyKinb06UV4CEwkg3QsafkxGpYYSlQZmEouZErupk k1FZLluEhS8J9WQ37C3NJaMv4S7DBwNIaIvz00WlAjLaWiSrJhhIQtGTK8FS JWTUsWEb5/cQEooZo1QyK8noYiNNbmcKCX3o/Caf2EVGIbnJrsxGEhp85/FX vJeMPIM6nH+1kpCM6cHK1/1kxLGh2y21k4TimlMvew6TEevzY4mWQySUUBvy xmiSjEYnS71+Mklo+KnB5eppMrJqqWn8sERC23R45dXpZLQ73mWP5xoJJVZ4 /9i9QNSzeMJ3E4cwGvU98ub9IlG/Gv3bYy5hJAdMS7F/ZOToOMTC4BFGNquZ ewJXycjWSvj0VX5h9LHk5sL6OhlJ/PYMaxAURv8DJCekKw== "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 1.}, PlotRange->{{-2, 5}, {0.693147183615263, 7.752100983379524}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Scaled[0.02]}]], "Output", CellChangeTimes->{{3.4467597962449923`*^9, 3.4467598013706408`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "ONE", " ", "OF", " ", "THE", " ", "NICE", " ", "NEW", " ", "FEATURES", " ", "IN", " ", "VERSIONS", " ", "6", " ", "AND", " ", "7", " ", "IS", " ", RowBox[{"MANIPULATE", ".", " ", "HERE"}], " ", "WE", " ", "SHOW", " ", "HOW", " ", "WE", " ", "CAN", " ", "MANIPULATE", " ", "THE", " ", "PLOT", " ", "OF", " ", "THE", " ", "SOLUTION", " ", "WITH", " ", "TWO", " ", "FREE", " ", "PARAMETERS"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "IMAGINE", " ", "THE", " ", "SOLUTION", " ", "TO", " ", "A", " ", "DIFFERENTIAL", " ", "EQUATION", " ", "IS", " ", "THE", " ", "FOLLOWING"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"solution", "[", RowBox[{"x_", ",", "c1_", ",", "c2_"}], "]"}], " ", ":=", " ", RowBox[{ RowBox[{"c1", " ", RowBox[{"(", RowBox[{ RowBox[{"x", "^", "2"}], " ", "-", " ", RowBox[{"2", "x"}], " ", "+", " ", "1"}], ")"}]}], " ", "+", " ", RowBox[{"c2", " ", RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "x"}], "/", "5"}], "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"solution", "[", RowBox[{"x", ",", "c1", ",", "c2"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", RowBox[{"xcenter", "-", "xrange"}], ",", RowBox[{"xcenter", "+", "xrange"}]}], "}"}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"xcenter", ",", RowBox[{"-", "100"}], ",", "100"}], "}"}], ",", RowBox[{"{", RowBox[{"xrange", ",", ".01", ",", "100"}], "}"}], ",", RowBox[{"{", RowBox[{"c1", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"c2", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.446759908404806*^9, 3.4467600622767935`*^9}, { 3.446760100533134*^9, 3.4467601906621313`*^9}, {3.4467602437827544`*^9, 3.4467602578756027`*^9}, {3.446760292425149*^9, 3.4467603328122125`*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`c1$$ = 0.6600000000000001, $CellContext`c2$$ = 0.9199999999999999, $CellContext`xcenter$$ = 1.5, $CellContext`xrange$$ = 3.1, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`xcenter$$], -100, 100}, { Hold[$CellContext`xrange$$], 0.01, 100}, { Hold[$CellContext`c1$$], -2, 2}, { Hold[$CellContext`c2$$], -2, 2}}, Typeset`size$$ = {360., {115., 118.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`xcenter$108664$$ = 0, $CellContext`xrange$108665$$ = 0, $CellContext`c1$108666$$ = 0, $CellContext`c2$108667$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`c1$$ = -2, $CellContext`c2$$ = -2, \ $CellContext`xcenter$$ = -100, $CellContext`xrange$$ = 0.01}, "ControllerVariables" :> { Hold[$CellContext`xcenter$$, $CellContext`xcenter$108664$$, 0], Hold[$CellContext`xrange$$, $CellContext`xrange$108665$$, 0], Hold[$CellContext`c1$$, $CellContext`c1$108666$$, 0], Hold[$CellContext`c2$$, $CellContext`c2$108667$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Plot[ $CellContext`solution[$CellContext`x, $CellContext`c1$$, \ $CellContext`c2$$], {$CellContext`x, $CellContext`xcenter$$ - \ $CellContext`xrange$$, $CellContext`xcenter$$ + $CellContext`xrange$$}], "Specifications" :> {{$CellContext`xcenter$$, -100, 100}, {$CellContext`xrange$$, 0.01, 100}, {$CellContext`c1$$, -2, 2}, {$CellContext`c2$$, -2, 2}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{409., {201., 206.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.4467600197629128`*^9, 3.446760062986709*^9}, { 3.4467601476189876`*^9, 3.446760191328104*^9}, 3.4467602638888893`*^9, 3.446760297012746*^9, 3.4467603336500497`*^9}] }, Open ]], Cell[BoxData[""], "Input", CellChangeTimes->{{3.446759855226569*^9, 3.446759855642558*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "HERE", " ", "WE", " ", "SOLVE", " ", "A", " ", "SECOND", " ", "ORDER", " ", "DIFFERENTIAL", " ", "EQUATION"}], ",", " ", RowBox[{ "FIRST", " ", "WITH", " ", "AND", " ", "THEN", " ", "WITHOUT", " ", "INITIAL", " ", RowBox[{"CONDITIONS", ".", " ", "NOTE"}], " ", "THE", " ", RowBox[{"SYNTAX", " ", "--"}], " ", "WE", " ", "HAVE", " ", "THE", " ", "DIFFERENTIAL", " ", "EQUATION", " ", "AND", " ", "THE", " ", "INITIAL", " ", "CONDITIONS", " ", "IN", " ", RowBox[{"{", "}"}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"AND", " ", "WE", " ", "USE"}], " ", "\[Equal]", " ", RowBox[{"AND", " ", "NOT"}]}], " ", "=", " ", RowBox[{"ALWAYS", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", "y", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "''"}], "[", "x", "]"}], " ", "+", " ", RowBox[{"y", "[", "x", "]"}]}], " ", "==", " ", "0"}], ",", " ", RowBox[{"y", "[", "x", "]"}], ",", " ", "x"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Clear", "[", "y", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "''"}], "[", "x", "]"}], " ", "+", " ", RowBox[{"y", "[", "x", "]"}]}], " ", "==", " ", "0"}], ",", " ", RowBox[{ RowBox[{"y", "[", "0", "]"}], " ", "==", " ", "2"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "0", "]"}], " ", "==", " ", "1"}]}], "}"}], ",", " ", RowBox[{"y", "[", "x", "]"}], ",", " ", "x"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.446680122062327*^9, 3.4466801592347527`*^9}, { 3.446680195502939*^9, 3.446680207010015*^9}, {3.4466802466910686`*^9, 3.4466802650502453`*^9}, {3.4466803217263055`*^9, 3.4466803719076643`*^9}, { 3.4467158449719687`*^9, 3.446715846237513*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], "\[Rule]", RowBox[{ RowBox[{ RowBox[{"C", "[", "1", "]"}], " ", RowBox[{"Cos", "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"C", "[", "2", "]"}], " ", RowBox[{"Sin", "[", "x", "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.4466801240787992`*^9, 3.44668015977671*^9}, { 3.446680198432439*^9, 3.4466802079513607`*^9}, {3.44668024787947*^9, 3.446680265441822*^9}, {3.44668033054703*^9, 3.446680372343183*^9}, { 3.446680406029504*^9, 3.4466804120047073`*^9}, 3.4467158508631935`*^9, 3.4467162537309628`*^9, 3.4467205118947277`*^9, 3.4467603852942047`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "x", "]"}], "\[Rule]", RowBox[{ RowBox[{"2", " ", RowBox[{"Cos", "[", "x", "]"}]}], "+", RowBox[{"Sin", "[", "x", "]"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.4466801240787992`*^9, 3.44668015977671*^9}, { 3.446680198432439*^9, 3.4466802079513607`*^9}, {3.44668024787947*^9, 3.446680265441822*^9}, {3.44668033054703*^9, 3.446680372343183*^9}, { 3.446680406029504*^9, 3.4466804120047073`*^9}, 3.4467158508631935`*^9, 3.4467162537309628`*^9, 3.4467205118947277`*^9, 3.446760385313735*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "ONE", " ", "OF", " ", "THE", " ", "NICE", " ", "NEW", " ", "FEATURES", " ", "IN", " ", "VERSIONS", " ", "6", " ", "AND", " ", "7", " ", "IS", " ", RowBox[{"MANIPULATE", ".", " ", "HERE"}], " ", "IS", " ", "SOME", " ", "CODE", " ", "TO", " ", "SHOW", " ", "HOW", " ", "ADDITIONAL", " ", "TERMS", " ", "IN", " ", "THE", " ", "TAYLOR", " ", "SERIES", " ", "EXPANSION", " ", "OF", " ", "SIN", RowBox[{"(", "X", ")"}], " ", "GIVE", " ", "A", " ", "BETTER", " ", "AND", " ", "BETTER", " ", RowBox[{"FIT", "."}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ "Note", " ", "that", " ", "once", " ", "we", " ", "have", " ", "51", " ", "or", " ", "so", " ", "terms"}], ",", " ", RowBox[{ "Mathematica", " ", "is", " ", "no", " ", "longer", " ", "evaluating", " ", "things", " ", RowBox[{"properly", ".", " ", "The"}], " ", "reason", " ", "is", " ", "we", " ", "are", " ", "now", " ", "raising", " ", "large", " ", "numbers", " ", "to", " ", "large", " ", "powers", " ", "and", " ", "then", " ", "dividing", " ", "by", " ", "large", " ", RowBox[{"factorials", " ", "--"}], " ", "if", " ", "we", " ", RowBox[{"don", "'"}], "t", " ", "group", " ", "the", " ", "terms", " ", "carefully"}], ",", " ", RowBox[{ "there", " ", "will", " ", "easily", " ", "be", " ", "overflow", " ", RowBox[{"problems", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Sin", "[", "x", "]"}], ",", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"-", "1."}], ")"}], "^", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}], " ", RowBox[{ RowBox[{"x", "^", RowBox[{"(", RowBox[{ RowBox[{"2", "n"}], "-", "1"}], ")"}]}], " ", "/", " ", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "n"}], " ", "-", " ", "1"}], ")"}], "!"}], ")"}]}]}], ",", RowBox[{"{", RowBox[{"n", ",", "1", ",", RowBox[{"Floor", "[", "num", "]"}]}], "}"}]}], "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", "0", ",", RowBox[{"xrange", " ", "2", " ", "Pi"}]}], "}"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "yrange"}], ",", "yrange"}], "}"}]}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"num", ",", "1", ",", "80"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"xrange", ",", " ", "1", ",", " ", "10"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"yrange", ",", "1.25", ",", "20"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4466804340745835`*^9, 3.446680578113216*^9}, { 3.446680618074526*^9, 3.4466806211563597`*^9}, {3.4466807212876463`*^9, 3.4466808005960464`*^9}, {3.4466808332892666`*^9, 3.4466809479674735`*^9}, {3.446715857294422*^9, 3.4467158887064743`*^9}, 3.4467205469696317`*^9, {3.4467205790398445`*^9, 3.446720626938146*^9}, { 3.4467206622688923`*^9, 3.44672067140112*^9}, {3.4467604056356764`*^9, 3.4467604417359047`*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`num$$ = 1, $CellContext`xrange$$ = 1, $CellContext`yrange$$ = 1.25, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`num$$], 1, 80}, { Hold[$CellContext`xrange$$], 1, 10}, { Hold[$CellContext`yrange$$], 1.25, 20}}, Typeset`size$$ = { 360., {105., 108.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`num$113281$$ = 0, $CellContext`xrange$113282$$ = 0, $CellContext`yrange$113283$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`num$$ = 1, $CellContext`xrange$$ = 1, $CellContext`yrange$$ = 1.25}, "ControllerVariables" :> { Hold[$CellContext`num$$, $CellContext`num$113281$$, 0], Hold[$CellContext`xrange$$, $CellContext`xrange$113282$$, 0], Hold[$CellContext`yrange$$, $CellContext`yrange$113283$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Plot[{ Sin[$CellContext`x], Sum[(-1.)^($CellContext`n + 1) ($CellContext`x^(2 $CellContext`n - 1)/Factorial[ 2 $CellContext`n - 1]), {$CellContext`n, 1, Floor[$CellContext`num$$]}]}, {$CellContext`x, 0, $CellContext`xrange$$ 2 Pi}, PlotRange -> {-$CellContext`yrange$$, $CellContext`yrange$$}], "Specifications" :> {{$CellContext`num$$, 1, 80}, {$CellContext`xrange$$, 1, 10}, {$CellContext`yrange$$, 1.25, 20}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{409., {177., 182.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{{3.446680550146256*^9, 3.4466805787723536`*^9}, 3.446680621989314*^9, {3.4466807217378125`*^9, 3.446680801567664*^9}, { 3.446680841813135*^9, 3.4466808521650114`*^9}, {3.4466808964024143`*^9, 3.446680948482089*^9}, 3.446715889356824*^9, 3.446716255828485*^9, 3.446720513089964*^9, 3.4467205479246483`*^9, 3.4467206278999987`*^9, { 3.4467206634348335`*^9, 3.4467206719235477`*^9}, {3.4467603885078664`*^9, 3.4467604424536324`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "HERE", " ", "IS", " ", "AN", " ", "EXAMPLE", " ", "OF", " ", "A", " ", "SYSTEM", " ", "OF", " ", "DIFFERENTIAL", " ", "EQUATIONS"}], ",", " ", RowBox[{"WITH", " ", "AND", " ", "WITHOUT", " ", "INITIAL", " ", RowBox[{"CONDITIONS", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", "x", "]"}], ";", " ", RowBox[{"Clear", "[", "y", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "t", "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"2", " ", RowBox[{"x", "[", "t", "]"}]}], " ", "+", " ", RowBox[{"3", " ", RowBox[{"y", "[", "t", "]"}]}]}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "t", "]"}], " ", "==", " ", RowBox[{ RowBox[{"3", " ", RowBox[{"x", "[", "t", "]"}]}], " ", "+", " ", RowBox[{"2", " ", RowBox[{"y", "[", "t", "]"}]}]}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"y", ",", " ", "x"}], "}"}], ",", " ", "t"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Clear", "[", "x", "]"}], ";", " ", RowBox[{"Clear", "[", "y", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "t", "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"2", " ", RowBox[{"x", "[", "t", "]"}]}], " ", "+", " ", RowBox[{"3", " ", RowBox[{"y", "[", "t", "]"}]}]}]}], ",", " ", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "t", "]"}], " ", "==", " ", RowBox[{ RowBox[{"3", " ", RowBox[{"x", "[", "t", "]"}]}], " ", "+", " ", RowBox[{"2", " ", RowBox[{"y", "[", "t", "]"}]}]}]}], ",", " ", RowBox[{ RowBox[{"x", "[", "0", "]"}], " ", "\[Equal]", " ", "4"}], ",", " ", RowBox[{ RowBox[{"y", "[", "0", "]"}], " ", "\[Equal]", " ", "2"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"y", ",", " ", "x"}], "}"}], ",", " ", "t"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.4467159192181935`*^9, 3.446716056808996*^9}, { 3.446716097972378*^9, 3.4467161005513144`*^9}, {3.446716151152568*^9, 3.4467161600152817`*^9}, {3.4467161997773848`*^9, 3.446716201431576*^9}, 3.4467162886554856`*^9, {3.4467164134072666`*^9, 3.446716461711792*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"x", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", "t", "}"}], ",", RowBox[{ RowBox[{ FractionBox["1", "2"], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], " ", RowBox[{"(", RowBox[{"1", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"6", " ", "t"}]]}], ")"}], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{ FractionBox["1", "2"], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"6", " ", "t"}]]}], ")"}], " ", RowBox[{"C", "[", "2", "]"}]}]}]}], "]"}]}], ",", RowBox[{"y", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", "t", "}"}], ",", RowBox[{ RowBox[{ FractionBox["1", "2"], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"6", " ", "t"}]]}], ")"}], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{ FractionBox["1", "2"], " ", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], " ", RowBox[{"(", RowBox[{"1", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"6", " ", "t"}]]}], ")"}], " ", RowBox[{"C", "[", "2", "]"}]}]}]}], "]"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446716154702145*^9, 3.4467161607789044`*^9}, 3.4467162021707864`*^9, 3.446716257859605*^9, {3.446716289383954*^9, 3.4467163002504463`*^9}, 3.446716334182845*^9, 3.4467164315047407`*^9, 3.4467164631765423`*^9, 3.4467604442874994`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"x", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", "t", "}"}], ",", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"6", " ", "t"}]]}]}], ")"}]}]}], "]"}]}], ",", RowBox[{"y", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", "t", "}"}], ",", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"3", " ", SuperscriptBox["\[ExponentialE]", RowBox[{"6", " ", "t"}]]}]}], ")"}]}]}], "]"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446716154702145*^9, 3.4467161607789044`*^9}, 3.4467162021707864`*^9, 3.446716257859605*^9, {3.446716289383954*^9, 3.4467163002504463`*^9}, 3.446716334182845*^9, 3.4467164315047407`*^9, 3.4467164631765423`*^9, 3.446760444309959*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "THIS", " ", "EXAMPLE", " ", "IS", " ", "TYPICAL", " ", "OF", " ", "WHAT", " ", "HAPPENS"}], ",", " ", RowBox[{ "NAMELY", " ", "THAT", " ", "MOST", " ", "OF", " ", "THESE", " ", "EQUATIONS", " ", "ARE", " ", "TOO", " ", "HARD", " ", "TO", " ", "BE", " ", "SOLVED", " ", RowBox[{"EXACTLY", ".", " ", "HERE"}], " ", "IS", " ", "AN", " ", "EXAMPLE"}], ",", " ", RowBox[{"WHERE", " ", RowBox[{"K", "[", "1", "]"}], " ", "IS", " ", "A", " ", "DUMMY", " ", "VARIABLE", " ", "OF", " ", RowBox[{"INTEGRATION", ".", " ", "IF"}], " ", "WE", " ", "CHANTE", " ", RowBox[{"t", "^", "2"}], " ", "TO", " ", "t", " ", "THEN", " ", "WE", " ", "CAN", " ", "SOLVE", " ", "IT", " ", "IN", " ", "TERMS", " ", "OF", " ", "SPECIAL", " ", "FUNCTIONS"}], ",", " ", RowBox[{ RowBox[{ "OR", " ", "IF", " ", "WE", " ", "CHANGE", " ", "THE", " ", "NON"}], "-", RowBox[{"HOMOGENOUS", " ", RowBox[{"TERM", "."}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "t", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t", "^", "2"}], " ", RowBox[{"y", "[", "t", "]"}]}]}], " ", "\[Equal]", " ", RowBox[{"Exp", "[", "t", "]"}]}], ",", " ", RowBox[{"y", "[", "t", "]"}], ",", " ", "t"}], "]"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "t", "]"}], " ", "+", " ", RowBox[{"t", " ", RowBox[{"y", "[", "t", "]"}]}]}], " ", "\[Equal]", " ", RowBox[{"Exp", "[", "t", "]"}]}], ",", " ", RowBox[{"y", "[", "t", "]"}], ",", " ", "t"}], "]"}], "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "t", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t", "^", "2"}], " ", RowBox[{"y", "[", "t", "]"}]}]}], " ", "\[Equal]", " ", "1"}], ",", " ", RowBox[{"y", "[", "t", "]"}], ",", " ", "t"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.4467164736719637`*^9, 3.446716561869444*^9}, { 3.4467165941935472`*^9, 3.446716664978079*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ SuperscriptBox["t", "3"], "3"]}]], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ SuperscriptBox["t", "3"], "3"]}]], " ", RowBox[{ SubsuperscriptBox["\[Integral]", "1", "t"], RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"K", "[", "1", "]"}], "+", FractionBox[ SuperscriptBox[ RowBox[{"K", "[", "1", "]"}], "3"], "3"]}]], RowBox[{"\[DifferentialD]", RowBox[{"K", "[", "1", "]"}]}]}]}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{ 3.446716511282838*^9, {3.4467165955206103`*^9, 3.446716602282873*^9}, 3.4467166661381607`*^9, 3.44676044737031*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ SuperscriptBox["t", "2"], "2"]}]], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", FractionBox["1", "2"]}], "-", FractionBox[ SuperscriptBox["t", "2"], "2"]}]], " ", SqrtBox[ FractionBox["\[Pi]", "2"]], " ", RowBox[{"Erfi", "[", FractionBox[ RowBox[{"1", "+", "t"}], SqrtBox["2"]], "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{ 3.446716511282838*^9, {3.4467165955206103`*^9, 3.446716602282873*^9}, 3.4467166661381607`*^9, 3.446760447395699*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"y", "[", "t", "]"}], "\[Rule]", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ SuperscriptBox["t", "3"], "3"]}]], " ", RowBox[{"C", "[", "1", "]"}]}], "-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", FractionBox[ SuperscriptBox["t", "3"], "3"]}]], " ", "t", " ", RowBox[{"Gamma", "[", RowBox[{ FractionBox["1", "3"], ",", RowBox[{"-", FractionBox[ SuperscriptBox["t", "3"], "3"]}]}], "]"}]}], RowBox[{ SuperscriptBox["3", RowBox[{"2", "/", "3"}]], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", SuperscriptBox["t", "3"]}], ")"}], RowBox[{"1", "/", "3"}]]}]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{ 3.446716511282838*^9, {3.4467165955206103`*^9, 3.446716602282873*^9}, 3.4467166661381607`*^9, 3.446760447439642*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "WHEN", " ", "WE", " ", "DO", " ", "GET", " ", "SOLUTIONS", " ", "LIKE", " ", "THE", " ", "ABOVE"}], ",", " ", RowBox[{ "IT", " ", "CAN", " ", "INVOLVE", " ", "FUNCTIONS", " ", "YOU", " ", RowBox[{"HAVEN", "'"}], "T", " ", "SEEN"}], ",", " ", RowBox[{"SUCH", " ", "AS", " ", RowBox[{"ERFI", ".", " ", "TO"}], " ", "FIND", " ", "OUT", " ", "MORE", " ", "ABOUT", " ", "THAT"}], ",", " ", RowBox[{ "EITHER", " ", "GO", " ", "TO", " ", "FIND", " ", "SELECTED", " ", "FUNCTION"}], ",", " ", RowBox[{ RowBox[{ "OR", " ", "TYPE", " ", "QUESTION", " ", "MARK", " ", "ERFI", " ", "AT", " ", "THE", " ", "NEXT", " ", "COMMAND", " ", RowBox[{"LINE", ".", " ", "YOU"}], " ", "CAN", " ", "THEN", " ", "CLICK", " ", "ON", " ", "THE"}], " ", ">>", " ", RowBox[{"FOR", " ", "MORE", " ", RowBox[{"INFO", "."}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"?", " ", "Erfi"}]}]], "Input", CellChangeTimes->{{3.4467166728652697`*^9, 3.44671673709651*^9}}], Cell[BoxData[ RowBox[{ StyleBox["\<\"\!\(\*RowBox[{\\\"Erfi\\\", \\\"[\\\", StyleBox[\\\"z\\\", \\\ \"TI\\\"], \\\"]\\\"}]\) gives the imaginary error function \ \!\(\*RowBox[{RowBox[{\\\"erf\\\", \\\"(\\\", RowBox[{\\\"i\\\", \\\"\ \[InvisibleSpace]\\\", \\\"z\\\"}], \\\")\\\"}], \\\"/\\\", \\\"i\\\"}]\). \"\ \>", "MSG"], "\[NonBreakingSpace]", ButtonBox[ StyleBox["\[RightSkeleton]", "SR"], Active->True, BaseStyle->"Link", ButtonData->"paclet:ref/Erfi"]}]], "Print", "PrintUsage", CellChangeTimes->{3.446760451431574*^9}, CellTags->"Info3446746051-8245676"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "You", " ", "can", " ", "also", " ", "use", " ", "Mathematica", " ", "to", " ", "solve", " ", "difference", " ", "equations"}], ",", " ", RowBox[{"which", " ", "it", " ", "calls", " ", "recurrence", " ", RowBox[{"equations", "."}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "Does", " ", "the", " ", "last", " ", "solution", " ", "mean", " ", "anything", " ", "to", " ", "you", RowBox[{"??", "?"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"RSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "2"}], "]"}], " ", "-", " ", RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}], " ", "-", " ", RowBox[{"a", "[", "n", "]"}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}], "\[IndentingNewLine]", RowBox[{"RSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "2"}], "]"}], " ", "-", " ", RowBox[{"3", " ", RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}]}], " ", "-", " ", RowBox[{"10", " ", RowBox[{"a", "[", "n", "]"}]}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}], "\[IndentingNewLine]", RowBox[{"RSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "2"}], "]"}], " ", "-", " ", RowBox[{"3", " ", RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}]}], " ", "-", " ", RowBox[{"10", " ", RowBox[{"a", "[", "n", "]"}]}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{"a", "[", "0", "]"}], " ", "\[Equal]", " ", "1"}], ",", " ", RowBox[{ RowBox[{"a", "[", "1", "]"}], " ", "\[Equal]", " ", "2"}]}], "}"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}], "\[IndentingNewLine]", RowBox[{"RSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "2"}], "]"}], " ", "-", " ", RowBox[{"6", " ", RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}]}], " ", "+", " ", RowBox[{"9", " ", RowBox[{"a", "[", "n", "]"}]}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}], "\[IndentingNewLine]", RowBox[{"RSolve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "2"}], "]"}], " ", "-", " ", RowBox[{"6", " ", "n", " ", RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}]}], " ", "+", " ", RowBox[{"6", " ", RowBox[{"a", "[", "n", "]"}]}]}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.446720176705244*^9, 3.4467203587551155`*^9}, { 3.446729976690366*^9, 3.4467300023693857`*^9}, {3.446730042866794*^9, 3.44673008695284*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{ RowBox[{ RowBox[{"C", "[", "1", "]"}], " ", RowBox[{"Fibonacci", "[", "n", "]"}]}], "+", RowBox[{ RowBox[{"C", "[", "2", "]"}], " ", RowBox[{"LucasL", "[", "n", "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446720283301937*^9, 3.4467203595714693`*^9}, 3.4467298016859417`*^9, 3.446730006718717*^9, {3.4467300543885174`*^9, 3.446730088242796*^9}, 3.446760455741844*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "2"}], ")"}], "n"], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["5", "n"], " ", RowBox[{"C", "[", "2", "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446720283301937*^9, 3.4467203595714693`*^9}, 3.4467298016859417`*^9, 3.446730006718717*^9, {3.4467300543885174`*^9, 3.446730088242796*^9}, 3.4467604557652807`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{ FractionBox["1", "7"], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "2"}], ")"}], "n"]}], "+", RowBox[{"4", " ", SuperscriptBox["5", "n"]}]}], ")"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446720283301937*^9, 3.4467203595714693`*^9}, 3.4467298016859417`*^9, 3.446730006718717*^9, {3.4467300543885174`*^9, 3.446730088242796*^9}, 3.4467604557691865`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{ RowBox[{ SuperscriptBox["3", "n"], " ", RowBox[{"C", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["3", "n"], " ", "n", " ", RowBox[{"C", "[", "2", "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446720283301937*^9, 3.4467203595714693`*^9}, 3.4467298016859417`*^9, 3.446730006718717*^9, {3.4467300543885174`*^9, 3.446730088242796*^9}, 3.446760455797505*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{ RowBox[{ RowBox[{"C", "[", "1", "]"}], " ", RowBox[{"Hypergeometric0F1Regularized", "[", RowBox[{"n", ",", RowBox[{"-", FractionBox["1", "6"]}]}], "]"}]}], "+", RowBox[{ RowBox[{"C", "[", "2", "]"}], " ", RowBox[{"Hypergeometric0F1Regularized", "[", RowBox[{"n", ",", RowBox[{"-", FractionBox["1", "6"]}]}], "]"}], " ", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{ RowBox[{"K", "[", "2", "]"}], "=", "0"}], RowBox[{ RowBox[{"-", "1"}], "+", "n"}]], RowBox[{ UnderoverscriptBox["\[Product]", RowBox[{ RowBox[{"K", "[", "1", "]"}], "=", "1"}], RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"K", "[", "2", "]"}]}]], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"Hypergeometric0F1Regularized", "[", RowBox[{ RowBox[{"2", "+", RowBox[{"K", "[", "1", "]"}]}], ",", RowBox[{"-", FractionBox["1", "6"]}]}], "]"}]}], "+", RowBox[{"6", " ", RowBox[{"Hypergeometric0F1Regularized", "[", RowBox[{ RowBox[{"1", "+", RowBox[{"K", "[", "1", "]"}]}], ",", RowBox[{"-", FractionBox["1", "6"]}]}], "]"}], " ", RowBox[{"K", "[", "1", "]"}]}]}], ")"}], "/", RowBox[{"Hypergeometric0F1Regularized", "[", RowBox[{ RowBox[{"2", "+", RowBox[{"K", "[", "1", "]"}]}], ",", RowBox[{"-", FractionBox["1", "6"]}]}], "]"}]}]}]}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.446720283301937*^9, 3.4467203595714693`*^9}, 3.4467298016859417`*^9, 3.446730006718717*^9, {3.4467300543885174`*^9, 3.446730088242796*^9}, 3.4467604561324444`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"?", " ", "Hypergeometric0F1Regularized"}]], "Input", CellChangeTimes->{{3.4467203622812567`*^9, 3.4467203706352143`*^9}}], Cell[BoxData[ RowBox[{ StyleBox["\<\"\!\(\*RowBox[{\\\"Hypergeometric0F1Regularized\\\", \ \\\"[\\\", RowBox[{StyleBox[\\\"a\\\", \\\"TI\\\"], \\\",\\\", StyleBox[\\\"z\ \\\", \\\"TI\\\"]}], \\\"]\\\"}]\) is the regularized confluent \ hypergeometric function \!\(\*RowBox[{RowBox[{SubscriptBox[\\\"\[Null]\\\", \ \\\"0\\\"], \\\"\[InvisibleSpace]\\\", SubscriptBox[\\\"F\\\", \\\"1\\\"], \\\ \"(\\\", RowBox[{\\\"a\\\", \\\";\\\", \\\"z\\\"}], \\\")\\\"}], \\\"/\\\", \ RowBox[{\\\"\[CapitalGamma]\\\", \\\"(\\\", \\\"a\\\", \\\")\\\"}]}]\). \"\>", "MSG"], "\[NonBreakingSpace]", ButtonBox[ StyleBox["\[RightSkeleton]", "SR"], Active->True, BaseStyle->"Link", ButtonData->"paclet:ref/Hypergeometric0F1Regularized"]}]], "Print", \ "PrintUsage", CellChangeTimes->{3.446760458214342*^9}, CellTags->"Info3446746058-4040371"] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "The", " ", "following", " ", "examples", " ", "illustrate", " ", "how", " ", "powerful", " ", "the", " ", "symbolic", " ", "algorithms", " ", "in", " ", "Mathematica", " ", RowBox[{"are", ".", " ", "We"}], " ", "proved", " ", "in", " ", "one", " ", "of", " ", "the", " ", "Green", " ", "Chicken", " ", "problem", " ", "sessions", " ", "that", " ", "the", " ", "sequence", " ", "below", " ", "has", " ", "limit", " ", "equal", " ", "to", " ", "log", RowBox[{"(", "2", ")"}]}], ";", " ", RowBox[{"Mathematica", " ", "is", " ", "able", " ", "to", " ", "see", " ", RowBox[{"this", "!"}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.446759511611843*^9, 3.446759553882575*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"RSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}], " ", "==", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", "n", "]"}], " ", "+", " ", RowBox[{"n", " ", RowBox[{"a", "[", RowBox[{"n", "-", "1"}], "]"}]}]}], ")"}], " ", "/", " ", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}]}], ",", " ", RowBox[{ RowBox[{"a", "[", "0", "]"}], " ", "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{"a", "[", "1", "]"}], " ", "\[Equal]", "1"}]}], "}"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}]], "Input", CellChangeTimes->{{3.4467298440074515`*^9, 3.4467299117306557`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{ RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"]}], " ", RowBox[{"LerchPhi", "[", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", RowBox[{"1", "+", "n"}]}], "]"}]}], "+", RowBox[{"Log", "[", "2", "]"}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4467299138838387`*^9, 3.4467604609758844`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"?", " ", "LerchPhi"}]], "Input"], Cell[BoxData[ RowBox[{ StyleBox["\<\"\!\(\*RowBox[{\\\"LerchPhi\\\", \\\"[\\\", \ RowBox[{StyleBox[\\\"z\\\", \\\"TI\\\"], \\\",\\\", StyleBox[\\\"s\\\", \ \\\"TI\\\"], \\\",\\\", StyleBox[\\\"a\\\", \\\"TI\\\"]}], \\\"]\\\"}]\) \ gives the Lerch transcendent \!\(\*RowBox[{StyleBox[\\\"\[CapitalPhi]\\\", \\\ \"TR\\\"], RowBox[{\\\"(\\\", RowBox[{StyleBox[\\\"z\\\", \\\"TI\\\"], \ \\\",\\\", StyleBox[\\\"s\\\", \\\"TI\\\"], \\\",\\\", StyleBox[\\\"a\\\", \\\ \"TI\\\"]}], \\\")\\\"}]}]\). \"\>", "MSG"], "\[NonBreakingSpace]", ButtonBox[ StyleBox["\[RightSkeleton]", "SR"], Active->True, BaseStyle->"Link", ButtonData->"paclet:ref/LerchPhi"]}]], "Print", "PrintUsage", CellChangeTimes->{3.44676046213206*^9}, CellTags->"Info3446746062-1890553"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Limit", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"]}], " ", RowBox[{"LerchPhi", "[", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", RowBox[{"1", "+", "n"}]}], "]"}]}], "+", RowBox[{"Log", "[", "2", "]"}]}], ",", RowBox[{"n", " ", "\[Rule]", " ", "Infinity"}]}], "]"}]], "Input", CellChangeTimes->{{3.446729922192877*^9, 3.446729945498026*^9}, { 3.4467594383616486`*^9, 3.4467594796480684`*^9}}], Cell[BoxData[ RowBox[{"Log", "[", "2", "]"}]], "Output", CellChangeTimes->{ 3.446730006548806*^9, {3.4467594684026947`*^9, 3.446759481080594*^9}, 3.446759558727968*^9, 3.446760463080242*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "the", " ", "following", " ", "is", " ", "a", " ", "slight", " ", "generalization", " ", "of", " ", "our", " ", "original", " ", RowBox[{"sequence", ".", " ", "It"}], " ", RowBox[{"isn", "'"}], "t", " ", "hard", " ", "to", " ", "show", " ", "that", " ", "if", " ", "the", " ", "first", " ", "equals", " ", "log", RowBox[{"(", "2", ")"}], " ", "then", " ", "this", " ", "equals", " ", "x"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"y", "-", "x"}], ")"}], " ", "log", RowBox[{"(", "2", ")"}]}]}], ";", " ", RowBox[{"Mathematica", " ", "sees", " ", "this", " ", "as", " ", RowBox[{"well", "!"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"RSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"n", "+", "1"}], "]"}], " ", "==", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", "n", "]"}], " ", "+", " ", RowBox[{"n", " ", RowBox[{"a", "[", RowBox[{"n", "-", "1"}], "]"}]}]}], ")"}], " ", "/", " ", RowBox[{"(", RowBox[{"n", "+", "1"}], ")"}]}]}], ",", " ", RowBox[{ RowBox[{"a", "[", "0", "]"}], " ", "\[Equal]", "x"}], ",", " ", RowBox[{ RowBox[{"a", "[", "1", "]"}], " ", "\[Equal]", "y"}]}], "}"}], ",", " ", RowBox[{"a", "[", "n", "]"}], ",", "n"}], "]"}]}]], "Input", CellChangeTimes->{{3.4467595694577503`*^9, 3.4467595818466053`*^9}, { 3.446759657894472*^9, 3.446759658961787*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"a", "[", "n", "]"}], "\[Rule]", RowBox[{"x", "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", "x", " ", RowBox[{"LerchPhi", "[", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", RowBox[{"1", "+", "n"}]}], "]"}]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", "y", " ", RowBox[{"LerchPhi", "[", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", RowBox[{"1", "+", "n"}]}], "]"}]}], "-", RowBox[{"x", " ", RowBox[{"Log", "[", "2", "]"}]}], "+", RowBox[{"y", " ", RowBox[{"Log", "[", "2", "]"}]}]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.4467595736742773`*^9, 3.446759582626829*^9}, 3.446759661144264*^9, 3.4467604641709924`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Limit", "[", RowBox[{ RowBox[{"x", "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", "x", " ", RowBox[{"LerchPhi", "[", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", RowBox[{"1", "+", "n"}]}], "]"}]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", "y", " ", RowBox[{"LerchPhi", "[", RowBox[{ RowBox[{"-", "1"}], ",", "1", ",", RowBox[{"1", "+", "n"}]}], "]"}]}], "-", RowBox[{"x", " ", RowBox[{"Log", "[", "2", "]"}]}], "+", RowBox[{"y", " ", RowBox[{"Log", "[", "2", "]"}]}]}], ",", " ", RowBox[{"n", " ", "\[Rule]", " ", "Infinity"}]}], "]"}]], "Input", CellChangeTimes->{{3.446759588317871*^9, 3.4467596529611945`*^9}}], Cell[BoxData[ RowBox[{"x", "-", RowBox[{"x", " ", RowBox[{"Log", "[", "2", "]"}]}], "+", RowBox[{"y", " ", RowBox[{"Log", "[", "2", "]"}]}]}]], "Output", CellChangeTimes->{3.4467596084464655`*^9, 3.446759662143224*^9, 3.446760466076144*^9}] }, Open ]] }, WindowSize->{673, 549}, WindowMargins->{{Automatic, 11}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{ "Info3446746051-8245676"->{ Cell[48398, 1209, 582, 13, 42, "Print", CellTags->"Info3446746051-8245676"]}, "Info3446746058-4040371"->{ Cell[56764, 1451, 845, 17, 62, "Print", CellTags->"Info3446746058-4040371"]}, "Info3446746062-1890553"->{ Cell[59854, 1537, 767, 15, 42, "Print", CellTags->"Info3446746062-1890553"]} } *) (*CellTagsIndex CellTagsIndex->{ {"Info3446746051-8245676", 65537, 1708}, {"Info3446746058-4040371", 65649, 1711}, {"Info3446746062-1890553", 65761, 1714} } *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 224, 5, 113, "Input"], Cell[772, 27, 1255, 26, 112, "Input"], Cell[2030, 55, 654, 16, 72, "Input"], Cell[2687, 73, 1115, 27, 112, "Input"], Cell[CellGroupData[{ Cell[3827, 104, 905, 22, 252, "Input"], Cell[CellGroupData[{ Cell[4757, 130, 1270, 28, 84, "Print"], Cell[6030, 160, 1439, 35, 62, "Print"] }, Open ]] }, Open ]], Cell[7496, 199, 94, 1, 31, "Input"], Cell[CellGroupData[{ Cell[7615, 204, 1743, 42, 132, "Input"], Cell[9361, 248, 596, 14, 33, "Output"], Cell[9960, 264, 687, 17, 33, "Output"] }, Open ]], Cell[10662, 284, 1781, 41, 174, "Input"], Cell[CellGroupData[{ Cell[12468, 329, 777, 17, 112, "Input"], Cell[13248, 348, 277, 7, 30, "Output"], Cell[13528, 357, 274, 7, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13839, 369, 486, 13, 52, "Input"], Cell[14328, 384, 6824, 117, 239, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21189, 506, 2220, 55, 152, "Input"], Cell[23412, 563, 2623, 49, 424, "Output"] }, Open ]], Cell[26050, 615, 92, 1, 31, "Input"], Cell[CellGroupData[{ Cell[26167, 620, 2136, 53, 172, "Input"], Cell[28306, 675, 707, 16, 30, "Output"], Cell[29016, 693, 614, 13, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[29667, 711, 3511, 81, 252, "Input"], Cell[33181, 794, 2790, 52, 376, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[36008, 851, 2664, 70, 152, "Input"], Cell[38675, 923, 2022, 56, 83, "Output"], Cell[40700, 981, 1174, 33, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[41911, 1019, 2378, 61, 192, "Input"], Cell[44292, 1082, 1006, 30, 51, "Output"], Cell[45301, 1114, 852, 27, 55, "Output"], Cell[46156, 1143, 1080, 34, 73, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[47273, 1182, 1122, 25, 152, "Input"], Cell[48398, 1209, 582, 13, 42, "Print", CellTags->"Info3446746051-8245676"] }, Open ]], Cell[CellGroupData[{ Cell[49017, 1227, 3226, 86, 172, "Input"], Cell[52246, 1315, 543, 14, 30, "Output"], Cell[52792, 1331, 576, 16, 33, "Output"], Cell[53371, 1349, 601, 17, 47, "Output"], Cell[53975, 1368, 532, 14, 33, "Output"], Cell[54510, 1384, 2067, 58, 206, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[56614, 1447, 147, 2, 31, "Input"], Cell[56764, 1451, 845, 17, 62, "Print", CellTags->"Info3446746058-4040371"] }, Open ]], Cell[57624, 1471, 782, 14, 92, "Input"], Cell[CellGroupData[{ Cell[58431, 1489, 804, 22, 52, "Input"], Cell[59238, 1513, 520, 16, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[59795, 1534, 56, 1, 31, "Input"], Cell[59854, 1537, 767, 15, 42, "Print", CellTags->"Info3446746062-1890553"] }, Open ]], Cell[CellGroupData[{ Cell[60658, 1557, 557, 16, 33, "Input"], Cell[61218, 1575, 198, 4, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61453, 1584, 1672, 42, 112, "Input"], Cell[63128, 1628, 918, 27, 55, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[64083, 1660, 847, 25, 55, "Input"], Cell[64933, 1687, 259, 7, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)