(************** Content-type: application/mathematica **************
                     CreatedBy='Mathematica 5.0'

                    Mathematica-Compatible Notebook

This notebook can be used with any Mathematica-compatible
application, such as Mathematica, MathReader or Publicon. 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[     39701,       1072]*)
(*NotebookOutlinePosition[     40614,       1103]*)
(*  CellTagsIndexPosition[     40525,       1097]*)
(*WindowFrame->Normal*)



Notebook[{
Cell[BoxData[{
    \(\(\(SetAttributes[Swap23, HoldAll]\)\(\n\)
    \)\), "\[IndentingNewLine]", 
    \(Swap23[L_] := 
      Module[{temp}, 
        temp = L[\([2]\)]; \[IndentingNewLine]L[\([2]\)] = 
          L[\([3]\)]; \[IndentingNewLine]L[\([3]\)] = 
          temp; \[IndentingNewLine]L]\)}], "Input"],

Cell[BoxData[{
    \(SetAttributes[InitializeMatrix, HoldFirst]\), "\[IndentingNewLine]", 
    \(\(InitializeMatrix[M_, n_] := \(M = 
          Table[If[i \[Equal] j, Infinity, \(-Infinity\)], {i, 1, n}, {j, 1, 
              n}]\);\)\), "\[IndentingNewLine]", 
    \(\(SetAttributes[SetToLambda, HoldFirst];\)\), "\[IndentingNewLine]", 
    \(\(SetToLambda[M_, i_, j_] := \ 
        If[i > j, SetToLambda[M, j, i], 
          M[\([i, j]\)]\  = \ \(M[\([j, 
                  i]\)]\  = \ \[Lambda]\_\(i\[InvisibleComma]j\)\)];\)\), "\
\[IndentingNewLine]", 
    \(\(SetAttributes[RecoverFromRowSum, 
        HoldFirst];\)\), "\[IndentingNewLine]", 
    \(\(\(RecoverFromRowSum[M_, i_, j_]\  := \ \((M[\([i, j]\)]\  = \ 0; \ 
          M[\([i, j]\)]\  = \ \(M[\([j, i]\)]\  = \ \(-Total[
                  M[\([i]\)]]\)\))\);\)\(\[IndentingNewLine]\)
    \)\), "\[IndentingNewLine]", 
    \(\(SetAttributes[SolveDet, HoldFirst];\)\), "\[IndentingNewLine]", 
    \(\(\(SolveDet[M_, i_, j_, rows_, cols_, \ 
        sieveFunction_: \((True &)\)]\  := \ 
      Module[{submat, x, soln, 
          eqn}, \ \[IndentingNewLine]M[\([i, j]\)] = \(M[\([j, i]\)] = 
            x\); \  (*\ 
          Double\ assignment\ important\ for\ the\ square - 
            root\ trick\ *) \[IndentingNewLine]submat = M[\([rows]\)]; \ 
        submat = 
          Transpose[\ \((Transpose[
                  submat])\)[\([cols]\)]]; \[IndentingNewLine] (*\ 
          Maxe\ soln\ the\ list\ of\ solutions*) \[IndentingNewLine]eqn = \
\((Det[submat] \[Equal] 0)\); \[IndentingNewLine]soln\  = \ 
          x /. Solve[eqn, x]; \[IndentingNewLine]soln = 
          Select[soln, \ sieveFunction]; \[IndentingNewLine]If[
          Length[soln] \[NotEqual] 1, \ Print[eqn]; 
          Message[recover::TooManySolutions, \ "\<Multiple or no solutions \
found where there should only be one\>"]]; \[IndentingNewLine] (*Curious\ \
\(thing : \ without\ this\ Simplify[]\), \ 
          the\ program\ is\ vastly\ \((100\ times)\)\ less\ efficient\ \
*) \[IndentingNewLine]M[\([j, i]\)]\  = \ \(M[\([i, j]\)] = 
            Simplify[
              First[soln]]\);\[IndentingNewLine]\[IndentingNewLine]]\)\(\
\[IndentingNewLine]\)
    \)\), "\[IndentingNewLine]", 
    \(\(SetAttributes[SqrtTrick, HoldFirst];\)\), "\[IndentingNewLine]", 
    \(\(SqrtTrick[M_, i_, 
          j_]\  :=  (*\(If[i > j, SqrtTrick[M, j, i], 
              SolveDet[M, i, j, {i, j}, {i, j}, 
                Function[x, \ x < 0]]];\)*) \[IndentingNewLine]\(M[\([i, 
              j]\)] = \(M[\([j, i]\)]\  = \ \(-Sqrt[
                M[\([i, i]\)]*
                  M[\([j, j]\)]]\)\)\);\)\), "\[IndentingNewLine]", 
    \(\(SetAttributes[ParametrizeIndex, 
        HoldFirst];\)\), "\[IndentingNewLine]", 
    \(\(ParametrizeIndex[M_, i_, 
          j_] := \(M[\([i, j]\)]\  = \ \(M[\([j, i]\)]\  = \ 
            Unique[]\)\);\)\), "\[IndentingNewLine]", 
    \(\)}], "Input"],

Cell[BoxData[{
    \(\(\(<< LinearAlgebra`MatrixManipulation`\)\(\[IndentingNewLine]\)
    \) (*\ \(Find\ the\ response\ matrix\ from\ Kirkoff' 
            s\ matrix\  &\)\ number\ of\ boundary\ nodes\ *) \), "\
\[IndentingNewLine]", 
    \(KtoL[K_, n_]\  := \ 
      Block[{A, B, C}, \ 
        If[Length[K] \[Equal] n, 
          K, \[IndentingNewLine]A = 
            TakeMatrix[K, {1, 1}, {n, n}]; \[IndentingNewLine]B = 
            TakeMatrix[
              K, {1, n + 1}, {n, Length[K]}]; \[IndentingNewLine]C\  = \ 
            TakeMatrix[
              K, \ {n + 1, n + 1}, \ {Length[K], 
                Length[K]}]; \[IndentingNewLine]A - 
            B . Inverse[C] . 
              Transpose[B]\[IndentingNewLine]]]\), "\[IndentingNewLine]", 
    \(LtoConditions[L_] := 
      Flatten[Table[\[Lambda]\_\(i, j\) = 
            0; \[Lambda]\_\(i, j\) =. ; \[Lambda]\_\(i, j\) -> 
            L[\([i, j]\)], {j, 1, Length[L]}, {i, 1, 
            j}]]\), "\[IndentingNewLine]", 
    \(\)}], "Input"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(K = 
        G /. {\(-Infinity\) \[Rule] 1, \ 
            Infinity \[Rule] 0};\)\), "\[IndentingNewLine]", 
    \(\(K[\([6, 9]\)]\  = \ \(K[\([9, 6]\)]\  = \ 
          2\);\)\), "\[IndentingNewLine]", 
    \(\(Map[Function[x, K[\([x, x]\)]\  = \ \(-Total[K[\([x]\)]]\)], 
        Range[Length[K]]];\)\), "\[IndentingNewLine]", 
    \(K\  // MatrixForm\)}], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {\(-2\), "0", "0", "0", "0", "0", "1", "0", "1"},
            {"0", \(-2\), "0", "0", "0", "0", "1", "1", "0"},
            {"0", "0", \(-2\), "0", "0", "0", "0", "1", "1"},
            {"0", "0", "0", \(-2\), "0", "0", "1", "0", "1"},
            {"0", "0", "0", "0", \(-2\), "0", "1", "1", "0"},
            {"0", "0", "0", "0", "0", \(-3\), "0", "1", "2"},
            {"1", "1", "0", "1", "1", "0", \(-4\), "0", "0"},
            {"0", "1", "1", "0", "1", "1", "0", \(-4\), "0"},
            {"1", "0", "1", "1", "0", "2", "0", "0", \(-5\)}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(KtoL[K, 6]\)\(\ \ \)\)\)], "Input"],

