(* 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[     94174,       2326]
NotebookOptionsPosition[     87505,       2107]
NotebookOutlinePosition[     88214,       2132]
CellTagsIndexPosition[     88171,       2129]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Lecture 2 - Quantum mechanics of identical bosons", "Title",
 CellChangeTimes->{{3.4165984729416*^9, 3.4165984729572*^9}}],

Cell["by Leonid Pryadko", "Subtitle"],

Cell[CellGroupData[{

Cell["Initialization", "Section",
 InitializationCell->True],

Cell[CellGroupData[{

Cell["Pauli Matrices", "Subsubsection",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SubscriptBox["\[Sigma]", "x"], "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "0"}], "}"}]}], "}"}]}], ";", 
  RowBox[{
   SubscriptBox["\[Sigma]", "y"], "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"-", "I"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"I", ",", "0"}], "}"}]}], "}"}]}], ";", 
  RowBox[{
   SubscriptBox["\[Sigma]", "z"], "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"1", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], ";", " ", 
  RowBox[{"e", "=", 
   RowBox[{"IdentityMatrix", "[", "2", "]"}]}], ";"}]], "Input",
 InitializationCell->True]
}, Open  ]],

Cell[CellGroupData[{

Cell["Direct product", "Subsubsection",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"Unprotect", "[", "NonCommutativeMultiply", "]"}], ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"NonCommutativeMultiply", "[", 
    RowBox[{
     RowBox[{"A_", "?", "MatrixQ"}], ",", 
     RowBox[{"B_", "?", "MatrixQ"}]}], "]"}], ":=", 
   RowBox[{"KroneckerProduct", "[", 
    RowBox[{"A", ",", "B"}], "]"}]}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4166025341726*^9, 3.416602538135*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Spin-1/2 spin matrices", "Subsubsection",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"sx", "=", 
   RowBox[{
    SubscriptBox["\[Sigma]", "x"], "/", "2"}]}], ";", 
  RowBox[{"sy", "=", 
   RowBox[{
    SubscriptBox["\[Sigma]", "y"], "/", "2"}]}], ";", " ", 
  RowBox[{"sz", "=", 
   RowBox[{
    SubscriptBox["\[Sigma]", "z"], "/", "2"}]}], ";"}]], "Input",
 InitializationCell->True]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Identical particles: bosons", "Section"],

Cell["Define truncated Hilbert space for bosons", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "7"}], ";"}]], "Input"],

Cell["Define boson creation and annihilation operators", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"bp", ":=", 
   RowBox[{"Block", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", "j"}], "}"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"j", "\[Equal]", 
          RowBox[{"i", "-", "1"}]}], ",", 
         RowBox[{"Sqrt", "[", "j", "]"}], ",", "0"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "NN"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "NN"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "5"}], ";", 
  RowBox[{"bp", "//", "MatrixForm"}]}]], "Input",
 CellChangeTimes->{{3.4166017214906*^9, 3.416601733487*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "0", "0", "0"},
     {"1", "0", "0", "0", "0"},
     {"0", 
      SqrtBox["2"], "0", "0", "0"},
     {"0", "0", 
      SqrtBox["3"], "0", "0"},
     {"0", "0", "0", "2", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.416601726311*^9, 3.4166017345478*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"bp", "//", "MatrixForm"}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "0", "0", "0"},
     {"1", "0", "0", "0", "0"},
     {"0", 
      SqrtBox["2"], "0", "0", "0"},
     {"0", "0", 
      SqrtBox["3"], "0", "0"},
     {"0", "0", "0", "2", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.4166018008322*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"bm", ":=", 
   RowBox[{"Transpose", "[", "bp", "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"bm", ".", "bp"}], "//", "MatrixForm"}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0", "0", "0"},
     {"0", "2", "0", "0", "0"},
     {"0", "0", "3", "0", "0"},
     {"0", "0", "0", "4", "0"},
     {"0", "0", "0", "0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.4166018206286*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"bp", ".", "bm"}], "//", "MatrixForm"}], 
  "\[IndentingNewLine]"}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "0", "0", "0"},
     {"0", "1", "0", "0", "0"},
     {"0", "0", "2", "0", "0"},
     {"0", "0", "0", "3", "0"},
     {"0", "0", "0", "0", "4"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.4166018441066*^9}]
}, Open  ]],

Cell["Commutator", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"bm", ".", "bp"}], "-", 
   RowBox[{"bp", ".", "bm"}]}], "//", "MatrixForm"}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0", "0", "0"},
     {"0", "1", "0", "0", "0"},
     {"0", "0", "1", "0", "0"},
     {"0", "0", "0", "1", "0"},
     {"0", "0", "0", "0", 
      RowBox[{"-", "4"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.4166018482406*^9}]
}, Open  ]],

Cell["This is the coordinate operator", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"xx", ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"bp", "+", "bm"}], ")"}], "/", 
    RowBox[{"Sqrt", "[", "2", "]"}]}]}], ";", 
  RowBox[{"xx", "//", "MatrixForm"}]}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", 
      FractionBox["1", 
       SqrtBox["2"]], "0", "0", "0"},
     {
      FractionBox["1", 
       SqrtBox["2"]], "0", "1", "0", "0"},
     {"0", "1", "0", 
      SqrtBox[
       FractionBox["3", "2"]], "0"},
     {"0", "0", 
      SqrtBox[
       FractionBox["3", "2"]], "0", 
      SqrtBox["2"]},
     {"0", "0", "0", 
      SqrtBox["2"], "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.4166019348206*^9}]
}, Open  ]],

Cell["Operator for non-linear oscillator", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"be", ":=", 
   RowBox[{"IdentityMatrix", "[", "NN", "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"be", "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.4166019604514*^9, 3.4166019627757998`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0", "0", "0"},
     {"0", "1", "0", "0", "0"},
     {"0", "0", "1", "0", "0"},
     {"0", "0", "0", "1", "0"},
     {"0", "0", "0", "0", "1"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.416601963197*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"H", ":=", 
   RowBox[{
    RowBox[{"bp", ".", "bm"}], "+", 
    RowBox[{"0.5", " ", "be"}], " ", "+", 
    RowBox[{"0.1", " ", 
     RowBox[{"(", 
      RowBox[{"xx", ".", "xx", ".", "xx", ".", "xx"}], ")"}]}]}]}], ";", " ", 
  
  RowBox[{"H", "//", "MatrixForm"}]}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.575`", "0", "0.21213203435596428`", "0", "0.1224744871391589`"},
     {"0", "1.875`", "0", "0.6123724356957946`", "0"},
     {"0.21213203435596428`", "0", "3.475`", "0", "0.7794228634059949`"},
     {"0", "0.6123724356957945`", "0", "4.875`", "0"},
     {"0.1224744871391589`", "0", "0.7794228634059949`", "0", "5.2`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.4166019967525997`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "17"}], ";", 
  RowBox[{"H", "//", "MatrixForm"}]}]], "Input"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.575`", "0", "0.21213203435596428`", "0", "0.1224744871391589`", "0", 
      "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"},
     {"0", "1.875`", "0", "0.6123724356957946`", "0", "0.2738612787525831`", 
      "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0"},
     {"0.21213203435596428`", "0", "3.475`", "0", "1.2124355652982142`", "0", 
      "0.474341649025257`", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "0"},
     {"0", "0.6123724356957945`", "0", "5.375`", "0", "2.012461179749811`", 
      "0", "0.724568837309472`", "0", "0", "0", "0", "0", "0", "0", "0", 
      "0"},
     {"0.1224744871391589`", "0", "1.2124355652982142`", "0", "7.575`", "0", 
      "3.012474066278414`", "0", "1.02469507659596`", "0", "0", "0", "0", "0",
       "0", "0", "0"},
     {"0", "0.2738612787525831`", "0", "2.012461179749811`", "0", "10.075`", 
      "0", "4.212481453965109`", "0", "1.3747727084867523`", "0", "0", "0", 
      "0", "0", "0", "0"},
     {"0", "0", "0.474341649025257`", "0", "3.012474066278414`", "0", 
      "12.875`", "0", "5.612486080160912`", "0", "1.774823934929885`", "0", 
      "0", "0", "0", "0", "0"},
     {"0", "0", "0", "0.724568837309472`", "0", "4.212481453965109`", "0", 
      "15.975`", "0", "7.212489168102786`", "0", "2.224859546128699`", "0", 
      "0", "0", "0", "0"},
     {"0", "0", "0", "0", "1.02469507659596`", "0", "5.612486080160912`", "0",
       "19.375`", "0", "9.012491331479882`", "0", "2.724885318687743`", "0", 
      "0", "0", "0"},
     {"0", "0", "0", "0", "0", "1.3747727084867523`", "0", 
      "7.212489168102785`", "0", "23.075000000000003`", "0", 
      "11.012492905786592`", "0", "3.2749045787625635`", "0", "0", "0"},
     {"0", "0", "0", "0", "0", "0", "1.774823934929885`", "0", 
      "9.012491331479882`", "0", "27.075`", "0", "13.212494087037465`", "0", 
      "3.8749193539995126`", "0", "0"},
     {"0", "0", "0", "0", "0", "0", "0", "2.224859546128699`", "0", 
      "11.012492905786592`", "0", "31.375`", "0", "15.612494995995995`", "0", 
      "4.524930938699508`", "0"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "2.724885318687743`", "0", 
      "13.212494087037467`", "0", "35.975`", "0", "18.212495710363257`", "0", 
      "5.224940191045253`"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "3.2749045787625635`", "0",
       "15.612494995995995`", "0", "40.875`", "0", "21.012496281974688`", 
      "0"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "3.8749193539995126`",
       "0", "18.212495710363253`", "0", "46.075`", "0", 
      "17.42842505793338`"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "4.524930938699508`", "0", "21.012496281974688`", "0", 
      "44.775000000000006`", "0"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "5.224940191045253`", "0", "17.428425057933378`", "0", "28.9`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.416602005551*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "7"}], ";", 
  RowBox[{
   RowBox[{"Eigenvalues", "[", "H", "]"}], "//", "Sort"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.5591062583247673`", ",", "1.7702171302575584`", ",", 
   "3.152810708294379`", ",", "4.557132332061639`", ",", "5.92289958870426`", 
   ",", "9.947650537680802`", ",", "10.140183444676596`"}], "}"}]], "Output",
 CellChangeTimes->{3.4166020402922*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "5"}], ";", 
  RowBox[{
   RowBox[{"Eigenvalues", "[", "H", "]"}], "//", "Sort"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.558593945429275`", ",", "1.7548148253980358`", ",", 
   "3.1840332804520743`", ",", "4.995185174601969`", ",", 
   "5.507372774118651`"}], "}"}]], "Output",
 CellChangeTimes->{3.4166020443482*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "35"}], ";", 
  RowBox[{"vals", "=", 
   RowBox[{
    RowBox[{"Eigenvalues", "[", "H", "]"}], "//", "Sort"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.5591463271835034`", ",", "1.7695026439490402`", ",", 
   "3.138624308497262`", ",", "4.628882808877491`", ",", "6.220300900045619`",
    ",", "7.899767229388651`", ",", "9.657839998004865`", ",", 
   "11.487315483715397`", ",", "13.382474187197824`", ",", 
   "15.33864496887862`", ",", "17.351935572372444`", ",", 
   "19.418825935853622`", ",", "21.53615176908419`", ",", 
   "23.70591852129439`", ",", "25.923383861078623`", ",", 
   "28.094993450046115`", ",", "30.552241678252503`", ",", 
   "33.04947996486148`", ",", "34.18322065159492`", ",", "40.95331747485684`",
    ",", "41.1479015545248`", ",", "53.060210641768116`", ",", 
   "53.13082771919421`", ",", "70.30718553948448`", ",", "70.34688055745814`",
    ",", "94.27442442632903`", ",", "94.30116880230621`", ",", 
   "127.44250124769455`", ",", "127.46219133491805`", ",", 
   "173.86849124343078`", ",", "173.88369097050193`", ",", 
   "241.1000939816919`", ",", "241.11210054183823`", ",", 
   "347.7254444378794`", ",", "347.73491926594664`"}], "}"}]], "Output",
 CellChangeTimes->{3.416602074191*^9}]
}, Open  ]],

Cell["Compare with the solution of the differential equation", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"solution", "=", 
   RowBox[{"NDSolve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"-", 
           RowBox[{
            SuperscriptBox["u", "\[DoublePrime]",
             MultilineFunction->None], "[", "x", "]"}]}], "/", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"x", "^", "2"}], "/", "2"}], "+", 
            RowBox[{"0.1", " ", 
             RowBox[{"x", "^", "4"}]}]}], ")"}], " ", 
          RowBox[{"u", "[", "x", "]"}]}]}], "\[Equal]", 
        RowBox[{"0.559", 
         RowBox[{"(*", " ", 
          RowBox[{"vals", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "*)"}], 
         RowBox[{"u", "[", "x", "]"}]}]}], ",", 
       RowBox[{
        RowBox[{
         SuperscriptBox["u", "\[Prime]"], "[", "0", "]"}], "==", "0"}], ",", 
       RowBox[{
        RowBox[{"u", "[", "0", "]"}], "\[Equal]", "1"}]}], "}"}], ",", "u", 
     ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "4.5"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.4166021997398*^9, 3.4166022304094*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"u", "[", "x", "]"}], "/.", "\[InvisibleSpace]", "solution"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "4.5`"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0", ",", "1"}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4166021461226*^9, 3.4166021504282*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwdl3k4ld/Xxo/hHA4HKUOokEhIiahvsVZEJIUiScbKmLk0IsMxZ0zGyJAh
