entropie3.nb

(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     76383,       2034]
NotebookOptionsPosition[     72753,       1907]
NotebookOutlinePosition[     73099,       1922]
CellTagsIndexPosition[     73056,       1919]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"k", "=", 
   RowBox[{"8.31451070", "*", 
    RowBox[{"10", "^", 
     RowBox[{"-", "3"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"m", "=", "18.0153"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"h", "=", "0.3990313224"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sym", "=", "2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ThAThBThC", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"h", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"8", 
        RowBox[{"Pi", "^", "2"}], "k"}], ")"}]}], ")"}], "^", "3"}], "/", 
   RowBox[{"(", 
    RowBox[{"0.00614517", "*", "0.01155002", "*", "0.01769519"}], 
    ")"}]}]}]}], "Input",
 CellChangeTimes->{{3.5098256238743076`*^9, 3.50982562531539*^9}, {
   3.509833466984908*^9, 3.509833484548912*^9}, 3.509833586208727*^9, {
   3.510085114781913*^9, 3.5100851150659294`*^9}, {3.510085291191003*^9, 
   3.5100852913900146`*^9}, {3.5111054983320823`*^9, 3.511105499291137*^9}, {
   3.5111055546043005`*^9, 3.5111055560323825`*^9}, {3.515828839284794*^9, 
   3.5158288615590677`*^9}, {3.5158288970961003`*^9, 3.515828939326516*^9}, {
   3.5158290505318766`*^9, 3.51582906685581*^9}, {3.515829922684761*^9, 
   3.5158299417938538`*^9}}],

