(* 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[     59153,       1451]
NotebookOptionsPosition[     57124,       1379]
NotebookOutlinePosition[     57466,       1394]
CellTagsIndexPosition[     57423,       1391]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "Rotation", " ", "of", " ", "general", " ", "triangle", " ", "into", " ", 
     "xy", " ", "plane", " ", "with", " ", "one", " ", "point"}], ",", " ", 
    RowBox[{"O", " ", "atom"}], ",", " ", 
    RowBox[{"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", " ", "7.3", 
       ".2013"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "T1", ",", "T2", ",", "U", ",", "V", ",", "T", ",", "P", ",", "p", ",", 
      "q", ",", "r"}], "]"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"rad", "=", "57.2957795"}], ";"}], "\[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`"}]}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"q", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", "P", "}"}], ",", 
       RowBox[{"{", "3", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"r", "=", "q"}], ";", " ", 
    RowBox[{"s", "=", "q"}], ";"}]}]}]], "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}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"rotate", " ", "p1", " ", "around", " ", "x"}], "-", 
    RowBox[{"axis", " ", "until", " ", "the", " ", "z"}], "-", "comp"}], " ", 
   "=", " ", "0"}], " ", "*)"}]], "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}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"T1", "=", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"ArcTan", "[", 
     RowBox[{
      RowBox[{"p", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "3"}], "]"}], "]"}], "/", 
      RowBox[{"p", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "]"}]}], "*", "rad", 
   "      ", 
   RowBox[{"(*", 
    RowBox[{"60.35405486", ";"}], "*)"}]}]}]], "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.571673026353874*^9}, {
   3.5716738009420347`*^9, 3.5716738278676825`*^9}}],

