(* 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[     51807,       1659]

NotebookOptionsPosition[     49044,       1566]

NotebookOutlinePosition[     49407,       1582]

CellTagsIndexPosition[     49364,       1579]

WindowFrame->Normal

ContainsDynamic->False*)

(* Beginning of Notebook Content *)

Notebook[{

Cell[BoxData[

 RowBox[{"Clear", "[", 

  RowBox[{"d", ",", "tm", ",", "m"}], "]"}]], "Input",

 CellChangeTimes->{{3.4241980742571344`*^9, 3.424198074928099*^9}, {

  3.424198111630875*^9, 3.424198112381955*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"tm", "[", "x_", "]"}], ":=", "x"}]], "Input",

 CellChangeTimes->{{3.424198093104235*^9, 3.4241981152961454`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"m", "[", "x_", "]"}], ":=", 

  RowBox[{"d", "[", "x", "]"}]}]], "Input",

 CellChangeTimes->{{3.424198093104235*^9, 3.4241981152961454`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[

 RowBox[{"d", "[", 

  RowBox[{"tm", "[", 

   RowBox[{"a", ",", "b"}], "]"}], "]"}]], "Input",

 CellChangeTimes->{{3.424197939102792*^9, 3.4241979434390273`*^9}, {

  3.424198082939619*^9, 3.424198084061232*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"m", "[", 

   RowBox[{"a", ",", "b"}], "]"}], "+", 

  RowBox[{"\[Sigma]", " ", 

   RowBox[{"t", "[", 

    RowBox[{"a", ",", "b"}], "]"}]}], "+", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], 

    RowBox[{"1", "+", 

     RowBox[{"deg", "[", "a", "]"}]}]], " ", 

   RowBox[{"tm", "[", 

    RowBox[{"a", ",", 

     RowBox[{"d", "[", "b", "]"}]}], "]"}]}], "-", 

  RowBox[{"tm", "[", 

   RowBox[{

    RowBox[{"d", "[", "a", "]"}], ",", "b"}], "]"}]}]], "Output",

 CellChangeTimes->{

  3.424197943959776*^9, {3.424198070712037*^9, 3.4241981177697024`*^9}, 

   3.4241981527399874`*^9, {3.424198220988123*^9, 3.4241982242928753`*^9}, {

   3.4241984914269953`*^9, 3.424198497575837*^9}, {3.4241986744701986`*^9, 

   3.424198684594757*^9}, 3.4241987900664177`*^9, {3.4241989830138626`*^9, 

   3.4241990466954317`*^9}}]

}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[

 RowBox[{"d", "[", 

  RowBox[{"tm", "[", 

   RowBox[{"a", ",", "b", ",", "c"}], "]"}], "]"}]], "Input",

 CellChangeTimes->{{3.424198928445397*^9, 3.4241989306886225`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"m", "[", 

   RowBox[{"a", ",", "b", ",", "c"}], "]"}], "+", 

  RowBox[{"\[Sigma]", " ", 

   RowBox[{"t", "[", 

    RowBox[{"a", ",", 

     RowBox[{"tm", "[", 

      RowBox[{"b", ",", "c"}], "]"}]}], "]"}]}], "+", 

  RowBox[{"\[Sigma]", " ", 

   RowBox[{"t", "[", 

    RowBox[{

     RowBox[{"tm", "[", 

      RowBox[{"a", ",", "b"}], "]"}], ",", "c"}], "]"}]}], "+", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], 

    RowBox[{"1", "+", 

     RowBox[{"deg", "[", "a", "]"}]}]], " ", 

   RowBox[{"tm", "[", 

    RowBox[{"a", ",", 

     RowBox[{"m", "[", 

      RowBox[{"b", ",", "c"}], "]"}]}], "]"}]}], "-", 

  RowBox[{"tm", "[", 

   RowBox[{

    RowBox[{"m", "[", 

     RowBox[{"a", ",", "b"}], "]"}], ",", "c"}], "]"}], "+", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], 

    RowBox[{"2", "+", 

     RowBox[{"deg", "[", "a", "]"}], "+", 

     RowBox[{"deg", "[", "b", "]"}]}]], " ", 

   RowBox[{"tm", "[", 

    RowBox[{"a", ",", "b", ",", 

     RowBox[{"d", "[", "c", "]"}]}], "]"}]}], "+", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], 

    RowBox[{"2", "+", 

     RowBox[{"deg", "[", "a", "]"}]}]], " ", 

   RowBox[{"tm", "[", 

    RowBox[{"a", ",", 

     RowBox[{"d", "[", "b", "]"}], ",", "c"}], "]"}]}], "+", 

  RowBox[{"tm", "[", 

   RowBox[{

    RowBox[{"d", "[", "a", "]"}], ",", "b", ",", "c"}], "]"}]}]], "Output",

 CellChangeTimes->{3.4241989309990687`*^9, 3.424199050130371*^9}]

}, Open  ]],

