(* 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[ 134802, 3919] NotebookOptionsPosition[ 128428, 3724] NotebookOutlinePosition[ 129025, 3748] CellTagsIndexPosition[ 128938, 3743] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "Set", " ", "the", " ", "location", " ", "and", " ", "load", " ", "the", " ", "package", " ", "DDESpecialSolutions"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"SetDirectory", "[", "\"\\"", "]"}], ";", RowBox[{"Get", "[", "\"\\"", "]"}]}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"Package DDESpecialSolutionsV2.m was successfully \ loaded.\"\>"], "Print", CellChangeTimes->{ 3.4419422547749147`*^9, 3.4419424052749147`*^9, {3.4419427680717897`*^9, 3.4419427841499147`*^9}}], Cell[BoxData["\<\"Version 2 was last updated: Monday, January 26, 2009 by \ Baldwin.\"\>"], "Print", CellChangeTimes->{ 3.4419422547749147`*^9, 3.4419424052749147`*^9, {3.4419427680717897`*^9, 3.4419427841499147`*^9}}], Cell[BoxData["\<\"Version 1 was first released on July 23, 2004.\"\>"], \ "Print", CellChangeTimes->{ 3.4419422547749147`*^9, 3.4419424052749147`*^9, {3.4419427680717897`*^9, 3.4419427841499147`*^9}}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"??", " ", "DDESpecialSolutions"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ StyleBox["\<\"DDESpecialSolutions[eqns, funcs, dvars, cvars, params, opts] \ solves a system of nonlinear differential-difference equations for funcs, \ with continuous independent variables cvars ({x,y,z,...}), discrete \ independent variables dvars ({n,m,...}) and non-zero parameters params in \ terms of hyperbolic tangent functions.\"\>", "MSG"]], "Print", "PrintUsage", CellChangeTimes->{3.4419428098061647`*^9}, CellTags->"Info3441917609-4281003"], Cell[BoxData[ InterpretationBox[GridBox[{ { RowBox[{ RowBox[{"Attributes", "[", "DDESpecialSolutions", "]"}], "=", RowBox[{"{", "ReadProtected", "}"}]}]}, {" "}, {GridBox[{ { RowBox[{ RowBox[{"Options", "[", "DDESpecialSolutions", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"Verbose", "\[Rule]", "False"}], ",", RowBox[{"\<\"InputForm\"\>", "\[Rule]", "False"}], ",", RowBox[{"NumericTest", "\[Rule]", "True"}], ",", RowBox[{"SymbolicTest", "\[Rule]", "False"}], ",", RowBox[{"DegreeOfThePolynomial", "\[Rule]", RowBox[{"{", "}"}]}], ",", RowBox[{"MinDegreeOfThePolynomial", "\[Rule]", "1"}], ",", RowBox[{"MaxDegreeOfThePolynomial", "\[Rule]", "3"}], ",", RowBox[{"SolveAlgebraicSystem", "\[Rule]", "True"}]}], "}"}]}]} }, BaselinePosition->{Baseline, {1, 1}}, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxItemSize->{"Columns" -> {{ Scaled[0.999]}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, "RowsIndexed" -> {}}]} }, BaselinePosition->{Baseline, {1, 1}}, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}], Definition[Calculus`DDESpecialSolutions`DDESpecialSolutions], Editable->False]], "Print", CellChangeTimes->{3.4419428098217897`*^9}, CellTags->"Info3441917609-4281003"] }, Open ]] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"In", " ", "difficult", " ", "cases"}], ",", " ", RowBox[{ "there", " ", "may", " ", "be", " ", "freedom", " ", "in", " ", "the", " ", "choice", " ", "of", " ", "the", " ", "degree", " ", "of", " ", "the", " ", "polynomial", " ", RowBox[{"solution", ".", " ", "Thus"}]}], ",", " ", RowBox[{ "the", " ", "software", " ", "investigates", " ", "all", " ", "the", " ", "possible", " ", "solutions", " ", "between", " ", "MinDegreeOfThePolynomial", " ", "and", " ", RowBox[{"MaxDegreeOfThePolynomial", "."}]}]}], " ", "*)"}]], "Input"], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Single", " ", "Toda", " ", "Lattice"}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", "2"}], "}"}]}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}], "-", RowBox[{"2", "*", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], "}"}], ",", RowBox[{"{", "u", "}"}], ",", " ", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", " ", RowBox[{"{", "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424197749147`*^9, 3.4419428102436647`*^9}], Cell["\<\ The following simplification rules are being used to test the solutions: \ {Sqrt[a^2]->a}\ \>", "Message", CellChangeTimes->{3.4419428102436647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "-", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "+", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424200874147`*^9, 3.4419428102592897`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Toda", " ", "as", " ", "System", " ", "of", " ", "First"}], "-", "order"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Use", " ", "the", " ", "Verbose"}], "\[Rule]", "True"}], ",", " ", RowBox[{ "option", " ", "to", " ", "see", " ", "a", " ", "trace", " ", "of", " ", "the", " ", "steps", " ", "of", " ", "the", " ", RowBox[{"algorithm", "."}]}]}], " ", "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], "\[Equal]", "0"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "-", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}], "-", RowBox[{"2", "*", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], "\[Equal]", "0"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"u", ",", " ", "v"}], "}"}], ",", " ", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", " ", RowBox[{"{", "}"}], ",", RowBox[{"Verbose", " ", "\[Rule]", " ", "True"}]}], "]"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData["\<\"The given system of differential-difference equations: \ \"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428102905397`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ InterpretationBox[ RowBox[{"u", "\[InvisibleSpace]", InterpretationBox[ SubscriptBox["\[Null]", InterpretationBox[ RowBox[{ "n", "\[InvisibleSpace]", "\<\";\"\>", "\[InvisibleSpace]", "t"}], SequenceForm[$CellContext`n, ";", $CellContext`t], Editable->False]], Subscript[ SequenceForm[$CellContext`n, ";", $CellContext`t]], Editable->False]}], SequenceForm[$CellContext`u, Subscript[ SequenceForm[$CellContext`n, ";", $CellContext`t]]], Editable->False], "-", SubscriptBox["v", "n"]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ InterpretationBox[ RowBox[{"v", "\[InvisibleSpace]", InterpretationBox[ SubscriptBox["\[Null]", InterpretationBox[ RowBox[{ "n", "\[InvisibleSpace]", "\<\";\"\>", "\[InvisibleSpace]", "t"}], SequenceForm[$CellContext`n, ";", $CellContext`t], Editable->False]], Subscript[ SequenceForm[$CellContext`n, ";", $CellContext`t]], Editable->False]}], SequenceForm[$CellContext`v, Subscript[ SequenceForm[$CellContext`n, ";", $CellContext`t]]], Editable->False], "-", RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["u", RowBox[{ RowBox[{"-", "1"}], "+", "n"}]], "-", RowBox[{"2", " ", SubscriptBox["u", "n"]}], "+", SubscriptBox["u", RowBox[{"1", "+", "n"}]]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", SubscriptBox["v", "n"]}], ")"}]}]}], "\[Equal]", "0"}]}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428102905397`*^9}], Cell[BoxData["\<\"Transform the system of differential-difference equations \ into nonlinear DDEs in T = Tanh.\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428102905397`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", "n", "]"}], "\[Prime]", MultilineFunction->None], "[", "T", "]"}]}], "-", RowBox[{ SuperscriptBox["T", "2"], " ", RowBox[{"c", "[", "2", "]"}], " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", "n", "]"}], "\[Prime]", MultilineFunction->None], "[", "T", "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "[", "T", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "\[Prime]", MultilineFunction->None], "[", "T", "]"}]}], "-", RowBox[{ SuperscriptBox["T", "2"], " ", RowBox[{"c", "[", "2", "]"}], " ", RowBox[{ SuperscriptBox[ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "\[Prime]", MultilineFunction->None], "[", "T", "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", RowBox[{ RowBox[{"-", "1"}], "+", "n"}], "]"}], "[", "T", "]"}], "+", RowBox[{"2", " ", RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", "n", "]"}], "[", "T", "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", RowBox[{"1", "+", "n"}], "]"}], "[", "T", "]"}], "-", RowBox[{ RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", RowBox[{ RowBox[{"-", "1"}], "+", "n"}], "]"}], "[", "T", "]"}], " ", RowBox[{ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "[", "T", "]"}]}], "+", RowBox[{"2", " ", RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", "n", "]"}], "[", "T", "]"}], " ", RowBox[{ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "[", "T", "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", RowBox[{"1", "+", "n"}], "]"}], "[", "T", "]"}], " ", RowBox[{ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "[", "T", "]"}]}]}]}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428102905397`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Time Used:\"\>", "\[InvisibleSpace]", "0.`"}], SequenceForm["Time Used:", 0.], Editable->False]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428102905397`*^9}], Cell[BoxData["\<\"Determine the maximal degree of the polynomial solutions in \ T.\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428103061647`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \ \[Rule] \\(\\(-1\\)\\)\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 0\\)\\)}\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" are being removed \ because they are (i) negative, (ii) contain freedom, (iii) fail to balance \ highest exponent terms from two different terms in the original system. If \ \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the transformation u -> 1/v \ may result in a system that DDESpecialSolutions can solve.\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424281030397`*^9, 3.4419428103530397`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \ \[Rule] \\(\\(-1\\)\\)\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 0\\)\\)}\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" are being removed \ because they are (i) negative, (ii) contain freedom, (iii) fail to balance \ highest exponent terms from two different terms in the original system. If \ \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the transformation u -> 1/v \ may result in a system that DDESpecialSolutions can solve.\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424281030397`*^9, 3.4419428103686647`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"m", "[", "1", "]"}], "\[Rule]", "1"}], ",", RowBox[{ RowBox[{"m", "[", "2", "]"}], "\[Rule]", "2"}]}], "}"}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428103842897`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Time Used:\"\>", "\[InvisibleSpace]", "0.046999999999999986`"}], SequenceForm["Time Used:", 0.046999999999999986`], Editable->False]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428103999147`*^9}], Cell[BoxData["\<\"Derive the nonlinear algebraic system for the coefficients.\ \"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428104155397`*^9}], Cell[BoxData["\<\"Seek polynomial solutions in T of the form:\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428104311647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"u", "[", "1", "]"}], "[", "n", "]"}], "[", "T", "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "+", RowBox[{"T", " ", RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}]}]}]}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"u", "[", "2", "]"}], "[", "n", "]"}], "[", "T", "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], "+", RowBox[{"T", " ", RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]}], "+", RowBox[{ SuperscriptBox["T", "2"], " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}]}]}]}]}], "}"}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428104467897`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"where T = Tanh[\"\>", "\[InvisibleSpace]", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "\[InvisibleSpace]", "\<\"].\"\>"}], SequenceForm[ "where T = Tanh[", $CellContext`phase + $CellContext`n $CellContext`c[ 1] + $CellContext`t $CellContext`c[2], "]."], Editable->False]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428104780397`*^9}], Cell[BoxData["\<\"The nonlinear algebraic system is: \"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428104936647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], "+", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], "-", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}]}], "-", RowBox[{"c", "[", "2", "]"}]}], ")"}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], "-", RowBox[{"c", "[", "2", "]"}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}], "+", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "+", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"c", "[", "2", "]"}], "-", RowBox[{"2", " ", RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "+", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], ")"}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}], "+", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "+", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "-", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "+", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], ",", RowBox[{"c", "[", "2", "]"}]}], "}"}], ",", RowBox[{"{", "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], ",", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], ",", RowBox[{"c", "[", "2", "]"}]}], "}"}]}], "}"}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428105092897`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Time Used:\"\>", "\[InvisibleSpace]", "0.016000000000000014`"}], SequenceForm["Time Used:", 0.016000000000000014`], Editable->False]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428105249147`*^9}], Cell[BoxData["\<\"Solve the nonlinear algebraic system.\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428105561647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], SqrtBox[ RowBox[{"1", "-", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]]]}]}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], "\[Rule]", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], "\[Rule]", RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], ",", RowBox[{ RowBox[{"c", "[", "2", "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], SqrtBox[ RowBox[{"1", "-", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]]]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], SqrtBox[ RowBox[{"1", "-", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]]]}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], "\[Rule]", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], ",", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], "\[Rule]", RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], ",", RowBox[{ RowBox[{"c", "[", "2", "]"}], "\[Rule]", FractionBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], SqrtBox[ RowBox[{"1", "-", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]]]}]}], "}"}]}], "}"}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428105717897`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Time Used:\"\>", "\[InvisibleSpace]", "0.062`"}], SequenceForm["Time Used:", 0.062], Editable->False]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428105874147`*^9}], Cell[BoxData["\<\"Build and test the travelling wave solutions.\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428106186647`*^9}], Cell[BoxData["\<\"The possible non-trivial solutions (before any testing): \"\ \>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428106342897`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "-", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "+", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428106655397`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428106967897`*^9}] }, Open ]], Cell["\<\ The following simplification rules are being used to test the solutions: \ {Sqrt[a^2]->a}\ \>", "Message", CellChangeTimes->{3.4419428107124147`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "-", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "+", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}], "}"}]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428107436647`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Time Used:\"\>", "\[InvisibleSpace]", "0.030999999999999972`"}], SequenceForm["Time Used:", 0.030999999999999972`], Editable->False]], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428107592897`*^9}], Cell[BoxData["\<\"\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428107905397`*^9}], Cell[BoxData["\<\"The (numerically and/or symbolically) tested \ solutions:\"\>"], "Print", CellChangeTimes->{3.4419424280249147`*^9, 3.4419428108217897`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "-", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], "+", RowBox[{ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424289624147`*^9, 3.4419428108374147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"Relativistic", " ", "Toda", " ", "Equation"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"To", " ", "use", " ", "the", " ", "symbolic", " ", "test"}], ",", " ", RowBox[{ RowBox[{"use", " ", "the", " ", "SymbolicTest"}], "\[Rule]", RowBox[{"True", " ", RowBox[{"option", "."}]}]}]}], " ", "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], " ", ",", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{"alpha", "*", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", " ", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "+", RowBox[{"alpha", "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], "-", RowBox[{"alpha", "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}]}], ")"}]}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"u", ",", "v"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", "alpha", "}"}], ",", RowBox[{"SymbolicTest", "\[Rule]", "True"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \ \[Rule] 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 2\\)\\)}\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" are being removed \ because they are (i) negative, (ii) contain freedom, (iii) fail to balance \ highest exponent terms from two different terms in the original system. If \ \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the transformation u -> 1/v \ may result in a system that DDESpecialSolutions can solve.\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424328686647`*^9, 3.4419428109155397`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \ \[Rule] 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 2\\)\\)}\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" are being removed \ because they are (i) negative, (ii) contain freedom, (iii) fail to balance \ highest exponent terms from two different terms in the original system. If \ \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the transformation u -> 1/v \ may result in a system that DDESpecialSolutions can solve.\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424328686647`*^9, 3.4419428109311647`*^9}], Cell[CellGroupData[{ Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424329467897`*^9, 3.4419428109936647`*^9}], Cell[BoxData["\<\"Symbolically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424329467897`*^9, 3.4419428110092897`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"1", "+", RowBox[{"alpha", " ", RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"alpha", " ", RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}]}], "alpha"]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "-", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}], "alpha"]}]}], "}"}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424330874147`*^9, 3.4419428111342897`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "To", " ", "force", " ", "the", " ", "function", " ", "to", " ", "use", " ", "your", " ", "choice", " ", "of", " ", "ansatz"}], ",", " ", RowBox[{"use", " ", "the", " ", "DegreeOfThePolynomial", " ", RowBox[{"option", "."}]}]}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"n", ",", "t"}], "]"}], "==", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{"alpha", " ", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"v", "[", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "+", "n"}], ",", "t"}], "]"}]}], "+", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{ SuperscriptBox["v", TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"n", ",", "t"}], "]"}], "==", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"1", "+", "n"}], ",", "t"}], "]"}], "-", RowBox[{"alpha", " ", RowBox[{"v", "[", RowBox[{ RowBox[{ RowBox[{"-", "1"}], "+", "n"}], ",", "t"}], "]"}]}], "+", RowBox[{"alpha", " ", RowBox[{"v", "[", RowBox[{ RowBox[{"1", "+", "n"}], ",", "t"}], "]"}]}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", "v"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", "alpha", "}"}], ",", " ", RowBox[{"DegreeOfThePolynomial", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"m", "[", "1", "]"}], "\[Rule]", "1"}], ",", RowBox[{ RowBox[{"m", "[", "2", "]"}], "\[Rule]", "1"}]}], "}"}]}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{ "DDESpecialSolutionsmSolver", "::", "\<\"DegreeOfThePolynomial\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Algorithm continuing with \ user-supplied exponent of the polynomial \\\\\\\"\\\", \\\"MT\\\"]\\)\ \[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \[Rule] \ 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 1\\)\\)}\\)}\\), \\\"MT\\\"]\\)\ \[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\".\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424363842897`*^9, 3.4419428112436647`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424364624147`*^9, 3.4419428113061647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"1", "+", RowBox[{"alpha", " ", RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"alpha", " ", RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}]}], "alpha"]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "-", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}], "alpha"]}]}], "}"}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424365874147`*^9, 3.4419428113374147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Volterra", " ", "Lattice"}], " ", "*)"}]], "Input"], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "To", " ", "get", " ", "the", " ", "output", " ", "in", " ", "a", " ", "copyable", " ", "form"}], ",", " ", RowBox[{ RowBox[{"use", " ", "the", " ", "InputForm"}], "\[Rule]", RowBox[{"True", " ", RowBox[{"option", "."}]}]}]}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", " ", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}]}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"u", ",", " ", "v"}], "}"}], ",", " ", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", " ", RowBox[{"{", "}"}], ",", " ", RowBox[{"InputForm", "\[Rule]", "True"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424391342897`*^9, 3.4419428114467897`*^9}], Cell["\<\ {{{u[n, t] -> -(c[2]*(Coth[c[1]] - Tanh[phase + n*c[1] + t*c[2]])), v[n, t] -> -(c[2]*(Coth[c[1]] + Tanh[phase + n*c[1] + t*c[2]]))}}}\ \>", "Output", CellChangeTimes->{3.4419424392280397`*^9, 3.4419428114467897`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "If", " ", "you", " ", "would", " ", "like", " ", "to", " ", "solve", " ", "the", " ", "algebraic", " ", "system", " ", "using", " ", "other", " ", "software"}], ",", " ", RowBox[{ RowBox[{"use", " ", "the", " ", "SolveAlgebraicSystem"}], "\[Rule]", RowBox[{"False", " ", RowBox[{"option", "."}]}]}]}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", " ", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}]}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"u", ",", " ", "v"}], "}"}], ",", " ", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", " ", RowBox[{"{", "}"}], ",", " ", RowBox[{"SolveAlgebraicSystem", "\[Rule]", "False"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}], "-", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}], "-", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], "+", RowBox[{"c", "[", "2", "]"}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], "+", RowBox[{"c", "[", "2", "]"}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}], "-", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"c", "[", "2", "]"}]}], "-", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], "-", RowBox[{"c", "[", "2", "]"}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], "-", RowBox[{"c", "[", "2", "]"}]}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"2", ",", "0"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], ",", RowBox[{"c", "[", "2", "]"}]}], "}"}], ",", RowBox[{"{", "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], ",", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], ",", RowBox[{"c", "[", "2", "]"}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424434780397`*^9, 3.4419428115249147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "A", " ", "new", " ", "integrable", " ", "system", " ", "related", " ", "to", " ", "the", " ", "relativistic", " ", "Toda", " ", "lattice"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Ref", ":", "Suris"}], ",", RowBox[{ RowBox[{"J", ".", "Phys", ".", "A", ".", "30"}], " ", RowBox[{"(", "1997", ")"}], " ", "1745"}]}], "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}], "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}]}]}], ",", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", " ", "v"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\({\\(\\(\\(\\(m[1]\\)\\)\ \[Rule] 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 0\\)\\)}\\)\\), \ \\(\\({\\(\\(\\(\\(m[1]\\)\\) \[Rule] 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] \ 2\\)\\)}\\)\\)}\\), \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" \ are being removed because they are (i) negative, (ii) contain freedom, (iii) \ fail to balance highest exponent terms from two different terms in the \ original system. If \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the \ transformation u -> 1/v may result in a system that DDESpecialSolutions can \ solve.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424499624147`*^9, 3.4419428116186647`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\({\\(\\(\\(\\(m[1]\\)\\)\ \[Rule] 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] \\(\\(-1\\)\\)\\)\\)}\\)\\), \ \\(\\({\\(\\(\\(\\(m[1]\\)\\) \[Rule] 1\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] \ 3\\)\\)}\\)\\)}\\), \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" \ are being removed because they are (i) negative, (ii) contain freedom, (iii) \ fail to balance highest exponent terms from two different terms in the \ original system. If \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the \ transformation u -> 1/v may result in a system that DDESpecialSolutions can \ solve.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419424499624147`*^9, 3.4419428116499147`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424500717897`*^9, 3.4419428117436647`*^9}], Cell["\<\ These solutions did not simplify to be less than 10^(-7.97729) when random \ numbers in [-1,1] were substituted for the unknowns. Please test these \ solutions by hand or interactively with Mathematica.\ \>", "Message", CellChangeTimes->{3.4419428117592897`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", RowBox[{"c", "[", "2", "]"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", RowBox[{"c", "[", "2", "]"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}]], "Print", CellChangeTimes->{3.4419424500717897`*^9, 3.4419428118061647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424503842897`*^9, 3.4419428118217897`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "On", " ", "some", " ", "integrable", " ", "systems", " ", "related", " ", "to", " ", "the", " ", "Toda", " ", "lattice"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Ref", ":", "Suris"}], ",", RowBox[{ RowBox[{"J", ".", "Phys", ".", "A", ".", "30"}], " ", RowBox[{"(", "1997", ")"}], " ", "2235"}]}], "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "-", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", "v"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424913530397`*^9, 3.4419428119311647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", RowBox[{"c", "[", "2", "]"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "-", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", RowBox[{"c", "[", "2", "]"}]}], " ", RowBox[{"(", RowBox[{ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}]}], ")"}]}]}]}], "}"}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419424914467897`*^9, 3.4419428119311647`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Schrodinger", " ", "System"}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "-", RowBox[{"alpha", "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"2", "*", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "+", RowBox[{"alpha", "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"2", "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "+", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", " ", "v"}], "}"}], ",", " ", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", " ", RowBox[{"{", "alpha", "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"AnalyzeAndSolve", "::", "\<\"soln\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The solution \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\(\\(\\(a[\\(\\(1, \ 0\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(4\\\\ alpha\\)\\/\\(\\((\\(\\(-1\\)\\) \ + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 0\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(a[\\(\\(1, 2\\)\\)]\\)\\) \[Rule] \\(8\\\\ alpha\ \\)\\/\\(\\((\\(\\(-1\\)\\) + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \ \\(\\(a[\\(\\(2, 2\\)\\)]\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(1, \ 1\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(2\\\\ \\(\\(c[2]\\)\\)\\)\\/\\(\\(\\((\ \\(\\(-1\\)\\) + \\@2)\\)\\)\\\\ \\(\\((1 + \\@2)\\)\\)\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 1\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(Coth[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] 0\\)\\), \ \\(\\(\\(\\(Tanh[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] \\(\\(-\\@2\\)\\)\\)\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is being removed \ because it is inconsistent.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", \ "MSG", CellChangeTimes->{3.4419424994467897`*^9, 3.4419428131342897`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"AnalyzeAndSolve", "::", "\<\"soln\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The solution \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\(\\(\\(a[\\(\\(1, \ 0\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(4\\\\ alpha\\)\\/\\(\\((\\(\\(-1\\)\\) \ + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 0\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(a[\\(\\(1, 2\\)\\)]\\)\\) \[Rule] \\(8\\\\ alpha\ \\)\\/\\(\\((\\(\\(-1\\)\\) + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \ \\(\\(a[\\(\\(2, 2\\)\\)]\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(1, \ 1\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(2\\\\ \\(\\(c[2]\\)\\)\\)\\/\\(\\(\\((\ \\(\\(-1\\)\\) + \\@2)\\)\\)\\\\ \\(\\((1 + \\@2)\\)\\)\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 1\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(Coth[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] 0\\)\\), \ \\(\\(\\(\\(Tanh[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] \\@2\\)\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is being removed \ because it is inconsistent.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", \ "MSG", CellChangeTimes->{3.4419424994467897`*^9, 3.4419428131655397`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419424994780397`*^9, 3.4419428131811647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"alpha", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "alpha", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}], RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "alpha", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"alpha", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", "alpha", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}], RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", "alpha", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419425009311647`*^9, 3.4419428132124147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Ablowitz", "-", RowBox[{ "Ladik", " ", "Discretization", " ", "of", " ", "NLS", " ", "Equation"}]}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"aa", "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"2", "*", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}], "+", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "==", RowBox[{ RowBox[{ RowBox[{"-", "aa"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"2", "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}], "-", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "+", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", " ", "v"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", "aa", "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"AnalyzeAndSolve", "::", "\<\"soln\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The solution \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\(\\(\\(a[\\(\\(1, \ 0\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(4\\\\ aa\\)\\/\\(\\((\\(\\(-1\\)\\) + \ \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 0\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(a[\\(\\(1, 2\\)\\)]\\)\\) \[Rule] \\(8\\\\ aa\\)\ \\/\\(\\((\\(\\(-1\\)\\) + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\ \\(\\(2, 2\\)\\)]\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(1, 1\\)\\)]\\)\\) \ \[Rule] \\(\\(-\\(\\(\\(2\\\\ \\(\\(c[2]\\)\\)\\)\\/\\(\\(\\((\\(\\(-1\\)\\) \ + \\@2)\\)\\)\\\\ \\(\\((1 + \\@2)\\)\\)\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 1\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(Coth[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] 0\\)\\), \ \\(\\(\\(\\(Tanh[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] \\(\\(-\\@2\\)\\)\\)\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is being removed \ because it is inconsistent.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", \ "MSG", CellChangeTimes->{3.4419425244780397`*^9, 3.4419428144467897`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"AnalyzeAndSolve", "::", "\<\"soln\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The solution \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\(\\(\\(a[\\(\\(1, \ 0\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(4\\\\ aa\\)\\/\\(\\((\\(\\(-1\\)\\) + \ \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 0\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(a[\\(\\(1, 2\\)\\)]\\)\\) \[Rule] \\(8\\\\ aa\\)\ \\/\\(\\((\\(\\(-1\\)\\) + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\ \\(\\(2, 2\\)\\)]\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(1, 1\\)\\)]\\)\\) \ \[Rule] \\(\\(-\\(\\(\\(2\\\\ \\(\\(c[2]\\)\\)\\)\\/\\(\\(\\((\\(\\(-1\\)\\) \ + \\@2)\\)\\)\\\\ \\(\\((1 + \\@2)\\)\\)\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 1\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(Coth[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] 0\\)\\), \ \\(\\(\\(\\(Tanh[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] \\@2\\)\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is being removed \ because it is inconsistent.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", \ "MSG", CellChangeTimes->{3.4419425244780397`*^9, 3.4419428144624147`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419425245092897`*^9, 3.4419428145092897`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"aa", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "aa", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}], RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "aa", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"aa", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", "aa", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}], RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", "aa", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419425259467897`*^9, 3.4419428145405397`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"NLS", " ", "Equation", " ", RowBox[{"(", RowBox[{"Ablowitz", "-", RowBox[{"Ladik", " ", "Discretization"}]}], ")"}]}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "==", RowBox[{ RowBox[{"aa", "*", RowBox[{"(", RowBox[{ RowBox[{"bb", "*", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], "-", RowBox[{"2", "*", "bb", "*", RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"bb", "*", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"bb", "*", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], "+", RowBox[{"bb", "*", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}]}], ")"}]}]}]}], ",", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "==", RowBox[{ RowBox[{ RowBox[{"-", "aa"}], "*", RowBox[{"(", RowBox[{ RowBox[{"bb", "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], "-", RowBox[{"2", "*", "bb", "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"bb", "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"bb", "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], "+", RowBox[{"bb", "*", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}]}], ")"}]}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", "v"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", RowBox[{"aa", ",", "bb"}], "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"AnalyzeAndSolve", "::", "\<\"soln\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The solution \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\(\\(\\(a[\\(\\(1, \ 0\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(4\\\\ aa\\)\\/\\(\\((\\(\\(-1\\)\\) + \ \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 0\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(a[\\(\\(1, 2\\)\\)]\\)\\) \[Rule] \\(8\\\\ aa\\)\ \\/\\(\\((\\(\\(-1\\)\\) + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\ \\(\\(2, 2\\)\\)]\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(1, 1\\)\\)]\\)\\) \ \[Rule] \\(\\(-\\(\\(\\(2\\\\ \\(\\(c[2]\\)\\)\\)\\/\\(\\(\\((\\(\\(-1\\)\\) \ + \\@2)\\)\\)\\\\ \\(\\((1 + \\@2)\\)\\)\\\\ bb\\\\ \\(\\(a[\\(\\(2, 2\\)\\)]\ \\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 1\\)\\)]\\)\\) \[Rule] \ 0\\)\\), \\(\\(\\(\\(Coth[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] 0\\)\\), \ \\(\\(\\(\\(Tanh[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] \\(\\(-\\@2\\)\\)\\)\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is being removed \ because it is inconsistent.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", \ "MSG", CellChangeTimes->{3.4419425326342897`*^9, 3.4419428159311647`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"AnalyzeAndSolve", "::", "\<\"soln\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The solution \\\\\\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\(\\(\\(\\(a[\\(\\(1, \ 0\\)\\)]\\)\\) \[Rule] \\(\\(-\\(\\(\\(4\\\\ aa\\)\\/\\(\\((\\(\\(-1\\)\\) + \ \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\\(\\(2, \ 2\\)\\)]\\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 0\\)\\)]\\)\\) \ \[Rule] 0\\)\\), \\(\\(\\(\\(a[\\(\\(1, 2\\)\\)]\\)\\) \[Rule] \\(8\\\\ aa\\)\ \\/\\(\\((\\(\\(-1\\)\\) + \\@2)\\)\\^2\\\\ \\((1 + \\@2)\\)\\^2\\\\ \\(\\(a[\ \\(\\(2, 2\\)\\)]\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(1, 1\\)\\)]\\)\\) \ \[Rule] \\(\\(-\\(\\(\\(2\\\\ \\(\\(c[2]\\)\\)\\)\\/\\(\\(\\((\\(\\(-1\\)\\) \ + \\@2)\\)\\)\\\\ \\(\\((1 + \\@2)\\)\\)\\\\ bb\\\\ \\(\\(a[\\(\\(2, 2\\)\\)]\ \\)\\)\\)\\)\\)\\)\\)\\)\\), \\(\\(\\(\\(a[\\(\\(2, 1\\)\\)]\\)\\) \[Rule] \ 0\\)\\), \\(\\(\\(\\(Coth[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] 0\\)\\), \ \\(\\(\\(\\(Tanh[\\(\\(c[1]\\)\\)]\\)\\) \[Rule] \\@2\\)\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" is being removed \ because it is inconsistent.\\\\\\\"\\\", \\\"MT\\\"]\\)\"\>"}]], "Message", \ "MSG", CellChangeTimes->{3.4419425326342897`*^9, 3.4419428159624147`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419425326655397`*^9, 3.4419428159936647`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"aa", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "aa", " ", "bb", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}], RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "aa", " ", "bb", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"aa", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", "aa", " ", "bb", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}], RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}]]}]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"a", "[", RowBox[{"2", ",", "1"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", "aa", " ", "bb", " ", "t", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], "]"}]}], ")"}]}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419425343217897`*^9, 3.4419428160249147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Modified", " ", "KdV", " ", "Equation", " ", RowBox[{"(", RowBox[{"with", " ", "parameter"}], ")"}]}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], "\[Equal]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"aa", "^", "2"}], "-", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "^", "2"}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", "u", ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", "aa", "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419425394311647`*^9, 3.4419428164155397`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", "aa"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", SuperscriptBox["aa", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"aa", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"2", " ", SuperscriptBox["aa", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419425764936647`*^9, 3.4419428164311647`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"1", "-", RowBox[{"D", " ", "Volterra", " ", "system"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "On", " ", "the", " ", "Dynamics", " ", "of", " ", "Rational", " ", "Solutions", " ", "for", " ", "1", "D", " ", "Volterra", " ", "System", " ", "by", " ", RowBox[{"A", ".", "S", ".", "Carstea"}], " ", "in", " ", "SIDE", " ", RowBox[{"III", "."}]}], " ", "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{"a", "+", RowBox[{"b", "*", RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}]}], "-", RowBox[{"c", "*", RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], "^", "2"}]}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"r", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], " ", "-", " ", RowBox[{"r", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{"a", "+", RowBox[{"b", "*", RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}]}], " ", "-", " ", RowBox[{"c", "*", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], "^", "2"}]}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"q", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}], "-", RowBox[{"q", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}]}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"q", ",", "r"}], "}"}], ",", RowBox[{"{", "n", "}"}], ",", RowBox[{"{", "t", "}"}], ",", RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"The potential solutions {{m[1] -> 1, m[2] -> 1}} are being \ removed because they are underdetermined. \"\>"], "Message", CellChangeTimes->{3.4419428165405397`*^9}], Cell[BoxData["\<\"The potential solutions {{m[1] -> 2, m[2] -> 2}} are being \ removed because they are underdetermined. \"\>"], "Message", CellChangeTimes->{3.4419428165561647`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419425825561647`*^9, 3.4419428175092897`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{ RowBox[{"-", "b"}], " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}]}], "-", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "+", RowBox[{"b", " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "+", RowBox[{ RowBox[{"c", "[", "2", "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "3"], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}], "2"]}]}], RowBox[{"2", " ", "c", " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], ")"}]}]]}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{ RowBox[{"-", "b"}], "+", RowBox[{"2", " ", "c", " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", SuperscriptBox[ RowBox[{"Coth", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], "-", RowBox[{"2", " ", "c", " ", RowBox[{"a", "[", RowBox[{"2", ",", "2"}], "]"}], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}], "2"]}]}], RowBox[{"2", " ", "c"}]]}]}], ",", RowBox[{"a", "\[Rule]", RowBox[{"-", FractionBox[ SuperscriptBox["b", "2"], RowBox[{"4", " ", "c"}]]}]}]}], "}"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "-", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "-", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "-", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "+", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "-", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "+", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"q", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "+", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}], ",", RowBox[{ RowBox[{"r", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{"b", "+", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], "+", RowBox[{"4", " ", "a", " ", "c"}]}]], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "+", RowBox[{"4", " ", "a", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.4419428177749147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"Self", "-", RowBox[{"dual", " ", "nonlinear", " ", "network", " ", "equation"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Time", "-", RowBox[{ "discretization", " ", "of", " ", "soliton", " ", "equations", " ", "by", " ", "Ryogo", " ", "Hirota", " ", "and", " ", "masataka", " ", "Iwao", " ", "in", " ", "SIDE", " ", RowBox[{"III", "."}]}]}], " ", "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "^", "2"}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}], " ", "-", " ", RowBox[{"v", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", "v", ",", "n", ",", "t", ",", RowBox[{"{", "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419428181499147`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], " ", RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "-", RowBox[{"2", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], "]"}]}]}], "}"}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419428181499147`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"Hybrid", " ", "Equation"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Time", "-", RowBox[{ "discretization", " ", "of", " ", "soliton", " ", "equations", " ", "by", " ", "Ryogo", " ", "Hirota", " ", "and", " ", "masataka", " ", "Iwao", " ", "in", " ", "SIDE", " ", RowBox[{"III", "."}]}]}], " ", "*)"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{"b", "*", RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}]}], "+", RowBox[{"c", "*", RowBox[{ RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}], "^", "2"}]}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"w", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "t"}], "]"}], " ", "-", " ", RowBox[{"w", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "t"}], "]"}]}], ")"}]}]}], ",", "w", ",", "n", ",", "t", ",", RowBox[{"{", RowBox[{"b", ",", "c"}], "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419428190874147`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{"-", "b"}], "+", RowBox[{ SqrtBox["c"], " ", SqrtBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]], " ", SqrtBox[ FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ SuperscriptBox["b", "2"], "-", RowBox[{"4", " ", "c"}]}], ")"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "c"]], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"4", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"b", "+", RowBox[{ SqrtBox["c"], " ", SqrtBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]], " ", SqrtBox[ FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ SuperscriptBox["b", "2"], "-", RowBox[{"4", " ", "c"}]}], ")"}], " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "c"]], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"2", " ", "c", " ", "phase"}], "+", RowBox[{"2", " ", "c", " ", "n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{ SuperscriptBox["b", "2"], " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}], "-", RowBox[{"4", " ", "c", " ", "t", " ", RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]}]}], RowBox[{"2", " ", "c"}]], "]"}]}]}], RowBox[{"2", " ", "c"}]]}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{ RowBox[{ SqrtBox["2"], " ", "b", " ", SqrtBox[ RowBox[{"c", "[", "2", "]"}]]}], "+", RowBox[{"2", " ", SqrtBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]], " ", SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], ")"}]}]]}], "-", RowBox[{ SqrtBox["2"], " ", "b", " ", SqrtBox[ RowBox[{"c", "[", "2", "]"}]], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}], "2"]}]}], RowBox[{"b", " ", SqrtBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]], " ", SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], ")"}]}]]}]]}]}], ",", RowBox[{"c", "\[Rule]", FractionBox[ SuperscriptBox["b", "2"], "4"]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"w", "[", RowBox[{"n", ",", "t"}], "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{ RowBox[{ RowBox[{"-", SqrtBox["2"]}], " ", "b", " ", SqrtBox[ RowBox[{"c", "[", "2", "]"}]]}], "+", RowBox[{"2", " ", SqrtBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]], " ", SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], ")"}]}]]}], "+", RowBox[{ SqrtBox["2"], " ", "b", " ", SqrtBox[ RowBox[{"c", "[", "2", "]"}]], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"phase", "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}]}], "+", RowBox[{"t", " ", RowBox[{"c", "[", "2", "]"}]}]}], "]"}], "2"]}]}], RowBox[{"b", " ", SqrtBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}]], " ", SqrtBox[ RowBox[{ SuperscriptBox["b", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox[ RowBox[{"Tanh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], ")"}]}]]}]]}]}], ",", RowBox[{"c", "\[Rule]", FractionBox[ SuperscriptBox["b", "2"], "4"]}]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.4419428191811647`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "2"}], ")"}], " ", "Toda", " ", "lattices"}], " ", "*)"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"s", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], ",", " ", "x", ",", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"s", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "x", ",", "t"}], "]"}], " ", "-", " ", RowBox[{"2", "*", RowBox[{"s", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}]}], "+", RowBox[{"s", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "x", ",", "t"}], "]"}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"s", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", "+", " ", "1"}], ")"}]}]}], ",", "s", ",", "n", ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", "t"}], "}"}], ",", RowBox[{"{", "}"}]}], "]"}]], "Input"], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419428193530397`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"s", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"phase", " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}], " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"t", " ", SuperscriptBox[ RowBox[{"c", "[", "3", "]"}], "2"]}], "+", RowBox[{"x", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], RowBox[{"c", "[", "3", "]"}]], "]"}]}]}], RowBox[{"c", "[", "3", "]"}]]}], "}"}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419428193530397`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"DDESpecialSolutions", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], ",", " ", "t"}], "]"}], " ", "\[Equal]", " ", RowBox[{"v", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}]}], ",", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], ",", " ", "x"}], "]"}], " ", "\[Equal]", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "x", ",", "t"}], "]"}], " ", "-", " ", RowBox[{"2", "*", RowBox[{"u", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}]}], "+", RowBox[{"u", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "x", ",", "t"}], "]"}]}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], " ", "+", " ", "1"}], ")"}]}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"u", ",", "v"}], "}"}], ",", "n", ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", "t"}], "}"}], ",", RowBox[{"{", "}"}]}], "]"}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \ \[Rule] \\(\\(-1\\)\\)\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 0\\)\\)}\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" are being removed \ because they are (i) negative, (ii) contain freedom, (iii) fail to balance \ highest exponent terms from two different terms in the original system. If \ \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the transformation u -> 1/v \ may result in a system that DDESpecialSolutions can solve.\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419428194311647`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"DDESpecialSolutionsmSolver", "::", "\<\"remove\"\>"}], RowBox[{ ":", " "}], "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"The potential solutions \\\\\ \\\"\\\", \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\!\\({\\({\\(\\(\\(\\(m[1]\\)\\) \ \[Rule] \\(\\(-1\\)\\)\\)\\), \\(\\(\\(\\(m[2]\\)\\) \[Rule] 0\\)\\)}\\)}\\), \ \\\"MT\\\"]\\)\[NoBreak]\\!\\(\\*StyleBox[\\\"\\\\\\\" are being removed \ because they are (i) negative, (ii) contain freedom, (iii) fail to balance \ highest exponent terms from two different terms in the original system. If \ \\\\\\\\!\\\\\\\\(M\\\\\\\\_i\\\\\\\\) < 0, then the transformation u -> 1/v \ may result in a system that DDESpecialSolutions can solve.\\\\\\\"\\\", \ \\\"MT\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.4419428194311647`*^9}], Cell[BoxData["\<\"Numerically testing the solutions, please wait.\"\>"], \ "Print", CellChangeTimes->{3.4419428195092897`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "0"}], "]"}], " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"], " ", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"phase", " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}], " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"t", " ", SuperscriptBox[ RowBox[{"c", "[", "3", "]"}], "2"]}], "+", RowBox[{"x", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], RowBox[{"c", "[", "3", "]"}]], "]"}]}]}], RowBox[{"c", "[", "3", "]"}]]}], ",", RowBox[{ RowBox[{"v", "[", RowBox[{"n", ",", "x", ",", "t"}], "]"}], "\[Rule]", RowBox[{ RowBox[{"-", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"phase", " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}], " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"t", " ", SuperscriptBox[ RowBox[{"c", "[", "3", "]"}], "2"]}], "+", RowBox[{"x", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], RowBox[{"c", "[", "3", "]"}]], "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"Tanh", "[", FractionBox[ RowBox[{ RowBox[{"phase", " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"n", " ", RowBox[{"c", "[", "1", "]"}], " ", RowBox[{"c", "[", "3", "]"}]}], "+", RowBox[{"t", " ", SuperscriptBox[ RowBox[{"c", "[", "3", "]"}], "2"]}], "+", RowBox[{"x", " ", SuperscriptBox[ RowBox[{"Sinh", "[", RowBox[{"c", "[", "1", "]"}], "]"}], "2"]}]}], RowBox[{"c", "[", "3", "]"}]], "]"}]}], ")"}]}]}]}], "}"}], "}"}], "}"}]], "Output", CellChangeTimes->{3.4419428195249147`*^9}] }, Open ]] }, WindowSize->{971, 899}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (November 10, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{ "Info3441917609-4281003"->{ Cell[1777, 62, 471, 7, 58, "Print", CellTags->"Info3441917609-4281003"], Cell[2251, 71, 1582, 38, 77, "Print", CellTags->"Info3441917609-4281003"]} } *) (*CellTagsIndex CellTagsIndex->{ {"Info3441917609-4281003", 128756, 3735} } *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 388, 9, 72, "Input"], Cell[CellGroupData[{ Cell[980, 35, 217, 4, 23, "Print"], Cell[1200, 41, 225, 4, 23, "Print"], Cell[1428, 47, 207, 4, 23, "Print"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[1684, 57, 68, 1, 31, "Input"], Cell[CellGroupData[{ Cell[1777, 62, 471, 7, 58, "Print", CellTags->"Info3441917609-4281003"], Cell[2251, 71, 1582, 38, 77, "Print", CellTags->"Info3441917609-4281003"] }, Open ]] }, Open ]], Cell[3860, 113, 633, 13, 72, "Input"], Cell[4496, 128, 108, 2, 31, "Input"], Cell[CellGroupData[{ Cell[4629, 134, 1049, 33, 52, "Input"], Cell[5681, 169, 151, 2, 23, "Print"], Cell[5835, 173, 160, 4, 22, "Message"], Cell[5998, 179, 1297, 39, 50, "Output"] }, Open ]], Cell[7310, 221, 530, 14, 52, "Input"], Cell[CellGroupData[{ Cell[7865, 239, 1346, 41, 52, "Input"], Cell[CellGroupData[{ Cell[9236, 284, 159, 2, 23, "Print"], Cell[9398, 288, 1816, 54, 23, "Print"], Cell[11217, 344, 194, 2, 23, "Print"], Cell[11414, 348, 2437, 76, 48, "Print"], Cell[13854, 426, 227, 5, 23, "Print"], Cell[14084, 433, 166, 2, 23, "Print"] }, Open ]], Cell[14265, 438, 843, 14, 57, "Message"], Cell[15111, 454, 843, 14, 57, "Message"], Cell[CellGroupData[{ Cell[15979, 472, 288, 8, 23, "Print"], Cell[16270, 482, 264, 5, 23, "Print"], Cell[16537, 489, 163, 2, 23, "Print"], Cell[16703, 493, 145, 1, 23, "Print"], Cell[16851, 496, 908, 30, 26, "Print"], Cell[17762, 528, 521, 12, 23, "Print"], Cell[18286, 542, 137, 1, 23, "Print"], Cell[18426, 545, 4910, 146, 92, "Print"], Cell[23339, 693, 264, 5, 23, "Print"], Cell[23606, 700, 139, 1, 23, "Print"], Cell[23748, 703, 2763, 85, 99, "Print"], Cell[26514, 790, 233, 5, 23, "Print"], Cell[26750, 797, 147, 1, 23, "Print"], Cell[26900, 800, 161, 2, 23, "Print"], Cell[27064, 804, 3292, 96, 92, "Print"], Cell[30359, 902, 151, 2, 23, "Print"] }, Open ]], Cell[30525, 907, 160, 4, 22, "Message"], Cell[CellGroupData[{ Cell[30710, 915, 3292, 96, 92, "Print"], Cell[34005, 1013, 264, 5, 23, "Print"], Cell[34272, 1020, 102, 1, 23, "Print"], Cell[34377, 1023, 159, 2, 23, "Print"] }, Open ]], Cell[34551, 1028, 3293, 96, 99, "Output"] }, Open ]], Cell[37859, 1127, 437, 12, 52, "Input"], Cell[CellGroupData[{ Cell[38321, 1143, 1728, 54, 72, "Input"], Cell[40052, 1199, 830, 14, 57, "Message"], Cell[40885, 1215, 830, 14, 57, "Message"], Cell[CellGroupData[{ Cell[41740, 1233, 151, 2, 23, "Print"], Cell[41894, 1237, 152, 2, 23, "Print"] }, Open ]], Cell[42061, 1242, 1390, 40, 87, "Output"] }, Open ]], Cell[43466, 1285, 326, 7, 31, "Input"], Cell[CellGroupData[{ Cell[43817, 1296, 2171, 71, 77, "Input"], Cell[45991, 1369, 579, 11, 21, "Message"], Cell[46573, 1382, 151, 2, 23, "Print"], Cell[46727, 1386, 1390, 40, 87, "Output"] }, Open ]], Cell[48132, 1429, 97, 2, 31, "Input"], Cell[48232, 1433, 334, 9, 31, "Input"], Cell[CellGroupData[{ Cell[48591, 1446, 1349, 42, 52, "Input"], Cell[49943, 1490, 151, 2, 23, "Print"], Cell[50097, 1494, 232, 4, 60, "Output"] }, Open ]], Cell[50344, 1501, 415, 10, 52, "Input"], Cell[CellGroupData[{ Cell[50784, 1515, 1361, 42, 52, "Input"], Cell[52148, 1559, 3467, 103, 69, "Output"] }, Open ]], Cell[55630, 1665, 451, 12, 72, "Input"], Cell[CellGroupData[{ Cell[56106, 1681, 1650, 52, 52, "Input"], Cell[57759, 1735, 932, 15, 57, "Message"], Cell[58694, 1752, 945, 15, 57, "Message"], Cell[59642, 1769, 151, 2, 23, "Print"], Cell[59796, 1773, 274, 5, 22, "Message"], Cell[60073, 1780, 2349, 72, 43, "Print"], Cell[62425, 1854, 1206, 36, 30, "Output"] }, Open ]], Cell[63646, 1893, 432, 12, 52, "Input"], Cell[CellGroupData[{ Cell[64103, 1909, 1224, 39, 52, "Input"], Cell[65330, 1950, 151, 2, 23, "Print"], Cell[65484, 1954, 1254, 38, 30, "Output"] }, Open ]], Cell[66753, 1995, 99, 2, 31, "Input"], Cell[CellGroupData[{ Cell[66877, 2001, 2386, 73, 72, "Input"], Cell[69266, 2076, 1299, 20, 95, "Message"], Cell[70568, 2098, 1286, 20, 95, "Message"], Cell[71857, 2120, 151, 2, 23, "Print"], Cell[72011, 2124, 3077, 88, 148, "Output"] }, Open ]], Cell[75103, 2215, 179, 5, 31, "Input"], Cell[CellGroupData[{ Cell[75307, 2224, 2375, 74, 72, "Input"], Cell[77685, 2300, 1293, 20, 94, "Message"], Cell[78981, 2322, 1280, 20, 94, "Message"], Cell[80264, 2344, 151, 2, 23, "Print"], Cell[80418, 2348, 3059, 88, 148, "Output"] }, Open ]], Cell[83492, 2439, 205, 6, 31, "Input"], Cell[CellGroupData[{ Cell[83722, 2449, 2675, 82, 92, "Input"], Cell[86400, 2533, 1300, 20, 94, "Message"], Cell[87703, 2555, 1287, 20, 94, "Message"], Cell[88993, 2577, 151, 2, 23, "Print"], Cell[89147, 2581, 3103, 88, 148, "Output"] }, Open ]], Cell[92265, 2672, 180, 4, 31, "Input"], Cell[CellGroupData[{ Cell[92470, 2680, 785, 26, 31, "Input"], Cell[93258, 2708, 151, 2, 23, "Print"], Cell[93412, 2712, 1249, 36, 55, "Output"] }, Open ]], Cell[94676, 2751, 492, 12, 52, "Input"], Cell[CellGroupData[{ Cell[95193, 2767, 1859, 59, 52, "Input"], Cell[97055, 2828, 184, 2, 22, "Message"], Cell[97242, 2832, 184, 2, 22, "Message"], Cell[97429, 2836, 151, 2, 23, "Print"], Cell[97583, 2840, 11514, 304, 540, "Output"] }, Open ]], Cell[109112, 3147, 493, 12, 52, "Input"], Cell[CellGroupData[{ Cell[109630, 3163, 727, 24, 31, "Input"], Cell[110360, 3189, 127, 2, 23, "Print"], Cell[110490, 3193, 1119, 34, 30, "Output"] }, Open ]], Cell[111624, 3230, 429, 10, 52, "Input"], Cell[CellGroupData[{ Cell[112078, 3244, 890, 29, 31, "Input"], Cell[112971, 3275, 127, 2, 23, "Print"], Cell[113101, 3279, 6780, 195, 303, "Output"] }, Open ]], Cell[119896, 3477, 156, 5, 31, "Input"], Cell[CellGroupData[{ Cell[120077, 3486, 1004, 30, 52, "Input"], Cell[121084, 3518, 127, 2, 23, "Print"], Cell[121214, 3522, 1178, 34, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[122429, 3561, 1339, 41, 72, "Input"], Cell[123771, 3604, 819, 14, 57, "Message"], Cell[124593, 3620, 819, 14, 57, "Message"], Cell[125415, 3636, 127, 2, 23, "Print"], Cell[125545, 3640, 2867, 81, 110, "Output"] }, Open ]] } ] *) (* End of internal cache information *)