(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[      5985,        180]
NotebookOptionsPosition[      5124,        153]
NotebookOutlinePosition[      5498,        169]
CellTagsIndexPosition[      5455,        166]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Time evolution of a non-stationary state", "Title",
 CellChangeTimes->{{3.716648179399791*^9, 
  3.716648189650872*^9}},ExpressionUUID->"a194c54d-d92b-40be-abbe-\
e55d4a1077a2"],

Cell["Infinite square well", "Text",
 CellChangeTimes->{{3.716648261283431*^9, 
  3.7166482665964913`*^9}},ExpressionUUID->"ce3e4edc-6f44-42ea-94c5-\
47c8b17663db"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Psi]", "[", 
   RowBox[{"x_", ",", "n_", ",", "a_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"Sqrt", "[", 
    RowBox[{"2", "/", "a"}], "]"}], 
   RowBox[{"Sin", "[", 
    RowBox[{"n", " ", "Pi", " ", 
     RowBox[{"x", "/", "a"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.7161170481948166`*^9, 
  3.7161170864958096`*^9}},ExpressionUUID->"73247eef-a971-407a-91ad-\
769f0f4db16f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Psi]Q", "[", 
   RowBox[{
   "x_", ",", "t_", ",", "c1_", ",", "\[Phi]_", ",", "n1_", ",", "n2_"}], 
   "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"\[Psi]1", ",", "\[Psi]2", ",", "c2"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"\[Psi]1", "=", 
      RowBox[{"\[Psi]", "[", 
       RowBox[{"x", ",", "n1", ",", "10"}], "]"}]}], ";", 
     RowBox[{"\[Psi]2", "=", 
      RowBox[{"\[Psi]", "[", 
       RowBox[{"x", ",", "n2", ",", "10"}], "]"}]}], ";", 
     RowBox[{"c2", "=", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"1", "-", 
        RowBox[{"c1", "^", "2"}]}], "]"}]}], ";", 
     RowBox[{
      RowBox[{
       RowBox[{"c1", "^", "2"}], " ", 
       RowBox[{"\[Psi]1", "^", "2"}]}], " ", "+", 
      RowBox[{
       RowBox[{"c2", "^", "2"}], " ", 
       RowBox[{"\[Psi]2", "^", "2"}]}], " ", "+", 
      RowBox[{"2", "c1", " ", "c2", " ", "\[Psi]1", " ", "\[Psi]2", " ", 
       RowBox[{"Cos", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"n1", "^", "2"}], "-", 
            RowBox[{"n2", "^", "2"}]}], ")"}], "t"}], "+", "\[Phi]"}], 
        "]"}]}]}]}]}], " ", "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.7166471070025063`*^9, 3.7166471471004877`*^9}, {
  3.7166471797445574`*^9, 3.7166472125136395`*^9}, {3.716647245642135*^9, 
  3.716647395001645*^9}, {3.7166475476583786`*^9, 
  3.7166475593627386`*^9}},ExpressionUUID->"783b2e8c-0fb6-4f6d-8a2a-\
42452185e191"],

Cell["Scaled for a nucleon in a nucleus", "Text",
 CellChangeTimes->{{3.716648234608741*^9, 
  3.7166482489696317`*^9}},ExpressionUUID->"8b23cb17-34e3-457d-b763-\
80a2a7c3a9e3"],

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"\[Psi]Q", "[", 
      RowBox[{"x", ",", "t", ",", "c", ",", "\[Phi]", ",", "n1", ",", "n2"}], 
      "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "10"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0.4"}], "}"}]}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<x(fm)\>\"", ",", "None"}], "}"}]}], ",", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"FontSize", "\[Rule]", "14"}]}], ",", 
     RowBox[{"PerformanceGoal", "\[Rule]", "\"\<Quality\>\""}], ",", 
     RowBox[{"Epilog", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Text", "[", 
        RowBox[{"\"\<time*\!\(\*SuperscriptBox[\(10\), \(22\)]\) s\>\"", ",", 
         RowBox[{"{", 
          RowBox[{"5", ",", "0.37"}], "}"}]}], "]"}], "}"}]}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"t", ",", "0"}], "}"}], ",", "0", ",", "5"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"c", ",", "0.707"}], "}"}], ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"\[Phi]", ",", "0"}], "}"}], ",", "0", ",", 
     RowBox[{"2", "Pi"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"n1", ",", "1"}], "}"}], ",", "1", ",", "20", ",", "1"}], "}"}],
    ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"n2", ",", "2"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
      "1", ",", "2", ",", "3", ",", "4", ",", "10", ",", "20", ",", "50"}], 
      "}"}]}], " ", "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7166477606524806`*^9, 3.716647858381542*^9}, {
  3.716647906074075*^9, 3.7166479242478676`*^9}, {3.7166480124603004`*^9, 
  3.7166480244303136`*^9}, {3.7166483546526575`*^9, 
  3.7166483751548214`*^9}},ExpressionUUID->"e6c563ed-4a85-4dd9-a7f7-\
c44f06d66907"]
}, Open  ]]
},
WindowSize->{759, 833},
WindowMargins->{{Automatic, 572}, {64, Automatic}},
Magnification:>0.9 Inherited,
FrontEndVersion->"11.1 for Microsoft Windows (64-bit) (April 18, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 183, 3, 91, "Title", "ExpressionUUID" -> \
"a194c54d-d92b-40be-abbe-e55d4a1077a2"],
Cell[766, 27, 164, 3, 28, "Text", "ExpressionUUID" -> \
"ce3e4edc-6f44-42ea-94c5-47c8b17663db"],
Cell[933, 32, 418, 12, 28, "Input", "ExpressionUUID" -> \
"73247eef-a971-407a-91ad-769f0f4db16f"],
Cell[1354, 46, 1558, 42, 119, "Input", "ExpressionUUID" -> \
"783b2e8c-0fb6-4f6d-8a2a-42452185e191"],
Cell[2915, 90, 177, 3, 28, "Text", "ExpressionUUID" -> \
"8b23cb17-34e3-457d-b763-80a2a7c3a9e3"],
Cell[3095, 95, 2013, 55, 97, "Input", "ExpressionUUID" -> \
"e6c563ed-4a85-4dd9-a7f7-c44f06d66907"]
}, Open  ]]
}
]
*)