Y4QUIULIVCFURGYHESW85/fuf57rc61nr+d+1r7XvvYWtXTSv8JIIBBeMhAI
/3suuohPrzy9qfJp43+DBpdDtpsKJRsCc4LtmzU6C5RIb5pMugrUX9aUZTqj
bsah9CVn8DEWGZ6i8+K7A7adSTegT+AhDtA5W7k2fkPtLvjkftjdRGfZL7tH
Y5e8oSMz3LaQzoG8Tmt7m32hsHnNPIrOQ7plvA1JVHCSKhx2pvN/oWt7TZyD
4MuE1sxJOsc0Htf4pRYK/AL93iJ0nmUMNQ3mDwdeNffk+XUa7EqxfRWxFAHN
rWJq1XTuPBp2wsQsClzurblQ6ez5pah7d3M0yNjwSZ2kcy/vylR1UiwsJhgr
1q/RwL9E6GYwKQ6KbpcH36bzAT1gNnSOB1FKiY0MnUNDqdtm1JJAdQ+xzucf
DQ5L5WZX5CfDK6kk5d10/tnYpuDHnwIXtc+qNa7SIPV+8NOwpVR4vPxs5s9f
GtzQTXJwy30C3XYPo0LprCNccOCCWRqMb9tXJ0jnv687q3c1pwNx+5a6PX9o
cHZla8+rpKdw+NkVrc3LNJBqlEp+opcFRhbhmvd/04AQe9QqgJQNv92/dk8t
0SBP3pym75wDEcVvw0sWaUByymKZVMsD867ONI4FGgwov2xrX8mDP2IG5Avz
NCihtES/yM+HyayBtZQ5Gpjnzgrf5y+EtWauKBEaDQYXBLeFLhXBwNhvx+kp
Gmh+EdtYP1MMDqqeOjx0Lq2VGXbJLYaSVIPXSpM0CI5QyT5v9hw8gqRPuYzT
QEnWUkGsuQSGi8OVC0doEGmbc6oyqQxsd+a9MBygAadzel/Tehl8aLKzF+6n
1/dG8tUe83JIsJWaH+6j+8M30ntxVwUsiTY7GvfQwDv5Vqls/ktIf9x0jL2L
Bo6dWtvSq16Dlv/lJvd3NJjuUcspFq6Cj+/21U3W08Duq7Jizf0q0M9wfnjx
LQ2uTsnpDqpXg1/S6vu9NTS4RBT04//wBhYn9PxjXtLg6balJXWnWnA/oNxj
/4wGo8KuR3VSa2FTUfPr8By6v8RoPuc6a6Fwt1NAQRYN0vZMcFnJ10HmA6X2
gXQaJCsO7vFergM+g1T4l0j3q27DpUqvelBiVbM8EUIDX7+HDbIRjbCnLnYL
71UaOHikaP+uaoQTB/dX1VrRwNAup7NqqhGMBZQ47Czo661bNXjqRBPYHCBl
5JvQoENodMluvQnE5PUaCGdpsKNEQSLL4T0UfyI+6VGhQcVQN1VYqw1qeQYs
T2+h6/s4SBm70Qau3+8/ddlEr2/jWFRBRhuYbDEpjOCggVn+6mNlwgfwZqPx
1LHQ/XJLrPxi+QfgiC7QWVidBanAl79aNneArVugfdCPWfAQJu9ivd8JwtH5
hgmFs3BVmnvp5qNOuPefs5tW3iwYKgm8m8jvhIMdU2mL2bOgcGaPTcuXTgha
qFeBtFlY8NTKe6DQBaoMp0IyY2bB4WuQAu94F3TcvkwtvzkL5sls6jvPfATl
Cv03f5Rn4aQQ5erR7T0wx8154VjVDKw7RojGHOyBoR/phWkVM/C8jmdwSqcH
eP2kWzdKZkDQbvu5RE963PuFaX7uDExWyKr+/dYD1slRwp8fzUCQod6OirRe
OPQoqN7MZQaaImM/y0t+gRLugCgL0RnQYBU9IS0/CAn/mQrGXJ8GnaUJSqTz
EMRSKmZFOKeAYtbaYOs3AnZm7m/dzk+A4FZuGAkagyiyTMfN+2PwsWVXseXF
SVgv9fzTWjwKYjdrXattZ8D/kLGxee4PCA+P5ePaOQey/K25LW+GoDR4rZsr
cB4MtnBZjNV/AzHGPYe9KL/A1kEo7Vp3P4z/fpUwd2MRfJHJ9o5OLyxbb8no
2FgC/Vba294XH6GAkBR2//4yUOcuicD1TnjFF+qmMLECGXzh08lX2mBzDjIY
Ov0FweZ3t/LYm+BNhUnIrf5V+O6uv/Ql+y2UH1euYr+6BiYLwRZKIm+AGuKv
JeCwBjpm/E7PGd/Aua62TxKua3CsvaDg2Gg1zJmZzh7zXAMXbqXipNxqkLrj
JXwzdg14Qm+pHFeshsfFtfdH3q1B0TFrRQO9KqDu0NCo2r0OByVN5oxjX4Hh
yul2x4l16PjvcEqNcQUIDgZa8dDWIS6Wd1ODWgV8ra1bfrm4DnafN33tk6mA
KyGKIkTCBqhrx91V2CgH9x07XBL5NiCS314uP6Mcoo/PbG5S24CoHR69WQtl
0BERYij6eAMkOFaOJiW8gEyekVnbPQTk/KZgriFQAoQI3dhbMgT0KNzkk8Rc
AsbsVcpB+whYMmilyDz3HLgYY0OzDxLwbmq4Emvjc7hFOyE1doyAouo8b4eu
Pwed98+uXDYm4M4sf5mPn4thycttwDSEgHxh9lu8MotAfYrx/dkZAhrovs42
GM+DhIM21/fNEXDqdFj/veY8mPNqE2H/RcAam+F9tbn0OE+cR90KAY+pbWeJ
v5YHs8p7xeWYGVCuK34pZPEZxEWc9+YUYsCbB7b/c2V5BpMH8xSbNRkw8ZFn
7U/MAfDmHs7QZsBNF46MGYnnQMz7G2HepxlQY6/f3REyPW52bETpHANO3TB6
c6I7G2ICP0c+NWPAJ6y2e5Sts0G5nzDje4MBcxikn/VEZ0GEt0G6SgYD/vMp
tXbYyASFvXfm2LMYsC9ukenaWCb09aUq9+UwoPbzJhe/9kzYKT/d41bIgP6p
RuHLKZlQOupDyalkwE6nVMeZY5nQe7Lo+pZOBiQblmZRgjNAmJddc3yNAbfq
WbyukUqH+tr9D18QGFEvx5fLlzcdbBwNh32YGPHn/OUnxoR0KGp8cnc7mRGf
3ik+o/U5DVRvHyrW52FEX+0jd8R80+DqtyuCVVKM2EB0Uan5/gQKcmqmI88z
4i15sWz1Z6mg/e0gf5Ix/f0nPlJMj1JhnOfZsaeXGNHL8uD3Hp9U2On9MLbS
ihGTq12z3xmnQqyhHf5wYsTPlQUMFeyp4MXEE6MQyIjp9yX3vvzxGHQvWR/p
qWDE56dWICg5CWaiBq4MvWLEs4MzXbG3kyC4SS9iqpoRAwUCtrSfT4J6+aOj
G/WMmHhPW+rN5iQ4xMYdLtnJiBe7xIMTghNBtLxy+NYEIx4uJ+4R90mAhU2c
wdsEmdDr7k2Bd+FxEMx7t2JhGxOWJQw9ornFwU7ByZ9NwkxoE2py5IhRHOiJ
vVO7Ls6Ez22ZEw+J0s8tBz3X2vYz4VzJW8vm0kfgaDzr7HWCPv81p9v+4ViY
TP9wftidCVWP8rT5GDwEn+yj1AoPJnw4tyzVpPwQBPNzSx/cZsIe4hj/XomH
oFVG5T7izYTfOsN1DZZjILtJpSU6hAmVPiRu254YA1dnClTU05hw1MDZ6NZ4
NAwphotntzOhefae9aH4KHCtPyX3qosJBQo8r530iwJGfTblD5+YsOE6l9In
xyjYec3/3GI/E/KdvX1f73gUXE6/44vjTLgY8PPh8Fwk/OSyGeolMKNmt7Pc
Od1ImBrHZPIBZhzkSNs6tjsC7t5Yz952kBmnz14vdOWLAArz69J9h5hR3QyE
dxAjQEZYqdVAhRlfRX30zW4IB0eDvatPTjKj/Mc6TQ8Ih7lagQv/WTFjQdmR
3HaVB7CUsMBjH8OM269nTgfZh8LHYOc6p0fMGOQvEPZJNxRKb886uScwo3Q+
pc5YMRTcjSff30tlxtQBRz9uxlBYEBj2jshjxqF+toq9SSFAi+ucflHPjESi
jiKlPxg+BJ5JqGxkxnSWpaHRt8FQcLPtxJv3zBifuFq4nBcMjkbNaU0dzJjx
Qd7irVcwTPPXGvUPMON3jmvTK7uDYSK2qJ6wxIxXFoIMWn2CoIkq60pcYUar
iUFFYacgyL6RJ8y2yoxRyXpbcy4GgY1h9u0tDEQ8n1r+2OJgEPzkTd0vwUHE
ZLf3Oy9PBcKPmPAkbXEiKuyv9Ui3CIQ6P46TupJELJmzmDfTC4Q095Dlc9JE
lLl1LFX/WCBYngvQN5UjokXGe48vooHwfYsXq4syEcdu9wq4jQbAYJSje6wB
ERN3aLbW3giACr1RcWsjIhZQbqhX2gVA1CaTHqWLRKzvDwgeMg0AjQda//WZ
E1FeJfza6xMBUBS0i7DNgYjVh6ZM0oUCIFgzsXjakYjNUfJ/uLkD4DLLZqsq
F/r3mOXPFZICQMB/o8HUg67357D6jQUq+Hh/CXniQ8SP4q4SB1upYAJ6yq7+
RJzlqKoXeEsFxfXGWdVAIlacFtqzt5IKU3de6I2EEZHbZknnezYVznlE8O+O
J+ILNWvzm4FUkFVkaV5OJCL/uJp/tzcVWJfu3W56TMTIhenQS7eoUOVqP2ib
Qc9XUHph3Y4KEo4a6XmFRKT8xx3ppksFwt6qc/eeE1HgUZOLxkkqfJmSJ51+
QWd9iq3mcSqE24ra0iqJqPZc83XfISrYScYJ1lQR8WVD+R5beSocH+Nsjagh
4oe2xI7dslRYufxP9sA7Il6tfD5/YBcVusTcvjM2E/G6rra3lzAV8oYnIrtb
6PnJEi6rglSwMO9ZdO+k6/eMfXp/MxWOCp9+qv6RiFxyyYyhnFTg+1p/nq+H
iMXVzbQGNiq8v/i8snyAiBob74qamKiQIbjHIfAbEQcS96dEEqjg2Zey/cIw
EVWLf50NXvMHozi+9j2jdL0XzQef//UH+fNh3n/HiMgjt/8k+4o/cPIxH2iZ
JOIuiYCc6CV/GP94+0fiDN1P7otsp375Q130fIzDHBFPebp5Ksz7Q5K+jYby
LyLOlznw6ND84Qb3t2WO30TU9i8bfDjjD3odBjlfV4hYJhU7xTntDzLhrcaF
q/R82TVnKib9gXRajeK9To+nX9wWNeEP3ymVVboMJBzQWbdKGveHypb9TqLM
JOQ+8fxwz5g/xARniSyQSCjhv1ikSmcnrR1ddWQSPlYtG/z60x9UI0fMOygk
JLtH1+bQmacvlzbIRULZ1gKRAjqPibh4Tm0moYiyOt80nSttlDj+8JIwnLKv
2oqeL6xoLZEkQMKOsCH6OdcfzFfeSvFsI+Fsv/MjAl2fPAa/FBUmYVl0JJPE
//QH6mru20nCN2EJj/yn/KGvna/nqDgJ78RRXYXp/5/HP3jlpCQJjx/4XfJr
1h+8zNIXz0uT0FJN7i4jvX76Wba+V2RJuKNXclqbXl9x2r7NbnIkbP1yU7yD
Xv8Vxd+p3gok7G2zPRlOX5/kdz7Vif+R8PqPlz+q1/3BmVNLJ0eZhHVhW1P3
MVJB1ZBroAxJ+O+ZYG8/kQrjo4l/OjVISHrToD/GQYVXey0Dv2mRUIVz6w8N
ur8eXJfknzlFQu2/VsHDfFRQIL5QYD1LQum3x2s/iFCBRefOW15DEqZwXs4Q
lqD3Q8wxfbELJDym+VA9X5oK3rs+OKmYkVDrRYqLlxIVWtTGct3tSfilzN+a
WZ8Kj0PyD/s4kvDPwxwZvEAF5263pnAXEvLO2+ysNKf72YrwM9eDhCXFZ61t
nalg6SMgOuRLwtJWoZXKKCocfP+tcJZKQk8qx9eZRHo/b36q8i+IhMGXg+LM
M6lQ8OSACX8ECTXLXziQXlLhT432I50kEsqpSc1c+EGFVtbNEhdTSPhbf7g1
cpauR7e31CaNhI5yv43Z/9L79fvlLt9sEhqu7/Cup+9XERv3OCpLSWi++/EW
w+MBsFulyFeijYSJFwz8D5QFwOI5Jz5CBwkpVOW6xPoAqLOXzenrImGVfRdZ
tTsATOPzPoT2kjCBQqvQmA+AmMVsgV8/SHjNnaFOaX8gEJ49KXzzl4TdpTbu
iy8D4UOtuWr8Ggkr3Y+v67QGQmKv8CdXAgtWb62gjX8NBEXS47/iJBbMVGxV
HmEOAgeLBPVQbhbM9WSX+HM2CHr5o/uNJFlQa5sB/8h6EDyV1Xc8IM2CPw3F
t97mDQY3dW4GiiwLMq3efGUoEwycbuESb+RZMHZC2HjUOBiOfwhxEQcWjE4R
UrN8FQzFfn4sCwb0+f5sjvuDQ8AzUS2hxYgFL7N9HLiUGQKnnjPuzbzIgi3F
9yR6a0Jg/Ju3vpEFC8Jr/vXSlRDYceRecrUDC/o63x8XcgiF4Pnr8iG+LJjz
OXWIyTIMsuyWrA2pLPjodIBl5P0wqP/hniQaxILzU5ktTk/CYP2TG7HiAQv2
KYcdFh0OA5dKl56ReBZ0ExxetbJ+AEa+1+6oFLHg60DBsd774XD933QhuYSu
Z4WZq+hpOERddxj5+IIFB07bHvneEg4t1vY69q9YUM7j0tRT9ghQOWUrHPeO
zmzfphQvRIA475W38wMsWPAiql50LQJUH4wsv/7Ggk3Hg0k3N0eCKctlmcBh
Fjxrw0zd2B0JcX8sH+4YZ8GxWhlq+tlIoHw1t9b+xYJlYYmlmvmRsPDUhD2T
zIraB5zrztpFQc2hc2eNFFmxP8yB8IcjBj7/J3bF/DArLhnYRjyQiYGpows3
bI6y4gDFR9JBOwb4jkUk3FRlRZnDstkMwTHgoNU6FHeaFXdWSZdnsj0EfuPj
zn1XWVFIrMnLQDAWZEy2+AzZsuKhwPTNxcqxoGo6HD3hwErfDzKUTC1iwcHS
u/yPKyveapvqLMyJhTq71+sCXqx4ed9or6oK/bx5RyHM+BErXoIKizqPOPC7
x5RimcCKVf8yMwNT4yDeq6vILpkVNb0s8541x0G9r/PH2+ms6O1+nIt9ezwI
huYJJRay4pu5gIXjTfHQkLQrt7+RFYVvHSAYyyXCtmqeRpMVVnSx95IQPfeY
vv8PDOxfZUUsEXzuRX0MR80zFpjXWdHR8nj1iZePwYRRYUc+Exmpt+Or+UVS
IFnj7PU1TjKe9ihKaV9IgR3tkWKPxckYZq08U/I5FWRTLxx2lSRjAfu9y5d/
pYKKi+gZDWky2shM/dDd9ARMtxTfnt1Pxr5XkovfTj6BlPMdnSpHyeibscV+
pO4JiHzn9P6mT0aNqD7t1eo02F/8+WGJARlNbTf9zfmWBujz+FmAERknPL6I
JdPvN+a7ZHv2mZJR7A8onFVLhyc2OrLetmSUWPRN4mxNh53zIQMi3mRkIGuP
ak1lgBgj+bBlPhnVngYcunQli75eIjJuRWT0CMxltqXf3xTKlIT9Ssj4keHK
6a+1WaBhd5X49CUZexbenyKLZINN19vOiQYyEgxCJO4NZUNBmqedy1cynlLP
anx3KxfK/z285DNERk22vZY/C3KhxjBfN3qEjHukPL29R3Khk21A8cUkGdWZ
DtAk9Z7Bouthpj+/yTgqNNtpvz8PDqstJnpzseFV4XmFDywFoJrMHh65mQ21
Rx2VJNUL4OTyTp80Xja0dqrt3ORbACbP9GzqhdhQSJtTO5pQCJ5bChVYJdnw
2fx2DVWGIno/27aFIxs6zJywLuQuhjN+39Yfu7BhytXayXWXEtDPpy3dcWdD
99WDMufiSsDw88a0kQcbDtnPb12vLgETSZH+zffY0G1IYstZjlKwazUvpway
YahSbDd/QSn48Qw7OaawYeLeiUcNhDIIVF64qp3Ghv7FVxwbZcog5CqjqWQm
Gw72ylySvlAGURU7Tw3nsqFiSIH09pIySDGxkjQsY0ORRbncOvtyqMwY+a7c
xoZ6/prCs3MVUN222CPYwYYHusXk/cReQu1v5vblLja06i7ZFmr4Epo0xauK
e9kwfsvuV0tVL+HT9JV48RE27JBUeX86qhLmFMb0OFbZ0CaWoY1P/zVINEzU
DUiy4/hm54CTQjVwU7fy4jNpdnRWcGXdf7QGWvqDF2/JsmP4v6RQs0s14Lwg
LcGvwI7+h0/apqTWwCthxyA9YEcF/bbLZVK1oHt74cw7A3q+JjauMZ06SCe+
HY8xYseeXZJm9W518Dsi+r7VRXbsl3esWI6vg8QshRKCBTs2qhe3WY3Vwc+P
N/iOOLDje/H8GRvqW7izb3Wg0IcdVaYGMiY+1cOHypbrnv7suHufmfMDhgYQ
0Uji1Alkx+Rrelce7W2AdybKOBnGjhKd8LSP2gCbgr3Sd8Wz48A9Zc4jKu/g
6QiTXXwhOw422Dga1jZCRzxl2WeAHZ8J9ZflirSAUIHMtQff2LGJ7ZRIv1kL
XK079SN+mK7nkvq2vJQWWJsMbS8aZ0dukExXFG0FqaMc2V8X2VH7FXF0h0wb
+H3luPAfBwWl5ciTHy3boX1hb7v6JgoO2nWf+pnXDoIsp9X1tlDwzDbT1RfL
7VC874GcjQAF3TzEFVtPdMBXb062WHEKrtYoiOqPdMAhMa5X88oUnL6PNd/3
d4Gv0j65NaQgVV6E75JJF3zQPpPFepyCLCFHLE4GdsEV9/Bo4ZMU3LU9oebu
UBdENXA56BhSMErpfrdVfDdMW2/anuNIQb1VwUMWIp8gNY/byzyFgvvYhwZ3
POyFn+FSf5TSKKhVJh14sa0XZNzU3LgyKWg9w+cjRuqDisPXratzKfjL/sKe
i7f6oONd75ltZRTMdZcOirr6BTa+p4j0tlLwgXeq761LA3C8viKhsJ2C6zcf
94mkDEBwVidPQBcFOWJ9PB2GBoDPkYms2EvPt8gTvWw9CLKrV+ejf1CwRj6L
wnnvK5jyytad+UvBQyVyDPaXv0PGisbR3WsUVPjVbdwW/B0m+83K1jcouNbm
UhNR/B1upEXm5hM58MX29kYVwhCE7VuKYufmQCWVeB6DtCF4rfXaqnE3B6a+
Pixg9ncYhDxPEsGAAyVSNJ8dYvwJxpPPTc4acaDm1JdmqyM/Ic5QqNT6Igc2
3vrYLXv9J/Dum7aMsODA1qGiHXmTP4Hre1jN0DUO3FhfY7s1OAZMqp13/Kgc
6LOTRW/byASoFhzujgviwLtMra4HJCfBWzBNKj+UA49abhcjX5uEtQXXvk9R
HBhWv/MYw59JWE7nUdqdwoGGHnzUS4LTMEU0Wmgu58AnjCUPxvxmQcq1Rutr
JQeeOSkhK/F5Fmy+Sj5ZqOLAlqwpfh1JGoyW/TkjVM+BcznHekzbacDc2SrW
08CBhP8fc/B/0Jq0fQ==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0, 4.5}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellChangeTimes->{{3.4166021427218*^9, 3.4166021516605997`*^9}, {
  3.4166022118297997`*^9, 3.4166022316106*^9}}]
}, Open  ]],

