(* 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[     53270,       1281]
NotebookOptionsPosition[     51296,       1214]
NotebookOutlinePosition[     51796,       1234]
CellTagsIndexPosition[     51753,       1231]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "Rotation", " ", "of", " ", "general", " ", "triangle", " ", "into", " ", 
     "xy", " ", "plane", " ", "with", " ", "one", " ", "point", " ", 
     RowBox[{"(", 
      RowBox[{"O", " ", "atom"}], ")"}], " ", "at", " ", "0", "rigine", " ", 
     "by", " ", "translation"}], ",", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"one", " ", "leg", " ", "on", " ", "x"}], "-", 
      RowBox[{"axis", " ", "O"}], " ", "-", " ", 
      RowBox[{"C", " ", "atoms"}]}], ";", " ", 
     RowBox[{"example", ":", " ", 
      RowBox[{
      "COH", " ", "in", " ", "Vitamine", " ", "A", "    ", "ES", " ", "15.3", 
       ".2013"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{"T1", ",", "U", ",", "T", ",", "P", ",", "p"}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"P", "=", "2"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"p", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0.8910058089999993`", ",", 
         RowBox[{"-", "1.2609385707599998`"}], ",", "2.215514762399999`"}], 
        "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1.3556894324000002`"}], ",", 
         RowBox[{"-", "1.02001740302`"}], ",", 
         RowBox[{"-", "0.6789785517999993`"}]}], "}"}]}], "}"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5716510170347843`*^9, 3.5716510826640997`*^9}, {
   3.5716511262193766`*^9, 3.5716512306927595`*^9}, {3.571651274435237*^9, 
   3.571651362856192*^9}, {3.5716514621659665`*^9, 3.5716515259076786`*^9}, {
   3.571651568979354*^9, 3.5716521347299476`*^9}, {3.571652232401719*^9, 
   3.571652378308776*^9}, {3.571652412753636*^9, 3.571652478757352*^9}, {
   3.571652545307069*^9, 3.5716525770999246`*^9}, 3.5716526669248824`*^9, 
   3.5716527375618067`*^9, {3.5716527820062847`*^9, 3.5716528145791416`*^9}, {
   3.571658921974269*^9, 3.5716589240022726`*^9}, {3.571658982455575*^9, 
   3.5716589854039803`*^9}, {3.571659202322361*^9, 3.571659204989966*^9}, {
   3.571659361723441*^9, 3.5716594007391095`*^9}, {3.5716608398884373`*^9, 
   3.571660843523244*^9}, {3.571667293318076*^9, 3.5716673002132883`*^9}, {
   3.571677324192623*^9, 3.5716774997709312`*^9}, 3.5716781072047987`*^9, {
   3.5716781806965275`*^9, 3.5716782034569674`*^9}, {3.5716783802208776`*^9, 
   3.571678450873402*^9}, {3.5717165445765185`*^9, 3.571716567336959*^9}, 
   3.5717197685873814`*^9, {3.571719848802722*^9, 3.571719857772738*^9}, {
   3.571720700829419*^9, 3.5717207082550316`*^9}, 3.571721050207632*^9, {
   3.5723323770483665`*^9, 3.5723323792947702`*^9}, {3.5723325671659*^9, 
   3.572332567727501*^9}, {3.5723332412898846`*^9, 3.5723332508683014`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"c1", "=", "20"}], ";", " ", 
  RowBox[{"c2", "=", "21"}], ";", " ", 
  RowBox[{"c3", "=", "51"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"v2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"u2", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.2557`", ",", "0.255`", ",", "0.255`", ",", "0.2557`", ",", "0.255`", 
     ",", "0.2557`", ",", "0.2547`", ",", "0.2547`", ",", "0.2557`", ",", 
     "0.2547`", ",", "0.2557`", ",", "0.2557`", ",", "0.2547`", ",", 
     "0.2557`", ",", "0.2557`", ",", "0.2557`", ",", "0.2557`", ",", 
     "0.2547`", ",", "0.2557`", ",", "0.255`", ",", "0.178`", ",", "0.5693`", 
     ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", ",", 
     "0.5241`"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plota", " ", "=", " ", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"p", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c1", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"p", "[", 
         RowBox[{"[", "2", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c3", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{"v2", ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c2", "]"}], "]"}]}], "]"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Plotc", "=", 
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", "v2"}], "}"}], ",", "0.02"}],
         "]"}], ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"p", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.02"}], "]"}], 
       ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", "v2"}], "}"}], ",", "0.02"}],
         "]"}]}], "}"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb0", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c3", ",", 
      RowBox[{"p", "[", 
       RowBox[{"[", "2", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb1", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c1", ",", 
      RowBox[{"p", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb2", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c2", ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
  "Plota", ",", "Plotb0", ",", "Plotb1", ",", "Plotb2", ",", "Plotc", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Axes", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"xl", ",", "yl", ",", "zl"}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "\[Rule]", "False"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "Automatic"}], ",", " ", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.5716510170347843`*^9, 3.5716510826640997`*^9}, {
   3.5716511262193766`*^9, 3.5716512306927595`*^9}, {3.571651274435237*^9, 
   3.571651362856192*^9}, {3.5716514621659665`*^9, 3.5716515259076786`*^9}, {
   3.571651568979354*^9, 3.5716521347299476`*^9}, {3.571652232401719*^9, 
   3.571652378308776*^9}, {3.571652412753636*^9, 3.571652478757352*^9}, {
   3.571652545307069*^9, 3.5716525770999246`*^9}, 3.5716526669248824`*^9, 
   3.5716527375618067`*^9, {3.5716527820062847`*^9, 3.5716528145791416`*^9}, {
   3.571658921974269*^9, 3.5716589240022726`*^9}, {3.571658982455575*^9, 
   3.5716589854039803`*^9}, {3.571659202322361*^9, 3.571659204989966*^9}, {
   3.571659361723441*^9, 3.5716594007391095`*^9}, {3.5716608398884373`*^9, 
   3.571660843523244*^9}, {3.571667293318076*^9, 3.5716673002132883`*^9}, {
   3.571677324192623*^9, 3.5716774997709312`*^9}, 3.5716781072047987`*^9, {
   3.5716781806965275`*^9, 3.5716782034569674`*^9}, {3.5716783802208776`*^9, 
   3.571678450873402*^9}, {3.5717165445765185`*^9, 3.571716567336959*^9}, 
   3.5717197685873814`*^9, {3.571719848802722*^9, 3.571719857772738*^9}, {
   3.571720700829419*^9, 3.5717207082550316`*^9}, 3.571721050207632*^9, {
   3.5723323770483665`*^9, 3.5723323806363726`*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.5], 
    SphereBox[{0.8910058089999993, -1.2609385707599998`, 2.215514762399999}, 
     0.255], SphereBox[{-1.3556894324000002`, -1.02001740302, \
-0.6789785517999993}, 0.5241], SphereBox[{0, 0, 0}, 0.178]}, 
   Text3DBox[
    "51", {-1.3556894324000002`, -1.02001740302, -0.6789785517999993}, {1, 
    0}], Text3DBox["20", {0.8910058089999993, -1.2609385707599998`, 
    2.215514762399999}, {1, 0}], 
   Text3DBox["21", {0, 0, 0}, {1, 0}], {
    TubeBox[{{0.8910058089999993, -1.2609385707599998`, 2.215514762399999}, {
      0, 0, 0}}, 0.02], 
    TubeBox[{{0.8910058089999993, -1.2609385707599998`, 
     2.215514762399999}, {-1.3556894324000002`, -1.02001740302, \
-0.6789785517999993}}, 0.02], 
    TubeBox[{{-1.3556894324000002`, -1.02001740302, -0.6789785517999993}, {0, 
      0, 0}}, 0.02]}},
  AspectRatio->Automatic,
  Axes->Automatic,
  AxesLabel->{
    FormBox["xl", TraditionalForm], 
    FormBox["yl", TraditionalForm], 
    FormBox["zl", TraditionalForm]},
  Boxed->False,
  PlotRange->Automatic,
  SphericalRegion->True]], "Output",
 CellChangeTimes->{3.5723354149197025`*^9, 3.572335447024558*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"rotate", " ", "p1", " ", "around", " ", "x"}], "-", 
     RowBox[{"axis", " ", "until", " ", "the", " ", "z"}], "-", "comp"}], " ",
     "=", " ", "0"}], " ", "*)"}], "\n", 
  RowBox[{
   RowBox[{"T1", "=", 
    RowBox[{"ArcTan", "[", 
     RowBox[{
      RowBox[{"p", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "3"}], "]"}], "]"}], "/", 
      RowBox[{"p", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "]"}]}], "\n", 
   RowBox[{
    RowBox[{"U", "=", "p"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"T", "=", 
     RowBox[{"RotationMatrix", "[", 
      RowBox[{"T1", ",", 
       RowBox[{"{", 
        RowBox[{"1", ",", "0", ",", "0"}], "}"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{"MatrixForm", "[", "T", "]"}], "\n", 
   RowBox[{"p", "=", 
    RowBox[{"U", ".", "T"}]}]}]}]], "Input",
 CellChangeTimes->{{3.5716510170347843`*^9, 3.5716510826640997`*^9}, {
   3.5716511262193766`*^9, 3.5716512306927595`*^9}, {3.571651274435237*^9, 
   3.571651362856192*^9}, {3.5716514621659665`*^9, 3.5716515259076786`*^9}, {
   3.571651568979354*^9, 3.5716521347299476`*^9}, {3.571652232401719*^9, 
   3.571652378308776*^9}, {3.571652412753636*^9, 3.571652478757352*^9}, {
   3.571652545307069*^9, 3.5716525770999246`*^9}, 3.5716526669248824`*^9, 
   3.5716527375618067`*^9, {3.5716527820062847`*^9, 3.5716528099459333`*^9}, {
   3.571653332843252*^9, 3.5716533405340652`*^9}, {3.571653579417285*^9, 
   3.5716535853140955`*^9}, {3.5716667591731386`*^9, 3.57166689118057*^9}, {
   3.571667308590503*^9, 3.571667311601309*^9}, 3.5716674264019103`*^9, 
   3.5716675428717146`*^9, 3.571673022609868*^9, {3.5723327300457864`*^9, 
   3.572332744663012*^9}}],

Cell[BoxData[
 RowBox[{"-", "1.053376974006295`"}]], "Output",
 CellChangeTimes->{3.5723354149353023`*^9, 3.5723354470713587`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0.`", "0.`"},
     {"0.`", "0.4946389506381359`", "0.8690985608730484`"},
     {"0.`", 
      RowBox[{"-", "0.8690985608730484`"}], "0.4946389506381359`"}
    },
    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.5723354149353023`*^9, 3.5723354470713587`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0.8910058089999993`", ",", 
     RowBox[{"-", "2.5492100230547097`"}], ",", "0.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.3556894324000002`"}], ",", "0.08555894437059663`", ",", 
     RowBox[{"-", "1.2223448953982992`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5723354149353023`*^9, 3.5723354470713587`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"Rotate", " ", "p1", " ", "around", " ", "z"}], "-", 
   RowBox[{"axis", " ", "to", " ", "make", " ", "y"}], "-", 
   RowBox[{"comp", " ", "to", " ", "0"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.571657528938622*^9, 3.5716575406698427`*^9}, {
   3.571658470103875*^9, 3.5716585072475405`*^9}, {3.57165863791337*^9, 
   3.5716586547769995`*^9}, {3.5716595085508995`*^9, 
   3.5716595093153005`*^9}, {3.5716595400317545`*^9, 
   3.5716595405465555`*^9}, {3.5716595715282097`*^9, 3.571659962823697*^9}, {
   3.5716669073265986`*^9, 3.571666964407099*^9}, {3.571667060846468*^9, 
   3.571667121109374*^9}, {3.571667254910809*^9, 3.5716672710880375`*^9}, 
   3.5716675639317517`*^9, {3.5716728084526916`*^9, 3.5716728130547*^9}, 
   3.5716730371958933`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"T2", "=", 
  RowBox[{"ArcTan", "[", 
   RowBox[{
    RowBox[{"p", "[", 
     RowBox[{"[", 
      RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", 
    RowBox[{"p", "[", 
     RowBox[{"[", 
      RowBox[{"1", ",", "1"}], "]"}], "]"}]}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"U", "=", "p"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"RotationMatrix", "[", 
    RowBox[{"T2", ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"MatrixForm", "[", "T", "]"}], "\n", 
 RowBox[{"p", "=", 
  RowBox[{"Chop", "[", 
   RowBox[{"U", ".", "T"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.571657528938622*^9, 3.5716575406698427`*^9}, {
   3.571658470103875*^9, 3.5716585072475405`*^9}, {3.57165863791337*^9, 
   3.5716586547769995`*^9}, {3.5716595085508995`*^9, 
   3.5716595093153005`*^9}, {3.5716595400317545`*^9, 
   3.5716595405465555`*^9}, {3.5716595715282097`*^9, 3.571659962823697*^9}, {
   3.5716669073265986`*^9, 3.571666964407099*^9}, {3.571667060846468*^9, 
   3.571667121109374*^9}, {3.571667254910809*^9, 3.5716672710880375`*^9}, 
   3.5716675639317517`*^9, {3.5716728084526916`*^9, 3.5716728130547*^9}, {
   3.5716730371958933`*^9, 3.5716730409399*^9}, 3.5716732176102104`*^9, 
   3.5716733652332697`*^9, {3.5716734033597364`*^9, 3.5716734637786427`*^9}, {
   3.5716740082975993`*^9, 3.5716740130400076`*^9}, 3.5716740558776827`*^9, {
   3.571675848820032*^9, 3.5716758536404405`*^9}, {3.571716495732833*^9, 
   3.571716500225641*^9}, 3.5717171313559494`*^9, 3.5717172354861317`*^9, {
   3.5717173286962957`*^9, 3.571717332175102*^9}, {3.5717173666043625`*^9, 
   3.5717173730627737`*^9}, {3.571717464400934*^9, 3.5717174651497355`*^9}, 
   3.571719768602981*^9, {3.571721118722953*^9, 3.571721158908623*^9}, {
   3.5717221478099604`*^9, 3.5717221537691708`*^9}, {3.5717221941108418`*^9, 
   3.571722237026517*^9}, 3.5717222852306013`*^9, 3.5723325444522605`*^9, 
   3.5723327037909403`*^9, 3.572332756909034*^9, {3.572332804504717*^9, 
   3.5723328097151265`*^9}, {3.57233319923221*^9, 3.572333202430216*^9}, {
   3.5723353622852097`*^9, 3.5723353707716246`*^9}, {3.5723354015036783`*^9, 
   3.5723354025644803`*^9}}],

Cell[BoxData[
 RowBox[{"-", "1.2345471182378898`"}]], "Output",
 CellChangeTimes->{3.572335414966502*^9, 3.572335447086959*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.32994867848144904`", "0.9439988715927289`", "0.`"},
     {
      RowBox[{"-", "0.9439988715927289`"}], "0.32994867848144904`", "0.`"},
     {"0.`", "0.`", "1.`"}
    },
    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.572335414966502*^9, 3.572335447102559*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2.700437574415364`", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "0.5280754835921541`"}], ",", 
     RowBox[{"-", "1.251539233788441`"}], ",", 
     RowBox[{"-", "1.2223448953982992`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.572335414966502*^9, 3.572335447102559*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"rotate", " ", "p2", " ", "around", " ", "x"}], "-", 
    RowBox[{"axis", " ", "to", " ", "make", " ", "z"}], "-", "comp"}], " ", 
   "=", " ", "0"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5716606820473604`*^9, 3.571660691968978*^9}, 
   3.571660959275447*^9, {3.5716669809899282`*^9, 3.571667011753182*^9}, {
   3.5716671327313943`*^9, 3.571667190357896*^9}, {3.5716672271895604`*^9, 
   3.5716672372047777`*^9}, {3.5716674632179747`*^9, 3.5716674704407873`*^9}, 
   3.5716676177986465`*^9, 3.5716728498551645`*^9, {3.5716729371217175`*^9, 
   3.571672942004526*^9}, {3.5716729911134124`*^9, 3.571672996667022*^9}, 
   3.5716730605959344`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"T3", "=", 
  RowBox[{"ArcTan", "[", 
   RowBox[{
    RowBox[{"p", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "3"}], "]"}], "]"}], "/", 
    RowBox[{"p", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "2"}], "]"}], "]"}]}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"U", "=", "p"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"RotationMatrix", "[", 
    RowBox[{"T3", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "0", ",", "0"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"MatrixForm", "[", "T", "]"}], "\n", 
 RowBox[{"p", "=", 
  RowBox[{"Chop", "[", 
   RowBox[{"U", ".", "T"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.5716606820473604`*^9, 3.571660691968978*^9}, 
   3.571660959275447*^9, {3.5716669809899282`*^9, 3.571667011753182*^9}, {
   3.5716671327313943`*^9, 3.571667190357896*^9}, {3.5716672271895604`*^9, 
   3.5716672372047777`*^9}, {3.5716674632179747`*^9, 3.5716674704407873`*^9}, 
   3.5716676177986465`*^9, 3.5716728498551645`*^9, {3.5716729371217175`*^9, 
   3.571672942004526*^9}, {3.5716729911134124`*^9, 3.571672996667022*^9}, 
   3.5716730605959344`*^9, 3.5716732321182356`*^9, 3.571673520406742*^9, {
   3.5717164705855885`*^9, 3.5717164745323954`*^9}, 3.5717171518855853`*^9, 
   3.5717197686185813`*^9, 3.5723320674970226`*^9, 3.572332553812277*^9, 
   3.5723326767248926`*^9, {3.5723327708554583`*^9, 3.572332777220269*^9}, {
   3.572332818185941*^9, 3.5723328223043485`*^9}, {3.5723332139586363`*^9, 
   3.57233321614264*^9}, {3.5723354316897316`*^9, 3.5723354400513463`*^9}}],

Cell[BoxData["0.773597698323206`"], "Output",
 CellChangeTimes->{3.5723354149821024`*^9, 3.572335447118159*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0.`", "0.`"},
     {"0.`", "0.7154015443235817`", 
      RowBox[{"-", "0.6987135538827298`"}]},
     {"0.`", "0.6987135538827298`", "0.7154015443235817`"}
    },
    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.5723354149821024`*^9, 3.572335447118159*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2.700437574415364`", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "0.5280754835921541`"}], ",", 
     RowBox[{"-", "1.7494220465679622`"}], ",", "0"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5723354149821024`*^9, 3.572335447118159*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Plota", " ", "=", " ", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"p", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c1", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"p", "[", 
         RowBox[{"[", "2", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c3", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{"v2", ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c2", "]"}], "]"}]}], "]"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Plotc", "=", 
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", "v2"}], "}"}], ",", "0.02"}],
         "]"}], ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"p", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.02"}], "]"}], 
       ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", "v2"}], "}"}], ",", "0.02"}],
         "]"}]}], "}"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb0", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c3", ",", 
      RowBox[{"p", "[", 
       RowBox[{"[", "2", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb1", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c1", ",", 
      RowBox[{"p", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb2", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c2", ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
  "Plota", ",", "Plotb0", ",", "Plotb1", ",", "Plotb2", ",", "Plotc", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Axes", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"xl", ",", "yl", ",", "zl"}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "\[Rule]", "False"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "Automatic"}], ",", " ", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.5716775916862926`*^9, 3.5716776153203344`*^9}, {
   3.57167792031647*^9, 3.571677980626176*^9}, {3.5716780600303154`*^9, 
   3.5716780606231165`*^9}, {3.571681187032473*^9, 3.5716812002300963`*^9}, {
   3.571716393895854*^9, 3.57171643149192*^9}, {3.5717200555342855`*^9, 
   3.5717201809585056`*^9}, 3.572332418092038*^9}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.5], SphereBox[{2.700437574415364, 0, 0}, 0.255], 
    SphereBox[{-0.5280754835921541, -1.7494220465679622`, 0}, 0.5241], 
    SphereBox[{0, 0, 0}, 0.178]}, 
   Text3DBox["51", {-0.5280754835921541, -1.7494220465679622`, 0}, {1, 0}], 
   Text3DBox["20", {2.700437574415364, 0, 0}, {1, 0}], 
   Text3DBox["21", {0, 0, 0}, {1, 0}], {
    TubeBox[{{2.700437574415364, 0, 0}, {0, 0, 0}}, 0.02], 
    TubeBox[{{
      2.700437574415364, 0, 0}, {-0.5280754835921541, -1.7494220465679622`, 
       0}}, 0.02], 
    TubeBox[{{-0.5280754835921541, -1.7494220465679622`, 0}, {0, 0, 0}}, 
     0.02]}},
  AspectRatio->Automatic,
  Axes->Automatic,
  AxesLabel->{
    FormBox["xl", TraditionalForm], 
    FormBox["yl", TraditionalForm], 
    FormBox["zl", TraditionalForm]},
  Boxed->False,
  PlotRange->Automatic,
  SphericalRegion->True]], "Output",
 CellChangeTimes->{3.5723354149977026`*^9, 
  3.5723354471337585`*^9},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztfHd0Y9d55y+aUbNkFTfZlh3LZbPZ3SSWLMeJEpfEiVtsb4pjeU+yycaR
XOQ4ciw3JY5saSLLGpWRNUUiORxyyBkOy3BmSKIRBAgCJArRSYAFYC/DDs6w
dxL7u+/dRz5iQEo+J2fP/mEcfsTDw3v33u93v3rv9/DFr37/m1/7569+/5GH
vnrPZx776r9885GHvnfPp7/zGE8d+DXg175J+rN7II7TgPYvzb/0I+KfOPOr
418d/ycev0P8OyDk7Iu6s+9SjpXTypGgP/slvn+n+HddZqt3Z7nr879Eq9nH
+p+FxLtl7ysrK+jp6cHCwgI+l8kTblC+d7lcePDBB/GZz3wGL7/8Mh566CH9
pXftcelTTz2FkpKS62RPPITZbFY+fUZ3+1uU45u3bz98+LByu7h2cHBQ+TYa
jSr3a23Nzs4q3YhR69t6kxzKxMSEcr245vHHH1fu57UHdLeTCeWT1qToWpwT
TX5a1+TtOu60JsXwZJM3ZGnyQBaOP6Vr8tY9ABPXssmbdU2eP39eN8qDyi3i
nHjxUn2rN+/TKkd2h/xawKGdHhkZ0Xi+TXx9UDT7IV2bN2TMTcZIlW/FMMXr
LVm4vi1DxjPbP5Ax5oy5v17eLoROU5E3ZYiAmIcbtwd/MItS3afrkJeKvoR8
iL4effRRBRAhZ4Kfd+sa14Th9ixN3qs22b7d7wf2l5L3vg5sdO12ZIHtvbrh
Z06FYKWtrQ3vz4LNLVJy9N3szEJnlkG8e3/x1BSmoKBA6VR8en8W3G7Ztor3
q30lsvT1jv370gRAs1Di0wfkOSnXyqdP7N3D2/YXsd9Wvr5+1wDES0D3Cd0s
X4vgW/eQJ/HS5Om3XsO2CbyEmRTXflw399ftKNSNWWVK09zfkWMTDQktkR1r
FkgzauLTx9T2j2+3L+2U9tK5mtece721+115ThM48emjassnssjxdVnm6HP7
96rZWE0CBH1YNxINQdnrK3JK9T1ovX92/57+OAtuf/jL8fIne8ib9Iaf0lkC
8bX49AfXzMz1u2ZGO/6kTpr2smCap/yUDh/NU8p+8rKM+uNy1K/hOBOJRIZH
29FD0cMDag+525x8bH+8vyCnKjNsEF8/oNe+PVv4c/l1pkeSPP++2khsH+Op
zcSXdTxq/Pyeenv09WmF+KjN6Zd186wh9hG1tdAeEiL19/p9WtMbdtlaYP+Z
uzFLaw9msdaqFsOXxY+JlxgjL3uD7k4tNvtKFk5V7YRbebtpz/BBDuprWbiT
Lbj2R/5AFu4eymLRNS6lJ3pVvB3I9NtZ2v8n+bUYWebAc3JyZHM5WRTqnv3d
zveyGM379m7t3TqeNM3XJvo6gfwPspjDfdq7W2dJsoQryreaJflhFpt4795N
v31/RA/qBqpJ0A+zgPHBvXt46/49PJkFig/q5l3nW69tQ7x4w9NZeFZdLSoy
ZDpbusKX2sLB7ctkTCEbKdXZo70suSbbonGN9ad01k7rVCAmrd3d2/yp0T5K
sqizbp6PZLApus2Wlr5Rba14f+CPZjEEd0sYMnyx0ryasKDwddjm41m82Q6z
B7IM+Sa18fxf3nJrfen9wDtfT1+5+6t8tlj2RBbI3p4Fsj1iDdnxK69Dmgqz
2Ni7trm6bvtIFxSqJkL5dDyLkxGsiXeeK9YJu8hOJWzv3nvIL+sM0F6JpTbR
RdL4ZQqGeL99W5resM2C6qRw5PUZIv3gTstb9IGmDrBsmaBqjPD8/p3dmOXO
s7oZ0eRMiXT3TGRlX8/uY7s43ZqP1ifNJTo5o9/SKXXWHn6ms3FCXDMFSmYf
xix2+9a9W/3b7Rn6jSzipItZsmUZhixo7dPX32RXovftP032LGZ/R67u2rfp
92RpWrxk0w6diuj8wfv3bfNd+5hGMS8Cs3rdXGVeI+ZJusD37dvRa6Th2oRo
IYH4pGEl5ELT2mxp53v37fg1svNsmV+jPCcu1aTunkwJ02KEPSQsnGENRfuy
je/oRpTpMslTaxYW1anHt1+HTRPjEeOIZTH8v64288jrs1v6AKo1Czuyta/v
31pyD/x3WXAZUmefJb565Neiey0MEkOQzMlG/kEHT6aX0sed4jWQwZBoUVUG
/P0+86ODR2/9+rPYjrt1ErM3PJNZpkldGcf/2uNO8d7W1jadwYEYnrzzy69D
TrTFtStZhq4qK/5qf9XRgBDwCrjFJ21Q4hYBvXAFsrU/15mavWKIG7I0mVL+
73bOGr9qJLNrTUk/YzLW1pa7NXx3Wr3WEUsIpCX+TfF2YMf47p1oSrXb1HWl
Ka9qgfDfsviSu/cXDv2iSDAYVD5pPehdonRJiiHc07e/RtqkXKRFEfrM5q17
Dz8zT9IERDRtt9uVi27JYs7esneTb9pHhzUp2V4XFtfemgWPN+/dvrZJse8y
1F3KV5rtGw4lcTnUhRFSOMeIFlKM1J5jQicpmWtCd54JvSdN6M83YbDAhOFC
Ey6fNmK02IjpOh8WYl1YjHdhKZ7EghRoDTHBnqZ5b9bJneJmbt9/1g5Kuz49
kkIvxxmt8cHFsbk5Lh8pkGPmmC0cs4VjrkVbbi0682qRzLOi+6QVvflW9BdY
MVhYi+GiWowU12L0bC3GSyyYOGfBZKkZqTIzpstNuFJhwtXzRiy4fFjt6ML6
ZOqqZEXYQs3CahMtWfnIPlaMk6Ay8AFMkYEkATZx8FYOvJ7k5KCbOGgvyZ9j
RSinDlFSLNeGtpM2dOTbkTxlR3eBHb0F9eg/bcdgkR3DZ+wYOWvH2Dk7xkvt
mCyzYaq8DqnzdbhSacXVC1bMXKzF7CUL5i6ZMV9lwkK1ICPWk13YmErN6Myy
lrbdqTJ0/2tYAulhVEl6H9o5K401XlSTMSNno5ZkI0MNZKiRzHhybGjOsSOY
W49IrgOtpLa8BjLXgMSpBnQXOtFb5EQ/afCME8NnnRg558RYqRPj5U5MVjRg
6nwDpi84cOWiAzOX6jFbZcdctQ3zNXVYqLFi0VB7HZaMFs7UsslEMiq00d2F
zVQKmZ5EcCGZvXcf8dt2AR/c9oVtZLY8x4BKMnuJjBo4gxYyWkdGHWTSlVNP
0XRQNBsomk6Ec11oIcXzGtF+sokMN6GrsAk9RU3oP+PG4Fk3hs+5cbm0CaPl
TRivaMJkZRNSF5owfakRV6saMVPtwmyNC/MGJxaMDVg0ObBkqsey2Y5liw0r
ljqs1NZitdaCVasZa1YTiTMd9GJzOiU9we4lCrm48ds6l6VlKYL7a13WfZig
ADs4z2VkvYLCe5GsV5N1E+fXSrbtZNtJtptyXBToRvhzmxDKcyN60o1Yvgft
+V50FniRPO1FT7EXfWd9GCjxYajUh8vlPoxW+DBe6cPkBR9Sl7yYrvLiarUX
MwYP5owezJvcWDC7sWhpwlJtI5atLqxYnVita8CarZ5kJ9mwbreSarFhN5NM
PCYUNgO2+pKYyhKMqAYI/2MPuyldrJp0f1giYFBQKCMClUSgiggYiICF3NvI
vYPcu3KaKPluSr4XwTwvIvk+tJ5qRrygGR2FfiSK/Og+E0BvSQD95wIYKgvg
ckUAo5UBjF0IYOJSAFNVAUzXBHDV4MeM0Y85UzPmLc1YqPVhyerDcp0XKzYP
Vm1urNmbsF7fiHWHCxsOJ8mBjYZ6bDbYSFZSLclMMpIM2Ix4sXUlNZFFM27X
meh9zIAMJlRgPqIDxkyyUkRsuEBQqgmKiaBYCYidgDhzPBQRH0WkmRriRzgv
gGh+ALGCINpPh5AoDqHrbBg958LoKw1jsDyM4fMRjFyIYOxSGBNVYUzVhDFt
COOKKYwZcxhzlhDma0NYsAaxZAti2R7ASr0fq/XNWHP4sN7gw4bTS3Jj09VE
asRWo5PUQKonEaRGK98tJBO2XATISZEZSI5JB5BtKfymbSemmgn8l9cX7eir
ElTtemBbs4RMldGYlBO684TuIqGrIXQWwmajLDkIm4uweXIClK0gZStE2Qqj
tSCCttMRdBRHkTwbRfe5FvSWtWCgogVDlS24fLEVo1WtGK9uxaShBSljC66Y
WzBjacFsbQvm6+jcbFEs2SNYro9g1RHGWkMY684QNlxBwhPAVpOf1IxNNw2L
24MttxtpD5NJD/2Hh1B6CKXHxmMr0m5ql9vI6wllI6G8msJl7CQ5Ws4rrIzU
RjU2uya/latpr7HucJuM13ZwrCXZiKODouiknjZRTz30UF56qGaKYoCiGKQo
hiiKYfjzogjnRxEtIG61PZhsHUcqNo7F8Xm5SPJOrPVNYL1vHOv949gYGMea
qwXrjVFsNEWJR4R4MCn2hJD2MpL1Ei9fM8mHdLOXRKyaiVWzk+QgMfv31ZEs
vNbE+4gTsdqaSY1lyXtu2BebbAsnuuTyDhkYqtiYJDZ2YtNAbFzExk0581LO
mmEmLnXExEFMXDkRuHOjGA6PY3Z0/u0S4q2ZeWwNj2OrnVlrMwPS5gjSfrLu
J+uBAMlPIttBD6mJRLaDZDtAtgMUD38t76F4+GiKyPbWUBJ6njX1un5fnt8t
zZMWZGoLRVpQK/MV1Zi9UfFeO3JhJznJfyN1zEMz1UxnHqAND1LPwtSzKPpD
43I5gwH73DzSIxxiiIyGyGiYjIbJaJiMhslomIyGyWiYehAikyHObZA6EODc
NpNJL+d2KKkuMdyqhLJ9jCMGSSL0DtJqRklxUifpsoFOL5zELGllNKUtRup9
1sF9oXmXTDEEFNpLgCpg0taetHUoFZ9bpOnewaaM+lJOmThPW3ORuFQTF19N
Nwc/L80cpWqeuIwRlxZi0kJMWohJSzOJ8t5CeY9y4qMNSEfqiQ1lPUjzGmBY
1kw8ZlJhCa8WQBo4BCupgeQhBSQubaRErgE9eQYMnDRguMCA0UIDxosMmCw2
YNbmxdp46sYsarM/Tjv2O9uWxM6qvnL+egmqyM201177fbfJa0W0qBp0DVQ3
yUdg/VS6EJpqumTcwTsmCGSMmMRoPOIUrBgFK0bBaiWIrdSeFoIYJYhhghg0
Yk0F0S/zOSHfO8GpgbbOQFk2MM0w0saZ6CvMiJDiuWZ05lrQc9KCgQILLheZ
MXbGjMlzJuZBRuZARsycJ6ikRZdX5EDILNAT7B7M9H17WO89ZiBb/n5dFoQF
msKqaQhrOz+3Se131DQr7rKMpkyIrEC3jKrcUJPE5Mjc70iVS08S3HZKaDsl
tI0S2kZg26ixcQLbSmBbCGyE8SKBXe/wNiu3vUPOn0GhCtIlJcUx0UCaaRxF
3iZyNhtTGjvTmXp0KTmag/lZPXOzekyW2zFdaWMuZsVcVS1TFTOWDMxNDAS3
xoAlOlINoIwAQ+7A6PaSVHHUZPI/AekbtuNbvTzTaNwhs4KRgWEpvvW7wG0L
jcrit1uQnhpFupMi20HQOqj37US1jajGiWorQyqiuhY0YH04CY9MhVVbbNCR
kXbGRBtTS5mto8zWU2YbmEG4EGHi1MYMIlHgRm+xB4PnPBgp92Ci0s0EyY2Z
mibMG11YMjcwCapnAlSH9ToLychjA1boe7O5RW3Ll+836wT44A7Au+uUDuy2
Bge3j6Rl1qDVll9afTEl/iijf9VwmxiZU2syOLmLc0gPJ5BOUIEThK2TsLUT
tja6jlYmL1ETlvwGtfrjniyIqV78Ag2LgUjV0ZM5aVy8SrbRzGjQj47TAWYY
QfSXhnD5fAjjF0NI1QRx1RRkFhHAUh0DZDuDY4ebmQIDYqeDgW8d35kzOQyY
ioazBa7aqtlNr0Ptr888vXPxLgCFo9IWIe6Qd7aFeqRWe0h+AZ6sJeJ8TY8g
3UXkkkQuQeQ66GTabNiKMQMOc+IvJ+Uex3uyQGdRJqZCCa6bmJt4GQz60Zgb
RHNumEFgBPHTLUiUtKKnvBWDlTEGzjFMGGKYtsQwW9eKRXsLVhqiDAAjDNpC
DHoZGHs5Fh979VHwvRRA5hCa4dQvDd74uqUtm+5r1b4iusnUXE7WnRLg4iNl
UmObBXAPSGuZXppFupu4dXGsSY61k8FKmxVrLZS2Tm+jlLbdOa6V5FBipcoc
scrjp46K+DAK98lWBAriiBa1o62kA10Vnei7kMBwdQJjpgSmahO4ak9g3pHA
sqsTa+4ObHjbsNkcY2zUosZSIQ4nRMsQalDipfn2uIaafsHzhkzUDm4fZWy3
62A6qNyue71ZtnzmyAUy5SU0s7I28UYVm14asW4OJklZ7LRjI27BQsggS6je
q5Mlk87ruClLzVTFEAPoFkYwMTTktsOT34lgURKt57rQUdGD7ou9GKjpxWVz
Hybq+jBd34dZZx8Wm/qw4u3Fur8Hm8EubIVpFaIddElx+v8WNR6IiYCqgXJY
ussZa7qoR2cvXbwhqy5uy9+1tZ4CO62eQTgK7fydcgBFzxRRUtQ6649qk6SA
yMClm0qZbMBmuxVLLUaXjPRUf6oJlVPqdoD4RaiLMfrWdtSdTMB5qhveol6E
SvrQ2zCGqY4ZXEnMYKZ7FivTywdkTL8+PIuNkRlsjZK66eFb++h3mJS0J2kS
OLJOylknZazDh77zFZgfGc6qkrJUVGn2pteN2I3XICaCE211RBO1iCuA4mfU
QoW3bdu8A9tH79Ej18dp7m7EJkVvJW7CYDDglF+XHi2TlqtBkdwyhuRlOS3q
XB/EYPQqZsaXNTm/VdP3xWV6GQLTTVB6KFc9lKnuMOJnKzERDWsJln5Dne8S
jszAQe1pl+rt1Ly8U14cdYURqq7B5tUhuS2WvYj4vRrwy0Jg3NjscmC1vRZD
TKucMmgrO1YhI2eP5LeV1E5TlFTaSI0sXTsruwOdW7bRJRDjdBoDSUx4mhEv
r8ZEvFWbZ7H9pN/wlythMrK9PkMatNf+cVe24bxjG6UIBj0WbIzF9kXpfbLp
9BiNZo8Lq511GI5IobgZZUcrCUijokIUBlIHqYtmamm7iT2Cxdvl6cn+acSr
GhG/aEZracFNGdqxU5FwTYq0e/5v2D7aSQzeoRlhZxSLCYZlvXQwyzNvl8KZ
jWHxfpEcWf782M3wPXjiTngfPHE9/x1nOy3fyMXg4SPoffrniHzrubeJ078F
75eP38V/x26H58vH3gPPXx/7MP8dvR6eLx0lrM9/8mekX5BySAXI//xZNiW+
U1nxfOllyYo4Orh9TtBcQeWtkgdHzQAMhdScGitilWe1DVd9Tn6dDqdsjum/
yqaG29qw2FmL9R7O41j8bbuR002VqXoVwaqNe5CwrbOtIXs/Rqx29FgSuF+c
+0skbesfRbJunTed+IGD5Ce1ouTwwPXKaZVFcaSyqF6qnRO03jv3Bi2koNQ4
aoYoQSJyDaH+nIOiUaZ/lErLwiSrN2WIhKDfkOK8NJbEYrsZ61312BITPzMs
d1N3X69i9IZd5wpcWzjXloZ/KI1Y/wZiLUNoqbajqTZ5N+KTyoaxh5lXJK58
+qz4/yO0KV8c/ZkZR59hmn0qzk/i3K1K208cqsHhQ7XKcf6xUeW91JhGeaM6
QyXtaeXq2XX5FNiN9FZXCITwWG2kiAzjCYqhQuafqsdVJu83JYpLo+Q7xgg6
YcNmN2e4zyOM3Fu3L9zRmQOyn/hwGpXxNFpG0lhJrdNldWPjXB5GXzy2bVfD
FWcw2KoWqjfI+8qeoU94lonrYAfi9T6mUAa88HAeHvvdF7cd9+n/mEdp0SbK
rentcyOLO0p3kzQlU5z8hpoR8peQ/EaVqFtNXyjzDS6N5/tk94sjScy3GLHW
UYvNJAPJnkbVg9Gkv0de0zmRhqU7jcRYGmuTqwxk2rB66iiWmaHdJVlz1LTD
WlSDgVa/NP234Oijecj57il87d5/x7/+ITX3f7qR+88jKH5xGWXlaZS50rAl
CdmUVi53C/ojV2HJ7UUVLWBFTieHHZfT5pexsF0J3EQQIln5kBzmXNiAtTjT
/07mEF0cRI9bDWPolO+R1/RMp+EaSKOfLG2MLymLgIsvP4f5GqPc7j+oxJPC
W9UWXcLVsWG1ruw2ZU3m2YeP4lsf+g88/sfFeO4rrch/PIWSV9dRXpPGOR9l
rzWFUNsVpCeu0lMxnplbVsO5O7A4tYxJRkBjbTPwn+1HU2Ev7HldMHOuLtEf
Vih+MSwZdct03aJEqr3xxC5eb8JMgElxqwmbbbXMnerVWFfEvCL2XZrVnpEb
nEnjVNMWRojw5ui8sn599ZlDmKow4M5d7PoUH2Q5XVknDUA8mJCymIvvUhaf
+bQJJ/6xm7I4h7wnWxRZLApQtWPDGGnrYsjGMK21nRFuAluRLkbBPYyGe7Hs
7cNCUz9mnP2Yqu/HmLUfQ6Z+9FX3IXGhF/GyHoTPdqO5kNleXidsxMJE43WR
sluhxAs+OfEKHtcJQG7Xg3EzVinDyxEjNpg4ptvr1IRIJEbdIhIjIFdGpBO+
EdbONEqjaUymqDqXOUuNNgz/+AmsJRPb7lSdAo9iJkynKyzK6d9UEBFrRw/f
+2P86x+8SnFuojhfzirOm81GypYd6QiHEeEQIiEmR8yCgjFs+tuw4evAqrsT
S64E5hqSuGJLYrI2iVFjEoNVSfRUJtBZlkDsbAdChe3wnozDSQmpIyoGmvUL
HGK5MkSXVAkVGWEtQ8ltw/BGPVBvwPpMCotiGSdqohJTctptatYtsm+RhXdx
pJcT75dYufrTMBCvq9NbNE+TSNuqsVJagKW2xHZMJhJ6R02nDGddSlLeFuo2
SnUWSvPzh19WlOanf1SCX3yl5VqlGU0jNKU5kXdhxNHIPJwjbLbJjRbKta9Z
2ZDZdIex1hjFckMrFuwxzNTFkTLHMW6I4/KlOPrPx9FVGkP7mVa0FLYgeDIK
D7WqgZIkNoyqidl5JdhyyNzPpOEm8mUNN6lt90sulju8WA0ZsRllkBCzqus7
HU51vadTrF5wdElmeKlRaOD5L9PG96Yxl9qkQR3BlrEMC6dyMRdNSHt9owKe
Gvb5ZAjoUBaE2kJdqJZ9Cx18/uEc6uDz+NmnjTjxVergoTmUnt7c1sE4rZnw
B6qteSdW6w3Ychmx1Vir7owyENx0urHu8GLF3ozFugDmLEFcMYYwWR3G6MUw
hirC6D0XQrI4hHhBEJG8APxErFFZWXHTTjVS+hs4ozZpl4wHJHIioc6K3Bux
QYlb9lPiwgI9i7KYqCwqisVFscgoFhvFomMnveJgUllR0xA869tC82AaS6kN
SuYANs8XIHX0BKaoh9fLaxw1SR2CbplZKfOKixJBobcXOdCH7/03PC709ouN
yPk29fYFqbfOHb1Vxfrt2JhOYcVqUEoORAnCus2KtTomdNYGLFlcmDc24WqN
G6lLHoxXenG5zIuBEi+6izzoOOVBa64bIVotL1XCqSzt1VNv66i3AjnTgSwy
J4s5PixhFEuta0RuI2SiPbWoi9tikTtGfY03ysVverc2P5Gk7PXTBi/MaXZO
hCMXGI60MhxZFeEIs+n14lex5rBuL9F3hMZoYkMZADrliqdFGd15CaIYr4Gi
+ciHDuGnnziLXzwYRf6PplDyyrWKrNrRt2EjlcKSwYBlo0Ep5FkymkWND+aq
RWGTDdOVdkyU12OEMepgcQN6CxuQyHegLdeBCI2an+NoYthSr6womxkTmIQ6
HMhi6GQdkPrw3u1Ya/fS+RjohEzKVks6QvBaxKIswWulq2x1q8sxMYpejKIX
C6nbNuNj0lPcsR31JMfTWJ9c4T2tWM55CeFDp1BJc1KheAmvdE5O6aBqlaUR
akSZBK5C6u+/fJj6+ykDjv9Dl6K/57Lor7rccJeyXbPY6MXcBYNCM5VGpawt
VWbCxDkzRs9YMHTagr58C7ryLGjPtaCVExYkSF6lIk3sbhiVvTgh+eVy1VQg
JuskPyJhW2/z0lnRXARMys6essMXEY6LPjRKpKJU0ihFLUpRi1LUokQrSrSi
EfrWbqRHx+UywBsZ+Szgkm8eA5MiuFpkewHMv/hzeP7tFR1KdrlwZFTkqySL
kv7ogVfw3BddyPmn4axKqqakb8baBI1LPInUWQMmz6hbiSOFBgzlG9CfZ0B3
rroVGyOFST5SI8lOMisbQQY4OQYBi6ykVdem34INAYuPsAg/HjCrG8Eh+qMQ
NTBE2xWmBobd6qZxyCcXQKmFIUITpCAFIkoGkL48jq3Z+XdKmZoZmcdgaJzu
KMpAJ6KswZaf6UQBY7RRJUabo5tzYerJJ5H37RNSnkzK5FGminXKWMP5fuS+
Q/jJJ87gpS9HtpWxrJrK6N1RxvfIyVkZTWEuou5H9/H+LlKHBCdC8pOaSP0E
Q5RopmjUZZ3272tGMe5VtvrFlr+y9d/MEMJPeQlQXgJi5VeUBxCYAIHxU7v8
BMZPYJopM80Ex0dgvBFseqPY8JDiPdgcHMfmAOnq/B0S+4WxeYy3jGMsOo6E
qQfN9OHuvAhcuWE4CFgdzZVZ8eU+2lMPXjoziLKWNKYY0k2EelH2+HGUPPoM
Sp4+LUA7vUsRX8V3Pvwcnv5UDRUxicKnZnGu8FpFlL/qoq64f1SUcyhlHaK8
I+0xKjsDaR8lwmdTdwp8ZN7nkjsHZN5LqfCQeQ/Ni9uPLXeA9wex0RjCmiuM
VWcEqw2MYhxRZRdiwdaCuTpGNIEeLHePY6lrHIvJcSyQ5hPjmOscx2zHOGba
x9Fb3oru0lYkSlrQXtyC+OkoWgoiCJ+MIJAXhjc3xHkMUsQCFPRmRj0+GIlS
FRWwkj68XPGQ9ThxvgeHnrbg9BEj888yTJSVYqk4n2lv4mSGTj50779SJ0/g
uS84VZ18fukanZQ/zaM+Wf8xpVBoq9XL6IOINRGxJpoXN1Fzi6IiouYmam6i
5iZqTYz3mpp4rQebtHqbLlHn1Yz1Bj/dlagBE7VgISzZQlioEzViYcxawrhq
jmDaGMGUIYKJalFXFsHIxQiGKyMYrIigryyCnnMRpQ6tsziMtsIwWk+FEDkZ
RDBXFF6JAixRiOWjXHkpVx7KVRM1y0W2nRQYh4x2hLkyKybrxedoUb0pzIuo
rvcytqrPYjbnVVwNJ+TThzva+c37nqJ2FlM7w8j/4STOZtHOv1Bh+6gmaANJ
iLq+LacM4FxmEoM4l5XEtNZVz0DOQXJig8n8RkOTsvO3Vu/Bqt2DFZuXWbmo
KxT1hX7MW/yYNQcwYwzgiiGIVHUQk1VBjF8MYrQyiMvngxgqF7uLQfSWBJVa
xkRRAO2FAcTy/WjJa0Y4V5T4+Wg9vRQsj7JDaacgWQmTicJURZguKHVIdTIw
VO2VsO5aZFjCUM6vhXLJPmyU5WPsyPGjOtV87uFXqJqHqZrVOP5/9lbNL6mI
qb/TcB+2+pPYYLi74aC3dzAhdph5XEuy8rwN6/V2rNsdWLOL+lJRZyrqTRux
XNvEQE7UooqaVA9mjaJG1Ysr1T6kqnyYvOjD+AUfxs77MFLejOGyZgyca1bq
XbuLfUie9qGzwIe2fC9a8zyI5HrofkWlqKgYFZWjooLUSaFqoJUXwZ+N/qaO
8UvttjBpLvlgRswWEzHb1BpDOyaJhcc5qxa8oNPJC4pOPo4fPnAch7/QgJxv
Dak6WbZbJ9UfJMAfKW+/jc1eRpR1BmVTfd1m4ruZny1Ys1qVDfdVqw0rtaIa
uR5LZgcWTQ1YMDkxbxCVy6KCuRFXqpowfbEJUxfcmDjvxliFGyNlavXzYIlb
qYbuKXKjq9CNzlNNSrV0LLcJ0dxGxsMuehknIxQnBclBCyXKje20UDYKkpX2
3EK9MytlyVp4oGVKIhirZTDWJYKxiWWGH7Scx17AgsGEn+/SO7Oid098vAgv
/XUYJ39AvTuxdo3eqU/hyR8S+e/Y7GF0TI+1ajFgxWLEitmEZbOZUWstlkxW
LBrr2JWoVBcV66JyvR5XL6mV7KlKUdnuxESFE2NlToyUOjFcIqrgXegvdqH3
tAvdBU5Gtk505DkRz21ACylMCAKEwKcUm4uicxsjXVGEXkszZKHEqBGvKFJv
D4n9gNQtulWsRqqTEmiNLSh1ejPPPo1UhQGHdunUCTx6/2H8x59W4djfJ1D4
ZHad+jsVDfWHTD6A9a4kVr1eLNYYSepDB/PVZsxVWTB7qVYphLl6oY7xex27
tGGy3Eb3YcdYqR0jJXalcGawWBTR1CsPPHTn1yN50k7u7Yzt7QxV7Qw5bNSX
OuqLlfoinjEQzxpYGBWJZw9MdFlGeGqoY6GkVsgiFrO0QGlrZJZ+w4HLT/xE
Wbr59wz1+Eeqxw9+/xgOf96xp3p8VeVa/SWcX8daglx3dmFWBtxXK024ct6M
6QrxbIl4xkQ8a1KLsZJaBuBWXC62Mgi3YqDAir5TVvSctCrPq3TmWpXnV8Rz
LFFSiFz5SV5y1iSDc1HP1cs57Q4ltV9PqutMK1FMSixMDdFs1huxWlaI5fYE
fqTjrprNfOO+JynhpynhoV0SXqJJOOXiYZW7L8hcTMTKS21JLMXVJ3zGz6pP
/YwUmTB82oTBQhP6C0zozVefEhJPC3WSxNNDMVILSTxVFCCJp4yGSIOh5N1y
9I39aRjJwUxqi5noGLZM57FQkIf5aAKP6daCVIl8lhJ5iRLZSYmcySqR31BH
/xV1mX1xJPUbumUVZ29adcDdQ9i8UITpY68oiwLfzhCCIw+fohAcpRDU49VH
BlH8nBCCrV1CoP6wyfbvSQhHFRhKY1lZwe/B+tlcrDVYv5kxAV+/76f4948X
4siXgjj5/YmsE/Co2vKhDOMeH5GL6PE2rOS9TCdklo+Vaggdxz/f/3Mc+tOL
OPZ3ROinRKhgE2UZCH1Xbf7prAZySVlwnD9ymPpqkOqtjl5UFb5IXL5PXJ79
vB2vflPgsngNLurvqOCZDIszSMY2x9SF5NShQ7T9BvkAqNq8MFhfv5fgfOwU
jvxVQAXnOMGp2g2OWoggn+BXVbuQqj2mqPYM0ox11i6coT4mpMNXwTn88DEV
nD+5QHA6UKCBU7sbnB+rzT+n061yRbfEQuYUtuqqsXz2lIh1ZeCqBycf3/u9
l/Hsn9kIzgCKDy+iNAOcn6rNv6ATflOCwj8tlvpGsWUoxVyuGhR+9hpwfoIf
fywfL/6VHye/N54VHFVktn+3IUCJcfWlsSCCp+4BbFQUYPzI8T/NwOXb9z9D
XCpx9H+3E5erWXFRf88Av8gm7R1dWCt6BSv22o/rxnxegeQkIfkFnv1cHV79
xgCKBCSluyE5rLb8sk7aL8ZVltaUxZMYlk4cUfz1A7rmLymQPIF/++hJQtKs
QnJsNyRBQqJCjWMZ0t4tpF2sNIQY4x5+BtPnDXIB7Zbt/Zp/Ii5P/UkFjv5t
Gwp+osPFv4PLS2rzJ2TzvZT2Jk3aZUo+8pMnFU/zQbmWru6NnMRjv/cSfv45
K4HpzwaMOmL1B3cyBV1sPzRYsFpepFh5zbsnGRKKyvIh0ihpnDRJohWUyxS3
bO9FfE3BLg8v/qUPeY+NEbvVa7BT8xL5qyyqMlS0pjEtlGFgAlu1F7FYeBLz
LQlZV/NGJSptUJ4GNDAeURcGgnKxICE2pHINGCQNnzQgxchZLCzcqcPdSI/+
rQ/9DE99slzF/Qnifupa3NUx4aROkcxUpFktu6o6gyvHX0GK3F4nh6ZVb5XL
tS3r9vCMiJO68owYyKdjO23ExBkOr8TAe6+UihXD2XqGeOMpOVD9PObhsY8c
wc8/W4tXvs55fHYBped2zeMpdaSn5Eh36SQTmvVzJ7HmrFNGp24o3EEXZGQg
pz4fKp4LFc+DiudAB07bMVqiPuspnvGcZ1glQq0lWU88Q5O9Nikr4bUs3Kju
onKqXxBT/d3sU62udcjfuLkL56jfQer3ypRYAO7ESv4xBJ7+hZKAVDP4sjE5
EclKJM+tPNTXd8arPMQ3eZHJUI0Hi+YmJksuJk/1TKZEQmXksQGzJoNMWnYW
hr/1oafx5CfL8PLfxDndV7JO9xl1eOoo36saibY02jUjwaC+5Scv4OK/FBI0
P/x5IcQKI8pTUEOV6lNP4mkn8YTTSkMEG41BbLmbseV1I+1zKMsyG00GXGlt
V8JfZevyIy9yWi145Wt9O9PawGlNqHZL/QEq9RedtHoP6xOU+r87IcstxOvN
GDF0YcrahTlHF5YbmS14O7EZaKflaaWvDWPCZELHyXw8/3mx152LF/7Cq87S
Uc7Spd2zVKHrE+9A6UvnlE5yXooo7xdeTsBzrB3R5xzK596iKMZsw7jiHMaC
Zwir/kFshvsRP2VhGmqG8Tu5aHi8BM1Pm2B5vgPn8q9F/oLaZZEUjFZ3CGUv
XkTZL5xKF4bcHlx6dRA5h4fhLVHrKoYoFCPH6zBVnVAyW+chMyzfuYjGH7sQ
eDGGltNJxKu7EHe0IWb3o87RiVFGHVW7u4q5oyh7nsnuEQ8qXmpRmm44O44O
1wy8HrX2SzVF43AcqiJwHhifmYT99DwclkXeP4mJjlG0NvagrSaKeJUb8UoD
xmPRzK7ijTGUHWZ+/bwf519sU5q+EhtFepIRxVW1rk4zMf5mhtuWHZCSQ+u4
3LEE7/kxRKt70WaJI25qJocWjLe3XsOVM6Y0V/ZsEBPBdrXwpX8Q6Z7LtAdj
yueN2BSWwinYfYuYaJrGUH0KXeZJtF4ah690FPWFQ6jO7YO/ph3xugDiZiYX
nRld3YixvnGUPlUH7ZVOxpFOJOiwe6R03oa14DAWfZdxpXEEY/YRDJgvI1E1
hGj5ALxFfbDldTH6aFcKQpqNzYhZazE7OValF8S3ou4lE3XQRmfdRBJPAUUp
2+wt1IGN5iRW3N2Yc/YgZevBiKkHfZe60VHehXBxAk157TTFrbRUIbkrUo+W
OpPcxtzuYq0/iU2xguq3kuTertev7OuuulqwWN+Kq7WiTjuO4QsxdJfGEC9q
QSAvojzoZ8rxyf1bm7JGMNbXu7uLN2GkO4UFq7oiJKrgxbOzG44GrNrEE8di
pceHqapmjJz3o7/Ej87CZkTzvMzPxLqWi4Zx13b61ZT6uLGgS2pPZ+Rq/kJ7
EvM1Ym/LgBXxsLrRjIUa8cC+eIhf5MNiX8uBrlMOGv965k42ZvlWOgWzyOqz
bdjLGVErTe9UsrZphtgzTOxnymnZeDxVwsytyIjBU2qm1kYXE1J+S8Go+OpL
GRs+6g9G4ZzydhsWWhhVnDJgJF913/256l5Am9wk8cgNkhrZjFzGUzf0oe5q
3Ywl+rkrYpdAxiYJXiWeghO7B13iHL9X6yb+//nN+18d/+r4V8e/Ov5/dIxf
+7+WhWo2\
\>"]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"now", " ", "reverse", " ", "the", " ", "three", " ", "rotations"}],
   " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5716782232378025`*^9, 3.571678237995428*^9}, {
  3.571678285138711*^9, 3.5716783310963917`*^9}, {3.5723331603881426`*^9, 
  3.572333180714978*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"U", "=", "p"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"RotationMatrix", "[", 
    RowBox[{
     RowBox[{"-", "T3"}], ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "0", ",", "0"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"MatrixForm", "[", "T", "]"}], "\n", 
 RowBox[{"p", "=", 
  RowBox[{"U", ".", "T"}]}]}], "Input",
 CellChangeTimes->{{3.572332884330057*^9, 3.5723329259509306`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0.`", "0.`"},
     {"0.`", "0.7154015443235817`", "0.6987135538827298`"},
     {"0.`", 
      RowBox[{"-", "0.6987135538827298`"}], "0.7154015443235817`"}
    },
    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.5723354150913024`*^9, 3.572335447242959*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2.700437574415364`", ",", "0.`", ",", "0.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "0.5280754835921541`"}], ",", 
     RowBox[{"-", "1.251539233788441`"}], ",", 
     RowBox[{"-", "1.2223448953982994`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5723354150913024`*^9, 3.572335447242959*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"U", "=", "p"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"RotationMatrix", "[", 
    RowBox[{
     RowBox[{"-", "T2"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"MatrixForm", "[", "T", "]"}], "\n", 
 RowBox[{"p", "=", 
  RowBox[{"U", ".", "T"}]}]}], "Input",
 CellChangeTimes->{{3.572332884330057*^9, 3.5723329607389917`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.32994867848144904`", 
      RowBox[{"-", "0.9439988715927289`"}], "0.`"},
     {"0.9439988715927289`", "0.32994867848144904`", "0.`"},
     {"0.`", "0.`", "1.`"}
    },
    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.572335415106902*^9, 3.572335447242959*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0.8910058089999991`", ",", 
     RowBox[{"-", "2.5492100230547097`"}], ",", "0.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.3556894324000002`"}], ",", "0.08555894437059663`", ",", 
     RowBox[{"-", "1.2223448953982994`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.572335415106902*^9, 3.572335447242959*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"U", "=", "p"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"RotationMatrix", "[", 
    RowBox[{
     RowBox[{"-", "T1"}], ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "0", ",", "0"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"MatrixForm", "[", "T", "]"}], "\n", 
 RowBox[{"p", "=", 
  RowBox[{"U", ".", "T"}]}]}], "Input",
 CellChangeTimes->{{3.572332884330057*^9, 3.5723329812218275`*^9}, {
  3.572333014730686*^9, 3.57233303396552*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0.`", "0.`"},
     {"0.`", "0.4946389506381359`", 
      RowBox[{"-", "0.8690985608730484`"}]},
     {"0.`", "0.8690985608730484`", "0.4946389506381359`"}
    },
    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.5723354151225023`*^9, 3.572335447258559*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0.8910058089999991`", ",", 
     RowBox[{"-", "1.2609385707599998`"}], ",", "2.2155147623999984`"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1.3556894324000002`"}], ",", 
     RowBox[{"-", "1.0200174030200002`"}], ",", 
     RowBox[{"-", "0.6789785517999996`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5723354151225023`*^9, 3.572335447258559*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Plota", " ", "=", " ", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"p", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c1", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"p", "[", 
         RowBox[{"[", "2", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c3", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{"v2", ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c2", "]"}], "]"}]}], "]"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Plotc", "=", 
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", "v2"}], "}"}], ",", "0.02"}],
         "]"}], ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"p", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.02"}], "]"}], 
       ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"p", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", "v2"}], "}"}], ",", "0.02"}],
         "]"}]}], "}"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb0", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c3", ",", 
      RowBox[{"p", "[", 
       RowBox[{"[", "2", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb1", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c1", ",", 
      RowBox[{"p", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb2", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c2", ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
  "Plota", ",", "Plotb0", ",", "Plotb1", ",", "Plotb2", ",", "Plotc", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Axes", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"xl", ",", "yl", ",", "zl"}], "}"}]}], ",", 
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "\[Rule]", "False"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "Automatic"}], ",", " ", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.572332884330057*^9, 3.5723329812218275`*^9}, {
  3.572333014730686*^9, 3.5723330352603226`*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.5], 
    SphereBox[{0.8910058089999991, -1.2609385707599998`, 2.2155147623999984`},
      0.255], SphereBox[{-1.3556894324000002`, -1.0200174030200002`, \
-0.6789785517999996}, 0.5241], SphereBox[{0, 0, 0}, 0.178]}, 
   Text3DBox[
    "51", {-1.3556894324000002`, -1.0200174030200002`, -0.6789785517999996}, {
    1, 0}], Text3DBox["20", {0.8910058089999991, -1.2609385707599998`, 
    2.2155147623999984`}, {1, 0}], 
   Text3DBox["21", {0, 0, 0}, {1, 0}], {
    TubeBox[{{0.8910058089999991, -1.2609385707599998`, 
      2.2155147623999984`}, {0, 0, 0}}, 0.02], 
    TubeBox[{{0.8910058089999991, -1.2609385707599998`, 
     2.2155147623999984`}, {-1.3556894324000002`, -1.0200174030200002`, \
-0.6789785517999996}}, 0.02], 
    TubeBox[{{-1.3556894324000002`, -1.0200174030200002`, \
-0.6789785517999996}, {0, 0, 0}}, 0.02]}},
  AspectRatio->Automatic,
  Axes->Automatic,
  AxesLabel->{
    FormBox["xl", TraditionalForm], 
    FormBox["yl", TraditionalForm], 
    FormBox["zl", TraditionalForm]},
  Boxed->False,
  PlotRange->Automatic,
  SphericalRegion->True]], "Output",
 CellChangeTimes->{3.5723354151381025`*^9, 3.572335447274159*^9}]
}, Open  ]]
},
WindowSize->{798, 857},
WindowMargins->{{484, Automatic}, {Automatic, 27}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
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, 2808, 57, 152, "Input"],
Cell[CellGroupData[{
Cell[3390, 81, 5674, 139, 412, "Input"],
Cell[9067, 222, 1160, 27, 375, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10264, 254, 1779, 40, 132, "Input"],
Cell[12046, 296, 130, 2, 31, "Output"],
Cell[12179, 300, 789, 20, 72, "Output"],
Cell[12971, 322, 415, 10, 31, "Output"]
}, Open  ]],
Cell[13401, 335, 817, 14, 31, "Input"],
Cell[CellGroupData[{
Cell[14243, 353, 2226, 43, 112, "Input"],
Cell[16472, 398, 127, 2, 31, "Output"],
Cell[16602, 402, 787, 20, 72, "Output"],
Cell[17392, 424, 388, 10, 31, "Output"]
}, Open  ]],
Cell[17795, 437, 721, 13, 31, "Input"],
Cell[CellGroupData[{
Cell[18541, 454, 1565, 33, 112, "Input"],
Cell[20109, 489, 111, 1, 31, "Output"],
Cell[20223, 492, 787, 20, 72, "Output"],
Cell[21013, 514, 355, 10, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21405, 529, 3496, 102, 232, "Input"],
Cell[24904, 633, 16391, 278, 345, 883, 22, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]],
Cell[41310, 914, 309, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[41644, 924, 460, 13, 92, "Input"],
Cell[42107, 939, 787, 20, 72, "Output"],
Cell[42897, 961, 394, 10, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43328, 976, 460, 13, 92, "Input"],
Cell[43791, 991, 787, 20, 72, "Output"],
Cell[44581, 1013, 411, 10, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45029, 1028, 508, 14, 92, "Input"],
Cell[45540, 1044, 787, 20, 72, "Output"],
Cell[46330, 1066, 454, 12, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[46821, 1083, 3269, 99, 232, "Input"],
Cell[50093, 1184, 1187, 27, 407, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)