(* 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[ 26684, 709] NotebookOptionsPosition[ 25968, 681] NotebookOutlinePosition[ 26313, 696] CellTagsIndexPosition[ 26270, 693] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"ClearAll", "[", " ", RowBox[{ "g", ",", " ", "ej", ",", " ", "ej", ",", " ", "a", ",", " ", "b", ",", " ", "acap", ",", " ", "bcap", ",", " ", "m", ",", " ", "display", ",", " ", "e"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"g", " ", "=", " ", RowBox[{"1", "/", RowBox[{"(", RowBox[{"1", " ", "-", " ", RowBox[{"rho", "^", "2"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ej", " ", "=", " ", RowBox[{"Exp", "[", RowBox[{"I", " ", "phij"}], " ", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ei", " ", "=", " ", RowBox[{"Exp", "[", RowBox[{"I", " ", "phii"}], " ", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"ej", " ", "-", RowBox[{ RowBox[{"rho", "^", "2"}], "/", "ej"}]}], ")"}], " ", "ei"}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"acap", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "ej"}], ")"}], " ", "-", " ", RowBox[{"ej", " ", RowBox[{"rho", "^", "2"}]}]}], ")"}], "/", "ei"}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", " ", "=", " ", RowBox[{ RowBox[{"-", "rho"}], RowBox[{ RowBox[{"(", RowBox[{"ej", " ", "-", " ", RowBox[{"1", "/", "ej"}]}], ")"}], "/", "ei"}]}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"bcap", " ", "=", " ", RowBox[{"rho", RowBox[{"(", RowBox[{"ej", " ", "-", " ", RowBox[{"1", "/", "ej"}]}], ")"}], " ", "ei"}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"m", " ", "=", " ", RowBox[{"g", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", " ", "b"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"bcap", ",", " ", "acap"}], "}"}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"display", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"phii", " ", "\[Rule]", " ", SubscriptBox["\[Phi]", "i"]}], ",", RowBox[{"phij", " ", "\[Rule]", " ", SubscriptBox["\[Phi]", "j"]}], ",", " ", RowBox[{"rho", " ", "\[Rule]", " ", SubscriptBox["\[Rho]", "ij"]}]}], "}"}]}], " ", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(", RowBox[{"m", " ", "/.", " ", "display"}], ")"}], "//", " ", "MatrixForm"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e", " ", "=", " ", RowBox[{"Eigensystem", "[", "m", "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.6898619928470206`*^9, 3.6898623786830053`*^9}, { 3.6898878920086713`*^9, 3.6898879432068367`*^9}, 3.6898879785726814`*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", RowBox[{"(", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SubscriptBox["\[Phi]", "j"]}]], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}]}], ")"}]}], RowBox[{"1", "-", SubsuperscriptBox["\[Rho]", "ij", "2"]}]], RowBox[{"-", FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SubscriptBox["\[Phi]", "i"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SubscriptBox["\[Phi]", "j"]}]]}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]]}], ")"}], " ", SubscriptBox["\[Rho]", "ij"]}], RowBox[{"1", "-", SubsuperscriptBox["\[Rho]", "ij", "2"]}]]}]}, { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{"-", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SubscriptBox["\[Phi]", "j"]}]]}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]]}], ")"}], " ", SubscriptBox["\[Rho]", "ij"]}], RowBox[{"1", "-", SubsuperscriptBox["\[Rho]", "ij", "2"]}]], FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SubscriptBox["\[Phi]", "i"]}]], " ", RowBox[{"(", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", SubscriptBox["\[Phi]", "j"]}]], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}]}], ")"}]}], RowBox[{"1", "-", SubsuperscriptBox["\[Rho]", "ij", "2"]}]]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "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.6898623049953475`*^9, 3.689862318000568*^9}, { 3.689862358520426*^9, 3.689862380651781*^9}, 3.6898657085857377`*^9, 3.689887922218553*^9}] }, Open ]], Cell[BoxData[{ RowBox[{"ClearAll", "[", " ", RowBox[{ "eigenvalues", ",", " ", "eigenvectors", ",", " ", "e1", ",", " ", "e2"}], " ", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eigenvalues", " ", "=", " ", RowBox[{ RowBox[{"e", "[", RowBox[{"[", "1", "]"}], "]"}], "//", " ", "FullSimplify"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eigenvectors", " ", "=", " ", RowBox[{ RowBox[{"e", "[", RowBox[{"[", "2", "]"}], "]"}], " ", "//", " ", "FullSimplify"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e1", " ", "=", " ", RowBox[{"eigenvectors", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e2", " ", "=", " ", RowBox[{"eigenvectors", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.6898623853566947`*^9, 3.6898624987354846`*^9}, { 3.68986254075554*^9, 3.6898625552289324`*^9}, {3.689862624228874*^9, 3.68986262546376*^9}, {3.6898879473015957`*^9, 3.6898879691956005`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"eigenvalues", " ", "//", " ", RowBox[{"FullSimplify", "/.", " ", "display"}]}], " ", "//", " ", "ColumnForm"}], " ", "//", " ", "TraditionalForm"}]], "Input", CellChangeTimes->{{3.689862576735078*^9, 3.6898625822361293`*^9}, { 3.6898656882594547`*^9, 3.689865691885331*^9}}], Cell[BoxData[ FormBox[ InterpretationBox[GridBox[{ { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]], " ", RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "-", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]]}], ")"}]}], "-", RowBox[{"2", " ", SqrtBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}]], " ", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}], "-", RowBox[{ RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "i"], "-", SubscriptBox["\[Phi]", "j"]}], ")"}], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}]}], ")"}], "2"], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "-", "1"}], ")"}], "2"]}], ")"}]}]]}], "-", "1"}], ")"}]}], RowBox[{"2", " ", RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "-", "1"}], ")"}]}]]}, { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}], "+", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]], " ", RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "-", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]]}], ")"}]}], "+", RowBox[{"2", " ", SqrtBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}]], " ", RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}], "-", RowBox[{ RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "i"], "-", SubscriptBox["\[Phi]", "j"]}], ")"}], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}]}], ")"}], "2"], "-", SuperscriptBox[ RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "-", "1"}], ")"}], "2"]}], ")"}]}]]}], "-", "1"}], ")"}]}], RowBox[{"2", " ", RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "-", "1"}], ")"}]}]]} }, BaselinePosition->{Baseline, {1, 1}}, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}], ColumnForm[{((Rational[1, 2] E^(Complex[0, -1] (Subscript[$CellContext`\[Phi], $CellContext`i] + Subscript[$CellContext`\[Phi], $CellContext`j])))/(-1 + Subscript[$CellContext`\[Rho], $CellContext`ij]^2)) (-1 + E^(Complex[0, 2] Subscript[$CellContext`\[Phi], $CellContext`i]) Subscript[$CellContext`\[Rho], $CellContext`ij]^2 + E^(Complex[0, 2] Subscript[$CellContext`\[Phi], $CellContext`j]) (- E^(Complex[0, 2] Subscript[$CellContext`\[Phi], $CellContext`i]) + Subscript[$CellContext`\[Rho], $CellContext`ij]^2) - 2 (E^(Complex[0, 2] (Subscript[$CellContext`\[Phi], $CellContext`i] + Subscript[$CellContext`\[Phi], $CellContext`j])) (-(-1 + Subscript[$CellContext`\[Rho], $CellContext`ij]^2)^2 + ( Cos[Subscript[$CellContext`\[Phi], $CellContext`i] + Subscript[$CellContext`\[Phi], $CellContext`j]] - Cos[Subscript[$CellContext`\[Phi], $CellContext`i] - Subscript[$CellContext`\[Phi], $CellContext`j]] Subscript[$CellContext`\[Rho], $CellContext`ij]^2)^2))^ Rational[1, 2]), ((Rational[1, 2] E^(Complex[0, -1] (Subscript[$CellContext`\[Phi], $CellContext`i] + Subscript[$CellContext`\[Phi], $CellContext`j])))/(-1 + Subscript[$CellContext`\[Rho], $CellContext`ij]^2)) (-1 + E^(Complex[0, 2] Subscript[$CellContext`\[Phi], $CellContext`i]) Subscript[$CellContext`\[Rho], $CellContext`ij]^2 + E^(Complex[0, 2] Subscript[$CellContext`\[Phi], $CellContext`j]) (- E^(Complex[0, 2] Subscript[$CellContext`\[Phi], $CellContext`i]) + Subscript[$CellContext`\[Rho], $CellContext`ij]^2) + 2 (E^(Complex[0, 2] (Subscript[$CellContext`\[Phi], $CellContext`i] + Subscript[$CellContext`\[Phi], $CellContext`j])) (-(-1 + Subscript[$CellContext`\[Rho], $CellContext`ij]^2)^2 + ( Cos[Subscript[$CellContext`\[Phi], $CellContext`i] + Subscript[$CellContext`\[Phi], $CellContext`j]] - Cos[Subscript[$CellContext`\[Phi], $CellContext`i] - Subscript[$CellContext`\[Phi], $CellContext`j]] Subscript[$CellContext`\[Rho], $CellContext`ij]^2)^2))^ Rational[1, 2])}], Editable->False], TraditionalForm]], "Output", CellChangeTimes->{{3.6898625669186225`*^9, 3.6898625829392805`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"e1", " ", "/.", " ", "display"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"e2", " ", "/.", " ", "display"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}]}]}]], "Input", CellChangeTimes->{{3.6898625972080593`*^9, 3.6898626548002996`*^9}, { 3.6898626945467043`*^9, 3.689862741884673*^9}}], Cell[BoxData[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]], " ", RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]]}], ")"}]}], "+", RowBox[{ SqrtBox["2"], " ", SqrtBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], " ", RowBox[{"cos", "(", RowBox[{"2", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "-", SubscriptBox["\[Phi]", "j"]}], ")"}]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"cos", "(", RowBox[{"2", " ", SubscriptBox["\[Phi]", "i"]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"cos", "(", RowBox[{"2", " ", SubscriptBox["\[Phi]", "j"]}], ")"}]}]}], ")"}]}], "+", RowBox[{"cos", "(", RowBox[{"2", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}], ")"}], "-", SubsuperscriptBox["\[Rho]", "ij", "4"], "+", RowBox[{"4", " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}], "-", "1"}], ")"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}], "-", "1"}], ")"}]}], RowBox[{"2", " ", SubscriptBox["\[Rho]", "ij"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]]}], ")"}]}]]}, {"1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]], TraditionalForm]], "Output", CellChangeTimes->{ 3.689862655956541*^9, {3.6898627088320856`*^9, 3.6898627428534083`*^9}, 3.689866415447497*^9}], Cell[BoxData[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { FractionBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]], " ", RowBox[{"(", RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]]}], ")"}]}], "-", RowBox[{ SqrtBox["2"], " ", SqrtBox[ RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}]], " ", RowBox[{"(", RowBox[{ RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{ SubsuperscriptBox["\[Rho]", "ij", "2"], " ", RowBox[{"cos", "(", RowBox[{"2", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "-", SubscriptBox["\[Phi]", "j"]}], ")"}]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"cos", "(", RowBox[{"2", " ", SubscriptBox["\[Phi]", "i"]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"cos", "(", RowBox[{"2", " ", SubscriptBox["\[Phi]", "j"]}], ")"}]}]}], ")"}]}], "+", RowBox[{"cos", "(", RowBox[{"2", " ", RowBox[{"(", RowBox[{ SubscriptBox["\[Phi]", "i"], "+", SubscriptBox["\[Phi]", "j"]}], ")"}]}], ")"}], "-", SubsuperscriptBox["\[Rho]", "ij", "4"], "+", RowBox[{"4", " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}], "-", "1"}], ")"}]}]]}], "-", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "i"]}]], " ", SubsuperscriptBox["\[Rho]", "ij", "2"]}], "-", "1"}], ")"}]}], RowBox[{"2", " ", SubscriptBox["\[Rho]", "ij"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "\[ImaginaryI]", " ", SubscriptBox["\[Phi]", "j"]}]]}], ")"}]}]]}, {"1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]], TraditionalForm]], "Output", CellChangeTimes->{ 3.689862655956541*^9, {3.6898627088320856`*^9, 3.6898627428534083`*^9}, 3.6898664155412493`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{"ClearAll", "[", "out", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"out", " ", "=", " ", RowBox[{"{", RowBox[{"0", ",", RowBox[{"psi", " ", "ei"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ RowBox[{"m", ".", " ", "out"}], " ", "/.", " ", "display"}], " ", "/.", " ", RowBox[{"psi", " ", "\[Rule]", "\[Psi]"}]}], "//", " ", "Matrixform"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"MatrixPower", "[", RowBox[{"m", ",", "n"}], "]"}], ".", "out"}], " ", "//", "FullSimplify"}], ")"}], " ", "/.", " ", "display"}], " ", "/.", " ", RowBox[{"psi", "\[Rule]", " ", "\[Psi]"}]}], "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}]}], "*)"}]}]}], "Input", CellChangeTimes->{{3.6898627669379125`*^9, 3.6898628415819063`*^9}, { 3.689862871936757*^9, 3.689862938357248*^9}, {3.6898629820458083`*^9, 3.6898629863113155`*^9}, {3.6898631386844544`*^9, 3.6898631793879185`*^9}, {3.6898647006355658`*^9, 3.6898647075582857`*^9}, 3.6898649479448223`*^9, {3.689887989277462*^9, 3.6898880003417625`*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ RowBox[{"e1", ".", "e1"}], " ", "//", " ", RowBox[{"FullSimplify", " ", "/.", " ", RowBox[{"display", "\[IndentingNewLine]", RowBox[{"e1", ".", "e2"}]}]}]}], " ", "//", " ", RowBox[{"FullSimplify", " ", "/.", " ", RowBox[{"display", "\[IndentingNewLine]", RowBox[{"e2", ".", "e2"}]}]}]}], " ", "//", " ", RowBox[{"FullSimplify", " ", "/.", " ", "display"}]}], "\[IndentingNewLine]", "*)"}]], "Input", CellChangeTimes->{{3.6898664204173503`*^9, 3.6898664629006233`*^9}, 3.689866575819744*^9}] }, WindowSize->{1413, 828}, WindowMargins->{{24, Automatic}, {-7, Automatic}}, 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[CellGroupData[{ Cell[579, 22, 2759, 83, 316, "Input"], Cell[3341, 107, 3505, 95, 143, "Output"] }, Open ]], Cell[6861, 205, 1071, 27, 166, "Input"], Cell[CellGroupData[{ Cell[7957, 236, 333, 7, 39, "Input"], Cell[8293, 245, 7344, 166, 199, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15674, 416, 516, 13, 91, "Input"], Cell[16193, 431, 3887, 96, 133, "Output"], Cell[20083, 529, 3889, 96, 133, "Output"] }, Open ]], Cell[23987, 628, 1361, 34, 116, "Input"], Cell[25351, 664, 613, 15, 116, "Input"] } ] *) (* End of internal cache information *)