C2tauSystem.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     92799,       2549]
NotebookOptionsPosition[     89415,       2490]
NotebookOutlinePosition[     89779,       2506]
CellTagsIndexPosition[     89736,       2503]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"coeffs", " ", "=", " ", 
   RowBox[{"Flatten", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "0", ",", "0"}], "]"}], "\[Rule]", "1"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"0", ",", "0", ",", "0"}], "]"}], "\[Rule]", "1"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"0", ",", "0", ",", "0"}], "]"}], "\[Rule]", "1"}]}], "}"}],
       ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "0", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"0", ",", "0", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"0", ",", "0", ",", "1"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          SuperscriptBox["\[Sigma]2", "2"]]}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "0", ",", "2"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"0", ",", "0", ",", "2"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"0", ",", "0", ",", "2"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{"4", " ", 
           SuperscriptBox["\[Sigma]2", "4"]}]]}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "1", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"0", ",", "1", ",", "0"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"]]}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"0", ",", "1", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "1", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"0", ",", "1", ",", "1"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{
           SuperscriptBox["\[Sigma]1", "2"], " ", 
           SuperscriptBox["\[Sigma]2", "2"]}]]}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"0", ",", "1", ",", "1"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{
           SuperscriptBox["\[Sigma]1", "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"]}]]}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "2", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"0", ",", "2", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"0", ",", "2", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"1", ",", "0", ",", "0"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"]]}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"1", ",", "0", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"1", ",", "0", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"1", ",", "0", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"1", ",", "0", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"1", ",", "0", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"1", ",", "1", ",", "0"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"]}]]}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"1", ",", "1", ",", "0"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"]}]]}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"1", ",", "1", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"2", ",", "0", ",", "0"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]]}], ",", 
        RowBox[{
         RowBox[{"t1", "[", 
          RowBox[{"2", ",", "0", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t2", "[", 
          RowBox[{"2", ",", "0", ",", "0"}], "]"}], "\[Rule]", "0"}]}], "}"}],
       ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "0", ",", "3"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "1", ",", "2"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "2", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"0", ",", "3", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"1", ",", "0", ",", "2"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"1", ",", "1", ",", "1"}], "]"}], "\[Rule]", 
         FractionBox["1", 
          RowBox[{
           SuperscriptBox["\[Sigma]1", "2"], " ", 
           SuperscriptBox["\[Sigma]2", "2"]}]]}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"1", ",", "2", ",", "0"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"2", ",", "0", ",", "1"}], "]"}], "\[Rule]", "0"}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"2", ",", "1", ",", "0"}], "]"}], "\[Rule]", 
         RowBox[{"-", 
          FractionBox[
           RowBox[{
            RowBox[{"-", "2"}], "+", 
            RowBox[{"2", " ", "\[Sigma]1"}], "-", 
            SuperscriptBox["\[Sigma]1", "2"], "+", 
            RowBox[{"2", " ", "\[Sigma]2"}], "-", 
            SuperscriptBox["\[Sigma]2", "2"]}], 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "2"}], "+", "\[Sigma]1", "+", "\[Sigma]2"}], 
              ")"}], "2"]}]]}]}], ",", 
        RowBox[{
         RowBox[{"t0", "[", 
          RowBox[{"3", ",", "0", ",", "0"}], "]"}], "\[Rule]", "0"}]}], 
       "}"}]}], "}"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.798267904677664*^9, 3.798267930583515*^9}, {
  3.798269159859371*^9, 3.798269160093379*^9}, {3.7982768826453943`*^9, 
  3.798276883235548*^9}, {3.7982770597772818`*^9, 3.798277060190648*^9}, {
  3.798277996353805*^9, 3.7982779968181477`*^9}, {3.79827830262808*^9, 
  3.798278302922859*^9}, {3.798337144812024*^9, 3.7983371464821053`*^9}, {
  3.798365545852171*^9, 3.7983655675830107`*^9}, {3.79836560976753*^9, 
  3.798365610183666*^9}, {3.798365690368574*^9, 3.798365690817889*^9}, {
  3.798365920648077*^9, 3.79836592201336*^9}, {3.798619647253272*^9, 
  3.798619692432066*^9}},
 CellLabel->"In[76]:=",ExpressionUUID->"45f0ad4a-754f-41f8-a483-c077c6705864"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", "=", "4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tau0", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"t0", "[", 
        RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}], " ", 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", 
   "coeffs"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tau1", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"t1", "[", 
        RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}], " ", 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", 
   "coeffs"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tau2", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"t2", "[", 
        RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}], " ", 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", 
   "coeffs"}]}]}], "Input",
 CellChangeTimes->{{3.7982677099707823`*^9, 3.798267739500039*^9}, {
   3.7982677798971043`*^9, 3.798267780007847*^9}, 3.7982678842385798`*^9, {
   3.798267932646373*^9, 3.7982679373269653`*^9}, 3.7982684190827427`*^9, {
   3.798269162551404*^9, 3.798269162707549*^9}, {3.798277876290277*^9, 
   3.798277876447708*^9}, {3.798278000608186*^9, 3.7982780197867737`*^9}, 
   3.798337153329358*^9, 3.7983654999648323`*^9, {3.798365550123809*^9, 
   3.798365572843627*^9}, {3.798365613852214*^9, 3.7983656141615868`*^9}, 
   3.7983657002262077`*^9, 3.798619670393963*^9, 3.7986197008262053`*^9},
 CellLabel->"In[77]:=",ExpressionUUID->"dc3effc5-390c-4b5c-9f6f-981959ddecaa"],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{"D_logt", "^", "2"}], " ", 
    RowBox[{"Log", " ", "\\", "tau"}]}], " ", "=", " ", 
   RowBox[{"product", " ", "of", " ", 
    RowBox[{"other", " ", "\\", "taus"}]}]}], "*)"}]], "Input",
 CellChangeTimes->{{3.7984382879635878`*^9, 
  3.798438308715067*^9}},ExpressionUUID->"c3677654-05fe-4766-9219-\
6586698020e9"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"taude0", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Sum", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], "^", 
             RowBox[{"(", 
              RowBox[{"k", "+", "1"}], ")"}]}], 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Tau0", "[", "t", "]"}], "-", "1"}], ")"}], "^", 
              "k"}], "/", "k"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
       RowBox[{
        RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
        RowBox[{
         RowBox[{"p", "^", "2"}], " ", 
         RowBox[{"t", "^", "p"}]}]}]}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        FractionBox[
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"-", "a1"}]], " ", 
          SuperscriptBox["t", 
           RowBox[{"1", "-", "\[Sigma]1"}]], " ", "\[Epsilon]"}], "\[Xi]1"], 
        RowBox[{
         RowBox[{"Tau1", "[", "t", "]"}], "/", 
         RowBox[{"Tau0", "[", "t", "]"}]}]}], "//", "Simplify"}], ")"}]}], 
    "\[Equal]", "0"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"taude1", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Sum", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], "^", 
             RowBox[{"(", 
              RowBox[{"k", "+", "1"}], ")"}]}], 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Tau1", "[", "t", "]"}], "-", "1"}], ")"}], "^", 
              "k"}], "/", "k"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
       RowBox[{
        RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
        RowBox[{
         RowBox[{"p", "^", "2"}], " ", 
         RowBox[{"t", "^", "p"}]}]}]}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{
       FractionBox[
        RowBox[{
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"a1", "-", "a2"}]], " ", 
         SuperscriptBox["t", 
          RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}]], " ", "\[Xi]1", " ", 
         "\[Epsilon]"}], "\[Xi]2"], 
       RowBox[{"Tau0", "[", "t", "]"}], 
       RowBox[{
        RowBox[{"Tau2", "[", "t", "]"}], "/", 
        RowBox[{
         RowBox[{"Tau1", "[", "t", "]"}], "^", "2"}]}]}], ")"}]}], "\[Equal]",
     "0"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"taude2", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Sum", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], "^", 
             RowBox[{"(", 
              RowBox[{"k", "+", "1"}], ")"}]}], 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Tau2", "[", "t", "]"}], "-", "1"}], ")"}], "^", 
              "k"}], "/", "k"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
       RowBox[{
        RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
        RowBox[{
         RowBox[{"p", "^", "2"}], " ", 
         RowBox[{"t", "^", "p"}]}]}]}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\[ExponentialE]", "a2"], " ", 
       SuperscriptBox["t", "\[Sigma]2"], " ", "\[Xi]2", " ", "\[Epsilon]", 
       " ", 
       RowBox[{
        RowBox[{"Tau1", "[", "t", "]"}], "/", 
        RowBox[{"Tau2", "[", "t", "]"}]}]}], ")"}]}], "\[Equal]", "0"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.798268465298442*^9, 3.798268481061843*^9}, {
  3.798268819234508*^9, 3.798268834118739*^9}, {3.798268998744268*^9, 
  3.7982690675998917`*^9}, {3.798364739560524*^9, 3.7983647471913757`*^9}},
 CellLabel->"In[81]:=",ExpressionUUID->"89f1c56e-4e87-497f-94bb-3925c7acf9c2"],

