(***********************************************************************

                    Mathematica-Compatible Notebook

This notebook can be used on any computer system with Mathematica 4.0,
MathReader 4.0, or any compatible application. The data for the notebook 
starts with the line containing 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[     15018,        480]*)
(*NotebookOutlinePosition[     15780,        507]*)
(*  CellTagsIndexPosition[     15736,        503]*)
(*WindowFrame->Normal*)



Notebook[{
Cell[TextData[{
  "\"http://mathworld.wolfram.com/SoddyCircles.html\"    ",
  ButtonBox["Soddy",
    ButtonData:>{
      URL[ "http://mathworld.wolfram.com/SoddyCircles.html"], None},
    ButtonStyle->"Hyperlink"],
  "    gives details"
}], "Text"],

Cell[TextData[{
  "cfr Gosper, R. W. \"Soddy's Theorem on Mutually Tangent Circles, \
Generalized to ",
  StyleBox["n",
    FontSlant->"Italic"],
  " Dimensions.\" ",
  Cell[BoxData[
      FormBox[
        StyleBox[
          
          ButtonBox[\(\(\(http\)\(:\)\) // \(www . ippi . com/rwg\)/
                Sodddy . htm\),
            ButtonData:>{
              URL[ "http://www.tweedledum.com/rwg/Sodddy.htm"], None},
            ButtonStyle->"Hyperlink"],
          "Text"], TraditionalForm]],
    TextJustification->0],
  ". "
}], "Text"],

Cell[CellGroupData[{

Cell[" inits", "Subsubsection"],

Cell[TextData[{
  "touch[{x1_, y1_, r1_}, {x2_, y2_, r2_}, r3_] := Block[{ },If[x1-x2==0.,\n\
",
  StyleBox["{x1 - Sqrt[-((r1 - r2 + y1 - y2)*(r1 + r2 + 2*r3 + y1 - y2)*(r1 - \
r2 - y1 + y2)*(r1 + r2 + 2*r3 - y1 + y2)*(y1-y2)^2)]/\n   (2*(y1 - y2)^2), \
(r1^2 - r2^2 + 2*r1*r3 - 2*r2*r3 - y1^2 + y2^2)/(-2*y1 + 2*y2),r3},",
    FontColor->RGBColor[0, 0, 1]],
  "\n {((x1 - x2)*(r2^2*(x1 - x2) + 2*r2*r3*(x1 - x2) + r1^2*(-x1 + x2) + \n  \
     2*r1*r3*(-x1 + x2) + (x1 + x2)*((x1 - x2)^2 + (y1 - y2)^2)) + \n     \
Sqrt[-(x1 - x2)^2*(-(r1 - r2)^2 + (x1 - x2)^2 + (y1 - y2)^2)*\n        (-(r1 \
+ r2 + 2*r3)^2 + (x1 - x2)^2 + (y1 - y2)^2)]*(-y1 + y2)",
  StyleBox["*Sign[-x1+x2]",
    FontColor->RGBColor[1, 0, 0]],
  ")/\n    (2*(x1 - x2)*((x1 - x2)^2 + (y1 - y2)^2)), \n   (",
  StyleBox["Sign[-x1+x2]*",
    FontColor->RGBColor[1, 0, 0]],
  "Sqrt[-(x1 - x2)^2*(-(r1 - r2)^2 + (x1 - x2)^2 + (y1 - y2)^2)*\n       \
(-(r1 + r2 + 2*r3)^2 + (x1 - x2)^2 + (y1 - y2)^2)] + r2^2*(y1 - y2) + \n     \
2*r2*r3*(y1 - y2) + r1^2*(-y1 + y2) + 2*r1*r3*(-y1 + y2) + \n     ((x1 - \
x2)^2 + (y1 - y2)^2)*(y1 + y2))/(2*((x1 - x2)^2 + (y1 - y2)^2)), r3}] ]"
}], "Input",
  InitializationCell->True,
  Background->GrayLevel[0.900008]],

Cell[BoxData[
    \(<< Statistics`NormalDistribution`\)], "Input",
  InitializationCell->True],

Cell["<<Statistics`DataManipulation`", "Input",
  InitializationCell->True],

Cell["<<Graphics`Arrow`", "Input",
  InitializationCell->True],

Cell["Clear[dist1]", "Input"],

Cell[BoxData[
    \(dist1 = NormalDistribution[0, 1]\)], "Input",
  InitializationCell->True],

Cell[BoxData[
    \(dist2 = NormalDistribution[\(-1\), 2^\(-4\)]\)], "Input",
  InitializationCell->True],

Cell[BoxData[
    \(rad := \((bidi = If[Random[] > 1\ , 0, 1])\) 
          2^Random[dist1] + \((1 - bidi)\) 2^Random[dist2]\)], "Input",
  InitializationCell->True],

Cell["(* rad:= 1/(0+Random[Integer,{1,2}]*1.) *)", "Input"],

Cell["{ 2/3 Sqrt[3.]-1,1/(2/3 Sqrt[3.]-1),Log[2.,2/3 Sqrt[3]-1]}", "Input"],

Cell[BoxData[
    \(\(chk = Table[rad, {10000}];\)\)], "Input"],

Cell[BoxData[
    \(\(chk = Part[Join[tes, inside], All, 1, 3];\)\)], "Input"],

Cell[BoxData[
    \(Length[chk]\)], "Input"],

Cell["\<\
yvals=Take[BinCounts[Log[2,chk], {u=  Min[Log[2,chk]],v= Max[Log[2,chk]], \
(v-u)/100}],100];\
\>", "Input"],

Cell["leftint=Table[x, {x,u,(v)-(v-u)/100, (v-u)/100}];", "Input"],

Cell["{Length[yvals],Length[leftint]}", "Input"],

Cell["\<\
ListPlot[Transpose[{leftint,yvals}],PlotJoined->True,PlotRange->All]\
\>", "Input"],

Cell["norm[li_List]:=Sqrt[li . li]", "Input",
  InitializationCell->True],

Cell["\<\
modsort2[li_List,n_]:=Block[{temp},RotateRight[Sort[li],Length[li]-Position[\
temp=Mod[ #2-#1,n,0]& @@@ Partition[Sort[li],2,1,{1,1}],Max[temp]][[-1,1]] ]]\
\
\>", "Input",
  InitializationCell->True],

Cell[BoxData[
    \(\(\( (*\ SeedRandom[2]\ *) \)\(\ \ \ \)\(SeedRandom[1]\)\)\)], "Input"],

Cell["Clear[foam]", "Input"],

Cell["the function:", "Text"],

Cell[TextData[{
  "foam[size_,spread_]:=Timing[Catch[si=108; \
tel=2;dist1=NormalDistribution[0,spread];ra1=rad;ra2=rad;\n\
tes={z[{-ra1,0,ra1},{ra2,0,ra2}],z[{ra2,0,ra2},{-ra1,0,ra1}]};prob={};inside={\
};mesh={};\n\
Do[rr=rad;candi=touch[p=First[tes[[1]]],s=Last[tes[[1]]],rr];lo=pos=1;hi=Mod[\
2,len=Length[tes],1];\nTimeConstrained[\n\
While[(prob=Position[First/@tes,hit_List/;\n     \
(Chop[0.+norm[Take[candi-hit,2]]-Last[candi+hit]]<0  ),1])=!={},\n     \
nice=",
  StyleBox["modsort2",
    FontColor->RGBColor[1, 0, 0]],
  "[Flatten[{1,prob}],len];   \n     rol=First[nice]-1;\n     \
tes=RotateLeft[tes,rol=First[nice]-1];\n     nice=Mod[nice-rol,len,1];\n     \
lo=First[nice];\n     hi=Last[nice];\n     \
candi=touch[Part[First/@tes,lo],Part[First/@tes,hi],rr];\n]; (* end while *)\n\
,4,Throw[Indeterminate]];\nnew={ \
z[First@tes[[lo]],candi],z[candi,First@tes[[hi]]]}; \ntail=Take[tes \
,{hi,-1}];\ntoinside=Take[ tes , hi-",
  StyleBox["2",
    FontColor->RGBColor[1, 0, 0]],
  "] ; tomesh=Take[ tes , hi-",
  StyleBox["1",
    FontColor->RGBColor[1, 0, 0]],
  "] ;\ninside=Join[inside,toinside];mesh=Join[mesh,tomesh]; try=tes;\n\
tes=RotateLeft[Join[ new, tail ], 2 ] ; tel++;\n(* If[Mod[tel,500]===0,",
  StyleBox["Show[Graphics[ tes/.w\[Rule]z/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  \
}&/@{p,q}]],AspectRatio\[Rule]Automatic,PlotRange->{{-si,si},{-si,si}}]",
    FontColor->RGBColor[1, 0, 1]],
  " ];",
  StyleBox[" ",
    FontColor->RGBColor[1, 0, 1]],
  "*)",
  "\n,{ size  -2 }];]]"
}], "Input",
  InitializationCell->True]
}, Open  ]],

Cell[CellGroupData[{

Cell["foam", "Subsubsection"],

Cell["\<\
still damn errors for x>1,    'cause of modsort2  goin' haywire (?)\
\>", "Text"],

Cell["foam[2^12,2^0.5]", "Input"],

Cell["\<\
Show[{Graphics[{Hue[0.8], mesh/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[0]],Graphics[ inside/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]]
,Graphics[Hue[.6]],Graphics[ tes/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
,Graphics[{Hue[0,1,0], try/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
},AspectRatio\[Rule]Automatic,PlotRange->All];\
\>", "Input"],

Cell["\<\
Show[{Graphics[{Hue[0], inside/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]]         }&/@{p,q}]}]
,Graphics[Hue[.6]],Graphics[ tes/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]]        }&/@{p,q}]] 
,Graphics[{Hue[0,1,0], try/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]]         }&/@{p,q}]}]
},AspectRatio\[Rule]Automatic,PlotRange->All];\
\>", "Input"],