Cell[BoxData[
    \({{\(-\(31\/20\)\), 1\/4, 1\/5, 9\/20, 1\/4, 
        2\/5}, {1\/4, \(-\(3\/2\)\), 1\/4, 1\/4, 1\/2, 1\/4}, {1\/5, 
        1\/4, \(-\(31\/20\)\), 1\/5, 1\/4, 13\/20}, {9\/20, 1\/4, 
        1\/5, \(-\(31\/20\)\), 1\/4, 2\/5}, {1\/4, 1\/2, 1\/4, 
        1\/4, \(-\(3\/2\)\), 1\/4}, {2\/5, 1\/4, 13\/20, 2\/5, 
        1\/4, \(-\(39\/20\)\)}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(LtoConditions[%]\)], "Input"],

Cell[BoxData[
    \({\[Lambda]\_\(1, 1\) \[Rule] \(-\(31\/20\)\), \[Lambda]\_\(1, 2\) \
\[Rule] 1\/4, \[Lambda]\_\(2, 2\) \[Rule] \(-\(3\/2\)\), \[Lambda]\_\(1, 3\) \
\[Rule] 1\/5, \[Lambda]\_\(2, 3\) \[Rule] 
        1\/4, \[Lambda]\_\(3, 3\) \[Rule] \(-\(31\/20\)\), \[Lambda]\_\(1, \
4\) \[Rule] 9\/20, \[Lambda]\_\(2, 4\) \[Rule] 
        1\/4, \[Lambda]\_\(3, 4\) \[Rule] 
        1\/5, \[Lambda]\_\(4, 4\) \[Rule] \(-\(31\/20\)\), \[Lambda]\_\(1, \
5\) \[Rule] 1\/4, \[Lambda]\_\(2, 5\) \[Rule] 
        1\/2, \[Lambda]\_\(3, 5\) \[Rule] 1\/4, \[Lambda]\_\(4, 5\) \[Rule] 
        1\/4, \[Lambda]\_\(5, 5\) \[Rule] \(-\(3\/2\)\), \[Lambda]\_\(1, 6\) \
\[Rule] 2\/5, \[Lambda]\_\(2, 6\) \[Rule] 1\/4, \[Lambda]\_\(3, 6\) \[Rule] 
        13\/20, \[Lambda]\_\(4, 6\) \[Rule] 2\/5, \[Lambda]\_\(5, 6\) \[Rule] 
        1\/4, \[Lambda]\_\(6, 6\) \[Rule] \(-\(39\/20\)\)}\)], "Output"]
}, Open  ]],

Cell[BoxData[
    \(Unique[]\)], "Input"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(M = {{\(-2\), 1, 1}, {1, \(-2\), 1}, \ {1, 
            1, \(-2\)}};\)\), "\[IndentingNewLine]", 
    \(M // MatrixForm\), "\[IndentingNewLine]", 
    \(LtoConditions[KtoL[M, \ 3]]\)}], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {\(-2\), "1", "1"},
            {"1", \(-2\), "1"},
            {"1", "1", \(-2\)}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"],

Cell[BoxData[
    \({\[Lambda]\_\(1, 1\) \[Rule] \(-2\), \[Lambda]\_\(1, 2\) \[Rule] 
        1, \[Lambda]\_\(2, 2\) \[Rule] \(-2\), \[Lambda]\_\(1, 3\) \[Rule] 
        1, \[Lambda]\_\(2, 3\) \[Rule] 
        1, \[Lambda]\_\(3, 3\) \[Rule] \(-2\)}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(With[Evaluate[%], \[Lambda]\_\(1, 1\)]\)], "Input"],

Cell[BoxData[
    RowBox[{\(With::"lvw"\), \(\(:\)\(\ \)\), "\<\"Local variable \
specification \\!\\({\\(\\(Hold[\\(\\(\[Lambda]\\_\\(i, j\\) = \
\\(\\(\\(\\({\\(\\({\\(\\(-2\\)\\), 1, 1}\\)\\), \\(\\({1, \\(\\(-2\\)\\), 1}\
\\)\\), \\(\\({1, 1, \\(\\(-2\\)\\)}\\)\\)}\\)\\) \[LeftDoubleBracket] \
\\(\\(i, j\\)\\) \[RightDoubleBracket]\\)\\)\\)\\)]\\)\\), \\(\\(\
\[LeftSkeleton] 4 \[RightSkeleton]\\)\\), \\(\\(Hold[\\(\\(\[Lambda]\\_\\(i, \
j\\) = \\(\\(\\(\\({\\(\\({\\(\\(-2\\)\\), 1, 1}\\)\\), \\(\\({1, \\(\\(-2\\)\
\\), 1}\\)\\), \\(\\({1, 1, \\(\\(-2\\)\\)}\\)\\)}\\)\\) \[LeftDoubleBracket] \
\\(\\(i, j\\)\\) \[RightDoubleBracket]\\)\\)\\)\\)]\\)\\)}\\) contains \
\\!\\(Hold[\\(\\(\[Lambda]\\_\\(i, j\\) = \\(\\(\\(\\({\\(\\({\\(\\(-2\\)\\), \
1, 1}\\)\\), \\(\\({1, \\(\\(-2\\)\\), 1}\\)\\), \\(\\({1, 1, \\(\\(-2\\)\\)}\
\\)\\)}\\)\\) \[LeftDoubleBracket] \\(\\(i, j\\)\\) \
\[RightDoubleBracket]\\)\\)\\)\\)]\\) which is not an assignment to a symbol. \
\\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\
\\\", ButtonFrame->None, ButtonData:>\\\"With::lvw\\\"]\\)\"\>"}]], "Message"],

Cell[BoxData[
    \(With[{Hold[\[Lambda]\_\(i, j\) = {{\(-2\), 1, 1}, {1, \(-2\), 1}, {1, 
                  1, \(-2\)}}\[LeftDoubleBracket]i, j\[RightDoubleBracket]], 
        Hold[\[Lambda]\_\(i, j\) = {{\(-2\), 1, 1}, {1, \(-2\), 1}, {1, 
                  1, \(-2\)}}\[LeftDoubleBracket]i, j\[RightDoubleBracket]], 
        Hold[\[Lambda]\_\(i, j\) = {{\(-2\), 1, 1}, {1, \(-2\), 1}, {1, 
                  1, \(-2\)}}\[LeftDoubleBracket]i, j\[RightDoubleBracket]], 
        Hold[\[Lambda]\_\(i, j\) = {{\(-2\), 1, 1}, {1, \(-2\), 1}, {1, 
                  1, \(-2\)}}\[LeftDoubleBracket]i, j\[RightDoubleBracket]], 
        Hold[\[Lambda]\_\(i, j\) = {{\(-2\), 1, 1}, {1, \(-2\), 1}, {1, 
                  1, \(-2\)}}\[LeftDoubleBracket]i, j\[RightDoubleBracket]], 
        Hold[\[Lambda]\_\(i, j\) = {{\(-2\), 1, 1}, {1, \(-2\), 1}, {1, 
                  1, \(-2\)}}\[LeftDoubleBracket]i, 
              j\[RightDoubleBracket]]}, \[Lambda]\_\(1, 1\)]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(%\)], "Input"],

Cell[BoxData[
    \({Hold[\[Lambda]\_\(i, j\)] \[Rule] \(-2\), 
      Hold[\[Lambda]\_\(i, j\)] \[Rule] 1, 
      Hold[\[Lambda]\_\(i, j\)] \[Rule] \(-2\), 
      Hold[\[Lambda]\_\(i, j\)] \[Rule] 1, 
      Hold[\[Lambda]\_\(i, j\)] \[Rule] 1, 
      Hold[\[Lambda]\_\(i, j\)] \[Rule] \(-2\)}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\[Lambda]\_\(1, 1\)\)], "Input"],

Cell[BoxData[
    \(\(-2\)\)], "Output"]
}, Open  ]],