Cell[BoxData[

 RowBox[{

  RowBox[{"d", "[", 

   RowBox[{"a", ":", 

    RowBox[{"t", "[", "___", "]"}]}], "]"}], ":=", 

  RowBox[{"With", "[", 

   RowBox[{

    RowBox[{"{", 

     RowBox[{"k", "=", 

      RowBox[{"Length", "[", "a", "]"}]}], "}"}], ",", 

    RowBox[{

     UnderoverscriptBox["\[Sum]", 

      RowBox[{"i", "=", "1"}], "k"], 

     RowBox[{

      SuperscriptBox[

       RowBox[{"(", 

        RowBox[{"-", "1"}], ")"}], 

       RowBox[{

        UnderoverscriptBox["\[Sum]", 

         RowBox[{"j", "=", "1"}], 

         RowBox[{"i", "-", "1"}]], 

        RowBox[{"deg", "[", 

         RowBox[{"a", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}], 

         "]"}]}]], 

      RowBox[{"MapAt", "[", 

       RowBox[{"d", ",", "a", ",", "i"}], "]"}]}]}]}], "]"}]}]], "Input",

 CellChangeTimes->{{3.4241985045058017`*^9, 3.4241986214940224`*^9}, {

  3.4241986537403903`*^9, 3.424198663494416*^9}}],

Cell[CellGroupData[{

Cell[BoxData[

 RowBox[{"d", "[", 

  RowBox[{"t", "[", 

   RowBox[{"a", ",", "b", ",", "c"}], "]"}], "]"}]], "Input",

 CellChangeTimes->{{3.4241986645559425`*^9, 3.4241986670695567`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], 

    RowBox[{

     RowBox[{"deg", "[", "a", "]"}], "+", 

     RowBox[{"deg", "[", "b", "]"}]}]], " ", 

   RowBox[{"t", "[", 

    RowBox[{"a", ",", "b", ",", 

     RowBox[{"d", "[", "c", "]"}]}], "]"}]}], "+", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], 

    RowBox[{"deg", "[", "a", "]"}]], " ", 

   RowBox[{"t", "[", 

    RowBox[{"a", ",", 

     RowBox[{"d", "[", "b", "]"}], ",", "c"}], "]"}]}], "+", 

  RowBox[{"t", "[", 

   RowBox[{

    RowBox[{"d", "[", "a", "]"}], ",", "b", ",", "c"}], "]"}]}]], "Output",

 CellChangeTimes->{3.424198667420061*^9, 3.4241987978676353`*^9}]

}, Open  ]],

Cell[BoxData[

 RowBox[{

  RowBox[{"d", "[", "X_Plus", "]"}], ":=", 

  RowBox[{"d", "/@", "X"}]}]], "Input",

 CellChangeTimes->{{3.4241988074414015`*^9, 3.424198812658904*^9}}],

Cell[BoxData[{

 RowBox[{

  RowBox[{"d", "[", 

   RowBox[{

    SuperscriptBox[

     RowBox[{"(", 

      RowBox[{"-", "1"}], ")"}], "n_."], " ", "X_"}], "]"}], ":=", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], "n"], " ", 

   RowBox[{"d", "[", "X", "]"}]}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"d", "[", 

   RowBox[{"s_\[Sigma]", " ", "X_"}], "]"}], ":=", 

  RowBox[{"s", " ", 

   RowBox[{"d", "[", "X", "]"}]}]}]}], "Input",

 CellChangeTimes->{{3.424199070639862*^9, 3.424199126960848*^9}, {

  3.424201311281747*^9, 3.4242013146165423`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"d", "[", 

   RowBox[{"d", "[", "_", "]"}], "]"}], ":=", "0"}]], "Input",

 CellChangeTimes->{{3.4241988199694157`*^9, 3.4241988233342543`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"tm", "[", 

   RowBox[{"___", ",", "0", ",", "___"}], "]"}], ":=", "0"}]], "Input",

 CellChangeTimes->{{3.424198831185544*^9, 3.424198836773579*^9}}],

