(************** Content-type: application/mathematica **************
                     CreatedBy='Mathematica 5.1'

                    Mathematica-Compatible Notebook

This notebook can be used with any Mathematica-compatible
application, such as Mathematica, MathReader or Publicon. The data
for the notebook starts with the line containing stars above.

To get the notebook into a Mathematica-compatible application, do
one of the following:

* Save the data starting with the line of stars above into a file
  with a name ending in .nb, then open the file inside the
  application;

* Copy the data starting with the line of stars above to the
  clipboard, then use the Paste menu command inside the application.

Data for notebooks contains only printable 7-bit ASCII and can be
sent directly in email or through ftp in text mode.  Newlines can be
CR, LF or CRLF (Unix, Macintosh or MS-DOS style).

NOTE: If you modify the data for this notebook not in a Mathematica-
compatible application, you must delete the line below containing
the word CacheID, otherwise Mathematica-compatible applications may
try to use invalid cache data.

For more information on notebooks and Mathematica-compatible 
applications, contact Wolfram Research:
  web: http://www.wolfram.com
  email: info@wolfram.com
  phone: +1-217-398-0700 (U.S.)

Notebook reader applications are available free of charge from 
Wolfram Research.
*******************************************************************)

(*CacheID: 232*)


(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[     15525,        489]*)
(*NotebookOutlinePosition[     16215,        513]*)
(*  CellTagsIndexPosition[     16171,        509]*)
(*WindowFrame->Normal*)