Cell[BoxData[
    \(\(With[{n = 1, K = M}, 
        A = TakeMatrix[K, {1, 1}, {n, n}]; \[IndentingNewLine]B = 
          TakeMatrix[
            K, {1, n + 1}, {n, Length[K]}]; \[IndentingNewLine]CC\  = \ 
          TakeMatrix[
            K, \ {n + 1, Length[K]}, \ {n + 1, Length[K]}]];\)\)], "Input"],

Cell[CellGroupData[{

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

Cell[BoxData[
    \(3\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(TakeMatrix[M, {2, 3}, {2, 3}]\)], "Input"],

Cell[BoxData[
    \({{1}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(MatrixForm\  /@ \ {A, B, CC}\)], "Input"],

Cell[BoxData[
    RowBox[{"{", 
      RowBox[{
        TagBox[
          RowBox[{"(", "\[NoBreak]", GridBox[{
                {\(-2\)}
                }], "\[NoBreak]", ")"}],
          Function[ BoxForm`e$, 
            MatrixForm[ BoxForm`e$]]], ",", 
        TagBox[
          RowBox[{"(", "\[NoBreak]", GridBox[{
                {"1", "1"}
                }], "\[NoBreak]", ")"}],
          Function[ BoxForm`e$, 
            MatrixForm[ BoxForm`e$]]], ",", 
        TagBox[
          RowBox[{"(", "\[NoBreak]", GridBox[{
                {"1"}
                }], "\[NoBreak]", ")"}],
          Function[ BoxForm`e$, 
            MatrixForm[ BoxForm`e$]]]}], "}"}]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(TakeMatrix[M, {1, n + 1}, {n, Length}]\  /. \ n \[Rule] 2\)], "Input"],

Cell[BoxData[
    \({{\(-2\), 1}, {1, \(-2\)}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(?? Solve\)], "Input"],

Cell[BoxData[
    RowBox[{"\<\"Solve[eqns, vars] attempts to solve an equation or set of \
equations for the variables vars. Solve[eqns, vars, elims] attempts to solve \
the equations for vars, eliminating the variables elims.\"\>", " ", 
      ButtonBox[
        StyleBox["More\[Ellipsis]",
          "SR"],
        ButtonData:>"Solve",
        Active->True,
        ButtonStyle->"RefGuideLink"]}]], "Print",
  CellTags->"Info3360852581-7681982"],

Cell[BoxData[
    InterpretationBox[GridBox[{
          {\(Attributes[Solve] = {Protected}\)},
          {" "},
          {GridBox[{
                {\(Options[Solve] = {InverseFunctions \[Rule] Automatic, 
                      MakeRules \[Rule] False, Method \[Rule] 3, 
                      Mode \[Rule] Generic, Sort \[Rule] True, 
                      VerifySolutions \[Rule] Automatic, 
                      WorkingPrecision \[Rule] \[Infinity]}\)}
                },
              GridBaseline->{Baseline, {1, 1}},
              ColumnWidths->0.999,
              ColumnAlignments->{Left}]}
          },
        GridBaseline->{Baseline, {1, 1}},
        ColumnAlignments->{Left}],
      Definition[ "Solve"],
      Editable->False]], "Print",
  CellTags->"Info3360852581-7681982"]
}, Open  ]],

Cell[BoxData[
    \(Assuming[x > 0, 
      Solve[{x^2 \[Equal] 1\ , x/Abs[x] == \(-1\)}, x]]\)], "Input"],

Cell[BoxData[
    \(Reduce\)], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
    \(x /. Solve[x^2 \[Equal] 1, x]\)], "Input"],

Cell[BoxData[
    \({\(-1\), 1}\)], "Output"]
}, Open  ]],

Cell[BoxData[
    \(SolveDet[A, 2, 3, {2, 3, 5}, {3, 4, 5}]\)], "Input"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(B = \(A = Table[Random[], {5}, {5}]\);\)\), "\[IndentingNewLine]", 
    \(A // MatrixForm\), "\[IndentingNewLine]", 
    \(\(SolveDet[A, 2, 3, {2, 3, 5}, {3, 4, 5}];\)\), "\[IndentingNewLine]", 
    \(A\  // MatrixForm\), "\[IndentingNewLine]", 
    \(A - B\)}], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {"0.4419362131396111`", "0.6857070562593893`", 
              "0.9083978020986319`", "0.40708396636658906`", 
              "0.7125074971033162`"},
            {"0.0742429021495689`", "0.9786521703832708`", 
              "0.9759604817134262`", "0.7808299449928588`", 
              "0.07440422219042807`"},
            {"0.5513787731780128`", "0.04364855733595576`", 
              "0.6682435292163451`", "0.6430648721115518`", 
              "0.32416321253467806`"},
            {"0.7756395038133836`", "0.15165370740666995`", 
              "0.2113998180126693`", "0.3267093949754902`", 
              "0.7594090101331057`"},
            {"0.9790966281033705`", "0.3598191016521737`", 
              "0.9852153923546223`", "0.12192238327085`", 
              "0.5371604149637594`"}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {"0.4419362131396111`", "0.6857070562593893`", 
              "0.9083978020986319`", "0.40708396636658906`", 
              "0.7125074971033162`"},
            {"0.0742429021495689`", "0.9786521703832708`", 
              "0.23531798519968328`", "0.7808299449928588`", 
              "0.07440422219042807`"},
            {"0.5513787731780128`", "0.23531798519968328`", 
              "0.6682435292163451`", "0.6430648721115518`", 
              "0.32416321253467806`"},
            {"0.7756395038133836`", "0.15165370740666995`", 
              "0.2113998180126693`", "0.3267093949754902`", 
              "0.7594090101331057`"},
            {"0.9790966281033705`", "0.3598191016521737`", 
              "0.9852153923546223`", "0.12192238327085`", 
              "0.5371604149637594`"}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"],