Cell["Now we can introduce many-body operators easily", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"solution", "=", 
   RowBox[{"NDSolve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"-", 
           RowBox[{
            SuperscriptBox["u", "\[DoublePrime]",
             MultilineFunction->None], "[", "x", "]"}]}], "/", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"x", "^", "2"}], "/", "2"}], "+", 
            RowBox[{"0.1", " ", 
             RowBox[{"x", "^", "4"}]}]}], ")"}], " ", 
          RowBox[{"u", "[", "x", "]"}]}]}], "\[Equal]", 
        RowBox[{
         RowBox[{"vals", "[", 
          RowBox[{"[", "2", "]"}], "]"}], 
         RowBox[{"u", "[", "x", "]"}]}]}], ",", 
       RowBox[{
        RowBox[{
         SuperscriptBox["u", "\[Prime]"], "[", "0", "]"}], "\[Equal]", "1"}], 
       ",", 
       RowBox[{
        RowBox[{"u", "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], ",", "u", 
     ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "4.5"}], "}"}]}], "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"u", "[", "x", "]"}], "/.", "\[InvisibleSpace]", "solution"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "4.5`"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0", ",", "1"}], "}"}]}]}], "]"}]], "Input"],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwll3c8Vu/7wO3N8zxWZURCg+wVOtdtplRCicyKUtm7FEJGZY+sZBXZRTIK
9YkWWdnzkfXYM8r6ne/r98+5Xu/Xfc59jfsa9zlw1dnIjoaKiioRf/xPrrqK
zW689MEUbv+/tH0saMX/zATKshZi9+ByX5kEcTrtOvhv2uqScbnaKHezLc0L
jsiahZfjkny+gqchLQTumd6Ay7js4dmYqU1LBN/8U7w8qYlgvLG3uybtJTyB
6TWuxJcQc/PVmeq0CrjQP5ovGlgBQcEJDVLRX+Bll6N24d4v4C3ELMr0oA1O
2TLzHaprhdP8bNfVBbuhNafG94RdJ+w4RR+IV+wGu4+d5zV8OuHNJ+7BmbPd
sKCmEaf6uBP4bgleSPXrhooYE4z9dSdMV0pp/hvuBr8/SlpyW50QbmK4vzKr
B65cim5xj+6CrzGJXfKH+2C10KonrbQbdJkOnJSQH4Qs7wJjcnMvtHTcFvl9
ahCyQ22ZFgd6wfR5xXaKzSBIsRdMr8z0wm2lM+UskYMgqHLlYD9zH0Tbeh+Y
nhqEkELvbnGdPuivb97MSx+CxoPWS8VVfVCpfsKUi34EWlgdE2xS+uFe1qiO
M3EEWr+LxrLm9YMGU5j8D/4R4NBvFit42w/fO9o5guRGgLNug+praz8M3rJv
WLYeAZLfQ5IPwwDQJMfJtVePgD4lnbnYeQDOrlHYYlzIMN6/oGtyfBBW/Zqd
R+6RQW90JTtGdxBSmF63S4eTIa5ljKXOeBCm+H2SfmaSwWUwe3zBcRCCNRlE
2TrI8KXf5fKrzEF4HyWiFqo4ClYlnz56MwyB5FHzm/f/jYJQNl++fcMQdJRh
TU0Mv+GlhkzXYOsQ3D0hIi3A9RtauOOStAeG4Ot5ymq1xG+4Oxn0a2J5CGy9
vQM2LH7Dfw/keu8LDUPa57in7nW/gbd2t3LFcxjYrJsabgaPgbXHP6Z11hFY
dH1d4h89BvLkj8Jv8Dj9Ck5MTkgbg4bX70Wu84zAs3wbx4/lY6DHrDJauX8E
pP6scu8bH4Nxq0ZVZtkROB8haPtVZxxChZRrL1wYAYUMmnNDhuNgKjkonm46
AnvLJpVXLcehXCXJecRiBMg9r1mFvcZht8T5uYHdCLiJ6pZ5vxwHO9Ychx7P
EYircaI5xDgB0zHUjnEJI+Dz03hWnWsCIkpoX91KHgELskqXkdAESOp11qk+
GwFRRtp8P+UJyDfH7jfkjEC50VPDrhsTEGeeLP+wbAS6KHUZD79OgJLYhwdF
LSPAt5cEY+GTcNVB+ngWDRko13tVjiZMwjGTDLOr9GSofJsp55IxCdYxLToC
TGQwMZYT36qYhPSspsD77GSIiTRm5x6bhI3gj8lre8jARJ/Yr4VNwaCBcup7
CTJ0G1t1hp+agvcZ0Xn7pMjwMku8peXCFFy63RbmIkMGbfTuk/ntKVDbY/CZ
VZEMAb49r9yTpuDg6wnljRNkWF/i88lenoKN0jNZr86RoRH9dqVsT4FNw/0j
OefJEB9VcFuamQL93CgyyYgMcsfUrWuEKCATkLriZEIGJ3tL3Y4zFJArvqrd
bkmGE+/E0D5TCgS9r/6Zak0GNob549bXKOC+oXnC+goZ8rP9js3cocADksf3
DlsyTAw956bNpUD8b4Y699tkeHvMnuPUGwp4lHsFszmSIfieDFPUBwo4V0Qq
Pncigwjfx02+XxRY0iqzeuNKBsuLo6Ny1NMANY6UKG8y/Poh+vqq+TQgjtjW
w4FkeIhJCKfdmAZJGpO+u0FkUHotG9XpPg3rB2raG4LJkJyIOeo9mQa23ogE
3VAyWF0xPSL1YRqqtI2jVh6TgfDLKunGt2l4s3f3OHsEGep17RgzO6eB3pdl
RCiSDAcl3ca55qeBqnFFVCaaDFN/Hmf+3T8DO5Tlj1tx+P43Y4nyEjPQG2Ll
PRRPhtMDSf4OyjPwiIlOqSqBDEUfX1gOG8xArcKJLrOnuH6Fwqa9FjOgv+Nb
I5REBmLuGzUj+xmQmzj9eghntyd1+xoCZuDj7v6ukym4/t2GsJ0nM9CyssS0
jPMv16Z15eQZ8ApZNE5IJYPypd7O/Ne4PlM11S9puH3fhrXHPsxASbl378Vn
ZEhRnygT/D4D58bqkgZx3jqwEhM9OgOHhgP9e9Jx++L/Un2fnwGbHkKB/nPc
PkYqF9rNGQjYatp5hzPxLsOwOuMs7FlrvieYQYaPs2znvLhm4evTo4fu4exm
zfWhRGgWyOSd3Q6cRdv3SVIkZkF56xazGN6nOrWFU0VUZiFYJU7HGeeQd+Is
Ftqz8Jgjo+INzspHj91JOD8LOvtqrRZwnkqTn/ppMQtU3PuRWBZuP0H1EtPN
WQjpXrW4gLN+IGrU8JyFk5d8393DeWtVV9H3wSz4ps8bpuNcfONsTnnELJxq
jpKqwtm6z5hrPnkWZnlyDJpxJp69HHjo5Sz0VXm868P5Y53Nks2bWZhhMLAn
4+wmd8MmpXYW7l8Itv0fi75wbOn4PgtPqW4V/+/9zj0eGHv3LOx1NNP7334h
j+4W6f6ehbatEoX/6VPZDhAIWJiFd8QWz//ZQ3EOfVy1OQvbdAIc9//nz2jE
v2XGOdCWoqK6+D9/LsbflOSeg6quqbPi//PnS0qPnfAcaKmd/LeIx6NINfPk
c8k5yAiLZnyLs1VRbkWPyhxE95zwccOZIFwsxqkzB9tSTeeO4FwfWx6vbzgH
fN7vI3rx8zjo89Gt9uYcpG31nhPHuVa8YOmQ5xzs9ZZ9+Q0/38ud8S4xAXMw
/9Lt5A2cY+XsnWyfzoFJ9l6XGDw/pMiGc81Zc0D7592mIM7fo9QclIvnQIGh
ojcHzyeaOY5bLA1zwKZ3LSETz7/01I0p99Y50BG09d6Ds+rp0RuD/XPgNnu9
PBTPV7fct3aly3NwjlVR7xKez6NXLK6YCM/Datxv7jK8PvyIuiN1EvMwF9rK
tpVIBr46aesjyvOg2xbhCjgbCdBabp2dh28pCY/e4vX2seuVaZbvPEhXHVez
jiGDxcO4btbQeeCOWev1xet1Q/6+iWfsPMTeo96IjSKDTMz5C3qv5qGJ635N
KV7vGfrrBvNd+Pf0n7Ly8Tmr/m/k56XRebh0Mt8vNYwMPXnfz36cmwcfLsJc
CN4/iAzp+vF0CxD0lffx2Yd4v63XPqkqtwAqopKklAA83xRjsZCIBSh3P+8l
50WGf799PywkLUD39b6cZg8yJMbaqZvlLICa3LHIa+5k+LmgoipZswDTHkmp
/vg9AMsfVmyfWoBXD7amfW6RISoqkZcgsghXr6SWuOH9957Ek325hxbBuUT3
3n/mZLD/8kAAji1CfbTJFcJlMmjuOog4qywCfwutfxLev1edtKRazi2CY1ye
oYsBGUzPLepE+S7CCiUqWwXh82N6XO/Qg0XwvcMhp46RQTakX78uZBGKOK/Q
qqmTgaW20XAhZhG0Oe3uHFUhw4djzywN8hZhcyhst1Ua79fspz0JnYuQdTtu
GNtPBo5X4JPbtwjR6r/H1/nJsKmt6AsjizCh7GJeuA/vX/eEHzjPLMKF/c+e
M/Lg9TP750kL9RL8kDs0EMpKhpkf2TlRUkswsD6V/3odn7+PtjsIYUtgSdsm
VNc0ArEqsTEHIpegQpB+p+jbCLhMiBvIxy+B13/XZ542joCE5vnvJhlLMP5K
k9ayfgSy/mbVP6tcggjR2vhCfH5H2+sVSVCWYG3yVKktPv8ddONCTp5eBiql
hMoLV0bgIM2R4/5sKxAbG7VX7usw/DlbsnOYcwWevT9ser1+GL4nK35u37MC
wqzVDXGVw+Aiq2kgfnAFeu9H6PblDcN7G3O75uMrEKP6NZM3fBgu1EXE8N9Y
gbymnl4jvWEIubdMefdxBfweRzmzfBrC50lNyqLXKpTlyofT5AzC+g2unNbd
NTjqEy01+6cPiqnSIh48WIePtqbCFwR6oIb3ibsCZQM0PUNvrJt3AucrRG3i
/A/4swpa0LV2qKu0eHynfxNoXhzRVHJsgXfaJz6wXt8GMLs59Y/qO5hsnGtx
ouyAWqPatKBkA7zgHpu/eYQK1Yg89Xkg8RF0Zmi+G89RIYFkJuGIjmpIUbT3
lF6kQi6T1z40V1XDon+zMOsKFerEGvWEMvB17iTvTxtUiEVu58CaQzXMnzgm
JktHjWqMz7+8zVgNSdGXAjj4qZGkvi5nIFTBtGKh0jc9apTNZ7hH6/07iA64
mI3lUKMxzv9G7CfLQeGY7yJrLjUaH3lSotZeDr29GSd6X1EjnUZBqv0fykFE
frbbvYQasVqx8JLiyqF8PJDtVTU1muRNeF4K5dBzutSTq40abdN8tm1LLQMh
Hla9qW1q9LzJk/Ou9RsoflU/G3OJBpVH/Uy5wFUK+sOKe9Iu06A/HTHR9lSl
MMVdoPHSkgZln/ZO5+svAZGAhMTqazRI2GqKizO6BBJNbqHfzjToyJm0MzGb
xeBPyx2vEEaD+k0fDf3XUwTnLW+odVfSIN/JbMXIggKYix2wI9fQoBxONaof
4QXw6Kth9EwtDUrTveN8yL4APsurj+9+pkGEhyc3MbECUGEhRR1uo0FB/kHU
kxn5cOBd9egdCg0K5IxyZ0t/BctEjkcCfLRIs+aoq0pJLjziuVe5LECLxpfK
v+vF5eL3p+mJr0K0qLA4dNTNOxcMDzZqeYrRIhK3EA83yoVSRb/tZhlapO9Q
vB3W/hKcLs+7+J+kRbrRdOEyWy9gOhtvfx60qEzlET/X1RwIzFMPqfSmRXqb
j7vgVA7wFeWXR96lRZ6/48KCZXLgVEUISS2AFvmLjz29uJsNeV+xH3GPaVFV
isRP7vRsuD5XjOlk0aIDHM0T/4azgKwUJZbXQouczW517XpkgtvnM7I17bTI
yKvg0XHLTKAxYjnxs5MW2T20QhE6mSDi+PDCaj8t+hZ9zy6QNxNss32D0BQt
onPgeJJZlQETBHtyDxUdMnr5ZS2LLgNmptAzZjk6pJBse8L68jO457WTJ6BI
h9psb1SwyD4DNrr35dIqdCimlxQ+zvgMJIWUmy5idOj4yy8vD1ekgdPFY5uZ
p+lQV35Q+Qh3Gix+3Gemeo0Onfb6XrzelwJrKcvct+Pp0Ju/AQu1YUnw65HL
J+endKh9xc2l3C4Jyu/OO3uk0CGulRSxDs0k8Lg8/f1+Bh26xL+S9Xj7KSzv
Gw2ILsT1hxseHvd4CgtJbbNvP9MhY/o5R4HbiUBJLP1MtUaHzqVynyr2iIev
IVJu9Bt0KHJA2DDMJB7yvAqFWDbpUMpJWqNQlXiwN8m7y0VNj0JKYqpWtuJg
gidDRpydHpEb9R8EhsbB7/ioNH0xesR6slDMPDsWBmOdPBIv0qNPEbGj3GvR
UGk4LnbDlB45vK0KYh6MhliiRbeyOT0S25f3ja8hGnQjT6n22tCjPCWzwdiE
aCgNF6UScKBHb7gHTbuUoyEwoO9xZiCu//GLzsO0USDupJtdWEKPjs+0Kvsa
PQGqYx8u3H9Dj1IUbKkljzyBvhl5hnNv6ZFOwpFuAtUTiLp54OZCNT3iyKsZ
dy5+DBu2W1JyjfTo4XOHU0/YH8N38zfV7wbo0dRwmgFbZzg4n9rf/omZAfGV
NmcOh4aCZsyYTSsbA8oXU5R+diUUuHvzFwYJDKhct68/Si0Uqu2V2f/yMKA3
i+qFXIshwBB2Xk9ahAEN6u+MRpqHwLPGwNpUVQZk6Gw0f0D9IfzQmsz3uM2A
aoKaDRREgiD9cdHxQCcGZL4T+P0wTRC4dLh/jXJlQMaPFLz0RwOB9xrVRL43
A+KVSnzClx0IVwP3HSAHMSBv5c2lXvFA+Fuv//RsGgOaOztBU6r0AA5hpUHi
zQzI6QxHbN6sH6xecOalamVAKk4ZT4er/eDTbalXve0M6PZpzr0W4X5glVz4
80kPA3qV+M9A9JAfxK/m7Vv5zYBI7Y94KHb3gaogs6TuHwPaeXg2lGrRF3r2
xPWbHmZEE9ckCy4p3IGXUkZOchKMKKHkQFA+yx1w1yFRs0kxok+teXeOk32A
wz1KvE6eERV6yL9kifIB7Z+PXcWAEU11Omm7zXnD6+BgxuWLjGiJjxy5760X
PFrylH8cxIgKgllEZx09IPfW2g2TEEZUe6ilaUbfAz7/9kg7EM6I5pTf7Bc8
6gE7ne70lZGM6CqndjPXpDu4Vrt2jyUzotsPBJtPX3MH0yBHX6yUETnGhNUV
2LqBGI/df0sDjOhEjqR1brgLaEaOrb8fZkTv8rkDbju64P8ntpJho4zouTP7
kcuGLpD092rC/ilGJBv9QL1inwuwDdnc0F9hRGbNgoe/mDrD8ksL1hfMTOjG
Qw32728coV7lgrGpEhMKnGx3Vs28BV2qB+1sjjMhbd0TVblOt2BGfdnLXp0J
/ao6IXdL/RbwakSn+GgyocEi2+bZnpvgcKqJnHSOCe0mzt/v57kJey5ru/Re
Z0IG9/eKrKffACdfhYjLT5nQO92DP+hmbCH4Pu3zqylMiG7yo73SJ1tI9m8v
vfWMCbVevLtbnGwLn4Ncft3NZkL3OKTc3E/ZAt+TQv7UEibk5Ty/c6LwGjSk
ieb3f2FCfWwmQ/5+V0GglvuLxQYTSk38tpfK2AYkowYGZDZx+yu6T91TsAF1
m5xluh0mVJx/PEGC1wYsaBT2F9Eyo/7A7Tztcmt4pmvsuc3BjHIChNZav1rB
/paYg+lizEh44QmTLJ0lCI9wBAwbMaNK1eQazXozkHndlVB2kRnVUXeeXow0
AxSYXhBqyozO7MT0LVuagY2oVLe0FTM6MMdItbZlCpn2Z6UCbjIjSaqt58OY
KYgsPR4QDmBGLo/OtJLbTPD7FfPxq0XMaM+HPbdNZC+ApIWwpHspM2InHA2/
wXABFCqUhYLLmJFOS9uBH/3GoHvrOv3LKmZE2mx/GhpqDPbt/7VRGphRqFZQ
WuOoERRn+d1yHWJGtxx4iidzDeHdVoJlIJkZma2+5SQ/MIR6k6LzcWPMyH/v
zeeiFobQxjKg9HaaGTE/KVAKIxnCqttx2r9/mJFFl3Jumf95OK61mhpAYEGL
cHvWuv0caD5jjYrhZEFON9TeHow6B6fXRQKzeFjQ8+Qf8Q5nzoFFgaH9Z34W
VK+t+if861nw4ypRYDrMghQah8T6Gs/g9XCzOQqxIO5aFzG2/tNgEDy8k+7K
gho4a7O6zE+CUdHCmq8HCzr/IWGcQeIkmHTtzpp6s6DIzBKe+H+6YHFYuJ/z
Pgs6dUXsx6cUXbjVZPMuJIwFlSR6RtoM60Aw96iz03N8/Wti5d472hB2Yvm6
fhYLansoeqD1rDY8vk5jdfgFC5KM2uLoF9GG2EqRM6P5LKjOukOd0qwFzy2u
HTapYEHykTK3uw9pQXXO2MiJZhYUqMy/N3BaA2qbV7v5WlnQvx0BF5nPGvDx
D13LejsLshZKrjySrgFf9cQ+vO5hQS9kEjk3jTSgc9YuWWyMBb17adHDUo9g
UWHSkH2TBb254aY0cg5AvIHyaeAwKxL2vtPFfF4NfM5XmxdIsKLkxNoOZx41
+NH/aPWOFCvKMH6ocqpfFVyWJcT3KLCiZqhgM7ZXhRohp3BDYEWJhWZdhaHH
4fzdZYPGi6yoVUn6TVOfMmTT/zcVb8qKHDymMbZcZfgTHffgmjkrEotLYy5z
V4bUXIUyqiusyPauv6c8hzJM/PLiVXNgRefLGJNnTiqBr/TmQEkgK+roH2N5
1qwAP6t/ePo9ZEVRrZaCW88VQFg3jeNsGCt6vAS279wUoNHiBJqOYEUrTlfn
YZ8CEB/5Z4sms6Keo0lfQm7Jw8sx2lvJJayIx4A60kpQDlqT2dYDB1jRxPfg
MaoAaeAvlnSMHGZFr5zVrnrKS8P1T2d+J4+yItGz7DsJk1KwPf2kpXSKFWXN
yHJ1GUrBUXX2vKFVVvRU2/WVhNQxCB5iN1NlZ0OdfAIPGekloGX5WIsOkQ1d
v/trIv3TUeBjPKdjyMWGFg05ypoDjsJr6UhZ+31sSC7ub/7szhEYCuBgSRRj
Q0H/qShw0BwBlYOEmqUTbOh79/OQXpFDEKQsLbuN2NCCTUWZGUUcfuob5DJp
syGh4LEOq1JxsPOIihM6zYZ6DO8lr4A4xDYQHM6asKHX038Fj9iJwewNouAr
JzZk/fMwb1/jQcgoJPnbPGdD9VjQK+8RYZiIOvpXOYsNrdFn0N2pFAZJdy13
wgs2xFtVLtcRLQyVxz1v1OazIeXlWJ0oTWFobewxEKhgQ2+Q5AnxC0KwO/Jc
uKeJDUUM99XakgRB+3NlSkkLvn+ejvfwDwF4lNvGHdrOhvT4Wl/+CREAXida
ZqUeNtSqZR65tsMPUpvXl+J+s6FmV9kNrj98YMUj9cngHxt6sXw+7cfuXsjZ
0FU/tM2GCLndBv7/7YXpfuuKnV02pHSFpeh96F7wyorJL6JnR5UWDcstpL0Q
Ib0Wy0piRyEr7WdUJffA+1Pvr305xI62NplOKtzhAX6/0/RwkR0lrhN9Gh04
4fL0GwtjU3ZUe9DaeEWNE5JM+MtvmLOjIqvbP9NYOYFHevZq9BV2pGZeffZK
IQkIIxH1ZEd2pP+t2S92hQi0mm2+wSHsyPKXE+0FaQJoFh/vSApnR0NXxLO1
BjgggC/raNETdhTno6TJ+ogDtpfdejtjcf2/fIf5pthhPZtb+dBzdrT5xsqR
sZgNZuhNl7+9Y0eu9eKvHlizwFG3+lND1ezIwJUpmWcfC9gPHc5c/sCOmFO2
a207mGG84q8B/2d2dPzY54Fjp5lh+EZqkUMrO2pMv8J/AjFBx7chewKFHUW0
JjIsnGGAmgi74Qt8HIhI/ejmySEa0LOI7LksyIG6bkgckS2iga6j79pshDmQ
nqSwZdY9Glj+wvTZQZwD/Ta+Ec0rSANHaQrzguU40OTFD6XittSQ5LXsWn6a
A8nXXjafZ6YCcR3+29VnOdDH14mKcR93sXIubdv68xxI9b2KlaLLLtZSmmDS
ZMKBnAsqwhRadzCGGRW1sascyE+jbWkiZRvztPan4/blQFG9MOxltonRSOVt
7/PjQOwdfp+peTax6K3WP0IPONDfFt7Uy23/sMJkEYpEKAfyjzI/6XzmHzbW
0dCsFceBHtmuK63p/8WM9Nieuhfg9tduG9A5r2MjvIrRd4o50H0TI608pXXM
adwy3P81B9JwFc2j2v2DhQeW+D5+x4HKuZOEv8X8werfG9lk/8eBtHSrVGPr
1jBp2eQjHf0cyCSPIXdAdRWr/Fkt/mKIA0XmnNORJ65iyGHgoDeZA9mHoQ/c
EyvY+Vyh/fyTHOjkuSBzj7gVzE3wJee1ZQ6kTG8u8GR1GatgLttcZiKg5yUW
jKYtSxiW+2vjMysBaeVIzy4ULmFftP+sJXIQEAPD1Ta6x0tYj//xRVVuAqo7
qxGdpbeE/VurGwsUIiDG2CeMyd8WMWy0qZlTiYCKHkUK5Py3gDX6z38fUyGg
4FMdVo0RC5iBIPFrhRoBSdxld20yXcBsTI0/XdYgoG2BCG3/pXks8GdvRdYZ
AkIn2p9aH5vHWBw2yzwMCKi2yFwkYGcOi2MWfK1rRECH8nbDfVvnsBxtm4Lp
SwTks1g0uOM1hzVWTzyXu0ZA64bUKgI/Z7FzpkzP6K4TUMczzUMbL2ax7rUj
KV32BOQpvPdirt8sNiXjGH/XiYCOhGjt+MrNYiy5K2Gf7hBQ6q2Q3MnsGSxW
myck/h4B1T8zeaYbOIPxjyoFXfcnoHmaQwluV2YwScG791keEpBskF6/lMgM
di6Oys0omoA+3HIfEi+axlqqXM9ciSMgG/tLrElR09j5kVFxl0QCsq/6aLfq
No0ZHWvoj0gjIN1lztXb6tNYh7FSRdpz3L9PHqMPRaaxC3dzowuyCMjp/nv3
MOZpzORLuM63PAISoKo8qdNHwbrn/gn1FODxmRsO3PmPgplyO/ybKCYgsX4w
Si+mYJevnCuhe4vH59XtzagQCtYfWhfOVUlApxw1OMnuFMyiWMZWpIaAYmub
lPdcpWDWm5z70EcCip69HCWrQcFGDgSvnPtMQArDFoY8chTsqt5as+UXAvqk
L6ozeJCC/Xa6nufwnYD+TLPqh/NSMNuE7kDfZgLSTjbT3cdCwcZr9CwfteL6
PtULRuxMYddHq5STOwhIpmmsZmxlCptkkuDM6yIg3qB+gvD0FGYvnTZb0UtA
zsEl7ECewigX2b80DOD+548Wa/ROYbfu+WX+GiYgW+/F5UPtU9hM1oLv71EC
knrQOrz0Ywpz+GZjsjxOQHEDPoHPG6ewuYU2GWoKARWH2SzIfJrCnHi1WImz
BHS+j1rrZe0UtqBePr5/gYBuTlpk7dZMYS7XxOqPLeO8fEftRPUUthSemKK+
RkAFnTM61lVTmFspo6f+Bp5PTka013Fe6fIxuLxJQOosg3WG+Pse25QjN3cI
6OLPpnaR91PYv5hjNO+piIj10Hp6L65PTey6sxYNEcnJ6sV54fbcq0wf+E5L
RKNiOhL/cHvf63efMqInor6GyAXbpilsa4jwrpeBiG52/vV+h/t7wk1P9AoT
Ee1asnks4/Hwo38QM8VMRKfSRFJ5R6ew2qSqHWdWIuJ+YZAjOjOF7Ugs315n
I6IHUconBdemMKg72uvHQURGHBvq1FQUrH48tSyCRERjimc5I/dSMKo7v4R5
uIiI8U2+u4o4BdNgY49M4yai9VStY20KFOyTnJ99wR4imv0ZNdV6gYLRNFZ0
yu0jon/vZdKOX6dgWmYLmtV8RBSTUMUe70PBPvvbCH4TJCLjfdLLBzIoGB1X
8qPzQkRkZ6dhe+ktBdN52bbeLUxEnkTT+MAfFKyxSbN94iARaU/lVdb+o2Df
9omF0h0lomR/XfEl22kshPvyH1cJIpLcd3iVOmAa0yZE2Q1LEpFtwQOZ5dRp
rJ7ur1a1NO7PtvMRr65prHKhicpFkYjYfh82kzWZwbymqV0GlYhIV7fn21Wf
GUxhXGn4tAoRMe11DfBLncFK+zI+iKkR0ZcZl9i7YzNYboPH3X5ERBJi5Rkl
D2Yx2/pXFD1NIjqoI+NkWTCLidQMmVZo4etUofPUXbNYeqmecowuEeXdjB1T
lJrDElMFVnXPEBFNye5f3+k57KHrZ8c3l4iolKWQQzV+AdN0+DsgZEZELf0+
oSo/FzCqG1JnIi4T0Z8H9+Vz6RaxexZJR+0tiWhoNmTHVmMR8zjpMCl4jYj6
s07zZH9cxOQ0M00e2+L5cvW/dJHfi9iielfDhh0R9fqH0KnQL2EOciinw56I
bBRHyUf0l7BrgtxXwp3w82xKyjAcWsKE955q/eNMRIUOQZm5dMvYEKcf2LoS
kay93b7nEsvYZaYpQfAgoklul3c+vsuY0UpN3+odPL+ym47THVzBNL5du2Dz
kIje/7yWTx23iomv+cQ3huDfXx+r0PhvFWM7EPlLMoyITkzKuqisrmI9PpUX
/j4iogNFOcmGpmuY02H2i7HRRETnTP/75pE/WEpoxcXPqUSk9+RJJPXSOhZQ
9iPh6DMiMk8lvs84toHZDY90RqcT0d0YupC9tzcwGWVWE8tMInpxzjg5dWoD
a5ywNvnzEq8Hf4caxfm/2LIO86XDb/DzdXKpIIpvYT2u+59GlhFRyd9RGXfX
LezDM/nu1XIi+r4VQN33YQsLW7O89PEdHt/Q+M56s21s/4s3ly5/IKKHI66x
p9N3sNN0FqZPvhJRk5Drr4RiKpCRcU1a/kZE9tPpyhKDVMBjEdJj+oOIyshn
cizZqGGkrNRU7CcRJdFd0cp0oAava/RmtR1EhOnzCvcr0ED2p2KzxSEiSvxT
e3B6nBboym4zeY4QUVWoRPIxUTq4nn343V8yERGqM2OMrtHB0eAsbtpxvD6W
70acHaODNzqJLTwzRHTt8YK89RI91H+9r6O2TkRh+l+4B1WZQKRKdbVug4h2
Zh31PoYxQfCr9Sztf0RU59AX+6GbCU4+cqU6t01EPMwxhRx3mPH7vV2NDS0J
fdjjl/zwGwsM/jwjG0IgIV1XzRCvF+yA1TGPsJBIiHN/zp9fDByQUdIYGcVJ
Qn160n0xNznANhrNJvGQ0MN4n2+q8gSYNVTILeAnoYmHRcdnbhPh3y9+wbZD
JHTXJvSwjyInWDT0/Lh4hISEDNJ5zdw4ofZtwt2+oyR0aleiiKeUEx4kEnp+
HyOhrRkDdSkpLmAypY37I09CdXe/3VRW5IY9/dNMAoiEXn4nUFW58sIfuZyq
EA0SOs9mJfu9khe6HlneWtIkoTmfC59/UO2BeNXWH406JEQh9XEXxe0BYsrb
SFd9EooR9Xql27AXWMwCuL+YkNBPfuzY6mV+oJQeb5A1JaFC7vLk4jf88JVp
xTPNjIQ+R3vXxbEKQEilXberBQnlsEv8mK8XANq9+ikCV0lIeT3mr5Xqftjq
4hF2cyQhG/+4L0wVwtAv1dI64ERCZYwTg1qzwlAdEvbgpAsJabfxr6scPAA+
SpujAu4k1M8zEmgedwDWEkZefPEhoSt8Km5N90Rg/kKBhGAwCVmMhSR3eYlC
c6HtQOhDElqJ/D66XSUKhXT7I5ZDSMjz5vJl9h1RuF0ePf8lnIQylvf8WgoT
gykur9duUSRkmqSWLJUvDuR2pPw1hYTMAic84wWPgPpCpolhGn7ec4k6M95H
IImV1qvvGQkFaZ8retpxBM5pN5TPZpDQQMGh4wxRR+H921NypFwSenXXrbaN
VxL2tucbpuSRkPesga7DPUlwn2d1PZhPQo8lfuSF/ZaEo4d+ligWkVBT0jc2
64pjkJRkdOxyGQkNW0fKyzhLw2p52Znf5SR0vfzx0xCyNBi0cTs4VJCQCFdg
/tejMsDA0p3vV4XH/+R9F+dMGfDwNT+cXYfnS5lVX2SyLBhY2R6c/UFC12g5
200r5SH/boOmZzMJmT+YDhmgUgCGp+JXd36SUPtOUYL3aQV43zKVQWwnISua
E8r0ZAU4quG4X7GHhBpCzm9s8StBiOXPE7W9ePw8z3pXOSoB+Y605cl+Enpa
lH36Tb0SJL1ZSjUbIiGfFo81h1vKwCDqtc9vjISoq69/TWxWgauoW4VxgoSo
/rqtNBw5Dh8sVEyjJ0lIQ080wzz0OHgk/EvMmiYh4W/3NNp0VWGUwZ/7yyLO
zM6kR91q8IESSiBuk5Cj8I2Hr3QA9h+UbeHHWUU3NFLsFICfRV/kIZxXmf12
M88AYC2S+E8eCSlmHtiXZgxQW97G5oSz/0iGZMIVgDp/fpYmnG99CBpOuQ/w
kaeELmwHj1ftxebVtwAiBqaf43Cuiqyh3K8CCAyjDn6OsxPliTPjBwDNTWPa
dzinG754x/8Z4NPIBvUEzoW/ONVOtQP8V6C5q7VLQlJNPHyV8wCi4zO1Bjgn
yyt0nFwGCN6f4GeOc7HG4z/dawDaMZPbbv9b7+W993cL4LPXk61MnMmhDJ8y
6BGIliq+L8I5fEBS7xYjgmDKkG8VznEHs2kVmBFoW8hutuI8djOZ7ysbgpyE
vuoBnIkW4l6xHAjoWoLuTuF8XVyQER+jYMt0THUV59fawc2HOBE0aHT93cW5
Q/7qt2UuBP8HrXQn9w==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0, 4.5}, {0, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellChangeTimes->{3.4166022850094*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"solution", "=", 
   RowBox[{"NDSolve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"-", 
           RowBox[{
            SuperscriptBox["u", "\[DoublePrime]",
             MultilineFunction->None], "[", "x", "]"}]}], "/", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"x", "^", "2"}], "/", "2"}], "+", 
            RowBox[{"0.1", " ", 
             RowBox[{"x", "^", "4"}]}]}], ")"}], " ", 
          RowBox[{"u", "[", "x", "]"}]}]}], "\[Equal]", 
        RowBox[{
         RowBox[{"vals", "[", 
          RowBox[{"[", "3", "]"}], "]"}], 
         RowBox[{"u", "[", "x", "]"}]}]}], ",", 
       RowBox[{
        RowBox[{
         SuperscriptBox["u", "\[Prime]"], "[", "0", "]"}], "\[Equal]", "0"}], 
       ",", 
       RowBox[{
        RowBox[{"u", "[", "0", "]"}], "\[Equal]", "1"}]}], "}"}], ",", "u", 
     ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "4.5"}], "}"}]}], "]"}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"u", "[", "x", "]"}], "/.", "\[InvisibleSpace]", "solution"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "4.5`"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1.2`"}], ",", "1"}], "}"}]}]}], "]"}]], "Input"],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV13k41F8XAHAx9sFMaLFnLWubSHJuloiSIiRbIUtCSH6KZKtEiCRblkIh