Cell[BoxData[
 RowBox[{
  RowBox[{"coeffeq", " ", "=", " ", 
   RowBox[{"LogicalExpand", "[", 
    RowBox[{
     RowBox[{"Series", "[", 
      RowBox[{"taude0", ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]1", ",", "0", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]2", ",", "0", ",", "n"}], "}"}]}], "]"}], " ", "&&", 
     " ", 
     RowBox[{"Series", "[", 
      RowBox[{"taude1", ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]1", ",", "0", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]2", ",", "0", ",", "n"}], "}"}]}], "]"}], "&&", " ", 
     RowBox[{"Series", "[", 
      RowBox[{"taude2", ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]1", ",", "0", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]2", ",", "0", ",", "n"}], "}"}]}], "]"}]}], " ", 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.798267886457643*^9, 3.798267894896447*^9}, {
   3.7982679905057364`*^9, 3.798267991711586*^9}, {3.798268391067833*^9, 
   3.7982684375380583`*^9}, {3.798268838293969*^9, 3.798268874829187*^9}, {
   3.7982690953300047`*^9, 3.7982691181354847`*^9}, 3.798364751676113*^9},
 CellLabel->"In[84]:=",ExpressionUUID->"c4dc56d4-0065-4260-bff8-81b8f2387be3"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{"coeffeq", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"t", "\[Rule]", "1"}], ",", 
      RowBox[{"a1", "\[Rule]", "0"}], ",", 
      RowBox[{"a2", "\[Rule]", "0"}]}], "}"}]}], ",", " ", 
   RowBox[{"Flatten", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"t0", "[", 
         RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}], ",", 
        RowBox[{"t1", "[", 
         RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}], ",", 
        RowBox[{"t2", "[", 
         RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.798269126436328*^9, 3.798269139706469*^9}, {
  3.79827798783261*^9, 3.79827799096078*^9}, {3.798278063754718*^9, 
  3.7982780683996077`*^9}},
 CellLabel->"In[85]:=",ExpressionUUID->"abf2448e-47fa-4ddf-9770-0706af56979c"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", "}"}], "}"}]], "Output",
 CellChangeTimes->{3.798619634137121*^9, 3.798619679637342*^9, 
  3.798619713174513*^9},
 CellLabel->"Out[85]=",ExpressionUUID->"78d3033e-65cb-42ac-bd88-f589928a8f24"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"t0", "[", 
       RowBox[{"k", ",", "k", ",", 
        RowBox[{"k", "+", "1"}]}], "]"}], "/.", "coeffs"}], "//", 
     "Simplify"}], ")"}], ",", 
   RowBox[{"{", 
    RowBox[{"k", ",", "0", ",", "2"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.798269886983911*^9, 3.798269928071876*^9}, {
   3.7982699774178457`*^9, 3.798270004778483*^9}, {3.7982766791960278`*^9, 
   3.798276680754891*^9}, {3.7982783238036737`*^9, 3.798278332686254*^9}, {
   3.7983659488714533`*^9, 3.798365974800027*^9}, {3.798366038897887*^9, 
   3.798366087363549*^9}, {3.798366183179206*^9, 3.7983662730815496`*^9}, {
   3.798366334436494*^9, 3.798366348092515*^9}, {3.798436159707526*^9, 
   3.798436172369213*^9}, 3.798619467115788*^9, {3.798619521818869*^9, 
   3.798619569746402*^9}},
 CellLabel->"In[56]:=",ExpressionUUID->"d333b561-ecda-412a-98dc-2a212e8a39bb"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0", ",", 
   FractionBox["1", 
    RowBox[{"4", " ", 
     SuperscriptBox["\[Sigma]2", "4"], " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "+", "\[Sigma]2"}], ")"}], "2"], " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"\[Sigma]1", "+", "\[Sigma]2"}], ")"}], "2"]}]], ",", 
   RowBox[{"t0", "[", 
    RowBox[{"2", ",", "2", ",", "3"}], "]"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.798365956516246*^9, 3.79836597518993*^9}, {
   3.798366039726715*^9, 3.7983660878208933`*^9}, {3.798366185084725*^9, 
   3.79836627370603*^9}, {3.798366338932425*^9, 3.7983663489681807`*^9}, 
   3.798436173430624*^9, {3.798619464442711*^9, 3.79861946854456*^9}, {
   3.798619536287662*^9, 3.798619570114386*^9}},
 CellLabel->"Out[56]=",ExpressionUUID->"b490ca9c-8c71-4499-8d9b-743b884e9457"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(", 
  RowBox[{
   RowBox[{
    RowBox[{"t0", "[", 
     RowBox[{"2", ",", "2", ",", "3"}], "]"}], "/.", "coeffs"}], "//", 
   "Simplify"}], ")"}]], "Input",
 CellChangeTimes->{{3.798619478918318*^9, 3.798619487521227*^9}},
 CellLabel->"In[53]:=",ExpressionUUID->"84a9803d-7eb4-4040-a97e-d089f6f4b672"],