Cell[BoxData[
    \({{0.`, 0.`, 0.`, 0.`, 0.`}, {0.`, 0.`, \(-0.7406424965137429`\), 0.`, 
        0.`}, {0.`, 0.1916694278637275`, 0.`, 0.`, 0.`}, {0.`, 0.`, 0.`, 0.`, 
        0.`}, {0.`, 0.`, 0.`, 0.`, 0.`}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(InitializeMatrix[A, 5]\)], "Input"],

Cell[BoxData[
    \({{\[Infinity], \[Infinity], \[Infinity], \[Infinity], \[Infinity]}, {\
\[Infinity], \[Infinity], \[Infinity], \[Infinity], \[Infinity]}, \
{\[Infinity], \[Infinity], \[Infinity], \[Infinity], \[Infinity]}, {\
\[Infinity], \[Infinity], \[Infinity], \[Infinity], \[Infinity]}, \
{\[Infinity], \[Infinity], \[Infinity], \[Infinity], \[Infinity]}}\)], \
"Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(TryUsingRMatrixRelations[lo0_, hi0_, lo1_, hi1_, lo2_, hi2_, i_, 
        j_]\  := \[IndentingNewLine]Module[{infinities, \ loToFind, \ 
          hiToFind, l, \ x, \ eqn, \ solutions, R1, \ R2, \ 
          R3}, \[IndentingNewLine]x =. ; \[IndentingNewLine]l = {{lo0, 
              hi0}, {lo1, hi1}, {lo2, hi2}}; \[IndentingNewLine]l = 
          Sort[l]; \ \[IndentingNewLine] (*\ Hopefully, \ 
          the\ relations\ are\ now\ \
R\_\(\(L[\([1]\)]\)[\([1]\)]\)\%\(\(L[\([1]\)]\)[\([2]\)]\) + 
              R\_\(\(L[\([3]\)]\)[\([1]\)]\)\%\(\(L[\([3]\)]\)[\([2]\)]\)\  \
== \ R\_\(\(L[\([2]\)]\)[\([1]\)]\)\%\(\(L[\([2]\)]\)[\([2]\)]\) . \ 
                  Imagine\ this\ as\ R\_12\%14 + R\_15\%16\  \[Equal] \ 
            R\_12\%16*) \[IndentingNewLine]infinities = 0; loToFind = 0; 
        hiToFind = 0; \[IndentingNewLine]Scan[
          Function[\[IndentingNewLine]\[IndentingNewLine]If[\(\(R[#[\([1]\)], \
#[\([2]\)]]\)[\([i]\)]\)[\([j]\)] \[Equal] Infinity, \ \((\(++infinities\); 
                loToFind\  = \ #[\([1]\)]; \ hiToFind\  = \ #[\([2]\)]; \ 
                Print[loToFind, "\< \>", \ 
                  hiToFind])\)]\[IndentingNewLine]], \ 
          l]; \[IndentingNewLine]Print[loToFind, \ "\< \>", \ 
          hiToFind, \ "\< \>", \ i, \ "\< \>", \ j]; \[IndentingNewLine]If[
          infinities\  \[Equal] \ 
            1, \[IndentingNewLine]\(\(R[loToFind, 
                  hiToFind]\)[\([i]\)]\)[\([j]\)]\ \  = 
            x; \[IndentingNewLine]R1\  := \ 
            R[\(l[\([1]\)]\)[\([1]\)], \ \(l[\([1]\)]\)[\([2]\)]]; \ 
          R2 := R[\(l[\([3]\)]\)[\([1]\)], \ \(l[\([3]\)]\)[\([2]\)]]; \ 
          R3 := R[\(l[\([2]\)]\)[\([1]\)], \ \(l[\([2]\)]\)[\([2]\)]]; \
\[IndentingNewLine]eqn\  = \ \((\ \(R1[\([i]\)]\)[\([j]\)]\  + \ \
\(R2[\([i]\)]\)[\([j]\)]\  \[Equal] \
\[IndentingNewLine]\(R3[\([i]\)]\)[\([j]\)])\); \[IndentingNewLine]solutions\ \
 = \ Solve[eqn, \ x]; \[IndentingNewLine]If[
            Length[solutions] \[NotEqual] 1, \ 
            Message[recover::TooManySolutions, \ "\<Multiple or no solutions \
found where there should only be one\>"]]; \[IndentingNewLine] (*\ 
            Substitute\ the\ solution\ in . \ 
                The\ weird\ syntax\ comes\ from\ the\ wierd\ output\ of\ \
Solve\ *) \[IndentingNewLine]\(\(R[loToFind, 
                  hiToFind]\)[\([i]\)]\)[\([j]\)]\  = \ 
            x\  /. \ First[
                solutions];\[IndentingNewLine]]\[IndentingNewLine]]\), "\
\[IndentingNewLine]", 
    \(True\), "\[IndentingNewLine]", 
    \(\)}], "Input"],

Cell[BoxData[
    \(True\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(R[1, 2]\), "\[IndentingNewLine]", 
    \(R[1, 2] = {{Infinity, \ 1}, \ {2, \ 3}}\), "\[IndentingNewLine]", 
    \(R[3, 
        5] = {{4, \ Infinity}, \ {5, \ Infinity}}\), "\[IndentingNewLine]", 
    \(R[1, 5] = {{7, 8}, \ {Infinity, \ Infinity}}\)}], "Input"],

Cell[BoxData[
    \({{\[Infinity], 1}, {2, 3}}\)], "Output"],

Cell[BoxData[
    \({{\[Infinity], 1}, {2, 3}}\)], "Output"],

Cell[BoxData[
    \({{4, \[Infinity]}, {5, \[Infinity]}}\)], "Output"],

Cell[BoxData[
    \({{7, 8}, {\[Infinity], \[Infinity]}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Attributes[Superscript]\)], "Input"],

Cell[BoxData[
    \({}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData[
    \({\[Infinity], 1}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Block[{i, j, l, h}\ , \ i = 1; j = 1; 
      Print\ ["\<|\>", \(\(R[1, 2]\)[\([i]\)]\)[\([j]\)]\  = \ 
          5, "\<|\>"]; \(\(R[1, 2]\)[\([i]\)]\)[\([j]\)]\ \  = x]\)], "Input"],

Cell[BoxData[
    RowBox[{\(Set::"setps"\), \(\(:\)\(\ \)\), "\<\"\\!\\(\\(\\(R[\\(\\(1, \
2\\)\\)]\\)\\) \[LeftDoubleBracket] i \[RightDoubleBracket]\\) in assignment \
of part is not a symbol. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::setps\\\"]\\)\"\>"}]], "Message"],

Cell[BoxData[
    InterpretationBox[\("|"\[InvisibleSpace]5\[InvisibleSpace]"|"\),
      SequenceForm[ "|", 5, "|"],
      Editable->False]], "Print"],

Cell[BoxData[
    RowBox[{\(Set::"setps"\), \(\(:\)\(\ \)\), "\<\"\\!\\(\\(\\(R[\\(\\(1, \
2\\)\\)]\\)\\) \[LeftDoubleBracket] i \[RightDoubleBracket]\\) in assignment \
of part is not a symbol. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::setps\\\"]\\)\"\>"}]], "Message"],

Cell[BoxData[
    \(x\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(Trace[
        TryUsingRMatrixRelations[1, 2, 3, 5, 1, 5, 1, 
          1]]\)[\([3]\)]\)], "Input"],

Cell[BoxData[
    InterpretationBox[\(1\[InvisibleSpace]" "\[InvisibleSpace]2\),
      SequenceForm[ 1, " ", 2],
      Editable->False]], "Print"],

Cell[BoxData[
    InterpretationBox[\(1\[InvisibleSpace]" "\[InvisibleSpace]2\
\[InvisibleSpace]" "\[InvisibleSpace]1\[InvisibleSpace]" \
"\[InvisibleSpace]1\),
      SequenceForm[ 1, " ", 2, " ", 1, " ", 1],
      Editable->False]], "Print"],

Cell[BoxData[
    RowBox[{\(Set::"setps"\), \(\(:\)\(\ \)\), \
"\<\"\\!\\(\\(\\(R[\\(\\(loToFind$291, hiToFind$291\\)\\)]\\)\\) \
\[LeftDoubleBracket] 1 \[RightDoubleBracket]\\) in assignment of part is not \
a symbol. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::setps\\\"]\\)\"\>"}]], "Message"],

Cell[BoxData[
    \(recover::"TooManySolutions" \(\(:\)\(\ \)\) 
      "-- Message text not found -- (\!\(\"Multiple or no solutions found \
where there should only be one\"\))"\)], "Message"],

Cell[BoxData[
    RowBox[{\(First::"first"\), \(\(:\)\(\ \)\), "\<\"\\!\\({}\\) has a \
length of zero and no first element. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\
\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"First::first\\\"]\\)\"\>"}]], "Message"],

Cell[BoxData[
    RowBox[{\(ReplaceAll::"reps"\), \(\(:\)\(\ \)\), "\<\"\\!\\({\\(First[\\(\
\\({}\\)\\)]\\)}\\) is neither a list of replacement rules nor a valid \
dispatch table, and so cannot be used for replacing. \
\\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\
\\\", ButtonFrame->None, ButtonData:>\\\"ReplaceAll::reps\\\"]\\)\"\>"}]], \
"Message"],

Cell[BoxData[
    RowBox[{\(Set::"setps"\), \(\(:\)\(\ \)\), \
"\<\"\\!\\(\\(\\(R[\\(\\(loToFind$291, hiToFind$291\\)\\)]\\)\\) \
\[LeftDoubleBracket] 1 \[RightDoubleBracket]\\) in assignment of part is not \
a symbol. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::setps\\\"]\\)\"\>"}]], "Message"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[""], "Input"],

Cell[BoxData[
    RowBox[{\(Syntax::"tsntxi"\), \(\(:\)\(\ \)\), "\<\"\\\"\\!\\(//TableForm\
\\)\\\" is incomplete; more input is needed.\\!\\(\\\"\\\"\\) \
\\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\
\\\", ButtonFrame->None, ButtonData:>\\\"General::tsntxi\\\"]\\)\"\>"}]], \
"Message"],

Cell[BoxData[
    StyleBox[
      ErrorBox[\(\(//\)\(TableForm\)\)],
      ShowStringCharacters->True]], "Message"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    RowBox[{"MatrixForm", "[", 
      InterpretationBox[\(R\_1\%2\),
        Subsuperscript[ R, 1, 2],
        Editable->False], "]"}], "\n", 
    RowBox[{"MatrixForm", "[", 
      InterpretationBox[\(R\_3\%5\),
        Subsuperscript[ R, 3, 5],
        Editable->False], "]"}], "\n", 
    RowBox[{"MatrixForm", "[", 
      InterpretationBox[\(R\_1\%5\),
        Subsuperscript[ R, 1, 5],
        Editable->False], "]"}]}], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {"\[Infinity]", "1"},
            {"2", "3"}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {"4", "\[Infinity]"},
            {"5", "\[Infinity]"}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {"7", "8"},
            {"\[Infinity]", "\[Infinity]"}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(R = {}\)\(\[IndentingNewLine]\)
    \)\)], "Input"],

Cell[BoxData[
    \({}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(R[\([5]\)]\  = \ 5;\)\)], "Input"],

Cell[BoxData[
    RowBox[{\(Set::"partw"\), \(\(:\)\(\ \)\), "\<\"Part \\!\\(5\\) of \
\\!\\({}\\) does not exist. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::partw\\\"]\\)\"\>"}]], "Message"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(R = SparseArray[{1 \[Rule] 0}]\)], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"SparseArray", "[", 
        RowBox[{
          InterpretationBox[\("<"\[InvisibleSpace]0\[InvisibleSpace]">"\),
            SequenceForm[ "<", 0, ">"],
            Editable->False], ",", \({1}\)}], "]"}],
      False,
      Editable->False]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(R[\([5]\)]\)], "Input"],

Cell[BoxData[
    RowBox[{\(Part::"partw"\), \(\(:\)\(\ \)\), "\<\"Part \\!\\(5\\) of \
\\!\\(\\(\\(\\*TagBox[\\(SparseArray[\\(\\(Automatic, \\(\\({1}\\)\\), 0, \\(\
\\({1, \\(\\({\\(\\({0, 0}\\)\\), \\(\\({}\\)\\)}\\)\\), \
\\(\\({}\\)\\)}\\)\\)\\)\\)]\\), False, Rule[Editable, False]]\\)\\) \
\[LeftDoubleBracket] 5 \[RightDoubleBracket]\\) does not exist. \
\\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\
\\\", ButtonFrame->None, ButtonData:>\\\"General::partw\\\"]\\)\"\>"}]], \
"Message"],

Cell[BoxData[
    RowBox[{\(Part::"partw"\), \(\(:\)\(\ \)\), "\<\"Part \\!\\(5\\) of \
\\!\\({0}\\) does not exist. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::partw\\\"]\\)\"\>"}]], "Message"],

Cell[BoxData[
    RowBox[{
      TagBox[
        RowBox[{"SparseArray", "[", 
          RowBox[{
            InterpretationBox[\("<"\[InvisibleSpace]0\[InvisibleSpace]">"\),
              SequenceForm[ "<", 0, ">"],
              Editable->False], ",", \({1}\)}], "]"}],
        False,
        Editable->False], "\[LeftDoubleBracket]", "5", 
      "\[RightDoubleBracket]"}]], "Output"]
}, Open  ]],

Cell[BoxData[
    TagBox[\(SparseArray[{}]\),
      False,
      Editable->False]], "Input"],

Cell[BoxData[
    TagBox[\(SparseArray[]\),
      False,
      Editable->False]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
    \(R$1$5 + 1\)], "Input"],

Cell[BoxData[
    \(1 + R$1$5\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(R$1$5 = 5\)], "Input"],

Cell[BoxData[
    \(5\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(R$1$5\)], "Input"],

Cell[BoxData[
    \(5\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Aargh[1, 2]\  = \ {{5, 6}, {7, 8}}\)], "Input"],

Cell[BoxData[
    \({{5, 6}, {7, 8}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData[
    \(6\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Remove[R]\)], "Input"],

Cell[BoxData[
    \({\[Lambda]\_\(1, 1\) \[Rule] \(-\(31\/20\)\), \[Lambda]\_\(1, 2\) \
\[Rule] 1\/4, \[Lambda]\_\(2, 2\) \[Rule] \(-\(3\/2\)\), \[Lambda]\_\(1, 3\) \
\[Rule] 1\/5, \[Lambda]\_\(2, 3\) \[Rule] 
        1\/4, \[Lambda]\_\(3, 3\) \[Rule] \(-\(31\/20\)\), \[Lambda]\_\(1, \
4\) \[Rule] 9\/20, \[Lambda]\_\(2, 4\) \[Rule] 
        1\/4, \[Lambda]\_\(3, 4\) \[Rule] 
        1\/5, \[Lambda]\_\(4, 4\) \[Rule] \(-\(31\/20\)\), \[Lambda]\_\(1, \
5\) \[Rule] 1\/4, \[Lambda]\_\(2, 5\) \[Rule] 
        1\/2, \[Lambda]\_\(3, 5\) \[Rule] 1\/4, \[Lambda]\_\(4, 5\) \[Rule] 
        1\/4, \[Lambda]\_\(5, 5\) \[Rule] \(-\(3\/2\)\), \[Lambda]\_\(1, 6\) \
\[Rule] 2\/5, \[Lambda]\_\(2, 6\) \[Rule] 1\/4, \[Lambda]\_\(3, 6\) \[Rule] 
        13\/20, \[Lambda]\_\(4, 6\) \[Rule] 2\/5, \[Lambda]\_\(5, 6\) \[Rule] 
        1\/4, \[Lambda]\_\(6, 6\) \[Rule] \(-\(39\/20\)\)}\)], "Output"]
}, Open  ]],

Cell[BoxData[""], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
    \(L1 = \(L2 = Table[0, {9}, {9}]\)\)], "Input"],

Cell[BoxData[
    \({{0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 
        0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 
        0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 
        0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\(\(With[\[IndentingNewLine]{i = 6, j = 9, 
        k = 23/13}, \[IndentingNewLine]L1[\([i, j]\)] = 1; 
      L2[\([i, j]\)] = k]\)\(\[IndentingNewLine]\)
    \)\)], "Input"],

Cell[BoxData[
    \(23\/13\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(L1 = {{0, 0, 0, 0, 0, 0, 1, 0, 1}, {0, 0, 0, 0, 0, 0, 1, 1, 0}, {0, 0, 
          0, 0, 0, 0, 0, 1, 1}, {0, 0, 0, 0, 0, 0, 1, 0, 1}, {0, 0, 0, 0, 0, 
          0, 1, 1, 0}, {0, 0, 0, 0, 0, 0, 0, 1, 2}, {0, 0, 0, 0, 0, 0, 0, 0, 
          0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 
          0}}\)], "Input"],

Cell[BoxData[
    \({{0, 0, 0, 0, 0, 0, 1, 0, 1}, {0, 0, 0, 0, 0, 0, 1, 1, 0}, {0, 0, 0, 0, 
        0, 0, 0, 1, 1}, {0, 0, 0, 0, 0, 0, 1, 0, 1}, {0, 0, 0, 0, 0, 0, 1, 1, 
        0}, {0, 0, 0, 0, 0, 0, 0, 1, 2}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 
        0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(L2 = {{0, 0, 0, 0, 0, 0, 19\/21, 0, 23\/21}, {0, 0, 0, 0, 0, 0, 19\/17, 
          15\/17, 0}, {0, 0, 0, 0, 0, 0, 0, 15\/13, 23\/26}, {0, 0, 0, 0, 0, 
          0, 19\/21, 0, 23\/21}, {0, 0, 0, 0, 0, 0, 19\/17, 15\/17, 0}, {0, 
          0, 0, 0, 0, 0, 0, 15\/13, 23\/13}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 
          0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}}\)], "Input"],

Cell[BoxData[
    \({{0, 0, 0, 0, 0, 0, 19\/21, 0, 23\/21}, {0, 0, 0, 0, 0, 0, 19\/17, 
        15\/17, 0}, {0, 0, 0, 0, 0, 0, 0, 15\/13, 23\/26}, {0, 0, 0, 0, 0, 0, 
        19\/21, 0, 23\/21}, {0, 0, 0, 0, 0, 0, 19\/17, 15\/17, 0}, {0, 0, 0, 
        0, 0, 0, 0, 15\/13, 23\/13}, {0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 
        0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(L1 = L1 + Transpose[L1];\)\), "\[IndentingNewLine]", 
    \(\(Function[i, L1[\([i, i]\)]\  = \ \(-Total[L1[\([i]\)]]\)]\  /@ \ 
        Range[9];\)\), "\[IndentingNewLine]", 
    \(L1\  // MatrixForm\)}], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {\(-2\), "0", "0", "0", "0", "0", "1", "0", "1"},
            {"0", \(-2\), "0", "0", "0", "0", "1", "1", "0"},
            {"0", "0", \(-2\), "0", "0", "0", "0", "1", "1"},
            {"0", "0", "0", \(-2\), "0", "0", "1", "0", "1"},
            {"0", "0", "0", "0", \(-2\), "0", "1", "1", "0"},
            {"0", "0", "0", "0", "0", \(-3\), "0", "1", "2"},
            {"1", "1", "0", "1", "1", "0", \(-4\), "0", "0"},
            {"0", "1", "1", "0", "1", "1", "0", \(-4\), "0"},
            {"1", "0", "1", "1", "0", "2", "0", "0", \(-5\)}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(L1\)], "Input"],

Cell[BoxData[
    \({{\(-2\), 0, 0, 0, 0, 0, 1, 0, 1}, {0, \(-2\), 0, 0, 0, 0, 1, 1, 
        0}, {0, 0, \(-2\), 0, 0, 0, 0, 1, 1}, {0, 0, 0, \(-2\), 0, 0, 1, 0, 
        1}, {0, 0, 0, 0, \(-2\), 0, 1, 1, 0}, {0, 0, 0, 0, 0, \(-3\), 0, 1, 
        2}, {1, 1, 0, 1, 1, 0, \(-4\), 0, 0}, {0, 1, 1, 0, 1, 1, 0, \(-4\), 
        0}, {1, 0, 1, 1, 0, 2, 0, 0, \(-5\)}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(L2 = L2 + Transpose[L2];\)\), "\[IndentingNewLine]", 
    \(\(Function[i, L2[\([i, i]\)]\  = \ \(-Total[L2[\([i]\)]]\)]\  /@ \ 
        Range[9];\)\), "\[IndentingNewLine]", 
    \(L2\  // MatrixForm\)}], "Input"],

Cell[BoxData[
    TagBox[
      RowBox[{"(", "\[NoBreak]", GridBox[{
            {\(-2\), "0", "0", "0", "0", "0", \(19\/21\), "0", \(23\/21\)},
            {"0", \(-2\), "0", "0", "0", "0", \(19\/17\), \(15\/17\), "0"},
            {"0", "0", \(-\(53\/26\)\), "0", "0", "0", 
              "0", \(15\/13\), \(23\/26\)},
            {"0", "0", "0", \(-2\), "0", "0", \(19\/21\), "0", \(23\/21\)},
            {"0", "0", "0", "0", \(-2\), "0", \(19\/17\), \(15\/17\), "0"},
            {"0", "0", "0", "0", "0", \(-\(38\/13\)\), 
              "0", \(15\/13\), \(23\/13\)},
            {\(19\/21\), \(19\/17\), "0", \(19\/21\), \(19\/17\), 
              "0", \(-\(1444\/357\)\), "0", "0"},
            {"0", \(15\/17\), \(15\/13\), "0", \(15\/17\), \(15\/13\), 
              "0", \(-\(900\/221\)\), "0"},
            {\(23\/21\), "0", \(23\/26\), \(23\/21\), "0", \(23\/13\), "0", 
              "0", \(-\(2645\/546\)\)}
            }], "\[NoBreak]", ")"}],
      Function[ BoxForm`e$, 
        MatrixForm[ BoxForm`e$]]]], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(N[Eigensystem[L2]]\)], "Input"],

Cell[BoxData[
    \({{\(-6.806824912338198`\), \(-5.244266937465268`\), \
\(-4.830726498152763`\), \(-2.4097516295719803`\), \(-2.`\), \(-2.`\), \
\(-1.3699838091374188`\), \(-1.2615231364112953`\), 
        0.`}, {{\(-0.28845616081343906`\), \(-0.16434426003286262`\), \
\(-0.3034719024108199`\), \(-0.28845616081343906`\), \
\(-0.16434426003286262`\), \(-0.6003684316110593`\), 0.3219854451881474`, 
          0.48745573053111`, 1.`}, {0.1015784901594543`, 0.7627808608517341`, 
          0.015567859990596844`, 0.1015784901594543`, 
          0.7627808608517341`, \(-0.35960341455620437`\), \
\(-1.5747632884613267`\), \(-0.8099198589934531`\), 
          1.`}, {\(-1.0067751750559637`\), 0.13314446300029203`, 
          0.8748131947986622`, \(-1.0067751750559637`\), 
          0.13314446300029203`, 0.8167615613743919`, 
          1.939368867367648`, \(-2.883682199428449`\), 
          1.`}, {\(-1.1122209815897222`\), 
          0.3833468329828884`, \(-4.611627271148372`\), \
\(-1.1122209815897222`\), 0.3833468329828884`, 
          5.05891084492697`, \(-0.7068199182808996`\), 0.7172846417157004`, 
          1.`}, {0.`, \(-1.`\), 0.`, 0.`, 1.`, 0.`, 0.`, 0.`, 
          0.`}, {\(-1.`\), 0.`, 0.`, 1.`, 0.`, 0.`, 0.`, 0.`, 
          0.`}, {4.223110187058179`, \(-1.348326910000551`\), \
\(-4.121204667001631`\), 
          4.223110187058179`, \(-1.348326910000551`\), \(-1.2042537156820738`\
\), 1.730167561395302`, \(-3.154275732826846`\), 
          1.`}, {0.31952085488999105`, \(-1.441223575207896`\), 
          1.1337985810477849`, 0.31952085488999105`, \(-1.441223575207896`\), 
          1.0625638824274404`, \(-0.9497297929379962`\), \
