NonlinearElasticity.nb

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

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

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    135580,       2950]
NotebookOptionsPosition[    131344,       2871]
NotebookOutlinePosition[    131709,       2887]
CellTagsIndexPosition[    131666,       2884]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Elasticity", "Section",
 CellChangeTimes->{{3.7712280935211678`*^9, 
  3.7712280972140417`*^9}},ExpressionUUID->"1db9d287-f7c6-42b4-a9d4-\
606839815e86"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Inv1", "[", "tensor_", "]"}], ":=", 
   RowBox[{"Tr", "[", "tensor", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Inv2", "[", "tensor_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"Tr", "[", "tensor", "]"}], ")"}], "2"], "-", 
      RowBox[{"Tr", "[", 
       RowBox[{"tensor", ".", 
        RowBox[{"Transpose", "[", "tensor", "]"}]}], "]"}]}], ")"}], "/", 
    "2"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Inv3", "[", "tensor_", "]"}], ":=", 
    RowBox[{"Det", "[", "tensor", "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"StrainGreen", "[", "mode_", "]"}], ":=", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", "SGr", "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"SGr", "=", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Transpose", "[", "Deform", "]"}], "+", "Deform"}], ")"}], 
         "/", "2"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"mode", "\[Equal]", "\"\<nonlinear2\>\""}], "||", 
          RowBox[{"mode", "\[Equal]", "\"\<nonlinear3\>\""}]}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"SGr", "=", 
          RowBox[{"SGr", "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"Transpose", "[", "Deform", "]"}], ".", "Deform"}], 
             ")"}], "/", "2"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", 
       RowBox[{"SGrTemplate", "=", 
        RowBox[{"{", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            SubscriptBox["SG", "xx"], ",", 
            SubscriptBox["SG", "xy"], ",", 
            SubscriptBox["SG", "xz"]}], "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            SubscriptBox["SG", "yx"], ",", 
            SubscriptBox["SG", "yy"], ",", 
            SubscriptBox["SG", "yz"]}], "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            SubscriptBox["SG", "zx"], ",", 
            SubscriptBox["SG", "zy"], ",", 
            SubscriptBox["SG", "zz"]}], "}"}]}], "}"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"SGrSubst", "=", 
        RowBox[{"{", " ", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["SG", "xx"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}], ",", 
          RowBox[{
           SubscriptBox["SG", "xy"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}], ",", 
          RowBox[{
           SubscriptBox["SG", "xz"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"1", ",", "3"}], "\[RightDoubleBracket]"}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{
           SubscriptBox["SG", "yx"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}], ",", 
          RowBox[{
           SubscriptBox["SG", "yy"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"2", ",", "2"}], "\[RightDoubleBracket]"}]}], ",", 
          RowBox[{
           SubscriptBox["SG", "yz"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"2", ",", "3"}], "\[RightDoubleBracket]"}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{
           SubscriptBox["SG", "zx"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"3", ",", "1"}], "\[RightDoubleBracket]"}]}], ",", 
          RowBox[{
           SubscriptBox["SG", "zy"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"3", ",", "2"}], "\[RightDoubleBracket]"}]}], ",", 
          RowBox[{
           SubscriptBox["SG", "zz"], "\[Rule]", 
           RowBox[{"SGr", "\[LeftDoubleBracket]", 
            RowBox[{"3", ",", "3"}], "\[RightDoubleBracket]"}]}]}], "}"}]}], 
       ";", "\[IndentingNewLine]", 
       RowBox[{"I1", " ", "=", " ", 
        RowBox[{"Inv1", "[", "SGrTemplate", "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"I2", " ", "=", " ", 
        RowBox[{"Inv2", "[", "SGrTemplate", "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"I3", " ", "=", " ", 
        RowBox[{"Inv3", "[", "SGrTemplate", "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"E1", " ", "=", " ", 
        RowBox[{"Inv1", "[", "SGrTemplate", "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"E2", " ", "=", " ", 
        RowBox[{"Inv1", "[", 
         RowBox[{"SGrTemplate", ".", "SGrTemplate"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"E3", " ", "=", " ", 
        RowBox[{"Inv1", "[", 
         RowBox[{"SGrTemplate", ".", "SGrTemplate", ".", "SGrTemplate"}], 
         "]"}]}], ";", "\[IndentingNewLine]", "SGr"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"PotEnergyTempl", "[", 
     RowBox[{"mode_", ",", "modules_"}], "]"}], ":=", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", "Pot", "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"Pot", "=", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"\[Lambda]", "+", 
             RowBox[{"2", " ", "\[Mu]"}]}], ")"}], "/", "2"}], " ", 
          SuperscriptBox["I1", "2"]}], "-", 
         RowBox[{"2", " ", "\[Mu]", " ", "I2"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"modules", "\[Equal]", "\"\<murnaghan\>\""}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"mode", "\[Equal]", "\"\<nonlinear2\>\""}], "||", 
            RowBox[{"mode", "\[Equal]", "\"\<nonlinear3\>\""}]}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"Pot", "=", 
             RowBox[{"Pot", "+", 
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"l", "+", 
                  RowBox[{"2", " ", "m"}]}], ")"}], "/", "3"}], " ", 
               SuperscriptBox["I1", "3"]}], "-", 
              RowBox[{"2", " ", "m", " ", "I1", " ", "I2"}], "+", 
              RowBox[{"n", " ", "I3"}]}]}], ";", "\[IndentingNewLine]", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{"mode", "\[Equal]", "\"\<nonlinear3\>\""}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{"Pot", "=", 
               RowBox[{"Pot", "+", 
                RowBox[{"g1", " ", 
                 SuperscriptBox["I1", "4"]}], "+", 
                RowBox[{"g2", " ", 
                 SuperscriptBox["I2", "2"]}], "+", 
                RowBox[{"g3", " ", "I1", " ", "I3"}], "+", 
                RowBox[{"g4", " ", 
                 SuperscriptBox["I2", "2"]}]}]}]}], "]"}], ";"}]}], "]"}], 
         ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"modules", "\[Equal]", "\"\<landau\>\""}], ",", 
            "\[IndentingNewLine]", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"mode", "\[Equal]", "\"\<nonlinear2\>\""}], "||", 
               RowBox[{"mode", "\[Equal]", "\"\<nonlinear3\>\""}]}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"Pot", "=", 
                RowBox[{"Pot", "+", 
                 RowBox[{
                  RowBox[{"CL", "/", "3"}], " ", 
                  SuperscriptBox["E1", "3"]}], "+", 
                 RowBox[{"BL", " ", "E1", " ", "E2"}], "+", 
                 RowBox[{
                  RowBox[{"AL", "/", "3"}], " ", "E3"}]}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{"If", "[", 
                RowBox[{
                 RowBox[{"mode", "\[Equal]", "\"\<nonlinear3\>\""}], ",", 
                 "\[IndentingNewLine]", 
                 RowBox[{"Pot", "=", 
                  RowBox[{"Pot", "+", 
                   RowBox[{"HL", " ", 
                    SuperscriptBox["E1", "4"]}], "+", 
                   RowBox[{"GL", " ", 
                    SuperscriptBox["E2", "2"]}], "+", 
                   RowBox[{"FL", " ", 
                    SuperscriptBox["E1", "2"], " ", "E2"}], "+", 
                   RowBox[{"EL", " ", "E1", " ", "E3"}]}]}]}], "]"}], ";"}]}],
              "\[IndentingNewLine]", "]"}], ",", 
            RowBox[{"Throw", "[", "\"\<Unknown modules\>\"", "]"}]}], "]"}], 
          ";"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       "Pot"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"PotEnergy", "[", 
     RowBox[{"mode_", ",", "modules_"}], "]"}], ":=", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"PotEnergyTempl", "[", 
        RowBox[{"mode", ",", "modules"}], "]"}], "/.", "SGrSubst"}]}], 
     "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Piola", "[", 
    RowBox[{"mode_", ",", "modules_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"P", ",", "Pot"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"mode", "\[Equal]", "\"\<linear\>\""}], ",", 
        RowBox[{"P", "=", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Lambda]", " ", 
            RowBox[{"Inv1", "[", "SGr", "]"}], 
            RowBox[{"IdentityMatrix", "[", "3", "]"}]}], "+", 
           RowBox[{"2", "\[Mu]", " ", "SGr"}]}], "//", "Simplify"}]}]}], 
       "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"mode", "\[Equal]", "\"\<nonlinear2\>\""}], "||", 
         RowBox[{"mode", "\[Equal]", "\"\<nonlinear3\>\""}]}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"Pot", "=", 
          RowBox[{"PotEnergyTempl", "[", 
           RowBox[{"mode", ",", "modules"}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"P", "=", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"IdentityMatrix", "[", "3", "]"}], "+", "Deform"}], 
             ")"}], ".", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "xx"]], "Pot"}], ",", 
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "xy"]], "Pot"}], ",", 
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "xz"]], "Pot"}]}], "}"}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "yx"]], "Pot"}], ",", 
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "yy"]], "Pot"}], ",", 
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "yz"]], "Pot"}]}], "}"}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "zx"]], "Pot"}], ",", 
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "zy"]], "Pot"}], ",", 
                RowBox[{
                 SubscriptBox["\[PartialD]", 
                  SubscriptBox["SG", "zz"]], "Pot"}]}], "}"}]}], "}"}]}], "/.",
            "SGrSubst"}]}]}]}], " ", "]"}], ";", "\[IndentingNewLine]", 
      "P"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.763533502188644*^9, 3.763533645609955*^9}, {
   3.763533675934937*^9, 3.763533697803871*^9}, {3.76353375402313*^9, 
   3.763533759890518*^9}, 3.7635338525281982`*^9, {3.763533900189554*^9, 
   3.7635339125773478`*^9}, {3.763533953085806*^9, 3.763533966554759*^9}, {
   3.7635340089303017`*^9, 3.763534144893711*^9}, {3.763534262260969*^9, 
   3.7635343005338383`*^9}, {3.76353440267231*^9, 3.763534420094318*^9}, {
   3.763534528815589*^9, 3.763534661779916*^9}, {3.763535083521009*^9, 
   3.763535095858418*^9}, {3.7643917518930893`*^9, 3.764391833009976*^9}, {
   3.764391916006719*^9, 3.764392015074498*^9}, {3.764392109122086*^9, 
   3.764392116003626*^9}, {3.766990195127303*^9, 3.766990221889421*^9}, {
   3.766990876832941*^9, 3.766990908343811*^9}, {3.7669909502903233`*^9, 
   3.766990959679263*^9}, {3.766990992128806*^9, 3.766991082245723*^9}, {
   3.766991932729607*^9, 3.766992021058549*^9}, {3.770386382492463*^9, 
   3.770386410150251*^9}, {3.7703865056235447`*^9, 3.770386508093234*^9}, {
   3.7703867564157457`*^9, 3.770386910553578*^9}, {3.770386940657357*^9, 
   3.770386965262755*^9}, {3.7703870392363377`*^9, 3.770387135754861*^9}, {
   3.770387194648162*^9, 3.7703871979483147`*^9}, {3.770387245749598*^9, 
   3.7703872499020157`*^9}, {3.770387319258977*^9, 3.770387418338616*^9}, {
   3.770387467879286*^9, 3.770387521327503*^9}, {3.770387588591833*^9, 
   3.770387594335718*^9}, {3.7703876300628223`*^9, 3.7703876614948893`*^9}, {
   3.770387822804311*^9, 3.770387875986969*^9}, {3.772335152773645*^9, 
   3.772335155207932*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"10bc8bc5-de5e-4ac8-abc4-9d37283af475"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Uniaxial tension", "Section",
 CellChangeTimes->{{3.772250066752474*^9, 3.772250072510826*^9}, {
  3.7722503029509087`*^9, 
  3.772250307135228*^9}},ExpressionUUID->"d669c3bd-62da-412e-88a6-\
05ff878fb8e6"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ClearAll", "[", 
   RowBox[{"U", ",", "V", ",", "W", ",", "\[Lambda]", ",", "\[Mu]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"U", "=", 
   RowBox[{"1", "+", "d"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Deform", "=", 
   RowBox[{
    RowBox[{"1", "/", "2"}], " ", 
    RowBox[{"DiagonalMatrix", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        SuperscriptBox["U", "2"], "-", "1"}], ",", 
       RowBox[{
        RowBox[{"1", "/", "U"}], "-", "1"}], ",", 
       RowBox[{
        RowBox[{"1", "/", "U"}], "-", "1"}]}], "}"}], "]"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.731858722647919*^9, 3.731858725748151*^9}, {
   3.7318589710473223`*^9, 3.731858971707371*^9}, {3.732815489305831*^9, 
   3.7328154913058605`*^9}, {3.761980579692844*^9, 3.76198058281301*^9}, {
   3.763532776472636*^9, 3.763532822881721*^9}, {3.7635328694929647`*^9, 
   3.763532901769594*^9}, {3.76353295740184*^9, 3.7635329599238462`*^9}, {
   3.7635330693416233`*^9, 3.7635330716142406`*^9}, {3.763534350899652*^9, 
   3.7635343516569033`*^9}, {3.763534409791201*^9, 3.7635344099618387`*^9}, 
   3.7635347029868*^9, {3.763535035488393*^9, 3.763535036165593*^9}, {
   3.763535080693055*^9, 3.763535081585973*^9}, {3.764437271436191*^9, 
   3.764437343746107*^9}, {3.7722501474044228`*^9, 3.772250204055682*^9}, {
   3.772250488205019*^9, 3.7722505259515667`*^9}, {3.772250567241746*^9, 
   3.7722505700247993`*^9}, {3.772251935730007*^9, 3.772251936425785*^9}, {
   3.772252010209832*^9, 3.772252019803437*^9}, {3.772252055827446*^9, 
   3.772252071600865*^9}, {3.77225239471063*^9, 3.7722524089597607`*^9}, {
   3.772252554184457*^9, 3.772252563588551*^9}, {3.772253290471219*^9, 
   3.77225329130687*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"118ca0f4-45d3-4e90-9b82-736b8a259fed"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "PMMA", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"AL", "=", 
     RowBox[{
      RowBox[{"-", "1.41"}], "*", 
      RowBox[{"10", "^", "9"}], " ", "a"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BL", "=", 
     RowBox[{
      RowBox[{"-", "7.02"}], "*", 
      RowBox[{"10", "^", "9"}], " ", "a"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"CL", "=", 
     RowBox[{
      RowBox[{"-", "3.91"}], "*", 
      RowBox[{"10", "^", "9"}], " ", "a"}]}], ";"}], "\[IndentingNewLine]", 
   StyleBox[
    RowBox[{
     RowBox[{"\[Lambda]", "=", 
      RowBox[{"\[CapitalEpsilon]", " ", 
       RowBox[{
        RowBox[{"\[Nu]", "/", 
         RowBox[{"(", 
          RowBox[{"1", "+", "\[Nu]"}], ")"}]}], "/", 
        RowBox[{"(", 
         RowBox[{"1", "-", 
          RowBox[{"2", "\[Nu]"}]}], ")"}]}]}]}], ";"}],
    FontFamily->"Courier New"], "\[IndentingNewLine]", 
   StyleBox[
    RowBox[{
     RowBox[{"\[Mu]", "=", 
      RowBox[{
       RowBox[{"\[CapitalEpsilon]", "/", "2"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", "\[Nu]"}], ")"}]}]}], ";"}],
    FontFamily->"Courier New"], 
   StyleBox["\[IndentingNewLine]",
    FontFamily->"Courier New"], 
   StyleBox[
    RowBox[{
     RowBox[{"\[CapitalEpsilon]", "=", 
      RowBox[{"4.92", "*", 
       RowBox[{"10", "^", "9"}]}]}], ";"}],
    FontFamily->"Courier New"], 
   StyleBox["\[IndentingNewLine]",
    FontFamily->"Courier New"], 
   StyleBox[
    RowBox[{
     RowBox[{"\[Nu]", "=", "0.34"}], ";"}],
    FontFamily->"Courier New"]}]}]], "Input",
 CellChangeTimes->{{3.772250791858264*^9, 3.7722508231358643`*^9}, {
  3.772250935288672*^9, 3.772250967258156*^9}, {3.772254921766569*^9, 
  3.77225492651581*^9}, {3.7722583137147293`*^9, 3.7722583354589453`*^9}, {
  3.772259592510427*^9, 3.772259594020467*^9}, {3.772335169103757*^9, 
  3.772335175403378*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"ae2879c8-059f-423b-9b96-ae967dadb5c8"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"mode", "=", "\"\<nonlinear3\>\""}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"modules", "=", "\"\<landau\>\""}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"SGr", "=", 
   RowBox[{"StrainGreen", "[", "mode", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Pot", "=", 
   RowBox[{"PotEnergy", "[", 
    RowBox[{"mode", ",", "modules"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"P", "=", 
   RowBox[{"Piola", "[", 
    RowBox[{"mode", ",", "modules"}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.772250110270763*^9, 3.772250110342271*^9}, {
  3.772253464127849*^9, 3.7722534652925653`*^9}, {3.772253758169621*^9, 
  3.772253761981372*^9}, {3.772254680482554*^9, 3.772254681427102*^9}, {
  3.772254746969104*^9, 3.772254749365308*^9}, {3.7723352055471354`*^9, 
  3.7723352057132673`*^9}},
 CellLabel->"In[18]:=",ExpressionUUID->"faa411b5-4659-4640-a167-295e5a203224"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Stress", "[", 
    RowBox[{"d_", ",", "a_", ",", "EL_", ",", "FL_", ",", "GL_", ",", "HL_"}],
     "]"}], "=", 
   RowBox[{
    RowBox[{"P", "[", 
     RowBox[{"[", 
      RowBox[{"1", ",", "1"}], "]"}], "]"}], "//", "Simplify"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Stress", "[", 
    RowBox[{"d_", ",", "a_"}], "]"}], "=", 
   RowBox[{"Stress", "[", 
    RowBox[{"d", ",", "a", ",", "0", ",", "0", ",", "0", ",", "0"}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.7722505727395*^9, 3.772250600126786*^9}, {
   3.772250998324006*^9, 3.772251001080229*^9}, {3.772252091694337*^9, 
   3.77225209281388*^9}, {3.772252414165577*^9, 3.772252414305666*^9}, {
   3.7722525462375937`*^9, 3.7722525465556173`*^9}, 3.772257729017083*^9, 
   3.7722577610978107`*^9, 3.7722592327297*^9, {3.772335192933798*^9, 
   3.7723352535565233`*^9}, {3.772350861428997*^9, 3.7723508876266813`*^9}},
 CellLabel->"In[23]:=",ExpressionUUID->"1e414e05-8461-40ee-917c-b603a07c0f19"],

Cell[CellGroupData[{

Cell["Extension", "Subsection",
 CellChangeTimes->{{3.7722768112339354`*^9, 
  3.772276812696094*^9}},ExpressionUUID->"7289dff9-9280-4a82-94d7-\
3c7a33b67c65"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Data", " ", "from", " ", "Okeke", " ", "2017"}], ":", " ", 
    RowBox[{"quasi", "-", 
     RowBox[{"static", " ", "tension"}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ExperDataStaticExt", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.005", ",", 
         RowBox[{"16.6", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.01", ",", 
         RowBox[{"29", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.015", ",", 
         RowBox[{"40", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"0.02", ",", 
         RowBox[{"48", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.025", ",", 
         RowBox[{"54", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.03", ",", 
         RowBox[{"57", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ExperDataStaticExtInterp", "=", 
     RowBox[{"Interpolation", "[", "ExperDataStaticExt", "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"Data", " ", "from", " ", "Wu", " ", "2004"}], ":", " ", 
     RowBox[{
     "intermediate", " ", "rate", " ", "of", " ", "tension", " ", "2.81"}]}], 
    " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ExperDataMedDynExt", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.00129", ",", 
         RowBox[{"4.95", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0025", ",", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"10", "-", 
            RowBox[{"5", "/", "12"}]}], ")"}], "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"0.00375", ",", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"15", "-", 
            RowBox[{"5", "/", "6"}]}], ")"}], "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.005", ",", 
         RowBox[{"18.75", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.00625", ",", 
         RowBox[{"23", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"0.0075", ",", 
         RowBox[{"27.5", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0086", ",", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"30", "+", 
            RowBox[{"5", "/", "3"}]}], ")"}], "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"0.01", "-", 
          RowBox[{"0.0025", "/", "8"}]}], ",", 
         RowBox[{"35.5", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"0.01", "+", 
          RowBox[{"0.0075", "/", "8"}]}], ",", 
         RowBox[{"40", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ExperDataMedDynExtInterp", "=", 
     RowBox[{"Interpolation", "[", "ExperDataMedDynExt", "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ExperDataMedDynExtSmooth", "[", "x_", "]"}], "=", 
     RowBox[{"Fit", "[", 
      RowBox[{"ExperDataMedDynExt", ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         SuperscriptBox["x", "2"]}], "}"}], ",", "x"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"Data", " ", "from", " ", "Wu", " ", "2004"}], ":", " ", 
     RowBox[{
     "intermediate", " ", "rate", " ", "of", " ", "tension", " ", "0.292"}]}],
     " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ExperDataMedDynExt2", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.00129", ",", 
         RowBox[{"4.8", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0025", ",", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"10", "-", 
            RowBox[{"5", "/", "3"}]}], ")"}], "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"0.00375", ",", 
         RowBox[{"12.5", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.005", ",", 
         RowBox[{"16", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0062", ",", 
         RowBox[{"19.5", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.0075", ",", 
         RowBox[{"23", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"0.00875", ",", 
         RowBox[{"27", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.01", ",", 
         RowBox[{"31", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ExperDataMedDynExtInterp2", "=", 
     RowBox[{"Interpolation", "[", "ExperDataMedDynExt2", "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"ExperDataMedDynExtSmooth2", "[", "x_", "]"}], "=", 
     RowBox[{"Fit", "[", 
      RowBox[{"ExperDataMedDynExt2", ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         SuperscriptBox["x", "2"]}], "}"}], ",", "x"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"ListLinePlot", "[", 
      RowBox[{"ExperDataMedDynExt", ",", 
       RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}], ",", 
       RowBox[{"GridLines", "\[Rule]", "Automatic"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListLinePlot", "[", 
      RowBox[{"ExperDataMedDynExt2", ",", 
       RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}], ",", 
       RowBox[{"GridLines", "\[Rule]", "Automatic"}]}], "]"}]}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.772254876400169*^9, 3.772254909437199*^9}, {
   3.7722552252851152`*^9, 3.772255232063011*^9}, {3.772255271654971*^9, 
   3.7722552808311663`*^9}, {3.7722553186072683`*^9, 3.772255395782542*^9}, 
   3.772255765392796*^9, {3.772255911047604*^9, 3.7722559148713303`*^9}, {
   3.772270826910883*^9, 3.7722710717138577`*^9}, {3.772272770566276*^9, 
   3.772272774095048*^9}, {3.7722730188803368`*^9, 3.772273034550399*^9}, {
   3.77227646194468*^9, 3.772276469047099*^9}, {3.77227657098955*^9, 
   3.772276582033478*^9}, {3.77227668035983*^9, 3.772276692862258*^9}, {
   3.772335271636469*^9, 3.7723352917166986`*^9}, {3.7723483827284107`*^9, 
   3.772348438771613*^9}, {3.772349674231357*^9, 3.772349931089978*^9}, {
   3.77234996706243*^9, 3.77234997740195*^9}, {3.772350012296361*^9, 
   3.772350012935498*^9}, 3.77235015957345*^9, {3.7723502544078913`*^9, 
   3.772350274562233*^9}, {3.772350360563099*^9, 3.77235036074614*^9}, {
   3.772357692414012*^9, 3.772357722145774*^9}, {3.7723578220138702`*^9, 
   3.772357822120489*^9}, {3.7723578523259907`*^9, 3.772357852803117*^9}, {
   3.772357883832855*^9, 3.7723579679492817`*^9}, {3.772358009296788*^9, 
   3.7723580949642677`*^9}, {3.772358130630684*^9, 3.7723583272019997`*^9}, {
   3.772358590314192*^9, 3.772358679535317*^9}, {3.7723587136665583`*^9, 
   3.772358716083259*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"6334fc9b-4e81-4604-a00b-737281ea1177"],

Cell[BoxData[
 GraphicsBox[{{{}, {{{}, {}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
        NCache[
         Rational[7, 360], 0.019444444444444445`]], AbsoluteThickness[1.6], 
       LineBox[{{0., 0.}, {0.00129, 4.95*^6}, {0.0025, 
        9.583333333333334*^6}, {0.00375, 1.4166666666666666`*^7}, {0.005, 
        1.875*^7}, {0.00625, 2.3*^7}, {0.0075, 2.75*^7}, {0.0086, 
        3.1666666666666668`*^7}, {0.0096875, 3.55*^7}, {0.0109375, 
        4.*^7}}]}}, {
      {RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[6], 
       AbsoluteThickness[1.6], GeometricTransformationBox[InsetBox[
         StyleBox["\<\"\[FilledCircle]\"\>",
          StripOnInput->False,
          FontSize->8.96], {0., 0.}], {{{0., 0.}}, {{0.00129, 4.95*^6}}, {{
        0.0025, 9.583333333333334*^6}}, {{0.00375, 
        1.4166666666666666`*^7}}, {{0.005, 1.875*^7}}, {{0.00625, 2.3*^7}}, {{
        0.0075, 2.75*^7}}, {{0.0086, 3.1666666666666668`*^7}}, {{0.0096875, 
        3.55*^7}}, {{0.0109375, 
        4.*^7}}}]}, {}}}, {}, {}, {}, {}}, {{}, {{{}, {}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
        NCache[
         Rational[7, 360], 0.019444444444444445`]], AbsoluteThickness[1.6], 
       LineBox[{{0., 0.}, {0.00129, 4.8*^6}, {0.0025, 8.333333333333333*^6}, {
        0.00375, 1.25*^7}, {0.005, 1.6*^7}, {0.0062, 1.95*^7}, {0.0075, 
        2.3*^7}, {0.00875, 2.7*^7}, {0.01, 3.1*^7}}]}}, {
      {RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[6], 
       AbsoluteThickness[1.6], GeometricTransformationBox[InsetBox[
         StyleBox["\<\"\[FilledCircle]\"\>",
          StripOnInput->False,
          FontSize->8.96], {0., 0.}], {{{0., 0.}}, {{0.00129, 4.8*^6}}, {{
        0.0025, 8.333333333333333*^6}}, {{0.00375, 1.25*^7}}, {{0.005, 
        1.6*^7}}, {{0.0062, 1.95*^7}}, {{0.0075, 2.3*^7}}, {{0.00875, 
        2.7*^7}}, {{0.01, 3.1*^7}}}]}, {}}}, {}, {}, {}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{Automatic, Automatic},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0, 0.0109375}, {0, 4.*^7}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.7723576977806387`*^9, 3.7723578560728683`*^9, {3.7723578958910923`*^9, 
   3.772357900978519*^9}, 3.772357938593742*^9, 3.772357968863305*^9, {
   3.772358012245631*^9, 3.772358024032151*^9}, 3.772358057882593*^9, {
   3.7723580905076303`*^9, 3.772358095583256*^9}, {3.772358134993271*^9, 
   3.7723582098761597`*^9}, {3.772358242819211*^9, 3.7723583069937687`*^9}, {
   3.77235867695652*^9, 3.772358680929736*^9}, 3.7723587181256*^9, 
   3.7723879331837683`*^9},
 CellLabel->"Out[33]=",ExpressionUUID->"d674ad7c-f351-4869-88a3-11d70336bd78"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"h", "=", "0.001"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"ExperDataStaticExtInterp", "[", "h", "]"}], "-", 
    RowBox[{"ExperDataStaticExtInterp", "[", "0", "]"}]}], ")"}], "/", 
  "h"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"ExperDataMedDynExtInterp", "[", "h", "]"}], "-", 
    RowBox[{"ExperDataMedDynExtInterp", "[", "0", "]"}]}], ")"}], "/", 
  "h"}]}], "Input",
 CellChangeTimes->{{3.7723354645031567`*^9, 3.7723355328114147`*^9}, {
  3.772349946076236*^9, 3.772349971376644*^9}},
 CellLabel->"In[34]:=",ExpressionUUID->"3227967f-a37c-49ed-9175-9fe4d431cc39"],

Cell[BoxData["3.790400000000002`*^9"], "Output",
 CellChangeTimes->{{3.7723354652954493`*^9, 3.7723354898880873`*^9}, {
   3.7723355257188377`*^9, 3.772335534689468*^9}, 3.772347811221644*^9, {
   3.7723499583349037`*^9, 3.772349978985915*^9}, 3.772350014550988*^9, 
   3.772387934917963*^9},
 CellLabel->"Out[35]=",ExpressionUUID->"872028df-9a71-42ab-9276-6fad752dd670"],

Cell[BoxData["3.8308456086675034`*^9"], "Output",
 CellChangeTimes->{{3.7723354652954493`*^9, 3.7723354898880873`*^9}, {
   3.7723355257188377`*^9, 3.772335534689468*^9}, 3.772347811221644*^9, {
   3.7723499583349037`*^9, 3.772349978985915*^9}, 3.772350014550988*^9, 
   3.7723879349199533`*^9},
 CellLabel->"Out[36]=",ExpressionUUID->"05d1abf5-4dd3-4657-b7f5-a7f53c114a71"]
}, Open  ]],

Cell[CellGroupData[{

Cell["3rd order", "Subsubsection",
 CellChangeTimes->{{3.772358814633583*^9, 
  3.772358816386402*^9}},ExpressionUUID->"bc2aeae8-2015-4b44-b1f7-\
cb7995fb441f"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dmax", "=", "0.03"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"f", "[", "x_", "]"}], "=", 
  RowBox[{"Fit", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"d", "/", "1000"}], ",", 
        RowBox[{"Stress", "[", 
         RowBox[{
          RowBox[{"d", "/", "1000"}], ",", "1"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"d", ",", "0", ",", 
        RowBox[{"0.001", "*", "1000"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"{", "x", "}"}], ",", "x"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Stress", "[", 
      RowBox[{"d", ",", "1"}], "]"}], ",", 
     RowBox[{"f", "[", "d", "]"}], ",", 
     RowBox[{"ExperDataStaticExtInterp", "[", "d", "]"}], ",", 
     RowBox[{"ExperDataMedDynExtSmooth", "[", "d", "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ExperDataMedDynExtSmooth2", "[", "d", "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"d", ",", "0", ",", "dmax"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", " ", "All"}], ",", 
   RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", "\[IndentingNewLine]", 
   
   RowBox[{"PlotLegends", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<Linear\>\"", ",", "\"\<Landau3*a\>\"", ",", 
      "\"\<Experim. low rate\>\"", ",", "\"\<Experim. medium rate 2.81\>\"", 
      ",", "\"\<Experim. medium rate 0.292\>\""}], "}"}]}]}], "]"}]}], "Input",\

 CellChangeTimes->{{3.7722508310115633`*^9, 3.7722508483839207`*^9}, {
   3.772250974625247*^9, 3.7722510045016403`*^9}, {3.772252103191379*^9, 
   3.77225211703337*^9}, {3.772252166064659*^9, 3.7722521709281397`*^9}, {
   3.7722522979244137`*^9, 3.772252298064801*^9}, {3.772252356742598*^9, 
   3.772252374784587*^9}, {3.772252423821867*^9, 3.772252439376793*^9}, {
   3.772252567994115*^9, 3.772252576904614*^9}, {3.77225261508735*^9, 
   3.772252615233014*^9}, {3.772252754599139*^9, 3.77225275727354*^9}, {
   3.77225375149601*^9, 3.772253752476163*^9}, {3.772253847489312*^9, 
   3.772253880323216*^9}, {3.772254230832884*^9, 3.772254236654265*^9}, {
   3.772254275536106*^9, 3.772254276123672*^9}, {3.772254541495772*^9, 
   3.772254567598219*^9}, {3.772254630671686*^9, 3.772254694618869*^9}, {
   3.7722547363253*^9, 3.772254742982931*^9}, {3.772254784779194*^9, 
   3.7722548196197653`*^9}, {3.772255435938899*^9, 3.772255517057938*^9}, {
   3.77225559746093*^9, 3.772255615778727*^9}, {3.772255647443912*^9, 
   3.772255681917776*^9}, {3.772255715553001*^9, 3.772255744147625*^9}, {
   3.7722559203252487`*^9, 3.772256018044723*^9}, {3.77226925623964*^9, 
   3.7722692753880663`*^9}, 3.772269356828538*^9, {3.772271087317779*^9, 
   3.772271111674114*^9}, {3.772276592175831*^9, 3.7722766270617847`*^9}, {
   3.7723351401361628`*^9, 3.772335143014654*^9}, {3.7723353099310217`*^9, 
   3.772335340389203*^9}, {3.772335563333227*^9, 3.77233556484482*^9}, {
   3.772348133837274*^9, 3.7723481396412*^9}, {3.772348196668665*^9, 
   3.772348232621614*^9}, {3.772349995503333*^9, 3.772349995875544*^9}, {
   3.772350028352704*^9, 3.772350031136403*^9}, {3.7723508994917994`*^9, 
   3.7723509057391233`*^9}, {3.772358377805985*^9, 3.7723584715494347`*^9}, {
   3.772358695715138*^9, 3.772358733859083*^9}},
 CellLabel->"In[37]:=",ExpressionUUID->"b6e0fb2b-21b9-4aa2-acd5-b4391d914e73"],

Cell[BoxData[
 RowBox[{"3.6757830939487143`*^9", " ", "x"}]], "Output",
 CellChangeTimes->{
  3.77225561617383*^9, {3.772255648999432*^9, 3.772255652475399*^9}, 
   3.77225568310441*^9, {3.772255715906597*^9, 3.7722557443801603`*^9}, 
   3.772255924093693*^9, {3.772255968248919*^9, 3.7722560187460003`*^9}, 
   3.772257731604315*^9, 3.772257763550974*^9, 3.772258318993701*^9, 
   3.772269234850355*^9, {3.772269275994589*^9, 3.7722692846993933`*^9}, {
   3.772269345431892*^9, 3.7722693574170313`*^9}, {3.772271090174115*^9, 
   3.772271112886817*^9}, {3.772276601996491*^9, 3.772276628044529*^9}, 
   3.772335143733492*^9, {3.772335327461673*^9, 3.772335347311558*^9}, 
   3.772335565216633*^9, 3.772347812067973*^9, 3.772348142517158*^9, {
   3.772348206374415*^9, 3.772348234340475*^9}, {3.772349996410104*^9, 
   3.772350032747957*^9}, 3.772350368879985*^9, 3.772350906096163*^9, {
   3.772358384410532*^9, 3.772358400673436*^9}, {3.772358450345291*^9, 
   3.772358471947193*^9}, {3.772358696823798*^9, 3.772358734104312*^9}, 
   3.7723879366675873`*^9},
 CellLabel->"Out[38]=",ExpressionUUID->"9b868566-770b-4b0c-80ed-f569c944dd08"],

Cell[BoxData[
 TemplateBox[{GraphicsBox[{{{{}, {}, 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.368417, 0.506779, 0.709798], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVxXk41HkcAGCbkWRzr5xpsq1KGi1jqc3HJx0YKROReDqFLZoeIoqZ3JOa
GcYcvtWPDvvkWlGOxDTl2Mpao9RjSdgOhDGrtVtPq3b3j/d56QeOsSPnaWlp
bfvP/6+4eXImXE/Hs2M6dOGjj194sWL8b+r5v/T0c9IJsucOe4VlFpQMsd55
TtKWVC/karxG8pmyRyxtYF5zDDHkamGpRY6JkmUKQzOTxIVrhLrPGaomlj1U
aNFz+Nyl2LiBxmhkuUByHatlPc8ZRXFxL+pZ3nD5twqdK2e80A9XHUm0DgLe
nIqZ5haI2d83lZcMRYDVytEkE9lu3JXXP7739kFgNoxzlo7tQ1HwgLz7bTTY
7O7Fw3sjMUpePp7dGgspgkt2TXkx6HIxbZVh+nHID4uY0pjE4eTVdCsL8wQw
5JvRPLyP4+xo7EBbaiLcbOydCAiJR1E/KdtvmQzuS8wrZuQn8JRetc+V3FPQ
k/jwfs+tJLTT2C4rmUiFYLOA2csdyVixMYq2FXlQHzjgbDn/NB5zdBkWJp6B
CwrxJGWXhqcniu+r2tPhQMAHU44LD6N0123UGGSC8uoDd132GUy+t7V8BLPA
9J5bvzAiHZ+y/7LoFGRDCu31Ay1eBoamW3Cy23Og/k3AvRBhJi7nD3R5mPPB
0WRNTHFdFnZPdD7o3HcWarcr87ueZWPJJKOSIcmDUpsWVderHLzT5irNe3kO
DtXEjnXo8XFuXVX0bWsBbKXH5+VYn0XbRTv72rcJgR1RWWi5Ig+9B/wdXvuL
gHe4w6yCdQ71P5ftUo2KwGdnATWz5zyeFITwG4T5YJrWp/bMEGBv54WFAtcC
GO6N9ggtFaJ3nyVv268FMJ0VlqQfJsKQHWPvF5wQQ82mV0HTkyLkyjSRVUaF
UNj5KZmemI/Sc0I6lhVC5KFwz5p/8pE9lTasCpCAut/tdYegAOm218W+4xJw
zjVvsV8sRkkj3U8hkoLLNOToVojR1v7pn9ZOMjDm2BWsYRaieL2dOrpDBt6B
TJ2q5kL86e7ttzVH5SBw3ZBRu12CgyaCgZe6RRDbXbVytk+CNT213cY3ikAa
Z/RCESvFeM84E0YwgUw6N/X5Zyn6phz0ZU4TOKimtazNlWEbvSzjO+kFeGzy
/OcvbeTIlq+oXfbtRfCRxR+x+VGO3c2ejLlHF+EWpR5LdirCP/xLq59wLkHl
tWK3PcoidKjc/LVYnwKz5h9UR30Jusvc7bYYUHDK2HBIxSLok+5o9cGIgqad
vVOuAQSjQ42N9ppToHRT6s+xCZbTBj+uolPQX2zqez6coFNEwpP7bhTw14c9
rDpG0MXgavrMfgqkX2kG1RKCmz5IUksPUbBFqK1hywkGvco9GRpFwWjt6nkN
hGBCU1yc4igF2kLFN9xigvWH14WdTaJg6s4zjmEZQY+7PWuXnadg+Vv+Audm
gn7lbaufCil4XNtlXaggGCZpcMgtoKCVR2e8VxJMOXLJVi2jYHHr38HKdoLN
i2P0mq5QoPgUdG1HN8Ff5oXTYkspWGpo2Hirh+DgVMBnu+sU1AU+67ToJfip
1XU2q4oCsXbeu+E+ggbVDhr3GxTMVyTpbh4guIRYTUzUUjDixbEuGyS4JmvR
G6qOgutrExiLhgkCR2sksJECNTPD+/jvBP8FsKdCqQ==
          "]]}, 
        Annotation[#, "Charting`Private`Tag$2753#1"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.880722, 0.611041, 0.142051], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVxXk41HkcAOB5rB5rPSYpuRbRipQktEg+PkliaDfJ9aSmsM5RUSGrMLlG
zQzm/Jm+OdY+TOnyNM6imnpyFB1a26xiJVdqUrQ2bbv7x/u8lgcOBsVo0Gi0
wP/8v21D+swe7UWeeU70vZ9oGl6M+IAG7YARTwMXRayKPeQVcaq04gXjved8
entMI1vtNVziIu5ifAWTJsroGjYNa4wK9DsYS2EkcVu0jK2HWn849LUwVgJT
oPVSzF6BTZs1HZoYTkBsx6uF7PXIT05+rmB4w/WUAbaA7YX+aJd4zDQYJG5z
z5guOzHfo0Ve8SISBKsymKYl4RhS/GxiX3MU0PxC45KHmcjfrZL0TsaBrFdT
4LIrBmMl8on82yyQ9ry5Rs+IRyfZCbvFuYchi0pkPPjIwtfVuSZGy49A9+MV
5++tP4yzYyyVMusYcAzkaxogFfnPqLr9xhnQ2r4wqpd5FDO1L22vKswE8xpW
eRo/DS3UZlYVU1lgaleswyzPwPNbYjV9MRs4hinXWh9l4sE1TkO8Yzlg75UW
seJlFv48de5W351cqByv+nZBLxtjtdy3qOmnYD6bOVFrl4MZN33lw5gHkxvr
nSIcc7E/aM6om5sPBpXPTXx2sDEs1+hQ/p0C+M323EnnsFNoXaS677a8CJ5E
b82JTM/D3qnue91MDhA1MH3L8rHitcMFB2Ex0KiuNpfyAmxVOouKR06D7KOb
IrK5ED+718c1m3KhL8XVnqMsQjPdXQN3AnnwaPOZnLQuDnqrAmxGA/hAv343
hbwqRp0vdSF9Y3z4M5w33Pn2NKZzQ4saeSXwVu4xvW4xF590l3/DdS6FbobX
0+DVPPQeMM4OfFAKFnnxBIz4GPrj+F9fHy2Do57WFe9+5eNJsTqmXk8Akvza
Ieu1JSg6zbPEOgFUzT3ltVwswaDpE0N9O4TQ2ZGmm+RRipZmtWV+E0JoMx7q
PNdRisImS/8bfBH09LQEq38oQ7OV/R9M7cVQS480v91fhmWbLN7E3RVDQro8
ZC5MgBfbmyevJEnAzyMhctuYAAf1uaoRLSnMct9VGSUJ8crDq71LLkvhxqYp
hdaCEFM9k/UddlPA0uO7hRaJ0O94lJ/LWwrGUxf+ttUVo9Kyjv29qBxyHFUc
Q5kYgyS2V602yGA0Kye80FKCvW2eDp+7ZLB+5pEo/BcJvguoufT40FmQdIiq
fDdI0eaCz3dlOgTum1lNNyqk6Cp2tdhGJ9B9XuSv1SzF7blrTOb1CNxz06kN
aZViXNgSvX3LCdzc/eHAbLsU5ZqDn+wsCSjO3B3Y0ClF+8gjj29tJEA+Jyjr
VVJ0olfnzuwncHCwobyaRuHWeWFWTTQBVqLt/IwGhcEvC9PDYgkkzMtCtiyi
8EhLcvKNJAJRBvlLhrUpVPzkHsFJIxAWGFZgvpRCt/aHjlZnCHhdX0iRWlPo
L1eu7ecR8GQc7hu3oTBC2GhTWEpg0++j9q52FB5PPGv2RkzAefbB+NN1FLYZ
xmu3VBFYbV+5d5krhT0aezRZNQRWtRq0RblTODi944tFLYGVfhzjBg8K/7nt
PJtXT8A8JvXJTqSQfslG7XqZgMn7McdKbwrNKZOpqasEDHP28NQ+FK7L031F
rhFYtvjha9hOIRyiDe9sIqB/1sef50/hv9LnQ80=
          "]]}, 
        Annotation[#, "Charting`Private`Tag$2753#2"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.560181, 0.691569, 0.194885], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwV0Hs0lAkYBnBJK4mkm7STlCKSU0wprdcbJUaUXIaaYiMKI/coGYNhvk8I
Y1xmvqRVKCkSqaMbXZRcyi5NhaUUxdBlS4v99o/3POf3z3Oe8+r+Huzsq6ig
oLCfvv/ToOrY2D6VGZYKCjyddIvpVqzDDlUqDn2WChzlC/KiXivPxMzCbtZn
S576sKnT2VGr3tNMcRNrOliVu6dvp6ZhsVay5h3WPOgZ37LXJXcuKr8yaa1j
rYBlUYWOyam6WPubkkktyxQMtobpa0Wswwwu9811ljUE2v568Ol2RHs0DIhc
4gKzJzWwN383CrbUlRV2cyAywKLn0BkPdCNffjhw4yC85OjbC0VemOEqy20Z
9Ifwh3Yu5Xwf9Mst+yC4HwQlqq/U1yr7o6nkpOEcfghMP94+2nYlAD+e42tr
LQyHjTXW2ZyVwfh1IEjWEBsJnc3P6zjLQjDjZX6p9+JoYC9gzqq6E4rHVSp2
FKUchzCbxa5i83DUkTOWFw7FQnWBx01LKgIvbvVTskUeNIzwWgIXRWGwkWlP
emQ8tBeUNLpzj+GJoTP3Whv50POBH5faFY1+ypu3ytUTwSa8OOHlhuMYfde2
rBeT4G9uvXlX4gnscP6m9SRNANV1IT6Oo7HI5msdFTQmg07jzn/2QxyuFMqa
Ny0UQrL3nOGvu3jYMvTk0RMvAuS9jkLZOA8LP5pcMhGRcDYl/raDKB5vNpjl
kH2pEOgqcBrbyMeJzeX+N5akgXuG6MhYEx8Zans6G3emA/Oj/BfjPQloLXPQ
f+uQAUzHsxylNwmoOlXq1jqQAW6rzfj5nol4LM1dWJN+GrZMzuCEvErEF08K
ZqWZZUJajPnPd25JaN25mLfzWSZEWN9jvW5IQvdd77/PjMgCnWvF16s2CjBO
LPct18gGPZlP7g9KgDmp6bpYmg3fPAdi1ygmo/Onkz2tjiLQspl+81RQMuoy
SrLsPoigKIz5Zaw5GUW1uvb1GTlgZR1tOm99CjJWdHxZYiyGfb738izJFMyy
0Bn2fyCGPuMH435vUvDy7RuDVwNzwaNdAlJTIb7WTJP1KedBp31lr8VJIV5t
q2yZeyUPJrpsmfzHQgyz5GqauObDNPP4vAQ1Au1iDtoxR/LBw8G3IcqWwAbd
0oSNOQUgaMtSX8Qj0DnXoHL5egmM9Q+K314msOWWpclEkwRsXPT4A38SqOFu
JcxulsD6fv+p5L8I3DWKfUatEkiM2a1q0Elg+8ptYs8OCdzitjw71EXgizTW
VE23BKx97it2ywjs9GK3hH6RgNJ4pUZdD4E9SqHB7xlSeMP/6Ko9SOCoQ3HF
86NSOEGdn8P7SSD3WVS5ZZgU/vj8c3ic9qCT/cXSCCnYjBj5R/xLYL/zyPm4
GCnEB77XPzxB97M3UUaJUtCJD5Cxpgis92lOTRBLYamep/KkIompJ74eNq2X
wi+yLq0Fs0hUmXzoR92RQo1CgXYCbUFcvq/KfSlQ3gXX5LR5fPDufiiFUr9V
V5tUSQxLEbLJNil0CSaaotVI9Mxm2Pb1S0G2/GhkgwaJ+pe26WWpUmC9189Z
dRGJ5mJzne3qFDzXiOdH0d7BN9L+oUEBx8jRpI+2P3uuxoGFFBTNuDlZq0Vi
mdLrn4a6FFQRg6H7tUk05oQ/v7eBgg3V44UZDBJN1c/xx7wpaI+9nHN6BYk2
P0SxxT4UWGinbf5E26U/5Rjbj4KS5n72Dj0Sw+u43PpAClpGKpQmaF8/tNmT
iKJgPufRQq9VJG663bZu+SkK7t7Yu0BtNYn2ZQ1rOtIpmFlkz/ai7Smq0U/J
pMA9JnN+Fe2YACljWEyBU8jb726GJN5adFilroiC4GvnuLlGJD5V3KcUVEzB
HbZvzADt158cp3RKKJB0ySeZa0icvG/2NamcAuYAI7SNtnqFvtz8CgWaFQVe
S41JXJqvPTRUScGtxwUPj9Bem6T2jqqm97nMvnCdNhxV6N1dS/8nSDY5bS2J
/wEUzMJ7
          "]]}, Annotation[#, "Charting`Private`Tag$2753#3"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.922526, 0.385626, 0.209179], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVjmk8lAkcgCXWWhnXVpimybUsRTkiNP/+OcJgJUXzY1WaRUop95XcVzPj
mHlf4l2rJFc52kLlCr925Gir3SK5OkQHlS2tY9sPz+/58nx4NA4dd+dKSkhI
uHzlf+s1RL73lpVmJSU1GJoVSu5gBzo3yDpPsCTC/9xlqzi2g5OcWzLC/sBK
COKK5xVmd4zlmBFi9ko4N65PFSmswDLVNOU2tgqMiaQ3LdGUUOaJ0UAzWwtC
xg+x+uU1sHG7lFEj2wRWB81yIuW2oCA4+Ok1tg2sD/us2PcNohPqB4XTPWDU
VxxBNOzGVOvmypIRH6huk+hdyeLgvqzBV75NftAXoBis5XQQBXuHyP6pAOCx
+G2pYi76k5WvUm8fg8zyOVmt+UA0KYrXV0gMgcXuI91KscH4+nyiuuqaUDju
LmPpcDIE514eG+qMC4erF/o+mcafQsFgYcVBtShQkFFRvNYQhjGyVxxK02Ng
VrR0xLA9ApkzDM2S6TgIVX8eJtkRhVU7/aV2YQLQq7orH0/H4HEDk1F++BkQ
+5mzzN7GYez0rx0DXYnw8Xa8O29VAvrLWO6coSXDsz/MVQroZzCqfVflGKYA
323GznBtIj50/0e1h5cKDB173dYfktArUfVEalcaMA3UXJc1k1EnY6h325oM
eGO71WTZOAX7p3vu9BzIBKLcQ7veOhVLXhtVGwmzgJZrNeNtkoY3Ok1FWRPZ
8Fk7eip2ezouWtYENNF5ADoXbjWZZiBDfs+jLhc+JG9bXfBBJxNthpx1nzsL
YOGy0dR+/SyUW67YN/BSADmqV9u11mVjJM8z4zo/B+hPTRyfbjiLD3rOfccz
zYXORZQ6zeShzSO1BJe+XOhIWyjvkeOjp9vk52/D8oA2SnfzlRbgaWKGW6OY
DwsZp4zbYwQoyuZrYEU+6ChzlK3eCdD9TfzogKsQmqq5gxNBOajBuJTn+EoI
B/Re2JEjOShs1HBqEYiA5Spz18c7FxlaDz/SNxFgfOd+LA7kYp4V821ANwH8
2NhVbMzDy61NU3VHSSj93sIJb+ThsDJvaEKmAKbStdSiNuZj3b36fqXaAlhr
y7h4+mI+nmIFKxvtLYTunbzSEoYQHaP9HM3eFUJ68eoQ3WwhdmpUJJmLzkGs
3fVWz5UidCf16jWNi8BDM5IjcVKE/TdZRoviItj9sd1rcliEs85lV+6fKIZh
gwfJe90J1K22086To8C7Vtrp5xYCLQgLpj2NgnamiANtBDokGqjPK1KghNG+
GzoIDPBSUvRdQ0F0wXb78S4CK6WG/9XXoODwb7Vp3F4CN/mE3u/YSkGhvMNg
4BMCTWjnE98fpCDP07ol8AuBtvPCuLLDFPClbPY7LhDo8Sw90sufggi/4km9
JQJDm4ODW45SQDN3efZyBYnXfrHkZEZQsM+s4CeuLInbWu9t0TxLgR1vfM5H
lUSnys6ND/lfe472X9bqJHKE13XTcyno+Lupkr6OxOigYsZbgoKZyZHNj5kk
3lwbKNtcSgH3kGrdHl0S70p6Sx0r+/p365OJ8Y8kDr9xXWZeoqBGrF+laEDi
0m3TuZQaCsq/ZIX3GpJIu6I7Y1FLQRqjWVy1mcT1herT0/UU7HlvpZJpTKJh
ivwL6ncKVkUpuQWYkggnJMZ2N1LQUGd5xn4rif8B8hAvSA==
          "]]}, 
        Annotation[#, "Charting`Private`Tag$2753#4"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.528488, 0.470624, 0.701351], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVkGs0lAkYgGdlV2qz0lkpZ7JuKVZT6YKY15siBsnmsk4lKpQxxv1Sg51B
4zYz3zAzoj5iCbl3cqskl1jaKNkjk9vRas1E0ratU7Ltj+c8v59HPyDM45wK
hUJx/cL/3nYrbvGE+tf01eFdQa1hFDvGeZdb6i7TdOkjdsvbreN2viniognG
O/rB4Wdxmlvn7aaIvbI+xirY7flRUmf02a5U57JWO2MDnMmK1TMw0EC1F7TB
VoYhXJCqjVKpW7DZVpXWzLAARty7Yv0NO1DEYo03MuzhvUH9KdqKLTqjaUiM
7nEY3Bf2Ruf7o5hm01pZNHESIrM6jZo+eKFX5uisX8sZ2LgdfB+YnUKRpzxv
QBEM9/mR2W8XAzAor3I2rTMUZCMZQ7m8QLS4mmj6HTcc3mlbKew9L+DrEu5m
He0oSOI7WBaqs/D9q1B5FycGYt1NIz61sFE0ml/hvykeFOF17WxBBF5Urz1S
zL8Iaydmrbj0KNRboBoUKTlg9YflQ+OpaLx5MEjVEZMhbf7ogdHrsRhmZjEp
jPkF2NTgup8hHi8pCzsGu7nA6bPsUywlYJCa9cEFjRRoVA2IqSq7hPEPHCun
MBXGwfsQa1ciDnv8o9MvSAO/nmoT1ZEk9OHqsNO6L8PEsSl5QU0yGqfLf7fS
Tgdbh09xhYm/4ICyv7f/dAaUE0/tPjhyseg1rYomyQTKc8EPVl/x8E7XHmnm
dBbUL+n4HO7h4bJ1dXCLrgCKbax7MzgpSF3300i3qxDsDeN/7TBPRXu5i8mf
LiLwnHNlbFek4tqVCq/BVyJ4ajvic1qWhnEC7/QmIQF6biIvB/fL+Ky/YI1g
jxhuyxMVcyp8tB/ZlOz6WAzi6aQWyg0+erv/9e/q6BzgSstCDL3SMUm2cK5a
MxdmZDuR8iEdpVlCfazIhdk2McdGnIEec4mTg24SsLSgJs/tz0R9anmO06wE
Ei+tBL4ZyERJs75zm0gKPcrnuTPhWUg1HP5b11wGAopEwfkmG3MO6M0HP5QB
v8BmdZ0kG2vutyjqmXnQq3DMP7BDgGNaAvm02hU42/Fb9+cWAdY/aRhYX3cF
OCMzHgVuQoyks7RonvnALku+e31MiE4JZ5z2vsmHsqaGENudIuzSr+DtlxbA
msh2Sk+0CD3ytjUY7L4KWkq1KP/bIhy4S6ct912FYNvGyaRFEb51Ka0dYl+D
l0tPS0ItCDSpOmyUs5aEomi/niYWgZYySz0HDRKsloVV98IIPMI127ykSYIR
rzeuk01gsM96TT9tEpTpxx4+jiCwUnXso6k+Cc68zr0vYwg0Pxk11LGPBLqx
SeCaJAItNEq4i/4kqGyd+/aQkMBDSxJO6VkS2mOpQw4iAo+/5Mf5BJEw0OOZ
7EQQGNXKYrUxSZgKGBe65RDYGGjtmxFLwlDKUZ63jECr+092GWSTMB9bf+ds
IYHOlV0/DgtJaKvxXzxXRKCvpMmELyZhcnqTZtB1AhNCrlHnZSRoOxVTzpcQ
eHfjefXWYhI8VrQaQ24Q+EjlhGpoKQlV21YZM8sJHJtzW9ErJ8HL/WMCs4LA
z5173qdWf/mTv+oF8yaBGrUmC5Z1X/rubXjFrCJwS/5mpbKBBPNxs+fMagJ3
pK6bIW+TwFpmNDBrCAQ2ZepYMwkU3agIZi2B/wGdgUaw
          "]]}, 
        Annotation[#, "Charting`Private`Tag$2753#5"]& ]}}, {}, {}}, {
    DisplayFunction -> Identity, Ticks -> {Automatic, Automatic}, 
     AxesOrigin -> {0, 0}, FrameTicks -> {{Automatic, 
        Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
        Charting`ScaledFrameTicks[{Identity, Identity}]}}, 
     GridLines -> {Automatic, Automatic}, DisplayFunction -> Identity, 
     PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.05], 
        Scaled[0.05]}}, PlotRangeClipping -> True, ImagePadding -> All, 
     DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> 
     Identity, Frame -> {{False, False}, {False, False}}, 
     FrameLabel -> {{None, None}, {None, None}}, 
     FrameTicks -> {{Automatic, Automatic}, {Automatic, Automatic}}, 
     GridLines -> {Automatic, Automatic}, GridLinesStyle -> Directive[
       GrayLevel[0.5, 0.4]], 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None, 
       "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
           (Identity[#]& )[
            Part[#, 1]], 
           (Identity[#]& )[
            Part[#, 2]]}& ), "CopiedValueFunction" -> ({
           (Identity[#]& )[
            Part[#, 1]], 
           (Identity[#]& )[
            Part[#, 2]]}& )}}, PlotRange -> {All, All}, PlotRangeClipping -> 
     True, PlotRangePadding -> {{Automatic, Automatic}, {
       Automatic, Automatic}}, Ticks -> {Automatic, Automatic}}],FormBox[
    FormBox[
     TemplateBox[{
      "\"Linear\"", "\"Landau3*a\"", "\"Experim. low rate\"", 
       "\"Experim. medium rate 2.81\"", "\"Experim. medium rate 0.292\""}, 
      "LineLegend", DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #3}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.922526, 0.385626, 0.209179], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.922526, 0.385626, 0.209179], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #4}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.528488, 0.470624, 0.701351], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.528488, 0.470624, 0.701351], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #5}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.24561133333333335`, 0.3378526666666667, 
                    0.4731986666666667], FrameTicks -> None, PlotRangePadding -> 
                    None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.368417`", ",", "0.506779`", ",", "0.709798`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.368417, 0.506779, 0.709798]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.368417, 0.506779, 0.709798], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.587148, 0.40736066666666665`, 0.09470066666666668], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.880722`", ",", "0.611041`", ",", "0.142051`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.880722, 0.611041, 0.142051]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.880722, 0.611041, 0.142051], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.37345400000000006`, 0.461046, 0.12992333333333334`], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.560181`", ",", "0.691569`", ",", "0.194885`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.560181, 0.691569, 0.194885]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.560181, 0.691569, 0.194885], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.922526, 0.385626, 0.209179], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.6150173333333333, 0.25708400000000003`, 
                    0.13945266666666667`], FrameTicks -> None, 
                    PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.922526`", ",", "0.385626`", ",", "0.209179`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.922526, 0.385626, 0.209179]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.922526, 0.385626, 0.209179], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.528488, 0.470624, 0.701351], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.3523253333333333, 0.3137493333333333, 
                    0.46756733333333333`], FrameTicks -> None, 
                    PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.528488`", ",", "0.470624`", ",", "0.701351`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.528488, 0.470624, 0.701351]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.528488, 0.470624, 0.701351], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}]}],
              "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2, ",", #3, ",", #4, ",", #5}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", "None"}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{
  3.77225561617383*^9, {3.772255648999432*^9, 3.772255652475399*^9}, 
   3.77225568310441*^9, {3.772255715906597*^9, 3.7722557443801603`*^9}, 
   3.772255924093693*^9, {3.772255968248919*^9, 3.7722560187460003`*^9}, 
   3.772257731604315*^9, 3.772257763550974*^9, 3.772258318993701*^9, 
   3.772269234850355*^9, {3.772269275994589*^9, 3.7722692846993933`*^9}, {
   3.772269345431892*^9, 3.7722693574170313`*^9}, {3.772271090174115*^9, 
   3.772271112886817*^9}, {3.772276601996491*^9, 3.772276628044529*^9}, 
   3.772335143733492*^9, {3.772335327461673*^9, 3.772335347311558*^9}, 
   3.772335565216633*^9, 3.772347812067973*^9, 3.772348142517158*^9, {
   3.772348206374415*^9, 3.772348234340475*^9}, {3.772349996410104*^9, 
   3.772350032747957*^9}, 3.772350368879985*^9, 3.772350906096163*^9, {
   3.772358384410532*^9, 3.772358400673436*^9}, {3.772358450345291*^9, 
   3.772358471947193*^9}, {3.772358696823798*^9, 3.772358734104312*^9}, 
   3.7723879368835773`*^9},
 CellLabel->"Out[39]=",ExpressionUUID->"1719ff9f-71ab-43a8-a205-d669b47b172b"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dmax", "=", "0.03"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"f", "[", "d", "]"}], ",", 
       RowBox[{"Stress", "[", 
        RowBox[{"d", ",", "a"}], "]"}], ",", 
       RowBox[{"ExperDataStaticExtInterp", "[", "d", "]"}], ",", 
       RowBox[{"ExperDataMedDynExtSmooth", "[", "d", "]"}], ",", 
       RowBox[{"ExperDataMedDynExtSmooth2", "[", "d", "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d", ",", "0", ",", "dmax"}], "}"}], ",", 
     RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Linear\>\"", ",", "\"\<Landau3*a\>\"", ",", 
        "\"\<Experim. low rate\>\"", ",", "\"\<Experim. medium rate 2.81\>\"",
         ",", "\"\<Experim. medium rate 0.292\>\""}], "}"}]}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"a", ",", "1", ",", "10", ",", "0.5"}], "}"}]}], "]"}]}], "Input",\

 CellChangeTimes->{{3.772258362208325*^9, 3.772258403530744*^9}, {
   3.772259320442834*^9, 3.7722593490827713`*^9}, {3.7722593915113907`*^9, 
   3.7722594035006533`*^9}, {3.772259441649233*^9, 3.772259473704123*^9}, {
   3.772259604056508*^9, 3.772259683275116*^9}, {3.772269395439508*^9, 
   3.772269409168303*^9}, {3.77226944871436*^9, 3.772269448835977*^9}, {
   3.772269511462352*^9, 3.7722695462823877`*^9}, {3.77227112177039*^9, 
   3.772271121883951*^9}, {3.7722711607953367`*^9, 3.772271177035329*^9}, 
   3.772271228797082*^9, 3.772276649581072*^9, 3.772276707017564*^9, {
   3.772335666337817*^9, 3.772335691055481*^9}, {3.7723357575664473`*^9, 
   3.772335757832345*^9}, {3.772347824949437*^9, 3.7723478355645723`*^9}, {
   3.7723481554768753`*^9, 3.772348157582074*^9}, {3.7723482399541483`*^9, 
   3.772348240086149*^9}, {3.772348270658585*^9, 3.772348335494632*^9}, {
   3.77234847120266*^9, 3.7723484726689157`*^9}, {3.7723503840311413`*^9, 
   3.7723503922735233`*^9}, {3.772350493664414*^9, 3.7723504939596043`*^9}, {
   3.772350594743127*^9, 3.772350602923962*^9}, {3.772350918783649*^9, 
   3.7723509227382927`*^9}, {3.7723576573175173`*^9, 
   3.7723576832168407`*^9}, {3.7723584848983927`*^9, 3.772358539922626*^9}, {
   3.77235876534759*^9, 3.772358767453516*^9}, {3.77235883525349*^9, 
   3.7723588353841743`*^9}},
 CellLabel->"In[40]:=",ExpressionUUID->"5beb9aa0-77e8-4bb4-a3b6-5f3301c24a15"],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`a$$ = 1, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`a$$], 1, 10, 0.5}}, Typeset`size$$ = {
    732., {125., 131.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`a$3011$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, "Variables" :> {$CellContext`a$$ = 1}, 
      "ControllerVariables" :> {
        Hold[$CellContext`a$$, $CellContext`a$3011$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> Plot[{
         $CellContext`f[$CellContext`d], 
         $CellContext`Stress[$CellContext`d, $CellContext`a$$], 
         $CellContext`ExperDataStaticExtInterp[$CellContext`d], 
         $CellContext`ExperDataMedDynExtSmooth[$CellContext`d], 
         $CellContext`ExperDataMedDynExtSmooth2[$CellContext`d]}, \
{$CellContext`d, 0, $CellContext`dmax}, GridLines -> Automatic, 
        PlotLegends -> {
         "Linear", "Landau3*a", "Experim. low rate", 
          "Experim. medium rate 2.81", "Experim. medium rate 0.292"}], 
      "Specifications" :> {{$CellContext`a$$, 1, 10, 0.5}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{789., {177., 185.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.77225840455881*^9, 3.7722584233880453`*^9}, 
   3.7722593496815577`*^9, {3.772259398112463*^9, 3.772259403986608*^9}, {
   3.772259453171679*^9, 3.772259481481946*^9}, {3.772259623347006*^9, 
   3.77225968559298*^9}, {3.772269385883868*^9, 3.772269409606758*^9}, 
   3.772269449246793*^9, {3.772269495848444*^9, 3.772269546565522*^9}, {
   3.772271118958425*^9, 3.7722711871455317`*^9}, 3.772271229317424*^9, {
   3.772276707773242*^9, 3.772276717736383*^9}, {3.77233569180835*^9, 
   3.7723357041657047`*^9}, {3.772335758268011*^9, 3.772335766668838*^9}, {
   3.772347818264567*^9, 3.772347836411436*^9}, 3.7723481098801403`*^9, {
   3.7723481586516438`*^9, 3.772348166423306*^9}, {3.772348209417955*^9, 
   3.7723482441297007`*^9}, {3.772348337856411*^9, 3.7723483432878437`*^9}, {
   3.772348473064625*^9, 3.7723484813836317`*^9}, {3.772350385102449*^9, 
   3.7723503957951117`*^9}, {3.772350494436432*^9, 3.772350496415592*^9}, {
   3.772350604192789*^9, 3.772350613139923*^9}, {3.772350915194627*^9, 
   3.77235092808322*^9}, {3.772357664161953*^9, 3.772357684454875*^9}, {
   3.772358500630064*^9, 3.772358550305295*^9}, {3.772358768342897*^9, 
   3.772358771034783*^9}, {3.772358835725505*^9, 3.772358838080752*^9}, 
   3.7723879403483133`*^9},
 CellLabel->"Out[41]=",ExpressionUUID->"727a5820-5a2d-4908-bffa-ada58eb8ae88"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["4th order", "Subsubsection",
 CellChangeTimes->{{3.772358820709601*^9, 
  3.7723588225068893`*^9}},ExpressionUUID->"a9047f7e-c4b0-42c2-9c13-\
8c3c44722a38"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"dmax", "=", "0.03"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"d", "/", "1000"}], ",", 
       RowBox[{"ExperDataMedDynExtSmooth2", "[", 
        RowBox[{"d", "/", "1000"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d", ",", "0", ",", 
       RowBox[{"dmax", "*", "1000"}]}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.772359142891779*^9, 3.772359230715601*^9}},
 CellLabel->"In[42]:=",ExpressionUUID->"7cfd2991-148b-4576-b2aa-5eedb925db1c"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fit", "=", 
  RowBox[{"FindFit", "[", 
   RowBox[{"data", ",", 
    RowBox[{"Stress", "[", 
     RowBox[{"d", ",", "a", ",", "EL", ",", "FL", ",", "GL", ",", "HL"}], 
     "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"a", ",", "EL", ",", "FL", ",", "GL", ",", "HL"}], "}"}], ",", 
    "d"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7723589317057657`*^9, 3.772358957632*^9}, {
   3.772359135510191*^9, 3.772359140131393*^9}, 3.7723592257596817`*^9},
 CellLabel->"In[44]:=",ExpressionUUID->"9821070a-224a-4597-957b-98b755efb9cc"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"a", "\[Rule]", "13.629953384447226`"}], ",", 
   RowBox[{"EL", "\[Rule]", 
    RowBox[{"-", "5.737089624261016`*^15"}]}], ",", 
   RowBox[{"FL", "\[Rule]", "7.982998685904592`*^15"}], ",", 
   RowBox[{"GL", "\[Rule]", "7.198289816347624`*^14"}], ",", 
   RowBox[{"HL", "\[Rule]", 
    RowBox[{"-", "8.686543720289861`*^15"}]}]}], "}"}]], "Output",
 CellChangeTimes->{3.7723589597670937`*^9, 3.772359232263002*^9, 
  3.772387947889003*^9},
 CellLabel->"Out[44]=",ExpressionUUID->"e85a73d4-0965-49b3-a358-822c65a35ea0"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fit", "=", 
  RowBox[{"FindFit", "[", 
   RowBox[{"data", ",", 
    RowBox[{"Stress", "[", 
     RowBox[{"d", ",", "4", ",", "EL", ",", "FL", ",", "GL", ",", "HL"}], 
     "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"EL", ",", "FL", ",", "GL", ",", "HL"}], "}"}], ",", "d"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.772358973508658*^9, 3.77235897661061*^9}, 
   3.772359226669353*^9},
 CellLabel->
  "In[396]:=",ExpressionUUID->"50bd4332-2385-43af-af52-7715da6e6b9d"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"EL", "\[Rule]", "8.133486894830394`*^15"}], ",", 
   RowBox[{"FL", "\[Rule]", 
    RowBox[{"-", "1.1366203039813472`*^16"}]}], ",", 
   RowBox[{"GL", "\[Rule]", 
    RowBox[{"-", "1.0182884423600852`*^15"}]}], ",", 
   RowBox[{"HL", "\[Rule]", "1.4227609644541294`*^16"}]}], "}"}]], "Output",
 CellChangeTimes->{3.772358979029627*^9, 3.772359233154496*^9},
 CellLabel->
  "Out[396]=",ExpressionUUID->"d9daccc6-b46e-4cf7-b024-d413425bc0f1"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"cons", "=", 
   SuperscriptBox["10", "12"]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"fit", "=", 
  RowBox[{"FindFit", "[", 
   RowBox[{"data", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Stress", "[", 
       RowBox[{"d", ",", "4", ",", "EL", ",", "FL", ",", "GL", ",", "HL"}], 
       "]"}], ",", 
      RowBox[{
       RowBox[{"Abs", "[", "EL", "]"}], "<", "cons"}], ",", 
      RowBox[{
       RowBox[{"Abs", "[", "FL", "]"}], "<", "cons"}], ",", 
      RowBox[{
       RowBox[{"Abs", "[", "GL", "]"}], "<", "cons"}], ",", 
      RowBox[{
       RowBox[{"Abs", "[", "HL", "]"}], "<", "cons"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"EL", ",", "FL", ",", "GL", ",", "HL"}], "}"}], ",", "d"}], 
   "]"}]}]}], "Input",
 CellChangeTimes->{{3.772359825138516*^9, 3.772359891092095*^9}},
 CellLabel->
  "In[421]:=",ExpressionUUID->"1bd97a10-0359-4014-a19e-427f5e9e34a9"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"EL", "\[Rule]", "1.`*^12"}], ",", 
   RowBox[{"FL", "\[Rule]", "1.`*^12"}], ",", 
   RowBox[{"GL", "\[Rule]", 
    RowBox[{"-", "2.2584896459158383`*^11"}]}], ",", 
   RowBox[{"HL", "\[Rule]", "1.`*^12"}]}], "}"}]], "Output",
 CellChangeTimes->{3.772359854946554*^9, 3.77235989601624*^9},
 CellLabel->
  "Out[422]=",ExpressionUUID->"700709f3-e9c9-411e-b595-b7790eb3dce8"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dmax", "=", "0.03"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"f", "[", "d", "]"}], ",", 
     RowBox[{
      RowBox[{"Stress", "[", 
       RowBox[{"d", ",", "a", ",", "EL", ",", "FL", ",", "GL", ",", "HL"}], 
       "]"}], "/.", "fit"}], ",", 
     RowBox[{"ExperDataStaticExtInterp", "[", "d", "]"}], ",", 
     RowBox[{"ExperDataMedDynExtSmooth", "[", "d", "]"}], ",", 
     RowBox[{"ExperDataMedDynExtSmooth2", "[", "d", "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"d", ",", "0", ",", "dmax"}], "}"}], ",", 
   RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<Linear\>\"", ",", "\"\<Landau4\>\"", ",", 
      "\"\<Experim. low rate\>\"", ",", "\"\<Experim. medium rate 2.81\>\"", 
      ",", "\"\<Experim. medium rate 0.292\>\""}], "}"}]}]}], "]"}]}], "Input",\

 CellChangeTimes->{{3.772258362208325*^9, 3.772258403530744*^9}, {
   3.772259320442834*^9, 3.7722593490827713`*^9}, {3.7722593915113907`*^9, 
   3.7722594035006533`*^9}, {3.772259441649233*^9, 3.772259473704123*^9}, {
   3.772259604056508*^9, 3.772259683275116*^9}, {3.772269395439508*^9, 
   3.772269409168303*^9}, {3.77226944871436*^9, 3.772269448835977*^9}, {
   3.772269511462352*^9, 3.7722695462823877`*^9}, {3.77227112177039*^9, 
   3.772271121883951*^9}, {3.7722711607953367`*^9, 3.772271177035329*^9}, 
   3.772271228797082*^9, 3.772276649581072*^9, 3.772276707017564*^9, {
   3.772335666337817*^9, 3.772335691055481*^9}, {3.7723357575664473`*^9, 
   3.772335757832345*^9}, {3.772347824949437*^9, 3.7723478355645723`*^9}, {
   3.7723481554768753`*^9, 3.772348157582074*^9}, {3.7723482399541483`*^9, 
   3.772348240086149*^9}, {3.772348270658585*^9, 3.772348335494632*^9}, {
   3.77234847120266*^9, 3.7723484726689157`*^9}, {3.7723503840311413`*^9, 
   3.7723503922735233`*^9}, {3.772350493664414*^9, 3.7723504939596043`*^9}, {
   3.772350594743127*^9, 3.772350602923962*^9}, {3.772350918783649*^9, 
   3.7723509227382927`*^9}, {3.7723576573175173`*^9, 
   3.7723576832168407`*^9}, {3.7723584848983927`*^9, 3.772358539922626*^9}, {
   3.77235876534759*^9, 3.772358767453516*^9}, {3.77235883525349*^9, 
   3.7723588353841743`*^9}, {3.7723590485874147`*^9, 3.772359092180375*^9}, 
   3.77235945892064*^9, {3.772387964781188*^9, 3.77238796502487*^9}},
 CellLabel->"In[45]:=",ExpressionUUID->"056f3927-ed5c-4691-a1d7-37b04defff55"],

Cell[BoxData[
 TemplateBox[{GraphicsBox[{{{{}, {}, 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.368417, 0.506779, 0.709798], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVxXk41HkcAOB5rB5rPSYpuRbRipQktEg+PkliaDfJ9aSmsM5RUSGrMLlG
zQzm/Jm+OdY+TOnyNM6imnpyFB1a26xiJVdqUrQ2bbv7x/u8lgcOBsVo0Gi0
wP/8v21D+swe7UWeeU70vZ9oGl6M+IAG7YARTwMXRayKPeQVcaq04gXjved8
entMI1vtNVziIu5ifAWTJsroGjYNa4wK9DsYS2EkcVu0jK2HWn849LUwVgJT
oPVSzF6BTZs1HZoYTkBsx6uF7PXIT05+rmB4w/WUAbaA7YX+aJd4zDQYJG5z
z5guOzHfo0Ve8SISBKsymKYl4RhS/GxiX3MU0PxC45KHmcjfrZL0TsaBrFdT
4LIrBmMl8on82yyQ9ry5Rs+IRyfZCbvFuYchi0pkPPjIwtfVuSZGy49A9+MV
5++tP4yzYyyVMusYcAzkaxogFfnPqLr9xhnQ2r4wqpd5FDO1L22vKswE8xpW
eRo/DS3UZlYVU1lgaleswyzPwPNbYjV9MRs4hinXWh9l4sE1TkO8Yzlg75UW
seJlFv48de5W351cqByv+nZBLxtjtdy3qOmnYD6bOVFrl4MZN33lw5gHkxvr
nSIcc7E/aM6om5sPBpXPTXx2sDEs1+hQ/p0C+M323EnnsFNoXaS677a8CJ5E
b82JTM/D3qnue91MDhA1MH3L8rHitcMFB2Ex0KiuNpfyAmxVOouKR06D7KOb
IrK5ED+718c1m3KhL8XVnqMsQjPdXQN3AnnwaPOZnLQuDnqrAmxGA/hAv343
hbwqRp0vdSF9Y3z4M5w33Pn2NKZzQ4saeSXwVu4xvW4xF590l3/DdS6FbobX
0+DVPPQeMM4OfFAKFnnxBIz4GPrj+F9fHy2Do57WFe9+5eNJsTqmXk8Akvza
Ieu1JSg6zbPEOgFUzT3ltVwswaDpE0N9O4TQ2ZGmm+RRipZmtWV+E0JoMx7q
PNdRisImS/8bfBH09LQEq38oQ7OV/R9M7cVQS480v91fhmWbLN7E3RVDQro8
ZC5MgBfbmyevJEnAzyMhctuYAAf1uaoRLSnMct9VGSUJ8crDq71LLkvhxqYp
hdaCEFM9k/UddlPA0uO7hRaJ0O94lJ/LWwrGUxf+ttUVo9Kyjv29qBxyHFUc
Q5kYgyS2V602yGA0Kye80FKCvW2eDp+7ZLB+5pEo/BcJvguoufT40FmQdIiq
fDdI0eaCz3dlOgTum1lNNyqk6Cp2tdhGJ9B9XuSv1SzF7blrTOb1CNxz06kN
aZViXNgSvX3LCdzc/eHAbLsU5ZqDn+wsCSjO3B3Y0ClF+8gjj29tJEA+Jyjr
VVJ0olfnzuwncHCwobyaRuHWeWFWTTQBVqLt/IwGhcEvC9PDYgkkzMtCtiyi
8EhLcvKNJAJRBvlLhrUpVPzkHsFJIxAWGFZgvpRCt/aHjlZnCHhdX0iRWlPo
L1eu7ecR8GQc7hu3oTBC2GhTWEpg0++j9q52FB5PPGv2RkzAefbB+NN1FLYZ
xmu3VBFYbV+5d5krhT0aezRZNQRWtRq0RblTODi944tFLYGVfhzjBg8K/7nt
PJtXT8A8JvXJTqSQfslG7XqZgMn7McdKbwrNKZOpqasEDHP28NQ+FK7L031F
rhFYtvjha9hOIRyiDe9sIqB/1sef50/hv9LnQ80=
          "]]}, 
        Annotation[#, "Charting`Private`Tag$3465#1"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.880722, 0.611041, 0.142051], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVj3s0lAkchm0pO1mzsmflNjsNi063WZdKJT8/umBQSbFWx7YKXQzJZZJL
xm0YZubDGMRnsLbGtSi3bG2TOssUM9vqyKxiVRJKJKlk7R/vec7z3/Myfgnz
OrZEQ0PDY3H/c00jZ8qfsswhqU6+pOvjF46s4+6NFPdhh5LgRFN+/KCjX0qO
9Alr2uHFQjpdfW7ScYjYJOliLQWVVO7XyNHASoN0vT9Y30Dm1sH36RG6qPUP
U9nGMoOm/sTVgUGrsWWHJrOFZQNaTXbl9Z4/oIjNftzEcgaTEKeJGIYjuuHa
k9HG3nA4SqmqjtiHafZtVdInhyEzotlctuCDh/j9owGtgRBXE6JaZxWAooPq
gp6XIaD0mFl+MDgQgwuqRtNuh8LHpcMWy4aD0KY4Ye3X3NOwv0ky1194Ascr
uEYG+pGgtf/h1DNrNs6MhKo74qNBXm+n5svDUdRfJDtieBaqd6a/i+ZH4DlK
vUs57xxYlAoPtFtHIn2SZiodi4dinwP6A91RWO0UrLkHz4PzaNXIbGYMhq2z
GRRGJ8Guy8KpUIuzGDdWKlfe4cLy9ensr4ZiMVhrm9MkNQUEjk/jHgrj8Oyt
PVVDmArDiiHOOC0Be73eGSgEaWCmcrfa3ZmIvlyD8LQ76ZCaxMvslp1H8wz1
/a36GbD8qM6MPC4Je8YUfyp+zoSHLp6PJl24KB1n1jDFfPCabthbppmM1zts
8/nDWaD5RFg+fi8Z57fVhrQaC+BNl2tjT0oK0nQO9N3xEELZm+r8bzenorPa
3fKZuwhO1t5o+zSVitoLskPKERFQ5u+5U6VpyBH4ZDQLCbjWbf/Cwjcd/1Zc
WCGwzYEyo1l6GIWHzn2G5z26c6BeeUqpqOahz74X77+MyoUVv19wa/fNwETJ
5LFa3Tx42buUFjCbgflZQgbK8sC6rM0/RJSJXhMJg0pPMWw/se/fxzZ8ZNAu
5bqOisGr3aTXpIuP4haG2w1RPoQ+V1+vO5GFNLPet8YbJGAt83Qt/JyFudvp
r0LuSuD4yqm+v7Kzse5m68srpwpA20Q0bW8pwAE9gXpYqxDmFHXehVcFeEXV
0LPyciF4+M0ZN7OEeMaBrcc8WAQbk3V2dQ4I0TU20HXT6yLg3BowU1uJsIMh
S96SfwEorIreGY4IvQrWNJhaF8Nd5m+zFa0i7Gl3YM53FcP9R4oaw1kRvnGv
rH8QXgKxnp1BhlsItKzZ9X2uNgnab4vqWsMItJPY0XdTSZC2zN+sCyfQhbvO
aE6XhPV0X17FaQJDfFfqBuiT8Elryj7rDIFVmgMf1zJIcBst3+gXQ+CGw5EP
5JtJQO+iDyMJBNpQK7hTR0j4oJ8f0JlN4M45cXzlURLsbP17rgoI9H7K4/gG
kxBdXj1RKiQwso3NvnGKhFU7FLZRBIFNQdv8MmNIUGY4PzMWE7j1psrKNJsE
5w8moh+LCXSr6ljfKyRBz0V9F0oI9BM3W/JyFvu2kxfNSQJjT5bQXklISJZR
giZLCWxfdZzSVk5C/vO93UkVBN5b4q8ZWknChDnV49ivBA5MeC7QL5HwU3vs
GZdKAj/ftp1JrSVhuIPZpHORQGq95aTdZRLijAZVrxf9uyKjsbEGEqazqXmq
SwRuTNV5Tl4joY/H/tQgIxDCNYb2tyz+cWhckVdF4H+FLVPA
          "]]}, 
        Annotation[#, "Charting`Private`Tag$3465#2"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.560181, 0.691569, 0.194885], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwV0Hs0lAkYBnBJK4mkm7STlCKSU0wprdcbJUaUXIaaYiMKI/coGYNhvk8I
Y1xmvqRVKCkSqaMbXZRcyi5NhaUUxdBlS4v99o/3POf3z3Oe8+r+Huzsq6ig
oLCfvv/ToOrY2D6VGZYKCjyddIvpVqzDDlUqDn2WChzlC/KiXivPxMzCbtZn
S576sKnT2VGr3tNMcRNrOliVu6dvp6ZhsVay5h3WPOgZ37LXJXcuKr8yaa1j
rYBlUYWOyam6WPubkkktyxQMtobpa0Wswwwu9811ljUE2v568Ol2RHs0DIhc
4gKzJzWwN383CrbUlRV2cyAywKLn0BkPdCNffjhw4yC85OjbC0VemOEqy20Z
9Ifwh3Yu5Xwf9Mst+yC4HwQlqq/U1yr7o6nkpOEcfghMP94+2nYlAD+e42tr
LQyHjTXW2ZyVwfh1IEjWEBsJnc3P6zjLQjDjZX6p9+JoYC9gzqq6E4rHVSp2
FKUchzCbxa5i83DUkTOWFw7FQnWBx01LKgIvbvVTskUeNIzwWgIXRWGwkWlP
emQ8tBeUNLpzj+GJoTP3Whv50POBH5faFY1+ypu3ytUTwSa8OOHlhuMYfde2
rBeT4G9uvXlX4gnscP6m9SRNANV1IT6Oo7HI5msdFTQmg07jzn/2QxyuFMqa
Ny0UQrL3nOGvu3jYMvTk0RMvAuS9jkLZOA8LP5pcMhGRcDYl/raDKB5vNpjl
kH2pEOgqcBrbyMeJzeX+N5akgXuG6MhYEx8Zans6G3emA/Oj/BfjPQloLXPQ
f+uQAUzHsxylNwmoOlXq1jqQAW6rzfj5nol4LM1dWJN+GrZMzuCEvErEF08K
ZqWZZUJajPnPd25JaN25mLfzWSZEWN9jvW5IQvdd77/PjMgCnWvF16s2CjBO
LPct18gGPZlP7g9KgDmp6bpYmg3fPAdi1ygmo/Onkz2tjiLQspl+81RQMuoy
SrLsPoigKIz5Zaw5GUW1uvb1GTlgZR1tOm99CjJWdHxZYiyGfb738izJFMyy
0Bn2fyCGPuMH435vUvDy7RuDVwNzwaNdAlJTIb7WTJP1KedBp31lr8VJIV5t
q2yZeyUPJrpsmfzHQgyz5GqauObDNPP4vAQ1Au1iDtoxR/LBw8G3IcqWwAbd
0oSNOQUgaMtSX8Qj0DnXoHL5egmM9Q+K314msOWWpclEkwRsXPT4A38SqOFu
JcxulsD6fv+p5L8I3DWKfUatEkiM2a1q0Elg+8ptYs8OCdzitjw71EXgizTW
VE23BKx97it2ywjs9GK3hH6RgNJ4pUZdD4E9SqHB7xlSeMP/6Ko9SOCoQ3HF
86NSOEGdn8P7SSD3WVS5ZZgU/vj8c3ic9qCT/cXSCCnYjBj5R/xLYL/zyPm4
GCnEB77XPzxB97M3UUaJUtCJD5Cxpgis92lOTRBLYamep/KkIompJ74eNq2X
wi+yLq0Fs0hUmXzoR92RQo1CgXYCbUFcvq/KfSlQ3gXX5LR5fPDufiiFUr9V
V5tUSQxLEbLJNil0CSaaotVI9Mxm2Pb1S0G2/GhkgwaJ+pe26WWpUmC9189Z
dRGJ5mJzne3qFDzXiOdH0d7BN9L+oUEBx8jRpI+2P3uuxoGFFBTNuDlZq0Vi
mdLrn4a6FFQRg6H7tUk05oQ/v7eBgg3V44UZDBJN1c/xx7wpaI+9nHN6BYk2
P0SxxT4UWGinbf5E26U/5Rjbj4KS5n72Dj0Sw+u43PpAClpGKpQmaF8/tNmT
iKJgPufRQq9VJG663bZu+SkK7t7Yu0BtNYn2ZQ1rOtIpmFlkz/ai7Smq0U/J
pMA9JnN+Fe2YACljWEyBU8jb726GJN5adFilroiC4GvnuLlGJD5V3KcUVEzB
HbZvzADt158cp3RKKJB0ySeZa0icvG/2NamcAuYAI7SNtnqFvtz8CgWaFQVe
S41JXJqvPTRUScGtxwUPj9Bem6T2jqqm97nMvnCdNhxV6N1dS/8nSDY5bS2J
/wEUzMJ7
          "]]}, Annotation[#, "Charting`Private`Tag$3465#3"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.922526, 0.385626, 0.209179], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVjmk8lAkcgCXWWhnXVpimybUsRTkiNP/+OcJgJUXzY1WaRUop95XcVzPj
mHlf4l2rJFc52kLlCr925Gir3SK5OkQHlS2tY9sPz+/58nx4NA4dd+dKSkhI
uHzlf+s1RL73lpVmJSU1GJoVSu5gBzo3yDpPsCTC/9xlqzi2g5OcWzLC/sBK
COKK5xVmd4zlmBFi9ko4N65PFSmswDLVNOU2tgqMiaQ3LdGUUOaJ0UAzWwtC
xg+x+uU1sHG7lFEj2wRWB81yIuW2oCA4+Ok1tg2sD/us2PcNohPqB4XTPWDU
VxxBNOzGVOvmypIRH6huk+hdyeLgvqzBV75NftAXoBis5XQQBXuHyP6pAOCx
+G2pYi76k5WvUm8fg8zyOVmt+UA0KYrXV0gMgcXuI91KscH4+nyiuuqaUDju
LmPpcDIE514eG+qMC4erF/o+mcafQsFgYcVBtShQkFFRvNYQhjGyVxxK02Ng
VrR0xLA9ApkzDM2S6TgIVX8eJtkRhVU7/aV2YQLQq7orH0/H4HEDk1F++BkQ
+5mzzN7GYez0rx0DXYnw8Xa8O29VAvrLWO6coSXDsz/MVQroZzCqfVflGKYA
323GznBtIj50/0e1h5cKDB173dYfktArUfVEalcaMA3UXJc1k1EnY6h325oM
eGO71WTZOAX7p3vu9BzIBKLcQ7veOhVLXhtVGwmzgJZrNeNtkoY3Ok1FWRPZ
8Fk7eip2ezouWtYENNF5ADoXbjWZZiBDfs+jLhc+JG9bXfBBJxNthpx1nzsL
YOGy0dR+/SyUW67YN/BSADmqV9u11mVjJM8z4zo/B+hPTRyfbjiLD3rOfccz
zYXORZQ6zeShzSO1BJe+XOhIWyjvkeOjp9vk52/D8oA2SnfzlRbgaWKGW6OY
DwsZp4zbYwQoyuZrYEU+6ChzlK3eCdD9TfzogKsQmqq5gxNBOajBuJTn+EoI
B/Re2JEjOShs1HBqEYiA5Spz18c7FxlaDz/SNxFgfOd+LA7kYp4V821ANwH8
2NhVbMzDy61NU3VHSSj93sIJb+ThsDJvaEKmAKbStdSiNuZj3b36fqXaAlhr
y7h4+mI+nmIFKxvtLYTunbzSEoYQHaP9HM3eFUJ68eoQ3WwhdmpUJJmLzkGs
3fVWz5UidCf16jWNi8BDM5IjcVKE/TdZRoviItj9sd1rcliEs85lV+6fKIZh
gwfJe90J1K22086To8C7Vtrp5xYCLQgLpj2NgnamiANtBDokGqjPK1KghNG+
GzoIDPBSUvRdQ0F0wXb78S4CK6WG/9XXoODwb7Vp3F4CN/mE3u/YSkGhvMNg
4BMCTWjnE98fpCDP07ol8AuBtvPCuLLDFPClbPY7LhDo8Sw90sufggi/4km9
JQJDm4ODW45SQDN3efZyBYnXfrHkZEZQsM+s4CeuLInbWu9t0TxLgR1vfM5H
lUSnys6ND/lfe472X9bqJHKE13XTcyno+Lupkr6OxOigYsZbgoKZyZHNj5kk
3lwbKNtcSgH3kGrdHl0S70p6Sx0r+/p365OJ8Y8kDr9xXWZeoqBGrF+laEDi
0m3TuZQaCsq/ZIX3GpJIu6I7Y1FLQRqjWVy1mcT1herT0/UU7HlvpZJpTKJh
ivwL6ncKVkUpuQWYkggnJMZ2N1LQUGd5xn4rif8B8hAvSA==
          "]]}, 
        Annotation[#, "Charting`Private`Tag$3465#4"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.528488, 0.470624, 0.701351], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJwVkGs0lAkYgGdlV2qz0lkpZ7JuKVZT6YKY15siBsnmsk4lKpQxxv1Sg51B
4zYz3zAzoj5iCbl3cqskl1jaKNkjk9vRas1E0ratU7Ltj+c8v59HPyDM45wK
hUJx/cL/3nYrbvGE+tf01eFdQa1hFDvGeZdb6i7TdOkjdsvbreN2viniognG
O/rB4Wdxmlvn7aaIvbI+xirY7flRUmf02a5U57JWO2MDnMmK1TMw0EC1F7TB
VoYhXJCqjVKpW7DZVpXWzLAARty7Yv0NO1DEYo03MuzhvUH9KdqKLTqjaUiM
7nEY3Bf2Ruf7o5hm01pZNHESIrM6jZo+eKFX5uisX8sZ2LgdfB+YnUKRpzxv
QBEM9/mR2W8XAzAor3I2rTMUZCMZQ7m8QLS4mmj6HTcc3mlbKew9L+DrEu5m
He0oSOI7WBaqs/D9q1B5FycGYt1NIz61sFE0ml/hvykeFOF17WxBBF5Urz1S
zL8Iaydmrbj0KNRboBoUKTlg9YflQ+OpaLx5MEjVEZMhbf7ogdHrsRhmZjEp
jPkF2NTgup8hHi8pCzsGu7nA6bPsUywlYJCa9cEFjRRoVA2IqSq7hPEPHCun
MBXGwfsQa1ciDnv8o9MvSAO/nmoT1ZEk9OHqsNO6L8PEsSl5QU0yGqfLf7fS
Tgdbh09xhYm/4ICyv7f/dAaUE0/tPjhyseg1rYomyQTKc8EPVl/x8E7XHmnm
dBbUL+n4HO7h4bJ1dXCLrgCKbax7MzgpSF3300i3qxDsDeN/7TBPRXu5i8mf
LiLwnHNlbFek4tqVCq/BVyJ4ajvic1qWhnEC7/QmIQF6biIvB/fL+Ky/YI1g
jxhuyxMVcyp8tB/ZlOz6WAzi6aQWyg0+erv/9e/q6BzgSstCDL3SMUm2cK5a
MxdmZDuR8iEdpVlCfazIhdk2McdGnIEec4mTg24SsLSgJs/tz0R9anmO06wE
Ei+tBL4ZyERJs75zm0gKPcrnuTPhWUg1HP5b11wGAopEwfkmG3MO6M0HP5QB
v8BmdZ0kG2vutyjqmXnQq3DMP7BDgGNaAvm02hU42/Fb9+cWAdY/aRhYX3cF
OCMzHgVuQoyks7RonvnALku+e31MiE4JZ5z2vsmHsqaGENudIuzSr+DtlxbA
msh2Sk+0CD3ytjUY7L4KWkq1KP/bIhy4S6ct912FYNvGyaRFEb51Ka0dYl+D
l0tPS0ItCDSpOmyUs5aEomi/niYWgZYySz0HDRKsloVV98IIPMI127ykSYIR
rzeuk01gsM96TT9tEpTpxx4+jiCwUnXso6k+Cc68zr0vYwg0Pxk11LGPBLqx
SeCaJAItNEq4i/4kqGyd+/aQkMBDSxJO6VkS2mOpQw4iAo+/5Mf5BJEw0OOZ
7EQQGNXKYrUxSZgKGBe65RDYGGjtmxFLwlDKUZ63jECr+092GWSTMB9bf+ds
IYHOlV0/DgtJaKvxXzxXRKCvpMmELyZhcnqTZtB1AhNCrlHnZSRoOxVTzpcQ
eHfjefXWYhI8VrQaQ24Q+EjlhGpoKQlV21YZM8sJHJtzW9ErJ8HL/WMCs4LA
z5173qdWf/mTv+oF8yaBGrUmC5Z1X/rubXjFrCJwS/5mpbKBBPNxs+fMagJ3
pK6bIW+TwFpmNDBrCAQ2ZepYMwkU3agIZi2B/wGdgUaw
          "]]}, 
        Annotation[#, "Charting`Private`Tag$3465#5"]& ]}}, {}, {}}, {
    DisplayFunction -> Identity, Ticks -> {Automatic, Automatic}, 
     AxesOrigin -> {0, 0}, FrameTicks -> {{Automatic, 
        Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
        Charting`ScaledFrameTicks[{Identity, Identity}]}}, 
     GridLines -> {Automatic, Automatic}, DisplayFunction -> Identity, 
     PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.05], 
        Scaled[0.05]}}, PlotRangeClipping -> True, ImagePadding -> All, 
     DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> 
     Identity, Frame -> {{False, False}, {False, False}}, 
     FrameLabel -> {{None, None}, {None, None}}, 
     FrameTicks -> {{Automatic, Automatic}, {Automatic, Automatic}}, 
     GridLines -> {Automatic, Automatic}, GridLinesStyle -> Directive[
       GrayLevel[0.5, 0.4]], 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None, 
       "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
           (Identity[#]& )[
            Part[#, 1]], 
           (Identity[#]& )[
            Part[#, 2]]}& ), "CopiedValueFunction" -> ({
           (Identity[#]& )[
            Part[#, 1]], 
           (Identity[#]& )[
            Part[#, 2]]}& )}}, 
     PlotRange -> {{0, 0.03}, {0., 1.1027349056798199`*^8}}, 
     PlotRangeClipping -> True, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.02], 
        Scaled[0.02]}}, Ticks -> {Automatic, Automatic}}],FormBox[
    FormBox[
     TemplateBox[{
      "\"Linear\"", "\"Landau4\"", "\"Experim. low rate\"", 
       "\"Experim. medium rate 2.81\"", "\"Experim. medium rate 0.292\""}, 
      "LineLegend", DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #3}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.922526, 0.385626, 0.209179], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.922526, 0.385626, 0.209179], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #4}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.528488, 0.470624, 0.701351], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.528488, 0.470624, 0.701351], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #5}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.24561133333333335`, 0.3378526666666667, 
                    0.4731986666666667], FrameTicks -> None, PlotRangePadding -> 
                    None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.368417`", ",", "0.506779`", ",", "0.709798`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.368417, 0.506779, 0.709798]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.368417, 0.506779, 0.709798], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.587148, 0.40736066666666665`, 0.09470066666666668], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.880722`", ",", "0.611041`", ",", "0.142051`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.880722, 0.611041, 0.142051]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.880722, 0.611041, 0.142051], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.37345400000000006`, 0.461046, 0.12992333333333334`], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.560181`", ",", "0.691569`", ",", "0.194885`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.560181, 0.691569, 0.194885]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.560181, 0.691569, 0.194885], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.922526, 0.385626, 0.209179], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.6150173333333333, 0.25708400000000003`, 
                    0.13945266666666667`], FrameTicks -> None, 
                    PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.922526`", ",", "0.385626`", ",", "0.209179`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.922526, 0.385626, 0.209179]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.922526, 0.385626, 0.209179], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.528488, 0.470624, 0.701351], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.3523253333333333, 0.3137493333333333, 
                    0.46756733333333333`], FrameTicks -> None, 
                    PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.528488`", ",", "0.470624`", ",", "0.701351`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.528488, 0.470624, 0.701351]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.528488, 0.470624, 0.701351], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}]}],
              "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2, ",", #3, ",", #4, ",", #5}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", "None"}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{{3.772359078102365*^9, 3.772359093552746*^9}, 
   3.772359238908114*^9, 3.772359460319889*^9, 3.7723598602237873`*^9, 
   3.7723599028132353`*^9, 3.77238796640283*^9},
 CellLabel->"Out[46]=",ExpressionUUID->"e3c621df-be3b-4603-b40a-7412f9c1b9c1"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"sbst4order", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"EL", "\[Rule]", 
      RowBox[{"8.133486894830394`*^15", " ", "b"}]}], ",", 
     RowBox[{"FL", "\[Rule]", 
      RowBox[{
       RowBox[{"-", "1.1366203039813472`*^16"}], "b"}]}], ",", 
     RowBox[{"GL", "\[Rule]", 
      RowBox[{
       RowBox[{"-", "1.0182884423600852`*^15"}], "b"}]}], ",", 
     RowBox[{"HL", "\[Rule]", 
      RowBox[{"1.4227609644541294`*^16", "b"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sbst4order", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"EL", "\[Rule]", 
      RowBox[{"1.`*^12", "b"}]}], ",", 
     RowBox[{"FL", "\[Rule]", 
      RowBox[{"1.`*^12", "b"}]}], ",", 
     RowBox[{"GL", "\[Rule]", 
      RowBox[{
       RowBox[{"-", "2.2584896459158383`*^11"}], "b"}]}], ",", 
     RowBox[{"HL", "\[Rule]", 
      RowBox[{"1.`*^12", "b"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Stress", "[", 
    RowBox[{"d_", ",", "a_", ",", "b_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"Stress", "[", 
     RowBox[{"d", ",", "a", ",", "EL", ",", "FL", ",", "GL", ",", "HL"}], 
     "]"}], "/.", "sbst4order"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.772359383331572*^9, 3.7723594212030354`*^9}, {
  3.772359650503766*^9, 3.7723596846334953`*^9}, {3.772359933431589*^9, 
  3.772359958382104*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"702e44d2-c727-4e8b-aaf1-6242b6d58b5a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sbst4order", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"EL", "\[Rule]", 
      RowBox[{"-", "3.3521917180042726`*^17"}]}], ",", 
     RowBox[{"FL", "\[Rule]", "4.693649596145145`*^17"}], ",", 
     RowBox[{"GL", "\[Rule]", "4.194849811670641`*^16"}], ",", 
     RowBox[{"HL", "\[Rule]", 
      RowBox[{"-", "6.142447703972079`*^17"}]}]}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.772388015623661*^9, 3.772388028865151*^9}},
 CellLabel->"In[47]:=",ExpressionUUID->"db643e60-c07e-40d8-b7e1-1a076806c001"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dmax", "=", "0.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"GraphicsRow", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Plot", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Stress", "[", 
            RowBox[{"d", ",", "a", ",", "b"}], "]"}], 
           RowBox[{"(*", 
            RowBox[{"/.", "sbst4order"}], "*)"}], ",", 
           RowBox[{"ExperDataStaticExtInterp", "[", "d", "]"}], ",", 
           RowBox[{"ExperDataMedDynExtSmooth", "[", "d", "]"}], ",", 
           RowBox[{"ExperDataMedDynExtSmooth2", "[", "d", "]"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"d", ",", "0", ",", "0.01"}], "}"}], ",", 
         RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
         RowBox[{"ImageSize", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{"100", ",", "100"}], "}"}]}]}], "]"}], ",", 
       RowBox[{"Plot", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Stress", "[", 
            RowBox[{"d", ",", "a", ",", "b"}], "]"}], 
           RowBox[{"(*", 
            RowBox[{"/.", "sbst4order"}], "*)"}], ",", 
           RowBox[{"ExperDataStaticExtInterp", "[", "d", "]"}], ",", 
           RowBox[{"ExperDataMedDynExtSmooth", "[", "d", "]"}], ",", 
           RowBox[{"ExperDataMedDynExtSmooth2", "[", "d", "]"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"d", ",", "0", ",", "dmax"}], "}"}], ",", 
         RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
         RowBox[{"PlotLegends", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
           "\"\<Landau4\>\"", ",", "\"\<Experim. low rate\>\"", ",", 
            "\"\<Experim. medium rate 2.81\>\"", ",", 
            "\"\<Experim. medium rate 0.292\>\""}], "}"}]}], ",", 
         RowBox[{"ImageSize", "\[Rule]", "Medium"}]}], "]"}]}], "}"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"a", ",", "1", ",", "20"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"b", ",", "0", ",", "1", ",", "0.01"}], "}"}]}], "]"}]}], "Input",\

 CellChangeTimes->{{3.772258362208325*^9, 3.772258403530744*^9}, {
   3.772259320442834*^9, 3.7722593490827713`*^9}, {3.7722593915113907`*^9, 
   3.7722594035006533`*^9}, {3.772259441649233*^9, 3.772259473704123*^9}, {
   3.772259604056508*^9, 3.772259683275116*^9}, {3.772269395439508*^9, 
   3.772269409168303*^9}, {3.77226944871436*^9, 3.772269448835977*^9}, {
   3.772269511462352*^9, 3.7722695462823877`*^9}, {3.77227112177039*^9, 
   3.772271121883951*^9}, {3.7722711607953367`*^9, 3.772271177035329*^9}, 
   3.772271228797082*^9, 3.772276649581072*^9, 3.772276707017564*^9, {
   3.772335666337817*^9, 3.772335691055481*^9}, {3.7723357575664473`*^9, 
   3.772335757832345*^9}, {3.772347824949437*^9, 3.7723478355645723`*^9}, {
   3.7723481554768753`*^9, 3.772348157582074*^9}, {3.7723482399541483`*^9, 
   3.772348240086149*^9}, {3.772348270658585*^9, 3.772348335494632*^9}, {
   3.77234847120266*^9, 3.7723484726689157`*^9}, {3.7723503840311413`*^9, 
   3.7723503922735233`*^9}, {3.772350493664414*^9, 3.7723504939596043`*^9}, {
   3.772350594743127*^9, 3.772350602923962*^9}, {3.772350918783649*^9, 
   3.7723509227382927`*^9}, {3.7723576573175173`*^9, 
   3.7723576832168407`*^9}, {3.7723584848983927`*^9, 3.772358539922626*^9}, {
   3.77235876534759*^9, 3.772358767453516*^9}, {3.77235883525349*^9, 
   3.7723588353841743`*^9}, {3.772359448658573*^9, 3.772359500867408*^9}, {
   3.7723596102355223`*^9, 3.772359610380136*^9}, {3.772359689726781*^9, 
   3.772359690063252*^9}, {3.772359964715962*^9, 3.7723599750853643`*^9}, {
   3.7723600072105*^9, 3.772360074409944*^9}, {3.772387994149927*^9, 
   3.77238800634783*^9}, {3.77238806657141*^9, 3.772388104747602*^9}, {
   3.77238815108239*^9, 3.7723881514788113`*^9}, {3.772388498034327*^9, 
   3.772388575711689*^9}, {3.7723887493940887`*^9, 3.772388790286489*^9}, 
   3.77238882217815*^9, {3.7723888580934353`*^9, 3.7723888654916487`*^9}, {
   3.7723889361238337`*^9, 3.7723890393924303`*^9}, {3.772389072563843*^9, 
   3.772389124783098*^9}, {3.772389163685472*^9, 3.7723891849516363`*^9}, {
   3.772389288399589*^9, 3.772389443959405*^9}, {3.772389492361594*^9, 
   3.772389506379928*^9}},
 CellLabel->
  "In[132]:=",ExpressionUUID->"554740c4-91f4-46cb-995a-dc2f477060ba"],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`a$$ = 5., $CellContext`b$$ = 1., 
    Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{
      Hold[$CellContext`a$$], 1, 20}, {
      Hold[$CellContext`b$$], 0, 1, 0.01}}, Typeset`size$$ = {
    720., {197., 204.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`a$44795$$ = 
    0, $CellContext`b$44796$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`a$$ = 1, $CellContext`b$$ = 0}, 
      "ControllerVariables" :> {
        Hold[$CellContext`a$$, $CellContext`a$44795$$, 0], 
        Hold[$CellContext`b$$, $CellContext`b$44796$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> GraphicsRow[{
         Plot[{
           $CellContext`Stress[$CellContext`d, $CellContext`a$$, \
$CellContext`b$$], 
           $CellContext`ExperDataStaticExtInterp[$CellContext`d], 
           $CellContext`ExperDataMedDynExtSmooth[$CellContext`d], 
           $CellContext`ExperDataMedDynExtSmooth2[$CellContext`d]}, \
{$CellContext`d, 0, 0.01}, GridLines -> Automatic, ImageSize -> {100, 100}], 
         Plot[{
           $CellContext`Stress[$CellContext`d, $CellContext`a$$, \
$CellContext`b$$], 
           $CellContext`ExperDataStaticExtInterp[$CellContext`d], 
           $CellContext`ExperDataMedDynExtSmooth[$CellContext`d], 
           $CellContext`ExperDataMedDynExtSmooth2[$CellContext`d]}, \
{$CellContext`d, 0, $CellContext`dmax}, GridLines -> Automatic, 
          PlotLegends -> {
           "Landau4", "Experim. low rate", "Experim. medium rate 2.81", 
            "Experim. medium rate 0.292"}, ImageSize -> Medium]}, ImageSize -> 
        Large], "Specifications" :> {{$CellContext`a$$, 1, 
         20}, {$CellContext`b$$, 0, 1, 0.01}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{777., {289., 297.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.772389414529623*^9, 3.772389444600236*^9}, {
  3.772389497293901*^9, 3.7723895068329897`*^9}, {3.772389565920397*^9, 
  3.7723895724813623`*^9}},
 CellLabel->
  "Out[133]=",ExpressionUUID->"b903bbf5-61b4-4e3b-9f52-3ab5ab86cc3f"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Compression [needs correction]", "Subsection",
 CellChangeTimes->{{3.772271247937764*^9, 3.772271256144603*^9}, {
  3.772276748457902*^9, 3.772276777668331*^9}, {3.772360150836989*^9, 
  3.772360156778953*^9}},ExpressionUUID->"999aa3f8-fd30-4a23-9ce8-\
fb12ac37d4d0"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Data", " ", "from", " ", "Zhang", " ", "2016"}], ":", " ", 
    RowBox[{"intermediate", " ", "compression", " ", "rate", " ", 
     RowBox[{"(", 
      RowBox[{"0.1", " ", 
       SuperscriptBox["s", 
        RowBox[{"-", "1"}]]}], ")"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
  
  RowBox[{
   RowBox[{"ExperDataSlowDynComp", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.03"}], ",", 
        RowBox[{
         RowBox[{"-", "20"}], "*", 
         SuperscriptBox["10", "6"]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.025"}], ",", 
        RowBox[{
         RowBox[{"-", "15.2"}], "*", 
         SuperscriptBox["10", "6"]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.02"}], ",", 
        RowBox[{
         RowBox[{"-", "10.8"}], "*", 
         SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.015"}], ",", 
        RowBox[{
         RowBox[{"-", "7.2"}], "*", 
         SuperscriptBox["10", "6"]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.01"}], ",", 
        RowBox[{
         RowBox[{"-", "4"}], "*", 
         SuperscriptBox["10", "6"]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.005"}], ",", 
        RowBox[{
         RowBox[{"-", "1.6"}], "*", 
         SuperscriptBox["10", "6"]}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}]}], "}"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.772276474470194*^9, 3.7722765656260977`*^9}, {
  3.772276784225371*^9, 3.772276797182976*^9}, {3.772276866572214*^9, 
  3.772276875991165*^9}, {3.77235889534264*^9, 
  3.772358900625572*^9}},ExpressionUUID->"47640f47-8f04-41d5-bd7b-\
352ac5d05a36"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Data", " ", "from", " ", "Li", " ", "2001"}], ":", " ", 
    RowBox[{"rapid", " ", "compression", " ", 
     RowBox[{"(", 
      RowBox[{"1500", " ", 
       SuperscriptBox["s", 
        RowBox[{"-", "1"}]]}], ")"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
  
  RowBox[{
   RowBox[{
    RowBox[{"ExperDataFastDynComp", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.002", ",", 
         RowBox[{"25", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.04", ",", 
         RowBox[{"50", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"0.01", ",", 
         RowBox[{"74", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.015", ",", 
         RowBox[{"80", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.02", ",", 
         RowBox[{"125", "*", 
          SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ExperDataFastDynCompInterp", "=", 
     RowBox[{"Interpolation", "[", "ExperDataFastDynComp", "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.7722768248115377`*^9, 3.7722768839753027`*^9}, {
  3.772276940417243*^9, 3.7722770007838507`*^9}, {3.772277085821553*^9, 
  3.772277122211185*^9}, {3.772277158841538*^9, 3.772277204713822*^9}, {
  3.772279199808434*^9, 3.7722792468293047`*^9}, {3.772279514640534*^9, 
  3.7722795175040827`*^9}, {3.772279570250833*^9, 3.7722795703807173`*^9}, {
  3.772358906453312*^9, 
  3.772358906586032*^9}},ExpressionUUID->"d695a0c2-6589-48a2-97ba-\
0695e9fa70f3"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{880, 911},
WindowMargins->{{79, Automatic}, {26, Automatic}},
Magnification:>1.25 Inherited,
FrontEndVersion->"11.3 for Linux x86 (64-bit) (March 6, 2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 159, 3, 85, "Section",ExpressionUUID->"1db9d287-f7c6-42b4-a9d4-606839815e86"],
Cell[742, 27, 14130, 337, 1658, "Input",ExpressionUUID->"10bc8bc5-de5e-4ac8-abc4-9d37283af475"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14909, 369, 212, 4, 85, "Section",ExpressionUUID->"d669c3bd-62da-412e-88a6-05ff878fb8e6"],
Cell[15124, 375, 1855, 38, 102, "Input",ExpressionUUID->"118ca0f4-45d3-4e90-9b82-736b8a259fed"],
Cell[16982, 415, 2018, 58, 243, "Input",ExpressionUUID->"ae2879c8-059f-423b-9b96-ae967dadb5c8"],
Cell[19003, 475, 960, 22, 156, "Input",ExpressionUUID->"faa411b5-4659-4640-a167-295e5a203224"],
Cell[19966, 499, 1044, 24, 69, "Input",ExpressionUUID->"1e414e05-8461-40ee-917c-b603a07c0f19"],
Cell[CellGroupData[{
Cell[21035, 527, 159, 3, 69, "Subsection",ExpressionUUID->"7289dff9-9280-4a82-94d7-3c7a33b67c65"],
Cell[CellGroupData[{
Cell[21219, 534, 8206, 211, 979, "Input",ExpressionUUID->"6334fc9b-4e81-4604-a00b-737281ea1177"],
Cell[29428, 747, 3433, 71, 290, "Output",ExpressionUUID->"d674ad7c-f351-4869-88a3-11d70336bd78"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32898, 823, 683, 17, 127, "Input",ExpressionUUID->"3227967f-a37c-49ed-9175-9fe4d431cc39"],
Cell[33584, 842, 371, 5, 44, "Output",ExpressionUUID->"872028df-9a71-42ab-9276-6fad752dd670"],
Cell[33958, 849, 374, 5, 44, "Output",ExpressionUUID->"05d1abf5-4dd3-4657-b7f5-a7f53c114a71"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34369, 859, 160, 3, 57, "Subsubsection",ExpressionUUID->"bc2aeae8-2015-4b44-b1f7-cb7995fb441f"],
Cell[CellGroupData[{
Cell[34554, 866, 3401, 67, 301, "Input",ExpressionUUID->"b6e0fb2b-21b9-4aa2-acd5-b4391d914e73"],
Cell[37958, 935, 1139, 17, 44, "Output",ExpressionUUID->"9b868566-770b-4b0c-80ed-f569c944dd08"],
Cell[39100, 954, 33308, 650, 279, "Output",ExpressionUUID->"1719ff9f-71ab-43a8-a205-d669b47b172b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72445, 1609, 2468, 46, 243, "Input",ExpressionUUID->"5beb9aa0-77e8-4bb4-a3b6-5f3301c24a15"],
Cell[74916, 1657, 3373, 59, 385, "Output",ExpressionUUID->"727a5820-5a2d-4908-bffa-ada58eb8ae88"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[78338, 1722, 162, 3, 57, "Subsubsection",ExpressionUUID->"a9047f7e-c4b0-42c2-9c13-8c3c44722a38"],
Cell[78503, 1727, 619, 16, 98, "Input",ExpressionUUID->"7cfd2991-148b-4576-b2aa-5eedb925db1c"],
Cell[CellGroupData[{
Cell[79147, 1747, 551, 12, 69, "Input",ExpressionUUID->"9821070a-224a-4597-957b-98b755efb9cc"],
Cell[79701, 1761, 569, 12, 81, "Output",ExpressionUUID->"e85a73d4-0965-49b3-a358-822c65a35ea0"]
}, Open  ]],
Cell[CellGroupData[{
Cell[80307, 1778, 496, 13, 39, "Input",ExpressionUUID->"50bd4332-2385-43af-af52-7715da6e6b9d"],
Cell[80806, 1793, 491, 11, 48, "Output",ExpressionUUID->"d9daccc6-b46e-4cf7-b024-d413425bc0f1"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81334, 1809, 922, 25, 98, "Input",ExpressionUUID->"1bd97a10-0359-4014-a19e-427f5e9e34a9"],
Cell[82259, 1836, 423, 10, 48, "Output",ExpressionUUID->"700709f3-e9c9-411e-b595-b7790eb3dce8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82719, 1851, 2498, 45, 185, "Input",ExpressionUUID->"056f3927-ed5c-4691-a1d7-37b04defff55"],
Cell[85220, 1898, 32585, 642, 279, "Output",ExpressionUUID->"e3c621df-be3b-4603-b40a-7412f9c1b9c1"]
}, Open  ]],
Cell[117820, 2543, 1469, 40, 127, "Input",ExpressionUUID->"702e44d2-c727-4e8b-aaf1-6242b6d58b5a"],
Cell[119292, 2585, 551, 12, 69, "Input",ExpressionUUID->"db643e60-c07e-40d8-b7e1-1a076806c001"],
Cell[CellGroupData[{
Cell[119868, 2601, 4412, 84, 359, "Input",ExpressionUUID->"554740c4-91f4-46cb-995a-dc2f477060ba"],
Cell[124283, 2687, 2929, 58, 609, "Output",ExpressionUUID->"b903bbf5-61b4-4e3b-9f52-3ab5ab86cc3f"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[127273, 2752, 273, 4, 69, "Subsection",ExpressionUUID->"999aa3f8-fd30-4a23-9ce8-fb12ac37d4d0"],
Cell[127549, 2758, 1902, 57, 109, "Input",ExpressionUUID->"47640f47-8f04-41d5-bd7b-352ac5d05a36"],
Cell[129454, 2817, 1862, 50, 139, "Input",ExpressionUUID->"d695a0c2-6589-48a2-97ba-0695e9fa70f3"]
}, Open  ]]
}, Open  ]]
}
]
*)