Cell[BoxData[
 RowBox[{"t0", "[", 
  RowBox[{"2", ",", "2", ",", "3"}], "]"}]], "Output",
 CellChangeTimes->{{3.798619482200333*^9, 3.798619487936987*^9}},
 CellLabel->"Out[53]=",ExpressionUUID->"281a22b1-6efe-46de-9acd-50699c53d2e7"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.798436174866186*^9, 
  3.798436175520501*^9}},ExpressionUUID->"14e70c01-de81-4b68-8c97-\
77c71c3eec30"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"t0", "[", 
    RowBox[{"2", ",", "2", ",", "2"}], "]"}], "/.", "coeffs"}], "//", 
  "Simplify"}]], "Input",
 CellChangeTimes->{{3.798366277892082*^9, 3.798366279217622*^9}},
 CellLabel->
  "In[628]:=",ExpressionUUID->"e659888a-46d1-4dcb-b09a-99b9bc05161b"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"3", " ", 
    SuperscriptBox["\[Sigma]2", "4"], " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       SuperscriptBox["\[Sigma]2", "2"]}], ")"}], "2"]}], "-", 
   RowBox[{
    SuperscriptBox["\[Sigma]1", "2"], " ", 
    SuperscriptBox["\[Sigma]2", "2"], " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       SuperscriptBox["\[Sigma]2", "2"]}], ")"}], "2"], " ", 
    RowBox[{"(", 
     RowBox[{"4", "+", 
      RowBox[{"3", " ", 
       SuperscriptBox["\[Sigma]2", "2"]}]}], ")"}]}], "+", 
   RowBox[{
    SuperscriptBox["\[Sigma]1", "8"], " ", 
    RowBox[{"(", 
     RowBox[{"3", "-", 
      RowBox[{"3", " ", 
       SuperscriptBox["\[Sigma]2", "2"]}], "+", 
      RowBox[{"4", " ", 
       SuperscriptBox["\[Sigma]2", "4"]}]}], ")"}]}], "+", 
   RowBox[{
    SuperscriptBox["\[Sigma]1", "6"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "6"}], "+", 
      RowBox[{"2", " ", 
       SuperscriptBox["\[Sigma]2", "2"]}], "+", 
      RowBox[{"4", " ", 
       SuperscriptBox["\[Sigma]2", "4"]}], "-", 
      RowBox[{"8", " ", 
       SuperscriptBox["\[Sigma]2", "6"]}]}], ")"}]}], "+", 
   RowBox[{
    SuperscriptBox["\[Sigma]1", "4"], " ", 
    RowBox[{"(", 
     RowBox[{"3", "+", 
      RowBox[{"5", " ", 
       SuperscriptBox["\[Sigma]2", "2"]}], "-", 
      RowBox[{"12", " ", 
       SuperscriptBox["\[Sigma]2", "4"]}], "+", 
      RowBox[{"4", " ", 
       SuperscriptBox["\[Sigma]2", "6"]}], "+", 
      RowBox[{"4", " ", 
       SuperscriptBox["\[Sigma]2", "8"]}]}], ")"}]}]}], 
  RowBox[{"8", " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"], " ", 
   SuperscriptBox["\[Sigma]1", "4"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"1", "+", "\[Sigma]1"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"], " ", 
   SuperscriptBox["\[Sigma]2", "4"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"1", "+", "\[Sigma]2"}], ")"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"\[Sigma]1", "+", "\[Sigma]2"}], ")"}], "2"]}]]], "Output",
 CellChangeTimes->{3.7983662796883717`*^9},
 CellLabel->
  "Out[628]=",ExpressionUUID->"6e780e73-6fdf-446d-be3b-c2b7f24c3c2a"]
}, Open  ]],

Cell[BoxData["t1"], "Input",
 CellChangeTimes->{{3.7983662333969316`*^9, 
  3.7983662341841087`*^9}},ExpressionUUID->"e8a00397-aef7-4303-94cc-\
3e1f0f0ca554"],

Cell[BoxData[
 RowBox[{
  RowBox[{"TauTest", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"t0", "[", 
        RowBox[{"n0", ",", "n1", ",", "n2"}], "]"}], " ", 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", 
        RowBox[{"-", "n"}], ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", 
        RowBox[{"-", "n"}], ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", 
        RowBox[{"-", "n"}], ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", "coeffs"}]}]], "Input",\

 CellChangeTimes->{{3.798280137308864*^9, 
  3.798280145442711*^9}},ExpressionUUID->"e3bc3cd1-d9c0-4245-bcce-\
faef4b39dce9"],

Cell[BoxData[
 RowBox[{"Table", "[", "]"}]], "Input",
 CellChangeTimes->{{3.798280162193467*^9, 
  3.798280163756844*^9}},ExpressionUUID->"540b1cc2-38c0-48db-9004-\
b39b156bb7d7"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "+", "1"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", "z", "]"}], 
   RowBox[{"G", "[", "z", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "+", "2"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "+", "1"}], "]"}], 
   RowBox[{"G", "[", 
    RowBox[{"z", "+", "1"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "+", "3"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "+", "2"}], "]"}], 
   RowBox[{"G", "[", 
    RowBox[{"z", "+", "2"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "+", "4"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "+", "3"}], "]"}], 
   RowBox[{"G", "[", 
    RowBox[{"z", "+", "3"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "-", "1"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"G", "[", "z", "]"}], "/", 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "1"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "-", "2"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"G", "[", 
    RowBox[{"z", "-", "1"}], "]"}], "/", 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "2"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "-", "3"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"G", "[", 
    RowBox[{"z", "-", "2"}], "]"}], "/", 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "3"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"G", "[", 
   RowBox[{"z_", "-", "4"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"G", "[", 
    RowBox[{"z", "-", "3"}], "]"}], "/", 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "4"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "+", "1"}], "]"}], ":=", 
  RowBox[{"z", " ", 
   RowBox[{"\[CapitalGamma]", "[", "z", "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "+", "2"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"z", "+", "1"}], ")"}], 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "+", "1"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "+", "3"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"z", "+", "2"}], ")"}], 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "+", "2"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "+", "4"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"z", "+", "3"}], ")"}], 
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "+", "3"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "-", "1"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", "z", "]"}], "/", 
   RowBox[{"(", 
    RowBox[{"z", "-", "1"}], ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "-", "2"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "1"}], "]"}], "/", 
   RowBox[{"(", 
    RowBox[{"z", "-", "2"}], ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "-", "3"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "2"}], "]"}], "/", 
   RowBox[{"(", 
    RowBox[{"z", "-", "3"}], ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]", "[", 
   RowBox[{"z_", "-", "4"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[CapitalGamma]", "[", 
    RowBox[{"z", "-", "3"}], "]"}], "/", 
   RowBox[{"(", 
    RowBox[{"z", "-", "4"}], ")"}]}]}]}], "Input",
 CellLabel->"In[90]:=",ExpressionUUID->"69c1f79b-4733-4927-a510-0c734b2ee0c5"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"F0", "[", 
       RowBox[{
        RowBox[{"\[Sigma]1", "+", "1"}], ",", "\[Sigma]2"}], "]"}], 
      RowBox[{
       RowBox[{"F0", "[", 
        RowBox[{
         RowBox[{"\[Sigma]1", "-", "1"}], ",", "\[Sigma]2"}], "]"}], "/", " ", 
       RowBox[{
        RowBox[{"F0", "[", 
         RowBox[{"\[Sigma]1", ",", "\[Sigma]2"}], "]"}], "^", "2"}]}]}], "/.", 
     RowBox[{
      RowBox[{"F0", "[", 
       RowBox[{"\[Sigma]1_", ",", "\[Sigma]2_"}], "]"}], "\[Rule]", 
      RowBox[{"1", "/", 
       RowBox[{"Abs", "[", 
        RowBox[{
         RowBox[{"G", "[", 
          RowBox[{"1", "-", "\[Sigma]1"}], "]"}], 
         RowBox[{"G", "[", 
          RowBox[{"1", "+", "\[Sigma]1"}], "]"}], 
         RowBox[{"G", "[", 
          RowBox[{"1", "-", "\[Sigma]2"}], "]"}], 
         RowBox[{"G", "[", 
          RowBox[{"1", "+", "\[Sigma]2"}], "]"}], 
         RowBox[{"G", "[", 
          RowBox[{"1", "+", "\[Sigma]1", "-", "\[Sigma]2"}], "]"}], 
         RowBox[{"G", "[", 
          RowBox[{"1", "+", "\[Sigma]1", "-", "\[Sigma]2"}], "]"}]}], 
        "]"}]}]}]}], "//.", 
    RowBox[{
     RowBox[{"Abs", "[", 
      RowBox[{"x_", " ", "y_"}], " ", "]"}], "\[Rule]", " ", 
     RowBox[{
      RowBox[{"Abs", "[", "x", "]"}], 
      RowBox[{"Abs", "[", "y", "]"}]}]}]}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Sigma]1", " ", ">", "0"}], ",", " ", 
     RowBox[{"\[Sigma]2", " ", ">", "0"}], ",", " ", 
     RowBox[{"\[Sigma]3", " ", ">", "0"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.798270076494328*^9, 3.798270144173478*^9}, {
  3.7982701901643457`*^9, 3.798270191691409*^9}, {3.798270456279217*^9, 
  3.798270617403667*^9}},
 CellLabel->
  "In[145]:=",ExpressionUUID->"5e004098-772f-467a-92d1-a6d9aecacf72"],

