(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 18300, 630] NotebookOptionsPosition[ 16708, 571] NotebookOutlinePosition[ 17043, 586] CellTagsIndexPosition[ 17000, 583] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Lax Pairs for the KdV equation", "Section", CellChangeTimes->{{3.497058189082504*^9, 3.497058200071642*^9}}], Cell[BoxData[ StyleBox[ RowBox[{ "First", " ", "we", " ", "define", " ", "the", " ", "Hamiltonian", " ", RowBox[{"operator", ":"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058205405012*^9, 3.497058216760372*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"H", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"4", RowBox[{"D", "[", RowBox[{"#", ",", RowBox[{"{", RowBox[{"x", ",", "3"}], "}"}]}], "]"}]}], "-", RowBox[{"3", RowBox[{"(", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", "x", "]"}], "*", "#"}], ",", "x"}], "]"}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], "*", RowBox[{"D", "[", RowBox[{"#", ",", "x"}], "]"}]}]}], ")"}]}]}], ")"}], "&"}]}]], "Input", CellChangeTimes->{{3.497053961396209*^9, 3.497054059394131*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"4", " ", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"{", RowBox[{"x", ",", "3"}], "}"}]], "#1"}]}], "-", RowBox[{"3", " ", RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "x"], RowBox[{"(", RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", "#1"}], ")"}]}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SubscriptBox["\[PartialD]", "x"], "#1"}]}]}], ")"}]}]}], "&"}]], "Output", CellChangeTimes->{3.497053984993398*^9, 3.497054060132504*^9, 3.497058229863216*^9}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{ "Example", " ", "for", " ", "the", " ", "use", " ", "of", " ", "the", " ", RowBox[{"Hamiltonian", ":"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058233135986*^9, 3.497058238992157*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"H", "[", RowBox[{"x", "^", "4"}], "]"}]], "Input", CellChangeTimes->{{3.497053987725083*^9, 3.497053995124348*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"96", " ", "x"}], "-", RowBox[{"3", " ", RowBox[{"(", RowBox[{ RowBox[{"8", " ", SuperscriptBox["x", "3"], " ", RowBox[{"u", "[", "x", "]"}]}], "+", RowBox[{ SuperscriptBox["x", "4"], " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}]}]], "Output", CellChangeTimes->{{3.497053991502734*^9, 3.497053995356415*^9}, 3.497054062204892*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[StyleBox["Next, define the Schroedinger operator :", \ "Subsection"]], "Subsubtitle", CellChangeTimes->{{3.497058241768033*^9, 3.497058249204941*^9}, { 3.497058438761869*^9, 3.49705845087325*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"L", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"D", "[", RowBox[{"#", ",", RowBox[{"{", RowBox[{"x", ",", "2"}], "}"}]}], "]"}]}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], "*", "#"}]}], ")"}], "&"}]}]], "Input", CellChangeTimes->{{3.497054083453381*^9, 3.497054127845465*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"{", RowBox[{"x", ",", "2"}], "}"}]], "#1"}]}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", "#1"}]}], "&"}]], "Output", CellChangeTimes->{{3.497054101764336*^9, 3.497054128388726*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{"Example", ":"}], "Subsection"]], "Input", CellChangeTimes->{{3.497058265179659*^9, 3.497058266824391*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"L", "[", RowBox[{"x", "^", "2"}], "]"}]], "Input", CellChangeTimes->{{3.497054103631685*^9, 3.497054119606182*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", "2"}], "+", RowBox[{ SuperscriptBox["x", "2"], " ", RowBox[{"u", "[", "x", "]"}]}]}]], "Output", CellChangeTimes->{{3.497054104803555*^9, 3.497054129474798*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{ "Applying", " ", "the", " ", "L", " ", "to", " ", "the", " ", "H", " ", "to", " ", "an", " ", "arbitrary", " ", RowBox[{"function", ":"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058285108891*^9, 3.497058299016789*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"L", "@", RowBox[{"H", "@", RowBox[{"f", "[", "x", "]"}]}]}]], "Input", CellChangeTimes->{{3.497054132923225*^9, 3.497054142319487*^9}, { 3.497054175872174*^9, 3.497054195932884*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "3"}], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "+", RowBox[{"3", " ", RowBox[{"(", RowBox[{ RowBox[{"5", " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{"2", " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "-", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "5", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}]], "Output", CellChangeTimes->{{3.497054165934281*^9, 3.497054196394765*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{ "Applying", " ", "the", " ", "H", " ", "to", " ", "the", " ", "L", " ", "and", " ", "to", " ", "an", " ", "arbitrary", " ", RowBox[{"function", ":"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058306486685*^9, 3.497058311838708*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"H", "@", RowBox[{"L", "@", RowBox[{"f", "[", "x", "]"}]}]}]], "Input", CellChangeTimes->{{3.4970542053409*^9, 3.497054208483954*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", "3"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{"u", "[", "x", "]"}]}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], ")"}]}], "+", RowBox[{"2", " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "-", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], ")"}]}]}], ")"}]}], "+", RowBox[{"4", " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], "-", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "5", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], ")"}]}]}]], "Output", CellChangeTimes->{3.497054210031022*^9}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{"The", " ", "commutator"}], ",", " ", RowBox[{"[", RowBox[{"L", ",", "H"}], "]"}], ",", " ", RowBox[{ "is", " ", "the", " ", "difference", " ", "between", " ", "the", " ", "two", " ", "operators", " ", "above"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058318063523*^9, 3.497058350286309*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{{3.49705421299067*^9, 3.49705422500985*^9}, { 3.497058355740602*^9, 3.497058356799909*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"L", "@", RowBox[{"H", "@", RowBox[{"f", "[", "x", "]"}]}]}], "-", RowBox[{"H", "@", RowBox[{"L", "@", RowBox[{"f", "[", "x", "]"}]}]}]}]], "Input", CellChangeTimes->{{3.497058361941091*^9, 3.497058372572323*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"3", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{"u", "[", "x", "]"}]}], "-", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], ")"}]}], "+", RowBox[{"2", " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "-", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], ")"}]}]}], ")"}]}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "3"}], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "+", RowBox[{"3", " ", RowBox[{"(", RowBox[{ RowBox[{"5", " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{"2", " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}], ")"}]}], "-", RowBox[{"4", " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], "+", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], "-", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "5", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}], ")"}]}], "-", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", TagBox[ RowBox[{"(", "5", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}]], "Output", CellChangeTimes->{{3.497054218997294*^9, 3.497054225401682*^9}, { 3.497058373701224*^9, 3.497058388118151*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{"Looks", " ", "bad"}], ",", " ", RowBox[{"but", ":"}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058376364333*^9, 3.497058379130985*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"%", "//", "Expand"}]], "Input", CellChangeTimes->{{3.49705838039754*^9, 3.497058390175203*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"6", " ", RowBox[{"f", "[", "x", "]"}], " ", RowBox[{"u", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", "\[Prime]", MultilineFunction->None], "[", "x", "]"}]}], "-", RowBox[{ RowBox[{"f", "[", "x", "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[ RowBox[{"(", "3", ")"}], Derivative], MultilineFunction->None], "[", "x", "]"}]}]}]], "Output", CellChangeTimes->{{3.497058384615912*^9, 3.497058390487699*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{"which", " ", "is", " ", RowBox[{"f", "[", "x", "]"}], " ", "times", " ", "the", " ", "KdV", " ", RowBox[{"equation", ".", " "}]}], "Subsection"]], "Input", CellChangeTimes->{{3.497058397771367*^9, 3.497058418804714*^9}}] }, Open ]] }, WindowSize->{640, 750}, WindowMargins->{{172, Automatic}, {Automatic, 77}}, FrontEndVersion->"6.0 for Linux x86 (32-bit) (April 20, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 115, 1, 73, "Section"], Cell[708, 26, 233, 5, 34, "Input"], Cell[CellGroupData[{ Cell[966, 35, 664, 22, 30, "Input"], Cell[1633, 59, 644, 22, 34, "Output"] }, Open ]], Cell[2292, 84, 239, 5, 34, "Input"], Cell[CellGroupData[{ Cell[2556, 93, 143, 3, 30, "Input"], Cell[2702, 98, 489, 15, 37, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[3240, 119, 213, 3, 29, "Subsubtitle"], Cell[CellGroupData[{ Cell[3478, 126, 380, 12, 30, "Input"], Cell[3861, 140, 319, 10, 34, "Output"] }, Open ]], Cell[4195, 153, 143, 3, 34, "Input"], Cell[CellGroupData[{ Cell[4363, 160, 143, 3, 30, "Input"], Cell[4509, 165, 211, 6, 33, "Output"] }, Open ]], Cell[4735, 174, 277, 6, 34, "Input"], Cell[CellGroupData[{ Cell[5037, 184, 216, 5, 30, "Input"], Cell[5256, 191, 2077, 67, 86, "Output"] }, Open ]], Cell[7348, 261, 289, 6, 34, "Input"], Cell[CellGroupData[{ Cell[7662, 271, 165, 4, 30, "Input"], Cell[7830, 277, 2427, 78, 114, "Output"] }, Open ]], Cell[10272, 358, 373, 9, 59, "Input"], Cell[10648, 369, 139, 2, 30, "Input"], Cell[CellGroupData[{ Cell[10812, 375, 269, 8, 30, "Input"], Cell[11084, 385, 4462, 143, 166, "Output"] }, Open ]], Cell[15561, 531, 198, 5, 34, "Input"], Cell[CellGroupData[{ Cell[15784, 540, 120, 2, 30, "Input"], Cell[15907, 544, 505, 16, 34, "Output"] }, Open ]], Cell[16427, 563, 265, 5, 34, "Input"] }, Open ]] } ] *) (* End of internal cache information *)