(* 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[ 20184, 570] NotebookOptionsPosition[ 19781, 552] NotebookOutlinePosition[ 20124, 567] CellTagsIndexPosition[ 20081, 564] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"ClearAll", "[", RowBox[{"theta", ",", " ", "phi"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"theta_", ",", " ", "phi_"}], "]"}], " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Cos", "[", "phi", "]"}], RowBox[{"Cos", "[", "theta", "]"}]}], ",", " ", RowBox[{ RowBox[{"Sin", "[", "phi", "]"}], RowBox[{"Cos", "[", "theta", "]"}]}], ",", " ", RowBox[{"-", RowBox[{"Sin", "[", "theta", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", "phi", "]"}]}], " ", RowBox[{"Sin", "[", "theta", "]"}]}], ",", " ", RowBox[{ RowBox[{"Cos", "[", "phi", "]"}], RowBox[{"Sin", "[", "theta", "]"}]}], ",", " ", "0"}], "}"}]}], "}"}]}], " ", ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"dtheta", "[", RowBox[{"theta_", ",", " ", "phi_"}], "]"}], " ", ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"theta", ",", " ", "phi"}], "]"}], ",", "theta"}], "]"}]}], " ", RowBox[{"(*", RowBox[{"/.", " ", RowBox[{"{", RowBox[{ RowBox[{"theta", " ", "\[Rule]", " ", "t"}], ",", " ", RowBox[{"phi", " ", "\[Rule]", " ", "p"}]}], "}"}]}], "*)"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"dphi", "[", RowBox[{"theta_", ",", " ", "phi_"}], "]"}], " ", ":=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"theta", ",", " ", "phi"}], "]"}], ",", "phi"}], "]"}]}], " ", RowBox[{"(*", RowBox[{"/.", " ", RowBox[{"{", RowBox[{ RowBox[{"theta", " ", "\[Rule]", " ", "t"}], ",", " ", RowBox[{"phi", " ", "\[Rule]", " ", "p"}]}], "}"}]}], "*)"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Grid", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ "\"\<\!\(\*SubscriptBox[\(A\), \(k\)]\) = \>\"", ",", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "\"\<\[PartialD]\!\(\*SubscriptBox[\(A\), \ \(k\)]\)/\[PartialD]\!\(\*SubscriptBox[\(\[Theta]\), \(k\)]\) = \>\"", ",", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"dtheta", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "\"\<\[PartialD]\!\(\*SubscriptBox[\(A\), \ \(k\)]\)/\[PartialD]\!\(\*SubscriptBox[\(\[Phi]\), \(k\)]\) = \>\"", ",", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"dphi", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"{", RowBox[{ "\"\<\!\(\*SubscriptBox[\(A\), \(k\)]\) \ \[PartialD]\!\(\*SuperscriptBox[SubscriptBox[\(A\), \(r\)], \(T\)]\)/\ \[PartialD]\!\(\*SubscriptBox[\(\[Theta]\), \(r\)]\) = \>\"", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], ".", " ", RowBox[{"(", RowBox[{ RowBox[{"dtheta", "[", RowBox[{ SubscriptBox["\[Theta]", "r"], ",", SubscriptBox["\[Phi]", "r"]}], "]"}], " ", "//", " ", "Transpose"}], ")"}]}], ")"}], " ", "//", " ", "Simplify"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "\"\<\!\(\*SubscriptBox[\(A\), \(k\)]\) \ \[PartialD]\!\(\*SuperscriptBox[SubscriptBox[\(A\), \(r\)], \(T\)]\)/\ \[PartialD]\!\(\*SubscriptBox[\(\[Phi]\), \(r\)]\) = \>\"", ",", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], ".", " ", RowBox[{"(", RowBox[{ RowBox[{"dphi", "[", RowBox[{ SubscriptBox["\[Theta]", "r"], ",", SubscriptBox["\[Phi]", "r"]}], "]"}], " ", "//", " ", "Transpose"}], ")"}]}], ")"}], " ", "//", " ", "Simplify"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"{", RowBox[{ "\"\<\!\(\*SubscriptBox[\(A\), \(k\)]\) \ \[PartialD]\!\(\*SuperscriptBox[SubscriptBox[\(A\), \(r\)], \(T\)]\)/\ \[PartialD]\!\(\*SubscriptBox[\(\[Theta]\), \(r\)]\) = \>\"", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"tk", ",", "pk"}], "]"}], ".", " ", RowBox[{"(", RowBox[{ RowBox[{"dtheta", "[", RowBox[{"tr", ",", "pr"}], "]"}], " ", "//", " ", "Transpose"}], ")"}]}], ")"}], " ", "//", " ", "Simplify"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ "\"\<\!\(\*SubscriptBox[\(A\), \(k\)]\) \ \[PartialD]\!\(\*SuperscriptBox[SubscriptBox[\(A\), \(r\)], \(T\)]\)/\ \[PartialD]\!\(\*SubscriptBox[\(\[Phi]\), \(r\)]\) = \>\"", ",", "\[IndentingNewLine]", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"tk", ",", "pk"}], "]"}], ".", " ", RowBox[{"(", RowBox[{ RowBox[{"dphi", "[", RowBox[{"tr", ",", "pr"}], "]"}], " ", "//", " ", "Transpose"}], ")"}]}], ")"}], " ", "//", " ", "Simplify"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], ")"}]}], "}"}]}], "*)"}], "\[IndentingNewLine]", "}"}], " ", "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], ".", " ", RowBox[{"(", RowBox[{ RowBox[{"dtheta", "[", RowBox[{ SubscriptBox["\[Theta]", "r"], ",", SubscriptBox["\[Phi]", "r"]}], "]"}], " ", "//", " ", "Transpose"}], ")"}]}], ")"}], " ", "//", " ", "Simplify"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{ SubscriptBox["\[Theta]", "k"], ",", SubscriptBox["\[Phi]", "k"]}], "]"}], ".", " ", RowBox[{"(", RowBox[{ RowBox[{"dphi", "[", RowBox[{ SubscriptBox["\[Theta]", "r"], ",", SubscriptBox["\[Phi]", "r"]}], "]"}], " ", "//", " ", "Transpose"}], ")"}]}], ")"}], " ", "//", " ", "Simplify"}], " ", "//", " ", "MatrixForm"}], " ", "//", " ", "TraditionalForm"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.662846447248891*^9, 3.6628464510012665`*^9}, 3.6628466738734217`*^9, {3.6628478583077135`*^9, 3.6628481216252766`*^9}, { 3.6628481521203594`*^9, 3.6628482086887665`*^9}, {3.662849517012821*^9, 3.662849597794898*^9}, {3.662849656412856*^9, 3.6628496608263364`*^9}, { 3.662849761852664*^9, 3.6628497664267764`*^9}, {3.662849939287547*^9, 3.662849959670011*^9}, {3.6628499917639413`*^9, 3.6628501674070606`*^9}, { 3.662850210460614*^9, 3.6628502374823694`*^9}, {3.6628510008942013`*^9, 3.662851047399355*^9}, {3.6628513735774064`*^9, 3.6628514663275137`*^9}}], Cell[BoxData[ TagBox[GridBox[{ {"\<\"\\!\\(\\*SubscriptBox[\\(A\\), \\(k\\)]\\) = \"\>", TagBox[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"sin", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], RowBox[{"-", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}]}]}, { RowBox[{ RowBox[{"-", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}]}], " ", RowBox[{"sin", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], "0"} }, 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$]]], TraditionalForm], TraditionalForm, Editable->True]}, {"\<\"\[PartialD]\\!\\(\\*SubscriptBox[\\(A\\), \ \\(k\\)]\\)/\[PartialD]\\!\\(\\*SubscriptBox[\\(\[Theta]\\), \\(k\\)]\\) = \"\ \>", TagBox[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"(", RowBox[{"-", RowBox[{"cos", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], ")"}]}], RowBox[{ RowBox[{"-", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}]}], " ", RowBox[{"sin", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], RowBox[{"-", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}]}]}, { RowBox[{ RowBox[{"-", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}]}], " ", RowBox[{"sin", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], "0"} }, 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$]]], TraditionalForm], TraditionalForm, Editable->True]}, {"\<\"\[PartialD]\\!\\(\\*SubscriptBox[\\(A\\), \ \\(k\\)]\\)/\[PartialD]\\!\\(\\*SubscriptBox[\\(\[Phi]\\), \\(k\\)]\\) = \ \"\>", TagBox[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{"-", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}]}], " ", RowBox[{"sin", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], "0"}, { RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"(", RowBox[{"-", RowBox[{"cos", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], ")"}]}], RowBox[{ RowBox[{"-", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}]}], " ", RowBox[{"sin", "(", SubscriptBox["\[Phi]", "k"], ")"}]}], "0"} }, 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$]]], TraditionalForm], TraditionalForm, Editable->True]} }, AutoDelete->False, GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}], "Grid"]], "Output", CellChangeTimes->{{3.662851428939417*^9, 3.662851467095209*^9}}], Cell[BoxData[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "r"], ")"}]}], "-", RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}]}], RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"sin", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}]}, { RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"sin", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}], RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}]} }, 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$]]], TraditionalForm]], "Output", CellChangeTimes->{{3.662851428939417*^9, 3.6628514671452503`*^9}}], Cell[BoxData[ FormBox[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"sin", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}], RowBox[{ RowBox[{"-", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "k"], ")"}]}], " ", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}]}, { RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"cos", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"cos", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}], RowBox[{ RowBox[{"sin", "(", SubscriptBox["\[Theta]", "k"], ")"}], " ", RowBox[{"sin", "(", SubscriptBox["\[Theta]", "r"], ")"}], " ", RowBox[{"sin", "(", RowBox[{ SubscriptBox["\[Phi]", "k"], "-", SubscriptBox["\[Phi]", "r"]}], ")"}]}]} }, 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$]]], TraditionalForm]], "Output", CellChangeTimes->{{3.662851428939417*^9, 3.662851467169272*^9}}] }, Open ]] }, WindowSize->{1492, 698}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, 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, 9235, 247, 673, "Input"], Cell[9817, 271, 5764, 157, 147, "Output"], Cell[15584, 430, 2184, 61, 79, "Output"], Cell[17771, 493, 1994, 56, 79, "Output"] }, Open ]] } ] *) (* End of internal cache information *)