Cell[BoxData[
 FractionBox["1", 
  RowBox[{
   SuperscriptBox["\[Sigma]1", "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"]}]]], "Output",
 CellChangeTimes->{
  3.798270097035798*^9, 3.798270148898827*^9, 3.7982701923889713`*^9, {
   3.798270452324078*^9, 3.798270617963896*^9}},
 CellLabel->
  "Out[145]=",ExpressionUUID->"86cd1d40-c155-4f35-9a8d-0a9483261ab1"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", "=", "2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vTau0", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"A0", "[", 
         RowBox[{
          RowBox[{"\[Sigma]1", "+", "n1", "-", "n0", "+", "1"}], ",", 
          RowBox[{"\[Sigma]2", "+", "n2", "-", "n1"}]}], "]"}], " ", "/", 
        RowBox[{"(", 
         RowBox[{"A0", "[", 
          RowBox[{
           RowBox[{"\[Sigma]1", "+", "1"}], ",", "\[Sigma]2"}], "]"}], 
         ")"}]}], 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", 
   "coeffs"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vTau1", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"A1", "[", 
         RowBox[{
          RowBox[{"\[Sigma]1", "+", "n1", "-", "n0", "+", "1"}], ",", 
          RowBox[{"\[Sigma]2", "+", "n2", "-", "n1", "-", "1"}]}], "]"}], " ",
         "/", 
        RowBox[{"(", 
         RowBox[{"A1", "[", 
          RowBox[{
           RowBox[{"\[Sigma]1", "+", "1"}], ",", 
           RowBox[{"\[Sigma]2", "-", "1"}]}], "]"}], ")"}]}], " ", 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", 
   "coeffs"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vTau2", "[", "t_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Sum", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"A2", "[", 
         RowBox[{
          RowBox[{"\[Sigma]1", "+", "n1", "-", "n0"}], ",", 
          RowBox[{"\[Sigma]2", "+", "n2", "-", "n1", "+", "1"}]}], "]"}], " ",
         "/", 
        RowBox[{"(", 
         RowBox[{"A2", "[", 
          RowBox[{"\[Sigma]1", ",", 
           RowBox[{"\[Sigma]2", "+", "1"}]}], "]"}], ")"}]}], 
       RowBox[{"t", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"1", "-", "\[Sigma]1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "\[Sigma]2"}]}], ")"}]}], 
       RowBox[{"E", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"n0", 
           RowBox[{"(", 
            RowBox[{"-", "a1"}], ")"}]}], "+", 
          RowBox[{"n1", 
           RowBox[{"(", 
            RowBox[{"a1", "-", "a2"}], ")"}]}], "+", 
          RowBox[{"n2", " ", "a2"}]}], ")"}]}], 
       RowBox[{"\[Epsilon]", "^", 
        RowBox[{"(", 
         RowBox[{"n0", "+", "n1", "+", "n2"}], ")"}]}], 
       RowBox[{"\[Xi]1", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n0"}], "+", "n1"}], ")"}]}], 
       RowBox[{"\[Xi]2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "n1"}], "+", "n2"}], ")"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"n0", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n1", ",", "0", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"n2", ",", "0", ",", "n"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{"n", "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], "/.", 
   "coeffs"}]}]}], "Input",
 CellChangeTimes->{{3.7982702376380043`*^9, 3.79827030829081*^9}, {
  3.798270371619322*^9, 3.798270391530945*^9}, {3.798270660314815*^9, 
  3.7982706795230103`*^9}, {3.798271264245088*^9, 3.798271346928906*^9}, {
  3.798280480104307*^9, 3.798280503017746*^9}, {3.7982806423717413`*^9, 
  3.798280652442116*^9}, {3.79828079996706*^9, 3.7982808411858807`*^9}},
 CellLabel->
  "In[446]:=",ExpressionUUID->"27366b93-7abd-4206-bdb0-2a3a75d81e3a"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"vtaude0", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"vTau0", "[", "t", "]"}], "//", "Normal"}], "]"}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
       RowBox[{
        RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
        RowBox[{
         RowBox[{"p", "^", "2"}], " ", 
         RowBox[{"t", "^", "p"}]}]}]}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        FractionBox[
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"-", "a1"}]], " ", 
          SuperscriptBox["t", 
           RowBox[{"1", "-", "\[Sigma]1"}]], " ", "\[Epsilon]"}], "\[Xi]1"], 
        RowBox[{
         RowBox[{"vTau1", "[", "t", "]"}], "/", 
         RowBox[{"vTau0", "[", "t", "]"}]}]}], "//", "Simplify"}], ")"}]}], 
    "\[Equal]", "0"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vtaude1", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"vTau1", "[", "t", "]"}], "//", "Normal"}], "]"}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
       RowBox[{
        RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
        RowBox[{
         RowBox[{"p", "^", "2"}], " ", 
         RowBox[{"t", "^", "p"}]}]}]}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{
       FractionBox[
        RowBox[{
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"a1", "-", "a2"}]], " ", 
         SuperscriptBox["t", 
          RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}]], " ", "\[Xi]1", " ", 
         "\[Epsilon]"}], "\[Xi]2"], 
       RowBox[{"vTau0", "[", "t", "]"}], 
       RowBox[{
        RowBox[{"vTau2", "[", "t", "]"}], "/", 
        RowBox[{
         RowBox[{"vTau1", "[", "t", "]"}], "^", "2"}]}]}], ")"}]}], 
    "\[Equal]", "0"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vtaude2", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"vTau2", "[", "t", "]"}], "//", "Normal"}], "]"}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
         RowBox[{
          RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
       RowBox[{
        RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
        RowBox[{
         RowBox[{"p", "^", "2"}], " ", 
         RowBox[{"t", "^", "p"}]}]}]}], ")"}], "-", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["\[ExponentialE]", "a2"], " ", 
       SuperscriptBox["t", "\[Sigma]2"], " ", "\[Xi]2", " ", "\[Epsilon]", 
       " ", 
       RowBox[{
        RowBox[{"vTau1", "[", "t", "]"}], "/", 
        RowBox[{"vTau2", "[", "t", "]"}]}]}], ")"}]}], "\[Equal]", "0"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.798270312251871*^9, 3.798270335347782*^9}, {
  3.79828077607124*^9, 3.79828077737891*^9}},
 CellLabel->
  "In[450]:=",ExpressionUUID->"3a0fe918-4fc7-4e0c-9e51-0767004f7db9"],

Cell[BoxData[
 RowBox[{
  RowBox[{"coeffeq", " ", "=", " ", 
   RowBox[{"LogicalExpand", "[", 
    RowBox[{
     RowBox[{"Series", "[", 
      RowBox[{"vtaude0", ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]1", ",", "0", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]2", ",", "0", ",", "n"}], "}"}]}], "]"}], " ", "&&", 
     " ", 
     RowBox[{"Series", "[", 
      RowBox[{"vtaude1", ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]1", ",", "0", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]2", ",", "0", ",", "n"}], "}"}]}], "]"}], "&&", " ", 
     RowBox[{"Series", "[", 
      RowBox[{"vtaude2", ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]1", ",", "0", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Xi]2", ",", "0", ",", "n"}], "}"}]}], "]"}]}], " ", 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7982703381434298`*^9, 3.798270341581685*^9}, {
  3.798280525622806*^9, 3.798280545529504*^9}},
 CellLabel->
  "In[453]:=",ExpressionUUID->"569e09a5-261a-4388-88a9-94f48502493f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"Log", "[", 
      RowBox[{
       RowBox[{"Tau0", "[", "t", "]"}], "//", "Normal"}], "]"}], "+", 
     RowBox[{
      RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
      RowBox[{"(", 
       RowBox[{"n", "+", "1"}], ")"}]}], "+", 
     RowBox[{
      RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "n"}], "+", "1"}], ")"}]}], "+", 
     RowBox[{
      RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "n"}], "+", "1"}], ")"}]}]}], ")"}], "/.", 
   RowBox[{
    RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
    RowBox[{
     RowBox[{"p", "^", "2"}], " ", 
     RowBox[{"t", "^", "p"}]}]}]}], ")"}]], "Input",
 CellLabel->
  "In[260]:=",ExpressionUUID->"77ae1bbe-bae5-469a-95ef-b7791b3752f1"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"(", 
    InterpretationBox[
     RowBox[{
      InterpretationBox[
       SuperscriptBox[
        RowBox[{"O", "[", "\[Xi]2", "]"}], "7"],
       SeriesData[$CellContext`\[Xi]2, 0, {}, 7, 7, 1],
       Editable->False], "+", 
      InterpretationBox[
       SuperscriptBox[
        RowBox[{"O", "[", "\[Xi]1", "]"}], "7"],
       SeriesData[$CellContext`\[Xi]1, 0, {}, 0, 7, 1],
       Editable->False]}],
     SeriesData[$CellContext`\[Xi]1, 0, {
       SeriesData[$CellContext`\[Xi]2, 0, {}, 7, 7, 1]}, 0, 7, 1],
     Editable->False], ")"}], "+", 
   FractionBox[
    RowBox[{
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{"-", "a1"}]], " ", 
     SuperscriptBox["t", 
      RowBox[{"1", "-", "\[Sigma]1"}]], " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "-", "\[Sigma]1"}], ")"}], "2"], " ", "\[Epsilon]"}], 
    RowBox[{"\[Xi]1", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"]}]], "+", 
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "a1"}]], " ", 
         SuperscriptBox["t", 
          RowBox[{"2", "-", 
           RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
         SuperscriptBox[
          RowBox[{"(", 
           RowBox[{"2", "-", 
            RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"]}], 
        RowBox[{
         SuperscriptBox["\[Xi]1", "2"], " ", 
         SuperscriptBox[
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]]}], "+", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"-", "2"}], " ", "a1"}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"2", "-", 
             RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"], " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "-", 
             RowBox[{"2", " ", "\[Sigma]1"}], "+", 
             SuperscriptBox["\[Sigma]1", "2"]}], ")"}]}], 
          RowBox[{"8", " ", 
           SuperscriptBox["\[Xi]1", "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{"-", "a2"}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"1", "-", "\[Sigma]2"}]], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"1", "-", "\[Sigma]2"}], ")"}], "2"]}], 
          RowBox[{"\[Xi]2", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"]}]]}], 
        ")"}]}]}], ")"}], " ", 
    SuperscriptBox["\[Epsilon]", "2"]}], "+", 
   RowBox[{
    FractionBox["1", "3"], " ", 
    RowBox[{"(", 
     RowBox[{
      FractionBox[
       RowBox[{
        SuperscriptBox["\[ExponentialE]", 
         RowBox[{"-", "a1"}]], " ", 
        SuperscriptBox["t", 
         RowBox[{"1", "-", "\[Sigma]1"}]], " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Sigma]1"}], ")"}], "2"], " ", 
        RowBox[{"(", 
         RowBox[{
          FractionBox[
           RowBox[{"9", " ", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{
              RowBox[{"-", "2"}], " ", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"]}], 
           RowBox[{"8", " ", 
            SuperscriptBox["\[Xi]1", "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]], "+", 
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{
              RowBox[{"-", "2"}], " ", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"], " ", 
            "\[Sigma]1"}], 
           RowBox[{"4", " ", 
            SuperscriptBox["\[Xi]1", "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]], "-", 
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{
              RowBox[{"-", "2"}], " ", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"], " ", 
            SuperscriptBox["\[Sigma]1", "2"]}], 
           RowBox[{"8", " ", 
            SuperscriptBox["\[Xi]1", "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]], "-", 
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"-", "a2"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"1", "-", "\[Sigma]2"}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"1", "-", "\[Sigma]2"}], ")"}], "2"]}], 
           RowBox[{"\[Xi]2", " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"]}]]}], 
         ")"}]}], 
       RowBox[{"\[Xi]1", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"]}]], "-", 
      FractionBox[
       RowBox[{"2", " ", 
        SuperscriptBox["\[ExponentialE]", 
         RowBox[{"-", "a1"}]], " ", 
        SuperscriptBox["t", 
         RowBox[{"1", "-", "\[Sigma]1"}]], " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Sigma]1"}], ")"}], "2"], " ", 
        RowBox[{"(", 
         RowBox[{
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{
              RowBox[{"-", "2"}], " ", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"], " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "-", 
              RowBox[{"2", " ", "\[Sigma]1"}], "+", 
              SuperscriptBox["\[Sigma]1", "2"]}], ")"}]}], 
           RowBox[{"8", " ", 
            SuperscriptBox["\[Xi]1", "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]], "+", 
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"-", "a2"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"1", "-", "\[Sigma]2"}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"1", "-", "\[Sigma]2"}], ")"}], "2"]}], 
           RowBox[{"\[Xi]2", " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"]}]]}], 
         ")"}]}], 
       RowBox[{"\[Xi]1", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"]}]], "+", 
      RowBox[{"3", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"9", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"3", " ", "a2"}]], " ", 
          SuperscriptBox["t", 
           RowBox[{"3", " ", "\[Sigma]2"}]], " ", 
          SuperscriptBox["\[Xi]2", "3"], " ", 
          SuperscriptBox["\[Sigma]2", "2"], " ", 
          RowBox[{"t0", "[", 
           RowBox[{"0", ",", "0", ",", "3"}], "]"}]}], "+", 
         RowBox[{
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"a1", "+", "a2"}]], " ", 
          SuperscriptBox["t", 
           RowBox[{"\[Sigma]1", "+", "\[Sigma]2"}]], " ", "\[Xi]1", " ", 
          "\[Xi]2", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"\[Sigma]1", "+", "\[Sigma]2"}], ")"}], "2"], " ", 
          RowBox[{"t0", "[", 
           RowBox[{"0", ",", "1", ",", "2"}], "]"}]}], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"2", " ", "a1"}], "-", "a2"}]], " ", 
           SuperscriptBox["t", 
            RowBox[{
             RowBox[{"2", " ", "\[Sigma]1"}], "-", "\[Sigma]2"}]], " ", 
           SuperscriptBox["\[Xi]1", "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"2", " ", "\[Sigma]1"}], "-", "\[Sigma]2"}], ")"}], 
            "2"], " ", 
           RowBox[{"t0", "[", 
            RowBox[{"0", ",", "2", ",", "1"}], "]"}]}], "\[Xi]2"], "+", 
         FractionBox[
          RowBox[{"9", " ", 
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{"3", " ", 
             RowBox[{"(", 
              RowBox[{"a1", "-", "a2"}], ")"}]}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"3", " ", 
             RowBox[{"(", 
              RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}]}]], " ", 
           SuperscriptBox["\[Xi]1", "3"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
           RowBox[{"t0", "[", 
            RowBox[{"0", ",", "3", ",", "0"}], "]"}]}], 
          SuperscriptBox["\[Xi]2", "3"]], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"-", "a1"}], "+", 
             RowBox[{"2", " ", "a2"}]}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"1", "-", "\[Sigma]1", "+", 
             RowBox[{"2", " ", "\[Sigma]2"}]}]], " ", 
           SuperscriptBox["\[Xi]2", "2"], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"1", "-", "\[Sigma]1", "+", 
              RowBox[{"2", " ", "\[Sigma]2"}]}], ")"}], "2"], " ", 
           RowBox[{"t0", "[", 
            RowBox[{"1", ",", "0", ",", "2"}], "]"}]}], "\[Xi]1"], "+", 
         RowBox[{"t", " ", 
          RowBox[{"t0", "[", 
           RowBox[{"1", ",", "1", ",", "1"}], "]"}]}], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{"a1", "-", 
             RowBox[{"2", " ", "a2"}]}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"1", "+", "\[Sigma]1", "-", 
             RowBox[{"2", " ", "\[Sigma]2"}]}]], " ", "\[Xi]1", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"1", "+", "\[Sigma]1", "-", 
              RowBox[{"2", " ", "\[Sigma]2"}]}], ")"}], "2"], " ", 
           RowBox[{"t0", "[", 
            RowBox[{"1", ",", "2", ",", "0"}], "]"}]}], 
          SuperscriptBox["\[Xi]2", "2"]], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "a1"}], "+", "a2"}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"2", "-", 
             RowBox[{"2", " ", "\[Sigma]1"}], "+", "\[Sigma]2"}]], " ", 
           "\[Xi]2", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}], "+", "\[Sigma]2"}], ")"}], 
            "2"], " ", 
           RowBox[{"t0", "[", 
            RowBox[{"2", ",", "0", ",", "1"}], "]"}]}], 
          SuperscriptBox["\[Xi]1", "2"]], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"-", "a1"}], "-", "a2"}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"2", "-", "\[Sigma]1", "-", "\[Sigma]2"}]], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"2", "-", "\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], 
           " ", 
           RowBox[{"t0", "[", 
            RowBox[{"2", ",", "1", ",", "0"}], "]"}]}], 
          RowBox[{"\[Xi]1", " ", "\[Xi]2"}]], "+", 
         FractionBox[
          RowBox[{
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"-", "3"}], " ", "a1"}]], " ", 
           SuperscriptBox["t", 
            RowBox[{"3", "-", 
             RowBox[{"3", " ", "\[Sigma]1"}]}]], " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"3", "-", 
              RowBox[{"3", " ", "\[Sigma]1"}]}], ")"}], "2"], " ", 
           RowBox[{"t0", "[", 
            RowBox[{"3", ",", "0", ",", "0"}], "]"}]}], 
          SuperscriptBox["\[Xi]1", "3"]]}], ")"}]}]}], ")"}], " ", 
    SuperscriptBox["\[Epsilon]", "3"]}], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "4"],
    SeriesData[$CellContext`\[Epsilon], 0, {}, 0, 4, 1],
    Editable->False]}],
  SeriesData[$CellContext`\[Epsilon], 0, {
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {}, 7, 7, 1]}, 0, 7, 1], 
    E^(-$CellContext`a1) $CellContext`t^(
      1 - $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-1) (
      1 - $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-2), 
    Rational[1, 2] (-
      E^((-2) $CellContext`a1) $CellContext`t^(2 - 
        2 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-2) (2 - 
        2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4) + 
     2 (Rational[1, 8] 
        E^((-2) $CellContext`a1) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-2) (2 - 
          2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4) (-1 - 
         2 $CellContext`\[Sigma]1 + $CellContext`\[Sigma]1^2) + 
       E^(-$CellContext`a2) $CellContext`t^(
          1 - $CellContext`\[Sigma]2) $CellContext`\[Xi]2^(-1) (
          1 - $CellContext`\[Sigma]2)^2 ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2))), 
    Rational[1, 3] (
     E^(-$CellContext`a1) $CellContext`t^(
        1 - $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-1) (
        1 - $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-2) (
       Rational[9, 8] 
        E^((-2) $CellContext`a1) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-2) (2 - 
          2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4) + 
       Rational[1, 4] 
        E^((-2) $CellContext`a1) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-2) (2 - 
          2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4) \
$CellContext`\[Sigma]1 + 
       Rational[-1, 8] 
        E^((-2) $CellContext`a1) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-2) (2 - 
          2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4) \
