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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     15747,        402]
NotebookOptionsPosition[     15076,        376]
NotebookOutlinePosition[     15421,        391]
CellTagsIndexPosition[     15378,        388]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
Problem 3.17(a) Sakurai.  See that l = 1.\
\>", "Text",
 CellChangeTimes->{{3.658799421873702*^9, 3.658799441795923*^9}, {
  3.6588014510974116`*^9, 3.658801459328019*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"ClearAll", "[", 
   RowBox[{
   "hbar", ",", " ", "theta", ",", " ", "phi", ",", " ", "r", ",", " ", "f", 
    ",", " ", "x", ",", " ", "y", ",", " ", "z", ",", " ", "Lsquared", ",", 
    " ", "Lz", ",", " ", "display", ",", " ", "braket", ",", " ", "psi", ",", 
    " ", "psiS"}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"$Assumptions", "=", 
    RowBox[{"(", 
     RowBox[{"And", "@@", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"#", ">", "0"}], "&"}], ")"}], "/@", 
       RowBox[{"{", 
        RowBox[{"hbar", ",", " ", "r"}], "}"}]}]}], ")"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"x", "[", 
    RowBox[{"r_", ",", " ", "theta_", ",", " ", "phi_"}], "]"}], " ", "=", 
   " ", 
   RowBox[{"r", " ", 
    RowBox[{"Sin", "[", "theta", "]"}], " ", 
    RowBox[{"Cos", "[", "phi", "]"}]}]}], " ", ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"y", "[", 
    RowBox[{"r_", ",", " ", "theta_", ",", " ", "phi_"}], "]"}], " ", "=", 
   " ", 
   RowBox[{"r", " ", 
    RowBox[{"Sin", "[", "theta", "]"}], " ", 
    RowBox[{"Sin", "[", "phi", "]"}]}]}], " ", ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"z", "[", 
     RowBox[{"r_", ",", " ", "theta_", ",", " ", "phi_"}], "]"}], " ", "=", 
    " ", 
    RowBox[{"r", " ", 
     RowBox[{"Cos", "[", "theta", "]"}]}]}], " ", ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"psiS", " ", ":=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"x", "[", 
       RowBox[{"#1", ",", "#2", ",", "#3"}], "]"}], " ", "+", " ", 
      RowBox[{"y", "[", 
       RowBox[{"#1", ",", "#2", ",", "#3"}], "]"}], " ", "+", " ", 
      RowBox[{"3", " ", 
       RowBox[{"z", "[", 
        RowBox[{"#1", ",", "#2", ",", "#3"}], "]"}]}]}], " ", "&"}]}], " ", 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"psi", " ", ":=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"psiS", "[", 
       RowBox[{"#1", ",", "#2", ",", "#3"}], "]"}], 
      RowBox[{"f", "[", "#1", "]"}]}], " ", "&"}]}], " ", ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Lsquared", "[", "f_", "]"}], " ", ":=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"hbar", "^", "2"}]}], " ", 
      RowBox[{"(", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"f", "[", 
            RowBox[{"r", ",", "theta", ",", "phi"}], "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"phi", ",", "2"}], "}"}]}], "]"}], "/", 
         RowBox[{
          RowBox[{"Sin", "[", "theta", "]"}], "^", "2"}]}], " ", "+", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "/", 
           RowBox[{"Sin", "[", "theta", "]"}]}], ")"}], " ", 
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Sin", "[", "theta", "]"}], 
            RowBox[{"D", "[", 
             RowBox[{
              RowBox[{"f", "[", 
               RowBox[{"r", ",", "theta", ",", "phi"}], "]"}], ",", "theta"}],
              "]"}]}], ",", "theta"}], "]"}]}]}], ")"}]}], " ", "//", " ", 
     "Simplify"}]}], " ", ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Lz", "[", "f_", "]"}], " ", ":=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "I"}], " ", "hbar", " ", 
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"f", "[", 
         RowBox[{"r", ",", " ", "theta", ",", " ", "phi"}], "]"}], ",", " ", 
        "phi"}], "]"}]}], "  ", "//", " ", "Simplify"}]}], " ", ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"unit", " ", "tests"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Lsquared", "[", " ", 
     RowBox[{
      RowBox[{"Sin", "[", "#3", "]"}], " ", "&"}], " ", "]"}], 
    "\[IndentingNewLine]", 
    RowBox[{"Lsquared", "[", " ", 
     RowBox[{
      RowBox[{"Cos", "[", "#3", "]"}], " ", "&"}], " ", "]"}], 
    "\[IndentingNewLine]", 
    RowBox[{"Lsquared", "[", " ", 
     RowBox[{
      RowBox[{"Cos", "[", "#2", "]"}], " ", "&"}], " ", "]"}], 
    "\[IndentingNewLine]", 
    RowBox[{"Lsquared", "[", " ", 
     RowBox[{
      RowBox[{"Sin", "[", "#2", "]"}], " ", "&"}], " ", "]"}]}], "*)"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"display", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"hbar", " ", "\[Rule]", " ", "\[HBar]"}], ",", " ", 
     RowBox[{"theta", " ", "\[Rule]", " ", "\[Theta]"}], ",", " ", 
     RowBox[{"phi", " ", "\[Rule]", " ", "\[Phi]"}]}], "}"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"psi", "[", 
    RowBox[{"r", ",", " ", "theta", ",", " ", "phi"}], "]"}], " ", "/.", " ", 
   "display"}], " ", "//", " ", "TraditionalForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Lsquared", "[", " ", "psi", " ", "]"}], "  ", "/.", " ", 
   "display"}], " ", "//", " ", "TraditionalForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Lz", "[", " ", "psi", " ", "]"}], "  ", "/.", " ", "display"}], 
   " ", "//", " ", "TraditionalForm"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.65879855514113*^9, 3.6587988694294415`*^9}, {
   3.6587990041210012`*^9, 3.658799007547123*^9}, {3.6587990571237307`*^9, 
   3.65879906950062*^9}, {3.658799189112285*^9, 3.6587993191292586`*^9}, {
   3.658835962125152*^9, 3.6588360796783648`*^9}, {3.658836144542302*^9, 
   3.658836163453373*^9}, {3.658836247352863*^9, 3.6588363382823734`*^9}, {
   3.6588363760223236`*^9, 3.65883641913488*^9}, {3.6588365138687887`*^9, 
   3.6588365290133677`*^9}, {3.658836566011056*^9, 3.6588366034799075`*^9}, {
   3.6588366454952745`*^9, 3.6588366549227295`*^9}, {3.658836702675954*^9, 
   3.658836705250261*^9}, 3.6588367449763975`*^9, {3.6588369010436296`*^9, 
   3.6588369020635386`*^9}, {3.6590321377937217`*^9, 3.659032139657397*^9}, {
   3.659034250782078*^9, 3.6590342513055573`*^9}, {3.659035969963339*^9, 
   3.6590359909504795`*^9}}],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"f", "(", "r", ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"r", " ", 
      RowBox[{"sin", "(", "\[Theta]", ")"}], " ", 
      RowBox[{"sin", "(", "\[Phi]", ")"}]}], "+", 
     RowBox[{"r", " ", 
      RowBox[{"sin", "(", "\[Theta]", ")"}], " ", 
      RowBox[{"cos", "(", "\[Phi]", ")"}]}], "+", 
     RowBox[{"3", " ", "r", " ", 
      RowBox[{"cos", "(", "\[Theta]", ")"}]}]}], ")"}]}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{
  3.6587988149017854`*^9, {3.6587988495613556`*^9, 3.6587988719597054`*^9}, 
   3.658799008294793*^9, {3.65879931012719*^9, 3.658799319724841*^9}, {
   3.65883601035765*^9, 3.658836039701683*^9}, {3.6588360735808926`*^9, 
   3.6588360805471373`*^9}, 3.6588361644072313`*^9, {3.6588363160334167`*^9, 
   3.6588363466316323`*^9}, 3.658836386024278*^9, 3.6588364202751513`*^9, 
   3.658836542850839*^9, 3.6588366062804203`*^9, 3.658836657041627*^9, 
   3.6588367094190006`*^9, {3.659023395053839*^9, 3.6590234114422946`*^9}, {
   3.6590310324699287`*^9, 3.6590310456636496`*^9}, 3.6590321404625006`*^9, 
   3.659032913972476*^9, 3.659034254584503*^9}],

Cell[BoxData[
 FormBox[
  RowBox[{"2", " ", "r", " ", 
   SuperscriptBox["\[HBar]", "2"], " ", 
   RowBox[{"f", "(", "r", ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"sin", "(", "\[Theta]", ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"sin", "(", "\[Phi]", ")"}], "+", 
        RowBox[{"cos", "(", "\[Phi]", ")"}]}], ")"}]}], "+", 
     RowBox[{"3", " ", 
      RowBox[{"cos", "(", "\[Theta]", ")"}]}]}], ")"}]}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{
  3.6587988149017854`*^9, {3.6587988495613556`*^9, 3.6587988719597054`*^9}, 
   3.658799008294793*^9, {3.65879931012719*^9, 3.658799319724841*^9}, {
   3.65883601035765*^9, 3.658836039701683*^9}, {3.6588360735808926`*^9, 
   3.6588360805471373`*^9}, 3.6588361644072313`*^9, {3.6588363160334167`*^9, 
   3.6588363466316323`*^9}, 3.658836386024278*^9, 3.6588364202751513`*^9, 
   3.658836542850839*^9, 3.6588366062804203`*^9, 3.658836657041627*^9, 
   3.6588367094190006`*^9, {3.659023395053839*^9, 3.6590234114422946`*^9}, {
   3.6590310324699287`*^9, 3.6590310456636496`*^9}, 3.6590321404625006`*^9, 
   3.659032913972476*^9, 3.6590342546755857`*^9}],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"-", "\[ImaginaryI]"}], " ", "r", " ", "\[HBar]", " ", 
   RowBox[{"f", "(", "r", ")"}], " ", 
   RowBox[{"sin", "(", "\[Theta]", ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"cos", "(", "\[Phi]", ")"}], "-", 
     RowBox[{"sin", "(", "\[Phi]", ")"}]}], ")"}]}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{
  3.6587988149017854`*^9, {3.6587988495613556`*^9, 3.6587988719597054`*^9}, 
   3.658799008294793*^9, {3.65879931012719*^9, 3.658799319724841*^9}, {
   3.65883601035765*^9, 3.658836039701683*^9}, {3.6588360735808926`*^9, 
   3.6588360805471373`*^9}, 3.6588361644072313`*^9, {3.6588363160334167`*^9, 
   3.6588363466316323`*^9}, 3.658836386024278*^9, 3.6588364202751513`*^9, 
   3.658836542850839*^9, 3.6588366062804203`*^9, 3.658836657041627*^9, 
   3.6588367094190006`*^9, {3.659023395053839*^9, 3.6590234114422946`*^9}, {
   3.6590310324699287`*^9, 3.6590310456636496`*^9}, 3.6590321404625006`*^9, 
   3.659032913972476*^9, 3.6590342547056127`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"braket", "[", 
    RowBox[{"f1_", ",", " ", "g1_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Integrate", "[", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"Sin", "[", "theta", "]"}], " ", 
      RowBox[{"Conjugate", "[", 
       RowBox[{"f1", "[", 
        RowBox[{"r", ",", " ", "theta", ",", " ", "phi"}], "]"}], "]"}], " ", 
      RowBox[{"g1", "[", 
       RowBox[{"r", ",", " ", "theta", ",", " ", "phi"}], "]"}]}], " ", ",", 
     " ", 
     RowBox[{"{", 
      RowBox[{"theta", ",", " ", "0", ",", " ", "Pi"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"phi", ",", " ", "0", ",", " ", 
       RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"LzPsi", "[", 
    RowBox[{"r_", ",", " ", "theta_", ",", " ", "phi_"}], "]"}], " ", ":=", 
   " ", 
   RowBox[{"Lz", "[", " ", "psiS", " ", "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"braket", "[", " ", 
   RowBox[{"psiS", ",", " ", "LzPsi"}], " ", "]"}], " ", "//", " ", 
  "Simplify"}]}], "Input",
 CellChangeTimes->{{3.6588364492961717`*^9, 3.65883645134501*^9}, {
  3.6588366104411497`*^9, 3.6588366185714397`*^9}, {3.658836747318379*^9, 
  3.658836767934866*^9}, {3.658836875818792*^9, 3.6588369390392017`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "\[Pi]"], 
    RowBox[{
     SubsuperscriptBox["\[Integral]", "0", 
      RowBox[{"2", " ", "\[Pi]"}]], 
     RowBox[{
      RowBox[{
       RowBox[{"Conjugate", "[", 
        RowBox[{"psiS", "[", 
         RowBox[{"r", ",", "theta", ",", "phi"}], "]"}], "]"}], " ", 
       RowBox[{"Sin", "[", "theta", "]"}]}], 
      RowBox[{"\[DifferentialD]", "phi"}], 
      RowBox[{"\[DifferentialD]", "theta"}]}]}]}], ")"}], " ", 
  RowBox[{"Lz", "[", "psiS", "]"}]}]], "Output",
 CellChangeTimes->{
  3.65883643284542*^9, 3.6588364697154803`*^9, 3.658836631526895*^9, 
   3.658836667343031*^9, 3.6588367216674995`*^9, 3.658836776197274*^9, 
   3.658836878887542*^9, {3.658836927565404*^9, 3.6588369404664783`*^9}, 
   3.6590233906915936`*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.659023226592938*^9, 3.659023368010092*^9}, {
   3.6590234003691206`*^9, 3.659023419013999*^9}, {3.6590234776848416`*^9, 
   3.659023640660867*^9}, {3.6590237053799567`*^9, 3.659023761357992*^9}, {
   3.6590239832340755`*^9, 3.6590240717250433`*^9}, {3.6590306910627613`*^9, 
   3.6590309315922036`*^9}, {3.659030973384384*^9, 3.6590309802986126`*^9}, {
   3.6590310134158497`*^9, 3.659031022875285*^9}, {3.659031055373348*^9, 
   3.6590312376597457`*^9}, {3.6590312854443264`*^9, 
   3.6590313934420643`*^9}, {3.6590314893323536`*^9, 
   3.6590316030393095`*^9}, {3.6590316446774144`*^9, 
   3.6590316778585453`*^9}, {3.6590317127675037`*^9, 3.65903172016536*^9}, {
   3.6590317528404455`*^9, 3.6590317976296654`*^9}, {3.6590318984855804`*^9, 
   3.6590319599441595`*^9}, {3.659032052984806*^9, 3.659032113629526*^9}, {
   3.659032179012823*^9, 3.6590321810968876`*^9}, {3.6590322765205684`*^9, 
   3.659032279337122*^9}, {3.6590323249736958`*^9, 3.6590324294806924`*^9}, {
   3.6590324825540843`*^9, 3.6590325116735992`*^9}, {3.659032613424*^9, 
   3.6590326147501397`*^9}, {3.6590326791064835`*^9, 
   3.6590326800413556`*^9}, {3.659032715710073*^9, 3.659032720106032*^9}, {
   3.659032750845251*^9, 3.659032781425968*^9}, {3.6590329617651834`*^9, 
   3.6590329621305065`*^9}, 3.6590333602193766`*^9, {3.6590338434724855`*^9, 
   3.6590339187597885`*^9}, {3.6590339489536633`*^9, 
   3.6590342152248716`*^9}, {3.6590342462229705`*^9, 
   3.6590342896390247`*^9}, {3.6590343480843997`*^9, 3.659034365599945*^9}, 
   3.659034443485089*^9, 3.6590359385123606`*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.6590309887949276`*^9, 3.6590309898038387`*^9}, {
   3.659032736999009*^9, 3.659032739265049*^9}, {3.6590328032678547`*^9, 
   3.6590329065598*^9}, {3.6590329384371624`*^9, 3.6590329571440706`*^9}, {
   3.6590330827398577`*^9, 3.659033161396392*^9}, {3.6590332049301767`*^9, 
   3.6590332687556515`*^9}, {3.6590333030683813`*^9, 
   3.6590333189907303`*^9}, {3.659033400205434*^9, 3.6590335029742002`*^9}, 
   3.6590335403675947`*^9, {3.65903357082638*^9, 3.6590336032022066`*^9}, {
   3.6590336374840927`*^9, 3.659033649879764*^9}, 3.6590336876986036`*^9, {
   3.659033720171196*^9, 3.659033832781952*^9}, 3.6590341263992105`*^9, 
   3.6590343218392725`*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.6590235937873507`*^9, 3.659023596086809*^9}}]
},
WindowSize->{1144, 676},
WindowMargins->{{Automatic, 16}, {Automatic, 22}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 182, 4, 40, "Text"],
Cell[CellGroupData[{
Cell[764, 28, 6350, 168, 816, "Input"],
Cell[7117, 198, 1148, 24, 57, "Output"],
Cell[8268, 224, 1163, 25, 60, "Output"],
Cell[9434, 251, 1022, 20, 57, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10493, 276, 1314, 34, 116, "Input"],
Cell[11810, 312, 828, 21, 95, "Output"]
}, Open  ]],
Cell[12653, 336, 1612, 23, 39, "Input"],
Cell[14268, 361, 707, 10, 39, "Input"],
Cell[14978, 373, 94, 1, 39, "Input"]
}
]
*)

(* End of internal cache information *)