(* 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[ 36957, 980] NotebookOptionsPosition[ 36002, 946] NotebookOutlinePosition[ 36579, 967] CellTagsIndexPosition[ 36536, 964] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"SiF4", " ", "computation"}], ";", " ", RowBox[{"coordinates", " ", "and", " ", "radii", " ", "from", " ", RowBox[{"SiF40zb7", ".", "exe"}]}]}], ",", " ", RowBox[{"13.02", RowBox[{".2012", "/", "29.05"}], ".2013"}]}], " ", "*)"}], "Text"], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"t", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.00000", ",", "0.00000", ",", "0.00000", ",", "0.59372"}], "}"}], ",", RowBox[{"{", RowBox[{"1.71129", ",", "0.00000", ",", "0.00000", ",", "1.11757"}], "}"}], ",", RowBox[{"{", RowBox[{"2.99074", ",", "0.00000", ",", "0.00000", ",", "0.16188"}], "}"}], ",", RowBox[{"{", RowBox[{"3.29787", ",", RowBox[{"-", "0.86867"}], ",", "0.00000", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{"3.29787", ",", "0.43434", ",", "0.75229", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{"3.29787", ",", "0.43434", ",", RowBox[{"-", "0.75229"}], ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.57043"}], ",", "1.61342", ",", "0.00000", ",", "1.11757"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.99691"}], ",", "2.81970", ",", "0.00000", ",", "0.16188"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.91828"}], ",", "2.81970", ",", "0.00000", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.68979"}], ",", "3.25404", ",", "0.75229", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.68979"}], ",", "3.25404", ",", RowBox[{"-", "0.75229"}], ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.57043"}], ",", RowBox[{"-", "0.80671"}], ",", "1.39726", ",", "1.11757"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.99691"}], ",", RowBox[{"-", "1.40985"}], ",", "2.44193", ",", "0.16188"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.91828"}], ",", RowBox[{"-", "1.40985"}], ",", "2.44193", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.68979"}], ",", RowBox[{"-", "2.27852"}], ",", "2.44193", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.68979"}], ",", RowBox[{"-", "0.97551"}], ",", "3.19422", ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.57043"}], ",", RowBox[{"-", "0.80671"}], ",", RowBox[{"-", "1.39726"}], ",", "1.11757"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.99691"}], ",", RowBox[{"-", "1.40985"}], ",", RowBox[{"-", "2.44193"}], ",", "0.16188"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.91828"}], ",", RowBox[{"-", "1.40985"}], ",", RowBox[{"-", "2.44193"}], ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.68979"}], ",", RowBox[{"-", "0.97551"}], ",", RowBox[{"-", "3.19422"}], ",", "0.75948"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.68979"}], ",", RowBox[{"-", "2.27852"}], ",", RowBox[{"-", "2.44193"}], ",", "0.75948"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"u", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.00000", ",", "0.00000", ",", "0.00000"}], "}"}], ",", RowBox[{"{", RowBox[{"2.99074", ",", "0.00000", ",", "0.00000"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.99691"}], ",", "2.81970", ",", "0.00000"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.99691"}], ",", RowBox[{"-", "1.40985"}], ",", "2.44193"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.99691"}], ",", RowBox[{"-", "1.40985"}], ",", RowBox[{"-", "2.44193"}]}], "}"}]}], "}"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.5788029926339245`*^9, 3.578803000948739*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "4"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{ 3.578802768773531*^9, {3.578803122925353*^9, 3.578803129212164*^9}, { 3.5788069613120947`*^9, 3.5788069621700964`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"xy", " ", "plane"}], ",", " ", RowBox[{"inspect", " ", "overlaps"}]}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"r", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"plot2a", "=", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "13", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "13", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "18", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "18", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.006", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], " ", "}"}]}], "}"}], " ", "]"}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.5787614756353836`*^9, 3.5787615518259172`*^9}, { 3.578761904339737*^9, 3.5787619435426054`*^9}, {3.5787619982519016`*^9, 3.5787620048351135`*^9}, {3.5787621307897344`*^9, 3.5787621663577967`*^9}, {3.5787622185086884`*^9, 3.578762243905533*^9}, { 3.578762354322527*^9, 3.5787623769269667`*^9}, {3.5788027139862347`*^9, 3.578802719649045*^9}, 3.578802912855384*^9, 3.5788031356705756`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Show", "[", RowBox[{"plot2", ",", "plot2a", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4.75"}], ",", "4.75"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "4.75"}], ",", "4.75"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.5787614756353836`*^9, 3.5787615518259172`*^9}, { 3.578761904339737*^9, 3.5787619435426054`*^9}, {3.5787619982519016`*^9, 3.5787620048351135`*^9}, {3.5787621307897344`*^9, 3.5787621663577967`*^9}, { 3.5787622185086884`*^9, 3.578762243905533*^9}, {3.578762354322527*^9, 3.5787623769269667`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0., 0.}, 0.59372], CircleBox[{1.71129, 0.}, 1.11757], CircleBox[{2.99074, 0.}, 0.16188], CircleBox[{3.29787, -0.86867}, 0.75948], CircleBox[{3.29787, 0.43434}, 0.75948], CircleBox[{3.29787, 0.43434}, 0.75948], CircleBox[{-0.57043, 1.61342}, 1.11757], CircleBox[{-0.99691, 2.8197}, 0.16188], CircleBox[{-1.91828, 2.8197}, 0.75948], CircleBox[{-0.68979, 3.25404}, 0.75948], CircleBox[{-0.68979, 3.25404}, 0.75948], CircleBox[{-0.57043, -0.80671}, 1.11757], CircleBox[{-0.99691, -1.40985}, 0.16188], CircleBox[{-1.91828, -1.40985}, 0.75948], CircleBox[{-0.68979, -2.27852}, 0.75948], CircleBox[{-0.68979, -0.97551}, 0.75948], CircleBox[{-0.57043, -0.80671}, 1.11757], CircleBox[{-0.99691, -1.40985}, 0.16188], CircleBox[{-1.91828, -1.40985}, 0.75948], CircleBox[{-0.68979, -0.97551}, 0.75948], CircleBox[{-0.68979, -2.27852}, 0.75948], {DiskBox[{2.99074, 0.}, 0.08], DiskBox[{-0.99691, 2.8197}, 0.08], DiskBox[{-0.99691, -1.40985}, 0.08], DiskBox[{-0.99691, -1.40985}, 0.08], DiskBox[{0, 0}, 0.1], {Thickness[0.006], LineBox[{{-0.99691, 2.8197}, {0, 0}}], LineBox[{{0, 0}, {2.99074, 0.}}], LineBox[{{0, 0}, {-0.99691, -1.40985}}], LineBox[{{0, 0}, {-0.99691, -1.40985}}]}}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-4.75, 4.75}, {-4.75, 4.75}}]], "Output", CellChangeTimes->{ 3.5788019637965174`*^9, 3.5788021387820244`*^9, 3.5788022815066752`*^9, 3.5788029209829984`*^9, {3.5788037190648003`*^9, 3.578803743478843*^9}, { 3.5788065111577044`*^9, 3.578806524230527*^9}, 3.5788066130726833`*^9, { 3.5788066884052153`*^9, 3.5788067299792886`*^9}, 3.5788067879489903`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Electron", " ", "domains"}], ",", " ", RowBox[{"turn", " ", "model", " ", "with", " ", "mouse"}]}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"plot2b", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"r", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot2b", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.5787614756353836`*^9, 3.5787615518259172`*^9}, { 3.578761904339737*^9, 3.5787619435426054`*^9}, {3.5787619982519016`*^9, 3.5787620048351135`*^9}, {3.5787621307897344`*^9, 3.5787621663577967`*^9}, {3.5787622185086884`*^9, 3.578762243905533*^9}, { 3.578762354322527*^9, 3.578762383744179*^9}, {3.578803685321941*^9, 3.578803705399176*^9}, 3.578803737207632*^9, {3.5788064873676624`*^9, 3.5788065187237177`*^9}, {3.5788066740375905`*^9, 3.5788067259544816`*^9}, 3.578806779696576*^9, {3.578806820334647*^9, 3.57880683329827*^9}, { 3.5788068689443326`*^9, 3.5788068871183643`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{0., 0., 0.}, 0.59372], SphereBox[{1.71129, 0., 0.}, 1.11757], SphereBox[{2.99074, 0., 0.}, 0.16188], SphereBox[{3.29787, -0.86867, 0.}, 0.75948], SphereBox[{3.29787, 0.43434, 0.75229}, 0.75948], SphereBox[{3.29787, 0.43434, -0.75229}, 0.75948], SphereBox[{-0.57043, 1.61342, 0.}, 1.11757], SphereBox[{-0.99691, 2.8197, 0.}, 0.16188], SphereBox[{-1.91828, 2.8197, 0.}, 0.75948], SphereBox[{-0.68979, 3.25404, 0.75229}, 0.75948], SphereBox[{-0.68979, 3.25404, -0.75229}, 0.75948], SphereBox[{-0.57043, -0.80671, 1.39726}, 1.11757], SphereBox[{-0.99691, -1.40985, 2.44193}, 0.16188], SphereBox[{-1.91828, -1.40985, 2.44193}, 0.75948], SphereBox[{-0.68979, -2.27852, 2.44193}, 0.75948], SphereBox[{-0.68979, -0.97551, 3.19422}, 0.75948], SphereBox[{-0.57043, -0.80671, -1.39726}, 1.11757], SphereBox[{-0.99691, -1.40985, -2.44193}, 0.16188], SphereBox[{-1.91828, -1.40985, -2.44193}, 0.75948], SphereBox[{-0.68979, -0.97551, -3.19422}, 0.75948], SphereBox[{-0.68979, -2.27852, -2.44193}, 0.75948]}, AspectRatio->Automatic, ImageSize->{569., 569.}, SphericalRegion->True, ViewPoint->{-2.5956483716307917`, -1.2764925231881108`, 1.7558975964148014`}, ViewVertical->{-0.05631351470372193, -0.11498386272248912`, 0.9935719770456726}]], "Output", CellChangeTimes->{ 3.578801963812117*^9, 3.5788021387976246`*^9, 3.578802281522275*^9, 3.578802920998598*^9, {3.578803719096*^9, 3.578803743494443*^9}, { 3.5788065111733046`*^9, 3.578806524230527*^9}, 3.5788066131038833`*^9, { 3.5788066884208155`*^9, 3.5788067299948883`*^9}, 3.57880678796459*^9}] }, Open ]], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Make", " ", "Stereopair"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"t1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"t", "=", RowBox[{"t1", ".", RowBox[{"RotationMatrix", "[", RowBox[{ RowBox[{"9", "*", RowBox[{"\[Pi]", "/", "40"}]}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"r", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3a", "=", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.008", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], " ", "}"}]}], "}"}], " ", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"t1", "=", "t"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"t", "=", RowBox[{"t1", ".", RowBox[{"RotationMatrix", "[", RowBox[{ RowBox[{"\[Pi]", "/", "20"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot4", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"r", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"plot4a", "=", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.008", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"18", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"13", ",", "2"}], "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], " ", "}"}]}], "}"}], " ", "]"}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.5787614756353836`*^9, 3.5787615518259172`*^9}, { 3.578761904339737*^9, 3.5787619435426054`*^9}, {3.5787619982519016`*^9, 3.5787620048351135`*^9}, {3.5787621307897344`*^9, 3.5787621663577967`*^9}, { 3.5787622185086884`*^9, 3.578762243905533*^9}, {3.578762354322527*^9, 3.5787623985486045`*^9}, {3.578762559338087*^9, 3.578762579820923*^9}, { 3.5787626112861786`*^9, 3.5787626633122697`*^9}, {3.578801945591285*^9, 3.5788019508328943`*^9}, {3.578802270586656*^9, 3.5788022755006647`*^9}, { 3.578802348399593*^9, 3.5788023767604427`*^9}, {3.5788066046954684`*^9, 3.5788066081586747`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Plota", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", "0", "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"Plotb", "=", "Plota"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Plota", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "plot3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Plota", "[", RowBox[{"[", "2", "]"}], "]"}], "=", "plot4"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Plotb", "[", RowBox[{"[", "1", "]"}], "]"}], "=", "plot3a"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Plotb", "[", RowBox[{"[", "2", "]"}], "]"}], "=", "plot4a"}], ";"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"Plota", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Plotb", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"ImageSize", "\[Rule]", RowBox[{"{", RowBox[{"300", ",", "300"}], "}"}]}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5787614756353836`*^9, 3.5787615518259172`*^9}, { 3.578761904339737*^9, 3.5787619435426054`*^9}, {3.5787619982519016`*^9, 3.5787620048351135`*^9}, {3.5787621307897344`*^9, 3.5787621663577967`*^9}, {3.5787622185086884`*^9, 3.578762243905533*^9}, { 3.578762354322527*^9, 3.578762394196197*^9}, {3.5787626822195024`*^9, 3.5787627730428624`*^9}, {3.5787631056150427`*^9, 3.5787631785919714`*^9}, {3.5787632443460865`*^9, 3.578763365199499*^9}, { 3.57876340019036*^9, 3.5787634045583677`*^9}, {3.5787634418580337`*^9, 3.5787634787052984`*^9}, {3.578798759611289*^9, 3.578798773464113*^9}, 3.5787988083613744`*^9, {3.5787989121639566`*^9, 3.5787989241759777`*^9}, { 3.578799193370051*^9, 3.578799222620102*^9}, {3.578799277984599*^9, 3.5787993174058685`*^9}, {3.578799668905686*^9, 3.5787997662186565`*^9}, { 3.578799848618002*^9, 3.5787999633718033`*^9}, {3.5788001185472765`*^9, 3.5788001221976824`*^9}, {3.578800191227804*^9, 3.578800204394227*^9}, { 3.578800282363164*^9, 3.578800282831165*^9}, {3.5788003601137004`*^9, 3.5788003641229076`*^9}, 3.5788004043241777`*^9, {3.5788007831396437`*^9, 3.5788007864780493`*^9}, {3.5788008270381203`*^9, 3.5788008329193306`*^9}, {3.5788009047574573`*^9, 3.5788009073158617`*^9}, {3.5788009897464066`*^9, 3.578801006064035*^9}, 3.5788010440033016`*^9, {3.578801094687791*^9, 3.5788010997733994`*^9}, { 3.57880117977034*^9, 3.578801214683202*^9}, {3.5788013492490377`*^9, 3.5788013736630807`*^9}, {3.578801466561244*^9, 3.5788015139697275`*^9}, { 3.578801580566244*^9, 3.578801600471879*^9}, {3.57880167499321*^9, 3.5788017153972807`*^9}, {3.578801762010163*^9, 3.578801799543829*^9}, { 3.578802061780289*^9, 3.5788021206079926`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ GraphicsBox[{CircleBox[{0., 0.}, 0.59372], CircleBox[{1.301275124871677, 0.}, 1.11757], CircleBox[{2.2741765375586365`, 0.}, 0.16188], CircleBox[{2.507720021773374, -0.86867}, 0.75948], CircleBox[{2.0191467494950603`, 0.43434}, 0.75948], CircleBox[{2.996293294051688, 0.43434}, 0.75948], CircleBox[{-0.43375837495722563`, 1.61342}, 1.11757], CircleBox[{-0.7580563111663268, 2.8197}, 0.16188], CircleBox[{-1.4586715556912273`, 2.8197}, 0.75948], CircleBox[{-1.0130937032895593`, 3.25404}, 0.75948], CircleBox[{-0.03594715873293147, 3.25404}, 0.75948], CircleBox[{-1.341206154967058, -0.80671}, 1.11757], CircleBox[{-2.3439629838252523`, -1.40985}, 0.16188], CircleBox[{-3.0445782283501526`, -1.40985}, 0.75948], CircleBox[{-2.1104271036701707`, -2.27852}, 0.75948], CircleBox[{-2.5990003759484845`, -0.97551}, 0.75948], CircleBox[{0.47368940505260676`, -0.80671}, 1.11757], CircleBox[{0.8278503614925987, -1.40985}, 0.16188], CircleBox[{0.12723511696769818`, -1.40985}, 0.75948], CircleBox[{1.549959513925994, -0.97551}, 0.75948], CircleBox[{1.0613862416476803`, -2.27852}, 0.75948], { DiskBox[{2.2741765375586365`, 0.}, 0.08], DiskBox[{-0.7580563111663268, 2.8197}, 0.08], DiskBox[{-2.3439629838252523`, -1.40985}, 0.08], DiskBox[{0.8278503614925987, -1.40985}, 0.08], DiskBox[{0, 0}, 0.1], {Thickness[0.008], LineBox[{{-0.7580563111663268, 2.8197}, {0, 0}}], LineBox[{{0, 0}, {2.2741765375586365`, 0.}}], LineBox[{{0, 0}, {0.8278503614925987, -1.40985}}], LineBox[{{0, 0}, {-2.3439629838252523`, -1.40985}}]}}}, AspectRatio->Automatic, ImageSize->{300, 300}], ",", GraphicsBox[{CircleBox[{0., 0.}, 0.59372], CircleBox[{1.1113939506269601`, 0.}, 1.11757], CircleBox[{1.9423302560630136`, 0.}, 0.16188], CircleBox[{2.141795235146663, -0.86867}, 0.75948], CircleBox[{1.5697494312854157`, 0.43434}, 0.75948], CircleBox[{2.71384103900791, 0.43434}, 0.75948], CircleBox[{-0.37046465020898667`, 1.61342}, 1.11757], CircleBox[{-0.6474412538608433, 2.8197}, 0.16188], CircleBox[{-1.2458232021508246`, 2.8197}, 0.75948], CircleBox[{-1.020028573118925, 3.25404}, 0.75948], CircleBox[{0.12406303460356989`, 3.25404}, 0.75948], CircleBox[{-1.432949489703286, -0.80671}, 1.11757], CircleBox[{-2.504299393438527, -1.40985}, 0.16188], CircleBox[{-3.102681341728508, -1.40985}, 0.75948], CircleBox[{-2.3048409088353607`, -2.27852}, 0.75948], CircleBox[{-2.8768867126966082`, -0.97551}, 0.75948], CircleBox[{0.6920201892853125, -0.80671}, 1.11757], CircleBox[{1.2094168857168404`, -1.40985}, 0.16188], CircleBox[{0.611034937426859, -1.40985}, 0.75948], CircleBox[{1.9809211741812538`, -0.97551}, 0.75948], CircleBox[{1.4088753703200065`, -2.27852}, 0.75948], { DiskBox[{1.9423302560630136`, 0.}, 0.08], DiskBox[{-0.6474412538608433, 2.8197}, 0.08], DiskBox[{-2.504299393438527, -1.40985}, 0.08], DiskBox[{1.2094168857168404`, -1.40985}, 0.08], DiskBox[{0, 0}, 0.1], {Thickness[0.008], LineBox[{{-0.6474412538608433, 2.8197}, {0, 0}}], LineBox[{{0, 0}, {1.9423302560630136`, 0.}}], LineBox[{{0, 0}, {1.2094168857168404`, -1.40985}}], LineBox[{{0, 0}, {-2.504299393438527, -1.40985}}]}}}, AspectRatio->Automatic, ImageSize->{300, 300}]}], "}"}]], "Output", CellChangeTimes->{ 3.578801963858917*^9, 3.5788021388444247`*^9, 3.5788022815534754`*^9, 3.5788029210297985`*^9, {3.5788037191584005`*^9, 3.5788037435412436`*^9}, { 3.578806511220104*^9, 3.5788065242773275`*^9}, 3.5788066131506834`*^9, { 3.5788066885768156`*^9, 3.578806730041689*^9}, 3.5788067880269904`*^9}] }, Open ]] }, WindowSize->{898, 824}, WindowMargins->{{433, Automatic}, {Automatic, 1}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.2, 841.8}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, 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, 4613, 128, 312, "Input"], Cell[5173, 150, 1259, 42, 92, "Input"], Cell[6435, 194, 4840, 135, 172, "Input"], Cell[CellGroupData[{ Cell[11300, 333, 1056, 24, 72, "Input"], Cell[12359, 359, 1780, 36, 375, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14176, 400, 1874, 45, 92, "Input"], Cell[16053, 447, 1683, 32, 584, "Output"] }, Open ]], Cell[17751, 482, 11197, 317, 452, "Input"], Cell[CellGroupData[{ Cell[28973, 803, 3131, 69, 152, "Input"], Cell[32107, 874, 3879, 69, 317, "Output"] }, Open ]] } ] *) (* End of internal cache information *)