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

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

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     24694,        729]
NotebookOptionsPosition[     23723,        691]
NotebookOutlinePosition[     24066,        706]
CellTagsIndexPosition[     24023,        703]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"Quit", "[", "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.6243788477743998`*^9, 3.6243788499814*^9}}],

Cell[CellGroupData[{

Cell[TextData[{
 "Para calcular funciones de partici\[OAcute]n de la red cuadrada. ",
 StyleBox["[Beale,  ",
  FontWeight->"Plain"],
 StyleBox["Exact Distribution of Energies in the Two-Dimensional Ising Model",
  FontWeight->"Plain",
  FontSlant->"Italic"],
 StyleBox[", PRL ",
  FontWeight->"Plain"],
 StyleBox["76",
  FontVariations->{"Underline"->True}],
 StyleBox[" 78]",
  FontWeight->"Plain"]
}], "Section",
 CellChangeTimes->{{3.5499764396844*^9, 3.5499764518792*^9}, {
  3.5499778603952*^9, 3.5499779322988*^9}, {3.6243633676266003`*^9, 
  3.6243633774046*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"n", " ", "=", " ", "8"}], ";"}], "  ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"dimensi\[OAcute]n", " ", "de", " ", "la", " ", 
     RowBox[{"red", ".", " ", "Para"}], " ", "una", " ", "red", " ", "de", 
     " ", "20", "x20", " ", "puede", " ", "tardar", " ", "unos", " ", "10", 
     " ", "segundos"}], ",", " ", 
    RowBox[{"dependiendo", " ", "de", " ", "la", " ", 
     RowBox[{"computadora", ".", " ", "Si"}], " ", "van", " ", "a", " ", 
     "trabajar", " ", "con", " ", "redes", " ", "muy", " ", "grandes"}], ",", 
    " ", 
    RowBox[{
    "convendr\[IAcute]a", " ", "que", " ", "guardaran", " ", "los", " ", 
     "resultados", " ", "en", " ", "un", " ", "archivo", " ", "para", " ", 
     "no", " ", "tener", " ", "que", " ", "calcularlos", " ", "cada", " ", 
     "vez"}]}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"m", "=", "n"}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"m", " ", "=", " ", 
    RowBox[{"n", " ", "para", " ", "la", " ", "red", " ", 
     RowBox[{"cuadrada", ".", " ", "Tambi\[EAcute]n"}], " ", "pueden", " ", 
     "calcularse", " ", "las", " ", "funciones", " ", "de", " ", 
     "partici\[OAcute]n", " ", "de", " ", "redes", " ", "rectangulares"}]}], 
   "*)"}], "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "a", ",", "b", ",", "c", ",", "s", ",", "c2", ",", "s2", ",", "x", ",", 
    "z", ",", "coef", ",", "z1", ",", "z2", ",", "z3", ",", "z4"}], "]"}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"acc", "=", 
   RowBox[{"Floor", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{"n", " ", "m", " ", 
       RowBox[{
        RowBox[{"Log", "[", "2", "]"}], "/", 
        RowBox[{"Log", "[", "10", "]"}]}]}], "]"}], 
     RowBox[{"(", " ", 
      RowBox[{"3", "/", "2"}], ")"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"Mod", "[", 
      RowBox[{"n", ",", "2"}], "]"}], "\[Equal]", "0"}], ",", 
    RowBox[{
     RowBox[{"z1", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"c2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", 
           RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z2", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"s2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", 
           RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z3", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"c", "[", "0", "]"}], " ", 
       RowBox[{"c", "[", "n", "]"}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"c2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "2", ",", 
           RowBox[{"n", "-", "2"}], ",", "2"}], "}"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z4", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"s", "[", "0", "]"}], " ", 
       RowBox[{"s", "[", "n", "]"}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"s2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "2", ",", 
           RowBox[{"n", "-", "2"}], ",", "2"}], "}"}]}], "]"}]}]}]}], ",", 
    RowBox[{
     RowBox[{"z1", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"c", "[", "n", "]"}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"c2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", 
           RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z2", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"s", "[", "n", "]"}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"s2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", 
           RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z3", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"c", "[", "0", "]"}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"c2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "2", ",", 
           RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"z4", "=", 
      RowBox[{
       RowBox[{"2", "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"m", " ", 
           RowBox[{"n", "/", "2"}]}], "-", "1"}], ")"}]}], " ", 
       RowBox[{"s", "[", "0", "]"}], " ", 
       RowBox[{"Product", "[", 
        RowBox[{
         RowBox[{"s2", "[", "k", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "2", ",", 
           RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}]}]}]}]}], "]"}],
   ";"}], "\n", 
 RowBox[{
  RowBox[{"z1", "+", "z2", "+", "z3", "+", "z4"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"c", "[", "0", "]"}], "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "x"}], ")"}], "^", "m"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x", 
         RowBox[{"(", 
          RowBox[{"1", "+", "x"}], ")"}]}], ")"}], "^", "m"}]}], ")"}], "/", 
    RowBox[{"2", "^", 
     RowBox[{"(", 
      RowBox[{"m", "/", "2"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"s", "[", "0", "]"}], "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "x"}], ")"}], "^", "m"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x", 
         RowBox[{"(", 
          RowBox[{"1", "+", "x"}], ")"}]}], ")"}], "^", "m"}]}], ")"}], "/", 
    RowBox[{"2", "^", 
     RowBox[{"(", 
      RowBox[{"m", "/", "2"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"c", "[", "n", "]"}], "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "x"}], ")"}], "^", "m"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x", 
         RowBox[{"(", 
          RowBox[{"1", "-", "x"}], ")"}]}], ")"}], "^", "m"}]}], ")"}], "/", 
    RowBox[{"2", "^", 
     RowBox[{"(", 
      RowBox[{"m", "/", "2"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"s", "[", "n", "]"}], "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "x"}], ")"}], "^", "m"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x", 
         RowBox[{"(", 
          RowBox[{"1", "-", "x"}], ")"}]}], ")"}], "^", "m"}]}], ")"}], "/", 
    RowBox[{"2", "^", 
     RowBox[{"(", 
      RowBox[{"m", "/", "2"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"b", "=", 
   RowBox[{"2", "x", 
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"x", "^", "2"}]}], ")"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"a", "[", "k_", "]"}], "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", 
       RowBox[{"x", "^", "2"}]}], ")"}], "^", "2"}], "-", 
    RowBox[{"b", " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"Pi", " ", 
       RowBox[{"k", "/", "n"}]}], "]"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"c2", "[", "k_", "]"}], "=", 
   RowBox[{
    RowBox[{"2", "^", 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       RowBox[{"2", "m"}]}], ")"}]}], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sum", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"m", "!"}], "/", 
           RowBox[{"j", "!"}]}], "/", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"m", "-", "j"}], ")"}], "!"}]}], 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"a", "[", "k", "]"}], "^", "2"}], "-", 
            RowBox[{"b", "^", "2"}]}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"j", "/", "2"}], ")"}]}], " ", 
         RowBox[{
          RowBox[{"a", "[", "k", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"m", "-", "j"}], ")"}]}]}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "m", ",", "2"}], "}"}]}], "]"}], "+", 
      RowBox[{"b", "^", "m"}]}], ")"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"s2", "[", "k_", "]"}], "=", 
   RowBox[{
    RowBox[{"2", "^", 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       RowBox[{"2", "m"}]}], ")"}]}], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sum", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"m", "!"}], "/", 
           RowBox[{"j", "!"}]}], "/", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"m", "-", "j"}], ")"}], "!"}]}], 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"a", "[", "k", "]"}], "^", "2"}], "-", 
            RowBox[{"b", "^", "2"}]}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"j", "/", "2"}], ")"}]}], " ", 
         RowBox[{
          RowBox[{"a", "[", "k", "]"}], "^", 
          RowBox[{"(", 
           RowBox[{"m", "-", "j"}], ")"}]}]}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "m", ",", "2"}], "}"}]}], "]"}], "-", 
      RowBox[{"b", "^", "m"}]}], ")"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"z", "=", 
   RowBox[{"Rationalize", "[", 
    RowBox[{"Chop", "[", 
     RowBox[{"Expand", "[", 
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"z1", "+", "z2", "+", "z3", "+", "z4"}], ",", "acc"}], "]"}], 
      "]"}], "]"}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"coef", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"m", " ", "n"}], "+", "1"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Do", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"coef", "[", 
       RowBox[{"[", 
        RowBox[{"j", "+", "1"}], "]"}], "]"}], "=", 
      RowBox[{"Coefficient", "[", 
       RowBox[{"z", ",", "x", ",", 
        RowBox[{"2", "j"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", 
       RowBox[{"m", " ", "n"}]}], "}"}]}], "]"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "el", " ", "resultado", " ", "se", " ", "guarda", " ", "en", " ", "la", 
    " ", "lista", " ", "coef"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5496260152834*^9, 3.5496261473912*^9}, {
   3.5496270953677998`*^9, 3.5496271182598*^9}, {3.5499764660352*^9, 
   3.5499765495746*^9}, {3.5499767751238003`*^9, 3.5499768627111998`*^9}, 
   3.6243593630892*^9, 3.6243598272594*^9, {3.624359877469*^9, 
   3.624359911534*^9}, {3.6243633829136*^9, 3.6243633895316*^9}, 
   3.6243634791594*^9}],

Cell[CellGroupData[{

Cell[BoxData["coef"], "Input",
 CellChangeTimes->{{3.6243593672441998`*^9, 3.6243593675912*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "2", ",", "0", ",", "128", ",", "256", ",", "4672", ",", "17920", ",", 
   "145408", ",", "712960", ",", "4274576", ",", "22128384", ",", "118551552",
    ",", "610683392", ",", "3150447680", ",", "16043381504", ",", 
   "80748258688", ",", "396915938304", ",", "1887270677624", ",", 
   "8582140066816", ",", "36967268348032", ",", "149536933509376", ",", 
   "564033837424064", ",", "1971511029384704", ",", "6350698012553216", ",", 
   "18752030727310592", ",", "50483110303426544", ",", "123229776338119424", 
   ",", "271209458049836032", ",", "535138987032308224", ",", 
   "941564975390477248", ",", "1469940812209435392", ",", 
   "2027486077172296064", ",", "2462494093546483712", ",", 
   "2627978003957146636", ",", "2462494093546483712", ",", 
   "2027486077172296064", ",", "1469940812209435392", ",", 
   "941564975390477248", ",", "535138987032308224", ",", "271209458049836032",
    ",", "123229776338119424", ",", "50483110303426544", ",", 
   "18752030727310592", ",", "6350698012553216", ",", "1971511029384704", ",",
    "564033837424064", ",", "149536933509376", ",", "36967268348032", ",", 
   "8582140066816", ",", "1887270677624", ",", "396915938304", ",", 
   "80748258688", ",", "16043381504", ",", "3150447680", ",", "610683392", 
   ",", "118551552", ",", "22128384", ",", "4274576", ",", "712960", ",", 
   "145408", ",", "17920", ",", "4672", ",", "256", ",", "128", ",", "0", ",",
    "2"}], "}"}]], "Output",
 CellChangeTimes->{
  3.6243593679392*^9, 3.6243598303234*^9, {3.6243634709863997`*^9, 
   3.6243634815604*^9}, 3.6243788661846*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Total", "[", "coef", "]"}], "\[IndentingNewLine]", 
 SuperscriptBox["2", 
  SuperscriptBox["n", "2"]]}], "Input",
 CellChangeTimes->{{3.6243634862204*^9, 3.6243634967194*^9}}],

Cell[BoxData["18446744073709551616"], "Output",
 CellChangeTimes->{{3.6243634885034*^9, 3.6243634971994*^9}, 
   3.6243788684505997`*^9}],

Cell[BoxData["18446744073709551616"], "Output",
 CellChangeTimes->{{3.6243634885034*^9, 3.6243634971994*^9}, 
   3.6243788684526*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
La funci\[OAcute]n de partici\[OAcute]n, la energ\[IAcute]a y el calor espec\
\[IAcute]fico:\
\>", "Subsection",
 CellChangeTimes->{{3.5499765681546*^9, 3.5499765798677998`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Z", "[", "x_", "]"}], ":=", " ", 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"coef", "[", 
       RowBox[{"[", 
        RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", 
      SuperscriptBox["x", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", 
          SuperscriptBox["n", "2"]}], "-", 
         RowBox[{"4", "i"}]}], ")"}]]}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"i", ",", " ", "0", ",", " ", 
       RowBox[{
        RowBox[{"Length", "[", "coef", "]"}], "-", "1"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"U", "[", "x_", "]"}], ":=", " ", 
   RowBox[{
    FractionBox["1", 
     RowBox[{
      SuperscriptBox["n", "2"], 
      RowBox[{"Z", "[", "x", "]"}]}]], " ", 
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", 
          SuperscriptBox["n", "2"]}], "-", 
         RowBox[{"4", "i"}]}], ")"}], " ", 
       RowBox[{"coef", "[", 
        RowBox[{"[", 
         RowBox[{"i", "+", "1"}], "]"}], "]"}], " ", 
       SuperscriptBox["x", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", 
           SuperscriptBox["n", "2"]}], "-", 
          RowBox[{"4", "i"}]}], ")"}]]}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"i", ",", " ", "0", ",", " ", 
        RowBox[{
         RowBox[{"Length", "[", "coef", "]"}], "-", "1"}]}], "}"}]}], 
     "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Cv", "[", "x_", "]"}], ":=", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"Log", "[", "x", "]"}], "2"], " ", 
    RowBox[{"(", "  ", 
     RowBox[{
      RowBox[{
       FractionBox["1", 
        RowBox[{
         SuperscriptBox["n", "2"], 
         RowBox[{"Z", "[", "x", "]"}]}]], " ", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", 
              SuperscriptBox["n", "2"]}], "-", 
             RowBox[{"4", "i"}]}], ")"}], "2"], " ", 
          RowBox[{"coef", "[", 
           RowBox[{"[", 
            RowBox[{"i", "+", "1"}], "]"}], "]"}], "  ", 
          SuperscriptBox["x", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", 
              SuperscriptBox["n", "2"]}], "-", 
             RowBox[{"4", "i"}]}], ")"}]]}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"i", ",", " ", "0", ",", " ", 
           RowBox[{
            RowBox[{"Length", "[", "coef", "]"}], "-", "1"}]}], "}"}]}], 
        "]"}]}], "-", 
      RowBox[{
       SuperscriptBox["n", "2"], " ", 
       SuperscriptBox[
        RowBox[{"U", "[", "x", "]"}], "2"]}]}], ")"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5499765869779997`*^9, 3.549976636307*^9}, {
  3.549976671808*^9, 3.5499767101776*^9}, {3.5499767448942003`*^9, 
  3.5499767594658003`*^9}, {3.6243788103278*^9, 3.6243788442514*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"T", ",", " ", 
      RowBox[{"Cv", "[", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"1", "/", "T"}]], "]"}]}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"T", ",", " ", "0.1", ",", " ", "5.", ",", " ", "0.05"}], 
     "}"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5499770694132*^9, 3.5499771610882*^9}, {
  3.6243788769546003`*^9, 3.6243788773845997`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000017`, 0.24, 0.6], PointBox[CompressedData["
1:eJw91QtMU2cUAOBbXqWlQEvBbjB1D5DhBoqI8gj9DwZQ48CpZMCcqEOBJdgK
ZhuOhGWbRTcmPiCj4zEdDLLRic4JBC0CXRnyKKKUQYEhlNcqrJTCkELV0fV0
J7m5+fOf850/N+fe+8r7wv3HrSiKSlq9TPeib01RTyiMrQGm+O3/tXm/Dden
w8z3blz/HGbO7yFjQXljmfSX+SnJpviDZDATY8KLTvLN9SpTvusZiYSv6DTF
ELHZW5dzNkTFNzsjpv13WOHPMH+UlOX6LXO73YnZHyNBeVl26ctbiLl+nPhE
HOTpy6KIud8k6T1bnfVCWRye6y/Smi5tXNmaTMyehlBHzosfvXkKvWmSWdue
m/04G70ZovXqcJs8fxY9LfnOIBbdXX8JPR2hmqzsu3WFmK8j/l9c0o1uvor+
HOnustO47anAej0RfJmW7RYnwX7zRPg4LVXl/wt6C0S7T7U2l9Sgt0Daq77J
qRmsR+8fsmPuDZ8Hjo3oLRJ+UHCUq1CG3hMyaLvmmq+kBb0lkpe/T/Uc2tBb
IpGhheK67Qr0DGQX48nOgM5u9JbJ2mG/c/u396C3QrKyFAHl473oGYlfyobF
03/2o2ckB2T1Ai/fIfSeEmOC71tJI4/Qe0YCPwwd7Kep0XtOpIfcA2f2jlvm
CMCDW8EbmMR8Cl4bcFenFWrQp2Bx8QfvXRkz5L/yABpUeFqnVkfOYj8aNN+P
LAqc1aFnBYsy+uBDSo+eFfytFWwMk+rRs4Lij0YT8p3m8TzW4KOs3Bml16Nn
DcH8VCV1Uo+eDbhd9i3+ungOPRuQFCQdP5SnQ88GLkxt2bxOrkXPFnhkMngp
fgY9WxAJWZv6VjTo2cFtcXjsdPMUenbwYmJZqaJ2Aj07SOQwDcbhMfTokPt9
gWBHiBo9OtgrZV2lLSPo2QMV0/FZjGEYPXtQK71p2ltD6K3uL1eFH701gB4D
PhaG1JVo+9FjwE8flIhyDvahxwSPT2M/YSz0oscEQ4v8IuumEj0mJI7yBt67
aJkPB3iakci5/NVD9BxgWtUl0pc8QI8FVNiZlni5Zb5YYHtsXXCQ8T56LMg8
1nBq6UIXeo6Qn33i3nqwzKcj6LMcdWmMTvScoDVC/CNP046eE8QpEmqkQ5b5
dgLnQNG70ep76DnDjGuHlLbUip4z0EO45dEvtaLHhohgb9++6N9xvthwhHHD
U5jbgj4bthlWSplKOdazwaOxQNKwQY792HDVMK+v/dzyXeTAtSiNIXRKhh4H
dqf3OlQekKHHgVov65LDrc3ocUB85bCzOrwZPQ5sdB8a/FXWhJ4LeJbTZpN3
N6HnArwq+RU/ieX9d4HX21ImehR30XOBJqN/Gm+5Ab3V+iom/eimBvS4cENU
6O51QorPhwvV2so9YzfvoMeFiZLkeC/qDnpcuD736prs2NvocUEgeDvv3HXL
f8EVUu0zt02x68m/lEaheg==
     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  Method->{},
  PlotRange->{{0, 5.}, {0, 1.1910584937062711`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{0.1, 0.1}, {0.023821169874125424`, 
   0.023821169874125424`}}]], "Output",
 CellChangeTimes->{
  3.5499770943171997`*^9, {3.5499771262872*^9, 3.5499771624472*^9}, 
   3.624363506799*^9, 3.6243788787446003`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"T", ",", " ", 
      RowBox[{"U", "[", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"1", "/", "T"}]], "]"}]}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"T", ",", " ", "0.1", ",", " ", "5.", ",", " ", "0.05"}], 
     "}"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{3.6243788809706*^9}],

Cell[BoxData[
 GraphicsBox[{{}, 
   {RGBColor[0.24720000000000017`, 0.24, 0.6], PointBox[CompressedData["
1:eJxF1A1MU3cQAPBSkM9S+gG4qYEZxmA4tmwwBojeCdMtgExxbBgVRQTEOT4U
CAxRMRnZMLAPGCCSxc1MM11YRGOCbhAGuPBlCxoG0tKW1hbkuxSklJZZeq9e
8vJy9+5+9/Lyz9t8NCs+lc1isVJeXOZ73SVzNAHLEhgcZI42a2553gmm1bWg
uhhCLQWw9D+CRxVrBUhPM8cAfDLRvJZb5odAc/dHgznv7TGHBBz3+/z30pMD
W+TJfdmvgO+1Sq3B6ithtK/+oNY6r4LfZiTSIes+NbzxdkXfVas3BvyojNow
qzcOmR3y2CIT401AemmOKxgZbxLO1t4Y2mNgvGnoHo5Tn1pivFmYyXolMknH
9M/CjQTR1bZpxp8DyfJBiNUw81ooE7WV1kmZffMgOe7UFt3HeDoob1fnz7cy
ng7eqZac9/2T8RYg6ufB+S9qGW8R9oz/VR9TzHjPIT/pol/kIcZbAqMiRdz+
AeMtgXf47IXtHMbTQ653A7g+MZG3DIM7SqvPXjGRZwBfbh6nO8lE3gq8v5d9
kuduIm8FwnaXZcW0GMkzwnZNoDD+sJE8E/BDGmL6tCvkrcKWYt7D+IIV6zma
SBkq3DttoH4Wtl987buUfQbyWSjsM3mfu7YMa+NBNphxXf9egkZP+2yw4fdb
C42eevLY+FiwtasrZIk8Nq5eL9wYGfOcPDZ2+v19eCpxkd7HFh/UeDdh6gJ5
ttimb3K+nKcjzw7jxLD+SsU8eXY4WPhHf/dtLXl2GMlhB8Rp5shbhwv7oyNe
D5gjbx1qc5JDPzszS549dmzj/1Nyb4Y8e2wJ2HCpsWaaPHtcDPt3qL18ijwH
lGWWJPrWTZLngPP7fPHzexPkOeItsUMCd/IZeY64q+BCY3jgM/IcUdR97qjL
mXHynLDSr/L2xMAYeU4o6fTUFEeMkeeMd9w9BqpuashzRvUDVXGgj4Y8Zyzl
7NYd/1VNnguWZ6cbS/3V5LngDH/sU8Wdp+RxkHM+9Zpm11PyOKiWZfpuGlGR
x8Fvp+waeV+pyHPFwOjsvGMbVOS5YsjXDarWZiV5XPTJb14vTVOSx0UHmyPp
pwRK8riYI/LYubN1lDw3jMiNiPI+PUqeG7IKej4U+Y2Sx0McfNLrL1PQ+eLh
gaqALNsX/xmLz0Nl7psfhycoaJ6HwrSyyzeFCtrHQ6+tR06HPpaTx8fK5C3x
7dVy8viYrWqdDj4gJ4+PXa8m6nM2y8njY9BbReLOfhl5fKx/uKN6LlVGngBt
e1s+Cl4ZIU+A1Xf9v/nypxHyBLjJLulQ0bsj5AmwJNtreZtYSp4AO34Iq6jJ
kZInRA/D0saTHlL6PkLUzLrF/3JfQp4QozyKM7yOScgTYuwJmX0PV0KeEAVp
/V5V94fJc8eKEJ0k+cQw/A/VkNCi
     "]]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 0.42836805894408136`},
  Method->{},
  PlotRange->{{0, 5.}, {0.42836805894408136`, 2.}},
  PlotRangeClipping->True,
  PlotRangePadding->{{0.1, 0.1}, {0.031432638821118374`, 
   0.031432638821118374`}}]], "Output",
 CellChangeTimes->{3.6243788816696*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1424, 779},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"9.0 for Microsoft Windows (32-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 136, 3, 31, "Input"],
Cell[CellGroupData[{
Cell[718, 27, 570, 16, 113, "Section"],
Cell[1291, 45, 11979, 372, 572, "Input"],
Cell[CellGroupData[{
Cell[13295, 421, 96, 1, 31, "Input"],
Cell[13394, 424, 1615, 26, 152, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15046, 455, 200, 4, 57, "Input"],
Cell[15249, 461, 137, 2, 31, "Output"],
Cell[15389, 465, 133, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15559, 472, 188, 4, 43, "Subsection"],
Cell[15750, 478, 2959, 93, 123, "Input"],
Cell[CellGroupData[{
Cell[18734, 575, 496, 13, 33, "Input"],
Cell[19233, 590, 2061, 40, 246, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21331, 635, 420, 12, 33, "Input"],
Cell[21754, 649, 1929, 37, 245, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)