\(-0.003227229901402451`\), 1.`}, {1.`, 1.`, 1.`, 1.`, 1.`, 1.`, 1.`, 1.`, 
          1.`}}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(N[Eigensystem[L1]]\)], "Input"],

Cell[BoxData[
    \({{\(-6.9777919682209655`\), \(-5.266482549405689`\), \
\(-4.732050807568877`\), \(-2.388667917140328`\), \(-2.`\), \(-2.`\), \
\(-1.3670575652330175`\), \(-1.2679491924311228`\), 
        0.`}, {{\(-0.2521074119496802`\), \(-0.1274691251661224`\), \
\(-0.2771462821112749`\), \(-0.2521074119496802`\), \(-0.1274691251661224`\), \
\(-0.5982154311051318`\), 0.2549382503322448`, 0.37957653711583816`, 
          1.`}, {0.23592142848814035`, 0.8853166145936804`, 
          0.03711577859305493`, 0.23592142848814035`, 
          0.8853166145936804`, \(-0.38772059248752555`\), \
\(-1.7706332291873608`\), \(-1.1212380430818456`\), 
          1.`}, {\(-1.3660254037844384`\), 0.3660254037844384`, 
          1.`, \(-1.3660254037844384`\), 0.3660254037844384`, 1.`, 
          2.732050807568877`, \(-3.732050807568877`\), 
          1.`}, {\(-0.8661304090090027`\), 
          0.33168144897928853`, \(-3.9479691849010985`\), \
