(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 20960, 770]*) (*NotebookOutlinePosition[ 21852, 801]*) (* CellTagsIndexPosition[ 21808, 797]*) (*WindowFrame->Normal*) Notebook[{ Cell[TextData["Zonohedrification"], "Title", Evaluatable->False, CellHorizontalScrolling->False, TextAlignment->Center], Cell["\<\ George W. Hart, August, 1996 http://www.georgehart.com \ \>", "Text", Evaluatable->False], Cell[TextData["Reference"], "Section", Evaluatable->False, FontFamily->"Helv", FontSize->14, FontWeight->"Bold"], Cell[TextData[ "\"Zonohedrification\" submitted to Mathematica Journal \nfor publication, by \ George W. Hart, August 1996."], "Input", PageWidth->Infinity], Cell[TextData["Discussion"], "Section", Evaluatable->False, FontFamily->"Helv", FontSize->14, FontWeight->"Bold"], Cell[TextData[ "Runs under Mathematica Kernel 2.1, Front end 2.2\nSee above reference for \ explanation of algorithm."], "Input", PageWidth->Infinity], Cell[TextData["Implementation"], "Section", Evaluatable->False], Cell[BoxData[ \(SetDirectory["\"]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(<< "\"\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(builtInPoly[polyhedron_] := Map[\(Vertices[polyhedron]\)\[LeftDoubleBracket]#1 \[RightDoubleBracket]&, Faces[polyhedron], {2}]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(vertices[poly_] := Union[Flatten[poly, 1]]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[poly_] := Show[Graphics3D[Polygon/@poly], Boxed \[Rule] False]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(cross[{ax_, ay_, az_}, {bx_, by_, bz_}] := {ay\ bz - az\ by, az\ bx - ax\ bz, ax\ by - ay\ bx}\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(\(tolerance = 0.000001; \)\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(collinear[v1_, v2_] := And@@\(\((Abs[#1] < tolerance&)\)/@cross[v1, v2]\)\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(setStar[vlist_] := Module[{selected = {}}, Scan[Function[v, If[v \[NotEqual] {0, 0, 0} && Select[selected, collinear[v, #1]&] == {}, AppendTo[selected, v]]], vlist]; Print[Length[selected], "\< zonal directions.\>"]; gStar = selected\/2]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(pToXYZ[p_] := Plus@@MapThread[Times, {gStar, p}]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(faceNormal[p_] := With[{indices = Position[p, 0]}, cross[gStar \[LeftDoubleBracket]indices\[LeftDoubleBracket]1, 1 \[RightDoubleBracket]\[RightDoubleBracket], gStar\[LeftDoubleBracket]indices\[LeftDoubleBracket]2, 1 \[RightDoubleBracket]\[RightDoubleBracket]]]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(endPoints[p_] := {\((If[#1 == 0, \(-1\), #1]&)\)/@p, \((If[#1 == 0, 1, #1]&)\)/@p}\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(allPairs[n_] := With[{indices = Table[i, {i, n}]}, Select[Flatten[Outer[List, indices, indices], 1], #1\[LeftDoubleBracket]1\[RightDoubleBracket] < #1 \[LeftDoubleBracket]2\[RightDoubleBracket]&]]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(makeFace[{i_, j_}] := With[{normal = cross[gStar\[LeftDoubleBracket]i\[RightDoubleBracket], gStar\[LeftDoubleBracket]j\[RightDoubleBracket]]}, \((approxSign[normal . #1]&)\)/@gStar]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(approxSign[x_] := If[Abs[x] < tolerance, 0, Sign[x]]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(allFaces := With[{halfTheFaces = makeFace/@allPairs[Length[gStar]]}, halfTheFaces \[Union] \(-halfTheFaces\)]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(makeEdge[i_, p_, faceNormal_] := With[{normal = cross[gStar\[LeftDoubleBracket]i\[RightDoubleBracket], faceNormal]}, MapThread[ If[#1 \[NotEqual] 0, #1, approxSign[#2 . normal]]&, {p, gStar}]] \)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(makeEdges[p_] := With[{normal = faceNormal[p]}, With[{halfEdges = \((makeEdge[#1, p, normal]&)\)/@Flatten[Position[p, 0]]}, halfEdges \[Union] \((2\ p - #1&)\)/@halfEdges]]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(makePolygon[p_] := \((pToXYZ[#1]&)\)/@ With[{edges = makeEdges[p]}, With[{vertices = Union[Flatten[endPoints/@edges, 1]]}, Select[cycleSort[Join[edges, vertices]], Count[#1, 0] == 0&]]] \)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(allEdges := Union[Flatten[makeEdges/@allFaces, 1]]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(\(lineDrawing := \((Line[pToXYZ/@endPoints[#1]]&)\)/@allEdges; \)\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(cycleSort[list_] := Module[{result = {list\[LeftDoubleBracket]1\[RightDoubleBracket]}, remaining = Rest[list]}, While[remaining \[NotEqual] {}, choices = Select[remaining, incident[#1, result\[LeftDoubleBracket]1\[RightDoubleBracket]]&]; If[choices \[NotEqual] {}, next = choices\[LeftDoubleBracket]1\[RightDoubleBracket], Print["\", list]; Break[]]; PrependTo[result, next]; remaining = Complement[remaining, {next}]]; result]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(incident[p1_, p2_] := Plus@@Abs[p1 - p2] == 1\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(Null\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(zonohedrify[poly_] := Module[{plist}, setStar[vertices[poly]]; plist = allFaces; printSummary[plist]; makePolygon/@plist]\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(printSummary[p2list_] := With[{maxSize = Max[\((Count[#1, 0]&)\)/@p2list]}, Do[With[{howMany = Length[Select[p2list, Count[#1, 0] == i&]]}, If[howMany > 0, Print["\< \>", 2\ i, "\<-gons: \>", howMany]]], {i, maxSize, 2, \(-1\)}]; Print["\", Length[p2list]]] \)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(\(polarStar[n_] := N[Table[{Sin[\(2\ \[Pi]\ i\)\/n], Cos[\(2\ \[Pi]\ i\)\/n], 1}, {i, 1, n}]]; \)\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(faceCenters[poly_] := \((Plus@@#1\/Length[#1]&)\)/@poly\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(edgeMidpoints[poly_] := Union[Flatten[\((1\/2\ \((#1 + RotateLeft[#1])\)&)\)/@poly, 1]]\)], "Input", PageWidth->Infinity, InitializationCell->True], Cell[BoxData[ \(unit[v_] := v\/\@\(v . v\)\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(\(viewStar[star_] := Show[Graphics3D[\((Line[{{0, 0, 0}, #1}]&)\)/@star, Boxed \[Rule] False]]; \)\)], "Input", PageWidth->Infinity, InitializationCell->True, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[TextData["Use the functions"], "Section", Evaluatable->False], Cell[BoxData[ \(view[builtInPoly[Icosahedron]]\)], "Input", PageWidth->Infinity], Cell[BoxData[ \(view[zonohedrify[builtInPoly[Octahedron]]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[{edgeMidpoints[builtInPoly[Cube]]}]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(\(truncatedCuboctahedron = zonohedrify[{ unit/@Join[faceCenters[builtInPoly[Cube]], edgeMidpoints[builtInPoly[Cube]]]}]; \)\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[truncatedCuboctahedron]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[{edgeMidpoints[builtInPoly[Icosahedron]]}]]\)], "Input",\ PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[builtInPoly[Cube]]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(star = N[{{1, 0, 0}, {0, 1, 0}, {\@2\/2, \@2\/2, 0}, {\@2\/2, \(-\(\@2\/2\)\), 0}, {0, 0, 1}}]; \)\), \(view[zonohedrify[{star}]]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[Nest[zonohedrify, builtInPoly[Octahedron], 3]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(star = N[{{1, 0, 1}, {\(-1\), 0, 1}, {0, 1, 1}, {0, \(-1\), 1}, {0, 0, 2}}]; \)\), \(view[zonohedrify[{star}]]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(rhombicEnneacontahedron = zonohedrify[builtInPoly[Dodecahedron]]; \)\), \(view[rhombicEnneacontahedron]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(rhombicTriacontahedron = zonohedrify[builtInPoly[Icosahedron]]; \)\), \(view[rhombicTriacontahedron]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(RTstar = gStar; \)\), \(view[zonohedrify[{Drop[RTstar, \(-1\)]}]]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[{Drop[RTstar, \(-2\)]}]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[{Drop[RTstar, \(-3\)]}]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(starFig3a = {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}, {\(-\(1\/2\)\), 1\/2, 1\/2}}; \)\), \(view[zonohedrify[{starFig3a}]]\)}], "Input", PageWidth->Infinity], Cell[BoxData[{ \(\(starFig3b = {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}, {1, 1, 0}}; \)\), \(view[zonohedrify[{starFig3b}]]\)}], "Input", PageWidth->Infinity], Cell[BoxData[{ \(\(setStar[starFig3a]; \)\), \(Show[Graphics3D[lineDrawing, Boxed \[Rule] False]]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[{ \(\(setStar[starFig3b]; \)\), \(Show[Graphics3D[lineDrawing, Boxed \[Rule] False]]\)}], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(viewStar[starFig3a]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(viewStar[starFig3b]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[rhombicTriacontahedron]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[ zonohedrify[{ unit/@Join[vertices[builtInPoly[Icosahedron]], faceCenters[builtInPoly[Icosahedron]], edgeMidpoints[builtInPoly[Icosahedron]]]}]]\)], "Input", PageWidth->Infinity], Cell[BoxData[ \(view[zonohedrify[truncatedCuboctahedron]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[{polarStar[10]}]]\)], "Input", PageWidth->Infinity], Cell[BoxData[ \(view[ zonohedrify[{ unit/@Join[vertices[builtInPoly[Icosahedron]], edgeMidpoints[builtInPoly[Icosahedron]]]}]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}], Cell[BoxData[ \(view[zonohedrify[rhombicEnneacontahedron]]\)], "Input", PageWidth->Infinity, FontFamily->"Courier", FontSize->12, FontWeight->"Plain", FontSlant->"Plain", FontTracking->"Plain", FontVariations->{"Underline"->False, "Outline"->False, "Shadow"->False}] }, FrontEndVersion->"Microsoft Windows 3.0", ScreenRectangle->{{0, 1600}, {0, 1134}}, AutoGeneratedPackage->None, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{492, 600}, WindowMargins->{{40, Automatic}, {Automatic, 5}}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, -1}}, ShowCellLabel->True, ShowCellTags->False, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False} ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1709, 49, 125, 3, 105, "Title", Evaluatable->False], Cell[1837, 54, 131, 5, 71, "Text", Evaluatable->False], Cell[1971, 61, 120, 4, 50, "Section", Evaluatable->False], Cell[2094, 67, 158, 3, 48, "Input"], Cell[2255, 72, 121, 4, 50, "Section", Evaluatable->False], Cell[2379, 78, 152, 3, 48, "Input"], Cell[2534, 83, 65, 1, 53, "Section", Evaluatable->False], Cell[2602, 86, 133, 3, 30, "Input", InitializationCell->True], Cell[2738, 91, 302, 11, 29, "Input", InitializationCell->True], Cell[3043, 104, 419, 13, 29, "Input", InitializationCell->True], Cell[3465, 119, 316, 11, 29, "Input", InitializationCell->True], Cell[3784, 132, 344, 12, 29, "Input", InitializationCell->True], Cell[4131, 146, 377, 12, 29, "Input", InitializationCell->True], Cell[4511, 160, 110, 3, 30, "Input", InitializationCell->True], Cell[4624, 165, 164, 4, 30, "Input", InitializationCell->True], Cell[4791, 171, 590, 18, 42, "Input", InitializationCell->True], Cell[5384, 191, 322, 11, 29, "Input", InitializationCell->True], Cell[5709, 204, 411, 9, 30, "Input", InitializationCell->True], Cell[6123, 215, 175, 4, 30, "Input", InitializationCell->True], Cell[6301, 221, 325, 7, 30, "Input", InitializationCell->True], Cell[6629, 230, 309, 7, 30, "Input", InitializationCell->True], Cell[6941, 239, 136, 3, 30, "Input", InitializationCell->True], Cell[7080, 244, 399, 13, 29, "Input", InitializationCell->True], Cell[7482, 259, 529, 17, 29, "Input", InitializationCell->True], Cell[8014, 278, 298, 7, 30, "Input", InitializationCell->True], Cell[8315, 287, 509, 16, 29, "Input", InitializationCell->True], Cell[8827, 305, 324, 11, 29, "Input", InitializationCell->True], Cell[9154, 318, 152, 4, 30, "Input", InitializationCell->True], Cell[9309, 324, 860, 23, 29, "Input", InitializationCell->True], Cell[10172, 349, 319, 11, 29, "Input", InitializationCell->True], Cell[10494, 362, 88, 3, 30, "Input", InitializationCell->True], Cell[10585, 367, 411, 13, 29, "Input", InitializationCell->True], Cell[10999, 382, 619, 17, 29, "Input", InitializationCell->True], Cell[11621, 401, 391, 13, 42, "Input", InitializationCell->True], Cell[12015, 416, 139, 3, 44, "Input", InitializationCell->True], Cell[12157, 421, 181, 5, 42, "Input", InitializationCell->True], Cell[12341, 428, 300, 11, 44, "Input", InitializationCell->True], Cell[12644, 441, 390, 13, 29, "Input", InitializationCell->True], Cell[13037, 456, 68, 1, 53, "Section", Evaluatable->False], Cell[13108, 459, 86, 2, 30, "Input"], Cell[13197, 463, 288, 10, 29, "Input"], Cell[13488, 475, 299, 10, 29, "Input"], Cell[13790, 487, 401, 13, 29, "Input"], Cell[14194, 502, 274, 10, 29, "Input"], Cell[14471, 514, 308, 11, 29, "Input"], Cell[14782, 527, 282, 10, 29, "Input"], Cell[15067, 539, 404, 13, 67, "Input"], Cell[15474, 554, 297, 10, 29, "Input"], Cell[15774, 566, 377, 13, 48, "Input"], Cell[16154, 581, 363, 12, 48, "Input"], Cell[16520, 595, 360, 12, 48, "Input"], Cell[16883, 609, 320, 11, 48, "Input"], Cell[17206, 622, 287, 10, 29, "Input"], Cell[17496, 634, 287, 10, 29, "Input"], Cell[17786, 646, 190, 4, 63, "Input"], Cell[17979, 652, 161, 3, 50, "Input"], Cell[18143, 657, 333, 11, 48, "Input"], Cell[18479, 670, 333, 11, 48, "Input"], Cell[18815, 683, 265, 10, 29, "Input"], Cell[19083, 695, 265, 10, 29, "Input"], Cell[19351, 707, 287, 10, 29, "Input"], Cell[19641, 719, 251, 6, 30, "Input"], Cell[19895, 727, 287, 10, 29, "Input"], Cell[20185, 739, 90, 2, 30, "Input"], Cell[20278, 743, 387, 13, 29, "Input"], Cell[20668, 758, 288, 10, 29, "Input"] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)