Cell[BoxData["11360.325359684379`"], "Output",
 CellChangeTimes->{3.515880291091324*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"densityF", " ", "=", 
   RowBox[{"SparseArray", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"300", "->", ".98"}], ",", 
      RowBox[{"360", "\[Rule]", ".96"}], ",", 
      RowBox[{"420", "\[Rule]", ".93"}], ",", 
      RowBox[{"480", "\[Rule]", ".9"}]}], "}"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5158780490060835`*^9, 3.5158780574935694`*^9}, {
  3.515878198649643*^9, 3.5158782658144846`*^9}, {3.515880272599266*^9, 
  3.5158802799536867`*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.515829820448913*^9, 3.515829824084121*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"len", "=", "10"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dt", "=", "0.004"}], ";"}]}], "Input",
 CellChangeTimes->{{3.509825630350678*^9, 3.509825638184126*^9}, 
   3.509833466987908*^9, 3.5098335824335113`*^9, {3.5098452653667364`*^9, 
   3.5098452658707657`*^9}, {3.5111055081296425`*^9, 3.511105510567782*^9}, {
   3.5111082558928056`*^9, 3.51110826349524*^9}, 3.5111157978811827`*^9, 
   3.511115915921934*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"S", "[", 
   RowBox[{"T_", ",", "nMol_", ",", " ", "file_", ",", "rho_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "c1", ",", " ", "Ct", ",", "gT", ",", "j", ",", "h1", ",", "De", ",", 
      "f1", ",", "f", ",", "gg", ",", "Dw", ",", "gs", ",", "Sig", ",", "fy", 
      ",", "z", ",", "SHs", ",", "Wgs", ",", "hb", ",", "beta", ",", "Who", 
      ",", "Svib", ",", "Sconf", ",", "Srot"}], "}"}], ",", 
    RowBox[{
     RowBox[{"c1", "=", 
      RowBox[{"k", "*", "T", "*", "3", " ", "*", "nMol"}]}], ";", " ", 
     RowBox[{"Ct", "=", 
      RowBox[{"Import", "[", 
       RowBox[{"file", ",", "\"\<Table\>\""}], "]"}]}], " ", ";", " ", 
     RowBox[{
      RowBox[{"Ct", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "2"}], "]"}], "]"}], "=", 
      RowBox[{
       RowBox[{"Ct", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "2"}], "]"}], "]"}], "*", "c1"}]}], ";", 
     RowBox[{"gT", "=", 
      RowBox[{
       RowBox[{
        RowBox[{"2", "/", "k"}], " ", "/", "T"}], "*", "dt", "*", 
       RowBox[{"(", 
        RowBox[{"Chop", "[", 
         RowBox[{
          RowBox[{"Fourier", "[", 
           RowBox[{
            RowBox[{"Ct", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", 
            RowBox[{"FourierParameters", "\[Rule]", 
             RowBox[{"{", 
              RowBox[{"1", ",", "1"}], "}"}]}]}], "]"}], "+", 
          RowBox[{"Fourier", "[", 
           RowBox[{
            RowBox[{"Ct", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", 
            RowBox[{"FourierParameters", "\[Rule]", 
             RowBox[{"{", 
              RowBox[{"1", ",", 
               RowBox[{"-", "1"}]}], "}"}]}]}], "]"}], " ", "-", " ", 
          RowBox[{"Ct", "[", 
           RowBox[{"[", 
            RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "]"}], ")"}]}]}], ";", 
     RowBox[{
      RowBox[{"j", "[", 
       RowBox[{"f_", ",", "y_"}], "]"}], ":=", 
      RowBox[{
       RowBox[{"2", "*", 
        RowBox[{"y", "^", "3"}], "*", 
        RowBox[{"f", "^", "3"}]}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"y", "+", 
          RowBox[{"6", 
           RowBox[{"y", "^", "2"}]}]}], ")"}], "*", 
        RowBox[{"f", "^", "2"}]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"6", "y"}]}], ")"}], "*", "f"}], "-", "2"}]}], ";", 
     RowBox[{
      RowBox[{"h1", "[", 
       RowBox[{"f_", ",", "D_"}], "]"}], ":=", 
      RowBox[{"j", "[", 
       RowBox[{"f", ",", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"f", "/", "D"}], ")"}], "^", 
         RowBox[{"(", 
          RowBox[{"3", "/", "2"}], ")"}]}]}], "]"}]}], ";", 
     RowBox[{"De", "=", 
      RowBox[{"2", "*", 
       RowBox[{
        RowBox[{
         RowBox[{"gT", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "/", "9"}], "/", "nMol"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Pi", "*", "k", "*", 
          RowBox[{"T", "/", "m"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"1", "/", "2"}], ")"}]}], "*", 
       RowBox[{"rho", "^", 
        RowBox[{"(", 
         RowBox[{"1", "/", "3"}], ")"}]}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"6", "/", "Pi"}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"2", "/", "3"}], ")"}]}]}]}], ";", 
     RowBox[{"f1", "=", 
      RowBox[{"f", "/.", 
       RowBox[{"ToRules", "[", 
        RowBox[{"Reduce", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"h1", "[", 
            RowBox[{"f", ",", "De"}], "]"}], "==", "0"}], ",", "f", ",", 
          "Reals"}], "]"}], "]"}]}]}], ";", 
     RowBox[{
      RowBox[{"gg", "[", "w_", "]"}], ":=", 
      RowBox[{
       RowBox[{"gT", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"gT", "[", 
             RowBox[{"[", "1", "]"}], "]"}], "*", 
            RowBox[{"w", "/", 
             RowBox[{"(", 
              RowBox[{"12", "*", "f1", "*", "nMol"}], ")"}]}]}], ")"}], "^", 
          "2"}]}], ")"}]}]}], ";", 
     RowBox[{"Dw", "=", 
      RowBox[{"2", "*", 
       RowBox[{"Pi", "/", "len"}]}]}], ";", 
     RowBox[{
      RowBox[{"gs", "[", "w_", "]"}], ":=", 
      RowBox[{
       RowBox[{"gT", "[", 
        RowBox[{"[", 
         RowBox[{"1", "+", 
          RowBox[{"w", "/", "Dw"}]}], "]"}], "]"}], "-", 
       RowBox[{"gg", "[", "w", "]"}]}]}], ";", 
     RowBox[{"Sig", "=", " ", 
      RowBox[{"k", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"1", "/", "rho"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"4", "*", "Pi", "*", "m", "*", 
              RowBox[{"3", "/", "2"}], "*", "nMol", "*", "k", "*", 
              RowBox[{"T", "/", 
               RowBox[{"(", 
                RowBox[{"3", "*", "nMol", "*", 
                 RowBox[{"h", "^", "2"}]}], ")"}]}]}], ")"}], "^", 
            RowBox[{"(", 
             RowBox[{"3", "/", "2"}], ")"}]}]}], "]"}], "+", 
         RowBox[{"5", "/", "2"}]}], ")"}]}]}], ";", 
     RowBox[{
      RowBox[{"z", "[", "y_", "]"}], ":=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "y", "+", 
         RowBox[{"y", "^", "2"}], "-", 
         RowBox[{"y", "^", "3"}]}], ")"}], "/", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", "y"}], ")"}], "^", "3"}]}]}], ";", 
     RowBox[{"fy", "=", 
      RowBox[{"f1", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"f1", "/", "De"}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"3", "/", "2"}], ")"}]}]}]}], ";", 
     RowBox[{"SHs", "=", 
      RowBox[{"Sig", "+", 
       RowBox[{"k", "*", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Log", "[", 
           RowBox[{"z", "[", "fy", "]"}], "]"}], "+", 
          RowBox[{"fy", "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"3", "*", "fy"}], "-", "4"}], ")"}], "/", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1", "-", "fy"}], ")"}], "^", "2"}]}]}]}], ")"}]}]}]}], 
     ";", 
     RowBox[{"Wgs", "=", 
      RowBox[{
       RowBox[{"1", "/", "3"}], "*", 
       RowBox[{"SHs", "/", "k"}]}]}], ";", 
     RowBox[{"hb", "=", 
      RowBox[{
       RowBox[{"h", "/", "2"}], "/", "Pi"}]}], ";", 
     RowBox[{"beta", "=", 
      RowBox[{
       RowBox[{"1", "/", "k"}], "/", "T"}]}], ";", 
     RowBox[{
      RowBox[{"Who", "[", "w_", "]"}], ":=", 
      RowBox[{
       RowBox[{"beta", "*", "hb", "*", 
        RowBox[{"w", "/", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Exp", "[", 
            RowBox[{"beta", "*", "hb", "*", "w"}], "]"}], "-", "1"}], 
          ")"}]}]}], "-", 
       RowBox[{"Log", "[", 
        RowBox[{"1", "-", 
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", "beta"}], "*", "hb", "*", "w"}], "]"}]}], "]"}]}]}], 
     ";", 
     RowBox[{"Svib", "=", 
      RowBox[{"k", "*", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"Dw", "/", "2"}], "/", "Pi"}], "*", 
          RowBox[{"gs", "[", "w", "]"}], "*", 
          RowBox[{"Who", "[", "w", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"w", ",", "Dw", ",", 
           RowBox[{
            RowBox[{"1", "/", "dt"}], "*", "Pi"}], ",", "Dw"}], "}"}]}], 
        "]"}]}]}], ";", 
     RowBox[{"Sconf", "=", 
      RowBox[{"k", "*", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"Dw", "/", "2"}], "/", "Pi"}], "*", 
          RowBox[{"gg", "[", "w", "]"}], "*", "Wgs"}], ",", 
         RowBox[{"{", 
          RowBox[{"w", ",", "0", ",", 
           RowBox[{
            RowBox[{"1", "/", "dt"}], "*", "Pi"}], ",", "Dw"}], "}"}]}], 
        "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Svib", "+", "Sconf"}], ")"}], "/", "nMol"}], "*", 
      "1000"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5098076621849575`*^9, 3.509807664636098*^9}, {
   3.509807712034809*^9, 3.509807739029353*^9}, {3.50980796101605*^9, 
   3.5098079633461833`*^9}, {3.509825585893135*^9, 3.50982558807326*^9}, {
   3.509825622365221*^9, 3.5098256351109505`*^9}, {3.5098316000531254`*^9, 
   3.5098316300778427`*^9}, 3.5098317074072657`*^9, {3.509838610274087*^9, 
   3.509838612351206*^9}, {3.5100847175281916`*^9, 3.510084787852214*^9}, {
   3.5100848235612564`*^9, 3.510084826357416*^9}, {3.5100848593523035`*^9, 
   3.510084933171526*^9}, {3.510084981875311*^9, 3.510085025441803*^9}, {
   3.5100850621109004`*^9, 3.5100850780978146`*^9}, {3.510085123248397*^9, 
   3.5100851820937634`*^9}, {3.5100852310265617`*^9, 
   3.5100854669290547`*^9}, {3.5100855589893203`*^9, 
   3.5100855637435923`*^9}, {3.510085654596789*^9, 3.510085654873804*^9}, {
   3.510087973580427*^9, 3.5100879855091095`*^9}, {3.510088569140491*^9, 
   3.510088619097348*^9}, {3.511108270363633*^9, 3.511108285472498*^9}, {
   3.511115627342428*^9, 3.511115639970151*^9}, {3.5152025783926644`*^9, 
   3.515202578971698*^9}, {3.5154842856584587`*^9, 3.5154842917698083`*^9}, 
   3.515828409660221*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Srot", "[", 
   RowBox[{"T_", ",", "nMol_", ",", " ", "file_", ",", "rho_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "c1", ",", " ", "Ct", ",", "gT", ",", "j", ",", "h1", ",", "De", ",", 
      "f1", ",", "f", ",", "gg", ",", "Dw", ",", "gs", ",", "Sig", ",", "fy", 
      ",", "z", ",", "SR", ",", "Wgs", ",", "hb", ",", "beta", ",", "Who", 
      ",", "Svib", ",", "Sconf", ",", "Srot"}], "}"}], ",", 
    RowBox[{
     RowBox[{"c1", "=", 
      RowBox[{"k", "*", "T", "*", "3", " ", "*", "nMol"}]}], ";", " ", 
     RowBox[{"Ct", "=", 
      RowBox[{"Import", "[", 
       RowBox[{"file", ",", "\"\<Table\>\""}], "]"}]}], " ", ";", " ", 
     RowBox[{
      RowBox[{"Ct", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "2"}], "]"}], "]"}], "=", 
      RowBox[{
       RowBox[{"Ct", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "2"}], "]"}], "]"}], "*", "c1"}]}], ";", 
     RowBox[{"gT", "=", 
      RowBox[{
       RowBox[{
        RowBox[{"2", "/", "k"}], " ", "/", "T"}], "*", "dt", "*", 
       RowBox[{"(", 
        RowBox[{"Chop", "[", 
         RowBox[{
          RowBox[{"Fourier", "[", 
           RowBox[{
            RowBox[{"Ct", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", 
            RowBox[{"FourierParameters", "\[Rule]", 
             RowBox[{"{", 
              RowBox[{"1", ",", "1"}], "}"}]}]}], "]"}], "+", 
          RowBox[{"Fourier", "[", 
           RowBox[{
            RowBox[{"Ct", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", 
            RowBox[{"FourierParameters", "\[Rule]", 
             RowBox[{"{", 
              RowBox[{"1", ",", 
               RowBox[{"-", "1"}]}], "}"}]}]}], "]"}], " ", "-", " ", 
          RowBox[{"Ct", "[", 
           RowBox[{"[", 
            RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "]"}], ")"}]}]}], ";", 
     RowBox[{
      RowBox[{"j", "[", 
       RowBox[{"f_", ",", "y_"}], "]"}], ":=", 
      RowBox[{
       RowBox[{"2", "*", 
        RowBox[{"y", "^", "3"}], "*", 
        RowBox[{"f", "^", "3"}]}], "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"y", "+", 
          RowBox[{"6", 
           RowBox[{"y", "^", "2"}]}]}], ")"}], "*", 
        RowBox[{"f", "^", "2"}]}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "+", 
          RowBox[{"6", "y"}]}], ")"}], "*", "f"}], "-", "2"}]}], ";", 
     RowBox[{
      RowBox[{"h1", "[", 
       RowBox[{"f_", ",", "D_"}], "]"}], ":=", 
      RowBox[{"j", "[", 
       RowBox[{"f", ",", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"f", "/", "D"}], ")"}], "^", 
         RowBox[{"(", 
          RowBox[{"3", "/", "2"}], ")"}]}]}], "]"}]}], ";", 
     RowBox[{"De", "=", 
      RowBox[{"2", "*", 
       RowBox[{
        RowBox[{
         RowBox[{"gT", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "/", "9"}], "/", "nMol"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Pi", "*", "k", "*", 
          RowBox[{"T", "/", "m"}]}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"1", "/", "2"}], ")"}]}], "*", 
       RowBox[{"rho", "^", 
        RowBox[{"(", 
         RowBox[{"1", "/", "3"}], ")"}]}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"6", "/", "Pi"}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"2", "/", "3"}], ")"}]}]}]}], ";", 
     RowBox[{"f1", "=", 
      RowBox[{"f", "/.", 
       RowBox[{"ToRules", "[", 
        RowBox[{"Reduce", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"h1", "[", 
            RowBox[{"f", ",", "De"}], "]"}], "==", "0"}], ",", "f", ",", 
          "Reals"}], "]"}], "]"}]}]}], ";", 
     RowBox[{
      RowBox[{"gg", "[", "w_", "]"}], ":=", 
      RowBox[{
       RowBox[{"gT", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "/", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"gT", "[", 
             RowBox[{"[", "1", "]"}], "]"}], "*", 
            RowBox[{"w", "/", 
             RowBox[{"(", 
              RowBox[{"12", "*", "f1", "*", "nMol"}], ")"}]}]}], ")"}], "^", 
          "2"}]}], ")"}]}]}], ";", 
     RowBox[{"Dw", "=", 
      RowBox[{"2", "*", 
       RowBox[{"Pi", "/", "len"}]}]}], ";", 
     RowBox[{
      RowBox[{"gs", "[", "w_", "]"}], ":=", 
      RowBox[{
       RowBox[{"gT", "[", 
        RowBox[{"[", 
         RowBox[{"1", "+", 
          RowBox[{"w", "/", "Dw"}]}], "]"}], "]"}], "-", 
       RowBox[{"gg", "[", "w", "]"}]}]}], ";", 
     RowBox[{"Sig", "=", " ", 
      RowBox[{"k", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"1", "/", "rho"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"4", "*", "Pi", "*", "m", "*", 
              RowBox[{"3", "/", "2"}], "*", "nMol", "*", "k", "*", 
              RowBox[{"T", "/", 
               RowBox[{"(", 
                RowBox[{"3", "*", "nMol", "*", 
                 RowBox[{"h", "^", "2"}]}], ")"}]}]}], ")"}], "^", 
            RowBox[{"(", 
             RowBox[{"3", "/", "2"}], ")"}]}]}], "]"}], "+", 
         RowBox[{"5", "/", "2"}]}], ")"}]}]}], ";", 
     RowBox[{
      RowBox[{"z", "[", "y_", "]"}], ":=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "y", "+", 
         RowBox[{"y", "^", "2"}], "-", 
         RowBox[{"y", "^", "3"}]}], ")"}], "/", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", "y"}], ")"}], "^", "3"}]}]}], ";", 
     RowBox[{"fy", "=", 
      RowBox[{"f1", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"f1", "/", "De"}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"3", "/", "2"}], ")"}]}]}]}], ";", 
     RowBox[{"SR", "=", 
      RowBox[{
       RowBox[{"Log", "[", 
        RowBox[{
         RowBox[{"Pi", "^", 
          RowBox[{"(", 
           RowBox[{"1", "/", "2"}], ")"}]}], "*", 
         RowBox[{
          RowBox[{"E", " ", "^", 
           RowBox[{"(", 
            RowBox[{"3", "/", "2"}], ")"}]}], "/", "sym"}], "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"T", "^", "3"}], "/", "ThAThBThC"}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"1", "/", "2"}], ")"}]}]}], "]"}], "*", "k"}]}], ";", " ", 
     RowBox[{"Wgs", "=", 
      RowBox[{
       RowBox[{"1", "/", "3"}], "*", 
       RowBox[{"SR", "/", "k"}]}]}], ";", 
     RowBox[{"hb", "=", 
      RowBox[{
       RowBox[{"h", "/", "2"}], "/", "Pi"}]}], ";", 
     RowBox[{"beta", "=", 
      RowBox[{
       RowBox[{"1", "/", "k"}], "/", "T"}]}], ";", 
     RowBox[{
      RowBox[{"Who", "[", "w_", "]"}], ":=", 
      RowBox[{
       RowBox[{"beta", "*", "hb", "*", 
        RowBox[{"w", "/", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Exp", "[", 
            RowBox[{"beta", "*", "hb", "*", "w"}], "]"}], "-", "1"}], 
          ")"}]}]}], "-", 
       RowBox[{"Log", "[", 
        RowBox[{"1", "-", 
         RowBox[{"Exp", "[", 
          RowBox[{
           RowBox[{"-", "beta"}], "*", "hb", "*", "w"}], "]"}]}], "]"}]}]}], 
     ";", 
     RowBox[{"Svib", "=", 
      RowBox[{"k", "*", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"Dw", "/", "2"}], "/", "Pi"}], "*", 
          RowBox[{"gs", "[", "w", "]"}], "*", 
          RowBox[{"Who", "[", "w", "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"w", ",", "Dw", ",", 
           RowBox[{
            RowBox[{"1", "/", "dt"}], "*", "Pi"}], ",", "Dw"}], "}"}]}], 
        "]"}]}]}], ";", 
     RowBox[{"Sconf", "=", 
      RowBox[{"k", "*", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"Dw", "/", "2"}], "/", "Pi"}], "*", 
          RowBox[{"gg", "[", "w", "]"}], "*", "Wgs"}], ",", 
         RowBox[{"{", 
          RowBox[{"w", ",", "0", ",", 
           RowBox[{
            RowBox[{"1", "/", "dt"}], "*", "Pi"}], ",", "Dw"}], "}"}]}], 
        "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Svib", "+", "Sconf"}], ")"}], "/", "nMol"}], "*", 
      "1000"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.515828392936264*^9, 3.5158284174226646`*^9}, {
   3.5158284552898307`*^9, 3.51582848936878*^9}, {3.5158285216556263`*^9, 
   3.515828670628147*^9}, {3.51582874095817*^9, 3.5158288044908037`*^9}, 
   3.5158293248025637`*^9, {3.5158301956463733`*^9, 3.515830196159403*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"nMol", "=", "988"}], ";", 
  RowBox[{
  "dir", "=", 
   "\"\<c:\\\\Users\\\\rschulz\\\\Documents\\\\work\\\\lignin_entropie\\\\\
water\\\\\>\""}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"T", "=", "300"}], ",", 
    RowBox[{"T", "\[LessEqual]", "480"}], ",", 
    RowBox[{"T", "+=", "60"}], ",", " ", 
    RowBox[{
     RowBox[{"l", "=", 
      RowBox[{"Map", "[", 
       RowBox[{
        RowBox[{"Function", "[", 
         RowBox[{"i", ",", 
          RowBox[{
           RowBox[{"fn", "=", 
            RowBox[{
             RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
             RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
           RowBox[{"S", "[", 
            RowBox[{"T", ",", "nMol", ",", 
             RowBox[{
             "dir", "<>", "\"\<vac\>\"", "<>", "fn", "<>", "\"\<.cvs\>\""}], 
             ",", 
             RowBox[{"nMol", "/", 
              RowBox[{
               RowBox[{"Import", "[", 
                RowBox[{
                 RowBox[{
                 "dir", "<>", "\"\<vol\>\"", "<>", "fn", "<>", 
                  "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
               RowBox[{"[", 
                RowBox[{"1", ",", "2"}], "]"}], "]"}]}]}], "]"}]}]}], "]"}], 
        ",", 
        RowBox[{"Range", "[", 
         RowBox[{"0", ",", "80", ",", "20"}], "]"}]}], "]"}]}], ";", 
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
       RowBox[{"ToString", "[", 
        RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
       RowBox[{"ToString", "[", 
        RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}]}], 
   "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.5100860650552654`*^9, 3.5100861296789618`*^9}, {
   3.510086308410185*^9, 3.510086337945874*^9}, {3.510086700679621*^9, 
   3.5100867663233757`*^9}, {3.5100869777054663`*^9, 
   3.5100869797955856`*^9}, {3.5100884918200684`*^9, 3.510088524836957*^9}, {
   3.510163788836868*^9, 3.51016379044096*^9}, {3.5111047416108003`*^9, 
   3.5111048369142513`*^9}, {3.511104869903138*^9, 3.511104882993887*^9}, {
   3.5111049816835318`*^9, 3.511104994450262*^9}, {3.5111050730027547`*^9, 
   3.511105138736514*^9}, {3.511105173610509*^9, 3.511105204840296*^9}, {
   3.5111052617035475`*^9, 3.511105324942165*^9}, {3.5111053647984447`*^9, 
   3.5111053824564543`*^9}, {3.511105477014863*^9, 3.51110547958401*^9}, {
   3.511108320218485*^9, 3.5111083323521786`*^9}, {3.511110125514742*^9, 
   3.511110132412136*^9}, {3.5111156595412703`*^9, 3.5111156617213945`*^9}, {
   3.5151951005739574`*^9, 3.515195132420779*^9}, 3.5158291536947775`*^9}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 54.3785 0.187403\"\>"], "Print",
 CellChangeTimes->{3.5151951338648615`*^9, 3.515202596745714*^9, 
  3.5154823154367685`*^9, 3.515483202506506*^9, 3.515484294664974*^9, 
  3.51551124541437*^9, 3.515825581194442*^9, 3.5158279240314445`*^9, 
  3.51582914482627*^9, 3.515830256285842*^9}],

Cell[BoxData["\<\"360 63.2357 0.0884352\"\>"], "Print",
 CellChangeTimes->{3.5151951338648615`*^9, 3.515202596745714*^9, 
  3.5154823154367685`*^9, 3.515483202506506*^9, 3.515484294664974*^9, 
  3.51551124541437*^9, 3.515825581194442*^9, 3.5158279240314445`*^9, 
  3.51582914482627*^9, 3.515830257323901*^9}],

Cell[BoxData["\<\"420 70.4708 0.0793841\"\>"], "Print",
 CellChangeTimes->{3.5151951338648615`*^9, 3.515202596745714*^9, 
  3.5154823154367685`*^9, 3.515483202506506*^9, 3.515484294664974*^9, 
  3.51551124541437*^9, 3.515825581194442*^9, 3.5158279240314445`*^9, 
  3.51582914482627*^9, 3.515830258376961*^9}],

Cell[BoxData["\<\"480 77.781 0.171635\"\>"], "Print",
 CellChangeTimes->{3.5151951338648615`*^9, 3.515202596745714*^9, 
  3.5154823154367685`*^9, 3.515483202506506*^9, 3.515484294664974*^9, 
  3.51551124541437*^9, 3.515825581194442*^9, 3.5158279240314445`*^9, 
  3.51582914482627*^9, 3.515830259417021*^9}]
}, Open  ]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"new", " ", "data"}], "*)"}]], "Input",
 CellChangeTimes->{{3.515827940962413*^9, 3.5158279441075926`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"T", "=", "300"}], ",", 
    RowBox[{"T", "\[LessEqual]", "480"}], ",", 
    RowBox[{"T", "+=", "60"}], ",", " ", 
    RowBox[{
     RowBox[{"l", "=", 
      RowBox[{"Map", "[", 
       RowBox[{
        RowBox[{"Function", "[", 
         RowBox[{"i", ",", 
          RowBox[{
           RowBox[{"fn", "=", 
            RowBox[{
             RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
             RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
           RowBox[{"S", "[", 
            RowBox[{"T", ",", "nMol", ",", 
             RowBox[{
             "dir", "<>", "\"\<corr\>\"", "<>", "fn", "<>", "\"\<.dat\>\""}], 
             ",", 
             RowBox[{"nMol", "/", 
              RowBox[{
               RowBox[{"Import", "[", 
                RowBox[{
                 RowBox[{
                 "dir", "<>", "\"\<vol\>\"", "<>", "fn", "<>", 
                  "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
               RowBox[{"[", 
                RowBox[{"1", ",", "2"}], "]"}], "]"}]}]}], "]"}]}]}], "]"}], 
        ",", 
        RowBox[{"Range", "[", 
         RowBox[{"0", ",", "80", ",", "20"}], "]"}]}], "]"}]}], ";", 
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
       RowBox[{"ToString", "[", 
        RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
       RowBox[{"ToString", "[", 
        RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}]}], 
   "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.5158257676481066`*^9, 3.5158257753245454`*^9}, 
   3.5158274861714*^9}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 54.4175 0.190288\"\>"], "Print",
 CellChangeTimes->{3.5158279295567603`*^9, 3.5158292473651347`*^9}],

Cell[BoxData["\<\"360 63.2237 0.206644\"\>"], "Print",
 CellChangeTimes->{3.5158279295567603`*^9, 3.5158292482881875`*^9}],

Cell[BoxData["\<\"420 70.5805 0.181005\"\>"], "Print",
 CellChangeTimes->{3.5158279295567603`*^9, 3.515829249217241*^9}],

Cell[BoxData["\<\"480 77.8311 0.187842\"\>"], "Print",
 CellChangeTimes->{3.5158279295567603`*^9, 3.5158292501552944`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"T", "=", "300"}], ",", 
    RowBox[{"T", "\[LessEqual]", "480"}], ",", 
    RowBox[{"T", "+=", "60"}], ",", " ", 
    RowBox[{
     RowBox[{"l", "=", 
      RowBox[{"Map", "[", 
       RowBox[{
        RowBox[{"Function", "[", 
         RowBox[{"i", ",", 
          RowBox[{
           RowBox[{"fn", "=", 
            RowBox[{
             RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
             RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
           RowBox[{"Srot", "[", 
            RowBox[{"T", ",", "nMol", ",", 
             RowBox[{
             "dir", "<>", "\"\<rot\>\"", "<>", "fn", "<>", "\"\<.dat\>\""}], 
             ",", 
             RowBox[{"nMol", "/", 
              RowBox[{
               RowBox[{"Import", "[", 
                RowBox[{
                 RowBox[{
                 "dir", "<>", "\"\<vol\>\"", "<>", "fn", "<>", 
                  "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
               RowBox[{"[", 
                RowBox[{"1", ",", "2"}], "]"}], "]"}]}]}], "]"}]}]}], "]"}], 
        ",", 
        RowBox[{"Range", "[", 
         RowBox[{"0", ",", "80", ",", "20"}], "]"}]}], "]"}]}], ";", 
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
       RowBox[{"ToString", "[", 
        RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
       RowBox[{"ToString", "[", 
        RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}]}], 
   "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.515827968094965*^9, 3.5158279701900845`*^9}, {
   3.5158291157306056`*^9, 3.5158291174667053`*^9}, 3.515830263733268*^9}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 13.1723 0.0280742\"\>"], "Print",
 CellChangeTimes->{{3.5158302259141045`*^9, 3.5158302649503374`*^9}}],

Cell[BoxData["\<\"360 18.3848 0.0397038\"\>"], "Print",
 CellChangeTimes->{{3.5158302259141045`*^9, 3.5158302658633895`*^9}}],

Cell[BoxData["\<\"420 23.4987 0.106856\"\>"], "Print",
 CellChangeTimes->{{3.5158302259141045`*^9, 3.5158302667884426`*^9}}],

Cell[BoxData["\<\"480 28.6772 0.134925\"\>"], "Print",
 CellChangeTimes->{{3.5158302259141045`*^9, 3.515830267722496*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"nMol", "=", "988"}], ";", 
  RowBox[{
  "dir", "=", 
   "\"\<c:\\\\Users\\\\rschulz\\\\Documents\\\\work\\\\lignin_entropie\\\\\
water\\\\\>\""}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Map", "[", 
  RowBox[{
   RowBox[{"Function", "[", 
    RowBox[{"T", ",", 
     RowBox[{"T", "\[Rule]", 
      RowBox[{
       RowBox[{"densityF", "[", 
        RowBox[{"[", "T", "]"}], "]"}], "*", 
       RowBox[{"nMol", "/", 
        RowBox[{"Mean", "[", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{"Function", "[", 
            RowBox[{"i", ",", 
             RowBox[{
              RowBox[{"fn", "=", 
               RowBox[{
                RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
                RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
              RowBox[{
               RowBox[{"Import", "[", 
                RowBox[{
                 RowBox[{
                 "dir", "<>", "\"\<vol\>\"", "<>", "fn", "<>", 
                  "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
               RowBox[{"[", 
                RowBox[{"1", ",", "2"}], "]"}], "]"}]}]}], "]"}], ",", 
           RowBox[{"Range", "[", 
            RowBox[{"0", ",", "80", ",", "20"}], "]"}]}], "]"}], 
         "]"}]}]}]}]}], "]"}], ",", 
   RowBox[{"Range", "[", 
    RowBox[{"300", ",", "480", ",", "60"}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"density", "=", 
   RowBox[{"SparseArray", "[", "%", "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5111087168991737`*^9, 3.511108971022709*^9}, {
   3.511109008829871*^9, 3.511109091014572*^9}, {3.5111091532001286`*^9, 
   3.5111091554012547`*^9}, {3.515193979853856*^9, 3.515193984920146*^9}, {
   3.5158778114044933`*^9, 3.51587782061002*^9}, {3.5158779919258184`*^9, 
   3.5158780035694847`*^9}, 3.515878036541371*^9, {3.5158780795958333`*^9, 
   3.5158780811919246`*^9}, {3.5158781185840635`*^9, 3.515878122653296*^9}, {
   3.5158803247432485`*^9, 3.5158803261403284`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"300", "\[Rule]", "32.043038224355534`"}], ",", 
   RowBox[{"360", "\[Rule]", "29.460475229072838`"}], ",", 
   RowBox[{"420", "\[Rule]", "26.018199487587246`"}], ",", 
   RowBox[{"480", "\[Rule]", "21.75889809729546`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.5111087884612665`*^9, {3.5111088316907396`*^9, 3.511108971874758*^9}, {
   3.511109015872274*^9, 3.511109044593917*^9}, {3.5111090845912046`*^9, 
   3.5111090915856047`*^9}, 3.5111091568863397`*^9, 3.5111158321311417`*^9, 
   3.51114881082753*^9, 3.5151939857051907`*^9, 3.515482319355993*^9, 
   3.5154832064637327`*^9, 3.515511248762562*^9, 3.5158255846356387`*^9, 
   3.5158280211930017`*^9, 3.5158303255208015`*^9, 3.515877815924752*^9, 
   3.5158780822909875`*^9, 3.5158782810143538`*^9, {3.515880319949974*^9, 
   3.515880326775365*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{"old", " ", "data"}], " ", "-", " ", 
    RowBox[{"difference", " ", "to", " ", "bulk", " ", "water"}]}], 
   ")"}]}]], "Input",
 CellChangeTimes->{{3.515827530983963*^9, 3.5158275431866612`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
  "dir", "=", 
   "\"\<c:\\\\Users\\\\rschulz\\\\Documents\\\\work\\\\lignin_entropie\\\\\>\"\
"}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"T", "=", "300"}], ",", 
    RowBox[{"T", "\[LessEqual]", "300"}], ",", 
    RowBox[{"T", "+=", "60"}], ",", " ", 
    RowBox[{"Do", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"l", "=", 
        RowBox[{"Map", "[", 
         RowBox[{
          RowBox[{"Function", "[", 
           RowBox[{"i", ",", 
            RowBox[{
             RowBox[{"fn", "=", 
              RowBox[{
               RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
               RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
             "\[IndentingNewLine]", 
             RowBox[{"nMol", "=", 
              RowBox[{
               RowBox[{
                RowBox[{"Import", "[", 
                 RowBox[{
                  RowBox[{
                  "dir", "<>", "\"\<size\>\"", "<>", "fn", "<>", 
                   "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
                RowBox[{"[", 
                 RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", "3"}]}], ";", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"54.3785", "-", 
                RowBox[{"S", "[", 
                 RowBox[{"T", ",", "nMol", ",", 
                  RowBox[{
                  "dir", "<>", "\"\<vac\>\"", "<>", "fn", "<>", 
                   "\"\<.cvs\>\""}], ",", 
                  RowBox[{"density", "[", 
                   RowBox[{"[", "T", "]"}], "]"}]}], "]"}]}], ")"}], "*", 
              "nMol"}]}]}], "]"}], ",", 
          RowBox[{"Range", "[", 
           RowBox[{"p", ",", 
            RowBox[{"p", "+", "2"}]}], "]"}]}], "]"}]}], ";", 
       RowBox[{"Print", "[", 
        RowBox[{
         RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
         RowBox[{"ToString", "[", 
          RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
         RowBox[{"ToString", "[", 
          RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}], ",", 
      RowBox[{"{", 
       RowBox[{"p", ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], "]"}], ";", "  ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "p", " ", "is", " ", "the", " ", "starting", " ", "number", " ", "of", 
     " ", "the", " ", "files", " ", "for", " ", "one", " ", "temp"}], " ", 
    "-", " ", 
    RowBox[{"to", " ", "average", " ", "seperate", " ", 
     RowBox[{"e", ".", "g", ".", " ", "over"}], " ", 
     RowBox[{"collapsed", "/", "extended"}]}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.511108350643225*^9, 3.511108411413701*^9}, {
   3.511109190600268*^9, 3.51110919867673*^9}, {3.5111093684484406`*^9, 
   3.511109368515444*^9}, {3.511109398921183*^9, 3.5111094857161474`*^9}, {
   3.5111096196878104`*^9, 3.511109620449854*^9}, {3.511109651801647*^9, 
   3.511109663705328*^9}, {3.5111097512113333`*^9, 3.511109822220394*^9}, {
   3.5111098715782175`*^9, 3.5111099093943806`*^9}, {3.5111099873388386`*^9, 
   3.5111100431580315`*^9}, {3.511115683021613*^9, 3.5111156838856626`*^9}, {
   3.515193900878339*^9, 3.5151939009623437`*^9}, {3.5151939589996634`*^9, 
   3.5151939612747936`*^9}, {3.515194063367633*^9, 3.5151940638536606`*^9}, {
   3.515194361634693*^9, 3.515194484552723*^9}, 3.5151948439332786`*^9, 
   3.51519488312652*^9, 3.5151951412722855`*^9, {3.5151951754892426`*^9, 
   3.5151952460802803`*^9}, {3.5155111918173046`*^9, 3.515511230166498*^9}, {
   3.515828002714945*^9, 3.5158280105083904`*^9}, 3.515830304534601*^9}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 1170.82 52.0112\"\>"], "Print",
 CellChangeTimes->{3.515193988710363*^9, 3.5151943886812396`*^9, 
  3.515194485859798*^9, 3.5151951424303513`*^9, 3.51519517806539*^9, 
  3.515195251024563*^9, 3.515202631184684*^9, 3.515482320044032*^9, 
  3.515483207900815*^9, 3.5154843047435503`*^9, 3.5155112187068424`*^9, 
  3.515511249365596*^9, 3.515825585285676*^9, 3.5158280345997686`*^9, 
  3.5158292800060015`*^9, 3.515830368649268*^9, 3.515877830382579*^9, 
  3.515880343903344*^9}],

Cell[BoxData["\<\"300 1206.59 63.289\"\>"], "Print",
 CellChangeTimes->{3.515193988710363*^9, 3.5151943886812396`*^9, 
  3.515194485859798*^9, 3.5151951424303513`*^9, 3.51519517806539*^9, 
  3.515195251024563*^9, 3.515202631184684*^9, 3.515482320044032*^9, 
  3.515483207900815*^9, 3.5154843047435503`*^9, 3.5155112187068424`*^9, 
  3.515511249365596*^9, 3.515825585285676*^9, 3.5158280345997686`*^9, 
  3.5158292800060015`*^9, 3.515830368649268*^9, 3.515877830382579*^9, 
  3.5158803445143795`*^9}]
}, Open  ]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{"new", " ", "data"}], " ", "-", " ", 
    RowBox[{"difference", " ", "to", " ", "bulk", " ", "water"}]}], 
   ")"}]}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"T", "=", "300"}], ",", 
     RowBox[{"T", "\[LessEqual]", "300"}], ",", 
     RowBox[{"T", "+=", "60"}], ",", " ", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"l", "=", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{"Function", "[", 
            RowBox[{"i", ",", 
             RowBox[{
              RowBox[{"fn", "=", 
               RowBox[{
                RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
                RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"nMol", "=", 
               RowBox[{
                RowBox[{
                 RowBox[{"Import", "[", 
                  RowBox[{
                   RowBox[{
                   "dir", "<>", "\"\<size\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
                 RowBox[{"[", 
                  RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", "3"}]}], ";", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"54.4175", "-", 
                 RowBox[{"S", "[", 
                  RowBox[{"T", ",", "nMol", ",", 
                   RowBox[{
                   "dir", "<>", "\"\<corr\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", 
                   RowBox[{"density", "[", 
                    RowBox[{"[", "T", "]"}], "]"}]}], "]"}]}], ")"}], "*", 
               "nMol"}]}]}], "]"}], ",", 
           RowBox[{"Range", "[", 
            RowBox[{"p", ",", 
             RowBox[{"p", "+", "2"}]}], "]"}]}], "]"}]}], ";", 
        RowBox[{"Print", "[", 
         RowBox[{
          RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"p", ",", 
         RowBox[{"{", 
          RowBox[{"1", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], "]"}], ";"}], 
  "  "}]], "Input",
 CellChangeTimes->{{3.515825564950513*^9, 3.515825568685726*^9}, 
   3.5158280139535875`*^9}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 1109.84 61.0493\"\>"], "Print",
 CellChangeTimes->{3.5158280461584296`*^9, 3.5158778359859*^9}],

Cell[BoxData["\<\"300 1349.98 104.095\"\>"], "Print",
 CellChangeTimes->{3.5158280461584296`*^9, 3.515877836538931*^9}]
}, Open  ]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{"old", " ", "data"}], " ", "-", " ", 
    RowBox[{"absolute", " ", "value"}]}], ")"}]}]], "Input",
 CellChangeTimes->{{3.5158275854830804`*^9, 3.5158275987208376`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"T", "=", "300"}], ",", 
     RowBox[{"T", "\[LessEqual]", "300"}], ",", 
     RowBox[{"T", "+=", "60"}], ",", " ", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"l", "=", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{"Function", "[", 
            RowBox[{"i", ",", 
             RowBox[{
              RowBox[{"fn", "=", 
               RowBox[{
                RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
                RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"nMol", "=", 
               RowBox[{
                RowBox[{
                 RowBox[{"Import", "[", 
                  RowBox[{
                   RowBox[{
                   "dir", "<>", "\"\<size\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
                 RowBox[{"[", 
                  RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", "3"}]}], ";", 
              RowBox[{"(", 
               RowBox[{"S", "[", 
                RowBox[{"T", ",", "nMol", ",", 
                 RowBox[{
                 "dir", "<>", "\"\<vac\>\"", "<>", "fn", "<>", 
                  "\"\<.cvs\>\""}], ",", 
                 RowBox[{"density", "[", 
                  RowBox[{"[", "T", "]"}], "]"}]}], "]"}], ")"}]}]}], "]"}], 
           ",", 
           RowBox[{"Range", "[", 
            RowBox[{"p", ",", 
             RowBox[{"p", "+", "2"}]}], "]"}]}], "]"}]}], ";", 
        RowBox[{"Print", "[", 
         RowBox[{
          RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"p", ",", 
         RowBox[{"{", 
          RowBox[{"1", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], "]"}], ";"}], 
  "  "}]], "Input",
 CellChangeTimes->{{3.5151952745329075`*^9, 3.5151952957561216`*^9}, {
  3.5151954875920935`*^9, 3.515195487908112*^9}, {3.5155116816063185`*^9, 
  3.5155116878176737`*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 52.5074 0.104896\"\>"], "Print",
 CellChangeTimes->{3.5151952622042027`*^9, 3.515195298672288*^9, 
  3.515482321510116*^9, 3.515483208587854*^9, 3.5155112506186676`*^9, 
  3.5155116892377553`*^9, 3.5158255876538115`*^9, 3.5158778475905633`*^9}],

Cell[BoxData["\<\"300 52.7171 0.0619765\"\>"], "Print",
 CellChangeTimes->{3.5151952622042027`*^9, 3.515195298672288*^9, 
  3.515482321510116*^9, 3.515483208587854*^9, 3.5155112506186676`*^9, 
  3.5155116892377553`*^9, 3.5158255876538115`*^9, 3.5158778481995983`*^9}]
}, Open  ]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{"new", " ", "data"}], " ", "-", " ", 
    RowBox[{"absolute", " ", "value"}]}], ")"}]}]], "Input",
 CellChangeTimes->{{3.5158276060562572`*^9, 3.5158276062352676`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"T", "=", "300"}], ",", 
     RowBox[{"T", "\[LessEqual]", "300"}], ",", 
     RowBox[{"T", "+=", "60"}], ",", " ", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"l", "=", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{"Function", "[", 
            RowBox[{"i", ",", 
             RowBox[{
              RowBox[{"fn", "=", 
               RowBox[{
                RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
                RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"nMol", "=", 
               RowBox[{
                RowBox[{
                 RowBox[{"Import", "[", 
                  RowBox[{
                   RowBox[{
                   "dir", "<>", "\"\<size\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
                 RowBox[{"[", 
                  RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", "3"}]}], ";", 
              RowBox[{"(", 
               RowBox[{"S", "[", 
                RowBox[{"T", ",", "nMol", ",", 
                 RowBox[{
                 "dir", "<>", "\"\<corr\>\"", "<>", "fn", "<>", 
                  "\"\<.dat\>\""}], ",", 
                 RowBox[{"density", "[", 
                  RowBox[{"[", "T", "]"}], "]"}]}], "]"}], ")"}]}]}], "]"}], 
           ",", 
           RowBox[{"Range", "[", 
            RowBox[{"p", ",", 
             RowBox[{"p", "+", "2"}]}], "]"}]}], "]"}]}], ";", 
        RowBox[{"Print", "[", 
         RowBox[{
          RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"p", ",", 
         RowBox[{"{", 
          RowBox[{"1", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], "]"}], ";"}], 
  "  "}]], "Input",
 CellChangeTimes->{{3.5158257370593567`*^9, 3.5158257417316236`*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 52.6447 0.0926025\"\>"], "Print",
 CellChangeTimes->{3.515830373901569*^9, 3.515877850722742*^9}],

Cell[BoxData["\<\"300 52.5484 0.266077\"\>"], "Print",
 CellChangeTimes->{3.515830373901569*^9, 3.515877851270774*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"rotational", " ", "-", " ", 
    RowBox[{"missing", " ", "correct", " ", "rotational", " ", "momentum"}]}],
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"T", "=", "300"}], ",", 
      RowBox[{"T", "\[LessEqual]", "300"}], ",", 
      RowBox[{"T", "+=", "60"}], ",", " ", 
      RowBox[{"Do", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"l", "=", 
          RowBox[{"Map", "[", 
           RowBox[{
            RowBox[{"Function", "[", 
             RowBox[{"i", ",", 
              RowBox[{
               RowBox[{"fn", "=", 
                RowBox[{
                 RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
                 RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
               "\[IndentingNewLine]", 
               RowBox[{"nMol", "=", 
                RowBox[{
                 RowBox[{
                  RowBox[{"Import", "[", 
                   RowBox[{
                    RowBox[{
                    "dir", "<>", "\"\<size\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
                  RowBox[{"[", 
                   RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", "3"}]}], ";", 
               RowBox[{"(", 
                RowBox[{"Srot", "[", 
                 RowBox[{"T", ",", "nMol", ",", 
                  RowBox[{
                  "dir", "<>", "\"\<rot\>\"", "<>", "fn", "<>", 
                   "\"\<.dat\>\""}], ",", 
                  RowBox[{"density", "[", 
                   RowBox[{"[", "T", "]"}], "]"}]}], "]"}], ")"}]}]}], "]"}], 
            ",", 
            RowBox[{"Range", "[", 
             RowBox[{"p", ",", 
              RowBox[{"p", "+", "2"}]}], "]"}]}], "]"}]}], ";", 
         RowBox[{"Print", "[", 
          RowBox[{
           RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
           RowBox[{"ToString", "[", 
            RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
           RowBox[{"ToString", "[", 
            RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}], 
        ",", 
        RowBox[{"{", 
         RowBox[{"p", ",", 
          RowBox[{"{", 
           RowBox[{"1", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], "]"}], ";"}], 
   "  ", "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.51582760942645*^9, 3.5158276692048693`*^9}, {
  3.5158302926979246`*^9, 3.515830292990941*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 13.1695 0.0291565\"\>"], "Print",
 CellChangeTimes->{3.5158803475995555`*^9}],

Cell[BoxData["\<\"300 13.0844 0.0214523\"\>"], "Print",
 CellChangeTimes->{3.515880348148587*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"T", "=", "300"}], ",", 
     RowBox[{"T", "\[LessEqual]", "300"}], ",", 
     RowBox[{"T", "+=", "60"}], ",", " ", 
     RowBox[{"Do", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"l", "=", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{"Function", "[", 
            RowBox[{"i", ",", 
             RowBox[{
              RowBox[{"fn", "=", 
               RowBox[{
                RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\<_\>\"", "<>", 
                RowBox[{"ToString", "[", "i", "]"}]}]}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"nMol", "=", 
               RowBox[{
                RowBox[{
                 RowBox[{"Import", "[", 
                  RowBox[{
                   RowBox[{
                   "dir", "<>", "\"\<size\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", "\"\<Table\>\""}], "]"}], "[", 
                 RowBox[{"[", 
                  RowBox[{"1", ",", "2"}], "]"}], "]"}], "/", "3"}]}], ";", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"13.1723", "-", 
                 RowBox[{"Srot", "[", 
                  RowBox[{"T", ",", "nMol", ",", 
                   RowBox[{
                   "dir", "<>", "\"\<rot\>\"", "<>", "fn", "<>", 
                    "\"\<.dat\>\""}], ",", 
                   RowBox[{"density", "[", 
                    RowBox[{"[", "T", "]"}], "]"}]}], "]"}]}], ")"}], "*", 
               "nMol"}]}]}], "]"}], ",", 
           RowBox[{"Range", "[", 
            RowBox[{"p", ",", 
             RowBox[{"p", "+", "2"}]}], "]"}]}], "]"}]}], ";", 
        RowBox[{"Print", "[", 
         RowBox[{
          RowBox[{"ToString", "[", "T", "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"Mean", "[", "l", "]"}], "]"}], "<>", "\"\< \>\"", " ", "<>", 
          RowBox[{"ToString", "[", 
           RowBox[{"StandardDeviation", "[", "l", "]"}], "]"}]}], "]"}]}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"p", ",", 
         RowBox[{"{", 
          RowBox[{"1", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], "]"}], ";"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.515828128268126*^9, 3.515828133791442*^9}, {
  3.5158304210142636`*^9, 3.515830426303566*^9}}],

Cell[CellGroupData[{

Cell[BoxData["\<\"300 1.44652 18.1918\"\>"], "Print",
 CellChangeTimes->{
  3.5158281348025*^9, {3.5158304277416487`*^9, 3.515830454261165*^9}, 
   3.51587785802216*^9, 3.5158803628424273`*^9}],

Cell[BoxData["\<\"300 63.1051 10.7848\"\>"], "Print",
 CellChangeTimes->{
  3.5158281348025*^9, {3.5158304277416487`*^9, 3.515830454261165*^9}, 
   3.51587785802216*^9, 3.5158803633794584`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"S", "[", 
  RowBox[{"300", ",", "637", ",", 
   RowBox[{"dir", "<>", "\"\<vac300_1.cvs\>\""}], ",", 
   RowBox[{"density", "[", 
    RowBox[{"[", "300", "]"}], "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.515193715820754*^9, 3.5151937197889814`*^9}, {
  3.5151938443341045`*^9, 3.5151938449371395`*^9}, {3.515194492834197*^9, 
  3.5151945828163433`*^9}, {3.5151946539574127`*^9, 3.5151946540394173`*^9}, {
  3.5151949005605173`*^9, 3.51519490286965*^9}, {3.5151974650992007`*^9, 
  3.515197465514224*^9}}],

Cell[BoxData[
 InterpretationBox[
  TagBox[
   PanelBox[GridBox[{
      {
       StyleBox[
        StyleBox[
         DynamicBox[ToBoxes[
           FEPrivate`FrontEndResource["FEStrings", "sizeExplanation"], 
           StandardForm],
          ImageSizeCache->{379., {3., 11.}}],
         StripOnInput->False,
         DynamicUpdating->True], "Panel",
        StripOnInput->False,
        Background->None]},
      {
       ItemBox[
        TagBox[
         RowBox[{
          FractionBox["1000", "637"], " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "0.23946363772105395`"}], "+", 
            RowBox[{"0.0083145107`", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"0.44422253881166995`", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "0.003207241080466569`"}], "+", 
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{"0.00320724425391222`", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"Fourier", "[", 
                    RowBox[{
                    RowBox[{"Symbol", "[", "]"}], ",", 
                    RowBox[{"FourierParameters", "\[Rule]", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", 
                    RowBox[{"-", "1"}]}], "}"}]}]}], "]"}], "+", 
                    RowBox[{"Fourier", "[", 
                    RowBox[{
                    RowBox[{"Symbol", "[", "]"}], ",", 
                    RowBox[{"FourierParameters", "\[Rule]", 
                    RowBox[{"{", 
                    RowBox[{"1", ",", "1"}], "}"}]}]}], "]"}], "-", 
                    RowBox[{"$Failed", "\[LeftDoubleBracket]", 
                    RowBox[{"1", ",", "2"}], "\[RightDoubleBracket]"}]}], 
                    ")"}]}], ")"}], "\[LeftDoubleBracket]", "3", 
                   "\[RightDoubleBracket]"}]}], ")"}]}], "+", 
               RowBox[{"\[LeftSkeleton]", "1249", "\[RightSkeleton]"}]}], 
              ")"}]}]}], ")"}]}],
         Short[#, 5]& ],
        Background->GrayLevel[1],
        BaseStyle->{Deployed -> False},
        Frame->True,
        FrameStyle->GrayLevel[0, 0.2],
        StripOnInput->False]},
      {
       RowBox[{
        ButtonBox[
         StyleBox[
          StyleBox[
           DynamicBox[ToBoxes[
             FEPrivate`FrontEndResource["FEStrings", "sizeShowLess"], 
             StandardForm],
            ImageSizeCache->{67., {1., 11.}}],
           StripOnInput->False,
           DynamicUpdating->True], "Panel",
          StripOnInput->False,
          Background->None],
         Appearance->Automatic,
         ButtonFunction:>OutputSizeLimit`ButtonFunction[
          Identity, 62, 23041219856762042107, 5/2],
         Enabled->True,
         Evaluator->Automatic,
         Method->"Queued"], "\[ThinSpace]", 
        ButtonBox[
         StyleBox[
          StyleBox[
           DynamicBox[ToBoxes[
             FEPrivate`FrontEndResource["FEStrings", "sizeShowMore"], 
             StandardForm],
            ImageSizeCache->{75., {1., 11.}}],
           StripOnInput->False,
           DynamicUpdating->True], "Panel",
          StripOnInput->False,
          Background->None],
         Appearance->Automatic,
         ButtonFunction:>OutputSizeLimit`ButtonFunction[
          Identity, 62, 23041219856762042107, 5 2],
         Enabled->True,
         Evaluator->Automatic,
         Method->"Queued"], "\[ThinSpace]", 
        ButtonBox[
         StyleBox[
          StyleBox[
           DynamicBox[ToBoxes[
             FEPrivate`FrontEndResource["FEStrings", "sizeShowAll"], 
             StandardForm],
            ImageSizeCache->{113., {3., 11.}}],
           StripOnInput->False,
           DynamicUpdating->True], "Panel",
          StripOnInput->False,
          Background->None],
         Appearance->Automatic,
         ButtonFunction:>OutputSizeLimit`ButtonFunction[
          Identity, 62, 23041219856762042107, Infinity],
         Enabled->True,
         Evaluator->Automatic,
         Method->"Queued"], "\[ThinSpace]", 
        ButtonBox[
         StyleBox[
          StyleBox[
           DynamicBox[ToBoxes[
             FEPrivate`FrontEndResource["FEStrings", "sizeChangeLimit"], 
             StandardForm],
            ImageSizeCache->{98., {1., 11.}}],
           StripOnInput->False,
           DynamicUpdating->True], "Panel",
          StripOnInput->False,
          Background->None],
         Appearance->Automatic,
         ButtonFunction:>FrontEndExecute[{
            FrontEnd`SetOptions[
            FrontEnd`$FrontEnd, 
             FrontEnd`PreferencesSettings -> {"Page" -> "Evaluation"}], 
            FrontEnd`FrontEndToken["PreferencesDialog"]}],
         Evaluator->None,
         Method->"Preemptive"]}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
       "RowsIndexed" -> {}},
     GridBoxDividers->{
      "Columns" -> {{False}}, "ColumnsIndexed" -> {}, "Rows" -> {{False}}, 
       "RowsIndexed" -> {}},
     GridBoxItemSize->{
      "Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
       "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[1.2]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    DefaultBaseStyle->{},
    FrameMargins->5],
   Deploy,
   DefaultBaseStyle->{Deployed -> True}],
  Out[62]]], "Output",
 CellChangeTimes->{{3.5151945334815216`*^9, 3.515194583907406*^9}, 
   3.515194663111936*^9, 3.5151949038227043`*^9, 3.5151974671173162`*^9, 
   3.515482322465171*^9, 3.5154832094569035`*^9, 3.515511251409713*^9, 
   3.51582558850286*^9, 3.515829192685007*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"S", "[", 
  RowBox[{"300", ",", "637", ",", 
   RowBox[{"dir", "<>", "\"\<corr300_1.dat\>\""}], ",", 
   RowBox[{"density", "[", 
    RowBox[{"[", "300", "]"}], "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5151974316742887`*^9, 3.515197434303439*^9}, {
  3.5151974690984297`*^9, 3.5151974694974527`*^9}, {3.515829286760388*^9, 
  3.515829288298476*^9}}],

Cell[BoxData["52.751866440645635`"], "Output",
 CellChangeTimes->{3.5151974352344923`*^9, 3.5151974700644846`*^9, 
  3.515482322675183*^9, 3.515483209634914*^9, 3.515511251585723*^9, 
  3.51582558867587*^9, 3.5158292888625083`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Srot", "[", 
  RowBox[{"300", ",", "637", ",", 
   RowBox[{"dir", "<>", "\"\<rot300_1.dat\>\""}], ",", 
   RowBox[{"density", "[", 
    RowBox[{"[", "300", "]"}], "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5158293023972826`*^9, 3.5158293354801745`*^9}}],

Cell[BoxData["13.209728528599326`"], "Output",
 CellChangeTimes->{{3.515829304140382*^9, 3.515829335940201*^9}, {
   3.5158298132645025`*^9, 3.5158298736019535`*^9}, 3.515829935185476*^9, 
   3.515830205780953*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "39255.6"}], "-", 
    RowBox[{"T", "*", 
     RowBox[{"-", "39255.6"}]}]}], ",", "T"}], "]"}]], "Input",
 CellChangeTimes->{{3.510088763752622*^9, 3.5100887863659153`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Plot", "::", "pllim"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Range specification \[NoBreak]\\!\\(T\\)\[NoBreak] is not \
of the form {x, xmin, xmax}. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/Plot/pllim\\\", ButtonNote -> \
\\\"Plot::pllim\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.510088787247966*^9, 3.5103673316704125`*^9, 
  3.5110601721738205`*^9, 3.515482322744187*^9, 3.5154832096519146`*^9, 
  3.515511251626725*^9, 3.515825588720872*^9}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "39255.6`"}], "-", 
    RowBox[{"T", " ", 
     RowBox[{"(", 
      RowBox[{"-", "39255.6`"}], ")"}]}]}], ",", "T"}], "]"}]], "Output",
 CellChangeTimes->{3.5103673316834135`*^9, 3.5110601721888213`*^9, 
  3.5154823227471867`*^9, 3.5154832096539154`*^9, 3.5155112516297255`*^9, 
  3.515825588723872*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "39255.6`"}], "-", 
    RowBox[{"T", " ", 
     RowBox[{"(", 
      RowBox[{"-", "39255.6`"}], ")"}]}]}], ",", "T", ",", "300", ",", 
   "480"}], "]"}]], "Input",
 CellChangeTimes->{{3.510088792051241*^9, 3.510088794421376*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Plot", "::", "nonopt"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Options expected (instead of \[NoBreak]\\!\\(480\\)\
\[NoBreak]) beyond position \[NoBreak]\\!\\(2\\)\[NoBreak] in \
\[NoBreak]\\!\\(Plot[\\(\\(\\(\\(\\(\\(-39255.6`\\)\\) - \\(\\(T\\\\ \
\\(\\((\\(-39255.6`\\))\\)\\)\\)\\)\\)\\), T, 300, 480\\)\\)]\\)\[NoBreak]. \
An option must be a rule or a list of rules. \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/nonopt\\\", ButtonNote -> \
\\\"Plot::nonopt\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.510088795097415*^9, 3.5103673317594175`*^9, 
  3.5110601723638315`*^9, 3.5154823228321915`*^9, 3.515483209672916*^9, 
  3.515511251669728*^9, 3.5158255887628746`*^9}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "39255.6`"}], "-", 
    RowBox[{"T", " ", 
     RowBox[{"(", 
      RowBox[{"-", "39255.6`"}], ")"}]}]}], ",", "T", ",", "300", ",", 
   "480"}], "]"}]], "Output",
 CellChangeTimes->{3.510367331762418*^9, 3.511060172366832*^9, 
  3.5154823228351917`*^9, 3.515483209676916*^9, 3.515511251672728*^9, 
  3.515825588765875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "39255.6"}], "+", 
        RowBox[{
         RowBox[{"3", "/", "4"}], "*", "nMol", "*", "k", "*", 
         RowBox[{"(", 
          RowBox[{"T", "-", "300"}], ")"}]}], "-", 
        RowBox[{"T", " ", "97.8589"}]}], ")"}], "/", "T"}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "36087.5"}], "+", 
        RowBox[{
         RowBox[{"3", "/", "4"}], "*", "nMol", "*", "k", "*", 
         RowBox[{"(", 
          RowBox[{"T", "-", "360"}], ")"}]}], "-", 
        RowBox[{"T", "*", "106.515"}]}], ")"}], "/", "T"}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "32818.5"}], "+", 
        RowBox[{
         RowBox[{"3", "/", "4"}], "*", "nMol", "*", "k", "*", 
         RowBox[{"(", 
          RowBox[{"T", "-", "420"}], ")"}]}], "-", 
        RowBox[{"T", "*", "113.627"}]}], ")"}], "/", "T"}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "29189.6"}], "+", 
        RowBox[{
         RowBox[{"3", "/", "4"}], "*", "nMol", "*", "k", "*", 
         RowBox[{"(", 
          RowBox[{"T", "-", "480"}], ")"}]}], "-", 
        RowBox[{"T", "*", "120.95"}]}], ")"}], "/", "T"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"T", ",", "300", ",", "480"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5100888143325152`*^9, 3.5100888159816093`*^9}, {
  3.510088879362235*^9, 3.5100889341373672`*^9}, {3.5100890044553895`*^9, 
  3.5100890149119873`*^9}, {3.510089056320356*^9, 3.5100890758884754`*^9}, {
  3.5100891129595957`*^9, 3.51008913243771*^9}, {3.5101349615400143`*^9, 
  3.51013498296924*^9}, {3.510135018190254*^9, 3.5101350518311787`*^9}, {
  3.5101485190869713`*^9, 3.51014854686656*^9}, {3.510148601147665*^9, 
  3.5101486347385864`*^9}, {3.5101486758809395`*^9, 3.5101486787281027`*^9}, {
  3.510148795579786*^9, 3.5101488119857244`*^9}, {3.5101488447745996`*^9, 
  3.5101488874780416`*^9}, {3.5101489305785074`*^9, 3.5101489385889654`*^9}, {
  3.5101642642800617`*^9, 3.5101642764757595`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVzndQkwccxnFGOW0UlOEAXCxFi2DdSrQ/gkLLqhADKK5y8GZqNCGBkLzW
lgZyKIcKpmjBo0JKPaUQTyoqwg8Z0uoFRBQh2rA8VoWAYkAD9O0fz33v89/j
Fi+MSrSysLAIp/Z/3yzutLZAEdzK7S2tUQvQu9Teo4ey4GqO5/tLAuQFfsOo
rxVBv5PwD2/KRmHlGVWdCGbEBZ8u5AjwY5N6euEjEawMOSdPPC9AOznT5KYT
gVoRvWBhpgC3Gp4MBb4SQWn5S7GYFKCqpOZpxqQIdosX5YsSBPjzsH6wZYsY
WJrO59v9BOgfx4mkCcUweWz14ooZPjrIzHGeJWIwe7Xu3NbCx/0hj8pne8VQ
cbTMa+wqH32CrViPlyUBsUfzj6WIj0KTODD2YBIMKiU3eAw+9rHetjWok0AW
apnc6MTHUcOLCo/mJDDlR552H+Zh4nJTF2O+BG4JV4abqnl4Kj0jLmqvBL5P
90kvVPOwrGCljEFKoMH5fnufkIcHOCNd7+5KwDX9rqEmiId5R9bnXRqXQNvr
tnFvNx4Gzi26POknBd9XcY2MKS6uOrZqN10ghYC0iJ7hVi5e7LXdeLJYChsn
Bw70lHHxpMvGk4ndUpg3os8szuRir6JxOsQlGfL6WnflcLj48OKF52MxydCi
bM5M2MPFZ9EJ+2S5ybAhyD9b7cbF1KASl7EnybA6+k/raksunhtc5beOlgJc
mWzWYODg5kpaDjMwBY77E4QeOWi3xXz5sDwFpMWKmNlrHGQmeg57V6bAV2P1
DZvSOOhbW3W2z5gCHmERjY4EB5Nazkq/85XBtup82yXBHHybFf9ZDU8G9Inx
GfpaDn6SXu+eXySDXNOJx722HJQ4rvl8WZcM9vcPtdca2fiT3umMq3Mq0Jo+
mGgv2NgQSbvyaV8q2CQ7u4RWsVEk3zqTn5UKLwaC9W2FbLzd0qr1bUyF23a8
Mk0GG+M7NGU/WMtBz9x8v1nARo3toLmeLoe8yv6O9fvZ6Grc9HYgSQ5+8QYm
n85Gs7JijU4rhzCfUJWDBxvFwSXaoiE5aCLqA3Tz2Miqalbu9FbAF7W6cst3
BHINs2ZtggLuZylP2+sJdPQtrHMuUECANn7Z0gYCvTcktzE6FfCtkafWlBL4
7HD15kOOJER1XCl3ySOwfIft9uhQEp7aXrsXnkbgVNNFpb2KhANFtQWLjhPY
nqZyf4AkkP5/aW/EEiiYGll/r5aE89ev+kdSthKydtx5SMIkX/r4QwyBPgfd
92nrSWhXrZsOoHxmQxWpaSIhLSH395cs6s/r0ZdZLST8+25gxRwmgbKtMdlH
ukjIcBUoOeEE2t18cCWum4SeGwk6O8rFbp6/xfaQcCI7bHlFGIG6+WNVUX0k
bJ/of2JJ2b03cyhogIT6Orr4lxAC/86u3utnJGH1aKe+OZjAozZekT5jJOzc
8mOalPL71LOH1o6TMMFf8OVyyiuIWLHHexKkrI4CXhCBIvp44ZJJynUu1TZ7
CZxzK/am0xQJDju6c2/uITB/Tc0d+48kKCZUp5iUGx3O6eaZSWCMyui/BhIY
pxrvmDtNgpP1w3VfUzZOx76xmSGhkDa4YpRBYLq4xmg1S0LEqZGlasqug17m
WcqF154u3UX5Pztxbus=
     "]]}, 
   {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJwVyXk41Akcx3GrbZNjnsTmIdsaaUflaCnHpu3bzKNze+hwZKKI3xxmjebA
zPg96bFqyrVIKHbtjqM8HZRyNDVfIaJEbAk9KZSjNK6QkZ3++Dzv5/V8qCH8
/WG6Ojo6e7X72oEVXYt0UADfpVio/xHx0Paq8eo3X33DiDIg5CGXsYteVyOA
zsmmLbZaq/mVcfJaARxrynl67TgPPzeen1/WIIDMaHqbMoKHFNmBaWqLACYl
qVY9bB66vHo0zOgRgN7tq1O0QB7Ki1Vtp2cEgAt26glPHmaNdA+1bhICRVZv
nWrKw81M9j59vhDsPA7Kz7wLx+USDdOmWAhOHrvPhCrD8eDuhtKFPiG43TTO
KkkJR7sduj7NliKoPusjnQ4NR/60kOEfIAID1u3Qx27h2O/zoaP+vAiyHlOD
BZRw/Pjq2a3VT0RQlUfvefyOi2E/TPfSDcVQr4qoaL3HxeOnTjP3e4oh+UYI
KziHi9fzfpTQSTHIG+yy+GIuHmKP9k5UiWELzXMUvLiYHWSfnTkuhuas+eYg
Oy4y9BQ5M45R8HCT2ZzfEi5aHbX61YMXBY7el27W9nEwvc/IKbIgCobc0ke2
1nIw0sIpMux1FLQLrsaW5nOwL/bB/G6LaPjz/egyk5McvJ+e9t+YXzTcfsZx
P3GEg+2+od6Sc9GQkZpb7rSVg9LtxRZjj6Lh8/oZs0IqB5OGrBzX6cfAGbKD
4arLwY2V+hkHGDEwm8L0ixhgI2WTJidQFgNm78suVj5k44EwmxHbyhjoTipf
6nqFjQ41ysR+dQx0ZorbHdPYKGpNjAp2kEBzmrNLhIiNH5JDvlVxJfBT3aBP
tT8b56IuvzZUSKCkId7+0FY2ik1oSy17JeCb4UrztmHjH92mcSvNpZCrUvyb
YMDG+n36F+a8pRBw+rpf/gQLBTKXL7nJUpA27Nxs08XC8tanZQ4PpLDKXe5a
f5+FIS8Kr59cJAPLuRPdLy+zsNBoSFPnIYPS8m2fvslg4Uq184dBkfYXV/lN
xbJQk3CL1lImgz09rJv8MBYKdxSXKYZl8Ft7+pZOLxb6KJ8k/GIbC456YXc6
3FjIebWgKQuNhQhqcHWeNQtNHPJrzfNiYe9Rka87hYW2G6I76F2xUBeaNBU1
Q2B74L2Nh01IkO4sMpf0E1jqbuTmu4cEfuBs3Po2Amcb0xOM5SRck0W+TVIS
+Dxebn0XSXj7LsLl4mUCebOj9tU1JORtjGfs0FqX7+NecZ+E7O/X+05cItAu
wNq7rI6Es0znpD1ax21QkoWNJJT8vmKzpohA25cfO5NbSTC5WKIIKiBQ4uKX
GtRLgkrskEv7m0DKlbsXmK9JMKcpAjr+IrCAalPk/4aEZw9zVp3UusVwTLm/
nwSdcv2KrjwCrfvODm8fJEH5dGRtSi6BTan3PB3VJNxx/tnsUw6BRxav2Wc3
RgIlLd5eofWkNPHw2nESnNd57fLWehXhL1w9ScKyAa+0kmwCBR7j+WYzJAwf
HAw5kkXgkhv+V0xnSaCypwoMtc6lqSqMP5NAPxb9vuo8gQ+WJ7UYaEi4NPg8
0VRrpnz8hd48CW1Bt/pqMglUz/sPLP5CQtG5xdv4Wp8SqtS6CyRs12lUWGq9
cmiNZkHrXjONQdM5Av8HAypudQ==
     "]]}, 
   {Hue[0.1421359549995791, 0.6, 0.6], LineBox[CompressedData["
1:eJwVzns0lAkYBnCpthrRRZcTJdelPSHpsnYd3h1FOCkkU6N08803l84wY5gx
85VLhk2yjbZNq7ZOMekk1IqVy0vjUiSbLhtaI9kyqCGbSyP77R/Pec7vnOeP
x+aQMCTK2MjIaDud/7tvWcdMIxSBq9WjovBlAnQqXGT3mvaarx58KlgqQJ6P
P1NTKwI7pU2QYYkA9cLyxPT7Ihhe0O18xVyAk03nphY2ikBTHFgysFCAZvLQ
MZtWEcyofVufYiLATd0tOp8uEaRlfL0MjQSYrq75M21cBAE/GD0lB/j4y0Bn
f9tGMYRw532YrOLj92wymCEUAzNKSarO8nGxzMC2V4uBsctL1yvg466AxuLp
XjEEtyenPmPyca2fcVjzylhoCuZ7yFfyUTgm9mHtjYW+jxtzuj7y8E3Y0NP6
c7GgXuCufdvCww/dz0vtHsdClaaFva+Ah1GrxrTM+RLoOnpdcCOFhzHKNHbI
VgkEDlQEeh7kYdHF1TImJYG7Hl0mzsDDPeR77cc/JHAroGrIbBUPz+93Pv/z
iAQivZWZEwYu+sy9mjPuGgfFDy6lD3Rw0fqAtZenIA7sm/N3rK/goqrXdH30
tThIOXDvUkAuF6Mt1kdH9cQBKb8clKTgYq+iYSrAIh5+q6hTHI3kYp3qzLPh
8Hhwi05WZQAX23cf2Sk7Gw9rJJGlz225mOCrthhuiYdB5k9h1+dw8VS/tes3
DCmEb+X6o47EDeWM7FAfKXCeRLT5tpFottGQs08uBa/oOzEP75IYGmU/4FQu
hekVqx3Uv5LoUluZ8UYvBTeF3r85mcTYtoy4gy4y6MvS/cgiSRzKPDSrhicD
jVCqfRVE4ue4gp75V2Xg57bYXbmZRIm547yVWhkUhi8tvG5F4onOJYmWKxLg
+Bbmqc45JNYHMy583pkA5AEL755hDorkm77kZiaA7ElyY3QHB39ve1Li0pAA
x7YnTQ5rOHjoZV5R0kw5bHn8IoRxi4N5pv0GjaccdlS4W83K4aCl3n3oXawc
tt0Kyqw8wUFDaqlja4kc9hzke1gIOSj2U5dc1clB3TzI9WJzMKzycep3Tgq4
MurfYO7LQW73tKHkiAJMAx955Llx0Nzl8v0VFxWgCap73mvFQad18U+ZHQpg
HTY5rWVwsH1f9YYIcwoKN/hLqXECiz1Mv90dSEHnurmJlf8QONGkSl2UTkE3
K0lV307gi5R02yqkYJ53jO2eOgIFE++dK2opsK2LGTShbSwM8yiro+DSlux7
1bUErt1ru7NEQ8Fpz81ce9qJ6yqpvCYKhgpHZ+hrCHR69eGvzDYKDPJ2p/Qq
AmWbwrP2a+l9Qva2ynICzW5WXWD3UNDSUFYqpH3Nxj6f9ZqCwZeDjna0W+cP
V4a8ocAzJM3yZBmBtr0ndb7vKGhVCL1Zdwl8mFW91VVPQe+xXNm/dwiMnO0Q
vHaYAvMFZzYX0B5NyIhYM0KB8vDpyQjaVgRLbDdKwau9V7I0twkUeY5cXj5O
gUpl8Sm7hMA5t1k3l0xQcCPUr8uPdq5jTdmiSQoYf0sbPxcT2LD4VKuJgYKo
kHc3D9Nmp4+8nDtFAd/LQb2ctn6K1Tf7C/1v7Eh+cxGBSnGN3niagurp/BvH
aVv2OximaSte6O640/4Pkh1zjQ==
     "]]}, 
   {Hue[0.37820393249936934`, 0.6, 0.6], LineBox[CompressedData["
1:eJwVz3s0lAkYBnC5bBFO0g1dFJaKMauas+TUewhdTrsukdxXfN83M/YMM2MY
4+tYoWklRYVqOl1cao+SEirxIkOJWG2bKNc6aGPIniEj++0fz3nO77/n2Rgp
8I3W1tLSOsjk//6w6q2OFgpBfRa8ePl8tLttYjXIuCHdNqMtj488931uT+uF
wJ7q7GIzVgmqU+SNQlj0nrqkPs/Hry0X5pc1CyG3MqonPZePxjI/9cZ2ITwI
aVxbdJqPnL4XY+69QmgT3rwwnsZHeUld54kZIbxLrJgujeNj3qee0Y4dInhj
pqtbtpePO4MpHwOBCDwUaR99rPi4XKoJti4RQVzG3sLjCzw8tL/57sKQCK73
mfae6eahvZe2f+taMSiiEs7+W8lDgVrkHhgkhrjvWq/9lMPDYf/Pr5ouiKFM
UbJgK+DhRN/rB1YvxZDqWcla9TMPo9ep+90M48HlpO2IiQMP4zJOBPt6xMM/
Q7n3tYx5WKbYIHWj40FLL9TcXsXFI9R4/5eH8WDtzjJL7uBifphD/vmpeGBf
Oa2W3+ei+5IbBTOOEijVeyFrP8dFywjLXa4xEmBnHk31SuBizpCRU2yhBLTF
nLwDIVyMNXeKjR6QgPUYrlbs5uJQsnJ+v3kCpFW+PBJhw8WGnLN/TR5OgM4M
y+JUfS52BUR5S88lwL0rbXkV4xQmeZaYT75IgJQwcsDhNYWnRi0dtxgkgumu
kutraijcXm2Q6+eeCMUr8m13FFJovENTECpLhJMcC2NOFoV+0daf7KoTQem0
KOChmEJWfU3msCoRhrXDlx8Io1DckSn5hSWFVm/Xm796Uvg5K1K3jieFlK1u
ZrGOFM5Jbg0Y3pCCQKJvFmFOYbyprf7afin4euw7pNKhMK1nRYqFWRLMNe8M
DJ0gscnH4OKcdxJQV23y5T0kCmWcb5ezkoArK3COVJJY0fFnOUuZBI3OmaOa
eyRGdheV/aYjgyL6y7zPFRKLjEY1T11lsFjy/NiRTBItVNs+j4hl4D0j0jWS
kqhJf2DbXi6DYrq36EQ0iSKvkvIbYzK4I04lLvuR6F/zMt3FLhm2/BC++DiQ
yO1b0JRHJcPKx9siTFkkmrKuNpopkuGN/SH19+tItGMnvHJ7mwxVDc/M1iwl
sSu0dnuIKQ2jyG5J+UrgXWejHwMO0GDpFFMrGyNwtiUn3UROQ+F0t154N4F/
H5dveoI0BF3Ldit+RmDM7LjDo3oaXLo4eRGMtQX+zlUNNKxvOqYyZ2wftMm7
/CkNw8sm7ma3EJjCrqGLWmjYe3QPV9pMoN27iTdZHTSo54OqDjYRKOUczg7r
p2HwjI5yFgk0Ln1yMXiAhgHou13BuHCjdXHgIA1zRgcvChi3G07W+A7TwHF8
l/KhjsBNQ7+PeY4wf943pHXWEvg8u9bDUUVD1Z1d5B81BIbr2fjYT9Kg/KiX
H814OikzZPMUDYpFgW2WjNcTgSKraRq2ml7yyHtMoNB16urqGRoKTM7EpD0i
cPG9wNIVszREt16r3834sm1dlclXGuQdLhZzDwlULj/VvlTD7JnS6o1jHCyf
6l4yT8OlzQ7uDoxV84Ef9L7RoD/4/s5INYEZojqV9gINy26t3FDI2GLURrPA
OMCwMzec8X/yg2rV
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{300., -228.},
  PlotRange->{{300, 480}, {-228.71089834641822`, -178.06995812363039`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{
  3.51008881651764*^9, {3.5100888848155465`*^9, 3.5100888979742994`*^9}, 
   3.510088935103423*^9, 3.510089015384015*^9, 3.510089076202493*^9, 
   3.5100891329937415`*^9, 3.5101349898666344`*^9, {3.5101350212324286`*^9, 
   3.5101350526022224`*^9}, {3.5101485277264657`*^9, 3.5101485473285866`*^9}, 
   3.5101486361786685`*^9, 3.5101486791371255`*^9, 3.5101488127087655`*^9, 
   3.5101488880110726`*^9, 3.5101489403560667`*^9, 3.510164277032791*^9, 
   3.510367332030433*^9, 3.5110601724038334`*^9, 3.5154823228951955`*^9, 
   3.515483209707918*^9, 3.5155112517137303`*^9, 3.515825588817878*^9}]
}, Open  ]]
},
WindowSize->{591, 749},
WindowMargins->{{Automatic, 308}, {Automatic, 20}},
FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 1279, 31, 112, "Input"],
Cell[1849, 55, 88, 1, 30, "Output"]
}, Open  ]],
Cell[1952, 59, 503, 12, 31, "Input"],
Cell[2458, 73, 92, 1, 31, InheritFromParent],
Cell[2553, 76, 464, 9, 52, "Input"],
Cell[3020, 87, 9416, 261, 532, "Input"],
Cell[12439, 350, 8594, 249, 532, "Input"],
Cell[21036, 601, 198, 6, 72, "Input"],
Cell[CellGroupData[{
Cell[21259, 611, 2551, 55, 152, "Input"],
Cell[CellGroupData[{
Cell[23835, 670, 307, 4, 23, "Print"],
Cell[24145, 676, 308, 4, 23, "Print"],
Cell[24456, 682, 308, 4, 23, "Print"],
Cell[24767, 688, 306, 4, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[25100, 696, 147, 3, 31, "Input"],
Cell[CellGroupData[{
Cell[25272, 703, 1694, 44, 152, "Input"],
Cell[CellGroupData[{
Cell[26991, 751, 122, 1, 23, "Print"],
Cell[27116, 754, 122, 1, 23, "Print"],
Cell[27241, 757, 120, 1, 23, "Print"],
Cell[27364, 760, 122, 1, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[27535, 767, 1746, 44, 152, "Input"],
Cell[CellGroupData[{
Cell[29306, 815, 125, 1, 23, "Print"],
Cell[29434, 818, 125, 1, 23, "Print"],
Cell[29562, 821, 124, 1, 23, "Print"],
Cell[29689, 824, 122, 1, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[29860, 831, 2022, 49, 232, "Input"],
Cell[31885, 882, 860, 15, 30, "Output"]
}, Open  ]],
Cell[32760, 900, 261, 7, 31, "Input"],
Cell[33024, 909, 151, 5, 31, "Input"],
Cell[CellGroupData[{
Cell[33200, 918, 3512, 77, 192, "Input"],
Cell[CellGroupData[{
Cell[36737, 999, 498, 7, 23, "Print"],
Cell[37238, 1008, 499, 7, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[37764, 1019, 193, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[37982, 1029, 2333, 60, 152, "Input"],
Cell[CellGroupData[{
Cell[40340, 1093, 117, 1, 23, "Print"],
Cell[40460, 1096, 119, 1, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[40606, 1101, 233, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[40864, 1111, 2312, 59, 152, "Input"],
Cell[CellGroupData[{
Cell[43201, 1174, 266, 3, 23, "Print"],
Cell[43470, 1179, 267, 3, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[43764, 1186, 233, 6, 31, "Input"],
Cell[CellGroupData[{
Cell[44022, 1196, 2209, 57, 152, "Input"],
Cell[CellGroupData[{
Cell[46256, 1257, 119, 1, 23, "Print"],
Cell[46378, 1260, 118, 1, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[46545, 1267, 2522, 63, 192, "Input"],
Cell[CellGroupData[{
Cell[49092, 1334, 99, 1, 23, "Print"],
Cell[49194, 1337, 97, 1, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[49340, 1344, 2357, 60, 152, "Input"],
Cell[CellGroupData[{
Cell[51722, 1408, 193, 3, 23, "Print"],
Cell[51918, 1413, 193, 3, 23, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[52160, 1422, 532, 10, 31, "Input"],
Cell[52695, 1434, 5897, 157, 206, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[58629, 1596, 380, 8, 31, "Input"],
Cell[59012, 1606, 231, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[59280, 1614, 282, 6, 31, "Input"],
Cell[59565, 1622, 214, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[59816, 1630, 250, 7, 31, "Input"],
Cell[60069, 1639, 606, 12, 21, "Message"],
Cell[60678, 1653, 380, 10, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[61095, 1668, 305, 9, 31, "Input"],
Cell[61403, 1679, 836, 15, 56, "Message"],
Cell[62242, 1696, 400, 11, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62679, 1712, 2143, 52, 112, "Input"],
Cell[64825, 1766, 7912, 138, 235, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)