\(-0.8661304090090027`\), 0.33168144897928853`, 
          4.14578104288938`, \(-0.6633628979585771`\), 0.5344489600297191`, 
          1.`}, {0.`, \(-1.`\), 0.`, 0.`, 1.`, 0.`, 0.`, 0.`, 
          0.`}, {\(-1.`\), 0.`, 0.`, 1.`, 0.`, 0.`, 0.`, 0.`, 
          0.`}, {12.132316392470566`, \(-3.3395289384068474`\), \
\(-12.312000311580675`\), 
          12.132316392470566`, \(-3.3395289384068474`\), \
\(-4.159845019296743`\), 6.679057876813695`, \(-8.792787454063724`\), 
          1.`}, {0.3660254037844386`, \(-1.3660254037844386`\), 1.`, 
          0.3660254037844386`, \(-1.3660254037844386`\), 
          1.`, \(-0.7320508075688772`\), \(-0.2679491924311228`\), 
          1.`}, {1.`, 1.`, 1.`, 1.`, 1.`, 1.`, 1.`, 1.`, 1.`}}}\)], "Output"]
}, Open  ]]
},
FrontEndVersion->"5.0 for Microsoft Windows",
ScreenRectangle->{{0, 1397}, {0, 918}},
WindowSize->{634, 878},
WindowMargins->{{234, Automatic}, {Automatic, 0}},
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->{
  "Info3360852581-7681982"->{
    Cell[13788, 362, 447, 10, 63, "Print",
      CellTags->"Info3360852581-7681982"],
    Cell[14238, 374, 790, 19, 99, "Print",
      CellTags->"Info3360852581-7681982"]}
  }
*)

(*CellTagsIndex
CellTagsIndex->{
  {"Info3360852581-7681982", 40327, 1088}
  }
*)

(*NotebookFileOutline
Notebook[{
Cell[1754, 51, 307, 7, 130, "Input"],
Cell[2064, 60, 2924, 55, 750, "Input"],
Cell[4991, 117, 1010, 22, 270, "Input"],

Cell[CellGroupData[{
Cell[6026, 143, 390, 8, 90, "Input"],
Cell[6419, 153, 737, 14, 153, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[7193, 172, 59, 1, 30, "Input"],
Cell[7255, 175, 375, 6, 110, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[7667, 186, 49, 1, 30, "Input"],
Cell[7719, 189, 886, 13, 110, "Output"]
}, Open  ]],
Cell[8620, 205, 41, 1, 30, "Input"],

Cell[CellGroupData[{
Cell[8686, 210, 217, 4, 70, "Input"],
Cell[8906, 216, 275, 8, 57, "Output"],
Cell[9184, 226, 262, 4, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[9483, 235, 71, 1, 30, "Input"],
Cell[9557, 238, 1127, 15, 88, "Message"],
Cell[10687, 255, 973, 13, 130, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[11697, 273, 34, 1, 30, "Input"],
Cell[11734, 276, 306, 6, 50, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[12077, 287, 52, 1, 30, "Input"],
Cell[12132, 290, 40, 1, 29, "Output"]
}, Open  ]],
Cell[12187, 294, 303, 6, 70, "Input"],

Cell[CellGroupData[{
Cell[12515, 304, 42, 1, 30, "Input"],
Cell[12560, 307, 35, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[12632, 313, 62, 1, 30, "Input"],
Cell[12697, 316, 39, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[12773, 322, 61, 1, 30, "Input"],
Cell[12837, 325, 680, 20, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[13554, 350, 90, 1, 30, "Input"],
Cell[13647, 353, 60, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[13744, 359, 41, 1, 30, "Input"],
Cell[13788, 362, 447, 10, 63, "Print",
  CellTags->"Info3360852581-7681982"],
Cell[14238, 374, 790, 19, 99, "Print",
  CellTags->"Info3360852581-7681982"]
}, Open  ]],
Cell[15043, 396, 105, 2, 30, "Input"],
Cell[15151, 400, 39, 1, 30, "Input"],

Cell[CellGroupData[{
Cell[15215, 405, 62, 1, 30, "Input"],
Cell[15280, 408, 45, 1, 29, "Output"]
}, Open  ]],
Cell[15340, 412, 72, 1, 30, "Input"],

Cell[CellGroupData[{
Cell[15437, 417, 294, 5, 110, "Input"],
Cell[15734, 424, 977, 20, 89, "Output"],
Cell[16714, 446, 978, 20, 89, "Output"],
Cell[17695, 468, 224, 3, 48, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[17956, 476, 55, 1, 30, "Input"],
Cell[18014, 479, 379, 6, 48, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[18430, 490, 2571, 43, 676, "Input"],
Cell[21004, 535, 38, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[21079, 541, 282, 5, 90, "Input"],
Cell[21364, 548, 60, 1, 29, "Output"],
Cell[21427, 551, 60, 1, 29, "Output"],
Cell[21490, 554, 70, 1, 29, "Output"],
Cell[21563, 557, 70, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[21670, 563, 56, 1, 30, "Input"],
Cell[21729, 566, 36, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[21802, 572, 53, 1, 30, "Input"],
Cell[21858, 575, 50, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[21945, 581, 201, 3, 50, "Input"],
Cell[22149, 586, 359, 5, 22, "Message"],
Cell[22511, 593, 150, 3, 25, "Print"],
Cell[22664, 598, 359, 5, 22, "Message"],
Cell[23026, 605, 35, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[23098, 611, 121, 3, 30, "Input"],
Cell[23222, 616, 146, 3, 25, "Print"],
Cell[23371, 621, 242, 5, 25, "Print"],
Cell[23616, 628, 383, 6, 37, "Message"],
Cell[24002, 636, 192, 3, 37, "Message"],
Cell[24197, 641, 287, 4, 22, "Message"],
Cell[24487, 647, 387, 6, 37, "Message"],
Cell[24877, 655, 383, 6, 37, "Message"]
}, Open  ]],

Cell[CellGroupData[{
Cell[25297, 666, 26, 0, 30, "Input"],
Cell[25326, 668, 322, 5, 22, "Message"],
Cell[25651, 675, 115, 3, 22, "Message"]
}, Open  ]],

Cell[CellGroupData[{
Cell[25803, 683, 447, 12, 73, "Input"],
Cell[26253, 697, 237, 7, 41, "Output"],
Cell[26493, 706, 247, 7, 41, "Output"],
Cell[26743, 715, 247, 7, 41, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[27027, 727, 75, 2, 50, "Input"],
Cell[27105, 731, 36, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[27178, 737, 56, 1, 30, "Input"],
Cell[27237, 740, 279, 4, 22, "Message"]
}, Open  ]],

Cell[CellGroupData[{
Cell[27553, 749, 63, 1, 30, "Input"],
Cell[27619, 752, 293, 8, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[27949, 765, 43, 1, 30, "Input"],
Cell[27995, 768, 527, 8, 37, "Message"],
Cell[28525, 778, 281, 4, 22, "Message"],
Cell[28809, 784, 385, 10, 29, "Output"]
}, Open  ]],
Cell[29209, 797, 92, 3, 30, "Input"],
Cell[29304, 802, 90, 3, 30, "Input"],

Cell[CellGroupData[{
Cell[29419, 809, 42, 1, 30, "Input"],
Cell[29464, 812, 43, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[29544, 818, 42, 1, 30, "Input"],
Cell[29589, 821, 35, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[29661, 827, 38, 1, 30, "Input"],
Cell[29702, 830, 35, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[29774, 836, 67, 1, 30, "Input"],
Cell[29844, 839, 50, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[29931, 845, 70, 1, 30, "Input"],
Cell[30004, 848, 35, 1, 29, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[30076, 854, 42, 1, 30, "Input"],
Cell[30121, 857, 886, 13, 110, "Output"]
}, Open  ]],
Cell[31022, 873, 26, 0, 30, "Input"],

Cell[CellGroupData[{
Cell[31073, 877, 65, 1, 30, "Input"],
Cell[31141, 880, 322, 4, 105, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[31500, 889, 194, 4, 90, "Input"],
Cell[31697, 895, 40, 1, 42, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[31774, 901, 343, 5, 110, "Input"],
Cell[32120, 908, 322, 4, 105, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[32479, 917, 403, 5, 152, "Input"],
Cell[32885, 924, 391, 5, 169, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[33313, 934, 236, 4, 70, "Input"],
Cell[33552, 940, 737, 14, 153, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[34326, 959, 35, 1, 30, "Input"],
Cell[34364, 962, 376, 5, 105, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[34777, 972, 236, 4, 70, "Input"],
Cell[35016, 978, 1033, 19, 233, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[36086, 1002, 51, 1, 30, "Input"],
Cell[36140, 1005, 1764, 29, 219, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[37941, 1039, 51, 1, 30, "Input"],
Cell[37995, 1042, 1690, 27, 219, "Output"]
}, Open  ]]
}
]
*)



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