Cell[BoxData["60.354054833041914`"], "Output",
 CellChangeTimes->{
  3.5716730966787977`*^9, 3.5716732718827057`*^9, 3.571673371613681*^9, {
   3.5716734123141527`*^9, 3.5716734688174515`*^9}, 3.571673525570351*^9, {
   3.571673814779259*^9, 3.57167383253209*^9}, 3.5716740197636194`*^9, 
   3.5716740620708933`*^9, 3.5716758611596537`*^9, 3.5716779926693974`*^9, 
   3.571678065271925*^9, 3.571681214301321*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"U", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "P", "}"}], ",", 
     RowBox[{"{", "4", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"V", "=", "U"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"P", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"p", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"p", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"p", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "4"}], "]"}], "]"}], "=", "1"}]}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "U", "]"}]}], "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}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.8910058089999993`", 
      RowBox[{"-", "1.2609385707599998`"}], "2.215514762399999`", "1"},
     {
      RowBox[{"-", "1.3556894324000002`"}], 
      RowBox[{"-", "1.02001740302`"}], 
      RowBox[{"-", "0.6789785517999993`"}], "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.571673096694398*^9, 3.571673271898306*^9, 3.571673371613681*^9, {
   3.5716734123297524`*^9, 3.5716734688330517`*^9}, 3.571673525585951*^9, {
   3.571673814794859*^9, 3.5716738325476904`*^9}, 3.571674019779219*^9, 
   3.5716740620864935`*^9, 3.5716758611596537`*^9, 3.5716779926693974`*^9, 
   3.5716780652875247`*^9, 3.571681214316921*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"T2", "=", 
   RowBox[{"T1", "/", "rad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "4", "}"}], ",", 
     RowBox[{"{", "4", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", "1.0"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"4", ",", "4"}], "]"}], "]"}], "=", 
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "2"}], "]"}], "]"}], "=", 
   RowBox[{"Cos", "[", "T2", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"3", ",", "3"}], "]"}], "]"}], "=", 
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "3"}], "]"}], "]"}], "=", 
   RowBox[{"Sin", "[", "T2", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"3", ",", "2"}], "]"}], "]"}], "=", 
   RowBox[{"-", 
    RowBox[{"T", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "3"}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "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.571652801459519*^9}, {
   3.5716528498352036`*^9, 3.571652858399619*^9}, {3.5716529333577504`*^9, 
   3.5716529342001514`*^9}, {3.571652969768214*^9, 3.5716532517387094`*^9}, {
   3.5716533648857083`*^9, 3.571653366430111*^9}, 3.571657597719143*^9, {
   3.571657731333378*^9, 3.571657878441636*^9}, {3.5716591457098618`*^9, 
   3.5716591461310625`*^9}, {3.5716592919757185`*^9, 3.571659294970924*^9}, {
   3.5716666093316755`*^9, 3.5716666223264985`*^9}, {3.5716673183873205`*^9, 
   3.5716673186837215`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0", "0", "0"},
     {"0", "0.4946389506381359`", "0.8690985608730484`", "0"},
     {"0", 
      RowBox[{"-", "0.8690985608730484`"}], "0.4946389506381359`", "0"},
     {"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.571673096709998*^9, 3.5716732719139056`*^9, 3.571673371629281*^9, {
   3.571673412345352*^9, 3.571673468848652*^9}, 3.571673525601551*^9, {
   3.571673814810459*^9, 3.5716738325632906`*^9}, 3.571674019794819*^9, 
   3.5716740621020937`*^9, 3.571675861175254*^9, 3.571677992684997*^9, 
   3.5716780652875247`*^9, 3.5716812143325214`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"V", "=", 
   RowBox[{"U", ".", "T"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "<", 
     RowBox[{"P", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"q", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"q", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"q", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", "q"}], "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.57165915589668*^9, 3.5716591635562935`*^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.5716730967255983`*^9, 3.5716732719295053`*^9, 3.5716733716604805`*^9, {
   3.5716734123609524`*^9, 3.5716734688642516`*^9}, 3.5716735256171513`*^9, {
   3.5716738148260593`*^9, 3.5716738325632906`*^9}, 3.571674019810419*^9, 
   3.571674062117694*^9, 3.5716758612064533`*^9, 3.5716779927005973`*^9, 
   3.571678065303125*^9, 3.5716812143325214`*^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[{
  RowBox[{"Clear", "[", 
   RowBox[{"T1", ",", "T2", ",", "U", ",", "V", ",", "T"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"T1", "=", 
  RowBox[{"180", "+", 
   RowBox[{
    RowBox[{"ArcTan", "[", 
     RowBox[{
      RowBox[{"q", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", 
      RowBox[{"q", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1"}], "]"}], "]"}]}], "]"}], "*", "rad", "     ", 
    RowBox[{"(*", " ", 
     RowBox[{"109.26548", ";"}], "*)"}]}]}]}]}], "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}}],

Cell[BoxData["109.26566053108144`"], "Output",
 CellChangeTimes->{
  3.571673096741198*^9, 3.5716732719451056`*^9, 3.5716733716760807`*^9, {
   3.5716734123765526`*^9, 3.5716734688798513`*^9}, 3.5716735256327515`*^9, {
   3.5716738148260593`*^9, 3.5716738325788903`*^9}, 3.5716740198260193`*^9, 
   3.5716740621332936`*^9, 3.5716758612220535`*^9, 3.571677992716197*^9, 
   3.5716780653187246`*^9, 3.571681214348121*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"U", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "P", "}"}], ",", 
     RowBox[{"{", "4", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"V", "=", "U"}], ";"}]}], "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[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"P", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"q", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"q", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"q", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "4"}], "]"}], "]"}], "=", "1"}]}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "U", "]"}]}], "Input",
 CellChangeTimes->{{3.5716591674719*^9, 3.571659171777508*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.8910058089999993`", 
      RowBox[{"-", "2.5492100230547097`"}], "0.`", "1"},
     {
      RowBox[{"-", "1.3556894324000002`"}], "0.08555894437059663`", 
      RowBox[{"-", "1.2223448953982992`"}], "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.571673096756798*^9, 3.5716732719607058`*^9, 3.571673371691681*^9, {
   3.5716734123921523`*^9, 3.5716734688954515`*^9}, 3.5716735256483517`*^9, {
   3.5716738148416595`*^9, 3.57167383259449*^9}, 3.5716740198416195`*^9, 
   3.5716740621488934`*^9, 3.5716758612376537`*^9, 3.571677992731797*^9, 
   3.571678065334325*^9, 3.571681214363721*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"T2", "=", 
   RowBox[{"T1", "/", "rad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "4", "}"}], ",", 
     RowBox[{"{", "4", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", 
   RowBox[{"Cos", "[", "T2", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "2"}], "]"}], "]"}], "=", 
   RowBox[{"-", 
    RowBox[{"Sin", "[", "T2", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
   RowBox[{"-", 
    RowBox[{"T", "[", 
     RowBox[{"[", 
      RowBox[{"1", ",", "2"}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"4", ",", "4"}], "]"}], "]"}], "=", "1"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "2"}], "]"}], "]"}], "=", 
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"3", ",", "3"}], "]"}], "]"}], "=", 
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"4", ",", "4"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "T", "]"}]}], "Input",
 CellChangeTimes->{{3.571657627093994*^9, 3.5716576537544413`*^9}, 
   3.571657963227785*^9, {3.57165800620586*^9, 3.5716581329092827`*^9}, {
   3.5716581674633436`*^9, 3.5716581993029995`*^9}, {3.571658274760332*^9, 
   3.571658296990371*^9}, {3.571659302536937*^9, 3.5716593120217543`*^9}, {
   3.571666638082526*^9, 3.5716666389249277`*^9}, {3.5716672827880583`*^9, 
   3.5716672849876623`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "0.32994867915859694`"}], 
      RowBox[{"-", "0.9439988713560505`"}], "0", "0"},
     {"0.9439988713560505`", 
      RowBox[{"-", "0.32994867915859694`"}], "0", "0"},
     {"0", "0", "1", "0"},
     {"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.571673096772398*^9, 3.571673271976306*^9, 3.571673371707281*^9, {
   3.5716734124077525`*^9, 3.571673468926652*^9}, 3.5716735256483517`*^9, {
   3.5716738148572593`*^9, 3.5716738326100903`*^9}, 3.5716740198572197`*^9, 
   3.5716740621644936`*^9, 3.571675861253254*^9, 3.571677992731797*^9, 
   3.571678065349925*^9, 3.571681214363721*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"V", "=", 
   RowBox[{"U", ".", "T"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "<", 
     RowBox[{"P", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"r", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"r", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"r", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", "r"}], "Input",
 CellChangeTimes->{{3.571657627093994*^9, 3.5716576537544413`*^9}, 
   3.571657963227785*^9, {3.57165800620586*^9, 3.5716581329092827`*^9}, {
   3.5716581674633436`*^9, 3.5716581993029995`*^9}, {3.571658274760332*^9, 
   3.571658296990371*^9}, {3.571659187424335*^9, 3.5716591941635466`*^9}, 
   3.5716730759307613`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.700437574415364`"}], ",", "1.9370738346680128`*^-9", ",",
      "0.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.5280754844899064`", ",", "1.2515392334096427`", ",", 
     RowBox[{"-", "1.2223448953982992`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.571673096787998*^9, 3.5716732719919057`*^9, 3.571673371707281*^9, {
   3.5716734124077525`*^9, 3.5716734689422517`*^9}, 3.5716735256639514`*^9, {
   3.571673814872859*^9, 3.5716738326256905`*^9}, 3.5716740198728194`*^9, 
   3.571674062180094*^9, 3.5716758612688537`*^9, 3.5716779927473974`*^9, 
   3.571678065349925*^9, 3.571681214379321*^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[{
  RowBox[{"Clear", "[", 
   RowBox[{"T1", ",", "T2", ",", "U", ",", "V", ",", "T"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"T1", "=", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"ArcTan", "[", 
     RowBox[{
      RowBox[{"r", "[", 
       RowBox[{"[", 
        RowBox[{"2", ",", "3"}], "]"}], "]"}], "/", 
      RowBox[{"r", "[", 
       RowBox[{"[", 
        RowBox[{"2", ",", "2"}], "]"}], "]"}]}], "]"}]}], "*", "rad", "  ", 
   RowBox[{"(*", 
    RowBox[{"44.323845", ";"}], "*)"}]}]}]}], "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}],

Cell[BoxData["44.32388315350226`"], "Output",
 CellChangeTimes->{
  3.571673096803598*^9, 3.571673272023106*^9, 3.5716733717228813`*^9, {
   3.5716734124389524`*^9, 3.571673468957852*^9}, 3.5716735256951513`*^9, {
   3.571673814888459*^9, 3.5716738326412907`*^9}, 3.571674019888419*^9, 
   3.5716740621956935`*^9, 3.5716758612688537`*^9, 3.5716779927629976`*^9, 
   3.571678065365525*^9, 3.5716812143949213`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"U", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "P", "}"}], ",", 
     RowBox[{"{", "4", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"V", "=", "U"}], ";"}]}], "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[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"P", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"r", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"r", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"r", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"U", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "4"}], "]"}], "]"}], "=", "1"}]}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "U", "]"}]}], "Input",
 CellChangeTimes->{{3.571660729689844*^9, 3.5716607386286592`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "2.700437574415364`"}], "1.9370738346680128`*^-9", "0.`", 
      "1"},
     {"0.5280754844899064`", "1.2515392334096427`", 
      RowBox[{"-", "1.2223448953982992`"}], "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.571673096819198*^9, 3.571673272038706*^9, 3.571673371738481*^9, {
   3.5716734124389524`*^9, 3.5716734689734516`*^9}, 3.5716735257107515`*^9, {
   3.5716738149040594`*^9, 3.571673832656891*^9}, 3.5716740199040194`*^9, 
   3.5716740622112937`*^9, 3.5716758612844534`*^9, 3.571677992778598*^9, 
   3.571678065381125*^9, 3.5716812143949213`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"T2", "=", 
   RowBox[{"T1", "/", "rad"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "4", "}"}], ",", 
     RowBox[{"{", "4", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", "1.0"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"4", ",", "4"}], "]"}], "]"}], "=", 
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"1", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "2"}], "]"}], "]"}], "=", 
   RowBox[{"Cos", "[", "T2", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"3", ",", "3"}], "]"}], "]"}], "=", 
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "3"}], "]"}], "]"}], "=", 
   RowBox[{"Sin", "[", "T2", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T", "[", 
    RowBox[{"[", 
     RowBox[{"3", ",", "2"}], "]"}], "]"}], "=", 
   RowBox[{"-", 
    RowBox[{"T", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "3"}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "T", "]"}]}], "Input",
 CellChangeTimes->{{3.571666647941743*^9, 3.571666651064948*^9}, {
  3.5716673391977572`*^9, 3.571667339494158*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1.`", "0", "0", "0"},
     {"0", "0.7154015442178727`", "0.6987135539909635`", "0"},
     {"0", 
      RowBox[{"-", "0.6987135539909635`"}], "0.7154015442178727`", "0"},
     {"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.571673096819198*^9, 3.571673272054306*^9, 3.5716733717540812`*^9, {
   3.5716734124545527`*^9, 3.571673468989052*^9}, 3.5716735257263517`*^9, {
   3.5716738149196596`*^9, 3.5716738326724906`*^9}, 3.5716740199040194`*^9, 
   3.571674062226894*^9, 3.5716758613000536`*^9, 3.571677992778598*^9, 
   3.5716780653967247`*^9, 3.5716812144105215`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"V", "=", 
   RowBox[{"U", ".", "T"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "<", 
     RowBox[{"P", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"s", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"s", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"s", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}], "=", 
      RowBox[{"V", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", "s", "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.571660822463207*^9, 3.5716608285472174`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2.700437574415364`"}], ",", "1.3857856125855326`*^-9", ",",
      "1.3534597433637913`*^-9"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "0.5280754844899064`", ",", "1.7494220462969692`", ",", 
     "1.1102230246251565`*^-16"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.571673096834798*^9, 3.571673272069906*^9, 3.571673371769681*^9, {
   3.571673412470153*^9, 3.571673468989052*^9}, 3.571673525741952*^9, {
   3.57167381493526*^9, 3.5716738326880903`*^9}, 3.5716740199196196`*^9, 
   3.5716740622424936`*^9, 3.571675861315654*^9, 3.5716779927941976`*^9, 
   3.571678065412325*^9, 3.571681214426121*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"c1", "=", "20"}], ";", " ", 
  RowBox[{"c2", "=", "21"}], ";", " ", 
  RowBox[{"c3", "=", "51"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"v2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0", ",", "1"}], "}"}]}], 
  ";"}], "\[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[{"s", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c1", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"s", "[", 
         RowBox[{"[", "2", "]"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c3", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"Sphere", "[", 
       RowBox[{
        RowBox[{"Take", "[", 
         RowBox[{"v2", ",", "3"}], "]"}], ",", 
        RowBox[{"u2", "[", 
         RowBox[{"[", "c2", "]"}], "]"}]}], "]"}]}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Plotc", "=", 
    RowBox[{"Graphics3D", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"s", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"Take", "[", 
            RowBox[{"v2", ",", "3"}], "]"}]}], "}"}], ",", "0.02"}], "]"}], 
       ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"s", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"s", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.02"}], "]"}], 
       ",", 
       RowBox[{"Tube", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"s", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"Take", "[", 
            RowBox[{"v2", ",", "3"}], "]"}]}], "}"}], ",", "0.02"}], "]"}]}], 
      "}"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb0", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c3", ",", 
      RowBox[{"s", "[", 
       RowBox[{"[", "2", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plotb1", "=", 
   RowBox[{"Graphics3D", "[", 
    RowBox[{"Text", "[", 
     RowBox[{"c1", ",", 
      RowBox[{"s", "[", 
       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}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.5], 
    SphereBox[{-2.700437574415364, 1.3857856125855326`*^-9, 
     1.3534597433637913`*^-9}, 0.255], 
    SphereBox[{0.5280754844899064, 1.7494220462969692`, 
     1.1102230246251565`*^-16}, 0.5241], SphereBox[{0, 0, 0}, 0.178]}, 
   Text3DBox["51", {0.5280754844899064, 1.7494220462969692`, 
    1.1102230246251565`*^-16}, {1, 0}], 
   Text3DBox["20", {-2.700437574415364, 1.3857856125855326`*^-9, 
    1.3534597433637913`*^-9}, {1, 0}], 
   Text3DBox["21", {0, 0, 0}, {1, 0}], {
    TubeBox[{{-2.700437574415364, 1.3857856125855326`*^-9, 
      1.3534597433637913`*^-9}, {0, 0, 0}}, 0.02], 
    TubeBox[{{-2.700437574415364, 1.3857856125855326`*^-9, 
     1.3534597433637913`*^-9}, {0.5280754844899064, 1.7494220462969692`, 
     1.1102230246251565`*^-16}}, 0.02], 
    TubeBox[{{0.5280754844899064, 1.7494220462969692`, 
      1.1102230246251565`*^-16}, {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.571678065427925*^9, 
  3.5716812146757216`*^9},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztfHd4ZVd17w/b44JNM810UxInD3gISEgIIaSSF/KS5xAehBJ4DmSwDcbG
BoMxGOPee5FmJI2k0UhX0qjM3Kp7Jd3eq3qvozZV06UZlfN+e599z70zPhqc
vO97X/7wfLM0R+eeu/daa6/yW3uvM/9ww09u/PdbbvjJTd+94dq/vf2G2268
6bt3XPvFW2/nrYtfB7zuRtKXroW41oD8D41/tZvED3HntevXrv8/X79V/LhY
2OI/F919m7yWt+WVoC//F/pccn2R+PQr/0F5328y8v8o+vy9xsh8SN5WDxY/
9C5DaRebDPe3RU++w3jyMpMnv2i6EFeaPPk3RU++2XjyIpMnr/gt336D8e0t
Jk9etpn4xUNcaTx0uckQZuwXhv2TonEuN1h5vcl3rjK5p7N8yTn3+OezRYNu
MQZ9739ugIsuqN73mNy75IKDblxQ4WZMXnzB8VaN8a77D383Z6zDdcbVRzcb
pXi9subDflhcyLsfM77x5xf8xrUXdJ3rjE/PM78vXHDQ9xlf+8SmKjn33ucv
ON67jfE+a/LdP9Ufe8HMVfTf+OfT5iO/0xj5Tzcf+cXzPtrMFj9tcPCOzdz2
8yZf+5z+xEvKIYo/yn/1SxdcJbNB/+Q/xvrfGBN8cfPBijS85RwN56//2hjl
SyajKJvfZvLRF4wv/k/jSk1aZkz6Z8ZHZhL8nfkCq1EGjVHMnP4fTO79kf7F
PmPSfzF56DP6Q93GQ9/c/KGM4ZvfNnnoD/WHkhdc6H81ufcH+hdjxhcvNXnI
bEbdIxA2vmiWqS7wxaCh0h9sZu3qSZ8xxXmfv0X+vBINpVbsIblJAVKclCX1
kYZJY2VWTAri9ZdNOPqkPmqnMc8lJg+9Uf68ChYOIiZsK5owQcqRBsVk26yY
3m7F7A4rFqusF+FAjZVD7a+yGkP94+Y8vGjw8CHj6p1q/vzMLSRnqQ2+Ugdi
pU7O7MLAtnaMlbsxXeXGwq52HLQ4caTJgWO77RfhRIvg4NhuKxZ2WjfloKQ4
jpz70fsNXt6lr5Vg5WrsJhM2Tt5Z6kG41It0mR9924IYrQphelcI8w0hHGwO
YaktgBM2H047OrHidOOMy44zTitO28lRvfV/mcyoB348ZRi92ZK8ZzPDuKbA
5Vv504nW0k64SgPwl4ZpHXF0VyQxVJPGRH0as7szOLAngyV7GidcKSx7Ejjb
FcWaL4QNvw/rfjc2fFac6bDiAqya6e1dhmWbsf+OApdv508PGkv92FsaRUdp
EqHtWaQqe9C3sxejDf2YbhnAgnUQh1wDON4xgNPefpwN9mItTBAQZXCIx0kh
Uhe0iB0rXnNu//vm3BbArVkMkJ/SFgW7ltIOUgjNZNRemkNXaR+iFYPI1o5g
oHEU423j2GefwH73JJa6JnAyMIGVyDjW4qNYTw1DywxAy/VA6ybj3Yw9Wfp4
wo7rTeb9+OYMFzD2B5W6a56skJq0cJ0tpWm6Si+c24bgqxxFfOckuhunMdS2
D5OOOcx55nDQN4+joXmcjs3hbGoWG9l90HqmoPVNQBsYgTZIRoeo4QFqt9eL
7EvleoqRNvmJIp+Vq3y1wdBHFLt1T9VITek2qv8M1OxDomFe/qaHvTdi3ncI
h8KHcTxxGMuZQ1jtOQit/wAnX4A2MgdtbBra+Di0iSH+243eOjd6q2rNFPZu
g5036/zVGVyZJT4ZZ4jY659qgOXZLnLbLz9vK5+W/wbr9yPZcgi99iXxq173
EOqPH4E2SSZnFqHtmyVNkobR2+xHb33LP23mrufde4PO4k6DRbP8qrN4KSxP
tsHyTJBs5tBUNiOfce86pPt+AflcX7IVrSU3wv6prcqSOegidbk4h15bD3r3
RsimFZvzeC4MODfBodrQ78fzfD1Ol3gqAsvTvWh4djwfinRtXyb5EX/u+0wT
nv/GIIaqTtGYerDx9M8Uf5egj57S6+pHrzNO/pxmWaqwrma5XS/dUGlw+wmV
Ii2P2slgDJYn+qHtpyXtPyC/+P2SH+HeP6jDc1/rk7/Xb9PQZNdg7dZw8hgH
nTsBzd2M7r/eqiS5FO27DqGtYpFZaIZ2MopeTxK9jvZ/vmDouwC32w1X0tPf
lbA8RG4fjcuH/73kJvz6kzXyuvLnJ1FXSg5tepSvHtQUV1fh0ZCG4YMc+Ahp
gFDt5V8jYDmA9qp5tJZNk9sxadW9XSn0tptye43Jsm/iL4r1UkPRZvLl65Z7
//FxKtqKp69PI49wmvZqaIxpqBnQ1L03YTlL108exuHIYSz4D+P6B7gQoxqO
i4VYWIbmtWH661sZz0akKJbSjAwrPW672jIxTS15EczA1RVFwUt+WYeQ78V3
aK53l5ThqX+I44mvN2DX8+S3lRTRsLNfw/3JswxNjI+909jI7cNqehan4/M4
Fp7HIf884+q8jK9DbbO4p1ND737ytEQaHoNW8bCMy5bSiIzT6fa9XzZs4AqD
las31/tz5/H7fnyL/N5V8gKe/FIE2289jF3PktdmmkpYQy35fTTHkr07SkqT
mHSyjOnpYawnRnEmOo6TwQkseSdxgC44y6Q13jbBJDYmk9n1Ny9h95CGw0fJ
xH6W6kE3Kv7+VpEEv2bCo9mm0JuNBdB3RfCMIYGOvT+Ir1OCn5U8jce/GMC2
Hx5A7VOUoIkSBDXs6tPwVK8+lMjpIrfLHJ8I6Tk/lsZGJIdVYoFlXz+Odw7g
sGsQi7ZBTLcOSOzQt7MPqYoehEqz8JSmcAfdPLOgYYNSWR6qheXmu/FqxdE3
OPDkBXPKJ5SF/W8K9tOSx/DYX3lRdtMCdj6hoaGBggU01FGoZ/vyWzjXSHgl
YNaGgFt+L9YJv1YJw1Y6EjjZnpLwTMC0ueYMJixpCd+6y5MSzvlpTgLetRDm
Cbj3i/usqHiGyhJRIerHiTu2ftWETTOYo+R77IIJ6WNKvq9QvjtKHsQjf+5B
6fdnsfPRdTTUUz4f5evR8MJA3m6vxoFGHfCecdlwtp0A2NUuAfFJAuOjewI4
1KIDZgGcBYDu2x6UgFoA6w66ipVAu7HUpuPFLfjWrVbE5jSsCsucIEhoKoWZ
jJtvUuFhQ0azfb/fVzJ+mTL+uORePPxnDrz8vWnUPLSGhl2UsUtDPXPFy4P5
vb03Y7HWKssMUXAIOtZsxxKrEFGNiKpkprod46xSBsraZdUiqhdRxYhqhlXN
xUK2ixTb5VkNsyKiHyIlaQe/utk01F1AwAfO87QrZCr+UcndePBze/HSDROo
vv8sLDUUpoOekNNQNpivKt+E/dVWHGS9JOq3RV7PVeqFnagmRaHXrQq/gCoE
21R5pq+Pzto3OXZgn4YzYpWmD0Fr22EqRGGDTt/Qxm9elX8JcX5Ycifu/+MW
vPDtUVT95gwsVRTHTXGov+2Dedx+JaZK9QJ4XBXF/apmjReJIIroRtIViv3n
ExqmDnNOQdk0Vh+6DWbpUwdo+LoRyj9kyHMh9xHs31TyYwKjRgmMdtyzjPpK
st/OMJHWUD6Ur1e24PqXNXRMaTgtNLnvKDR7nV7fyPl0uI1vnDeZ+PMBQ5Of
Lpp2a8kPce+ndfSz4+7TqN/OaZ2cNqWhcii/9bAFjzOJjB5S2KK3l5jtTvy9
IZxZVaeXIrqaXkXl+gcquX2v5Pu4h2Dn2a/kUHkXwU6ZDscakhp2DGp/pALW
9Q9psI9pOCER2kmJ0LJ/sdWsDrhWZ+RbJh8Vti3PAwWfVNzcQBX9qqRC4paK
O0+g7iVyY2VWimsCev2p0s69jAEDB1R6H2Bl8vI9pkXnBVh5m8GKmZ4+rvj5
Nvn5RUmpxCTldxzFrhfIyx4DR/15Xju3HEfLsIYlYSQLK9B8dkx+davZbscH
ikz2lcCvUFnmI8c3yMHPS57XUcZtRBnPnYuKfhM/o+9fUi0/pSHlFjmYYGN4
QqIeMx7eX2S10p7e8oppv8Zp7yx5SocGtxAaPM0pdxfAzSOZdbV1uQV1/P2g
WIn9a9BCHuz/3la9WDWfdtNtt08WZbeflDzK7N2FspsXUPskzbGxAEue7Mlv
fW7BzVRFYl7DuhB5jEV07TNFfqKUfaMhnhlm/ria9p857e0lDzCpuvWk+tgG
Giyc1q+Dhmf68oeHW1DJmD0vRBaJPsZEf/tWdbL3iiVVUv/7K8LBP3G+20p+
jYc/rxLcw0xwdZzPW0jiX1Xz/WuthvCshrNCzMlFmXXN9qpVAvo309hzS8kv
CgnogVU01Bay6UuM2N9Uk71I15/JZ8BUDCu/vFmttVnNiu8YEedTRZP9oOSn
Mj28mE8P1YVsV8rJvpOPr4w3XdMaloVkM4dlpjKb7N1Frvxbl1IwcGPJbbjv
DxsY4IdQ9esVWHao/JTRsI0MfE8x8CQdaTyfa3IZrD54m1pKk8Ky2GnyEXRr
yc0M6bv0kP5LhvRyTuQyMsmNeUEf1eCcUDXu7DFojnp1xGs6z9fOW8LL8V0K
dU9JtYrUp/TC2WGkjlvUPPfRfIbyVWlfn8wcZocjeumpo4Ki1bsc/4fz/LKk
vBCDVfnbmJAZ4cf5mPfjU2gbUaXi/ClonhakvrD1L43xlCD/9CpW7HJZTP2i
5GU90P7kGHa9qOkFKwN/1aB2pxLuZ9RrT76uGxyG9tI9KvqZJBrowW/TncX/
pub+ugyxz+HJvwtj+4+P6CG2jRSVQf5uNXfDoKrHFs8wwDsw/pWt+CuTYXXg
o0ehTQ8orlOW+lUZZ5/U4+yPDqL2Gc7bosI7A8DPfYu4RzHwQ/KUUuWTNjIp
w7ue4d6IJNFTRiGrHqLaftIgabjMhtFtNoxvt2GynMV8hQ37dtgwW2XDfLX1
7Sas6fd0J7y42K31uPyIisuLMi7LcjGk54SH0mu4X/FaxQC2XyzSgXXmhA4s
fnerMoh3EXkL5GdjhWEjErcTzDqQJuWIyntIfWVOInUnhre5yLuLyN2FyQoX
pne4MFPlwmyNE/M7BbJ3Yn+dAwfqHThkseNwgw1HGm1E/lbC6E0FKxxwfNC4
MgO7+WDyZZkX7ld5YQ47H9/Qi8midPSwEvs7LFCi+epofFamI70J5VrYKbaL
InsospcUoMghihtlQRInpUrbkS11o7vMjd7tHgyUezBEGqnwYKyiAxM7OjBV
1YGZGg9md3owv8uDxTqmW4sbBxvcONTYjiNNLiztduFosxPHWhw43mrHiVYb
TrZZ32Aios4afmez4H7O8ougemvJr/BQPl89oorOovz4hNLDy2lVRIlAFA/g
+G1b1ZL8rjxG2iNPcOwswBwsAJzopPw+yh+k/BEWnTFSsqwDmbJO5Ei9ZV3o
396FwQovhiu9GKvyYoI0VePFzE4vZnd5MVfnxYLFi/0NXhxs7MKh3Z040tyJ
pZYOHGvtwPE2D07scePkXhdOWZ04bXOQ7Fi22UhWs36Tt+r6ea8R08yqimtM
NXVLyc/x4J/s0ZPtg6t66aoy+/N0lmeUpr7MAsA3oyq1qQMys+usfFRWdKIs
EuddrdSWldpyUFPt1FQHNeWllgKlHXSkLlpRFx3Ji1SZD1lSjzwTC9CKAhiq
CGB0RwDj1UFM7gxiujaIfXVBzNUHsdAQxP6mAA7uDuBQcwBHWv1YavPj2B4f
ju/14oTVi5O2Tpyyd2DZ4cGKw40VpwtnnE55nnZWbCu4rHopL/OicrMPGzeu
MdT3ySIF3VxyBwFCswQI1fepcrizgEZeVAp6mlF48rBCI+k4lu++Kd8HJPTz
cf60oYnaaaEd7aF2bNSMk1pxl3bSrrpoV17alY925aefBWhXQWS2BZHbHqKf
hdBfHsZgZRgjVRGM1UQwURvB1K4IZuojmG2IYL4xgsXdERxojuBQawSH28JY
2hPG0b0hHLeFcMIexClHAKedfiy7fFhp9+JMexfOujtIHqx62klOrHns/NdG
ytfpr8CKyh8/spkpbppGiytbscH+GwF8vqmAjyjMPXphLpBXWR6QPK6hfVJV
tjNHiLyq9FpIV2yJVKyFim2gYnfT5Fqp2L1UrJ1KdVGpHiq1i0r1U6mh0iBN
MIREWRip7RFkK6LoJvVVxjCwI4bh6jhGa+IYr41jsi6OaUsCs40JzDUlsNAc
x/6WOA62xXF4TxxL1jiO2mI4bo/hhCOKU64ITrdHsOwOY8UdwhlPEKsdAax2
+rDW6cVaVyfJg3VvO9a7nCQbyWq2ufSWzQNeYb80D7jEOcC9n6rVgd2vClsE
CkFWKiD0IOPfSL5U785KAGlo8dNKi05qsZ3O7EEzNdhGDVqpPQe156bmOqk5
b2mYzhyhM0cZ/mI00zjS2+PIVSTQsyOB/qokhqpTGNmZwtiuFCbqU5hqSGNf
Yxpzu9NYaEljsVXsVaZxaG8aR2wpHLWncMyRxAlXEifbEzjtTmDZE8dKRwxn
O6NY7YpgzRuWR85rvgDW/X65B7oR6OS/1KjfxX/tcn90zWs1RVFvOT+rFvoK
84jy3yRyrdKR6y9OFTY9dIRcIx/7HK7/6Qr2jqo9htnjEiHrp74OkosadJM6
qMUuOruPzu6nswfp7CHaYoS2GKUtxujscTp7AvGyJFLbUsiUp9BdmUZvVQYD
NRkM1WYxWpfFuCWLyYYcpptymG3OYb41h8W2HPbvzeGgNYfD9hyWHDkcc2Vx
vD2Lk+4MTnkyWO5MY6UrjbPeFFZ9Saz5E1gPxKm1KNaDEVIIG6EgcY8fWpja
JP7ZCLn5uxNa0EqtmmtSHZt+6rxwqe9C/LJku47If1a0R6Mj/3qFq+6ibfbl
d0X6+rHx1J3yxLyB+mqkvprop63U115amp36aqeuOqgrL3UVKE3S8lL04TQS
2zJIl2eRrcyipyqHvppuDNZ2Y7iuR+24vh0nPD041dmN013dWPHmcMaXw6o/
i7VAhvKnKW+KsiehRaiXSAwb0Qi0aJgUIPlIXSQPP6dOwvTXoLlOdNCgbzz+
luh3Ob4pi4iXVBFxvLB7pFcuDQW3vFoaVdOQ2rKZP83ypRVPf+YnVFKQwS5M
F43QRWN00ThdNEkXTTGTpGlcGYTLsohtzyFZ0X2NWrm1cDfWIznKmYUWy8jj
ES1O+eMJaIkYKaIfnSRoEwkv71P2mIs6sGEtaDXDFUr0j12whskDV7EpnL/3
EaWOf5F1zbN4Il/XPK82r/Saqkm5r+5iNuVm7fJ87EcsvbKLfPSZAJoe2Ivd
332QLhejyyXockm6XIpaydDlsnS5nCpDaBzJHInipyh+KikTtpaOkrj0abpE
iqInO6gCih6zmUl0qWH++sacfiAgn7lM6aWZHDtInaSg2koWxU+vaoUaFT1X
oh2qworZSivmSYs7rLIX6r1KOfoR0RNG4SVP71qNoq9ZTVVQjlMpR3SgiG4F
PymIu+8Lo/qZKJqfjKL1rnpU/cW99K0M82OWUFeM8Xpo2YzcxdayVEhWdJ3Q
FrJURoZ+kKEykk5sxK2Fvfjz9sXOPw4r0s11m31ni1pccTAg1GU31GWjuuws
NOyss+wYKnNgrNyBqUo7a0I7Fmps2F9rw6E6Kw7XW7FksdLGlhpEI5e4Prdy
uFxumt1R8hAe+0vWhD9Y1I8Qm0lhWb+2mqrRrdToJQXkebKFobuRnvatHybw
3OMptD5DlY1OyRI3fyCk9fBeD8NKD22qR5yrUoU5qjBLV0o5sBI2DyFm9y4u
AlqbFoAXF8cLoUSrrN4cDJZOpmgXba6dCNItq7RxVmnT1R7M1bIqqxcVmV6N
iUrseKuDFZidZLsIp/YIXZ5otbJAsxbFfX1D8Mclv2Gh2Y7SG+fkqWVjY6G+
fiC1ijb57NWb2GVeoWGpUAudtYHO2sSwJb5VRtwyny/Kwp1YuGHr7+cNrJ+P
9DFU9VGvvdRrD/0058ZG0maqP7MTWl2FRUX16wv6u0YcITJ6OOgY7YyjHUzQ
XbJWEXXKIGuUsSq9LhE1yWKjXosssRYRNchJWxdOs/7Qe/icrDf0ukP08i07
rDhpVdVt/qDtlyxS7SxSZ1AjTkYt+smvKNYf796AVfllAVsU6y+o9BeT+rOU
pi9RDvWVcg3BfWoPdmJO1va/n7ewQepvgPoboP76qb9eprecy0xNm0fxonbH
qwqqe7fE4M10GoG7Oxl3BNZObIsQGUYltha4erJOx9OLzQkcatMxtMDPp1xR
iZslZu70Eyd3EScLrOwiMfd0WXG2w4rTLusVRSr8YcnP8ICsXif1g9c6dbjM
Ov9pkl09q6vQLuGZ3v+WN8GIUmFSHfJuwbNMO9P5beV4EEdv3ar0x2A1RN0N
UncDTI+9Haz5X7iAjb2iblLR8H1G6tBhyvsUf05aXyczeoCQJ0zrixHmJAmu
dUg4XCugYA4zu3NYIPw7QOh3hLBPQD4B9yTU86UkzNsIxrAeYi4LB3VoF/YY
sG6NsG65I1/bFU4Z7/+j3XqRe/9ZY8td7AI8R5d2qTB6bmws1mFc6PD1Sk3X
P6nBM6U2y6cPQmsqUyqkBkYYFofI2gBDYrfTTHuFDKP64AtFuxmsKvTxFBoS
irT7AaXddgYZn8RMNnLbITDS9iwSFd3IVfViYFcfRhsHMNUygLk9g9hvH5R9
Isc8AzjZ1Y9lv+gn7cNquAfr0ZwOnuKM8gmKk2SiTBIwJulN8XYdLIWsONVl
VZZwuTxQ/c0fWvRq994V45xBbCO8OKC5DQXnDVX4emdRnJRG+gal4IeDRYcB
aRZLv7gJv5fX3ShT9zAhXL/HtHFBNTZddiET3eS2iZHLWzpbv6u49tOMiTII
vZyEXN7SXkS29yNVPYieumEMNY5gvHUUM9YxLDjGcNA9jqXOMZzwjeF0cBRn
IiNYi4n2oyFspEQrUh9TZw/DFBFrjgrPUeG5sN4Tm6FZJ+3ycPjX+Qr4nmX9
DKWwldCpAlrBekUE8OUV++a81f58FbYxDadEXbdvCdqeqoJSx+j0Q968DYv3
rmB0LBVhGt0UjUT9/6zeomk+gopHShXfMdpxhgVSD9rLBuAvH0a0ahSZunH0
NU1ipHUKk7ZpzDpnsOiZwWHvDI4FZnAyPIOV2DRWE1NYT09iIztOFxxjFh1h
NBuixfSTqOp+qro/idZb7kLXD57V6+G7T+knRoWNBZ/SW8Feu1ReiqodyS24
m6sweEBtOvSw9nrgVoU3qbMJBtPRYBHQMa8Ziu9dXmS8l5gBIj1InBMQLjk3
XXGM6seqlUdlqcg+AqYheMrHEKyaRHzXtP4qw5WM9JNcdOpoeBQttz+Oztur
Eb2PeOqlKHp3RQt7A4Oa/roGp9KV0aGAd+Qd6nbzsDqBEl2jDkteBb0N9chV
P18kmhl+ubCvbjn/duHhguSmzY1bsPMxC+qfFEvXR0WMwFY2oSuuoFtpzToD
V8lT/rtKXpSVa4aW10Nn7m8fQF+oB5k4scWocBMf1nud8D9X7lfT1D78gnI8
r4LRsaKNQsXbbayAuxeVBvJbi0KNPbtKi5g2Qymvxh7MnPQS42v5XfVdj7bA
8oSIuyNGk1LxUPl/9Wneeg4vgr6qWhVFLVt++5JsxJDng3ppz7gyCm1mkFmR
bjZFhY2HsTFEmFlWDr9aOF1b7rwJvVNZR/FUtX2KHdHE0dGG8Ge3fsjQmBc9
lvJXr7GCEs4zkoKKipyHKnpkLyyP0/kPDBjt2edZ4EVqdvHnK6rH8fEv+rHt
1kOFU8uYrO755dIS0YeevOdpJH/1GJK/fADJu36FxE9+Cv8Pvv9ZvFyy9c/w
UsnWL+DFkq1X4QU+fTFuL7mNdKf87Z3ypz6buMqT+JOcV20mQ2PyEPZaQ0sh
Ga96GneYvsG75VU43qWmjmdig+cmy6vUmHUP2WgKDIHzNIX9Q+fj7cJA+jxl
P9ivd7e0GMX/FZh8eZYjLVlj8h2iI03tmNvZ9CFxm5+9NHsLf+x7E3/MvAeT
L85w9NLv7iE55W/vlD9fL0cXV3nK/y46DBfz/TJh+wcM1QVkhOvdXWX60rIC
udKWLv9PaustahHqHnRgY5Loah9dZY5AYGH4FWdb54+hSnqCEm+/4GaytxeT
2SQGIp3izn3w9WtfgL9fPP/o/X55Jejt4qeyefG77vj5TwWJzpYPGHE+rHKx
cFS9KJF2ovzr/Iyuj3uOuxXCa36X4kwfxR1lHJgghJtmEt5HkWeHSeNmmtpM
e3kLE439vc4gevY6ELNU+/MO/ADjyyPyrRfVBGC26WDWobZF1etBletdEq+a
ndiqSLLlvKXP/znHj17VawVvVSpaztJrBgj7Rphxx+k9E1TTJNU0RfwyNVYo
EuTVGxXbLfIVh0FSN7ptdoQslciro/a+p2F5WFfHO41vFhbnQ4boYRWSPSrL
W9X4pkXy+Vo8d9BLjatLjCtdxCtxImbFRjeLtX4C1yEC+RGBk8RGUy+JiHiU
FjHCeDY8qTZfL0O4Zgqd5ROsvUdYgw+QuR5Z1GYdNhgIS1pWVP72ebUwevEW
KFpPAZGtv00us3TxQeW5R3ws1ZMUIefS9xf6Q/r7YQNi70G8K8bV6mcK7KMY
vUyH3eNEopNEpFM4m5jGcnRG2c2VGNs7jYHmKeQsk0jsHEewYhQegjVbaT/x
Sk6KKEJ50t5miLkxxDnH5Hu36g3zSxXYCymP9ajtFJs6xzRDXh9UZn+addx6
wsoyhNJkuSC5gHpbIqG/npfr5v1e+creRmqQzw5jNTaKlfAoTgXGcMw7jsMd
49jfPo5Z+zgm94xhpHkUfZYRZHcOI1Y5CD+laadl7iGmb5KbOSFV/NkMoQZY
IQ77ZR3yObXmulB5o+xUGFw3zLcV5S9pXPKYm9FryWPFmQBXJ0qZ4k59lzvF
bJ4KySY6LZlkTSsOBcQhQTfWwr04E+zHad8ATnQNYMk9iIPOQSzYBjHTNojx
3QMYqu9Hj3iHorIXYVnyZVj6Jel1UcaggGKOXqfcfb2H8w5ylUbColr9Y0Oe
ZJFMeZt0y8WSe0xKKL3h5f044aYoXooSpCghO7QIS28GeS0izrM4WSjKz+JY
8ydx1pvGcmcGJz1ZHHPpZ2b793ZjrqUbU405jNbl0F+TRa4ig0RZCkEuQwfZ
t3Ipmgyjkb4RUvHgbIaRqI/LMsjphsNyY+MzCvi8UhSvGsJpLJHqP9AbuN+N
E079FdrVTkokpPIyF/jcWPeKU2Mf7wdwxhPCsjuCk84ojtljOLI3jgOtCczv
TmDGksB4bQJDVXH0lMeQLosgwum9nN5FCVplUeKSiowot5Au2seFGKBRD0b0
7a2h1KcMKVLQNxfPl6TDMDY5nhIl/1bXCasVpx3iVQ6x/3G23cbQpr/S4cKK
04PT9k6ctHpl74bo4xA9HaK/Q/R6TNYEMVIZQP+2ALJlohPCiwDtx0PtiXep
d3POiwsyXIa1NFe+262Hmr6gvrUpw01KMXQVjTCjhImrDGYmjDO/wqoz7pN5
cdqscjf89F4hzqk9NlzNe3Ycb3HIbfPDhF0HLG4s1Hmwr8aDyR0ejJZ7MFDm
QTftN8nRwxy9S7bF2LkSNtkuw5liStNn41zwDL0i59H3tMXedi9F7KMgfSlS
Rv3HJW9geM9QCykOkSwSKGQikCMvUL7PUP0PGkcbxUnB0QarPDARxycH62xY
3Cla7eyY2WHH5HYHRrc5KIKDItiRIkXkK/U2RqrCixAcOyGHfj0Dn00FSQaU
HEXoFoGSIvTEVLBMiYB5nZJiL0vv1tK0fEW6iVI0SBuLKHDhV85/riTS1pQo
uqe9UZ6VLezQz872lVsxtc2K8TL9FY4B9RZKihQh+dRrHHtLrWnFxtkwFR+n
BaXIeZqcZ8h5hpxnxQ4TLSnLBcgwJGb4jTShRjqnUsObaJI5elaGZpmmQlIM
dwkuTIyKiRQtiE+J4nmlKKq1Nf9fo/Sq/5ZB/Nar5liP8PeYXT+ITIoDJFp5
yq/CdUSF7ERR2M7ooTvWrfqe34LQthz8ZVnaX4bip2VotpHXNvLaTONpojc0
GGr3qnDtVnigiF/VKSsT5NsUKltnftSidv2QOObR372L+2Q3nf6+PZmMkcko
FRkVB+0pbITTWAtlsBrM4mwghzP+bqz4etQB09vRV9OD7qpuZCpzSJXnEN+e
RWRbBsGyNHwMzp0UoJ0246AAVrnVGKbiQ1R8QO7xFmzHbURcKcB7igTQY9Zn
VPag7YbpfhHGkQitIOLVs0gkqGeScETPJqE41oMJrAWSWPWnmH1Eh0VGZpZT
HVlmJNGBkcNRl9ia78YhWzcOWLuxuKcb863dmGW2md7djUlmHLGVP8KsM1Sr
Zx6xxS+yj9juT1LQGC0qTEEDsvMhJrORi5Zlp7B7uFotFFZkpoZzhDVWTMr7
3iJ59aD2adnVsUEAoAUdJC5akIsWZNYMUuagn58FsOEP8TnRbxPFqjeKs10x
rHTGsdwhenOS8n3D4y7Rs5OW7x0epj8d3JvG/rYMFlr1dxD3NWVk34/4LyTG
dqUxvDONweoU+nak0FNBZLo9gVRZAgnKF6V8IdkVI7pjwrTSEK00SCv100p9
tNIueTTSYHiRvqgyEbyvSEg9OJTI46J1kUR9NpKT1E7yYN3XyfterHl9fEa0
ZgVxtiOEFU9YJtbT7VGcdIl2rhgTbBxLtjiTbAKH9iRkol1s0ZOtOMISCXey
nkl3VwKjO5l4q+MY2BFHb0UMuW0xZMqiXMgIF1I0S4mmKdE8FeRCivMdP63X
S+vtkol5NwVrKAgm5fpAkVy6530MawQIa52C7CQnqR2rHR5SB856unDG7cVK
uw/LLj9OOwM45QjihD2EYzbRBxjGkT2iLzCCgy0R7N8dxUJTFHMNUeyzRDFV
F8VEbVT2Fg5XRTBYGUF/eQQ928LMxCECihAXS7TPiTY60U4n2upEe51os+uk
YXbQMN00zHYapmgkc1wkY4eUxgDWxWLp/vhRggPiHrf+3ufZdqcECmdcbgkW
lh2dBAxdOGkTTZ6i2VO8D+qXGx6HWwI42CwaQ4NYaAxizqKDiOldegOpaCQd
3RHEUEUQA9sD6CsLoJugIkNKkn0BLiKy7VK0X4o2TNGOKdoy2+X7pG0Uo4V2
1iTTqJ5KP1zEv25612GFYGeFlday3U4S7bqidbcdJ/e6cWKPB8fbOnC0VbT4
duJwcxcONnXhQKMXixYv5uv1tuCZWh+manyYqPJhrNKL4XIvBrd5ybMXPWVd
zAxdSJPHBHmMymZa0VTr5hK0cwlcXAInl8DB5GYn3zb6jI4yfqeI32tVaSPQ
zEniGIliWh041uIkkhEt2e0SzRxsdEtEs1jvwTxRzWytaOfuwFR1B9FNh2zz
HtnegaFtHvQT5fRKpONBRqGdGClCnoLy/VXRNi3ap+3yPVabPGC3qpMVna8P
Kwx8dDchSZPojhdd8nYctOhd84u7nFiodWJupwv7qvUO+8lKveNedN6LDvyB
MpfsyO/mfFlSinOKjv2oKDxIfs7bJRsjrB8rmvr3VKJd3EkMUWPFXLUN+6ps
mK60YbLCxinEiwl2DJfZMVhmQx+HEC8u5Ehpknih4VNF4/1X+E8uX7t+7fq1
69euX7t+7fq3XON1/xcXVjEE\
\>"]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "now", " ", "add", " ", "O", " ", "lone", " ", "pairs", " ", 
   "perpendicular", " ", "to", " ", "COH", " ", "plane", " ", "and", " ", 
   "retransform", " ", "into", " ", "molecular", " ", "skeleton"}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.5716782232378025`*^9, 3.571678237995428*^9}, {
  3.571678285138711*^9, 3.5716783310963917`*^9}}]
},
WindowSize->{882, 857},
WindowMargins->{{1, Automatic}, {Automatic, 2}},
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, 2899, 65, 232, "Input"],
Cell[3459, 87, 1046, 17, 31, "Input"],
Cell[CellGroupData[{
Cell[4530, 108, 1287, 26, 31, "Input"],
Cell[5820, 136, 412, 6, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6269, 147, 2156, 53, 172, "Input"],
Cell[8428, 202, 1160, 27, 56, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9625, 234, 2580, 67, 192, "Input"],
Cell[12208, 303, 1118, 26, 86, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13363, 334, 1555, 40, 132, "Input"],
Cell[14921, 376, 720, 15, 31, "Output"]
}, Open  ]],
Cell[15656, 394, 817, 14, 31, "Input"],
Cell[CellGroupData[{
Cell[16498, 412, 1386, 29, 52, "Input"],
Cell[17887, 443, 419, 6, 31, "Output"]
}, Open  ]],
Cell[18321, 452, 836, 17, 52, "Input"],
Cell[CellGroupData[{
Cell[19182, 473, 1163, 35, 132, "Input"],
Cell[20348, 510, 1130, 26, 56, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21515, 541, 1985, 60, 192, "Input"],
Cell[23503, 603, 1162, 28, 86, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24702, 636, 1377, 38, 132, "Input"],
Cell[26082, 676, 701, 15, 33, "Output"]
}, Open  ]],
Cell[26798, 694, 721, 13, 31, "Input"],
Cell[CellGroupData[{
Cell[27544, 711, 1060, 24, 52, "Input"],
Cell[28607, 737, 412, 6, 31, "Output"]
}, Open  ]],
Cell[29034, 746, 735, 15, 52, "Input"],
Cell[CellGroupData[{
Cell[29794, 765, 1167, 35, 132, "Input"],
Cell[30964, 802, 1118, 26, 60, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32119, 833, 1685, 55, 192, "Input"],
Cell[33807, 890, 1126, 26, 86, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34970, 921, 1147, 34, 152, "Input"],
Cell[36120, 957, 707, 16, 33, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[36864, 978, 4853, 131, 392, "Input"],
Cell[41720, 1111, 14988, 255, 345, 1174, 27, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]],
Cell[56723, 1369, 397, 8, 52, "Input"]
}
]
*)

(* End of internal cache information *)