Notebook[{

Cell[CellGroupData[{
Cell["Problem Set 4", "Subtitle"],

Cell[CellGroupData[{

Cell["Prob. 2", "Subsection"],

Cell["\<\
First, express the radial wave functions P and Q in terms of F1 and \
F2\
\>", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(x\  = \ 2\ \[Lambda]\ r\)], "Input",
  CellLabel->"In[1]:="],

Cell[BoxData[
    \(2\ r\ \[Lambda]\)], "Output",
  CellLabel->"Out[1]="]
}, Open  ]],

Cell[BoxData[
    \(P[x_]\  := \ 
      Sqrt[1 + \[Epsilon]]\ Exp[\(-x\)/2] \((F1[x] + F2[x])\)\)], "Input",
  CellLabel->"In[2]:="],

Cell[BoxData[
    \(Q[x_]\  := \ 
      Sqrt[1 - \[Epsilon]]\ Exp[\(-x\)/2] \((F1[x] - F2[x])\)\)], "Input",
  CellLabel->"In[3]:="],

Cell["Next,  write down the radial Dirac equations", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(eq1\  = \ \((\(-Z\)/r + c^2)\)\ P[x]\  + \ 
          c \((D[Q[x], r] - \[Kappa]/r\ Q[x])\)\  \[Equal] 
        c^2\ \[Epsilon]\ P[x]\)], "Input",
  CellLabel->"In[4]:="],

Cell[BoxData[
    RowBox[{
      RowBox[{\(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \((c\^2 - 
              Z\/r)\)\ \@\(1 + \[Epsilon]\)\ \((F1[2\ r\ \[Lambda]] + 
              F2[2\ r\ \[Lambda]])\)\), "+", 
        RowBox[{"c", " ", 
          RowBox[{"(", 
            
            RowBox[{\(-\(\(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \@\(1 - \
\[Epsilon]\)\ \[Kappa]\ \((F1[2\ r\ \[Lambda]] - 
                          F2[2\ r\ \[Lambda]])\)\)\/r\)\), 
              "-", \(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \@\(1 - \
\[Epsilon]\)\ \[Lambda]\ \((F1[2\ r\ \[Lambda]] - F2[2\ r\ \[Lambda]])\)\), 
              "+", 
              
              RowBox[{\(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\), 
                " ", \(\@\(1 - \[Epsilon]\)\), " ", 
                RowBox[{"(", 
                  RowBox[{
                    RowBox[{"2", " ", "\[Lambda]", " ", 
                      RowBox[{
                        SuperscriptBox["F1", "\[Prime]",
                          MultilineFunction->None], "[", \(2\ r\ \[Lambda]\), 
                        "]"}]}], "-", 
                    RowBox[{"2", " ", "\[Lambda]", " ", 
                      RowBox[{
                        SuperscriptBox["F2", "\[Prime]",
                          MultilineFunction->None], "[", \(2\ r\ \[Lambda]\), 
                        "]"}]}]}], ")"}]}]}], ")"}]}]}], 
      "\[Equal]", \(c\^2\ \[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \[Epsilon]\ \
\@\(1 + \[Epsilon]\)\ \((F1[2\ r\ \[Lambda]] + 
            F2[2\ r\ \[Lambda]])\)\)}]], "Output",
  CellLabel->"Out[4]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(eq2\  = \ \((\(-Z\)/r - c^2)\)\ Q[x]\  - \ 
          c \((D[P[x], r] + \[Kappa]/r\ P[x])\)\  \[Equal] 
        c^2\ \[Epsilon]\ Q[x]\)], "Input",
  CellLabel->"In[5]:="],

Cell[BoxData[
    RowBox[{
      RowBox[{\(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \((\(-c\^2\) - 
              Z\/r)\)\ \@\(1 - \[Epsilon]\)\ \((F1[2\ r\ \[Lambda]] - 
              F2[2\ r\ \[Lambda]])\)\), "-", 
        RowBox[{"c", " ", 
          RowBox[{"(", 
            
            RowBox[{\(\(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \@\(1 + \
\[Epsilon]\)\ \[Kappa]\ \((F1[2\ r\ \[Lambda]] + 
                        F2[2\ r\ \[Lambda]])\)\)\/r\), 
              "-", \(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \@\(1 + \
\[Epsilon]\)\ \[Lambda]\ \((F1[2\ r\ \[Lambda]] + F2[2\ r\ \[Lambda]])\)\), 
              "+", 
              
              RowBox[{\(\[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\), 
                " ", \(\@\(1 + \[Epsilon]\)\), " ", 
                RowBox[{"(", 
                  RowBox[{
                    RowBox[{"2", " ", "\[Lambda]", " ", 
                      RowBox[{
                        SuperscriptBox["F1", "\[Prime]",
                          MultilineFunction->None], "[", \(2\ r\ \[Lambda]\), 
                        "]"}]}], "+", 
                    RowBox[{"2", " ", "\[Lambda]", " ", 
                      RowBox[{
                        SuperscriptBox["F2", "\[Prime]",
                          MultilineFunction->None], "[", \(2\ r\ \[Lambda]\), 
                        "]"}]}]}], ")"}]}]}], ")"}]}]}], 
      "\[Equal]", \(c\^2\ \[ExponentialE]\^\(\(-r\)\ \[Lambda]\)\ \@\(1 - \
\[Epsilon]\)\ \[Epsilon]\ \((F1[2\ r\ \[Lambda]] - 
            F2[2\ r\ \[Lambda]])\)\)}]], "Output",
  CellLabel->"Out[5]="]
}, Open  ]],

Cell["Now solve the equations for dF1/dx and dF2/dx", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Solve[{eq1, 
        eq2}, {\(F1'\)[2\ \[Lambda]\ r], \(F2'\)[
          2\ \[Lambda]\ r]}]\)], "Input",
  CellLabel->"In[6]:="],