QirbRZZkixQJkZQtM4Oxhnfe+WPm+TzP3Pnee+6559zZccHztBMzExOTI+Pt
/5/zV2Sml3L8tQgb/39RwPGeqK1wmjkkSvBcnWN4e6kCaTL1IhwxbO9tZxiZ
PFXPpnvBU6F+2RiG55v2unal+sG3p08c1BjOO1z3eEPnBvSqf4xuWqeAcr/c
WCI9GEJXyp9pMnxH0HNNqSUU9jj15j9eo8CISblgY2oE3Ml5lTH4jwIaUWtK
1l53YSNuOIyD4YRm3aNzOlGgkMR/VmiVAjPMUbaRW2NgfwWvGP8KBaSfuFbG
0mNBJPtFz+wSBbo0o/Wt7R4AV35P4LtFCgT1F3+Wa4kHjc/3hZwXKNAnuDRV
k5oIgTWS+67NUSC8VNg/ki0JWAudX32jUWDvKSCYez2GLC7+nVJUCkRFRYj8
1UmFgKkzhEvTFDgo/yLvbWEafERr7q6TFPjd3L4/bOsTqEiV7zEep0DGrcic
aHoGsFlse9YxSgE/k1R3nxeZ0NEbtdl1hAInxF/uPWuXBdEVlPDJIQqsVHXV
SLdkA697W0TWNwqYLm3rrUzNgcF33ip/2ikg3yyflnkqF76nWnFVfKQAU6Km
w222PJDQGFzzaaZAwT57ymmv58C+JepgWi0F2Dxz2Sd1CsAzh2DW+IoCA4ff
tXcuFcBZ/hydrkIKlBJb418XFoLu3gKHhucUsH8xI35raxFs32yibZ9JgcFZ
IZEoejHQX+sOeMRQwKBfamP95Cs4Et9WXBZJgbI6xZ9XXrwC1xPJ3SPhFIiM
1cqzsCsBsb8rPLPXKaCmfGG/VEspHM6Z7J28SIE41+fHK1LLIXWvScoRdQrw
emV/+7BeDk0a1NN+exjx9Uu72Gv/BiJYC8/el2fkR2hc8Lz0W/C2MPB0EaFA
cNp/ZcqF76ByNu5G4/oMeHQdE8murgL9h1Xhu2pmYLpX5/kr8WoIf1N2L/H1
DLgNHT5Qe6savr2ztfxTMAMXp/aYDOrVQH+LtfP+lBmwYRUK29qBYSCA7USg
3wzkiNDpep51MKX2cHu33AyMiXtrnsioA3mn8pKbojMgLUUJMeuqA+S37xmZ
fwaydk3wOeyrh9Hy+zU/1/5C2oHBXcGL9XCnNuinRvdfSDBptKm42QBmlONH
2q7+hdCwh43Ksc1QEbbF+H3eNLhfe2K0UN0MT7RTsldSpsHc7XlX9VQznJ2r
kRCOmQZ5k+rB4/of4Cgzb/XWq9PwSXiM7rb+AUp8dxUzoWkQK90vm+v+Ear1
601Suqbg7cjnCPFj7eCF1Ir1fk1CVs8g8Y9fO/A8i9d62zMJUc1/Hrx82g6F
w3Vego2TYFe4mn6YqQMU5ZvmY55NAtt/Um/OvemAtKDr0xFOkyB/591c6+ZP
MLZr2vLOyARcE+eU5rjVBbTKj+EaH8fhogKZ7v+oC8K6xLSuvR0Hc7XtTROF
XeCVsUUsM2cc9p/c5dLa3wWKbbXRVSHjMBt0rOD+/m5Y+fazil9jHNyH7u4X
HO+Go8lWZJNnf8A+jUtP8mQP1Ctbi7t6/AZDYeJFTdFeyG6lhN7rGIV1j9gd
Caq9QE7vsS4oH4WSeoHBqRO9sH9a91Vl+igIuYmapQT1QsRVkMnxGIXJt8ra
Kz964UvPWG4BzyjcNT8l9jarD6ZuDpbHG/yED3GJX/ft7Ae7rS+rXUuG4SjH
Dn2FfYOQfV0sory8H07QJ4hxXiOwjST+PKPoIxDt2hpdw36B1iO9LHWzNCy0
jQy/7v6Bnn1sHziffcI9rdKvLpybBL1fmzp63g5iKf867xrXv3Coy9FA5tAY
rpHNp8kx8qChQMZY7tQYtvqS4BUX/BcCHOnrQs5j+MFeFw/HR3/BKcZmd++D
Mcz8l9eNq/EvjP/cciNzYgz/PG993lxiBtoXjmn4JPzGGUaLJ2e+zkAUkae+
YugP1lwZ7rD4OQMjt8I9v8/+wX15H0/U/Z2BiPmlc9Ns45jElm6UQKCAQsbT
D0PK4zi4VldfYy8FbHLv9RCCxrGd6gOtiGgKrN9V5DqxfQLHxCRu4ZOkgvam
D5xROpP4hkLU9lw5KmwroawePTOJXZpviYASFSxW+AppFyex9oa7pKc6FfKj
QpW2RE7ieQ8d5U5jKtjy1XxO7pzElsZUvZjrVIhR/GAcZTGFJXkMr/J9ocLY
GIgMW01j3ufgn9tPBS7zFK1wt2m8qqt6HYapUGqy7dj2gGncc0PilucUFZ78
GNjgeTyNI6YXojo30eCr//TZXV+n8VRr9tMYZRqw3XguIXXiLy6LXPvMd4cG
D9NPWW5SnsEP1B/E7bhPAyvb4c3pmjPY67fsyX0JNOhO8FuTMZrBCtomH80z
aJAZYGay6DyDs5azatPe0uB574bXk4wZHOtiUKgwQQOOFBlOZl4Kdj8aH6Fv
OAuZv3fdXuqjYMN5Od2zJrNAKFjbNP+LgndmVW26ZD4LbgY+uv1UCh5dG7tx
/8IsHEoLj9pHoGKr1+o+XwJmweHIowkkTMUG0kN2DvmzUP76/l4bfSqWYt51
8CZxDv7JG2Y9fkTFCyeK1ndunoOakK2OcylU/PGxakP31jn4okQd0cmgYq89
2idlpeZADlT5W/KouMr+nFP7wTkwD/hj4fiOis1wdJyw8xzI7iz/OdRHxXLc
AhaN7nPglLGQ8XSAilfMk0U8vefA4/gT6fPDVJwxk5NbHzgHbofib9b8oeK/
orXVrglzwLqYlsa1QMURN2Yn3tTNAbnH3H6CRMNnP/gXnW+eg/mLEVPCAjSs
KLDhy93OmN9d3RN6W2m4O5+bya5vDg7fqHt8Q5SGxb5Lb2GjzEHC+zeBgbto
+M1BC20zkXkIr2CubwEajgwfZF/fMQ9+iulObto0bNvl0J4rNw8K0sFLm/Ro
mODqZbm6dx6uOenPCRrSsEnSXY+sY/Pw1bDJq8uUhscXKpOpfvPgoiOnneRE
wwVmI0ICgQzHn3vb40zDXiVsKWqh88BWFFfI5kbDy5dPpQTFzEPoTh8JIw8a
5hr7nULMm4d/HXsDzP1ouEObKLL75Tx4hp/VAH8afpCxJ9W0bB5ubrWKEw2g
YWHrG6nJtfNQSr1DqgxkrP8zOW3nt3kw5OB++V8YDVN3q4ke/zEPypyKN/ki
aLjsvnWa59g8ZNDObkq5TcOHDXPTymnz8I07RzY+koaZ89pE+xfnIawHf9q4
R8NNrLNpa2vz8PDG4MEL0TRsXKeZrstFh6GAIAuOWBreLH5BzIVEh90aevzH
42j4643b6fe20CE0OTE1/AEN26l3p3dL0kHlS51sfwJj/UXOT+wP0eGPbmgo
52Ma3s8TLR52hA6LXc8F/zG85FbyJFefDg83kxJ+JdPwLdm1JzOmdFjfluL0
KJWGj4ZJSmy2ogPhX0/t+TRGvH7qZ6ja0+FKlf82yXRGfNIeZNxwp4NbVHh7
2BMaNl99I5HhTQerKP19chk0LHR2MOO9Px16D8rkYIaHypl3/AmiQ+o/O3nj
TBrOFtiZyRVOB53oIw3dDLt4n9ihfI8OC9PC/ieyGPH95J15Ko4OojbHDGsY
piol7bj6iA5bLRRBJpsR33vVmUlpdPgdqnIuhGH/iZ87qrLp8Ku0MesLw5r6
HFk/njPm+9VVSPwpDTM9U5JkKaZDPi29yZbhRmbTLNlyOmwXouYnMBxp7y9p
WEWH827UjlqGjWvSsi7X02FlRWP/KMObRd5Lxn2gw7dv3sOrDH/9bzyrrIMO
Csw6fVzPaDi5l0eqr4cRP1NbCRLDtqr7slf76bDp3s1qboYl4y2lxEfooCXu
WrbGGP+bGpit/YcO13aPcf5m+IVxttTFv3T4FNlY38CwR8GH7LtzdCCmj44+
Zngv14xU4TIddvhKeDsxvOjM//TTBh3qcx295Ri2inzn3Mu6AMJudsYjjPVX
FdgpDHEvgOrz/9bjGBbvZKX8Ii/AfYfLtw8yHELLL5naugCt11f/9jLiO8Z/
2m9WdAGEiqcOXGbY4MDSwWWpBYiNIzovM/aH97puHdvuBdh7TzN2gbGfV9Im
w3gOLEBSIle8M8M9ONZAQHMBLkjfju1k5MNjwmDnDoMFkJcwjQlj5Ms/2dD4
ncYLIFebmNzFyCfbY7ssVMwWwED1Spkgw9L3rw5p2i/A7RWhnaEpNHy7WDhL
5+IC1OZujipg5Odkd52TofsCbP9M29zOyN/irbx/LfwXgNa/ljb3iIb5NcqK
bYMWgFkxomA5kYb9rK18ncIWQORzwezCQxo+lJmz6h27AKK8XRq98TSc/v44
/i9xAZyErPUrGOeH6fdsSHDqAmhbZ2QkMM5Xkzxw3c9bAJalEmOVGBo+Xda7
7XntAmx9Mq6ydpeGX38NHChqWoBlsxX1lDs0vG1ZKqO8bQHSzohV7Gac96HD
XnINfQtw3PWdhjajPrg2cxz4QV0ATdkjXRuM+tE28XJ5bGEBKqvqDljdoGFl
4pnq6X8LwOMrcaSAUW/mTTJ1VzgWwYXf4JzaNRoO7j9oKrhjERrNxI5OetLw
6L8fW0TkFuFYgL4PM6OeHRWP6JdUWoSwQwE+m91pmOjYfX73wUUIzPMjCLnQ
8KO/bl5GpxYhR+gI7wc7Gn7JlBp969Yi5B6FtacnaJitbt9ZiduL4KC+e+yi
EWN/brVK46hFeC7GRZA8RsN8zKuV/x4tgvCzmyuhuoz6wmI1ee3lItjFXrBY
1KDhPWzb9D2+L0L0VX+dDFnGeWoq2swzsgjfByqWF6Vo+GeE/lD+70UgzIeo
H9vBqB/s165O0Bbh/XyAY58wDc9yfM125FiCyaY4jheMflPGncB0TnWJcV+U
FBdfpmJim0LbssYScBH/9ZAZ/ckx6v2jJLQEK6OsLuuzVCzIM6f81WgJ2l7t
O9Q4TcV+vKdtTl1YgvdSxkpcjH6nTiZV6McsAX/XJRa/Riqu3BLls39iCdaZ
1dn/3aNiwkT98nHKEmzVODvz8zYVG1cu33SiL8GNqi0e70OpeNjOOSpx0zL4
iOuZXblOxWzPtZ8tbV8GYzXmiFtuVGxyaPlrleEyUC9/6f+P0d9/213U0CtY
hqYCD/HVZQre/BxtMvdcgeovfYeHi2bwqe2uMWxXV+B1qNfu588Y94+7caJv
AlaA72/gPreUGczj9vPgtogVkH9XndITMYM5FMK8+1NX4G20uOFxmxm8VtA8
avtxBTSykpu7OGbweLFxk7PMKnD7BrpGWv3F+K31vf++r8LWlQ88FlNT+I3u
4Wrui2sQbRX6b+7BODZfMu70mFgH441/vhMyY/iZwK8Z111MiFV63bzu3wjW
m2L+aPqXCYn6q07NBwzi2OAz2VpPN6EQT/rqOK0Xv3xeOx1nwYzEO/4pvZX4
jI1+qG5NtWJGNwg23Vo8n/G4QP6RHBtmdM7p4WTXcjeWDH6YWOHAjB4OCtwU
/dyNE83d0KgnM9K66OxlGNGNb7IIJOy/w4wkF3cwrc10YRMb50O9b5nR9xvU
+faWT3iWxBspIsSCajrtnl42bMeRgjfezoqwoBzH4i9UhXYsKTT5+4M4CzKz
LPp6l6cdn5Jq0rkqw4LmTHi01j614WLVoLX23Szo5a5gqwnLNuxhNeN1U58F
3esgnOxxb8WT2YzrtS8LWtX5lTz4tAWPHIiRyetkQSOfLnF8NGvC3g3H91R2
s6Brly5Y3znYhJlPcx3u+MKCkq/4/bMVa8KSl8PN5r+zoJJ8JaYz443YMft6
KBpnQQor7ooOgY34N5/LSB8TAVHjLsUzFTTgqXGUxrmXgPhMazmZhN/jG37r
eSKqBBQh6FldzvweEwlVZSrqBNQzRnoUNlmPFcXV2s5oEdDUsEO8U0U99jij
tJppSECn79TGHDpXj6l1289qOBBQy7iTs3JmHaYnzwpcSiCgcKLOrdHDtbgn
0qve8xEBvdelXOzYWYvLAmY8fZMJ6KI5Nuzgr8W+VpMfAzMIqNhIsm57D8az
238GxxYQkAm5Q4hDBWNKUtf06wYCuh85Jp05XY0nEosbmOgEFKfmWCgVUYk/
RCh7sy4RUNvRqini+Uqc51cgzrVKQI/s6nRENSuxi3leAP8mViRN7Bwuma3A
vwUzdsvysKJGJxPFovMVeDQhJtVIhhXt+ubb+fXoOzz4wMM38Qwraine+d8B
pTf47akxGWdLVpRmH2DCxPcGPyBZ96qdY0USce3BE9RyfPT+MY1v9qwovWQt
efvrclx8V5pJxJ0VPW/aGzCjVY5DgvvvZYawohe6B/5ZWr7Gsh5HswuKWNGH
pGHZzMxSzKRUbRZYwoqif328oHy7FPdP7WMzfs2K+JyynfrcS3GM6w5XSgUr
cjY2vByiXoqXHP8p721iRZbOpUixqwR/PFdS8WaAFRncaulvYC/BnsfEuus5
2dBBKuc//7vFWDvul/0nIhuaVjSuj/AqxgLfXlAG+djQZuN/pqUWxbjCRY1n
WZANHR0P93KRLcZsd0wMVCTZEDch3O5mZBFOawqpSdFgQ+/8kqq2kF/iVp0/
L3wvsSHBHkOz+mP5OP1e4cEQDzZkGbcpUUE8H3t99vkQc4UNFX8Sp7TOv8Bb
HJh+v7jGhlpS3wm9yniBL4Rs3zESyoYMU57Wsq8+x8u1Ro9OpLKhXz3579+9
y8NyWsWhsu1s6Ey+zVSnTQ6eN/PcwvSJDU3NhK+Kaebg+kvKz791s6F124BH
SUI52PZxQUdUHxtie6n1mavvGU6Yz9s+N8oY7zrC9MfsGWbKzyzCK2yodndg
+jnzp7hva/x3y53sSO8eLV36WhbOUT7tsVeBHclpH936n1UW9tEjbyIqsyPe
sH6tucNZmNcnRhbvY0evXHoIlqxZWLfj3hUZYEfLoaUnfzzMxK/Cwthnz7Aj
gYMihxuqM3Ak7eq+e6HsyDFAvKTtZDrOdaM7m0ewI04/iyE9iXTcMOqbuuMu
O3pyifRqnZqG17/4sL69z47OPLqLlOPT8JWKK72/HrMjusrV9sD+VGwZevm6
VjE7ik//z+CZbwqWEXR6TxtgR8Z+mqamX5Kw9v1fi1U/2JGRdGaSTmEStmV3
VLzzkx3xNH0NvRSehJOWLzwUG2dHFSU7zM8fSMLEIXtnozl2tPUpXXMo+RGe
zbHmfsbJgaROF6iPX07EtepmppYHOFD9jZKsdI0E/FVDysn+IAf6/C9m/cCW
BDylOevnosmBtES+7OaYjcdbjsQm+2tzIL7cAyumL+Kx+7G2kSRjDpTRe93y
rmg83mql6/XtIgdicdXCSrwPsMf1/dFWjziQbm6eWuCWWBwWyPLkQjIHctwS
NW9Ki8GPb3YXu6VxILdj3hK85TG4IdSrJyCbA018dnDLQjFYKKpAOKWIAw32
6PZoWt3HjanSL743c6CiUNX5j+lRWKRGoNl6iTGf2A75ede7WDFmYGD3KgeK
V07PKtK8izXtn84S1jlQxWAWTwnfXWzNvF+skIUTBXYJGUWW38FpR02vrvFy
ok3b9S2C2O9gsc44qXQZTmQ7XXNFpjwCSwzzBv84zYlY/Jh62HXD8O5XXx+W
nuFE4SqPB+bFwjAKSc+/bcmJ/O35/oquhGJ7aeVeFVtO9NHRZkjyVSjOdDmh
HOzKif4dcJak7AjFkrR7AxLBnOiZQ6+2EW8I4/8558ELhZxIx7dAZ44vGCta
Syj6FHOintA+ubChm3h/uZp4WCknuufLrZGUdhMfdbvImvOOE31dKP8mIXoT
u3S/75po5ETJso0oY0cQfpkV5HZliBOdVuLVq4Ib+KDOfEowHxfKS3zM/K7V
H2unccfEbeZCji/VuD4/8seGi5IhWYJciM2kSN/Y0R9b559yaRDmQhnIc9Zp
/RoO4i/az7GTCw22u2zhVLvGyGfX9hjEhSQHChIfv7uKT4b9WE+/woWS+h+y
jW5449OFFPp1Xy5U3no2+fU3b2z+dWPa8hoXGud8pdhe6o2td0p83xzIhR42
x6a2uHhjtzb7NxF3uFDWOMxNfbmCwwR+eno84UI+glWRAVVeuOLpr+HD7Vzo
17I984LnZVzTPt8r9IkLSQjciDLadxnXLRA6F7u5UARrPmHzojv+YCBT/aqP
C3lftBueuemOv0w7PZb5xYV2tOy40JN4CVP3/znFs8qFXm1XiS7pc8WyjRP1
Azu5Ufuo+SR33kXsb1JxLl+BGwWaftwr8t9F3Po9cv4/ZW40P7iSG2Z4EXvN
Kshu3c+NJH8Kdp2bccKV4h53TwE32h2QfIZfwwmbBMyebDrDjeqIze4KIw44
m/X9eIIlN0q/fqPy+hsHvBAbf8vhHDdKm43wFr7vgFNy95cynedGaoHLzPqa
Dvh3j9+WQ+7cqN6rUg2lXsDXVVYHikK4Ebf7maxjbudxR0Xr1aBwbqRXTTU2
1TuPJY6m8p64w42WvtKsn0qcx03Wh9FkNDcyqNMoE+qzx6TIm9nSj7kR39vz
24wM7XHOLxa3x0XcaKOtzTnHzxZ/ekxcDBngRlUnnBQ7Hlth4ZeKl+//4Eb4
o+uXbees8MX646OPf3KjF9KSTpMiVnhtMqqzeJwb0RYPKYplncXymjx5Q/Pc
iEOsyKagxBKHDfGc1eAhonNFMc13psxx56xSpx6JiPh1LRpby8yxELux3il+
Iuo5/HlfRJA5fqVyf4/LdiJyObSvyobfHA8F83IlyhCR9Dm2Y0HaZ7C6FF8l
7TARJXzaHVv72hSHqqnsWUNEtCpQq0WKMMUdRidzOXSJqNU4zKfF3BQ7+cbE
ixsSkbbrBbmry6fxg0Y+9xPmRGR6YvtGss5pPO1MEn3uQUSJp48XeU6a4IwC
8k37J0SkmvY6aGn5OP4dI7+slkVEUu4g11R+HCv66PjwPSOiLa8PO0v6Hsdv
D151rnlBRJF9HXYhVCP8qanvpEg5EXXrO8bVUwzxxvATib42IvLJrX1vzX4M
6za8TS7qJKJdzcdaE9sMcGRul8DtbiJq0lzjPvnAAG/xYOE80Md4XtmizQFx
A6y8epEWP0pE2RVRUY1a+thWULn+5AoRzXQ0jAom6OGnS0c15daI6I5/9Hy9
vR6e/G5Xvr5BRG8jZ0J/KOlhv6y4F4WsPOj7ru5chxZdHK1Cf8BN5kHbZQdF
g1l1cdWxKodmOR70+NvnduM4bSwcZMgKZ3jQgvSyuusoYKvJEmtTSx70zIfp
08Z9wEnmwmXO53jQl5gA6/OHAAuqTF+IPc+Drp8uPhKUqIX5hqNrRy7zoFVy
+OjfM4cxi3bX9bAIHjQZMJstvaSBtV8e/Jx0lwfFjX/ITynSwMFCWfKFUTzo
Zf2tj/XOGnht1vvblwc86FXn6Bbh/oN4MVtATe4Jwxp/Wz40qeMpVsvZljc8
KP3N1pcNNQewvHftsaEKHlQ1F/kkIvgAdhnamTlbzYMcHK4ODGsfwGPlyyeF
G3hQOb9YjWWLKv7hnFLo/okHaZy1O7p7cD/+3DLkwjfBWO8vbGomvg9XRjv9
MBPiRZt5/m7JIezGBtb3+6xEeVGHnT4hpkgFf5V/02UvwYv4Ag0O5Jur4Nlm
jgZ3WV7EUr0Y7JyrjOWZC/LC9vIiNjP2mjenlXCS3+yVMkNelKbuXk/olMey
esKXKk7wIgfSjkDaLXlcxq/rWGvCi8wk5jKSVOVxZ/FD8zZzXkR8/OYD75Nd
mG1K/dCvC7yoaVkk0CVgJ75qd5MgcJ0XhfbfijE5JYuZlfPWtgfxorXhA4tk
Plkc++/TgvgtXvTl/K3L3u0yuOCx5ITCbV70o9FnmfO4DP71ubFdJ54Xrfur
W140kcanDYiPfPJ50VbDfT7XfSWxyp7Huz5/Z/zexUDbURDDbzsqZJ8N8aJN
Ur2xlFFRjNwHpK6N8DL6sdCd2jui2CRXXEz4Dy/6qsmTktcjgr1FczY7zPKi
c6N+/1leFcblnKWrsxx86PsTEWOfvm1Y62db++YDfKhZW8lAfZsAbro58/GX
Oh8yVZGYFejixydFSR/KD/EhUiy1779IfmxvaVpvdYQPrYpmPfqyvhmHdHwr
zzrOhxZu16nk/CXjporfT/Y68KGB85751xr4sHE8k/fpWD70RKxF+bcaB+58
d+X4+Xg+pL4yV99QyY5Nhn/KeiXyoV815aOXEDs+rdT4PTqVD/2IijUiHGfD
5s139Vry+NDO2W0bc5cI2G5183ZUx4eeRhyPS25lwl4OMrVKs3zo2fG2gV0j
czUt22VuE+RJiOl2eAjxda9WhIDVwhUFEjrnfZSFo7RfS5cvxumHIglZWLQ8
HXg1qFVLWNapUCGhtRdHWtW/jWi9pbQxeamSUGz6Huul4T9auY2+Ad8RCe2K
U1XgIc5qhV9puFxiQUK7+9SmvWvXtbTdlwfEz5LQtC+pSjZ+Q4vJWfl4tBUJ
NbdfaPTexwQ3rJPkXWxIjPvVjhefODaBr777H1EHElqtNFLzrGAGB1GB83c9
SGjTt67L4ZpscKTFwcw+nIQoY+ldE308IEv3T2iKYNhNX4psyQvEHfd7FO+Q
kJCCsMHBPl7o839rthxJQuyxNp3p3/nAYyfPmQexJHSb9VdRxiIJkm+Xn2lI
IaHnJicprHL8EFza+lA+jYTe/Ph0stqTH5x+DH+JTSchG++DtiLv+GG3Gre5
TSYJVcvnfBk4IQBNv+3MF3JI6G1+1J/XNwVhVo/TYmcJCR0w3LxTm30b9F0R
e3S/lIQkzlqdaDm3DarT9vXOl5GQjoLaO9lX2+AO3cai7g0JdXPbFF2x3Q5i
z0osrKpJqE0r2zHnvRAYEqwtoz6QkByXUsW9SlHYvftK0mwLCdkVzOeQxcVA
0Dqiz7KVhDQupu3JCRWD4dJiS5kOEkLiFYrcJuLg58B6tuYzCbkN1TcJXZeA
7PqXZ6lDJKTswBnSfVgSCKWXOK4OM75vdryBL0gSLmbvfLM8QkK3sl+/EcCS
IB+WJcAyRkLr359u3NKRghK9xE7BKRI6b3DPF5lLQ+2HQL1DiySkSiuKlciR
Bcl3GvN4iYQ+fN/z+uy8LIQ9X8zSXSEho8tb1GN15EA/8gqT8RoJtcgIqRHG
5KDDyKnSnoWMdhZ3XPHevwsGO47vieAjo8MzXWeWORRBC3MOc5HJKLZ/cP63
iyJkFDXdj9lMRifmS54LfFQEx1g0nSRIRrmlW7foxCrB9Kn9ufnCZPQ343uW
r6IKrPQIi3bJkVHm2T1ES7M9YN3Y13pmFxkRTLXsg2L2QM3rhwH98mQU5L2j
z6N1D9xK5OsbVSKjxPlU5wq9vcBhyRK/sI+M2uWFvp3S3Qdbv09yiCAyUnzf
fvm0uyos7H36LuIIGVnK+bHnF6vC10gbN5o2GaX5V1SO0FUhQeNTa5MeGXXz
+1yjhxwAUvLr+1eMyGjx2lM9njw14DobLNBsTkYh8crOkxIaMFF8sHGPJRn5
Bpsqv/DVgA8cc1dTz5LRDL8ld1CLBkS8deq9Yk1G/11ot7nidwhYthkli1wg
I81M82n6oCb8+yoo4X2ZjF4MN+sL/gH4rtz5acCDjCRUfBGTHIKKiDu39L3I
6Pvl+N87ziHwP7D6U8SHjCpiJ+q8GxHQHw4/a/YnI/XbBb/vZx2BGbN8BdEw
Rvykr9xdCtaB9gLHgdvhZGRdxB1pV6sDBQSx6NkIMuryu7uQyqQLl8piZ5rv
kpEkXfVEbqgujPP7vfKOIaOwwnUfkwQ9GOlGah+SyaiTp4PY+lMfNCmZ5qdS
Gc+7h9LeqBlAEjeLX38aGXE/q5f8EW0AxrqNZdMZZLTyVVJm9vAxqHp9bC85
l4yyC5VUygoNISnptJJVKRnJOR+yO/rpBMyXlR4fLSMjoWPM/caHjOFkl4C7
ezkZte4+/nM81xjYuHpfBL0joyXBL2qT4SfB9/q5ndmYjK42avpNvjCBk7aO
UtOtZPToGU5Pe3YaXgQ0al9tJ6P4ETrr4YnTwPZI9sJ6B2O/440HBpRMoapz
PIPUTUZ5KgcKlSpMQf7IZTHVPjLSZd1bd/uHGUTYdByu+UZGDd7qCWkKZ2Dk
PxUb/e9kdK1bxeuL/xlIKqGlnB0io5KGghPVW8yBTdpve9AvMhL5MXaPZm8B
F1CvOvtvMjJ+P+mWXG4B1dbqlrF/yGikqizkHo8l+D5cScyaJKNC3zGqFbaE
n2w3BZqpZOR+gXBId58VVE/c5iOtMeIfbrp5/3kbEJPa0ynM8BfvVDUPRxsI
su6/L8fwjZCayjxnG9DqVOQFhrWvVQ6LedhATVkX0YNhwr+v6XzXbQDfFOZq
Y7i9OJLO/NAG6gSLCHfWycjAoNd6T4sNSJ60bIhn2HbV2udamw2E3NkU9oRh
lcRdk9WdNqC9asryhuFb0Y1LRl9toH54adNvhts4jFLcR23gfb72hs4GGdVq
Pd3AazYgPTZVc5LhfAKTBHGTLYSJPQw6x7Dz+XOZlgRb0I37s+bN8Lv90DzL
ZQsNflH/MhnW/LyRobDNFqSLVasKGV7/lvoxQJgxfmLo+juGB7USz3wUY4y3
3rP6ieHdqW0xbjK28PRhf8UAw0I9ElC50xYInaEB4wzzTmycJCragiOHksY8
w6mBV+ttVGyh8cjX5Q2Ga1SjE4v22sL/AJ3X1pM=
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0, 4.5}, {-1.2, 1}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellChangeTimes->{3.4166022916706*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Two decoupled oscillators", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "5"}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"HH", ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"bp", ".", "bm"}], "+", 
        RowBox[{"0.5", " ", "be"}]}], ")"}], "**", "be"}], "+", 
     RowBox[{"be", "**", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"bp", ".", "bm"}], "+", 
        RowBox[{"0.5", " ", "be"}]}], ")"}]}]}]}], ";"}], " "}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"vals", "=."}], ";", 
  RowBox[{"vecs", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.4166024484974003`*^9, 3.4166024536766*^9}}],

Cell[BoxData[
 RowBox[{"bp", "**", "be"}]], "Input",
 CellChangeTimes->{{3.4166025568862*^9, 3.4166025651386003`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"HH", "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.4166024930198*^9, 3.4166025193058*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "2.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "3.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "4.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "5.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "2.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "3.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "4.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "5.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "6.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "3.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "4.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "5.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "6.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "7.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "4.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "5.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "6.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "7.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "8.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "5.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "6.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "7.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "8.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "9.`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.416602494455*^9, 3.416602519571*^9}, 3.4166025749666*^9}
  ]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"vals", ",", "vecs"}], "}"}], "=", 
   RowBox[{"Eigensystem", "[", "HH", "]"}]}], ";", "vals"}]], "Input",
 CellChangeTimes->{{3.4166025988814*^9, 3.4166026001294003`*^9}, {
  3.4166027018258*^9, 3.416602702715*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "9.`", ",", "8.`", ",", "8.`", ",", "7.`", ",", "7.`", ",", "7.`", ",", 
   "6.`", ",", "6.`", ",", "6.`", ",", "6.`", ",", "5.`", ",", "5.`", ",", 
   "5.`", ",", "5.`", ",", "5.`", ",", "4.`", ",", "4.`", ",", "4.`", ",", 
   "4.`", ",", "3.`", ",", "3.`", ",", "3.`", ",", "2.`", ",", "2.`", ",", 
   "1.`"}], "}"}]], "Output",
 CellChangeTimes->{{3.416602437281*^9, 3.416602454987*^9}, {
  3.4166025960578003`*^9, 3.416602600535*^9}, {3.416602699517*^9, 
  3.4166027029802*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"vals", "=", 
   RowBox[{"Reverse", "[", "vals", "]"}]}], ";", 
  RowBox[{"vecs", "=", 
   RowBox[{"Reverse", "[", "vecs", "]"}]}], ";", "vals"}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1.`", ",", "2.`", ",", "2.`", ",", "3.`", ",", "3.`", ",", "3.`", ",", 
   "4.`", ",", "4.`", ",", "4.`", ",", "4.`", ",", "5.`", ",", "5.`", ",", 
   "5.`", ",", "5.`", ",", "5.`", ",", "6.`", ",", "6.`", ",", "6.`", ",", 
   "6.`", ",", "7.`", ",", "7.`", ",", "7.`", ",", "8.`", ",", "8.`", ",", 
   "9.`"}], "}"}]], "Output",
 CellChangeTimes->{{3.4166026809062*^9, 3.416602709657*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"N1", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"bp", ".", "bm"}], ")"}], "**", "be"}]}], ";", 
  RowBox[{"N2", "=", 
   RowBox[{"be", "**", 
    RowBox[{"(", 
     RowBox[{"bp", ".", "bm"}], ")"}]}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "\[LessEqual]", 
    RowBox[{"Length", "[", "vals", "]"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Print", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"vecs", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ".", "HH", ".", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], ",", " ", "\"\<   \>\"", ",", 
     RowBox[{
      RowBox[{"vecs", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ".", "N1", ".", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], ",", " ", "\"\<   \>\"", ",", 
     RowBox[{
      RowBox[{"vecs", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ".", "N2", ".", 
      RowBox[{"vecs", "[", 
       RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}], "]"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "1.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`"}],
  SequenceForm[1., "   ", 0., "   ", 0.],
  Editable->False]], "Print",
 CellChangeTimes->{3.416602729391*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "2.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`"}],
  SequenceForm[2., "   ", 1., "   ", 0.],
  Editable->False]], "Print",
 CellChangeTimes->{3.416602729391*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "2.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`"}],
  SequenceForm[2., "   ", 0., "   ", 1.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027294066*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "3.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`"}],
  SequenceForm[3., "   ", 2., "   ", 0.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027294378*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "3.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`"}],
  SequenceForm[3., "   ", 1., "   ", 1.],
  Editable->False]], "Print",
 CellChangeTimes->{3.416602729469*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "3.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`"}],
  SequenceForm[3., "   ", 0., "   ", 2.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027295002003`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "4.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`"}],
  SequenceForm[4., "   ", 3., "   ", 0.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027295313997`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "4.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`"}],
  SequenceForm[4., "   ", 2., "   ", 1.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027295626*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "4.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`"}],
  SequenceForm[4., "   ", 0., "   ", 3.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027295938*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "4.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`"}],
  SequenceForm[4., "   ", 1., "   ", 2.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027296094*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "5.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`"}],
  SequenceForm[5., "   ", 4., "   ", 0.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027296406*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "5.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`"}],
  SequenceForm[5., "   ", 3., "   ", 1.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027296717997`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "5.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`"}],
  SequenceForm[5., "   ", 2., "   ", 2.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027296874*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "5.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`"}],
  SequenceForm[5., "   ", 1., "   ", 3.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027297186003`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "5.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "0.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`"}],
  SequenceForm[5., "   ", 0., "   ", 4.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027297342*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "6.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`"}],
  SequenceForm[6., "   ", 4., "   ", 1.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027297654*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "6.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`"}],
  SequenceForm[6., "   ", 3., "   ", 2.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027297966003`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "6.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`"}],
  SequenceForm[6., "   ", 2., "   ", 3.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027298122*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "6.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "1.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`"}],
  SequenceForm[6., "   ", 1., "   ", 4.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027298434*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "7.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`"}],
  SequenceForm[7., "   ", 4., "   ", 2.],
  Editable->False]], "Print",
 CellChangeTimes->{3.416602729859*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "7.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`"}],
  SequenceForm[7., "   ", 3., "   ", 3.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027298901997`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "7.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "2.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`"}],
  SequenceForm[7., "   ", 2., "   ", 4.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027299214*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "8.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`"}],
  SequenceForm[8., "   ", 4., "   ", 3.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027299370003`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "8.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "3.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`"}],
  SequenceForm[8., "   ", 3., "   ", 4.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027299681997`*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "9.`", "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`", 
   "\[InvisibleSpace]", "\<\"   \"\>", "\[InvisibleSpace]", "4.`"}],
  SequenceForm[9., "   ", 4., "   ", 4.],
  Editable->False]], "Print",
 CellChangeTimes->{3.4166027299838*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Two linearly coupled oscillators", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"HH", ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"bp", ".", "bm"}], "+", 
        RowBox[{"0.5", " ", "be"}]}], ")"}], "**", "be"}], "+", 
     RowBox[{"be", "**", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"bp", ".", "bm"}], "+", 
        RowBox[{"0.5", " ", "be"}]}], ")"}]}], "+", 
     RowBox[{"0.2", " ", 
      RowBox[{"xx", "**", "xx"}]}]}]}], ";"}], " "}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Eigenvalues", "[", "HH", "]"}], "//", "Sort"}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.994936153078953`", ",", "1.8893633682901467`", ",", 
   "2.0903812873466294`", ",", "2.783793648352695`", ",", 
   "2.9848133028390476`", ",", "3.185828277032211`", ",", 
   "3.6784083006863653`", ",", "3.879664192672512`", ",", 
   "4.080577889706295`", ",", "4.281354888560959`", ",", "4.578235946818142`",
    ",", "4.79004371067676`", ",", "4.993260906326349`", ",", 
   "5.1853955991935985`", ",", "5.37862879618715`", ",", "5.517149799942078`",
    ",", "5.840462755951073`", ",", "6.17130910016814`", ",", 
   "6.451263772983425`", ",", "6.522980903783384`", ",", "7.039747387290595`",
    ",", "7.479984099056089`", ",", "7.6631915412069675`", ",", 
   "8.456873102485401`", ",", "9.082351269365029`"}], "}"}]], "Output",
 CellChangeTimes->{3.4166027921809998`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NN", "=", "10"}], ";", 
  RowBox[{
   RowBox[{"Eigenvalues", "[", "HH", "]"}], "//", "Sort"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.9949361530051244`", ",", "1.889363344005036`", ",", 
   "2.0903812680154537`", ",", "2.7837905350049574`", ",", 
   "2.9848084590153694`", ",", "3.1858263830257867`", ",", 
   "3.6782177260049798`", ",", "3.8792356500155494`", ",", 
   "4.080253574025908`", ",", "4.281271498036175`", ",", "4.572644917017692`",
    ",", "4.773662841055878`", ",", "4.974680765074097`", ",", 
   "5.175698689061935`", ",", "5.3767166130516815`", ",", 
   "5.467072109050465`", ",", "5.668090036090773`", ",", 
   "5.8691079624409905`", ",", "6.070125885134564`", ",", 
   "6.271143806092159`", ",", "6.361499357496935`", ",", "6.472161728388311`",
    ",", "6.5625174924063945`", ",", "6.763535669840103`", ",", 
   "6.964553618336755`", ",", "7.165571320973971`", ",", "7.25592872147692`", 
   ",", "7.366589023732574`", ",", "7.456956135609446`", ",", 
   "7.56760685722762`", ",", "7.657989087121777`", ",", "7.859014646769058`", 
   ",", "8.060025076825303`", ",", "8.150409080964627`", ",", 
   "8.26102890303782`", ",", "8.35169116927618`", ",", "8.462037460893082`", 
   ",", "8.553209307591167`", ",", "8.663052356532308`", ",", 
   "8.754629147438829`", ",", "8.955580930024846`", ",", "9.045595821505835`",
    ",", "9.156152263621577`", ",", "9.250724609985896`", ",", 
   "9.356736074541192`", ",", "9.46058329161425`", ",", "9.557540554295796`", 
   ",", "9.671416739877829`", ",", "9.758503095076817`", ",", 
   "9.87603346575972`", ",", "9.944843051130105`", ",", "10.071373808513522`",
    ",", "10.171706302194757`", ",", "10.262366125767288`", ",", 
   "10.427859197953074`", ",", "10.455873369118015`", ",", 
   "10.653737765685493`", ",", "10.705828141308826`", ",", 
   "10.85334372501691`", ",", "10.854018440990187`", ",", 
   "10.989321564524749`", ",", "11.12808225018947`", ",", 
   "11.262421071080757`", ",", "11.460902899519677`", ",", 
   "11.513440306829336`", ",", "11.73958201021266`", ",", 
   "11.777746622964939`", ",", "11.824337414477695`", ",", 
   "11.948342451114085`", ",", "12.133704653107543`", ",", 
   "12.186221693984407`", ",", "12.519369227254682`", ",", 
   "12.570250805132506`", ",", "12.727471634928964`", ",", 
   "12.804635736839078`", ",", "13.030823530220221`", ",", 
   "13.038403799518218`", ",", "13.210296092582027`", ",", 
   "13.463924202733317`", ",", "13.725972955360103`", ",", 
   "13.774508815969916`", ",", "13.836107208425428`", ",", 
   "14.121200399696711`", ",", "14.337948329618737`", ",", 
   "14.392732425982903`", ",", "14.76907311986579`", ",", 
   "14.858414800870971`", ",", "15.095134348799153`", ",", 
   "15.190187734417588`", ",", "15.736135345050474`", ",", 
   "15.756532680414207`", ",", "16.016432530142957`", ",", 
   "16.25077618922392`", ",", "16.54951139315629`", ",", 
   "17.268393244028946`", ",", "17.322718718811423`", ",", 
   "17.52646084157879`", ",", "18.148877334690525`", ",", 
   "19.242121486747262`", ",", "19.484630979787855`"}], "}"}]], "Output",
 CellChangeTimes->{3.4166028065486*^9}]
}, Open  ]],