Cell[BoxData[{

 RowBox[{

  RowBox[{"tm", "[", 

   RowBox[{"X___", ",", "Y_Plus", ",", "Z__"}], "]"}], ":=", 

  RowBox[{

   RowBox[{

    RowBox[{"tm", "[", 

     RowBox[{"X", ",", "#", ",", "Z"}], "]"}], "&"}], "/@", 

   "Y"}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"tm", "[", 

   RowBox[{"X__", ",", "Y_Plus", ",", "Z___"}], "]"}], ":=", 

  RowBox[{

   RowBox[{

    RowBox[{"tm", "[", 

     RowBox[{"X", ",", "#", ",", "Z"}], "]"}], "&"}], "/@", 

   "Y"}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"tm", "[", 

   RowBox[{"X___", ",", 

    RowBox[{

     SuperscriptBox[

      RowBox[{"(", 

       RowBox[{"-", "1"}], ")"}], "n_."], "Y_"}], ",", "Z___"}], "]"}], ":=", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], "n"], 

   RowBox[{"tm", "[", 

    RowBox[{"X", ",", "Y", ",", "Z"}], "]"}]}]}]}], "Input",

 CellChangeTimes->{{3.4242002106090574`*^9, 3.4242002710960336`*^9}, {

  3.4242009022035217`*^9, 3.4242009251565266`*^9}}],