Cell[BoxData[
    RowBox[{"{", 
      RowBox[{"{", 
        RowBox[{
          RowBox[{
            RowBox[{
              SuperscriptBox["F1", "\[Prime]",
                MultilineFunction->None], "[", \(2\ r\ \[Lambda]\), "]"}], 
            "\[Rule]", \(\(-\(\(1\/\(2\ c\ r\ \@\(1 - \[Epsilon]\)\ \[Lambda]\
\)\)\((c\^2\ r\ \@\(1 + \[Epsilon]\)\ F1[2\ r\ \[Lambda]] - 
                    Z\ \@\(1 + \[Epsilon]\)\ F1[2\ r\ \[Lambda]] - 
                    c\^2\ r\ \[Epsilon]\ \@\(1 + \[Epsilon]\)\ F1[
                        2\ r\ \[Lambda]] - 
                    c\ \@\(1 - \[Epsilon]\)\ \[Kappa]\ F1[2\ r\ \[Lambda]] - 
                    c\ r\ \@\(1 - \[Epsilon]\)\ \[Lambda]\ F1[
                        2\ r\ \[Lambda]] + 
                    c\^2\ r\ \@\(1 + \[Epsilon]\)\ F2[2\ r\ \[Lambda]] - 
                    Z\ \@\(1 + \[Epsilon]\)\ F2[2\ r\ \[Lambda]] - 
                    c\^2\ r\ \[Epsilon]\ \@\(1 + \[Epsilon]\)\ F2[
                        2\ r\ \[Lambda]] + 
                    c\ \@\(1 - \[Epsilon]\)\ \[Kappa]\ F2[2\ r\ \[Lambda]] + 
                    c\ r\ \@\(1 - \[Epsilon]\)\ \[Lambda]\ F2[
                        2\ r\ \[Lambda]])\)\)\) - \((Z\ F1[2\ r\ \[Lambda]] + 
                    c\ \@\(1 - \[Epsilon]\)\ \@\(1 + \[Epsilon]\)\ \[Kappa]\ \
F1[2\ r\ \[Lambda]] - c\^2\ r\ F2[2\ r\ \[Lambda]] + 
                    Z\ \[Epsilon]\ F2[2\ r\ \[Lambda]] + 
                    c\^2\ r\ \[Epsilon]\^2\ F2[2\ r\ \[Lambda]] - 
                    c\ r\ \@\(1 - \[Epsilon]\)\ \@\(1 + \[Epsilon]\)\ \
\[Lambda]\ F2[
                        2\ r\ \[Lambda]])\)/\((2\ c\ r\ \@\(1 - \[Epsilon]\)\ \
\@\(1 + \[Epsilon]\)\ \[Lambda])\)\)}], ",", 
          RowBox[{
            RowBox[{
              SuperscriptBox["F2", "\[Prime]",
                MultilineFunction->None], "[", \(2\ r\ \[Lambda]\), "]"}], 
            "\[Rule]", \(-\(\((Z\ F1[2\ r\ \[Lambda]] + 
                    c\ \@\(1 - \[Epsilon]\)\ \@\(1 + \[Epsilon]\)\ \[Kappa]\ \
F1[2\ r\ \[Lambda]] - c\^2\ r\ F2[2\ r\ \[Lambda]] + 
                    Z\ \[Epsilon]\ F2[2\ r\ \[Lambda]] + 
                    c\^2\ r\ \[Epsilon]\^2\ F2[2\ r\ \[Lambda]] - 
                    c\ r\ \@\(1 - \[Epsilon]\)\ \@\(1 + \[Epsilon]\)\ \
\[Lambda]\ F2[
                        2\ r\ \[Lambda]])\)/\((2\ c\ r\ \@\(1 - \[Epsilon]\)\ \
\@\(1 + \[Epsilon]\)\ \[Lambda])\)\)\)}]}], "}"}], "}"}]], "Output",
  CellLabel->"Out[6]="]
}, Open  ]],

Cell["Simplify the result;  dF1/dx -> f1p  dF2/dx -> f2p", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(f1p\  = \ 
      Expand[FullSimplify[\(F1'\)[
              2\ \[Lambda]\ r]\  /. \ %6]]\)\(\[IndentingNewLine]\)
    \)\)], "Input",
  CellLabel->"In[7]:="],

Cell[BoxData[
    \({1\/2\ F1[
            2\ r\ \[Lambda]] - \(c\ F1[2\ r\ \[Lambda]]\)\/\(2\ \@\(1 - \
\[Epsilon]\^2\)\ \[Lambda]\) + \(Z\ \[Epsilon]\ F1[2\ r\ \[Lambda]]\)\/\(2\ c\
\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) + \(c\ \[Epsilon]\^2\ F1[2\ r\ \
\[Lambda]]\)\/\(2\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) + \(Z\ F2[2\ r\ \
\[Lambda]]\)\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) - \(\[Kappa]\ \
F2[2\ r\ \[Lambda]]\)\/\(2\ r\ \[Lambda]\)}\)], "Output",
  CellLabel->"Out[7]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(f2p\  = \ 
      Expand[FullSimplify[\(F2'\)[2\ \[Lambda]\ r]\  /. \ %6]]\)], "Input",
  CellLabel->"In[8]:="],

Cell[BoxData[
    \({\(-\(\(Z\ F1[
                  2\ r\ \[Lambda]]\)\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \
\[Lambda]\)\)\) - \(\[Kappa]\ F1[2\ r\ \[Lambda]]\)\/\(2\ r\ \[Lambda]\) + 
        1\/2\ F2[
            2\ r\ \[Lambda]] + \(c\ F2[2\ r\ \[Lambda]]\)\/\(2\ \@\(1 - \
\[Epsilon]\^2\)\ \[Lambda]\) - \(Z\ \[Epsilon]\ F2[2\ r\ \[Lambda]]\)\/\(2\ c\
\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) - \(c\ \[Epsilon]\^2\ F2[2\ r\ \
\[Lambda]]\)\/\(2\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\)}\)], "Output",
  CellLabel->"Out[8]="]
}, Open  ]],