$CellContext`\[Sigma]1^2 - 
       E^(-$CellContext`a2) $CellContext`t^(
         1 - $CellContext`\[Sigma]2) $CellContext`\[Xi]2^(-1) (
         1 - $CellContext`\[Sigma]2)^2 ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2)) - 2 
     E^(-$CellContext`a1) $CellContext`t^(
       1 - $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-1) (
       1 - $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-2) (
      Rational[1, 8] 
       E^((-2) $CellContext`a1) $CellContext`t^(2 - 
         2 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-2) (2 - 
         2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4) (-1 - 
        2 $CellContext`\[Sigma]1 + $CellContext`\[Sigma]1^2) + 
      E^(-$CellContext`a2) $CellContext`t^(
         1 - $CellContext`\[Sigma]2) $CellContext`\[Xi]2^(-1) (
         1 - $CellContext`\[Sigma]2)^2 ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2)) + 
     3 (9 E^(3 $CellContext`a2) $CellContext`t^(
          3 $CellContext`\[Sigma]2) $CellContext`\[Xi]2^3 $CellContext`\
\[Sigma]2^2 $CellContext`t0[0, 0, 3] + 
       E^($CellContext`a1 + $CellContext`a2) $CellContext`t^($CellContext`\
\[Sigma]1 + $CellContext`\[Sigma]2) $CellContext`\[Xi]1 $CellContext`\[Xi]2 \
($CellContext`\[Sigma]1 + $CellContext`\[Sigma]2)^2 $CellContext`t0[0, 1, 2] + 
       E^(2 $CellContext`a1 - $CellContext`a2) $CellContext`t^(
          2 $CellContext`\[Sigma]1 - $CellContext`\[Sigma]2) $CellContext`\
\[Xi]1^2 $CellContext`\[Xi]2^(-1) (
          2 $CellContext`\[Sigma]1 - $CellContext`\[Sigma]2)^2 \
$CellContext`t0[0, 2, 1] + 
       9 E^(3 ($CellContext`a1 - $CellContext`a2)) $CellContext`t^(
          3 ($CellContext`\[Sigma]1 - $CellContext`\[Sigma]2)) $CellContext`\
\[Xi]1^3 $CellContext`\[Xi]2^(-3) ($CellContext`\[Sigma]1 - $CellContext`\
\[Sigma]2)^2 $CellContext`t0[0, 3, 0] + 
       E^(-$CellContext`a1 + 2 $CellContext`a2) $CellContext`t^(
          1 - $CellContext`\[Sigma]1 + 
          2 $CellContext`\[Sigma]2) $CellContext`\[Xi]1^(-1) $CellContext`\
\[Xi]2^2 (1 - $CellContext`\[Sigma]1 + 
          2 $CellContext`\[Sigma]2)^2 $CellContext`t0[
         1, 0, 2] + $CellContext`t $CellContext`t0[1, 1, 1] + 
       E^($CellContext`a1 - 2 $CellContext`a2) $CellContext`t^(
          1 + $CellContext`\[Sigma]1 - 
          2 $CellContext`\[Sigma]2) $CellContext`\[Xi]1 \
$CellContext`\[Xi]2^(-2) (1 + $CellContext`\[Sigma]1 - 
          2 $CellContext`\[Sigma]2)^2 $CellContext`t0[1, 2, 0] + 
       E^((-2) $CellContext`a1 + $CellContext`a2) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1 + $CellContext`\[Sigma]2) $CellContext`\
\[Xi]1^(-2) $CellContext`\[Xi]2 (2 - 
          2 $CellContext`\[Sigma]1 + $CellContext`\[Sigma]2)^2 \
$CellContext`t0[2, 0, 1] + 
       E^(-$CellContext`a1 - $CellContext`a2) $CellContext`t^(
          2 - $CellContext`\[Sigma]1 - $CellContext`\[Sigma]2) $CellContext`\
\[Xi]1^(-1) $CellContext`\[Xi]2^(-1) (
          2 - $CellContext`\[Sigma]1 - $CellContext`\[Sigma]2)^2 \
$CellContext`t0[2, 1, 0] + 
       E^((-3) $CellContext`a1) $CellContext`t^(3 - 
          3 $CellContext`\[Sigma]1) $CellContext`\[Xi]1^(-3) (3 - 
          3 $CellContext`\[Sigma]1)^2 $CellContext`t0[3, 0, 0]))}, 0, 4, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.798275389338586*^9},
 CellLabel->
  "Out[260]=",ExpressionUUID->"2d8b7efa-8a49-4fe1-869e-dbecb0fdb317"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"Tau0", "[", "t", "]"}], 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Tau0", "[", "t", "]"}], "/.", 
         RowBox[{
          RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
          RowBox[{
           RowBox[{"p", "^", "2"}], " ", 
           RowBox[{"t", "^", "p"}]}]}]}], ")"}]}], "-", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Tau0", "[", "t", "]"}], "/.", 
          RowBox[{
           RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
           RowBox[{"p", " ", 
            RowBox[{"t", "^", "p"}]}]}]}], ")"}], "^", "2"}], ")"}]}], ")"}], 
    "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
     RowBox[{"(", "3", ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
     RowBox[{"(", "3", ")"}]}], "+", 
    RowBox[{
     RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
     RowBox[{"(", "3", ")"}]}]}], "//", "Simplify"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.798366831029063*^9, 3.798366890817655*^9}, {
  3.79836745060773*^9, 3.798367500218089*^9}, {3.798367614779303*^9, 
  3.7983676321802883`*^9}},
 CellLabel->
  "In[657]:=",ExpressionUUID->"7f672b3c-c1d5-4f5a-8ed9-6d7b8984f4c2"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"(", 
    InterpretationBox[
     RowBox[{
      InterpretationBox[
       SuperscriptBox[
        RowBox[{"O", "[", "\[Xi]2", "]"}], "3"],
       SeriesData[$CellContext`\[Xi]2, 0, {}, 3, 3, 1],
       Editable->False], "+", 
      InterpretationBox[
       SuperscriptBox[
        RowBox[{"O", "[", "\[Xi]1", "]"}], "3"],
       SeriesData[$CellContext`\[Xi]1, 0, {}, 0, 3, 1],
       Editable->False]}],
     SeriesData[$CellContext`\[Xi]1, 0, {
       SeriesData[$CellContext`\[Xi]2, 0, {}, 3, 3, 1]}, 0, 3, 1],
     Editable->False], ")"}], "+", 
   RowBox[{
    RowBox[{"(", 
     InterpretationBox[
      RowBox[{
       FractionBox[
        InterpretationBox[
         RowBox[{
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"-", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"1", "-", "\[Sigma]1"}]], " ", 
            SuperscriptBox["\[Sigma]1", "2"]}], 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"]], "+", 
          InterpretationBox[
           SuperscriptBox[
            RowBox[{"O", "[", "\[Xi]2", "]"}], "13"],
           SeriesData[$CellContext`\[Xi]2, 0, {}, 0, 13, 1],
           Editable->False]}],
         SeriesData[$CellContext`\[Xi]2, 0, {
          E^(-$CellContext`a1) $CellContext`t^(
             1 - $CellContext`\[Sigma]1) (-1 + $CellContext`\[Sigma]1)^(-2) \
$CellContext`\[Sigma]1^2}, 0, 13, 1],
         Editable->False], "\[Xi]1"], "+", 
       InterpretationBox[
        SuperscriptBox[
         RowBox[{"O", "[", "\[Xi]1", "]"}], "12"],
        SeriesData[$CellContext`\[Xi]1, 0, {}, -1, 12, 1],
        Editable->False]}],
      SeriesData[$CellContext`\[Xi]1, 0, {
        SeriesData[$CellContext`\[Xi]2, 0, {
         E^(-$CellContext`a1) $CellContext`t^(
            1 - $CellContext`\[Sigma]1) (-1 + $CellContext`\[Sigma]1)^(-2) \
$CellContext`\[Sigma]1^2}, 0, 13, 1]}, -1, 12, 1],
      Editable->False], ")"}], " ", "\[Epsilon]"}], "+", 
   RowBox[{
    RowBox[{"(", 
     InterpretationBox[
      RowBox[{
       FractionBox[
        InterpretationBox[
         RowBox[{
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{
              RowBox[{"-", "2"}], " ", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"2", "-", 
              RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"1", "-", 
               RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"]}], 
           RowBox[{"4", " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]], "+", 
          InterpretationBox[
           SuperscriptBox[
            RowBox[{"O", "[", "\[Xi]2", "]"}], "13"],
           SeriesData[$CellContext`\[Xi]2, 0, {}, 0, 13, 1],
           Editable->False]}],
         SeriesData[$CellContext`\[Xi]2, 0, {
          Rational[1, 4] 
           E^((-2) $CellContext`a1) $CellContext`t^(2 - 
             2 $CellContext`\[Sigma]1) (1 - 
             2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4)}, 
          0, 13, 1],
         Editable->False], 
        SuperscriptBox["\[Xi]1", "2"]], "+", 
       RowBox[{"(", 
        InterpretationBox[
         RowBox[{
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"-", "a2"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"1", "-", "\[Sigma]2"}]], " ", 
            SuperscriptBox["\[Sigma]2", "2"]}], 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"], " ", 
            "\[Xi]2"}]], "+", 
          InterpretationBox[
           SuperscriptBox[
            RowBox[{"O", "[", "\[Xi]2", "]"}], "12"],
           SeriesData[$CellContext`\[Xi]2, 0, {}, -1, 12, 1],
           Editable->False]}],
         SeriesData[$CellContext`\[Xi]2, 0, {
          E^(-$CellContext`a2) $CellContext`t^(
             1 - $CellContext`\[Sigma]2) ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2) $CellContext`\
\[Sigma]2^2}, -1, 12, 1],
         Editable->False], ")"}], "+", 
       InterpretationBox[
        SuperscriptBox[
         RowBox[{"O", "[", "\[Xi]1", "]"}], "11"],
        SeriesData[$CellContext`\[Xi]1, 0, {}, -2, 11, 1],
        Editable->False]}],
      SeriesData[$CellContext`\[Xi]1, 0, {
        SeriesData[$CellContext`\[Xi]2, 0, {
         Rational[1, 4] 
          E^((-2) $CellContext`a1) $CellContext`t^(2 - 
            2 $CellContext`\[Sigma]1) (1 - 
            2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4)}, 
         0, 13, 1], 0, 
        SeriesData[$CellContext`\[Xi]2, 0, {
         E^(-$CellContext`a2) $CellContext`t^(
            1 - $CellContext`\[Sigma]2) ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2) $CellContext`\
\[Sigma]2^2}, -1, 12, 1]}, -2, 11, 1],
      Editable->False], ")"}], " ", 
    SuperscriptBox["\[Epsilon]", "2"]}], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "3"],
    SeriesData[$CellContext`\[Epsilon], 0, {}, 0, 3, 1],
    Editable->False]}],
  SeriesData[$CellContext`\[Epsilon], 0, {
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {}, 3, 3, 1]}, 0, 3, 1], 
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {
       E^(-$CellContext`a1) $CellContext`t^(
          1 - $CellContext`\[Sigma]1) (-1 + $CellContext`\[Sigma]1)^(-2) \
$CellContext`\[Sigma]1^2}, 0, 13, 1]}, -1, 12, 1], 
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {
       Rational[1, 4] 
        E^((-2) $CellContext`a1) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1) (1 - 
          2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4)}, 0, 
       13, 1], 0, 
      SeriesData[$CellContext`\[Xi]2, 0, {
       E^(-$CellContext`a2) $CellContext`t^(
          1 - $CellContext`\[Sigma]2) ($CellContext`\[Sigma]1 - $CellContext`\
\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2) $CellContext`\[Sigma]2^2}, \
-1, 12, 1]}, -2, 11, 1]}, 0, 3, 1],
  Editable->False]], "Output",
 CellChangeTimes->{{3.798366893391548*^9, 3.798366919228387*^9}, 
   3.798367461545299*^9, {3.798367492765407*^9, 3.798367510782126*^9}, {
   3.798367623680216*^9, 3.798367634196685*^9}},
 CellLabel->
  "Out[657]=",ExpressionUUID->"f0f9a4dd-9b07-4eff-a862-21da1388e6e4"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"-", "1"}], ")"}], "^", 
       RowBox[{"(", 
        RowBox[{"k", "+", "1"}], ")"}]}], 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Tau0", "[", "t", "]"}], "-", "1", "+", 
          RowBox[{
           RowBox[{"O", "[", "\[Epsilon]", "]"}], "^", 
           RowBox[{"(", "3", ")"}]}], "+", 
          RowBox[{
           RowBox[{"O", "[", "\[Xi]1", "]"}], "^", 
           RowBox[{"(", "3", ")"}]}], "+", 
          RowBox[{
           RowBox[{"O", "[", "\[Xi]2", "]"}], "^", 
           RowBox[{"(", "3", ")"}]}]}], ")"}], "^", "k"}], "/", "k"}]}], ",", 
     
     RowBox[{"{", 
      RowBox[{"k", ",", "1", ",", "2"}], "}"}]}], "]"}], ")"}], "/.", 
  RowBox[{
   RowBox[{"t", "^", "p_"}], "\[RuleDelayed]", " ", 
   RowBox[{
    RowBox[{"p", "^", "2"}], " ", 
    RowBox[{"t", "^", "p"}]}]}]}]], "Input",
 CellChangeTimes->{{3.798366983673348*^9, 3.798367022315482*^9}, {
   3.7983672173005667`*^9, 3.7983672267007303`*^9}, {3.798367399073*^9, 
   3.7983674138481007`*^9}, 3.798367540576016*^9},
 CellLabel->
  "In[655]:=",ExpressionUUID->"7f69390b-1e35-4163-997d-41456c8142b3"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   RowBox[{"(", 
    InterpretationBox[
     RowBox[{
      InterpretationBox[
       SuperscriptBox[
        RowBox[{"O", "[", "\[Xi]2", "]"}], "3"],
       SeriesData[$CellContext`\[Xi]2, 0, {}, 3, 3, 1],
       Editable->False], "+", 
      InterpretationBox[
       SuperscriptBox[
        RowBox[{"O", "[", "\[Xi]1", "]"}], "3"],
       SeriesData[$CellContext`\[Xi]1, 0, {}, 0, 3, 1],
       Editable->False]}],
     SeriesData[$CellContext`\[Xi]1, 0, {
       SeriesData[$CellContext`\[Xi]2, 0, {}, 3, 3, 1]}, 0, 3, 1],
     Editable->False], ")"}], "+", 
   RowBox[{
    RowBox[{"(", 
     InterpretationBox[
      RowBox[{
       FractionBox[
        InterpretationBox[
         RowBox[{
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"-", "a1"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"1", "-", "\[Sigma]1"}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"1", "-", "\[Sigma]1"}], ")"}], "2"]}], 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "2"]], "+", 
          InterpretationBox[
           SuperscriptBox[
            RowBox[{"O", "[", "\[Xi]2", "]"}], "3"],
           SeriesData[$CellContext`\[Xi]2, 0, {}, 0, 3, 1],
           Editable->False]}],
         SeriesData[$CellContext`\[Xi]2, 0, {
          E^(-$CellContext`a1) $CellContext`t^(1 - $CellContext`\[Sigma]1) (
             1 - $CellContext`\[Sigma]1)^2 (-1 + \
$CellContext`\[Sigma]1)^(-2)}, 0, 3, 1],
         Editable->False], "\[Xi]1"], "+", 
       InterpretationBox[
        SuperscriptBox[
         RowBox[{"O", "[", "\[Xi]1", "]"}], "2"],
        SeriesData[$CellContext`\[Xi]1, 0, {}, -1, 2, 1],
        Editable->False]}],
      SeriesData[$CellContext`\[Xi]1, 0, {
        SeriesData[$CellContext`\[Xi]2, 0, {
         E^(-$CellContext`a1) $CellContext`t^(1 - $CellContext`\[Sigma]1) (
            1 - $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-2)},
          0, 3, 1]}, -1, 2, 1],
      Editable->False], ")"}], " ", "\[Epsilon]"}], "+", 
   RowBox[{
    RowBox[{"(", 
     InterpretationBox[
      RowBox[{
       FractionBox[
        InterpretationBox[
         RowBox[{
          RowBox[{"-", 
           FractionBox[
            RowBox[{
             SuperscriptBox["\[ExponentialE]", 
              RowBox[{
               RowBox[{"-", "2"}], " ", "a1"}]], " ", 
             SuperscriptBox["t", 
              RowBox[{"2", "-", 
               RowBox[{"2", " ", "\[Sigma]1"}]}]], " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"2", "-", 
                RowBox[{"2", " ", "\[Sigma]1"}]}], ")"}], "2"]}], 
            RowBox[{"4", " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "1"}], "+", "\[Sigma]1"}], ")"}], "4"]}]]}], "+", 
          InterpretationBox[
           SuperscriptBox[
            RowBox[{"O", "[", "\[Xi]2", "]"}], "3"],
           SeriesData[$CellContext`\[Xi]2, 0, {}, 0, 3, 1],
           Editable->False]}],
         SeriesData[$CellContext`\[Xi]2, 0, {
          Rational[-1, 4] 
           E^((-2) $CellContext`a1) $CellContext`t^(2 - 
             2 $CellContext`\[Sigma]1) (2 - 
             2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4)}, 
          0, 3, 1],
         Editable->False], 
        SuperscriptBox["\[Xi]1", "2"]], "+", 
       RowBox[{"(", 
        InterpretationBox[
         RowBox[{
          FractionBox[
           RowBox[{
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"-", "a2"}]], " ", 
            SuperscriptBox["t", 
             RowBox[{"1", "-", "\[Sigma]2"}]], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"1", "-", "\[Sigma]2"}], ")"}], "2"]}], 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"\[Sigma]1", "-", "\[Sigma]2"}], ")"}], "2"], " ", 
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", "\[Sigma]2"}], ")"}], "2"], " ", 
            "\[Xi]2"}]], "+", 
          InterpretationBox[
           SuperscriptBox[
            RowBox[{"O", "[", "\[Xi]2", "]"}], "2"],
           SeriesData[$CellContext`\[Xi]2, 0, {}, -1, 2, 1],
           Editable->False]}],
         SeriesData[$CellContext`\[Xi]2, 0, {
          E^(-$CellContext`a2) $CellContext`t^(1 - $CellContext`\[Sigma]2) (
             1 - $CellContext`\[Sigma]2)^2 ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2)}, -1, 2, 1],
         Editable->False], ")"}], "+", 
       InterpretationBox[
        SuperscriptBox[
         RowBox[{"O", "[", "\[Xi]1", "]"}], "1"],
        SeriesData[$CellContext`\[Xi]1, 0, {}, -2, 1, 1],
        Editable->False]}],
      SeriesData[$CellContext`\[Xi]1, 0, {
        SeriesData[$CellContext`\[Xi]2, 0, {
         Rational[-1, 4] 
          E^((-2) $CellContext`a1) $CellContext`t^(2 - 
            2 $CellContext`\[Sigma]1) (2 - 
            2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4)}, 
         0, 3, 1], 0, 
        SeriesData[$CellContext`\[Xi]2, 0, {
         E^(-$CellContext`a2) $CellContext`t^(1 - $CellContext`\[Sigma]2) (
            1 - $CellContext`\[Sigma]2)^2 ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2)}, -1, 2, 
         1]}, -2, 1, 1],
      Editable->False], ")"}], " ", 
    SuperscriptBox["\[Epsilon]", "2"]}], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "\[Epsilon]", "]"}], "3"],
    SeriesData[$CellContext`\[Epsilon], 0, {}, 0, 3, 1],
    Editable->False]}],
  SeriesData[$CellContext`\[Epsilon], 0, {
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {}, 3, 3, 1]}, 0, 3, 1], 
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {
       E^(-$CellContext`a1) $CellContext`t^(1 - $CellContext`\[Sigma]1) (
          1 - $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-2)}, 
       0, 3, 1]}, -1, 2, 1], 
    SeriesData[$CellContext`\[Xi]1, 0, {
      SeriesData[$CellContext`\[Xi]2, 0, {
       Rational[-1, 4] 
        E^((-2) $CellContext`a1) $CellContext`t^(2 - 
          2 $CellContext`\[Sigma]1) (2 - 
          2 $CellContext`\[Sigma]1)^2 (-1 + $CellContext`\[Sigma]1)^(-4)}, 0, 
       3, 1], 0, 
      SeriesData[$CellContext`\[Xi]2, 0, {
       E^(-$CellContext`a2) $CellContext`t^(1 - $CellContext`\[Sigma]2) (
          1 - $CellContext`\[Sigma]2)^2 ($CellContext`\[Sigma]1 - \
$CellContext`\[Sigma]2)^(-2) (-1 + $CellContext`\[Sigma]2)^(-2)}, -1, 2, 
       1]}, -2, 1, 1]}, 0, 3, 1],
  Editable->False]], "Output",
 CellChangeTimes->{
  3.79836722932701*^9, 3.7983672645253477`*^9, {3.798367399973819*^9, 
   3.798367414471375*^9}, 3.798367541253806*^9},
 CellLabel->
  "Out[655]=",ExpressionUUID->"a6ad2b89-0c37-4bb0-8ae0-107d50f5d5c3"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(", 
  RowBox[{
   RowBox[{"Tau0", "[", "t", "]"}], "-", "1"}], ")"}]], "Input",
 CellChangeTimes->{{3.798367231995789*^9, 
  3.798367257329693*^9}},ExpressionUUID->"e02caf23-695b-4ca1-a31c-\
30938ec4117b"]
},
WindowSize->{1920, 1055},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
Magnification:>1.5 Inherited,
FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 8621, 222, 490, "Input",ExpressionUUID->"45f0ad4a-754f-41f8-a483-c077c6705864"],
Cell[9182, 244, 6609, 194, 354, "Input",ExpressionUUID->"dc3effc5-390c-4b5c-9f6f-981959ddecaa"],
Cell[15794, 440, 380, 10, 47, "Input",ExpressionUUID->"c3677654-05fe-4766-9219-6586698020e9"],
Cell[16177, 452, 5775, 170, 306, "Input",ExpressionUUID->"89f1c56e-4e87-497f-94bb-3925c7acf9c2"],
Cell[21955, 624, 1229, 29, 82, "Input",ExpressionUUID->"c4dc56d4-0065-4260-bff8-81b8f2387be3"],
Cell[CellGroupData[{
Cell[23209, 657, 1137, 30, 47, "Input",ExpressionUUID->"abf2448e-47fa-4ddf-9770-0706af56979c"],
Cell[24349, 689, 236, 5, 53, "Output",ExpressionUUID->"78d3033e-65cb-42ac-bd88-f589928a8f24"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24622, 699, 954, 20, 47, "Input",ExpressionUUID->"d333b561-ecda-412a-98dc-2a212e8a39bb"],
Cell[25579, 721, 852, 19, 81, "Output",ExpressionUUID->"b490ca9c-8c71-4499-8d9b-743b884e9457"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26468, 745, 326, 8, 47, "Input",ExpressionUUID->"84a9803d-7eb4-4040-a97e-d089f6f4b672"],
Cell[26797, 755, 234, 4, 53, "Output",ExpressionUUID->"281a22b1-6efe-46de-9acd-50699c53d2e7"]
}, Open  ]],
Cell[27046, 762, 152, 3, 47, "Input",ExpressionUUID->"14e70c01-de81-4b68-8c97-77c71c3eec30"],
Cell[CellGroupData[{
Cell[27223, 769, 303, 8, 47, "Input",ExpressionUUID->"e659888a-46d1-4dcb-b09a-99b9bc05161b"],
Cell[27529, 779, 2510, 78, 89, "Output",ExpressionUUID->"6e780e73-6fdf-446d-be3b-c2b7f24c3c2a"]
}, Open  ]],
Cell[30054, 860, 158, 3, 47, "Input",ExpressionUUID->"e8a00397-aef7-4303-94cc-3e1f0f0ca554"],
Cell[30215, 865, 2155, 67, 116, "Input",ExpressionUUID->"e3bc3cd1-d9c0-4245-bcce-faef4b39dce9"],
Cell[32373, 934, 179, 4, 47, "Input",ExpressionUUID->"540b1cc2-38c0-48db-9004-b39b156bb7d7"],
Cell[32555, 940, 3927, 122, 558, "Input",ExpressionUUID->"69c1f79b-4733-4927-a510-0c734b2ee0c5"],
Cell[CellGroupData[{
Cell[36507, 1066, 1851, 50, 116, "Input",ExpressionUUID->"5e004098-772f-467a-92d1-a6d9aecacf72"],
Cell[38361, 1118, 420, 11, 81, "Output",ExpressionUUID->"86cd1d40-c155-4f35-9a8d-0a9483261ab1"]
}, Open  ]],
Cell[38796, 1132, 7295, 217, 354, "Input",ExpressionUUID->"27366b93-7abd-4206-bdb0-2a3a75d81e3a"],
Cell[46094, 1351, 4431, 131, 272, "Input",ExpressionUUID->"3a0fe918-4fc7-4e0c-9e51-0767004f7db9"],
Cell[50528, 1484, 1061, 28, 82, "Input",ExpressionUUID->"569e09a5-261a-4388-88a9-94f48502493f"],
Cell[CellGroupData[{
Cell[51614, 1516, 878, 28, 47, "Input",ExpressionUUID->"77ae1bbe-bae5-469a-95ef-b7791b3752f1"],
Cell[52495, 1546, 19886, 490, 637, "Output",ExpressionUUID->"2d8b7efa-8a49-4fe1-869e-dbecb0fdb317"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72418, 2041, 1385, 40, 150, "Input",ExpressionUUID->"7f672b3c-c1d5-4f5a-8ed9-6d7b8984f4c2"],
Cell[73806, 2083, 6937, 174, 122, "Output",ExpressionUUID->"f0f9a4dd-9b07-4eff-a862-21da1388e6e4"]
}, Open  ]],
Cell[CellGroupData[{
Cell[80780, 2262, 1271, 37, 47, "Input",ExpressionUUID->"7f69390b-1e35-4163-997d-41456c8142b3"],
Cell[82054, 2301, 7112, 178, 122, "Output",ExpressionUUID->"a6ad2b89-0c37-4bb0-8ae0-107d50f5d5c3"]
}, Open  ]],
Cell[89181, 2482, 230, 6, 47, "Input",ExpressionUUID->"e02caf23-695b-4ca1-a31c-30938ec4117b"]
}
]
*)