Cell[BoxData[{

 RowBox[{

  RowBox[{"t", "[", 

   RowBox[{"X___", ",", "Y_Plus", ",", "Z__"}], "]"}], ":=", 

  RowBox[{

   RowBox[{

    RowBox[{"t", "[", 

     RowBox[{"X", ",", "#", ",", "Z"}], "]"}], "&"}], "/@", 

   "Y"}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"t", "[", 

   RowBox[{"X__", ",", "Y_Plus", ",", "Z___"}], "]"}], ":=", 

  RowBox[{

   RowBox[{

    RowBox[{"t", "[", 

     RowBox[{"X", ",", "#", ",", "Z"}], "]"}], "&"}], "/@", 

   "Y"}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"t", "[", 

   RowBox[{"X___", ",", 

    RowBox[{

     SuperscriptBox[

      RowBox[{"(", 

       RowBox[{"-", "1"}], ")"}], "n_."], "Y_"}], ",", "Z___"}], "]"}], ":=", 

  RowBox[{

   SuperscriptBox[

    RowBox[{"(", 

     RowBox[{"-", "1"}], ")"}], "n"], 

   RowBox[{"t", "[", 

    RowBox[{"X", ",", "Y", ",", "Z"}], "]"}]}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"t", "[", 

   RowBox[{"X___", ",", 

    RowBox[{"s_\[Sigma]", " ", "Y_"}], ",", "Z___"}], "]"}], ":=", 

  RowBox[{"s", " ", 

   RowBox[{"t", "[", 

    RowBox[{"X", ",", "Y", ",", "Z"}], "]"}]}]}]}], "Input",

 CellChangeTimes->{{3.424201175196066*^9, 3.4242012093451695`*^9}, {

  3.424201303861077*^9, 3.424201306564965*^9}}],

Cell[BoxData[{

 RowBox[{

  RowBox[{"t", "[", 

   RowBox[{"X__", ",", 

    RowBox[{"t", "[", "Y__", "]"}], ",", "Z___"}], "]"}], ":=", 

  RowBox[{"t", "[", 

   RowBox[{"X", ",", "Y", ",", "Z"}], "]"}]}], "\[IndentingNewLine]", 

 RowBox[{

  RowBox[{"t", "[", 

   RowBox[{"X___", ",", 

    RowBox[{"t", "[", "Y__", "]"}], ",", "Z__"}], "]"}], ":=", 

  RowBox[{"t", "[", 

   RowBox[{"X", ",", "Y", ",", "Z"}], "]"}]}]}], "Input",

 CellChangeTimes->{{3.42420197886168*^9, 3.424201999771747*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"deg", "[", 

   RowBox[{"d", "[", "X_", "]"}], "]"}], ":=", 

  RowBox[{

   RowBox[{"deg", "[", "X", "]"}], "-", "1"}]}]], "Input",

 CellChangeTimes->{{3.4241991779742017`*^9, 3.4241991852346416`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"deg", "[", "a_m", "]"}], ":=", 

  RowBox[{

   RowBox[{"Length", "[", "a", "]"}], "-", "2", "+", 

   RowBox[{"Plus", "@@", 

    RowBox[{"(", 

     RowBox[{"deg", "/@", 

      RowBox[{"List", "@@", "a"}]}], ")"}]}]}]}]], "Input",

 CellChangeTimes->{{3.4242003314828653`*^9, 3.4242003740440655`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"deg", "[", "a_tm", "]"}], ":=", 

  RowBox[{

   RowBox[{"Length", "[", "a", "]"}], "-", "1", "+", 

   RowBox[{"Plus", "@@", 

    RowBox[{"(", 

     RowBox[{"deg", "/@", 

      RowBox[{"List", "@@", "a"}]}], ")"}]}]}]}]], "Input",

 CellChangeTimes->{{3.4242003910885744`*^9, 3.4242003932016125`*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"d", "[", 

   RowBox[{"a", ":", 

    RowBox[{"HoldPattern", "[", 

     RowBox[{"tm", "[", "___", "]"}], "]"}]}], "]"}], ":=", 

  RowBox[{"With", "[", 

   RowBox[{

    RowBox[{"{", 

     RowBox[{"k", "=", 

      RowBox[{"Length", "[", "a", "]"}]}], "}"}], ",", "\[IndentingNewLine]", 

    RowBox[{

     RowBox[{

      UnderoverscriptBox["\[Sum]", 

       RowBox[{"l", "=", "1"}], 

       RowBox[{"k", "-", "1"}]], 

      RowBox[{

       RowBox[{"\[Sigma]", "[", 

        RowBox[{"l", ",", 

         RowBox[{"List", "@@", "a"}]}], "]"}], " ", 

       RowBox[{"t", "[", 

        RowBox[{

         RowBox[{"Take", "[", 

          RowBox[{"a", ",", "l"}], "]"}], ",", 

         RowBox[{"Drop", "[", 

          RowBox[{"a", ",", "l"}], "]"}]}], "]"}]}]}], "+", 

     "\[IndentingNewLine]", 

     RowBox[{

      UnderoverscriptBox["\[Sum]", 

       RowBox[{"l", "=", "1"}], "k"], 

      RowBox[{

       UnderoverscriptBox["\[Sum]", 

        RowBox[{"lp", "=", "0"}], 

        RowBox[{"l", "-", "1"}]], 

       RowBox[{"With", "[", 

        RowBox[{

         RowBox[{"{", 

          RowBox[{

           RowBox[{"l1", "=", "lp"}], ",", 

           RowBox[{"l2", "=", 

            RowBox[{"k", "-", "l", "+", "1"}]}], ",", 

           RowBox[{"l3", "=", 

            RowBox[{"l", "-", "1", "-", "lp"}]}]}], "}"}], ",", 

         RowBox[{

          SuperscriptBox[

           RowBox[{"(", 

            RowBox[{"-", "1"}], ")"}], 

           RowBox[{

            RowBox[{"l2", 

             RowBox[{"(", 

              RowBox[{"k", "-", "1", "+", 

               RowBox[{

                UnderoverscriptBox["\[Sum]", 

                 RowBox[{"j", "=", "1"}], "lp"], 

                RowBox[{"deg", "[", 

                 RowBox[{

                 "a", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}], 

                 "]"}]}]}], ")"}]}], "+", 

            RowBox[{

             RowBox[{"(", 

              RowBox[{"k", "-", "l"}], ")"}], 

             RowBox[{"(", "lp", ")"}]}]}]], 

          RowBox[{"tm", "@@", 

           RowBox[{"(", 

            RowBox[{

             RowBox[{"Take", "[", 

              RowBox[{

               RowBox[{"List", "@@", "a"}], ",", "lp"}], "]"}], "~", "Join", 

             "~", 

             RowBox[{"{", 

              RowBox[{"m", "@@", 

               RowBox[{"Take", "[", 

                RowBox[{

                 RowBox[{"List", "@@", "a"}], ",", 

                 RowBox[{"{", 

                  RowBox[{

                   RowBox[{"lp", "+", "1"}], ",", 

                   RowBox[{"lp", "+", "k", "-", "l", "+", "1"}]}], "}"}]}], 

                "]"}]}], "}"}], "~", "Join", "~", 

             RowBox[{"Drop", "[", 

              RowBox[{

               RowBox[{"List", "@@", "a"}], ",", 

               RowBox[{"lp", "+", "k", "-", "l", "+", "1"}]}], "]"}]}], 

            ")"}]}]}]}], "]"}]}]}]}]}], "\[IndentingNewLine]", 

   "]"}]}]], "Input",

 CellChangeTimes->{{3.4241976743520994`*^9, 3.4241979358681407`*^9}, {

  3.4241980205799503`*^9, 3.4241980693100204`*^9}, {3.4241981831937776`*^9, 

  3.424198219375805*^9}, {3.4241984876916237`*^9, 3.424198488572891*^9}, {

  3.4241987823553295`*^9, 3.424198785980542*^9}, {3.4241989378989906`*^9, 

  3.424198966930736*^9}, {3.4241990027121873`*^9, 3.424199045193272*^9}, {

  3.424200498352813*^9, 3.4242005144659824`*^9}, {3.4242006132580385`*^9, 

  3.424200613448312*^9}, {3.42420072089281*^9, 3.4242007272219105`*^9}, {

  3.424200772507027*^9, 3.424200782811845*^9}, {3.4242012813687344`*^9, 

  3.424201292494733*^9}, {3.4242071329028273`*^9, 3.424207333311*^9}, {

  3.4242075895995245`*^9, 3.424207624329464*^9}}],