Cell["Find coefficients of F1 and F2 on RHS of equations", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(f11\  = \ Coefficient[f1p, \ F1[2\ \[Lambda]\ r]]\)], "Input",
  CellLabel->"In[9]:="],

Cell[BoxData[
    \({1\/2 - 
        c\/\(2\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) + \(Z\ \[Epsilon]\)\/\(2\
\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) + \(c\ \[Epsilon]\^2\)\/\(2\ \
\@\(1 - \[Epsilon]\^2\)\ \[Lambda]\)}\)], "Output",
  CellLabel->"Out[9]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(f12\  = \ \ \ Coefficient[f1p, \ 
        F2[2\ \[Lambda]\ r]]\)\(\[IndentingNewLine]\)
    \)\)], "Input",
  CellLabel->"In[10]:="],

Cell[BoxData[
    \({Z\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) - \[Kappa]\/\(2\ r\ \
\[Lambda]\)}\)], "Output",
  CellLabel->"Out[10]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(f21\  = \ Coefficient[f2p, \ F1[2\ \[Lambda]\ r]]\)], "Input",
  CellLabel->"In[11]:="],

Cell[BoxData[
    \({\(-\(Z\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\)\)\) - \[Kappa]\
\/\(2\ r\ \[Lambda]\)}\)], "Output",
  CellLabel->"Out[11]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(f22\  = \ Coefficient[f2p, \ F2[2\ \[Lambda]\ r]]\)], "Input",
  CellLabel->"In[12]:="],

Cell[BoxData[
    \({1\/2 + 
        c\/\(2\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) - \(Z\ \[Epsilon]\)\/\(2\
\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) - \(c\ \[Epsilon]\^2\)\/\(2\ \
\@\(1 - \[Epsilon]\^2\)\ \[Lambda]\)}\)], "Output",
  CellLabel->"Out[12]="]
}, Open  ]],

Cell["\<\
Combine 1st, 2nd, and 4th terms of f11 and of f22 to further \
simplify\
\>", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Simplify[
      f11[\([1, 1]\)] + f11[\([1, 2]\)] + \ 
          f11[\([1, 4]\)] /. \ \[Lambda]\  -> \ 
          c\ Sqrt[1 - \[Epsilon]^2]]\)], "Input",
  CellLabel->"In[13]:="],

Cell[BoxData[
    \(0\)], "Output",
  CellLabel->"Out[13]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(Simplify[
      f22[\([1, 1]\)] + f22[\([1, 2]\)] + \ 
          f22[\([1, 4]\)] /. \ \[Lambda]\  -> \ 
          c\ Sqrt[1 - \[Epsilon]^2]]\)\(\[IndentingNewLine]\)
    \)\)], "Input",
  CellLabel->"In[14]:="],

Cell[BoxData[
    \(1\)], "Output",
  CellLabel->"Out[14]="]
}, Open  ]],

Cell["\<\
Recombine the simplified coefficients and write out the resulting \
RHS terms\
\>", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(g11\  = \ \ f11[\([1, 3]\)]\ F1[x]\)], "Input",
  CellLabel->"In[15]:="],

Cell[BoxData[
    \(\(Z\ \[Epsilon]\ F1[2\ r\ \[Lambda]]\)\/\(2\ c\ r\ \@\(1 - \
\[Epsilon]\^2\)\ \[Lambda]\)\)], "Output",
  CellLabel->"Out[15]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(g12\  = f12[\([1]\)]\ F2[x]\)], "Input",
  CellLabel->"In[16]:="],