Cell["Try the conventional method for the ground state", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Apply", "[", 
   RowBox[{"Plus", ",", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      RowBox[{"Eigenvalues", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"1", ",", "0.2"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"0.2", ",", "1"}], "}"}]}], "}"}], "]"}], "//", "Sort"}], 
     "]"}]}], "]"}], "/", "2"}]], "Input"],

Cell[BoxData["0.994936153005124`"], "Output",
 CellChangeTimes->{3.4166028439105997`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{"E1", ",", "E2"}], "}"}], "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"Eigenvalues", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"1", ",", "0.2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.2", ",", "1"}], "}"}]}], "}"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4166030776298*^9, 3.4166030822318*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1.0954451150103321`", ",", "0.8944271909999159`"}], 
  "}"}]], "Output",
 CellChangeTimes->{3.4166030837138*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"0.5", " ", "E1"}], "+", 
  RowBox[{"1.5", "E2"}]}]], "Input",
 CellChangeTimes->{{3.4166030881442003`*^9, 3.4166030935573997`*^9}}],

Cell[BoxData["1.88936334400504`"], "Output",
 CellChangeTimes->{3.416603094197*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Bosons with strong repulsion U", "Subsection"],

Cell[TextData[{
 "Introduce operator ",
 StyleBox["n",
  FontWeight->"Bold"]
}], "Text"],

Cell[BoxData[
 RowBox[{"bn", ":=", 
  RowBox[{"bp", ".", "bm"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"U", ":=", "20"}]], "Input"],

Cell[BoxData[
 RowBox[{"H", ":=", 
  RowBox[{"bn", "+", 
   RowBox[{
    RowBox[{"U", "/", "2"}], " ", 
    RowBox[{"bn", ".", 
     RowBox[{"(", 
      RowBox[{"bn", "-", "be"}], ")"}]}]}]}]}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Eigenvalues", "[", "H", "]"}], "//", "Sort"}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0", ",", "1", ",", "22", ",", "63", ",", "124", ",", "205", ",", "306", 
   ",", "427", ",", "568", ",", "729"}], "}"}]], "Output",
 CellChangeTimes->{3.4166031442262*^9}]
}, Open  ]],

