(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 48146, 1359] NotebookOptionsPosition[ 44617, 1229] NotebookOutlinePosition[ 45045, 1247] CellTagsIndexPosition[ 45002, 1244] WindowFrame->Normal ContainsDynamic->True *) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Two - Layer Geodesic Dome", "Section", CellChangeTimes->{{3.424127614465*^9, 3.4241276194*^9}}], Cell["\<\ G.W. Hart, 3 July 2008, Visualization of one triangle of a two-layer dome \ (Montreal Expo67 style, with dual triangular and hexagonal layers) and edge \ length calculation for MathCamp workshop. Start with P1, P2, P3 as three \ vertices of one face of a regular icosahedron. Point p[i,j,n] is the (i,j)th \ one of a n-frequency division of the triangle. Vertex v[i,j,n] is the \ projection of that point to a sphere of radius Rv; these are the vertices of \ the triangle layer. In the hexagon layer, on sphere of radius Rab, are points \ a[i,j,n] (centers of the the \"up\" triangles\") and points b[i,j,n] (centers \ of the \"down\" triangles. In addition, mirrored across the big triangle edge \ from point a[0,j,n] is point aa[0,j,n]. The d[,] function gives the distance \ between two points. \ \>", "Text", CellChangeTimes->{{3.4241804760786*^9, 3.4241808297496*^9}, { 3.4241808629175997`*^9, 3.4241809409156*^9}, {3.4241809715716*^9, 3.4241810710216*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"tau", "=", RowBox[{ RowBox[{"(", RowBox[{"1", "+", RowBox[{"Sqrt", "[", "5", "]"}]}], ")"}], "/", "2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"unit", "[", "vec_", "]"}], " ", ":=", " ", RowBox[{"vec", "/", RowBox[{"Sqrt", "[", RowBox[{"vec", ".", "vec"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"d", "[", RowBox[{"v1_", ",", "v2_"}], "]"}], ":=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"(", RowBox[{"v1", "-", "v2"}], ")"}], ".", RowBox[{"(", RowBox[{"v1", "-", "v2"}], ")"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.4241212341949997`*^9, 3.424121253012*^9}, { 3.424124259953*^9, 3.4241242915109997`*^9}, {3.4241811003526*^9, 3.4241811029416*^9}, 3.4241812938606*^9}], Cell[BoxData[ RowBox[{ RowBox[{"P1", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "tau"}], "}"}]}], ";", " ", RowBox[{"P2", "=", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "tau"}], ",", "1"}], "}"}]}], ";", " ", RowBox[{"P3", "=", RowBox[{"{", RowBox[{"1", ",", "0", ",", "tau"}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.4241212581940002`*^9, 3.424121260449*^9}, { 3.4241212996219997`*^9, 3.424121346251*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"p", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{ RowBox[{"P1", " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"i", "/", "n"}], ")"}], " ", RowBox[{"(", RowBox[{"P2", "-", "P1"}], ")"}]}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{"j", "/", "n"}], ")"}], " ", RowBox[{"(", RowBox[{"P3", "-", "P1"}], ")"}]}]}], "//", "N"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"v", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{"Rv", " ", RowBox[{"unit", "[", RowBox[{"p", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{"Rab", " ", RowBox[{"unit", "[", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "3"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"p", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}], ",", "n"}], "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"b", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{"Rab", " ", RowBox[{"unit", "[", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "3"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"p", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{ RowBox[{"i", "-", "1"}], ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{"i", ",", RowBox[{"j", "-", "1"}], ",", "n"}], "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"aa", "[", RowBox[{"0", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ".", RowBox[{"a", "[", RowBox[{"0", ",", "j", ",", "n"}], "]"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.4241213563570004`*^9, 3.424121419979*^9}, { 3.424121530586*^9, 3.4241215310699997`*^9}, {3.424122284855*^9, 3.424122291459*^9}, {3.424124275324*^9, 3.424124281925*^9}, { 3.4241811583956003`*^9, 3.4241811915556*^9}}], Cell[CellGroupData[{ Cell["Edge Lengths for MathCamp model:", "Subsection", CellChangeTimes->{{3.4241815524266*^9, 3.4241815778456*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"n", "=", "4"}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Rv", "=", "295"}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Rab", "=", "265"}], ";"}], " "}]}], "Input", CellChangeTimes->{{3.424121495288*^9, 3.42412152316*^9}, {3.424121609177*^9, 3.42412161055*^9}, {3.424121754525*^9, 3.424121782152*^9}, 3.4241218143929996`*^9, {3.424122057724*^9, 3.4241220734370003`*^9}, { 3.424122171683*^9, 3.4241222433710003`*^9}, {3.42412235792*^9, 3.424122374708*^9}, {3.4241224842060003`*^9, 3.4241224843310003`*^9}, { 3.4241225262279997`*^9, 3.4241225758780003`*^9}, {3.42412260794*^9, 3.4241227340290003`*^9}, {3.424122840849*^9, 3.424122843008*^9}, { 3.424122919572*^9, 3.424123011676*^9}, {3.424123192273*^9, 3.424123268769*^9}, {3.424123319487*^9, 3.424123322645*^9}, { 3.42412336474*^9, 3.42412340021*^9}, {3.424123432013*^9, 3.4241234356800003`*^9}, {3.424123466285*^9, 3.4241234811289997`*^9}, { 3.424123524457*^9, 3.424123685777*^9}, {3.4241237546940002`*^9, 3.4241238758129997`*^9}, {3.424127485718*^9, 3.4241275066280003`*^9}, { 3.424130320594*^9, 3.424130324212*^9}, {3.4241812597546*^9, 3.4241813048196*^9}, {3.4241813446526003`*^9, 3.4241813818906*^9}, { 3.4241816011506*^9, 3.4241816064955997`*^9}}], Cell[CellGroupData[{ Cell["Triangle layer :", "Subsubsection", CellChangeTimes->{{3.4241816764326*^9, 3.4241816799896*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241814597946*^9, 3.4241814774456*^9}}], Cell[BoxData["74.68945575625698`"], "Output", CellChangeTimes->{{3.4241814682076*^9, 3.4241814783786*^9}, 3.42421339655*^9, 3.4242149538274*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"0", ",", "2", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241814928255997`*^9, 3.4241815171466*^9}}], Cell[BoxData["88.08349942819889`"], "Output", CellChangeTimes->{{3.4241814955456*^9, 3.4241815180976*^9}, 3.424213396567*^9, 3.4242149538464003`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"1", ",", "0", ",", "4"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{3.4241815245986*^9}], Cell[BoxData["87.09633360907624`"], "Output", CellChangeTimes->{3.4241815265316*^9, 3.424213396585*^9, 3.4242149538654003`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241816192496*^9, 3.4241816209166*^9}}], Cell[BoxData["86.88659595300764`"], "Output", CellChangeTimes->{3.4241816220136003`*^9, 3.4242133966029997`*^9, 3.4242149538864*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "2", ",", "4"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241816419006*^9, 3.4241816431816*^9}}], Cell[BoxData["92.2963343737362`"], "Output", CellChangeTimes->{3.4241816442236*^9, 3.4242133966210003`*^9, 3.4242149539054003`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"1", ",", "2", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241816602386*^9, 3.4241816650296*^9}}], Cell[BoxData["95.85131038870736`"], "Output", CellChangeTimes->{3.4241816660716*^9, 3.424213396641*^9, 3.4242149539264*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Hexagon layer :", "Subsubsection", CellChangeTimes->{{3.4241816764326*^9, 3.4241816799896*^9}, { 3.4241817142106*^9, 3.4241817153445997`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241814597946*^9, 3.4241814774456*^9}, { 3.4241817328636*^9, 3.4241817420796003`*^9}}], Cell[BoxData["41.14689850485721`"], "Output", CellChangeTimes->{3.4241817437525997`*^9, 3.4242133966800003`*^9, 3.4242149539463997`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{3.4241817922076*^9}], Cell[BoxData["47.46971075366911`"], "Output", CellChangeTimes->{3.4241817934765997`*^9, 3.4242133967060003`*^9, 3.4242149539684*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "2", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241818092566*^9, 3.4241818112586*^9}}], Cell[BoxData["48.035703921066045`"], "Output", CellChangeTimes->{3.4241818123456*^9, 3.424213396729*^9, 3.4242149539884*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "2", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241818329396*^9, 3.4241818482036*^9}}], Cell[BoxData["49.9341204125193`"], "Output", CellChangeTimes->{3.4241818490446*^9, 3.4242133967539997`*^9, 3.4242149540094*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}], ",", RowBox[{"aa", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.424127690744*^9, 3.42412772649*^9}, 3.424127795519*^9, { 3.424127857132*^9, 3.424127950242*^9}, {3.424128083012*^9, 3.424128123384*^9}, {3.424128181267*^9, 3.424128217164*^9}, { 3.424128299683*^9, 3.424128310049*^9}, 3.424128489365*^9, { 3.4241286623859997`*^9, 3.424128664984*^9}, {3.4241287833129997`*^9, 3.424128784973*^9}, {3.4241288374300003`*^9, 3.424128844308*^9}, { 3.4241288818190002`*^9, 3.424128940454*^9}, {3.424129376367*^9, 3.424129387873*^9}, {3.424129420144*^9, 3.42412942223*^9}, { 3.424129790242*^9, 3.424129815797*^9}, {3.4241814205296*^9, 3.4241814221935997`*^9}}], Cell[BoxData["41.022708835594216`"], "Output", CellChangeTimes->{{3.424127707269*^9, 3.424127727625*^9}, 3.424127764729*^9, 3.4241277964960003`*^9, {3.4241278635290003`*^9, 3.424127950972*^9}, { 3.424128098007*^9, 3.42412812413*^9}, {3.424128182502*^9, 3.424128218324*^9}, 3.424128311358*^9, 3.4241284906689997`*^9, 3.424128665933*^9, 3.424128786593*^9, {3.4241288387320004`*^9, 3.424128845987*^9}, {3.42412890341*^9, 3.424128942139*^9}, 3.4241293890959997`*^9, 3.424129423724*^9, 3.424129623329*^9, 3.424129691549*^9, {3.42412979356*^9, 3.4241298174849997`*^9}, { 3.4241813969676*^9, 3.4241814237676*^9}, 3.42421339679*^9, 3.4242149540854*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"a", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"aa", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input"], Cell[BoxData["44.88682891689037`"], "Output", CellChangeTimes->{3.4241814275346003`*^9, 3.424213396815*^9, 3.4242149541394*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Middle layer :", "Subsubsection", CellChangeTimes->{{3.4241816764326*^9, 3.4241816799896*^9}, { 3.4241817599886*^9, 3.4241817608066*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.424127690744*^9, 3.42412772649*^9}, 3.424127795519*^9, { 3.424127857132*^9, 3.424127950242*^9}, {3.424128083012*^9, 3.424128123384*^9}, {3.424128181267*^9, 3.424128217164*^9}, { 3.424128299683*^9, 3.424128310049*^9}, 3.424128489365*^9, { 3.4241286623859997`*^9, 3.424128664984*^9}, {3.4241287833129997`*^9, 3.424128784973*^9}, {3.4241288374300003`*^9, 3.424128844308*^9}, { 3.4241288818190002`*^9, 3.424128940454*^9}, {3.424129376367*^9, 3.424129387873*^9}, {3.424129420144*^9, 3.42412942223*^9}, { 3.424129790242*^9, 3.424129815797*^9}, {3.4241814205296*^9, 3.4241814221935997`*^9}, {3.4241819249206*^9, 3.4241819267766*^9}}], Cell[BoxData["47.5553918160656`"], "Output", CellChangeTimes->{{3.424127707269*^9, 3.424127727625*^9}, 3.424127764729*^9, 3.4241277964960003`*^9, {3.4241278635290003`*^9, 3.424127950972*^9}, { 3.424128098007*^9, 3.42412812413*^9}, {3.424128182502*^9, 3.424128218324*^9}, 3.424128311358*^9, 3.4241284906689997`*^9, 3.424128665933*^9, 3.424128786593*^9, {3.4241288387320004`*^9, 3.424128845987*^9}, {3.42412890341*^9, 3.424128942139*^9}, 3.4241293890959997`*^9, 3.424129423724*^9, 3.424129623329*^9, 3.424129691549*^9, {3.42412979356*^9, 3.4241298174849997`*^9}, { 3.4241813969676*^9, 3.4241814237676*^9}, 3.4241819295406*^9, 3.4242133968389997`*^9, 3.4242149541954*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"0", ",", "0", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241819444186*^9, 3.4241819742656*^9}}], Cell[BoxData["55.189028397276815`"], "Output", CellChangeTimes->{3.4241814275346003`*^9, 3.4241819752455997`*^9, 3.424213396866*^9, 3.4242149542454*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241820020136003`*^9, 3.4241820119556*^9}}], Cell[BoxData["55.81067905004486`"], "Output", CellChangeTimes->{{3.4241820079306*^9, 3.4241820189996*^9}, 3.4242133968900003`*^9, 3.4242149542904*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241820365526*^9, 3.4241820422056*^9}}], Cell[BoxData["57.61512704968028`"], "Output", CellChangeTimes->{3.4241820428746*^9, 3.424213396914*^9, 3.4242149543274*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{3.4241820660875998`*^9}], Cell[BoxData["54.9341608418462`"], "Output", CellChangeTimes->{3.4241820668056*^9, 3.424213396941*^9, 3.4242149543704*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "2", ",", "4"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241820826275997`*^9, 3.4241820855016003`*^9}}], Cell[BoxData["58.38212328211421`"], "Output", CellChangeTimes->{3.4241820861516*^9, 3.424213396965*^9, 3.4242149544074*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"0", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241821000976*^9, 3.4241821027656*^9}}], Cell[BoxData["58.89522231645497`"], "Output", CellChangeTimes->{3.4241821034516*^9, 3.424213397*^9, 3.4242149544544*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"0", ",", "2", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "2", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241821204226*^9, 3.4241821263546*^9}}], Cell[BoxData["57.61512704968023`"], "Output", CellChangeTimes->{3.4241821273296003`*^9, 3.424213397024*^9, 3.4242149544974003`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"b", "[", RowBox[{"1", ",", "2", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241821407216*^9, 3.4241821459896*^9}}], Cell[BoxData["60.42393108579186`"], "Output", CellChangeTimes->{3.4241821469586*^9, 3.424213397049*^9, 3.4242149545284*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "[", RowBox[{ RowBox[{"v", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"1", ",", "1", ",", "4"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241821627116003`*^9, 3.4241821660416*^9}}], Cell[BoxData["60.627468477416436`"], "Output", CellChangeTimes->{3.4241821669396*^9, 3.424213397073*^9, 3.4242149545624*^9}] }, Open ]] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Dynamic Visualization", "Section", CellChangeTimes->{{3.4241275862939997`*^9, 3.424127592626*^9}, { 3.4241826488986*^9, 3.4241826566286*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"a", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{"R2", " ", RowBox[{"unit", "[", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "3"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"p", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}], ",", "n"}], "]"}]}], ")"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{"R2", " ", RowBox[{"unit", "[", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "3"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"p", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{ RowBox[{"i", "-", "1"}], ",", "j", ",", "n"}], "]"}], "+", RowBox[{"p", "[", RowBox[{"i", ",", RowBox[{"j", "-", "1"}], ",", "n"}], "]"}]}], ")"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"aa", "[", RowBox[{"0", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ".", RowBox[{"a", "[", RowBox[{"0", ",", "j", ",", "n"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"v", "[", RowBox[{"i_", ",", "j_", ",", "n_"}], "]"}], ":=", RowBox[{"R1", " ", RowBox[{"unit", "[", RowBox[{"p", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "]"}]}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"Graphics3D", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{"Blue", ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"v", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "}"}], ",", "sph"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "0", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"v", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"v", "[", RowBox[{"j", ",", "i", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"j", ",", RowBox[{"i", "+", "1"}], ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"v", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{ RowBox[{"i", "-", "1"}], ",", RowBox[{"j", "+", "1"}], ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "Red", ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"a", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "}"}], ",", "sph"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], "}"}], ",", "sph"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"a", "[", RowBox[{ RowBox[{"i", "-", "1"}], ",", RowBox[{"j", "-", "1"}], ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"a", "[", RowBox[{ RowBox[{"i", "-", "1"}], ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"a", "[", RowBox[{"i", ",", RowBox[{"j", "-", "1"}], ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "Green", ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"aa", "[", RowBox[{"0", ",", "j", ",", "n"}], "]"}], ",", "sph"}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"0", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"aa", "[", RowBox[{"0", ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", "Yellow", ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}], ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"n", "-", "i", "-", "1"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{ RowBox[{"i", "-", "1"}], ",", "j", ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b", "[", RowBox[{"i", ",", "j", ",", "n"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"i", ",", RowBox[{"j", "-", "1"}], ",", "n"}], "]"}]}], "}"}], ",", "cyl"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"n", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"n", "-", "i"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", "}"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}]}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"n", ",", "4"}], "}"}], ",", "2", ",", "20", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"R1", ",", "295"}], "}"}], ",", "1", ",", "500"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"R2", ",", "265"}], "}"}], ",", "1", ",", "500"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"sph", ",", "5"}], "}"}], ",", ".01", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"cyl", ",", "3"}], "}"}], ",", ".01", ",", "10"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.4241230781879997`*^9, 3.424123116626*^9}, 3.424123167165*^9, {3.424123915243*^9, 3.424123983026*^9}, { 3.424124035014*^9, 3.424124065796*^9}, {3.4241241041210003`*^9, 3.424124119507*^9}, {3.42412432241*^9, 3.424124361152*^9}, { 3.424124671078*^9, 3.424124675612*^9}, {3.4241256740290003`*^9, 3.4241256744449997`*^9}, {3.4241822184556*^9, 3.4241822843515997`*^9}, { 3.4241823758276*^9, 3.4241823946426*^9}, {3.4241824526576*^9, 3.4241825008915997`*^9}, {3.4241825573366003`*^9, 3.4241825890546*^9}, 3.4241826710636*^9}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`cyl$$ = 3, $CellContext`n$$ = 4, $CellContext`R1$$ = 295, $CellContext`R2$$ = 265, $CellContext`sph$$ = 5, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{{ Hold[$CellContext`n$$], 4}, 2, 20, 1}, {{ Hold[$CellContext`R1$$], 295}, 1, 500}, {{ Hold[$CellContext`R2$$], 265}, 1, 500}, {{ Hold[$CellContext`sph$$], 5}, 0.01, 10}, {{ Hold[$CellContext`cyl$$], 3}, 0.01, 10}}, Typeset`size$$ = { 630., {261., 263.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`n$188$$ = 0, $CellContext`R1$189$$ = 0, $CellContext`R2$190$$ = 0, $CellContext`sph$191$$ = 0, $CellContext`cyl$192$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`cyl$$ = 3, $CellContext`n$$ = 4, $CellContext`R1$$ = 295, $CellContext`R2$$ = 265, $CellContext`sph$$ = 5}, "ControllerVariables" :> { Hold[$CellContext`n$$, $CellContext`n$188$$, 0], Hold[$CellContext`R1$$, $CellContext`R1$189$$, 0], Hold[$CellContext`R2$$, $CellContext`R2$190$$, 0], Hold[$CellContext`sph$$, $CellContext`sph$191$$, 0], Hold[$CellContext`cyl$$, $CellContext`cyl$192$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> ($CellContext`a[ Pattern[$CellContext`i$, Blank[]], Pattern[$CellContext`j$, Blank[]], Pattern[$CellContext`n$, Blank[]]] := $CellContext`R2$$ $CellContext`unit[(1/ 3) ($CellContext`p[$CellContext`i$, $CellContext`j$, \ $CellContext`n$] + $CellContext`p[$CellContext`i$ + 1, $CellContext`j$, $CellContext`n$] + \ $CellContext`p[$CellContext`i$, $CellContext`j$ + 1, $CellContext`n$])]; $CellContext`b[ Pattern[$CellContext`i$, Blank[]], Pattern[$CellContext`j$, Blank[]], Pattern[$CellContext`n$, Blank[]]] := $CellContext`R2$$ $CellContext`unit[(1/ 3) ($CellContext`p[$CellContext`i$, $CellContext`j$, \ $CellContext`n$] + $CellContext`p[$CellContext`i$ - 1, $CellContext`j$, $CellContext`n$] + \ $CellContext`p[$CellContext`i$, $CellContext`j$ - 1, $CellContext`n$])]; $CellContext`aa[0, Pattern[$CellContext`j, Blank[]], Pattern[$CellContext`n, Blank[]]] := Dot[{{1, 0, 0}, {0, -1, 0}, {0, 0, 1}}, $CellContext`a[0, $CellContext`j, $CellContext`n]]; $CellContext`v[ Pattern[$CellContext`i$, Blank[]], Pattern[$CellContext`j$, Blank[]], Pattern[$CellContext`n$, Blank[]]] := $CellContext`R1$$ $CellContext`unit[ $CellContext`p[$CellContext`i$, $CellContext`j$, $CellContext`n$]]; Graphics3D[{Blue, Table[ Sphere[{ $CellContext`v[$CellContext`i, $CellContext`j, \ $CellContext`n$$]}, $CellContext`sph$$], {$CellContext`i, 0, $CellContext`n$$}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i}], Table[ Cylinder[{ $CellContext`v[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i + 1, $CellContext`j, $CellContext`n$$]}, $CellContext`cyl$$], \ {$CellContext`i, 0, $CellContext`n$$ - 1}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i - 1}], Table[ Cylinder[{ $CellContext`v[$CellContext`j, $CellContext`i, \ $CellContext`n$$], $CellContext`v[$CellContext`j, $CellContext`i + 1, $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 0, $CellContext`n$$ - 1}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i - 1}], Table[ Cylinder[{ $CellContext`v[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i - 1, $CellContext`j + 1, $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 1, $CellContext`n$$}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i}], Red, Table[ Sphere[{ $CellContext`a[$CellContext`i, $CellContext`j, \ $CellContext`n$$]}, $CellContext`sph$$], {$CellContext`i, 0, $CellContext`n$$ - 1}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i - 1}], Table[ Sphere[{ $CellContext`b[$CellContext`i, $CellContext`j, \ $CellContext`n$$]}, $CellContext`sph$$], {$CellContext`i, 1, $CellContext`n$$}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}], Table[ Cylinder[{ $CellContext`b[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`a[$CellContext`i - 1, $CellContext`j - 1, $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 1, $CellContext`n$$ - 1}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}], Table[ Cylinder[{ $CellContext`b[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`a[$CellContext`i - 1, $CellContext`j, $CellContext`n$$]}, $CellContext`cyl$$], \ {$CellContext`i, 1, $CellContext`n$$ - 1}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}], Table[ Cylinder[{ $CellContext`b[$CellContext`i, $CellContext`j, $CellContext`n$$], $CellContext`a[$CellContext`i, $CellContext`j - 1, $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 1, $CellContext`n$$ - 1}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}], Green, Table[ Sphere[ $CellContext`aa[ 0, $CellContext`j, $CellContext`n$$], $CellContext`sph$$], \ {$CellContext`j, 0, $CellContext`n$$ - 1}], Table[ Cylinder[{ $CellContext`a[0, $CellContext`j, $CellContext`n$$], $CellContext`aa[ 0, $CellContext`j, $CellContext`n$$]}, $CellContext`cyl$$], \ {$CellContext`j, 0, $CellContext`n$$ - 1}], Yellow, Table[ Cylinder[{ $CellContext`a[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i, $CellContext`j, \ $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 0, $CellContext`n$$ - 1}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i - 1}], Table[ Cylinder[{ $CellContext`a[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i + 1, $CellContext`j, $CellContext`n$$]}, $CellContext`cyl$$], \ {$CellContext`i, 0, $CellContext`n$$ - 1}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i - 1}], Table[ Cylinder[{ $CellContext`a[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i, $CellContext`j + 1, $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 0, $CellContext`n$$ - 1}, {$CellContext`j, 0, $CellContext`n$$ - $CellContext`i - 1}], Table[ Cylinder[{ $CellContext`b[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i, $CellContext`j, \ $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 1, $CellContext`n$$ - 1}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}], Table[ Cylinder[{ $CellContext`b[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i - 1, $CellContext`j, $CellContext`n$$]}, $CellContext`cyl$$], \ {$CellContext`i, 1, $CellContext`n$$ - 1}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}], Table[ Cylinder[{ $CellContext`b[$CellContext`i, $CellContext`j, \ $CellContext`n$$], $CellContext`v[$CellContext`i, $CellContext`j - 1, $CellContext`n$$]}, $CellContext`cyl$$], {$CellContext`i, 1, $CellContext`n$$ - 1}, {$CellContext`j, 1, $CellContext`n$$ - $CellContext`i}]}, Boxed -> False]), "Specifications" :> {{{$CellContext`n$$, 4}, 2, 20, 1}, {{$CellContext`R1$$, 295}, 1, 500}, {{$CellContext`R2$$, 265}, 1, 500}, {{$CellContext`sph$$, 5}, 0.01, 10}, {{$CellContext`cyl$$, 3}, 0.01, 10}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{675., {346.5, 349.5}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{ 3.424123118243*^9, 3.424124066888*^9, 3.424124120743*^9, {3.424124341192*^9, 3.4241243619519997`*^9}, 3.4241246768120003`*^9, 3.424124867343*^9, 3.424125678314*^9, 3.4241822890776*^9, 3.4241823958086*^9, 3.4241825015956*^9, {3.4241825655636*^9, 3.4241825911345997`*^9}, 3.424213424565*^9, {3.4242149287644*^9, 3.4242149614274*^9}}] }, Open ]] }, Open ]] }, WindowSize->{860, 1006}, WindowMargins->{{-888, Automatic}, {-49, Automatic}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, Magnification->0.75, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 102, 1, 53, "Section"], Cell[695, 26, 978, 15, 65, "Text"], Cell[1676, 43, 861, 26, 54, "Input"], Cell[2540, 71, 484, 14, 23, "Input"], Cell[3027, 87, 2734, 85, 84, "Input"], Cell[CellGroupData[{ Cell[5786, 176, 116, 1, 29, "Subsection"], Cell[5905, 179, 1361, 25, 54, "Input"], Cell[CellGroupData[{ Cell[7291, 208, 103, 1, 21, "Subsubsection"], Cell[CellGroupData[{ Cell[7419, 213, 267, 7, 23, "Input"], Cell[7689, 222, 149, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7875, 229, 271, 7, 23, "Input"], Cell[8149, 238, 154, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8340, 245, 245, 7, 23, "Input"], Cell[8588, 254, 131, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8756, 261, 267, 7, 23, "Input"], Cell[9026, 270, 136, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9199, 277, 267, 7, 23, "Input"], Cell[9469, 286, 135, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9641, 293, 267, 7, 23, "Input"], Cell[9911, 302, 124, 1, 23, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[10084, 309, 151, 2, 21, "Subsubsection"], Cell[CellGroupData[{ Cell[10260, 315, 316, 8, 23, "Input"], Cell[10579, 325, 140, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10756, 332, 245, 7, 23, "Input"], Cell[11004, 341, 136, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11177, 348, 267, 7, 23, "Input"], Cell[11447, 357, 125, 1, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11609, 363, 267, 7, 23, "Input"], Cell[11879, 372, 131, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12047, 379, 844, 16, 23, "Input"], Cell[12894, 397, 679, 10, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13610, 412, 206, 6, 23, "Input"], Cell[13819, 420, 131, 2, 23, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[13999, 428, 146, 2, 21, "Subsubsection"], Cell[CellGroupData[{ Cell[14170, 434, 885, 16, 23, "Input"], Cell[15058, 452, 703, 10, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15798, 467, 267, 7, 23, "Input"], Cell[16068, 476, 156, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16261, 483, 271, 7, 23, "Input"], Cell[16535, 492, 155, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16727, 499, 267, 7, 23, "Input"], Cell[16997, 508, 124, 1, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17158, 514, 249, 7, 23, "Input"], Cell[17410, 523, 123, 1, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17570, 529, 275, 7, 23, "Input"], Cell[17848, 538, 124, 1, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18009, 544, 267, 7, 23, "Input"], Cell[18279, 553, 121, 1, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18437, 559, 267, 7, 23, "Input"], Cell[18707, 568, 135, 2, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[18879, 575, 267, 7, 23, "Input"], Cell[19149, 584, 124, 1, 23, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19310, 590, 271, 7, 23, "Input"], Cell[19584, 599, 125, 1, 23, "Output"] }, Open ]] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[19782, 608, 150, 2, 53, "Section"], Cell[CellGroupData[{ Cell[19957, 614, 14662, 398, 444, "Input"], Cell[34622, 1014, 9967, 211, 710, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)