Cell[BoxData[
    \(\((Z\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\) - \[Kappa]\/\(2\ \
r\ \[Lambda]\))\)\ F2[2\ r\ \[Lambda]]\)], "Output",
  CellLabel->"Out[16]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(g21\  = \ f21[\([1]\)]\ F1[x]\)], "Input",
  CellLabel->"In[17]:="],

Cell[BoxData[
    \(\((\(-\(Z\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\)\)\) - \
\[Kappa]\/\(2\ r\ \[Lambda]\))\)\ F1[2\ r\ \[Lambda]]\)], "Output",
  CellLabel->"Out[17]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(g22\  = \ \((1 + \ f22[\([1, 3]\)])\)\ F2[x]\)], "Input",
  CellLabel->"In[18]:="],

Cell[BoxData[
    \(\((1 - \(Z\ \[Epsilon]\)\/\(2\ c\ r\ \@\(1 - \[Epsilon]\^2\)\ \[Lambda]\
\))\)\ F2[2\ r\ \[Lambda]]\)], "Output",
  CellLabel->"Out[18]="]
}, Open  ]],

Cell[BoxData[
    \(Clear[x]\)], "Input",
  CellLabel->"In[19]:="],

Cell[CellGroupData[{

Cell[BoxData[
    \(r\  = \ x/\((2\ \[Lambda])\)\)], "Input",
  CellLabel->"In[20]:="],

Cell[BoxData[
    \(x\/\(2\ \[Lambda]\)\)], "Output",
  CellLabel->"Out[20]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(g11 + g12\)], "Input",
  CellLabel->"In[21]:="],

Cell[BoxData[
    \(\(Z\ \[Epsilon]\ F1[x]\)\/\(c\ x\ \@\(1 - \[Epsilon]\^2\)\) + \((Z\/\(c\
\ x\ \@\(1 - \[Epsilon]\^2\)\) - \[Kappa]\/x)\)\ F2[x]\)], "Output",
  CellLabel->"Out[21]="]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(g21 + g22\)], "Input",
  CellLabel->"In[22]:="],

Cell[BoxData[
    \(\((\(-\(Z\/\(c\ x\ \@\(1 - \[Epsilon]\^2\)\)\)\) - \[Kappa]\/x)\)\ F1[
          x] + \((1 - \(Z\ \[Epsilon]\)\/\(c\ x\ \@\(1 - \[Epsilon]\^2\)\))\)\
\ F2[x]\)], "Output",
  CellLabel->"Out[22]="]
}, Open  ]],