Cell[TextData[{
 "In the limit of ",
 StyleBox["U\[Rule]\[Infinity]",
  FontWeight->"Bold"],
 ", only two states remain.  Many bosons with hopping -- Bose Hubbard model \
(possibly, Infinite-",
 StyleBox["U)",
  FontWeight->"Bold"]
}], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"H", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"U", "/", "2"}], " ", 
       RowBox[{"bn", ".", 
        RowBox[{"(", 
         RowBox[{"bn", "-", "be"}], ")"}]}]}], ")"}], "**", "be"}], "+", 
    RowBox[{"be", "**", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"U", "/", "2"}], " ", 
       RowBox[{"bn", ".", 
        RowBox[{"(", 
         RowBox[{"bn", "-", "be"}], ")"}]}]}], ")"}]}], "-", 
    RowBox[{"bp", "**", "bm"}], "-", 
    RowBox[{"bm", "**", "bp"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.4166032647986*^9, 3.4166033417222*^9}, 
   3.4166033946062*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Eigenvalues", "[", 
   RowBox[{"N", "[", "H", "]"}], "]"}], "//", "Sort"}]], "Input",
 CellChangeTimes->{{3.416603396135*^9, 3.4166034276158*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.9999999999999956`"}], ",", 
   RowBox[{"-", "0.19803902718557276`"}], ",", "0.`", ",", 
   "1.0000000000000042`", ",", "17.928692494294513`", ",", 
   "19.999999999999854`", ",", "20.198039027185423`", ",", 
   "21.921215971661052`", ",", "39.415624397556066`", ",", 
   "59.933407243254194`", ",", "60.07130750570548`", ",", "60.0787840283389`",
    ",", "60.517616456796986`", ",", "76.81457802732554`", ",", 
   "82.78478746499121`", ",", "118.85815186158717`", ",", 
   "120.06659275674582`", ",", "120.06675914564693`", ",", 
   "120.12292001506674`", ",", "120.15270797608211`", ",", 
   "139.83369047876025`", ",", "140.97470797574448`", ",", 
   "175.66095954517968`", ",", "183.58689068552047`", ",", 
   "200.0625019576079`", ",", "200.06250455892692`", ",", 
   "200.10630803615152`", ",", "200.107138646355`", ",", 
   "220.18902947089825`", ",", "220.26308273562913`", ",", 
   "238.15065044287752`", ",", "259.70105425414494`", ",", 
   "261.5479171044182`", ",", "300.0600014850882`", ",", "300.0600015163132`",
    ",", "300.09167684033`", ",", "300.0916924349549`", ",", 
   "314.4709193215408`", ",", "320.15893768210987`", ",", 
   "320.16142417049576`", ",", "324.3238391510418`", ",", 
   "360.2665503679223`", ",", "360.41357601460874`", ",", 
   "397.3149418490047`", ",", "419.53576754055814`", ",", "420.058334143592`",
    ",", "420.05833414389565`", ",", "420.08286472060536`", ",", 
   "420.08286493906616`", ",", "422.21503861298515`", ",", 
   "440.12919274377475`", ",", "440.1292472651593`", ",", 
   "480.22420888622264`", ",", "480.2299950914972`", ",", 
   "493.24731198588063`", ",", "504.99220470488206`", ",", 
   "540.3526269664429`", ",", "540.6075890415315`", ",", "560.057143343139`", 
   ",", "560.0571433431418`", ",", "560.0770872528259`", ",", 
   "560.077087255254`", ",", "580.1114783677705`", ",", "580.1114792410418`", 
   ",", "596.3713657161935`", ",", "619.3381621543324`", ",", 
   "620.175182752659`", ",", "620.1753279476858`", ",", "622.9552590399929`", 
   ",", "680.3022340013201`", ",", "680.3137687989301`", ",", 
   "711.9927748900222`", ",", "720.0562503139361`", ",", "720.0562503139362`",
    ",", "720.128545205449`", ",", "720.1285452054714`", ",", 
   "725.5889703145564`", ",", "740.2248782950174`", ",", "740.2248783059005`",
    ",", "760.4460222802942`", ",", "760.8495042623504`", ",", 
   "780.3596038443477`", ",", "780.3596064448834`", ",", "835.3380236245796`",
    ",", "840.5612028296831`", ",", "840.5615254230931`", ",", 
   "859.1131095770388`", ",", "863.7547626945623`", ",", "920.886890422961`", 
   ",", "920.9072136808581`", ",", "970.7215506804707`", ",", 
   "986.1394289005079`", ",", "1021.278449319529`", ",", 
   "1021.8605710994918`", ",", "1114.3795006481864`", ",", "1140.`", ",", 
   "1145.6204993518134`", ",", "1270.9999999999998`", ",", 
   "1288.9999999999998`", ",", "1440.`"}], "}"}]], "Output",
 CellChangeTimes->{{3.4166034127802*^9, 3.4166034285674*^9}}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