Cell[BoxData[

 SuperscriptBox[

  RowBox[{"(", 

   RowBox[{"-", "1"}], ")"}], 

  RowBox[{

   RowBox[{

    RowBox[{"(", 

     RowBox[{"k", "-", "l", "+", "1"}], ")"}], 

    RowBox[{"(", 

     RowBox[{"k", "-", "1", "+", 

      RowBox[{

       UnderoverscriptBox["\[Sum]", 

        RowBox[{"j", "=", "1"}], "lp"], 

       RowBox[{"deg", "[", 

        RowBox[{"a", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}], 

        "]"}]}]}], ")"}]}], "+", 

   RowBox[{

    RowBox[{"(", 

     RowBox[{"k", "-", "l"}], ")"}], 

    RowBox[{"(", "lp", ")"}]}]}]]], "Input"],

Cell[BoxData[

 SuperscriptBox[

  RowBox[{"(", 

   RowBox[{"-", "1"}], ")"}], 

  RowBox[{

   RowBox[{

    RowBox[{"(", 

     RowBox[{"l1", "+", "l2"}], ")"}], "l3"}], "+", 

   RowBox[{"l2", " ", 

    RowBox[{"(", 

     RowBox[{"Plus", "@@", 

      RowBox[{"(", 

       RowBox[{"deg", "/@", 

        RowBox[{"Take", "[", 

         RowBox[{

          RowBox[{"List", "@@", "a"}], ",", "l1"}], "]"}]}], ")"}]}], 

     ")"}]}]}]]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[

 RowBox[{"Collect", "[", 

  RowBox[{

   RowBox[{"d", "[", 

    RowBox[{"d", "[", 

     RowBox[{"tm", "[", 

      RowBox[{"a", ",", "b"}], "]"}], "]"}], "]"}], ",", 

   RowBox[{"(", 

    RowBox[{"_tm", "|", "_d", "|", "_t"}], ")"}], ",", "Together"}], 

  "]"}]], "Input",

 CellChangeTimes->{{3.424198799299694*^9, 3.424198802804734*^9}, {

  3.4241991375360546`*^9, 3.4241991665077133`*^9}, {3.424199203320648*^9, 

  3.4241992055338306`*^9}}],

Cell[BoxData["0"], "Output",

 CellChangeTimes->{{3.424198803315469*^9, 3.424198837724947*^9}, {

   3.4241990626684*^9, 3.424199205944421*^9}, 3.424200074182886*^9, 

   3.4242001944958878`*^9, 3.424200409054408*^9, {3.4242004997448144`*^9, 

   3.4242005193029375`*^9}, 3.424200728914344*^9, 3.424200784414149*^9, 

   3.424200932707384*^9, 3.4242013241802945`*^9, 3.4242014646823263`*^9, {

   3.4242073102378225`*^9, 3.424207336445507*^9}, {3.42420760386003*^9, 

   3.4242076281149073`*^9}}]

}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[

 RowBox[{"d", "[", 

  RowBox[{"tm", "[", 

   RowBox[{"a", ",", "b"}], "]"}], "]"}]], "Input",

 CellChangeTimes->{{3.424200732269168*^9, 3.424200734181918*^9}}],

Cell[BoxData[

 RowBox[{

  RowBox[{"m", "[",