Cell[TextData[{
  "Therefore, noting that   \[Lambda] = c ",
  Cell[BoxData[
      \(TraditionalForm\`\@\(1 - \[Epsilon]\^2\)\)]],
  ", we find\n\n",
  Cell[BoxData[
      \(TraditionalForm\`dF1\/dx\)]],
  "=",
  StyleBox[" g11 + g12 ",
    FontFamily->"Courier"],
  "=  ",
  Cell[BoxData[
      \(TraditionalForm\`Z\[Epsilon]\/\(\[Lambda]\ x\)\)]],
  " ",
  StyleBox["F1[x]",
    FontFamily->"Courier"],
  " +",
  Cell[BoxData[
      \(\((Z\/\(\(\ \)\(\[Lambda]\ x\)\) - \[Kappa]\/x)\)\ F2[x]\)]],
  "\n",
  Cell[BoxData[
      \(TraditionalForm\`dF2\/dx\)]],
  "=  ",
  StyleBox["g21 + g22",
    FontFamily->"Courier"],
  "  = ",
  Cell[BoxData[
      \(\((\(-\(Z\/\(\(\ \)\(\[Lambda]x\)\)\)\) - \[Kappa]\/x)\)\ F1[x]\)]],
  "+ ",
  Cell[BoxData[
      \(\((1 - \(Z\ \[Epsilon]\)\/\(\(\ \)\(\[Lambda]\ x\)\))\)\ F2[x]\)]]
}], "Text"]
}, Closed]]
}, Open  ]]
},
FrontEndVersion->"5.1 for X",
ScreenRectangle->{{0, 1280}, {0, 1024}},
WindowSize->{788, 987},
WindowMargins->{{120, Automatic}, {0, Automatic}},
Magnification->1.5,
StyleDefinitions -> "ArticleModern.nb"
]

(*******************************************************************
Cached data follows.  If you edit this Notebook file directly, not
using Mathematica, you must remove the line containing CacheID at
the top of  the file.  The cache data will then be recreated when
you save this file from within Mathematica.
*******************************************************************)

(*CellTagsOutline
CellTagsIndex->{}
*)

(*CellTagsIndex
CellTagsIndex->{}
*)

(*NotebookFileOutline
Notebook[{

Cell[CellGroupData[{
Cell[1776, 53, 33, 0, 59, "Subtitle"],

Cell[CellGroupData[{
Cell[1834, 57, 29, 0, 82, "Subsection"],
Cell[1866, 59, 96, 3, 37, "Text"],

Cell[CellGroupData[{
Cell[1987, 66, 80, 2, 61, "Input"],
Cell[2070, 70, 73, 2, 60, "Output"]
}, Open  ]],
Cell[2158, 75, 132, 3, 61, "Input"],
Cell[2293, 80, 132, 3, 61, "Input"],
Cell[2428, 85, 60, 0, 37, "Text"],

Cell[CellGroupData[{
Cell[2513, 89, 190, 4, 84, "Input"],
Cell[2706, 95, 1569, 32, 267, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[4312, 132, 190, 4, 84, "Input"],
Cell[4505, 138, 1567, 32, 267, "Output"]
}, Open  ]],
Cell[6087, 173, 61, 0, 37, "Text"],

Cell[CellGroupData[{
Cell[6173, 177, 148, 4, 61, "Input"],
Cell[6324, 183, 2417, 44, 415, "Output"]
}, Open  ]],
Cell[8756, 230, 66, 0, 37, "Text"],

Cell[CellGroupData[{
Cell[8847, 234, 181, 5, 84, "Input"],
Cell[9031, 241, 492, 8, 149, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[9560, 254, 130, 3, 61, "Input"],
Cell[9693, 259, 527, 9, 149, "Output"]
}, Open  ]],
Cell[10235, 271, 66, 0, 37, "Text"],

Cell[CellGroupData[{
Cell[10326, 275, 106, 2, 61, "Input"],
Cell[10435, 279, 259, 5, 93, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[10731, 289, 156, 4, 84, "Input"],
Cell[10890, 295, 145, 3, 87, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[11072, 303, 107, 2, 61, "Input"],
Cell[11182, 307, 154, 3, 87, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[11373, 315, 107, 2, 61, "Input"],
Cell[11483, 319, 260, 5, 93, "Output"]
}, Open  ]],
Cell[11758, 327, 95, 3, 37, "Text"],

Cell[CellGroupData[{
Cell[11878, 334, 198, 5, 84, "Input"],
Cell[12079, 341, 60, 2, 60, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[12176, 348, 234, 6, 107, "Input"],
Cell[12413, 356, 60, 2, 60, "Output"]
}, Open  ]],
Cell[12488, 361, 101, 3, 37, "Text"],

Cell[CellGroupData[{
Cell[12614, 368, 92, 2, 61, "Input"],
Cell[12709, 372, 148, 3, 88, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[12894, 380, 85, 2, 61, "Input"],
Cell[12982, 384, 170, 3, 90, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[13189, 392, 87, 2, 61, "Input"],
Cell[13279, 396, 179, 3, 90, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[13495, 404, 102, 2, 61, "Input"],
Cell[13600, 408, 158, 3, 90, "Output"]
}, Open  ]],
Cell[13773, 414, 66, 2, 61, "Input"],

Cell[CellGroupData[{
Cell[13864, 420, 86, 2, 61, "Input"],
Cell[13953, 424, 78, 2, 78, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[14068, 431, 67, 2, 61, "Input"],
Cell[14138, 435, 186, 3, 90, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[14361, 443, 67, 2, 61, "Input"],
Cell[14431, 447, 216, 4, 90, "Output"]
}, Open  ]],
Cell[14662, 454, 835, 31, 131, "Text"]
}, Closed]]
}, Open  ]]
}
]
*)



(*******************************************************************
End of Mathematica Notebook file.
*******************************************************************)