AutoGeneratedPackage->None,
WindowSize->{1016, 696},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"Magnification"->1,
"PaperOrientation"->"Portrait",
"PaperSize"->{612, 792},
"PostScriptOutputFile":>FrontEnd`FileName[{$RootDirectory, "home", "leonid", 
   "text", "class", "seminar"}, "lecB3x.ps", CharacterEncoding -> 
  "iso8859-1"]},
ShowSelection->True,
Magnification->1.5,
FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 128, 1, 203, "Title"],
Cell[721, 26, 37, 0, 74, "Subtitle"],
Cell[CellGroupData[{
Cell[783, 30, 60, 1, 105, "Section",
 InitializationCell->True],
Cell[CellGroupData[{
Cell[868, 35, 66, 1, 42, "Subsubsection",
 InitializationCell->True],
Cell[937, 38, 841, 30, 72, "Input",
 InitializationCell->True]
}, Open  ]],
Cell[CellGroupData[{
Cell[1815, 73, 66, 1, 42, "Subsubsection",
 InitializationCell->True],
Cell[1884, 76, 128, 3, 43, "Input",
 InitializationCell->True],
Cell[2015, 81, 366, 10, 43, "Input",
 InitializationCell->True]
}, Open  ]],
Cell[CellGroupData[{
Cell[2418, 96, 74, 1, 42, "Subsubsection",
 InitializationCell->True],
Cell[2495, 99, 333, 11, 46, "Input",
 InitializationCell->True]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[2877, 116, 46, 0, 105, "Section"],
Cell[2926, 118, 57, 0, 41, "Text"],
Cell[2986, 120, 68, 2, 43, "Input"],
Cell[3057, 124, 64, 0, 41, "Text"],
Cell[3124, 126, 570, 18, 43, "Input"],
Cell[CellGroupData[{
Cell[3719, 148, 168, 4, 43, "Input"],
Cell[3890, 154, 802, 23, 150, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4729, 182, 60, 1, 43, "Input"],
Cell[4792, 185, 781, 23, 150, "Output"]
}, Open  ]],
Cell[5588, 211, 107, 3, 43, "Input"],
Cell[CellGroupData[{
Cell[5720, 218, 84, 2, 43, "Input"],
Cell[5807, 222, 749, 21, 134, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6593, 248, 124, 4, 71, "Input"],
Cell[6720, 254, 749, 21, 134, "Output"]
}, Open  ]],
Cell[7484, 278, 26, 0, 41, "Text"],
Cell[CellGroupData[{
Cell[7535, 282, 134, 4, 43, "Input"],
Cell[7672, 288, 771, 22, 134, "Output"]
}, Open  ]],
Cell[8458, 313, 47, 0, 41, "Text"],
Cell[CellGroupData[{
Cell[8530, 317, 218, 7, 43, "Input"],
Cell[8751, 326, 949, 31, 208, "Output"]
}, Open  ]],
Cell[9715, 360, 50, 0, 41, "Text"],
Cell[9768, 362, 112, 3, 43, "Input"],
Cell[CellGroupData[{
Cell[9905, 369, 126, 2, 43, "Input"],
Cell[10034, 373, 748, 21, 134, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10819, 399, 312, 10, 43, "Input"],
Cell[11134, 411, 922, 21, 134, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12093, 437, 107, 3, 43, "Input"],
Cell[12203, 442, 3467, 64, 362, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15707, 511, 139, 4, 43, "Input"],
Cell[15849, 517, 295, 6, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16181, 528, 139, 4, 43, "Input"],
Cell[16323, 534, 241, 6, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16601, 545, 168, 5, 43, "Input"],
Cell[16772, 552, 1115, 19, 118, "Output"]
}, Open  ]],
Cell[17902, 574, 70, 0, 41, "Text"],
Cell[17975, 576, 1193, 35, 98, "Input"],
Cell[CellGroupData[{
Cell[19193, 615, 390, 11, 43, "Input"],
Cell[19586, 628, 9607, 163, 378, "Output"]
}, Open  ]],
Cell[29208, 794, 63, 0, 41, "Text"],
Cell[29274, 796, 1097, 34, 71, "Input"],
Cell[CellGroupData[{
Cell[30396, 834, 328, 10, 43, "Input"],
Cell[30727, 846, 14107, 237, 378, "Output"]
}, Open  ]],
Cell[44849, 1086, 1097, 34, 71, "Input"],
Cell[CellGroupData[{
Cell[45971, 1124, 353, 11, 43, "Input"],
Cell[46327, 1137, 13154, 221, 350, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[59518, 1363, 47, 0, 54, "Subsection"],
Cell[59568, 1365, 68, 2, 43, "Input"],
Cell[59639, 1369, 405, 14, 43, "Input"],
Cell[60047, 1385, 164, 4, 43, "Input"],
Cell[60214, 1391, 118, 2, 43, "Input"],
Cell[CellGroupData[{
Cell[60357, 1397, 122, 2, 43, "Input"],
Cell[60482, 1401, 5532, 92, 514, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66051, 1498, 280, 7, 43, "Input"],
Cell[66334, 1507, 524, 10, 68, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66895, 1522, 190, 5, 43, "Input"],
Cell[67088, 1529, 433, 8, 68, "Output"]
}, Open  ]],
Cell[67536, 1540, 260, 9, 43, "Input"],
Cell[CellGroupData[{
Cell[67821, 1553, 823, 23, 98, "Input"],
Cell[CellGroupData[{
Cell[68669, 1580, 298, 7, 32, "Print"],
Cell[68970, 1589, 298, 7, 32, "Print"],
Cell[69271, 1598, 299, 7, 32, "Print"],
Cell[69573, 1607, 299, 7, 32, "Print"],
Cell[69875, 1616, 298, 7, 32, "Print"],
Cell[70176, 1625, 303, 7, 32, "Print"],
Cell[70482, 1634, 303, 7, 32, "Print"],
Cell[70788, 1643, 299, 7, 32, "Print"],
Cell[71090, 1652, 299, 7, 32, "Print"],
Cell[71392, 1661, 299, 7, 32, "Print"],
Cell[71694, 1670, 299, 7, 32, "Print"],
Cell[71996, 1679, 303, 7, 32, "Print"],
Cell[72302, 1688, 299, 7, 32, "Print"],
Cell[72604, 1697, 303, 7, 32, "Print"],
Cell[72910, 1706, 299, 7, 32, "Print"],
Cell[73212, 1715, 299, 7, 32, "Print"],
Cell[73514, 1724, 303, 7, 32, "Print"],
Cell[73820, 1733, 299, 7, 32, "Print"],
Cell[74122, 1742, 299, 7, 32, "Print"],
Cell[74424, 1751, 298, 7, 32, "Print"],
Cell[74725, 1760, 303, 7, 32, "Print"],
Cell[75031, 1769, 299, 7, 32, "Print"],
Cell[75333, 1778, 303, 7, 32, "Print"],
Cell[75639, 1787, 303, 7, 32, "Print"],
Cell[75945, 1796, 299, 7, 32, "Print"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[76305, 1810, 54, 0, 54, "Subsection"],
Cell[76362, 1812, 473, 16, 43, "Input"],
Cell[CellGroupData[{
Cell[76860, 1832, 92, 2, 43, "Input"],
Cell[76955, 1836, 818, 14, 118, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[77810, 1855, 141, 4, 43, "Input"],
Cell[77954, 1861, 3022, 49, 343, "Output"]
}, Open  ]],
Cell[80991, 1913, 64, 0, 41, "Text"],
Cell[CellGroupData[{
Cell[81080, 1917, 402, 13, 43, "Input"],
Cell[81485, 1932, 89, 1, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81611, 1938, 395, 12, 43, "Input"],
Cell[82009, 1952, 152, 4, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82198, 1961, 166, 4, 43, "Input"],
Cell[82367, 1967, 83, 1, 42, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[82499, 1974, 52, 0, 54, "Subsection"],
Cell[82554, 1976, 88, 4, 41, "Text"],
Cell[82645, 1982, 76, 2, 43, "Input"],
Cell[82724, 1986, 51, 1, 43, "Input"],
Cell[82778, 1989, 204, 7, 43, "Input"],
Cell[CellGroupData[{
Cell[83007, 2000, 91, 2, 43, "Input"],
Cell[83101, 2004, 214, 5, 42, "Output"]
}, Open  ]],
Cell[83330, 2012, 243, 8, 41, "Text"],
Cell[83576, 2022, 649, 21, 43, "Input"],
Cell[CellGroupData[{
Cell[84250, 2047, 181, 4, 43, "Input"],
Cell[84434, 2053, 3019, 48, 318, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)