Cell["\<\
Show[{Graphics[{Hue[0.8], mesh/. z[p_,q_]\[RuleDelayed]List[{              \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[{Hue[0], inside/. z[p_,q_]\[RuleDelayed]List[{              \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[.6]],Graphics[ tes/. z[p_,q_]\[RuleDelayed]List[{           \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
,Graphics[{Hue[0,1,0], try/. z[p_,q_]\[RuleDelayed]List[{           \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
} (*  /.Line[{aa:__}]:>Arrow[aa,HeadScaling->Relative] *)   ,AspectRatio\
\[Rule]Automatic,PlotRange->All];\
\>", "Input"],

Cell["\<\
Show[{
Graphics[{Hue[0], inside/. z[p_,q_]\[RuleDelayed]List[{              \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[0,1,0]],Graphics[ tes/. z[p_,q_]\[RuleDelayed]List[{           \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
} /.Line[{aa:__}]:>Arrow[aa,HeadScaling->Relative],AspectRatio\[Rule]\
Automatic,PlotRange->All];\
\>", "Input"],

Cell[BoxData[
    \(Length[Join[tes, inside]]\)], "Input"],

Cell[BoxData[
    \({Length@tes, Length@inside, Length@mesh}\)], "Input"],

Cell["{whole,side,whole+side,meadow,achieved,ratio}", "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["calculate areas", "Subsubsection"],

Cell["surface inside fence :", "Text"],

Cell[BoxData[
    \(\(meadow := \(Last[
              Plus @@ \ \((Take[tes, All\ ]\  /. \ 
                    z[p_, q_] \[RuleDelayed] 
                      Cross[p\ {1, 1, 0}, q\ {1, 1, 0}])\)]/\(-2\)\)/
          Pi;\)\)], "Input",
  InitializationCell->True],

Cell["5 Sqrt[3.]/Pi", "Input"],

Cell["surface complete inside circles", "Text"],

Cell[BoxData[
    \(\(whole := 
        Plus @@ \((Flatten[
                Take[inside, All, \(-1\), \(-1\)]\  /. \ 
                  z \[Rule] List]^2)\);\)\)], "Input",
  InitializationCell->True],

Cell["partial surfaces of fence posts", "Text"],

Cell[BoxData[
    \(side := 
      Plus @@ Chop[\((\(Partition[Take[tes, All], 2, 
                  1, {1, 1}] /. {z[p_, q_], 
                    z[q_, r_]} \[RuleDelayed] \((\((1 - 
                            w[Take[q - p, 2], Take[r - q, 2]]/Pi)\)/2 
                       Last[q]^2)\)\) /. 
              w[a_, b_] \[RuleDelayed] \((\(-Sign[
                        Last[Cross[Append[a, 0], Append[b, 0]]]]\) 
                    ArcCos[\(a . b/norm[a]\)/
                        norm[b]])\))\), \((10^\(-8\))\)]\)], "Input",
  InitializationCell->True],

Cell[BoxData[
    \(\(achieved := 
        Chop[\((whole + side)\)/meadow, \((10^\(-9\))\)];\)\)], "Input",
  InitializationCell->True],

Cell[BoxData[
    \(perfect = \(Pi/2\)/Sqrt[3. ]\)], "Input",
  InitializationCell->True],

Cell[BoxData[
    \(ratio := achieved/perfect\)], "Input",
  InitializationCell->True]
}, Open  ]],

Cell[CellGroupData[{

Cell["taken to bits", "Subsubsection"],

Cell["the dirty underside of things:", "Text"],

Cell["detecting problems?", "Text"],

Cell["rr=rad", "Input"],

Cell["\<\
candi=touch[p=First[tes[[1]]],s=Last[tes[[1]]],rr]
lo=pos=1;hi=Mod[2,len=Length[tes],1]\
\>", "Input"],

Cell[" da capo", "Text"],

Cell["\<\
(prob=Position[First/@tes,hit_List/;
     (Chop[norm[Take[candi-hit,2]]-Last[candi+hit]]<0  ),1])\
\>", "Input"],

Cell[BoxData[
    \(tes[\([{18, 19}]\)]\)], "Input"],

Cell[BoxData[
    \(tes[\([{2, 4}]\)]\)], "Input"],

Cell[BoxData[
    \(tes[\([{2}]\)]\)], "Input"],

Cell["nice=modsort[Flatten[{1,prob}],len]", "Input"],

Cell["nice=modsort2[Flatten[{1,prob}],len]", "Input"],

Cell["{lo,hi,pos,len}", "Input"],

Cell[BoxData[
    \(rol = First[nice] - 1\)], "Input"],

Cell[BoxData[
    \(\(tes = RotateLeft[tes, rol = First[nice] - 1];\)\)], "Input"],

Cell["nice=Mod[nice-rol,len,1]", "Input"],

Cell[BoxData[{
    \(lo = First[nice]\), "\[IndentingNewLine]", 
    \(hi = Last[nice]\), "\[IndentingNewLine]", 
    \(candi = 
      touch[Part[First /@ tes, lo], Part[First /@ tes, hi], rr]\)}], "Input"],

Cell["\<\
(prob=Position[First/@tes,hit_List/;
     (Chop[norm[Take[candi-hit,2]]-Last[candi+hit]]<0  ),1])\
\>", "Input"],

Cell["and back up again again untill {}", "Text"],

Cell[BoxData[
    \(\(\(\ \)\({norm[Take[candi - \((First[tes[\([2]\)]])\), 2]], 
      Last[candi + \((First[tes[\([2]\)]])\)]\ }\)\)\)], "Input"],

Cell[BoxData[
    \(\(\(\ \ \)\({Chop[
        norm[Take[candi - \((First[tes[\([2]\)]])\), 2]] - 
          Last[candi + \((First[tes[\([2]\)]])\)]\ \ ]}\)\)\)], "Input"],

Cell["after looping, it continues :", "Text"],

Cell["new={ z[First@tes[[lo]],candi],z[candi,First@tes[[hi]]]}", "Input"],

Cell["\<\
Show[{Graphics[{Hue[0], inside/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[.6]],Graphics[ tes/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
,Graphics[{Hue[0,1,0], try/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[.4]],Graphics[ new/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
},AspectRatio\[Rule]Automatic,PlotRange->All];\
\>", "Input"],

Cell["\<\
Show[{Graphics[{Hue[0], inside/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[.6]],Graphics[ tes/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
,Graphics[{Hue[0,1,0], try/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]}]
,Graphics[Hue[.4]],Graphics[ z[candi,candi]/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  }&/@{p,q}]] 
},AspectRatio\[Rule]Automatic,PlotRange->All];\
\>", "Input"],

Cell["tail=Take[tes ,{hi,-1}]", "Input"],

Cell["toinside=Take[ tes , hi-1] ", "Input"],

Cell["inside=Join[inside,toinside];", "Input"],

Cell["tes=RotateRight[Join[ new, tail ], 2 ] ", "Input"],

Cell["verder", "Text"],

Cell["new={ z[First@tes[[lo]],candi],z[candi,First@tes[[hi]]]}", "Input"],

Cell["tail=Take[tes ,{hi,-1}]", "Input"],

Cell[TextData[{
  "toinside=Take[ tes , hi-",
  StyleBox["2",
    FontColor->RGBColor[1, 0, 0]],
  "] "
}], "Input"],

Cell["\<\
inside=Join[inside,toinside]
 try=tes\
\>", "Input"],

Cell["tes=RotateLeft[Join[ new, tail ], 2 ] ", "Input"],

Cell[TextData[{
  "If[Mod[tel,1]===0,",
  StyleBox["Show[Graphics[ tes/.w\[Rule]z/. \
z[p_,q_]\[RuleDelayed]List[{Circle[Take[#,2],Last[#]], \
Line[Take[#,2]&/@{p,q}]  \
}&/@{p,q}]],AspectRatio\[Rule]Automatic,PlotRange->{{-si,si},{-si,si}}]",
    FontColor->RGBColor[1, 0, 1]],
  " ];",
  StyleBox[" ",
    FontColor->RGBColor[1, 0, 1]]
}], "Input"]
}, Closed]]
},
FrontEndVersion->"4.0 for Microsoft Windows",
ScreenRectangle->{{0, 1024}, {0, 695}},
AutoGeneratedPackage->None,
WindowSize->{1016, 668},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
Magnification->1
]


(***********************************************************************
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[1717, 49, 248, 7, 33, "Text"],
Cell[1968, 58, 547, 18, 36, "Text"],

Cell[CellGroupData[{
Cell[2540, 80, 31, 0, 43, "Subsubsection"],
Cell[2574, 82, 1218, 22, 244, "Input",
  InitializationCell->True],
Cell[3795, 106, 94, 2, 30, "Input",
  InitializationCell->True],
Cell[3892, 110, 75, 1, 30, "Input",
  InitializationCell->True],
Cell[3970, 113, 62, 1, 30, "Input",
  InitializationCell->True],
Cell[4035, 116, 29, 0, 30, "Input"],
Cell[4067, 118, 93, 2, 30, "Input",
  InitializationCell->True],
Cell[4163, 122, 105, 2, 30, "Input",
  InitializationCell->True],
Cell[4271, 126, 165, 3, 30, "Input",
  InitializationCell->True],
Cell[4439, 131, 59, 0, 30, "Input"],
Cell[4501, 133, 75, 0, 30, "Input"],
Cell[4579, 135, 63, 1, 30, "Input"],
Cell[4645, 138, 78, 1, 30, "Input"],
Cell[4726, 141, 44, 1, 30, "Input"],
Cell[4773, 144, 118, 3, 30, "Input"],
Cell[4894, 149, 66, 0, 30, "Input"],
Cell[4963, 151, 48, 0, 30, "Input"],
Cell[5014, 153, 93, 2, 30, "Input"],
Cell[5110, 157, 73, 1, 30, "Input",
  InitializationCell->True],
Cell[5186, 160, 210, 5, 48, "Input",
  InitializationCell->True],
Cell[5399, 167, 91, 1, 30, "Input"],
Cell[5493, 170, 28, 0, 30, "Input"],
Cell[5524, 172, 29, 0, 33, "Text"],
Cell[5556, 174, 1594, 37, 426, "Input",
  InitializationCell->True]
}, Open  ]],

Cell[CellGroupData[{
Cell[7187, 216, 29, 0, 43, "Subsubsection"],
Cell[7219, 218, 91, 2, 33, "Text"],
Cell[7313, 222, 33, 0, 30, "Input"],
Cell[7349, 224, 592, 14, 102, "Input"],
Cell[7944, 240, 398, 8, 84, "Input"],
Cell[8345, 250, 579, 11, 102, "Input"],
Cell[8927, 263, 355, 8, 84, "Input"],
Cell[9285, 273, 58, 1, 30, "Input"],
Cell[9346, 276, 73, 1, 30, "Input"],
Cell[9422, 279, 62, 0, 30, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[9521, 284, 40, 0, 43, "Subsubsection"],
Cell[9564, 286, 38, 0, 33, "Text"],
Cell[9605, 288, 264, 6, 30, "Input",
  InitializationCell->True],
Cell[9872, 296, 30, 0, 30, "Input"],
Cell[9905, 298, 47, 0, 33, "Text"],
Cell[9955, 300, 200, 5, 30, "Input",
  InitializationCell->True],
Cell[10158, 307, 47, 0, 33, "Text"],
Cell[10208, 309, 560, 11, 90, "Input",
  InitializationCell->True],
Cell[10771, 322, 135, 3, 30, "Input",
  InitializationCell->True],
Cell[10909, 327, 89, 2, 30, "Input",
  InitializationCell->True],
Cell[11001, 331, 86, 2, 30, "Input",
  InitializationCell->True]
}, Open  ]],

Cell[CellGroupData[{
Cell[11124, 338, 38, 0, 43, "Subsubsection"],
Cell[11165, 340, 46, 0, 33, "Text"],
Cell[11214, 342, 35, 0, 33, "Text"],
Cell[11252, 344, 23, 0, 30, "Input"],
Cell[11278, 346, 112, 3, 48, "Input"],
Cell[11393, 351, 24, 0, 33, "Text"],
Cell[11420, 353, 122, 3, 48, "Input"],
Cell[11545, 358, 52, 1, 30, "Input"],
Cell[11600, 361, 50, 1, 30, "Input"],
Cell[11653, 364, 47, 1, 30, "Input"],
Cell[11703, 367, 52, 0, 30, "Input"],
Cell[11758, 369, 53, 0, 30, "Input"],
Cell[11814, 371, 32, 0, 30, "Input"],
Cell[11849, 373, 54, 1, 30, "Input"],
Cell[11906, 376, 82, 1, 30, "Input"],
Cell[11991, 379, 41, 0, 30, "Input"],
Cell[12035, 381, 206, 4, 70, "Input"],
Cell[12244, 387, 122, 3, 48, "Input"],
Cell[12369, 392, 49, 0, 33, "Text"],
Cell[12421, 394, 147, 2, 30, "Input"],
Cell[12571, 398, 171, 3, 30, "Input"],
Cell[12745, 403, 45, 0, 33, "Text"],
Cell[12793, 405, 73, 0, 30, "Input"],
Cell[12869, 407, 591, 14, 102, "Input"],
Cell[13463, 423, 602, 14, 102, "Input"],
Cell[14068, 439, 40, 0, 30, "Input"],
Cell[14111, 441, 44, 0, 30, "Input"],
Cell[14158, 443, 46, 0, 30, "Input"],
Cell[14207, 445, 56, 0, 30, "Input"],
Cell[14266, 447, 22, 0, 33, "Text"],
Cell[14291, 449, 73, 0, 30, "Input"],
Cell[14367, 451, 40, 0, 30, "Input"],
Cell[14410, 453, 116, 5, 30, "Input"],
Cell[14529, 460, 62, 3, 48, "Input"],
Cell[14594, 465, 55, 0, 30, "Input"],
Cell[14652, 467, 350, 10, 48, "Input"]
}, Closed]]
}
]
*)




(***********************************************************************
End of Mathematica Notebook file.
***********************************************************************)